地球物理学报  2014, Vol. 57 Issue (6): 1977-1989   PDF    
不同孔隙流体下延长组砂岩波速各向异性的实验研究
乔二伟1, 龙长兴1, 马洪涛2    
1. 国土资源部新构造运动与地质灾害重点实验室, 中国地质科学院地质力学研究所, 北京 100081;
2. 中国石油化工集团公司西北油田分公司勘探开发研究院, 乌鲁木齐 830011
摘要:利用Autolab 2000岩石物性测试设备,在0~200 MPa及三种不同孔隙流体条件下(干燥、水饱和及油饱和)研究了三种来自延长油田延长组的砂岩沿平行和垂直层面两个方向(即H和V)的纵波、横波速度,并分别计算了这些弹性波速的各向异性系数.结果表明:三种砂岩Y1、Y2和Y3两个方向上的VPVSHVSV均随围压增加而基本呈对数函数形式变化,且纵波饱油曲线的斜率远大于饱水;Y1、Y2砂岩H方向的VPVSHVSV在干燥、饱水和饱油条件下基本上均大于对应的V方向的波速;干燥、饱水和饱油条件下Y1、Y2的波速各向异性系数εγζ和围压之间呈指数或二次函数降低;Y1、Y2、Y3三个岩样之间的速度各向异性差异明显,以Y2的各向异性最显著,其次是Y1,而Y3基本没有各向异性.另外,Y1、Y2的各向异性系数εγζ之间表现出明显的与压力和流体状态相关性,且在实验压力范围内,饱油条件下的εγ均大于对应的饱水条件.实验结果可以为该地区地震资料的解释、油水层划分及与声波测井之间的对比提供重要的基础数据.
关键词干燥     水饱和     油饱和     延长组     砂岩     各向异性    
An experimental study of P- and S-wave velocities anisotropy in sandstones filled with differential pore fluids from the Yanchang Formation
QIAO Er-Wei1, LONG Chang-Xing1, MA Hong-Tao2    
1. Key Laboratory of Neotectonic Movement and Geohazard, Ministry of Land and Resources, Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, China;
2. Exploration and Development Research Institute, Northwest Oilfield Company, SINOPEC, Vrümqi 830011, China
Abstract: This work applied AutoLab 2000 computer-controlled servo hydraulic triaxial test system to measure accurate P- and S-wave velocities at pressures up to 180 MPa on three kinds of core samples of sandstones that were cut parallel and perpendicular to the bedding planes (i.e., H and V) in dry, water and oil-saturated conditions from the Yanchang Formation in the Yanchang oilfield, and the wave velocity anisotropy parameters were also calculated. The results show the wave velocities VP. VSH, and VSV for all H and V samples cut from three kinds of samples of sandstone Y1,Y2 and Y3 change basically logarithmically with increasing confining pressure p, and for all VP, the slopes of the oil-saturated curves are greater than those of the corresponding water-saturated curves;for all H samples, in dry, water- and oil-saturated conditions the VP, VSH, and VSV are larger than the corresponding ones for V samples;for two kinds of sandstones Y1 and Y2, the wave velocity anisotropy parameters ε, γ, and ζ change exponentially or quadratically with increasing confining pressure, and Y2 has a greater correlation coefficient of anisotropy parameters than Y1, but there is no wave velocity anisotropy in Y3. In addition, in dry, water- and oil-saturated conditions the relationships between ε,γ and p or γ, ζ and p are quite distinct from each other, that is, the anisotropy factors ε, γ and ζ of Y1 and Y2 exhibit a clear correlation with pressure and the fluid property, and the ε and γ in oil-saturated condition are larger than that in the corresponding water-saturated condition at the experimental pressure. The experimental results can provide an important basis for the interpretation of seismic data, identifying oil and water zone and comparing with well log data in that region.
Key words: Dry     Water saturated     Oil saturated     Yanchang Formation     Sandstone     Anisotropy    

1 引言

岩石的弹性波速VPVS和波速各向异性是建立地壳和上地幔物质组成和结构模型的基础,对了解地壳、地幔的演化起着重要的作用,是目前认识地球内部物质成分、结构和物理状态最有效的方法,被广泛地应用于区域构造研究和浅部地震勘探中(Christensen and Mooney, 1995Rudnick and Fountain, 1995Weiss et al., 1999Gao et al., 2000Park and Levin, 2002张中杰,2002嵇少丞等, 20062007Barberini et al., 2007Chen et al., 2013彭艳菊等,2007王勤,2007Zhang et al., 2008).

岩石弹性波速及其各向异性的影响因素很多,除了矿物组分外,还有压力、温度、孔隙度、流体类型及微裂隙等(Potter and Mast, 1963Nur and Simmons, 1969Timur,1977Witt and Braun, 1983刘斌等,1998张中杰等,1999Koesoemadinata and McMechan, 2004刘丽娟等,2005苏晓捷等,2005;王玉梅等2007;Crampin and Peacock, 2008Chichinina et al., 2009).在低围压下,某些因素对波速的影响程度甚至远远大于矿物成分的影响.比如在上地壳,有些学者认为孔隙和流体是控制岩石波速的主要因素(王勤,2007;王玉梅等2007).尽管也有学者对此进行了一些局部区域的研究(Batzle and Wang, 1992李爱兵和马丽云,1994Best,1997刘斌等,1998Khaksar et al., 1999史謌等, 20012003卢琳等,2005李洁和李书光,2009肖昆等,2010),然而同一岩性由于其来源地不同,所得到的结果也是不同的,也就是说没有统一的公式适用于任何地区.

陕西延长油田是我国西部重要的油气基地之一,然而目前很少有地层岩石弹性性质方面的报道,仅有作者对部分岩石的纵波、横波速度进行了一些介绍(乔二伟等,2012),但还不够全面,也没有涉及到波速的各向异性等.我们知道,地震资料的解释离不开实验数据的内插和外延.因此,确定各种典型岩石的弹性波性质随压力(围压)、温度等的变化规律,对于地震资料的解释更为重要.所以很有必要对此区域进行更全面和深入的实验研究,从而可以为该地区地震资料的解释、油水层划分及与声波测井之间的对比提供重要的参考依据.鉴于温度对地震波速的影响很小(Timur,1977李爱兵和马丽云,1994卢琳等,2005孙建国,2006肖昆等,2010),本文着重从压力(围压)和含不同孔隙流体方面,通过实验来探讨其对砂岩波速及其各向异性的影响.

2 实验方法 2.1 实验样品

实验所用样品均取自延长油田.延长油田位于鄂尔多斯盆地陕北斜坡东部,区域构造为一平缓的西倾单斜,局部具有差异压实形成的低幅鼻状隆起(翟光明等,1996;李玉宏等2005).延长组是鄂尔多斯盆地在晚三叠世持续拗陷和稳定沉降过程中,堆积的一套河流—三角洲—湖泊相陆源碎屑沉积岩系(武富礼等,2004屈红军等,2011).实验所用的三种砂岩Y1、Y2和Y3,都是来自该地层,它们是三口不同深度的钻井取芯,其取芯深度分别为2007.65~ 2007.78 m、2232.71~2232.86 m和1403.8~1405.67 m. Y1、Y2和Y3的矿物组成和其他有关参数见表 1.

表 1 岩石样品的矿物组成、颗粒大小、密度和孔隙度Table 1 Composition,grain size,density and porosity of the samples

根据实验要求,将岩芯分别沿平行和垂直层面方向(记为H和V)(Best,2007)加工成直径25 mm、高度50 mm 的圆柱体一个或两个,表面抛光,放入温度为90 ℃的真空干燥箱中真空干燥48 h以上,或者真空下用蒸馏水饱和或煤油饱和24 h以上.

2.2 实验装置及测量方法

实验所用装置为New Engl and Research(NER)公司生产的AutoLab 系列(邓继新等,2005He et al., 2010Purcell et al., 2010Mur et al., 2011)2000岩石物性测试系统.其最大围压为200MPa,误差为±0.5%.该仪器波速测量用的是超声波脉冲穿透法,可以同时实现一个纵波和两个相互正交的S波即SH和SV波(He et al., 2010Purcell et al., 2010Mur et al., 2011)的输出结果.对于本次研究的V方向样品,其实只有一个SH波(Best,2007).根据Hornby(1998)的方法,计算出P波和S波的相对误差分别在1.0%和0.77%以下.

实验岩样饱水用的是去离子水、饱油用的是煤油.在室温条件下,实验对每种岩样(包括H和V方向的岩样)按干燥、饱水和饱油的先后顺序分别测量其在围压作用下的纵波、横波速度.一般从5 MPa或10 MPa开始,逐步增加围压直至180 MPa.在每一围压值进行测量前至少需30 min的时间进行稳定,以确保样品内部的应力分布均匀.

3 结果与讨论 3.1 波速VPVSHVSV

图 1是Y1砂岩H、V方向样品即Y1(H)和Y1(V)在干燥、饱水及饱油条件下其纵波、横波速度随围压变化的情况.可以看出,两种样品的三种波速即VPVSHVSV均随围压增加而基本呈对数曲线形式变化(图中给出了纵波的拟合公式及相关系数;另外,根据压力与深度之间的关系(吴泰然等,2011),也给出了参考深度坐标).对于纵波VP:无论是Y1(H)还是Y1(V),除个别点外,均为饱油的波速最大,饱水的次之,干燥的为最小,即饱油>饱水>干燥.这与施行觉等(1998)对大庆砂岩的实验结果基本一致.另外,从图中给出的拟合结果也可看出饱油曲线的斜率远大于饱水,其差值分别为54(图 1a)和123(图 1b).对于横波VSHVSV:则是:饱油≥干燥>饱水.

图 1 Y1岩样在干燥、饱水及饱油情况下VPVSHVSV与围压p的关系(a)H方向;(b)V方向.Fig. 1 P- and S-wave velocity in dry,water and oil saturated conditions changing with confining pressure for(a)sample Y1(H) and (b)sample Y1(V)

上述结果可用以下理论进行解释.

根据弹性波速度公式(马沃克,2008)可知,样品中的饱和水或油会从两方面影响其中传播的纵波、横波波速:一是岩石内部孔隙充满水或油会增大样品的有效体积模量和剪切模量,使波速增大;二是 孔隙内充水或油会加大样品的密度值,使波速变小.

VP来说,低围压下饱和水或油会使样品有效弹性模量的增加起主要作用.根据弹性波速度公式(马沃克,2008),这时的VP均比干燥的大,它们所对应的曲线也在干燥之上.随着围压增加,密度增大的效应逐渐显著.同样由弹性波速度公式(马沃克,2008)可知,VP增大的幅度会逐渐缩小,与干燥之间的差距也在缩小,三条波速曲线逐渐靠近.当围压增加到约130 MPa时,饱水与干燥曲线相交;150 MPa后,干燥曲线则明显位于饱水之上,并向上继续逼近饱油曲线或与之汇合(图 1a,b).由于油的黏滞系数比水大,密度比水小,相对而言,油饱和使样品的有效弹性模量增加更显著,密度增大效应则较弱,所以饱油条件下的VP均比饱水的大.

VSHVSV的情况相对复杂.水饱和使样品密度增大起主要作用,由弹性波速度公式(马沃克,2008)及油水物理性质的差异可知,此时的VSHVSV和饱油及干燥状态相比是最小的.而饱油条件在低围压下是弹性模量增加效应起主要作用,其VSHVSV均比相对应的干燥的大;随着压力上升,密度增大效应开始逐渐显著,两者之间的差距也越来越小,两种曲线逐渐接近并近似平行(图 1a),或者干燥曲线超越饱油曲线,位于其之上(图 1b).另外,高围压尤其是100 MPa以上饱水条件下的VSHVSV基本没有差异,也是密度增大效应作用的结果.

另外,从图 1a图 1b的对比可以看出,H方向不同孔隙流体条件下的三种波速(除饱水VSH外)基本上均大于对应的V方向,这从侧面表明Y1的波速存在各向异性.从图 2还可以看出,除饱水VSH外,Y1(H)其他条件下的VPVSHVSV曲线近似平行于对应的Y1(V),说明尽管Y1存在波速各向异性,但两个方向上波速随围压变化的规律性却基本相同.

图 2 不同围压及不同孔隙流体条件下Y1(H)和Y1(V)岩样VPVSHVSV的对比(a)VP;(b)VSHVSV.Fig. 1 Comparison of measured(a)P- and (b)S-wave velocity for sample Y1(H) and corresponding to those for YI(V)with different pore fluids under different confining pressures

图 3是Y2砂岩在干燥、饱水及饱油条件下其纵波、横波速度随围压变化的情况.不难看出,它们 均随压力增加而以对数曲线形式变化,且纵波饱油曲线的斜率远大于饱水(其差值分别为87(图 3a)和116(图 3b)),这和Y1相似;但在三种波速间的关系方面,却有不尽相同的表现形式.对VP来说,其表现形式基本是:饱水>干燥>饱油,而VSHVSV则是:干燥>饱水>饱油.很明显,由于饱油波速的显著下降而导致了三者间不同于Y1的特征变化.其主要原因可能有以下两个:

图 3 Y2岩样在干燥、饱水及饱油情况下VPVSHVSV与围压p的关系(a)H方向;(b)V方向.Fig. 3 P- and S-wave velocity in dry,water and oil saturated conditions changing with confining pressure for(a)sample Y2(H) and (b)sample Y2(V)

一是Y2和Y1类砂岩的粒度不同.从表 1图 4可以看出,Y2以细沙为主,区别于以中砂为主的Y1.较细颗粒比粗粒所形成的孔隙更狭小,也更容易使充填的油被压实,从而促使密度增大的效果更明显,进而导致Y2饱油波速随压力增加而上升的幅度低于饱水和干燥.但饱水和干燥波速并不受粒 度变化的影响,因此,它们之间的关系和对应的Y1相比,基本保持不变.

图 4 实验样品显微照片(a)Y1样品;(b)Y2样品. 标有“Qtz、Fsp、Bt和Deb”的矿物颗粒分别是石英、长石、黑云母和岩屑.Fig. 4 Microscopic photographs for sample Y1(a) and Y2(b)The mineral grains marked with Qtz,Fsp,Bt and Deb are quartz,feldspar,biotite and debris,respectively.

二是含油饱和度可能较低.较低的含油饱和度使样品有效弹性模量增加的效果不明显,进而影响到其波速的增加.而这一切归根到底是由Y2的细粒结构所形成的狭小孔隙引起的,因为狭小孔隙不利于黏度相对较大的油进入,因此同样的岩样,含油和含水饱和度相比会较低.

另外,同Y1的情况类似,H方向的三种波速均大于对应的V方向,这从侧面表明Y2的波速也存在各向异性.从图 5还可以看出,Y2(H)的VPVSHVSV曲线近似平行于对应的Y2(V),这也说明Y2尽管存在波速各向异性,但两个方向上波速随围压变化的规律性却基本相同.

图 5 不同围压及不同孔隙流体条件下Y2(H)和Y2(V)岩样VPVSHVSV的对比(a)VP;(b)VSHVSV.Fig. 5 Comparison of measured(a)P- and (b)S-wave velocity for sample Y2(H) and corresponding to those for Y2(V)with different pore fluids under different confining pressures

图 6是Y3砂岩在干燥、饱水及饱油条件下VPVSHVSV随压力变化的情况.很明显,它们均随压力增加而以对数曲线形式变化,这点同Y2和Y1相似.另外,还有一个重要的相似点或者说普遍存在于Y1、Y2和Y3的规律是纵波饱油曲线的斜率远大于饱水(图 6),这可作为区分油水性质即识别油水层的重要标志.

图 6 Y3岩样在干燥、饱水及饱油情况下VPVSHVSV与围压p的关系(a)H方向;(b)V方向.Fig. 6 P- and S-wave velocity in dry,water and oil saturated conditions changing with confining pressure for(a)sample Y3(H) and (b)sample Y3(V)

不过,在三种波速的关系方面,Y3则有自己的独特之处. 对VP来说,其表现特征为:饱油≈干燥>饱水;而VSHVSV则是:干燥>饱油>饱水.其内在原因可能与Y3具有较大的孔隙度(表 1)有关.因为岩样充水或充油后,较大孔隙度会更容易使其密度增大的效应起主导作用.这样,随着压力上升,其波速增加趋势也会大幅减缓.由于水比油密度大,因此它对波速的影响也会更突出.所以无论对VPVSHVSV来说,饱水是三种波速中最小的(图 6).另外,Y3具有的较大孔隙度也影响到了其H和V方向波速间的关系,使得这两个方向的VPVSHVSV曲线并不像Y1和Y2那样近似平行,而是相互交织在一起(图 7).这从侧面表明Y3的波速已基本不存在各向异性,下面对此还要做进一步说明.

图 7 不同围压及不同孔隙流体下Y3(H)和Y3(V)岩样VPVSHVSV的对比(a)VP;(b)VSHVSV.Fig. 7 Comparison of measured(a)P- and (b)S-wave velocity for sample Y3(H) and corresponding to those for Y3(V)with different pore fluids under different confining pressures

此外,从图 1—3图 5—7也可以看出,所有的纵波、横波在0~50 MPa段,波速随压力增加较快,明显呈非线性变化,这是由于裂隙和孔隙的闭合引起的(Birch,1961Carlson and Gangi, 1985郝晋昇等,1985刘永祥等,1995杨树锋等,1997Kern et al., 2001李华等,2005Higo et al., 2006马中高等,2006David and Zimmerman, 2012Sun et al., 2012).虽然不同类型砂岩增加的幅度不尽相同,但这一变化却是十分显著的;在压力超过50 MPa后,由于岩石的裂隙和孔隙基本闭合,波速增加呈斜率不同的线性曲线,且变化非常小.由于总体上波速随压力呈对数曲线变化,为方便说明,本文没有按压力进行分段拟合.

波速在0~50 MPa段的非线性特征也直接地影响了下面要介绍的波速各向异性的变化.

3.2 波速各向异性系数ε、γ和ζ

根据Thomsen(1986)Best等(2007),P波和S波各向异性系数可以表示如下:

其中,参数ζ仅作为补充,它只是对水平岩石样品的测量.

这里的P、SH和SV 分别指P波、水平极化S波和垂直极化S波.

图 8是Y1的波速各向异性系数ε、γ和ζ在不同孔隙流体条件下随围压变化的情况.可以看出,它 们基本上随围压的增加而呈下降趋势,但9条曲线拟合的情况并不相同.其中饱水及饱油条件下的6 条曲线均采用二次曲线拟合,相关系数都在90%以下;尤其是饱水条件下的γ,不仅相关系数较低(52.6%), 且其值几乎全部为负值,说明饱水条件下的S波在Y1的两个方向上的各向异性较弱.其主要原因是由 于S波受到了密度增大效应的影响,这在前面已分 析过,此处不再赘述.相对而言,干燥条件下的三条曲线规律性较好,拟合程度也较高(式(4)—(6)).另 外,除ζ外,饱油条件的ε、γ 在较高压力下(>50 MPa)均大于对应的干燥和饱水条件(0~50 MPa段主要受裂隙和孔隙的影响),这点类似于衰减系数的各向异性(史謌和邓继新,2005),且在实验压力范围内,饱油条件的ε、γ 均大于对应的饱水条件,这可以作为区分流体性质即识别油水层的判定标准;但同一孔隙流体下,ε、γ和ζ 的关系比较复杂,其中:干燥为ζ>γ>ε,饱水为ε>ζ>γ,饱油为γ>ε>ζ.

图 8 不同孔隙流体条件下Y1岩样的各向异性系数ε、γ和ζ随不同围压的变化Fig. 8 The wave velocity anisotropy parameters ε、γ and ζ for sample Y1 changing with different pore fluids under different confining pressures

为了更清楚地展现ε、γ和ζ在干燥、饱水及饱油三种条件下随围压变化的规律性,绘制了如图 9的三维图,并可以获得如下的关系式:

图 9 不同孔隙流体条件下Y1岩样的各向异性系数ε、γ和ζ随围压变化的关系图Fig. 9 The relation between wave velocity anisotropy parameters ε、γ and ζ for sample Y1 changing with different pore fluids under different confining pressures

干燥条件

饱水条件

饱油条件

可以看出,ε、γ和围压之间以及γ、ζ和围压之 间呈明显的线性关系;其中以干燥条件下的拟合关 系最好,相关系数都在95%以上;较高压力下(>50 MPa)饱油、干燥及饱水曲线基本上完全分离.很显然,各向异性系数ε、γζ之间表现出明显的与压力和流体状态相关性.

图 10是Y2的波速各向异性系数ε、γζ在不同孔隙流体条件下随围压变化的情况.可以看出,它们基本上也是随着围压增加而呈下降趋势.这同Y1相似,但也有明显不同之处,即9条曲线都为幂函数 关系,且相关系数基本都在90%以上(除干燥下的ζ外).另外,饱油条件的ε、γζ在较高压力下(>50 MPa)均大于对应的饱水和干燥条件,这点也类似于衰减系数的各向异性,且在实验压力范围内,饱油条件的ε、γ均大于对应的饱水条件,这可以作为区分流体性质即识别油水层的判定标准;同一孔隙流体下,ε、γζ则均表现为ε>γ>ζ.这不同于Y1,可能跟Y2具有较细的粒度有关,从而影响到其P波和S波的波速(前面已详细分析过),并进而影响到波速的各向异性系数.具体拟合的情况如下:

图 10 不同孔隙流体条件下Y2砂岩的各向异性系数ε、γ和ζ随不同围压的变化 Fig. 10 The wave velocity anisotropy parameters ε、γ and ζ for sample Y2 changing with different pore fluids under different confining pressures

干燥条件

饱水条件

饱油条件

图 11是Y2的ε、γ和ζ在干燥、饱水及饱油三 种条件下随围压变化的三维图.可以看出,各向异性系数ε、γ和ζ之间表现出明显的与压力和流体状态相关性.

图 11 不同孔隙流体条件下Y2岩样各向异性系数ε、γ和ζ随围压变化的关系图Fig. 11 The relation between wave velocity anisotropy parameters ε、γ and ζ for sample Y2 changing with different pore fluid under different confining pressures

拟合后,得到如式(22)—(27)的关系式.可以看出,其总体规律性要好于Y1,这同样可能与Y2具有较细的粒度有关.

干燥条件

饱水条件

饱油条件

图 12是Y3的波速各向异性系数ε、γ、ζ在不同孔隙流体下随围压变化的情况.可以看出,它们基本上没有规律性,且明显不同于Y1和Y2.除极个别点外,几乎所有数据都为负值且其绝对值在1左右.负值的原因,一方面跟计算方法有关,因为并不是采用常规的最大与最小波速之差来计算(Birch,1961Babuška,1984);另一方面也可能跟所处的沉积环境有关,因为Y3位于河流相而不同于Y1、Y2的湖泊和三角洲相(李玉宏等,2005),但目前并没有这方面的相关报道.

图 12 不同孔隙流体条件下Y3砂岩的各向异性系数ε、γ和ζ随不同围压的变化Fig. 12 The wave velocity anisotropy parameters ε、γ and ζ for sample Y3 changing with different pore fluids under different confining pressures

综合以上Y3的波速及各向异性系数的分析,说明Y3的波速各向异性较弱或已基本没有各向异性.

4 结论

由本实验可以得出如下认识:

(1)三种砂岩Y1、Y2和Y3两个方向即H和V上的VPVSHVSV均随围压增加而基本呈对数曲线形式变化,其相关系数都在90%以上(除Y3的干燥VP外),且纵波饱油曲线的斜率远大于饱水,这可作为区分油水性质即识别油水层的重要标志.

(2)干燥、饱水和饱油波速之间的关系因砂岩种类的不同而不同,这主要取决于砂岩的力学参数(有效弹性模量)、内部结构(粒度、孔隙度)及流体性质(密度、黏滞系数)等.

(3)Y1、Y2砂岩H方向的VPVSHVSV在干燥、饱水和饱油条件下基本上均大于对应的V方向的波速,且对应的波速具有基本相同的规律性.而Y3则没有上述特征.

(4)干燥、饱水和饱油条件下Y1、Y2的波速各向异性系数ε、γ和ζ和围压之间呈指数或二次函数降低.以Y2的各向异性规律最明显,其次是Y1,而Y3基本上没有各向异性;ε、γ和围压之间以及γ、ζ和围压之间呈明显的线性关系.其中以干燥条件下的拟合关系最好,相关系数都在95%以上.

(5)Y1、Y2的各向异性系数ε、γ和ζ之间表现出明显的与压力和流体状态相关性,且在实验压力范围内,饱油条件下的ε、γ均大于对应的饱水条件,这也可以作为区分流体性质即识别油水层的判定标准.

(6)VPVSHVSVε、γ和ζ随围压和不同孔隙流体条件的变化规律,可以为该地区地震资料的解释、油水层划分及与声波测井之间的对比提供基础数据和参考依据,同时也为深部地球物理探测数据的正确解译提供依据.

致谢 感谢地质力学所王红才研究员、赵卫华和孙东生助理研究员等在实验过程中给予的帮助,感谢匿名审稿专家对本文提出的宝贵意见及建议!

参考文献
[1] Babuška V. 1984. P-wave velocity anisotropy in crystalline rocks. Geophysical Journal International, 76(1): 113-119.
[2] Barberini V, Burlini L, Zappone A. 2007. Elastic properties, fabric and seismic anisotropy of amphibolites and their contribution to the lower crust reflectivity. Tectonophysics, 445(3-4): 227-244.
[3] Batzle M, Wang Z J. 1992. Seismic properties of pore fluids. Geophysics, 57(11): 1396-1408.
[4] Best A I. 1997. The effect of pressure on ultrasonic velocity and attenuation in near-surface sedimentary rocks. Geophysical Prospecting, 45(2): 345-364.
[5] Best A I, Sothcott J, McCann C. 2007. A laboratory study of seismic velocity and attenuation anisotropy in near-surface sedimentary rocks. Geophysical Prospecting, 55(5): 609-625.
[6] Birch F. 1961. The velocity of compressional waves in rocks to 10 kilobars, part 2. Journal of Geophysical Research, 66(7): 2199-2224.
[7] Carlson R L, Gangi A F. 1985. Effect of cracks on the pressure dependence of p wave velocities in crystalline rocks. Journal of Geophysical Research, 90(B10): 8675-8684.
[8] Chen Y, Zhang Z J, Sun C Q, et al. 2013. Crustal anisotropy from Moho converted Ps wave splitting analysis and geodynamic implications beneath the eastern margin of Tibet and surrounding regions. Gondwana Research, 24(3-4): 946-957.
[9] Chichinina T I, Obolentseva I R, Ronquillo-Jarillo G. 2009. Anisotropy of seismic attenuation in fractured media: theory and ultrasonic experiment. Transport in Porous Media, 79(1): 1-14.
[10] Christensen N I, Mooney W D. 1995. Seismic velocity structure and composition of the continental crust: A global view. J. Geophys. Res., 100(B6): 9761-9788.
[11] David E C, Zimmerman R W. 2012. Pore structure model for elastic wave velocities in fluid-saturated sandstones. Journal of Geophysical Research, 117(B7), doi: 10.1029/2012JB009195.
[12] Deng J X, Wang S X, Yu J. 2005. Experimental results in partially saturated sandstone under condition of different pressure and their theoretical interpretation. Oil Geophysical Prospecting (in Chinese), 40(5): 530-534.
[13] Crampin S, Peacock S. 2008. A review of the current understanding of seismic shear-wave splitting in the Earth's crust and common fallacies in interpretation. Wave Motion, 45(6): 675-722.
[14] Gao S, Kern H, Liu Y S, et al. 2000. Measured and calculated seismic velocities and densities for granulites from xenolith occurrences and adjacent exposed lower crustal sections: a comparative study from the North China craton. J. Geophys. Res., 105(B8): 18965-18976.
[15] Hao J S, Lu X H, Li J H. 1985. An experimental study on the density and elastic wave velocity of typical rocks of north china under high pressure up to 10kb. Acta Seismologica Sinica (in Chinese), 1(3): 276-284.
[16] He T, Zou C C, Pei F G, et al. 2010. Laboratory study of fluid viscosity induced ultrasonic velocity dispersion in reservoir sandstones. Applied Geophysics, 7(2): 114-126.
[17] Higo Y, Inoue T, Li B S, et al. 2006. The effect of iron on the elastic properties of ringwoodite at high pressure. Physics of the Earth and Planetary Interiors, 159(3-4): 276-285.
[18] Hornby B E. 1998. Experimental laboratory determination of the dynamic elastic properties of wet, drained shales. Journal of Geophysical Research, 103(B12): 29945-29964.
[19] Ji S C, Wang Q, Wang Q, et al. 2007. Correlation between elastic properties and density for Sulu-Dabie UHP metamorphic rocks. Acta Petrologica Sinica (in Chinese), 23(12): 3054-3064.
[20] Ji S C, Wang Q, Denis M, et al. 2006. Pressure-dependence of P-wave velocities in Sulu UHP metamorphic rocks. Acta Geologica Sinica (in Chinese), 80(12): 1807-1812.
[21] Kern H, Popp T, Gorbatsevich F, et al. 2001. Pressure and temperature dependence of VP and VS in rocks from the superdeep well and from surface analogues at Kola and the nature of velocity anisotropy. Tectonophysics, 338(2): 113-134.
[22] Khaksar A, Griffiths C M, McCann C. 1999. Compressional- and shear-wave velocities as a function of confining stress in dry sandstones. Geophysical Prospecting, 47(4): 487-508.
[23] Koesoemadinata A P, McMechan G A. 2004. Effects of diagenetic processes on seismic velocity anisotropy in near-surface sandstone and carbonate rocks. Journal of Applied Geophysics, 56(3): 165-176.
[24] Li A B, Ma L Y. 1994. Effects of frequency and temperature on shear velocity and Q-s in sandstone. Chinese Journal of Geophysics (in Chinese), 37(5): 653-658.
[25] Li H, Wang L S, Li C, et al. 2005. Experimental study on wave velocity of rocks from tarim basin at high PT conditions and its geological implications. Geological Journal of China Universities (in Chinese), 11(4): 601-607.
[26] Li J, Li S G. 2009. Characteristics and changes of rock acoustic velocity under the temperature and pressure of natural reservoir. Chinese Journal of Engineering Geophysics (in Chinese), 6(6): 765-768.
[27] Li L J, Yan G J, Chen J W. 2005. Influences of fluid on velocities of seismic wave. Marine Geology Letters (in Chinese), 21(9): 8-12.
[28] Li Y H, Zhang H A, Wei X Y, et al. 2005. Sedimentary characteristics and distribution of oil and gas in Chang 6 in the Yanchang oilfield. Journal of Northwest University (Natural Science Edition) (in Chinese), 35(5): 606-610, 616.
[29] Liu Y X, Wu F Y, Zhang S H. 1995. The characteristics of elastic wave velocity of rocks at high t and p: implications for the study of deepseated geology. Progress in Geophys. (in Chinese), 10(3): 82-94.
[30] Liu B, Kern H, Popp T. 1998. Velocities and attenuation of P- and S-waves in dry and wet rocks with different porosities under different confining pressures. Chinese J. Geophys.(in Chinese), 41(4): 537-546.
[31] Lu L, Yan G J, Chen J W. 2005. Influences of temperature and pressure on velocities of seismic wave. Marine Geology Letters (in Chinese), 21(9): 13-16, 21.
[32] Ma Z G, Wu X Y, Wang Z H. 2006. Effect of effective pressure on compressional and shear wave velocities. Progress in Exploration Geophysics (in Chinese), 29(3): 183-186.
[33] Mavko G, Mukerji T, Dvorkin J. 2008. The Rock Physics Handbook: Tools for Seismic Analysis of Porous Media (in Chinese). Translated by Xu H B, Dai J C. Hefei: Press of University of Science and Technology of China, 45.
[34] Mur A, Purcell C, Soong Y, et al. 2011. Integration of core sample velocity measurements into a 4D seismic survey and analysis of SEM and CT images to obtain pore scale properties. Energy Procedia, 4: 3676-3683.
[35] Nur A, Simmons G. 1969. The effect of saturation on velocity in low porosity rocks. Earth Planet. Sci. Lett., 7(2): 183-193.
[36] Park J, Levin V. 2002. Seismic anisotropy: tracing plate dynamics in the mantle. Science, 296(5567): 485-489.
[37] Potter P E, Mast R F. 1963. Sedimentary structures, sand shape fabrics, and permeability. The Journal of Geology, 71(4): 441-471.
[38] Purcell C, Mur A, Soong Y, et al. 2010. Integrating velocity measurements in a reservoir rock sample from the SACROC unit with an AVO proxy for subsurface supercritical CO2. The Leading Edge, 29(2): 192-195.
[39] Qiao E W, Zhao W H, Long C L. 2012. A laboratory study of seismic wave velocity in sandstone bearing differential pore fluids at different pressures. Chinese J. Geophys.(in Chinese), 2012, 55(12): 4257-4265.
[40] Qu H J, Yang X C, Cao J Z, et al. 2011. Oil accumulation rules in deep zones of Upper Triassic Yanchang Formation in Ordos Basin. Acta Petrolei Sinica (in Chinese), 32(2): 243-248.
[41] Rudnick R L, Fountain D M. 1995. Nature and composition of the continental crust: A lower crustal perspective. Rev. Geophys., 33(3): 267-309.
[42] Shi G, Deng J X. 2005. Experimental study on the attenuation anisotropy of shale and mudstone under reservoir condition. Science in China (Series D) (in Chinese), 35(3): 268-275.
[43] Shi G, Shen W L, Yang D Q. 2003. The relationship of wave velocities with saturation and fluid distribution in pore space. Chinese J. Geophys.(in Chinese), 46(1): 138-142.
[44] Shi G, Yang D Q. 2001. The regression analysis study on velocity and porisity, and clay content of rocks. Acta Scientiarum Naturalium Universitatis Pekinensis (in Chinese), 37(3): 379-384.
[45] Shi X J, Xia C J, Wu Y G. 1998. The laboratory study on wave velocity under reservoir condition and its affection factors. Chinese J. Geophys.(in Chinese), 41(2): 234-241.
[46] Su X J, Yan G J, Chen J W, et al. 2005. Influences of rock pore and clay content on velocities of seismic. Marine Geology Letters (in Chinese), 21(9): 3-7.
[47] Sun J G. 2006. Rock Physics (in Chinese). Beijing: Geology Press, 112.
[48] Sun S S, Ji S C, Wang Q, et al. 2012. Seismic velocities and anisotropy of core samples from the Chinese Continental Scientific Drilling borehole in the Sulu UHP terrane, eastern China. Journal of Geophysical Research, 117, doi: 10.1029/2011JS008672.
[49] Thomsen L. 1986. Weak elastic anisotropy. Geophysics, 51(10): 1954-1966.
[50] Timur A. 1977. Temperature dependence of compressional and shear wave velocities in rocks. Geophysics, 42(5): 950-956.
[51] Wang Q. 2007. Experimental and theoretical study on seismic properties of rocks. Bulletin of Mineralogy Petrology and Geochemistry (in Chinese), 26(2): 118-126.
[52] Wang Y M, Miao Y K, Meng X J, et al. 2007. The impact of the fluid properties on the seismic wave velocities in rocks. Petroleum Geophysics (in Chinese), 5(3): 36-39.
[53] Weiss T, Siegesmund S, Rabbel W, et al. 1999. Seismic velocities and anisotropy of the lower continental crust: A review. Pure Appl. Geophys., 156(1-2): 97-122.
[54] Witt K J, Braun S J. 1983. Permeability-anisotropy due to particle-shape. Journal of Geotechnical Engineering, 109(9): 1181-1187.
[55] Wu F L, Li W H, Li Y H, et al. 2004. Delta sediments and evolution of the Yanchang Formation of Upper Triassic in Ordos Basin. Journal of Palaeogeography (in Chinese), 6(3): 307-315.
[56] Wu T R, He G Q. 2011. Physical Geology (in Chinese). 2nd ed. Beijing: Peking University Press, 27.
[57] Xiao K, Zou C C, Nie X, et al. 2010. Experimental research of S-wave and P-wave velocities for rocks under the reservoir conditions in G sag, Jiangsu Oilfield. Chinese Journal of Engineering Geophysics (in Chinese), 7(6): 688-693.
[58] Yang S F, Cheng H L, Jiang J S, et al. 1997. The velocity characteristics and its geological significance of I and S type granite from Southern China under high temperature and high pressure. Science in China (Series D) (in Chinese), 27(1): 33-38.
[59] Zhai G M. 1996. Petroleum geology of China (in Chinese). Vol. 12. Beijing: Petroleum Industry Press, 67.
[60] Zhang Z J. 2002. A review of the seismic anisotropy and its applications. Progress in Geophysics (in Chinese), 17(2): 281-293.
[61] Zhang Z J, Teng J W, He Z H. 1999. Azimuth anisotropic studies of Seismic velocity, attenuation and quality factor in EDA medium. Science in China (Series E) (in Chinese), 29(6): 569-574.
[62] Zhang Z J, Zhang X, Badal J. 2008. Composition of the curst beneath southeastern China derived from an integrated geophysical data set. Journal of Geophysical Research-Solid Earth, 113(B4), doi: 10.1029/2006JB004503.
[63] 邓继新, 王尚旭, 俞军. 2005. 不同压力条件下部分饱和砂岩速度实验结果及理论解释. 石油地球物理勘探, 40(5): 530-534.
[64] 葛瑞·马沃克, 塔基·木克基, 杰克·德沃金. 2008. 岩石物理手册: 孔隙介质中地震分析工具. 徐海滨, 戴建春译. 合肥: 中国科学技术大学出版社, 45.
[65] 郝晋昇, 刘晓红, 李纪汉. 1985. 华北地区典型岩石在1万巴压力下的密度和弹性波速度的实验研究. 地震学报, 1(3): 276-284.
[66] 嵇少丞, 王茜, 王勤等. 2007. 苏鲁—大别超高压变质岩的弹性力学性质与密度的关系. 岩石学报, 23(12): 3054-3064.
[67] 嵇少丞, 王茜, Denis M等2006. 苏鲁超高压变质岩中地震波速随围压的变化规律. 地质学报, 80(12): 1807-1812.
[68] 李爱兵, 马丽云. 1994. 频率、温度对砂岩的横波速度和Qs值的影响. 地球物理学报, 37(5): 653-658.
[69] 李华, 王良书, 李成等. 2005. 塔里木盆地岩石高温高压下波速的实验研究及其地质意义. 高校地质学报, 11(4): 601-607.
[70] 李洁, 李书光. 2009. 储层温压条件下岩石的波速特性及变化规律. 工程地球物理学报, 6(6): 765-768.
[71] 刘丽娟, 闫桂京, 陈建文. 2005. 流体对地震波速度的影响. 海洋地质动态, 21(9): 8-12.
[72] 李玉宏, 张化安, 魏仙样等. 2005. 延长油田长6油层组沉积特征与油气分布. 西北大学学报(自然科学版), 35(5): 606-610, 616.
[73] 刘永祥, 吴福元, 张世红. 1995. 高温高压下岩石弹性波速特征及其在深部地质研究中的意义. 地球物理学进展, 10(3): 82-94.
[74] 刘斌, Kern H, Popp T. 1998. 不同围压下孔隙度不同的干燥及水饱和岩样中的纵波、横波速度及衰减. 地球物理学报, 41(4): 537-546.
[75] 卢琳, 闫桂京, 陈建文. 2005. 地层温度和压力对地震波度的影响. 海洋地质动态, 21(9): 13-16, 21.
[76] 马中高, 伍向阳, 王中海. 2006. 有效压力对岩石纵横波速度的影响. 勘探地球物理进展, 29(3): 183-186.
[77] 彭艳菊, 黄忠贤, 苏伟等. 2007. 中国大陆及邻区海域地壳上地幔各向异性研究. 地球物理学报, 50(3): 752-759. (本条文献无对应的英文信息)
[78] 乔二伟, 赵卫华, 龙长兴. 2012. 含不同孔隙流体的砂岩地震波速度随压力变化的实验研究. 地球物理学报, 55(12): 4257-4265.
[79] 屈红军, 杨县超, 曹金舟等. 2011. 鄂尔多斯盆地上三叠统延长组深层油气聚集规律. 石油学报, 32(2): 243-248.
[80] 史謌, 邓继新. 2005. 地层条件下泥、页岩衰减各向异性实验研究. 中国科学D辑, 35(3): 268-275.
[81] 史謌, 沈文略, 杨东全. 2003. 岩石弹性波速度和饱和度、孔隙流体分布的关系. 地球物理学报, 46(1): 138-142.
[82] 史謌, 杨东全. 2001. 岩石波速和孔隙度、泥质含量之间的关系研究. 北京大学学报(自然科学版), 37(3): 379-384.
[83] 施行觉, 夏从俊, 吴永钢. 1998. 储层条件下波速的变化规律及其影响因素的实验研究. 地球物理学报, 41(2): 234-241.
[84] 苏晓捷, 闫桂京, 陈建文等. 2005. 岩石孔隙及泥质含量对地震波速度的影响. 海洋地质动态, 21(9): 3-7.
[85] 孙建国. 2006. 岩石物理学基础. 北京: 地质出版社, 112.
[86] 王勤. 2007. 岩石地震波性质的实验与理论研究. 矿物岩石地球化学通报, 26(2): 118-126.
[87] 王玉梅, 苗永康, 孟宪军等. 2007. 流体性质对岩石地震波速度的影响. 油气地球物理, 5(3): 36-39.
[88] 武富礼, 李文厚, 李玉宏等. 2004. 鄂尔多斯盆地上三叠统延长组三角洲沉积及演化. 古地理学报, 6(3): 307-315.
[89] 吴泰然, 何国琦. 2011. 普通地质学(第二版). 北京: 北京大学出版社, 27.
[90] 肖昆, 邹长春, 聂昕等. 2010. 江苏油田G凹陷储层条件下岩石纵波、横波速度实验研究. 工程地球物理学报, 7(6): 688-693.
[91] 杨树锋, 陈汉林, 姜继双等. 1997. 高温高压下华南Ⅰ型和S型花岗岩的波速特征及其地质意义. 中国科学(D辑), 27(1): 33-38.
[92] 翟光明. 1996. 中国石油地质志. 卷12. 北京: 石油工业出版社, 67.
[93] 张中杰. 2002. 地震各向异性研究进展. 地球物理学进展, 17(2): 281-293.
[94] 张中杰, 滕吉文, 贺振华. 1999. EDA介质中地震波速度、衰减与品质因子方位异性研究. 中国科学(E辑), 29(6): 569-574.