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

Формулы: Корреляция и карты

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

Нормировка гамма-каротажа по опорной скважине
GRн=GRclоп+IGR(GRshоп−GRclоп),IGR=GR−GRclGRsh−GRclGR_{\text{н}} = GR_{cl}^{\text{оп}} + I_{GR}\left(GR_{sh}^{\text{оп}} - GR_{cl}^{\text{оп}}\right),\qquad I_{GR} = \frac{GR - GR_{cl}}{GR_{sh} - GR_{cl}}

единицы API; индекс безразмерен

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

Корреляция скважин →
Толщины интервала в наклонной скважине
h=Δlcos⁡γ,hв=hcos⁡δ,ΔH=Δlcos⁡θh = \Delta l\cos\gamma,\qquad h_{\text{в}} = \frac{h}{\cos\delta},\qquad \Delta H = \Delta l\cos\theta

длины в одних единицах; углы в градусах

Справедливо, если: пласт плоский, кровля и подошва параллельны; ствол на интервале прямолинейный — иначе интервал делят на части; азимуты пластов и ствола отсчитаны от одного севера; формула для cos⁡γ\cos\gamma — из темы векторный расчёт.

Корреляция скважин →
Линия корреляции по методу наименьших квадратов
zБ=zˉБ+b (zА−zˉА),b=∑k(zА,k−zˉА)(zБ,k−zˉБ)∑k(zА,k−zˉА)2z_{\text{Б}} = \bar z_{\text{Б}} + b\,(z_{\text{А}} - \bar z_{\text{А}}),\qquad b = \frac{\sum_k (z_{\text{А},k} - \bar z_{\text{А}})(z_{\text{Б},k} - \bar z_{\text{Б}})}{\sum_k (z_{\text{А},k} - \bar z_{\text{А}})^2}

метры; одинаковая привязка глубин в обеих скважинах (по вертикали)

Справедливо, если: события одновозрастны в обеих скважинах; в подбор берут надёжные события одного отрезка без перерывов; отношение скоростей накопления на отрезке постоянно; где оно меняется или есть перерыв, линию строят по частям; ошибки глубин событий в Б больше, чем в А (например, шлам против керна): поэтому подбирают zБz_{\text{Б}} по zАz_{\text{А}}.

Корреляция скважин →
Вероятность, что тело доходит до соседней скважины
P=1−d sin⁡αW(P≥0)P = 1 - \frac{d\,\sin\alpha}{W}\quad (P \ge 0)

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

Справедливо, если: тело прямолинейное, постоянной ширины, его положение относительно скважин случайно; формула грубеет при малых α\alpha, когда тело вытянуто почти вдоль линии скважин: тогда решают извилистость и длина тела, а не ширина; при dsin⁡α>Wd\sin\alpha > W вероятность нулевая: тело уже, чем расстояние между скважинами вдоль его поперечника.

Корреляция скважин →
Коэффициенты макронеоднородности
Kпес=∑hэф∑hобщ,Kр=1N∑i=1Nni,Kсв=SсвS,Kлв=SкSK_{\text{пес}} = \frac{\sum h_{\text{эф}}}{\sum h_{\text{общ}}},\qquad K_{\text{р}} = \frac{1}{N}\sum_{i=1}^{N} n_i,\qquad K_{\text{св}} = \frac{S_{\text{св}}}{S},\qquad K_{\text{лв}} = \frac{S_{\text{к}}}{S}

безразмерные; толщины и площади в одинаковых единицах

Справедливо, если: коэффициенты считают для пласта или горизонта, выделенного корреляцией; KпесK_{\text{пес}} зависит от граничных значений, по которым выделен коллектор (граничные значения); здесь он посчитан с весом толщин; KсвK_{\text{св}} и KлвK_{\text{лв}} берут с карт, построенных по скважинам, поэтому их точность ограничена плотностью сетки; в части методик площадь слияния относят не ко всей площади залежи, а к площади распространения коллекторов пласта, — сравнивая объекты, проверяют, как посчитан KсвK_{\text{св}}.

Строение коллектора и поведение пласта при разработке →
Меры связности
C=Vmax⁡Vк,Cw=VwVкC = \frac{V_{\max}}{V_{\text{к}}},\qquad C_w = \frac{V_w}{V_{\text{к}}}

доли единицы

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

Строение коллектора и поведение пласта при разработке →
Постоянная времени блока за перегородкой
τ=Vp ctT,T=kb Aμ b\tau = \frac{V_p\,c_t}{T},\qquad T = \frac{k_b\,A}{\mu\,b}

СИ

Справедливо, если: давление внутри блока одинаково — блок много проницаемее перегородки; давление по другую сторону постоянно: соседний пласт много больше блока или поддерживается закачкой; если блоки соизмеримы, VpV_p заменяют на V1V2/(V1+V2)V_1V_2/(V_1 + V_2); течение через перегородку однофазное и почти установившееся; сжимаемостью самой перегородки пренебрегают.

Строение коллектора и поведение пласта при разработке →
Где изолиния пересекает сторону треугольника
sc=L c−zazb−zas_c = L\,\frac{c - z_a}{z_b - z_a}

любые одинаковые единицы длины; значения — в одних единицах

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

Структурные карты: данные и построение изолиний →
Метод обратных расстояний
z^(x,y)=∑i=1nwi zi∑i=1nwi,wi=1di p\hat z(x, y) = \frac{\sum_{i=1}^{n} w_i\,z_i}{\sum_{i=1}^{n} w_i},\qquad w_i = \frac{1}{d_i^{\,p}}

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

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

Структурные карты: данные и построение изолиний →
Условие минимальной кривизны
∬(zxx2+2 zxy2+zyy2)dx dy→min⁡,z(xi,yi)=zi\iint \left(z_{xx}^2 + 2\,z_{xy}^2 + z_{yy}^2\right) dx\,dy \to \min,\qquad z(x_i, y_i) = z_i

согласованные единицы длины

Справедливо, если: поверхность гладкая и непрерывная; разломы в ней не предусмотрены — их задают отдельно, иначе программа «сгладит» уступ в изгиб; между точками данных такая поверхность подчиняется бигармоническому уравнению ∇4z=0\nabla^4 z = 0; программа решает его на сетке итерациями.

Структурные карты: данные и построение изолиний →
Изохора как разность поверхностей и изопахита
hв(x,y)=zп(x,y)−zк(x,y),h=hвcos⁡δbh_\text{в}(x, y) = z_\text{п}(x, y) - z_\text{к}(x, y),\qquad h = h_\text{в}\cos\delta_b

длины в одних единицах; углы в градусах

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

Разрезы по картам, совмещение карт и разломы на картах →
Амплитуды и угол сместителя по многоугольнику разлома
w=Δ⊥,t=zв−zл,tg δ=tww = \Delta_\perp,\qquad t = z_\text{в} - z_\text{л},\qquad \mathrm{tg}\,\delta = \frac{t}{w}

длины в метрах (на местности, а не на карте); углы в градусах

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

Разрезы по картам, совмещение карт и разломы на картах →
Сдвиг многоугольника подошвы относительно многоугольника кровли
Δq=hвtg δ±tg δb\Delta q = \frac{h_\text{в}}{\mathrm{tg}\,\delta \pm \mathrm{tg}\,\delta_b}

м; углы в градусах

Справедливо, если: «+», если пласт падает навстречу сместителю, «−», если в ту же сторону — так же, как в формуле вертикального смещения пласта в теме о разломах; пласт и сместитель плоские на отрезке Δq\Delta q, толщина постоянна; одинаково для обоих крыльев.

Разрезы по картам, совмещение карт и разломы на картах →
Смещение одноимённых изогипс вдоль разлома
Δs=ttg δ∥\Delta s = \frac{t}{\mathrm{tg}\,\delta_\parallel}

м

Справедливо, если: смещение по падению: точка висячего крыла и её «пара» в лежачем лежат друг против друга по нормали к разлому и различаются по глубине на tt; пласт наклонён вдоль разлома; если он простирается параллельно разлому (δ∥=0\delta_\parallel = 0), изогипсы не смещаются, а идут вдоль полосы, и амплитуду читают по числу изогипс, «пропущенных» поперёк неё.

Разрезы по картам, совмещение карт и разломы на картах →
Объём пород выше контакта по кривой «площадь — глубина»
GRV=∫zсвzВНКAк(z) dz\mathrm{GRV} = \int_{z_\text{св}}^{z_\text{ВНК}} A_\text{к}(z)\,dz

СИ: площадь в м², глубина в м, объём в м³; удобно: км² × м = млн м³

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

Общий объём пород по картам →
Площадь многоугольника по координатам вершин
A=12∣∑i=1n(xi yi+1−xi+1 yi)∣,(xn+1,yn+1)=(x1,y1)A = \frac12\left|\sum_{i=1}^{n} \left(x_i\,y_{i+1} - x_{i+1}\,y_i\right)\right|,\qquad (x_{n+1}, y_{n+1}) = (x_1, y_1)

координаты в метрах — площадь в м²

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

Общий объём пород по картам →
Объём слоя между соседними изогипсами: правила трапеций и пирамид
ΔVтр=Δz2 (Ai+Ai+1),ΔVпир=Δz3(Ai+Ai+1+Ai Ai+1)\Delta V_\text{тр} = \frac{\Delta z}{2}\,(A_i + A_{i+1}),\qquad \Delta V_\text{пир} = \frac{\Delta z}{3}\left(A_i + A_{i+1} + \sqrt{A_i\,A_{i+1}}\right)

площади и сечение в согласованных единицах

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

Общий объём пород по картам →
Шапка над первой изогипсой и формула Симпсона
Vш=A0hш3 (конус) … A0hш2 (параболоид),V=Δz3(A0+4A1+2A2+…+4An−1+An)V_\text{ш} = \frac{A_0 h_\text{ш}}{3}\ \text{(конус)}\ \ldots\ \frac{A_0 h_\text{ш}}{2}\ \text{(параболоид)},\qquad V = \frac{\Delta z}{3}\left(A_0 + 4A_1 + 2A_2 + \ldots + 4A_{n-1} + A_n\right)

согласованные единицы

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

Общий объём пород по картам →
Объём пород тонкого пласта по кривым кровли и подошвы
GRV=∫zсвzВНК[Aк(z)−Aп(z)]dz\mathrm{GRV} = \int_{z_\text{св}}^{z_\text{ВНК}} \left[A_\text{к}(z) - A_\text{п}(z)\right] dz

м², м, м³

Справедливо, если: площадь на подошве считается от того же свода: выше подошвы её кривая равна нулю; если толщина пласта по вертикали постоянна и равна hвh_\text{в}, подошва — кровля, опущенная на hвh_\text{в}, и Aп(z)=Aк(z−hв)A_\text{п}(z) = A_\text{к}(z - h_\text{в}); при переменной толщине нужна своя карта подошвы.

Общий объём пород по картам →
Поправочный коэффициент геометрии
GCF=GRVAВНК hв;параболоид: GCF=1−hв2H  (hв≤H),H2hв  (hв≥H)\mathrm{GCF} = \frac{\mathrm{GRV}}{A_\text{ВНК}\,h_\text{в}};\qquad \text{параболоид: } \mathrm{GCF} = 1 - \frac{h_\text{в}}{2H}\ \ (h_\text{в} \le H),\qquad \frac{H}{2h_\text{в}}\ \ (h_\text{в} \ge H)

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

Справедливо, если: формулы справа — для купола-параболоида, у которого площадь растёт линейно с глубиной; для конуса GCF=1−(1−hв/H)33hв/H\mathrm{GCF} = \dfrac{1 - (1 - h_\text{в}/H)^3}{3h_\text{в}/H} при hв≤Hh_\text{в} \le H и H/(3hв)H/(3h_\text{в}) при hв≥Hh_\text{в} \ge H; толщина постоянна, контакт горизонтален.

Общий объём пород по картам →

Обозначения

СимволЗначение и где встречается
Δl\Delta l
  • длина ствола между кровлей и подошвой интервала (разность глубин по стволу), мКорреляция скважин
ΔH\Delta H
z^\hat z
Δq\Delta q
Δs\Delta s
Δz\Delta z
ΔV\Delta V
(xi,yi), zi(x_i, y_i),\ z_i
zˉА, zˉБ\bar z_{\text{А}},\ \bar z_{\text{Б}}
AA
A0…AnA_0 \ldots A_n
A0, hшA_0,\ h_\text{ш}
AВНКA_\text{ВНК}
Aк(z)A_\text{к}(z)
Aк(z), Aп(z)A_\text{к}(z),\ A_\text{п}(z)
Ai, Ai+1A_i,\ A_{i+1}
bb
cc
CC
ctc_t
CwC_w
dd
did_i
GRGR
GRнGR_{\text{н}}
GRcl, GRshGR_{cl},\ GR_{sh}
  • показания против опорного чистого пласта и опорной глинистой пачки в нормируемой скважине, APIКорреляция скважин
GRclоп, GRshопGR_{cl}^{\text{оп}},\ GR_{sh}^{\text{оп}}
GRV\mathrm{GRV}
hh
HH
hвh_\text{в}
hэф, hобщh_{\text{эф}},\ h_{\text{общ}}
h, hвh,\ h_{\text{в}}
IGRI_{GR}
KлвK_{\text{лв}}
KпесK_{\text{пес}}
KрK_{\text{р}}
KсвK_{\text{св}}
kbk_b
LL
nn
ni, Nn_i,\ N
pp
PP
Sсв, Sк, SS_{\text{св}},\ S_{\text{к}},\ S
scs_c
tt
TT
tg δ∥\mathrm{tg}\,\delta_\parallel
VкV_{\text{к}}
Vmax⁡V_{\max}
VpV_p
VwV_w
ww
WW
xi, yix_i,\ y_i
zА,k, zБ,kz_{\text{А},k},\ z_{\text{Б},k}
zВНКz_\text{ВНК}
zк, zпz_\text{к},\ z_\text{п}
zл, zвz_\text{л},\ z_\text{в}
zсвz_\text{св}
zсв, zВНКz_\text{св},\ z_\text{ВНК}
za, zbz_a,\ z_b
ziz_i
zxx, zxy, zyyz_{xx},\ z_{xy},\ z_{yy}
z(x,y)z(x, y)
α\alpha
γ\gamma
  • угол между стволом и нормалью к напластованию: cos⁡γ=∣sin⁡δ sin⁡θ cos⁡(a−φ)−cos⁡δ cos⁡θ∣\cos\gamma = |\sin\delta\,\sin\theta\,\cos(a - \varphi) - \cos\delta\,\cos\theta|, °Корреляция скважин
δ\delta
δ, a\delta,\ a
Δ⊥\Delta_\perp
δb\delta_b
θ, φ\theta,\ \varphi
μ\mu
τ\tau