[发明专利]基于均炮源波场照明的逆时偏移空间振幅补偿方法

摘要

本发明是保持地震波动力学特征的逆时偏移空间振幅补偿方法。采用深度域层速度模型和正演炮集数据，选取震源子波在深度域进行单炮波动方程正演，存储所有时刻单炮源波场并计算单炮源波场照明，用单炮记录作为初始边界沿时间轴逆向正演得到接收波场，计算未补偿的单炮逆时偏移成像和均炮源波场照明，得到均炮源波场照明补偿后的单炮逆时偏移成像，将所有单炮逆时偏移成像进行叠加得到均炮源波场照明补偿后的叠加逆时偏移成像。本发明从构造成像同相轴相对保持、背景噪音压制、波形和沿层地震属性空变相对保持方面好于单炮源波场照明的振幅补偿方法。

主权项

基于均炮源波场照明的逆时偏移空间振幅补偿方法，特点是通过如下步骤实现：1）采集地震勘探数据，处理得到深度域层速度模型和正演炮集数据；2）利用下式在深度域进行单炮波动方程正演： 1 v 2 ∂ 2 P s ( X , t ) ∂ t 2 - Δ P s ( X , t ) = δ ( X - X s ) S ( t ) , 并存储所有时刻的单炮源波场；同时利用下式计算单炮源波场照明：ISource(X)＝∫Ps(X,t)Ps(X,t)dt；其中：v为地震波速度，Ps(X,t)为单炮源波场，Δ为Laplace算子，δ(X‑Xs)为单位脉冲函数，S(t)为震源子波，X=(x,z)为波场所在二维空间位置坐标，Xs=(xs,zs)为震源所在二维空间位置坐标，x和xs分别为波场和震源横坐标，z和zs分别为波场和震源纵坐标，t为时间；3）利用单炮地震记录作为初始边界条件，用下式从最大记录时刻开始，沿时间轴逆向进行波动方程正演，得到接收波场： 1 v 2 ∂ 2 P g ( X , t ) ∂ t 2 - Δ P g ( X , t ) = G ( X g , t ) ; 其中：Pg(X,t)为接收波场，G(Xg,t)为单炮地震记录，Xg=(xg,zg)为检波器所在二维空间位置坐标，xg为检波器横坐标，zg为检波器纵坐标；4）在接收波场逆推时，读取相应时刻单炮源波场，利用下式计算未补偿的单炮逆时偏移成像：I(X)=∫Ps(X,t)Pg(X,t)dt；5）利用下式计算均炮源波场照明： I ‾ Source ( X ) = ( Σ Shot I Source ( X ) ) / N Shot ; 其中：NShot为单炮成像孔径范围内的总炮数；6）利用下式计算均炮源波场照明补偿后的单炮逆时偏移成像： I SC ( X ) = S { I ( X ) } I ‾ Source ( X ) ; 其中：S为满足二阶导数连续的平滑算子以保证成像波形的连续性（所述的平滑算子采用三次样条函数插值平滑算子）；7）将所有单炮逆时偏移成像进行叠加，得到均炮源波场照明补偿后的叠加逆时偏移成像。