Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90...

19
Stima dei parametri di anisotropia tipo VTI mediante la tomografia dei tempi d’arrivo GNGTS - 37° Convegno Nazionale Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - Trieste Gualtiero Böhm 19 – 21 novembre 2018 Bologna

Transcript of Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90...

Page 1: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Stima dei parametri di anisotropia tipo VTImediante la tomografia dei tempi d’arrivo

GNGTS - 37° Convegno Nazionale

Istituto Nazionale di Oceanografiae di Geofisica Sperimentale - Trieste

Gualtiero Böhm

19 – 21 novembre 2018Bologna

Page 2: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Sfruttare i segmenti dei raggi (ed i loro angoli)definiti nel ray-tracing della tomografia dei tempi d’arrivo

L’idea

Page 3: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

La teoria

( )θεθθδθ 422o 1VV sincossin)( ++~

θ è l’angolo con la verticaleVo è la componente verticale della velocitàε e δ sono I parametri di Thomsen

Se conosciamo V0 e V90 …

Caso VTI (con anisotropia debole)

Dall’approssimazione della matrice di elasticità:

( ) 0090 VVV /−=ε

?=δ

Page 4: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

La teoria

( )θεθθδθ 422o 1VV sincossin)( ++~

θ è l’angolo con la verticaleVo è la componente verticale della velocitàε e δ sono I parametri di Thomsen

Se conosciamo V0 e V90 …

Dall’approssimazione della matrice di elasticità:

( ) 0090 VVV /−=ε

?=δ )(cos)(sin))(sin()(

θθθεθδ 22

0

40

V1VV +−

=

Caso VTI (con anisotropia debole)

Page 5: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Il metodo

Se non conosciamo V0 and V90 … 1. Modello iniziale (V0, ε e δ)2. Stima di V90 dai raggi orizzontali

3. Ray tracing anisotropo usando il modello attuale

4. Aggiornamento della velocità con l’algoritmo di inversione (SIRT)

5. Calcolo di dV0 per ciascun segmento di raggio

6. Ottengo V0 dalla media dei dV

7. Stima di ε

8. Aggiornamento di ε e V0

9. Ritorno al passo 3

)(ϑdV

( )θεsinθθcosδsin1dV

dV 4220 ++=

)(ϑ

∑=

+=n

1i

ioinio dVVV

n1

( ) 0090 VVV /−=ε

LOO

P

Page 6: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Se non conosciamo V0 and V90 …

Il metodo

Page 7: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Se non conosciamo V0 and V90 … 1. Modello iniziale (V0, ε e δ)2. Stima di V90 dai raggi orizzontali

3. Ray tracing anisotropo usando il modello attuale

4. Aggiornamento della velocità con l’algoritmo di inversione (SIRT)

5. Calcolo di dV0 per ciascun segmento di raggio

6. Ottengo V0 dalla media dei dV

7. Stima di ε

8. Aggiornamento di ε e V0

9. Ritorno al passo 3

)(ϑdV

( )θεsinθθcosδsin1dV

dV 4220 ++=

)(ϑ

∑=

+=n

1i

ioinio dVVV

n1

( ) 0090 VVV /−=ε

LOO

P

Il metodo

Page 8: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Esperimento cross-hole

Page 9: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Esperimento cross-hole

Page 10: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Esperimento cross-hole

Page 11: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Esperimento cross-hole

Page 12: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Esperimento cross-hole

Page 13: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Esperimento cross-hole

Page 14: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Esperimento cross-hole

Page 15: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Tomografia a riflessione

Risultati finali dalla tomografia:

Vvert=3.55 km/s Vor=4.34 km/s rms time residuals (arrivi riflessi) = 0.1 ms (0.03%) epsilon=0.235 H2=-0.303 km

Arrivi riflessi e rifratti

Page 16: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Arrivi riflessi e rifrattiTTI

Tomografia a riflessione

Page 17: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Tomografia da VSP (walk-away)

Arrivi trasmessi e riflessi

TTIVTI

Page 18: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Conclusioni

Metodo per determinare i parametri di anisotropia (VTI)usando la procedura iterativa della tomografia dei tempi d’arrivo.

Il metodo è stato testato con due esempi sintetici di tomografiain trasmissione e tomografia a riflessione/rifrazione.

Lavoro futuro:Analisi della sensitivitàIncludere la stima di δ nella proceduraApplicare il metodo a dati reali;

Problemi/difficoltà: Valore fisso di δ

Variazioni laterali di velocità (dati cross-hole)

Interfacce complesse (tomografia a riflessione)

Page 19: Stima dei parametri di anisotropia tipo VTI mediante la ...Il metodo Se non conosciamo V 0 and V 90 … 1. Modello iniziale (V 0, ε e δ) 2. Stima di V90 dai raggi orizzontali 3.

Ringraziamenti

Grazie per l’attenzione!


C×R type SSOP package 60 V, 80 V and 100 V ... - Panasonic

C×R type SSOP package 60 V, 80 V and 100 V ... - Panasonic

V(r) r rmrm εmεm r=σ Regione attrattiva Regione repulsiva V(r m )=-ε, F attr =F rep V(σ)=0, V attr =V rep.

V(r) r rmrm εmεm r=σ Regione attrattiva Regione repulsiva V(r m )=-ε, F attr =F rep V(σ)=0, V attr =V rep.

INA12x Precision, Low-Power Instrumentation · PDF filepurpose instrumentation amplifiers offering excellent ... (V+) - 2 (V+) - 1.4 V (V) + 2 ... Slew rate VO = ±10 V,

INA12x Precision, Low-Power Instrumentation · PDF filepurpose instrumentation amplifiers offering excellent ... (V+) - 2 (V+) - 1.4 V (V) + 2 ... Slew rate VO = ±10 V,

Gebruiksaanwijzing FAKRO FTS - FTS-V - FTT - FTP-V - FTU-V

Gebruiksaanwijzing FAKRO FTS - FTS-V - FTT - FTP-V - FTU-V

Ω V V Ω V V Ω V V Ω V V Resistance range 0.01 Ω V V Ω V V ... · The “C’DIAG TOUCH” continuity tester is powered by two 1.5 V batteries and its innovative design provides

Ω V V Ω V V Ω V V Ω V V Resistance range 0.01 Ω V V Ω V V ... · The “C’DIAG TOUCH” continuity tester is powered by two 1.5 V batteries and its innovative design provides

PIBETA: A new precise measurement of rare pion and muon …V cd V cs V cb V td V ts V tb ... −0.0017(1) In the meantime we: • have reduced region C background, • have refined

PIBETA: A new precise measurement of rare pion and muon …V cd V cs V cb V td V ts V tb ... −0.0017(1) In the meantime we: • have reduced region C background, • have refined

Análisis y Diseño de Edificaciones de Mampostería · 1 2 V 3 V 4 V 5 V u Cortante sísmico V R,piso = V i V i = F AE F R (0.5v m*A T + 0.3P) i v m* = esfuerzo cortante de diseño

Análisis y Diseño de Edificaciones de Mampostería · 1 2 V 3 V 4 V 5 V u Cortante sísmico V R,piso = V i V i = F AE F R (0.5v m*A T + 0.3P) i v m* = esfuerzo cortante de diseño

Chapitre Les solides équivalents 33 V petit cône ≈ 154,28 m3 On détermine le volume de l’aquarium : V aquarium = V grand cône – V petit cône V aquarium ≈ 509,2 – 154,28

Chapitre Les solides équivalents 33 V petit cône ≈ 154,28 m3 On détermine le volume de l’aquarium : V aquarium = V grand cône – V petit cône V aquarium ≈ 509,2 – 154,28

VSVS L C R At every instant, the generator / supply voltage: (V S = V m sinωt) is given by: V S = V L + V C + V R IZ = IX L + IX C + IR These relationships.

VSVS L C R At every instant, the generator / supply voltage: (V S = V m sinωt) is given by: V S = V L + V C + V R IZ = IX L + IX C + IR These relationships.

Teoria degli Errori e Fondamenti di Statistica11.4.2 Stima a posteriori degli errori di misura . . . . . . . . 183 ... senza il quale questo libro non sarebbe probabilmente mai venuto

Teoria degli Errori e Fondamenti di Statistica11.4.2 Stima a posteriori degli errori di misura . . . . . . . . 183 ... senza il quale questo libro non sarebbe probabilmente mai venuto

Ç o v^ ] } v · î ô &ODVV

Ç o v^ ] } v · î ô &ODVV

CHARAKTERYSTYKI STAŁOPRĄDOWE … · dsp =β p V in −V DD −V tp] 2 [( ) 2 1 2 out dsn n in tn out V I =βV −V V ...

CHARAKTERYSTYKI STAŁOPRĄDOWE … · dsp =β p V in −V DD −V tp] 2 [( ) 2 1 2 out dsn n in tn out V I =βV −V V ...

Review of Resonance - rfic.eecs.berkeley.edurfic.eecs.berkeley.edu/142/pdf/module5.pdf · Series Resonance v R v C v L v s!

Review of Resonance - rfic.eecs.berkeley.edurfic.eecs.berkeley.edu/142/pdf/module5.pdf · Series Resonance v R v C v L v s!

Features Description - Analog, Discrete, Logic, Mixed-Signal · CC Operating Voltage - 1.65 - 5.5 V V IN Control Input Voltage - 0 - V CC V V S Switch Input Voltage - 0 - V CC V V

Features Description - Analog, Discrete, Logic, Mixed-Signal · CC Operating Voltage - 1.65 - 5.5 V V IN Control Input Voltage - 0 - V CC V V S Switch Input Voltage - 0 - V CC V V

Laser LAX 300 - produktinfo.conrad.com · Δ 3 - Δ 2 < ± 2mm Nel calcolo bisogna tenere conto del segno iniziale ! it: H3 H1

Laser LAX 300 - produktinfo.conrad.com · Δ 3 - Δ 2 < ± 2mm Nel calcolo bisogna tenere conto del segno iniziale ! it: H3 H1

Lavoro di un campo elettrico uniforme Consideriamo una carica q che si sposta da una posizione iniziale A ad una posizione finale B in una regione in cui.

Lavoro di un campo elettrico uniforme Consideriamo una carica q che si sposta da una posizione iniziale A ad una posizione finale B in una regione in cui.

Tema 8: Amplificación. leyes de Kirchoff que permiten expresar cualquiera 2 variables del transistor en función de otras dos: 2 1 2 1 1 2, ,, , i i i i v g i v v v v i v f i v C

Tema 8: Amplificación. leyes de Kirchoff que permiten expresar cualquiera 2 variables del transistor en función de otras dos: 2 1 2 1 1 2, ,, , i i i i v g i v v v v i v f i v C

(Progetto Lingua Italiana Dante Alighieri) · PDF fileitaliano iniziale Επίπεδο Α2 ... p.l.i.d.a. b1/plida juniores b1 livello soglia italiano fondamentale Επίπεδο Β2

(Progetto Lingua Italiana Dante Alighieri) · PDF fileitaliano iniziale Επίπεδο Α2 ... p.l.i.d.a. b1/plida juniores b1 livello soglia italiano fondamentale Επίπεδο Β2

Gebruiksaanwijzing FYP-V FYU-V tuimeldakraam

Gebruiksaanwijzing FYP-V FYU-V tuimeldakraam

AlvaZnøm Jnem-^-Øns‚ kp`{ZX°v th≠n XymK-߃°v Xøm-dm-Ip-hn≥ · v v I

AlvaZnøm Jnem-^-Øns‚ kp`{ZX°v th≠n XymK-߃°v Xøm-dm-Ip-hn≥ · v v I