微动勘探与高密度电法在小尺度老采空区探测中的应用

贾煦， 孙建伟， 刘向东， 程贤达， 黄虹霖

（1.中国地质调查局西安矿产资源调查中心，陕西 西安 710100；2.中国地质大学（北京）地球物理与信息技术学院，北京 100083）

0 引言

地下油页岩开采后会形成一定规模的采空区，采空区沉降或塌陷会对地面建筑和环境产生严重影响。目前，采空区探测主要采用钻探和物探方法[1]。钻探方法精度高，准确直观，但成本较高，对大范围的采空区探测或复杂采空区探测，尤其是缺乏历史资料的采空区探测可行性不高。物探方法经济快捷，不仅可以确定采空区的平面位置，还能直接显现采空区的空间形态特征，近年来得到了广泛应用[2-3]。但是每种方法都有其适用性和局限性，需要根据研究区的地质特征，选择合理的方法或方法组合来进行探测。

栾川县赵庄村一带曾于20世纪有过多期次油页岩开采活动，随着时间推移，采空区地面沉降问题逐渐凸显，多处房屋出现变形开裂现象，对人居安全产生威胁，急需查明采空区分布范围和空间形态特征，以便有效预测地面沉降的发展趋势，控制和减轻地面沉降的危害。然而，由于该区采矿活动历史久远，采掘相关资料缺失，给采空区探测带来一定难度。

微动勘探具有简便、经济、不受电磁干扰影响、对环境及地形适应能力强等特点，对探测村庄覆盖区之下的速度异常区域具有一定的技术优势[4]。前人多是将其用于大规模采空区的探测[5-8]，对于小尺度老采空区的探测报道相对较少。本文针对赵庄村一带的采空区问题，尝试将微动与高密度电法组合进行探测实验，研究其用于小尺度油页岩老采空区探测的可行性。

1 采空区地质背景及物探方法选择

1.1 采空区地质背景

栾川县赵庄村一带属于潭头盆地油页岩分布区，位于秦岭地轴北侧的三级构造中。区内主要出露古近系高峪沟组、潭头组以及第四系洪冲积层（图1）。高峪沟组为一套红色岩系，下部为砾岩与钙质长石粉砂岩不等厚互层，中部为砾岩、含砾砂岩夹粉砂质泥岩，上部为砾岩与钙质长石粉砂岩不等厚互层夹砂岩。潭头组主要为一套灰色岩系，下部以泥岩为主，夹细粒岩屑石英砂岩，中部和上部为薄-中层泥晶灰岩夹泥岩，部分区域夹油页岩及煤线，该套地层是本次勘探对比解释的主要目的层。第四系冲洪积层为含砾粉砂质黏土、含砾砂质黏土夹砾石层。

图1 研究区地质背景与工作部署

研究区所属潭头盆地受扬子板块与华北板块碰撞挤压及西太平洋板块俯冲控制，南部以马超营大断裂为界，东侧受NE向断裂控制，总体呈现向SE倾斜的掀斜式单向断陷盆地，区内发育2组裂隙，分别为69°～90°和340°～360°。

1.2 采空区特征

通过实地调查和走访群众，赵庄村一带采空区有以下4个特点: ①规模较小，高度普遍在2 m左右，宽度约3 m； ②埋深较浅，受限于民采的开采设备，埋深普遍在30 m以浅； ③分布杂乱无序，因开采活动无统一规划，且从20世纪初至20世纪70年代有过不同期次的开采历史，致使采空区分布杂乱； ④相关资料缺乏。

1.3 物探方法选择

目前，采空区探测常用的物探方法可分为3大类，即地震类方法、电磁类方法和放射性探测[9-17]。地震类方法主要根据采空区相对于围岩通常表现为低速异常来识别,电磁类方法主要根据采空区充水与否，表现为高阻或低阻异常来识别,放射性探测则主要根据采空区通常为氡气汇集区，相对于周围表现为高氡气异常来识别。人工源地震因震源激发存在安全隐患和环境噪声污染，考虑到研究区人员聚集的实际情况，该方法不具备实施条件。电磁类方法因受人文电磁干扰影响，效果难以保障。微动勘探观测方法便捷，无需人工震源，不受电磁干扰影响，对环境噪声无特殊要求，因而对工作区的适用性较其他方法更强，本次研究将其作为主要探测方法。为了减少断裂构造或其他地质条件引起的低速异常对微动探测的干扰，同时开展了高密度电法勘探，以期与微动勘探相互验证，提升推断解释的准确性。

由于研究区资料缺乏，因此本区地层电性特征参考的是邻区的以往资料。由表1[18]可见，区内主要岩性层之间电阻率差异均大于20 Ω·m，差异较为明显，因此可首先选择在地表变形区及其邻区部署高密度电法探测，利用该方法工作效率高、成本低的优点快速摸清地下地层结构情况，再结合地质资料，初步判断油页岩分布区域和深度，识别断层等可能对微动推测产生干扰的区域。此外，若采空区未充水则表现为高阻异常，若采空区充水则表现为低阻异常[19]，将其与微动结果相互验证，可降低单一物性参数产生的多解性。

针对采空区规模小、分布杂乱的特点，在由高密度电法圈定的油页岩分布区域开展了微动勘探工作。不论采空区是否充水，其相对围岩都表现为低速特征，再结合采空区埋深浅的特点，因此微动探测将重点关注浅部低速异常区，据此圈定采空区的分布位置和形态。

表1 潭头盆地主要岩性电阻率测试结果[18]

2 物探工作方法与数据处理

2.1 高密度电法

高密度电法是直流电阻率法的一个分支，与常规电阻率法相同，通过发射电极A、B向地下供入电流I，测量M、N极之间的电位差ΔU，从而求得一定深度范围内地层的视电阻率值。为了提升探测效率，高密度电法通过一次布设多个电极，利用仪器控制改变各电极在发射电极和测量电极间的“角色”转换与组合，完成电阻率剖面测量和测深工作，获得相对较多的地层信息，所测数据经过处理与反演获得电阻率断面图。

本次研究使用重庆奔腾WDA-1B（120道）高密度电法仪开展数据采集工作，一次布设电极120根，电极间距为5 m，选用温纳装置，供电时间200 ms，采用滚动测量的方式延长剖面长度，测量层数33层。在实地调查基础上，根据地面房屋变形分布情况，结合地形条件，布设高密度电法剖面2条（图1），点距均为5m，其中G1剖面长940 m，测量点数 4 553个，G2剖面长890 m，测量点数4 224个。

对实测数据进行畸变点剔除、地形校正，利用阻尼最小二乘反演算法进行反演计算，获得相应测线下方的电阻率断面。G1剖面经过4次迭代，拟合误差1.7%，G2剖面经过5次迭代，拟合误差6.45%。图2给出了2条剖面的原始视电阻率拟断面与由反演结果经正演计算得到的视电阻率拟断面的对比，可以看出二者相似性较好，说明反演电阻率模型可靠，可用于后续的地质解释。

图2 高密度电法视电阻率拟断面

2.2 微动勘探

地球上人类活动和各种自然现象引发的微弱震动称为微动，它是一种由体波和面波组成的复杂震动，其中面波能量占微动总能量的2/3以上，信噪比高，在分层的地层介质会携带很多地层介质信息。因面波在一定时空范围内满足统计稳定性，可从观测到的微动信号中提取面波的频散曲线，通过对频散曲线的反演来推断地壳浅部的横波速度结构[20]。其理论基础由Aki[21]和Capon[22-23]于20世纪中叶提出[21-23]，Okada等[24-25]在此基础上于20世纪90年代将其系统化，20世纪80年代，王振东[26-27]将该方法引入我国。微动勘探通常是利用拾震器记录面波垂直分量信号（面波质点位移的垂直分量约为水平分量的1.5倍），再通过数据处理从记录的信号中提取面波频散曲线。

根据高密度电法反演结果和推断解释，在G1剖面选择泥岩、砂岩地层埋深较浅的445～685区段布设微动测线W1，包含测点9个，点距30 m，测线长度240 m; 在G2剖面选择160～550区段布设微动测线W2，兼顾电性模型推断的采空区和泥岩、砂岩埋深较浅区，包含测点14个，点距同为30 m，测线长度390 m。为了降低物探结果多解性，微动剖面与高密度电法剖面相应区段重合。根据前人工作经验，观测系统选用三层嵌套式等边三角形台阵（图3），结合实验结果，确定最小边长15 m，最大边长60 m，每个测点共布设10个拾振器。

图3 三层嵌套式等边三角形台阵

使用北京水电物探研究所WD-1（有线）型设备，动态范围120 dB，2 Hz拾震器进行数据采集。采样过程中，采样率设置为5 ms，因所用仪器可在采集现场实时显示频散曲线，采集过程中频散曲线形态稳定即可终止采集，因此各点连续记录时长略有差异，总体约13 min。采集过程中实时关注现场频散曲线形态，及时调整迭代频率范围，直至频散曲线形态正常，共获得高质量微动测点23个。

目前常用的频散曲线提取方法有频率-波数法（F-K法）、空间自相关法（SPAC法）和扩展的空间自相关法（ESPAC法）[28]。其中，空间自相关法因其精度高、方法成熟，在实践中应用更为普遍，本次野外采集的地震记录即选用该方法获得面波频散相速度数据。图4给出了部分测点频散曲线，图中横坐标vc为相速度，纵坐标L/2为半波长。

图4 微动部分测点频散曲线

对于采空区探测而言，只需了解剖面上横波速度的相对变化，无需精确获得横波速度绝对值[29]，故未进行反演工作，仅对相速度按周期作提高峰度计算，获得视速度

式中:vx,i为周期i的视速度，m/s；vc,i为周期i的相速度,m/s；ti为周期。再结合地形高程数据，即可得到二维视速度映像图。

3 物探成果解释与验证

3.1 高密度电法成果解释

图5、图6分别为G1、G2测线反演电阻率断面和相应的地质解释。可以看出，该区域地下电性结构总体呈现出较好的层状特征，由浅至深电阻率逐渐升高，电阻率范围为3～500 Ω·m，这与该区地质情况相符。浅部电阻率小于10 Ω·m的薄层推断为第四系冲洪积层，厚度约10 m，局部有增厚现象； 浅覆盖层之下，电阻率10～100 Ω·m电性层推断为新生界古近系砂、泥岩互层，也是油页岩的可能分布层位； 深部电阻率大于100 Ω·m的电性体推断为古近系砂岩层。

G1剖面中，440～520 m和640～720 m的砂泥岩层中出现电性梯度带，表现为相对低阻特征，且体现出一定的规模，结合民采活动采矿深度和规模有限（30 m以浅范围）的情况，推断为2条隐伏断裂破碎带。G2剖面160～200、580～600、620～760 m区段出现电性梯度带，表现为相对低阻特征，且体现出一定的规模，与G1剖面同理，推断其为3条规模、产状不同的隐伏断裂破碎带； G2剖面260～300 m区段的低阻异常可能为采空区充水引起。2条剖面的其他区段未见明显的电性异常，无法判断采空区的分布范围，推断可能与采空区规模小、塌落、碎石充填以及围岩电性差异不明显有关。

图5 G1剖面高密度电法反演结果（上）与地质解释（下）

图6 G2剖面高密度电法反演结果（上）与地质解释（下）

3.2 微动勘探成果解释

图7为W1、W2剖面的二维视速度（vx）映像图。可以看出，探测区域地下结构总体呈现层状结构特点，垂向上随着深度增加，视速度逐渐增高，横向波速变化不明显，等值线的起伏基本反映沉积地层界面形态，这与高密度电法的探测结果有着较好的对应关系。

图7 微动剖面二维视速度映像

W1剖面二维视速度映像图中，445～545区段，标高550 m下方出现一系列“串珠状”低速异常，这与根据电性特征推断的隐伏断层F1的空间位置有较好的对应关系，推断低速异常带为断层破碎带的反映； 645～685区段，标高560 m下方的“串珠状”低速异常带与推断的隐伏断层F2的空间位置有较好的对应关系，推断此低速异常带为断层破碎带的反映； 在535～580区段，标高560～540 m深度范围内横向上相对于周边区域出现明显的“圈闭状”低速异常，该低速异常位于根据电性结构推断的砂、泥岩互层（油页岩赋存地层）内，且该区段砂泥岩层埋深较浅（约5 m），便于民采活动，结合该段地面房屋变形严重的实际情况，该区为油页岩采空区的可能性较大。综合高密度电法和微动勘探结果，在W1剖面535～580区段，标高540～560 m 深度范围内可能存在采空区。

W2剖面上，根据低阻异常推断的采空区在视速度映像图上表现为横向延伸较大的宽缓低速异常，推断该低速异常可能为古近系泥岩、砂岩地层出现凹陷，使得第四系冲洪积层增厚所引起。在220～270 m区段，标高560 m下方出现近直立的低速异常，与根据电性特征推断的产状较陡的砂岩地层在空间位置上有较好的对应关系，分析认为出现该低速异常是因为该区域地下地层产状较陡，不满足微动勘探适用于水平层状介质的理论假设，造成对相速度的低估而出现的低速假异常。综合高密度电法和微动勘探结果，W2剖面覆盖区域存在采空区的可能性较小。

3.3 钻探验证

根据物探结果，推断G1剖面535～580区段，标高540～560 m深度范围内可能存在采空区。为了验证这一推断的可靠性，结合地面房屋分布的实际情况，在G1测线575 m处实施了GZK-01，终孔深度51 m。孔内0～6 m为第四系粉质黏土，其下为泥岩、粉砂岩互层，局部深度见煤线及油页岩层（图8）。钻孔施工至标高561 m时，出现钻进突然加速，钻井液大量漏失的现象，标高558～561 m段的岩心呈黑色碎块状，采取率不足。钻进情况和岩心特征证实该深度范围为油页岩采空区，且采空区出现塌落情况，被碎石充填，但未充水。推断结果得到验证，证明了该组合方法探测小尺度老采空区的有效性。

图8 GZK-01钻孔柱状图Fig.8 Columnar section of Bore Hole GZK-01

3.4 讨论

在利用组合方法对采空区进行探测的过程中，高密度电法和微动各自发挥着不同的作用。

高密度电法快速获得了较高精度的地下电性结构信息，为地层划分、地层产状判断和构造识别提供帮助。但同时也注意到，由于小尺度采空区规模较小，尤其是老采空区发生塌落充填后，与围岩电性差异降低，其在电性结构上难以引起可识别的异常,如G1剖面中的采空区段在电性结构上毫无反映。此外，仅仅依靠电性异常容易引起误判，例如G2剖面中局部地段中的低阻异常，后经钻孔揭露实为砂、泥岩层下凹，第四系盖层增厚的反映。

微动勘探较好地识别了小尺度采空区引起的低速异常，但其同时容易受构造和陡立地层引起的低速异常的干扰，如W1剖面445～545 m区段的隐伏断层产生一系列“串珠状”低速异常，W2剖面220～270 m区段陡立地层引起的低速异常。

此外，还有一个问题需要注意，2种方法综合应用虽然可有效识别出小规模老采空区，但在采空区底界面的推断上还不够准确。如GZK-01揭露的采空区规模在纵向上小于W1剖面低速异常推断的采空区规模，存在“下拉”现象，无法准确圈定采空区规模，仍需借助钻探工作验证，这可能与其“体积效应”有关。解决这一问题需开展进一步的探索研究工作。

4 结论

（1）利用微动和高密度电法组合勘探，可有效探测小尺度油页岩老采空区的分布位置，采空区在面波二维视速度映像图上表现为低速异常。

（2）高密度电法在组合方法中主要用于地层划分，缩小微动勘探的目标区域，排除构造和陡立地层引起的低速异常干扰。

（3）该组合方法虽能推断油页岩采空区的分布位置，但在采空区底界面的判断上还不够准确，相较于实际情况，底界面下延，需借助钻探工作确定底界面深度，才能准确圈定采空区规模。

致谢：本项目在微动数据采集和数据处理方面得到北京水电物探研究所刘云祯所长、金荣杰老师和该所技术人员的大力支持和无私帮助，此外，项目组成员崔鑫、王延州、李宏等同志在整个施工过程中付出了艰辛的努力，在此一并表示感谢！