首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Natural evaporite dissolution in the subsurface can lead to cavities having critical dimensions in the sense of mechanical stability. Geomechanical effects may be significant for people and infrastructures because the underground dissolution may lead to subsidence or collapse (sinkholes). The knowledge of the cavity evolution in space and time is thus crucial in many cases. In this paper, we describe the use of a local nonequilibrium diffuse interface model for solving dissolution problems involving multimoving interfaces within three phases, that is, solid–liquid–gas as found in superficial aquifers and karsts. This paper generalizes developments achieved in the fluid–solid case, that is, the saturated case [1]. On one hand, a local nonequilibrium dissolution porous medium theory allows to describe the solid–liquid interface as a diffuse layer characterized by the evolution of a phase indicator (e.g., porosity). On the other hand, the liquid–gas interface evolution is computed using a classical porous medium two‐phase flow model involving a phase saturation, that is, generalized Darcy's laws. Such a diffuse interface model formulation is suitable for the implementation of a finite element or finite volume numerical model on a fixed grid without an explicit treatment of the interface movement. A numerical model has been implemented using a finite volume formulation with adaptive meshing (e.g., adaptive mesh refinement), which improves significantly the computational efficiency and accuracy because fine gridding may be attached to the dissolution front. Finally, some examples of three‐phase dissolution problems including density effects are also provided to illustrate the interest of the proposed theoretical and numerical framework. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

2.
The aging of gypsum in underground mines   总被引：1，自引：0，他引：1  
The aging properties of gypsum extracted from two underground mines in France (Livry–Gargan and Grozon) have been highlighted by means of observation performed using scanning electron microscopy on samples taken along horizontal boreholes drilled through to the middle of several pillars. The aging process is exhibited by the presence of traces of dissolution (edges of the dissolved crystals, corroded crystalline surfaces, important intra- and intercrystalline space). These dissolution figures decrease in number and in intensity from the wall heading towards the middle of the pillar. The “older” pillars display a greater number of dissolution traces than the “more recent” pillars. Then, the role of the aging has been assessed by determining various physical and mechanical parameters on cylindrical samples of 38 mm in diameter and 76 mm in height, machined from the set of drilled boreholes. The following parameters have been quantified: density, grains density, total porosity, porosity accessible to water, intrinsic permeability to nitrogen, velocity of ultrasonic compression and shear waves, dynamic Young's modulus, dynamic Poisson's ratio, uniaxial compressive strength and static Young's modulus. Significant variations in parameter values between the wall and the middle of the pillar were recorded and recognized as being strongly correlated with the intensity of dissolution traces, and hence with gypsum aging.  相似文献   

3.
Authigenic gypsum was found in a gravity core, retrieved from the top of Mound Perseverance, a giant cold‐water coral mound in the Porcupine Basin, off Ireland. The occurrence of gypsum in such an environment is intriguing, because gypsum, a classic evaporitic mineral, is undersaturated with respect to sea water. Sedimentological, petrographic and isotopic evidence point to diagenetic formation of the gypsum, tied to oxidation of sedimentary sulphide minerals (i.e. pyrite). This oxidation is attributed to a phase of increased bottom currents which caused erosion and enhanced inflow of oxidizing fluids into the mound sediments. The oxidation of pyrite produced acidity, causing carbonate dissolution and subsequently leading to pore‐water oversaturation with respect to gypsum and dolomite. Calculations based on the isotopic compositions of gypsum and pyrite reveal that between 21·6% and 28·6% of the sulphate incorporated into the gypsum derived from pyrite oxidation. The dissolution of carbonate increased the porosity in the affected sediment layer but promoted lithification of the sediments at the sediment‐water interface. Thus, authigenic gypsum can serve as a signature for diagenetic oxidation events in carbonate‐rich sediments. These observations demonstrate that fluid flow, steered by environmental factors, has an important effect on the diagenesis of coral mounds.  相似文献   

4.
Quartz dissolution in melts in the KCAS and KCMAS systems results in the formation of a silica- and potassium-enriched boundary layer next to the dissolving crystals. The presence of potassium in CAS melts has no discernible effect on dissolution rate compared with that in K-free melts with otherwise similar composition despite a small decrease in the diffusivity of silica in the potassium-bearing melts. The decrease in silica diffusivity is offset by an increase in the solubility of silica in the K-bearing melts. Addition of potassium to CMAS melts results in a large decrease in the dissolution rate of quartz. Even though the solubility of silica is enhanced, the addition of potassium leads to large changes in the structure of the melt in the boundary layer (as measured by NBO/T), which results in a large decrease in the diffusivity of silica and thus slower dissolution. There is significant diffusive coupling of Al2O3, CaO and MgO during dissolution, which leads to local uphill diffusion of these components. K2O is decoupled from the other components, as shown by its much thicker diffusion zone. Potassium moves through the boundary layer as a result of two homogeneous reactions: uphill diffusion in which potassium diffuses into the silica-enriched melt adjacent to the dissolving quartz crystal and downhill diffusion in the region furthest from the crystal–melt interface where SiO2 and K2O diffuse away from the interface together.  相似文献   

5.
Convective mixing of dissolved carbon dioxide (CO2) with formation brine has been shown to be a significant factor for the rate of dissolution of CO2 and thus for determining the viability of geological CO2 storage sites. In most previous convection investigations, a no-flow boundary condition was used to represent the interface between an upper region with CO2 and brine and the single-phase brine region beneath. However, due to interfacial tension between the phases, the water phase is partly mobile in the upper region and advection may occur. Based on linear stability analysis and numerical simulations, we show that advection across the interface leads to considerable destabilization of the system. In particular, the time of onset of instability is reduced by a factor of two and the rate of dissolution is enhanced by a factor of two for three of four formations we consider, and by 40 % for the fourth formation. It is found that exponential decay of the relative permeability away from the interface provides a useful approximation to the real system. In addition, the exponential decay also simplifies the linear stability analysis. Interestingly, formations with large absolute permeability and small porosity have the largest impact from the transition zone, despite the fact that the relative permeability decays quickly above the interface in these formations. This is because the length-scale of instability is smallest in these formations.  相似文献   

6.
In this paper, a local non-equilibrium diffuse interface model is introduced for describing solid–liquid dissolution problems. The model is developed based on the analysis of Golfier et al. (J Fluid Mech 457:213–254, 2002) upon the dissolution of a porous domain, with the additional requirement that density variations with the mass fraction are taken into account. The control equations are generated by the upscaling of the balance equations for a solid–liquid dissolution using a volume averaging theory. This results into a diffuse interface model (DIM) that does not require an explicit treatment of the dissolving interface, e.g., the use of arbitrary Lagrangian–Eulerian (ALE) methods, for instance. Test cases were performed to study the features and influences of the effective coefficients inside the DIM. In particular, an optimum expression for the solid–liquid exchange coefficient is obtained from a comparison with the referenced solution by ALE simulations. Finally, a Ra–Pe diagram illustrates the interaction of natural convection and forced convection in the dissolution problem.  相似文献   

7.
Scallops and flutes are common dissolution rock forms encountered in karst caves and surface streams. Their evolution is only partially understood and no numerical model that simulates their formation has been presented. This work at least partially fills the gap by introducing a numerical approach to simulate the evolution of different initial forms of soluble surfaces embedded in a turbulent fluid. The aim is to analyze wall dissolution phenomena from basic principles and to identify stable profiles. The analysis is based on a finite volume moving boundary method. The underlying mathematical model is a \(k-\epsilon \) turbulent model for fluid flow coupled with turbulent scalar transport. The rock wall is treated as a moving boundary, where the normal wall retreat velocity is proportional to the under-saturation of the boundary fluid cells with respect to the mineral comprising the wall. As the flow time scale is several orders of magnitude smaller than the dissolution time scale, stationary flow field, concentration field and wall propagation velocity are calculated for each iteration. The boundary at all points is then moved by distracting minimal velocity along the entire boundary from the actual velocity at a certain location, and then normalized to the maximum allowed shift, which is equal to half the height of the boundary cell. In this way only deformation of the initial wall is calculated. The method was applied to several different initial profiles. During the evolution, the profiles progressively converged towards stable forms. In this work, a framework is proposed for a computation of the moving boundary problem related to slow dissolution of a soluble surface.  相似文献   

8.
A Lagrangian numerical approach for the simulation of rapid landslide runouts is presented and discussed. The simulation approach is based on the so‐called Particle Finite Element Method. The moving soil mass is assumed to obey a rigid‐viscoplastic, non‐dilatant Drucker–Prager constitutive law, which is cast in the form of a regularized, pressure‐sensitive Bingham model. Unlike in classical formulations of computational fluid mechanics, where no‐slip boundary conditions are assumed, basal slip boundary conditions are introduced to account for the specific nature of the landslide‐basal surface interface. The basal slip conditions are formulated in the form of modified Navier boundary conditions, with a pressure‐sensitive threshold. A special mixed Eulerian–Lagrangian formulation is used for the elements on the basal interface to accommodate the new slip conditions into the Particle Finite Element Method framework. To avoid inconsistencies in the presence of complex shapes of the basal surface, the no‐flux condition through the basal surface is relaxed using a penalty approach. The proposed model is validated by simulating both laboratory tests and a real large‐scale problem, and the critical role of the basal slip is elucidated. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

9.
A numerical model of coupled saturated/unsaturated water flow, heat transfer and multi-component reactive solute transport is presented to evaluate the long-term geochemical evolution in bentonite, concrete and clay formation for a potential geological radioactive waste repository. Changes in formation porosity caused by mineral dissolution/precipitation reactions are taken into account. Simulations were carried out with a general-purpose multicomponent reactive transport code, CORE2D V4. Numerical results show that pH in the bentonite porewater can vary from neutral to up to 13 over a time scale of 1 Ma although dissolution of silica minerals and precipitation of secondary calcium silicate hydrate (CSH) minerals in bentonite buffer the effect of the hyperalkaline plume. Mineral precipitation reduces the volume of pore space in bentonite close to the bentonite–concrete interface due to the precipitation of CSH minerals. Model results indicate that bentonite porosity decreases less than 25%. The hyperalkaline plume from the concrete only extends to a distance of 0.7 m in the clay formation over the time range of 1 Ma.  相似文献   

10.
The widespread and dissected nature of the Angolan gypsiferous salt residuals offers a uniquely detailed view of the lateral and vertical relations inherent to secondary evaporite textures, which typify exhumed salt masses worldwide. Such secondary textures are sometimes misinterpreted as primary evaporite textures. Thin, metre‐scale and patchy, dome‐like gypsum accumulations are well‐exposed within strongly incised present‐day river valleys along the eastern margin of the Namibe and Benguela basins (south‐west Angola). These sections are time equivalent to the main basinward subsurface evaporites (Aptian Loeme Formation) which mostly consist of halite. The gypsum (here called the Bambata Formation) is interpreted to represent the final residual product of fractional dissolution and recrystallization of the halite mass that occurred during Late Cretaceous margin uplift and continues today. This halite underwent multiple episodes of diagenetic alteration between its deposition and its final exhumation, leading to the formation of various secondary gypsum fabrics and solution‐related karst and breccia textures that typify the current evaporite outcrop. Four different diagenetic gypsum fabrics are defined: thinly bedded alabastrine, nodular alabastrine, displacive selenite rosettes and fibrous satin‐spar gypsum. Current arid conditions are responsible for a thin weathered crust developed at the top of the outcropping gypsum, but the fabrics in the main core of the current at‐surface evaporite unit mostly formed during the telogenetic stage of uplift prior to complete subaerial exposure. Alteration occurred as various dissolving and rehydrating saline minerals encountered shallow aquifers in the active phreatic and vadose zones. Geomorphological and petrographic analyses, mostly based on the cross‐cutting relations and crystallographic patterns in the outcrop, are used to propose a sequence of formation of these different fabrics.  相似文献   

11.
This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

12.
Mineral dissolution is an important factor in many magmatic processes such as melting, assimilation and magma mixing. Since it is not possible to determine dissolution rates or mechanisms from natural samples, experimental measurements are very useful. However, the geometry of the crystal–melt system can have a large effect on the measured rate, depending on whether the contaminated melt formed during dissolution is gravitationally stable or unstable. This study examines the effects of the crystal–melt geometry on the dissolution rate and mechanism. The experiments were performed using basanite melt and cylinders and spheres prepared from a single crystal of natural quartz. All of the experiments were performed in the piston cylinder apparatus at 0.5 GPa and 1350 °C. Four crystal–melt geometries were used: (1) quartz cylinders on top of a column of melt; (2) quartz cylinders beneath a column of basanite melt; (3) quartz cylinders in the middle of column of melt; (4) quartz spheres on top of a column of basanite melt. These geometries allow an examination of non-convective, convective and mixed non-convective/convective dissolution. Sphere experiments were included, as this has been the most commonly used geometry in previous experimental studies. In all of the experiments quartz dissolves directly into the basanite without formation of cristobalite or tridymite. Quartz on top of a column of melt dissolves at a rate almost proportional to the square root of time and forms a silica-rich compositional boundary layer that is gravitationally stable. All of the samples show well-defined compositional gradients in the boundary layer; however, the melt at the interface varies in composition with time and plots of concentration as a function of distance normalized to time show that the diffusion rate of SiO2 increases with time. These data suggest that the rate-controlling step during quartz dissolution is interface reaction rather than cation diffusion. Quartz on the bottom of a column of basanite dissolves much more quickly than in the quartz-on-top experiments and the dissolution rate is linear, due to the periodic gravitational instability and resultant convection of the boundary layer. Even though interface kinetics are the rate-controlling step in quartz dissolution, convection causes an increase in dissolution rate because it replenishes the boundary layer with new, silica-undersaturated melt, which dissolves the quartz more quickly than the contaminated melt. These data suggest that the interface reaction rate is controlled by the degree of undersaturation of the solvent melt in the dissolving component. Both quartz-in-middle and quartz sphere experiments dissolve at a rate intermediate between the two extremes and both show a power law rate. Both dissolve by a combination of convective and non-convective dissolution but the sphere experiments are affected by an additional factor. During the experiment the sphere can sink through the capsule causing forced convection which adds another complication to the interpretation of the dissolution rate data. The results of this study indicate that the choice of experiment geometry plays a major role in determining the observed dissolution rate. Mineral spheres, which have been widely used in the past, are not ideal for dissolution studies. Instead, dissolution rates and mechanisms are best determined in the absence of convection. These experiments have an additional advantage in that for diffusion-controlled dissolution, they allow determination of cation diffusivity. Received: 2 March 2000 / Accepted: 11 April 2000  相似文献   

13.
何明倩  黄文辉  久博 《地学前缘》2021,28(4):327-336
膏质白云岩是鄂尔多斯盆地中东部下古生界马家沟组M51最有利储层,影响膏质白云岩储层性质的最大因素在于溶蚀作用。本文在已有的研究基础上,应用岩心资料及阴极发光、地球化学等方法,测定M51膏质白云岩及其特殊的示底构造的化学成分,并分析其成因模式。岩心、薄片以及物性资料显示:石膏溶蚀形成的大量膏模孔及伴生的次生孔是储层的主要储集空间,示底构造规模性出现于膏质云岩中。阴极发光、常微量元素分析结果显示:示底构造上部充填埋藏期方解石;下部充填碎屑白云石,与示底构造外围的准同生期泥-粉晶白云石具有相近的地球化学特征,与深部埋藏期形成的细-中晶白云石差异明显。示底构造的演化是膏质白云岩有利储层形成的重要依据,本文依据化学离子平衡原理,对不同成岩阶段发生的溶蚀作用、填充作用进行离子间相互置换的重演,分析得知不稳定的化学组成是石膏易被溶蚀的内在原因,石膏的规模性溶蚀是膏质云岩有利储层形成的关键。  相似文献   

14.
This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

15.
The stability of underground mines represents a key issue for active and abandoned mines. Over the last few years, several collapses of underground mines in France have affected existing buildings and infrastructures. Many factors are generally identified as the cause of failures: pillar ageing, fractures, and pillars’ height to width ratio, etc. Among the treatment techniques available to prevent instability and reduce the deterioration of pillars, backfill is the most frequently used. A research programme, supported by the French Ministry of the Environment, was developed to study the operability of partial and total backfill using waste material in the Livry-Gargan gypsum mine (near Paris, France), where pillar height is 17 m. The paper focuses on: (1) the characterisation of the gypsum and fill material (laboratory and in situ tests), (2) the in situ measurements, involving 5 pillars equipped with 19 pressure cells, since 1999, (3) and numerical modelling of fractured pillars performed in order to improve understanding of the effects of backfill on the stability of room-and-pillar mines. The study clearly shows the operability and the advantages of partial and total backfill for short-term pillar stability. The induced horizontal pressure generated by backfill can reach 200 kPa. The use of numerical modelling also shows the effect of backfill on fractures and that backfill reduces indicatively the shear displacement and the opening of fractures. Numerical modelling helps in identifying the mechanisms of backfill and in a better understanding of the behaviour of backfilled mines.  相似文献   

16.
露天地下开采隔离层稳定性研究   总被引：2，自引：0，他引：2  
岩小明  李夕兵  郭雷  高峰 《岩土力学》2007,28(8):1682-1686
露天地下开采隔离层稳定性分析是矿山中经常遇到的问题。以大宝山矿露天地下开采的工程实例,采用数值分析方法,对隔离层稳定性进行分析计算,发现空区顶板的拉应力是关系隔离层稳定性的关键因素。通过对隔离层的安全厚度数值计算,与5种理论计算方法结果进行了比较,综合求和归一法数据处理和多项式数值逼近得到了不同空区跨度与隔离层安全厚度关系。结果为露天地下开采隔离层稳定性分析的方法进行了验证和补充,也为矿山设计隔离层厚度提供了参考,对指导露天地下开采安全生产施工具有重要意义。  相似文献   

17.
A large-scale evaporate series is developed in Paleogene-Neogene strata in the Kuqa basin. The series is composed mainly of evaporate with thin beds of clastic rock (mainly mudstone and siltstone). In grayish white medium- and coarse-grained sandstone in Miocene strata, the formation of copper minerals is in close connection with brine. In joint planes, which are developed in vertical strata, are filled with gypsum. Gypsum and copper-mineralized sandstone contains enormous copper minerals, mainly atacamite. According to the SEM analysis for salt rock, gypsum rock, limestone, grayish green siltstone, grayish white medium-coarse-grained sandstone, some minerals are composed of metallic elements including Au, Ag, Cu, Zn, Pb, Co, Ni and U etc., in which Au occurs in a native form, Cu occurs in a native form or as atacamite in salt rock, gypsum rock and limestone, Ag occurs as silver sulfide in gypsum, and Zn, Pb, Co, Ni, U occur as compounds along with the above metallic ions in evaporate or clastic rock. From SEM images, we can see that metallic elements or their compounds (oxides or sulfides) “take root” as grains in salt or gypsum crystals, which belong to primary chemical sedimentation along with evaporate, while some grains “float” on surface of salt or gypsum. In the former case, mineral grains were formed together with salt (gypsum) crystals; while in the latter case, minerals were enriched from internal metallic ions (Paleogene evaporate samples) or external metallic ions (Neogene gypsum samples) in the late stage of evaporate formation. The metallic ions in Paleogene evaporate samples might originate from weathered or denudated materials in the south Tianshan Mountains. The metallic ions in the Neogene evaporate samples might be from metal-bearing brine, which migrated upward to surface along fractures and leached into evaporate (gypsum). Occurrence of metallic minerals and their compounds (elementary substance) in Paleogene evaporate proves that diversified metallic minerals exist in evaporate. The source of metallic ions in the Neogene evaporate series shows that evaporate could provide materials for late-stage metallic mineralization.  相似文献   

18.
Subsurface dissolution (subrosion) of evaporites such as halite and gypsum can lead to extensive land subsidence. Recent land subsidences have been surveyed at six separate locations in northwestern Switzerland. The diameters of the affected surface areas range from 100 to 1,500 m, and corresponding subsidence rates reached more than 100 mm/year. Based on a geometrical model, three sites could be outlined where land subsidence can likely be attributed to salt solution mining. The effects of increased hydrostatic gradient due to both groundwater withdrawal and fluid density contrasts were evaluated in more detail for the remaining sites with a series of 2D density-coupled solute-transport simulations along an approximately 1,000-m-long and 150-m-deep 2D cross section. Simulation results indicate that the upconing process of saline groundwater into the main aquifer occurs under different distributions of subsurface parameters and hydraulic boundary conditions. For the presented setup, the simulations also revealed that the most sensitive factor for the dissolution rate is the structure or dip of the halite formation, which leads to an increase of dissolution rate with increasing dip. Due to the increased density of the brine, an intrinsic flow dynamic develops which follows the direction of the dip.  相似文献   

19.
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.  相似文献   

20.
石膏对白云岩溶解影响的实验模拟研究   总被引：25，自引：2，他引：25  
黄思静  杨俊杰 《沉积学报》1996,14(1):103-109
表生到埋藏成岩作用的温度与压力（４０-１３０℃、常压－３０ＭＰａ）条件下，含膏与不含膏白云岩的溶解实验证明:在表生与相对浅埋藏的温压条件（低于７５℃，２０ＭＰａ）下，石膏（或硬石膏）的存在可不同程度地加速白云岩的溶解，随着实验温度和压力的升高，石膏（或硬石膏）对白云岩溶解的这种积极作用逐渐降低。在相对深埋藏的温压条件（高于７５℃、２０ＭＰａ）下，石膏（或硬石膏）的存在显著阻止白云岩的溶解，随着实验温度和压力的继续升高，石膏（或硬石膏）对白云岩溶解的这种消极作用也逐渐增加。从实验的这种结果可以预测，在近地表条件下和埋藏成岩作用的早期阶段，由溶解作用造成的含膏白云岩地层的次生孔隙将比不含膏的白云岩地层更为发育，因而在经历了古风化作用的地层中，含膏白云岩层更易形成良好的储层；与之相反，在相对高温高压的深埋藏成岩阶段，不含膏的白云岩地层中将更容易因酸性水的溶解作用而形成次生孔隙。因而在非蒸发沉积环境中形成的白云岩体（如正常海沉积环境的灰岩中的白云岩透镜体）更易因深埋藏溶蚀作用而形成良好的储层。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号