О моделировании подтопления и дренирования земель с помощью фортрана
Автор: Сологаев В.И.
Журнал: Вестник Омского государственного аграрного университета @vestnik-omgau
Рубрика: Науки о земле
Статья в выпуске: 1 (25), 2017 года.
Бесплатный доступ
Рассмотрена методология моделирования фильтрации подземных вод в почвах и грунтах при защите от подтопления в городском строительстве и сельском хозяйстве. Целью исследования является использование численного метода конечных разностей для моделирования движения подземных вод с помощью алгоритмов, написанных на языке программирования фортран. Научная новизна статьи в особых приемах и алгоритмах, содержащихся в разработке автора на языке программирования фортран. Проведена привязка теоретических выкладок к практике гидрогеологических работ в условиях подтопляемых земель. Исследование проведено с теоретическими выкладками на примере города Калуги с применением авторской программы на фортран, специально разработанной для данной темы. Рассмотрена напорная нестационарная фильтрация подземной воды в песке при откачке из 13 скважин по схеме неограниченного в плане водоносного пласта. Для моделирования принят участок размером в плане 500 на 500 метров. Компьютерная модель водопонижения методом конечных разностей использует явную схему. Шаги пространственной сетки модели по горизонтали и вертикали приняты равными. Шаг времени на модели принят максимально возможный с учетом критерия устойчивости счета по явной схеме. Расчетное время модели на компьютере с тактовой частотой 2.13 GHz получалось от нескольких секунд до 10 минут. Результаты расчета записаны в выходные текстовые файлы. По ним построены карты напоров подземных вод при работе системы из 13 водопонизительных скважин. В отличие от традиционного алгоритма вычислительной математики с циклом for-end, для решения разностной схемы применен вложенный цикл фортрана do-end. Тесты показали увеличение скорости расчета модели на три порядка, то есть примерно в тысячу раз.
Фильтрация, подземная вода, грунты, защита от подтопления, моделирование, алгоритмы, фортран
Короткий адрес: https://sciup.org/142199328
IDR: 142199328
Текст научной статьи О моделировании подтопления и дренирования земель с помощью фортрана
Фильтрация подземных вод количественно может быть рассчитана с помощью закона Дарси. Необходимо совершенствовать методологию моделирования движения воды в почвах и грунтах, отталкиваясь от изначальных теоретических основ с опорой на опытные данные. Наиболее современный способ моделирования – использование численных методов. Автор отдает предпочтение методу конечных разностей (МКР). Далее показано использование языка программирования фортран для решения поставленных практикой задач моделирования подтопления и дренирования.
Объекты и методы
Компьютерное моделирование водопонижения из 13 скважин выполнено для стадии строительства многофункционального спортивного комплекса «Дворец спорта» в городе Калуге. Научно-исследовательская работа проведена на основании исходных данных, изложенных в техническом отчете об инженерно-геологических изысканиях 11267-ИГИ ОАО «КалугаТИСИЗ» 2014 г. Все необходимые данные предоставил советник генерального директора территориального проектного института ТПИ Омскграж-данпроект, кандидат технических наук Раис Шакирович Абжалимов.
Площадка изысканий расположена в городе Калуге по улице Ленина, на территории существующего стадиона. На момент проведения изысканий площадка занята под футбольное поле, трибуны, административное помещение и торговые павильоны. Рельеф ровный, спланированный. Отметки рельефа 193,00–195,60 м. Перепад – 2,6 м. Климат умеренно континентальный. Инженерно-гидрогеологические данные схематизированы по рекомендациям справочного пособия к СНиП [1]. За условный (относительный) ноль принята средняя отметка низа водонасыщенного песка с абсолютной отметкой 180,25 м. Коэффициент фильтрации песка k = 2 м/сут. Его средняя мощность M = 6 м. Выше песка до поверхности земли расположен суглинистый грунт толщиной 7 м с коэффициентом фильтрации ничтожно малой величины 0,00003 м/сут. Уровень подземных вод (УПВ) залегает на 2 м ниже уровня земли. Водонасыщенный песок подстилается водоупорными глинистыми грунтами.
Средняя глубина водопонизительной скважины принята 20 м. Она проходит с поверхности земли 7 м суглинистого грунта, ниже – 6 м водонасыщенного песка и заглубляется на 5 м в водоупорные глинистые грунты. Общая глубина скважины около 20 м, в том числе 2 м отстойник. Напор в скважине при водопонижении принят постоянным на отметке 180,25 м. В схематизированной схеме водоносного пласта это соответствует низу песка при контакте с водоупором. Скважина гидродинамически совершенная, ее фильтр принят на всю мощность М = 6 м водонасыщенного песка.
Компьютерная модель построена методом конечных разностей (МКР). Рассмотрена напорная нестационарная фильтрация подземной воды в песке при откачке из 13 скважин по схеме неограниченного в плане водоносного пласта. Напор в скважине Нс = 0 м. Начальный естественный напор в подземных водах Не = 13 м. Мощность во- доносного песчаного пласта М = 6 м. Коэффициент фильтрации песка, по данным изысканий, k = 2 м/сут. Коэффициент водоотдачи песка принят по Бецинскому [2]
ц = 0 ,117 7 k = 0 , 117^2 = 0 , 13 . (1)
Для моделирования принят участок размером в плане 500 на 500 метров.
Схема расположения водопонизительных скважин, предложенная Р.Ш. Абжали-мовым, показана на рис. 1. На расстоянии около 15 м от скважин 1-го и 5-го верхнего ряда расположены, соответственно, два 5-этажных жилых дома (5КЖ) – левый и правый. Для них на компьютерной модели необходимо было получить информацию о понижении УПВ при водопонижении.
Первая контрольная точка 1КТ по определению пониженного уровня подземных вод (УПВ) находится в центре котлована здания бассейна (рис. 1). Вторая контрольная точка 2КТ – в центре котлована здания дворца.
Компьютерная модель водопонижения методом конечных разностей (МКР) использует явную схему: подробности изложены в нашей монографии «Фильтрационные расчеты и компьютерное моделирование при защите подтопления в городском строительстве» (Омск, СибАДИ, 2002).

Рис. 1. Схема расположения водопонизительных скважин
Радиус водопонизительной скважины принят R с = 0,1 м. Шаги пространственной МКР-сетки по оси х (по горизонтали) и оси у (по вертикали) приняты равными, соответственно, dx = dy = dL = 0,5 м. Это связано с тем, что при задании постоянного напора в скважине на уровне низа водонасыщенного песка Н с = 0 м (относительная отметка) R c и dL связаны соотношением вида
R c = 0,2 ∙ dL = 0,2 ∙ 0,5 = 0,1 м. (2)
Шаг времени на модели dt принят максимально возможный с учетом критерия устойчивости счета МКР-модели по явной схеме и вычислен:
, ц - dL 0 , 13 • 0 , 5
dt =-------=---------= 0 , 000677 сут = 58,5 сек ~ 1 мин.
4 • к • M 4 • 2 • 6
При размерах участка моделирования в плане 500 х 500 м при пространственном шаге dx = dy = dL = 0,5 м получается плановая модель фильтрации подземных вод к 13 водопонизительным скважинам, работающим в режиме постоянной откачки с фильтрационным напором Н с = 0 м. Количество узлов модели превышает миллион, так как 1001 ∙ 1001 = 1 002 001.
Заметим, что известная компьютерная программа геологической службы США ModFlow PMwin, версия 5.3 которой официально сертифицирована для применения в России, способна применять МКР-модели с максимальным количеством узлов не более 250 тыс. Это для целей нашего моделирования водопонижения с количеством узлов более 1 млн непригодно.
Поэтому в плане инновации решено построить МКР-модель водопонижения с оригинальной разработкой на языке программирования Fortran 95 [3]. С целью независимости от иностранных технологий (импортозамещение) использовано свободное программное обеспечение с 64-битной операционной системой Linux. Компилятор фортрана gfortran. Исходный код программы записан в файле Bass13cKB.f95 и скомпилирован в исполняемый файл Bass13cKB командой gfortran–std=f95 Bass13cKB.f95 Bass13cKB
Время водопонижения на модели рассчитано с такими значениями (с учетом организации строительства):
0,1 сут; 1 сут; 5 сут; 10 сут; 15 сут; 30 сут; 45 сут; 60 сут.
Время расчета модели на компьютере с тактовой частотой 2,13 GHz получалось от нескольких секунд до 10 мин. Результаты расчета записаны в выходные текстовые файлы. По ним построены карты напоров подземных вод при работе системы из 13 водопонизительных скважин, что выборочно показано на рис. 2–4.
1-е сутки

Не = 13 м, Нскв = 0 м, М = 6 м, Rскв = 0,1 м
K = 2 м/сут, mu = 0,13, t = 1 сут, 500х500 м2, шаг сетки 0,5 м 13 скважин (рис. 1)
Понижение УПВ = 0,139 м (контр. точка 1КТ)
Понижение УПВ = 0,012 м (контр. точка 2КТ) Понижение УПВ под левым 5КЖ = 0,157 м
Понижение УПВ под правым 5КЖ = 0,670 м
Рис. 2. Прогнозная карта водопонижения на 1-е сутки
5-е сутки

Не = 13 м, Нскв = 0 м, М = 6 м, Rскв = 0,1 м
K = 2 м/сут, mu = 0,13, t = 5 сут, 500х500 м2, шаг сетки 0,5 м
13 скважин (рис. 1)
Понижение УПВ = 2,746 м (контр. точка 1КТ)
Понижение УПВ = 1,495 м (контр. точка 2КТ)
Понижение УПВ под левым 5КЖ = 1,303 м
Понижение УПВ под правым 5КЖ = 2,423 м
Рис. 3. Прогнозная карта водопонижения на 5-е сутки
30-е сутки

Не = 13 м, Нскв = 0 м, М = 6 м, Rскв = 0,1 м
K = 2 м/сут, mu = 0,13, t = 30 сут, 500х500 м2, шаг сетки 0,5 м
13 скважин (рис. 1)
Понижение УПВ = 9,056 м (контр. точка 1КТ)
Понижение УПВ = 7,366 м (контр. точка 2КТ)
Понижение УПВ под левым 5КЖ = 4,439 м
Понижение УПВ под правым 5КЖ = 5,934 м
Рис. 4. Прогнозная карта водопонижения на 30-е сутки
Водопритоки в скважины в течение водопонижения сведены в таблицу по временным периодам. С течением времени водопонижения водопритоки в скважины уменьшаются. Понижение уровня подземных вод (УПВ) увеличивается. Разумеется, при условии постоянной откачки воды из скважин, которая должна быть автоматизирована. Уровень воды в скважинах рекомендовано поддерживать постоянным на уровне низа водонасыщенного песка. Фильтры скважин расположить в пределах водонасыщенного слоя грунта, то есть скважины должны быть гидродинамически совершенными.
Прогнозные притоки подземных вод в скважины 1–13 (рис. 1), м3/сут
№ скв. |
Время от начала откачки при водопонижении, сут |
|||||||
0,1 |
1 |
5 |
10 |
15 |
30 |
45 |
60 |
|
1 |
104,0 |
63,0 |
41,4 |
30,7 |
24,6 |
16,0 |
12,4 |
10,4 |
2 |
104,0 |
62,7 |
36,8 |
24,0 |
17,5 |
9,6 |
6,9 |
5,6 |
3 |
104,0 |
62,7 |
36,5 |
23,0 |
16,1 |
8,2 |
5,7 |
4,5 |
4 |
104,0 |
62,7 |
36,8 |
23,9 |
17,3 |
9,2 |
6,5 |
5,2 |
5 |
104,0 |
63,0 |
41,3 |
30,5 |
24,2 |
15,0 |
11,3 |
9,3 |
6 |
104,0 |
63,0 |
41,3 |
30,5 |
24,1 |
14,9 |
11,2 |
9,2 |
7 |
104,0 |
62,7 |
36,6 |
23,3 |
16,2 |
8,0 |
5,4 |
4,2 |
8 |
104,0 |
62,7 |
36,1 |
21,6 |
14,0 |
5,9 |
3,7 |
2,8 |
9 |
104,0 |
62,7 |
36,0 |
21,4 |
14,0 |
6,0 |
3,8 |
2,9 |
10 |
104,0 |
62,9 |
38,7 |
25,6 |
18,5 |
9,8 |
6,8 |
5,4 |
11 |
104,0 |
62,4 |
38,1 |
27,4 |
21,5 |
13,2 |
9,8 |
8,1 |
12 |
104,0 |
62,4 |
40,3 |
30,5 |
24,8 |
16,1 |
12,3 |
10,2 |
13 |
104,0 |
63,3 |
46,1 |
36,4 |
29,4 |
18,1 |
13,2 |
10,7 |
Анализируя прогнозные карты водопонижения (рис. 2–4), нетрудно обнаружить: с течением времени работы скважин, из которых непрерывно ведется откачка, воронки кривых депрессий пониженного УПВ начинают сливаться в общую воронку водопонижения. Это особенно заметно на 30-е сут работы водопонизительной системы (рис. 4). То есть подтверждается гипотеза «большого колодца», принятая для системы близко расположенных скважин [4]. Водопритоки скважин (таблица) в начале водонижения практически одинаковые, скважины вначале не влияют друг на друга. Затем, в ходе непрерывной откачки воды, скважины начинают влиять все больше и больше, причем у внутренних скважин типа 8 и 9 (рис. 1) водоприток уменьшается в десятки раз.
Рассмотрим более подробно приемы программирования и алгоритмы, примененные в ходе данной научно-исследовательской работы.
В самом начале кода программы было применено новое замечательное ключевое слово фортрана allocatable [5–7] для объявления массива напоров, а именно:
real, allocatable :: H0(:,:) real, allocatable :: H(:,:)
где Н0 – начальные напоры непониженного УПВ; Н – напоры УПВ при водопонижении.
Это позволило записывать программный код компактно, без применения громоздких циклов, как, например, обычно принято в ФОРТРАН 77.
МКР-схема была принята явная, для нестационарной фильтрации подземных вод. В терминах вычислительной математики мы применили явную разностную схему для решения дифференциальных уравнений фильтрации параболического типа, то есть для нестационарного движения подземных вод в плане в рамках гипотезы Буссинеска.
В отличие от традиционного алгоритма вычислительной математики с циклом for-end [8, с. 185], для решения упомянутой разностной схемы применен вложенный цикл фортрана do-end. Суть данного приема в том, что параметр времени моделирования вынесен из тела массива. Этот особый алгоритмический прием позволил избавиться от ограничения памяти персонального компьютера (ПК), т.е. фактически предохраниться от overflow-переполнения оперативной памяти ПК. В результате алгоритм стал быстрым. Тесты показали увеличение скорости расчета модели на три порядка, т.е. примерно в тысячу раз.
Заключение
В представленной работе предпринято исследование по совершенствованию методологии моделирования движения подземных вод. Показано, что с помощью фортрана с учетом особых алгоритмов можно существенно ускорить расчет компьютерных моделей водопонижения при защите от подтопления земель. Это должно повысить надежность проектирования защиты от подтопления в городском строительстве и сельском хозяйстве.
V.I. Sologaev
Omsk State Agrarian University named after P.A. Stolypin, Omsk
About the modeling of underflooding and land drainage with the fortran
Список литературы О моделировании подтопления и дренирования земель с помощью фортрана
- Прогнозы подтопления и расчет дренажных систем на застраиваемых и застроенных территориях: справ. пособ. к СНиП/ВНИИ ВОДГЕО. -М.: Стройиздат, 1991. -272 с.
- Климентов П.П. Динамика подземных вод/П.П. Климентов, В.М. Кононов. -М.: Высшая школа, 1986. -384 с.
- Немнюгин С.А. Современный Фортран/С.А. Немнюгин, О.Л. Стесик. -СПб.: БХВ-Петербург, 2004. -496 с.
- Развитие исследований по теории фильтрации в СССР/под ред. П.Я. Полубариновой-Кочиной. -М.: Наука, 1969. -546 с.
- Фортран 90. Международный стандарт/пер. с англ. -М.: Финансы и статистика, 1998. -416 с.
- Артемов И.Л. FORTRAN: основы программирования/И.Л. Артемов. -М.: Диалог-МИФИ, 2007. -304 с.
- Алгазин С.Д. Программирование на Visual Fortran/С.Д. Алгазин, В.В. Кондратьев. -М.: Диалог-МИФИ, 2008. -472 c.
- Поршнев С.В. Вычислительная математика/С.В. Поршнев. -СПб.: БХВ-Петербург, 2004. -320 с.