基于广义弹性阻抗的流体识别因子反演方法研究与应用

贾凌云,李　琳,王千遥,马劲风,王浩璠,王大兴

(1.西北大学地质学系,二氧化碳捕集与封存技术国家地方联合工程研究中心,陕西西安710069;2.中煤科工集团西安研究院有限公司,陕西西安710077;3.中国石油长庆油田公司勘探开发研究院,陕西西安710018)

对于常规储层,采用一般的流体识别因子可以达到识别流体的目的;而对于致密储层,由于其孔隙结构复杂、非均质强,采用常规的流体识别因子会产生一些流体识别假象[1]。国内外学者基于弹性参数提出了不同的流体识别因子,且其求取的方法也不同。GOODWAY等[2]提出了拉梅参数反演(Lambda-Mu-Rho,LMR)流体识别技术,即在纵、横波阻抗的平方运算和差运算后提取出λρ和μρ来识别储层中的流体。RUSSELL[3]使用BIOT[4]和GASSMANN[5]孔隙介质弹性理论,引入了岩石的体变模量K和切变模量μ,化简了GOODWAY等[2]的公式,进一步估算出了储层的流体识别因子。QUAKENBUSH等[6]提出了泊松阻抗的概念,提高了对岩性及流体的识别能力。目前致密储层识别流体的主要方法是通过叠前弹性参数反演获得纵横波速度、密度、纵横波阻抗等数据体,间接推算出Russell流体因子[7-8]或其它流体敏感参数进行识别,获得有利孔隙砂岩的分布[9]。地震属性间接计算流体识别敏感因子,会产生较大的累积误差。许多学者提出流体敏感因子的直接提取方法,避免了间接计算带来的累积误差。王保丽等[10]用拉梅系数λ和剪切模量μ表示弹性阻抗公式,认为直接求取比通过纵横波速度和密度间接求取精度高,对流体的反映较为清楚。宗兆云等[11]基于叠前弹性阻抗反演的方法,实现了拉梅参数的直接提取。桂金咏等[12]提出了弹性阻抗直接反演泊松比的方法。孙瑞莹等[13]基于随机地震反演,提出了Russell流体因子的直接估算方法。印兴耀等[14-15]认为Gassmann流体项识别流体较Russell流体因子或其它弹性参数敏感,推导出了含Gassmann流体项的弹性阻抗公式并研究了流体项的直接提取方法。杨培杰等[16]建立了一种敏感因子定量分析方法,认为Gassmann流体项对研究区流体识别最为敏感,并通过反演方法直接将其提取出来。

Zoeppritz方程[17]表示了P波入射到弹性介质界面时P波和S波的反射和透射系数的数学方程,由于Zoeppritz方程的非线性形式使得振幅估算的反演结果不稳定,因此地震AVO分析或叠前反演常用Zoeppritz方程的线性近似公式[18-22],并且这些近似公式直接引入了容易理解的地震参数。其中WANG[20]用流体项和刚性项表达了Zoeppritz方程,马劲风等人对其进一步化解[23-25],提出了广义弹性阻抗(GEI)和射线弹性阻抗(REI)计算方法,该方法对Zoeppritz方程的近似具有较高的精度,并可以分解为流体项和刚性项。

埋深较大或致密储层受强烈的压实和后期成岩作用的改造,孔隙结构复杂[26-27],与常规储层地质特征差异大,地震资料存在信噪比低、高频信号弱等问题。常规叠前弹性阻抗反演方法在识别流体中虽然被广泛应用[28-30],但应用于致密砂岩储层流体识别时效果不理想。本文将广义弹性阻抗方程中的流体项进一步化解为含Gassmann方程流体项f的公式,可以直接反演获得参数f/vS。该参数是一种新的流体识别因子,避免了参数间接计算的累积误差,且该流体因子对流体的识别较为敏感,提高了致密储层的流体识别能力。将该方法应用于苏里格致密砂岩储层的流体识别与预测中,取得了良好效果。

1　理论公式

1.1　流体因子敏感指示系数

DILLON等[7]提出了用流体因子敏感指示系数(S)来分析流体敏感性的方法,其值等于含烃流体因子与含水流体因子平均值之差的绝对值与含烃流体因子标准差的比值,针对气层可表示为:

(1)

1.2　基于GEI的敏感流体项弹性阻抗方程

弹性阻抗在油气储层预测中得到了较好的应用,但随着地震勘探难度的增加,弹性阻抗已经不能完全满足油气预测的要求。马劲风等人对弹性阻抗进行深入研究后发现弹性阻抗本身存在一些缺点和不足[23-25],即弹性阻抗随入射角的增大,误差逐渐增大,容易导致岩性识别的错误,针对此问题,将上下岩层密度的变化率考虑为横波速度变化率的线性关系,并推导出了含有流体项和刚性项的广义弹性阻抗值IGE公式,具体为:

(2)

基于Gassmann理论,多孔流体饱和岩石的纵横波速度可以表示为:

式中:Kdry为干岩石体变模量;μ为岩石的切变模量,μ=μdry=μsat;f为流体、孔隙项。

(6)

进一步可以表示为:

(7)

将公式(7)代入公式(2)中,化简后,有:

(8)

(9)

1.3　精度分析

为验证广义弹性阻抗方程的合理性和精确度,以模型数据为基础分析公式(2)得到反射系数的近似精度,依据苏里格气田苏46井测井数据建立两层介质模型,盖层为泥岩层,储层为砂岩层,具体的弹性参数取值如表1所示。

假设最大入射角为30°,分别采用精确Zoeppritz方程、常规弹性阻抗(EI)方程和广义弹性阻抗(GEI)方程计算两层模型的反射系数,利用AVO曲线定量对比GEI方程和EI方程与精确Zoeppritz方程之间的误差(图1,图2)。入射角较小时,EI方程和GEI

表1　两层模型参数

图1　反射系数计算结果

图2　误差分析

方程计算的反射系数与精确Zoeppritz方程很接近,随入射角度增加,二者逐渐偏离精确解,GEI方程的精度高于EI方程的精度,最大误差为-0.001左右。因此,广义弹性阻抗方程应用于致密砂岩储层可以提高储层预测精度。

2　参数求取

干岩石纵横波速度比γdry的计算方法[31-34]主要分为以下两种:一种是实验室测量方法,一种是基于实际测井资料并结合Gassmann方程的估算方法。实际应用中需根据研究区的具体条件选择合理的求取方法,只有得到适应工区的γdry,才能保证反演流体因子的精度。本文采用Gassmann方程求出测井资料目的层段所有采样点的γdry,求平均值后作为研究区的干岩石纵横波速度比γdry。

反演时将饱和岩石的纵横波速度比γsat作为反演结果的一个参数,即可以通过公式(8)进行广义弹性阻抗反演求得饱和岩石纵横波速度比γsat。

2.2　岩性和流体参数求取

当入射角度不变时,本文推导的公式(9)中每个采样点的广义弹性阻抗值IGE是未知量γsat和f/vS的函数,选用两个角度的广义弹性阻抗值IGE(θ),可得方程组

(10)

解方程(10)可求出参数γsat和f/vS。进一步可由饱和岩石的纵横波速度比γsat通过公式(11)计算出泊松比,用来识别岩性或流体;横波速度vS对剪切模量变化不敏感,其随流体变化或含水饱和度变化几乎不变,而流体项f对流体较敏感,因此f/vS的流体敏感指示系数与f相近,可用来识别流体变化。在苏里格致密砂岩储层中,广义弹性阻抗方法较常规弹性阻抗方法更精确,利用广义弹性阻抗反演方法求取的流体识别因子具有较高的流体识别能力。

(11)

3　实际应用及效果分析

鄂尔多斯盆地苏里格气田上古生界二叠系石盒子组盒8段属于特低渗致密砂岩储层,分为盒8上和盒8下两段,储层孔隙类型多样、孔隙结构复杂、非均质性强,因而气、水分布复杂。含气砂岩和围岩的弹性参数存在较大范围的叠置,储层预测存在明显的多解性,很难预测致密砂岩储层中的流体有利区,需要寻找致密砂岩储层有效的岩性识别参数和流体识别因子。

流体检测的方法一般是采用弹性数据直接求出流体参数来预测流体[35],本文采用的方法是:先识别岩性,再在目的层中求出流体参数,最后进行流体预测。

3.1　岩性识别

致密砂岩储层的围岩对流体识别影响较大[36]。首先分析研究区内岩性识别的敏感参数,将砂岩储层从岩层中识别出来。根据研究区盒8段测井资料,选用纵、横波阻抗,纵横波速度比,泊松比和广义弹性阻抗的数据运算与自然伽马数据进行交会分析(图3),对比可知,纵、横波阻抗难以区分砂、泥岩,纵横波速度比虽有一定的区分能力但没有泊松比清晰,广义弹性阻抗的简单数据运算也可区分岩性,但在实际应用中不可避免地会产生累积误差。因此,选择泊松比作为最优岩性识别参数来识别致密砂岩储层最合理。

3.2　流体敏感因子分析

不同的流体因子在不同地质条件下表现出不同的敏感特征。本文选用几项流体因子对致密砂岩储层进行流体识别,采用公式(1)计算各项流体因子的流体指示系数(表2),并绘制成直方图(图4),可直观看到致密砂岩储层中常规属性如纵波速度、纵横波速度比、泊松比等流体指示系数均较小且差别不大,难以识别流体,而Gassmann流体项f及f/vS对流体的指示较为敏感,其中f/vS为最敏感的流体识别因子。

Gassmann流体项f主要是含水饱和度和孔隙度的函数,苏里格气田盒8段致密砂岩储层的孔隙度范围为2.2%～14.1%,平均值为7.4%,即孔隙度变化不大,对流体项的影响较小,因此流体项f或f/vS主要受含水饱和度影响,是流体识别的敏感因子。致密砂岩储层干岩石的各项弹性参数与流体项相比大得多[37],包含干岩石信息的常规参数识别流体的能力较低,流体识别指示系数差异不大,而流体项f或f/vS去除了干岩石的影响,对流体的识别效果尤为明显。Gassmann方程中假设流体替换前后剪切模量不变,地层含气较含水密度小,根据剪切模量与密度和横波速度的关系((3b)式)可知,含气横波速度较含水横波速度有轻微的增大,流体识别因子f/vS的流体指示系数有可能比f略大。

图3　岩性识别敏感参数分析a 纵波阻抗与伽马交会结果; b 横波阻抗与伽马交会结果; c 纵横波速度比与伽马交会结果; d 泊松比与伽马交会结果; e IGE0-IGE30°与伽马交会结果; f(IGE0-IGE30°)/IGE0与伽马交会结果

IP/(g·cm-3·km·s-1)vP/vSσ(k-μ)/GPaμρ/(GPa·g·cm-3)f/GPafvS/(GPa·km-1·s)(I2P-γ2dryI2S)/(GPa·g·cm-3)含气砂岩平均值11.2841.5800.1623.42251.1960.1330.0450.062标准差0.8550.0600.2103.8606.6790.1610.0509.651含水砂岩平均值11.7341.6260.1946.31052.3163.7421.3017.369流体指示系数0.5260.7680.8160.7480.16822.42925.2190.757

图4　不同流体因子的流体指示系数

图5对比了完全含水和完全含气情况下流体因子f/vS和泊松比σ与孔隙度的交会结果,可以看出,f/vS参数对流体指示的效果较为明显,是流体指示的最佳敏感参数。

3.3　广义弹性阻抗正反演

3.3.1正演分析

采用主频为30Hz的Ricker子波,1ms的时间采样率,分别由精确Zoeppritz方程、EI方程和GEI方程对致密砂岩储层测井数据进行正演模拟,并求出EI和GEI的正演误差(图6)。对比可知,随入射角增大,GEI与EI计算的误差逐渐增大,但GEI误差相对较小,因此GEI比EI更适合本文研究区实际地质情况。盒8段致密砂岩储层深度范围为3562.8～3621.3m(表3),时间范围为1420～1446ms(图6)。

图5　流体因子敏感性分析a f/vS与孔隙度交会结果; b 泊松比与孔隙度交会结果

表3　测井储层段解释信息

3.3.2反演效果分析

对精确Zoeppritz方程得到的正演模型进行GEI反演,同样采用主频为30Hz的Ricker子波,获得不同入射角度的广义弹性阻抗GEI数据。图7为实际地层在入射角0,15°,25°的GEI反演结果与测井曲线计算相对应入射角的GEI值对比结果,可以看出,反演的GEI值与测井计算的GEI值相似度较高,说明本方法对参数GEI值的估算合理。选用入射角为0和15°的GEI反演数据代入公式(10)计算出纵横波速度比、泊松比及流体因子f/vS数据(图8)。不同流体因子之间的量纲虽有所不同,但可以看到f/vS对流体的反映最为清晰。

由前文可知,在研究区中,泊松比为最优岩性识别参数,f/vS为最优流体识别因子。因此,本文首先采用反演的泊松比识别砂岩,识别出储层段,然后来圈定反演的f/vS流体因子。与测井计算的储层f/vS流体因子对比可知,两者相似性较高,结果如图9所示。图9中,第一列为反演获得的f/vS流体因子,第二列为测井计算的f/vS流体因子,加入10%的噪声后反演流体因子产生了一定程度的误差,但仍然可以有效识别出流体。同时,图9还显示,孔隙度越大,f/vS随含水饱和度的变化越明显,孔隙度在0.05以下时,f/vS随含水饱和度的变化难以识别。可见,GEI直接反演致密砂岩储层流体因子的方法具有一定的有效性和较好的抗噪性。

图6　正演模拟

图7　GEI反演

图8　GEI反演参数

3.3.3实际资料应用效果

将广义弹性阻抗直接反演纵横波速度比(或泊松比)和流体识别因子f/vS的方法应用于苏里格致密砂岩储层预测中,获得泊松比剖面和流体识别因子f/vS剖面。图10为泊松比反演剖面,图中红色层位为盒8下段的底界面,储层显示为低泊松比值。剖面中的A井、B井和C井在盒8段测井解释的砂岩厚度分别为32.5,16.0,15.5m;气层总厚度分别为6.6,2.2,5.1m;含气饱和度分别为39.3%,56.0%,63.0%。从泊松比反演剖面中识别出砂岩储层,圈定流体识别因子f/vS的范围,并估算出储层的含气饱和度剖面,结果如图11所示,红色为高含气饱和度,蓝色为低含气饱和度,与测井解释结果一致。利用广义弹性阻抗反演获得的f/vS参数可以有效预测出致密砂岩储层的含气饱和度,有利于目标区和目标层的识别。

图9　GEI反演流体因子与测井计算值对比

图10　泊松比反演剖面

图11　流体饱和度剖面

4　结论

由于致密砂岩储层孔隙度小,渗透率低,导致其与围岩的密度、速度等弹性参数差异小,常规弹性阻抗反演对致密砂岩储层流体识别难度大,本文采用广义弹性阻抗反演的方法,直接求取流体识别因子f/vS,对苏里格致密砂岩储层进行流体识别,取得了较好的效果,得到如下结论:

1) 在苏里格致密砂岩气藏中,广义弹性阻抗方法较弹性阻抗方法精度更高,能更好地为储层预测提供技术服务;

2) 将广义弹性阻抗公式表示为纵横波速度比vP/vS(可计算泊松比σ)和f/vS的函数,直接反演出岩性识别参数σ和流体识别因子f/vS,避免了间接计算的累积误差,提高了各个参数的估算精度。泊松比σ是苏里格致密砂岩岩性识别的最优参数,f/vS流体识别因子能够较好地满足致密砂岩储层中流体的识别需要,并具有一定的抗噪性;

3) 致密砂岩储层的围岩对流体的识别影响较大,在实际应用中需要首先识别出致密砂岩储层,再针对储层采用流体敏感因子识别出流体,才能取得理想的流体识别效果。

[1]张世鑫.基于地震信息的流体识别方法研究与应用[D].东营:中国石油大学,2012

ZHANG S X.Methodology and application of fluid identification with seismic information[D].Dongying:China University of Petroleum,2012

[2]GOODWAY W N,CHEN T,DOWNTON J.Improved AVO fluid detection and lithology discrimination using Lamé petrophysical parameters; “kk”,“ll”,and “k/μ fluid stack,” from P and S inversions[J].Expanded Abstracts of 67thAnnual Internat SEG Mtg,1997:183-186

[3]RUSSELL B,HEDLIN K,HILTERMAN F,et al.Fluid-property discrimination with AVO:a Biot-Gassmann perspective[J].Geophysics,2003,68(1):29-39

[4]BIOT M A.General theory of three-dimensional consolidation[J].Journal of Applied Physics,1941,12(2):155-164

[5]GASSMANN F.Elastic waves through a packing of spheres[J].Geophysics,1951,16(3):673-682

[6]QUAKENBUSH M,SHANG B,TUTTLE C.Poisson impedance[J].The Leading Edge,2006,25(2):128-138

[7]DILLON L,SCHWEDERSKY G,VSQUEZ G,et al.A multiscale DHI elastic attributes evaluation[J].The Leading Edge,2003,22(10):1024-1029

[8]姜仁,欧阳永林,曾庆才,等.Russell流体因子在致密砂岩气层检测中的应用[J].天然气工业,2017,37(1):76-81

JIANG R,OUYANG Y L,ZENG Q C,et al.Application of the Russell fluid factor in tight sandstone gas detection[J].Natural Gas Industry,2017,37(1):76-81

[9]张建坤,吴吉忠,徐文会,等.基于岩石物理的致密砂岩油藏叠前优势储层预测——以渤海湾盆地南堡凹陷4号构造东营组二段为例[J].石油学报,2017,38(7):793-803

ZHANG J K,WU J Z,XU W H,et al.Pre-stack favorable reservoir prediction in tight sandstone reservoirs based on rock physics:a case study of the 2ndMember of Dongying Formation,Structure 4,Nanpu sag,Bohai Bay Basin[J].Acta Petrolei Sinica,2017,38(7):793-803

[10]王保丽,印兴耀,张繁昌.基于Gray 近似的弹性波阻抗方程及反演[J].石油地球物理勘探,2007,42(4):435-439

WANG B L,YIN X Y,ZHANG F C.Gray approximation based elastic wave impedance equation and inversion[J].Oil Geophysical Prospecting,2007,42(4):435-439

[11]宗兆云,印兴耀,张繁昌.基于弹性阻抗贝叶斯反演的拉梅参数提取方法研究[J].石油地球物理勘探,2011,46(4):598-604

ZONG Z Y,YIN X Y,ZHANG F C.Elastic impedance Bayesian inversion for lame parameters extracting[J].Oil Geophysical Prospecting,2011,46(4):598-604

[12]桂金咏,印兴耀,曹丹平.基于弹性阻抗反演理论的泊松比反演方法研究[J].石油物探,2011,50(5):463-469

GUI J Y,YIN X Y,CAO D P.Poisson’s ratio inversion method based on elastic impedance inversion theory[J].Geophysical Prospecting for Petroleum,2011,50(5):463-469

[13]孙瑞莹,印兴耀,王保丽,等.基于随机地震反演的Russell流体因子直接估算方法[J].地球物理学报,2016,59(3):1143-1150

SUN R Y,YIN X Y,WANG B L,et al.A direct estimation method for the Russell fluid factor based on stochastic seismic inversion[J].Chinese Journal of Geophysics,2016,59(3):1143-1150

[14]印兴耀,张世鑫,张繁昌,等.利用基于Russell近似的弹性波阻抗反演进行储层描述和流体识别[J].石油地球物理勘探,2010,45(3):373-380

YIN X Y,ZHANG S X,ZHANG F C,et al.Utilizing Russell approximation-based elastic wave impedance inversion to conduct reservoir description and fluid identification[J].Oil Geophysical Prospecting,2010,45(3):373-380

[15]印兴耀,张世鑫,张峰.针对深层流体识别的两项弹性阻抗反演与Russell流体因子直接估算方法研究[J].地球物理学报,2013,56(7):2378-2390

YIN X Y,ZHANG S X,ZHANG F.Two-term elastic impedance inversion and Russell fluid factor direct estimation method for deep reservoir fluid identification[J].Chinese Journal of Geophysics,2013,56(7):2378-2390

[16]杨培杰,董兆丽,刘昌毅,等.敏感流体因子定量分析与直接提取[J].石油地球物理勘探,2016,51(1):158-164

YANG P J,DONG Z L,LIU C Y,et al.Sensitive fluid factor extraction and analysis[J].Oil Geophysical Prospecting,2016,51(1):158-164

[17]ZOEPPRITZ K,ERDBEBENWELLEN Ü.On the reflection and propagation of seismic waves[J].Göttinger Nachrichten,1919(1):66-84

[18]BORTFELD R.Approximations to the reflection and transmission coefficients of plane longitudinal and transverse waves[J].Geophysical Prospecting,1961,9(4):485-502

[19]AKI K,RICHARDS P G.Quantitative seismology:theory and methods[M].San Francisco:Freeman W H and Company,1980:1-557

[20]WANG Y H.Approximations to the Zoeppritz equations and their use in AVO analysis[J].Geophysics,1999,64(6):1920-1927

[21]CONNOLLY P.Elastic impedance[J].The Leading Edge,1999,18(4):438-452

[22]GRAY F D.The application of AVO and inversion to the estimation of rock properties[J].Expanded Abstracts of 70thAnnual Internat SEG Mtg,2000:549-552

[23]马劲风.地震勘探中广义弹性阻抗的正反演[J].地球物理学报,2003,46(1):118-124

MA J F.Forward modeling and inversion method of generalized elastic impedance in seismic exploration[J].Chinese Journal of Geophysics,2003,46(1):118-124

[24]MA J F.Generalized elastic impedance[J].Expanded Abstracts of 73rdAnnual Internat SEG Mtg,2003:254-257

[25]MA J F,Morozov I B.Ray-path elastic impedance:2004 Canadian Society of Exploration Geophysicists National Convention,AVO[R].Calgary,Canada:CSEG,2004

[26]赵海波,唐晓花,李奎周,等.基于地震岩石物理分析与叠前地质统计学反演技术的齐家地区致密薄储层预测[J].石油物探,2017,56(6):853-862

ZHAO H B,TANG X H,LI K Z,et al.Tight thin-bed reservoir prediction using rock physics analysis and prestack geostatistical inversion in the Qijia area[J].Geophysical Prospecting for Petroleum,2017,56(6):853-862

[27]李瑞磊,杨勤林,田建华,等.松辽盆地龙凤山气田致密砂岩含气性预测研究[J].石油物探,2017,56(6):874-881

LI R L,YANG Q L,TIAN J H,et al.Tight sandstone gas prediction in the Longfeng Mountain gas field of Songliao Basin,China[J].Geophysical Prospecting for Petroleum,2017,56(6):874-881

[28]郎晓玲,彭仕宓,康洪全,等.叠前同时反演方法在流体识别中的应用[J].石油物探,2010,49(2):164-169

LANG X L,PENG S M,KANG H Q,et al.Application of prestack simultaneous inversion in fluid identification[J].Geophysical Prospecting for Petroleum,2010,49(2):164-169

XI X Y,LING Y,SUN D S,et al.Application of prestack seismic inversion[J].Geophysical Prospecting for Petroleum,2010,49(4):384-389

[30]龚幸林,牛刚,白文成,等.叠前弹性反演在Y研究区盐下碳酸盐岩气藏预测中的应用[J].石油物探,2011,50(1):65-68

GONG X L,NIU G,BAI W C,et al.Application of pre-stack elastic inversion for sub-salt gas prediction in Block Y[J].Geophysical Prospecting for Petroleum,2011,50(1):65-68

[31]PICKETT G R.Acoustic character log and their application in formation evaluation[J].Journal of Petroleum Technology,1963,15(6):659-667

[32]MURPHY W F,SCHWARTZ L M,HORNBY B.Interpretation physics ofvPandvSin sedimentary rocks[J].Transactions SPWLA 32ndAnnual Logging Symposium,1991:1-24

[33]BRIE A,PAMPURI F,MARSALA A F,et al.Shear sonic interpretation in gas-bearing sands[J].SPE Annual Technical Conference and Exhibition,1995:SPE-30595-MS

[34]LI L,MA J F.Shear wave velocity prediction during CO2-EOR and sequestration in Gao89 area of Shengli Oilfield[J].Applied Geophysics,2017,14(3):372-380

[35]桂金咏,高建虎,李胜军,等.面向实际储层的流体因子优选方法[J].石油地球物理勘探,2015,50(1):129-135

GUI J Y,GAO J H,LI S J,et al.Reservoir-oriented fluid factor optimization method[J].Oil Geophysical Prospecting,2015,50(1):129-135

[36]GULIYEV E,DAVIS T L.Interpretation ofvP/vSvelocity ratio for improved tight gas sandstone reservoir characterization,Rulison field,Colorado[J].Expanded Abstracts of 77thAnnual Internat SEG Mtg,2007 Annual Meeting,2007:1451-1455

[37]MAVKO G,MUKERJI T,DVORKIN J P.The rock physics handbook[M].Cambridge:Cambridge University Press,1998:1-524