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}}}
  где:
  m_k - мощность k-ого слоя;
  (k_z)_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 г.