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

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

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

рости потока увеличиваются, захват материала становится все большим. В результате возможно предположить наличие момента, когда поток всю свою энергию будет тратить только на перенос частиц. Смыва здесь уже не происходит (область транзита). Ниже по склону, в связи с уменьшением уклонов и соответственно потерей скоростей, материал уже начинает выпадать в осадок. Здесь также возможно предпо­ ложить наличие точки склона, где мощность отложений до­ стигает максимума. Подобный процесс смыва по поверхности склона можно показать графически (гр. 61). Здесь знак плюс показывает величину сноса, знак минус — величину мощности отложений. Картина перемещения материала, изображенная на графике 61 фактически показывает нам характер изменения коэффициента денудации склона kd по функции kd = kd(x). Подобное распределение можно аппроксимировать синусои­ дальной функцией, типа (Трофимов, Переведенцев, 1969)

= а sin wx,

(2.4.3—1)

причем с достаточной точностью. Здесь: а — мощность (или глубина) вымытого материала в точке х с; ш = 2^/xc. Урав­

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

kd = a sin шл:-ср (л),

(2.4.3—2)

где: <?(х) — показывает интенсивность денудации

и аккумуля­

ции пропорционального уклонам склона.

склона; чем

Величина а зависит от среднего уклона

меньше средний уклон (а), тем меньше величина а и наобо­ рот. Характер функции a = jc(a) найден Р. Е. Хортоном (1948)

а — i4-sina/tg0,3a,

где А — коэффициент пропорциональности.

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

ния развития

склона

 

^

-

3 7 1 f (s) (sl" *■*•» <*> e~ “ 1Иг +

 

 

+

</<*) (s

<2.

где: / (s) — функция Хортона.

190

2. М етод реш ен ия уравнения развития делю виального склона

Для

простоты

решения

запишем

уравнение

(2.4.3—3)

в виде

 

 

 

 

 

 

 

— = — {а sin *хе~к‘} —

+ sin wxe~ Al}

(2.4.3—4) •

d t

д х 1

 

' д х

'

д х 2

 

Полагая

 

 

 

y = X (x)-T(t)

 

 

 

 

 

 

 

 

запишем

из (2.4.3—4)

 

 

 

или

X T ' — awe x/cosoxxX 'T + ae

keslnwxX"T

 

Т'

—\t

 

X'

х ш

 

 

 

 

 

e

 

— aw cos wx----- h a sin «ox — — c

 

 

T

 

 

X

 

X

 

где: штрихи обозначают производные, разделяющие перемен­ ные в уравнении (2.4.3—4),

г — const величина.

 

 

Г

- с 7£Гх' = 0;

(2.4.3—5)

sin «охX" + cocos «охX' — — ^ = 0.

(2.4.3—6)

Из (2.4.3—5) получим

 

 

dt

или Г = с0е х р [ - |- < Г х' ] .

(2.4.3—7)

- x t

Последнее (2.4.3—7) представляет собой рёшение уравнения

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

времени.

 

 

 

 

 

Для решения уравнения относительно х (2.4.3—7) сделаем

замену

 

 

 

 

 

 

v

у /

f ZdZ

 

 

Z = — ; X = b ]

;

 

 

 

 

X

 

 

 

 

 

X ' = Z X \ X " = Z 'X + Z X ' = Z 'X + Z 2X = X ( Z ' + Z 2)

 

 

 

 

 

 

(2.4.3—8)

Используя (2.4.3—8) запишем в окончательном виде

Z' + Z 2+ «о ctg «ох-Z •

a s i n

=

0 ,

(2.4.3—9)

 

 

 

а х

 

 

и получим уравнение Риккати,

которое

не

интегрируется

в квадратурах.

Чтобы определить общее решение, надо найти

три частных.

решение

уравнения

(2.4.3—4)

привело нас

Фактически

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

(2.4.3—3).

3. У равнение развития делю виального склона

Попробуем логически разобрать вариант построения урав­ нения. Начальная форма склона при / = 0 может быть запи­ сана функцией

У= <Р(•*),

укоторой (см. гр. 61) отметки в начале снижаются (процесс денудации), достигая максимума а в точке л:с , а затем повы­

шаются (аккумуляция); оба эти процесса обозначим симво­ лом ± о. На графике 61 пунктирной линией показано изме­

нение о по синусоидальному закону, т. е. o =

sinwx, а общее

уравнение у = tp(x) + .<о. В точке

наибольшего

сноса

о = + 1 ,

в точке наименьшего о =

— 1.

Первое

положение

(а = +

1)

показывает, что функцию

sin<o;c мы

можем

умножить

на

величину а, таким образом, показав истинную мощность сносимого материала в точке х с, т. е. a = asinu>c. Отсюда

получим, что мощность денудации в точке л;с равная а должна

в основании склона соответствовать мощности отложений, максимум которых также будет равен а, что вряд ли прием­ лемо. Очевидно, в природе этот процесс контролируется конкретными условиями склона, т. е. функция а sinwx: про­ порциональна функции <р(х), тогда о = а sin шхер (л:), а общее уравнение запишется

у — (х) + а sin «шр (х).

Рассмотрим теперь как изменяется о во времени. Исходя из теории затухающих геоморфологических процессов А. С. Девдариани (1967), можно показать, что а уменьшается по экспо­ ненциальному закону

 

« =

ао(1 - О .

где п — показатель степени.

сущности теории должна быть

Показателем

степени

по

величина, увеличивающаяся

во времени. Таким показателем

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

работой

А. М. Трофимова и Ю. П. Пере-

веденцева (1969) может быть функция

 

п =

1 — ехр[— X/],

которая при t = 0 дает ехр[— W] = 0; л = 1 .

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

у — ?(х) + а<р(х)\ 1 — exp [1 — е~и}} sin шх (2.4.3—10)

192

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

у = (f(x) + akl .

(2.4.3—Юа)

Подставляя вместо <?(х) ее значение,

например <р (х:) = h ё~(тх)

и полагая л: = 0, получим у = Л, т. е. наиболее высокую точку склона (водораздел). Далее, полагая / = 0, из (2.4.3—10) можно получить у = <р(л;) т. е. начальную конфигурацию склона. Оба эти положения показывают, что уравнение (2.4.3—10) спра­

ведливо

относительно наших граничных условий.

При t — 0

и х ф О

получим начальную конфигурацию

склона, т. е.

У = <?(х),

а при і ф 0 и х ф О — представление

о

характере

развития

конфигурации склона во времени.

 

 

4. Исследование уравнения развития делювиального склона

Проведем анализ уравнения (2.4.3—10)

у = <р (*) + ср (л;) а sin сох {1 — exp [1 — е~Хі)]}.

Здесь <р (л) — функция профиля склона, описывающая его кон­ фигурацию. Величина asinowc показывает величину удаляю­ щегося или аккумулирующего материала в каждой точке склона. Поскольку ш = 2тс/л;с, эта функция в точке л;с прини­

мает максимальное значение, в точке 2хс совпадает с функ­

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

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

Д-316.-13

193

§ 4 . О С Н О В Н Ы Е О С О Б Е Н Н О С Т И В Р А З В И Т И И Д Е Л Ю В И А Л Ь Н Ы Х С К Л О Н О В

Полевые наблюдения показывают, что положение области максимального сноса склона (д:с) зависит от конфигурации

склона, его склона, интенсивности сноса и общей физикогеографической обстановки. Изменение одного из факторов может привести к смещению ее вниз либо вверх по склону. Таким образом, положение области х с не является устойчи­

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

профиля динамического равновесия.

 

По

мере перехода склона

к равновесному состоянию все

более

четко

прослеживается

о с н о в н а я

тенденция — преиму­

щественное

перемещение области хс

в сторону вершины

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

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

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

особенностями ее перемещения.

Как показывают многочисленные исследования (Дедков, 1966; Гравис, 1969 и др.) современные процессы перестраи­ вают склоны, развитие которых происходило в отличных климатических условиях. Соотношение хода современной переработки (но не интенсивности процесса!) с переработкой, которую испытывали склоны ранее, можно проиллюстриро­ вать графиком 62 хода изменения уклонов склонов во вре­ мени. Строился график на основании ряда данных, в том

194

График 62. Ход изменения уклонов склонов во времени. Индек­ сами обозначены процессы: gr — обвально-осыпные; dp — оползне­ вые; sf — солифлюкционные: dl — делювиальные; p r— пролювиаль­ ные.

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

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

З А К Л Ю Ч Е Н И Е

Итак, использование математического моделирования в гео­ морфологии при изучении отдельных форм, сводится к сле­ дующим этапам:

к оценке относительной роли факторов, установлению связи между ними, нахождению вида этих связей, оценке временного фактора и т. п.

13*

195

График 63. Инерционно-геоморфо­ логический прогноз развития скло­ нов. Обозначения: А — основная тен­ денция повторяемости во времени конфигурации поперечного профиля крутого подмываемого склона; Б —

то же

для

осыпного склона; В — то

же

для

делювиального склона.

Использование, в свою очередь этих положений поз­ воляет перейти к изложению теоретических основ разви­ тия форм.

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

цессами

перемещения точки

перелома

склона

к бровке, а

у делювиальных

склонов —

с процессом перемещения об­ ласти максимального сноса (хс) к бровке. Установивший­

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

 

лювиального склона).

 

Найдено также, что особенность развития осыпных скло­

нов связывается с изменением

во времени

соотношения

вы­

соты склона и его уклона, а

особенность

развития делюви­

альных склонов связана с положением на склоне области

и

ее перемещения во времени.

 

 

 

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

По всей видимости вогнутый профиль склона отвечает понятию профиля динамического равновесия.

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

196

тых подрезаемых (А), осыпных (Б) и делювиальных (В) скло­ нов. Поскольку rt не уменьшается сразу до 0, он, стало быть, обладает прогностической информативностью. Попе­ речный профиль крутого подрезаемого склона так или иначе в конечном счете повторяет прежнюю конфигурацию в до­ вольно сходном виде (*4; г, = 0,98), хотя отклонения были значительны. Направленность в развитии осыпных и делюви­ альных склонов, заключающаяся в перестройке поперечного профиля, выражена более четко.

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

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

Л И Т Е Р А Т У Р А

А г е е в К. С., Д и т м а р А. В. Некоторые особенности рельефа высо­ когорского района Корякского нагорья.— Уч. зап. НИИГ Арктики. Регион,

геол., вып. 4,

1964.

 

 

Е. Д. Инструментальное изуче­

А н а н ь е в

Г. С., С ы р о е ч к о в с к а я

ние делювиальных

процессов

в Забайкалье.— Вопросы

географии, сб. 85,

Изд. „Мысль“, 1971.

Изучение

геоморфологических процессов эксперимен­

А р м а н д Д. Л.

тальным методом.— Тр. ин-та геогр., т. 47,

Материалы

по геом. и палео­

граф. СССР, вып. 4, М.—Л., 1950.

 

проектирования сети

А р м а н д

Д. Л. Физико-географические основы

полезащитных лесных полос. М., Изд-во АН СССР, 1961.

развитие берегов

А р ч и к о в

Е. И. Влияние склоновых процессов на

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

(Опыт зонального

исследования).

Автор, дисс. на соискание уч. степени

кандидата геогр.

наук. Владивосток, 1970.

 

П р е о б р а ж е н с к и й В. С.„

А р х и п о в Ю. Р., Б л а ж к о Н. И.,

С т у п и ш и н А. В., Т р о ф и м о в А. М.

Принципиальные вопросы исполь­

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

сы советской географической науки. М.,

.Наука“, 1972.

Б а б а н о в

Ю.

В., Т р о ф и м о в

А. М.

Возможные вариации развития

крутой части склона.—Геогр. сб., вып. 2, Казань, 1967а.

Б а б а н о в

Ю.

В., Т р о ф и м о в

А. М. Формы цоколей склонов при их

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

геологическом строении.— Вопросы

географии и геологии,

сб. IV, Казань, 1967 б.

А. М.

Динамика и морфология скло­

Б а б а н о в

Ю.

В., Т р о ф и м о в

нов (на примере Среднего Поволжья).— Сб.

Вопр.

геогр. и геол. Казань,

Изд. КГПИ, 1971.

 

 

 

Основы грунтоведения

Б а б к о в

В. Ф., Г е р б р у т - Г е й б о в и ч А. В.

и механики грунтов. М., Автотрансиздат, 1956.

руды. Металлургиздат,

Б а р о н

Л.

И.

Вторичное дробление и

выпуск

М., 1950.

Л. И.

Характеристика трения горных пород. М., .Наука*, 1967.

Б а р о н

Б а т ы р

В.

В.

Выступление по

докладу

Н. И. Соколова .О соотноше­

нии карста и явления отседания склонов“.— Тр. лаборат. гидро-геол. проб­ лем им. Саваренского. Общие вопросы карстоведения, М., 1961.

Б л а г о в о л и н Н. С., М у р а т о в

В.

М., Т и м о ф е е в

Д. А. Некото­

рые вопросы формирования склонов в условиях различных

морфоструктур.

Изв. АН СССР, сер. геогр., № 3, 1963.

Г.,

З а б о т и н Я.

И. Математико-

Б л а ж к о Н. И., Г р и г о р ь е в С.

географические методы исследования городских поселений. Казань, Изд-во КГУ, 1970.

Б ол и г А. Очерки по геоморфологии.

М., Изд-во ин. лит., 1956.

Б о р с у к

О. А.,

С и м о н о в

Ю. Г. О

механизме движения обломков

на солифл«опционных

склонах.— „Вести.,

научн. информ. Забайкальского

отд. географ, об-ва СССР“,

№ 3, Чита, 1965.

Б о р с у к

О. А.,

С и м

о н о в

Ю. Г.

Анализ кривых распределения

в геоморфологии.— „Вести, научн. информ. Забайкальского отдела ГО СССР",

№ 9, 1968.

198

Б о р с у к О. А., С и м о н о в Ю. Г. Статистические и вероятностные методы в геоморфологии.— Математические методы в географии. М., Изд-во

МГУ, 1968а.

 

В. П. Вопросы прогноза переработки берегов

водохрани­

 

Б о ч к а р е в

лищ Капчагайской ГЭС на р. Или.— Иссл. берегов водохр.

и

оз.

Байкал.

М.,

1964.

 

Ю. С. О причинах перемещения обломочного материала

 

Б у д и л и н

по склонам.— Тр. Центр, н.-и. горноразвед.

ин-та, вып. 56,

1963.

 

 

для

Б у т а к о в

Г.

П. Использование гранулометрических

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

определения

генезиса крупнообломочного

материала.— Вести,

студ.

научн. об-ва.

Естествен, науки, вып. 2, Казань, 1962.

 

 

 

 

 

Б у т а к о в

Г. П. Влияние размеров галек на

их форму.— Географиче­

ский сборник,

вып. 2, Казань, 1967.

 

 

сыпучих

осно­

 

Б у ц ь к о

3. Н. Об определении несущей способности

ваний.— Инж. сб.,

т. 26,

1958.

 

 

 

 

 

 

 

В е л и к а н о в

М. А. Моделирование руслового процесса.— ДАН СССР,

т. 70, № 3, 1950.

Л.

Проблемы формирования

рельефа

чащи

(берегов

 

В е н д р о в

С.

н дна) больших водохранилищ. Изд-во АН СССР, 1959.

 

 

 

 

 

В с т е л и у с

А. Б. Морфометрия обломочных частиц.— Тр. лаборатории

•аэрофотометодов, т. 9, 1960.

 

Саратовского

ун-та,

1967.

 

Вопросы морфометрии (сб.). Вып. 2. Изд-во

 

В о с к р е с е н с к и й

С. С. Изучение

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

эндогенных и

экзогенных сил в историческом аспекте — основная задача геоморфологии.— Вести. Моек, ун-та. География, № 2, 1968.

В о с к р е с е н с к и й

С.

С. Склоны, их формирование и строение.—

Вестник МГУ, сер. геогр. № 3, 1969.

В о с к р е с е н с к и й

С.

С. Динамическая геоморфология. Формирование

склонов. Изд-во МГУ, 1971.

С., З о р и н Л. В., С и м о н о в Ю. Г. Законо­

В о с к р е с е н с к и й

С.

мерности формирования склонов в Восточной Сибири.— Вестник МГУ, сер.

геогр., 5, № 5, 1960.

 

С. С., Т и м о ф е е в

Д. А.

Формирование скло­

В о с к р е с е н с к и й

нов.— Современные экзогенные процессы рельефообразования. М.,

.Наука",

1970.

 

Н. А. Тектонические трещины и сила бортового отпо­

Г в о з д е ц к и й

ра.— Вести. Моек, ун-та. География, № \ t

1966.

 

 

 

 

 

Г е р а с и м о в

И. П. О движении почвенно-грунтовых масс на скло­

нах.— „Почвоведение", № 7/8, 1941.

 

 

 

 

 

интерпретации общей

Г е р а с и м о в

И. П. Опыт геоморфологической

схемы геологического

строения СССР.— „Проблемы

геоморфологии", XII,

1946.

 

И. П. Конструктивная география; цели, методы, резуль­

Г е р а с и м о в

таты. Изв. ВГО, т. 98, вып. 5, 1966.

К.

К.

 

Четвертичная

геология. М.,

Г е р а с и м о в

И.

П.,

М а р к о в

 

Учпедгиз, 1939.

Н. Ф., Л у к а ш о в А.

А.,

С и м о н о в

Ю. Г. Некото­

Г л а з о в с к и й

рые особенности размещения крурумов на

ЮВ

Читинской

области.— Зап.

Забайк. отд. геогр. об-ва СССР, вып. 21, 1963.

Устойчивость откосов при

Г о в я д и н о в

А.

И.,

Фа ль ко в ич

С.

 

В.

предельном состоянии равновесия.— Инж. сб.,

т.

14, 1952.

 

накопления

Г р а в и с

Г. Ф. Склоновые отложения

Якутии

(условия

и промерзания. Криогенное строение). М., „Наука", 1969.

 

 

 

Г р а д ш т е й н

И.

С.,

Р ы ж и к И.

М.

Таблицы интегралов сумм, рядов

и произведений. М., 1962.

 

 

 

 

 

 

 

 

 

 

Г р а ч е в

А. Ф.

К теории педипленизации. Мат. по геоморфологии и

неотектонике Урала и Поволжья. Вып. 1. Уфа, 1962.

размыва

берегов на

Г р е ч и щ е в К. К. Метод расчета

ширины

зоны

примере Братского водохранилища. Иркутск, 1961.

и плотность

вероят­

Г у р е в и ч

Б. Л.

Плотность населения

города

199

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