19/11/2021

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

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

17/11/2021

GFLOW бесплатно

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

20/10/2021

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

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

27/09/2021

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

Иногда данных о рельефе территории исследования так мало, что приходится прибегать к помощи космических технологий. Точность у этой информации, прямо скажем, не ахти, особенно на  застроенных или лесистых территориях, но уж что есть —  дарёному коню, как говорится. Знаю, уже есть алгоритмы, которые с помощью нейросетей умеют удалять с этих массивов дома и небольшие лесочки, но пока не встречал их в открытом доступе, пригодном к тому же для использования неквалифицированному пользователю.
Конечной целью предполагается получение CSV файла вида X, Y, Z.  Потребуется следующий инструментарий:
  1. QGIS
  2. Плагины к нему:
    • Какой-нибудь плагин для подгрузки подложек типа QuickMapServices. Ну или загрузить их вручную через XYZ Connections.
    • SRTM-Downloader
    • Point sampling tool
  3. Аккаунты как минимум на одном из сайтов с рельефом со спутников:
Сразу предупрежу, что выудить данные с японского сайта — тот еще квест, но в детали вдаваться не буду. В сети есть инструкции. С 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. Редактировать, разумеется, нельзя. Но можно скачивать и использовать в свое удовольствие.

14/09/2021

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

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

11/08/2021

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

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

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

10/08/2021

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

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