一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法
【專利摘要】一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，采集測(cè)井?dāng)?shù)據(jù)，建立隨壓力變化的孔隙介質(zhì)橫波速度預(yù)測(cè)模型；將校正過(guò)的各個(gè)參數(shù)代入隨壓力變化的孔隙介質(zhì)橫波速度預(yù)測(cè)模型得到孔隙介質(zhì)的與壓力相關(guān)的縱橫波速度；根據(jù)孔隙介質(zhì)下的與壓力相關(guān)的縱波速度、橫波速度得到裂隙介質(zhì)對(duì)稱方向的與壓力相關(guān)的縱波速度、橫波速度并代入Thomson公式中，得到平行于裂隙方向的與壓力相關(guān)的橫波速度和垂直于裂隙方向的與壓力相關(guān)的橫波速度。本發(fā)明既考慮了壓力對(duì)于介質(zhì)橫波速度的影響，也考慮到了儲(chǔ)層中含有裂隙，更加符合實(shí)際。根據(jù)本發(fā)明得到的裂隙介質(zhì)橫波速度，可以更好的評(píng)估油氣藏巖石屬性，包括孔隙度、密度、巖性與流體含量，可以建立AVO模。
【專利說(shuō)明】一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于地震勘探中巖石物理領(lǐng)域，具體涉及一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法。
【背景技術(shù)】
[0002]人們通常用縱橫波速度比為常數(shù)來(lái)代替橫波速度，但是對(duì)于不同的介質(zhì)縱橫波速度比應(yīng)該是不同的，這樣的假設(shè)不合理。因此很多橫波速度的預(yù)測(cè)方法隨之而出，有經(jīng)驗(yàn)公式，也有基于巖石物理理論的。
[0003]比較著名的經(jīng)驗(yàn)公式包括Castagna等人在1985年提出的著名的泥巖線為Vs=0.862Vp-l.172。Gardner于1974年給出了不同的巖性之間的速度與密度的關(guān)系，其中它的平均變換式為P =0.23V°_25，這個(gè)平均變換式是對(duì)所有巖性的速度與密度關(guān)系的最佳擬合，它適合于所有巖性，不僅僅適用用某種巖性。而Castagna于1993年又對(duì)Gardner的
公式進(jìn)行了擴(kuò)充，得到了不同巖性的速度與密度之間的關(guān)系:對(duì)于砂巖有
【權(quán)利要求】
1.一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，包括以下步驟: I)采集測(cè)井?dāng)?shù)據(jù):采集總孔隙度Φζ、巖石的體積密度P以及實(shí)際縱波速度Vp、流體飽和度；采集巖石骨架的體變模量Kma，巖石骨架的切變模量μΜ，差異壓力P，測(cè)量配位數(shù)Cp’，巖石顆粒變形之前接觸區(qū)域的半徑a與巖石顆粒的半徑R ； 分別對(duì)總孔隙度Φζ、巖石的體積密度P進(jìn)行校正得到孔隙介質(zhì)的孔隙孔隙度Φρ、孔隙介質(zhì)的體積密度P ρ，利用流體飽和度得到混合流體的體變模量Kf ;計(jì)算巖石骨架的泊松比 Vma; 對(duì)測(cè)量配位數(shù)Cp’進(jìn)行線性擬合并且加權(quán)，得到含有加權(quán)系數(shù)W的配位數(shù)Cp，然后利用Digby> Mindlin公式以及Gassmann方程的變形公式建立隨壓力變化的孔隙介質(zhì)橫波速度預(yù)測(cè)模型； 2)對(duì)巖心進(jìn)行巖石物理測(cè)試，得到巖心的縱波速度、巖心的橫波速度、巖心的密度以及巖心的孔隙度，然后利用隨壓力變化的孔隙介質(zhì)橫波速度預(yù)測(cè)模型對(duì)巖心的縱波速度、巖心的橫波速度、巖心的密度以及巖心的孔隙度進(jìn)行計(jì)算得到巖心的加權(quán)系數(shù)Wc^，然后利用巖心的加權(quán)系數(shù)對(duì)加權(quán)系數(shù)W進(jìn)行校正，得到校正后的各向同性介質(zhì)的加權(quán)系數(shù)Wisotropy ，如公式(14)所示:


2.一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟I)中若取巖心的井與進(jìn)行采集測(cè)井?dāng)?shù)據(jù)的井是兩口井，則標(biāo)定巖心在進(jìn)行采集數(shù)據(jù)的井中的層位。
3.根據(jù)權(quán)利要求1所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟I)中隨壓力變化的孔隙介質(zhì)橫波速度預(yù)測(cè)模型建立的具體過(guò)程為: a)對(duì)測(cè)量配位數(shù)Cp’進(jìn)行線性擬合，得到測(cè)量配位數(shù)Cp’與孔隙度Φ的關(guān)系為:CP’=1 1.759θ1_φ-12.748(I) 對(duì)測(cè)量配位數(shù)Cp’進(jìn)行加權(quán)，得到配位數(shù)Cp，如公式(2)所示， Cp=ff*C；(2) b)在均勻的孔隙介質(zhì)的條件下，利用Gassmann方程的變形進(jìn)行縱橫波速度的預(yù)測(cè)；公式(3)- (5)為Gassmann方程的變形為:
4.根據(jù)權(quán)利要求3所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟b)中Digby公式為:
5.根據(jù)權(quán)利要求3所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟b)和步驟c)中Mindlin公式為:
6.根據(jù)權(quán)利要求3所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟b)中
7.根據(jù)權(quán)利要求1所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟3)中將總孔隙度Φζ校正為孔隙介質(zhì)下的孔隙孔隙度Φρ的具體過(guò)程為:將總孔隙度02減去裂隙孔隙度Φ。得到孔隙孔隙度Φρ。
8.根據(jù)權(quán)利要求1所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟I)中對(duì)巖石的體積密度P進(jìn)行校正，得到孔隙介質(zhì)情況下的孔隙介質(zhì)的密度P ρ的具體過(guò)程如下: 首先利用公式Pd=Pi(1-Ctz)計(jì)算得到干巖石的密度，其中Pd為干巖石的密度，然后利用公式ρρ=ρ/(1-φρ)得到孔隙介質(zhì)的密度ρρ。
9.根據(jù)權(quán)利要求1所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟4)中對(duì)稱方向的彈性介質(zhì)的縱波速度Citl、橫波速度β C1是通過(guò)以下過(guò)程得到的: 根據(jù)公式(18)和(19)，計(jì)算得到對(duì)稱方向的彈性介質(zhì)的縱波速度Citl、橫波速度
10.根據(jù)權(quán)利要求1所述的一種隨壓力變化的裂隙介質(zhì)橫波速度預(yù)測(cè)方法，其特征在于，所述步驟4)中各向異性參數(shù)ε、Y和δ的表達(dá)式分別為:
【文檔編號(hào)】G01V1/28GK103576195SQ201310518774
【公開(kāi)日】2014年2月12日 申請(qǐng)日期:2013年10月28日 優(yōu)先權(quán)日:2013年10月28日
【發(fā)明者】李琳, 馬勁風(fēng) 申請(qǐng)人:西北大學(xué)