海洋天然气水合物的地球物理研究(III):似海底反射

对天然气水合物研究中与似海底反射有关的一些观点进行讨论 ,以推动天然气水合物地震研究的认识 .30年的似海底反射研究表明 ,似海底反射仍然是指示天然气水合物沉积存在的最好手段之一 .有似海底反射 ,是可以认为存在天然气水合物的 .虽然存在“游离气带速度模型”与“水合物楔速度模型” ,但似海底反射主要由天然气水合物稳定带底界下方的游离气引起 .BSR上方的天然气水合物、下方的游离气与天然气再循环和含甲烷的流体流有关 .由于天然气水合物稳定带计算控制因素难以准确确定等因素 ,似海底反射与天然气水合物稳定带底界只是近似的对应关系 .需从动态的观点考虑天然气水合物 天然气体系及其与似海底反射的关系 .     

水合物层下伏游离气渗漏过程的数值模拟及实例分析

海洋环境中天然气水合物层是理想的毛细管封闭层,游离气被抑制在水合物层下,游离气层的气体压力随气体聚集和气层厚度的增加而升高,当气压超过封闭层的毛细管力时,游离气会克服毛细管进入压力、刺入上伏封闭层孔隙空间,毛细管封闭作用随之消失,从而形成水合物下伏游离气向海底的渗漏.通过对该过程进行的数值模拟计算表明：渗漏气体是以活塞式驱动上伏沉积层中的孔隙水向海底排出,水合物稳定带内流体渗漏速度随水流柱高度的减小而增加,当水流阻抗大于相应沉积层段的静岩压力时,沉积层将转变为流沙,流沙沉积被海流移除后便在海底留下凹陷麻坑.麻坑形成后流体运移通道演化为气体通道,气体快速排放.麻坑深度主要取决于游离气层的厚度和水合物封闭层（底界）的深度,而与沉积层的渗透率无关.麻坑深度一定程度上指示了渗漏前水合物层下伏游离气层的资源量.对布莱克海台海底麻坑深度的数值模拟计算表明,形成4 m深的海底麻坑需要至少22 m厚的游离气层.     

Quantitative assessment of gas washing of oils in the Tazhong area of the Tarim Basin,Northwest China

Gas washing has been known in the Tazhong area of the Tarim Basin, but its quantitative assessment has not yet been reported. Here the influence of gas washing fractionation in the area was discussed based on the gas chromatogram data of 68 oils and the results of the mixing experiments of a black oil and a condensate. The results show that the intensity of gas washing fractionation decreased generally from northern to southern part and vertically from deep reservoirs to shallow reservoirs. The gas washing fractionation was mainly controlled by fault systems in this area, with the increase of n-alkane mass depletion positively correlated to the number and scale of faults. Gas washing fractionation appears to have affected the hydrocarbon property, and as a result the diversity of the crude oils is markedly controlled by gas washing. In addition, the occurrence of waxy oil in this area may be resulted from multiple factors including gas washing, mixed filling and migration fractionation. Supported by National Basic Research Program of China (Grant No. 2006CB202303) and the National Natural Science Foundation of China (Grant No. 40672091)     

南海北部琼东南海域活动冷泉特征及形成模式

近年来,活动冷泉的研究越来越受到关注.本文利用多波束数据、多道地震数据以及底质取样结果研究琼东南海域活动冷泉系统,分析活动冷泉的羽状流特征、海底地貌与底质特征以及流体活动构造特征.多波束水体数据上,观测到多个延伸高度超过750 m的气泡羽状流,海底流体活动非常强烈;多道地震上识别出麻坑、流体运移通道、气烟囱等流体渗漏相关的构造,与其他海域观测到的反射特征不同,羽状流的下方流体运移通道呈强振幅"串珠"反射;重力活塞取样在两个站位上获得浅表层块状天然气水合物.其中一个站位位于活动冷泉附近,天然气水合物赋存于海底以下8 m左右.基于以上三方面的数据,笔者提出了一个用于描述活动冷泉系统的形成模式,游离气通过气烟囱向上运移到达浅层,一部分在天然气水合物稳定带内形成天然气水合物,另一部分穿透天然气水合物稳定带到达海底,形成活动冷泉的羽状流.     

基于双相介质理论的AVO属性交汇图识别天然气水合物

天然气水合物(gas hydrates)作为一种新的潜在能源 是目前世界范围内研究的热点和难点之一。目前 普遍认为拟海底反射(简称BSR)是天然气水合物存 在的重要标志之一,但有关水合物的识别仍有很 多问题需要解决。例如世界各地已发现水合物的 例证表明在没有BSR的情况下仍然存在水合物,在 没有BSR的情况下如何识别水合物仍是没有解决 的研究课题。本文以天然气水合物的识别为研究 目的,以AVO属性交汇图为研究手段,通过双相介 质正演理论模型研究了当地层含有天然气水合物 或游离气时的AVO属性特征,并与由双相介质的 交错网格有限差分法得到的合成地震记录所反演 得出的AVO交汇图进行了对比,得到高度一致的 对比结果,表明利用AVO属性交汇图是识别天然 气水合物和游离气并能估测其含量的有效技术。     

Simulation experiments on gas production from hydrate-bearing sediments

Experiments were made on 58 sediment samples from four sites(1244,1245,1250 and 1251) of ODP204 at five temperature points(25,35,45,55 and 65℃) to simulate methane production from hydrate-bearing sediments.Simulation results from site 1244 show that the gas components consist mainly of methane and carbon dioxide,and heavy hydrocarbons more than C2+ cannot be detected.This site also gives results,similar to those from the other three,that the methane production is controlled by experimental temperatures,gene...     

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

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

An explicit four-dimensional variational data assimilation method based on the proper orthogonal decomposition: Theoretics and evaluation

The proper orthogonal decomposition (POD) method is used to construct a set of basis functions for spanning the ensemble of data in a certain least squares optimal sense. Compared with the singular value decomposition (SVD), the POD basis functions can capture more energy in the forecast ensemble space and can represent its spatial structure and temporal evolution more effectively. After the analysis variables are expressed by a truncated expansion of the POD basis vectors in the ensemble space, the control variables appear explicitly in the cost function, so that the adjoint model, which is used to derive the gradient of the cost function with respect to the control variables, is no longer needed. The application of this new technique significantly simplifies the data assimilation process. Several assimilation experiments show that this POD-based explicit four-dimensional variational data assimilation method performs much better than the usual ensemble Kalman filter method on both enhancing the assimilation precision and reducing the computation cost. It is also better than the SVD-based explicit four-dimensional assimilation method, especially when the forecast model is not perfect and the forecast error comes from both the noise of the initial filed and the uncertainty of the forecast model. Supported by the National Natural Science Foundation of China (Grant No. 40705035), National High Technology Research and Development Program of China (Grant No. 2007AA12Z144), Knowledge Innovation Project of Chinese Academy of Sciences (Grant Nos. KZCX2-YW-217 and KZCX2-YW-126-2), and National Basic Research Program of China (Grant No. 2005CB321704)     

Check of Earth’s free oscillations excited by Sumatra-Andaman Large Earthquake and discussions on the anisotropy of inner core

Sumatra-Andaman Large Earthquake on Dec. 26, 2004 generated not only the Indian Ocean Tsunami but also the Earth’s free oscillations (EFO). The signals of Earth’s free oscillations were perfectly recorded by the superconducting gravimeter C0-32 at Wuhan station in China. After the pre-treatment and spectral analysis on the observational data from Wuhan station, we obtained more than ninety EFO modes including 42 fundamental modes, 2 radial modes and 49 harmonic modes. On the basis of the discussions on some observed harmonic modes and abnormal splitting phenomena, we considered that the real rigidity might be lower than the theoretical prediction of PREM model in the inner core and however the anisotropy of compressive wave was brightly higher than the present estimations in the inner core. This suggested that the anisotropy of the inner core could be much more complicated than our present understanding, and there might be some new geophysical phenomena in the formation process of the inner core. Supported jointly by the National Natural Science Foundation of China (Grant Nos. 40404005 and 40374029) and the Excellent Prize of President Scholarship and Hundred Talents Program of the Chinese Academy of Sciences     

海洋环境甲烷水合物溶解度及其对水合物发育特征的控制

除合适的温度和压力条件外,甲烷水合物的形成还需要有充足的甲烷供给,沉积物孔隙水中的甲烷浓度必须大于甲烷水合物的溶解度.本文建立了水合物－水－游离气三相体系、水合物－水二相体系、气－水二相体系的甲烷溶解度计算优选方法,计算确定了水合物系统的甲烷溶解度-深度相图,依此划分出游离气、溶解气、水-水合物、水-水合物-游离气四个甲烷不同相态分布区.对水合物脊ODP1249和1250钻位、布莱克海台ODP997钻位稳定带甲烷水合物含量和稳定带之下游离气含量进行了计算.ODP1249浅部13.5～72.4 mbsf(mbsf表示海底以下深度)的甲烷水合物是沉积物孔隙体积的10%～61%,ODP1250钻位35～1065 mbsf的甲烷水合物约为孔隙体积的0.7%～1.9%,水合物层之下游离气层厚约22 m,游离气含量约占孔隙的4%.布莱克海台ODP997钻位的浅部146.9 mbsf处无水合物发育,202.4～433.3 mbsf之间水合物占孔隙体积的约5%～7%,水合物层之下游离气层厚约80 m,游离甲烷含量为孔隙的0.2%～28%.     

Gas hydrate saturation at Site C0002, IODP Expeditions 314 and 315, in the Kumano Basin,Nankai trough

The degree of gas hydrate saturation at Integrated Ocean Drilling Program (IODP) Site C0002 in the Kumano Basin, Nankai Trough, was estimated from logging‐while‐drilling logs and core samples obtained during IODP Expeditions 314 and 315. Sediment porosity data necessary for the calculation of saturation were obtained from both core samples and density logs. Two forms of the Archie equation (‘quick‐look’ and ‘standard’) were used to calculate gas hydrate saturation from two types of electrical resistivity log data (ring resistivity and bit resistivity), and a three‐phase Biot‐type equation was used to calculate gas hydrate saturation from P‐wave velocity log data. The gas hydrate saturation baseline calculated from both resistivity logs ranges from 0% to 35%, and that calculated from the P‐wave velocity log ranges from 0% to 30%. High levels of gas hydrate saturation (>60%) are present as spikes in the ring resistivity log and correspond to the presence of gas hydrate concentrations within sandy layers. At several depths, saturation values obtained from P‐wave velocity data are lower than those obtained from bit resistivity data; this discrepancy is related to the presence of free gas at these depths. Previous research has suggested that gas from deep levels in the Kumano Basin has migrated up‐dip towards the southern and seaward edge of the basin near Site C0002. The high saturation values and presence of free gas at site C0002 suggest that a large gas flux is flowing to the southern and seaward edge of the basin from a deeper and/or more landward part of the Kumano Basin, with the southern edge of the Kumano Basin (the location of site C0002) being the main area of fluid accumulation.     

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

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

尼日尔三角洲天然气水合物体系地震特征研究

天然气水合物体系一般由地震上的似海底反射层BSR和下伏强振幅带(解释为游离气FGZs)所指示,但并非所有BSR及下伏强振幅带都与水合物和游离气有关.本文通过提取多种属性剖面,优选出视极性属性来辅助判断与水合物有关的BSR和下伏FGZs.—般来说,BSR和FGZs的顶部反射层表现为负视极性、高振幅的强反射特征.研究发现,尼日尔三角洲南部的水合物主要分布在与重力作用有关的生长断层及伴生的滚动背斜地区,广泛发育的断层、气烟囱、不整合面,以及砂岩层都可作为流体运移通道.除为水合物体系提供气源,这些流体运移通道还在水合物分解或FGZs超压时发挥作用,使游离气沿通道在FGZs和含水合物地层间循环,部分游离气可运移至海底进入海水甚至大气中·总之,视极性是判断水合物相关地震指示标志的有效属性,流体运移通道对尼日尔三角洲的水合物成藏具有重要作用.     

The Pengguan tectonic dome of Longmen Mountains, Sichuan Province: Mesozoic denudation of a Neoproterozoic magmatic arc-basin system

Neoproterozoic igneous and metamorphic complexes occur as tectonic domes in the Longmen Mountains of the western margin of the Yangtze Block, and are important in reconstructing the Rodinian supercontinent and constraining the timing and mechanism of tectonic denudational processes. The Pengguan dome consists of granitic intrusions and metamorphic rocks of the Huangshuihe Group and is tectonically overlain by ductilly deformed Sinian to Paleozoic strata. The plutonic intrusions consist of granites with abundant amphibolite enclaves. New LA-ICP-MS zircon U-Pb dating yielded an emplacement age of 809±3 Ma and a protolith age of 844±6 Ma for the granite. The granitic rocks have geochemical signatures typical of A-type granites, indicating their formation under an extensional environment, by melting of newly formed tonalite-trondhjemite-granodiorite (TTG) rocks. A detachment fault, characterized by variable ductile shear deformation of S-C fabric and ESE-ward kinematics, separates the Pengguan dome from the Sinian-Paleozoic cover. 40Ar/39Ar dating of muscovite from the mylonite in the detachment fault of the dome demonstrates that ductile deformation occurred at ~160 Ma. This study indicates the existence of a Neoproterozoic magmatic arc-basin system, which was denudated by a Jurassic middle crustal ductile channel flow along the Longmenshan thrust belt.     

Passing gas through the hydrate stability zone at southern Hydrate Ridge,offshore Oregon

We present an equilibrium model of methane venting through the hydrate stability zone at southern Hydrate Ridge, offshore Oregon. Free gas supplied from below forms hydrate, depletes water, and elevates salinity until pore water is too saline for further hydrate formation. This system self-generates local three-phase equilibrium and allows free gas migration to the seafloor. Log and core data from Ocean Drilling Program (ODP) Site 1249 show that from the seafloor to 50 m below seafloor (mbsf), pore water salinity is elevated to the point where liquid water, hydrate and free gas coexist. The elevated pore water salinity provides a mechanism for vertical migration of free gas through the regional hydrate stability zone (RHSZ). This process may drive gas venting through hydrate stability zones around the world. Significant amount of gaseous methane can bypass the RHSZ by shifting local thermodynamic conditions.     

外界条件变化对天然气水合物相平衡曲线及稳定带厚度的影响

目前，有关天然气水合物的相关研究越来越多，而天然气水合物相平衡曲线和稳定带厚度的研究也变得越来越重要.本文利用Sloan的CSMHYD程序研究了外界条件变化对天然气水合物相平衡曲线及稳定带厚度的影响.研究结果表明：当天然气水合物中含有其他气体时，除了氮气会使水合物稳定存在的区域变小外；其他气体都会使稳定区域变大，且甲烷含量越少，水合物越容易形成；对于本文中所提到的几种气体，丙烷和硫化氢对相平衡曲线的影响最大；另外，水合物稳定存在的区域会随着盐度增加而变小.地温梯度、水深、海底温度、气体组成和孔隙水盐度对稳定带厚度的影响不同，其中稳定带厚度与地温梯度呈指数相关关系，与水深呈对数相关关系，与海底温度、水合物中甲烷含量及气体组成呈线性相关关系.水深从1 000 m增加到4 000 m时，稳定带厚度增加了大约400 m；水深2 000 m情况下，地温梯度从0.02℃/m到0.1℃/m变化时，稳定带厚度减薄了大约660 m；底水温度从0~17℃的变化过程中，稳定带厚度减薄了大约1000m；在水合物中气体组成从纯甲烷到含20%乙烷时，稳定带厚度增加了大约170m；盐度在0~4.5 wt％的变化中，稳定带厚度减薄了大约130 m.由此可见，在这几种因素中，地温梯度和底水温度对稳定带厚度的影响较大.     

Carbon isotopic composition and its implications on paleoclimate of the underground ancient forest ecosystem in Sihui,Guangdong

We present the carbon isotopic composition of the total organic carbon (TOC) and fine roots in the sedimentary profile from the underground ancient forest in Sihui to study the climatic and environmental changes from 4.5 ka BP to 0.6 ka BP. Results show that C₃ plant was the main vegetation from 4.5 ka BP to 0.6 ka BP in this region. The ancient forest began to develop in the wetland at around 4 ka BP and disappeared together with the wetland at about 3.0 ka BP, implying that the climate had changed greatly at around 3.0 ka BP. As indicated by the simulation results, the content of atmospheric CO₂ increased slightly during 3.5 ka BP to 3.0 ka BP, implying climate warming during that period. The interval of radiocarbon age between 3.0 ka BP to 1.2 ka BP was possibly caused by the strong erosion when the block was lifted in the neotectonic movement. From 1.2 ka BP to 0.6 ka BP, the region remained in terrestrial sedimentary environment, and the surface plant biomass declined gradually. Drought caused by the climate change was the likely cause for the disappearance of the ancient forest. South transition of Intertropical Convergence Zone (ITCZ) was probably the main mechanism for the climate change. Supported by National Natural Science Foundation of China (Grant Nos. 40231015 and 40473002), National Basic Research Program of China (Grant No. 2005CB422004), the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KSCX2-SW-133) and Open Funds of State Key Laboratory of Organic Geochemistry (Grant No. OGL-200607)     

Full waveform inversion of gas hydrate reflectors in Northern South China Sea

Bottom Simulating Reflectors (BSRs) are considered to be the bottom of gas hydrate bearing sediments; hence, BSRs are used to identify gas hydrate and free gas. In order to obtain accurate velocity structure of BSRs, this paper presents a full waveform inversion strategy based on Genetic Algorithm. Synthetic seismograms are calculated using the slowness technique. Through numerical experiments made with noisy synthetic data, the inversion algorithm shows stable performance, and genetic operators are defined. This method was applied to field data from the northern South China Sea. Inversion results show that obvious velocity anomaly of BSRs can be detected, which indicates the existence of gas hydrate and free gas.     

Variability of the currents in the Luzon Strait during spring of 2002 obtained from observations and satellite geostrophic currents and spectral analyses

The structure and variability of the currents in the Luzon Strait during spring of 2002 are studied, based on the current measurements at the average position of the mooring station (20°49′57"N, 120°48′12"E) from March 17 to April 15, 2002, satellite geostrophic currents in the Luzon Strait, and the spectral analyses, using the maximum entropy method. The subtidal currents at the mooring station show de-creased amplitudes downward with an anti-cyclonic rotation, suggesting that the currents enter and exit t...     

天然气水合物形成条件与含量影响因素的半定量分析

天然气水合物形成条件及其含量的研究是天然气水合物调查研究中最为关心的两个问题.从天然气水合物形成过程的热动力学理论模型出发,半定量地探讨了不同因素（温度、压力、气体组成、孔隙水盐度、沉积物孔隙大小等）对天然气水合物形成作用及含量的影响程度.结果显示,气体组成特别是丙烷的加入对天然气水合物形成的温度和压力条件影响最大,孔隙水盐度也会对天然气水合物形成的温度和压力条件产生重要影响,沉积物孔隙在一定范围内（1×10^-6m～4×10^-8m）对天然气水合物形成的温度和压力条件影响有限.天然气水合物含量受孔隙大小和盐度影响较小,主要只与气体的供应大小有关.