Skip to content
Brent$117,6Газ TTF$1 192USD453,87 ₸
MunaiHub

Формулы: Приток законтурной воды

Формулы уроков темы «Приток законтурной воды» раздела «Физика пласта и разработка» с условиями применимости и единицами, а также обозначения: что означает каждый символ в этих уроках.

Упругий запас водоносного горизонта
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

согласованные единицы: м³, МПа и 1/МПа

Справедливо, если: горизонт насыщен водой, сжимаемости постоянны; радиальная геометрия — сектор кольца с углом θ\theta; для полосообразного горизонта шириной ww и длиной LL объём равен wLhϕwLh\phi; снижение давления успело дойти до всех частей горизонта.

Модели притока законтурной воды →
Приток из малого водоносного горизонта
We=(cw+cf) W (pi−p)W_e = (c_w + c_f)\,W\,(p_i - p)

м³, МПа, 1/МПа (или баррели, psi и 1/psi)

Справедливо, если: давление во всём горизонте равно давлению на контакте: горизонт выравнивается быстрее, чем меняется давление в залежи; питания извне нет; давление на контакте принимают равным среднему давлению залежи.

Модели притока законтурной воды →
Безразмерное время водоносного горизонта и время выхода на его границу
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}

СИ

Справедливо, если: однородный горизонт постоянной толщины, радиальная геометрия; tat_a — оценка порядка величины: заметное снижение давления у границы начинается раньше, а выравнивание всего горизонта наступает в несколько раз позже.

Модели притока законтурной воды →
Стационарная модель притока воды
 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

м³ в пластовых условиях, сутки, МПа; в промысловых единицах — баррели в сутки на psi

Справедливо, если: давление на внешней границе горизонта постоянно и равно начальному; течение в горизонте установившееся: переходный период короток по сравнению с историей разработки; давление на контакте принимают равным среднему давлению залежи.

Модели притока законтурной воды →
Коэффициент продуктивности радиального водоносного горизонта (установившееся течение)
Ja=2π f k hμw ln⁡(ra/ro)J_a = \frac{2\pi\,f\,k\,h}{\mu_w\,\ln(r_a/r_o)}

СИ

Справедливо, если: радиальный горизонт постоянной толщины и проницаемости; давление на радиусе rar_a постоянно; течение установившееся и однофазное.

Модели притока законтурной воды →
Полустационарная модель водоносного горизонта
 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)

м³ в пластовых условиях, сутки, МПа

Справедливо, если: горизонт замкнут; после переходного периода давление в нём падает всюду с одной скоростью; давление на контакте принимают равным среднему давлению залежи.

Модели притока законтурной воды →
Приток при постоянном снижении давления на контакте (нестационарная модель)
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

м³, МПа, 1/МПа, м

Справедливо, если: однородный радиальный горизонт; давление на контакте снизилось на Δp\Delta p скачком и дальше не меняется; для меняющегося давления историю разбивают на ступени и складывают отклики (Принцип суперпозиции).

Модели притока законтурной воды →
Отклик полустационарного горизонта на постоянное снижение давления
 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}

м³, сутки, МПа

Справедливо, если: снижение давления на контакте постоянно; при t=0t = 0 горизонт не тронут: We=0W_e = 0.

Модели притока законтурной воды →
Задача ван Эвердингена–Херста в безразмерном виде
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

безразмерная форма

Справедливо, если: однородный горизонт постоянной толщины, насыщенный водой; сжимаемости малы и постоянны, поэтому уравнение линейно; течение радиальное, без перетока через кровлю и подошву; внешняя граница непроницаема; для бесконечного горизонта условие на ней заменяют условием P→0P \to 0 при rD→∞r_D \to \infty.

Нестационарный приток воды: метод ван Эвердингена–Херста →
Приток при постоянном снижении давления на контакте
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 и UU — в одной системе объёма и давления (м³ и МПа); tDt_D — в СИ

Справедливо, если: задача выше: однородный радиальный горизонт, давление на контакте снизилось скачком и дальше не меняется; давление на контакте принимают равным среднему давлению залежи — так же, как в материальном балансе; горизонт насыщен водой, ct=cw+cfc_t = c_w + c_f, сжимаемости постоянны.

Нестационарный приток воды: метод ван Эвердингена–Херста →
Приток как расширение воронки снижения давления
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

м³, МПа, 1/МПа

Справедливо, если: в горизонте одна вода; ctc_t постоянна; питания извне нет.

Нестационарный приток воды: метод ван Эвердингена–Херста →
Приток при ступенчатом снижении давления на контакте
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)

м³, МПа

Справедливо, если: условия линейности: свойства воды и горизонта не зависят от давления; давление на контакте между моментами ступеней постоянно; насколько хорошо ступени заменяют плавную историю, зависит от выбора перепадов и шага (следующий раздел); все ступени отсчитываются от одного начального состояния: до разработки давление в горизонте всюду равно pip_i.

Нестационарный приток воды: метод ван Эвердингена–Херста →
Перепады ступеней по средним давлениям интервалов
Δ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}

МПа

Справедливо, если: при разных интервалах между замерами формула та же — ступенчатая кривая по-прежнему идёт по средним давлениям интервалов, — но самые длинные интервалы вносят наибольшую ошибку; давление на интервале меняется плавно, без резких скачков (пуск или остановка множества скважин).

Нестационарный приток воды: метод ван Эвердингена–Херста →
Темп притока и безразмерный дебит
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}

СИ

Справедливо, если: постоянное снижение давления Δp\Delta p на контакте с момента t=0t = 0.

Нестационарный приток воды: метод ван Эвердингена–Херста →
От темпа притока к накопленному притоку
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

СИ

Справедливо, если: те же, что для задачи в безразмерном виде.

Нестационарный приток воды: метод ван Эвердингена–Херста →
Изображение безразмерного притока
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)}

безразмерные

Справедливо, если: замкнутая внешняя граница на raDr_{aD}; справа — бесконечный горизонт.

Нестационарный приток воды: метод ван Эвердингена–Херста →
Приток из линейного водоносного горизонта
We=UL Δp WD(tD),UL=w h L ϕ ct,tD=k tϕ μwct L2;WD=2tD/π  (бесконечная полоса)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{(бесконечная полоса)}

м³, МПа, 1/МПа, м; tDt_D — СИ (метрическая форма — с множителем 3,553⋅10−33{,}553\cdot10^{-3} при времени в часах)

Справедливо, если: однородная полоса постоянного сечения; течение перпендикулярно контакту. для бесконечной полосы LL сокращается: We=2whϕctΔpηt/πW_e = 2wh\phi c_t\Delta p\sqrt{\eta t/\pi}.

Нестационарный приток воды: метод ван Эвердингена–Херста →
Приток модели нестационарного горизонта
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}

м³, МПа; tDt_D — в СИ или с множителем 3,553⋅10−33{,}553\cdot10^{-3} при kk в мД, tt в часах, μw\mu_w в мПа·с, ctc_t в 1/МПа

Справедливо, если: однородный радиальный горизонт, давление на контакте равно среднему давлению залежи; функцию ww считают с короткими ступенями или точно: ошибка ступенчатого представления давления иначе переходит в подобранные параметры.

Уточнение параметров водоносного горизонта по истории разработки →
Кажущаяся постоянная водоносного горизонта
Un=We,nwn=Fn−NEt,nwnU_n = \frac{W_{e,n}}{w_n} = \frac{F_n - N E_{t,n}}{w_n}

м³/МПа

Справедливо, если: запасы известны независимо; газовой шапки нет или её размер известен (EtE_t включает mEgmE_g); ранние даты, где снижение давления сравнимо с погрешностью его замера, получают меньший вес.

Уточнение параметров водоносного горизонта по истории разработки →
Прямая Хавлены–Одеха с нестационарной функцией притока
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}

м³, МПа, м³/МПа

Справедливо, если: газовой шапки нет или её размер mm известен; ранние точки, где EtE_t мало, ненадёжны и не должны определять наклон.

Уточнение параметров водоносного горизонта по истории разработки →
Мера согласия с историей давления
M=1n∑k=1n(pkзамер−pkрасчσ)2,M≤1 ⟺ СКО=1n∑k(pkзамер−pkрасч)2≤σM = \frac1n\sum_{k=1}^{n}\left(\frac{p_k^{\text{замер}} - p_k^{\text{расч}}}{\sigma}\right)^{2},\qquad M \le 1\ \Longleftrightarrow\ \text{СКО} = \sqrt{\frac1n\sum_{k}\big(p_k^{\text{замер}} - p_k^{\text{расч}}\big)^2} \le \sigma

безразмерная MM; давления в МПа

Справедливо, если: погрешности замеров независимы и одинаковы по величине; каждый вариант горизонта подобран так, чтобы MM было наименьшим при его raDr_{aD} и проницаемости, — постоянная UU подбирается при этом автоматически.

Уточнение параметров водоносного горизонта по истории разработки →
Полустационарная модель горизонта в обозначениях метода Фетковича
 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

м³, сутки, МПа

Справедливо, если: горизонт замкнут и после переходного периода истощается так, что давление в нём падает всюду с одной скоростью (псевдоустановившийся режим); сжимаемость и проницаемость постоянны; давление на контакте равно среднему давлению залежи.

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Шаг метода Фетковича
Δ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}

м³, сутки, МПа

Справедливо, если: внутри шага давление на контакте заменено его средним значением; полустационарная модель горизонта (уравнения выше).

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Коэффициент продуктивности водоносного горизонта для метода Фетковича
Jрад=2πf k hμw(ln⁡raD−34),Jлин=3 k w hμw LJ_{\text{рад}} = \frac{2\pi f\,k\,h}{\mu_w\left(\ln r_{aD} - \dfrac34\right)},\qquad J_{\text{лин}} = \frac{3\,k\,w\,h}{\mu_w\,L}

СИ

Справедливо, если: псевдоустановившийся режим: напор — разность среднего давления горизонта и давления на контакте; формула для радиального горизонта верна при raD≫1r_{aD} \gg 1; при малых raDr_{aD} — точное выражение для кольца (дополнительный раздел).

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Функция влияния радиального водоносного горизонта
pD(tD)=2πfkh (pi−po(t))qw μw;pD≈2tD/π  (tD≪1),pD≈12(ln⁡tD+0,809)  (бесконечный),pD≈2tDraD2−1+ln⁡raD−34  (замкнутый, поздно)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{(бесконечный)},\qquad p_D \approx \frac{2t_D}{r_{aD}^2 - 1} + \ln r_{aD} - \frac34\ \ \text{(замкнутый, поздно)}

СИ; безразмерная функция

Справедливо, если: однородный радиальный горизонт; темп притока постоянен с момента t=0t = 0; скина на контакте нет; логарифмическая форма — при tDt_D больше примерно 25, пока горизонт ведёт себя как бесконечный; последняя — после выхода на замкнутую границу, при raD≫1r_{aD} \gg 1 (для малых raDr_{aD} вместо ln⁡raD−3/4\ln r_{aD} - 3/4 — точное выражение для кольца).

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Шаг метода Картера–Трейси
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

м³, МПа; время — безразмерное

Справедливо, если: темп притока внутри шага постоянен; вся история до шага заменена эквивалентной: мгновенный приток в начальный момент плюс постоянный темп текущего шага с самого начала (дополнительный раздел); функция влияния отвечает горизонту: замкнутому с его raDr_{aD} или бесконечному.

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Давление на контакте как свёртка истории притока
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

м³ (обе части)

Справедливо, если: линейная задача, до начала разработки давление в горизонте всюду равно pip_i; мгновенный приток объёма VV в момент λ\lambda даёт отклик V pD′(tD−λ)V\,p_D'(t_D - \lambda) — производную отклика на постоянный темп.

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Снижение давления от эквивалентной истории притока
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})

м³ (обе части)

Справедливо, если: эквивалентная история вместо настоящей; функция влияния — для данного горизонта.

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Предел формулы Картера–Трейси при очень коротких шагах
(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)

м³ (обе части); время безразмерное

Справедливо, если: шаги стремятся к нулю; горизонт и функция влияния те же, что в формуле метода.

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Безразмерное сопротивление кольцевого горизонта в псевдоустановившемся режиме
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})}

безразмерная; JJ — в СИ (в метрических промысловых единицах 2π→0,5362\pi \to 0{,}536)

Справедливо, если: однородное кольцо; внешняя граница непроницаема; течение установилось в псевдоустановившемся режиме.

Упрощённые модели притока воды: Фетковича и Картера–Трейси →
Шаг метода Фетковича при линейном изменении давления на контакте
Δ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]

м³, МПа, сутки

Справедливо, если: полустационарная модель горизонта; давление на контакте внутри шага меняется по прямой.

Упрощённые модели притока воды: Фетковича и Картера–Трейси →

Обозначения

СимволЗначение и где встречается
Δp\Delta p
СКО\text{СКО}
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
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
F/EtF/E_t
h, ϕh,\ \phi
Iν, KνI_\nu,\ K_\nu
Δpj\Delta p_j
JJ
JлинJ_{\text{лин}}
JрадJ_{\text{рад}}
JaJ_a
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
pkзамер, pkрасчp_k^{\text{замер}},\ p_k^{\text{расч}}
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
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