[發(fā)明專利]一種基于初至子波提取的可控震源整形濾波方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201410316271.8
• 申請(qǐng)日: 2014-07-04
• 公開(kāi)（公告）號(hào): CN105319587A
• 公開(kāi)（公告）日: 2016-02-10
• 發(fā)明（設(shè)計(jì)）人: 關(guān)鍵;梁鴻賢;秦寧;石林光;李凌云;徐輝;傅金榮;高麗;孫朋朋;張傳強(qiáng);孫治國(guó);武天祥;王延光;王興謀;王靜軒;王敬閣;王磊;陳震林;孫淑琴;馬季;揭景榮
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號(hào): G01V1/30
• 分類號(hào): G01V1/30;G01V1/36;G01V1/28;G06F19/00
• 代理公司: 東營(yíng)雙橋?qū)＠碛邢挢?zé)任公司 37107
• 代理人: 侯華頌

1.一種基于初至子波提取的可控震源整形濾波方法，其特征是包括：

（1）選取多組野外同炮點(diǎn)位置激發(fā)，同排列接收得到的炸藥震源單炮和可控震源單炮；

（2）對(duì)步驟（1）選取的炸藥震源單炮和可控震源單炮分別進(jìn)行初至拾取；

（3）初至校平，根據(jù)拾取到的初至，將地震單炮沿初至校到同一時(shí)間；

（4）對(duì)校平后的炸藥震源及可控震源單炮直接進(jìn)行水平疊加，提取各自的初至子波；

（5）利用得到的多組子波統(tǒng)一求取可控震源整形濾波算子；

（6）將求得的可控震源整形濾波算子作用于可控震源的所有地震道，完成可控震源子波整形濾波的處理。

2.根據(jù)權(quán)利要求1所述的基于初至子波提取的可控震源整形濾波方法，其特征是：所述可控震源整形濾波算子的求取按下列方式進(jìn)行：

設(shè)可控震源與炸藥震源單炮的初至子波記錄分別為可控震源子波x_i(t)（i=1,2，…N）和炸藥震源子波y_i（t）（i=1,2，…N），其中i為炮號(hào)；N為參與運(yùn)算的總炮數(shù)；

設(shè)一個(gè)整形濾波算子m_i（t）作用于x_i(t)，使可控震源子波x_i(t)經(jīng)過(guò)褶積濾波后逼近炸藥震源子波y_i（t），設(shè)整形濾波后實(shí)際輸出x_i(t)*m_i（t）與期望輸出y_i（t）的誤差為e_i(t)，那么則有

（1）

誤差能量總和為

（2）

應(yīng)用最小二乘法原理，令總誤差能量E對(duì)m_i（t）的偏導(dǎo)數(shù)等于零，即

（3）

得到求解整形濾波算子的托布里茲矩陣方程

（4）

式中：R_xx表示輸入道x_i(t)的自相關(guān)函數(shù)矩陣；M為整形濾波算子向量；R_zr表示期望輸出道y_i(t)與輸入道x_i（t）的互相關(guān)函數(shù)向量，托布里茲矩陣方程的矩陣形式為

（5）

求解公式（5）的托布里茲矩陣方程，得到可控震源第i道的整形濾波算子m_i（t），對(duì)每組炸藥震源和可控震源初至子波得到的N個(gè)整形濾波算子進(jìn)行平均，就得到可控震源整形濾波算子

（6）

再將求得的可控震源整形濾波算子作用于可控震源的所有地震道，完成可控震源子波整形濾波的處理，

（7）。