The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore st...The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.展开更多
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.展开更多
This article reports recent developments and advances in the simulation of the CO2-formation fluid displacement behaviour at the pore scale of subsurface porous media. Roughly, there are three effective visualization ...This article reports recent developments and advances in the simulation of the CO2-formation fluid displacement behaviour at the pore scale of subsurface porous media. Roughly, there are three effective visualization approaches to detect and observe the CO2-formation fluid displacement mechanism at the micro-scale, namely, magnetic resonance imaging, X-ray computed tomography and fabricated micromodels, but they are not capable of investigating the dis- placement process at the nano-scale. Though a lab-on-chip approach for the direct visualization of the fluid flow behaviour in nanoscale channels has been developed using an advanced epi-fluorescence microscopy method combined with a nanofluidic chip, it is still a qualitative analysis method. The lattice Boltzmann method （LBM） can simulate the CO2 displacement processes in a two-dimensional or three-dimensional （3D） pore structure, but until now, the CO2 displace- ment mechanisms had not been thoroughly investigated and the 3D pore structure of real rock had not been directly taken into account in the simulation of the CO2 displacement process. The status of research on the applications of CO2 displacement to enhance shale gas recovery is also analyzed in this paper. The coupling of molecular dynamics and LBM in tandem is proposed to simulate the CO2-shale gas displacement process based on the 3D digital model of shale obtained from focused ion beams and scanning electron microscopy.展开更多
In this paper, the mechanical response of saturated geological rock under tidal force is explored by poroelastic theory. First, we use the free energy formula of saturated rock under a tidal force to study the relatio...In this paper, the mechanical response of saturated geological rock under tidal force is explored by poroelastic theory. First, we use the free energy formula of saturated rock under a tidal force to study the relationships of pore pressure with stress, and stress with strain. Then we analyze the relationship between rock strain and tidal potential by the equilibrium differential equations of saturated rock under tidal force. Finally, we derive the physical relationship between the two parameters （pore pressure and tidal mean stress） of saturated rock and tidal potential. The relationship shows that：pore pressure is directly proportional with tidal potential, but tidal mean stress of saturated rock is inversely proportional with tidal potential. The ratio coefficient is related not only to the Lame coefficients of rock skeletons, but also to the Blot modulus. By using this model to analyze observational well water level of C-18 well which locates in Huili, Sichuan Province, the well level response coefficient （D） was estimated. This way, we derive the Skempton coefficient （B）, the coefficient A and C which refer to the response coefficients of pore pressure and tidal stress to tidal potential respectively. Then we compare the differences among each coefficient in coupling and uncoupling conditions. It shows that for saturated rocks, the response of stress and pore pressure to earth tides is a product of coupling, and it is necessary to take into account the coupling effect when we study the mechanical response. The model will provide the basis not only for the study of mechanics and hydrodynamics of well-confined aquifer systems, and the mechanics of faulting under tidal force, but also for quantitative research of the triggering mechanism of tidal forces.展开更多
基金Projects (50934002, 51074013, 51104100) supported by the National Natural Science Foundation of China
文摘The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.
基金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.
基金The authors gratefully acknowledge the financial support of the National Natural Science Foundation of China （Grant Nos. 51374213 and 51674251）, the State Key Research Development Program of China （Grant No. 2016YFC0600705）, the National Natural Science Fund for Distinguished Young Scholars of China （Grant No. 51125017）, the Fund for Innovative Research and Development Group Program of Jiangsu Province （Grant No. 2014- 27）, the Science Fund for Creative Research Groups of the National Natural Science Foundation of China （Grant No. 51421003）, and the Priority Academic Program Development of the Jiangsu Higher Education Institutions （PAPD 2014）.
文摘This article reports recent developments and advances in the simulation of the CO2-formation fluid displacement behaviour at the pore scale of subsurface porous media. Roughly, there are three effective visualization approaches to detect and observe the CO2-formation fluid displacement mechanism at the micro-scale, namely, magnetic resonance imaging, X-ray computed tomography and fabricated micromodels, but they are not capable of investigating the dis- placement process at the nano-scale. Though a lab-on-chip approach for the direct visualization of the fluid flow behaviour in nanoscale channels has been developed using an advanced epi-fluorescence microscopy method combined with a nanofluidic chip, it is still a qualitative analysis method. The lattice Boltzmann method （LBM） can simulate the CO2 displacement processes in a two-dimensional or three-dimensional （3D） pore structure, but until now, the CO2 displace- ment mechanisms had not been thoroughly investigated and the 3D pore structure of real rock had not been directly taken into account in the simulation of the CO2 displacement process. The status of research on the applications of CO2 displacement to enhance shale gas recovery is also analyzed in this paper. The coupling of molecular dynamics and LBM in tandem is proposed to simulate the CO2-shale gas displacement process based on the 3D digital model of shale obtained from focused ion beams and scanning electron microscopy.
基金supported by R&D Special Fund for Public Welfare Industry of Ministry of Science and Technology(200808055 and 200808079),the People's Republic of ChinaScience Research Plan Project of Hebei Province(Z2009104)
文摘In this paper, the mechanical response of saturated geological rock under tidal force is explored by poroelastic theory. First, we use the free energy formula of saturated rock under a tidal force to study the relationships of pore pressure with stress, and stress with strain. Then we analyze the relationship between rock strain and tidal potential by the equilibrium differential equations of saturated rock under tidal force. Finally, we derive the physical relationship between the two parameters （pore pressure and tidal mean stress） of saturated rock and tidal potential. The relationship shows that：pore pressure is directly proportional with tidal potential, but tidal mean stress of saturated rock is inversely proportional with tidal potential. The ratio coefficient is related not only to the Lame coefficients of rock skeletons, but also to the Blot modulus. By using this model to analyze observational well water level of C-18 well which locates in Huili, Sichuan Province, the well level response coefficient （D） was estimated. This way, we derive the Skempton coefficient （B）, the coefficient A and C which refer to the response coefficients of pore pressure and tidal stress to tidal potential respectively. Then we compare the differences among each coefficient in coupling and uncoupling conditions. It shows that for saturated rocks, the response of stress and pore pressure to earth tides is a product of coupling, and it is necessary to take into account the coupling effect when we study the mechanical response. The model will provide the basis not only for the study of mechanics and hydrodynamics of well-confined aquifer systems, and the mechanics of faulting under tidal force, but also for quantitative research of the triggering mechanism of tidal forces.
