Размещено на http://www.allbest.ru/

ИДЕНТИФИКАЦИЯ ВОДОПРОВОДИМОСТИ ПЛАСТА В СТАЦИОНАРНОМ ПОТОКЕ ПОДЗЕМНЫХ ВОД

М.У. Мурзакматов

Рассматривается применение метода регуляризации А.Н. Тихонова к идентификации водопроводимости пористой среды в двумерном стационарном уравнении фильтрации подземных вод.

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

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

Установившееся движение подземных вод в водоносном пласте, ограниченном сверху и снизу непроницаемыми прослойками, описывается уравнением [1]:

(1)

с граничным условием

,(2)

Где Т=Т(х,у) водопроводимость пласта (м²/сут); H=H(x,y) напорная функция (м); f(x,y) функция источников и стоков (м/сут); и заданные функции; D область фильтрации в плане, Г ее граница.

В данной задаче, кроме условия (2), задаются так называемые внутренние граничные условия

, (3)

, (4)

о которых говорилось выше.

Задача заключается в определении функции Т(x,y) из уравнения (1) при соблюдении условий (2) (4). Поскольку значения напоров задаются с определенной погрешностью и в недостаточном объеме, то задача нахождения коэффициента уравнения (1) является некорректной, поэтому для ее решения применяем метод регуляризации А.Н. Тихонова [2]. Задача сводится к нахождению функции Т(x,y), сообщающей в области D минимум функционалу [3]

,(5)

где вариация функции Т(х,у); параметр регуляризации; H_i(T)расчетные значения напоров, которые находятся как решение задачи (1), (2).

Задача (1), (2) решается методом конечных элементов [4]. Область фильтрации D разбивается на треугольные элементы таким образом, чтобы точки, в которых заданы экспериментальные значения и , совпали с вершинами элементов. В этих точках используются условия (3) и (4), а в остальных узлах начальные приближения функций Н(х,у) и T(х,у) должны удовлетворять условиям и соответственно. В элементе (е) с вершинами i, j,k функция Н(х,у) выражается формулой [5]

, (6)

где

,

.

Суммируя равенство (6) по всем элементам (е), получаем формулу

,(7)

где n число всех узлов сетки.

Применяем к задаче (1), (2) принцип Галеркина:

Используя в двойном интеграле формулу Грина, приходим к системе уравнений

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

(8)

где

,

. (9)

Систему (8) решаем одним из точных или приближенных методов.

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

Теперь займемся минимизацией функционала (5). При каждом наборе значений функции Т(х,у) получаем вполне определенные значения функции Н(х,у), т. е. имеем оператор Н(Т), определенный алгоритмически по формулам (1), (2), (8), (9). Этот оператор в общем случае является нелинейным. Линеаризуем его следующим образом

, (10)

где значение функции Т в точке (x_s,y_s), полученное в предыдущей итерации; R₂() остаточный член разложения. Подставляя (10) в (5) и используя необходимое условие минимума функции многих переменных, получаем систему линейных алгебраических уравнений

или

(11)

с коэффициентами

,

и с правыми частями

Здесь

Матрицы систем (8) и (11) являются симметричными и имеют диагональное преобладание и они легко решаются методом Гаусса. Производная является операторной, т. е.

,

а частные производные аппроксимируются разностными отношениями

,

где номер итерации.

Вычислительная процедура осуществляется в следующем порядке. Используя начальные значения Н и Т в качестве нулевого приближения, решается задача (1), (2) и определяется первое приближение Н⁽¹⁾ . Затем, придавая приращение функции T(х,у), находим второе приближение Н⁽²⁾ . Это дает возможность приближенно определить производную и решить систему (11) при некотором значении параметра регуляризации . Итерация по проводится до установления фильтрационного процесса. Если при этом полученные значения напоров в пределах ошибок не совпадут с данными экспериментальными значениями , то итерация проводится по параметру , выбор которого может быть осуществлен методом невязок [2].

Работа алгоритма проверена на решении следующей тестовой задачи, рассмотренной в работах [6,7]. Областью фильтрации D является круг , который разбит на 54 треугольника (элемента) с максимальной длиной сторон . Число узлов сетки (вершин треугольников)37, из них 18 граничных. В этой области заданы функции Н(х,у)=х²+у² +5, Q(x,y)=0, f(x,y)=40(2x²+2y²+1). Искомой функцией является T(x,y)=10(x²+y²+1). Число узлов, в которых задаются экспериментальные (точные) значения функций H(x,y), равно: p= 37, 22, 5. В каждом из этих случаев задаются точные значения искомой функции в q=17, 9, 5 точках.

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

Таблица 1

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

В табл. 2 решение данной задачи сравнивается с соответствующими результатами, полученными другими методами [6,7]. Следует отметить, что идентификация параметров водоносных горизонтов с использованием методов малых возмущений [6] и регуляризации является устойчивой процедурой, что очень важно при проведении гидрогеологических расчетов в реальных условиях

Таблица 2

Сравнение с результатами, полученными другими методами

Примечание: aпогрешности метода малых возмущений [6];

bпогрешности, полученные в работе [7];

спогрешности метода регуляризации, полученные при p=5.

регуляризация водопроводимость пористый фильтрация

Литература

Полубаринова Кочина П.Я. Теория движения грунтовых вод. М.: Наука, 1977.664 с.

Тихонов А. Н., Арсенин В. Я. Методы решения некорректных задач. М.: Наука, 1974 233 с.

Мурзакматов М. У., Мамыров Ж. М. Об идентификации коэффициента фильтрации в неоднородном водоносном горизонте //Вестник ИссыкКульского университета, -№ 3, 1999. с. 7377.

Джаныбеков Ч. Математическое моделирование движения грунтовых вод в многослойных средах. Фрунзе: Илим, -1982. 280 с.

Сегерлинд Л. Применение метода конечных элементов. М.: Мир, 1979. 392 с.

Мурзакматов М.У., Байболотов Б. А. Об идентификации водопроводимости напорного потока методом теории возмущений //Вестник ИссыкКульского университета, № 9 Каракол, 2003.С. 4954.

Мурзакматов М. У., Исабеков К.А. Об идентификации параметров планового фильтрационного потока //Вестник ИссыкКульского университета, № 9 Каракол, 2003.С. 2633.

Размещено на Allbest.ru