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

Кітапхана әзірге орыс тілінде. Ағылшын және қазақ тілдеріндегі нұсқалары кейінірек қосылады.

Формулалар: Повреждение пласта и интенсификация

«Добыча нефти и газа» бөліміндегі «Повреждение пласта и интенсификация» тақырыбы сабақтарының формулалары қолданылу шарттарымен және бірліктерімен, сондай-ақ белгілердің осы сабақтардағы мағынасы.

Скин зоны пониженной проницаемости (формула Хокинса)
Sd=(kkd−1)ln⁡rdrwS_d = \left(\frac{k}{k_d} - 1\right)\ln\frac{r_d}{r_w}

любые согласованные: входят только отношения k/kdk/k_d и rd/rwr_d/r_w; логарифм натуральный

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

Повреждение пласта и скин повреждения →
Скин повреждения для нескольких концентрических зон
Sd=∑i=1n(kki−1)ln⁡riri−1,r0=rwS_d = \sum_{i=1}^{n} \left(\frac{k}{k_i} - 1\right)\ln\frac{r_i}{r_{i-1}},\qquad r_0 = r_w

любые согласованные

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

Повреждение пласта и скин повреждения →
Скин повреждения при непрерывном профиле проницаемости
Sd=∫rwrd(kk(r)−1) drrS_d = \int_{r_w}^{r_d} \left(\frac{k}{k(r)} - 1\right)\frac{\dd r}{r}

любые согласованные

Мына жағдайда орындалады: радиальный установившийся поток в призабойной зоне, однофазное течение по закону Дарси.

Повреждение пласта и скин повреждения →
Инерционная добавка к перепаду давления в призабойной зоне
Δpин=β ρ q24π2h2(1rw−1r)\Delta p_{\text{ин}} = \frac{\beta\,\rho\,q^2}{4\pi^2 h^2}\left(\frac{1}{r_w} - \frac{1}{r}\right)

СИ

Мына жағдайда орындалады: однородная по β\beta зона, радиальный поток по всей толщине; в повреждённой зоне вместо β\beta подставляют её собственный коэффициент βd\beta_d.

Повреждение пласта и скин повреждения →
Коэффициент восстановления образца из нескольких участков
LR=∑iLiRi,Ri=kiпослеkiдо\frac{L}{R} = \sum_i \frac{L_i}{R_i},\qquad R_i = \frac{k_i^{\text{после}}}{k_i^{\text{до}}}

любые согласованные: длины в одних единицах

Мына жағдайда орындалады: поток одномерный, участки работают последовательно; исходная проницаемость одинакова по всей длине образца; иначе участки взвешивают не по длине LiL_i, а по Li/kiдоL_i/k_i^{\text{до}}; проницаемость до и после измерена одним флюидом при той же насыщенности.

Источники повреждения пласта →
Критический расход воды на метр интервала
qкр=uкр Aпр,Aпр=2πrw  или  n π dк lкq_{\text{кр}} = u_{\text{кр}}\,A_{\text{пр}},\qquad A_{\text{пр}} = 2\pi r_w\ \ \text{или}\ \ n\,\pi\,d_{\text{к}}\,l_{\text{к}}

СИ; для расхода в м³/сут умножают на 86 400

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

Источники повреждения пласта →
Индекс насыщения (на примере сульфата бария)
SI=lg⁡γBa [Ba2+]  γSO4 [SO42−]Ksp\mathrm{SI} = \lg\frac{\gamma_{\text{Ba}}\,[\mathrm{Ba^{2+}}]\;\gamma_{\text{SO}_4}\,[\mathrm{SO_4^{2-}}]}{K_{sp}}

концентрации в моль/л; логарифм десятичный

Мына жағдайда орындалады: SI > 0 — вода пересыщена, соль может выпасть; SI < 0 — недонасыщена; SI = 0 — равновесие; в рассолах коэффициенты активности много меньше единицы; их и KspK_{sp} при пластовых температуре и давлении рассчитывают по моделям рассолов в специальных программах; индекс термодинамический: он говорит, может ли соль выпасть, но не как быстро.

Источники повреждения пласта →
Сколько сульфата бария может выпасть из смеси
x=(a+b)−(a−b)2+4Ksp′2x = \frac{(a + b) - \sqrt{(a - b)^2 + 4K'_{sp}}}{2}

моль/л; масса осадка — xx, умноженное на молярную массу сульфата бария (233,4 г/моль)

Мына жағдайда орындалады: смесь дошла до равновесия; другие ионы и осадки не учтены; коэффициенты активности при выпадении не меняются.

Источники повреждения пласта →
График Холла
∑j(pw,j−pˉ)Δtj=μB(ln⁡rerw+S)2πkh W\sum_j \left(p_{w,j} - \bar p\right)\Delta t_j = \frac{\mu B\left(\ln\dfrac{r_e}{r_w} + S\right)}{2\pi k h}\,W

СИ; в метрических промысловых единицах 2π2\pi заменяют на 0,536, и наклон получается в МПа·сут/м³ при kk в мД, hh в м, μ\mu в мПа·с

Мына жағдайда орындалады: закачка близка к установившейся: давление на расстоянии rer_e равно пластовому pˉ\bar p и почти не меняется; подвижности воды и вытесняемой нефти близки; часто вместо забойного давления берут устьевое: при постоянном расходе гидростатика и трение дают постоянную поправку; наклон меняется и по другим причинам — рост пластового давления, обводнение зоны вокруг скважины; точнее всего скин измеряет исследование на падение давления.

Источники повреждения пласта →
Кратность увеличения дебита после снижения скина
q2q1=X+S1X+S2,X=ln⁡rerw−34\frac{q_2}{q_1} = \frac{X + S_1}{X + S_2},\qquad X = \ln\frac{r_e}{r_w} - \frac34

безразмерная; логарифм натуральный

Мына жағдайда орындалады: среднее пластовое давление и депрессия до и после обработки одинаковы: подъёмник удерживает прежнее забойное давление при большем дебите; псевдоустановившийся однофазный приток в круговой области дренирования со скважиной в центре; обработка меняет только скин; проницаемость пласта и его толщина прежние.

Кислотные обработки: принцип, экономика и выбор скважин →
Абсолютный прирост дебита после снижения скина
Δq=2πkh ΔpμB(1X+S2−1X+S1)\Delta q = \frac{2\pi kh\,\Delta p}{\mu B}\left(\frac{1}{X + S_2} - \frac{1}{X + S_1}\right)

СИ (любая согласованная система единиц)

Мына жағдайда орындалады: те же, что для формулы кратности.

Кислотные обработки: принцип, экономика и выбор скважин →
Срок окупаемости при экспоненциальном спаде прироста
tок=−1aln⁡(1−aCΔq0 v)t_{\text{ок}} = -\frac{1}{a}\ln\left(1 - \frac{aC}{\Delta q_0\,v}\right)

любые согласованные: здесь сутки, м³/сут и доллары

Мына жағдайда орындалады: прирост спадает экспоненциально с постоянным темпом aa, чистый доход с кубометра постоянен; дисконтирование не учитывается: при сроках окупаемости до года оно мало меняет ответ; решение существует, только если aC<Δq0vaC < \Delta q_0 v. Иначе весь будущий доход от прироста, Δq0v/a\Delta q_0 v/a, меньше затрат, и обработка не окупится никогда.

Кислотные обработки: принцип, экономика и выбор скважин →
Добавка к добыче до предела рентабельности после обработки
ΔN=qпрD1(1−1F)\Delta N = \frac{q_{\text{пр}}}{D_1}\left(1 - \frac{1}{F}\right)

любые согласованные: дебит и темп падения в одних единицах времени

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

Кислотные обработки: принцип, экономика и выбор скважин →
Реакции соляной кислоты с карбонатами
2 HCl+CaCO3→CaCl2+CO2+H2O4 HCl+CaMg(CO3)2→CaCl2+MgCl2+2 CO2+2 H2O\begin{aligned} 2\,\mathrm{HCl} + \mathrm{CaCO_3} &\rightarrow \mathrm{CaCl_2} + \mathrm{CO_2} + \mathrm{H_2O}\\ 4\,\mathrm{HCl} + \mathrm{CaMg(CO_3)_2} &\rightarrow \mathrm{CaCl_2} + \mathrm{MgCl_2} + 2\,\mathrm{CO_2} + 2\,\mathrm{H_2O}\end{aligned}

коэффициенты — числа молей

Мына жағдайда орындалады: кислота расходуется полностью; реакция с кальцитом очень быстрая, с доломитом медленнее, особенно при низкой температуре; сидерит FeCO3\mathrm{FeCO_3} реагирует так же, как кальцит, но переводит в раствор железо.

Кислотные составы и проектирование обработки →
Реакции плавиковой кислоты
SiO2+4 HF→SiF4+2 H2O,SiF4+2 HF→H2SiF6Al4Si4O10(OH)8+32 HF→4 AlF2++4 SiF62−+18 H2O+4 H+CaCO3+2 HF→CaF2↓+CO2+H2O\begin{aligned} \mathrm{SiO_2} + 4\,\mathrm{HF} &\rightarrow \mathrm{SiF_4} + 2\,\mathrm{H_2O},\qquad \mathrm{SiF_4} + 2\,\mathrm{HF} \rightarrow \mathrm{H_2SiF_6}\\ \mathrm{Al_4Si_4O_{10}(OH)_8} + 32\,\mathrm{HF} &\rightarrow 4\,\mathrm{AlF_2^{+}} + 4\,\mathrm{SiF_6^{2-}} + 18\,\mathrm{H_2O} + 4\,\mathrm{H^{+}}\\ \mathrm{CaCO_3} + 2\,\mathrm{HF} &\rightarrow \mathrm{CaF_2}\downarrow + \mathrm{CO_2} + \mathrm{H_2O}\end{aligned}

коэффициенты — числа молей

Мына жағдайда орындалады: кислота в избытке: кремний уходит во фторсиликат-ион, алюминий — в комплексы с двумя фторидами; при нехватке фтора продукты другие, и HF на моль минерала уходит меньше; с глинами и полевыми шпатами реакции быстрые, с кварцем медленные; скорость растёт с температурой; с кальциевыми минералами HF даёт нерастворимый фторид кальция — поэтому их удаляют соляной кислотой до плавиковой.

Кислотные составы и проектирование обработки →
Растворяющая способность кислоты
β=νмMмνкMк,X=C β ρкρм\beta = \frac{\nu_{\text{м}} M_{\text{м}}}{\nu_{\text{к}} M_{\text{к}}},\qquad X = \frac{C\,\beta\,\rho_{\text{к}}}{\rho_{\text{м}}}

отношения масс и объёмов; единицы плотностей и молярных масс одинаковы в числителе и знаменателе

Мына жағдайда орындалады: кислота расходуется полностью и только на этот минерал; для HF коэффициент νк\nu_{\text{к}} зависит от того, какие фториды образуются: здесь кремний переходит во фторсиликат-ион; в таблицах, где продуктом считают SiF4\mathrm{SiF_4}, значения в 1,3–1,5 раза больше; плотность раствора берут при температуре его приготовления.

Кислотные составы и проектирование обработки →
Объём кислоты на метр интервала для фронта растворения на радиусе r
Vh=π (r2−rw2)[ϕ+(1−ϕ) xX]\frac{V}{h} = \pi\,(r^2 - r_w^2)\left[\phi + (1 - \phi)\,\frac{x}{X}\right]

СИ (м, м³/м)

Мына жағдайда орындалады: кислота поступает в пласт равномерно по толщине и по окружности; минерал растворяется за резким фронтом полностью (быстрая реакция); медленная реакция HF с кварцем, перемешивание и рост пористости не учтены — объём занижен, поэтому к расчёту добавляют запас; для глинокислоты xx — доля быстро реагирующих минералов, доступных кислоте: глин на стенках пор, мелочи, частиц бурового раствора.

Кислотные составы и проектирование обработки →
Скин кольца с осадком из отработанной кислоты
Sп=(kkп−1)ln⁡r2r1=12(kkп−1)ln⁡(1+Vπϕh r12)S_{\text{п}} = \left(\frac{k}{k_{\text{п}}} - 1\right)\ln\frac{r_2}{r_1} = \frac12\left(\frac{k}{k_{\text{п}}} - 1\right)\ln\left(1 + \frac{V}{\pi\phi h\,r_1^2}\right)

любые согласованные; логарифм натуральный

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

Кислотные составы и проектирование обработки →
Предельный расход матричной закачки
qmax⁡=2πkh (pf−Δpзап−pˉ)μ(ln⁡rerw+S)q_{\max} = \frac{2\pi k h\,(p_f - \Delta p_{\text{зап}} - \bar p)}{\mu\left(\ln\dfrac{r_e}{r_w} + S\right)}

СИ; ln⁡\ln — натуральный логарифм

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

Кислотные составы и проектирование обработки →
Давление на устье при закачке
pу=pw−ρgH+Δpтрp_{\text{у}} = p_w - \rho g H + \Delta p_{\text{тр}}

СИ; в метрических единицах ρgH=0,00981 ρH\rho g H = 0{,}00981\,\rho H (МПа при ρ\rho в г/см³ и HH в м)

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

Кислотные составы и проектирование обработки →
Кажущийся скин во время закачки
Sкаж=2πkh (pw−pˉ)q μ−ln⁡rerwS_{\text{каж}} = \frac{2\pi k h\,(p_w - \bar p)}{q\,\mu} - \ln\frac{r_e}{r_w}

СИ; в метрических промысловых единицах 2π2\pi заменяют на 0,536 (расход в м³/сут)

Мына жағдайда орындалады: поток установившийся; неустановившиеся эффекты в начале закачки и после смены расхода искажают значение; в кажущийся скин входят и псевдоскины, и вязкий вал кислоты, и действие отклонителя.

Кислотные составы и проектирование обработки →
Число кислотной ёмкости и объём этапа
Nк=ϕ X(1−ϕ) x,Vh=π (r2−rw2) ϕ(1+1Nк)N_{\text{к}} = \frac{\phi\,X}{(1 - \phi)\,x},\qquad \frac{V}{h} = \pi\,(r^2 - r_w^2)\,\phi\left(1 + \frac{1}{N_{\text{к}}}\right)

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

Мына жағдайда орындалады: те же, что для расчёта объёма по резкому фронту.

Кислотные составы и проектирование обработки →
Скин зоны вязкого отклонителя
Sв=(μвμ−1)ln⁡rвrwS_{\text{в}} = \left(\frac{\mu_{\text{в}}}{\mu} - 1\right)\ln\frac{r_{\text{в}}}{r_w}

любые согласованные; логарифм натуральный

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

Кислотные составы и проектирование обработки →
Число Дамкёлера для червоточины
Da=π d l κqwhDa = \frac{\pi\,d\,l\,\kappa}{q_{wh}}

СИ (безразмерное)

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

Кислотные обработки карбонатов и особых скважин →
Скорость роста червоточин и кривая эффективности кислоты
vwh=vopt1/3PVbt,opt vi2/3[1−e−4vi2/vopt2]2,PVbt=vivwh=PVbt,opt (vi/vopt)1/3[1−e−4vi2/vopt2]2v_{wh} = \frac{v_{opt}^{1/3}}{PV_{bt,opt}}\,v_i^{2/3}\left[1 - e^{-4v_i^2/v_{opt}^2}\right]^2,\qquad PV_{bt} = \frac{v_i}{v_{wh}} = PV_{bt,opt}\,\frac{(v_i/v_{opt})^{1/3}}{\left[1 - e^{-4v_i^2/v_{opt}^2}\right]^2}

любые согласованные единицы скорости

Мына жағдайда орындалады: постоянные voptv_{opt} и PVbt,optPV_{bt,opt} берут из опытов на керне той же породы с той же кислотой при пластовой температуре; формула сглаживает минимум: при vi=voptv_i = v_{opt} она даёт PVbt≈1,04 PVbt,optPV_{bt} \approx 1{,}04\,PV_{bt,opt}; при больших скоростях PVbtPV_{bt} растёт как vi1/3v_i^{1/3}, при малых — очень круто; на образцах большего размера оптимальная скорость и PVbtPV_{bt} меньше, чем на стандартных, поэтому лабораторные значения к пласту переносят с поправкой или калибруют по результатам обработок.

Кислотные обработки карбонатов и особых скважин →
Скорость в порах на фронте червоточин при радиальной закачке
vi=q2πr h ϕ,ropt=q2πhϕ voptv_i = \frac{q}{2\pi r\,h\,\phi},\qquad r_{opt} = \frac{q}{2\pi h\phi\,v_{opt}}

СИ

Мына жағдайда орындалады: кислота входит в интервал равномерно по толщине и по окружности; фронт червоточин считается цилиндрическим; сами червоточины занимают ничтожную долю порового объёма.

Кислотные обработки карбонатов и особых скважин →
Объём кислоты на метр для червоточин до радиуса r_wh
Vh=∫rwrwh2πr ϕ PVbt(vi(r)) dr  ≈  π ϕ PVbt (rwh2−rw2)\frac{V}{h} = \int_{r_w}^{r_{wh}} 2\pi r\,\phi\,PV_{bt}\big(v_i(r)\big)\,dr \;\approx\; \pi\,\phi\,PV_{bt}\,\left(r_{wh}^2 - r_w^2\right)

СИ (м, м³/м)

Мына жағдайда орындалады: фронт червоточин продвигается со скоростью vwh(vi)v_{wh}(v_i), как в линейном опыте при той же скорости в порах; расход постоянен; вторая форма — «объёмная модель» с одним значением PVbtPV_{bt}: годится, если скорость на всём пути близка к оптимальной; иначе интеграл считают по кольцам; порода однородна; естественные трещины и каверны принимают кислоту быстрее и меняют картину.

Кислотные обработки карбонатов и особых скважин →
Доля расхода прослоя при закачке с общим забойным давлением
qjq=kjhj/(ln⁡(re/rw)+Sj)∑ikihi/(ln⁡(re/rw)+Si)\frac{q_j}{q} = \frac{k_jh_j\big/\left(\ln(r_e/r_w) + S_j\right)}{\sum_i k_ih_i\big/\left(\ln(r_e/r_w) + S_i\right)}

доли; единицы kk и hh одинаковы для всех прослоев

Мына жағдайда орындалады: установившаяся закачка, как в формуле предельного расхода; пластовое давление и вязкость жидкости в прослоях одинаковы, перетока между прослоями у ствола нет; вал отклонителя добавляет к скину прослоя (μв/μ−1)ln⁡(rв/rw)(\mu_{\text{в}}/\mu - 1)\ln(r_{\text{в}}/r_w), где μв\mu_{\text{в}} и rвr_{\text{в}} — его кажущаяся вязкость и радиус (скин вязкого отклонителя).

Кислотные обработки карбонатов и особых скважин →
Идеальная ширина травления кислотной трещины
wi=X V(1−ϕ)⋅2xfhfw_i = \frac{X\,V}{(1 - \phi)\cdot 2x_f h_f}

СИ

Мына жағдайда орындалады: вся кислота прореагировала на стенках в пределах трещины — без утечек в пласт и без расхода у ствола; на деле ширина меньше; травление равномерно; проводимость же создаёт как раз неравномерность, поэтому wiw_i — мера объёма работы кислоты, а не ширина канала после закрытия.

Кислотные обработки карбонатов и особых скважин →
Дополнительная добыча за период против падающей базы
ΔNT=q2(1−e−a2T)a2−q1(1−e−a1T)a1,a2=a1 X+S1X+S2\Delta N_T = \frac{q_2\left(1 - e^{-a_2T}\right)}{a_2} - \frac{q_1\left(1 - e^{-a_1T}\right)}{a_1},\qquad a_2 = a_1\,\frac{X + S_1}{X + S_2}

согласованные единицы: дебит и темп падения в одних единицах времени

Мына жағдайда орындалады: замкнутая область дренирования без поддержки давления, постоянное забойное давление (депрессия убывает вместе с пластовым давлением): падение экспоненциальное, а его темп пропорционален продуктивности; при поддержке давления закачкой прирост держится дольше, и формула его занижает; при обводнении после обработки — завышает.

Кислотные обработки карбонатов и особых скважин →
Скин скважины с червоточинами внутри повреждённого кольца
S={kkdln⁡rdrwh−ln⁡rdrw,rwh<rd−ln⁡rwhrw,rwh≥rdS = \begin{cases} \dfrac{k}{k_d}\ln\dfrac{r_d}{r_{wh}} - \ln\dfrac{r_d}{r_w}, & r_{wh} < r_d\\[2mm] -\ln\dfrac{r_{wh}}{r_w}, & r_{wh} \ge r_d \end{cases}

безразмерная; логарифм натуральный

Мына жағдайда орындалады: зона червоточин проводит бесконечно хорошо, давление в ней равно забойному: скин получается предельным, на деле он несколько выше; повреждение равномерно в кольце rw…rdr_w \ldots r_d; при rwh=rwr_{wh} = r_w формула переходит в формулу Хокинса (Скин-фактор).

Кислотные обработки карбонатов и особых скважин →
Безразмерная проводимость трещины
FCD=kf wk xfF_{CD} = \frac{k_f\,w}{k\,x_f}

любые согласованные: kfk_f и kk в одних единицах, ww и xfx_f — тоже

Мына жағдайда орындалады: вертикальная трещина на всю толщину пласта с двумя одинаковыми крыльями; проводимость трещины и проницаемость пласта одинаковы по длине; однофазное течение по закону Дарси; инерционные потери и многофазный поток в трещине уменьшают её эффективную проводимость.

Гидроразрыв пласта: назначение и выбор скважин →
Эквивалентный скин трещины высокой проводимости
Sf=−ln⁡rw′rw,rw′≈xf2S_f = -\ln\frac{r_w'}{r_w},\qquad r_w' \approx \frac{x_f}{2}

любые единицы длины, одинаковые для всех радиусов

Мына жағдайда орындалады: трещина высокой проводимости: FCDF_{CD} от 10 и выше; при меньшей проводимости приведённый радиус меньше (следующий раздел); трещина мала по сравнению с областью дренирования: полудлина заметно меньше rer_e; псевдоустановившийся приток; в первые месяцы, пока течение к трещине линейное, продуктивность выше; трещина обходит повреждение, поэтому скин повреждения в сумму не входит.

Гидроразрыв пласта: назначение и выбор скважин →
Кратность увеличения дебита после гидроразрыва
q2q1=X+S1X+Sf,X=ln⁡rerw−34\frac{q_2}{q_1} = \frac{X + S_1}{X + S_f},\qquad X = \ln\frac{r_e}{r_w} - \frac34

безразмерная; логарифм натуральный

Мына жағдайда орындалады: одинаковые среднее пластовое давление и депрессия; подъёмник и наземное оборудование пропускают возросший дебит; псевдоустановившийся приток в круговой области дренирования.

Гидроразрыв пласта: назначение и выбор скважин →
Приведённый радиус трещины конечной проводимости (аппроксимация)
ln⁡xfrw′=1,65−0,328 u+0,116 u21+0,180 u+0,064 u2+0,005 u3,u=ln⁡FCD\ln\frac{x_f}{r_w'} = \frac{1{,}65 - 0{,}328\,u + 0{,}116\,u^2}{1 + 0{,}180\,u + 0{,}064\,u^2 + 0{,}005\,u^3},\qquad u = \ln F_{CD}

безразмерная; логарифмы натуральные

Мына жағдайда орындалады: вертикальная трещина на всю толщину пласта, скважина в центре области дренирования, трещина мала по сравнению с областью; аппроксимация графика с погрешностью в несколько процентов в практически важном диапазоне — примерно от 0,3 до сотен; при меньшей проводимости она завышает rw′r_w' (при FCD=0,1F_{CD} = 0{,}1 — примерно на 12 %, подробнее в уроке Приток к скважине после гидроразрыва); при больших FCDF_{CD} правая часть приближается к ln⁡2\ln 2, то есть rw′r_w' — к xf/2x_f/2.

Гидроразрыв пласта: назначение и выбор скважин →
Безразмерное время по полудлине трещины
tDxf=k tϕ μ ct xf2t_{Dxf} = \frac{k\,t}{\phi\,\mu\,c_t\,x_f^2}

СИ (любая согласованная система единиц)

Мына жағдайда орындалады: однофазный слабосжимаемый флюид; та же пьезопроводность k/(ϕμct)k/(\phi\mu c_t), что в уравнении пьезопроводности (Уравнение пьезопроводности).

Приток к скважине после гидроразрыва →
Эквивалентный скин трещины через универсальную кривую
Sf=ln⁡xfrw′∣FCD−ln⁡xfrwS_f = \ln\frac{x_f}{r_w'}\bigg|_{F_{CD}} - \ln\frac{x_f}{r_w}

безразмерная; логарифмы натуральные

Мына жағдайда орындалады: псевдорадиальный приток: возмущение давления ушло от трещины на расстояние в несколько её полудлин; трещина мала по сравнению с областью дренирования (коэффициент проникновения меньше примерно 0,3, см. следующий раздел); однородная проводимость по длине; пережатие у ствола и повреждение стенок учитывают отдельными скинами (последний раздел).

Приток к скважине после гидроразрыва →
Приведённый радиус трещины низкой проводимости
rw′≈12γ kfwk=0,2807 kfwkr_w' \approx \frac{1}{2\gamma}\,\frac{k_f w}{k} = 0{,}2807\,\frac{k_f w}{k}

kfw/kk_f w/k имеет размерность длины: при kfwk_f w в мД·м и kk в мД радиус получается в метрах

Мына жағдайда орындалады: FCD≲0,1F_{CD} \lesssim 0{,}1: при FCDF_{CD} = 0,1 точная кривая лежит на 3,3 % ниже асимптоты, при 0,05 — на 1,6 %; псевдорадиальный приток; трещина длиннее участка, который реально работает (длины порядка kfw/kk_f w/k).

Приток к скважине после гидроразрыва →
Коэффициент проникновения трещины
Ix=2xfxeI_x = \frac{2x_f}{x_e}

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

Мына жағдайда орындалады: трещина проходит через центр области параллельно её стороне; при IxI_x = 1 концы трещины достигают границ.

Приток к скважине после гидроразрыва →
Безразмерный коэффициент продуктивности
JD=qBμ2πkh (pˉ−pwf),J=2πkhBμ JDJ_D = \frac{qB\mu}{2\pi kh\,(\bar p - p_{wf})},\qquad J = \frac{2\pi kh}{B\mu}\,J_D

СИ (любая согласованная система единиц)

Мына жағдайда орындалады: псевдоустановившийся режим в замкнутой области (Псевдоустановившийся приток и коэффициент формы); в неустановившемся режиме JDJ_D зависит от времени (раздел о первых дебитах); однофазный приток; у вертикальной скважины в квадрате 1/JD1/J_D равно ln⁡(xe/rw)−1,311+S\ln(x_e/r_w) - 1{,}311 + S, где −1,311-1{,}311 — это 12ln⁡(4/(γCA))\tfrac12\ln\big(4/(\gamma C_A)\big) при CAC_A = 30,88.

Приток к скважине после гидроразрыва →
Предел безразмерного коэффициента продуктивности
JD≤JD,max⁡=6π⋅xeyeJ_D \le J_{D,\max} = \frac{6}{\pi}\cdot\frac{x_e}{y_e}

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

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

Приток к скважине после гидроразрыва →
Билинейный приток к трещине конечной проводимости
pwD=πΓ(5/4)2 tDxf1/4FCD=2,45 tDxf1/4FCDp_{wD} = \frac{\pi}{\Gamma(5/4)\sqrt2}\,\frac{t_{Dxf}^{1/4}}{\sqrt{F_{CD}}} = 2{,}45\,\frac{t_{Dxf}^{1/4}}{\sqrt{F_{CD}}}

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

Мына жағдайда орындалады: ранние времена: при FCD≥3F_{CD} \ge 3 режим длится примерно до tDxf≈0,1/FCD2t_{Dxf} \approx 0{,}1/F_{CD}^2; при FCDF_{CD} от 1,6 до 3 — дольше, чем по этой формуле (при FCDF_{CD} = 1,6 — до tDxf≈0,7t_{Dxf} \approx 0{,}7); при FCD<1F_{CD} < 1 — меньше: работает лишь короткий участок трещины у ствола, и уже при tDxft_{Dxf} порядка FCD2/400F_{CD}^2/400 приток начинает переходить в псевдорадиальный, минуя линейный режим; трещина однородна и не пережата у ствола; влияние ствола закончилось.

Приток к скважине после гидроразрыва →
Линейный приток к трещине высокой проводимости
pwD=π tDxf⟺Δp=qBμ2hxftπkϕμctp_{wD} = \sqrt{\pi\,t_{Dxf}}\quad\Longleftrightarrow\quad \Delta p = \frac{qB\mu}{2hx_f}\sqrt{\frac{t}{\pi k\phi\mu c_t}}

СИ (любая согласованная система единиц)

Мына жағдайда орындалады: FCDF_{CD} больше примерно 300; режим длится до tDxf≈0,016t_{Dxf} \approx 0{,}016, при меньшей проводимости его вытесняет билинейный приток; трещина на всю толщину пласта; у трещины, которая ниже пласта, начальный участок искажён.

Приток к скважине после гидроразрыва →
Псевдорадиальный приток к трещине
pwD=12(ln⁡tDxf+0,809)+ln⁡xfrw′p_{wD} = \frac12\left(\ln t_{Dxf} + 0{,}809\right) + \ln\frac{x_f}{r_w'}

безразмерная; логарифмы натуральные

Мына жағдайда орындалады: производная выходит на полку 0,5 с точностью 7–8 % при tDxf=3t_{Dxf} = 3 и 5 % при tDxf=5t_{Dxf} = 5 при любых FCDF_{CD} от 1 до бесконечности; бесконечный пласт: возмущение ещё не дошло до границ области.

Приток к скважине после гидроразрыва →
Неустановившийся коэффициент продуктивности и кажущаяся кратность прироста
JD(t)=1pwD(t)−2πtDA,nкаж(t)=JD(t)JD,доJ_D(t) = \frac{1}{p_{wD}(t) - 2\pi t_{DA}},\qquad n_{\text{каж}}(t) = \frac{J_D(t)}{J_{D,\text{до}}}

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

Мына жағдайда орындалады: замкнутая однородная область, однофазный приток; до обработки скважина работала в псевдоустановившемся режиме; при постоянном забойном давлении дебит меняется во времени похоже, но не точно так же: эта формула — мера, а не прогноз дебита по дням.

Приток к скважине после гидроразрыва →
Скин пережатия трещины
sch≈πk xs(1kf,chw−1kfw)≈πk xskf,chws_{ch} \approx \pi k\,x_s\left(\frac{1}{k_{f,ch}w} - \frac{1}{k_f w}\right) \approx \frac{\pi k\,x_s}{k_{f,ch}w}

безразмерная; проницаемости и проводимости в согласованных единицах (мД и мД·м)

Мына жағдайда орындалады: весь поток крыла проходит пережатый участок вдоль трещины: на каждое крыло — половина дебита; это верхняя оценка: если kf,chwk_{f,ch}w не намного больше k xsk\,x_s, часть флюида обходит пережатый участок по пласту, и фактический скин меньше (в примере ниже — примерно на 40 %).

Приток к скважине после гидроразрыва →
Скин стенок трещины
sfs=πbs2xf(kks−1)s_{fs} = \frac{\pi b_s}{2x_f}\left(\frac{k}{k_s} - 1\right)

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

Мына жағдайда орындалады: тонкий слой, bs≪xfb_s \ll x_f, одинаковый вдоль всей трещины; поток через него перпендикулярен стенке; в линейном режиме sfss_{fs} прибавляется к pwDp_{wD} точно; в псевдорадиальном итоговый скин растёт примерно на sfss_{fs}, но не меньше: по модели — на 5–20 % больше при sfss_{fs} от 0,5 до 1 и относительно ещё больше при малых sfss_{fs}, потому что приток перераспределяется вдоль трещины.

Приток к скважине после гидроразрыва →
Горизонтальные напряжения пороупругого слоя с тектоническими деформациями
σh=ν1−ν(σv−αpp)+αpp+E1−ν2(εh+ν εH)\sigma_h = \frac{\nu}{1-\nu}\left(\sigma_v - \alpha p_p\right) + \alpha p_p + \frac{E}{1-\nu^2}\left(\varepsilon_h + \nu\,\varepsilon_H\right)

любые одинаковые единицы напряжения и модуля (МПа или psi); деформации безразмерны

Мына жағдайда орындалады: слой линейно упругий и изотропный, широкий по сравнению с толщиной; вертикальное напряжение — главное; деформации отсчитаны от состояния, в котором горизонтальные напряжения создавал только вес пород; εh\varepsilon_h и εH\varepsilon_H одинаковы для всех слоёв разреза, но напряжения от них в каждом слое свои; модули статические — измеренные при медленном нагружении керна или пересчитанные из динамических по зависимости, проверенной на керне (Акустический каротаж); σH\sigma_H получают той же формулой, поменяв местами εh\varepsilon_h и εH\varepsilon_H; при εh=εH=0\varepsilon_h = \varepsilon_H = 0 оба горизонтальных напряжения равны, и при α=1\alpha = 1 получается формула Итона.

Механика горных пород для гидроразрыва →
Изменение горизонтального напряжения при изменении порового давления
Δσh=A Δpp,A=α (1−2ν)1−ν\Delta\sigma_h = A\,\Delta p_p,\qquad A = \frac{\alpha\,(1 - 2\nu)}{1 - \nu}

любые одинаковые единицы давления

Мына жағдайда орындалады: давление меняется во всём широком пласте, и пласт деформируется только по вертикали (одноосная деформация, uniaxial strain); порода упругая; горное давление при разработке не меняется; измеренные на месторождениях отношения Δσh/Δpp\Delta\sigma_h/\Delta p_p разбросаны и могут заметно отличаться от расчётного; его проверяют повторными замерами давления закрытия.

Механика горных пород для гидроразрыва →
Напряжения на стенке вертикальной скважины (решение Кирша)
σθ=σH+σh−2 (σH−σh)cos⁡2θ−pw,σr=pw\sigma_\theta = \sigma_H + \sigma_h - 2\,(\sigma_H - \sigma_h)\cos 2\theta - p_w,\qquad \sigma_r = p_w

любые одинаковые единицы напряжения; сжатие положительно

Мына жағдайда орындалады: порода линейно упругая, однородная и изотропная; скважина вертикальна и вертикальное напряжение — главное; поровое давление у стенки не отличается от пластового: стенку отделяет от скважины фильтрационная корка или давление меняется так быстро, что жидкость не успевает войти в породу (влияние фильтрации — ниже); температура стенки равна пластовой; плоская деформация вдоль оси скважины; с удалением от стенки возмущение быстро затухает, и уже в нескольких радиусах от оси напряжения почти не отличаются от σH\sigma_H и σh\sigma_h (раздел для углублённого изучения).

Механика горных пород для гидроразрыва →
Давление разрыва стенки при проникновении жидкости в породу
pb=3σh−σH+T−A pp2−A,A=α (1−2ν)1−νp_b = \frac{3\sigma_h - \sigma_H + T - A\,p_p}{2 - A},\qquad A = \frac{\alpha\,(1 - 2\nu)}{1 - \nu}

любые одинаковые единицы давления

Мына жағдайда орындалады: те же, что у решения Кирша, но поровое давление на стенке равно давлению в скважине, а изменение порового давления вокруг скважины осесимметрично; разрыв — когда эффективное окружное напряжение σθ−pw\sigma_\theta - p_w достигает −T-T; медленное нагружение: в каждый момент порода успевает прийти в равновесие с давлением в порах; реальное давление разрыва лежит между этим значением и значением без проникновения; какое ближе, зависит от скорости закачки, вязкости жидкости, корки и проницаемости.

Механика горных пород для гидроразрыва →
Забойное давление при закачке и после остановки
pwf=pc+pnet+Δpперф+Δpпз,pISIP=pc+pnetp_{wf} = p_c + p_{net} + \Delta p_{\text{перф}} + \Delta p_{\text{пз}},\qquad p_{ISIP} = p_c + p_{net}

любые одинаковые единицы давления; все давления — на глубине пласта

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

Механика горных пород для гидроразрыва →
Функция G (для трещины, площадь которой росла пропорционально времени)
G(ΔtD)=163π[(1+ΔtD)1,5−ΔtD1,5−1],ΔtD=ΔttpG(\Delta t_D) = \frac{16}{3\pi}\left[(1 + \Delta t_D)^{1{,}5} - \Delta t_D^{1{,}5} - 1\right],\qquad \Delta t_D = \frac{\Delta t}{t_p}

безразмерная; Δt\Delta t и tpt_p в одних единицах

Мына жағдайда орындалады: закачка с постоянным расходом; после остановки площадь трещины не меняется; фильтрация из трещины в пласт подчиняется закону Картера: скорость утечки обратно пропорциональна корню из времени, прошедшего с момента, когда данный участок стенки открылся; жёсткость трещины постоянна: чистое давление пропорционально объёму жидкости в ней; существует и вторая предельная форма функции — для трещины, площадь которой росла как корень из времени; две формы различаются не больше чем примерно на 20 % (тем меньше, чем позже закрытие), а на найденное давление закрытия выбор не влияет.

Механика горных пород для гидроразрыва →
Наибольшая ширина трещины при заданном чистом давлении
wmax⁡=2 pnet hfE′,E′=E1−ν2w_{\max} = \frac{2\,p_{net}\,h_f}{E'},\qquad E' = \frac{E}{1 - \nu^2}

согласованные: pnetp_{net} и E′E' в одних единицах, ww в единицах hfh_f

Мына жағдайда орындалады: трещина много длиннее своей высоты, её вертикальное сечение — эллипс (плоская деформация); чистое давление одинаково по высоте; порода однородна и упруга; это оценка порядка: на самом деле давление падает к кончику трещины, а слои разной жёсткости раскрываются по-разному.

Механика горных пород для гидроразрыва →
Равновесная высота трещины между двумя барьерами
hf=hcos⁡(π pnet2 Δσ),pnet<Δσh_f = \frac{h}{\cos\left(\dfrac{\pi\,p_{net}}{2\,\Delta\sigma}\right)},\qquad p_{net} < \Delta\sigma

pnetp_{net} и Δσ\Delta\sigma в одних единицах; hfh_f в единицах hh; угол в радианах

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

Механика горных пород для гидроразрыва →
Решение Кирша для вертикальной скважины
σr=σH+σh2(1−rw2r2)+σH−σh2(1−4rw2r2+3rw4r4)cos⁡2θ+pwrw2r2σθ=σH+σh2(1+rw2r2)−σH−σh2(1+3rw4r4)cos⁡2θ−pwrw2r2τrθ=−σH−σh2(1+2rw2r2−3rw4r4)sin⁡2θ\begin{aligned} \sigma_r &= \tfrac{\sigma_H+\sigma_h}{2}\left(1-\tfrac{r_w^2}{r^2}\right) + \tfrac{\sigma_H-\sigma_h}{2}\left(1-\tfrac{4r_w^2}{r^2}+\tfrac{3r_w^4}{r^4}\right)\cos 2\theta + p_w\tfrac{r_w^2}{r^2}\\ \sigma_\theta &= \tfrac{\sigma_H+\sigma_h}{2}\left(1+\tfrac{r_w^2}{r^2}\right) - \tfrac{\sigma_H-\sigma_h}{2}\left(1+\tfrac{3r_w^4}{r^4}\right)\cos 2\theta - p_w\tfrac{r_w^2}{r^2}\\ \tau_{r\theta} &= -\tfrac{\sigma_H-\sigma_h}{2}\left(1+\tfrac{2r_w^2}{r^2}-\tfrac{3r_w^4}{r^4}\right)\sin 2\theta \end{aligned}

любые одинаковые единицы напряжения; сжатие положительно

Мына жағдайда орындалады: те же, что у формулы для стенки; при r=rwr = r_w получаются σr=pw\sigma_r = p_w, формула для σθ\sigma_\theta на стенке и τrθ=0\tau_{r\theta} = 0.

Механика горных пород для гидроразрыва →
Предельное отношение эффективных напряжений при трении на разломах
σ1−ppσ3−pp≤(μ2+1+μ)2\frac{\sigma_1 - p_p}{\sigma_3 - p_p} \le \left(\sqrt{\mu^2 + 1} + \mu\right)^2

любые одинаковые единицы давления

Мына жағдайда орындалады: в массиве есть разломы любой ориентации без сцепления; эффективное напряжение — по Терцаги (коэффициент Био для разломов принимают равным 1).

Механика горных пород для гидроразрыва →
Эффективное напряжение на проппанте
σэф=pc−pwf\sigma_{\text{эф}} = p_c - p_{wf}

любые одинаковые единицы давления

Мына жағдайда орындалады: у ствола давление в трещине равно забойному; дальше по трещине оно выше, и нагрузка меньше; при истощении давление закрытия снижается вместе с пластовым — на A ΔpˉA\,\Delta\bar p, где AA — пороупругий коэффициент (Механика горных пород для гидроразрыва); обычно это меньше, чем снижение забойного давления, поэтому проектную нагрузку считают по начальному давлению закрытия и наименьшему планируемому забойному — с запасом.

Размеры и проводимость трещины: проектирование →
Ширина и проводимость закреплённой трещины
wp=Cp(1−ϕp) ρp,(kfw)пл=kf,лаб(σэф) R (wp−2e)w_p = \frac{C_p}{(1 - \phi_p)\,\rho_p},\qquad (k_f w)_{\text{пл}} = k_{f,\text{лаб}}(\sigma_{\text{эф}})\,R\,(w_p - 2e)

проводимость — в мД·м при kk в мД и ширине в м

Мына жағдайда орындалады: проппант лёг равномерно, без оседания к низу трещины; множители RR и ee определяют испытаниями с жидкостью разрыва и с образцами породы пласта; без испытаний их принимают с запасом.

Размеры и проводимость трещины: проектирование →
Эффективная проницаемость упаковки с инерционными потерями
kf,эф=kf1+23 Fow,Fow=βf kf ρ vwμ,vw=qB/2hp wpk_{f,\text{эф}} = \frac{k_f}{1 + \tfrac23\,\text{Fo}_w},\qquad \text{Fo}_w = \frac{\beta_f\,k_f\,\rho\,v_w}{\mu},\qquad v_w = \frac{qB/2}{h_p\,w_p}

СИ

Мына жағдайда орындалады: приток в трещину равномерен по длине — так ведёт себя трещина высокой проводимости; при низкой проводимости приток сосредоточен у ствола, и потери ближе к значению у ствола kf/(1+Fow)k_f/(1 + \text{Fo}_w); βf\beta_f упаковки растёт при дроблении зёрен, остатках геля и при наличии второй фазы.

Размеры и проводимость трещины: проектирование →
Число проппанта
Nprop=Ix2 FCD=4xf2xe2⋅kfwk xf=2 kf Vpropk VresN_{prop} = I_x^2\,F_{CD} = \frac{4x_f^2}{x_e^2}\cdot\frac{k_f w}{k\,x_f} = \frac{2\,k_f\,V_{prop}}{k\,V_{res}}

безразмерная; kfk_f и kk в одних единицах

Мына жағдайда орындалады: вертикальная трещина в центре квадратной области дренирования, параллельная её стороне; в VpropV_{prop} идёт только проппант в пределах пласта: если трещина выше пласта, а проппант распределён по высоте равномерно, — доля h/hfh/h_f от всего проппанта; kfk_f — проницаемость упаковки в пласте с учётом нагрузки, остатков полимера и вдавливания (предыдущий раздел).

Размеры и проводимость трещины: проектирование →
Наибольшая продуктивность при малом числе проппанта
1JD,max⁡≈0,99−12ln⁡Nprop,FCD,opt≈1,6–1,7\frac{1}{J_{D,\max}} \approx 0{,}99 - \frac12\ln N_{prop},\qquad F_{CD,opt} \approx 1{,}6\text{–}1{,}7

безразмерная; логарифм натуральный

Мына жағдайда орындалады: NpropN_{prop} до 0,1: при 0,1 формула завышает JD,max⁡J_{D,\max} на 0,7 %, при 0,3 — уже на 3 %; квадратная область, псевдоустановившийся режим; вывод — в разделе для углублённого изучения.

Размеры и проводимость трещины: проектирование →
Размеры трещины по числу проппанта и выбранной безразмерной проводимости
xf=kf Vprop2FCD k h,w=FCD k Vprop2kf hx_f = \sqrt{\frac{k_f\,V_{prop}}{2F_{CD}\,k\,h}},\qquad w = \sqrt{\frac{F_{CD}\,k\,V_{prop}}{2k_f\,h}}

согласованные; kfk_f и kk в одних единицах

Мына жағдайда орындалады: следует из Vprop=2xfwhV_{prop} = 2x_f w h и FCD=kfw/(kxf)F_{CD} = k_f w/(k x_f); эквивалентная запись — xf=Ixxe/2x_f = I_x x_e/2 при Ix=Nprop/FCDI_x = \sqrt{N_{prop}/F_{CD}}.

Размеры и проводимость трещины: проектирование →
Ширина трещины у ствола и чистое давление по модели PKN
ww≈2,75(μ qi xfE′)1/4,wˉ≈0,58 ww,pnet=E′ ww2hfw_w \approx 2{,}75\left(\frac{\mu\,q_i\,x_f}{E'}\right)^{1/4},\qquad \bar w \approx 0{,}58\,w_w,\qquad p_{net} = \frac{E'\,w_w}{2h_f}

СИ (любая согласованная система единиц)

Мына жағдайда орындалады: трещина много длиннее своей высоты, высота постоянна, каждое вертикальное сечение — эллипс; ньютоновская жидкость, ламинарное течение; множитель 2,75 и отношение 0,58 — середина между решением без утечек (2,82 и 0,59) и предельным случаем сильных утечек (2,68 и 0,57); вывод — в разделе для углублённого изучения. Для жидкостей, вязкость которых падает с ростом скорости сдвига, берут кажущуюся вязкость при скорости сдвига в трещине; множитель отнесён к полному расходу в оба крыла. Если подставлять расход в одно крыло, множитель в 21/4≈1,192^{1/4} \approx 1{,}19 раза больше — около 3,27; в таком виде формулу часто и приводят, поэтому, сравнивая формулы, проверяйте, к какому расходу отнесён множитель.

Размеры и проводимость трещины: проектирование →
Баланс объёмов при закачке с утечками
qit=Af wˉ+2Al(κ CLt+Sp),η=Af wˉqitq_i t = A_f\,\bar w + 2A_l\left(\kappa\,C_L\sqrt t + S_p\right),\qquad \eta = \frac{A_f\,\bar w}{q_i t}

СИ (любая согласованная система единиц)

Мына жағдайда орындалады: закон Картера: скорость утечки обратно пропорциональна корню из времени, прошедшего с открытия участка стенки; через обе стенки, отсюда множитель 2; утечки идут только в проницаемый пласт; проппант входит в закачанный объём, но не утекает; κ\kappa получается усреднением 2t−τ2\sqrt{t-\tau} по площади при росте площади по степенному закону; в примерах урока κ\kappa = 1,5, мгновенная утечка не учитывается.

Размеры и проводимость трещины: проектирование →
Эффективность и коэффициент утечки по калибровочной закачке
η=3π16Gc1+3π16Gc,CL=m hf2E′ hltp\eta = \frac{\frac{3\pi}{16}G_c}{1 + \frac{3\pi}{16}G_c},\qquad C_L = \frac{m\,h_f^2}{E'\,h_l\sqrt{t_p}}

СИ; CLC_L в м/с1/2^{1/2}

Мына жағдайда орындалады: функция GG той формы, что в уроке о механике пород: площадь трещины при закачке росла пропорционально времени; со второй предельной формой эффективность получается на несколько процентов меньше; после остановки площадь не меняется, мгновенной утечки нет, к моменту закрытия в трещине не осталось жидкости; податливость трещины PKN; для коротких трещин нужна податливость модели KGD, и коэффициент утечки получится другим.

Размеры и проводимость трещины: проектирование →
Доля подушки и закон роста концентрации проппанта
ε=1−η1+η,Vpad=εVi,c(t)=ce(t−tpadte−tpad)ε,ce=MηVi\varepsilon = \frac{1 - \eta}{1 + \eta},\qquad V_{pad} = \varepsilon V_i,\qquad c(t) = c_e\left(\frac{t - t_{pad}}{t_e - t_{pad}}\right)^{\varepsilon},\qquad c_e = \frac{M}{\eta V_i}

согласованные

Мына жағдайда орындалады: утечки по закону Картера, площадь трещины растёт во время закачки; проппант не оседает и движется вместе с жидкостью; концентрация — масса проппанта на единицу объёма смеси; в промысловой практике её задают на объём чистой жидкости: cж=c/(1−c/ρp)c_{\text{ж}} = c/(1 - c/\rho_p); эффективность берут из расчёта основной закачки, откалиброванного по калибровочной; при сомнениях подушку немного увеличивают — стоп обходится дороже лишней жидкости.

Размеры и проводимость трещины: проектирование →
Продуктивность короткой трещины через число проппанта
1JD=ln⁡2−1,311−12ln⁡Nprop+[12ln⁡FCD+ln⁡xfrw′(FCD)]\frac{1}{J_D} = \ln2 - 1{,}311 - \frac12\ln N_{prop} + \left[\frac12\ln F_{CD} + \ln\frac{x_f}{r_w'}\left(F_{CD}\right)\right]

безразмерная; логарифмы натуральные

Мына жағдайда орындалады: трещина мала по сравнению с областью (IxI_x меньше примерно 0,3), псевдоустановившийся режим, квадратная область.

Размеры и проводимость трещины: проектирование →
Показатели степенной модели по двум точкам и кажущаяся вязкость
n=ln⁡(τ2/τ1)ln⁡(γ˙2/γ˙1),K=τ1γ˙1 n,μa=Kγ˙ n−1n = \frac{\ln(\tau_2/\tau_1)}{\ln(\dot\gamma_2/\dot\gamma_1)},\qquad K = \frac{\tau_1}{\dot\gamma_1^{\,n}},\qquad \mu_a = K\dot\gamma^{\,n-1}

СИ; KK в Па·сⁿ (в промысловых единицах — lbf·sⁿ/ft², 1 lbf·sⁿ/ft² = 47,88 Па·сⁿ)

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

Жидкости разрыва и гидроразрыв с остановкой роста трещины →
Скорость сдвига у стенки щели и трубы для степенной жидкости
γ˙щ=2n+13n⋅6vw,γ˙тр=3n+14n⋅8vd\dot\gamma_{\text{щ}} = \frac{2n+1}{3n}\cdot\frac{6v}{w},\qquad \dot\gamma_{\text{тр}} = \frac{3n+1}{4n}\cdot\frac{8v}{d}

СИ

Мына жағдайда орындалады: ламинарное течение; щель — две параллельные стенки (сечение трещины на деле эллиптическое, и формула даёт характерную скорость сдвига); при n=1n = 1 получаются ньютоновские 6v/w6v/w и 8v/d8v/d (скорость сдвига в трубе); множители учитывают более плоский профиль скорости разжижающейся жидкости.

Жидкости разрыва и гидроразрыв с остановкой роста трещины →
Опыт на фильтрацию: коэффициент корки и мгновенная утечка
VA=Sp+2Cwt,Cw=m2A,Sp=bA\frac{V}{A} = S_p + 2C_w\sqrt t,\qquad C_w = \frac{m}{2A},\qquad S_p = \frac{b}{A}

любые согласованные: CwC_w получается в единицах длины на корень из единицы времени

Мына жағдайда орындалады: скорость утечки через корку Cw/tC_w/\sqrt t — отсюда множитель 2 при интегрировании: ∫0tCw dt′/t′=2Cwt\int_0^t C_w\,dt'/\sqrt{t'} = 2C_w\sqrt t; перепад давления в опыте близок к пластовому (иначе коэффициент пересчитывают; у полимерных корок он растёт с перепадом слабее, чем как корень из него); образец той же породы и проницаемости.

Жидкости разрыва и гидроразрыв с остановкой роста трещины →
Коэффициенты утечки, ограниченной фильтратом и сжатием пласта, и их сочетание
Cv=k ϕ Δp2μf,Cc=Δpk ϕ ctπ μr,CL=2CcCwCvCwCv+Cw2Cv2+4Cc2(Cw2+Cv2)C_v = \sqrt{\frac{k\,\phi\,\Delta p}{2\mu_f}},\qquad C_c = \Delta p\sqrt{\frac{k\,\phi\,c_t}{\pi\,\mu_r}},\qquad C_L = \frac{2C_cC_wC_v}{C_wC_v + \sqrt{C_w^2C_v^2 + 4C_c^2\left(C_w^2 + C_v^2\right)}}

СИ (kk в м²: 1 мД = 0,987⋅10−150{,}987\cdot10^{-15} м²); в промысловых единицах коэффициенты чаще дают в ft/мин1/2^{1/2}: 1 м/с1/2^{1/2} = 25,4 ft/мин1/2^{1/2}

Мына жағдайда орындалады: линейная (одномерная) фильтрация в пласт, перепад давления постоянен; фильтрат вытесняет пластовый флюид поршнем; общий коэффициент — точное сочетание трёх сопротивлений в ряд при одном общем перепаде, если коэффициенты корки и фильтрата растут как корень из своей доли перепада, а коэффициент сжатия — пропорционально ей; часто пользуются и приближением 1/CL=1/Cw+1/Cv+1/Cc1/C_L = 1/C_w + 1/C_v + 1/C_c, которое всегда даёт меньшее значение — в примере ниже на 15 %; вывод формул — в разделе для углублённого изучения.

Жидкости разрыва и гидроразрыв с остановкой роста трещины →
Скорость оседания зерна в ньютоновской и степенной жидкости
vs=g Δρ d218 μ,vs=[g Δρ d n+118⋅3 n−1K]1/nv_s = \frac{g\,\Delta\rho\,d^2}{18\,\mu},\qquad v_s = \left[\frac{g\,\Delta\rho\,d^{\,n+1}}{18\cdot 3^{\,n-1}K}\right]^{1/n}

СИ

Мына жағдайда орындалады: одиночное сферическое зерно, малое число Рейнольдса ρfvsd/μ≪1\rho_f v_s d/\mu \ll 1; стенки трещины далеко, жидкость неподвижна; вторая формула — закон Стокса с кажущейся вязкостью при скорости сдвига 3vs/d3v_s/d; при n=1n = 1 совпадает с первой. Это оценка: при очень малых скоростях сдвига реальная вязкость перестаёт расти, и формула занижает скорость оседания; в концентрированной смеси зёрна мешают друг другу и оседают медленнее — по эмпирическому закону Ричардсона и Заки при малых числах Рейнольдса скорость умножается примерно на (1−cv)4,65(1 - c_v)^{4{,}65}, где cvc_v — объёмная доля проппанта.

Жидкости разрыва и гидроразрыв с остановкой роста трещины →
Утечки и объём трещины после остановки её роста
ΔVL=2AlCLtост[g(δ)−π2],g(δ)=δ+(1+δ)arcsin⁡11+δ,δ=t−tостtост\Delta V_L = 2A_lC_L\sqrt{t_{\text{ост}}}\left[g(\delta) - \frac{\pi}{2}\right],\quad g(\delta) = \sqrt\delta + (1 + \delta)\arcsin\frac{1}{\sqrt{1 + \delta}},\quad \delta = \frac{t - t_{\text{ост}}}{t_{\text{ост}}}

СИ

Мына жағдайда орындалады: до остановки площадь трещины росла как корень из времени (предел сильных утечек), каждый участок стенки фильтрует по закону Картера с момента своего открытия; после остановки площадь не меняется; мгновенная утечка вся случилась до остановки; та же функция описывает утечки после остановки насосов при калибровочной закачке: g(δ)−π/2=π4G(δ)g(\delta) - \pi/2 = \tfrac{\pi}{4}G(\delta), где GG — вторая предельная форма функции Нолте, для трещины, площадь которой росла как корень из времени (функция G там записана для роста пропорционально времени); вывод — в разделе для углублённого изучения.

Жидкости разрыва и гидроразрыв с остановкой роста трещины →
Концентрация проппанта по плотности смеси
cv=ρсм−ρжρp−ρж,c=cv ρp,cж=c1−cvc_v = \frac{\rho_{\text{см}} - \rho_{\text{ж}}}{\rho_p - \rho_{\text{ж}}},\qquad c = c_v\,\rho_p,\qquad c_{\text{ж}} = \frac{c}{1 - c_v}

СИ; концентрацию на кубометр жидкости в промысловой практике дают в фунтах на галлон: 1 фунт/галлон = 119,8 кг/м³

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

Жидкости разрыва и гидроразрыв с остановкой роста трещины →

Белгілер

БелгіМағынасы
pˉ\bar p
Δq\Delta q
Δp\Delta p
ΔN\Delta N
Δt\Delta t
wˉ\bar w
pˉ, pwf\bar p,\ p_{wf}
[Ba2+], [SO42−][\mathrm{Ba^{2+}}],\ [\mathrm{SO_4^{2-}}]
Δq0\Delta q_0
0,8090{,}809
Δpзап\Delta p_{\text{зап}}
Δpперф\Delta p_{\text{перф}}
Δpпз\Delta p_{\text{пз}}
Δpтр\Delta p_{\text{тр}}
aa
AA
a1, a2a_1,\ a_2
AпрA_{\text{пр}}
AfA_f
AlA_l
a, ba,\ b
Al4Si4O10(OH)8\mathrm{Al_4Si_4O_{10}(OH)_8}
AlF2+\mathrm{AlF_2^{+}}
bsb_s
CC
cec_e
CLC_L
CpC_p
ctc_t
cvc_v
Cv, Cc, CLC_v,\ C_c,\ C_L
CwC_w
c, cжc,\ c_{\text{ж}}
c(t)c(t)
CaCl2, MgCl2\mathrm{CaCl_2},\ \mathrm{MgCl_2}
CaCO3\mathrm{CaCO_3}
CaF2\mathrm{CaF_2}
CaMg(CO3)2\mathrm{CaMg(CO_3)_2}
dd
D1D_1
dк, lкd_{\text{к}},\ l_{\text{к}}
d, ld,\ l
DaDa
ee
EE
E′E'
FF
FCDF_{CD}
Fow\text{Fo}_w
gg
GG
GcG_c
g(δ)g(\delta)
hh
HH
H2SiF6, SiF62−\mathrm{H_2SiF_6},\ \mathrm{SiF_6^{2-}}
hfh_f
hf, hlh_f,\ h_l
hp, wph_p,\ w_p
h, xfh,\ x_f
h, ϕh,\ \phi
IxI_x
Δtj\Delta t_j
JJ
JDJ_D
JD,доJ_{D,\text{до}}
JD,max⁡J_{D,\max}
JD(t)J_D(t)
kk
kdk_d
kfk_f
kfwk_f w
kf, kk_f,\ k
kf, k, hk_f,\ k,\ h
kf,лабk_{f,\text{лаб}}
kf,эфk_{f,\text{эф}}
kf,chw, kfwk_{f,ch}w,\ k_f w
kik_i
kiдо, kiпослеk_i^{\text{до}},\ k_i^{\text{после}}
kj, hjk_j,\ h_j
ksk_s
KspK_{sp}
k, hk,\ h
k, h, re, rw, Sk,\ h,\ r_e,\ r_w,\ S
k, h, μ, pˉ, re, rwk,\ h,\ \mu,\ \bar p,\ r_e,\ r_w
k, kпk,\ k_{\text{п}}
k, kdk,\ k_d
k, ϕk,\ \phi
k, ϕ, ct, tk,\ \phi,\ c_t,\ t
k, rwk,\ r_w
Ksp′K'_{sp}
k(r)k(r)
ΔVL\Delta V_L
L, LiL,\ L_i
ln⁡(xf/rw′)\ln(x_f/r_w')
ln⁡(xf/rw′)∣FCD\ln(x_f/r_w')\big|_{F_{CD}}
  • значение универсальной кривой при данной безразмерной проводимости: от ln⁡2≈0,693\ln 2 \approx 0{,}693 (бесконечная проводимость) до нескольких единиц (низкая проводимость), —Приток к скважине после гидроразрыва
ln⁡(xf/rw′)(FCD)\ln(x_f/r_w')(F_{CD})
mm
MM
Mм, MкM_{\text{м}},\ M_{\text{к}}
m, bm,\ b
nn
NкN_{\text{к}}
nкаж(t)n_{\text{каж}}(t)
NpropN_{prop}
Nprop, FCDN_{prop},\ F_{CD}
n, Kn,\ K
Δpp\Delta p_p
pуp_{\text{у}}
pbp_b
pcp_c
pfp_f
pISIPp_{ISIP}
pnetp_{net}
pnet, hfp_{net},\ h_f
ppp_p
pwp_w
pw,jp_{w,j}
pwDp_{wD}
pwD, tDxfp_{wD},\ t_{Dxf}
pwD(t)p_{wD}(t)
pwfp_{wf}
PVbtPV_{bt}
PVbt(vi)PV_{bt}(v_i)
qq
q1, q2q_1,\ q_2
qкрq_{\text{кр}}
  • критический расход воды на метр интервала (в пластовых условиях; для воды почти равен расходу на поверхности), м³/с на мИсточники повреждения пласта
qпрq_{\text{пр}}
qiq_i
qi, tq_i,\ t
qj, qq_j,\ q
qmax⁡q_{\max}
qwhq_{wh}
q, B, μq,\ B,\ \mu
qBqB
rr
RR
r1, r2r_1,\ r_2
rв, rwr_{\text{в}},\ r_w
rdr_d
re, rwr_e,\ r_w
re, rw, Sr_e,\ r_w,\ S
RiR_i
ri−1, rir_{i-1},\ r_i
roptr_{opt}
rwr_w
rw, rr_w,\ r
rw′r_w'
rwh, rwr_{wh},\ r_w
S1S_1
S1, S2S_1,\ S_2
SвS_{\text{в}}
SкажS_{\text{каж}}
SпS_{\text{п}}
schs_{ch}
SdS_d
SfS_f
sfss_{fs}
SjS_j
SpS_p
SI\mathrm{SI}
SiO2\mathrm{SiO_2}
tt
ΔNT\Delta N_T
TT
tокt_{\text{ок}}
tостt_{\text{ост}}
tDAt_{DA}
tDxft_{Dxf}
tpt_p
tpad, tet_{pad},\ t_e
uu
uкрu_{\text{кр}}
  • критическая скорость фильтрации по испытанию на керне: расход, при котором проницаемость начала падать, делённый на площадь сечения образца, м/сИсточники повреждения пласта
vv
VV
viv_i
vopt, PVbt,optv_{opt},\ PV_{bt,opt}
Vpad, ViV_{pad},\ V_i
VpropV_{prop}
VresV_{res}
vsv_s
vwv_w
vwhv_{wh}
V/hV/h
ww
WW
wiw_i
wmax⁡w_{\max}
wpw_p
www_w
w, dw,\ d
xx
XX
xex_e
xfx_f
xf, kx_f,\ k
xf, wx_f,\ w
xf, hfx_f,\ h_f
xf, rw, rw′x_f,\ r_w,\ r_w'
xsx_s
X, S1, S2X,\ S_1,\ S_2
yey_e
α\alpha
α, ν\alpha,\ \nu
β\beta
βf\beta_f
γ\gamma
γ˙\dot\gamma
Γ(5/4)\Gamma(5/4)
γ˙щ, γ˙тр\dot\gamma_{\text{щ}},\ \dot\gamma_{\text{тр}}
γBa, γSO4\gamma_{\text{Ba}},\ \gamma_{\text{SO}_4}
Δρ\Delta\rho
Δσ\Delta\sigma
Δσh\Delta\sigma_h
ε\varepsilon
εh, εH\varepsilon_h,\ \varepsilon_H
η\eta
θ\theta
κ\kappa
μ\mu
μ, B\mu,\ B
μв, μ\mu_{\text{в}},\ \mu
μa\mu_a
μf, μr\mu_f,\ \mu_r
ν\nu
νм, νк\nu_{\text{м}},\ \nu_{\text{к}}
ρ\rho
ρ, μ\rho,\ \mu
ρк, ρм\rho_{\text{к}},\ \rho_{\text{м}}
ρсм, ρж, ρp\rho_{\text{см}},\ \rho_{\text{ж}},\ \rho_p
σ1, σ3\sigma_1,\ \sigma_3
σэф\sigma_{\text{эф}}
σh\sigma_h
σh, σH\sigma_h,\ \sigma_H
σH, σh\sigma_H,\ \sigma_h
σr\sigma_r
σr, σθ, τrθ\sigma_r,\ \sigma_\theta,\ \tau_{r\theta}
σv\sigma_v
σθ\sigma_\theta
τ1, τ2\tau_1,\ \tau_2
ϕ\phi
ϕ, X, x, r, rw\phi,\ X,\ x,\ r,\ r_w
ϕ, μ, ct\phi,\ \mu,\ c_t
ϕp, ρp\phi_p,\ \rho_p