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

Equations: Aquifer influx

The equations of the “Aquifer influx” lessons in “Reservoir engineering”, with their conditions and units, and what each symbol means in them.

Elastic store of an aquifer
We,max⁡=(cw+cf) W Δp,W=π (ra2−ro2) h ϕ fW_{e,\max} = (c_w + c_f)\,W\,\Delta p,\qquad W = \pi\,(r_a^2 - r_o^2)\,h\,\phi\,f

consistent units: m³, MPa and 1/MPa

Holds when: the aquifer is full of water, and the compressibilities are constant the geometry is radial, a sector of a ring with the angle θ\theta; for a strip-shaped aquifer of width ww and length LL the volume is wLhϕwLh\phi the pressure drop has had time to reach all parts of the aquifer

Aquifer models: water influx →
Influx from a small aquifer
We=(cw+cf) W (pi−p)W_e = (c_w + c_f)\,W\,(p_i - p)

m³, MPa, 1/MPa (or barrels, psi and 1/psi)

Holds when: the pressure throughout the aquifer equals the pressure at the contact: the aquifer equalises faster than the pressure in the accumulation changes there is no recharge from outside the pressure at the contact is taken as equal to the average pressure of the accumulation

Aquifer models: water influx →
Dimensionless time of the aquifer and the time to reach its boundary
tD=η tro2=k tϕ μw ct ro2,ta≈(ra−ro)24ηt_D = \frac{\eta\,t}{r_o^2} = \frac{k\,t}{\phi\,\mu_w\,c_t\,r_o^2},\qquad t_a \approx \frac{(r_a - r_o)^2}{4\eta}

SI

Holds when: a uniform aquifer of constant thickness, with radial geometry tat_a is an order-of-magnitude estimate: a noticeable pressure drop at the boundary starts earlier, and the equalisation of the whole aquifer comes several times later

Aquifer models: water influx →
Steady-state model of water influx
 dWe dt=Ja (pi−p),We=Ja∫0t(pi−p) dt\frac{\dd W_e}{\dd t} = J_a\,(p_i - p),\qquad W_e = J_a\int_0^t (p_i - p)\dd t

m³ at reservoir conditions, days, MPa; in oilfield units, barrels per day per psi

Holds when: the pressure at the outer boundary of the aquifer is constant and equal to the initial pressure the flow in the aquifer is steady: the transient period is short compared with the development history the pressure at the contact is taken as equal to the average pressure of the accumulation

Aquifer models: water influx →
Productivity index of a radial aquifer (steady flow)
Ja=2π f k hμw ln⁡(ra/ro)J_a = \frac{2\pi\,f\,k\,h}{\mu_w\,\ln(r_a/r_o)}

SI

Holds when: a radial aquifer of constant thickness and permeability the pressure at the radius rar_a is constant; the flow is steady and single-phase

Aquifer models: water influx →
Pseudo-steady-state aquifer model
 dWe dt=J (pˉa−p),We=(cw+cf) W (pi−pˉa)\frac{\dd W_e}{\dd t} = J\,(\bar p_a - p),\qquad W_e = (c_w + c_f)\,W\,(p_i - \bar p_a)

m³ at reservoir conditions, days, MPa

Holds when: the aquifer is closed; after a transient period the pressure in it falls everywhere at the same rate the pressure at the contact is taken as equal to the average pressure of the accumulation

Aquifer models: water influx →
Influx for a constant pressure drop at the contact (unsteady-state model)
We=U Δp WD(tD),U=2π f ϕ ct h ro2W_e = U\,\Delta p\,W_D(t_D),\qquad U = 2\pi\,f\,\phi\,c_t\,h\,r_o^2

m³, MPa, 1/MPa, m

Holds when: a uniform radial aquifer; the pressure at the contact fell by Δp\Delta p in one step and does not change afterwards for a changing pressure the history is split into steps and the responses are added (Superposition principle)

Aquifer models: water influx →
Response of a pseudo-steady-state aquifer to a constant pressure drop
 dWe dt=J[Δp−We(cw+cf)W]  ⇒  We=(cw+cf) W Δp(1−e−t/τ),τ=(cw+cf) WJ\frac{\dd W_e}{\dd t} = J\left[\Delta p - \frac{W_e}{(c_w + c_f)W}\right]\;\Rightarrow\; W_e = (c_w + c_f)\,W\,\Delta p\left(1 - \e^{-t/\tau}\right),\quad \tau = \frac{(c_w + c_f)\,W}{J}

m³, days, MPa

Holds when: the pressure drop at the contact is constant; at t=0t = 0 the aquifer is untouched: We=0W_e = 0

Aquifer models: water influx →
The van Everdingen–Hurst problem in dimensionless form
1rD∂∂rD ⁣(rD ∂P∂rD)=∂P∂tD,P(1,tD)=1,P(rD,0)=0,∂P∂rD∣rD=raD=0\frac{1}{r_D}\frac{\partial}{\partial r_D}\!\left(r_D\,\frac{\partial P}{\partial r_D}\right) = \frac{\partial P}{\partial t_D},\qquad P(1, t_D) = 1,\quad P(r_D, 0) = 0,\quad \left.\frac{\partial P}{\partial r_D}\right|_{r_D = r_{aD}} = 0

dimensionless form

Holds when: a uniform aquifer of constant thickness, full of water; the compressibilities are small and constant, so the equation is linear the flow is radial, with no crossflow through the top and bottom the outer boundary is impermeable; for an infinite aquifer the condition on it is replaced by P→0P \to 0 as rD→∞r_D \to \infty

Unsteady-state water influx: the van Everdingen–Hurst method →
Influx for a constant pressure drop at the contact
We=U Δp WD(tD, raD),U=2πfϕcthro2,tD=k tϕ μwct ro2W_e = U\,\Delta p\,W_D(t_D,\,r_{aD}),\qquad U = 2\pi f\phi c_t h r_o^2,\qquad t_D = \frac{k\,t}{\phi\,\mu_w c_t\,r_o^2}

WeW_e and UU in one system of volume and pressure (m³ and MPa); tDt_D in SI

Holds when: the problem above: a uniform radial aquifer, with the contact pressure falling in one step and not changing afterwards the pressure at the contact is taken as equal to the average pressure of the accumulation, as in the material balance the aquifer is full of water, ct=cw+cfc_t = c_w + c_f, and the compressibilities are constant

Unsteady-state water influx: the van Everdingen–Hurst method →
Influx as the expansion of the pressure-drop funnel
We=ct∫Vp(pi−p) dVp=U Δp∫1raDP rD drD,WD=∫1raDP rD drDW_e = c_t\int_{V_p}(p_i - p)\dd V_p = U\,\Delta p\int_1^{r_{aD}} P\,r_D\dd r_D,\qquad W_D = \int_1^{r_{aD}} P\,r_D\dd r_D

m³, MPa, 1/MPa

Holds when: there is only water in the aquifer; ctc_t is constant; there is no recharge from outside

Unsteady-state water influx: the van Everdingen–Hurst method →
Influx for a stepwise pressure drop at the contact
We(tn)=U∑j=0n−1Δpj WD(tDn−tDj)W_e(t_n) = U\sum_{j=0}^{n-1}\Delta p_j\,W_D\big(t_{Dn} - t_{Dj}\big)

m³, MPa

Holds when: the conditions of linearity: the properties of the water and the aquifer do not depend on pressure the contact pressure is constant between the moments of the steps; how well the steps replace a smooth history depends on the choice of the drops and of the step (the next section) all the steps are counted from one initial state: before development the pressure everywhere in the aquifer equals pip_i

Unsteady-state water influx: the van Everdingen–Hurst method →
Step drops from the average pressures of the intervals
Δp0=p0−p12,Δpj=pj−1−pj+12(j≥1),∑j=0n−1Δpj=p0−pn−1+pn2\Delta p_0 = \frac{p_0 - p_1}{2},\qquad \Delta p_j = \frac{p_{j-1} - p_{j+1}}{2}\quad (j \ge 1),\qquad \sum_{j=0}^{n-1}\Delta p_j = p_0 - \frac{p_{n-1} + p_n}{2}

MPa

Holds when: if the intervals between measurements differ, the formula is the same, since the stepped curve still follows the average pressures of the intervals, but the longest intervals bring the largest error the pressure in an interval changes smoothly, without sharp jumps (the start-up or shutdown of many wells)

Unsteady-state water influx: the van Everdingen–Hurst method →
Rate of influx and dimensionless rate
qw=2πfkh roμw∂p∂r∣ro=2πfkh Δpμw qD(tD),qD=−∂P∂rD∣rD=1q_w = \frac{2\pi f k h\,r_o}{\mu_w}\left.\frac{\partial p}{\partial r}\right|_{r_o} = \frac{2\pi f k h\,\Delta p}{\mu_w}\,q_D(t_D),\qquad q_D = -\left.\frac{\partial P}{\partial r_D}\right|_{r_D = 1}

SI

Holds when: a constant pressure drop Δp\Delta p at the contact from the moment t=0t = 0

Unsteady-state water influx: the van Everdingen–Hurst method →
From the rate of influx to the cumulative influx
We=∫0tqw dt=2πfkh Δpμw⋅ϕμwctro2k∫0tDqD dtD=2πfϕcthro2 Δp WDW_e = \int_0^t q_w\dd t = \frac{2\pi fkh\,\Delta p}{\mu_w}\cdot\frac{\phi\mu_wc_tr_o^2}{k}\int_0^{t_D}q_D\dd t_D = 2\pi f\phi c_thr_o^2\,\Delta p\,W_D

SI

Holds when: the same as for the problem in dimensionless form

Unsteady-state water influx: the van Everdingen–Hurst method →
Laplace transform of the dimensionless influx
WˉD(s)=K1(s) I1(raDs)−I1(s) K1(raDs)s3/2[K0(s) I1(raDs)+I0(s) K1(raDs)]  →  raD→∞    K1(s)s3/2 K0(s)\bar W_D(s) = \frac{K_1(\sqrt s)\,I_1(r_{aD}\sqrt s) - I_1(\sqrt s)\,K_1(r_{aD}\sqrt s)}{s^{3/2}\big[K_0(\sqrt s)\,I_1(r_{aD}\sqrt s) + I_0(\sqrt s)\,K_1(r_{aD}\sqrt s)\big]}\;\xrightarrow{\;r_{aD}\to\infty\;}\;\frac{K_1(\sqrt s)}{s^{3/2}\,K_0(\sqrt s)}

dimensionless

Holds when: a closed outer boundary at raDr_{aD}; on the right is the infinite aquifer

Unsteady-state water influx: the van Everdingen–Hurst method →
Influx from a linear aquifer
We=UL Δp WD(tD),UL=w h L ϕ ct,tD=k tϕ μwct L2;WD=2tD/π  (infinite strip)W_e = U_L\,\Delta p\,W_D(t_D),\qquad U_L = w\,h\,L\,\phi\,c_t,\qquad t_D = \frac{k\,t}{\phi\,\mu_wc_t\,L^2};\qquad W_D = 2\sqrt{t_D/\pi}\ \ \text{(infinite strip)}

m³, MPa, 1/MPa, m; tDt_D in SI (the metric form has the multiplier 3.553⋅10−33.553\cdot10^{-3} for time in hours)

Holds when: a uniform strip of constant cross-section; the flow is perpendicular to the contact for an infinite strip LL cancels: We=2whϕctΔpηt/πW_e = 2wh\phi c_t\Delta p\sqrt{\eta t/\pi}

Unsteady-state water influx: the van Everdingen–Hurst method →
Influx of the unsteady-state aquifer model
We(tn)=U wn,wn=∑j=0n−1Δpj WD(tDn−tDj; raD),tD=k tϕ μwct ro2W_e(t_n) = U\,w_n,\qquad w_n = \sum_{j=0}^{n-1}\Delta p_j\,W_D\big(t_{Dn} - t_{Dj};\ r_{aD}\big),\qquad t_D = \frac{k\,t}{\phi\,\mu_wc_t\,r_o^2}

m³, MPa; tDt_D in SI, or with the multiplier 3.553⋅10−33.553\cdot10^{-3} for kk in mD, tt in hours, μw\mu_w in mPa·s and ctc_t in 1/MPa

Holds when: a uniform radial aquifer; the pressure at the contact equals the average pressure of the accumulation the function ww is calculated with short steps, or exactly; otherwise the error of the stepped representation of the pressure passes into the fitted parameters

History matching of aquifer parameters →
Apparent aquifer constant
Un=We,nwn=Fn−NEt,nwnU_n = \frac{W_{e,n}}{w_n} = \frac{F_n - N E_{t,n}}{w_n}

m³/MPa

Holds when: the oil in place is known independently; there is no gas cap, or its size is known (EtE_t includes mEgmE_g) the early dates, where the pressure drop is about the same size as the error of its measurement, get a smaller weight

History matching of aquifer parameters →
Havlena–Odeh line with an unsteady-state influx function
FEt=N+U wEt,Et=Eo+mEg+Efw\frac{F}{E_t} = N + U\,\frac{w}{E_t},\qquad E_t = E_o + mE_g + E_{fw}

m³, MPa, m³/MPa

Holds when: there is no gas cap, or its size mm is known the early points, where EtE_t is small, are unreliable and must not decide the slope

History matching of aquifer parameters →
Measure of agreement with the pressure history
M=1n∑k=1n(pkmeas−pkcalcσ)2,M≤1 ⟺ RMS=1n∑k(pkmeas−pkcalc)2≤σM = \frac1n\sum_{k=1}^{n}\left(\frac{p_k^{\text{meas}} - p_k^{\text{calc}}}{\sigma}\right)^{2},\qquad M \le 1\ \Longleftrightarrow\ \text{RMS} = \sqrt{\frac1n\sum_{k}\big(p_k^{\text{meas}} - p_k^{\text{calc}}\big)^2} \le \sigma

MM is dimensionless; pressures in MPa

Holds when: the measurement errors are independent and equal in size each candidate aquifer is fitted so that MM is the smallest for its raDr_{aD} and permeability; the constant UU is then fitted automatically

History matching of aquifer parameters →
The pseudo-steady-state aquifer model in the notation of the Fetkovich method
 dWe dt=J (pˉa−p),We=ctW (pi−pˉa)=Wei(1−pˉapi),Wei=ctW pi\frac{\dd W_e}{\dd t} = J\,\big(\bar p_a - p\big),\qquad W_e = c_tW\,\big(p_i - \bar p_a\big) = W_{ei}\left(1 - \frac{\bar p_a}{p_i}\right),\qquad W_{ei} = c_tW\,p_i

m³, days, MPa

Holds when: The aquifer is closed. After the transient period it drains so that the pressure falls at the same rate everywhere in it (the pseudo-steady-state regime). Compressibility and permeability are constant. The pressure at the contact equals the average pressure of the accumulation.

Simplified water influx models: Fetkovich and Carter–Tracy →
The Fetkovich step
ΔWe,n=ctW (pˉa,n−1−pˉn)(1−e−Δtn/τ),pˉa,n=pi−We,n−1+ΔWe,nctW,pˉn=pn−1+pn2\Delta W_{e,n} = c_tW\,\big(\bar p_{a,n-1} - \bar p_n\big)\left(1 - \e^{-\Delta t_n/\tau}\right),\qquad \bar p_{a,n} = p_i - \frac{W_{e,n-1} + \Delta W_{e,n}}{c_tW},\qquad \bar p_n = \frac{p_{n-1} + p_n}{2}

m³, days, MPa

Holds when: Within the step, the pressure at the contact is replaced by its average value. The pseudo-steady-state aquifer model (the equations above).

Simplified water influx models: Fetkovich and Carter–Tracy →
Productivity index of the aquifer for the Fetkovich method
Jrad=2πf k hμw(ln⁡raD−34),Jlin=3 k w hμw LJ_{\text{rad}} = \frac{2\pi f\,k\,h}{\mu_w\left(\ln r_{aD} - \dfrac34\right)},\qquad J_{\text{lin}} = \frac{3\,k\,w\,h}{\mu_w\,L}

SI

Holds when: Pseudo-steady-state regime: the pressure difference that drives the flow is the average aquifer pressure minus the pressure at the contact. The formula for a radial aquifer is correct for raD≫1r_{aD} \gg 1. For small raDr_{aD} use the exact expression for a ring (see the advanced section).

Simplified water influx models: Fetkovich and Carter–Tracy →
Influence function of a radial aquifer
pD(tD)=2πfkh (pi−po(t))qw μw;pD≈2tD/π  (tD≪1),pD≈12(ln⁡tD+0.809)  (infinite),pD≈2tDraD2−1+ln⁡raD−34  (closed, late)p_D(t_D) = \frac{2\pi f k h\,\big(p_i - p_o(t)\big)}{q_w\,\mu_w};\qquad p_D \approx 2\sqrt{t_D/\pi}\ \ (t_D \ll 1),\qquad p_D \approx \tfrac12\big(\ln t_D + 0.809\big)\ \ \text{(infinite)},\qquad p_D \approx \frac{2t_D}{r_{aD}^2 - 1} + \ln r_{aD} - \frac34\ \ \text{(closed, late)}

SI; a dimensionless function

Holds when: A uniform radial aquifer. The influx rate is constant from the moment t=0t = 0. There is no skin at the contact. The logarithmic form holds for tDt_D above about 25, while the aquifer still behaves as an infinite one. The last form holds after the pressure drop has reached the closed boundary, for raD≫1r_{aD} \gg 1 (for small raDr_{aD}, use the exact expression for a ring instead of ln⁡raD−3/4\ln r_{aD} - 3/4).

Simplified water influx models: Fetkovich and Carter–Tracy →
The Carter–Tracy step
We,n=We,n−1+(tDn−tD,n−1) U Δpn−We,n−1 pD′(tDn)pD(tDn)−tD,n−1 pD′(tDn),Δpn=pi−pnW_{e,n} = W_{e,n-1} + \big(t_{Dn} - t_{D,n-1}\big)\,\frac{U\,\Delta p_n - W_{e,n-1}\,p_D'(t_{Dn})}{p_D(t_{Dn}) - t_{D,n-1}\,p_D'(t_{Dn})},\qquad \Delta p_n = p_i - p_n

m³, MPa; time is dimensionless

Holds when: The influx rate is constant within the step. The whole history before the step is replaced by an equivalent one: an instant influx at the initial moment plus the constant rate of the current step from the very beginning (see the advanced section). The influence function matches the aquifer: a closed one with its own raDr_{aD}, or an infinite one.

Simplified water influx models: Fetkovich and Carter–Tracy →
Pressure at the contact as a convolution of the influx history
U (pi−p(tD))=∫0tD dWe dλ  pD′(tD−λ) dλU\,\big(p_i - p(t_D)\big) = \int_0^{t_D}\frac{\dd W_e}{\dd\lambda}\;p_D'(t_D - \lambda)\dd\lambda

m³ (both sides)

Holds when: The problem is linear. Before the start of development the pressure in the aquifer equals pip_i everywhere. An instant influx of volume VV at the moment λ\lambda gives the response V pD′(tD−λ)V\,p_D'(t_D - \lambda), the derivative of the response to a constant rate.

Simplified water influx models: Fetkovich and Carter–Tracy →
Pressure decline from the equivalent influx history
U Δpn=an pD(tDn)+(We,n−1−an tD,n−1) pD′(tDn)U\,\Delta p_n = a_n\,p_D(t_{Dn}) + \big(W_{e,n-1} - a_n\,t_{D,n-1}\big)\,p_D'(t_{Dn})

m³ (both sides)

Holds when: The equivalent history instead of the real one; the influence function is for the given aquifer.

Simplified water influx models: Fetkovich and Carter–Tracy →
The limit of the Carter–Tracy formula for very short steps
(pD−tD pD′)  dWe dtD+pD′ We=U (pi−p(tD))\big(p_D - t_D\,p_D'\big)\,\frac{\dd W_e}{\dd t_D} + p_D'\,W_e = U\,\big(p_i - p(t_D)\big)

m³ (both sides); time is dimensionless

Holds when: The steps tend to zero. The aquifer and the influence function are the same as in the formula of the method.

Simplified water influx models: Fetkovich and Carter–Tracy →
Dimensionless resistance of a ring-shaped aquifer in the pseudo-steady-state regime
b(raD)=4 raD4ln⁡raD−3 raD4+4 raD2−14 (raD2−1)2  →  raD≫1    ln⁡raD−34,J=2πfkhμw b(raD)b(r_{aD}) = \frac{4\,r_{aD}^4\ln r_{aD} - 3\,r_{aD}^4 + 4\,r_{aD}^2 - 1}{4\,\big(r_{aD}^2 - 1\big)^2}\;\xrightarrow{\;r_{aD}\gg1\;}\;\ln r_{aD} - \frac34,\qquad J = \frac{2\pi f k h}{\mu_w\,b(r_{aD})}

dimensionless; JJ in SI (in metric oilfield units 2π→0.5362\pi \to 0.536)

Holds when: A uniform ring. The outer boundary is impermeable. The flow has settled into the pseudo-steady-state regime.

Simplified water influx models: Fetkovich and Carter–Tracy →
The Fetkovich step when the pressure at the contact changes linearly
ΔWe,n=ctW[(pˉa,n−1−pn−1)(1−e−Δtn/τ)+(pn−1−pn)(1−τΔtn(1−e−Δtn/τ))]\Delta W_{e,n} = c_tW\left[\big(\bar p_{a,n-1} - p_{n-1}\big)\big(1 - \e^{-\Delta t_n/\tau}\big) + \big(p_{n-1} - p_n\big)\left(1 - \frac{\tau}{\Delta t_n}\big(1 - \e^{-\Delta t_n/\tau}\big)\right)\right]

m³, MPa, days

Holds when: The pseudo-steady-state aquifer model. The pressure at the contact changes along a straight line inside the step.

Simplified water influx models: Fetkovich and Carter–Tracy →

Symbols

SymbolMeaning and context
Δp\Delta p
pˉa\bar p_a
ana_n
pˉa,n−1\bar p_{a,n-1}
pˉa,n−1, pˉa,n\bar p_{a,n-1},\ \bar p_{a,n}
bb
ctc_t
ctWc_tW
cw, cfc_w,\ c_f
  • compressibility of the water and of the pore volume of the aquifer (the fraction by which a volume changes for a unit change in pressure), 1/MPaAquifer models: water influx
cw+cfc_w + c_f
WˉD(s)\bar W_D(s)
ΔWe,n\Delta W_{e,n}
 dWe/ dt\dd W_e/\dd t
ff
Fn, Et,nF_n,\ E_{t,n}
f, k, h, μwf,\ k,\ h,\ \mu_w
f, ϕ, ct, h, rof,\ \phi,\ c_t,\ h,\ r_o
  • fraction of the perimeter, porosity, compressibility cw+cfc_w + c_f and thickness of the aquifer, radius of the accumulation, —, —, 1/MPa, m, mAquifer models: water influx
F/EtF/E_t
h, ϕh,\ \phi
Iν, KνI_\nu,\ K_\nu
Δpj\Delta p_j
JJ
JaJ_a
JlinJ_{\text{lin}}
JradJ_{\text{rad}}
k, hk,\ h
k, h, μwk,\ h,\ \mu_w
k, h, μw, fk,\ h,\ \mu_w,\ f
k, μwk,\ \mu_w
k, ϕk,\ \phi
MM
pˉn\bar p_n
Δtn\Delta t_n
Δpn\Delta p_n
NN
pp
PP
p0=pip_0 = p_i
pD, pD′p_D,\ p_D'
pD′p_D'
pip_i
pi−pp_i - p
pi−p(tD)p_i - p(t_D)
pi, pp_i,\ p
pjp_j
pkmeas, pkcalcp_k^{\text{meas}},\ p_k^{\text{calc}}
pn−1, pnp_{n-1},\ p_n
po(t)p_o(t)
P, rDP,\ r_D
qDq_D
qwq_w
rar_a
raDr_{aD}
rDr_D
ror_o
ro, rar_o,\ r_a
RMS\text{RMS}
ss
tt
tat_a
tDt_D
tDn−tDjt_{Dn} - t_{Dj}
tDn, tD,n−1t_{Dn},\ t_{D,n-1}
UU
ULU_L
UnU_n
U, WDU,\ W_D
VpV_p
WW
WDW_D
WD, raD, tDW_D,\ r_{aD},\ t_D
WD=∫0tDqD dtDW_D = \int_0^{t_D}q_D\dd t_D
WeW_e
We,max⁡W_{e,\max}
We,n, We,n−1W_{e,n},\ W_{e,n-1}
We,n−1W_{e,n-1}
We,n−1−an tD,n−1W_{e,n-1} - a_n\,t_{D,n-1}
We(λ)W_e(\lambda)
We(tn)W_e(t_n)
WeiW_{ei}
wnw_n
w, h, Lw,\ h,\ L
W, cw+cfW,\ c_w + c_f
W, cw+cf, ΔpW,\ c_w + c_f,\ \Delta p
w/Etw/E_t
η\eta
μw\mu_w
σ\sigma
τ\tau
τ, Δtn\tau,\ \Delta t_n
ϕ, ct, h\phi,\ c_t,\ h
ϕ, ct, ro\phi,\ c_t,\ r_o