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 ...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.展开更多
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 directi...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.展开更多
Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the...Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the concept of large-scale stimulation by fracture network,balanced propagation and effective support of fracture network in fracturing design and developed the extreme massive hydraulic fracturing technique for deep coalbed methane(CBM)horizontal wells.This technique involves massive injection with high pumping rate+high-intensity proppant injection+perforation with equal apertures and limited flow+temporary plugging and diverting fractures+slick water with integrated variable viscosity+graded proppants with multiple sizes.The technique was applied in the pioneering test of a multi-stage fracturing horizontal well in deep CBM of Linxing Block,eastern margin of the Ordos Basin.The injection flow rate is 18 m^(3)/min,proppant intensity is 2.1 m^(3)/m,and fracturing fluid intensity is 16.5 m^(3)/m.After fracturing,a complex fracture network was formed,with an average fracture length of 205 m.The stimulated reservoir volume was 1987×10^(4)m^(3),and the peak gas production rate reached 6.0×10^(4)m^(3)/d,which achieved efficient development of deep CBM.展开更多
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 m...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.展开更多
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...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.展开更多
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 meth...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.展开更多
Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based ...Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based on engineering properties to simulate the gradual collapse of the roof during longwall top coal caving(LTCC).A numerical model is established using the material point method(MPM)and the strain-softening damage constitutive model according to the structure of the physical model.Numerical simulations are conducted to analyze the LTCC process under different hard roofs for ground hydraulic fracturing.The results show that ground hydraulic fracturing releases the energy and stress of the target stratum,resulting in a substantial lag in the fracturing of the overburden before collapse occurs in the hydraulic fracturing stratum.Ground hydraulic fracturing of a low hard roof reduces the lag effect of hydraulic fractures,dissipates the energy consumed by the fracture of the hard roof,and reduces the abutment stress.Therefore,it is advisable to prioritize the selection of the lower hard roof as the target stratum.展开更多
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 conne...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.展开更多
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 productiv...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.展开更多
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 mechani...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.展开更多
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 vit...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.展开更多
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), fact...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.展开更多
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-se...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.展开更多
This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation...This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.展开更多
In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution a...In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.展开更多
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 wit...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.展开更多
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 con...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.展开更多
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 postfr...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.展开更多
After volume fracturing of horizontal wells in shale gas reservoir, an extremely complex fracture system is formed. The space area of the fracture system is the reservoir reconstruction volume of shale gas reservoir. ...After volume fracturing of horizontal wells in shale gas reservoir, an extremely complex fracture system is formed. The space area of the fracture system is the reservoir reconstruction volume of shale gas reservoir. The geometric parameters such as crack length, crack width, crack height, and characteristic parameters such as crack permeability and fracture conductivity proposed for a single crack in conventional fracturing are insufficient to describe and characterize the complex network fracture system after volume fracturing. In this paper, the discrete fracture modeling method is used to establish the volume fracturing network fracture model of horizontal wells in shale gas reservoir by using the random modeling method within the determined reservoir space. The model is random and selective, and can fully provide different forms of volume fracturing fracture expansion, such as conventional fracture morphology, line network model and arbitrarily distributed network fractures. The research results provide a theoretical basis for the development plan and stimulation plan of shale gas reservoir, and have important reference value and significance for other unconventional gas reservoir fracturing.展开更多
In-situ conversion presents a promising technique for exploiting continental oil shale formations,characterized by highly fractured organic-rich rock.A 3D in-situ conversion model,which incorporates a discrete fractur...In-situ conversion presents a promising technique for exploiting continental oil shale formations,characterized by highly fractured organic-rich rock.A 3D in-situ conversion model,which incorporates a discrete fracture network,is developed using a self-developed thermal-flow-chemical(TFC)simulator.Analysis of the model elucidates the in-situ conversion process in three stages and defines the transformation of fluids into three distinct outcomes according to their end stages.The findings indicate that kerogen decomposition increases fluid pressure,activating fractures and subsequently enhancing permeability.A comprehensive analysis of activated fracture permeability and heating power reveals four distinct production modes,highlighting that increasing heating power correlates with higher cumulative fluid production.Activated fractures,with heightened permeability,facilitate the mobility of heavy oil toward production wells but hinder its cracking,thereby limiting light hydrocarbon production.Additionally,energy efficiency research demonstrates the feasibility of the in-situ conversion in terms of energy utilization,especially when considering the surplus energy from high-fluctuation energy sources such as wind and solar power to provide heating.展开更多
基金the China Research and Pilot Test on Key Technology of Efficient Production of Changqing Tight Oil(Grant No.2021DJ2202).
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.42077243,52209148,and 52079062).
文摘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.
基金Supported by the National Natural Science Foundation of China Project(52274014)Comprehensive Scientific Research Project of China National Offshore Oil Corporation(KJZH-2023-2303)。
文摘Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the concept of large-scale stimulation by fracture network,balanced propagation and effective support of fracture network in fracturing design and developed the extreme massive hydraulic fracturing technique for deep coalbed methane(CBM)horizontal wells.This technique involves massive injection with high pumping rate+high-intensity proppant injection+perforation with equal apertures and limited flow+temporary plugging and diverting fractures+slick water with integrated variable viscosity+graded proppants with multiple sizes.The technique was applied in the pioneering test of a multi-stage fracturing horizontal well in deep CBM of Linxing Block,eastern margin of the Ordos Basin.The injection flow rate is 18 m^(3)/min,proppant intensity is 2.1 m^(3)/m,and fracturing fluid intensity is 16.5 m^(3)/m.After fracturing,a complex fracture network was formed,with an average fracture length of 205 m.The stimulated reservoir volume was 1987×10^(4)m^(3),and the peak gas production rate reached 6.0×10^(4)m^(3)/d,which achieved efficient development of deep CBM.
基金supported by the National Natural Science Foundation of China(Grant No.52090081)the State Key Laboratory of Hydroscience and Engineering(Grant No.2022-KY-02).
文摘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.
基金the support from the National Key Research and Development Program of China(Grant No.2022YFE0137200)the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2024JC-YBQN-0381,2023JC-QN-0403)+2 种基金the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2022JC-37)the Innovation Capability Support Program of Shaanxi(Program No.2023-CX-TD31)the Funded by Open Foundation of Shaanxi Key Laboratory of Carbon Dioxide Sequestration and Enhanced Oil Recovery,and the Youth Innovation Team of Shaanxi Universities。
文摘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.
基金funded by the National Natural Science Foundation of China(No.51974268)Open Fund of Key Laboratory of Ministry of Education for Improving Oil and Gas Recovery(NEPUEOR-2022-03)Research and Innovation Fund for Graduate Students of Southwest Petroleum University(No.2022KYCX005)。
文摘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.
基金the National Natural Science Foundation of China(No.51974042)National Key Research and Development Program of China(No.2023YFC3009005).
文摘Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based on engineering properties to simulate the gradual collapse of the roof during longwall top coal caving(LTCC).A numerical model is established using the material point method(MPM)and the strain-softening damage constitutive model according to the structure of the physical model.Numerical simulations are conducted to analyze the LTCC process under different hard roofs for ground hydraulic fracturing.The results show that ground hydraulic fracturing releases the energy and stress of the target stratum,resulting in a substantial lag in the fracturing of the overburden before collapse occurs in the hydraulic fracturing stratum.Ground hydraulic fracturing of a low hard roof reduces the lag effect of hydraulic fractures,dissipates the energy consumed by the fracture of the hard roof,and reduces the abutment stress.Therefore,it is advisable to prioritize the selection of the lower hard roof as the target stratum.
基金sponsored by the General Program of the National Natural Science Foundation of China(Grant Nos.52079129 and 52209148)the Hubei Provincial General Fund,China(Grant No.2023AFB567)。
文摘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.
基金SINOPEC's Scientific and Technological Research Project:Research on effective production strategies of Jurassic continental shale oil and gas(No.P21078-5).
文摘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.
基金financially supported by the National Key R&D Program of China(Grant No.2018YFB1501802)。
文摘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.
基金The authors are grateful to the financial support from China Postdoctoral Science Foundation(2022M712645)Opening Fund of Key Laboratory of Enhanced Oil Recovery(Northeast Petroleum University),Ministry of Education(NEPU-EOR-2021-03).
文摘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.
基金supported by the Natural Science Foundation of China under (Nos. 42172293, 4190020747, and 41472268)。
文摘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.
基金Supported by the National Natural Science Foundation of China (51974332)。
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.52274038,5203401042174143)+1 种基金the Taishan Scholars Project(No.tsqnz20221140)the Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University)of China(No.PLN2020-5)。
文摘This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.
基金Supported by the National Natural Science Foundation of China(52192622,52204005,U20A20265)Sichuan Outstanding Young Scientific and Technological Talents Project(2022JDJQ0007).
文摘In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.
文摘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.
基金This work was financially supported by the National Key R&D Program of China(No.2021YFC2900500)the National Natural Science Foundation of China(Nos.52074020 and 42202306)+2 种基金the Open Fund of State Key Laboratory of Water Resource Protection and Utilization in Coal Mining(No.WPUKFJJ2019-06)the Interdisciplinary Research Project for Young Teachers of USTB(Fundamental Research Funds for the Central Universities)(No.FRF-IDRY-21001)the Natural Science Foundation of Jiangsu Province,China(No.BK20200993).
文摘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.
基金This work is supported by Open Fund Project“Study on Multiphase Flow Semi-Analytical Method for Horizontal Wells of Continental Shale Condensate Gas”of Sinopec Key Laboratory of Shale Oil/Gas Exploration and Production Technology.
文摘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.
文摘After volume fracturing of horizontal wells in shale gas reservoir, an extremely complex fracture system is formed. The space area of the fracture system is the reservoir reconstruction volume of shale gas reservoir. The geometric parameters such as crack length, crack width, crack height, and characteristic parameters such as crack permeability and fracture conductivity proposed for a single crack in conventional fracturing are insufficient to describe and characterize the complex network fracture system after volume fracturing. In this paper, the discrete fracture modeling method is used to establish the volume fracturing network fracture model of horizontal wells in shale gas reservoir by using the random modeling method within the determined reservoir space. The model is random and selective, and can fully provide different forms of volume fracturing fracture expansion, such as conventional fracture morphology, line network model and arbitrarily distributed network fractures. The research results provide a theoretical basis for the development plan and stimulation plan of shale gas reservoir, and have important reference value and significance for other unconventional gas reservoir fracturing.
基金supported by the National Natural Science Foundation of China (Grant No.42090023)the Alliance of International Science Organization (ANSO)Scholarship for Young Talents+3 种基金the Key Deployment Program of Chinese Academy of Sciences (YJKYYQ20190043,ZDBS-LY-DQC003,KFZD-SW-422,ZDRW-ZS-2021-3-1)the Scientific Research and Technology Development Project of China National Petroleum Corpo ration (2022DJ5503)the CAS Key Technology Talent ProgramSupercomputing Laboratory,IGGCAS。
文摘In-situ conversion presents a promising technique for exploiting continental oil shale formations,characterized by highly fractured organic-rich rock.A 3D in-situ conversion model,which incorporates a discrete fracture network,is developed using a self-developed thermal-flow-chemical(TFC)simulator.Analysis of the model elucidates the in-situ conversion process in three stages and defines the transformation of fluids into three distinct outcomes according to their end stages.The findings indicate that kerogen decomposition increases fluid pressure,activating fractures and subsequently enhancing permeability.A comprehensive analysis of activated fracture permeability and heating power reveals four distinct production modes,highlighting that increasing heating power correlates with higher cumulative fluid production.Activated fractures,with heightened permeability,facilitate the mobility of heavy oil toward production wells but hinder its cracking,thereby limiting light hydrocarbon production.Additionally,energy efficiency research demonstrates the feasibility of the in-situ conversion in terms of energy utilization,especially when considering the surplus energy from high-fluctuation energy sources such as wind and solar power to provide heating.
