Natural gas leakage of Mizhi gas reservoir in Ordos Basin,recorded by natural gas fluid inclusion

Abundant natural gas inclusions were found in calcite veins filled in fractures of Central Fault Belt across the centre of Ordos Basin. Time of the calcite veins and characteristics of natural gas fluid inclusion were investigated by means of dating of thermolum luminescence (TL) and analyzing stable isotope of fluid inclusion. Results show that natural gas inclusion formed at 130–140°C with salinity of 5.5 wt%–6.0 wt% NaCl. It indicates that natural gas inclusion is a kind of thermal hydrocarbon fluid formed within the basin. Method of opening inclusion by heating was used to analyze composition of fluid inclusion online, of which the maximal hydrocarbon gas content of fluid inclusion contained in veins is 2.4219 m³/t rock and the maximal C₁/ΣC _i ratio is 91%. Laser Raman spectroscopy (LRS) was used to analyze chemistry of individual fluid inclusion in which the maximal hydrocarbon gas content is 91.6% compared with little inorganic composition. Isotope analysis results of calcite veins show that they were deposited in fresh water, in which the δ ¹³C_PDB of calcite veins is from ?5.75‰ to 15.23‰ and δ ¹⁸O_SMOW of calcite veins is from 21.33‰ to 21.67‰. Isotope results show that δ ¹³C_{1 PDB} of natural gas fluid inclusion is from ?21.36‰ to ?29.06‰ and δD_SMOW of that is from ?70.89‰ to ?111.03‰. It indicates that the gas of fluid inclusion formed from coal source rocks and it is the same as that of natural gas of Mizhi gas reservoir. Results of TL dating show that time of calcite vein is (32.4±3.42)×10⁴ a, which is thought to be formation time of gas inclusion. It indicated that natural gas inclusion contained in calcite veins recorded natural gas leakage from Mizhi gas reservoir through the Central Fault Belt due to Himalayan tectonic movement.     

Natural gas leakage of Mizhi gas reservoir in Ordos Basin,recorded by natural gas fluid inclusion

Abundant natural gas inclusions were found in calcite veins filled in fractures of Central Fault Belt across the centre of Ordos Basin. Time of the calcite veins and characteristics of natural gas fluid inclusion were investigated by means of dating of thermolum luminescence (TL) and analyzing stable isotope of fluid inclusion. Results show that natural gas inclusion formed at 130–140°C with salinity of 5.5 wt%–6.0 wt% NaCl. It indicates that natural gas inclusion is a kind of thermal hydrocarbon fluid formed within the basin. Method of opening inclusion by heating was used to analyze composition of fluid inclusion online, of which the maximal hydrocarbon gas content of fluid inclusion contained in veins is 2.4219 m³/t rock and the maximal C₁/ΣC _i ratio is 91%. Laser Raman spectroscopy (LRS) was used to analyze chemistry of individual fluid inclusion in which the maximal hydrocarbon gas content is 91.6% compared with little inorganic composition. Isotope analysis results of calcite veins show that they were deposited in fresh water, in which the δ ¹³C_PDB of calcite veins is from −5.75‰ to 15.23‰ and δ ¹⁸O_SMOW of calcite veins is from 21.33‰ to 21.67‰. Isotope results show that δ ¹³C_{1 PDB} of natural gas fluid inclusion is from −21.36‰ to −29.06‰ and δD_SMOW of that is from −70.89‰ to −111.03‰. It indicates that the gas of fluid inclusion formed from coal source rocks and it is the same as that of natural gas of Mizhi gas reservoir. Results of TL dating show that time of calcite vein is (32.4±3.42)×10⁴ a, which is thought to be formation time of gas inclusion. It indicated that natural gas inclusion contained in calcite veins recorded natural gas leakage from Mizhi gas reservoir through the Central Fault Belt due to Himalayan tectonic movement.

     

New indicators from bedding-parallel beef veins for the fault valve mechanism

Beef structures(bedding-parallel veins of fibrous calcite)are widespread within the Lower Triassic carbonate rocks in the Sichuan Basin of China,especially within clay-rich strata of low permeability.In the veins,fibrous calcite occurs in the outer zones,and coarse equant calcite in the inner zones.At least two generations of calcite crystallization took place during aqueous alteration,at the same time as deformation recorded by the calcite.The first-generation calcite fibers are at steep angles to the hydrocarbon-bearing host beds,and they grew vertically against the force of gravity at a time when the source rocks were maturing.Second-generation calcite occurs as coarse equant grains that sealed pores via localized fluid flow during horizontal tectonic compression,so that shear stresses acted at the fracture margins.Shearing of the host rock was accommodated in part by dissolution-precipitation creep(DPC),grain rotation,and grain slippage,recorded in crystallographic preferred orientations(CPOs)of the host calcite grains beside the crack walls.The beef veins formed during high pore-fluid overpressures along hydrofractures.We propose that the bedding-parallel veins with beef structures are evidence of a"crack-seal slip"fault valve process during hydrocarbon generation.The hydrocarbon-bearing calcite beef structure may be a good indicator of oil or gas migration,and of the flow direction of aqueous solutions.     

Ore-forming fluid and mineralization of Caijiaying and Dajing polymetallic ore deposits

Fluid inclusions from the Dajing and Caijiaying deposits have nearly the same ho-mogenization temperature. Correlation between temperature and salinity shows that both Sn-and Cu-bearing fluids Sn and Cu were present in the Dajing deposit but only one kind of fluids continuously evolved in the Caijiaying deposit. Study on rare earth elements (REE) in ancient fluid from the inclusions indicates that the fluid of Sn mineralizing stage in Dajing was derived from remelting magma of the continental crust, and the fluid of Cu-Pb-Zn mineralizing stage in the Dajing deposit and the fluid of Cu-Pb-Zn mineralization in the Caijiaying deposit were derived from the mantle. It is concluded that the Cu-Pb-Zn mineralizations in the Dajing and Caijiaying deposits resulted from the identical tectono-thermal event of magma-fluids induced by Mesozoic tectonic transition and extension in the eastern part of North China Craton.     

Effect of oil composition on fluid substitution in heavy oil reservoirs

Unlike light oils, heavy oils do not have a well‐established scheme for modelling elastic moduli from dynamic reservoir properties. One of the main challenges in the fluid substitution of heavy oils is their viscoelastic nature, which is controlled by temperature, pressure, and fluid composition. Here, we develop a framework for fluid substitution modelling that is reliable yet practical for a wide range of cold and thermal recovery scenarios in producing heavy oils and that takes into account the reservoir fluid composition, grounded on the effective‐medium theories for estimating elastic moduli of an oil–rock system. We investigate the effect of fluid composition variations on oil–rock elastic moduli with temperature changes. The fluid compositional behaviour is determined by flash calculations. Elastic moduli are then determined using the double‐porosity coherent potential approximation method and the calculated viscosity based on the fluid composition. An increase in temperature imposes two opposing mechanisms on the viscosity behaviour of a heavy‐oil sample: gas liberation, which tends to increase the viscosity, and melting, which decreases the viscosity. We demonstrate that melting dominates gas liberation, and as a result, the viscosity and, consequently, the shear modulus of the heavy oils always decrease with increasing temperature. Furthermore, it turns out that one can disregard the effects of gas in the solution when modelling the elastic moduli of heavy oils. Here, we compare oil–rock elastic moduli when the rock is saturated with fluids that have different viscosity levels. The objective is to characterize a unique relation between the temperature, the frequency, and the elastic moduli of an oil–rock system. We have proposed an approach that takes advantage of this relation to find the temperature and, consequently, the viscosity in different regions of the reservoir.     

渤海新构造运动及其对晚期油气成藏的影响

自中新世晚期（大约5.2 Ma B.P.）以来,渤海海域盆地进入裂后构造活动最活跃的新构造运动时期,致使渤海海域沉降沉积中心迁移、玄武岩喷发、地震频繁活动、深大断裂继承性活动和大量晚期断裂的生成．该期构造运动海域强度明显高于相邻陆域,呈现出幕式活动的特征．在新构造运动的影响和控制作用下,渤海海域含油气盆地形成了一批大型浅层背斜圈闭,油气输导体系得到优化,烃源岩晚期快速生烃,油气幕式充注成藏．本文以PL19-3油气藏为例,初步总结了渤海新构造运动控制油气晚期快速动态成藏的模式．     

Threshold pressure gradient of fluid flow through multi-porous media in low and extra-low permeability reservoirs

After analyzing many studies of fluid flow theory of multi-porous media in low and extra-low permeability reservoirs and the numerical simulation of non-Darcy flow, we found that a negative flow rate occurs in the existing non-Darcy flow equation, which is unreasonable. We believe that the existing equation can only be considered as a discriminant to judging Darcy flow or non-Darcy flow, and cannot be taken as a fluid flow governing equation of multi-porous media. Our analysis of the experimental results shows that the threshold pressure gradient (TPG) of low and extra-low permeability reservoirs is excessively high, and does not conform to fluid flow through multi-porous media in the actual reservoir situation. Therefore, we present a reasonable TPG ranging from 0.006 to 0.04 MPa/m at the well depth of 1500 m and oil drainage distance of 500 m. The results of our study also indicate that the non-Darcy flow phenomenon will disappear when the TPG reaches a certain value. In addition, the TPG or non-Darcy flow in low and extra-low permeability reservoirs does not need to be considered in the productivity prediction and reservoir numerical simulation. At present, the black oil model or dual-porous media is suitable for simulating low and extra-low permeability reservoirs.     

Extreme flood events and climate change around 3500 aBP in the Central Plains of China

The occurrence of prehistoric extreme flood events has been the common interest of geologists, geographers and archaeologists. It is recognized that from 5000 aBP to 3000 aBP was a period of extreme floods frequently occurring around the world. For in-stance, flood events have been found in North Amer-ica, rapid sea level rise has been found in both the Black Sea and the Mediterranean Sea, and flooded ancient cities have been found in a number of coast areas[1—6]. In China, records of pre…     

热释光年代测定方法在鱼洞水库地区断层活动性研究中的应用

用作者研制的RS 81 1型热释光仪对云南昭通鱼洞水库地区两条断层的破碎带内发育的方解石脉进行了热释光年代测定和显微构造分析 .结果显示 :2条断层破碎带内发育 4期方解石脉 ,其中有 3期发生了变形 .F1 断层距今 5× 1 0 4 a以来 ,F2 断层距今 2 5× 1 0 4 a以来没有强烈活动的迹象 .断层泥中石英颗粒表面形态分析结果与热释光法相同 .     

泥质砂岩液相渗透率计算新方法

液相渗透率描述了岩石的渗流特性,是评价储层与预测油气产能的重要参数.液相渗透率是指盐水溶液在岩石孔隙中流动且与岩石孔隙表面黏土矿物发生物理化学作用时所测得的渗透率;液相渗透率的实验测量条件更加接近实际地层泥质砂岩的条件,使得液相渗透率更能反映地层条件下泥质砂岩的渗流特性;然而,现有的液相渗透率评价模型较少,且模型未能揭示液相渗透率与溶液矿化度之间的关系.基于此,开展了液相渗透模型推导与计算方法研究;文中首先将岩石等效为毛管束模型,推导建立了液相渗透率与比表面、喉道曲折度、总孔隙度、黏土束缚水孔隙度等参数之间的关系;其次,根据岩石物理体积模型,推导建立了黏土束缚水孔隙度与阳离子交换容量、溶液矿化度等参数的关系;最终,将黏土束缚水孔隙度引入液相渗透率计算公式,建立了基于总孔隙度、阳离子交换容量、溶液矿化度、比表面、喉道曲折度等参数的液相渗透率理论计算模型.液相渗透率计算模型与两组实验数据均表明,液相渗透率随阳离子交换容量的增大而降低,随溶液矿化度的增大而增大.然而,液相渗透率理论计算模型的实际应用中喉道曲折度、比表面等参数求取困难,直接利用理论模型计算液相渗透率受到限制.在分析液相渗透率与孔隙渗透率模型的基础上,建立了液相渗透率与空气渗透率之间的转换模型,形成了利用转化模型计算液相渗透率的新方法.为进一步验证液相渗透率与空气渗透率转化模型的准确性,基于两组实验数据,利用转换模型计算了液相渗透率;液相渗透率计算结果与岩心测量液相渗透率实验结果对比显示,液相渗透率计算结果与实际岩心测量结果吻合较好,文中建立的液相渗透率与空气渗透率转化模型合理可靠.     

Fluxes of fluid and heat from the oceanic crustal reservoir

Recent discoveries define a global scale fluid reservoir residing within the uppermost igneous oceanic crust, a region of seafloor that is both warm and may harbor a substantial biosphere. This hydrothermal fluid reservoir formed initially within volcanic rocks newly erupted at mid-ocean ridges, but extends to the vastly larger and older ridge flanks. Upper oceanic crust is porous and permeable due to the presence of lava drainbacks, fissuring, and inter-unit voids, and this porosity and permeability allows active fluid circulation to advect measurable quantities of lithospheric heat from the crust to an average age of 65 Myr. A compilation of crustal porosities shows that this fluid reservoir contains nearly 2% of the total volume of global seawater. Heat flow and sediment thickness data allow calculation of reservoir temperatures, predicting 40°C mean temperatures in Cretaceous crust. Utilizing these temperature estimates, heat flow measurements and models for the thermal structure and evolution of the oceanic lithosphere, we have computed mean hydrothermal fluxes into the deep ocean as a function of plate age. The total hydrothermal volume flux into the oceans approaches 20% of the total riverine input and may contribute to the global seawater mass balance.     

陆相薄互层稠油油藏热采时移地震监测研究

目前,针对陆相薄互层油藏实施时移地震监测尚没有明显的突破,加强时移地震方法、理论的研究和应用,对于我国时移地震监测技术的发展和应用具有积极的意义.根据油藏工程理论,在详细分析火烧油层对储层物性影响的基础上,结合已有的实验室岩石物理测量结果,重点就稠油热采过程中温度的剧烈变化对储层岩石和孔隙流体弹性特性的影响进行了较为详细的分析.并开展了薄互层油藏火烧油层地震正演模拟研究.结果表明,在稠油热采过程中,高温高压不仅会使得孔隙流体特性发生显著的改变,而且储层岩石本身弹性特性的变化也非常明显.即使是对于薄层或薄互层,两方面共同的变化将使得稠油热采前后储层弹性特性产生巨大的差异,进而引起显著的地震异常.因此,开展陆相薄互层稠油热采时移地震监测,不仅具有坚实的岩石物理基础,而且具有切实的可行性.     

Simulation of two-phase fluid flow through compactible reservoirs

Fluid-flow simulators used in the oil industry model the movement of fluids through a porous reservoir rock. These simulators either ignore coupling between the flow and concurring deformation of the solid rock frame or take it into account approximately, in the so-called loose or staggered-in-time mode. In contrast to existing simulators, the one we describe here fully couples two-phase (oil and water) flow to subsurface deformation and simultaneously accounts for all relevant physical phenomena. As such, our flow simulator inherently links time-dependent fluid pressures, saturations, permeabilities and flow velocities to stresses in the whole subsurface. These stresses relate to strains through the non-linear theory of elasticity, allowing us to model time-lapse changes in seismic velocities and anisotropy. The velocity variations manifest themselves in time shifts and reflection amplitudes that are conventionally measured from 4D seismic data. Changes in anisotropy produce time-dependent shear-wave splitting that can be used for monitoring the horizontal stresses.     

利用振幅谱差识别水淹层的变化

储层在长期水驱前后,其中的流体及储层本身会发生变化,导致速度和密度的变化,这一变化为应用四维地震监测油藏的变化提供了依据.本文基于平面波入射水平层状介质,推导出了用振幅谱差（四维振幅谱差）预测油藏变化的公式,并给出数值模拟的结果.理论计算表明,不同程度水驱和不同层位发生水淹其振幅谱差呈现明显的不同,从而在理论上证实了用振幅谱差预测油藏油水变化的可能.这为实际工作中利用振幅谱差预测地下油水变化奠定了理论基础.     

Pore fluid pressure effects in anisotropic rocks: Mechanisms of induced seismicity and weak faults

Many observations and studies indicate that pore fluid pressure in the crustal rocks plays an important role in deformation, faulting, and earthquake processes. Conventional models of pore pressure effects often assume isotropic porous rocks and yield the nondeviatoric pressure effects which seem insufficient to explain diverse phenomena related to pore pressure variation, such as fluid-extraction induced seismicity and crustal weak faults. We derive the anisotropic effective stress law especially for transversely-isotropic and orthotropic rocks, and propose that the deviatoric effects of pore fluid pressure in anisotropic rocks not only affect rock effective strength but also cause variation of shear stresses. Such shear stress variations induced by either pore pressure buildup or pore pressure decline may lead to faulting instability and trigger earthquakes, and provide mechanisms for the failure of crustal weak faults with low level of shear stresses. We believe that the deviatoric effects of pore fluid pressure in anisotropic rocks are of wide application in studies of earthquake precursors and aftershocks, oil and gas reservoir characterization, enhanced oil recovery, and hydraulic fracturing.     

Influence of frequency-dependent anisotropy on seismic amplitude-versus-offset signatures for fractured poroelastic rocks

Frequency-dependent amplitude variation with offset offers an effective method for hydrocarbon detections and analysis of fluid flow during production of oil and natural gas within a fractured reservoir. An appropriate representation for the frequency dependency of seismic amplitude variation with offset signatures should incorporate influences of dispersive and attenuating properties of a reservoir and the layered structure for either isotropic or anisotropic dispersion analysis. In this study, we use an equivalent medium permeated with aligned fractures that simulates frequency-dependent anisotropy, which is sensitive to the filled fluid of fractures. The model, where pores and fractures are filled with two different fluids, considers velocity dispersion and attenuation due to mesoscopic wave-induced fluid flow. We have introduced an improved scheme seamlessly linking rock physics modelling and calculations for frequency-dependent reflection coefficients based on the propagator matrix technique. The modelling scheme is performed in the frequency-slowness domain and can properly incorporate effects of both bedded structure of the reservoir and velocity dispersion quantified with frequency-dependent stiffness. Therefore, for a dispersive and attenuated layered model, seismic signatures represent a combined contribution of impedance contrast, layer thickness, anisotropic dispersion of the fractured media and tuning and interference of thin layers, which has been avoided by current conventional methods. Frequency-dependent amplitude variation with offset responses was studied via considering the influences of fracture fills, layer thicknesses and fracture weaknesses for three classes amplitude variation with offset reservoirs. Modelling results show the applicability of the introduced procedure for interpretations of frequency-dependent seismic anomalies associated with both layered structure and velocity dispersion of an equivalent anisotropic medium. The implications indicate that anisotropic velocity dispersion should be incorporated accurately to obtain enhanced amplitude variation with offset interpretations. The presented frequency-dependent amplitude variation with offset modelling procedure offers a useful tool for fracture fluid detections in an anisotropic dispersive reservoir with layered structures.     

Hot fluid migration: an efficient source of ground deformation: application to the 1982–1985 crisis at Campi Flegrei-Italy

Some solutions of the forced heat advection problem in compressible media are worked out employing a perturbative approach and their implication for thermoelastic deformation are discussed. A sharp temperature front, which migrates at a speed in the order of Darcy flow rate, develops in the medium, giving rise to significant deformation via thermal expansion. A thermally induced pressure source accompanies the temperature front, which may be significant only in very high temperature cases. Results are applied to interpreting the uplift episode of 1982–1984 at Campi Flegrei (C.F.), near Naples, Italy. A mechanism is envisaged for uplift at Campi Flegrei in which a sudden connection is established between a deep, hot, high-pressure fluid reservoir and a shallow, relatively cold, low-pressure aquifer. The inclusion of fluid migration in the deformation model allows simple explanations of several geochemical and geophysical observations made during the bradyseismic crisis. It appears that the proposed mechanism may explain the large observed uplift, without requiring unreasonable pressure increase within the magma chamber. Furthermore, the deformation source may be allowed to be shallower than the magma chamber, as required by any reasonable deformation model at C.F.     

Mn content of reservoir calcite cement: A novel inorganic geotracer of secondary petroleum migration in the tectonically complex Junggar Basin (NW China)

Electronic probe microanalysis(EPMA) results of reservoir calcite cement from fourteen core samples in the Junggar Basin show that Mn-content varies largely between different samples from below the detect limitation to 4.14%,while it displays a generally good correlation with oil-gas shows.This,therefore,likely indicates that concentration of the Mn-content of the calcite cement has a close rela-tion to the intensity of petroleum fluid charging during hydrocarbon secondary migration.In order to assess this hypothesis,oxygen and strontium isotopic measurements on sixteen calcite veins host in source sequences were carried out to investigate the feature of the oil-source petroleum fluid.Analytical results imply that during hydrocarbon generation and migration,deep hot fluid has dissolved volcanic minerals interlined between mudstone source rocks.As Mn is a kind of typical trace element enriched in volcanic rocks,it is reasonable to conclude that the petroleum fluid formed in the source sequences would be Mn-rich.Consequently,calcite cements precipitated from such Mn-rich petroleum fluid would be Mn-rich accordingly.Due to the geologic chromatographic effect during migration along reservoir rocks,the decreasing of the Mn-content of the reservoir calcite cements indicates the migration direction.Then,this novel geotracer was further successfully applied in the study of hydrocarbon migration in the Junggar Basin in combination with organic geochemical analyses during the hydrocarbon migration.The Mn content of the reservoir calcite cement appears promising as a novel inorganic geotracer for the petroleum migration.This paper represents a search for novel indicators of secondary petroleum migration in tectonically complex basins based on fundamentals of the reservoir fluid-rock interactions.     

Research on seismic fluid identification driven by rock physics

Seismic fluid identification works as an effective approach to characterize the fluid feature and distribution of the reservoir underground with seismic data. Rock physics which builds bridge between the elastic parameters and reservoir parameters sets the foundation of seismic fluid identification, which is also a hot topic on the study of quantitative characterization of oil/gas reservoirs. Study on seismic fluid identification driven by rock physics has proved to be rewarding in recognizing the fluid feature and distributed regularity of the oil/gas reservoirs. This paper summarizes the key scientific problems immersed in seismic fluid identification, and emphatically reviews the main progress of seismic fluid identification driven by rock physics domestic and overseas, as well as discusses the opportunities, challenges and future research direction related to seismic fluid identification. Theoretical study and practical application indicate that we should incorporate rock physics, numerical simulation, seismic data processing and seismic inversion together to enhance the precision of seismic fluid identification.     

海底冷泉流体渗漏的原位观测技术及冷泉活动特征

海底冷泉渗漏是海洋环境中广泛分布的自然现象.本文主要介绍海底冷泉原住观测技术及全球几个典型冷来渗漏活动地区的原位观测结果.冷泉渗漏活动主要发育在海底裂隙和断层发育的地区,海底常有泥火山、麻坑、水合物、冷泉生物群落和冷泉碳酸盐岩的发育.冷泉渗漏主要包括液体和气体两种形式,流体流动包括上升流和下降流两种方式,流量变化大,主要与海底环境及构造作用有关.烃类渗漏是以气泡、油滴的形式沿断层向海底喷溢,释放量大,运移过程主要受涌浪、静水压力和裂隙压力等因素控制.冷泉渗漏活动是复杂的水文地质过程,受水合物形成和分解产生的气藏的充注和排空、深部常规油气藏的破坏及海底环境变化等多种作用的控制.