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


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