REFLEKSIJSKA SEIZMIKA
REFLEKSIJSKA SEIZMIKA
REFLEKSIJSKA SEIZMIKA
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).