南海北部神狐海域天然气水合物分解的测井异常

南海北部神狐海域GMGS-1钻探揭示SH3井天然气水合物位于稳定带上部,厚度约为10 m.氯离子异常计算的水合物饱和度最高达26%,高水合物饱和度层出现高电阻率和低纵波速度.为分析该低纵波速度异常,本文基于简化的三相介质理论计算了饱和水纵波速度,在深度195 m附近,测量的纵波速度小于饱和水纵波速度.利用阿尔奇公式,基于原位温度、盐度、密度孔隙度和测量的电阻率,利用交会分析确定了该井的阿尔奇常数为a=1.1和m=2.3.基于该参数,利用阿尔奇方程计算的水合物饱和度占孔隙空间5%~20%,局部地层水合物饱和度达26.8%,在垂向上分布不均匀.由于钻探可能导致水合物发生分解而产生游离气,原位游离气和水合物分解产生的气体都能造成低纵波速度异常.由于地震资料采集在测井之前完成,利用不同速度制作合成地震记录并与地震资料进行对比,能够确定水合物稳定带上部的低速异常形成原因.     

天然气水合物勘探开发新技术进展

天然气水合物是一种非常规能源,因其巨大的资源量成为目前能源研究的热点.水合物的研究目前仍处于初期阶段,人们在水合物勘探、开发模拟方面做了大量尝试,总结了许多有利于水合物勘探和开发的关键技术.本文首先介绍了天然气水合物形成的地质条件和储层特征,然后介绍了针对海底天然气水合物的地震勘探和电磁、地化等非地震勘探最新技术进展,分析了天然气水合物矿藏的测井响应特征,结合岩石物理模型对储层孔隙度、饱和度进行估算用于储量预测.最后,论述了天然气水合物开发过程中的深水区的钻井问题,水合物的开采方式总结以及商业化生产所面临的挑战.     

南海北部神狐海域天然气水合物成藏模式研究

南海北部陆坡神狐海域是我国海洋天然气水合物勘探开发研究的重点靶区,独特的水合物成藏特征,难以利用当前观测到的沉积速率和流体流动条件对其成藏机理进行解释和量化说明,对其形成演化模式和控制因素尚不明确.本文构建了海洋天然气水合物形成演化过程的动力学模型,模型的主控参量为海底沉积速率和水流速率,以此计算了神狐海域天然气水合物聚集演化过程,并与饱和度的盐度测试值进行对比.最后,在研究神狐海域地质构造活动和水合物成藏动力学基础上建立了神狐天然气水合物形成演化模式.认为神狐海域当前的天然气水合物是在上新世末—更新世早期断裂体系水合物基础上继承演化而来的,神狐海域天然气水合物形成演化具有典型的二元模式.第一阶段水合物形成发生在距今1.5 Ma之前构造活动形成的断裂体系中,高达50 m/ka的孔隙水流动携带了大量的甲烷进入水合物稳定带,导致了水合物的快速生成,在4万年内形成了饱和度达20%的甲烷水合物;第二阶段发生在1.5 Ma以来,泥质粉砂沉积使沉积体渗透率骤减,0.7 m/ka的低速率水流使甲烷供给不足,在海底浅层新沉积体中无法生成水合物,仅在水合物稳定带底部有缓慢的水合物继承增长,并因此形成了神狐海域当前观测到的水合物产出特征,而且水合物资源量仍在减少.     

天然气水合物和游离气饱和度估算的影响因素

讨论了不同水合物胶结类型的流体饱和多孔隙固体中地震波的衰减情况,分析了估算天然气水合物和游离气饱和度影响因素.结果表明, 地层孔隙度、纵波速度模型和弹性模量的计算方法是影响反演水合物和游离气饱和度的关键因素.含水合物地层的吸收与水合物胶结类型密切相关,当水合物远离固体颗粒,像流体一样充填在孔隙时,品质因子出现负异常,而当水合物胶结固体颗粒影响骨架的弹性性质,其品质因子出现正异常.根据布莱克海台地区164航次995井的测井资料,分别应用低频和高频速度模型估算了水合物和游离气饱和度.由低频速度模型得到的水合物饱和度(占孔隙空间的)10%～20%,游离气饱和度(占孔隙空间的)05%～1%;而由高频速度模型得到的水合物饱和度(占孔隙空间的)5%～10%,游离气饱和度(占孔隙空间的)1%～2%.     

天然气水合物饱和度与声学参数响应关系的实验研究

天然气水合物饱和度是评估天然气水合物资源量的重要参数,而用来估算饱和度的速度模型则是有限的几个,它们有的为经验公式,有的是以实验数据或野外资料为基础建立起来的,需要进行实验验证和参数分析以确定其使用的适用性.本文首次利用超声和时域反射联合探测技术,研究了沉积物中水合物饱和度与声学特性的关系,并对时间平均方程、伍德及其修正方程、李权重方程和BGTL（Biot-Gassmann Theory by Lee）理论等常用的水合物饱和度估算模型进行了验证.实验结果表明,超声和时域反射联合探测技术能有效地实时获得水合物饱和度和纵、横波速的实验数据,李权重方程和BGTL理论的速度预测值与实验值比较吻合,有广泛的适用性.     

含裂隙储层中的天然气水合物分布:来自地震各向异性测量的启示

不同的天然气水合物分布对储层地震性质的影响存在显著差异,深入认识储层中的水合物分布,对基于地震性质定量评价水合物资源具有重要意义.然而,目前对含裂隙储层中水合物分布的认识仍然严重不足.为此,文章设计并开展了针对性岩石物理实验,测量了含硬币状定向裂隙岩石在水合物生成过程中的各向异性速度,并进一步分析了不同水合物饱和度条件下裂隙中的水合物分布,及其对各向异性速度和各向异性参数的影响.实验结果表明:随着水合物饱和度的增加,各向异性速度均呈现出整体增加的趋势;当水合物饱和度大于10%时,各向异性速度的增加幅度变得更大,且各向异性速度间的增加速率也表现出强烈差异,具体表现为,垂直于层理面与裂隙面传播的纵波速度和偏振方向垂直于层理面的横波速度分别比其他纵波和横波速度增加得更快.对各向异性结果的深入解释表明:在水合物饱和度小于10%时,水合物可能倾向于黏结在裂隙表面的矿物颗粒上,而当水合物饱和度大于10%时,以黏结矿物颗粒分布的水合物演化为连接裂隙上下表面的桥接分布.研究结果提出了对含裂隙储层水合物分布及其对各向异性地震性质影响的新认识,为含裂隙水合物储层的定量评价提供了重要依据.     

裂隙充填型天然气水合物的各向异性分析及饱和度估算——以印度东海岸NGHP01-10D井为例

沿裂隙发育的天然气水合物是印度深水盆地细粒沉积物中水合物的重要产出方式,水合物以结核状或脉状充填在高角度裂隙中.天然气水合物主要沿着构造主应力方向生成,由于裂隙的存在,含水合物的沉积物层呈现各向异性.利用孔隙介质中水合物呈均匀分布的速度模型计算的NGHP01-10D井水合物饱和度高达40%,而压力取芯表明水合物饱和度占孔隙空间的20%左右.为了研究水合物饱和度差异,基于层状介质的各向异性模型计算了裂隙充填型水合物的饱和度.在垂直井孔中,由于波入射角与裂隙倾角有关,考虑裂隙倾角变化,利用纵波和横波速度同时反演水合物饱和度和裂隙倾角.利用层状介质模型计算的水合物占孔隙空间的15%～25%,裂隙的倾角在60°～90°,多为高角度裂隙.在NGHP01-10D井中,纵横波速度联合计算的饱和度与压力取芯结果吻合更好.     

南海北部神狐海域天然气水合物成藏演化分析研究

南海北部陆坡具备天然气水合物成藏的基本地质条件,神狐海域天然气水合物是当前我国海洋天然气水合物勘探开发研究的重点靶区.然而,神狐海域水合物集中分布在水合物稳定带的底部薄层中,饱和度高,其水合物特征与典型的低甲烷通量控制的水合物分布有很大差异,对其成藏机理和控制因素尚不明确.本文构建了针对神狐水合物成藏过程的一维动力学模型,模型包括沉积压实作用、甲烷溶解度、以及水合物生成和沉积体渗透率,模拟计算的主控参量为海底沉积速率和水流通量,在孔隙水中甲烷浓度一定的情况下,水流速率决定了溶解甲烷的迁移速率和稳定带中甲烷的供给速率,并以此模型计算了神狐海域水合物聚集成藏的动力学过程.模型讨论了特定沉积速率和水流通量条件下水合物成藏与分布特征,并与实际观测数据进行比较研究.研究发现,基于当前沉积速率和水流通量条件模拟的水合物形成演化过程,与神狐海域实际水合物分布特征存在很大差异,但在假定系统中水合物饱和度初值达16~20%时,模拟的水合物饱和度分布特征与观测数据吻合,并因此推测在早期地质历史上,神狐海域存在更加丰富的甲烷水合物,当前的水合物分布特征是在对早期水合物继承基础上发展而成的,而且神狐海域水合物含量正逐渐减少.     

南海北部无明显BSR地区天然气水合物识别研究

似海底反射(BSR)作为海域天然气水合物的重要地震识别标志之一,已得到广泛认同.然而,科学钻探证实,天然气水合物和BSR之间并不具有严格对应关系,即存在水合物的地区却并不一定存在BSR.本文在分析BSR与天然气水合物非充分必要对应关系及其原因的基础上,着重对无明显BSR地区天然气水合物的地震识别方法和应用进行研究.对经钻探证实的存在天然气水合物的神狐海域地震资料进行处理分析后发现,含水合物沉积层具有层速度相对较高、高波阻抗、瞬时高频等特征,且层速度反演、波阻抗反演和瞬时属性分析等方法能有效识别无明显BSR区域的天然气水合物.最后,综合利用这些识别方法,应用于琼东南盆地无明显BSR地区的天然气水合物地震识别,取得了较好的效果.     

南中国海神狐海域天然气水合物 地震识别及分布特征

南中国海北部陆坡神狐海域是天然气水合物发育的有利地区,钻井取芯表明此区域的天然气水合物以细颗粒状分布于水合物稳定带内.以三维地震数据为基础,讨论海洋天然气水合物调查的高分辨率地震成像关键技术,综合利用不同的地球物理方法(阻抗反演、属性聚类分析及神经网络方法等),以钻井结果为约束分析不同方法对天然气水合物识别的适用性,确定利用属性聚类方法来获得神狐海域天然气水合物与游离气的三维空间展布特征.将钻井取芯结果与地质条件相结合,综合多方面地球物理信息,分析产生这一分布特征的地质原因,指出断层与裂缝是南海北部陆坡神狐海域水合物、游离气的分布模式及天然气水合物形成的主要地质控制因素.     

Magnetotelluric data inversion with seismic data constraint

Introduction The research on the structure and physical property of ancient hidden hill, igneous rocks and basement is relatively difficult by using seismic data only. If we combine seismic data, magneto-telluric (MT) data and geophysical data together, better results can be obtained for the above problem. A number of geophysicists at home and abroad, such as CHEN and WANG (1990), Siri-punvarapor and Egbert (2000) have tried many methods to solve the problem by the inversion of seismic da…     

数据自相关多次波偏移成像

在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义.     

Super-resolution data assimilation

Ocean Dynamics - Increasing model resolution can improve the performance of a data assimilation system because it reduces model error, the system can more optimally use high-resolution...     

Nonlinear data processing method for the signal enhancement of GPR data

An alternative data processing procedure is proposed in this paper for the purpose of enhancing the signal/noise (S/N) ratio of ground penetrating radar (GPR) data. The processing methodology is achieved by performing the logarithmic transform in conjunction with the ensemble empirical mode decomposition (EEMD), a new nonlinear data analysis method in signal processing. The synthetic model study and field example indicate that the logarithmic transform is effective in alleviating the attenuation problem. Additionally, the spectrogram obtained from Hilbert-Huang transform (HHT) shows that the decomposition sensitivity of the EEMD method is greatly improved with the aid of the logarithmic transform. This new method allows us to extract the signal components from noisy GPR data efficiently. The success of this study suggests a possible nonlinear analysis application in future GPR investigation, particularly in the filter design and gain correction.     

An automatic method for data processing of seismic data in tunneling

We present a new method for the prediction of the discontinuities and lithological variations ahead of the tunnel face. The automatic procedure is applied to data collected by seismic reflection surveys, with the sources and sensors located along the tunnel. The method allows: i) to estimate an average value of the wave velocity; ii) to detect the discontinuities for each source point; and iii) to analyze and plot the number of superposing estimates for each node of the domain. The final result can be interpreted as the probability to detect a discontinuity at a certain distance from the tunnel face. The method automatically estimates the peaks in the seismograms that can be related to a reflection. On the base of this process, the method only requires the source–receiver geometry and the data acquisition parameters. The procedure has been tested on synthetic and real data coming from a seismic survey on a tunnel under construction. The results indicate that the method runs very fast and it is reliable in the identification of lithological changes and discontinuities, up to more than 100 m ahead of the tunnel face.     

基于大数据技术的地震数据处理新思路

现代地震观测积累了大量、种类繁多的地震前兆数据和震后数据,但数据中隐含的地震规律及趋势无法用经典的算例、数学公式、物理公式进行定量解释、分析和预测,需要探索新方法、新技术。案例表明,大数据分析具有预测事物发展趋势、改变传统观念和发现新事物的功能,有助于从错误信息中挖掘有价值信息。本文由此提出基于大数据技术处理地震数据的新思路,对地震监测预报新思路进行探索;通过在大数据平台上基于地震目录的余震预测研究,给出大数据技术在地震数据处理方面的一个应用实例,验证该思路的可行性。