CONVENZIONE ENEL-INGV Possibile modello di velocità da dati sismici a riflessione Unità di...

CONVENZIONE ENEL-INGVCONVENZIONE ENEL-INGV

Possibile modello di velocità da dati sismici a Possibile modello di velocità da dati sismici a riflessioneriflessione

Unità di Business Produzione Geotermica

Ingegneria Mineraria

Arezzo, 11 Febbraio 2005

2 / 10

Uso: Riservato aziendaleMappa delle linee sismiche e delle postazioni rilevamento sismico

3 / 10

Uso: Riservato aziendale

I dati sismici forniscono una misura indiretta della velocità di propagazione delle onde sismiche. Le analisi di velocità vengono effettuate in dominio Common Depth Point (in media ogni 20 CDP)

Da queste analisi si ottengono:

• velocità di stack: è una velocità quadratica media utilizzata per le correzioni di NMO

• velocità intervallare: velocità media in un intervallo tra due riflettori

Analisi di velocità su dati sismici

4 / 10


LINE LAR-44 SP 461 CDP 780

X 1663656 Y 4779415 T 33 V 2293 T 200 V 2974 T 433 V 3854 T 574 V 4206 T 1359 V 5163 T 2078 V 5536 T 4493 V 5816 T 5903 V 6985

TT1 1 – V– V11

TT2 2 – V– V22

V V intint

1

12

12

int

ii

iiiii tt

tVtVV

Calcolo delle velocità intervallari Dall’analisi di velocità, si ottengono coppie tempi/velocità che caratterizzano le riflessioni identificate (TWT dell’evento a offset zero e velocità di stack, espressi rispettivamente in ms e m/s)

La formula di Dix permette di trasformare le velocità di stack in velocità intervallari caratteristiche dell’intervallo compreso tra due riflettori di individuati .

5 / 10


LINE LAR-44 SP 461 CDP 780

X 1663656 Y 4779415 T 33 V 2293 T 200 V 2974 T 433 V 3854 T 574 V 4206 T 1359 V 5163 T 2078 V 5536 T 4493 V 5816 T 5903 V 6985

Dallo spessore in tempi doppi e dalla velocità intervallare si è potuti risalire allo spessore in metri dei singoli intervalli

Cumulando gli spessori calcolati è possibile ricavare l’andamento delle velocità intervallari in funzione delle profondità

( ) 3091=33200

33×2293200×2974=

22

int -

- V m/

s

Spessore

singolo intervall

o

( )258=3091×

200033200

=- m

Dalle Vint al Calcolo delle distanze tra due riflettori

6 / 10


0

200

400

600

800

1000

1200

1400

1600

1800

2000

2200

2400

2600

2800

3000

3200

3400

0 1000

2000

3000

4000

5000

6000

7000

Velocità (m/s)

Prof. v

ert G

.L/fl

oat

ing (m

)

CDP 260

TWTV.

STACK.V.

INT.PROF

.

0 0 0

2293

33 2293 38

3091

200 2974 296

4473

433 3854 817

5138

574 4206 1179

5834

807 4734 1859

5733

1359 5163 3441

Ricostruzione della funzione “Profondità – Velocità”

7 / 10


• Identificati i pozzi con VSP in prossimità di linee sismiche e selezionati

quelli in cui il rilievo ha raggiunto la profondità maggiore

• Estratte le analisi di velocità relative ai CDP più vicini ai pozzi (le distanze

variano tra un minimo di 130 ad un massimo di 412 m)

• Con il metodo precedentemente descritto, sono stati ricavati per questi

CDP le velocità in funzione delle profondità e confrontate con le velocità

intervallari e formazionali ottenute dai dati di pozzo

Valutazione dell’attendibilità del metodo

8 / 10


Pozzo BADIA_1A – Linea sismica 2D LAR_9

0

200

400

600

800

1000

1200

1400

1600

1800

2000

2200

2400

2600

2800

3000

0 1000

2000

3000

4000

5000

6000

7000

8000

Velocità (m/s)Pr

of. v

ert. G

.L./floa

ting

(m)

DATI VSP

VEL. FORM. DA VSP

CDP 260

CDP 280

FLYSCH

SCAGLIE

TETTONICHE

FILLADI

MICASCISTI

GNEISS

Distanza CDP 260 dal pozzo: 328 m


Confronto delle funzioni di velocità ottenute (1)

9 / 10


Pozzo MONTIERI_4 – Linea sismica 2D LAR_44

0

500

1000

1500

2000

2500

3000

3500

4000

0 1000

2000

3000

4000

5000

6000

7000

8000

9000

Velocità (m/s)

Prof. v

ert. G

.L./float

ing (m

)

DATI VSP FILTRATI (5)

VEL FORM DA VSP

CDP_760

CDP 780

FLYSCH

SCAGLIE

TETTONICHE

FILLADI

ROCCETERMOMET.

ROCCE

INTRUSIVE




10 / 10


Pozzo CHIUSDINO_4 – Linea sismica 2D LAR_45

0

500

1000

1500

2000

2500

3000

3500

4000

0 1000

2000

3000

4000

5000

6000

7000

8000

9000

1000

0

1100

0

1200

0

1300

0

Velocità (m/s)

Prof

. ver

t. G

.L./

floa

ting

(m)

DATI VSP

VEL. FORM. DA VSP

CDP 360

CDP 380

SCAGLIE

TETTONICHE

FILLADI

ROCCETERMOMET.

ROCCE

INTRUSIVE




11 / 10


Pozzo PALAZZACCIO_2 – Linea sismica 2D LAR_36

0

500

1000

1500

2000

2500

3000

3500

4000

0 1000

2000

3000

4000

5000

6000

7000

8000

9000

1000

0

1100

0

Velocità (m/s)

Prof. v

ert. G

.L./float

ing (m

)

DATI VSPVEL. FORM. DA VSPCDP 1160CDP 1180

FLYSCH

ANIDRITI

BURANO

FILLADI

MICASCISTI

CarotaGNEISS

CarotaGNEISS

PTC

PTC




12 / 10


Commento confronto dei dati

1. I rapporti giaciturali e l’assetto geologico strutturale dell’area di Larderello-Travale sono caratterizzati da forti variazioni laterali che si verificano principalmente nelle unità sedimentarie (depositi marino-continentali del Neogene, Unità in facies di Flysch del Cretaceo, Unità a Serie Toscana di età mesozoica e dal complesso a scaglie tettoniche del Trias.

2. I pozzi con VSP di cui è disponibile la funzione velocità non coincidono come ubicazione con i CDP analizzati, ma le distanze variano da 130 a 412 m.

Per quanto sopra possiamo pensare che un modello di velocità ricavato con i dati disponibili di sismica a riflessione possa avere una buona attendibilità.

Buona corrispondenza tra i profili di velocità ricostruiti dai dati di pozzo e quelli ottenuti dai CDP ricavati dai dati di sismica a riflessione.

Le principali differenze sono evidenti soprattutto nelle porzioni superficiali, mentre in profondità le differenze si fanno sempre meno marcate. Questo è dovuto principalmente a due fattori:

A) Neoautoctono

B) Insieme Ligure/Subligure

C) Insieme Toscano

D) Complesso a Scaglie Tettoniche

E) Basamento Metamorfico

F) Termometamorfiche e Intrusive

13 / 10


• Le analisi di velocità sono riferite al floating datum

• Il floating datum è una superficie che si utilizza nel processing pre-NMO e approssima la topografia

• Maggiormente accidentata è la topografia, maggiore è la differenza tra le due curve

• NON E’ DIRETTAMENTE DISPONIBILE PER CIASCUN CDP IL VALORE DELLE QUOTE DEL FLOATING DATUM

• Valutiamo un possibile incertezza della quota di riferimento di ciascuna analisi di velocità non superiore a 300 m

• Questa incertezza può essere eliminata calcolando per ogni CDP la quota del floating datum

Affidabilità

14 / 10


Per calcolare la differenza di quota tra floating e datum di riferimento (+200) occorre :

•Estrarre i valori C_STATIC dalla header di traccia di ciascun CDP (con adeguato software)

•Identificare i valori delle velocità superficiali utilizzate per le correzioni statiche

•Calcolare la quota del floating datum rispetto al piano di riferimento (+ 200 m slm)

Non avendo a disposizione il software adeguato, occorre

rivolgersi a terzi

Calcolo della quota del floating

15 / 10


Dai 1024 CDP possono essere ricavati circa 9500 valori di velocità intervallare riferibili a:• un volume di crosta di circa 16.800 Km3

• su una superficie di 800 Km2

• per una profondità media di 21 Km

In tale volume ricadono il 94% di tutti i terremoti registrati dalla Rete Sismica di LARDERELLO-TRAVALE

Quantità di dati disponibili dalla sismica a riflessione 2DPer tutta l’ara di Larderello-Travale sono disponibili 1024 CDP con analisi di velocità.Ulteriori A.V. potrebbero essere ricavate da recenti rilievi 3D

16 / 10


Bruciano 101 A.V. con griglia di 500 x 500 m

Montieri-Chiusdino 130 A.V. con griglia di 500x800m

Gabbro-Sesta 140 A.V. con griglia 500x800 m

Lagoni-Rossi Val di Cornia ancora da elaborare

Quantità di dati disponibili dalla sismica a riflessione 3DUlteriori Analisi di Velocità potrebbero essere ricavate da rilievi 3D

17 / 10


Scelta del Modello apriori di Velocità

• Valore di velocità omogeneo per tutta l’area

• Valori di velocità dei modelli attualmente utilizzati per le stazioni della rete sismica

• Modello unico per tutte le stazioni (Tipo Velest) già utilizzato come test (Tesi di Righi)

• Dati di velocità del modello 3D (Tomografia Zollo) già utilizzato per la sequenza Bulera

• Modello ricostruibile dalle A.V. da Rilievi sismici a riflessione (1024 CDP)

CONVENZIONE ENEL-INGV Possibile modello di velocità da dati sismici a riflessione Unità di...

Documents

Transcript of CONVENZIONE ENEL-INGV Possibile modello di velocità da dati sismici a riflessione Unità di...