О платных консультациях

Есть такой анекдот:

- Саша, хотел у тебя спросить. Задача вот какая...
- Погоди, погоди, тебе совет или консультацию?
- А в чем разница?
- Совет бесплатный, консультация за деньги.
- Совет, конечно!
- Мой тебе совет: запишись на консультацию.

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

GFLOW бесплатно

Еще одна крайне приятная новость: с недавних пор создатель очень хорошей программы GFLOW, предназначенной для фильтрационного моделирования с помощью методики аналитических элементов, Dr. Henk Haitjema сделал свою программу бесплатной. Однако техническая поддержка оказывается только для купивших программу пользователей, имейте это в виду.
Кстати, насчет аналитических элементов: есть еще одна бесплатная программа, поддерживающая сразу несколько АЭ «солверов» — Visual AEM. Жалко только, что она не поддерживает современные версии солвера TimML — там ведь даже трехмерные задачи можно решать.

КиберЛенинка

Не прошло и 10 лет, как индеец Зоркий Глаз обнаружил, что на сайте КиберЛенинки доступно немало статей по гидрогеологии на русском азыке.

Как выудить в точки рельеф из SRTM или AW3D30

Иногда данных о рельефе территории исследования так мало, что приходится прибегать к помощи космических технологий. Точность у этой информации, прямо скажем, не ахти, особенно на  застроенных или лесистых территориях, но уж что есть —  дарёному коню, как говорится. Знаю, уже есть алгоритмы, которые с помощью нейросетей умеют удалять с этих массивов дома и небольшие лесочки, но пока не встречал их в открытом доступе, пригодном к тому же для использования неквалифицированному пользователю.

Конечной целью предполагается получение CSV файла вида X, Y, Z.  Потребуется следующий инструментарий:

1. QGIS
2. Плагины к нему:
• Какой-нибудь плагин для подгрузки подложек типа QuickMapServices. Ну или загрузить их вручную через XYZ Connections.
• SRTM-Downloader
• Point sampling tool
3. Аккаунты как минимум на одном из сайтов с рельефом со спутников:
• https://search.earthdata.nasa.gov/ - регистрация по ссылке https://urs.earthdata.nasa.gov//users/new
  
• https://www.eorc.jaxa.jp/ALOS/a/en/dataset/aw3d30/aw3d30_e.htm - регистрация по ссылке https://www.eorc.jaxa.jp/ALOS/a/en/dataset/aw3d30/registration.htm

Сразу предупрежу, что выудить данные с японского сайта — тот еще квест, но в детали вдаваться не буду. В сети есть инструкции. С NASA возни гораздо меньше, особенно если вам повезло и SRTM-Downloader работает как надо.

Итак, алгоритм следующий:

1. В QGIS любым известным способом открываем картографическую подложку. Это не обязательно, но как правило сильно облегчает работу.
2. Создаем временный слой типа Polygon. Рисуем в нем полигон по контуру территории, для которой нам нужен рельеф. Зуммируемся так, чтоб нарисованный полигон занимал большую часть экрана (Zoom to Layer).
3. С помощью строенного в QGIS инструмента Vector/Random points in polygons генерируем достаточное количество точек в пределах нарисованного в предыдущем пункте полигона. При этом создается новый временный слой с точками.
4. С помощью SRTM-Downloader скачиваем растры с рельефом для области, показанной на экране. Для чего жмём последовательно Set Canvas Extend и Download. Для данных ALOS придется скачивать данные через сайт и подгружать в QGIS вручную.
5. Запускаем Point sampling tool, указываем слой с точками из пункта 3 в качестве Layer containing sampling points, а в качестве Layers with fields/bands to get values from — растры (можно указать сразу несколько) из пункта 4. Тут важно отметить, что точки и растры должны быть в одной системе координат (WGS 84 -  EPSG 4326 в случае если мы работаем с SRTM). 
6. На выходе получаем векторный слой с точками с атрибутами в виде значений рельефа. Если область интереса попадает сразу на несколько растров, то столбцов с атрибутами будет несколько. Надеюсь, не надо рассказывать, как из нескольких столбцов получить один — это можно сделать как через Field Calculator в самом QGIS, так и в Excel или любом другом табличном редакторе.
7. Пересохраняем полученный точечный слой в виде CSV, не забывая выбрать нужную нам систему координат (если не выбирать, то по умолчанию  координаты будут представлены в виде градусов) и указать, что GEOMETRY сохраняется в виде AS_XY.

Аналитический расчет барражного эффекта

Давно чесались руки реализовать в электронных таблицах расчет, представленный в монографии В.И. Сологаева «Фильтрационные расчеты и компьютерное моделирование при защите от подтопления в городском строительстве» и реализованный в Ansdimat.
Наконец-то появилось на это время. Расчет доступен на Google Drive. Редактировать, разумеется, нельзя. Но можно скачивать и использовать в свое удовольствие.

Processing Modflow 8 бесплатно

К слову сказать, с июля этого года компания Simcore начала бесплатно раздавать Processing Modflow 8. В дополнение к уже давно бесплатной Processing Modflow 7. А версию 5.3 так и вовсе найти почти невозможно.

Шаблон обработки кустовой откачки в Excel

Вообще автоматически обрабатывать откачки в Excel довольно проблематично. В моем шаблоне сделана довольно кривая попытка, но она нормально работает только на рисованных откачках и, очевидно, сбоит на реальных. В качестве исходных данных там использована как-раз таки «рисованная» по такому случаю в Modflow откачка.

Скачать его можно по ссылке: Gdrive.

Processing MODFLOW — почему не считается модель

Опытные коллеги конечно же об этом знают, а вот начинающим будет полезно знать, что найти причину ошибки можно двумя путями:

1. При запуске модели на счет (Models/MODLFOW (Run Simulation)/Run...) надо поставить галочку у пункта Check Model Data. Запустить счет, а потом в простом текстовом редакторе посмотреть содержимое файла CHECK.LIS, который появится в каталоге с моделью.
2. Некоторые неочевидные ошибки проверка из предыдущего пункта к сожалению пропускает. Тогда придется после неудачного запуска модели на счет воспользоваться «тяжелой артиллерией» — открыть файл output.dat (опять же в простом текстовом редакторе). В нем хранится вся служебная информация о последнем запуске, в т.ч. и ошибки. Файл пишется последовательно по мере выполнения MODFLOW, поэтому указание на ошибку чаще всего надо искать в самом конце файла. 

Кстати, похожие инструменты доступны и в других препроцессорах.