To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model ...To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model coupled with flow,stress and damage is proposed.A traction-separation law is used to describe the mixed-mode response of the damaged adhesive fractures,and the cubic law is used to describe the fluid flow within the fractures.The rock deformation is controlled by the in-situ stress,fracture cohesion and fluid pressure on the hydraulic fracture surface.The coupled finite element equations are solved by the explicit time difference method.The effects of the fracturing treatment parameters including fluid viscosity,pumping rate and cluster spacing on the geometries of multifractures are investigated.The results show that variable fluid-viscosity injection can improve the complexity of the fracture network and height of the main fractures simultaneously.The pumping rate of15 m^(3)/min,variable fluid-viscosity of 3-9-21-36-45 mPa s with a cluster spacing of 7.5 m is the ideal treatment strategy.The field application shows that the peak daily production of the application well with the optimized injection procedu re of variable fluid-viscosity fracturing is 171 tons(about 2.85 times that of the adjacent well),which is the highest daily production record of a single shale oil well in China,marking a strategic breakthrough of commercial shale oil production in the Jiyang Depression,Shengli Oilfield.The variable fluid-viscosity fracturing technique is proved to be very effective for improving shale oil production.展开更多
It is well known that the complicated channeling of fluid flow and heat transfer is strongly related with the intricate natural fracture system.However,it is still challenging to set up the fracture network model whic...It is well known that the complicated channeling of fluid flow and heat transfer is strongly related with the intricate natural fracture system.However,it is still challenging to set up the fracture network model which is strong heterogeneous.Compared with other methods(e.g.equivalent continuum model(ECM),discrete fracture model(DFM),and ECM-DFM),the fracture flow module in the COMSOL Multiphysics simulator is powerful in definition of fractures as the inner flow boundary existing in the porous media.Thus it is selected to simulate the fluid flow and heat transfer in the geothermal-developed fractured granite of Sanguliu area located at Liaodong Peninsula,Eastern China.The natural faults/fractures based on field investigation combined with the discrete fracture network(DFN)generated by the MATLAB are used to represent the two-dimensional geological model.Numerical results show that early thermal breakthrough occurs at the production well caused by quick flow of cold water along the highly connected fractures.Suitable hydraulic fracturing treatments with proper injection rates,locations,etc.can efficiently hinder the thermal breakthrough time in the natural fracture system.Large well spacing helps the long-term operation of geothermal production,but it is highly dependent on the geometrical morphology of the fracture network.The enhancement of reservoir properties at the near-well regions can also increase the geothermal production efficiency.The results in this study can provide references to achieve a sustainable geothermal exploitation in fractured granitic geothermal reservoirs or hot dry rocks at depth.展开更多
Characterizing natural fractures has a decisive effect on production forecasts in fractured oil and gas reservoirs.Discrete Fracture Networks(DFN)constitutes the main modeling framework for fractured geosystems.Howeve...Characterizing natural fractures has a decisive effect on production forecasts in fractured oil and gas reservoirs.Discrete Fracture Networks(DFN)constitutes the main modeling framework for fractured geosystems.However,myriads of uncertainties are enclosed prior modeling representative stochastic or deterministic DFN ensembles.This paper presents a novel methodology for DFN calibration and an efficacious field application,which incorporatesWell-testing interpretation,Embedded Discrete Fracture Model(EDFM)framework,and numerical reservoir simulation.The proposed workflow starts with the DFN generation from seismic data,imaging logging data and core data.After multiple DFNs are modeled,well-test analysis is employed to calibrate the intrinsic properties of fractures at different locations.Then,these fracture networks are characterized dynamically by EDFM,which promotes capturing the optimal fracture model quickly screened.Finally,pressure and production history match are reached for the DFN realization that honors the optimal fracture model.展开更多
In this paper, the effect of pre-existing discrete fracture network(DFN) connectivity on hydraulic fracturing is numerically investigated in a rock mass subjected to in-situ stress. The simulation results show that DF...In this paper, the effect of pre-existing discrete fracture network(DFN) connectivity on hydraulic fracturing is numerically investigated in a rock mass subjected to in-situ stress. The simulation results show that DFN connectivity has a significant influence on the hydraulic fracture(HF) & DFN interaction and hydraulic fracturing effectiveness, which can be characterized by the total interaction area, stimulated DFN length, stimulated HF length, leak-off ratio, and stimulated total length. In addition, even at the same fluid injection rate, simulation models exhibit different responses that are strongly affected by the DFN connectivity. At a low injection rate, total interaction area decreases with increasing DFN connectivity; at a high injection rate, total interaction area increases with the increase of DFN connectivity. However, for any injection rate, the stimulated DFN length increases and stimulated HF length decreases with the increase of connectivity. Generally, this work shows that the DFN connectivity plays a crucial role in the interaction between hydraulic fractures, the pre-existing natural fractures and hydraulic fracturing effectiveness; in return, these three factors affect treating pressure, created microseismicity and corresponding stimulated volume. This work strongly relates to the production technology and the evaluation of hydraulic fracturing effectiveness. It is helpful for the optimization of hydraulic fracturing simulations in naturally fractured formations.展开更多
Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-ex...Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-existing joints in the damage evolution around the underground opening is of critical importance as they govern the fracturing mechanisms and influence the brittle responses of these hard rock masses under highly anisotropic in situ stresses.In this study,the main focus is the impact of joint network geometry,joint strength and applied field stresses on the rock mass behaviours and the evolution of excavation induced damage due to the loss of confinement as a tunnel face advances.Analysis of such a phenomenon was conducted using the finite-discrete element method (FDEM).The numerical model is initially calibrated in order to match the behaviour of the fracture-free,massive Lac du Bonnet granite during the excavation of the Underground Research Laboratory (URL) Test Tunnel,Canada.The influence of the pre-existing joints on the rock mass response during excavation is investigated by integrating discrete fracture networks (DFNs) of various characteristics into the numerical models under varying in situ stresses.The numerical results obtained highlight the significance of the pre-existing joints on the reduction of in situ rock mass strength and its capacity for extension with both factors controlling the brittle response of the material.Furthermore,the impact of spatial distribution of natural joints on the stability of an underground excavation is discussed,as well as the potentially minor influence of joint strength on the stress induced damage within joint systems of a non-persistent nature under specific conditions.Additionally,the in situ stress-joint network interaction is examined,revealing the complex fracturing mechanisms that may lead to uncontrolled fracture propagation that compromises the overall stability of an underground excavation.展开更多
The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks（well interlocked undisturbed rock mass with blocks formed by three or ...The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks（well interlocked undisturbed rock mass with blocks formed by three or less intersecting joints）. A synthetic rock mass modelling（SRM） approach is employed to determine the mechanical properties of the rock mass. In this approach, the intact body of rock is represented by the discrete element method（DEM）-Voronoi grains with the ability of simulating the initiation and propagation of microcracks within the intact part of the model. The geometry of the preexisting joints is generated by employing discrete fracture network（DFN） modelling based on field joint data collected from the Brockville Tunnel using LiDAR scanning. The geometrical characteristics of the simulated joints at a representative sample size are first validated against the field data, and then used to measure the rock quality designation（RQD）, joint spacing, areal fracture intensity（P21）, and block volumes. These geometrical quantities are used to quantitatively determine a representative range of the geological strength index（GSI）. The results show that estimating the GSI using the RQD tends to make a closer estimate of the degree of blockiness that leads to GSI values corresponding to those obtained from direct visual observations of the rock mass conditions in the field. The use of joint spacing and block volume in order to quantify the GSI value range for the studied rock mass suggests a lower range compared to that evaluated in situ. Based on numerical modelling results and laboratory data of rock testing reported in the literature, a semi-empirical equation is proposed that relates the rock mass modulus to confinement as a function of the areal fracture intensity and joint stiffness.展开更多
Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlight...Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.展开更多
Rockbursting in deep tunnelling is a complex phenomenon posing significant challenges both at the design and construction stages of an underground excavation within hard rock masses and under high in situ stresses. Wh...Rockbursting in deep tunnelling is a complex phenomenon posing significant challenges both at the design and construction stages of an underground excavation within hard rock masses and under high in situ stresses. While local experience, field monitoring, and informed data-rich analysis are some of the tools commonly used to manage the hazards and the associated risks, advanced numerical techniques based on discontinuum modelling have also shown potential in assisting in the assessment of rockbursting. In this study, the hybrid finite-discrete element method(FDEM) is employed to investigate the failure and fracturing processes, and the mechanisms of energy storage and rapid release resulting in bursting, as well as to assess its utility as part of the design process of underground excavations.Following the calibration of the numerical model to simulate a deep excavation in a hard, massive rock mass, discrete fracture network(DFN) geometries are integrated into the model in order to examine the impact of rock structure on rockbursting under high in situ stresses. The obtained analysis results not only highlight the importance of explicitly simulating pre-existing joints within the model, as they affect the mobilised failure mechanisms and the intensity of strain bursting phenomena, but also show how the employed joint network geometry, the field stress conditions, and their interaction influence the extent and depth of the excavation induced damage. Furthermore, a rigorous analysis of the mass and velocity of the ejected rock blocks and comparison of the obtained data with well-established semi-empirical approaches demonstrate the potential of the method to provide realistic estimates of the kinetic energy released during bursting for determining the energy support demand.展开更多
基金funded by the National Natural Science Foundation of China(Nos.52192622,51874253,U20A202)
文摘To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model coupled with flow,stress and damage is proposed.A traction-separation law is used to describe the mixed-mode response of the damaged adhesive fractures,and the cubic law is used to describe the fluid flow within the fractures.The rock deformation is controlled by the in-situ stress,fracture cohesion and fluid pressure on the hydraulic fracture surface.The coupled finite element equations are solved by the explicit time difference method.The effects of the fracturing treatment parameters including fluid viscosity,pumping rate and cluster spacing on the geometries of multifractures are investigated.The results show that variable fluid-viscosity injection can improve the complexity of the fracture network and height of the main fractures simultaneously.The pumping rate of15 m^(3)/min,variable fluid-viscosity of 3-9-21-36-45 mPa s with a cluster spacing of 7.5 m is the ideal treatment strategy.The field application shows that the peak daily production of the application well with the optimized injection procedu re of variable fluid-viscosity fracturing is 171 tons(about 2.85 times that of the adjacent well),which is the highest daily production record of a single shale oil well in China,marking a strategic breakthrough of commercial shale oil production in the Jiyang Depression,Shengli Oilfield.The variable fluid-viscosity fracturing technique is proved to be very effective for improving shale oil production.
基金financial support from the projects of the National Natural Science Foundation of China(NSFC)(Grant Nos.51809259,51774056,and 51774095)the CAS Pioneer Hundred Talents Program in China。
文摘It is well known that the complicated channeling of fluid flow and heat transfer is strongly related with the intricate natural fracture system.However,it is still challenging to set up the fracture network model which is strong heterogeneous.Compared with other methods(e.g.equivalent continuum model(ECM),discrete fracture model(DFM),and ECM-DFM),the fracture flow module in the COMSOL Multiphysics simulator is powerful in definition of fractures as the inner flow boundary existing in the porous media.Thus it is selected to simulate the fluid flow and heat transfer in the geothermal-developed fractured granite of Sanguliu area located at Liaodong Peninsula,Eastern China.The natural faults/fractures based on field investigation combined with the discrete fracture network(DFN)generated by the MATLAB are used to represent the two-dimensional geological model.Numerical results show that early thermal breakthrough occurs at the production well caused by quick flow of cold water along the highly connected fractures.Suitable hydraulic fracturing treatments with proper injection rates,locations,etc.can efficiently hinder the thermal breakthrough time in the natural fracture system.Large well spacing helps the long-term operation of geothermal production,but it is highly dependent on the geometrical morphology of the fracture network.The enhancement of reservoir properties at the near-well regions can also increase the geothermal production efficiency.The results in this study can provide references to achieve a sustainable geothermal exploitation in fractured granitic geothermal reservoirs or hot dry rocks at depth.
基金This researchwas funded by“CNPC Science and technology project:Fine evaluation and prediction technology for complex reservoirs in overseas natural gas reservoirs,grant number 2018D-4305”.
文摘Characterizing natural fractures has a decisive effect on production forecasts in fractured oil and gas reservoirs.Discrete Fracture Networks(DFN)constitutes the main modeling framework for fractured geosystems.However,myriads of uncertainties are enclosed prior modeling representative stochastic or deterministic DFN ensembles.This paper presents a novel methodology for DFN calibration and an efficacious field application,which incorporatesWell-testing interpretation,Embedded Discrete Fracture Model(EDFM)framework,and numerical reservoir simulation.The proposed workflow starts with the DFN generation from seismic data,imaging logging data and core data.After multiple DFNs are modeled,well-test analysis is employed to calibrate the intrinsic properties of fractures at different locations.Then,these fracture networks are characterized dynamically by EDFM,which promotes capturing the optimal fracture model quickly screened.Finally,pressure and production history match are reached for the DFN realization that honors the optimal fracture model.
基金the National Natural Science Foundation of China(Grant Nos.41227901,41502294&41330643)the Strategic Priority Research Program of the Chinese Academy of Sciences(Grants Nos.XDB10030000,XDB10030300&XDB10050400)
文摘In this paper, the effect of pre-existing discrete fracture network(DFN) connectivity on hydraulic fracturing is numerically investigated in a rock mass subjected to in-situ stress. The simulation results show that DFN connectivity has a significant influence on the hydraulic fracture(HF) & DFN interaction and hydraulic fracturing effectiveness, which can be characterized by the total interaction area, stimulated DFN length, stimulated HF length, leak-off ratio, and stimulated total length. In addition, even at the same fluid injection rate, simulation models exhibit different responses that are strongly affected by the DFN connectivity. At a low injection rate, total interaction area decreases with increasing DFN connectivity; at a high injection rate, total interaction area increases with the increase of DFN connectivity. However, for any injection rate, the stimulated DFN length increases and stimulated HF length decreases with the increase of connectivity. Generally, this work shows that the DFN connectivity plays a crucial role in the interaction between hydraulic fractures, the pre-existing natural fractures and hydraulic fracturing effectiveness; in return, these three factors affect treating pressure, created microseismicity and corresponding stimulated volume. This work strongly relates to the production technology and the evaluation of hydraulic fracturing effectiveness. It is helpful for the optimization of hydraulic fracturing simulations in naturally fractured formations.
基金the Natural Sciences and Engineering Research Council of Canadathe Ministry of National Defensethe RMC Green Team for providing the funding and the resources
文摘Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-existing joints in the damage evolution around the underground opening is of critical importance as they govern the fracturing mechanisms and influence the brittle responses of these hard rock masses under highly anisotropic in situ stresses.In this study,the main focus is the impact of joint network geometry,joint strength and applied field stresses on the rock mass behaviours and the evolution of excavation induced damage due to the loss of confinement as a tunnel face advances.Analysis of such a phenomenon was conducted using the finite-discrete element method (FDEM).The numerical model is initially calibrated in order to match the behaviour of the fracture-free,massive Lac du Bonnet granite during the excavation of the Underground Research Laboratory (URL) Test Tunnel,Canada.The influence of the pre-existing joints on the rock mass response during excavation is investigated by integrating discrete fracture networks (DFNs) of various characteristics into the numerical models under varying in situ stresses.The numerical results obtained highlight the significance of the pre-existing joints on the reduction of in situ rock mass strength and its capacity for extension with both factors controlling the brittle response of the material.Furthermore,the impact of spatial distribution of natural joints on the stability of an underground excavation is discussed,as well as the potentially minor influence of joint strength on the stress induced damage within joint systems of a non-persistent nature under specific conditions.Additionally,the in situ stress-joint network interaction is examined,revealing the complex fracturing mechanisms that may lead to uncontrolled fracture propagation that compromises the overall stability of an underground excavation.
基金the Nuclear Waste Management Organization (NWMO) of Canadathe National Science and Engineering Research Council (NSERC)+1 种基金the Canadian Ministry of National Defence (DND)the RMC Green Team for funding this research
文摘The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks（well interlocked undisturbed rock mass with blocks formed by three or less intersecting joints）. A synthetic rock mass modelling（SRM） approach is employed to determine the mechanical properties of the rock mass. In this approach, the intact body of rock is represented by the discrete element method（DEM）-Voronoi grains with the ability of simulating the initiation and propagation of microcracks within the intact part of the model. The geometry of the preexisting joints is generated by employing discrete fracture network（DFN） modelling based on field joint data collected from the Brockville Tunnel using LiDAR scanning. The geometrical characteristics of the simulated joints at a representative sample size are first validated against the field data, and then used to measure the rock quality designation（RQD）, joint spacing, areal fracture intensity（P21）, and block volumes. These geometrical quantities are used to quantitatively determine a representative range of the geological strength index（GSI）. The results show that estimating the GSI using the RQD tends to make a closer estimate of the degree of blockiness that leads to GSI values corresponding to those obtained from direct visual observations of the rock mass conditions in the field. The use of joint spacing and block volume in order to quantify the GSI value range for the studied rock mass suggests a lower range compared to that evaluated in situ. Based on numerical modelling results and laboratory data of rock testing reported in the literature, a semi-empirical equation is proposed that relates the rock mass modulus to confinement as a function of the areal fracture intensity and joint stiffness.
文摘Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.
文摘Rockbursting in deep tunnelling is a complex phenomenon posing significant challenges both at the design and construction stages of an underground excavation within hard rock masses and under high in situ stresses. While local experience, field monitoring, and informed data-rich analysis are some of the tools commonly used to manage the hazards and the associated risks, advanced numerical techniques based on discontinuum modelling have also shown potential in assisting in the assessment of rockbursting. In this study, the hybrid finite-discrete element method(FDEM) is employed to investigate the failure and fracturing processes, and the mechanisms of energy storage and rapid release resulting in bursting, as well as to assess its utility as part of the design process of underground excavations.Following the calibration of the numerical model to simulate a deep excavation in a hard, massive rock mass, discrete fracture network(DFN) geometries are integrated into the model in order to examine the impact of rock structure on rockbursting under high in situ stresses. The obtained analysis results not only highlight the importance of explicitly simulating pre-existing joints within the model, as they affect the mobilised failure mechanisms and the intensity of strain bursting phenomena, but also show how the employed joint network geometry, the field stress conditions, and their interaction influence the extent and depth of the excavation induced damage. Furthermore, a rigorous analysis of the mass and velocity of the ejected rock blocks and comparison of the obtained data with well-established semi-empirical approaches demonstrate the potential of the method to provide realistic estimates of the kinetic energy released during bursting for determining the energy support demand.
