One kind of 3D coupled thermo-hydro-mechanical-migratory model for saturated-unsaturated dual-porosity medium was established,in which the stress field and the temperature field are single,but the seepage field and th...One kind of 3D coupled thermo-hydro-mechanical-migratory model for saturated-unsaturated dual-porosity medium was established,in which the stress field and the temperature field are single,but the seepage field and the concentration field are double,and the influences of sets,spaces,angles,continuity ratios,stiffness of fractures on the constitutive relationship of the medium can be considered.The relative three-dimensional program of finite element method was also developed.By comparing with the existing computation example,reliability of the model and the program were verified.Taking a hypothetical nuclear waste repository as a calculation example,the radioactive nuclide leak was simulated numerically with both the rock mass and the buffer being unsaturated media,and the temperatures,negative pore pressures,flow velocities,nuclide concentrations and normal stresses in the rock mass were investigated.The results showed that the temperatures,negative pore pressures and nuclide concentrations in the buffer all present nonlinear changes and distributions that even though the saturation degree in porosity is only about 1/9 of that in fracture,the flow velocity of underground water in fracture is about 6 times of that in porosity because the permeability coefficient of fracture is almost four orders higher than that of porosity,and that the regions of stress concentration occur at the vicinity of two sides of the boundary between buffer and disposal pit wall.展开更多
Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressur...Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.展开更多
This paper compares numerical modeling of the effect of stress on solute transport （advection and matrix diffusion） in fractured rocks in which fracture apertures are correlated with fracture lengths. It is mainly m...This paper compares numerical modeling of the effect of stress on solute transport （advection and matrix diffusion） in fractured rocks in which fracture apertures are correlated with fracture lengths. It is mainly motivated by the performance and safety assessments of underground radioactive waste repositories. Five research teams used different approaches to model stress/deformation, flow and transport pro- cesses, based on either discrete fracture network or equivalent continuum models. The simulation results derived by various teams generally demonstrated that rock stresses could significantly influence solute transport processes through stress-induced changes in fracture apertures and associated changes in per- meability. Reasonably good agreement was achieved regarding advection and matrix diffusion given the same fracture network, while some observed discrepancies could be explained by different mechanical or transport modeling approaches.展开更多
Stress-dependent permeability models are developed for the organic pores and inorganic cleats/fractures in unconventional gas reservoirs,which are modeled as Biot’s porous media of dual-porosity.Further considering m...Stress-dependent permeability models are developed for the organic pores and inorganic cleats/fractures in unconventional gas reservoirs,which are modeled as Biot’s porous media of dual-porosity.Further considering multiple flow mechanisms such as dynamic effects of gas flow and surface diffusion,apparent permeability models are obtained to investigate the characteristics of unconventional gas migration.Compared to the gas transfer in single-porosity reservoirs,the gas migration ability of cleats in dual-porosity stratums rarely changes while that of organic pores is greatly improved because cleats sustain major geomechanical shrinkage deformation when the pore pressure drops.Further,the mass flux of reservoirs is dominated by the mass flux of cleats,which has a lower peak value,but a much longer production term than those in single-porosity reservoirs due to the interaction between organic pores and cleats.Parametric analysis is conducted to identify key factors significantly impacting mass flux in unconventional reservoirs.Reasons for the mass flux variation are also explored in terms of gas migration ability and pore pressure distribution.展开更多
This article presents a new well test model for stress-sensitive composite dual-porosity reservoirs based on the concept of permeability modulus, where the rock and fluid properties as well as the formation thickness ...This article presents a new well test model for stress-sensitive composite dual-porosity reservoirs based on the concept of permeability modulus, where the rock and fluid properties as well as the formation thickness vary in the radial direction. An analytical solution in the Laplace space for the pressure-transient behavior for a line-source, constant-rate well of this type of reservoir is obtained with the Laplace transformation and the perturbation technique. The pressure and its derivative in the reservoir and the effects of relevant parameters on the pressure-transient response are obtained. The model as well as the corresponding type curves may be used in predicting the production performance or analyzing the production data for this type of reservoir.展开更多
A coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturated-unsaturated ubiquitous-joint rockmass was established,in which the stress field and the temperature field were single,but the see...A coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturated-unsaturated ubiquitous-joint rockmass was established,in which the stress field and the temperature field were single,but the seepage field and the concentration field were double,and the influences of sets,spaces,angles,continuity ratios,stiffnesses of fractures on the constitutive relationship of the medium were considered.Also,the relative two-dimensional program of finite element method was developed.Taking a hypothetical nuclear waste repository as a calculation example,the case in which the rockmass was unsaturated dual-porosity medium and radioactive nuclide leak was simulated numerically,and the temperatures,negative pore pressures,saturations,flow velocities,nuclide concentrations and principal stresses in the rockmass were investigated.The results show that the negative pore pressures and nuclide concentrations in the porosity and fracture present different changes and distributions.Even though the saturation degree in porosity is only about 1/10 that in fracture,the flow velocity of underground water in fracture is about three times that in porosity because the permeability coefficient of fracture is almost four orders higher than that of porosity.The value of nuclide concentration in fracture is close to that in porosity.展开更多
With the models of stress corrosion and pressure solution, by Yasuhara et al., two computation conditions were designed for a hypothetical nuclear waste repository in an unsaturated dual-porosity rock mass: (i) the fr...With the models of stress corrosion and pressure solution, by Yasuhara et al., two computation conditions were designed for a hypothetical nuclear waste repository in an unsaturated dual-porosity rock mass: (i) the fracture apertures are changed with the stress corrosion and pressure solution (the porosity of rock matrix is also a function of stress); (ii) the fracture apertures and the porosity of rock matrix are constants. Then the corresponding two-dimensional FEM analyses for the coupled thermohydro-mechanical processes were carried out, and the states of temperatures, rates and magnitudes of aperture closure, pore pressures, flow velocities and stresses in the rock mass were investigated. The results show the followings: The aperture closure rates caused by stress corrosion are almost six orders higher than those caused by pressure solution, and the two kinds of closure rates climb up, then decline, furthermore tend towards stability; the fracture apertures decrease from the original value and tend towards the residual value, and the contact-area ratios of asperities increase also from the original value and tend towards the nominal value; when the effects of stress corrosion and pressure solution are considered, the negative fracture pressures in near field rise very highly; the differences between the magnitudes and distributions of stresses within the rock mass in two calculation cases are very small.展开更多
The Asmari Formation in the G oilfield on the Iran-Iraq border is a fractured-porous multi-lithology mixed reservoir, for which fracture is an important factor affecting oil productivity and water cut. The characteriz...The Asmari Formation in the G oilfield on the Iran-Iraq border is a fractured-porous multi-lithology mixed reservoir, for which fracture is an important factor affecting oil productivity and water cut. The characterization and modeling of fractures in the carbonate reservoir of G oilfield are challenging due to weak conventional well log responses of fractures and a lack of specific logs, such as image logs. This study proposes an integrated approach for characterizing and modeling fractures in the carbonate reservoir. The features, formation mechanism, influencing factors, and prediction methods of fractures in the Asmari Formation carbonate reservoirs of G oilfield were studied using core observation, thin section, image log, cross-dipole acoustic log (CDAL), geomechanics numerical simulation (GNS), and production data. According to CDAL-based fracture density interpretation, GNS-based fracture intensity prediction between wells, and DFN-based rock fracture properties modeling, the quantitative fracture characterization for G oilfield was realized. This research shows that the fractures in the Asamri Formation are mainly medium-to high-angle shear fractures. The substantial compression stress during the Miocene played a major role in the formation of the prominent fractures and determined their trend in the region, with primary trends of NNW-SSE and NNE-SSW. The fracture distribution has regularity, and the fractures in zone A dolomites are more highly developed than that in zone B limestones vertically. Horizontally, fractures intensity is mainly controlled by faults and structural location. The results of this study may benefit the optimization of well design during field development. From 2019 to 2021, three horizontal wells pilot tests were deployed in the fractures belt in zone A, and these fractures prominently increased the permeability of tight dolomite reservoirs. The initial production of the wells is four to five times the average production of other wells in the area, showing a good development effe展开更多
The models of stress corrosion, pressure solution and flee-face dissolution/precipitation were introduced. Taking a hypothetical nuclear waste repository in an unsaturated dual-porosity rock mass as the calculation ob...The models of stress corrosion, pressure solution and flee-face dissolution/precipitation were introduced. Taking a hypothetical nuclear waste repository in an unsaturated dual-porosity rock mass as the calculation objective, four cases were designed 1） the fracture aperture is a function of stress corrosion, pressure solution and free-face dissolution/precipitation; 2） the fracture aperture changes with stress corrosion and pressure solution; 3） the fracture aperture changes with pressure solution and free-face dissolution/precipitation; 4） the fracture aperture is only a function of pressure solution, and the matrix porosity is also a function of stress in these four cases. Then, the corresponding two-dimensional FEM analyses for the coupled thermo-hydro-mechanical processes were carried out. The results show that the effects of stress corrosion are more prominent than those of pressure solution and free-face dissolution/precipitation, and the fracture aperture and relevant permeability caused by the stress corrosion arc only about 1/5 and 1/1000 of the corresponding values created by the pressure solution and free-face dissolution/precipitation, respectively Under the action of temperature field from released heat, the negative pore and fracture pressures in the computation domain rise continuously, and are inversely proportional to the sealing of fracture aperture. The vector fields of flow velocity of fracture water in the cases with and without considering stress corrosion are obviously different. The differences between the magnitudes and distributions of stresses within the rock mass are very small in all cases.展开更多
基金supported by the National Key Basic Research and Development Program of China("973"Project)(Grant No.2010CB732101)the Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering(Grant No.SKLQ008)
文摘One kind of 3D coupled thermo-hydro-mechanical-migratory model for saturated-unsaturated dual-porosity medium was established,in which the stress field and the temperature field are single,but the seepage field and the concentration field are double,and the influences of sets,spaces,angles,continuity ratios,stiffness of fractures on the constitutive relationship of the medium can be considered.The relative three-dimensional program of finite element method was also developed.By comparing with the existing computation example,reliability of the model and the program were verified.Taking a hypothetical nuclear waste repository as a calculation example,the radioactive nuclide leak was simulated numerically with both the rock mass and the buffer being unsaturated media,and the temperatures,negative pore pressures,flow velocities,nuclide concentrations and normal stresses in the rock mass were investigated.The results showed that the temperatures,negative pore pressures and nuclide concentrations in the buffer all present nonlinear changes and distributions that even though the saturation degree in porosity is only about 1/9 of that in fracture,the flow velocity of underground water in fracture is about 6 times of that in porosity because the permeability coefficient of fracture is almost four orders higher than that of porosity,and that the regions of stress concentration occur at the vicinity of two sides of the boundary between buffer and disposal pit wall.
基金Project supported by the National Natural Science Foundation of China(No.5140232)the National Science and Technology Major Project(No.2011ZX05038003)the China Postdoctoral Science Foundation(No.2014M561074)
文摘Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.
基金the context of the international DECOVALEX-2011 ProjectLBNL from NDA via SERCO TAS was provided through the U.S. Department of Energy Contract No. DE-AC02-05CH11231supported by the Ministry of Education of the Czech Republic within the SGS project No. 7822/115 on the TUL
文摘This paper compares numerical modeling of the effect of stress on solute transport （advection and matrix diffusion） in fractured rocks in which fracture apertures are correlated with fracture lengths. It is mainly motivated by the performance and safety assessments of underground radioactive waste repositories. Five research teams used different approaches to model stress/deformation, flow and transport pro- cesses, based on either discrete fracture network or equivalent continuum models. The simulation results derived by various teams generally demonstrated that rock stresses could significantly influence solute transport processes through stress-induced changes in fracture apertures and associated changes in per- meability. Reasonably good agreement was achieved regarding advection and matrix diffusion given the same fracture network, while some observed discrepancies could be explained by different mechanical or transport modeling approaches.
基金the National Natural Science Foundation of China(Grant Nos.12102372 and 11872324)the Open Research Fund of CNMGE Platform&NSCC-TJ(Grant No.CNMGE2023011).
文摘Stress-dependent permeability models are developed for the organic pores and inorganic cleats/fractures in unconventional gas reservoirs,which are modeled as Biot’s porous media of dual-porosity.Further considering multiple flow mechanisms such as dynamic effects of gas flow and surface diffusion,apparent permeability models are obtained to investigate the characteristics of unconventional gas migration.Compared to the gas transfer in single-porosity reservoirs,the gas migration ability of cleats in dual-porosity stratums rarely changes while that of organic pores is greatly improved because cleats sustain major geomechanical shrinkage deformation when the pore pressure drops.Further,the mass flux of reservoirs is dominated by the mass flux of cleats,which has a lower peak value,but a much longer production term than those in single-porosity reservoirs due to the interaction between organic pores and cleats.Parametric analysis is conducted to identify key factors significantly impacting mass flux in unconventional reservoirs.Reasons for the mass flux variation are also explored in terms of gas migration ability and pore pressure distribution.
基金supported by the National Key Basic Research Project of China (973 Program, Grant No. 2011CB201005)the Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20105121110006)the National Science Fund for Distinguished Young Scho-lars of China(Grant No. 51125019)
文摘This article presents a new well test model for stress-sensitive composite dual-porosity reservoirs based on the concept of permeability modulus, where the rock and fluid properties as well as the formation thickness vary in the radial direction. An analytical solution in the Laplace space for the pressure-transient behavior for a line-source, constant-rate well of this type of reservoir is obtained with the Laplace transformation and the perturbation technique. The pressure and its derivative in the reservoir and the effects of relevant parameters on the pressure-transient response are obtained. The model as well as the corresponding type curves may be used in predicting the production performance or analyzing the production data for this type of reservoir.
基金Project(2010CB732101) supported by the National Basic Research Program of ChinaProject(51079145) supported by the National Natural Science Foundation of China
文摘A coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturated-unsaturated ubiquitous-joint rockmass was established,in which the stress field and the temperature field were single,but the seepage field and the concentration field were double,and the influences of sets,spaces,angles,continuity ratios,stiffnesses of fractures on the constitutive relationship of the medium were considered.Also,the relative two-dimensional program of finite element method was developed.Taking a hypothetical nuclear waste repository as a calculation example,the case in which the rockmass was unsaturated dual-porosity medium and radioactive nuclide leak was simulated numerically,and the temperatures,negative pore pressures,saturations,flow velocities,nuclide concentrations and principal stresses in the rockmass were investigated.The results show that the negative pore pressures and nuclide concentrations in the porosity and fracture present different changes and distributions.Even though the saturation degree in porosity is only about 1/10 that in fracture,the flow velocity of underground water in fracture is about three times that in porosity because the permeability coefficient of fracture is almost four orders higher than that of porosity.The value of nuclide concentration in fracture is close to that in porosity.
基金supported by the National Key Basic Research and Development Program of China ("973"Project) (Grant No. 2010CB732101)the National Natural Science Foundation of China (Grant No. 51079145)
文摘With the models of stress corrosion and pressure solution, by Yasuhara et al., two computation conditions were designed for a hypothetical nuclear waste repository in an unsaturated dual-porosity rock mass: (i) the fracture apertures are changed with the stress corrosion and pressure solution (the porosity of rock matrix is also a function of stress); (ii) the fracture apertures and the porosity of rock matrix are constants. Then the corresponding two-dimensional FEM analyses for the coupled thermohydro-mechanical processes were carried out, and the states of temperatures, rates and magnitudes of aperture closure, pore pressures, flow velocities and stresses in the rock mass were investigated. The results show the followings: The aperture closure rates caused by stress corrosion are almost six orders higher than those caused by pressure solution, and the two kinds of closure rates climb up, then decline, furthermore tend towards stability; the fracture apertures decrease from the original value and tend towards the residual value, and the contact-area ratios of asperities increase also from the original value and tend towards the nominal value; when the effects of stress corrosion and pressure solution are considered, the negative fracture pressures in near field rise very highly; the differences between the magnitudes and distributions of stresses within the rock mass in two calculation cases are very small.
基金supported by the National Science and Technology Major Project“Reservoir Characterization of Typical Thick Carbonate Reservoirs in the Middle East”(Grant No.2017ZX05032004-001).
文摘The Asmari Formation in the G oilfield on the Iran-Iraq border is a fractured-porous multi-lithology mixed reservoir, for which fracture is an important factor affecting oil productivity and water cut. The characterization and modeling of fractures in the carbonate reservoir of G oilfield are challenging due to weak conventional well log responses of fractures and a lack of specific logs, such as image logs. This study proposes an integrated approach for characterizing and modeling fractures in the carbonate reservoir. The features, formation mechanism, influencing factors, and prediction methods of fractures in the Asmari Formation carbonate reservoirs of G oilfield were studied using core observation, thin section, image log, cross-dipole acoustic log (CDAL), geomechanics numerical simulation (GNS), and production data. According to CDAL-based fracture density interpretation, GNS-based fracture intensity prediction between wells, and DFN-based rock fracture properties modeling, the quantitative fracture characterization for G oilfield was realized. This research shows that the fractures in the Asamri Formation are mainly medium-to high-angle shear fractures. The substantial compression stress during the Miocene played a major role in the formation of the prominent fractures and determined their trend in the region, with primary trends of NNW-SSE and NNE-SSW. The fracture distribution has regularity, and the fractures in zone A dolomites are more highly developed than that in zone B limestones vertically. Horizontally, fractures intensity is mainly controlled by faults and structural location. The results of this study may benefit the optimization of well design during field development. From 2019 to 2021, three horizontal wells pilot tests were deployed in the fractures belt in zone A, and these fractures prominently increased the permeability of tight dolomite reservoirs. The initial production of the wells is four to five times the average production of other wells in the area, showing a good development effe
基金Project(2010CB732101) supported by the National Basic Research Program of ChinaProject(51079145) supported by the National Natural Science Foundation of ChinaProject(2009BAK53B03) supported by the National Key Technology R&D Program of China
文摘The models of stress corrosion, pressure solution and flee-face dissolution/precipitation were introduced. Taking a hypothetical nuclear waste repository in an unsaturated dual-porosity rock mass as the calculation objective, four cases were designed 1） the fracture aperture is a function of stress corrosion, pressure solution and free-face dissolution/precipitation; 2） the fracture aperture changes with stress corrosion and pressure solution; 3） the fracture aperture changes with pressure solution and free-face dissolution/precipitation; 4） the fracture aperture is only a function of pressure solution, and the matrix porosity is also a function of stress in these four cases. Then, the corresponding two-dimensional FEM analyses for the coupled thermo-hydro-mechanical processes were carried out. The results show that the effects of stress corrosion are more prominent than those of pressure solution and free-face dissolution/precipitation, and the fracture aperture and relevant permeability caused by the stress corrosion arc only about 1/5 and 1/1000 of the corresponding values created by the pressure solution and free-face dissolution/precipitation, respectively Under the action of temperature field from released heat, the negative pore and fracture pressures in the computation domain rise continuously, and are inversely proportional to the sealing of fracture aperture. The vector fields of flow velocity of fracture water in the cases with and without considering stress corrosion are obviously different. The differences between the magnitudes and distributions of stresses within the rock mass are very small in all cases.
