The paper presents a novel hydraulic fracturing model for the characterization and simulation of the complex fracture network in shale gas reservoirs. We go beyond the existing method that uses planar or orthogonal co...The paper presents a novel hydraulic fracturing model for the characterization and simulation of the complex fracture network in shale gas reservoirs. We go beyond the existing method that uses planar or orthogonal conjugate fractures for representing the ''complexity'' of the network. Bifurcation of fractures is performed utilizing the Lindenmayer system based on fractal geometry to describe the fracture propagation pattern, density and network connectivity. Four controlling parameters are proposed to describe the details of complex fractures and stimulated reservoir volume(SRV). The results show that due to the multilevel feature of fractal fractures, the model could provide a simple method for contributing reservoir volume calibration. The primary-and second-stage fracture networks across the overall SRV are the main contributions to the production, while the induced fracture network just contributes another 20% in the late producing period. We also conduct simulation with respect to different refracturing cases and find that increasing the complexity of the fracture network provides better performance than only enhancing the fracture conductivity.展开更多
The generation method of three-dimensional fractal discrete fracture network(FDFN)based on multiplicative cascade process was developed.The complex multi-scale fracture system in shale after fracturing was characteriz...The generation method of three-dimensional fractal discrete fracture network(FDFN)based on multiplicative cascade process was developed.The complex multi-scale fracture system in shale after fracturing was characterized by coupling the artificial fracture model and the natural fracture model.Based on an assisted history matching(AHM)using multiple-proxy-based Markov chain Monte Carlo algorithm(MCMC),an embedded discrete fracture modeling(EDFM)incorporated with reservoir simulator was used to predict productivity of shale gas well.When using the natural fracture generation method,the distribution of natural fracture network can be controlled by fractal parameters,and the natural fracture network generated coupling with artificial fractures can characterize the complex system of different-scale fractures in shale after fracturing.The EDFM,with fewer grids and less computation time consumption,can characterize the attributes of natural fractures and artificial fractures flexibly,and simulate the details of mass transfer between matrix cells and fractures while reducing computation significantly.The combination of AMH and EDFM can lower the uncertainty of reservoir and fracture parameters,and realize effective inversion of key reservoir and fracture parameters and the productivity forecast of shale gas wells.Application demonstrates the results from the proposed productivity prediction model integrating FDFN,EDFM and AHM have high credibility.展开更多
After calculation on the fracture angles under various conditions of specific surface energies with different symmetry operations of rotation, the complicated behavior of dependence of fractal dimension on the structu...After calculation on the fracture angles under various conditions of specific surface energies with different symmetry operations of rotation, the complicated behavior of dependence of fractal dimension on the structure of crystal is shown. It is found that the crack propagates along the weakest crystal plane no matter what the direction of the maximum stress is if the anisotropy is sufficiently strong; and then, the fractal dimension of the fractured surfaces might be determined by the approximate fractal structure already existed in the material. Specificity of the fractal dimension of fractured surfaces would be easy to appear in this case. Reversely, the crack propagates along the direction of the maximum stress no matter what direction of the weakest crystal plane is if the anisotropy is sufficiently weak. Universality of the fractal dimension of fractured surfaces would be possible to appear in this case. In many real materials, universality and specificity of the materials are associated. The fractal dimension measured may more or less be influenced by the structure of materials and it shows its universality through the specificity of materials.展开更多
In this paper, simulated experiment device of coal and gas outburst was employed to perform the experiment on gas-containing coal extrusion. In the experiment, coal surface cracks were observed with a high-speed camer...In this paper, simulated experiment device of coal and gas outburst was employed to perform the experiment on gas-containing coal extrusion. In the experiment, coal surface cracks were observed with a high-speed camera and then the images were processed by sketch. Based on the above description, the paper studied the fractal dimension values from different positions of coal surface as well as their changing laws with time. The results show that there is a growing parabola trend of crack dimension value in the process of coal extrusion. Accordingly, we drew the conclusion that extruded coal crack evolution is a process of fractal dimension value increase. On the basis of fractal dimension values taken from different parts of coal masses, a fractal dimension of the contour map was drawn. Thus, it is clear that the contour map involves different crack fractal dimension values from different positions. To be specific, where there are complicated force and violent movement in coal mass, there are higher fractal dimension values, i.e., the further the middle of observation surface is from the exit of coal mass, and the lower the fractal dimension value is. In line with fractal geometry and energy theory of coal and gas outburst, this study presents the relation between fractal dimension and energy in the process of extruding. In conclusion, the evolution of crack fractal dimension value can signify that of energy, which has laid a solid foundation for the quantification research on the mechanism of gas-containing coal extrusion.展开更多
In this paper a two dimensional (2 D) model of earthquake fault rupturing was presented. It was estabilished on the basis of 1 D spring block model. Using this model, we studied the dynamical plane strain fractur...In this paper a two dimensional (2 D) model of earthquake fault rupturing was presented. It was estabilished on the basis of 1 D spring block model. Using this model, we studied the dynamical plane strain fracture problem, modeled the whole dynamical process of nucleating, expanding and propagating of fracture on a 2 D fault with homogeneous or inhomogeneous rupture strength distribution. Our studies show that under homogeneous prestress condition, the fault will gain enough momentum to tear strong obstacles in their propagating path. The rupturing fronts can also propagate forth around the isolated barriers. It is shown that the stopping conditions for rupturing processes play an important role in modeling whole earthquake process. We also studied the dynamical rupturing problems of the fault on which the rupture strength distribution is inhomogeneous, and modeled the earthquake sequence generated on a 2 D fault with the strength distribution of fractal structure. It possesses some similar features as a seismic sequence in the nature. These features mainly depend on the distribution of rupture strength on the fault plane and the level of initial stress drop. The modeling studies which were established on the basis of experiments and observations provided the physical basis for explaining some statistical rules of seismicity.展开更多
基金supported by National Natural Science Foundation of China(No.51674279)China Postdoctoral Science Foundation(No.2016M602227)a grant from National Science and Technology Major Project(No.2017ZX05049-006)
文摘The paper presents a novel hydraulic fracturing model for the characterization and simulation of the complex fracture network in shale gas reservoirs. We go beyond the existing method that uses planar or orthogonal conjugate fractures for representing the ''complexity'' of the network. Bifurcation of fractures is performed utilizing the Lindenmayer system based on fractal geometry to describe the fracture propagation pattern, density and network connectivity. Four controlling parameters are proposed to describe the details of complex fractures and stimulated reservoir volume(SRV). The results show that due to the multilevel feature of fractal fractures, the model could provide a simple method for contributing reservoir volume calibration. The primary-and second-stage fracture networks across the overall SRV are the main contributions to the production, while the induced fracture network just contributes another 20% in the late producing period. We also conduct simulation with respect to different refracturing cases and find that increasing the complexity of the fracture network provides better performance than only enhancing the fracture conductivity.
基金Supported by the National Science and Technology Major Project(2017ZX05063-005)Science and Technology Development Project of PetroChina Research Institute of Petroleum Exploration and Development(YGJ2019-12-04)。
文摘The generation method of three-dimensional fractal discrete fracture network(FDFN)based on multiplicative cascade process was developed.The complex multi-scale fracture system in shale after fracturing was characterized by coupling the artificial fracture model and the natural fracture model.Based on an assisted history matching(AHM)using multiple-proxy-based Markov chain Monte Carlo algorithm(MCMC),an embedded discrete fracture modeling(EDFM)incorporated with reservoir simulator was used to predict productivity of shale gas well.When using the natural fracture generation method,the distribution of natural fracture network can be controlled by fractal parameters,and the natural fracture network generated coupling with artificial fractures can characterize the complex system of different-scale fractures in shale after fracturing.The EDFM,with fewer grids and less computation time consumption,can characterize the attributes of natural fractures and artificial fractures flexibly,and simulate the details of mass transfer between matrix cells and fractures while reducing computation significantly.The combination of AMH and EDFM can lower the uncertainty of reservoir and fracture parameters,and realize effective inversion of key reservoir and fracture parameters and the productivity forecast of shale gas wells.Application demonstrates the results from the proposed productivity prediction model integrating FDFN,EDFM and AHM have high credibility.
基金National Natural Science Foundation of China!59671093 National Natural Science Foundation of China !19874064
文摘After calculation on the fracture angles under various conditions of specific surface energies with different symmetry operations of rotation, the complicated behavior of dependence of fractal dimension on the structure of crystal is shown. It is found that the crack propagates along the weakest crystal plane no matter what the direction of the maximum stress is if the anisotropy is sufficiently strong; and then, the fractal dimension of the fractured surfaces might be determined by the approximate fractal structure already existed in the material. Specificity of the fractal dimension of fractured surfaces would be easy to appear in this case. Reversely, the crack propagates along the direction of the maximum stress no matter what direction of the weakest crystal plane is if the anisotropy is sufficiently weak. Universality of the fractal dimension of fractured surfaces would be possible to appear in this case. In many real materials, universality and specificity of the materials are associated. The fractal dimension measured may more or less be influenced by the structure of materials and it shows its universality through the specificity of materials.
基金the National Natural Science Foundation of China (Nos. 50904067 and 51104156)the New Century Excellent Talents in University (No. NCET-10-0768) for their support of this project
文摘In this paper, simulated experiment device of coal and gas outburst was employed to perform the experiment on gas-containing coal extrusion. In the experiment, coal surface cracks were observed with a high-speed camera and then the images were processed by sketch. Based on the above description, the paper studied the fractal dimension values from different positions of coal surface as well as their changing laws with time. The results show that there is a growing parabola trend of crack dimension value in the process of coal extrusion. Accordingly, we drew the conclusion that extruded coal crack evolution is a process of fractal dimension value increase. On the basis of fractal dimension values taken from different parts of coal masses, a fractal dimension of the contour map was drawn. Thus, it is clear that the contour map involves different crack fractal dimension values from different positions. To be specific, where there are complicated force and violent movement in coal mass, there are higher fractal dimension values, i.e., the further the middle of observation surface is from the exit of coal mass, and the lower the fractal dimension value is. In line with fractal geometry and energy theory of coal and gas outburst, this study presents the relation between fractal dimension and energy in the process of extruding. In conclusion, the evolution of crack fractal dimension value can signify that of energy, which has laid a solid foundation for the quantification research on the mechanism of gas-containing coal extrusion.
文摘In this paper a two dimensional (2 D) model of earthquake fault rupturing was presented. It was estabilished on the basis of 1 D spring block model. Using this model, we studied the dynamical plane strain fracture problem, modeled the whole dynamical process of nucleating, expanding and propagating of fracture on a 2 D fault with homogeneous or inhomogeneous rupture strength distribution. Our studies show that under homogeneous prestress condition, the fault will gain enough momentum to tear strong obstacles in their propagating path. The rupturing fronts can also propagate forth around the isolated barriers. It is shown that the stopping conditions for rupturing processes play an important role in modeling whole earthquake process. We also studied the dynamical rupturing problems of the fault on which the rupture strength distribution is inhomogeneous, and modeled the earthquake sequence generated on a 2 D fault with the strength distribution of fractal structure. It possesses some similar features as a seismic sequence in the nature. These features mainly depend on the distribution of rupture strength on the fault plane and the level of initial stress drop. The modeling studies which were established on the basis of experiments and observations provided the physical basis for explaining some statistical rules of seismicity.