Overview of numerical models for interactions between hydraulic fractures and natural fractures: Challenges and limitations

The intersection of natural fractures with hydraulic fractures results in formation of complex fracture networks, including non-planar fractures or multi-stranded fractures. On one hand, opening of these natural fractures improves productivity of the formation; on the other hand, coalescence of these fractures into a hydraulic fracture makes pressure analysis and prediction of fracture growth very complicated. Overall, interactions between natural fractures and hydraulic fractures pose more challenges in the fracturing design and its execution. Investigation and understanding of their interaction are crucial in achieving successful fracture treatments in formations with pre-existing natural fracture network. In this paper, we will review the numerical works that have been done in the last decade to model opening of natural fractures during hydraulic fracturing, focusing especially on mechanical models that address propagation of hydraulic fractures in naturally fractures reservoirs. Linear elastic fracture mechanics, cohesive element methods and continuum damage mechanics techniques utilized to understand interaction of hydraulic fractures with natural fractures are discussed here based on their capability to reproduce experimental results and field observations.     

Coupled Numerical Evaluations of the Geomechanical Interactions Between a Hydraulic Fracture Stimulation and a Natural Fracture System in Shale Formations

Due to the low permeability of many shale reservoirs, multi-stage hydraulic fracturing in horizontal wells is used to increase the productive, stimulated reservoir volume. However, each created hydraulic fracture alters the stress field around it, and subsequent fractures are affected by the stress field from previous fractures. The results of a numerical evaluation of the effect of stress field changes (stress shadowing), as a function of natural fracture and geomechanical properties, are presented, including a detailed evaluation of natural fracture shear failure (and, by analogy, the generated microseismicity) due to a created hydraulic fracture. The numerical simulations were performed using continuum and discrete element modeling approaches in both mechanical-only and fully coupled, hydro-mechanical modes. The results show the critical impacts that the stress field changes from a created hydraulic fracture have on the shear of the natural fracture system, which in-turn, significantly affects the success of the hydraulic fracture stimulation. Furthermore, the results provide important insight into: the role of completion design (stage spacing) and operational parameters (rate, viscosity, etc.) on the possibility of enhancing the stimulation of the natural fracture network (‘complexity’); the mechanisms that generate the microseismicity that occurs during a hydraulic fracture stimulation; and the interpretation of the generated microseismicity in relation to the volume of stimulated reservoir formation.     

Injection-Sensitive Mechanics of Hydraulic Fracture Interaction with Discontinuities

We develop a new analytical model, called OpenT, that solves the elasticity problem of a hydraulic fracture (HF) contact with a pre-existing discontinuity natural fracture (NF) and the condition for HF re-initiation at the NF. The model also accounts for fluid penetration into the permeable NFs. For any angle of fracture intersection, the elastic problem of a blunted dislocation discontinuity is solved for the opening and sliding generated at the discontinuity. The sites and orientations of a new tensile crack nucleation are determined based on a mixed stress- and energy-criterion. In the case of tilted fracture intersection, the finite offset of the new crack initiation point along the discontinuity is computed. We show that aside from known controlling parameters such stress contrast, cohesional and frictional properties of the NFs and angle of intersection, the fluid injection parameters such as the injection rate and the fluid viscosity are of first-order in the crossing behavior. The model is compared to three independent laboratory experiments, analytical criteria of Blanton, extended Renshaw?Pollard, as well as fully coupled numerical simulations. The relative computational efficiency of OpenT model (compared to the numerical models) makes the model attractive for implementation in modern engineering tools simulating hydraulic fracture propagation in naturally fractured environments.     

离散裂隙网络对页岩气井产能影响的数值模拟

页岩作为典型的非常规储层,基质孔隙小,渗透率极低,水平井多级水力压裂为其商业开发的主要手段。准确模拟页岩气产能,应同时考虑水力裂隙和天然裂隙的渗流。基于离散裂隙模型和等效连续模型建立页岩气渗流数学模型,利用有限元分析方法进行数值求解,研究不同走向裂隙组对页岩气井产能的影响。研究认为,页岩基质为气体的生产提供了主要气源,天然裂隙作为渗流的主要通道,将气体输送到水力裂缝,进而到达井筒。模拟结果表征,离散裂隙的渗流特征对于页岩气井的产能有重要影响。根据页岩储层的天然裂隙走向,可以优化相应的水平井方位。对于二维离散裂隙网络模型,水平井沿着2个裂隙组夹角的平分线更有利于生产。      

Application of Stochastic Fracture Network with Numerical Fluid Flow Simulations to Groundwater Flow Modeling in Fractured Rocks

ＩＮＴＲＯＤＵＣＴＩＯＮＧｒｏｕｎｄｗａｔｅｒｏｒｆｌｕｉｄｆｌｏｗｍｏｄｅｌｉｎｇｉｎｆｒａｃｔｕｒｅｄｒｏｃｋｓｉｓａｃｏｍｐｌｉｃａｔｅｄｔｈｅｏｒｅｔｉｃａｌａｎｄａｐｐｌｉｅｄｔｏｐｉｃ.Ｂｏｔｈｔｈｅｏｒｅｔｉｃａｌｌｙａｎｄｏｐｅｒａｔｉｏｎａｌｌｙ ,ｉｔｉｓｉｍｐｏｒｔａｎｔｉｎｍａｎｙｆｉｅｌｄｓｓｕｃｈａｓｇｅｏｌｏｇｉｃａｌａｎｄｈｙｄｒｏｇｅｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ ,ｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｅｔｒｏｌｅｕｍｅｎｇｉｎｅｅｒｉｎ…     

Simulation of groundwater flow in fractured rocks using a coupled model based on the method of domain decomposition

Groundwater flow in fractured rocks is modeled using a coupled model based on the domain decomposition method. In the model, the fractured porous medium is divided into two non-overlapping sub-domains. One is the rock matrix, in which the medium is described using a continuum model. The other consists of deep fractures and fissure zones, where the medium is described using a discrete fracture network (DFN) model. The two models are coupled through the continuity of the hydraulic heads and fluxes on the common boundaries. The coupled model is used to simulate groundwater flow in a hydropower station. The results show that the model simulates groundwater levels that are in agreement with the measured groundwater levels. Furthermore, the model’s parameters relating to deep fractures and fissure zones are verified by comparing three different models (the continuum model, coupled model, and DFN model). The results show that the coupled model can capture and duplicate the hydrogeological conditions in the study domain, whereas the continuum model overestimates and the DFN model underestimates the measured hydraulic heads. A sensitivity analysis shows that fracture aperture has a considerable effect on the groundwater level. So, when the fracture aperture is large, the coupled model or DFN model is more appropriate than the continuum model in the fracture domain.     

An explicitly coupled hydro‐geomechanical model for simulating hydraulic fracturing in arbitrary discrete fracture networks

Modeling hydraulic fracturing in the presence of a natural fracture network is a challenging task, owing to the complex interactions between fluid, rock matrix, and rock interfaces, as well as the interactions between propagating fractures and existing natural interfaces. Understanding these complex interactions through numerical modeling is critical to the design of optimum stimulation strategies. In this paper, we present an explicitly integrated, fully coupled discrete‐finite element approach for the simulation of hydraulic fracturing in arbitrary fracture networks. The individual physical processes involved in hydraulic fracturing are identified and addressed as separate modules: a finite element approach for geomechanics in the rock matrix, a finite volume approach for resolving hydrodynamics, a geomechanical joint model for interfacial resolution, and an adaptive remeshing module. The model is verified against the Khristianovich–Geertsma–DeKlerk closed‐form solution for the propagation of a single hydraulic fracture and validated against laboratory testing results on the interaction between a propagating hydraulic fracture and an existing fracture. Preliminary results of simulating hydraulic fracturing in a natural fracture system consisting of multiple fractures are also presented. Copyright © 2012 John Wiley & Sons, Ltd.     

A fracture mapping and extended finite element scheme for coupled deformation and fluid flow in fractured porous media

This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

Development of a full 3D numerical model to investigate the hydraulic fracture propagation under the impact of orthogonal natural fractures

Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.

     

Numerical Modeling of Natural Fracture Distributions in Shale

The production efficiency of shale gas is affected by the interaction between hydraulic and natural fractures. This study presents a simulation of natural fractures in shale reservoirs, based on a discrete fracture network (DFN) method for hydraulic fracturing engineering. Fracture properties of the model are calculated from core fracture data, according to statistical mathematical analysis. The calculation results make full use of the quantitative information of core fracture orientation, density, opening and length, which constitute the direct and extensive data of mining engineering. The reliability and applicability of the model are analyzed with regard to model size and density, a calculation method for dominant size and density being proposed. Then, finite element analysis is applied to a hydraulic fracturing numerical simulation of a shale fractured reservoir in southeastern Chongqing. The hydraulic pressure distribution, fracture propagation, acoustic emission information and in situ stress changes during fracturing are analyzed. The results show the application of fracture statistics in fracture modeling and the influence of fracture distribution on hydraulic fracturing engineering. The present analysis may provide a reference for shale gas exploitation.     

Numerical modeling of stress-permeability coupling in rough fractures

A numerical model is described for coupled flow and mechanical deformation in fractured rock. The mechanical response of rock joints to changes in hydraulic pressure is strongly influenced by the geometric characteristics of the joint surfaces. The concept of this work is to combine straightforward finite element solutions with complex and realistic fracture surface geometry in order to reproduce the non-linear stress-deformation-permeability coupling that is commonly observed in fractures. Building on the numerous studies that have expanded the understanding of the key parameters needed to describe natural rough-walled fractures, new methods have been developed to generate a finite element mesh representing discrete fractures with realistic rough surface geometries embedded in a rock matrix. The finite element code GeoSys/Rockflow was then used to simulate the coupled effects of hydraulic stress, mechanical stress, and surface geometry on the evolving permeability of a single discrete fracture. The modeling concept was experimentally verified against examples from the literature. Modeling results were also compared to a simple interpenetration model.     

Flow Deflection in Intersected Tubes

ＩＮＴＲＯＤＵＣＴＩＯＮＴｈｅｔｈｅｏｒｅｔｉｃａｌｃａｌｃｕｌａｔｉｏｎｓｏｆｗａｔｅｒｍｉｇｒａｔｉｏｎ ,ｗｈｅｔｈｅｒｉｎａｎｉｎｄｉｖｉｄｕａｌｆｒａｃｔｕｒｅｏｒｉｎａｆｉｓｓｕｒｅｎｅｔｗｏｒｋ ,ｈａｖｅａｌｌｂｅｅｎｂａｓｅｄｏｎｔｈｅｅｑｕａｔｉｏｎｆｏｒｖｉｓｃｏｕｓｆｌｏｗｂｅｔｗｅｅｎｐａｒａｌｌｅｌｐｌａｔｅｓ.ＷｉｌｓｏｎａｎｄＷｉｔｈｅｒ ｓｐｏｏｎ (1976 )ｆｉｒｓｔａｃｃｏｍｐｌｉｓｈｅｄａｈｙｄｒａｕｌｉｃｓｉｍｕｌａｔｉｏｎｏｆｆｒａｃｔｕｒｅｗａｔｅｒｃｒｏｓｓ…     

Numerical simulation of particle plugging in hydraulic fracture by element partition method

Hydraulic fracturing (HF) treatment often involves particle migration and is applied for propping or plugging fractures. Particle migration behaviors, e.g., bridging, packing, and plugging, significantly affect the HF process. Hence, it is crucial to effectively simulate particle migration. In this study, a new numerical approach is developed based on a coupled element partition method (EPM). The EPM is used to model natural and hydraulic fractures, in which a fracture is allowed to propagate across an element, thereby avoiding remeshing in fracture simulations. To characterize the water flow process in a fracture, a fully hydromechanical coupled equation is adopted in the EPM. To model particle transportation in fractures with water flow, each particle is treated as a discrete element. The particles move in the fracture as a result of being dragged by fluid. Their movement, contact, and packing behaviors are simulated using the discrete element method. To reflect the plugging effect, an equivalent aperture approach is proposed. Using this method, the particle migration and its effect on water flow are well simulated. The simulation results show that this method can effectively reproduce particle bridging, plugging, and unblocking in a hydraulic fracture. Furthermore, it is demonstrated that particle plugging significantly affects water flow in a fracture and hence the propagation of hydraulic fracture. This method provides a simple and feasible approach for the simulation of particle migration in a hydraulic fracture.     

煤层水力压裂应力与裂隙演化的细观规律

水力压裂作为煤层强化增透技术的一种,其应力演化特征及裂隙形态与扩展范围的判断尤为重要。采用离散元数值方法,以导向压裂为背景,建立水力压裂流固耦合模型;通过应力路径、裂纹热点图等手段,探究水力压裂过程中压裂排量、泊松比、天然裂隙密度对应力演化和裂隙演化的影响及其细观规律。结果表明:不同压裂排量下的应力演化方向及最终应力路径曲线形状有着明显的不同,低排量下裂隙附近的应力比值逐渐增大,而在高排量下先增大后减小;煤层泊松比越大,平均压裂半径越低,但对起裂时间及裂隙的扩展形态影响不明显;天然裂隙的发育情况对水力裂隙的扩展起着关键性作用,高裂隙发育煤层水力裂隙扩展的方向性无法预测,应力演化方向会出现反转现象;压裂过程中不同区域的应力演化特征能够反映出裂隙的扩展状态,现场可通过监测压裂区域附近应力变化,判断水力压裂缝网的扩展范围。      

Mechanical and hydraulic properties of rocks related to induced seismicity

Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23–55.The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass.In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid.However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now available that clearly demonstrates this stress-flow behavior. Two approaches have been used in attempting to numerically model such behavior: (1) continuum models, and (2) discrete models. The continuum approach only needs information as to average values of fracture spacing and material properties. But because of the inherent complexity of fractured rock masses and the corresponding decrease in symmetry, it is difficult to develop an equivalent continuum that will simulate the behavior of the entire system. The discrete approach, on the other hand, requires details of the fracture geometry and material properties of both fractures and rock matrix. The difficulty in obtaining such information has been considered a serious limitation of discrete models, but improved borehole techniques can enable one to obtain the necessary data, at least in shallow systems. The possibility of extending these methods to deeper fracture systems needs more investigation. Such data must be considered when deciding whether to use a continuum or discrete model to represent the interaction of rock and fluid forces in a fractured rock system, especially with regard to the problem of induced seismicity. When one is attempting to alter the pressure distribution in a fault zone by injection or withdrawal of fluids, the extent to which this can be achieved will be controlled in large measure by the behavior of the fractures that communicate with the borehole. Since this is essentially a point phenomenon, i.e., the changes will propagate from a relatively small region around the borehole, the use of a discrete model would appear to be preferable.     

Simulation of hydraulic fracture propagation near a natural fracture using virtual multidimensional internal bonds

A virtual multidimensional internal bond (VMIB) model developed to simulate the propagation of hydraulic fractures using the finite‐element method is formulated within the framework of the virtual internal bond theory (VIB) that considers a solid as randomized material particles in the micro scale, and derives the macro constitutive relation from the cohesive law between the material particles with an implicit fracture criterion. Hydraulic pressure is applied using a new scheme that enables simulation of hydraulically driven cracks. When the model is applied to study hydraulic fracture propagation in the presence of a natural fracture, the results show the method to be very effective. It shows that although the in situ stress ratio is the dominant factor governing the propagation direction, a natural fault can also strongly influence the hydraulic fracture behavior. This influence is conditioned by the shear stiffness of the fault and the distance to the original hydraulic fracture. The model results show that when the fault is strong in shear, its impact on hydraulic fracture trajectory is weak and the hydraulic fracture will likely penetrate the fault. For a weak fault, however, the fracture tends to be arrested at the natural fault. The distance between the fault and the hydraulic fracture is also important; the fault influence increases with decreasing distance. The VMIB does not require selection of a fracture criterion and remeshing when the fracture propagates. Therefore, it is advantageous for modeling fracture initiation and propagation in naturally fractured rock. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.     

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr–Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils.     

A Comprehensive Numerical Study of Hydraulic Fracturing Process and Its Affecting Parameters

Production of hydrocarbons from shale formations has been made possible mainly due to hydraulic fracture (HF) technology. It increases the permeability of reservoir rocks around a well by connecting fractures and improving conductivity. HF behavior especially in presence of natural fractures have recently given much attention in studies. This paper investigates HF propagation and its associated parameters in various conditions. A higher order displacement discontinuity method is used to achieve higher accuracy in the results. First, behavior of crack opening displacement (COD) of an HF i.e. HF width in various conditions is studied. COD is a key parameter in determination of an HF operation success. It is proportional to production rate of oil and gas wells and provides a path for proppant entrance into the fractures. An equation considering many important parameters, based on numerous numerical modellings of various mechanical and geometrical effects on COD is proposed with coefficient of determination and standard error of 94.35% and 4.37 × 10^?4 respectively. The next part of the paper studies the HF propagation in a naturally fractured reservoir. These natural fractures alter the stress fields and hence affect the propagation of a hydraulic fracture. In fact, it is shown that in certain orientations of hydraulic fractures and natural fractures, the effect of natural fractures disappear or completely changes propagation path. Using a combination of several interaction criteria, a new modeling of HFs and NFs interaction is presented. The modellings showed that spacing and angle of intersection can significantly affect HFs propagation. The results of COD and HF propagation in presence of natural fractures may be considered in HG design and primary orientation of perforated fractures.