渗流-化学溶解耦合作用下岩石单裂隙渗透特性研究

为揭示在渗流-化学溶解耦合作用下单裂隙渗透特性的变化规律,建立了描述二维渗流-化学溶解耦合作用的偏微分方程组,并利用COMSOL Multiphysics软件成功地求解该方程组。首先,模拟了文献[1]中的盐岩渗流-溶解耦合渗流试验结果,数值模拟结果与试验结果较为吻合,验证了数学模型的正确性和有效性。然后,利用分形理论生成了一个粗糙的裂隙面数字模型,着重分析了二维石灰岩粗糙裂隙面在水流、矿物溶解和输运过程中其渗透特性的变化规律。数值分析显示,（1）溶质浓度对裂隙面的溶解具有非常重要的作用,从而水流进口端的溶解厚度比出口端大得多。（2）裂隙的整体渗透性在初始时刻增加较慢,随着裂隙开度的增大和贯通,溶解速度会逐渐增大,是一个加速的过程。     

Physical and chemical steady-state compaction in deep-sea sediments: Role of mineral reactions

Marine sediments typically exhibit steep porosity gradients in their uppermost centimeters. Although the decrease in porosity with depth below the sediment-water interface is primarily due to compression arising from the accumulation of overlying sediment, early diagenetic mineral dissolution and precipitation reactions may potentially also affect the porosity gradient. Here, we present a steady state compaction model, based on the mass and momentum conservation of total fluid and solid phases, in order to quantify the relative contributions of mineral reactions and physical compaction on porosity changes. The compaction model is applied to estimate hydraulic conductivity and compressive response coefficients of deep-sea sediments from measured porosity depth profiles. The results suggest an inverse relation between the compressive response coefficient and the lithogenic content of marine sediments. For deep-sea sediments exhibiting high rates of dissolution of siliceous shell fragments, the compaction model ignoring mineral reactions overestimates the hydraulic conductivity and compressive response coefficients. In contrast to non-compacting porous media, mineral dissolution in surficial sediments can lead to lower porosity. However, as illustrated for a deep-sea sediment in the equatorial Atlantic characterized by extensive dissolution of calcareous shell fragments, the effect of mineral dissolution and precipitation reactions on porosity gradients is, in most cases, negligible.     

A numerical study of mechanical behavior of a cement paste under mechanical loading and chemical leaching

Based on relevant experimental data of a petroleum cement paste under mechanical loading and chemical leaching, an elastic‐plastic model is first proposed by taking into account plastic shearing and pore collapse. The degradation of mechanical properties induced by the chemical leaching is characterized by a chemical damage variable which is defined as the increase of porosity. Both elastic and plastic properties of the cement paste are affected by the chemical damage. The proposed model is calibrated from and applied to describe mechanical responses in triaxial compression tests respectively on sound and fully leached samples. In the second part, a phenomenological chemical model is defined to establish the relationship between porosity change and calcium dissolution process. The dissolution kinetics is governed by a diffusion law taking into account the variation of diffusion coefficient with calcium concentration. The chemical model is coupled with the mechanical model, and both are applied to describe mechanical response of cement paste samples subjected to progressive chemical leaching and compressive stresses. Comparisons between experimental data and numerical results are presented.     

Diagenesis and Porosity Evolution Model of the Quartzose Sandstone Reservoir with Moderate Palaeogeothermal,Long Period of Large Buried-depth:A Case Form the Lower Shihezi Reservoir in Yan’an,Ordos Basin

This research took Lower Shihezi Reservoir in Yan’an, Ordos Basin as research object. Based on the achievement of former scholars, the article first analyzed burial history and palaeogeothermal history. Furthermore, it analyzed the diagenesis characteristics and divided the diagenetic stages through the application of test data, such as thin section, SEM, cathode luminescence, x-ray diffraction, homogenization temperature of fluid inclusion, etc. Finally, the article analyzed and identified the initial porosity model and various kinds of porosity evolution model in the process of diagenetic, optimized the suitable method of Lower Shihezi Reservoir, established a calculation method that the porosity of reservoir increased in the process of dissolution by applying the effect simulation, simulated porosity evolution of reservoir in research area and analyzed the porosity evolution model. The results showed that the reservoir of quartzose sandstone with moderate palaeogeothermal, long period of large buried-depth, carbonate cement was absent in early phase in the coal measure strata, and the porosity loss by compaction was substantial. Though secondary porosity was generated by dissolution, the quartz and clay mineral were precipitated from dissolution mass, which filled intergranular pores and dissolved pores. In addition to substantial porosity loss by compaction, the filling of massive quartz overgrowth, clay mineral, and ferrocarbonatite, always made such reservoirs to be low or ultra-low permeability reservoirs.     

Temporal alteration of fracture permeability in granite under hydrothermal conditions and its interpretation by coupled chemo-mechanical model

Examining the evolution of fracture permeability under stressed and temperature-elevated conditions, a series of flow-through experiments on a single rock fracture in granite has been conducted under confining pressures of 5 and 10 MPa, under differential water pressures ranging from 0.04 to 0.5 MPa, and at temperatures of 20–90 °C, for several hundred hours in each experiment. Measurements of fluid and dissolved mass fluxes, and post-experimental microscopy, were conducted to constrain the progress of mineral dissolution and/or precipitation and to examine its effect on transport properties. Generally, the fracture aperture monotonically decreased with time at room temperature, and reached a steady state in relatively short periods (i.e., <400 h). However, once the temperature was elevated to 90 °C, the aperture resumed decreasing and kept decreasing throughout the rest of the experimental periods. This reduction may result from the removal of the mineral mass from the bridging asperities within the fracture. Post-experimental observations by scanning electron microscopy, coupled with energy dispersive X-ray spectroscopy (SEM-EDX), revealed the formation of several kinds of secondary minerals such as silica and calcite. However, the precipitated minerals seemed to have had little influence on the flow characteristics within the fracture, because the precipitation was limited to quite local and small areas. The evolving rates and ultimate magnitudes of the fracture aperture are likely to be controlled by the stress exerted over the contacting asperities and temperatures, and by the prescribed flow conditions. Thus, this complex behavior should be attributed to the coupled chemically- and mechanically-induced effect. A coupled chemo–mechano conceptual model, accounting for pressure and free-face dissolutions, is presented in this paper to follow the evolution of the fracture permeability observed in the flow-through experiments. This model addresses the two dissolution processes at the contacting asperities and the free walls within the fractures, and is also capable of describing multi-mineral dissolution behavior. The model shows that the evolution of a fracture aperture (or related permeability) and of element concentrations may be followed with time under arbitrary temperature and pressure conditions. The model predictions for the evolving fracture aperture and elements concentrations show a relatively good agreement with the experimental measurements, although it is not possible to replicate the abrupt reduction observed in the early periods of the experiments, which is likely to be due to an unaccounted mechanism of more stress-mediated fracture compaction driven by the fracturing of the propping asperities.     

OpenGeoSys-ChemApp: a coupled simulator for reactive transport in multiphase systems and application to CO2 storage formation in Northern Germany

Sequestration of CO₂ into a deep geological reservoir causes a complex interaction of different processes such as multiphase flow, phase transition, multicomponent transport, and geochemical reactions between dissolved CO₂ and the mineral matrix of the porous medium. A prognosis of the reservoir behaviour and the feedback from large-scale geochemical alterations require efficient process-based numerical models. For this purpose, the multiphase flow and multicomponent transport code OpenGeoSys-Eclipse have been coupled to the geochemical model ChemApp. The newly developed coupled simulator was successfully verified for correctness and accuracy of the implemented reaction module by benchmarking tests. The code was then applied to assess the impact of geochemical reactions during CO₂ sequestration at a hypothetical but typical Bunter sandstone formation in the Northern German Basin. Injection and spreading of 1.48 × 10⁷ t of CO₂ in an anticline structure of the reservoir were simulated over a period of 20 years of injection plus 80 years of post-injection time. Equilibrium geochemical calculations performed by ChemApp show only a low reactivity to the geochemical system. The increased acidity of the aqueous solution results in dissolution of small amounts of calcite, anhydrite, and quartz. Geochemical alterations of the mineral phase composition result in slight increases in porosity and permeability, which locally may reach up to +0.02 and 0.1 %, respectively.     

Evolution of fluid chemistry in quartz compaction systems: Experimental investigations and numerical modeling

The evolution of fluid chemistry in compacting rock is controlled by coupled chemical processes and rock deformation. In order to characterize this evolution, we conducted water-rock interaction experiments using quartz aggregates at 150 °C and effective pressure of 34.5 MPa. A coupled fluid flow, chemical reaction, and creep compaction model is developed, in which both free-surface reaction and grain-contact dissolution are considered as system volume and porosity evolve.The direct experimental measurement and numerical modeling indicate that effective pressure has significant effects on pore-fluid chemistry. At the early stages of compaction, pore fluids are supersaturated with respect to bulk quartz. With increasing compaction and time, solute concentrations gradually decrease to saturated conditions. Supersaturation is caused mainly by dissolution of ultrafines and high-energy, unstable surfaces which are produced by stress concentrations at grain contacts during the very early stages of compaction. Grain-contact dissolution also contributes to the solute increase in pore fluid in the early stage of compaction, but the effect is small compared to that of ultrafines and unstable surfaces and only slight supersaturation can be produced by it. The gradual decrease in pore-fluid concentration is related to the mechanical removal of ultrafines by pore-fluid flow and the dissolution of ultrafines and unstable surfaces. It also results from the lessening of grain-contact dissolution.Pore fluids in compacting sedimentary basins of quartz sandstone are nearly saturated throughout most of diagenetic processes. Ultrafines and unstable surfaces produced by stress appear not to be the major sources of quartz cement.     

砂岩三轴循环加卸载条件下的渗透率研究

渗透率是地下工程的流-固耦合分析中的一个关键因素。对多孔红砂岩进行了三轴压缩试验,在不同变形阶段实施了轴向应力循环加卸载,并在试验全过程中测量轴向渗透率,得到了试样破坏全过程的渗透率演化规律。从平均应力和循环加卸载对渗透率的影响等两方面进行了深入分析,结果表明,（1）随着轴向变形的增加,初始压密阶段和弹性变形试样渗透率均匀减小;进入塑性变形阶段,渗透率与轴向变形的曲线逐渐趋于水平,低围压条件下渗透率略有增加。（2）轴向加载使骨架颗粒被压缩,引起孔隙减小,造成渗透率减小;采用经验公式定量描述了渗透率和平均应力之间的关系。（3）轴向应力循环加卸载过程中,骨架颗粒的不可恢复变形引起渗透率产生不可恢复现象。（4）峰值后渗透率只发生少许突跳,说明对于多孔砂岩,孔隙和裂隙对渗透率的影响相当,且渗透率的突跳程度随着围压的升高而降低。     

水-岩作用对储层渗透性影响的数值模拟研究——以鄂尔多斯盆地东北部上古生界砂岩储层为例

油气储层的质量及演化往往受沉积作用和成岩作用的共同制约,成岩作用是影响储层孔渗最直接的因素,水-岩反应贯穿于整个成岩作用阶段,决定了孔渗的演化过程。本次研究以鄂尔多斯盆地东北部上古生界陆源碎屑岩储层下石盒子组-山西组为对象,运用多相流反应溶质运移数值模拟的方法,定量研究水-岩化学作用对含油气盆地储层渗透性能的影响,再现水-岩反应所造成的储层矿物溶解、沉淀以及储层渗透性的变化。结果表明：在成岩作用初期,矿物溶解过程占主导地位,模型运行50 Ma时,储层孔隙度由初始的0.30增至0.36,渗透率由初始的5.00×10^-15 m²增至8.70×10^-15 m²;在油气充注阶段（气体CO₂为主）,整个成岩作用矿物的沉淀量大于矿物的溶解量,导致固相体积增加,模型运行20 Ma时,储层孔隙度由0.20降低至约0.15,渗透率由5.00×10^-13 m²下降至约2.00×10^-13 m²。     

Mechanical and transport constitutive models for fractures subject to dissolution and precipitation

Transient changes in the permeability of fractures in systems driven far‐from‐equilibrium are described in terms of proxy roles of stress, temperature and chemistry. The combined effects of stress and temperature are accommodated in the response of asperity bridges where mineral mass is mobilized from the bridge to the surrounding fluid. Mass balance within the fluid accommodates mineral mass either removed from the flow system by precipitation or advection, or augmented by either dissolution or advection. Where the system is hydraulically closed and initially at equilibrium, reduction in aperture driven by the effects of applied stresses and temperatures will be augmented by precipitation on the fracture walls. Where the system is open, the initial drop in aperture may continue, and accelerate, where the influent fluid is oversaturated with respect to the equilibrium mineral concentration within the fluid, or may reverse, if undersaturated. This simple zero‐dimensional model is capable of representing the intricate behavior observed in experiments where the feasibility of fracture sealing concurrent with net dissolution is observed. This zero‐order model is developed as a constitutive model capable of representing key aspects of changes in the transport parameters of the continuum response of fractured media to changes in stress, temperature and chemistry. Copyright © 2009 John Wiley & Sons, Ltd.     

Coupled flow and geomechanical processes during enhanced coal seam methane recovery through CO2 sequestration

The sensitivity of coal permeability to the effective stress means that changes in stress as well as pore pressure within a coal seam lead to changes in permeability. In addition coal swells with gas adsorption and shrinks with desorption; these sorption strains impact on the coal stress state and thus the permeability. Therefore the consideration of gas migration in coal requires an appreciation of the coupled geomechanical behaviour. A number of approaches to representing coal permeability incorporate the geomechanical response and have found widespread use in reservoir simulation. However these approaches are based on two simplifying assumptions; uniaxial strain (i.e. zero strain in the horizontal plane) and constant vertical stress. This paper investigates the accuracy of these assumptions for reservoir simulation of enhanced coalbed methane through CO₂ sequestration. A coupled simulation approach is used where the coalbed methane simulator SIMED II is coupled with the geomechanical model FLAC3D. This model is applied to three simulation case studies assembled from information presented in the literature. Two of these are for 100% CO₂ injection, while the final example is where a flue gas (12.5% CO₂ and 87.5% N₂) is injected. It was found that the horizontal contrast in sorption strain within the coal seam caused by spatial differences in the total gas content leads to vertical stress variation. Thus the permeability calculated from the coupled simulation and that using an existing coal permeability model, the Shi–Durucan model, are significantly different; for the region in the vicinity of the production well the coupled permeability is greater than the Shi–Durucan model. In the vicinity of the injection well the permeability is less than that calculated using the Shi–Durucan model. This response is a function of the magnitude of the strain contrast within the seam and dissipates as these contrasts diminish.     

A coupled compressible flow and geomechanics model for dynamic fracture aperture during carbon sequestration in coal

This paper presents the development of a discrete fracture model of fully coupled compressible fluid flow, adsorption and geomechanics to investigate the dynamic behaviour of fractures in coal. The model is applied in the study of geological carbon dioxide sequestration and differs from the dual porosity model developed in our previous work, with fractures now represented explicitly using lower-dimensional interface elements. The model consists of the fracture-matrix fluid transport model, the matrix deformation model and the stress-strain model for fracture deformation. A sequential implicit numerical method based on Galerkin finite element is employed to numerically solve the coupled governing equations, and verification is completed using published solutions as benchmarks. To explore the dynamic behaviour of fractures for understanding the process of carbon sequestration in coal, the model is used to investigate the effects of gas injection pressure and composition, adsorption and matrix permeability on the dynamic behaviour of fractures. The numerical results indicate that injecting nonadsorbing gas causes a monotonic increase in fracture aperture; however, the evolution of fracture aperture due to gas adsorption is complex due to the swelling-induced transition from local swelling to macro swelling. The change of fracture aperture is mainly controlled by the normal stress acting on the fracture surface. The fracture aperture initially increases for smaller matrix permeability and then declines after reaching a maximum value. When the local swelling becomes global, fracture aperture starts to rebound. However, when the matrix permeability is larger, the fracture aperture decreases before recovering to a higher value and remaining constant. Gas mixtures containing more carbon dioxide lead to larger closure of fracture aperture compared with those containing more nitrogen.     

Prediction of dilation and permeability changes in rock salt

A numerical procedure has been developed for predicting dilation (porosity) and gas permeability changes in rock salt. The hierarchical single-surface constitutive model of Desai and co-workers is used a finite element program to calculate the state of stress and strain surrounding excavations in rock salt. The elastoplastic constitutive model accounts for strain hardening, a non-associative volumetric response and stress-path-dependent behaviour. The calculated stress and strain fields are used in a flow model based on the equivalent channel concept to predict permeability. Parameters for both the mechanical and permeability models are developed from laboratory test results. Two field experiments adjacent to underground excavations are modelled. The extent of the dilated rock zone around the excavation is predicted well, but the magnitude of the porosity and gas permeability is underpredicted very near the excavations. This discrepancy is attributed to model parameters derived from loading-only laboratory tests, whereas significant unloading occurs in the field. The shape of the yield surface was found to be an important factor in dilation and permeability predictions. Similar stress, strain and permeability fields were obtained with different model types (plane strain or axisymmetric) and initial stress states, and with instantaneous and progressive excavation.     

Effects of medium and pore‐fluid compressibility on chemical‐dissolution front instability in fluid‐saturated porous media

Theoretical analysis and computational simulations have been carried out to investigate how medium and pore‐fluid compressibility affects the chemical‐dissolution front propagation, which is associated with a fully‐coupled nonlinear problem between porosity, pore‐fluid pressure, pore‐fluid density and reactive chemical‐species transport within a deformable and fluid‐saturated porous medium. When the fully‐coupled nonlinear system is in a subcritical state, some analytical solutions have been derived for a special case, in which the ratio of the equilibrium concentration to the solid molar density of the chemical species is approaching zero. To investigate the effect of either medium compressibility or pore‐fluid compressibility on the evolutions of chemical dissolution fronts in supercritical chemical dissolution systems, numerical algorithms and procedures have been also proposed. The related theoretical and numerical results have demonstrated that: (i) not only can pore‐fluid compressibility affect the propagating speeds of chemical dissolution fronts in both subcritical and supercritical systems, but also it can affect the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; (ii) medium compressibility may have a little influence on the propagating speeds of chemical dissolution fronts, but it can have significant effects on the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; and (iii) both medium and pore‐fluid compressibility may stabilize irregular chemical‐dissolution‐fronts in supercritical chemical dissolution systems. Copyright © 2011 John Wiley & Sons, Ltd.     

The impact of diagenetic fluid–rock reactions on Rotliegend sandstone composition and petrophysical properties (Altmark area, central Germany)

In the framework of the German R&D joint project CLEAN (CO₂ large-scale enhanced gas recovery in the Altmark natural gas field), Rotliegend reservoir sandstones of the Altensalzwedel block in the Altmark area (Saxony-Anhalt, central Germany) have been studied to characterise litho- and diagenetic facies, mineral content, geochemical composition, and petrophysical properties. These sands have been deposited in a playa environment dominated by aeolian dunes, dry to wet sand flats and fluvial channel fills. The sediments exhibit distinct mineralogical, geochemical, and petrophysical features related to litho- and diagenetic facies types. In sandstones of the damp to wet sandflats, their pristine red colours are preserved and porosity and permeability are only low. Rocks of the aeolian environment and most of the channel fill deposits are preferentially bleached and exhibit moderate to high porosity and permeability. Although geochemical element whole rock content in these rocks is very similar, element correlations are different. Variations in porosity and permeability are mainly due to calcite and anhydrite dissolution and differences in clay coatings with Fe-bearing illitic-chloritic composition exposed to the pore space. Moreover, mineral dissolution patterns as well as compositions (of clays and carbonate) and morphotypes of authigenic minerals (chlorite, illite) are different in red and bleached rocks. Comparison of the geochemical composition and mineralogical features of diagenetically altered sandstones and samples exposed to CO₂-bearing fluids in laboratory batch experiments exhibit similar character. Experiments prove an increase in wettability and water binding capacity during CO₂ impact.     

Effects of groundwater flow on mineral diagenesis, with emphasis on carbonate aquifers

 This article provides a critical synopsis of the effects of groundwater flow on mineral diagenesis. Emphasis is placed on those aspects and processes that change porosity and permeability in carbonate aquifers, because they are of particular importance to human societies as sources of supplies of water for human consumption (drinking, irrigation) and of crude oil and natural gas. Diagenetic settings in carbonates as well as clastics are generally ill defined. This paper proposes a new comprehensive classification of diagenetic settings into near-surface, shallow-, intermediate-, and deep-burial diagenetic settings; hydrocarbon-contaminated plumes; and fractures. These settings are defined on the basis of mineralogy, petroleum, hydrogeochemistry, and hydrogeology. This classification is applicable to all sedimentary basins. Diagenesis is governed by various intrinsic and extrinsic factors that include thermodynamic and kinetic constraints, as well as microstructural factors that may override the others. These factors govern diagenetic processes, such as dissolution, compaction, recrystallization, replacement, and sulfate–hydrocarbon redox-reactions. Processes such as cementation, dissolution, and dolomitization require significant flow of groundwater driven by an externally imposed hydraulic gradient. Other processes, such as stylolitization and thermochemical sulfate reduction, commonly take place without significant groundwater flow in hydrologically nearly or completely stagnant systems that are geochemically "closed." Two major effects of groundwater flow on mineral diagenesis are enhancement and reduction of porosity and permeability, although groundwater flow can also leave these rock properties essentially unchanged. In extreme cases, an aquifer or hydrocarbon reservoir rock can have highly enhanced porosity and permeability due to extensive mineral dissolution, or it can be plugged up due to extensive mineral precipitation. Received, April 1998 · Revised, July 1998 · Accepted, September 1998     

煤层瓦斯卸压抽放动态过程的气-固耦合模型研究

根据瓦斯渗流与煤体变形的基本理论,引入煤体孔隙变形与透气性演化的耦合作用方程,建立了考虑煤层吸附、解吸作用的含瓦斯煤岩固-气耦合作用模型。应用该模型模拟研究了不同压力影响下瓦斯抽放过程中煤层透气性的演化和抽放孔周围瓦斯压力的变化规律,其结果对于深入理解瓦斯抽放作用机制并采取相应的瓦斯预防和控制措施等具有重要的理论和实践意义。     

Numerical simulation of carbon dioxide migration in a sandstone aquifer considering the fluid–solid coupling

Very limited investigations have been done on the numerical simulation of carbon dioxide (CO₂) migration in sandstone aquifers taking consideration of the interactions between fluid flow and rock stress. Based on the poroelasticity theory and multiphase flow theory, this study establishes a mathematical model to describe CO₂ migration, coupling the flow and stress fields. Both finite difference method (FDM) and finite element method (FEM) were used to discretize the mathematical model and generate a numerical model. A case study was carried out using the numerical model on the Jiangling sandstone aquifer in the Jianghan basin, China. The rock mechanics parameters of reservoir and overlying strata of Jiangling depression were obtained by triaxial tests. A two-dimensional model was then built to simulate carbon dioxide migration in the sandstone aquifer. The numerical simulation analyzes the carbon dioxide migration distribution rule with and without considering capillary pressure. Time-dependent migration of CO₂ in the sandstone aquifer was analyzed, and the result from the coupled model was compared with that from a traditional non-coupled model. The calculation result indicates a good consistency between the coupled model and the non-coupled model. At the injection point, the CO₂ saturation given by the coupled model is 15.39 % higher than that given by the non-coupled model; while the pore pressure given by the coupled model is 4.8 % lower than that given by the non-coupled model. Therefore, it is necessary to consider the coupling of flow and stress fields while simulating CO₂ migration for CO₂ disposal in sandstone aquifers. The result from the coupled model was also sensitized to several parameters including reservoir permeability, porosity, and CO₂ injection rate. Sensitivity analyses show that CO₂ saturation is increased non-linearly with CO₂ injection rate and decreased non-linearly with reservoir porosity. Pore pressure is decreased non-linearly with reservoir porosity and permeability, and increased non-linearly with CO₂ injection rate. When the capillary pressure was considered, the computed gas saturation of carbon dioxide was increased by 10.75 % and the pore pressure was reduced by 0.615 %.     

Stress-Dependent Fluid Flow and Permeability in Fractured Media: from Lab Experiments to Engineering Applications

Summary. Permeability is a physical property in rocks of extreme importance in energy engineering, civil and environmental engineering, and various areas of geology. Early on, fractures in fluid flow models were assumed to be rigid. However, experimental research and field data confirmed that stress-deformation behavior in fractures is a key factor governing their permeability tensor. Although extensive research was conducted in the past, the three-dimensional stress-permeability relationships, particularly in the inelastic deformation stage, still remain unclear. In this paper, laboratory experiments conducted on large concrete blocks with randomly distributed fractures and rock core samples are reported to investigate fluid flow and permeability variations under uniaxial, biaxial and triaxial complete stress-strain process. Experimental relationships among flowrate, permeability and fracture aperture in the fractured media are investigated. Results show that the flowrate and stress/aperture exhibit “cubic law” relationship for the randomly distributed fractures. A permeability-aperture relationship is proposed according to the experimental results. Based on this relationship, stress-dependent permeability in a set of fractures is derived in a three-dimensional domain by using a coupled stress and matrix-fracture interactive model. A double porosity finite element model is extended by incorporating such stress-dependent permeability effects. The proposed model is applied to examine permeability variations induced by stress redistributions for an inclined borehole excavated in a naturally fractured formation. The results indicate that permeability around underground openings depends strongly on stress changes and orientations of the natural fractures.     

渝东南下志留统龙马溪组页岩储集层成岩作用及其对孔隙发育的影响*

以重庆綦江、南川和涪陵地区的下志留统龙马溪组下部页岩样品为例,通过X射线衍射(XRD)、场发射扫描电镜、TOC测试、R_O测试、高压压汞和低温氮吸附等手段,对研究区龙马溪组下部页岩储集层的矿物成分、有机组分、孔隙特征及成岩作用进行研究。实验结果表明,龙马溪组页岩成岩作用类型主要包括机械压实作用、化学压实作用、溶蚀作用、交代作用、破裂作用及矿物转化等,演化阶段介于晚成岩作用与极低级变质作用之间;成岩作用对页岩储集层的孔渗特征、孔隙成因类型、非均质性和形貌特征均有显著影响。对成岩作用与孔隙演化进行了匹配,认为机械/化学压实作用、溶蚀作用及过程中的有机生烃、矿物转化对孔隙发育影响最为显著,并将储集层成岩作用划分为溶蚀作用等建设性成岩作用与压实作用等破坏性成岩作用。