三维非常规被动源面波勘探在露天矿采空区探测中的应用研究.pdf

第 卷第期 年月山西煤炭 文章编号：（）犇 犗 犐：三维非常规被动源面波勘探在露天矿采空区探测中的应用研究王立熊（中煤平朔集团有限公司 地质测量中心，山西 朔州 ）摘要：露天煤矿采空区的常规物探方法易受场地环境、地质条件等的影响，探测效果欠佳。以安家岭露天煤矿为例，通过改进传统的被动源面波勘探方法，采取呈三条平行线展布的检波点阵列记录被动源面波信号，并利用平面波基函数直接扫描不同时间窗口的原始信号，来获取与地下介质有关的面波频散信息，然后反演频散曲线，得到测区的横波速度剖面图，根据剖面图低速特征圈定出采空区范围与边界。通过钻孔验证发现，探测异常与实际采空区吻合，说明该方法在煤矿采空区探测中有效，为今后煤矿采空区的探测提供了指导。关键词：露天煤矿；采空区；面波勘探；基函数；频散曲线中图分类号：文献标识码：开放科学（资源服务）标识码（犗 犛 犐 犇）：犃 狆 狆 犾 犻 犮 犪 狋 犻 狅 狀狅 犳 犇犝 狀 犮 狅 狀 狏 犲 狀 狋 犻 狅 狀 犪 犾犘 犪 狊 狊 犻 狏 犲犛 狅 狌 狉 犮 犲犛 狌 狉 犳 犪 犮 犲犠 犪 狏 犲犈 狓 狆 犾 狅 狉 犪 狋 犻 狅 狀 犻 狀犌 狅 犪 犳犇 犲 狋 犲 犮 狋 犻 狅 狀犻 狀犗 狆 犲 狀 狆 犻 狋犕 犻 狀 犲犠犃 犖 犌犔 犻 狓 犻 狅 狀 犵（犌 犲 狅 犾 狅 犵 犻 犮 犪 犾 犛 狌 狉 狏 犲 狔犆 犲 狀 狋 犲 狉，犆 犺 犻 狀 犪犆 狅 犪 犾犘 犻 狀 犵 狊 犺 狌 狅犌 狉 狅 狌 狆犆 狅，犔 狋 犱，犛 犺 狌 狅 狕 犺 狅 狌 ，犆 犺 犻 狀 犪）犃 犫 狊 狋 狉 犪 犮 狋：，；，；，犓 犲 狔狑 狅 狉 犱 狊：；收稿日期：作者简介：王立熊（），男，山西运城人，本科，工程师，从事煤矿地测技术研究工作，：引文格式：王立熊三维非常规被动源面波勘探在露天矿采空区探测中的应用研究山西煤炭，（）：，（）：工程概况平朔矿区煤炭资源丰富，埋藏相对较浅。经多年开采，平朔矿区形成了分布众多、形状各异、大小不一的小窑采空区，给后续规模化开采带来了严重的安全隐患。安家岭露天煤矿界内有多座已关闭的小窑，关闭时间久远，现有的采空区资料可靠性较差，无法准确指导生产。因此，在开采施工之前查明采空区的赋存情况是十分重要的。目前探测煤矿采空区主要依靠地面瞬变电磁法、地面高密度电法、浅层反射波地震法等物探方法。安家岭露天矿采剥区内，原始地表已被剥离成台阶状，且现场布置有大量大型设备。因此，地面瞬变电磁法和地面高密度电法受安家岭露天煤矿现有开采与场地环境条件限制，无法开展探测；安家岭露天煤矿采空区埋深为 ，地震反射波与面波在时间窗口上相互重叠并形成强烈干扰，同时场地表层剥离施工作业与运煤车均产生强烈的振动噪声干扰，极大地降低了反射波的成像效果，所以浅层反射波地震法也不适用。为解决安家岭露天煤矿浅部煤层采空区探测难题，结合场地环境噪声条件，在此采用被动源面波法进行勘探。研究现状面波勘探也称弹性波频率测深，是近 年发展起来的一种地球物理勘探技术，具有简便、快速等优点。在许多领域得到应用，并取得了良好的应用效果。根据源的不同，面波勘探可分为天然源面波勘探和人工源面波勘探 。主动源面波勘探由稳态发展到瞬态，再由瞬态的表面波谱法发展到多道瞬态面波勘探法，勘探精度及勘探效率也在逐渐提高，应用领域也越来越广泛。近些年，随着计算机的发展和面波研究的不断深入，面波勘探技术在工程实践中取得了较好的成果及显著的社会效益和经济效益 。与主动源不同，被动源法在工程应用方面起步较晚，但发展速度较快，方法也越来越成熟。天然源面波勘探在工程安全评价、施工规范设计中发挥了不可替代的作用，它突破了源的限制，使得在嘈杂的城市环境中也能取得较好的勘探效果，但是将天然源面波勘探技术用于煤层采空区探测，国内还处于尝试和研究阶段。年，天然源勘探方法首次应用于煤矿采区勘探，结果表明，与其他物探方法相比，该方法具有明显的技术优势和应用前景。年，王建文等将双源面波勘探技术应用于煤层采空区探测，在国内属于首次尝试，在国外尚无先例，勘探结果表明，该方法是可行、有效的。勘探原理 面波勘探的基本原理地表及地下微震动会产生横波和纵波，横波和纵波相互叠加干涉将使地下介质按一定轨迹运动，这种新形成的波称为瑞雷波。瑞雷波的能量基本限制在一个波长范围内，强能量主要集中在地表附近，且频率低、速度低。瑞雷波的特点为：在自由表面附近沿波传播方向的垂直平面内，瑞雷面波质点运动的轨迹是椭圆的，椭圆的水平轴和垂直轴长度比为 ，且质点的垂直位移比水平位移相位超前。质点的震动轨迹和振幅随着单位波长深度的变化而变化，当质点震动波长与深度狕的比值小于 时，质点运动为逆时针转动的椭圆；当质点震动波长与深度狕的比值大于 时，质点运动轨迹为顺时针转动的椭圆，且振幅随着距自由表面距离的增加而衰减。面波具有频散特性，不同的频率具有不同的面波相速度，这是面波勘探的理论基础。依据面波的特点，通过测定不同频率面波的相速度狏，即可了解地下构造的有关性质并计算相应地层的动力学特征参数，达到地质勘探的目的 。改进后的被动源面波勘探原理传统的被动源面波方法包含两项基本假设：噪声源是围绕检波点呈 均匀、随机、非相干分布的，即空间域平衡假设；长时间的采集策略会使从被动源面波记录中恢复出的经验格林函数趋于平稳，即时间域平稳假设。一般来说，这两项假设是很难被满足的，尤其是当利用环境中频率大于 的噪声源时，基本假设与真实世界中的噪声源条件不相容的程度会更剧烈，而绝大多数的煤矿采空区探测正是利用这一频带的噪声源。为了降低传统被动源面波方法的基本假设与现实条件的不相容程度，即在利用被动源面波信号探测和识别煤矿采空区的过程中，充分考虑噪声源随时间变化的非平稳特征，在此对传统方法作以下改进：采取呈三条平行线展布的检波点阵列记录被动源面波信号；利用平面波基函数直接扫描不同时间窗口的原始信号，从而揭露噪声源方向的时变特征，以及与地下介质结构有关的面波频散信息。改进的被动源面波方法面波传播的过程中存在着两种扩散模式，分别第期王立熊：三维非常规被动源面波勘探在露天矿采空区探测中的应用研究为柱面波与平面波模式。当远场条件被满足时，则可以认为面波以平面波模式扩散。若忽略地下介质的衰减性质，此时，某检波点在圆频率处记录的频率域面波信号狌可以表达为：狌犃（）犖犼犉 狊犼，（）犲犻 狊犼犽狅狉狊（）犼狉狊槡犼（）式中：犉（狊犼，）为第犼个震源狊犼的振幅，；狊 犼为震源狊犼的初始相位，；狉狊 犼为震源狊犼到检波点的距离，；犽狅为与地下介质结构有关的面波波数，犽狅犮；犮是面波的相速度，。从式（）可以看出，检波点记录的信号是犖个震源产生信号的叠加。被动源面波方法本质上就是通过某种波场变换方法从信号中提取与地下介质结构有关的面波频散信息，即面波频散曲线。上文中提到，以互相关为核心的经典被动源面波频散分析方法包含两项基本假设，而这两项基本假设在现实中是很难被满足的。此时，噪声源的方位角和时变特性必须在数据分析时予以充分考虑。和 （）指出，二维的检波点阵列能够感知波源的方位角。因此，本文采用呈三条平行线展布的检波点阵列记录测区中的被动源面波信号，且需考虑噪声源方位角的时变性，采用对原始信号分时窗处理的方法获得噪声源随时间的变化规律。假设平面波扩散是被满足的，则对于某一时窗信号，使用平面波基函数对其进行扫描，就可以获得某一方位角的能量谱矩阵，其表达式为：犈，（）狏犕犕犼犲犻 犼，（）狏狌犼（）狌犼（）（）式中：为方位角，（）；狏为面波相速度，；犕为检波点数量；狌犼（）为某一检波点的频谱；犼（，狏）为以某一点为原点至检波点处的时移，其表达式为：犼，（）狏狓犼 （）狏狔犼 （）狏（）式中：狓犼与狔犼是第犼个检波点的横坐标与纵坐标。从式（）可知，矩阵犈（，狏）是三维矩阵，需要关于相速度和频率作叠加才能得到该时窗信号的方位角能量分布向量犘（），这个叠加过程可以描述为：犘（）犉犽犖犺犈 犽，狏（）犺（）式中：犉为信号有效带宽内的频点数量；犖为相速度扫描的数量。方位角能量分布向量犘（）中最大值所在的方位角就是面波源的方位角。面波源方位角被确定后，就可以按公式（）计算信号的频散谱犈（，狏），表达式为：犈，（）狏犈，（）狏（）对原始信号所有时窗的频散谱作叠加就能够得到最终的叠加频散谱犈（，狏），然后按最大值原则从叠加频散谱的图像中拾取与阵列下方结构有关的面波频散曲线。模型模拟试算为了验证本文提出的噪声源方位角时变分析与面波频散测量方法的可靠性，以两层递增型地质模型为例，基于平面波扩散假设，分别合成了噪声源方位角时不变与时变情况下的噪声记录。模型的地质参数如表所示。表两层递增型地质模型参数 层序号纵波速度狏（）横波速度狏（）密度（）层厚度犺 噪声源方位角时不变条件下的检波点与噪声源配置示意图如图所示。噪声源的数量为 ，随机分布在距离原点 ，方位角为 的范围内；以三条平行线展布的方式布设 个检波点，每条测线都有 个检波点，点距为，测线间距为。图右上角给出了局部放大后的检波点阵列示意图。图噪声源方位角时不变条件下的噪声源与检波点配置示意图 噪声源子波的主频随机分布在 之间，振幅随机分布在 之间，模拟获得时间长度为 的噪声记录。按公式（）（）对合成的噪声记录分时窗扫描即可得到噪声源方位角能量谱的时变图，如图（）所示。从图中可以看出，在 方位角处一直存在明显的能量聚焦，这与模拟中的噪声源条件是高度一致的，说明本文提出的方法可以有效地山西煤炭第 卷从信号中提取出噪声源方位角的时变特征。最后，通过在图（）中拾取每个时窗的噪声面波信号的方位角，利用公式（）就得到了如图（）所示的面波叠加频散谱。图中的黑色点线是由多模式快速矢量传递算法计算得出的两层递增型地质模型的基阶面波频散曲线的解析解。叠加频散谱能量的最大值已与基阶面波频散曲线的解析解吻合良好（图（），这说明本文提出的噪声记录处理与分析方法对于测量与地下横波速度结构有关的面波频散信息是有效的。图噪声源方位角时不变条件下的噪声源方位角能量谱时变图与面波叠加频散谱 噪声源方位角时变条下的噪声源与检波点配置示意图，如图所示。合成记录的过程中，噪声源的数量与检波点的坐标参数保持与上一测试相同；时，噪声源随机分布在距离原点 ，方位角为 范围内；时，噪声源随机分布在距离原点 ，方位角为 范围内。按公式（）（）处理合成记录得到了如图（）所示的噪声源方位角能量谱时变图。从图中可以发现，时间段和 时间段的噪声源方位角均已被精确地估计，对于噪声源方位角发生改变的时间点的描述非常准确。从图（）中拾取图噪声源方位角时变条件下的噪声源与检波点配置示意图 每个时间窗口的噪声源方位角，按照公式（）就得到了合成记录的面波叠加频散谱，如图（）所示。黑色的点线是两层递增型地质模型的基阶面波频散曲线解析解，已与合成记录的面波叠加频散谱的能量最大值精确对应。这进一步说明了本文提出方法的稳定性与精确性，且该方法对于噪声源方位角时不变与时变两种情况均是适用的。图噪声源方位角时变条件下的噪声源方位角能量谱时变图与面波叠加频散谱 第期王立熊：三维非常规被动源面波勘探在露天矿采空区探测中的应用研究实例分析为了验证该方法的效果，在安家岭露天煤矿中进行了实践应用。实验参数为：接收道距，线距（根据现场地形条件布置），接收道数 道（道条线），采样率 ，记录长度，观测时间，检波器频率。测线布置示意图如图所示。图测线布置示意图 首先对特定子接收排列的时域记录和频谱进行分析，以确定被动源面波的带宽，从而进行频散测量。图和图分别显示了 的子阵记录及其振幅频谱，振幅被归一化为轨迹。从图中可以看出，记录的相对振幅随着时间的推移有一定的变化，这表明本次调查的噪声源具有一定程度的非稳态和时变的方位特性。同时，位于图底部的该记录的第一条轨迹在振幅强度上有明显的变化，这也表明噪声源具有很强的非稳态性。图中被动源面波的能量主要集中在 以下，低频大约达到。然而，在 范围内仍存在微弱的能量，同时为了避免在接收间隔的情况下出现空间混叠，所以将频散测量的带宽确定为 。图子排列记录的 被动源面波 图 被动源源面波记录振幅谱 以作为分割子阵记录的窗口，计算每个时间段的方位角能量分布向量，叠加后得到图（）所示的方位角能量谱。从图中可以看出，本次调查的噪声源具有很强的方位角时间变化特征，动态变化范围从 ，这也说明在非稳态现实条件下处理被动源面波图时，考虑噪声源的方位角非常重要。从方位角能谱中提取噪声源的时变方位角，计算每个时间段的频散能谱，然后根据公式得到最终的叠加频散能谱，如图（）所示。如上所述，图（）中的两条白线限制了频散能的可靠范围。然后，使用 算法对提取的频散曲线进行反演。另外，算法可以对基模式和多模式频散曲线进行联山西煤炭第 卷合反演，但由于后者在本调查区域不稳定，所以只选择前者进行反演。初始模型被定义为 层均匀的半空间模型，横波速度为 ，基于频散曲线。图（）给出了由初始模型和反演模型计算出的预测值与子阵记录的实测值之间的对比，可以看出反演后的频散曲线与实测频散曲线拟合得非常好。图（）给出了 横波速度结构满足 探测深度的反演结果，显示在 深度范围内存在明显的低速异常，这可能是采煤区的标志。图子阵列记录的频散测量和反演结果 将获得的频散曲线进行反演得到了测区的横波速度剖面图，如图所示。通过分析发现，在 的水平范围内，推测煤层被部分开采，形成小规模的采煤区，有一定的低速异常，对上覆岩层有一定的扰动。在 的水平范围内，假定煤层完全被采空，形成大规模的采空区，有明显的低速异常（大部分地区速度低于 ），对上覆岩层有较大影响。大面积的采空区很容易造成上覆岩层的应力平衡被破坏，在重力作用下发生断裂，导致剪切波速度下降，分层不良。这种现象在 水平范围内的覆土层中得到了充分体现。在 的水平范围内，推测煤层未被开采，根据低速异常的向上倾斜的几何形状，解释为一条显示为黑色铁路线的巷道，并与煤矿开采区相连。上述结果与收集到的地质和钻孔数据基本一致，说明改进的被动面波法对煤矿采空区的探测是有效的。此外，受限于表土剥离作业的实际情况，在 处布置了一个验证钻孔。利用验证钻孔，在 深处发现了高的空洞，这与该处开采的解释结果一致。图横波速度剖面图 结论采用三维非常规被动源面波探测法，通过改进传统的被动源面波勘探方法，采取呈三条平行线展布的检波点阵列记录被动源面波信号，并利用平面波基函数直接扫描不同时间窗口的原始信号，来获第期王立熊：三维非常规被动源面波勘探在露天矿采空区探测中的应用研究取与地下介质有关的面波频散信息，然后反演频散曲线，得到测区的横波速度剖面图，根据剖面图低速特征圈定出采空区范围与边界。通过模型模拟试算及安家岭露天煤矿现场实验验证发现，在场地环境、地质条件复杂的区域，采用改进后的被动源面波勘探方法探测露天煤矿浅层采空区，效果明显、可靠，该方法可在同类型的采空区探测中推广应用。下一步研究将主要结合露天煤矿采空区特点，充分利用场地强烈的振动噪声，进一步提高探测成果的验证率。参考文献：王振东面波勘探技术要点与最新进展物探与化探，（）：，（）：，王建文，孙秀容，王宏科，等双源面波地震勘探在煤层采空区探测中的应用工程地球物理学报 ，（）：，（）：王光文，王海燕，李洪强，等面波勘探技术的研究现状及进展地球科学前言，（）：，（）：朱鑫面波勘探在工程勘察中的应用世界有色金属 （）：，（）：辛永祺多道瞬态面波技术在岩土工程勘察中的应用山西建筑 ，（）：，（）：毋光荣，余凯，马若龙天然源面波勘探技术在工程中的应用研究工程地球物理学报 ，（）：，（）：肖柏勋，李长征瑞雷面波勘探技术研究评述工程地球物理学报 ，（）：，（）：刘庆华，鲁来玉，王凯明主动源和被动源面波浅勘方法综述地球物理学进展 ，（）：，（），（）：黄真萍，朱鹏超，胡艳主动源与被动源面波勘探方法对比分析与应用路基工程 （）：，（）：赵东被动源面波勘探方法与应用物探与化探，（）：，（）：（编辑：安娜）山西煤炭第 卷