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


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