PO.VFP.MukherjeeBrill.Pin.SI
Description
Calculates inlet pipe pressure using Mukherjee-Brill (1985), [kPa]. For inclined wells. Gas properties (ρg, μg, Bg) use DAK (Z-factor) and LGE (viscosity) (SI units).
Formula
=LAMBDA(Ql, Qg, P_out, T, Rho_l, Ul, SG_gas, IFT_gl, pipe_ID, pipe_length, pipe_roughness, pipe_angle, LET(
Ql_bbld, PO.UnitConverter(Ql, "m3/d", "bbl/d"),
Qg_mmscfd, PO.UnitConverter(Qg, "sm3/d", "mmscf/d"),
P_out_psi, PO.UnitConverter(P_out, "kPa", "psi"),
T_degF, PO.UnitConverter(T, "degC", "degF"),
Rho_l_lbft3, PO.UnitConverter(Rho_l, "kg/m3", "lb/ft3"),
IFT_gl_dynecm, PO.UnitConverter(IFT_gl, "mN/m", "dyne/cm"),
pipe_ID_in, PO.UnitConverter(pipe_ID, "cm", "in"),
pipe_length_ft, PO.UnitConverter(pipe_length, "m", "ft"),
result_psi, PO.VFP.MukherjeeBrill.Pin(Ql_bbld, Qg_mmscfd, P_out_psi, T_degF, Rho_l_lbft3, Ul, SG_gas, IFT_gl_dynecm, pipe_ID_in, pipe_length_ft, pipe_roughness, pipe_angle),
PO.UnitConverter(result_psi, "psi", "kPa")))Parameters
| Parameter | Unit | Description |
|---|---|---|
Ql |
m3/d | Liquid rate, [m3/d] |
Qg |
sm3/d | Gas rate, [sm3/d] |
P_out |
kPa | Outlet pipe pressure, [kPa] |
T |
degC | Temperature at middle point of the pipe, [degC] |
Rho_l |
kg/m3 | Liquid density, [kg/m3] |
Ul |
- | Liquid viscosity, [cP] |
SG_gas |
- | Gas specific gravity, [dimensionless] |
IFT_gl |
mN/m | Liquid-gas interfacial (surface) tension, [mN/m] |
pipe_ID |
cm | Pipe inner diameter, [cm] |
pipe_length |
m | Pipe length, [m] |
pipe_roughness |
- | Pipe relative roughness, [dimensionless] |
pipe_angle |
- | Pipe angle from horizontal, [degrees]. Horizontal flow: angle = 0; Vertical well producer: angle = 90; Vertical well injector: angle = -90 |
Examples
| Invocation | Expected |
|---|---|
=PO.VFP.MukherjeeBrill.Pin.SI(159, 28317, 6895, 93.33, 800, 2, 0.65, 20, 6.35, 914.4, 0.0006, 90) |
10755.57 |