Optimal Designs of Well Pattern for the Development of Tight Gas Reservoirs (TGRs)
Jiao Yuwei^{1, *}, Xia Jing^{1}, Yan Jianye^{2}, Xu Daicai^{2}
Article Information
Identifiers and Pagination:
Year: 2016Volume: 9
First Page: 77
Last Page: 90
Publisher Id: TOEFJ-9-77
DOI: 10.2174/1876973X01609010077
Article History:
Received Date: 16/6/2015Revision Received Date: 21/09/2016
Acceptance Date: 26/09/2016
Electronic publication date: 31/10/2016
Collection year: 2016
open-access license: This is an open access article licensed under the terms of the Creative Commons Attribution-Non-Commercial 4.0 International Public License (CC BY-NC 4.0) (https://creativecommons.org/licenses/by-nc/4.0/legalcode), which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.
Abstract
Both horizontal well and fractured-horizontal well have been widely used to develop TGRs. However, the costs of horizontal well and fractured-horizontal well are much higher than the vertical well. Therefore, it is necessary to consider the reservoir conditions for evaluating the potential benefit when choosing well pattern or designing well parameters. In this paper, a simulator of simulating the development of TGRs including slippage flow and stress dependence in matrix, and high-velocity non-Darcy flow and stress effect in hydraulic fractures was firstly developed. Then, it was used to study the development effects of different TGRs using different well patterns and well parameters. Based on the simulation results, the incremental ratio models of horizontal well to vertical well and fractured-horizontal well to horizontal well were achieved. These models can be used to predict the incremental production using horizontal well or fractured-horizontal well. We also obtained the plates of choosing well pattern and designing the corresponding parameters to achieve a good profit in the field.
INTRODUCTION
Tight gas as one of the main unconventional resources is playing more and more important role in recent years throughout the world. In USA, tight gas comprised 26% by the end of 2010, and it will still keep about 21% by the end of 2035 [1] (Fig. 1). 70% of the current worldwide annual production of tight gas is coming from USA, and it will further expand in 2030. However, the production of tight gas outside of North America is an emerging business, and it appears that China and Russia have significant stakes in this domain [2].
In 1970s, the U.S. government decided that the definition of TGR is one in which the expected value of permeability to gas flow would be less than 0.1 md [3]. Holditch [3] reported that the best definition of TGR is “a reservoir that cannot be produced at economic flow rates nor recover economic volumes of natural gas unless the well is stimulated by a large hydraulic fracture treatment or produced by use of a horizontal wellbore or multilateral wellbores”. Some experts thought that the permeability should be between 0.1md and 0.0001md [4, 5].
From the above definitions or discussions, we can see that economic development is the ultimate objective of TGR. Actually, the production of gas is not only a function of permeability, and to define TGR is of limited significance by choosing only a single value of permeability [3]. The reservoir thickness, porosity, permeability and other factors have significant impacts on gas production [2, 6-8]. Moreover, the well pattern also significantly affects the gas recovery, and horizontal and fractured-horizontal wells are common for tight gas development. Comparing with that of the vertical well, the cost of horizontal well is 1.5 to 2.5 times of vertical well, and the overall current commercial success rate of horizontal wells appears to be 65% [9]. Moreover, the cost will sharply increase as the length of horizontal segment increases. Hydraulic fracturing is also a high cost treatment for increasing production rate contrast to horizontal well. The cost of proppant is approximately 90% of the total cost, and if larger and more fractures are created, it may result in a totally uneconomical project [10]. Thus, a project for horizontal well or fractured-horizontal well should consider the reservoir condition to evaluate the potential benefit.
Fig. (1). US natural gas production, 1990-2035, EIA Annual Energy Outlook 2012 Early Release. |
The objective of this paper is to achieve the models for optimizing well pattern and parameters in different TGRs. We first developed a simulator of simulating the development of TGR including slippage flow and stress dependence in matrix, and high-velocity non-Darcy flow and stress effect in hydraulic fractures. Then, we studied the development effects of different well patterns and parameters in TGRs with different petro-physical properties. The incremental production-ratio models of horizontal well to vertical well and fractured-horizontal well to horizontal well were achieved for the TGRs with different features, and two TGRs in China are used to perform the design of well pattern. These models can be effectively used to optimize well patterns and corresponding parameters.
MODELING OF GAS TRANSPORT
Based on the mass balance equation and extended Darcy's law, gas transport in tight reservoirs should be
(1) |
where K_{a} is the apparent permeability; ρ is the gas density; μ is the gas viscosity; Φ is the potential; q is the sink/source term; φ_{e} is the effective porosity; t is time.
Incorporation of Non-darcy Flow Regimes
In TGRs, the pore size is in nano-scale. Non-Darcy flow has a significant impact on the apparent permeability [11]. The apparent permeability of TGRs significantly deviates from the intrinsic permeability because gas flow in such formation occurs by various mechanisms different from the viscous flow regime represented by Darcy’s law [12-14]. Gas slippage effect or Klinkenberg effect is incorporated into this model by modifying the apparent permeability of gas flow [15]:
(2) |
in which,
(3) |
where K_{∞} is the absolute permeability; c is a proportional coefficient; λ is the molecular mean free path; is the mean gas pressure; r is the hydraulic radius.
Incorporation of Stress Dependence
The impacts of effective stress on gas permeability in tight formation were studied by Soeder, Bustin et al. and Wang and Reed [16-18]. The stress-dependent relationship between permeability and effective mean stress for the first type of rock (Type-I rock) called by Raghavan and Chin [19] is in well agreement with the experimental data,
(4) |
Thus, the absolute permeability at the initial scenario should be:
(5) |
where K_{0} is the permeability at effective mean stress σ_{m}=0; η is a parameter; σ_{m} is the mean effective stress, which is mainly related to the vertical overburden load and reservoir pressure by [19]:
(6) |
and the initial mean effective stress σ_{m,i} by:
(7) |
where is the total mean stress pressure, which is a function of σ_{ob,v} and is nearly a constant value for the system studied. Substituting Eqs. (6) and (7) into Eqs. (4) and (5) for σ_{m} and σ_{m,i}, respectively, the absolute permeability considering stress dependence effect becomes [20].
(8) |
Stress Dependence and High-velocity Non-darcy Flow in Hydraulic Fractures
According to the theory of equivalent fracture conductivity, the relationship between the hydraulic fracture and its equivalent refined grid used in simulation is:
(9) |
where K_{hf} is the permeability of the hydraulic fracture; w_{hf} is the fracture aperture; is the equivalent permeability of the hydraulic fracture because it is usually handled by a refined grid; Δx_{hf} is the length of the refined grid along the well.
According to experimental results [21], the relationship between the effective permeability and gas pressure can be written as:
(10) |
Similarly handling Eq. (10) with Eq. (8) yields:
(11) |
where α_{s} is stress sensitivity coefficient for hydraulic fracture; K_{hf}(p_{g}) is the permeability of hydraulic fracture considering stress effect. Applying Eqs. (9) and (11) yields:
(12) |
In order to consider the high-velocity non-Darcy effect, the Forchheimer equation is usually used [22-24]:
(13) |
By introducing the equivalent permeability and applying the Darcy's law, the following equation is obtained:
(14) |
Applying Eqs. (12), (13), and (14) yields:
(15) |
where β is the high-velocity non-darcy coefficient. The correction of Evans and Civan [25] is used to determine the high-velocity non-Darcy coefficient,
(16) |
This correlation for β matched with more than 180 data points including those for propped fractures and was found to match the data very well with the correlation coefficient of 0.974 [26].
No. | Permeability, md | Porosity | Thickness, m |
---|---|---|---|
1 | 0.1 | 0.05 | 90 |
2 | 0.1 | 0.075 | 45 |
3 | 0.1 | 0.1 | 30 |
4 | 0.1 | 0.125 | 22.5 |
5 | 0.1 | 0.15 | 18 |
6 | 0.05 | 0.05 | 45 |
7 | 0.05 | 0.075 | 30 |
8 | 0.05 | 0.1 | 22.5 |
9 | 0.05 | 0.125 | 18 |
10 | 0.05 | 0.15 | 90 |
11 | 0.01 | 0.05 | 30 |
12 | 0.01 | 0.075 | 22.5 |
13 | 0.01 | 0.1 | 18 |
14 | 0.01 | 0.125 | 90 |
15 | 0.01 | 0.15 | 45 |
16 | 0.005 | 0.05 | 22.5 |
17 | 0.005 | 0.075 | 18 |
18 | 0.005 | 0.1 | 90 |
19 | 0.005 | 0.125 | 45 |
20 | 0.005 | 0.15 | 30 |
21 | 0.001 | 0.05 | 18 |
22 | 0.001 | 0.075 | 90 |
23 | 0.001 | 0.1 | 45 |
24 | 0.001 | 0.125 | 30 |
25 | 0.001 | 0.15 | 22.5 |
Fig. (2). The fitting results of incremental ratios for horizontal wells with different lengths. |
RESULTS AND DISCUSSION
Based on the above model, a simulator of simulating tight gas production was developed. After that, it was used to simulate the development of TGR using vertical well, horizontal well, or fractured-horizontal well.
Modeling of the Incremental Ratio (IR) of Horizontal Well to Vertical Well in Different TGRs
Based on the features of horizontal well and sensitivity analysis, the permeability, porosity, and reservoir thickness are the dominant factors of the IR of horizontal well to vertical well. Moreover, the length of the horizontal well is the key parameter of drilling cost and IRs. Thus, we first choose permeability, porosity, and reservoir thickness to set up simulation cases using Orthogonal Experimental Design (OED), as shown in Table 1. Then, the horizontal wells with different well lengths (L_{H}) for developing TGRs are performed to obtain the IRs.
The model of the IRs of horizontal well to vertical well is obtained based on the simulation results using data regression. The fitting results of IRs for horizontal wells with different lengths are shown in Fig. (2). We can see that as the permeability and reservoir thickness increase, the IRs decreases; as the porosity increases, the IRs increases. In order to extend the application of Eq. (17), we further obtain the relations between the coefficients a, b, c and L_{D} based on Fig. (3), which are shown as follows:
(17) |
(18) |
(19) |
(20) |
where R_{H/V} is the IR of horizontal well to vertical well; K is the intrinsic permeability; H is reservoir thickness; φ is the intrinsic porosity; L_{D} is the dimensionless length of horizontal well, which equals to the ratio of well length to reservoir length; a, b, and c are the coefficients. Then, the above models are used to predict the IRs of different TGRs, and the comparison of predicted value with simulated value is shown in Fig. (4). Thus, the above predicted model is validated.
Fig. (3). Coefficients a and b vs. Dimensionless length L_{D}. |
Fig. (4). Comparisons between predicted value and simulated value. |
Modeling of The IR of Fractured-horizontal Well to Horizontal Well in Different TGRs
Based on the features of fractured-horizontal well and sensitivity analysis, the permeability, porosity, reservoir thickness, and hydraulic fracture parameters are the dominant factors of the IRs of fractured-horizontal well to horizontal well. But the hydraulic fracture parameters also significantly affect the cost. Hence, we firstly choose permeability, porosity and reservoir thickness to set up simulation cases using Orthogonal Experimental Design (OED), as shown in Table 2. Then, we compare the IRs of fractured-horizontal well to horizontal well with different products of fracture number (n_{f}) and half-length (h_{f}) for developing TGR.
Fig. (5). The fitting results of IRs for horizontal wells with different hydraulic fractures n_{f}h_{f}.. |
The model of the IRs of fractured-horizontal well to horizontal well is obtained based on the simulation results. The fitting results of IRs with different products of n_{f} and h_{f} are shown in Fig. (5). We can see that as the permeability increases, the IRs decreases; as the porosity and reservoir thickness increase, the IRs increases. In order to extend the application of Eq. (21), we obtained the relation between the coefficients a, b, c and n_{f}h_{f} based on Fig. (6), which are shown as follows:
No. | Permeability, md | Porosity | Thickness, m |
---|---|---|---|
1 | 0.1 | 0.05 | 90 |
2 | 0.1 | 0.075 | 45 |
3 | 0.1 | 0.1 | 30 |
4 | 0.1 | 0.125 | 22.5 |
5 | 0.1 | 0.15 | 18 |
6 | 0.05 | 0.05 | 45 |
7 | 0.05 | 0.075 | 30 |
8 | 0.05 | 0.1 | 22.5 |
9 | 0.05 | 0.125 | 18 |
10 | 0.05 | 0.15 | 90 |
11 | 0.01 | 0.05 | 30 |
12 | 0.01 | 0.075 | 22.5 |
13 | 0.01 | 0.1 | 18 |
14 | 0.01 | 0.125 | 90 |
15 | 0.01 | 0.15 | 45 |
16 | 0.005 | 0.05 | 22.5 |
17 | 0.005 | 0.075 | 18 |
18 | 0.005 | 0.1 | 90 |
19 | 0.005 | 0.125 | 45 |
20 | 0.005 | 0.15 | 30 |
21 | 0.001 | 0.05 | 18 |
22 | 0.001 | 0.075 | 90 |
23 | 0.001 | 0.1 | 45 |
24 | 0.001 | 0.125 | 30 |
25 | 0.001 | 0.15 | 22.5 |
Fig. (6). Comparisons between calculated value and simulated value. |
(21) |
(22) |
(23) |
(24) |
where R_{Hf/Hw} is the IR of fractured-horizontal well to horizontal well; n_{f} is the fracture number; h_{f} is the half-length of the hydraulic fracture; a, b, and c are the coefficients. Then, the above models are used to predict the IRs of different TGRs with fractured-horizontal well, and the comparison of predicted value with simulated value is shown in Fig. (7). We can see that the agreement is pretty well.
Fig. (7). Coefficients a, b, and c vs. n_{f}h_{f}.. |
Model Application for Designing Projects
It is an essential problem to determine the well pattern or horizontal segment length for developing TGR. High profit is the ultimate objective of TGR development. From the economic perspective, horizontal well is not always better than vertical well, and a long horizontal well is not always better than a short horizontal well as well; Moreover, hydraulic fracturing is not always necessary due to its high cost.
Firstly, we study how to optimize the length of horizontal well. The cost of a horizontal well includes initial and production investment. In this paper, the initial investment which is the main difference between horizontal well and vertical well is taken into account. The initial investment is the sum of vertical segment cost, inclined segment cost, horizontal segment cost, and the well completion cost of horizontal segment [27]:
(25) |
Thus, the cost ratio (CR) of horizontal well to vertical well is as follows:
(26) |
where C_{tH} is the initial investment; C_{V}, C_{I}, C_{H} are the vertical segment cost, inclined segment cost, horizontal segment cost per meter, respectively; r_{h} is incremental coefficient of well completion cost of horizontal segment per square meter.
A higher cost must require a higher gas production. We introduce the coefficient ξ as the ratio of minimum incremental ratio of production (PIR) to the incremental ratio of cost (CIR). In fact, ξ is related to the gas price, and a higher gas price permits a smaller ξ. ξ can be described as:
(27) |
Hence, the following relation must be met for achieving an ideal profit with horizontal well,
(28) |
We first use the above models to study the feasibility of horizontal well in different reservoirs with different scenarios. The related parameters are shown in Table 3. Then, we can obtain the plate of different reservoir thicknesses with different ξ, as shown in Fig. (8). We can see that the plate can be used to determine the reasonable well length for different reservoirs. Thin reservoir is more suitable to use a horizontal well, and longer well does not mean better profit.
Parameters | Value |
---|---|
Reservoir length, m | 1000 |
Porosity | 0.1 |
Permeability, md | 0.01 |
Reservoir thickness, m | 20, 40, 60, 80 |
C_{V}, $/m | 500 |
C_{I}, $/m | 2000 |
C_{H}, $/m | 2250 |
L_{V}, m | 3000 |
L_{I}, m | 300 |
L_{H}, m | 200-1000 |
r_{h}, m^{-2} | 0.01 |
Fig. (8). Dimensionless well length vs. Production ratios of different thicknesses and ξ. |
As it is discussed above, hydraulic fractures also need be designed. The cost of hydraulic fracturing is higher than 5,000,000 $ in USA, which is twice of well drilling. Similarly, we can compare the initial investment of fractured-horizontal well with horizontal well. The initial investment is the sum of vertical segment cost, inclined segment cost, horizontal segment cost, and hydraulic fracturing:
(29) |
Thus, the cost ratio of fractured-horizontal well to horizontal well is as follows,
(30) |
where C_{F} is the cost of hydraulic fracturing per meter. Similarly, the following relation must be met for achieving an ideal profit with fractured-horizontal well,
(31) |
We first use the above models to study the feasibility of fractured-horizontal well in different reservoirs with different scenarios. The related parameters are shown in Table 4. Then, we can obtain the plate of different permeability with different ξ, as shown in Fig. (9). We can find that the plate can be used to determine the hydraulic fracture parameters; and the lower the permeability is, the more necessary to carry out hydraulic fracturing.
Parameters | Value |
---|---|
Reservoir length, m | 1000 |
Well length, m | 800 |
Porosity | 0.1 |
Reservoir thickness, m | 70 |
C_{V}, $/m | 500 |
C_{I}, $/m | 2000 |
C_{H}, $/m | 2250 |
C_{F}, $/m | 5000 |
L_{V}, m | 3000 |
L_{I}, m | 300 |
L_{H}, m | 200-1000 |
Fig. (9). Product of fracture number and half-length vs. Production ratios of different permeability and ξ. |
Cases Study
Two TGRs in China are used to perform the design of different well patterns. The parameter is shown in Table 5.
Fig. (10). Production ratios vs. Dimenssionless well length in TGR A. |
Fig. (11). Production ratios vs. Dimenssionless well length in TGR B. |
Paramters | TGR A | TGR B |
---|---|---|
Permeability, md | 0.01 | 0.002 |
Reservoir length, m | 1500 | 1500 |
Reservoir thickness, m | 30 | 35 |
Reservoir depth, m | 3000 | 3200 |
Porosity | 0.08 | 0.1 |
Well length, m | 1200 |
In TGR A, the permeability is 0.01md, which is relatively high for gas flow into well, so hydraulic fracturing may be not necessary, and vertical well and horizontal well are compared. Fig. (10) shows the production ratios of horizontal well to vertical well with different well lengths. We can see that when L_{D} is less than 0.2, the production ratio is less than 1, vertical well is apparently better. For ξ=0.3, L_{D} must be larger than 0.4 to satisfy the requirement that R_{H/V}> R_{H/V}^{min}; Similarly, L_{D} must be larger than 0.6 for ξ=0.4; whereas horizontal well is always unreasonable for ξ=0.5.
In TGR B, the permeability is 0.002md. We first evaluate the feasibility of horizontal well. Fig. (11) shows the production ratios of horizontal well to vertical well with different well lengths. We can see that R_{H/V} is much greater than R_{H/V}^{min} even for ξ=1.0, so horizontal well is significantly necessary. Then, we evaluate the feasibility of fractured-horizontal well. Fig. (12) shows the production ratios of fractured-horizontal well to horizontal well with different products of fracture number and half-length. We can see that when ξ=0.3, the product should be less than 900 to satisfy the requirement that R_{Hf/Hw}> R_{Hf/Hw}^{min}; Similarly, it must be less than 600 and 400 and for ξ=0.4 and ξ=0.5, respectively.
Fig. (12). Production ratios vs. Product of fracture number and half-length in TGR B. |
CONCLUSION
In this paper, we first developed a simulator for simulating the development of TGRs including slippage flow and stress dependence in matrix, and high-velocity non-Darcy flow in hydraulic fractures. Then, we studied the development effects of different well patterns and well parameters for developing TGRs. The incremental ratio models of horizontal well to vertical well and fractured-horizontal well to horizontal well were achieved. Moreover, some conclusions were achieved:
(1) The models can be used to obtain the plates of choosing a suitable well type and corresponding parameters;
(2) The effect of horizontal well is significant in the tight gas reservoirs with low permeability, high porosity and small thickness;
(3) The effect of fractured-horizontal well is significant in the tight gas reservoirs with low permeability, high porosity and large thickness;
(4) Longer horizontal well or larger n_{f}h_{f} does not mean to achieve a better economic benefit.
CONFLICT OF INTEREST
The authors confirm that this article content has no conflict of interest.
ACKNOWLEDGEMENTS
This study was supported by Science Foundation of CNPC Research Institute of Petroleum Exploration & Development, Research on Key Technologies of Natural Gas Development (No. 2014B-1507).