Файл: Трофимов, А. М. Основы аналитической теории развития склонов (на примере осыпных и делювиальных).pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 122

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

лингом

(Culling, 1963)

и

Суше (Souchez, 1964) уравнению

развития

склонов

 

 

 

 

Ь

. и

Ѵу

(2.4.1-5)

 

b t

 

»x* ’

 

 

 

представляющему собой уравнение теплопроводимости с по­ стоянным коэффициентом.

А. Е. Шайдеггер (1964), детально анализирующий про­ цессы переработки склона, указывал, что „в первом при­ ближении процесс формирования склона определяется зави­ симостью между скоростью течения и транспортирующей

способностью потока“

[стр. 119].

Смысл

этой зависимости

можно показать на двух примерах:

а) Если плоскостной или

ручейковый поток несет наносов

меньше, чем

позволяет

его транспортирующая

способность, то

он будет

добирать

его и тем самым эродирует склон; б) Если способность потока переносить материал уменьшается, то последний выпадает в осадок. Исходя из этих двух примеров, А. Е. Шай­ деггер (1964) делает вывод, что „существует динамическое равновесие между потоком и материалом, слагающим его ложе. При отклонении от этого равновёсия склон эроди­ руется или становится ареной аккумуляции осадков“ [стр. 120]. Отсюда заключаем, что в каждой точке склона процесс деформации различен и определяется соотношением факто­ ров перемещения материала. Интенсивность выветривания материала определяется соотношением величины прихода

материала и

его

удаления

в каждой

точке

склона.

В вер­

шинных частях склонов

поступление

материала извне менее

значительно, а отсюда и большая

интенсивность выветри­

вания. Обратная

картина

наблюдается в областях

аккуму­

ляции

склона.

 

 

 

 

 

показывает, что

Таким образом, сама сущность процесса

коэффициент k

изменяется

одновременно

в зависимости и

от изменений х и от изменений t,

и,

следовательно,

не мо­

жет

быть

величиной

постоянной.

 

В соответствии

с этим

форма профиля склона наиболее точно может быть записана также уравнением теплопроводности, однако в отличие от (2.4.1—5) оно сильно изменится (Трофимов, Переведенцев, 1969)

{kd{ x , t ) \ ^ - + {kd( x , t ) \ - ^

(2.4.1-6)

»X

 

 

 

В этом случае коэффициент

kd, которому может

быть при­

писан смысл коэффициента

денудации склона, является

не

постоянной величиной, а зависит одновременно от х и t, т.

е.

kd {x,t) = kd{x)-kd{t)

(2.4.1—7)

184


где

 

kd(t) = e~u .

(2.4.1- 8а)

Последняя запись (2.4.1—8а) наиболее очевидна, ибо, как было нами показано ранее, интенсивность развития склонов во времени носит замедленный характер. Из (2.4.1—5) за­ ключаем, что в уравнении (2.4.1—6) неизвестным остается зависимость kd(;x).

Коэффициент денудации склона kd определяется,

вообще

говоря, множеством факторов.

Однако

наиболее значимым

из них можно

считать

эрозию: степень

или

интенсивность

ее воздействия. Как показали исследования

Р. Е. Хотона

(1948), интенсивность эрозии различна,

во-первых,

на

скло­

нах различной крутизны и, во-вторых, в

каждой

точке

склона

(т.

е.

зависимость

нелинейная) и записывается в виде

функции

 

 

X (s) =

sin <*/tg0" \

 

 

 

(2.4.1 —9)

 

 

 

 

 

 

где x(s)

интенсивность

воздействия

эрозии

на

склон

в за­

висимости

от

уклона (а)

в каждой

точке склона

(гр. 44).

Таким образом, согласно приведенной зависимости,

которая

получила

название функции

Хортона,

эффект

денудации

в каждой точке склона не представляет собой линейной зависимости до незначительных по величине уклонов (Fardley, Vivant, 1967). Распределение же уклонов по склону

выражается в виде параболической

кривой (т. е. уклоны

от минимальных у основания склона

переходят к более зна­

чительным и выше бровки вновь уменьшаются); причем мак­ симальные уклоны делювиальных склонов не превосходят значения критических по схеме Хортона (1948). Учитывая это положение, распределение интенсивности эрозии может быть подчинено распределению уклонов и, следовательно, выразится кривой, повторяющей кривую распределения

уклонов (гр. 44, 45, 46).

 

что

 

 

Из сказанного полагаем,

 

 

 

 

kd(л;) =

ах2 + Ьх +

с,

 

где а, Ь, с — коэффициенты уравнения,

и на

основании этого,

переписав

(2.4.1—7),

в виде

 

 

 

 

kd(x,t) = (ах2 + bx + с) е~и

(2.4.1-10)

запишем

окончательный

вид дифференциального уравнения

(2 .4 .1 -6)

 

 

 

 

 

 

 

— = ~ {(ад:2 + Ьх + с) е~х<) — +

 

öf

дх

'

 

 

' » X

 

+ {(ах2 +

Ьх + с ) е ' ХІ\ ^ .

(2.4.1-11)

185


Для того, чтобы решение (2.4.1 — 11) было однозначным, поставим начальное и граничное условие

Уi t - о) = ?(■*).

т. е. первичная конфигурация склона может быть принята произвольной:

У = ч(х),

(2.4.1-12)

УІХ- О) = ? і(0 ,

(2.4.1-13)

Уі, - ,) = ? а ( 0 .

(2.4.1-14)

т. е. в каждый последующий момент времени конфигурация

склона изменяется

в соответствии с уравнением

(2.4.1 — 12).

Найдем решение

уравнения

(2.4.1—12) методом Фурье.

Записав

у = ХТ,

 

 

(2.4.1-15)

 

 

 

где Х = Х(х); T — T(t), и подставив

(2.4.1 — 15)

в уравнение

(2.4.1—11), получим

 

 

 

 

X T ' = <ГХГ(2ах + Ь )Х 'Т + (ах2 +

Ьх +

с) е~ХіТХ",

(2 .4 .1 -15a)

где штрихи обозначают производные, разделяющие пере­

менные в уравнении

(2.4.1—11).

 

 

 

 

 

Наконец, запишем

 

 

 

 

 

 

 

 

 

 

 

Г

+ те~Хі Г = 0,

 

 

(2.4.1-16)

(ах2 + Ьх +с) X" +

(2ах +

Ь) X' +

т Х = 0,

(2.4.1-17)

где т — параметр

разделения,

определяющийся

из

гранич­

ных условий.

 

 

 

(2.4.1—17) в общем случае

известно

Решение уравнения

(Камке, 1965). Подстановкой

 

 

 

 

 

 

 

Х = уі (5);

X = p + ( q - p ) t ,

 

(2.4.1-18)

где р, <? — неравные

корни уравнения

ах2 + Ьх + с = 0 ,

сво­

дим уравнение к гипергеометрическому виду

 

 

 

5(5— 1 )Ѵ '+

Г25+

2ар + Ь-1 V + —

= 0.

(2.4.1-19)

 

L

 

 

a {q

р ) J

а

 

 

 

 

Стандартное же

гипергеометрическое

уравнение

имеет вид

х (х — 1 )у" + [(« + ß + 1)х — і \ у г + а$у = 0

(2.4.1—20)

решение которого

дается

в

виде

сходящегося

ряда

при

| х |< 1, а в частном

случае

и при х — 1, когда т — а — ß > О

/=Ч*. Р. Т,

і +

^

+ а (я + 1) Р(3 + 1) r 2 I

(2.4.1-31)

186

 

 

1*7

1-2.Т(Т+1)

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 


Сравнивая коэффициенты в уравнениях (2.4.1—19) и (2.4.1—20), можно записать

а == 2; р.

1.

т

2ар + Ь .

- L = ap = - 2 (2.4.1-22)

* *

і

/

V*

 

 

 

a(q - р )

а

откуда можно заключить, что параметр разделения т равен

 

 

 

 

т

2 а .

 

 

Для нашего случая решение (2.4.1—19) будет

иметь вид

 

 

 

F = 1 +

-M frzig L I

 

(2.4.1-23)

 

 

 

 

2ap + b

 

 

 

а далее ряд обрывается, так

как ß = — 1.

получим

Заменяя $,

А, х истинными значениями,

 

 

F = 1 +■

- - ( х - р ) .

 

(2.4.1-24)

 

 

 

 

2ар + b

 

 

 

Уравнение (2.4.1—16) имеет

решение

 

 

 

 

Г (0 = ^ , е х р ^ е - х/] .

 

(2.4.1-25)

где т — 2 а ;

с, — начальное

условие.

примет вид

И, наконец, общее решение (2.4.1—11)

у ( х , 0 - [ і

+

^ _ г ( л - Л ] с , е х р [ - - ^ - в - м] . (2.4.1-26)

С помощью

начального

условия

(2.4.1—12) определим зна­

чение с х в уравнении (2.4.1—26) в виде

 

 

 

 

 

_

 

< ? ( х ) е 2аІХ

 

 

 

 

 

сі

 

2а

( х — р )

 

 

 

 

 

 

1 + 2ар + b

 

 

 

 

 

 

 

 

 

и подставляя

его

в уравнение (2.4.1—17)

запишем

у ( х ,

0 =

? ( * ) е х р [ - ^ ( 1 - < Г х') ] ,

 

(2.4.1-27)

где <р(х) — начальная форма склона.

 

характеризует

Учитывая,

что квадратическое

уравнение

параболу, обращенную выпуклостью вверх, следует при­ писать коэффициенту а отрицательный знак. Итак, (2.4.1—27)

перепишется в окончательном

виде

 

 

 

 

 

у ( X,

t ) =

cp ( X ) ехр[—

(1 -

e~xt) ] .

(2.4.1-28)

Исследуя уравнение (2.4.1—28),

можно

привести два гра­

ничных

варианта,

когда

( = 0

и

t

оо.

В первом случае,

полагая

t=* 0

из

(2.4.1—28)

получим

у

=ф (х) — начальное

условие

(как и следовало

полагать).

Во втором случае, при

t — о о , т. е. время,

в течение которого

склон

превращается

187


разрушаясь в слабонаклонную поверхность. Из (2.4.1—28) получаем

у (t = оо) = <р (л:) ехр [— 2а/Ц.

§ 2. ИССЛЕДОВАНИЕ УРАВНЕНИЯ РАЗВИТИЯ ДЕЛЮВИАЛЬНОГО СКЛОНА

Рассмотрим смысл каждого члена уравнения (2.4.1—28)

У (х, 0 = ? (■*) ехр [(—2aß) (1 — e~lt)\.

Неявная функция ср(л:)— есть начальная форма склона, кото­ рая в соответствии с (2.1.1—5) может быть выражена урав­ нением

y==he~(kx)\

(2.4.2—1)

Функцией (2.4.2—1) можно характеризовать профили склонов, имеющие три основные сегмента: выпуклый в верхней части» прямой в средней и вогнутый в нижней.

Член уравнения (2.4.1—28)—(—2aß) определяет интенсив­ ность денудации склона в каждой точке. Поскольку со вре­ менем (при t —»п) степень (—№), под которой стоят основа­ ния натуральных логарифмов, увеличивается. В связи с этим

уменьшается выражение е~и. При t —>оо, lim (1 — ё~и) = 1. Из последнего следует, что деформация профиля во времени определяется величиной (— 2 aß). В целом, из степени

[(— 2а/Х) (1 — следует, что в конечном счете величина интенсивности денудации связана с (—2aß). Коэффициент (X), стоящий перед t, определяет интенсивность денудации во времени, из выражения (—2а/Х) следует тот же вывод, однако из первого полагаем, что он в прямой зависимости от t и в обратной от а. Впрочем, это положение справедливо только при Х > 1. При Х <1 следует обратное соотношение.

Из степени [(— 2а/Х)(\ é~lt)\ учитывая сказанное, заклю­ чаем, что интенсивность денудации склона определяется коэффициентом а, которому может быть придан смысл коэф­ фициента денудации, как и в исходном уравнении (2.4.1 — 10) коэффициенту kd, однако после замены его через х и t функ­

ций kä(x, t) (2.4.1 — 11) его сущность выражается через а. Из уравнения (2.4.1—28), особенно из анализа степени

[(—2aß) (1е~х‘)] можно заключить, что все точки профиля <р(х) со временем (t —*n) понижаются в соответствии с вели­ чиной а и пропорционально коэффициенту >. во времени. Отсюда можно сделать вывод, что если склоны развиваются по механизму, описанному моделью (2.4.1—28), то они в своем развитии выполаживаются. Параллельного отступания какихлибо точек здесь нет.

188