[發明專利]一種反射地震記錄頻變速度分析方法

有效

• 申請號: 201811546192.0
• 申請日: 2018-12-18
• 公開（公告）號: CN109581500B
• 公開（公告）日: 2020-06-30
• 發明（設計）人: 李紅星;倪然;陶春輝;劉財;韓立國;徐文斌
• 申請（專利權）人: 東華理工大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 南昌新天下專利商標代理有限公司 36115
• 代理人: 李炳生

1.一種反射地震記錄頻變速度分析方法，其特征在于，包括如下步驟：

a、利用廣義S變換將地震記錄由空間-時間域變換到空間-時間-頻率域，在頻率切片求取時頻譜極大值對應時間，提取不同頻率地震波的反射時距曲線；

b、進行速度掃描，選取計算的理論時距曲線和地震記錄時頻譜中提取的時距曲線擬合精度最高的速度，即使用速度掃描，計算理論時距曲線與提取的時距曲線差的2范數，范數的極小值對應的速度即為某頻率地震波傳播速度；

理論時距曲線表示為：

式中：Tt_Fj(x_i，F)表示頻率為F的地震波在炮檢距為x_i的第i道檢波器的旅行時，t_0Fj表示頻率為F的地震波在第j個反射層的自激自收旅行時，V_0Fj為頻率為F的地震波在第j個反射層處的傳播速度；

理論時距曲線與提取的時距曲線差的2范數表示為：||Tt_Fj-Tp_Fj||₂，范數的極小值對應的速度即為某頻率地震波傳播速度：

V_jF＝find(||Tt_Fj-Tp_Fj||₂＝min(||Tt_Fj-Tp_Fj||₂))

上式表示通過速度掃描，找出使得||Tt_Fj-Tp_Fj||₂取最小值的速度，即為某頻率地震波傳播速度。

2.根據權利要求1所述的一種反射地震記錄頻變速度分析方法，其特征在于，采用廣義S變換時頻分析技術，反射地震記錄u(x，t)的廣義S變換為：

這里u(x，t)為地震記錄，x表示地震記錄道空間，t表示地震記錄時間，U(x，τ，f)為空間-時間-頻率域地震記錄，τ為地震記錄時頻譜時間，f為頻率，λ、p 為廣義S變換調節參數，exp表示e指數。

3.根據權利要求1或2所述的一種反射地震記錄頻變速度分析方法，其特征在于，通過空間-時間-頻率譜的分析，提取能譜極大值位置，獲得不同頻率的反射波時距關系，具體步驟如下：

步驟1：在時間域中輸入空間-時間域二維地震數據，然后采用廣義S變換，將重數據從空間-時間域變換到空間-時間-頻率域；

步驟2：基于空間-時間-頻率譜的頻率切片，提取能譜極大值對應的不同地震道的該頻率的地震波履行時間；

[x_i，Tp_F]＝find(U(x_i，t，F)＝max(U(x_i，t，F)))

[x_i，Tp_F]表示頻率為F的地震波在第i地震道的反射到達時為Tp_F，也即頻率為F Hz的地震波的反射時距關系，x_i為第i道記錄的炮檢距，U(x_i，t，F)為地震記錄空間-時間-頻率譜在F Hz頻率處的切片。