Scientific journal

Modern problems of science and education

ISSN 2070-7428

"Перечень" ВАК

ИФ РИНЦ = 1,006

Тенденция математизации характерна для большинства научных направлений исследований. Применение математического аппарата уже не ограничивается использованием статистических методов для обработки экспериментальных данных, все чаще акцент делается на математическое моделирование. Математическое моделирование эрозионно-аккумулятивных процессов почв - относительно молодое научное направление, основанное на методах моделирования сложных динамических систем, которые в значительной степени зависят от принципа иерархической организации или принципа интегративных уровней [9]. Этот принцип утверждает, что для предсказания поведения сложной системы не обязательно точно знать, как ее компоненты построены из более простых компонентов. В зависимости от степени упрощения для одной и той же системы-оригинала можно получить несколько различных моделей. Степень детализации модели, форма ее представления в первую очередь определяются целями исследования.

Областью наших исследований является эрозионный процесс почв на склонах горных и предгорных ландшафтов, который обычно делят на два этапа [6]: площадной и линейный. Первый этап связан с поверхностным стокообразованием, эрозионным размывом почв, слагающих поверхность водосбора, и поступлением в русловую сеть наносов. В качестве объекта нашего исследования была выбрана склоновая (плоскостная) водная  эрозия [6, 7, 8]. Активность эрозии определяется различными факторами, такими как: интенсивность поверхностного стокообразования, податливость почвы или рыхлой горной породы размыву, характеристики рельефа и подстилающей поверхности. Интенсивная эрозия в исследуемых районах наблюдается в период ливней и паводков [3, 4].

Положительная сторона методов, основанных на математическом моделировании данных по эрозионно-аккумулятивному процессу, состоит в том, что при их использовании:

• есть возможность получения предельно точной информации о характеристиках смыва (это направление потока веществ, места их отложения, скорость смыва и т.д.);
• ускоряется анализ и доступ к данным большого объема;
• появляется возможность точного прогноза развития почвенной эрозии.

Математическое моделирование эрозионно-аккумулятивных процессов насчитывает десятки моделей, представленных в виде алгебраических или дифференциальных уравнений, основанных на различных подходах и опирающихся на общие и частные физические законы. Можно выделить три типа моделей для прогнозирования эрозии: физические, эмпирические и полуэмпирические. Так, значительный объем исследований в данной области основан на теории взвесенесущих потоков (Маккавеев В.М., Гончаров В.Н., Мирцхулава Ц.Е., Швебс Г.И., Кузнецов М.С. и др.) [1, 12], которая составляет основу для различных полуэмпирических зависимостей поверхностной эрозии почв.

Существующие методики количественной оценки смыва почвы можно свести к следующим:

-        методы балльной оценки;

-        сравнительные методы;

-        эмпирические зависимости;

-        методы, основанные на физике гидромеханического процесса.

Суть метода балльной оценки сводится к количественной оценке качественных явлений. Каждому из существенных с точки зрения изучения моделируемого процесса факторов назначается балл. В некоторых моделях предполагается возможность ранжировать основополагающие факторы с указанием их весового значения. Интегральная оценка смыва почвы представляет собой линейную свертку назначенных баллов. В основу данной методики может быть положен математический аппарат нечеткой логики, предусматривающий возможность количественной оценки качественных показателей, которым присуща субъективность [11].

Метод сравнительной оценки - это разновидность метода балльной оценки. Особенностью данной методики является необходимость сравнения количественных данных, полученных в результате вычислительного эксперимента, с эталоном. В результате имеется объективное обоснование полученных данных и выводов. Основной недостаток данного метода заключается в разрозненности оценок взаимосвязанных величин.

Одним из более совершенных приемов считается метод эмпирической оценки смыва, в основу которого положен аналитический аппарат обобщения количественных законов изучаемого процесса. В качестве проблем данного подхода Швебс Г.И. указывает отсутствие надежного обоснования структуры расчетной схемы и формально-статистический подход к определению многих критериев расчетной формулы.

Гидромеханический метод в силу возможности раскрытия механизма явления имеет ряд преимуществ, но зачастую это делается без должного изучения динамики процесса. В данном случае необходимо детально рассматривать все исходные условия, что налагает определенные ограничения на механический перенос уже известных выражений на изучаемые процессы и явления.

В настоящее время имеется довольно обширный ряд методик, позволяющих прогнозировать эрозию, а также большое число противоэрозионных мероприятий. Однако  в научной литературе не приводится описание обобщенных методик, связующих их в единую последовательность логических или расчетных процедур. 

Целью нашего исследования стала разработка обобщенной математической модели эрозионного процесса склоновой подсистемы горного и предгорного ландшафта, алгоритм которой приводится на рисунке.

Алгоритм построения обобщённой математической модели эрозионного процесса склоновой подсистемы горного и предгорного ландшафта

Одним из самых распространенных методов определения почвенных потерь является метод, основанный на универсальном уравнении потери почвы [13].

,

где  - потеря почвы в кг/м2,  - фактор эродирующей способности дождя,  - фактор подверженности почв эрозии,  - фактор длины склона,  - фактор уклона,  - фактор системы ведения растениеводства,  - фактор борьбы с эрозией.

Факторы длины  и крутизны  склона часто рассматриваются как единый топографический фактор  или, иначе именуемый, «фактор рельефа». Для склонов со сложным профилем Уишмейером было предложено упрощенное решение, предусматривающее разбиение сложного склона на несколько сегментов, каждый из которых должен иметь ровный уклон и однородный почвенный покров. Значение фактора  находится путем деления суммы значений каждого сегмента на полную длину склона.

Данная методика легла в основу действующего государственного стандарта ГОСТ 17.4.4.03-86, введенного в 01.07.1987 г. Настоящий стандарт устанавливает метод определения потенциальной опасности эрозии почв под воздействием дождей с разбивкой по классам интенсивности для составления карт эрозионной опасности. Метод основан на определении факторов атмосферных осадков, устойчивости почв, длины и крутизны склона, севооборота, агротехники и зависимостей между ними [2].

Модели, в основе которых лежит «универсальное уравнение потерь почвы», обладают большими возможностями в области прогнозирования и планирования противоэрозионных мероприятий, а также довольно просты в применении на практике. В силу этого они получили большое распространение во многих странах и легли в основу исследований, направленных на разработку региональных вариантов этой модели.

В общем виде влияние рельефа на величину смыва почвы со склонов может быть выражено зависимостью [10]:

,

где  - среднее значение модуля смыва для склона длиною L, а, n и p - эмпирические показатели. По результатам экспериментальных данных было установлено, что величина p  колеблется в диапазоне [0,4; 0,5], а n - в диапазоне [0,7; 2,2]. Величина а получена обратным расчетом для набора опытов и осреднена [10].

Экспериментальные исследования и обобщение обширных материалов дают общее выражение для среднего модуля смыва [10]:

.

Показатель степени при уклоне зависит от эрозионной устойчивости подстилающей поверхности и приводится в таблице.

Зависимость показателя степени крутизны склона от эрозионной устойчивости  поверхности

Для осредненных условий, водосбора прямоугольной формы  и прямого профиля была получена следующая зависимость:

 [10].

Широко применяются модели смыва почв, где в качестве основных  факторов эрозионного процесса учитывают длину и крутизну склона либо используют их как интегральный показатель - «фактор рельефа» (Сурмач Г.П., Мирцхулава Ц.Е., Швебс Г.И., Ларионов Г.А. и др.). По заключению Гаршинева Е.А., фактор рельефа  (I - крутизна, L - длина склона) является значимым фактором, объединяющим в себе гидрологические и топографические характеристики.

Реализация предложенного выше алгоритма является начальным этапом в построении  системы оценки экологической стабильности склоновой подсистемы на основе полученной обобщенной математической модели эрозионного процесса склоновой подсистемы.

Анахаев К.Н., д.т.н., профессор, зам. директора по селевой проблематике ФГБУ «Высокогорный геофизический институт», г. Нальчик;

Ламердонов З.Г., д.т.н., профессор, Кабардино-Балкарский аграрный университет,                   г. Нальчик.