Klimatyzacja - praca, Referat z Klimatyzacji, ANALIZA ZNANYCH Z LITERATURY METOD PROGNOZOWANIA TEMPERATURY POWIETRZA


METODY PROGNOZOWANIA TEMPERATURY ŚREDNIEJ POWIETRZA W PRĄDACH OPŁYWOWYCH

  1. Klasyczna metoda Vossa.

  2. Zmodyfikowana metoda Vossa.

  3. Graficzna metoda prognozowania temperatury powietrza kopalnianego (metoda wykresów Ψ i Ω ).

  4. Metoda J. Wacławika.

  5. Metoda Szczerbania, Kremniewa, Żurawlenki.

  6. Metoda PTO-2.

POLITECHNIKA WROCŁAWSKA WROCŁAW 18.01.2001.

WYDZIAŁ GÓRNICZY

REFERAT Z KLIMATYZACJI

TEMAT: METODY PROGNOZOWANIA TEMPERATURY ŚREDNIEJ POWIETRZA W PRĄDACH OPŁYWOWYCH.

WYKONALI :

GRZEGORZ CIELICZKA

MICHAŁ DZIEDZIC

METODY PROGNOZOWANIA TEMPERATURY ŚREDNIEJ POWIETRZA W PRĄDACH OPŁYWOWYCH.

Metody prognozowania warunków klimatycznych w kopalniach podziemnych rozwijane są od kilkudziesięciu lat [11]. Początkowo były one oparte tylko na prawie zachowania energii i teorii ruchu ciepła. Z kolei wobec dużej rozbieżności wyników prognozy z rzeczywistością kopalnianą, konieczne okazało się uwzględnienie zasady zachowania energii i zasady zachowania masy (wilgoci), jak również łącznej teorii ruchu ciepła i masy.

Starsze metody, tj. metody opracowane w bieżącym stuleciu do końca lat sześćdziesiątych oparte są na teorii przewodzenia ciepła w skałach suchych i stosowane są przeważnie do prognozowania wyłącznie temperatury powietrza [12]. W niektórych jednak z tych metod wykorzystuje się pomiary przyrostów wilgotności i temperatury powietrza w wybranych wyrobiskach kopalnianych, by na ich podstawie przewidywać wilgotność powietrza w wyrobiskach badanych, a także uwzględniać wpływ parowania wody na temperaturę powietrza. Metody te jednak okazały się zawodne.

W wyrobiskach górniczych i w otaczających skałach zachodzi wiele złożonych procesów wymiany ciepła i masy, które powodują zmiany parametrów przepływającego powietrza [4]. Transport ciepła odbywa się w sposób jawny drogą przewodnictwa cieplnego, konwekcji i promieniowania oraz w sposób niejawny drogą parowania wody, w szczególności wilgoci. Temperatura powietrza w wyrobiskach górniczych różni się zwykle od temperatury skał. W związku z tą różnicą zachodzi wymiana ciepła pomiędzy skałami a przepływającym powietrzem. Temperatury obu ośrodków zależą m. in. od czasu, proces przepływu powietrza jest więc nieustalony.

W szybach i w wyrobiskach nachylonych następuje sprężanie powietrza w polu siły ciężkości przy ruchu w dół lub rozprężanie przy wznoszącym ruchu powietrza. Wzrostowi ciśnienia towarzyszy zwiększenie temperatury. Przy obniżaniu ciśnienia temperatura powietrza maleje.

W wyrobiskach kopalnianych wydziela się ciepło od źródeł miejscowych, takich jak m. in. procesy utleniania, pracujące maszyny i urządzenia. Następują zmiany wilgotności powietrza.

W celu określenia wpływu temperatury skał na temperaturę powietrza wyznacza się pole temperatury skał wokół wyrobiska kopalnianego. W ogólności temperatura skał jest funkcją położenia i czasu:

0x01 graphic
(1)

Zakłada się, że skały są ciałem stałym, w którym transport energii odbywa się drogą przewodnictwa cieplnego. Przyjmując, że skały stanowią ośrodek izotropowy ze względu na zjawisko przewodnictwa cieplnego, to wektor gęstości strumienia cieplnego q jest prostopadły do powierzchni izotermicznej, a w kierunku stycznym do izoterm nie ma przepływu ciepła.

Z równania bilansu energii przeprowadzonego dla dowolnego obszaru przestrzennego w jednorodnym i izotropowym masywie skalnym otrzymuje się równanie różniczkowe przewodnictwa cieplnego, które ma postać:

0x01 graphic
(2)

W obliczeniach dotyczących temperatury skał korzysta się z cylindrycznego układu współrzędnych r, ϕ, s, gdzie:

r - promień wodzący,

ϕ - kąt kierunkowy,

s - oś wyrobiska.

Równanie przewodnictwa cieplnego ma postać:

0x01 graphic
(3)

Na skutek wymiany ciepła między skałami a przepływającym powietrzem następuje zmiana temperatury skał bezpośrednio otaczających wyrobisko. Wartość bezwzględna przyrostu temperatury skał na jednostkę długości w kierunku promienia r ma rząd wielkości 1°C/m i jest znacznie większa od wartości bezwzględnej zmiany temperatury na jednostkę długości w kierunku równoległym do osi wyrobiska, która na ogół nie przekracza 10-2 K/m.

Wobec tego w równaniu przewodnictwa cieplnego pomija się ciepło przewodzone w kierunku równoległym do osi wyrobiska. Po tych uproszczeniach równanie bilansu ciepła w skałach wyrażone we współrzędnych cylindrycznych przyjmuje postać:

0x01 graphic
(4)

W rozważaniach nad zmianami temperatury powietrza kopalnianego korzysta się z pojęcia uogólnionej entalpii właściwej powietrza:

0x01 graphic
, J/kg (5)

gdzie:

h - entalpia uogólniona właściwa, J/kg,

cp - ciepło właściwe przy stałym ciśnieniu, J/(kg K),

z - wysokość niwelacyjna, m.

W wyrażeniu (5) prócz klasycznej termodynamicznej entalpii występują składniki związane z energią potencjalną i kinetyczną. W rozważanym przypadku zmiana uogólnionej entalpii przypadająca na jednostkę długości wyrobiska równa jest strumieniowi ciepła dopływającego do jednostki masy powietrza, na tym odcinku:

0x01 graphic
, J/(kg m) (6)

gdzie:

qs - strumień ciepła przekazywany od skał do powietrza przez powierzchnię calizny skalnej na odcinku wyrobiska o jednostkowej długości, W/m,

qd - strumień ciepła przenoszony od innych źródeł ciepła (prócz ciepła od skał) do powietrza, na odcinku wyrobiska o jednostkowej długości, W/m.

Pomijając w zależności (6) zmiany energii kinetycznej jako małe w porównaniu z innymi składnikami, po przekształceniu otrzymuje się równanie energii przepływającego powietrza:

0x01 graphic
, K/m (7)

Wyrażenia (4) i (7) stanowią układ równań różniczkowych na temperatury skał i powietrza. W celu wyznaczenia funkcji spełniających te równania i odpowiadających konkretnemu zagadnieniu fizycznemu należy sformułować warunki graniczne (tj. warunek początkowy i brzegowy) dla równania przewodnictwa cieplnego (4) oraz warunek początkowy dla równania energii (7).

Temperaturę skał 0x01 graphic
i temperaturę powietrza T(s,t) przedstawia się w postaci dwóch składników:

0x01 graphic
(8)

Pierwsze składniki prawych stron wzorów (8) noszą nazwy średniej temperatury skał 0x01 graphic
i średniej temperatury powietrza Ts. Wyznacza się je przy założeniu, że temperatura powietrza na wlocie do rozważanego wyrobiska lub do kopalni jest stała w czasie.

Zagadnienie brzegowo- początkowe do wyznaczenia funkcji 0x01 graphic
i Ts obejmuje:

0x01 graphic
(9)

0x01 graphic
dla t = 0 (10)

0x01 graphic
dla r = ro (11)

0x01 graphic
(12)

To = Tos dla s = 0 (13)

Drugie składniki prawych stron wyrażeń (8) noszą nazwę okresowych zmian temperatur skał 0x01 graphic
i powietrza Tr. Zagadnienie brzegowo - początkowe do wyznaczenia funkcji 0x01 graphic
i Tr obejmuje:

0x01 graphic
(14)

0x01 graphic
dla t = 0 (15)

0x01 graphic
(16)

0x01 graphic
(17)

0x01 graphic
(18)

Jak widać rozwiązanie wymaga wielu założeń upraszczających i danych początkowych. W związku z tym wiele metod jest często przystosowanych do konkretnych warunków wybranych zagłębi czy nawet kopalń.

W niniejszym rozdziale zostaną przedstawione i omówione wybrane metody prognozowania temperatury powietrza.

1. Klasyczna Metoda Vossa

Metoda J. Vossa została opracowana w latach 1965÷1973 przez Zespół Badawczy ds Klimatyzacji Kopalń w RFN [5].

Żeby móc ocenić wpływ poszczególnych czynników na kształtowanie się warunków klimatycznych jak również przewidzieć jak będą się one kształtowały w projektowanym wyrobisku górniczym, konieczne jest przeprowadzenie analizy obliczeniowej temperatury powietrza w tym wyrobisku [6].

W wyrobiskach górniczych bez dodatkowych źródeł ciepła, całkowita ilość ciepła, którą powietrze przyjmuje jest odebrana z masywu skalnego. Ciepło to przedostaje się z głębi masywu skalnego na jego powierzchnię w skałach litych w wyniku przewodnictwa, a w występujących porach i szczelinach skał przez dyfuzję lub przepływ wody albo pary wodnej.

Całkowity strumień ciepła przepływający przez powierzchnię masywu skalnego przedstawia się następująco:

0x01 graphic
, kJ/(m2 h K) (19)

gdzie:

λeq - ekwiwalentny współczynnik przewodnictwa cieplnego, W/m K,

tgo - temperatura nienaruszonego robotami górniczymi masywu skalnego, °C,

tp - temperatura przepływającego powietrza, °C

ro - promień równoważny poprzecznego przekroju wyrobiska, m,

Ki - wielkość bezwymiarowa, przedstawiająca zmiany ilości przyjmowanego ciepła od skał będąca funkcją liczby Fouriera F0 i liczby Biota Bi.

Z całkowitego strumienia ciepła, jak wykazują pomiary, tylko jego część powoduje przyrost temperatury powietrza, mierzonej termometrem suchym, a mianowicie:

0x01 graphic
, kJ/(m2 h K) (20)

przy czym 0x01 graphic
- współczynnik ciepła konwekcyjnego, wyrażający stosunek suchego do całkowitego przejmowania ciepła przez powietrze, a zatem:

0x01 graphic
(21)

Wartości współczynnika ciepła konwekcyjnego εs przedstawiono w tabeli 5.1.

Tabela 5.1.Wartości współczynnika ciepła konwekcyjnego

Współczynnik ciepła konwekcyjnego εs

Nazwa wyrobiska

0x01 graphic

Chodniki kamienne

0,35

Ściany prowadzone na zawał gdy są zainstalowane MK

0,25

Ściany prowadzone na zawał gdy nie są zainstalowane MK

0,15

Temperaturę powietrza mierzoną termometrem suchym wyznacza się z zależności:

0x01 graphic
, °C (22)

gdzie:

B - obwód poprzecznego przekroju wyrobiska, m

0x01 graphic
- wydatek masowy powietrza, kg/h

L - długość wyrobiska, m

tpso - temperatura powietrza na wlocie do wyrobiska mierzona termometrem suchym, °C.

cpp - ciepło właściwe powietrza przy stałym ciśnieniu, równe 1,005 kJ/(kg K),

Wartości liczbowe ekwiwalentnego współczynnika przewodnictwa cieplnego λeq zostały określone poprzez staranne i obszerne pomiary w drogach powietrznych wyrobisk górniczych. Okazało się, że z powodu dość dużego rozrzutu wartości pomiarowych celowym działaniem jest rozłączenie wartości λeq dla chodników kamiennych, chodników przyścianowych oraz dla wyrobisk ścianowych.

O ile tylko przy określaniu warunków klimatycznych nie są wymagane bardzo dokładne wartości dla ekwiwalentnego współczynnika przewodnictwa cieplnego można przyjmować dla kopalni węgla następujące średnie wartości:

  1. chodniki kamienne, λeq = 5,8 W/mK

  2. chodniki przyścianowe:

w przypadku urabiania, λeq = 8,1 W/mK

w przypadku nieurabiania, λeq = 2,3 W/mK

  1. wyrobiska ścianowe:

z zawałem, λeq = 7 W/mK

z podsadzką, λeq = 3 W/mK

Wielkość bezwymiarowa Ki = f(F0, Bi) jest wielkością charakteryzującą ochłodzenie się górotworu. Ochłodzenie to zależy naturalnie od całkowitego przepływu ciepła z górotworu q a nie tylko od jego strumienia częściowego qs. Z tego też powodu musi się liczbę Fouriera F0 związać z ekwiwalentnym współczynnikiem przewodnictwa temperaturowego aeq w następujący sposób:

0x01 graphic
(23)

i równocześnie

0x01 graphic
(24)

gdzie:

cs - ciepło właściwe skał, J/(kg K),

ρs - gęstość skał, kg/m3.

Tabela 5.2. Wartości ciepła właściwego i gęstości skał.

Ciepło właściwe skał i ich gęstość cs i ρs

Rodzaj skał

0x01 graphic
,

kg/m3

cs ,

J/(kg K)

Zlepieńce i żwirowce

2200

961

Piaskowiec gruboziarnisty

2400

696

Piaskowiec drobnoziarnisty

2500

705

Łupek piaszczysty

2550

850

Łupek ilasty

2600

1000

Łupek miedzionośny

2650

954

Dolomit

2520

840

Węgiel kamienny

1300

439

Liczba Biota jest określona następującym związkiem:

0x01 graphic
(25)

a współczynnik przejmowania ciepła αk można określić z następującego empirycznego wzoru

0x01 graphic
, W/m2 K (26)

gdzie:

w - prędkość przepływu powietrza, m/s,

ξ/ξo - stosunek spadków ciśnienia w rozpatrywanym wyrobisku do spadku ciśnienia w rurowym przewodzie gładkim, którego promień jest równy hydraulicznemu promieniowi wyrobiska.

Dla wyrobisk chodnikowych wartość tego stosunku zawarta jest średnio w granicach 3,5 ÷ 5,8, natomiast dla wyrobisk ścianowych średnio w granicach 4,7 ÷ 9,3.

Wartości dla Ki można odczytać z diagramu Ki = f (F0, Bi) przedstawionego na rys. 5.1 Przy określaniu bezwymiarowego współczynnika Ki można korzystać również z uproszczonego wzoru wg Batzel'a w następującej postaci:

0x01 graphic
(27)

0x08 graphic

Rys. 5.1 Nomogram dla funkcji Ki = f (F0, Bi)

Ciśnienie barometryczne panujące w chodniku można określić z tzw. barometrycznej formuły wysokości w postaci:

0x01 graphic
, Pa (28)

po - ciśnienie powietrza na zrębie szybu wdechowego, Pa,

Strumień masy powietrza określony jest wzorem:

0x01 graphic
, kg/h (29)

przy czym gęstość powietrza suchego można wyznaczyć z zależności:

0x01 graphic
, kg/m3 (30)

gdzie:

Ra - indywidualna stała gazowa równa 287,04 J/(kg K)

T - temperatura powietrza mierzona termometrem suchym wyrażona w stopniach °K

T = 273,15 + tpso (31)

gdzie tpso oznacza temperaturę powietrza na wlocie do wyrobiska mierzona termometrem suchym.

Natomiast prędkość przepływu powietrza przez wyrobisko wynosi:

0x01 graphic
, m/s (32)

gdzie A - pole przekroju poprzecznego wyrobiska, m2.

Równoważny promień poprzecznego przekroju wyrobiska określony jest wzorem

0x01 graphic
(33)

gdzie B oznacza obwód poprzecznego przekroju wyrobiska.

Temperaturę nienaruszonego robotami górniczymi masywu skalnego można wyznaczyć za pomocą następującego wzoru:

0x01 graphic
(34)

gdzie:

t0 - temperatura skał na głębokości 0x01 graphic
,°C,

0x01 graphic
- stopień geotermiczny, m/°C,

H - głębokość zalegania pokładu, m.

Przyrost temperatury powietrza w chodniku mierzony termometrem suchym wynosi:

0x01 graphic
, °C (35)

Przyrost entalpii powietrza wynosi:

0x01 graphic
, kJ/kg (36)

Przyrost zawartości pary wodnej w powietrzu płynącym wyrobiskiem określa się z następującego wzoru:

0x01 graphic
, kg/kg (37)

gdzie:

rb - ciepło parowania wody, równe 2486 kJ/kg

Zawartość pary wodnej w powietrzu wynosi:

0x01 graphic
, kg/kg (38)

zaś temperaturę powietrza na termometrze wilgotnym tpwz odczytujemy z diagramu h - x w sposób pokazany na rys. 5.2 , p. 2.

0x08 graphic

Rys. 5.2 Wykres h - x.

Entalpię powietrza na wlocie do wyrobiska można odczytać z wykresu h - x przedstawionego na rys. 5.2 --> [Author:FR] lub obliczyć na podstawie następującego wzoru:

0x01 graphic
, kJ/kg (39)

gdzie:

Cpp - ciepło właściwe powietrza przy stałym ciśnieniu, równe 1,005 kJ/(kg K),

xo - zawartość pary wodnej w powietrzu, kg/kg

0x01 graphic
(40)

Prężność pary wodnej pp otrzymuje się z empirycznego wzoru Sprung'a

0x01 graphic
, hPa (41)

gdzie:

tpwo - temperatura powietrza na wlocie do wyrobiska mierzona termometrem

wilgotnym, °C,

tpso - temperatura powietrza na wlocie do wyrobiska mierzona termometrem suchym,

°C,

p - bezwzględne ciśnienie statyczne powietrza, hPa,

pw - prężność pary wodnej nasyconej w powietrzu (określana na podstawie tablicy

5.3), hPa.

Wartość entalpii powietrza na wylocie z wyrobiska wyznacza się z następującego wzoru:

0x01 graphic
, kJ/kg (42)

Jak wykazały przeprowadzone badania, tylko od 25% do 35% z wszystkich zainstalowanych urządzeń elektrycznych w wyrobiskach górniczych działa równocześnie, ze względu na to, że nie wszystkie urządzenia działają cały czas.

W rezultacie tarcia występującego w maszynach energia zamieniona na pracę przechodzi z powrotem na ciepło. Tylko w przypadku zmiany jej na inne postacie energii (np. potencjalnej przy transporcie w górę) część energii pobranej przez urządzenie nie przechodzi w ciepło.

Ponieważ metoda ta nie uwzględnia nachylenia wyrobiska w związku z tym nie jest wykonywana żadna praca podnoszenia. Dlatego wartość oddawanego ciepła przez urządzenie elektryczne można przyjąć jako:

0x01 graphic
, kJ/h (43)

gdzie:

Nz --> [Author:FR] - moc zainstalowanych silników elektrycznych, kW,

Zatem dodatkowy przyrost entalpii powietrza w wyrobisku powstały w wyniku pracujących tam urządzeń elektrycznych wyraża się następująco:

0x01 graphic
, kJ/kg (44)

Liczne pomiary wykazały również, że tylko 90% uprzednio określonej ilości ciepła będzie oddziaływać, ponieważ ten dodatkowy przyrost entalpii między górotworem a prądem przepływającego powietrza powoduje, że zostaje przytamowany strumień cieplny z górotworu .

Dlatego entalpia powietrza na wylocie z wyrobiska przy uwzględnieniu pracy urządzeń elektrycznych wynosi:

0x01 graphic
, kJ/kg (45)

Z tego, spowodowanego przez pracę urządzeń elektrycznych przyrostu entalpii, tylko około 20% występuje w postaci suchego przekazania ciepła, a z tego, z powodu przyhamowania strumienia ciepła z górotworu, oddziałuje na temperaturę suchą tylko 80%.

Zatem przyrost temperatury powietrza na termometrze suchym wynikający z umieszczenia pracujących urządzeń elektrycznych wynosi:

0x01 graphic
, °C (46)

Wartość temperatury powietrza na termometrze suchym na wylocie wyrobiska, uwzględniająca ciepło pochodzące od pracujących urządzeń elektrycznych, będzie wynosić:

0x01 graphic
, °C (47)

Natomiast wartość temperatury powietrza na termometrze wilgotnym na wylocie z wyrobiska uwzględniająca ciepło pochodzące z pracujących urządzeń elektrycznych określa się z diagramu h - x w sposób pokazany na rys. 5.2, p.7

Tabela 5.3 Prężność pary wodnej nasyconej nad wodą w zależności od temperatury

Tempe-ratura

oC

Prężność pary wodnej nasyconej w powietrzu pw, hPa

0,0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

-4

-3

-2

-1

-0

0

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

4,55

4,89

5,28

5,68

6,11

6,11

6,56

7,05

7,57

8,13

8,72

9,35

10,01

10,72

11,48

12,27

13,12

14,01

14,97

15,97

17,04

18,17

19,37

20,63

21,96

23,37

24,87

26,43

28,09

29,84

31,68

33,61

35,65

37,80

40,05

42,43

44,93

47,56

50,31

53,20

56,24

59,43

62,77

66,27

69,95

4,51

4,87

5,24

5,64

6,07

6,16

6,61

7,11

7,63

8,19

8,79

9,41

10,08

10,80

11,55

12,36

13,21

14,11

15,07

16,08

17,16

18,29

19,49

20,76

22,11

23,52

25,01

26,60

28,25

30,01

31,87

33,81

35,87

38,03

40,29

42,68

45,19

47,83

50,60

53,51

56,55

59,76

63,11

66,63

70,32

4,48

4,83

5,20

5,60

6,03

6,19

6,67

7,16

7,68

8,24

8,84

9,48

10,15

10,87

11,63

12,44

13,29

14,20

15,17

16,19

17,27

18,41

19,61

20,89

22,24

23,67

25,17

26,76

28,43

30,19

32,05

34,01

36,08

38,24

40,52

42,92

45,44

48,09

50,88

53,80

56,87

60,08

63,45

66,99

70,99

4,44

4,79

5,16

5,56

5,97

6,24

6,71

7,21

7,73

8,31

8,91

9,53

10,23

10,95

11,71

12,52

13,39

14,29

15,27

16,29

17,37

18,52

19,73

21,03

22,39

23,81

25,32

26,92

28,60

30,37

32,24

34,21

36,28

38,47

40,76

43,17

45,71

48,37

51,16

54,11

57,17

60,41

63,80

67,35

71,07

4,41

4,76

5,12

5,52

5,93

6,29

9,76

7,25

7,79

8,36

8,96

9,61

10,29

11,01

11,79

12,61

13,47

14,39

15,36

16,40

17,49

18,64

19,87

21,16

22,52

23,96

25,48

27,08

28,77

30,56

32,44

34,41

36,49

38,69

41,00

43,41

45,96

48,64

51,45

54,40

57,49

60,75

64,15

67,72

71,45

4,37

4,72

5,08

5,47

5,89

6,33

6,80

7,31

7,85

8,43

9,03

9,68

10,36

11,09

11,87

12,69

13,56

14,48

15,47

16,51

17,61

18,76

19,99

21,29

22,67

24,11

25,64

27,25

28,95

30,75

32,63

34,61

36,71

38,92

41,23

43,67

46,23

48,92

51,73

54,71

57,81

61,08

64,49

68,08

71,84

4,35

4,68

5,04

5,44

5,84

6,37

6,85

7,36

7,91

8,48

9,09

9,75

10,43

11,17

11,95

12,77

13,65

14,59

15,57

16,61

17,72

18,88

20,12

21,43

22,80

24,25

25,08

27,41

29,12

30,92

32,83

34,83

36,93

39,15

41,47

43,92

46,49

49,19

52,03

55,01

58,13

61,41

64,85

68,45

72,23

4,31

4,65

5,01

5,39

5,80

6,43

6,91

7,41

7,96

8,53

9,16

9,81

10,51

11,24

12,03

12,87

13,75

14,68

15,67

16,72

17,83

19,00

20,24

21,56

22,95

24,41

25,95

27,59

29,31

31,11

33,01

35,03

37,15

39,37

41,71

44,17

46,75

49,47

52,32

55,32

58,45

61,75

65,20

68,83

72,61

4,28

4,61

4,97

5,36

5,76

6,47

6,96

7,47

8,01

8,60

9,21

9,88

10,57

11,32

12,11

12,95

13,84

14,77

15,77

16,83

17,95

19,12

20,37

21,69

23,09

24,56

26,11

27,75

29,48

31,29

33,21

35,24

37,36

39,60

41,95

44,43

47,01

49,75

52,61

55,63

58,77

62,08

65,56

69,19

73,00

4,24

4,59

4,93

5,32

5,72

6,52

7,00

7,52

8,08

8,65

9,28

9,95

10,65

11,40

12,19

13,04

13,93

14,87

15,88

16,93

18,05

19,25

20,15

21,83

23,23

24,71

26,27

27,92

29,65

31,48

33,41

35,44

37,57

39,83

42,19

44,68

47,28

50,03

52,91

55,93

59,11

62,43

65,91

69,56

73,39

40

41

42

43

44

45

46

47

48

49

73,79

77,81

82,03

86,43

91,04

95,87

100,89

106,17

111,67

117,41

74,17

78,23

82,45

86,88

91,52

96,36

101,41

106,71

112,23

118,00

74,57

78,64

82,89

87,33

91,99

96,85

101,93

107,25

112,80

118,59

74,97

79,05

83,32

87,79

92,47

97,35

102,45

107,79

113,37

119,17

75,37

79,47

83,76

88,25

92,95

97,85

102,97

108,33

113,93

119,79

75,77

79,89

84,20

88,71

93,43

98,36

103,51

108,89

114,51

120,37

76,17

80,32

84,64

89,17

93,91

98,85

104,04

109,44

115,08

120,99

76,59

80,73

85,08

89,64

94,40

99,36

104,57

109,99

115,65

121,57

76,99

81,16

85,53

90,11

94,88

99,88

105,09

110,55

116,24

122,19

77,40

81,59

85,97

90,57

95,37

100,39

105,63

111,11

116,83

122,80

2. Zmodyfikowana metoda Vossa

W oryginalnej metodzie Vossa wychodzi się z podstawowych równań ruchu ciepła i masy. Wyprowadzona zależność nie uwzględnia jednak zmian temperatury powietrza wynikających z nachylenia wyrobiska. Dlatego też zmodyfikowano tę metodę tak, aby mogła służyć do prognozowania temperatury w wyrobiskach nachylonych [7].

Zależności opisujące metodę prognozowania temperatury powietrza kopalnianego

mierzonej termometrem suchym

Temperaturę suchą powietrza na końcu wyrobiska wyznacza się ze wzoru:

0x01 graphic
(48)

przy czym :

0x01 graphic
(49)

gdzie:

0x01 graphic
- średnia temperatura pierwotna skał otaczających wyrobisko, °C, przy czym

0x01 graphic
(50)

0x01 graphic
- odpowiednio temperatura pierwotna skał w przekroju dopływowym i

wypływu, oC, liczona ze wzoru:

0x01 graphic
(51)

gdzie:

t0 - temperatura skał na głębokości 0x01 graphic
,°C,

0x01 graphic
- stopień geotermiczny, m/°C,

ro - promień równoważny wyrobiska, m, przy czym

ro = 0x01 graphic
(52)

0x01 graphic
- współczynnik ciepła konwekcyjnego ogrzewającego powietrze suche,

równy ilorazowi ciepła konwekcyjnego do ciepła całkowitego,

0x01 graphic
-ekwiwalentny współczynnik przewodzenia ciepła w wilgotnym masywie skalnym, W/(m K)

A - pole przekroju wyrobiska, m2,

B - obwód wyrobiska, m.,

L - długość wyrobiska, m.,

cpa- pojemność cieplna powietrza suchego,

cpa = 1005 J/(kgK)

0x01 graphic
- strumień masy powietrza suchego, kg/s,

0x01 graphic
(53)

0x01 graphic
- prędkość średnia powietrza w przekroju, dla którego wyznacza się strumień

masy powietrza, m/s,

0x01 graphic
- odpowiednio wysokości niwelacyjne przekroju dopływu i wypływu

wyrobiska, dla którego wykonuje się prognozę temperatury powietrza, m,

0x01 graphic
- zagęszczenie strumienia od zewnętrznych źródeł ciepła,

W/m, przy czym:

0x01 graphic
(54)

Nz - moc zewnętrznych źródeł ciepła, kW,

czs- współczynnik określający jaka część energii z zewnętrznych źródeł wpływa

na podwyższenie temperatury mierzonej termometrem suchym.

Liczbę Fouriera wyznacza się ze wzoru:

0x01 graphic
(55)

gdzie:

0x01 graphic
- czas przewietrzania wyrobiska, s,

0x01 graphic
- ekwiwalentny współczynnik wyrównywania temperatury, m2/s,

równy:

0x01 graphic
(56)

cs - pojemność cieplna skał, J/(kg K),

0x01 graphic
- gęstość pozorna skał otaczających wyrobisko, kg/m3,

Liczba Biota dana wzorem:

0x01 graphic
(57)

0x01 graphic
- współczynnik przejmowania ciepła z górotworu, W/(m2K),

0x01 graphic
(58)

Liczba Kirpiczewa charakteryzująca ochłodzenie się górotworu wyznaczana jest z monogramów lub wzorów empirycznych jako funkcja liczb Fouriera i Biota, przy czym

jeżeli: Fo ≥ 1

0x01 graphic
(59)

dla <1

0x01 graphic
(60)

natomiast jeżeli: 0,001 ≤ Fo <1

0x01 graphic
(61)

gdzie:

0x01 graphic
(62)

0x01 graphic
(63)

Prognozowanie temperatury powietrza mierzonej termometrem wilgotnym w wyrobisku górniczym

Prognozowanie temperatury powietrza mierzonej termometrem wilgotnym w wyrobisku górniczym przeprowadzić możemy dopiero po prognostycznym wyznaczeniu temperatury mierzonej termometrem suchym.

Temperaturę wilgotną w wyrobisku górniczym wyznaczamy ze wzoru:

0x01 graphic
(64)

gdzie:

xd - stopień zawilżenia powietrza na dopływie wyrobiska, kg/kg,

xw - stopień zawilżenia powietrza na wypływie z wyrobiska,

kg/kg,

rb - ciepło parowania wody w termometrze 0°C;

rb = 2500000 J/kg,

cpw - pojemność cieplna pary wodnej przy stałym ciśnieniu,

cpw = 1927 J/(kg K).

Znając temperaturę mierzoną termometrem suchym 0x01 graphic
i stopień zawilżenia 0x01 graphic
na końcu wyrobiska możemy obliczyć temperaturę mierzoną termometrem wilgotnym na końcu wyrobiska, która jest niewiadomą w następującym równaniu:

0x01 graphic
(65)

Występującą w tym równaniu wartość ciśnienia statycznego, bezwzględnego 0x01 graphic
na końcu wyrobiska można określić z przybliżonego wzoru

0x01 graphic
(66)

gdzie:

0x01 graphic
- liczba oporu wyrobiska,

0x01 graphic
- gęstość powietrza dla warunków normalnych; 0x01 graphic
=1.20 kg/m3.

Równanie (34) rozwiązuje się w sposób przybliżony uzyskując w wyniku temperaturę powietrza mierzoną termometrem wilgotnym 0x01 graphic
na końcu prognozowanego wyrobiska.

3. Graficzna metoda prognozowania temperatury powietrza kopalnianego

W dotychczasowych metodach prognozowania temperatury powietrza kopalnianego stosowane są skomplikowane, trudne do stosowania w praktyce wzory matematyczne.

W 1969 r. opracowano w Głównym Instytucie Górnictwa podstawy teoretyczne metody wykreślnej prognostycznego wyznaczania temperatury powietrza w wyrobiskach górniczych. Na podstawie tego opracowania wykonano w 1970 r. w Głównym Biurze Studiów i Projektów Górniczych „Nomogramy do prognozowania warunków klimatycznych i obliczania parametrów klimatu”. Metoda wykreślna jest łatwiejsza i zarazem mniej pracochłonna.

Model bocznicy sieci wentylacyjnej stosowany w metodzie wykreślnej

Dla prognozowania temperatury powietrza w kopalni dzieli się bocznice sieci wentylacyjnej na odcinki. Zakłada się, że znana jest temperatura powietrza wpływającego do pierwszego z takich odcinków, a przekrój wypływu ostatniego z nich pokrywa się z miejscem w kopalni, którego temperatura podlega prognozie. Jeżeli ciąg takich odcinków zawiera węzły, do którego wpływa powietrze z co najmniej dwu bocznic, to temperaturę powietrza wypływającego z niego oblicza się jako średnią ważoną względem mas strumieni powietrza wpływającego.

0x08 graphic

λ = idem

a = idem

q = idem

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
(d) (w)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

Zi Qi ro

0x08 graphic
0x08 graphic
0x08 graphic
β

0x08 graphic
x=0 x xw

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
x

0x08 graphic
0x01 graphic
p

0x08 graphic
0x08 graphic
td 0x01 graphic
tw

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
cp=idem

0x08 graphic
0x08 graphic

0x08 graphic

k(x) 1m2

0x08 graphic
0x08 graphic
L

Rys. 5.3 Model odcinka bocznicy kopalnianej sieci wentylacyjnej.

d,w - przekroje odpowiednio - dopływu i wypływu powietrza; td, tw - temperatury powietrza w przekrojach - odpowiednio dopływu i wypływu; 0x01 graphic
- strumień masy powietrza; x - współrzędna odległościowa; xw - współrzędna odległościowa przekroju wypływu, równa długości L odcinka bocznicy; ro -promień równoważny; cp - pojemność cieplna właściwa powietrza (pod stałym ciśnieniem); Zi-i - to dodatkowe źródło ciepła; Qi - strumień ciepła ze źródła Zi; 0x01 graphic
- zagęszczenie strumienia ciepła dopływającego z powierzchni górotworu; k(x) - płaszczyzna prostopadła do osi x w punkcie osi o współrzędnej x; λ - współczynnik przewodzenia ciepła skał; a - współczynnik wyrównania temperatury skał; q - gęstość masy skał; p-p - poziom

Stosowany w metodzie wykreślnej model jednego ze wspomnianych odcinków przedstawiono na rys. 5.3. Odcinek ten ograniczony jest przekrojami dopływu d i wypływu w powietrza. Odcinek modelu jest prostoliniowy, a jego oś x zwrócona jest w stronę strzałki wyznaczającej kierunek przepływu powietrza. Oś x jest nachylona do poziomu p-p pod kątem β, równym średniemu kątowi nachylenia odcinka rzeczywistego. Znak β jest dodatni dla wyrobisk opadających w kierunku osi x. Z osią x wiąże się współrzędna odległościowa x. Współrzędna x przekroju dopływu d jest równa zeru, a współrzędna x przekroju wypływu w, oznaczona przez xw, równa jest długości L odcinka rzeczywistego.

Odsłonięta powierzchnia górotworu ma w modelu kształt powierzchni walca. Promień tego walca, zwany równoważnym promieniem ro, wyznacza się za pomocą średniego obwodu B przekroju poprzecznego w odcinku rzeczywistym wyrobiska stosując wzór:

0x01 graphic
(67)

Skały górotworu, otaczające model odcinka wyrobiska (rys.5.3), są jednorodne oraz izotropowe i rozciągają się do nieskończoności w kierunkach prostopadłych do osi x. Charakteryzują je niezależne od czasu: współczynnik λ przewodzenia ciepła, współczynnik a wyrównywania temperatury oraz gęstość (masy) ρs. Współczynnik λ, a oraz gęstość ρs są to średnie parametry analogicznych wielkości skał rzeczywistych, zalegających bliskie otoczenie rozpatrywanego odcinka wyrobiska górniczego w kopalni (obszar skalny sięgający pod powierzchnię odkrytą górotworu w wyrobisku na głębokość ok. 30 m, położony między dwiema płaszczyznami przekrojów: dopływu i wypływu omawianego odcinka wyrobiska).

Niech symbol τ oznacza przedział czasu, którego początek odnosi się do chwili uzyskania wentylacji opływowej w odcinku bocznicy, koniec zaś do chwili, w której temperatura powietrza podlega prognozie.

Zakłada się, że w tym przedziale czasu przez model odcinka przepływa powietrze o niezmiennym strumieniu masy 0x01 graphic
równym średniemu (w tym przedziale) strumieniowi masy powietrza w rzeczywistym odcinku. Zakłada się również, że w modelu właściwa pojemność cieplna cp powietrza (pod stałym ciśnieniem) jest niezmienna (cp=idem).

Z dowolnym punktem osi x (mającym współrzędną x) modelu wiąże się płaszczyznę k(x), przechodzącą przez ten punkt oraz prostopadłą do osi x. W chwili początkowej, w której uruchomiono wentylację opływową (dla τ = 0), można skałom modelu położonym w płaszczyźnie k(x) przypisać temperaturę 0x01 graphic
, którą nazwano temperaturą pierwotną

0x01 graphic
(68)

gdzie

0x01 graphic
- temperatura pierwotna skał w płaszczyźnie przekroju dopływu powietrza,

0x01 graphic
- gradient geotermiczny dla obszaru, pod którym znajduje się odcinek wyrobiska

górniczego.

Temperatura tw powietrza w przekroju wypływu w zależy od temperatury td powietrza w przekroju dopływu (d) oraz od ciepła dostarczonego z górotworu oraz źródeł dodatkowych Zi, do powietrza w rozpatrywanym odcinku wyrobiska, gdzie i = 1, 2...n.

Zagęszczenie 0x01 graphic
strumienia ciepła dopływającego z powierzchni górotworu, związanej ze współrzędną x, tzn. ciepło dopływające w jednostce czasu do powietrza z powierzchni odkrytej górotworu o długości Δx = 1m, odległej o x od przekroju dopływu (d) można wyznaczyć ze wzoru

0x01 graphic
(69)

gdzie

Ki - liczba Kirpiczewa,

tx - temperatura powietrza w przekroju położonym w płaszczyźnie k(x).

Zagęszczenie 0x01 graphic
można również wyrazić za pomocą wzoru Newtona.

0x01 graphic
(70)

gdzie

α - współczynnik wnikania ciepła z powierzchni górotworu do powietrza; w modelu

przyjmuje się α = idem,

0x01 graphic
- temperatura odkrytej powierzchni górotworu w płaszczyźnie k(x).

Zaznaczone w modelu dodatkowe źródła ciepła Zi, dla i = 1, 2 ... n, ilustrują rzeczywiste dodatkowe źródła czynne w wyrobisku. Do takich źródeł ciepła zalicza się np. maszyny i urządzenia górnicze, urobek zwałowany i transportowany lub utleniające się powierzchnie węgla. Zakłada się, że niektóre źródła Zi nie zależą od temperatury powietrza w wyrobisku, inne natomiast zależą tylko od średniej - szacunkowej temperatury tego powietrza. Strumień ciepła Qi, czyli ciepło oddane w jednostce czasu z i - tego dodatkowego źródła ciepła Zi, równe jest strumieniowi i - tego rzeczywistego źródła ciepła, czynnego w wyrobisku.

Korzystając z zależności

0x01 graphic
(71)

w modelu i metodzie wykreślnej posługujemy się zagęszczeniem strumienia ciepła z i - tego dodatkowego źródła ciepła, odniesionym do 1 m długości odcinka, które jest określone wzorem:

0x01 graphic
(72)

W metodzie wykreślnej korzysta się również ze stosowanych w teorii ruchu ciepła bezwymiarowych liczb Fouriera Fo

0x01 graphic
(73)

i Biota Bi

0x01 graphic
(74)

Wprowadzono także bezwymiarową współrzędną odległościową

0x01 graphic
(75)

Liczba Kirpiczewa Ki we wzorze (69) jest funkcją liczb Fouriera Fo i Biota Bi. Można ją wyznaczyć na podstawie diagramu przedstawionego na rys. 5.1 lub za pomocą wzorów empirycznych (59 - 63)

K = K(Fo, Bi)

Podstawowy wzór metody wykreślnej prognozowania temperatury powietrza

Podstawowy wzór metody wykreślnej prognozowania dotyczy przyrostu temperatury powietrza

Δt = tw - td (76)

w modelu odcinka wyrobiska górniczego.

Wzór ten ma postać

0x01 graphic
(77)

W wyrażeniu (72) uczestniczą dwa bezwymiarowe współczynniki

0x01 graphic
(78)

0x01 graphic
(79)

będące funkcjami liczb Fouriera Fo i Biota Bi oraz bezwymiarowej współrzędnej odległościowej 0x01 graphic
. W GIG przy korzystaniu z danych dotyczących funkcji Kirpiczewa Ki i stosowaniu maszyny matematycznej wykonano tablicę funkcji (78) i (79), a na rys.5.4 i 5.5 przedstawiono ich wykresy

0x01 graphic
(80)

0x01 graphic
(81)

dla przypadku Bi = ∞. Z wykresów funkcji (rys. 5.4 i 5.5) można korzystać, gdy tylko Bi ≥ 25.

Dalsze symbole wzoru (72), tj. 0x01 graphic
oraz 0x01 graphic
, mają następującą treść fizyczną:

0x01 graphic
jest różnicą miedzy temperaturą 0x01 graphic
pierwotną skał, w płaszczyźnie wyznaczonej przez przekrój dopływu d, a temperaturą powietrza td w tym przekroju, czyli

0x01 graphic
(82)

0x01 graphic
natomiast jest nieskorygowanym przyrostem temperatury powietrza w odcinku bocznicy sieci kopalnianej wynikłym z działania dodatkowych źródeł Zi ciepła, wyznaczonym za pomocą wzoru

0x01 graphic
(83)

Aby uzyskać rzeczywisty przyrost temperatury powietrza w odcinku bocznicy jako skutek działania źródeł dodatkowych Zi, gdzie i = 1, 2 ... n, należy w myśl wzoru (72) pomnożyć nieskorygowany przyrost Δtz przez współczynnik korygujący Ω. Współczynnik Ω uwzględnia oddziaływanie górotworu na przyrost temperatury powietrza, pochodzący z dodatkowych źródeł Zi.

Efekty cieplne związane z nachyleniem odcinka wyrobiska (β ≠ 0) - a w tym efekty sprężania lub rozprężania się powietrza - jak również zmianę temperatury pierwotnej skał zgodnie z wzorem (68), uwzględniono w metodzie wykreślnej w postaci dodatkowego źródła ciepła Zi zawartego we wzorach (72) i (83). Źródło to ma zagęszczenie 0x01 graphic
strumienia powietrza równe

0x01 graphic
(84)

gdzie g - przyspieszenie ziemskie,

0x01 graphic
(85)

W literaturze rozpowszechniony jest pogląd na efekty cieplne związane ze zmianą wilgotności powietrza kopalnianego, według którego ciepło parowania wody pobierane jest z powietrza. Zgodnie z tym, w znanych prognostycznych metodach obliczeniowych uwzględniano te efekty jako dodatkowe (ujemne) źródło ciepła. Prowadziło to do dużych błędów w prognozie temperatury powietrza kopalnianego. Wobec tego w opracowanej wykreślnej metodzie oraz w podstawowym wzorze (72) nie uwzględniono dodatkowego (ujemnego) źródła ciepła związanego ze wzrostem wilgotności powietrza.

Przy stosowaniu wzoru (72) obliczania przyrostu temperatury powietrza Δt w odcinku bocznicy sprowadza się do wyznaczenia różnicy temperatur ΔΘd (wzór (82)), przyrostu Δtz(wzór (83)), liczby Fouriera Fo (wzór (73)) i bezwymiarowej współrzędnej0x01 graphic
(wzór (75)), odczytania wartości funkcji 0x01 graphic
i 0x01 graphic
z wykresów (rys. 5.4 i 5.5 ) oraz wykonania działania zgodnie ze wzorem (72). Tak obliczony przyrost Δt służy do wyznaczenia temperatury tw powietrza w przekroju wypływu w, w myśl przekształconego wzoru (76)

tw = td + Δt (76')

Jeśli przekrój wypływu (w) odcinka nie leży w węźle, do którego dopływa więcej niż jeden strumień powietrza, to temperaturę tw przyjmuję się równą temperaturze przekroju dopływu td odcinka następnego, w którym wyznacza się ponownie przyrost temperatury powietrza zgodnie z wzorem (72).

0x08 graphic

Rysunek 5.4 Wykres bezwymiarowego współczynnika ψ jako funkcji liczby Fouriera F0 dla różnych wartości parametru 0x01 graphic
wyznaczony ze wzoru (60) przy założeniu, że liczba Biota Bi .

0x08 graphic

Rysunek 5.5 Wykres bezwymiarowego współczynnika Ω jako funkcji liczby Fouriera F0 dla różnych wartości parametrów 0x01 graphic
, wyznaczony ze wzoru (61), przy założeniu, że liczba Biota Bi → ∞.

4. Metoda J. Wacławika

Istnieje wiele sposobów określania temperatury powietrza w kopalni, podawanych przez różnych autorów zarówno polskich jak i zagranicznych. W latach 1971 - 74 J. Wacławik podał rozwiązanie, które uwzględnia także wpływ sezonowych zmian temperatury powietrza na wlocie do kopalni.

Temperatura średnia

Dla uproszczenia przyjmuje się, że przekrój wyrobiska ma kształt kołowy o promieniu hydraulicznym

0x01 graphic
(86)

gdzie

A - powierzchnia przekroju wyrobiska, m2,

B - obwód przekroju wyrobiska, m.

Ponadto zakłada się, że pole temperatur skał jest symetryczne względem osi wyrobiska. Pomijając przewodnictwo w skałach w kierunku równoległym do wyrobiska, bilans ciepła w skałach ujmuje równanie.

0x01 graphic
(87)

gdzie

r - promień wodzący współrzędnych cylindrycznych r0 ≤ r < ∞, m,

t - czas przewietrzania wyrobiska, lata,

0x01 graphic
- temperatura górotworu, °C

a - współczynnik wyrównania temperatury skał (tabela 5.4)

0x01 graphic
(88)

przy czym

λs - współczynnik przewodnictwa cieplnego skał (tabela 5.5),

ρs - gęstość skał, kg/m3,

cs - ciepło właściwe górotworu, J/(kg K).

Tabela 5.4. Współczynnik wyrównania temperatury skał a, (m2/s) 106

Rodzaj skał

Wartość współczynnika a, (m2/s) 106

Górnośląskie Zagłębie Węglowe

Dolnośląskie

Zagłębie Węglowe

Lubelskie Zagłębie Węglowe

Lubińsko - Głogowskie Zagłębie Miedzi

Zlepieńce i żwirowce

Piaskowce gruboziarniste

Piaskowce drobnoziarniste

Łupki piaszczyste

Łupki ilaste

Węgle kamienne

Wapień

Kreda pisząca

Dolomit

Margle miedzionośne

1,65

1,5

1,3

0,9

0,8

0,32

0,9

1,45

1,26

1,0

0,76

0,47

1,6

1,2

1,02

0,35

1,25

0,95

1,1

1,05

1,2

1,2

Tabela 5.5 Współczynnik przewodnictwa ciepła skał λs, W/m °C

Rodzaj skał

Wartość współczynnika λs, W/m oC

Górnośląskie Zagłębie Węglowe

Dolnośląskie

Zagłębie Węglowe

Lubelskie Zagłębie Węglowe

Lubińsko - Głogowskie Zagłębie Miedzi

Zlepieńce i żwirowce

Piaskowce gruboziarniste

Piaskowce drobnoziarniste

Łupki piaszczyste

Łupki ilaste

Węgle kamienne

Wapień

Kreda pisząca

Dolomit

Margle miedzionośne

3,4

3,5

3,1

2,2

2,1

0,60

2,3

3,4

2,9

2,5

1,9

0,66

4,04

2,5

2,5

0,4

2,8

1,8

2,3

3,0

2,5

2,5

Jeżeli przyjmie się stały wydatek powietrza 0x01 graphic
i pominie zmiany energii kinetycznej, to równanie energii dla powietrza przyjmie postać

0x01 graphic
(89)

gdzie

ρ - gęstość powietrza, kg/m3,

qs - ilość ciepła dopływającego do powietrza od skał na jednostkowej długości

przewodu i w jednostce czasu, J/(kg K),

T - temperatura powietrza, °C,

g - przyspieszenie ziemskie, m/s2,

s - współrzędna przestrzenna, mierzona wzdłuż osi wyrobiska, m,

qd - dodatkowa ilość ciepła (np. od maszyn), przypadająca na jednostkę długości

wyrobiska i jednostkę czasu, J/(kg K)

z - współrzędna wysokościowa, m.

Jeżeli temperatura pierwotna skał na wlocie do wyrobiska wynosi 0x01 graphic
, to temperaturę pierwotną w innym punkcie oblicza się ze wzoru

0x01 graphic
(90)

gdzie:

0x01 graphic
- temperatura pierwotna skał na wlocie do wyrobiska, °C,

σ - gradient geotermiczny, °C/m,

0x01 graphic
- przyrost temperatury pierwotnej skał na jednostkę długości wyrobiska.

Dla równania (87) przyjmuje się warunek brzegowy w postaci

0x01 graphic
(91)

wyrażający, że ilość ciepła wymienionego pomiędzy skałami a powietrzem jest proporcjonalna do różnicy temperatury ścianek wyrobiska i przepływającego powietrza. Współczynniki przejmowania ciepła dla przewodów gładkich zestawiono w tabeli 5.6.

Carslaw i Jeager podali rozwiązanie równania różniczkowego (87) przy założeniu, że temperatura pierwotna skał wokół wyrobiska 0x01 graphic
jest stała i stała jest także temperatura powietrza na wlocie T(s0). Podali oni także wyrażenie na ilość ciepła odbieranego przez powietrze z jednostkowej długości przewodu i w jednostkowym czasie

0x01 graphic
(92)

gdzie

u - zmienna całkowania,

0x01 graphic
- liczba Biota,

0x01 graphic
- liczba Fouriera.

Wyrażenie w nawiasie {} będzie dalej oznaczone jako q*.

Tabela 5.6 Współczynnik α przejmowania ciepła dla przewodów gładkich w zależności od prędkości przepływu w i promienia hydraulicznego ro.

Promień hydrauliczny r0, m

Prędkość przepływu w, m/s

0,25

0,5

1,0

2,0

4,0

8,0

16,0

0,5

1,0

1,5

2,0

2,5

3,0

3,5

4,0

1,10

0,96

0,92

0,83

0,80

0,77

0,74

0,73

1,91

1,67

1,61

1,45

1,39

1,34

1,30

1,26

3,33

2,90

2,80

2,50

2,41

2,33

2,26

2,20

5,80

5,05

4,88

4,40

4,20

4,05

3,93

3,83

10,10

8,79

8,49

7,65

7,32

7,06

6,84

6,66

17,58

15,31

14,79

13,33

15,74

12,29

11,91

11,60

30,11

26,21

25,33

22,82

21,82

21,04

20,40

19,87

Tabela 5.7 Bezwymiarowy strumień cieplny q* dla Bi = jako funkcja liczby Fouriera Fo.

Fo

q*

Fo

q*

Fo

q*

Fo

q*

Fo

q*

Fo

q*

0,01

0,02

0,03

0,04

0,05

0,06

0,08

38,51

28,10

23,47

20,72

18,83

17,43

15,47

0,1

0,2

0,3

0,4

0,5

0,6

0,8

14,13

10,78

9,28

8,37

7,75

7,29

6,63

1

2

3

4

5

6

8

6,18

5,03

4,50

4,17

3,59

3,78

3,53

10

20

30

40

50

60

80

3,35

2,90

2,68

2,54

2,44

2,36

2,25

100

200

300

400

500

600

800

2,17

1,96

1,85

1,78

1,72

1,68

1,62

1000

2000 3000

4000

5000

6000

8000

1,58

1,45

1,39

1,35

1,32

1,29

1,26

Równanie energii (89) przyjmuje postać

0x01 graphic
(93)

gdzie qd oznacza ciepło doprowadzone od maszyn w jednostce czasu i na jednostkę długości wyrobiska.

Gdy dla s = 0 Ts = T0 śr, wówczas rozwiązaniem równania jest funkcja

0x01 graphic
(94)

Z analizy wzoru (94), przy pominięciu dopływu ciepła do maszyn (qd = 0), wynika, że gdy czas przewietrzania t jest mały (t → 0), wówczas temperatura powietrza zmierza do wartości temperatury pierwotnej skał. Wtedy bowiem dla nieskończenie dużej wartości współczynnika przejmowania ciepła q* jest

0x01 graphic
(95)

dla

t → 0, q*

Gdy czas przewietrzania rośnie, wówczas q* → 0. Wtedy ze wzoru (94) otrzymuje się

0x01 graphic
(96)

t, q* → 0

skąd wynika, że temperatura powietrza ze wzrostem okresu przewietrzania maleje do wielkości temperatury przemiany adiabatycznej (temperatura wzrasta w stosunku do wlotowej o wielkość 0x01 graphic
) na 1 m zmiany głębokości.

Średnia temperatura powietrza na podszybiach szybów wdechowych

Temperatura na wlocie do szybu wdechowego ulega wahaniom w zależności od pory dnia i pory roku. Ponieważ rozwiązanie (94) podane jest przy założeniu niezmiennej temperatury na wlocie do szybu, w wyniku tego do obliczenia średniej temperatury powietrza na podszybiu należy podstawić za T0 śr średnią roczną temperaturę powietrza wlotowego. Zależy ona od miejscowości i jest w przybliżeniu równa temperaturze skał na głębokości 25 do 35 m od powierzchni ziemi.

Tabela 5.8 Wartość średnia, amplituda i przesunięcie fazowe temperatury powietrza atmosferycznego

Miejscowość

Temperatura średnia

T0 śr, oC

Amplituda

A0, oC

Przesunięcie fazowe

ε0

Katowice

Bytom

Gliwice

Legnica

7,5

8,0

7,8

8,4

10,3

10,3

10,2

10,2

-1,89

-1,89

-1,89

-1,89

Wtedy 0x01 graphic
- T0 śr ≈ 0. Ponieważ dla szybu wdechowego 0x01 graphic
, więc wzór (76) na średnią temperaturę Tp na podszybiu można uprościć do postaci

0x01 graphic
(97)

gdzie s w tym przypadku równa się głębokości szybu H.

Gdy ilość przepływającego powietrza jest tak duża, że

0x01 graphic

wówczas temperatura na podszybiu wyniesie

0x01 graphic
(98)

Temperatura średnia powietrza w wyrobisku poziomym

Dla wyrobiska poziomego 0x01 graphic
, więc przy pominięciu ciepła dopływającego od maszyn (qd = 0), otrzymuje się na podstawie wzoru (63)

0x01 graphic
(99)

Temperatura powietrza w wyrobiskach poziomych zależy przede wszystkim od temperatury powietrza na wlocie, temperatury pierwotnej skał, długości wyrobiska i ilości przepływającego powietrza. Ponieważ jedyną wielkością, którą można zmieniać w szerszych granicach jest ilość powietrza 0x01 graphic
, więc najbardziej interesujący jest wpływ natężenia przepływu na temperaturę powietrza.

5. Metoda Szczerbania, Kremniewa, Żurawlenki

Przy prognozowaniu temperatury powietrza tą metodą zależność na prognozowaną temperaturę powietrza uzależniono od czasu przewietrzania wyrobiska.

Wyrobiska przewietrzane powyżej roku

  1. Poziome

Temperaturę powietrza w poziomym wyrobisku określa się ze wzoru

0x01 graphic

(100)

przy czym dla uproszczenia przyjmuje się:

0x01 graphic
(101)

0x01 graphic
(102)

gdzie

ϕ1, ϕ2 - odpowiednia wilgotność względna powietrza na początku (1) i na końcu (2) wyrobiska, %, przy czym:

ϕ2 = ϕ1 ± ξL

ξ - zmiana wilgotności względnej powietrza na drodze od szybów do przodków wybierkowych (do miejsca eksploatacji), odnoszący się do 1 m tej drogi,

L - długość wyrobiska, m,

t1 - temperatura powietrza na początku wyrobiska, °C

tp - temperatura skał na początku wyrobiska, °C,

ΣQM - ilość ciepła wydzielana w miejscowych źródłach w wyrobisku (na danym odcinku wyrobiska), kcal/h,

G - strumień masy powietrza przepływającego przez wyrobisko, kg/h,

qo - ciepło utleniania z 1 m2 powierzchni wyrobiska w jednostce czasu, kcal/m2 h,

kτ - współczynnik niestacjonarnej wymiany ciepła

m, n, l - współczynniki uwzględniające zależność zawartości ciepła od temperatury równe:

m = 3,872, n = 0,169, l = 0,0189

cp - pojemność cieplna powietrza przy stałym ciśnieniu (można przyjąć 0,25 kcal/kg °C)

r - ciepło parowania wody (można przyjąć 0.59 kcal/g),

U - obwód wyrobiska, m.

Jeśli wilgotność względna strumienia powietrza od początku do końca wyrobiska zwiększa się to wartość ξ we wzorze (102) przyjmuje znak dodatni, a jeśli się zmniejsza to znak ujemny.

Jeśli wiemy jaka jest oczekiwana temperatura na wylocie z wyrobiska, możemy obliczyć jaka powinna być temperatura wlotowa na podstawie następującego wzoru:

0x01 graphic
(103)

przy czym dla uproszczenia przyjmuje się:

0x01 graphic
(104)

0x01 graphic
(105)

0x01 graphic
(106)

gdzie t2 oznacza temperaturę na końcu wyrobiska.

Przybliżone obliczenia dla wyrobisk poziomych można przeprowadzić wg wzoru zestawionego dla niezmiennej wilgotności względnej. Wzór ten ma postać

0x01 graphic
(107)

Wartości K i D przy różnych wilgotnościach względnych przedstawione są w tabeli 5.9.

Tabela 5.9 Wartości K i D przy różnych wilgotnościach względnych ϕ.

Wilgotność względna ϕ, %

K

D

75

80

85

90

95

-19,1

-18,3

-17,4

-16,7

-16,1

29,3

27,5

25,9

24,4

23,2

Wartość B zależną od temperatury i wilgotności względnej można odczytać z tablic podanych w [3].

  1. wyrobiska nachylone

Temperaturę powietrza w nachylonym wyrobisku oblicza się ze wzoru

0x01 graphic
(108)

przy czym dla uproszczenia przyjęto:

0x01 graphic
(109)

gdzie

ψ - kąt nachylenia wyrobiska, °

σ - gradient geotermiczny

0x01 graphic
, °C/m

Uproszczona formuła dla wyrobisk nachylonych ma postać

0x01 graphic
(110)

Wzór na temperaturę powietrza płynącego w kierunku odwrotnym do normalnego przepływu w nachylonych wyrobiskach ma postać:

0x01 graphic
(111)

przy czym dla uproszczenia przyjęto:

0x01 graphic
(112)

Wyrobiska przewietrzane do roku

  1. poziome

Temperatura na końcu wyrobiska lub rejonu (oddziału) określamy ze wzoru

0x01 graphic
(113)

gdzie

φ - współczynnik wilgotności (przy zmiennej wilgotności względnej ϕ)

0x01 graphic
; 0x01 graphic
(114)

Przy niezmiennej wilgotności względnej

0x01 graphic
(115)

Γ - współczynnik wymiany ciepła, równy przy występowaniu grzejących i chłodzących rurociągów

0x01 graphic
(116)

natomiast przy braku grzejących i chłodzących rurociągów:

0x01 graphic
(117)

T - współczynnik temperaturowy równy w przypadku występowania grzejących lub chłodzących rurociągów

0x01 graphic
(118)

KT, Kx - współczynnik przekazywania ciepła od grzejącego (chłodzącego) rurociągu,

kcal/(m2 h °C),

Ut, UX - obwód przekroju grzejącego (chłodzącego) rurociągu, m.

przy ich braku (KT = 0, KX = 0)

0x01 graphic
(119)

δ - współczynnik wilgotności określony według wzoru

0x01 graphic
(120)

gdzie B' oznacza ciśnienie barometryczne mierzone w mm Hg, natomiast wartości n' i p'cp są zestawione w tabeli 5.10 dla różnych przedziałów temperatur.

Tabela 5.10 Wartości n' i p'cp dla różnych temperatur

t, °C

n'

p'cp

t, °C

n'

p'cp

0-10

2-12

4-14

6-16

8-18

10-20

12-22

14-24

705

795

885

990

1105

1240

1410

1575

5,7

6,5

7,5

8,5

9,5

11,0

12,5

14,5

16-26

18-28

20-30

22-32

24-34

26-36

28-38

30-40

1725

1915

2120

2370

2630

2880

3200

3480

16,3

18,5

20,5

23,0

26,0

29,0

32,0

36,0

Dla ułatwienia obliczeń wg wzoru (113) wykonano nomogram pokazany na rys. 5.6 Bezwymiarową temperaturę powietrza 0x01 graphic
w wyrobisku na tym monogramie określa się jako funkcje bezwymiarowych współczynników Γ, φ, B1. Korzystanie z nomogramu pokazane jest za pomocą strzałek. Na jego podstawie określa się bezwymiarową temperaturę 0x01 graphic
a potem temperaturę na końcu wyrobiska, przy czym

0x01 graphic
; 0x01 graphic
(121)

0x01 graphic
(122)

Nomogram jest wykonany dla przypadku gdy nie ma grzejących i chłodzących rurociągów

  1. wyrobiska nachylone

Temperaturę na końcu nachylonego wyrobiska określa się wg wzoru

0x01 graphic
(123)

gdzie

0x01 graphic
(124)

Wszystkie pozostałe oznaczenia jak w poprzednich wzorach dla poziomych wyrobisk. Przy określaniu T, tp okazuje się temperaturą w górnej części nachylonego wyrobiska.

0x08 graphic
Przy obliczaniu różnicy wysokości ΔH = (H2 - H1) większej od 100 m w ΣQM należy uwzględnić ilość ciepła potrzebną na ochłodzenie powietrza przy jego rozszerzeniu, równym ΔQpac ≅ 0,0025GΔH, kcal/h, gdzie G - strumień masy powietrza, kg/h.

Rysunek 5.6 Nomogram do obliczeń temperatury w wyrobiskach górniczych.

6. Metoda PTO-2

Temperaturę powietrza kopalnianego w danym wyrobisku z wentylacją opływową oblicza się z wzoru:

tw = td + ψΔΘd + ΩΔtz + ΦΔtx

td - temperatura powietrza w przekroju dopływu wyrobiska , °C,

ψ - bezwymiarowy współczynnik określony wzorem

ψ = 1- exp[- K(Fo,Bi) x

Fo- liczba Fouriera

0x01 graphic

Bi- liczba Biota

Bi = 0x01 graphic

X - bezwymiarowa współrzędna odległościowa:

X = 0x01 graphic

Gdzie:

K - K(Fo,Bi) - liczba Kirpiczewa wyznaczana z wykresu (rys.1)

a - współczynnik wyrównywania temperatury skał wyznaczany ze wzoru:

a = 0,35(3-2kw)10-6

kw - oznacza jaką część całkowitego obwodu wyrobiska B stanowi węgiel

τ - czas przewietrzania wyrobiska [s]

rs - promień równoważny wyrobiska:

rs = 0,1592B

λ - współczynnik przewodzenia ciepła skał:

λ = 0,395(6-5kw)

l - długość wyrobiska, m,

m - strumień masy powietrza w wyrobisku [kg/s]

cp - właściwa pojemność cieplna powietrza przy stałym ciśnieniu [J/(kg K)]

α - współczynnik wnikania ciepła z pobocznicy wyrobiska [W/(m2 K)]

α = 0x01 graphic

ΔΘd - różnica między temperaturą pierwotną skał tpg , a temperaturą powietrza w płaszczyźnie przekroju dopływu wyrobiska:

ΔΘd = tpg - td

Dalsze symbole oznaczają:

Ω - bezwymiarowy współczynnik :

Ω = 0x01 graphic

Δtz - przyrost temperatury powietrza w wyrobisku , wynikający z działania w nim dodatkowych źródeł ciepła:

Δtz = Σ 0x01 graphic

Qi - ciepło wydzielone w jednostce czasu z i-tego źródła ciepła [W]

Ciepło wydzielone z dodatkowych źródeł

Q1 - ciepło wydzielone z maszyn lub urządzeń elektrycznych czynnych w wyrobisku wyznacza się ze wzoru:

Q1 = 0,1*N

N - moc znamionowa silnika maszyny lub urządzenia [W]

Q2 - ciepło wydzielone z procesu utleniania węgla :

Q2 = qo * kw * B * l

qo - 4,36 [W/m2] dla wyrobisk korytarzowych przewietrzanych powyżej 3 lat

qo - 7,56 [W/m2] dla pozostałych wyrobisk korytarzowych

qo - 16,28 [W/m2] dla wyrobisk eksploatacyjnych

kw - oznacza jaką część całkowitego obwodu wyrobiska B stanowi węgiel

l - długość wyrobiska [m.]

Q3 - związane jest z wymianą ciepła między rurociągiem, a powietrzem w wyrobisku wyznacza się :

Qi=r = 7,719*(wd*Dr)0,62*(tr-td)* lr

wd - prędkość powietrza w przekroju /d/ wyrobiska [m/s]

Dr - średnica rurociągu [m.]

tr - temperatura powierzchni rurociągu [°C]

td - temperatura powietrza w przekroju dopływu wyrobiska [°C]

lr - długość rurociągu w wyrobisku

Q4 - ciepło związane z wymianą ciepła między transportowanym na taśmociągu urobkiem a powietrzem w wyrobisku wyznacza się:

Q4 = 1,1*mw*L0,8*[tpg-(tn+b)]

mw - masa urobku przenoszona w jednostce czasu taśmociągiem w wyrobisku [kg/s]

L - długość taśmociągu w wyrobisku [m.]

tpg - temperatura węgla w caliźnie przodka wydobywczego z którego węgiel jest przenoszony taśmociągiem [°C]

tn - szacowana średnia temperatura na termometrze wilgotnym w części wyrobiska , gdzie znajduje się taśmociąg [°C]

b - empiryczna wielkość, która zależy od temperatury węgla w caliźnie w przodku wydobywczym oraz od masy urobku transportowanego taśmociągiem.

Dla tp 35 °C:

b = 7 °C, dla 11 kg/s mw 33 kg/s

b = 4 °C, dla 33 kg/s mw 55 kg/s

b = 2,5 °C, dla 55 kg/s mw 84 kg/s

Dla tp = 30 °C ± 4° C :

b = 5 °C, dla 11 kg/s mw 33 kg/s

b = 2 °C, dla 33 kg/s mw 55 kg/s

b = 0,5 °C, dla 55 kg/s mw 84 kg/s

Φ - bezwymiarowy współczynnik:

Φ = ε'*Ω*(1- 0x01 graphic
)

ε' - liczba bezwymiarowa , która jest wielkością empiryczną zależną od czasu przewietrzania τ:

0,8 - dla bardzo mokrych szybów wdechowych niezależnie od czasu przewietrzania

0,5 - dla normalnie mokrych szybów wdechowych niezależnie od czasu przewietrzania

0,5 - dla wyrobisk poziomych i nachylonych o czasie przewietrzania do 3 m-cy

0,32 - dla wyrobisk poziomych i nachylonych o czasie przewietrzania do 6 m-cy

0,16 - dla wyrobisk poziomych i nachylonych o czasie przewietrzania powyżej 6 m- cy

Δtx - przyrost temperatury powietrza w wyrobisku związany ze zmianą stanu skupienia wody przez dany wzór :

Δtx = - 0x01 graphic
* Δx

a' - ciepło parowania wody [J/kg]

Δx - przyrost stopnia zawilżenia powietrza w wyrobisku

Stopień zawilżenia powietrza X w danym wyrobisku

wyznacza się z wzoru:

xw = xd + Δx

xw - stopień zawilżenia powietrza w przekroju wypływu wyrobiska [kg/kg]

xd - stopień zawilżenia powietrza w przekroju dopływu wyrobiska [kg/kg]

Δx - przyrost stopnia zawilżenia w wyrobisku [kg/kg]

Δx = l * χ

l - długość wyrobiska [m.]

χ - jednostkowy stopień zawilżenia powietrza [kg/(kg m.)]

Jednostkowy stopień zawilżenia powietrza wyznacza się z bezpośrednich pomiarów , jakie wykonuje się w głębokich kopalniach. Na podstawie tych pomiarów dla każdego rodzaju wyrobiska w zależności od przeznaczenia można wyznaczyć jednostkowy przyrost stopnia zawilżenia korzystając ze wzoru:

χn = 0x01 graphic

Przy prognozowaniu stopnia zawilżenia w wyrobisku należy do obliczeń wykonywanych zgodnie ze wzorem na Δx wprowadzić korektę jednostkowego stopnia zawilżenia stosując wzór:

χ = χm = 0x01 graphic

gdzie:

m - strumień masy powietrza w projektowanym wyrobisku [kg/s]

mm - strumień masy powietrza w istniejącym wyrobisku dla którego wyznaczono χn

φd - szacowana wilgotność względna powietrza w projektowanym wyrobisku

φdm - wilgotność względna powietrza w istniejącym wyrobisku dla którego wyznaczono χm.

Θ = Θ(φ) - eksperymentalny potencjał ruchu wilgoci wyznaczony z tablicy:

0,50

0,51

0,53

0,55

0,57

0,59

0,60

0,61

0,63

0,65

0,67

0,69

0,70

0,71

26,2

26,6

27,5

28,4

29,4

30,3

30,9

31,4

32,5

33,6

34,7

36,1

36,8

37,4

0,73

0,75

0,77

0,80

0,83

0,85

0,87

0,90

0,91

0,93

0,95

0,97

0,99

1,00

38,9

40,6

42,6

45,9

49,8

53,2

56,9

63,7

66,1

71,5

78,4

86,8

95,5

100

Znając zaprognozowaną temperaturę powietrza tw (mierzoną termometrem suchym ) i stopień zawilżenia xw wyznaczamy z wykresem i-x temperaturę powietrza t'w (mierzonym termometrem mokrym). Następnie korzystając z niżej podanych wzorów wyznaczamy katastopnie wilgotne:

W 1 K' = (0,35+0,850x01 graphic
)*(36,5-t')

W> 1 K' = (0,10+1,10x01 graphic
)*(36,5-t')

w - prędkość powietrza w wyrobisku [m/s]

t' - temperatura powietrza na termometrze mokrym [°C]

Znając temperaturę powietrza i katastopnie wilgotne można dać odpowiedź czy podjęte środki profilaktyczne są wystarczające dla zapewnienia prawidłowych warunków klimatycznych w danym wyrobisku.

?

,może N, Nz lub podobnie



Wyszukiwarka

Podobne podstrony:
Klimatyzacja - praca, 5. Analiza znanych z literatury metod, ANALIZA ZNANYCH Z LITERATURY METOD PROG
Klimatyzacja - praca, referat PTO-2, Politechnika Wroclawska
Praca domowa 2a Analiza Matematyczna
Motyw Bohatera, Notatki, Filologia polska i specjalizacja nauczycielska, Analiza dzieła literackiego
Parafraza, Polonistyka, Poetyka i analiza dzieła literackiego
S. Balbus, poetyka z el. teorii literatury, analiza dzieła literackiego
Stylizacja, filmoznawstwo, Semestr 1, Poetyka z Analizą Dzieła Literackiego
WIERSZ WOLNY, filmoznawstwo, Semestr 1, Poetyka z Analizą Dzieła Literackiego
Chwilowe temperatury powietrza zewnetrznego w okresie letnim, Pomoce naukowe, Wentylacja i klimatyza
praca magisterska licencjacka Analiza wybranych możliwości finansowania projektów innowacyjnychx
metafora, filmoznawstwo, Semestr 1, Poetyka z Analizą Dzieła Literackiego
34a. Metafora i metonimia. Fabuła i akcja w dramacie. Kasia Kabat, poetyka z el. teorii literatury,
praca, Motyw pracy, MOTYW PRACY W LITERATURZE,SPOJRZENIE NA PRACĘ
Praca domowa 1a Analiza Matematyczna
sylabotonizm, filmoznawstwo, Semestr 1, Poetyka z Analizą Dzieła Literackiego
odpowiedzi na ćwiczenia z 10.10, filmoznawstwo, Semestr 1, Poetyka z Analizą Dzieła Literackiego

więcej podobnych podstron