塌陷采空区的数值模拟及属性解释

图２ 模型响应的时间剖面图
间剖面外，还分析了塌陷采空区时间剖面上的频谱 特征。通常认为频谱分析中包括振幅、相位以及频率 等地震动力学特征，而多数地震动力学特征是构造、 地层、岩性与油气或煤层等综合因素的反映。对煤层 以及塌陷采空区延迟反射波进行频谱分析见图3，4： 煤层记录道的主频为54 Hz，平均振幅为0.001 8，塌陷 采空区延迟反射波的主频变为48 Hz，平均振幅为 0.000 4。对比分析发现塌陷采空区平均振幅明显减 弱，主频也减小。 ３ 实例分析
山西焦煤科技
２０１４ 年第 １２ 期
图１ 塌陷采空区地质模型图
表１ 模型参数表
地层
纵波速度 ／ｍ／ｓ
横波速度 ／ｍ／ｓ
砂岩
３ ６００
２ ０７８
裂隙带压实区 ２ ４５０／３ ２７０～３ ６００
／离层区
１ ７３２
冒落带压实区 ２ ３１０／２ ６４０～３ ２７０
／离层区
１ ６１６
煤层
２ ２００
１ ２７０
灰岩
通过沉积岩的地球物理特征可知：设定砂岩、裂 隙带、冒落带、煤层以及灰岩的纵波速度分别为3 600 m/s、2 450/3 270～3 600 m/s、2 310/2 640～3 270 m/s、2 200 m/s、4 500 m/s。其中，裂隙带与冒落带的地层速度采 用梯度变化，且符合采动裂隙“O”形圈理论。品质因 子是描述岩石弹性的重要参数，描述了介质对能量的 吸收与散射，并依据经验公式Q=14VP2.2获得品质因子 Q，品质因子依次为235、156、134、79以及380，其余参 数见表1。 １．３ 地球物理特征
关键ຫໍສະໝຸດ Baidu 塌陷采空区；数值模拟；频谱分析；地震属性 中国分类号：TD166 文献标识码：A 文章编号：1672-0652（2014）12-0011-04
近年来，随着国家对煤矿整合、兼并力度的不断 加大，塌陷采空区给国有煤矿的安全生产留下了隐 患，直接影响着矿井机械化采煤工作面的布置，更是 矿井透水事故的主要来源。因此，提前查明老窑采空 区的位置是各大型矿井亟需解决的问题。
三维地震探测采空区主要依据的是不同介质间 的波阻抗差异。波阻抗差异越大，反射波的能量越强，
收稿日期：２０１４－１０－２８ 作者简介：徐书恩（１９６９—），男，山西五台人，１９９１年毕业于山东矿业学院，工程师，主要从事地测防治水技术管理工作
（Ｅ－ｍａｉｌ）ｘｕｓｈｕｅｎ１０９＠１６３．ｃｏｍ
· １２ ·
本文根据采空区的实际赋存形态建立了地质模 型，在时间剖面上较系统地分析了塌陷采空区的形态 特征，并分析了时间剖面上反射波的频谱特征。结合 一些典型的工程实例，认为煤矿塌陷采空区在时间剖 面上的形态特征以及地震属性特征均比较明显，可为 煤矿生产与安全提供指导意义。 １ 地震地质模型 １．１ 模型建立
当煤层被开采或被采空时，上覆地层的原有地应 力平衡遭到破坏。刘天泉院士等对我国煤矿开采覆 岩破坏与导水裂隙分布作了大量的实测和理论研 究[1-3]，对采场上覆岩层移动破断与采动裂隙分布规 律提出了“横三区”、“竖三带”的总体认识，即沿工作 面推进方向覆岩将分别经历煤壁支承影响区、离层 区以及重新压实区；由下往上岩层移动分为垮落带、 断裂带、整体弯曲下沉带，采空区中部的离层裂隙基
２０１４ 年第 １２ 期
Numerical Simulation and Attribute Analysis of Collapse Goaf
Xu Shu-en，Guo Zi-hua Abstract Collapse goaf in old coal mine seriously threatens the production and safety of large scale mine ， detecting the collapse goaf has become an important research topic in current 3D earthquake. Geological model is established according to the actual occurrence of the goaf，using ray tracing theory carries out numerical simulation. Systematically analyzes the morphological characteristics of collapse goaf and analyzes the spectrum characteristics of reflected wave on the time section. Combines with some typical engineering examples，sums up the seismic response characteristics and seismic attributes characteristics of the collapse goaf in the seismic data ，it provides a guiding significance for coal mine safety production. Key words Collapse goaf；Numerical simulation；Spectrum analysis；Seismic attributes
本被压实，而在采空区四周存在连通的离层裂隙发育 区，其形状与老顶岩板破断的“O-X”型相似，称之为 采动裂隙“O”形圈［4－6］。由于“竖三带”中弯曲带只产生 弯曲沉降，并不产生破裂，带内的岩层仍为层状结构， 因此，引起的物性差异普遍比冒落带以及裂隙带微弱 的多，此带可忽略不计。为此建立如下地质模型，见图 1，共建立了5层地层，分别为砂岩、冒落带、裂缝带、煤 层以及灰岩。 １．２ 地层物理参数
图７ 优势频率属性图
４结论 综上所述，时间剖面上的塌陷采空区内会产生许
多不连续的界面反射波，类似于“蚯蚓状”排列，波形 杂乱，塌陷采空区的煤层底板同相轴也会出现明显的 时间延迟现象，振幅减弱且频率降低；在地震属性解 释上其特征也比较明显，平均振幅以及优势频率也明 显地降低。上述特征为地震勘探方法在塌陷采空区 上的成功运用提供了现实与指导意义。
· １４ ·
山西焦煤科技
［11］ 程建远，孙红星，赵庆彪，等.老窑采空区的探测技术与实例研究［J］.煤炭学报，2008，33（3）:251-255. ［12］ 陈相府，安西峰，王高伟.浅层高分辨地震勘探在采空区勘测中的应用［J］.地球物理学进展，2005，20（2）:381-386. ［13］ 杨双安，宁书年.老窑采空区的地震探测与研究［J］.中国煤田地质，2004，16（1）:44-47. ［14］ 张延玲，杨长春，贾曙光.地震属性技术的研究和应用［J］.地球物理学进展，2005，20（4）:1129-1133. ［15］ 魏 艳，尹 成，丁 峰，等.地震多属性综合分析的应用研究［J］.石油物探，2007，46（1）:42-47.
本次模型是采用射线追踪理论进行正演模拟。模 型宽度710 m，深度600 m，选用自激自收的放炮方 式，雷克子波主频为60 Hz，道间距10 m。
依此获得模型的响应时间（图2）：煤层反射波同 相轴下凹、中断，以中断为主，且同相轴变粗、频率变 低，在塌陷采空区上部同相轴发生抖动、连续性变差、 “蚯蚓状”排列，波形杂乱。除了分析塌陷采空区的时
参考文献 ［1］ 刘天泉.矿山岩体采动影响与控制工程学及其应用［J］.煤炭学报，1995，20（1）:1-5. ［2］ 钱鸣高，刘听成.矿山压力及其控制［M］.北京:煤炭工业出版社，1991：99-100. ［3］ 国家煤炭工业局制定.建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程［M］.北京:煤炭工业出版社，2000：157-159. ［4］ 钱鸣高，许家林.覆岩采动裂隙分布的“O-X”形圈特征研究［J］.煤炭学报，1998，2003（5）：466-469. ［5］ 李 智，王汉鹏，李术才，等.煤层开采过程中上覆岩层裂隙演化规律研究［J］.山东大学报，2011，41（3）:142-147. ［6］ 姜永东，鲜学福，李晓红，等.采场上覆岩层塑性断裂区域和移动规律研究［J］.中国矿业，2009，18（3）:67-71. ［7］ 裔传标，朱书阶，孙永亮，等.三维地震勘探技术在老窑采空区探测中的应用［J］.煤炭技术，2008，27（4）:88-91. ［8］ 裴文春，王德民，程增庆，等.三维地震资料解释技术分析煤层冲刷及采空区［J］.煤炭科学技术，2007，35（8）:32-34. ［9］ 王 晶，张建州，范庆荣，等.综合物探技术在整合煤矿复杂采空区勘探中的应用［J］.煤矿开采，2012，17（4）:21-24. ［10］ 刘 岩.采空区上的地震勘探方法［J］.煤炭技术，2008，27（4）:100-101.
４ ５００
２ ５９８
密度 ／ｋｇ／ｍ３ ２ ４０１
２ ２９４
２ ２５５ ２ １２３ ２ ５３９
品质 因子 ２３５ １５６
１３４ ７９ ３８０
勘探效果越好。煤层未开采时，呈现低密度、低频、低 速度的特征，与周围围岩存在着明显的波阻抗，能形 成能量较强的反射波［7-9］。当煤层被开采形成采空区 时，上覆地层变得疏松，密度降低，传播于其中的地震 波的能量、速度、频率均发生不同程度的变化。在时间 剖面上主要表现为煤层反射波同相轴及辅助相位 不 连 续 、波 形 异 常 、无 规 律 分 布 、频 率 较 低 、延 时 现 象明显 。 ［10-13］ ２ 模型正演
工区位于晋西黄土高原的最南端，属强烈侵蚀 的低、中山地貌，勘探区地形复杂，相对高差120 m。区 内坡度较陡，黄土大面积覆盖，浅表层地球物理条件 较差。2#煤位于山西组下部，埋深250～430 m，平均厚 度6 m，其下50 m左右为10#煤层，厚度3.5 m，与其顶 底板岩层形成明显的波阻抗，故会形成两个良好的 反射界面，适合三维地震勘探。
第 １２ 期 ２０１４ 年 １２ 月
·试验研究·
山西焦煤科技 Ｓｈａｎｘｉ Ｃｏｋｉｎｇ Ｃｏａｌ Ｓｃｉｅｎｃｅ ＆ Ｔｅｃｈｎｏｌｏｇｙ
Ｎｏ．１２ Ｄｅｃ．２０１４
塌陷采空区的数值模拟及属性解释
徐书恩，郭子华
（山西煤炭运销集团，山西 太原 030013）
摘 要 老窑塌陷采空区严重威胁着大型矿井的生产与安全，探测出塌陷采空区已成为当前三 维地震的一项重要研究课题。根据采空区的实际赋存形态建立了地质模型，并采用射线追踪理论进 行数值模拟，在时间剖面上较系统地分析了塌陷采空区的形态特征，并分析了反射波的频谱特征。结 合一些典型的工程实例，总结出塌陷采空区在地震资料上的地震响应特征以及地震属性特征，对煤 矿安全生产具有指导意义。
图３ 煤层记录道振幅谱图 图４ 塌陷采空区记录道振幅谱图
2014 年第 12 期
徐书恩等：塌陷采空区的数值模拟及属性解释
本次地震勘探新发现采空区: 采空区I （CK1）、采 空 区 II（CK2）、采 空 区 III（CK1）、采 空 区 IV（CK1）共 4 块，仅以其中某个采空区为例，从地震时间剖面上看， 采空区均表现为T2 波同相轴中断，对应区域煤层反 射波能量变的很弱；波形异常，辅助相位排列杂乱、 无规律、断断续续；频率变低，有绕射波、多次波等多 种波混杂叠加，见图5。
· １３ ·
图６ 平均振幅属性图
图５ 采空区在时间剖面上反映示意图
根据采空区发育地区的地震属性可知，利用沿煤 储层的层间地震统计特征更有利于对煤层采空区等 异常地质体的识别。地震属性可划分为振幅统计类、 频（能）谱统计类、相位统计类、复地震道类、层序统计 类以及相关统计类6大类［14］。振幅类属性反映波阻抗 差异、地层厚度、孔隙度及含流体成分变化；频谱类属 性特征反映储藏岩性变化、裂缝分布范围、油气水分 布变化和储藏上覆地层对地震波的衰减 ［15］。经过对 比分析，分别选用了平均振幅以及优势频率属性。从 地震属性分析：平均振幅强度在采空区内明显减弱， 由175减弱到75；优势频率也明显减小，由20 000减小 到5 000，见图6，7。