薄互层条件下围岩变化对煤层反射波的影响研究

常锁亮，张 生，刘 晶 ，刘最亮，陈 强，刘 波

薄互层条件下围岩变化对煤层反射波的影响研究

常锁亮1,2，张 生1,2，刘 晶1,2，刘最亮3，陈 强1,2，刘 波1,2

(1. 太原理工大学 地球科学与工程系，山西 太原 030024；2. 煤与煤系气地质山西省重点实验室，山西 太原 030024；3. 华阳新材料科技集团有限公司，山西 阳泉 045000 )

为研究薄互层条件下围岩变化对煤层反射波的影响，以地震勘探中/4薄层范围内含煤地层为研究对象，建立围岩岩性、厚度与结构变化三类模型，基于薄层反射系数谱理论中Brekhovskikh方程，计算并总结围岩变化对煤层AVO曲线、属性及道集的影响。研究结果表明：/4范围内围岩岩性变化会对煤层AVO响应产生显著影响，其中顶板为砂岩会使得煤层AVO截距和梯度属性明显增大，顶板为泥岩会使得煤层AVO截距和梯度属性增大，且顶板岩性不同，对应的煤层AVO道集特征也会发生变化；/4范围内围岩互层结构和厚度变化会对煤层AVO响应产生一定影响，但是影响较小，其中围岩互层结构的变化会使得煤层AVO道集特征产生变化，围岩厚度的变化会使得煤层AVO截距属性产生变化；基于界面型的Zoeppritz方程不适用于薄互层含煤地层的正演模拟，应选取更适用于薄互层的Brekhovskikh正演方程或者其他模拟方法。

薄互层；围岩变化；反射系数谱；煤层反射波

我国含煤盆地多形成于陆相或海陆过渡相的沉积环境[1-2]，含煤地层多以薄互层形式存在[3]，且岩性和厚度横向变化较大，使得煤层的反射波响应中既包含煤层与煤体结构的响应，也包含煤厚变化与围岩结构的影响。已有研究表明，对于单个薄煤层反射波来说，煤层厚度[4-6]、互层围岩的岩性[7-9]、互层结构[10-12]以及入射波主频[13-14]等都会对薄煤层的地震响应造成影响。

已有对煤层反射波的影响研究可归纳为两种：第一种是通过建立单一薄煤层模型讨论各类影响要素对薄煤层纵波与横波AVA特征的影响；张铁强等[6]基于反射率法从AVO振幅、截距、梯度以及频谱等方面区分了不同厚度煤层模型的特征与差异；Zou Guangui等[7]通过建立楔形薄煤层模型，研究了煤层厚度与地震波振幅的关系以及影响这种关系的因素；Chen Tiansheng等[8]基于反射系数谱理论讨论了单薄层的多波AVO响应；Pan Wenyong等[9]基于反射系数谱理论分析了薄层弹性介质中AVO效应与频率的关系；蔡希玲等[10]讨论了地震子波主频和厚度对薄层AVO振幅和频谱的影响。但是该类模型没有考虑含煤地层中煤层与围岩互层分布、且横向厚度不均一的实际情况，忽略了薄互层中围岩变化对煤层反射波造成的影响，模拟结果存在偏差。

第二种是通过建立薄互层含煤地层模型讨论整个含煤地层的反射波特征及影响因素。安莹等[11]基于各向异性与长波长理论，根据煤厚与含煤地层总厚的比值建立模型，研究含煤地层的多波AVO地震响应；王增玉等[12]基于反射系数谱理论分析了薄煤层中夹矸和构造煤对煤层AVO特征的影响。这种模型符合实际地层展布，但目前研究仅限于整个含煤地层的综合地震响应及影响因素，实际煤田地震勘探更关注含煤地层中单一煤层反射波所受影响，已有研究[15-18]对煤层自身厚度及结构变化对煤层反射波响应产生的影响进行了探讨，然而实际含煤地层为不同岩性的薄互层组成，煤层上下围岩岩性、结构以及厚度均呈现较强非均一性。

因此，针对薄互层含煤地层中围岩变化对煤层反射波造成的影响展开研究。分别建立围岩岩性、结构以及厚度变化的含煤地层模型，模拟并总结围岩变化对煤层反射波AVO曲线、属性以及道集特征造成的影响，同时研究正演模拟采用的薄层反射系数谱法与常规基于Zoeppritz方程法的差异，证明研究方法的合理性。

1 薄层反射系数谱原理

薄层反射系数谱理论是从弹性波动力学理论的角度讨论平面简谐波在层状介质中的反射和透射，其理论基础是Brekhovskikh方程。如果地层厚度小于等于一个波长，则这种地层上下界面的波动将会发生明显的干涉和叠加，使每一层面所形成的波动在强度变化及波形特点各方面都会具有与厚层不同的特点。对比图1中的单界面模型和交互薄层模型，可以看出交互薄层的反射和透射是由薄互层中经受不同次数反射后形成的单波叠加结果，每个单波之间的时差及其干涉条件取决于薄层的厚度、薄层中波的传播速度及入射波的入射角[17]。

如图1a所示，设平面波从顶层入射到下面各层中，纵横波到各层入射角为，则顶层中各波的位移函数表达式分别为：

图1 波的传播

Fig.1 Schematic diagram of wave propagation

位移-应力关系为：

反射系数和透射系数计算公式为：

如图1b所示，对于一系列薄层来说，来自叠加薄层的反射波可以被认为是来自每个薄层的反射波和透射波的复合波。因此，可以利用递推算法来计算薄层叠加层顶底之间的位移和应力关系，从而可以求取每一层的反射透射系数。

其中，

式中：k为水平方向慢度；s、d分别为纵波和SV波慢度；sz、dz分别为纵波和SV波垂向慢度；为地层厚度。

2 围岩对煤层反射波的影响

为分析围岩对煤层地震响应的影响，将含煤地层视为VTI型薄互层[15-18]，根据表1中的物性参数，分别构建围岩岩性、厚度与结构变化地层模型，基于Brekhovskikh方程计算不同模型各个层位反射系数，提取煤层反射波的AVO曲线[19-20]、梯度和截距属性[21-23]以及AVO道集[24-26]，讨论围岩变化对其造成的影响。

表1 模型参数

2.1 模型构建

由于含煤地层多以薄互层形式存在，为研究煤层薄层顶底板存在对不同厚度煤层顶板反射波是否存在影响，分别建立不包含煤层顶底板围岩的一层模型(图2a)和包含顶底板围岩交互的三层模型(图2b)进行对比，每类模型都包含1～8 m不同厚度的煤层模型，围岩为2 m厚的泥岩。

2.2 围岩对煤层反射波AVO曲线特征影响分析

基于反射系数谱方法选取40 Hz雷克子波以及表1中的模型参数计算煤层顶板反射系数，两类模型对应的AVO曲线如图3所示。由图3可以发现：不同模型的煤层顶板AVO曲线中，总体表现为反射系数值均为负值，且绝对值随着偏移距的增大先减小后增大的趋势(由于40°以后的数据超过了临界角所以不予考虑)。与一层模型(不含顶底板围岩)相比，可以发现三层模型(含顶底板模型)由于顶底板薄层围岩存在调谐效应，所以对煤层反射产生了不可忽略影响，当煤厚小于2 m时，顶底板存在使反射系数的绝对值减小，当煤厚大于2 m时，顶底板存在使反射系数的绝对值增大。所以三层结构的建模方式更加符合实际情况，不能忽略顶底板岩性存在对煤层反射波所造成的影响。

图3 煤厚变化单煤层及煤层与顶底板交互模型的AVO曲线

3 围岩变化对煤层反射波的影响

为分析围岩变化对煤层地震响应的影响，分别构建围岩岩性、结构以及厚度变化3类模型，讨论煤层反射波的AVO曲线、属性以及道集特征的变化。

3.1 围岩岩性变化影响

3.1.1 围岩岩性变化模型构建

为研究围岩岩性变化对煤层反射波的影响，选择顶、底板和砂、泥岩2项参数变化作为变量分别组成4组模型，模拟实际含煤地层中围岩岩性变化情况。如图4所示，模型1中煤层顶底板均为泥岩；模型2中煤层顶底板均为砂岩；模型3中煤层顶板为泥岩，底板为砂岩；模型4中煤层顶板为砂岩底板为泥岩。模型中煤厚为/12，围岩厚度为/12。

图4 围岩岩性变化模型

3.1.2 围岩岩性变化对煤层反射波AVO曲线影响特征分析

图5为4类模型的AVO曲线，由图5可以看出4类模型的煤层顶板AVO曲线总体表现为反射系数值均为负值，且绝对值随着偏移距的增大先减小后增大。由直接顶板的岩性不同可将AVO曲线分为两类：顶板为泥岩类和顶板为砂岩类。第一类模型反射系数变化范围为-0.625～-0.450；第二类反射系数变化范围为-0.7～-0.2。研究结果表明，煤层顶板围岩的岩性变化会对煤层反射波产生不可忽略的影响，且可由顶板岩性将煤层AVO曲线分为两大类：泥岩类和砂岩类。

图5 围岩岩性变化模型AVO曲线

3.1.3 围岩岩性变化对煤层反射波AVO属性特征影响分析

为进一步刻画不同模型的截距和梯度属性特征，对模型1—模型4进行梯度和截距值统计，得到截距、梯度交会图(图6)。顶板为砂岩类模型的梯度明显大于顶板为泥岩类模型，同时截距绝对值大于顶板为泥岩类模型，其交会点位于图中左上角，泥岩类模型属性交会点位于右下角。由此可区分出不同顶板岩性的两大类模型，说明围岩岩性变化对煤层反射波的AVO属性特征产生明显影响。

3.1.4 围岩岩性变化对煤层反射波的AVO道集特征影响分析

为研究围岩岩性变化对煤层道集特征的影响，分别对顶板为泥岩和砂岩两类模型计算各个反射层的反射系数，得到模型对应的反射系数谱，接着用主频为40 Hz雷克子波与该反射系数谱进行褶积运算，得到模型AVO道集，如图7所示，其中图7a为泥岩围岩模型道集，图7b为砂岩围岩模型道集。为突出道集特征，进一步提取地震道集的均方根振幅属性。从图中可以看出，围岩岩性变化会对含煤地层AVO道集产生明显影响，砂岩类振幅随入射角增大而减小的拐点角度小于泥岩类，前者为20°，后者为30°。

图6 围岩岩性变化模型AVO梯度–截距属性

3.2 围岩结构变化的影响

3.2.1 围岩结构变化模型构建

实际含煤地层中围岩并不是单一结构，而是由砂泥岩互层结构组成，为研究围岩结构变化对煤层反射波的影响，固定煤厚为/12，含煤地层总厚为/4，逐步增加围岩内部砂泥岩互层数建立模型，如图8所示，砂泥岩互层数从模型1到模型5依次从1层增加到5层。

图7 围岩岩性变化模型AVO道集

图8 围岩结构变化模型

3.2.2 围岩结构变化对煤层反射波的AVO曲线影响分析

围岩结构变化模型对应的AVO曲线如图9所示，由图可知不同围岩结构模型AVO 曲线形态基本一致。随着互层数由一层递增到三层，反射系数的模稍有减小，当互层数增加到三层以上时，不再产生影响，说明围岩结构变化对煤层反射波的AVO曲线影响较小，所以对其梯度和截距影响特征不再分析。

图9 围岩结构变化模型AVO曲线

3.2.3 围岩结构变化对煤层反射波的AVO道集影响分析

为研究含煤地层AVO道集是否受到围岩互层结构影响，分别对一层(图10a)、两层(图10b)和三层(图10c)结构的模型计算每个反射层的反射系数，得到模型对应的反射系数谱，接着用主频为40 Hz雷克子波与该反射系数谱进行褶积运算，得到不同围岩结构模型对应的AVO道集，为凸显道集特征进一步提取了地震道集的均方根振幅属性。从图10中可以看出随着互层结构的增多，煤层反射波振幅明显增强，且表现为三层强反射轴。以上研究表明围岩互层结构对整个含煤地层AVO道集产生显著影响。

3.3 围岩厚度变化影响

3.3.1 围岩厚度变化模型构建

为研究围岩厚度变化对煤层反射波的影响，首先固定模型中煤厚为/12，逐步增加围岩厚度，如图11所示，厚度依次由/120增加到/12，一共包含10组模型，所有模型的含煤地层总厚保持在/4(属地震勘探薄层)范围内，其中煤层直接顶底板为泥岩。

3.3.2 围岩厚度变化对煤层反射波AVO曲线影响分析

不同围岩厚度模型的AVO曲线如图12所示，可以看出随着围岩厚度增加，反射系数的绝对值逐渐增大，但是增大范围基本小于0.05，变化程度不大，曲线形态保持不变，说明围岩厚度变化对煤层反射波的AVO曲线具有一定程度影响，实际应用中不能忽略围岩厚度变化的影响。

3.3.3 围岩厚度变化对煤层反射波AVO属性影响分析

为进一步刻画AVO属性特征，对围岩厚度变化模型进行截距值统计，得到截距随围岩厚度变化的相关图，如图13所示，分别统计了5组入射角的截距–围岩厚度变化情况。从图中可以看出，不论入射角怎么变化，截距均随围岩厚度增加而变大，且变化范围在0～0.05。说明围岩厚度变化会对煤层反射波的AVO截距属性产生一定影响。

3.3.4 围岩厚度变化对煤层反射波的AVO道集影响分析

为研究含煤地层AVO道集是否受到围岩厚度影响，选取围岩厚度为/120和8/120两个厚度模型计算每个反射层的反射系数，得到模型对应的反射系数谱，接着用主频为40 Hz雷克子波与该反射系数谱进行褶积运算，得到模型对应的AVO道集，为凸显道集特征进一步提取均方根振幅属性，如图14所示。由图14可以看出，厚度对含煤地层AVO道集整体影响较小，随厚度增加，AVO道集出现两层强反射轴。说明围岩厚度变化对含煤地层的AVO道集影响较小。

图10 围岩结构变化模型AVO道集

图11 围岩厚度变化模型

图12 围岩厚度变化模型AVO曲线

图13 围岩厚度变化模型的AVO截距–厚度交会图

4 实例验证

为验证本文基于反射系数谱计算方法的合理性，通过阳泉某矿区叠前地震道集资料(图15a)提取XJ-6钻孔处3号煤层不同偏移距对应的振幅值，由钻孔信息可知，井点3号煤埋深为516 m，所以利用埋深和偏移距计算得到入射角，利用振幅值以及入射角拟合3号煤AVO曲线。同时以钻孔岩性及厚度分布情况建模，分别采用本文基于反射系数谱法中Brekhovskikh方程和常规Zoeppritz方程计算3号煤AVO曲线。可以看出，本次基于反射系数谱法计算结果更接近实际AVO拟合曲线，常规Zoeppritz方程模拟结果截距和梯度均较小。为进一步验证该结果，通过2种方法计算过研究区3-166井的AVO梯度属性，图16a和16b分别为Brekhovskikh方程和Zoeppritz方程计算的AVO梯度属性结果，从3号煤反射波梯度属性可以看出，Brekhovskikh方程计算梯度衰减结果要明显强于Zoeppritz方程计算结果，计算结果与过井处道集特征一致，进一步证明薄互层条件下围岩变化对煤层反射波的影响不同于常规储层，在实际资料处理与解释中需要进一步关注和研究。

图14 围岩厚度变化模型与AVO道集

图15 某矿区XJ-6井叠前AVO道集与AVO曲线

5 结论

a./4范围内围岩岩性变化会对煤层AVO响应产生显著影响，当顶板岩性为砂岩时，AVO截距和梯度属性明显增大，当顶板岩性为泥岩时，AVO截距和梯度属性增大，且当煤层顶板岩性变化时，对应的AVO道集特征也会发生变化。

b. λ/4范围内围岩互层结构和厚度变化对煤层AVO响应影响较小，其中随着围岩互层结构增多，煤层AVO道集特征产生变化，随着单层围岩厚度增加，煤层AVO截距属性略有增加。

c.基于界面型的Zoeppritz方程不适用于薄互层含煤地层的正演模拟，应选取更适用于薄互层的Brekhovskikh正演方程或者其他模拟方法。

图16 过研究区3-166井2种方法计算AVO梯度剖面

[1] KENNETT B L N. Seismic wave propagation in stratified media[J]. Geophysical Journal International，1986，86(1)： 219–220.

[2] KENNETT B L N. Seismic waves in a stratified half space[J]. Geophysical Journal of the Royal Astronomical Society，1980，61(1)：1–10.

[3] 师素珍，郭家成，谷剑英，等. 调谐效应影响下的AVO属性在构造煤预测中的应用[J]. 煤田地质与勘探，2018，46(3)：133–138.

SHI Suzhen，GUO Jiacheng，GU Jianying，et al. Application of tectonic coal prediction by AVO attributes under tuning effect[J]. Coal Geology & Exploration，2018，46(3)：133–138.

[4] YANG Chun，WANG Yun，WANG Yanghua. Reflection and transmission coefficients of a thin bed[J]. Geophysics，2016，81(5)：N31–N39.

[5] 杨震，芦俊，孟星浑，等. 薄煤层 PP 波与 PS 波 AVA 地震响应特征[J]. 煤炭学报，2015，40(6)：1435–1441.

YANG Zhen，LU Jun，MENG Xinghun，et al. PP-and PS-wave AVA response characteristics for thin coal seam[J]. Journal of China Coal Society，2015，40(6)：1435–1441.

[6] 张铁强，孙鹏远，钱忠平，等. 薄煤层 AVO 响应特征分析[J]. 石油地球物理勘探，2013，48(4)：597–603.

ZHANG Tieqiang，SUN Pengyuan，QIAN Zhongping，et al. AVO analysis on thin coalbed[J]. Oil Geophysical Prospecting，2013，48(4)：597–603.

[7] ZOU Guangui，XU Zhiliang，PENG Suping，et al. Analysis of coal seam thickness and seismic wave amplitude：A wedge model[J]. Journal of Applied Geophysics，2018，148：245–255.

[8] CHEN Tiansheng，LIU Yang. Multi-component AVO response of thin beds based on reflectance spectrum theory[J]. Applied Geophysics，2006，3(1)：27–36.

[9] PAN Wenyong，INNANEN K A. AVO/AVF analysis of thin beds in elastic media[J]. SEG Technical Program Expanded Abstracts，2013：373–377.

[10] 蔡希玲，刁文川，周兴元，等. 薄层反射波非零炮检距的属性特征[J]. 石油地球物理勘探，2007，42(3)：277–282.

CAI Xiling，DIAO Wenchuan，ZHOU Xingyuan，et al. Attribute features of reflection on non-zero offsets in thin layers[J]. Oil Geophysical Prospecting，2007，42(3)：277–282.

[11] 安莹，芦俊，杨春. 含煤地层多波AVO地震响应[J]. 煤炭学报，2018，43(3)：793–800.

AN Ying，LU Jun，YANG Chun. Multi-wave AVO seismic response of coal bearing strata[J]. Journal of China Coal Society，2018，43(3)：793–800.

[12] 王增玉，杨德义，曹志勇，等. 构造煤及夹矸对煤层AVO正演模拟结果影响分析[J]. 地球物理学进展，2018，33(2)：754–759.

WANG Zengyu，YANG Deyi，CAO Zhiyong，et al. Analysis on the influence of tectonic coal and parting on AVO forward modeling of coal seam[J]. Progress in Geophysics，2018，33(2)：754–759.

[13] LIU Yinbin，SCHMITT D R. Amplitude and AVO responses of a single thin bed[J]. Geophysics，2003，68(4)：1126–1422.

[14] SWAN H. Amplitude-versus-offset measurement errors in a finely layered medium[J]. Geophysics，1991，56(1)：9–159.

[15] 赵晨光. 薄互层地震反射波的特征分析[J]. 石油地球物理勘探，1986，21(1)：32–46.

ZHAO Chenguang. The analysis of seismic reflection from thin interbedding[J]. Oil Geophysical Prospeting，1986，21(1)：32–46.

[16] BACKUS G E．Long-wave elastic anisotropy produced by horizontal layering[J]．Journal of Geophysical Research，1962，67(11)：4427–4440．

[17] 吴飞勇，马培元，黄明亮. 薄(互)储层的反射系数谱AVO分析[J]. 地球物理学进展，2010，25(1)：71–75.

WU Feiyong，MA Peiyuan，HUANG Mingliang. Thin reservoir AVO analysis of spectral reflection coefficient[J]. Progress in Geophysics，2010，25(1)：71–75.

[18] POSTMA G W.Wave propagation in a stratified medium[J]. Geophysics，1955，20(4)：745–961．

[19] SAYERS C. Long-wave seismic anisotropy of heterogeneous reservoirs[J]. Geophysical Journal International，1998，132(3)：667–673．

[20] SIDLER R，HOLLIGER K. Seismic reflectivity of the sediment-covered seafloor：Effects of velocity gradients and fine-scale layering[J].Geophysical Journal International，2010，181(1)：521–531．

[21] 汪恩华，贺振华，李庆忠. 基于薄层的反射系数谱理论与模型正演[J]. 成都理工学院学报，2001，28(1)：70–74.

WANG Enhua，HE Zhenhua，LI Qingzhong. Reflectance spectrum theory and model computation based on thin beds[J]. Journal of Chengdu University of Technology，2001，28(1)：70–74.

[22] 李雪英，陈树民，王建民，等. 薄层时频特征的正演模拟[J]. 地球物理学报，2012，55(10)：3410–3419.

LI Xueying，CHEN Shumin，WANG Jianmin，et al. Forward modeling studies on the time-frequency characteristics of thin layers[J]. Chinese Journal of Geophysics(in Chinese)，2012，55(10)：3410 –3419.

[23] 徐果明，李跃，倪四道，等. 薄互层等效横向各向同性的研究[J]. 石油地球物理勘探，1996，31(6)：792–805.

XU Guoming，LI Yue，NI Sidao，et al. Research on equivalent azimuthal isotropy of thin interbed medium[J]. Oil Geophysical Prospecting，1996，31(6)：792–805.

[24] ZHANG Sheng，HUANG Handong，ZHU Baoheng，et al. Seismic facies-controlled prestack simultaneous inversion of elastic and petrophysical parameters for favourable reservoir prediction[J]. Exploration Geophysics，2017，49(5)：655–668.

[25] ZHANG Sheng，HUANG Handong，LUO Yaneng. Study on qp-wave elastic impedance characteristic in fractured equivalent TTI media[C]//79thEAGE Conference & Exhibition 2017. Pairs，France，June 12-15，2017.

[26] ZHANG Sheng，HUANG Handong，DONG Yinping，et al. Direct estimation of the fluid properties and brittleness via elastic impedance inversion for predicting sweet spots and the fracturing area in the unconventional reservoir[J]. Journal of Natural Gas Science & Engineering，2017，45：415–427.

Influence of surrounding rock changes on the coal seam reflected wave under thin interbed condition

CHANG Suoliang1,2, ZHANG Sheng1,2, LIU Jing1,2, LIU Zuiliang3, CHEN Qiang1,2, LIU Bo1,2

(1. Department of Geosciences and Engineering, Taiyuan University of Technology, Taiyuan 030024, China; 2. Key Lab of Geology of Coal and Coal Measure Gas, Taiyuan 030024, China; 3. Huayang New Material Technology Group Co., Ltd.,Yangquan 045000, China)

In order to study the influence of the changes of surrounding rock on the reflected wave of coal seam under the condition of thin interbeds, the coal-bearing strata in the/4 thin layer range in seismic exploration is taken as the research object, and three types of models of surrounding rock lithology, thickness and structure change are established. By Brekhovskikh equation in the theory of thin-bed reflection coefficient spectrum, the influence of surrounding rock changes on the AVO curve, attributes and gathers of coal seams are calculated and summarized. The research results show that the changes of surrounding rock lithology in the range of/4 have a significant impact on the AVO response of the coal seam. The sandstone roof will significantly increase the AVO intercept and gradient properties of the coal seam. The mudstone roof makes the AVO intercept and gradient properties of the coal seam increase; because of the different roof lithology, the corresponding coal seam AVO gather characteristics will also change. The surrounding rock interbed structure and thickness changes in the range of/4 will have a certain impact on the coal seam AVO response, but the impact is small. Among them, the changes of the surrounding rock interbed structure will cause changes in the AVO gather characteristics of the coal seam, and the surrounding rock thickness will cause changes in the AVO intercept attribute of the coal seam. The Zoeppritz equation based on the interface type is not suitable for the normalization of the thin interbedded coal-bearing strata. For simulation, the Brekhovskikh forward equation or other simulation methods that are more suitable for thin interlayers should be selected.

thin interbed; surrounding rock change; reflection coefficient spectrum; coal seam reflection wave

移动阅读

语音讲解

P65；P631

A

1001-1986(2021)05-0220-10

2020-10-19；

2021-05-28

阳泉煤业(集团)有限责任公司研究开发项目(GY18027)；山西省面上青年基金项目(201901D211005)

常锁亮，1972年生，男，山西灵石人，博士，副教授，研究方向为煤与煤层气地震勘探. E-mail：changsuoliang@tyut.edu.cn

张生，1986年生，男，安徽合肥人，博士，讲师，研究方向为非常规油气地震预测.E-mail：zhangsheng@tyut.edu.cn

常锁亮，张生，刘晶，等. 薄互层条件下围岩变化对煤层反射波的影响研究[J]. 煤田地质与勘探，2021，49(5)：220–229. doi: 10.3969/j.issn.1001-1986.2021.05.024

CHANG Suoliang，ZHANG Sheng，LIU Jing，et al. Influence of surrounding rock changes on the coal seam reflected wave under thin interbed condition[J]. Coal Geology & Exploration，2021，49(5)：220–229. doi: 10.3969/j.issn. 1001-1986. 2021.05.024

(责任编辑 聂爱兰)