CC BY
36
ГАЗОВАЯ ПРОМЫШЛЕННОСТЬ
УДК 622.691 105
Об уравнениях движения газоконденсатных систем в пористых средах и трубах вблизи давления конденсатообразования
Р.М. Саттаров
доктор технических наук, профессор, научный консультант1 r.sattarov@yahoo.com
Б.З. Тухфатов
кандидат технических наук, генеральный директор2 tuxfatov@mail.ru
1НИИ «Геотехнологических Проблем Нефти, Газа и Химии», Баку, Азербайджан 2КГП «ОблТрансГаз», Атырау, Казахстан
Эффективность технологических процессов разработки и эксплуатации газовых и газоконденсатных месторождений неразрывно связана с выбором и установлением в системе пласт-скважина-трубопроводные коммуникации термогидрогазодинамически оптимальных режимов течения газожидкостной смеси, сопровождаемые, как правило, неравновесными явлениями и фазовыми переходами. Неравновесные процессы и фазовые переходы в газожидкостных системах проявляются в достаточно широком диапазоне давлений, включая и области выше и ниже критических давлений, при этом, как правило, релаксационные свойства отмеченных смесей зависят как от уровня давления, так и от темпа изменения давления в системе.
Поскольку при движении газожидкостных систем происходит снижение давления, то естественно предположить, что по мере движения может происходить образование и рост новой фазы (на нано и микро уровнях)конденсации, часть из которых, как правило, выносится газовым потоком, а часть скапливается на стенках пор или труб, изменяя гидравлические характеристики среды.
Известно, при разработке и эксплуатации газовых и газоконденсатных месторождений возникает необходимость исследования те рмогидро газодинамического состояния флюидов, как в свободном объеме, так и в пористых средах, поскольку контроль и управление технико-технологическим процессом газо-конденсатодобычи требует детального теоретического анализа процессов переноса углеводородного сырья на всех этапах добычи.
При этом моделирование технико-технологических процессов и выбор оптимальных режимов добычи газа и конденсата представляет важный практический интерес не только на этапе проектирования технологических схем, но и на любой стадии разработки месторождения.
Причем особый интерес представляют исследования фазового состояния и кинетики фазовых превращений флюидов в различных скважинно-трубопроводных и пластовых условиях, поскольку кинетика восстановления термодинамического равновесия, как правило, может нести уникальную информацию, как о самих флюидах, так и о свойствах свободного объема и порового пространства.
Построение рациональной системы эксплуатации газовых и газоконден-сатных месторождений на основе установления термогидрогазодинамически оптимальных режимов течения газожидкостной смеси в системе разработки месторождения, системах транспортировки, отличается тем, что добываемые природные углеводороды представляет собой состав компонентов, при движении которых, как правило, происходят неравновесные процессы, сопровождаемые фазовыми переходами.
В принципе, можно считать установленным, что неравновесные процессы в газоконденсатных системах могут проявляться как при давлениях несколько выше давления начала конденсации, когда происходят процессы образования зародышей (на нано и микро уровнях) новой фазы, так и в дальнейшем, по мере снижения давления, когда происходит кластеризация новых фаз на более высоком масштабном уровне.
Кроме того, как показывают экспериментальные исследования, в областях критических термодинамических параметров, при фазовых переходах, наряду с процессами конденсации, происходят процессы испарения (обратный фазовый переход из газоконденсатного в однофазное газовое состояние) конденсата, что также предопределяет особый характер неравновесного поведения системы.
В такой ситуации, принятия адекватных инженерных решений для повышения
эффективности разработки и эксплуатации систем добычи природного газа и конденсата требует, прежде всего, надежного моделирования процессов течения газожидкостных углеводородных смесей в пластовых и трубопроводных системах, разработки достаточно упрощенных методов расчета, с учетом термогидрогазоди-намически неравновесных особенностей многофазных систем.
Особенностям неравновесных процессов газожидкостных систем в достаточно широком диапазоне давлений, включая и области выше и ниже критических давлений, посвящены ряд исследований, в которых показано, что релаксационные свойства отмеченных смесей зависят как от уровня давления, так и от темпа изменения давления в системе [1, 2, 3, 4, 5].
Следует отметить, что результаты отмеченных исследований находятся в соответствии с кинетической теории фазовых превращений по Я.И. Френкелю [6], согласно которой возникновение новой фазы (на нано и микро уровнях) начинается задолго до критической точки под действием локальных флуктуаций.
В то же время, для оценки релаксационных свойств неравновесных газожидкостных систем, как правило, используются методы идентификационного моделирования [1, 3, 4].
Такой подход, на основе полученных идентификационных моделей не всегда позволяет выявлять физическую сущность исследуемых процессов.
Поэтому, в этом случае, сочетание идентификационных и термогидрогазоди-намических методов моделирования для получения некоторых уравнений состояния и течений в пористых средах и трубах газожидкостных систем, в предпереход-ных областях с учетом неравновесных фазовых переходов, является наиболее целесообразным.
Поскольку при движении газожидкостных систем происходит снижение давления, то естественно предположить, что по мере движения может происходить образование и рост новой фазы (на нано и микро уровнях) конденсации, часть из которых, как правило, выносится газовым потоком, а часть скапливается на стенках пор или труб, изменяя гидравлические характеристики среды.
В соответствии с изложенным, можно записать зависимость плотности газожидкостной системы от концентрации вновь образованной новой фазы (на нано и микро уровнях) конденсации
Р^ =Рё (1 - С)+Р/С (1)
рд1 — плотность газожидкостной системы; р — плотность газа;
^ д
рг — плотность жидкости;
С — концентрация вновь образованной новой фазы (на нано и микро уровнях) конденсации.
Предполагается, что при образовании новой фазы конденсации на наноуровне, возникающие центры жидкостных нанов-ключений, с уменьшением давления могут кооперироваться первоначально в нано-кластеры, а затем на микроуровне в жидкостные микровключения, равномерно перераспределенные в газожидкостном потоке.
Предполагается также, что термодинамическая зависимость плотности газа от давления Р и температуры Т может быть записана как рд = р(Р,Т).
Однако если допустить, что в среде отсутствуют источники выделения или поглощения энергии, а изменение температуры в процессе движения газожидкостной системы крайне малы и остается приблизительно постоянной и равной Т0, то уравнение термодинамического состояния для газа примет вид [7]
(2)
а плотность жидкости, при этом в силу также своей слабой сжимаемости по сравнению с плотностью газа, практически принимается постоянной pf = pfo = const.
Поскольку рассматриваемые процессы, из-за времени образования, кооперирования, роста и перераспределения на нано и микро уровне, носят неравновесный характер, то зависимость концентрации жидкостных включений C от уровня давлений P может быть записана как
(3)
или после некоторых преобразований как
О
dt
е,=е
Ф, dp
dt
(5)
Полученная зависимость (5), описывающая процессы конденсации в газожидкостной системе в интервале от давления вновь образованной новой фазы (на нано и микро уровнях) конденсации, является нелинейным неравновесным термодинамическим уравнением состояния газожидкостной смеси.
Нелинейное неравновесное термодинамическое уравнение состояние (5), при определенных принятых предположениях, может быть существенно упрощено.
Если предположить, что газ подчиняется закону термодинамически идеальному газу
(6)
где:
в — характерное время образования жидкостных включений в газожидкостной системе; вс — коэффициент вновь образованной новой фазы
(на нано и микро уровнях) конденсации; Рс — давление вновь образованной новой фазы (на нано и микро уровнях) конденсации; t— время.
Предложенное уравнение (3) может описывать процессы конденсации в газожидкостной системе в интервале от давления вновь образованной новой фазы (на нано и микро уровнях) конденсации Рс (выше отмеченного давления в газе полностью отсутствуют жидкие включения), до уровня давления, когда система не может рассматриваться как гомогенная газожидкостная среда.
Уравнение (1) относительно С можно записать в следующем виде
(4)
С учетом соотношений (2) и (4), а также предполагая, что плотность жидкой фазы постоянна и равна р0, дифференциальная зависимость (3) запишется
(р и Р — постоянные), а также полагая
^ до о "
для параметра вс правомерность замены, в первом приближении, плотности газожидкостной смеси плотностью газа рд! =рд(Р), а для параметра всс правомерность замены величины давления Р ее осред-ненным значением Рт, то уравнение (5), в упрощенном линеаризованном виде, запишется
ndpgf +р -й dP + в P
o.„=o
к'
Материалы и методы
1. Зависимость плотности газожидкостной системы от концентрации вновь образованной новой фазы (на нано и микроуровнях) конденсации принята на основе термодинамических методов. Использовано также, уравнение термодинамического состояния для газа, а в силу своей слабой сжимаемости по сравнению с плотностью газа, принято допущение о постоянстве плотности жидкости.
2. Предложена неравновесная зависимость концентрации жидкостных включений от уровня давлений, описывающая процессы конденсации в газожидкостной системе в интервале от давления вновь образованной новой фазы (на нано и микроуровнях) конденсации Рс (выше отмеченного давления в газе полностью отсутствуют жидкие включения), до уровня давления, когда система не может рассматриваться как гомогенная газожидкостная среда.
3.Построено нелинейное неравновесное термодинамическое уравнение состояния газожидкостной смеси, описывающее процессы конденсации в газожидкостной системе в интервале от давления вновь образованной новой фазы (на нано и микроуровнях) конденсации, до уровня давления, когда система не может рассматриваться как гомогенная газожидкостная среда.
4. Получено дифференциальное уравнение изотермической фильтрации газоконден-сатных систем, с учетом эффектов неравновесности в области образований новых фаз (на нано и микроуровнях) конденсации.
5. Получена система нелинейных уравнений, описывающая изотермические процессы течения газоконденсатных систем в скважинах и трубах, с учетом образований новых фаз (на нано и микроуровнях) конденсации.
Ключевые слова
термогидрогазодинамика, флюид, газодобыча, газоконденсатная система, нано, микро, дифференциальное уравнение
(7)
Для описания фильтрационных процессов газоконденсатных систем в пластах, с учетом эффектов неравновесности в области образований новых фаз (на нано и микроуровнях) конденсации, можно воспользоваться условием, что
сжимаемость газожидкостнои смеси на несколько порядков превышает сжимаемость пористоИ среды. Тогда уравнение неразрывности запишется в виде
:■■■. - - ■... :■ (8)
где:
m — пористость пласта;
U, — скорость фильтрации газожидкостной смеси.
Предполагается, что уравнение фильтрации газожидкостноИ смеси описывается законом Дарси
где:
k — коэффициент проницаемости; ц — вязкость газожидкостной смеси.
Поскольку, как отмечалось выше, часть образовавшегося конденсата (на нано и микроуровнях), как правило, выносится газовым потоком, а часть скапливается на стенках пор, то коэффициент проницаемости и вязкость газожидкостноИ смеси в общем случае, с учетом изотермической фильтрации, могут быть записаны в следующем виде:
(10)
Теперь вводятся функции
ц
(и)
Уравнение (8), с учетом (7), (9) и (11), после некоторых преобразованиИ (в уравнении (7) субстанциональные производные, в силу пренебрежения конвективными составляющими, заменяются локальными производными), может быть записано как
= О
at
ij^i
(12)
— оператор Лапласа, х, у, z — пространственные координаты.
В простейшем случае, можно газо-конденсатную смесь (с конденсатными включениями на нано и микроуровнях) считать термодинамически идеальной, с постоянной, не зависящей от давления, вязкостью цд)(Р) = ц°д1 const, а коэффициент проницаемости пласта, также считаются постоянным k(P) = k° = const.
При этом плотность рассматриваемой смеси в первом приближении принимается в соответствии зависимость- И-также является постоянной.
Тогда, с учетом зависимостей функций
Р(Р)= }
(13)
уравнение (12) преобразуется к виду
Mv а )
м.
(14)
(9)
<
(15)
О—+1 Iii11 dl
Поскольку, уравнения (14) и (15) получены в предположении постоянства температуры газожидкостных смесей Т, то они являются дифференциальными уравнениями изотермической фильтрации газоконденсатных систем, с учетом эффектов неравновесности в области образований новых фаз (на нано и микроуровнях) конденсации.
Предложенные дифференциальные уравнения (14) и (15), в частном случае, когда не имеют место фазовые образования конденсаций (6 = о, вс =0, вс0 = р0^ / Ро и = ^08), полностью совпадают с дифференциальными уравнениями изотермической фильтрации газа [7].
Для описания процессов течения га-зоконденсатных систем в скважинах и трубах, с учетом эффектов неравновесности в области образований новых фаз (на нано и микроуровнях) конденсации, можно воспользоваться одномерным уравнением неразрывности
(16)
— осредненная скорость газожидкостной смеси в трубе.
Уравнение движения газожидкостной смеси в трубах, при этом, может быть записан в следующем виде
(17)
г у£ра
(11)
di dx
a§L-v
dx P*gdx 2D
где:
§ — ускорение свободного падения; й — диаметр трубы;
Л — коэффициент гидравлического сопротивления; х — координата, совпадающая с осью трубы и направленная по течению газожидкостной смеси; г — превышение центра тяжести сечения трубы над горизонтальной плоскостью.
Система уравнений (7), (16) и (17), описывающая изотермические процессы течения газоконденсатных систем в скважинах и трубах, с учетом образований новых фаз (на нано и микроуровнях) конденсации, является системой нелинейных уравнение, решение которой относительно P, рд/ и w сопряжено с определенными трудностями. Поэтому для их решения, можно прибегнуть к некоторым упрощениям и различного типа линеаризациям, в зависимости от конкретных технологических задач.
Итоги
Предложенная система дифференциальных уравнений (7), (16) и (17), в частном случае, когда не имеют место фазовые образования конденсаций (в = о, ßc =0, ßс = Р0д/ /Po, ц0дf = цод, pgo = рд и wgo = wg), полностью совпадают с системой дифференциальных уравнений течения газа в скважинах и трубах.
Выводы
Полученные уравнения движения газоконденсатных систем в пористых средах и трубах вблизи давления конденсато-образования, могут быть использованы при решении конкретных инженерных задач контроля и регулирования образований новых фаз (на нано и микроуровнях) конденсации для повышения эффективности различных технологических процессов газодобычи и газоконденсатоотдачи пласта.
Список используемой литературы:
1. Саттаров Р.М., Фарзане П.Я. Исследование движения газожидкостных систем с учетом образования микрозародышей //
ИФЖ. 1987. т. 52. № 5. С. 765-771.
2. Болотов А.А., Мирзаджанзаде А.Х., Нестеров И.И. Реологические свойств растворов газов
в жидкости в области давления насыщения // МЖГ. 1988. № 1. С. 172-175.
3. Саттаров Р.М. Неустановившееся движение реологически сложных жидкостей в трубах. Баку: Элм, 1999. 412 с.
4. Мирзаджанзаде А.Х., Хасанов М.М., Бахтизин Р.Н. Этюды о моделировании сложных систем нефтдобычи. Уфа: Гилем, 1999. 464 с.
5. Мирзаджанзаде А.Х., Кузнецов О.Л., Басниев К.С.,
Алиев З.С.
Основы технологии добычи газа. М.: Недра, 2003. 88о с.
6. Френкель Я.И. Кинетическая теория жидкостей. М.: Наука, 1975. 592 с.
7. Баренблатт Г.И., Ентов В.М., Рыжик В.М.
Теория нестационарной фильтрации
жидкости и газа.
М.: Недра, 1972. 288 с.
About the movement equations of gas-condensate system in porous udc 622.691
mediums and pipes nearby a pressure of condensate formation
Authors:
Rakiz M. Sattarov — doctor of science, scientific adviser1; r.sattarov@yahoo.com Bauyrzhan Z. Tukhfatov — candidate of science, general director2; tuxfatov@mail.ru
1Research Institute "Geotechnological Problems of Oil, Gas and Chemistry", Baku, Azerbaijan 2MSE "OblTransGas", Atyrau, Kazakhstan
ENGLISH GAS INDUSTRY
Abstract
Efficiency of technological processes
of development and exploitation of gas
and gas-condensate fields inseparably
linked with a choice and an establishment
thermohydrogasdynamic optimal a
flow regime of liquid-gas mixtures
accompanied, as a rule,
by the nonequilibrium phenomena and
phase transitions in
system of stratum-well-pipeline
communications.
Nonequilibrium processes and phase transitions in liquid-gas mixtures are revealed in sufficient wide range of pressure, including and areas above and below critical pressure, thus, as a rule, relaxational properties of noted mixes depend both as on pressure level so and from change rate of pressure in system. Since under a movement of gas-condensate systems is a pressure decreasing, then in process of movement there can be a formation and growth of a new phase (on nano and micro levels) condensation, the part from which, as a rule, is taken out by a gas stream, and the part accumulates on pore walls or pipes, changing hydraulic characteristics of medium.
Materials and methods
1. Dependence of density of gas-liquid system on concentration of repeatedly formed new phase (on nano and micro levels) condensation is accepted on the basis of thermodynamic methods. The equation of a thermodynamic condition for gas,and owing to the weak compressibility in comparison with gas density is also used,
the assumption about constant density of liquid is accepted.
2. No equilibrium is dependent of concentration of liquid inclusions on level of pressure, describing condensation processes in gas-liquid
system in the range from pressure of again formed new phase (on nano and micro levels) is offered Pc (at constant pressure) condensation (above noted pressure in gas completely there are no liquid inclusions), to pressure level when the system can't be considered as the homogeneous gas-liquid environment.
3. The nonlinearity of no equilibrium thermodynamic equation of a condition of the gas-liquid mix,
describing condensation processes
in gas-liquid system in the range from
pressure of again formed
new phase (on nano and micro levels)
condensation,
to a pressure
level when the system can't be considered as the homogeneous gasliquid environment is formed.
4. The differential equation of an isothermal filtration of gas-condensate systems, taking into account effects in the field of formations of new phases (on nano and micro levels) condensation is received.
5. The system of the nonlinear equations describing isothermal processes of a current of gas-condensate systems in wells and pipes, taking into account formations of new phases (on nano and micro levels) condensation is received.
Oonclusions
The equations of movement of gas-condensate systems in porous environments and pipes near the condensate formation pressure which can be used at the solution of specific engineering objectives of control and regulation of formations of new phases (on nano and micro levels) condensation; for increase of efficiency of various technological processes of gas production and gas condensate output layer are achieved.
Keywords
thermohydrogasdynamics, a liquid, gas-extraction, gas-condensate system, nano, micro, the differential equation
References a liquid in the saturation pressure Ufa: Gilem, 1999, 464 p.
1. Sattarov R.M., Farzane P.Ya. region]. 5. Mirzadjanzadeh A.Kh., Kuznetsov
Issledovanie dvijeniya gazojidkostnikh MJG, 1988, issue 1, pp. 172-175. O.L, Basniev K.S., Aliev Z.S. Osnovi
system s uchetom obrazovaniya 3. Sattarov R.M. Neustanovivsheesya tekhnologii dobichi gaza [Basic
mikrozarodishey [Motion of gas- dvijenie reologicheski slojnikh foundation of technology of a gas
liquid systems, taking account of jidkostey v trubakh production]. Moscow: Nedra, 2003,
micronucleus formation]. Journal [Nonsteady movement of rheological 880 p.
of Engineering Physics and complex fluids in pipes]. 6. Frenkel Ya.I. Kineticheskaya teoriya
Thermophysics, 1987, Vol. 52, issue 5, Baku: Elm, 1999, 412 p. jidkostey [The kinetic theory of fluids].
pp. 555-561. 4. Mirzadjanzadeh A.Kh., Khasanov M.M., Moscow: Nauka, 1975, 592 p.
2. Bolotov A.A., Mirzadjanzadeh A.Kh., Bakhtizin R.N. 7. Barenblatt G.I., Entov V.M., Rijik V.M.
Nesterov I.I. Reologicheskie svoystva Etudi o modelirovanii slojnikh system Teoriya nestatsionarnoy filtratsii
rastvorov gazov v jidkosti v oblasti neftedobichi jidkosti i gaza [The theory of a
davleniya nasisheniya [Rheological [Etudes about modelling of complex nonstationary filtration of a fluid and
properties of solutions of gases in systems of oil production]. gas]. Moscow: Nedra, 1972, 288 p.
