一地球物理特征

1表层地震地质条件

测区内沟谷纵横，多呈“V”形地貌，地形高差较大。地表出露的二叠系红色泥质砂岩地层常形成较陡的滑坡，粉砂岩区多形成深沟或者悬崖峭壁以及崩塌堆积场，攀越困难。区内村庄和花椒树地较多，影响激发接收效果和地震测线的布置。表层地震地质条件较差。

2浅层地震地质条件

测区内基岩和黄土塬均有分布，基岩出露较多。基岩能够较好地传播地震波，能量损失较小，激发效果较好。黄土塬的黄土对地震波有较强的吸收散射和低通滤波作用，降低了地震波的能量和频率，激发效果不佳。因此浅层地震地质条件较好。

3深层地震地质条件

测区内地层自下而上有中奥陶统下马家沟组、上马家沟组、峰峰组，石炭系中上统本溪组及太原组，二叠系下统山西组、下石盒子组，二叠系上统上石盒子组、石千峰组及第四系更新统和全新统。含煤地层为石炭—二叠系，煤层顶、底板为泥岩、砂质泥岩、粉砂岩等，其饱和抗压强度均小于30MPa，属于软弱岩石，二者之间波阻抗差异较大，具有形成反射波的良好条件。但测区局部的3煤层已经开采，可能对11煤层、奥灰的反射波形成存在较大影响。综上所述，本区表层地震地质条件较差，浅、深层地震地质条件相对较好。

二地震数据野外采集和处理

1地震数据野外采集

通过试验工作确定了地震数据野外采集参数。野外观测系统采用8线8炮制，单线96道中点放炮，道间距10m，线距40m，炮排距80m，CDP网格为5m（纵向）×10m（横向），叠加次数24次(横向4次，纵向6次)，最小非纵距10m；激发条件为单井激发，厚黄土覆盖区，钻孔打到粘土层中，药量3.0kg；薄黄土覆盖区，钻孔打到基岩面或粘土层，药量2.0kg；基岩出露区，井深4m，药量1.5kg；检波器采用60Hz二串二并堆放插置；仪器使用加拿大产ARIES遥测数字地震仪，记录的数据格式为SEG-Y，仪器前放增益用24dB，0~500Hz全频带接收，采样间隔1ms，记录长度1.5s。线束方向布置遵循一般垂直地层走向和主要构造走向的原则，呈东西向布置。

2地震数据处理

本次地震资料处理是在我公司SunBIade2000工作站上进行，使用了多套大型处理软件，本着“高信噪比、高分辨率、高保真度”的原则，加强对处理模块和参数的反复试验，确定了本次地震资料处理流程。

2.1初至折射静校正

由于地表高程及地表低（降）速带厚度、速度的横向变化使得地震波旅行时差会对信号的叠加效果产生一定的不利影响，致使反射波同相轴信噪比下降、频率降低。应用合适的静校正模块和参数，可以消除这种时差，确保叠加剖面的质量。测区属典型复杂山区，最低高程481.99m，最高高程919.61m，最大高差为437.62m，地表低、降速带地层的横向变化较大。经过对比试验，初至折射校正法可以较好解决该地区的静校正问题。

2.2反褶积技术

为了消除大地的滤波作用，拓宽频带，压缩地震子波，提高地震资料的纵向分辨率，经大量的测试对比后，选择了地表一致性预测反褶积。该方法是基于地震子波可以被分解为共炮点、共接收点、共偏移距、共反射点等多种成份的思想，它不仅能压缩地震子波，而且能进一步消除地表条件的变化对地震波的振幅特性和相位特性的影响，同时对多次波也有压制作用。由于反褶积在提高分辨率的同时将会降低资料信噪比，所以处理时在保证资料信噪比的情况下再提高分辨率。经对比分析，最终选定的处理参数为：预测步长8ms，算子长度150ms。

2.3速度分析和剩余静校正

速度分析的精度和剩余静校正的准确程度是相互影响的，为求取更准确的叠加速度场，必须求取准确的剩余静校正量。剩余静校正的求取是在给定的条带、时窗、倾角范围内，在纵横两个方向进行倾角调查形成模型道，在共炮点道集和共检波点道集做互相关，求取每个炮点、检波点的剩余静校正量，因此在选择时窗时要对层位的拾取进行试验。剩余静校正和速度分析是一个反复迭代的过程，迭代的次数在一定程度上影响着处理的精度。在本次资料处理过程中，进行了3次迭代分析，剩余静校正后剖面目的层同相轴连续性明显提高，剖面质量得到了明显改善。

2.4叠后去噪

叠后去噪是利用相邻地震道之间的反射波具有相似性而干扰波不具有相似性的这种特点使相干波得以加强的处理方法。利用相关函数测定2个以上记录道时间序列的相似程度，按相似程度的大小对记录进行加权，使得相似性好的波得以加强。本次资料处理选择多项式拟合方法来提高叠后信噪比。提高信噪比的同时影响剖面的分辨率，在叠后采用谱白化方法进一步提高频率，得到高分辨高信噪比资料。

2.5三维偏移

三维偏移的主要目的是消除地下倾斜界面对反射波的影响，使之成像归位到真实的反射界面位置上去，从而正确地反映地下形态和构造变化情况。有限差分偏移的效果主要决定于偏移速度，我们选用叠加速度经过转换建立偏移速度模型，并进行了反复测试和调整。用人工剔除奇异值，采用机器自身平滑，对速度的百分比进行偏移试验。经对比，95%的偏移效果较好。

三采空区的时间剖面特征和切片属性特征分析

1采空区的时间剖面特征

三维地震勘探成果包含了丰富的地质信息，其中包含采空区的地质信息，用来解释煤层的采空区范围。在时间剖面上采空区表现有3种现象，其一是煤层反射波变弱，在采空区边界处反射波同相轴频率和产状发生突变，在采空区内部反射波同相轴不连续且杂乱无章；其二是煤层反射波同相轴的变弱，频率和产状变化较大，与周围非采空区煤层反射波存在明显差别，而且煤层反射波之下的层位反射波同相轴增强，频率和产状突变，从整张剖面来看无采空区的地段煤层反射波较强，其下伏层位反射波很弱，而有采空区的地段煤层反射波同相轴表现微弱，其下伏层位反射波同相轴较强，形成明显反差；其三在地震时间剖面上表现为煤层反射波缺失。以韩城桑树坪井田为例，根据前人资料和矿方提供的采掘工程资料得知，井田内部，3号煤层正常采掘范围比较集中，2号煤层和3号煤层的老窑采空区比较分散。本次三维地震勘探圈定了采空区主要位于测区的东部以及中南部。根据采空区在三维地震勘探资料上的显示特征，解释了采空区的范围，面积约0.98km2。

2采空区的切片属性特征

地震属性技术的关键在于属性提取，提取方式包括同相轴属性提取和数据体属性提取。同相轴属性是与某个界面有关的地震属性，具体提取方法包括瞬时提取法、单道分时窗提取法和多道分时窗提取法。利用地震属性分析，有利于采空区边界的界定，利用Geoframe解释软件，提取了26种地震属性，选择了对采空区反映较好的弧长、均方根振幅、最大振幅和能量和四种地震属性参与了解释，对采空区范围的判别和划分起到了很好的辅助解释作用。

四结论

三维地震勘探技术应用于韩城桑树坪井田采空区探测获得了很好的效果，诸多成果得到了钻孔和实际采掘工程的验证，进而在三维地震勘探技术应用于复杂地形采空区勘探领域取得了重要的技术和生产经验。尤其值得注意的是，在数据解释方面，时间剖面解释和切片属性分析在采空区判断和范围划分具有同等重要的地位，二者相辅相成，兼顾彼此是提高采空区解释精度的关键。

作者：刘芳晓赵禄顺李元杰王晓亮单位：中国煤炭科工集团西安研究院有限公司