一種地球物理勘探橫波估算方法和系統(tǒng)的制作方法
【專利摘要】本發(fā)明公開了一種地球物理勘探橫波估算方法和系統(tǒng)。該方法包括根據(jù)測(cè)井曲線分析儲(chǔ)層的巖性特征和流體特征；根據(jù)儲(chǔ)層的巖性特征和流體特征分別確定相應(yīng)的巖性物理模量參數(shù)和流體物理模量參數(shù)；建立儲(chǔ)層的巖石物理模型；根據(jù)巖石物理模型估算儲(chǔ)層的縱波速度，分析估算的縱波速度與實(shí)測(cè)的縱波速度之間的誤差，以最小化此誤差為目標(biāo)，求取相應(yīng)的虛擬孔隙度；基于虛擬孔隙度求取儲(chǔ)層含水巖石的橫波速度；通過流體替換求取儲(chǔ)層含油氣巖石的橫波速度；輸出橫波速度曲線。本發(fā)明方法通過引入虛擬孔隙度的概念，重構(gòu)一條虛擬孔隙度曲線用來描述儲(chǔ)層實(shí)際地質(zhì)與理論模型之間的差異，從而避免了現(xiàn)有技術(shù)中必須基于縱橫速度比進(jìn)行地層分析的難題，大大簡化了預(yù)測(cè)流程，并且提高了估算結(jié)果的精度，使得計(jì)算出的橫波速度更為合理和準(zhǔn)確。
【專利說明】一種地球物理勘探橫波估算方法和系統(tǒng)

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及石油地球物理勘探【技術(shù)領(lǐng)域】，尤其涉及一種能夠應(yīng)用于石油地球物理勘探中的儲(chǔ)層精細(xì)描述和流體識(shí)別的橫波估算方法和系統(tǒng)。

【背景技術(shù)】
[0002]隨著油氣勘探開發(fā)的不斷深入，隱蔽性油氣藏的比例越來越高，儲(chǔ)層描述難度也越來越大，因此迫切需要新的技術(shù)對(duì)儲(chǔ)層巖性、物性和含油氣性的準(zhǔn)確描述給予支撐。這其中對(duì)于包含很多信息的疊前地震信息處理解釋方法的研究尤其重要。而疊前地震信息處理解釋方法的關(guān)鍵技術(shù)則是疊前地震反演。通過疊前地震反演技術(shù)，在已知縱波震源地震資料的情況下，利用疊前地震時(shí)間偏移道集資料，不僅可以得到縱波速度，而且可以得到橫波速度和密度參數(shù)，從而可以通過多種反演參數(shù)和彈性參數(shù)對(duì)儲(chǔ)層巖性、物性和流體性質(zhì)進(jìn)行有效的分析解釋。例如，縱橫波速度比是儲(chǔ)層重要的油氣指示，利用它計(jì)算出的泊松比是儲(chǔ)層巖性劃分的重要指標(biāo)；同時(shí)對(duì)于疊前地震彈性反演和AVO(振幅隨偏移距變化)處理，需要已知井的縱、橫波測(cè)井速度作為約束條件才有意義。但是在實(shí)際勘探中，考慮到橫波測(cè)井成本和數(shù)量的限制，絕大多數(shù)井并不做橫波測(cè)井。這就導(dǎo)致疊前地震反演缺失某些約束條件，進(jìn)而在一定程度上影響疊前地震反演結(jié)果的準(zhǔn)確度。因此如何基于有限的測(cè)井資料進(jìn)行巖石物理建模以及在此基礎(chǔ)上完成準(zhǔn)確的橫波估算是疊前地震反演技術(shù)中一個(gè)重要的課題。
[0003]巖石物理建模的主要任務(wù)是研究地球物理勘探所獲得的物理量與地下儲(chǔ)層參數(shù)的對(duì)應(yīng)關(guān)系，以了解巖石物性與地震波特性的關(guān)系，從而獲得各種用于儲(chǔ)層識(shí)別及含油氣性分析的敏感巖石物理參數(shù)，為準(zhǔn)確地預(yù)測(cè)儲(chǔ)層性質(zhì)(巖性、物性)及其狀態(tài)(飽和度、孔隙壓力等)提供有力的支持。
[0004]目前實(shí)際生產(chǎn)應(yīng)用中較為流行的巖石物理建模主要為以下幾種:
[0005]I) Gassmann 方程:
[0006]在研究巖石孔隙以及含流體性質(zhì)對(duì)地震性質(zhì)的影響時(shí)，普遍采用Gassmann方程，基于固體礦物、干骨架及孔隙流體的體積模量計(jì)算流體飽和多孔介質(zhì)的體積模量。在Gassmann方程中，干骨架的彈性模量及密度是指在束縛水狀態(tài)下的彈性模量和體密度。在此基礎(chǔ)上，B1t、Gassmann又發(fā)展了適用于地震勘探頻率范圍的縱橫波速度計(jì)算方法(簡稱B-G方程)。Gassmann和B-G方程中都要涉及介質(zhì)的總孔隙度參數(shù)。除了巖石骨架的彈性性質(zhì)、巖石中流體性質(zhì)以及孔隙度大小以外，孔隙的形狀(即縱橫比)也是影響縱橫波速度的一個(gè)重要因素?？紫兜男螤顚?duì)縱橫波速度的影響主要是通過干骨架體積模量體現(xiàn)。在使用Gassmann和B-G方程模型時(shí)均假設(shè)飽和流體不影響巖石的剪切模量，而這種假設(shè)條件只有在球形孔隙條件下才成立，因此Gassmann和B-G方程實(shí)際上隱含了孔隙形狀為球形的應(yīng)用條件。但是由于實(shí)際的孔隙形狀通常非常復(fù)雜，涵蓋了球形、橢球形、扁球形至裂隙所有譜系，因此基于該方法建立的模型與實(shí)際的地層模型有一定的偏差。
[0007]2)利用縱波速度預(yù)測(cè)橫波速度:
[0008]如利用Castagna泥巖趨勢(shì)線和Greenberg — Castanga砂泥巖公式等。
[0009]Vs = 0.862Vp-l.172
[0010]以及Castagna (1993)的拋物線關(guān)系法
[0011]Vs = AVp2+BVp+C
[0012]Vp為縱波速度；VS為橫波速度；A、B、C為常數(shù)。
[0013]由于使用該方法時(shí)需要憑經(jīng)驗(yàn)估計(jì)很多參數(shù)，因此基于該方法獲得的結(jié)果精度較低。
[0014]3)理論模型法:
[0015]目前主要流行的理論模型法是基于Kuster-TokWz模型和XU-WHITE模型的方法。
[0016]Kuster- ToksSz基于散射理論，導(dǎo)出了長波長一階近似條件下流體飽和孔隙介質(zhì)的彈性模量方程，即Kuster- ToksSz方程(簡稱K-T方程)。與Gassmann方程類似，K-T方程中隱含孔隙形狀為橢球和扁球形的應(yīng)用條件，同時(shí)考慮介質(zhì)的孔隙度及其縱橫比對(duì)縱波速度和橫波速度的影響。因此基于該方法建立的模型與實(shí)際的地層模型也有一定的偏差。
[0017]XU-WHITE提出了一種根據(jù)巖石孔隙度和粘土含量來預(yù)測(cè)聲波速度的方法。該方法從礦物的組成入手，綜合了幾種理論方法，主要分析泥、砂混合狀態(tài)下的縱橫波的速度、密度與孔隙度、泥質(zhì)含量的關(guān)系，并把粘土成分對(duì)聲波的影響歸因于泥巖和砂巖的孔隙幾何形狀和面孔率的差異。其中，由于泥頁巖孔隙的面孔率比砂巖小，因此泥頁巖孔隙更易變形，對(duì)聲波速度有明顯的影響。相比于前面幾種方法，基于該方法獲得的結(jié)果更加準(zhǔn)確。但是由于該方法需要輸入的參數(shù)繁多，計(jì)算過程復(fù)雜，因此整體上效率偏低，并且非測(cè)井專業(yè)人士很難掌握。
[0018]綜上所述，雖然上述幾種方法能夠在一定程度上解決特定地層段的巖石物理建模問題，但是由于所建立的模型與實(shí)際的地層模型有一定的偏差，因此最終獲得的橫波速度估算結(jié)果精度不高。


【發(fā)明內(nèi)容】

[0019]針對(duì)上述問題，本發(fā)明提出了一種基于虛擬孔隙度曲線的橫波估算方法和系統(tǒng)。
[0020]該方法包括以下步驟:
[0021]步驟一，根據(jù)測(cè)井曲線分析儲(chǔ)層的巖性特征和流體特征；
[0022]步驟二，根據(jù)儲(chǔ)層的巖性特征和流體特征分別確定相應(yīng)的巖性物理模量參數(shù)和流體物理模量參數(shù)；
[0023]步驟三，建立儲(chǔ)層的巖石物理模型；
[0024]步驟四，根據(jù)巖石物理模型估算儲(chǔ)層的縱波速度，分析估算的縱波速度與實(shí)測(cè)的縱波速度之間的誤差，以最小化此誤差為目標(biāo)，求取相應(yīng)的虛擬孔隙度；
[0025]步驟五，基于虛擬孔隙度求取儲(chǔ)層含水巖石的橫波速度；
[0026]步驟六，通過流體替換求取儲(chǔ)層含油氣巖石的橫波速度；
[0027]步驟七，輸出橫波速度曲線。
[0028]根據(jù)本發(fā)明的實(shí)施例，上述步驟一中，測(cè)井曲線可以包括伽瑪曲線、密度曲線、縱波速度曲線、橫波速度曲線、孔隙度曲線、滲透率曲線、電阻率曲線和巖性曲線，以及泥質(zhì)含量曲線。
[0029]根據(jù)本發(fā)明的實(shí)施例，上述步驟二中，巖性物理模量參數(shù)可以包括砂巖和泥巖的體積模量、剪切模量和密度，以及孔隙扁度。
[0030]根據(jù)本發(fā)明的實(shí)施例，上述步驟二中，流體物理模量參數(shù)可以包括流體的體積模量、剪切模量和密度。
[0031]根據(jù)本發(fā)明的實(shí)施例，上述步驟三中，根據(jù)以下Xu-White方程模型，計(jì)算儲(chǔ)層的孔隙度，建立巖石物理模型:
[0032]

【權(quán)利要求】
1.一種地球物理勘探橫波估算方法，包括以下步驟: 步驟一，根據(jù)測(cè)井曲線分析儲(chǔ)層的巖性特征和流體特征； 步驟二，根據(jù)儲(chǔ)層的巖性特征和流體特征分別確定相應(yīng)的巖性物理模量參數(shù)和流體物理模量參數(shù)； 步驟三，建立儲(chǔ)層的巖石物理模型； 步驟四，根據(jù)巖石物理模型估算儲(chǔ)層的縱波速度，分析估算的縱波速度與實(shí)測(cè)的縱波速度之間的誤差，以最小化此誤差為目標(biāo)，求取相應(yīng)的虛擬孔隙度； 步驟五，基于虛擬孔隙度求取儲(chǔ)層含水巖石的橫波速度； 步驟六，通過流體替換求取儲(chǔ)層含油氣巖石的橫波速度； 步驟七，輸出橫波速度曲線。
2.根據(jù)權(quán)利要求1所述的橫波估算方法，其特征在于，所述步驟一中，測(cè)井曲線包括伽瑪曲線、密度曲線、縱波速度曲線、橫波速度曲線、孔隙度曲線、滲透率曲線、電阻率曲線和巖性曲線，以及泥質(zhì)含量曲線。
3.根據(jù)權(quán)利要求1所述的橫波估算方法，其特征在于，所述步驟二中，巖性物理模量參數(shù)包括砂巖和泥巖的體積模量、剪切模量和密度，以及孔隙扁度。
4.根據(jù)權(quán)利要求1所述的橫波估算方法，其特征在于，所述步驟二中，流體物理模量參數(shù)包括流體的體積模量、剪切模量和密度。
5.根據(jù)權(quán)利要求1所述的橫波估算方法，其特征在于，所述步驟三中，根據(jù)以下Xu-White方程模型，計(jì)算儲(chǔ)層的孔隙度，建立巖石物理模型:
其中，Vsh、9、a #、a $分別表示儲(chǔ)層的泥質(zhì)含量、孔隙度以及砂和泥的孔隙扁度， VP、VS、P分別表示實(shí)測(cè)的縱波速度、橫波速度和密度， other parameters代表計(jì)算時(shí)所使用的其他的已知的巖性物理模量參數(shù)和流體物理模量參數(shù)。
6.根據(jù)權(quán)利要求1所述的橫波估算方法，其特征在于，所述步驟五中，首先根據(jù)以下Gassmann經(jīng)驗(yàn)方程求取含水巖石的體積模量和剪切模量:
上式中，K7jc為含水巖石體積模量；Kd為干巖石骨架的體積模量；kf7jc為水的體積模量；Km為基質(zhì)的體積模量；Φ為虛擬孔隙度；μ為剪切模量；μ d為干巖石的剪切模量， 然后再根據(jù)以下公式求取含水巖石的縱波速度vp#和橫波速度Vs#:
上式中，P 7jc為含水巖石的密度。
7.根據(jù)權(quán)利要求6所述的橫波估算方法，其特征在于，所述步驟六中，根據(jù)目的層的含水飽和度曲線，在其中含水飽和度Sw〈l的區(qū)域，通過以下Gassmann經(jīng)驗(yàn)方程求取含油氣巖石的體積模量和剪切模量:
上式中，為含油氣巖石體積模量；Kd為干巖石骨架的體積模量；kf_n為油氣的體積模量；Km為基質(zhì)的體積模量；9為虛擬孔隙度；μ為剪切模量；μ d為干巖石的剪切模量，然后再根據(jù)以下公式求取含油氣巖石的縱波速度和橫波速度
上式中，P 為含油氣巖石的密度。
8.—種地球物理勘探橫波估算系統(tǒng)，包括: 第一單元，其根據(jù)測(cè)井曲線分析儲(chǔ)層的巖性特征和流體特征； 第二單元，其根據(jù)儲(chǔ)層的巖性特征和流體特征分別確定相應(yīng)的巖性物理模量參數(shù)和流體物理模量參數(shù)； 第三單元，其建立儲(chǔ)層的巖石物理模型； 第四單元，其根據(jù)巖石物理模型估算儲(chǔ)層的縱波速度，分析估算的縱波速度與實(shí)測(cè)的縱波速度之間的誤差，以最小化此誤差為目標(biāo)，求取相應(yīng)的虛擬孔隙度； 第五單元，其基于虛擬孔隙度求取儲(chǔ)層含水巖石的橫波速度； 第六單元，其通過流體替換求取儲(chǔ)層含油氣巖石的橫波速度； 第七單元，其輸出橫波速度曲線。
【文檔編號(hào)】G01V1/28GK104181585SQ201410302430
【公開日】2014年12月3日 申請(qǐng)日期:2014年6月27日 優(yōu)先權(quán)日:2014年6月27日
【發(fā)明者】孫振濤, 馬永強(qiáng) 申請(qǐng)人:中國石油化工股份有限公司, 中國石油化工股份有限公司石油物探技術(shù)研究院