Идентификация параметров газоносного пласта на основе решения обратной задачи теории фильтрации

УДК 517.958:532 ББК 22.181 Ш 96

М.М. Шумафов, Р. Цей

Идентификация параметров газоносного пласта на основе решения обратной задачи теории фильтрации

Аннотация

В работе определяются фильтрационные и емкостные параметры газоносного пласта на основе решения обратной задачи теории фильтрации методом модулирующих функций.

Ключевые слова: обратная задача, математическое моделирование, идентификация, модулирующие функции, фильтрационные и ёмкостные параметры.

M.M. Shumafov, R. Tsey

Determination of gas reservoir parameters on the basis of a solution of an inverse problem of the filtration theory

Abstract

In this paper, the filtrational and capacity parameters of gas reservoir are determined. The method of solving this problem is based on the solution of inverse problem of filtration theory. In this case the modulating functions method is used.

Key words: inverse problem, mathematical simulation, identification, modulating functions, filtrational and capacity parameters.

Введение

О той большой роли, которую играют обратные задачи при математическом моделировании в различных областях науки, отмечено в работах [1,2]. Напомним, что к прямым задачам математической физики относят задачи нахождения следствий заданных причин, а к обратным - задачи отыскания неизвестных причин заданных следствий. Обратные задачи возникают в приложениях и имеют важное значение при решении вопросов математического моделирования в сложных системах. Это, в частности, относится и к процессам, связанным с добычей нефти и газа - одним из самых приоритетных направлений в энергетике.

В настоящей статье нами рассматривается обратная задача теории фильтрации, а именно, задача определения обобщенных фильтрационно-емкостных параметров газоносного пласта. Для решения этой задачи применяется метод модулирующих функций (далее, М-метод). Отметим, что идея применения М-метода для решения обратных задач восходит к работам Дж. Лоэба и Г. Кахена (J. Loeb, G. Cahen) [3, 4]. Возможность применения М-метода для решения задач нефтегазовой науки впервые была высказана В.Б. Георгиевским и им были разработаны унифицированные алгоритмы для решения обратных задач подземной гидрогазодинамики [1]. В работах [5,6] приведены некоторые способы программной реализации алгоритмов, предложенных в работах В.Б. Георгиевского (см., например, [1]). В работе [7] сделана попытка обобщить М-метод на случай любой степени полиномов разложения неизвестных параметров газоносного пласта.

Математическая постановка задачи

Рассмотрим нелинейное эволюционное уравнение параболического типа, описывающее процесс нестационарной фильтрации реального газа [8]

_3_ 3 х

k(х, у^(х, у) 3р2

т (р) г (р) 3 х

3 у

k(х, y)h(х, у) 3р2

3

= 2а (х, у)т(х, у^(х, у)—

3 ^

т (р)^(р) 3у

+ ^(х, у, ^)ра

,2( Р\

(х, у) е О , t е (0,Т), Ос я2 с начально-граничными условиями

Р(х у,0) = Ро(х у), (х у) е 0 ,

Lpг = ю

Здесь р( х, у, t) - давление в точке пласта с координатами (х,у) в момент времени t, k(х,у) - коэффициент проницаемости пласта, т(х,у) - коэффициент пористости пласта, ^х,у) - эффективная толщина пласта, ц(р) и г(р) - соответственно коэффициенты динамической вязкости и сверхсжимаемости газа при давлении р и пластовой температуре Гпд, а(х,у) - коэффициент газонасыщенности, Q(х,y,t) - объемный расход газа, отнесенный к единице площади пласта в точке (х,у) в момент времени t, приведенный к атмосферному давлению рат и пластовой температуре Тпл. Далее, L - это единичный оператор в случае первой краевой задачи; в случае второй краевой задачи L =3 /3 п - производная по внешней нормали к границе Г области О ; в случае третьей краевой задачи Lp = 3р / 3 п + Q(х, у, 0р (Q(х, у, 0, ю (х, у, 0 - известные функции).

Отметим, что оператор L в условии (3) может иметь и более сложную структуру.

Как известно, если выполнены условия согласованности начальных и краевых условий при заданных k, h, т , г,а , т, Q, рат, р0, ю , то задача нахождения функции р

, удовлетворяющей условиям (1)-(3) в классическом или обобщенном смысле, называется прямой задачей. При решении прямой задачи возникает ряд вопросов, в частности, вопросы однозначной разрешимости задачи (1)-(3); корректности в том или ином смысле и т.д. (см., например, [9, 10]).

Предположим теперь, что решение задачи (1)-(3) - р( х, у, t) - известно и требуется найти функции а( х, у) := k (х, у^( х, у) и Ь( х, у) := 2а т( х, y)h( х, у). Функции а( х, у) и Ь( х, у) будем называть обобщенными фильтрационным и емкостным параметрами.

Сформулированная выше задача явялется обратной задачей. Возникает естественный вопрос: существует ли решение обратной задачи, если существует, то единственно ли оно?

Рассмотрим сначала случай, когда а( х, у) и Ь( х у) являются однородными многочленами степени п р:

а( х у) = X X

а

кт

хк - тут

к = 0 т= 0

Чх, у) = X X Ьт,хк- тут

к = 0 т= 0

Требуется найти коэффициенты акт , ькт в

являющаяся решением уравнения

3 а( х, у) 3 р2 3 + у,) X 3 р 2 = Ь( х, у) -3- / \ р

3 х т(р)г(р) 3 х _ 3 у . т(р) г(р) 3 у . 31 г (р) V 1 /

+ с

3

в предположении, что т (р), г(р), с(х, у, t) - известные функции.

Для решения сформулированной задачи нами применяется метод модулирующих функций. Как было ранее отмечено, при этом не требуется знать конкретный вид начально-граничных условий (2), (3).

Введем в рассмотрение функцию

У

і т

где р = р( X ^t).

Частная производная по х от У (р(х, у, t)) равна

3у (р) ёу (р) 3р_ 2р 3р_

1

3 р2

Аналогично

3 х ёр з х т (р)2 (р) з х т (р)2 (р) з х

3 у (р) _ 1 3р2

3 у т (р)2 (р) 3 у

Учитывая соотношения (7)-(9), перепишем уравнение (6) следующим образом:

_3_ 3 х

а(x, у) ^У (р) 3х

_3_ 3 у

а

3 У

= Ъ( x, У) 3^7

3 7

2

+ с( х, у, 7).

Вместо уравнения (10) рассмотрим несколько более общее уравнение:

_3_ 3 х

а

3У і

3

+ -----

3 у

а

3У

3 у

Ъ( х у) Ч7Г +у

3 7

3х

где у 1 = У I(х,у,^р), у 2 = у 2(х,у,t;р) , у 3 = у 3(х,y,t;р), у 4 = у 4(х,у,^р) - известные функции, функция р = р( х, у, t) задана, а а(х, у) и Ь(х, у) имеют вид (4), (5).

Требуется найти акт.

Ъ

кш

Теоретическое обоснование применения метода модулирующих функций

Пусть известна функция р = р(х, у, t) ((х, у) е G , t е [0,Т]), обладающая требуемой гладкостью и удовлетворяющая уравнению (11), где а( х-, у) и Ь( х, у) имеют вид (4), (5). Требуется найти а( х, у), Ь( х, у).

Сначала задачу (11) решим в случае, когда G есть прямоугольник: {(х, у) е Я2: х е [х1, х2], у е [у1, у2]}. Умножим уравнение (11) на функцию ^(х,у,0 = /(х)/(у)/(t), где /(х)е С2[х^,х2], /(у)е С2[у 1,у2], /(t)е С'[0,Т] и

/ (х1) = / (х2) = 0, / (х1) = /'(х2) = 0, I (у1) = / (у 2 ) = 0, I' (у1) = I' (у2) = 0,

I (0) = I (Т) = 0.

+

2

Затем полученное соотношение проинтегрируем по параллелепипеду V = {[х1,х2]х [у1,у2]х [0,Т]} и произведем интегрирование по частям с учетом условий (12)-(16). Для первого слагаемого в левой части уравнения (11) имеем

_3_ 3 х

а

3У 1

3х

Т у 2 х2

Г (х, у, t)dхdydt = [ [[ —

3х

0 у1 х1

а

3У 1

у2 х2 3

/ ()dt( / (у)ф\

а

3У 1

3х

3х / (х)Лх.

0 у1

Для внутреннего интеграла правой части (17) с учетом условий (12), (13) получаем

_3_ 3 х

а( х, у)3у 1

3х

/ (х)Лх = а( х, у) -У-1 У (х) 3х

( а( х, у) 3-У-1 /' (х)Лх =

3х

х2

( а(х, у) у-1 /(х)Лх =- a(х, у^' (х)у 1 х + ( у 1[а( х у)/' (х)1А =

У \ [а(х,у)/'(х)1хЛх.

Из (17) и (18) следует

V

_3_ 3 х

а

3У 1

_3_ V 3 у

а

3х

ос

3 у

Т у2 х2

Г(х, y, t= (/^)Л(/(у)Лу( У1[а(х, у^'(х%&.

0 у1 х1

Аналогично, для остальных слагаемых (11) имеем:

Т х2 у 2

Г (х, y, ^ = (У (0^ (У (х)Лх (у 2[а( х у)/'( х)]Улу ,

0 х1 у1

3у

Ь(х, у) (х, у, t)dV = - ( /(х)Лх( /(у)Ь(х, у)Лу( /'(t)у ъdt,

у2

Из (11), (19)-(22) получаем следующее интегральное тождество

у2

/^)dt \ /(у)Лу I у 1[а(х, у)/'(х)]!А + \ /^)dt IУ(х)Лх I у 2[а(х у)/'(х)]УЛУ+

у2

у2

Т

+ ( /(х)Лх( /(у)Ь(х,у)Лу( /'^)у зdt = ( у 4у(х)/(у)/-^)dV.

х1 У1 0 V

Учитывая вид (4), (5) функций а( х, у) и Ь( х, у), тождество (23), после некоторых элементарных преобразований, примет вид:

р к

1 1 (а ктакт + Ь ктЬкт) = Г ,

к = 0 т= 0

где

кт

(/(){У ,/(У)[хк-тут/"(х) + (к - т)хк-т'ут/(х)] +

V

+ у /(х)[хк-тут/• (у) + тх‘-тут-'/(у)]Л',

Ькт = (У 3/(х)/(У)/'(')х‘-"У"ЛУ, 7 = (у 4/(х)/(у)/(г)<^,

V V

(V = [ х1, х2] х [yl, у2] х [0,Т]).

В соотношении (24) акт , Ькт - неизвестные, а а кт, Ь кт, 7 - известные числа.

хх

х

х

х

Т

Т

х

х

Нетрудно установить, что при k - пр, т - 0,1,..., k число неизвестных {акт}

и

/а } (пр + 1)(пр + 2) п - п б

{ькт} соответственно равно па - —--------------- и пь = па, а общее число неизвестных

равно пе - (пр + 1)(пр + 2).

Для однозначного определения пе неизвестных необходимо иметь пе линейно независимых алгебраических уравнений, причем детерминант матрицы этой системы должен быть отличным от нуля. С этой целью умножим уравнение (11) на функции

Р, (х У, *) - I,(х)/, (У)/, (*) , где /, (х) /, (У), /, (*) удовлетворяют условиям (12)-(16) при ] - 0,1,...,пе - 1, и проинтегрируем полученное уравнение по параллелепипеду V. Имеем:

а т = ¡ /,(0{У / №‘-тут/,(х) + (к - т)хк'тУЦ(х)]+

'

+ У г/,(х)[хк-тут/'(у) + тхк-тут-'/'(у)])^,

Ь,' - (у 31,(X), 1,(у)Xі--у-]](I)Л’,

V

11) - 1 у 4 ]. (X)], (у)], (I)Л'

Г

(і - 0,1,...,

п.

(27) т - 0,1,...,і

где V - {[x1, Х2 ] X [^1а У2] X ^ ґ2]},

] - аи, п -1).

Введем в рассмотрение матрицу W размера пе X пе

W - }, :- {а ь кт)},

и пе - мерный вектор

0 - {акш, Ьт }Г .

Вместо двойных индексов кт мы иногда будем использовать одинарный индекс і. Таким образом, индекс і будет обозначать номер столбца, а ] - номер строки матрицы W. Тогда соотношение (24) можно переписать в векторно-матричной форме

W • 0 - У .

Нетрудно показать, что, выбирая специальным образом систему модулирующих функций ],(х) ],(у) ],(|) , можно добиться выполнения условия невырожденности матрицы W: detW Ф 0.

Теперь предположим, что функции а( х, у), Ь( х, у) являются достаточно гладкими, а именно

а

Представим а( х, у), Ь( х, у) по формуле Маклорена:

п 1

а(х у) - а(0,0) + X 77 к-1 к!

п 1

Ь(x, у) - Ь(0,0) + X к!

к- 1 к!

3 3

х—+ у— 3 х 3 у

к

а(0,0) + е 1,

33

х—+ у — 3 х 3 у

к

Ь(0,0)+ е 2,

где

(п + 1)!

33

х—+ у— 3 х 3 у

\ п+ 1

а(вх,ву),

1

е

(п + 1)!

З З

х—+ у— З х З у

\ п+ 1

Ь(8х,8у),

(0 < 8 < 1). Выражение

ЗЗ

х—+ у — З х З у

\ k

а(0,0) есть однородный многочлен по х и у степени k.

Поэтому соотношения (29) и (30) можно переписать в виде

а(ху) = І І а^хк~тут + еі,

k = 0 т= 0

Ь( х, у) =І І ЬктХк- тут + Є 2.

к = 0 т= 0

В силу условий (28) на основании (31) и (32) для £ 1 и £ 2 справедливы оценки

т

(п + 1)! ’

(п + 1)!'

Здесь т 1 (соответственно т 2) константа, зависящая от чисел

/1 = х2 - х1, 12 = у2 - у1, тах

З

к, + кл

(х> у"> ° З хк1 З ук'

З к1 + к:

(соответственно тах

а( х, у),

Ь( х, У) ),

(х-у>- о 3 хkl 3 yk2 k1 + k2 = и + 1, ^,k2 = 0,1,...,пр.

Следовательно, при выполнении условий гладкости (28) мы можем определить коэффициенты а( х, у) и Ь( х, у) с точностью, зависящей от степени их гладкости.

Если предположить, что а( х, у), Ь( х, у) являются аналитическими функциями своих переменных, то можно их восстановить по давлению р = р( х, у, t) с наперед заданной точностью.

Решение обратной задачи теории фильтрации с использованием метода модулирующих функций

Вернемся к нашей основной задаче. Уравнение

_З_ З х

З

а

Зх

З

+ -----

З у

З

а

Зу

Ь( х, у) Зі

З і

*

+ с( х, у, і),

где У (р):

2Ї4

о т

Как следует из вышеполученных результатов, в предположении, что функции а( х, у), Ь( х, у) являются достаточно гладкими, аналогично (24), мы получим следующее соотношение:

пР к

І І (а к"акт + Ь2Х) = 1<Л, j = °Д.....п. - 1

к = 0 т= 0

е

2

р

где значения а к—", Ь ш , 1(1) равны соответственно (25), (26), (27), а функции а(х, у) и Ь( х, у) представлены по формуле Тейлора (см. (33), (34)).

Проводя рассуждения, вполне аналогичные вышеприведенным, для определения неизвестных коэффициентов акт, Ькт получим систему алгебраических уравнений (35), где

а к— = (/,(гУ (Р(х,у,г)){/,(у)[х‘-"у/х) + (к- —)хк-—1 у-/’(х)]+

К

+ /, (х)[хк-"у"/'(у) + тхк - -у-- /'(у)]^,

ЬГ = (у (Р(х,у,г))хк-ту"/;(х)/;(у)/;(гук, 1';) = (с(х,у,г)/,(х)/,(у)/,(г)dF,

К V

(V = {[х1,х2]х [у1,у2]х [г1,г2]}, dV = dхdydг, к = 0,1,...,ир ; — = 0,1,...,к ; ; = 0,1,...,пе - 1).

Здесь функции /;(х), /;(y), /;(г), 1 = 0,1,...,пе - 1 являются гладкими (/,(х)- С2, /;(у)- С2, /;(г)- С1 ) линейно-независимыми функциями, удовлетворяющими условиям

/; (х1) = /1 (х2) = 0,

/] (х1) = /; (х2) = 0,

/; (у1) = /1 (у 2 ) = 0,

/; (у1) = / (у2) = 0,

/, (0) = /, (Т) = 0.

Ясно, что вид модулирующих функций определяется геометрией области О, а их минимальная гладкость - порядком дифференциальных операторов, входящих в рассматриваемую математическую модель. Выбор модулирующих функций зависит от специфики каждой конкретной задачи. Кроме того, их выбор влияет на обусловленность матрицы полученной системы алгебраических уравнений. Эти вопросы будут предметом более подробного рассмотрения в другой работе.

Некоторые обобщения математической модели

Выше вкратце был изложен М-метод для решения обратной задачи теории фильтрации на примере дифференциального уравнения, описывающего процесс неустановившейся фильтрации газа в пористой среде.

Случай, когда коэффициенты а(х, у, р) , Ь(х, у, р) имеют специальную структуру

Рассмотрим теперь случай, когда коэффициенты а

_з_

3 х

а

3У 1

3х

_з_ 3 у

а

3 у

3у

7/ ч3у 3

= Ь(х ^ р)~т^ +У 3г

У 3(х,y,р(хy,г)),

параллелепипеде

Здесь р = р(х,y,г), у 1 = у 1 (у,г,р(х,у,г)), у 2 = у 2(х,г,р(х,у,г)), у 3 У 4 = у 4(х, у, г, р(х, у, г)) - заданные функции в

[х1 < х < х2] х [у1 < у < у2] х [0 < г < Т] (на поверхности параллелепипеда не обязательно требуется знать значения этих функций). Требуется восстановить функции а(х, у, р) ,

2

+

4

Ъ( х, у, р) по этим данным с требуемой точностью в зависимости от степени их гладкости. Сначала, мы выделим один случай, когда не требуются разложения по степеням р функций а(х, у, р), Ъ(х, у, р) .

Пусть функции а(х, у, р) и Ъ( х, у, р) представимы в следующем виде: а(х, у, р) = А(х, у) (р), Ъ(х, у, р) = В(х, у)§ (р) .

В таком случае, после подстановок выражений

♦ 1 (г,у,р) = ((X))у ^у,Х) 4 ,

'0

♦ ,

У

♦ 3

3Х

3у 3

3Х

dX

♦ 4(х У, г, р) = у 4(х, ^ г, р), в уравнение (36), получим следующее соотношение

_3_ 3 х

3 ♦

А( х, у)3 1

3х

3

+ -----

3 у

3 ♦

А( х, у)3 2

3у

3 ♦

В( х, у)^ + ♦ v 'У) 3 г

Уравнение (37) принадлежит классу уравнений, рассмотренных выше.

Случай, когда коэффициенты а(х, у, р) , Ъ(х, у, р) являются произвольными гладкими функциями

Предположим, что функции а(х, у, р) , Ъ(х, у, р) являются достаточно гладкими функциями своих аргументов х, у, р (или аналитическими функциями). Представим их по формуле Маклорена по аргументу р. Отбросывая остаточный член, будем иметь:

а

к = 0 к = 0

Подставляя (38) в уравнение (36) после замены {у г} на {♦ г }г= 1,2,3 в силу следующих интегральных соотношений

♦ и- (г, у, р) = ( к dx,

0

♦ (г,х,р) =((Xк 3у 2(X1 dx

4

♦ 4(х, У, г, р) ° у 4(х У, г, р),

получаем следующее дифференциальное уравнение в частных производных эволюционного типа

_3_ 3 х

г 3

1 ак(х, У) 3- ф и

3х

к = 1

3

+ -----

3 У

г 3

1 ак (^ У) 3“ Ф 2,к 3у

к = 1

г 3

1 Ък(х у) 37Ф

3,к + ♦ 4.

к = 1

Уравнение (39) принадлежит рассмотренному выше классу уравнений. Следовательно, мы можем применить М-метод к уравнению (39) для определения неизвестных коэффициентов ак(х,у), Ък(х,у), предварительно представив функции ак (х, у) , Ък (х, у) по формуле Маклорена и оставляя лишь конечную сумму, обеспечивающую нужную точность.

Ясно, что аналогичные рассуждения, проведенные относительно уравнения (36) для определения коэффициентов а(х, у, р) , Ъ(х, у, р) М-методом применимы также и к более широкому классу уравнений, чем (36):

п 3

1 г

Г1 3 х

а1

3у

3 хг

1 Ъг

Г 1

3_ь

3 хг

+ с(хи х2,..., х„ , г, р)^ 1( хl, х2,..., х„ , г, р) =

3 ♦

= d (хl, х2,..., х„ , г, р^^-^ + ♦ 3(х^ х2,..., хп , г, р).

3 гг

Здесь функции у г, ф г предполагаются явно не зависящими от хг и зависящими от остальных независимых переменных {х;-}, от г и р (т.е. у зависит косвенно от хг через р = р(х1,х2,...,хп,г)); функция ♦ 2 явно не зависит от г, т.е. ♦ 2 =♦ 2(х1,х2,...,хп,р) . Предполагается, что р = р(х1, х2,..., хп, г), у , ф , Г = 1,2,..., п, ♦ 1, ♦ 2, ♦ 3 - заданные функции, а аг, Ъг, с и d подлежат определению.

Замечание. М-метод применим и в случае, когда функции у г-, ф явно зависят и от

переменного хг, а ♦ 2 - от г. Но в этом случае интегральные формулы для определения неизвестных коэффициентов разложения неизвестных функций в многомерный ряд Маклорена будут содержать под интегралом частные производные по хг, (Г = 1,2,..., п) и

по г от функции р(х1,х2,...,хп,г), что не очень желательно в случае, когда р - таблично заданная функция. Последнее объясняется тем, что, операция дифференцирования от таблично заданной функции является, как хорошо известно, некорректной операцией. Основная суть М-метода и заключается в том, чтобы избежать операцию дифференцирования от экспериментально заданной функции р.

Заключение

В работе приводится математическая постановка обратной задачи теории фильтрации. Задача сводится к определению коэффициентов нелинейного эволюционного уравнения параболического типа. Для ее решения применяется метод модулирующих функций. В работе приводятся также некоторые обобщения рассматриваемой математической модели. На основе приведенного в работе решения рассматриваемой задачи можно разработать алгоритм для численного решения поставленной задачи.

1. Георгиевский В.Б. Унифицированные алгоритмы для определения фильтрационных параметров. Справочник. Киев, 1971. 328 с.

2. Мирзаджанзаде А.Х., Хасанов М.М., Бахтизин

Р.Н. Моделирование процессов

нефтегазодобычи. Нелинейность,

неравновесность, неопределенность. М.; Ижевск, 2004. 368 с.

3. Loeb J., Cahen G. Extraction, a partik des

enregistrements de mesures, des paramétrés dynamiques d um system // Automatisme. 1963. № 12. P. 17-28.

4. Loeb J., Cahen G. More about process identification

// Trans. on Automatic Control. 1965. P. 359-361.

5. Лукнер Л., Шестаков В.М. Моделирование

геофильтрации. М.,1976. 407 с.

6. Трофимов В.В., Батищева Г.А. Реализация на

ЭВМ унифицированных алгоритмов В.Б. Георгиевского // Сборник научных трудов ЮжНИИгидротехники и мелиорации. 1976. Вып. 9. С. 111-114.

7. Юдин А.И., Юдина О.К. Расчет фильтрационно-

ёмкостных параметров по промысловым данным эксплуатации газового

месторождения // Термодинамика

кооперативных процессов в гетерогенных средах. Тюмень, 1985. С. 80-85.

8. Закиров С.Н., Лапук Б.Б. Проектирование и

разработка газовых месторождений. М., 1974. С.39.

9. Тихонов А.Н., Арсенин В.Я. Методы решения

некорректных задач. М., 1986. 287 с.

10. Hadamard J. Le probleme de Cauchy et les

equations aux derivees partielles lineaires hyperboliques. Paris, 1932.

1. Georgievsky V.B. Unified algorithms for

definition of filtrational parameters. Directory. Kiev, 1971. 328 pp.

2. Mirzadzhanzade A.H., Khasanov M.M., Bakhtizin

R.N. Simulation of petroleum-gas production processes. Nonlinearity, nonequilibrium, uncertainty. Moscow - Izhevsk, 2004. 368 pp.

3. Loeb J., Cahen G. Extraction, a partik des

enregistrements de mesures, des parametres dynamiques d um system // Automatisme. 1963. No. 12. P. 17-28.

4. Loeb J., Cahen G. More about process

identification // Trans. on Automatic Control. 1965. P. 359-361.

5. Lukner L., Shestakov V.M. Modelling of a

geofiltration. M., 1976. 407 pp.

6. Trofimov V.V., Batishcheva G.A. Realization with the help of computer of the V.B. Georgievsky unified algorithms. Proc. Southern Research Hydrotechnical and Melioration Institute, 1976, Issue 9. P. 111-114.

7. Yudin A.I., Yudina O.K. Calculation of filtrational-

capacity parameters from trade data on exploitation of a gas deposit // Thermodynamics of cooperative processes in heterogeneous environments. Tyumen, 1985. P. 80-85.

8. Zakirov S.N., Lapuk V.V. Designing and

development of gas deposits. M., 1974. P. 39.

9. Tikhonov A.N., Arsenin V.Ya. Methods of

solution of incorrect problems. M., 1986. 287 pp.

10. Hadamard J. Le probleme de Cauchy et les

equations aux derivees partielles lineaires hyperboliques. Paris, 1932.