Научная статья на тему 'Оценка изменчивости элементов водного баланса половодья и межени в бассейне реки Волги'

Оценка изменчивости элементов водного баланса половодья и межени в бассейне реки Волги Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
213
51
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭЛЕМЕНТЫ ВОДНОГО БАЛАНСА / ELEMENTS OF WATER BALANCE / БАССЕЙН РЕКИ ВОЛГИ / RIVER BASIN OF VOLGA / МЕЖГОДОВАЯ ИЗМЕНЧИВОСТЬ / VARIABILITY BETWEEN SEVERAL YEARS / ПОЛОВОДЬЕ / HIGH WATER / МЕЖЕНЬ / LOW WATER / КЛИМАТ / CLIMATE / АВТОКОРРЕЛЯЦИОННАЯ ФУНКЦИЯ / AUTOCORRELATION FUNCTION

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Исмайылов Габил Худуш Оглы, Муращенкова Наталья Владимировна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Исмайылов Габил Худуш Оглы, Муращенкова Наталья Владимировна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Assessment of the changeability of water balance elements of high water and low water in the Volga river basin

Results of the analysis and estimation of variability of elements of water balance for the periods of a spring high water and low water on 11 private reservoirs and as a whole on a river basin of Volga are resulted. In a river basin of Volga four water-economic areas on conditions of uniformity of the basic statistical parameters of elements of water balance, and also their communications inside of some and between numbers are allocated. Dynamics of variability of elements of water balance of the river of Volga in space and in time is considered. The equations of communication defining values of a river drain for the periods of a high water and low water from climatic factors defining them are received.

Текст научной работы на тему «Оценка изменчивости элементов водного баланса половодья и межени в бассейне реки Волги»

Гидравлика, гидрология, водные ресурсы

УДК 502/504 : 556.55

Г. Х. ИСМАЙЫЛОВ, Н. В. МУРАЩЕНКОВА

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Московский государственный университет природообустройства»

ОЦЕНКА ИЗМЕНЧИВОСТИ ЭЛЕМЕНТОВ ВОДНОГО БАЛАНСА ПОЛОВОДЬЯ И МЕЖЕНИ БАССЕЙНА РЕКИ ВОЛГИ*

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

Элементы водного баланса, бассейн реки Волги, межгодовая изменчивость,

половодье, межень, климат, автокорреляционная функция.

Results of the analysis and estimation of variability of elements of water balance for the periods of a spring high water and low water on 11 private reservoirs and as a whole on a river basin of Volga are resulted. In a river basin of Volga four water-economic areas on conditions of uniformity of the basic statistical parameters of elements of water balance, and also their communications inside of some and between numbers are allocated. Dynamics of variability of elements of water balance of the river of Volga in space and in time is considered. The equations of communication defining values of a river drain for the periods of a high water and low water from climatic factors defining them are received.

Elements of water balance, river basin of Volga, variability between several years,

high water, low water, a climate, autocorrelation function.

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

условиям формирования и амплитуде колебаний ЭВБ зависит во многом от внутризональной дифференциации, обусловленной азональными факторами и ландшафт-

ландшафтных комплексов определяет общий характер влагооборота в бассейне.

Пространственную изменчивость ЭВБ можно рассматривать в разных аспектах:

как изменчивость отдельных характери-

ной структурой бассейна. Речной бассейн __________ „„„„„„„ „„,„„,„

__ _ стик элементов водного баланса (их стати-

* Работа выполнена при финансовой поддержке стических параметров, годовых, сезонных, РФФИ (проект № 12-05-00193а). месячных значений и тд.) по территории

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

Различное сочетание климатических и неклиматических факторов, включая и антропогенные, и определяет многообразие реакции речных водосборов на одни и те же осадки. Поэтому при рассмотрении пространственной изменчивости элементов водного баланса в первую очередь выделяются районы, однородные по структуре водного баланса (ВБ) и многолетней динамике его элементов. Формальной основой такого районирования может служить матрица взаимной корреляции временных рядов ЭВБ частных водосборов, принятых к рассмотрению.

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

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

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

Для выявления пространственно-временных закономерностей изменения ЭВБ речных бассейнов авторами получены временные ряды за период с 1914/1915 по 2000/2001 годы N = 87 лет) основных элементов водного баланса (речной сток, атмосферные осадки, суммарное испарение и бассейновые влагозапасы) по 11 частным водосборам бассейна реки Волги за период половодья, межени и года в целом [1, 2]. Наличие столь длительных временных рядов позволило осуществить оценку изменения основных характеристик ЭВБ бассейна реки Волги как в пространстве, так и во времени. Дополнительно к традиционно рассматриваемым элементам водного баланса были определены также разности «осадки - сток», «осадки - испарение». Кроме того, для оценки изменчивости гидротермического режима территории бассейна реки Волги авторы рассмотрели также изменение температуры воздуха в указанных 11 частных водосборах бассейна реки Волги, включая бассейн в целом.

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

Таблица 1

Выборочные оценки основных статистических параметров временных рядов годовых и сезонных величин ЭВБ зоны формирования стока в бассейне реки Волги за 1914/1915 - 2000/2001 годы N = 87 лет), мм

Статистический параметр Весеннее половодье (IV - VI) Межень (VII - III) Год (IV - III)

РВ ИВ ЕВ ±VB РМ ИМ ЕМ ±VM РГ КГ ЕГ ±VГ

Среднее значение 329 87 141 -101 341 56 388 0 670 143 529 0

Минимальное значение 218 51 -22 -233 217 29 218 -351 481 85 312 -351

Максимальное значение 464 134 226 13 481 107 517 150 924 216 665 150

Амплитуда 246 83 248 246 264 78 299 501 443 131 353 501

Стандарт 47 18 43 50 58 16 61 87 84 28 71 87

Коэффициент вариации Слг 0,14 0,21 0,31 -0,49 0,17 0,28 0,16 - 0,13 0,20 0,13 -

Коэффициент асимметрии Сч 0,29 0,42 0,61 -0,98 0,34 0,57 0,31 - 0,26 0,40 0,26 -

г[1] -0,08 0,18 0,18 0,17 0,08 0,44 0,36 0,48 0,11 0,51 0,53 0,48

Кб)

№ 2' 2013

Анализ временных рядов ЭВБ позво- ^ лил разделить всю территорию бассейна реки Волги на четыре района по условиям ц

б

стационарности основных статистических й параметров элементов водного баланса, а также по их внутрирядной и межряд- М

ной связи. По изменчивости структуры ®

и динамике элементов водного баланса 8 §

территория бассейна Волги разделена на В

четыре водохозяйственных района: ®

Верхняя Волга (частные водосборы Ивань- ^ ^

ковского, Углического, Рыбинского и Ни- £ оо

жегородского водохранилищ), Средняя § I I

Волга (частные водосборы Чебоксарского ® ^

и Куйбышевского водохранилищ), река

о В й

Кама (частные водосборы Камского, Во- * !?

нн ^

ткинского и Нижнекамского водохрани- д

лищ) и Нижняя Волга (частные водосборы § ®

Саратовского и Волгоградского водохра- 2

нилищ). По климатическим условиям об- § § ширную территорию бассейна реки Волги ® 3

климатологи разделяют на три зоны: Верх- ^ I ей

нюю Волгу, Среднее Поволжье и Нижнее 3 ^ ^

Поволжье [3]. Для временных рядов вы- д ® й

® \ о

деленных водохозяйственных районов § ^ й

определены выборочные оценки основ- аЙ о

М эЫ

ных статистических параметров (табл. 2). м ^ ей

В период половодья среднемноголет- ^ м ^

ние значения осадков колеблются от <2 ® 3

261 мм (354 км3) - для водосбора Волгоград- Ц <5 И

ского водохранилища до 344 мм (468 км3) & £

н

- для водосбора Камского водохранили- И

я 2

Волги и Камы, достаточного увлажнения ®

ни

- район Средней Волги и недостаточного ^

увлажнения - район Нижней Волги.

Прежде чем характеризовать изменение речного стока за период половодья и межени, отметим, что для бассейна реки Волги характерно преимущественно

2 м

е

н 8 Ь

В 5 о

й

ща. Характерной особенностью распреде- В а §

ления осадков по территории является то, ® и х

что наибольшие значения приходятся на 5® ч

районы Верхней и Средней Волги и район ^ ^

Камы, а наименьшие значения - на район ® °

ей Сй

Нижней Волги - 261...270 мм (355...367 км3). I? <5

Аналогичная тенденция свойственна ^ <5

и для меженного периода. Наиболее ^

высокие значения осадков относятся к § ^

районам Верхней и Средней Волги и рай- § Д

ону Камы - 336.383 мм (457.521 км3), а ® и

наименьшие значения осадков относятся И Ц2

к району Нижней Волги - 222.249 мм ®

(302.339 км3). По степени валовой ув- м

лажненности территорию бассейна реки § н

но

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Волги можно разделить на три зоны: из- ^ « быточного увлажнения - районы Верхней

1-Н 1-Н 1-Н 1 > 1—1 Р-И о 0,18 -0,05

Р-Е со ю 0,42 0,48

и > 1 1 0,31

ч о и и н о 0,12 0,11

и Й со ю 0,22 0,50

со о со 0,15 0,22

Межень (VII - III) Р-И 1 со ю со 0,23 | о*

Р-Е 00! 0,86 0,38

±VM 1 1 0,31

ЕМ ^ ю со 0,14 0,03

ИМ 00! 0,31 0,42

§ Рч со ю ^ 0,21 0,12

Весеннее половодье (IV - VI) Р-И о 00 0,27 -0,03

Р-Е со ^ 0,27 0,06

±VB со о - -0,52 -0,03

м н ь- 00 0,17 0,26

м й со ю 0,23 0,21

РВ со со 0,17 | 0,02 |

Статистический параметр | Среднее значение | Коэффициент вариации Сл,

и г л о В

й е н

и Н

н о й

а р

й

ы

н

н

е

в

т

с

й

я

з

о

х

о

д

о В

Год (IV - III) Р-И 4 со 4 0,20 с с

Р-Е 00 со -2,06 а о; с

Г > 1 1 ^ о; с

и н 4 со 4 0,25 ^ с

и Й 00 со 0,46 с

-И 2 о 5 0,18 с- с с

Межень (VII - III) Р-И 3 3 0,31 а с

Р-Е о 5,59 со; с

±VM | 1 1 ^ о; с

ЕМ 3 3 0,24 с

ИМ о 0,57 ОС 1С с

РМ 4 со 0,30 о; о; с

Весеннее половодье (IV - VI) Р-И 7 2 0,28 Ь'О 0

Р-Е 3 ю 1 0,64 с- с

±VB -104 -0,76 с с

м н 3 1 0,60 се с

ИВ 00 ^ 0,50 сс се с

РВ 6 со 2 0,24 с с 1

Статистический параметр | Среднее значение | Коэффициент вариации С„

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

Анализ изменения среднемноголет-него стока за период половодья по частным водосборам бассейна реки Волги показал следующее: для района Верхней Волги среднемноголетнее значение стока колеблется от 110 до 149 мм (от 150 до 203 км3), для района Средней Волги - от 46 до 79 мм (от 62 до 107 км3), для района Камы - от 121 до 191 мм (от 164 до 260 км3) и для района Нижней Волги - от 33 до 62 мм (от 45 до 84 км3). Аналогичная картина характерна и для периода межени. Так, например, для района Камы среднемноголетний сток за период межени составляет от 75 до 130 мм (от 102 до 177 км3), а для района Нижней Волги - от 10 до 29 мм (от 14 до 39 км3).

В районах Верхней и Средней Волги, реки Камы сток за период половодья хорошо коррелирует с атмосферными осадками, выпадающими на водосборе за зимний и весенний периоды. Коэффициент корреляции изменяется от 0,45 до 0,60, чего нельзя сказать про район Нижней Волги (г = 0,18...0,27). Анализ связи между стоком и эффективными осадками (разность «осадки - испарение») показал, что существует более тесная связь в годовом и сезонном разрезе. Для весеннего половодья эта связь изменяется от 0,45 (для частного водосбора Воткинского водохранилища) до 0,99 (для частного водосбора Куйбышевского водохранилища), а для межени - от 0,61 до 0,99 (для тех же частных водосборов). Это свидетельствует о том, что для оценки величин речного стока можно выбрать прежде всего «эффективные осадки», а также собственно атмосферные осадки и суммарное испарение.

В отношении изменения испарения по территории обнаруживается несколько иная ситуация. Так, например, в период половодья, если для зоны избыточного увлажнения испарение изменяется от 65 до 123 мм (от 88 до 167 км3), то для зоны достаточного увлажнения - от 153 до 191 мм (от 208 до 260 км3), для зоны недостаточного увлажнения - от 109 до 132 мм (от 148 до 180 км3). Периоду межени свойственна достаточно высокая степень испарения для всех частных водосборов,

и она изменяется от 325 до 441 мм (соответственно от 442 и до 600 км3).

Изменению бассейновых влагозапа-сов свойственна следующая закономерность. В период половодья они имеет накопительное свойство, а в период межени происходит сработка, и в результате годовой баланс стремится к нулю. Для района Верхней Волги в период половодья накопление бассейновых влагозапасов изменяется от 88 до 105 мм (от 120 до 143 км3), для района Средней Волги - от 96 до 97 мм (от 130 до 132 км3), для района Камы -от 101 до 106 мм (от 137 до 144 км3), для района Нижней Волги - от 102 до 106 мм (от 139 до 144 км3). Таким образом, бассейновые влагозапасы, накопленные в период половодья, участвуют в формировании дополнительного испарения и поверхностного стока в период межени, благодаря чему повышается степень водо-обеспеченности территории частного водосбора речного бассейна в период межени.

Анализ тенденции изменения ЭВБ показывает, что для района Верхней Волги в период межени характерно увеличение речного стока за исследуемый период (1914/1915) на 14 мм за 87 лет. В то же время элементам водного баланса -атмосферным осадкам и испарению свойственно умеренное снижение. За рассматриваемый период осадки уменьшаются на 18 мм за 87 лет, а испарение - на 15 мм за тот же период. Бассейновые влагозапасы, которые участвуют в компенсации как речного стока, так и испарения, увеличиваются на 17 мм за 87 лет. К сожалению, исследование изменения ЭВБ за период половодья показало незначительное изменение этих элементов, особенно малозаметное изменение свойственно речному стоку (3 мм/87 лет) и испарению (5 мм/87 лет).

Наиболее заметные изменения ЭВБ происходят в районах Средней и Нижней Волги и районе Камы за период межени (лето - осень и зима). Как видно из рисунка, за период межени в целом по бассейну реки Волги атмосферные осадки имеют тенденцию к увеличению на 42 мм за 87 лет. Аналогичная тенденция наблюдается в динамике речного стока, увеличение которого составляет 25 мм за 87 лет. Суммарное испарение и изменение бассейновых влагозапасов за рассматриваемый период уменьшаются соответственно на 113 мм/87 лет и 130 мм/87 лет.

RМ, мм

10 20 30 40 50 60 70 80 90 б

г VМ, мм

10 20 30 40 50 60 70 80 90

10 20 30 40 50 60 70 80 90

в г

Изменение элементов водного баланса в бассейне реки Волги за период межени с 1914/1915 по 2000/2001 годы: а - атмосферные осадки; б - речой сток; в - суммарное испарение; г - изменение бассейновых влагозапасов.

Выявленное увеличение и снижение ЭВБ наиболее заметно за последнее 30-летие. При этом изменения ЭВБ за период весеннего половодья не столь существенны по сравнению с периодом межени. Одна из причин увеличения речного стока объясняется увеличением атмосферных осадков за рассматриваемый 87-летний период.

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

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

Средней Волги и Камы и колеблются от 0,15 до 0,19. Изменчивость суммарных осадков за период половодья и межени выше, чем изменчивость годовых осадков.

Коэффициент вариации речного стока за период половодья, начиная с частного водосбора Иваньковского водохранилища, постепенно уменьшается с 0,33 и на частном водосборе Куйбышевского водохранилища достигает 0,27, а затем, начиная с частного водосбора Саратовского водохранилища, повышается с 0,48 до 0,74 (частный водосбор Волгоградского водохранилища). Для района Камы коэффициент вариации речного стока в период половодья достаточно высок - 0,57 (частный водосбор Воткинского водохранилища). Для периода межени изменчивость речного стока также имеет высокое значение и колеблется от 0,33 (частные водосборы Чебоксарского и Камского водохранилищ) до 0,70 (частный водосбор Волгоградского водохранилища).

Наибольшего значения коэффициент вариации испарения за период половодья достигает в частных водосборах бассейна реки Камы и в районе Нижней Волги, а в районе Верхней Волги коэффициент вариации снижается и колеблется в пределах от 0,17 до 0,21.

Меженному периоду свойственно понижение коэффициента вариации до 0,06...0,10 для района Верхней Волги,

тогда как для района Камы и Нижней Волги коэффициент вариации повышается и достигает 0,32 (частный водосбор Волгоградского водохранилища).

Выводы

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

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

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

1. Исмайылов Г. Х., Федоров В. М.

Межгодовая изменчивость и взаимосвязь элементов водного баланса бассейна реки Волги // Водные ресурсы. - 2008. - Т. 35. - № 3. - С. 259-276.

2. Исмайылов Г. Х., Муращенкова Н. В. Пространственно-временные закономерности изменчивости и взаимосвязи ЭВБ половодья и межени в бассейне реки Волги в ХХ веке: Современные проблемы стохастической гидрологии и регулирования стока: Труды Всероссийской научной конференции, посвященной памяти выдающегося ученого-гидролога, профессора А. В. Рождественского. - М.: РосНИИВХ, 2012. - С. 354-363.

3. Гидрометеорологические условия Волжского региона и современные изменения климата / О. А. Анисимов [и др.] // Метеорология и гидрология. - № 5. -2011. - С. 33-42.

Материал поступил в редакцию 28.04.12. Исмайылов Габил Худуш оглы, доктор технических наук, профессор, зав. кафедрой «Гидрология,метеорология и регулирование стока» Тел. 8 (499) 976-23-68 E-mail: [email protected] Муращенкова Наталья Владимировна, кандидат технических наук, доцент Тел. 8 (499) 976-17-45 E-mail: [email protected]

УДК 502/504 : 551.311.2 : 551.33 О. Я. МАСЛИКОВА

Институт водных проблем РАН

ЗИМНИЕ НАВОДНЕНИЯ В УСЛОВИЯХ КРИОЛИТОЗОНЫ, ВЫЗВАННЫЕ ЛЕДОВЫМИ ЯВЛЕНИЯМИ НА РЕКАХ

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

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

There are given investigation results of the disturbance waves influence downstream of hydraulic works on bottom and banks deformation including under the conditions of catastrophic floods caused by the impact of release waves or dams breakage and ice jams formation. A special attention is given to the investigation of parameters which characterize permanently frozen soils.

Channel flow,downstream, waves of releases,ice,jams,permanently frozen soil,math-ematical models, laboratory simulation.

№ 2' 2013

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

и

i Надоели баннеры? Вы всегда можете отключить рекламу.