Programma Operativo Nazionale “Ricerca e Competitività 2007-2010”
Sistemi Integrati per il Monitoraggio, l’Early
Warning e la Mitigazione del Rischio Idrogeologico
lungo le grandi vie di comunicazione
OR 4 - MODELLI
25 Luglio 2013
OR 4 – MODELLI
Modelli di previsione degli spostamenti
Previsione (Forecasting): capacità di valutare
probabilisticamente tempo, ubicazione e
dimensione di un prossimo evento disastroso.
OR 4 – MODELLI - Composizione
OR.4
WP 4.1
Modello areale - GEOTOP
WP 4.2
Modello completo di versante - SUSHI
WP 4.3
Modello di propagazione – SCIDDICA d1
Modello Idrologico FLaIR
(Forecasting of Landslides Induced by Rainfall)
OR 4 – MODELLI - Composizione
Il modello FLaIR  operativo
Il modello SCIDDICA  operativo
Il modello GETOP  [I]
Il modello SUSHI  [I]
STRUTTURA SOFTWARE
Pre-Processing
Processing
Modulo
Idraulico
[C]
Post-Processing
Modulo
Geotecnico
[C]
Pre-Processing
Processing
1. Carica Dxf Strati
[C]
Lettura file Dxf, Individuazione e gestione Polilinee
2. Carica Dati Stratigrafia
Post-Processing
GEOMETRIA E
STRATIGRAFIA
[C]
Lettura file testo, Associazione dati
fisico-meccanici alle polilinee
3. Carica Dati Pioggia
[C]
- Lettura e memorizzazione in Struttura dati
4.  polilinea, combinazioni di Macroelementi
[C]
5. Def. combinazione ottimale Macroelementi
[C]
 Criterio di Valutazione Regolarità
6. Da strutt. Polilinea a Strutt. Macroelementi
INPUT DATI
• Rainfall data
• Physical and mechanical param.
• Boundary and initial conditions
CREAZIONE MESH
PRIMARIA
[C]
7. Visualizzazione a schermo della mesh generata [N]
8. Numerazione nodi nel Riferimento Globale
[I]
9. Infittimento mesh
[N]
10. Algoritmo di rinumerazione dei nodi
[N]
11.Gestione Interfaccia Grafica
[I]
INFITTIMENTO
MESH
Pre-Processing
Processing
Post-Processing
Modulo Idraulico
Mag. Lug
2013 2013
1. Implementazione e gestione elemento finito
[C]
[C]
[C]
[C]
3. Formulazione isoparametrica
[C]
[C]
4. Traduzione e passaggio al Linguaggio C++
[N]
[C]
5. Implementazione matrici dell’elemento
[C]
[C]
6. Assemblaggio matrici globali
[C]
[C]
7. Definizione forzante
[C]
[C]
8. Definizione condizioni al contorno
[C]
[C]
9. Discretizzazione temporale
[C]
[C]
10.Procedura iterativa non lineare
[C]
[C]
2. Scrittura Equazione di Richards
in forma matriciale
EXIT
Modulo
Geotecnico
Procedura iterativa
implementata in
MAPLE
CONFRONTO
Procedura iterativa
implementata in
C++
Pre-Processing
Processing
Post-Processing
Modulo
Idraulico
Mag. Lug
2013 2013
…
Calcolo Fattore di sicurezza
Modulo
Geotecnico
[C]
[C]
1. Inizializzazione e allocazione memoria
[C]
[C]
2. Acquisizione modello di simulazione
[C]
[C]
3. Def. Sistema algebrico risolutivo
[C]
[C]
[C]
[C]
5. Calcolo spostamenti e stato tenso-deformativo
[C]
[C]
6. Verifica violazione criterio snervamento
[C]
[C]
[I]
[C]
Risposta tenso-deformativa
4. Generazione stato tensionale iniziale
con procedura gravity loading
7. Analisi step-by-step
in funzione dell’output del modulo idraulico
Pre-Processing
Processing
Post-Processing
MODEL VALIDATION
Software
Output
COMPARISON
Instruments
in Place
DISPLAYING RESULTS
Modello areale - GEOTOP
Fino ad Aprile 2013
GEOtop-OMS3 Integration
Sarno Application
Modello areale - GEOTOP
Work in progress:
Applicazione a nuovi casi
Sviluppo del modulo per il calcolo probabilistico del fattore di
sicurezza
Sviluppo della procedura per la visualizzazione delle mappe FS
Work in progress (1/3): Drake river basin calibration
Input data: Radiation
Rainfall
Air Temperature
Relative Humidity
Simulation time step: daily
Work in progress (2/3): Drake river basin raster sim.
Input data: Radiation
Rainfall
Air Temperature
Relative Humidity
Simulation time step: daily
Soil moisture raster maps at 30 cm depth
Work in progress (3/3): frana A3 Sa-Rc Rogliano
Soil Water Content Simulation
Km 260
25-01-2009 ore
22.00 circa
Simulation period: 01-02-2008 – 01-02-2009
Time step: hourly
Where the landslides are?
3D-(Multi)Point simulations
Measurement
data
MODEL OUTPUT:
Model
Setup
I/O GIS
integration
Parameter
estimation
•
Soil moisture maps
for each layer (SM)
•
Suction maps for
each layer (WP)
Raster simulations
SAFETY FACTOR MAP
COMPUTATION:
Model
Setup
I/O GIS
integration
17
Formetta
G. Formetta,
G., ARS-USDA-Fort
Leipzig
University
Trento 24
05June
of
July
Trento
2011
2012
Collins
(Italy)
(CO)
GEOtop
Model
FS(WP, SM, C, φ,
SLOPE, etc.)
15
Prove e parametri necessari

Definizione della geometria – dominio di indagine

Definizione delle coperture

Caratterizzazione dei parametri idraulici e geotecnici

Definizione delle Curve Caratteristiche (SWRC)
Validazione degli spostamenti calcolati:
• Spostamenti superficiali
• Spostamenti profondi
Sensori superficiali
Misure inclinometriche
Programma Operativo Nazionale “Ricerca e Competitività 2007-2010”
Sistemi Integrati per il Monitoraggio, l’Early
Warning e la Mitigazione del Rischio Idrogeologico
lungo le grandi vie di comunicazione
OR 4 - MODELLI
25 Luglio 2013
Modello elasto-plastico perfetto alla Mohr-Coulomb
I parametri che descrivono il modello sono 6:
• Peso dell’unità di volume γ
• Modulo elastico di Young E
• Coefficiente di Poisson ν
• Coesione c’
• Angolo di attrito interno φ’
• Angolo di dilatanza ψ
Identificazione del
terreno
Prove triassiali,
edometriche
Prove triassiali, taglio
diretto
Validazione
• Valori della pressione neutra e suzione (piezometro e tensiometro)
• Valori del contenuto d’acqua (TDR)
• Spostamenti superficiali
• Spostamenti profondi
• Stato tensionale (celle di pressione)
Sensori superficiali
Misure inclinometriche
Modello elasto-plastico perfetto alla Mohr-Coulomb
I parametri che descrivono il modello sono 6:
• Peso dell’unità di volume γ
Identificazione del terreno
• Modulo elastico di Young E
• Coefficiente di Poisson ν
Prove triassiali, edometriche
• Coesione c’
• Angolo di attrito interno φ’
• Angolo di dilatanza ψ
Prove triassiali, taglio diretto
OR 4 – MODELLI
1. Collegamenti con le altre OR
2. Dati e misure in ingresso ai modelli
3. Uscite dei modelli
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Uscite dei modelli
OR 1
OR 2
OR 6
OR4
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Uscite dei modelli
OR 1
Scelta appropriata del sito di indagine/applicazione dei modelli
OR 2
Scelta dei sensori/strumentazioni a supporto dei modelli
OR 6
Definizione dei flussi di scambio tra il CAED ed i modelli
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Uscite dei modelli
Modello puntuale – SUSHI
• Contenuti Idrici
• Valori delle suzioni
e pressioni neutre
Modulo
Idraulico
Modulo
Geotecnico
• Stati tensionali
• Stati deformativi
• Spostamenti
• LIVELLO 1
• FS Equilibrio Limite (rigido
plastico)
• LIVELLO II
• FS pendio (elasto-plastico)
• Tensioni e Deformazioni
USCITA
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Uscite dei modelli
Modello puntuale – SUSHI
• Contenuti Idrici
• Valori delle suzioni
e pressioni neutre
Modulo
Idraulico
Modulo
Geotecnico
• Stati tensionali
• Stati deformativi
• Spostamenti
• LIVELLO 1
• FS Equilibrio Limite (rigido
plastico)
• LIVELLO II
• FS pendio (elasto-plastico)
• Tensioni e Deformazioni
USCITA
File txt con contenuti alfa-numerici
Immagini a bande di colori (s, e, s, u…)
Immagini con vettori direzione spostamento……
E.W.S (???)
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
E.W.S (???)
Uscite dei modelli
FS pendio
FSpendio
>FSsoglia
<FSsoglia
FSpendio≥FS1
FS1<FSpendio≤FS2
FSpendio>FS2
Procedura 1
Procedura 2
Procedura 3
Nessuna
procedura
E.W.S (???)
Cumulative Rainfal
Collegamenti con le altre OR
150
H15 - Stanley
Uscite dei modelli
N49 - Tuen Mun
Dati e misure in ingresso ai modelli
Risposta tenso-deformativa
100s
50
0
u
19/5/2007 19/5/2007 19/5/2007 19/5/
12:00
15:00
18:00
21:
PEt    PEt | Y  p Y   dY
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Uscite dei modelli
Y(t)
Valutare la Criticità
Y (t )
Trasmissione dei dati di
pioggia IN TEMPO REALE
Y3
Y2
Y1
t t
t
t
Attivazionedella
dellafase
fase132
attivazione
t
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Uscite dei modelli
Modello propagazione – Sciddica
Ad ogni Step e per ogni
Cella è possibile avere le
seguenti informazioni:
• Quota
• Spessore di Detrito in
Transito o Accumulato
• Velocità ed Energia del
Detrito
• Profondità della
Copertura Detritica
Alla fine della
Simulazione è possibile
avere GRID-ASCII
contenete:
• Area invasa dalla Frana
• Morfologia Modificata
• Depositi di Detrito
Finali
• Profondità della
Copertura Detritica
Finale
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Uscite dei modelli
Modello propagazione – Sciddica
Step 0
Collegamenti con le altre OR
Dati e misure in ingresso ai modelli
Modello propagazione – Sciddica
Step finale
Uscite dei modelli
Scarica

Uscite dei modelli