Дайджест ссылок

Hydrodynamic Porosity: A New Perspective on Flow through Porous Media, Part I и Hydrodynamic Porosity: A New Perspective on Flow through Porous Media, Part II. Крайне любопытные статьи о моделировании проницаемости поровых сред. Там, как оказалось, не «всё так однозначно».

Pore-scale flow velocity is an essential parameter in determining transport through porous media, but it is often miscalculated. Researchers use a static porosity value to relate volumetric or superficial velocities to pore-scale flow velocities. We know this modeling assumption to be an oversimplification.

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

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

The World Wide HYdrogeological Parameters DAtabase. Матёрая база данных гидрогеологических параметров и свойств пород. Пополняется всем миром. Если кто подскажет, как её открыть без плясок с бубном - буду премного благодарен, т.к. с наскоку никакие вьюверы sqlite баз данных это дело открыть не смогли. Пришлось заново прогонять скрипт mysql2sqlite для конвертации дампа MySQL в sqlite.

This is the landing page of the wwhypda project, a collaborative project aimed at collecting into a unique repository all the measurements of hydrogeological parameters spreaded into many technical reports and pubblications.

Accurate Simulation of Flow through Dipping Aquifers with MODFLOW 6 Using Enhanced Cell Connectivity. Оказывается, MODFLOW6 научился учитывать вертикальное смещение смежных модельных блоков. В классической постановке ячейки и вовсе могут быть разорваны в пространстве (т.е. отметка подошвы одной ячейки выше отметки кровли в соседней) и это никак на влияет на результаты расчетов потока между этими ячейками.

In simulations of groundwater flow through dipping aquifers, layers of model cells are often “deformed” to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells within the same model layer are vertically offset from one another, and the standard conductance-based (two-point) formulation for flow between cells does not rigorously account for these offsets. The XT3D multi-point flow formulation in MODFLOW 6 is designed to account for geometric irregularities in the grid, including vertical offsets, and to provide accurate results for both isotropic and anisotropic groundwater flow

Проницаемость трещины: теоретическая зависимость и ключевые логические следствия. Замечательная статья от уважаемого Александра Лукьянова о трещинной проницаемости скальных пород. Как нельзя по пористости судить о проницаемости (см. первую статью подборки), так и количество трещин само по себе ничего о ней не говорит. Куда теперь бедному крестьянину податься? Не откачки же делать, в самом деле. Это древнее искусство считай забыто.

В моей практике был случай, когда рудник получил весьма неприятный аварийный водоприток с расходом 2500 м3/час (примерно 700 л/с) из одной-единственной горной выработки. Зная перепад напоров (который был весьма существенным), я примерно оценил ширину раскрытия трещины протяженностью 10 м, которая, условно, могла бы опоясывать эту горную выработку по кругу и обеспечивать в неё этот приток. Расчётное раскрытие трещины получилось равным 4,8 мм. Это вовсе не означает, что в данном конкретном случае приток шёл через одну-единственную трещину. Но он мог бы!

Уходит эпоха

Сегодня в Москве в 8 утра скончался Борис Владимирович Боревский.

Дайджест ссылок

Горнопромышленная гидрогеология. Что-то расплодилось окологидрогеологических телеграмм-каналов и сообществ. Вот еще одно. Хорошее.
QGIS для начинающих. Видеоуроки от Василия Швалева. Вот прямо азы, зато бесплатно, а за настоящим «Qgis кунг-фу» — добро пожаловать на его платные курсы. Рекомендую — специалист хороший, берет недорого.

Открытый курс на котором научитесь основам работы в QGIS, даже если только сегодня открыли программу

Getting Started with MODFLOW. Книга под авторством Richard B. Winston, посвященная моделированию в Modflow с помощью ModelMuse. Надо отметить, что благодаря этой книге я почти перестал бояться интерфейса ModelMuse.

Numerical modeling of groundwater flow systems was once accessible only to modeling specialists in the hydrogeological community. Software such as MODFLOW—the most frequently used groundwater modeling program in the world—and associated graphical user interfaces (GUIs) have made modeling possible for most groundwater scientists. This book provides the bridge from understanding to implementing models by introducing the basics of MODFLOW version 6 and providing readers who have a working knowledge of groundwater flow with a guide through construction of their first groundwater model.

Несколько полезных для гидрогеологов библиотек на python: pastas, ttim и timml.

Лекции о гидрогеологии и гидрологии океанских и морских побережий

2023 Darcy Lecture 1: Subseafloor Hydrogeology: Moving Beyond Watersheds

2023 Darcy Lecture 2: Please Check Your Assumptions at the Coastline

Мощности горизонтов

На заметку коллегам, которые, к моему удивлению, умудряются путаться даже в таком, казалось бы, простом вопросе: 

1. Мощность напорного водоносного горизонта считается как разница в отметках (глубинах) залегания кровли и подошвы слоя. Кровля — соответствует низу верхнего водоупора, подошва — верху нижнего водоупора. 
2. Мощность безнапорного водоносного горизонта считается как разница в отметках (глубинах) залегания уровня воды в горизонте и подошвы слоя.

Установка последней версии SagaGIS в QGIS

В профильном сообществе попросили написать гайд по установке актуальной версии SagaGIS в QGIS. Исполняю обещание. В принципе, в составе QGIS SagaGIS уже включён, но во-первых, это временное явление (вроде как в новых версиях решили исключить SagaGIS из ядра), а во-вторых, в составе установочного пакета QGIS поставляется какая-то уж очень древняя версия Saga — кажется 7-я, тогда как актуальная версия на момент написания заметки — 9.2. 

По пунктам, кратенько. Потом постараюсь добавить подробностей, когда будет больше свободного времени. 

1. Скачиваем последнюю версию SagaGIS. 
2. Распаковываем скачанный архив в какую-нибудь папку. Я предпочитаю OSGeo4W64\apps\, но можно любую другую.
3. В QGIS отключаем плагин ядра SAGA GIS provider и скачиваем Processing Saga NextGen Provider. В Настройках последнего указываем путь к SagaGIS — куда мы его распаковали на предыдущем шаге. Закрываем QGIS. 
4. В папке с SagaGIS лежит архив saga4qgis.zip, внутри которого есть инструкция, что надо делать с содержимым этого архива. Не буду повторять её, но суть такова: содержимое архива надо распаковать (удалив предварительно уже существующее) в папку     C:\Users\...\AppData\Roaming\QGIS\QGIS3\profiles\default\python\plugins\processing_saga_nextgen 
   Только не делайте это «в лоб», сперва найдите где конкретно лежит файл SagaNameDecorator.py и папка \description\ — удалите их и скопируйте на их место содержимое архива.
5. В принципе всё готово, можно запускать QGIS и радоваться возможности работать с последней версией SagaGIS. Но есть еще одна маленькая хитрость, сильно облегчающая работу. Надо открыть файл SagaAlgorithm.py и заменить

elif isinstance(param, (QgsProcessingParameterString, QgsProcessingParameterField)):
                command += ' -{} "{}"'.format(param.name(), self.parameterAsString(parameters, param.name(), context))

на

elif isinstance(param, (QgsProcessingParameterString, QgsProcessingParameterField)):
                command += ' -{} "{:.10}"'.format(param.name(), self.parameterAsString(parameters, param.name(), context))

Для чего нужен последний пункт: дело в том, что инструменты SagaGIS, работающие с векторными исходниками, на входе принимают только shape файлы SHP. А те, в свою очередь, согласно стандарту не поддерживают имена полей длиннее 10 символов. Поэтому если вы попробуете подсунуть в SagaGIS на вход какой-нибудь Geopackage с именем поля 11 символов и больше, то в исходном варианте SagaAlgorithm.py вы получите ошибку. Мой же код автоматически обрезает имя поля по 10 символу и в результате формируется правильный shape-файл. 
К слову сказать, похожие улучшения можно вставлять и в другие Processing Provider плагины — они все сделаны более или менее похожим образом. Однако важно понимать, что если вы ухитрились создать в Geopackage несколько столбцов (полей) длинными именами, первые 10 символов которых побуквенно совпадают, то готовьтесь к удивительным глюкам.

Фриланс и гидрогеология

По мотивам увлекательной беседы в одном профильном сообществе в Telegram родилась идея поделиться советами и мыслями с начинающими фрилансерами на ниве гидрогеологии. Если кто не знает, то сам я «фрилансю» уже более 10 лет, так что опыт какой-никакой, да имеется.

1. Переходить на 100% фриланс имеет смысл только тогда, когда доход от «леваков и шабашек» начнет существенно превышать зарплату. Причем не в месяц, а в среднем за год, а лучше полтора-два.
2. Очень сложно не психануть в периоды кризисов, когда по несколько месяцев вообще нет работы. Соответственно надо иметь подушку безопасности минимум на полгода привычной жизни.
   
3. Из п. 2 следует несколько выводов:
• Перед уходом в свободное плавание эту подушку следует сперва наполнить. Т.е. сначала подушка, потом увольнение из офиса.
• Если у вас ипотека или кредиты, то вероятно оптимальным вариантов будет сперва досрочно закрыть их за счет дохода с «леваков», а уж потом увольняться. Однако это плохо сочетается с актуальным рынком ипотеки, которую внезапно стало невыгодно досрочно погашать (но это отдельный большой вопрос, в котором я не сильно разбираюсь, просто краем уха слышал).
• Хочешь-не-хочешь, но придется вникать в тему инвестиций. С нынешним уровнем инфляции хранить деньги на депозите — равно терять по 2-3% накоплений в год.
4. Крайне рекомендую перед уходом в свободное плавание обзавестись как минимум одним, а лучше 2-3 постоянными клиентами с регулярными заказами. Будьте готовы к тому, что для этого костяка либо придется надолго зафиксировать цены, либо изначально обговорить постепенное повышение тарифа на уровень инфляции.
5. Поток заказов имеет довольно чёткую сезонность с двумя пиками (догадайтесь с трёх раз, с чем они связаны): конец весны и конец года. Соответствующим образом следует подходить к тайм-менеджменту: к началу мая и ноября крайне желательно подходить с минимумом текучки, иначе новые заказы не переварите.
6. Не стесняйтесь отказываться от работы. Чувствуете, что не сдюжите — скажите, чтоб писали мне. Я постараюсь что-нибудь придумать. Не хапайте всё подряд, особенно если чувствуете, что заказчик «мутный». Да, таких в последнее время тоже полно. И ладно бы просто не платили. Бывает и гораздо хуже — можно за 100 рублей и «срок с пола поднять».
7. Не беритесь за срочные заказы в области, в которой плохо разбираетесь, но если чувствуете  себе желание расти и развиваться, то при гуманных сроках не стоит огульно отказываться от всего нового. 
8. Держите связь с коллегами и если уж отказываете заказчику, то обязательно посоветуйте пару человек, к которым можно обратиться.

О новостройках

Недавно довелось принять участие в комплексных инженерно-геологических изысканиях под строительство весьма пафосного жилого комплекса. Ну очень пафосного, в очень недешевом месте. При этом наблюдалось довольно редкое явление: к изыскателям и проектировщикам прикрепили куратора из весьма уважаемого института. Так вот по совершенно непонятной мне причине, этот самый куратор с усердием достойным лучшего применения занимался кроиловом буквально на карандашах. И не только на них: «насыпухи много набурили — убирайте, а то мы на сваях разоримся», «вода слишком солёная в пробах — поправьте, а то вода слишком к арматуре агрессивна» (а там реально рассол из-за того, что рядом находится склад поваренной соли), ну и вишенка уже по моей части: «что-то у нас зона влияния от водопонижения слишком большая — сделайте коэффициент фильтрации в 2 раза меньше и пересчитайте прогноз, а то нам придется половину района на аварийность обследовать». Причем редкий случай, когда изыскатели таки соизволили сделать реальные откачки. Не очень представительные, но всё-же. Похожие запросы у него и к проектировщикам, но тут я не специалист — просто по тону их ответов во время совместной телеконференции понял, что требуют от них подобной залипухи.
В общем, если уж в уважаемых институтах такой бардак и развал, то всё, приехали. Меня иногда спрашивают, каким застройщикам можно доверять при покупке недвижимости в столице — по всей видимости, никаким. Беда в том, что с такими строителями нельзя успокаиваться и тем, кто живёт в старом жилом фонде — вот построят такие вот молодцы свою башню у тебя во дворе по программе реновации, а у тебя по дому трещины из-за неравномерных осадок.