CC BY
16
УДК 532.546
МНОГОДАТЧИКОВЫЕ ТЕХНОЛОГИИ ИССЛЕДОВАНИЯ НЕФТЯНЫХ ГОРИЗОНТАЛЬНЫХ СКВАЖИН НА МЕСТОРОЖДЕНИЯХ РЕСПУБЛИКИ ТАТАРСТАН
М.Х. Хайруллин1, Е.Р. Бадертдинова2*, Н.А. Назимов3
'Институт механики и машиностроения Казанского научного центра РАН, Казань, Россия 2Казанский государственный технологический университет, Казань, Россия 3ПАО «Татнефть», Альметьевск, Россия
В Республике Татарстан в разработку широко вовлекаются трудноизвлекаемые запасы углеводородов, приуроченные к низкопроницаемым, неоднородным и расчлененным коллекторам. Повышение эффективности выработки трудноизвлекаемых запасов в значительной степени зависит от оперативного регулирования системы разработки на основе сведений о фильтрационных и теплофизических свойствах нефтяного пласта. Проблемы, связанные с интерпретацией геолого-промысловой информации, приводят к некорректным, в смысле Адамара, математическим задачам. Численное решение таких задач требует разработки специальных методов. В работе строится математическая модель термогидродинамических процессов, происходящих в системе «пласт-горизонтальная скважина». На основе этой модели и методов регуляризации А. Н. Тихонова предлагается вычислительный алгоритм для интерпретации результатов термогидродинамических исследований горизонтальных скважин и пластов. В качестве исходной информации берутся кривые изменения температуры, снятые одновременно несколькими глубинными приборами, установленными на разных участках горизонтальной части ствола скважины. Этот подход позволяет оценить фильтрационные параметры неоднородного пласта и строить профиль притока вдоль ствола горизонтальной скважины.
Ключевые слова: горизонтальная скважина, давление, температура, проницаемость, многодат-чиковые технологии, обратная задача
DOI: https://doi.Org/10.18599/grs.19.3.14
Для цитирования: Хайруллин М.Х., Бадертдинова Е.Р., Назимов Н.А. Многодатчиковые технологии исследования нефтяных горизонтальных скважин на месторождениях Республики Татарстан. Георесурсы. 2017. Т. 19. № 3. Ч. 1. С. 234-238. DOI: https://doi.Org/10.18599/grs.19.3.14
Проведение и обработка результатов термогидродинамических исследований при разработке месторождений с трудноизвлекаемыми запасами углеводородов, как правило, сопряжены со значительными трудностями. К ним относятся: механизированность фонда добывающих скважин, затрудняющая доставку глубинных измерительных приборов на забои скважин; низкие дебиты, приводящие к малой информативности дебитометрических исследований; термодинамические процессы, сопровождающиеся небольшими изменениями температуры; большая продолжительность гидродинамических исследований, затрудняющая использование классических методов интерпретации кривых восстановления давления.
Измерения температуры и давления в стволе горизонтальной скважины на основе многодатчиковой технологии (Фархуллин и др., 2003; Хайруллин и др., 2006) дают достаточно полную информацию о термогидродинамических процессах, происходящих в пласте и стволе. Изменение температуры в стволе горизонтальной скважины является интегральным показателем процессов тепломассопереноса, происходящих как в самой скважине, так и в пласте.
Будем считать, что ствол горизонтальной скважины параллелен кровле и подошве пласта, движение жидкости
*Ответственный автор: Бадертдинова Елена Радитовна E-mail: badertdinova@yandex. ru
в стволе - одномерное. Процесс распределения давления в стволе является квазистационарным, приток флюида к стволу при пуске скважины - радиальным. При этих предположениях из законов сохранения массы, импульса и энергии следует (Васильев, Воеводин, 1968; Чарный, 1975):
8v дх
2 w
, W--
k др2
ц дг
0<x<L,
(1)
dp, d(v 2) vi/ i i — = P^— + — PW , 0 <x< L, ox ox 4 r -
(2)
dTx
+ v
8pl
—L + e —
дх дх
2(a.
■ woC ) H p-( T,
P c.re
- Tr),
0 <x< L, 0 <t < t
(3)
р«ф2 = 1 д dt rdr
k_rd[>2 M 8r
0 <x< L, r < r <R , 0 <t < t ,
— c k — exp'
k dp2 f 8T2 dp.
+ e
dt ' Р Ц дг v дг дг
0 <х< L, r < r <R , 0 <t < t ,
— 7 c k — exp'
с начальными условиями:
(4)
(5)
p2(x, r, 0)=p0(x, r), T2(x, r, 0)=T0(x, r), 0<x<L, rc <r<R,, (6)
Í—^ПМНШ SCIENTIFIC AND TECHNICAL JOURNAL
^АШ GEDRESURSY
и граничными условиями:
= 0 </ < / .1 ^ ~ *
■ Ц дг
p(x, Rv t)=pv T2(x, Rv t)=Tk.
(7)
(8)
Здесь р1=р1(х), Т1=Т1(х, /) - давление и температура в стволе горизонтальной скважины, р2=р2(х, г, /), Т2=Т2(х, г, /), - давление и температура в пласте, рк - пластовое давление, Тк - пластовая температура, q - дебит горизонтальной скважины, - поверхность ствола горизонтальной скважины, гс - радиус ствола скважины, Rk - радиус контура питания, в* - упругоемкость пласта, у(х) - скорость флюида в стволе ГС, р - плотность флюида, е - коэффициент Джоуля-Томсона, у - коэффициент гидравлического сопротивления, а - коэффициент теплопередачи ствола горизонтальной скважины, Ср- удельная теплоемкость флюида, w - скорость фильтрации, Ь - длина ствола горизонтальной скважины, /ехр - время работы скважины.
Метод решения краевой задачи (1)-(8) основывается на сопряжении внешней (пласт) и внутренней (ствол горизонтальной скважины) задач. Система (1)-(8) решается численно методом конечных разностей. Область фильтрации покрывается неравномерной сеткой, которая сгущается к скважине. Полученная нелинейная система разностных уравнений решается итерационно.
На модельной задаче с данными, соответствующими реальным месторождениям Республики Татарстан, исследуются термогидродинамические процессы, происходящие в системе «пласт-горизонтальная скважина».
Рассматривается модельный нефтяной пласт, который разрабатывается горизонтальной скважиной. Горизонтальная скважина запускается в эксплуатацию с постоянным отбором жидкости из пласта. Исходные данные: Сп = 1.48 106Дж/(м3К), Ср = 1929 Дж/(кг К), Тк = 300 К, рк = 15 МПа, в* = 10-4МПа-\ ц = 25 мПа с, р = 800 кг/м3, е = 0.4 К/МПа, Ь = 100 м, г = 0.1 м, Rl = 5 м,
' ' 'с ' к '
q = 30 м3/сут. Предполагается, что проницаемость пласта является кусочно-постоянной функцией. В каждой зоне однородности пласта расположен глубинный прибор (Рис. 1). Рассматриваются следующие варианты зон неоднородности пласта:
к>к2, к1 = 0.05 мкм2, к2 = 0.01 мкм2, к<к2, к1 = 0.01 мкм2, к2 = 0.05 мкм2. В ствол горизонтальной скважины жидкость поступает из зон неоднородности нефтяного пласта с различной температурой за счёт эффекта Джоуля-Томсона (Рис. 1). Изменение температуры в стволе горизонтальной скважины о бъясняется калориметрическим эффектом.
Результаты численных расчетов показали, что для однородного нефтяного пласта температура по стволу является постоянной и со временем возрастает. Скорость жидкости по стволу горизонтальной
kl Пр.1 k2 Пр.2
1 <=> =» 1
50 м 50 м
скважины изменяется линейно. На рис. 2-3 представлены результаты расчетов изменения температуры и скорости потока для вариантов 1, 2 на момент времени t = 120 ч. Температура (Рис. 2), скорость потока (Рис. 3) в стволе горизонтальной скважины меняются нелинейно. Интенсивность притока флюида к стволу горизонтальной скважины (Рис. 4-5) имеет разрыв в точке, соответствующей границе зон однородности пласта по проницаемости. Распределение температуры, скорости потока по стволу горизонтальной скважины и интенсивность притока флюида к стволу горизонтальной скважины зависят от значений проницаемостей и размеров зон однородностей.
Отличительной чертой обратных задач нефтегазовой гидромеханики, связанных с исследованием математических моделей реальных процессов фильтрации в нефтяных пластах, является то, что характер дополнительной информации определяется возможностями промыслового эксперимента. На основе предложенной математической модели и методов регуляризации А.Н. Тихонова предлагается вычислительный алгоритм для интерпретации результатов термогидродинамических исследований горизонтальных скважин. Результаты измерений изменения температуры в период пуска скважины на разных участках горизонтального ствола скважины используются в качестве исходной информации. Участки размещения глубинной измерительной аппаратуры определяются на основе геофизических исследований скважины. В работах (Фархуллин и др., 2003; Хайруллин и др., 2006) описана технология проведения термогидродинамических
300.44
300.40
300.36
300.32
300.28
T, C
300.24
L, м
Рис. 1. Расположения приборов в неоднородном пласте
0 20 40 60 80 100 Рис. 2. Распределение температуры по стволу ГаС. ■ - Вар.1 к1>к2, ♦ - Вар. 2 k1<k2
0-012-1 v, м/с
0.008
0.004
0.000
0 20 40 60 80 100 Рис. 3. Распределение скорости по стволу ГС. ■ - Вар.1 k>k2, ♦ - Вар. 2 k<k2
НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ
ГЕОРЕСУРСЫ
0.5
0.4-
0.3
0.3
0.2
0.1
q м3/(сут м)
L, м
0 20 40 60 80 100 Рис. 5. Распределение притока по стволу ГС- Вар. 2 k¡<k2
0.4
0.3
0.3
0.2
0.1
Ь, м
1 | I | I | I | I |
0 20 40 60 80 100 Рис. 4. Распределение притока по стволу ГС.ш - Вар.1 k>k2
исследований горизонтальной скважины с помощью нескольких глубинных автономных приборов.
Пусть в местах расположения глубинных измерительных приборов в стволе горизонтальной скважины с координатами х ,, i=l,N, снимаются кривые изменения температуры:
T, (t)= T(xp t)=^(t), i=\,N, 0 <t < te,
(9)
Обратная коэффициентная задача формулируется следующим образом: определить коэффициент проницаемости k = к(х,г), когда термогидродинамические процессы в нефтяном пласте и стволе горизонтальной скважины описываются уравнениями (1)-(8). В качестве исходной информации используются замеренные значения температуры глубинными автономными приборами.
Оценка коэффициента проницаемости ищется в классе кусочно-постоянных функций к(х,г) = к , (х,г)е V ,
N _
-V, где Уп, п=\,Ы - области однородности (Рис. 1).
и=1
Приближенное решение обратной коэффициентной задачи (1)-(9) ищется из минимизации среднеквадратичного отклонения между наблюдаемыми и вычисленными величинами:
(10)
и=1 о
где ф() - наблюдаемые значения температуры, Т1 и(0 - вычисленные значения температуры, полученные из численного решения уравнений (1)-(8), а = (к1, к2, ..., к- искомый
параметр, 0 < m< k < M (m, M = const).
A n — n — nvn'n '
Исследование ГС № 18326 (Назимов, 2007). Скважина расположена в залежи N° 665 Ромашкинского месторождения Республики Татарстан. Скважина имеет 313 метровый открытый горизонтальный участок в данково-лебедянском горизонте в интервале от 1475 до 1788 м. В 2004 году в скважине проведены глубинные термогидродинамические исследования на основе многодатчиковых технологий. На рис. 6 указаны места расположения приборов. После окончания подземного ремонта скважина была пущена в работу с дебитом 7,8 м3/сут.
Предложенный вычислительный алгоритм применяется для интерпретации кривых изменения температуры, снятых глубинными приборами №№ 1879, 1721, 1726 и 1885. Для этого пласт разбивается на четыре зоны однородности, в которых расположены приборы .№№ 1879, 1721, 1726 и 1885.
На рисунках 7-8 приводятся результаты расчетов распределения температуры и скорости течения жидкости в
V" 1 „' . .=. h.....: Ь LTl : я г fe я ? С £ я
1 ÍÍÍ3
ms ?s
-ÍSJS
j
-síi-í
гт^убина. м
Рис. 6. Схема траектории ствола ГС № 18326 и места расположения приборов
Т, с
26.760
26.755
26.750
26.745
26.740
0 100 200 300
Рис. 7. Распределение температуры по стволу скважины
0.003 п v, м/с
0.002
0.001
0.000
0 100 200 300
Рис. 8. Распределение скорости жидкости по стволу скважины
0.5 -п q , м3/(сут м)
SCIENTIFIC AND TECHNICAL JOURNAL
GEDRESURSY
стволе горизонтальной скважины. Распределение притока жидкости вдоль ствола горизонтальной скважины приводится на рис. 9.
На рисунках 10-12 и в таблице приводятся результаты интерпретации кривых изменения температуры по приборам №№ 1879, 1721 и 1885. На рис. 10-12 приводятся вычисленные и наблюдаемые кривые изменения температуры. В таблице приводятся оценки проницаемости в зонах расположения приборов.
Из полученных результатов следует, что зоны расположения приборов № 1721, 1726 имеют низкую проницаемость - в этих участках приток к стволу горизонтальной скважины наименьший (Рис. 9). Геофизические исследования показали, что в районе 1680-1721 м (прибор № 1721) ствол скважины проходит по низкопроницаемому включению, а в районе 1620-1670 м (прибор №2 1726) - по слабонефтенасыщенному включению.
Результаты исследований, проведенных ОАО «Перм-нефтегеофизика» в середине 2006 года с технологическим комплексом «ЛАТЕРАЛЬ-2005», показали, что в зонах расположения этих приборов наблюдается низкий приток к стволу горизонтальной скважины.
Оценки проводимости в районах расположения приборов №№ 1885 и 1879, полученные по кривым изменения температуры и давления хорошо согласуются. При интерпретации результатов гидродинамических исследований
26.76 — т к
26.77
26.76 -
26.75
26.74 -
26.73
26.72
T, K
t, час
0 1 10 100 Рис. 9. Распределение притока жидкости по стволу скважины 0.05-, q*, м 3/(сут м)
0.04
0.03
0.02
0.01
0.00
L, м
100
200
300
Рис. 10. Прибор №1879. Кривая изменения температуры. • наблюдаемая, ♦ - вычисленная
26.75 -
26.74 -
26.73 -
26.72
♦-♦ ♦♦♦
26.71
0
1
10
t, час
мм
100
Рис. 11. Прибор №1721. Кривая изменения температуры. • наблюдаемая, ♦ - вычисленная
26.76 -и т К
26.75 -
26.74 -
26.73
26.72
10
t, час
мм
100
Рис. 12. Прибор №1885. Кривая изменения температуры. • -наблюдаемая, ♦ - вычисленная
Зона пр. №1879 Зона пр. №1721 Зона пр. №1726 Зона пр. №1885
Интерпретация кривых изменения температуры
k/\í (мю^/мПа-с) 1.04-10"3 2.13-10^ 1.93-КГ1 1.07-10"3
Интерпретация кривых изменения давления
к/Ц (мкм'/мПа-с) 3,42-Ю"3 4,46-Ю"3 3,6310"3 6,3410"3
Таблица. ГС № 18326. Оценки фильтрационных параметров.
скин-эффект не учитывался, поэтому имеется расхождение оценок проводимости в районах расположения приборов №№ 1721, 1726.
Литература
Васильев О.Ф., Воеводин А.Ф. О газотермодинамическом расчете потоков в простых и сложных трубопроводах (постановка задачи). Изв. Сиб. отд. АН СССР. 1968. №13. Вып 3. С. 53.
Назимов Н.А. Особенности характера течения флюидов в горизонтальных скважинах по данным глубинных исследований. Дисс. канд. тех. наук. Бугульма. 2007. 162 с.
Фархуллин Р.Г., Никашев О.А., Хайруллин М.Х., Муслимов Р.Х., Хисамов Р.С., Полушин В.И. Способ гидродинамических исследований горизонтальных скважин. Патент на изобретение RUS 2243372. 2003.
0
НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ OTWk^^Hi
ШШШШшжЬ'
Хайруллин М.Х., Хисамов Р.С., Шамсиев М.Н., Фархуллин Р.Г. Интерпретация результатов гидродинамических исследований скважин методами регуляризации. М.-Ижевск: НИЦ «Регулярная и хаотическая динамика»; Институт компьютерных исследований. 2006. 172 с.
Чарный И.А. Неустановившееся движение реальной жидкости в трубах. М.: Недра. 1975. 296 с.
Сведения об авторах
Мухамед Хильмиевич Хайруллин - доктор тех. наук, профессор, заведующий лабораторией подземной гидродинамики, Институт механики и машиностроения Казанского научного центра РАН
Россия, 420111, Казань, ул. Лобачевского, 2/31
Елена Радитовна Бадертдинова - доктор тех. наук, профессор кафедры информатики и прикладной математики, Казанский государственный технологический университет
Россия, 420015, Казань, ул. К. Маркса, д. 68
Нафис Анасович Назимов - канд. тех. наук, начальник службы геологии и опытно-промышленных работ по трудноизвлекаемым запасам, ПАО «Татнефть» Россия, 423450, Альметьевск, ул. Ленина, д. 75
Статья поступила в редакцию 17.06.2017;
Принята к публикации 27.07.2017; Опубликована 30.08.2017
IN ENGLISH
Multisensor Research Technologies of Oil Horizontal Wells on Fields of the Republic of Tatarstan
M.Kh. Khairullin1, E.R Badertdinova2, N.A. Nazimov3
'Institute of Mechanics and Engineering, Kazan Science Center, Russian Academy of Sciences, Kazan, Russia 2Kazan State Technological University, Kazan, Russia 3PJSC Tatneft, Almetyevsk, Russia
*Corresponding author: Elena R. Badertdinova, e-mail: [email protected]
Abstract. In the Republic of Tatarstan the development widely involves hard-to-recover reserves of hydrocarbons, confined to low-permeable, heterogeneous and dissected reservoirs. Improvement of recovery efficiency of such reserves largely depends on the operational control system based on information about the filtration and thermal properties of the oil reservoir. Issues, related to the interpretation of geological field information, lead to incorrect mathematical problems in terms of Hadamard. The numerical solution of such problems requires the development of special methods. A mathematical model of thermohydrodynamic processes occurring in the 'horizontal well' system is constructed in the paper.
Based on this model and regularization methods of A.N. Tikhonov, a computational algorithm is proposed for interpreting the results of thermohydrodynamic studies of horizontal wells and layers. Curves of the temperature changes are taken as the initial information, taken simultaneously by several deep instruments installed at different parts of the horizontal part of the wellbore. This approach makes it possible to evaluate the filtration parameters of an inhomogeneous reservoir and to construct an inflow profile along the trunk of a horizontal well.
Key words: horizontal well, pressure, temperature, permeability, multi-sensor technologies, inverse task
For citation: Khairullin M.Kh., Badertdinova E.R., Nazimov N.A. Multisensor Research Technologies of Oil Horizontal Wells on Fields of the Republic of Tatarstan.
Georesursy = Georesources. 2017. V. 19. No. 3. Part 1. Pp. 234-238. DOI: https://doi.org/10.18599/grs.19.3.14
References
Charny I.A. Unsteady flow of a real kuids in pipes. Moscow: Nedra. 1975. 296 p. (In Russ.)
Farkhullin R.G., Nikashev O.A., Khairullin, M.Kh. Muslimov R. Kh., Khisamov R. S., Polushin, V. I. Method of hydrodynamic investigations of horizontal wells. Patent RF 2243372. 2003. (In Russ.)
Khairullin M.Kh., Khisamov R.S., Shamsiev M.N., Farkhullin R.G. Interpretation of the results of well hydrodynamic research by regularization methods. Moscow-Izhevsk: SRC "Regular and chaotic dynamics"; Institute of Computer Science. 2006. 172 p. (In Russ.)
Nazimov N.A. Features of flow of fluids in horizontal wellsusing deep investigations. Avtoref. Diss. kand. tech. nauk [Abstract Cand. engineer. sci. diss.]. Bugulma. 2007. 162 p. (In Russ.)
Vasiliev O.F., Voevodin A.F. On the gas-thermodynamic calculation of flows in simple and complex pipelines (statement of the problem). Izv. Sib. Otd. AN USSR. 1968. V. 3. No. 13. Pp. 53. (In Russ.)
About the Authors
MuhamedKh. Khairullin - DSc (Engineering), Professor, Head of the Laboratory of Underground Hydrodynamics, Institute of Mechanics and Engineering, Kazan Science Center, Russian Academy of Sciences
420111, Russia, Kazan, Lobachevsky str, 2/31
Elena R. Badertdinova - DSc (Engineering), Professor, Department of Informatics and Applied Mathematics
Kazan State Technological University
420015, Russia, Kazan, K. Marks str, 68
Nafis A. Nazimov - PhD (Engineering), Head of the Service of Geology and Pilot Projects for Hard to Recover Reserves
PJSC Tatneft
423400, Russia, Almetyevsk, Lenin str., 75
Manuscript received 17 June 2017; Accepted 27 July 2017;
Published 30 August 2017
SCIENTIFIC AND TECHNICAL JOURNAL
GEDRESURSY
