Estimated gas hydrate saturation from the reservoir of clayey silt with sandy interlayers at Site U1517, Tuaheni landslide complex on the Hikurangi margin, New Zealand
-
摘要: 准确评估新西兰Hikurangi边缘Tuaheni滑坡复合体（TLC）区域的天然气水合物含量与储层分布对TLC慢滑移现象与产生机制的解释有重要作用。本文分析了IODP372航次U1517站位测井和取心数据，发现在局部地层纵波速度增加（＞1.7 km/s）和电阻率升高（＞1.5 Ω·m）的104～160 mbsf层段存在天然气水合物，其中112～114、130～145和150～160 mbsf层段饱和度相对较高。根据岩性划分了不同井段对应的矿物成分含量，用于纵波速度模型计算，并利用简化三相介质（STPE）和改进的Biot-Gassmann模型（BGTL）分别估算了104～160 mbsf层段的天然气水合物饱和度，平均饱和度分别为5.2%和6.0%，最高饱和度分别为22.7%和21.6%。同时，与阿尔奇公式估算的水合物饱和度比较，在104～160 mbsf层段3种方法估算的饱和度值随深度变化相似，天然气水合物平均饱和度相近（约6.0%），在130～145 mbsf层段的水合物平均饱和度最高（约8.5%）。本研究使用两种声速模型和更为精细的参数估算饱和度，其估算结果更为可靠，可为Tuaheni滑坡复合体慢滑移现象研究提供良好的基础数据支撑。
-
关键词:
- 天然气水合物 /
- 水合物饱和度 /
- STPE和BGTL模型 /
- 新西兰Hikurangi边缘 /
- IODP372航次
Abstract: Accurate assessment of natural gas hydrate saturation and reservoir distribution in the Tuaheni landslide complex (TLC) area on the Hikurangi margin, New Zealand plays a critical role in explaining the creeping TLC phenomenon and its forming mechanisms. Gas hydrates are discovered in the interval from 104 mbsf to the BGHS at 160 mbsf based on the logging- while-drilling (LWD) and coring data from the site U1517 located on the extensional and creeping part of the TLC. Elevated P-wave velocity (＞1.7 km/s), electrical resistivity (＞1.5 Ω·m), and the high saturation intervals occur in 112～114 mbsf, 130～145 mbsf and 150～160 mbsf respectively. The mineral components of different well intervals are used to calculate the P-wave velocity, to estimate the gas hydrate saturation of the 104～160 mbsf interval with the simplified three-phase equation (STPE) and the biot- gassmann theory by Lee (BGTL) models, and the average saturation was 5.2% and 6.0% and the highest was 22.7% and 21.6%, respectively. Adoption of the new method made it more efficient to get BGTL model parameters. Compared with the hydrate saturation estimated by the Expedition 372 from the Archie equation was used, the average gas hydrate saturation is similar in the 104～160 mbsf interval (about 6.0%) and about 8.5% in the interval of 130～145 mbsf. In this study, two sound velocity models were used to estimate the hydrate saturation of U1517, making the results more reliable. The accurate estimation of the gas hydrate saturation and distribution at site U1517 will provide basic data for simulation of the creeping TLC phenomenon on the Hikurangi margin.-
Keywords:
- gas hydrate /
- hydrate saturations /
- STPE and BGTL /
- the Hikurangi margin /
- IODP 372
-
天然气水合物是一种笼状结构的类冰状结晶化合物，主要是由甲烷和水分子结合而成，因其在冻土地区和海洋大陆边缘广泛分布、与海底稳定性相关，以及可能对全球气候具有潜在影响而广受关注[1-2]。新西兰Hikurangi大陆边缘每两年左右发生一次慢滑移事件[3]，有关证据显示，多期次水合物形成分解可能是造成该区产生蠕变的重要原因之一[4]。
2017年11月—2018年1月执行了“蠕变中的天然气水合物滑动和Hikurangi随钻测井”为主旨的IODP372航次。该航次的主要目的之一是调查天然气水合物和海底滑坡的关系，因此，在新西兰Hikurangi边缘Tuaheni滑坡复合体（Tuaheni Landslide Complex，TLC）的U1517站位进行了随钻测井工作（图1）。该站位钻井的主要任务是通过在滑坡体和天然气水合物稳定区进行测井和采样，研究水合物与蠕变的关系。
![]()
20世纪末，研究人员在新西兰Hikurangi边缘发现地震高振幅异常和似海底反射（BSR）标志[4-7]。该区域反射地震[8-11]、电磁[12-14]、甲烷渗漏、海水甲烷浓度、与渗漏相关的沉积坍塌和冷泉等证据均指示有天然气水合物存在[15-19]。通过反射地震发现，Hikurangi大陆边缘的Tuaheni滑坡复合体显示了活动蠕变变形的特征，且蠕变中的近陆边缘与海底天然气水合物稳定带底部的尖灭相一致[4]，因此，科学家认为水合物分解—形成过程可能与新西兰Hikurangi边缘的多期次慢滑移密切相关[4, 20-21]。TLC地区水合物分布和含量估算是研究蠕变与水合物关系的必要环节。Mountjoy提出3种机制解释浅层天然气水合物是如何导致慢滑移，主要认为是由于水合物的分解导致沉积物液化失稳、水合物分解对地层孔隙压力的影响和水合物含量对地层提供的不同支撑模式的影响[4]。不同饱和度的水合物对沉积物的支撑模式不同[22]，因此，准确估算慢滑移区域水合物的饱和度可以进一步分析天然气水合物导致TLC慢滑移的原因，IODP372航次测井和取心为水合物饱和度估算提供了可靠资料。
由于声速和电阻率对水合物储层最为敏感，常用声学、电学模型来估算天然气水合物饱和度[23-24]。基于电阻率的模型有阿尔奇方程[25]和连通性方程[26]；基于声速的模型有权重方程（Weighted equation，WE）[27]、等效介质理论（Effective Media Theory，EMT）模型[28]、改进的Biot-Gassmann理论（Biot-Gassmann Theory by Lee，BGTL）模型[29-30]和简化三相介质方程（Simplified Three-Phase Equation，STPE）[31]等，其中，常用于测井应用的模型主要是STPE和等效介质理论[32]两种。由于STPE模型参数较易获取，所以多次被实际应用于估算天然气水合物饱和度，均得到理想的预测效果[33-34]；Hu等采用超声探测技术和时域反射技术实时探测了沉积物的纵横波速度和水合物饱和度的变化情况，检验了多种理论模型，发现BGTL理论预测的纵、横波速度更接近实测值[35]，但BGTL模型中与岩石固结程度相关参数难以通过实际地层数据计算，从而导致较少被应用于测井数据估算水合物饱和度。本文主要通过STPE与BGTL模型对U1517站位水合物饱和度进行研究，在BGTL模型使用新的参数选取方法，使参数获取更为简易，计算过程中，根据岩性划分不同层段对应的矿物成分含量，用于纵波速度模型计算，以精确模型判断水合物储层深度分布和天然气水合物饱和度计算。
1. 测井数据分析
IODP372航次U1517站位测井位于38°S、178°E（图1），井深约205 mbsf。该航次通过随钻测井采集了井径、声波速度、伽马密度、孔隙度、自然伽马和电阻率等数据，其中纵波数据在160～168 mbsf层段内未获取。通过对LDEO（Lamont Doherty Earth Observatory）数据库提供的U1517站位原始测井数据进行解释分析，并拟合背景趋势线，结果如图2所示，实测纵波声速与背景拟合声速相比，速度增加出现在94～160 mbsf层段，氯离子浓度异常出现在104～160 mbsf层段；在94～104 mbsf层段纵波速度和孔隙度等明显增加，而电阻率与密度减小，氯离子浓度并无异常，该层段的异常与21～28 mbsf层段相似，可能为井孔和局部岩性变化造成（图中黄色区域）。天然气水合物稳定区大概位于104～160 mbsf，并且在130～145 mbsf层段（图中绿色区域），纵波速度、电阻率和井径明显增加，密度减小。该航次从U1517站位获得多个柱状样品，使用红外热像仪进行扫描，温度异常表明天然气水合物的存在，并发现在岩心的上层沉积物或其岩心采集器中有甲烷释放[37]。图3所示为地层因子与纵波速度交会图，由于含天然气水合物的沉积物具有较高的纵波速度和地层因子，所以含天然气水合物的沉积地层的交会图显示高于饱和水沉积地层[34]，在130～145 mbsf层段显示较高的地层因子和纵波速度。
![]()
图 2 U1517站位测井数据蓝色实线为实测数据，红色实线为拟合背景趋势，黑色虚线为BSR。Figure 2. Well logs at Site U1517The blue lines are measured logging data and the red lines are fitted background trends; the black dashed line is BSR.![]()
图 3 U1517站位井地层因子与纵波速度交会图黑色实线为0～90 mbsf拟合背景趋势线。Figure 3. Cross plot of formation factor versus the measured P-wave velocities at Site U1517The black line is fitted background trends of 0～90 mbsf.2. 储层水合物饱和度估算
依据国际大洋发现计划出版物（International Ocean Discovery Program Publications）[36]获取矿物类型及含量数据，基于该数据，根据岩性划分不同层段对应的矿物成分含量，结果如图4所示，用于纵波速度模型计算，以精确天然气水合物饱和度计算值。岩石骨架的不同矿物类型的物性参数如表1所示。
表 1 骨架组分及物性参数Table 1. Constants used for the modeling![]()
图 4 U1517站位所取岩心的岩性和岩石矿物成分相对含量数据黑线为依据岩性平均矿物成分相对含量。Figure 4. Simplified lithostratigraphic column with bulk powder XRD results, Site U1517Black line is average mineral composition based on core data.2.1 简化三相介质模型（STPE）
含天然气水合物储层具有相对较高的纵横波速度。本次研究中使用STPE模拟U1517站位井的纵波速度，其中用于纵波速度（Vp）建模的STPE[31, 33]使用等式（1）对水合物储层的纵波速度建模：
$${V_{\rm{p}}} = \sqrt {\frac{{k + {\rm{4}}\mu /{\rm{3}}}}{\rho }} , \; {V_{\rm{s}}} = \sqrt {\frac{\mu }{\rho }} $$ (1) 式中ρ是含天然气水合物沉积模型的体积密度，k是体积模量，μ是剪切模量。建模参数ε是解释在加强主体沉积物骨架方面天然气水合物形成相对于压实的影响减小，Lee和Waite推荐使用ε=0.12为建模数值[30]。在Vp和Vs建模中使用的参数α使用等式（2）计算：
$${\alpha _{\rm{i}}} = {\alpha _0}{({p_0}/{p_i})^{\rm{n}}} \approx {\alpha _0}{({d_0}/{d_{\rm{i}}})^{\rm{n}}}$$ (2) 式中α0是有效压力p0和深度d0的固结参数，αi是有效压力pi和深度di的固结参数，固结参数可以使用饱和水沉积物的速度来估算[40]。固结参数取决于固结程度和该区域的有效压力，Mindlin认为体积模和剪切模量为有效压力的1/3幂[41]，因此不同位置，根据研究区域的主要岩性，α的值随深度而变化[40]。通过建模速度基线和实测纵波速度之间的最佳拟合选定α的值[40]，本次研究用于U1517站位井的固结参数αi=42（60/di）1/3。使用上述参数，获得了井下剖面背景纵波速度和天然气水合物饱和度，饱和水沉积地层VP符合程度较高（图5，图6），黑色实线为航次实测纵波速度，红色实线为本文利用STPE模型计算结果，如图5所示，104～160 mbsf层段内测井实测VP大于理论基线速度，可能属于天然气水合物储层区，利用航次实测数据和模型结果计算出水合物饱和度（图6）。结果显示，在104～160 mbsf的深度区间内平均饱和度约为5.2%，最高饱和度达到22.7%，其中130～145 mbsf层段内水合物饱和度较高，平均饱和度为7.9%。
![]()
图 5 使用STPE在U1517站位井测量的纵波速度和计算的基线速度的比较Figure 5. Measured and calculated baseline P-wave velocities with STPE at the Site U1517![]()
图 6 使用STPE计算水合物储层区的背景纵波速度及饱和度Figure 6. Measured and calculated baseline resistivities and P-wave velocities at the Site U15172.2 改进的Biot-Gassmann模型（BGTL）
BGTL理论建立在经典的BGT理论（Biot-Gassmann Theory）上，在预测速度时不仅考虑了分压的影响，而且还考虑了岩石的孔隙度、固结度等因素的影响[35]。在BGTL模型计算中，将天然气水合物作为基质中的一种矿物成分。本研究用于VP建模的BGTL模型使用等式（1）进行。
公式（1）中沉积介质的剪切模量μ可由下式计算：
$$ \mu=\frac{\mu_{\rm{ma}}{ G}^2\left(1-\varPhi\right)^{2{ n}}{ k}}{{{k_{\rm{ma}}}}+4\mu_{\rm{ma}}\left[1-{ G}^2\left(1-\varPhi\right)^{2{ n}}\right]/3} $$ (3) 式中，kma为岩石骨架的体积模量；μma为岩石骨架的剪切模量；Φ为孔隙度；常数G主要用来校正由基质中的黏土引起的差异。Han等通过实验室数据表明G = 1对清洁砂岩有利[42]，随着黏土体积增加，G将按下式（4）计算减少：
$${{G}} = 0.955\;2 + 0.044\;8{{\rm{e}}^{ - {C_{\rm{v}}}{\rm{/}}0.067\;14}}$$ (4) 式中，泥质含量Cv可使用来自U1517A井的伽马射线测井数据通过公式（5）[34]估算：
$${C_{\rm{v}}} = 0.083({2^{{\rm{GCUR}} \times {{\rm{I}}_{{\rm{GR}}}}}} - 1)$$ (5) 式中，GCUR是与地层有关的经验系数，新地层(古新近系地层)GCUR=3.7[43]，IGR为通过伽马测井数据计算的伽马射线指数，可由公式（6）计算：
$${{\rm{I}}_{{\rm{GR}}}} = \frac{{{\rm{G}}{{\rm{R}}_{{\rm{log}}}} - {\rm{G}}{{\rm{R}}_{{\rm{min}}}}}}{{{\rm{G}}{{\rm{R}}_{{\rm{max}}}} - {\rm{G}}{{\rm{R}}_{{\rm{min}}}}}}$$ (6) 式中，GRlog为测井伽马值，GRmin为砂岩层伽马值，GRmax为泥岩层伽马值。
公式（3）中参数n取决于分压大小及岩石的固结程度，可由公式（7）得到：
$${{n}} = [{\rm{1}}{{\rm{0}}^{({\rm{0}}.{\rm{426}} - {\rm{0}}.{\rm{235Log10}}p)}}]{{/m}}$$ (7) 测量数据表明m≈5适合于固结沉积物，m≈1.5适用于疏松沉积物[29]；如图7所示利用BGTL预测饱和水沉积地层段（0～90 mbsf）速度和实测纵波速度对比，其中使用P=8.0 MPa和CV=58%，改变m值预测速度，在高孔隙度低纵波速度时m=2.5预测速度拟合程度高，而在低孔隙度高声速时，m值应小于2.5，大于1。在可能含天然气水合物地层（104～160 mbsf）的中子孔隙度主要为45%～65%，因此，本次研究建模使用m=2.5。
使用上述参数，获得了井下剖面背景纵波速度和天然气水合物饱和度（图8，图9），如图所示，黑色实线为航次实测纵波速度，红色实线为本文使用BGTL模型计算结果，绿色区域为利用航次实测与模型速度计算的水合物饱和度，104～160 mbsf层段内测井实测VP大于理论饱和水背景VP，可能属于水合物储层区，因此导出水合物饱和度。结果显示，在104～160 mbsf的深度区间内平均饱和度约为6.0%，最高饱和度达到21.6%，其中130～145 mbsf层段内平均水合物饱和度为8.5%。
![]()
图 8 使用BGTL模型在U1517站位井测量的纵波速度和计算的基线速度比较Figure 8. Measured and calculated baseline P-wave velocities with BGTL at the Site U1517![]()
图 9 使用BGTL模型计算水合物储层区的背景纵波速度及饱和度Figure 9. Background P-wave velocities and gas hydrate saturations at Site U1517 with BGTL3. 讨论
图10为U1517站位井在104～160 mbsf层段的测井曲线，通过测井数据和背景基线看出存在3层纵波速度和电阻率明显异常的含水合物层段，同时，井径、密度和伽马测井数据均有不同程度的异常。在112～114 mbsf层段，环电阻率和声波速度明显增加，最高峰值分别为9.25 Ω·m和1.79 km/s，可能为含天然气水合物的薄层；在130～145 mbsf层段，环电阻率和声波速度最高峰值分别为2.88 Ω·m和1.90 km/s，属于较厚层的含水合物区域；在150～160 mbsf层段，密度与自然伽马降低较为明显，环电阻率和声波速度最高峰值分别为1.97 Ω·m和1.83 km/s。
![]()
图 10 U1517站位井井径、纵波速度、电阻率、密度和伽马测井曲线Vp为实测声波速度，Vpw为背景速度值；Rt-Ring为环电阻率；Rt-P40L为低频随钻相移电阻；R0为背景电阻率值。Figure 10. The well logs from site U1517A showing the caliper, P-wave velocity, resistivity, density and gamma rayVp is measured velocity, Vpw is calculated velocity baseline, Rt-Ring is ring resistivity, Rt-P40L is 400 kHz phasor resistivity, R0 is calculated resistivity baseline.如图11所示，由纵波速度数据通过STPE和BGTL模型估算了U1517站位井104～160 mbsf层段的天然气水合物饱和度，并与IODP372航次科学家利用阿尔奇公式和氯离子浓度两种方法计算结果相比较。STPE、BGTL和航次科学家利用阿尔奇公式3种模型在104～160 mbsf层段计算的平均饱和度分别为5.2%、6.0%和6.5%，130～145 mbsf层段的平均饱和度分别为7.9%、8.5%和9.6%；130～145 mbsf层段符合航次科学家使用氯离子浓度含量估算的高饱和度层段。112～114 mbsf和130～145 mbsf层段，阿尔奇公式估算的最高饱和度分别为56%和49%，大于BGTL和STPE计算结果，但是电阻率识别高饱和度薄层水合物约为2～5 cm，而声波测井分辨率约为15 cm，该薄层饱和度异常可能由于声波测井无法探测到而引起的，同时井径也发生变化，可能影响随钻测井速度与电阻率。氯离子异常在局部地层出现异常高值，从岩心分析看，异常高值与薄砂层相对应。因此，在104～160 mbsf层段3种方法估算的饱和度随深度变化相似，表明不同测井数据之间差异不大，且天然气水合物平均饱和度最高的层段为130～145 mbsf。使用STPE和BGTL模型计算的饱和水地层（0～90 mbsf）的纵波速度与实测纵波速度比较见图12所示，对于U1517井BGTL模型比STPE模型更适用于该站位水合物饱和度估算。
![]()
图 11 根据BGTL、STPE与电阻率、氯离子估算的天然气水合物饱和度的对比Figure 11. The gas hydrate saturation calculated by BGTL, STPE compared with the hydrate saturation calculated by the Expedition 372 scientists used the Archie equation and the chloride concentration![]()
图 12 STPE与BGTL在饱和水地层（0～90 mbsf）预测纵波速度与实测纵波速度对比Figure 12. The measured P-wave velocity compared with the calculated results of the STPE and BGTL in water-saturated sediments (0～90 mbsf)4. 结论
（1）通过U1517站位随钻测井和岩心数据综合分析，证实了该站位黏土质粉砂岩性不同层位存在天然气水合物，水合物呈层状分布。天然气水合物储层区域在104～160 mbsf层段，其中存在3层纵波速度和电阻率明显异常的含水合物层段（112～114、130～145和150～160 mbsf），112～114 mbsf层段可能为薄的天然气水合物层，而130～145 mbsf层段相较于其他层段水合物饱和度相对较高。其中112～114 mbsf层段天然气水合物饱和度最高，130～145 mbsf层段为主要天然气水合物赋存区域。
（2）依据LWD和取心数据，在计算过程中，根据岩性划分不同井段对应的矿物成分含量，用于纵波速度模型计算，并使用饱和水地层孔隙度与纵波速度拟合得到BGTL模型参数的方法，使BGTL模型更便于应用到测井资料估算水合物饱和度，通过STPE和BGTL模型计算出了U1517站位的水合物饱和度，并比较分析两种模型在饱和水地层的预测与实测纵波速度表明BGTL拟合度高于STPE；计算结果与航次科学家估算的饱和度相比，平均饱和度相近，3种方法计算的水合物饱和度值随深度变化相似，表明计算结果的合理性。
致谢：本研究所用样品和数据由IODP372航次提供，中国IODP办公室提供了胡高伟参加航次的旅费资助，在此一并致谢！
-
![]()
![]()
图 2 U1517站位测井数据
蓝色实线为实测数据，红色实线为拟合背景趋势，黑色虚线为BSR。
Figure 2. Well logs at Site U1517
The blue lines are measured logging data and the red lines are fitted background trends; the black dashed line is BSR.
![]()
图 3 U1517站位井地层因子与纵波速度交会图
黑色实线为0～90 mbsf拟合背景趋势线。
Figure 3. Cross plot of formation factor versus the measured P-wave velocities at Site U1517
The black line is fitted background trends of 0～90 mbsf.
![]()
图 4 U1517站位所取岩心的岩性和岩石矿物成分相对含量数据
黑线为依据岩性平均矿物成分相对含量。
Figure 4. Simplified lithostratigraphic column with bulk powder XRD results, Site U1517
Black line is average mineral composition based on core data.
![]()
图 5 使用STPE在U1517站位井测量的纵波速度和计算的基线速度的比较
Figure 5. Measured and calculated baseline P-wave velocities with STPE at the Site U1517
![]()
图 6 使用STPE计算水合物储层区的背景纵波速度及饱和度
Figure 6. Measured and calculated baseline resistivities and P-wave velocities at the Site U1517
![]()
图 8 使用BGTL模型在U1517站位井测量的纵波速度和计算的基线速度比较
Figure 8. Measured and calculated baseline P-wave velocities with BGTL at the Site U1517
![]()
图 9 使用BGTL模型计算水合物储层区的背景纵波速度及饱和度
Figure 9. Background P-wave velocities and gas hydrate saturations at Site U1517 with BGTL
![]()
图 10 U1517站位井井径、纵波速度、电阻率、密度和伽马测井曲线
Vp为实测声波速度，Vpw为背景速度值；Rt-Ring为环电阻率；Rt-P40L为低频随钻相移电阻；R0为背景电阻率值。
Figure 10. The well logs from site U1517A showing the caliper, P-wave velocity, resistivity, density and gamma ray
Vp is measured velocity, Vpw is calculated velocity baseline, Rt-Ring is ring resistivity, Rt-P40L is 400 kHz phasor resistivity, R0 is calculated resistivity baseline.
![]()
图 11 根据BGTL、STPE与电阻率、氯离子估算的天然气水合物饱和度的对比
Figure 11. The gas hydrate saturation calculated by BGTL, STPE compared with the hydrate saturation calculated by the Expedition 372 scientists used the Archie equation and the chloride concentration
![]()
图 12 STPE与BGTL在饱和水地层（0～90 mbsf）预测纵波速度与实测纵波速度对比
Figure 12. The measured P-wave velocity compared with the calculated results of the STPE and BGTL in water-saturated sediments (0～90 mbsf)
-
[1] Nixon M F, Grozic J L H. Submarine slope failure due to gas hydrate dissociation: a preliminary quantification [J]. Canadian Geotechnical Journal, 2007, 44(3): 314-325. doi: 10.1139/t06-121
[2] Brewer P G. Gas hydrates and global climate change [J]. Annals of the New York Academy of Sciences, 2000, 912(1): 195-199.
[3] Wallace L M, Webb S C, Ito Y, et al. Slow slip near the trench at the Hikurangi subduction zone, New Zealand [J]. Science, 2016, 352(6286): 701-704. doi: 10.1126/science.aaf2349
[4] Mountjoy J J, Pecher I, Henrys S, et al. Shallow methane hydrate system controls ongoing, downslope sediment transport in a low-velocity active submarine landslide complex, Hikurangi Margin, New Zealand [J]. Geochemistry, Geophysics, Geosystems, 2014, 15(11): 4137-4156. doi: 10.1002/2014GC005379
[5] Singh C. Multiphase fluid flow in presence of Transverse Magnetic Field [J]. Eos Transactions American Geophysical Union, 1987(53).
[6] Townend J. Estimates of conductive heat flow through bottom-simulating reflectors on the Hikurangi and southwest Fiordland continental margins, New Zealand [J]. Marine Geology, 1997, 141(1-4): 209-220. doi: 10.1016/S0025-3227(97)00073-X
[7] Henrys S A, Ellis S, Uruski C. Conductive heat flow variations from bottom-simulating reflectors on the Hikurangi margin, New Zealand [J]. Geophysical Research Letters, 2003, 30(2): 1065.
[8] Pecher I A, Henrys S A, Zhu H. Seismic images of gas conduits beneath vents and gas hydrates on Ritchie Ridge, Hikurangi margin, New Zealand [J]. New Zealand Journal of Geology and Geophysics, 2004, 47(2): 275-279. doi: 10.1080/00288306.2004.9515054
[9] Crutchley G J, Pecher I A, Gorman A R, et al. Seismic imaging of gas conduits beneath seafloor seep sites in a shallow marine gas hydrate province, Hikurangi Margin, New Zealand [J]. Marine Geology, 2010, 272(1-4): 114-126. doi: 10.1016/j.margeo.2009.03.007
[10] Navalpakam R S, Pecher I A, Stern T. Weak and segmented bottom simulating reflections on the Hikurangi Margin, New Zealand—Implications for gas hydrate reservoir rocks [J]. Journal of Petroleum Science & Engineering, 2012, 88-89: 29-40.
[11] Gorman A R, Fletcher P T, Baker D, et al. Characterisation of focused gas hydrate accumulations from the Pegasus Basin, New Zealand, using high-resolution and conventional seismic data [J]. ASEG Extended Abstracts, 2018, 2018(1): 1-4.
[12] Schwalenberg K, Haeckel M, Poort J, et al. Evaluation of gas hydrate deposits in an active seep area using marine controlled source electromagnetics: Results from Opouawe Bank, Hikurangi Margin, New Zealand [J]. Marine Geology, 2010, 272(1-4): 79-88. doi: 10.1016/j.margeo.2009.07.006
[13] Schwalenberg K, Wood W, Pecher I, et al. Preliminary interpretation of electromagnetic, heat flow, seismic, and geochemical data for gas hydrate distribution across the Porangahau Ridge, New Zealand [J]. Marine Geology, 2010, 272(1-4): 89-98. doi: 10.1016/j.margeo.2009.10.024
[14] Schwalenberg K, Rippe D, Koch S, et al. Marine-controlled source electromagnetic study of methane seeps and gas hydrates at Opouawe Bank, Hikurangi Margin, New Zealand [J]. Journal of Geophysical Research: Solid Earth, 2017, 122(5): 3334-3350. doi: 10.1002/2016JB013702
[15] Faure K, Greinert J, Pecher I A, et al. Methane seepage and its relation to slumping and gas hydrate at the Hikurangi margin, New Zealand [J]. New Zealand Journal of Geology and Geophysics, 2006, 49(4): 503-516. doi: 10.1080/00288306.2006.9515184
[16] Naudts L, Greinert J, Poort J, et al. Active venting sites on the gas-hydrate-bearing Hikurangi Margin, off New Zealand: Diffusive-versus bubble-released methane [J]. Marine Geology, 2010, 272(1-4): 233-250. doi: 10.1016/j.margeo.2009.08.002
[17] Greinert J, Lewis K B, Bialas J, et al. Methane seepage along the Hikurangi Margin, New Zealand: Overview of studies in 2006 and 2007 and new evidence from visual, bathymetric and hydroacoustic investigations [J]. Marine Geology, 2010, 272(1-4): 6-25. doi: 10.1016/j.margeo.2010.01.017
[18] Fraser D R A, Gorman A R, Pecher I A, et al. Gas hydrate accumulations related to focused fluid flow in the Pegasus Basin, southern Hikurangi Margin, New Zealand [J]. Marine and Petroleum Geology, 2016, 77: 399-408. doi: 10.1016/j.marpetgeo.2016.06.025
[19] Crutchley G J, Kroeger K F, Pecher I A, et al. Gas hydrate formation amid submarine canyon incision: Investigations from New Zealand’s Hikurangi subduction margin [J]. Geochemistry, Geophysics, Geosystems, 2017, 18(12): 4299-4316. doi: 10.1002/2017GC007021
[20] Crutchley G J, Kroeger K F, Pecher I A, et al. How tectonic folding influences gas hydrate formation: New Zealand's Hikurangi subduction margin [J]. Geology, 2019, 47(1): 39-42. doi: 10.1130/G45151.1
[21] Gross F, Mountjoy J J, Crutchley G J, et al. Free gas distribution and basal shear zone development in a subaqueous landslide-Insight from 3D seismic imaging of the Tuaheni Landslide Complex, New Zealand [J]. Earth and Planetary Science Letters, 2018, 502: 231-243. doi: 10.1016/j.jpgl.2018.09.002
[22] 胡高伟, 李承峰, 业渝光, 等. 沉积物孔隙空间天然气水合物微观分布观测[J]. 地球物理学报, 2014, 57(5):1675-1682. [HU Gaowei, LI Chengfeng, YE Yuguang, et al. Observation of gas hydrate distribution in sediment pore space [J]. Chinese Journal of Geophysics, 2014, 57(5): 1675-1682. doi: 10.6038/cjg20140530 [23] Hu G W, Ye Y G, Zhang J, et al. Acoustic response of gas hydrate formation in sediments from South China Sea [J]. Marine and Petroleum Geology, 2014, 52: 1-8. doi: 10.1016/j.marpetgeo.2014.01.007
[24] Bu Q T, Hu G W, Ye Y G, et al. The elastic wave velocity response of methane gas hydrate formation in vertical gas migration systems [J]. Journal of Geophysics and Engineering, 2017, 14(3): 555-569. doi: 10.1088/1742-2140/aa6493
[25] Archie G E. The electrical resistivity log as an aid in determining some reservoir characteristics [J]. Transactions of AIME, 1942, 146(1): 54-62. doi: 10.2118/942054-G
[26] Montaron B. Connectivity theory-a new approach to modeling non-Archie rocks [J]. Petrophysics, 2009, 50(2): 102-115.
[27] Lee M W, Hutchinson D R, Collett T S, et al. Seismic velocities for hydrate-bearing sediments using weighted equation [J]. Journal of Geophysical Research: Solid Earth, 1996, 101(B9): 20347-20358. doi: 10.1029/96JB01886
[28] 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(13): 2021-2024. doi: 10.1029/1999GL900421
[29] Lee M W. Modified Biot-Gassmann theory for calculating elastic velocities for unconsolidated and consolidated sediments [J]. Marine Geophysical Researches, 2002, 23(5-6): 403-412.
[30] Lee M W, Waite W F. Estimating pore-space gas hydrate saturations from well log acoustic data [J]. Geochemistry, Geophysics, Geosystems, 2008, 9(7): Q07008.
[31] Lee M W. Models for gas hydrate-bearing sediments inferred from hydraulic permeability and elastic velocities[R]. Scientific Investigations Report 2008-5219. U. S. Geological Survey, 2008.
[32] 王秀娟, 钱进, Lee M. 天然气水合物和游离气饱和度评价方法及其在南海北部的应用[J]. 海洋地质与第四纪地质, 2017, 37(5):35-47. [WANG Xiujuan, QIAN Jin, Lee M. Methods for estimation of gas hydrate and free gas saturations and application to the northern slope of South China Sea [J]. Marine Geology and Quaternary Geology, 2017, 37(5): 35-47. [33] Lee M W, Collett T S. Scale-dependent gas hydrate saturation estimates in sand reservoirs in the Ulleung Basin, East Sea of Korea [J]. Marine and Petroleum Geology, 2013, 47: 195-203. doi: 10.1016/j.marpetgeo.2012.09.004
[34] Wang X J, Lee M, Collett T, et al. Gas hydrate identified in sand-rich inferred sedimentary section using downhole logging and seismic data in Shenhu area, South China Sea [J]. Marine and Petroleum Geology, 2014, 51: 298-306. doi: 10.1016/j.marpetgeo.2014.01.002
[35] Hu G W, Ye Y G, Zhang J, et al. Acoustic properties of gas hydrate-bearing consolidated sediments and experimental testing of elastic velocity models [J]. Journal of Geophysical Research, 2010, 115(B2): B02102.
[36] Pecher I A, Barnes P M, LeVay L J, et al. International ocean discovery program expedition 372 preliminary report creeping gas hydrate slides and Hikurangi LWD [J]. Integrated Ocean Drilling Program: Preliminary Reports, 2018, 372: 1-35.
[37] Pecher I A, Barnes P M, LeVay L J, et al. Proceedings of the International Ocean Discovery Program, 372A[Z]. College Station, TX: International Ocean Discovery Program, 2019.
[38] Waite W F, Santamarina J C, Cortes D D, et al. Physical properties of hydrate-bearing sediments [J]. Reviews of Geophysics, 2009, 47(4): RG4003.
[39] Mavko G, Mukerji T, Dvorkin J. The Rock Physics Handbook[M]. 2nd ed. Cambridge: Cambridge University Press, 2009.
[40] Lee M W. Proposed moduli of dry rock and their application to predicting elastic velocities of sandstones[R]. U. S. Geological Survey, 2005.
[41] Mindlin R D. Compliance of elastic bodies in contact [J]. Journal of Applied Mechanics, 1949, 16(3): 259-268.
[42] Han D H, Nur A, Morgan D. Effects of porosity and clay content on wave velocities in sandstones [J]. Geophysics, 1986, 51(11): 2093-2107. doi: 10.1190/1.1442062
[43] Tiab D, Donaldson E C. Petrophysics: Theory and Practice of Measuring Reservoir Rock and Fluid Transport Properties[M]. 3rd ed. Houston: Gulf Professional Publishing, 2012.
-
期刊类型引用(2)
1. 聂鑫，杜文波，高红芳，胡小三，杨楚鹏，鞠东. 南海北部神狐峡谷区沉积相特征及BSR地质灾害指示. 地质学报. 2024(09): 2737-2752 . 
百度学术
2. 龚智，王秀娟，周吉林，邓炜，匡增桂，苏丕波，闫伟超，李三忠. 基于随钻电阻率测井的天然气水合物赋存类型的识别方法. 地质学报. 2024(09): 2709-2722 . 百度学术
其他类型引用(3)
下载:
计量
- 文章访问数: 2624
- HTML全文浏览量: 593
- PDF下载量: 51
- 被引次数: 5
