共查询到20条相似文献，搜索用时 62 毫秒
1.
A numerical model is developed for describing the transport of virus in a fracture-matrix coupled system with fracture-skin.An advective dispersive virus transport equation,including firstorder sorption and inactivation constant is used for simulating the movement of viruses.Implicit finite-difference numerical technique is used to solve the coupled non-linear governing equations for the triple continuum model consisting of fracture,fracture-skin and the rock-matrix.A varying grid is adopted at the fracture and fracture-skin interface to capture the mass transfer.Sensitivity analysis was performed to investigate the effect of various properties of the fracture-skin as well as viruses on the virus concentration in the fractured formation with fracture-skin.Simulation results suggest that the virus concentration in the fracture decreases with increment in the fracture-skin porosity,fracture-skin diffusion coefficient,mass transfer coefficient,inactivation constant and sorption distribution coefficient, and with reduction in the fracture aperture. 相似文献
2.
In this study, the behavior of thermal fronts along the fracture is studied in the presence of fracture-skin in a coupled
fracture-matrix system. Cold water is injected into the fracture, which advances gradually towards production well, while
extracting heat from the surrounding reservoir matrix. The heat conduction into the fracture-skin and the rock-matrix from
the high permeability fracture is assumed to be one dimensional perpendicular to the axis of fluid flow along the fracture.
Constant temperature cold water is injected through an injection well at the fracture inlet. The fluid flow takes place along
the horizontal fracture which ensures connectivity between the injection and production wells. Since the rock-matrix is assumed
to be tight, the permeability of fracture-skin as well as the rock-matrix is neglected. The present study focuses on the heat
flux transfer at the fracture-skin interface as against the earlier studies on fracture-matrix interface, and the sensitivity
of additional heterogeneity in the form of fracture skin in a conventional fracture-matrix coupled system is studied. The
behavior of thermal fronts for various thermal conductivity values of the fracture-skin and rock-matrix is analyzed. Spatial
moment analysis is performed on the thermal distribution profiles resulting from numerical studies in order to investigate
the impact on mobility and dispersion behavior of the fluid in the presence of fracture-skin. The presence of fracture skin
affects the heat transfer significantly in the coupled fracture-matrix system. The lower order spatial moments indicate that
the effective thermal velocity increases with increase in skin thermal conductivity and a significant thermal dispersion is
observed at the inlet of the fracture owing to the high thermal conductivity of the fracture-skin at the early stages. Furthermore
the higher spatial moments indicate that the asymmetricity increases with decrease in skin thermal conductivity unlike the
case with half fracture aperture and fluid velocity and the kurtosis is maximum with higher skin thermal conductivity which
implies enhanced heat extraction from the fracture-skin into the fracture. Results suggest that the amount of heat extraction
by the circulating fluid within the fracture from the reservoir not only depends on the rock-matrix module of the reservoir
characteristics but also the fracture-skin characteristics of the system and subsequently influence the reservoir efficiency. 相似文献
3.
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. 相似文献
4.
This study deals with transport of solutes through a saturated sub-surface rock formation with well-defined horizontal parallel fractures. For this purpose, a simplified conceptual model consisting of a single fracture and its associated rock-matrix is considered in the presence of a fracture-skin in order to study the mobility and mixing of solutes along the fracture. In this paper, a coupled fracture-skin-matrix system is modeled numerically using finite difference method in a pseudo two-dimensional domain with a constant continuous source at fracture inlet. Flow and transport processes are considered parallel to the fracture axis, while the transport processes in fracture-skin as well as in rock-matrix are considered perpendicular to the fracture axis. Having obtained the concentration distribution along the fracture, method of spatial moments is employed to study the mobility and spreading of solutes. Sensitivity analyses have been done to understand the effect of various fracture-skin parameters like porosity, thickness, and diffusion coefficient. Further, the influence of non-linear sorption and radioactive decaying of solutes are carried out for different sorption intensities and decay constants. Results suggest that the presence of fracture-skin significantly influences the mobility and spreading of solutes along the fracture in comparison with a coupled fracture-matrix system without fracture-skin. 相似文献
5.
Hide Yasuhara Naoki Kinoshita Hiroaki Ohfuji Dae Sung Lee Shinichiro Nakashima Kiyoshi Kishida 《Applied Geochemistry》2011,26(12):2074-2088
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. 相似文献
6.
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. 相似文献
7.
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. 相似文献
8.
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. 相似文献
9.
构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟 总被引：12，自引：0，他引：12
建立了一个综合的构造流体成矿体系的反应输运力学耦合动力学模型。利用有限元方法求解岩石变形、断裂作用和断裂网络统计动力学、流体流动、有机和无机地球化学反应及成岩成矿作用、压力溶液和其它压实力学、热迁移的方程组 ,可以对构造流体成矿体系的动力学演化过程进行 1～ 3维数值模拟。模拟的主要内容是在各种过程耦合作用下描述构造流体成矿体系的主要变量的时空演化 :( 1)与成矿流体的形成和性质有关的变量 ,如地层中矿物 (包括成矿物质 )的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等 ;( 2 )与构造变形和流体运移有关的各变量 ,如应力与变形速率、岩石孔隙度、构造 (断裂 )渗透率等 ;( 3 )与沉淀成矿有关的变量 ,如矿物 (金属矿物和脉石矿物 )的成核速率、各矿物的沉淀量等 ;( 4 )上述各有关变量间的时空耦合关系 ,如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。以湖南沃溪金锑钨矿床为例 ,应用该模型和方法对成矿动力学过程和动力学机制进行了初步的模拟与分析。 相似文献
10.
The effects of stress/deformation on fluid flow and contaminant transport in fractured rocks is one of the major concerns for performance and safety assessments of many subsurface engineering problems, especially radioactive waste disposal and oil/gas reservoir fields. However, very little progress has been made to study this issue due to difficulties in both experiments and numerical modeling. The objective of this study is to systematically investigate the influence of stress on solute transport in fractured rocks for the first time, considering different stress and hydraulic pressure conditions. A hybrid approach combining discrete element method (DEM) for stress-flow simulations and a particle tracking algorithm is developed. The impact of matrix diffusion (diffusion of molecular size solutes in and out of the rock matrix, and sorption onto the surface of micropores in rock matrix) is also included. The numerical results show that stress not only significantly changes the solute residence time through the fracture networks, but also changes the solute travel paths. Matrix diffusion plays a dominant role in solute transport when the hydraulic gradient is small, which is often encountered in practice. 相似文献
11.
A numerical model using a hybrid formulation of a finite element method (FEM) coupled with the volume of fluid (VOF) technique to simulate the fracture grouting processes in soils is described. The numerical model considered the couplings of the stress distribution, with two-phase fluid flows, and the mesh element damage. The hardening of grout in soil is described by a time-dependent Young’s modulus and viscosity. Crack initiation, branching, propagation, and grout vein growth in homogeneous and heterogeneous soils can be numerically reproduced. Although the method is developed particularly for simulating fracture grouting, the processes of compaction grouting and permeation grouting can also be numerically simulated. Some grouting cases have been simulated with results similar to the experimental results. This further confirms the adequacy and the power of the numerical approach. 相似文献
12.
Effects of linking up of discontinuities on fracture growth and groundwater transport 总被引：3，自引：1，他引：2
It is proposed that the growth of fractures is the basic process for generating and maintaining permeability in solid rock
(bedrock). Many extension fractures grow as hydrofractures, whereas many shear (and extension) fractures grow through the
formation of transverse fractures that connect the adjacent tips of existing fractures. In a boundary-element analysis, the
hydrofractures are modeled as being driven open by a fluid overpressure that varies linearly from 10 MPa at the fracture centre
to 0 MPa at the fracture tip. The host rock has a uniform Young's modulus of 10 GPa, a Poisson's ratio of 0.25, and is dissected
by vertical joints and horizontal contacts, each of which is modeled as an internal spring of stiffness 6 MPa m−1. The number of joints and contacts, and their location with respect to the hydrofracture tip are varied in different model
runs. The results of the analyses indicate that the tensile stresses generated by overpressured hydrofractures open up joints
and contacts out to considerable distances from the fracture tip, so that they tend to link up to form a hydraulic pathway.
Using the same Young's modulus, Poisson's ratio, and internal spring constant for joints as in the hydrofracture models, boundary-element
models were made to study the interaction stresses that cause neighbouring joints to become interconnected through the growth
of linking transverse fractures that, ultimately, may evolve into shear fractures. The models were subjected to tensile stress
of 6 MPa acting normal to the joint planes as the only loading. The offset (horizontal distance) and underlap (vertical distance)
between the adjacent tips of the joints were varied between model runs. The results show a concentration of tensile and shear
stresses in the regions between the neighbouring tips of the joints, but these regions become smaller when the underlap of
the joints decreases and changes to overlap. These stress-concentration regions favour the development of transverse (mostly
shear) fractures that link up the nearby tips of the joints, so as to form a segmented shear or extension fracture. Analytical
results on aperture variation of a hydrofracture in a homogeneous, isotropic rock are compared with boundary-element results
for a hydrofracture dissecting layered rocks. The aperture is larger where the hydrofracture dissects soft (low Young's modulus)
layers than where it dissects stiff layers. Aperture variation may encourage subsequent groundwater-flow channeling along
a pathway generated by a hydrofracture in layered rocks.
Electronic Publication 相似文献
13.
渗流-化学溶解耦合作用下岩石单裂隙渗透特性研究 总被引：3，自引：0，他引：3
为揭示在渗流-化学溶解耦合作用下单裂隙渗透特性的变化规律,建立了描述二维渗流-化学溶解耦合作用的偏微分方程组,并利用COMSOL Multiphysics软件成功地求解该方程组。首先,模拟了文献[1]中的盐岩渗流-溶解耦合渗流试验结果,数值模拟结果与试验结果较为吻合,验证了数学模型的正确性和有效性。然后,利用分形理论生成了一个粗糙的裂隙面数字模型,着重分析了二维石灰岩粗糙裂隙面在水流、矿物溶解和输运过程中其渗透特性的变化规律。数值分析显示,（1）溶质浓度对裂隙面的溶解具有非常重要的作用,从而水流进口端的溶解厚度比出口端大得多。（2）裂隙的整体渗透性在初始时刻增加较慢,随着裂隙开度的增大和贯通,溶解速度会逐渐增大,是一个加速的过程。 相似文献
14.
This paper presents numerical investigation on the ore-forming fluid migration driven by
tectonic deformation and thermally-induced buoyancy force in the Chanziping ore district in South
China. A series of numerical scenarios are considered to examine the effect of meteoric water precipitation,
the dip angle of the faults, unconformity surface, and thermal input on the ore genesis. Our computations
reveal that the downward basinal fluid flow driven by extensional stress mixes with the upward
basal fluid driven by the thermal input from depth at the junction of two faults at a temperature of about
200 C, triggering the precipitation of the Chanziping uranium deposit. 相似文献
15.
构造成矿非线性动力学:1.递增应力流变学模型 总被引：1，自引：1，他引：1
本文提出了构造成矿作用动力学研究的递增应力流变学方法。利用该方法将构造变形、应力、流体流动、地球化学反应及成岩成矿作用等多过程耦合起来，可以从多种地质过程的耦合与反馈作用对构造成矿的动力学演化过程进行１～３维数值模拟。模拟的主要内容是在各种过程耦合作用下，以下描述构造成矿体系的主要变量的时空演化：（１）与成矿流体的形成和性质有关的变量，如地层中矿物（包括成矿物质）的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等；（２）与构造变形和流体运移有关的各变量，如应力与变形速率、岩石孔隙度、构造（断裂）渗透率等；（３）与沉淀成矿有关的变量，如矿物（金瞩矿物和脉石矿物）的成核速率、各矿物的沉淀量等；（４）上述各有关变量间的时空耦合关系，如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。 相似文献
16.
碳酸盐岩单裂隙渗流-溶蚀耦合模型及其参数敏感性分析 总被引：1，自引：1，他引：0
岩溶地下水系统是由碳酸盐岩裂隙含水介质演化形成的,系统初始的裂隙网络介质特征及边界条件决定了其演化过程。为揭示岩溶系统演化过程中裂隙介质特征和边界条件的影响程度,建立了裂隙溶蚀扩展的渗流-溶蚀耦合模型,并对不同边界条件下不同隙宽的单裂隙溶蚀扩展特征进行了模拟分析。结果表明:裂隙溶蚀扩展受水的侵蚀性(CO2分压)、水动力条件(水力梯度)、裂隙介质特征(裂隙初始隙宽)等综合作用影响,Ca2+的平衡浓度、水力梯度以及裂隙初始隙宽等参数的增加均能促进裂隙的快速扩展。在这些参数中,初始隙宽B0对岩溶发育的影响最为敏感,水力梯度J和Ca2+平衡浓度Ceq对岩溶发育具有相同的敏感性;此外,随着各参数值的不断增大,参数变化对岩溶发育的敏感程度越来越低。 相似文献
17.
A coupled hydro-chemo-mechanical numerical model is developed for these coupled phenomena in many engineering fields. The
model has been applied to predicting the response of a stressed rockmass column to an injected reactive fluid (reagent) flow.
The response includes evolutions of porosity, permeability, reagent and mineral concentrations during dissolution. In the
model, the progress of dissolution is defined by the change in porosity ratio and the porosity increases with dissolution
assuming there is no precipitation. The numerical evolutions of porosity, permeability, reagent and mineral concentrations
during dissolution are validated against steady state solutions. The model results show that these evolutions are regulated
to a certain extent by the applied external loadings: an applied extensional stress enhances the progress of the dissolution
process while an applied compression stress slows the progress of the dissolution process.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
18.
Permeability and stress in crystalline rocks 总被引：2，自引：0，他引：2
Groundwater from crystalline rocks is a significant resource in many areas of the world. It is also an important medium for contaminant transport from, for example, deep nuclear waste repositories. Stress distributions in fractured rocks are important in controlling groundwater flow in several ways: (i) palaeostress fields are responsible for the evolution of fracture systems which transmit groundwater; (ii) current in situ stress fields will influence the shape and aperture of fractures; (iii) humans can influence the natural stress field in a rock mass to enhance fracture flows. The significance of stresses for groundwater flow can be investigated by field techniques (hydraulic fracturing), laboratory techniques (stress cells) or by numerical modelling. 相似文献
19.
This paper proposes a three-dimensional coupled hydrothermal model for fractured rock based on the finite-discrete element method to simulate fluid flow and heat transport. The 3D coupled hydrothermal model is composed of three main parts: a heat conduction model for the rock matrix, a heat transfer model for the fluid in the fractures (including heat conduction and heat convection), and a heat exchange model between the rock matrix and the fluid in the fractures. Four examples with analytical solutions are provided to verify the model. A heat exchange experiment of circulating water in a cylindrical granite sample with one fracture is simulated. The simulation results agree well with the experimental results. The effects of the fracture aperture, fluid viscosity, and pressure difference on the heat exchange between the fluid and rock are studied. Finally, an application concerned with heat transport and fluid flow in fractured rock is presented. The simulation results indicate that the 3D fully coupled hydrothermal model can capture the fluid flow and temperature evolution of rocks and fluids. 相似文献
20.
W. Zhou 《Environmental Geology》2001,40(7):797-808
Two-phase flow in fractured rock is an important phenomenon related to a range of practical problems, including non-aqueous phase liquid contamination of groundwater. Although fractured rocks consist of fracture networks, the study of two-phase flow in a single fracture is a pre-requisite. This paper presents a conceptual and numerical model of two-phase flow in a variable fracture. The void space of the fracture is conceptualized as a system of independent channels with position-dependent apertures. Fundamental equations, governing two-phase displacement in each channel, are derived to represent the interface positions and fractional flows in the fracture. For lognormal aperture distributions, simple approximations to fractional flows are obtained in analytical form by assuming void occupancy based on a local capillary allowability criterion. The model is verified by analytical solutions including two-phase flow in a parallel-plate fracture, and used to study the impacts of aperture variation, mobility ratio and fracture orientation on properties of two-phase flow. Illustrative examples indicate that aperture variation may control the distribution of wetting and non-wetting fluids within the fracture plane and hence the ability of the fracture to transmit these fluids. The presence of wetting fluid does little to hinder non-wetting fluid flow in fractures with large aperture variations, whereas a small volume of non-wetting fluid present in the fracture can significantly reduce wetting fluid flow. Large mobility ratios and high fracture slope angles facilitates the migration of non-wetting fluid through fractures. 相似文献