Полтавський національний технічний університет
імені Юрія Кондратюка
Кафедра екології
Курсова робота з дисципліни:
Метеорологія і кліматологія
Виконав: студент групи 201-СЕ
Степаненков Г.В.
№ залікової книжки 07113
Перевірив: Ільяш О.Е.
Полтава 2009
ВИЗНАЧЕННЯ МЕТЕОРОЛОГІЧНИХ І КЛІМАТИЧНИХ ФАКТОРІВ
1. Визначення радіаційного балансу діяльного шару землі
1.1 Розрахунок радіаційного балансу
Радіаційний баланс діяльного шару землі R являє собою суму короткохвильової радіації RК та довгохвильової радіації RД і залежить від складових його потоків
R = RK + RД = (SГ + D) · (1 - A) - BЕФ , (1.1)
де SГ - потік сонячної радіації на горизонтальну земну поверхню, Вт/м2, [1,2] ;
D - потік розсіяної радіації на земну поверхню, Вт/м2, [1,2] ;
А - альбедо земної поверхні (див. завдання);
ВЕФ- ефективне випромінювання, узяте зі зворотнім знаком, Вт/м2 .
Величина ВЕФ у свою чергу визначається як
ВЕФ = σ · δз ·[ T34 - TA4· ( 0.61 + 0.05 √е )] , (1.2)
де σ - постійна Стефана-Больцмана, що дорівнює 5.7·10-11 кВт/(м2 К4);
δз - відносна випромінююча властивість земної поверхні чи коефіцієнт випромінювання, δз=0,86;
ТЗ - температура земної поверхні ,0К, що дорівнює ТЗ = 273 + t3 . Температура земної поверхні t3 ,0С, приймається за [1];
Т3 = 273+ 23=296° К;
ТА - середнє значення температури повітря найбільш теплого місяця року, 0К(див анотацію);
ТА=273 + ta=290°К;
е - парціальний тиск водяної пари (пружність водяної пари), гПа, [1,2].
е - потенційний тиск водяної пари (пружність водяної пари), гПа,
Пружність водяної пари визначається за формулою
де φ - відносна вологість повітря φ = 65%; Е - тиск насиченої пари, гПа, Е=20,24 гПа.
(гПа)0,099(Вт/м2)
(Вт/м2) (Вт/м2) (Вт/м2) (Вт/м2) (Вт/м2) (Вт/м2) (Вт/м2) (Вт/м2)Розрахунок величини радіаційного балансу R зводиться в табличну форму - табл.1.1.
Таблиця 1.1
Величинапотоків | Значення потоків, Вт/м2 | |||||||||
Часи доби | ||||||||||
2-3 | 3-4 | 4-5 | 5-6 | 6-7 | 7-8 | 8-9 | 9-10 | 10-11 | 11-12 | |
I | - | - | 17 | 98 | 213 | 360 | 495 | 590 | 677 | 730 |
Sr | - | - | 15 | 56 | 84 | 99 | 112 | 126 | 129 | 133 |
Q | - | - | 32 | 154 | 297 | 459 | 607 | 716 | 806 | 863 |
R | - | - | 28 | 135 | 261 | 403 | 534 | 630 | 709 | 759 |
21-22 | 20-21 | 19-20 | 18-19 | 17-18 | 16-17 | 15-16 | 14-15 | 13-14 | 12-13 | |
Часи доби |
1.2 Побудова діаграми добового розподілу радіаційного балансу
На основі даних таблиці 1.1 будують сумісну діаграму добового розподілу сумарної сонячної радіації та радіаційного балансу у теплий період (липень).По горизонтальній осі відкладають часи доби, а по вертикальній - значення величини Q,(Вт/м2) та величини R , (Вт/м2) - рис1.1.
2.Аналіз теплового режиму атмосфери
2.1 Визначення розподілу температури атмосферного повітря по висоті
Розподіл температури атмосферного повітря по висоті характеризується вертикальним температурним градієнтом, 0С/м ,
γ = - (∆ t /∆ z) ·100 , (2.1.)
де ∆ t = tВРn - tВР n-1 - різниця температур повітря на верхньому (заданому) та нижньому рівнях, 0С;
∆ z - різниця рівнів, м ( див. завдання);
γ - вертикальний температурний градієнт відповідно до кожного рівня висоти ( див. завдання).
Визначення розподілу температури по висоті tВР1¸ tВР5 проводять за формулою
tВРn = - [ ( γ · ∆ z) / 100 ] + tВР n-1 (2.2.)
s
1)
Z1=50м2)
Z2=100м3)
Z3=700м4)
Z4=1300м5)
Z5=2000м2.2 Побудова графіка кривої стратифікації
Користуючись визначеними даними зміни температури повітря з висотою у заданий період року, будують графік кривої стану атмосферного повітря (рис.2.1) – кривої стратифікації. На горизонтальній осі координат відкладають значення температур (0С) ,а на вертикальній – висоту (м).
2.3 Визначення характеру стану атмосфери
Визначення характеру стану атмосфери на заданих рівнях можна проводити, користуючись двома методами.
Ι метод передбачає порівняння динаміки зміни кривої стратифікації (1) з кривою стану суміші повітря, що адіабатично підіймається (2), яка додатково будується на рис.2.1.
Крива (2) характеризується величиною адіабатичного вертикального градієнта
γа = (∆ t /∆ z) ·100 = 1 0С/ 100м (2.3)
1)
Z1=50м2)
Z2=100м3)
Z3=700м4)
Z4=1300м5)
Z5=2000мЗа даними рис.2.1 проводиться аналіз характеру стану атмосфери на кожному рівні:
Z1=50м γ1=0,3 γ1< γa-сильно усталена ;
Z2=100м γ2=0,9 γ2< γa - слабо усталена
Z3=700м γ3=-0,5 γ3< 0- інверсія;
Z4=1300м γ3=-1 γ4<<0 - інверсія;
Z5=2000м γ3=-1,3 γ5>γa- неусталена.
ІІ метод передбачає дослідження стану атмосфери шляхом визначення зміни з висотою величини потенційної температури, ОС,
Θ = tВРп + γа· ( Z / 100 ) . (2.4)
Розподіл потенційної температури з висотою зображений на мал. 2.2., де по горизонтальній вісі координат відкладається температура 0, °С, а по вертикальній - висота z,, м.
До висоти Z4 = 1300 м значення Θ збільшується , шо характеризує стан атмосфери як усталений, вище від 1300м до 2000м значення Θ зменшується, тому стан атмосфери неусталений.
3 Аналіз розподілу водяної пари в атмосфері
3.1 Визначення характеристик вологості повітря
Пружність (парциальний тиск) водяної пари визначається за формулою гПа,
е= φ · Е , (3.1)
де φ – відносна вологість повітря, %, φ = 88%;
Е – тиск насиченої пари, гПа, приймається за додатком 1 відповідно до значення tа(Е=4,15 гПа).
Абсолютна вологість повітря визначається за формулою , г/м3,
а = 217 · е / Т , (3.2)
де Т = 273 + tа , оК .
г/м3Питома вологість повітря , г/кг, визначається як
s = 622 · е / p , (3.3)
де p – атмосферний барометричний тиск приймається рівним 1013.3 гПа.
(г/кг)3.2 Визначення розподілу вологості повітря по висоті
Розподіл вологості по висоті характеризується величинами відносної φz або абсолютної аzвологості на відповідних рівнях
φz = еz / Еz· 100 % , (3.4)
аz = 217 · еz / ТВР , (3.4а)
де еz - парціальний тиск ненасиченого вологого повітря, гПа, на висоті Z визначається за формулою
еz = ео · 10 -z/6300 , (3.5)