共查询到20条相似文献，搜索用时 281 毫秒
1.
用格子波耳兹曼方法模拟双重孔隙介质中的流体迁移 总被引：3，自引：0，他引：3
作者在本文中介绍了基于格子波耳兹曼模型的双重孔隙介质中流体运移的数值模拟计算方法。我们从格子波耳兹曼碰撞模型出发，利用格子波耳兹曼方程、Chapman-Enskog展开，以及多尺度技术，得到了描述双重孔隙介质中流体迁移的二维扩散方程。利用格子气自动机方法计算该扩散方程，实现了对双重孔隙介质中流体运移过程的数值模拟仿真。数值实验表明，我们所使用的方法正确、有效。 相似文献
2.
Donald J. DePaolo 《Geochimica et cosmochimica acta》2006,70(5):1077-1096
A mathematical model is presented that describes the effects of pore fluid aqueous diffusion and reaction rate on the isotopic exchange between fluids and rocks in reactive geo-hydrological systems where flow is primarily through fractures. The model describes a simple system with parallel equidistant fractures, and chemical transport in the matrix slabs between fractures by aqueous diffusion through a stagnant pore fluid. The solid matrix exchanges isotopes with pore fluid by solution-precipitation at a rate characterized by a time constant, R (yr−1), which is an adjustable parameter. The effects of reaction on the isotopes of a particular element in the fracture fluid are shown to depend on the ratio of the diffusive reaction length for that element (L) to the fracture spacing (b). The reaction length depends on the solid-fluid exchange rate within the matrix, the partitioning of the element between the matrix pore fluid and the matrix solid phase, the porosity and density of the matrix, and the aqueous diffusivity. For L/b < 0.3, fluid-rock isotopic exchange is effectively reduced by a factor of 2L/b relative to a standard porous flow (single porosity) model. For L/b > 1, the parallel fracture model is no different from a porous flow model. If isotopic data are available for two or more elements with different L values, it may be possible to use the model with appropriate isotopic measurements to estimate the spacing of the primary fluid-carrying fractures in natural fluid-rock systems. Examples are given using Sr and O isotopic data from mid-ocean ridge (MOR) hydrothermal vent fluids and Sr isotopes in groundwater aquifers hosted by fractured basalt. The available data for MOR systems are consistent with average fracture spacing of 1-4 m. The groundwater data suggest larger effective fracture spacing, in the range 50-500 m. In general, for fractured rock systems, the effects of fracture-matrix diffusive exchange must be considered when comparing isotopic exchange effects for different elements, as well as for estimating water age using radioactive and cosmogenic isotopes. 相似文献
3.
Thermal recovery from a hot dry rock (HDR) reservoir viewed as a deformable fractured medium is investigated with a focus on the assumption of local thermal non‐equilibrium (LTNE). Hydraulic diffusion, thermal diffusion, forced convection and deformation are considered in a two‐phase framework, the solid phase being made by impermeable solid blocks separated by saturated fractures. The finite element approximation of the constitutive and field equations is formulated and applied to obtain the response of a generic HDR reservoir to circulation tests. A change of time profile of the outlet fluid temperature is observed as the fracture spacing increases, switching from a single‐step pattern to a double‐step pattern, a feature which is viewed as characteristic of established LTNE. A dimensionless number is proposed to delineate between local thermal equilibrium (LTE) and non‐equilibrium. This number embodies local physical properties of the mixture, elements of the geometry of the reservoir and the production flow rate. All the above properties being fixed, the resulting fracture spacing threshold between LTNE and LTE is found to decrease with increasing porosity or fluid velocity. The thermally induced effective stress is tensile near the injection well, illustrating the thermal contraction of the rock, while the pressure contribution of the fracture fluid is negligible during the late period. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
4.
引言岩石和岩体的本构关系(包括适用于可压缩流体的状态方程)是岩石介质的应力(或压力)、应变(或比容)、比内能(或温度)、及加载路径、应变率等之间的普遍关系式.它是研究应力波在岩石介质中的传播规律及其与地下结构相互作用的基础;在估价地下结构对各种爆炸作用的抵抗能力以设计经得起这些作用的防护结构时,在预报周围环境地运动及确定结构--介质相互作用时,都需要岩石和岩体介质的动力响应的精确描述; 相似文献
5.
Wang Wu-ling Feng Yi-xing 《岩土力学》1980,2(1):15-36
The constitutive relations of rock and rock mass are essential to the study of propagation of stress wave in rock medium and its interaction with structure .in this paper ,the development of the research on constitutive relations of rock and rock mass is briefly reviewed ;the equations of state of rock in fluid state at high pressures and the con equations of rock in solid state are mainly described .The cop model as the constitutive equation of rock and its generaligations especially discussed . 相似文献
6.
Double Porosity Finite Element Method for Borehole Modeling 总被引：2，自引：0，他引：2
Summary. This paper considers the mechanical and hydraulic response around an arbitrary oriented borehole drilled in a naturally fractured formation. The formation is treated as a double porosity medium consisting of the primary rock matrix as well as the fractured systems, which are each distinctly different in porosity and permeability. The poro-mechanical formulations that couple matrix and fracture deformations as well as fluid flow aspects are presented. A double porosity and double permeability finite element solution for any directional borehole drilled in the fractured porous medium is given. Compared with conventional single-porosity analyses, the proposed double-porosity solution has a larger pore pressure in the matrix and a smaller tensile stress in the near-wellbore region. The effects of time, fracture, mud weight, and borehole inclination in the double-porosity solution are parametrically studied to develop a better understanding of physical characteristics governing borehole problems. 相似文献
7.
水泥基速凝浆液广泛应用于富水裂隙围岩注浆工程。为研究水泥基速凝浆液在裂隙注浆中的扩散机制与压力特性,选取水泥-水玻璃浆液及高聚物改性水泥浆液,采用不同水灰比及混合体积,开展了静水条件下单一平板裂隙的注浆扩散模型试验,通过监测浆液的注浆扩散形态和注浆过程中不同测点的压力变化,研究了水泥基速凝浆液的静水扩散规律和注浆扩散压力分布规律。试验结果表明,水泥-水玻璃浆液和高聚物改性水泥浆液具有相似的扩散形态和压力分布变化,两者的显著差别在于高聚物改性水泥浆液扩散所需的注浆压力明显小于水泥-水玻璃浆液,且凝结固化时间相对较长,有利于延长注浆扩散距离。研究结果应用于张马屯铁矿注浆堵水工程,解决了裂隙岩体渗漏水问题,实现了围岩稳定性控制。研究成果具有一定的理论和实践意义,对类似工程具有参考价值。 相似文献
8.
This paper presents the development of a mathematical model for methane gas migration in coal seams. The major focus of this model is the coupling between the gas flow and deformation of solid coal. The effect of diffusion of adsorbed methane gas from the solid matrix to the voids has been taken into account. The adsorption of gas in the coal seam causes a two-phase state of gas flow. The governing equation for the two-phase gas flow is a non-linear partial differential equation with non-linear boundary conditions. A finite element model has been developed for simulation of the distribution of pressure and concentration of methane gas due to gas migration in coal seams. 相似文献
9.
Numerical Simulation of 3D Hydraulic Fracturing Based on an Improved Flow-Stress-Damage Model and a Parallel FEM Technique 总被引：4，自引：4，他引：0
L. C. Li C. A. Tang G. Li S. Y. Wang Z. Z. Liang Y. B. Zhang 《Rock Mechanics and Rock Engineering》2012,45(5):801-818
The failure mechanism of hydraulic fractures in heterogeneous geological materials is an important topic in mining and petroleum engineering. A three-dimensional (3D) finite element model that considers the coupled effects of seepage, damage, and the stress field is introduced. This model is based on a previously developed two-dimensional (2D) version of the model (RFPA2D-Rock Failure Process Analysis). The RFPA3D-Parallel model is developed using a parallel finite element method with a message-passing interface library. The constitutive law of this model considers strength and stiffness degradation, stress-dependent permeability for the pre-peak stage, and deformation-dependent permeability for the post-peak stage. Using this model, 3D modelling of progressive failure and associated fluid flow in rock are conducted and used to investigate the hydro-mechanical response of rock samples at laboratory scale. The responses investigated are the axial stress–axial strain together with permeability evolution and fracture patterns at various stages of loading. Then, the hydraulic fracturing process inside a rock specimen is numerically simulated. Three coupled processes are considered: (1) mechanical deformation of the solid medium induced by the fluid pressure acting on the fracture surfaces and the rock skeleton, (2) fluid flow within the fracture, and (3) propagation of the fracture. The numerically simulated results show that the fractures from a vertical wellbore propagate in the maximum principal stress direction without branching, turning, and twisting in the case of a large difference in the magnitude of the far-field stresses. Otherwise, the fracture initiates in a non-preferred direction and plane then turns and twists during propagation to become aligned with the preferred direction and plane. This pattern of fracturing is common when the rock formation contains multiple layers with different material properties. In addition, local heterogeneity of the rock matrix and macro-scale stress fluctuations due to the variability of material properties can cause the branching, turning, and twisting of fractures. 相似文献
10.
John M. Eiler Lukas P. Baumgartner John W. Valley 《Contributions to Mineralogy and Petrology》1992,112(4):543-557
We formulated a numerical model for stable isotope interdiffusion which predicts the temperatures recorded between two or more minerals, and the intragranular distribution of stable isotopes in each mineral, as functions of mineral grain sizes and shapes, diffusivities, modes, equilibrium isotopic fractionations, and the cooling rate of a rock. One of the principal assumptions of the model is that grain boundaries are regions of rapid transport of stable isotopes. This Fast Grain Boundary (FGB) model describes interdiffusion between any number of mineral grains, assuming that local equilibrium and mass balance restrictions apply on the grain boundaries throughout the volume modeled. The model can be used for a rock containing any number of minerals, and number of grain sizes of each mineral, several grain shapes, and any thermal history or domain size desired. Previous models describing stable isotope interdiffusion upon cooling have been based on Dodson's equation or an equivalent numerical analogue. The closure temperature of Dodson is the average, bulk temperature recorded between a mineral and an infinite reservoir. By using Dodson's equation, these models have treated the closure temperature as an innate characteristic of a given mineral, independent of the amounts and diffusion rates of other minerals. Such models do not accurately describe the mass balance of many stable isotope interdiffusion problems. Existing models for cation interdiffusion could be applied to stable isotopes with some modifications, but only describe exchange between two minerals under specific conditions. The results of FGB calculations differ considerably from the predictions of Dodson's equation in many rock types of interest. Actual calculations using the FGB model indicate that closure temperature and diffusion profiles are as strongly functions of modal abundance and relative differences in diffusion coefficient as they are functions of grain size and cooling rate. Closure temperatures recorded between two minerals which exchanged stable isotopes by diffusion are a function of modal abundance and differences in diffusion coefficient, and may differ from that predicted by Dodson's equation by hundreds of degrees C. Either or both of two minerals may preserve detectable zonation, which may in some instances be larger in the faster diffusing mineral. Rocks containing three or more minerals can record a large span of fractionations resulting from closed system processes alone. The results of FGB diffusion modeling indicate that the effects of diffusive exchange must be evaluated before interpreting mineral fractionations, concordant or discordant, recorded within any rock in which diffusion could have acted over observable scales. The predictions of this model are applicable to thermometry, evaluation of open or closed system retrogression, and determination of cooling rates or diffusion coefficients. 相似文献
11.
The development of outgrowths or pits of various shapes on coated grains is explained via a quantitative model of grain growth/dissolution kinetics coupled to evolving grain geometry (morphological dynamics). Grain-coating thinning or fracturing occurring due to nonplanar growth (and consequent grain surface area increase) is shown to underlie an instability to the formation of bumps, or in the case of undersaturated systems, to pitting. Examples of diagenetic outgrowth phenomena on clay-coated quartz are presented. A quantitative model of coupled quartz growth and coating dynamics is shown to imply many features observed in natural systems. Crystal growth anisotropy is shown to strongly influence the morphology of the outgrowths. The creation of inclusions is shown to be closely related to the present morphological instability. These morphological instability phenomena are interesting examples of geochemical self-organization.A steady-state model of the diffusion of solutes across the grain coating is shown to yield a novel nonlinear equation to be solved for the rate of growth of coated grains. This equation leads to a complex dependence of the growth/dissolution rate on saturation (or more generally on the composition of the fluid) in the medium surrounding the coated grain. The feedback between the dynamics of the coating thickness and morphology changes makes the phenomenon of interest here distinct from that arising from the coupling of grain growth and diffusion in the surrounding medium. This makes pits and needles possible even in a well-stirred surrounding medium, a fact of interest in interpreting the geologic record. For example, the present model can explain the development of spike, mushroom, and other outgrowths on clay-coated quartz grains in a sedimentary rock, whereas the classic Mullins and Sekerke diffusion model cannot, i.e., the time scale for eliminating concentration gradients in a pore is much shorter than that for grain growth. Predictions of the model are consistent with observations on quartz when typical values of diffusion, growth rate coefficients, and other parameters are used. In this paper, we emphasize that pitting can be closely related to outgrowth instabilities. If this be the case, then pitted dissolution of feldspars is a likely example. As feldspars dissolve, they commonly surround themselves with a clay coating, often leaving only a ghost remnant of the original growth. 相似文献
12.
We have calculated the coefficient of turbulent diffusion in a random flow with time restoration, describing the interstellar medium. Such a flow abruptly loses its memory at random times, forming a Poisson flow of events. The coefficient of turbulent diffusion in the flow is determined by the rms velocity and correlation time, as in mixing-length theory, but the numerical coefficient differs from that predicted by this theory. The closure equation derived by us for the transport of the mean concentration of a passive scalar takes a more complicated form than obtained in standard mean-field theory, but the main properties of the equation retain their validity. The possibility of extending the results of this exactly solved problem to the problem of transport in the turbulent interstellar medium is discussed. 相似文献
13.
Clarisse Alboin Jérôme Jaffré Patrick Joly Jean E. Roberts Christophe Serres 《Computational Geosciences》2002,6(3-4):523-543
Contaminant transport in a fractured porous medium can be modeled, under appropriate conditions, with a double porosity model. Such a model consists of a parabolic equation with a coupling term describing contaminant exchange between the fractures, which have high permeability, and the matrix block, which has low permeability. A locally conservative method based on mixed finite elements is used to solve the parabolic problem, and the calculation of the coupling term, which involves the solution of diffusion equations in the matrix blocks, is based on an analytic expression. Numerical experiments show that this semi-analytic method for the coupling term is accurate and faster than several other methods but at a small expense of computer memory. 相似文献
14.
变形介质流固耦合渗流的数值模型及其应用 总被引：7，自引：0，他引：7
在地下流体的开采过程中,随着地下流体的不断采出,必然造成孔隙流体压力的逐渐降低,由此导致储层岩石骨架的有效应力增大,地层产生变形或压实,甚至造成地表沉陷.针对饱和土固结问题,利用可变形多孔介质中流体渗流的流固耦合有限元数值模型,对该问题进行了数值模拟.同时,对油井开采过程进行了模拟. 相似文献
15.
Summary Temperature increase in saturated porous materials under undrained conditions leads to thermal pressurization of the pore
fluid due to the discrepancy between the thermal expansion coefficients of the pore fluid and of the solid matrix. This increase
in the pore fluid pressure induces a reduction of the effective mean stress and can lead to shear failure or hydraulic fracturing.
The equations governing the phenomenon of thermal pressurization are presented and this phenomenon is studied experimentally
for a saturated granular rock in an undrained heating test under constant isotropic stress. Careful analysis of the effect
of mechanical and thermal deformations of the drainage and pressure measurement system is performed and a correction of the
measured pore pressure is introduced. The test results are modelled using a non-linear thermo-poro-elastic constitutive model
of the granular rock with emphasis on the stress-dependent character of the rock compressibility. The effects of stress and
temperature on thermal pressurization observed in the test are correctly reproduced by the model.
Correspondence: Dr. S. Ghabezloo, CERMES, Ecole Nationale des Ponts et Chaussées, 6-8 avenue Blaise Pascal, Cité Descartes,
77455 Champs-sur-Marne, Marne la Vallée cedex 2, France 相似文献
16.
17.
深冰层取心钻探是极地科学考察的重要组成部分。在极地勘探取样过程中做到在获得高质量冰岩心样品的同时,将钻井液对极地环境的污染降低到最低程度,对保护极地原生环境具有重要意义。依据南极科学钻探的地层条件,应用Visual Modflow建立了污染物在雪层中的运移模型,在污染物浓度保持恒定时,对钻孔内污染物进行了数值模拟,预测了10年后污染物在雪中污染羽扩散范围。在此基础上,对比分析了不同弥散系数时,雪对污染物弥散的能力。极地科学钻探使用的钻井液粘度通常是水的10倍,随着污染物粘度的增大,其扩散的距离也将相应减少。研究结果为深入开展极地冰层取心钻探钻井液对极地环境影响评价研究奠定了重要的理论与技术基础。 相似文献
18.
一直以来,由Biot孔隙弹性动力方程得到的饱和土地下源Green函数都是u-w形式（u为固相介质位移,w为流相相对于固相的平均位移）。应用两相介质纵波解耦理论,得到了饱和土半空间地下点源荷载的u-P形式（P为孔压）Green函数频域解答;克服了u-w形式Green函数在边界元（BEM）积分时的增根影响。再由Hankel反演,结合Somigliana表象积分,完成BEM计算。并以计算结果分析了地下集中力作用时,饱和土位移、孔压、排水量等动力特性,这对地铁等交通工程、地震工程、土-结构动力相互作用（SSI）的响应计算都具有较重要应用价值。 相似文献
19.
层状岩体由于层理面或者结构面的存在,在力学上具有横观各向同性的特点,现有的各向同性蠕变模型难以全面反映横观各向同性岩体的蠕变力学特性。为构建能够反映横观各向同性岩体的三维蠕变模型,以能反映瞬时应变、减速蠕变和稳态蠕变特征的黏弹性Burgers模型为基础,基于常泊松比假定,在三维各向同性蠕变本构方程的基础上,按照算子替换的方法,将横观各向同性柔度矩阵代替各向同性柔度矩阵,并考虑了平行和垂直层理方向岩体蠕变力学行为的差异性,推导了横观各向同性岩体的三维蠕变本构方程。根据本构方程的特点,提出了根据平行和垂直方向岩体蠕变试验结果进行三维蠕变本构模型中蠕变参数的辨识方法。将提出的模型应用于三轴蠕变试验参数辨识,从而获得了一套完整的三维蠕变参数,并与试验结果进行对比分析,从而验证了所提模型与方法的合理性和有效性。进一步,指出了传统蠕变试验设计方案的局限性,给出了横观各向同性材料蠕变试验设计建议。研究成果为研究岩体三维蠕变机制提供了新思路,可对岩体蠕变试验设计提供相应的科研支撑。 相似文献