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

有效

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

本發明涉及一種反射地震記錄頻變速度分析方法，包括如下步驟：a、利用廣義S變換將地震記錄由空間?時間域變換到空間?時間?頻率域，在頻率切片求取時頻譜極大值對應時間，提取不同頻率地震波的反射時距曲線；b、進行速度掃描，選取計算的理論時距曲線和地震記錄時頻譜中提取的時距曲線擬合精度最高的速度，即認為是某頻率地震波的傳播速度。本發明提出了新的反射地震記錄頻變速度分析方法，克服了傳統速度分析方法無法獲得速度頻散屬性的缺點，頻變速度分析技術可以利用速度頻散屬性進行油氣儲層預測，基于頻變速度分析技術可以提高地震處理速度分析、動校正和疊加的精度。

技術領域

本發明涉及地震記錄的速度分析方法，尤其是一種反射地震記錄頻變速度分析方法。

背景技術

在油氣儲層、礦產資源、工程地質及環境地質反射地震勘探中，地震速度分析技術是反射地震數據處理的重要和必須環節。速度分析質量的優劣直接影響反射地震成像精度。同時，地震波傳播速度作為描述地質體的重要地震屬性，也被用來對巖性，尤其時儲層性質進行定性和定量解釋。

儲層和海底介質等是典型的孔隙多相介質，地震波在儲層和海底介質中的傳播速度從在頻散已經得到實踐證明(基于儲層砂巖微觀孔隙結構特征的彈性波頻散響應分析[J].地球物理學報,2015,58(9):3389-3400.)，地震波速度頻散成為儲層預測越來越被重視的屬性。但是，目前的反射地震速度分析技術研究尚未有報道頻變速度分析技術，學者們都是在不同復雜程度的地質模型基礎上或不同噪聲條件下探討將所有頻率地震波傳播速度認為是同一速度的速度分析方法。不依賴頻率的速度分析方法獲得的速度參數不隨頻率變化，會導致反射地震資料處理中的動校正和疊加不能取得更好的效果，會導致無法利用速度分析技術獲得速度參數的頻散屬性，從而預測油氣儲層。

發明內容

本發明的目的是為了解決現有技術存在的缺陷，提供一種反射地震記錄頻變速度分析方法。

為了實現上述目的，本發明采用的技術方案是：

一種反射地震記錄頻變速度分析方法，包括如下步驟：

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

b、進行速度掃描，選取計算的理論時距曲線和地震記錄時頻譜中提取的時距曲線擬合精度最高的速度，即認為是某頻率地震波的傳播速度。

進一步，采用廣義S變換時頻分析技術，反射地震記錄u(x,t)的廣義S變換為：

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

進一步，通過空間-時間-頻率譜的分析，提取能譜極大值位置，獲得不同頻率的反射波時距關系，具體步驟如下：

步驟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，也即頻率為FHz的地震波的反射時距關系，x_i為第i道記錄的炮檢距，U(x_i,t,F)為地震記錄空間-時間-頻率譜在FHz頻率處的切片。