[發明專利]碳酸鹽巖溶洞儲集體發育規模地震識別方法及系統

有效

• 申請號: 201710725084.9
• 申請日: 2017-08-22
• 公開（公告）號: CN109425890B
• 公開（公告）日: 2020-09-15
• 發明（設計）人: 胡華鋒;馬靈偉;楊麗;呂慧;胡起
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種碳酸鹽巖溶洞儲集體發育規模地震識別方法，其特征在于，該方法包括：

基于疊后地震數據提取各分頻能量地震數據體；

基于各分頻能量地震數據體提取頻率梯度屬性；

基于頻率梯度屬性實現碳酸鹽巖溶洞儲集體發育規模預測；

其中，基于疊后時間域地震數據Seis(x,y,t)，獲得n個不同頻率地震數據體

Seis_f(f₁,x,y,t)、Seis_f(f₂,x,y,t)、…、 Seis_f(f_n,x,y,t)，

其中，x,y代表每個地震道對應的空間位置，t代表時間域，Seis_f(f₁,x,y,t)為頻率為f₁的分頻能量地震數據體，Seis_f(f₂,x,y,t)為頻率為f₂的分頻能量地震數據體，Seis_f(f_n,x,y,t)為頻率為f_n的分頻能量地震數據體，并且滿足頻率f₁至頻率f_n均分布在主頻帶范圍內，即f_n＝f₁+(n-1)*Δf，Δf為頻率步長；

采用最小二乘擬合計算頻率梯度屬性，擬合各單頻能量隨頻率的變化的規律，擬合公式為：

y＝a*x+b

其中，x＝{f₁,f₂,…,f_n}，y＝{Seis_f(f₁,x,y,t),Seis_f(f₂,x,y,t),…,Seis_f(f_n,x,y,t)}，擬合系數a即為頻率的能量隨頻率變化的梯度屬性，b為常數。

2.根據權利要求1所述的碳酸鹽巖溶洞儲集體發育規模地震識別方法，其特征在于，利用傅里葉變換或者匹配追蹤方法提取各分頻能量地震數據體。

3.根據權利要求1所述的碳酸鹽巖溶洞儲集體發育規模地震識別方法，其特征在于，通過已知井建立溶洞規模同頻率梯度的擬合關系，在獲得頻率梯度屬性體中，頻率梯度與溶洞規模呈負相關關系。

4.一種碳酸鹽巖溶洞儲集體發育規模地震識別系統，其特征在于，該系統包括：

存儲器，存儲有計算機可執行指令；

處理器，所述處理器運行所述存儲器中的計算機可執行指令，執行以下步驟：

基于疊后地震數據提取各分頻能量地震數據體；

基于各分頻能量地震數據體提取頻率梯度屬性；

基于頻率梯度屬性實現碳酸鹽巖溶洞儲集體發育規模預測；

其中，基于疊后時間域地震數據Seis(x,y,t)，獲得n個不同頻率地震數據體

Seis_f(f₁,x,y,t)、Seis_f(f₂,x,y,t)、…、 Seis_f(f_n,x,y,t)，

其中，x,y代表每個地震道對應的空間位置，t代表時間域，Seis_f(f₁,x,y,t)為頻率為f₁的分頻能量地震數據體，Seis_f(f₂,x,y,t)為頻率為f₂的分頻能量地震數據體，Seis_f(f_n,x,y,t)為頻率為f_n的分頻能量地震數據體，并且滿足頻率f₁至頻率f_n均分布在主頻帶范圍內，即f_n＝f₁+(n-1)*Δf，Δf為頻率步長；

采用最小二乘擬合計算頻率梯度屬性，擬合各單頻能量隨頻率的變化的規律，擬合公式為：

y＝a*x+b

其中，x＝{f₁,f₂,…,f_n}，y＝{Seis_f(f₁,x,y,t),Seis_f(f₂,x,y,t),…,Seis_f(f_n,x,y,t)}，擬合系數a即為頻率的能量隨頻率變化的梯度屬性，b為常數。