Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual conne...Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.展开更多
Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recogni...Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.展开更多
The effect of roughness on flow in fractures was investigated using lattice Boltzmann method(LBM).Simulations were conducted for both statistically generated hypothetical fractures and a natural dolomite fracture. The...The effect of roughness on flow in fractures was investigated using lattice Boltzmann method(LBM).Simulations were conducted for both statistically generated hypothetical fractures and a natural dolomite fracture. The effect of increasing roughness on effective hydraulic aperture, Izbash and Forchheimer parameters with increasing Reynolds number(Re) ranging from 0.01 to 500 was examined. The growth of complex flow features, such as eddies arising near the fracture surface, was directly associated with changes in surface roughness. Rapid eddy growth above Re values of 1, followed by less rapid growth at higher Re values, suggested a three-zone nonlinear model for flow in rough fractures. This three-zone model, relating effective hydraulic conductivity to Re, was also found to be appropriate for the simulation of water flow in the natural dolomite fracture. Increasing fracture roughness led to greater eddy volumes and lower effective hydraulic conductivities for the same Re values.展开更多
Simulation of flow and transport through rough walled rock fractures is investigated using the latticeBoltzmann method （LBM） and random walk （RW）, respectively. The numerical implementation isdeveloped and validat...Simulation of flow and transport through rough walled rock fractures is investigated using the latticeBoltzmann method （LBM） and random walk （RW）, respectively. The numerical implementation isdeveloped and validated on general purpose graphic processing units （GPGPUs）. Both the LBM and RWmethod are well suited to parallel implementation on GPGPUs because they require only next-neighbourcommunication and thus can reduce expenses. The LBM model is an order of magnitude faster onGPGPUs than published results for LBM simulations run on modern CPUs. The fluid model is verified forparallel plate flow, backward facing step and single fracture flow; and the RWmodel is verified for pointsourcediffusion, Taylor-Aris dispersion and breakthrough behaviour in a single fracture. Both algorithmsplace limitations on the discrete displacement of fluid or particle transport per time step to minimise thenumerical error that must be considered during implementation. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely dis...The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.展开更多
Taking 91105 working face as the research object, the observation method of water flowing fracture<span style="font-family:Verdana;">d</span><span style="font-family:Verdana;"> zo...Taking 91105 working face as the research object, the observation method of water flowing fracture<span style="font-family:Verdana;">d</span><span style="font-family:Verdana;"> zone and the layout of mining holes were determined by analyzing the field geological structure</span><span style="font-family:Verdana;">. </span><span style="font-family:Verdana;">It was shown that the fractured zone height and the ratio given by the measured method were 52.33 and 12.46, respectively. By the numerical simulation method with the software of UDEC, the fractured zone height and the ratio were 42.5 and 10.12. By comparison of measured height data and UDEC numerical simulation, there were some differences between the measured height and the calculated results of UDEC numerical simulation method. The method of simulation can be used as the technical basis for the design of waterproof coal pillar in the future.</span>展开更多
Fluid flow in fractures controls subsurface heat and mass transport,which is essential for developing enhanced geothermal systems and radioactive waste disposal.Fracture permeability is controlled by fracture microstr...Fluid flow in fractures controls subsurface heat and mass transport,which is essential for developing enhanced geothermal systems and radioactive waste disposal.Fracture permeability is controlled by fracture microstructure(e.g.aperture,roughness,and tortuosity),but in situ values and their anisotropy have not yet been estimated.Recent advances in geophysical techniques allow the detection of changes in electrical conductivity due to changes in crustal stress and these techniques can be used to predict subsurface fluid flow.However,the paucity of data on fractured rocks hinders the quantitative interpretation of geophysical monitoring data in the field.Therefore,considering different shear displacements and chemical erosions,an investigation was conducted into the hydraulic-electric relationship as an elevated stress change in fractures.The simulation of fracture flows was achieved using the lattice Boltzmann method,while the electrical properties were calculated through the finite element method,based on synthetic faults incorporating elastic-plastic deformation.Numerical results show that the hydraulic and electrical properties depend on the rock's geometric properties(i.e.fracture length,roughness,and shear displacement).The permeability anisotropy in the direction parallel or perpendicular to the shear displacement is also notable in high stress conditions.Conversely,the permeability econductivity(i.e.,formation factor)relationship is unique under all conditions and follows a linear trend in logarithmic coordinates.However,both matrix porosity and fracture spacing alter this relationship.Both increase the slope of the linear trend,thereby changing the sensitivity of electrical observations to permeability changes.展开更多
Multiple fractured horizontal well(MFHW) is widely applied in the development of shale gas. To investigate the gas flow characteristics in shale, based on a new dual mechanism triple continuum model, an analytical sol...Multiple fractured horizontal well(MFHW) is widely applied in the development of shale gas. To investigate the gas flow characteristics in shale, based on a new dual mechanism triple continuum model, an analytical solution for MFHW surrounded by stimulated reservoir volume(SRV) was presented. Pressure and pressure derivative curves were used to identify the characteristics of flow regimes in shale. Blasingame type curves were established to evaluate the effects of sensitive parameters on rate decline curves, which indicates that the whole flow regimes could be divided into transient flow, feeding flow, and pseudo steady state flow. In feeding flow regime, the production of gas well is gradually fed by adsorbed gases in sub matrix, and free gases in matrix. The proportion of different gas sources to well production is determined by such parameters as storability ratios of triple continuum, transmissibility coefficients controlled by dual flow mechanism and fracture conductivity.展开更多
Xin’an coal mine, Henan Province, faces the risk of water inrush because 40% of the area of the coal mine is under the surface water of the Xiaolangdi reservoir. To forecast water disaster, an effective aquifuge and ...Xin’an coal mine, Henan Province, faces the risk of water inrush because 40% of the area of the coal mine is under the surface water of the Xiaolangdi reservoir. To forecast water disaster, an effective aquifuge and a limit of water infiltration were determined by rock-phase analysis and long term observations of surface water and groundwater. By field monitoring, as well as physical and numerical simulation experiments, we obtained data reflecting different heights of a water flow fractured zone (WFFZ) under different mining conditions, derived a formula to calculate this height and built a forecasting model with the aid of GIS. On the basis of these activities, the coal mine area was classified into three sub-areas with different potential of water inrush. In the end, our research results have been applied in and verified by industrial mining experiments at three working faces and we were able to present a successful example of coal mining under a large reservoir.展开更多
Even though a large number of large-scale arch dams with height larger than 200 m have been built in the world, the transient groundwater flow behaviors and the seepage control effects in the dam foundations under dif...Even though a large number of large-scale arch dams with height larger than 200 m have been built in the world, the transient groundwater flow behaviors and the seepage control effects in the dam foundations under difficult geological conditions are rarely reported. This paper presents a case study on the transient groundwater flow behaviors in the rock foundation of Jinping I double-curvature arch dam, the world＇s highest dam of this type to date that has been completed. Taking into account the geological settings at the site, an inverse modeling technique utilizing the time series measurements of both hydraulic head and discharge was adopted to back-calculate the permeability of the foundation rocks,which effectively improves the uniqueness and reliability of the inverse modeling results. The transient seepage flow in the dam foundation during the reservoir impounding was then modeled with a parabolic variational inequality（PVI） method. The distribution of pore water pressure, the amount of leakage, and the performance of the seepage control system in the dam foundation during the entire impounding process were finally illustrated with the numerical results.展开更多
The validity of Local Cubic Law （LCL） is an important issue to study groundwater flow and transport in fractured media. According to laboratory simulaion tests, the average velocity with a lower gradient in a single...The validity of Local Cubic Law （LCL） is an important issue to study groundwater flow and transport in fractured media. According to laboratory simulaion tests, the average velocity with a lower gradient in a single fracture is calculated by the LCL, which is compared with the measured average velocity. Then dye tracer test is designed and completed. The evidence for non-LCL, is drawn from the results of the simulation tests and the dye tracer tests. Then the Reynolds number of groundwater is calculated, the critical value of Re for laminar flow is discussed in a single fracture under different conditions. The motion types for groundwater flow have been discussed.展开更多
In the mining process of coal mine,waterproof coal pillars should be set between the weathered oxidation zone and the first mining face.In order to determine the reasonable upper limit of the first mining face of Hong...In the mining process of coal mine,waterproof coal pillars should be set between the weathered oxidation zone and the first mining face.In order to determine the reasonable upper limit of the first mining face of Hongyi Coal Mine,the waterproof coal pillar needs to be wide enough to resist the lateral hydrostatic pressure of the oxidation zone,and to ensure that the top plate aquifer does not run through the water guide crack zone,while also liberating as much stagnant coal as possible.In this paper,the first coal mine face’s waterproof coal pillar was calculated using conventional and optimized calculation methods,taking into account the dip angle of the coal seam,the height of the water-guide crack zone,and the rock shift influence zone,and finally the upper limit of the first face of Hongyi Coal Mine was determined.展开更多
Analyzed the rule of the Water Flowing Fractured (WFF) zone's development during the fully mechanized top coal caving.Six influence factors of WFF's height were selected,viz.mining thickness,base rock thicknes...Analyzed the rule of the Water Flowing Fractured (WFF) zone's development during the fully mechanized top coal caving.Six influence factors of WFF's height were selected,viz.mining thickness,base rock thickness,dip angle,uniaxial compressing strength of roof,mudstone proportion in overlying rock,and structure of overlying rock. The height-forecasting model of WFF was established based on the Artificial Neural Net-work techniques,and was applied in the first fully mechanized top coal caving face under sea in China.展开更多
Taking advantage of the water sensitivity of the transient electromagnetic method(TEM),this study assesses the effectiveness of the water curtain system for underground LPG storage caverns during the excavation period...Taking advantage of the water sensitivity of the transient electromagnetic method(TEM),this study assesses the effectiveness of the water curtain system for underground LPG storage caverns during the excavation period.It also detects fracture water flow during the excavation process in light of the practice of two pilot large underground LPG storage caverns in China.Comparative maps of apparent resistivity derived from TEM measurements before and after water-filling during the excavation process have been discussed to improve the quality of the water curtain system.This is the first case to apply TEM to detect the quality of the water curtain system during the construction of underground LPG storage cavern projects,and it is found to be practical,more visualized and worth popularizing.展开更多
The western Sichuan hydrothermal area is located at the northeastern margin of the eastern syntaxis of the Qinghai-Tibet Plateau, which is also the eastern end of the Mediterranean-Himalayan geothermal activity zone. ...The western Sichuan hydrothermal area is located at the northeastern margin of the eastern syntaxis of the Qinghai-Tibet Plateau, which is also the eastern end of the Mediterranean-Himalayan geothermal activity zone. There are 248 warm or hot springs in this area, and 11 have temperatures beyond the local boiling temperature. Most of these hot springs are distributed along the Jinshajiang, Dege-Xiangcheng, Ganzi-Litang, and Xianshuihe faults, forming a NW-SE hydrothermal belt. A geothermal analysis of this high-temperature hydrothermal area is an important basis for understanding the deep geodynamic process of the eastern syntaxis of the Qinghai-Tibet Plateau. In addition, this study offers an a priori view to utilize geothermal resources, which is important in both scientific research and application. We use gravity, magnetic, seismic, and helium isotope data to analyze the crust-mantle heat flow ratio and deep geothermal structure. The results show that the background terrestrial heat flow descends from southwest to northeast. The crustal heat ratio is not more than 60%. The high temperature hydrothermal active is related to crustal dynamics processes. Along the Batang-Litang-Kangding line, the Moho depth increases eastward, which is consistent with the changing Qc/Qm(crustal/mantle heat flow) ratio trend. The geoid in the hydrothermal zone is 4–6 km higher than the surroundings, forming a local "platform". The NW-SE striking local tensile stress zone and uplift structure in the upper and middle crust corresponds with the surface hydrothermal active zone. There is an average Curie Point Depth(CPD) of 19.5–22.5 km in Batang, Litang, and Kangding. The local shear-wave(S-wave) velocity is relatively low in the middle and lower crust. The S-wave shows a low velocity trap(Vs<3.2 km s.1) at 15–30 km, which is considered a high-temperature partial melting magma, the crustal source of the hydrothermal active zone. We conclude that the hydrothermal system in this area can be divided into Batang-type and Kangding-type, both of which rely on a crustal heating cycle of atmospheric precipitation and surface water along the fracture zone. The heat is derived from the middle and lower crust: groundwater penetrates the deep faults bringing geothermal energy back to the surface and forming high-temperature springs.展开更多
基金sponsored by the General Program of the National Natural Science Foundation of China(Grant Nos.52079129 and 52209148)the Hubei Provincial General Fund,China(Grant No.2023AFB567)。
文摘Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.
基金supported by the National Natural Science Foundation of China(No.52174045)。
文摘Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.
文摘The effect of roughness on flow in fractures was investigated using lattice Boltzmann method(LBM).Simulations were conducted for both statistically generated hypothetical fractures and a natural dolomite fracture. The effect of increasing roughness on effective hydraulic aperture, Izbash and Forchheimer parameters with increasing Reynolds number(Re) ranging from 0.01 to 500 was examined. The growth of complex flow features, such as eddies arising near the fracture surface, was directly associated with changes in surface roughness. Rapid eddy growth above Re values of 1, followed by less rapid growth at higher Re values, suggested a three-zone nonlinear model for flow in rough fractures. This three-zone model, relating effective hydraulic conductivity to Re, was also found to be appropriate for the simulation of water flow in the natural dolomite fracture. Increasing fracture roughness led to greater eddy volumes and lower effective hydraulic conductivities for the same Re values.
文摘Simulation of flow and transport through rough walled rock fractures is investigated using the latticeBoltzmann method （LBM） and random walk （RW）, respectively. The numerical implementation isdeveloped and validated on general purpose graphic processing units （GPGPUs）. Both the LBM and RWmethod are well suited to parallel implementation on GPGPUs because they require only next-neighbourcommunication and thus can reduce expenses. The LBM model is an order of magnitude faster onGPGPUs than published results for LBM simulations run on modern CPUs. The fluid model is verified forparallel plate flow, backward facing step and single fracture flow; and the RWmodel is verified for pointsourcediffusion, Taylor-Aris dispersion and breakthrough behaviour in a single fracture. Both algorithmsplace limitations on the discrete displacement of fluid or particle transport per time step to minimise thenumerical error that must be considered during implementation. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金ChinaCommitteeofEducation theUniver sityofArizonaandtheMetropolitanWaterDistrictofSouthernCaliforni a.
文摘The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.
文摘Taking 91105 working face as the research object, the observation method of water flowing fracture<span style="font-family:Verdana;">d</span><span style="font-family:Verdana;"> zone and the layout of mining holes were determined by analyzing the field geological structure</span><span style="font-family:Verdana;">. </span><span style="font-family:Verdana;">It was shown that the fractured zone height and the ratio given by the measured method were 52.33 and 12.46, respectively. By the numerical simulation method with the software of UDEC, the fractured zone height and the ratio were 42.5 and 10.12. By comparison of measured height data and UDEC numerical simulation, there were some differences between the measured height and the calculated results of UDEC numerical simulation method. The method of simulation can be used as the technical basis for the design of waterproof coal pillar in the future.</span>
基金supported in part by the Japan Society for the Promotion of Science (JSPS)under JSPS KAKENHI (Grant Nos.JP22K14635 and JP22H05303)a supporting program titled“Program to Support Research and Investigation on Important Basic Technologies Related to Radioactive Waste (2023 FY)”under the contract with the Ministry of Economy,Trade and Industry,Japan.
文摘Fluid flow in fractures controls subsurface heat and mass transport,which is essential for developing enhanced geothermal systems and radioactive waste disposal.Fracture permeability is controlled by fracture microstructure(e.g.aperture,roughness,and tortuosity),but in situ values and their anisotropy have not yet been estimated.Recent advances in geophysical techniques allow the detection of changes in electrical conductivity due to changes in crustal stress and these techniques can be used to predict subsurface fluid flow.However,the paucity of data on fractured rocks hinders the quantitative interpretation of geophysical monitoring data in the field.Therefore,considering different shear displacements and chemical erosions,an investigation was conducted into the hydraulic-electric relationship as an elevated stress change in fractures.The simulation of fracture flows was achieved using the lattice Boltzmann method,while the electrical properties were calculated through the finite element method,based on synthetic faults incorporating elastic-plastic deformation.Numerical results show that the hydraulic and electrical properties depend on the rock's geometric properties(i.e.fracture length,roughness,and shear displacement).The permeability anisotropy in the direction parallel or perpendicular to the shear displacement is also notable in high stress conditions.Conversely,the permeability econductivity(i.e.,formation factor)relationship is unique under all conditions and follows a linear trend in logarithmic coordinates.However,both matrix porosity and fracture spacing alter this relationship.Both increase the slope of the linear trend,thereby changing the sensitivity of electrical observations to permeability changes.
基金Project(2011ZX05015)supported by Important National Science and Technology Specific Projects of the "Twelfth Five-years" Plan Period,China
文摘Multiple fractured horizontal well(MFHW) is widely applied in the development of shale gas. To investigate the gas flow characteristics in shale, based on a new dual mechanism triple continuum model, an analytical solution for MFHW surrounded by stimulated reservoir volume(SRV) was presented. Pressure and pressure derivative curves were used to identify the characteristics of flow regimes in shale. Blasingame type curves were established to evaluate the effects of sensitive parameters on rate decline curves, which indicates that the whole flow regimes could be divided into transient flow, feeding flow, and pseudo steady state flow. In feeding flow regime, the production of gas well is gradually fed by adsorbed gases in sub matrix, and free gases in matrix. The proportion of different gas sources to well production is determined by such parameters as storability ratios of triple continuum, transmissibility coefficients controlled by dual flow mechanism and fracture conductivity.
基金Project 2007CB209400 supported by the National Basic Research Program of China
文摘Xin’an coal mine, Henan Province, faces the risk of water inrush because 40% of the area of the coal mine is under the surface water of the Xiaolangdi reservoir. To forecast water disaster, an effective aquifuge and a limit of water infiltration were determined by rock-phase analysis and long term observations of surface water and groundwater. By field monitoring, as well as physical and numerical simulation experiments, we obtained data reflecting different heights of a water flow fractured zone (WFFZ) under different mining conditions, derived a formula to calculate this height and built a forecasting model with the aid of GIS. On the basis of these activities, the coal mine area was classified into three sub-areas with different potential of water inrush. In the end, our research results have been applied in and verified by industrial mining experiments at three working faces and we were able to present a successful example of coal mining under a large reservoir.
基金financially supported through NSERC Discovery Grant(RGPIN/4994-2014)
文摘Even though a large number of large-scale arch dams with height larger than 200 m have been built in the world, the transient groundwater flow behaviors and the seepage control effects in the dam foundations under difficult geological conditions are rarely reported. This paper presents a case study on the transient groundwater flow behaviors in the rock foundation of Jinping I double-curvature arch dam, the world＇s highest dam of this type to date that has been completed. Taking into account the geological settings at the site, an inverse modeling technique utilizing the time series measurements of both hydraulic head and discharge was adopted to back-calculate the permeability of the foundation rocks,which effectively improves the uniqueness and reliability of the inverse modeling results. The transient seepage flow in the dam foundation during the reservoir impounding was then modeled with a parabolic variational inequality（PVI） method. The distribution of pore water pressure, the amount of leakage, and the performance of the seepage control system in the dam foundation during the entire impounding process were finally illustrated with the numerical results.
基金Project supported by the National Natural Science Foundation of China (Grant No :40202027) and Fok Ying Tung Ed-ucation Foundation (Grant No :91079)
文摘The validity of Local Cubic Law （LCL） is an important issue to study groundwater flow and transport in fractured media. According to laboratory simulaion tests, the average velocity with a lower gradient in a single fracture is calculated by the LCL, which is compared with the measured average velocity. Then dye tracer test is designed and completed. The evidence for non-LCL, is drawn from the results of the simulation tests and the dye tracer tests. Then the Reynolds number of groundwater is calculated, the critical value of Re for laminar flow is discussed in a single fracture under different conditions. The motion types for groundwater flow have been discussed.
基金the Open Fund of State Key Laboratory of Water Resource Protection and Utilization in Coal Mining(Grant No.GJNY-18-73.15)the National Natural Science Foundation of China(NSFC)under grant No.41602254+2 种基金the National Natural Science Foundation of China(NSFC)under grant No.41807221the Science and Technology Innovation Project of Northwest A&F University,the Double-Class Discipline Group Dry Area Hydrology and Water Resources Regulation Research funding project(Z102022011)the Special Funding Project for Basic Scientific Research Business Fees of Central Universities(2452016179).
文摘In the mining process of coal mine,waterproof coal pillars should be set between the weathered oxidation zone and the first mining face.In order to determine the reasonable upper limit of the first mining face of Hongyi Coal Mine,the waterproof coal pillar needs to be wide enough to resist the lateral hydrostatic pressure of the oxidation zone,and to ensure that the top plate aquifer does not run through the water guide crack zone,while also liberating as much stagnant coal as possible.In this paper,the first coal mine face’s waterproof coal pillar was calculated using conventional and optimized calculation methods,taking into account the dip angle of the coal seam,the height of the water-guide crack zone,and the rock shift influence zone,and finally the upper limit of the first face of Hongyi Coal Mine was determined.
基金National Science Support Plan of China(2006BAB16B04)
文摘Analyzed the rule of the Water Flowing Fractured (WFF) zone's development during the fully mechanized top coal caving.Six influence factors of WFF's height were selected,viz.mining thickness,base rock thickness,dip angle,uniaxial compressing strength of roof,mudstone proportion in overlying rock,and structure of overlying rock. The height-forecasting model of WFF was established based on the Artificial Neural Net-work techniques,and was applied in the first fully mechanized top coal caving face under sea in China.
基金We thank CREGC(20130503)Innovative Research Team in SWJTU–China(IRT0955)+1 种基金the Fundamental Research Funds for the Central Universities(SWJTU11ZT33)China Scholarship Council(201407000040)for their sponsorships of this project,and we also thank the consultants and engineers contacted,especially Jinzhou and Huizhou project departments of China Railway Erju Second Engineering Co.Ltd.for providing their data during this study.
文摘Taking advantage of the water sensitivity of the transient electromagnetic method(TEM),this study assesses the effectiveness of the water curtain system for underground LPG storage caverns during the excavation period.It also detects fracture water flow during the excavation process in light of the practice of two pilot large underground LPG storage caverns in China.Comparative maps of apparent resistivity derived from TEM measurements before and after water-filling during the excavation process have been discussed to improve the quality of the water curtain system.This is the first case to apply TEM to detect the quality of the water curtain system during the construction of underground LPG storage cavern projects,and it is found to be practical,more visualized and worth popularizing.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41574074, 41174085, 41430319)the Innovation Team Project of Chinese Academy of Sciences (Grant No. KZZD-EW-TZ-19)the Strategic Pilot Technology of Chinese Academy of Sciences (Grant No. XDA1103010102)
文摘The western Sichuan hydrothermal area is located at the northeastern margin of the eastern syntaxis of the Qinghai-Tibet Plateau, which is also the eastern end of the Mediterranean-Himalayan geothermal activity zone. There are 248 warm or hot springs in this area, and 11 have temperatures beyond the local boiling temperature. Most of these hot springs are distributed along the Jinshajiang, Dege-Xiangcheng, Ganzi-Litang, and Xianshuihe faults, forming a NW-SE hydrothermal belt. A geothermal analysis of this high-temperature hydrothermal area is an important basis for understanding the deep geodynamic process of the eastern syntaxis of the Qinghai-Tibet Plateau. In addition, this study offers an a priori view to utilize geothermal resources, which is important in both scientific research and application. We use gravity, magnetic, seismic, and helium isotope data to analyze the crust-mantle heat flow ratio and deep geothermal structure. The results show that the background terrestrial heat flow descends from southwest to northeast. The crustal heat ratio is not more than 60%. The high temperature hydrothermal active is related to crustal dynamics processes. Along the Batang-Litang-Kangding line, the Moho depth increases eastward, which is consistent with the changing Qc/Qm(crustal/mantle heat flow) ratio trend. The geoid in the hydrothermal zone is 4–6 km higher than the surroundings, forming a local "platform". The NW-SE striking local tensile stress zone and uplift structure in the upper and middle crust corresponds with the surface hydrothermal active zone. There is an average Curie Point Depth(CPD) of 19.5–22.5 km in Batang, Litang, and Kangding. The local shear-wave(S-wave) velocity is relatively low in the middle and lower crust. The S-wave shows a low velocity trap(Vs<3.2 km s.1) at 15–30 km, which is considered a high-temperature partial melting magma, the crustal source of the hydrothermal active zone. We conclude that the hydrothermal system in this area can be divided into Batang-type and Kangding-type, both of which rely on a crustal heating cycle of atmospheric precipitation and surface water along the fracture zone. The heat is derived from the middle and lower crust: groundwater penetrates the deep faults bringing geothermal energy back to the surface and forming high-temperature springs.
