05.05.2013 Views

REFLEKSIJSKA SEIZMIKA

REFLEKSIJSKA SEIZMIKA

REFLEKSIJSKA SEIZMIKA

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

1. Uvod<br />

Refleksijska seizmika<br />

Refleksijska seizmika je glavna geofizička metoda u istraživanju ugljikovodika. Metoda registrira<br />

djelovanje s niza diskontinuiteta u podzemlju (za razliku od ostalih metoda koje sumiraju djelovanje raznih<br />

diskontinuiteta u jednu anomaliju) i stoga daje najpreciznije podatke o strukturama u podzemlju. Seizmička<br />

slika je odraz petrofizikalnih svojstava stijena zbog čijih se razlika i stvaraju refleksi. Ako su dva medija s<br />

jednakim brzinama širenja seizmičkog vala, tada neće biti niti refleksije, niti refrakcije. Mediji mogu biti<br />

litološki različiti, ali sa istim prolaznim brzinama, pa se neće seizmički razlikovati. Isto tako, dva medija<br />

mogu biti generalno jednaki, ali vrlo mala varijacija u brzinama (zbog poroziteta, fluida, pukotina,...) može<br />

dati refleks.<br />

2. Snimanje seizmičkih podataka<br />

Snimanje seizmičkih podataka obavlja se po profilima (2D ili 3D). Kod 2D prikaza, izvor energije i<br />

prijemnici nalaze se na liniji snimanja. Reflektirani signal dolazi od reflektora do geofona i snima se u<br />

digitalnom obliku.<br />

Mjesto izvora seizmičke energije naziva se točka paljenja (SP-shot point), a mjesto prijema signala su<br />

geofoni. Sam odziv odnosi se na CDP (common depth point – zajednička dubinska točka) ili CMP<br />

(common midd point – zajednička središnja točka). CMP uvijek predstavlja točku na polovici razmaka<br />

između točke paljenja i točke prijema. CDP je zajednička točka u dubini na mjestu refleksije. Kako položaj<br />

CDP, zbog zakona refleksije i kosine reflektora, nije uvijek na polovici između dvije točke paljenja, u<br />

obradi se koristi CMP kao zajednička točka na površini. Samo u slučaju horizontalnog reflektora CDP i<br />

CMP imaju potpuno jednak položaj.<br />

Slika1. Odnos CDP-CMP<br />

SP CMP G<br />

CDP<br />

CDP<br />

CMP<br />

Ranije su u refleksijskim istraživanjima mjerenja izvođena metodom kontinuiranog prekrivanja granice.<br />

Unaprijeđenjem tehnologije omogućen je razvoj i uporaba metode višestrukog prekrivanja (CDP).<br />

Metodom višestrukog prekrivanja se postiže da se signal od jedne točke reflektora dobije nekoliko puta<br />

ovisno o broju prekrivanja kako bi se zbrajanjem dobio što čišći signal sa što manje smetnji. Prekrivanje<br />

G


može biti 12, 24, 48, 64 puta i više. Geofonski raspored se sastoji od puno grupa geofona na duljini od više<br />

kilometara. Točka paljenja se postavi na odgovarajućoj udaljenosti od prve grupe geofona. Izvode se<br />

uzastopna paljenja. Nakon prvog paljenja iskopča se prva grupa geofona (grupa uz točku paljenja), a spoji<br />

grupa na drugom kraju. Upotrijebljenu grupu geofona prebaci se na novi položaj na početku profila, izvodi<br />

se novo paljenje i cijeli postupak se ponavlja. Broj prekrivanja (broj tragova od iste dubinske točke)<br />

1 n ⋅d<br />

određen je geometrijom sustava mjerenja i iznosi f = ⋅<br />

gdje je f = broj prekrivanja, d = udaljenost geofona ili grupe geofona, s = udaljenost točaka paljenja i n=<br />

broj kanala ili aktivnih geofona.<br />

TP 1<br />

Slika 2. Metoda višestrukog prekrivanja<br />

3. Obrada podataka<br />

Podaci refleksijskih mjerenja se moraju dodatno obraditi. Obrada podataka obuhvaća tri osnovne vrste<br />

korekcija: statičke, dinamičke i preostale korekcije.<br />

1. Statičke korekcije<br />

TP 2<br />

TP 3<br />

CDP<br />

Provode se zbog:<br />

- varijacije u nadmorskoj visini duž profila. Seizmički profili su postavljeni na neravnom terenu tj. svi<br />

geofoni se ne nalaze na istoj nadmorskoj visini. (Do geofona smještenih u dolini reflektirani val stigne<br />

prije nego do geofona na uzvišenju).<br />

- varijacije sloja niske brzine (low-velocity layer LVL, kora trošenja, trošna zona). Pošto je taj sloj<br />

anizotropno sredstvo vrlo promjenjivih svojstava, poželjno je točku paljenja postaviti ispod sloja niske<br />

brzine.<br />

Korekcija uzima u obzir dubinu trošne zone i promjene brzine u njoj, koje se određuju refrakcijskim<br />

mjerenjima.<br />

- korekcijom se sva mjerenja svode na odgovarajuću horizontalnu ravninu → referentna ravnina (nivo<br />

svođenja; “datum plane”) – slika 4.<br />

Budući da je iznos ovih korekcija jednak za cijeli trag zovu se statičke korekcije.<br />

G 1<br />

2<br />

G 2<br />

s<br />

G 3


površina<br />

Slika 4. Primjena statičkih korekcija<br />

2. Dinamičke korekcije<br />

Provode se zbog povećanja vremena putovanja vala duž geofonskog rasporeda tj. povećanja udaljenosti<br />

geofona i točke paljenja.<br />

T 0<br />

TP<br />

Slika 5. Primjena dinamičke korekcije<br />

grupa geofona varijacije u<br />

nadm. vis.<br />

točka paljenja smještena<br />

ispod LVL<br />

G 1<br />

G 1 G 2<br />

G 2<br />

baza<br />

LVL<br />

nivo svođenja<br />

G 3<br />

ΔT<br />

G 3


Dinamička korekcija se naziva još i NMO (Normal Move Out).<br />

Zbog varijacija u offsetu (udaljenosti točke paljenja od geofona), dolazi do porasta vremena nailazaka.<br />

NMO korigira svaki energetski put vala prema prema onom putu koji je okomit (normalan) na površinu<br />

reflektiranja. Sva vremena nailazaka ( TX ) svode se na nulto vrijeme ( T0 ) tj. najkraće vrijeme putovanja<br />

vala od točke paljenja do reflektora i natrag – okomiti put zrake vala.<br />

Ako nam je poznata brzina vala do reflektora, možemo izračunati dinamičku korekciju za svaki offset.<br />

Nakon toga možemo pomaknuti refleks za vrijednost izračunate korekcije kako bi njegov nailazak<br />

odgovarao vremenu potrebnom da refleks stigne do geofona na nultom offsetu. Na taj način vrijeme<br />

nailaska korigirano je za dodatno vrijeme (∆T) potrebno valu da prijeđe kosi put do udaljenog geofona.<br />

(slika 5.)<br />

Slika 5: Proračun NMO korekcije<br />

gdje je: x – udaljenost TP i geofona,V – brzina u sloju iznad granice refleksije (pretpostavlja se da je<br />

konst.), T0 – vrijeme za okomiti put zrake<br />

2h<br />

T0<br />

V<br />

= TX = T0<br />

+ ΔT<br />

2<br />

2 2<br />

2<br />

2 2 2 2<br />

→ ( ) = ( VT ) + x → V ( T + ΔT<br />

) = V T + X<br />

VT X<br />

2<br />

2 2<br />

Vrijeme putovanja reflektiranog vala je: TX = ( T0<br />

+ ΔT<br />

) = T0<br />

+ 2<br />

0<br />

2<br />

2 x<br />

Povećanje vremena putovanja vala duž geofonskog rasporeda je: ΔT<br />

= T0<br />

+ −T<br />

2 0<br />

V<br />

Brzina reflektiranog vala<br />

2<br />

2 2 2<br />

ili: V ( T ΔT<br />

) = V T + X<br />

→<br />

TP<br />

VT<br />

0<br />

2<br />

x G<br />

0<br />

x<br />

V<br />

2<br />

2 2<br />

V ( TX<br />

2 2<br />

− T0<br />

) = X → V 2 2<br />

TX<br />

T0<br />

−<br />

=<br />

2<br />

2<br />

2<br />

0 + → 0<br />

( 0 + 2T0ΔT<br />

+ ΔT<br />

)<br />

X<br />

2 X<br />

T = T0<br />

+<br />

V<br />

2<br />

2<br />

2 2<br />

V ( 2T0Δ<br />

T + ( ΔT<br />

) ) = X → V =<br />

2<br />

2T0ΔT + ( ΔT<br />

)<br />

x<br />

VT X<br />

2<br />

NMO korekcija ovisi o vremenu putovanja vala za svaki pojedini refleks ( tj. ovisi o brzini) i zato se naziva<br />

dinamička korekcija.<br />

Za provođenje dinamičke korekcije potrebno je poznavati:<br />

2<br />

2<br />

2<br />

0


Vrijeme putovanja vala ( T ),<br />

Udaljenost geofona od točke paljenja ( x ),<br />

Seizmičku brzinu u sloju iznad granice refleksije ( V ).<br />

NMO korekcijom vremena nailazaka su izjednačena i refleksi se svrstavaju u jednu horizontalnu liniju, za<br />

razliku od drugih nailazaka tj. šumova koji zadržavaju oblik krivulje. Tako korigirani refleksi se sumiraju.<br />

Sumiranje tragova (''stekiranje'') ujedno je i kontrola NMO korekcije. Ako je sumirano n tragova,<br />

amplituda refleksa mora biti povećana n puta. Ako je dobro određena brzina V i dobiven dobar ∆T, onda<br />

amplituda mora biti maksimalna. Povećanje amplitude je: = ∑<br />

i=<br />

t A(<br />

)<br />

Slika 6. Stekiranje<br />

Obrađeni seizmički podaci se obično prikazuju vremenskim seizmičkim profilom. Na vodoravnoj osi<br />

unešeni su položaji točaka na profilu, to jest udaljenosti u metrima, a na okomitoj osi unešena su vremena<br />

nailazaka odbijenih valova svedena na vrijednost kada točka paljenja i geofon padaju u istu točku, tj.<br />

vrijeme za okomiti put zrake. Slika podzemlja u vremenskom mjerilu vrlo je slična geološkom presjeku<br />

terena.<br />

Podaci se preračunavaju u dubine. Proračun seizmičkih brzina najvažniji je za određivanje dubine i<br />

geološku interpretaciju – određivanje strukturne građe i litološkog sastava stijena.<br />

Seizmička brzina se dobiva:<br />

- mjerenjem geofonom u bušotini<br />

- zvučnom karotažom<br />

- bušotinskom seizmikom<br />

- iz refrakcijske seizmike<br />

- iz analize brzina u obradi refleksijske seizmike<br />

n<br />

1<br />

a<br />

i t (<br />

)


ZADATAK<br />

(ponijeti paus A4 i milimetarski papir A4)<br />

1. Na paus papir (A4) precrtati sa dobivenog seizmograma:<br />

• prvi nailazak (N1),<br />

• tri najbolje izražena refleksa (R1, R2, R3).<br />

Precrtavati nailaske između 1. i 24. geofona ili 25. i 48. geofona, i to za svaki drugi geofon.<br />

Raspored duž profila:<br />

1 2 3 24 TP 25 48 TP<br />

795m 0m<br />

Udaljenost između susjednih geofona je Δx = 30 m, a udaljenost između zadnjeg geofona i TP je 105<br />

m.<br />

Ukupna duljina od najudaljenijeg geofona do točke paljenja je x = 795 m.<br />

2. Izračunati brzine širenja valova za pojedine reflekse:<br />

Nailasci ΔT (s) V (m/s)<br />

N1<br />

R1<br />

R2<br />

R3<br />

T Tx<br />

T0<br />

− = Δ<br />

T0, Tx – očitati iz seizmograma<br />

(T0 je jednak za sve geofone, a Tx očitati za posljednji geofon).<br />

Brzinu izračunati prema formuli:<br />

2<br />

x<br />

V = <br />

T T<br />

2<br />

x −<br />

2<br />

0<br />

3. Zakon brzina:<br />

V<br />

=<br />

x<br />

2 ( ΔT<br />

) + 2T<br />

ΔT<br />

Na apscisu nanijeti brzine za pojedine reflekse (izračunate u drugom dijelu zadatka), a na ordinatu<br />

pripadna očitana vremena (T0). Graf napraviti na milimetarskom papiru.<br />

2<br />

0


2T(s)<br />

N1<br />

R2<br />

R1<br />

V (m/s)<br />

R3<br />

4. Proračun NMO-a:<br />

• NMO korekciju napraviti za jedan refleks ( R1, R2 ili R3 ).<br />

Br. geofona x (m) Tx (s) ΔT (s) NMO(s)<br />

25 (26)<br />

27 (28)<br />

29 (30)<br />

.<br />

.<br />

.<br />

.<br />

.<br />

.<br />

47 (48)<br />

795 (765)<br />

735 (705)<br />

675 (645)<br />

.<br />

.<br />

.<br />

.<br />

.<br />

.<br />

.<br />

135 (105)<br />

(Ako se radi između 1. i 24. geofona položaji x su isti)<br />

Tx – očitati iz seizmograma za svaki precrtani refleks<br />

T0 – očitani iz sezmograma, jednak je za sve geofone<br />

V – brzina izračunata u prvom dijelu zadatka.<br />

2<br />

2 x<br />

ΔT<br />

= T0<br />

+ −T<br />

2<br />

V<br />

NMO = Tx – ΔT<br />

0<br />

• Vremena dobivena za NMO ucrtati na paus iz prvog dijela zadatka (ako je NMO korekcija<br />

dobro napravljena vrijednosti moraju biti približno jednake T0).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!