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

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

ГеологияСтатистика, объёмы и запасы Орташа 18 мин оқу

Распределения и достаточность выборки

Распределение (distribution) — закон, который описывает, как часто величина принимает те или иные значения. Гистограмма показывает распределение выборки; распределение всей совокупности описывают математической кривой — нормальной, логнормальной или другой. Достаточность выборки (sample sufficiency) — ответ на вопрос, сколько измерений нужно, чтобы среднему по ним можно было доверять с заданной точностью.

Не үшін қажет

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

Қысқаша негізгісі

  • Гистограмма выборки — приближение к плотности вероятности совокупности. Площадь под кривой плотности между двумя значениями — доля породы с такими значениями.
  • Нормальное распределение задают среднее μ\mu и стандартное отклонение σ\sigma: в пределах μ±σ\mu \pm \sigma лежит 68 % значений, в пределах μ±2σ\mu \pm 2\sigma — 95 %. К нему близки пористость одной породы и погрешности измерений.
  • Логнормально распределена величина, логарифм которой распределён нормально: проницаемость и вообще всё, что получается перемножением многих факторов. У такого распределения мода меньше медианы, а медиана меньше среднего.
  • На вероятностном графике нормальное распределение выглядит прямой. Изгибы и изломы выдают неподходящий масштаб или смесь пород.
  • Стандартная ошибка среднего s/Ns/\sqrt N убывает медленно: чтобы уточнить среднее вдвое, нужно вчетверо больше образцов. Для точности ±20 % с вероятностью 95 % нужно около (10 Cv)2(10\,C_v)^2 образцов.
  • Смещение — систематическая ошибка отбора — от числа образцов не зависит. Образцы только из крепкой, хорошей породы дают точное, но неверное среднее.
Негізгі қатынастар
z=x−μσz = \dfrac{x - \mu}{\sigma}

стандартизованное значение; доля значений меньше xx равна Φ(z)\Phi(z)

x50=eμ,xˉ=eμ+σ2/2x_{50} = \e^{\mu},\quad \bar x = \e^{\mu + \sigma^2/2}

медиана и среднее логнормального распределения (μ\mu и σ\sigma — для ln⁡x\ln x)

SE=sNSE = \dfrac{s}{\sqrt{N}}

стандартная ошибка среднего

xˉ±t sN\bar x \pm t\,\dfrac{s}{\sqrt{N}}

доверительный интервал среднего; t≈2t \approx 2 при 95 % и N>30N > 30

N0=(10 Cv)2N_0 = (10\,C_v)^2

число образцов для точности ±20 % (95 %)

Ps=200 CvNsP_s = \dfrac{200\,C_v}{\sqrt{N_s}}

точность среднего, %, достигнутая при NsN_s образцах

Жиі қателіктер
  • Принимать стандартное отклонение данных за погрешность среднего — или наоборот.
  • Применять правило «± 2σ» к проницаемости в миллидарси, а не к её логарифму.
  • Считать P50 ожидаемым значением: у скошенного вправо распределения среднее больше медианы.
  • Считать сто образцов из одной скважины сотней независимых сведений о пласте.
  • Думать, что большое число образцов исправляет неудачный отбор.

Сабақтан кейін сіз

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

От гистограммы к распределению

Гистограмма описывает выборку: сорок образцов, сто, тысячу (как её строят). Если бы образцов становилось всё больше, а классы — всё уже, ступенчатый контур гистограммы превратился бы в плавную кривую. Она описывает уже не выборку, а совокупность и называется плотностью вероятности (probability density function).

Высота кривой плотности сама вероятностью не является. Вероятность — это площадь под кривой: доля значений между aa и bb равна площади над отрезком от aa до bb, а вся площадь под кривой равна единице. Вероятность здесь — число от 0 (невозможно) до 1 (достоверно); для породы её удобно понимать как долю объёма с такими свойствами.

Накопленной кривой выборки соответствует функция распределения (cumulative distribution function) F(x)F(x) — вероятность того, что величина не превысит xx. Она растёт от 0 до 1. Вероятность превышения равна 1−F(x)1 - F(x); именно её имеют в виду, когда оценку запасов называют P90.

Выборка и совокупность: что чему соответствует
По выборке (оценка)
гистограмма
У совокупности (параметр)
плотность вероятности f(x)f(x)
По выборке (оценка)
накопленная кривая
У совокупности (параметр)
функция распределения F(x)F(x)
По выборке (оценка)
среднее xˉ\bar x
У совокупности (параметр)
среднее μ\mu
По выборке (оценка)
стандартное отклонение ss
У совокупности (параметр)
стандартное отклонение σ\sigma

Функцией распределения пользуются и в обратную сторону: по заданной вероятности находят значение — квантиль. На этом построен метод Монте-Карло: случайное число от 0 до 1 принимают за накопленную вероятность, по функции распределения находят отвечающую ей пористость, толщину или площадь и повторяют так тысячи раз (Вероятностная оценка запасов методом Монте-Карло).

Нормальное распределение

Нормальное распределение (normal, Gaussian distribution) — симметричная колоколообразная кривая. Оно возникает там, где величина складывается из многих независимых небольших вкладов, ни один из которых не преобладает: случайные погрешности измерений, средние по большому числу образцов.

Плотность нормального распределения
f(x)=1σ2π exp⁡[−(x−μ)22σ2]f(x) = \frac{1}{\sigma\sqrt{2\pi}}\,\exp\left[-\frac{(x - \mu)^2}{2\sigma^2}\right]
мұндағы:
  • f(x)f(x)плотность вероятности, 1/единицы величины
  • xxзначение величины, единицы величины
  • μ\muсреднее; у нормального распределения оно же медиана и мода, единицы величины
  • σ\sigmaстандартное отклонение, единицы величины

Бірліктер: xx, μ\mu и σ\sigma — в единицах измеряемой величины

Мына жағдайда орындалады
  • величина может принимать любые значения, в том числе отрицательные; для пористости или толщины это лишь приближение, пригодное, пока μ\mu больше примерно 3σ3\sigma;
  • распределение одновершинное и симметричное.

Кривую полностью задают два числа: μ\mu сдвигает колокол вдоль оси, σ\sigma растягивает его. По выборке их оценивают средним арифметическим xˉ\bar x и стандартным отклонением ss.

Стандартизованное значение
z=x−μσz = \frac{x - \mu}{\sigma}
мұндағы:
  • zzрасстояние от среднего в стандартных отклонениях, —
  • xxзначение величины, единицы величины
  • μ, σ\mu,\ \sigmaсреднее и стандартное отклонение распределения, единицы величины

Бірліктер: безразмерное

Мына жағдайда орындалады
  • распределение нормальное; тогда zz подчиняется стандартному нормальному распределению со средним 0 и отклонением 1.

Любая задача о нормальном распределении сводится к одной кривой: доля значений меньше xx равна Φ(z)\Phi(z) — функции стандартного нормального распределения, которую берут из таблиц или вычисляют в электронной таблице. Обратно: значение с накопленной вероятностью pp равно μ+zpσ\mu + z_p\sigma.

Какая доля значений нормального распределения лежит в интервале
Интервал
μ±0,674σ\mu \pm 0{,}674\sigma
Доля значений
50 %
Примечание
между квартилями
Интервал
μ±σ\mu \pm \sigma
Доля значений
68,3 %
Примечание
около двух третей
Интервал
μ±1,282σ\mu \pm 1{,}282\sigma
Доля значений
80 %
Примечание
между 10-м и 90-м процентилями
Интервал
μ±1,645σ\mu \pm 1{,}645\sigma
Доля значений
90 %
Примечание
между 5-м и 95-м процентилями
Интервал
μ±1,960σ\mu \pm 1{,}960\sigma
Доля значений
95 %
Примечание
основа 95-процентных интервалов
Интервал
μ±2σ\mu \pm 2\sigma
Доля значений
95,4 %
Примечание
удобное округление
Интервал
μ±3σ\mu \pm 3\sigma
Доля значений
99,7 %
Примечание
почти все значения
Нормальное распределение пористости: доли значений в пределах одного и двух стандартных отклонений12,715,819,022,225,300,050,10Пористость, %Плотность вероятности, 1/%μ − 2σμ − σμμ + σμ + 2σ68 %95 %

Колоколообразная кривая плотности вероятности пористости с вершиной при 19,0 процента. Вертикальные штриховые линии проведены при 12,7; 15,8; 19,0; 22,2 и 25,3 процента и подписаны как среднее минус два сигма, минус сигма, среднее, плюс сигма, плюс два сигма. Площадь под кривой между 15,8 и 22,2 закрашена темнее и подписана 68 процентов, между 12,7 и 25,3 — светлее и подписана 95 процентов.

ЕсептеуНормальная кривая со средним 19,0 % и стандартным отклонением 3,16 % — параметрами пористости десяти образцов из примера ниже. Тёмная полоса — интервал μ±σ\mu \pm \sigma, светлая — μ±2σ\mu \pm 2\sigma.

В данных о пласте к нормальному распределению близки четыре вида величин.

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

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

Есептеу мысалыДоля породы ниже граничного значения пористости
Бастапқы деректер (шартты)
Пористость пласта
близка к нормальной; по десяти образцам ϕˉ=19,0\bar\phi = 19{,}0 %, s=3,162s = 3{,}162 %
Граничное значение (cut-off)
15 %
Шешім
  1. Стандартизуем границу: z=(15−19,0)/3,162=−1,265z = (15 - 19{,}0)/3{,}162 = -1{,}265.
  2. По таблице нормального распределения Φ(−1,265)=0,103\Phi(-1{,}265) = 0{,}103: около 10 % породы имеет пористость ниже 15 %.
  3. 10-й и 90-й процентили: 19,0−1,282⋅3,162=14,919{,}0 - 1{,}282\cdot3{,}162 = 14{,}9 % и 19,0+1,282⋅3,162=23,119{,}0 + 1{,}282\cdot3{,}162 = 23{,}1 %.
Жауап: Ниже граничного значения — около 10 % породы; 80 % значений лежат между 14,9 и 23,1 %.

Бұл нені білдіреді. В самой выборке ниже 15 % один образец из десяти, но по одному образцу о доле судить нельзя; кривая распределения использует все десять значений сразу. Независимая проверка: по сорока образцам того же пласта 10-й и 90-й процентили получились 14,4 и 23,0 % (пример с гистограммой). Оценка годится, пока распределение действительно близко к нормальному, а образцы представляют весь пласт. Данные примера условные.

Логнормальное распределение

Величина распределена логнормально (log-normal), если нормально распределён её логарифм. Нормальное распределение возникает при сложении многих независимых влияний, логнормальное — при их перемножении: логарифм произведения равен сумме логарифмов.

Логнормально — точно или приближённо — распределены многие величины, с которыми работают геологи и разработчики.

  • Проницаемость — главный пример: она зависит от произведения факторов (размер зёрен, пористость, сортировка, цементация).
  • Размеры зёрен и поровых каналов — поэтому в гранулометрии размеры откладывают в логарифмической шкале.
  • Запасы — результат перемножения объёма пород, доли коллектора, пористости, насыщенности и пересчётного коэффициента (Подсчёт запасов объёмным методом). Даже если каждый сомножитель распределён симметрично, распределение произведения скошено вправо.
  • Размеры залежей в нефтегазоносном бассейне — приближённо: мелких много, крупных мало.
Медиана, среднее и мода логнормального распределения
x50=eμ,xˉ=eμ+σ2/2,xмод=eμ−σ2x_{50} = \e^{\mu}, \qquad \bar x = \e^{\mu + \sigma^2/2}, \qquad x_{\text{мод}} = \e^{\mu - \sigma^2}
мұндағы:
  • μ\muсреднее натурального логарифма величины, —
  • σ\sigmaстандартное отклонение натурального логарифма величины, —
  • x50x_{50}медиана; она же среднее геометрическое, единицы величины
  • xˉ\bar xсреднее арифметическое, единицы величины
  • xмодx_{\text{мод}}мода — наиболее вероятное значение, единицы величины

Бірліктер: xx — в единицах величины (например, мД); μ\mu и σ\sigma относятся к ln⁡x\ln x и безразмерны

Мына жағдайда орындалады
  • логарифм величины распределён нормально;
  • μ\mu и σ\sigma вычислены по натуральным логарифмам. Для десятичных: μln⁡=2,303 μlg⁡\mu_{\ln} = 2{,}303\,\mu_{\lg}, σln⁡=2,303 σlg⁡\sigma_{\ln} = 2{,}303\,\sigma_{\lg}.

Всегда мода меньше медианы, а медиана меньше среднего, и разрыв между ними растёт с σ\sigma. При σ=1\sigma = 1 среднее в e0,5≈1,65\e^{0{,}5} \approx 1{,}65 раза больше медианы, а мода в e≈2,7\e \approx 2{,}7 раза меньше её. По выборке μ\mu и σ\sigma оценивают средним и стандартным отклонением логарифмов.

Процентили и коэффициент вариации логнормального распределения
xp=eμ+zpσ,Cv=eσ2−1x_p = \e^{\mu + z_p\sigma}, \qquad C_v = \sqrt{\e^{\sigma^2} - 1}
мұндағы:
  • xpx_pзначение, ниже которого лежит доля pp значений, единицы величины
  • zpz_pквантиль стандартного нормального распределения: −1,282 для 10 %, 0 для 50 %, +1,282 для 90 %, —
  • CvC_vкоэффициент вариации самой величины (не её логарифма), —
  • μ, σ\mu,\ \sigmaсреднее и стандартное отклонение натурального логарифма, —

Бірліктер: xpx_p — в единицах величины; остальное безразмерно

Мына жағдайда орындалады
  • логарифм величины распределён нормально.

Процентили расположены симметрично не по разности, а по отношению: 90-й во столько же раз больше медианы, во сколько 10-й меньше её, — в e1,282σ\e^{1{,}282\sigma} раз. Коэффициент вариации зависит только от σ\sigma.

Логнормальное распределение проницаемости: мода, медиана и среднее010020030040000,0040,0080,012Проницаемость, мДПлотность вероятности, 1/мДмода 21 мДмедиана 61 мДсреднее 105 мД

Кривая плотности вероятности проницаемости от 0 до 400 миллидарси: крутой подъём до вершины при 21 миллидарси и длинный пологий спад вправо. Отмечены три точки: мода 21 миллидарси на вершине, медиана 61 миллидарси на склоне и среднее 105 миллидарси ещё правее и ниже.

ЕсептеуЛогнормальная кривая с параметрами проницаемости десяти образцов из примера ниже (μ=4,116\mu = 4{,}116, σ=1,041\sigma = 1{,}041 для натурального логарифма проницаемости в мД). Мода, медиана и среднее различаются в разы.
Есептеу мысалыЛогнормальная модель проницаемости по десяти образцам
Бастапқы деректер (шартты)
Натуральные логарифмы проницаемости (мД) десяти образцов
среднее m=4,1155m = 4{,}1155, стандартное отклонение s=1,041s = 1{,}041
По самим образцам
среднее арифметическое 100 мД, медиана 60 мД
Шешім
  1. Медиана: e4,1155=61,3\e^{4{,}1155} = 61{,}3 мД.
  2. Среднее: показатель 4,1155+1,0412/2=4,6574{,}1155 + 1{,}041^2/2 = 4{,}657; e4,657=105\e^{4{,}657} = 105 мД.
  3. Мода: показатель 4,1155−1,0412=3,0324{,}1155 - 1{,}041^2 = 3{,}032; e3,032=20,7\e^{3{,}032} = 20{,}7 мД.
  4. 10-й и 90-й процентили: e4,1155−1,282⋅1,041=16,1\e^{4{,}1155 - 1{,}282\cdot1{,}041} = 16{,}1 мД и e4,1155+1,282⋅1,041=233\e^{4{,}1155 + 1{,}282\cdot1{,}041} = 233 мД; их отношение — около 14.
  5. Доля породы с проницаемостью ниже 10 мД: z=(ln⁡10−4,1155)/1,041=−1,74z = (\ln 10 - 4{,}1155)/1{,}041 = -1{,}74, Φ(z)=0,041\Phi(z) = 0{,}041 — около 4 %.
Жауап: Медиана 61,3 мД, среднее 105 мД, мода 20,7 мД; 80 % значений — от 16,1 до 233 мД.

Бұл нені білдіреді. Модель согласуется с выборкой: медиана 61 против 60 мД, среднее 105 против 100 мД. Наиболее вероятное значение (мода, около 21 мД) впятеро меньше среднего: «чаще всего встречается» и «в среднем» для такой породы — совсем разные числа. Десятая часть породы хуже 16 мД, десятая — лучше 230 мД; именно эти «хвосты» определяют застойные зоны и пути прорыва воды. Данные примера условные.

Вероятностный график: проверка распределения

Гистограмма десяти или двадцати образцов слишком груба, чтобы судить о форме распределения. Надёжнее вероятностный график (probability plot): его ось накопленной вероятности растянута так, что нормальное распределение изображается прямой линией. Раньше для этого печатали специальную вероятностную бумагу, теперь шкалу вычисляют.

Как построить
  1. 1
    Упорядочить

    Расположить NN значений по возрастанию и пронумеровать: i=1,2,…,Ni = 1, 2, \dots, N.

  2. 2
    Приписать вероятности

    Каждому значению приписать накопленную вероятность pi=(i−0,5)/Np_i = (i - 0{,}5)/N — середину «своей» доли 1/N1/N. Встречаются и другие формулы, например i/(N+1)i/(N + 1); различие заметно только у крайних точек.

  3. 3
    Найти квантили

    Для каждого pip_i найти ziz_i — квантиль стандартного нормального распределения (по таблице или функцией электронной таблицы).

  4. 4
    Построить

    Отложить значения против ziz_i: сами значения — для проверки на нормальность, значения на логарифмической оси — на логнормальность.

  5. 5
    Прочитать

    Если точки ложатся на прямую, распределение подходит. Значение при z=0z = 0 — медиана; изменение на единицу zz — стандартное отклонение.

Данные для вероятностного графика проницаемости десяти образцов (условные данные)
ii
1
kik_i, мД
12
pip_i
0,05
ziz_i
−1,645
ii
2
kik_i, мД
20
pip_i
0,15
ziz_i
−1,036
ii
3
kik_i, мД
30
pip_i
0,25
ziz_i
−0,674
ii
4
kik_i, мД
40
pip_i
0,35
ziz_i
−0,385
ii
5
kik_i, мД
55
pip_i
0,45
ziz_i
−0,126
ii
6
kik_i, мД
65
pip_i
0,55
ziz_i
0,126
ii
7
kik_i, мД
90
pip_i
0,65
ziz_i
0,385
ii
8
kik_i, мД
120
pip_i
0,75
ziz_i
0,674
ii
9
kik_i, мД
168
pip_i
0,85
ziz_i
1,036
ii
10
kik_i, мД
400
pip_i
0,95
ziz_i
1,645
Есептеу мысалыВероятностный график проницаемости десяти образцов
Бастапқы деректер (шартты)
Проницаемость по возрастанию, мД
12; 20; 30; 40; 55; 65; 90; 120; 168; 400
Среднее и отклонение логарифмов
m=4,1155m = 4{,}1155, s=1,041s = 1{,}041
Шешім
  1. Вероятности pi=(i−0,5)/10p_i = (i - 0{,}5)/10: 0,05; 0,15; 0,25 и так далее до 0,95.
  2. Квантили для первой половины ряда: −1,645; −1,036; −0,674; −0,385; −0,126; для второй половины — те же числа с плюсом в обратном порядке.
  3. Прямая, отвечающая среднему и отклонению логарифмов: ln⁡k=4,1155+1,041 z\ln k = 4{,}1155 + 1{,}041\,z. При z=0z = 0 она даёт 61 мД, при z=−1z = -1 — 22 мД, при z=+1z = +1 — 174 мД.
  4. Сильнее всех от прямой отклоняется самый проницаемый образец: он в 1,2 раза выше неё. Остальные девять точек лежат в пределах 9 % от прямой.
Жауап: Точки ложатся на прямую: логнормальное распределение подходит. Медиана 61 мД; интервал «плюс-минус одно отклонение» — от 22 до 174 мД.

Бұл нені білдіреді. На обычной, не логарифмической оси те же точки легли бы на крутую дугу — признак того, что для самой проницаемости нормальное распределение не годится. Десяти точек мало, чтобы уверенно различать близкие распределения, но достаточно, чтобы отвергнуть заведомо неподходящее. Данные примера условные.

Вероятностный график проницаемости десяти образцов25165084959851020501002005001000Накопленная доля образцов, % (нормальная шкала)Проницаемость, мД22 мД61 мД174 мД

График: по вертикали проницаемость от 5 до 1000 миллидарси в логарифмической шкале, по горизонтали накопленная доля образцов от 2 до 98 процентов в шкале нормального распределения. Десять точек лежат вдоль восходящей прямой; только крайняя правая точка, 400 миллидарси, заметно выше неё. Штриховые линии показывают 22 миллидарси при 16 процентах, 61 миллидарси при 50 процентах и 174 миллидарси при 84 процентах.

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

Отклонения точек от прямой тоже несут информацию.

  • Плавный изгиб по всей длине — масштаб выбран неверно: например, проницаемость отложена без логарифма.
  • Излом, две прямые с разным наклоном — в выборке смешаны две совокупности: две породы, матрица и трещины. Их разделяют и описывают по отдельности.
  • Горизонтальная «полка» на краю — значения обрезаны пределом измерения прибора.
  • Одна-две точки в стороне на краях — возможные выбросы. Прежде чем отбрасывать, проверяют образец и замер; крайние точки вообще наименее надёжны.

Насколько точно известно среднее

Среднее по образцам — оценка. Если бы мы взяли другие десять образцов из того же пласта, среднее получилось бы другим. Вопрос в том, как далеко такие средние разбросаны вокруг истинного.

Ответ даёт центральная предельная теорема (central limit theorem): средние по независимым выборкам объёмом NN распределены приблизительно нормально вокруг истинного среднего μ\mu — тем теснее, чем больше NN, и почти при любом распределении самих данных. Стандартное отклонение таких средних называют стандартной ошибкой (standard error, SE).

Стандартная ошибка среднего
SE=sNSE = \frac{s}{\sqrt{N}}
мұндағы:
  • SESEстандартная ошибка среднего, единицы величины
  • ssстандартное отклонение данных, единицы величины
  • NNчисло независимых измерений, —

Бірліктер: те же, что у измеряемой величины

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

Стандартное отклонение ss описывает разброс данных и с ростом числа образцов не убывает. Стандартная ошибка описывает неопределённость среднего и убывает как 1/N1/\sqrt N: чтобы уменьшить её вдвое, образцов нужно вчетверо больше.

Доверительный интервал среднего
xˉ−t sN  ≤  μ  ≤  xˉ+t sN\bar x - t\,\frac{s}{\sqrt{N}} \;\le\; \mu \;\le\; \bar x + t\,\frac{s}{\sqrt{N}}
мұндағы:
  • μ\muистинное среднее совокупности, единицы величины
  • xˉ\bar xсреднее по выборке, единицы величины
  • ttкоэффициент Стьюдента для выбранной доверительной вероятности и N−1N - 1 степеней свободы, —
  • s, Ns,\ Nстандартное отклонение данных и число измерений, единицы величины, —

Бірліктер: те же, что у измеряемой величины

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

При большом числе образцов и вероятности 95 % коэффициент tt близок к 2 (точнее, к 1,96). При малом числе он больше, потому что само ss оценено неточно.

Коэффициент Стьюдента для доверительной вероятности 95 %
Число образцов NN
5
tt
2,78
Число образцов NN
10
tt
2,26
Число образцов NN
15
tt
2,14
Число образцов NN
20
tt
2,09
Число образцов NN
30
tt
2,05
Число образцов NN
60
tt
2,00
Число образцов NN
очень много
tt
1,96

Двусторонний интервал, N−1N - 1 степеней свободы.

Есептеу мысалыДоверительный интервал средней пористости и средней проницаемости
Бастапқы деректер (шартты)
Пористость, 10 образцов
ϕˉ=19,0\bar\phi = 19{,}0 %, s=3,162s = 3{,}162 %
Проницаемость, те же образцы
kˉ=100\bar k = 100 мД, s=116s = 116 мД
Коэффициент Стьюдента
t=2,262t = 2{,}262 (95 %, 9 степеней свободы)
Шешім
  1. Пористость: SE=3,162/10=1,00SE = 3{,}162/\sqrt{10} = 1{,}00 %. Полуширина интервала: 2,262⋅1,00=2,262{,}262\cdot1{,}00 = 2{,}26 %.
  2. Интервал для средней пористости — от 16,7 до 21,3 %, то есть ±12 % от самого среднего.
  3. Проницаемость: SE=116/10=36,7SE = 116/\sqrt{10} = 36{,}7 мД. Полуширина: 2,262⋅36,7=832{,}262\cdot36{,}7 = 83 мД — это ±83 % от среднего.
Жауап: Средняя пористость 19,0±2,2619{,}0 \pm 2{,}26 %; средняя проницаемость 100±83100 \pm 83 мД.

Бұл нені білдіреді. Одни и те же десять образцов определяют среднюю пористость с точностью около 12 %, а среднюю проницаемость — лишь с точностью ±83 %. Причём для проницаемости формальный интервал (от 17 до 183 мД) ещё и ненадёжен: при таком скошенном распределении среднее по десяти образцам распределено не нормально. Среднюю пористость интервала по этим образцам оценить можно; среднюю проницаемость — практически нет. Данные примера условные.

Сколько образцов достаточно

Формулу доверительного интервала можно обратить: задать желаемую точность и найти число образцов. Точность удобно задавать допуском PP в процентах от среднего: «среднее по образцам отличается от истинного не более чем на PP % с вероятностью 95 %».

Необходимое число образцов для заданной точности среднего
N0=(100 t CvP)2N_0 = \left(\frac{100\,t\,C_v}{P}\right)^2
мұндағы:
  • N0N_0необходимое число образцов; результат округляют вверх, —
  • ttкоэффициент Стьюдента; около 2 при вероятности 95 % и числе образцов больше тридцати, —
  • CvC_vкоэффициент вариации s/xˉs/\bar x, —
  • PPдопуск — допустимое отклонение среднего по образцам от истинного, % от среднего

Бірліктер: PP — в процентах; остальные величины безразмерны

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

Формула получается из условия, что полуширина доверительного интервала t s/N0t\,s/\sqrt{N_0} равна допуску P xˉ/100P\,\bar x/100. Число образцов растёт как квадрат коэффициента вариации и как квадрат обратного допуска: вдвое более неоднородная порода или вдвое более жёсткий допуск требуют вчетверо больше образцов.

Правило для числа образцов: точность среднего ±20 % с вероятностью 95 %
N0=(10 Cv)2,Ps=200 CvNsN_0 = (10\,C_v)^2, \qquad P_s = \frac{200\,C_v}{\sqrt{N_s}}
мұндағы:
  • N0N_0число образцов для допуска ±20 %, —
  • CvC_vкоэффициент вариации, —
  • NsN_sфактическое число образцов, —
  • PsP_sточность среднего, достигнутая при NsN_s образцах, %

Бірліктер: PsP_s — в процентах от среднего

Мына жағдайда орындалады
  • t=2t = 2: вероятность 95 %, число образцов больше тридцати;
  • первая формула относится к допуску 20 %, вторая — к любому числу образцов при t=2t = 2;
  • для очень неоднородных пород (CvC_v заметно больше 1, например сложные карбонаты) допуск ±20 % часто недостижим; тогда принимают ±50 %, и N0=(4 Cv)2N_0 = (4\,C_v)^2.

Правило позволяет за минуту ответить на два вопроса: сколько измерений нужно и чего стоят имеющиеся. Оно получено и проверено главным образом для проницаемости.

Число образцов для разных пород и допусков
CvC_v
0,2
±50 %
1
±20 %
4
±10 %
16
CvC_v
0,5
±50 %
4
±20 %
25
±10 %
100
CvC_v
1,0
±50 %
16
±20 %
100
±10 %
400
CvC_v
1,5
±50 %
36
±20 %
225
±10 %
900
CvC_v
2,0
±50 %
64
±20 %
400
±10 %
1600
CvC_v
3,0
±50 %
144
±20 %
900
±10 %
3600

Рассчитано при t=2t = 2. Числа меньше 30 занижены: при малых выборках tt больше двух, а оценить сам разброс по двум-трём образцам невозможно.

Точность среднего в зависимости от числа образцов при разных коэффициентах вариации410251004001000020406080100Число образцов NТочность среднего, ± %Cv = 0,2Cv = 0,5Cv = 1,0Cv = 2,0

График: по горизонтали число образцов от 3 до 1000 в логарифмической шкале, по вертикали точность среднего от 0 до 100 процентов. Четыре убывающие кривые для коэффициентов вариации 0,2; 0,5; 1,0 и 2,0. Горизонтальная штриховая линия на уровне 20 процентов пересекает их при 4, 25, 100 и 400 образцах.

ЕсептеуТочность среднего Ps=200 Cv/NsP_s = 200\,C_v/\sqrt{N_s} в зависимости от числа образцов при четырёх значениях коэффициента вариации. Штриховая линия — допуск ±20 %; точки на ней — необходимое число образцов N0=(10 Cv)2N_0 = (10\,C_v)^2.
Есептеу мысалыХватает ли десяти образцов
Бастапқы деректер (шартты)
Интервал
3,0 м, 10 образцов через 0,3 м
Коэффициенты вариации
пористость 0,166; проницаемость 1,16
Шешім
  1. Пористость: N0=(10⋅0,166)2=2,8N_0 = (10\cdot0{,}166)^2 = 2{,}8 — формально хватило бы 3 образцов. Десяти достаточно с запасом; они дают точность 200⋅0,166/10=10,5200\cdot0{,}166/\sqrt{10} = 10{,}5 %.
  2. Проницаемость: N0=(10⋅1,16)2=134,6N_0 = (10\cdot1{,}16)^2 = 134{,}6, то есть 135 образцов. Десять образцов дают 200⋅1,16/10=73200\cdot1{,}16/\sqrt{10} = 73 %.
  3. Шаг отбора для 135 измерений в интервале 3,0 м: 3,0/135=0,0223{,}0/135 = 0{,}022 м, или 2,2 см.
Жауап: Для средней пористости десяти образцов хватает (±10,5 %). Для средней проницаемости с допуском ±20 % нужно около 135 измерений — в тринадцать с лишним раз больше.

Бұл нені білдіреді. Шаг в 2 см меньше диаметра самого цилиндрического образца: столько цилиндров из трёх метров керна не выбурить. Плотные измерения проницаемости делают точечным пермеаметром (probe permeameter; его называют также профильным или мини-пермеаметром): наконечник прижимают к поверхности распиленного керна и по расходу газа определяют проницаемость в точке. Замеры быстрые и не разрушают керн, поэтому сотня точек на интервал — посильная программа. Точности ±10,5 и ±73 % получены по правилу с t=2t = 2; с точным t=2,26t = 2{,}26 для десяти образцов они хуже — ±12 и ±83 %, как в предыдущем примере. Данные примера условные.

Как планировать отбор
  1. 1
    Пробная серия

    Сделать около 25 равномерно распределённых измерений в изучаемом интервале — пласте, пачке, отдельной породе.

  2. 2
    Оценить разброс

    Вычислить CvC_v. Если он не больше 0,5, двадцати пяти измерений уже достаточно: (10⋅0,5)2=25(10\cdot0{,}5)^2 = 25.

  3. 3
    Сгустить сеть

    Если CvC_v больше, поставить новые точки между уже сделанными — столько, чтобы общее число измерений достигло N0N_0. Шаг равномерной сети равен длине интервала, делённой на N0N_0.

  4. 4
    Проверить

    Пересчитать CvC_v и достигнутую точность PsP_s. Если нужное число измерений недостижимо, прямо указать более широкий допуск.

Смещение: когда больше образцов не помогает

Всё сказанное о точности относилось к случайной ошибке: она убывает с ростом числа образцов. Но есть и другая ошибка — смещение (bias), систематическая ошибка отбора. Если образцы взяты не из всей совокупности, а из её удобной части, их среднее сходится к неверному значению — тем увереннее, чем образцов больше.

В данных керна смещение накапливается на каждом шаге — от выбора скважины до измерения в лаборатории.

  • В каких скважинах отбирают керн. Это разведочные и первые эксплуатационные скважины, обычно в сводовой части залежи. На крыльях и в водоносной зоне керна почти нет, хотя порода там может быть иной: вторичные изменения в нефтенасыщенной и водонасыщенной частях пласта могут идти по-разному.
  • Что удаётся вынести. Рыхлые, трещиноватые и кавернозные интервалы разрушаются при бурении и подъёме; вынос керна из лучших коллекторов нередко самый низкий.
  • Из чего выбуривают образцы. Цилиндр получается из крепкой однородной породы. Рыхлый песок рассыпается — теряются самые высокие значения проницаемости; тонкие плотные прослои обходят — теряются самые низкие.
  • Где сгущают отбор. Продуктивные песчаники опробуют гуще, чем алевролиты и глины. Простое среднее по образцам тогда описывает коллектор, а не разрез в целом.
  • Условия измерения. Образцы очищены, высушены и разгружены от горного давления, поэтому пористость и особенно проницаемость по керну обычно несколько выше пластовых.
  • Пределы прибора. Значения ниже и выше диапазона измерения обрезаны.
Есептеу мысалыГусто в хорошей породе, пусто в лучшей
Бастапқы деректер (шартты)
Крепкий чистый песчаник
40 % толщины пласта; 24 образца; средняя проницаемость 300 мД
Глинистый песчаник
45 % толщины; 6 образцов; 40 мД
Рыхлый крупнозернистый песок
15 % толщины; образцы выбурить не удалось; по другим данным около 2000 мД
Шешім
  1. Простое среднее по 30 образцам: (24⋅300+6⋅40)/30=248(24\cdot300 + 6\cdot40)/30 = 248 мД.
  2. Среднее по тем же двум породам, но с весом их толщин: (0,40⋅300+0,45⋅40)/0,85=162(0{,}40\cdot300 + 0{,}45\cdot40)/0{,}85 = 162 мД. Простое среднее завышено в 1,5 раза, потому что хороший песчаник опробован гуще.
  3. Среднее по всему пласту, включая рыхлый песок: 0,40⋅300+0,45⋅40+0,15⋅2000=4380{,}40\cdot300 + 0{,}45\cdot40 + 0{,}15\cdot2000 = 438 мД.
  4. Если удвоить число образцов при том же способе отбора (48 и 12), простое среднее останется равным 248 мД; уменьшится только его стандартная ошибка.
Жауап: 248 мД по образцам; 162 мД для опробованной части разреза; около 438 мД для пласта в целом.

Бұл нені білдіреді. Смещений два, и направлены они в разные стороны. Первое — неравномерную густоту отбора — исправляют весами: средние по породам умножают на их доли в разрезе, известные по каротажу. Второе — породу, которая в образцы не попала вовсе, — никакими весами не исправить, а здесь 15 % толщины дают 68 % проводимости вдоль напластования. Нужны другие данные: керн, отобранный с сохранением рыхлых интервалов, боковой керн, каротаж, гидродинамические исследования скважин (ГДИС). Данные примера условные.

Доли пород в разрезе пласта и в образцах кернаКрепкий чистый песчаник300 мДв разрезе40 %в образцах80 %Глинистый песчаник40 мДв разрезе45 %в образцах20 %Рыхлый крупнозернистый песок2000 мДв разрезе15 %в образцах0 %

Три пары горизонтальных столбиков. Крепкий чистый песчаник, 300 миллидарси: 40 процентов разреза и 80 процентов образцов. Глинистый песчаник, 40 миллидарси: 45 процентов разреза и 20 процентов образцов. Рыхлый крупнозернистый песок, 2000 миллидарси: 15 процентов разреза и ни одного образца.

ЕсептеуДоли трёх пород примера в толщине пласта (по каротажу) и в общем числе образцов керна. Справа — средняя проницаемость породы.

Полностью избавиться от смещения нельзя, но его можно уменьшить и, главное, вовремя заметить.

  • Описывать весь керн, считать средние по породам и взвешивать их по долям пород в разрезе, взятым с каротажа.
  • Отмечать интервалы без выноса керна и без образцов и выяснять, чем они сложены.
  • Сопоставлять керн с каротажем и с результатами ГДИС: систематическое расхождение — признак смещения.
  • Планировать отбор заранее: равномерная сеть и дополнительные точки в изменчивых прослоях, а не «где удобно выбурить».

Терминдер мен белгілер

Плотность вероятности f(x)f(x)
кривая, площадь под которой между двумя значениями равна доле совокупности с такими значениями.
Функция распределения F(x)F(x)
вероятность того, что величина не превысит xx; растёт от 0 до 1.
μ\mu, σ\sigma
среднее и стандартное отклонение совокупности; для логнормального распределения — логарифма величины.
zz, Φ(z)\Phi(z)
стандартизованное значение и накопленная вероятность стандартного нормального распределения.
Вероятностный график
график значений против квантилей нормального распределения; нормальные данные ложатся на прямую.
SESE
стандартная ошибка среднего, s/Ns/\sqrt N: мера неопределённости среднего, а не разброса данных.
tt
коэффициент Стьюдента: число стандартных ошибок в полуширине доверительного интервала; около 2 при вероятности 95 % и большом числе образцов.
N0N_0, PsP_s
число образцов, необходимое для заданного допуска, и точность среднего, достигнутая при имеющемся числе образцов.
Смещение
систематическое отличие выборки от совокупности; не уменьшается с ростом числа образцов.
VDPV_{DP}
коэффициент Дикстры–Парсонса: мера неоднородности проницаемости от 0 до 1, определяемая по вероятностному графику.

Жиі қателіктер

Қате

«Если взять больше образцов, стандартное отклонение уменьшится».

Шын мәнінде

Стандартное отклонение — свойство самих данных: с ростом выборки оно уточняется, но не убывает. Убывает стандартная ошибка среднего, s/Ns/\sqrt N.

Қате

«95-процентный доверительный интервал содержит 95 % образцов».

Шын мәнінде

Он относится к среднему, а не к отдельным значениям: это диапазон, который с вероятностью 95 % накрывает истинное среднее. 95 % отдельных значений нормального распределения занимают гораздо более широкий интервал xˉ±2s\bar x \pm 2s.

Қате

«Образцов много — значит, среднее надёжно».

Шын мәнінде

Число образцов уменьшает только случайную ошибку. Если отбор смещён — керн только из свода, образцы только из крепкой породы, — среднее точно сходится к неверному значению.

Қате

«Проницаемость распределена нормально, как всё в природе».

Шын мәнінде

Приблизительно нормально распределён её логарифм. В миллидарси распределение резко скошено: мода, медиана и среднее различаются в разы.

Қате

«P50 — это ожидаемое, то есть среднее, значение».

Шын мәнінде

P50 — медиана. У скошенного вправо распределения среднее больше медианы: у логнормального — в eσ2/2\e^{\sigma^2/2} раз. Для запасов это означает, что ожидаемое значение выше оценки P50.

Толық: тереңдетілген

Коэффициент Дикстры–Парсонса

В нефтяной литературе неоднородность проницаемости часто выражают коэффициентом Дикстры–Парсонса (Dykstra–Parsons coefficient; встречается и написание «Дайкстра»), предложенным в середине XX века для расчётов заводнения слоистых пластов. Он опирается на логнормальное распределение и читается прямо с вероятностного графика.

Коэффициент Дикстры–Парсонса
VDP=k50−k84,1k50=1−e−σln⁡kV_{DP} = \frac{k_{50} - k_{84{,}1}}{k_{50}} = 1 - \e^{-\sigma_{\ln k}}
мұндағы:
  • VDPV_{DP}коэффициент Дикстры–Парсонса, —
  • k50k_{50}медиана проницаемости, мД
  • k84,1k_{84{,}1}проницаемость, которую превышают 84,1 % образцов, — на одно стандартное отклонение логарифма ниже медианы, мД
  • σln⁡k\sigma_{\ln k}стандартное отклонение натурального логарифма проницаемости, —

Бірліктер: безразмерный, от 0 до 1

Мына жағдайда орындалады
  • проницаемость распределена логнормально; для других распределений коэффициент может вводить в заблуждение;
  • k50k_{50} и k84,1k_{84{,}1} берут по прямой, проведённой через точки вероятностного графика, а не по отдельным образцам;
  • процент 84,1 отсчитан как доля образцов с большей проницаемостью — это вероятность превышения, а не накопленная доля. На графике с осью накопленной доли, как в этом уроке, k84,1k_{84{,}1} читают у отметки 16 % (z=−1z = -1).

VDP=0V_{DP} = 0 — однородный пласт; чем ближе коэффициент к 1, тем сильнее неоднородность. Вторая форма следует из того, что k84,1=k50 e−σk_{84{,}1} = k_{50}\,\e^{-\sigma}. Через десятичные логарифмы: VDP=1−10−σlg⁡kV_{DP} = 1 - 10^{-\sigma_{\lg k}}.

Есептеу мысалыКоэффициент Дикстры–Парсонса по десяти образцам
Бастапқы деректер (шартты)
Прямая на вероятностном графике
ln⁡k=4,1155+1,041 z\ln k = 4{,}1155 + 1{,}041\,z
Шешім
  1. Медиана: k50=e4,1155=61,3k_{50} = \e^{4{,}1155} = 61{,}3 мД. На одно отклонение ниже: k84,1=e4,1155−1,041=21,6k_{84{,}1} = \e^{4{,}1155 - 1{,}041} = 21{,}6 мД.
  2. VDP=(61,3−21,6)/61,3=0,65V_{DP} = (61{,}3 - 21{,}6)/61{,}3 = 0{,}65. Проверка: 1−e−1,041=0,651 - \e^{-1{,}041} = 0{,}65.
  3. Коэффициент вариации той же логнормальной модели: e1,0412−1=1,40\sqrt{\e^{1{,}041^2} - 1} = 1{,}40. По самим десяти образцам получалось 1,16: в малых выборках из скошенных распределений коэффициент вариации обычно занижен.
Жауап: VDP≈0,65V_{DP} \approx 0{,}65.

Бұл нені білдіреді. У большинства коллекторов коэффициент лежит между 0,5 и 0,9, так что наш пласт — в умеренной части диапазона. Шкала сжата у верхнего края: значениям VDPV_{DP} 0,6; 0,8 и 0,9 отвечают σln⁡k\sigma_{\ln k} около 0,9; 1,6 и 2,3, то есть пласты различаются гораздо сильнее, чем кажется по самим числам. Как коэффициент используют в расчётах заводнения — в уроке «Неоднородность пласта и осреднение проницаемости». Данные примера условные.

Меры неоднородности логнормального распределения
σln⁡k\sigma_{\ln k}
0,5
Среднее / медиана
1,13
CvC_v
0,53
VDPV_{DP}
0,39
x90/x10x_{90}/x_{10}
3,6
σln⁡k\sigma_{\ln k}
1,0
Среднее / медиана
1,65
CvC_v
1,31
VDPV_{DP}
0,63
x90/x10x_{90}/x_{10}
13
σln⁡k\sigma_{\ln k}
1,5
Среднее / медиана
3,08
CvC_v
2,91
VDPV_{DP}
0,78
x90/x10x_{90}/x_{10}
47
σln⁡k\sigma_{\ln k}
2,0
Среднее / медиана
7,39
CvC_v
7,32
VDPV_{DP}
0,86
x90/x10x_{90}/x_{10}
168

Здесь x10x_{10} и x90x_{90} — 10-й и 90-й процентили по накопленной доле. Все четыре меры однозначно связаны с σln⁡k\sigma_{\ln k}, поэтому для логнормальных данных достаточно знать любую из них.

Доверительный интервал для логнормальных данных

Для сильно скошенных данных интервал строят в логарифмах, где распределение близко к нормальному. Получается интервал для среднего логарифма, то есть для геометрического среднего (медианы).

Доверительный интервал геометрического среднего
exp⁡(m−t sln⁡N)  ≤  x50  ≤  exp⁡(m+t sln⁡N)\exp\left(m - t\,\frac{s_{\ln}}{\sqrt N}\right) \;\le\; x_{50} \;\le\; \exp\left(m + t\,\frac{s_{\ln}}{\sqrt N}\right)
мұндағы:
  • x50x_{50}медиана (геометрическое среднее) совокупности, единицы величины
  • mmсреднее натуральных логарифмов выборки, —
  • sln⁡s_{\ln}стандартное отклонение натуральных логарифмов, —
  • t, Nt,\ Nкоэффициент Стьюдента и число измерений, —

Бірліктер: границы — в единицах величины

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

В обычном масштабе интервал несимметричен, а «в разах» симметричен: оценку умножают и делят на один и тот же множитель exp⁡(t sln⁡/N)\exp(t\,s_{\ln}/\sqrt N).

Есептеу мысалыИнтервал для геометрического среднего проницаемости
Бастапқы деректер (шартты)
Логарифмы проницаемости десяти образцов
m=4,1155m = 4{,}1155, sln⁡=1,041s_{\ln} = 1{,}041
Коэффициент Стьюдента
t=2,262t = 2{,}262
Шешім
  1. Полуширина в логарифмах: 2,262⋅1,041/10=0,7452{,}262\cdot1{,}041/\sqrt{10} = 0{,}745. Множитель: e0,745=2,11\e^{0{,}745} = 2{,}11.
  2. Интервал: от 61,3/2,11=2961{,}3/2{,}11 = 29 мД до 61,3⋅2,11=12961{,}3\cdot2{,}11 = 129 мД.
Жауап: С вероятностью 95 % геометрическое среднее лежит между 29 и 129 мД.

Бұл нені білдіреді. Геометрическое среднее по десяти образцам известно с точностью до множителя 2,1. Для среднего арифметического логнормальных данных интервал ещё шире и тоже несимметричен: к неопределённости медианы добавляется неопределённость разброса, от которого среднее зависит через eσ2/2\e^{\sigma^2/2}. Вывод тот же, что и в основном тексте: десятью образцами среднюю проницаемость неоднородного пласта не определить. Данные примера условные.

Зависимые образцы

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

Стандартная ошибка среднего при зависимых измерениях
SE=sN 1+(N−1) ρˉ,Nэф=N1+(N−1) ρˉSE = \frac{s}{\sqrt{N}}\,\sqrt{1 + (N - 1)\,\bar\rho}, \qquad N_{\text{эф}} = \frac{N}{1 + (N - 1)\,\bar\rho}
мұндағы:
  • SESEстандартная ошибка среднего, единицы величины
  • ρˉ\bar\rhoсредний коэффициент корреляции между парами измерений: 0 — независимы, 1 — повторяют друг друга, —
  • NэфN_{\text{эф}}эффективное число независимых измерений, —
  • s, Ns,\ Nстандартное отклонение данных и число измерений, единицы величины, —

Бірліктер: SESE — в единицах измеряемой величины

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

При ρˉ=0\bar\rho = 0 формула переходит в обычную. При положительной корреляции стандартная ошибка с ростом NN не стремится к нулю, а останавливается на уровне sρˉs\sqrt{\bar\rho}: сколько бы образцов ни взяли в одной скважине, о пласте в целом они говорят ограниченно.

Есептеу мысалыЧего стоят десять похожих образцов
Бастапқы деректер (шартты)
Число образцов
10
Средняя корреляция между образцами ρˉ\bar\rho
0,3
Шешім
  1. Эффективное число независимых измерений: 10/(1+9⋅0,3)=2,710/(1 + 9\cdot0{,}3) = 2{,}7.
  2. Стандартная ошибка больше, чем по формуле для независимых образцов, в 1+9⋅0,3=1,9\sqrt{1 + 9\cdot0{,}3} = 1{,}9 раза.
Жауап: Десять образцов со средней корреляцией 0,3 равноценны 2,7 независимого образца.

Бұл нені білдіреді. Как быстро сходство убывает с расстоянием, описывает вариограмма (Пространственная корреляция и вариограмма). Практический вывод простой: для суждения о пласте в целом число скважин с керном важнее числа образцов в каждой из них. Данные примера условные.

Басқа тақырыптармен байланыс

Өзіңізді тексеріңіз

  1. 1

    Пористость пласта распределена нормально со средним 20 % и стандартным отклонением 3 %. Какая доля породы имеет пористость от 17 до 23 %?

    1. А.Около 50 %
    2. Б.Около 68 %
    3. В.Около 95 %
    4. Г.Около 99,7 %
    Жауабын көрсету

    Дұрыс жауап: Б. Около 68 %

    Интервал от 17 до 23 % — это среднее плюс-минус одно стандартное отклонение, а в пределах μ±σ\mu \pm \sigma лежит 68 % значений нормального распределения. 95 % лежит в пределах двух отклонений — от 14 до 26 %.

  2. 2

    По 25 образцам коэффициент вариации проницаемости равен 0,8. Сколько образцов нужно для точности ±20 % и какую точность дают имеющиеся?

    Жауабын көрсету

    Нужно N0=(10⋅0,8)2=64N_0 = (10\cdot0{,}8)^2 = 64 образца. Имеющиеся 25 дают Ps=200⋅0,8/25=32P_s = 200\cdot0{,}8/\sqrt{25} = 32 %: среднее известно с точностью примерно до трети. Чтобы дойти до ±20 %, сеть нужно сгустить в два с половиной раза.

  3. 3

    На вероятностном графике логарифмов проницаемости точки образуют две прямые с изломом между ними. Что это означает?

    1. А.Образцов слишком много
    2. Б.В выборке смешаны две породы с разными распределениями
    3. В.Нужно вернуться от логарифмов к самим значениям
    4. Г.Распределение нормальное
    Жауабын көрсету

    Дұрыс жауап: Б. В выборке смешаны две породы с разными распределениями

    Одна прямая — одна совокупность с одним средним и одним разбросом. Излом означает смесь двух совокупностей, например двух пород или матрицы и трещин; их нужно разделить и описывать по отдельности. Число образцов на форму графика не влияет.

  4. 4

    Логарифм проницаемости распределён нормально: μ=ln⁡50\mu = \ln 50, σ=1,2\sigma = 1{,}2 (натуральные логарифмы, проницаемость в мД). Найдите медиану и среднее арифметическое.

    Жауабын көрсету

    Медиана: eμ=50\e^{\mu} = 50 мД. Среднее: eμ+σ2/2=50⋅e0,72≈103\e^{\mu + \sigma^2/2} = 50\cdot\e^{0{,}72} \approx 103 мД — вдвое больше медианы. Чем шире разброс, тем дальше среднее уходит от «типичного» значения.

  5. 5

    Керн отобран только в сводовых скважинах, образцы — только из крепкого песчаника. Число образцов увеличили со 100 до 400. Что изменилось?

    1. А.Стандартная ошибка уменьшилась вдвое, смещение осталось прежним
    2. Б.Ничего не изменилось
    3. В.Смещение уменьшилось вдвое
    4. Г.И стандартная ошибка, и смещение уменьшились вдвое
    Жауабын көрсету

    Дұрыс жауап: А. Стандартная ошибка уменьшилась вдвое, смещение осталось прежним

    Стандартная ошибка убывает как 1/N1/\sqrt N: вчетверо больше образцов — вдвое меньше ошибка. Но это точность оценки среднего опробованной породы. Рыхлые прослои и крылья залежи в выборку по-прежнему не попали, и смещение осталось тем же.

  6. 6

    Средняя пористость по 16 образцам — 18 %, стандартное отклонение — 2 %. Найдите стандартную ошибку и 95-процентный доверительный интервал среднего (t=2,13t = 2{,}13).

    Жауабын көрсету

    SE=2/16=0,5SE = 2/\sqrt{16} = 0{,}5 %. Полуширина интервала: 2,13⋅0,5≈1,12{,}13\cdot0{,}5 \approx 1{,}1 %. Истинная средняя пористость с вероятностью 95 % лежит между 16,9 и 19,1 %. Отдельные образцы при этом разбросаны гораздо шире — примерно от 14 до 22 %.

Жаңартылған: 2 октября 2026 г.