Skip to content
Brent$117.6Gas TTF$1,192USD453.87 ₸
MunaiHub

Equations: Distributed pressure measurements

The equations of the “Distributed pressure measurements” lessons in “Well testing and pressure transient analysis”, with their conditions and units, and what each symbol means in them.

Drawdown mobility at steady flow to the probe
(kμ)d=C q2π rpe Δpd\left(\frac{k}{\mu}\right)_d = \frac{C\,q}{2\pi\,r_{pe}\,\Delta p_d}

SI (any consistent system)

Holds when: The flow is steady: the piston works for longer than the pressure needs to level off at a constant value (advanced section). One phase (one fluid) of known viscosity flows, usually the filtrate in the invaded zone. The rock at the probe is uniform, and the flow follows Darcy's law. In an anisotropic rock (one with different permeabilities in different directions), the result reflects a mix of horizontal and vertical permeability. It usually lies between the spherical and the horizontal permeability. There is no extra resistance at the opening of the probe from leftover mudcake, damaged rock or fine particles. If there is, the mobility comes out too low.

Pressure measurement with the wireline formation tester →
Random error of a point on the pressure–depth plot
σp2=σgauge2+(G σz)2\sigma_p^2 = \sigma_{\text{gauge}}^2 + (G\,\sigma_z)^2

SI (in practice kPa, kPa/m and m)

Holds when: The errors of pressure and depth are random and independent. Systematic errors (errors that always push the same way, unlike random ones) are not included: a common shift of the depths, or a constant error of the gauge. They shift the line but do not change its slope. An error that grows with depth (cable stretch, survey error) changes the slope too.

Pressure measurement with the wireline formation tester →
Pressure gradient from the points by the least-squares method
G=∑i(zi−zˉ) pi∑i(zi−zˉ)2,sG=s∑i(zi−zˉ)2G = \frac{\sum_i (z_i - \bar z)\,p_i}{\sum_i (z_i - \bar z)^2},\qquad s_G = \frac{s}{\sqrt{\sum_i (z_i - \bar z)^2}}

any consistent system (in practice MPa and m; the gradient in kPa/m)

Holds when: All the points belong to one continuous phase (for example, one column of oil) in hydrostatic equilibrium: the fluid is at rest, so the pressure changes only with depth. They also belong to one hydraulically connected system, in which fluid can flow between all the parts. The errors of the points are random and of equal size. Points that we know are shifted in one direction, such as a supercharged point or an unfinished build-up, have been excluded beforehand. The density of the fluid in the interval is constant. For tall gas columns this is not so (pressure in a gas column).

Pressure measurement with the wireline formation tester →
Pressure at the probe with tool storage
S  dp dt=−qpist+k/μK (pres−p),τ=S Kk/μS\,\frac{\dd p}{\dd t} = -q_{\text{pist}} + \frac{k/\mu}{K}\,(p_{\text{res}} - p),\qquad \tau = \frac{S\,K}{k/\mu}

cm³, MPa, s; mobility in mD/(mPa·s)

Holds when: The flow to the probe is quasi-steady: in the rock at the probe the pressure adjusts much faster than in the tool. Tool storage is constant. The transient spherical flow in the formation (flow that has not yet become steady) is not included. It gives a slow “tail” to the build-up: the pressure gains the last few kilopascals not exponentially but noticeably more slowly. So the stabilisation time from this estimate is a lower limit. With a strict criterion (tenths of a kilopascal), the real time is several times longer.

Pressure measurement with the wireline formation tester →
Supercharging as a pressure divider
Δps=ps−pres=Δpover RresRmc+Rres,Rmc=μf hmckmc\Delta p_s = p_s - p_{\text{res}} = \Delta p_{\text{over}}\,\frac{R_{\text{res}}}{R_{\text{mc}} + R_{\text{res}}},\qquad R_{\text{mc}} = \frac{\mu_f\,h_{mc}}{k_{mc}}

SI

Holds when: The flow through the cake is steady: its thickness and permeability are constant (dynamic equilibrium, or slow growth in the static period). The resistance of the rock changes slowly with time (logarithmically), so the formula is applied to the moment of the measurement.

Supercharging and permeability from formation-tester pretests →
Supercharging by the model of filtrate injection
Δps=u2 rw μf2kf[u1u2ln⁡Td+tsts+ln⁡kf tsϕ μf ct rw2+0.809]\Delta p_s = \frac{u_2\,r_w\,\mu_f}{2k_f}\left[\frac{u_1}{u_2}\ln\frac{T_d + t_s}{t_s} + \ln\frac{k_f\,t_s}{\phi\,\mu_f\,c_t\,r_w^2} + 0.809\right]

SI; natural logarithm

Holds when: The cake sets the filtration rate: the supercharging is small compared with the overbalance (otherwise use the pressure divider above). In each stage the rate is constant. The flow of filtrate is single-phase (one fluid), in a uniform rock. The logarithmic approximation at the wall holds a few minutes after injection starts, even in a layer of 0.03 mD, long before logging. The filtrate spreads away from the wall horizontally (radially), as in a uniform formation. A thin tight layer between permeable parts partly drains into its neighbours in the vertical direction, and the supercharging in it is smaller than calculated.

Supercharging and permeability from formation-tester pretests →
Distortion of the gradient and of the contact by supercharging
Gapp=G+Δ(Δps)Δz,δz∗=Δps,o−Δps,wGw−GoG_{\text{app}} = G + \frac{\Delta(\Delta p_s)}{\Delta z},\qquad \delta z^* = \frac{\Delta p_{s,o} - \Delta p_{s,w}}{G_w - G_o}

SI (in practice kPa, kPa/m and m)

Holds when: The lines are fitted to the supercharged points as to valid ones. The supercharging is small compared with the pressure differences on which the classification of the fluids is based.

Supercharging and permeability from formation-tester pretests →
Pressure build-up after a pretest under spherical flow
p(Δt)=p∗−ms fs(Δt),fs=1Δt−1tp+Δt,ms=q μϕ μ ct4π3/2 ks3/2p(\Delta t) = p^* - m_s\,f_s(\Delta t),\qquad f_s = \frac{1}{\sqrt{\Delta t}} - \frac{1}{\sqrt{t_p + \Delta t}},\qquad m_s = \frac{q\,\mu\sqrt{\phi\,\mu\,c_t}}{4\pi^{3/2}\,k_s^{3/2}}

SI

Holds when: The rock is uniform (possibly anisotropic). The disturbance has moved away from the probe by more than a few hole radii. The flow to the probe is then a full sphere, because the hole is a thin impermeable tube. (For flow to a hemisphere at a flat wall, the slope would be twice as large.) The influence of tool storage has ended (see the section on tool storage). The disturbance has not yet reached the boundaries of the formation. One phase of known viscosity flows. With two rates of draw, fsf_s is built by superposition, by the same rules.

Supercharging and permeability from formation-tester pretests →
Response of an observation probe at a constant draw rate
Δpobs(t)=qμ4πkh Δz erfc⁡Δz4ηvt  ≈  qμ4πkh Δz−mst,ηv=kvϕμct\Delta p_{\text{obs}}(t) = \frac{q\mu}{4\pi k_h\,\Delta z}\,\operatorname{erfc}\frac{\Delta z}{\sqrt{4\eta_v t}} \;\approx\; \frac{q\mu}{4\pi k_h\,\Delta z} - \frac{m_s}{\sqrt t},\qquad \eta_v = \frac{k_v}{\phi\mu c_t}

SI; the approximate form holds when Δz/4ηvt<0.15\Delta z/\sqrt{4\eta_v t} < 0.15

Holds when: A point sink in a uniform anisotropic medium; the hole is not allowed for. The hole slightly increases the response (by no more than a factor of two, the limit of a flat wall), and then the first formula underestimates khk_h. In processing real data, the constant is taken from a model of the tool with a hole. The rate at the sink is constant from the start; the flow is single-phase.

Supercharging and permeability from formation-tester pretests →
Relation between the permeability of core to gas and to formation water at reservoir conditions
kb=c ka n⟺lg⁡kb=lg⁡c+nlg⁡kak_b = c\,k_a^{\,n}\qquad\Longleftrightarrow\qquad \lg k_b = \lg c + n\lg k_a

mD; the constant cc refers specifically to millidarcies

Holds when: The relation is found on samples of the same rock type, and it is applied within the permeabilities of these pairs. It already includes gas slippage (gas slipping along the pore walls, which makes gas permeability read too high; Klinkenberg correction), the stress and the swelling of clays. There is no need to separate these corrections if the pairs were measured in this way.

Comparing formation-tester, core and well-test data →
Horizontal absolute permeability from the drawdown mobility
kh=(k/μ)d μfkrf′ fA,fA=kdkh=π2K(1−1/A)k_h = \frac{(k/\mu)_d\,\mu_f}{k_{rf}'\,f_A},\qquad f_A = \frac{k_d}{k_h} = \frac{\pi}{2K(1 - 1/A)}

khk_h in mD, for (k/μ)d(k/\mu)_d in mD/(mPa·s) and μf\mu_f in mPa·s

Holds when: The flushed zone is at the residual saturation of the displaced phase, and there is no extra resistance at the opening of the probe (cake, fine particles). Otherwise khk_h comes out too low. fAf_A comes from the model of a round probe on a flat impermeable wall (derivation). Solutions with an ellipsoidal sink give values about 10% smaller. AA is taken from pairs of horizontal and vertical samples of the same layer. At the scale of the probe, lamination (thin layering) can give a larger anisotropy than that of the samples.

Comparing formation-tester, core and well-test data →
Flow capacity from the Horner line, converted to absolute permeability
(kh)abs=(kh)okro′=2.149 qBμm kro′(kh)_{\text{abs}} = \frac{(kh)_o}{k_{ro}'} = \frac{2.149\,qB\mu}{m\,k_{ro}'}

metric oilfield units; 2.149=ln⁡10/(4π⋅0.0853)2.149 = \ln 10/(4\pi\cdot 0.0853)

Holds when: Only oil flows in the formation around the well, with irreducible water (the pressure is above the bubble point pressure, and the water is immobile). The slope is taken from the period of radial flow in an infinite reservoir (Diagnosis with the pressure derivative). The whole interval is open, and the layers are continuous within the radius of investigation.

Comparing formation-tester, core and well-test data →
Systematic bias and scatter of the ratios
di=lg⁡kpt,ikcore,i,dˉ=1n∑idi,s=∑i(di−dˉ)2n−1,dˉ±t snd_i = \lg\frac{k_{\text{pt},i}}{k_{\text{core},i}},\qquad \bar d = \frac{1}{n}\sum_i d_i,\qquad s = \sqrt{\frac{\sum_i (d_i - \bar d)^2}{n - 1}},\qquad \bar d \pm t\,\frac{s}{\sqrt n}

dimensionless (decimal logarithms)

Holds when: The deviations of the points are independent and close to lognormal (their logarithms follow the usual bell-shaped normal distribution). Points with a known defect of the measurement are excluded beforehand, not because of the size of their deviation. 10dˉ10^{\bar d} is the geometric mean of the ratios, a systematic multiplier. 10s10^{s} is the typical scatter of one point, as “by how many times”.

Comparing formation-tester, core and well-test data →
Number of points needed to tell bias from scatter
n≈(z slg⁡F)2n \approx \left(\frac{z\,s}{\lg F}\right)^2

dimensionless

Holds when: The condition is ∣dˉ∣>z s/n|\bar d| > z\,s/\sqrt n at dˉ=lg⁡F\bar d = \lg F. For small nn, instead of zz we take Student's coefficient for n−1n - 1 degrees of freedom. It is noticeably larger than two (3.18 for four points), and more points are needed. This condition is for a mean equal to lg⁡F\lg F. Because dˉ\bar d itself scatters, a bias of exactly a factor of FF will be significant in only about half of the sets. To notice it in four sets out of five, we take zz as about 2.8, and about twice as many points are needed. The points are independent: two points in one core window are in effect one.

Comparing formation-tester, core and well-test data →
The pressure that the formation tester records at the oil–water contact
pmeas(z)={po(z)=pFWL+Go (z−z∗),z<zOWCpw(z)=pFWL+Gw (z−z∗),z>zOWCp_{\text{meas}}(z) = \begin{cases} p_o(z) = p_{\text{FWL}} + G_o\,(z - z^*), & z < z_{\text{OWC}} \\ p_w(z) = p_{\text{FWL}} + G_w\,(z - z^*), & z > z_{\text{OWC}} \end{cases}

SI (in practice MPa, kPa/m and m)

Holds when: The accumulation is in hydrostatic and capillary equilibrium, and the formation water does not flow. The rock is water-wet, and the oil in the transition zone is continuous. In an oil-wet rock the picture is the mirror image (advanced section). The measurements are valid ones, without supercharging and without an unfinished build-up (classes of pretests). The well is drilled with water-based mud.

Pressure gradients, fluid contacts and capillary equilibrium →
Depth of the crossing of two pressure lines
z∗=z0+pa(z0)−pb(z0)Gb−Gaz^* = z_0 + \frac{p_a(z_0) - p_b(z_0)}{G_b - G_a}

any consistent system (MPa, MPa/m and m, or kPa, kPa/m and m)

Holds when: Both lines belong to one hydraulically connected system in equilibrium. The densities of the phases are constant with height. If the composition changes noticeably with depth, the line of the oil or gas is curved (section on fluid systems).

Pressure gradients, fluid contacts and capillary equilibrium →
Standard error of a contact from two lines
σz∗=σa2(z∗)+σb2(z∗)∣Gb−Ga∣,σ2(z)=s2[1n+(z−zˉ)2∑i(zi−zˉ)2]\sigma_{z^*} = \frac{\sqrt{\sigma_a^2(z^*) + \sigma_b^2(z^*)}}{|G_b - G_a|},\qquad \sigma^2(z) = s^2\left[\frac1n + \frac{(z - \bar z)^2}{\sum_i (z_i - \bar z)^2}\right]

any consistent system (in practice kPa, kPa/m and m)

Holds when: The errors of the points are random, independent and of equal size. The points within each line belong to one phase. Systematic errors are not included. A constant error of the gauge in one survey shifts both lines in the same way and does not change the crossing. A common shift of depths moves the crossing together with the logs of the same well. For lines from different wells, these errors are added separately (section on the regional water line). With few points, ss is estimated only roughly. The confidence interval (the range in which the true value is expected to lie) is taken with Student's coefficient for n−2n - 2 degrees of freedom, and not just ±2σ (about the standard error).

Pressure gradients, fluid contacts and capillary equilibrium →
Density and gradient of formation water from analysis
ρw,res=ρw,scBw,Gw=ρw,res g\rho_{w,\text{res}} = \frac{\rho_{w,\text{sc}}}{B_w},\qquad G_w = \rho_{w,\text{res}}\,g

SI; in metric form GwG_w [kPa/m] =9.81 ρw,res= 9.81\,\rho_{w,\text{res}} [g/cm³]

Holds when: Dissolved gas is ignored (otherwise the density is calculated from the mass of water and gas). Salinity and temperature are the same over the interval through which the line is extended.

Pressure gradients, fluid contacts and capillary equilibrium →
GOC from the oil line and the bubble point pressure
zGOC=zˉo+pb−pˉoGo,σzGOC≈σpbGoz_{\text{GOC}} = \bar z_o + \frac{p_{\text{b}} - \bar p_o}{G_o},\qquad \sigma_{z_{\text{GOC}}} \approx \frac{\sigma_{p_{\text{b}}}}{G_o}

any consistent system

Holds when: The accumulation is in equilibrium, and the gas cap and the oil are connected. The composition of the oil is constant with height. With compositional grading, the bubble point pressure of a sample taken below the GOC is smaller than that of the oil at the contact, and the GOC from this formula comes out above the true one.

Pressure gradients, fluid contacts and capillary equilibrium →
A pressure line with a set gradient and its uncertainty
p(z)=pˉ+ρg (z−zˉ),σ2(z)=s2n+(g σρ)2 (z−zˉ)2p(z) = \bar p + \rho g\,(z - \bar z),\qquad \sigma^2(z) = \frac{s^2}{n} + (g\,\sigma_\rho)^2\,(z - \bar z)^2

SI (in practice MPa, kPa/m and m); gg = 9.81 m/s²

Holds when: All the chosen points belong to one continuous phase, are valid, and have no noticeable supercharging. The density of the sample refers to the fluid of this line, and it is constant over the interval of extrapolation. The errors of the points and the error of the density are independent. With one point, the line passes through it, and ss is taken from other measurements.

Forced-gradient technique and multi-well analysis →
What the shifts of the lines of a well are made of
δo=δ+Go Δz+Δo,δw=δ+Gw Δz+Δw\delta_o = \delta + G_o\,\Delta z + \Delta_o,\qquad \delta_w = \delta + G_w\,\Delta z + \Delta_w

SI (in practice kPa, kPa/m and m)

Holds when: The lines of the field and of the well are single-phase, in equilibrium, with the same gradients. The errors are constant over the well. A depth error that grows with depth changes the slope of the lines too.

Forced-gradient technique and multi-well analysis →
How shifts of the lines move the FWL of a well
δz∗=δo−δwGw−Go,δpFWL=δo+Go δz∗\delta z^* = \frac{\delta_o - \delta_w}{G_w - G_o},\qquad \delta p_{\text{FWL}} = \delta_o + G_o\,\delta z^*

SI (in practice kPa, kPa/m and m)

Holds when: The gradients of the phases at the well are the same as in the field. A depth error Δz\Delta z shifts the level vertically by −Δz-\Delta z at the same pressure. It is allowed for separately.

Forced-gradient technique and multi-well analysis →
Weighted mean level and a check of consistency
zˉ∗=∑izi∗/σi2∑i1/σi2,σzˉ∗=1∑i1/σi2,χ2=∑i(zi∗−zˉ∗σi)2\bar z^* = \frac{\sum_i z_i^*/\sigma_i^2}{\sum_i 1/\sigma_i^2},\qquad \sigma_{\bar z^*} = \frac{1}{\sqrt{\sum_i 1/\sigma_i^2}},\qquad \chi^2 = \sum_i\left(\frac{z_i^* - \bar z^*}{\sigma_i}\right)^2

any consistent system

Holds when: The errors of the wells are independent and correctly estimated. Errors common to all wells (the fluid density, systematic errors of one type of tool) are not reduced by this method. All the wells belong to one system. A well that does not fit the others is excluded only after the cause is explained.

Forced-gradient technique and multi-well analysis →
Free water level with a common oil line
z∗=zref+Φo−ΦwGw−Go,zB∗−zA∗=Φw,A−Φw,BGw−Goz^* = z_{\text{ref}} + \frac{\Phi_o - \Phi_w}{G_w - G_o},\qquad z^*_{\text{B}} - z^*_{\text{A}} = \frac{\Phi_{w,\text{A}} - \Phi_{w,\text{B}}}{G_w - G_o}

any consistent system (MPa, MPa/m and m)

Holds when: The oil is at rest and connected between the wells, so its datum pressure is one. At the contact in each well the water is continuous. The densities of the phases are constant. For gas, the gas gradient is taken instead of GoG_o. The depths are vertical, from one level. The pressures are absolute, from calibrated gauges.

Perched and tilted contacts, trapped water →
Excess of the pressure of perched water over the main aquifer
Δpperch=(Gw−Go) (zmain∗−zperch∗)\Delta p_{\text{perch}} = (G_w - G_o)\,(z^*_{\text{main}} - z^*_{\text{perch}})

any consistent system

Holds when: The water of the block is closed on all sides, and since it was separated it has not leaked through seals and faults. The densities of the phases in the blocks are the same.

Perched and tilted contacts, trapped water →
Tilt of the contact from the gradient of the datum water pressure
tan⁡β=∣∇Φw∣Gw−Go=ρwρw−ρo ∣∇h∣\tan\beta = \frac{|\nabla\Phi_w|}{G_w - G_o} = \frac{\rho_w}{\rho_w - \rho_o}\,|\nabla h|

a gradient in kPa/km and pressure gradients in kPa/m give tan⁡β\tan\beta in m/km

Holds when: The oil is at rest, the flow of water is steady, and the densities are constant. The datum water pressures are compared on one datum and at one water density (otherwise, see the advanced section). Capillary forces are ignored: this is about the FWL. The OWC from logs is parallel to it, at the height hdh_d (transition zone).

Perched and tilted contacts, trapped water →
Height of the accumulation in a dome with a tilted contact
h(x,y)=Heff−ax (x−x0)2−ay (y−y0)2,x0=tx2ax,y0=ty2ay,Heff=H+tx24ax+ty24ay,V=π Heff22axay\begin{aligned} h(x, y) &= H_{\text{eff}} - a_x\,(x - x_0)^2 - a_y\,(y - y_0)^2, \\ x_0 &= \frac{t_x}{2a_x},\quad y_0 = \frac{t_y}{2a_y},\quad H_{\text{eff}} = H + \frac{t_x^2}{4a_x} + \frac{t_y^2}{4a_y}, \\ V &= \frac{\pi\,H_{\text{eff}}^2}{2\sqrt{a_x a_y}} \end{aligned}

coordinates in km, depths and heights in m

Holds when: The top is an elliptic paraboloid with axes along xx and yy. The formation is thicker than the accumulation (there is bottom water), so the rock volume is the volume between the top and the contact. The contact is a plane. The trap does not spill: the outline of the accumulation does not reach the spill point (the lowest point through which oil can leave).

Perched and tilted contacts, trapped water →
Pressure gradient with vertical flow
 dp dz=ρg+μ uzkz\frac{\dd p}{\dd z} = \rho g + \frac{\mu\,u_z}{k_z}

SI

Holds when: single-phase flow; the velocity uzu_z and the permeability kzk_z are constant over the measured interval, so the profile is a straight line if kzk_z changes with depth, the viscous component changes too: a thin tight layer gives a step (see the section on barriers) depths are true vertical depths, and the fluid density is known independently (from the analysis of a sample)

Pressure measurements in producing reservoirs: barriers and accuracy →
Significance of the deviation from the hydrostatic gradient
Gvisc=G−ρg,σGvisc=sG2+(g σρ)2G_{\text{visc}} = G - \rho g,\qquad \sigma_{G_{\text{visc}}} = \sqrt{s_G^2 + (g\,\sigma_\rho)^2}

SI (in practice kPa/m)

Holds when: the errors of the slope and of the density are random and independent; the points of the line are valid and belong to one unit a systematic scale error of the depths (a wrong directional survey) is not in the formula and is checked separately we call the flow detected if ∣Gvisc∣|G_{\text{visc}}| is larger than two or three of its errors. With few points the error of the slope is only a rough estimate, so the threshold is taken from Student's t coefficient (a factor from statistical tables that widens the threshold when there are few points)

Pressure measurements in producing reservoirs: barriers and accuracy →
Potential step across a barrier
δΦb=μ uz bkb\delta\Phi_b = \frac{\mu\,u_z\,b}{k_b}

SI

Holds when: the barrier is uniform and continuous; its storage capacity (how much fluid it gives up as the pressure falls) is small, so the flow through it is almost steady-state the units above and below the barrier are in vertical equilibrium and are filled with one fluid

Pressure measurements in producing reservoirs: barriers and accuracy →
Material balance of two units with crossflow
ctV1  dΦ1 dt=−q1−q12,ctV2  dΦ2 dt=−q2+q12,q12=T (Φ1−Φ2)\begin{aligned} c_tV_1\,\frac{\dd\Phi_1}{\dd t} &= -q_1 - q_{12}, \\ c_tV_2\,\frac{\dd\Phi_2}{\dd t} &= -q_2 + q_{12}, \qquad q_{12} = T\,(\Phi_1 - \Phi_2) \end{aligned}

metric oilfield: MPa, m³, m³/d, 1/MPa, days

Holds when: the potential is the same everywhere inside each unit: the unit is much more permeable than the barrier, and the drawdown around the wells (the extra pressure drop that drives fluid into a well) is ignored the withdrawals are constant and the compressibility is constant; inflow from outside (edge water) is counted as a negative withdrawal the storage capacity of the barrier itself is small; the flow is single-phase

Pressure measurements in producing reservoirs: barriers and accuracy →
Steady-state step and time constant of the barrier
Δss=(r2−r1) τ=q2V1−q1V2T (V1+V2),ri=qictVi,τ=ct V1V2(V1+V2) T\begin{aligned} \Delta_{\text{ss}} &= (r_2 - r_1)\,\tau = \frac{q_2V_1 - q_1V_2}{T\,(V_1 + V_2)}, \\ r_i &= \frac{q_i}{c_t V_i}, \qquad \tau = \frac{c_t\,V_1V_2}{(V_1 + V_2)\,T} \end{aligned}

metric oilfield, as above

Holds when: constant withdrawals and equal initial potentials of the units the compressibility is not in the steady-state step: it cancels in the right-hand form

Pressure measurements in producing reservoirs: barriers and accuracy →
Potential profile in a tight unit that flows into its neighbour
Φ(ζ)−Φ1=ΔΦmax⁡[2ζh−(ζh)2],ΔΦmax⁡=μ u h2 kz,u=qA⋅V2V1+V2\begin{aligned} \Phi(\zeta) - \Phi_1 &= \Delta\Phi_{\max}\left[\frac{2\zeta}{h} - \left(\frac{\zeta}{h}\right)^2\right], \\ \Delta\Phi_{\max} &= \frac{\mu\,u\,h}{2\,k_z}, \qquad u = \frac{q}{A}\cdot\frac{V_2}{V_1 + V_2} \end{aligned}

SI

Holds when: the top of the tight unit is impermeable, the flow in it is vertical, and no fluid flows to the wells along the tight unit itself common pseudo-steady-state regime: both units lose pressure at the same rate (the time criterion is given below) single-phase flow, and kzk_z is the same over the thickness of the unit

Pressure measurements in producing reservoirs: barriers and accuracy →
When the profile of the tight unit settles
tDh=kz tϕ μ ct h2>2t_{Dh} = \frac{k_z\,t}{\phi\,\mu\,c_t\,h^2} > 2

SI

Holds when: the withdrawal from the permeable unit is roughly constant; all of it goes through that unit, and the crossflow from the tight unit builds up from zero the limit holds for any ratio of the pore volumes of the units. If a constant rate were taken from the tight unit itself from the start, tDh>0.5t_{Dh} > 0.5 would be enough (see the advanced section)

Pressure measurements in producing reservoirs: barriers and accuracy →
Lag of the unit behind the barrier during a pulse
τloc=ϕ1h1 ct μ bkb,Δp2−Δp1≈τloc  d(Δp2) dt\tau_{\text{loc}} = \frac{\phi_1 h_1\,c_t\,\mu\,b}{k_b},\qquad \Delta p_2 - \Delta p_1 \approx \tau_{\text{loc}}\,\frac{\dd(\Delta p_2)}{\dd t}

SI

Holds when: the well that creates the pulse is open only to unit 2, and unit 1 receives the disturbance only through the barrier lateral flow inside unit 1 is ignored; this overestimates the lag if unit 1 is permeable across the area Δp2\Delta p_2 is the pressure of unit 2 at the observation well, including the crossflow. The barrier feeds fluid into unit 2 over the whole area that the pulse has reached. The leakage factor is B=k2h2b/kbB = \sqrt{k_2h_2b/k_b} (the distance over which the shale can feed unit 2; k2k_2 and h2h_2 are the permeability and thickness of unit 2). If BB is smaller than the distance to the well that creates the pulse, the line-source solution for unit 2 alone strongly overestimates its response (see the example below). The line-source solution is the standard formula for the pressure around a well in an infinite reservoir the approximate equality holds for an almost linear change of the pressure of unit 2 and for times much longer than τloc\tau_{\text{loc}}; the exact relation is Δp2−Δp1=τloc  d(Δp1)/ dt\Delta p_2 - \Delta p_1 = \tau_{\text{loc}}\,\dd(\Delta p_1)/\dd t

Pressure measurements in producing reservoirs: barriers and accuracy →
Potentials of the two units in time
Φˉ(t)=Φ0−(q1+q2) tct (V1+V2),Δ(t)=Δss(1−e−t/τ),Φ1=Φˉ+V2V1+V2 Δ,Φ2=Φˉ−V1V1+V2 Δ\begin{aligned} \bar\Phi(t) &= \Phi_0 - \frac{(q_1 + q_2)\,t}{c_t\,(V_1 + V_2)}, \qquad \Delta(t) = \Delta_{\text{ss}}\left(1 - e^{-t/\tau}\right), \\ \Phi_1 &= \bar\Phi + \frac{V_2}{V_1 + V_2}\,\Delta, \qquad \Phi_2 = \bar\Phi - \frac{V_1}{V_1 + V_2}\,\Delta \end{aligned}

metric oilfield, as in the main section

Holds when: constant withdrawals, equal initial potentials, and the other assumptions of the two-unit model.

Pressure measurements in producing reservoirs: barriers and accuracy →
Response of two units with crossflow to a change in the rate of a well in unit 2
pˉ2(s)=2⋅0.933 qμk2h2 s K0 ⁣(rBτ2s+τloc s1+τloc s),pˉ1(s)=pˉ2(s)1+τloc s\bar p_2(s) = \frac{2\cdot 0.933\,q\mu}{k_2h_2\,s}\,K_0\!\left(\frac{r}{B}\sqrt{\tau_2 s + \frac{\tau_{\text{loc}}\,s}{1 + \tau_{\text{loc}}\,s}}\right),\qquad \bar p_1(s) = \frac{\bar p_2(s)}{1 + \tau_{\text{loc}}\,s}

metric oilfield: MPa, m³/d, mD, m, mPa·s, days

Holds when: unit 1 does not carry flow across the area (lateral flow in it is ignored), the storage capacity of the shale is small, and the formation is uniform across the area and unbounded the pressures in time are found by numerical inversion of the Laplace transform (turning the result in ss back into time; Stehfest algorithm) without crossflow (kb→0k_b \to 0) the formula turns into the line-source solution (Line-source solution and the pressure function)

Pressure measurements in producing reservoirs: barriers and accuracy →
Pressure gradient with vertical flow of two phases
 dp dz=gρˉ+utλt,ρˉ=λ1ρ1+λ2ρ2λ1+λ2,λj=kz krjμj\frac{\dd p}{\dd z} = g\bar\rho + \frac{u_t}{\lambda_t},\qquad \bar\rho = \frac{\lambda_1\rho_1 + \lambda_2\rho_2}{\lambda_1 + \lambda_2},\qquad \lambda_j = \frac{k_z\,k_{rj}}{\mu_j}

SI

Holds when: one-dimensional vertical flow in a uniform unit; the phases are incompressible over the measured interval the pressures of the phases are equal, so capillary pressure is ignored (the correction is in the section on capillary pressure) the formula is local: the mobilities are taken at the saturation of the point in question. If the saturation changes with depth, the gradient changes with it, and the profile bends

Two-phase flow and reservoir description from formation-tester data →
Density from the gradient in countercurrent flow of gas and oil
ρˉ=ρo+Mρg1+M,M=krg/μgkro/μo\bar\rho = \frac{\rho_o + M\rho_g}{1 + M},\qquad M = \frac{k_{rg}/\mu_g}{k_{ro}/\mu_o}

densities in any identical units; the gradient is gρˉg\bar\rho

Holds when: there is no total vertical flow: the gas rises and the oil drains down at the same volumetric rate the gas–oil capillary pressure is negligible (the interfacial tension between the two fluids is low); the water is immobile at irreducible saturation (the water that stays in the pores and cannot be displaced)

Two-phase flow and reservoir description from formation-tester data →
Gradient of the oil pressure in a two-phase zone with the capillary term
Goil= dpo dz=G+fw  dpc dz,fw=λwλtG_{\text{oil}} = \frac{\dd p_o}{\dd z} = G + f_w\,\frac{\dd p_c}{\dd z},\qquad f_w = \frac{\lambda_w}{\lambda_t}

SI (in practice kPa/m)

Holds when: a water-wet rock, with continuous oil in the zone; the capillary pressure curve is the imbibition curve (the curve for water entering the rock) if the zone is being flooded (imbibition profile) the other assumptions of the two-phase gradient; the derivation is in the advanced section

Two-phase flow and reservoir description from formation-tester data →
Shift of the contact if the water line is set to the hydrostatic gradient
δz∗=−Gvisc LGw−Go\delta z^* = -\frac{G_{\text{visc}}\,L}{G_w - G_o}

any consistent units (kPa/m and m)

Holds when: the line with the set gradient passes through the middle point of the water readings, and the true oil and water lines are straight. Then the formula is exact. With the measured, steeper water gradient in the denominator, it would underestimate the shift the oil line is found from its own points, without a shift

Two-phase flow and reservoir description from formation-tester data →
Water that has entered the accumulation, from the rise of the contact
We−WpBw≈A ϕ (1−Swi−Sor) ΔzW_e - W_pB_w \approx A\,\phi\,(1 - S_{wi} - S_{or})\,\Delta z

SI (m, m², m³)

Holds when: the water rises evenly over the whole area AA: bottom water, vertical saturation equilibrium. With cones near wells or inflow from the flanks, the rise in one well is not representative the swept zone stays at residual oil saturation; the change in the shape of the transition zone and the compressibility of the pores and the water are ignored the initial and the current FWL are both found from pressure, with common depth matching of the wells

Two-phase flow and reservoir description from formation-tester data →
Material balance of a tank in the multi-tank model
Cj  dpj dt=−qj+∑kTjk (pk−pj),Cj=ct Vp,jC_j\,\frac{\dd p_j}{\dd t} = -q_j + \sum_k T_{jk}\,(p_k - p_j),\qquad C_j = c_t\,V_{p,j}

metric oilfield: MPa, m³, m³/d, days

Holds when: the pressure inside a tank is the same: it equalises much faster than between tanks single-phase undersaturated oil or water, with constant compressibility and viscosity; below the bubble point pressure the storage capacity depends on pressure, and the equations are nonlinear (see the advanced section) the crossflow through a link is proportional to the pressure difference of the tanks: the flow through it is steady-state, and the storage capacity of the link itself is ignored

Lumped multi-tank material balance model →
Rate of a well in a small block behind a fault
q(t)=q∞+(q0−q∞) e−t/τ,τ=CJ+T,q0=J (pnb−pwf),q∞=J TJ+T (pnb−pwf)\begin{aligned} q(t) &= q_\infty + (q_0 - q_\infty)\,\e^{-t/\tau},\qquad \tau = \frac{C}{J + T}, \\ q_0 &= J\,(p_{\text{nb}} - p_{wf}),\qquad q_\infty = \frac{J\,T}{J + T}\,(p_{\text{nb}} - p_{wf}) \end{aligned}

metric oilfield

Holds when: the block is a tank; the pressure of the neighbouring block is constant; at the start the pressure of the block equals that of its neighbour the bottomhole pressure and the productivity are constant; the flow is single-phase, above the bubble point pressure

Lumped multi-tank material balance model →
Time constants of three tanks in a chain
λ2−s1λ+s2=0,τ=1/λ,s1=TABCA+TAB+TBCCB+TBCCC,s2=TAB TBC (CA+CB+CC)CA CB CC\begin{aligned} &\lambda^2 - s_1\lambda + s_2 = 0,\qquad \tau = 1/\lambda, \\ &s_1 = \frac{T_{AB}}{C_A} + \frac{T_{AB} + T_{BC}}{C_B} + \frac{T_{BC}}{C_C},\qquad s_2 = \frac{T_{AB}\,T_{BC}\,(C_A + C_B + C_C)}{C_A\,C_B\,C_C} \end{aligned}

metric oilfield

Holds when: a chain A–B–C with no direct A–C link; the storage capacities and transmissibilities are constant.

Lumped multi-tank material balance model →

Symbols

SymbolMeaning and context
zˉ\bar z
Δz\Delta z
∇h\nabla h
dˉ, s\bar d,\ s
zˉ, pˉ\bar z,\ \bar p
(k/μ)d(k/\mu)_d
(kh)o(kh)_o
δz∗\delta z^*
zˉ∗, σzˉ∗\bar z^*,\ \sigma_{\bar z^*}
 dp/ dz\dd p/\dd z
∇Φw\nabla\Phi_w
0.8090.809
pˉ2, pˉ1\bar p_2,\ \bar p_1
Δp2, Δp1\Delta p_2,\ \Delta p_1
AA
ax, aya_x,\ a_y
bb
BB
BwB_w
b, kbb,\ k_b
CC
CA, CB, CCC_A,\ C_B,\ C_C
CjC_j
ctc_t
ct, Vp,jc_t,\ V_{p,j}
c, nc,\ n
 dpc/ dz\dd p_c/\dd z
Δpd\Delta p_d
did_i
FF
fAf_A
fwf_w
δpFWL\delta p_{\text{FWL}}
gg
GG
Ga, GbG_a,\ G_b
Gapp, GG_{\text{app}},\ G
GoG_o
Go, GwG_o,\ G_w
GoilG_{\text{oil}}
GviscG_{\text{visc}}
GwG_w
Gw−GoG_w - G_o
Gw, GoG_w,\ G_o
G, sGG,\ s_G
hh
HH
HeffH_{\text{eff}}
h(x,y)h(x, y)
JJ
kk
KK
K0K_0
k2, h2k_2,\ h_2
kak_a
kbk_b
kfk_f
kpt,i, kcore,ik_{\text{pt},i},\ k_{\text{core},i}
krf′k_{rf}'
  • end-point relative permeability of the filtrate in the flushed zone: krw′k_{rw}' for a water-based filtrate in the oil-bearing part, 1 in the water-bearing part, kro′k_{ro}' for a hydrocarbon filtrate, —Comparing formation-tester, core and well-test data
krg, krok_{rg},\ k_{ro}
kro′k_{ro}'
ksk_s
kzk_z
kz, krj, μjk_z,\ k_{rj},\ \mu_j
K(m)K(m)
k/μk/\mu
LL
mm
MM
msm_s
nn
zˉo, pˉo\bar z_o,\ \bar p_o
Δpobs\Delta p_{\text{obs}}
Δpover\Delta p_{\text{over}}
pp
p∗p^*
pa(z0), pb(z0)p_a(z_0),\ p_b(z_0)
pbp_{\text{b}}
pjp_j
pmeas(z)p_{\text{meas}}(z)
pnbp_{\text{nb}}
pwfp_{wf}
p(Δt)p(\Delta t)
Δpperch\Delta p_{\text{perch}}
qq
q1, q2q_1,\ q_2
q12q_{12}
qjq_j
qpistq_{\text{pist}}
q, Aq,\ A
q, B, μq,\ B,\ \mu
q, q0, q∞q,\ q_0,\ q_\infty
q, μ, ϕ, ct, kh, kvq,\ \mu,\ \phi,\ c_t,\ k_h,\ k_v
rr
r1, r2r_1,\ r_2
RmcR_{\text{mc}}
rper_{pe}
RresR_{\text{res}}
rwr_w
ss
Δps\Delta p_s
SS
s1s_1
s2s_2
sGs_G
Swi, SorS_{wi},\ S_{or}
s, n, zˉs,\ n,\ \bar z
Δps,o, Δps,w\Delta p_{s,o},\ \Delta p_{s,w}
tt
TT
TAB, TBCT_{AB},\ T_{BC}
TdT_d
tDht_{Dh}
TjkT_{jk}
  • transmissibility of the link between tanks jj and kk: the crossflow per 1 MPa of pressure difference, Tjk=TkjT_{jk} = T_{kj}. For a section of a fault it is 0.0853 (kf/lf) Af/μ0.0853\,(k_f/l_f)\,A_f/\mu, m³/(d·MPa)Lumped multi-tank material balance model
tpt_p
tst_s
tx, tyt_x,\ t_y
uu
u1, u2u_1,\ u_2
utu_t
uzu_z
VV
V1, V2V_1,\ V_2
WeW_e
WpBwW_pB_w
x0, y0x_0,\ y_0
zz
z∗z^*
zmain∗z^*_{\text{main}}
zperch∗z^*_{\text{perch}}
z∗, zA∗, zB∗z^*,\ z^*_{\text{A}},\ z^*_{\text{B}}
z∗, pFWLz^*,\ p_{\text{FWL}}
z0z_0
ziz_i
zi∗, σiz_i^*,\ \sigma_i
zi, piz_i,\ p_i
zOWCz_{\text{OWC}}
zrefz_{\text{ref}}
β\beta
δ\delta
Δ, Δss, τ\Delta,\ \Delta_{\text{ss}},\ \tau
Δ(Δps)/Δz\Delta(\Delta p_s)/\Delta z
δo, δw\delta_o,\ \delta_w
Δo, Δw\Delta_o,\ \Delta_w
Δss\Delta_{\text{ss}}
δΦb\delta\Phi_b
ζ\zeta
ηv\eta_v
λ\lambda
λj, λt\lambda_j,\ \lambda_t
μ\mu
μ, ϕ, ct\mu,\ \phi,\ c_t
μf\mu_f
μf, hmc, kmc\mu_f,\ h_{mc},\ k_{mc}
μf, ϕ, ct\mu_f,\ \phi,\ c_t
μg, μo\mu_g,\ \mu_o
ρ\rho
ρˉ\bar\rho
ρ, σρ\rho,\ \sigma_\rho
ρ1, ρ2\rho_1,\ \rho_2
ρo, ρg\rho_o,\ \rho_g
ρw, ρo\rho_w,\ \rho_o
ρw,sc\rho_{w,\text{sc}}
σ(z)\sigma(z)
σa(z∗), σb(z∗)\sigma_a(z^*),\ \sigma_b(z^*)
σGvisc\sigma_{G_{\text{visc}}}
σgauge\sigma_{\text{gauge}}
σp\sigma_p
σpb\sigma_{p_{\text{b}}}
σz\sigma_z
σz∗\sigma_{z^*}
σzGOC\sigma_{z_{\text{GOC}}}
τ\tau
τ2, τloc\tau_2,\ \tau_{\text{loc}}
τloc\tau_{\text{loc}}
ϕ\phi
Φˉ\bar\Phi
ϕ, ct, h\phi,\ c_t,\ h
Φ0\Phi_0
Φ1\Phi_1
ϕ1, h1\phi_1,\ h_1
Φ1, Φ2\Phi_1,\ \Phi_2
Φo\Phi_o
Φw, Φw,A, Φw,B\Phi_w,\ \Phi_{w,\text{A}},\ \Phi_{w,\text{B}}
χ2\chi^2