Extrém alacsony frekvenciájú mágneses terek mérése környezetünkben TDK dolgozat
Zsigmond Anna Julia ELTE TTK Fizika BSc 2. évfolyam
Témavezető: dr. Veres Gábor ELTE TTK Atomfizikai Tanszék
ELTE TTK 2008
Tartalomjegyzék Bevezetés
2
1. Mérőműszerek 1.1. A kis tekercs induktivitásának kiszámolása 1.2. A kis tekercs induktivitásának kimérése . . 1.3. A tekercsek érzékenysége . . . . . . . . . . 1.4. Térerősségmérő műszer . . . . . . . . . . . 1.5. A műszer kalibrálása . . . . . . . . . . . . 1.6. Mérési tartományok és módszerek . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
4 4 6 7 8 9 10
2. Helmholtz-tekercsek, biológiai vizsgálatok segítése 12 2.1. Helmholtz-tekercsek terének kiszámítása . . . . . . . . . . . . . . . . 12 2.2. Helmholtz-tekercsek terének mérése . . . . . . . . . . . . . . . . . . . 14 3. Elektromágneses terek környezetünkben 3.1. Katódsugárcsöves monitor . . . . . . . . 3.2. Erősáramú földvezeték . . . . . . . . . . 3.3. Nagyfeszültségű távvezeték . . . . . . . . 3.4. Közlekedési eszközök . . . . . . . . . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
17 17 17 21 24
4. Összefoglalás
27
Irodalomjegyzék
28
1
Bevezetés Az elektromosság felfedezése óta egyre több villamos berendezés kerül környezetünkbe, melyek nem-ionizáló sugárzásokat és kvázisztatikus elektromos és mágneses tereket keltenek maguk körül. A nem-ionizáló sugárzások hullámhossza 100 nm-nél nagyobb, frekvenciája 3 PHz-nél kisebb, ezeket foglalja össze az 1. táblázat. Mi ezek közül az LF, VLF, ELF frekvenciatartományokkal foglalkozunk. Sugárzás típusa Ionizáló Ultraibolya Látható fény Infravörös Extrém magas frekvencia (EHF) Szuper-magas frekvencia (SHF) Ultra-magas frekvencia (UHF) Nagyon magas frekvencia (VHF) Magas frekvencia (HF) Közép frekvencia (MF) Alacsony frekvencia (LF) Nagyon alacsony frekvencia (VLF) Extrém alacsony frekvencia (ELF) Sztatikus terek
Frekvencia-tartomány >3 PHz 3-0,75 PHz 750-350 THz 350-0,3 THz 300-30 GHz 30-3 GHz 3-0,3 GHz 300-30 MHz 30-3 MHz 3-0,3 MHz 300-30 kHz 30-0,3 kHz 16-300 Hz 0 Hz
Hullámhossz <100 nm 100-400 nm 400-800 nm 0,8-1000 m 1-10 mm 1-10 cm 10-100 cm 1-10 m 10-100 m 100-1000 m 1-10 km 10-1000 km >1000 km Végtelen
1. táblázat. Elektromágneses sugárzások frekvencia- és hullámhossztartományai A nem-ionizáló elektromágneses sugárzások biológiai hatásai és mérése az utóbbi évtizedekben a közérdeklődés középpontjába kerültek. Az egészségügyi határértékeket az ICNIRP [1] illetve az EU ajánlásában lehet megtalálni, ahol a lakosságra és a munkahelyre határoznak meg maximális értékeket. Ezek az értékek lakosságra (24 órára) 100 µT, munkahelyre (8 órára) 500 µT. Az ekkoránál kisebb terek az eddigi tapasztalatok alapján nem veszélyesek az egészségre. A különböző frekvenciáknak különböző a hatásuk. A nagyobb frekvenciák esetében a sugárzás hőhatását kell figyelembe venni (pl. mobiltelefon), alacsony frekvenciák esetében (50 Hz-es elektromos hálózat) pedig a testben indukálódott áramsűrűségek hatását. Az orvosi gyakorlatban különböző frekvenciájú változó tereket és sztatikus tereket is használnak. Mai tudásunk szerint az orvosi MRI vizsgálat során alkalmazott nagy sztatikus mágneses térnek (0,2-2 Tesla) és a rádiófrekvenciás sugárzásnak nincs káros biológiai hatása. Mi most főleg a hálózati frekvenciájú terekkel foglalkozunk [2]. Ezek biológiai hatásának vizsgálati eredményei között sok az ellentmondás. Vannak, akikre jótékony hatással van a mágneses térrel való kezelés, és vannak, akik meg vannak győződve róla, hogy az alvóhelyükön lévő elektromos berendezéseknek káros hatásuk van szervezetükre. Bizonyított tény azonban, hogy ha az indukálódott áramsűrűség meghalad egy bizonyos határértéket, az indukált áram képes közvetlenül ingerelni az ideg- és izomszövetet. 2-20 mA/m2 áramsűrűségeknek hatásuk van a sejtmemb2
rán működésére, 10-100 mA/m2 áramsűrűség-tartományban szövetszintű hatások figyelhetők meg. Minél nagyobb az indukált áramsűrűség, annál veszélyesebb lehet a hatás az idegrendszerre illetve a szívműködésre. A vizsgálatok során többféle mérési módszer van. Egyrészt lehet mérni a terek pillanatnyi nagyságát, másrészt pedig lehet mérni az időbeli lefolyást, amely sokkal összetettebb, mivel az expozíció nagysága és időbeni eloszlása nagyon változatos lehet. A következőkben elsősorban a mágneses terek nagyságának mérését tárgyaljuk. Az extrém alacsony frekvenciájú terek mérésére különböző eszközöket lehet vásárolni, melyek sokszor jelentős beruházást jelentenek a lakosság számára, pedig házilag is elkészíthető eszközzel is megbízhatóan és megfelelő pontossággal mérhetünk változó mágneses teret bizonyos frekvenciatartományokban. Mi egyszerű tekercseket - amik tönkrement elektromos eszközökből kivehetők, illetve mi is tekercselhetünk - használunk változó mágneses terek mérésére. A következő fejezetben részletezzük, hogyan mérhetjük a mágneses teret a tekercseinkkel. Ezután szerepel, hogyan segítettük az egyetemünkön akvítan művelt neurobiológiai kutatást azzal, hogy jellemeztük a patkányokat érő mágneses teret, amivel a saját tekercseinkről is pontosabbak lettek az adataink. Az utolsó fejezetben pedig bemutatjuk néhány a környezetünkben előforduló mágneses tér mérését földvezeték, távvezeték, közlekedési eszközök esetében, illetve néhány további érdekességet.
3
1.
Mérőműszerek
A változó mágneses teret tekercsekkel mérjük az indukció törvénye alapján. Ehhez ismernünk kell a tekercsek geometriai adatait, és a menetszámukat. A geometriai adatokat le lehet mérni, de ez nem mindig egyszerű, ha nem szabályos a tekercs alakja. Két tekercsünk van, amiket az 1. ábrán láthatunk együtt: a karika (2. ábra) minden adatát tudjuk, mivel mi készítettük: a menetszáma 200, az átlagos átmérője (27, 6±0, 5) cm. A kis tekercs (3. ábra) geometriai adatait lemértük, de a menetszámot nem tudjuk, mert egy szétszedett transzformátorból van. A felületének oldalai (2, 20 ± 0, 05) cm és (1, 55 ± 0, 05) cm, a hossza pedig (2, 00 ± 0, 05) cm.
1. ábra. Tekercseink: a "karika" és a "kis tekercs"
2. ábra. A karika adatai
1.1.
A kis tekercs induktivitásának kiszámolása
A kis tekercs menetszámának kiszámításához a tekercs induktivitását határozzuk meg. Párhuzamos rezgőkört alakítunk ki (4. ábra) egy kapacitással, és oszcilloszkóppal vizsgáljuk az áramkör válaszát a jelgenerátor különböző frekvenciájú szinuszos jelére. A mért feszültség amplitúdója ekkor U0 U0mert = q 1 + R R+2ReZ |Z|2 4
(1)
3. ábra. A kis tekercs adatai
4. ábra. Párhuzamos rezgőkör
lesz, ahol Z a párhuzamos LC kör impedanciája, amit a következő módon számolhatunk ki: 1 ZC = iωC ZL = iωL 1 1 1 1 1 r ωL = + + = + 2 + i(ωC − 2 ) 2 2 Z ZC ZL + r RC RC r +ω L r + ω 2 L2 1 + R2rC − 2ω 2LC 1 1 2 2 = 2 +ω C + |Z|2 RC r 2 + ω 2 L2
(2)
ahol RC a kondenzátor ellenállása és r a tekercs belső ellenállása. Az (1) és a (2) egyenletből kiszámolható, hogy a maximális feszültséget az alábbi körfrekvenciánál mérhetjük: vq u u 1 + 2r + 2r2 C + 2r − r2 C t RC L R L ω= (3) LC Ez azonban jó közelítéssel átalakítható a rezonanciafrekvenciánál jól ismert képletté (4), mert RC végtelennek tekinthető (kívül esik a multiméter méréshatárán), és r
5
általában nagyon kicsi. A mi tekercseink ohmikus ellenállása: a karikáé r◦ = 19, 8 Ω, a kis tekercsé r = 323, 4 Ω. 1 (4) ν= √ 2π LC A közelítés mértékét behelyettesítéssel ellenőriztük, ahol RC -t végtelennek tekintettük, mivel kívül esett a multiméter méréshatárán, és a kis tekercs ellenállásával számoltunk. Ekkor azt kaptuk, hogy a számolt frekvenciák között kevesebb, mint 0,2% a különbség, ami jóval kisebb a méréseinknél adódó egyéb hibáknál, tehát számolhatunk a (4) képlettel.
1.2.
A kis tekercs induktivitásának kimérése
A rezgőkör jelgenerátorra adott válaszából kimértük az amplitúdó frekvenciafüggését, amire a (1) és (2) képletek alapján illesztettük a rezonanciagörbét. Ez látható a 5. ábrán. 110 mert adatok illesztett gorbe
Amplitúdó (mV)
100 90 80 70 60 50 40 30 20 10 0 100
1000
10000
Frekvencia (Hz) 5. ábra. Rezonanciagörbe a kis tekercsből és a 94,6 nF-os kapacitásból összeállított rezgőkörnél A rezonanciagörbe maximuma szimmetrikusnak tekinthető néhány százalékos hibával, ezért rezonanciafrekvenciát olymódon keressük meg, hogy a maximális amplitúdó kis környezetében mindkét oldalán megkeressük ugyanazt az amplitúdót, és a két frekvencia átlaga a rezonanciafrekvencia. A kapacitást multiméterrel mérjük meg, és a rezgőkör rezonanciafrekvenciára vonatkozó képlete (4) alapján kiszámoljuk a tekercs induktivitását. Ez a kis tekercsre a 2. táblázatban lévő értékeket adta, egy nagyobb és egy kisebb kondenzátorral összekötve illetve egy kábellel, mint kapacitással összekötve. A kábel egy körülbelül 84 cm hosszú koaxális BNC kábel. Ennél a mérésnél észrevettük, hogy a kapacitásokhoz hozzá kell adnunk a kábel kapacitását is, így az egyes értékek az induktivitásra közelebb kerültek egymáshoz, de így is nagy a különbség 6
közöttük. Ennek oka valószínűleg, hogy a kapacitások mérése nem volt tökéletes a kis kondenzátornál illetve a kábelnél. kábel nagy kondenzátor kis kondenzátor
C(nF) ν (kHz) L(mH) 0,38±0,05 39,055 43,7±5,7 94,6±0,1 2,3015 50,5±0,1 1,12±0,1 19,100 61,0±5,4
2. táblázat. A kis tekercs induktivitása a rezonanciafrekvencia és a kapacitás alapján A menetszámot a hosszú szolenoidra vonatkozó képletből számoljuk ki: s lL N= ≈ 1553 ± 63 µ0 A
(5)
Ez a szám csak becslés, mivel a tekercsünk csak nagy közelítéssel tekinthető hosszú szolenoidnak. Ehelyett a képlet helyett használhatunk egy pontosabb képletet a többrétegű, kör keresztmetszetű tekercsre [3]: L=
0, 8r 2 N 2 6r + 9l + 10d
(6)
ahol L a tekercs induktivitása µH-ben, N a tekercs menetszáma, r a tekercs átlagos sugara, d a tekercs külső és belső sugarának különbsége, l a tekercs hossza (mindegyik inch-ben). Ebben a képletben a kis tekercs téglalap alakú felületét körrel közelítjük, tehát ez a képlet is csak közelítés. A képletbe való behelyettesítés után a tekercs menetszámára N ≈ 2230 ± 100 adódik. Itt látszik, hogy a menetszám és az induktivitás kapcsolata az egyes tekercseknél bonyolult kérdés. Az ω körfrekvenciájú, B0 amplitúdójú szinuszosan változó mágneses tér az A felületű tekercsben U = −N
∂(BA) ∂ = −NAB0 sin(ωt) = −NAB0 ω cos(ωt) = −U0 cos(ωt) ∂t ∂t
(7)
változó feszültséget hoz létre. Így a tekercsben indukált feszültség mérésével meghatározható a mágneses tér amplitúdója: B0 =
U0 NA2πν
Az effektív érték szinuszos jel esetén ennek mindig időbeli átlaga ennek π2 ≈ 0, 637-szerese.
1.3.
(8) √1 -szerese, 2
az abszolut érték
A tekercsek érzékenysége
A cél minden tekercsnél megállapítani a váltószámot az indukált feszültség és a tekercs síkjára merőleges mágneses indukció amplitúdója között, ami a tekercs érzékenysége. Ez minden frekvenciára más, de mi az 50 Hz-es hálózati frekvenciával foglalkozunk a legtöbbet, ezért ezeket az értékeket részletezzük. 7
A kis tekercsen hálózati frekvencia esetén 1 mV mért feszültség (6,01±0,48) µT nagyságú mágneses teret jelent az (5) képlet alapján számolva a menetszámot. Ha a (6) képlet alapján számoljuk a menetszámot, akkor 1 mV feszültség (4,19±0,21) µT nagyságú mágneses térnek feleltethető meg. Ezzel a tekerccsel 30-40 kHz-nél nagyobb frekvenciájú változó teret nem tudunk mérni, mivel ez közel esik a rezonanciafrekvenciájához, amit a 2. táblázatban is láthatunk, az első sorban, ahol a kapacitást csak a kábel kapacitása adta. Ennek következményeiről még lesz szó a 3. fejezetben. A karikára pontosan meg tudjuk mondani az egyes frekvenciákra ezt a váltószámot, mivel ismerjük a menetszámát. 50 Hz-es hálózati frekvencia esetén a karikán mért 1 mV feszültség 0,266±0,009 µT nagyságú változó mágneses teret jelent. Ezt az értéket ezentúl pontosnak tekintjük, és ezzel ellenőrizzük majd ki a kis tekercs pontosabb érzékenységét 50 Hz esetén az (5) és (6) képlet közelítő jellege miatt.
1.4.
Térerősségmérő műszer
A tekercseink mellé vásároltunk egy érzékenyebb műszert: egy ME 3030B elektromos térerősséget és mágneses indukciót mérő műszert (6. ábra). Ez elektromos és mágneses váltakozó terek mérésére alkalmas 2000 V/m illetve 2000 nT alatt, 16 Hz és 2 kHz közötti frekvenciatartományban. Mi mágneses teret mérünk vele, főleg a külső helyszínen folytatott mérésekkor.
6. ábra. ME 3030B elektromos és mágneses tér mérő Méréssel ellenőriztük, hogy a frekvencitartományban milyen a műszer karakterisztikája. Ismert változó mágneses térbe helyezve a frekvenciát változtatva vizsgáltuk, hogy mekkora értéket ír ki az ismert amplitúdóhoz képest. A mérés eredményét a 7. ábrán rögzítettük, amin látható, hogy azokon a frekvenciákon, ahol mi használjuk megbízhatóan tudunk vele mérni. A műszer egy komponenst mér, és érzékeny a mozgatásra csakúgy, mint a tekercsek, például ha forgatjuk, akkor a földmágnesség 30-60 µT nagyságú terét is érezni, ami a méréshatár sokszorosa. Tehát mérések közben figyelni kell, hogy nyugodtan tartsuk a műszert, és figyelembe vegyük az irányokat. Általában a mérésekkor annak az iránynak értékeit írjuk fel, ahol a maximális az érték, és csak ha ettől eltérünk, akkor beszélünk az irányokról. 8
1.1
1
Áteresztés
0.9
0.8
0.7
0.6
0.5 10
100
Frekvencia (Hz)
1000
7. ábra. A műszer karakterisztikája: a különböző frekvenciáknál a relatív áteresztése
1.5.
A műszer kalibrálása
A karikánkat, mivel ismerjük a geometriáját, használhatjuk ismert mágneses tér létrehozására is. Így egy ellenálláson keresztül megtápláljuk a karikát ismert amplitúdójú váltóárammal, amiből kiszámolhatjuk a karika tengelyén a mágneses tér amplitúdóját: NI0 µ0 R2 B0 = (9) 2 (R2 + z 2 ) 32 A műszert ráhelyezzük a karika forgástengelyére ismert z magasságban, és figyeljük a kiírt értéket, miközben a tápfeszültség amplitúdóját változtatjuk. A számolt és a kiírt értékpárokat ábrázolhatjuk, és az egyenes meredekségéből lehet tudni, hogy mennyivel kell szorozni a kiírt értékeket, hogy megkapjuk a mágneses indukció amplitúdóját. Az egyenesek meredeksége nem pontosan 1. Ez azért van, mert a műszer az időben váltakozó mágneses teret időben átlagolja. Ebben esetben szinuszos jelet mérünk, a szinuszfüggvény abszolutértékének átlaga az amplitúdó π2 ≈ 0, 637szerese. A szinuszos jel effektív értéke pedig az amplitúdó √12 -szerese. A különböző mérések során erre az időátlagolásra is figyelni kell. A 8. ábrán látható az 50 Hz-es frekvenciájú áram esete. Az egyenes meredeksége: 0, 666 ± 0, 001, tehát ha tiszta szinuszos 50 Hz-es frekvenciájú teret mérünk, akkor a műszer által kiírt értéket szorozni kell körülbelül 1, 5-tel. A 9. ábrán látható 200 Hz-es frekvencia mellett a kiírt és a számolt mágneses tér viszonya. Az egyenes meredeksége ebben az esetben 0, 6532 ± 0, 0005, tehát ekkor a szorzófaktor 1, 53, de ritkán mérünk tiszta szinuszos 200 Hz-es jeleket.
9
2000 1800 1600
Kiírt B (nT)
1400 1200 1000 800 600 400 200 0 0
500
1000
1500
2000
2500
3000
Számolt B0 (nT) 8. ábra. Mágneses indukció ME 3030B műszeren kiírt értéke a számolt amplitúdó függvényében 50 Hz-es frekvencia esetén
2000 1800
Kiírt B (nT)
1600 1400 1200 1000 800 600 400 200 0 0
500
1000
1500
2000
Számolt B0 (nT)
2500
3000
3500
9. ábra. Mágneses indukció ME 3030B műszeren kiírt értéke a számolt amplitúdó függvényében 200 Hz-es frekvencia esetén
1.6.
Mérési tartományok és módszerek
Ezekkel a mérőműszerekkel a mérhető amplitúdó- és frekvencia-tartományok a következők: a tekercsekkel 1 µT fölött és 50 kHz alatt tudunk mérni, az ME 3030B sugárzásmérővel 1 nT és 2 µT között illetve 16 Hz és 2 kHz között. A terepi méréseken a tekercseken fellépő feszültséget általában digitális multiméterrel mértük, de ezenkívül a feszültség időbeni változását vagyis a mágneses 10
10. ábra. Az ME3030B (piros) és a tekercsek (kék) amplitúdó- és frekvenciatartományai
tér jelalakját hordozható számítógépen futó digitális oszcilloszkóp-programmal [4] vizsgáltuk. Ez a program a számítógép hangkártyájával van közvetlen összeköttetésben, így a hangkártya bemenetére kötjük rá a tekercs két végét. A program képes megállapítani a jelben megjelenő frekvenciákat Fourier-transzformációval, ami jellemző a jel forrására.
11
2.
Helmholtz-tekercsek, biológiai vizsgálatok segítése
Az ELTE TTK Élettani és Neurobiológiai Tanszékén folytatnak kutatásokat azzal kapcsolatban, hogy milyen hatással vannak az alacsony frekvenciájú elektromágneses terek a patkányok agyára [5] illetve viselkedésére [6]. Mágneses térbe helyezett patkányok vizsgálatához a mágneses teret két áramjárta tekercs között hozzák létre egy egyszerű szerkezettel, melyen állítható nagyságú az 50 Hz-es váltakozó mágneses tér. Mi a karikák közti térrész bármely pontjában szeretnénk tudni a mágneses indukciót, ezt kiszámoljuk egységnyi áramerősség esetében, majd az ismert tekercsünkkel, a karikával számszerűsítjük.
11. ábra. A mágneses teret létrehozó tekercspár, és a szabályozója
2.1.
Helmholtz-tekercsek terének kiszámítása
A két tekercs Helmholtz-tekercsekhez hasonló módon van elhelyezve, de míg a valódi Helmholtz-tekercsek távolsága egyenlő a sugarukkal, itt egy kicsit távolabb van egymástól a két tekercs (12. ábra). A Helmholtz-tekercsekre az igaz, hogy a forgástengely mentén a mágneses térnek csak függőleges komponense van, és a nagysága a tekercspár között körülbelül homogén. Itt majd látni fogjuk, hogy ez nem pontosan teljesül, mert a forgástengely mentén a középsíktól a karikákhoz közeledve kis mértékben nő a mágneses tér amplitúdója. A tekercsek közötti mágneses tér nagyságának kiszámításához a Biot-Savarttörvényt alkalmazzuk: ~ = Iµ0 d~s × ~r (10) dB 4π r 3 A mágneses tér komponenseinek amplitúdóját a középponttól mért távolság függvényében számoljuk ki a középsíktól mért különböző távolságokban. A koordinátarendszert úgy vesszük fel hasonlóan 13. ábrához, hogy a z tengely megy át a tekercsek középpontján, és az origó a két tekercs közti felezősíkban a középpontban van. 12
12. ábra. A mágneses teret létrehozó tekercspár elhelyezkedése és adatai
13. ábra. A tekercsek terének kiszámításához használt koordináták [7]
A két karikából létrejött tér függőleges komponensét összeadhatjuk, így a (11) képletet kell kiszámolni a mágneses tér függőleges komponenséhez, ahol I az áramerősség, R a karikák sugara és h a karikák középsíktól mért távolsága, vagyis a távolságuk fele. NIRµ0 ( Bz (x, y, z) = 4π
Z
2π
0
+
(R − y sin ϕ − x cos ϕ)dϕ
3
((x − R cos ϕ)2 + (y − R sin ϕ)2 + (z − h)2 ) 2 Z 2π (R − y sin ϕ − x cos ϕ)dϕ 0
+
3
((x − R cos ϕ)2 + (y − R sin ϕ)2 + (z + h)2 ) 2
) (11)
A rendszer forgásszimmetriája miatt az x és az y tengely egyforma helyzetben van, ezért az y = 0 síkban az x függvényében kell kiszámolni a mágneses teret. Írtunk egy algoritmust, amiben paraméterként szerepel a z magasság, és 100 különböző x értékre (−2R < x < 2R-re) kiintegrálja a (11) kifejezést ϕ szerint. Ezt több magasságban kiszámoljuk, és ábrázoljuk az eredményt, így alakul ki a 14. ábra alakja egy szorzófaktor erejéig. A mágneses indukció függőleges komponensének nagyságát még számszerűsítenünk kell. Ezt a nagy karikával tettük meg, mert annak ismerjük minden adatát. 13
2.2.
Helmholtz-tekercsek terének mérése
Az összes mérést úgy végezzük, hogy a tekercsekhez tartozó szabályozógombot 1 mTra állítjuk be, amit a berendezés kiír, és majd a végén ellenőrizzük, hogy ez a kiírás mennyire igaz a szabályozó állásának többi értékénél. A karikát berakjuk a mágneses térbe a tekercsek közötti felezősíkba, így mérhetjük a függőleges komponens felületi integrálját, ami a karikán átmenő fluxus. Oszcilloszkóppal mérjük a karikán indukálódott feszültséget. A mért feszültség frekvenciája a hálózati 50 Hz, így mérhetjük a feszültség effektív értékét multiméterrel is, amit azonnal átszámolunk feszültség amplitúdóvá.
Mágneses indukció (mT)
2 z=0 cm 2 cm 4 cm 6 cm 8 cm
1.5
1
0.5
0
-0.5 -0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
Középponttól mért távolság (m) 14. ábra. A tekercsek közötti térben a mágneses indukció amplitúdója a forgástengelytől mért távolság függvényében a középsíktól különböző magasságokban A karikán a középsíkban 4,54 V-ot mértünk. A karikán átmenő fluxust a számolt 14. ábra alapján ki kell számolni. Ezt úgy számoljuk, hogy az ismert függvényértékeket megszorozzuk az adott középponttól mért távolságból számolt körkerülettel, majd ezeket összeadjuk a középpont és a karika sugara közötti részen és megszorozzuk a lépésközzel, vagyis két távolság különbségével. Ezt az összeget kell összehasonlítanunk a középpontban (x = y = z = 0) a mágneses indukció értékének és karika felületének a szorzatával, ami akkor állna fent, ha homogén lenne a tér. Ezt a számítást Octave-ban [8] elvégeztük, és azt kaptuk, hogy a homogén térrel számolva 1,16-szor nagyobb az amplitúdó, mint a karikán vett átlag, ezért a karikán mért feszültséget meg kell szorozni, ezzel az 1,16-tal, hogy a középpontban kapjuk meg a mágneses indukció függőleges komponensének nagyságát. Ezzel a számítással a középpontban (1,34±0,04) mT a mágneses indukció függőleges komponensének amplitúdója, ami az indukcióvektor amplitúdója, mert más irányú komponens nincs jelen. A tér többi pontjában a függőleges komponens amplitúdóját a számolt görbékkel határoztuk meg úgy, hogy megszoroztuk az összes görbét a fent meghatározott közös szorzófaktorral. Így kialakult a 14. ábra. 14
A számítások ellenőrzéséhez a kis tekercsünkkel mértünk a tekercsek közötti felezősíkban feszülséget a középponttól távolodva centiméterenként. Ez arra is jó volt, hogy az 1. fejezetben csak közelítőleg kiszámolt érzékenységet pontosítsuk. A feszültségek amplitúdóknak a felezősíkra kiszámolt görbére (14. ábra fekete görbe) kell illeszkedniük, ha megszorozzuk őket a kis tekercs érzékenységével. Erre az érzékenységre a középponti értékek alapján 4,13 adódik, vagyis a kis tekerccsel mért 1 V (vagy mV) feszültség 4,13±0,13 mT (vagy µT) mágneses indukciónak felel meg 50 Hz-es váltakozó tér esetén. Ez majdnem megegyezik a (6) képletből kapott menetszám alapján számolt érzékenységgel. Ezzel az értékkel beszorozva a mért feszültségeket ábrázoljuk a számolt görbénkkel együtt (15. ábra), és látjuk, hogy illeszkednek az adatsorok. Ugyanezt elvégeztük egy másik magasságban, a középsíktól 15,3 cm-re, aminek az eredménye a 16. ábrán látható. Ezek alapján az illeszkedések alapján elfogadjuk a számításainkat a tekercsek közötti mágneses térről. Ezzel egyúttal a kis tekercsünk érzékenységét is megismertük, amit a 3. táblázatban rögzítettünk.
Mágneses indukció (mT)
1.4 mert adatok z=0 cm
1.2 1 0.8 0.6 0.4 0.2 0 -0.2 -0.4
-0.2
0
0.2
0.4
Középponttól mért távolság (m) 15. ábra. A számolt (görbe) és a kis tekercsen mért (adatpontok) mágneses indukciók összehasonlítása a középsíkban
Karika Kis tekercs
0, 266 ± 0, 009 4, 13 ± 0, 13
µT mV µT mV
3. táblázat. A két tekercsünk érzékenysége Látható, hogy a tekercsek beállítóján lévő 1 mT-hoz képest a mágneses indukció amplitúdója nagyobb. Ezért ezt még tovább vizsgáltuk. Megmértük, hogy a szabályozó néhány különböző kiírása mellett a középpontban mekkora a mágneses tér, vagyis a kis tekercsen mekkora feszültség indukálódik, amit átszámítottunk mágneses térré, és rögzítettük a 17. ábrán. Az illesztett egyenes meredeksége 1, 03 ± 0, 01 és tengelymetszete −0, 37 ± 0, 03 mT, tehát a kiírás 15
Mágneses indukció (mT)
1.6 mert adatok z=15.3 cm
1.4 1.2 1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
Középponttól mért távolság (m) 16. ábra. A számolt (görbe) és a kis tekercsen mért (adatpontok) mágneses indukciók összehasonlítása a középsík felett 15,3 cm magasságban, vízszintes síkban
közelítőleg arányos a térrel a középpontban, de hozzá kell adni ezt a 0,37-et, vagyis a kiírás maga az amplitúdó csak eltolva. Az ábrán látszik, hogy kis terek esetén nem illeszkednek a pontok az egyenesre, tehát kis tereknél bizonytalan, hogy pontosan mekkora a mágneses tér a tekercsek között, ha csak a berendezés által kiírt értékre hagyatkozunk. 4 3.5
Kiírt (mT)
3 2.5 2 1.5 1 0.5 0 -0.5 0
0.5
1
1.5
2
2.5
3
3.5
4
4.5
Mért amplitúdó (mT) 17. ábra. A tekercsek szabályozóján kiírt érték a középpontban lévő mágneses indukció amplitúdójának függvényében
16
3. 3.1.
Elektromágneses terek környezetünkben Katódsugárcsöves monitor
Mint azt már említettük, oda kell figyelnünk, hogy mekkora frekvenciákat mérhetünk a tekercsekkel, mivel félrevezető értékeket kaphatunk. Egy mindennapjainkban gyakran előforduló sugárforrás, a monitor, aminek jele 30 kHz-es, és kis súllyal szerepel benne a felharmonikus is. A mágneses terének mérésekor például előfordult, hogy kétszer nagyobb frekvenciát láttunk az oszcilloszkópon, mint a monitor frekvenciája. Ez abban az esetben történik, ha a tekercsünkből és a kábel kapacitásából alkotott rezgőkör rezonanciafrekvenciájához közeli a monitor frekvenciájának első felharmonikusa. Oda kell figyelnünk, hogy milyen kábelellel kötjük össze a tekercset az oszcilloszkóppal, mert különböző kábelek esetén különböző a rezonanciafrekvenciája a tekercsből és a kábelből álló rezgőkörnek. Például egy körülbelül 3,5 m-es kábellel a rezonanciafrekvencia a monitor 30 kHz-es frekvenciájához esik közel, amitől csak a 30 kHz-es jelet látjuk az oszcilloszkópon, mert a 60 kHz-es felharmonikus már a rezonanciagörbe lecsengő ágán van. Egy rövidebb 80 cm-es kábellel a rezonanciafrekvencia a 60 kHz-es felharmonikushoz esik közel, ezért az oszcilloszkópon olyan, mintha frekvenciakétszereződést látnánk. Ez látható a 18. ábrán. Természetesen, ha jobban megnézzük az oszcilloszkópot, nem tökéletesen szinuszos a kétszeres frekvenciájú jel, de félrevezető.
18. ábra. A két szinuszgörbe összegének megjelenése az oszcilloszkópon
3.2.
Erősáramú földvezeték
A katódsugárcsöves monitorok nagy mágneses teret keltenek maguk körül, és működési elvükből következően nagyon érzékenyek a külső elektromágneses terekre is. Például, ha egy mágnesrudat túl közel rakunk a monitorhoz vagy tv-hez, akkor megváltoznak a színei. Illetve ha egy monitorhoz túl közel megy egy földvezeték, akkor 17
az is bezavarhat a képnek, és remegni kezd tőle. Egy ilyen probléma miatt kerültünk egy elég erős földvezetékhez Kispesten, aminek mérési eredményeit tárgyaljuk a következőkben. A Zalaegerszeg utcában húzódik egy földvezeték, aminek a földtől 1,2 m magasságban mérhető terét a 19. ábrán látható térképen rögzítettük. A vezeték tere kimagaslóan nagy a két oszlop között, amiket a 20. ábrán láthatunk. Tehát a két oszlop között a járda alatt húzódik a vezeték.
19. ábra. A mérés helyszíne és az utcában mért mágneses indukció értékek nT-ban A 20. ábrán látható oszlopnál digitális oszcilloszkóppal figyeltük meg a tekercsben indukálódott feszültséget. A jelalak a 21. ábrán látható. A digitális oszcilloszkóp frekvenciaanalízise alapján a jelalak főfrekvenciája 150 Hz. A jelalak nagyon nem szabályos, ami elrontja a mérések pontosságát. Meg tudjuk állapítani, hogy a földvezeték milyen mélyen van az aszfalthoz képest. A mágneses teret mérjük az aszfalttól mért magasság függvényében úgy, hogy a kis tekercsen mérjük az indukálódott feszültséget a földvezeték fölött. Mivel a mágneses indukció a vezetéktől mért távolság reciprokával arányos, ezért a feszültség reciprokát a magasság függvényében ábrázolva egyenest kapunk. A 22. ábrán látható ez az egyenes, amiről leolvasható az egyenes és az x tengely metszéspontja, ami a földvezeték távolsága az aszfalttól. Ez a mélység 44 ± 2 cm-nek adódott. A Kossuth utcában mérve a mágneses teret a vezetéktől mért merőleges távolság 18
20. ábra. Az oszlop és a tetején az összekötések
21. ábra. Az oszcilloszkópon megjelenő jelalak a földvezeték végénél lévő oszlopnál
függvényében a 23. ábrát kapjuk. Erre a görbére illeszthető egy elméleti görbe, vagyis hogy a mágneses tér a távolság reciprokával arányos, amiben már szerepel, hogy 44 cm mélyen van a vezeték. Azzal, hogy a távolság reciprokával arányos a tér azt igazoltuk, hogy egy áramjárta vezeték okozza a mágneses teret, aminek áramerősségére becslést tudunk adni az illesztés alapján. A 23. ábrán látható illesztésben 0I szereplő szorzófaktor felel meg a µ2π hányadosnak, amiből az I áramerősség amp19
2 1/U 1.8
1/U (1/mV)
1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0 -50
0
50
100
150
200
Magassaág (cm) 22. ábra. Földvezeték mélységének megállapítása
litúdójára körülbelül 10,3 A adódik. Ennek a földvezetéknek a néhány µT nagyságú mágneses tere már elég ahhoz, hogy a közelében lévő monitor képe remegjen. 1600
Mágneses indukció (nT)
B 1400 1200 1000 800 600 400 200 0 -120
-100
-80
-60
-40
-20
0
20
40
60
Távolság (m) 23. ábra. Földvezeték mágneses tere a távolság függvényében Az 23. ábrán látható, hogy a mérési pontoknak nagy hibájuk van, és nem túl pontosan illeszkednek az elméleti görbére. Ennek az az oka, hogy időben változik a műszeren kiírt mágneses tér. A kijelzőt felvettük videóra egy kis ideig, és a videóról felvettük a t-B adatpárokat, amiket a 24. ábrán láthatunk. Jól látható, hogy nagyon változik a mért tér egy rövid idő alatt, tehát nagyon nehéz pontos mérést végezni ilyen nem szabályos jelek esetén. 20
800
700
B (nT)
600
500
400
300
200
100
0 0
20
40
60
80
100
t (s)
120
140
160
180
24. ábra. Földvezeték mágneses terének időbeni változása az ME-3030B műszerrel mérve
3.3.
Nagyfeszültségű távvezeték
A közérdeklődés középpontjában leginkább a távvezetékek vannak, ezért végeztünk méréseket a Savoya Park mellett lévő távvezeték alatt, amiről 25. ábrán láthatunk egy képet.
25. ábra. A mért távvezeték egy oszlopa
21
Megnéztük az oszcilloszkópon a karikán indukálódott feszültség jelalakját. Ez látható a 26. ábrán. A főfrekvencia jól láthatóan a hálózati 50 Hz, ami mellett fellép egy erős zaj, ami elrontja a mérések pontosságát.
26. ábra. Az oszcilloszkópon megjelenő jelalak a távvezeték alatt A távvezeték vonalára merőlegesen 5 m-enként mérjük a mágneses tér függőleges illetve vízszintes komponensét. A távvezetékkel párhuzamos irányban nem lép fel a távvezetékből adódó mágneses tér. A mért értékeket összehasonlítjuk az elvi számításainkkal, amit a 27. ábrán láthatunk. Az számítások alapjául azt vesszük, hogy a 6 vezetékből adódó mágneses tér szuperpozícióját mérjük. A távvezetékben 3 fázisban folyik az áram 120◦-os fáziskülöbségekkel, ezeket a fázisokat jelölik R, S, T betűkkel. Feltételezzük, hogy a vezetékekben egyforma amplitúdójú áram folyik. A következő módon számoljuk a mágneses tér nagyságát: a különböző fázisok amplitúdóját elnevezzük a, b, c betűkkel, így a következő módon adódnak össze a fázisok: a sin ϕ + b sin(ϕ + 120◦ ) + c sin(ϕ − 120◦) = √ √ 3 3 1 1 b cos ϕ − c sin ϕ − c cos ϕ = = a sin ϕ − b sin ϕ + 2 2 2 √ 2 √ 3 b c 3 = (a − − ) sin ϕ + ( b− c) cos ϕ = 2 2 2 2 = A cos δ sin ϕ + A sin δ cos ϕ (12)
22
Minket az A amplitúdó érdekel: A2 = a2 +
b2 c2 bc 3b2 3c2 3bc + − ab − ac + + + − = 4 4 2 4 4 2 = a2 + b2 + c2 − ab − ac − bc (13)
A különböző fázisok amplitúdói a vezetékek elhelyezkedésétől függenek. Egyenes vezető mágneses tere a vezetéktől való távolsággal fordítottan arányos a következő képlet alapján, ahol I az áramerősség. B= A
µ0 I 2π r
(14)
µ0 I 2π
kifejezést elnevezzük k konstansnak. A mágneses tér vízszintes (x) komponensének amplitúdója a középvonaltól mért távolság függvényében a következőképpen számolható, ahol h a megfelelő vezeték magassága, l pedig a megfelelő vezeték középvonaltól mért távolsága a megfelelő előjellel. h Bx (x) = k (15) (x + l)2 + h2
A függőleges (y) komponens amplitúdója az előzőhöz hasonló módon: By (x) = k
x+l (x + l)2 + h2
(16)
A hat vezetékből adódó mágneses tér amplitúdóját végül úgy számoljuk, hogy a két-két egyforma fázisú vezeték adott irányú terét összeadjuk, majd a (13)-nek megfelelően számoljuk az amplitúdót az adott irányban. A számítások eredménye a 27. ábrán látható a mérési pontokkal együtt. A k konstansra körülbelül 70000 nT m adódik, vagyis ezzel a számmal illeszkednek legjobban az adatpontok a görbékre. Ebből a konstansból számolható körülbelül a vezetékek árama, amire 350 A adódik. A mért adatpontok alakjából a fáziskiosztásra is következtethetünk, mert ha máshogy vannak kiosztva a fázisok nagyon eltérő eredményeket kaphatunk. Ha úgy osztjuk ki a fázisokat a számolásban, hogy az azonos fázisú vezetékek mind párhuzamosak, akkor a 28. ábrához jutunk, amin jól láthatóan az origóban a függőleges komponensnek nullához kéne tartania, de a mérési eredmények szerint nem így van, tehát ez a fáziskiosztás rossz. Az adatpontokra legjobban illeszkedő görbékhez (27. ábra) úgy jutunk, ha a felső két vezeték fázisa megegyezik, az alsó kettő pedig fel van cserélve.
23
1400 By Bx |B|
1200
B (nT)
1000 800 600 400 200 0 -60
-40
-20
0
20
40
60
x (m) 27. ábra. Távvezeték mágneses terének két komponensének és abszolut értékének számolt görbéje és mért pontjai a legjobban illeszkedő fáziskiosztással
1400 By Bx |B|
1200
B (nT)
1000 800 600 400 200 0 -60
-40
-20
0
20
40
60
x (m) 28. ábra. Távvezeték mágneses terének két komponensének és abszolut értékének számolt görbéje és mért pontjai rossz fáziskiosztással
3.4.
Közlekedési eszközök
Végeztünk méréseket villamosokon és HÉV-en. Ezeket a méréseket egyelőre nem tudjuk maximális pontossággal elvégezni, mert ezek a közlekedési eszközök egyenárammal működnek, és mi csak a változó tereket érzékeljük a jelenlegi műszereinkkel. Viszont a villamos motorjának közelében mérhető a motor működésénél keletkező változó mágneses tér. 24
A közlekedési eszközök, és a felsővezetékük változó mágneses tere általában kicsi, néhány µT, ezért az me3030B műszerrel mértünk legtöbbször, de a tekercsekkel és az oszcilloszkóppal megnéztük, hogy milyen jeleket mérünk. A HÉV és a villamos esetében is a jelek frekvenciája 300 Hz-es főfrekvenciából és a felharmonikusaiból áll. Erre látható egy példa a 29. ábrán, ahol a mozgó villamos mellett a megállóban mértük ezt a 300 Hz-es jelet. Ez a jel valószínűleg a feszültség egyenirányítása után marad rajta az egyenfeszültségen az egyenirányítás jellegéből adódóan.
29. ábra. Az oszcilloszkópon megjelenő jelalak a villamoson Az oszcilloszkóppal való mérésekkor a villamosnál és a HÉV-nél is jól látható volt, hogy a gyorsuló jármű mellett megnő az indukálódott feszültség. Ez abból adódik, hogy nő a jármű teljesítménye, ezzel a felvett áramerősség, amivel arányos a mágneses tér nagysága, és ez a növekedés feszültséget indukál a tekercsben. A másik műszerrel egyszerűbben mérhetünk a járműveken, úgyhogy ezzel több mérést tudtunk végezni. Az ilyen méréseket úgy végeztük, hogy a műszert olyan helyzetben tartottuk, hogy a maximális legyen a tér erőssége, és így felvettük videóra a kijelzőt. A videókról utólag felvettük az idő-mágneses indukció adatpárokat, amiket ábrázolhatunk. Az első mérést a 18-as villamoson végeztük a Savoya Park felé. A 30. ábrán látható, hogy hogyan változik a tér gyorsításkor, és hogy menynyire lecsökken, amikor a villamos a megállóban áll. A 31. ábrán kinagyítva látszik két jellemző indulástól megállásig tartó szakasz. Jól látható, hogy a gyorsítás alatt megnő a mágneses tér, a lassítás elején egy visszacsatolásból származó csúcs jelenik meg és a fékezésnél több kisebb csúcs látható.
25
30. ábra. A 18-as villamos útja során mérhető mágneses indukció
31. ábra. A 18-as villamoson mért jellemző csúcsok
26
4.
Összefoglalás
A törvényben szabályozott határértékeket (100-500 µT) ritkán éri el bármi a hétköznapi életben, de sokkal alacsonyabb, nem feltétlenül hivatalos ajánlások is megjelennek az irodalomban (20 nT, [9]). Az egészségügyi hatások tisztázása és a megengedett határértékek ismerete fontos ahhoz, hogy a lakosság félelmeit feleslegesen ne növeljük, melyeket a nagyon erzékeny műszerek által már kimérhető elektromágneses terek esetleg előidéznek. Ezek az érzékeny eszközök mindemellett nagyon hasznosak a feleslegesen jelen lévő, könnyen megszüntethető elektromágneses tér forrásainak felderítésében, amíg az egészségügyi hatások egyértelműen nem tisztázódnak. Fontos, hogy az elektromágneses tereknek műszaki hatását érzékeljük, de biológiai hatásukról még nagyon keveset tudunk, ezért is hasznos, hogy tudjuk ezeket a tereket mérni egyszerű, hitelesített eszközeinkkel. Ezzel tudjuk segíteni a biológiai kutatásokat, melyek az elektromágneses terek viselkedésre gyakorolt hatását vizsgálják. Ezeknek az eredményei között kevés a szignifikáns eltérés a 500 µT-ás mágneses térrel kezelt illetve a nem kezelt csoportok viselkedése között [6]. A kísérletek tovább folynak, nagyobb amplitúdójú mágneses terek alkalmazásával. A mérési tartományok, főleg a mérhető frekvenciatartomány kiszélesítésével lehetne folytatni a mérőeszközök fejlesztését. Ezen kívül minél több helyen, minél több mérést kell elvégezni, hogy jobb képet kapjunk a környezetünkben előforduló alacsonyfrekvenciás elekromágneses sugárzásokról. A mérések folytatásának egy másik iránya, hogy az egyes emberekre ható expozícióját mérjük a különböző frekvenciájú elektromágneses sugárzásoknak, vagyis hogy az egyes emberre pl. egy nap alatt mekkora terhelés jut. Ilyen mérésekhez más jellegű eszközök kellenek, amik rögzítik a különböző frekvenciájú sugárzások intenzitását bizonyos időközönként. A rögzített adatokat egy a mérés során írt naplóból lehet kiértékelni. Ilyen jellegű mérések is folynak az Atomfizika Tanszéken [10].
27
Hivatkozások [1] Guidelines for Limiting Exposure to Time-Varying Electric, Magnetic, and Electromagnetic Fields (up to 300 GHz). Health Physics 74 (4): 494-522; 1998. (http://www.icnirp.de/PubEMF.htm) [2] Thuróczy György, Szabó Judit, Bakos József: Hálózati frekvenciájú elektromágneses terek környezetünkben (Széchenyi Füzetek - 5. számú útmutató az egészség megőrzéséhez), OSSKI, 2004. [3] http://en.wikipedia.org/wiki/Inductor [4] http://www.zeitnitz.de/Christian/Scope/Scope_en.html [5] Petra Varró, Renáta Szemerszky, György Bárdos, Ildikó Világi: Changes in synaptic efficacy and seizure susceptibility in rat brain slices following extremely low-frequency electromagnetic filed exposure, Bioelectromagnetics 2008. (revízió alatt) [6] Szemerszky Renáta: Hosszú-tartamú alacsony frekvenciás elektromágneses expozíció viselkedésélettani hatásai patkányban, ELTE Élettani és Neurobiológiai Tanszék, Szakdolgozat 2007. (kézirat) [7] Kovács István, Párkányi László: Fizikai példatár II. (438. feladat) [8] http://www.gnu.org/software/octave [9] http://eszmog.hu [10] Finta Viktória: Milyen hatásai vannak a környezetünkben lévő nem ionizáló elektromágneses sugárzásoknak? Fizikai Szemle 2007/11.
28