表层Q值确定性求取与空变补偿方法

2019-09-28曾华会田彦灿徐兴荣肖明图

石油地球物理勘探 2019年5期

苏勤曾华会田彦灿徐兴荣肖明图吴杰

(中国石油勘探开发研究院西北分院，甘肃兰州 730020)

0 引言

陆上地震勘探中，疏松表层介质对地震波高频成分有强烈的吸收衰减作用，导致地震资料垂向分辨率降低；由于表层岩性、速度、厚度的横向变化，吸收衰减程度会不同，造成地震波道间能量、频率和相位的不一致，影响同相叠加结果的横向分辨率。地震波的吸收和频散特性与介质的成分、饱和度、孔隙度等密切相关，研究衰减补偿时一般用品质因子Q描述这种吸收与衰减的总体效应。

目前改善地表一致性、提高分辨率的常规方法有基于地表一致性(或单道)的反褶积(含振幅补偿)和基于深层地层的Q补偿方法。前者基于统计方法，通常难以满足假设条件，因此对消除表层影响能力有限。而且反褶积存在不同程度的保真性问题，甚至产生副作用[1-5]。后者忽略了表层严重的吸收衰减作用，因此难以获得理想的处理成果。

近年来，对地震信号的吸收衰减研究日趋重视[6-7]，但大多数的地层Q值求取及吸收补偿方法针对的是深层成岩介质[8-13]，Q模型的估算难度大，存在不确定性；反Q滤波振幅补偿的稳定性、补偿结果的评价等问题都亟待解决。针对反Q滤波振幅补偿的稳定性问题，Wang[14]基于Kolsky-Futterman模型，从一维波动方程出发，考虑能量的吸收和频散，导出能量补偿和相位校正公式，实现了稳健有效的频率域反Q滤波方法。

随着油气勘探领域对地震数据精度要求的提高，学者们开始重视表层吸收补偿问题[15-20]。尹喜玲等[18]就沙漠地区的近地表吸收规律进行了初步探讨，认为多次反射和透射是造成近地表衰减的主要原因。于承业等[19]提出利用双井微测井资料，通过频移法估算近地表Q值，强调了近地表衰减补偿的重要性，指出现有处理系统的地表一致性振幅补偿和地表一致性反褶积是定性的，不能在时间、频率、空间三个域有效消除近地表影响。

为了有效解决表层吸收补偿问题，本文提出了一种确定性表层Q求取和空变吸收补偿方法。该方法基于地表一致性原理，实现了表层相对衰减系数统计估算、表层Q模型建立、表层Q标定以及稳定的空变Q补偿，有效解决表层横向多变的吸收问题。经模型计算验证了方法的正确性和保真性，目前已应用于新疆、青海、大庆等多个地区实际生产。生产实践表明：该方法能合理拓宽有效频带并大幅提高地震数据分辨率，且不降低数据信噪比。经表层吸收补偿后，数据一致性更高，波组关系更加明确，地质现象和构造特征更加清楚。井标定和属性分析等均证明了补偿结果的可靠性，相比常规处理方法效果明显。

1 方法原理

首先利用野外表层吸收调查数据的直达波计算调查点的表层实测Q值。在没有表层实测Q值的情况下，也可利用表层调查的速度模型计算近似Q值：基于共炮检域统计迭代求出炮检点相对衰减系数，并结合表层旅行时得到表层相对Q值；再利用实测Q值或近似Q值对相对Q值的大小和空间变化趋势进行标定，得到表层Q场；最后利用表层Q场和表层旅行时，通过有效稳定的Q补偿算法对地震数据进行表层Q空变补偿。表层吸收补偿即相当于将炮点和检波点从地表延拓至高速顶界的过程。该方法流程如图1示。

图1 表层Q空变补偿流程

1.1 求表层实测Q值

实测Q值的求取是表层吸收补偿方法研究的基础，其结果对计算整个工区的Q值起着约束标定作用。利用微测井数据井底道和地面道直达波的频率属性和振幅属性，通过改进的峰值频移法或谱比法计算实测Q值[13,19-20]。

改进峰值频率频移法：表层对地震波的吸收会导致峰值频率由高变低，通过分析井底道峰值频率fw到地面道峰值频率fg的移动可反演调查点的表层Q值，如图2所示。改进峰值频率频移法结合了质心频率频移法和峰值频率频移法的优点，目的是先求质心频率fwc，再通过式(1)将质心频率换算成峰值频率fw，提高算法的稳定性。

图2 峰值频率频移法示意图

(1)

峰值频率频移法的计算公式为

(2)

频谱比率法(改进的谱比法)：表层对地震波的吸收还会造成振幅减弱，谱比法通过分析井中道振幅A0与地面道振幅A的比值计算Q值

(3)

图3 谱比法示意图

实际应用中，应先评价吸收调查记录是否可用。由于表层的吸收作用，地面道信号相对井底道信号其振幅明显更弱，主频明显降低。有时由于采集中井中道的耦合不好，所求出的Q值比较大，这样的值不可用。野外表层Q值调查点的设计密度应以能控制区域的表层Q的长波长变化为宜。

对于表层较厚的沙漠区，也可以利用多项式拟合质心频移法估算近地表Q值[19]。

1.2 表层相对Q值

表层吸收会造成反射信号的振幅减弱。若表层速度低或旅行时长，则表现出较强的吸收特征，这种吸收衰减具有地表一致性特点，把这种比例关系定义为相对衰减系数，作为求取表层相对Q值的重要属性参数。具体方法为：先对原始数据去噪，选取给定时窗内信噪比较高的反射同相轴数据作为输入，计算每道的均方根振幅，然后统计求取共炮检点的平均振幅，迭代求出炮检点相对衰减系数R，并结合表层旅行时、根据谱比法得到表层相对Q值。

R相当于谱比法中衰减前、后的能量比值。检波点位置的R与表层Q值之间的关系式为

(4)

式中：t表示旅行时，根据表层模型求得；f表示地震波主频，可根据去噪后频谱分析计算得到；scale表示调节系数，该系数可改变Q的趋势大小，逼近于实测Q值。

实际工作中对R做异常值过滤处理，由式(4)可以推导出表层Q值为

(5)

1.3 表层Q场的建立

基于野外实测数据求取的实测Q值及通过大炮数据统计求取炮检点表层的相对Q值，采用图4给出的表层Q约束标定流程可实现相对Q值的标定约束，并建立表层空变Q场，标定结果不改变表层Q值的短波长，但校正了表层Q值的长波长问题。

除表层Q值外，表层旅行时也是表层Q补偿的重要参数，可通过静校正计算得到。具体过程如图4所示。

图4 表层Q约束标定流程图

求得了每个炮检点的Q和t，也就完整地建立了表层Q模型和补偿参数，为叠前数据表层Q空变补偿奠定了基础。

1.4 表层Q补偿

建立表层Q场后，用稳定的表层补偿算法就可以对叠前数据在频率域进行表层Q空变补偿

(6)

Wang[14]提出了稳定的Q补偿算法，对振幅项做了改进，提出稳定的振幅补偿量

上述Q补偿算法的稳定性体现在可以通过增益限制控制补偿的频带范围，防止高频噪声的过分补偿。由于补偿量随频率和时间的增大逐渐增大，通常的补偿算法都会对高频端过分补偿，造成噪声过量、信号失真，而Wang[14]采用对完全吸收的频率成分不再进行补偿的策略，避免了对高频噪声的过量补偿。该算法的另一优点是补偿振幅的同时也可调整相位，可以同时解决地层吸收造成的地震波能量损失和频散问题。

算法中的参数scale对补偿起着重要作用，在实际工作中需要根据具体情况给定合理的值。图5是参数scale中心频率对相位调整项的影响。由图可见，中心频率发生变化时，相位调整量也随之变化，仅仅在频率为ωh时相位不做调整，其他情况都需要调整。

图5 相位调整量随中心频率的变化曲线

图6是在固定Q、t的情况下，给定不同Glim时的振幅补偿曲线。由图可见，Glim控制了拐点频率(即截止频率fh)，在实际数据补偿中控制高频端的补偿量，因此合理设定Glim能降低补偿结果中的噪声水平。图7和图8分别为不同表层Q值、相同旅行时的振幅和相位补偿曲线，也进一步说明了表层Q补偿算法的自适应性。

图6 增益限制Glim对振幅补偿量影响曲线

图7 相同旅行时(t=60ms)、不同Q值时的振幅补偿曲线

图8 相同旅行时(t=60ms)、不同Q值时的相位补偿曲线

2 模型试验

2.1 模型1

为了测试方法的正确性及适应性，利用理论模型进行正、反演论证。参照一口实际井的数据建立黏弹模型。图9是表层速度模型和Q模型，表层Q是横向变化的。通过波动方程正演获得二维地震数据，利用本文方法对数据求取相对衰减系数、建立表层Q模型并验证表层Q补偿。图10为正演Q值和反演相对Q值曲线，可见两者吻合较好。图11是模型正演、反褶积、表层Q补偿道集及频谱对比图。需要说明的是试验采用的是地表一致性反褶积，参数是经扫描后优选得到的。模型1的计算结果证明了本文法求取Q值及表层Q补偿是准确且稳定的。

图9 表层速度理论模型(上)和Q值理论模型(下)

图10 理论模型正演Q值(蓝色)和反演相对Q值(红色)

图12为模型正演、反褶积、表层Q补偿叠加剖面及模型合成记录对比图，剖面中内嵌的蓝色信号为合成记录(模型)。可以看出表层Q补偿的数据与合成记录吻合性更好，分辨率更高。进一步验证了本文方法的保真性。

2.2 模型2

模型2来自于实际的表层吸收调查数据。图13为吸收调查的现场观测方式示意图，通过激发井和接收井构成双微测井调查，分别取井底道与地面道信号，用频移法或谱比法求取野外实测Q值。用得到的实测Q值补偿地面道，与井底道对比如图14所示。可以看出：地面道相对与井底道的衰减更显著、振幅更弱、主频更低，经反演Q值补偿的信号与井底道的信号很相似。模型2的计算结果进一步证实，表层吸收严重的信号经衰减处理后，通过表层补偿可以近似得到恢复。