Численное моделирование температуры и давления в газовой скважине Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

УДК 550.3

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРЫ И ДАВЛЕНИЯ В ГАЗОВОЙ СКВАЖИНЕ

© Н. Ф. Файзуллин*, А. Ш. Рамазанов

Башкирский государственный университет Россия, Республика Башкортостан, 450074 г. Уфа, ул. Заки Валиди 32.

Тел./факс: +7 (347) 272 60 56. *Email: dr.faizmeister@gmail. com

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

Ключевые слова: моделирование, газ, конвекция.

Введение

Технологический режим работы газовой скважины определяется термодинамическими условиями фильтрации газа в пласте и условиями движения газа в скважине при теплообмене с окружающими горными породами [8, 9]. При этом условия тепло-переноса при фильтрации газа в пласте и движении потока газа по скважине различны, однако их параметры определяются аналогичными термодинамическими процессами: дроссельным и адиабатическими эффектами, конвективным и кондуктивным переносом тепла.

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

Постановка задачи

Для моделирования термогидродинамических процессов, происходящих в скважине с произвольной траекторией и произвольным числом интервалов перфорации решается следующая система уравнений неизотермического установившегося течения газа со следующими допущениями: течение однофазное (фазовые переходы отсутствуют), изменением значения величин по поперечному сечению скважины пренебрегаем, в пределах интервалов перфорации задаются объемный дебит и температура поступающего из пласта газа. Задается состав газа (здесь для расчетов взят следующий состав газа для Уренгойского месторождения: 97.8% метан, 0.1% этан, 0.03% пропан, 0.01% группа пентана и выше, 0.18% углекислый газ, 1.7% азот). Течение газа по колонне, свободной от НКТ.

скважина, давление, теплообмен, теплоперенос,

d(pv) 2pjvj

д I

dp л

дТ dp ,

— = — £— +

dl dl Rpv

R

d(pv2) dl

2pivi

2

pv2 r 4R '

+

(T, — T) + ( Tg — T) — ?cose

(1)

ЯрУСр - - - ср

Здесь р - плотность газа в скважине, кг/м3; и -средняя скорость течения газа по трубе, м/с; I - расстояние от забоя по стволу скважины, м; р1 - плотность газа, втекающего в скважину, кг/м3, и1 - скорость притекающего газа, м/с; р - давление, Па; в -угол наклона участка скважины от вертикали, град.; / - коэффициент трения, безр.; ц - вязкость газа, Па с; X - теплопроводность газа, Вт/(мК); Т - температура, К; а - коэффициент теплопередачи, Вт/(м2 К); Я -внутренний радиус обсадной колонны, м; - невозмущенная температура горных пород, К.

Первое уравнение описывает баланс масс, второе уравнение является законом распределения давления [2], третье - уравнением баланса энергии. В правой части уравнения энергии первое слагаемое совместно с четвертым отвечает за изменение температуры в результате охлаждения за счет расширения и диссипативный нагрев, второе - за эффект калориметрического смешивания в продуктивных интервалах, третье - учитывает конвективный теплообмен с горными породами. К ним добавляется уравнение состояния газа 2(р, Т) для нахождения р(р, Т), е(р,Т), ср(р,Т) и корреляционные соотношения для ¡(р,Т), Х(р,Т).

Коэффициент теплопередачи, учитывающий тепловое сопротивление материалов скважины и горных пород равен [3]:

1

& = ~ 1 гс 1 (2) --+---+—с---

кс Ьсш к^

где Ис - коэффициент теплоотдачи от газа к стенке скважины, Вт/(м2 К), Нс^ - коэффициент теплопередачи от внутренней стенки обсадной колонны к

2а

стенке скважины, Вт/(м2К), к/ - коэффициент теплопередачи между стенкой скважины и горными породами, с невозмущенной температурой, Вт/(м2К).

Уравнение состояния реального газа

В качестве уравнения состояния реального газа

возьмем уравнение Пенга-Робинсона: _ пт

Р Т7_Ь 1

v-ь v(v+b)+b(b-v)

(3)

где

а - 0.45724'

Рс

Ъ - 0.077Q —

(4)

1+m(l

m - 0.37464 + 1.54226м - 0.26992м2 (5)

(6)

Рс

приводится к уравнению третьего порядка относительно коэффициента сверхсжимаемости [4]:

г3 +г2(в -1) + г (а -2В - зв2) +

+ (АВ - В2 - В3) = 0 (7)

где коэффициенты А и В определяются следующим образом:

А=^,В = Р± (8)

п2т2 пт у '

где Тс - критическая температура, К, рс - критическое давление, Па, ю - ацентрический фактор.

Расчет параметров на основе уравнения состояния

рут = гят (9)

Р = РТ (1°)

где Я - универсальная газовая постоянная, Дж/(моль К), М - молярная масса газа, кг/моль, Ъ -коэффициент сверхсжимаемости, безразм., Ут=М/р - молярный объем, м3/моль.

Коэффициент Джоуля-Томсона и адиабатический коэффициент

=

1-ßT

СрР ßT

v - —

СрР

(11) (12)

Здесь в - коэффициент температурного расширения, К-1,

ет = 1+Шр (13)

Теплопроводность

Теплопроводность X может быть рассчитана по формулам при известных ?т1х (теплопроводность смеси при р = 0.1 МПа), критической величине сверсжимаемости 2с, приведенной плотности газа ра и параметра зависящего от молекулярной массы смеси и критических давлении и температуры [5].

С точностью до 5% значение Щт^ определяем по формуле

_ Т,?=оУ1ЛаЦМ1/з

Ant --„ 1 /

(14)

£?=оУ1м1/з

где Хац - теплопроводность /-го компонента при р=0.1 МПа.

Для метана, нафтеновых и ароматических углеводородов при р=0.1 МПа

- 2.047 • 10-3 • ер»/!-

и 1/ М /2Т1Ь

5 = ^-

(15)

(16)

A -

Ant +

Рс

14^10-8(e-0535Pd-l)

Jz5 '

Pd < 0.5

An+ +

An+ +

13.b10-s(e-0(i7Pd-1,069) JZ5 '

0.5 <Pd<2

2.976^10-s(e-1155pd + 2.016)

(17)

где Ра

-P/Pd =

ZcTcP

Pd> 2

/ р =--приведенная плотность, без-

га гтрс

разм., 1С = г(рс,Тс) - коэффициент сверхсжимаемости при критических параметрах, безразм.

Вязкость

Вязкость газа при данном давлении и температуре, согласно компонентному составу, определяется по формуле

1/

_ ZtoWiM/2 V — 1/

(18)

где - вязкость /-го компонента при заданной температуре и давлении.

Согласно методу Голубева, основанному на преобразовании Трауца [6]

Vc • Td0965,Td < 1

Vi =

Wc • Td

9/т

d,Td > 1

1/ 2/3 1 /2 •pj

Vc

(19)

(20)

где Та = Т/т - приведенная температура, безразм.,

' 'с

цс - вязкость при критических параметрах, мкПа с.

Теплоемкость

Теплоемкость природных газов зависит от давления, температуры и состава газа, а также от процесса передачи тепла. Теплоемкость можно вычислить через идеально-газовую теплоемкость и фуги-тивность /, характеризующую отклонение газа от идеального:

ОТ21

_ ,±_

Ср = Ср0 + дт

дт ).

RT2

(21)

Фугитивность определяется из соотношения фугитивности и энтальпии:

fainA _ н0(т)-н(р,т)

д T Р

(22)

ыг-ыР + !Р(2-1)Р

В общем виде фугитивность определяется из уравнения состояния:

(23)

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

Ln^ = z-1-ln(z-В)-

V.

/б

T

с

Рис. 1. Зависимость коэффициента сверхсжимаемости (слева) и вязкости (справа) от давления при заданной температуре. 1 - рассчитанные в симуляторе, 2 - данные МБТ.

1

2^2 b

z+(^2 + I)B dT ") z-(fi-\)B

Н-Н0 = RT(z -1) +

1 ( da z+{^2+i)b +--F" ( T--a ) In —77=—■—

2\2b \ dT ) z—yy2 — 1)B

(24)

(25)

Граничные условия

1. На забое скважины задается давление:

Р1=0 = Ро (26)

2. Распределение температуры горных пород подчиняется геотермическому распределению.

То(г) = Тд(г) (27)

Решение и тестирование

Задача (1-27) решалась численно методом Эйлера, с использованием разностной схемы по направлению потока на равномерной сетке.

1. Проверка расчета параметров газа

Как было сказано выше, параметры газа рассчитываются в зависимости от давления и температуры. Тестовые расчеты проводились для одного компонента - метана. Приведены примеры сравнения значений коэффициента сверхсжимаемости и вязкости с их значениями, взятыми из базы NIST [7] при постоянной температуре 100 °С.

Для сверхсжимаемости СКО составило 0.019, а для вязкости СКО составило 0.0029 сПз.

2. Сравнение расчетов с постоянными параметрами и параметрами, зависящими от давления и температуры

Рассмотрим вертикальную скважину глубиной 1200 м, радиусом скважины по долоту 12 см, внутренним радиусом обсадной колонны 10 см. Забойное давление составляет 40 атм, геотермический градиент равен 0.02 К/м, температура горных пород на забое 1=0 равна 45 °С, зумпф скважины заполнен водой, скважина имеет единственный продуктивный интервал [1000-1030] м, дебит с этого интер-

вала 50000 ст. м3/сут, температурный эффект Джо-уля-Томсона в пласте, составляет -20 К. В модели с постоянными коэффициентами примем значения параметров согласно табл. 1.

Таблица 1

Свойства газа для модели с постоянными параметрами | Плотность р | 15 кг/м3 |

Удельная теплоемкость ср Коэффициент Джоуля-Томсона е

Вязкость ^ Теплопроводность Л

3 кДж/(кгК)

-0.3 К/атм

0.015 сПз 0.1 Вт/м/К

На рис. 2 показаны смоделированные профили температуры и давления для случая с постоянными параметрами газа и с параметрами, зависящими от температуры и давления. Для температуры СКО составляет 0.55 К, для давления - 0.13 атм, максимальное отклонение, соответственно, 0.76 К и 0.22 атм.

Выводы

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

ЛИТЕРАТУРА

1. Чекалюк Э. Б. Термодинамика нефтяного пласта. М.: Недра, 1965. 238 с.

2. Sui W. Determining multilayer formation properties from transient temperature and pressure measurements: PhD dissertation, Texas A&M University, 2009.

3. Hasan, A. R. and C. S. Kabir, 2002. Fluid Flow and Heat Transfer in Wellbores. Society of Petroleum Engineers, Richardson, Tex, pp: 181.

4. Садретдинов А. А. Неизотермическая фильтрация сжимаемого флюида в системе скважина-пласт: дис. канд. физ.-мат. наук. БашГУ, Уфа, 2012. 125 с.

5. Зотов Г. А. Руководство по исследованию скважин. М.: Наука, 1995.

6. Глумов Д. Н., Стрекалов А. В. Способ расчета динамической вязкости газов в широком диапазоне давлений // Нефтегазовое дело, 2011, №1.

7. National Institute of Standards and Technology. Chemistry WebBook. URL: http://webbook.nist.gov/chemistry.

8. Рамазанов А. Ш., Акчурин Р. З. Моделирование распределения температуры в бурящейся скважине // Вестник Башкирского университета. 2016. Т. 21. №2. С. 269-273.

9. Исламов Д. Ф., Рамазанов А. Ш. Нестационарное температурное поле при фильтрации жидкости в неоднородном пласте // Вестник Башкирского университета. 2016. Т. 21. №1. С. 4-8.

Рис. 2. Сравнение распределений температуры (слева) и давления (справа) в стволе газовой скважины. 1 - модель с постоянными параметрами, 2 - модель с параметрами, зависящими от температуры и давления.

Поступила в редакцию 27.02.2017 г.

NUMERICAL MODELING OF TEMPERATURE AND PRESSURE IN A GAS WELL

© N. F. Fayzullin*, A. Sh. Ramazanov

Bashkir State University 32 Zaki Validi Street, 450076 Ufa, Republic of Bashkortostan, Russia.

Phone: +7 (347) 272 60 56.

*Email: dr.faizmeister@gmail. com

The authors of the article consider the problem of numerical modeling of stationary distribution of pressure and temperature in the gas well. During operation, heat exchange processes occurs in the well resulting in changing the distribution of temperature and pressure in the wellbore. Mathematical model that describes steady-state non-isothermal gas flow in the wellbore has been developed. The mathematical model is solved numerically by Euler formula in direction of flow. The simulator designed on the basis of the numerical model is able to simulate the distribution of temperature, pressure, and velocity of the gas in the wellbore taking into account the thermodynamic effects and convective heat transfer from the rocks. The program is tested for the case of the constant gas parameters and variable gas parameters depending on pressure and temperature. The model can be used to solve the inverse problem of determining the formation parameters and the flow according to the distribution of temperature, pressure, and gas flow rate in the wellbore.

Keywords: modeling, gas, well, pressure, heat exchange, heat transfer, convection.

Published in Russian. Do not hesitate to contact us at bulletin_bsu@mail.ru if you need translation of the article.

REFERENCES

1. Chekalyuk E. B. Termodinamika neftyanogo plasta [Thermodynamics of oil reservoir]. Moscow: Nedra, 1965.

2. Sui W. Determining multilayer formation properties from transient temperature and pressure measurements: PhD dissertation, Texas A&M University, 2009.

3. Hasan, A. R. and C. S. Kabir, 2002. Fluid Flow and Heat Transfer in Wellbores. Society of Petroleum Engineers, Richardson, Tex, pp: 181.

4. Sadretdinov A. A. Neizotermicheskaya fil'tratsiya szhimaemogo flyuida v sisteme skvazhina-plast: dis. kand. fiz.-mat. nauk. BashGU, Ufa, 2012.

5. Zotov G. A. Rukovodstvo po issledovaniyu skvazhin [Guidebook on studying oil wells]. Moscow: Nauka, 1995.

6. Glumov D. N., Strekalov A. V. Neftegazovoe delo, 2011, No. 1.

7. National Institute of Standards and Technology. Chemistry WebBook. URL: http://webbook.nist.gov/chemistry.

8. Ramazanov A. Sh., Akchurin R. Z. Vestnik Bashkirskogo universiteta. 2016. Vol. 21. No. 2. Pp. 269-273.

9. Islamov D. F., Ramazanov A. Sh. Vestnik Bashkirskogo universiteta. 2016. Vol. 21. No. 1. Pp. 4-8.

Received 27.02.2017.