| THEORY & FORMULAE |
The M.B.E. remains an important analytical tool in reservoir engineering, employed mainly in the preliminary design of reservoir development and operations, prelude/complement to fully-fledged reservoir simulation, and in economic/planning studies.
The modified Schilthuis procedure is implemented here to compute oil and gas production and ultimate recovery from a combination of solution gas drive/gas-cap drive reservoir, producing at or below the bubble-point pressure. Direct production/injection of gas fro/to the gas cap is explicitly included here. It is assumed that the gas-oil contact does not move into the oil zone, but gas expansion from the gas cap will diffuse through the oil zone. Also water encroachment is considered negligible.
User Input: The user is expected to enter PVT and Relative-Permeability data, each consisting of 10 to 15 data points. The PVT pressure values should be in descending order and pressure intervals must be constant. The input may include the fraction of produced gas (solution or diffused gas-cap gas) withdrawn from/injected into the gas cap. Injection into the gas cap is indicated by a negative sign, and production is assumed positive. The liquid saturation values in the Rel. Perm table, should be in ascending order rising to 1.0.
The main MBE equation and three auxillary equations are given below:
    
where:
     SDI = solution gas drive index
     GDI = gas-cap gas drive index
     N = initial oil in place, bbl
     NPn = cumulative oil produced at pressure step n, fraction of N
     GPn = cumulative gas produced at pressure step n, fraction of N
     Gin = cumulative gas injected to pressure step n, fraction of N
     Mi = initial ratio of gas-cap to oil reservoir pore space
     Bt = two-phase formation vol. factor at pressure step n
     Bo = oil formation vol. factor at pressure step n
     Boi = initial oil formation volume factor, bbl/scf
     Bgn = gas formation vol. factor at pressure step n
     Bgi = initial gas formation volume factor, res.bbl/scf
     Rsn = solution gas-oil ratio at pressure step n
     Rsi = initial solution gas-oil ratio, scf/bbl
     RIn = instataneous gas-oil ratio
     SLn = total liquid saturation at pressure step n
     Sw = water saturation
     μ = viscosity, cp
     k = permeability, mD
The cubic spline interpolation method is used to generate continuous values of relative perm. data from the entered set of discrete values.