基于SC-Gassmann模型的縱橫波速度預(yù)測

張祖豪， 林 凱， 賀錫雷

(成都理工大學(xué) 地球物理學(xué)院，成都 610059)

基于SC-Gassmann模型的縱橫波速度預(yù)測

張祖豪， 林 凱， 賀錫雷

(成都理工大學(xué) 地球物理學(xué)院，成都 610059)

通過自相容近似模型和Biot理論相結(jié)合建立巖石物理模型，提出基于SC-Gassmann模型的速度預(yù)測方法。從實(shí)際測井資料中提取巖石物性參數(shù)，利用測井資料基于V-R-H模型來計(jì)算巖石基質(zhì)彈性模量，然后通過自相容近似模型和Biot理論結(jié)合建立的巖石物理模型來計(jì)算干巖樣骨架彈性模量，最后通過Gassmann方程加入流體，求得飽和巖石彈性模量，并預(yù)測儲(chǔ)層巖石縱橫波速度。實(shí)際資料計(jì)算表明，通過巖石物理模型計(jì)算的飽和巖石縱橫波速度與實(shí)際數(shù)據(jù)吻合較好，符合實(shí)際地質(zhì)情況，證明該方法可靠有效。

Gassmann方程； 自相容近似模型； 縱橫波速度預(yù)測

0 引言

儲(chǔ)層含流體雙相介質(zhì)的縱橫波速度在當(dāng)前油氣儲(chǔ)層研究中有重要的作用，通過縱橫波速度和屬性參數(shù)來研究儲(chǔ)層的物性和充填其中的流體影響是重要的技術(shù)方法。

基于巖石物理模型是速度預(yù)測的主要手段。許多學(xué)者通過巖石物理測試，試圖建立巖石物性參數(shù)與縱橫波速度之間的經(jīng)驗(yàn)關(guān)系；或者基于巖石物理模型來進(jìn)行速度預(yù)測。Castagna等[4]給出了含水碎屑硅酸巖中縱波速度與橫波速度的經(jīng)驗(yàn)關(guān)系；Han[5]在經(jīng)過巖石物理測試后，得到縱橫波速度間的線性回歸方程；Xu等[6]建立了針對含泥砂巖的等效介質(zhì)模型，也被稱為Xu-White模型；Xu等[7]把Xu-White模型中的砂巖孔隙和泥巖孔隙替換為砂巖孔、泥巖孔、粒間孔和微裂隙四種類型，成功將該模型引入碳酸鹽巖速度預(yù)測中。國內(nèi)也有一些學(xué)者做了相關(guān)研究，李宏兵等[8]提出了基于非線性模擬退火發(fā)反演巖石等效孔隙縱橫比,進(jìn)行儲(chǔ)層孔隙結(jié)構(gòu)評價(jià)和橫波速度預(yù)測的方法。

筆者將Gassmann理論、Biot理論、自相容近似模型與彈性波基礎(chǔ)理論相結(jié)合，同時(shí)考慮巖石基質(zhì)的彈性性質(zhì)、泥質(zhì)含量、孔隙度和孔隙扁率等影響。基于SC-Gassmann模型的縱橫波速度預(yù)測中，巖石物理模型具有堅(jiān)實(shí)的理論基礎(chǔ)，所涉及的參數(shù)具有明確的物理含義，符合地層的實(shí)際情況，有較高的預(yù)測精度。在測井資料的基礎(chǔ)上，基于自相容近似模型來分析巖石在含流體狀態(tài)下的彈性性質(zhì)，并且根據(jù)速度與體積模量和剪切模量的關(guān)系求出含流體狀態(tài)下的巖石的速度，通過速度來對儲(chǔ)層的流體識別提供依據(jù)。首先以Voigt-Reuss-Hill平均模型來計(jì)算巖石基質(zhì)的彈性模量K0；然后根據(jù)Biot理論計(jì)算出干巖樣的體積模量Kdar；最后通過Gassmann理論來計(jì)算含流體飽和巖石的彈性模量Ksat，進(jìn)而得到含流體飽和巖石速度。

1 方法原理

1.1 等效模量的自相容近似

對于含孔隙彈性介質(zhì)的等效彈性模量的計(jì)算來說一般依賴于：①各個(gè)巖石基質(zhì)礦物的彈性模量；②多各個(gè)巖石基質(zhì)礦物的體積含量；③各個(gè)巖石基質(zhì)礦物相互組合在一起的幾何細(xì)節(jié)。大多數(shù)等效彈性模量的計(jì)算方法，都是在一種巖石礦物基質(zhì)的無限背景介質(zhì)中加入另一種巖石礦物基質(zhì), 再用不同的方法來計(jì)算當(dāng)包含物有一定分布時(shí)的等效彈性模量[1-12]。

自相容近似模型含N相混合物的一般公式：

(1)

(2)

式中：i是指第i種基質(zhì)；xi是體積含量；P和Q是與孔隙扁率有關(guān)的幾何參數(shù)；P和Q的上標(biāo)i指的是此幾何參數(shù)是針對具有自相容等效模量Ksat和μsat的背景介質(zhì)中包含物基質(zhì)i[1-12]。

1 .2 Gassmann方程

在干巖石(不等同于空氣飽和巖石，定義為孔隙壓縮只導(dǎo)致骨架體積形變而不誘發(fā)孔隙壓力的變化)的體積模量和孔隙度之間存在一個(gè)關(guān)系[13-24]：

(3)

式中：Kφ是干燥孔隙空間剛度；φ是孔隙度。

由Bettie功能互等定理，可得飽和巖石體積模量和孔隙度之間的類似關(guān)系[8-19]：

(4)

(5)

式中：Ksat是以有效體積模量為Kf1的孔隙流體所飽和巖石的有效體積模量。

Gassmann方程的轉(zhuǎn)換為：

(6)

(7)

1.3 Biot理論

Biot理論是研究飽和多孔介質(zhì)中固體和流體之間變形關(guān)系的基本理論。在低頻時(shí)，Biot方程變?yōu)镚assmann關(guān)系式，可以說，Biot理論是對Gassmann理論的推廣。由于孔隙流體的粘滯性對不同頻率波的衰減影響不同，因此，Biot理論根據(jù)孔隙流體的相對運(yùn)動(dòng)將孔隙巖石彈性性質(zhì)分為兩個(gè)頻率范圍來描述，并定義了由“低”頻域躍遷到“高”頻域的“臨界頻率”fc為[13-24]：

(8)

式中：η為孔隙流體粘滯系數(shù)；φ為孔隙度；ρf1為孔隙流體密度；k為滲透率。

Geertsma和Smit(1961)根據(jù)巖石骨架性質(zhì)預(yù)測飽和巖石的頻相關(guān)速度可以表示為[8-19]：

(9)

其中：Vp0、Vp∞是零頻率和無窮大頻率的速度值。零頻率飽和流體巖石的縱波速度與Gassmann方程一致[13-24]，即

(10)

無窮大頻率飽和流體巖石的縱波速度的近似表達(dá)式為：

(11)

橫波速度公式：

(12)

其中：β有時(shí)被稱為Biot系數(shù)，β=1-Kdry/K0；ρ0為巖石基質(zhì)礦物顆粒的平均密度；ρsat為飽和巖石的低頻密度，ρsat=(1-φ)ρ0+φρf1；α為曲折度(也稱構(gòu)造因子)，是一個(gè)與固體或流體密度無關(guān)的純幾何參數(shù)，總大于1。而對于低頻的下限速度，Biot模型的零頻率Vp0表達(dá)式與Gassmann預(yù)言的一樣[13-24]。

1.4 基于自相容近似模型的縱橫波速度計(jì)算流程

筆者將Gassmann理論、Biot理論、自相容近似模型與彈性波基礎(chǔ)理論相結(jié)合，同時(shí)考慮巖石基質(zhì)的彈性性質(zhì)、孔隙度和孔隙扁率等影響。提出基于SC-Gassmann模型的縱橫波速度計(jì)算流程(圖1)。

圖1 基于SC-Gassmann模型的縱橫波速度計(jì)算流程Fig.1 The calculation of the waves velocity based on the SC-Gassmann model

該方法包括以下步驟：

1)利用測井資料中的各種礦物的彈性模量和體積含量，采用V-R-H模型計(jì)算巖石基質(zhì)彈性模量。

2)利用自相容近似模型和Biot理論計(jì)算干巖樣的彈性模量。

3)利用Wood方程計(jì)算孔隙流體的彈性模量。

4)向孔隙系統(tǒng)中加入孔隙流體，利用Gassmann方程計(jì)算飽和流體巖石的彈性模量。

5)利用巖石彈性性質(zhì)和速度之間的關(guān)系計(jì)算巖石縱橫波速度。

2 實(shí)際資料計(jì)算與分析

利用某工區(qū)的測井?dāng)?shù)據(jù)，進(jìn)行了縱橫波速度預(yù)測，該地區(qū)巖石特性為低孔低滲高含水飽和度的砂泥巖。圖2中的VP、VS是實(shí)際測井資料中的巖石速度曲線，VP-FL、VS-FL是基于自相容近似模型的速度預(yù)測方法計(jì)算得到的巖石速度。

由圖2可以看出，速度預(yù)測結(jié)果與實(shí)際測井?dāng)?shù)據(jù)的巖石速度基本吻合。在深度為2 418 m～2 420 m這一段地層的泥質(zhì)含量、孔隙度降低，密度升高，地層的巖性發(fā)生了變化，通過修正巖石的物性參數(shù)，最后得到計(jì)算結(jié)果與測井速度曲線基本吻合。由圖3可以看出，速度預(yù)測結(jié)果與實(shí)際測井?dāng)?shù)據(jù)的巖石速度比較吻合，且其中有一套薄泥層，速度預(yù)測計(jì)算結(jié)果與測井速度曲線吻合，符合實(shí)際地質(zhì)情況。由圖4可以看出，速度預(yù)測結(jié)果與實(shí)際測井?dāng)?shù)據(jù)的巖石速度比較吻合，并且其中有一套薄泥層，速度預(yù)測計(jì)算結(jié)果與測井速度曲線吻合，符合實(shí)際地質(zhì)情況。

圖2 某9井速度預(yù)測計(jì)算結(jié)果與測井速度曲線Fig.2 Calculation results of X9 well velocity prediction and well log velocity curve

圖3 某107井速度預(yù)測計(jì)算結(jié)果與測井速度曲線Fig.3 Calculation results of X107 well velocity prediction and well log velocity curve

圖4 某111井速度預(yù)測計(jì)算結(jié)果與測井速度曲線Fig.4 Calculation results of X111 well velocity prediction and well log velocity curve

3 結(jié)論

筆者基于實(shí)際的地層資料，利用SC-Gassmann模型來進(jìn)行速度預(yù)測，在某工區(qū)計(jì)算得到的速度預(yù)測結(jié)果與實(shí)際的測井速度曲線吻合較好，地層中的泥巖薄層能夠有效的適用。而且在地層巖性變化時(shí)，修正巖性的物性參數(shù)，可以得到與實(shí)際測井曲線相吻合的計(jì)算結(jié)果。結(jié)果表明，本方法在砂泥巖地層應(yīng)用效果良好，對砂泥巖薄互層有敏感的反映。因此筆者提出的方法是可靠有效的，是符合實(shí)際地質(zhì)情況的。

[1] MAVKO,G.,MUKERJI,T.,DVORKIN,J. The rock physics handbook-tools for seismic in porous media[M]. Cambridge: Cambridge University Press,2003.

[2] LIU, W.L. Reservoir exploitation seismic techniques[M].Beijing: Petroleum Industrial Press,1996.

[3] CHEN, R. HUANG, T.F. Petrophysics[M].Beijing:Peking University Press, 2001.

[4] CASTAGNA J P, BATZLE M L, EASTWOOD R L. Relationships between compressional-wave and shear-wave velocities in clastic silicate rocks[J]. Geophysics, 1985,50(4):571-581.

[5] HAN D H, NUR A, MORGAN D. Effects of porosity and clay content on wave velocities in sandstones[J]. Geophysics, 1986,51(11):2093-2107.

[6] XU S, WHITE R E. A new velocity model for clay-sand mixtures[J]. Geophysical Prospecting,1995,43(1):91-118.

[7] XU S, PAYNE M A. Modeling elastic properties in carbonate rocks[J]. The Leading Edge,2009,28(1):66-74.

[8] 李宏兵,張佳佳,姚逢昌.巖石的等效孔隙縱橫比反演及其應(yīng)用[J].地球物理學(xué)報(bào),2013,56(2):608-615. LI H B, ZHANG J J, YAO F C. Inversion of effective pore aspect ratios for porous rocks and its applications[J]. Chinese Journal of Geophysics,2013,56(2):608-615.(In Chinese)

[9] 董寧,霍志周,孫贊東,等.泥頁巖巖石物理建模研究[J].地球物理學(xué)報(bào),2014,57(6):1990-1998. DONG N, HUO Z, SUN Z D, et al. An investigation of a new rock physics model for shale[J].Geophys,2014,57(6):1990-1998.(In Chinese)

[10]胡起,陳小宏,李景葉.基于各向異性巖石物理模型的頁巖氣儲(chǔ)層橫波速度預(yù)測[J].石油物探,2014,53(3):254-261. HU Q, CHEN X H, LI J Y. Shear wave velocity prediction for shale gas reservoirs based on anisotropic rock physics model[J].Geophysical Prospecting for Petroleum,2014,53(3):254-261. (In Chinese)

[11]林凱,賀振華,熊曉軍,等. 基于Gassmann方程的流體替換流程[J].長江大學(xué)學(xué)報(bào),2009,(6)1:180-181. LIN K, HE Z H, XIONG X J, et al. Based on the Gassmann equation fluid substitution flow[J].Journal of Yangtze University(Nat Sci Edit),2009,(6)1:180-181. (In Chinese)

[12]史燕紅. 基于Gassmann方程的流體替換[D].成都:成都理工大學(xué),2009. SHI Y H. The fluid substitution based on Gassmann equation[D].Chengdu:Chengdu University of Technology,2009. (In Chinese)

[13]王炳章,朱曄,王丹.多孔介質(zhì)的流體機(jī)制模型及其頻散機(jī)理 [J].勘探地球物理進(jìn)展 ,2008,31(6):405-413. WANG B Z, ZHU Y, WANG D. Fluid mechanism models and their velocity dispersions in porous media [J].Progress in Exploration Geophysics ,2008,31(6):405-413. (In Chinese)

[14]王炳章.地震巖石物理學(xué)及其應(yīng)用研究[D].成都:成都理工大學(xué),2008. WANG B Z. Seismic rock physics and its applied research[D]. Chengdu:Chengdu University of Technology,2008. (In Chinese)

[15]陳潔.蘇里格氣田儲(chǔ)層的AVO分析及流體識別[D].荊州：長江大學(xué),2012. CHEN J. AVO analysis and fluid identification in Su lige gas field reservoir[D]. Jingzhou:Yangtze University,2012. (In Chinese)

[16]薛濤.時(shí)移地震技術(shù)在油田開發(fā)區(qū)塊中的應(yīng)用研究[D].西安：長安大學(xué),2008. XUE T. The application of research time-lapse seismic technology in the development block of oil field[D]. Xi'an:Chang’an University,2008. (In Chinese)

[17]江偉.多參數(shù)流體識別方法及應(yīng)用研究 [D].成都:成都理工大學(xué),2010. JIANG W. Multi-parameter fluid identification and applied resarch [D]. Chengdu:Chengdu University of Technology,2010. (In Chinese)

[18]李國福.多參數(shù)儲(chǔ)層預(yù)測及流體識別方法研究[D].成都:成都理工大學(xué),2011. LI G F. Research on multi-parameters reservoir prediction and fluid identification methods[D]. Chengdu:Chengdu University of Technology,2011. (In Chinese)

[19]撒利明.儲(chǔ)層反演油氣檢測理論方法研究及其應(yīng)用[D].北京：中國科學(xué)院研究生院,2003. SA L M. Research and application of reservoir inversion theory and method of oil and gas detection[D]. Beijing:Graduate School of Chinese Academy of Sciences,2003. (In Chinese)

[20]鎮(zhèn)晶晶,劉洋.裂縫介質(zhì)巖石物理模型研究綜述[J].地球物理學(xué)進(jìn)展,2011,26(5):1708-1716. ZHEN J J, LIU Y. Review over physical model of fractured rock medium[J].Progress in Geophys,2011,26(5):1708-1716. (In Chinese)

[21]鄒文.基于地震資料的流體識別技術(shù)研究[D].成都:成都理工大學(xué),2008. ZOU W. Fluid identification technique research based on seismic data[D].Chengdu:Chengdu University of Technology,2008. (In Chinese)

[22]林凱.基于Gassmann方程的彈性參數(shù)數(shù)值模擬研究 [D].成都:成都理工大學(xué),2010. LIN K. Based on Gassmann equation Elastic parameters Numerical simulation[D]. Chengdu:Chengdu University of Technology,2010. (In Chinese)

[23]鄔龍.碳酸鹽巖儲(chǔ)層流體替代及其AVO分析[D].上海：同濟(jì)大學(xué),2014. WU L. Fluid substitution and its AVO analysis in carbonate reservoir[D]. Shanghai:Tongji University,2014. (In Chinese)

[24]馬妮.多參數(shù)儲(chǔ)層流體識別方法研究[D].成都:成都理工大學(xué),2014. MA N. Multi-parameters reservoir fluid identification methods[D].Chengdu:Chengdu University of Technology,2014. (In Chinese)

Waves velocity prediction based on self-consistent model

ZHANG Zuhao, LIN Kai, HE Xilei

(College of Geophysics, Chengdu University of Technology，Chengdu 610059, China)

In this paper, the Gassmann theory, Biot theory and elastic wave theory are combined to quantitatively analyze the relationship between the physical parameters of the rock and the seismic dynamic parameters. The physical parameters of rocks are extracted from the actual well logging data. Based on the self-consistent model, the log data is used to predict the wave velocity. The elastic properties of rock with fluid and the waves velocity are calculated through the rock physical model, which provides a basis for reservoir prediction and fluid identification in this paper. The predicted results obtained by this method are in agreement with the log data. This is consistent with the actual geological conditions, which proves that the method is reliable and effective.

Gassmann equation; self-consistent model; waves velocity prediction

2016-11-18 改回日期：2017-03-16

國家自然科學(xué)基金項(xiàng)目(41374111);國家自然科學(xué)基金青年基金項(xiàng)目(41404101)

張祖豪(1992-)，男，碩士，研究方向?yàn)橛蜌獾厍蛭锢恚?E-mail：296023193@qq.com。

1001-1749(2017)02-0248-05

P 631.4

A

10.3969/j.issn.1001-1749.2017.02.15