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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 st...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.展开更多
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.展开更多
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.展开更多
Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining indus-tries.However,accurately simulating the fluid flow in realistic fracture networks(compared to t...Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining indus-tries.However,accurately simulating the fluid flow in realistic fracture networks(compared to the statistical fracture net-works)is still challenging due to the fracture complexity and computational burden.This work proposes a simple yet efficient numerical framework for the flow simulation in fractured porous media obtained by 3D high-resolution images,aiming at both computational accuracy and efficiency.The fractured rock with complex fracture geometries is numerically constructed with a cell-based discrete fracture-matrix model(DFM)having implicit fracture apertures.The flow in the complex fractured porous media(including matrix flow,fracture flow,as well as exchange flow)is simulated with a pipe-based cell-centered finite volume method.The performance of this model is validated against analytical/numerical solutions.Then a lab-scale true triaxial hydraulically fractured shale sample is reconstructed,and the fluid flow in this realistic fracture network is simu-lated.Results suggest that the proposed method achieves a good balance between computational efficiency and accuracy.The complex fracture networks control the fluid flow process,and the opened natural fractures behave as primary fuid pathways.Heterogeneous and anisotropic features of fluid flow are well captured with the present model.展开更多
Fracture networks within hot dry rock(HDR)geothermal reservoirs are complex,and heat extraction via water injection is thus a coupled process of heat-fluid-solid multifield.In this paper,utilizing the theory of normal...Fracture networks within hot dry rock(HDR)geothermal reservoirs are complex,and heat extraction via water injection is thus a coupled process of heat-fluid-solid multifield.In this paper,utilizing the theory of normally distributed random functions,we develop a corresponding pre-processing subprogram to establish a discrete network model of complex fracture distribution in HDR reservoirs;then construct a heat-fluid-solid finite element model for heat extraction via water injection and compare the numerical solution with the analytical solution of the one-dimensional non-isothermal consolidation problem for verification.The numerical simulation results show that the main factors affecting the heat extraction efficiency of HDR reservoirs include fracture width,fracture density,fracture permeability,and matrix permeability.When a HDR reservoir is injected with water for heat extraction,there is a certain threshold value of these influential parameters,beyond which the outlet temperature drops significantly,resulting in an obvious thermal breakthrough.When injecting water for heat extraction,the values of these parameters should be controlled and kept at a reasonable level,otherwise,the HDR reservoir may enter a thermal breakthrough stage in advance,which is not conducive for long-period heat extraction.Influenced by the random distribution of complex fractures,the leading edge of the cold front may present an irregular distribution.During the process of heat extraction,the stress gradually changes from a compressional state to a tensile state,which induces further damage to the HDR reservoir.展开更多
Based on the comprehensive understanding on microfractures and matrix pores in reservoir rocks,numerical algorithms are used to construct fractured porous media and fracture-pore media models.Connectivity coefficient ...Based on the comprehensive understanding on microfractures and matrix pores in reservoir rocks,numerical algorithms are used to construct fractured porous media and fracture-pore media models.Connectivity coefficient and strike factor are introduced into the models to quantitatively characterize the connectivity and strike of fracture network,respectively.The influences of fracture aperture,fracture strike and fracture connectivity on the permeability of porous media are studied by using multi-relaxation-time lattice Boltzmann model to simulate fluid flow in them.The greater the strike factor and the smaller the tortuosity of the fractured porous media,the greater the permeability of the fractured porous media.The greater the connectivity coefficient of the fracture network is,the greater the permeability of the fracture-pore media is,and the more likely dominant channel effect occurs.The fracture network connectivity has stronger influence on seepage ability of fracture-pore media than fracture aperture and fracture strike.The tortuosity and strike factor of fracture network in fractured porous media are in polynomial relation,while the permeability and fracture network connectivity coefficient of the fracture-pore media meet an exponential relation.展开更多
Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the con...Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.展开更多
Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed ...Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.展开更多
Stress sensitivity is a very important index to understand the seepage characteristics of a reservoir.In this study,dedicated experiments and theoretical arguments based on the visualization of porous media are used t...Stress sensitivity is a very important index to understand the seepage characteristics of a reservoir.In this study,dedicated experiments and theoretical arguments based on the visualization of porous media are used to assess the effects of the fracture angle,spacing,and relevant elastic parameters on the principal value of the permeability tensor.The fracture apertures at different angles show different change rates,which influence the relative permeability for different sets of fractures.Furthermore,under the same pressure condition,the fractures with different angles show different degrees of deformation so that the principal value direction of permeability rotates.This phenomenon leads to a variation in the water seepage direction in typical water-injection applications,thereby hindering the expected exploitation effect of the original well network.Overall,the research findings in this paper can be used as guidance to improve the effectiveness of water injection exploitation in the oil field industry.展开更多
基金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(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.
基金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.
基金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.
基金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.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFC3080200)the National Natural Science Foundation of China(Grant No.42022053)the China Postdoctoral Science Foundation(Grant No.2023M731264).
文摘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.
基金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.
文摘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.
基金supported by the Natural Sciences and Engineering Research Council of Canada(NSERC)with NSERC/Energi Simulation Industrial Research Chair program,NSERC Discovery 341275,and CRDPJ 54389419.
文摘Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining indus-tries.However,accurately simulating the fluid flow in realistic fracture networks(compared to the statistical fracture net-works)is still challenging due to the fracture complexity and computational burden.This work proposes a simple yet efficient numerical framework for the flow simulation in fractured porous media obtained by 3D high-resolution images,aiming at both computational accuracy and efficiency.The fractured rock with complex fracture geometries is numerically constructed with a cell-based discrete fracture-matrix model(DFM)having implicit fracture apertures.The flow in the complex fractured porous media(including matrix flow,fracture flow,as well as exchange flow)is simulated with a pipe-based cell-centered finite volume method.The performance of this model is validated against analytical/numerical solutions.Then a lab-scale true triaxial hydraulically fractured shale sample is reconstructed,and the fluid flow in this realistic fracture network is simu-lated.Results suggest that the proposed method achieves a good balance between computational efficiency and accuracy.The complex fracture networks control the fluid flow process,and the opened natural fractures behave as primary fuid pathways.Heterogeneous and anisotropic features of fluid flow are well captured with the present model.
基金This work is financially supported by the National Science Foundation of China(Grant No.52192622,No.51936001,No.52274002,No.51804033 and No.U20A20265)Beijing Natural Science Foundation(Grant No.3222030)+2 种基金the PetroChina Science and Technology Innovation Foundation Project(2021DQ02e0201)the Award Cultivation Foundation from Beijing Institute of Petrochemical Technology(Grant No.BIPTACF-002)the Fund of the Beijing Municipal Education Commission(Grant No.22019821001).
文摘Fracture networks within hot dry rock(HDR)geothermal reservoirs are complex,and heat extraction via water injection is thus a coupled process of heat-fluid-solid multifield.In this paper,utilizing the theory of normally distributed random functions,we develop a corresponding pre-processing subprogram to establish a discrete network model of complex fracture distribution in HDR reservoirs;then construct a heat-fluid-solid finite element model for heat extraction via water injection and compare the numerical solution with the analytical solution of the one-dimensional non-isothermal consolidation problem for verification.The numerical simulation results show that the main factors affecting the heat extraction efficiency of HDR reservoirs include fracture width,fracture density,fracture permeability,and matrix permeability.When a HDR reservoir is injected with water for heat extraction,there is a certain threshold value of these influential parameters,beyond which the outlet temperature drops significantly,resulting in an obvious thermal breakthrough.When injecting water for heat extraction,the values of these parameters should be controlled and kept at a reasonable level,otherwise,the HDR reservoir may enter a thermal breakthrough stage in advance,which is not conducive for long-period heat extraction.Influenced by the random distribution of complex fractures,the leading edge of the cold front may present an irregular distribution.During the process of heat extraction,the stress gradually changes from a compressional state to a tensile state,which induces further damage to the HDR reservoir.
基金Supported by the Science and Technology Major Project of PetroChina(2016E-06)National Natural Science Foundation of China(U1562217)。
文摘Based on the comprehensive understanding on microfractures and matrix pores in reservoir rocks,numerical algorithms are used to construct fractured porous media and fracture-pore media models.Connectivity coefficient and strike factor are introduced into the models to quantitatively characterize the connectivity and strike of fracture network,respectively.The influences of fracture aperture,fracture strike and fracture connectivity on the permeability of porous media are studied by using multi-relaxation-time lattice Boltzmann model to simulate fluid flow in them.The greater the strike factor and the smaller the tortuosity of the fractured porous media,the greater the permeability of the fractured porous media.The greater the connectivity coefficient of the fracture network is,the greater the permeability of the fracture-pore media is,and the more likely dominant channel effect occurs.The fracture network connectivity has stronger influence on seepage ability of fracture-pore media than fracture aperture and fracture strike.The tortuosity and strike factor of fracture network in fractured porous media are in polynomial relation,while the permeability and fracture network connectivity coefficient of the fracture-pore media meet an exponential relation.
基金co-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project(GrantNo.2011ZX05004003)
文摘Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.
文摘Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.
基金This work is financially supported by the National Natural Science Foundation Project(No.51374222)National Major Project(No.2017ZX05032004-002)+2 种基金the National Key Basic Research&Development Program(No.2015CB250905)CNPC’s Major Scientific and Technological Project(No.2017E-0405)SINOPEC Major Scientific Research Project(No.P18049-1).
文摘Stress sensitivity is a very important index to understand the seepage characteristics of a reservoir.In this study,dedicated experiments and theoretical arguments based on the visualization of porous media are used to assess the effects of the fracture angle,spacing,and relevant elastic parameters on the principal value of the permeability tensor.The fracture apertures at different angles show different change rates,which influence the relative permeability for different sets of fractures.Furthermore,under the same pressure condition,the fractures with different angles show different degrees of deformation so that the principal value direction of permeability rotates.This phenomenon leads to a variation in the water seepage direction in typical water-injection applications,thereby hindering the expected exploitation effect of the original well network.Overall,the research findings in this paper can be used as guidance to improve the effectiveness of water injection exploitation in the oil field industry.
