CC BY
17
БОБРЕНЁВА Ю. О., ГУБАЙДУЛЛИН И. М., ЖАЛНИН Р. В.
НЕЯВНАЯ СХЕМА НА ОСНОВЕ РАЗРЫВНОГО МЕТОДА ГАЛЁРКИНА ДЛЯ МОДЕЛИРОВАНИЯ ТЕМПЕРАТУРНЫХ ПОЛЕЙ В НЕФТЯНОМ ПЛАСТЕ1
Аннотация. В работе для моделирования температурного поля в системе «скважина-трещина-пласт» предложена неявная схема, основанная на разрывном методе Галёркина. Рассматривается двумерное уравнение теплопроводности с постоянными коэффициентами, равными единице. Получены выражения для определения элементов матрицы системы линейных алгебраических уравнений.
Ключевые слова: разрывный метода Галёркина, неявная схема, термометрия, РМГ, DG.
BOBRENEVA YU. O., GUBAYDULLIN I. M., ZHALNIN R. V.
IMPLICIT SCHEME BASED ON DISCONTINUOUS GALERKIN METHOD FOR MODELING TEMPERATURE IN OIL RESERVOIR
Abstract. The implicit scheme is developed to simulate the temperature in the "well-fracture-layer" system. The scheme is based on the discontinuous Galerkin method. The two-dimensional heat conduction equation with constant coefficients equal to unity is considered. The expressions to determine matrix elements of the system of linear equations are obtained.
Keywords: discontinuous Galerkin method, implicit scheme, thermometry, DG.
При гидродинамических исследованиях скважин, особенно для скважин с гидравлическим разрывом пласта [1], необходимо совместное рассмотрение гидродинамического состояния системы «скважина-пласт» с температурным полем, или термометрией. Термометрия исторически является первым методом исследования скважин. Основным параметром, который несет информационную нагрузку в данном методе, является температура. Профиль температуры зависит от скорости фильтрации, градиента давления, а также свойств жидкости и породы [2; 3]. Для детального изучения термогидродинамических процессов в пласте со скважиной необходима разработка эффективного вычислительного алгоритма для численного моделирования совместной работы системы «скважина-трещина-пласт» [4]. Как известно, процессы теплопередачи описываются уравнениями параболического типа. При этом для использования алгоритмов сквозного счета необходимо строить сетки, имеющие достаточно мелкие ячейки в области, соответствующей трещине. Это накладывает существенные ограничения на шаг дискретизации по времени и, тем
1 Работа выполнена при финансовой поддержке РФФИ (грант 16-31-50018 мол_нр).
самым, требует большого числа итераций. Поэтому возникает необходимость построения неявных схем, для которых такое ограничение отсутствует. Рассмотрим уравнение: ди д2и д2и
= °,(х,у) е в, г > °,и = и^,х,у); (1)
д1 дх2 ду2
и = р(х,у),(х,у) едв,г> °;и\ = =^(х,у). (2)
Для дискретизации разрывным методом Галеркина запишем (1) в виде системы уравнений [5]:
ди дцх дцу & дх ду д и
д и
Покроем рассматриваемую область треугольной сеткой и введем в каждом треугольнике Д ^ систему базисных функций
{ р1т}с Рк(х,у), (4)
~ . , л (К + 2)(К + 1) т-^К /- N ТГ
т = °,..., N — 1, N =---, Рл (х, у) - пространство полиномов степени не выше К.
Решение системы (3) будем искать в виде:
N-1
и\х,у,г)=} щущ
1(х,уЛ) = ^ и\(Ор11(х,у),
1 = 0 N-1
Ч1х(х,у,0 = ^ Ч1Х1^)р11(х,у), (5)
=0 N-1
чУ(х,уЛ) = ^ ц1у1(Ор\(х,у). =0
Найдем проекцию уравнений (3) на пространство кусочно-непрерывных полиномов с базисом (4). Для этого подставим (5) в (3) и умножим скалярно каждое уравнение на базисную функцию:
Г ди1 Г дц1х дц1у
] -д^Ртйхйу + 1(— + —)(р}пйхйу = 0 А1 А1
д и
¡4Хср^йхйу + \—р^йхйу = 0,
А1 АI
д и
I 4у Рт йх йу + I рт йх йу = 0. А А д у
После преобразований получим
N-1
I Р11 Рт йх йу + + пу)рт йу —
1 = 0 АI дА1
N-1 £ N-1 £
^4x1 1 р\ддР-йхйу -^Ч1у1 1 Р11дрр^йхйу = 0,
дх " Аи^Ч ^ ду
=0 А =0 А
N-1 N-1
^ 4х1 I р\ рт йхйу + £ (и1 • пхрт йу — ^ и\ I р\ рт йхйу = 0,
^ 41у1 I р11 рт йхйу + £ (и1 • пурт йу — ^ и\ I р\ рт йхйу = 0.
1 = 0 А1 дАь 1 = 0 А[
Дискретизацию по времени выполним по неявной схеме. Значения на границе ячеек положим равным полусумме значений из ячеек, которые разделяет данное ребро. Окончательно получим:
N-1
—^ и I р\ рт йх йу +
1 = 0 а I
(N-1 N-1
^41х1Пх I р\ркйу + ^4х1Пх I Р1 Ртйу +
1 = 0 Г1к 1 = 0 Пк
N-1 N-1 ^
+ ^ 41у1Пу I р\ р1тйу + ^ 4уП I рI р1тйу\ —
1 = 0 Пк 1 = 0 Пк '
N-1 д I N-13 I ^N-1
— ^41Х1 | р\-р-йхйу — ^4у1 р11-рр^йхйу |Р¡РlmйхйУ,
=0 А =0 А =0 А
N-1
X 4lxi f Pi Pin dx dy + 1=0 м
+ X ^[X Ulnx f V ^dY + X и?Пх f pk VlmdY ) -
kenelgh(Af) \l=0 yik
N-1
-Iulf
d<pm
1 i i Pi^rdxdy = 0
1=0 Ai
N-1
Z4lfpipmdxdy+
1 = 0
(N-1 N-1
^ufoy f p\ pm.dY + X ufny f pk pkdY
1 = 0 yik l = 0 Yik
kenelgh(A{) \ 1=0 yik
N-1
v l f ldpm,, n
X U f pl——dxdy = 0.
1=0 д. У
(8)
Здесь neigh (Д/) - множество номеров соседних ячеек ячейки Д i; Yl? - ребро между ячейками с номерами i и к; и1 - значение на предыдущем слое по времени.
Таким образом, получили систему линейных алгебраических уравнений:
Ay = Ь, (9)
где матрица A - состоит из блоков размера 3N x 3N, столбец y состоит из блоков yl, а столбец Ь - из блоков bl, соответствующих i -ой ячейке сетки:
\
y =
/ и0 \
uNl -1
4x0
4xN-1 4y0
WyN-V
, ь1 =
N-1
(
1 = 0 д i
ll f p\ p0
d x d y
^ N-1
lXu fPllPN-1dxdy
1 = 0 д1
0 0
( ° )
где Ь = °, ...М — 1,М - количество ячеек расчетной сетки.
Для решения системы (9) можно воспользоваться библиотекой HYPRE [6,7]. Таким образом, предложен вычислительный алгоритм, основанный на неявной схеме для разрывного метода Галёркина, предназначенный для моделирования температурных
0
полей в системе «скважина-трещина-пласт». Неявная схема позволяет проводить расчеты с временным шагом, ограниченным только условием требуемой точности.
ЛИТЕРАТУРА
1. Малышев А. Г. и др. Анализ влияния технологических факторов и механических свойств горных пород на эффективность ГРП // Нефть Сургута. - М.: Нефтяное хозяйство, 1997. - С. 224-237.
2. Чекалюк Э. Б. Термодинамика нефтяного пласта. - М.: Недра, 1965. - 238 с.
3. Рамазанов А. Ш., Нагимов В. М. Аналитическая модель для расчета температурного поля в нефтяном пласте при нестационарном притоке жидкости // Электронный научный журнал «Нефтегазовое дело». - 2007. - № 1. - С. 1-9.
4. Губайдуллин И. М., Линд Ю. Б., Коледина К. Ф. Методология распараллеливания при решении многопараметрических обратных задач химической кинетики // Вычислительные методы и программирование: новые вычислительные технологии. -2012. - Т. 13, № 2 (26). - С. 28-36.
5. Жалнин Р. В., Масягин В. Ф., Панюшкина Е. Н. О применении разрывного метода Галеркина для численного решения двумерных уравнений диффузионного типа на неструктурированных разнесенных сетках // Современные проблемы науки и образования. - 2013. - № 6. - С. 874.
6. Falgout К. В., Yang U. M. Hypre: a Library of High Performance Preconditioners // Computational Science / P.M.A. Sloot, C.J.K. Tan. JJ Dongarra, A.G. Hoekstra (Eds.). -Springer-Verlag, 2002. - Vol. 2331. - P. 632-641.
7. Falgout R. D., Jones J. E., Yang U. M. The Design and Implementation of hypre, a Library of Parallel High Performance Preconditioners // Numerical Solution of Partial Differential Equations on Parallel Computers / A.M. Bruaset, A. Tveito (Eds.). - SpringerVerlag, 2006. - Vol. 51. - P. 267-294.
