Geoestatistikaren oinarrizko kontzeptuak geofisikan
Geoestatistika estatistikaren adarra da, espazialki erlazionatutako datuak aztertzeko bereziki diseinatua, hau da, kokapenak eragina duen balioa duten datuak. Geofisikan, datuak ia beti dira espazialak: grabitatearen, magnetismoaren, sismikoen eta erresistentziaren neurketak, baita geokimikoak eta putzuen erregistroak ere, koordenatu espezifikoei lotuta daude. Beraz, geoestatistikak oinarri ezinbestekoa eskaintzen du lurpeko parametroen banaketa-ereduak ulertzeko, neurtu gabeko kokapenetako balioak kalkulatzeko eta interpretazioen ziurgabetasuna ebaluatzeko. Artikulu honek geoestatistikaren oinarrizko kontzeptuak aztertzen ditu, testuinguru geofisikoetan gehien erabiltzen direnak.
Zergatik da garrantzitsua geoestatistika geofisikan?
Geofisikako azterketak etengabe aurkitzen dituzte laginketa-mugak. Ezinezkoa da gainazaleko edo lurpeko puntu guztiak neurtzea kostu, denbora eta sarbide-mugak direla eta. Ondorioz, interpretazio geofisikoak interpolazioa eta modelizazioa behar ditu. Alderantzizko distantzia-ponderazioa (IDW) bezalako interpolazio sinpleak komenigarriak dira, baina askotan korrelazio espazialeko egiturak alde batera uzten dituzte eta ez dute ziurgabetasunaren neurririk ematen. Geoestatistikak bi behar nagusi jorratzen ditu: (1) datuen korrelazio espazialeko ereduak ustiatzea estimazio errealistagoak lortzeko, eta (2) ziurgabetasunaren esparru kuantitatibo bat eskaintzea esplorazio-erabakiak informatzeko.
Geofisikako praktikan, geoestatistika anomalia-mapa finduagoak eta koherenteagoak sortzeko, arroken propietateak modelatzeko (porositatea, iragazkortasuna, uhinen abiadura), datu iturri anitzekoak integratzeko (adibidez, sismikoak eta putzuak) eta arrisku-eszenarioetarako lurpeko ereduak simulatzeko erabiltzen da.
Datu espazialak eta ausazko eremuen kontzeptua
Geoestatistikaren funtsa fenomeno geologiko edo geofisikoak funtzio edo eremu aleatorio gisa ikustea da. Horrek esan nahi du dentsitatea, suszeptibilitate magnetikoa edo P abiadura bezalako parametro-balioak kokapen bakoitzean egitura espezifiko bat duen prozesu aleatorio baten gauzatze gisa hartzen direla. Ikuspegi honetan, gure helburua ez da "mapa bat marraztea" soilik, baizik eta balioen banaketa eta kokapenen arteko korrelazioak kalkulatzea.
Bi osagai garrantzitsu daude eredu geoestatistiko batean: joera eta aldakortasun lokala. Joerak eskala handiko aldaketak deskribatzen ditu (adibidez, dentsitatearen handitzea iparraldetik hegoaldera eskualdeko litologia-aldaketak direla eta). Aldakortasun lokalak eskala txikiko gorabeherak deskribatzen ditu, askotan arroken heterogeneotasunarekin, hausturekin edo fazies-aldaketekin lotuta. Joeraren eta aldakortasunaren arteko bereizketak metodo egokia aukeratzen laguntzen digu: baldintza geldikorrak suposatu besterik ez egin, edo jitoa/joera esplizituki sartu.
Gelkortasuna: Oinarrizko hipotesi erabiliena
Geoestatistika metodo klasiko askok geldiunearen hipotesian oinarritzen dira, hau da, datuen propietate estatistikoak aldatu gabe mantentzen dira kokapena aldatzean. Forma ohikoena bigarren mailako geldiunea da: batez bestekoa konstantea da, eta kobariantza bereizketaren distantziaren eta norabidearen (atzerapenaren) araberakoa da soilik, ez posizio absolutuaren araberakoa.
Geofisikan, hipotesi hau ez da beti zuzena, baldintza geologikoak pixkanaka aldatzen baitira askotan. Hala ere, eskala jakin batzuetan (adibidez, litologia-domeinu bakar baten barruan), geldikortasunaren hipotesia nahiko arrazoizkoa da askotan. Datuek joera sendoa erakusten badute, joera kentzeko prozesua egiten da normalean edo joera hori egokitzen duen metodo-aldaera bat erabiltzen da, hala nola kriging unibertsala.
Variograma: Geoestatistikaren bihotza
Geoestatistikan kontzepturik ikonikoena bariograma (edo erdibariograma) da. Bariogramak datuen balioen antzekotasuna distantziarekin nola aldatzen den erakusten du. Intuitiboki, elkarrengandik hurbil dauden bi puntuk balio antzekoak izateko joera dute, eta urrun dauden puntuek, berriz, balio desberdinak izateko joera dute. Bariogramak printzipio hau kuantifikatzen du.
Erdibariograma enpirikoa, oro har, honela kalkulatzen da:
\[
\gamma(h)=\frac{1}{2N(h)}\sum_{i=1}^{N(h)}[Z(x_i)-Z(x_i+h)]^2
\]
non \(h\) atzerapena (distantzia eta norabidea) den, \(N(h)\) atzerapen horretan dauden datu-bikoteen kopurua, eta \(Z(x)\) datuen balioa.
Hiru bariograma-parametro garrantzitsu:
1. Nugget-a: zerora hurbiltzen den atzerapenarekin erdibariogramaren balioa. Nugget-ek neurketa-errorea, zarata edo heterogeneotasuna islatzen dute laginketa-distantzia baino eskala txikiagoan.
2. Sill: erdibariogramaren balioa goi-ordoki batera iristen denean. Hau domeinu horretako datuen bariantza osoarekin erlazionatuta dago.
3. Ibiltartea: bariograma ertzara hurbiltzen den distantzia. Ibiltarte horren azpitik, datuak oraindik korrelazionatuta daude; ibiltarte horren gainetik, korrelazioa ahula/galduta dago.
Bariogramek anisotropia ere erakuts dezakete, hau da, norabidearen arabera aldatzen den korrelazioa. Geofisikan, anisotropia askotan egitura geologikoetatik sortzen da, hala nola geruzatzetik, failatik edo sedimentuen fluxuaren norabidetik. Bariograma norabidedun batek jarraitutasun handiena eta txikiena duten norabideak zehazten laguntzen du, eta hori funtsezkoa da urtegien modelizazioan edo egitura-interpretazioan.
Bariograma Ereduak: Funtzio Enpirikoetatik Matematikoetaraino
Bariograma enpirikoak askotan zaratatsuak dira eta ez dituzte beti kriging-ean erabiltzeko baldintza matematikoak betetzen. Beraz, honelako eredu teoriko batekin egokitu behar dira:
– Esferikoa
– Esponentziala
– Gaussiarra
– Matern (malguagoa, baina konplexuagoa)
Modeloaren hautaketa bariograma enpirikoaren formak eta ulermen geologikoak gidatzen dute. Adibidez, modelo gaussiarrek trantsizio oso leunak sortzen dituzte distantzia laburretan, parametro oso jarraituetarako egokiak. Modelo esponentzialak zakarragoak dira distantzia laburretan, aldaketa azkarrak dituzten fenomenoetarako egokiak.
Kriging: Variograman Oinarritutako Estimazio Optimoa
Kriginga geoestatistikako interpolazio metodo bat da, bariograma bat erabiltzen duena estimatzaile lineal inpartzial onena (BLUE) emateko. Interpolazio metodo deterministak ez bezala, krigingak:
1. Kontuan hartu distantzia eta espazio-korrelazioa.
2. Estimazio-mapa bat eta kriging bariantza-mapa (ziurgabetasun) bat sortzen ditu.
Kriging mota ohikoenak:
– Kriging Sinplea (SK): batez bestekoa ezaguna eta konstantea da.
– Kriging arrunta (OK): batez bestekoa ezezaguna da, baina tokiko auzoan konstantea dela suposatzen da.
– Kriging Unibertsala (Erresuma Batua): joera/noraeza barne hartzen ditu (adibidez, koordenatuekiko funtzio polinomiala).
– Ko-kriginga: bigarren mailako aldagaiak erabiliz (adibidez, porositatea kalkulatzea inpedantzia sismikoaren laguntzaz).
– Adierazleen Kriging-a: datu kategorikoetarako/gertaeren araberakoetarako (adibidez, litologia jakin baten probabilitatea).
Geofisika aplikatuan, kriging arrunta hasierako aukera sendoa da askotan, malgua delako eta ez duelako batez besteko global ezagun baten hipotesirik behar.
Simulazio geoestatistikoa: mapa bat baino gehiago
Interpolazioak "onena" den eredu bakarra sortzen du, baina lurpeko azalera ez da inoiz ziurra. Arriskua ebaluatu eta eszenarioak sortzeko, simulazio geoestatistikoak erabiltzen dira, hala nola:
– Simulazio Gaussiar Sekuentziala (SGS) aldagai jarraituetarako.
– Adierazle Sekuentzialen Simulazioa (SIS) aldagai kategorikoetarako.
– Entrenamendu-irudietan oinarritutako puntu anitzeko estatistikak eredu geologiko konplexuetarako.
Simulazioek datuekin eta bariogramarekin koherenteak diren hainbat gauzatze sortzen dituzte, eta horrek aukera ematen digu probabilitate-tarteak, kuantilak eta probabilitateak kalkulatzeko. Testuinguru geofisiko batean, simulazioak ezinbestekoak dira zulaketak planifikatzeko, bolumen-ziurgabetasuna ebaluatzeko eta fluxu-modelizazioarekin integratzeko.
Modeloen Balidazioa: Gurutzatutako Balidazioa eta Diagnostikoak
Geoestatistika ona ez da bariogramak eta kriging-ak sortzearekin amaitzen. Ebaluazioak egin behar dira, adibidez:
– Bat kanpoan utzitako gurutzadura-balioztatzea: puntu bakoitza beste puntu batzuk erabiliz iragartzen da, eta gero benetako balioarekin alderatzen da.
– Hondarren analisia: hondarrak ausaz banatuta dauden ala ez, alborapen sistematikorik dagoen ala ez.
– Kriging bariantzaren egiaztapena: ziurgabetasuna arrazoizkoa den ala ez (handia datu urriko eremuetan, baxua datu-dentsitate handiko eremuetan).
Balidazioak bariograma "leunegia" den, tartea luzeegia/laburregia den edo jaso ez den anisotropiarik dagoen zehazten laguntzen du.
Geofisika eta Geoestatistika Aplikazioetako Ohiko Erronkak
Maiz sortzen diren erronka batzuk:
1. Zarata eta ez-gaussiarrak: Datu geofisikoek askotan muturreko balioak eta banaketa ez-normalak izaten dituzte. Batzuetan eraldaketak (adibidez, puntuazio normala) beharrezkoak dira.
2. Laginketa irregularra: neurketa-bideak (ibilbideak) datuak norabide batean trinkoak eta bestean urriak izatea eragiten du.
3. Ez-egonkortasuna: litologian edo eskualde-egituran izandako aldaketek joera sendoak sortzen dituzte.
4. Eskala anitzeko integrazioa: sismika eskala handikoa da, baina bereizmena desberdina da putzuko datu oso zehatzekin alderatuta.
Erronka hauei aurre egiteko, ulermen estatistikoaren eta intuizio geologiko/geofisikoaren konbinazioa behar da.
Itxiera
Geoestatistikaren oinarrizko kontzeptuek —egonkortasuna, bariograma, kriging-a eta simulazioa— kokapenean oinarritutako datu geofisikoak prozesatzeko esparru indartsua eskaintzen dute. Bariograma batekin, korrelazio espazialaren egitura mapatzen dugu; kriging-arekin, estimazio optimoak eta haien ziurgabetasunak lortzen ditugu; eta simulazioarekin, lurpeko ziurgabetasunak modu errealistan irudikatzen dituzten hainbat eszenatoki eraikitzen ditugu. Azken finean, geoestatistika ez da mapatzeko teknika bat soilik, baizik eta geofisikako esplorazio eta interpretazio erabaki informatuagoak hartzeko ikuspegi kuantitatibo bat.
Nahi izanez gero, aplikazio kasuen adibideak ere gehi ditzaket (adibidez, anomalia magnetikoen mapaketa edo abiadura sismikoaren estimazioa), edo software espezifiko batean bariogramaren kalkulutik kriging-era arteko lan-fluxu praktiko bat sar dezaket.
