05.05.2013 Views

REFLEKSIJSKA SEIZMIKA

REFLEKSIJSKA SEIZMIKA

REFLEKSIJSKA SEIZMIKA

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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!