Динамика Подземных Вод Реферат

Динамика Подземных Вод Реферат



➡➡➡ ПОДРОБНЕЕ ЖМИТЕ ЗДЕСЬ!






























Динамика Подземных Вод Реферат
1. Движение воды в зонах аэрации и насыщения
2. Движение подземных вод в водоносных пластах. Определение скорости движения подземных вод
3. Установившееся и неустановившееся движение подземных вод. Методы моделирования фильтрации
4. Приток воды к водозаборным сооружениям
4.1 Приток безнапорных вод в совершенную горизонтальную дрену (канаву)
4.2 Расчет притока грунтовых вод в скважину
4.3 Расчет притока напорных вод в совершенную дрену
4.4 Расчет притока артезианских вод в скважину
5. Методы определения коэффициента фильтрации горных пород. Определение радиуса влияния
Гидрогеология — наука, изучающая подземные воды Земли, их историю, происхождение, формирование, состав, режим, геологическую и геохимическую деятельность. Но главной целью гидрогеологических исследований остается выявление новых резервов питьевой воды. В связи с быстрым ростом численности населения нашей планеты проблема запасов природных вод становится особенно острой. В ряде районов уже сейчас испытывается большой недостаток пресной воды, некоторые страны импортируют воду из других государств.
По количеству водных ресурсов Россия занимает одно из первых мест в мире. Однако промышленность и население нашей страны распределены неравномерно, в связи с чем и у нас в отдельных районах возникают очень большие потребности в питьевой воде. Такое положение сложилось в отдельных районах Урала. В РФ, к тому же, есть ряд засушливых районов с отсутствием водных ресурсов или острым их недостатком.
Роль подземных вод в жизни человека, животных и растений исключительно велика, в связи с этим значение гидрогеологии для народного хозяйства переоценить трудно.
1. Движение воды в зонах аэрации и насыщения

В зоне аэрации, т. е. в толще пород, расположенной между дневной поверхностью и зеркалом грунтовых вод, находятся:
а) водяной пар, заполняющий поры породы;
б) гигроскопическая влага, обусловливающая гигроскопическую влажность пород;
в) пленочная вода, обволакивающая зерна пород в виде пленок различной толщины, и
г) капиллярная вода, располагающаяся в виде капиллярной каймы над зеркалом грунтовых вод.
Движение подземных вод в зоне аэрации может происходить в виде передвижения пара, в виде пленочного движения, свободного просачивания и капиллярного движения.
Движение парообразной и гигроскопической влаги. А. Ф. Лебедевым было экспериментально доказано, что влага в парообразном состоянии передвигается от участка с большей упругостью водяного пара к участку с меньшей его упругостью. Упругость же зависит от температуры и влажности пород. Таким образом, если между различными участками горных пород появляется разница в температуре или влажности, возникает движение водяных паров. При одинаковой температуре движение направлено от более влажных частиц к менее влажным; при одинаковой влажности — от более к менее нагретым. Поэтому летом парообразная влага движется сверху вниз, а зимой — снизу вверх.
Гигроскопическая влага также передвигается в порах пород в виде водяного пара.
Движение воды в пленочном состоянии. По А. Ф. Лебедеву, движение воды в пленочном состоянии происходит под действием молекулярных сил и не подчиняется влиянию силы тяжести.
Рассмотрим движение пленочной воды на примере. Допустим, что мы имеем две одинаковые по диаметру частицы породы, соприкасающиеся между собой. Частица с центром О 1

покрыта пленкой воды толщиной Р 1
, а
вторая частица — более тонкой пленкой, толщиной Р 2
.
Рассмотрим влияние частиц породы на частицу воды, расположенную в точке С.
Легко убедиться, что расстояние О 1
С=R+P 1
и оно больше, чем О 2
С=R+P 2
т. е. частица 2 будет оказывать большее притяжение на частицу воды в точке С, чем частица породы с центром О 1
,
В результате частица воды С
перейдет на пленку, обволакивающую частицу породы 2. Движение частиц воды происходит до тех пор, пока толщина пленок на обеих частицах породы станет одинаковой.
Движение воды в виде просачивания. Просачивание в породах может происходить в виде отдельных струек и в виде сплошной массы воды. В первом случае отдельные струйки воды движутся самостоятельно, разрозненно. Вначале происходит смачивание частиц грунта, после чего под действием сил тяжести избыточная вода в виде гравитационной просачивается вниз.
Такой вид движения Г. Н. Каменский назвал свободным просачиванием.
Второй вид движения наблюдается в случае, если породы насыщены водой полностью. Движение воды здесь происходит сверху вниз под действием силы тяжести. Этот вид движения влаги назван инфильтрацией.

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

В породах, насыщенных водой, т. е. в зоне насыщения, движение воды может происходить в двух формах:
1) ламинарного,
при котором струйки воды текут параллельно, без перемешивания и
2) турбулентного,
при котором происходит хаотическое движение частиц жидкости и интенсивное перемешивание ее слоев. Переход от ламинарного движения к турбулентному и обратно происходит при достижении определенной скорости частиц жидкости, называемой критической скоростью.
Движение подземных вод в нескальных породах происходит по типу ламинарного.
Чтобы установить закономерности движения жидкости в породах, французский ученый X. Дарси в 1856 г. поставил несложный опыт, который заключался в следующем. В цилиндр, наполненный песком, наливали слой воды, поддерживая ее уровень постоянным. Вода после просачивания через песок выливалась через кран в нижней части цилиндра. В цилиндр были вставлены изогнутые трубки, так называемые пьезометры.
Вода в них устанавливалась на различных уровнях (в верхнем пьезометре — выше) в связи с тем, что в процессе фильтрации через поры грунта вода преодолевала сопротивление и на это терялась часть напора.
В результате проведенных исследований Дарси установил, что количество воды, профильтровавшейся через песок в единицу времени (расход, О),
прямо пропорционально разности уровней воды в пьезометрических трубках (∆Н=Н 2

— Н 1
),
площади поперечного сечения цилиндра (F)
и некоторому коэффициенту пропорциональности (К)
и обратно пропорционально высоте слоя песка (I).
Оказалось, что коэффициент К
зависит от свойств песка и его стали называть коэффициентом фильтрации (Кф).
Эта зависимость получила название закона Дарен и обычно записывается в следующем виде (1):
обозначают буквой / и называют напорным, градиентом
или гидравлическим уклоном.
Тогда можно записать
Если разделить обе части уравнения на F, то получим скорость фильтрации (υ) (2):
Таким образом, скорость фильтрации прямо пропорциональна коэффициенту фильтрации и напорному градиенту. Формула (3) представляет собой уравнение прямой линии, в связи с чем закон Дарси называют линейным законом фильтрации.

Если в выражении (3) принять I=1 ,
что имеет место при уклоне, равном 45°, получим
т. е. коэффициент фильтрации — это та скорость просачивания, которую имел бы поток при уклоне, равном единице.
Не следует при этом смешивать скорость фильтрации со скоростью движения частиц воды. Дело в том, что Дарси при расчетах принимал площадь поперечного сечения потока (F)равной сечению цилиндра, тогда как в действительности вода передвигалась в породе только по порам. Чтобы получить действительную скорость (и)движения подземных вод в порах грунта, необходимо расход воды разделить на площадь поперечного сечения и пористость грунта (n).
Это выражение показывает, что действительная скорость движения подземных вод больше скорости фильтрации, так как величина пористости всегда меньше единицы.
Необходимо заметить, что коэффициент фильтрации выражают в м/сут, хотя в некоторых случаях применяют см/с и км/год.
Если движение подземных вод происходит в крупных пустотах горных пород, то оно становится турбулентным и подчиняется нелинейному закону фильтрации, который выражается уравнением Шези
— Краснопольского

Таким образом, скорость фильтрации при турбулентном движении пропорциональна коэффициенту фильтрации и напорному градиенту в степени ½
2. Движение подземных вод в водоносных пластах. Определение скорости движения подземных вод

Для определения направления движения подземных вод используют карты гидроизогипс,
на которых в виде изолиний показан «рельеф» зеркала грунтовых вод. Перпендикуляры к гидроизогипсам, направленные в сторону снижения отметок, называются линиями тока,
показывающими направление движения грунтовых вод.
По взаимному расположению гидроизогипс и линий тока потоки грунтовых вод разделяют на плоские и радиальные (рис. 3).В плоском потоке гидроизогипсы в плане имеют вид параллельных прямых и линии тока при пересечении с ними образуют сеть прямоугольников. Плоский поток может иметь место в междуречьях; между рекой и дреной, текущими параллельно; в случае дренирования грунтовых вод горизонтальными выработками (канавами, штольнями).
В радиальном потоке гидроизогипсы представляют соб»й систему кривых линий, а линии тока имеют вид радиусов. Наиболее наглядным примером радиального потока может быть приток воды в колодец или скважину во время интенсивного водоотбора. Радиальный поток может быть расходящимся
(например, возле излучины реки) и сходящимся
(к водозабору). При расходящемся потоке ширина его по
направлению движения увеличивается, а при сходящемся, наоборот, уменьшается.
График изменения содержания ионов хлора в подземных водах при определении действительной скорости потока
Скорость движения подземных вод можно определить несколькими способами. Один из них основан на введении в воду поваренной
соли. На некотором расстоянии от опытной скважины (шурфа или колодца) проходят наблюдательную скважину, которую закладывают ниже по направлению движения подземных вод. Перед началом опыта определяют содержание хлора в опытной и наблюдательной выработках. Затем в опытную выработку вводят раствор поваренной соли, в котором концентрация ионов хлора в 2000 раз выше, чем в подземных водах. Естественно, время ввода соли (t 1
) необходимо отметить. Через каждые 10 мин из наблюдательной скважины отбирают пробы воды и при помощи азотнокислого серебра определяют содержание хлора. Данные анализов наносят на график (рис, 3) и находят время прохождения пика (t 2
). Действительная скорость
Где l - расстояние между выработками, м.
Этот способ очень удобен, но применение его невозможно при естественном содержании хлора в воде свыше 500—600 мг/л и при резких неровностях водоупорного слоя. В первом случае анализами трудно определить изменения содержания хлора, во втором — более тяжелый, чем вода, раствор поваренной соли может задержаться в понижениях водоупора.
Можно также применять органические красители,
присутствие которых в воде обнаруживается при ничтожно малых концентрациях (до 10- 6
%). Для этого применяют флуоресцеин, имеющий при слабых концентрациях зеленовато-желтый цвет, метиленовый синий краситель и др. Для определения содержания красителя в воде используют флюороскоп — набор стеклянных трубок с разной концентрацией красителя. Сравнивая цвет воды в отобранных пробах с цветом трубок-эталонов, легко и быстро можно определить содержание красителя в пробе воды. Затем строят график изменения во времени содержания красителя в воде и аналогично вышеописанному способу определяют скорость движения подземных вод.
Скорость движения подземных вод можно определять и электролитическим
способом. Для этого в опытную скважину вводят электролит (обычно хлористый аммоний) и следят за изменением электропроводимости между опытной и наблюдательной скважинами. Для этой цели используют миллиамперметр, по данным которого строят график изменения силы тока во времени.
Новейшие достижения физики и химии позволяют использовать «меченные атомы»
— изотопные индикаторы. Высокая чувствительность и простота радиоактивных измерений позволяют фиксировать минимальное количество изотопов в подземных водах.
3. Установившееся и неустановившееся движение подземных вод. Методы моделирования фильтрации

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

Для моделирования фильтрации в основном используются гидравлическая и электрическая аналогии, реализуемые на сплошных и сеточных моделях.
Сплошные гидравлические
модели, представленные фильтрационными лотками различных видов, в гидрогеологических расчетах применяются редко.
В развитии методов моделирования фильтрации подземных вод основная роль принадлежит сплошным и сеточным электрическим моделям, основанным на использовании метода электрогидродинамических аналогий
(ЭГДА), сущность которого наглядно представляется сопоставлением основных законов движения фильтрационного потока и электрического тока:
где Q — расход; F — площадь поперечного сечения потока; Н — напор; х-—расстояние;I— сила тока; с — удельная проводимость, ; р — удельное сопротивление; площадь поперечного сечения проводника; U— электрический потенциал, l — длина проводника.
Приведенная формула закона Ома получена путем несложных преобразований
Идентичность записи законов Дарси и Ома очевидна. В них соответствуют физические характеристики — коэффициент фильтрации Кф
и удельная проводимость с
(физическое подобие), силовые характеристики — напор Н
и потенциал U
(динамическое подобие) и, наконец, расход потока Q и сила тока (кинематическое подобие).
На сплошных
моделях ЭГДА фильтрационный поток моделируется сплошным электрическим полем, геометрически подобным. Для этого применяются электропроводная бумага и электролиты. Электропроводная бумага изготавливается с удельным сопротивлением от 100 до 100 000 Ом/см, в зависимости от количества содержащихся в ней сажи и графита.
Участки поля с различной проницаемостью пород моделируются кусками бумаги различной удельной проводимости. Между собой участки модели скрепляются специальным электропроводным клеем.
Электролиты также широко используются в качестве материала модели и обычно представляют собой растворы солей, причем наибольшее распространение получили водные растворы поваренной соли и медного купороса. Кроме того, можно использовать электропроводные краски, клеи, электропроводный картон, гипс и т. д.
Определение приведенного потенциала на моделях ЭГДА производится с помощью мостовой измерительной схемы.
При составлении сеточных
моделей поток разбивается на отдельные блоки, центры которых связываются электрическими резисторами. В таких моделях геометрическое подобие модели и объекта не сохраняется.
4. Приток воды к водозаборным сооружениям

Среди водозаборных сооружений мы будем рассматривать такие горные выработки, как дрены (канавы) и скважины. В гидрогеологии горные выработки разделяют на совершенные и несовершенные.
Гидродинамически совершенной называется горная выработка, вскрывающая водоносный горизонт от кровли до подошвы.
Рис. 1.Схемы совершенной (а)
и несовершенной (б, в)
выработок
Приток воды к
ней происходит по всей поверхности соприкосновения стенок выработки с водоносным горизонтом (рис. 1, а).
Если же выработка не доходит до водоупора, она называется несовершенной по степени вскрытия
водоносного горизонта (рис. 1, б). Зачастую выработки закрепляются от обрушения, цементируются скважины оборудуются обсадными трубами, фильтрами и т. п. Естественно, что приток воды в такие выработки затруднен и их называют несовершенными по характеру вскрытия
водоносного горизонта.Основные уравнения притока воды к водозаборам (скважинам и дренам) будем выводить при условии совершенства выработок.
Представим себе плоский поток грунтовых вод. Гидравлический градиент I в данном случае равен
где х —
расстояние между сечениями h 1
и h 2


Если мы будем сближать сечения h 1
и h 2

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

После устройства дрены скорость движения воды в ней увеличивается и уровень воды понижается на величину S, которую в гидрогеологии принято называть величиной понижения.
Иными словами, величина понижения представляет собой разницу между статическим и динамическим
уровнями. Мощность водоносного горизонта до понижения обозначим через H, глубину воды в дрене — через h o
. В результате понижения уровня в дрене в водоносном горизонте образуется депрессионная воронка,
показанная на рис. 50 сплошной жирной линией. Расстояние R,
на которое сказывается влияние понижения, называют радиусом влияния.

Для расчета притока воды в дрену Q выбираем на расстоянии х
от стенки дрены сечение с напором к,
которое находится в интервале от нуля до R.

В общем виде приток воды в дрену будет равен выражению (13). Подставим сюда величину площади фильтрации
При расчете притока воды в дрену удобно пользоваться понятием единичного притока д,
т. е. притока воды на единицу длины дрены
Отсюда элементарная формула для расчета притока воды
Разделим переменные в выражении (18), т. е. умножим обе его части на dх
и проинтегрируем
Формула (22) выражает величину единичного притока с одной стороны дрены. Для получения полного притока воды в дрену необходимо умножить единичный приток на два, а затем — на длину дрены. Приток воды в торцы дрены обычно не учитывают, так как он при большой длине дрены составляет ничтожную долю.
По формуле (22) можно рассчитать расход плоского грунтового потока. Подставив вместо радиуса влияния расстояние между сечениями, равное I
, получим
Исследуя выражение (22), мы сможем решить одну из весьма важных задач в гидрогеологических расчетах — вывести уравнение депрессионной кривой. Построение депрессионной кривой необходимо при возникновении угрозы затопления подземными водами котлованов, подвалов зданий и т. п..
Изменив пределы интегрирования в выражении (20) по X
от 0 до х,
а по У
от h 0
до h К
получим
Естественно, что приток воды в выражениях (22) и (27) одинаков, т. е.
Для построения депрессионной кривой мы задаемся величиной hо в зависимости от 5, мощность водоносного горизонта H легко получить по данным бурения, величину радиуса влияния можно найти по эмпирическим формулам (об этом мы поговорим позже).
На миллиметровой бумаге строим разрез через дрену и котлован (рис. 52) и, задаваясь разными значениями х(х и
x 2
,.
.., х п
),
например 10, 20, 30 и т. д. метров, получаем величины h(h,
h 2
,..., h п
).
Соединив полученные точки плавной линией, получим кривую депрессии. Если она проходит через котлован, строят новую кривую, задавшись большей величиной понижения и, естественно, меньшим значением глубины воды в дрене. Построение производят до тех пор, пока депрессионная кривая не опустится ниже дна котлована.
4.2 Расчет притока грунтовых вод в скважину

Здесь мы имеем дело не с плоским потоком, как в предыдущем примере, а с радиальным. На рис. 10 показаны все обозначения, которые нам ясны из предыдущей задачи, кроме г
— радиуса скважины.
Расчет начинаем с уравнения (13). Площадь притока воды равна площади боковой поверхности цилиндра, радиус которого равен х,
т. е.
Разделяем переменные (т. е. умножаем обе части уравнения на ) и интегрируем
Интегрирование по х
производим не от нуля, а от r
— стенки скважины, в результате получаем формулу Дюпюи
Рис. 2. Схема притока грунтовых вод в совершенную скважину
В таком виде использовать выражение (2) не совсем удобно, так как в нем присутствует натуральный логарифм. Подставим вместо него десятичный (1nх = 2,30 lgх), а вместо л его значение и получим более удобное выражение для расчета притока безнапорных вод в скважину:
И, подставив его в зависимость 54), получим
Для построения депрессионной кривой возвращаемся к (35) и изменяем пределы интегрирования: по X
от r
до х
а по Y от h 0
до h:

Решаем равенство относительно h
и получаем уравнение кривой депрессии
4.3 Расчет притока напорных вод в совершенную дрену

Площадь фильтрации в сечении h, расположенном на расстоянии х
от стенки дрены, будет равна
Мы здесь снова не учитываем приток воды через торцы дрены. Подставляем площадь в (43):
Разделяя переменные и интегрируя, получим
Выражение (47) представляет собой единичный приток артезианских вод в один из бортов канавы. Полный приток составит
Если нам необходимо получить уравнение депрессионной кривой, то (43) нужно проинтегрировать по Х от 0 до ж, а по Y
от h о
до h:

Решаем уравнение (51) относительно h:

Анализируя выражение (52), мы видим, что это уравнение прямой линии. На самом деле депрессионная кривая криволинейна.
4.4 Расчет притока артезианских вод в скважину

В выражение (14) подставляем величину площади фильтрации, которая равна
Разделяя переменные и интегрируя по X
от r
до Н,
а по Y
от h о

до Н,
получим
Переведем натуральный логарифм в десятичный и подставим значение я. Получим выражение для расчета притока артезианских вод в совершенную скважину:
Для расчета кривой депрессии возвратимся к (55) и сменим пределы интегрирования: по X от r до х,
а по Y
от hо
до h:

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

Дебит (по-фр.— сбыт, расход) — количество воды, нефти, газа, даваемое источником, колодцем, скважиной за единицу времени. Единицами измерения дебита для подземных вод являются м 3
/с или м 3
/сут, для нефти — т/сут. Удельный дебит получают при делении величины дебита на величину понижения уровня
5. Методы определения коэффициента фильтрации горных пород. Определение радиуса влияния

Как видно из приведенных в предыдущем параграфе формул, в большинстве из них присутствуют величины коэффициента фильтрации и радиуса влияния, на методах определения которых мы и остановимся.
Определение коэффициента фильтрации горных пород можно проводить на основании использования эмпирических формул, лабораторных данных и полевых опытов.
Эмпирические формулы позволяют быстро определить коэффициент фильтрации горных пород на основании данных об их пористости и механического состава. Однако эти формулы дают лишь приблизительные представления о водопроницаемости пород и могут быть использованы только при предварительных расчетах.
Для определения коэффициента фильтрации песков с эффективным диаметром частиц, равным от 0,1 до 3,0 мм, и при коэффициенте неоднородности менее 5,0 можно применять формулу Хазена
где С—эмпирический коэффициент, равный, по О. К. Ла-нге,
п
— пористость породы, %; d э
~
эффективный диаметр частиц, мм; t°
— температура воды, °С.
Для определения эффективного диаметра и коэффициента-неоднородности необходимо построить интегральную (суммарную) кривую механического состава. Обычно результаты анализов механического состава лаборатория выдает в виде стандартной таблицы (табл. 1). Для построения суммарной кривой необходимо знать суммарные содержания фракций диаметром менее 0,005; 0,01; 0,05; 0,1; 0,25 и т. д. мм. Для этого в таблице справа налево мы суммируем содержания фракций и записываем ниже (под содержанием фракций, %).
Банк рефератов содержит более 364 тысяч рефератов , курсовых и дипломных работ, шпаргалок и докладов по различным дисциплинам: истории, психологии, экономике, менеджменту, философии, праву, экологии. А также изложения, сочинения по литературе, отчеты по практике, топики по английскому.










Название: Расчет динамики подземных вод
Раздел: Рефераты по геологии
Тип: курсовая работа
Добавлен 16:47:40 21 октября 2009 Похожие работы
Просмотров: 2842
Комментариев: 15
Оценило: 3 человек
Средний балл: 4.7
Оценка: неизвестно     Скачать


Динамика подземных вод
Курсовая работа: Расчет динамики подземных вод - BestReferat.ru
Динамика подземных вод | реферат [24,7 K], добавлен 09.05.2007
Реферат на тему: Геология. Динамика подземных вод
Динамика подземных вод - Геология - KazEdu.kz
Реферат : Динамика подземных вод - Xreferat.com - Банк...
Основы динамики подземных вод
ОСНОВЫ ДИНАМИКИ ПОДЗЕМНЫХ ВОД — Студопедия
Динамика подземных вод - реферат - скачать бесплатно
Динамика подземных вод .
Реферат : Динамика подземных вод - All-Referats.com - Сайт...
23 - Динамика подземных вод » СтудИзба
Роль Дружбы В Жизни Человека Декабрьское Сочинение
Личность Ивана Третьего Эссе
Сочинение На Тему Мой Домашний Питомец
Контрольная Работа За 1 Четверть 3 Класс
Эссе На Тему Творчество

Report Page