ПОИСК Статьи Рисунки Таблицы Особенности численного моделирования миграции подземных вод переменной плотности и некоторые частные результаты из "Проблемы гидрогеоэкологии Том 1" В ряде работ доказывают преимущества первого подхода, когда значения давлений (или приведенных напоров и концентраций) рассматриваются как независимые переменные [19]. Для решения могут использоваться любые из рассмотренных выше методов — МКР, МКЭ, МХ. [c.406] Важным физическим результатом моделирования явилось доказательство динамического характера наблюдаемого (квазистационарного) предельного равновесного состояния в зоне соленых вод образуется конвективная ячейка — в нижней части языка вода течет от моря, а в верхней — в обратном направлении. [c.407] Схема решения такого рода задач в наиболее общей — трехмерной — постановке предложена в работе [15 ], Уравнения для приведенных напоров и концентраций решаются на каждом шаге последовательно, методом конечных элементов. Соответствующие матричные уравнения связаны друг с другом через скорости фильтрации и эффект всплывания 7 (задача нелинейная). Для устранения трудностей, связанных с исследованием трехмерного процесса, решение ищется на ряде последовательных слоев по координате I. [c.408] НОЙ И вынужденной конвекции, коща на верхней границе купола возникает весьма чувствительный баланс между поперечной профильной дисперсией (стремящейся обеспечить отток растворенных солей вверх) и нисходящей плотностной конвекцией. Численная программа TOUGH, основанная на методе сопряженных градиентов, обеспечивала расчет скоростей и массовых градиентов на границах блоков не только в ортогональном, но и в параллельном им направлении, что давало повыщенную точность реализации поперечной дисперсии. Были также включены утонченные процедуры для проверки текущих результатов счета на ложную сходимость, для чего итерации на каждом щаге продолжались до достижения весьма малых невязок по балансу массы, аккумулируемой в каждом блоке. Результаты моделирования оказались даже качественно отличными от полученных ранее с помощью других численных программ и исключительно чувствительными к численной поперечной дисперсии, особенно в периоды времени, близкие к стабилизации процесса. Для достижения стационара требуются времена порядка сотен лет, причем при сравнительно слабой поперечной дисперсии вообще не возникает конвективной ячейки и растворенная соль выносится на поверхность в области разгрузки инфильтрующихся вод. Как характер процесса, так и время достижения стационара сильно зависят от начальных условий при задании исходного ореола рассолов над соляным куполом, их вовлечение в конвективный процесс происходит крайне медленно и общее время ста-билизации увеличивается в несколько раз (по сравнению с вариантом моделирования, при котором в начальный момент естественные рассолы над куполом отсутствуют). [c.409] Подбор оптимальной дискретизации модели практически невозможен без специальных тестов, которые следует начинать со стационарных задач без переходной зоны. Временной шаг предпочтительнее принимать постепенно возрастающим во времени, начиная с величин, отвечающих занижению критерия Куранта (7.9) примерно на порядок. Это требование оказывается особенно существенным для блоков, прилежащих к скважинам. Необходима и проверка качества моделирования расходов субвертикальных потоков, обусловленных плотностной конвекцией, что требует специальных тестов, которые лучше всего ориентировать на типовые условия естественной конвекции. [c.411] Вернуться к основной статье