[發明專利]基于瞬時地震子波的疊前域估算地層Q值的方法

在審

• 申請號: 201710451042.0
• 申請日: 2017-06-15
• 公開（公告）號: CN107462923A
• 公開（公告）日: 2017-12-12
• 發明（設計）人: 王延光;劉浩杰;陳雨茂;李民龍;鄭靜靜;趙茂強;王蓬
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 東營雙橋專利代理有限責任公司37107
• 代理人: 方圓

1.基于瞬時地震子波的疊前域估算地層Q值的方法，其特征在于，所述方法包括如下步驟：

步驟1：將子波假定為脈沖子波，將振幅譜假定為光滑譜，設定初始地震子波模型；

步驟2：將測井數據帶入地震頻帶的巖石物理模型；計算得到初始Q值的范圍；然后將初始Q值帶入步驟1得到的初始地震子波模型，獲得初始的瞬時地震子波集；

步驟3：根據步驟2計算所得的初始地震子波集，采用四階累積量匹配法將地震子波譜從地震記錄振幅譜中估計出來，獲得最優的地層擬Q值；

步驟4：利用地層擬Q值與偏移距平方之間的線性關系，將不同偏移距處的地層Q值線性回歸，獲得零偏移距處的地層Q值，即地層常Q值。

2.根據權利要求1所述的基于瞬時地震子波的疊前域估算地層Q值的方法，其特征在于，所述步驟1設定初始地震子波模型的具體過程可描述為：

首先，在常Q介質模型中，確定地震波U(ω,z)在一維黏彈性介質中沿z軸的增加方向以角頻率ω傳播時的傳播方程:

其中式(1)中，ω是角頻率；z是傳播距離；c(ω)是相速度，Q為地層衰減參數，U(ω,0)是深度為0時的地震波；

然后，將震源設定為理想脈沖源，即|U(ω,0)|＝1，并且忽略由散射引起的衰減，因此將式(1)變為：

其中式(2)中，相速度c(ω)滿足：

其中式(3)中，c₀為角頻率ω₀處的相速度；

令t＝z/c(ω)，則將式(2)變為：

然后，假定子波為脈沖子波，假定振幅譜為光滑譜，匹配生成初始地震子波模型：

其中式(5)中，A為幅值參數，Wa(w,t)表示不同時間的地震子波模型。

3.根據權利要求2所述的基于瞬時地震子波的疊前域估算地層Q值的方法，其特征在于，所述步驟2的地震頻帶巖石物理模型為：

其中：E_D為骨架的平面波模量，f為Gassmann方程流體項，為含氣飽和度，f_m為衰減主頻；

f、f_m、E_D均由測井數據確定獲得，將其帶入公式(6)計算初始Q值的范圍，再將計算結果帶入公式(5)獲得初始的地震子波。