Новая работа

Немного неожиданно, в большой степени формально (впрочем, не без перспектив), но у меня теперь новая строчка в резюме: Ведущий гидрогеолог-эксперт в ФГУП «Геолэкпертиза».

Где искать гидрогеологов

Если вам нужен хороший гидрогеолог и вы готовы принять в штат человека с небольшим опытом, то основным критерием выбора сотрудника после личных качеств становится ВУЗ, который заканчивал соискатель. Я не могу говорить за все университеты, но про столичные ВУЗы у меня скопилась приличная статистика. В Москве более или менее приличных гидрогеологов есть шанс выпустить только у двух заведений: МГУ им М.В.Ломоносова и МГРИ. Из прочих геологических заведений столицы нынче выпускают студентов, верующих в «лозоходство». Буквально.
Если говорить за всю страну, то я слышал о хороших гидрогеологических школах в университетах Санкт-Петербурга и Новосибирска (аберрация памяти — не слышал я ничего про гидрогеологов из Новосибирска, а мой научный руководитель, хоть и заканчивал НГУ, но во-первых он не гидрогеолог, а во-вторых, то был математический факультет), но лично со специалистами недавних выпусков оттуда не знаком, к сожалению.
P.S.: Простите, если кого обидел - ничего личного. Более того, буду рад ошибаться. Кстати, эту заметку вполне можно считать моей рекомендацией - куда стоит идти учиться на гидрогеолога в столице.

How Rivers Work - The Role of Groundwater. Отличное общеобразовательное видео, в котором, что называется, на пальцах показана роль подземных вод в круговороте воды.

This video, produced by the UK Groundwater Forum, helps explain the important role that groundwater plays in the environment. The video shows the contribution groundwater makes to the flow in Britain’s rivers and in doing so provides some basic information on groundwater, subsurface dynamics and the hydrological cycle as a whole.

Как моделировали в 80-х

What did groundwater modeling look like in 1980?
Охренеть.

Гидрогеологический LinkedIn

Если вы владеете языком «вероятного противника» хотя бы на уровне чтения специализированной литературы и у вас есть аккаунт на LinkedIn, то я настоятельно рекомендую вступить в тамошнюю группу, посвященную геофильтрационному моделированию. Там что ни тема — то заголовок для диссертации. Очень сильный состав участников и глубина затрагиваемых вопросов.

Как задавать граничные условия

В комментарии к моей недавней записи о схематизации граничных условий мне совершенно логично заметили, что я написал много слов по теме, но так и не рассказал, как задавать эти граничные условия. Исправляю допущенную, не скрою — умышленно, оплошность.
Граничное условие I-рода в MODFLOW может быть задано двумя способами. Первый способ: в массиве IBOUND (меню Grid/Cell Status) в соответствующие ячейки надо занести любое целое отрицательное значение (не обязательно -1, можно и -23 — разницы никакой), затем в массив Initial & Prescribed Hydraulic Head (в меню Parameters) в этих же ячейках надо задать величину напора в этих ячейках. Тут важно не допустить самую главную ошибку новичка — не задать безнапорном слое I-род с напором ниже отметки подошвы слоя. Второй способ: Time-Variant Specified-Head (MODFLOW/Flow Packages), где можно задать линейно меняющийся напор на границе. Напор, разумеется, может и не меняться. В принципе, почти любая граница III-рода может быть превращена в I-род, для этого достаточно задать очень большую проницаемость границы, но об этом позже.
Границы II-рода задаются так-же двумя способами: пакет Well (MODFLOW/Flow Packages), где отрицательными значениями задается постоянный отток, а положительными — приток; и пакет Recharge (MODFLOW/Flow Packages), с помощью которого задается инфильтрация (положительное значение) или испарение (отрицательное).
Для границ III-рода придумано наибольшее количество пакетов, которые постоянно пополняются по мере выхода новых версий MODFLOW. Все нижеперечисленные пакеты расположены в меню MODFLOW/Flow Packages. General Head Boundary — наиболее общий способ задания г.у. III-рода, там все очевидно — задается напор и проницаемость границы (величина, обратная общепринятому у нас фильтрационному сопротивлению). Drain — так задаются дрены и участки высачивания, принцип аналогичен GHB, но если уровень воды в ячейке с границей оказывается ниже уровня воды, заданного на дрене, то расход через границу становится нулевым (GHB в аналогичном случае станет питающей границей, Drain питающей границей быть не может). River — пакет реки позволяет учесть эффект возникновения режима дождевания при отрыве уровня в пласте от дна реки (для этого в дополнение к уровню на границе и ее проницаемости в этом пакете задается отметка дна реки). Evapotranspiration — пакет для задания испарения с поверхности грунтовых вод, тоже формально может быть отнесен к границам III-рода.

Осушающиеся блоки модели

MODFLOW очень «не любит» маломощные безнапорные горизонты. Когда уровень воды в ячейке модели становится очень близок к отметке ее подошвы, MODFLOW ее делает неактивной и больше в расчете не использует. Даже если вокруг потоп и все окружающие эту ячейку блоки затоплены «под завязку», MODFLOW ее не включит в расчет и не обводнит. Следствием такой особенности являются проблемы с точностью и сходимостью расчета. Если отключенных ячеек оказывается относительно много, то модель почти гарантированно не сойдется при умолчальных настройках расчетных модулей (Solvers).
Для борьбы с этой напастью придумано масса способов:

• Можно поиграться с настройками солверов. Особенно помогает уменьшить Damping Parameter в настройках модуля PCG2.
• Иногда помогает просто сменить солвер. Неплохо сходится солвер GMG, но он поддерживает только MODFLOW-2000 и не все препроцессоры с ним совместимы.
• Полезно воспользоваться пакетом «обводнение» (Wetting Capability), который специально создан для борьбы с этим эффектом. К сожалению, само по себе использование этого пакета может привести к проблемам со сходимостью. За что боролись, как говорится. Хотя иногда он очень помогает, особенно при нестационарных расчетах процессов обводнения изначально сухих горизонтов.
• Самый радикальный способ — отказаться от расчета в безнапорной постановке и считать все слои напорными. В этом случае MODFLOW не будет отключать никакие сухие ячейки, в силу того, что проводимость в них не зависит от уровня воды. Такой подход приводит к заведомо большим расходам потока, но зато позволяет быстро получить хоть какой-нибудь результат. К тому же, при расчете тех же дренажей, небольшое увеличение расхода никакой беды не представляет, создавая некий запас прочности (главное помнить, что таких «запасов» часто оказывается не один и не два и в сумме они нарисуют такой  «запасище», что проектанты на стенку полезут).

О прикладных моментах workflow

Занятный психологический момент заметил тут на днях. Раньше, когда я работал в большой гидрогеологической компании, я сначала проводил все расчеты, рисовал графику, а уж потом начинал писать текст записки. Это несомненно наиболее правильная постановка workflow. Однако, на вольных хлебах, в силу особенностей распределения рабочего времени, мне часто приходится  «плясать от двери» — начинать работу с написания «рыбы» записки. Это не очень удобно и не совсем правильно с методической точки зрения, но так уж сложилось.
И тут возникает забавный момент: поскольку раньше на написании записки работа заканчивалась, то сейчас, дописав «рыбу», я подсознательно переключаюсь из состояния «а-а-а, все попало, ничего не успеваю» в состояние «все нормально, еще куча времени». Хотя объективной причины к этому совершенно не наблюдается, даже наоборот.