0 引言
天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内。地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] 。
速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多。实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层。相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式。沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] 。因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] 。水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂。孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] 。
正是由于影响水合物速度的因素多,单纯利用速度识别水合物可能会带来一定误差[19 ] ,特别是在海底浅层的MTD、河道沉积往往产生与水合物类似的反射特征,此时则需要去除由于孔隙度变化引起的振幅异常,岩石物理建模分析可以有效地指导开展高精度的水合物识别,其关键在于构建高敏感度的水合物指示因子。事实上,在实际生产应用中,针对不同储层应使用不同的弹性参数或者构建相应的指示因子[20 ] ,叠后阻抗、纵横波速度、泊松比、Gassman流体因子、杨氏模量、泊松比、流体模量、组合流体因子等相继被提出并应用到相关研究中,且国内外众多学者相继提出了不同弹性参数的反演方法,但是对于水合物储层的流体因子构建和反演研究尚少。
叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似。为了提高反演的稳定性,本文使用弹性阻抗反演的方法。弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述。但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法。
笔者基于SCA-DEM岩石物理建模方法[30 ] ,针对水合物未固结特性合理设置建模参数,定量地描述物性参数对弹性参数以及AVO 特征的影响,并在此基础上构建了水合物指示因子。相比于传统弹性参数,该指示因子与水合物饱和度具有更高相关性。然后在平面波近似下,推导了与该指示因子直接相关的AVO反射特征方程,基于该方程则可直接获取水合物指示因子,提高水合物识别精度。
1 速度影响因素分析
1.1 建模参数优选
首先,根据矿物组分含量,选取VRH公式估算矿物(水合物、方解石、白云石和泥质)混合后的模量,VRH计算公式为:
(1) M VRH = M V + M R 2
M V = ∑ i = 1 N f i M i , 1 M R = ∑ i = 1 N f i M i ,
式中:M VRH 、MV 和MR 可以是任何模量,如体积模量、剪切模量或杨氏模量等;fi 是第i 个组分的体积分量;Mi 是第i 个组分的弹性模量。
然后,利用SCA-DEM模型向基质中加入有一定连通性的干燥微孔隙,假设孔隙为固定纵横比的理想椭球体,计算干岩石骨架的弹性模量,由于水合物储层未固结,因此给定的纵横比小于 0.01,表现为颗粒接触边界模型。各向异性的SCA模型如下所示:
(2) C SCA = ∑ n = 1 N v n C n [ I + G ˙ ( C n - C SCA ) ] - 1 × ∑ n = 1 N v n [ I + G ˙ ( C n - C SCA ) ] - 1 - 1
式中: G ˙ I 是单位张量;Cn 是第n 种成分的刚度张量;vn 是第n 种矿物成分的体积含量。最后向孔隙中加入流体等,则可以获得饱和岩石的刚度张量C sat 。纵波速度、横波速度则可以表达为:
(3) v p = C 33 sat ρ v s = C 55 sat ρ
特别地,在SCA-DEM模型中,我们需要首先确认模型中的建模常数,包括矿物组分的模量、孔隙纵横比、孔隙内流体类型等,其中,孔隙纵横比较难确定。海底浅层物性参数变化较小,假设孔隙结构变化也应较小,水合物富集在浅层未固结岩层中,软孔隙较发育,纵横比小于0.01。图1 展示了南海琼东南海域一口含水合物井的部分测井曲线,水合物层(1 800~1 860 m)为高电阻率、高纵波速度特征。通过对比不同孔隙纵横比建模结果与实测数据,可以获取研究区的孔隙纵横比,如图2 所示,当硬孔隙含量较低(平均值约0.3)时,模型计算出的纵波速度与实测的纵波速度几乎完全一致,表明此时模型参数设置合理,可以为下一步分析提供基础。
图1
图1
南海琼东南海域一口含水合物井的部分测井曲线
Fig.1
Log data of a typical hydrate-bearing well in Qiongdongnan area of South China Sea
图2
图2
建模结果与对应的硬孔隙度
Fig.2
Modeling results with corresponding hard porosity
1.2 速度的影响因素
我们的目的是研究速度对哪些因素最敏感,固定孔隙纵横比后,改变孔隙度、饱和度等参数则可以研究不同物性参数与速度的关系。
图3 为纵横波速度随着水合物饱和度和含气饱和度的变化趋势,图4 为纵横波速度随着孔隙度和含气饱和度的变化趋势。不难发现,含气饱和度对纵波速度影响显著,即使少量含气也会使得纵波速度明显减小。孔隙度与速度负相关,水合物饱和度与速度正相关(图3 a,图4 a),图5 展示了不同含气饱和度下速度随着孔隙度理论上的变化,即纵波速度与孔隙度之间呈指数型关系,横波速度与孔隙度之间几乎呈线性关系,水合物饱和度与孔隙度对速度影响都很大。
图3
图3
孔隙度为45%时纵横波速度与水合物饱和度、含气饱和度的关系
a—纵波速度与水合物饱和度、含气饱和度的关系;b—横波速度与水合物饱和度、含气饱和度的关系
Fig.3
Velocity corresponding to different gas saturation and hydrate saturation when porosity is 45%
a—P-velocity corresponding to different gas saturation and hydrate saturation;b—S-velocity corresponding to different gas saturation and hydrate saturation
图4
图4
纵横波速度与孔隙度、含气饱和度的关系
a—纵波速度与孔隙度、含气饱和度的关系;b—横波速度与孔隙度、含气饱和度的关系
Fig.4
Velocity corresponding to different gas saturation and porosity
a—P-velocity corresponding to different gas saturation and porosity;b—S-velocity corresponding to different gas saturation and porosity
图5
图5
不同含气饱和度下纵横波速度与孔隙度的关系
a—纵波速度与孔隙度的关系;b—横波速度与孔隙度的关系
Fig.5
Velocity corresponding to porosity at different gas saturations
a—P-velocity corresponding to different porosity;b—S-velocity corresponding to different porosity
w1和w2分别为琼东南海域不含水合物井与高饱和度水合物井,图6 a、b分别为w1(低饱和度)与w2(高饱和度)孔隙度与速度交会图,可以看到,w1纵波速度与孔隙度呈指数型负相关,与模型分析结果一致。而w2中孔隙度与速度负相关关系不明显,这是水合物赋存在不同孔隙度的储层造成的,且高饱和度水合物层未表现出高速异常(图6 b中黑色椭圆),推测是孔隙度较大使得速度减小造成的;高速异常区与水合物饱和度关系不明显,速度无法区分开高饱和度水合物与较高饱和度的水合物(图6 b中红色椭圆),反而与孔隙度关系较明显。那么,这种情况下若要高精度识别水合物,需要去除孔隙度的影响。
图6
图6
南海琼东南海域典型高饱和度水合物与低饱和度水合物井速度—孔隙度测井交会分析
a—w1纵波速度与孔隙度的关系;b—w2横波速度与孔隙度的关系
Fig.6
Typical high saturation hydrate and P-velocity corresponding to porosity with a no hydrate-bearing well in Qiongdongnan area of South China Sea
a—P-velocity corresponding to porosity in w1 well;b—S-velocity corresponding to porosity in w2 well
当速度、孔隙度、水合物饱和度之间保持上述的关系时,我们需要确定此时其AVO特征的正确性,即需要保证模型正演出的AVO属性与实际道集分析结果一致。利用w1井确定不含水合物层的参数,如泥质含量、孔隙度等,含水合物地层的弹性参数通过改变水合物饱和度再建模获取。基于上述的岩石物理模型,从理论上分析了反射系数与水合物饱和度、入射角之间的关系,如图7 所示。可以看到,水合物顶为第一类AVO,AVO截距P 属性为正值,振幅随着角度的增加而逐渐减小,但梯度G 属性较小,而且G 随着水合物饱和度增大而减小,即高饱和度水合物的顶AVO特征不明显。而对于水合物底则呈第三类AVO,P 为负值,G 为负值,振幅绝对值不断增大,且G 绝对值随着水合物饱和度增大而减小,相比于水合物顶,水合物底的P 、G 属性更加明显。
图7
图7
水合物顶底AVO分析
a—水合物顶AVO特征;b—水合物底AVO特征
Fig.7
Modeled AVO at the hydrate top and bottom
a—AVO at the hydrate top;b—AVO at the hydrate bottom
利用实际井对水合物顶底进行了AVO分析,如图8 所示。在顶底附近10个采样点分别计算了其上下两点的AVO特征,发现在水合物顶有明显的一类AVO现象,底层有明显的三类AVO现象,与实际道集分析结果(图中黑色虚线)吻合,其余层位均无此现象(图中蓝色曲线)。这一特征与建模结果吻合,说明上述模型很好地刻画了水合物的AVO异常,进一步证实了模型可靠性。
图8
图8
w3井水合物顶底AVO特征曲线
a—水合物底AVO特征;b—水合物顶AVO特征
Fig.8
AVO at the hydrate top and bottom in w3 well
a—AVO at the hydrate bottom;b—AVO at the hydrate top
2 水合物指示因子构建与反射特征方程
2.1 水合物指示因子构建:孔隙度弹性解耦
由上述分析可知,对于研究区而言,纵波速度与孔隙度呈指数型,与横波速度呈线性。构建水合物指示因子的关键在于:从速度中去除孔隙度影响,突出水合物以及下浮游离气作用。水合物富集使得纵、横波速度增大,即:
(4) v p = v p 0 + v p gh v s = v s 0 + v s gh
式中:v p 、v s 为实测纵波速度、横波速度; v p 0 v s 0 v p gh v s gh
(5) v p = a 1 ϕ b 1 + v p gh v s = a 2 ϕ + b 2 + v s gh ' = a 2 ϕ + v s gh
式中: v s gh v p gh a 1 、a 2 、b 1 可以通过不含水合物井段数据拟合得到。然而,式(5)中 v s gh v p gh 图9 所示,水合物饱和度与孔隙度呈分段线性正相关,在低饱和度或者高饱和度水合物层,水合物饱和度随着孔隙度的变化梯度是几乎不变的,超高饱和度水合物则对应最大的孔隙度,因此 v p gh ϕn 、 v s gh ' ϕn 更能够表示水合物富集特征。那么式(5)可以写成:
(6) v p · ϕ n = a 1 ϕ b 1 + n + v p gh · ϕ n v s · ϕ n = a 2 ϕ 1 + n + v s gh · ϕ n
图9
图9
水合物饱和度与孔隙度的关系
a—w2水合物饱和度与孔隙度的关系;b—w3水合物饱和度与孔隙度的关系
Fig.9
Relation between hydrate saturation and porosity
a—relation between hydrate stauration and porosity in w2 well;b—relation between hydrate stauration and porosity in w3 well
由于水合物富集使得纵横波速度均增大,而游离气对横波几乎无影响,那么可构建如式(7)所示的水合物指示因子:
(7) F p = v p gh · ϕ n + v s gh · ϕ n = v p · ϕ n - a 1 ϕ b 1 + n + v s · ϕ n - a 2 ϕ 1 + n
式中:Fp 为水合物指示因子。取n =1,利用w2井进行交会分析,图10 为Fp 与纵波速度、纵波阻抗的交会图,颜色轴为水合物饱和度,高饱和度水合物层Fp 较大,高饱和度水合物层的Fp 位于100~ 300之间,水合物饱和度越高,Fp 越大,不含水合物层的Fp 低于50,而纵波速度、纵波阻抗等参数未能区分不同饱和度水合物,Fp 识别效果明显优于纵波阻抗等常规参数。下一步则需要考虑如何从地震数据中得到Fp ,本文通过建立以Fp 表示的AVO反射系数方程进行直接反演。
图10
图10
w2井Fp 与纵波速度以及纵波阻抗相对于水合物饱和度的交会分析
a—Fp 与纵波速度的交会;b—Fp 与纵波阻抗的交会
Fig.10
Crossplots of Fp and P-velocity and P-impedance in w2 well
a—crossplot of Fp and P-velocity;b—crossplot of Fp and P-impedance
2.2 水合物指示因子表示的AVO反射特征方程
下一步需要求取Fp 。首先需要建立起与之有关的反射特征方程,设 v p ϕ v p ·ϕn ; v s ϕ v s ·ϕn ,推导可得到(推导过程见附录A):
(8) R pp ϕ ( θ ̅ ) = R pp ( θ ̅ ) + n Δ ϕ ϕ n se c 2 θ ̅ 2 + 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅
(9) R pp ϕ ( θ ̅ ) ≈ se c 2 θ ̅ 2 Δ v p ϕ v p ϕ - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ v s ϕ v s ϕ + 1 2 1 - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ ρ ρ ̅
式(9)为孔隙耦合方程,与孔隙度无关,可用于反演 v p ϕ v s ϕ Fp 。孔隙度的反演方法见附录B。式(8)本质上是对原始反射系数(原始道集)进行预处理,孔隙度引起的速度变化则在等式左边进行了补偿或者删减。若孔隙度较大引起速度减小,则会使得 R pp ϕ
3 基于地震岩石物理的的水合物识别
从式(8)可以看到,反演的精度与孔隙度准确程度密切相关。若n 取值太大,则反演结果对于孔隙度预测结果依赖性较大, 因此反演结果受到更多方面的干扰。n =0不能很好地达到构建水合物指示因子的目的。综合权衡后,本文取n =1,此时Fp 仍具有较好的识别效果(图10 )。建立起如下的反演流程,如图11 所示。首先利用Kumar 等提出的方法,基于纵波阻抗求取孔隙度;然后利用背景速度与孔隙度的关系,求取式(2)中的常数;进一步利用孔隙耦合方程式(8)获得,最后利用式(7)计算出水合物指示因子。
图11
图11
水合物指示因子反演流程
Fig.11
Flow chart of hydrate indicator inversion
实际资料来源于中国南海琼东南海域某二维测线,海底深度约为1 700 m。对实际资料进行叠前反演之前,需要对地震数据进行保幅处理,包括精细的波前扩散补偿、震源组合与检波器组合效应的校正、反Q滤波、地表一致性处理、叠前去噪处理、去除多次波等,并假设处理后的层间多次波、各向异性的影响可以忽略不计。
通过钻探取心已经证实了该区分布大量渗透型高饱和度水合物。图12 展示了该测线不同中心入射角的部分角度叠加剖面, 大约在CDP 380 处, 可以看到从海底到2 500 ms范围内发育明显的渗透通道,底部存在气体模糊带,同相轴下拉,疑似为气体运移所致,气烟囱顶部AVO异常明显。
图12
图12
琼东南研究区某测线地震剖面
a—小角度部分角度叠加数据(1°~10°);b—中角度部分角度叠加数据(11°~20°);c—大角度部分角度叠加数据(21°~30°)
Fig.12
The target seismic section in Qiongdongnan area
a—small angle partial angle stack seismic data(1°~10°);b—middle angle partial angle stack seismic data(11°~20°);c—big angle partial angle stack seismic data(21°~30°)
依据本文提出的方法,首先求取孔隙度。通过孔隙度反演,得到了图13 所示的反演结果,最终的孔隙度结果与测井孔隙度保持很好的一致性。
图13
图13
孔隙度反演结果
Fig.13
Inverted results of porosity
然后进行孔隙耦合方程的反演,图14 展示了 v p ϕ v s ϕ 图 14 已经可以大致看出水合物与游离气分布,黑色区域为疑似游离气分布区,水合物层则分布在游离气顶部。图15 为最终的Fp 反演结果。图中红色井曲线为水合物饱和度曲线,可以看到Fp 剖面中高值区域与井位置处水合物高饱和度保持了很好的一致性,由前文分析可知,高饱和度水合物Fp 介于100~ 300,即图中黄色和红色区域,测井显示5处含水合物层段均已识别出来,下方游离气层位低值异常,水合物和游离气区分性较好,并在一定程度上预测了水合物的空间展布。
图14
图14
v p ϕ ( a )与v ϕ s ( b)反演结果
Fig.14
Inverted results of v ϕ p ( a) and v ϕ s ( b)
图15
图15
Fp 反演结果
Fig.15
Inverted results of Fp
4 结论与认识
多种因素共同影响水合物储层速度,使得水合物预测多解性增强。对于浅层水合物储层而言,水合物饱和度、孔隙度为影响最为显著的参数。利用本文提出的水合物指示因子,可以消除孔隙度的干扰。实际钻探数据分析表明,相比于纵波阻抗等参数,水合物指示因子与水合物饱和度存在更高相关性。高饱和度水合物的指示因子位于100~300内,水合物指示因子反演结果与钻探认识吻合较好,可以有效地预测水合物空间分布,验证了方法的实用性与有效性。
附录A
Aki-Richards近似是AVO反演的基础,其描述了地震反射界面两侧能量分布,Aki近似方程表达为:
(10) R pp ( θ ̅ ) ≈ se c 2 θ ̅ 2 Δ v p v ̅ p - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ v s v ̅ s + 1 2 1 - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ ρ ρ ̅
式中:Δv p 、Δv s 、Δρ 分别表示纵波速度、横波速度和密度的变化量; θ ̅ ρ ̅ v ̅ s v ̅ p v p ϕ v p ·ϕ n ; v s ϕ v s ·ϕ n ,则有:
(11) Δ v p ϕ v p ϕ = Δ v p v p + n Δ ϕ ϕ n Δ v s ϕ v s ϕ = Δ v s v s + n Δ ϕ ϕ n
将式(11)代入式(10),Aki近似方程可以写成:
(12) R pp ( θ ̅ ) ≈ se c 2 θ ̅ 2 Δ v p ϕ v p ϕ - n Δ ϕ ϕ n - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ v s ϕ v s ϕ - n Δ ϕ ϕ n + 1 2 1 - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ ρ ρ ̅
(13) R pp ( θ ̅ ) ≈ se c 2 θ ̅ 2 Δ v p ϕ v p ϕ - se c 2 θ ̅ 2 n Δ ϕ ϕ n - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ v s ϕ v s ϕ + 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ n Δ ϕ ϕ n + 1 2 1 - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ ρ ρ ̅
(14) R pp ϕ ( θ ̅ ) ≈ se c 2 θ ̅ 2 Δ v p ϕ v p ϕ - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ v s ϕ v s ϕ + 1 2 1 - 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ Δ ρ ρ ̅
(15) R pp ( θ ̅ ) ≈ R pp ϕ ( θ ̅ ) - se c 2 θ ̅ 2 + 4 v ̅ s 2 v ̅ p 2 si n 2 θ ̅ n Δ ϕ ϕ n
式(15)表示了Aki近似方程R pp 与去除孔隙度影响后的反射系数方程 R pp ϕ
附录B
Fp 求取需要获取孔隙度,本文是Kumar等[18 ] 提出的方法的优化,孔隙度反射系数与纵波阻抗反射系数的线性关系求解孔隙度:
(16) R pro = a R I p
其中, R I p I p 2 - I p 1 ) / ( I p 2 I p 1 ) 。
利用w1井进行拟合(图16 a),w1线性关系明显,拟合关式为式(17),但是w2由于水合物存在而线性关系较差(图16 b)。
(17) R pro = - 1.14 R I p
图16
图16
南海琼东南海域阻抗与孔隙反射系数关系拟合
a—w1井纵波阻抗与孔隙反射系数关系;b—w2井纵波阻抗与孔隙反射系数关系
Fig.16
Relation between P-impedance and porosity reflection coefficient in Qiongdongnan area of South China Sea
a—relation between P-impedance and porosity in w1 well;b—relation between P-impedance and porosity in w2 well
由于道积分获得的是相对孔隙度变化,积分后结果低频部分不准确,因此本文中最终孔隙度预测结果为积分结果减去低频然后加上实测孔隙度低频部分:
(18) ϕ = ϕ INT - ϕ INT LFM + ϕ well LFM
式中:ϕ INT 、 ϕ INT LFM ϕ well LFM 图17 、18分别为w1和w2孔隙度预测结果,可以看到由于w1井阻抗与孔隙度拟合关系较好,预测结果与实测结果几乎一致,低频趋势差别也较小。而w2积分计算孔隙度相对变化趋势与实测孔隙度几乎一致,但是低频部分的差距在1 790 m处开始增大,局部极值引起的误差需要实测低频来补偿,因此按照式(18)得到的w2最终孔隙度与实测孔隙度差别较小,满足后续反演要求精度。
图17
图17
w1井孔隙度预测结果
Fig.17
Inverted porosity in w1 well
图18
图18
w2井孔隙度预测结果
Fig.18
Inverted porosity in w2 well
参考文献
View Option
[1]
Dickens G R Quinby H M S . Methane hydrate stability in seawater
[J]. Geophysical Research Letters , 1994 ,21 (19 ):2115 -2118 .
DOI:10.1029/94GL01858
URL
[本文引用: 1]
[2]
Kvenvolden K A . Gas hydrates:Geological perspective and global change
[J]. Review in Geophysics , 1993 ,31 :173 -187 .
[本文引用: 1]
[3]
Collett T S . Gas hydrate resources of the United States
[G]// Gautier D L,Dolton G L,Takahashi K I,et al.National Assessment of United States Oil and Gas Resources on CD-ROM,U.S. , 1995 .
[本文引用: 1]
[4]
梁金强 , 付少英 , 陈芳 , 等 . 南海东北部陆坡海底甲烷渗漏及水合物成藏特征
[J]. 天然气地球科学 , 2017 ,28 (5 ):761 -770 .
[本文引用: 1]
Liang J Q Fu S Y Chen F , et al . Characteristics of methane leakage and hydrate accumulation in the continental slope of the northeastern South China Sea
[J]. Natural Gas Geoscience , 2017 ,28 (5 ):761 -770 .
[本文引用: 1]
[5]
杨胜雄 , 梁金强 , 陆敬安 , 等 . 南海北部神狐海域天然气水合物成藏特征及主控因素新认识
[J]. 地学前缘 , 2017 ,24 (4 ):1 -14 .
[本文引用: 1]
Yang S X Liang J Q Lu J A , et al . New understanding of gas hydrate accumulation characteristics and main controlling factors in Shenhu sea area in the northern South China Sea
[J]. Earth Science Frontiers , 2017 ,24 (4 ):1 -14 .
[本文引用: 1]
[6]
郭依群 , 杨胜雄 , 梁金强 , 等 . 南海北部神狐海域高饱和度天然气水合物分布特征
[J]. 地学前缘 , 2017 ,24 (4 ):24 -31 .
[本文引用: 1]
Guo Y Q Yang S X Liang J Q , et al . Distribution characteristics of high saturation gas hydrate in Shenhu sea area in the northern South China Sea
[J]. Geoscience , 2017 ,24 (4 ):24 -31 .
[本文引用: 1]
[7]
Waite W Helgerud M B Nur A , et al . Laboratory measurements of compressional and shear wave speeds through methane hydrate
[G]// Holder G D,Bishnoi P R. Gas hydrates:Challenges for the future:Annals of the New York Academy of Science , 2000 : 1003 -1010 .
[本文引用: 1]
[8]
Hamilton E L . Geoacoustic modeling of the sea floor
[J]. Journal of the Acoustical Society of America , 1980 ,68 :1313 -1340 .
DOI:10.1121/1.385100
URL
[本文引用: 1]
[9]
胡高伟 , 业渝光 , 张剑 , 等 . 沉积物中天然气水合物微观分布模式及其声学响应特征
[J]. 天然气工业 , 2010 ,30 (3 ):120 -124 .
DOI:10.3787/j.issn.10000976.2010.03.031
URL
[本文引用: 1]
天然气水合物（以下简称水合物）的微观分布与其形成条件、流体运移通道等密切相关，对天然气水合物资源勘探与评价具有重要意义。为了解水合物在固结沉积物和松散沉积物中的微观分布及其声学响应特征，采用超声和时域反射联合探测技术实时测量了沉积物中水合物生成、分解过程中声速等参数的变化情况。结果表明，在固结沉积物中，水合物先在孔隙流体中形成，随后逐渐向骨架靠拢，当饱和度大于30%后水合物开始胶结沉积物颗粒生成，这种胶结模式会圈闭部分孔隙流体，使之因得不到气源的补充而难以形成水合物，因此固结沉积物中水合物饱和度最终为65.5%左右；在松散沉积物中，少量的水合物（饱和度1%左右）胶结沉积物颗粒生成，当饱和度大于1%后水合物开始在孔隙流体中以悬浮状形态生成，由于水合物与沉积物颗粒间尚有流体运移通道，水合物能进一步生成，最终几乎完全充填沉积物孔隙。不同的水合物微观分布特征对沉积物的声速具有不同影响：水合物在孔隙流体中生成时，10%饱和度的水合物对固结沉积物的声速影响不明显；当水合物胶结沉积物颗粒生成时，约1%饱和度的水合物可使松散沉积物的纵波速度增长200~300 m/s。
Hu G W Ye Y G Zhang J , et al . Microscopic distribution patterns and acoustic response characteristics of natural gas hydrates in sediments
[J]. Natural Gas Industry , 2010 ,30 (3 ):120 -124 .
[本文引用: 1]
[10]
胡高伟 , 李承峰 , 业渝光 , 等 . 沉积物孔隙空间天然气水合物微观分布观测
[J]. 地球物理学报 , 2014 ,57 (5 ):1675 -1682 .
DOI:10.6038/cjg20140530
URL
[本文引用: 1]
在高分辨率工业型X射线层析扫描仪上（X-CT，X-Ray Computerized Tomography），研制了适用于沉积物中水合物微观赋存状态研究的CT原位探测装置，可通过模拟天然气水合物生成/分解过程，直接观测水合物在沉积孔隙中的微观分布特征.对粒径为0.425~0.85 mm沉积物中水合物形成过程进行了CT观测研究，结果表明，沉积孔隙中水合物呈混合分布模式，但在水合物不同形成阶段，以某种分布模式为主导：在水合物形成初期，仅有极少量水合物悬浮在流体中，水合物主要以接触或胶结模式为主；在水合物形成中期（如饱和度为24.6%、35.0%时），水合物倾向于在孔隙流体中以悬浮状形态生成；在水合物形成后期（如饱和度为51.4%之后），悬浮状的水合物慢慢生长聚拢在一起，水合物又重新胶结沉积物颗粒.
Hu G W Li C F Ye Y G , et al . Microscopic distribution observation of gas hydrate in sediment pore space
[J]. Chinese Journal of Geophysics , 2014 ,57 (5 ):1675 -1682 .
[本文引用: 1]
[11]
潘豪杰 , 刘堂晏 , 邬龙 , 等 . 天然气水合物赋存形态识别及其饱和度预测
[J]. 地球物理学进展 , 2014 ,29 (4 ):1735 -1740 .
DOI:10.6038/pg20140432
URL
[本文引用: 1]
含水合物储层纵横波速度不仅与天然气水合物饱和度有关，还取决于水合物的赋存形态.本文基于考虑水合物微观形态的岩石物理模型预测了水合物海脊1247B和1250F井位的水合物饱和度，并与核磁共振——密度孔隙度评价的饱和度进行比较，以确定水合物在沉积物中的存在形态.而且对比了有效介质模型、改进的Boit-Gassmann模型和简化的三相方程在同一赋存形态下预测的饱和度，以此探究三种模型在含水合物储层定量评价中的适用性.对比预测结果显示，1247B和1250F井位的天然气水合物主要以骨架支撑形态存在.虽然各岩石物理模型在含水合物层段预测趋势一致，但是，相对于简化三相方程预测结果而言，有效介质模型和改进的Biot-Gassmann模型能更准确预测出海底沉积物中水合物饱和度，并且在同一种模型中纵横波速度联合反演比单纵波预测的水合物饱和度更合理.
Pan H J Liu T Y Wu L , et al . Gas hydrate occurrence pattern recognition and saturation prediction
[J]. Progress in Geophysics , 2014 ,29 (4 ):1735 -1740 .
[本文引用: 1]
[12]
Helgerud M B Dvorkin J Nur A , et al . Elastic-wave velocity in marine sediments with gas hydrates:Effective medium modeling
[J]. Geophysical Research Letters , 1999 ,26 :2021 -2024 .
DOI:10.1029/1999GL900421
URL
[本文引用: 1]
[13]
Sava D Hardage B . Rock-physics characterization of hydrate-bearing deepwater sediments
[J]. The Leading Edge , 2006 ,25 :616 -619 .
DOI:10.1190/1.2202666
URL
[本文引用: 1]
[14]
Domenico S N . Effect of brine-gas mixture on velocity in an unconsolidated sand reservoir
[J]. Geophysics , 1976 ,41 :882 -894 .
DOI:10.1190/1.1440670
URL
[本文引用: 1]
[15]
王秀娟 , 吴时国 , 王吉亮 , 等 . 南海北部神狐海域天然气水合物分解的测井异常
[J]. 地球物理学报 , 2013 ,56 (8 ):2799 -2807 .
DOI:10.6038/cjg20130828
URL
[本文引用: 1]
a=1.1和m=2.3.基于该参数,利用阿尔奇方程计算的水合物饱和度占孔隙空间5%~20%,局部地层水合物饱和度达26.8%,在垂向上分布不均匀.由于钻探可能导致水合物发生分解而产生游离气,原位游离气和水合物分解产生的气体都能造成低纵波速度异常.由于地震资料采集在测井之前完成,利用不同速度制作合成地震记录并与地震资料进行对比,能够确定水合物稳定带上部的低速异常形成原因.]]>
Wang X J Wu S G Wang J L , et al . Logging anomalies of natural gas hydrate decomposition in the Shenhu sea area in the northern South China Sea
[J]. Chinese Journal of Geophysics , 2013 ,56 (8 ):2799 -2807 .
[本文引用: 1]
[17]
Chatterjee R Singha D K Ojha M , et al . Porosity estimation from pre-stack seismic data in gas-hydrate bearing sediments, Krishna-Godavari basin,India
[J]. Journal of Natural Gas Science and Engineering , 2016 ,33 :562 -572 .
DOI:10.1016/j.jngse.2016.05.066
URL
[本文引用: 1]
[18]
Kumar R Das B Chatterjee R , et al . A methodology of porosity estimation from inversion of post-stack seismic data
[J]. Journal of Natural Gas Science and Engineering , 2016 ,28 (1 ):356 -364 .
DOI:10.1016/j.jngse.2015.12.028
URL
[本文引用: 2]
[19]
Lee M W Collett T S . In-situ gas hydrate saturations estimated from various well logs at Mount Elbert well, North Slope of Alaska
[J]. Marine and Petroleum Geology , 2011 ,28 :439 -449 .
DOI:10.1016/j.marpetgeo.2009.06.007
URL
[本文引用: 1]
Gas hydrate saturations were estimated from various well logs such as nuclear magnetic resonance (NMR), P- and S-wave velocity, and electrical resistivity logs along with pore-water salinity. Gas hydrate saturations from the NMR log agree well with those estimated from P- and S-wave velocity data. Because of the low salinity of the connate water and the low formation temperature, the resistivity of connate water is comparable to that of shale. Therefore, the effect of clay should be accounted for to accurately estimate gas hydrate saturations from the resistivity data. Two highly gas hydrate-saturated intervals are identified - an upper similar to 43 ft zone with an average gas hydrate saturation of 54% and a lower similar to 53 ft zone with an average gas hydrate saturation of 50%; both zones reach a maximum of about 75% saturation. Published by Elsevier Ltd.]]>
[20]
Yin X Y Zong Z Y Wu G C . Research on seismic fluid identification driven by rock physics
[J]. Science China Earth Sciences , 2015 ,58 (2 ):159 -171 .
DOI:10.1007/s11430-014-4992-3
URL
[本文引用: 1]
[21]
Aki K Richards P G . Quantitative seismology
[M]. California:University Science Books , 2002 .
[本文引用: 1]
[23]
Smith G C Gidlow P M . Wighted stacking for rock property estimation and detection of gas
[J]. Geophysical Prospecting , 1987 ,35 (9 ):992 -1014 .
[本文引用: 1]
[24]
Fatti J L Smith G C Vail P J , et al . Detection of gas in sandstone reservoirs using AVO analysis:A 3-D seismic case history using the Geostack technique
[J]. Geophysics , 1994 ,59 (9 ):1362 -1376 .
DOI:10.1190/1.1443695
URL
[本文引用: 1]
[25]
Russell B H Gray D Hampson D P . Linearized AVO inversion and poroelasticity
[J]. Geophysics , 2011 ,76 (3 ):C19 -C29 .
DOI:10.1190/1.3555082
URL
[本文引用: 1]
The technique of amplitude variation with offset (AVO) allows geoscientists to extract fluid and lithology information from the analysis of prestack seismic amplitudes. Various AVO parameterizations exist, all of which involve the sum of three weighted elastic-constant terms. In present-day AVO approaches, the weighting terms involve either knowledge of the incidence angle only, or knowledge of both the incidence angle and the in situ V(P)/V(S) ratio. We have used the theory of poroelasticity to derive a generalized AVO approximation that provides the estimation of fluid, rigidity, and density parameters. We have combined two previously independent AVO formulations, thus reducing, instead of adding to, the total number of formulations. This new approach requires knowledge of a third parameter to compute the weights: the dry-rock V(P)/V(S) ratio. We have derived a new equation and applied it to model and real data sets. The new formulation has allowed us to estimate fluid properties of the reservoir in a more direct manner than previous formulations.
[26]
Yin X Y Zhang S X . Bayesian inversion for effective pore-fluid bulk modulus based on fluid-matrix decoupled amplitude variation with offset approximation
[J]. Geophysics , 2014 ,79 (5 ):R221 -R232 .
DOI:10.1190/GEO2013-0372.1
URL
[本文引用: 1]
Fluid indicators estimated from seismic data play important roles in reservoir characterization and prospect identification. Traditionally, there are a variety of fluid indicators proposed, but they are very likely to provide ambiguous results for fluid identification due to the fact that their sensitivity is dependent upon the mixed effect of pore fluid and rock porosity. To raise the sensitivity of fluid indication, we used the effective pore-fluid bulk modulus as a fluid indicator. Starting with the poroelastic amplitude variation with offset (AVO) theory and the corresponding rock-physics model with the homogeneous sorting trend, we derived a new AVO approximation that allowed us to estimate the effective pore-fluid bulk modulus in a direct fashion. The inversion for the fluid indicator is formulated in Bayesian framework with the Cauchy distribution as a prior constraint. We tested the method on synthetic data and analyzed the feasibility and stability of the inversion. A field data example shows that the effective pore-fluid bulk modulus can reduce the ambiguity caused by the rock porosity and improve the quality of fluid discrimination in a clastic reservoir. Further research needs to be done on the reservoirs that do not fit the rock-physics model without a sorting trend.
[28]
邓炜 , 印兴耀 , 宗兆云 . 基于等效流体体积模量直接反演的流体识别方法
[J]. 石油地球物理勘探 , 2017 ,52 (2 ):315 -325 .
[本文引用: 1]
Deng W Yin X Y Zong Z Y . Fluid identification method based on direct inversion of equivalent fluid volume modulus
[J]. Oil Geophysical Prospecting , 2017 ,52 (2 ):315 -325 .
[本文引用: 1]
[29]
沙志彬 , 万晓明 , 赵忠泉 , 等 . 叠前同时反演技术在珠江口盆地西部海域天然气水合物储层预测中的应用
[J]. 物探与化探 , 2019 ,43 (3 ):476 -485 .
[本文引用: 1]
Sha Z B Wan X M Zhao Z Q , et al . Application of pre-stack simultaneous inversion technique in gas hydrate reservoir prediction in the western waters of the Pearl River Mouth Basin
[J]. Geophysical and Geochemical Exploration , 2019 ,43 (3 ):476 -485 .
[本文引用: 1]
[30]
Hill R . A self-consistent mechanics of composite materials
[J]. Journal of the Mechanics and Physics of Solids , 1965 ,13 :213 -222 .
DOI:10.1016/0022-5096(65)90010-4
URL
[本文引用: 1]
Methane hydrate stability in seawater
1
1994
... 天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内.地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] . ...
Gas hydrates:Geological perspective and global change
1
1993
... 天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内.地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] . ...
Gas hydrate resources of the United States
1
1995
... 天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内.地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] . ...
南海东北部陆坡海底甲烷渗漏及水合物成藏特征
1
2017
... 天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内.地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] . ...
南海东北部陆坡海底甲烷渗漏及水合物成藏特征
1
2017
... 天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内.地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] . ...
南海北部神狐海域天然气水合物成藏特征及主控因素新认识
1
2017
... 天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内.地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] . ...
南海北部神狐海域天然气水合物成藏特征及主控因素新认识
1
2017
... 天然气水合物是在一定的温度和压力下由水和天然气形成的类似冰的笼型结晶化合物[1 ,2 ,3 ] , 主要蕴藏在浅部海底盆地、陆架坡折带和陆地永久冻土区等的天然气水合物稳定域内.地震识别是检测海底水合物以及游离气最有效的方法之一[4 ,5 ] . ...
南海北部神狐海域高饱和度天然气水合物分布特征
1
2017
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
南海北部神狐海域高饱和度天然气水合物分布特征
1
2017
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
Laboratory measurements of compressional and shear wave speeds through methane hydrate
1
2000
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
Geoacoustic modeling of the sea floor
1
1980
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
沉积物中天然气水合物微观分布模式及其声学响应特征
1
2010
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
沉积物中天然气水合物微观分布模式及其声学响应特征
1
2010
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
沉积物孔隙空间天然气水合物微观分布观测
1
2014
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
沉积物孔隙空间天然气水合物微观分布观测
1
2014
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
天然气水合物赋存形态识别及其饱和度预测
1
2014
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
天然气水合物赋存形态识别及其饱和度预测
1
2014
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
Elastic-wave velocity in marine sediments with gas hydrates:Effective medium modeling
1
1999
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
Rock-physics characterization of hydrate-bearing deepwater sediments
1
2006
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
Effect of brine-gas mixture on velocity in an unconsolidated sand reservoir
1
1976
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
南海北部神狐海域天然气水合物分解的测井异常
1
2013
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
南海北部神狐海域天然气水合物分解的测井异常
1
2013
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
Fundamentals of seismic rock physics
1
2012
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
Porosity estimation from pre-stack seismic data in gas-hydrate bearing sediments, Krishna-Godavari basin,India
1
2016
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
A methodology of porosity estimation from inversion of post-stack seismic data
2
2016
... 速度在水合物勘查中具有重要意义[6 ] ,但速度的影响因素较多.实验室已测出纯合成水合物速度[7 ] ,即纵波速度3 650 m/s,横波速度1 890 m/s,密度900 kg/m3 ,速度高于海底未成岩的水合物储层.相比于不含水合物的岩石,大量水合物的聚集会使得岩石的纵横波速度增加,泊松比减小[8 ] ,然而少量的水合物对于速度的影响则依赖于水合物赋集方式.沉积物孔隙空间中水合物形成并非单一模式,在同一时刻水合物可能出现“悬浮”、“接触”和“胶结”的混合模式,在形成过程中水合物先以胶结沉积物颗粒的微观分布为主,形成中期以悬浮状形态或接触模式为主,形成后期又以胶结模式为主[9 ,10 ,11 ] .因此高饱和度的水合物应以胶结模式为主,在建模中主要考虑胶结模式的水合物,如果水合物以胶结态赋集,则纵横波速度增大[12 ,13 ] .水合物一般与游离气伴生,游离气对速度的影响更加复杂,即使孔隙中存在微量游离气(<1%),纵波速度也会发生明显变化[14 ] ,含气对横波速度的影响不明显,水合物储层中可能存在分解游离气与水合物共存的现象[15 ] ,此时储层弹性特征较为复杂.孔隙度对于弹性特征影响非常明显[16 ] ,水合物储层孔隙度可以基于孔隙度反射率以及波阻抗反射率、速度反射率之间的关系进行反演获取[17 ,18 ] . ...
... Fp 求取需要获取孔隙度,本文是Kumar等[18 ] 提出的方法的优化,孔隙度反射系数与纵波阻抗反射系数的线性关系求解孔隙度: ...
In-situ gas hydrate saturations estimated from various well logs at Mount Elbert well, North Slope of Alaska
1
2011
... 正是由于影响水合物速度的因素多,单纯利用速度识别水合物可能会带来一定误差[19 ] ,特别是在海底浅层的MTD、河道沉积往往产生与水合物类似的反射特征,此时则需要去除由于孔隙度变化引起的振幅异常,岩石物理建模分析可以有效地指导开展高精度的水合物识别,其关键在于构建高敏感度的水合物指示因子.事实上,在实际生产应用中,针对不同储层应使用不同的弹性参数或者构建相应的指示因子[20 ] ,叠后阻抗、纵横波速度、泊松比、Gassman流体因子、杨氏模量、泊松比、流体模量、组合流体因子等相继被提出并应用到相关研究中,且国内外众多学者相继提出了不同弹性参数的反演方法,但是对于水合物储层的流体因子构建和反演研究尚少. ...
Research on seismic fluid identification driven by rock physics
1
2015
... 正是由于影响水合物速度的因素多,单纯利用速度识别水合物可能会带来一定误差[19 ] ,特别是在海底浅层的MTD、河道沉积往往产生与水合物类似的反射特征,此时则需要去除由于孔隙度变化引起的振幅异常,岩石物理建模分析可以有效地指导开展高精度的水合物识别,其关键在于构建高敏感度的水合物指示因子.事实上,在实际生产应用中,针对不同储层应使用不同的弹性参数或者构建相应的指示因子[20 ] ,叠后阻抗、纵横波速度、泊松比、Gassman流体因子、杨氏模量、泊松比、流体模量、组合流体因子等相继被提出并应用到相关研究中,且国内外众多学者相继提出了不同弹性参数的反演方法,但是对于水合物储层的流体因子构建和反演研究尚少. ...
Quantitative seismology
1
2002
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
A simplification of the Zoeppritz equations
1
1985
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
Wighted stacking for rock property estimation and detection of gas
1
1987
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
Detection of gas in sandstone reservoirs using AVO analysis:A 3-D seismic case history using the Geostack technique
1
1994
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
Linearized AVO inversion and poroelasticity
1
2011
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
Bayesian inversion for effective pore-fluid bulk modulus based on fluid-matrix decoupled amplitude variation with offset approximation
1
2014
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
Elastic impedance
1
1999
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
基于等效流体体积模量直接反演的流体识别方法
1
2017
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
基于等效流体体积模量直接反演的流体识别方法
1
2017
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
叠前同时反演技术在珠江口盆地西部海域天然气水合物储层预测中的应用
1
2019
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
叠前同时反演技术在珠江口盆地西部海域天然气水合物储层预测中的应用
1
2019
... 叠前地震反演的基础是Zoeppritz方程,Aki等[21 ] 导出了在界面两侧弹性差异较小的情况下纵波反射系数的近似表达式,该方程广泛用于反演纵横波速度;Shuey[22 ] 导出了与泊松比有关的纵波相对反射系数近似式,用于含气性检测;Smith等[23 ] 将Gardner公式代入Aki近似式中,消除了密度项,提高反演的稳定性;Fatti等推导了基于纵横波阻抗和密度的反射系数近似方程,该方程提高了纵波阻抗反演稳定性,被广泛应用于工业生产[24 ] ;Russell等[25 ] 基于Boit-Gassmann理论对饱和流体多孔介质的AVO理论进行研究,提出了包含Gassmann流体项的反射系数近似公式,一定程度消除了岩石骨架的影响;Yin和Zhang[26 ] 针对复杂孔隙结构储层推导了基于流体模量Kf 表示的固液解耦四项近似近似.为了提高反演的稳定性,本文使用弹性阻抗反演的方法.弹性阻抗反演由于其固有的优秀抗噪性得到广泛使用,国内外众多学者在上述AVO近似式基础上提出了不同的弹性阻抗反演方法[27 ,28 ] ,在此不一一赘述.但浅层水合物储层尚未成岩,其弹性特征受到孔隙度的影响较大,水合物饱和度与孔隙度对速度等弹性参数的影响程度不明确,需要构建对水合物敏感性更好的参数[29 ] ,并发展相应的反演方法. ...
A self-consistent mechanics of composite materials
1
1965
... 笔者基于SCA-DEM岩石物理建模方法[30 ] ,针对水合物未固结特性合理设置建模参数,定量地描述物性参数对弹性参数以及AVO 特征的影响,并在此基础上构建了水合物指示因子.相比于传统弹性参数,该指示因子与水合物饱和度具有更高相关性.然后在平面波近似下,推导了与该指示因子直接相关的AVO反射特征方程,基于该方程则可直接获取水合物指示因子,提高水合物识别精度. ...
