Счетные осцилляции при расчете взаимодействия ППВ могут воз- никать также и при «критическом» состоянии поверхностного стока, когда его расход близок к нулю [Гриневский, 1994]. Если в этот пери- од взаимосвязь ППВ происходит при свободном режиме фильтрации под рекой (см. 1.2), то, согласно (1.5), расход фильтрации из реки оп- ределяется разностью напоров (Н, - Но) = Й, + о (рис. 1.2). При обме- лении реки Л, — 0, и (Н, - Но) — то, а уДеЛЬНЫИ расход х —› ат Т.е. / — К. Когда участок реки полностью обмелел (Р = 0 и й, = 0), фильт- рация из реки прекращается: г = (, что означает разрывность расхода взаимодействия 1В (аналогично функции Хэвисайда) — скачок на величину, пропорциональную коэффициенту фильтрации экрани- рующих ложе реки отложений &) Учитывая, что такой резкий скачок удельного расхода приводит к адекватному изменению расхода взаимодействия ПШВ в блоке & — на всем расчетном участке реки пропорционально /Р, это может выра- жаться как в осцилляциях уровня подземных вод Й, так и в «мерца- нии» расхода реки Р. Физическая причина таких осцилляций связана с недоучетом в расчетной модели емкостных свойств экранирующих русло отложений, когда они не представлены самостоятельным слоем модели (см. 1.6). Для преодоления таких осцилляций рекомендуется использование сглаживающих функций [Гриневский, 1994], посте- пенно «выводящих» тр из разности напоров ППВ при приближении к полному обмелению реки: если И, —› 0, то (Н, - Но) — 0. Условие полного исчезновения: речного стока (Р = О ий, = 0) тре- бует балансового контроля расчетного расхода взаимодействия ППВ К, чтобы он не оказался больше суммарного расхода поверхностных вод, поступивших в блок [Гриневский, 1994]; если В, > Р_| + Р' + (М ‚1), + О,., то В, = Р_| + Р +(М,‚);+0, и Р =0. { Это также может приводить к счетным осцилляциям, которые ис- правляются уменьшением расчетного шага по времени. Такой кон- троль производится на каждом шаге по времени, что позволяет моде- лировать восстановление речного стока, при невыполнении условия (4.19). В практике гидрогеологического моделирования широко исполь- зуется расчетный прием, основанный на применении метода суперпо- зиции [Лукнер, Шестаков, 1976], согласно которому решение уравне- ния геофильтрации ведется относительно изменения (понижения) на- поров $, по сравнению с начальными (на момент времени / = 0): $ = й — /. Такой прием позволяет «напрямую» оценивать балансово- гидрогеодинамическое воздействие, связанное с техногенными пони- жениями уровня подземных вод, что особенно эффективно при реше- нии геоэкологических задач — в частности, при оценке воздействия эксплуатации подземных вод на речной сток [Гриневский, Штенгелов, 1988]. (4.19) 86