[發(fā)明專利]采集腳印噪聲去除方法、電子設(shè)備及計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)

在審

• 申請(qǐng)?zhí)枺?/td>: 202110913600.7
• 申請(qǐng)日: 2021-08-10
• 公開（公告）號(hào): CN113640882A
• 公開（公告）日: 2021-11-12
• 發(fā)明（設(shè)計(jì)）人: 李緒深;霍守東;張海榮;黃亮;石太昆;周旭暉;鄒佳儒;穆盛強(qiáng);舒國旭
• 申請(qǐng)（專利權(quán)）人: 南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室（湛江）;中國科學(xué)院地質(zhì)與地球物理研究所
• 主分類號(hào): G01V1/36
• 分類號(hào): G01V1/36
• 代理公司: 北京眾達(dá)德權(quán)知識(shí)產(chǎn)權(quán)代理有限公司 11570
• 代理人: 梁凱

本發(fā)明公開了一種采集腳印噪聲去除方法、電子設(shè)備及計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，其中所述方法包括：采集多個(gè)地震道對(duì)應(yīng)的初始地震數(shù)據(jù)；對(duì)所述初始地震數(shù)據(jù)進(jìn)行頻率波數(shù)徑向變換，得到所述初始地震數(shù)據(jù)的主能量?jī)A角；根據(jù)所述初始地震數(shù)據(jù)的主能量?jī)A角，采用多方向矢量中值濾波對(duì)所述初始地震數(shù)據(jù)進(jìn)行去噪處理，從而獲得去除所述采集腳印噪聲后的目標(biāo)地震數(shù)據(jù)。采用本申請(qǐng)，能精準(zhǔn)、有效地去除地震數(shù)據(jù)中的采集腳印噪聲，進(jìn)而提升地震數(shù)據(jù)質(zhì)量。

技術(shù)領(lǐng)域

本發(fā)明涉及地球物理勘探技術(shù)領(lǐng)域，尤其涉及一種采集腳印噪聲去除方法、電子設(shè)備及計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)。

背景技術(shù)

隨著巖性油氣藏和隱蔽性油氣藏勘探的不斷深入，對(duì)地震成像分辨率和地震資料的細(xì)節(jié)刻畫提出了更高的要求。采集腳印是在地震數(shù)據(jù)采集過程中，由于采集或檢波點(diǎn)空間采樣不規(guī)則方式，導(dǎo)致地下面元之間地震屬性分布不均勻而產(chǎn)生的一種噪聲。

其中，陸地采集到的地震資料(地震數(shù)據(jù))，因?yàn)槊嬖g的方位角和炮間距在空間上的分布特征，使得從地下發(fā)射回來的噪聲泄露到疊加數(shù)據(jù)中產(chǎn)生采集腳印噪聲。而海洋采集到的地震數(shù)據(jù)，因?yàn)橛馉铍娎|隨海水搖擺漂移、洋流潮汐等采集因素影響，造成檢波點(diǎn)空間采樣不規(guī)則，導(dǎo)致產(chǎn)生采集腳印噪聲。由上可知：由于地震數(shù)據(jù)采取的處理手段導(dǎo)致剩余時(shí)差、偏移距或振幅的變化，往往導(dǎo)致采集腳印噪聲進(jìn)一步加強(qiáng)，影響后續(xù)的地震數(shù)據(jù)處理與解釋工作。因此，亟需提出一種地震數(shù)據(jù)除噪方案。

發(fā)明內(nèi)容

本申請(qǐng)實(shí)施例通過提供一種采集腳印噪聲去除方法，能精準(zhǔn)、有效地去除地震數(shù)據(jù)中的采集腳印噪聲，進(jìn)而提升地震數(shù)據(jù)質(zhì)量。

一方面，本申請(qǐng)通過本申請(qǐng)的一實(shí)施例提供一種采集腳印噪聲去除方法，所述方法包括：

采集多個(gè)地震道對(duì)應(yīng)的初始地震數(shù)據(jù)，所述初始地震數(shù)據(jù)為所述地震道在不同采樣點(diǎn)處采集的時(shí)間空間域的地震數(shù)據(jù)，且包括采集腳印噪聲；

對(duì)所述初始地震數(shù)據(jù)進(jìn)行頻率波數(shù)徑向變換，得到所述初始地震數(shù)據(jù)的主能量?jī)A角；

根據(jù)所述初始地震數(shù)據(jù)的主能量?jī)A角，采用多方向矢量中值濾波對(duì)所述初始地震數(shù)據(jù)進(jìn)行去噪處理，從而獲得去除所述采集腳印噪聲后的目標(biāo)地震數(shù)據(jù)。

可選地，所述對(duì)所述初始地震數(shù)據(jù)進(jìn)行頻率波數(shù)徑向FKR變換，得到所述初始地震數(shù)據(jù)的主能量?jī)A角包括：

對(duì)所述初始地震數(shù)據(jù)進(jìn)行頻率波數(shù)變換，以將所述時(shí)間空間域的初始地震數(shù)據(jù)變換為頻率波數(shù)域的中間地震數(shù)據(jù)；

在預(yù)設(shè)傾角范圍內(nèi)對(duì)所述中間地震數(shù)據(jù)進(jìn)行角度掃描，以在徑向方向上識(shí)別得到所述初始地震數(shù)據(jù)對(duì)應(yīng)的主能量?jī)A角。

可選地，所述對(duì)所述初始地震數(shù)據(jù)進(jìn)行頻率波數(shù)變換，以將時(shí)間空間域的所述初始地震數(shù)據(jù)變換為頻率波數(shù)域的中間地震數(shù)據(jù)包括：

對(duì)時(shí)間空間域的所述初始地震數(shù)據(jù)進(jìn)行非均勻傅里葉變換，得到頻率空間域的變換地震數(shù)據(jù)；

對(duì)所述頻率空間域的變換地震數(shù)據(jù)進(jìn)行非均勻傅里葉反變換，得到頻率波數(shù)域的中間地震數(shù)據(jù)。

可選地，所述對(duì)所述頻率空間域的變換地震數(shù)據(jù)進(jìn)行非均勻傅里葉反變換，得到頻率波數(shù)域的中間地震數(shù)據(jù)包括：

計(jì)算每個(gè)地震道對(duì)應(yīng)的所述初始地震數(shù)據(jù)的傅里葉系數(shù)；

從多個(gè)所述初始地震數(shù)據(jù)的傅里葉系數(shù)中，選擇最大傅里葉系數(shù)；

利用所述最大傅里葉系數(shù)，對(duì)所述頻率空間域的變換地震數(shù)據(jù)進(jìn)行傅里葉反變換重計(jì)算，以獲得所述頻率波數(shù)域的中間地震數(shù)據(jù)。

可選地，所述得到所述初始地震數(shù)據(jù)的主能量?jī)A角之后，所述方法還包括：

對(duì)多個(gè)所述初始地震數(shù)據(jù)的主能量?jī)A角進(jìn)行排序，以找出傾角大于預(yù)設(shè)閾值的所述初始地震數(shù)據(jù)的主能量?jī)A角；