Nowy tydzień - stary problem -Geoelektryka

Dziś obiecane dokończenie ostatniego posta. Zacznę może jeszcze raz.Po co to robimy?? 

Celem jest określenie przestrzeni zalegania bariery zrobionej z DSM o raz na ile to możliwe wskazanie (jeżeli takie istnieją) rozluźnień wewnątrz kolumny. 

Proste?? - chyba tak, zwłaszcza że znam geologię z przed wykonania kolumn. 

Jak widzicie na profilu geologicznym występują gliny rzeczne. Pod względem litologicznym wykształcone są one jako pyły, pyły piaszczyste i gliny pylaste, gliny oraz gliny pylaste zwięzłe.
Gliny rzeczne są zróżnicowane pod względem stanu konsystencji, dominuje jednak stan plastyczny. Lokalnie stwierdza się występowanie gruntów miękkoplastycznych.
Wśród warstwy glin rzecznych, na głębokości 10,5 - 12,0 m ppt  stwierdzono warstwę piasków średnich przewarstwionych gliną pylastą. 

Znając budowę geologiczną i wiedząc że pod wpływem mieszania cementu z gruntem powstanie ośrodek o stosunkowo odmiennej oporności elektrycznej (od otaczającego tła jakim są gliny rzeczne), potrafie przewidzieć jaki otrzymam obraz z pomiarów terenowych. 
Zaprezentuje to na przykładzie prostego modelowania.

Zakładam oporność tła(glin) na poziomie 60 Ohm, a cementogruntu na poziomie około 150 Ohm, przy odpowiednim rostawie rzędu 120m, powinienem osiągnąć max głębokość pomiaru około 15-20m. Dodaję delikatne szumy, i otrzymuję taki oto obraz:

Nic trudnego, i mając trochę obycia z geoelektryką łatwe do przewidzenia. Ale co można powiedzieć o wyniku ? Po pierwsze primo, na pierwszy rzut oka wygląda bardzo ładnie. Interacja zatrzymana po trzecim obliczeniu przy błędzie na poziomie 2.4% - co oznacza że zarejestrowane(modelowane) dane są bardzo dobrej jakości.  Po drugie primo, znając mniej więcej położenie bariery zastosowany układ nie jest wstanie "dokladnie" określić granic pionowych bariery. I po trzecie primo - ultimo ! przy zastosowanym rozstawie elektrod, rozdzielczość metody w tym wypadku jest zbyt słaba.
W teorii tyle, teraz czas na pomiar terenowy ;)
Badania przeprowadzone zostały przy + temperaturze i miłej atmosferze ;) Obeszło się bez problemów i tzw rzucania błędami. Okład zastosowany to "pełny wenner" co oznacza dwie rzeczy:
-pomiar będzie trwał długo
-osiągniemy maksymalną możliwą głębokość przy takim rozstawie elektrod (2m)

 Ta czerwona kreska na rysunku to nasza linia pomiarowa, a te kółka to kolumny DSM. Środek rozstawu przypada dokładnie na środku kolumny.
Tyle o terenie, teraz trochę o wynikach.

Po odpaleniu programu i szybkiej obróbce, otrzymałem ładny kolorowy obraz, ale ku mojemu zaskoczeniu, zupełnie inny niż zakładałem - wtf?? Otrzymałem co prawda zarys kolumny ale spodziewałem się dokładniejszego odwzorowania, zwłaszcza dna kolumny. A co najważniejsze spodziewałem otrzymać wysokooporowy ośrodek pośród niskooporowego tła!
Nasuwa się pytanie co tak naprawdę się stało ? Czy otrzymany wynik jest poprawny, czy może należy powtórzyć badanie? Jak to interpretować?
Zajmijmy się samym kształtem.
Mozolne przetwarzanie w programie, próby zastosowań innych metod obliczeniowych i nadal to samo. Nagle, zaświeciło mi w głowie - odpalę starą demówkę, jeszcze ze studiów, res2dinv - zobaczymy jak ten program sobie poradzi. Program zatrzymuje się po 3 interacjach i wynik jest jeszcze gorszy - szkoda, ale nie poddaję się, szukam starych projektów i znajduję ;) - "modelowanie i inwersja w res2dinv". Eureka !! -opcja "robust inversion" - zastosowanie do pionowych granic. Tak tak, to jest to. Jaka różnica ;) Szkoda że to tylko demo. Program res2dinv poradził sobie znacznie lepiej, możemy już z dużym przybliżeniem wskazać granicę kolumny, a zwłaszcza spąg.

Pozostaje jeszcze kwestia niskiej oporności bariery. Jak to zinterpretować ? - w tej chwili ciężko cokolwiek o tym powiedzieć, a nie chce wysuwać pochopnych wniosków. Jasne jest że kolumna nie jest jednorodna. Pozostaje tylko pojechać tam jeszcze raz i zastosować inną metodę geofizyczną.
Taki o to wydawało by się prosty temat, a dał nam więcej pytań niż odpowiedzi.
Teraz pozostaje tylko czekać na dobre warunki atmosferyczne, jechać w teren i rozpracować ten obiekt za pomocą metod sejsmicznych.