Horizontal versus vertical wells interference in hydraulically fractured shale reservoirs

The impact of well interference on in-situ stresses, drainage area, and pressure response in hydraulically fractured shale reservoirs is examined. In-situ stress distribution in the reservoir resulting from fracture propagation and poro-elasticity and its influence on hydraulic fracture orientation and well spacing configurations are studied using iterative numerical methods. The results of the simulation indicated that drainage distance (XDL) from the well centre is restricted to the immediate environment of the well and with little effect on the external reservoir. In contrast, in-situ stress change has a wider and more complex reservoir reach away from the well with stress orthogonal reorientation occurring from a distance-of-stress-orthogonality (Ⱶσ), while pressure response has the farthest reach (XPT). A new approach utilised in this study, which considers in-situ stress, drainage area and pressure interference (such that XDL< Ⱶσ< XPT), suggests that a spacing range of 450ft to 750ft, with an optimum of 600ft for minimal interference will be adequate. Furthermore, parallel orientation of infill wells within this range is less feasible due to complex stress reorientation over the productive years. Wells drilled and fractured perpendicular to the parent well showed incremental cumulative production.


Introduction
Well interference induced by hydraulic fracturing can affect the development of shale reservoirs. Hydraulic fracture stimulation has been employed successfully to enhance production in very low permeability reservoirs that are in themicro or nano midrange (usually less than 1md). Initially implemented on vertical wells, it was uneconomical to implement in low permeability shale reservoirs because of the high cost of drilling many vertical wells. 1 Today, many reports have shown that the use of multistage hydraulically fractured horizontal wells proven to be a better option. Although, certain reservoir geometry and in-situ stress condition require vertical wells to be drilled; further advancement in this technology, specifically, in hydraulic fracturing design in horizontal wells and well spacing, increased its applicability in various shale plays.
Unlike fracture spacing optimisation, well spacing has no definite formula, but has been modelled to suit pressure profile and cumulative hydrocarbon output. In Eagle Ford shale, well spacing averages at 700ft, but recently, it has become tighter. 2 Reduced hydraulic fracture spacing with high fracture density and well density have been revealed to improverecovery. [3][4][5] However, these studies have shown that closer well spacing in hydraulically fractured wells results in well-to-well interference. This may affect well performance overall development strategy.
In-situ stress, but not pressure, in the stress shadow and fracture drainage volume was the focus of earlier performance studies and resulted in the Zipper and Texas Two Step well and fracture patterns. In contrast, well spacing studies have centred on well drainage volume and pressure, but not in-situ stress. 5 The implication of this is unfavourable hydraulic fracture stimulation of infill well leading to under optimised production. 6 The modified zipper configuration attempted stress inclusion in its design to address this gap, but it is not widely implemented for interference reasons. 7 The inclusion of in-situ stress distribution in addition to drainage volume and pressure transient will provide a holistic understanding of well-towell interference in shale reservoirs. This will help in placement and stimulation of wells.
In this study, the effect of in-situ stresses and thermo-poroelasticity on production from shale reservoirs was investigated. A coupled geomechanics, heat and fluid flow model for hydrocarbon production in shale reservoirs is adopted. A simple reservoir model with two wells (one horizontal and one vertical)having several planar fractures was developed. We investigated drainage volume, pressure response and in-situ stress distribution of well interference. Sensitivity analysis was performed on the fracture orientation to determine the impact on fluid drainage. Our study provides key insights into well spacing and fracture orientation when developing infill wells.

Model formulation
The governing equations for hydrocarbon production in unconventional reservoirs include porous medium deformation, pore and fracture fluid flow, and fracture propagation. These equations are expressions of mass, momentum and energy equilibrium. Thermal energy changes due to production or injection have shown significant effect on in-situ stresses, thus heat has been accounted for in the model. 8

Geomechanical model
Typical mechanical properties: Young's modulus (M), Shear modulus (G), Poisson's ratio(v), Bulk Modulus (K) and Biot's coefficient (α), all makeup the fundamental component of a geomechanical model. Also, rock strength, which is the amount of resistance a body of rock must wedge against external energy, plays a significant role in fracture initiation and propagation. Unconfined compressive strength (UCS), Tensile strength (T o ), cohesion (C o ), and internal friction (θ), all make up different forms for rock strength that come into force to resist deformation in rocks

Principal stresses
These stresses play a major role in determining direction of wellbore during drilling and the direction of hydraulic fracture propagation. Zhao et al. 9 demonstrated that hydraulic fracture reroutes its path towards the direction of maximum horizontal stress direction, irrespective of its direction of initiation. Implication of this travel path is the formation of simple non-planar fractures at early time with low fracture density, and complex fractures at late time with high fracture density. Thus, it is generally recommended that fracture stimulation be initiated in the direction perpendicular to the minimum horizontal stress.
Principal stresses acting on the subsurface; minimum and maximum horizontal stresses, are a translation of the effect of vertical stress (loading) at the depth of interest. Typical change in stress for anisotropic formation can be written as 1 : (2)

Porous medium deformation
Porous medium deformation is modelled as a poroelastic material undergoing quasi-static deformation. At equilibrium, the initial stress in the medium is zero when body forces are neglected: where, . ∇ is the divergence operator, σ is the Cauchy total stress tensor and b F g ρ = is the body force, comprising of, g the gravity vector and, b ρ the bulk density., ρ is the density of the solid phase, and φ is the true porosity. The effective stress and rock constitutive relation for thermo-elasticity is given as 10 : where, ' σ is effective stress, α is the Biot's coefficient, s β is the coefficient of liner solid thermal expansion, dr C is the rank-4 drained elasticity tensor, dr K is the drained bulk modulus, 1 is the rank-2 identity tensor, p is the pore pressure and T is temperature. The incremental stress form that includes pore pressure and temperature effects for the geomechanical model can be written as: When subjected to small strain ε , the medium deforms to a new stress state: The strain tensor, ε due to the infinitesimal transformation assumption, is the symmetric gradient of displacement : u ( ) and the resulting volumetric strain v ε is given as: .

Flow model
For single-phase flow of a slightly compressible fluid in a poroelastic medium, the fluid mass conservation equation can be written as:

Pore fluid flow
Under the assumption of small volumetric strains and Darcy flow, the thermoporoelasticity equation can be written by substituting Biot's poroelastic constitutive equation 10 : Where the velocity of fluid through the network of interconnected pores is given as:

Geothermal effect
In-situ temperature change during production or injection may account for additional change in in-situ stress. We have included reservoir temperature effect in the porous deformation and pore fluid flow equations above. The solution for stepwise change in temperature can be obtained by using Fourier's law as the heat conduction term and a thermoporoelastic constitutive equation 10 in the conservation equation (Eq. 10). Thus, the thermal diffusion equation is obtained as 8 : Where the total volumetric heat capacity is ( ) 1 .

Fracture initiation and propagation
In response to change in stress, fracture propagation is modelled in accordance with linear elastic fracture mechanics (LEFM). 11 The model predicts the amount of stress required to initiate and propagate fracture in terms of fracture geometry (r,) θ and stress intensity factor K I .

cos
1 sin sin 2 2 2 2 K I is proportional to applied stress. The rock's maximum tensile strength corresponds to critical intensity factor, K IC, (Fracture toughness) obtained at r = 0 and I K dampens asymptotically with increasing r until normal far field stress is achieved. At equilibrium, K IC is equal to applied intensity factor K I . Thus, fracture propagation occurs when applied stress is greater than sum of the rock's maximum tensile strength and in-situ stress parallel to the direction of fracture.
Stress intensity at Equilibrium: Equation (15) and (16) represents initiation and propagation criteria of fractures. Additional information on the equations for the corresponding fracture height growth, width profile, fluid pressure and stress intensity factor can be found in Yang et al. 12 and Weng et al. 13

Interference modelling
The effect of communication of reservoir fluids between wells is significant. Developed models in this regard measure pressure effect which is directly impactful on fluid recovery in conventional reservoirs. This is not the case for shale reservoirs where fluid flow is primarily dependent on matrix permeability. Unconventional reservoirs may have low flow tendency, but interference is still possible for varying low permeability. 14 Variable flow rate, amongst others, is characteristic of shale reservoirs and makes interpretation of well test difficult. Although, constant production rate is possible in vertical wells, it is uneconomical and inapplicable in gas well test. In horizontal wells, constant rate is practically impossible. 15 Complexities of shale reservoirs influences fluid flow behaviour and accounts for varying production rates and pressure. The proposed correlation for multi-rate flow Q n , and normalised pressure, (∆P/ Q n ) used to study and analyse production from shale formations was derived from instantaneous source function. The semi analytical model is presented in Appendix B. Pressure alone cannot determine well spacing in unconventional reservoirs because of its indirect impact on fluid flow. Thus, pressure interference in unconventional reservoirs would be assumed to occur in most cases. Consequently, interference as a result of other determinants like drainage volume and in-situ stress may be investigated and considered for optimum hydrocarbon recovery.

Iterative solution approach using CMG-GEM
In the iterative approach, the reservoir porosity Φ , is first imputed from initial predetermined value, then subsequently updated by solving the thermoporoelasticity model: The coefficients of reservoir porosity c i and a j are defined in appendix A Pressure and saturation change is obtained from the fluid model, while reservoir porosity and its derivatives from the geomechanical model. The sequence of execution commences with validating initial reservoir properties at equilibrium. A dynamic process simulates pressure change due to production/injection and stress change due to pore pressure change (internal process). It also measures change in stress due to fracture stimulation. Porosity is updated, and the cycle is repeated. An adaptive implicit method in discretised system of fluid flow and geomechanical equations is used to model pressure.

Model setup and verification
The model is a synthetic shale system in a reservoir with dimensions 3600ft x 3600ft x 100ft represented by a 24 x 24 x 1 grid system and block length and width of 150ft 2 each. The fluid system is atypical black oil system with low gas oil ratio (GOR) of Eagle Fordshale formation extracted from publicly available data contained in Table  1&2 obtained from Yu et al. 16

Strategy model-output focus
In order to achieve our objectives, the following outcomes from the simulation were selected for investigation: Production rate and Cumulative production; Pressure and pressure drop (Field pressure, Bottom Hole Pressure, BHP); Oil saturation (Drainage volume per unit area); Minimum stress and Minimum principal effective stress vector; and, Maximum stress and Maximum principal effective stress vector. The usual output focus for interference and well spacing studies is production rate, cumulative production, and pressure. CMG-GEM outputs drainage volume and principal stress vector, which is beneficial for an improved optimization study in unconventional reservoir.

Results and discussions
Oil production from tight sands has a characteristic steep drop in production rate after fluid in the immediate fracture is produced. Applying a minimum bottom-hole pressure of 3200 psi and maximum stock tank oil rate of 2,000 BOPD at surface, Figures 5 and 6 show oil and gas rate respectively from horizontal well with 22 fractures at 150 fracture spacing.  Spikes seen in the rate are indicative of variable non-uniform production rates. At early times, fluid collected in the pores of the fractures within the stimulated reservoir volume (SRV) are produced in large volume. This is represented by high rates recorded until the stimulated volume is depleted. Further on, production depends on fluid transmissibility from matrix into fractures. More so, subsequent variable rates are due to the number and conductivities of secondary fractures present and generated in the process. Comparatively, a similar output for history matched model is shown in Figure 7.

Case study
To properly evaluate well placement and interference, three key physical properties are be examined for impact of the reservoir on regions away from the well location. This study examines drainage volume, pressure transient, and effective minimum principal stress.

Drainage volume
Single fracture -vertical well: The 2D simulation model for this test consists of a vertical well with planar fractures propagating outwards and penetrating total reservoir thickness. The drainage for a simulated productive period of 30 years is shown in Figure 8a assuming that the well is in a virgin pressure region. Drainage coverage stretches beyond the stimulated reservoir volume (SRV) by about 75ft over its fracture half length, and more than 150ft from the fracture face (west -east direction). Maximum volumetric output is possible without any form of nearby connection or interference. This is not expected in multi-fracture system. Instead, reduced fracture performance and far reaching effect into the reservoir. Multiple fracture -horizontal well: In Figure 8b, a single horizontal well has more drainage coverage than the vertical well. This is because of its lateral length and number of fractures. Drainage advancement is consequent on merging of single regions by two or more fractures.
The combined influence results in extended reservoir reach beyond fracture half length (a total distance X DL );although, a little longer than seen in the vertical well. But Drainage length, at both ends of the well, opposite a fracture face, is the same as in single fracture vertical well. Basically, number of fractures is inversely proportional to fracture spacing. Table 3 shows drainage length gain beyond SRV for corresponding fracture spacing. Considers interference in pressure, as well as drainage area and in-situ principal stresses.
Spacing of interfering infill well is independent of production effect of parent well.
Infill well spacing and hydraulic fracture orientation depends on the duration of production of parent well as it affects in-situ stress orientation

Infill wells always have same orientation with parent well
This is achievable only if all the wells are drilled and stimulated at same time.
An infill well cannot have same orientation with parent well beyond 450ft well spacing after years of production.
Volumetric gain from the hydraulically fractured horizontal well is about 1400% more than that of the vertical well. This gain is subject to incremental number of fractures resulting from reduced fracture spacing shown in Figure 9. Overall production performance is good, but the volumetric performance of each fracture is less than recorded for vertical well.

Pressure transient
Pressure drop at the well throughout the field life is up to 30% and 54% for the vertical and horizontal well respectively. But pressure at the outer reservoir remains unchanged at 8000psig. Distance between these extreme pressures make up distance travelled by pressure transient disturbance (X PT ) and connotes limit of interference in the adjacent well. In the single fracture vertical well, X PT is 470ft.But in a multiple hydraulic fracture system, transient travel distance varies with fracture spacing as shown in Figure 10.
The lowest fracture spacing, 150ft, accounts for maximum pressure transient travel of 640ft. A distance marking minimum well spacing for limited interference effect. This distance, X PT, is higher than the corresponding drainage length (X DL ) for same fracture spacing, but less than 700ft well spacing commonly employed in practice. At maximum drainage length, 340ft, corresponding pressure is 5614psi and is higher than BHP. The pressure gradient can still support production. Hence, drainage distance can be used as the yard stick for placing adjacent well for optimum reservoir drainage, but with pressure interference. In shared relation, both X DL, and X PT are affected by fracture half length, and serve as a key determinant of well spacing. However, minimum stress is another criterion that can give a more meaningful conclusion since it affects fracture orientation and geometry.

In-situ stress distribution
Minimum principal stress and its orientation play critical role in hydraulic fracture stimulation. At initial condition, principal stress is unidirectional, and its vector is shown in Figure11(a,b)with horizontal well running from left to right and hydraulic fractures propagating perpendicular to the minimum principal stress in the north-south direction. Altered by mechanical (hydraulic fracture) and poro-elastic effects (pore pressure changes), the stress vector changes direction from 0 o to 90 o as its magnitude changes with defined gradient.
Separated by 'axes-of-rotation' running horizontally through fracture/ well centre and perpendicular through the centre of fracture/mid-length of horizontal well, two kinds of stresses can be defined. Backslash (\), where North-West and South-East stresses rotate anticlockwise, and Forward slash (/), where North-East and South-West stress rotate clock wise. At maximum rotation, minimum effective principal stress changes orientation to vertical direction pointing southward. Around the fractures, opposite the fracture faces, stress rotates in response to fracture opening and they rotate more during production.
Due to fluid continuity, sustained production alters pressure further away from the well, consequently, causing stress rotation between the wells and a maximum distance into the reservoir (far-field). Figure  12 shows stress response to vertical and horizontal wells. In both wells, only stress vectors acting along axes-of-rotation are aligned in extreme rotation (0 o and 90 o ) for all distance from the well into farfield reservoir while others stress vectors take any acute inclination. Furthermore, a different peculiar pattern was observed in horizontal well. Stress acting in the direction of fracture face, over its full length, align in horizontal direction (0 o ) after a short period of production. The complex orientation of stresses seen around fracture at start of production all tend toward horizontal alignment over the observed time. The initial distortion in stress could be attributed to transient response of stress to transition from static to dynamic flow state, whereas the stress alignment afterwards could be attributed to a steady flow condition. This horizontal tendency in stress vectors, spreads beyond the fracture length (or SRV) into external reservoir at distance of stress orthogonality Ⱶ σ , measured from horizontal axis of the well. Beyond this distance, vertical acting stress vectors become dominant. At first, vertical oriented stress vectors act along the vertical axis of horizontal well, subsequently, this pattern spreads to the left and to the right simultaneously as production continues. At very long time, the whole stress vectors align in such a pattern that makes the whole horizontal well's length act like a single fracture as seen in vertical well. In both wells, orthogonal stress vectors along the axis of rotation assume their extreme anglesas production commences, while surrounding vectors take form later.

Stress orthogonality
In vertical well, two distances of stress orthogonality (Ⱶ σ ) can be seen; one at the centre of the well (Ⱶ σ1 ) with zero magnitude and the second (Ⱶ σ2 ) at 750ft. In horizontal well, Ⱶ σ1 measured from centre of the well varies with fracture spacing by observation. Figure 13 shows this relationship. This behaviour is most likely the effect of combination of stress characteristics of individual fractures. The early time post production horizontally aligned stresses take shape within Ⱶ σ1. However, Ⱶ σ2 could not be determined due to the size of the model. Even if made visible by increasing the model, the magnitude is large enough to validate well placement conclusions as reported in this study.

Figure 13 Distance of stress orthonality versus fractue spacing
Following the previous discussion, the model with lowest fracture spacing (150ft) gave highest production and corresponds to Ⱶ σ1 of 450ft. Though, at very late time, this distance increases to 600ft. In summary, by combining and comparing the effect of drainage, pressure and stress on near and far field reservoir, X DL is less than Ⱶ σ , and both are less than X PT in the relation: This correlation constitutes recommended region of placement of in fill well following period of production from parent wells. In very low permeability reservoir, it is very tempting to assume placement of infill well outside the stimulated drainage region X DL . The consequence would be poor hydraulic fracture job with high pressure interference. However, placing the well beyond Ⱶ σ2 could be the better guess. This is because it leaves pressure interference effect to latter time when a good portion of the reservoir must have been produced.

Infill model: new well placement
To optimise the suggested region of well placement for minimal interference, the first infill model was a vertical well located at 450ft next to the parent horizontal well with 22 hydraulic fractures at 150ft fracture spacing. Hydraulic fractures of both wells are parallel. The resulting drainage, stress, and pressure profile are shown in Figures 14  and 15. Hypothetically, the vertical well can be placed as such and yet have good stimulation, although complex northward effect of fracture hit is expected to take prominence. Nonetheless, hydrocarbon output would determine the effectiveness of the idea.

Drainage profile
Overall, there seems to be a merger in drainage volume per unit area as shown in Figure 14. Drainage length in the horizontal well closely matches its single state, but that of the vertical well shows less drainage compared to its standalone state (Figure 8a). The effect is evident in cumulative oil production by a reduction of 2%. This is a result of drainage area interference which is consequent on created fracture network and fracture hit. The region of intersection of both fractures is highly depleted, probably the positive effect of fracture hit. The overall output reduction could be because of direction of flow of the reservoir fluids; away from the vertical well and towards the horizontal well.

Stress and field pressure
The merger effect is demonstrated in the distribution of minimum effective principal stress vector as in Figure 15. The Far field stress vectors are in conformity with standalone horizontal well. Stress vector around the vertical well has fair semblance with its standalone profile. The combined stress vector act like a single horizontal well. However, difference lies in the distance of orthogonality. Asymmetric in the horizontal well, Ⱶ σ1 increased in the north by 150ft but remained zero in the vertical well. Because of the location of both wells, their vertical axis aligns into a single vertical axis. It is expected that if the vertical well is relocated further left or right of its current position, the overall effect will be a complete spread of vertical orientation of minimum effective principal stress over the length of the horizontal well above Ⱶ σ1 .

Reverse infill model
In this model, the infill vertical well was located within an area where parent horizontal well has influence. The objective was to locate the infill well within a region of horizontal oriented minimum effective principal stress so that vertical plan arhydraulic fractures can be achieved. The consequence was interference resulting in drainage area and pressure. In reverse, horizontal infill well located close to a vertical well within same spacing is possible. But absolute transverse fractures could not be achieved on its north side due to altered stress state around the vertical well. Non-planar and complex fractures resulted. Overall drainage is not expected to be affected and stress response to the horizontal infill well still dominates at later time. To avoid early conclusion on this interference effect, this study examined other possible well spacing, where optimum fracture propagation could be achieved.

Sensitivity
To evaluate interference against well spacing within permissible range expressed by equation (18), the distance between 450ft and 750ft was demarcated for parallel vertical fractures, and 300ft to 750ft for horizontal fractures in the vertical infill well.

Horizontal and vertical wells -parallel fracture
Interference is a function of well spacing. At wider well spacing, drainage area interference is absent. Closer well spacing results in increased interference. In Figure 16a, drainage section is independent for each well. But, in b and c convergence builds due to some central tendency harbouring potential force. Pressure, which travels farthest at the boundary front, is the potential causing this interference. Volumetric implication is seen in Figure 17.

Horizontal and vertical wells -perpendicular fracture
Like parallel vertical fracture model, oil drainage has nodrainage area interference at large spacingat 600ft and beyond. However, below this, drainage area convergence exists. BHP at the horizontal well was constant for most times but decreases at the vertical well. Also, productivity shows inverse proportionality to well spacing, and the highest output was recorded for the farthest well spacing.

Inference: well spacing -hydraulic fracture orientation -Stress orientation
For both parallel and perpendicular fractures of infill vertical well examined, pressure interference was present in all cases. Larger well spacing resulted in minimum drainage interference and increased production. At such distance, only one of both fracture orientations can be applicable. Usually, infill wells are added at early to mid-period of the field life; 5 years and over. In the model demonstration in a previous section, stress reorientation in horizontal parent well was examined at five years post production in parent well, and minimum distance of orthogonality was approximately 450ft from the horizontal well axis. Beyond this distance into the reservoir, minimum principal effective stress vector acts vertical. Consequently, only horizontal oriented fracture will propagate perpendicular to the minimum stress direction. Although vertical oriented fracture of infill well at farther spacing is not feasible due to vertical stress orientation. By comparison, at same spacing of 750ft, horizontal oriented fracture gave higher oil production as a justification for the stress direction in the location.

Validation of the model: comparative approach
Notable well spacing applied in unconventional reservoir models give various values ranging from 330ft to 1320ft, with a mode spacing value of 700ft applied in field of practice in Eagle Ford shale. Metrics for the varying distances have been on pressure interference and cumulative production only. In this study, a combination of stress distribution and hydrocarbon drainage gives spacing range of 450ft and 750ft with 600ft as minimum recommended. The premises of this inference are in concordance with existing models, whose findings show that pressure transient is a determinant of well spacing and travels beyond SRV into the reservoir for single well. Also, well spacing affects interference and cumulative production inversely.
However, this study offers new insights, particularly in the orientation of infill well and its hydraulic fracture. Table 4 bears these findings.

Conclusion
Hydrocarbon production in shale reservoir was modelled using coupled geomechanics and fluid flow model. Pressure, drainage volume and stress with respect to time are the main variables considered in this work. Pressure interference models in horizontal and vertical wells were developed as a tool to study well spacing. Using a commercial numerical simulator, drainage volume of a well limited by low permeability formation was delineated, and stress vector profile over the reservoir was obtained. These, including pressure, were used in the well-to-well interference study to examine well spacing and orientation. The following conclusions were drawn: • Interference effect resulted in reduced cumulative production.
Reduction of interference effect is dependent on well spacing; and well spacing is dependent on drainage area of a well, pressure, and distance of stress orthogonality.
• The drainage area of a well stretches beyond stimulated reservoir volume into external reservoir by about 750ft.Locating infill well with drainage area next to a parent well (Zipper pattern) is not economically profitable.
• In-situ stress distribution is uniformly oriented initially but rotates with poroelastic changes in the formation and has tendency to achieve extreme rotation at 90 o . Thus, infill wells and their hydraulic fracture cannot have same orientation with parent well.
• Two distances of stress orthogonality (Ⱶ σ ) exists. One toward the well(Ⱶ σ1 ), and the other, at further distance into the formation (Ⱶ σ2 ). In single fracture (vertical well), Ⱶ σ1 is 0ft, and Ⱶ σ2 is 750ft.
In multiple fractures (horizontal well), Ⱶ σ1 is 450ft, and Ⱶ σ2 is unknown (beyond know wides well spacing).
• Separation between Ⱶ σ andⱵ σ2 indicates regions where infill well and its hydraulic fractures must be propagated perpendicular to the parent well. Outside this separation, wells and hydraulic fractures parallel to parent well can be achieved.
• Fracture spacing lower than 230ft should be implemented in very tight reservoir because of increase production recorded. Texas Two Step (Alternate Fracturing) pattern can be employed to achieve this.
• Optimum well spacing for minimum inter-well interference recorded is 600ft • If horizontal infill-well is considered after productive years, it should be drilled perpendicular to the parent well. This is valid for well spacing of 600ft and above.
• For irregular reservoir geometry, vertical wells can be drilled, for improved recovery, in locations where horizontal well is not a good fit. The minimum spacing is valid for both horizontal and vertical wells.

Appendix B
The proposed correlation for multi-rate flow and normalised pressure distribution in the reservoir system is given as follows: Pressure drop at observation well: The corresponding discretised pressure drop is given as: The dimensionless pressure drop for horizontal well in anisotropic system is: ( ) where the pressure drop is approximated as; where,