Profil georadarowy
![](http://upload.wikimedia.org/wikipedia/commons/thumb/7/7f/Profil_georadarowy_rzeki.gif/360px-Profil_georadarowy_rzeki.gif)
Profil georadarowy (także radargram, falogram) w badaniach georadarowych jest graficznym obrazem badanego podłoża. Profil georadarowy złożony jest z hodografów odbitej w podłożu fali elektromagnetycznej, która w przypadku większości sprzętu spolaryzowana jest liniowo[1]. Poszczególne wydarzenia są tu graficznie przedstawiane jako amplituda i średnia częstotliwość odbitego sygnału w poziomej skali odległości i pionowej skali czasu (pomiędzy emisją a rejestracją powracającego odbitego sygnału). Przy znajomości prędkości fali elektromagnetycznej w badanym ośrodku skala czasu może być przekształcona na głębokość[2].
Odbicie sygnału
![](http://upload.wikimedia.org/wikipedia/commons/thumb/3/30/Partial_transmittance.gif/360px-Partial_transmittance.gif)
![](http://upload.wikimedia.org/wikipedia/commons/thumb/0/04/Georadar-odbicie_fali.jpg/360px-Georadar-odbicie_fali.jpg)
Rozchodzenie się fal elektromagnetycznych w gruncie, który ma właściwości dielektryczne następuje w wyniku indukcji i polaryzacji elektrycznej ośrodka. Gdy fala dociera do granicy pomiędzy warstwami o różnej względnej przenikalności elektrycznej r część energii ulega odbiciu a część zostaje przekazana w dół. Odbicie następuje tu zgodnie z zasadą Huygensa, tzn. poprzez konstruktywną interferencję pomiędzy falami, które uległy dyfrakcji na sąsiednich punktach położonych na granicy pomiędzy ośrodkami o różnych właściwościach dielektrycznych. Podobnie jak w przypadku fali sejsmicznej, amplituda sygnału na hodografie odbitej fali elektromagnetycznej jest zależna od fresnelowskiego (amplitudowego) współczynnika odbicia R[3]:
W przypadku fali elektromagnetycznej, rozchodzącej się w ośrodku niepochłaniającym fali i niemagnetycznym, współczynnik załamania Fresnela zależy od przenikalności elektrycznej ośrodków[4]:
Wartość R będzie więc zawarta pomiędzy −1 a +1. Przy czym wartości ujemne oznaczają zmianę fazy sygnału odbitego w stosunku do padającego. Współczynnik odbicia określa więc na profilu georadarowym fazę odbitej od granicy litologicznej fali. Przykładowo granica pomiędzy suchym a nawodnionym piaskiem charakteryzuje się ujemną wartością współczynnika odbicia Fresnela (co także wiąże się ze wzrostem przenikalności elektrycznej poniżej tej granicy), podobna sytuacja będzie miała miejsce w przypadku wzrostu zagęszczenia suchego piasku. Jeżeli jednak piasek jest nawodniony, nagły wzrost zagęszczenia cechuje się dodatnim współczynnikiem odbicia Fresnela (pozytywną fazą odbitej i zarejestrowanej fali). Wynika to ze zmniejszenia zawartości płynów porowych (czyli wzrostu przenikalności elektrycznej na granicy ośrodków).
Fala odbita od powierzchni terenu zawsze charakteryzuje się negatywnym współczynnikiem R (powietrze ma zawsze niższą wartość przenikalności elektrycznej niż grunt), a więc faza zarejestrowanej fali jest zawsze negatywna[5]. Podobnie górna granica (strop) warstwy torfu zalegającej pomiędzy dwoma warstwami nawodnionego piasku będzie zaznaczona poprzez dodatnią fazę fali (sygnał przechodzi w dół przez granicę, która charakteryzuje się spadkiem przenikalności elektrycznej), natomiast dolna granica (spąg) poprzez negatywną fazę fali (wzrost przenikalności elektrycznej na granicy)[6][7].
Istnieją różnice w prezentowaniu fazy fali na profilach georadarowych. Dla przykładu przy interpretacji profili radarowych na potrzeby drogownictwa powszechnie uznawane od lat 90. standardy[8] zalecają by pozytywna faza fali odbitej prezentowana była w skali czarno-białej za pomocą koloru białego. Dlatego powierzchnia terenu będzie reprezentowana na profilu jako czarna linia. W innych jednak zastosowaniach (badaniach naukowych, archeologii i innych) kolor czarny oznacza pozytywną fazę i dlatego powierzchnia terenu oznaczana jest kolorem białym. Podobnie też linia o kolorze białym reprezentuje zwierciadło wód gruntowych.
Prędkość fali elektromagnetycznej
W geofizyce prędkość fali w danym ośrodku ma bardziej charakter wartości skalarnej niż wektorowej, to znaczy opisuje ona charakter ośrodka. Jest więc ona niezależna od częstotliwości wysyłanego sygnału elektromagnetycznego, czyli zastosowanej anteny (pozwala to m.in. na stosowanie kilku anten o różnych częstotliwościach jednocześnie). Niemniej w ośrodkach anizotropowych prędkość fali elektromagnetycznej może w pewnym stopniu przybierać różne wartości w różnych punktach i kierunkach. W skałach ogólnie istnieje tendencja spadku średniej prędkości fali elektromagnetycznej wraz z głębokością (w przeciwieństwie do prędkości fali sejsmicznej) co wiąże się ze wzrostem średniej przenikalności elektrycznej.
Różnica czasu powrotu (czas od wysłania do zarejestrowania) fali odbitej w podłożu i fali powierzchniowej t od jej prędkości V jest opisywana wzorem[9]:
gdzie:
- – głębokość przeszkody, na której zachodzi odbicie,
- – rozstaw anten.
Przy badaniach, w których stosuje się jedną antenę nadawczo-odbiorczą, lub przy większych głębokościach, gdy odległość anten nie mają znaczenia, równanie ma więc postać:
Równanie to może być wykorzystywane do obliczania średniej prędkości fali w określonym ośrodku, gdy znana jest głębokość jakiegoś wydarzenia na profilu georadarowym (na przykład znana jest głębokość zwierciadła wody gruntowej).
W badaniach common-offset przy zmianie prędkości w ośrodku z na (np. przy przejściu przez zwierciadło wody gruntowej), ze względu na zniekształconą geometrię przebiegu fali, ta zależność jest bardziej złożona. Zmiana własności ośrodka działa podobnie jak soczewka w optyce. Prędkość fali elektromagnetycznej w danym ośrodku o niskiej utracie energii (opisanej poniżej) jest opisywana zależnością[4]:
gdzie:
- – prędkość światła w próżni (około 0,30 m/ns),
- – względna przenikalność elektryczna warstwy.
Równanie to może być wykorzystywane do obliczeń prędkości fali w ośrodku, którego właściwości elektryczne możemy zbadać laboratoryjnie (np. próbki betonu, skały itp.).
Częstotliwość a rozdzielczość
Podobnie jak w badaniach sejsmicznych, pionowa rozdzielczość w metodzie georadarowej jest ograniczona do przedmiotów o pionowej średnicy większych od jednej czwartej długości fali:
gdzie: jest długością fali elektromagnetycznej i jest wyrażana jako:
gdzie:
- – dominująca częstotliwość,
- – prędkość fal w danym ośrodku, czyli:
Dlatego też wyższe częstotliwości wykorzystuje się do wykrywania mniejszych obiektów. Fala o częstotliwości 250 MHz, przy założeniu prędkości 0,1 m/ns jest w stanie prawidłowo przedstawić obiekt o miąższości 10 centymetrów. Jednak fale o wyższej częstotliwości są silniej tłumione, w wyniku czego zakres częstotliwości fal powracających jest zubożony o fale o wyższej częstotliwości i średnia częstotliwość i rozdzielczość falogramu jest obniżona.
Znajomość spektrum częstotliwości powracającej fali jest także przydatna przy przetwarzaniu danych georadarowych, tzn. do usuwania tych jej zakresów, które są związane z szumem i wzmacnianiu sygnału. Jest to możliwe dzięki zastosowaniu transformacji Fouriera, która przekształca wykres amplitudy w skali czasu w wykres amplitudy w skali częstotliwości.
Tłumienie sygnału
Rozchodzenie się sygnału w ośrodku pochłaniającym falę można opisać przez użycie zespolonego współczynnika załamania, w którym składowa urojona określa pochłanianie fali i jest interpretowana jako opór elektryczny lub hamowanie ruchu wzbudzonych falą spolaryzowanych ładunków (dipoli). Składowa rzeczywista, jest związana z przenikalnością elektryczną ośrodka jest określana jako zdolność materiału do magazynowania ładunku elektrycznego. Jest to ilość energii jaka przy zadanym polu elektrycznym jest magazynowana w postaci spolaryzowanych ładunków (polaryzacja orbit elektronów, cząsteczkowa, dipolowa oraz polaryzacja jonowa).
Przy niskiej częstotliwości, ładunki w ośrodku mają czas w pełni zareagować na pole elektryczne, największa ilość energii jest magazynowana w postaci spolaryzowanych cząstek a straty są niewielkie. Przy wyższych częstotliwościach czas zmiany pola elektrycznego jest mniejszy, polaryzacja następuje jedynie częściowo i większa część energii jest tracona wskutek oporu ośrodka. Dlatego też fale o niskich częstotliwościach ulegają mniejszemu tłumieniu. Utrata energii może być także związana z wysoką przewodnością elektryczną ośrodka na skutek transportu ładunków (np. w wodzie morskiej) lub wysoką magnetycznością ośrodka (ferromagnetyki: żelazo, kobalt, nikiel)[4].
Przypisy
- ↑ G.S. Baker, H.M. Jol „Stratigraphic analyses using GPR” Geological Society of America, 2007.
- ↑ Practical processing of GPR data, Sensors & Software Inc., 1999.
- ↑ K. Żukowska, T. Wiktorczyk: Wyznaczanie współczynnika załamania cienkich warstw dielektrycznych metodą spektrofotometryczną. [dostęp 2009-06-01].
- ↑ a b c Neal, A. „Ground Penetrating Radar and its use in sedimentology: principles, problems and progress” Earth Science Reviews 66, 2004.
- ↑ Subsurface Views Vol 3, No. 10 April 2007 Sensors & Software Inc. [dostęp 11:30 24.09.2009].
- ↑ Unsworth, M., „Ground-penetrating radar (GPR)” Geophysics 223, Marzec 2009.
- ↑ Agron J., Gogineni S.P., „One-Dimensional Inverse Scattering: Reconstruction of Permittivity Profiles for Stratified Dielectric Media” Technical Report, Online.
- ↑ Saarenketo, T., Electrical properties of road materials and subgrade soils and the use of GPR in traffic infrastructure surveys online [dostęp 11:10 24.09.2009].
- ↑ Fisher, S.C., Stewart, R.R., Jolt, H.M. Processing Ground Penetrating Radar (GPR) Data, CREWES Research Report Volume 4, 1992.
Bibliografia
- Daniels, J.J., Ground Penetrating Radar Fundamentals, Online: Appendix to a Report to the U.S.EPA, Region V, 2000.
- Fisher, E., McMechan, G.A., Annant, P., Cosways, S. Examples of reverse-time migration of single-channel, ground-penetrating radar profiles, Geophysics 1992.