Estimation of the anisotropy of hydraulic conductivity through 3D fracture networks using the directional geological entropy

With an extension of the geological entropy concept in porous media,the approach called directional entrogram is applied to link hydraulic behavior to the anisotropy of the 3D fracture networks.A metric called directional entropic scale is used to measure the anisotropy of spatial order in different directions.Compared with the traditional connectivity indexes based on the statistics of fracture geometry,the directional entropic scale is capable to quantify the anisotropy of connectivity and hydraulic conductivity in heterogeneous 3D fracture networks.According to the numerical analysis of directional entrogram and fluid flow in a number of the 3D fracture networks,the hydraulic conductivities and entropic scales in different directions both increase with spatial order(i.e.,trace length decreasing and spacing increasing)and are independent of the dip angle.As a result,the nonlinear correlation between the hydraulic conductivities and entropic scales from different directions can be unified as quadratic polynomial function,which can shed light on the anisotropic effect of spatial order and global entropy on the heterogeneous hydraulic behaviors.

Propagation and sealing efficiency of chemical grouting in a two-dimensional fracture network with flowing water

In this study, an orthogonal array experiment is conducted by using a transparent fracture network replica. Image processing and theoretical analysis are performed to investigate the model sealing efficiency(SE), factors influencing SE, and the effect of flowing water on propagation. The results show that grout propagation can be classified into three patterns in the fracture network: sealing off, partial sealing,and major erosion. The factors controlling the SE in a descending order of the amount of influence are the initial water flow speed, fracture aperture, grout take, and gel time. An optimal value for the combination of the gel time and grout take(artificial factors) can result in a good SE. The grouting and seepage pressures are measured, and the results reveal that their variations can indicate the SE to some extent. The SE is good when the seepage pressure at each point increases overall;the frequent fluctuations in the seepage pressure indicate a moderately poor SE, and an overall decline in the seepage pressure indicates a major erosion type. The deflection effect of grouting shows an approximately elliptical propagation with the long axis expanding along the wider fracture opening, demonstrating further application in grouting design.

Development of an improved three-dimensional rough discrete fracture network model:Method and application

Structure plane is one of the important factors affecting the stability and failure mode of rock mass engineering.Rock mass structure characterization is the basic work of rock mechanics research and the important content of numerical simulation.A new 3-dimensional rough discrete fracture network(RDFN3D)model and its modeling method based on the Weierstrass-Mandelbrot(W-M)function were presented in this paper.The RDFN3D model,which improves and unifies the modelling methods for the complex structural planes,has been realized.The influence of fractal dimension,amplitude,and surface precision on the modeling parameters of RDFN3D was discussed.The reasonable W-M parameters suitable for the roughness coefficient of JRC were proposed,and the relationship between the mathematical model and the joint characterization was established.The RDFN3D together with the smooth 3-dimensional discrete fracture network(DFN3D)models were successfully exported to the drawing exchange format,which will provide a wide application in numerous numerical simulation codes including both the continuous and discontinuous methods.The numerical models were discussed using the COMSOL Multiphysics code and the 3-dimensional particle flow code,respectively.The reliability of the RDFN3D model was preliminarily discussed and analyzed.The roughness and spatial connectivity of the fracture networks have a dominant effect on the fluid flow patterns.The research results can provide a new geological model and analysis model for numerical simulation and engineering analysis of jointed rock mass.

Earthquake-induced fracture displacements and transmissivity changes in a 3D fracture network of crystalline rock for spent nuclear fuel disposal

During the long service period of a nuclear waste repository in crystalline rock,large earthquake(s)may occur nearby the repository site and coseismically alter the local stress field around pre-existing fractures within the geological formation.The resulting fracture normal/shear displacements may lead to fracture opening and further promote the transport of leaked radionuclides into the groundwater system.Thus,it is of central importance to analyze the consequences of potential future earthquake(s)on the hydrogeological properties of a repository site for spent nuclear fuel disposal.Based on the detailed site characterization data of the repository site at Forsmark,Sweden,we conduct a three-dimensional(3D)seismo-hydro-mechanical simulation using the 3Dimensional Distinct Element Code(3DEC).We explicitly represent a primary seismogenic fault zone and its surrounding secondary fracture network associated with a power-law size scaling and a Fisher orientation distribution.An earthquake with a magnitude of M_(w)=5.6 caused by the reactivation of the primary fault zone is modeled by simulating its transient rupture propagating radially outwards from a predefined hypocenter at a specified rupture speed,with the faulting dynamics controlled by a strength weakening law.We model the coseismic response of the off-fault fracture network subject to both static and dynamic triggering effects.We further diagnose the distribution of fracture hydro-mechanical properties(e.g.mechanical/hydraulic aperture,hydraulic transmissivity)before and after the earthquake in order to quantify earthquakeinduced hydraulic changes in the fracture network.It is found that earthquake-induced fracture transmissivity changes tend to follow a power-law decay with the distance to the earthquake fault.Our simulation results and insights obtained have important implications for the long-term performance assessment of nuclear waste repositories in fractured crystalline rocks.

Well interference evaluation considering complex fracture networks through pressure and rate transient analysis in unconventional reservoirs

Severe well interference through complex fracture networks(CFNs)can be observed among multi-well pads in low permeability reservoirs.The well interference analysis between multi-fractured horizontal wells(MFHWs)is vitally important for reservoir effective development.Well interference has been historically investigated by pressure transient analysis,while it has shown that rate transient analysis has great potential in well interference diagnosis.However,the impact of complex fracture networks(CFNs)on rate transient behavior of parent well and child well in unconventional reservoirs is still not clear.To further investigate,this paper develops an integrated approach combining pressure and rate transient analysis for well interference diagnosis considering CFNs.To perform multi-well simulation considering CFNs,non-intrusive embedded discrete fracture model approach was applied for coupling fracture with reservoir models.The impact of CFN including natural fractures and frac-hits on pressure and rate transient behavior in multi-well system was investigated.On a logelog plot,interference flow and compound linear flow are two new flow regimes caused by nearby producers.When both NFs and frac-hits are present in the reservoir,frac-hits have a greater impact on well#1 which contains frac-hits,and NFs have greater impact on well#3 which does not have frac-hits.For all well producing circumstances,it might be challenging to see divergence during pseudosteady state flow brought on by frac-hits on the logelog plot.Besides,when NFs occur,reservoir depletion becomes noticeable in comparison to frac-hits in pressure distribution.Application of this integrated approach demonstrates that it works well to characterize the well interference among different multi-fractured horizontal wells in a well pad.Better reservoir evaluation can be acquired based on the new features observed in the novel model,demonstrating the practicability of the proposed approach.The findings of this study can help for better evaluating well interference degree in multi-well systems combing PTA and RTA,which can reduce the uncertainty and improve the accuracy of the well interference analysis based on both field pressure and rate data.

Pump-stopping pressure drop model considering transient leak-off of fracture network

By introducing the coupling flow expressions of main fracture-matrix, secondary fracture-matrix and main fracture-secondary fracture into the traditional main fracture material balance equation, the “main fracture-secondary fracture-matrix” leak-off coupling flow model is established. The pressure-dependent fracture width equation and the wellbore injection volume equation are coupled to solve the pressure-rate continuity problem. The simulation and calculation of the bottomhole pressure drop and fracture network closure after the pump stopping in slickwater volumetric fracturing treatment are realized. The research results show that the log-log curve of pump-stopping bottomhole pressure drop derivative presents five characteristic slope segments, reflecting four dominant stages, i.e. inter-fracture crossflow, fracture network leak-off, fracture network closure and residual leak-off, after pump shutdown. At the initial time of pump shutdown for volumetric fracturing treatment of horizontal well, the crossflow between main and secondary fractures is obvious, and then the leak-off becomes dominant. The leak-off of main and secondary fractures shows a non-uniform decreasing trend. Specifically, the leak-off of main fractures is slow, while that of secondary fractures is fast;the fracture network as a whole presents the leak-off law of fast first, then slow, until close to zero. The influence of fracture network conductivity on the shape of pressure decline curve is relatively weaker than that of fracture network size. The fracture network conductivity is positively correlated with leak-off volume and fracture closure. The secondary fracture size is positively correlated with leakoff volume and closure of the secondary fracture, but negatively correlated with closure of the main fracture. Field data validation proves that the proposed model and simulation results can effectively reflect the closure characteristics of the fracture network, and the interpretation results are reliable and can reflect the non-uniform stimulation performance of each fracturing stage of an actual horizontal well.

Horizontal Well Interference Performance and Water Injection Huff and Puff Effect on Well Groups with Complex Fracture Networks:A Numerical Study

Well interference has become a common phenomenon with the increasing scale of horizontal well fracturing.Recent studies on well interference in horizontal wells do not properly reflect the physical model of the postfracturing well groups and the realistic fracturing process of infill wells.Establishing the correspondence between well interference causative factors and manifestations is of great significance for infill well deployment and secondary oil recovery.In this work,we develop a numerical model that considers low velocity non-Darcy seepage inshale reservoirs to study the inter-well interferencephenomenon that occurs in theSantanghufield,andconstruct an explicit hydraulic fracture and complex natural fracture network model with an embedded discrete fracture model,focusing on the effect of fracture network morphology on well interactions.The model also considers a multi-segment wellbore model to accommodate the effect of inter-well crossflow on wellbore tubular flow.The changes in formation pressure and water saturation during fracturing are performed by controlling the injection pressure and water injection rate.The result shows that the shape of the fracture network generated by the infill well with the old well determines the subsequent fluid and oil-increasing performance of the disturbed well.The synergistic production or competitive relationship formed by fractures with different connectivity between the two wells determines the positive and negative effects of the interference.The paper also investigates the adaptation study of water injection huff and puff schemes for well groups with different connectivity,and demonstrated a potential yield increase of up to 10.85%under adaptation injection.This method of identifying well interference based on the production dynamics of affected wells and the subsequent corresponding water injection method provides valuable references for the selection of secondary oil recovery measures.

Quantitative characterization of irregular microfracture network and its effect on the permeability of porous media

Based on the comprehensive understanding on microfractures and matrix pores in reservoir rocks,numerical algorithms are used to construct fractured porous media and fracture-pore media models.Connectivity coefficient and strike factor are introduced into the models to quantitatively characterize the connectivity and strike of fracture network,respectively.The influences of fracture aperture,fracture strike and fracture connectivity on the permeability of porous media are studied by using multi-relaxation-time lattice Boltzmann model to simulate fluid flow in them.The greater the strike factor and the smaller the tortuosity of the fractured porous media,the greater the permeability of the fractured porous media.The greater the connectivity coefficient of the fracture network is,the greater the permeability of the fracture-pore media is,and the more likely dominant channel effect occurs.The fracture network connectivity has stronger influence on seepage ability of fracture-pore media than fracture aperture and fracture strike.The tortuosity and strike factor of fracture network in fractured porous media are in polynomial relation,while the permeability and fracture network connectivity coefficient of the fracture-pore media meet an exponential relation.

3D reconstruction method and connectivity rules of fracture networks generated under different mining layouts

In current research,a series of triaxial tests,which were employed to simulate three typical mining layouts(i.e.,top-coal caving,non-pillar mining and protected coal seam mining),were conducted on coal by using MTS815 Flex Test GT rock mechanics test system,and the fracture networks in the broken coal samples were qualitatively and quantitatively investigated by employing CT scanning and 3D reconstruction techniques.This work aimed at providing a detail description on the micro-structure and fractureconnectivity characteristics of rupture coal samples under different mining layouts.The results show that:(i)for protected coal seam mining layout,the coal specimens failure is in a compression-shear manner and oppositely,(ii)the tension-shear failure phenomenon is observed for top-coal caving and non-pillar mining layouts.By investigating the connectivity features of the generated fractures in the direction ofσ_1,under different mining layouts,it is found that the connectivity level of the fractures of the samples corresponding to non-pillar mining layout was the highest.

Fractured reservoir modeling by discrete fracture network and seismic modeling in the Tarim Basin,China

Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.

Characterizing the influence of stress-induced microcracks on the laboratory strength and fracture development in brittle rocks using a finite-discrete element method-micro discrete fracture network FDEM-μDFN approach

Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.

An Elevated Perspective: Dyke-Related Fracture Networks Analysed with Uav Photogrammetry

An abundance of data from seismic and geodetic monitoring has provided new insight into dyke propagation and emplacement mechanisms.These studies show that faulting and fracturing is part of the magma

Simulation of the Pressure-Sensitive Seepage Fracture Network in Oil Reservoirs with Multi-Group Fractures

Stress sensitivity is a very important index to understand the seepage characteristics of a reservoir.In this study,dedicated experiments and theoretical arguments based on the visualization of porous media are used to assess the effects of the fracture angle,spacing,and relevant elastic parameters on the principal value of the permeability tensor.The fracture apertures at different angles show different change rates,which influence the relative permeability for different sets of fractures.Furthermore,under the same pressure condition,the fractures with different angles show different degrees of deformation so that the principal value direction of permeability rotates.This phenomenon leads to a variation in the water seepage direction in typical water-injection applications,thereby hindering the expected exploitation effect of the original well network.Overall,the research findings in this paper can be used as guidance to improve the effectiveness of water injection exploitation in the oil field industry.

Nonlinear fluid flow through three-dimensional rough fracture networks:Insights from 3D-printing,CT-scanning,and high-resolution numerical simulations

Nonlinear flow behavior of fluids through three-dimensional(3D)discrete fracture networks(DFNs)considering effects of fracture number,surface roughness and fracture aperture was experimentally and numerically investigated.Three physical models of DFNs were 3D-printed and then computed tomography(CT)-scanned to obtain the specific geometry of fractures.The validity of numerically simulating the fluid flow through DFNs was verified via comparison with flow tests on the 3D-printed models.A parametric study was then implemented to establish quantitative relations between the coefficients/parameters in Forchheimer’s law and geometrical parameters.The results showed that the 3D-printing technique can well reproduce the geometry of single fractures with less precision when preparing complex fracture networks,numerical modeling precision of which can be improved via CT-scanning as evidenced by the well fitted results between fluid flow tests and numerical simulations using CT-scanned digital models.Streamlines in DFNs become increasingly tortuous as the fracture number and roughness increase,resulting in stronger inertial effects and greater curvatures of hydraulic pressure-low rate relations,which can be well characterized by the Forchheimer’s law.The critical hydraulic gradient for the onset of nonlinear flow decreases with the increasing aperture,fracture number and roughness,following a power function.The increases in fracture aperture and number provide more paths for fluid flow,increasing both the viscous and inertial permeabilities.The value of the inertial permeability is approximately four orders of magnitude greater than the viscous permeability,following a power function with an exponent a of 3,and a proportional coefficient b mathematically correlated with the geometrical parameters.

Performance-based fractal fracture model for complex fracture network simulation

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.

Fracture Network Volume Fracturing Technology in High-temperature Hard Formation of Hot Dry Rock

It is more difficult for a hot dry rock to form a fracture network system than shale due to its special lithology, physical and mechanical properties under high temperature. The essential characteristics, rock mechanics and in-situ stress characteristics of a hot rock mass have been systematically studied by means of laboratory tests and true tri-axial physical simulation. The fracture initiation and propagation characteristics under different geological and engineering conditions are physically simulated, and the main controlling factors for the formation of a complex fracture network are revealed. The technology of low displacement for enhancing thermal cracking, gel fluid for expanding fracture and variable displacement cyclic injection for increasing a fracture network has been applied in the field, and good results have been achieved. Microseismic monitoring results demonstrate that complex fractures were formed in the field test, and the stimulation volume for heat exchanging reaches more than 3 million cubic meters. The research results play an important role in the stimulation technology of an enhanced geothermal system(EGS) and realize a breakthrough for power generation.

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.

A smart productivity evaluation method for shale gas wells based on 3D fractal fracture network model

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.

An optimization model for conductivity of hydraulic fracture networks in the Longmaxi shale,Sichuan basin,Southwest China

Shale gas is an important unconventional resource.The economic recovery of shale gas is only possible when a fracture network with sufficient conductivity is created by hydraulic fracturing,that,if effectively propped,connects fracturing fractures and natural fractures.Focusing on the Longmaxi shale in the Sichuan Basin,Southwest China,we built an optimization model for conductivity of multi-grade fractures based on equivalent seepage theory.We then experimentally analyzed the conductivity of self-propped and sand-propped fractures,and optimized the propping patterns of multi-grade hydraulic fractures in shale gas reservoirs.We concluded that the propping effectiveness of fracture networks could be improved by using low concentrations of small-sized sands and by focusing on creating a large number of self-propped fractures.By applying this understanding to the optimization of fracturing designs for the Longmaxi shale,we successfully created networks of well-propped fractures.

Study on the Impact of Massive Refracturing on the Fracture Network in Tight Oil Reservoir Horizontal Wells

Class III tight oil reservoirs have low porosity and permeability,which are often responsible for low production rates and limited recovery.Extensive repeated fracturing is a well-known technique to fix some of these issues.With such methods,existing fractures are refractured,and/or new fractures are created to facilitate communication with natural fractures.This study explored how different refracturing methods affect horizontal well fracture networks,with a special focus on morphology and related fluid flow changes.In particular,the study relied on the unconventional fracture model(UFM).The evolution of fracture morphology and flow field after the initial fracturing were analyzed accordingly.The simulation results indicated that increased formation energy and reduced reservoir stress differences can promote fracture expansion.It was shown that the length of the fracture network,the width of the fracture network,and the complexity of the fracture can be improved,the oil drainage area can be increased,the distance of oil and gas seepage can be reduced,and the production of a single well can be significantly increased.