[發明專利]非穩態地震規則噪聲壓制方法和系統

在審

• 申請號: 202310247121.5
• 申請日: 2023-03-14
• 公開（公告）號: CN116299712A
• 公開（公告）日: 2023-06-23
• 發明（設計）人: 李皓;邢孟江;楊圓圓;谷巍巍;趙東鳳;王文慶;邢孟道;錢江;呂海濤;王峣鈞;王磊;高鵬
• 申請（專利權）人: 電子科技大學長三角研究院（湖州）
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 成都虹盛匯泉專利代理有限公司 51268
• 代理人: 陳婷

1.一種非穩態地震規則噪聲壓制方法，其特征在于，包括：

S1、求取原始地震數據傾角時差；

S2、利用面波傾角時差對局部地震數據進行拉平處理，再進行數據能量與道間時差調整；

S3、求取f-k濾波器時間域算子，并對經步驟S2處理后的地震數據進行時間域f-k濾波；

S4、對經步驟S3得到的地震數據進行反道間時差調整、反能量調整以及反拉平處理，得到純面波的地震數據；

S5、將步驟S4得到的純面波的地震數據與原始地震數據相減，得到有效信號。

2.根據權利要求1所述的一種非穩態地震規則噪聲壓制方法，其特征在于，步驟S1具體包括以下分步驟：

S11、記炮點和檢波點在同一個排列的地震數據為d⁰(t,x)，從d⁰(t,x)上拾取規則干擾的速度v⁰，并計算炮點和檢波點不在同一個排列上的視速度v^*：

其中，t表示時間點，x表示地震道，θ為方位角；

S12、計算地震數據面波傾角時差：

其中，dt、dx分別為時間采樣間隔和空間采樣間隔。

3.根據權利要求2所述的一種非穩態地震規則噪聲壓制方法，其特征在于，步驟S2具體包括以下分步驟：

S21、記窗口內拉平的地震數據為并選取中間道作為參考道x_r；

S22、對窗口內地震道進行均方根能量調整，使之與參考道能量一致：

其中，q為時間采樣點數；

S23、計算窗口內每道與參考道的相關，取地震道與參考道最大相關處，作為道間靜態時差進行調整，使橫向一致性更好；窗口內每道與參考道的相關的計算式為：

其中，x_r為參考道，l為地震道x與參考道x_r之間需要調整的靜態時差；

S24、求取相位濾波器h₁(t,x)，使窗口中各道向參考道進行動態時差調整即相位調整，來消弱面波固有的頻散特性。

其中，obj為目標函數，x_r為參考道，F()為傅立葉變換。

4.根據權利要求3所述的一種非穩態地震規則噪聲壓制方法，其特征在于，步驟S3具體包括以下分步驟：

S31、濾波器時間域表達式為：

其中，f_N為截止頻率，v為視速度；

S32、對地震數據進行濾波：

其中，表示二維褶積算子。

5.一種非穩態地震規則噪聲壓制系統，其特征在于，包括：地震數據傾角時差模塊、地震數據正向調整模塊、濾波模塊、地震數據反向調整模塊以及減法模塊；地震數據正向調整模塊包括依次串聯的：拉平處理單元、數據能量調整單元、道間時差調整單元；地震數據反向調整模塊包括依次串聯的：反道間時差調整單元、反數據能量調整單元、反拉平處理單元；

地震數據傾角時差模塊用于計算原始地震數據面波傾角時差；

地震數據正向調整模塊基于傾角時差對地震數據進行拉平處理，然后進行數據能量調整、道間時差調整；

經地震數據正向調整模塊處理后的地震數據作為濾波模塊的輸入，濾波模塊的輸出作為地震數據反向調整模塊的輸入；

地震數據反向調整模塊的輸出為純面波的地震數據；

原始地震數據與純面波的地震數據作為減法模塊的輸入，減法模塊的輸出為有效信號。