[發明專利]子波處理二步法反褶積方法

無效

• 申請號: 201010139220.4
• 申請日: 2010-04-02
• 公開（公告）號: CN101813786A
• 公開（公告）日: 2010-08-25
• 發明（設計）人: 王西文;王宇超;胡自多;王述江;田顏燦;肖明圖
• 申請（專利權）人: 中國石油集團西北地質研究所
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 甘肅省知識產權事務中心 62100
• 代理人: 劉繼春

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)與需要調整相位區塊的疊前地震道進行褶積， 即可完成相位校正子波處理。