Метод модулирующих функций и его применение при решении обратных задач Текст научной статьи по специальности «Математика»

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

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

Метод модулирующих функций и его применение при решении обратных задач

(Рецензирована)

Аннотация

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

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

Введение

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

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

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

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

обзор которых имеется, например, в [1]. При этом в разных работах использовались различные методы решения обратных задач.

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

Однако унифицированные алгоритмы В. Б. Георгиевского не получили широкого применения для решения прикладных задач, за исключением работ [6-8].

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

Основные недостатки существующих методов определения фильтрационных параметров

Подавляющее большинство существующих методов определения параметров систем (например, метод стохастической аппроксимации и методы теории чувствительности, используемые в регрессионном анализе, методы, основанные на применении преобразований Лапласа, метод детерминированных моментов и др.), описываемых дифференциальными уравнениями, основано на использовании решений уравнений. Разработке методов посвящено весьма большое число работ. Однако всем этим методам присущи серьезные недостатки [1, 5].

К недостаткам существующих методов определения фильтрационных параметров следует отнести следующее:

- низкая точность определяемых параметров,

- трудности обработки опытных данных,

- необходимость организации опытных работ по специальной схеме,

- невозможность интерпретации сложных случаев фильтрационного движения,

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

Главный из этих недостатков - низкая точность определяемых параметров. Причиной тому является абстрагирование от реальных условий и схематизация натурного движения. Тем самым, априори в основу математической модели закладываются параметры, определяемые подчас с заведомой большой погрешностью.

Имеются также и другие недостатки этих методов. Например:

- число расчетных схем и формул велико,

- структура формул во многих случаях сложна,

- расчетные формулы не универсальны и не охватывают сложные случаи фильтрации.

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

Отметим, что в ряде работ [9-12 и др.] в теории разработки газовых месторождений исследовались задачи однофазной и двухфазной фильтрации. В этих работах использовались и развивались методы определения фильтрационных параметров в вариационной постановке. Показателем качества идентификации параметров пласта выступал некото-

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

Таким образом, можно сказать, что причины всех существующих недостатков применяемых методов коренятся в исходном принципе построения этих методов, а именно

- в использовании решений краевых задач.

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

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

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

Метод модулирующих функций (М-метод)

Метод модулирующих функций (далее, М-метод) впервые был описан в работах [3, 4] в связи с решением задачи идентификации объектов автоматического управления. Это весьма гибкий и удобный своей компактностью метод. Этот метод позволяет определить параметры модели без использования решения краевых задач.

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

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

Проиллюстрируем силу М-метода на примерах.

Пример 1. Начнем с простого примера. Пусть задано простейшее дифференциальное уравнение

У = /(х) (У= Ф1 dx), (1)

где_Дх) - заданная на некотором интервале (а, Ь) непрерывная функция.

Умножим обе части уравнения (1) на гладкую (непрерывно дифференцируемую) на интервале (а, Ь) функцию р(х). Далее, проинтегрируем полученное уравнение по отрезку [х0, х1] с (а, Ь):

•ч -ч

|у р(х)ёх = |/(х)р(х)ёх . (2)

х0 х0

Применяя к левой части равенства (2) формулу интегрирования по частям, получим:

х1 х1

у( х)р(х)| хо -1у( хр( х11(х)р( х№. (3)

х0 х0

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

р( х0) = р( х1) = 0. Тогда из (3) имеем:

•Ч -Ч

|y(x)p(x)dx = |f (x)p(x)dx .

(4)

Соотношение (4) является интегральным аналогом дифференциального уравнения (1). Как видим, в (4) операция дифференцирования «перешла» от решения у(х) к «модулирующей» функции р(х) .

Пример 2. Пусть задано уравнение

/ = f (x) (y" = d2 y / dx2),

(5)

гдеfx) - заданная на некотором интервале (a, b) непрерывная функция.

Как в примере 1, умножим обе части уравнения (5) на дважды непрерывно дифференцируемую функцию р(x) и проинтегрируем по отрезку [x0, x1] с (a, b) :

xi xi

Iy"(x)p(x)dx = If (x)p(x)dx . (6)

x0 x0

Применим к левой части равенства (6) формулу интегрирования по частям два раза:

(7)

У'(x)P(x)| ,о - jy'(x)p'(x)dx =y'(x)p(x)| *

y(x)p(x)| x0 - jy(x)p(x)dx

Выбирая функцию р( x) такой, что р( x0) = р( x0) = 0 и р' (x0) = р( x0) = 0, из (6), (7) получим:

11 |y(x)p(x)dx = |f (x)p(x)dx .

(8)

Соотношение (8) есть интегральный аналог дифференциального уравнения (5). Пример 3. Пусть задано уравнение

У(n) = f (x) (y (n) = dny / dxn).

(n) _ An,

(9)

x

x

0

0

x

x

0

x

x

0

0

Поступая аналогично как в примерах 1 и 2, умножим обе части уравнения (9) на достаточно гладкую функцию р(х) (имеющей непрерывно производные до п порядка включительно) и проинтегрируем левую часть полученного соотношения по отрезку [х0, х1] с (а, Ь) . Выбирая «модулирующую» функцию р(х), удовлетворяющей условиям

Лр* = 0, = 0, (*=0,1,...,п-1)

(хк (хк

и применяя формулу интегрирования по частям п раз, получим интегральный аналог уравнения (9):

(-1) п | у( х) (х = 1 /(х)р( х)(х.

(хп

х0 х0

Определение постоянных коэффициентов в дифференциальных уравнениях

Теперь применим использованную в примерах 1-3 методику для решения простейших обратных коэффициентных задач.

Пример 4. Рассмотрим задачу определения коэффициента к в уравнении

У = ку (у = (у / (X). (10)

Уравнение (10) - это уравнение нормального размножения, если к>0, и - уравнение радиоактивного распада, если к<0.

Умножим обе части уравнения (10) на гладкую (непрерывно дифференцируемую) на интервале (а, Ь) функцию р(х) и проинтегрируем полученное уравнение по отрезку [х0, х1] с (а, Ь) . Имеем:

х1 х1

|у р(х)(х = |кур(х)(х . (11)

х0 х0

Применяя к левой части равенства (11) формулу интегрирования по частям, получим:

х1 х1

у(х)р(х)|х‘ - |у(х)р'(х)(х = |ку(х)р(х)(х . (12)

|х0

х0

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

р( х0) = р( х1) = 0. Тогда из (12) имеем:

х1 х1

- |ур'(х)(х = к |ур(х)(х. (13)

х0 х0

Из равенства (13) определяем коэффициент к:

х1

|ур'( х)(х

к = -х------------.

х1

|ур( х)(х

X

X

Пример 5. Определить коэффициент к в уравнении взрыва

У = ку2 (/ = ёу / ёх). (14)

Как и в примере 4, умножим обе части уравнения (14) на гладкую (непрерывно

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

уравнение по отрезку [х0, х1] с (а, Ь). Имеем:

х1 х1

|ур(х)ёх = |ку2р(х)ёх . (15)

х0 х0

Применяя к левой части равенства (15) формулу интегрирования по частям, и выбирая «модулирующую» функцию р( х) такой, что р( х0) = р( х1) = 0, получим

Л1 Л1

|ур(х)йх = к |у2р(х)ёх. (16)

Из равенства (16) определяем коэффициент к:

х1

|ур'( х)ёх

к=

|у 2р( х)ёх

х0

Пример 6. Определить коэффициент о2 в уравнении колебаний маятника

у" + о2 у = 0 (у" = й2у / ёх2). (17)

Поступая аналогично, как и в примере 2, получим

х1 х1

|урр(х)ёх = -о2 |ур(х)ёх . (18)

х0 х0

Из равенства (18) определяем коэффициент о2:

х1

|урр ( х)йх

2 х0

о =—0---------------.

х1

|ур( х)ёх

х0

Теперь продемонстрируем применение М-метода для определения коэффициентов в дифференциальных уравнениях с частными производными.

Пример 7. Определить коэффициент а2 в уравнении свободных малых колебаний струны

д2и 2 д2и гтт ттг чч

-гг- = а (и = и (х, г)). (19)

дг дх

Умножим обе части уравнения (19) на дважды непрерывно дифференцируемые функции р1 (х) и р2 (г) и проинтегрируем соответственно по отрезкам х0 < х < х1, г0 < г < г1:

х

х

0

х

0

х

г1 х1 д 2 и г1 х д 2и

—— р, (х)р2 (г)ёхёг = а2 [ [—— р1 (х)р2 (г)ёхёг.

2 1 2 2 1 2

г1 х1 2 дг2

00

дх

00

Уравнение (20) можно записать еще так:

г1 д2 дг

|р1 (х) |р2 (г)ёгёх = а2 р2 (г) |дд ^ р1 (х)ёхйг.

х0 г0

00

(21)

Выберем модулирующие функции р1 (х) и р2 (г) такими, чтобы р1 (х0 ) = р1 (х1 ) = 0 , Р2 (г0 ) = Р2 (г1 ) = 0 , р1 (х0 ) = р1 (х1 ) = 0 , Р2 (г0 ) = Р2 (г1 ) = 0 . (22)

гГ д 2и хг д 2и

Применим к интегралам I—— р2(г)йг и I—— р1(х)ёх по два раза формулу интег-

•’ дг * дх2

дх2

рирования по частям. Имеем:

г1_ д 2гг г1 ^2

дг

1 д 2и 1~дгГ Р2(г )ёг = Iи

Ьг д 2Р2 (г)

х1 2

дг2

ёг.

}р1( х)ёх = |и

д р1(х) дх2

(23)

Тогда из соотношений (21)-(23) получим

|р1 (х) |и д Р22(г) ёгёх = а2 |р2 (г) |и ° х) ёхёг, или

х1 2

д р1( х)

х0 г0 г1 х1

дг2 у24 дх2

г0 х0

г1 х1

| |ир1 (х)р2(г )ёхёг = а21 |ир!( х)р2 (г )ёхёг.

г0 х0 г0 х0

Из равенства (24) теперь определяем коэффициент а2 :

г1 х1

| |ир1 (х)р2(г )ёхёг

(24)

а

2 г0 х0

11

| |ир!( х)р2 (г )ёхёг

г0 х0

Пример 8. Определить коэффициент а в уравнении теплопроводности

а

д2и д2и +

дх2 ду2

ди

дг

(25)

Умножим обе части уравнения (25) на дважды непрерывно дифференцируемые функции р1 (х) и р2 (у) и непрерывно дифференцируемую функцию р3 (г). Полученное уравнение проинтегрируем соответственно по отрезкам х0 < х < х1, у0 < у < у1, г 0 < г < г1:

а

г1 у1 х д 2 и г‘ х у1 д 2 и

рз(г) |р2(у) I дТ р1(х)ёхёуёг + рз(г) |р(х) I -у р2 (у )ёуёхёг

г0 у0

у1 х1

г0 х0 у0

ду2

= |р2 (у) |р1 (х) 1 ^ Рз (г)ёгёхёу.

(26)

у0

00

х

0

0

х

х

0

Выбирая модулирующие функции (рх (х), р2 (у), р3 (‘) такими, что

р1(хо) = р1(х1) = 0, Р 2 (У 0 ) = р2(У1) = 0, Рэ(^0) = р3(Ь) = 0, Р!(Х0) = Р1(Х1) = 0, Р2(У0) = Р2(У1) = 0 ,

и применяя в (26) формулу интегрирования по частям к интегралам [д ^ р1(х)йХ и

•’ Эх

• ох2

х0

У\ д 2и . .. \ ди ...

I—— р2 (у)ау два раза, а к интегралу I-р3 (‘)а‘ - один раз, получим:

•’ ду2 ; д‘

у0 ^ ‘0

‘1 У1 х1

а

‘1_У1_Ж1_ *1.У1х1.

‘0 У0 х0 ‘0 У0 х0

= - jjjир1р2р3ёхёуЛ . (27)

‘0 у0 х0

Из равенства (27) определяем коэффициент а:

‘1 у1 х1

11 \ир1р2р3 ёхёуЛ

‘0 у0 х0

а =

‘0 у0 х0 ‘0 у0 х0

ч М-^-1 ч М-^-1

11 |ирр2р3 ёхёуЛ +11 \ирхр\ръ ёхёуЛ

Определение переменных коэффициентов в дифференциальных уравнениях

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

Рассмотрим примеры применения М-метода для определения переменных коэффициентов дифференциальных уравнений.

Пример 9. Определить коэффициент к(х) в уравнении

у = к(х)у х е (а, Ь) с Я . (28)

Предполагая функцию к(х) достаточно гладкой, разложим ее по формуле Тейлора: к (х) =к 0+к1 х + к 2 х2 +--------------+ кпхп + Я( х), (29)

где Я(х) - остаточный член.

Ограничиваясь первыми п+1 членами разложения и подставляя (29) в уравнение (28), получаем:

у' =к0у + к!ху + к2х2у + ••• + кпхпу . (30)

Умножим уравнение (30) на непрерывно дифференцируемые модулирующие функции рi (х) (I = 0,1,..., п) и проинтегрируем полученные уравнения по отрезку [х0, х1] с (а, Ь) :

І у (х)& = І (к к (() + ( хур0 (х) + к 2 х2 Уф0 (х) + • • • + кпхпур0 (х) );

х0 хо

х1 хі

І у' р (х)<іх = І (к оУ<Рі (х) + к ху^1 (х) + к 2 х2 урі (х) + • + кпхпурі (х));

х0 х0

хі хі

І у' Рп(х)іїх = Ік оуРп(х) + кі хурп(х) + к 2х 2 урп (х) + • + кпхПуРп(х))

х0 х0

Полученную систему уравнений (3і) можно переписать так:

хі п хі

Іу'рг (х)йх = ^кі Іхгург (х^ (і = 0,і,..., п).

(3і)

(32)

Выбирая модулирующие функции р (х) такими, что

р (х0) = р (х1) = 0 0' = n),

и применяя к левым частям уравнений системы (32) формулу интегрирования по частям, получим

- Іур0(х)^х = к0 Іур0 (х)дх + кі |хур0(х^х + к2 |х2ур0(х)$х +------к |хпур0(х)<іх;

х0 х0 х0 х0 х0

х хі хі хі хі

|урі(х)ёх = к0 |урі (х)ёх + кі |хурі (х)ёх + к2 |х2урі (х)ёх +-+кп |хпурі (х)ёх;

(33)

Л^ Л^ Л-1

Iур'п (х)Ох = к0 |урп (х)Ох + к1 |хурп (х)Ох + к2 |х2урп (х)Ох + ••• +кп |хпурп (х)ах.

х0 х0 х0 х0 х0

Система (33) есть система линейных алгебраических уравнений относительно неиз-

вестных к0, к1, к, кп.

Решая эту систему одним из известных методов (например, методом исключения Гаусса), найдем искомые параметры - коэффициенты к (' = 0,1,...,п) в разложении (29).

Отметим, что независимость уравнений системы (33) можно достичь, как было отмечено выше, например, варьируя вид модулирующих функций р (х) .

Пример 10. Определить коэффициент к(х) в уравнении

у' + к (х) у' + £ (у) = / ( х).

(34)

Здесь предполагается, что функции У(х) и g(y) непрерывны, а функция к(х) достаточно гладкая. Уравнение (34) - это уравнение движения маятника в сопротивляющейся среде с переменным коэффициентом трения к(х) и с восстанавливающей силой g(y) под действием силы _Дх).

Как и в примере 9, разлагая функцию к(х) по формуле Тейлора и умножая уравнение (34) на дважды непрерывно дифференцируемые модулирующие функции р (х) (' = 0,1,..., п), а затем интегрируя полученные уравнения, получим:

г =0

х

х

0

х

х

х

х

х

0

0

к) IуР(х)^х + кі |ху'р(х)ёх +-----------+ кп |хпу'р()(х)ёх =

х0 х0 х0

хі хі хі

= I/(х)р0(х)ёх - I£(х)р0(х)ёх -1ур0(х)ёх; х0 х0 х0

хі хі хі к0 |урі(х)ёх + кі |хурі(х)ёх +-+ кп |хпурі(х)ёх =

х0 х0 х0

x х х

= I/(х)Рі(х)ёх - I£(у)Рі(х)ёх - Iу'Рі(х)ёх;

х0 х0 х0

xi хі хі

к0 Iу'рп (х)ёх + кі Iху'рп (х)ёх + ••• + кп Iхпу рп (х)ёх =

х0 х0 х0

хі хі хі

= I/(х)Рп(х)ёх - I£(у)Рп (х)ёх - Iу'Рп(х)ёх

х0 х0 х0

Выбирая модулирующие функции рі (х) такими, что

Р (х0) = Р (хі) = ^ р (х0) = р (хі) = 0 (і = п)

хі

и применяя к интегралам |у0рп (х)ёх формулу интегрирования по частям два раза, а к

х0

хі

интегралам кі |хгу'рёх (і = 0,і,...,п) - один раз, получим

х0

хі хі хі

к0 Iур0(х)ёх + кі Iхур0(х)ёх +-----+ кп Iхпур0(х)ёх =

х0 х0 х0

х х х

Iур0(х)ёх - If (х)Р0(х)ёх + I£(ур0(х)ёх; х0 х0 х0

х х х

к0 Iурі(х)ёх + кі Iхурі(х)ёх +-----------+ кп Iхпур[(х)ёх =

х0 х0 х0

x х х

= Iурі(х)ёх- I/(х)Рі(х)ёх + I£(у)Рі(х)ёх; (35)

х0 х0 х0

xi хі хі

к0 Iурп (х)ёх + к Iхурп (х)ёх + ••• + кп Iхпурп (х)ёх =

х0 х0 х0

хі хі хі

= Iур"п(х)ёх- I/(х)Рп(х)ёх + I£(у)Рп(х)ёх

х

х

х

0

0

Получили, как и в примере 9, систему алгебраических уравнений относительно неизвестных к (/ = 0,1,...,п). Решая ее, находим искомые коэффициенты к разложения функции к(х).

Рассмотрим теперь более содержательный пример.

Пример 11. Определить коэффициент к(ху) в уравнении фильтрации

д_

дх

к (х, у)

дР2 дх

+ -

д_

дУ

к(х, у)

дР2 дУ

др

ді

(36)

Здесь Р - давление в точке пласта с координатами (х, у) в момент времени і, к(ху) -коэффициент проницаемости, /I - коэффициент динамической вязкости газа.

Как и выше, разложим функцию к(х,у) по формуле Тейлора

к (х, у) =к 0+к1 х + к2 у + к3 ху + к4 х2 + к5 у2 + • После подстановки (37) в уравнение (36) получим

(37)

[о( Р !)'х + [ х( Р !)'х + [ у( Р !)'х + к, ху (Р !)'х + к4 х 2( Р !)'х + к5 у 2( Р 2 )'х + ...] +

_3_

дх

+ д- [о( Р 2)'у + к1х(Р 2)у. + [ у( Р 2)у. + к3 ху(Р 2)у + к4х ’( Р ’Гу + к5у ЧР 2)у + ••.] =

о дР = 21—.

^ ді

Далее имеем

{)' + к,(х(Р% )', + к 2 у( Р 2)1 +...}+

+ {ко(Рк + к(Р% + к,(у(Р2)у)У + ...}= 2}2)і '

(38)

(39)

Умножим обе части уравнения (39) на гладкие функции р(х),р(у),р^)

2 1 (/ = 0,1,...,п) (причем р(х) ир(у) - функции класса С , р(^) - функции класса С ) и

проинтегрируем полученное уравнение соответственно по отрезкам х0 < х < х1, у0 < у < у1, 10 < t < ^ . Имеем:

іі уі

іо уо

+

+ к

||р0 )Фг(у)ФЖ • к0 |(Р 2) ІР( ^ + к1 |(х(Р 2)'х )Р( х№ + к 2 |у(Р 2) 1 Р (х)^

’о |_ х0 х0 х0

х1 уі уі уі

I \фі (і Р (х)^х^і • ко I(Р 2)"уу Рг (У¥У + к1 Iх(Р 2) Уу Р (у¥у + к 2 I (у(Р 2 )'у )'у Р (у¥у

іо хо _ уо уо

у1 х1 і1

= 2і 11р(у)р(х)ёхёу |Рр(і)Л (і = о,1,...,п).

+

Уо

+ к

уо хо

Выбирая модулирующие функции Рі (х), Рі (у), Рі (і), удовлетворяющие условиям

х

х

х

о

Рг (хо) = Рг (х1) = ^ Рг (уо) = Рг СЮ = ^ Рг (іо) = Рг (і1 ) = ^

Р (хо) = Р (х1) = 0, Р (Уо) = Р (У1) = ^ (і = п)

и применяя формулу интегрирования по частям по два раза к тем интегралам, в которых присутствуют (Р2)X., (Р2), а к интегралам, в которых присутствуют (Р2)X,

(Р2) у, (Р2)" - один раз, получим

'і Уі

'o У0 h_ Хі

+

хі хі , хі

I p (t)p (У)dydt ■ k0 IP2p"(x)dx + kl IP2 (xp'(x)) dx + k2 |yP2p"(x)dx + ...

’o L x0 x0 x0

h Хі Уі Уі Уі ,

I p(t)P(x)dxdt ■ k0 IP2p(УУу + kl IxP2p(У )dy + k2 IP2 (yp(У )) dy

'o Xo L У0 У0 У0

'і Уі xl

-2p 11 IPp (x)p (y)p(t)dxdydt, (i = 0,1,..., n).

'o У0 Xo

Систему уравнений (40) можно переписать так:

'і Уі xl

-0 ff [P2[(x)p(y) + p"(y)p(x)"|p(t)dxdydt +

+

(40)

ko\\]P 2 \( x)P (y) + p( yp (x)] (t )dxdydt

to У0 Xo h Уі xl

+kl UfP 2 [(xp( x))p (y)+xp( y )p(x )](t )dxdydt+

to У0 Xo h Уі Хі

+ k2 ]]\P 2 [(y p (y ))p ( x) + yp( x)P (y)]p (t )dxdydt +

to У0 Xo

h Уі xl

+ ••• = -2^ ff [Pp(x)p(y)p(t)dxdydt, (i = 0,1,..., n).

(41)

Определение коэффициента k(xy) заключается в нахождении постоянных разложения к 0, к,..., кп,.... Последние находятся из системы линейных алгебраических уравнений (41). Необходимую точность аппроксимации функции k(xy) соответствующим многочленом Тейлора можно обеспечить увеличением числа членов разложения в (37).

Замечание. В приведенных выше примерах в качестве модулирующих функций можно использовать, например, функции вида <р( x) = (x - x0)(x1 - x) в уравнениях, в которых присутствуют производные первого порядка и (р( x) = (x - x0)2( x1 - x)2 - в уравнениях с производными второго порядка.

Примеры модулирующих функций

Универсальной модулирующей функцией является, например, функция:

(р(x) = а(x) exp ( - "У - ) , (42)

(x xo )(xi x)

где a(x), ^(x) - функции, непрерывные вместе со своими производными вплоть до граничных точек xo и xj.

Функция (42) вместе со своими производными любого порядка равна нулю на концах интервала (x0, x1) (в этом легко убедиться непосредственным вычислением). С помощью функций a(x), ^(x) можно придать любой «вес» экспериментальным функциям в области наблюдения. Схематически график функции (42) изображен на рис. 1.

Рис. 1. Типичный график модулирующей функции.

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

Для дифференциальных уравнений фильтрации (т.е. уравнений, как правило, второго порядка) модулирующие функции для промежутка [х0, х1] могут быть взяты весьма простого вида, например:

п

к( х) = а( х)( х - х0 )2 (х1 - х)2, к( х) = а( х) sin2------(х - х0),

(х1 - х0)

где а(х) - любая непрерывная вплоть до точек х0, х1 функция. Множители

п

(х - х0)2( х1 - х)2 и sin2-(х - х0) обеспечивают выполнение условий (22). Мно-

(х1 - х0)

житель же а(х) обеспечивает необходимую форму модулирующей функции на отрезке [х0, х1]. В частности, можно взять а(х)=1.

Для аргумента, по которому уравнение имеет первый порядок, модулирующие функции могут быть выбраны следующими:

к( х) = а( х)( х - х0)(х1 - х), к( х) = а( х) sin —х°).

(х1 - х0)

Заключение

Исходя из вышесказанного, можно сделать следующие выводы.

1. М-метод позволяет перейти от некорректной операции дифференцирования экспериментальной функции к корректной операции интегрирования.

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

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

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

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

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

зин Р.Н. Моделирование процессов нефтегазодобычи. Нелинейность, неравновес-ность, неопределенность. - Москва-

: -

ваний, 2004. - 368 с.

2. ., . . моделирование и обратные задачи // Вестник Адыгейского государственного университета. - 2008. - № 4(32). - С.18-24.

3. Loeb J., Cahen G. Extraction, a partik des enregistrements de mesures, des parametres dynamiques d um system. // Automatisme. -1963. - № 12. - PP. 17-28.

4. Loeb J., Cahen G. More about process identification. // Trans. on Automatic Control. - 1965. - PP. 359-361.

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

6. ., . . геофильтрации. - М.: Недра,1976. - 407 с.

7. . ., . .

на ЭВМ унифицированных алгоритмов . . . - . . . -НИИгидротех. и мелиор., 1976, вып. 9. - С. 111-114.

8. Юдин АЛ., Юдина ОХ. Расчет фильтрационно-ёмкостных параметров по промысловым данным эксплуатации газового месторождения // Термодинамика кооперативных процессов в гетерогенных средах. Тюмень, 1985. - С. 80-85.

9. Закиров С.Н., Коршунова Л.Г., Нанивский

..

теории разработки газовых месторожде-

- , . -ка и эксплуатация газовых и газоконден-

. . 1 2. - .,

1975.

10. -

ботки газовых месторождений. / Закиров С.Н., Васильев В И., Гутников АЛ. и др. -М.: Недра, 1984. - 295 с.

11. . .

разработки газовых и газоконденсатных месторождений. М.: Недра, 1989. - 334 с. 12. . ., . . -кация параметров газовых залежей при газовом и водонапорном режимах разработки. - Обз. инофрм. ВНИИЭгазпрома. Сер. Разработка и эксплуатация газовых и газоконденсатных месторождений. - М., 1990.