Showing posts with label полезняшки. Show all posts
Showing posts with label полезняшки. Show all posts

01/12/2023

Установка последней версии 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 символов которых побуквенно совпадают, то готовьтесь к удивительным глюкам.

04/01/2023

Генерация пользовательской СК QGIS

В рамках семинара для коллег по основам работы в QGIS «родилась» вот такая утилитка для генерации пользовательской СК в формате строки "proj.4" на основе уже существующей СК и высчитанных сдвижек по осям X и Y.
Задача традиционная: нам прислали какой-то там DXF в «левых» метровых координатах (далее пусть будет СК-X) с нулём координат у местной водокачки и надо худо-бедно положить его на спутниковую  (или любую другую в общеупотребимой СК) подложку в QGIS. Сразу хочу предупредить, что если вы не можете похвастаться отменным везением, то результат будет именно что «худо и бедно», но хоть как-то.
Простейший вариант решения проблемы: сместить карту прямо в автокаде так, что б общие точки на на карте и подложке более или менее совпадали. Вариант хороший и рабочий, но если вам начнут присылать другие данные в этой «левой» СК-X, то придется и их смещать вручную.
Я предлагаю подойти к решению этой проблемы с другого бока: взять одну из общеупотребимых  спроецированных метровых СК (далее СК-А) и на её основе сгенерировать новую пользовательскую СК (СК-Б - иными словами, это та же СК-X, но с известными параметрами), подобрав для неё смещения по X и Y так, чтоб DXF с присвоенной СК-Б совмещался с подложкой. Тут самое сложное: найти такую СК-А, чтоб угол совпадал с СК-X. По опыту, в 90% случаев подходит либо одна из 3-градусных зон Pulkovo 42, либо соответствующая региональная МСК (частенько эта таинственная СК-X и вовсе является какой-нибудь «МСК-59 зона 2»). Ну а дальше всё просто: «скалываем» координаты общей точки в СК-X (в автокаде, к примеру) и в СК-А (для этого можно использовать один из множества плагинов в QGIS, не забыв сделать СК-А системой координат проекта), копируем их в соответствующие ячейки моей таблицы, потом вставляем строку proj.4 от СК-А и строкой ниже получаем новую строку proj.4 уже для нашей новой кастомной СК-Б. Теперь эту строку можно добавить в качестве пользовательской проекции к QGIS и открывать в ней наш DXF, каталоги скважин и прочие данные.
Отлично понимаю, что с точки зрения тру топографов, географов и картографов все эти ужимки - жуткая ересь, прошу отнестись снисходительно.
P.S.: Добавил усреднение по нескольким общим точкам с проверкой на несовпадение угла поворота.

27/09/2021

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

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

11/08/2021

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

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

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

28/09/2020

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

The Simple Online Groundwater Model. Простая Онлайновая Геофильтрационная Модель. В общем-то и добавить нечего. Для простых прикидочных расчетов очень даже удобно.
The Simple Online Groundwater Model (SOGWaM) is built for budding hydrogeology enthusiasts. It allows you to simulate and visualize groundwater flow, and to understand how system changes affect the flow field by developing and testing hypotheses.
Cadmapper. Онлайновый транслятор открытых данных OpenStreetMap, NASA и USGS в формат Autocad.
Architects and urban planners use Cadmapper to save hours of routine drawing. It transforms data from public sources such as OpenStreetMap, NASA, and USGS into neatly organized CAD files.
FloPyArcade. На базе MODFLOW уже игры пишут. Дожили до светлого дня.
FloPyArcade is a MODFLOW-powered groundwater arcade-type game. It builds on functionality of the library FloPy, which is a wrapper adding pre- and postprocessing options around MODFLOW as well as its related software.
Groundwater Learning Free of Charge. Очень полезное начинание от канадских коллег. Фактически, это онлайн-библиотека с бесплатными книгами и учебными пособиями по гидрогеологии. Насколько я знаю, к проекту уже присоединились русскоязычные авторы и переводчики. Искренне желаю им всяческих успехов.
The Groundwater Project, a non-profit organization, registered in Canada, is committed to contribute to advancement in education by creating and making available online free high-quality groundwater educational material for all. The Groundwater Project is led by Dr. John Cherry, recipient of the 2020 Stockholm Water Prize, and managed by a 11-member international Steering Committee and a larger and globally diverse Advisory Committee.
HatariChem. Онлайн утилита для рисования диаграмм химического состава воды.
We present our own webapp for the representation of the Piper Diagram, Stiff Diagram and Scholler Diagram. The webapp was developed in Python Django and it is entirely free for everyone. The main objective behind this webapp was to develop a user friendly and minimum requirement tool to create these water quality / hidrogeochemical diagrams.
Graphical Construction of Groundwater Flow Nets. Когда я учился в МГУ, нас тоже учили рисовать сетку потока подземных вод. Не думал, что этой теме можно аж целую книгу посвятить.
A groundwater flow net consists of two families of intersecting lines: equipotential lines, which connect locations of equal hydraulic head, and flow lines, which show paths of groundwater flow. Together, these two sets of lines provide a visual, two-dimensional representation of the groundwater conditions under steady state (that is, neither set of lines change with time).
Managing Interference Between Geothermal Installations. Расчет эффекта интерференции между геотермальными системами отопления/охлаждения.
There is an increasing interest in utilizing shallow ground and groundwater as a source for geothermal heating and cooling. Either open- or closed-loop systems can be used for heat exchange with the subsurface to supply heat pumps of buildings.
Environmental Modeling Solutions - Webinars. Серия вебинаров по актуальным гидрогеологическим темам.
How agriculture software for nutrient optimization can help to protect groundwater
According to a new report by the European Environment Agency (EEA), only 40% of the continent’s water bodies are compliant with the agency’s environmental standards set by the European directive on water, which was adopted in 2000. The directive offers a framework for the evaluation, management and improvement of water resources throughout Europe.

05/06/2020

Прореживатель данных с логгеров

В последнее время при осуществлении опытно-фильтрационных работ для замера уровней (или еще каких параметров — минерализации, к примеру) в скважинах все чаще используются т.н. логгеры. Производят их много кто, вроде бы даже есть отечественные производители, но у меня есть опыт работы только с устройствами от фирмы Solinst. Не буду цитировать рекламные буклеты, просто скажу, что штука это невероятно удобная, настоятельно рекомендую попробовать — потом за уши не оттащишь. Стоят они однако весьма негуманных денег.
Однако у этих устройств кроме цены есть еще один недостаток, хотя и весьма условный: по-умолчанию замеры там сохраняются с постоянным шагом по времени, который хоть и задается в настройках, но при ОФР, очевидно, специалистами выставляется не очень большим. Как следствие на выходе мы получаем огромный CSV файл с кучей замеров с шагом в 5-10 секунд на протяжении 2-3-10 суток. Очевидно, что уже к концу первого часа такая частота замеров становится избыточной, особенно если возникнет необходимость в оформлении журнала откачки. 
Хочу поделиться своим инструментом для прореживания данных с логгера. Надеюсь, никому не потребуется инструкция к применению — там все предельно просто: вставляем исходные данные во вкладку Input, выставляем настройки во вкладке Options, копируем прореженный результат из вкладки Output.
Убедительная просьба: не надо присылать мне запросы на разрешение редактирования документа — все-равно не дам. Либо сделайте копию на свой Drive, либо сохраните в виде XLS на локальный диск. Все это можно сделать через меню File.

05/12/2019

Расчет количества блоков модели

Сначала о грустном: ввиду претензий правообладателей пришлось удалить рубрику «Книжная полка гидрогеолога». Жалко конечно, но ничего не поделаешь.
Дабы подсластить пилюлю хочу поделиться одним расчетом, который будет полезен при обосновании смет на моделирование. В последнее время, в экспертизе довольно остро встал вопрос обоснованности размеров модели, размеров блоков и их количества. К счастью, в СП 250.1325800.2016 есть конкретные (ну не совсем уж конкретные, но хоть что-то) указания, какие надо брать величины.
По ссылке можно скачать Расчет количества блоков геофильтрационной модели согласно рекомендациям СП 250.1325800.2016. Документ открыт только для чтения, поэтому для ввода своих значений вам надо либо его скачать на локальный диск (если сохраните в формате Excel, то надо будет поиском и заменой заменить в документе pow на power), либо сохранить на свой персональный Google Drive. Но можно оставлять комментарии, особенно если обнаружились ошибки. Для сметного расчета я использовал такие значения параметров, дабы цена получалась наименьшая (но можно ввести любые, только будьте готовы к проблемам с их обоснованием в сметной экспертизе).

19/03/2019

Учитесь программировать

На днях привалило мне «счастье» - 100 с гаком гигабайтов базы данных. Очень важных и очень интересных. Однако сразу всплыло несколько "но":
  1. База создана в Borland Interbase 7.5.
  2. Сами важные данные хранятся в виде BLOB-записей в графическом формате.
Сначала пришлось изрядно повозиться только для того, чтоб просто открыть базу, т.к. общедоступный Firebird этого формата не понимает. Потом выяснилось, что ни одна из популярных программ для работы с базами Interbase не умеет в пакетном режиме выуживать из оттуда BLOB-записи и сохранять их в виде файлов. Задача еще усложнялась тем, что мало выудить файл, надо еще его сохранить под правильным именем, которое хранится в отдельном столбце в таблице.
Делать нечего, пришлось экспортировать базу в формат SQLite. В принципе, на этом можно было бы и остановиться, SQLite вполне себе общеупотребимый и удобный формат, но опять проблемы:
  1. BLOB-записи из Interbase почему-то сконвертировались в HEX-строки SQLite (хотя тоже в формате BLOB). Т.е. из бинарного вида они превратились в ASCII: "04В458A4C1...". Т.е. напрямую с этими данными работать не получится — придется извлекать, конвертировать в бинарный вид и уже после этого открывать в графическом редакторе. Либо, как вариант, можно пройтись по базе и перекодировать данные.
  2. Сама идея хранить файлы в базе данных мне кажется сугубо порочной. Для файлов уже создана отличная база данных и называется она — файловая система.
  3. Программы, которые я подобрал для извлечения BLOB-записей из SQLite, отлично работают на маленьких базах, но на больших либо денег просят, либо падают от недостатка памяти при экспорте >1000-й строки (а всего там их около 80 тысяч), а часто и то и другое.
Помог скрипт sqlite-blob-dumper.py. Пришлось его немного модифицировать (о слава тем далеким годам, когда я занимался программированием), чтобы имена экспортируемых файлов брались из базы, а не генерировались на основе номера строки, названия базы и столбца с данными и все заработало. Причем довольно быстро и без проблем с памятью.
Скрипт уже выгрузил 80% базы, когда я понял, что могу сэкономить еще одну итерацию: ведь файлы сохранялись все в том же HEX-виде и для перекодирования в бинарный вид я собирался использовать отдельную программу xxd.exe — хорошую, но медленную под Windows. Тормознул скрипт, модифицировал, добавив в него функцию перекодирования. Теперь файлы сразу сохраняются в нужном бинарном виде — существенная экономия времени и дискового пространства.

P.S.: Для особо прозорливых, кто догадался, что это за база: извините, не моё — не дам.

08/06/2014

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

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

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

21/02/2013

Сертификат для Modflow дубль два

К моей давнишней записи появилось существенное дополнение.
На сайте сертификационного центра выложен список сертифицированных программ: csert.ru/list.php. Вот цитата из него:
Обозначение программной продукции Название программной продукции № сертификата соответствия и срок его действия. Нормативные документы, на соответствие которым выдан сертификат соответствия
Подкласс 2.41. Защита от опасных геофизических воздействий.
Программа PMWIN.
Версия 5.3
Геофильтрационные и геомиграционные расчеты на основе численных методов конечных элементов. РОСС RU.СП15.Н00551
(28.11.2012 – 27.11.2014)
СНиП 2.04.02-84*, СНиП 22-02-2003, СП 116.13330.2012 (СНиП 22-02-2003), СНиП 2.06.15-85
Таким образом, у одной из версий PMWIN (древней, но к счастью бесплатной) таки есть сертификат. Осталось раздобыть где-нибудь его копию, дабы от экспертов отвязаться раз и навсегда. Хотя бы по этому вопросу.

17/05/2012

Ссылки

Groundwater Science Textbook. Еще один вариант обработки откачек с помощью Excel.
EPA on-line tools for deliniating plumes and site assessment. Замечательная подборка онлайновых (даже так!) модулей с миграционными расчетами:

Formulas:

Models

14/04/2012

Русско-английский гидрогеоэкологический словарь

Тут у меня недавно спрашивали, где достать упомянутый мною гидргеологический англо-русский словарь. Я в ответ написал, что про существование такого не слышал. Однако, добрые люди указали мне на оплошность и прислали таки русско-английский гидрогеоэкологический словарь. Не англо-русский, конечно, но всяко лучше, чем ничего. Он немного неудобен в использовании — там все внутренние ссылки ведут на внешний сервер, который сейчас недоступен, я в ближайшее время это исправлю, а пока качайте как есть.

05/04/2012

Сертификат соответствия для Modflow

Есть у наших доблестных экспертов и представителей СРО такой пунктик: требовать от гидрогеологов сертификат соответствия и лицензию на расчетные программы. Вообще, согласно нашему законодательству, сертификация программных средств является добровольной и экспертам рекомендовано интересоваться наличием этих самых сертификатов соответствия. Насколько я помню, из всей массы программ для геофильтрационного моделирования сертификат есть только у Geolink ModTech, да и то не факт, что он до сих пор действителен. Modflow, являясь некоммерческой разработкой Геологической Службы США, такой сертификации не проходил. Кроме того, почему-то сертифицируют у нас не сам расчетный модуль (что было бы логично), а т.н. пре- и постпроцессоры. Вот, две полезные ссылки по теме:
  1. http://csert.ru/doc/r3.htm
  2. http://csert.ru/doc/letter.htm
Так что же делать, если эксперт требует этот несчастный сертификат. Первым делом — насторожиться, поскольку по моему личному опыту такой вопрос чаще всего возникает в случаях, когда эксперту поручено «зарезать» отчет. В качестве отписки могу лишь посоветовать отправить эксперту нижеследующий текст — иногда срабатывает. А вот что делать, если в уставе СРО прописано обязательное использование исключительно сертифицированных программных средств, я не знаю.
Справка об используемых программных средствах
В настоящее время наиболее разработанными и широко используемыми во всем мире для решения задач связанных с моделированием движения подземных вод являются программа MODFLOW и связанные с ней программы MODPATH, MT3DMS,PEST и некоторые другие. Начало разработки этого пакета программ положено более 20 лет и в настоящее время он является наиболее оттестированным и надежным. MODFLOW и связанные с ней программы фактически являются мировым стандартом для решения задач фильтрации и массопереноса. Большая часть задач, связанных с поземной гидросферой, решается на основе использования этих моделей. При порождении международных экспертиз часто требованием экспертов является именно использование этих программ.
В США MT3D и MODFLOW рекомендуются к использованию Агентством окружающей среды и Геологической службой. На сайтах этих организаций представлены тестовые примеры, подтверждающие правильность результатов расчетов, выполненных с помощью MODFLOW и связанных с ней программ.
Согласно Федеральному закону от 27 декабря 2002 г. N 184-ФЗ «О техническом регулировании»: сертификация делится на обязательную и добровольную. Программные средства для гидрогеологических расчетов не включены «Перечень продукции, подлежащей обязательной сертификации». Таким образом, сертификация таких программ является добровольной.
В России MODFLOW и связанные с ней программы широко применяются в таких уважаемых организациях как: ГИДЭК, МГУ, СПГУ, НИИ ВОДГЕО, Гидропроект и многих других. Для подготовки данных и визуализации результатов расчетов при использовании MODFLOW и связанных с ней программ в нашей организации используется следующее программное обеспечение:
  • PMWIN 5.1 — распространяется бесплатно.
  • ModelMuse 2.13.0.0 — распространяется бесплатно.
Что же касается лицензий на дорогостоящие программные комплексы типа PMWIN 8.0, Visual MODFLOW, GMS и т.п., то все что я могу посоветовать — это купить их. Других вариантов у меня нет.

Важное дополнение! На сайте сертификационного центра выложен список сертифицированных программ: http://csert.ru/list.php. Вот цитата из него:
Обозначение программной продукции Название программной продукции № сертификата соответствия и срок его действия. Нормативные документы, на соответствие которым выдан сертификат соответствия
Подкласс 2.41. Защита от опасных геофизических воздействий.
Программа PMWIN.
Версия 5.3
Геофильтрационные и геомиграционные расчеты на основе численных методов конечных элементов. РОСС RU.СП15.Н00551
(28.11.2012 – 27.11.2014)
СНиП 2.04.02-84*, СНиП 22-02-2003, СП 116.13330.2012 (СНиП 22-02-2003), СНиП 2.06.15-85
Таким образом, у одной из версий PMWIN (древней, но к счастью бесплатной) таки есть сертификат. Осталось раздобыть где-нибудь его копию, дабы от экспертов отвязаться раз и навсегда.

04/04/2012

Расчет самозатопления карьера

Решил поделиться одной своей наработкой из области горной гидрогеологии. Заранее оговорюсь: этот подход совершенно неуниверсален и крайне маловероятно, что получится его применить еще где-нибудь.

Постановка задачи:

Есть карьер, расположенный склоне небольшой горы (высотой примерно 800 метров). Глубина карьера составляет 550 м (абсолютная отметка дна карьера равна 410 м). Под горой протекает река с абсолютной отметкой уреза — 450 м. Расстояние от карьера до реки в среднем составляет 80 м. Периметр подгорной части карьера равен 100 м. Особенности рельефа нагорной части позволяют утверждать, что все атмосферные осадки в пределах области водосбора за вычетом испарения будут поступать только в карьер при любом уровне воды в нем. Области подземного и поверхностного водосбора совпадают. Водоупор находится на глубине 550 м и что важно — рельеф его кровли полностью повторяет рельеф поверхности.

Задача:

Надо посчитать за какой срок после окончания разработки и процедур осушения карьер наполнится водой и до какой отметки.

Решение:

Ну, с водопритоком с нагорной части все понятно:

Q_n=\left (\omega -\varphi \right )\cdot S...(1)
где:
  • ω – среднегодовые атмосферные осадки;
  • φ – испарение;
  • S – площадь водосбора.

Воде на склоне некуда деваться и она либо поверхностным, либо подземным стоком неминуемо притечет в карьер. Чтоб увязать Qn и q1 из рисунка — надо разделить Qn на периметр нагорной части карьера, но для расчета это не нужно.

Водоприток с подгорной части происходит за счет фильтрации подземных вод в/из реки. Приток (отток) подземных вод связан с уровнем воды в реке следующим выражением (при отсутствии поверхностнго стока, см. рисунок):

q_2=k\frac{h_{0}+h_{L}}{2}\cdot\frac{H_0-H_L}{L}-\frac{b(\omega -\varphi )}{2}L...(2)
где:
  • q2 – удельный расход потока;
  • k – коэффициент фильтрации;
  • h0 – мощность потока вблизи карьера;
  • hL – мощность потока вблизи ручья;
  • H0 – уровень воды в карьере;
  • HL – уровень воды в ручье;
  • L – расстояние между карьером и ручьем;
  • b – коэффициент стока.

В данном выражении я увязал вместе две известных зависимости для фильтрации между двумя разноуровенными водоемами: при наклонном водоупоре и при наличии инфильтрации.

Формула приведена для расчета удельного расхода потока, для получения величины водопритока (Qp) в карьер удельный расход должен быть умножен на периметр подгорной части карьера (P):

Q_p=q \cdot P...(3)

Дальше, казалось бы, задача превращается в очевидную: приравнять выражения (1) и (2) с учетом (3), но увы, тут нас ждет подвох. Дело в том, что выражении (2) присутствует величина H0 (уровень воды в карьере), которая вообще-то зависит от объема притекшей на заданный момент времени воды в карьер, т.е. эта величина динамически изменяется в зависимости от текущего соотношения Qp и Qn. К тому же h0 ( мощность потока вблизи карьера) равна в нашем безнапорном случае:

h_0 = H_0 - H_v...(4)
где:
  • Hv – абсолютная отметка кровли водоупора на краю карьера, обращенном к реке;

Таким образом, мы приходим к классической ситуации: расход зависит от напора, а напор от расхода. Решение я опущу в силу его громоздкости, но поясню суть (если кому интересно, xls-ку можно скачать по ссылке: drive.google.com..., только не просите меня объяснить, как оно работает). Основные положения:

  1. Зависимость уровня от объема поступившей воды — ступенчатая, в силу ступенчатости бортов карьера (на рисунке этого не видно, я поленился рисовать уступы).
  2. Между «ступеньками» эта зависимость — линейная. Ну, тоже понятно, борта уступов почти линейны и имеют небольшой наклон. Тут я просто вытряс из проектировщиков объемы карьера на заданных отметках (кромках уступов) — предложил было аппроксимировать борт карьера усеченным конусом, но получил по рогам.
  3. Выбираем шаг по времени Δt.
  4. По формулам (1) и (2,3) считаем объем воды (Qn + Qp), поступивший в карьер к моменту времени Δt.
  5. По выявленной ранее зависимости уровня воды в карьере от объема поступившей в него воды получаем новый уровень H0.
  6. Считаем водоприток, подставляя в выражение (2) уже новое значение H0. Объем поступившей воды в карьер соответствует моменту времени 2Δt.
  7. Повторяем последние два пункта пока не достигнем «стационара», т.е. до того момента, пока водоприток в нагорной части не уравняется с оттоком в сторону реки. Получаем искомую равновесную отметку воды в карьере и момент времени nΔt, когда произошло это радостное событие.
  8. Попутно отмечаем момент, когда q2 в выражении (2) меняет знак, т.е. подземный приток от реки сменяется на отток.

Краткий итог:

  • Насколько я знаю, такая схема называется явной. Очевидно, она чувствительна к выбору шага по времени Δt.
  • Важно помнить, что равновесный уровень воды в карьере теоретически может оказаться выше нижней кромки карьера, т.е. вода просто начнет переливаться через край и уйдет в реку поверхностным стоком (у меня на одном из карьеров так и случилось), и тут возникает большая неприятность с величиной ω-φ, но учитывать еще и это мне показалось совсем уж крохоборством.
  • А еще там периметр P зависит от H0 (при подъеме уровня в образующемся озере, очевидно, растет и его периметр) — это я тоже не стал учитывать, а вообще следовало бы т.к. результат оказался очень чувствителен к этому параметру.
  • Этот нехитрый расчет вроде бы не очень сложно запихнуть в тот же MODFLOW, но я пока не обдумывал детали — наверняка всплывет какая-нибудь заковыка.

Вот и все. Если что непонятно — спрашивайте. Заметили ошибки — ругайте.

22/03/2012

Excel для гидрогеологов

Наткнулся на еще одну очень интересную подборку полезных для гидрогеолога табличек в Excel. Что характерно: опять американцы и опять USGS. А остальным должно быть стыдно, да.
Кстати, вот еще одна интересная программа (это даже не программа, а большой такой макрос для MS Office Access), предназначенная для оптимизации расположения наблюдательных скважин при мониторинге распространения загрязнения — MAROS. Бесплатная, если не врут.
The Monitoring and Remediation Optimization System (MAROS) software has been developed by GSI and UH for the Air Force Center for Environmental Excellence (AFCEE) in accordance with the AFCEE Long-Term Monitoring Optimization guide. The software provides site managers with a strategy for formulating appropriate long-term groundwater monitoring programs that can be implemented at lower costs. MAROS is a decision support tool based on statistical methods applied to site-specific data that accounts for relevant current and historical site data as well as hydrogeologic factors (e.g. seepage velocity) and the location of potential receptors (e.g., wells, discharge points, or property boundaries). Based on this site-specific information the software suggests an optimization plan for the current monitoring system in order to efficiently achieve the termination of the monitoring program. This public domain software was developed for the Air Force Center for Environmental Excellence through a Broad Agency Announcement (BAA) contract with GSI Environmental Inc. Inc. of Houston, TX and the University of Houston. MAROS is programmed in Microsoft Access.

24/02/2012

Обработка откачки-восстановления из одиночной скважины по схеме Тейса-Джейкоба

Вот, как и обещал, добавил обработку данных откачки-восстановления с использованием схемы Тейса-Джейкоба (однородный изолированный неограниченный пласт).
Скачать можно тут: pump_Theis_Jacob.7z

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

14/02/2012

Шаблон журнала откачки

Вот, решил поделиться «наработкой». Простая табличка в формате MS Excel. Без каких-либо расчетов (есть мыслишка адаптировать мериканские таблицы обработки откачек для наших реалий), а пока только шаблон для заполнения «в полях». Вроде удобно, сам проверял. pumping.7z

02/02/2012

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

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

04/12/2011

Несколько интересных ссылок

A groundwater model is a hypothesis; A hypothesis can be disproved by never proven. Как-раз то, что я всегда ощущал на уровне понимания, но так и не додумался сформулировать. Геофильтрационная модель - не истина в последней инстанции, а всего-лишь предположение ее создателя о том, как все там под землей устроено. И это предположение нельзя доказать - только опровергнуть.
“The word model has so many meanings and is so overused that it is sometimes difficult to know what one is referring to. For this discussion, we define a model as a representation of a real system or process…We define a conceptual model as a hypothesis for how a system or process operates. The idea can be expressed quantitatively as a mathematical model.”


Sharing MODFLOW Models on ArcGIS Online. Что творят, черти. Кажется, действительно пора изучать ArcGIS. А то я со своей мапинфой совсем от жизни отстану.
Groundwater simulation models are mathematical models that help us understand and predict how water behaves in subsurface aquifers. Groundwater models, such as MODFLOW, are commonly used in the industry to predict different aspects of groundwater management including water availability, effects of pumping, contaminant transport, and more.

One of the challenges facing organizations is the communication of information regarding groundwater systems from professional modelers to others within the organization or to the public. Usually, visualization of groundwater models is achieved with specialized software packages. The model input and output files are a set of text or binary files that are difficult to interpret and map without specialized tools.

The workflow presented enables organizations to quickly create maps and GIS layers of MODFLOW models and to publish those on ArcGIS Online. By posting on ArcGIS Online you can share MODFLOW inputs and outputs with anyone within your organization, or share the model results with the public.

17/11/2011

Несколько полезняшек для гидрогеолога

Carlos Molano опубликовал на своем сайте несколько полезных таблиц с гидрогеологическими расчетами:

1D analytical transient flow in confined aquifers
Automatic fit for pump test data to the Theis equation
Capture type curves for wells Hantush leaky aquifer test analysis
Inverse particle tracking
Neuman unconfined aquifer test analysis
Piper diagram
Superposition Thiem well field
The Theis test for confined aquifers
The Theis fit for variable pumping rates
The Theis-Jacob Test