Next Article in Journal
Bio-Inspired Optimization Algorithms Applied to the GAPID Control of a Buck Converter
Previous Article in Journal
Heterogeneous Effect of “Eco-Friendly” Dwellings on Transaction Prices in Real Estate Market in Portugal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Fracture Propagation in a Dual-Porosity System: Pseudo-3D-Carter-Dual-Porosity Model

Petroleum Engineering Department, Khalifa University of Science and Technology, Abu Dhabi P.O. Box 127788, United Arab Emirates
*
Author to whom correspondence should be addressed.
Energies 2022, 15(18), 6779; https://doi.org/10.3390/en15186779
Submission received: 21 August 2022 / Revised: 11 September 2022 / Accepted: 12 September 2022 / Published: 16 September 2022
(This article belongs to the Section H1: Petroleum Engineering)

Abstract

:
Despite the significant advancements in geomodelling techniques over the past few decades, it is still quite challenging to obtain accurate assessments of hydraulic fracture propagation. This work investigates the effect of fluid leak-off in a dual-porosity system on the hydraulic fracture propagation geometry, which, in turn, affects hydrocarbon recovery from tight and unconventional reservoirs. Fracture propagation within tight reservoirs was analyzed using the Pseudo Three-Dimensional-Carter II model for single- (P3D-C) and dual-porosity systems (P3D-C-DP). Previous studies have accounted for leak-off in single-porosity models; however, studies within dual-porosity systems are still quite limited. We present a novel approach to coupling fluid leak-off in a dual-porosity system along with a fracture-height growth mechanism. Our findings provide important insights into the complexities within hydraulic fracturing treatment design using our new and pragmatic modeling approach. The simulation results illustrate that fluid leak-off in dual-porosity systems contributes to a confined fracture half-length ( x f ) , that is 31% smaller using the P3D-C-DP model as opposed to the single-porosity model (P3D-C). As for the fracture height growth ( h f ) , the P3D-C-DP model resulted in a 40% shorter fracture height compared to the single-porosity model.

1. Introduction

Hydraulic fracturing is considered the most successful well stimulation technique [1]. The first fracturing experiment took place in 1947. In 1950, hydraulic fracturing was deemed commercial. Furthermore, 2.5 million fracs have been estimated worldwide by the Society of Petroleum Engineers, and more than 1 million of these fracs have taken place in the USA [2]. Fracture propagation predictions are non-trivial and require a careful and well-thought-out plan of the hydraulic fracturing design elements. It is essential to have ample information on the reservoir rock properties, fracturing fluids, and the magnitude and direction of the in situ stresses in order to reach an accurate prediction of fracture propagation in terms of the opening width, length, and height [3]. Furthermore, a set of parameters for the hydraulic fracturing treatment, such as time, injection rate, proppant type and concentration, and fracturing fluid viscosity with power law parameters, must be considered at the design stage [4].
Tight reservoirs are challenging in terms of their heterogeneity and extreme low permeabilities. In 1992, Rahim and Holditch [5] emphasized the importance of using 3D hydraulic fracture propagation models for obtaining accurate results; however, most of the fracture treatments in the past were designed using 2D fracture propagation models, especially for optimization processes [6]. Several introductory and key papers were published between late 1950s and early 1970s that set the foundational development for hydraulic fracture modeling approaches. In 1957, Carter [7] developed a model based on a volumetric balance, which assumes a constant fracture width. Then, in 1961 and 1972, the Perkins–Kern–Nordgren (PKN) fracture geometry model was developed by combining the fracture width model of Perkins and Kern [8] with the closed-form equations that account for fluid leak-off and storage [9]. In addition, the P3D model for fracture height propagation was developed for two categories; the first was developed for a symmetric three-layer case where the stress barrier contrasts above and below the pay-zone are assumed to be constant [10], while the second one was established in 1987 by Fung et al. [11], where they derived a generic solution for the P3D model for non-symmetric multilayer cases. In 2010, Rahman and Rahman [12] coupled fluid leak-off equations with the P3D model formulation. In 2018, Wang et al. [13] extended the application of the PKN-C model by coupling fluid leak-off equations with dual-porosity systems, where an equivalent Carter leak-off coefficient was introduced by coupling the matrix and fracture leak-off coefficients. In this work, we build a P3D-C-DP model that involves computing the equivalent leak-off coefficient for the matrix and the fracture systems.
Hydraulic fracturing design requires a full treatment analysis to carefully predict the fracture geometry parameters such as the length, height, and width of the fracture. Numerous analytical models developed over the years have provided several important aspects of hydraulic fracturing processes. Paluszny et al. [14] showed how these models are very useful in terms of highlighting the important factors influencing the growth of the hydraulic fracture. However, they also argued that such models exhibit a limitation when accounting for multiple factors regarding the matrix’s permeability and thermal effects. Wang [15] added that the existence of pre-existing natural fractures may influence hydraulic fracture propagation and the related flow capacity depending on their density, size, direction, and cementation as well as the area where the hydraulic fracturing is carried out. A number of factors, such as the in situ stresses, treatment parameters (such as injection rate and fluid viscosity), and the geometric and mechanical features of the natural fractures, all have the potential to have a major impact on the hydraulic fracture and natural fracture interactions [16,17,18,19]. A recent study by Oyarhossein and Dusseault [20] presented how geomechanical properties and injection parameters affect the fracture height and aperture, while taking into account all technical and environmental constraints. Interesting findings from a pulse-decay experiment performed by Iia and Xian [21] on a fractured shale core and a history matching of the matrix-fracture system demonstrate that, for matrix-only cores, the permeability estimated from the pulse-decay experiment precisely follows Darcy’s law, which varies proportionately with the fluid viscosity. However, the flow in the fracture and the matrix interacts with one another in the fracture–matrix system. Jia et al.’s [22] investigation into the CO2 storage capacity has shed more light into the flow behavior in tight formations. In this study, the fractures serve as the primary flow conduit for fluid transport in a tight porous medium, which has been stimulated using a dual porosity and permeability model.
Limitations of analytical solutions in terms of handling geometries and heterogeneities can be addressed using robust numerical approaches. Kim and Um [23] described the various analytical and numerical approaches used in hydraulic fracturing propagation. Analytical hydraulic fracturing models are classified into two main categories with regards to their dimensional form: 2D and 3D models. These models include the 2D PKN and KGD (Khristianovic and Zeltova-Geertsma and de Klerk) models, where different assumptions are applied for each model. For 3D, the P3D model is mostly used. This model was introduced back in the 1980s, and it is considered an extension to the classical 2D PKN model. The 2D PKN model assumes a fixed fracture height, and the fracture half-length is a variable that propagates with time to capture the extended fracture length ( x f h f ). In contrast to the PKN model, the P3D model allows the fracture height to vary and grow vertically into adjacent layers, and it is not limited to the thickness of the reservoir. It is important to note that both PKN and P3D models are used when the fracture height remains relatively small compared to the fracture half-length, which are incorporated in modeling the production of oil and gas from a hydraulically fractured well [24]. Figure 1 depicts the conventional P3D fracture propagation model, and the symbols are described in the relevant equations in Section 2.

2. Mathematical Formulation

2.1. Fracture Height Growth at the Wellbore

The P3D-C model allows for a varying fracture height, and the letter “C” in the name denotes the Carter Equation II as a solution for the material balance [12]. Moreover, it accounts for the fluid leak-off test as a result of the fracture propagation. The term pseudo refers to the model’s notion of not accounting for fracture geometry change in a fully three-dimensional space, while the P3D model alters the two-dimensional PKN model by including the height variation concept [25]. In other words, the height is considered to change over time and is computed as a function of the lateral position at each time step. The expansion of the aforementioned height is governed by pressure, which functions solely in response to the lateral position [26]. In this model, the fracture half-length, x f , is fixed with a preassumed target length value, and the fracture height, h f , is allowed to vary and propagate along the fracture length; the shape can be linear or parabolic [12].
Using a non-Newtonian fracturing fluid, the power law parameters that are correlated with the fluid’s viscosity are utilized to derive the fracture geometry parameters from the following mathematical correlations [4]:
n = 0.1756 μ 0.1233 ,
and
k = ( 0.457 μ 0.0159 ) 47.880 ,
where n is a dimensionless exponent, k is the consistency index in (Pa.sn), and μ is the fluid viscosity (Pa.s). The model is initialized with a set of treatment variables, which consist of the injection rate, Newtonian fluid viscosity, target length, and fracture height. The fracture height growth at the wellbore is then computed by neglecting the material property variations in each stress layer in order to reach an equilibrium height [12]. The governing equations required for the P3D model begin with computing the stress intensity factors at the top and bottom layers, K l c t and K l c b (psi. inch ), respectively [27]:
K l c t = 1 π a a a p ( y ) a + y a y d y ,
and
a = h f 2 ,
where K l c t is the stress-intensity factor or fracture toughness at the top of the fracture (psi. inch ), a is the fracture half-height (ft), p ( y ) is the net fracture pressure distribution opening the fracture (psi), and y is the vertical position from the center of the fracture.
A geometry constraint term (in ft) is added for the heights ( b 2 and b 3 ) straddling the payzone height (h):
b 3 = h b 2 .
The fracture height at the wellbore can be solved for by integrating Equation (3) and using a similar equation for the bottom layer giving the following two sets of equations:
π 2 a ( K I c b + K I c t ) = ( σ 2 σ 1 ) sin 1 ( b 2 a ) + ( σ 3 σ 1 ) × sin 1 ( b 3 a ) ( σ 2 + σ 3 2 p w ) π 2 ,
and
π 2 ( K I c b K I c t ) = ( σ 2 σ 1 ) a 2 b 2 2 ( σ 3 σ 1 ) a 2 b 2 2 ,
where p w is the treatment pressure at the wellbore (psi), σ 1 is the minimum horizontal in situ stress (psi), σ 2 is the intermediate principal stress (psi), and σ 3 is the minimum principal stress (psi). The treatment pressure at the wellbore is estimated by rearranging Equation (6). The height growth of the fracture can be computed using an iterative algorithm, which will be described in the following section after coupling the P3D model with the Carter Equation II.

2.2. Pseudo 3D Model Coupled with Carter Equation II

Rahim and Holditch [5] and Rahim [6] emphasized the importance of using 3D hydraulic propagation models. However, the majority of the fracture treatments are designed using 2D fracture propagation models. To be able to determine the fracture width, it is required to have an equilibrium condition between the pressure in the wellbore and the fracture using closed-form equations [12]. It is recommended to use the 2D PKN-C equations to compute for the maximum and average fracture width at the wellbore [25]:
w f = 9.15 ( 1 2 n + 2 ) 3.98 ( n 2 n + 2 ) [ 1 + 2.14 n n ] ( n 2 n + 2 ) k ( 1 2 n + 2 ) ( ( q i 2 ) n h f 1 n t l E ) ( 1 2 n + 2 ) ,
and
w ¯ = π 5 w f ,
where E is the plane strain modulus (Pa), q i is the injection rate (m3/s), and t l is the target length (m). The fracture height growth estimation requires knowing the in situ stresses, fracture toughness above and below the payzone height, and thickness of the layer. In addition, net pressure can be estimated either from previous fracture treatments or from the output obtained from the P3D model [28]. The net fracture pressure is the difference between the fracturing fluid pressure and the closure pressure, and it refers to the driving mechanism behind the fracture growth [29]. The net fracture pressure is obtained from the following equation [26]:
p n , w = E 2 h f w f .
Within the proposed implementation of the P3D-C model in this study, the default fracture height was set to 101 ft. The initial estimate of the fracture half-height at the wellbore can be computed from Equation (4). b 2 can be computed numerically from Equation (7) using a root-finding method, such as the bisection method. In Equation (10), the net pressure was computed with the changes of the fracture height growth. Rahim and Holditch [28] presented typical net pressure values from field experience ranging between 300 to 2000 psi. It is important to note that the fracture height at any point in the fracture acts as a function of the net pressure. Following Rahman and Rahman [12], a user-defined tolerance can be set for the fracture height computations as:
T o l e r a n c e = A b s ( p t r e a t p w ) p t r e a t × 100 ,
where the treatment pressure of the fracturing fluid is computed from the following equation:
P t r e a t = σ 1 + P n e t .  
An iterative scheme is then implemented based on the following incremental change:
Δ h f = 0.1 × h f ( i )
where Δ h f is the deferential height change (ft). In addition, if P w < P t r e a t then the new fracture height, h f ( i + 1 ) , is expected to increase by:
h f ( i + 1 ) = h f ( i ) + Δ h f ,
and if P w > P t r e a t it would decrease by:
h f ( i + 1 ) = h f ( i ) Δ h f .
Assuming a linear height variation along the fracture half-length, the fracture area is calculated as:
A p = x f ( h f + h 2 ) .
For two fracture wings, the propped fracture area, A p , is given by:
A p = x f ( h f + h ) ,
and the total surface area for two fracture wings, A f , would be:
A f = 4 A p .
The Carter Equation II gives a relationship between the fracture-half length, height, average fracture width, and injection time for a given injection rate [4]. This relationship can be extended further by adding the summation of the changing fracture height and the payzone height, leading to the following expression for the fracture half-length:
x f = ( w ¯ + 2 S p ) 4 C l 2 π q i ( h f + h ) [ e x p ( β 2 ) e r f c ( β ) + 2 β π 1 ] ,
and
β = 2 C L π t i w ¯ + 2 S p ,
where S p is the spurt loss (m), C l is the leak-off coefficient (m/ s ), β is the auxiliary variable for the Carter II Equation, t i is the injection time (seconds), and w ¯ is the average fracture width (m). Finally, the fracture volume, V f (m3), is calculated from:
V f = π 2 γ ( h f + h 2 ) x f w f

2.3. Governing Equations for Fluid Leak-Off in a Dual-Porosity System

Fluid leak-off is investigated in this work for a dual-porosity system. In a hydraulic fracturing process, leak-off occurs when the influx of the fracturing fluid propagates to the adjacent natural fracture path layers. Ghaderi et al. [30] stated that the propagating dimension would be affected significantly by the leak-off rate. Here, we compare our P3D-C-DP model in a dual-porosity system against a single-porosity system (P3D-C model) and analyze the results in terms of the induced fracture aperture, width, fracture half-length, and height. Ghaderi et al. [30] further mentioned that the fracture and matrix permeabilities play an important role in linking the leak-off phenomena in a naturally fractured reservoir. There are three main concepts to consider for a hydraulic fracture treatment in naturally fractured reservoirs:
  • fluids penetrating into the matrix;
  • fluids flowing into natural fractures;
  • the transfer of fluids from natural fractures into the matrix.
Wang et al. [13] presented an equivalent Carter leak-off coefficient for a dual medium implemented in the PKN-C model. Following Wang et al. [13], the equivalent Carter leak-off approach was used in this work and implemented in a P3D-C model for dual-porosity systems. The P3D-C model applied in this study accounts for fracturing fluids that exhibit a non-Newtonian behavior. The 1-D single-porosity leak-off perpendicular to the hydraulic fracture was compared to the dual-porosity leak-off. The governing equations follow the same set of equations in Wang et al. [13], with the exception of the extension to P3D-C models. During the first stage of initiating the hydraulic fracture treatment, it is assumed that the natural fractures remain sealed and unchanged. However, when the treatment pressure increases, it would lead to dilating and reactivating the natural fractures. When the natural fractures open, it is required to calculate the critical fissure opening pressure, p c , given by:
p c = 1 α f [ ( 2 v p + σ m a x ) sin 2 θ + ( p + σ m i n ) cos 2 θ ] ,
where p is the net pressure in the hydraulic fracture, and θ is the angle between the hydraulic fracture and the natural fracture. Furthermore, the value of the critical fissure opening pressure and the treatment pressure determines which equation to use for the hydraulic aperture of the natural fracture (m),   b , using the following criterion:
b = { b 0 p f < p c b 0 e x p ( χ σ n e ) p f > p c ,
where b 0 is the initial natural fracture aperture, χ is a relatively insignificant constant that describes the compliance of a natural fracture with regard to a change in pressure, and σ n e is the effective normal stress. The dynamic porosity,   ϕ f , and permeability, k f , of the natural fracture system are calculated, respectively, from [31,32]:
ϕ f = ϕ f 0 b b 0 ,
and
k f = k f 0 ( b b 0 ) 3 [ cos 2 θ sin θ cos θ sin θ cos θ sin 2 θ ] .
Assuming a constant leak-off coefficient, the fluid leak-off velocity directly penetrating into the matrix can be written as follows:
u 1 ( x , t ) = 2 C l m t τ ( x ) ,
and
C l m = Δ p k m ϕ m μ π K w ,
where C l m is the matrix leak-off coefficient (m/ s e c ), t is the total injection time ( s e c ), τ ( x ) is the arrival time ( s e c ) of the fracture tip at location x, k m is the matrix permeability (m2), ϕ m is the matrix porosity, K w is the fluid bulk modulus (Pa), and Δ p is the pressure differential between the hydraulic fracture and the formation (Pa). The fluid bulk modulus is computed from the following mathematical expression:
K w = E 3 ( 1 2 v ) .
The differential pressure, Δ p , is computed as:
Δ p = p f l u i d ( x , t ) p 0
where p f l u i d ( x , t ) is the treatment pressure (Pa) and p 0 is the initial formation pressure (Pa). In addition, the dimensionless time scaling is given by [13]:
τ = t t * ,
and the characteristic quantity, t * is computed as:
t * = π 2 H 6 p * 5 4 E 4 μ π Q o 2   ,
where p * is a characteristic net pressure (Pa). Assuming that the fluid transfer from the natural fractures to the matrix is ignored, then the fluid leak-off velocity flowing into the natural fractures can be approximated by [13]:
u 2 ( x , t ) = 2 C l f t τ ( x ) ,
and
C l f = Δ p k f ( p ) ϕ f μ π K w .
The fracture leak-off coefficient, C l f , is considered a varying parameter due to the fracture permeability, k f ( p ) , being a function of net pressure, as shown in Equation (25). In addition, the net pressure varies slowly with time, t , and location, x . Therefore, the use of an average net pressure, p ¯ , is adopted for simplicity [13]. Furthermore, when the fluid transfer from the natural fractures to the matrix is not ignored, then the fluid leak-off velocity can be approximated by:
u 3 ( x , t ) = 8 v t r C l f t τ ( x ) b ( p ¯ ) ,
where the fluid leak-off velocity from the natural fractures to the matrix, v t r (m/s), using the fracture spacing, s (m), is given by:
v t r = k m μ Δ p s / 2   .
The total fluid leak-off velocity introduced by Wang et al. [13] can then be approximated by:
u ( x , t ) = u 1 ( x , t ) + u 2 ( x , t ) + u 3 ( x , t ) = 2 C l e t τ ( x ) ,
where the equivalent Carter leak-off coefficient, C l e , is expressed by:
C l e = C l m + C l f + 8 k m Δ p C l f [ t τ ( x ) ] s μ b ( p ¯ ) .
The Carter Equation II, Equation (19), with the equivalent leak-off coefficient from Equation (37) predicts the fracture half-length as the fracture propagates.

3. Numerical Results

3.1. Model Validation

Wang et al. [13] presented an equivalent Carter leak-off coefficient into a pressure-sensitive dual medium. The model was used to study the interaction between the natural fractures and the induced fracture, and it was treated as a dual-porosity medium. The 1-D single-porosity leak-off perpendicular to the hydraulic fracture was compared with the dual-porosity leak-off. In our work, we have developed two codes: (1) PKN-C-DP and (2) P3D-C-DP. Our work extends Wang et al.’s [13] work to the P3D model. The PKN-C-DP code was validated with the same input parameters used in Wang et al. [13]. The algorithm for the P3D-C-DP model is shown in Figure 2 and Table 1 presents the input parameters used for the model validation [12].
The pressure difference between the hydraulic fracture and the formation is needed to compute the fluid leak-off coefficient in both the matrix and the fracture. Using a non-dimensional time, Wang et al. [13] ensured that the simulation progressed from a small-time to a large-time solution. Comparing the small-time asymptote at the beginning of the simulation, Wang et al. [13] found the storage dominated regime occurs at the beginning of the simulation. As a result, the fracture’s development stretches to the end of the leak-off dominated regime at a large time asymptote. Equation (31) was utilized as a basis for validation, but since the resultant dimensionless time value was larger than 1, it could only be employed in the leak-off-dominated regime. This validation neglected the storage-dominated regime. The model did not analyze fracture geometry parameters after shut-in for a 3000-second injection time. Table 2 shows the results after a 3000-second injection time.
Figure 3a,b provide a comparison of model validation in the dual-porosity medium with the Wang et al. [13] model. Figure 3a showcases the average width (mm) changes as a function of total pumping time. The validation model and the observed average width changes at the inlet are very similar. Figure 3a also represents the fugitive fluid leak-off from the fracture into the surrounding formation, which eventually leads to a reduction in the maximum width. This is due to the interaction between the induced hydraulic fracture and the natural fracture. The average fracture width at the middle region does not exactly match Wang et al.’s [13] model, which could be due to a variety of factors including: differences in the way the time function is constructed, use of non-Newtonian equations, or variation of the target length value used. Leak-off in the dual-porosity medium is expected to impact fracture geometry during treatment by reducing the fracture length. Figure 3b shows a validated example of fracture length (m), and it closely matches Wang et al.’s [13] model.

3.2. Numerical Examples

Table 3 presents a reservoir well data for a single-porosity system. Here, the P3D-C and PKN-C models were used to compare single and dual-porosity systems (PKN-C-DP and P3D-C-DP) utilizing previously published numerical example data [12].
Both single- and dual-porosity system models used the same initial set of input parameters. Table 4 presents additional input parameters for computing the equivalent Carter leak-off coefficient [13].

3.2.1. Fracture Geometry Propagation

The fracture geometry parameters include the fracture half-length, maximum fracture width at the wellbore, and average fracture width. Figure 4a,b present the fracture half-length evolution for the single- and dual-porosity systems using PKN-C and P3D-C, respectively. The fracture half-length, in Figure 4a, for the dual-porosity system, is approximately 653 ft less (21% shorter) than it in the single-porosity system. This is expected since the fluid leak-off is much larger in the dual medium. The fracture half-length predicted by the P3D-C model for a single-porosity system in Figure 4b gives a larger fracture half-length when compared to Figure 4a. However, the fracture half-length in the P3D-C-DP is 150 ft less than that of the PKN-C-DP model and exhibits a 31% reduction when compared to the P3D-C model. This is because the fracture half-length and fracture height obtained from the P3D-C and P3D-C-DP models are inversely proportional to each other [12]. Both Figure 4a,b depict a common significant decrease when the dual-porosity-medium system is introduced. This is due to the fracture permeability, where it exerts a significant amount of fluid pressure, which in turn results in a much reduced length of the propagating fracture.
Figure 5a,b indicate that the PKN-C model in the single-porosity system predicts a slightly greater average fracture width and maximum fracture width at the wellbore. In contrast, lower values are obtained from the PKN-C-DP model. This is because a higher leak-off coefficient results in a lower build-up pressure, leading to a narrower fracture width [33]. Similarly, when comparing the single-porosity system to the dual-porosity system, we find that more fluid leak-off occurs in the dual porosity medium resulting in a smaller maximum and average fracture width. Thus, the fracture for the dual-porosity leak-off shrinks and becomes more confined in size compared to the fracture for the single-porosity leak-off.
Figure 6a,b depict the average fracture width and the maximum fracture width at the wellbore for the single and dual mediums, respectively. The average fracture width in the single porosity medium is 12% greater than that in the dual medium. Similarly, the maximum fracture width at the wellbore in the single medium is 11% greater than in the dual medium. In the dual-medium setting, the fluid leak-off has a particularly strong impact and as fluid leak-off increases, the fracture width decreases. The average and maximum values for fracture width in the P3D-C and P3D-C-DP models are similar to the findings obtained from the PKN-C and PKN-DP models. It is important to note that the intersection angle between the natural fractures and the hydraulic fracture affects the fracture length and width; a larger angle causes the fracture width and length to be reduced [13]. In this work, however, the angle was kept constant at π / 2 .

3.2.2. Fracture Height Growth

The fracture height in the single-porosity system in Figure 7 tends to grow linearly with time as it reaches 188 ft. However, there is a significant contrast in the P3D-C dual-porosity model, where the fracture height is 40% shorter compared to the single-porosity system. In contrast to the P3D-C model, where the fracture height propagation grows more slowly due to the constant leak-off coefficient used in the single-porosity system, the P3D-C-DP coupled with fluid leak-off equation for the dual-porosity system displays a steady increase in terms of the fracture height propagation. However, in the dual-porosity system, when the equivalent Carter leak-off equation is introduced, both the matrix and fracture leak-off coefficients influence fracture height growth and since the equivalent Carter leak-off velocity is substantial, the resulting fracture height growth is limited. This is because a higher fluid leak-off contributes to more confined fracture height and limits its growth.

3.2.3. Net Fracture Pressure

The net fracture pressure is required to provide the energy needed for fracture propagation. Figure 8a,b present the net fracture pressure changes for the single and dual mediums using the PKN-C and P3D-C models, respectively. As a consequence of fugitive fluid leak-off into the surrounding formation, the net fracture pressure in the dual media decreases, which, in turn, decreases the maximum fracture width [13]. This is because the maximal fracture width at the wellbore is directly proportional to the net fracture pressure. When compared to the 2D PKN model, the P3D model’s net fracture pressure values are lower, since they are based on the average fracture height. In all models, the net fracture pressure in the single-porosity system is greater than that in the dual-porosity medium.

3.3. Sensitivity Analysis

A series of sensitivity analyses were performed to examine the effect of varying the fracture height and half-length. Figure 9a,b illustrate the effect of changing the fracture height on the fracture half-length in the PKN-C and PKN-C-DP models, respectively. Naturally, a larger fracture height contributes to a reduced fracture half-length.
Figure 10a,b present the impact of changing the fracture half-length on the fracture height growth; a higher fracture height growth is achieved in the single porosity P3D-C model in Figure 10a, whereas the fracture height growth is more confined in the dual-porosity system due to the fugitive fluid leak-off that limits the fracture height growth and results in a more confined fracture height in Figure 10b.

4. Conclusions

In this work, a semi-analytical model to evaluate and interpret the hydraulic fracture propagation in dual-porosity reservoirs was developed. The model entailed the coupling of the so-called equivalent Carter leak-off coefficient for dual-porosity models in a P3D approach. The specific findings from this work are listed below.
  • The model was successfully validated against a case study conducted on fluid leak-off in a dual-porosity system.
  • Generally, the fracture geometry parameters ( x f , h f , w f , w ¯ ) in the single-porosity models are significantly larger when compared to the dual-porosity models, and this is primarily due to the effect of fluid leak-off into the natural fractures.
  • The equivalent Carter leak-off coefficient in a dual-porosity system is directly related to the permeability of the natural fractures; higher permeability values accelerate the fluid leak-off velocity, which ultimately affects the fracture geometry propagation and limits the fracture growth.
  • For a given injection time in the PKN-C-DP model, the effect of increasing the fracture height yields shorter fracture half-lengths. Furthermore, the confinement of fracture height in the P3D-C-DP model increases the fracture half-lengths.
  • The simulation findings show that fluid leak-off in a dual-porosity system contributes to a constrained fracture half-length ( x f ), with x f being 21% and 31% smaller in PKN-C-DP and P3D-C-DP models, respectively, compared to the single-porosity model. The P3D-C-DP model shows a 40% reduction in fracture height growth ( h f ) compared to the single-porosity model (P3D-C).

Author Contributions

Conceptualization, M.M.R.; Data curation, F.A.H.; Formal analysis, F.A.H. and M.A.K.; Investigation, A.S., M.A.K., M.M.R. and M.H.; Methodology, F.A.H. and A.S.; Project administration, M.M.R.; Resources, A.S., M.A.K., M.M.R. and M.H.; Supervision, M.A.K. and M.M.R.; Validation, F.A.H.; Writing—original draft, F.A.H.; Writing—review & editing, M.A.K., M.M.R. and M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not Applicable.

Acknowledgments

We would like to thank Khalifa University of Science and Technology for the scholarship provided to the first author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ghalambor, A.; Ali, S.A.; Norman, W.D. Frac Packing Handbook; Society of Petroleum Engineers: Richardson, TX, USA, 2009; 610p. [Google Scholar]
  2. King, G.E. Hydraulic fracturing 101: What every representative, environmentalist, regulator, reporter, investor, university researcher, neighbor and engineer should know about estimating frac risk and improving frac performance in unconventional gas and oil wells. In Proceedings of the SPE Hydraulic Fracturing Technology Conference, The Woodlands, TX, USA, 6–8 February 2012. [Google Scholar] [CrossRef]
  3. Nasirisavadkouhi, A. A Comparison Study of KGD, PKN and a Modified P3D Model. 30 July 2015. Available online: https://www.Researchgate.net/publication/280568251 (accessed on 1 August 2021).
  4. Rahman, M.M. Constrained hydraulic fracture optimization improves recovery from low permeable oil reservoirs. Energy Sources Part A 2008, 30, 536–551. [Google Scholar] [CrossRef]
  5. Rahim, Z.; Holditch, S.A. The effects of mechanical properties and selection of completion interval upon the created and propped fracture dimensions in layered reservoirs. In Proceedings of the SPE 24349, SPE Rocky Mountain Reg. Meet., Denver, CO, USA, 18–21 May 1992. [Google Scholar]
  6. Rahim, Z. Proppant Transport Down a Three-Dimensional Planar Fracture. Ph.D. Thesis, Texas A&M University, College Station, TX, USA, 1988. [Google Scholar]
  7. Carter, R.D. Derivation of the general equation for estimating the extent of the fractured area. In Optimum Fluid Characteristics for Fracture Extension; Appendix I; Drilling Production Practice; Howard, G.C., Fast, C.R., Eds.; American Petroleum Institute: New York, NY, USA, 1957; p. 261. [Google Scholar]
  8. Perkins, T.K.; Kern, L.R. Widths of hydraulic fractures. J. Pet. Technol. 1961, 13, 937–949. [Google Scholar] [CrossRef]
  9. Nordgren, R.P. Propagation of a Vertical Hydraulic Fracture. SPE J. 1972, 12, 306–314. [Google Scholar] [CrossRef]
  10. Simonson, E.R.; Abou-Sayed, A.S.; Clifton, R.J. Containment of Massive Hydraulic Fractures. SPE J. 1978, 18, 27–32. [Google Scholar] [CrossRef]
  11. Fung, R.L.; Vijayakumar, S.; Cormack, D.E. Calculation of vertical fracture containment in layered formations. SPE Form. Eval. 1987, 2, 518–522. [Google Scholar] [CrossRef]
  12. Rahman, M.M.; Rahman, M.K. A review of hydraulic fracture models and development of an improved pseudo-3D model for stimulating tight oil/gas sand. Energy Sources Part A Recovery Util. Environ. Eff. 2010, 32, 1416–1436. [Google Scholar] [CrossRef]
  13. Wang, J.; Elsworth, D.; Denison, M.K. Hydraulic fracturing with leakoff in a pressure-sensitive dual porosity medium. Int. J. Rock Mech. Min. Sci. 2018, 107, 55–68. [Google Scholar] [CrossRef]
  14. Paluszny, A.; Salimzadeh, S.; Zimmerman, R.W. Finite-element modeling of the growth and interaction of hydraulic fractures in poroelastic rock formations. In Hydraulic Fracture Modeling; Gulf Professional Publishing: London, UK, 2018; pp. 1–19. [Google Scholar]
  15. Wang, H. Hydraulic fracture propagation in naturally fractured reservoirs: Complex fracture or fracture networks. J. Nat. Gas Sci. Eng. 2019, 68, 102911. [Google Scholar] [CrossRef]
  16. Rahman, M.M.; Rahman, S.S. Fully coupled finite-element–based numerical model for investigation of interaction between an induced and a preexisting fracture in naturally fractured poroelastic reservoirs: Fracture diversion, arrest, and breakout. Int. J. Geomech. 2013, 13, 390–401. [Google Scholar] [CrossRef]
  17. Zhao, K.; Jiang, P.; Feng, Y.; Sun, X.; Cheng, L.; Zheng, J. Numerical investigation of hydraulic fracture propagation in naturally fractured reservoirs based on lattice spring model. Geofluids 2020, 2020, 8845990. [Google Scholar] [CrossRef]
  18. Bakhshi, E.; Golsanami, N.; Chen, L. Numerical modeling and lattice method for characterizing hydraulic fracture propagation: A review of the numerical, experimental, and field studies. Arch. Comput. Methods Eng. 2021, 28, 3329–3360. [Google Scholar] [CrossRef]
  19. Qin, M.; Yang, D.; Chen, W. Three-dimensional hydraulic fracturing modeling based on peridynamics. Eng. Anal. Bound. Elem. 2022, 141, 153–166. [Google Scholar] [CrossRef]
  20. Oyarhossein, M.; Dusseault, M.B. Sensitivity Analysis of Factors Affecting Fracture Height and Aperture. Upstream Oil Gas Technol. 2022, 9, 100079. [Google Scholar] [CrossRef]
  21. Jia, B.; Xian, C.G. Permeability measurement of the fracture-matrix system with 3D embedded discrete fracture model. Pet. Sci. 2022, 19, 1757–1765. [Google Scholar] [CrossRef]
  22. Jia, B.; Chen, Z.; Xian, C. Investigations of CO2 storage capacity and flow behavior in shale formation. J. Pet. Sci. Eng. 2022, 208, 109659. [Google Scholar] [CrossRef]
  23. Kim, J.; Um, E.S. A Framework of Integrated Flow–Geomechanics–Geophysics Simulation for Planar Hydraulic Fracture Propagation. In Hydraulic Fracture Modeling; Gulf Professional Publishing: London, UK, 2018; pp. 21–39. [Google Scholar]
  24. Rahman, M.M. Productivity Prediction for Fractured Wells in Tight Sand Gas Reservoirs Accounting for Non-Darcy Effects. In Proceedings of the SPE 115611, SPE Russian Oil & Gas Technical Conference and Exhibition, Moscow, Russia, 28–30 October 2008. [Google Scholar]
  25. Settari, A.; Cleary, M.P. Three-dimensional simulation of hydraulic fracturing. J. Pet. Technol. 1984, 36, 1177–1190. [Google Scholar] [CrossRef]
  26. Valkó, P.; Economides, M.J. Hydraulic Fracture Mechanics; Wiley: Chichester, UK, 1995; Volume 28, 206p. [Google Scholar]
  27. Rice, J.R. Mathematical analysis in the mechanics of fracture. Fract. Adv. Treatise 1968, 2, 191–311. [Google Scholar]
  28. Rahim, Z.; Holditch, S.A. Using a three-dimensional concept in a two-dimensional model to predict accurate hydraulic fracture dimensions. J. Pet. Sci. Eng. 1995, 13, 15–27. [Google Scholar] [CrossRef]
  29. Belyadi, H.; Fathi, E.; Belyadi, F. Hydraulic Fracturing in Unconventional Reservoirs: Theories, Operations, and Economic Analysis; Gulf Professional Publishing: London, UK, 2019. [Google Scholar]
  30. Ghaderi, A.; Taheri-Shakib, J.; Sharifnik, M.A. The effect of natural fracture on the fluid leak-off in hydraulic fracturing treatment. Petroleum 2019, 5, 85–89. [Google Scholar] [CrossRef]
  31. Robertson, E.P.; Christiansen, R.L. A permeability model for coal and other fractured, sorptive-elastic media. SPE J. 2008, 13, 314–324. [Google Scholar] [CrossRef] [Green Version]
  32. Lei, G.; Dong, P.C.; Mo, S.Y.; Yang, S.; Wu, Z.S.; Gai, S.H. Calculation of full permeability tensor for fractured anisotropic media. J. Pet. Explor. Prod. Technol. 2015, 5, 167–176. [Google Scholar] [CrossRef] [Green Version]
  33. Newman, M.S.; Pavloudis, M.; Rahman, M.M. Importance of Fracture Geometry and Conductivity in Improving Efficiency of Acid Fracturing in Carbonates. In Proceedings of the Canadian International Petroleum Conference, Calgary, AB, Canada, 16–18 June 2009. [Google Scholar]
Figure 1. The Pseudo 3D fracture propagation model.
Figure 1. The Pseudo 3D fracture propagation model.
Energies 15 06779 g001
Figure 2. Algorithm for the P3D-C-DP model.
Figure 2. Algorithm for the P3D-C-DP model.
Energies 15 06779 g002
Figure 3. (a) Evolution of the average fracture width versus pumping time for leak-off into a dual-porosity system: Model validation versus Wang et al.’s [13] model; (b) Evolution of the fracture half-length versus pumping time for leak-off into a dual-porosity system: Model validation versus Wang et al.’s [13] model.
Figure 3. (a) Evolution of the average fracture width versus pumping time for leak-off into a dual-porosity system: Model validation versus Wang et al.’s [13] model; (b) Evolution of the fracture half-length versus pumping time for leak-off into a dual-porosity system: Model validation versus Wang et al.’s [13] model.
Energies 15 06779 g003
Figure 4. (a) Fracture half-length versus total injection time for PKN-C and PKN-C-DP models; (b) fracture half-length versus total injection time for P3D-C and P3D-C-DP models.
Figure 4. (a) Fracture half-length versus total injection time for PKN-C and PKN-C-DP models; (b) fracture half-length versus total injection time for P3D-C and P3D-C-DP models.
Energies 15 06779 g004aEnergies 15 06779 g004b
Figure 5. (a) Average fracture width versus total injection time for PKN-C and PKN-C-DP models; (b) maximum fracture width at the wellbore versus total injection time for PKN-C and PKN-C-DP models.
Figure 5. (a) Average fracture width versus total injection time for PKN-C and PKN-C-DP models; (b) maximum fracture width at the wellbore versus total injection time for PKN-C and PKN-C-DP models.
Energies 15 06779 g005aEnergies 15 06779 g005b
Figure 6. (a) Average fracture width versus total injection time for P3D-C and P3D-C-DP models; (b) maximum fracture width at the wellbore versus total injection time for P3D-C and P3D-C-DP models.
Figure 6. (a) Average fracture width versus total injection time for P3D-C and P3D-C-DP models; (b) maximum fracture width at the wellbore versus total injection time for P3D-C and P3D-C-DP models.
Energies 15 06779 g006aEnergies 15 06779 g006b
Figure 7. Fracture height growth versus total injection time for P3D-C and P3D-C-DP models.
Figure 7. Fracture height growth versus total injection time for P3D-C and P3D-C-DP models.
Energies 15 06779 g007
Figure 8. (a) Net fracture pressure versus total injection time for PKN-C and PKN-C-DP models; (b) net fracture pressure versus total injection time for P3D-C and P3D-C-DP models.
Figure 8. (a) Net fracture pressure versus total injection time for PKN-C and PKN-C-DP models; (b) net fracture pressure versus total injection time for P3D-C and P3D-C-DP models.
Energies 15 06779 g008aEnergies 15 06779 g008b
Figure 9. (a) Impact of fracture height variation on fracture half-length versus total injection time for the PKN-C model; (b) impact of fracture height variation on fracture half-length versus total injection time for the PKN-C-DP model.
Figure 9. (a) Impact of fracture height variation on fracture half-length versus total injection time for the PKN-C model; (b) impact of fracture height variation on fracture half-length versus total injection time for the PKN-C-DP model.
Energies 15 06779 g009aEnergies 15 06779 g009b
Figure 10. (a) Impact of fracture half-length variation on fracture height versus total injection time for the P3D-C model; (b) impact of fracture half-length variation on fracture height versus total injection time for the P3D-C-DP model.
Figure 10. (a) Impact of fracture half-length variation on fracture height versus total injection time for the P3D-C model; (b) impact of fracture half-length variation on fracture height versus total injection time for the P3D-C-DP model.
Energies 15 06779 g010aEnergies 15 06779 g010b
Table 1. Input parameters for model validation.
Table 1. Input parameters for model validation.
Input Design ParametersValue
Injection rate, Q i , m3/s 0.004
Fracture height, h f , m 10
Plain strain Young’s modulus, E , GPa 25
Poisson’s ratio, v 0.2
Fluid dynamic viscosity, μ , Pa.s 0.2
Fluid bulk modulus, K w , GPa 2.2
Permeability of the natural fracture system, k f , m2 1.5 × 10 17
Initial permeability of the natural fracture system, k f 0 , m2 3.25 × 10 15
Initial porosity of the natural fracture system, ϕ f 0 0.002
Fracture spacing, s , m 0.2
Biot coefficient, α f 0.8
Permeability of the matrix, k m , m2 1.48 × 10 17
Porosity of the matrix, ϕ m 0.06
Initial formation pore pressure, p 0 , MPa 6
Minimum in situ principal stress, σ m i n , MPa 15
Maximum in situ principal stress, σ m a x , MPa 16
Angle between NFs and HF, θ π / 2
Table 2. Model validation results for the PKN-C-DP code.
Table 2. Model validation results for the PKN-C-DP code.
ParametersPKN-C-DP Code Results
Fracture half-length, x f , m64.75
Average fracture width, w a v g , m0.0025
Maximum fracture width, w m a x , m0.0039
Total pumping time, t e , seconds 3000
Net pressure, P n e t , MPa 4.90
Table 3. Reservoir and well data for a single-porosity system.
Table 3. Reservoir and well data for a single-porosity system.
Input Design ParametersValue
Spurt loss, S p , in 0.0
Fluid viscosity, μ , cP 250
Pumping rate, q i , bbl/min 20
Wellbore radius, r w , ft. 0.35
Fracture height, h f , ft. 100
Target length, t l , ft. 1968.5
Young’s modulus, E , psi 5.075 × 10 6
Fluid leak-off coefficient, C L , ft/min0.5 0.00025
Max. horizontal stress, σ H , psi 7000
Min. horizontal stress, σ 1 , psi 6000
Min. horizontal stress (shale), σ 2 , psi 6700
Min. horizontal stress (shale), σ 3 , psi 7200
Fracture toughness, K I c b , psi-in0.5 1700
Fracture toughness, K I c t , psi-in0.5 1500
Porosity 0.1
Permeability, k , mD 0.1
Spurt loss, S p , in 0.0
Table 4. Additional input parameters for the dual-porosity system.
Table 4. Additional input parameters for the dual-porosity system.
Input Design ParametersValue
Fluid bulk modulus, K w , psi 4.1 × 10 11
Initial permeability of the natural fracture system, k f 0 , mD 3.3
Initial porosity of the natural fracture system, ϕ f 0 0.002
Fracture spacing, s , ft. 0.66
Biot coefficient, α f 0.8
Permeability of the matrix, k m , mD 0.1
Porosity of the matrix, ϕ m 0.1
Angle between NFs and HF, θ π / 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al Hameli, F.; Suboyin, A.; Al Kobaisi, M.; Rahman, M.M.; Haroun, M. Modeling Fracture Propagation in a Dual-Porosity System: Pseudo-3D-Carter-Dual-Porosity Model. Energies 2022, 15, 6779. https://doi.org/10.3390/en15186779

AMA Style

Al Hameli F, Suboyin A, Al Kobaisi M, Rahman MM, Haroun M. Modeling Fracture Propagation in a Dual-Porosity System: Pseudo-3D-Carter-Dual-Porosity Model. Energies. 2022; 15(18):6779. https://doi.org/10.3390/en15186779

Chicago/Turabian Style

Al Hameli, Fatima, Abhijith Suboyin, Mohammed Al Kobaisi, Md Motiur Rahman, and Mohammed Haroun. 2022. "Modeling Fracture Propagation in a Dual-Porosity System: Pseudo-3D-Carter-Dual-Porosity Model" Energies 15, no. 18: 6779. https://doi.org/10.3390/en15186779

APA Style

Al Hameli, F., Suboyin, A., Al Kobaisi, M., Rahman, M. M., & Haroun, M. (2022). Modeling Fracture Propagation in a Dual-Porosity System: Pseudo-3D-Carter-Dual-Porosity Model. Energies, 15(18), 6779. https://doi.org/10.3390/en15186779

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop