基于深水区宽频地震的天然气水合物识别方法

仅利用地震似海底反射（BSR）识别琼东南盆地深水区天然气水合物存在一定的局限性,从而影响天然气水合物的勘探成效。笔者利用天然气水合物已钻井数据,分析该盆地深水区天然气水合物岩石弹性参数特征,用以查明天然气水合物的岩石物理规律;同时,利用地震正演模拟,明确了研究区发育的孔隙型、烟囱型水合物的地震反射特征。在此基础上,利用AVO正演判识真假BSR：天然气水合物底界面反射具有Ⅲ类AVO且存在AVO异常,此为真BSR反射;而块体流（MTD）底界面虽类似BSR反射,但其AVO为Ⅳ类且AVO无异常特征。利用宽频地震数据和三维地震速度体进行速度模型下的宽频确定性反演,并通过高速异常、高阻抗异常描述天然气水合物发育情况。总之,利用地震反射特征、AVO特征、无井宽频地震反演等手段,实现了琼东南盆地深水区多种类型天然气水合物的地震识别,判识圈定了水合物矿藏。     

东海天然气水合物的地震特征

使用中国科学院海洋研究所"科学一号"调查船于2001年以及20世纪80年代在东海地区采集的多道地震资料,以海域天然气水合物研究为目的,对这些资料进行了数据处理并获得了偏移地震剖面。通过对地震剖面的解释,在6条剖面上确定了6段异常反射为BSR,均有振幅强、与海底相位相反的特点。6段BSR基本上都没有出现和沉积地层相交的现象。分析认为,这与东海地区第四纪以来的沉积特征有关,并不能由此否认这些异常反射是BSR。6段BSR出现的水深为750~2 000 m,埋深在0.1~0.5 s (双程时间)之间。随着海底深度的增大,BSR埋深有增大的趋势。计算结果显示,6段BSR所处的温度和压力条件都满足水合物稳定赋存所需要的温度和压力条件。本文的BSR主要与北卡斯凯迪亚盆地以及智利海域水合物的温度、压力条件相似,而与日本南海海槽、美国布莱克海台等海域水合物的温度、压力条件相差比较大。在地震剖面上,6段BSR所处的局部构造位置都和挤压、断层有关,有利于水合物的发育;在空间上,它们主要分布在东海陆坡近槽底的位置以及与陆坡相近的槽底。在南北方向上,除分布在吐噶喇断裂和宫古断裂附近外,还与南奄西、伊平屋和八重山热液活动区相邻。热液活动和水合物虽然没有直接的成因关系,但岩浆活动为水合物气源的形成提供了热源条件,为流体和气体的运移、聚集提供了通道条件,从而有利于水合物的发育与赋存。根据地震剖面反射特征推断,剖面A1A2和A14A23发育BSR的位置应该有气体或者流体从海底流出,可能是海底冷泉发育的位置。剖面A14A23上BSR发育处,振幅比的异常增大和BSR埋深的降低是相关联的。这种关联支持该处发育海底冷泉的推测。     

南海北部琼东南海域气烟囱发育特征及其对水合物形成与分布的影响

琼东南海域地震资料解释发现了BSR(似海底反射)、BSR下伏强反射及烟囱体等天然气水合物的地震响应特征。研究发现,区内气烟囱的分布与BSR的分布存在明显的相关性,气烟囱是气体垂直运移的主要通道,气体垂直向上运移至水合物稳定带大量聚集,从而形成水合物。因此,精细刻画研究区气烟囱发育特征对于区内天然气水合物成藏及分布研究具有重要意义。传统气烟囱识别方法只通过地震剖面上的弱反射或相关属性分析,笔者利用基于MLP算法的神经网络,高效地分析了本区域气烟囱的分布,并根据烟囱与BSR分布之间的关系,分析了气烟囱对天然气水合物形成及分布的影响。     

用于研究东海天然气水合物的地震资料处理方法

通过对东海海域二维地震、单道地震、浅层剖面等资料进行的综合研究表明 :用于研究天然气水合物的地震资料处理应提高速度分析精度和分辨能力 ,进行子波估算 ,压制多次波 ,相对保持振幅 ,DMO,AVO及波阻抗特殊处理等。在地震剖面上天然气水合物主要特征有 :BSR、振幅异常、速度异常、AVO异常等标志特征。据此 ,可对天然气水合物进行识别和预测。东海海域是天然气水合物可能赋存的有利部位 ,其中冲绳海槽是天然气水合物成藏的目标区域。     

基于声波阻抗估算印度安达曼海天然气水合物饱和度

印度国家天然气水合物01航次NGHP01-17井测井资料表明,安达曼海的薄层火山灰中存在天然气水合物,饱和度高达孔隙空间的20%。BSR上厚度为20 m地层的水合物饱和度约为10%。利用稀疏脉冲反演了含水合物层的声波阻抗剖面,基于测井获得声波阻抗与饱和水孔隙度之间的关系,计算了地层饱和水孔隙度剖面,结合阿尔奇公式估算了地震剖面的水合物饱和度。从地震数据计算的天然气水合物饱和度占孔隙空间的6%~15%,局部饱和度高达20%,在空间分布上具有不均匀性。     

海洋天然气水合物的地球物理勘探技术

天然气水合物将成为21世纪的替代能源，地球物理方法是勘探天然气水合物的重要手段。本文比较全面地分析和总结了天然气水合物的各种地球物理识别技术以及地震资料的特殊处理和分析方法，详细地介绍了水平地震剖面、垂直地震剖面、测井以及旁侧声纳剖面上天然气水合物的表现和识别方法。特别地，针对海洋地震资料的特点以及天然气水合物在地震剖面上的识别标志BSR、振幅空白带等特征，文章引入了真振幅处理、子波处理以及多项式拟合等处理方法来提高天然气水合物识别标志在地震剖面上的显示效果。最后，为了全面了解海底天然气水合物的分布 以及微细结构，文章介绍了AVO分析、全波形反演速率分析、叠加速度分析和走时反演等正、反演技术。     

南海北部揭阳凹陷天然气水合物的地震异常特征分析

大量研究表明南海北部珠江口盆地是天然气水合物发育区,但是该盆地东部揭阳凹陷水合物研究较少。本文利用揭阳凹陷新采集三维地震资料,对该三维地震资料进行成像道集优化和叠前时间偏移处理,得到针对水合物的新处理地震数据体,并通过高精度网格层析反演得到层速度数据体。利用该数据开展叠后约束稀疏脉冲反演,获得含天然气水合物地层波阻抗异常,综合分析反演与地震属性识别水合物。从新处理地震资料看,该区域似海底反射(bottom simulation reflection,BSR)反射呈连续、不连续与地层斜交等特征,BSR发育在一个继承性小型水道上,且下部断裂和气烟囱发育。通过分析BSR特征及BSR上下地层的速度、波阻抗、振幅、频率、相干等属性异常,结合水合物成藏条件,发现了南海北部新的天然气水合物有利富集区,为该区域水合物勘探提供基础。     

利用地震谱反演技术圈定南海神狐海域天然气水合物储集层

地震谱反演是近年发展起来的一种新的地震方法,可利用地震数据构建反射系数剖面,极大地提高地震资料的分辨率。利用谱反演方法对南海神狐海域天然气水合物赋存区的地震资料进行了处理,获得该区海底地层的反射系数剖面,结合BSR和测井资料,确定了该区天然气水合物储集层的顶界和底界,从而圈定了该区天然气水合物储集层。     

地震多参数约束下的真假BSR识别

地震勘探是探测海底天然气水合物的重要手段,利用地震资料的诸多特征可以较好地识别海底天然气水合物,尤其是在识别似海底反射(BSR)方面发挥着重要作用。由于多次波等特征与BSR有很多相似之处,如果辨别不当就很容易被误认为是BSR,将会得出错误的结论。以我国某海域实际资料为例,从研究BSR的地震特征出发,指明多次波、气泡效应等多种假BSR现象,提出了利用精细速度分析、AVO特征分析、多次波压制等多种地震参数约束以识别真假BSR,进而提高海洋地震勘探精度,为寻找更多的海底天然气水合物提供技术保障。     

台湾增生楔天然气水合物的地震特征

通过对南海973航次在该区域的多道地震剖面的分析,结合该区域的地质背景,认为台湾增生楔具有天然气水合物存在和分布的地球物理特征,在地震剖面上出观海底反射层(BSR)、振幅空白(BZ)、极性反转等地震识别标志。BSR所在区域位于南海向菲律宾海板块俯冲的增生楔上,南海新生代沉积不仅提供了富含有机质的丰富物源,而且类似于活动大陆边缘的构造体系又为天然气水合物的形成提供了良好的条件。     

冲绳海槽天然气水合物稳定带特征及资源量评价

根据冲绳海槽多道地震资料的处理解释,在16条地震剖面上发现了水合物似海底反射层BSR,经过AVO、波形反演等特殊的处理技术,首次直接利用BSR圈定了冲绳海槽天然气水合物的具体分布范围,直接利用数据得出了天然气水合物稳定带厚度在冲绳海槽的分布趋势,认为海槽南部最厚,中部次之,北部最薄,并通过计算得出了冲绳海槽水合物稳定带的厚度和水合物资源量,对今后海槽水合物勘查和资源量评价具有一定的指导意义.     

Source and accumulation of gas hydrate in the northern margin of the South China Sea

The northern South China Sea (NSCS) experienced continuous evolution from an active continental margin in the late Mesozoic to a stable passive continental margin in the Cenozoic. It is generally believed that the basins in the NSCS evolved as a result of Paleocene–Oligocene crustal extension and associated rifting processes. This type of sedimentary environment provides a highly favourable prerequisite for formation of large-scale oil- and gas–fields as well as gas hydrate accumulation. Based on numerous collected data, combined with the tectonic and sedimentary evolution, a preliminary summary is that primitive coal-derived gas and reworked deep gas provided an ample gas source for thermogenic gas hydrate, but the gas source in the superficial layers is derived from humic genesis. In recent years, the exploration and development of the NSCS oil, gas and gas hydrate region has provided a basis for further study. A number of 2D and 3D seismic profiles, the synthetic comparison among bottom simulating reflector (BSR) coverage characteristics, the oil-gas area, the gas maturity and the favourable hydrate-related active structural zones have provided opportunities to study more closely the accumulation and distribution of gas hydrate. The BSR has a high amplitude, with high amplitude reflections below it, which is associated with gas chimneys and pockmarks. The high amplitude reflections immediately beneath the BSR are interpreted to indicate the presence of free gas and gas hydrate. The geological and geochemical data reveal that the Cenozoic northern margin of the NSCS has developed coal-derived gas which forms an abundant supply of thermogenic gas hydrate. Deep-seated faults and active tectonic structures facilitate the gas migration and release. The thermogenic gas hydrate and biogenic gas are located at different depths, have a different gas source genesis and should be separately exploited. Based on the proven gas hydrate distribution zone, we have encircled and predicted the potential hydrate zones. Finally, we propose a simple model for the gas hydrate accumulation system in the NSCS Basin.     

天然气水合物BSR的识别与地震勘探频率

地震勘探是调查天然气水合物广泛使用而有效的方法,而BSR是水合物赋存的主要标志.通过对实际调查资料的分析对比,结合国外的调查研究成果,探讨了地震勘探频率在BSR识别中的影响和作用,提出了在我国海洋天然气水合物的地震调查中有利于BSR识别的合适的频率范围.     

The enigma of double BSRs: indicators for changes in the hydrate stability field?

 A classical bottom simulating reflector (BSR) and a presently unknown double BSR pattern are detectable in reflection seismic profiles from the Storegga Slide area west of Norway. Pressure and temperature modeling schemes lead to the assumption that the strong BSR marks the base of a hydrate stability zone with a typical methane gas composition of 99%. The upper double BSR may mark the top of gas hydrates and the lower double BSR may represent a relict of former changes of the hydrate stability field from glacial to interglacial times or the base of gas hydrates with a gas composition including heavier hydrocarbons.     

Mass-transport deposits and gas hydrate occurrences in the Ulleung Basin,East Sea – Part 2: Gas hydrate content and fracture-induced anisotropy

Mass-transport-deposits (MTDs) and hemipelagic mud interbedded with sandy turbidites are the main sedimentary facies in the Ulleung Basin, East Sea, offshore Korea. The MTDs show similar seismic reflection characteristics to gas-hydrate-bearing sediments such as regional seismic blanking (absence of internal reflectivity) and a polarity reversed base-reflection identical to the bottom-simulating reflector (BSR). Drilling in 2007 in the Ulleung Basin recovered sediments within the MTDs that exhibit elevated electrical resistivity and P-wave velocity, similar to gas hydrate-bearing sediments. In contrast, hemipelagic mud intercalated with sandy turbidites has much higher porosity and correspondingly lower electrical resistivity and P-wave velocity.At drill-site UBGH1-4 the bottom half of one prominent MTD unit shows two bands of parallel fractures on the resistivity log-images indicating a common dip-azimuth direction of about ∼230° (strike of ∼140°). This strike-direction is perpendicular to the seismically defined flow-path of the MTD to the north-east. At Site UBGH1-14, the log-data suggest two zones with preferred fracture orientations (top: ∼250°, bottom: ∼130°), indicating flow-directions to the north-east for the top zone, and north-west for the bottom zone. The fracture patterns may indicate post-depositional sedimentation that gave rise to a preferred fracturing possibly linked to dewatering pathways. Alternatively, fractures may be related to the formation of pressure-ridges common within MTD units.For the interval of observed MTD units, the resistivity and P-wave velocity log-data yield gas hydrate concentrations up to ∼10% at Site UBGH1-4 and ∼25% at Site UBGH1-14 calculated using traditional isotropic theories such as Archie's law or effective medium modeling. However, accounting for anisotropic effects in the calculation to honor observed fracture patterns, the gas hydrate concentration is overall reduced to less than 5%. In contrast, gas hydrate was recovered at Site UBGH1-4 near the base of gas hydrate stability zone (GHSZ). Log-data predict gas hydrate concentrations of 10–15% over an interval of 25 m above the base of GHSZ. The sediments of this interval are comprised of the hemipelagic mud and interbedded thin sandy turbidites, which did contain pore-filling gas hydrate as identified from pore-water freshening and core infra-red imaging. Seismically, this unit reveals a coherent parallel bedding character but has overall faint reflection amplitude. This gas-hydrate-bearing interval can be best mapped using a combination of regular seismic amplitude and seismic attributes such as Shale indicator, Parallel-bedding indicator, and Thin-bed indicator.     

AVO analysis of BSR to assess free gas within fine-grained sediments in the Shenhu area,South China Sea

Gas hydrates have been identified from two-dimensional (2D) seismic data and logging data above bottom simulating reflector (BSR) during China’s first gas hydrate drilling expedition in 2007. The multichannel reflection seismic data were processed to be preserved amplitudes for quantitatively analyzing amplitude variation with offset (AVO) at BSRs. Low P-wave velocity anomaly below BSR, coinciding with high amplitude reflections in 2D seismic data, indicates the presence of free gas. The absolute values of reflection coefficient versus incidence angles for BSR range from 0 to 0.12 at different CMPs near Site SH2. According to logging data and gas hydrate saturations estimated from resistivity of Site SH2, P-wave velocities calculated from effective media theory (EMT) fit the measured sonic velocities well and we choose EMT to calculate elastic velocities for AVO. The rock-physics modeling and AVO analysis were combined to quantitatively assess free gas saturations and distribution by the reflection coefficients variation of the BSRs in Shenhu area, South China Sea. AVO estimation indicates that free gas saturations immediately beneath BSRs may be about 0.2 % (uniform distribution) and up to about 10 % (patchy distribution) at Site SH2.     

冲绳海槽天然气水合物稳定带特征和资源量评价

根据冲绳海槽多道地震资料的处理解释,在16条地震剖面上发现了水合物拟海底反射层BSR,经过AVO、波形反演等特殊的处理技术,首次直接利用BSR圈定了冲绳海槽天然气水合物具体分布范围,直接利用数据得出了天然气水合物稳定带厚度在冲绳海槽的分布趋势,认为海槽南部最厚,中部次之,北部最薄,并通过计算得出了冲绳海槽水合物稳定带的厚度和水合物资源量,对今后海槽水合物勘查和资源量评价具有一定的指导意义。     

Geothermal modeling of the gas hydrate stability zone along the Krishna Godavari Basin

A wide-spread bottom simulating reflector (BSR), interpreted to mark the thermally controlled base of the gas hydrate stability zone, is observed over a close grid of multichannel seismic profiles in the Krishna Godavari Basin of the eastern continental margin of India. The seismic data reveal that gas hydrate occurs in the Krishna Godavari Basin at places where water depths exceed 850 m. The thickness of the gas hydrate stability zone inferred from the BSR ranges up to 250 m. A conductive model was used to determine geothermal gradients and heat flow. Ground truth for the assessment and constraints on the model were provided by downhole measurements obtained during the National Gas Hydrate Program Expedition 01 of India at various sites in the Krishna Godavari Basin. Measured downhole temperature gradients and seafloor-temperatures, sediment thermal conductivities, and seismic velocity are utilized to generate regression functions for these parameters as function of overall water depth. In the first approach the base of gas hydrate stability is predicted from seafloor bathymetry using these regression functions and heat flow and geothermal gradient are calculated. In a second approach the observed BSR depth from the seismic profiles (measured in two-way travel time) is converted into heat flow and geothermal gradient using the same ground-truth data. The geothermal gradient estimated from the BSR varies from 27 to 67°C/km. Corresponding heat flow values range from 24 to 60 mW/m². The geothermal modeling shows a close match of the predicted base of the gas hydrate stability zone with the observed BSR depths.