[發(fā)明專利]一種確定地層巖性和孔隙流體的方法

有效

• 申請(qǐng)?zhí)枺?/td>: 200710119461.0
• 申請(qǐng)日: 2007-07-25
• 公開（公告）號(hào): CN101354444A
• 公開（公告）日: 2009-01-28
• 發(fā)明（設(shè)計(jì)）人: 張紅英;侯愛源
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油天然氣集團(tuán)公司;中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司
• 主分類號(hào): G01V1/40
• 分類號(hào): G01V1/40;G01V1/28;G01V1/30;G01V1/36
• 代理公司: 北京市中實(shí)友知識(shí)產(chǎn)權(quán)代理有限責(zé)任公司
• 代理人: 劉天語(yǔ)

1.一種確定地層巖性和孔隙流體的方法，其特征在于：采用以下具體步驟實(shí)現(xiàn)：

1)在地面人工激發(fā)地震，在測(cè)井井中不同深度采樣點(diǎn)放置三分量地震波接收器，記錄三分量VSP地震波波場(chǎng)；

2)對(duì)三分量VSP地震數(shù)據(jù)進(jìn)行波場(chǎng)分離、走廊疊加和頻帶拓寬，得到零井源距地震波場(chǎng)；

3)利用步驟2)處理得到的零井源距地震波場(chǎng)，用下式計(jì)算不同深度采樣點(diǎn)的地震波波場(chǎng)的透射縱波：

透射縱波＝入射縱波+反射縱波+轉(zhuǎn)換橫波-透射轉(zhuǎn)換橫波

4)根據(jù)下式計(jì)算縱橫波速度比：

Φ(d)=(Σk,l=13|mk,l(d)|2+γ1Σm=13|r0-dmrm|2)---(1)]]>

上式中：d是速度比向量，v_p1、v_p2、v_s1、v_s2分別是界面上、下兩側(cè)介質(zhì)的縱、橫波速度；m_k，l是位移梯度矩陣的二階子式，r₀是界面的法線方向矢量；r_m是各類型波的傳播方向矢量和界面的法線方向矢量的內(nèi)積，當(dāng)m＝0時(shí)為下行P波，當(dāng)m＝1時(shí)為上行P波，當(dāng)m＝2時(shí)為上行S波，當(dāng)m＝3時(shí)為下行S波；γ₁是一權(quán)系數(shù)；

令Ф(d)→min，通過Newton迭代法我們可求得界面上、下兩側(cè)介質(zhì)的縱、橫波速度比

5)根據(jù)以下公式計(jì)算縱、橫波反射系數(shù)：

ρ=ρ2ρ1---(2)]]>

Rpp=d1-ρd1+ρ---(3)]]>

Rss=d3-ρd2d3+ρd2---(4)]]>

上式中，為步驟4)中求得縱橫波速度比，R_pp為縱波反射系數(shù)、R_ss為橫波反射系數(shù)，ρ₁為反射界面上覆層的密度，ρ₂為反射界面下伏層的密度；

6)利用以下公式計(jì)算縱、橫波波阻抗：

Imppp=1+Rpp1-Rpp---(5)]]>

Impss=1+Rss1-Rss---(6)]]>

上式中，Im?p_pp為反射界面下伏層的縱波阻抗，Im?p_ss為反射界面下伏層的橫波波阻抗，R_pp為縱波反射系數(shù)、R_ss為橫波反射系數(shù)，

7)根據(jù)步驟6)得到的地層波阻抗變化曲線，對(duì)比已知的測(cè)井資料就可確定地層的巖性和孔隙流體性質(zhì)。