Pressure data for field example (taken from [14]).

## Abstract

This chapter focuses on the application of Tiab’s direct synthesis (TDS) technique for practical and accurate interpretation of pressure tests on vertical wells in conventional reservoirs, so bilinear, linear, and elliptical flow regimes can be used for fracture characterization. Most fractured well interpretation tests are conducted using nonlinear regression analysis if the pressure model is available. This method has some drawbacks associated with the nonuniqueness of the solution. Also, the conventional straight-line method requires one plot for each individual flow regime observed in the pressure tests, and the estimated parameters cannot be verified. Tiab’s direct synthesis (TDS) methodology, which uses specific lines and intersection points found on the pressure and pressure derivative plot, is used in some direct equations which are obtained from the solution of the diffusivity equation for a given flow regime. It has been proven to provide accurate results, and its power allows verification of most results which is not possible from any other technique. The methodology has been successfully explained and tested by its application in two examples, although there exists more than a hundred articles that provide many useful applications.

### Keywords

- bilinear flow
- linear flow
- elliptical flow
- half-length fracture
- fracture conductivity
- hydraulic fracturing

## 1. Introduction

Throughout their history, well test analyses for fractured wells have received many contributions. For practical purposes, let us name the most important ones for this chapter. A good place to start is by mentioning the work in [1], which described the pressure behavior for infinite-conductivity and uniform-flux fractured wells, so people started conducting interpretation tests on such wells by using type-curve matching. Later, [2] introduced the concept of finite-conductivity fractures and established the onset value of dimensionless conductivity as 300. Values lower than that are considered finite-conductivity values, and those above 300 are classified as infinite conductivity. In [2], a fine semi-analytical solution was introduced for describing the well-pressure behavior in hydraulically fractured wells. This solution was then applied in [3] to provide a well interpretation method using type-curve matching. Since then, other mathematical solutions have been presented for finite-conductivity fractures. Among them, the work in [4] using fractal theory is worth mentioning.

The way of conducting well test interpretation was changed by the introduction of Tiab’s direct synthesis (TDS) technique by [5]. This revolutionary and modern technique focuses on the different flow regimes seen on the pressure derivative curve. Defined lines are drawn through each individual flow regime, and the intersection points found among them are read and used for reservoir characterization. Additionally, reading arbitrary points on the pressure and pressure derivative of each flow regime also serve for reservoir parameter determination. A great number of applications of the * TDS* technique are given in [6]. The second work [7], by the same author of [5], presented

*technique for infinite-conductivity and uniform-flux fractures in vertical wells. In [7], the elliptical or biradial flow regime was introduced and characterized. This elliptical flow is also seen in horizontal wells and was characterized in [8, 9, 10]. Because of the similarity between the mathematical models of hydraulic fractures and horizontal wells, this concept was applied by [11] to determine the average reservoir pressure in formations drained by horizontal wells using the*TDS

*technique. The infinite-conductivity model in [7] also included the late-time pseudosteady-state period as well as some equations involved in the drainage area (conventional analysis for this case was included in [12]). This may be disadvantageous for inexperienced users of*TDS

*technique when interpreting pressure tests without reaching reservoir boundaries because the equations involved the use of the unknown reservoir drainage area, although it can be still applied by using the intersection points. To overcome this drawback, [13] presented a new mathematical model excluding the late-time pseudosteady-state period.*TDS

* TDS* technique for finite-conductivity fractured wells is given in [14], with practical field applications to demonstrate the usefulness of the technique. The fracture parameters can be readily obtained by using an arbitrary point on the flow regimes.

*technique plays an important role when analyzing short pressure tests because a user can “make up” nonexisting flow regimes since, for instance, the radial flow horizontal line can be obtained from the reservoir permeability even though radial flow regime is absent. [15, 16] extended the works of finite- and infinite-conductivity fractures in naturally fractured reservoirs. The equations provided by these works can also be applied to either homogeneous or naturally fractured formations since they involve a dummy variable that takes the value of one for the homogeneous case or the value of the dimensionless storativity coefficient for the case of a naturally fractured formation.*TDS

* TDS* technique has also been extended to several scenarios related to hydraulically fractured wells. For instance, when a finite-conductivity fracture intersects with a fault, the pressure trace changes; then, the equations developed in [17] apply for this case. There are cases where a threshold pressure is required to start the flow. The work in [18] includes this concept in uniform-fractured vertical wells, and the work in [19] includes the concept for horizontal wells. Also, when the fractured face is damaged, a pseudolinear flow regime develops along the fracture. [1] included

*technique to characterize such systems. [16] presented*TDS

*technique for fractured wells in gas composite reservoirs.*TDS

*technique can also be usefully applied to transient-rate analysis, as seen in [20]. Application of*TDS

*technique to horizontally isolated fractured wells was presented and characterized in [21] and in conventional analysis in [22]. The works in [23, 24] use*TDS

*technique for shale reservoirs. Other applications of*TDS

*technique to these systems are given by [25] under transient-rate analysis and [26] for pressure-transient analysis conditions. Other important applications of TDS Technique to fractured wells are given by [29, 30].*TDS

This chapter is devoted to the application of * TDS* technique to hydraulically fractured wells in either homogeneous or naturally fractured formations. Without given detailed derivations, the expressions for characterizing the hydraulic fracture parameters are presented along with the way they should be used. Important relationships and practical exercises are included.

## 2. TDS basis

The pioneer publication on the * TDS* technique, [5], explains in detail the derivation of the equations. The Laplace space solution of the arithmetic pressure derivative for a homogeneous and infinite reservoir with skin and wellbore storage is also presented in [5] and given by

However, we know that the pressure derivative is a horizontal line during radial flow regime. The dimensionless pressure derivative during radial line is easier represented by

Then, to obtain practical equations, dimensionless parameters must be used. The dimensionless time, based upon half-fracture length and reservoir drainage area, is given below:

and

The dimensionless pressure and pressure derivative parameters for oil reservoirs are given by

and

Finally, the dimensionless fracture conductivity introduced in [3] is defined as

It is observed from Eq. (5) that the two key parameters of a hydraulic fracture are the half-fracture length, * x*, and the fracture conductivity,

_{f}

*. The total length of the fracture is given by 2*k

_{f}w

_{f}

*.*x

_{f}

The easiest application of * TDS* technique is given by replacing the dimensionless pressure derivative defined by Eqs. (6) and (2), to provide an expression to readily determine formation permeability:

where (* t**∆

*′)*P

*is the pressure derivative value during radial flow regime. The equations for the*

_{R}

*technique are derived in the same manner Eq. (8) was obtained.*TDS

## 3. Biradial flow regime

Biradial or elliptical flow normally results in a hydraulically fractured well when areal anisotropy is present. This is recognized on the pressure derivative versus time log-log plot by a straight line with a slope of 0.36. In hydraulic fractures, the flow from the formation to the fracture is described by parallel flow lines resulting in a linear flow geometry better known as linear flow regime and characterized by a slope of 1/2 on the pressure derivative versus time log-log plot.

Both linear flow and biradial/elliptical flow regimes are seen on the plot of dimensionless pressure and pressure derivative versus dimensionless time based on half-fracture length for a naturally fractured formation. New expressions for the elliptical flow regime introduced in [13] excluding reservoir drainage area are given by.

and

being * ξ* a dummy variable that defines either a homogeneous or naturally fractured formation. When

*= 1, a homogeneous reservoir is considered. For the case of naturally fractured formations,*ξ

*=*ξ

*, the dimensionless storativity coefficient.*ω

Once dimensionless parameters given by Eqs. (3), (5), and (6) are replaced into Eqs. (9) and (10), respectively, and solve for the half-fracture length, which yields

and

TDS technique is based on drawing a straight line throughout a given flow regime; then, the user is expected to read the pressure, Δ* P*, and pressure derivative, (

_{BR}

**Δ*t

*’)*P

*, at a given time,*

_{BR}

*. A better way to reduce noise effects consists of extrapolating the mentioned straight line (biradial for this case) to the time of 1 h and read the pressure derivative value, (*t

_{BR}

**Δ*t

*’)*P

_{BR1}, at 1 h. For this case, the pressure and pressure derivative set in Eqs. (11) and (12) is changed to Δ

P

_{BR1}and (

**Δ*t

*’)*P

_{BR1}, respectively.

When bilinear flow is unseen, fracture conductivity can be found with an expression presented in [27]

[5] also provided an equation for the determination of the skin factor using an arbitrary point read during radial flow regime:

The pseudosteady-state regime governing the pressure derivative equation is given by

[7] used the point of intersection, * t*, of Eqs. (2) and (15) to derive an equation for the estimation of the drainage area:

_{RPi}

The derivation of Eq. (16) follows a similar idea as that presented later in Section 4 for the use of the points of intersection.

## 4. Bilinear and linear flow regimes

Bilinear flow regime takes place when two linear flows, normal one flowing into the other, take place simultaneously. This situation occurs in low conductivity fractures where linear flow along the fracture and linear flow from the formation to the fracture are observed. Bilinear flow is recognized in the pressure derivative curve by a slope of 0.25. However, this is not shown in Figure 1 since bilinear flow is absent. The governing expressions for early bilinear and linear flow regimes for vertical fractures in naturally fractured systems were, respectively, presented in [16]

and

Linear flow regime can be used to find the half-fracture length, and bilinear flow regime allows finding the fracture conductivity. Once the dimensionless quantities of Eqs. (1) and (3)–(5) are replaced in Eqs. (16)–(19), the fracture conductivity is solved for then

Once the fracture conductivity is found, Eq. (7) applies to find the dimensionless fracture conductivity if reservoir permeability and the half-fracture length are known. When bilinear flow is absent, the fracture conductivity may be found from Eq. (13), or the dimensionless fracture conductivity can be read from Figure 2:

and

## 5. Points of intersection

If bilinear flow also takes place, then the point of intersection between the pressure derivatives of the bilinear and biradial flow lines, * t*, given by Eqs. (10) and (18), respectively, allows the development of an equation to find the half-fracture as follows:

_{BLBRi}

Simplifying,

Replacing the dimensionless quantities, Eqs. (3) and (7) in Eq. (26) lead to

Solving for the half-fracture from Eq. (27), we readily obtain

By the same token, the intercept of Eq. (20) with Eq. (18), * t*, provides another expression to find the half-fracture length:

_{LBRi}

Bilinear flow regime is absent in the plot of Figure 1. Linear, biradial, and radial flow regimes along with the late pseudosteady-state period are seen. The interception points formed by the possible combinations of such periods can be represented schematically in this plot.

Another way to find the half-fracture length comes from the intersection of Eqs. (2) and (10), * t*, and Eq. (10) with Eq. (15), so that

_{RBRi}

and

The intercept point resulting between linear flow and bilinear flow lines given by the governing pressure derivative solutions, Eqs. (18) and (19), can be used to find either half-fracture length or permeability:

* t* is the intersection of the bilinear pressure derivative line given by Eq. (18) with the radial flow regime line (Eq. (2)). This intersection point serves as the estimation of either permeability or fracture conductivity:

_{BLRi}

## 6. Other estimations

The expressions for determination of the naturally fractured reservoir parameters cannot be included in this chapter for space reasons. However, they can be found in [15, 16], which also used intersection points and maximum and minimum data read from the pressure and pressure derivative curve.

Radial flow regime may be absent in short tests run in fractured wells with the sole purpose of determining fractured parameters. For these cases, the skin factor can be estimated from any of the two empirical correlations presented by [27]

and

where

Additionally, fracture conductivity can be read from the plot given in Figure 2.

Finally, space reasons prevent including * TDS* technique for fractured wells in unconventional shale formations. The reader is referred to [23, 24, 25, 26].

## 7. Examples

### 7.1 Field example

[14] presented a field example of a fractured well test. Pressure and pressure derivative data are given in Table 1 and Figure 3. Other relevant data are provided below:

, h | Δ, psia | *Δ’, psia | , h | Δ, psia | *Δ’, psia |
---|---|---|---|---|---|

0.23 | 102 | 26.3 | 15 | 390 | 117 |

0.39 | 115 | 30 | 20 | 423 | 112 |

0.6 | 130 | 35.8 | 25 | 446 | 120 |

1 | 145 | 40.8 | 30 | 471 | 141 |

1.8 | 183 | 57.2 | 35 | 493 | 136.5 |

2.4 | 195 | 67 | 40 | 510 | 132 |

3.8 | 260 | 83.3 | 45 | 526 | 135 |

4.1 | 265 | 69.2 | 50 | 540 | 150 |

4.96 | 280 | 96.9 | 55 | 556 | 137.5 |

6.2 | 308 | 102.3 | 60 | 565 | 144 |

8.5 | 320 | 103.3 | 65 | 580 | 121.1 |

10 | 345 | 149 | 71 | 583 |

Using a commercial well test software, the following parameters were estimated by nonlinear regression analysis:

The objective is to compute the hydraulic fracture parameters using the * TDS* technique and compare results obtained from the regression analysis.

#### 7.1.1 Solution

7.1.1.1 Step 1: Obtain the characteristic points

Once the pressure and pressure derivative versus time log-log plot is built and reported in Figure 3, the characteristic points are read from such plot as follows:

7.1.1.2 Step 2: Estimate permeability and skin factor

Permeability and skin factor are found in Eqs. (8) and (14) to be 0.76 md and −4.68, respectively.

7.1.1.3 Step 3: Estimate fracture conductivity

Fracture conductivity is estimated using Eqs. (21) and (22):

From Figure 3, * t* = 200 hr. A very close value is obtained from Eq. (33):

_{BLRi}

which indicates that the calculation of the fracture conductivity is accurate. Notice that instead of estimating * t* the fracture conductivity can be found instead to obtain another value of fracture conductivity; then, Eq. (33) can also be expressed as

_{BLRi}

7.1.1.4 Step 4: Half-fractured length and dimensionless fracture conductivity estimation

Find half-fracture length with Eqs. (23) and (24):

Solve for half-fracture length from Eq. (13) and find this:

Find the dimensionless fracture conductivity using Eq. (5):

The above value confirms that the fracture has finite conductivity.

### 7.2 Synthetic example

[13] presented a synthetic example of a pressure test run in a bounded homogeneous reservoir with the information given below:

Estimate the half-fracture length by the * TDS* technique, and compare the answer with the value used for generating the test.

#### 7.2.1 Solution

7.2.1.1 Step 1: Obtain the characteristic points.

A pressure and pressure derivative versus time log–log plot is presented in Figure 4, from which the following characteristic points are read:

## 8. Comments on the results

The two given examples show three aspects of the * TDS* technique: (1) practical use, (2) accuracy, and (3) self-confirmation.

As shown in the exercises, the process includes defining flow regimes, drawing a few lines, and finally computing the necessary parameters. Contrary to the conventional straight-line method, which requires a plot for each flow regime, TDS technique uses only the pressure and pressure derivative versus time log-log plot. Computations are straight forward.

Table 2 summarizes the main parameters obtained in the two worked examples. The results show a good agreement between the calculated results by * TDS* technique and the results obtained from commercial software packages, for the field case. The results of the half-fracture length for the synthetic case using

*technique are even better compared to the input value use to simulate the test. This demonstrates that*TDS

*technique is an accurate methodology which has been also presented in many publications, not only in the list reference but also in others not mentioned here.*TDS

The last aspect dealt with is self-confirmation. In the field example, three values of half-fracture length and three values of fractured conductivity were found, and for the synthetic example, two values of half-fracture length were estimated from different equations. All the estimations match with the reference values.

## 9. Conclusion

It has been shown that * TDS* technique is a powerful, practical, and accurate tool for well test interpretation because manipulations are easy to do and parameters can be confirmed from different sources from the same pressure test. Compared to reference values, the worked examples provided accurate results of both half-fracture length and hydraulic fracture conductivity. Besides being accurate,

*technique has the great advantage of being able to estimate a given parameter, such as half-fracture length or fracture conductivity, from more than one source or equation. This provides a means of verifying that the estimated parameter is in a good range.*TDS

## Acknowledgments

The author wishes to express his gratitude to Universidad Surcolombiana for providing him the time to write this chapter.

## Nomenclature

A | Draining area (ft2) |

B | Oil volume factor (rb/STB) |

CfD | Dimensionless fracture conductivity |

ct | Compressibility (1/psi) |

h | Formation thickness (ft) |

k | Formation permeability (md) |

kfwf | Fracture conductivity (md-ft) |

P | Pressure (psi) |

Pwf | Well-flowing pressure (psi) |

q | Oil flow rate (STB/D) |

qg | Gas flow rate (MSCF/D) |

rw | Wellbore radius (ft) |

xf | Half-fracture length (ft) |

s | Skin factor |

t | Test time (h) |

tp | Production time (h) |

t*∆P′ | Pressure derivative (psi) |

tD*PD’ | Dimensionless pressure derivative |

Greek symbols | |

∆ | Change |

ϕ | Porosity (fraction) |

λ | Interporosity flow parameter |

μ | Viscosity (cp) |

ξ | Variable to identify homogeneous (ξ = 1) or heterogeneous (ξ = ω) reservoirs |

ω | Dimensionless storativity coefficient |

Suffixes | |

BL | Bilinear |

BL1 | Bilinear at 1 h |

BLL | Bilinear-linear intersection |

BR | Birradial |

BR1 | Birradial at 1 h |

BRBLi | Birradial-bilinear intersection |

BRPi | Birradial-pseudosteady intersection |

D | Dimensionless |

DA | Dimensionless based on area |

Dxf | Dimensionless based on half-fractured length |

DLBRi | Dual linear-birradial intersection |

LBRi | Linear-birradial intersection |

R | Radial |

RBRi | Radial-birradial intersection |

RPi | Intersect of radial-pseudosteady-state lines |

w | Well |

t | Time |

P | Pseudosteady state |

## References

- 1.
Gringarten AC, Ramey HJ, Raghavan R. Applied pressure analysis for fractured wells. Journal of Petroleum Technology. Richardson, Texas, USA. 1975:793-800. DOI: 10.2118/5496-PA - 2.
Cinco LH, Samaniego VF, Dominguez AN. Transient pressure behavior for a well with a finite-conductivity vertical fracture. SPE Journal. 1978; 18 (4):253-264. DOI: 10.2118/6014-PA - 3.
Cinco-Ley H, Samaniego VF. Transient pressure analysis for fractured wells. Journal of Petroleum Technology. 1981; Sept. :1479-1766. DOI: 10.2118/7490-PA - 4.
Cossio M. A Semi-Analytic Solution for Flow in Finite-Conductivity Vertical Fractures Using Fractal Theory. Richardson, Texas, USA: Society of Petroleum Engineers; 2012. DOI: 10.2118/163057-STU - 5.
Tiab D. AAnalysis of pressure and pressure derivative without type-curve Matching: 1-skin and wellbore storage. Journal of Petroleum Science and Engineering. 1995; 12 :171-181. Also Tiab D. Analysis of pressure and pressure derivatives without type-curve matching: i-skin and wellbore storage. Society of Petroleum Engineers. 1993. DOI: 10.2118/25426-MS - 6.
Escobar FH. Recent Advances in Practical Applied Well Test Analysis. New York: Nova Science Publishers, Inc.; 2015. 422p - 7.
Tiab D. Analysis of pressure derivative without type-curve matching: Vertically fractured wells in closed systems. Journal of Petroleum Science and Engineering. Richardson, Texas, USA. 1994; 11 :323-333. DOI: 10.2118/26138-MS. This paper was originally presented as Society of Petroleum Engineers - 8.
Escobar FH, Muñoz OF, Sepulveda JA. Horizontal permeability determination from the elliptical flow regime for horizontal wells. CT&F – Ciencia, Tecnología y Futuro. 2004; 2 (5):83-95. ISSN 0122-5383 - 9.
Escobar FH, Montealegre M. Conventional Analysis for the Determination of the Horizontal Permeability from the Elliptical Flow of Horizontal Wells. Richardson, Texas, USA: Society of Petroleum Engineers; 2007. DOI: 10.2118/105928-MS - 10.
Escobar FH, Montealegre-M M. Determination of horizontal permeability from the elliptical flow of horizontal wells using conventional analysis. Journal of Petroleum Science and Engineering. 2008; 61 :15-20. ISSN: 0920-4105 - 11.
Escobar FH, Cantillo JH, Santos N. A practical approach for the estimation of the average reservoir pressure from multi-rate tests in long horizontal wells. Fuentes: El Reventón Energético. 2011; 9 (1):13-20 - 12.
Escobar FH, Montealegre MM, Cantillo JH. Conventional analysis for characterization of bi-radial (elliptical) flow in infinite-conductivity vertical fractured wells. CT&F – Ciencia, Tecnología y Futuro. 2006; 3 (2):141-147 - 13.
Escobar FH, Ghisays-Ruiz A, Bonilla LF. New model for elliptical flow regime in hydraulically-fractured vertical wells in homogeneous and naturally-fractured systems. Journal of Engineering and Applied Sciences. 2014; 9 (9):1629-1636 - 14.
Tiab D, Azzougen A, Escobar FH, Berumen S. Analysis of Pressure Derivative Data of Finite-Conductivity Fractures by the “Direct Synthesis” Technique. Society of Petroleum Engineers; 1999. DOI: 10.2118/52201-MS - 15.
Escobar FH, Zhao YL, Fahes M. Characterization of the naturally fractured reservoir parameters in infinite-conductivity hydraulically-fractured vertical wells by transient pressure analysis. Journal of Engineering and Applied Science. 2015; 10 (12):5352-5362 - 16.
Tiab D, Bettam Y. Practical Interpretation of Pressure Tests of Hydraulically Fractured Wells in a Naturally Fractured Reservoir. Richardson, Texas, USA: Society of Petroleum Engineers; 2007. DOI: 10.2118/107013-MS - 17.
Escobar FH, Tiab D, Berumen-Campos S. Well Pressure Behavior of a Finite-Conductivity Fracture Intersecting a Finite Sealing-Fault. Richardson, Texas, USA: Society of Petroleum Engineers; 2003. DOI: 10.2118/80547-MS - 18.
Escobar FH, Zhao YL, Pournik M, Liu QG, Olaya-Marin G. Interpretation of pressure tests in uniform-flux fractured vertical wells with threshold pressure gradient in low permeability reservoirs. Journal of Engineering and Applied Science. 2015; 10 (20):9364-9372 - 19.
Escobar FH, Zhao YL, Zhang LH. Interpretation of pressure tests in horizontal wells in homogeneous and heterogeneous reservoirs with threshold pressure gradient. Journal of Engineering and Applied Sciences. 2014; 9 (11):2220-2228 - 20.
Escobar FH, Castro JR, Mosquera JS. Rate-transient analysis for hydraulically fractured vertical oil and gas wells. Journal of Engineering and Applied Science. 2014; 9 (5):739-749. ISSN: 1819-6608 - 21.
Al Rbeawi SJH, Tiab D. Effect of the Number and Length of Zonal Isolations on Pressure Behavior of Horizontal Wells. Society of Petroleum Engineers; 2011. DOI: 10.2118/142177-MS - 22.
Escobar FH, Meneses AR, Losada LM. Straight-line conventional transient pressure analysis for horizontal wells with isolated zones. Dynamis. 2014; 81 (185):78-85 - 23.
Bernal KM, Escobar FH, Ghisays-Ruiz A. Pressure and pressure derivative analysis for hydraulically-fractured shale formations using the concept of induced permeability field. Journal of Engineering and Applied Science. 2014; 9 (10):1952-1958 - 24.
Escobar FH, Montenegro LM, Bernal KM. Transient-rate analysis for hydraulically-fractured gas shale wells using the concept of induced permeability field. Journal of Engineering and Applied Science. 2014; 9 (8):1244-1254 - 25.
Escobar FH, Rojas JD, Ghisays-Ruiz A. Transient-rate analysis for hydraulically-fractured horizontal wells in naturally-fractured shale gas reservoirs. Journal of Engineering and Applied Science. 2015; 10 (1):102-114 - 26.
Escobar FH, Bernal KM, Olaya-Marin G. Pressure and pressure derivative analysis for fractured horizontal wells in unconventional shale reservoirs using dual-porosity models in the stimulated reservoir volume. Journal of Engineering and Applied Sciences. 2014; 9 (12):2650-2669 - 27.
Tiab D. Advances in Pressure Transient Analysis — TDS Technique. Lecture Notes Manual. Norman, Oklahoma, USA: The University of Oklahoma; 2003. 577p - 28.
Economides MJ, Watters, Dunn-Norman S. Petroleum Well Construction. New York: John Wiley & Sons; 1988. 622p - 29.
Escobar FH, Gonzalez RA, Hernandez LM, Hernandez CM. Pressure and pressure derivative analysis for hydraulically fractured vertical wells with face skin. Journal of Engineering and Applied Science. 2016; 11 (13):8268-8273 - 30.
Zhao YL, Escobar FH, Hernandez CM, Zhang CP. Performance analysis of a vertical well with a finite-conductivity fracture in gas composite reservoirs. Journal of Engineering and Applied Sciences. 2016; 11 (15):8992-9003