Механіка підземних споруд

[ виправити ] текст може містити помилки, будь ласка перевіряйте перш ніж використовувати.

скачати


Введення

Гірська кріплення - штучна споруда, що споруджується в виробках для запобігання обвалення навколишніх порід і збереження необхідних площ перерізів виробок у нормальному стані впродовж всього терміну експлуатації.

У загальному комплексі робіт, пов'язаних з видобутком корисних копалин, важливе місце займає кріплення та підтримання капітальних виробок. Для подальшого вдосконалення гірничо-підготовчих робіт у виробках особливу роль займає проблема розвитку методів прогнозування проявів гірського тиску і розрахунку кріплень. Із застосуванням чисельних методів і ПЕОМ методи вибору способів охорони виробок і проектування кріплень отримали достатньо достовірну наукову основу, що призводить до поліпшення стану виробок і сприяє підвищенню техніко-економічних показників підземної розробки.

Мета даного курсового проектування - придбання навичок проектування кріплення гірничих виробок з використанням сучасних ПЕОМ.

Вихідними даними для проектування є гірничо-геологічні і гірничотехнічні умови проведення виробок в реальних умовах з наступними характеристиками:

Перетин вироблення S св - 12,8 м2

Висота вироблення-3400 мм

Ширина вироблення - 5300 мм

Глибина залягання - 950м

Міцність порід, МПа / потужність, м: 1 аргіліти - 32 / 2,0;

алевроліти - 45 / 2,8; пісковики-58 / 6,9; вугілля - 18 / 2,6.

Обводненість: невелика

Відстань від підошви пласта - 4 / 0

Кут падіння порід - 9 градусів.

Напрямок проходки - під кутом

Термін служби - 4 роки, одиночна

Тип вироблення - розкривають

Умови проходки - БПР

Розрахунок внутрішніх зусиль у кріплення проводиться за програмою, розробленою авторами під керівництвом проф. Колоколова С.Б.

Курсовий проект складається з двох частин: розрахунку металевої арочної податливою кріплення та розрахунку монолітної бетонної кріпи.

1. Розрахунок гірського тиску на кріплення виробки

1.1. Визначення розрахункових опорів оточуючих вироблення порід

По заданому перетину вироблення визначаються основні геометричні розміри виробітку (висота, ширина, радіус). Дані за геометричними розмірами кріплення наведені в довідковій літературі (довідник по кріпленню). Дані в завданні наведені для перетину виробки у світлі, а для розрахунків необхідне використання геометричних розмірів для перетину виробки у проходці. Для цього розміри по ширині виробки повинні бути збільшені на 0,4 м, а по висоті на 0,3 м.

Нпр = h +300 мм (1.1) Нпр = 3400 +300 = 3700мм

Впр = b +400 мм. (1.2)

Впр = 5300 +400 = 5700мм

Розрахункова глибина Нр розміщення підземної гірничої виробки визначається за формулою:

Нр = Н * k, м (1.3)

де Н-проектна глибина закладення вироблення, м; k - коефіцієнт враховує відміну чинного напруженого стану масиву від власної ваги товщі порід до поверхні. Приймається в зонах тектонічних порушень рівним 1.5. і рівним 1.0 у всіх інших випадках.

Нр = 900 * 1 = 900м

Розрахунковий опір порід стисненню Rc визначається за формулою:

Rc = R * kc * kH2O, М па

R - середнє значення опору порід у зразку на одноосьовий стиск, МПа;

kc - коефіцієнт враховує додаткову нарушенность масиву порід.

kH 2 O - коефіцієнт обводнення

Rc аргіл = 1 * 32 * 0,75 = 24 МПа

Rc алеврит = 1 * 45 * 0,8 = 36 МПа

Rc песча = 1 * 58 * 0,9 = 52 МПа

Rc вугілля = 1 * 18 * 1 = 18 МПа

У разі обводнення проектованої вироблення зниження міцності порід за рахунок впливу вологи враховують шляхом зменшення міцності пісковиків до 20%, алевролітів до 40% і аргілітів до 50%.

Для оцінки раніше виконаної проектної документації і при використанні наявних нормативних документів допускається визначати опір порід стисненню у зразку R за коефіцієнтом міцності порід f відповідно до шкали М.М. Протодьяконова.

R = 10 * f (1.5)

де f - коефіцієнт міцності порід, що приймається за шкалою М.М. Протодьяконова.

Розрахунковий опір порід стисненню Rc по контуру поперечного перерізу виробки повинен визначатися з урахуванням всіх пересічних виробленням шарів (пластів) потужністю більше 0,1 м, що залягають на відстанях від контуру перерізу виробки: у покрівлі 1,5 B, в грунті і боках виробки - 1 B, де B - ширина виробки (див. рис.1). При наявності п шарів, що залягають по контуру поперечного перерізу виробки, з мінливістю Rc в межах 30% приймається для всієї вироблення усереднене значення розрахункового опору порід стиску, що визначається за формулою:

Rcy = (m1R1 + m2R2 ... + mnR) / (m1 + m2 + ... mn)

де R - розрахунковий опір порід стиску, МПа;

n - число породних шарів;

m - потужність i-го шару.

R покрівля = МПа

R грунт = МПа

R-пліч = МПа

Визначення необхідних характеристик кріплення

За методикою Слухаючи, для прогнозу стійкості горизонтальних виробок вугільних шахт як критерію стійкості приймається величина очікуваних зміщень на контурі перетину незакріпленої виробітку за весь термін її служби, що визначаються за формулою:

u = К a До q До S КВК tU т, мм

де, u т - зміщення (типові), що визначаються за графіками (рис.1.2);

Малюнок 1.2. - Графіки для визначення типових зміщень порід.

U т кр = 280; U т пліч = 1050; U т поч = 500

До q - коефіцієнт, що враховує напрямок зміщень порід і приймає значення 1 при визначенні зсувів з боку покрівлі або грунту (у вертикальному напрямку). При визначенні бічних зсувів порід (у горизонтальному напрямку) До q приймається за табл.1.1.

Таблиця 1.1

Напрямок

проходки

вироблення

Коефіцієнти До a і К q в залежності від кута падіння порід a, град


до 20

21-30

31 - 40

41-50

більше 50


До a

До q

До a

До q

До a

До q

До a

До q

До a

До q

За простиранню

1

0,35

0,95

0,55

0,8

0,8

0,65

1,2

0,6

1,5

Навхрест простягання

0,7

0,55

0,6

0,8

0,45

0,95

0,25

0,95

0,2

0,8

Під кутом до простяганню

0,85

0,45

0,8

0,65

0,65

0,9

0,45

1,05

0,35

1,1

До α кр = 0,85; До α-пліч = 0,85; До α поч = 0,85

До q - коефіцієнт, що враховує напрямок зміщень порід і приймає значення 1 при визначенні зсувів з боку покрівлі або грунту (у вертикальному напрямку). При визначенні бічних зсувів порід (у горизонтальному напрямку) До q приймається за табл.1.1.

До θ кр = 0,45; До θ поч = 0,45; До θ пліч = 0,45

До S - коефіцієнт впливу прольоту виробітку для покрівлі і грунту приймається рівним К S = 0,2 (b -1), для боків виробки До S = 0,2 (h -1) (b - проліт вироблення, h - висота виробки) ;

До s кр. пч. = 0,2 (5,7-1) = 0,94 (1.8)

До s пліч = 0,2 (3,7-1) = 0,54

КВ - коефіцієнт впливу інших виробок, рівний для одиночних виробок 1;

До t - коефіцієнт впливу часу, що дорівнює 1 для виробок, термін служби яких більше 15 лет.д.ля виробок термін служби яких більше року До t визначається за рис.1.3. Зсуви визначаються окремо для покрівлі, грунту і боків виробки.

Малюнок 1.3. - Графіки для визначення коефіцієнта Кt при терміні служби вироблення більше року.

U кр = 0,85 * 0,45 * 0,94 * 1 * 1 * 280 = 101 мм.

U поч = 0,85 * 0,45 * 0,94 * 1 * 1 * 500 = 180 мм.

U-пліч = 0,85 * 0,45 * 0,54 * 1 * 0,5 * 1050 = 108 мм.

1.3. Визначення розрахункових навантажень на кріплення

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

Нормативні навантаження на замкнену кріплення також визначаються за графіками рис.1.4. окремо для покрівлі, грунту і боків. Усереднена вертикальне навантаження визначається за значеннями нормативних навантажень з боку покрівлі і грунту, усереднена горизонтальне навантаження - за значеннями нормативних навантажень з боку боків.



Малюнок 1.4. - Графіки визначення нормативного навантаження на кріплення

Розрахункова вертикальне навантаження на 1 м довжини виробітку визначається за формулою:

Рв = k п * k н * k пр * Рнв, кПа, (1.9)

де k п - коефіцієнт перевантаження, що враховує мінливість навантаження (табл.1.2);

k н - коефіцієнт, який приймається для головних розкриваючих виробок дорівнює 1,1, для решти-1;

k пр - коефіцієнт умов проведення виробок, що дорівнює 1 при проведенні виробок буровибуховим способом.

Рнв - нормативна вертикальне навантаження.

Таблиця 1.2. - Значення коефіцієнта k п


U, мм

Значення k п для виробок


Розкриваючих

Підготовляють

до 50

1,25

1,1

51-200

1,1

1,05

201-500

1,05

1

більше 500

1

1

Рв = 1,1 * 1 * 1 * 45 = 50 кПа

Розрахункове навантаження на 1 м виробки з боку боків визначається за формулою:

Рб = kп * kн * Kпр * рнб, кПа, (1.9)

де РНБ - нормативна горизонтальне навантаження, кПа.

Рб = 1,1 * 1 * 1 * 29 = 32 кПа

2. Розрахунок металевої арочної податливою кріплення

За отриманими даними розрахункових опорів оточуючих вироблення порід, визначення зміщень контуру виробки і визначення розрахункових навантажень на кріплення. Вибір даних параметрів виробляється згідно розробленої методики Слухаючи.

Отримані дані повинні бути введені наступні дані:

1) Радіус криволінійної частини кріплення R і висота прямолінійної частини (стійок) h, (рис.2.1).

R = В / 2 (2.1)

R = 5700 / 2 = 2850 мм

h = HR (2.2)

h = 3700-2850 = 850 мм

Малюнок 2.1. - Розрахункова схема для металевої аркового кріплення.

2) Число частин розподілу кожної з вертикальних стійок 4 і число поділу криволінійної частини 9.

3) Вертикальна і горизонтальна навантаження на кріплення рівні відповідно 50 кПа і 32 кПа.

4) Кут нахилу стійкий до вертикалі дорівнює нулю.

Отримані результати розрахунку містять значення згинальних моментів і поздовжніх сил в кожній заданій точці кріплення (таблиця 2.1). На підставі даних значень знаходиться максимальний згинальний момент (за модулем) і відповідне цьому моменту нормальна сила.

Таблиця 2.1. - Розрахунок металевого кріплення на ЕОМ

i = х [i] = y [i] = М [i] = N [i] =

1-2.850.00 0.0000-142.5000

2-2.850.21 1.0002-142.5000

3-2.850.43 0.5554-142.5000

4-2.850.64 - 1.3344-142.5000

5-2.850.85 - 4.6692-143.0609

6-2.811.34 - 11.9141-143.7196

7-2.681.82 - 14.7291-142.4850

8-2.472.28 - 13.6014-138.6918

9-2.182.68 - 9.4427-133.0686

10-1.833.03 - 3.4556-126.5388

11-1.423.32 3.0324-120.1014

12-0.973.53 8.7480-114.7045

13-0.493.66 12.6401-111.1255

140.003.70 14.0177-109.8735

150.493.66 12.6401-111.1255

160.973.53 8.7480-114.7045

171.423.32 3.0324-120.1014

181.833.03 - 3.4556-126.5388

192.182.68 - 9.4427-133.0686

202.472.28 - 13.6014-138.6918

212.681.82 - 14.7291-142.4850

222.811.34 - 11.9141-143.7196

232.850.85 - 4.6692-143.0609

242.850.64 - 1.3344-142.5000

252.850.43 0.5554-142.5000

262.850.21 1.0002-142.5000

272.850.00 0.0000-142.5000

У нашому випадку після розрахунків М max = 14,01 кН * м, і N = 109,87 кН * м

Перевірка міцності здійснюється за формулою:

(2.3)

де W пр - момент опору з (табл.2.2);

F пр - площа поперечного перерізу з (табл.2.2)

Таблиця 2.2. - Характеристики спецпрофілю.

Найменування профілю

Вага 1 п. р. м, кг

Площа поперечно. перерізу, F см2

Моменти опору, см3

Моменти інерції, СМ4




Wx

Wy

Jx

Jy

СВП-17

17,06

21,73

50,3

57,9

243,4

382,3

СВП-22

21,9

27,91

81,3

77,8

428,6

566,3

СВП-27

27,0

34,3

100,2

101,5

639,5

763,1

СВП-33

33,8

42,46

138,5

148,0

1000

1228

R - розрахунковий опір вигину прокатної сталі може бути прийнято 230 МПа (2300кг/см2)

L - крок кріплення приймається 0,33; 0,5; 0,75; 1.

Таблиця 2.3.

Тип профілю


0,33


0,5


0,75


1


СВП-22


15 * 107

15 * 107

15 * 107

15 * 107


69,7 * 107

46 * 107

30,6 * 107

23 * 107


СВП-27


14,8 * 107

14,8 * 107

14,8 * 107

14,8 * 107


69,7 * 107

46 * 107

30,6 * 107

23 * 107


СВП-33

11 * 107

11 * 107

11 * 107

11 * 107


69,7 * 107

46 * 107

30,6 * 107

23 * 107

На підставі отриманого значення вибирається тип спецпрофілю.

СВП-22 з моментом опору W пр = 93,36 і площею поперечного перерізу F пр = 27,91 см2 з кроком кріплення 1 м. Ми обрали цей тип спецпрофілю і крок кріплення виходячи з рівняння.

15 * 107 ≤ 23 * 107

Арочні податливі триланковий металеві кріплення є найбільш поширеним видом кріплення горизонтальних і похилих гірничих виробок. Арка складається з трьох ланок: двох стійок і верхняка.

3. Розрахунок бетонного кріплення.

Розрахунок параметрів бетонної крепі, проводиться за тими ж самим параметрами що і для розрахунку арочної металевого кріплення, проте сама бетонне кріплення має ряд конструктивних відмінностей від металевого кріплення розрахованої вище. Схема бетонної крепі показана на малюнку 3.1

Отримані вище дані про вертикальному і горизонтальному тиску повинні бути використані при розрахунку міцності бетонної кріпи з використанням ПЕОМ. Для проведення розрахунку в програму повинні бути введені наступні дані:

Висота склепіння у просвіті може прийматися за СНиП:

h 0 = В / 3

де В - ширина виробки у проходці.

h 0 = 5700 / 3 = 1900мм

Для визначення товщини зводу в замку можна скористатися формулою С.С. Давидова:

де L 0 - полупролет вироблення, м;

f - коефіцієнт міцності порід по Протодьяконова.

м

Товщина зводу в п'яті:

d п = (1,25-1,50) d 0

d п = 1,35 * 0,25 = 0,325 м

Товщина стін:

з = (1,00-2,00) d п

з = 1,5 * 0,325 = 0,49 м

Для проведення розрахунку в програму повинні бути введені наступні дані:

Кількість елементів кріплення (зазвичай прямолінійний елемент, малий радіус, великий радіус). Для спрощення можна прийняти прямолінійний елемент і великий радіус.

Вертикальну і горизонтальну навантаження на кріплення.

Координату Х крайньої нижньої (лівої) точки кожного елемента

Координату У крайньої нижньої (лівої) точки кожного елемента

Х1 = В / 2; У1 = 1,8 + В / 3;

Х1 = 5,7 / 2 = 2,85 м; У1 = 1,8 +5,7 / 3 = 3,7 м

Х2 = В / 2; У2 = В / 3;

Х2 = 5,7 / 2 = 2,85 м; У2 = 5,7 / 3 = 1,9 м

Х3 = 0,383 * В; У3 = 0,115 * В;

Х3 = 0,383 * 5,7 = 2,18 м; У3 = 0,115 * 5,7 = 0,65 м

Радіус кривизни в метрах (для прямолінійного ділянки = 100)

Для малого радіуса r = 0,262 * В; R = 0,692 * У

r = 0,262 * 5,7 = 1,49 м; R = 0,692 * 5,7 = 3,9 м.

Кількість ділянок розбиття елементів:

- Для прямолінійної частини (парне число) 4-8

- Для криволінійної частини (парне число) 4-10, приймається з умови зручності розбиття криволінійної частини на кути по яких відкладаються отримані навантаження.

Відносну жорсткість елемента (у разі однакової товщини кріплення дорівнює 1)

Тому що при розрахунку приймаються симетричні бічні та вертикальне навантаження щодо осі Y, то для визначення внутрішніх зусиль у кріплення досить ввести координати X і Y тільки однієї половини контуру виробки.

Малюнок 3.1. - Розрахункова схема для бетонного кріплення.

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

Виділена із зведення арка одиничної ширини є бесшарнірной.

Отримані результати розрахунку містять значення згинальних моментів і поздовжніх сил в кожній заданій точці кріплення. На підставі даних значень проводиться розрахунок бетонної крепі в наступній послідовності.

Таблиця 3.1. - Розрахунок бетонного кріплення на ЕОМ

Nел Х YR NJ EI

1 2,85 3,70 100 4 1

2 2,85 1,90 1.49 4 1

3 2,18 0,66 3.94 4 1

i = x [i] = y [i] = M [i] = N1 [i] = N2 [i] =

1

2.85

3.70

19.80

-142.50

-142.50

2

2.85

3.25

19.92

-142.50

-142.50

3

2.85

2.80

13.56

-142.50

-142.50

4

2.85

2.35

0.72

-142.50

-142.50

5

2.85

1.90

-18.60

-147.57

-146.77

6

2.81

1.54

-32.61

-153.10

-150.88

7

2.67

1. 19

-37.50

-151.11.

-148.00

8

2.46

0.89

-33.55

-142.61

-139.35

9

2.38

0.66

-21.82

-131.88

-127.44

10

1.68

0.37

0.09

-122.46

-138.99

11

1.14

0.17

17.01

-115.45

-113.25

12

0.58

0.04

27.72

-111.40

-130.65

13

0.00

0.00

31.44

-130.65

-111.40

14

-0.58

0.04

27.72

-113.25

-115.45

15

-3.14

0.17

17.01

-138.99

-122.46

Продовження.

16

-1.68

0.37

0.09

-127.44

-131.88

17

-2.18

0.66

-21.82

-139.35

-142.61

18

-2.46

0.89

-33.55

-148.00

-151.11

19

-2.67

1. 19

-37.50

-150.88

-153.10

20

-2.81

1.54

-32.61

-146.77

-147.57

21

-2.85

1.90

-18.60

-142.50

-142.50

22

-2.85

2.35

0.72

-142.50

-142.50

23

-2.85

2.80

13.56

-142.50

-142.50

24

-2.85

3.25

19.92

-142.50

-142.50

25

-2.81

3.70

19.80

0.00

0.00

Після розрахунку внутрішніх зусиль на ПЕОМ ми визначаємо точку з координатами X і Y, де є максимальний згинальний момент М max = 37,50 і відповідна поздовжня сила N = 151,11 в кріплення. Стискуюче напруга від сили N визначиться за формулою:

де b - розмір по довжині виробки, м (b = 1 м).

h - товщина кріплення, м (попередньо вибирається за емпіричною формулою)

σ сж = 151,11 / 1 * 0,25 = 604, кН/м2

Крім того, у перетині кріплення діє згинальний момент Mmax, тоді:

де W-момент опору перерізу, для прямокутного перерізу

.

σ виг = 37,50 / 0,0104 = 3606 кН/м2

Тоді на одній грані перерізу бетонної крепі стискаючі напруги складуть:

На іншій грані кріплення розтягують напруги складуть:

σ сж = 3606 +604 = 4210кПа = 4,21 мПа

σ р = 604-3606 = кПа = 3мПа

Після вибирається марка бетону (В 20) для якої визначається межа міцності на стиск Rb і розтягнення Rbt (див. таблицю 3.2)

Таблиця 3.1. - Опір бетону.

Бетон

В15

В20

В25

В30

В35

Межа міцності на стиск Rb, мПа

8,5

11,5

14,5

17,0

19,5

Межа міцності на розтяг Rbt, мПа

0,75

0,9

1,05

1,2

1,3

Дані напруги порівнюються з граничними відповідно до маркою бетону. Наші отримані значення не перевищують граничні Rb (4,21 <11.5), але перевищують Rbt (3> 0.9)

У нас не виконується умова по граничним розтягуючим напруженням в кріплення, тому найбільш доцільним можна вважати застосування армування бетонної крепі.

Армування бетонного кріплення може проводиться одиночної арматурою (сітка з розмірами 0,2 х0, 2м). Приймається арматура класу А-III, з межею міцності на розтяг (Rs) рівною 340 мПа. Причому, переріз арматури на 1м заздалегідь не відомо (As), тому орієнтовно можна прийняти її рівною 1% від площі.

As = 0,25 * 1 * 0,01 = 0,0025 м2

Після чого знаходимо висоту стиснутої зони бетону за формулою:

Х = 0,0025 * 340/11, 5 * 1 = 0,07, м

Несуча спроможність перетину визначається за формулою:

h 0 - відстань від розтягнутої арматури до найбільш віддаленої точки стиснутої зони h 0 = h - а, м (а - товщина захисного шару розтягнутої арматури а = 0,02-0,03 м)

h 0 = 0,25-0,03 = 0,22 м

кН * м

У випадку якщо отриманий згинальний момент перевищує М max більше, ніж на 10-15% несучу спроможність перетину М1, слід зменшити площу застосовуваної арматури.

Зменшуємо площа перерізу арматури:

As = 0,25 * 1 * 0,01 / 5 = 0,0005

Потім знаходимо висоту стиснутої зони бетону:

Х = 0,0005 * 340/11, 5 * 1 = 0,015

Несуча спроможність у цьому випадку визначається за формулою

М1 = 11,5 * 1 * 0,015 (0,22 / 0,015 / 2) = 37 кН * м

Отриману підсумкову площа арматури слід розділити на п'ять одиниць (так як сітка має розміри 0,2 х0, 2м) тобто п'ять прутків на 1 метр. Потім від площі арматури необхідно перейти до її діаметра і округлити до найближчого цілого числа.

Перетин прутка визначається за формулою:

S = As / 5

S = 0,0005 / 5 = 0,0001

Визначаємо діаметр прутка за формулою:

0,011 м або d = 1,1 см

Приймаємо діаметр прутка 1,2 см.

Ми визначили, що при наших гірничо-геологічних умовах при будівництві гірничої виробки із застосуванням бетонного кріплення необхідно дотримуватися наступних умов:

- Товщина кріплення дорівнює 0,25 м;

- Марка застосовуваного бетону дорівнює В20;

- Армування прутком діаметром 1,2 см.

Монолітна залізобетонна кріплення застосовується для кріплення капітальних горизонтальних і похилих виробок з великим терміном служби і знаходяться поза зоною очисних робіт. Найбільш широко вона застосовується для кріплення виробок приствольних дворів, камер, вузлів сполучення, капітальних виробок і стовбурів.

СПИСОК ісползовать ДЖЕРЕЛ

1. Каретников В.Н., Клейменов В.Б., Нуждіхін А.Г. Кріплення капітальних і підготовчих гірських виробок. М.: Недра, 1989р

2. Буличов Н.С., Фотієва М.М., Стрільців Є. В, Проектування і розрахунок кріплення капітальних виробок. М.: Надра 1986р

3. Судариков А.Е. методические указания к курсовому по дисциплине "Расчет конструкций подземных сооружений".


Додати в блог або на сайт

Цей текст може містити помилки.

Геологія, гідрологія та геодезія | Курсова
74.6кб. | скачати


Схожі роботи:
Розробка методів аналізу деформацій підземних споруд
Технологія будівництва підземних споруд методом опускного колодязя
Неоптолемеевская механіка як механіка ери космосу
Видобуток підземних вод
Розрахунок динаміки підземних вод
Облаштування підземних переходів в м Суми
Ліцензування видобутку підземних вод
Практичне значення та охорона підземних вод
Охорона і раціональне використання поверхневих та підземних вод
© Усі права захищені
написати до нас