Способ определения потерь газа при эксплуатации подземных хранилищ газа
METHOD FOR DETERMINING GAS LOSSES IN THE OPERATION OF UNDERGROUND GAS STORAGE FACILITIES

Изобретение относится к газодобывающей промышленности и может использоваться при эксплуатации подземных хранилищ газа (ПХГ).
Известен способ исследования динамических процессов газовой среды ПХГ (патент РФ №2167288, Е21В 47/00, опубл. 20.05.2001), включающий введение в пласт через разные нагнетательные скважины индикаторов в газовом носителе, отбор проб газа из добывающих скважин и определение концентраций индикаторов во времени в продукции добывающих скважин. В период максимального давления газа выбирают центральные нагнетательные скважины, расположенные в одном или нескольких эксплуатационных горизонтах, исходя из системы расположения добывающих скважин по площади, при этом используют индикаторы нескольких цветов, а закачивают индикатор одного цвета в виде газонаполненных микрогранул со степенью дисперсности 0,5-0,6 мкм, состоящих из смеси поликонденсационной смолы и органического люминесцирующего вещества в расчетном количестве. В период снижения давления до минимальной средневзвешенной по площади величины одновременно отбирают пробы газа из добывающих скважин, расположенных в одном или нескольких эксплуатационных горизонтах, и определяют изменения во времени концентрации индикаторов каждого цвета и объемной скорости газа всех добывающих скважин, находят суммарное количество индикатора каждого цвета, поступившего в каждую добывающую скважину, по заданной формуле. Строят карты и по величине долей мигрирующего газа выявляют направления внутрипластовых и межпластовых перетоков и оконтуривают газодинамически различные зоны. Недостатком известного способа является возможность только качественной оценки потерь газа и отсутствие возможности их количественной оценки.
t
t
t
t
y
y
ф
Р
P
ф
Наиболее близким к предложенному способу (прототипом) является способ определения герметичности подземных хранилищ газа (патент РФ №2526434, Е21В 47/10, опубл. 20.08.2014). В известном способе осуществляют циклическое воздействие на пласт, при котором каждый цикл включает закачку газа в пласт с последующим отбором газа. Воздействие на пласт осуществляют, по меньшей мере, в течение 10 циклов. В каждом цикле периодически одновременно измеряют текущее пластовое давление (Р) и объем отбора (или закачки) газа. С учетом измеренных параметров определяют расчетное давление в ПХГ (P) для режима эксплуатации хранилища без утечек газа и для режима эксплуатации хранилища с утечками газа. Затем определяют функцию (F) как среднеарифметическое значение отклонений (P) от (Р), полученных при каждом i-м измерении, для режима эксплуатации хранилища без утечек газа и функцию (F) для режима эксплуатации хранилища с утечками газа и при выполнении неравенства F<F делают вывод о наличии утечек газа в хранилище. Указанный способ не позволяет определить величину потерь газа при утечках из ПХГ при его эксплуатации.
Задачей, на решение которой направлено предлагаемое изобретение, является разработка способа количественной оценки потерь газа при эксплуатации ПХГ, позволяющего определять величину потерь газа из ПХГ.
Техническим результатом, на достижение которого направлено предлагаемое изобретение, является повышение точности учета газа в хранилище, надежности ПХГ и обеспечение проектных показателей при эксплуатации ПХГ.
t
t
t
t
t
t
0
п
в
ф
Р
P
ф
ф
Указанный технический результат достигается за счет того, что в предлагаемом способе определения потерь газа при эксплуатации подземных хранилищ газа (ПХГ), включающем циклическое воздействие на пласт, при котором каждый цикл включает закачку газа через эксплуатационные скважины в пласт до достижения величины пластового давления, не превышающего максимально допустимого проектного значения, с последующим отбором газа до достижения величины не ниже минимально допустимого проектного значения, в каждом цикле периодически одновременно измеряют текущее пластовое давление (Р) и объем отбора/закачки газа (q). Затем, с учетом измеренных параметров, определяют расчетное давление (P) на момент времени t, после чего определяют функцию (F), как среднеарифметическое значение отклонений (P) от (Р), полученных при каждом i-м измерении. При этом текущее пластовое давление (Р) измеряют в купольной части пласта. Расчетное давление определяют путем итерации при различных значениях газонасыщенного порового объема пласта (Ω), коэффициента потерь газа (С) и коэффициента активности пластовой воды (С) для эксплуатации ПХГ с водонапорным режимом из соотношения
о
где Р - начальное пластовое давление,
о
Z - начальный коэффициент сверхсжимаемости газа,
t
Z - коэффициент сверхсжимаемости газа на момент времени t,
t
q - объем закачки (или отбора) газа на момент времени t,
п
С - коэффициент потерь газа,
в
Q - объем подвижной пластовой воды, причем
в
где С - коэффициент активности пластовой воды,
t
P - текущее давление в залежи,
ст
Р - начальное гидростатическое давление,
0
п
или при различных значениях (Ω) и (С) для эксплуатации ПХГ с газовым режимом из соотношения
0
п
в
при этом итерацию с изменением параметров (Ω), (С), (С) выполняют до достижения минимального значения функции (F)
где n - количество замеров пластового давления,
п
i - порядковый номер замера пластового давления, после чего, с учетом полученного в результате итерации коэффициента потерь газа (С), определяют потери газа  в ПХГ на момент времени t
При эксплуатации ПХГ без количественного учета пластовых потерь газа нарушается учет объема газа в пласте, что может привести к серьезным осложнениям при создании и выходе хранилища на циклический режим эксплуатации.
Для ПХГ изменение объема газа в пласте во времени определяется уравнением
t
где V - объем газа в пласте в момент времени t,
Переходя к интегральному виду, получаем
o
где V - объем газа в начальный момент времени.
Из уравнения материального баланса (Закиров С.Н. «Проектирование и разработка газовых месторождений». М.: Недра, 1974 г., с. 28-35) известно
t
где Ω - газонасыщенный поровый объем пласта в момент времени t,
t
P - пластовое давление газа в момент времени t,
t
Z - коэффициент сверхсжимаемости газа в момент времени t.
Коэффициент сверхсжимаемости (Z) зависит от состава газа, температуры, давления и является справочным показателем (Требин Ф.А. «Добыча природного газа». М.: Недра, 1976 г., с. 78-85). Значения Z можно аппроксимировать полиномом вида
где а, b, с - коэффициенты полинома.
При водонапорном режиме эксплуатации ПХГ изменение величины Q, за время t происходит за счет движения пластовой воды при закачке/отборе газа в ПХГ, тогда уравнение (8) для ПХГ с водонапорным режимом примет вид
Эксплуатацию ПХГ с водонапорным режимом описывают через измеряемые параметры отбора (закачки) газа и пластового давления системой уравнений
t
t
При газовом режиме эксплуатации ПХГ изменения порового объема Ω за время t не происходит, т.к. Ω - постоянная величина и уравнение (8) для ПХГ с газовым режимом примет вид
Эксплуатацию ПХГ с газовым режимом описывают через измеряемые параметры отбора (закачки) газа и пластового давления системой уравнений
При наличии потерь газа при эксплуатации ПХГ уравнение (6) примет вид
- объем потерь газа при эксплуатации ПХГ в момент времени t.
Объем потерь газа при эксплуатации ПХГ можно описать уравнением вида (Закиров С.Н. «Проектирование и разработка газовых месторождений». М.: Недра, 1974 г., с. 220-226)
п
где С - коэффициент потерь газа при эксплуатации ПХГ.
При эксплуатации ПХГ с водонапорным режимом с потерями газа уравнение (11) примет вид
Уравнение (17) с учетом движения пластовой воды при закачке (отборе) газа в ПХГ примет вид
При эксплуатации ПХГ с газовым режимом с потерями газа уравнение (13) примет вид
t
P
Значение расчетного давления в ПХГ определяют путем решения формул (18) и (19) относительно (P).
t
t
t
t
P
Ф
P
Ф
Для оценки отклонения расчетного пластового давления (P) от фактического (Р) используют функцию (F), определяемую как среднеарифметическое значение абсолютных отклонений (Р) от фактического (Р), полученных при каждом i-м измерении
t
0
п
в
P
При эксплуатации ПХГ с водонапорным режимом функцию (F) определяют путем решения уравнения (18) относительно пластового давления (Р) при различных значениях параметров (Ω, C, С).
t
0
п
P
При эксплуатации ПХГ с газовым режимом функцию (F) определяют путем решения уравнения (19) относительно пластового давления (Р) при различных значениях параметров (Ω, C).
Способ осуществляют следующим образом.
t
t
t
t
0
п
в
0
п
min
п
Ф
Ф
P
В процессе эксплуатации ПХГ осуществляют циклическое воздействие на продуктивный пласт. В каждом цикле через эксплуатационные скважины проводят закачку газа в продуктивный пласт с последующим отбором газа. Закачку газа проводят до достижения пластового давления в ПХГ, не превышающего максимально допустимого проектного значения. Отбор газа проводят до достижения пластового давления не ниже минимально допустимого проектного значения. В течение каждого цикла замеряют текущее пластовое давление (Р) в хранилище и объем отбора (закачки) газа (q). Текущее пластовое давление (Р) измеряют в купольной части пласта, что позволяет получить наиболее достоверные значения измеряемых параметров и обеспечивает повышение точности учета газа. Затем определяют расчетное давление (Р) по формуле (18) для ПХГ с водонапорным режимом и по формуле (19) для ПХГ с газовым режимом. После чего определяют значения функции (F) по формуле (20) путем итерации при различных значениях параметров (Ω, C, С) для ПХГ с водонапорным режимом и при различных значениях (Ω, C) для ПХГ с газовым режимом. Итерацию выполняют до момента, когда значение функции (F) перестает изменяться и определяют минимальное значение функции (F). По найденному в результате итерации значению коэффициента потерь газа (С) определяют по формуле (16) объем потерь газа, количественное значение которого позволяет оценить риски при эксплуатации ПХГ и обеспечить проектные показатели и надежность ПХГ. Пример 1.
t
t
t
Ф
Р
С помощью предлагаемого способа исследовали Калужское ПХГ с водонапорным режимом эксплуатации. Полученные в процессе исследования фактические значения пластового давления (P) и объема закачки (отбора) газа (q), а также расчетные значения пластовых давлений (P) и потерь газа  в различные моменты времени приведены в таблице 1.
п
min
п
По результатам сравнения фактических и расчетных параметров и найденному в результате итерации значению коэффициента (С) были определены потери газа в различные моменты времени в указанном ПХГ при F=2,823 и С=0,00001.
Пример 2.
t
t
t
ф
P
С помощью предлагаемого способа исследовали Осиповичское ПХГ с газовым режимом эксплуатации. Полученные в процессе исследования фактические значения пластового давления (Р) и объема закачки (отбора) газа (q), а также расчетные значения пластовых давлений (P) и потерь газа  в различные моменты времени приведены в таблице 2.
П
min
П
По результатам сравнения фактических и расчетных параметров и найденному в результате итерации значению коэффициента (С) были определены потери газа в различные моменты времени в указанном ПХГ при F=1,92 и С=0,037.
Таким образом, определение реальных количественных потерь газа с использованием наиболее достоверных, измеренных в купольной части пласта, значений текущего пластового давления в различные моменты времени позволяет повысить точность учета газа в хранилище, надежность ПХГ и обеспечить проектные показатели при эксплуатации ПХГ.