Просто две ссылки и предложение поразмыслить над сложившейся в отрасли ситуации.
Прикладная гидрогеология, геофильтрационное моделирование, инженерная геология. Блог специалиста по гидрогеологическому моделированию Копылова Александра.
14/11/2012
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.
- Wen-Hsing Chiang: “Processing Modflow PRO”, April 6, 2006;
- Simcore Software: “Processing Modflow. An Integrated Modeling Environment for the Simulation of Groundwater Flow, Transport and Reactive Processes”, July 5, 2012;
- Шестаков В.М.: «Гидрогеодинамика», 1995 г.
Метки:
Modflow,
PMWin,
моделирование,
руководство,
схематизация
30/08/2012
Гидрогеологический дайджест
Прошу прощения за длительное отсутствие. Работы много.
- Hydrogeology - A Global Perspective. Сборник статей о гидрогеологии, доступный для бесплатного скачивания.
- MODFLOW-USG (UnStructured Grids): A Game Changer for Groundwater Modeling. Теперь MODFLOW умеет не только ортогональные сетки. Последние преимущества у конечно-элементных моделей отнимают.
- Aquifer Testing Reference List. Подборка ссылок на статьи и книги, посвященные опытно-фильтрационным работам.
16/07/2012
Моделирование в понижениях
Коллега в телефонной беседе напомнил про один очень полезный и неочевидный трюк.
Как известно, нестационарные модели очень чувствительны к начальному распределению напоров. Получить это распределение можно двумя путями: интерполяцией фактически наблюденных напоров в наблюдательных точках (скважинах); либо с помощью откалиброванной стационарной модели. Очевидны недостатки обоих подходов: для первого не всегда хватает исходных данных, а второй чреват ошибками (все упирается в качество решения обратной задачи).
Однако, существует массивный пласт задач, для которых совершенно не обязательно знание о распределении напоров. К примеру: определение водопритока к котловану. Нам совершенно неважно какой там напор, в качестве результата нас устроит достигнутое понижение и расход.
В этом случае допустим такой приём: заменяем в нашем расчете термин «напор» на термин «понижение». Таким образом, «начальный напор» становится «начальным понижением», которое, очевидно, равно нулю. Все модельные слои должны быть напорными (Confined), даже если они не являются таковыми. Важно помнить: изолинии и массивы результатов, получаемые при таком расчете, будут показывать не напоры, а понижения. Расходы будут такими же, как и при расчете в термине «напор».
Совершенно ясно, что и этот подход не лишен недостатков и ограничений. Так, без дополнительных ухищрений невозможно учесть эффект осушения пластов, а величины расходов могут оказаться существенно завышенными для безнапорных горизонтов, как это всегда бывает, если их задать напорными.
Однако, существует массивный пласт задач, для которых совершенно не обязательно знание о распределении напоров. К примеру: определение водопритока к котловану. Нам совершенно неважно какой там напор, в качестве результата нас устроит достигнутое понижение и расход.
В этом случае допустим такой приём: заменяем в нашем расчете термин «напор» на термин «понижение». Таким образом, «начальный напор» становится «начальным понижением», которое, очевидно, равно нулю. Все модельные слои должны быть напорными (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
Недавно обнаружил, что я не один такой «хитрый»: 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, или заполните форму на сайте.
Если кому нужен почтовый ящик типа ваше_имя@modflow.ru, пишите на water+blog@alick.ru, или заполните форму на сайте.
21/06/2012
Гидрогеологический дайджест
An Alternative to Low Flow Groundwater Sampling: HydraSleeves. Оказывается, даже в такой консервативной области, как отбор проб воды из скважин, могут быть новые разработки. Как вспомню, как мы желонкой 50-метровую скважину «прокачивали» — в дрожь бросает.
Benefits of the HydraSleeve:Groundwater Sampling in a Snap. Еще одна ссылочка на вышеобозначенную тему.Additional benefits of the HydraSleeve can be found here: http://www.hydrasleeve.com/learn-about-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
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.
Метки:
excel,
Modflow,
гидрогеохимия,
карты,
моделирование,
ссылка,
устройства,
эвапотранспирация
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.
Subscribe to:
Posts (Atom)