CC BY
3
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
DOI -10.32743/UniTech.2023.117.12.16592
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ВРЕМЕНИ ЗАПАЗДЫВАНИЯ ФИЛЬТРАЦИИ ЖИДКОСТИ В НЕОДНОРОДНЫХ СРЕДАХ
Эрназаров Мирзохид Юлдаш угли
преподаватель,
Термезский университет экономики и сервиса, Республика Узбекистан, г. Термез E-mail: ernazarov. mirzohid@mail.ru
Холияров Эркин Чоршанбиевич
канд. физ.-мат. наук, доц., Термезский университет экономики и сервиса, Республика Узбекистан, г. Термез E-mail: e.kholiyarov@mail.ru
Тойиров Акбар Хасанович
канд. физ.-мат. наук, Термезский университет экономики и сервиса, Республика Узбекистан, г. Термез E-mail: akbartoyirov@gmail.com
Тураев Дилмурод Шокир угли
преподаватель,
Термезский университет экономики и сервиса, Республика Узбекистан, Сурхандарьинская область, г. Термез
E-mail: turaevdilmurod8@gmail.com
DETERMINATION OF THE DELAY TIME COEFFICIENT OF LIQUID FILTRATION
IN HETEROGENEOUS ENVIRONMENTS
Mirzokhid Ernazarov
Teacher,
Termez University of Economics and Service, Uzbekistan, Termez
Erkin Kholiyarov
Ph.D. physics and mathematics Sciences, Associate Professor, Termez University of Economics and Service, Uzbekistan, Termez
Akbar Toirov
Ph.D. physics and mathematics Sciences, Termez University of Economics and Service, Uzbekistan, Termez
Dilmurod Turaev
Teacher,
Termez University of Economics and Service, Uzbekistan, Termez
АННОТАЦИЯ
В работе поставлена и численно решена обратная задача по определению коэффициента запаздывания времени фильтрации жидкости в неоднородных пористых средах. Для решения задачи были использованы численные методы первого и второго порядка. Решение обратной задачи сводится к задаче нахождения минимума функционала, т.е. к задаче безусловной оптимизации. По дополнительным невозмущенным исходным данным коэффициент практически восстанавливается за три-четыре итерации в различных начальных приближениях. Показано, что метод первого порядка точнее, чем метод второго порядка. По этой причине для решения обратной задачи с возмущенным исходным данным был использован метод первого порядка.
Библиографическое описание: ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ВРЕМЕНИ ЗАПАЗДЫВАНИЯ ФИЛЬТРАЦИИ ЖИДКОСТИ В НЕОДНОРОДНЫХ СРЕДАХ // Universum: технические науки : электрон. научн. журн. Эрназаров М.Ю. [и др.]. 2023. 12(117). URL: https://7universum. com/ru/tech/archive/item/16592
¿к UNIVERSUM:
№11X117}_ТЕХНИЧЕСКИЕ НАУКИ_декабрь. 2023 г.
ABSTRACT
In this work, the inverse problem of determining the delay coefficient of liquid filtration time in inhomogeneous porous media is posed and numerically solved. To solve the problem, first and second order methods were used. Solving the inverse problem reduces to the problem of finding the minimum of the functional, i.e. to the unconstrained optimization problem. Using additional unperturbed initial data, the coefficient is practically restored in three to four iterations in various initial approximations. It is shown that the first order method is better than the second order method. For this reason, a first-order method was used to solve the inverse problem with perturbed input data.
Ключевые слова: фильтрация, неоднородная пористая среда, обратная задача, трещиновато-пористая среда, уравнение просачивания, жидкость, конечные разности, метод движения, численное решение.
Keywords: filtration, inhomogeneous porous medium, inverse problem, fractured-porous medium, seepage equation, fluid, iteration, finite differences, motion method, numerical solution.
Введение. Среди неоднородных пористых сред чаще встречаются трещиновато-пористые и пористые среды. В таких средах неоднородности (трещины, каверны, плохо проницаемые породы) распределены по объему равномерно и часто [1,2]. Такие породы называются средами с двумя различными проницаемостями. При этом окружающая среда рассматривается как сложная система, состоящая из двух взаимосвязанных сред. Первая среда состоит из трещины или пор с высокой проницаемостью, а вторая среда состоит из пористых блоков с низкой проницаемостью. Математическая модель истечения жидкости в трещинно-пористых средах описывается в работах [1-3]. Из размерного анализа получено следующее выражение для перетока жидкости между трещинами и пористыми блоками [3]:
Ч = —(р2 -Рх)> И
где р, р2 - давление в трещинах и блоках соответственно, р - плотность жидкости, и - вязкость жидкости, га - безразмерная величина, описывающая флюидообмен между блоками и трещинами.
В случае, когда жидкость слабосжимаема, вязкость постоянна, обе среды упругие и проницаемость обеих сред постоянна, движение жидкости в трещиновато-пористой среде записывается в виде следующей системы уравнений:
„* dp к d2p а / ч
= ТТ" + Р2 -Pi)'
dt u dx u
(1)
„* dp2 k2 д2p2 а / ч
P2^ = + -(P2 - Pi h
dt u dx u
где р* = рс/ + т1 Р7 , (I = 1,2), к1 - проницаемость,
т{ = т01 + Рс/ (р - р0), (I = 1,2), т1 - поPистость, Р* - упругоемкость среды, р - коэффициент сжимаемости среды, р^. - коэффициент сжимаемости
жидкости, х - координата, I - время, индекс I = 1 -соответствует трещинам, I = 2 - пористом блокам.
Если в системе (1) имеем тх << т2 , Рс1 « Рс2, к2 << к тогда мы приходим к следующей системе упрощенных уравнений:
к, d2p, а / ч _
"¡г + (P2 -Pi)=
dx
dp
(2)
dp2 + а(P2 - Pi)= dt u
Вычитая второе из первого уравнения системы (2), получаем следующее уравнение:
dP2 dt
ki d 2Pi up* dx2
(3)
Если взять производную по времени I из первого уравнения системы (2), то уравнение принимает вид:
к d3 P а (dp dP
+
up* dx2dt up*
V
dt dt
= (4)
у
Подставив (3) в (4), получим следующее уравнение относительно p :
dPi = ki d 3 Pi , ki d 2 Pi dt а dx2 dt uP * dx2
к
Введя замены _L = ^ а
Pi = P
в уравнение (5), получим следующее:
(5)
о* ^i
up2
dp _ d3 p d ~ Л dx2dt +Xi
d 2 n
—rr, ° < x < L, ° < t < T, (6)
dx.2
где ц - коэффициент время запаздывания,
Ь - длина одномерной области, Т - максимальное время, в течение которого исследуется процесс.
Эффективные численные методы решения обратных задач математической физики приведены в [4-9]. Некоторые обратные задачи фильтрации жидкости в пористых и трещиновато-пористых средах
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
решены в [10-19], а обратные задачи массообмена в пористых средах решались в работах [20-23].
В данной работе на основе (6) рассматривается обратная задача по определению коэффициента времени запаздывания ц. Сначала дается постановка задачи и затем её решение. Для минимизации функционала невязки используется итерационная процедура по неизвестному коэффициенту. В качестве критерия остановки итерационной процедуры используются как минимизация значения функционала, так и сходимость итерационных значений коэффициента к предельному. Исходные данные для решения обратной задачи подготовлены из решения соответствующей прямой задачи, т.е. проводится "квазиреальный" вычислительный эксперимент. Получены решения обратной задачи и для искусственно возмущенных исходных данных.
1 Постановка обратной задачи
Для решения уравнения (6) введем начальные и граничные условия в следующем виде:
p(0, х) = p0 = const, 0 < х < L,
(7)
Q = Q
M
dp ] d1 p
dx X\ dxdt p (t, L) = p0, Q = const, 0 < t < T.
(8)
Кроме этих условий для решения обратной задачи необходимо дополнительное условие, т.е. известно изменение давления на нефтяной скважине (в точке х=0):
p(t ,0) = z(t).
(9)
Обратная задача ставится следующим образом: коэффициент ц уравнения (6) определяется из минимизации следующего функционала невязки:
т
—(ц)=1\р^ ,о)-Ж )]2 л. (10)
2 Решение обратной задачи методом первого порядка. Условие стационарности (10) функционала будет иметь такой вид:
— = 2\\р{( ,0)- 2 (/ )И' ,0& = ¥ (ц) = 0. (11)
& 0
Здесь ™ = - функция чувствительности.
д]
Разложим в ряд функцию p в окрестности Л с
точностью до членов второго порядка следующим образом [4, 5]
s+1 s f s+1 As
p (t, x)« p(t, x ) + l л-Л I w(t, x). (12)
Поставим разложение (12) в соотношение (11), получаем следующее выражение:
p (t ,0)+f T-ttI w (t ,0)-z(t)
w( t ,0)dt = 0.
Если функции p(t,0) и w(t,0) известны, тогда из последнего выражения легко вычислить прибли-
s+1
жение л [4]:
T =
Ts s s
J T1 w(t,0)- p(t,0)+ z(t) w(t,0)dt
(13)
J Wt,0)
dt
Дифференцируем задачи (6)-(8) по параметру Л получим следующую задачу:
dw д3 w д3 p д2 w
dt T дх2 dt + dx2 dt + % dx2
(14)
0 < x < L, 0 < t < T, w(0, x) = 0, 0 < x < L,
M
/dw T d1 w 1 d1 p^ — + —-+--—
dx % dxdt % dxdt
(15)
0, w(t,L) = 0. (16)
Численный алгоритм определения коэффициента ц можно построить следующим образом:
1. Приведем начальное приближение ц (^ = 0);
2. Решаем задачу (6)-(8) от t = 0 до t = т и определяем функцию р. Решим также задачу (14)-(16)
от t = 0 до t = т и определяем функцию W;
5 +1
3. Вычисляем приближению ц по формуле (13);
4. Этапы 2 и 3 продолжаются до тех пор, пока не будут выполнены следующие условия:
s+1 s
J - J
<e,,
s+1 s
t-T
<F.
— л
где 81 и С 2 - достаточно малые величины.
3. Решение обратной задачи методом второго порядка. Взяв производную функции ¥(ц) в формуле (11) по ц, получим следующее соотношение:
F\ T 1 + dF[]]-TT ! = 0,
(17)
T
0
с,
0
k
n
с
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
то есть
Л=Л- FI Л
dF_ dл
(18)
где
dF
^ 0
j|[ p (t ,0)-z (t )]с( t ,0) + [ w (t ,0)]2 J dt = 0,
C =
dw _ д2 p дл дл
2
(19)
Взяв производную от (14)-(16) по ц, получим следующее задачу:
дс _ д3с дC д3w
~dt ~Л~дХх2дг + X HÖF + ^Föt
0 < x < L, 0 < t < T, w(0,x)= 0, 0<x<L,
(20) (21)
дс л д2с „ д3w — + ---- + 2-
дxдt
дх Xi дх x
= 0,
(22)
ю( I, Ь ) = 0.
Численный алгоритм определения коэффициента ц методом второго порядка состоит из следующих шагов:
1. Приведем начальное приближение ц (5 = 0);
2. Решаем задачи (6)-(8), (14)-(16) и (20)-(22) от| = 0 до i = т и определяем функции р, ж, ю.
Вычисляется функционал (10) и интегралы (11), (19);
5 +1
3. Вычисляем приближению ц по формуле (18);
4. Этапы 2 и 3 продолжаются до тех пор, пока не будут выполнены следующие условия:
5 + 1 5
J - J
<8„
5 + 1 5
Л - Л
<8,,
Q,
j(х,tj), Xi = h i = 0,1,...,N, h
где е и & - достаточно малые величины.
4. Разностная задача. Задачи (6)-(8), (14)-(16) и (20)-(22) численно решаются методом разностных конечностей [24]. Для этого введем сетку в области Б = {0 < х < Ь,0 < I < Т}:
1 Л
М
—, tj = jr, j = 0,1,...,M, т = — }.
Аппроксимируя задачу (6)-(8) с помощью чисто неявной схемы на сетке ПА , получается следующую
разностную задачу:
pj+1 - p', = л 1 r r
j - 2 pj+1 + j p- - 2 pj + pj
i+1
+
pj+1 - 2pj+1 + p
j+1
+ X
j+1
i = 1,2,..., N -1, j = 0,1,...,M -1,
(23)
pi = p(), i = 0,1,...,N,
k
-1
5+1 5
5
1
h
h
4-1
h
а =-4
Ц
p*1
J+1
+
V1
./+1
pj
p 2
h
j = 0,1,..., M -1.
, pN+1 = p
h
h
\
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
Аналогично, аппроксимируя задачи (14)-(16), приходим к этой разностной задаче:
wJ+1 - wJ
1
+ -X
5 1
= Л"
X
j+i
W+1 - 2wJ+1 + W+1 wJ, - 2wJ + wJ,
г-1 г г + 1 г-1 г г + 1
+1 +1
h
h
+
pj+11 - 2 p/+1 + pJ+11 pJ-1 - 2 p/ + p/
+1
h
h
+
+ X1
wJ+1 - 2 wJ+1 + WJ+1
h
, г = 1,2,..., n -1, j = 0,1,...,m -1,
w° = о, г = 0,1,...,n,
+1
w0+1 - w1+1 л
%1X
w0+1 - wJ+1
+
+
%1X
+1
p0 - p1 p0 - p1
1L v
h h y
= 0, wf = 0, j = 0,1,...,M -1.
После аппроксимации задачи (19)-(22) принимают
вид:
ю/+1 -Ю/
2
+ -
X
ю
£} - 2ю/+1 +«J+11 fflJ-1 - 2ю/ +«j+1
= Л
wZ-11 - 2w/+1 + w/++1 w^1 - 2wj + w+
+
+
+ X1
J - 2ю/+1 + юг+11
, г = 1,2,...,N -1, j = 0,1,...,M -1,
ю0 = 0, г = 0,1,...,n,
Ю0+1 -ю1+1 , _Л_ h X1X
Ю0+1 - ю/+1 ю0 - ю2
+
X1X
+1 - w
j +1
J ^J
w0 - w
= 0
, ю J+1 = 0, j = 0,1,...,M -1
Разностные задачи (23), (24) и (25) можно решать методом прогонки. В этом случае запишем разностную задачу (23) в следующем виде:
Лр^ - Ср\+1 + Бр^ =- , / = 1,2,...,#-1, у = 0,1,...,М-1, р0 = Р0, / = 0,1,..., #,
где
A = B = -т-h2
(24)
(25)
Л+ХХ
с = 2A +1, F = P/(p-1 -2pj + p+i)
а=- ^
pn +1 - pj Л
h
p+1 - pj pJ - pj
1 X v
h
pj;1 = p, j = 0,1,...,m -1
X
k
w J - w J
0
h
V
1
h
h
1
5
2
2
h
h
X
X
2
2
h
h
2
h
2
h
h
h
h
h
j
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
Алгоритм решения разностной задачи (26) методом прогонки принимает вид:
а-1, + л
5 (pj - Г! ),
а
i+1
Хг + Л| Х1г + Л
B , ß„ = ^, i = 0,1.....N -1,
C-a A C-а A
ii
pr =aMpÜ +ßi+1, i = N -1,...,1,0, pN+ = p0.
Разностная задача (24) представляется следующим образом:
Aw^1 - Cwj+1 + Bwj+1 = -W/, i = 0,1,..., N -1, j = 0,1,..., M -1
w0 = 0, i = 0,1,...,N,
wQ+1 - wj+1 | h Xr
/ wj +1 - wj+1 wj - wj
+
+
'pJ0 +1 - p1 +1 p0 - P1
Xr
h
= 0, wj+1 = 0, j = 0,1,...,M -1,
где
W = wj (<1 -2wj + wj+1)+ ^[(p-11 -2pr + pi+11 )-(p/_1 -2pj + pj+1)].
Решения разностная задача (27) методом прогонки принимает следующий вид:
(27)
5
а« = 1, ß« (w0 -w1)--((poj+1 -рГ)-(pj -p1)),
Xr + Л
XГ + Л
h
h
1
h
V
а (1) = а i+1
B
C -а(1)A
ß(+ =
Aß (1) + Wj C - а(1) A
, i = 1,2,..., N -1,
w/+1 = а М +ß (+>, i = N -1.....1,0, wN+1 = 0.
Разностную задачу (25) запишем следующим образом:
Ac j+ - Cc j+1 + Bc j+ =-V/, i = 0,1,..., N -1, j = 0,1,..., M -1 c0 = 0, i = 0,1,...,N,
c0+1 -c1+1 Л
4 +1 -cj+1 c0 -c2
h
Xr
1 r V
h
+
+
V+1
Xr
w!+1 wj - w!
= 0, с-г = 0, j = 0,1,...,m -1
k
h
2
h
h
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
где
V/ = ой j - -Л (<1 - 2oj + <1) + A [(wj+11 - 2w/+1 + wj+11) - (w-1 - 2 wj + wj+1)] = 0.
Алгоритм решения разностной задачи (28) имеет
вид:
а
(1)
= 1, ß(1) =
а (1) = а /+1
1г к )--(W+1 - W/+1 )-(w0 - w()),
%1 X + T 5
C
(1)
ß (+i =
%1 X + T
Aß (1) + V/
ш/+1 = а (+1 Ш/+1 +ß (+1
i = 1,1,...,N -1.
C -а(1) A i = N -1,...,1,0, ш/+1 = 0
5 Результаты численного расчета. Для проведения квазиреального эксперимента [6] сначала рассмотрим прямая задача (6)-(8) с заданным зна-
__exact г\ 1
чением л = 0,1 решается численно методом конечных разностей [24]. По результатам численных расчетов определяется сеточная функция zj = z (t ),
j = 0,1,...,M. Функция z(t ) служит исходными
данными для решения обратной задачи. Для моделирования погрешностей функция z(t) зашумляется случайными погрешностями следующим образом:
zj = zj +15
ö
a / --
(29)
где ау - случайная функция, равномерно распределенная на интервале (0,1).
Для численных расчетов использованы следующие значения исходных данных: Т = 2000 с, Ь = 10 м, ^ = 1. 10-13 м2, ро = 10 МПа, ц = 1. 10-8 МПас, ^ = 1 • 10-3 м2/с, <2о = 1. 10-5 м/с. Разностная сетка делит координатный участок \0, Ь] на 200 интервалов, а временной участок \0, т] на 4000 интервалов.
был подготовлен в 200 точках «времени» на основе сеточных решений «данных измерений» (29) (Рис. 1).
2
V
10
го 9
о_
10
500 1000 1500 2000 t, s
11 10 9 8 7
7 Г
— 5=0.0 ■" 5=0.02
500 1000 1500 2000 t, s
7
V
------5=0.0
......... 5=0.01
0 500 1000 1500 2000 t, s
0 500 1000 1500 2000 t, s
Рисунок 1. График сеточной функции z* (t).
q
Я
Я
7
Г)
Г)
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
На рис. 2 показаны численные расчеты (чего?) методом идентификации первого порядка при 5 = 0 . Как видно из рис.2 (а-ф, различные начальные приближения ц восстанавливаются за две-три
0
итерации. Начальные приближения ц дает хорошие результаты, когда в 10 раз меньше и 6,5 раз
больше чем ц^^ (Рис.2, a-d). Коэффициент ц
восстанавливаются за две-три итерации. При
0
ц = 0,687 процесс итерации становится несходящимся (Рис.2, d). Численные расчёты показывают,
0
что когда начальные приближение ц изменяются
0
в диапазоне 0,0001 < ц < 0,687 метод первого порядка дает хорошие результаты.
1.2 1
0.8
® 0.6
^ 0.4 0.2 0
0
4
3 I
2 1 0
3 2.5
1
1
0.5
0
12 s
Л0=0.15
Л0=0.2
-ЕЪ- л0=0.41
Л0=0.53
Л0=0.65
Рисунок 2. Восстановление коэффициента ц по невозмущенными данными (при 5 = 0)
методом первого порядка
На рис.3. показаны численные расчеты методом идентификации второго порядка при 5 = 0. Как видно из рис.3 Ь), различные начальные приближения ц восстанавливаются за три-четыри итерации. Численные расчеты (рис.3, a-b) показывают, что когда начальное приближение ц0 изменяется в диа-
пазоне 0,04 < ц < 0,34 коэффициент ц восстанавливается за три-четыре итерации. При начальных 0 0 приближениях ц = 0,001 (рис.3, а) и ц = 0,34
(рис.3, Ь) процесс итерации становится несходящимся, то есть не удаётся восстановит коэффициент ц .
b
a
1
2
3
1
2
3
5
0
3
n
1
5
с
s
№ 12 (117)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
декабрь, 2023 г.
1.5
0.5
-0.5
3.5
л0=0.001 л0=0.006
^°=0.02
„0_
л0=0.04
0
23 s
2.5
1.5
со
1
0.5
-0.5
0
23 s
Рисунок 3. Восстановление коэффициента Ц заданными невозмущенными данными (при 5 = 0)
методом второго порядка
Численные расчеты рис. 2 и рис. 3 показывают, что метод первого порядка дает наилучшые результаты, по сравнению с методом второго порядка. Такие покозатели могут быть результатом ошибки вычисления, поскольку в методе первого порядка в
каждом итерации вычесляется две прямые задачи, а в методе второго порядка в каждом итерации вычисляется три прямые задачи.
Результаты расчетов с возмущенными исходными данными приведены в табл. 1.
Таблица 1.
Восстановления коэффициента Ц при заданными возмущенными данными
b
a
1
4
1
4
S = 0,01 м2 = 0,4 м2
5 , м2 Относительная погрешность, 1», „exact I 5 , м2 Относительная погрешность, 1», „exact I J^-100, %
0,0 3 0,09999980 0,000200 3 0,10000000 0,000000
0,005 4 0,09957584 0,424160 4 0,10033725 0,337250
0,01 5 0,09630044 3,699559 5 0,10455201 4,552010
0,02 7 0,11177254 11,772542 6 0,10578710 5,787100
0,05 8 0,08595013 14,049868 8 0,13199806 31,998060
Из таблицы 1 видно, что относительные погрешности восстановления коэффициента Ц изменяются
от 0,3373% до 31,9981%. Из табл. 1 видно, что при малых погрешностях восстановления коэффициента Ц дает наилучшие результаты.
Вывод. Поставлена и решена коэффициентная обратная задача фильтрации жидкости в неоднородных пористых средах. Определено коэффициенты времени запаздывания фильтрации жидкости. Для решения задачи использованы два метода - методы
первого и второго порядка. Численные результаты, полученные на невозмущенных данных, показывают, что метод первого порядка точнее, чем метод второго порядка. Исходя из этого, для расчетов восстановления коэффициента Ц при возмущенных исходных
данных применены метод первого порядка. Численные результаты расчеты показывают, что в случае с возмущенными исходными данными коэффициент времени запаздывания восстанавливается с ошибками, соответствующими порядкам возмущения.
¿к UNiVERSUM:
№ 12 П171_ НАУК1_декабрь. 2023 г.
Список литературы:
1. Баренблатт Г.И., Желтов Ю.П. Об основных уравнениях фильтрации однородных жидкостей в трещиноватых породах // ДАН СССР. 1960. Т.132, №3. С. 545-548.
2. Баренблатт Г.И., Желтов Ю.П., Кочина И.Н. Об основных представлениях теории фильтрации однородных жидкостей в трещиноватых породах // ПММ. 1960. Т.24, вып. 5. С. 852-864.
3. Баренблатт Г.И., Ентов В.М., Рыжик В.М. Движение жидкостей и газов в природных пластах. Москва.: «Недра», 1984.
4. Бабе Г.Д., Бондарев Э.А., Воеводин А.Ф., Каниболотский М.А. Идентификация моделей гидравлики. Новосибирск: Наука, 1980.
5. Алифанов О.М., Артюхин Е.А., Румянцев С.В. Экстремальные методы решения некорректных задач. Москва: Наука, 1988.
6. Samarskii A.A., Vabishchevich P.N. Numerical Methods for Solving Inverse Problems of Mathematical Physics. Berlin: Walter de Gruyter, 2007. - 438 p.
7. Алифанов О.М. Обратные задачи теплообмена. Москва: Машиностроение, 1988.
8. Hao D. N. Methods for inverse heat conduction problems. - Lang. pub. Inc., Peter, 1998.
9. Beck J.V., Blackwell B., Clair C.R. Inverse Heat Conduction: Ill-Posed Problems. Wiley, 1985.
10. Хайруллин М.Х., Хисамов Р.С., Шамсиев М.Н., Фархуллин Р.Г. Интерпретация результатов гидродинамических исследований скважин методами регуляризации. Регулярная и хаотическая динамика, Ижевский институт компьютерных исследований, 2006.
11. 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. No. 11, pp. 35-40.
12. Khairullin M. H, et.al. Thermohydrodynamic studies of vertical wells with hydraulic fracturing of a reservoir. High Temperature, 2011. Vol. 49. No. 5, pp. 769-772.
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., No.3, pp. 517-525.
14. Khuzhayorov B.K., Kholiyarov E.C., Mamatov S. Identification of the coefficients of relaxation filtration of a homogeneous liquid in porous media. AIP Conference Proceedings. - AIP Publishing LLC, 2022. - Vol. 2637. No. 1. 040009.
15. Хужаёров Б.Х., Холияров Э.Ч. Определение коэффициента перетока и проницаемости при фильтрации однородной жидкости в трещиновато-пористых средах // Проблемы вычислительной и прикладной математики.
2022. №1(38). С. 66-76.
16. Хужаёров Б.Х., Холияров Э.Ч., Эрназаров М.Ю., Тураев М. Обратная задача по определению коэффициента перетока в модели фильтрации Уоррена-Рута // Научный вестник СамГУ. 2022. №1(131). С. 115-123.
17. Холияров Э.Ч., Шадмонов И.Э., Шерпулатов Ш.Ш. Определение коэффициента перетока в модели фильтрации Уоррена-Рута на основе решения обратной задаче // Universum: технические науки : электрон. научн. журн. 2022. №5(98). C. 51-58.
18. Kholiyarov E.Ch., Ernazarov M.Y., Jurayev O.A., et al. Coefficient inverse problem for a simplified model of filtration of a homogeneous fluid in fractured-porous medium. - AIP Conference Proceedings 2637, 040021 (2022)
19. Kholiyarov E.Ch., Ernazarov M.Y. Determination of Relaxation and Flow Coefficients During Filtration of a Homogeneous Liquid in Fractured-Porous Media // Current Problems in Applied Mathematics and Computer Science and Systems. - Cham : Springer Nature Switzerland, 2023. .- P. 54-67.
20. 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.
21. Khuzhayorov B., Begmatov T., Kholiyarov E.Ch., Fayziev B. The inverse problem of determining the kinetic coefficients in the model of suspension filtration in a porous medium // Eurasian Journal of Mathematical and Computer Applications. 2022. Volume 10, Issue 4 (2022) P. 96-106.
22. Холияров Э.Ч., Бегматов Т.И., Файзиев Б.М. Идентификация коэффициента кинетики в модели фильтрации суспензии в пористой среде // Проблемы вычислительной и прикладной математики. 2022. №1(38). С. 9-17.
23. Khuzhayorov B., Kholiyarov E.Ch., Khaydarov O.Sh. Inverse Problem of Contaminant Transport in Porous Media // Current Problems in Applied Mathematics and Computer Science and Systems. - Cham : Springer Nature Switzerland,
2023. - P. 87-97.
24. Самарский А.А. Теория разносттных схем. - М.: Наука. 1989. - 616 с.
