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

有效

• 申請(qǐng)?zhí)枺?/td>: 201811546192.0
• 申請(qǐng)日: 2018-12-18
• 公開(kāi)（公告）號(hào): CN109581500B
• 公開(kāi)（公告）日: 2020-06-30
• 發(fā)明（設(shè)計(jì)）人: 李紅星;倪然;陶春輝;劉財(cái);韓立國(guó);徐文斌
• 申請(qǐng)（專利權(quán)）人: 東華理工大學(xué)
• 主分類號(hào): G01V1/30
• 分類號(hào): G01V1/30
• 代理公司: 南昌新天下專利商標(biāo)代理有限公司 36115
• 代理人: 李炳生

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

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

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

理論時(shí)距曲線表示為：

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

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

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

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

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

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

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

步驟1：在時(shí)間域中輸入空間-時(shí)間域二維地震數(shù)據(jù)，然后采用廣義S變換，將重?cái)?shù)據(jù)從空間-時(shí)間域變換到空間-時(shí)間-頻率域；

步驟2：基于空間-時(shí)間-頻率譜的頻率切片，提取能譜極大值對(duì)應(yīng)的不同地震道的該頻率的地震波履行時(shí)間；

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

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