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.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.

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.

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.

Graph theoretical analysis of limestone fracture network damage patterns based on uniaxial compression test

The topological attributes of fracture networks in limestone,subject to intense hydrodynamics and intricate geological discontinuities,substantially influence the mechanical and hydraulic characteristics of the rock mass.The dynamical evolution of fracture networks under stress is crucial for unveiling the interaction patterns among fractures.However,existing models are undirected graphs focused on stationary topology,which need optimization to depict fractures'dynamic development and rupture process.To compensate for the time and destruction terms,we propose the damage network model,which defines the physical interpretation of fractures through the ternary motif.We focus primarily on the evolution of node types,topological attributes,and motifs of the fracture network in limestone under uniaxial stress.Observations expose the varying behavior of the nodes'self-dynamics and neighbors'adjacent dynamics in the fracture network.This approach elucidates the impact of micro-crack behaviors on large brittle shear fractures from a topological perspective and further subdivides the progressive failure stage into four distinct phases(isolated crack growth phase,crack splay phase,damage coalescence phase,and mechanical failure phase)based on the significance profile of the motif.Regression analysis reveals a positive linear and negative power correlation between fracture network density and branch number to the rock damage resistance,respectively.The damage network model introduces a novel methodology for depicting the interaction of two-dimensional(2D)projected fractures,considering the dynamic spatiotemporal development characteristics and fracture geometric variation.It helps dynamically characterize properties such as connectivity,permeability,and damage factors while comprehensively assessing damage in rock mass fracture networks.

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.

Proppant transport in rough fracture networks using supercritical CO_(2)

Proppant transport within fractures is one of the most critical tasks in oil,gas and geothermal reservoir stimulation,as it largely determines the ultimate performance of the operating well.Proppant transport in rough fracture networks is still a relatively new area of research and the associated transport mechanisms are still unclear.In this study,representative parameters of rough fracture surfaces formed by supercritical CO_(2) fracturing were used to generate a rough fracture network model based on a spectral synthesis method.Computational fluid dynamics(CFD)coupled with the discrete element method(DEM)was used to study proppant transport in this rough fracture network.To reveal the turning transport mechanism of proppants into branching fractures at the intersections of rough fracture networks,a comparison was made with the behavior within smooth fracture networks,and the effect of key pumping parameters on the proppant placement in a secondary fracture was analyzed.The results show that the transport behavior of proppant in rough fracture networks is very different from that of the one in the smooth fracture networks.The turning transport mechanisms of proppant into secondary fractures in rough fracture networks are gravity-driven sliding,high velocity fluid suspension,and fracture structure induction.Under the same injection conditions,supercritical CO_(2)with high flow Reynolds number still has a weaker ability to transport proppant into secondary fractures than water.Thickening of the supercritical CO_(2)needs to be increased beyond a certain value to have a significant effect on proppant carrying,and under the temperature and pressure conditions of this paper,it needs to be increased more than 20 times(about 0.94 m Pa s).Increasing the injection velocity and decreasing the proppant concentration facilitates the entry of proppant into the branching fractures,which in turn results in a larger stimulated reservoir volume.The results help to understand the proppant transport and placement process in rough fracture networks formed by reservoir stimulation,and provide a theoretical reference for the optimization of proppant pumping parameters in hydraulic fracturing.

Fracture network types revealed by well test curves for shale reservoirs in the Sichuan Basin,China

Pressure buildup testing can be used to analyze fracture network characteristics and conduct quantitative interpretation of relevant parameters for shale gas wells,thus providing bases for assessing the well productivity and formulating proper development strategies.This study establishes a new well test interpretation model for fractured horizontal wells based on seepage mechanisms of shale reservoirs and proposes a method for identifying fracturing patterns based on the characteristic slopes of pressure buildup curves and curve combination patterns.The pressure buildup curve patterns are identified to represent three types of shale reservoirs in the Sichuan Basin,namely the moderately deep shale reservoirs with high pressure,deep shale reservoirs with ultra-high pressure,and moderately deep shale reservoirs with normal pressure.Based on this,the relationship between the typical pressure buildup curve patterns and the fracture network types are put forward.Fracturing effects of three types of shale gas reservoir are compared and analyzed.The results show that typical flow patterns of shale reservoirs include bilinear flow in primary and secondary fractures,linear flow in secondary fractures,bilinear flow in secondary fractures and matrix,and linear flow in matrix.The fracture network characteristics can be determined using the characteristic slopes of pressure buildup curves and curve combinations.The linear flow in early secondary fractures is increasingly distinct with an increase in primary fracture conductivity.Moreover,the bilinear flow in secondary fractures and matrix and the subsequent linear flow in the matrix occur as the propping and density of secondary fractures increase.The increase in the burial depth,in-situ stress,and stress difference corresponds to a decrease in the propping of primary fractures that expand along different directions in the shale gas wells in the Sichuan Basin.Four pressure buildup curve patterns exist in the Sichuan Basin and its periphery.The pattern of pressure buildup curves of shale reservoirs in the Yongchuan area can be described as 1/2/→1/4,indicating limited stimulated reservoir volume,poorly propped secondary fractures,and the forming of primary fractures that extend only to certain directions.The pressure buildup curves of shale reservoirs in the main block of the Fuling area show a pattern of 1/4/→1/2 or 1/2,indicating greater stimulated reservoir volume,well propped secondary fractures,and the forming of complex fracture networks.The pattern of pressure buildup curves of shale reservoirs in the Pingqiao area is 1/2/→1/4→/1/2,indicating a fracturing effect somewhere between that of the Fuling and Yongchuan areas.For reservoirs with normal pressure,it is difficult to determine fracture network characteristics from pressure buildup curves due to insufficient formation energy and limited liquid drainage.

Hydrocarbon gas huff-n-puff optimization of multiple horizontal wells with complex fracture networks in the M unconventional reservoir

The oil production of the multi-fractured horizontal wells(MFHWs) declines quickly in unconventional oil reservoirs due to the fast depletion of natural energy. Gas injection has been acknowledged as an effective method to improve oil recovery factor from unconventional oil reservoirs. Hydrocarbon gas huff-n-puff becomes preferable when the CO_(2) source is limited. However, the impact of complex fracture networks and well interference on the EOR performance of multiple MFHWs is still unclear. The optimal gas huff-n-puff parameters are significant for enhancing oil recovery. This work aims to optimize the hydrocarbon gas injection and production parameters for multiple MFHWs with complex fracture networks in unconventional oil reservoirs. Firstly, the numerical model based on unstructured grids is developed to characterize the complex fracture networks and capture the dynamic fracture features.Secondly, the PVT phase behavior simulation was carried out to provide the fluid model for numerical simulation. Thirdly, the optimal parameters for hydrocarbon gas huff-n-puff were obtained. Finally, the dominant factors of hydrocarbon gas huff-n-puff under complex fracture networks are obtained by fuzzy mathematical method. Results reveal that the current pressure of hydrocarbon gas injection can achieve miscible displacement. The optimal injection and production parameters are obtained by single-factor analysis to analyze the effect of individual parameter. Gas injection time is the dominant factor of hydrocarbon gas huff-n-puff in unconventional oil reservoirs with complex fracture networks. This work can offer engineers guidance for hydrocarbon gas huff-n-puff of multiple MFHWs considering the complex fracture networks.

Prediction of three-dimensional ocean temperature in the South China Sea based on time series gridded data and a dynamic spatiotemporal graph neural network

Ocean temperature is an important physical variable in marine ecosystems,and ocean temperature prediction is an important research objective in ocean-related fields.Currently,one of the commonly used methods for ocean temperature prediction is based on data-driven,but research on this method is mostly limited to the sea surface,with few studies on the prediction of internal ocean temperature.Existing graph neural network-based methods usually use predefined graphs or learned static graphs,which cannot capture the dynamic associations among data.In this study,we propose a novel dynamic spatiotemporal graph neural network(DSTGN)to predict threedimensional ocean temperature(3D-OT),which combines static graph learning and dynamic graph learning to automatically mine two unknown dependencies between sequences based on the original 3D-OT data without prior knowledge.Temporal and spatial dependencies in the time series were then captured using temporal and graph convolutions.We also integrated dynamic graph learning,static graph learning,graph convolution,and temporal convolution into an end-to-end framework for 3D-OT prediction using time-series grid data.In this study,we conducted prediction experiments using high-resolution 3D-OT from the Copernicus global ocean physical reanalysis,with data covering the vertical variation of temperature from the sea surface to 1000 m below the sea surface.We compared five mainstream models that are commonly used for ocean temperature prediction,and the results showed that the method achieved the best prediction results at all prediction scales.

Numerical simulation of three-dimensional fracturing fracture propagation in radial wells

A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.

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.

3D DEM simulation of hard rock fracture in deep tunnel excavation induced by changes in principal stress magnitude and orientation

To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.

Borehole stability in naturally fractured rocks with drilling mud intrusion and associated fracture strength weakening:A coupled DFN-DEM approach

Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.

Effect of navigation endoscopy combined with threedimensional printing technology in the treatment of orbital blowout fractures

●AIM:To explore the combined application of surgical navigation nasal endoscopy(NNE)and three-dimensional printing technology(3DPT)for the adjunctive treatment of orbital blowout fractures(OBF).●METHODS:Retrospective analysis was conducted on the data of patients with OBF who underwent surgical treatment at the Affiliated Eye Hospital of Nanchang University between July 2012 and November 2022.The control group consisted of patients who received traditional surgical treatment(n=43),while the new surgical group(n=52)consisted of patients who received NNE with 3DPT.The difference in therapeutic effects between the two groups was evaluated by comparing the duration of the operation,best corrected visual acuity(BCVA),enophthalmos difference,recovery rate of eye movement disorder,recovery rate of diplopia,and incidence of postoperative complications.●RESULTS:The study included 95 cases(95 eyes),with 63 men and 32 women.The patients’age ranged from 5 to 67y(35.21±15.75y).The new surgical group and the control group exhibited no statistically significant differences in the duration of the operation,BCVA and enophthalmos difference.The recovery rates of diplopia in the new surgical group were significantly higher than those in the control group at 1mo[OR=0.03,95%CI(0.01–0.15),P<0.0000]and 3mo[OR=0.11,95%CI(0.03–0.36),P<0.0000]postoperation.Additionally,the recovery rates of eye movement disorders at 1 and 3mo after surgery were OR=0.08,95%CI(0.03–0.24),P<0.0000;and OR=0.01,95%CI(0.00–0.18),P<0.0000.The incidence of postoperative complications was lower in the new surgical group compared to the control group[OR=4.86,95%CI(0.95–24.78),P<0.05].●CONCLUSION:The combination of NNE and 3DPT can shorten the recovery time of diplopia and eye movement disorder in patients with OBF.

Supramolecular polymer-based gel fracturing fluid with a double network applied in ultra-deep hydraulic fracturing

A gel based on polyacrylamide,exhibiting delayed crosslinking characteristics,emerges as the preferred solution for mitigating degradation under conditions of high temperature and extended shear in ultralong wellbores.High viscosity/viscoelasticity of the fracturing fluid was required to maintain excellent proppant suspension properties before gelling.Taking into account both the cost and the potential damage to reservoirs,polymers with lower concentrations and molecular weights are generally preferred.In this work,the supramolecular action was integrated into the polymer,resulting in significant increases in the viscosity and viscoelasticity of the synthesized supramolecular polymer system.The double network gel,which is formed by the combination of the supramolecular polymer system and a small quantity of Zr-crosslinker,effectively resists temperature while minimizing permeability damage to the reservoir.The results indicate that the supramolecular polymer system with a molecular weight of(268—380)×10^(4)g/mol can achieve the same viscosity and viscoelasticity at 0.4 wt%due to the supramolecular interaction between polymers,compared to the 0.6 wt%traditional polymer(hydrolyzed polyacrylamide,molecular weight of 1078×10^(4)g/mol).The supramolecular polymer system possessed excellent proppant suspension properties with a 0.55 cm/min sedimentation rate at 0.4 wt%,whereas the0.6 wt%traditional polymer had a rate of 0.57 cm/min.In comparison to the traditional gel with a Zrcrosslinker concentration of 0.6 wt%and an elastic modulus of 7.77 Pa,the double network gel with a higher elastic modulus(9.00 Pa)could be formed only at 0.1 wt%Zr-crosslinker,which greatly reduced the amount of residue of the fluid after gel-breaking.The viscosity of the double network gel was66 m Pa s after 2 h shearing,whereas the traditional gel only reached 27 m Pa s.

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.

A three-dimensional feature extraction-based method for coal cleat characterization using X-ray μCT and its application to a Bowen Basin coal specimen

Cleats are the dominant micro-fracture network controlling the macro-mechanical behavior of coal.Improved understanding of the spatial characteristics of cleat networks is therefore important to the coal mining industry.Discrete fracture networks(DFNs)are increasingly used in engineering analyses to spatially model fractures at various scales.The reliability of coal DFNs largely depends on the confidence in the input cleat statistics.Estimates of these parameters can be made from image-based three-dimensional(3D)characterization of coal cleats using X-ray micro-computed tomography(m CT).One key step in this process,after cleat extraction,is the separation of individual cleats,without which the cleats are a connected network and statistics for different cleat sets cannot be measured.In this paper,a feature extraction-based image processing method is introduced to identify and separate distinct cleat groups from 3D X-ray m CT images.Kernels(filters)representing explicit cleat features of coal are built and cleat separation is successfully achieved by convolutional operations on 3D coal images.The new method is applied to a coal specimen with 80 mm in diameter and 100 mm in length acquired from an Anglo American Steelmaking Coal mine in the Bowen Basin,Queensland,Australia.It is demonstrated that the new method produces reliable cleat separation capable of defining individual cleats and preserving 3D topology after separation.Bedding-parallel fractures are also identified and separated,which has his-torically been challenging to delineate and rarely reported.A variety of cleat/fracture statistics is measured which not only can quantitatively characterize the cleat/fracture system but also can be used for DFN modeling.Finally,variability and heterogeneity with respect to the core axis are investigated.Significant heterogeneity is observed and suggests that the representative elementary volume(REV)of the cleat groups for engineering purposes may be a complex problem requiring careful consideration.

Research of the molecular mechanism of Hai Honghua medicinal liquor in the treatment of fracture based on network pharmacology combined with molecular docking

Background:Hai Honghua medicinal liquor(HHML),a famous hospital formula composed of 19 traditional Chinese medicines,has been successfully applied in treating soft tissue injury,fresh closed fracture,limb dysfunction after fracture healing,shoulder,neck and leg pain,knee joint pain and other clinical multiple diseases for 30 years in clinical.However,research on the material basis of HHML for the treatment of fracture healing-related disorders is still in a gap.Therefore,it is particularly important to explore the active ingredients,core targets and potential pharmacological mechanisms of HHML to promote fracture healing.Methods:We screened the core active components of each traditional Chinese medicine in formula and its action targets through the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform and Encyclopedia of Traditional Chinese Medicine database.The fracture related targets were retrieved from several different public databases,including GeneCards,Online Mendelian Inheritance in Man,DisGeNET and Therapeutic Target Database.Bioinformatics analysis to obtain key bioactive components,underlying targets and signaling pathways,containing the Venn diagram of the intersection with components and diseases gene targets,protein–protein interaction,as well as the Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes analysis,and finally molecular docking.Results:A total of 249 bioactive ingredients of HHML and 325 HHML-fracture-related targets were screened.The network analysis revealed that quercetin,luteolin,kaempferol,Licochalcone A,naringenin and 8-Isopentenyl-kaempferol may be potential candidate agents.Multiple targets are involved including TP53,MAPK3,STAT3,AKT1,MAPK1,HSP90AA1,ESR1 and PIK3CA may be closely linked targets.PI3K-AKT signaling pathway may play a significant role of HHML in treatment of fracture.What’s more,molecular docking suggested that 8-isopentenyl kaempferol,glycyrrhiza chalcone A,and naringenin bound to AKT1,PIK3CA,and ESR1,respectively,exhibiting lower energy and more stable characteristics.Conclusions:The findings indicate the potential active ingredients,target proteins and molecular mechanisms of HHML for the treatment of fractures to provide the exact idea for the next research on the mechanism of action of HHML formula for fracture treatment.

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.