СПОСОБ ОПРЕДЕЛЕНИЯ ПЛАСТОВОГО ДАВЛЕНИЯ В НЕФТЯНОЙ СКВАЖИНЕ Российский патент 2001 года по МПК E21B47/06 

Описание патента на изобретение RU2167289C2

Изобретение относится к области добычи нефти и может быть использовано для определения пластового давления в нефтяных скважинах.

Известны способы определения пластового давления и коэффициента продуктивности скважин, основанные на экспериментальных методах восстановления давления и установившихся отборов [1], [2].

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

В качестве прототипа выбран способ определения пластового давления в добывающих и нагнетательных скважинах [3], включающий остановку (закрытие) скважины, снятие с помощью глубинного манометра кривой восстановления давления, измерение забойного давления до остановки скважины, а также спустя некоторое время после ее остановки. Недостатком этого способа является то, что требуется полное снятие кривой восстановления давления и, как следствие, длительное время остановки скважины и присутствие на ней обслуживающего персонала. Кроме того, полное снятие кривых восстановления давления требует большого количества времени, ведет к значительным потерям добычи нефти и большим эксплуатационным затратам.

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

Отличие предлагаемого способа состоит в том, что измеряют давление на начальном участке кривой восстановления давления и по формулам вычисляют значение пластового давления p0, для чего измеряют глубинным манометром забойное давление pз, затем, спустя 2-3 мин, останавливают скважину и через время tк = 10-15 мин вновь измеряют забойное давление, принимая эту точку отсчета за начало координат расчетного участка кривой восстановления давления с приращением забойного давления от момента остановки скважины на величину Δ pн, и затем через каждые 5-10 мин в течение 100 мин измеряют приращение текущего значения забойного давления Δ p и по формуле

методом наименьших квадратов вычисляют текущее расчетное приращение пластового давления Δpор = Δpo-Δpн, а по нему значение конечного приращения пластового давления Δpo = Δpор+Δpн, постоянные времени T1, T2, T3 и коэффициенты A1, A2 путем сравнения в определенный момент времени давлений расчетных и действительных, снятых на начальном участке кривой восстановления давления, и далее экстраполируют кривую восстановления давления до момента, при котором разность приращений пластового давления за заданный промежуток времени будет равно нулю, и по найденным значениям вычисляют пластовое давление по формуле
p0 = Δ p0 + pз.

На фиг. 1 показаны кривые восстановления давления в глубинно-насосных скважинах.

На фиг. 2 показаны кривые восстановления давления в фонтанных скважинах.

На фиг. 3 приведен график зависимости между экспериментальными и расчетными значениями пластовых давлений по группам скважин.

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

где ω - средняя истинная скорость перемещения контура нефтеносности, м/с;
k - проницаемость пласта на данном участке, м;
m - динамическая пористость пласта, доли единиц;
p1-p2 - перепад давлений между соседними изобарами;
μ - динамическая вязкость жидкости, Па•с;
Δx - расстояние между соседними изобарами, м.

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

Рассматриваемый способ определения пластового давления основан на вышеприведенной математической модели пласта, но в ином представлении, описывающемся дифференциальным уравнением вида

где ω - постоянный коэффициент пьезопроводности пласта;
p - давление в момент t в любой точке пласта, удаленной на расстояние r от его центра.

Для плоско-радиального движения жидкости в ограниченном открытом круговом пласте известно решение этого уравнения, которое можно представить в виде более простого:

где p - пластовое давление;
p0 - текущее забойное давление;
к = 1, 2, 3,..., ∞;
Bk - коэффициент разложения решения дифференциального уравнения по экспоненциальным функциям от времени с постоянными Tк, которые однозначно находятся из известного решения уравнения (2).

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

В действительности большинство из указанных условий не соблюдается и уравнение не отражает полностью тех физических процессов, которые происходят в пласте. Поэтому известное решение уравнения (2) не отвечает действительным КВД и при соответствующем подборе коэффициентов Bк в уравнении (3) можно достичь адекватности обеих кривых. В таком случае Bк будут уже отвечать реальным граничным и начальным условиям исходного дифференциального уравнения (2). Однако - для нахождения всех коэффициентов необходимо знать всю КВД.

Для практических целей не обязательно иметь разложение решения уравнения (2) в виде бесконечной суммы. Как показывает опыт, достаточно иметь в разложении (3) суммы двух - трех слагаемых, где Bк и Tк компенсируют недостающие члены разложения. Эти коэффициенты находят из кривой КВД путем сравнения реальных давлений с давлениями, рассчитанными по аппроксимирующей кривой на начальном участке экспериментальной кривой. Характерные формы кривых для фонтанных и насосных скважин представлены на фиг. 1 и 2.

Экспериментально была опробована базовая аппроксимирующая кривая

где к = 1, 2, 3 с условием A1 + A2 + A3 = 1, или

т.е. разложения включают только три члена.

В разложении (5) неизвестно шесть параметров p0, A1, A2, T1, T2, T3.

Неизвестные параметры определяют методом наименьших квадратов путем сравнения в определенный момент времени давлений из уравнения (4) и действительных давлений.

Таким образом, находится минимум функционала:

где к = 1, 2, 3;
n - число точек, взятых на КВД в момент времени tm.

Минимум функционала (6) отыскивается методом случайного поиска, например, методом Монте-Карло для нахождения экстремума функций.

С учетом того, что перед остановкой скважины в ней имеется забойное давление, а начальный участок КВД искажен из-за притока жидкости в скважину после ее остановки, начало координат аппроксимирующей кривой сдвигается относительно основных координат на некоторое время tк (см. фиг. 1), которое принимается равным 10-15 мин после остановки насосной скважины или закрытия выкидной задвижки на фонтанной скважине.

Тогда формула (5) принимает вид:

где Δp, Δpo, Δpн - соответственно, приращения текущего забойного, пластового (депрессии) и начального (до времени tк) давлений;
Δpo-Δpн= Δpop - текущее расчетное приращение пластового давления.

Тогда конечное приращение пластового давления равно:
(8)
Пластовое же давление равно (см. фиг. 1)
p0 = Δ p0 + pз,
где pз - забойное давление на момент остановки скважины.

Рассмотрим пример расчета пластового давления по данным, приведенным в источнике - /Справочная книга по добыче нефти. Под ред. д-ра тех. наук Ш.К. Гиматудинова, М. : "Недра", 1974, 704/ - на странице 179. Эти данные - результаты исследования скважины со снятием кривой восстановления давления на забое - приведены в таблице 1.

Для расчета возьмем следующие данные:
Забойное давление pз = 120,3 кгс/см2. Время восстановления давления t = 240 мин после остановки скважины или t = 220 мин после начала расчетного участка кривой восстановления давления.

Дебит скважины - 100 т/сут.

Расчетные данные после остановки скважины приведены в таблице 2.

В результате расчета получим следующие значения (см. формулу (7)):
A1 = 0,301; A2 = 0,223; A3 = (1 - A1 - A2) = 0,476;
A1+A2+A3= 0,301+0,223+0,476 = 1.

Текущее расчетное приращение пластового давления:
Δpop= Δpo-Δpн = 3,01397 кгс/см2.

Тогда конечное приращение пластового давления равно:
Δpo= Δpop+Δpн = 3,01397 + 5,11 = 8,12397 кгс/см2.

Пластовое же давление равно:
p0 = Δ p0 + pз = 8,12397 + 120,3 = 128,42397 кгс/см2.

По экспериментальным данным, конечное приращение пластового давления равно Δ p0эксп = 8,10 кгс/см2. Следовательно, пластовое давление для экспериментальных данных равно
p0эксп = Δ p0эксп+pз = 8,10+120,3 = 128,4 кгс/см2.

Отсюда следует, что разница между расчетными и экспериментальными данными по измерению пластового давления составляет 0,01866%.

Для насосных скважин, КВД которых даны на фиг. 1, значения приращения текущего забойного давления по каждой скважине равны (см. формулу (7)):
Для первой кривой:

A1+A2+A3 = 0,847+0,175-0,028 ≈ 1.

Для второй кривой:

A1+A2+A3 = 0,893+0,218-0,11 ≈ 1.

Для фонтанной скважины (см. фиг. 2, кривая 1) зависимость приращения текущего забойного давления имеет следующий вид:

A1+A2+A3 = 0,79+0,439-0,229 ≈ 1.

В вышеприведенных формулах время измерения t берется в начале снятия расчетных точек, то есть уже учитывая время сдвига tк. Конечное приращение пластового давления и пластовое давление для этих формул рассчитываются по формулам (8) и (9).

При известном массовом дебите скважины (Q, т/сут) можно определить коэффициент продуктивности скважины как: η = Q/Δp.
Расчетный участок КВД необходимо выбирать после графического построения экспериментальной части кривой, на которой нужно взять точку на явно выраженном выпуклом участке, определяемом стягивающей прямой.

Решение базового уравнения (5) было реализовано на ЭВМ, при этом обработано более 70 кривых, снятых в основном на фонтанных скважинах ПО "Куйбышевнефть" и на глубинно-насосных скважинах ПО "Татнефть". Результаты сравнения значений приращений пластовых давлений, полученных расчетным путем, с экспериментальными значениями даны на фиг. 3.

Анализ полученных результатов показал, что расхождения расчетных данных с экспериментальными составили: для 50% скважин - до ±2,5%; для 20% скважин - до ±5%; для 15% скважин - до ±10%; для остальных скважин - более ±10%.

Дальнейшие исследования показали, что в скважинах с большим отклонением расчетных данных от экспериментальных КВД были сняты не полностью. При анализе полученных данных не учитывались и возможные различия характеристик скважин.

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

Источники информации
[1] Щелкачев В.В. Разработка нефтеводоносных пластов при упругом режиме, М.: Гостоптехиздат, 1959.

[2] Справочная книга по добыче нефти /Под ред. д.т.н. Ш.К. Гиматудинова, -М.: Недра, 1974.

[3] А.С. СССР N 1265303 A1, кл. E 21 B 47/06, опубл. 23.10.1986.

Похожие патенты RU2167289C2

название год авторы номер документа
СПОСОБ ИЗМЕРЕНИЯ КОЛИЧЕСТВА ВОДЫ, ЗАКАЧИВАЕМОЙ ЦЕНТРОБЕЖНЫМ ЭЛЕКТРОНАСОСОМ В НЕФТЯНЫЕ ПЛАСТЫ 1998
  • Кричке В.О.
  • Громан А.О.
  • Кричке В.В.
RU2176732C2
СПОСОБ НЕПРЕРЫВНОГО КОНТРОЛЯ ЗА РАБОТОЙ НАСОСНО-ТРУБОПРОВОДНЫХ СИСТЕМ ДЛЯ ПЕРЕКАЧКИ ВОДЫ И НЕФТЕПРОДУКТОВ 1997
  • Кричке В.О.
  • Громан А.О.
  • Кричке В.В.
RU2114325C1
АВТОМАТИЗИРОВАННАЯ ИНФОРМАЦИОННАЯ СИСТЕМА ДЛЯ НЕПРЕРЫВНОГО КОНТРОЛЯ ЗА РАБОТОЙ НАСОСНО-ТРУБОПРОВОДНОГО КОМПЛЕКСА ДЛЯ ПЕРЕКАЧКИ ВОДЫ И НЕФТЕПРОДУКТОВ 1997
  • Кричке В.О.
  • Громан А.О.
  • Кричке В.В.
RU2165642C2
СПОСОБ ИЗМЕРЕНИЯ ПРОИЗВОДИТЕЛЬНОСТИ НАГНЕТАТЕЛЬНОЙ УСТАНОВКИ, ОБОРУДОВАННОЙ ЦЕНТРОБЕЖНЫМ КОМПРЕССОРОМ С ПРИВОДОМ ОТ ЭЛЕКТРОДВИГАТЕЛЯ 1996
  • Кричке В.О.
RU2126526C1
СПОСОБ ИСПРАВЛЕНИЯ УГЛА I ОДНОВРЕМЕННО У ТРЕХ НИВЕЛИРОВ 1997
  • Русков А.М.
RU2138779C1
СПОСОБ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ РЕЖИМА РАБОТЫ СКВАЖИНЫ, ОБОРУДОВАННОЙ ПОГРУЖНЫМ ЦЕНТРОБЕЖНЫМ ЭЛЕКТРОНАСОСОМ 1997
  • Кричке В.О.
  • Кричке В.В.
RU2140523C1
СПОСОБ ЭКСПЛУАТАЦИИ СКВАЖИНЫ С ПОГРУЖНЫМ ЦЕНТРОБЕЖНЫМ НАСОСОМ И УСТРОЙСТВО ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ 1998
  • Галицков С.Я.
  • Люстрицкий В.М.
  • Масляницын А.П.
  • Шишков С.А.
RU2165037C2
СПОСОБ ОПРЕДЕЛЕНИЯ ТОЛЩИНЫ СЛОЯ БЕТОНА С НАРУШЕННОЙ СТРУКТУРОЙ 1993
  • Ильин Н.А.
RU2072520C1
СПОСОБ ВОССТАНОВЛЕНИЯ СЖАТЫХ ЖЕЛЕЗОБЕТОННЫХ ЭЛЕМЕНТОВ 1996
  • Ильин Н.А.
RU2116416C1
СПОСОБ УПРАВЛЕНИЯ ПИТАТЕЛЬНЫМИ НАСОСАМИ КОТЕЛЬНЫХ УСТАНОВОК 1999
  • Кричке В.О.
  • Громан А.О.
  • Кричке В.В.
RU2163702C2

Иллюстрации к изобретению RU 2 167 289 C2

Реферат патента 2001 года СПОСОБ ОПРЕДЕЛЕНИЯ ПЛАСТОВОГО ДАВЛЕНИЯ В НЕФТЯНОЙ СКВАЖИНЕ

Изобретение относится к нефтедобывающей промышленности и может быть использовано на нефтяных месторождениях для измерения пластового давления для контроля и управления процессом добычи нефти. Техническим результатом является уменьшение времени простоя скважин при исследовании, повышение точности и надежности результатов измерения пластового давления. Измеряют глубинным манометром забойное давление р3. Спустя 2-3 мин останавливают скважину и через время tk, равное 10-15 мин, вновь измеряют забойное давление, принимая эту точку отсчета за начало координат расчетного участка кривой восстановления давления с приращением забойного давления от момента остановки скважины на величину Δpн. Затем каждые 8-10 мин в течение 100 мин измеряют приращение забойного давления Δp. По приведенной формуле методом наименьших квадратов вычисляют текущее расчетное приращение пластового давления. По нему вычисляют значение конечного приращения пластового давления, постоянные времени и коэффициенты путем сравнения в определенный момент времени давлений расчетных и действительных, снятых на начальном участке кривой восстановления давления. Экстраполируют кривую восстановления давления до момента, при котором приращение пластового давления за заданный промежуток времени будет равно нулю. По найденным значениям вычисляют пластовое давление. 3 ил. , 2 табл.

Формула изобретения RU 2 167 289 C2

Способ определения пластового давления в нефтяной скважине, включающий остановку скважины, снятие с помощью глубинного манометра кривой восстановления давления, измерение забойного давления до остановки скважины, а также спустя некоторое время после ее остановки, отличающийся тем, что измеряют давление на начальном участке кривой восстановления давления и по формулам вычисляют значение пластового давления ро, для чего измеряют глубинным манометром забойное давление рз, затем спустя 2 - 3 мин останавливают скважину и через время tк = 10 - 15 мин вновь измеряют забойное давление, принимая эту точку отсчета за начало координат расчетного участка кривой восстановления давления с приращением забойного давления от момента остановки скважины на величину Δрн и затем через каждые 8-10 мин в течении 100 мин измеряют приращение текущего значения забойного давления Δр, и по формуле

методом наименьших квадратов вычисляют текущее расчетное приращение пластового давления Δp = Δp0-Δpн, а по нему значение конечного приращения пластового давления Δp0 = Δp+Δpн постоянные времени T1, Т2, Т3 и коэффициенты А1, А2 путем сравнения в определенный момент времени давлений расчетных и действительных, снятых на начальном участке кривой восстановления давления, и далее экстраполируют кривую восстановления давления до момента, при котором разность значений приращений пластового давления за заданный промежуток времени будет равна нулю, и по найденным значениям вычисляют пластовое давление по формуле р0 = Δр0 + рз.

Документы, цитированные в отчете о поиске Патент 2001 года RU2167289C2

Способ определения пластового давления в добывающих и нагнетательных скважинах 1984
  • Бучковский Станислав Степанович
SU1265303A1
Способ определения пластового давления в процессе бурения 1990
  • Бабаян Эдуард Вартанович
  • Баринов Валентин Николаевич
SU1714108A1
Способ определения пластового давления в обсаженных скважинах 1985
  • Бучковский Станислав Степанович
  • Свягла Василий Михайлович
  • Юрасик Любомир Мирославович
SU1283367A1
Способ определения пластового давления в процессе бурения 1985
  • Бабаян Эдуард Вартанович
  • Бугаенко Анатолий Александрович
  • Зуева Ирина Юрьевна
  • Суханов Владислав Борисович
SU1296717A1
СПОСОБ УСТАНОВЛЕНИЯ ПЛАСТОВОГО ДАВЛЕНИЯ НА НЕФТЯНОЙ ЗАЛЕЖИ 1997
  • Хисамов Р.С.
  • Лапицкий В.И.
  • Файзуллин И.Н.
RU2108460C1
RU 2055182 C1, 04.10.1985
Справочная книга по добыче нефти, / Под ред
д.т.н
Ш.К.ГИМАТУДИНОВА
- М.: Недра, 1974, с.169.

RU 2 167 289 C2

Авторы

Кричке В.О.

Даты

2001-05-20Публикация

1999-01-19Подача