煤田岩性地震勘探技术进展

1、2020/8/23,中国矿业大学,1,煤田岩性地震勘探技术进展,崔若飞 中国矿业大学 二零一一年十月,2020/8/23,中国矿业大学,2,汇报提纲,煤矿安全开采与岩性地震勘探 地震属性技术 地震反演技术 地震相分析技术 谱分解技术,2020/8/23,中国矿业大学,3,1 煤矿安全开采与岩性地震勘探,提高地质成果的勘探精度、实现安全开采是我国煤炭行业今后的主要奋斗目标。影响煤矿安全开采的地质因素很多，如构造、水文、瓦斯、煤层顶底板条件等。目前，主要成熟的勘探手段之一是三维地震勘探，但也仅限于解决构造问题。而煤矿安全开采中的主要地质问题，包括煤层顶底板水文地质条件、瓦斯突出条件与力学性质均属岩

2、性勘探范畴。,2020/8/23,中国矿业大学,4,1 煤矿安全开采与岩性地震勘探,构造地震勘探主要利用地震波的运动学特征，即利用地震波的旅行时和波速，计算出地层分界面上各点的埋藏深度，从而确定出地层的构造形态。岩性地震勘探除了利用地震波的运动学特征外，还利用地震波的动力学特征来研究地层的岩性。地震属性技术、相关/方差体技术、地震反演技术、地震相分析技术和谱分解技术是煤田岩性地震勘探的重要手段。,2020/8/23,中国矿业大学,5,2 地震属性技术2.1 地震属性,地震属性是指从叠前和叠后地震数据中提取出来的运动学、动力学和统计学地震特殊测量值。地震属性技术是指提取、显示、分析和评价地震属性

3、的技术，在煤田地震勘探中包括地震属性的提取、地震属性的分析、利用地震属性区分构造、岩性并进行目的层预测。 地震属性的分类没有统一的标准，结合煤田地震勘探的特点，可以根据运动学/动力学特征把地震属性分成八个类别：时间、振幅、频率、相位、波形、相关、吸收衰减、速度。地震属性的类型很多，要根据解决的地质问题来选择相应的地震属性。,2020/8/23,中国矿业大学,6,2.1 地震属性,地震属性技术的关键在于属性提取，提取方式包括同相轴属性提取和数据体属性提取。同相轴属性是与某个界面有关的地震属性，具体提取方法包括瞬时提取法、单道分时窗提取法和多道分时窗提取法。基于数据体的地震属性将产生一个完整的属性

4、体，其最大优点是能产生相关型的数据，从而提供逐道之间地震信号相似性和连续性的有用信息。将固定的三维数据体转化为能反映一定地球物理特征的新三维数据体。最常见的是相干数据体和方差数据体。,2020/8/23,中国矿业大学,7,2.2 应用实例1构造解释,阳煤集团新景煤矿保安区东部完成了三维地震勘探工作，对3煤层、8煤层和15煤层进行了构造解释。图1是该区的典型地震剖面，可以看出地震剖面的信噪比较高，对应三个煤层的T3波、T8波和T15波同相轴连续，可以连续追踪。 图2是3煤层构造解释方案，3煤层的构造相对比较简单，共解释断层18条，陷落柱8个。,2020/8/23,中国矿业大学,8,图1 典型地震

5、剖面,2020/8/23,中国矿业大学,9,图2 3煤层构造解释方案,2020/8/23,中国矿业大学,10,提取了T3波的多个地震属性，根据本区煤层反射波的波形特征，选择时窗长度为25ms；主频为45Hz；主频带为4050Hz。图3是T3波的地震属性，分别是振幅、方差、相似性、主频带能量和相位。 从属性切片中均可发现，在勘探区的东北部有一个构造异常区域(用黑色圆圈表示)，在原构造方案中均未进行解释。图4为异常区域内的地震剖面。结合地质资料进行综合分析，该构造异常区域为煤层变薄冲刷带。,2.2 应用实例1构造解释,2020/8/23,中国矿业大学,11,振 幅,图3 T3波地震属性,方 差,2

6、020/8/23,中国矿业大学,12,图3 T3波地震属性,主 频 带 能 量,相 位,2020/8/23,中国矿业大学,13,图4 3煤层构造异常区域的地震剖面,2020/8/23,中国矿业大学,14,图5是8煤层构造解释方案，8煤层的构造相对比较简单，共解释断层16条，陷落柱8个。 图6是T8波的地震属性，分别是振幅、相似性、主频带能量和相位。可以发现，所有地震属性在勘探区北部和南部差异明显，而原构造方案中均未进行解释。结合地质资料进行综合分析，勘探区南部为煤层变薄区域(0.8m)。,2.2 应用实例1构造解释,2020/8/23,中国矿业大学,15,图5 8煤层构造解释方案,2020/8

7、/23,中国矿业大学,16,振 幅,图6 T8波地震属性,相 似 性,2020/8/23,中国矿业大学,17,图6 T8波地震属性,主 频 带 能 量,相 位,2020/8/23,中国矿业大学,18,新景煤矿芦南二区北三正、副巷在掘进过程中发生瓦斯突出，7202综采工作面在回采过程中也发生瓦斯突出，瓦斯突出点(红星)位置见图7。图7中的红线范围是利用三维地震资料解释的3煤层变薄冲刷缺失带。北三正、副巷瓦斯突出点正在这个范围内，7202综采工作面瓦斯突出点在这个范围的边缘。显然，瓦斯突出与煤层的变薄缺失密切相关。,2.4应用实例2岩性解释,2020/8/23,中国矿业大学,19,图7 3煤层瓦斯

8、突出点位置,2020/8/23,中国矿业大学,20,2.4应用实例2岩性解释,图8是北三正巷的地震剖面，图9是北三副巷的地震剖面，图10是通过7202综采工作面瓦斯突出点的地震剖面，黄线代表T3波，瓦斯突出区域用红色箭头表示。从地震剖面中可以看出，对应3煤层的T3波有两个相位，且信噪比较高。原解释方案根据第一相位的振幅变弱，相位转移等波形特征确定了3煤层的变薄缺失带范围是基本可信的。,2020/8/23,中国矿业大学,21,图8 北三副巷地震剖面,2020/8/23,中国矿业大学,22,图9 北三正巷地震剖面,2020/8/23,中国矿业大学,23,图10 7202综采工作面瓦斯突出点地震剖面

9、,2020/8/23,中国矿业大学,24,2.4应用实例2岩性解释,提取了T3波的多个地震属性，根据本区煤层反射波的波形特征，选择时窗长度为17ms；主频为65Hz；主频带为6070Hz。图11是T3波的地震属性，分别是振幅、相似性、主频、主频带能量和平均频率相位。 从T3波地震属性图中可以看出3煤层的变薄缺失情况，在振幅切片上能量明显变弱(红色)，在相似性切片上异常突出(红色)，在主频切片上表现为高值(黑色)，在主频带能量切片上表现为低值(蓝色)，在相位切片上也存在较大差异。,2020/8/23,中国矿业大学,25,振 幅,图11 T3波地震属性,主 频 带 能 量,相 似 性,相 位,主频

10、,2020/8/23,中国矿业大学,26,2.4应用实例2岩性解释,相对于常规地震解释方法，利用地震属性能够更准确、更细致地划分3煤层的变薄缺失带。在这个范围内，煤层是部分变薄、部分缺失，北三正、副巷瓦斯突出点(图11中异常)的所有地震属性值接近煤层，与7202综采工作面瓦斯突出点的地震属性值基本一致。 基于上述事实，本区的瓦斯富集带位于煤层变薄缺失带附近，更准确地说应该是3煤层的变薄可能会导致瓦斯突出。因此，准确划分3煤层的变薄缺失带在就显得十分重要。,2020/8/23,中国矿业大学,27,2.5 结论,地震属性技术的主要特点是利用更多的地震信息(主要是动力学信息)，并且借助计算机的高分辨

11、能力对地震属性参数定量化，进行综合分析，避免了人为因素的影响。,2020/8/23,中国矿业大学,28,3 地震反演技术3.1 地震反演技术,地震反演利用地表观测地震资料，以已知地质规律和钻井、测井资料为约束，对地下岩层空间结构和物理性质进行成像(求解)的过程，是反演地层波阻抗(或速度)的地震特殊处理解释技术。地震反演方法具有明确的物理意义，是预测岩性的确定性方法，在油气勘探中取得了显著的地质效果。 煤田地震反演工作起步较晚，处在叠后地震反演的研究和初步应用阶段。近年来，我们把地震反演技术应用于多家煤矿，其关注的重点是煤矿安全开采的有关地质问题，获得了丰富的地质成果，主要包括：,2020/8/

12、23,中国矿业大学,29,3.1 地震反演技术,(1) 利用反演剖面提供的岩性信息来划分地层，研究煤层顶板的稳定性； (2) 确定奥陶系灰岩顶部岩层中的含隔水性，查明含、隔水层的空间分布和厚度分布； (3) 圈定岩浆岩侵入煤层的范围； (4) 预测煤层厚度； (5) 预测煤层瓦斯富集带； (6) 划分新生界下部地层岩性和含隔水性。,2020/8/23,中国矿业大学,30,3.2 应用实例1确定奥陶系灰岩顶部含、隔水层的空间分布和厚度分布,在临沂矿业集团新驿煤矿下组煤水文地质补充勘探工作中，利用地震反演技术查明部分水文地质条件，主要包括奥陶系顶部含、隔水层的空间分布和厚度分布。研究成果突破了传统

13、煤田地质学理论，表明奥陶系灰岩不完全是一个含水层，其顶部就有隔水层存在，可以为建立矿井突水的水量预测模型提供基础资料。,2020/8/23,中国矿业大学,31,3.2 应用实例1确定奥陶系灰岩顶部含、隔水层的空间分布和厚度分布,在常规地震剖面上，无法对奥陶系内部进行分层，原因是信噪比低和反射系数小(相对煤层而言)。但是利用地震反演剖面，借助相对波阻抗值的差异，便能够确定奥陶系顶部灰岩中的含、隔水层的空间分布和厚度分布。图12是水1群2联井反演地震剖面，从中可以看出奥陶系顶界面、隔水层和含水层的分层位置。图13是隔水层底界面的沿层波阻抗切片；图14是含水层底界面的沿层波阻抗切片。沿层波阻抗切片能

14、够提供整个研究区域内的含、隔水层岩性信息，隔水层的平均波阻抗明显高于含水层的平均波阻抗，利用该信息能够比较准确地确定奥陶系灰岩中的含、隔水层的空间分布。,2020/8/23,中国矿业大学,32,图12 联井反演地震剖面,16煤,奥陶系顶界面,隔水层底界面,含水层底界面,2020/8/23,中国矿业大学,33,图13 隔水层底界面沿层波阻抗切片,2020/8/23,中国矿业大学,34,图14 含水层底界面沿层波阻抗切片,2020/8/23,中国矿业大学,35,3.2 应用实例1确定奥陶系灰岩顶部含、隔水层的空间分布和厚度分布,奥陶系顶部隔水层时间间隔为15ms，换算为厚度30m；奥陶系顶部含水层

15、时间间隔为5ms，换算为厚度10m；二者之间的过渡带厚度是变化的。利用上述信息，可以确定奥陶系顶部的含、隔水层的空间分布和厚度分布。,2020/8/23,中国矿业大学,36,3.3 应用实例2预测煤层厚度,彬长煤田大佛寺井田含煤4层，其中4煤层可采，厚度变化范围为1.9317.04m，平均12.37m，属特厚煤层。在常规地震剖面上，4煤层对应的T4波振幅主要是由煤层顶、底板的反射系数决定的，无法辨别煤层厚度的变化。而反演地震剖面反映了地层的波阻抗变化，煤层在煤系地层中特征明显，其波阻抗值很低，与围岩存在巨大的阻抗值差异，因此能够比较准确地分辨出4煤层厚度的变化。,2020/8/23,中国矿业大

16、学,37,3.3 应用实例2预测煤层厚度,图15为D18井与D19井的联井反演地震剖面，可以清楚地看出4煤层厚度与波阻抗的变化关系。表1为本区9个钻孔处4煤层的预测厚度与实际揭露厚度对比结果。通过表1可以看出，4煤层预测的厚度与实际厚度基本吻合，根据该方法预测了全区4煤层厚度。,2020/8/23,中国矿业大学,38,图15 联井反演地震剖面,4煤,2020/8/23,中国矿业大学,39,表1 4煤层实际厚度与预测厚度对比表,2020/8/23,中国矿业大学,40,3.4 结论,地震反演技术是岩性地震勘探的重要手段之一，是一门集地震、测井、地质、计算机等多学科的综合地球物理勘探技术。根据钻孔测

17、井数据纵向分辨率很高的有利条件，对井旁地震资料进行约束反演，并在此基础上对孔间地震资料进行反演，推断煤系地层岩性在平面上的变化情况，这样就把具有高纵向分辨率的已知测井资料与连续观测的地震资料联系起来，实行优势互补，大大提高三维地震资料的纵、横向分辨率和对地下地质情况的勘探研究程度。,2020/8/23,中国矿业大学,41,4 地震相分析技术4.1 地震相分析技术,在进行煤系地层岩性解释过程中，普遍采用波阻抗反演和地震属性技术。随着煤矿开发对地层岩性信息的要求不断增加，这两种技术在某些程度上无法满足实际生产的需要，因为它们均无法评价地震信号的总体变化程度。 任何与地震波传播有关的物理参数变化都反

18、映在地震道波形的变化中，可以使用样点值随时间的变化来刻画和衡量地震道波形变化，这就是基于波形的地震相分析技术。基于波形的地震相分析技术综合利用了地震波的频率、相位、速度、能量等各种信息，是基于地震信号整体差异的分类，克服了上述缺陷，具有独特解决问题的能力。,2020/8/23,中国矿业大学,42,4.2 应用实例确定煤层中冲刷带范围,恒源煤矿62采区三维地震资料动态解释过程中，在6煤层发现了两个冲刷变薄带。时间剖面上的反射波形特征明显，主要表现为振幅削弱、频率增高、相位增多、反射波连性变差等特点；如果煤层缺失，则表现为反射波的消失，见图16中的红圈。,中国矿业大学,43,图16 冲刷变薄带在地

19、震剖面上的显示,2020/8/23,中国矿业大学,44,4.2 应用实例确定煤层中冲刷带范围,鉴于正常煤层和冲刷变薄带的波形特征存在明显差异，可以借助地震属性和地震相分析技术确定冲刷变薄带的范围。沿6煤层位时间向上、向下各10ms创建一个21ms等厚层段，在这个层段上提取多个地震属性并划分了5类地震相，充分考虑了地震波的振幅、频率、相位和波形特征。图17是6煤层的地震属性切片和分类数为5的地震相分布图，利用这些地震属性能够准确确定6煤层的冲刷变薄带范围。,2020/8/23,中国矿业大学,45,能 量,图17 6煤层的地震属性切片,5 类 地 震 相 分 布,相 似 性,相 位,2020/8/

20、23,中国矿业大学,46,4.3 结论,基于波形的地震相分析技术综合利用了地震波的频率、相位、速度、能量等各种信息，即基于地震信号的整体差异进行分类，因此可以作为弥补常规地震属性技术无法把握地震信号总体变化的一种手段。,2020/8/23,中国矿业大学,47,5 谱分解技术5.1 谱分解技术,谱分解技术的理论基础是薄层反射系统可产生复杂的谐振反射，利用薄层反射的频率属性可指示时间厚度变化。谱分解技术就是利用薄层调谐体离散频率特性，通过分析复杂岩层内频谱变化和局部相位的不稳定性，识别薄地层的横向分布特征。,2020/8/23,中国矿业大学,48,5.2 应用实例圈定岩浆岩侵入煤层的范围,皖北煤电

21、集团公司卧龙湖煤矿岩浆岩侵入现象严重，利用利用谱分解技术和地震相分析技术解释北一采区6煤层的岩浆岩侵入，确定了6煤层的岩浆岩侵入边界。 对于6煤层反射波(T6波)，正常煤层与岩浆岩侵入煤层之间存在着明显的频率差异，见图18。图18(a)是正常煤层T6波的频谱，其主频为65Hz；图18(b)是岩浆岩侵入T6波的频谱，其主频为89Hz。通过对比可知，当煤层被岩浆岩侵入后其反射波的主频增高，频谱宽度明显增加。因此，对T6波进行谱分解处理，选择低频、中频和高频谱分解切片作为解释的依据，见图19。,中国矿业大学,49,图18 T6波频谱特征,浆岩侵入反射波频谱征,正常煤层反射波频谱,2020/8/23,

22、中国矿业大学,50,25 Hz,图19 T6波谱分解切片,65 Hz,45 Hz,85 Hz,2020/8/23,中国矿业大学,51,5.2 应用实例圈定岩浆岩侵入煤层的范围,通过对比T6波的谱分解处理结果，发现岩浆岩侵入6煤层在谱分解切片上有以下主要特征： (1) 在低频段谱分解切片上表现为低能量，随着频率的提高，其能量越来越强。形成这一现象的主要原因在于煤层被岩浆岩侵入后，残存的煤层较薄，反射波能量主要集中在高频段。在25Hz切片上其能量极低；在45Hz切片上其能量有所提高，但明显低于正常煤层的能量；在65Hz切片上其能量与正常煤层的能量同样强；在85Hz切片上其能量略强于正常煤层的能量。,2020/8/23,中国矿业大学,52,5.2 应用实例圈定岩浆岩侵入煤层的范围,(2) 在所有频段的谱分解切片上都有一定的