Владимир Останин
Добрый день, коллеги!
Прошу совета по вопросу по гидрогеологическому моделированию.
Задача: сделать модель для для подсчета запасов подземных (грунтовых) вод для водозабора в долине небольшой горной реки. Водоносный горизонт один в четвертичных отложениях, гидродинамически связан с рекой. Уровни по 30 скважинам линейно зависят от гипсометрии с коэффициентом корреляции, близким 0.9. Уровни реки и искусственного водоема гипсометрически выше водозабора укладываются в тенденцию.
Решение: для моделирования фильтрации в стационарном режиме строю карту начальных уровней от гипсометрии беру изопьезу выше искусственного водоема в качестве верхней границы первого рода, уровень реки в качестве нижней границы первого рода, решаю обратную задачу (калибрую) с определением коэффициента фильтрации, водопроводимости ложа реки и водного баланса. Затем собираюсь задать с помощью модуля recharge удаленную область питания и посчитать модель с работающими скважинами водозабора.
Вопроса два – 1) содержательный и 2)технический.
1) Корректно ли вообще такое решение, поскольку при сопоставлении результатов моделирования стационарной фильтрации, уровня осадков и модуля стоков получил область питания для осадков (без испарения) 10 км2, а по модулю стока – 20,4 км2, что превышает площадь области питания, если ее выделять по водоразделу.
2) В качестве программного средства для калибровки и решения обратной задачи пытаюсь использовать последние версии ModelMuse и ModelMate, которые якобы интегрированы между собой согласно хелпам. При попытке запустить из ModelMate оценку параметров в приложении UCODE не читаются названия параметров из модели, пробовал задавать вручную – задаются, но их UCODE в модели нe видит. Чтобы посоветовали делать в таком случае? Можно ли как-то экспортировать модель в понятный для PMWIN формат? Или проще вообще использовать PEST из консоли? Как вообще можно узнать названия параметров модели Modflow (Kx, Conductance для реки, K упругоемкости) по умолчанию? Пробовал использовать названия из файла lst – не получаетсяL(
Заранее благодарен за любой ответ.