1、开究页目2023May年5 月GEOLOGREVIEW质地评论No.Vol.69第3 期第6 9 卷反射槽波探测煤层断层构造的应用研究杨顺1，2），孙文斌3），朱广辉3），韩建光1)1)中国地质科学院，北京，1 0 0 0 37；2）中国地质大学，武汉，430 0 7 4；3)北京探创资源科技有限公司，北京，1 0 0 0 6 7WWW 0 1 6）。随着煤田开采越来越趋于高效化和经济化，煤层的开采效率很大程度上依赖于地下地质条件的勘查程度。目前阻碍开采进度或造成事故的煤层地质问题多是由于探测时未查明的小构造所造成的（乐勇等，2 0 1 3），如小断层、陷落柱和褶曲等都会影响煤层工作面回采，因
2、此查明地下煤层中小构造，对施工设计、储量估计、地质灾害预防以及煤矿开采工作至关重要。三维地震勘探是目前煤田采区最有效的勘探技术之一，在煤田勘探中已经得到广泛应用，并取得了良好的成效（武磊彬等，2 0 0 9），但地面地震勘探技术在精细探查煤层构造时存在明显缺点，如高频信号在地层的吸收衰减、煤层间的多次波干扰现象、能量扩散等问题，同时三维地震偏移速度精度不够、小断层反应模糊以及地震时一深转换速度不准等，导致地面地震勘探在勘查煤层时分辨率不够高（程建远等，2009），地震解释成果无法完全满足生产部门的生产要求。相比于地面地震勘探，槽波地震勘探是利用在井下煤层开采工作面中激发、传播和接收的导波来探查
3、煤层地质构造以及岩性变化的地球物理方法（姬广忠，2 0 1 7），具有高效、准确、分辨率高的优点，其在煤矿勘探中探测距离大、抗干扰能力强、波形明显（胡国泽等，2 0 1 3），槽波地震勘探能够准确探明煤层空间中的不连续性、断层、陷落柱、冲刷带等异常体，还可以获得工作面内部的其他地质信息（韩德品等，2 0 0 9），在探测小型构造时具有显著的优势。Evision（1 955）在新西兰煤矿的煤层中首次发现了通过煤层制导的槽波，阐述了槽波在煤层勘探中应用前景；Krey（1 96 3）从理论上计算了槽波的频散关系，奠定了槽波地震勘探的基础；随后越来越多的学者相继对槽波进行了研究，推动了槽波勘探的发展。
4、刘天放等（1 994)在其编著的槽波地震勘探中详细的论述了槽波相关理论，提出了提取频散曲线的方法；李天元等（1 997）研究了我国槽波的赋存状态，并将其分为四类，还研究了通过钻孔实现槽波地震勘探的技术；杨真等（2 0 1 0)研究了SH型槽波的频散特征及波形模式和薄煤层中透射槽波的时频域特征、速度特征；姬广忠等（2 0 1 2）对煤矿井下槽波进行了三维数值模拟，发现由于巷道的影响，巷道壁上产生很强的巷道振型槽波，煤层中则出现了以Love型为主的槽波,并分析了实际槽波的形成机理；皮娇龙等（2 0 1 8）开展了煤层内槽波地震勘探注：本文为国家自然科学基金资助项目（编号：42 1 7 41 57）
5、、中国地质科学院基本科研业务费项目（编号：JKY202216）和企业合作开发（编号：56 0 0 B2021000004）的成果收稿日期：2 0 2 2-0 6-2 8；改回日期：2 0 2 2-1 1-2 1；网络首发：2 0 2 2-1 2-2 0；责任编辑：刘志强。Doi：1 0.1 6 50 9/j.g e o r e v i e w.2 0 2 2.1 2.0 2 1作者简介：杨顺，男，1 996 年生，硕士研究生，地球物理专业；Email：1 6 0 37 58 344 q q.c o m。通讯作者：韩建光，男，1 98 7 年生，博士，副石员，主要从事地震波偏移成像及深部探测方法
6、研究；Email：h a n j i a n g u a n g 6 1 3 1 6 3.c o m。103220233年质论地评的数值模拟和物理模拟的初步研究工作。目前，槽波地震勘探技术已经较为成熟，可以成功的探测出煤层中的不连续地质，如断层或断裂带、泥床分布和冲刷带等,是一种有效的小构造探测方法（Mechie etal.，2 0 1 9）。槽波勘探主要有反射法、透射法以及反射一透射联合勘探法。槽波透射法主要应用于煤层厚度变化情况、夹研厚度分布情况、空区、空巷探查等（程久龙等，2 0 1 4）。槽波反射法主要探查煤层中的各种大小断层、侵人体以及岩体等异常体反射槽波勘探精度高，抗干扰能力强，同
7、时还具有强振幅、波形易于识别、探测结果直观等优点，在多个矿区已经得到了应用，适合高精度槽波地震勘探（杨思通等，2 0 1 2）。槽波探测的精度不但需要高质量的采集数据，而且取决于数据处理阶段参数的确定、资料分析和处理手段。目前反褶积压缩波列、径向道变换压制直达波、槽波绕射偏移成像以及包络叠加方法均能有效提高反射槽波的成像精度。本文利用反射槽波勘探方法对我国山东某矿区煤层小断层构造进行应用研究。首先简单介绍反射槽波地震勘探的基本原理，给出了反射槽波数据精细处理流程，并将其应用于实际煤层小断层勘探，对实际资料中的反射槽波进行频散分析、单炮记录分析，采用基于包络叠加的方法来计算槽波波列的能量，提高成
8、像精度，最后有效识别出煤层中的断层构造，验证了方法的有效性和实用性1反射槽波的基本原理在实际煤系地层中，煤层的速度和密度低于顶底板的围岩，所以煤层是夹在顶底层间的一个低速层，形成了波阻抗差异，当在煤层中激发震源产生地震波，就会在顶底板之间形成多次反射波，不同类型的层间多次波相互干涉叠加形成槽波（刘天放等，1994；乔勇虎等，2 0 1 8）。当煤层工作面中有断层、夹研、陷落柱等不连续的地质构造体时，会形成强烈的波阻抗差异界面，震源激发的槽波在传播过程中遇到异常构造后形成反射波，进而传回巷道内的检波器被接收，即反射槽波（王季，2 0 1 5；杨辉，2018），反射槽波形成原理如图1 所示。若煤层
9、中有较为明显的异常时,则反射槽波信号较强，若煤层中异常构造体小或不突出，则反射能量弱或没有，因此可以根据反射槽波的有无或者强弱来判断煤层内部构造情况。2反射槽波数据处理流程及关键技术运输巷道煤层工作面异常构造体检波点炮点通风巷道图1 反射槽波的形成原理示意图Fig.1 The formation principle of the reflectedin-seam waves原始数据1、数据和记录格式转换2、剔除坏道3、建立观测系统预处理4、初至校正5、能量补偿1、频散分析单炮记录分析2、频谱分析3、时频分析4、速度分析滤波5、育能量分析包络计算叠加处理偏移成像图2 反射槽波数据处理流程图Fig
10、.2 Reflected in-seam waves data processing flowchart2.1处理流程数据的分析、处理是槽波地震探测技术的关键。槽波能量主要集中在煤层中并沿着煤层以柱面波的形式传播，频率越高，在煤层中的能量越集中。当遇到断层或其他异常构造时，槽波的能量将发生显著的变化。因此,提出针对槽波能量处理的一系列措施，从而反演出煤层中的断层或破碎带等构造异常体。反射槽波数据处理流程如图2 所示，包括预处理、单炮记录分析、滤波、包络计算、叠加处理、偏移成像等。2.2关键技术2.2.1频散分析槽波最显著的特点就是频散，即槽波的速度随频率而变化（DavidandStewart，
11、1 97 9；任亚平，2015）。在常规地震勘探中，地震波的速度与传播介质的密度、弹性参数、传播方向有关，但与地震波1033第3 期杨顺等：反射槽波探测煤层断层构造的应用研究的频率无关,而槽波在煤层中的波速是由本身的频率以及煤厚变化所决定的。槽波由不同的层间多次波干涉叠加而形成，在传播过程中，不同频率的简谐波波速不同，因此在同一时间的传播距离各有差异。当传播一定时间后,短脉冲简谐波会变为长波列,最后形成槽波的波形在传播过程中不断变化的状态，即槽波的频散（SchottandWaclawik，2 0 1 5；郭银景等，2 0 2 0）。槽波的频散特征对求取煤层速度、围岩信息以及煤层速度有着关键作用
12、，在槽波数据解释时，频散特征不仅能够证明槽波的存在，也可以用来解释煤层地质信息。频散分析实际上就是求取槽波记录中的群速度与相速度曲线,Love 型槽波的频散方程为：-inN+1hZhikBN+1ei+1二N-inN+IkZhi(-inN+IkuN+IB+IeikA。avan-ian-2a,二(inokuoAo)(mllm12ik:Ao(1)m21m22(inokuoAocosO,in,u,sino,其中,a,=(n=1,(in.u,sing,coso,2,3N），A,B为任意系数；N为地层数目;h：层厚；u，：第n层的剪切模量。2.2.2槽波包络叠加成像煤矿井下槽波数据通常噪音大、信噪比低，需
13、要压制噪音，包括工频噪音压制、声波噪音压制、随机噪音衰减等；由于槽波波列长，需要研究反射槽波特征函数，对槽波波列进行压缩；槽波波场含有纵波、直达槽波、反射槽波等多种波，采用频率滤波、极化滤波等方法，对槽波波场进行分离；反射槽波能量较弱，需要消除几何扩散能量衰减影响，增强反射槽波。(1)包络计算。槽波的整体频率偏高，尤其是在埃里相频率周围；此外,由于频散效应的影响,波列之间相隔甚远，相邻地震道的槽波相位不能对比，进行相应的叠加处理时会引起反向消除。因此，需要在对槽波波列做适当处理之后再进行叠加和速度分析。信号包络就是常用的处理方法之一，该方法是将波形序列转换为能量序列，转换之后的信号在外形与之前
14、的信号一致,转换使得信号频率降低（赵立明等，2 0 1 4）。实测的地震道可看成一个复地震道在实平面上的投影，称为实地震道；复地震道在虚平面上的投影称为虚地震道，它与实地震道互为正交地震道，复地震道可表示为：c(t)=x(t)+ix;(t)(2)式中，c(t）表示复地震道；x(t）为实地震道；x（t）为虚地震道。根据希尔伯特变换，虚地震道可直接从实地震道求取，表达式如下：1(3)一*TTTtTT-8t-T式中，x(t）称为x(t）的希尔伯特变换,实际上相当于对信号道进行一次90 纯相位滤波。其滤波因子及频率响应分别是：1h(t)Tt(4)一f0if0根据以上推导，可以得到实地震道x（t）和虚地
15、震道x（t），则地震道的包络E（t）可可按下式计算：E(t)=/x(t)+x;(t)2(5)地震道包络可理解为地震波的瞬时振幅，它把地震道的概念向复数域扩展。求取反射槽波包络的过程实现了从波形序列向能量序列的变换，变换之后信号仍然与之前保持一致。反射槽波包络频率较低，相干性强，在叠加之后效果明显，能够实现反射槽波的精确成像。(2)叠加处理成像。共深度点（CDP）或共反射点叠加技术是地面地震勘探的核心方法，它以不同的炮检距对同一反射点进行多次重复观测，然后把同一反射点的CDP道集经过正常时差校正或动校正后进行叠加，这种方法对于提高信噪比、压制多次波及突出反射波有明显效果。反射槽波信号一般较弱，记
16、录的信噪比低，严重影响到反射槽波成像。因此，本文采用多次叠加技术处理反射槽波。槽波是频散波,相速度与群速度存在明显差异且依赖于频率,但是槽波埃里相具有速度变化小的特点,因此可以通过估算出埃里相速度来进行动校正叠加处理，如图3所示，运用槽波埃里相的特性，估算出速度值进行多次叠加后，反射槽波同相轴有了明显的增强，3实际煤田应用分析3.1工区概况我国山东某煤矿530 5工作面煤层厚度最大103420233年质论地评水平距离（m）6012018024030036050-100送1 50200-250图3反射槽波的叠加同相轴Fig.3 The superposition event of the reflectedin-seam waves运输顺槽工作面轨道顺槽图4工作面测点布置图Fig.4 Arrangement measurement pointlayout of the work-plane7.4m,最小3.4m，平均7.7 m，属结构简单、稳定的中厚煤层。顶板、底板岩性为泥岩、细粒砂岩、粉砂岩、砂质泥岩。顶底板结构稳定，与煤层形成了良好的低速夹层，为槽波的形成创造了有利条件，工作面内断层
