Файл: Количественные методы в мелиорации засоленных почв..pdf

При инфильтрации (поливе):

У—Ср01 при h = H + h k.

3. Уравнение распределения концентрации солей в во­ доносном слое (0<.х<.1, О <Л<Я, 0<£=S^koh)

Условия согласования на поверхности потока грунтовых вод:

■Dill = D2ff- ,У= С при h = H(x, t),

h — вертикальная, х — горизонтальная, t — временная ко­ ординаты.

Приведенная система уравнений решается методом ко­

АЛГОЛ-60, использованы следующие обозначения:

I — длина расчетного участка, [/и]; N — число горизон­

п, i — целые числа для индексации массивов значений ин­ тервалов сетки, концентраций и временных интервалов; П — число вертикальных интервалов сетки в области пото­

верхности потока грунтовых вод (пъезометрический напор)

дН

в узлах сетки, [лг], (г = 0,1,. . . , N); d H [ i ] = k j —— скорость

потока грунтовых вод в расчетных узлах (с обратным зна­

159

тп1 = тп\, m2 —m2— активные пористости зоны аэрации и во­

ент уровнепроводности, [м2/сут] ; Ш = Нт1П— наименьшая толщина потока грунтовых вод на участке, [ж]; у п = у п — критическое значение концентрации солей в почвенном растворе, выше которого процессы сорбции и кристаллиза­ ции преобладают над процессами растворения и десорбции, [г/л];

L[i] =L(iAx) — высота дневной поверхности в расчетных уз­

условии уравнения (2), [ ж ] ; б — характерная толщина слоя интенсивного испарения и поглощения влаги корнями ра­ стений в верхней части капиллярной каймы, (б ~ 0,01 ж).

160

Следующие обозначения используются для выполнения ЭВМ расчетов, связанных с функциями времени,

а. Параметры инфильтрации или поливов.

Число периодов инфильтрации предусмотрено не более шести в год, учитывая наиболее значительные по объему (не промывной режим!). Схема задания параметров полива приведена на рис. 15 со следующими обозначениями: W \n\=W n— время начала очередного периода инфильтра­ ции в сутках от начала года, л = 1, 2 , ... , 6; а [ п \ = а п— про­ должительность 7г-го периода инфильтрации (полива) в сут­ ках, 71 = 1, 2 , .. . , 6; qa\ii\=qan— средняя скорость инфиль­ трации В 7 7 -Н Ы Йпериод, [м/сут] .

Рис. 15. Схема задания параметров полива.

б. У ровни грунтовых в о д на вх оде и вы ходе расчетного участка (пьезометрический напор).

Задаются в виде линий из шести отрезков (рис. 16) по наблюдениям за уровнем грунтовых вод в течение года. W1 [тг] = W l n— дни замеров уровня, считая от начала года, 7г=0, 1 , ... , 6. Всегда W lo=0, a W l6 = 365; АН[п] =АНп и АН1[п]=АН1„— соответственно пьезометрические напоры в

эвапотранспирацию в формуле С. Ф, Аверьянова.

Задается в виде ломаной линии, состоящей из шести от­ резков (рис. 17) по данным наблюдений за испарением в те­ чение года. W2[n\ = W 2n— дни замеров испарения, начиная

пяти-шести суток после заданного времени £=7кон = 2кон. Выдача на печать промежуточных значений уровня грунто­ вых вод и концентрации солей по желанию может произво­ диться периодически через время 7=Шеч=^аВ,. Например, чтобы прогнозировать ВСР на 320 дней, начиная с 1 апреля,

182

и получать промежуточные прогнозы через каждые 30 дней,

— переключатель метода счета граничного условия урав­ нения 2 в зависимости от режима испарения или инфильт­ рации.

Рис. 19. Объем почвы, в котором рассчитывается изменение содержания солей.

Кроме того, программа вычисляет приращение количе­ ства солей — Sal [кг] на участке капиллярной каймы шири­ ной 1 м, длиной I и высотой Ай (рис. 19), к конечному моменту времени по формуле

S a l= Sytx—Syodv + SBy,

где

S y tx = m l ^y(x, A, t K0H)dv}

V

S yodv= mi J y ( x r A, tnil4)d v ; v

163