-
砂砾岩体,是指主要由粒径大于2 mm的砾石和砂岩一并构成的粗粒碎屑沉积体,其内部砾岩、砂岩和泥质杂基含量不一[1⁃2],常见冲积扇、近岸水下扇、扇三角洲、辫状河三角洲、深水浊积扇、滑塌浊积扇等不同成因的扇体在断陷湖盆陡坡带形成快速堆积、有序发育的扇体群[3⁃4]。作为重力流沉积的主要研究对象,砂砾岩体受到全球沉积学界的重点关注[5⁃7],同时也在油气勘探和水文地质勘探领域占据重要地位[1,8⁃10],尤其是准噶尔盆地玛湖凹陷特大型砾岩油藏和渤海湾盆地渤中19-6大型整装凝析气田的连续发现,使得砂砾岩体迅速成为我国油气资源规模增储的接替阵地[1,11⁃14]。由于砂砾岩体多靠近生烃中心,常与深湖相烃源岩直接接触或被其包裹,因而具有优越的成藏条件[1,4,8]。然而,在构造活动强烈、古地形起伏多变、物源近、坡度陡的沉积背景下,砂砾岩体往往表现出沉积动力复杂、扇体迁移叠置、岩性岩相多样、内部非均质性强等突出特征,导致砂体成层性差、储层相变快且预测风险较高[1,15⁃18]。另外,埋藏后的差异性流体―岩石作用过程进一步加剧了储层孔喉结构和油水关系的复杂性[13,19⁃20]。
岩石作为不同沉积、成岩过程的产物,其物理性质必然携带不同地质过程的地震属性信息[13,21]。显然,借助岩石物理实验建立储层岩性、物性、流体与地震属性之间的关系,能够有效约束地震反演,进而提高储层预测精度[22]。在众多属性中,由于地震波频散特征与流体流动密切相关而显得极为重要。当地震波经过岩石储层时,诱发流体流动,导致能量损耗,从而产生地震波频散现象,称为波致流[23⁃26]。Gregory[27]对孔隙度为4%~41%的岩石样品进行了超声纵波速度测量,发现在低孔隙度岩石中,流体饱和度对纵波速度的影响要比在高孔隙度岩石中更加显著。Domenico[28]测量了盐水饱和度对疏松含气砂岩纵波速度的影响,发现在低盐水饱和度时,速度测量值与理论值基本一致,而在高盐水饱和度时,速度测量值与理论结果存在差异,并认为这种差异与波致流现象相关。Toksöz et al.[29]在实验室中使用超声波法测量了盐水饱和的砂岩和灰岩,发现存在纵波衰减,且随着有效压力增加,纵波衰减显著减小。King et al.[30]使用超声透射法测量了砂岩储层样品的纵、横波速度,对于盐水饱和样品,研究发现实验测量结果与Gassmann[31]预测结果不相符,孔隙尺度上的流体流动即喷射流理论可为其测量结果提供合理的解释。Lebedev et al.[32]使用超声波研究了饱和度对砂岩的影响,在低饱和度时的测量结果接近Gassmann-Wood（GW）边界,而在高饱和度时,测量结果接近Gassmann-Hill（GH）边界,在中等饱和度时位于GW和GH边界之间,同样认为这种现象由波致流现象引起。
实验研究表明,频散和衰减普遍存在,但测量频段主要集中在超声频段,难以确定地震频段相应的频散和衰减机制。为实施地震频段测量,Spencer[33]首先利用应力—应变法在4~400 Hz测量了水饱和砂岩的衰减,并使用应力松弛理论解释实验数据。Batzle et al.[34]在地震频段测量了砂岩速度频散和衰减,并定义“流体流动性”控制频散和衰减特征。Mikhaltsevitch et al.[35]对完全水饱和状态下的低渗透砂岩开展地震频段衰减测量,但未解释测量结果。Subramaniyan et al.[36]测量了枫丹白露砂岩的频散和衰减,其实验结果表明喷射流存在。Pimienta et al.[37]在对饱和枫丹白露砂岩的宽频带（0.005~100 Hz,1 MHz）测量中观察到两个衰减峰值,并将其解释为与喷射流和排水/不排水过渡的组合。未晛等[38]利用应力—应变低频岩石物理测量装置研究了孔隙流体对不同渗透率岩石地震波速度的影响,揭示在低饱和度下,致密砂岩在地震和超声频段下没有明显的频散,但在高饱和度下纵波速度的频散变得明显。另外,针对砂岩储层[39]和灰岩储层[40],笔者曾分别测量了二者在地震频段与超声频段的频散和衰减特征并分析了其机理,认为流体在不同尺度的流动是引发地震波频散和衰减的本质原因,为使用地震波速度频散和衰减属性预测流体性质提供了可能。
上述研究获得了大量的频散、衰减数据和重要的理论认识,有效地指导了砂岩和灰岩储层的地球物理勘探工作,但目前缺少对断陷湖盆砂砾岩体岩石物理响应特征的系统总结,其复杂的沉积特征、孔隙结构和非均质性对地震波弹性响应的影响在根本上限制着流体反演方法的可靠性和流体识别结果。因此,本文针对东营凹陷陡坡带典型砂砾岩样品开展跨频段（地震频段和超声频段）弹性参数测量与分析,通过实验结果探究敏感弹性参数,总结频散规律及机制,为砂砾岩储层的岩性和烃类检测提供岩石物理依据。
-
本文测量装置分为地震频带测量单元和超声频带测量单元。地震频带测量单元的有效测量频率范围为1~100 Hz,设备配置与Borgomano et al.[41]研发的设备基本一致。图1显示了该装置的设计图,包括A：轴压加载,用以加载轴向压力；B：压电陶瓷震源,用以激励震源产生正弦信号；C：胡克腔,用以加载围压,可加载到60 MPa；D：加温套,用来加载温度,可达150 ℃；E：底座,用于调整高度；F：应变片,用以测量岩石样品和铝制样品的轴向和径向应变；G：超声探头,用以测量岩石的超声波速度；H：参考铝样；I：加压流体泵；J：惠斯通电桥,用以识别应变片的微弱信号；K：气瓶。
图 1 多频段岩石物理设备图
Figure 1. Multi⁃frequency band apparatus diagram
地震频带测量单元的工作原理为应力—应变法。在测量过程中,测试样品和铝制参考样品被黏合在一起并沿其轴线垂直对齐,以使应力能够在横向和轴向方向均匀分布。当在测试样品上施加特定频率的正弦应力时,应变片可以检测轴向和径向应变,据此计算岩石模量和泊松比：
(1) (2) 式中：Eal 是参考铝样的杨氏模量,值为72 GPa。εal 为参考铝样的轴向应变,εax 为测试样品轴向应变,εra 为参考铝样的径向应变。基于各向同性假设,可以利用杨氏模量,泊松比和密度计算岩石样品的纵、横波速度,具体参考岩石物理手册第二章相关内容[42]。
超声频带测量单元的工作原理为行波法。P波和S波探头镶嵌在参考铝样内部（图1,H）,上部探头激发频率为1 MHz的信号,下部探头接收该信号,假设样品长度为
(3) -
济阳坳陷东营凹陷北部陡坡带因坡度陡、物源近、构造活动强烈,在沙河街组三段、四段时期发育冲积扇、扇三角洲、近岸水下扇、浊积扇等多种类型的砂砾岩沉积体[2,8,15]。本次实验共钻取沙三段40块岩心样品,颗粒粒度涵盖细砾岩、砂质砾岩、砾质砂岩、细砂岩等各粒级代表性样品。孔隙度的测试结果（图2）显示,大部分样品孔隙度集中在0~14%,少量样品孔隙度大于20%,其中不同岩性孔隙度分布范围重叠；渗透率测试结果（图3）显示,大部分样品的渗透率集中在10-5~10-3 μm2,少量样品渗透率高于10-2 μm2,不同岩性渗透率分布范围重叠交叉且渗透率和孔隙度呈弱正相关关系。整体来看,以扇中分支水道含砾砂岩、细砂岩孔渗条件最优,其次为扇根主水道砾岩,泥岩物性最差。
图 2 沙三段40块样品孔隙度统计结果
Figure 2. Porosity of 40 samples in the Third member of the Shahejie Formation (Es3)
图 3 沙三段40块样品渗透率统计结果
Figure 3. Permeability of 40 samples in in the Third member of the Shahejie Formation (Es3)
在此基础上,选取21块代表性样品开展超声测量,依据岩性分组：砾岩6块,砂岩8块,含砾砂岩5块,泥质砂岩1块,灰岩1块。为揭示砂砾岩体的频散特征,优选出4块典型岩石样品进行地震频带测量,样品信息为：16号砂质砾岩（图4a）、8号细砾岩（图4b）、10号细砂岩（图4c）和13号砾质砂岩（图4d）。样品所在深度、物性参数、矿物组成和颗粒体积模量等信息如表1~3所示,颗粒体积模量K0 通过Voigt-Reuss-Hill平均法[42]计算得出。因4块样品的矿物组分类似（表3）,对应体积模量差异不大。从渗透率的角度分析,细砂岩和砾质砂岩的渗透率较低,而砂质砾岩和细砾岩的渗透率较高；从孔隙度的角度分析,砾质砂岩和砂质砾岩的孔隙度较高,而细砂岩和细砾岩的孔隙度较低。
图 4 用于低频测试的代表性岩心样品
Figure 4. Representative core samples for low⁃frequency measurement
图 5 饱和气、水和油的21块样品的纵横波速度随压力变化曲线
Figure 5. Evolution of P⁃and S⁃wave velocity versus increasing pressure for 21 samples of saturated gas, water, and oil
图 6 砂砾岩饱和水条件下的纵、横波速度频散
Figure 6. P⁃and S⁃wave velocity versus frequency under different effective pressures for glutenite saturated with water
图 7 岩性指示因子
Figure 7. Lithology indicator factors
图 8 流体指示因子
Figure 8. Fluid indicator factors
图 9 裂隙密度随压力变化情况
Figure 9. Crack density versus pressure
图 10 软孔隙随纵横比变化的分布特征
Figure 10. Soft⁃porosity distribution versus aspect ratio
表 1 16号、8号、10号和13号样品岩性
岩心序号 岩性 深度/m 16 砂质砾岩 2 991.5 8 细砾岩 2 713.8 10 细砂岩 2 822.8 13 砾质砂岩 2 744.7 表 2 样品物性参数
岩心序号 密度/（g/cm3） 孔隙度/% 渗透率/×10-3 μm2 16 2.33 13.88 22.23 8 2.40 8.62 17.70 10 2.42 9.91 0.15 13 2.36 12.84 3.52 表 3 16号、8号、10号和13号样品的矿物组分
岩心序号 黏土总量/% 石膏/% 石英/% 钾长石/% 斜长石/% 方解石/% 铁白云石/% 白云石/% 菱铁矿/% K0/GPa 16 6.3 3.3 27.4 0 34.2 15.4 3.7 9.7 0 62.2 8 5.0 3.8 26.6 0 29.4 23.8 3.8 7.7 0 63.2 10 4.7 0 25.8 4.4 16.7 33.3 4.7 7.9 2.5 66.2 13 5.0 3.8 26.6 0 29.4 23.8 3.8 7.7 0 63.2 -
对21块样品分别饱和气、水和油（表4）后进行超声测量,获得纵横波速度（图5）。
表 4 流体物性参数
流体 密度/（g/cm3） 黏度/cP 体积模量/GPa 空气 10-3 10-5 10-3 水 1 1 2.2 油 0.89 68 1.8 对于干燥样品,纵、横波速度特征如下：灰岩具有最高的纵波速度,超过了6 560.0 m/s,而泥岩的纵波速度次之,介于5 500.0~5 700.0 m/s；砾岩和砂岩的纵波速度差异较大,分布范围分别为3 900.0~6 400.0 m/s和2 170.0~6 600.0 m/s。横波速度特征类似,灰岩同样具有最高的横波速度,超过4 050.0 m/s,而泥岩的横波速度最低,基本介于3 260.0~3 350.0 m/s；砾岩和砂岩的横波速度也存在较大的差异,分别为1 990.0~4 000.0 m/s和1 200.0~4 060.0 m/s。
对于饱和水的样品,其纵、横波速度具有以下特点：灰岩的纵波速度最高,超过6 660.0 m/s,泥岩次之,速度介于5 800.0~5 890.0 m/s。砾岩、含砾砂岩和砂岩的纵波速度跨度较大,速度介于3 270.0~6 600.0 m/s。对比干燥样品,饱和水纵波速度整体提升。另一方面,灰岩的横波速度最高,超过4 000.0 m/s,而砾岩的横波速度介于2 000.0~4 000.0 m/s,砂岩的横波速度介于1 200.0~3 700.0 m/s。对比干燥样品可知,流体对横波存在一定影响,但比对纵波的影响幅度小。
对于饱和油的样品,其纵、横波速度具有以下特点：灰岩纵波速度最高,大于6 640.0 m/s；泥岩次之,纵波速度在5 800.0 m/s以上；砾岩的纵波速度介于4 820.0~6 640.0 m/s；砂岩的纵波速度介于3 200.0~6 500.0 m/s。从横波速度看,灰岩横波速度同样最高,超过4 000.0 m/s；泥岩横波速度为3 300.0 m/s左右；砾岩横波速度介于2 000.0~4 000.0 m/s；砂岩横波速度介于1 200.0~4 000.0 m/s。对比饱和水样品,饱和油与饱和水的纵波速度基本相似。
综上,灰岩的速度高于砾岩、泥质砂岩和含砾砂岩,表明当砂岩中存在泥质或砾石时,其非均质性更强,速度差异也更大。
除了流体和岩性影响速度外,压力同样影响速度。由测量结果可知,无论是干燥还是饱和状态,五种岩石的纵、横波速度随着压力的增大而增大,总体呈现非线性增长,其变化率在较高压力处小,在较低压力处大。引起这种变化规律的原因是岩石中存在裂缝等纵横比较小的软（易压缩）孔隙,其随有效压力增加逐渐闭合,导致岩石骨架刚度变大,进而引发相应的纵、横波速度以及弹性模量的增大。由于砂岩和砂砾岩具有较高的孔隙度,需要40~50 MPa的压力才足以使得岩石中软孔隙完全闭合,即速度压力曲线呈现轻微的变缓的趋势。相对而言,泥岩和灰岩孔隙度较低,随着压力的变化规律并不敏感。将岩石完全饱和水或油后,砂岩的纵波速度压力敏感性降低,其原因是可压缩性更小的水替代了孔隙中原有的空气,降低了岩石的整体可压缩性（即岩石刚度提升）。
-
通过杨氏模量和泊松比计算得出16号样品在水饱和状态下的纵波速度和横波速度（图6a,b）。设备所测杨氏模量E的不确定度约为5％,泊松比的不确定度约为4％。基于这些不确定度推导出的纵波速度不确定度约为6%[40]。在1 MPa压力下时,随着频率从1 Hz增加到1 MHz,16号样品纵波速度从2 913.6 m/s增加到3 430.9 m/s,频散幅度约为17%。当压力增加到5 MPa时,随着频率从1 Hz增加到1 MHz,纵波速度从3 066.7 m/s增加到3 546.8 m/s,频散幅度约为15.7%。当压力增加到20 MPa时,纵波速度从3 644.8 m/s增加到3 968.8 m/s,频散幅度约为9%。因此可以明确,随着压力增加,频散幅度减小。通过表1可知,16号样品所在深度为2 991.0 m,根据Avseth et al.[43]中提到的关于有效压力与深度间的换算关系,当孔隙压力为静水压力时,30 MPa的有效压力能够近似岩心所处地层的压力条件。在30 MPa压力下,随着频率从1 Hz到1 MHz,地震波速度从3 860.0 m/s增加到4 059.0 m/s,频散幅度降低到6%,表明在地下原位储层中,16号样品的频散变低。
在水饱和状态下,8号样品的纵波速度和横波速度（图6c,d）的不确定度约为6%[40]。在1 MPa压力下,随着频率从1 Hz增加到1 MHz,纵波速度从4 292.2 m/s增加到4 849.1 m/s,频散幅度约为12%。当压力增加到5 MPa时,随着频率从1 Hz增加到1 MHz,纵波速度从4 514.0 m/s增加到4 985.0 m/s,频散幅度约为10.4%。当压力增加到20 MPa时,纵波速度从4 928.0 m/s增加到5 125.0 m/s,频散幅度约为4%。与16号样品所得规律一致,随着压力增加,频散幅度减小。通过表1可知,8号样品所在深度为2 713.0 m,根据Avseth et al.[43]中提到的关于有效压力与深度间的换算关系,当孔隙压力为静水压力时,30 MPa的有效压力能够近似岩心所处深度地层的压力条件。在30 MPa压力下,随着频率从1 Hz增加到1 MHz,纵波速度从5 103.0 m/s增加到5 202.0 m/s,频散幅度降低到2%,这表明在地下原位储层压力下,8号样品的频散变得很低,基本可以忽略不计。
在水饱和状态下,10号样品的纵波速度和横波速度（图6e,f）的不确定度约为6%[40]。当压力为1 MPa时,随着频率从1 Hz增加到1 MHz,纵波速度从3 960.0 m/s增加到4 369.0 m/s,频散幅度约为10%。当压力增加到5 MPa时,随着频率从1 Hz增加到1 MHz,纵波速度从4 017.0 m/s增加到4 395.0 m/s,频散幅度约为9%。当压力增加到20 MPa时,纵波速度从4 499.0 m/s增加到4 646.0 m/s,频散幅度约为3%。随着压力增加,频散幅度减小。通过表1可知,10号样品所在深度为2 822.8 m,根据Avseth et al.[43]中提到的关于有效压力与深度间的换算关系,当孔隙压力为静水压力时,30 MPa的有效压力能够近似岩心所处2 800.0 m深度地层的压力条件。在30 MPa压力下,随着频率从1 Hz增加到1 MHz,地震波速度从4 645.0 m/s增加到4 723.3 m/s,频散幅度降低到2%,这表明在地下原位储层中,10号样品的频散相对较低,可以忽略不计。
在水饱和状态下,13号样品的纵波速度和横波速度（图6g,h）的不确定度约为6%[40]。在1 MPa压力下,随着频率从1 Hz增加到1 MHz,纵波速度从3 104.1 m/s增加到3 677.7 m/s,频散幅度约为18%。当压力增加到5 MPa时,随着频率从1 Hz增加到1 MHz,纵波速度从3 227.4 m/s增加到3 744.9 m/s,频散幅度约为16%。当压力增加到20 MPa时,纵波速度从3 752.1 m/s增加到4 009.3 m/s,频散幅度约为7%。随着压力增加,频散幅度减小。通过表1可知,13号样品所在深度为2 744.7 m,根据Avseth et al.[43]中提到的关于有效压力与深度间的换算关系,当孔隙压力为静水压力时,30 MPa的有效压力能够近似岩心所处2 744.7 m 深度地层的压力条件。在30 MPa压力下,随着频率从1 Hz增加到1 MHz,地震波速度从3 959.0 m/s增加到4 081.2 m/s,频散幅度降低到3%,这表明在地下原位储层中,13号样品的频散相对较低。
-
在常规岩石物理分析中,寻找敏感参数区分岩性及流体是常用的分析手段。本文采用Dillon et al.[44]提出的统计学方法寻找敏感参数,即构建指示因子区分两组数据,从而获得敏感属性参数：
(4) 式中：I为指示因子,X为特定物相的特定属性的平均值,X′为参考物相的特定属性的平均值。σX′ 为参考物相的特定属性的标准差。当I大于2时,表明两组物相数据的中心点距离大于参考相数据的范围,因此两组物相数据可以区分。在区分砂岩和砾岩时,X为砾岩特定属性平均值,X′为砂岩特定属性平均值,σX′ 为砂岩特定属性标准差。在区分流体时,例如水和气,则
以水为参考相,可计算得到砂岩、砾岩及砂砾岩的流体指示因子（图8a~c）。综合三种岩性,气—水差异明显,主要的敏感弹性参数为v,ρVp/Vs,ρv,λv。然而,油—水差异很小,没有敏感弹性参数可以区分。选取v和ρVp/Vs作为敏感参数时的交会图（图8d）分析砂砾岩储层,可见气—水可明显区分,但油—水无法区分。
综合上述敏感因子分析可知：从岩性上看,砂砾岩与其他岩性较容易区分；从流体上看,砂岩、砾岩及砂砾岩中饱和流体为水—气和油—气容易区分,但油—水之间不易区分。
-
参照邓继新等[45]的研究认识,利用超声频带纵、横波速度计算表1中各样品的裂隙密度：
(5) (6) 式中：Kh,vh 和Gh 是岩石骨架无微裂隙体积模量,泊松比和剪切模量,可通过纵、横波速度高压极限测量结果获得,ε是裂隙密度。Kd 和Gd 是由纵、横波计算所得干燥样品体积模量和剪切模量。图9是裂隙密度随压力的演化结果：在1 MPa压力下,砂质砾岩（图9a）和砾质砂岩（图9d）裂隙密度约为0.9,而砂岩（图9c）和砾岩（图9b）的裂隙密度约为0.2；随着压力增加,所有样品的裂隙密度均降低。上述四块样品所在的深度为2 700.0~2 900.0 m,根据Avseth et al.[43]中提到的关于有效压力与深度间的换算关系,当孔隙压力为静水压力时,20~30 MPa的有效压力能够近似岩心所处深度地层的压力条件。因此,对于砂砾岩（图9a,d）而言,原位压力裂隙密度约0.2,而砾岩（图9b）和砂岩（图9c）的原位压力裂隙密度约0.04,是砂砾岩裂隙密度的1/5。
参照Duan et al.[46]和Sun et al.[47]的研究,基于裂隙密度,可计算岩石中未闭合软孔隙的最小初始纵横比：
(7) 式中：εp 为特定压力的裂隙密度。Kd (εp )为特定压力情况含气体积模量,可通过特定压力的纵、横波速度（图6）计算。ε0 和εp 分别为压力为0和p时的裂隙密度。则样品中软孔隙的纵横比值随压力的变化关系为：
(8) 通过裂隙密度和裂隙纵横比计算软孔隙度Φc 为：
(9) 软孔隙随纵横比分布的演化关系（图10）表明：在1 MPa压力下,砂砾岩（图10a,d）样品中裂隙纵横比分布的最大值为2.2×10-3,而砂岩（图10c）和砾岩（图10d）的裂隙纵横比分布的最大值为1.2×10-3。随着压力增加,裂隙纵横比分布范围逐渐减小,这是因为压力增加使小纵横比的裂隙闭合,同时使大纵横比的微裂隙纵横比变小。当压力达到储层原位压力时（20~30 MPa）,砂砾岩的纵横比分布的最大值为1.2×10-3,砂岩和砾岩的纵横比分布的最大值为0.6×10-3,与此同时,软孔隙度变少,但砂砾岩的软孔隙度约为细砾岩和细砂岩软孔隙度的10倍。
裂隙密度及纵横比与岩石的频散特征密切相关[48⁃49],参照邓继新等[45]计算速度的方法,可计算表1中样品饱和水时纵、横波速度的频散特征（图6,实线）。该预测结果与实验结果完美匹配,表明砂砾岩体频散特征满足“喷射流”机制。基于喷射流理论,相对于水,油的高黏度使流体交换需要更多的时间来达到平衡,特征频率向低频移动,降低频散出现频率范围,从而使地震频带的频散范围增大；反之,可以利用频散范围区分油和水,从而解决传统敏感参数无法区分油水的问题。图6中的加号为20 MPa下饱和油样品的速度,可以看到,与上文分析一致,频散发生的频率范围更低,频散幅度更大。
压力增加将导致裂隙密度（图9）减小,进而改变速度频散特征。因此,探究储层压力下岩石样品的裂隙密度及其对应的频散特征十分必要。不同岩性下速度频散范围（图11a,b）、裂隙密度（图11c）及裂隙闭合压力（图11d）之间的关系显示：砂砾岩纵波速度频散最大,最大值为16%,砂岩次之,最高值为12%,砾岩更小,最大值为10%,泥灰岩最小,最大值为4%；各岩性横波速度频散范围均小于5%,其中砂砾岩横波频散范围偏大,泥灰岩频散最小,几乎可以忽略；砂砾岩体裂隙密度最大,砂岩次之,砾岩和泥灰岩裂隙密度较小；砂砾岩微裂隙闭合压力介于40~45 MPa,砂岩的闭合压力介于20~40 MPa,泥灰岩需要的压力最小,介于6~10 MPa。
图 11 砂岩、砾岩、砂砾岩及泥灰岩综合分析
Figure 11. Comprehensive analysis of sandstone, conglomerate, glutenite, and mudstone⁃limestone
总体而言,砂砾岩裂隙发育,其速度频散特征与泥灰岩存在明显差异且具有流体敏感性。因此,传统流体反演方法,如常用的Gassmann流体项反演方法[49]等,均需要考虑频率的影响才可对砂砾岩储层内部流体进行有效识别和区分。研究中尚未考虑砂砾岩储层的强非均质性,其导致的强各向异性对于地震波弹性响应的影响同样值得深入探究。此外,由于喷射流机制对应的特征频率出现在较高频段,而在103~106 Hz之间缺少相应的实验数据与预测结果进行比对,后续研究工作还需在高频段为其寻找更加有力的证据。
-
（1） 超声频带测量结果表明,区分砂砾岩与泥灰岩的主要敏感弹性参数（优选属性）为ρ,K,λ,λρ,λΦ,区分气—水的主要敏感弹性参数为v,ρVp/Vs,ρv,λv,无敏感参数可区分油—水。
（2） 地震频带测量结果表明,储层压力下砂砾岩纵、横波速度频散较泥、灰岩更强,反演所得裂隙密度更大。
（3） 联合超声和地震频带纵、横波速度测量结果,确定了砂砾岩中的重要频散物理机制为喷射流机制,预测其频散特征的裂隙纵横比为一组分布。基于喷射流理论,饱和油和水的岩石样品的频散所在的频带不同,因此可作为区分油—水的敏感参数。
Investigation on Elastic Parameters of Glutenite by Multi-Frequency-Band Measurements: A case study for rocks from the steep slope zone in Dongying Sag
-
摘要: 砂砾岩体不仅是重力流沉积的主要研究对象,在油气勘探和水文地质勘探领域同样占据重要地位。选用东营凹陷北部陡坡带典型砂砾岩样品,分析其物性,测量其超声频段弹性波速度,依据统计学理论计算指示因子,寻找岩性及流体敏感参数；基于超声频段弹性波速度反演裂隙密度和软孔隙特征,计算砂砾岩频散特征并与地震频段弹性参数测量结果对比,分析跨频段（地震频段和超声频段）岩石物理响应特征。结果表明：（1）依据敏感参数,砂砾岩与泥岩、砾岩差异明显,但砂砾岩与砂岩之间不易区分；对于流体,砂砾岩储层中气和水易区分,但油和水不易区分；（2）在原位储层压力下,相较于砂岩、砾岩、泥岩和灰岩,砂砾岩纵波速度频散范围更大,裂隙密度更大,相应的闭合压力范围也较大；（3）频散特征满足“喷射流”机制,其中微裂隙纵横比为一组分布,频散发生的频率范围及幅度可作为区分油和水的敏感参数。研究成果为区分砂砾岩储层中不同频段地震波响应特征及进一步识别流体类型提供了理论依据,有助于实现储层与非储层的地震精准预测。Abstract: Glutenites, which significantly contribute to gravity flow sedimentation studies, have become prime targets for oil, gas, and groundwater exploration, owing to their advantageous reservoir conditions. However, their complex depositional attributes, heterogeneous composition, and intricate pore structures pose challenges, affecting the elastic response of seismic waves. This complexity can compromise the reliability of fluid inversion methodologies and the outcomes of fluid identification processes. Consequently, a comprehensive understanding of the interrelationship between reservoir lithology, its physical properties, involved fluids, and seismic attributes is vital. Such understanding, attainable through petrophysical experiments, is critical for accurately deciphering the petrophysical response characteristics of glutenites. We undertook multi-frequency-band measurements （employing seismic frequencies at 1⁃100 Hz and ultrasonic frequency at 1 MHz） on the elastic parameters of the glutenite from the Dongying Sag in the Bohai Bay Basin. Initially, we selected a series of samples from each grain level for porosity and permeability tests, followed by classification and analysis of the rock types and physical properties. We then measured the P- and S-wave ultrasonic velocities in samples saturated with varying fluids, assessing the effects of fluid, lithology, and pressure on velocity. This led to the determination of the relationship between pressure and velocity by evaluating different water-saturated samples, followed by the calculation of indicative factors to identify rock types and fluid sensitivity parameters. Furthermore, we measured the Young’s modulus and Poisson’s ratio to calculate the P- and S-wave velocities of the seismic band. Lastly, we computed the crack density and soft porosity using ultrasonic velocity and compared the velocity dispersion with the seismic frequency band measurements. Our results suggest that： （1） sensitive parameters ρ, K, λ, λρ, and λϕ help distinguish glutenite from shale and conglomerate, but not sandstone. For fluids, gas-water can be separated through υ, ρVP/VS, ρυ, and λυ, whereas oil and water cannot； （2） at in-situ reservoir pressure, glutenite demonstrates larger P-wave velocity dispersion, crack density, and corresponding closure pressure compared to those of sandstone, conglomerate, and mudstone； （3） the dispersion can be explained by the "squirt flow" mechanism, aligning with a crack aspect ratio distribution that mirrors the actual scenario. The velocity dispersion range of oil-saturated samples is found at lower frequencies, with a larger dispersion span, thus positioning seismic band dispersion range as a sensitive parameter for differentiating between oil and water, overcoming the limitations of traditional sensitive parameters. This study offers insights into lithology and fluid-sensitive parameters via ultrasonic experiments and explores the pronounced velocity dispersion characteristics of glutenite under in-situ reservoir conditions through seismic frequency band experimental results, thereby identifying the key dispersion physical mechanism in glutenite. We assert that the fractured glutenite possesses excellent reservoir properties, exhibiting velocity dispersion characteristics markedly different from marl with fluid sensitivity. This prompts the consideration of frequency effects in fluid identification for traditional inversion methods. Our work systematically encapsulates the petrophysical response characteristics of glutenite bodies, establishing a theoretical groundwork for studying seismic wave response and identifying fluids in glutenite reservoirs. This could enhance reservoir prediction accuracy and effectively guide the geophysical exploration of glutenite reservoirs.
-
Key words:
- glutenite /
- reservoir sensitivity factor /
- low-frequency rock physics /
- dispersion /
- crack density
-
图 1 多频段岩石物理设备图
A.轴压加载；B.压电陶瓷震源；C.胡克腔；D.加温套；E.底座；F.应变片；G.超声探头；H.铝合金堵头；I.加压流体泵；J.惠斯通电桥；K.气瓶
Figure 1. Multi⁃frequency band apparatus diagram
A. piston for axial pressure loading; B. piezoelectric source (PIZO) for seismic source; C. tri⁃axial cell for confining pressure; D. heating sleeve; E. aluminum base; F. strain gauge; G, ultrasonic transducers; H. aluminum alloy plug head; I. pressurized fluid pump; J. Wheatstone bridge; K. gas cylinder
图 2 沙三段40块样品孔隙度统计结果
Figure 2. Porosity of 40 samples in the Third member of the Shahejie Formation (Es3)
Fig.2
图 3 沙三段40块样品渗透率统计结果
Figure 3. Permeability of 40 samples in in the Third member of the Shahejie Formation (Es3)
Fig.3
图 4 用于低频测试的代表性岩心样品
（a）16号,砂质砾岩；（b）8号,细砾岩；（c）10号,细砂岩；（d）13号,砾质砂岩
Figure 4. Representative core samples for low⁃frequency measurement
(a) No. 16, sandy anglomerate; (b) No. 8, fine⁃grained conglomerate; (c) No. 10, fine⁃grained sandstone; (d) No. 13, grave sandstone
图 5 饱和气、水和油的21块样品的纵横波速度随压力变化曲线
Figure 5. Evolution of P⁃and S⁃wave velocity versus increasing pressure for 21 samples of saturated gas, water, and oil
Fig.5
图 6 砂砾岩饱和水条件下的纵、横波速度频散
（a）样品16,砂质砾岩；（b）样品16,砂质砾岩；（c）样品8,细砾岩；（d）样品8,细砾岩；（e）样品10,细砂岩；（f）样品10,细砂岩；（g）样品13,砾质砂岩；（h）样品13,砾质砂岩
Figure 6. P⁃and S⁃wave velocity versus frequency under different effective pressures for glutenite saturated with water
Fig.6
图 7 岩性指示因子
（a）气饱和样品指示因子；（b）水饱和样品指示因子；（c）油饱和样品指示因子；（d）油饱和样品敏感参数交会图
Figure 7. Lithology indicator factors
(a) gas⁃saturated samples; (b) water⁃saturated samples; (c) oil⁃saturated samples; and (d) intersection diagram of sensitive parameters of oil⁃saturated samples
图 8 流体指示因子
（a）砂岩指示因子；（b）砾岩指示因子；（c）砂砾岩指示因子；（d）砂砾岩敏感参数交会图
Figure 8. Fluid indicator factors
(a) sandstone indicator factor; (b) conglomerate indicator factor; (c) conglomerate indicator factor; (d) intersection map of glutenite sensitive parameters
图 9 裂隙密度随压力变化情况
（a）16号样品,砂质砾岩；（b）8号样品,细砾岩；（c）10号样品,细砂岩；（d）13号样品,砾质砂岩
Figure 9. Crack density versus pressure
(a) sample No. 16, sandy conglomerate; (b) sample No. 8, fine⁃grained conglomerate; (c) sample No. 10, fine⁃grained sandstone; (d) sample No. 13, gravel sandstone
图 10 软孔隙随纵横比变化的分布特征
（a）16号样品,砂质砾岩；（b）8号样品,细砾岩；（c）10号样品,细砂岩；（d）13号样品,砾质砂岩
Figure 10. Soft⁃porosity distribution versus aspect ratio
(a) No.16, sandy conglomerate; (b) No.8, fine⁃grained conglomerate; (c) No.10, fine⁃grained sandstone; (d) No.13, grave sandstone
图 11 砂岩、砾岩、砂砾岩及泥灰岩综合分析
（a）纵波频散范围；（b）横波频散范围；（c）裂隙密度；（d）裂隙闭合压力
Figure 11. Comprehensive analysis of sandstone, conglomerate, glutenite, and mudstone⁃limestone
(a) P⁃wave dispersion; (b) S⁃wave dispersion; (c) crack density; (d) crack closure pressure
表 2 样品物性参数
岩心序号 密度/（g/cm3） 孔隙度/% 渗透率/×10-3 μm2 16 2.33 13.88 22.23 8 2.40 8.62 17.70 10 2.42 9.91 0.15 13 2.36 12.84 3.52 
下载: 导出CSV
表 3 16号、8号、10号和13号样品的矿物组分
岩心序号 黏土总量/% 石膏/% 石英/% 钾长石/% 斜长石/% 方解石/% 铁白云石/% 白云石/% 菱铁矿/% K0/GPa 16 6.3 3.3 27.4 0 34.2 15.4 3.7 9.7 0 62.2 8 5.0 3.8 26.6 0 29.4 23.8 3.8 7.7 0 63.2 10 4.7 0 25.8 4.4 16.7 33.3 4.7 7.9 2.5 66.2 13 5.0 3.8 26.6 0 29.4 23.8 3.8 7.7 0 63.2
下载: 导出CSV
-
[1] 刘强虎,朱红涛,杜晓峰,等. 渤海海域砂砾岩体沉积响应进展及热点[J]. 中国科学（D辑）：地球科学,2020,45（5）：1676-1705. Liu Qianghu, Zhu Hongtao, Du Xiaofeng, et al. Development and hotspots of sedimentary response of glutenite in the offshore Bohai Bay Basin[J]. Science China (Seri. D): Earth Science, 2020, 45(5): 1676-1705. [2] 王鑫,林承焰,马存飞,等. 东营凹陷北部陡坡带利563区块沙四上亚段砂砾岩扇体沉积特征及沉积模式[J]. 吉林大学学报（地球科学版）,2020,50（3）：705-720. Wang Xin, Lin Chengyan, Ma Cunfei, et al. Sedimentary characteristics and sedimentary model of glutenite fans in upper Es4 in L563 area, north steep slope of Dongying Depression[J]. Journal of Jilin University (Earth Science Edition), 2020, 50(3): 705-720. [3] 孔凡仙. 东营凹陷北带砂砾岩扇体勘探技术与实践[J]. 石油学报,2000,21（5）：27-31. Kong Fanxian. Exploration technique and practice of sandy-conglomeratic fans in the northern part of Dongying Depression[J]. Acta Petrolei Sinica, 2000, 21(5): 27-31. [4] 胡阳,刘惠民,郝雪峰. 断陷湖盆陡坡带砂砾岩油藏特征及控制因素：以东营凹陷古近系为例[J]. 地质论评,2019,65（增刊1）：151-152. Hu Yang, Liu Huimin, Hao Xuefeng. Characteristics and controlling factors of glutenite reservoir in steep slope zone of faulted lacustrine basin: A case study of Paleogene in Dongying Depression[J]. Geological Review, 2019, 65(Suppl. 1): 151-152. [5] Dodd T J H, Graham Leslie A, Gillespie M R, et al. Deep to shallow-marine sedimentology and impact of volcanism within the Middle Triassic palaeo-tethyan Semantan Basin, Singapore[J]. Journal of Asian Earth Sciences, 2020, 196: 104371. [6] Niu X B, Yang T, Cao Y C, et al. Characteristics and formation mechanisms of gravity-flow deposits in a lacustrine depression basin: Examples from the Late Triassic chang 7 oil member of the Yanchang Formation, Ordos Basin, central China[J]. Marine and Petroleum Geology, 2023, 148: 106048. [7] 朱筱敏,陈贺贺,谈明轩,等. 从太平洋到喜马拉雅的沉积学新航程：21届国际沉积学大会研究热点分析[J]. 沉积学报,2023,41（1）：126-149. Zhu Xiaomin, Chen Hehe, Tan Mingxuan, et al. A new journey in sedimentology from the Pacific to the Himalayas: Analysis of research hotpots from the 21st International Sedimentological Congress[J]. Acta Sedimentologica Sinica, 2023, 41(1): 126-149. [8] 解强旺,王艳忠,操应长,等. 东营凹陷陡坡带盐斜229地区沙四上亚段砂砾岩油藏成藏控制因素[J]. 中南大学学报（自然科学版）,2019,50（7）：1626-1636. Xie Qiangwang, Wang Yanzhong, Cao Yingchang, et al. Control factors on the hydrocarbon accumulation of the Es4s reservoirs in Yanxie 229 area, Dongying Sag[J]. Journal of Central South University (Science and Technology), 2019, 50(7): 1626-1636. [9] Kra K L, Qiu L W, Yang Y Q, et al. Sedimentological and diagenetic impacts on sublacustrine fan glutenites reservoir quality: An example of the Paleogene Shahejie Formation (Es4s member) in the Dongying Depression, Bohai Bay Basin (East China)[J]. Sedimentary Geology, 2022, 427: 106047. [10] Fawad N, Liu T X, Fan D D, et al. Sequence stratigraphic divisions and correlation of the middle sub-member of Eocene Shahejie Formation in the Bonan Sag of Bohai Bay Basin (China): Implication for facies and reservoir heterogeneities[J]. Geoenergy Science and Engineering, 2023, 225: 211622. [11] 徐长贵,于海波,王军,等. 渤海海域渤中19-6大型凝析气田形成条件与成藏特征[J]. 石油勘探与开发,2019,46（1）：25-38. Xu Changgui, Yu Haibo, Wang Jun, et al. Formation conditions and accumulation characteristics of Bozhong 19-6 large condensate gas field in offshore Bohai Bay Basin[J]. Petroleum Exploration and Development, 2019, 46(1): 25-38. [12] 李国欣,覃建华,鲜成钢,等. 致密砾岩油田高效开发理论认识、关键技术与实践：以准噶尔盆地玛湖油田为例[J]. 石油勘探与开发,2020,47（6）：1185-1197. Li Guoxin, Qin Jianhua, Xian Chenggang, et al. Theoretical understandings, key technologies and practices of tight conglomerate oilfield efficient development: A case study of the Mahu oilfield, Junggar Basin, NW China[J]. Petroleum Exploration and Development, 2020, 47(6): 1185-1197. [13] 邓继新,柴康伟,宋连腾,等. 差异性成岩过程对百口泉组砂砾岩岩石物理特征的影响[J]. 地球物理学报,2022,65（11）：4448-4459. Deng Jixin, Chai Kangwei, Song Lianteng, et al. The influence of diagenetic evolution on rock physical properties of glutenite of Baikouquan Formation[J]. Chinese Journal of Geophysics, 2022, 65(11): 4448-4459. [14] 高阳. 准噶尔盆地玛湖凹陷砂砾岩储层物性分类及控制因素[J]. 成都理工大学学报（自然科学版）,2022,49（5）：542-551,560. Gao Yang. Physical property classification and controlling factors of glutenite reservoir in Mahu Sag, Junggar Basin, China[J]. Journal of Chengdu University of Technology (Science & Technology Edition), 2022, 49(5): 542-551, 560. [15] 韩宏伟,崔红庄,林松辉,等. 东营凹陷北部陡坡带砂砾岩扇体地震地质特征[J]. 特种油气藏,2003,10（4）：28-30. Han Hongwei, Cui Hongzhuang, Lin Songhui, et al. Seismic geology of glutenite fan in the north actic region of Dongying Sag[J]. Special Oil and Gas Reservoirs, 2003, 10(4): 28-30. [16] Hori K, Nagasawa S, Sato Y, et al. Response of a coarse grained, fluvial to coastal depositional system to glacioeustatic sealevel fluctuation since the Last Glacial Maximum: An Example from the Tenryu River, Japan[J]. Journal of Sedimentary Research, 2017, 87(2): 133-151. [17] 周晓光,黄晓波,王启明,等. 渤海海域石南陡坡带多期砂砾岩扇体识别与展布特征[J]. 东北石油大学学报,2020,44（2）：46-55. Zhou Xiaoguang, Huang Xiaobo, Wang Qiming, et al. Identification and description of the multi-stage glutenite fan body in Shinan steep slope zone, Bohai Sea[J]. Journal of Northeast Petroleum University, 2020, 44(2): 46-55. [18] 胡鑫,邹红亮,胡正舟,等. 扇三角洲砂砾岩储层特征及主控因素：以准噶尔盆地东道海子凹陷东斜坡二叠系上乌尔禾组为例[J]. 东北石油大学学报,2021,45（6）：15-26. Hu Xin, Zou Hongliang, Hu Zhengzhou, et al. Reservoir characteristics and main controlling factors of glutenite reservoir in fan dalta glutenite: A case study of the Upper Urho Formation of Permian in the east slope of Dongdaohaizi Sag, Junggar Basin[J]. Journal of Northeast Petroleum University, 2021, 45(6): 15-26. [19] 雷蕾,韩宏伟,于景强. 近岸水下扇沉积样式及地震响应特征新认识[J]. 石油地球物理勘探,2019,54（5）：1151-1158. Lei Lei, Han Hongwei, Yu Jingqiang. New understandings of near-shore subaqueous fan sedimentary styles and its seismic responses[J]. Oil Geophysical Prospecting, 2019, 54(5): 1151-1158. [20] 杜猛,向勇,贾宁洪,等. 玛湖凹陷百口泉组致密砂砾岩储层孔隙结构特征[J]. 岩性油气藏,2021,33（5）：120-131. Du Meng, Xiang Yong, Jia Ninghong, et al. Pore structure characteristics of tight glutenite reservoirs of Baikouquan Formation in Mahu Sag[J]. Lithologic Reservoirs, 2021, 33(5): 120-131. [21] 吕复苏,黄小平,任涛. 地震属性信息在砂砾岩油藏开发中的应用：以克拉玛依油田上二叠统上乌尔禾组油藏为例[J]. 新疆石油地质,2003,24（4）：310-312. Fusu Lü, Huang Xiaoping, Ren Tao. Application of seismic attribute information in glutenite reservoir development-An example of Wuerhe reservoir of Upper Permian in Karamay oilfield[J]. Xinjiang Petroleum Geology, 2003, 24(4): 310-312. [22] 邢文军,吴开龙,吴鑫,等. 储层砂岩宽频段地震岩石物理特征的实验研究[J]. 地球物理学进展,2018,33（4）：1609-1616. Xing Wenjun, Wu Kailong, Wu Xin, et al. Multi-frequency laboratory measurement of rock physics property on sandstone[J]. Progress in Geophysics, 2018, 33(4): 1609-1616. [23] Müller T M, Gurevich B, Lebedev M. Seismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks: A review[J]. Geophysics, 2010, 75(5): 75A147-75A164. [24] Pride S R, Berryman J G, Harris J M. Seismic attenuation due to wave-induced flow[J]. Journal of Geophysical Research: Solid Earth, 2004, 109(B1): B01201. [25] Zhang L, Ba J, Carcione J M. Wave propagation in Infinituple-Porosity media[J]. Journal of Geophysical Research: Solid Earth, 2021, 126(4): e2020JB021266. [26] Zhang L, Ba J, Carcione J M, et al. Seismic wave propagation in partially saturated rocks with a fractal distribution of fluid-patch size[J]. Journal of Geophysical Research: Solid Earth, 2022, 127(2): e2021JB023809. [27] Gregory A R. Fluid saturation effects on dynamic elastic properties of sedimentary rocks[J]. Geophysics, 1976, 41(5): 895-921. [28] Domenico S N. Effect of brine‐gas mixture on velocity in an unconsolidated sand reservoir[J]. Geophysics, 1976, 41(5): 882-894. [29] Toksöz M N, Johnston D H, Timur A. Attenuation of seismic waves in dry and saturated rocks: I. Laboratory measurements[J]. Geophysics, 1979, 44(4): 681-690. [30] King M S, Marsden J R, Dennis J W. Biot dispersion for P- and S-wave velocities in partially and fully saturated sandstones[J]. Geophysical Prospecting, 2000, 48(6): 1075-1089. [31] Gassmann F. Über die elastizität poröser medien, Vier der Natur[J]. Gesellshaft in Zurich, 1951, 96: 1-23. [32] Lebedev M, Toms-Stewart J, Clennell B, et al. Direct laboratory observation of patchy saturation and its effects on ultrasonic velocities[J]. The Leading Edge, 2009, 28(1): 24-27. [33] Spencer J W. Stress relaxations at low frequencies in fluid-saturated rocks: Attenuation and modulus dispersion[J]. Journal of Geophysical Research: Solid Earth, 1981, 86(B3): 1803-1812. [34] Batzle M L, Han D H, Hofmann R. Fluid mobility and frequency-dependent seismic velocity:Direct measurements[J]. Geophysics, 2006, 71(1): N1-N9. [35] Mikhaltsevitch V, Lebedev M, Gurevich B. A laboratory study of elastic and anelastic properties of Savonnieres limestone[C]//Proceedings of the 76th EAGE Conference & Exhibition. Amsterdam, Netherlands: European Association of Geoscientists and Engineers, 2014: 1-5. [36] Subramaniyan S, Quintal B, Madonna C, et al. Laboratory-based seismic attenuation in Fontainebleau sandstone: Evidence of squirt flow[J]. Journal of Geophysical Research: Solid Earth, 2015, 120(11): 7526-7535. [37] Pimienta L, Fortin J, Guéguen Y. Experimental study of Young’s modulus dispersion and attenuation in fully saturated sandstones[J]. Geophysics, 2015, 80(5): L57-L72. [38] 未晛,王尚旭,赵建国,等. 含流体砂岩地震波频散实验研究[J]. 地球物理学报,2015,58（9）：3380-3388. Wei Xian, Wang Shangxu, Zhao Jianguo, et al. Laboratory study of velocity dispersion of the seismic wave in fluid-saturated sandstones[J]. Chinese Journal of Geophysics, 2015, 58(9): 3380-3388. [39] Sun C, Tang G Y, Fortin J, et al. Dispersion and attenuation of elastic wave velocities: Impact of microstructure heterogeneity and local measurements[J]. Journal of Geophysical Research: Solid Earth, 2020, 125(12): e2020JB020132. [40] Sun C, Fortin J, Borgomano J V M, et al. Influence of fluid distribution on seismic dispersion and attenuation in partially saturated limestone[J]. Journal of Geophysical Research: Solid Earth, 2022, 127(5): e2021JB023867. [41] Borgomano J V M, Gallagher A, Sun C, et al. An apparatus to measure elastic dispersion and attenuation using hydrostatic- and axial-stress oscillations under undrained conditions[J]. Review of Scientific Instruments, 2020, 91(3): 034502. [42] Mavko G, Mukerji T, Dvorkin J. The rock physics handbook: Tools for seismic analysis of porous media[M]. Cambridge: Cambridge University Press, 2009. [43] Avseth P, Jørstad A, van Wijngaarden A J, et al. Rock physics estimation of cement volume, sorting, and net-to-gross in North Sea sandstones[J]. The Leading Edge, 2009, 28(1): 98-108. [44] Dillon L, Schwedersky G, Vásquez G, et al. A multiscale DHI elastic attributes evaluation[J]. The Leading Edge, 2003, 22(10): 1024-1029. [45] 邓继新,周浩,王欢,等. 基于储层砂岩微观孔隙结构特征的弹性波频散响应分析[J]. 地球物理学报,2015,58（9）：3389-3400. Deng Jixin, Zhou Hao, Wang Huan, et al. The influence of pore structure in reservoir sandstone on dispersion properties of elastic waves[J]. Chinese Journal of Geophysics, 2015, 58(9): 3389-3400. [46] Duan C S, Deng J X, Li Y, et al. Effect of pore structure on the dispersion and attenuation of fluid-saturated tight sandstone[J]. Journal of Geophysics and Engineering, 2018, 15(2): 449-460. [47] Sun Y Y, Gurevich B. Modeling the effect of pressure on the moduli dispersion in fluid-saturated rocks[J]. Journal of Geophysical Research: Solid Earth, 2020, 125(8): e2019JB019297. [48] Gurevich B, Makarynska D, de Paula O B, et al. A simple model for squirt-flow dispersion and attenuation in fluid-saturated granular rocks[J]. Geophysics, 2010, 75(6): N109-N120. [49] 印兴耀,张世鑫,张峰. 针对深层流体识别的两项弹性阻抗反演与Russell流体因子直接估算方法研究[J]. 地球物理学报,2013,56（7）：2378-2390. Yin Xingyao, Zhang Shixin, Zhang Feng. 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. -
点击查看大图
计量
- 文章访问数: 110
- HTML全文浏览量: 21
- PDF下载量: 46
- 被引次数: 0
