[發明專利]一種利用地震流體阻抗進行流體確定的方法

有效

• 申請號: 200710178116.4
• 申請日: 2007-11-27
• 公開（公告）號: CN101446645A
• 公開（公告）日: 2009-06-03
• 發明（設計）人: 李紅兵;崔興福
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/40
• 代理公司: 北京市中實友知識產權代理有限責任公司
• 代理人: 劉天語

1.一種利用地震流體阻抗進行流體確定的方法，其特征在于通過以下具體步驟實現：

1)激發并記錄地震波，按常規地震資料進行高保真處理，形成用于振幅隨偏移距變化分析的正常時差校正后的共中心點CMP道集；

2)利用下式將步驟1)中的共中心點CMP道集轉換到角度域中，得到角度道集，

sinθ=xvpvrms2(t02+x2vrms2)1/2---(1)]]>

式中，x為炮檢距，V_p為地層的層速度，t₀為零偏移距傳播時間，V_rms為地層的均方根速度；

3)對角度道集作頻率補償以消除由于動校正拉伸所導致的信號頻率畸變；以下公式校正不同偏移距的振幅變化，使遠道的統計平均振幅譜匹配零偏移距道，

αx=∂t∂t0≈t0t(1-2(t-t0)VNMO∂VNMO∂t0)---(2)]]>

式中α_x為壓縮系數，t₀為零偏移距傳播時間，t為傳播時間，V_NMO為動校正均方根速度；

以下公式消除遠偏移距和由于動校正拉伸所造成的降頻，

wθ(t)=w(tcosθ)→FTw(ω·cosθ)---(3)]]>

式中w_θ為拉伸后子波，w為零偏移距子波，θ為反射角，t為傳播時間，ω為角頻；

4)利用下列的Fatti公式從補償后的角度道集中提取縱波反射率屬性和加權橫波反射率屬性

(4)式中的兩個屬性的提取過程求解下面的方程組：

ΔIpI‾p4Vs2Vp2ΔIsI‾s1+tan2θ12-sin2θ11+tan2θ22-sin2θ2······1+tan2θn2-sin2θn-1R(θ1)R(θ2)···R(θn)---(5)]]>

式中θ_i，i＝1，2，…，n為步驟3)中得到的角度道集地震數據的入射角，R(θ_i)為步驟3)中得到的角度道集P波反射地震數據，V_p、V_s及分別是反射層面上下層的P波和S波的平均速度值及平均阻抗值，ΔI_p及ΔI_s分別是反射層面上下層的P波、S波的變化量；

5)對以上縱波反射率屬性和加權橫波反射率屬性利用下式計算改進后的流體因子：

Rf=ΔIpI‾p-c4Vs2Vp2ΔIsI‾s---(6)]]>

式中R_f為流體反射率，也稱為流體因子，c是常數為2-3；

6)采用常規方法取得測井資料，處理后從縱、橫波和密度曲線中取得流體阻抗曲線；

7)用步驟5)提取的流體因子和步驟6)提取的流體阻抗曲線進行迭代反演，得到流體阻抗剖面；

迭代反演算法是：

(1)利用流體阻抗曲線計算流體因子R_f(t)，再把它與地震子波w(t)褶積，得到合成流體因子記錄d(t)＝R_f(t)*w(t)；

(2)把合成流體因子記錄與從地震資料中提取的流體因子的井旁道進行標定，

(3)采用通常距離加權的插值方法由測井流體阻抗建立初始的流體阻抗模型；

(4)從流體因子剖面出發，得到最終的流體阻抗剖面，在流體阻抗剖面上進行流體分布的橫向追蹤；

8)利用流體因子屬性定性分析孔隙中的流體異常的變化，利用流體阻抗定量描述流體的空間分布。