天然气水合物储层吸收衰减机制及岩石物理理论研究进展

武存志 张 峰 李向阳

(中国石油大学(北京)，北京102249)

0 引言

反射地震是水合物地球物理勘探的重要技术。赋存天然气水合物的海底沉积物的反射地震响应往往存在似海底反射(BSR)特征[1]。BSR是由上覆高速的含水合物层与下伏低速的含游离气层之间的阻抗差异产生的，大体平行于海底，指示天然气水合物稳定带的底部，具有负极性、强振幅特征。但天然气水合物和游离气与BSR特征的关系仍欠明确，基于BSR的天然气水合物识别具有多解性[2]。此外，BSR仅能指明天然气水合物底部位置，却无法用于水合物含量的定量解释。

地震波的吸收衰减是天然气水合物地层地震响应的一个显著特征。研究发现，在加拿大Mallik冻土带、日本南海海槽等水合物探区BSR位置的上方存在与天然气水合物有关的低地震振幅或振幅“空白”带[3]。虽然空白的成因尚不完全清楚，但人们普遍认为空白的“强度”与地震波吸收衰减直接相关，可作为天然气水合物分布的一项指标。多个水合物探区的声波测井波形也显示，天然气水合物的存在导致了异常的衰减值[4]。国内外学者针对性地开展了岩石物理实验和理论研究，目的是探索地震波衰减机制及其与天然气水合物含量的定量关系，提高水合物储量的预测能力。

本文概述了不同探区天然气水合物储层的地震波吸收衰减特征，探讨并总结了水合物储层可能存在的衰减机制及相关的岩石物理理论模型研究现状，展望了水合物储层吸收衰减特性的应用及岩石物理模型构建的未来发展方向。

1 不同地区天然气水合物储层的地震波吸收衰减特征

1.1 加拿大Mallik冻土带

Mallik冻土带位于加拿大西北部的Mackenzie三角洲，是世界上已经证实的冻土带中含有天然气水合物的典型地区。该区2L38井的阵列声波测井波形记录结果显示，含水合物地层的波形振幅发生明显减弱[4](图1)，表明声波在含水合物地层中传播时可能发生了强烈吸收衰减。Guerin等[4]利用质心频率偏移法分别计算了纵波和横波衰减，发现在含水合物地层二者都表现为高衰减值。文中用品质因子Q的倒数(逆品质因子Q-1)表示衰减强度，则高衰减对应高Q-1值或低Q值。进而对含水合物饱和度和纵、横波衰减(QP-1、QS-1)进行交会分析，得知纵、横波衰减与水合物饱和度大致呈线性关系，且相同含水合物饱和度条件下，横波衰减大于纵波衰减(图2)。这也是在含水合物地层吸收衰减特性分析研究中，首次将纵、横波衰减特性与含水合物饱和度进行定量分析。

图2 Mallik-2L38井中利用电阻率法估算的天然气水合物饱和度与声波衰减之间的关系[4]

Pratt等[5]和Bauer等[6]利用基于射线的层析成像技术对井间地震资料进行了衰减剖面的成像。结果显示，在含水合物地层地震波衰减也表现出明显的高值异常。由于声波测井频率一般在100～25kHz，而井间地震数据频率是在100～1000Hz(如采用美国Tomoseis公司的压电震源、10单元有源阵列接收系统，扫描频率为100～2000Hz，处理后频率为100～1000Hz)，因此Pratt等[5]和Bauer等[6]的研究是对Guerin等[4]分析结果的重要补充。

1.2 日本南海海槽

日本南海海槽是上新世以来菲律宾海板块向欧亚板块俯冲形成的年轻海沟，在水深800～2000m的内陆坡，分布数个新近系沉积充填的弧前沉积盆地。在弧前盆地的外侧，较深的陆坡处发育增生楔[7]。

由于该区在早期地震调查中识别出广泛分布的BSR，为了进一步确认天然气水合物藏的存在，JMITI(Japanese Ministry of International Trade and Industry)在日本南海海槽展开钻探，采集到含有大量

天然气水合物的砂岩岩心。与加拿大Mallik冻土带相似，在该区的阵列声波测井中，声波波形振幅也发生明显减弱。Matsushima[8]使用质心频率偏移法对阵列声波测井波形记录数据分别计算了纵波和横波衰减，虽然该区含水合物地层较薄，且水合物饱和度较低，但纵、横波衰减仍然都表现为高值异常，并且纵、横波衰减与含水合物饱和度也大致呈线性关系。Matsushima[9]又利用VSP数据进行了吸收衰减参数的提取，在地震频带范围内进一步研究了日本南海海槽含水合物地层的地震波吸收衰减特性。结果显示在含游离气区域，纵波衰减表现出明显高值异常，而在含水合物地层位置，VSP提取的衰减参数并没有明显的异常(图3)。因此Matsushima认为，含水合物地层的吸收具有频率选择性，在30～110Hz的频率范围内，含水合物地层并未对地震波产生强烈的吸收作用。

图3 对日本南海海槽PSW1现场VSP数据利用谱比法(a)和质心频率偏移法(b)估算的纵波逆品质因子[9]

1.3 阿拉伯海Makran增生楔

Makran增生楔位于西印度洋阿拉伯海海域，该区虽无井资料直接证明天然气水合物的存在，但Sain等[10]、Ojha等[11]基于该区地震资料通过多种方法证明了天然气水合物和游离气藏的存在。为了分析该区含水合物地层的地震波吸收特征，Sain等[12]在叠前道集上利用谱比法针对三个典型CDP处的水合物上覆地层、水合物稳定带及游离气层分别求取层间品质因子(图4)，这三个CDP分别对应强BSR、弱BSR及无BSR三个典型位置。分析这三个CDP处的地层品质因子变化规律可知，Makran增生楔地区含水合物地层吸收效应较弱，而BSR下部游离气藏的吸收效应较强。Sain等[12]指出，在Makran增生楔地区，含水合物地层的地震吸收特征与Mallik地区和日本南海海槽区的特征不一致，可能缘于该区含水合物地层中只有水合物和岩石骨架两相，而不存在流体，水合物与骨架的胶结作用增加了岩石的硬度从而造成地层吸收效应减弱。

图4 Makran增生楔地区 CDP 4286(a)、4372(b)和4524(c)处的地震道集分别与估计的层间品质因子Q叠加[12]

1.4 墨西哥湾

墨西哥湾西北陆缘是油气十分富集的被动大陆边缘，中新生代沉积盆地发育了巨厚海相碳酸盐岩，形成大规模的盐底辟构造和断裂构造。目前在众多区块识别出BSR[13-15]，并采集到大量天然气水合物样品。在该区阵列声波测井中，声波振幅也发生明显减弱。王吉亮等[16]使用质心频率偏移法对墨西哥湾GC955H井的阵列声波测井波形数据计算了纵波衰减(图5)。该区含有两种不同类型的含水合物地层，裂隙填充型的水合物层Ⅰ及孔隙填充型的砂岩水合物层Ⅱ。前者衰减大小与上覆背景泥岩层段相比并无明显变化，推测原因有三：①泥岩层水合物对衰减无影响；②泥岩层段内水合物饱和度太低；③可能是因为水合物充填于垂直裂隙中，对进入地层垂向传播的声波作用不大。而后者表现出明显的衰减高值异常，并且纵波衰减与含水合物饱和度也大致呈线性关系，随着水合物饱和度增加，纵波衰减逐渐增强。

图5 墨西哥湾GC955H井声波衰减计算结果[16]

1.5 中国珠江口盆地南海神狐海域

位于南海北部陆缘的珠江口盆地珠二坳陷白云凹陷南侧的神狐海域是中国近年来天然气水合物勘探的重点目标。广州海洋地质调查局于2007年及2015年先后在神狐海域进行了GMGS1、GMGS3两航次水合物钻探，均获得了水合物实物样品。钻探资料研究结果表明，神狐海域水合物储集层以黏土质粉砂细粒沉积物为主，局部存在具有较高孔渗性的粉砂细粒沉积物，水合物以分散式或者胶结式充填在其孔隙中[17]。GMGS1钻探区的水合物分布在海底以下约200m，厚度处于10～47m，平均饱和度约20%，最高达48%。GMGS3钻探区的水合物主要分布于海底以下100～250m，饱和度介于0～76%，较GMGS1钻探区厚度大、饱和度高、展布规模及储量大[18]。目前，针对中国南海水合物地震吸收衰减特征的研究较少。Li等[19]对神狐海域天然气水合物的反射地震数据分析表明，含水合物层的吸收衰减明显低于相邻地层。Wang等[20]利用神狐海域垂直电缆地震(VCS)资料，采用谱比法估算了海底各层品质因子Q值，研究发现，含水合物地层的Q值比不含水合物或游离气地层的Q值大，即吸收衰减弱，含游离气地层的Q值比不含水合物或游离气地层的Q值小，即吸收衰减强，但造成这一特征的主要原因和影响因素有待进一步研究。

1.6 人造岩样实验室测量结果

Priest等[21]在实验室合成水合物饱和岩石样品，利用共振法在频率50～400Hz下对13个含不同体积水合物的砂岩样本进行了衰减值的测量。结果(图6)显示：在水合物饱和度3%～5%范围内，纵横波衰减随饱和度增加迅速增至最大值，尔后纵波衰减迅速降低，接着又呈现缓慢增加的趋势，而横波衰减迅速减小且逐渐保持平稳，但仍比不含水合物地层的衰减值高2～3倍；当水合物饱和度增至3%～5%的临界饱和度时，衰减达到最强。

图6 Priest实验试样H4-1及H5-1在500kPa围压下纵波衰减及横波衰减随水合物饱和度变化[21]

Priest等[21]认为，实验室合成的含水合物砂岩样品为部分饱和，包含岩石骨架、水合物、吸附水和游离气等多个相态。在岩石中的水合物形成之前，地震波传播会引起吸附水和游离气在岩石裂隙和孔隙中产生喷射流，从而造成地震波能量的衰减。当水合物形成之后，水合物与岩石骨架的胶结会增加裂缝和孔隙的长度，从而加剧衰减作用。随着水合物饱和度进一步增加，水合物与岩石的胶结会逐渐完全堵塞裂隙和孔隙，吸附水和游离气被束缚，喷射流效应减弱，因而吸收作用逐渐减弱，但由于水合物本身也是多孔隙结构，吸附水和游离气可能在水合物内部结构中也会产生一定的喷射流效应，所以此时含水合物地层的吸收衰减效应仍强于不含水合物地层。对于横波的衰减只在水合物饱和度小于3%～5%时才会有敏感变化，该结论与Mallik地区的测井分析结果并不一致。在Mallik地区，随着水合物饱和度增加，横波衰减要比纵波衰减严重。Priest等[21]认为这主要是因为Mallik地区为完全饱和岩石，只含有岩石骨架、水合物和流体水三种相态，而实验样本含气体，气体的膨胀压缩作用，导致纵波衰减的变化更明显。

1.7 小结

为便于对比，将五个水合物探区及人工岩样的纵、横波衰减特征进行归纳总结(表1)。可以发现，对于不同的水合物探区或样品、使用不同的资料，地震波表现出不同的衰减特征，目前人们对此现象并未达成共识，这也是后续开展水合物地层吸收衰减机制及岩石物理研究的主要论题之一。

表1 各水合物探区及人工岩样的地震波衰减特征对比

2 天然气水合物储层吸收衰减机制及岩石物理理论研究

地下油气储层为多孔隙介质，孔隙结构的不规则性和孔隙介质的多相性导致多孔介质的宏观物理性质对孔隙流体含量、流体性质或孔隙结构的细小变化十分敏感。因此，当地震波在多孔隙介质中传播时，受其影响，介质孔隙发生闭合或张开、流体产生相对运动，致使多孔隙岩石的宏观物理性质甚至化学性质发生变化，从而引起地震波传播速度的频散、能量和振幅的频率依赖性衰减。可见，研究地震波在多孔隙介质中传播特性随孔隙流体及孔隙结构动态变化的关系，认识地震波在其中的传播规律，对于解决譬如储层流体预测等问题有重要的理论与实际意义。

国内外学者对含流体孔隙介质中的地震波传播与衰减问题进行了大量研究，建立了多种双相孔隙介质波动理论模型，如经典的Biot理论[22]、White模型[23]、BISQ理论[24]等。这些理论模型广泛应用于常规的油气储层研究，并取得了显著的应用效果[25]。而天然气水合物因其自身物理状态(冰态、结晶状化合物)及赋存方式(常见孔隙填充型、颗粒支撑性、胶结型)的特殊性，对储层流变性质及孔隙结构的影响远比油气层复杂得多，基于传统的理论来模拟天然气水合物与其宿主沉积物之间的相互作用、了解能量耗散过程，结果并不准确。因此学者们基于以上理论，考虑了各种水合物存在时可能带来的附加衰减机制，用来描述含水合物储层的地震波吸收衰减现象，由此形成了一系列岩石物理模型，如Leclaire模型[26]、CT模型[27]、改进的Leclaire模型[28]、HEG模型[29]、HBES模型[30]等，并结合水合物探区实际数据取得了较好的理论模拟及应用效果，为利用地震波衰减属性进行水合物储层识别及定量评价奠定了理论基础。

2.1 Leclaire模型(三相Biot理论)

Biot[22]最早(1956年)提出的描述含饱和流体孔隙介质的弹性波理论——Biot双相介质理论，奠定了孔隙介质弹性波理论的基础。该理论假设介质为各向同性、孔隙空间连续，且孔隙流体是相对于固体骨架可以进行Poiseuille流动的牛顿流体，并将孔隙介质受到的力分解为作用于固体骨架和孔隙流体两部分的应力，当地震波传播经过孔隙介质时，导致了固体骨架和孔隙流体之间的相互运动，其能量由于黏性摩擦而发生衰减，这就是Biot耗散机制。Biot理论不仅首次提出了流固耦合作用及其对地震波响应的影响，还预测了慢纵波的存在。

Leclaire等[26]认为在0℃以下，孔隙空间中的水和冰晶共存，并将结冰过程中的冰晶从连续态转化为不连续态，将适用于饱和流体多孔介质的Biot理论拓展到包含岩石颗粒、冰晶和孔隙流体三相的多孔介质，称其为Leclaire模型或三相Biot理论(图7)。该模型有两个重要假设，一是岩石颗粒与冰晶之间无直接接触；二是最小波长大于孔径，使骨架在宏观上可看成有效介质。模型基于连续介质力学给出了介质的动能、势能和耗散能量，推导出在岩石颗粒、冰晶、孔隙流体三相介质中波的传播方程

图7 含岩石颗粒、冰晶和孔隙流体的Leclaire模型[26]

(1)

式中：φ和ψ表示波的纵向和横向位移，上标点表示对时间求导；R、μ、ρ和F分别表示刚度、剪切模量、质量密度和摩擦系数矩阵，且有

(2)

Guerin等[4]将Leclaire模型中的冰晶替换为水合物，分析了不同水合物饱和度下储层衰减变化(图8)。模型孔隙度设置为35%，纵波频率采用12kHz，横波频率为2.5kHz，图中点线、虚线、实线分别对应地层渗透率κ为10-8m2、10-9m2、10-10m2，黑点为加拿大Mallik冻土区2L38含水合物井的实测数据。分析发现虽然模型与实测数据κ=10-10m2的条件下能匹配，但模型曲线并未很好地描述实际数据的趋势，其主要原因是该模型未考虑岩石颗粒与冰晶之间的相互作用。

图8 纵(左)、横(右)波逆品质因子随水合物饱和度的变化[4]

2.2 CT模型

Carcione等[27]考虑了岩石颗粒与冰晶之间接触而产生的势能与动能，发展了Leclaire理论，即CT(Carcione Tinivella)模型，并将CT模型中的冰晶替换为水合物颗粒，应用于计算水合物沉积地层的弹性波速度。但CT模型忽略了岩石颗粒与水合物之间的相互摩擦。Guerin等[28]利用Leclaire模型和CT模型分析了水合物饱和度及岩性对含水合物沉积层纵、横波速度及衰减的影响(图9)。相比Leclaire模型，CT模型中纵、横波速度随水合物饱和度增大的趋势更明显。此外，横波速度曲线表明：当水合物饱和度小于30%时，横波速度的增量可忽略；而当水合物饱和度较高时，横波速度的增加非常显著。但布莱克海脊数据表明，横波速度在低水合物饱和度时也有明显变化，与模型不相符[31]。Leclaire模型和CT模型均考虑了黏性流体摩擦，但衰减特征并不显著，这与Mallik冻土区(图8)数据明显矛盾，其原因是Guerin等[4]使用水合物渗透率和剪切模量的极值。表明Leclaire模型和CT模型均不能准确描述该含水合物地层的衰减现象。

2.3 改进的Leclaire模型

2.3.1 考虑沉积物颗粒与水合物间胶结作用

Guerin等[28]假设天然气水合物和沉积物颗粒之间的胶结作用可用耦合剪切模量来描述。使用类似于Johnson等[32]固结模型中剪切模量的表达式，该模量的简单表达式为岩石颗粒和天然气水合物体积分数的二阶多项式函数

μsh=μsh0(φhφs)2=μsh0[φh(1-φ)]2

(3)

图9 Leclaire模型和CT模型的模拟结果[28]

2.3.2 考虑沉积物颗粒与水合物之间的摩擦作用

岩石矿物颗粒与天然气水合物之间的胶结作用是速度随水合物饱和度增加而增加的必要条件[32]，这两相之间可能发生接触和摩擦，而Leclaire模型和CT模型都忽略了摩擦。Berryman等[33]描述了双孔隙、双渗透率介质中固相之间的摩擦。Guerin等[28]假设固体颗粒与水合物之间的摩擦力是其间速度差的线性函数，类似于对孔隙流体的黏性摩擦力，并定义新的摩擦系数为天然气水合物和沉积物组分的二阶多项式函数，类似于Berryman等[33]计算的阻力系数，从而给出摩擦矩阵，将沉积物颗粒与水合物之间的摩擦损耗引入Leclaire模型：

(4)

2.3.3 考虑沉积物颗粒、水合物、孔隙水之间的惯性耦合

2.3.4 考虑喷射流

图10 各Leclaire模型模拟结果(线)与Mallik-2L38井数据(点)对比[28]

将以上模型结果与Mallik-2L38井实际数据对比发现，在三相Leclaire模型中，摩擦、胶结和喷射流动机制的组合(无惯性耦合)，能够较合理地再现Mallik探区速度和衰减随水合物饱和度增加而增加的情况。

2.4 HEG(Hydrate Effective Grains)模型

2.4.1 BSM模型与存在问题

Stoll等[36]假设松散海洋沉积物骨架中存在全局摩擦损耗。在Biot理论中引入常数复骨架剪切模量和体积模量，这种对Biot理论的扩展，称为Biot-Stoll模型(BSM)。Stoll[37]的实验室数据表明，在松散的细粒海相沉积物中，纵波与横波的对数衰减率非常接近。因此，如果给定了复剪切模量，则可以通过假设骨架的合理泊松比估计无法通过测量直接获得的复骨架体积模量；Stoll[38]实验室数据进一步解释了此结果，认为细粒海洋沉积物中，存在由于骨架中的局部黏弹性损失引起的附加频率依赖性衰减，即全局骨架摩擦损耗，还使用Cole-Cole模型[39]中的公式引入复剪切模量。通过与实际数据对比，发现Biot理论明显低估了衰减，而Biot-Stoll模型也给出了低估的衰减值。Wood等[40]认为为了匹配这些数据而假设一个更高的剪切模量是不合理的，因为会导致很高的纵波速度。

2.4.2 EGM模型

Leurer[41]认为需一个附加衰减机制来增加纵波能量的衰减，该机制在加载期间分别通过纯压缩和膨胀激活，即使剪切波对其吸收应非常高，也对纵波速度和衰减的影响很小，因为细颗粒近地表海洋沉积物的剪切模量很低[42]。因此Leurer[41]提出了等效颗粒模型(Effective-Grain Model，EGM模型)，认为通常假设的体积模量对沉积物颗粒的描述是不充分的，因而假设沉积物具有较多的黏土组分，并且存在膨胀性的黏土矿物(以蒙脱石为代表)，因为在几乎所有含黏土的海洋沉积物中都有不同数量的膨胀性黏土矿物。由于膨胀，即吸水，蒙脱石变成一个两相体系。为了描述等效颗粒材料的动态行为，膨胀蒙脱石被认为是均匀各向同性的平面柱状包裹体，假设包裹体被晶内水饱和，包裹体边缘与自由孔隙水所占据的孔隙空间相连。因此，颗粒材料被视为对应力具有黏弹性响应的等效介质，引入复体积模量来描述。衰减机制为应力诱导挤压引起的流体运动，因此可视为喷射流动过程，这一现象的一般特征是从低纵横比包裹体到纵横比接近统一的孔隙的局部流体流动。

2.4.3 EGM的拓展——含水合物层HEG模型

Priest等[43-44]使用过量气(Excess gas)和过量水(Excess water)两种方法在砂岩样本中形成甲烷水合物。过量气法包括将已知质量的水分布在整个砂样中，用甲烷气体使部分水饱和的砂饱和，水倾向于覆盖湿砂粒，甲烷与水反应形成颗粒覆盖型水合物，其中一部分作用于在颗粒接触处胶结颗粒[45]。相比之下，过量水法需要向砂样中注入甲烷气体，然后用水冲洗砂样，这被认为产生了一种孔隙填充型水合物形态，因为甲烷气体倾向于分布在砂岩的孔隙空间中，周围是自然附着在湿砂颗粒上的水。该解释分别得到了文献[44]和文献[46]中报告的速度和弹性模量结果的支持。Best等[29]提出水合物颗粒可认为是含有气体或水的柔性材料，而不是单纯的固体矿物。包裹体是在Excess gas或Excess gas条件下水合物形成过程中捕获的孤立气或水的结果。在此情况下，水合物可能以与沉积物和岩石中其他微孔、柔性矿物类似的方式，产生局部黏性流体(喷射流)。Ecker等[47-48]不仅给出关于颗粒胶结水合物的骨架模量表达式，还给出承载型和孔隙填充型水合物的骨架体积和剪切模量的表达式，并将这些骨架模量值输入到Gassman[49]方程，以计算流体饱和沉积物中的地震速度。Best等[29]将其引入Biot理论中，以预测频率相关的纵波和横波速度和衰减。为了解释与微孔水合物颗粒相关的喷射流，Best等[29]使用Leurer等[41，50]的EGM方法，使骨架弹性模量是与频率相关的复模量。这里，Ecker等[47-48]表达式中的固体水合物相的体积模量和剪切模量被含有流体包裹体的等效水合物颗粒的体积模量和剪切模量代替。并遵循Johnston等[51]的方案，利用对应性原理引入频率相关的复包裹材料弹性模量，使水合物的等效颗粒模量为频率相关的复模量。由此提出水合物等效颗粒模型(Hydrate Effective Grains，HEG，图11)。

图11 水合物等效颗粒模型示意图[47]

Best等[29]模拟发现，胶结水合物HEG模型估算的衰减可忽略不计；孔隙填充水合物HEG模型只能解释饱和水条件下纵波衰减的增加(预测横波衰减为零)；承载HEG模型预测纵波和横波衰减的增加。

图12 气体饱和时承载型HEG模型模拟结果[29]

图13显示含水饱和(Sw=1)、承载型水合物的类似结果，图13c和图13d显示通过固定Si=0.5和不同的包裹体纵横比预测Sh=0.15时衰减峰值。

图13 水饱和时承载型HEG模型模拟结果(同图12)

图14 承载型HEG模型结果与Best等[29]实验数据对比

2.5 HBES模型

Marín-Moreno等[30]给出含水合物沉积物微观结构的理想状态(图15)。他们将HEG模型[29]引入的水合物中流体包裹体与孔隙之间的流动称为亚微米喷射流，并以HEG模型和Biot-Stoll模型[36]为基础，将水合物生成孔隙和原始孔隙(如砂岩骨架主孔隙)之间的微米喷射流以及气泡阻尼引起的衰减加入其中，提出含水合物等效沉积物模型(Hydrate-Bearing Effective Sediment，HBES)，模拟多重衰减机制下纵波衰减随水合物饱和度的变化。图16模型示意图展示了不同的单独衰减机制，含水合物沉积物中的总能量损失由每种衰减机制的贡献给出。

图15 含水合物沉积物微观结构的理想状态[30]

图16 含水合物有效沉积物(HBES)模型流程图[30]

图17为孔隙填充型水合物层的的纵波衰减随频率及饱和度的变化，地震频率下显示的衰减峰值是由水合物(亚微米喷射流)中的流体包裹体引起的，而在超声波频率下，衰减峰值是由孔隙填充水合物存在导致的初始孔隙度降低和不同纵横比孔隙的引起的。水合物中流体包裹体产生的亚微喷射流的衰减最大值出现在整个频率范围内，这取决于包裹体中流体的纵横比和类型(图18a和图18b)，而不同纵横比孔隙产生的微喷射流的峰值出现在声波频率和超声波频率下，与孔隙的纵横比无关(图18c)。本质上，水合物中流体包裹体的类型控制着衰减变化的幅度和衰减峰值对应频率。相比之下，水合物中流体包裹体的纵横比仅控制衰减峰值的频率依赖性，而包裹体的浓度仅控制衰减幅度。

图17 孔隙填充型HBES模型结果[30]

图18 孔隙填充型HBES模型变化模拟结果(相对于不含水合物)[53]

将HBES模型结果(图19)与Priest等[21]实验数据对比，可发现模型预测的纵横波速度与实验数据吻合较好，也同样预测了较低饱和度下的纵波衰减峰值，但结果偏低，而预测的横波衰减为零，与数据不符。

图19 孔隙填充型HBES模型结果(相对于不含水合物)与Priest等[21]实验数据对比[43]

2.6 小结

为方便读者对比，本文将各种岩石物理模型的假设条件与优缺点、衰减机制、模拟结果及对实际探区或实验数据的适用性进行了归纳总结(表2)。

表2 各模型假设条件与优缺点、衰减机制、模拟结果及对实际探区或实验数据适用性分析对比

3 结论与展望

本文回顾了国内外多探区天然气水合物地震波吸收衰减特征，综述了水合物储层可能存在的衰减机制及相关岩石物理理论，得到以下几点认识：

(1)很多地区的水合物地层的地震数据呈现较明显的吸收衰减特征，但是相互之间存在较大的差异，造成差异的主要原因是水合物的赋存状态不同，以及储层是否含流体。

(2)地震衰减具频散特征，但现有的现场观测和实验测试数据有限，难以全面体现衰减随频率变化的性质。

(3)水合物对衰减的作用还不明确，现有的吸收衰减机制和岩石物理理论在一定程度上能够描述地震波衰减随水合物饱和度的关系，例如改进的Leclaire模型与井数据相匹配，HEG模型与Best的实验数据的趋势相一致，但是适用性和精度还需要进一步结合实际数据进行分析和验证。

(4)由于实验室内较难合成与原位地层条件相同的水合物样品，基于人工合成岩心的岩石物理实验得到的水合物与各种岩石物理参数之间的关系，并不能准确描述实际地层的情况。

(5)与天然气水合物存在有关的测井声波波形异常及低地震振幅或振幅“空白”带，与地震波吸收衰减直接相关，可作为天然气水合物量化的一个指标。但由于我们对沉积物中的衰减机制，特别是水合物的作用认识有限，因此对声波测井或者地震衰减测量的解释是不确定的，很难定量地将吸收衰减与特定的水合物或沉积物条件联系起来，但实验室研究可为其提供理论基础。

基于以上认识，本文认为：

(1)基于孔隙介质地震波传播理论，开展一系列赋存方式和分布特征各异的含水合物储层的吸收衰减机制、岩石物理理论和衰减特征研究，对实验和现场数据的更可靠解释具有决定性的作用，可提高将地震振幅与天然气水合物含量联系起来的可能性，提高估算水合物储量的能力。

(2)对于含水合物沉积物，研究水合物对岩石微观结构改变的作用机理，分离固有损耗机制并准确预测其依赖性(如测量频率、等效压力、温度、渗透率、孔隙流体类型和饱和度、水合物饱和度和形态)，从而建立描述这种衰减机理的最佳岩石物理模型，仍是目前的主要挑战。