Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Трофимов, А. М. Основы аналитической теории развития склонов (на примере осыпных и делювиальных)

.pdf
Скачиваний:
7
Добавлен:
20.10.2023
Размер:
8.31 Mб
Скачать

V

График 39. Кривые, иллюстрирующие различ­ ный ход отступания крутой части склона во времени: равномерный, замедляющейся и коле­ бательный (см. график 38).

Кривая этой функции показана на графике 39 (линия Б). Находя производную от функции (1.3.5—8)

~ == ( К е Р° sin ® t y — bx e ~ p>t (іо cos (o t p 0 sin W)

запишем уравнение скорости перемещения точки перелома склона по уравнению (1.3.4—10) в виде

d y

_( \xhbx (1—х)

е х р [ - /» о О (ш c o s ü )( — Po sin м()2| 112. (1.3.5—9)

d t

(2(1— m)sinwf

 

Уравнение (1.3.5—9) уже с большей точностью может

описать скорость

развития осыпных склонов во времени.

Глава 4. ЛИНЕЙНАЯ И НЕЛИНЕЙНАЯ ТЕОРИИ

ИНТЕНСИВНОСТИ ДЕНУДАЦИИ СКЛОНА А. Е. ШАЙДЕГГЕРА. НЕКОТОРЫЕ ОБОБЩЕНИЯ И КРИТИКА

ТЕОРИИ. ПОДБОР КОЭФФИЦИЕНТОВ ДЛЯ НАИБОЛЕЕ РЕАЛЬНЫХ МОДЕЛЕЙ

Одной из наиболее завершенных теорий развития крутых склонов следует считать линейную и нелинейную теории ин­ тенсивности денудации А. Е. Шайдеггера (1964).

§1. ЛИНЕЙНАЯ ТЕОРИЯ ИНТЕНСИВНОСТИ ДЕНУДАЦИИ

Врассматриваемых выше моделях денудационных процес­ сов в результате отложения материала у основания склона формируется цоколи коренных пород. Подобные модели бы­ ли Даны Беккером и его последователями*. Теория Беккера,

однако, применима и к незащищенным осыпным материалом

* Обзор работ дан у Стралера (Strahler, 1952) и Шайдеггера А. Е (1964).

120

склонам. Использование основных факторов денудации скло­ на без учета отложения материала у его основания — основ­ ная йдея теории интенсивности денудации А. Е. Шайдеггера

(1964).

Учитывая, что модели .развития крутой части склона по Беккеру происходит пѵтем параллельного отступания склона и выполаживанием (радиальная схема развития склона), А. Е. Шайдеггер систематически исследовал, какую форму склона можно получить в результате применения названных гипотез интенсивности денудации. Он выдвинул предположе­ ние, что могут существовать три основных характерных слу­ чая развития склонов:

Случай 1. Процессы денудации происходят с одинаковой интенсивностью в каждой точке склона.

Случай 2. Интенсивность денудации пропорциональна вы­ соте расположения точки над некоторым базисным уровнем.

Случай 3. Интенсивность денудации пропорциональна кру­ тизне склона.

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

Обозначая высоту склона через у, длину через х, он вы­

водит

 

 

 

 

 

4 ^ = - const-Ф,

 

(1.4.1 — 1)

где Ф — интенсивность воздействия на склон со

случаями

Ф =

1; у; Ьу/Ьх.

(1 .4 .1 -2 );

(1.4.1 -

3; 1.4.1 - 4)

Очевидно,

всегда есть

возможность изменить время так,

чтобы постоянство (1.4.1 — 1) было равно

единице. Другими

словами, изменяя масштаб времени, мы можем произвести замену

^' = const-^,

(1.4.1 — 5)

и уравнение (1.4.1 — 1) приводится

к виду,

где const = 1 . Та­

ким образом, получается частное

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

ние.

Начальная форма склона представляет собой произволь­

ную функцию, входящую

в решение каждого

частного диф­

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

 

Случай 1 (гр. 40 А-1).

 

 

Полученное дифференциальное уравнение

 

~

= - 1

(1.4.1. — 6)

121

с начальным условием

 

У = /о (■*)

(1 .4 .1 -7 )

имеет решение

 

y = M x ) - t .

( 1 .4 .1 - 8)

Данное уравнение решается А. Е. Шайдеггером графи­ чески; отступание склона идет параллельно с равномерным понижением высот всех точек склона.

Случай 2 (гр. 40 А-2). Дифференциальное уравнение

^ = ~ У

(1.4.1 - 9 )

с начальным условием (1.4.1 — 7) имеет решение

J '= / o W - ехр[ — і].

(1 .4 .1 -10)

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

Случай 3 (гр. 40, А-3).

Дифференциальное уравнение

8_у_

»V

(1 .4 .1 -11)

Ы

8д:

 

с начальным условием (1.4.1— 7) имеет решение

Ѵ =/о { x - t ) .

(1 .4 .1 -12)

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

§ 2. НЕЛИНЕЙНАЯ ТЕОРИЯ ИНТЕНСИВНОСТИ ДЕНУДАЦИИ

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

В соответствии с этим и в отличие от случая линейной теории, А. Е. Шайдеггер дает другую, нелинейную формулу, более близко описывающую естественный процесс (гр. 40 Б).

T - - V 1 + ( і 7 ) ‘ -Ф

(1А.2— 1)

.122

График 40. Графическое выражение развития склонов по линейной (А) и нелинейной (Б) тео­ риям интенсивности денудации А. Е. Шайдеггера. А — модели склонов, показывающие характер развития с учетом коэффициентов а , ß, у в со­

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

рассматривая в качестве функции воздействия на склон

Ф (^; 9у/Ъх)

(1.4.2 — Іа)

те же частные виды (1.4.1 — 2; 1.4.1 — 3; 1.4.1—4).

Формула (1.4.2— 1) показывает,

что вертикальный эффект

воздействия (8у/№) интенсивности

денудации не просто про­

порционален величине самого действия (Ф) для любой точки, но зависит от наклона склона в каждой точке.

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

§ 3. КРИТИКА ТЕОРИЙ

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

123

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

Второй закон А. Е. Шайдеггера дает возможность учета одного существенного фактора интенсивности разрушения склона: высоту расположения рассматриваемой точки на склоне. При этом, однако, получается, что разрушение водо­ раздельной горизонтальной поверхности происходит быстрее, чем сам склон. Это также нереально. Этим законам верно отражается только то положение, что интенсивность разру­ шения склоновой поверхности при отсутствии подмыва тем больше, чем выше область на склоне (см., например, Спири­ донов, 1951; Маккавеев, 1955)*.

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

В нелинейных формулах законов А. Е. Шайдеггера эти резкие границы между тремя случаями несколько сглажива­ ются, так как всюду вводится коэффициент, отражающий крутизну склона.

Таким образом, приходим к выводу, что все три формулы А. Е. Шайдеггера, взятые по отдельности, дефектны. Наиболее существенный фактор учитывает третья зависимость, менее существенный — вторая, и почти нереальная — первая зависи­ мость.

Учитывая критику теории, наиболее правильным будет обобщение этих законов, где реальное значение каждого из природных факторов будет учтено величиной соответствую­ щего коэффициента. При этом, обобщенное уравнение будет описывать наиболее естественные формы развития склонов при условии равенства нулю коэффициента, соответствую­ щего первому случаю А. Е. Шайдеггера и когда второй коэффициент значительно меньше третьего.

* Существует, однако, область отсутствия эрозии близ водораздельной полосы (Хортон, 1948; Козьменко, 1949).

124

§ 4. ОБОБЩЕНИЕ ТЕОРИЙ ИНТЕНСИВНОСТИ ДЕНУДАЦИИ

Линейный случай

Исходя из вышеизложенного и полагая, что природный процесс включает в себя весь комплекс факторов, соотноше­ ние которых, правда, может меняться в силу тех или иных обстоятельств, целесообразно рассмотреть развитие склонов по данной теории в обобщенном виде (Митин, Трофимов, 1964).

Соответственно этому вводятся коэффициенты і, Р и і (для каждого из трех случаев). Введение коэффициентов не­ сколько упорядочит общую теорию, ибо учтет, что в общем развитии склона процент пропорциональности интенсивности

денудации высоте, к примеру, может

колебаться в зависи­

мости от колебания других соотношений.

 

С учетом последнего, обобщенная

формула будем иметь

вид

 

 

 

-âL. = a +

Р3, + Т -^ -.

(1.4.4 — 1)

Найдем решение уравнения

(1.4.4 — 1) при

начальных усло­

виях

 

 

 

Уы.о“ Я *)-

 

(1 .4 .4 -2 )

Уравнение (1.4.4— 1) является линейным неоднородным диф­ ференциальным уравнением в частных производных. Решение его распадается на два случая:

1. М О ,

2 . р = 0 .

1. ß =f=0 .

Введением неявной зависимости

{у, t, х) = const

(1.1.4 — 3)

сводим уравнение (1.4.4 — 1) к однородному. Но

в этом урав­

нении функция у является независимой переменной.

Искомое уравнение есть

 

(а + Ру) ср^ + <р, — те, = о.

(1.4.4 — 4)

Последнее легко решается с помощью сведения задачи к ре­ шению системы обыкновенных дифференциальных уравнений *, общий интеграл которых дает решение уравнения (1.4.4 —4) в виде

? = <р[* + • у ІП ( у +

— *] = const. (1.4.4—5)

* См., например: В. И. Смирнов (1954).

125

Отсюда

^ = ехр р - Ф ^ + ^ + Р^] —j ,

(1.4.4 6>

где Ф — произвольная функция от ( t + —^ .

Используя начальное условие (1.4.4 — 2), получим решение задачи Коши для уравнения (1.4.4 — 1) в виде

.У = [/(Т* + *)1ехр[ОД — -7 О — ехрІОД).

(1 .4.4-7)

р

 

2. ß = 0. Решается аналогично предыдущему. Общее ре­ шение записывается следующей формулой:

y = at + f ( x + it).

(1.4.4 —8)

Анализируя формулы (1.4.4 — 7; 1.4.4 —8) приходим к вы­ воду, что изменение формы склона сводится к замене, вопервых, аргумента на (л: + yt). Таким образом, график будет смещаться вправо или влево по оси х на величину |?| А Так как склон отступает в сторону водораздела, то естественно принять 7 с отрицательным знаком.

Во-вторых, множителем перед функцией f ( x + 7О является величина ехр[3/]. Очевидно, что высота склона не может увеличиваться в результате денудации, так что ß < 0 .

Второй член справа в формуле (1.4.4 — 7) должен быть отрицательным так как смещение профиля склона может происходить только в сторону уменьшения у. Отсюда сле­ дует принять а < 0 .

Итак, все коэффициенты а, ß, 7 являются отрицательными величинами. Заметим, что в случаях линейной теории они принимали значения (соответственно)

1)

- 1,

0 ,

0

2)

0 ,

- 1,

0

3)0 , 0 , - 1 .

Нелинейный случай

Так же, как и в линейном случае наше обобщение будет заключаться в том, чтобы рассмотреть взамен (1.4.2-- 1) уравнение

 

 

 

(1.4.4 --9)

 

 

)2-(а+Ру + 71 ь Г ).

Уравнение (1.4.4 — 9) является нелинейным

дифференциаль­

ным уравнением в частных производных.

 

Решим задачу Коши

для начальных условий (1.4.4 —2),

предварительно

(как и в

линейном случае)

разбив решение

на два случая:

ß =£0 и ß = 0 .

 

126 .

l . p ^

приступить к решению,

произведем

замену

Прежде чем

 

* =

Ѵ=а± . + у

 

 

 

 

Р

 

 

тогда уравнение

(1.4.4 —9)

перепишется

следующим

образом

% - Ѵ 1 + і т ; Т { ѵ + щ

<‘-4-4- 9а>

где X= a/ß.

 

При решении уравнения (1.4.4 —9а) наиболее удобно при­

менить

преобразование Лежандра

(Курант,

Гильберт, 1951).

і

 

 

9t» _,

 

9t»

 

 

 

 

 

Ь х

 

 

Эо»

 

 

 

X

8о»

t'

 

(1 .4 .4 -10)

 

 

«П ’

 

8т]

 

 

 

 

 

 

 

 

W + v = ix + rjt'.

 

В новых

переменных

уравнения

(1.4.4 —9а)

будет иметь вид

 

9о»

ч

9а» = w +

Ѵ \

X?.

(1 .4 .4 -11)

 

+

ц

 

+ 52

 

Последнее является линейным неоднородным уравнением в частных производных, метод решения которых был нами

использован

выше.

 

 

 

 

 

Общее решение (1.4,4— 11) дается формулой

 

® = $ / ^ _ ( X i ; - T ))ln& -7jln

(1 -К К Г + Т 2),

(1 .4 .4 -12)

где / — является производной функцией от (S/т»j).

и подста­

Из этой формулы методом

дифференцирования

новки можно

найти

 

 

 

 

 

x = f $ h ) + / Ф/ ч) —^inS( ь —

5 /

--------- - Kl + S2;

 

П

 

\

I + V 1 + £2

 

 

 

 

 

 

(1 .4 .4 -13)

У =

^ Г Ш

+ \ п 1 - \ п ( \

+ \/Т+1*);

(1 .4 .4 -1 4 )

 

т

 

 

 

 

 

® =

------—

3 ^ -—

г - - ( «

- ч).

(1 .4 .4 -15)

 

(1 + Ѵ і + «*) Ѵ Т Т І 2

 

 

 

Решение задачи Коши получим следующим образом: зная

начальное условие (1.4.4 —2)

и уравнение (1.4.4 — 9а), можно

найти D; bvßt-, b v ß x

в точке

t' = 0.

 

 

127

Таким образом, начальные условия

(1.4.4 — 2) при t' = О

запишутся

 

y = f o ( x ) ,

(1.4.4— 16а)

^ ~ / і ( * ) ,

(1.4.4 — 166)

" / üW ,

(1 .4 .4 - 16в)

где / 0; f x\ / 2 — получаются с помощью

(1.4.4 —2) дифферен­

цированием и подстановкой.

Из этих уравнений можно х, £, у выразить через функцию

от ЬҢ.

 

 

Принимая

в уравнения (1.4.4— 13;

1.4.4— 14; 1.4.4 — 15)

j' = 0;

<о = Л ( ^ ) ; W e

(to); -ч- h m

(последние четыре равенства получены с помощью системы уравнений (1.4.4 — 16)) можно найти функцию /, а с помощью ее и решение нашей задачи

 

 

 

 

 

 

 

У ) ,

 

 

 

где функция

<р получается исключением

переменных S и у из

уравнений

(1.4.4— 13; 1.4.4 — 14;

1.4.4— 15).

 

П р и м е р .

Возьмем

случай

начальной линейной

зависи­

мости

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у = ах;

t = 0.

 

 

 

Произведя

несложные вычисления, получим решение

 

У =

— In а + — + — ln ak

+ — ln

/

2é-exp [ßf]

\)

\

£2 — exp 12ßf]

Л

 

 

ß

 

ße

ß

 

ß

 

 

 

/

2&-exp[ßf]

\

_'a_

 

(1 .4 .4 -1 7 )

 

 

 

\ k * - e x p [ 2 ßf]

J

ß ’

 

 

 

 

 

 

 

где k = (1 4- V \

+ a?)ja.

 

 

 

 

 

 

Анализируя

данный пример, а также

 

формулу (1.4.4 — 7),

можно сказать,

что

величина —

показывает время

„жизни*

склона. Другими словами, тангенс угла наклона по истечении

этого времени уменьшается

по крайней мере в „ет раз.

2. ß = 0.

 

Уравнение (1.4.4 — 9) для

этого случая запишется

(1 .4 .4 -1 8 )

S - j / ‘ +(£)'(■ + ’ с)-

Последнее будем решать с помощью теории характеристик.

128

Характеристическая система для данного уравнения будет

иметь

вид

 

 

 

 

 

 

 

 

 

 

 

 

d x : d t : d y : d p : dq = ■

\р(л + т ] + (1 + />2)т ]: - 1 :

 

 

 

 

 

 

Ѵ і + Р г

 

 

 

 

 

 

 

 

 

 

■ — Р— - (Г + + 2тр -

д : 0:0,

 

(1.4.4— 19)

 

 

 

Ѵ і + р2

 

р

 

 

 

 

 

 

 

где р — Ьу/Ьх; q = Ь ѵ / Ы .

 

 

 

 

 

 

 

 

 

Решение

этой

системы представится так

 

 

 

 

 

 

 

 

 

Р =

Po't

 

 

 

 

 

(1.4.4 — 20а)

 

 

 

 

 

 

Я =

Яо\

 

 

 

 

 

(1 .4 .4 -206)

х =

- ■1_

 

0 (< 1 + Wo) +

(1 + />§)]t +

х0;

(1.4.4 - 2 0 в )

 

 

 

1^1 +

Р2

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

а - яЬт

+ У о-

 

 

 

(1.4.4 — 20е)

 

 

 

 

 

 

і

 

 

 

 

 

 

 

 

- Ѵ Т 7 7 о

 

 

 

 

 

 

 

Используя начальное условие

(1.4.4 —2),

зададим

началь­

ную полоску

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(0, *о) «=/(■*<>);

 

 

 

(1.4.4- 2 1 а )

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 .4 .4 -216)

 

 

 

 

 

 

Р о = Г ( Р о ) \

 

 

 

 

(1.4.4 -2 1 в )

 

 

 

Яо *=

К і

+

[ / '

( * о )] 2 - [*т+/

W

] ;

 

(1.4.4 — 21г)

 

 

 

 

 

 

/ =

о.

 

 

 

 

 

(1.4.4 -2 1 д )

Тогда

решение можно представить

системой

 

 

 

 

 

 

У=*

а \ Р

(-Ур)!3т

 

 

 

 

 

(1 .4 .4 -2 2 )

 

 

 

 

 

 

} /+ /( jc 0);

 

 

 

 

 

 

К і + [ / ' (*о)12

 

 

 

 

 

 

 

X * * —

 

 

1- ■-

{/' (JC0) la + f (*o)T l - d + l / '

С*Ь)]*Т)} t J r X 0.

 

Kl + [/'Uo)l2

 

 

 

 

 

 

 

 

(1 .4 .4 -23)

 

 

 

 

 

 

 

 

 

 

 

 

 

Для того, чтобы решение нашей задачи Коши было един­

ственным,

необходимо,

чтобы определитель

 

 

 

 

 

 

 

 

 

н _

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8t

8 х0

8 х0

8t

 

 

 

 

был равен

нулю (при у ф 0).

 

 

случае Д = — 1.

 

Легко подсчитать, что в нашем

и под­

Разрешая

уравнение

(1.4.4 — 23)

относительно х

ставляя

в уравнение

(1.4.4 — 22),

получим

искомое решение.

д-316.—9

129

Соседние файлы в папке книги из ГПНТБ