14/11/2012

О качестве инженерно-геологических изысканий

Просто две ссылки и предложение поразмыслить над сложившейся в отрасли ситуации.

  1. Ссылка нумеро уно.
  2. Ссылка нумеро дуо.

01/11/2012

Processing MODFLOW. Параметры.

К параметрам модели относятся:
  • Горизонтальный коэффициент фильтрации и Проводимость Horizontal Hydraulic Conductivity and Transmissivity.
    Коэффициент фильтрации должен быть задан для безнапорных и напорно/безнапорных слоев с переменной проводимостью (тип 1 и 3). Проводимость может быть задана для напорных и напорно/безнапорных слоев с постоянной проводимостью (тип 0 и 2).
    Горизонтальный коэффициент фильтрации в общем случае характеризует проницаемость вдоль модельных строк. Проницаемость вдоль модельных столбцов получается путем умножения этой величины на коэффициент горизонтальной анизотропии (anisotropy factor, задается в диалоговом окне Layer Property dialog box).
    Для безнапорного слоя проводимость считается как произведение коэффициента фильтрации на разность уровня воды и подошвы слоя (если в том же Layer Property dialog box явно не задана опция user-specified transmissivity).
    Важно понимать, что «внутре» MODLFOW всегда коэффициент фильтрации пересчитывает в проводимость. Просто в случае с безнапорным слоем проводимость уточняется на каждой итерации (точно не помню, но кажется этот момент тоже настраивается — можно пересчитывать проводимость не на каждом расчетном шаге).
  • Перетекание и вертикальный коэффициент фильтрации Vertical Leakance and Vertical Hydraulic Conductivity.
    Параметр перетекания расчитывается MODFLOW по формуле:
    VCONT=\frac{2}{\frac{m_k}{(k_z)_{k,i,j}}+\frac{m_{k+1}}{(k_{z})_{k+1,i,j}}}
    где:
    mk - мощность k-ого слоя;
    (kz)k,i,j - его вертикальный коэффициент фильтрации.
    Вы можете задать VCONT явно (Vertical Leakance), либо в виде коэффициента фильтрации (Vertical Hydraulic Conductivity). Казалось бы, зачем в здравом уме заниматься ручным расчетом перетекания, если можно не мучаться и задать коэффициент фильтрации — однако, есть ряд задач, где это необходимо.
  • Коэффициент вертикальной анизотропииVertical Anisotropy.
    Ну, тут все просто — это отношение горизонтального коэффициента фильтрации к вертикальному. Поддерживается версиями MODFLOW-2000 и старше. Сильно упрощает процесс решения обратных задач геофильтрации — подбирая проницаемость толщи не приходится синхронизировать значения коэффициентов фильтрации по разным осям.
  • Эффективная пористость Effective Porosity.
    Не буду напоминать определение эффективной пористости — она есть в любом учебнике. Отмечу лишь, что на решение задач геофильтрации этот параметр не влияет, в отличие от миграционных задач.
  • Упругая емкость (породы) Specific Storage, Упругая емкость пласта Storage Coefficient и Гравитационная емкость Specific Yield.
    Эти параметры используются при работе с нестационарными моделями. Упругая емкость пласта — это упругая емкость, умноженная на мощность пласта, и вы ее можете задать явно, либо предоставить выполнение операции умножения программе. Для безнапорных и напорно/безнапорных слоев надо еще задать гравитационную емкость Specific Yield.
За подробностями и определениями отправляю вас к литературе:
  1. Wen-Hsing Chiang: “Processing Modflow PRO”, April 6, 2006;
  2. Simcore Software: “Processing Modflow. An Integrated Modeling Environment for the Simulation of Groundwater Flow, Transport and Reactive Processes”, July 5, 2012;
  3. Шестаков В.М.: «Гидрогеодинамика», 1995 г.

30/08/2012

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

Прошу прощения за длительное отсутствие. Работы много.

16/07/2012

Моделирование в понижениях

Коллега в телефонной беседе напомнил про один очень полезный и неочевидный трюк. Как известно, нестационарные модели очень чувствительны к начальному распределению напоров. Получить это распределение можно двумя путями: интерполяцией фактически наблюденных напоров в наблюдательных точках (скважинах); либо с помощью откалиброванной стационарной модели. Очевидны недостатки обоих подходов: для первого не всегда хватает исходных данных, а второй чреват ошибками (все упирается в качество решения обратной задачи).
Однако, существует массивный пласт задач, для которых совершенно не обязательно знание о распределении напоров. К примеру: определение водопритока к котловану. Нам совершенно неважно какой там напор, в качестве результата нас устроит достигнутое понижение и расход.
В этом случае допустим такой приём: заменяем в нашем расчете термин «напор» на термин «понижение». Таким образом, «начальный напор» становится «начальным понижением», которое, очевидно, равно нулю. Все модельные слои должны быть напорными (Confined), даже если они не являются таковыми. Важно помнить: изолинии и массивы результатов, получаемые при таком расчете, будут показывать не напоры, а понижения. Расходы будут такими же, как и при расчете в термине «напор».
Совершенно ясно, что и этот подход не лишен недостатков и ограничений. Так, без дополнительных ухищрений невозможно учесть эффект осушения пластов, а величины расходов могут оказаться существенно завышенными для безнапорных горизонтов, как это всегда бывает, если их задать напорными.

11/07/2012

Симуляция откачек на сеточных численных моделях

Когда я только начинал карьеру гидрогеолога, мне попалась достаточно интересная задача — надо было более или менее автоматизировать процесс обработки экспресс-откачек. Ясное дело, что там все элементарно можно было сделать с помощью аналитических формул и «экселя», но я не стал искать легких путей и зарядил все это дело в MODFLOW. Результатом трудов стала статья в сборник работ института «ВОДГЕО», ознакомиться с которой вы можете по ссылке: Определение фильтрационных параметров по данным экспресс-откачек с использованием численного моделирования.pdf.
Недавно обнаружил, что я не один такой «хитрый»: Rushton, K. R. & Chan, Y. K., 1976, A numerical model for pumping test analysis, Proc. Instn Civ. Engrs, Part 2, 61, pp. 281-296 and Barrash, W. and M.E. Dougherty, 1997, Modeling axially symmetric and nonsymmetrical flow to a well with MODFLOW, and application to Goddard2 well test, Boise, Idaho, Ground Water ,Vol. 35, No. 4: pp. 602-611 — http://info.ngwa.org/gwol/pdf/972062742.PDF

25/06/2012

Modflow.ru

Скиберсквотил вкусный домен: modflow.ru. Что делать с ним - пока не решил. Пока поставил заглушку. Со временем планирую там сделать более или менее полноценный сайт-визитку, но это дело не ближайших месяцев.
Если кому нужен почтовый ящик типа ваше_имя@modflow.ru, пишите на water+blog@alick.ru, или заполните форму на сайте.

21/06/2012

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

An Alternative to Low Flow Groundwater Sampling: HydraSleeves. Оказывается, даже в такой консервативной области, как отбор проб воды из скважин, могут быть новые разработки. Как вспомню, как мы желонкой 50-метровую скважину «прокачивали» — в дрожь бросает.
Benefits of the HydraSleeve:
  • You can sample for all compounds (not just VOCs)
  • You can sample a well in less than 15 minutes
  • You can cut costs by 50% to 80%
  • You can reduce purge water to a big fat goose egg
Additional benefits of the HydraSleeve can be found here: http://www.hydrasleeve.com/learn-about-the-hydrasleeve.
Groundwater Sampling in a Snap. Еще одна ссылочка на вышеобозначенную тему.
Back in March, I wrote a post on the Joys of Low Flow Groundwater Sampling (or not). Instead, I reviewed an alternative to this procedure — the Hydrasleeve. At the same time, I came across another option for those of you out there that get as much enjoyment from low flow sampling as you do paper cuts, thistles, and flaming dog poo: The Snap Sampler.
Modeling Salt Water Intrusion: the Freshwater/Salt Water Interface. Меня тут недавно спрашивали об особенностях моделирования потока соленых и пресных вод. так вот, MODFLOW в чистом виде для этого не годится. Нужны другие модели. Типа той, о которой идет речь в этой заметке.
Due to the difference in density, there is a fairly sharp interface between freshwater and sea water in the subsurface. Thus, as freshwater flows out to the sea, it floats on top of a sea water wedge. The region of mixing at this interface is quite thin, leading to fairly sharp devide between the two fluids. In general, the brackish water that develops in this mixing zone flows outward to sea, rather than contaminating the freshwater. Under steady-state conditions, this situation is fairly stable, but pumping freshwater from the aquifer can change the pressure regime of the subsurface. This may lead to seawater being drawn into the freshwater aquifer and permanently degrading the aquifer’s water quality.
World average annual evapotranspiration web map now available! В то время, пока наши ФГУПы и ТИСИЗЫ секретят и ДСПсят все что только можно и нельзя, мериканцы выложили во всеобщий доступ интерактивную карту величин среднегодовой эвапотранспирации. Посмотрел цифры для Москвы и области — вроде правдоподобные.
Most of us understand the hydrologic cycle in terms of the visible paths that water can take: rainstorms, rivers, waterfalls, swamps, etc.  However, an even larger volume of water flows through the air all around us in two invisible paths: evaporation and transpiration. They claim 61% of all terrestrial precipitation, and together are referred to as evapotranpsiration. In order to better understand this important process, Esri’s Mapping Center has produced a web map showing the world’s average annual evapotranspiration.
Carlos Molano’s YouTube Channel on Groundwater Modeling Spreadsheets. Я давно слежу за творчеством Карлоса Молано. Когда-то я наивно считал себя чуть ли не гуру расчетов в Excel-е. Как же я ошибался: по сравнению с Карлосом я вообще ничегео не умею.
Last year, I wrote about the groundwater modeling spreadsheets developed by hydrologist Carlos Molano.  Molano has published a series of these spreadsheets on his personal website.

I recently learned that Molano also has a youtube channel, where he discusses using spreadsheets to solve groundwater modeling problems.

15/06/2012

Ссылки

SUTRA, an open-source option for modeling groundwater flow on a mesh. Вот, еще одна моделька появилась.
UTRA is a groundwater modeling code published by the United States Geological Survey. In many ways, it is similar to MODFLOW. SUTRA solves equations that simulate groundwater flow, solute transport with adsorption and decay, and thermal energy transport in groundwater. SUTRA can handle both steady-state or transient simulations. Physical properties can be projected on nodes or elements, or can be represented as a functions. Fluid density and viscosity are functions of temperature and solute concentration. Unsaturated flow properties are handled by a subroutine, UNSAT. Each node and element passes properties to and from UNSAT, which models flow through a series of unsaturated flow functions. Like MODFLOW, SUTRA communicates with input and output text files and runs in Fortran.
MODFLOW-USG, coming soon. А скоро и еще одна будет.

MODFLOW-USG (unstructured grid) is a new version of MODFLOW being developed by Sorab Panday of GMS.

According to the Aquaveo forum, MODFLOW-USG will probably be fully released at some point in 2012.

Водопровод. Интересный сайт с заметками о водоснабжении, в том числе и о гидрогеологическом моделировании. Вроде как блог, но логики в навигации сильно недостает.