CC BY
5
DOI - 10.32 743/UniTech.2022.98.5.13606
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ПЕРЕТОКА В МОДЕЛИ ФИЛЬТРАЦИИ УОРРЕНА-РУТА НА ОСНОВЕ РЕШЕНИЯ ОБРАТНОЙ ЗАДАЧЕ
Холияров Эркин Чоршанбиевич
канд. физ.-мат. наук, доц. кафедры Информатики и информационной технологии, Педагогический институт Термезского государственного университета,
Республики Узбекистан, г. Термез E-mail: e. kholiyarov@mail. ru
Шадманов Илхом Эргашевич
канд. техн. наук, доц. кафедры Высшей математики, Самаркандский институт экономика и сервиса, Республики Узбекистан, г. Самарканд E-mail: shadmanov_ilkhom@mail. ru
Шерпулатов Шерзод Шухратович
магистрант кафедры математическое моделирование, Самаркандский государственный университета, Республики Узбекистан, г. Самарканд E-mail: sherzodsherpulatov@mail.ru
DETERMINATION OF THE FLOW COEFFICIENT IN THE WARREN-ROOT FILTRATION MODEL BASED ON THE SOLUTION OF THE INVERSE PROBLEM
Erkin Kholiyarov
Cand. phys.-math. Sci., Associate Profeccor of the Department of Informatics and Information Technology, Pedagogical Institute of the Termez State University,
Uzbekistan, Termez
Ilkhom Shadmanov
Cand. of the tech. Sci., Associate Professor of the Department of Higher Mathematics,
Samarkand Institute of Economics and Service, Uzbekistan, Samarkand
Sherzod Sherpulatov
Master student of the Department of Mathematical Modeling,
Samarkand State University, Uzbekistan, Samarkand
АННОТАЦИЯ
В работе поставлена и численно решена обратная задача по определению коэффициента перетока в модели Уоррена-Рута для фильтрации однородной жидкости в трещиновато-пористых средах. Для решения задачи использованы методы идентификации первого порядка. Установлено, что коэффициент перетока при различных нулевых приближениях с невозмущенными исходными данными восстанавливаются достаточно хорошо при небольших числах итераций.
ABSTRACT
In this paper, the inverse problem of determining the flow coefficient in the Warren-Root model for filtering a homogeneous fluid in fractured-porous media is posed and numerically solved. To solve the problem, first-order identification methods were used. It has been established that the flow coefficient at various zero approximations with unperturbed initial data is restored quite well at a small number of iterations.
Библиографическое описание: Холияров Э.Ч., Шадманов И.Э., Шерпулатов Ш.Ш. ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ПЕРЕТОКА В МОДЕЛИ ФИЛЬТРАЦИИ УОРРЕНА-РУТА НА ОСНОВЕ РЕШЕНИЯ ОБРАТНОЙ ЗАДАЧЕ // Universum: технические науки : электрон. научн. журн. 2022. 5(98). URL: https://7universum.com/ru/tech/archive/item/13606
A UNiVERSUM:
№ 5 (98)_- * ♦ > г ■■■ - ._май. 2022 г.
Ключевые слова: метод идентификация, обратная задача, фильтрация, трещиновато-пористая среда, устойчивость решения.
Keywords: identification method, inverse problem, filtration, fractured-porous medium, solution stability.
Введение
Теория фильтрации однородных жидкостей в трещиновато-пористых средах (ТПС) рассмотрены в [3, 4]. Согласно по этой теории ТПС описывается как взаимодействующая среда, состоящая из системы трещин и пористых блоков. Уравнения движения и сохранения массы записываются независимо для каждой среды, т. е. в каждой точке вводятся две проницаемости, пористости, скорости фильтрации и два давления. Переток жидкости из одной среды в другую учитывается введением функции источника-стока в уравнениях сохранения массы. Предполагается, что пласт однороден, изотропен и течение в обеих средах происходит в пределах справедливости закона Дарси. Жидкость сла-босжимаема, обе среды - упругие, происходит обмен жидкостью между трещинами и пористыми блоками и масса протекающей из блоков в трещины жидкости q подчиняется соотношению
q = а0 Р (p2 - pl),
м
(1)
где а о - безразмерный коэффициент, зависящий от геометрических характеристик пористых блоков; Ро - плотность при первоначальном давлении Ро , Ц - динамическая вязкость жидкости, р1, р2 - давления в трещинах и пористых блоках,
соответственно .
При этих предпосылках уравнения фильтрации в одномерном случае принимают вид:
р*&=^ %+ао (p2 - p),
И dt м dx2 м 2 1
* dp2 k2 d2p2 а0 / ч
p2^7" =--;т2---(p2 -pi),
dt м dx м
(2)
где
p* =pci + moipf; p = p0f1+pf(Pi-Po)];
m
i = moi + pci (Pi- Po); v = -—, l =1,2
м dx
к - проницаемости, щ - пористости, шы -пористости при р = р , V - скорость фильтрации, - коэффициент сжимаемости жидкости, Рс1 -коэффициент сжимаемости среды, р - плотность жидкости, индекс I = 1 соответствует трещинам, I = 2 - пористым блокам.
Из обобщенной модели (2) в условиях щ << щ, рс1 << Рс2, к << к, получается упрощенная система уравнений
ап
ki d 2 Pi мР2 dx2 мР
+ ■
■(P2 - Pi) = °
2
(3)
р2 %+а0(P2 - Pi )=o. dt м
Уоррен и Рут [17] учитывали сжимаемости трещин, но пренебрегали движением жидкости в пористых блоках. В этих соображениях из (2) получается следующую систему уравнений
* dP1 к1 d2 P1 а o ( ч
+ ( P2 - P1)' dt м dx м
(4)
Р
dP2 dt
+
а
^(P2 - P1 ) = o.
м
В некоторых источниках система уравнений (4) называется «усеченной».
Модели (2), (3), (4) широкое применяется в анализе процессов разработки залежей нефти с трещинными и трещиновато-пористыми коллекторами [5, 9, 16].
Некоторый анализ моделей движения жидкостей в ТПС можно найти в работах [11, 16]. Постановка задач для упрощенных и «усеченных» систем уравнений (3), (4) имеет ряд особенностей, этими сведениями можно ознакомиться в [10, 11].
Некоторые коэффициентные обратные задачи фильтрации жидкости в пористых и трещиновато-пористых средах решались в [6, 8, 12, 13, 14].
В этой работе используя систему уравнений (4) решается коэффициентная обратная задача по определению параметра перетока а о . В качестве исходных данных для обратной задачи принимается решение прямой задачи в определенных точках области исследования. Следовательно, проводиться квазиреальный вычислительный эксперимент.
Постановка обратной задачи
Система уравнений фильтрации однородной жидкости в ТПС (4), называемой уравнениями Уор-рена-Рута, записываем в виде
^=А %+А (P2 - P1 )=o, dt мР* dx2 мР*
dP2
+ —Т(P2 - Pi) = o, o < x < L, o < t < T.
dt мР:
2
Для решения системы уравнений (5) задается началные и граничные условия в виде
Pi (0,х) = p2 (0,х) = p0, p0 = const, 0 < х < L,
А d^
ц дх
= v0 = const,
(6)
(7)
р (X,Ь) = р0, 0 < X < Т.
Задача (5), (6), (7) соответствует прямой постановке. Для решения обратной задачи по определению а0 необходимо задать дополнительные условия, в качестве которых используем решение задачи (5), (6), (7) при известном а0 в заданных точках области. Пусть такая информация дается в точке х = 0 как решение р (/, 0), которого обозначим как
г (X) .
Обратная задача ставится следующим образом: определить коэффициент перетока а0 из условия минимума следующего функционала
J(а0) = {[p (t,0)-z(t)]2 dt,
(8)
где р (X, 0) - решение задачи (5) - (7) при заданном
а0 .
Решение обратной задачи
Условие стационарности функционала (8) и (®о)/ Са о запишется в виде
и (а ) Т
-а = 21 [р1 (X, 0) - г (/)] щ (/,0) С = ^ (а0) = 0, (9)
са
др.
где щ = —1— функция чувствительности [1, 2] по
да0
коэффициенту а .
Разложим в ряд функцию р в окрестности а5 с точностью до членов второго порядка [2]
5+1 5 5
p (t, х) « Pi (t, х) + (аЦ+1 - а05) w, (t, х). (10)
где а5 - 5 -ое приблежение а0, р1, р1 - приближения, соответствующие а5 и а15+1.
Подставляем разложение (10) в соотношение (9) получим следующие соотношение:
T Г
J Pi(t,0) + (а5+1 -а5)W(t,0)-z(t)
щ (t,0) dt = 0,
из последнего выражения можно вычислить прибли-
5 +1 5 5
жение а 0 , если функции р1 (X ,0) и щ (X ,0) известны
J
а5+1 = 0 а0
а05 wi (t,0)- TP, (t,0) + z (t)
W,
\ (t ,0) dt
J Wi (t,0)
dt
(11)
Дифференцировав систему уравнений (5) по коэффициенту а получим следующую систему уравнений
dwj k д2 w а
dt [ß* дх [ß*
+ —(W2 - Wi) + (P2 - Pi ) = 0,
ß
dW2 + 00*(W2 - Wi ) + 4*(P2 - Pi ) = 0, 0 < х < L, 0 < t < T, dt [ß* [ß*
(12)
др др2 где щ =-, щ =- - функции чувствительно-
да0 да
сти [1, 2] по коэффициенту а.
Начальные и граничные условия для системы уравнений (12) получаются из условий (6), (7) дифференцированием по а :
1^(0, х) = w (0, х) = 0, 0 < х < L, (13)
= 0, w (t, L ) = 0, 0 < t < T. (14)
k dWj
ц дх
Численный алгоритм определения коэффициента а можно построить следующим образом:
1) Выбирая начальное приближение а 0 (полагая 5 = 0);
2) Решается задачу (5) - (7) от X = 0 до X = Т и определяется функцию р . Вычисляется значение функционала (8). После этого решается задачу (12) -(14) от X = 0 до X = Т и определяется функцию щ ;
3) По формуле (11) вычисляется приближение
аг
,5+i
4) Повторяется этапы 2), 3) до выполнения усло-
вия
J5+ - J5
J5
■<S,
а 5+i - а5 Гл0 0
а„
где Sj, е2 - допустимые погрешности.
2
х=0
0
0
х=0
5 5 + i
В рамках квазиреального эксперимента [15] сначала рассматривается прямая задача (5) - (7) с известным а™304 = 3,6 • 10-16. Эта задача решается численно методом конечных разностей [7]. По результатам численных расчетов определяется сеточная функция 23 = 2), 3 = 0,1,...,М. Кроме
того, при решении обратной задачи сеточная функция 2 (() зашумляется случайными погрешностями [15]
следующим образом:
■ 251 tfj --
(15)
где ст3 - случайная функция, равномерно распределенная на интервале [0,1], 8 - уровень
погрешности.
Для численного решения задачи (5) - (7) использованы следующие исходные значения параметров: Т = 2000 с, Ь = 60 м, к = 1-10-12 м2, ро = 10
МПа, ц = 2,5-10-8 МПа-с, 0* = 1-10-5 МПа-1,
V = 2 • 10-6 м/с. График ^ (/) представлен на рис.1.
Рисунок 1. График функции Zg (t)
Разностные задачи
Задачи (5) - (7), (12) - (14) при а0 = а^ решаются численно с помощью метода конечных разностей [7].
В области Б = {0 < х < Ь, 0 < t < Т} введем равномерную сетку
О =
[(х,.,tj), х,. = ih, i = 0,1,...,N, h = LN, tj = jz, j = 0,1,...,M, т = T/MJ;
где к - шаг сетки по координате х, т - шаг сетки по времени t. Вводим обозначения р13 = р (tJ, х),
Р23 = Р2 , X ) , ^ = (tj , X ) , 3 = (tj , X ) , соЦ = о (г, х), о 23 = а2 (tJ, х).
Для разностной аппроксимации задачи (5) - (7) используем чисто неявную разностную схему:
Pij+i - Pij
= *N>1{+i + а(P2j -Pij+i) = 0, [ß2
i = i,2,..., N -1, j = 0,i,...,M-i,
P 2j+i - P 2j
+ —^r (P2j+i - Pj ) = 0, i = 0,i,..., N, j = 0,i,..., M -1,
P1,0 = Po, P2,0 = Po, i = 0,1,...,N,
ki p1,+i - Pii+i
= Vo, PiN+1 = Po, j = 0,1,...,M -1,
(16)
где
APij = (Pij-i - 2Pij + Pi;f+1),
h [ß2
вид
к'
Функционал (8) после аппроксимации принимает
м-1 2
и = £*( рГ - г'+1) , ^ = г (г]).
] = 0
Аналогично аппроксимируем задачу (12) - (14)
Щ]+1 - Щ] . „ а'п
= %Awij+1 +а5*(w2j -wij+1 ) + -i-(P2j+1 -Pij+1) = 0,
w2 j+1 - w2 j
i = 1,2,...,N -1, j = 0,1,...,M -1,
+-ß (w2j+'-—4 (P2j+'-Pij+1 )=0
i = 0,1,...,N, j = 0,1,...,M-1, wi0 = 0, w2, = 0, i = 0,1,...,N, — w1j+1 - w1j+1 ,,
--1 —0-^ = 0, —1д+1 = 0, j = 0,1,...,M-1.
[i h
(17)
После аппроксимации приближение а5+1 (11) вычисляется следующим образом
а5+1 = а 0
m -i
Я—о—+1 - P10+1 + zj+1 ] —10+1
M -1 0
Я wij+1 ]2
j =0
Системы разностных уравнений (16), (17) решаются методом прогонки [7], и находятся разностные решения р1], р2], щ1], щ2].
Результаты численных расчетов
Сетка разбивала координатный отрезок [0,60] на 120 интервалов, временный отрезок [0,2000] -
на 4000 интервалов. «Данные измерения» (15)
подготовлены на основе этого решения в 200 точках «время». Результаты расчетов восстановления коэффициента а0 с невозмущенными исходными данными при различных начальных приближениях представлены на Рис. 2-3.
В случае, когда начальное приближение до три раза больше (или четыре раза меньше) чем точное значение искомого коэффициента, для восстановления параметра а0 достаточно 5-7 итераций (Рис.2.а.,
Рис.2.б.). В случае, когда начальное приближение в 3-5 раза больше чем точное значение искомого коэффициента, требуется 10-23 итераций (Рис.2.в., Рис.2.г.). А в случае, когда начальное приближение от пяти до шести раза больше чем точное значение искомого коэффициента, требуется 28-75 итераций (Рис.З.а., Рис.З.б.). Чем больше начальное приближение удаляется от точки равновесия, требуется тем большее число итерации.
Результаты расчетов с возмущенными исходными данными представлены в табл. 1. Численные расчеты проведены с начальным приближением а° = 7,2 -10-16. Относительные погрешности восстановления коэффициента а0 изменяется в пределах 0,000028 % до 7,912444 %. Как видно из табл. 1 относительная погрешность определения а с
увеличением погрешности исходных данных увеличивается.
г
h
т
т
I
♦
_I
5
Рисунок 2. Восстановление коэффициента а0 с невозмущенными исходными данными (при 8 = 0),
ап
■ использованное при подготовке исходных данных значение параметра а0 •
Рисунок 3. Восстановление коэффициента а0 с невозмущенными исходными данными (при 8 = 0),
а
Таблица 1.
Восстановление коэффициента ап с возмущенными исходными данными
S Ц0 = 7,2 • 10-16
s < , с а - aexact 0 0 Относительная погрешность--100, % А exact а0
0,0 28 3,599999 -10-16 0,000028
0,005 23 3,589153-10-16 0,301306
0,01 22 3,450903 -10-16 4,141583
0,02 24 3,405337 -10-16 5,407306
0,05 18 3,884848 -10-16 7,912444
Заключение
В данной работе рассмотрена обратная коэффициентная задача фильтрации однородной жидкости для модели Уоррена-Рута в трещиновато-пористых средах. Искомым параметрам коэффициент перетока находится из решения обратной задачи.
Для того чтобы подготовить дополнительную информацию для решения обратной задачи рассматривалась соответствующая прямая задача с известным значением коэффициента перетока. Таким образом подготавливается «исходные данные» для решения обратной задачи. Проводились также расчеты с возмущенными исходными данными, которые подготовлены путем зашумления данных случайными погрешностями.
Список литературы:
1. Алифанов О.М., Артюхин Е.А., Румянцев С.В. Экстремальные методы решения некорректных задач. - М.: Наука, 1988. - 288 с.
2. Бабе Г.Д., Бондарев Э.А., Воеводин А.Ф., Каниболотский М.А. Идентификация моделей гидравлики. -Новосибирск: Наука, 1980. - 161 с.
3. Баренблатт Г.И., Желтов Ю.П. Об основных уравнениях фильтрации однородных жидкостей в трещиноватых породах. // ДАН СССР, - 1960. - Т. 132, - №3. - С. 545-548.
4. Баренблатт Г.И., Желтов Ю.П., Кочина И.Н. Об основных представлениях теории фильтрации однородных жидкостей в трещиноватых породах // ПММ. - 1960. - Т. 24, вып. 5. - С. 852-864.
5. Майдебор В.Н. Особенности разработки нефтяных месторождений с трещиноватыми коллекторами. - М.: Недра, 1980. - 288 c.
6. Нармурадов Ч.Б., Холияров Э.Ч., Гуломкодиров К.А. Численное моделирование обратной задачи релаксационной фильтрации одно однородной жидкости в пористой среде // Проблемы вычислительной и прикладной математики. - 2017. - №2. - С. 12-19.
7. Самарский А.А. Теория разностных схем. - М.: Наука. 1989. - 616 с.
8. Хайруллин М.Х., Хисамов Р.С., Шамсиев М.Н., Фархуллин Р.Г. Интерпретация результатов гидродинамических исследований скважин методами регуляризации. М. - Ижевск: НИЦ «Регулярная и хаотическая динамика»; Институт компьютерных исследований, 2006. - 172 с.
9. Шаймуратов Р.В. Гидродинамика нефтяного трещиноватого пласта. - М.: Недра, 1980. - 223 с.
10. Barenblatt G.I., Entov V.M., Ryzhic V.M. Theory of Fluid Flows through Natural Rocks. London, 1990.
11. Chen Z.-X. Transient flow of slightly compressible fluids through double-porosity, double-permeability systems // Transport in Porous Media. 1989. Vol. 4. P. 147-184.
12. Khairullin M.H, Abdullin A.I., Morozov P.E., Shamsiev M.N. The numerical solution of the inverse problem for the deformable porous fractured reservoir // Matem. Mod. - 2008. Vol. 20. № 11, - P. 35-40.
13. Khuzhayorov B., Kholiyarov E. Inverse problems of elastoplastic filtration of liquid in a porous medium // Journal of Engineering Physics and Thermophysics. - 2007. Vol. 8. № 3. - P. 517-525.
Для решения обратной задачи использован метод идентификации первого порядка. Метод идентификации основан путем минимизации функционала невязки. Минимум функционала определяется из условия стационарности по искомому коэффициенту. Результаты расчетов показывают, что если в итерационном процессе начальные приближения близки к точке равновесия, коэффициент перетока восстанавливается достаточно быстро. При удалением начальные приближения от точки равновесия, при восстановление коэффициента перетока требуется больше итераций.
14. Khuzhayorov B., Ali Md. F., Sulaymonov F., Kholiyarov E. Inverse coefficient problem for mass transfer in two-zone cylindrical porous medium // AIP Conference Proceedings. - 2016. Vol. 1739. 020028.
15. Samarskii A.A., Vabishchevich P.N. Numerical Methods for Solving Inverse Problems of Mathematical Physics. Berlin: Walter de Gruyter, 2007. - 438 p.
16. Van Golf-Racht, T.D., Fundamentals of Fractured Reservoir Engineering. Developments in Petroleum Science, Elsevier Scientific, Amsterdam, Oxford, New York. 1982. - 365 p.
17. Warren, J.E. and Root, P.J. The behavior of naturally fractured reservoirs // Soc. Petrol. Eng.J. 1963, Sept., - P. 245-255.
