Linear Material Balance Equation (MBE) in Reservoir Engineering and Excel Application
Hi 🙋♂️ I am Franklin Obika, a reservoir engineer and data scientist passionate about solving oil and gas challenges using Excel, data analytics, and machine learning. Beyond my technical skills, I am recognized as a valuable team member, always eager to collaborate, share knowledge, and contribute to a dynamic work environment. I am committed to continuous learning and always seek new challenges where I can apply my skills to make a meaningful difference. Let's connect if you are exploring opportunities in AI, machine learning, or Excel, or if you are interested in collaborating on impactful research projects! 😊🤝
1. Introduction
The material balance equation (MBE) is one of the most fundamental tools in reservoir engineering. It is based on the principle of conservation of mass and states that the total volume of fluids produced from a reservoir is equal to the total expansion of fluids and rock within the reservoir plus any fluid entering the reservoir from external sources (such as an aquifer or injection). Through MBE, engineers can estimate original hydrocarbons in place, identify dominant drive mechanisms, and evaluate reservoir performance using production and pressure data. Despite its usefulness, the general form of the material balance equation is mathematically complex. It contains multiple expansion terms for oil, gas, water, rock, and sometimes injected fluids. This complexity makes direct interpretation difficult, especially when trying to diagnose reservoir drive mechanisms or validate reservoir models.
To overcome this challenge, Havlena and Odeh reformulated the material balance equation into a linear (straight-line) form. This transformation allows engineers to plot carefully selected variable groups against one another. When the correct reservoir model is assumed, the data align along a straight line. The slope and intercept of this line then provide key reservoir parameters such as stock-tank oil initially in place (STOIIP), gas-cap size, or aquifer strength. Deviations from linearity indicate the presence of unaccounted energy mechanisms or data inconsistencies.
2. Linear Form of the Material Balance Equation
The linear form of the MBE expresses total underground withdrawal as a linear combination of fluid expansion and external energy terms. In its general form, it can be written as:
$$F = N E_o + m E_g + E_{f,w} + W_e B_w$$
where:
- F is the total underground withdrawal,
$$F = N_p[B_o+(Rp-Rs)B_g]+WpBw$$
- Eo represents oil and dissolved gas expansion,
$$E_o=(B_o-B_{oi})-(R_{si}-Rs)B_g$$
- Eg represents gas-cap expansion,
$$E_g=B_{oi}\left(\frac{B_g}{B_{gi}}-1\right)$$
- \(E_{f,w}\) accounts for rock and connate water expansion,
$$(1+m)B_{oi}\left(\frac{c_wS_{wi}+c_f}{1-S{wi}}\right)\Delta P$$
WeBw represents water influx,
N is the stock-tank oil initially in place,
m is the gas-cap size ratio.
Depending on reservoir conditions, some terms become negligible, allowing the equation to reduce to simpler linear forms. The following sections describe the most common reservoir scenarios.
2.1 Undersaturated Reservoir Without Water Influx
An undersaturated reservoir is one in which the reservoir pressure remains above the bubble-point pressure, so no free gas is produced. If the reservoir is also volumetric, meaning there is no water influx, the only sources of energy are oil expansion, rock compaction, and connate water expansion.
Under these assumptions, the material balance equation reduces to:
$$F = N (E_o + E_{f,w})$$
In this case:
\(F = N_p B_o\),
\(E_o = B_o - B_{oi}\),
\(E_{f,w}=(1+m)B_{oi}\left(\frac{c_wS_{wi}+c_f}{1-S{wi}}\right)\Delta P\)
When F is plotted against \((E_o + E_{f,w})\), the data form a straight line passing through the origin. The slope of this line gives the value of N, the stock-tank oil initially in place. A straight line confirms purely volumetric depletion, while curvature suggests additional energy support or data errors.
2.2 Undersaturated Reservoir with Water Influx
In many reservoirs, pressure support is provided by an aquifer even before the bubble point is reached. In this case, the reservoir is still undersaturated, but water influx contributes additional energy.
The linearized material balance equation becomes:
$$F = N (E_o + E_{f,w}) + W_e$$
Dividing through by \((E_o + E_{f,w})\) gives:
$$\frac{F}{E_o + E_{f,w}} = N + \frac{W_e}{E_o + E_{f,w}}$$
where
$$F=N_pB_o+W_pB_w$$
$$E_o=B_o-B_{oi}$$
$$E_{f,w}=(1+m)B_{oi}\left(\frac{c_wS_{wi}+c_f}{1-S{wi}}\right)\Delta P$$
A plot of \((\frac{F}{E_o + E_{f,w}})\) versus \((\frac{W_e}{E_o + E_{f,w}})\) yields a straight line with:
an intercept equal to N,
a slope of unity.
2.3 Saturated Reservoir with Water Influx
A saturated (or a two-phase reservoir) reservoir operates below the bubble-point pressure, meaning free gas evolves from solution and contributes to production. If water influx is also present, the reservoir benefits from both internal gas expansion and external aquifer support.
For this case, the linear MBE reduces to:
$$F = N E_o + W_e$$
where:
\(F = N_p B_o + (R_p - R_s) B_g + W_p B_w\),
\( E_o=(B_o-B_{oi})-(R_{si}-Rs)B_g\)
Dividing through by Eo:
$$\frac{F}{E_o} = N + \frac{W_e}{E_o}$$
A straight-line plot of \((\frac{F}{E_o})\) versus \((\frac{W_e}{E_o})\) gives:
intercept = N,
slope = 1.
The position and trend of the data indicate whether the aquifer is weak, moderate, or strong. Excessive curvature may suggest an improperly modeled aquifer or inaccurate water influx estimates.
2.4 Saturated Reservoir Without Water Influx
This scenario represents a solution-gas drive reservoir with no external pressure support. All reservoir energy comes from oil expansion and gas liberation.
The material balance equation simplifies to:
$$F = NE_o$$
where:
\(F = N_p [B_o + (R_p - R_s) B_g]\),
\(E_o = (B_o - B_{oi}) + (R_{si} - R_s) B_g\).
A plot of F versus Eo produces a straight line passing through the origin. The slope of this line is the STOIIP, N.
2.5 Gas Cap Drive Reservoir
A gas cap drive reservoir is characterized by the presence of a free gas zone overlying the oil column at initial conditions. As reservoir pressure declines, gas in the cap expands and provides additional energy that supports oil production. This mechanism often results in slower pressure decline compared to purely solution-gas drive reservoirs.
In a gas-cap drive reservoir, the effects of natural water influx and the compressibility of water and rock pores are typically considered negligible due to the overwhelming presence of highly compressible gas. The linear form of the material balance equation is written as:
$$F = N(E_o + mE_g)$$
where:
\(F = N_p [B_o + (R_p - R_s) B_g]\)
\(E_o = (B_o - B_{oi}) + (R_{si} - R_s) B_g\)
\(E_g=B_{oi}\left(\frac{B_g}{B_{gi}}-1\right)\)
Known Gas-Cap Size
When the gas-cap size (m) is known from geological or volumetric studies, (Eo + mEg) is evaluated at each pressure point and plotted against F. A straight line passing through the origin confirms the dominance of gas-cap drive. The slope of the line gives the value of N.
Unknown Gas-Cap Size
When the gas-cap size is unknown, the equation is rearranged to a linear form:
$$\frac{F}{E_o} = N + mN\left(\frac{E_g}{E_o}\right)$$
A plot of \((\frac{F}{E_o})\) versus \((\frac{E_g}{E_o})\) yields:
an intercept equal to the STOIIP, N,
a slope equal to mN.
From this plot, the gas-cap size (m) is obtained by dividing the slope by the intercept. Once m and N are known, the gas initially in place (GIIP) can be calculated.
2.6 Combination Drive Reservoir
A combination drive reservoir is one in which multiple drive mechanisms contribute simultaneously to production. These mechanisms may include solution-gas drive, gas-cap expansion, rock and connate water expansion, and water influx from an aquifer. Most real reservoirs fall into this category rather than exhibiting a single, idealized drive mechanism.
For a combination drive system, the linear material balance equation takes the form:
$$F = N E_t + W_e$$
where:
\(F = N_p [B_o + (R_p - R_s) B_g]+W_pB_w\)
\(E_t = E_o + mE_g + E_{f,w}\) is the total expansion term
\(E_o = (B_o - B_{oi}) + (R_{si} - R_s) B_g\)
\(E_g=B_{oi}\left(\frac{B_g}{B_{gi}}-1\right)\)
\( E_{f,w}=(1+m)B_{oi}\left(\frac{c_wS_{wi}+c_f}{1-S{wi}}\right)\Delta P\)
Dividing through by Et gives:
$$\frac{F}{E_t} = N + \frac{W_e}{E_t}$$
A plot of \((\frac{F}{E_t})\) versus \((\frac{W_e}{E_t})\) produces a straight line with:
an intercept equal to the STOIIP, N,
a slope of unity.
2.7 Solving Examples 6.3 – 6.5
This section demonstrates how the linear form of the material balance equation is applied in practice using representative reservoir scenarios. The examples (examples 6.3, 6.4, and 6.5) from Sylvester Okotie’s textbook, “Reservoir Engineering Fundamentals and Applications,” illustrate how straight-line plots are constructed and interpreted to estimate key reservoir parameters such as stock-tank oil initially in place (STOIIP), gas-cap size, and gas initially in place (GIIP).
Example 6.3: Saturated Reservoir Without Water Influx
This example considers a saturated oil reservoir operating below the bubble-point pressure with no aquifer support. Under these conditions, reservoir energy comes solely from oil expansion and solution gas liberation.
The appropriate linear material balance form, as discussed in section 2.4 of this article, is:
$$F = N E_o$$
The calculation procedure involves:
Computing F at each pressure point using cumulative oil production and produced gas in excess of solution gas.
Evaluating Eo from PVT data by combining oil expansion and gas liberation effects.
Plotting F versus Eo.
The solution to this example in Excel is shown in Figure 1.
Figure 1: Solution to Example 6.3 in Excel
Because there is no water influx, the plot should pass through the origin. The slope of the straight line obtained from the data represents the STOIIP, gotten as 776.56 MMstb. In Example 6.3, the data align well along a straight line with an R-squared score of 0.993 (Figure 1), confirming a solution-gas drive mechanism and validating the absence of aquifer support.
Example 6.4: Gas-Cap Drive Reservoir with Known Gas-Cap Size
In this case, the reservoir contains a gas cap of known size and is supported by both gas-cap expansion and limited water influx. The objective is to estimate the STOIIP using the linear material balance equation.
For a gas-cap drive reservoir with known gas-cap size (m), the linear form, as discussed in section 2.5 of this article, becomes:
$$F = N (E_o + mE_g)$$
The solution steps are:
Calculate F using cumulative oil production and produced gas.
Determine Eo from PVT data by combining oil expansion and gas liberation effects.
Compute Eg
Combine the expansion terms as (Eo + mEg).
Plot F versus (Eo + mEg).
Figure 2: Solution to Example 6.4 in Excel
A straight-line relationship, as shown in Figure 1, confirms that the assumed gas-cap size is reasonable. The slope of the line gives the STOIIP. In Example 6.4, the data form a linear trend with an R-squared score of 0.9998, and the resulting slope yields the estimated stock-tank oil initially in place for the reservoir as 3440 MMstb.
Example 6.5: Determination of Gas-Cap Size, STOIIP, and GIIP
Example 6.5 extends the previous analysis by treating the gas-cap size as unknown. Here, the goal is to determine:
the gas-cap size (m),
the STOIIP (N),
and the gas initially in place (G).
The linearized equation is rearranged as:
$$\frac{F}{E_o} = N + mN \left( \frac{E_g}{E_o} \right)$$
This form is analogous to a straight-line equation where:
the intercept equals (N),
the slope equals (mN).
The procedure involves:
Computing F, Eo, and Eg at each pressure point.
Computing the ratios \((\frac{F}{E_o})\) and \((\frac{E_g}{E_o})\).
Plotting \((\frac{F}{E_o})\) versus \((\frac{E_g}{E_o})\).
Figure 3: Solution to Example 6.5 in Excel
From the resulting straight line in Figure 3:
the intercept directly gives the STOIIP as 37.377 MMstb,
the slope divided by the intercept yields the gas-cap size (m) as 1.679.
Once m and N are known, the gas initially in place is calculated using the relationship between gas-cap size, oil volume, and formation volume factors, which yields 90487.309 MMscf. In Example 6.5, this method provides consistent estimates of all three parameters, demonstrating the robustness of the linear MBE approach for gas-cap reservoirs.
Conclusion
The linear form of the material balance equation provides a powerful and practical framework for understanding reservoir performance and estimating key reservoir parameters. By transforming the general material balance equation into straight-line relationships, complex interactions between fluid expansion, production, and external energy support become easier to interpret and diagnose. Through the various reservoir scenarios discussed, ranging from undersaturated and saturated systems to cases with and without water influx or gas-cap drive, the linear MBE demonstrates its flexibility and robustness. The worked examples further show how simple plots can be used to estimate stock-tank oil initially in place, gas-cap size, and gas initially in place, while simultaneously validating reservoir drive mechanisms. When production and pressure data align linearly, confidence in the assumed reservoir model is strengthened; when they do not, the deviation itself provides valuable diagnostic insight. Ultimately, the strength of the linear material balance approach lies not only in its ability to quantify hydrocarbons in place but also in its role as a reservoir diagnostic and history-matching tool. When applied with reliable PVT and production data, it remains an essential technique for reservoir engineers, bridging the gap between theoretical reservoir behavior and real-field performance.
