Introduction

Hydraulic fracturing is a very common way to stimulate hydrocarbon wells. Also, well test analysis is the most used practice to evaluate the results of a fracturing job. Well test analysis is mostly based upon analytical solutions which allow us to understand the formation and well flow dynamics. For well pressure tests, analytical models can be grouped into two types according to their solution: semi-analytic models and asymptotic analytical models. A semi-analytic model was first developed by Gringarten, Ramey and Raghavan (1974) for hydraulically fractured wells with fractures of infinite conductivity, and later, this semi-analytic model was adapted for hydraulically fractured wells with finite-conductivity fractures by Cinco-Ley, Samaniego and Dominguez (1978). They developed a mathematical model to study transient behavior in wells with vertical finite-conductivity fractures for dimensional times greater than 0.001, which may be correlated with the dimensional finite conductivity of fractures. Their provided a type-curve matching interpretation technique.

Subsequently, Cinco-Ley and Samaniego (1981) developed an asymptotic analytical solution for the interpretation of transient pressure data for wells intercepted by a vertical fracture. With this, they demonstrated that for dimensional fracture conductivity less than 300, the fracture conductivity cannot be considered as infinite. They also provided type curves for the determination of the reservoir and fracture parameters.

Later, Cinco-Ley (1982) developed a method for fractured wells considering different flow geometries. This work recognized that wells intercepted by a fracture had different flow behaviors than non-fractured wells, and application of traditional pressure analysis methods based upon radial flow theory, could yield erroneous results. His method was based on the bilinear flow rate, which is present when the linear flow of the fracture and the linear flow within the formation occur simultaneously.

Several other models were also introduced. For instance, Cossio, Moridis, and Blasingame (2013) presented a semi-analytic solution for a finite-conductivity fractured using fractal geometry theory. Wong, Harrington and Cinco-Ley (1984) presented a mathematical solution for pressure behavior under pseudolinear conditions which has useful applications when the fracture has a skin factor that allows having a linear flow from the fracture tip to the well.

As a final flow model mention, there are more, Lee and Brockenbrough (1986) introduced a new analytical solution, based on a trilinear flow model as shown in Fig. 1. This sketch shows how three linear flows make up the trilinear flow model. Flow in region 1 renders the fracture with a dominant flow in the x direction. If the fracture is considered as a bending line, the flow in region 2 will be controlled by the y direction. With similar reasoning, the flow in region 3 would be approximate in the x direction for a very short period. These reasoning of flow behavior is reasonable for fluid flow before the pseudo-radial flow regime shows up. In addition to the formulation of their trilinear flow model, they also conducted a study comparing this trilinear flow solution with that proposed by Cinco-Ley et al. (1981). They found that the trilinear model fits perfectly with the model proposed by Cinco-Ley et al. (1981) up to a value of tD = 1.0. After this value, the deviation becomes more appreciable.

Fig. 1
figure 1

Taken from Lee and Brockenbrough (1986)

A quadrant of the top view of the fractured well system and trilinear flow approximation.

The authors are focused on well test interpretation by the TDS Technique, Tiab (1995), which is a novel, accurate and revolutionary methodology based upon characteristic points and lines found on the pressure and pressure derivative versus time log–log plot -also to the reciprocal rate and reciprocal rate derivative versus time log–log plot. Details of this technique is widely exposed on the textbooks of Escobar (2015, 2019) and a very comprehensive state-of-the-art is presented by Escobar, Jongkittnarukorn and Hernandez (2018). Classical articles on the TDS Technique for hydraulic fracture wells are presented by Tiab (1994) and Tiab et al (1999). Also, Escobar, Gonzalez and Hernandez (2016) extended the TDS Technique for well test interpretation under pseudolinear conditions using the already mentioned model by Wong et al (1984).

In this work, TDS Technique and straight-line conventional analysis are extended for the practical test data interpretation of finite-conductivity vertical fractured wells in hydrocarbons reservoirs using the trilinear flow model of Lee and Brockenbrough (1986). Both pressure transient and rate transient analysis are provided for the TDS Technique and conventional analysis, as well. The proposed equations were successfully tested with synthetic and actual field data.

Mathematical formulation

Lee and Brockenbrough (1986) provided a Laplace space pressure solution of a vertical fractured well where the flow is dominated by trilinear flow regime,

$$\overline{P}_{wD} (s) = \frac{b}{{s\left[ {sbC_{D} - \psi \;\tanh \,\psi } \right]}}$$
(1)

where;

$$a \equiv \frac{2}{{C_{fD} }};\quad b \equiv \frac{\pi }{{C_{fD} }}$$
(2)
$$\psi = \sqrt {\frac{a(s + \sqrt s )}{{1 + \sqrt {(s + \sqrt s )} S}} + \frac{s}{{\eta_{f} }}}$$
(3)
$$\frac{1}{{\eta_{f} }} = \frac{{k\phi_{f} c_{tf} }}{{k_{f} \phi c_{t} }}$$
(4)

The dimensionless quantities are defined below as:

$$C_{D} = \frac{0.894C}{{\phi c_{t} hx_{f}^{2} }}$$
(5)
$$C_{fD} = \frac{{k_{f} w_{f} }}{{kx_{f} }}$$
(6)

The dimensionless parameters for an oil reservoir are given below:

$$t_{{Dx_{f} }} = \frac{0.0002637kt}{{\phi \mu c_{t} x_{f}^{2} }}$$
(7)
$$P_{D} = \frac{{kh(P_{i} - P_{wf} )}}{141.2q\mu B}$$
(8)
$$t_{D} *P_{D} ^{\prime} = \frac{kh(t*\Delta P^{\prime})}{{141.2q\mu B}}$$
(9)

The solution of the flow equation for the constant-pressure case;

$$\overline{q}_{D} (s) = \frac{1}{{s^{2} \overline{P}_{wD} (s)}}$$
(10)

Then, the dimensionless parameters are given by:

$$\frac{1}{{q_{D} }} = \frac{{kh\left( {P_{i} - P_{wf} } \right)}}{141.2q\mu B}$$
(11)
$$t_{D} *(1/q_{D} )^{\prime} = \frac{kh\Delta P}{{141.2\mu B}}[t*(1/q)^{\prime}]$$
(12)

The dimensionless quantities for gas reservoirs, constant-rate case, are given as follows:

$$t_{{Dx_{f} }} = \frac{0.0002637kt}{{\phi (\mu c_{t} )_{i} x_{f}^{2} }}$$
(13)
$$m(P)_{D} = \frac{{kh\left[ {m(P_{i} ) - m(P_{wf} )} \right]}}{{1424.26q_{sc} T}}$$
(14)
$$t_{D} *m(P)_{D} ^{\prime} = \frac{{kh\left[ {t*\Delta m(P)^{\prime}} \right]}}{{1424.26Tq_{sc} }}$$
(15)

The dimensionless parameters for the constant-pressure case, which are used in transient rate analysis, are provided below:

$$\frac{1}{{q_{D} }} = \frac{{kh\left[ {m\left( {P_{i} } \right) - m(P_{wf} )} \right]}}{1424.26T}\frac{1}{{q_{g} }}$$
(16)
$$\left[ {t_{D} *(1/q_{D} )^{\prime}} \right] = \frac{{kh\left[ {\Delta m\left( P \right)} \right]}}{1424.26T}\left[ {t*(1/q)^{\prime}} \right]$$
(17)

Wellbore rate and pressure behavior

Both dimensionless pressure and dimensionless pressure derivative versus dimensionless time behavior provided in Fig. 2 was generated with Eq. 1 while the reciprocal rate and reciprocal rate derivative behavior was obtained by using Eq. 10. In both cases the transient behaviors are about the same. As indicated in both plots, at early time the linear flow regime dominates the test. This is followed by a unit-slope tendency observed on both derivatives, although is shorter for the reciprocal rate case and finally a bilinear flow regime is developed. During the intermediate time when reciprocal rate derivative displays a unit-slope behavior, see Fig. 3, the reciprocal rate does not display a unit-slope behavior but a variety of slopes going to values higher than 2. This, therefore, difficult the application of conventional analysis.

Fig. 2
figure 2

Log–log plot of dimensionless pressure and dimensionless pressure derivative versus dimensionless time for different dimensionless conductivities

Fig. 3
figure 3

Log–log plot of dimensionless reciprocal rate and dimensionless reciprocal rate derivative versus dimensionless time for different dimensionless conductivities

For a better understanding of the flow dynamics let us refer to Fig. 1 which sketches the trilinear model proposed by Lee and Brockenbrough (1986). They consider three flow zone. The first zone provides the pressure response of the fracture. The second zone is responsible for the simultaneous answer of two linear flows: one occurring inside the fracture and the other one from the formation to the fracture. Finally, region 3 consists of linear flow to the well-fracture system. As reflected in both pressure derivative and reciprocal rate derivative, the first flow regime occurs corresponds to a linear flow from formation to fracture. This means that the flow inside the fracture is imperceptible. It looks that the unit-slope behavior results from the transition between the early linear flow regime to the subsequent bilinear flow regime which takes place at later time. This bilinear flow may be formed by the simultaneous flow regime from the formation to the fracture and from zone 3 to the fracture-well system.

Pressure transient analysis by the TDS technique for oil reservoirs

TDS Technique, Tiab (1995), is a methodology based upon finding unique characteristic and points on the pressure and pressure derivative versus time log–log plot (this is extended to reciprocal rate and reciprocal rate derivative). Escobar (2015, 2019) presented two textbooks which clearly illustrate this technique.

For extension of the TDS Technique is necessary to find unified pressure/rate and pressure/reciprocal rate derivative behaviors so unique equations can be developed. For instance, Fig. 4 presents a unified behavior of the pressure derivative multiplied by the dimensionless fracture conductivity. The following equation is.

Fig. 4
figure 4

Unified pressure derivative times dimensionless conductivity versus dimensionless time during early linear flow regime

then obtained:

$$(t_{D} *P_{D} ^{\prime})_{L} C_{fD} = 0.6\pi (t_{D} )_{L}^{0.5}$$
(18)

Suffix L in Eq. (18) stands for “Linear”. Plugging Eqs. (6), (7) and (9) into Eq. (18) and solving for the fracture conductivity, it yields:

$$k_{f} w_{f} = \frac{4.3q\mu B}{{h(t*\Delta P^{\prime})_{L} }}\sqrt {\frac{{kt_{L} }}{{\phi \mu c_{t} }}}$$
(19)

Integration of Eq. (18) leads to find the governing pressure equation:

$$(P_{D} )_{L} = \frac{{0.6\pi (t_{D} )_{L}^{0.5} }}{{0.5C_{fD} }} + c$$
(20)

where c is a given constant of no interest for this work. Again, replacing Eqs. (6, 7 and 8) into Eq. (20), another equation is found to find the fracture conductivity.

$$k_{f} w_{f} = \frac{8.6q\mu B}{{h(\Delta P)_{L} }}\sqrt {\frac{{kt_{L} }}{{\phi \mu c_{t} }}}$$
(21)

Another unified behavior is found during the unit-slope behavior as shown in Fig. 5. The following expression governs such behavior:

$$(t_{D} *P_{D} ^{\prime})_{{{\text{US}}}} = 0.9\pi \frac{{(t_{D} )_{{{\text{US}}}} }}{{C_{fD} }}$$
(22)
Fig. 5
figure 5

Unified pressure derivative times dimensionless conductivity versus dimensionless time during the unit-slope flow behavior

Again, we found by integration that the unified dimensionless pressure equation is governed by:

$$(P_{D} )_{{{\text{US}}}} = 0.9\pi \frac{{(t_{D} )_{{{\text{US}}}} }}{{C_{fD} }}$$
(23)

Replacement of the dimensionless quantities given by Eqs. 7, 8, 9 into Eqs. (22) and (23) allows to find two expressions for the determination of the half-fracture length:

$$x_{f} = \frac{0.1qBk}{{h\phi c_{t} k_{f} w_{f} }}\left( {\frac{{t_{{{\text{US}}}} }}{{(t*\Delta P^{\prime})_{{{\text{US}}}} }}} \right)$$
(24)
$$x_{f} = \frac{0.1qBk}{{h\phi c_{t} k_{f} w_{f} }}\left( {\frac{t}{\Delta P}} \right)_{{{\text{US}}}}$$
(25)

As seen in Fig. 2 a unified behavior exists by itself during linear flow regime. The pressure and pressure derivative governing equations during the bilinear flow regime are found to be:

$$(t_{D} *P_{D} ^{\prime})_{L} = 0.43414(t_{{{\text{Dxf}}}}^{0.25} )_{L}$$
(26)
$$(P_{D} )_{BL} = 1.736554(t_{Dxf}^{0.25} )_{BL}$$
(27)

Replacing the dimensionless parameters in the above equations, it yields:

$$x_{f} = 60.254\left( {\frac{q\mu B}{{{\text{kh}}(t*\Delta P^{\prime})_{{{\text{BL}}}} }}} \right)^{2} \sqrt {\frac{{{\text{kt}}_{{{\text{BL}}}} }}{{\phi \mu c_{t} }}}$$
(28)
$$x_{f} = 976.35\left( {\frac{q\mu B}{{{\text{kh}}(\Delta P)_{{{\text{BL}}}} }}} \right)^{2} \sqrt {\frac{{kt_{{{\text{BL}}}} }}{{\phi \mu c_{t} }}}$$
(29)

We may also solve for formation permeability,

$$k = 15.3694\left( {\frac{q\mu B}{{h(t*\Delta P^{\prime})_{{{\text{BL}}}} }}} \right)^{4/3} \left( {\frac{{t_{{{\text{BL}}}} }}{{\phi \mu c_{t} x_{f}^{2} }}} \right)^{1/3}$$
(30)
$$k = 98.417\left( {\frac{q\mu B}{{h(\Delta P)_{{{\text{BL}}}} }}} \right)^{4/3} \left( {\frac{{t_{{{\text{BL}}}} }}{{\phi \mu c_{t} x_{f}^{2} }}} \right)^{1/3}$$
(31)

Pressure transient analysis by the conventional method for oil reservoirs

Once the dimensionless parameters are replaced in Eq. 20, the following is obtained:

$$\Delta P = P_{i} - P_{{{\text{wf}}}} = \frac{8.6q\mu B}{{hk_{f} w_{f} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}} \sqrt t$$
(32)

Which can be expressed as:

$$P_{{{\text{wf}}}} = P_{i} - m_{L} \sqrt t$$
(33)

The above implies that a Cartesian plot of \(\Delta P\, = \,\left( {P_{i} - P_{{{\text{ws}}}} } \right){\text{ vs}}.\;t^{0.5}\) (for drawdown) or \(\Delta P\, = \,\left( {P_{{{\text{ws}}}} - P_{{{\text{wf}}}} } \right){\text{ vs}}. \, \left[ {\left( {t_{p} \, + \,\Delta t} \right)^{0.5} {-\!\!-}\Delta t^{0.5} } \right]\) (for buildup) yields a straight line which slope, mL, provides the fracture conductivity. Such slope, \(m_{L}\) is given by:

$$m_{L} = \frac{8.6q\mu B}{{hk_{f} w_{f} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}}$$
(34)

From which;

$$k_{f} w_{f} = \frac{8.6q\mu B}{{{\text{hm}}_{L} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}}$$
(35)

The unit-slope pressure equation obtained from Eq. 23 allows finding:

$$P_{{{\text{wf}}}} = P_{i} - \frac{{0.1q{\text{Bk}}}}{{h\phi c_{t} k_{f} w_{f} }}t$$
(36)
$$x_{f} = \frac{{0.1q{\text{Bk}}}}{{h\phi c_{t} k_{f} w_{f} }}\left( {\frac{t}{\Delta P}} \right)_{{{\text{US}}}}$$
(37)

Which can be expressed as:

$$P_{{{\text{wf}}}} = P_{i} - m_{{{\text{US}}}} t$$
(38)

A Cartesian plot of ΔP vs. t yields a straight line which slope, mUS, provides the half-fracture length. This slope is given by:

$$m_{{{\text{US}}}} = \frac{{0.1q{\text{Bk}}}}{{h\phi c_{t} k_{f} w_{f} x_{f} }}$$
(39)

From which;

$$x_{f} = \frac{0.1qBk}{{h\phi c_{t} k_{f} w_{f} m_{US} }}$$
(40)

The bilinear pressure behavior obtained from Eq. 27 allows to develop:

$$P_{{{\text{wf}}}} = P_{i} - 31.247\frac{q\mu B}{{{\text{kh}}}}\sqrt[4]{{\frac{k}{{\phi \mu c_{t} x_{f}^{2} }}}}\sqrt[4]{t}$$
(41)

Which can be expressed as:

$$P_{{{\text{wf}}}} = P_{i} - m_{{{\text{BL}}}} \sqrt[4]{t}$$
(42)

Which implies that a Cartesian plot of \(\Delta P\, = \,\left( {P_{{{\text{ws}}}} - P_{{{\text{wf}}}} } \right){\text{ vs}}. \, \left[ {\left( {t_{p} \, + \,\Delta t} \right)^{0.25} - \Delta t^{0.25} } \right]\) (for drawdown) or \(\Delta P\, = \,\left( {P_{{{\text{ws}}}} - P_{{{\text{wf}}}} } \right){\text{ vs}}. \, \left[ {\left( {t_{p} \, + \,\Delta t} \right)^{0.25} - \Delta t^{0.25} } \right]\) (for buildup) provides a straight line which slope, mBL, provides the half-fracture length:

$$x_{f} = 976.35\left( {\frac{q\mu B}{{{\text{kh}}}}} \right)^{2} \sqrt[{}]{{\frac{k}{{\phi \mu c_{t} }}}}$$
(43)

Pressure transient analysis by the TDS technique for gas reservoirs

Equations 18 and 20, 22 and 23 and 26 and 27 apply to gas reservoirs. Dimensionless pressure and dimensionless pressure derivative are changed by dimensionless pseudopressure and dimensionless pseudopressure derivative. Then, after replacing the dimensionless time expression given by Eq. 13 and the dimensionless pseudopressure and pseudopressure derivative given by Eqs. 14 and 15 in the above-mentioned equations, it yields:

$$k_{f} w_{f} = \frac{{43.594q_{{{\text{sc}}}} T}}{{{\text{kh}}\left[ {t*\Delta m(P)^{\prime}} \right]}}\sqrt {\frac{{{\text{kt}}_{L} }}{{\phi (\mu c_{t} )_{i} }}}$$
(44)
$$k_{f} w_{f} = \frac{{87.2q_{{{\text{sc}}}} T}}{{h[\Delta m(P)]_{L} }}\sqrt {\frac{{{\text{kt}}_{L} }}{{\phi \mu c_{t} }}}$$
(45)
$$x_{f} = \frac{{1.062q_{{{\text{sc}}}} kT}}{{\phi \mu c_{t} {\text{hk}}_{f} w_{f} }}\left[ {\frac{t}{t*\Delta m(P)^{\prime}}} \right]_{US}$$
(46)
$$x_{f} = \frac{{1.062q_{{{\text{sc}}}} kT}}{{\phi \mu c_{t} {\text{hk}}_{f} w_{f} }}\left[ {\frac{t}{\Delta m(P)}} \right]_{{{\text{US}}}}$$
(47)
$$x_{f} = 6208.565\left( {\frac{qT}{{kh[t*\Delta m(P^{\prime})]_{{{\text{BL}}}} }}} \right)^{2} \sqrt {\frac{{kt_{{{\text{BL}}}} }}{{\phi (\mu c_{t} )_{i} }}}$$
(48)
$$k = 337.8\left( {\frac{qT}{{h[t*\Delta m(P^{\prime})]_{{{\text{BL}}}} }}} \right)^{4/3} \left( {\frac{{t_{{{\text{BL}}}} }}{{\phi (\mu c_{t} )_{i} x_{f}^{2} }}} \right)^{1/3}$$
(49)
$$x_{f} = 6208.565\left( {\frac{qT}{{{\text{kh}}[t*\Delta m(P^{\prime})]_{BL} }}} \right)^{2} \sqrt {\frac{{kt_{{{\text{BL}}}} }}{{\phi (\mu c_{t} )_{i} }}}$$
(50)
$$k = 2144.915\left( {\frac{qT}{{h[\Delta m(P)]_{{{\text{BL}}}} }}} \right)^{4/3} \left( {\frac{{t_{{{\text{BL}}}} }}{{\phi (\mu c_{t} )_{i} x_{f}^{2} }}} \right)^{1/3}$$
(51)

The points of intersection between the different pressure derivative equations are very important for the verification of the reservoir parameters. These developed expressions are valid for both oil and gas cases, with the exception that for gas case the product of compressibility and viscosity is given at initial conditions. Thereby, the intersection points between Eqs. 18 and 22, 18 and 27, 22 and 27 and 11 and 27 provides the following expressions:

$$x_{f} = \frac{1}{41.054}\sqrt {\frac{{kt_{LUSi} }}{{\phi \mu c_{t} }}}$$
(52)
$$k = 1.606\sqrt[5]{{\frac{{\phi \mu c_{t} (k_{f} w_{f} )^{4} }}{{x_{f}^{2} t_{LBLi} }}}}$$
(53)
$$k = 11.7166\sqrt[7]{{\left( {\frac{{\phi \mu c_{t} }}{{t_{USBLi} }}} \right)^{3} x_{f}^{2} (k_{f} w_{f} )^{4} }}$$
(54)
$$k = \frac{{670.606\phi \mu c_{t} k_{f} w_{f} x_{f} }}{{t_{USRi} }}$$
(55)

If radial flow is seen, Tiab (1995) found that during the radial flow regime the pressure derivative is governed by the following simple equation:

$$(t_{D} *P_{D} ^{\prime})_{r} = 0.5$$
(56)

Then, the intersection between Eqs. 18, 22 and 27 with Eq. 53 yields the following:

$$k = 6.438\sqrt[3]{{\frac{{\phi \mu c_{t} (k_{f} w_{f} )^{2} }}{{t_{LRi} }}}}$$
(57)
$$k = \frac{{670.606\phi \mu c_{t} k_{f} w_{f} x_{f} }}{{t_{USRi} }}$$
(58)
$$x_{f} = \frac{1}{81.6822}\sqrt {\frac{{kt_{BLRi} }}{{\phi \mu c_{t} }}}$$
(59)

Transient rate analysis by the TDS technique for oil reservoirs

As shown in Fig. 6, the unified reciprocal rate derivative behavior obeys the following equation:

$$[t_{D} *(1/q_{D} )^{\prime}]_{L} = 1.77\pi (t_{D} /C_{fD}^{2} )_{L}^{0.5}$$
(60)
Fig. 6
figure 6

Unified reciprocal rate derivative versus dimensionless time divided by the dimensionless fracture conductivity squared observed during early linear flow regime

Replacing Eqs. 7 and 12 allows to solve for an expression to find fracture conductivity,

$$k_{f} w_{f} = \frac{4.034q\mu B}{{h[t*(1/q)^{\prime}]_{L} }}\sqrt {\frac{{{\text{kt}}_{L} }}{{\phi \mu c_{t} }}}$$
(61)

Integration of Eq. 56 leads to find the governing expression for the reciprocal rate,

$$(1/q_{D} )_{L} = 1.12\pi (t_{D} /C_{fD}^{2} )_{L}^{0.5}$$
(62)

After replacing Eqs. 7 and 11 and solving for the fracture conductivity from the resulting replacement leads to,

$$k_{f} w_{f} = \frac{12.75\mu B}{{h\Delta P(1/q)_{L} }}\sqrt {\frac{{{\text{kt}}_{L} }}{{\phi \mu c_{t} }}}$$
(63)

Also, another unified reciprocal rate derivative behavior is given in Figs. 7 and 8 from which can be written:

$$[t_{D} *(1/q_{D} )^{\prime}]_{US} = 0.9\pi (t_{D} /C_{fD} )$$
(64)
Fig. 7
figure 7

Unified reciprocal rate derivative versus dimensionless time divided by dimensionless fracture conductivity observed unit-slope behavior

Fig. 8
figure 8

Unified reciprocal rate versus dimensionless time divided by the square root of dimensionless fractured conductivity both raised squared, observed during the unit-slope behavior of the reciprocal rate derivative

Which also provides an expression for the estimation of the half-fracture length once Eqs. 7 and 12 are replaced in the above equation. Therefore,

$$x_{f} = \frac{{0.1053q{\text{Bk}}}}{{\phi c_{t} h\,k_{f} w_{f} }}\left( {\frac{t}{[t*(1/q)^{\prime}]}} \right)_{{{\text{US}}}}$$
(65)

It is observed during the unit-slope behavior that the reciprocal rate does not follow a unit slope, then the Equation for this was found to be,

$$(1/q_{D} )]_{{{\text{US}}}} = \frac{{1.8772(t_{D}^{2} )_{{{\text{US}}}} }}{{\sqrt {C_{fD} } }}$$
(66)

It should be used only if the log–log slope of the reciprocal rate is about 2. From where an expression to find the half-fracture length is found once the dimensionless quantities have been replaced in the above equation,

$$x_{f} = \frac{1}{37.86}\sqrt[3]{{\left( {\frac{\mu B}{{kh(\Delta P)k_{f} w_{f} (1/q)_{{{\text{US}}}} }}} \right)\left( {\frac{{kt_{{{\text{US}}}} }}{{\phi \mu c_{t} }}} \right)^{2} }}$$
(67)

A unified behavior is also found during the bilinear flow regime in Fig. 3. The governing expression for this flow regime is given by:

$$[t_{D} *(1/q_{D} )^{\prime}]_{{{\text{BL}}}} = 0.43414(t_{Dxf}^{0.25} )_{{{\text{BL}}}}$$
(68)

An expression to find the half-fractured length is obtained once the dimensionless time, Eq. 7, and dimensionless reciprocal rate derivative, Eq. 12, are plugged in Eq. 64, therefore:

$$x_{f} = 61.02\left( {\frac{\mu B}{{{\text{kh}}\Delta P[t*(1/q)^{\prime}]_{{{\text{BL}}}} }}} \right)^{2} \sqrt {\frac{{kt_{{{\text{BL}}}} }}{{\phi \mu c_{t} }}}$$
(69)

Or,

$$k = 15.5\left( {\frac{\mu B}{{{\text{kh}}\Delta P[t*(1/q)^{\prime}]_{{{\text{BL}}}} }}} \right)^{4/3} \left( {\frac{{t_{{{\text{BL}}}} }}{{\phi \mu c_{t} }}} \right)^{1/3}$$
(70)

Integrating Eq. 64 results in,

$$(1/q_{D} )_{BL} = 1.736554(t_{Dxf}^{0.25} )_{BL} + c$$
(71)

Again, the constant, c, has no effect in this work. Replacing Eqs. 7 and 11 in the above expression leads to find either:

$$x_{f} = 976.35\left( {\frac{\mu B}{{kh\Delta P(1/q)_{BL} }}} \right)^{2} \sqrt {\frac{{kt_{BL} }}{{\phi \mu c_{t} }}}$$
(72)

Or,

$$k = 98.417\left( {\frac{\mu B}{{kh\Delta P(1/q_{{{\text{BL}}}} )}}} \right)^{4/3} \left( {\frac{{t_{{{\text{BL}}}} }}{{\phi \mu c_{t} }}} \right)^{1/3}$$
(73)

Rate transient analysis by the conventional method for oil reservoirs

Equation 58 can be rewritten as:

$$1/q = \frac{12.75\mu B}{{h\Delta Pk_{f} w_{f} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}} \sqrt t$$
(74)

Which can be expressed as:

$$\frac{1}{q} = m_{L} \sqrt t$$
(75)

A Cartesian plot of 1/q vs. t0.5 yields a straight line which slope, mL, provides the fracture conductivity,

$$m_{L} = \frac{12.75\mu B}{{h\Delta P\;k_{f} w_{f} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}}$$
(76)

Therefore,

$$k_{f} w_{f} = \frac{12.75\mu B}{{h\Delta Pm_{L} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}}$$
(77)

Once the dimensionless parameters have been replaced into Eq. 62, this be written as:

$$(1/q) = \frac{1}{54255.27}\left( {\frac{\mu B}{{\phi c_{t} h(\Delta P)k_{f} w_{f} x_{f}^{3} }}} \right)\left( {\frac{k}{{\phi \mu c_{t} }}} \right)^{2} t^{2}$$
(78)

Or,

$$\frac{1}{q} = m_{{{\text{US}}}} \;t^{2}$$
(79)

A Cartesian plot of 1/q vs. t2 yields a straight line which slope, mUS, provides the half fracture length,

$$x_{f} = \frac{1}{37.86}\sqrt[3]{{\left( {\frac{\mu B}{{\phi c_{t} h(\Delta P)k_{f} w_{f} m_{{{\text{US}}}} }}} \right)\left( {\frac{k}{{\phi \mu c_{t} }}} \right)^{2} }}$$
(80)

Equation 67 is rewritten as:

$$\frac{1}{q} = 31.247\frac{\mu B}{{{\text{kh}}\Delta P}}\sqrt[4]{{\frac{k}{{\phi \mu c_{t} x_{f}^{2} }}}}\sqrt[4]{t}$$
(81)

Or,

$$\frac{1}{q} = m_{{{\text{BL}}}} \sqrt[4]{t}$$
(82)

A Cartesian plot of 1/q vs. t0.25 yields a straight line which slope, mBL, provides the half fracture length,

$$x_{f} = 31.247\left( {\frac{\mu B}{{kh\Delta Pm_{{{\text{BL}}}} }}} \right)^{2} \sqrt[{}]{{\frac{k}{{\phi \mu c_{t} }}}}$$
(83)

Transient rate analysis by the TDS technique for gas reservoirs

Equations 56, 58, 60, 62,64 y 67also work for gas reservoirs. Now, the reciprocal rate and reciprocal rate derivative are affected by pseudopressure instead of pressure and Eqs. 16 and 17 are employed instead of Eqs. 11 and 12. Under these circumstances the resulting equations are:

$$k_{f} w_{f} = \frac{40.69T}{{h[\Delta m(P)][t*(1/q)^{\prime}]_{L} }}\sqrt {\frac{{{\text{kt}}_{L} }}{{\phi \mu c_{t} }}}$$
(84)
$$k_{f} w_{f} = \frac{81.36T}{{{\text{kh}}[\Delta m(P)](1/q)_{L} }}\sqrt {\frac{{{\text{kt}}_{L} }}{{\phi \mu c_{t} }}}$$
(85)
$$x_{f} = \frac{{1.062kTt_{{{\text{US}}}} }}{{\phi \mu c_{t} h[\Delta m(P)]k_{f} w_{f} [t*(1/q)^{\prime}]_{{{\text{US}}}} }}$$
(86)
$$k_{f} w_{f} = \frac{T}{{34460.6kh[\Delta m(P)])(1/q)_{{{\text{US}}}} }}\left( {\frac{{kt_{{{\text{US}}}} }}{{\phi \mu c_{t} }}} \right)^{2}$$
(87)
$$x_{f} = 6208.565\left( {\frac{T}{{kh\Delta m(P)[t*(1/q)^{\prime}]_{BL} }}} \right)^{2} \sqrt {\frac{{kt_{BL} }}{{\phi (\mu c_{t} )_{i} }}}$$
(88)

Or,

$$k = 337.8\left( {\frac{T}{{{\text{kh}}\Delta m(P)[t*(1/q)^{\prime}]_{{{\text{BL}}}} }}} \right)^{4/3} \left( {\frac{{{\text{kt}}_{{{\text{BL}}}} }}{{\phi (\mu c_{t} )_{i} x_{f}^{2} }}} \right)^{1/3}$$
(89)

Or,

$$k = 2144.915\left( {\frac{T}{{{\text{kh}}\Delta m(P)(1/q)_{{{\text{BL}}}} }}} \right)^{4/3} \left( {\frac{{kt_{{{\text{BL}}}} }}{{\phi (\mu c_{t} )_{i} x_{f}^{2} }}} \right)^{1/3}$$
(90)

Intersection points also play an important role in rate transient analysis when the TDS Technique is applied. The following expressions resulted from the reciprocal rate derivative intersection points between the different flow regimes:

$$x_{f} = \frac{1}{48.576}\sqrt {\frac{{{\text{kt}}_{LUSi} }}{{\phi \mu c_{t} }}}$$
(91)
$$k = 1.697\sqrt[5]{{\frac{{\phi \mu c_{t} (k_{f} w_{f} )^{4} }}{{x_{f}^{2} t_{LBLi} }}}}$$
(92)
$$k = 11.7166\sqrt[7]{{\left( {\frac{{\phi \mu c_{t} }}{{t_{USBLi} }}} \right)^{3} x_{f}^{2} (k_{f} w_{f} )^{4} }}$$
(93)
$$k = 6.741\sqrt[3]{{\frac{{\phi \mu c_{t} (k_{f} w_{f} )^{2} }}{{t_{LRi} }}}}$$
(94)
$$k = \frac{{670.606\phi \mu c_{t} k_{f} w_{f} x_{f} }}{{t_{USRi} }}$$
(95)
$$x_{f} = \frac{1}{81.6822}\sqrt {\frac{{kt_{BLRi} }}{{\phi \mu c_{t} }}}$$
(96)

Equations 8792 apply to either oil or gas wells.

Transient rate analysis by the conventional method for gas reservoirs

Equation 58 can be rewritten as:

$$1/q_{D} = 1.12\pi \sqrt {\frac{{t_{D} }}{{C_{fD}^{2} }}}$$
(97)
$$1/q = \frac{81.36T}{{kh[\Delta m(P)]k_{f} w_{f} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}} \sqrt t$$
(98)

Which can be expressed as Eq. 71 . A Cartesian plot of 1/q vs. t0.5 yields a straight line which slope, mL, provides the fracture conductivity,

$$k_{f} w_{f} = \frac{81.36T}{{{\text{kh}}[\Delta m(P)]m_{{{\text{BL}}}} }}\sqrt {\frac{k}{{\phi \mu c_{t} }}}$$
(99)

Once the dimensionless parameters have been replaced into Eq. 62, this be written as:

$$\frac{1}{{q_{{{\text{sc}}}} }} = \frac{T}{{5378.824{\text{kh}}[\Delta m(P)]k_{f} w_{f} }}\left( {\frac{k}{{\phi \mu c_{t} }}} \right)^{2} t^{2}$$
(100)

Or,

$$\frac{1}{q} = m_{{{\text{US}}}} \;t$$
(101)

A Cartesian plot of 1/q vs. t2 yields a straight line which slope, mUS, provides the half fracture length,

$$k_{f} w_{f} = \frac{T}{{5378.824{\text{kh}}[\Delta m(P)]m_{{{\text{US}}}} }}\left( {\frac{k}{{\phi \mu c_{t} }}} \right)^{2}$$
(102)

Equation 67 is rewritten as:

$$\frac{1}{q} = \frac{315.18T}{{{\text{kh}}\left[ {\Delta m\left( P \right)} \right]}}\sqrt[4]{{\frac{k}{{\phi (\mu c_{t} )_{i} x_{f}^{2} }}}}\sqrt[4]{{t_{{{\text{BL}}}} }}$$
(103)

Which can be written as Eq. 78. A Cartesian plot of 1/q vs. t0.25 yields a straight line which slope, mBL, provides the half fracture length,

$$x_{f} = 99337.95\left( {\frac{T}{{kh\left[ {\Delta m\left( P \right)} \right]m_{BL} }}} \right)^{2} \sqrt[{}]{{\frac{k}{{\phi (\mu c_{t} )_{i} }}}}$$
(104)

From Eq. 53, Tiab (1995), found an expression to find reservoir permeability,

$$k = \frac{70.6q\mu B}{{h(t*\Delta P^{\prime})_{r} }}$$
(105)

For gas wells,

$$k = \frac{711,26qT}{{h\left[ {t_{a} (p)*\Delta m^{\prime}(p)} \right]_{r} }}$$
(106)

According to Tiab (1995) the mechanical skin factor can be found from:

$$S_{m} = 0.5\left( {\frac{{\Delta P_{r} }}{{(t*\Delta P^{\prime})_{r} }} - \ln \left[ {\frac{{k\;t_{r} }}{{\phi \mu c_{t} r_{w}^{2} }}} \right] + 7.43} \right)$$
(107)

Or for gas wells, the apparent skin factor is given by Escobar,

$$S_{m} ^{\prime} = 0.5\left[ {\frac{{(\Delta m(P))_{r} }}{{\left[ {t*\Delta m(P)^{\prime}} \right]_{r} }} - \ln \left( {\frac{{{\text{kt}}_{r} }}{{\phi (\mu c_{t} )_{i} r_{w}^{2} }}} \right) + 7.4316} \right]$$
(108)

Escobar et al. (2016) presented a couple of expressions to find the mechanical skin factor when radial flow regime is absent. Also, a correlation to find half-fractured length as a function of fracture conductivity of vice versa.

Data interpretation

Example 1: Synthetic example

The second column of Table 1 presents the input data used to generate the simulated drawdown reported in Fig. 9.

Table 1 Input data for worked examples
Fig. 9
figure 9

Pressure and pressure derivative versus time log–log plot for synthetic example 1

TDS Technique

From Fig. 9 the following information was read:

$$t_{L} \, = \,0.091 h\Delta P_{L} = 0.811 \, psi(t*\Delta P^{^{\prime}} )_{L} = 0.405 \, psi$$
$$t_{US} \, = \,552.84 h\Delta P_{US} \, = \,250.03 \, psi \, (t*\Delta P^{^{\prime}} )_{US} \, = \,217.46 \, psi$$
$$t_{LUSi} \, = \,9 h$$

The linear flow regime was used to calculate the fracture conductivity by means of Eq. 19 and 21 using pressure derivative and pressure data, respectively:

$$\begin{gathered} k_{f} w_{f} = \frac{4.3(101)(0.45)(1.507)}{{(42)(0.405)}} \hfill \\ \sqrt {\frac{0.76(0.091)}{{0.08(0.45)(1.77 \times 10^{ - 5} )}}} = 5704.3\;{\text{md - ft}} \hfill \\ \end{gathered}$$
$$\begin{gathered} k_{f} w_{f} = \frac{8.6(101)(0.45)(1.507)}{{(42)(0.811)}} \hfill \\ \sqrt {\frac{0.76(0.091)}{{0.08(0.45)(1.77 \times 10^{ - 5} )}}} = 5697.25\;{\text{md - ft}} \hfill \\ \end{gathered}$$

The half-fracture length is then found from the unit-slope trend on the pressure derivative curve by using Eq. 24 and 25,

$$x_{f} = \frac{0.1(101)(0.76)(1.507)}{{(42)(0.08)(1.77 \times 10^{ - 5} )(6209.7)}}\left( {\frac{552.84}{{217.46}}} \right) = 86.7\;{\text{ft}}$$

Use of Eq. 49 leads to verify the half-fracture length:

$$x_{f} = \frac{1}{41.054}\sqrt {\frac{0.76(9)}{{0.08(0.45)(1.77 \times 10^{ - 5} )}}} = 79.8\;{\text{ft}}$$

Conventional analysis

Figures 10 and 11 present Cartesian plots from which the slope of the linear flow regime, mL, and the unit-slope behavior, mUS, were found to be:

$$x_{f} = \frac{0.1(101)(0.76)(1.507)}{{(42)(0.08)(1.77 \times 10^{ - 5} )(6209.7)}}\left( {\frac{552.84}{{250.03}}} \right) = 75.5\;{\text{ft}}$$
Fig. 10
figure 10

Cartesian plot of pressure drop versus the square root of time for synthetic example 1

Fig. 11
figure 11

Cartesian plot of pressure drop versus time for synthetic example 1

mL = 2.677 hr0.5/psi and mUS = 0.3229 h/psi. Use of Eqs. 33 and 37 allows finding the fracture conductivity,

$$\begin{gathered} k_{f} w_{f} = \frac{8.6(101)(0.45)(1.507)}{{42(2.677)}} \hfill \\ \sqrt {\frac{0.76}{{0.08(0.76)1.77 \times 10^{ - 5} )}}} = 5721.52\;{\text{md - ft}} \hfill \\ \end{gathered}$$
$$x_{f} = \frac{0.1(101)(1.507)(0.76)}{{(42)(0.08)(1.77 \times 10^{ - 5} )(5403)(0.3286)}} = 94\;{\text{ft}}$$

Example 2: Synthetic example

The third column of Table 1 were used to generate the reciprocal rate and reciprocal rate derivative data presented in Fig. 12.

Fig. 12
figure 12

Reciprocal rate and reciprocal rate derivative versus time log–log plot for synthetic example 2

Solution to this problem has two limitations which can be drawn from Fig. 12. During the unit-slope behavior the slope of the reciprocal rate is different than 2.0, therefore it is not recommended to use conventional analysis, Eq. 67 nor TDS Technique with the reciprocal rate value Eq. 63. Conventional analysis will only be applied to the linear flow regime.

TDS Technique

The following information was read from Fig. 12.

$$t_{L} \, = \,0.0871 h(1/q)_{L} = 0.00276 \, 1/STB$$
$$[t*(1/q)^{^{\prime}} ]_{L} = 0.00008797 \, 1/STBt_{LUSi} \, = \,9 h$$
$$t_{US} \, = \,66.082 h[t*(1/q)^{^{\prime}} ]_{US} \, = \,0.00845 \, 1/STB$$

From the linear flow regime fracture conductivity was estimated using Eq. 57 and 59 to be 1646.7 md-ft and 1657.92 md-ft respectively, as indicated below:

$$\begin{gathered} k_{f} w_{f} = \frac{4.034(250)(0.74)(1.1)}{{(36)(0.00008797)}} \hfill \\ \sqrt {\frac{0.15(0.0871)}{{(0.07)(0.74)(1.6 \times 10^{ - 6} )}}} = 1646.7\;{\text{md - ft}} \hfill \\ \end{gathered}$$
$$\begin{gathered} k_{f} w_{f} = \frac{12.75(250)(0.74)(1.1)}{{(36)(0.00845)}} \hfill \\ \sqrt {\frac{0.15(0.0871)}{{(0.07)(0.74)(1.6 \times 10^{ - 6} )}}} = 1657.9\;{\text{md - ft}} \hfill \\ \end{gathered}$$

The reciprocal rate derivative during the unit-slope behavior was used to find the half-fracture length with Eq. 61 to be 81.9 ft.

$$x_{f} = \frac{0.1053(250)(0.15)(1.1)}{{(36)(0.07)(1.6 \times 10^{ - 6} )(1616.8)}}\left( {\frac{66.082}{{0.00845}}} \right) = 81.9\;{\text{ft}}$$

Also, Eq. 87 allows us to verify the half-fracture length:

$$x_{f} = \frac{1}{48.576}\sqrt {\frac{0.15(9)}{{(0.07)(0.74)(1.6 \times 10^{ - 6} )}}} = 83.1\;{\text{ft}}$$

Conventional analysis

As commented previously, conventional analysis was applied only to the linear flow regime. A Cartesian plot of the reciprocal rate versus the square root of time is shown in Fig. 13 which slope mL = 0.0009336 hr0.5/STB. With this the fracture conductivity is computed with Eq. 73 to be 1661.76 md-ft.

$$\begin{gathered} k_{f} w_{f} = \frac{12.75(250)(0.74)(1.1)}{{(36)(0.0009336)}} \hfill \\ \sqrt {\frac{0.15}{{(0.07)(0.74)(1.6 \times 10^{ - 6} )}}} = 1661.8\;{\text{md - ft}} \hfill \\ \end{gathered}$$
Fig. 13
figure 13

Cartesian plot of the reciprocal rate versus the square root of time for synthetic example 2

Example 3: Field example

The fourth column of Table 1 presents information of a massive hydraulically fractured gas well reported by Agarwal, Carter, and Pollock (1979) which time and reciprocal rate are given in Table 2. The reciprocal rate derivative was estimated in this work and it is also reported in Table 2.

Table 2 Input data for worked examples

TDS technique

The reciprocal rate and reciprocal rate derivative versus time log–log plot is provided in Fig. 14. Only bilinear flow is observed and the following information is obtained from there:

$$t_{{{\text{BL}}}} = \, 6000{\text{ hr }}[t*(1/q)^{\prime } ]_{{{\text{BL}}}} = \, 0.002014 \, 1/{\text{STB}}$$
Fig. 14
figure 14

Reciprocal rate and reciprocal rate derivative versus time for field example

Reciprocal rate value during bilinear time is not read since the difference between the reciprocal rate derivative and the reciprocal rate derivative has a difference less than 4 units. Use of Eq. 84 allows to find the half-fracture length:

$$\begin{gathered} x_{f} = 6208.565\left( {\frac{720}{{0.0081(32)(3.96 \times 10^{8} )[0.002015]}}} \right)^{2} \hfill \\ \sqrt {\frac{0.0081(6000)}{{0.107(0.0176)(2.34 \times 10^{ - 4} )}}} = 790.15\;{\text{ft}} \hfill \\ \end{gathered}$$

Conventional analysis

The Cartesian plot of the reciprocal rate versus the fourth root of time are given versus time is given in Fig. 15. The slope of this plot:

$$m_{{{\text{BL}}}} = \, 0.00086{\text{ hr}}^{0.25} /{\text{STB}}$$
Fig. 15
figure 15

Reciprocal rate versus the fourth root of time for field example (conventional analysis)

The half-fracture length is found from Eq. 100 to be,

$$\begin{gathered} x_{f} = 99337.95\left( {\frac{720}{{0.0081(32)(3.96 \times 10^{8} )(0.00086)}}} \right)^{2} \hfill \\ \sqrt {\frac{0.0081}{{0.107(0.0176)(2.34 \times 10^{ - 4} )}}} = 896\;{\text{ft}} \hfill \\ \end{gathered}$$

Comments on the results

Table 3 summarizes the results of the worked examples. As expected, the results obtained from both TDS Technique and conventional analysis are satisfactory with errors ranging from 0.03% up to almost 10% for the synthetic examples and error of 8% for the estimation of the half-fracture length in the actual field case example. However, notice that the value of xf = 727 ft reported in Table 1 was obtained by using type-curve matching which is not the best well test interpretation method.

Table 3 Summary of results

It is customary to obtain fracture conductivity from the classical bilinear flow regime and half-fracture length from the classical linear flow regime. For our case, the trilinear model, the situation is inverted.

Conclusions

  1. a)

    TDS technique has been proven to be accurate and also to be a practical tool for well test analysis. In this case, it is extended to a trilinear flow model in oil and gas hydraulically finite-conductivity fractured wells under both constant-rate and constant-pressure cases. Several equations were developed and successfully tested against synthetic examples and a field case. Also, equations for the conventional analysis method are developed and verified with the mentioned examples. However, for the case of reciprocal rate during the intermediate unit-slope behavior, conventional analysis is not recommended to be applied if the slope of the reciprocal rate versus time log–log plot is different than two.

  2. b)

    The bilinear flow regime and the linear flow regime allow to obtained either half-fracture length and fracture conductivity, respectively. This situation is contrary to the classical bilinear and linear flow regimes.

Ethical Statement

This manuscript has not be submitted to other journal for simultaneous consideration.

The submitted work is original and has not been published elsewhere in any form or language (partially or in full.

This study is not split up into several parts to increase the quantity of submissions and submitted to various journals or to one journal over time.

Results are presented clearly, honestly, and without fabrication, falsification or inappropriate data manipulation (including image based manipulation). Authors adhere ourselves to discipline-specific rules for acquiring, selecting and processing data.

No data, text, or theories by others are presented as if they were the author’s own (‘plagiarism’). Proper acknowledgements to other works must be given (this includes material that is closely copied (near verbatim), summarized and/or paraphrased), quotation marks (to indicate words taken from another source) are used for verbatim copying of material, and permissions secured for material that is copyrighted.