Таблички для обработки и анализа откачек

Spreadsheets for the Analysis of Aquifer-Test and Slug-Test Data. И как я их раньше не замечал. Настоятельно рекомендую почитать сперва FAQ — сразу отвалится куча глупых вопросов.

И еще раз об осушенных блоках модели

Возвращаясь к напечатанному. Довольно типичная ситуация: «стена в грунте», внутри стены задан скважинный дренаж. Откачивающие скважины заданы как граничное условие второго рода. Запускаем расчет и получаем вот такую картину:

Такое впечатление, что из шести скважин работают только верхние две. Самое забавное, что так оно и есть, но в чем же дело?
Вот как вся эта ситуация это выглядит при итеративном расчете:

1. Работает скважина, снижается уровень воды.
2. Вдруг уровень стал ниже подошвы слоя — отключили ячейку со скважиной.
3. О, скважина исчезла (ячейка то с ней отключена) — уровень воды начинается повышаться.

А проблема в том, что отключенная из-за осушения ячейка не включится, а следовательно и скважина не заработает.
Как бороться с такой напастью. Да очень просто, почти все рецепты я уже описал в своей предыдущей записи на эту тему, к которым можно добавить еще один, а именно: задать скважины не расходом, а напором (понижением) — кстати, вовсе не обязательно использовать для этого первый род, можно попытаться учесть сопротивление прискважинной зоны и задать скважины третьим родом (пакетами General Head или даже лучше Drain). А расход скважин можно получить уже из Water Budget.

Новая работа

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

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

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