Автор И.В.Авилина
Программы NASA предназначена для обработки данных откачек в одно- и двухпластовых водоносных системах с использованием компьютерных методов. В программе реализованы схемы упругого гомогенного и упругого гетерогенного перетекания в разделяющих и покровных пластах. Решения прямых задач получены Шестаковым В.М. (Шестаков, Невечеря, 1998; Шестаков, 2011) с использованием преобразования Лапласа - Карсона и обратным переходом к оригиналам по методу Л.К. Гохберга. В покровном пласте и в разделяющем прослое в двухпластовой системе принимается вертикальное перетекание в упругой гетерогенной среде, состоящей из проводящих каналов и слабопроницаемых блоков. Для схемы упругого гомогенного перетекания предусмотрена процедура автоматизированного поиска параметров моделей путем минимизации функции качества. Для схемы гетерогенного упругого режима автоматизированный поиск параметров модели не предусматривается, возможно имитационное моделирование. Понижением в покрове при откачке пренебрегаем.
Критерием нахождения искомых геофильтрационных параметров в программе служит минимум функции качества, определяемой как сумма квадратов разности между модельными и измеренными понижениями по всем расчетным моментам времени и всем наблюдательным скважинам. Минимум функции качества находится путем многократных решений прямой задачи при изменении начальных величин искомых параметров. В программе NASA реализован поиск минимума функции качества методом Розенброка.
Остановка расчета в программе происходить по заданному количеству вычислений функции качества, поэтому рекомендуется проводить минимизацию несколько раз до получения устойчивого минимума.
Для схемы упругого гомогенного перетекания искомые параметры обратной задачи в системе -1 пласт с покровом : T1 - проводимость 1 пласта, μ1 - -упругая емкость 1 пласта, k0 - коэффициент фильтрации покрова, ηk (etak) - упругая емкость проводящих каналов. Количество параметров может меняться от 1 до 4.
Для схемы упругого гомогенного перетекания искомые параметры обратной задачи в системе -2 пласта со слабопроницаемым прослоем: T1 - проводимость 1 пласта, μ1 - упругая емкость 1 пласта, k0 - коэффициент фильтрации слабопроницаемого слоя, ηk (etak) - характеризует упругий гомогенный режим в прослое, T2 - проводимость 2 пласта, μ2 - упругая емкость 2 пласта. Количество параметров может меняться от 1 до 6.
Параметры ищутся из условия минимума функции качества:
 ,
где xi - координаты наблюдательных скважин, tj - моменты времени (от начала откачки). Начало прямоугольных координат x=0 -местоположение центральной скважины, Sm - численно- аналитическое решение, Se -натурные наблюдения. Метод минимизации - Розенброка.
Для схемы упругого гетерогенного перетекания решение обратной задачи по поиску параметров не предусмотрено. В этом случае в исходных данных вместо параметра ηk (etak) задается параметр ηb (etab)- упругая емкость слабопроницаемых блоков ( другой параметр режима строго фиксирован в программе и равен 1).
Скачать программу. Содержимое архива: anno3.doc (аннотация), ins3.doc (инструкция), nasa.bat, nasa.exe(сама программа), nasa.bat(nasa.exe+vgamma.com, при несовпадение русификатора), upr1.dat (пример файла данных), upr1.res (пример файла результатов), другие примеры - upr1a.dat , upr1a.res, upr2.dat , upr2.res, upr3.dat , upr3.res, upr4.dat , upr4.res
Литература.
1. Шестаков В.М., Невечеря И.К. Теория и методы интерпретации опытных откачек. - М.: МГУ. 1998. 160 с.
2. Шестаков В.М., 2011 - в рукописи издательского проекта РФФИ 2011 года "Методы расчета опытных откачек в водоносных пластах с перетеканием",
авторы - Шестаков В.М., Невечеря И.К., Авилина И.В.
|
