A numerical procedure for the evaluation of equivalent permeability tensor for fractured vuggy porous media is presented.At first we proposed a new conceptual model,i.e.,discrete fracture-vug network model,to model th...A numerical procedure for the evaluation of equivalent permeability tensor for fractured vuggy porous media is presented.At first we proposed a new conceptual model,i.e.,discrete fracture-vug network model,to model the realistic fluid flow in fractured vuggy porous medium on fine scale.This new model consists of three systems:rock matrix system,fractures system,and vugs system.The fractures and vugs are embedded in porous rock,and the isolated vugs could be connected via discrete fracture network.The flow in porous rock and fractures follows Darcy’s law,and the vugs system is free fluid region.Based on two-scale homogenization theory,we obtained an equivalent macroscopic Darcy’s law on coarse scale from fine-scale discrete fracture-vug network model.A finite element numerical formulation for homogenization equations is developed.The method is verified through application to a periodic model problem and then is applied to the calculation of equivalent permeability tensor of porous media with complex fracture-vug networks.The applicability and validity of the method for these more general fractured vuggy systems are assessed through a simple test of the coarse-scale model.展开更多
The traditional Green-Ampt model does not accurately represent the infiltration behavior of clay soils.Infiltration in clay is influenced by low hydraulic conductivity,strong capillary forces,and a gradual transition ...The traditional Green-Ampt model does not accurately represent the infiltration behavior of clay soils.Infiltration in clay is influenced by low hydraulic conductivity,strong capillary forces,and a gradual transition zone between saturated and unsaturated zones.These factors often lead to overestimated infiltration rates and underestimated infiltration durations.Therefore,it is necessary to improve the model to better reflect the characteristics of clay infiltration and enhance its predictive accuracy and practical applicability.This study conducts hydraulic characterization tests,one-dimensional soil column rainfall infiltration experiments,and numerical analysis on a representative clay sampled from Wuhan,China,to investigate infiltration behaviors under varying rainfall intensities and initial moisture conditions.The study reveals that the proportion of the transition layer within the wetting layer decreases with increasing wetting front depth,following a power-law function.Under the same initial moisture content,this proportion tends to converge to a stable value regardless of rainfall intensity.In contrast,under the same rainfall intensity,a higher initial moisture content leads to a larger proportion of the transition layer at a given wetting front depth.Based on the NMR curve,the unsaturated permeability coefficients corresponding to different volumetric water contents of clay can be obtained quickly,accurately,and at low cost.By utilizing the unsaturated permeability coefficient prediction model based on the nuclear magnetic resonance(NMR)curve,the study refines the computational method for the equivalent permeability coefficient in the wetting layer during clay rainfall infiltration,and proposes an improved clay Green-Ampt infiltration model that considers the saturated-unsaturated differentiation layer and the dynamic variation of its equivalent permeability coefficient under continuous rainfall conditions.The computational results of the improved model were compared with measured infiltration data,numerical simulations,and predictions from the traditional GA model.The results indicate that the improved model effectively captures the dynamic variation between the transition layer and wetting layer and provides more accurate predictions of wetting front depth in clay,with an accuracy approximately 68.36%higher than that of the traditional GA model.展开更多
The ability to capture permeability of fractured porous media plays a significant role in several engineering applications, including reservoir, mining, petroleum and geotechnical engineering. In order to solve fluid ...The ability to capture permeability of fractured porous media plays a significant role in several engineering applications, including reservoir, mining, petroleum and geotechnical engineering. In order to solve fluid flow and coupled flow-deformation problems encountered in these engineering applications,both empirical and theoretical models had been proposed in the past few decades. Some of them are simple but still work in certain circumstances; others are complex but also need some modifications to be applicable. Thus, the understanding of state-of-the-art permeability evolution model would help researchers and engineers solve engineering problems through an appropriate approach. This paper summarizes permeability evolution models proposed by earlier and recent researchers with emphasis on their characteristics and limitations.展开更多
Based on the characteristics of fractured vuggy porous media,a novel mathematical model was proposed to model fluid flow in such media on fine scale,i.e.,the discrete fracture-vug network model.The new model consists ...Based on the characteristics of fractured vuggy porous media,a novel mathematical model was proposed to model fluid flow in such media on fine scale,i.e.,the discrete fracture-vug network model.The new model consists of three systems:porous rock system,fracture system,and vug system.The fractures and vugs are embedded in porous rock,and the isolated vugs could be connected via the discrete fracture network.The flow in porous rock and fractures follows Darcy's law,and the vugs system is free fluid region.Using a two-scale homogenization limit theory,we obtained a macroscopic Darcy's law governing the media on coarse scale.The theoretical formula of the equivalent permeability of the fractured vuggy porous media was derived.The model and method of this paper were verified by some numerical examples.At the end the permeability of some fractured vuggy porous media with typical fracture-vug structures was analyzed.展开更多
The applicability of cement grout （or cement-based grout） has been considered as an alternative to bentonite grout commonly used to backfill closed-loop vertical ground heat exchangers. In a geothermal heat pump sys...The applicability of cement grout （or cement-based grout） has been considered as an alternative to bentonite grout commonly used to backfill closed-loop vertical ground heat exchangers. In a geothermal heat pump system, repeated heating-cooling cycles may cause adverse effects on the integrity of cement grout in the ground heat exchanger. To account for the temperature cycling effect, the strength degradation of cement grout due to temperature cycling has been examined by measuring the unconfined compression strength of cured specimens in a humidity-temperature controlling chamber with applying temperature cycles between -5℃ and 50℃. There is a tendency that the unconfined compression strength decreases with an increase in the number of temperature cycles. On the other hand, an equivalent hydraulic conductivity of a pipe-embedded cement grout specimen was evaluated by carrying out a modified flexible wall permeameter test equipped with a water circulating system to control temperature inside the pipe section. The applied operating temperature range was from 5 to 35℃. After three cycles of heating-cooling circulation, the equivalent hydraulic conductivity becomes asymptotic to a constant value, which implies there is no severe detachment of the pipe from the cement grout.展开更多
During the long-time operation of salt rock storage cavern,between its formations,damaged interfaces induced by discontinuous creep deformations between adjacent layers will possibly lead to serious gas leakage.In thi...During the long-time operation of salt rock storage cavern,between its formations,damaged interfaces induced by discontinuous creep deformations between adjacent layers will possibly lead to serious gas leakage.In this paper,damaged interfaces are considered as main potential leakage path:firstly in meso-level,gas flow rule along the interface is analyzed and the calculation of equivalent permeability is discussed.Then based on porous media seepage theory,gas leakage simulation model including salt rock,cavity interlayers and interface is built.With this strategy,it is possible to overcome the disadvantage of simulation burden with porous-fractured double medium.It also can provide the details of gas flowing along the damaged zones.Finally this proposal is applied to the salt cavern in Qianjian mines(East China).Under different operation pressures,gas distributions around two adjacent cavities are simulated;the evolvement of gas in the interlayers and salt rock is compared.From the results it is demonstrated that the domain of creep damage area has great influence on leakage range.And also the leakage in the interface will accelerate the development of leakage in salt rock.It is concluded that compared with observations,this new strategy provides closer answers.The simulation result proves its validity for the design and reasonable control of operating pressure and tightness evaluation of group bedded salt rock storage caverns.展开更多
文摘A numerical procedure for the evaluation of equivalent permeability tensor for fractured vuggy porous media is presented.At first we proposed a new conceptual model,i.e.,discrete fracture-vug network model,to model the realistic fluid flow in fractured vuggy porous medium on fine scale.This new model consists of three systems:rock matrix system,fractures system,and vugs system.The fractures and vugs are embedded in porous rock,and the isolated vugs could be connected via discrete fracture network.The flow in porous rock and fractures follows Darcy’s law,and the vugs system is free fluid region.Based on two-scale homogenization theory,we obtained an equivalent macroscopic Darcy’s law on coarse scale from fine-scale discrete fracture-vug network model.A finite element numerical formulation for homogenization equations is developed.The method is verified through application to a periodic model problem and then is applied to the calculation of equivalent permeability tensor of porous media with complex fracture-vug networks.The applicability and validity of the method for these more general fractured vuggy systems are assessed through a simple test of the coarse-scale model.
基金financial support from the Joint Funds of the National Nature Science Foundation of China(No.U22A20232)Supported by Open Project Funding of Key Laboratory of Intelligent Health Perception and Ecological Restoration of Rivers and Lakes,Ministry of Education(HGKFZ07)+2 种基金the National Natural Science Foundation of China(No.51978249)Innovation Research Team Project of the Hubei Provincial Department of Science and Technology(2025AFA020)the International Collaborative Research Fund for Young Scholars in the Innovation Demonstration Base of Ecological Environment Geotechnical and Ecological Restoration of Rivers and Lakes.
文摘The traditional Green-Ampt model does not accurately represent the infiltration behavior of clay soils.Infiltration in clay is influenced by low hydraulic conductivity,strong capillary forces,and a gradual transition zone between saturated and unsaturated zones.These factors often lead to overestimated infiltration rates and underestimated infiltration durations.Therefore,it is necessary to improve the model to better reflect the characteristics of clay infiltration and enhance its predictive accuracy and practical applicability.This study conducts hydraulic characterization tests,one-dimensional soil column rainfall infiltration experiments,and numerical analysis on a representative clay sampled from Wuhan,China,to investigate infiltration behaviors under varying rainfall intensities and initial moisture conditions.The study reveals that the proportion of the transition layer within the wetting layer decreases with increasing wetting front depth,following a power-law function.Under the same initial moisture content,this proportion tends to converge to a stable value regardless of rainfall intensity.In contrast,under the same rainfall intensity,a higher initial moisture content leads to a larger proportion of the transition layer at a given wetting front depth.Based on the NMR curve,the unsaturated permeability coefficients corresponding to different volumetric water contents of clay can be obtained quickly,accurately,and at low cost.By utilizing the unsaturated permeability coefficient prediction model based on the nuclear magnetic resonance(NMR)curve,the study refines the computational method for the equivalent permeability coefficient in the wetting layer during clay rainfall infiltration,and proposes an improved clay Green-Ampt infiltration model that considers the saturated-unsaturated differentiation layer and the dynamic variation of its equivalent permeability coefficient under continuous rainfall conditions.The computational results of the improved model were compared with measured infiltration data,numerical simulations,and predictions from the traditional GA model.The results indicate that the improved model effectively captures the dynamic variation between the transition layer and wetting layer and provides more accurate predictions of wetting front depth in clay,with an accuracy approximately 68.36%higher than that of the traditional GA model.
基金supported by the National Nature Science Foundation of China(No.51278383,No.51238009 and No.51025827)Key Scientific and Technological Innovation Team of Zhejiang Province(No.2011R50020)Key Scientific and Technological Innovation Team of Wenzhou(No.C20120006)
文摘The ability to capture permeability of fractured porous media plays a significant role in several engineering applications, including reservoir, mining, petroleum and geotechnical engineering. In order to solve fluid flow and coupled flow-deformation problems encountered in these engineering applications,both empirical and theoretical models had been proposed in the past few decades. Some of them are simple but still work in certain circumstances; others are complex but also need some modifications to be applicable. Thus, the understanding of state-of-the-art permeability evolution model would help researchers and engineers solve engineering problems through an appropriate approach. This paper summarizes permeability evolution models proposed by earlier and recent researchers with emphasis on their characteristics and limitations.
基金supported by the National Basic Research Program of China("973"Program)(Grant No.2006CB202404)
文摘Based on the characteristics of fractured vuggy porous media,a novel mathematical model was proposed to model fluid flow in such media on fine scale,i.e.,the discrete fracture-vug network model.The new model consists of three systems:porous rock system,fracture system,and vug system.The fractures and vugs are embedded in porous rock,and the isolated vugs could be connected via the discrete fracture network.The flow in porous rock and fractures follows Darcy's law,and the vugs system is free fluid region.Using a two-scale homogenization limit theory,we obtained a macroscopic Darcy's law governing the media on coarse scale.The theoretical formula of the equivalent permeability of the fractured vuggy porous media was derived.The model and method of this paper were verified by some numerical examples.At the end the permeability of some fractured vuggy porous media with typical fracture-vug structures was analyzed.
基金supported by the Fundamental Research and Development Program of the Center of New and Renewable Energy of the Ministry of Knowledge and Economy (Grant No. 2008-N-GE08-R-01)the National Research Foundation of Korea Grant funded by the Korean Government (Grant No. 2010-0011159)
文摘The applicability of cement grout （or cement-based grout） has been considered as an alternative to bentonite grout commonly used to backfill closed-loop vertical ground heat exchangers. In a geothermal heat pump system, repeated heating-cooling cycles may cause adverse effects on the integrity of cement grout in the ground heat exchanger. To account for the temperature cycling effect, the strength degradation of cement grout due to temperature cycling has been examined by measuring the unconfined compression strength of cured specimens in a humidity-temperature controlling chamber with applying temperature cycles between -5℃ and 50℃. There is a tendency that the unconfined compression strength decreases with an increase in the number of temperature cycles. On the other hand, an equivalent hydraulic conductivity of a pipe-embedded cement grout specimen was evaluated by carrying out a modified flexible wall permeameter test equipped with a water circulating system to control temperature inside the pipe section. The applied operating temperature range was from 5 to 35℃. After three cycles of heating-cooling circulation, the equivalent hydraulic conductivity becomes asymptotic to a constant value, which implies there is no severe detachment of the pipe from the cement grout.
基金We acknowledge the following funds to give financial supports.They are China National Program on National natural sciences foundation of China Grant no.51104108 and 41172284,Key Basic Research Project(973 Program)Grant no.2009CB724603.
文摘During the long-time operation of salt rock storage cavern,between its formations,damaged interfaces induced by discontinuous creep deformations between adjacent layers will possibly lead to serious gas leakage.In this paper,damaged interfaces are considered as main potential leakage path:firstly in meso-level,gas flow rule along the interface is analyzed and the calculation of equivalent permeability is discussed.Then based on porous media seepage theory,gas leakage simulation model including salt rock,cavity interlayers and interface is built.With this strategy,it is possible to overcome the disadvantage of simulation burden with porous-fractured double medium.It also can provide the details of gas flowing along the damaged zones.Finally this proposal is applied to the salt cavern in Qianjian mines(East China).Under different operation pressures,gas distributions around two adjacent cavities are simulated;the evolvement of gas in the interlayers and salt rock is compared.From the results it is demonstrated that the domain of creep damage area has great influence on leakage range.And also the leakage in the interface will accelerate the development of leakage in salt rock.It is concluded that compared with observations,this new strategy provides closer answers.The simulation result proves its validity for the design and reasonable control of operating pressure and tightness evaluation of group bedded salt rock storage caverns.
