О разделении труда

На днях пришло письмо от проектировщика ВК (водоснабжение и канализация) с просьбой помочь в расчете водопритока к дренажу с «верховодки». У человека явно серьезные проблемы с понятийным аппаратом (не потому что дурак, а потому что не учат у нас проектировщиков гидрогеологии, и правильно делают), но надо формулу найти и надо как-то ее применить. Нашли формулу, но из формулы стало ясно, что ничего не ясно — водоприток то с «верховодки» целиком и полностью определяется формой и размерами линзы водоупора, на котором образовалась эта «верховодка». Отсюда возникает необходимость в дополнительных инженерно-геологических изысканиях. В общем, все грустно.
И вот так оно везде. Для меня навсегда останется загадкой: почему проектировщиков регулярно заставляют заниматься гидрогеологическими расчетами, тогда как гидрогеологов довольно редко просят что-то проектировать. Хотя и такое случается, особенно в области лицензирования месторождений подземных вод. С довольно закономерным результатом, надо заметить. Уверен, что любой настоящий проектировщик будет со слезами от смеха смотреть на все эти «писульки», что называются у нас «проектом» водозабора или разработки месторождения.

Обновления FAQ

Вопрос #14:
Для моделирования работы деталей элементов мелиоративных систем мне понадобится существенно сгущать сетку и тут есть ограничения на на размеры соседних ячеек. В текстах про Visual MODFLOW вычитал, что не более чем 1,5 раза и соотношение сторон ячейки не более чем 1:10. В некоторых примерах сетки видел нарушения этих ограничений, может не такое строгое оно?

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

Вопрос #15:
Где можно раздобыть примеры решения задач на Modflow для самостоятельного обучения моделированию?

Ответ:
В составе программного комплекса Processing Modflow (где скачать — см. выше) есть примеры и пошаговое руководство по созданию модели. Что насчет моих собственных примеров, то тут все довольно сложно. Специального набора примеров у меня нет, а «рабочие» модели либо слишком сложны, либо являются коммерческой тайной моих заказчиков, а часто и то и другое.

Граничное условие I-рода

Что такое граничное условие I-рода применительно к геофильтрационной схематизации — это граница с заданной на ней функцией изменения напора во времени H(t). В частном случае, когда напор постоянен и неизменен H(t)=const (если расчет ведется «в понижениях», то на такой границе понижение dH(t) = 0).
В природе границ I-рода не бывает. Вообще никогда. II-рода, кстати, тоже. Есть только III-род (т.е. граница с заданной на ней зависимостью расхода от напора Q=(Hг-H)/С), а все остальные граничные условия — это лишь частные случаи и упрощения. Ноги у этих упрощений растут из тех времен, когда расчеты велись с помощью аналитических формул. Конечно, третий род тоже считался аналитически, но вспомните как усложняются расчеты, когда возникает сопротивление на границе, а если оно нелинейно, то тогда вообще «тушите свет». Для академических упражнений в математике это может и семечки, но когда дело доходит до практики — не до изысков.
С появлением численных методов дело существенно упростилось, но все-равно исследователи продолжали по-возможности упрощать третий род до второго или первого, поскольку это позволяло иногда довольно существенно ускорить расчет.
Современный инструментарий гидрогеолога-модельера позволяет вовсе не использовать границы первого рода (со вторым чуть сложнее) — General Head Boundary с высоким значением проницаемости (Conductance) в принципе ничем не отличается от границы с заданным напором, а в случае чего — эту Conductance можно и скорректировать (только не забудьте это обосновать).

Гидрогеологический дайджест

Depletion of Central Valley's groundwater may be causing earthquakes. Снижение уровня подземных вод может служить причиной землетрясений. И наоборот, кстати.

For years, scientists have wondered about the forces that keep pushing up California's mighty Sierra Nevada and Coast Ranges, causing an increase in the number of earthquakes in one part of Central California.

Можно ли использовать воду из скважины в коттеджном поселке для питья? Добротная порция «страшилок» от уважаемого коллеги из Питера. Жаль, что в статье не приведена статистика о том, сколько коттеджных поселков водоснабжаются из нелегальных скважин. Подозреваю, что эта цифирь должна вселять железобетонную уверенность в завтрашнем дне у специалистов, занимающихся оценкой запасов подземных вод — без работы они точно не останутся. Разумеется, если у нас вдруг законодательство не изменится.

В соответствии с действующим законодательством, эксплуатирующая организация должна самостоятельно подавать отчетность об объемах добычи подземных вод, положении уровней в питьевом водоносном горизонте и качестве подземных вод в контролирующие органы. Однако, если управляющая организация не получила лицензию на право добычи подземных вод, надзорные органы могут не подозревать о существовании водозаборной скважины в том или ином коттеджном поселке.
То есть, подобного рода скважины используются нелегально, без лицензии, что при проведении плановой проверки сотрудниками Роспотребнадзора или Росприроднадзора может повлечь за собой штраф до 1 миллиона рублей (статья 7.3 КОАП РФ) или полный запрет на эксплуатацию скважины и требование об ее ликвидации.

MOD16 for actual evapotranspiration measurement. Карта величины эвапотранспирации, полученных со спутниковых датчиков. Карты лежат на ftp-сервере в форматах GeoTIFF и hdf (т.е. легко подгружается в любой более или менее современных GIS-пакет).

This project is part of NASA/EOS project to estimate global terrestrial evapotranspiration from earth land surface by using satellite remote sensing data. MOD16 global evapotranspiration product can be used to calculate regional water and energy balance, soil water status; hence, it provides key information for water resource management.

Australian groundwater modelling guidelines. Кратенькое (всего-то 200 страниц) руководство по гидрогеологическому моделированию из далекой Австралии. Очень толковое, должен заметить: мне особенно понравились советы по построению диаграммы сопоставления наблюденных и расчетных значений — узнал несколько полезных приемов.

Groundwater models are computational methods that simulate the profile of an underground water system. They have proven to be useful tools to provide insight into complex water system behaviour; to address a range of groundwater problems; and to support decision-making processes.

Решение обратных задач с помощью Visual Modflow

На днях со мной связался автор одной очень любопытной программы, предназначенной для автоматизации решения обратной задачи геофильтрации с помощью Visual Modflow. С любезного согласия автора я выложил саму программу и её описание у себя на соответствующей странице (она всегда доступна в меню справа).

От себя добавлю, что лично мне кажется такой подход к решению ОЗ не отличается «физичностью» результата. Да, в результате получается очень похожая на фактическую модельная карта гидроизогипс, но это крайне редко бывает самоцелью. Считать какие-либо прогнозы на такой модели следует крайне осторожно. Тем не менее, в качестве первого приближения к решению обратной задачи такой подход имеет все права на существование.

Гидрогеологический дайджест

Прошу прощения за долгое молчание. Дел навалилось — не почать и край.

Уникальные геологические объекты России. Интересный ресурс для любителей побродить по геологическим достопримечательностям родного края.

Здесь представлена информация о геологических достопримечательностях России: сейчас в системе более 2600 уникальных геологических объектов (геологических памятников природы) и около 3800 стратиграфических объектов, т.е. стратотипов (типовых разрезов) и опорных стратиграфических разрезов.

Хочу представить интересный гидрогеологический блог: Gidahatari. Насколько я понял, это корпоративный блог перуанской фирмы, занимающейся гидрогеологическими работами. Вот, примеры последних записей:

• How to do a groundwater flow model with your company logo. Как сделать геофильтрационную модель с логотипом вашей фирмы. Не бином Ньютона, но результат вполне симпатичный — сойдет для украшения презентации.
• Feflow or Modflow? A recurrent question. А вот это уже более серьезная тема. Кстати, а в нашей стране есть гидрогеологи, активно использующие Feflow?

О том, почему у нас аквапарки разваливаются

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

P.S.: Справедливости ради стоит отметить, что по результатам состоявшегося совещания и написанных гневных писем (не только и не столько моих), строители таки сделали шпунтовую завесу и нормальное водопонижение.

Получение карты гидроизогипс по архивным данным инженерных изысканий

В ряду бредовых «хотелок» некоторых экспертов МГЭ достойное место занимает карта гидроизогипс, построенная по архивным скважинам. Желают они её видеть, не отдавая себе отчет в том, что смысла такая карта несет ровно 0 целых, 0 десятых. И вот почему. Во-первых, не секрет, что инженеры-геологи не особенно стараются при гидрогеологическом обследовании скважин — пошел мокрый песочек на шнеке, значит типа водоносный горизонт. Точность такого определения — плюс-минус полтора метра, если не хуже. А во-вторых, архивные данные обладают таким неприятным свойством, как разновременностью замеров. Какие-то изыскания делались в период зимней межени, а какие-то в период весеннего паводка. Все это приводит к тому, что попытка нарисовать по таким точкам гидроизогипсы или гидроизопьезы — бестолковое, а главное — методологически вредное занятие.
Особенно веселые картинки получаются, если не мудрствовать и запихнуть всё в какой-нибудь автоматический интерполятор типа Golden Software Surfer. Мало кто утруждает себя глубокими настройками параметров интерполяции, а в результате получаются «гидроизогипсы», говорящие о наличии источников/стоков в самих наблюдательных скважинах, что в свою очередь является свидетельством вопиющей безграмотности специалиста (но экспертов такие картинки радуют, что характерно).
Справедливости ради стоит заметить, что есть способ сделать всё методологически чисто: нужно снормировать архивные показатели с учетом колебаний уровня в ближайшей режимной скважине. Довольно кропотливая работа, даже при наличии этой самой режимной скважины, в которой к тому же имеются замеры на весь период, охваченный архивными данными. В противном случае, когда архивы у вас имеются, к примеру, с 1956 по 2002 гг, а режим есть только с 1980 по 1993, то задачка становится еще более нетривиальной. Анализ среднегодовых и среднемноголетних колебаний — это, поверьте, то еще развлечение. И вот только по таким нормированным показаниям можно строить вожделенную карту гидроизогипс. Смысла от которой все-равно немного — разве что показать плохому специалисту, куда и откуда вода течет, но ему она все-равно не поможет, а хорошему не нужна.
Однако, существует один частный случай, когда такая карта имеет хоть какой-то смысл — это когда одна или несколько внешних границ модели задаются I-родом с постоянным напором, полученным на основе интерполяции уровней воды в этих самых архивных скважинах. Но об этом грязном приемчике как-нибудь в другой раз.
И да, забыл уточнить: все вышесказанное имеет отношение к задачкам типа расчета барражного эффекта при освоении подземного пространства и прочих т.н. инженерно-гидрогеологических расчетах. Допустим, при поиске и разведке подземных вод такая карта строится обязательно и разумеется с применением множества методологических приемов и хитростей.