ПРОБЛЕМЫ НЕФТЕДОБЫЧИ, НЕФТЕХИМИИ, НЕФТЕПЕРЕРАБОТКИ И ПРИМЕНЕНИЯ НЕФТЕПРОДУКТОВ
УДК 517: 532, 536
М. Х. Хайруллин, М. Н. Шамсиев, П. Е. Морозов, А. И. Абдуллин, В. Р. Гадильшина, И. Т. Салимьянов
ЧИСЛЕННОЕ РЕШЕНИЕ ПРЯМЫХ И ОБРАТНЫХ ЗАДАЧ ТЕПЛОМАССОПЕРЕНОСА
В НЕФТЯНЫХ ПЛАСТАХ
Ключевые слова: неизотермическая фильтрация, термогидродинамические исследования, тепловое воздействие.
В данной работе строится математическая модель для описания термогидродинамических процессов в нефтяном пласте с учетом влияния ствола скважины. Предлагается численный метод к математической обработке результатов измерений, полученных при эксплуатации вертикальных нефтяных скважин, основанный на применении методов регуляризации. Предложена математическая модель процесса прогревания призабой-ной зоны скважины. Показано, что применение скважинного нагревателя при разработке высоковязких неф-тей приводит к увеличению производительности скважины.
Key words: _nonisothermal filtration thermohydrodynamic study, the thermal effects.
In this paper, a mathematical model to describe the thermodynamic processes in the oil reservoir with the influence of the well bore. The numerical approach to the mathematical treatment of the results of measurements obtained in the operation of vertical oil wells, based on the theory of inverse problems of mathematical physics. A mathematical model of the process of warming up well bottom zone . It is shown that the use of a downhole heater in the development of high-viscosity leads to an increase in productivity of the well.
Введение
Изучение температурных явлений (экспериментальные термодинамические исследования явления дросселирования), связанных с фильтрацией флюида в пористых средах были начаты Б.Б. Ла-пуком в 30-х годах прошлого столетия. При изучении термодинамических эффектов фильтрации рассматривается система дифференциальных уравнений с частными производными с соответствующими начальными и граничными условиями, которая выводится из законов сохранения, количества движения и энергии [2, 3].
Исследования термодинамических режимов работы скважин необходимы для оценки теплофи-зических и фильтрационно-емкостных параметров пластов. Пластовая температура наряду с давлением определяет термодинамические процессы в пласте. В работах Э.Б. Чекалюка [4] решается линеанеари-зованная задача неизотермической фильтрации. По стационарному распределению давления вычисляется скорость фильтрации, которая затем подставляется в уравнение энергии. На основе этого подхода был предложен метод термодинамического зондирования для оценки фильтрационных параметров пласта. При этом пренебрегаются коэффициентами теплопроводности и адиабатического расширения. Это предположение обосновывается тем, что скорость изменения давления гораздо больше скорости изменения температуры. Это действительно имеет место при фильтрации упругой жидкости в пластах с хорошими коллекторскими свойствами, но не обязательно верно при фильтрации с низкими коэффициентами подвижности.
1. Процесс тепломассопереноса в круговом нефтяном пласте, вскрытый вертикальной скважиной, описывается следующей системой уравнений:
в* д-Р - 1-±[ ИгЁР\, Г e(rc,Rk), t > 0. dt г дг {у дг ' Vc k'
с T -1АЬ rТ
°П ~dt ~ r ~дг i П Г ~дг
(1)
др
+ РтсрЦ— +
k др [ дТ др"] t „ Ч , „ + рср-+ ед.\, г е(г,c,Rk), t > 0;
н у дг i дг дг )
с начальными
р(г,0) - Ро(г), Т(г,0)- То(), г е ]fc,Rk) и граничными условиями
2яИ
к
М дг
др дг
- q + C-др- или р(с^)- р3
- 0 или Т([с,t)- Т3
(2)
(3)
(4)
(5)
= Рк, Т^к,^ =тк, (6)
где р = р(г^), Т = Т(г^) - давление и температура пласте, Н - толщина пласта, рк - пластовое давление, Тк - пластовая температура, рз - забойное давление, Тз - температура на забое скважины, р = р(Т) - зависимость вязкости от температуры, гс - радиус ствола скважины, - радиус контура питания, С - коэффициент влияния ствола скважины, р - упругоемкость пласта, сп = трОр + (1 - т)с3р3 - объемная теплоемкость пласта, Ср, с3 - удельная теплоемкость жидкости и
с
среды соответственно, р , ps - плотность жидкости и среды соответственно, m - пористость пласта, Яп = (l - m)s + ml - теплопроводность пласта, Я ,
ls - теплопроводность жидкости и среды соответственно, / - коэффициент адиабатического расширения, s - коэффициент Джоуля-Томсона. Для численного решения системы (1) - (6) применяется метод конечных разностей.
В первом примере предполагается, что у = const. Приводится сравнение с аналитическим решением, полученным Э.Б. Чекалюком.
Пример 1. Рассматривается пуск вертикальной скважины с постоянным дебитом. Расчеты проводились при следующих параметрах: cp =1800
Дж/кг К, , cs = 800 Дж/кг К, р = 800 кг/м3, ps =2700 кг/м3, m = 0.2, Я = 0.1385 Вт/м К, ls = 2.42 Вт/м К, / = 0.014 К/МПа, s = 0.4 К/МПа, Tk = 300 K, pk = 10 МПа, гс = 0.1 м, Rk = 100 м, H = 10 м, C = 0.7 м3/МПа, и = 5 мПа с, k = 0.01 мкм2, j* = 10-4 МПа-1, q = 5 м3/сут, Гэкс = 100 сут.
На рис. 1 приводится результаты сравнения изменения температуры, полученные из решения системы уравнений (1)-(5) и по аналитической формуле, предложенной Э.Б. Чекалюком [4].
0.8 —| л т.С
0.6 -
0.4 -
-TTTTTTTf
0.001 0.010
1.000
"I ........I
10.000 100.000
Рис. 1 - Изменение температуры на забое скважины. + - решение по формуле Э.Б. Чекалюка; • - численное решение
Кривые изменения температуры практически совпадают. Такое совпадение получено при условии, что на скважине мгновенно устанавливается дебит, т.е. не учитывается влияние ствола скважины, а также ) = 0 , Х = 0 . Температурное поле, связанное с дроссельным эффектом, при постоянном отборе совпадает с таким же полем для несжимаемой жидкости. Результаты расчетов с учетом влияния ствола скважины приведены на рис.2. При пуске скважины в начальный момент жидкость отбирается в основном из затрубного пространства и через некоторое время в работу включается пласт, что ска-
зывается на изменениях температуры на забое скважины. Это особенно сильно влияет на изменения забойной температуры для малодебитных скважин.
Результаты расчетов показывают, что при вариации коэффициентов теплоемкости и теплопроводности на 25% изменение температуры на забое скважины составляет около 5%. Коэффициент адиабатического расширения влияет только на начальный участок кривой. Наибольшое влияние на изменение кривой температуры имеют коэффициенты Джоуля - Томсона и проницаемости. Коэффициент адиабатического расширения практически не влияет на кривую изменения температуры.
0.8 —| а т\с
0.001
Рис. 2 - Изменение температуры на забое скважины (учитываются все эффекты). Численное решение с учетом (♦) и без учета (•) влияния ствола скважины
2. Естественное температурное поле вокруг действующей скважины бывает нарушено за счет теплообмена движущегося в ней потока жидкости с окружающей средой. После остановки скважины начинается процесс восстановления давления температуры. Он зависит от многих факторов, прежде всего от предыстории эксплуатации скважины. При постановке обратной задачи об интерпретации кривой восстановления температуры необходимо знать распределение температуры в пласте перед остановкой скважины. Поэтому в данной работе рассматривается задача интерпретации результатов термогидродинамических исследований после пуска скважины. В этом случае формирование происходит на фоне независящего от времени первоначального распределения параметров. Для этого продолжительность простоя скважины должно быть больше, чем продолжительность цикла измерений.
Обратная задача определения коэффициентов объемной теплоемкости пласта сп, теплопроводности пласта Лп, Джоуля-Томсона е и проводимости к сводится к минимизации функционала И
[1, 5]:
F(a)-JHf)-T(rc,t)dt, 0
где ф(() - наблюдаемые значения температуры, Т(rc,t) - вычисленные значения температуры на
забое
скважины,
где
a -I сп, Лп, e,k
М
0 < a; < а, < bj (a,, bj - const).
Итерационная последовательность для минимизации функционала-невязки (7) строится на основе метода Левенберга-Марквардта. Значения переменных минимизации на j-ой итерации вычисляются по формуле: aj+1 - aj - (hj + ve) VFj, где H - приближенная матрица вторых производных, H - АТA , A - матрица чувствительности, v - параметр регуляризации, VF - градиент функционала-невязки.
Критериями остановки итерационного процесса служит выполнение хотя бы одного из неравенств:
F а+1)- f а)
< si,
а +1 - а' \ <s2 .
Модельные кривые изменения температуры, приведенные в примере 1, используются в качестве исходной информации при решении обратной задачи. Для исследования устойчивости в модельную кривую изменения температуры вводили случайным образом погрешности. Результаты интерпретации возмущенной кривой приведены на рис. 3. Из результатов видно, что предложенный алгоритм устойчив относительно погрешности входной информации.
300.8 -, AT, С
300.6 -
300.4 -
300.2 -
300.0
I, сут
299.8 --1||| ии|—1||| ии|—1||| ии|—г"гтттттт|—г"гтттттт|
0.001 0.010 0.100 1.000 10.000 100.000 Рис. 3 - • - возмущенная кривая,--вычисленная кривая температуры
Предложенный метод интерпретации результатов термогидродинамических исследований вертикальных скважин позволяет оценивать фильтрационные и теплофизические параметры нефтяного пласта.
3. Тепловые методы воздействия на пласт для увеличения нефтеотдачи, помимо методов за-
воднения, рассматриваются как единственная реализуемая на промышленном уровне альтернатива.
Нагревание ствола скважины и призабойной зоны пласта применяется для уменьшения вязкости нефти и облегчения процесса ее извлечения, а также для расплавления и удаления парафиносмолистых отложений.
В этом пункте предложена математическая модель процесса прогревания призабойной зоны скважины. Исследуется скорость распространения теплового фронта в призабойной зоне в зависимости от мощности нагревателя и проницаемости пласта. Показано, что применение скважинного нагревателя при разработке высоковязких нефтей приводит к увеличению производительности скважины. Пример 2: Рассматривается пуск вертикальной скважины с постоянным забойным давлением при постоянном прогреве призабойной зоны. Предполагается, что вязкость зависит от температуры. Расчеты проводились при следующих параметрах: Ср =1800
Дж/кг К, , С8 = 800 Дж/кг К, р = 996 кг/м3, р3 =2700 кг/м3, т = 0.2, 2 = 0.1385 Вт/м К, Лв = 2.42 Вт/м К, 77 = 0.014 К/МПа, £ = 0.4 К/МПа, Тк = 293 К, Тз = 320 К, рк = 10 МПа, рз = 6 МПа, гс = 0.1 м, Ик = 200 м, Н = 10 м, к = 0.1 мкм2, р* = 10-4 МПа-1. Экспериментальная зависимость вязкости от температуры приведена в табл.1. Результаты расчетов приведены на рис. 4 - 5.
Таблица 1 - Экспериментальная зависимость вязкости от температуры
T, K М, мПа-с
293 340
303 210
323 120
343 80
340 - \ T, K p, МПа / Г ю
330 - - 9
320 - \T Р/ - 8
310 - -
300 - - 7
290 - ....... ! ! - 6
0.1
10
100 r, м
Рис. 4 - Распределение температуры и давления в пласте. Т3 = 340 K
Рис. 5 - Изменение дебита скважины по времени. • - температура прогрева 320 К, ■ - температура прогрева 340 К
На рис. 4 приводятся распределения температуры и давления в пласте на момент времени 7=360 сут. Распределение давления в полулогарифмических координатах имеет нелинейный характер. Радиус прогретой зоны около 10 м.
Изменение дебита скважины по времени в зависимости от температуры нагревателя приводится на рис. 5. Для данного примера, после 30 сут прогрева наблюдается рост дебита скважины.
Вывод
Построена математическая модель для описания термодинамических процессов в нефтяном пласте с учетом влияния ствола скважины. Разработан метод интерпретации результатов термогидродинамических исследований вертикальных скважин на основе методов регуляризации. Предложена математическая модель процесса прогревания призабойной зоны скважины. Показано, что применение скважинного нагревателя при разработке высоковязких нефтей приводит к увеличению производительности скважины.
Работа выполнена при поддержке гранта РФФИ №12-01-97005-р_поволжье_а.
Литература
1. A.I. Abdullin, M.N. Shamsiev, E.R. Badertdinova, M.H. Khairullin. Interpretation of temperature curves measured in a vertical well // XVI International Conference on the methods of aerophysical research, Kazan, 19-25 August 2012. -Composing, ITAM SB RAS. P.5-6.
2. Басниев К.С., Дмитриев Н.М., Каневская Р.Д., Максимов В.М. Подземная гидромеханика. М. - Ижевск: Институт компьютерных исследований, 2006. - 488 с.
3. Закиров С. Н., Индрупский И. М., Закиров Э. С., Аникеев Д. П. Новый подход к исследованию скважин и пластов // Нефтяное хозяйство, 2002. № 6, С. 113-116.
4. Чекалюк Э.Б. Термодинамика нефтяного пласта. М.: Недра, 1965. - 238 с.
5. Шамсиев М.Н., Бадертдинова Е.Р. Оценка фильтрационных и теплофизических параметров нефтяного пласта по результатам измерений температуры на забое скважины // Измерительная техника, 2012. №3. С. 45-47.
© М. Х. Хайруллин - д.т.н., проф., зав. лаб. подземной гидродинамики Института механики и машиностроения КазНЦ РАН, khairullin@mail.knc.ru; М. Н. Шамсиев - д.т.н., вед. науч. сотр. Института механики и машиностроения КазНЦ РАН, МБИаш-siev@yandex.ru; П. Е. Морозов - к.т.н., ст. науч. сотр. Института механики и машиностроения КазНЦ РАН, Morozov@mail.knc.ru; А. И. Абдуллин - к. ф.-м. н., ст. науч. сотр. Института механики и машиностроения КазНЦ РАН, adel777@rambler.ru; В. Р. Гадильшина - науч. сотр. Института механики и машиностроения КазНЦ РАН; И. Т. Салимьянов - к.т.н., доц. каф. информатики и прикладной математики КНИТУ, inisal@yandex.ru.