碳酸盐岩单裂隙渗流-溶蚀耦合模型及其参数敏感性分析

岩溶地下水系统是由碳酸盐岩裂隙含水介质演化形成的,系统初始的裂隙网络介质特征及边界条件决定了其演化过程。为揭示岩溶系统演化过程中裂隙介质特征和边界条件的影响程度,建立了裂隙溶蚀扩展的渗流-溶蚀耦合模型,并对不同边界条件下不同隙宽的单裂隙溶蚀扩展特征进行了模拟分析。结果表明:裂隙溶蚀扩展受水的侵蚀性(CO2分压)、水动力条件(水力梯度)、裂隙介质特征(裂隙初始隙宽)等综合作用影响,Ca2+的平衡浓度、水力梯度以及裂隙初始隙宽等参数的增加均能促进裂隙的快速扩展。在这些参数中,初始隙宽B0对岩溶发育的影响最为敏感,水力梯度J和Ca2+平衡浓度Ceq对岩溶发育具有相同的敏感性;此外,随着各参数值的不断增大,参数变化对岩溶发育的敏感程度越来越低。     

The motion of sediment-water mixtures during intense bedload transport: computer simulations

A new model, which couples fluid and particle dynamics, has been developed to study the motion of the sediment-water mixture during intense bedload transport, including the velocity profiles of both sediment and water, the roughness length of an upper plane bed and the thickness of moving sediment layers. Standard mixing length theory is used to model the motion of water above the boundary between the overlying water and the sediment-water mixture. The turbulent flow within the moving sediment layers is described by a shear stress model, in which the effective viscosity of the flowing water is proportional to the velocity difference between the fluid and the sediment. The particle dynamics method, in which the equations of motion of each of many particles are solved directly, is applied to model the movement of sediment particles. The particle-fluid interaction is expressed by a velocity-squared fluid drag force exerted on each sediment particle. Both computer simulation results and theoretical analysis have shown that the velocities of both sediment and fluid during intense sediment transport decrease exponentially with depth in the top layers of a fast-moving sediment—water mixture. The thickness of the moving sediment layers, obtained from the computer simulation results, is proportional to the shear stress, which agrees with previous experimental observations.     

河间隐伏型岩溶裂隙含水系统演化的数值模拟

为揭示隐伏型岩溶含水系统在不同边界条件与介质特征下的发育状况,基于地下水渗流和碳酸盐岩溶蚀理论构建耦合裂隙网络渗流与碳酸盐岩溶蚀的岩溶裂隙网络溶蚀扩展模型,模拟再现了河间隐伏型岩溶含水系统平面上发育演化过程。结果显示:岩溶系统发育过程中流量变化呈现三个阶段:第Ⅰ阶段为缓慢层流阶段,流量增速为2.2×10 -6 mL·(s·a) -1,持续时间与初始水头差呈对数关系;第Ⅱ阶段为极快速紊流阶段,流量增速为5.5×10 2 mL·(s·a) -1,持续时间与初始水头差呈线性关系;第Ⅲ阶段为快速紊流阶段,流量增速为26 mL·(s·a) -1。优势裂隙的存在使得岩溶裂隙网络溶蚀演化加快,非均质系统进入极快速紊流阶段所用时间较均质系统缩短了53%,裂隙溶蚀扩展的速度加快了12.5%。      

构造成矿非线性动力学:2.湘西金矿研究实例

应用递增应力流变学方法对湘西沃溪金矿成矿过程中构造变形、应力、流体流动、地球化学反应及成岩成矿作用等多过程间的耦合作用和非线性特征进行了模拟和分析。构造压溶作用对矿源层中成矿物质的溶解与迁移起了重要作用。断裂作用可以产生比原岩渗透率高近4个数量级的断裂渗透率,从而促进流体的流动、组分扩散和在断裂中的汇聚成矿。不同岩性对断裂作用和断裂渗透率的演化有重要影响,并且构造、流体、岩石之间存在强烈的非线性耦合作用,马底驿组中段是最有利的构造成矿层位。     

Numerical simulation of rock salt dissolution in dynamic water

The dissolution process of the rock salt is actually a complicated process in fluid dynamics and chemical kinetics; therefore, the mathematical model for the solute-fluid flow during the cavern building may be established to research the change rule of the form of the cavern boundary based on the kinetic analysis of the convection–diffusion process, with integration of the basic theories of fluid mechanics and the principles of chemical kinetics. In order to research the dissolution mechanism of the rock salt under the condition of dynamic water, the method of numerical simulation has been employed to simulate and study the factors with effect on the fluid transport mechanism in solution mining cavern of the rock salt, and its dissolution in dynamic water has been simulated and researched. The result shows that the model established in this paper may well describe the dissolution process of rock salt under the condition of dynamic water, providing basic analyzing means for the research on the rock salt dissolution in dynamic water.     

The interaction of migrating grain boundaries with fluid inclusions in rock analogues: the effect of wetting angle and fluid inclusion velocity

The distribution of fluids in grain boundaries, fluid inclusion morphology and kinetics have important effects on the evolution of microstructure and transport properties and should be understood for correct interpretations for studies of thermobarometry and paleorheology. We compare results of in situ annealing experiments on rock analogues in the presence of different pore fluids in transmitted light: bischofite with saturated brine, camphor with ethanol, and camphor with ethylene glycol. The solid–liquid systems vary in terms of wettability and solubility, while homologous temperatures, strain rates, annealing times, and the initial textures are similar. In agreement with earlier work and theory, we observe different types of grain boundary–fluid inclusion interaction at sufficiently low grain boundary velocity such as drag and drop, necking, and the break up into arrays of smaller inclusions. In all three systems the maximum possible velocity of a fluid inclusion being dragged by a moving grain boundary is dependent on the fluid inclusion size. We interpret this to reflect the fluid inclusion mobility, and compare the trend with theoretical models which suggest that for all three systems the rate-limiting process is bulk diffusion and not surface diffusion or solution-precipitation.     

Some Investigations on Early Organization of Karst System

Karstsystemsaredynamicsystemsinspontaneoussuc-cessiveevolutionarisenfromtheinteractionbetweenwaterflowandsolub1emedium-Thecircu1ationofwaterprogres-sivelyandcontinuouslymodifiesthestructureofmediumfield,consequentlymodifyingthehydraulicbehaviorofthesystemitself'Someauthors(e.g.,FordandWilliams,l989)believedthatmostcarbonaterockswillbecomeconduitaquifersifsurficialerosionisnegligibleorsufficienttimeisgiven.Karstsystemsfailtodevelopconnectedconduiteitherduetoincompleteextension(becauseinsuffici…     

The transport of heat and matter by fluids during metamorphism

Non-dimensional solutions to the equations for the combined advective and diffusive one-dimensional transport of heat and solute in a layer are derived for fixed temperature/concentration on the boundaries and initial conditions of a linear gradient across the layer or a step function at the lower boundary. The solutions allow distinction of regimes in which advective or diffusive transport of either heat or solute predominate as a function of fluid flux, time and a length scale. The much lower diffusive coefficients for solute than heat results in a significant range of length scales and fluid flux rates characterised by advection of matter and diffusion of heat. The advective velocity of a component is a function of its fluid:rock partition coefficient. The most rapidly transported tracers which partition largely into the fluid phase, such as He, will travel orders of magnitude faster than heat or compatible solutes such as oxygen. Geochemical profiles in boundary layer regions where both advective and diffusive transport are significant are shown to be particularly informative as to properties of the rocks related to fluid flow such as porosity, permeability, time scales and fluid flux rates. The importance of advection can be directly estimated from the asymmetry of the geochemical profiles across individual layers.     

钻井液流态对孔壁稳定性的影响

沉积地层中孔壁的稳定性与钻井液流态有着密切关系。结合河南叶舞凹陷盆地2 386 m盐矿钻探工程,采用不同的泵量和上返流速在不同孔段进行了钻探试验。试验和理论计算表明:在水敏性地层（黏土、泥岩、砂质泥岩等）浅层岩石孔隙率和塑性指数较高时,在钻井液层流流态情况下,主要以缩径为主,最大缩径率达26.6%;中部和深部同等地层,在紊流流态情况下,缩径和超径同时存在。其中,随着孔深的增加,缩径率呈减小趋势,超径率则有增大趋势,泥岩地层最大超径率达157.4%,盐岩地层超径率高达225.7%。      

A numerical model simulating reactive transport and evolution of fracture permeability

A numerical model is presented to describe the evolution of fracture aperture (and related permeability) mediated by the competing chemical processes of pressure solution and free‐face dissolution/precipitation; pressure (dis)solution and precipitation effect net‐reduction in aperture and free‐face dissolution effects net‐increase. These processes are incorporated to examine coupled thermo‐hydro‐mechano‐chemo responses during a flow‐through experiment, and applied to reckon the effect of forced fluid injection within rock fractures at geothermal and petroleum sites. The model accommodates advection‐dominant transport systems by employing the Lagrangian–Eulerian method. This enables changes in aperture and solute concentration within a fracture to be followed with time for arbitrary driving effective stresses, fluid and rock temperatures, and fluid flow rates. This allows a systematic evaluation of evolving linked mechanical and chemical processes. Changes in fracture aperture and solute concentration tracked within a well‐constrained flow‐through test completed on a natural fracture in novaculite (Earth Planet. Sci. Lett. 2006, in press) are compared with the distributed parameter model. These results show relatively good agreement, excepting an enigmatic abrupt reduction in fracture aperture in the early experimental period, suggesting that other mechanisms such as mechanical creep and clogging induced by unanticipated local precipitation need to be quantified and incorporated. The model is applied to examine the evolution in fracture permeability for different inlet conditions, including localized (rather than distributed) injection. Predictions show the evolution of preferential flow paths driven by dissolution, and also define the sense of permeability evolution at field scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Two-dimensional reactive transport modeling of the alteration of a fractured limestone by hyperalkaline solutions at Maqarin (Jordan)

Two-dimensional reactive transport modeling of the Maqarin Eastern Springs site, a natural analogue for the alteration of a fractured limestone by high-pH Portland cement waters, has been performed using the CrunchFlow code. These 2D calculations included transport by advection–dispersion–diffusion along a single fracture and diffusion in the wall rock. Solute transport was coupled to mineral dissolution and precipitation. A limited sensitivity analysis evaluated the effect of different values of primary mineral surface areas, flow velocity and sulfate concentration of the inflowing high-pH solution.Major secondary minerals include ettringite–thaumasite, C–S–H/C–A–S–H and calcite. C–S–H/C–A–S–H precipitation is controlled by the dissolution of primary silicates. Ettringite precipitation is controlled by diffusion of sulfate and aluminum from the wall rock to the fracture, with aluminum provided by the dissolution of albite. Calcite precipitation is controlled by diffusion of carbonate from the wall rock. Extents of porosity sealing along the fracture and in the fracture-wall rock interface depend on assumptions regarding flow velocity and composition of the high-pH solution. The multiple episodes of fracture sealing and reactivation evidenced in the fracture infills were not included in the simulations. Results can qualitatively reproduce the reported decrease in porosity in the fractures and in the wall rock next to the fractures. Instances of porosity increase next to fractures caused by carbonate dissolution were not reproduced by the calculations.     

Grain boundaries as microreactors during reactive fluid flow: experimental dolomitization of a calcite marble

Limestone dolomitization is an example of a fluid-induced mineralogical transformation that commonly affects extensive rock volumes. To understand the mechanisms enabling these efficient replacement reactions, we investigated experimentally the dolomitization of a fractured calcite marble under flow-through conditions at mild hydrothermal conditions. Contrary to most earlier studies of coupled dissolution reprecipitation reactions that were conducted using small, individual grains, in this study, the integrity of the rock was preserved, so that the experiment explored the links between flow in a fracture and fluid–rock interaction. In these experiments, grain boundaries acted as microreactors, in which a Mg-poor ‘protodolomite’ formed initially, and then transformed into dolomite. The difficulty in nucleating dolomite played a key role in controlling the evolution of the porosity, by allowing for (1) initial dissolution along grain boundaries, and (2) formation of coarse porosity at the reaction interface. This porosity evolution not only enabled the reaction to progress efficiently, but also controlled the mineralogy of the system, as shown by brucite replacing calcite near the fracture once the fluid along calcite grain boundaries became sufficiently connected to the fluid flowing through the fracture. This study illustrates the role of grain boundaries, porosity evolution and nucleation in controlling reaction progress as well as the nature and textures of the products in pervasive mineralogical transformations.     

初始裂隙对岩溶水紊流形成的影响

岩溶地区地下发育着大量的溶洞和地下河管道,地下水流状态既有层流也有紊流,而紊流是溶洞管道形成的重要条件。紊流的形成受到岩石初始裂隙的影响,初始裂隙的张开度、分布、走向、迹长、密度等因素都影响着裂隙发育过程中水流状态的变化。通过对不同统计特征的初始裂隙网络进行水流和溶蚀的数值模拟发现,以张开度标准差反映的裂隙网络非均匀性越强,模拟紊流出现的时间就越早;主要裂隙的存在使裂隙网络的非均性增强,主要裂隙与水力梯度总方向的角度越小,紊流出现的时间就越早;当裂隙平均迹长过小时会导致裂隙连通性较差,影响裂隙水流和溶蚀作用;裂隙密度,尤其是主要裂隙密度,对岩溶发育的影响较大。相对于次要裂隙,如果主要裂隙密度偏小,紊流形成时间会大大增加,甚至很难形成紊流。当初始裂隙张开度小于0.001 cm,增大水力梯度仍没有紊流发生,岩溶几乎不发育。      

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

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

A theoretical and numerical framework for modeling gypsum cavity dissolution

In this paper, a large‐scale diffuse interface model is used to describe the evolution of a gypsum cavity formation induced by dissolution. The method is based upon the assumption of a pseudo‐component dissolving with a thermodynamic equilibrium boundary condition. A methodology is proposed based on numerical computations with fixed boundaries in order to choose suitable parameters for the diffuse interface model, and hence predict the correct dissolution fluxes and surface recession velocity. Additional simulations were performed to check which type of momentum balance equation should be used. The numerical results did not show a strong impact of this choice for the typical initial boundary value problems under consideration. Calculations with a variable density and Boussinesq approximation were also performed to evaluate the potential for natural convection. The results showed that the impact of density driven flows was negligible in the cases under investigation. The potential of the methodology is illustrated on two large‐scale configurations: one corresponding to a gypsum lens located strictly within a porous rock formation and the other to an isolated pillar in a flooded gypsum room and pillar quarry. Copyright © 2016 John Wiley & Sons, Ltd.     

Prediction of rock alteration patterns: A potential tool in mineral exploration

Hydrothermal alteration of a quartz‐K‐feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO₂ gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300°C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass‐action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass‐balance equations are solved sequentially using an implicit scheme in a finite‐element code. The pore‐fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K‐feldspar. Our model simulates, in a simplified way, the acid‐induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin.     

层流—紊流共存流场中岩溶裂隙网络演化过程的数值模拟方法

研究岩溶水系统的演化过程对许多资源与环境问题有着重要意义。岩溶地区水资源、油气资源的预测与开发、水土流失成因及防治等课题都与此密切相关。岩溶管道是地下水对裂隙的逐渐溶蚀扩展形成的,岩溶演化初期,所有裂隙宽度不大,流场整体成层流状态。随着溶蚀的进行,一部分裂隙优先扩大使其中的水流进入紊流状态,而另一部分裂隙中水流仍呈层流状态。文章提出一种数值方法,能够模拟层流-紊流共存流场的岩溶演化过程。利用蒙特卡洛(Monte-Carlo)方法模拟初始裂隙网络,通过非连续介质方法模拟裂隙网络中的渗流。裂隙扩展的速度通过岩石表面溶蚀速度经验公式计算,使用迭代方法求解层流-紊流共存条件下的水头非线性方程。构建了能够利用解析法求解的模型,把数值解和解析解的结果进行对比,验证了本研究的数值法及软件的可行性。      

构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟

建立了一个综合的构造流体成矿体系的反应输运力学耦合动力学模型。利用有限元方法求解岩石变形、断裂作用和断裂网络统计动力学、流体流动、有机和无机地球化学反应及成岩成矿作用、压力溶液和其它压实力学、热迁移的方程组 ,可以对构造流体成矿体系的动力学演化过程进行 1～ 3维数值模拟。模拟的主要内容是在各种过程耦合作用下描述构造流体成矿体系的主要变量的时空演化 :( 1)与成矿流体的形成和性质有关的变量 ,如地层中矿物 (包括成矿物质 )的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等 ;( 2 )与构造变形和流体运移有关的各变量 ,如应力与变形速率、岩石孔隙度、构造 (断裂 )渗透率等 ;( 3 )与沉淀成矿有关的变量 ,如矿物 (金属矿物和脉石矿物 )的成核速率、各矿物的沉淀量等 ;( 4 )上述各有关变量间的时空耦合关系 ,如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。以湖南沃溪金锑钨矿床为例 ,应用该模型和方法对成矿动力学过程和动力学机制进行了初步的模拟与分析。     

单轴压缩条件下岩溶化裂隙岩体损伤破坏特征研究

岩溶化裂隙岩体是普遍发育于自然界中具有初始损伤的岩体。为了研究岩溶化裂隙岩体损伤破坏特征,本文以贵州某地赋存的溶蚀岩体为研究对象,运用损伤力学理论构建岩溶化裂隙岩体在单轴压缩条件下的损伤演化模型,并建立岩溶化裂隙岩体损伤演化方程。采用颗粒流数值软件进行单轴压缩数值试验,进一步研究岩溶化裂隙岩体试件在单轴压缩条件下的损伤演化特征,分析岩溶化裂隙岩体的微观损伤特征。结果表明:岩溶化裂隙岩体的初始损伤主要包括溶蚀损伤和裂隙损伤。岩溶化裂隙岩体的初始损伤随着溶蚀率的增加而增加,最终增加速率趋于平缓;岩溶化裂隙岩体的损伤演化曲线均呈“S”型变化,先缓慢增加,再迅速增加,最后缓慢增加至损伤值1;岩溶化裂隙岩体存在异构特征,导致破坏裂隙起源于具有初始损伤的溶蚀孔洞和裂隙处,随后裂隙经历萌发、扩张和剪切作用、数量和长度增加以及裂隙贯通4个阶段后发生宏观破坏。     

横流环境湍射流涡动力学特性数值模拟

采用三维湍流模型及其混合有限分析解法进行横流中湍射流这一复杂三维流动的研究。在利用实验数据对模型及其计算方法验证的基础上,对多种喷口形式和流速比工况下的流速场和涡量场进行了数值计算,模拟得到了其旋涡结构发展演化特性。在射流初始阶段,横流在射流背流面形成绕流分离旋涡,其结构与射流喷口形式和流速比有关。射流主轴沿流向布置的窄缝射流产生的旋涡最为特殊和复杂,在喷口侧面存在4个分离点,从而形成4个旋涡。在射流远区,流动主要由反向旋转涡对所控制,并诱导出二次涡对。分析得出了射流喷口形式及流速比对反向旋转涡对涡核位置和旋转强度的影响。