There is a huge demand in the industry to forecast production in shale gas reservoirs accurately. There are many methods including volumetric, Decline Curve Analysis (DCA), analytical simulation and numerical simulation. Each one of these methods has its advantages and disadvantages, but only the DCA technique can use readily available production data to forecast rapidly and to an extent accurately. The DCA methods in use in the industry such as the Arps method had originally been developed for Boundary dominated flow (BDF) wells but it has been observed in shale reservoirs the predominant flow regime is transient flow. Therefore it was imperative to develop newer models to match and forecast transient flow regimes. The SEDM/SEPD, the Duong model and the Arps with a minimum decline rate are models that have the ability to match and forecast wells with transient flow followed by boundary flow. I have revised the Duong model to forecast better than the original model. I have also observed a certain variation of the Duong model proves to be a robust model for most of the well cases and flow regimes. The modified Duong has been shown to work best compared to other deterministic models in most cases. For grouped datasets the SPED & Duong models forecast accurately while the Modified Arps does a poor job.

Production forecasting in shale (ultra-low permeability) gas reservoirs is of great interest due to the advent of multi-stage fracturing and horizontal drilling. The well renowned production forecasting model, Arps? Hyperbolic Decline Model, is widely used in industry to forecast shale gas wells. Left unconstrained, the model often overestimates reserves by a great deal. A minimum decline rate is imposed to prevent overestimation of reserves but with less than ten years of production history available to analyze, an accurate minimum decline rate is currently unknown; an educated guess of 5% minimum decline is often imposed. Other decline curve models have been proposed with the theoretical advantage of being able to match linear flow followed by a transition to boundary dominated flow. This thesis investigates the applicability of the Stretched Exponential Production Decline Model (SEPD) and compares it to the industry standard, Arps' with a minimum decline rate. When possible, we investigate an SEPD type curve. Simulated data is analyzed to show advantages of the SEPD model and provide a comparison to Arps' model with an imposed minimum decline rate of 5% where the full production history is known. Long-term production behavior is provided by an analytical solution for a homogenous reservoir with homogenous hydraulic fractures. Various simulations from short-term linear flow (~1 year) to long-term linear flow (~20 years) show the ability of the models to handle onset of boundary dominated flow at various times during production history. SEPD provides more accurate reserves estimates when linear flow ends at 5 years or earlier. Both models provide sufficient reserves estimates for longer-term linear flow scenarios. Barnett Shale production data demonstrates the ability of the models to forecast field data. Denton and Tarrant County wells are analyzed as groups and individually. SEPD type curves generated with 2004 well groups provide forecasts for wells drilled in subsequent years. This study suggests a type curve is most useful when 24 months or less is available to forecast. The SEPD model generally provides more conservative forecasts and EUR estimates than Arps' model with a minimum decline rate of 5%.

This study compares four recently developed decline curve methods and the traditional Arps or Fetkovich approach. The four methods which are empirically formulated for shale and tight gas wells are: 1. Power Law Exponential Decline (PLE). 2. Stretched Exponential Decline (SEPD). 3. Duong Method. 4. Logistic Growth Model (LGM). Each method has different tuning parameters and equation forms. The main objective of this work is to determine the best method(s) in terms of Estimated Ultimate Recovery (EUR) accuracy, goodness of fit, and ease of matching. In addition, these methods are compared against each other at different production times in order to understand the effect of production time on forecasts. As a part of validation process, all methods are benchmarked against simulation. This study compares the decline methods to four simulation cases which represent the common shale declines observed in the field. Shale wells, which are completed with horizontal wells and multiple traverse highly-conductive hydraulic fractures, exhibit long transient linear flow. Based on certain models, linear flow is preceded by bilinear flow if natural fractures are present. In addition to this, linear flow is succeeded by Boundary Dominated Flow (BDF) decline when pressure wave reaches boundary. This means four declines are possible, hence four simulation cases are required for comparison. To facilitate automatic data fitting, a non-linear regression program was developed using excel VBA. The program optimizes the Least-Square (LS) objective function to find the best fit. The used optimization algorithm is the Levenberg-Marquardt Algorithm (LMA) and it is used because of its robustness and ease of use. This work shows that all methods forecast different EURs and some fit certain simulation cases better than others. In addition, no method can forecast EUR accurately without reaching BDF. Using this work, engineers can choose the best method to forecast EUR after identifying the simulation case that is most analogous to their field wells. The VBA program and the matching procedure presented here can help engineers automate these methods into their forecasting sheets.