WPROWADZENIE do Cwiczenie 4, Geologia GZMiW UAM 2010-2013, II rok, Geofizyka, CD-2 2012, Projekt - 2


Materiały uzupełniające do Ćwiczenia 4

Ilościowa interpretacja sondowań elektrooporowych w oparciu o model horyzontalnie warstwowanej przewodzącej półprzestrzeni (tzw. model 1D) - z wykorzystaniem programu INT_PSE lub IPI2win

WSTĘP - Na czym polega interpretacja ilościowa sondowań elektrooporowych ?

Oporność pozorna ρa jest parametrem geoelektrycznym, który charakteryzuje niejednorodny ośrodek geologiczny znajdujący się w „polu widzenia” układu pomiarowego. Można traktować ją jak pewną wypadkową oporność elektryczną niejednorodnego ośrodka, w którym przestrzenny rozkład oporności „rzeczywistej” narzucony jest przez budowę geologiczną. Oporność pozorna jest zależna od:

W przypadku horyzontalnie warstwowanego przekroju geoelektrycznego można

ten rozkład opisać podając oporności i miąższości warstw występujących w profi-

lu: ρ1, h1, ρ2, h2,..., ρN-1, hN-1 , ρN (ostatnia N-ta „warstwa” na nieskończenie dużą

miąższość)

symetrycznego układu pomiarowego AMNB).

ρa = f {ρ(x, y, z) , Pp(xp, yp, zp) , AB/2}

Sondowanie elektrooporowe polega na określaniu oporności pozornej ρa w stałym punkcie pomiarowym (Pp const.) przy zwiększających się rozstawach AB układu pomiarowego. Im większy rozstaw układu tym większe jest jego „pole widzenia”(większy zasięg penetracji).

ρa = f {ρ1, h1, ρ2, h2, ... , ρN-1, hN-1 , ρN, (AB/2)i } przy Pp const.

Krzywą sondowania - wykres zależności ρa = f (AB/2) można otrzymywać w wyniku:

(rozkład oporności: ρ1, h1, ρ2, h2, ... jest nieznany, mamy go ustalić)

(rozkład oporności jest znany , krzywa sondowania obliczana jest dla założo-

nego rozkładu oporności: [ρ1,h1,ρ2,h2,...,ρN-1,hN-1,ρN] - model 1D )

Używane są terminy: polowa krzywa i teoretyczna (lub modelowa) krzywa sondowania elektrooporowego.

Interpretacja ilościowa sondowania (w oparciu o model 1D) polega na określeniu profilu geoelektrycznego, w miejscu gdzie wykonano sondowanie, na podstawie pomierzonej krzywej polowej ρaP = f (AB/2). Ten etap badań nazywany jest interpretacją geofizyczną, a wyinterpretowany profil geoelektryczny opisany jest przez podanie oporności i miąższości wyodrębnionych warstw: ρ1,h1,ρ2,h2,...,ρN-1,hN-1,ρN. Kolejny etap badań to interpretacja geologiczna - poszczególnym warstwom przypisywany jest charakter litologiczny na podstawie ich oporności. Efektem końcowym interpretacji pojedynczego sondowania jest zgeneralizowany profil geologiczny (należy zaznaczyć, że jest to profil otrzymany z sondowania geoelektrycznego a nie np. z wiercenia).

W większości programów komputerowych (np.: INT_PSE lub IPI2win), który stanowi „narzędzie” do interpretacji ilościowej sondowań elektrooporowych, zastosowano tzw. metodę doboru. Istotą tej metody jest porównywanie pomierzonej polowej krzywej sondowania z krzywą teoretyczną obliczaną dla zadawanego przez interpretatora modelu (deklarowana jest ilość warstw: N oraz ich oporności i miąższości: ρ1,h1,ρ2,h2,...,ρN-1,hN-1,ρN). Przyjęte są też pewne kryteria - wskaźniki oceniające podobieństwo („identyczność”) krzywych: teoretycznej i polowej. Jeżeli te kryteria nie są spełnione wprowadzana jest korekta modelu, obliczana jest nowa krzywa teoretyczna i porównywana z krzywą polową. Czynność ta jest powtarzana (procedura iteracyjna wykonuje kolejne kroki) aż do „dopasowania” krzywej teoretycznej do krzywej polowej. Jeśli osiągniemy zadawalające dopasowanie to można uznać, że ostatni z użytych modeli odpowiada profilowi geoelektrycznemu w miejscu, gdzie pomierzono krzywą polową.

A więc szukamy „identyczności” funkcji:

ρaP = f (AB/2)i ρaT = f {ρ1, h1, ρ2, h2, ... , ρN-1, hN-1 , ρN, (AB/2)i }

Za wskaźnik dopasowania krzywej teoretycznej do krzywej polowej, w programie INT_PSE, przyjęto średnią kwadratową różnicę względną E między tymi krzywymi:

0x08 graphic

gdzie: k - ilość punktów na krzywej = ilość rozstawów

W procedurze dopasowywania szukane są takie wartości zmiennych niezależnych (ρ1,h1,ρ2,h2,..., ρN-1, hN-1, ρN ), przy których wartość funkcji E jest minimalna.

E = f(ρ1, h1, ρ2, h2, ... , ρN-1, hN-1 , ρN) minimum

Jest to problem matematyczny, badania minimum funkcji wielu zmiennych i w programie komputerowym da się „zautomatyzować”. Znaczy to, że interpretator po analizie wizualnej kształtu krzywej polowej (wyświetlanej na monitorze) określa tylko:

ich oporności i miąższości: ρ1,h1,ρ2,h2,...,ρN-1,hN-1,ρN),

( średnia kwadratowa różnica względna mniejsza od 0.1%),

natomiast dalsze poprawianie modelu aż do dopasowania krzywych przejmuje komputer (procedura gradientowa określa, którą zmienną, o ile i w jakim kierunku zmieniać by funkcja

E(ρ1, h1, ρ2, h2, ... , ρN-1, hN-1 , ρN) dążyła do minimum). Powyższe uwagi odnoszą się do programu o nazwie: INT_PSE.

W IPI2win, innym z udostępnianych w ramach ćwiczeń programie na interpretację ilościową krzywych sondowań elektrooporowych w oparciu o model 1D, określenie tzw. MODELU STARTOWEGO można „scedować” na sam program, wybierając odpowiednią opcję w menu (patrz plik Guide_IPI2win.pdf na dysku CD-2_2010, sam program jest na CD-1_2010)

Przykład interpretacji krzywej sondowania pokazany jest w pliku: Instrukcja - opis menu - przykład.doc w folderze: INT_PSE na CD-2_2010

Naukę posługiwania się programem IPI2win proponuję przeprowadzić na tym samym przykładzie. Metodyka posługiwania się tym „narzędziem” - programem IPI2win - przedstawiona została na ćwiczeniach. Opis menu programu i podstawowe zasady jego używania podane są w pliku: Guide IPI2win.pdf folder: IPI2win na CD-2_2010 (oraz na CD-1_2010).

Janusz Antoniuk 1, 2, e-mail: antoniuk@geol.agh.edu.pl

1. Zakład Geofizyki, Wydział Geologii, Geofizyki i Ochrony Środowiska, AGH w Krakowie

2. Zakład Geologii Środowiskowej, Instytut Geologii, UAM w Poznaniu

3

0x01 graphic



Wyszukiwarka

Podobne podstrony:
Cwiczenie 4 dla gr 1, Geologia GZMiW UAM 2010-2013, II rok, Geofizyka, CD-2 2012, Projekt - 2
Cwiczenie 4 dla gr 1, Geologia GZMiW UAM 2010-2013, II rok, Geofizyka, CD-2 2012, Projekt - 2
UAM-wyniki pomiarów podatności magnetycznej, Geologia GZMiW UAM 2010-2013, II rok, Geofizyka, CD-1
Projekt 2B, Geologia GZMiW UAM 2010-2013, II rok, Geofizyka, Projekt 2, B
Projekt 4, Geologia GZMiW UAM 2010-2013, II rok, Geofizyka, Projekt 4
sedymentologia3 sko 02.03, Geologia GZMiW UAM 2010-2013, II rok, Sedymentologia, Wykłady, wykłady 1-
EGZAMIN - pojęcia 02, Geologia GZMiW UAM 2010-2013, II rok, Hydrogeologia, EGZAMIN
EGZAMIN - pojęcia 01, Geologia GZMiW UAM 2010-2013, II rok, Hydrogeologia, EGZAMIN
pochodzenie wód podziemnych, Geologia GZMiW UAM 2010-2013, II rok, Hydrogeologia, Materiały
sedymentologia1 rot 23.02, Geologia GZMiW UAM 2010-2013, II rok, Sedymentologia, Wykłady, wykłady 1-
sedymentologia4 sko 03.03, Geologia GZMiW UAM 2010-2013, II rok, Sedymentologia, Wykłady, wykłady 1-
Omówienie prezentacji RZGW Warszawa, Geologia GZMiW UAM 2010-2013, III rok, Gospodarowanie zasobami
QUATERNARY EVOLUTION OF A CARPATHIAN FOOTHILLS - OSADY CZWARORZĘDOWE, Geologia GZMiW UAM 2010-2013,
GF w3 2.03, Geologia GZMiW UAM 2010-2013, I rok, Geologia fizyczna, Geologia fizyczna - wykłady, 01,

więcej podobnych podstron