Методика расчета интенсивности истечения газа при фонтанировании скважин
- коэффициент гидравлического сопротивления, б/р;
g - ускорение свободного падения, м/с2;
, Tн, Pн - плотность, температура, давление газа при нормальных условиях;
Tср - средняя температура в скважине, К;
Zср - среднее значение коэффициента сжимаемости;
Pз - давление в скважине напротив работающего интервала;
a, b - коэффициенты фильтрационного сопротивления;
Pа - атмосферное давление (рассматривается открытый фонтан);
a - коэффициент линейного фильтрационного сопротивления скважины, кг·см2
b - коэффициент квадратичного фильтрационного сопротивления скважины, кг·см2;
dнар i - наружный диаметр кольцевого пространства, м;
dвн i - внутренний диаметр кольцевого пространства, м;
- вязкость газа в пластовых условиях, сп.;
K - проницаемость пласта в окрестности скважины, дарси;
m - пористость пласта в окрестности скважины, б/р;
hпл - эффективная толщина пласта, м;
Zср - среднее значение коэффициента сжимаемости.
Gl(t) - массовый расход из скважины, кг/с.
Модель установившегося истечения из скважины. Предполагается, что по скважине и пласту движется установившийся поток газа. Массовый расход через любое поперечное сечение скважины одинаков:
где: Q - объемный расход через сечение;
- средняя плотность газа в сечении.
Предполагается, что в пределах скважины канал, по которому движется газ, составлен из N-прямолинейных равнопроходных секций, каждая из которых имеет поперечное сечение кольцевой формы. Таким образом, геометрия канала определяется набором следующих параметров:
li, dнар i, dвн i, , i = 1, 2, 3,..., N, (10-2)
dнар i - наружный диаметр кольцевого сечения;
- зенитный угол (угол между направлением оси секции и вертикалью).
В пределах секции справедливо уравнение количества движения:
где: di - гидравлический диаметр, вычисляемый по формуле (10-4):
l - расстояние от устья (при фонтанировании скорость в формуле (10-3) отрицательна);
Предполагается, что в местах соединения секций давление изменяется непрерывным образом (потери, вызванные изменением поперечного сечения и направления потока, не учитываются).
Уравнение состояния газа записывается в обычном виде:
где: , Tн, Pн - плотность, температура, давление газа при нормальных условиях.
Потери давления в пласте при стационарной фильтрации описываются следующим уравнением:
где: Pз - давление в скважине напротив работающего интервала;
a, b - коэффициенты фильтрационного сопротивления;
Уравнение (10-6) можно рассматривать как граничное условие для системы (формулы (10-1) - (10-3)). Условие на устье имеет вид:
где: Py, Wy, C - давление, скорость газа и скорость звука на устье;
Pа - атмосферное давление (рассматривается открытый фонтан).
Задача решается в предположении постоянства температуры и коэффициента сжимаемости:
Задача состоит в определении дебита фонтана при заданных параметрах пласта (пластовое давление, коэффициенты фильтрационного сопротивления), геометрии ствола и параметров уравнений состояния (10-1), (10-8). Рекомендуется использовать в качестве средних значений в формуле (10-8) среднеарифметические значения температуры и коэффициента сжимаемости для пластовых и устьевых условий. Решение проводится методом деления отрезка пополам. В качестве нижней границы корня принимается нулевое значение дебита. Верхняя граница определяется путем расчета забойного давления для нескольких последовательно возрастающих значений Qн.
Настоящая модель залпового выброса из скважины предназначена для расчета залпового выброса, который возникнет при мгновенной разгерметизации устья закрытой скважины. Максимальный объем поступивших в атмосферу продуктов достигается в ситуации, при которой авария происходит на скважине, заполненной неподвижным газом (например, в процессе исследований скважины).
Предполагается, что скважина вертикальна и канал, по которому происходит выброс, имеет постоянное поперечное сечение. Нестационарное течение газа описывается системой уравнений, выражающих законы сохранения массы и количества движения:
Уравнение состояния принимается в виде:
g - ускорение свободного падения;
P, W, - давление, скорость и плотность газа;
- коэффициент гидравлического сопротивления (принимается постоянным).
В случае газоконденсатных смесей при получении уравнения состояния принимается модель гомогенного потока. Предположение о равенстве скоростей фаз позволяет определить зависимость плотности от давления и температуры по данным о контактной конденсации, которое можно аппроксимировать выражением, совпадающим по форме с уравнением состояния (10-12).
В начальный момент устье скважины закрыто, и распределение давления P0 в неподвижном столбе газа описывается следующим уравнением:
На забое давление в скважине совпадает с пластовым давлением:
где: lb - координата кровли проявляющего интервала;
Pпл - пластовое давление в местах расположения скважин.
Пусть в момент t = 0 происходит мгновенная разгерметизация устья. Скорость в выходном сечении будет равна местной скорости звука, а текущий дебит вычисляется по следующей формуле:
Fh - площадь выходного сечения;
- плотность флюида в выходном сечении;
CD - коэффициент расхода, зависящий от формы выходного сечения.
Начиная с этого момента вниз по столбу газа будет передвигаться волна разрежения. Дойдя до забоя, волна, частично отразившись, перейдет в пласт, где сформируется возрастающая во времени депрессионная воронка. Для расчета выброса применяется метод смены стационарных состояний, в соответствии с которым область течения разделяется на два участка. На нижнем участке находится покоящийся столб газа, а на верхнем движется стационарный поток.
Таким образом, в приустьевой части ствола выполняются уравнения:
На подвижной границе выполняется условие:
где lf - текущее положение фронта.
Зная распределение давления по стволу, можно определить массу газа, находящегося в момент t в скважине. Из сказанного следует, что масса полностью определяется положением фронта: M(t) = M(lf). Из условия материального баланса, примененного ко всему стволу, следует уравнение перемещения фронта:
После того как волна достигла забоя, приходит в движение флюид в пласте. Предполагая течение симметричным относительно оси скважины, обозначим через Rf радиус границы (радиус депрессионной воронки), отделяющей область неподвижного флюида от прискважинной области, в которой поток стационарен и его дебит равен мгновенному дебиту фонтана. Для вычисления Rf используется уравнение, аналогичное (10-19), в котором под M понимается масса газа в стволе и круговой области пласта, радиус которой R0 выбран так, что на рассматриваемом интервале времени Rf < R0. Чтобы определить M(Rf), рассматривается задача о стационарном течении в системе "скважина-пласт", удовлетворяющем условию на устье (10-15) и условию на подвижном контуре:
Кроме того, выполняются условие сопряжения (непрерывное изменение давления) и условие массового расхода при переходе от пласта к скважине.
- Гражданский кодекс (ГК РФ)
- Жилищный кодекс (ЖК РФ)
- Налоговый кодекс (НК РФ)
- Трудовой кодекс (ТК РФ)
- Уголовный кодекс (УК РФ)
- Бюджетный кодекс (БК РФ)
- Арбитражный процессуальный кодекс
- Конституция РФ
- Земельный кодекс (ЗК РФ)
- Лесной кодекс (ЛК РФ)
- Семейный кодекс (СК РФ)
- Уголовно-исполнительный кодекс
- Уголовно-процессуальный кодекс
- Производственный календарь на 2025 год
- МРОТ 2024
- ФЗ «О банкротстве»
- О защите прав потребителей (ЗОЗПП)
- Об исполнительном производстве
- О персональных данных
- О налогах на имущество физических лиц
- О средствах массовой информации
- Производственный календарь на 2024 год
- Федеральный закон "О полиции" N 3-ФЗ
- Расходы организации ПБУ 10/99
- Минимальный размер оплаты труда (МРОТ)
- Календарь бухгалтера на 2024 год
- Частичная мобилизация: обзор новостей