[發明專利]一種射線參數域轉換波靜校正方法及系統

在審

• 申請號: 201811033015.2
• 申請日: 2018-09-05
• 公開（公告）號: CN110879413A
• 公開（公告）日: 2020-03-13
• 發明（設計）人: 朱童;雷朝陽;李京南
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種射線參數域轉換波靜校正方法，其特征在于，包括：

1)輸入共檢波點道集；

2)將步驟1)中輸入的共將波點道集進行τ-p變換，得到射線參數域道集，并分選為共射線參數道集；

3)對步驟2)中獲得的共射線參數道集進行沿構造趨勢平滑處理獲得平滑后的共射線參數道集；

4)將步驟2)中獲得的共射線參數道集和步驟3)中獲得的平滑后的共射線參數道集互相關，獲得靜校正量道集；

5)將步驟2)中的共射線參數道集與步驟4)中的靜校正量道集進行褶積，獲得靜校正后的共射線參數域道集，并分選為常規射線參數域道集；

6)將步驟5)中獲得的常規射線參數域道集進行反τ-p變換，得到靜校正后的共檢波點道集，并分選為共炮點道集。

2.根據權利要求1所述的射線參數域轉換波靜校正方法，其中，在步驟2)中，利用如下公式進行τ-p變換：

其中，p表示射線參數，τ表示截距時間，u(t,x)表示地震信號，表示τ-p正變換，表示希爾伯特變換的導數。

3.根據權利要求2所述的射線參數域轉換波靜校正方法，其中，所述截距時間τ通過如下公式表達：

其中，q_i為第i層的垂直慢度，q_i＝cos(θ_i)/v_i，θ為地震波傳播入射角度，ΔZ_i為第i層的厚度，ΔZ_i＝Z_i+1-Z_i，d表示下行波射線路徑，u表示上行波的射線路徑，其中，針對轉換波的下行波為縱波，上行波為橫波，τ^d表示檢波器端的縱波近地表效應，τ^u表示檢波器端的橫波近地表效應，i表示第i個與第i+1個地層的分界面，n表示第n個地層的層數。

4.根據權利要求3所述的射線參數域轉換波靜校正方法，其中，通過如下公式表示檢波器端的橫波近地表效應τ^u和檢波器端的橫波近地表靜校正量Δτ^u：

其中，v₀表示第0層橫波速度，v₁表示第1層橫波速度，θ為入射角，ΔZ_i為第i層的厚度，ΔZ_i＝Z_i+1-Z_i表示地層厚度，Z₂表示第二層地層的深度，Z₁表示第一層地層的深度；

其中，垂直慢度q_i與射線參數p的關系表達式為：

將公式(5)代入公式(4)中得到其中，射線參數p可以通過如下公式表示：

式中，v_p,i與v_s,i為第i層的P波與S波速度，θ_p,i與θ_s,i為P波與S波的的傳播角度。

5.根據權利要求1所述的射線參數域轉換波靜校正方法，其中，在步驟4)中，對互相關結果進行譜白化與邊界處理，獲得靜校正量道集。

6.根據權利要求1所述的射線參數域轉換波靜校正方法，其中，在步驟3)中，利用高斯低通濾波對共射線參數道集進行平滑處理。