где Р (0) — коэффициент . диффу‘зии влаги, или влагопроводности, аналогичный, как видно из выражения ([8.1]1), по физическому `смыслу коэффициенту уровнепроводности 56]. Уравнения (8.8), (8.10) и (8.12) могут содержать, как и дифферен- циальные уравнения движения грунтовых вод, свободный член, ха- рактеризующий интенсивность истока & в виде отбора влаги корнями растений. Величина ё связана с удельной поверхностью корневых во- лосков г на некоторой глубине и суммарной транспирацией Е, (в — = Е,г) и зависит от глубины и времени, поскольку корневая система развивается со временем. Величина ё может характеризовать также интенсивность водообмена крупных трещин, по которым в первую очередь движется инфильтрующаяся вода, с агрегатами почвы и бло- ками пород зоны аэрации при гетерогенно-блоковой ее структуре. Дифференциальные уравнения (8.10) и (8.12) являются нелинейными и их коэффициенты зависят от Н, } или 0. Решение их связано с серь- езными трудностями. Существуют различные способы линеаризации этих уравнений, ‘сводящиеся к усреднению параметров влагопереноса [40, 47 ]. Аналитических решений таких уравнений мало. Нелинейные уравнения решаются численным методом с примене- нием АВМ или ЭЦВМ [56]. Однако и такое решение весьма трудо- емко, так как требуются большой объем исходной информации и ВЫ- полнение итерационных приближений в решении нелинейного урав- нения. 8.6. РАСЧЕТНЫЕ СХЕМЫ ПРОЦЕССОВ ПИТАНИЯ, НАСЫЩЕНИЯ И ОСУШЕНИЯ Рассмотрим схемы процессов питания, осушения и насыщения и их динамические характеристики. Составим балансовое уравнение, хХа- рактеризующее приток @ к свободной поверхности со стороны зоны аэрации, отток 0 от нее в грунтовые воды и приращение уровня й: о - = @—. (8.13) ° @ Полный водный баланс элементарного объема высотой /, (от поверх- ности земли до первого регионального водоупора или до некоторой границы в зоне насыщения) включает поступление воды интенсив- ностью / через поверхность земли, отток интенсивностью © через нижнюю границу и изменение влагозапасов И в элементе за @. Учитывая переменную по высоте влажность @, получим [56] 1 [--о=- " д [ 6а. @й дЕ 0 Рассматривая приращение влагозапасов раздельно. в зоне насы- щения (от О до іъ) и в зоне аэрации (от Й до /.) е° 1 ' 1 баг == —- | 04г - — | баг Г д ё Т д [ 145