Well test, rate  transient analysis and reservoir simulation for characterizing multi-fractured unconventional oil and gas reservoirs

Nick Bahrami1, David Pena1, Ian Lusted1


Aurora Oil and Gas, Houston, USA





Unconventional  reservoirs  such  as  shale  gas and shale oil have become an increasingly important source of energy in the USA with potential reservoirs identified worldwide. Due to the insufficient permeability of the shale reservoirs, they require efficient stimulation using multistage  hydraulic fractures to  produce gas  in  commercial quantities. A critical challenge in the reservoirs is performance evaluation of the fracturing and characterization of the  stimulated  reservoir volume  (SRV) for permeability and hydraulic fracture size. Conventional well test analysis in multi-stage fractured shale reservoirs may not provide reliable results due to the extensive wellbore storage effect, fracture complexities, and heterogeneity of the low-permeability reservoir. To overcome such issues, advanced well test analysis techniques integrated with rate transient analysis can be used to reduce uncertainties associated with estimation   of   the   reservoir   and   hydraulic   fractures’ dynamic parameters. This paper proposes a practical methodology and  workflow for  characterizing  the  SRV parameters in multi-fractured wells in unconventional oil and gas reservoirs using well test and rate transient data analysis based on diffusivity equation solution for linear and elliptical flow regimes integrated with numerical reservoir simulation. A reservoir simulation model is built and run for a typical fractured shale reservoir to verify the reliability of the  proposed simplified approach. Furthermore, multi-fractured unconventional reservoir field examples of well test analysis, reservoir simulation and history matching are presented to show how the stimulated reservoir volume can be characterized to perform a more

reliable  production  forecast  in  shale  oil  and  shale  gas reservoirs.




Unconventional reservoirs such as shale gas and shale oil have become an increasingly important source of energy worldwide. Due to the insufficient permeability of the shale reservoirs, they  require  efficient stimulation  using hydraulic fractures to produce hydrocarbon in commercial quantities (Bagherian et al. 2010).

The propagation and direction of hydraulic fractures are influenced by many factors but are primarily controlled by in situ stresses. Where there is high contrast between minimum  and  maximum  horizontal stresses, the  stimulation creates a narrow or linear fracture fairway, and where the stress contrast is low, wide or complex fracture geometry is created during the treatment as shown in Fig. 1 (Fan et al. 2010). Considering a horizontal well in a normal faulting stress regime, if the horizontal well is drilled in the direction of maximum horizontal stress, the longitudinal hydraulic fractures are likely to be initiated along or parallel to the wellbore, and if the horizontal well is drilled in the direction of minimum horizontal stress, then the transverse hydraulic fractures are initiated perpendicular to the wellbore axis as shown in Fig. 2 (Hossain and Rahman 2008).

organic content (non-productive shale intervals). Also at the well location, although the well may penetrate a productive shale, but in the same interval, away from the well location there might be non-productive shale sections. In other words, in unconventional shale reservoirs, productive shale zones and non-productive shale sections have randomly been distributed as illustrated in Fig. 3.

Therefore, a common practice in unconventional reservoirs is to drill a horizontal well, complete the well with series of perforated intervals with an specific spacing and length  to control where the fracturing stages to be regardless of where the productive shale zones are at wellbore location, and then hydraulically fracture the whole formation along the wellbore in multi-stages to create a stimulated reservoir volume  (SRV) with  an  enhanced  effective  permeability around each hydraulic fracture that has gone through productive and non-productive shale sections.

In the case of a hydraulic fracture network around the wellbore where fracturing creates  a stimulated reservoir volume as shown in Fig. 3 (the typical scenario in shale reservoirs), the  main  SRV  parameters  are  the  reservoir thickness (h), SRV permeability (KSRV), fracture half- length size (Xf), fracture spacing (2 x Ye),  drainage area (A), and cross-section area (Ac). In the simplified analytical methods used in studying shale gas and shale oil reservoirs, it is assumed that the whole SRV is homogeneous productive shale. For the SRV shown in Fig. 3, the drainage area (A) and the cross-section area perpendicular to the flow (Ac) can be defined as follows, considering ‘n’ stages of hydraulic fractures (Arevalo et al. 2002):

The SRV is a network of hydraulic fractures and untreated  matrix  which  are  randomly  distributed in  the dual-porosity  and  dual-permeability  system.  The  main SRV parameters are SRV permeability, hydraulic fracture spacing  (shale  matrix  block  size)  and  SRV  porosity (volume  of  the  open  fractures  inside  the  stimulated reservoir volume). Characterizing the SRV parameters is a critical challenge as it indicates the effectiveness of the stimulation of the low-permeability shale, the efficiency of the drained volume due to well spacing, and the well future production performance.


One of the common methods for characterizing the dynamic reservoir parameters is well testing, in which the pressure transient data are recorded and analyzed using the plot of transient pressure and its derivative versus time on Log–Log scale. The test data are normally analyzed using a diagnostic plot to identify the radial flow regime and calculate the average permeability of the reservoir accordingly. In a conventional reservoir, wellbore storage effect may last few hours, and then in few days the reservoir will exhibit its average properties upon reaching infinite acting radial flow regime. However, for multi-fractured unconventional reservoirs, things are radically different because of the large drop in mobility, fracture complexities, and heterogeneity of the low-permeability formation that cause an extensive wellbore storage effect and slow propagating pressure pulse into the formation. As a result, we may have to wait months or years to detect the complicated flow regimes in the drainage area around the wellbore and fractures, and maybe centuries to eventually detect the equivalent infinite acting radial flow regime (KAPPA Engineering 2015). Therefore, using the conventional well test analysis to characterize the SRV parameters in multi-fractured shale reservoirs may fail to provide reliable results, and advanced techniques may be required.

Reservoir  flow regimes in hydraulically fractured wells



A pressure transient test breaks into several flow regimes, each seeing deeper in the reservoir than the last. Depending on the well completion type and the reservoir geological and geometric attributes, different flow regime might be revealed on pressure transient diagnostic plots. In conventional well test analysis, diagnosis of the radial flow regime is critical in quantitative well test interpretation, since  a  reliable  value for reservoir permeability  can  be estimated when the late-time radial flow regime is established in the reservoir (Badazhkov 2008; Bourdarot 1998).

In the case of a multi-fractured well as illustrated in ‘‘Appendix A’’ and shown in Fig. 4, the main reservoir flow regimes after wellbore storage effect, are the early-time  linear  flow regime  towards  the  hydraulic  fracture wings in SRV in the vicinity of the fractures, the early-time elliptical flow regime in the SRV towards the drainage area of linear flow perpendicular to each hydraulic fracture, and the early-time boundary-dominated flow when stimulated reservoir volume, SRV, around each hydraulic fracture is depleted. After the SRV depletion, hydrocarbon flow is provided from untreated shale around SRV (the untreated reservoir rock surrounding the stimulated reservoir volume acts as boundary).

After the reservoir flow regimes in the SRV, the SRV boundary-dominated flow effect is then followed by linear flow regime inside the untreated reservoir zone towards the drainage area of the multi-stage fractured well (Fig. 4c), then elliptical flow regime, and finally at late time when pressure disturbance propagates deep enough into the reservoir,  a  pseudo  radial  flow  regime  is  established (Fig. 4d), with slope of zero on pressure derivative. On the pressure derivative curve based on solution of radial flow diffusivity equation, the slope of +1 shows wellbore storage effect, the slopes of -0.5, +0.5, +0.25 and +0.36 indicate spherical, linear, bi-linear and elliptical flow regimes, respectively, and the slope of zero indicates radial flow regime.  In  typical  testing  time  duration  in  multifractured shale reservoirs, only the early-time flow regimes might be the ones that can be detected on the diagnostic plots, but not the late-time flow regimes.


In  multi-fractured  shale  reservoirs, field observations indicate  that  the  early-time  linear  flow regime  and  the early-time  elliptical  flow regime  are  established  in  the SRV, but radial flow regime cannot be established in the SRV. The reservoir flow regimes in SRV are followed by the SRV boundary-dominated effect. Well test analysis methods such as rate transient analysis (RTA) and pressure buildup analysis (PBA) in the multi-fractured wells require analyzing the reservoir flow regimes that are detected on well test analysis diagnostic plots, to characterize the SRV dynamic parameters.

Rate transient analysis (RTA) for multi-fractured shale reservoirs



In rate transient analysis for oil wells, the pressure data are normalized using the rate data, and then the RTA plots are made to characterize SRV parameters. For gas wells, bottom-hole pressure, P, should be converted to pseudo pressure, m(p), and then the pseudo pressure data are normalized using gas production rate. The normalized pressure data are used to make the following RTA plots for SRV characterization (Malallah et al. 2007):


  • A plot of the normalized pressure (d[Pwf]/d[Q]) versus time (d[t]) on the Log–Log scales as shown in Fig. 5a can identify the data points related the linear flow regime based on the line slope of  1/2 on the diagnostic plot. When the slope increases from  1/2   to higher values, this indicates the time at which the boundary-dominated flow effect starts (tBDF).  For gas wells, normalized pseudo pressure should be plotted versus time as shown in Fig. 5c.
  • A plot of the normalized pressure (d[Pwf]/d[Q]) versus (d[t1/2]) as shown in Fig. 5b can provide slope of the linear flow regime straight line (mLF). For gas wells, normalized pseudo pressure should be plotted versus the time function as shown in Fig. 5d.

The RTA equations use tBDF  and mLF to characterize the SRV parameters. For oil wells with constant production rate, the SRV permeability and hydraulic fracture size can be estimated as follows (Rasdi and Chu 2012):

For constant production rate in gas wells, the following equations have been proposed (Ibrahim et al. 2006):

By combining the Eqs. 1 and 2 with the Eqs. 5 and 6, the SRV permeability and hydraulic fracture size can be estimated as follows:

In the above equations, P is pressure (psia), t is time (days), Qo   is oil flow rate (STBD), Qg   is gas flow rate (MSCFD), B is formation volume factor, h is reservoir thickness (ft), μ is viscosity (cp), φSRV  is the effective average porosity of the porous media in SRV that contributes  to  the  flow  (fraction),  Ct   is  total compressibility, Ye is distance to the boundary (optimistically, Ye  is half distance between each two adjacent hydraulic fractures), L is horizontal well length (ft), KSRV  is permeability of the stimulated reservoir volume, T is reservoir temperature (R), and n is number of bottom-hole pressure data, then WHP data should be converted to BHP using multi-phase flow well model correlations, and then be input into the RTA models. RTA can practically be used in multi-fractured shale reservoirs to characterize the SRV parameters.

Pressure buildup  analysis (PBA) for multifractured shale reservoirs



In well test analysis of shale reservoirs, the pressure transient data and the derivatives should be plotted on the Log–Log diagnostic plot: d(p)/d(t0.5) as linear flow derivative, d(p)/d(t0.36) as elliptical flow derivative and d(p)/d(ln[t]) as radial flow derivative (see the ‘‘Appendix A’’ for more details). The plot is used to identify the pressure derivative values for infinite acting linear flow, elliptical flow and radial flow regimes (mRF) from zero slope line on radial flow derivative, mE11 from zero slope line on elliptical flow derivative and mLF  from zero slope line on linear derivative) as shown in Fig. 6 (Arevalo et al. 2002; Martinez et al. 2012). If the test duration is short, then on the Log–Log diagnostic plot only the linear flow regime might be observed, and if the testing time is long enough, then both the linear and elliptical flow might be observed. For the PBA method, the pressure transient data before the effect of the early-time boundary-dominated flow effect are used for SRV characterization depending on the detected flow regimes.

The  linear  flow and  elliptical  flow regimes  are  both function of KSRV and Xf, and combining the solutions of the diffusivity equation for linear and elliptical flow regimes and solving the two non-linear equations simultaneously can provide the two unknowns KSRV and Xf. The solution of the diffusivity equation for infinite acting elliptical flow as function of elliptical flow regime slope, mE11, has been proposed for horizontal oil and gas wells (Martinez et al. 2012). The modified and re-derived solution of elliptical the hydraulic fracturing stages.


The  RTA  method  requires  recording of  surface  production rates and bottom-hole pressure data. However, if only the well-head pressure data are available instead of flow equation and linear flow equation for multi-fractured horizontal wells are as follows:

Integrating the above two Eqs. 9 and 10 can provide the two unknowns Xf and KSRV for an oil well in shale reservoirs (assuming that the effective average porosity of SRV to be known).

For gas wells and using the pseudo pressure in well test analysis, the following equations can provide fracture size and SRV permeability:

In the above equations, P is pressure (psia), t is time (h), q is oil flow rate (STBD), B is formation volume factor, h is reservoir thickness (ft), μ is viscosity (cp), φSRV is the effective average porosity of SRV (fraction) that may be dominated by matrix porosity effect since the fractures may have very small porosity, Ct is total compressibility, n is number of hydraulic fractures, and KSRV  is SRV average permeability that may be mainly controlled by permeability of the fractures. Based on the analytical derivations, the value of mLF  can be a good indicator of well deliverability. The production data for the first few months can be displayed on RTA plots to evaluate well deliverability of producing wells in a field.


For multi-fractured wells in unconventional low-permeability gas reservoirs, only the linear flow regime may be observed on the diagnostic plots (testing time is not long enough to have radial flow regime detected). To get some estimates of the SRV parameters, one can predict a theoretical late-time radial flow regime inside the SRV using advanced methods such as the second derivative of transient pressure, to integrate the theoretical radial flow equation with the linear flow equation, and solve the two equations for the two unknowns with some uncertainties (Bahrami and Siavoshi 2013).


It should be noted that the above simplified equations do not take into account the changing gas viscosity and gas compressibility with the use of pseudo time, and they are based on single-phase flow. To get more accurate well test analysis results, gas slippage effect with pseudo time and pseudo pressure and also multi-phase parameters can be considered  in  the  analytical  approach  of  well  test  data analysis.

Uncertainties of SRV characterization using analytical  methods



The main advantage of RTA method is that in RTA, using just WHP and Q data which can always be recorded and available during production, the SRV can be characterized. However, the RTA equations may have some uncertainties as well. The disadvantage of RTA method is that the early- time  data  specially  the  linear  flow time  period may  be affected by well clean-up period, when oil and/or gas are produced with water and fracturing fluids. This causes the early-time linear flow data points on RTA plot to be scattered, and therefore mLF cannot be identified very accurately. In other words, mLF  on RTA plot may not be the actual representative of the fracture performance. On RTA plots, the data points may also be scattered because of the changes in well conditions during production (well head, choke size, liquid loading in wellbore, etc.). These affect RTA plots and calculation of slope of the lines. Also if the effective number of producing fractures (n) is not known, it makes uncertainties to be more. Running a production log is required to identify the effective number of hydraulic fractures, and improve reliability of RTA results.


The main advantage of pressure buildup test data analysis is that the flow rate is constant during the test (Q = 0); therefore, the data points are significantly less scattered compared to the data on RTA plots, and mLF can be identified more accurately using PBA plots. The accurate estimation of mLF  is very important, as it is an indicator of the well’s future production performance. Another advantage of PBA is that KSRV  is calculated independent from Ye (that maybe unknown), and therefore the calculation of Xf will be more reliable. The uncertainties associated with the PBA  method  are  when  only  a  linear  flow  regime  is observed on the diagnostic plots (one equation with two unknowns). It should be noted that the units of mLF in RTA and  PBA  equations  are  different.  Use  mLF,RTA = 4.9/ Q x mLF,PBA to convert mLF  from RTA units to mLF  from PBA units.


RTA equations are based on assuming Ye  = L/n/2 (i.e., the reservoir volume around the horizontal well between each two fracturing stages, being fully stimulated). If the actual Ye  is smaller, then RTA underestimates the Xf value. Also the RTA method is based on assuming the whole SRV being productive shale with 100 % net volume, whereas in SRV there is productive net shale volume, as well as non- productive shale volume, and in reality, NTG volume may not be 100 %.

Also SRV permeability is not uniform and it is higher near the wellbore (larger fractures), and it is lower away from wellbore (smaller fractures). Using the analytical methods to estimate SRV parameters, it assumes a homogeneous single-porosity system and a fully stimulated pore volume with 100 % productive shale (optimistic assumption), which may result in under-estimation of hydraulic fracture size.

Therefore, in designing the well spacing based on the calculated Xf, the actual well spacing should be considered larger than the one determined from RTA.

The RTA and PBA methods both have different advantages and also uncertainties and, therefore, integrating the methods can provide more accurate characterization of SRV parameters. The proposed methodology is as follows:


From RTA diagnostic plot, estimate tBDF.

From PBA Log–Log diagnostic plot, estimate mLF  and mEll.

To estimate Xf and KSRV  (2 unknowns, 2 equations) for oil wells, use the values of mLF  and mEll, and integrate Eqs. 9 and 10. For gas wells, integrate Eqs. 11 and 12. To estimate Xf, KSRV  and Ye  (3 unknowns, 3 equations) for oil wells, use the values of tBDF,  mLF  and mEll, and integrate  Eqs. 3, 9 and 10. For gas wells, integrate Eqs. 7, 11 and 12.


To examine accuracy of the proposed method for SRV characterization, a simple reservoir simulation model with one hydraulic fracture in a closed system (equivalent to one fractured stage in a multi-stage fractured well) is built and run using typical fractured shale reservoir characteristics. The input parameters into the SRV model (a single-porosity system) have been shown in Fig. 7 and Table 1. The model was run for production period of 1000 h with flow rate of 300 MSCFD, followed by 1000 h pressure buildup. Then RTA plots and equations were used to analyze the pressure drawdown data, and PBA plots and equations were used to analyze the pressure buildup data.

The RTA results for the numerically generated data are shown in Fig. 8, and PBA results in Fig. 9. The results indicate that the mLF  value calculated from both RTA and PBA methods are in good agreement, confirming reliability of the proposed equations. The calculation results for SRV parameters are reported in Table 2, which indicate that the proposed equations provide reasonably good results (average error around 10 %) for the SRV characteristics (considering the fact that the simplified analytical equations have been used to analyze the numerically generated data). In a real case, since production rate during pressure buildup test is constant (Q = 0), PBA may provide slightly more  accurate  results than  RTA;  however, due  to  economical reasons, pressure buildup tests are not common as it required the wells not to produce, but production data are always available to apply RTA method.

It should be noted that SRV is a dual-porosity and dual-permeability system, but the RTA and PBA equations are based on single-porosity system. Using the matrix porosity in RTA calculations as the effective SRV porosity, it assumes a single-porosity system and a fully stimulated pore volume, in which the volume occupied by hydraulic fractures is equal to porosity of matrix (this is an optimistic assumption, since fracture porosity is significantly less than matrix porosity). The assumption may result in under-estimation of hydraulic fracture size, and therefore RTA results may be optimistic. In this paper, reservoir simulation  and  history  matching  approach  in  a  dual-porosity model are presented to characterize SRV more accurately.

Characterizing SRV using dual-porosity simulation  model



The RTA and PBA equations are based on single-porosity models, whereas the SRV in multi-fractured shale reservoirs is actually a dual-porosity system, and the stimulated rock volume is a random network of untreated shale matrix blocks and fracture planes. Characterization of a dual- porosity SRV in multi-fractured shale reservoirs generally includes estimating the dynamic parameters such as average SRV permeability and hydraulic fracture size.


Well test analysis of multi-fractured shale reservoirs using dual-porosity concept is challenging and may not be practical due to the extensive wellbore storage effect in a long fractured wellbore, heterogeneity of the reservoir, and complexities in the dual-porosity stimulated rock volume (Restrepo and Tiab 2009), and therefore conventional well test analysis for dual-porosity behavior using interporosity flow coefficient (λ) and fracture storativity ratio (ω) may not be meaningful.


The main input parameters that are required to model fluid flow in a fractured system are fracture permeability, fracture  compressibility,  fracture  porosity,  shape  factor, and the productive shale net volume in SRV (NTGshale). Shape factor (d) can be defined for Kazemi model, as a function of matrix block size (a) that can be simplified to δ= 4/a2  (Saeedi 2012).


Fracture compressibility (Cf) in dual-porosity systems might be 1–100-folds higher than matrix compressibility. Fracture compressibility can  be  estimated  from the  following equation (Tiab et al. 2006):

where Ki is well test permeability at initial reservoir pres- sure, K is well test permeability at the current reservoir pressure (permeability from RTA or PBA), and ΔP is pressure depletion in the system.


For characterizing fracture porosity (φf), fracture opening (b) and matrix block size (or fracture spacing: a), processing of image log data if it is an open-hole well can be  used  (Dashti et  al.  2009;  Luthi  1990). However, in multi-fractured cased-hole wells, logging applications for this purpose are very limited and running fracture characterization  image  logs may  not be  practical  in  fractured cased-hole wells. The alternative method is to use history matching approach and tuning the values to estimate the fracture parameters φf and δ.


Permeability of fractures in SRV can be estimated from the following equation based on Kazemi model (Bahrami et al. 2013):

In dual-porosity reservoir simulation approach, to reduce the simulation convergence issues and improve the run time, the actual size of fracture grids (hf,actual = b) are increased to be equal to matrix grid size (hf,adjusted = L), and then permeability of the enlarged fractures are adjusted (reduced) to the new equivalent size of   the   fractures   that   results   in   the   same   fracture productivity  (Kf,actual x hf,actual = Kf,DP  x hf,adjusted).  For a  dual-porosity  simulation  model,  Eq. 14  for  fracture permeability can be written as follows:

In the above equations, Kf,actual  is the actual intrinsic permeability of fractures, Kf,DP   is the equivalent fracture permeability for dual-porosity reservoir simulation model (permeability for the fracture grids), Kwell   test   is average well test permeability (from RTA or PBA), b is average fracture aperture, a is average matrix block size, and L is average  grid size  in  the  SRV section of the  simulation model.


The downhole stresses applied from different directions to the formation rock can control aperture (and, therefore, permeability) of the hydraulically created fractures. The fractures that are perpendicular to the minimum stress direction may have significantly better permeability than the fractures perpendicular to the maximum stress direction. The permeability estimated from RTA or PBA is equivalent to SRV Kaverage,  the average permeability in the SRV model

To characterize SRV, reservoir simulation and history matching  can  be  done using dual-porosity dual-permeability model. The reservoir model should be built considering that the given dual-porosity SRV volume is surrounded by the untreated single-porosity rock. Then the cumulative volume of injected water, oil production, water production, gas production, and flowing bottom-hole  pressure  should  all  be  matched  during  history matched  using  the  matching  parameters  such  as  SRV size   and   SRV   permeability.   The   history   match   is achieved when the  correct values for matrix  and fractures  parameters in  the  SRV are  used in  the  reservoir model.

Field example:  Eagle Ford  Shale



The static and dynamic data from some multi-fractured wells in Eagle Ford shale are analyzed to characterize SRV in the reservoir. The typical petrophysical data for the shale reservoir are shown in Fig. 10: Gamma Ray versus CGR, log water saturation, neutron porosity, density porosity and corrected  porosity.  In  shale  intervals  due  to  the  bound water in the shale, porosity from Neutron Log is affected and  it  significantly overestimates the  effective  porosity, compared to density porosity. Over a shale interval with negligible TOC, the shale lithology effect on density log is not significant and density porosity should be the corrected porosity, but when TOC is significant, then density log readings are affected by presence of Kerogen. To calculate the effective porosity, TOC should be calculated first, and then density porosity log should be corrected based on that accordingly (the actual porosity in an organic-rich shale intervals is less than density porosity as shown in Fig. 10).

Typical core data have shown matrix permeability of 1E-6 md. For well test analysis, some wells with pressure and rate data available are selected to guess-estimate reservoir characteristics, and then reservoir simulation and history matching are used to characterize the reservoir as dual-porosity system, and therefore perform more accurate production forecast from the  multi-fractured  shale reservoir.

PBA for a shale gas multi-fractured horizontal well



Pressure buildup test data from a hydraulically fractured horizontal well in a shale gas field were analyzed to estimate SRV permeability and the average hydraulic fracture size. As shown in Fig. 11, the test duration is not long enough to reach radial flow regime, and analysis of the well test data have uncertainties using the conventional methods.

Using the diagnostic plot shown in Fig. 11, the linear and elliptical flow regimes are both detected, and therefore the proposed well test analysis method can be used based on  linear  flow (LF)  derivative  and  elliptical  flow (EF) derivative. By integrating the diffusivity equation solution for  the  two  reservoir  flow regimes,  the  two  unknowns K and Xf are calculated (Eqs. 11, 12): SRV permeability of 0.0005 mD and hydraulic fracture half-length of 48 ft.

RTA for shale gas and shale oil multi-fractured horizontal wells



The rate and pressure transient data are studied for two multi-fractured wells, one in a shale oil reservoir (well O1) and one in a shale gas reservoir (well G2), to show applications of RTA in characterizing SRV in multi-fractured wells.

The plots and results of RTA in the gas well G-2 are shown in Fig. 12, which shows determination of mLF  and tBDF   from the RTA plots. The values are input in Eqs. 7 and 8, and the results showed average KSRV  = 0.002 md and Xf = 94 ft.

The plots and results of RTA in the oil well O-2 are shown in Fig. 13, which shows determination of mLF  and tBDF   from the RTA plots. The values are input in Eqs. 7 and 8, and the results showed average KSRV  = 0.006 md and Xf = 70 ft.

Characterizing multi-fractured shale oil reservoir using dual-porosity simulation  model



The main objective of shale gas reservoir simulation and history matching  is to  build  a  representative simulation model that can provide a reliable production forecast. For a multi-fractured shale oil well in Eagle Ford Shale, the dual-porosity reservoir simulation model was built (to model the fracturing process in the simulation software, transmissibility of the grids that their pressure goes above rock breakdown pressure during water injection are increased to allow water to flow into those simulation cells). Oil production rate was selected as simulation control parameter, and then the well injectivity and SRV parameters were tuned to achieve a good history match for gas production rate, water production rate, well flowing bottom-hole pressure, and cumulative injected water during fracturing. The input data into the simulation model are reported in Table 3, and the history match reservoir simulation results are presented in Fig. 14. A satisfactory history match could be achieved by considering fracturing fluid loss of 54 % (history match of the other wells in this field also showed fluid loss), which could probably be loss of fluid into some non-productive  shale  zones  or  fluid loss  behind  casing through some micro-channels between cement, casing and formation (in other words, the SRV could ideally–theoretically be larger and production rates be higher, if there was no frac fluid loss). Table 4 shows the SRV parameters that provided a good history match, as a reliable model to be used in production forecasting for the multi-fractured shale oil well.




Rate transient analysis (RTA) and pressure buildup analysis (PBA) can practically be used to have estimations of the stimulated reservoir volume (SRV) parameters in shale reservoirs.


In RTA method, the slope of the linear flow regime data (mLF) and the time at which boundary-dominated flow starts (tBDF)  are used for KSRV  and Xf estimation.


In  PBA  method, integration  of  diffusivity equations solutions  for  linear  and  elliptical  flow regimes  can provide estimation of SRV characteristics. The values of mEll  and mLF  can be determined, respectively, from zero slope line on elliptical flow derivative and zero slope line on linear flow derivative.


The  value  of  mLF can be a good  indicator  of  well deliverability. The mLF from PBA may be more reliable than  the  value  of  mLF   from  RTA,  since  during  a pressure  buildup  test,  production  rate  is  constant (Q = 0)  and,  therefore,  the  data  quality  is  better. However, in  terms of  data  availability  in  producing wells, most wells have production data for RTA, but few wells may have the pressure buildup data required for PBA.


Integration  of  rate  transient  analysis  and  pressure buildup analysis methods can reduce uncertainties in SRV characterization.


Reservoir simulation and history matching using dual-porosity SRV model is more reliable method for SRV characterization, compared  with the analytical methods.






The authors would like to thank Aurora Oil and Gas management for supporting the technical research study, Chet Ozgen and Basar Basbug from Nitec LLC for the technical discussions and support on simulation of multi-fractured shale reservoirs, and Baharak Ghaffari Nik (Schlumberger) for providing useful technical information on this topic. We acknowledge KAPPA-Engineering for use of KAPPA-ECRIN, Coat’s Engineering for use of SENSOR reservoir simulator,  and also  LYNX pre- and post-processing reservoir simulation tool in this research study.


Open  Access   This  article  is  distributed under  the  terms  of  the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.


Author: Nick Bahrami, David Pena, Ian Lusted, Aurora Oil and Gas, Houston, USA


Appendix A: General  forms of fluid flow equations for different  flow regimes


(Bahrami et al. 2013; Martinez et al. 2012; Kappa Engineering 2012)

In the above equations, P is pressure (psia), t is time (h), q is production rate (oil: STBD, gas: MSCFD), h is reservoir thickness (ft), B is formation volume factor, µ is viscosity (cp),  Ct   is  total  compressibility, u is  effective  average porosity (fraction), φ is permeability of the stimulated rock volume, T is reservoir temperature (R), mRF   is pressure derivative value for infinite acting radial flow, mLF  is pressure derivative value for infinite acting linear flow, and mEF is pressure derivative value for infinite acting elliptical flow. For multi-fractured horizontal wells, n is total number of fractures and h is total reservoir thickness, but for multi- fractured vertical wells, n = l for h as total reservoir thickness of all the fractured layers.Time functions for pressure derivative calculation in a pressure buildup test:




Arevalo JA, Wattenbarger R, Samaniego F (2002) Analysis of performance in tight gas wells. In: SPE annual technical conference, New Orleans, USA


Badazhkov D (2008) Analysis of production data with elliptical flow regime in tight gas reservoirs. In: SPE 117023, SPE Russian oil and gas technical conference and exhibition, Moscow, Russia


Bagherian B, Sarmadivaleh M, Ghalambor A, Nabipour A, Rasouli V, Mahmoudi M (2010) Optimization of multiple-fractured hori- zontal tight gas well. In: SPE 127899, symposium and exhibition on formation damage control. Society of Petroleum Engineers, Lafayette, Louisiana


Bahrami H, Rezaee R, Asadi M, Minaeian V (2013) Effect of stress anisotropy on well productivity in unconventional gas reservoir. In: SPE 164521, SPE production and operations symposium. Oklahoma, USA


Bahrami H, Siavoshi J (2013) Interpretation of reservoir flow regimes and analysis of welltest data in unconventional reservoirs. In: SPE-164033, Middle East unconventional gas conference and exhibition held in Muscat, Oman Bourdarot G  (1998) Well  testing  interpretation  methods. Editions Technip


Dashti R, Bagheri MB (2009) Fracture characterization of a porous fractured carbonate reservoir. In: SPE 125329, SPE reservoir characterization conference, UAE Fan L, Thompson JW, Robinson JR (2010)


Understanding gas production mechanism and effectiveness of well stimulation in the  Haynesville  shale  through  reservoir  simulation.  In:  SPE 136696, Canadian  unconventional  resources and international petroleum conference, Calgary, Alberta


Hossain  MMM,  Rahman  MK  (2008)  Numerical  simulation  of complex fracture growth during tight reservoir stimulation by hydraulic fracturing. J Pet Sci Eng 60:86–104


Ibrahim  M, Canal  SU, Wattenbarger  RA (2006) Analysis of rate dependence in transient linear flow in tight gas wells. In: International petroleum exhibition and conference, Abu Dhabi, UAE


KAPPA engineering (2012) Dynamic Data Analysis

Kappa Engineering (2015) Kappa DDA book unconventional resources chapter


Luthi SM (1990) Fracture apertures from electrical borehole scans. J Geophys 55(7):821–833


Malallah A, Nashawi I, Algharaib M (2007) Constant pressure analysis   of   oil   wells   intercepted   by   infinite  conductivity hydraulic  fracture  using  rate  derivative  functions.  In:  SPE Middle East oil and gas conference, Bahrain


Martinez JA, Escobar FH, Bonilla LF (2012) Reformulation of the elliptical flow governing equation for a more complete welltest data interpretation in horizontal wells. ARPN J Eng Appl Sci 7(3)


Rasdi F, Chu L (2012) Diagnosing fracture network pattern and flow regimes aids production performance analysis in unconventional oil reservoirs. In: SPE European unconventional resources conference, Vienna, Austria


Restrepo DP, Tiab D (2009) Multiple fractures transient response. In: SPE 121594, Latin American and Caribbean petroleum engi- neering conference. Cartagena, Colombia


Saeedi A (2012) Experimental study of multi-phase flow in porous media during CO2 geo-sequestration processes, Springer Theses. Springer, Berlin


Tiab  D,  Restrepo  DP,  Lgbokoyi  A  (2006)  Fracture  porosity  of naturally fractured reservoirs. In: SPE 104056, First international oil conference and exhibition in Mexico, Mexico