南海神狐海域水合物一维速度结构研究： 基于中—低饱和度水合物区的约束

高精度OBS探测作为目前研究水合物的常用技术手段,可以获取水合物矿体精细速度结构,在研究天然气水合物饱和度、水合物资源的预测与开发等方面具有重要指导意义。然而正演模拟OBS速度结构是繁琐漫长的过程,构建合理的初始模型是后续精细结构快速成像的重要前提。本研究基于全球18个已探明地震波速度的水合物区,分析了水合物矿体内纵波速度的共性特征和影响因素,拟合了水合物矿体带纵波速度与海水深度、沉积物厚度的经验公式。综合经验公式、OBS数据与多道地震剖面,建立了神狐水合物区横向均匀初始模型,并通过射线追踪与走时拟合模拟了神狐海域的一维纵波速度结构。结果表明,神狐水合物矿体带具有高纵波速度(1.83~1.92 km/s）,游离气层具有低P波速度（1.60~1.70 km/s）,此外,基于全球的水合物速度经验公式对神狐海域速度结构模拟具有重要参考意义,有望为获取神狐海域二维/三维精确速度结构提供可靠的初始模型,进而指导水合物精细勘察与资源评估。     

Assessment of natural gas hydrate reservoirs at Site GMGS3-W19 in the Shenhu area,South China Sea based on various well logs

To obtain the characteristics of the gas hydrate reservoirs at GMGS3-W19, extensive geophysical logging data and cores were analyzed to assess the reservoir properties. Sediment porosities were estimated from density, neutron, and nuclear magnetic resonance (NMR) logs. Both the resistivity and NMR logs were used to calculate gas hydrate saturations, the Simandoux model was employed to eliminate the effects of high clay content determined based on the ECS and core data. The density porosity was closely in agreement with the core-derived porosity, and the neutron porosity was higher while the NMR porosity was lower than the density porosity of sediments without hydrates. The resistivity log has higher vertical resolution than the NMR log and thus is more favorable for assessing gas hydrate saturation with strong heterogeneity. For the gas hydrate reservoirs at GMGS3-W19, the porosity, gas hydrate saturation and free gas saturation was 52.7%, 42.7% and 10%, on average, respectively. The various logs provide different methods for the comprehensive evaluation of hydrate reservoir, which supports the selection of candidate site for gas hydrate production testing.©2022 China Geology Editorial Office.     

Petrophysical parameters and modelling of the Eocene reservoirs in the Qadirpur area,Central Indus Basin,Pakistan: implications from well log analysis

The Eocene rock units of the Qadirpur field, Central Indus Basin (Pakistan), are investigated petrophysically for their detailed reservoir characterization. The different petrophysical parameters determined include the following: true resistivity, shale volume, total porosity, effective porosity, density and neutron porosity, water and hydrocarbon saturation, bulk volume of water, lithology, gas effect, P-wave velocity, movable hydrocarbon index and irreducible water saturation and integrated with different cross-plots. The Eocene reservoirs are excellent with high effective porosity (2–32 %) and hydrocarbon saturation (10–93 %). Among these, the Sui Upper Limestone is an overall a poor reservoir; however, it has some hydrocarbon-rich intervals with high effective porosity and better net pay. All the net pay zones identified show low and variable shale volume (5–30 %). The secondary porosity has added to the total and effective porosities in these reservoirs. The main contributors to the porosity are the chalky, intercrystalline and vuggy/fracture types. The thickness of the reservoirs zones ranges from 4.5 to 62 m. These reservoirs are gas-producing carbonates with almost irreducible water saturation (0.002–0.01) and are likely to produce water-free hydrocarbons. The lower values of moveable hydrocarbon index (0.07–0.9) show that the hydrocarbons are moveable spontaneously to the well bore. The proposed correlation model shows that the reservoirs have an inclined geometry and are a part of an anticlinal trap.     

天然气水合物地震响应研究——中国南海HD152测线应用实例

目前识别水合物的主要依据是BSR和其上具有高速特征的空白带。当BSR不存在, 或由于水合物饱和度较高导致空白带特征不明显时, 水合物的识别是未解决的问题。针对这种情况, 提出了一种更实用的水合物识别模式: 低速背景中近似平行于海底的高速地质体是水合物带的特征, 水合物带泊松比降低, 利用纵横波速度信息易于识别水合物带下是否存在游离气。基于上述水合物的识别模式, 以中国南海HD152测线为实例进行了天然气水合物的识别研究。     

The Sensitive Properties of Hydrate Reservoirs Based on Seismic Stereoscopic Detection Technology

Higher-precision determinations of hydrate reservoirs, hydrate saturation levels and storage estimations are important for guaranteeing the ability to continuously research, develop and utilize natural gas hydrate resources in China. With seismic stereoscopic detection technology, which fully combines the advantages of different seismic detection models, hydrate formation layers can be observed with multiangle, wide-azimuth, wide-band data with a high precision. This technique provides more reliable data for analyzing the distribution characteristics of gas hydrate reservoirs, establishing velocity models, and studying the hydrate-sensitive properties of petrophysical parameters;these data are of great significance for the exploration and development of natural gas hydrate resources. Based on a velocity model obtained from the analysis of horizontal streamer velocity data in the hydrate-bearing area of the Shenhu Sea, this paper uses three VCs(longitudinal spacing of 25 m) and four OBSs(transverse spacing of 200 m) to jointly detect seismic datasets consisting of wave points based on an inversion of traveltime imaging sections. Accordingly, by comparing the differences between the seismic phases in the original data and the forward-modeled seismic phases, multiangle coverage constraint corrections are applied to the initial velocity model, and the initial model is further optimized, thereby improving the imaging quality of the streamer data. Petrophysical elastic parameters are the physical parameters that are most directly and closely related to rock formations and reservoir physical properties. Based on the optimized velocity model, the rock elastic hydrate-sensitive parameters of the hydrate reservoirs in the study area are inverted, and the sensitivities of the petrophysical parameters to natural gas hydrates are investigated. According to an analysis of the inversion results obtained from these sensitive parameters, λρ, Vp and λμ are simultaneously controlled by the bulk modulus and shear modulus, while Vs and μρ are controlled only by the shear modulus, and the latter two parameters are less sensitive to hydrate-bearing layers. The bulk modulus is speculated to be more sensitive than the shear modulus to hydrates. In other words, estimating the specific gravity of the shear modulus among the combined parameters can affect the results from the combined elastic parameters regarding hydrate reservoirs.     

基于电阻率测井的天然气水合物饱和度估算及估算精度分析

利用电阻率测井基于阿尔奇公式能够估算孔隙空间中天然气水合物饱和度。基于神狐海域的水合物钻井(SH2站位)资料,利用密度反演的地层孔隙度与地层因子交会图,计算出估算水合物饱和度的阿尔奇常数。假定地层中水合物呈均匀分布,利用电阻率资料计算的该站位水合物饱和度平均为24%,最高达44%,水合物饱和度在垂向上具有明显的不均匀性;利用孔隙水氯离子异常估算的SH2站位水合物饱和度平均值为25%,最高值达48%,厚度为25 m。这两种方法估算的水合物饱和度基本吻合,利用电阻率资料计算的饱和度与阿尔奇常数和饱和度指数及孔隙度有关。     

Research Progress on the Electrical Properties of Gas Hydrate-bearing Sediments

Electrical properties are important physical parameters of natural gas hydrate, and, specifically, resistivity has been widely used in the quantitative estimation of hydrate saturation. There are three main methods to study the electrical properties of gas hydrate-bearing sediments: experimental laboratory measurements, numerical simulation, and resistivity logging. Experimental measurements can be divided into three categories: normal electrical measurement, complex resistivity measurement, and electrical resistivity tomography. Experimental measurements show that the resistivity of hydrate-bearing sediment is affected by many factors, and its distribution as well as the hydrate saturation is not uniform; there is a distinct non-Archie phenomenon. The numerical method can simulate the resistivity of sediments by changing the hydrate occurrence state, saturation, distribution, etc. However, it needs to be combined with X-ray CT, nuclear magnetic resonance, and other imaging techniques to characterize the porous characteristics of the hydrate-bearing sediments. Resistivity well logging can easily identify hydrate layers based on their significantly higher resistivity than the background, but the field data of the hydrate layer also has a serious non-Archie phenomenon. Therefore, more experimental measurements and numerical simulation studies are needed to correct the parameters of Archie''s formula.     

南海北部天然气水合物富集特征及定量评价

为寻找有资源前景的高富集天然气水合物及水合物储层的精细刻画方法,利用南海6次钻探发现的高饱和度水合物层的测井、岩心和三维地震数据,分析水合物富集层测井与地震异常特征.发现:（1）不同饱和度的孔隙与裂隙充填型水合物层的测井和地震异常不同,裂隙充填型水合物层具有各向异性;（2）受高通量流体运移的影响,在粉砂沉积物的水合物稳定带底界附近能形成中等饱和度的水合物层,识别标志为稳定带内极性与海底一致的强振幅反射,而非BSR和振幅空白;（3）裂隙充填型中等饱和度水合物层在地震剖面上表现为地层上拱和弱-中等强度振幅反射.储层-疏导-气源的耦合控制着水合物的富集特征和分布,断层与流体运移控制着细粒粉砂质沉积物中水合物的富集与厚度.基于饱和度岩相的统计学反演,能识别3 m非水合物和低饱和度水合物层及空间分布.      

基于声波测井的冻土区孔隙型水合物储层水合物饱和度估算方法

水合物饱和度参数的准确计算对于水合物资源量的评价至关重要。本文提出利用超声波测井资料与等效介质模型相结合的方法,可有效评价祁连山冻土区孔隙型水合物储层水合物饱和度变化特征,并在典型孔隙型水合物钻孔DKXX 13进行了应用。基于等效介质理论的弹性波速度模型正演模拟的纵波速度相比基于双相介质理论的弹性波速度模型更加吻合实际测井纵波速度,可用于分析孔隙型水合物储层的纵波速度特征;通过正演模拟的纵波速度与实际测井纵波速度的对比,识别出X30. 0~X30. 2m、X30. 3~X30. 4m、X31. 1~X31. 6m、X31. 7~X31. 9m、X32. 0~X32. 2m井段存在水合物,水合物赋存井段地层的水合物饱和度变化范围为13. 0%~85. 0%,平均值为61. 9%,与标准阿尔奇公式估算结果和现场岩芯测试结果基本一致。研究结果可为祁连山冻土区水合物地层测井评价与地震勘探提供理论依据和技术支撑。     

南海神狐海域天然气水合物声波测井速度与饱和度关系分析

利用南海北部神狐海域A站位的地震和测井资料综合分析神狐海域含天然气水合物沉积层的声波测井速度及水合物饱和度的分布特征和变化规律,并对水合物饱和度的理论计算值和实测值进行对比分析,同时对水合物稳定带的纵波速度特征与饱和度的关系进行了综合研究。结果表明：神狐海域A站位的水合物层厚度约20 m,纵波速度在1 873～2 226m/s之间,水合物饱和度在15.0%～47.3%之间变化,水合物饱和度值相对较高;受海底复杂地质因素的影响,根据岩心孔隙水的氯离子淡化程度实测的水合物饱和度随声波速度的变化并不是单一的正比例关系,而是随声波速度的升高而上下波动,波动幅度在10%～20%之间,总体趋势上随声波速度的升高而升高,并集中分布在理论曲线附近;利用热弹性理论速度模型计算并校正后的水合物饱和度随声波速度的增加而有规律地增加,水合物饱和度的理论计算值与实测数据比较吻合,说明所建立的岩石物理模型正确,模型参数选取合理。根据声波速度计算水合物饱和度这一方法可扩展到整个研究区域,并为研究区的水合物资源量评价提供基础数据。     

南海北部陆坡天然气水合物的地震速度研究

采用D ix公式法和速度反演法分别对南海北部陆坡测线A和B的层速度及纵波速度进行计算,结合BSR、振幅空白带以及波形极性反转等多种水合物赋存信息的分析,对水合物成矿带的速度特征进行了综合研究。结果表明:低速背景中的高速异常,是天然气水合物赋存的重要特征;高速异常体一般呈平行于海底的带状分布;在高速异常体的内部,速度也是不断变化的,一般在异常体的中心速度最高,由中心到边缘速度逐渐降低,反映在水合物矿带内部,水合物饱和度由矿体中心向边缘逐渐降低的特征。研究结果表明高精度速度分析不仅可以帮助寻找水合物矿点,还可以进一步判定水合物的富集层位。     

超声探测技术在天然气水合物模拟实验中的应用

为了解不同介质中天然气水合物的声学特性,在特制的高压反应釜中分别进行了纯水、松散沉积物和岩心中甲烷水合物的生成和分解的模拟实验,同时应用超声技术进行了探测。在纯水-甲烷体系中,声波速度的变化主要受温度的制约,水中生成的絮状水合物并没有使声波速度发生明显变化;在纯水-松散沉积物-甲烷体系中,声波速度和系统主频的变化灵敏地反映出体系内水合物的生成和分解;在纯水-岩心-甲烷体系中,随着水合物的生成,纵波速度、横波速度以及纵波幅度均增大,这说明纵波和横波的速度随着孔隙度的减小而增大,而纵波幅度的衰减则随着孔隙度的减小而减小。实验结果显示,超声探测是天然气水合物模拟实验中的一项有效的探测技术。     

南海神狐海域含水合物地层测井响应特征

分析了南海北部神狐海域含天然气水合物沉积层声波速度及密度的分布特征和变化规律,并通过对比DSDP 84航次570号钻孔含天然气水合物层段测井资料,总结出神狐海域含水合物地层的测井响应规律特征:神狐海域含水合物地层存在着明显的高声波速度、低密度特征,地层密度随声波速度的变化并不是单一的反比例关系,总体趋势上随声波速度的升高而降低;含水合物地层高声波速度值主要集中在197~220 m段,饱和度值在15%~47%之间,低密度值集中在200~212 m段,分布在水合物饱和度大于20%的地层内;含水合物地层声波速度平均值为2 076 m/s,其上覆和下伏地层的声波速度平均值为1 903 m/s和1 892 m/s,所对应的地层密度值分别为1.89 g/cm3、1.98 g/cm3和2.03 g/cm3,声波速度受孔隙度和饱和度的共同影响,地层密度受水合物饱和度影响较大;从水合物上覆地层到声波速度最高值段,声波速度值增加了9.1%,相对应的地层密度值减少了4.55%,从水合物声波速度最高值段到下伏地层,声波速度值减少了8.86%,相对应的地层密度值增加了7.41%。这些测井响应特征,可用来识别地层中天然气水合物,并可以用来计算水合物的饱和度,同时结合其他地质和地球物理资料,确定水合物层的厚度、分布范围,计算天然气水合物的资源量。     

粘弹介质叠前四参数同步反演及应用

叠前地震资料中包含纵波速度、横波速度、密度和吸收衰减系数等弹性及粘弹性信息,其中吸收衰减属性对储层物性及流体性质尤为敏感.本文从粘弹介质精确Zoeppritz方程出发,基于介质分解理论通过公式推导及近似,建立了包含纵波速度、横波速度、密度和吸收衰减系数的纵波反射系数特征方程,利用贝叶斯反演框架,实现了叠前纵波速度、横波...     

射线追踪法在南海天然气水合物速度分析中的应用

采用射线追踪法对南海北部陆坡A测线层速度进行计算,结合BSR(Bottom Simulating Reflector)、振幅空白带,以及波形极性反转等多种水合物赋存信息的分析,对水合物成矿带的速度特征进行了综合研究。结果表明:低速背景中的高速异常,是天然气水合物赋存的重要特征;高速异常体一般呈平行于海底的带状分布;在高速异常体的内部,速度也是不断变化的,一般在异常体的中心速度最高,由中心到边缘速度逐渐降低,反映在水合物矿带内部,水合物饱和度由矿体中心向边缘逐渐降低的特征。研究结果还表明,高精度速度分析不仅可以帮助寻找水合物矿点,还可以进一步判定水合物的富集层位。     

海底天然气水合物储集层饱和度的估算方法

据当前国内外海底天然气水合物饱和度预测技术的文献分析,讨论了岩芯孔隙水的氯离子浓度、地球物理测井和地震波速度等方面在估算海底天然气水合物饱和度方面的应用和计算方法,并对它们的优劣进行简单的评述。     

Coexistence of natural gas hydrate,free gas and water in the gas hydrate system in the Shenhu Area,South China Sea

Shenhu Area is located in the Baiyun Sag of Pearl River Mouth Basin, which is on the northern continental slope of the South China Sea. Gas hydrates in this area have been intensively investigated, achieving a wide coverage of the three-dimensional seismic survey, a large number of boreholes, and detailed data of the seismic survey, logging, and core analysis. In the beginning of 2020, China has successfully conducted the second offshore production test of gas hydrates in this area. In this paper, studies were made on the structure of the hydrate system for the production test, based on detailed logging data and core analysis of this area. As to the results of nuclear magnetic resonance (NMR) logging and sonic logging of Well GMGS6-SH02 drilled during the GMGS6 Expedition, the hydrate system on which the production well located can be divided into three layers: (1) 207.8–253.4 mbsf, 45.6 m thick, gas hydrate layer, with gas hydrate saturation of 0–54.5% (31% av.); (2) 253.4–278 mbsf, 24.6 m thick, mixing layer consisting of gas hydrates, free gas, and water, with gas hydrate saturation of 0–22% (10% av.) and free gas saturation of 0–32% (13% av.); (3) 278–297 mbsf, 19 m thick, with free gas saturation of less than 7%. Moreover, the pore water freshening identified in the sediment cores, taken from the depth below the theoretically calculated base of methane hydrate stability zone, indicates the occurrence of gas hydrate. All these data reveal that gas hydrates, free gas, and water coexist in the mixing layer from different aspects.     

Estimation of the 3-D velocity structure of the Sorachi–Yezo and Oshima Belts in the western part of Hokkaido, Japan

We adopted the seismic tomography technique to refine the three-dimensional velocity structure model of the western part of Hokkaido, Japan. Using the P-wave first arrival data listed by Japan Meteorological Agency from 2002 to 2005, we could estimate a 3-D inhomogeneous velocity structure model with a low velocity at a depth of 14 km beneath Asahikawa. The crustal structure near Sapporo was characterized by lateral velocity change toward the southern seaside. The low-velocity zone near Urakawa, proposed by previous research, was also clarified. In general, the present model showed lower-velocity values for most of the crustal layers in the area concerned. The results of this study were affected by less number of higher magnitude events (M?≥?0.5) in the central part of the area of interest. However, the perturbation results for comparatively shallow layers (6–50 km) were good in resolution. It was found that the source region of the Rumoi–Nanbu earthquake of December 14, 2004 was characterized by a low-velocity zone, located between high velocity zones. Such an inhomogeneous crustal structure might play an important role in the relatively high seismic activity in the Rumoi–Nanbu earthquake source region.     

Modeling of the Shale Volume in the Hendijan Oil Field Using Seismic Attributes and Artificial Neural Networks

Petrophysical properties have played an important and definitive role in the study of oil and gas reservoirs, necessitating that diverse kinds of information are used to infer these properties. In this study, the seismic data related to the Hendijan oil field were utilised, along with the available logs of 7 wells of this field, in order to use the extracted relationships between seismic attributes and the values of the shale volume in the wells to estimate the shale volume in wells intervals. After the overall survey of data, a seismic line was selected and seismic inversion methods (model-based, band limited and sparse spike inversion) were applied to it. Amongst all of these techniques, the model-based method presented the better results. By using seismic attributes and artificial neural networks, the shale volume was then estimated using three types of neural networks, namely the probabilistic neural network (PNN), multi-layer feed-forward network (MLFN) and radial basic function network (RBFN).     

Crustal structure in Tengchong Volcano-Geothermal Area, western Yunnan, China

Based upon the deep seismic sounding profiles carried out in the Tengchong Volcano-Geothermal Area (TVGA), western Yunnan Province of China, a 2-D crustal P velocity structure is obtained by use of finite-difference inversion and forward travel-time fitting method. The crustal model shows that a low-velocity anomaly zone exists in the upper crust, which is related to geothermal activity. Two faults, the Longling–Ruili Fault and Tengchong Fault, on the profile extend from surface to the lower crust and the Tengchong Fault likely penetrates the Moho. Moreover, based on teleseismic receiver functions on a temporary seismic network, S-wave velocity structures beneath the geothermal field show low S-wave velocity in the upper crust. From results of geophysical survey, the crust of TVGA is characterized by low P-wave and S-wave velocities, low resistivity, high heat-flow value and low Q. The upper mantle P-wave velocity is also low. This suggests presence of magma in the crust derived from the upper mantle. The low-velocity anomaly in upper crust may be related to the magma differentiation. The Tengchong volcanic area is located on the northeast edge of the Indian–Eurasian plate collision zone, away from the eastern boundary of the Indian plate by about 450 km. Based on the results of this paper and related studies, the Tengchong volcanoes can be classified as plate boundary volcanoes.