[发明专利]基于立体成像条件的一次波和多次波联合成像方法及系统

权利要求书

1.一种基于立体成像条件的一次波和多次波联合成像方法，其特征在于，所述成像方法包括：

获取地表观测点的地震数据，所述地震数据包括一次波和多次波；

利用表层多次波去除方法从所述地震数据中提取震源子波；

根据所述地震数据和所述震源子波，采用单程波波动方程沿深度轴进行正向波场延拓，确定深度为z处的波场，记为延拓后震源波场；

根据所述地震数据，采用单程波波动方程沿深度轴进行反向波场延拓，确定深度为z处的波场，记为延拓后检波波场；

利用局部倾斜叠加方法分别对所述延拓后震源波场和所述延拓后检波波场进行分解，得到不同局部斜率的震源波场和不同局部斜率的检波波场；

根据设定的速度场，利用平面波解构滤波器确定震源波场和检波波场局部斜率的取值范围；

在所述局部斜率的取值范围内，将分解的震源波场和检波波场进行互相关确定成像结果。

2.根据权利要求1所述的基于立体成像条件的一次波和多次波联合成像方法，其特征在于，所述利用表层多次波去除方法从所述地震数据中提取震源子波，具体包括：

根据所述地震数据，基于反馈迭代的表层多次波去除方法确定地表加权算子；

根据所述地表加权算子，采用傅里叶反变换确定所述震源子波。

3.根据权利要求2所述的基于立体成像条件的一次波和多次波联合成像方法，其特征在于，所述根据所述地震数据，基于反馈迭代的表层多次波去除方法确定地表加权算子，具体公式为：

ΔP⁽ⁿ⁺¹⁾(z₀)＝P(z₀)-A⁽ⁿ⁺¹⁾(w)ΔP⁽ⁿ⁾(z₀)P(z₀)

其中，n表示迭代次数，n＝0,1,2...，z₀表示地表观测点的深度，ΔP⁽ⁿ⁺¹⁾(z₀)表示第n次迭代后一次波数据，ΔP⁽⁰⁾(z₀)＝P(z₀)，P(z₀)表示地震数据，A⁽ⁿ⁺¹⁾(w)表示第n次迭代后地表加权算子，w表示角频率。

4.根据权利要求1所述的基于立体成像条件的一次波和多次波联合成像方法，其特征在于，所述单程波波动方程具体为：

其中，u(x,z,t)表示波场，v(x,z,w)表示地震波速度，+符号表示上行波场，－符号表示下行波场，x表示水平方向坐标，z表示垂直方向坐标，t表示地震波传播时间，w表示角频率。

5.根据权利要求1所述的基于立体成像条件的一次波和多次波联合成像方法，其特征在于，所述在所述局部斜率的取值范围内，将分解的震源波场和检波波场进行互相关确定成像结果，具体成像公式为：

I(x,z)＝∫∫W_S(x,z,p,t)W_D(x,z,p,t)dpdt

其中，p表示空间和时间的局部斜率三维矢量，W_S(x,z,p,t)表示分解的震源波场，W_D(x,z,p,t)表示分解的检波波场，x表示水平方向坐标，z表示垂直方向坐标，t表示地震波传播时间。