Mechanisms of fracture propagation from multi-cluster using a phase field based HMD coupling model in fractured reservoir

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.

Numerical modeling of the effects of roughness on flow and eddy formation in 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.

Numerical modelling of flow and transport in rough fractures

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.

Application of Stochastic Fracture Network with Numerical Fluid Flow Simulations to Groundwater Flow Modeling in Fractured Rocks

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.

Height Detection and Analysis of Water Flowing Fractured Zone of Coal Face

Taking 91105 working face as the research object, the observation method of water flowing fractured zone and the layout of mining holes were determined by analyzing the field geological structure. 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.

Algorithmic approach to discrete fracture network flow modeling in consideration of realistic connections in large-scale fracture networks

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.

Rate decline analysis of multiple fractured horizontal well in shale reservoir with triple continuum

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.

Forecasting water disaster for a coal mine under the Xiaolangdi reservoir

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.

Direct shear tests on cemented paste backfill-rock wall and cemented paste backfill-backfill interfaces

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.

NON-LCL AND TRACER TEST FOR GROUNDWATER FLOW IN A SINGLE FRACTURE

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.

Optimum Calculation of Coal Pillars in Inclined Weathered Oxidation Zone

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.

Forecasting of destroyed height of overlying rock with the top coal caving based on ANN

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.

Transient electromagnetic detection method in water-sealed underground storage caverns

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.

Geothermal data analysis at the high-temperature hydrothermal area in Western Sichuan

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.