УДК 553.98:532.311.8 |
Определение аномально высоких поровых давлений в глинистых породах методом компрессионной кривой
В.М. ДОБРЫНИН, В.А. СЕРЕБРЯКОВ, А.Д. СРЕБРОДОЛЬСКИЙ (МИНХиГП)
В статье описана новая методика определения аномально высоких поровых давлений (АВПоД) в глинистых породах. Методика основана на использовании свойств компрессионной кривой - закономерности, характеризующей уплотнения глинистых пород в зависимости от разности между давлением массы вышележащих осадочных пород (горным давлением) и давлением насыщающих глины флюидов (s - р).
Используемые в настоящее время методики оценки АВПоД основаны на изучении графических зависимостей физических свойств глинистых пород от глубины их залегания. Эти зависимости закономерны в зоне нормально уплотненных глинистых пород. В зоне АВПоД они нарушаются, что и является качественным и количественным признаком аномально поровых давлений.
Компрессионные кривые более универсальны. Они едины в зоне залегания глинистых пород с нормальным и аномальным давлением поровой жидкости. Это свойство предоставляет широкие возможности для аналитического определения параметров компрессионной кривой и использования вычислительной техники. Кроме того, нет необходимости прибегать к трудоемким графическим построениям.
Изучение уплотнения глинистых пород показывает, что при неизменных свойствах насыщающего флюида и неизменном минералогическом составе глин между различными физическими свойствами глинистых пород (х) и разностью давлений (s - р) имеют место экспоненциальные зависимости, которые в полулогарифмическом масштабе становятся прямыми (см. рисунок ). В природе этот факт находит подтверждение в существовании широко известных экспоненциальных зависимостей между пористостью, плотностью или удельным электрическим сопротивлением нормально уплотненных глинистых пород и глубиной их залегания; в экспериментальных исследованиях грунтоведов при уплотнении глинистых пород также отмечаются экспоненциальные формы компрессионных кривых. Более широкое теоретическое толкование этому явлению дано в работе [1].
Отклонения от прямой линии могут возникнуть в случае непостоянства минерализации погребенных в глинах поровых вод, а также в случае неравномерного увеличения температуры глин с глубиной. Первое явление наблюдается обычно в глинах, залегающих на глубинах менее 800-1000 м, второе определяется путем аналитического приведения всей кривой к одной и той же температуре, как это описано в работе [1].
Выделим на прямой, отображенной на рисунке , точки 1, 2 и 3. Все измерения параметра X на прямой приведем к температуре точки 1. В точках 1 и 2 давление поровых флюидов в глинистых поводах известно и равно p1 и р2. В точке 3 это давление требуется определить (Pa).
Запишем с учетом сказанного координаты точек 1, 2, 3:
где a(x) - значение температурного коэффициента для физического свойства х [1]; Г - среднее значение геотермического градиента в интервале h1-h2; h, h1, h2 - глубины точки соответственно 3, 1 и 2; s, s1, s2 - средние нормальные напряжения соответственно на глубинах h, h1 и h2.
Уравнение прямой, проходящей через точки 1 и 2, имеет вид
Это уравнение можно привести к более удобному для вычисления уравнению прямой - уравнению с угловым коэффициентом:
где введены следующие обозначения для параметров этой прямой:
Таким образом, для определения аномального порового давления в глинистых породах по методу компрессионной кривой необходимо знать значения параметров b (х) и К(х), а также возможное их изменение по площади изучаемого региона.
В районах с относительно простым геологическим строением, где имеются мощные пачки нормально уплотненных чистых глинистых пород, коэффициенты b(х) и К(х) могут быть определены по формулам (3) с использованием геофизических диаграмм и кривой изменения плотности глинистых пород с глубиной. В районах с более сложным строением, включающим подсолевые отложения, коэффициенты b(х) и К(х) и их изменения определяются путем использования геофизических материалов и прямых измерений пластовых давлений (нормальных или аномальных) в коллекторах, вмещаемых глинистыми породами.
Этот второй прием основывается на подтвержденном положении о том, что пластовое давление в коллекторах равно поровому давлению во вмещающих эти коллекторы глинистых породах. В работе [2] это положение было подтверждено результатами исследования АВПД в скважинах девяти нефтегазоносных районов Советского Союза. Исключение в данных исследованиях составила Азербайджанская нефтегазоносная область, где в ряде нефтяных месторождений пластовые давления в коллекторах оказались существенно ниже, чем во вмещающих глинах.
В таблице приведены результаты контрольных определений поровых давлений по предлагаемой методике в 17 скважинах месторождений различных нефтегазоносных областей, в которых ранее была успешно применена известная методика “кривых нормально уплотненных глин”. Другими словами, для определения параметров b(х) и К(х), необходимых для оценки АВПоД по методу компрессионной кривой, использовались ранее полученные зависимости измеренных геофизических параметров от глубины залегания пластов. Как и ожидалось, по обеим методикам в одних и тех же изучаемых интервалах получены практически не отличающиеся результаты. Среднее отклонение значений аномальных поровых давлений от данных прямых замеров составило ±5,58 % (см. таблицу ).
Таким образом, использование методики компрессионной кривой дает надежные результаты, которые хорошо согласуются с данными других методик определения АВПоД и с данными, полученными прямыми измерениями давлений в коллекторах глубинными манометрами.
Методику компрессионной кривой целесообразно использовать в тех районах, где трудно установить положение “кривой нормально уплотненных глин”. Однако при этом следует учитывать, что при нахождении параметров b(х) и К(х) в одиночных скважинах по ограниченному числу известных значений s и р точность определения поровых давлений в связи со случайными погрешностями может существенно снизиться. Поэтому рекомендуется применять статистический подход по изучаемым или даже по соседним скважинам, находящимся в сходных геологических условиях. Например, коэффициенты b(х) и К(х), определенные в четырех скважинах месторождения Кенкияк, изменяются: b(х) от -0,26 до -0,53, К(х) от 0,0032 до 0,0036. Среднее арифметическое значение этих параметров составляет: К(х) = 0,0034, b(х) = 0,38.
Значения АВПоД, рассчитанные по средним значениям параметров К(х) и b(х), в данном конкретном случае отличаются от значений поровых давлений, определенных по параметрам каждой конкретной скважины, не более чем на 7 %. Таким образом, можно говорить о возможности использования средних значений параметров К(х) и b(х) для прогнозирования АВПД на конкретной площади или на территории со сходным геологическим строением.
Рассмотрим более сложный случай, когда значения физического свойства глинистых пород х в разрезе какого-либо месторождения или целого нефтегазового района не могут быть оценены в условиях нормального уплотнения (например, в результате влияния искажающих факторов). В этом случае общепринятые методики не могут быть использованы для прогнозирования АВПД. В то же время предлагаемый метод компрессионной кривой может найти применение и в данном случае, если имеются результаты двух замеров пластового давления в одной или нескольких имеющихся скважинах. Эти замеры давлений совместно с данными о физических свойствах глин в близлежащих интервалах позволяют с помощью уравнений (3) оценить необходимые параметры К(х) и b(х) и с их помощью провести расчет АВПоД по данным промысловой геофизики в любой соседней скважине этого региона.
Подобный пример рассмотрен нами также на месторождении Кенкияк, где в качестве исходных материалов были приняты два замера пластового давления в скв. Г-91 на глубине 4230 м (ра = 75,6 МПа, rп.гл =3.5 Ом-м) и на глубине 4330 м (Ра = 67.3 МПа, rп.гл=5,5 Ом-м). По этим данным были рассчитаны коэффициенты К(х) = 0.00183 и b(х) = 0,0374. которые в дальнейшем использовались для оценки АВПоД (43 определения) в других восьми скважинах этого месторождения. Сравнение полученных значений давлений с результатами пяти манометрических замеров в этих скважинах показало хорошую сходимость результатов, расхождение ±8 %.
Таким образом, исследования, выполненные по опробованию предлагаемой методики "компрессионной кривой”, свидетельствуют о возможности ее применения для прогнозирования АВПД в различных геологических условиях, в том числе в районах со сложным геологическим разрезом, где из-за отсутствия достаточно мощной толщи глин в верхних интервалах разреза, наличия мощной толщи соли или значительного влияния других искажающих факторов трудно установить закономерность изменения физических свойств глинистых пород, залегающих в условиях нормальных давлений.
Вышеизложенное позволяет сделать следующие выводы.
Предлагается новая методика прогнозирования АВПоД, основанная на использовании уравнения “компрессионной кривой” - зависимости изменения различных петрофизических свойств чистых глинистых пород от разности между давлением массы вышележащих пород и поровым давлением насыщающего флюида.
Методика “компрессионной кривой” дает возможность оценивать АВПоД в районах со сложным геологическим строением, а также значительно упрощает применение ЭВМ для прогнозирования АВПД по сравнению с известными методиками.
СПИСОК ЛИТЕРАТУРЫ
Поступила 13/XI 1981 с.
Площадь |
Скважина |
Глубина кровли коллектора, м |
Поровое давление, определенное по методике компрессионной кривой и приведенное к глубине коллектора, Рпор, МПа |
К(х) |
b(x) |
Пластовое давление, измеренное в коллекторе, р, МПа |
Относительное расхождение определения р и Рпор, % |
Кенкияк |
Г-86 |
4080 |
67,4 |
0,0032 |
-0,26 |
75,2 |
-10,4 |
Г-86 |
4150 |
69,1 |
0,0034 |
-0,36 |
75,2 |
-8,1 |
|
Г-90 |
3950 |
69,0 |
0,0035 |
-0,36 |
61,9 |
+11,1 |
|
Г-90 |
4110 |
77,8 |
0,0035 |
-0,33 |
74,0 |
+5,1 |
|
Г-91 |
4210 |
72,9 |
0,0034 |
-0,38 |
75,6 |
-3,6 |
|
Г-91 |
4230 |
73,3 |
0,0034 |
-0,40 |
75,6 |
-3,1 |
|
Г-91 |
4330 |
66,8 |
0,0036 |
-0,49 |
67,3 |
-0,9 |
|
Г-92 |
4570 |
71,9 |
0,0036 |
-0,53 |
78,5 |
-8,4 |
|
Каратюбе |
Г-25 |
4450 |
80,6 |
0,0032 |
-0,40 |
82,6 |
-2.4 |
Г-33 |
4440 |
78,0 |
0,0037 |
-0,68 |
80,5 |
-3.2 |
|
Г-34 |
4330 |
77,1 |
0,0038 |
-0,69 |
79,5 |
-3.0 |
|
Г-34 |
4560 |
80,7 |
0,0037 |
-0,68 |
81,1 |
-0,5 |
|
Г-35 |
4540 |
77,4 |
0,0037 |
-0,68 |
81,1 |
-4,6 |
|
Фонтановская |
Ф-5 |
3230 |
45,3 |
0,0031 |
-0,62 |
52,9 |
-14,5 |
Ф-7 |
3320 |
43,9 |
0,0031 |
-0,71 |
47,0 |
-6,7 |
|
Ф-8 |
3150 |
45,0 |
0,0036 |
-0,50 |
51,3 |
-12,3 |
|
Новодмитриевская |
59 |
2600 |
41,8 |
0,0028 |
-0,73 |
41,7 |
+1,5 |
275 |
2500 |
32,2 |
0,0027 |
-0,72 |
34,0 |
-5,4 |
|
505 |
2410 |
39,3 |
0,0027 |
-0,71 |
38,5 |
+2..2 |
|
510 |
2700 |
40,1 |
0,0028 |
-0,67 |
42,5 |
-5,6 |
|
Горомай |
20 |
2350 |
31,6 |
0,0028 |
-0,25 |
33,3 |
-5,1 |
20 |
2660 |
37,7 |
0,0026 |
-0,54 |
36,8 |
+2,4 |
|
Паромай |
202 |
850 |
13,4 |
0,0027 |
-0,61 |
14,5 |
-7,8 |
Погрешности: Dcр = ±5,58; Dcр - от -5,87 до +4,54
Схематический график изменения физических свойств глинистых пород в зависимости от эффективного напряжения