[發明專利]下伏灰質地層干涉的薄層厚度定量預測方法

在審

• 申請號: 202011420611.3
• 申請日: 2020-12-07
• 公開（公告）號: CN114594514A
• 公開（公告）日: 2022-06-07
• 發明（設計）人: 宋亮;蘇朝光;巴素玉;張盼盼;王楠;張秀娟;張鵬
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/36
• 代理公司: 東營雙橋專利代理有限責任公司 37107
• 代理人: 侯華頌

1.一種下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于包括以下步驟：

步驟110，地震正演模擬，定量分析下伏灰質地層對薄層砂巖的地震干涉特征；

步驟120，通過實鉆砂體厚度統計，確定所述薄層砂巖的優勢調諧頻率；

步驟130，地震分頻處理，提取所述優勢調諧頻率下的最佳地層切片組合；

步驟140，實鉆井厚度約束，建立分頻切片融合公式；

步驟150，利用所述分頻切片融合公式，對灰質干涉下的薄層厚度進行定量預測。

2.根據權利要求1所述的一種下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于，步驟110中所述地震正演模擬，定量分析下伏灰質地層對薄層的地震干涉特征，具體包括：

通過實鉆井分析，明確下伏灰質地層及薄層砂巖的厚度、速度及密度的巖石物理特征，建立下伏灰質地層的薄層砂巖地震正演模型；

通過波動方程地震正演模擬，確定下伏灰質地層對薄層砂巖開始產生干涉時的地震響應特征，并定量統計薄層砂巖厚度與不同位置地層切片振幅的變化關系。

3.根據權利要求2所述的一種下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于，進一步包括：選擇能較好反映薄層砂巖厚度變化的地層切片組合，通過相鄰波峰振幅比法處理，即下伏灰質地層波峰切片振幅與砂巖波峰切片振幅之比，統計薄層砂巖厚度與不同位置地層切片振幅的變化關系。

4.根據權利要求2所述的下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于，在步驟120中所述通過實鉆砂體厚度統計，確定所述薄層砂巖的優勢調諧頻率，具體包括：

按照厚度調諧原理，在地層速度確定的前提下，當砂巖厚度值等于四分之一波長時，根據地震縱波傳播速度公式，求取對應的頻率：

f＝v/λ＝v/4h

其中，v為縱波速度，f為地震資料主頻，λ為地震縱波波長，h為調諧厚度，令h等于目標砂體厚度，利用上式計算得到f即為反映目標砂體的優勢調諧頻率。

5.根據權利要求3所述的下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于，在步驟130所述地震分頻處理，提取所述優勢調諧頻率下的最佳地層切片組合，具體包括：在對地震數據進行分頻處理，選擇S變換作為時頻分析算法，根據如下公式進行計算：

其中，t表示時間，f表示頻率，τ為高斯窗在時間軸上的位置，時窗寬度隨頻率的變化而改變，頻率的倒數決定時窗的大小，s(t)表示時間域輸入信號，exp為自然對數函數，π為圓周率，i表示虛數單位，dt表示積分函數。

對地震數據進行分頻處理后，生成優勢的調諧頻率體，并提取目標砂體地層切片的振幅屬性，獲得所述優勢調諧頻率下的最佳地層切片組合。

6.根據權利要求3所述的下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于，在步驟130所述地震分頻處理具體包括：采用傅里葉變換、小波變換時頻分析算法，對地震數據進行分頻處理。

7.根據權利要求5所述的下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于，所述提取目標砂體地層切片的振幅屬性具體包括：提取最佳地層切片組合的相鄰波峰振幅比，即下伏灰質地層波峰切片振幅與砂巖波峰切片振幅之比。

8.根據權利要求4所述的下伏灰質地層干涉的薄層厚度定量預測方法，其特征在于，在步驟140中所述實鉆井厚度約束，建立分頻切片融合公式，具體包括以下步驟：

對優勢頻率范圍內的地層切片進行所述相鄰波峰振幅比法處理，依據實鉆井砂體厚度與分頻切片的統計擬合關系，將比值后的分頻切片進行加權線性回歸，建立分頻切片融合公式，以求得目標砂體的厚度值：

Thickness＝A₁k₁+A₂k₂+A₃k₃+……+A_nk_n+C

其中，Thickness為預測的薄層砂巖厚度；A為優勢調諧頻率的最佳地層切片組合即相鄰波峰振幅比，A₁為第1個優勢頻率的地層切片組合，A_n為第n個優勢頻率的地層切片組合；k₁為第1個地層切片組合的加權系數值，k_n為第n個地層切片組合的加權系數值；C為加權常數。