[发明专利]子波处理二步法反褶积方法

权利要求书

1.一种子波处理二步法反褶积方法，其特征在于按以下步骤进行：

第一步、将采集于不同施工年代、不同施工参数的三维地震资料按常规方法进行分块预 处理，包括静校正处理、振幅衰减补偿处理、噪音压制处理；

第二步、地震子波频宽处理

假设同一位置的不同时间采集的不同信号，包含震源不同信号，是反射系数与地震子波 的褶积，

第一块三维采集的地震信号为：s₁(t)＝w₁(t)*r(t)

第二块三维采集的地震信号为：s₂(t)＝w₂(t)*r(t)

……

第N块三维采集的地震信号为：s_n(t)＝w_n(t)*r(t)       (1)

其中s(t)是用检波器接收到的地震子波，w(t)为激发的震源子波，r(t)是地层反射系 数，假定同一位置的地震反射系数r(t)是不变的，只是子波受到不同采集参数的影响；首先， 不考虑子波相位的差异，认为子波只是频宽不同；

根据子波的特征，采用反褶积算子a(t)，对地震记录s₁(t)和s₂(t)进行处理，相当于对地 震子波w₁(t)、w₂(t)……w_n(t)进行处理，得到处理后的w_a1(t)、w_a2(t)……w_an(t)：

s_a1(t)＝a(t)*s₁(t)＝a(t)*(r(t)*w₁(t))

＝r(t)*(a(t)*w₁(t))

＝r(t)*w_a1(t)

s_a2(t)＝a(t)*s₂(t)＝a(t)*(r(t)*w₂(t))

＝r(t)*(a(t)*w₂(t))

＝r(t)*w_a2(t)

……

s_an(t)＝a(t)*s_n(t)＝a(t)*(r(t)*w_n(t))

＝r(t)*(a(t)*w_n(t))              (2)

＝r(t)*w_an(t)

将反褶积算子a(t)设计成零相位反褶积算子，进行地震子波频宽处理；

第三步、将地震子波频宽处理后的资料进行常规的速度分析与剩余静校正、叠加处理；

第四步、求取相位校正反褶积算子

设s_a1(t)，s_a2(t)为经过第二步地震子波频宽处理后的任意地理位置上相邻的两个区块的 叠前或叠后地震记录信号，

s_a1(t)＝m(t)*s_a2(t)          (3)

式中，m(t)为相位校正反褶积算子；再利用式(4)求相位校正反褶积算子m(t)：

Q=Σt=0K+M[sa1(t)-Στ=0mm(t)sa2(t-τ)]2---(4)]]>

Q为趋于零的极小值，根据函数求极值的原理，

令s＝0，1，2…，M    (5)可列出M+1个方程组，

组成托布利兹矩阵方程组：

Στ=0mraa(τ-s)sa1(τ)=rab(s)---(6)]]>

式中：r_aa(τ-s)是时间延迟为τ-s的地震信号s_a1(t)自相关数列；r_ab(s)是时间延迟 为s的地震信号s_a1(t)与s_a2(t)的互相关数列；用递推法解托布利兹矩阵方程(6)，得到相 位校正反褶积算子m(t)；

所述s_a1(t)、s_a2(t)采用叠加后地震道记录；叠加后的地震道记录和叠加前的地震道记录 与相位校正反褶积算子m(t)的关系如下式：

sa2(t)*m(t)={sa21(t)+sa22(t)+...+sa2n(t)}*m(t)]]>

=sa21(t)*m(t)+sa22(t)*m(t)+...sa2n(t)*m(t)]]>

式中s_a2(t)为叠加后的地震道记录，为叠加前的地震道记录；

第五步、用求取的相位校正反褶积算子m(t)与需要调整相位区块的叠前地震道进行褶积， 即可完成相位校正子波处理。