[發明專利]一種基伸縮調頻同步提取地震儲層預測方法

在審

• 申請號: 202110558646.1
• 申請日: 2021-05-21
• 公開（公告）號: CN113296155A
• 公開（公告）日: 2021-08-24
• 發明（設計）人: 胡英;朱冰雪;陳輝;陳旭平;李蕊;周懷來;方玉霞
• 申請（專利權）人: 成都理工大學
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/36
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種基伸縮調頻同步提取地震儲層預測方法，其特征在于，包括以下步驟：

S1、輸入待分析的原始二維地震剖面信號s(x；t)；其中x表示地震道，t表示時間；

S2、對S1輸入的二維地震剖面信號s(x；t)進行基伸縮自適應調頻變換，獲取時頻變換結果SBCT(x；f,t)，計算方法如下：

其中，t表示時間中心，f表示頻率中心，h(x；τ)表示高斯函數，表示相位函數，τ為時間變量；

S3、根據S2得到的時頻變換結果SBCT(x；f,t)估計第x道地震信號在各個時頻位置(t,f)處的瞬時頻率

S4、根據同步提取原理，在該時頻域上構建一種以信號瞬時頻率曲線為中心的同步提取算子SEO(x；f,t)用于提取原時頻譜能量，獲得新的時頻系數T(x；f,t)，剔除掉模糊的時頻能量；

S5、對S4獲得的時頻系數取模，得到基伸縮調頻同步提取變換后的時頻譜；

S6、采用傅氏變換計算原始地震信號s(x；t)的主頻范圍，將主頻范圍中的頻率最大值的85％設置為高頻值，將主頻范圍中的頻率最小值設置為低頻值，然后從S5得到的時頻譜中抽取高頻值對應的頻譜，得到高頻值的共頻率剖面；從S5得到的時頻譜中抽取低頻值對應的頻譜，得到低頻值的共頻率剖面。通過對比兩個共頻率剖面中地震信號的衰減情況來確定該地震剖面中是否存在儲層。

2.根據權利要求1所述的基伸縮調頻同步提取地震儲層預測方法，其特征在于，所述步驟S2中二維地震剖面信號s(x；t)的基伸縮調頻自適應變換的時頻變換結果為：

其中，t表示時間中心，f表示頻率中心，h(x；τ)表示高斯函數，表示相位函數，τ為時間變量，所述相位函數為：

其中，k＝1,2,…,n，n為相位函數的階數，(a₁,a₂,…,a_n)是相位函數的自適應參數，可根據峰度理論確定，即：

其中，argmax(·)表示函數取到最大值時的參數集。

3.根據權利要求1所述的基伸縮調頻同步提取地震儲層預測方法，其特征在于，所述步驟S3中瞬時頻率的計算包括根據S2得到的時頻變換結果SBCT(x；f,t)估計信號在各個時頻位置(t,f)處的瞬時頻率

其中，表示窗函數h(x；t)對時間t的偏導，表示在窗函數下的SBCT結果；表示在函數t^kh(x；t)下的SBCT結果。

4.根據權利要求1所述的基伸縮調頻同步提取地震儲層預測方法，其特征在于，所述步驟S4中獲取新的時頻系數T(x；f,t)的具體方法為：根據同步提取原理，在該時頻域上構建一種以信號瞬時頻率為曲線中心的同步提取算子SEO(x；f,t)提取原時頻譜能量獲得新的時頻系數T(x；f,t)，剔除掉模糊的時頻能量，同步提取算子SEO(x；f,t)滿足：

其中δ(·)為單位脈沖函數，新的時頻系數T(x；f,t)為：

T(x；f,t)＝SBCT(x；f,t)SEO(x；f,t) (7) 。

5.根據權利要求1所述的基伸縮同步提取地震儲層預測方法，其特征在于，所述步驟S5可采用下式對(7)中的T(x；f,t)進行逆變換，重構地震信號s(x；t)；

其中，本發明適用于重構地震信號。