Formation of X-shaped hydraulic fractures in deep thick glutenite reservoirs:A case study in Bohai Bay Basin,East China

Tight glutenite reservoirs are widely developed in Bohai Bay Basin,East China.They are mostly huge thick and rely on hydraulic fracturing treatment for commercial exploitation.To investigate the propagation behavior of hydraulic fractures in these glutenite reservoirs,the geological feature of reservoirs in Bohai Bay Basin is studied firstly,including the reservoir vertical distribution feature and the heterogeneous lithology.Then,hydraulic fracturing treatments in block Yan 222 are carried out and the fracturing processes are monitored by the microseismic system.Results show the hydraulic fractures generated in the reservoirs are mostly in X shape.The cause of X-shaped hydraulic fractures in this study is mainly ascribed to(I)the reservoir heterogeneity and(II)the stress shadow effect of two close hydraulic fractures propagating in the same orientation,which is confirmed by the following numerical simulation and related research in detail.This study can provide a reference for the research on the fracturing behavior of the deep thick glutenite reservoirs.

Experimental Investigation on Condensate Revaporization During Gas Injection Development in Fractured Gas Condensate Reservoirs

The gas field in the Bohai Bay Basin is a fractured metamorphic buried-hill reservoir with dual-media characteristics. The retrograde vaporization mechanism observed in this type of gas condensate reservoir differs significantly from that observed in sand gas condensate reservoirs. However, studies on improving the recovery of fractured gas condensate reservoirs are limited;thus, the impact of retrograde vaporization on condensate within fractured metamorphic buried-hill reservoirs remains unclear. To address this gap, a series of gas injection experiments are conducted in pressure-volume-temperature(PVT) cells and long-cores to investigate the retrograde vaporization effect of condensate using different gas injection media in fractured gas condensate reservoirs. We analyze the variation in condensate volume, gas-to-oil ratio, and condensate recovery during gas injection and examine the influence of various gas injection media(CO_(2), N_(2), and dry gas) under different reservoir properties and varying gas injection times. The results demonstrate that the exchange of components between injected gas and condensate significantly influences condensate retrograde vaporization in the formation. Compared with dry gas injection and N_(2) injection,CO_(2) injection exhibits a superior retrograde vaporization effect. At a CO_(2) injection volume of 1 PV, the percentage shrinkage volume of condensate is 13.82%. Additionally, at the maximum retrograde condensation pressure, CO_(2) injection can increase the recovery of condensate by 22.4%. However, the condensate recovery is notably lower in fractured gas condensate reservoirs than in homogeneous reservoirs, owing to the creation of dominant gas channeling by fractures, which leads to decreased condensate recovery. Regarding gas injection timing, the effect of gas injection at reservoir pressure on improving condensate recovery is superior to that of gas injection at the maximum retrograde condensation pressure. This research provides valuable guidance for designing gas injection development plans and dynamic tracking adjustments for fractured gas condensate reservoirs.

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

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

A multi-mechanism numerical simulation model for CO_(2)-EOR and storage in fractured shale oil reservoirs

Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.

Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

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

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

Evaluation of Well Spacing for Primary Development of Fractured Horizontal Wells in Tight Sandstone Gas Reservoirs

Methods for horizontal well spacing calculation in tight gas reservoirs are still adversely affected by the complexity of related control factors,such as strong reservoir heterogeneity and seepage mechanisms.In this study,the stress sensitivity and threshold pressure gradient of various types of reservoirs are quantitatively evaluated through reservoir seepage experiments.On the basis of these experiments,a numerical simulation model(based on the special seepage mechanism)and an inverse dynamic reserve algorithm(with different equivalent drainage areas)were developed.The well spacing ranges of Classes I,II,and III wells in the Q gas field are determined to be 802–1,000,600–662,and 285–400 m,respectively,with their average ranges as 901,631,and 342.5 m,respectively.By considering both the pairs of parallel well groups and series well groups as examples,the reliability of the calculation results is verified.It is shown that the combination of the two models can reduce errors and provide accurate results.

Naturally fractured hydrocarbon reservoir simulation by elastic fracture modeling

Accurate fluid flow simulation in geologically complex reservoirs is of particular importance in construction of reservoir simulators.General approaches in naturally fractured reservoir simulation involve use of unstructured grids or a structured grid coupled with locally unstructured grids and discrete fracture models.These methods suffer from drawbacks such as lack of flexibility and of ease of updating.In this study,I combined fracture modeling by elastic gridding which improves flexibility,especially in complex reservoirs.The proposed model revises conventional modeling fractures by hard rigid planes that do not change through production.This is a dubious assumption,especially in reservoirs with a high production rate in the beginning.The proposed elastic fracture modeling considers changes in fracture properties,shape and aperture through the simulation.This strategy is only reliable for naturally fractured reservoirs with high fracture permeability and less permeable matrix and parallel fractures with less cross-connections.Comparison of elastic fracture modeling results with conventional modeling showed that these assumptions will cause production pressure to enlarge fracture apertures and change fracture shapes,which consequently results in lower production compared with what was previously assumed.It is concluded that an elastic gridded model could better simulate reservoir performance.

Pressure transient analysis of a finite-conductivity multiple fractured horizontal well in linear composite gas reservoirs

Faulted gas reservoirs are very common in reality,where some linear leaky faults divide the gas reservoir into several reservoir regions with distinct physical properties.This kind of gas reservoirs is also known as linear composite(LC)gas reservoirs.Although some analytical/semi-analytical models have been proposed to investigate pressure behaviors of producing wells in LC reservoirs based on the linear composite ideas,almost all of them focus on vertical wells and studies on multiple fractured horizontal wells are rare.After the pressure wave arrives at the leaky fault,pressure behaviors of multiple fractured horizontal wells will be affected by the leaky faults.Understanding the effect of leaky faults on pressure behaviors of multiple fractured horizontal wells is critical to the development design.Therefore,a semi-analytical model of finite-conductivity multiple fractured horizontal(FCMFH)wells in LC gas reservoirs is established based on Laplace-space superposition principle and fracture discrete method.The proposed model is validated against commercial numerical simulator.Type curves are obtained to study pressure characteristics and identify flow regimes.The effects of some parameters on type curves are discussed.The proposed model will have a profound effect on developing analytical/semi-analytical models for other complex well types in LC gas reservoirs.

An intelligent prediction method of fractures in tight carbonate reservoirs

An intelligent prediction method for fractures in tight carbonate reservoir has been established by upgrading single-well fracture identification and interwell fracture trend prediction with artificial intelligence,modifying construction of interwell fracture density model,and modeling fracture network and making fracture property equivalence.This method deeply mines fracture information in multi-source isomerous data of different scales to reduce uncertainties of fracture prediction.Based on conventional fracture indicating parameter method,a prediction method of single-well fractures has been worked out by using 3 kinds of artificial intelligence methods to improve fracture identification accuracy from 3 aspects,small sample classification,multi-scale nonlinear feature extraction,and decreasing variance of the prediction model.Fracture prediction by artificial intelligence using seismic attributes provides many details of inter-well fractures.It is combined with fault-related fracture information predicted by numerical simulation of reservoir geomechanics to improve inter-well fracture trend prediction.An interwell fracture density model for fracture network modeling is built by coupling single-well fracture identification and interwell fracture trend through co-sequential simulation.By taking the tight carbonate reservoir of Oligocene-Miocene AS Formation of A Oilfield in Zagros Basin of the Middle East as an example,the proposed prediction method was applied and verified.The single-well fracture identification improves over 15%compared with the conventional fracture indication parameter method in accuracy rate,and the inter-well fracture prediction improves over 25%compared with the composite seismic attribute prediction.The established fracture network model is well consistent with the fluid production index.

Efficient production optimization for naturally fractured reservoir using EDFM

Naturally fractured reservoirs make important contributions to global oil and gas reserves and production.The modeling and simulation of naturally fractured reservoirs are different from conventional reservoirs as the existence of natural fractures.To address the development optimization problem of naturally fractured reservoirs,we propose an optimization workflow by coupling the optimization methods with the embedded discrete fracture model(EDFM).Firstly,the effective and superior performance of the workflow is verified based on the conceptual model.The stochastic simplex approximate gradient(StoSAG)algorithm,the ensemble optimization(EnOpt)algorithm,and the particle swarm optimization(PSO)algorithm are implemented for the production optimization of naturally fractured reservoirs based on the improved versions of the Egg model and the PUNQ-S3 model.The results of the two cases demonstrate the effectiveness of this optimization workflow by finding the optimal well controls which yield the maximum net present value(NPV).Compared to the initial well control guess,the final NPV obtained from the production optimization of fractured reservoirs based on all three optimization algorithms is significantly enhanced.Compared with the optimization results of the PSO algorithm,StoSAG and EnOpt have significant advantages in terms of final NPV and computational efficiency.The results also show that fractures have a significant impact on reservoir production.The economic efficiency of fractured reservoir development can be significantly improved by the optimization workflow.

A Review of the Dynamic Modeling Approaches for Characterizing Fluid Flow in Naturally Fractured Reservoirs

Fluid flow in fractured media has been studied for decades and received considerable attention in the oil and gas industry because of the high productivity of naturally fractured reservoirs.Due to formation complexity and reservoir heterogeneity,characterizing fluid flow with an appropriate reservoir model presents a challenging task that differs relatively from homogeneous conventional reservoirs in many aspects of view,including geological,petrophysical,production,and economics.In most fractured reservoirs,fracture networks create complex pathways that affect hydrocarbon flow,well performance,hence reservoir characterization.A better and comprehensive understanding of the available reservoir modeling approaches is much needed to accurately characterize fluid flow behavior in NFRs.Therefore,in this paper,a perspective review of the available modeling approaches was presented for fluid flow characterization in naturally fractured medium.Modeling methods were evaluated in terms of their description,application,advantages,and disadvantages.This study has also included the applications of these reservoir models in fluid flow characterizing studies and governing equations for fluid flow.Dual continuum models were proved to be better than single continuum models in the presence of large scale fractures.In comparison,discrete models were more appropriate for reservoirs that contain a smaller number of fractures.However,hybrid modeling was the best method to provide accurate and scalable fluid flow modeling.It is our understanding that this paper will bridge the gap between the fundamental understanding and application of NFRs modeling approaches and serve as a useful reference for engineers and researchers for present and future applications.

The application study on the multi-scales integrated prediction method to fractured reservoir description

In this paper,we implement three scales of fracture integrated prediction study by classifying it to macro-（ 1/4/λ）,meso-（ 1/100λ and 1/4λ） and micro-（ 1/100λ） scales.Based on the multi-scales rock physics modelling technique,the seismic azimuthal anisotropy characteristic is analyzed for distinguishing the fractures of meso-scale.Furthermore,by integrating geological core fracture description,image well-logging fracture interpretation,seismic attributes macro-scale fracture prediction and core slice micro-scale fracture characterization,an comprehensive multi-scale fracture prediction methodology and technique workflow are proposed by using geology,well-logging and seismic multi-attributes.Firstly,utilizing the geology core slice observation（Fractures description） and image well-logging data interpretation results,the main governing factors of fracture development are obtained,and then the control factors of the development of regional macro-scale fractures are carried out via modelling of the tectonic stress field.For the meso-scale fracture description,the poststack geometric attributes are used to describe the macro-scale fracture as well,the prestack attenuation seismic attribute is used to predict the meso-scale fracture.Finally,by combining lithological statistic inversion with superposed results of faults,the relationship of the meso-scale fractures,lithology and faults can be reasonably interpreted and the cause of meso-scale fractures can be verified.The micro-scale fracture description is mainly implemented by using the electron microscope scanning of cores.Therefore,the development of fractures in reservoirs is assessed by valuating three classes of fracture prediction results.An integrated fracture prediction application to a real field in Sichuan basin,where limestone reservoir fractures developed,is implemented.The application results in the study area indicates that the proposed multi-scales integrated fracture prediction method and the technique procedureare able to deal with the strong heterogeneity and multi-scales problems in fracture prediction.Moreover,the multi-scale fracture prediction technique integrated with geology,well-logging and seismic multi-information can help improve the reservoir characterization and sweet-spots prediction for the fractured hydrocarbon reservoirs.

Borehole stability in naturally fractured reservoirs luring production tests

Based on the plane of weakness theory, a model for predicting borehole stability in fractured reservoirs under different stress states was estiblisted and the equations for solving borehole stability were developed. The minimum downhole pressures required to maintain borehole stability under different natural fracture occurrences were calculated by using the data from a well in the Tazhong （central Tarim） area, Tarim Basin, west China. Several conclusions were drawn for naturally fractured reservoirs with a dip angle from less than 10° to greater than 30°. Application in three wells in the Tazhong area indicates that this model is practically useful.

Experimental investigation of propagation mechanisms and fracture morphology for coalbed methane reservoirs

Fracture propagation mechanisms in coalbed methane(CBM) reservoirs are very complex due to the development of the internal cleat system. In this paper, the characteristics of initiation and propagation of hydraulic fractures in coal specimens at different angles between the face cleat and the maximum horizontal principal stress were investigated with hydraulic fracturing tests. The results indicate that the interactions between the hydraulic fractures and the cleat system have a major effect on fracture networks. "Step-like’’ fractures were formed in most experiments due to the existence of discontinuous butt cleats. The hydraulic fractures were more likely to divert or propagate along the butt cleat with an increase in the angles and a decrease in the horizontal principal stress difference. An increase in the injection rate and a decrease in the fracturing fluid viscosity were more conducive to fracture networks. In addition, the influence on fracture propagation of the residual coal fines in the wellbore was also studied. The existence of coal fines was an obstacle in fracturing, and no effective connection can be formed between fractures. The experimental investigation revealed the fracture propagation mechanisms and can provide guidance for hydraulic fracturing design of CBM reservoirs.

Modeling flow in naturally fractured reservoirs： effect of fracture aperture distribution on dominant sub-network for flow

Fracture network connectivity and aperture （or conductivity） distribution are two crucial features controlling flow behavior of naturally fractured reservoirs. The effect of connectivity on flow properties is well documented. In this paper, however, we focus here on the influence of fracture aperture distribution. We model a two dimensional fractured reservoir in which the matrix is impermeable and the fractures are well connected. The fractures obey a power-law length distribution, as observed in natural fracture networks. For the aperture distribution, since the information from subsurface fracture networks is limited, we test a number of cases： log-normal distributions （from narrow to broad）, power-law distributions （from narrow to broad）, and one case where the aperture is pro- portional to the fracture length. We find that even a well- connected fracture network can behave like a much sparser network when the aperture distribution is broad enough （c~ 〈 2 for power-law aperture distributions and σ ≥ 0.4 for log-normal aperture distributions）. Specifically, most fractures can be eliminated leaving the remaining dominant sub-network with 90% of the permeability of the original fracture network. We determine how broad the aperture distribution must be to approach this behavior and the dependence of the dominant sub-network on the parameters of the aperture distribution. We also explore whether one can identify the dominant sub-network without doing flow calculations.

The tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy: A case study in the Wangyao Oilfield of Ordos Basin, China

Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantitative information on a large number of fractures, based on which the accuracy of subsurface fracture characterization can be improved. Here we take the tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy, a case study of the Chang6t~ Formation of the Upper Triassic Yanchang Group of the Wangyao Oilfield in the Ordos Basin of China. An outcrop at the edge of the basin is a suitable analog for the reservoir, but the prerequisite is that they must have equivalent previous stress fields, similar final structural characteristics, relative timing and an identical depositional environment and diagenesis. The relationship among fracture density, rock type and bed thickness based on the outcrop is one of the most important fracture distribution models, and can be used to interpret fracture density in individual wells quantitatively. Fracture orientation, dip, geometry and scale, also should be described and measured in the outcrop, and can be used together with structure restoration and single well fracture density interpretation to guide fracture intensity prediction on bed surfaces and to constrain the construction of the 3D fracture geometry model of the subsurface reservoir. The application of the above principles shows the outcrop-based tectonic fracture models of the target ultra-low permeability sandstone reservoir are consistent with fractures inferred from microseismic interpretation and tracer tests. This illustrated that the fracture modeling based on the outcrop analogy is reliable and can reduce the uncertainty in stochastic fracture modeling.

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

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

Anisotropic characteristics of electrical responses of fractured reservoir with multiple sets of fractures

In fractured reservoirs, the fractures not but also form the main flow channels which connect productivity of reservoirs. However, because of the only provide the storage space for hydrocarbons, the pores of the matrix, so fractures dominate the heterogeneity and randomness of the distribution of fractures, exploration and evaluation of fractured reservoirs is still one of the most difficult problems in the oil industry. In recent years, seismic anisotropy has been applied to the assessment of fractured formations, whereas electrical anisotropy which is more intense in fractured formations than seismic anisotropy has not been studied or used so extensively. In this study, fractured reservoir models which considered multiple sets of fractures with smooth and partly closed, rough surfaces were established based on the fractures and pore network, and the vertical and horizontal electrical resistivities were derived as a function of the matrix and fracture porosities according to Ohm＇s law. By using the anisotropic resistivity equations, variations of the electrical anisotropy of three types of fractured models under the conditions of free pressure and confining pressure were analyzed through the variations of the exerted pressure, matrix porosity, fracture aperture and formation water resistivity. The differences of the vertical and horizontal resistivities and the anisotropy between the connected and non-connected fractures were also analyzed. It is known from the simulated results that an increase of the confining pressure causes a decrease of electrical anisotropy because of the elasticity of the closed fractures and the decrease of the fracture aperture. For a fixed fracture porosity, the higher the matrix porosity, the weaker the electrical anisotropy in the rock formation.

The role of natural fracture activation in hydraulic fracturing for deep unconventional geo-energy reservoir stimulation

The presence of sealed or semi-sealed,multiscale natural fracture systems appears to be crucial for the successful stimulation of deep reservoirs.To explore the reaction of such systems to reservoir stimulation,a new numerical simulation approach for hydraulic stimulation has been developed,trying to establish a realistic model of the physics involved.Our new model successfully reproduces dynamic fracture activation,network generation,and overall reservoir permeability enhancement.Its outputs indicate that natural fractures facilitate stimulation far beyond the near-wellbore area,and can significantly improve the hydraulic conductivity of unconventional geo-energy reservoirs.According to our model,the fracture activation patterns are jointly determined by the occurrence of natural fractures and the in situ stress.High-density natural fractures,high-fluid pressure,and low effective stress environments promote the formation of complex fracture networks during stimulation.Multistage or multicluster fracturing treatments with an appropriate spacing also increase the stimulated reservoir area(SRA).The simulation scheme demonstrated in this work offers the possibility to elucidate the complex multiphysical couplings seen in the field through detailed site-specific modeling.