Несколько полезняшек для гидрогеолога

Carlos Molano опубликовал на своем сайте несколько полезных таблиц с гидрогеологическими расчетами:

1D analytical transient flow in confined aquifers
Automatic fit for pump test data to the Theis equation
Capture type curves for wells Hantush leaky aquifer test analysis
Inverse particle tracking
Neuman unconfined aquifer test analysis
Piper diagram
Superposition Thiem well field
The Theis test for confined aquifers
The Theis fit for variable pumping rates
The Theis-Jacob Test

Оценка и переоценка запасов подземных вод

Намедни встречался с однокурсниками с кафедры гидрогеологии МГУ. Был приятно удивлён тем, что многие по-прежнему работают по-специальности. В частности, оказалось, что у меня есть кому передавать работы по оценке запасов (сам я этим, к сожалению, не занимаюсь уже занимаюсь, хотя и не беру на себя весь объем работ) — в последнее время участились звонки от заказчиков на эту тему. Кризис, видать, прижал и переплата «крупняку» даже каких-то 200 т.р. — уже серьезный повод обратиться к частникам.
Надо чаще встречаться, да.

Немного о приятном

Не знаю как у тебя, мой читатель, а у меня самый приятный момент в работе — это промежуток времени между нажатием кнопки send в почтовом клиенте и сообщением «письмо отправлено». В письме — оформленный отчет о проделанной научной работе. Именно в этот момент происходит смена состояния «я должен» в положение «мне должны».
Я даже момент получения гонорара ценю не на столько.
Это я к чему: закончился изматывающий почти двухмесячный рабочий марафон, а следовательно у меня освободится немного времени на ведение блога. Следите за обновлениями.

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

How to write a plan for a numerical groundwater model

Хороший план, жаль только, что на английском. Надо будет перевести.

Шикарный подход

«Анализ данных по физико-механическим свойствам горных пород месторождения …, рекомендованных для геомеханических расчетов институтом … показал, что значения угла внутреннего трения и особенно сцепления являются завышенными, соответственно при расчетах устойчивости дадут необоснованно высокие значения углов откосов карьера. Скорректированные в соответствии с методиками … показатели приведены в таблице».

Показатели нам не нравятся, поэтому возьмем те, которые нам хочется. Я уже давно не удивляюсь, что у нас самолеты падают. Для чего тогда вообще изыскания делать, если в расчетах все равно литературные данные используются?

Новый способ решения обратной задачи

A new way to estimate hydraulic conductivity.
Для тех, у кого проблемы с языком вероятного противника, в двух словах. Суть метода заключается в том, что на модель накладывается дополнительная элементная сетка, триангулируется, значение коэффициента фильтрации для каждого треугольника усредняется по величинам в вершинах. Итеративно меняем величину Кф в вершинах до тех пор, пока эти изменения не перестанут влиять на невязку.
Изящно, но лично у меня есть большие сомнения в сходимости.

Processing MODFLOW. Создание модельной сетки.

Итак, наконец-то у меня дошли руки до давно обещанного краткого (подчеркиваю!) руководства по геофильтрационному моделированию в программе Processing Modflow. Сразу предупрежу: у меня нет под рукой бесплатной версии этой программы (5.3.1), поэтому пункты меню я буду приводить на примере седьмой версии. Принципиальной разницы между ними нет, но стоит иметь в виду, когда будете сталкиваться с несоответствиями.

Начинается процесс с создания модели: File/New Model – настоятельно рекомендую создавать для каждой модели отдельную папку и назначать им (модели и папке) – «говорящие» названия.

Затем мы переходим к заданию сетки: Grid/Mesh Size. В открывшемся диалоговом окне предоставляется возможность задать следующие параметры:

• количество слоев (Number of Layers) – тут главное не переборщить, т.к. добавить в случае чего новый слой достаточно легко, а вот удалить практически невозможно (правда, можно отключить, сделав все ячейки слоя неактивными);
• толщину модели (Model Thickness) и отметку кровли верхнего слоя (Model Top Elevation) – опциональные параметры т.к. толщина модели все-равно определяется суммарной мощностью слоев, геометрия которых задается позже, где задается и рельеф кровли верхнего слоя;
• количество ячеек по оси Y или количество строк (Number of Rows);
• размеры модели по оси Y (Model Extent) – очевидно, что размеры ячеек по оси Y будут определяться как отношение Model Extent к Number of Rows;
• количество ячеек по оси X или количество столбцов (Number of Columns);
• размеры модели по оси X (Model Extent);
• вертикальный масштаб (Vertical Exaggeration).

Переключившись во вкладку Coordinate System, мы можем задать координаты сетки и просмотрового окна. Важно заметить, что в качестве начала координат используется верхний левый угол модели. Кроме того, там же задается угол поворота сетки.

После нажатия на “Ok” мы попадаем в главное рабочее окно программы (мы его еще ни раз увидим, когда будем задавать параметры и граничные условия). Здесь мы можем подправить небольшие огрехи, допущенные на этапе задания сетки: подвинуть сетку, изменить угол поворота, изменить размеры ячеек. Для разбиения сетки в областях наибольшего интереса рекомендую сначала выйти из главного окна с сохранением, а потом зайти обратно в Grid/Mesh Size – в этом случае в меню по правому клику на ячейке появятся пункты, отвечающие за разбиение ячейки.

Кстати, уже на этом этапе появляется возможность подгрузить векторные и растровые подложки (Options/Maps…). PM достаточно привередливая программа в плане поддерживаемых форматов: векторный DXF, к примеру, она «кушает» только допотопных версий 2000 года, а растровый файл должен заведомо превышать размеры просмотрового окна. В общем, тут им много чего надо доработать.

Злобное

На протяжении моей без малого (полгода до юбилея осталось, кстати) десятилетней карьеры гидрогеолога, я бесчисленное количество раз отвечал на разной степени идиотизма вопросы экспертов и заказчиков. Сейчас я сам выступаю в качестве въедливого заказчика и пишу злобные замечания на отчеты приснопамятного Гидэка. Получаю ни с чем не сравнимое удовольствие. Вот такая я какашка, как говорилось в одном известном анекдоте.