[發(fā)明專利]一種地震數(shù)據(jù)處理方法及裝置

有效

• 申請(qǐng)?zhí)枺?/td>: 201810460300.6
• 申請(qǐng)日: 2018-05-15
• 公開（公告）號(hào): CN108897036B
• 公開（公告）日: 2020-04-10
• 發(fā)明（設(shè)計(jì)）人: 王偉;趙玉敏;李皓;周振曉;唐博文;張文波
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油天然氣集團(tuán)有限公司;中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司;中國(guó)石油大學(xué)(北京)
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28;G01V1/30
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司 11127
• 代理人: 李輝

1.一種地震數(shù)據(jù)處理方法，其特征在于，包括：

根據(jù)采集到的地震數(shù)據(jù)，獲得反映地層結(jié)構(gòu)特征的投影誤差算子；

根據(jù)所述投影誤差算子，構(gòu)建地震數(shù)據(jù)反演的目標(biāo)函數(shù)，所述目標(biāo)函數(shù)中包括：所述地震數(shù)據(jù)、去除噪聲的地震數(shù)據(jù)、所述投影誤差算子之間的函數(shù)關(guān)系；

對(duì)所述目標(biāo)函數(shù)中進(jìn)行迭代反演，獲得去除噪聲的目標(biāo)地震數(shù)據(jù)。

2.如權(quán)利要求1所述的一種地震數(shù)據(jù)處理方法，其特征在于，所述投影誤差算子包括：

上式中，

h(x,t)表示所述投影誤差算子，x表示空間方向上的參數(shù)，x＝-px,-(px-1),…,(px-1),px，t表示時(shí)間方向上的參數(shù)，t＝-pt,-(pt-1),…,(pt-1),pt，px是空間方向的算子長(zhǎng)度，pt是時(shí)間方向的算子長(zhǎng)度，D表示地震數(shù)據(jù)組成的列向量，D＝d(nx,nt)，nx＝1,2,…N,nt＝1,2,…T，N表示所述地震數(shù)據(jù)的總道數(shù)，T表示每道采樣點(diǎn)數(shù)，I表示單位矩陣，λ表示權(quán)重系數(shù)，A表示所述地震數(shù)據(jù)重新排布后的矩陣，A的第i行的元素為d(nx,nt),m≤nx≤m+2px且nx≠m+px，n≤nt≤n+2pt，m＝(i-1)|(N-2px)+1，n＝(i-1)\(N-2px)+1。

3.如權(quán)利要求1所述的一種地震數(shù)據(jù)處理方法，其特征在于，所述目標(biāo)函數(shù)包括：

J＝||d-s||²+μF(hs)

上式中，J表示所述目標(biāo)函數(shù)，s表示所述去除噪聲的地震數(shù)據(jù)，d表示所述地震數(shù)據(jù)，μ表示反演權(quán)重系數(shù)，h表示所述投影誤差算子的褶積矩陣，F(xiàn)表示稀疏約束函數(shù)。

4.如權(quán)利要求3所述的一種地震數(shù)據(jù)處理方法，其特征在于，所述稀疏約束函數(shù)包括：

上式中，s表示所述去除噪聲的地震數(shù)據(jù)，h表示所述投影誤差算子的褶積矩陣，L表示(hs)的個(gè)數(shù)，δ表示尺度因子。

5.如權(quán)利要求4所述的一種地震數(shù)據(jù)處理方法，其特征在于，所述對(duì)所述目標(biāo)函數(shù)中進(jìn)行迭代反演，獲得去除噪聲的目標(biāo)地震數(shù)據(jù)，包括：

設(shè)置所述去除噪聲的地震數(shù)據(jù)的初始值、所述反演權(quán)重系數(shù)的取值、所述尺度因子的取值；

利用如下公式循環(huán)迭代求解所述去除噪聲的地震數(shù)據(jù)，直至所述目標(biāo)函數(shù)的取值滿足迭代終止條件，將滿足迭代終止條件時(shí)對(duì)應(yīng)的去除噪聲的地震數(shù)據(jù)作為所述目標(biāo)地震數(shù)據(jù)；

上式中，s表示所述去除噪聲的地震數(shù)據(jù)，I表示單位矩陣，λ表示所述權(quán)重系數(shù)，Q^k表示反演對(duì)角距陣，表示所述反演對(duì)角距陣中的對(duì)角元素，k表示迭代次數(shù)，h表示所述投影誤差算子的褶積矩陣，d表示所述地震數(shù)據(jù)。

6.如權(quán)利要求5所述的一種地震數(shù)據(jù)處理方法，其特征在于，所述目標(biāo)函數(shù)的取值滿足迭代終止條件，包括：

當(dāng)所述目標(biāo)函數(shù)的取值滿足時(shí)，確定滿足所述迭代終止條件；

上式中，J^k表示第k次迭代時(shí)所述目標(biāo)函數(shù)的取值，Δ表示誤差標(biāo)準(zhǔn)。

7.一種地震數(shù)據(jù)處理裝置，其特征在于，包括：

誤差算子獲取模塊，用于根據(jù)采集到的地震數(shù)據(jù)，獲得反映地層結(jié)構(gòu)特征的投影誤差算子；

反演目標(biāo)函數(shù)構(gòu)建模塊，用于根據(jù)所述投影誤差算子，構(gòu)建地震數(shù)據(jù)反演的目標(biāo)函數(shù)，所述目標(biāo)函數(shù)中包括：所述地震數(shù)據(jù)、去除噪聲的地震數(shù)據(jù)、所述投影誤差算子之間的函數(shù)關(guān)系；

迭代反演模塊，用于對(duì)所述目標(biāo)函數(shù)中進(jìn)行迭代反演，獲得去除噪聲的目標(biāo)地震數(shù)據(jù)。