山区煤层气三维地震成像关键技术研究

黄素荷 张 凯 袁兴赋

（安徽省勘查技术院，安徽 合肥 230031）

郑庄D区块位于沁水煤层气田西部的主体部位，地震地质条件十分复杂，地形起伏大，表层岩性横向变化，部分地段煤层埋深浅，致使有效覆盖次数低，小断层及陷落柱发育。种种因素造成全区资料存在严重静校正问题，局部地段信噪比低以及偏移成像难度大，加上煤层气开发对地震成像精度要求高，需要进行煤储层的裂缝预测和含气性检测等[1]，给该区三维地震资料处理带来严峻的挑战。因此，需要采用不同的地震数据处理措施，形成有效的针对性的配套技术，才能达到提高地震成像品质的目的。

1 高精度静校正联合

D区块地形高差大，以中低山为主，总体呈两山夹一河的地势，海拔高程在600～1150m之间，平均海拔标高在900m左右。表层岩性变化大，主要激发岩性有砂岩、泥岩、砾石、胶泥、黄土等五类。地形起伏及地表岩性变化等因素造成原始资料静校正问题非常严重，静校正问题是需要优先解决的问题。本次解决静校正问题的配套技术：初至波拾取修改→初至波折射静校正→地表一次性剩余静校正→非地表一次性静校正。首先应用折射静校正解决长波长及中波长静校正问题，其次应用地表一次性剩余静校正解决短波长静校正问题，最后应用非地表一次性静校正进一步解决局部残余短波长静校正问题。

1.1 初至波折射静校正

D区块大多地段为基岩出露，虽然地形起伏大及表层岩性变化，但基本能够拾取同一折射层的初至波。初至波拾取的质量直接影响静校正的计算结果，由于本区块覆盖次数不高，设计为30次，拾取初至波按照程序计算、手工逐道修改的办法，并约束初至波的拾取：

（1）在全区进行初至波极性调查，判断采用何种极性进行初至波拾取。

（2）在初至波清楚的单炮上拾取初至参考线，影射到其他单炮，指导初至拾取。

（3）调出地形线，判断拐曲位置的初至是否与地形起伏对应。

（4）加滤波、动平衡等处理使得初至波变得比较清楚，利于正确拾取。

（5）同时调出多个单炮，互相比较，保证所拾取的初至波是同一个相位。

（6）对实在不好辨认的初至舍弃，保证参与计算初至波时间的准确率。

（7）拾取多个初至文件进行静校正，比较静校正效果，并及时修改初至波。

利用拾取的初至波时间信息反演近地表模型，然后求解静校正量。静校正量需要分解成低频分量和高频分量，先进行高频分量的校正，叠前时间偏移后的CRP道集在叠加后再进行低频分量的校正。

D区块应用初至波折射静校正后，图1（a）显示单炮记录畸变的同相轴连续性得到改善，地形和低速层引起的同相轴错动现象基本消除。图1（ b）显示叠加剖面信噪比明显提高，煤层反射波组特征突出，影响构造的中长波长静校正问题得到了较好的解决。

图1 折射静校正前后单炮记录及叠加剖面对比

1.2 地表一致性剩余静校正

D区块采用地表一致性最大能量剩余静校正来消除残余短波长静态时差的影响。该方法的主要参数是相关时窗和最大校正量，相关时窗应选反射波能量强、波形稳定、连续性好的层段；最大校正量选取小于相关时窗内主周期的一半，剩余静校正后的叠加剖面不应有串相位现象。将剩余静校正与速度分析迭代进行，直到90%以上的校正量小于一个采样间隔。图2（a）显示通过剩余静校正后，速度分析的监控道集和叠加段有效波相位更加连续，速度谱的能量团更集中，速度分析的精度得到改善，从而拾取的速度更真实准确。图2（b）显示通过剩余静校正后，剖面同相轴的连续性明显变好。

图2 剩余静校正前后速度分析及叠加剖面对比

1.3 非地表一致性静校正

D区块由于基岩出露、地形起伏、低降速带速度及厚度变化等因素的影响，造成实际情况与地表一致性假设的条件存在偏差。经过折射静校正和多次地表一致性剩余静校正后，虽然静校正问题已基本解决，但局部资料还有残余短波长时差的影响，非地表一致性静校正能够较好地解决这类残余静态时差的影响。

在一些地表条件复杂地区，平滑静校正能够作为进一步解决静校正问题的补充。使用平滑静校正一定要对计算时窗、相关道数和校正时差进行充分的试验，一定要对叠加结果进行严格的质量控制，确保叠加剖面只是相位的连续性有所改善而构造形态没有明显的异常变化，只能在静校正问题基本解决的情况下使用[2]。

在D区块处理中，平滑静校正采用大时窗、小相关道数、小校正时差。时窗为500ms，相关道数为11道，校正时差很小为1ms，并对叠加剖面进行了严格的质量控制。图3显示叠加剖面除波组的连续性略有改善外，构造形态没有明显异常变化，说明参数选取合理。

图3 平滑静校正前（左）后（右）叠加剖面对比

2 叠前多域“递进”去噪

针对干扰波类型与特点，采用先强后弱、先规则后随机的去噪顺序在多域“递进”去噪。噪声压制过程为：强能量干扰波衰减→线性干扰波衰减→残余强能量干扰波衰减→随机干扰波衰减→叠后残余干扰波衰减。单炮上的强振幅干扰如声波干扰、猝发脉动、异常扰动、工业干扰、环境干扰及低速低频面波等强能量干扰采用地表一致性分频压噪技术在炮域、检波点域、CMP域进行压制。

针对具有线性特征的噪声，在炮域采用相干噪音衰减进行压制。基本思路是根据线性干扰波与有效波之间在视速度、频率、位置和能量上的差异，在T-X域采用线性预测方法确定线性干扰的视速度、分布范围及规律，将规则干扰自动识别出来，从原始数据中减去，实现线性干扰波的滤除，被滤波的部分主要集中在被干扰波覆盖的区域，其他部分不受影响。

随机干扰噪音在记录上普遍存在，根据随机性特点与有效波相干性特点的区别，选用随机噪音衰减技术有效地压制。反褶积之前在炮域采用“噪音自动识别与压制技术”进行一次随机噪音衰减，以保证反褶积的处理效果；反褶积后在炮域和共偏移距域再次进行随机噪音衰减，有效地提高叠前资料的信噪比，为后续的叠前时间偏移奠定良好的基础。

图4显示，单炮记录经过综合去噪后，各种干扰波被压制，反射信息突出。

图4 干扰波综合去除前（上）后（下）单炮记录对比

3 叠前时间偏移

叠前时间偏移输入为水平叠加之前的非零炮检距CMP道集数据，比叠后偏移更适用于复杂构造区，在地层速度横向变化不是很显著的构造复杂区能够取得很好的成像效果[3]。D区块两套含煤地层基本为层状介质，速度场横向变化不大，采用“渐进+迭代”方式速度建模，提高了小断层和陷落柱等地质构造成像的准确性。

3.1 “渐进+迭代”偏移速度建模

（1）在资料预处理的主要阶段都要进行速度分析，利用速度分析与剩余静校正迭代来逐步提高速度场质量。

（2）将用于叠前时间偏移的道集进行DMO速度分析，进一步改善速度场质量。

（3）DMO叠加速度作为初始速度场进行偏移处理，输出CRP道集。

（4）根据CRP道集校平与否来判断速度的合理性，采用目标线剩余速度分析与偏移的迭代求出较正确的偏移速度，直到获得满意的偏移效果，才得到最终的偏移速度场。

3.2 偏移参数的选择

GEODEPTH系统中克希霍夫积分叠前时间偏移的重要参数即为偏移孔径、拉伸因子及反假频因子。

（1）偏移孔径的选择。偏移孔径大小主要与地层倾角、目的层深度和速度有关。过小会出现假频，大倾角地层不能正确成像；过大深部会带来大量的偏移噪声，延长偏移时间。

（2）拉伸因子的选择。过大产生弧状噪声，浅部地层成像效果信噪比变差；过小产生假频，模糊构造，倾斜地层不能成像。

（3）反假频因子的选择。过大倾斜构造不能成像，产生假频；过小产生弧状噪声，成像效果变差，浅部地层反映特别明显。

4 处理效果分析

配套处理技术对关键环节进行处理，较好地解决了静校正问题，提高了剖面的信噪比，资料的成像品质获得了不错的效果。最终成果剖面信噪比高，保真度高，分辨率适中；反射波组层次清晰，波组特征明显，易于追踪对比；断裂刻画清晰，断点干脆，断裂的交互关系和结构清楚，如图5（a）；主要勘探目的层内幕反射清晰，剖面揭示的构造形态可靠，陷落柱特征明显，如图5（b）。

图5 叠前时间偏移效果

5 结论

（1）初至时间的拾取质量是影响初至折射波静校正效果的关键因素，多步骤多方法的有机联合能够较好地解决山区煤层气三维地震勘探中静校正的影响。

（2）小断层及陷落柱等构造成像需要叠前资料有较好的保幅性，叠前噪声压制要分类别、分阶段“递进”式进行。

（3）在叠前时间偏移中，采用“渐进+迭代”模式建立偏移速度模型，能够改善速度模型质量，提高成像精度，对地质构造有较好的成像效果。