[發明專利]計算復雜地表區的浮動基準面靜校正量的方法

在審

• 申請號: 201810930300.8
• 申請日: 2018-08-15
• 公開（公告）號: CN110837120A
• 公開（公告）日: 2020-02-25
• 發明（設計）人: 余青露;居興國;李進;鄒少峰;肖盈;高艷霞;祝媛媛;劉思思
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;張杰

本發明涉及用于計算復雜地表區的浮動基準面靜校正量的方法，所述浮動基準面靜校正量為三維地震數據中每個CMP點相對固定基準面的靜校正量，所述方法包括：通過平均靜校正量法計算出所述每個CMP點的浮動基準面靜校正量；利用非局部均值濾波法在空間域平滑所述每個CMP點的浮動基準面靜校正量。

技術領域

本發明總體涉及地震信號處理技術領域，尤其涉及一種基于非局部均值濾波的用于計算復雜地表區的浮動基準面靜校正量的方法。

背景技術

在復雜地表區以及高精度三維地震資料的地震數據處理時，會因近地表變化引起靜校正量劇烈變化，而造成反射波雙曲線型時距曲線的畸變。通常采用浮動基準面校正法來消除這種影響。目前常用的浮動基準面靜校正量計算方法有平均靜校正量法和平滑地表法，不同的計算方法得到的浮動基準面靜校正量有所差別。對于不同的地區，由于近地表起伏不同，采用的浮動基準面靜校正量計算方法也不一樣。目前，比較常用的浮動基準面計算方法是平均靜校正量法。在我國南方和中西部地區，由于地表起伏劇烈，計算出的平均靜校正量會產生突變。

在復雜地表區，采用平均靜校正量法計算浮動基準面來進行地震資料疊前時間偏移處理時，剖面上局部會存在鋸齒現象，這是由于靜校正量的突變引起的。因此，需要采用更加平滑的浮動基準面來進行計算。也就是說，采用平均靜校正量法計算時，需要在空間域平滑各CMP(Common Middle Point，共中心點)點的靜校正量。

發明內容

為了解決上述技術問題，本發明提供了一種基于非局部均值濾波的浮動基準面靜校正量的計算方法，該方法可以消除靜校正量的異常突變點，更好的保護信號的結構信息，并且算法效果較好，對于變化較大的離散點不產生影響。

根據本發明的一個方面，提供了一種用于計算復雜地表區的浮動基準面靜校正量的方法，所述浮動基準面靜校正量為三維地震數據中每個CMP點相對固定基準面的靜校正量，所述方法包括：

通過平均靜校正量法計算出所述每個CMP點的浮動基準面靜校正量；

利用非局部均值濾波法在空間域平滑所述每個CMP點的浮動基準面靜校正量。

在一實施例中，通過平均靜校正量法計算出所述每個CMP點的浮動基準面靜校正量包括：

對每個CMP道集的各道的炮點和檢波點相對固定基準面的靜校正量進行求和平均，得到每個CMP點的浮動基準面靜校正量。

在一實施例中，所述炮點和檢波點相對固定基準面的靜校正量T_m的計算公式為：

其中，E_g(m)是第m個炮點或檢波點的高程；h_nm為第m個炮點或者檢波點對應的第n層低速帶厚度，V_n是低速帶各層速度；N為低速帶層數；v_c是靜校正替換速度；E_d是固定基準面的高程。

在一實施例中，在計算所述每個CMP點的浮動基準面靜校正量時加入所述每個CMP點的權重值。

在一實施例中，利用非局部均值濾波法在空間域平滑所述每個CMP點的浮動基準面靜校正量包括：

選取控制濾波強度的參數，使得經平滑的所述每個CMP點的浮動基準面靜校正量與通過平均靜校正量法得到的所述每個CMP點的浮動基準面靜校正量相比，二者形態相似。

在一實施例中，利用非局部均值濾波法在空間域平滑所述每個CMP點的浮動基準面靜校正量包括：

假設p點是濾波器要處理的點，當濾波器遍歷到p點時，需要對p點鄰域內的像素點進行遍歷，假設q點是p點鄰域內的一個像素，算法遍歷到q點，則需要計算p點和q點之間的權重值，其包括：