Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Sedimentary characteristics of tide-dominated estuary in Donghetang Formation（Upper Devonian）, central Tarim Basin

The Donghetang Formation （Upper Devonian） in central Tarim Basin has been thought an important oil and gas reservoir since the abundant oil and gas resources were found in the wells W16, W20, W34, and other fields. However, the sedimentary environment of the Donghetang Formation has been disputed because it suffered from both tidal and fluvial actions and there were not rich fossils in the sandstone. After the authors analyzed sedimentary features by means of drill cores, well logging data, paleosols, and with SEM obseruations, three kinds of sedimentary environments were distinguished： alluvial fan, tide-dominated estuary, and shelf. Particularly, the sedimentary features of tide-dominated estuary were studied in detail. Besides, the authors discussed sedimentary characteristics of the Donghetang Formation which was divided into two fourth-order sequences and five system tracts. At the same time, according to the forming process of five system tracts, the whole vertical evolution and lateral transition of tide-dominated estuary were illustrated clearly. Finally, the reservoir quality was evaluated based on porosity and permeability.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

鄂尔多斯盆地西峰地区长8砂岩微观非均质性的实验分析

通过真实砂岩微观模型实验,客观地分析了鄂尔多斯盆地西峰地区长8砂岩的微观非均质性特征,发现该区长8储层为致密型砂岩,直观地显示出长8砂岩储层具有很强的微观非均质性,砂岩的成岩作用和孔隙结构是影响鄂尔多斯盆地西峰地区长8砂岩储层微观非均质性的主要原因。在低渗透砂岩中,引起砂岩产生强烈的非均质性,而这些砂岩则形成长8砂岩储层的高渗带。     

新疆阿尔泰蒙库铁矿床的成矿流体及成矿作用

蒙库大型铁矿床赋存于上志留统—下泥盆统康布铁堡组变质火山-沉积岩系中,容矿岩石为石榴子石矽卡岩、变粒岩、浅粒岩和大理岩。矿体总体顺层分布,空间上与矽卡岩密切相关。研究表明,矽卡岩期石榴子石以发育玻璃质熔融包裹体、流体熔融包裹体和流体包裹体为特征,晚期矽卡岩阶段矿物中发育液相包裹体,变质期矿物中主要发育液相包裹体和含子矿物包裹体。矽卡岩期熔融包裹体的均一温度为1100℃,早期矽卡岩阶段流体包裹体均一温度变化于193～499℃,在450℃、350℃和230℃出现峰值。中期矽卡岩阶段均一温度变化于236～550℃,峰值为350℃。区域变质期均一温度介于132～513℃,在350℃、230℃和190℃出现峰值。流体包裹体的盐度w(NaCleq)介于1.23%～60.31%,流体密度变化于0.60～1.16g/cm3。石榴子石、石英和方解石的δ18OSMOW变化于0.2‰～8.4‰,δ18OH2O介于-5.1‰～5.33‰,δD为-127‰～-81‰,表明矽卡岩期成矿流体主要是岩浆水,混合少量大气降水;变质期流体主要为大气降水,为混合变质水。方解石δ13CPDB变化于-6.1‰～-2.3‰,表明流体中碳来自深部或地幔。成矿时代为早泥盆世早期(略晚于404～400Ma),成矿作用与矽卡岩的退化变质作用有关。     

高、低煤阶煤层气藏地质特征及控气作用差异性研究

高、低煤阶煤层气地质特征及控气作用差异性是研究煤层气富集成藏的重要组成部分,是煤层气勘探开发理论研究过程中重要的基础性研究领域之一。本文以中国沁水、阜新盆地和美国粉河盆地等典型的含气盆地为例,探讨了高、低煤阶煤层气的储层物性差异,分析了构造控气和水文地质控气作用的差异性。研究表明,高煤阶气藏含气量高,CH4百分含量高,δ13C1值大于-38.75‰,储层渗透率变化小,储层改造难,构造热事件对煤层气的生成、富集贡献大,持续的水动力使气藏遭到破坏,且破坏幅度大,现今地下水格局对气藏的形成具有一定的影响;低煤阶气藏含气量低,CH4百分含量低,δ13C1值大于-49.11‰,储层渗透率变化大,储层易改造,煤热演化史及煤阶影响着煤层气的生成、富集,在煤层气生成过程中活跃的水动力是甲烷生成的主要的水文地质条件之一,但持续的水动力使气藏遭到破坏,且破坏幅度小,而合适的地层水矿化度则是低煤阶煤层气生成的重要条件,地下水格局对气藏的调整和改造起到决定性的影响。