химических процессов, в математической модели эти процессы идеализируются, то есть учитываются наиболее существенные из них, связи между ними © отвлечением от несущественных процессов и связей. Предлагаем — схему — моделирования (В виде алгоритма) со следующим оптимальным порядком выполнения моделирующих операций: 1) расчет среднегодо- вой скорости испарения; 2) расчет ионного состава посту- пающих на поверхность почвогрунта вод; 3) расчет ион- ного состава гипотетических растворов в подошве зоны транспирации; 4) расчет ионного состава гипотетических растворов в горизонте максимального соленакопления зоны транспирации; 5) учет явления отрицательной адсорбции; 6) учет явления положительной адсорбции; 7) осаждение из раствора избытка СаСО,; 8) осаждение избытка МеСО,; 9) осаждение избытка СаЗО,; 10) приведение в равновесное состояние катионообменной системы Ма* * [Са?*] == Са?* == + [No']; 1!) приведение в равновесное состояние катионооб- — — менной системы К*=[Ме**] => Ме?* == [К*]; 12) оценка явле- — р ННУННННН } ний содового засоления, осолонцевания и суммы растворимых и токсичных солей в поровом растворе. При составлении модели учитывалось, что, представляя реальные процессы с достаточной точностью, она должна быть достаточно простой. Каждая операция представляет со- бой численный алгоритм, который сводит описание сложного процесса к арифметическим и алгебраическим действиям над числами. Последовательное выполнение всех операций в пред- ложенном порядке приводит к требуемому результату; при этом точность конечного результата зависит лишь от точности промежуточных вычислений. Ниже приведены численные ал- горитмы для каждой операции; концентрации ионов Во всех уравнениях выражены в киломолях на кубометр раствора. Операция 1. _ С, ‚1__—————">хг>0; оо ехр —- = 0,, шаг == 0,05 к где л;— прогнозируемая глубина залегания уровня иррига- ционно-грунтовых вод, м; @ — находится по формуле (10.9). Операция 2. 10 > 8 > 0, д , * * та —0г° 0С* = Со; Со == С, при 9 << 1; при @0>1 из раствора следует осадить избыток СаСО, и СабО,: 0,4888 + 0,00087* == А*; 0,3241 - 0,000177* = В*; и* < [НСО, ; — 0,001) при [НСО, ; < 2 [Са*, , 241