Gas Flowing Material Balance
This methodology should only be used for dry/wet gas wells.
1. Theory
This module is designed based on the work presented by Palacio and Blasingame (1993) known as "Conventional Gas Flowing Material Balance". The work involves plotting rate-normalized pseudopressure \(\left(\frac{\Delta m_{wf}}{q_g}\right)\) vs. material balance pseudotime (\(t_{ma}\)). During pseudosteady state flow, this results in a straight line with slope inversely proportional to the contacted pore volume \(\left(V_p\right)\).
Conventional gas FMB is given by
Where \(\Delta m_{wf}\) is real gas pseudo-pressure drop defined as
in which subscript \(i\) and \(wf\) reflect initial and well flowing bottomhole condition, respectively. Gas viscosity (cp) is represented by \(\mu_g\), while Z is gas compressibility factor and \(p\) is pressure (psia).
Palacio and Blasingame (1993) introduced material balance pseudotime, \(t_{ma}\) (days) to account for the changes in gas viscosity and compressibility with the decline of the average pressure.
Where \(\bar{p}\) denotes the average reservoir pressure in psia, \(q_g\) is gas production rate in Mcf/d and \(c_t\) is total compressibility in 1/psia.
Total compressibility accounts for the compressibility of the fluid phases present in the system as well as the rock compressibility and is written as
where \(c_g\), \(c_o\), \(c_w\) and \(c_r\) denote gas, oil, water and formation (rock) compressibility, respectively. The default value of \(c_r\) is 4x10\(^{-6}\) 1/psia and \(c_w\) is 2.8x10\(^{-6}\) 1/psia in whitson+.
According to Eq. \eqref{eq:gas-FMB}, the slope of the straight line during pseudo-steady-state flow is denoted by \(m_{pss}\) in psia\(^2\)/cp-Mscf,
Original gas in place is calculated by,
with gas formation volume factor given by,
Rearranging Eq. \eqref{eq:slope-of-straight-line} with Eqs. \eqref{eq:OGIP} and \eqref{eq:Bg}, then contacted pore volume \(\left(V_p\right)\) in Mscf can be determined.
Subscript \(sc\) refers to standard condition and \(T\) is temperature in Rankine.
From simple material balance (ignoring rock and water compressibility),
Please note that this method requires the user to input first guess of \(G_i\) to determine \(t_{ma}\). It uses an iterative procedure to correctly estimate the \(G_i\).
2. Workflow
Gas FMB Workflow
- From provided PVT data, calculate initial properties (\(B_{gi}\), \(c_{gi}\), \(\mu_{gi}\) and \(m(p_{i})\)).
- From the production data, calculate the following properties for every time step:
- Cumulative gas production (\(G_p(t)\))
- From initial guess of gas in place (\(G_{i-guess}\)), calculate average gas formation volume factor (\(\bar{B_g}(t)\)) using Eq. \eqref{eq:Average-Bg}
- From PVT data, find average pressure (\(\bar{p}(t)\)) that gives \(\bar{B_g}(t)\)
- Do the same interpolation for average gas viscosity and compressibility (\(\bar{\mu_g}(t)\) and \(\bar{c}(t)\)) at \(\bar{p}(t)\)
- Material balance pseudo time using Eq. \eqref{eq:MB-Pseudotime}
- From PVT data, calculate \(m(p_{wf})\)
- Rate-normalized pseudopressure \(\left(\frac{\Delta m_{wf}(t)}{q_g(t)}=\frac{m(p_i)-m(p_{wf})}{q_g(t)}\right)\)
- Plot material balance pseudotime \(\left(t_{ma}\right)\) in x-axis and rate-normalized pseudopressure \(\left(\frac{\Delta m_{wf}(t)}{q_g(t)}\right)\) in y-axis.
- Find slope (\(m_{pss}\)) and intercept (\(b_{pss}\)).
- Calculate contacted pore volume (\(V_p\)) using Eq. \eqref{eq:pore-volume}.
- Calculate original gas in place (\(G_i\)) using Eq. \eqref{eq:OGIP}.
3. Gas Flowing Material Balance (FMB) accounting for Rock Compressibility
Consider a tank with initial volume \(V_{pi}\) containing gas at initial pressure \(p_{i}\) and water at initial saturation \(S_{wi}\). The compressibility of the tank (rock) and water are \(c_r (p)\) and \(c_w (p)\). Water is assumed to be immobile in the system (if it is mobile we can use the multiphase FMB instead). If the initial volume of gas in place is \(G_{i}\), then we can express this in terms of the tank volume as
After removal of a certain volume of gas (\(G_p\)) the tank pressure falls to p, the insitu volume of gas is given by
In which (\(Υ_V\)) and (\(Υ_w\)) are the volume multiplier for rock and water irrespectively. We assume that (\(Υ_w (p)\)) is 1 in this module and that water compressibility, \(c_w\), is 2.8x10\(^{-6}\) 1/psia.
Substituting for \(V_{pi}\) from Eq. \eqref{eq:GiBgdi}, we have
or
4. Computing Original Gas in Place with Adsorption
Ambrose et al.[2] argue that the free OGIP must be corrected for the adsorbed gas that occupies some of the hydrocarbon pore volume (HCPV). If porosity is estimated from core plugs where core preparation has removed the adsorbed gas, then the OGIP will be overestimated.
The convention of reporting OGIP when adsorption is included, is to report it in units of scf/ton (standard cubic feet per short ton). For a fluid system having only dry/wet gas and water, the total OGIP is
where
and with \(C_1\approx 5.7060\) and \(C_2\approx1.318\cdot10^{-6}\) are unit conversion factors.
The OGIP resolved in GFMB is the free OGIP. We can use this volume to compute the rock mass, which in turn can be used to calculate the adsorbed OGIP. Multiplying equation \eqref{eq:freegasperton} with the rock mass (\(G=\hat{G}_f m_r\)) gives the free OGIP in units of scf. Solving for the rock mass yields
This rock mass can then be multiplied by equation \eqref{eq:adsorbedgasperton} to get the adsorbed OGIP in scf, i.e.,
For more information regarding adsorption, please refer to this page
Total Compressibility when Adsorption is ON
In order to honor mass balance, it becomes necessary to modify total system compressibility when accounting for adsorption. Total system compressibility can be determined using the following equation:
The derivation of this equation is shown here and is contructed upon Ambrose et al.[2].