PROGETTO SPOT - SISMICITÀ POTENZIALMENTE INNESCABILE OFFSHORE E TSUNAMI Report integrato di fine progetto Una collaborazione del MiSE, con il supporto tecnico del DPC, con CNR-ISMAR, INGV, EUCENTRE e ReLUIS DGISSEG Antoncecchi I., Ciccone F., Dialuce G., Grandi S., Terlizzese F. Di Bucci D., Dolce M. Argnani A., Mercorella A., Pellegrini C., Rovere M., con la partecipazione dell’Università di Bologna: Armigliato A., Pagnoni, G., Paparo M.A., Tinti, S., Zaniboni F. Basili R., Cavallaro D., Coltelli M., Firetto Carlino M., Lipparini L., Lorito S., Maesano F.E., Romano F., Scarfì L., Tiberti M.M., Volpe M., con la partecipazione dell’Università di Pavia: Fedorik J., Toscani G. Borzi B., Faravelli M., Bozzoni F., Pascale V., Quaroni D., Germagnoli F. Belliazzi S., Del Zoppo M., Di Ludovico M., Lignola G.P., Prota A. I contenuti di questo report sono tutelati dalla licenza Creative Commons 4.0 International (CC BY 4.0). Il permesso di utilizzare tali contenuti viene concesso con il pieno riconoscimento della fonte, utilizzando la seguente citazione raccomandata: Antoncecchi I., Ciccone F., Dialuce G., Grandi S., Terlizzese F., Di Bucci D., Dolce M., Argnani A., Mercorella A., Pellegrini C., Rovere M., Armigliato A., Pagnoni, G., Paparo M.A., Tinti, S., Zaniboni F., Basili R., Cavallaro D., Coltelli M., Firetto Carlino M., Lipparini L., Lorito S., Maesano F.E., Romano F., Scarfì L., Tiberti M.M., Volpe M., Fedorik J., Toscani G., Borzi B., Faravelli M., Bozzoni F., Pascale V., Quaroni D., Germagnoli F., Belliazzi S., Del Zoppo M., Di Ludovico M., Lignola G.P., Prota A. (2020) - Progetto SPOT - Sismicità Potenzialmente innescabile Offshore e Tsunami. Report integrato di fine progetto. DOI: 10.5281/zenodo.3732887 ISBN 9788894366945 31 marzo 2020
45
Embed
Report integrato di fine progetto...PROGETTO SPOT - Report integrato di fine progetto 3 Nome Affiliazione ORCID e-mail Antoncecchi Ilaria 1,9* 0000-0002-3164-0794 [email protected]
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
PROGETTO SPOT - SISMICITÀ POTENZIALMENTE INNESCABILE OFFSHORE E TSUNAMI
Report integrato di fine progetto
Una collaborazione del MiSE, con il supporto tecnico del DPC, con CNR-ISMAR, INGV, EUCENTRE e ReLUIS
DGISSEG
Antoncecchi I., Ciccone F., Dialuce G., Grandi S., Terlizzese F.
Di Bucci D., Dolce M.
Argnani A., Mercorella A., Pellegrini C., Rovere M., con la partecipazione dell’Università di Bologna: Armigliato A., Pagnoni, G., Paparo M.A., Tinti, S., Zaniboni F.
Basili R., Cavallaro D., Coltelli M., Firetto Carlino M., Lipparini L., Lorito S., Maesano F.E., Romano F., Scarfì L., Tiberti M.M., Volpe M., con la partecipazione dell’Università di Pavia: Fedorik J., Toscani G.
Borzi B., Faravelli M., Bozzoni F., Pascale V., Quaroni D., Germagnoli F.
Belliazzi S., Del Zoppo M., Di Ludovico M., Lignola G.P., Prota A.
I contenuti di questo report sono tutelati dalla licenza Creative Commons 4.0 International (CC BY 4.0). Il permesso di utilizzare tali contenuti viene concesso con il pieno riconoscimento della fonte, utilizzando la seguente citazione raccomandata:
Antoncecchi I., Ciccone F., Dialuce G., Grandi S., Terlizzese F., Di Bucci D., Dolce M., Argnani A., Mercorella A., Pellegrini C., Rovere M., Armigliato A., Pagnoni, G., Paparo M.A., Tinti, S., Zaniboni F., Basili R., Cavallaro D., Coltelli M., Firetto Carlino M., Lipparini L., Lorito S., Maesano F.E., Romano F., Scarfì L., Tiberti M.M., Volpe M., Fedorik J., Toscani G., Borzi B., Faravelli M., Bozzoni F., Pascale V., Quaroni D., Germagnoli F., Belliazzi S., Del Zoppo M., Di Ludovico M., Lignola G.P., Prota A. (2020) - Progetto SPOT - Sismicità Potenzialmente innescabile Offshore e Tsunami. Report integrato di fine progetto.
DOI: 10.5281/zenodo.3732887 ISBN 9788894366945
31 marzo 2020
PROGETTO SPOT - Report integrato di fine progetto
2
Ministero dello sviluppo economico Direzione generale per la sicurezza anche ambientale delle attività minerarie ed energetiche
Ufficio nazionale minerario per gli idrocarburi e le georisorse Palazzo Piacentini - Via Veneto, 33 - 00187 Roma
1. Ministero dello Sviluppo Economico – DGISSEG 2. Presidenza del Consiglio dei Ministri – Dipartimento della Protezione Civile 3. Consiglio Nazionale delle Ricerche – Istituto di Scienze Marine 4. Università di Bologna – Dipartimento di Fisica e Astronomia 5. Istituto Nazionale di Geofisica e Vulcanologia 6. Università di Pavia – Dipartimento di Scienze della Terra e dell’Ambiente 7. EUCENTRE Centro Europeo di Formazione e Ricerca in Ingegneria Sismica 8. Rete dei Laboratori Universitari di Ingegneria Sismica 9. Ricerca sul Sistema Energetico S.p.A.
* Affiliazione secondaria
Attività svolta nell'ambito del programma CLYPEA - Innovation Network for Future Energy, finanziato dal Ministero dello Sviluppo Economico su fondi per la "sicurezza offshore" ai sensi dell'art.35 D.L. 83/2012
PROGETTO SPOT - Report integrato di fine progetto
4
Introduzione
Il progetto SPOT (Sismicità Potenzialmente innescabile Offshore e Tsunami) è stato sviluppato con lo scopo di supportare le Autorità italiane nell'applicazione della Direttiva Europea sulla sicurezza delle operazioni in mare nel settore degli idrocarburi (2013/30/EU), su fondi di cui art. 35 del Decreto Legge 83/2012, e dei decreti italiani che ne derivano (Antoncecchi et al., 2019). Il progetto, della durata di 21 mesi, è stato ideato e finanziato dal Ministero dello Sviluppo Economico italiano – Direzione Generale per la sicurezza delle attività minerarie ed energetiche (DGS-UNMIG) nell’ambito della rete di ricerca CLYPEA, con il supporto tecnico del Dipartimento della Protezione Civile nazionale.
Finalità
Il progetto ha avuto come obiettivo generale l’implementazione di alcuni degli approfondimenti suggeriti dalla Commissione Grandi Rischi dopo il terremoto dell’Emilia del 2012. È stata condotta un’ampia analisi delle strutture geologiche presenti nell’offshore delle coste italiane, volta a valutare l’esistenza di faglie potenzialmente sismogeniche in prossimità delle piattaforme di estrazione olio e gas in produzione, e di corpi sedimentari potenzialmente instabili dal punto di vista gravitativo e suscettibili di (ri-) attivazione a seguito di terremoti generati dalle suddette faglie. A questa analisi ha fatto seguito una valutazione generale, alla scala nazionale, del potenziale impatto e delle conseguenti perdite causate da terremoti e tsunami generati dalle strutture riconosciute, sia tettoniche sia gravitative. Il progetto, dunque, riguarda alcune pericolosità naturali, e deve essere considerato come una fase propedeutica a un’eventuale valutazione, successiva, della sismicità potenzialmente innescata connessa a operazioni svolte sulle piattaforme.
Enti di ricerca coinvolti
Sono stati coinvolti i seguenti quattro enti di ricerca, in possesso di competenze riconosciute sulle tematiche trattate: l’Istituto di Scienze Marine del Consiglio Nazionale delle Ricerche (ISMAR), con la collaborazione dell’Università di Bologna (Dipartimento di Fisica ed Astronomia); l’Istituto Nazionale di Geofisica e Vulcanologia (INGV); il Centro Europeo per la Formazione e la Ricerca nell’Ingegneria dei Terremoti (EUCENTRE); la Rete Universitaria dei Laboratori di Ingegneria Sismica (ReLUIS).
Linee di ricerca
Le linee di ricerca (L) sviluppate nel corso del progetto sono le seguenti:
1. identificazione di faglie potenzialmente attive nei pressi delle piattaforme offshore in produzione; 2. identificazione di faglie potenzialmente sismogeniche nei pressi delle piattaforme offshore in
produzione e caratterizzazione della loro attività; 3. identificazione di frane sottomarine e instabilità di sedimenti lungo i margini continentali, che possano
essere attivati da scuotimento sismico generato dalle faglie potenzialmente sismogeniche; 4. modellazione dell’onda di maremoto potenzialmente prodotta dalle frane sottomarine; 5. modellazione dell’onda di maremoto potenzialmente prodotta dall’attivazione delle faglie selezionate
per gli scenari di danno; 6. modellazione di scenari di danno su strutture e infrastrutture costiere, e conseguenti perdite dovute ai
terremoti generati dalle faglie di scenario (incluso webGIS dedicato); 7. caratterizzazione della vulnerabilità delle strutture residenziali costiere rispetto all’impatto di un’onda di
maremoto; 8. modellazione di scenari di danno, dovuti ai maremoti sopra menzionati, su strutture residenziali
costiere, e conseguenti perdite.
Le indagini relative alle diverse linee di ricerca sono riferite alla linea di costa italiana che fronteggia, in senso lato, le piattaforme offshore in produzione considerate. Tale linea di costa è stata opportunamente
PROGETTO SPOT - Report integrato di fine progetto
5
divisa in cinque diverse aree di studio: Adriatico Settentrionale, Adriatico Centrale, Adriatico Meridionale, Mar Ionio (offshore della Calabria ionica), Canale di Sicilia (Figura 1).
Fig. 1 – Aree di studio del progetto: Adriatico Settentrionale, Adriatico Centrale, Adriatico Meridionale, Mar Ionio (offshore della Calabria ionica), Canale di Sicilia.
Risultati raggiunti
Linea di ricerca 1. Identificazione di faglie potenzialmente attive nei pressi delle piattaforme offshore in produzione
Nell’ambito di questa linea di ricerca, al fine di ricostruire l’assetto stratigrafico-strutturale del sottofondo marino, il CNR-ISMAR ha analizzato circa: i) 90.000 km di profili di sismica a riflessione monocanale di tipo CHIRP (penetrazione 30–50 m e risoluzione verticale decimetrica); ii) circa 45.000 km2 di dati batimetrici da ecoscandagli multifascio ad alta risoluzione (risoluzione spaziale variabile da pochi metri a 90 m); iii) circa 25.000 km di sismica a riflessione multicanale provenienti da acquisizioni ministeriali, progetto CROP e proprie (penetrazione 1–4 km e risoluzione verticale metrica o decametrica); iv) svariati campioni di sedimento (disponibili solo lungo la scarpata del Bacino di Gela e nell’Adriatico Centrale) e pozzi per l’esplorazione petrolifera, disponibili sul database ViDEPI accessibile alla pagina web https://www.videpi.com/videpi/videpi.asp#) migliaia di chilometri di sismica multicanale messi a disposizione da Eni S.p.A. durante due sessioni di lavoro in data room presso gli uffici di San Donato
PROGETTO SPOT - Report integrato di fine progetto
6
Milanese (MI). Sono state considerate tutte le eventuali faglie attive, le deformazioni del fondo mare e le evidenze di instabilità di scarpata nelle aree di progetto.
L’area attorno ai pozzi Gela-Panda, nell’area di studio del Canale di Sicilia, ospita parzialmente la Falda di Gela (Argnani, 1987), il più esterno e giovane sovrascorrimento del sistema di catena appenninico in Sicilia, che si è imbricato sopra l’avampaese pelagiano tra il tardo Miocene e il Pleistocene (Figura L1-1A). Le faglie più importanti nell’area delimiterebbero dei bacini estensionali del tardo Miocene, e sarebbero state riattivate in compressione durante il Pliocene-Pleistocene (Ghisetti et al., 2009). Tale sistema di faglie è ben evidente in Figura L1-1B, ma non risulta visibile nei profili CHIRP CNR-ISMAR, che mostrano la parte recente successione sedimentaria e il fondo mare (Figura L1-1C). Nella parte orientale del bacino, in prossimità della scarpata che lo borda e dove sono ubicate le frane gemelle (Minisini et al., 2007; Twin Slides in Figura L1-1A), è stata individuata una faglia del sistema flessurale che risulta interessare il fondo mare con una piccola scarpata. Va tuttavia notato che nicchie di frana di dimensioni simili sono distribuite lungo tutta la scarpata del bacino, anche lontano dall'emersione delle faglie (Figura L1-2B), e non risultano dunque rappresentare una condizione necessaria alla destabilizzazione.
Figura L1-1 – A: mappa batimetrica e strutturale dell’area di studio attorno ai pozzi Gela-Panda. Sono indicate le tracce dei profili mostrati in B e C. I cerchi rossi indicano l’area di interesse attorno ai pozzi. B: profilo sismico a riflessione multicanale GS53 acquisito nel 1990 dal CNR-ISMAR, che mostra il sistema di faglie principale che caratterizza l’area. C: profilo sismico a riflessione monocanale CHIRP, che attraversa tutta l’area di studio e mette in evidenza l’assenza di dislocazioni del fondo mare. TWT = Two Way Time, indica i tempi di andata e ritorno del segnale acustico dalla sorgente al sottofondo e viceversa.
PROGETTO SPOT - Report integrato di fine progetto
7
L’area di studio nel Canale di Sicilia, che include l’impianto Vega, si trova nell’immediato offshore di un sistema di faglie denominato Scicli-Ragusa (Figura L1-2A), prevalentemente considerate trascorrenti destre (Grasso et al., 2000), con una possibile riattivazione tardo quaternaria in trascorrenza sinistra (Catalano et al., 2008). La successione plio-quaternaria della piattaforma continentale, che è presente con spessori molto ridotti e talora è assente, non sembra essere interessata da deformazioni importanti, nemmeno sulla prosecuzione a mare delle strutture di Scicli-Ragusa e di Comiso, mentre mostra una complessa relazione fra erosione e sedimentazione (Figura L1-2B). In nessuno dei profili sismici CHIRP di CNR-ISMAR è stato possibile rintracciare evidenza a fondo mare della prosecuzione a mare delle strutture di Scicli-Ragusa (Fig. L1-2B a titolo esemplificativo) e del sistema Comiso-Chiaramonte. Il substrato carbonatico mostra invece la presenza di ampie pieghe chilometriche ad asse subverticale che hanno direzione da N–S a NNE–SSO, che sono probabilmente il risultato di un'inversione paleogenica e che non mostrano attività successiva (Argnani e Torelli, 2001).
L’area attorno ai campi Luna, nell’offshore crotonese, è piuttosto complessa da un punto di vista morfo-batimetrico e strutturale. La piattaforma continentale è larga al massimo 7 km, la scarpata superiore ha una pendenza molto accentuata che può raggiungere localmente 20°, ed è caratterizzata dalla fitta presenza di canyon e canalizzazioni erosive (Figura L1-3A). La scarpata intermedia è caratterizzata dalla presenza a fondo mare di numerose evidenze di espulsioni di fluidi (pockmark; Figura L1-3A), che sono collegate a numerose faglie estensionali, presenti nella successione plio-quaternaria; queste ultime hanno favorito la migrazione di gas verso il fondo mare a partire dalla successione miocenica, che rappresenta la roccia serbatoio ed è posta a 1262 m di profondità sotto il livello del mare (Figura L1-3B). Sono abbondanti le
Figura L1-2 – A: mappa batimetrica e strutturale dell’area di studio attorno ai pozzi Vega. B: profilo sismico a riflessione monocanale CHIRP, la cui traccia è mostrata in A, che attraversa tutta la piattaforma e mette in evidenza l’assenza di dislocazioni del fondo mare. Sono indicate le superfici erosive ES1, ES2 ed ES3 riferite agli ultimi stazionamenti bassi del livello del mare di circa 29.000, 123.000 e 250.000 anni fa.
PROGETTO SPOT - Report integrato di fine progetto
8
frane (Figura L1-3C), probabilmente favorite nel loro innesco dalla presenza di fluidi, dagli alti gradienti di scarpata e da una morfologia caratterizzata da zone maggiormente rilevate e, quindi, più facilmente soggette a scivolamento gravitativo. L'analisi dei profili dell'Eni S.p.A. ha consentito di individuare strutture con direzione NNO-SSE parallele ad allineamenti di scarpate che hanno espressione morfologica (Figura L1-3A) e hanno confermato la presenza di una diffusa deformazione dovuta a faglie normali, di probabile origine gravitativa, che spesso sembrano raccordarsi alla base della successione pliocenico-quaternaria e talvolta a quella pre-pliocenica.
Gli impianti del pozzo Aquila sono ubicati al piede della ripida scarpata erosionale che caratterizza il margine orientale della Piattaforma Apula (Figura L1-4A), nell’area di studio dell’Adriatico Meridionale. Nell'area non sono presenti strutture tettoniche importanti di età neogenica (Figura L1-4B, C). Al margine della Piattaforma Apula si evidenziano fenomeni di instabilità gravitativa nei sedimenti sovrastanti (Figura L1-4C), legati a compattazione differenziale fra i sedimenti clastici di spessore ridotto, posti sopra ai carbonati mesozoici di piattaforma, e la spessa successione clastica terziaria, presente sopra all'adiacente bacino pelagico mesozoico (Figura L1-4B). Nella parte più orientale dell'area sono presenti strutture di inversione che interessano faglie estensionali, probabilmente originatesi per la flessurazione dell'avampaese, sottoposto al carico della catena ellenico-albanese (Figura L1-4B). Tali strutture di inversione deformano la superficie di erosione messiniana, ma non interessano i sedimenti quaternari, come suggerisce la mappa (non mostrata) computata per il riflettore corrispondente allo stazionamento basso di 800.000 anni fa (800 ka in Figura L1-4B), che conferma l'assenza di strutture tettoniche attive nel Quaternario.
Figura L1-3 – A: mappa batimetrica e strutturale dell’area di studio attorno al campo Luna. B: profilo sismico multicanale ViDEPI. C: profilo sismico monocanale CHIRP che mette in evidenza la presenza di frane a fondo mare. Le tracce dei profili sono mostrate in A.
PROGETTO SPOT - Report integrato di fine progetto
9
Nell’area di studio dell’Adriatico Centrale, gli impianti di produzione attorno al pozzo Rospo Mare sono ubicati sopra la Piattaforma Apula di età mesozoica, la quale rappresenta l’antico margine passivo del dominio adriatico, che si estende sulla maggior parte dell'area (Figura L1-5). Qui il fronte appenninico si trova a terra e i sedimenti dell'avanfossa plio-quaternaria hanno uno spessore ridotto. Al margine della Piattaforma Apula si osserva una limitata instabilità gravitativa dovuta a compattazione differenziale, che talora si manifesta con faglie estensionali a limitato rigetto che si radicano al ciglio della piattaforma. Nell'insieme, non sono state rilevate strutture significative che risultino attive durante il Quaternario e che interessino il fondo mare.
Figura L1-4 – A: mappa batimetrica e strutturale dell’area di studio attorno al pozzo Aquila. B. profilo sismico multicanale ministeriale. C: profilo sismico monocanale CHIRP che mostra la presenza di una frana sepolta da circa 50 m di sedimento.
PROGETTO SPOT - Report integrato di fine progetto
10
Figura L1-5 – profondità della superficie di base del Pliocene mappata nell’Adriatico Centrale e Settentrionale, con sovraimposte le strutture tettoniche principali ritenute attive alla base del Pliocene (in rosso) e alla base del Quaternario (in verde).
Nell'area di studio dell’Adriatico Centrale sono presenti varie strutture tettoniche neogeniche (Figura L1-5). Nella parte meridionale dell’area, il fronte appenninico è ubicato all’incirca sulla costa (campo Rospo Mare in Figura L1-5), per poi spostarsi verso mare nel settore settentrionale. Uno spesso bacino plio-quaternario di avanfossa (Bacino di Pescara, Figura L1-5) è presente nella parte meridionale, mentre si osserva una forte riduzione degli spessori dei sedimenti plio-quaternari verso nord (Figura L1-5). I fronti di sovrascorrimento esterni appaiono sigillati dalla sedimentazione quaternaria. Nel complesso, non si rilevano nell'area marina importanti strutture tettoniche che mostrino attività dal tardo Pleistocene a oggi, con la possibile eccezione della struttura costiera, che risulta però difficile da valutare con i soli dati disponibili a mare.
PROGETTO SPOT - Report integrato di fine progetto
11
Una caratteristica particolarmente rilevante nell'Adriatico centro-settentrionale è data dalla presenza di sistemi di progradazione plio-quaternari estesi su tutta l'area, che rappresentano un utile elemento stratigrafico per stimare l'attività recente delle strutture tettoniche. La mappatura della base e del tetto dell’inviluppo di questi sistemi progradazionali, che hanno la geometria di clinoformi, è stata utilizzata proprio per evidenziare lo stile di sedimentazione, ed eventuale deformazione o dislocazione indotta da attività tettonica, in un intervallo temporale compreso tra Pleistocene inferiore (base) e superiore (tetto). La superficie di base delle unità mostra limitate evidenze di attività deformativa o dislocativa (Figura L1-6A), mentre sulla superficie di tetto non sono visibili perturbazioni imputabili ad attività tettonica (Figura L1-6B).
Figura L1-6 – A: mappatura della base dell’inviluppo dei clinoformi quaternari. B: tetto dei clinoformi (profondità in millisecondi). In blu sono rappresentate le strutture tettoniche ritenute attive fino alla base del Pliocene, in verde quelle attive fino alla base del Quaternario. I cerchi rossi indicano le aree di progetto SPOT.
PROGETTO SPOT - Report integrato di fine progetto
12
Nell'area di studio dell’Adriatico Settentrionale sono presenti varie strutture neogeniche (Figura L1-5). Nella parte più settentrionale dell’area, in prossimità della costa, sono presenti i sovrascorrimenti esterni appenninici (Fantoni e Franciosi, 2010), che appaiono sigillati dal sistema di progradazione quaternaria (es.: Ghielmi et al., 2013). Nel settore adiacente la costa di Cattolica (Figura L1-6C), i sistemi progradanti e le sovrastanti unità in aggradazione mostrano, a grande scala, superfici progressivamente inclinate verso est, che indicano una possibile crescita delle strutture ubicate in prossimità della costa, in maniera simile a quanto già osservato in precedenza nell'Adriatico Centrale. Nella parte sud-orientale dell'area di studio (attorno ai pozzi del campo Clara) sono presenti blande pieghe che interessano l'avampaese e sono assimilabili al sistema deformativo presente nell'Adriatico Centrale (es.: Argnani e Frugoni, 1987). In questo caso, tuttavia, i sedimenti quaternari non sono interessati da deformazione, e questo si vede sia nei dati di sismica multicanale (es.: Pace, 2018), sia nei dati di sismica monocanale ad alta risoluzione di CNR-ISMAR. Inoltre, anche nell’Adriatico Settentrionale sono stati mappati la base e il tetto dei clinoformi (Figura L1-6C e D) e, nel complesso, non sono presenti nell'area strutture tettoniche che mostrino attività nel Pleistocene.
Linea di ricerca 2. Identificazione di faglie potenzialmente sismogeniche nei pressi delle piattaforme offshore in produzione e caratterizzazione della loro attività
Identificazione di potenziali faglie sismogeniche
Le cinque Aree di Studio riportate nella Figura 1 sono state oggetto di uno studio approfondito, condotto dall’INGV, finalizzato all’individuazione di faglie potenzialmente sismogeniche (Figura L2-1). Le aree sono caratterizzate da regimi tettonici e stili strutturali molto diversi tra di loro, come emerso anche dall’analisi di dettaglio della letteratura scientifica disponibile; tale analisi è iniziata con uno studio bibliografico approfondito, a cui ha fatto seguito la consultazione di log di pozzo e dati geologici ricavabili da relazioni redatte nell’ambito di permessi di ricerca e concessioni di coltivazione, disponibili sul database ViDEPI. Sono state, inoltre, interpretate diverse migliaia di km di profili di sismica a riflessione, disponibili in formato SEGY nel database SNAP (Seismic data Network Access Point; Diviacco et al., 2018), sotto licenza dell’Istituto Nazionale di Oceanografia e di Geofisica Sperimentale (OGS, Trieste). Quanto emerso dallo studio effettuato è stato ulteriormente vincolato attraverso la consultazione di dati confidenziali (linee sismiche ad elevata risoluzione e log di pozzo), avvenuta durante due sessioni di data room presso gli uffici di ENI S.p.A. di San Donato Milanese (MI).
PROGETTO SPOT - Report integrato di fine progetto
13
Figura L2-1 – Mappe dei maggiori lineamenti strutturali riferiti alla Base del Pliocene/tetto del Messiniano identificati nelle cinque aree di studio: AS Adriatico Settentrionale, AC Adriatico Centrale, BR Adriatico Meridionale, KR Mar Ionio (offshore della Calabria ionica), CS Canale di Sicilia.
Sono stati individuati e caratterizzati da un punto di vista geometrico e, ove possibile, cinematico, i principali piani di faglia riconosciuti all’interno delle aree di studio, e mappati alcuni orizzonti stratigrafici di importanza regionale; in particolare, la base del Pliocene/tetto del Messiniano è stata individuata su tutte le aree e a tale orizzonte sono state riferite le carte strutturali relative alle faglie principali (Figura L2-1). I piani di faglia, inizialmente individuati nei profili sismici in tempi doppi (TWT, si veda ad esempio la Figura L1-1B), sono stati convertiti in profondità utilizzando modelli di velocità derivati dai log di pozzo e dai profili sismici multicanale, o ricavati dalla letteratura scientifica. Sono stati infine analizzati la sismicità storica (CPTI15, Rovida et al., 2016; Guidoboni et al., 2018) e strumentale, consultabili nei principali database INGV, e i dati batimetrici disponibili sul portale “EMODnet Bathymetry” (http://www.emodnet-bathymetry.eu; EMODnet Bathymetry Consortium, 2018).
L’Adriatico Settentrionale (AS in Figura L2-1) include i fronti esterni della catena Nord-Appenninica, caratterizzati da uno stile strutturale a pieghe e faglie noto in letteratura. Sono stati ricostruiti 19 piani di faglia principali per una lunghezza complessiva di 389 km e, su gran parte dell’area di studio, l’orizzonte stratigrafico relativo alla base del Pliocene, rispetto al quale sono stati stimati i rigetti. In generale, le faglie ricostruite si radicano all’interno di livelli stratigrafici pre-pliocenici, raggiungendo, in alcuni casi, i livelli evaporitici triassici, ad eccezione di due piani di faglia minori scollati alla base della successione pliocenica.
PROGETTO SPOT - Report integrato di fine progetto
14
Tutte le faglie mappate mostrano una cinematica inversa (thrust), sono sepolte e confinate all’interno della successione pliocenico-quaternaria e, ad eccezione di una, dislocano la base del Pliocene.
La sismicità strumentale registrata nell’area tra il 2005 ed il 2015 è maggiormente concentrata a ridosso della linea di costa, e presenta magnitudo comprese tra 1.2 e 4.5, con una preponderanza di eventi di magnitudo 1.6–1.7. La profondità degli eventi varia tra 0 e 40 km, ma la maggior parte è confinata nei primi 20 km. L’evento di magnitudo più alta è stato registrato nel giugno 2012, è localizzato nell’immediato offshore di Ravenna a una profondità di 26 km e mostra un meccanismo focale con cinematica trascorrente (http://cnt.rm.ingv.it/event/925601).
L’Adriatico Centrale (AC in Figura L2-1) mostra uno stile strutturale paragonabile a quello dell’Adriatico Settentrionale. L’analisi delle linee sismiche ha evidenziato 25 strutture compressive, per una lunghezza complessiva di 366 km, che si scollano su due differenti livelli stratigrafici: uno superficiale, localizzato alla base della successione pliocenica, e uno più profondo, all’interno della successione evaporitica triassica. Sono state inoltre riconosciute 4 strutture con cinematica trascorrente e geometria del piano di faglia sub-verticale, per una lunghezza di 56 km. Tutte le strutture compressive sono blind thrust (sovrascorrimenti ciechi), in quanto la loro terminazione superficiale si colloca all’interno della successione pliocenica e non raggiunge il fondo mare. Tuttavia, gli orizzonti pliocenici e la base del Quaternario mostrano deformazioni a pieghe legate all’attività delle strutture compressive sottostanti. Le strutture compressive sono state quindi suddivise in due categorie: quelle che scollano alla base delle successioni silicoclastiche e quelle che si radicano nei carbonati. I thrust radicati nei carbonati mesozoici mostrano occasionalmente l’associazione con strutture tipo back-thrust (retroscorrimenti) a vergenza occidentale.
La sismicità strumentale registrata nell’area tra il 2005 ed il 2015 è sparsa e mostra magnitudo < 3.4. Le profondità variano tra 0 e 46 km; tuttavia, oltre il 30% della sismicità è collocato nei primi 10 km, mentre una seconda distribuzione importante (circa il 20%) si rinviene a circa 30 km di profondità.
L’Adriatico Meridionale (BR in Figura L2-1) si colloca all’interno di un’area di avampaese (Piattaforma Apula) del sistema orogenico appenninico. Sono stati interpretati e mappati l’orizzonte del tetto del Messiniano e quello al tetto dei depositi carbonatici cretacei; quest’ultimo rappresenta l’orizzonte più profondo individuabile dai dati sismici a disposizione. I rigetti che dislocano questo orizzonte rappresentano elementi di natura tettonica o paleo-morfologica legati all’evoluzione mesozoica della stessa formazione, del tutto suturati dalla sedimentazione successiva e, ad oggi, non attivi nell’attuale contesto di avampaese. Pertanto, il margine della Piattaforma Apula appare come elemento passivo già a partire dal Giurassico medio. L’orizzonte tetto del Messiniano mostra una generale continuità su tutta l’area di interesse ed è dislocato da due faglie dirette coniugate, orientate ONO-ESE, a cui è associato un rigetto minimo. Non sono pertanto state riconosciute faglie potenzialmente sismogeniche all’interno di tale area di studio.
Inoltre, dall’analisi della sismicità storica e strumentale è emersa l’assenza di eventi rilevanti.
La Calabria Ionica (KR in Figura L2-1) si colloca in un settore dell’orogene appenninico caratterizzato dal passaggio tra la zona di subduzione calabra e le aree di collisione continentale con il margine apulo. L’area risulta interessata da una tettonica polifasica, che riflette i vari stadi di evoluzione della subduzione della crosta ionica e la successiva collisione continentale con il margine apulo. L'elemento strutturale dominante consiste in un cuneo di accrezione che si estende nel Mar Ionio fino alla latitudine di Malta. L’analisi delle linee sismiche interpretate ha permesso di riconoscere strutture compressive databili al Miocene e al Pliocene, e distensive/transtensive ascrivibili al Pleistocene. Le strutture compressive sono sia legate all’impilamento delle unità tettoniche del cuneo di accrezione (sovrascorrimenti) sia a strutture di tipo pop-up (combinazione di scorrimenti e retroscorrimenti), probabilmente connesse alle fasi tettoniche transpressive. Le strutture principali consistono in thrust che non dislocano la discordanza zancleana (Pliocene inferiore). Tali strutture si inseriscono in un contesto più generale in cui si evidenziano sollevamenti attivi durante il Pleistocene superiore e l’Olocene. Le strutture che mostrano evidenze di attività nel Plio-Pleistocene sono legate ad elementi estensionali/transtensivi che tagliano i precedenti elementi compressivi. Complessivamente, le strutture di maggior interesse che sono state mappate sono: 2 thrust superficiali, per una lunghezza totale di 53 km, 5 thrust al di sotto della discordanza zancleana, per una lunghezza totale di 102 km, e un elemento trascorrente di 35 km di lunghezza.
PROGETTO SPOT - Report integrato di fine progetto
15
L’analisi della sismicità, basata sui dati offshore dal 2005 al 2015 rilocalizzati da Anselmi et al. (2018), ha permesso l’identificazione di pochi eventi (23) con magnitudo compresa tra 1.3 e 3.6, con un picco per eventi di magnitudo 3, e profondità che varia tra 6 e 30 km. Complessivamente il 40% della sismicità è collocato tra 5 e 10 km di profondità, mentre il restante 60% è collocato tra i 20 ed i 30 km.
Il Canale di Sicilia (CS in Figura L2-1) ricade, nel quadro geodinamico del Mediterraneo Centrale, all’interno della porzione più settentrionale dell’avampaese africano, il cui assetto strutturale attuale è dominato dal processo di collisione continentale neogenica con la placca europea. All’interno di questa area di studio sono state mappate 24 faglie principali, per una lunghezza complessiva di 844 km.
In particolare, il settore più orientale (offshore di Pozzallo) è dominato dalla prosecuzione offshore della Linea di Scicli, qui caratterizzata a terra da diversi segmenti di faglia con cinematica di tipo trascorrente (Grasso et al., 2000; Figura L1-2A); complessivamente, essa si estende in direzione NNE-SSO per oltre 65 km ed è stata indentificata fino a oltre 10 km di profondità, producendo strutture a fiore che si estendono lateralmente anche per oltre 10 km. Alcuni segmenti di faglia orientati NE-SO, lunghi poche decine di km e identificati fino a profondità di circa 8 km, convergono verso la prosecuzione offshore della Linea di Scicli e si estendono nelle aree emerse, dove sono stati definiti come probabili sorgenti sismogeniche (Azzaro e Barbano, 2000).
Il settore occidentale (offshore di Gela) è invece contraddistinto dal fronte più avanzato della Falda di Gela e da un sistema di faglie generalmente distensive, localmente riattivate in compressione. Queste sono orientate circa NO-SE nelle porzioni occidentali e ENE-OSO in quelle orientali, e interessano l’intera successione mesozoico-miocenica fino all’orizzonte corrispondente al tetto del Messiniano. La successione plio-quaternaria non mostra, in genere, evidenze di tettonica attiva ascrivibile alle suddette faglie, ma solo rare strutture superficiali di minore entità. Le uniche deformazioni a fondo mare, oltre a quelle legate a locali fenomeni franosi e/o erosivo-deposizionali, sono riconducibili a fenomeni di compattazione differenziale dei sedimenti in corrispondenza dei rigetti maggiori dell’orizzonte tetto del Messiniano.
L’analisi della sismicità locale, relativa al periodo 1981-2017, ha permesso di identificare circa 1600 eventi con magnitudo compresa tra 0.4 e 4.5. La sismicità risulta confinata nei primi 30 km di profondità, con pochi eventi localizzati al di sotto di tale limite. In particolare, sono state evidenziate due sequenze sismiche che hanno interessato l’offshore di Pozzallo nell’ottobre del 1990 e nel dicembre del 2017, che indicano movimenti prevalentemente di tipo trascorrente sinistro lungo piani sub-paralleli alla prosecuzione offshore della Linea di Scicli. Il Catalogo dei Forti Terremoti in Italia (Guidoboni et al., 2018) non riporta alcun evento nell’intera area di studio.
Caratterizzazione dell’attività delle faglie potenzialmente sismogeniche identificate
La caratterizzazione della possibile attività delle strutture tettoniche identificate è stata condotta dall’INGV utilizzando tre diversi metodi: calcolo della propensione al movimento (slip tendency; Morris et al., 1996; Neves et al., 2009), calcolo dei tassi di movimento di lungo termine, calcolo dei tassi di sismicità.
La propensione al movimento delle faglie all’interno dei campi di sforzi a cui sono sottoposte attualmente è rappresentata da un valore compreso tra 0 e 1; al crescere di tale valore aumenta la probabilità che si abbia movimento lungo la struttura. L’inizio del movimento avviene in genere per valori di slip tendency uguali o superiori al coefficiente di attrito statico della roccia (mediamente pari a 0.6; Byerlee, 1978).
Il calcolo di slip tendency è fortemente influenzato dal regime tettonico in cui si trova la faglia, dai valori degli assi degli sforzi e dal loro rapporto reciproco. Considerata l’estensione dell’area mappata, il calcolo di slip tendency è stato applicato in maniera semplificata a scala regionale, adottando delle assunzioni che permettessero un confronto dei risultati per tutte le aree di studio. I risultati ottenuti costituiscono, quindi, un livello base di conoscenze che rende possibile valutare quali strutture possano essere di interesse per ulteriori approfondimenti e studi specifici. Le aree di studio sono state suddivise in settori più piccoli, all’interno dei quali si può assumere che il regime tettonico e l’orientazione degli sforzi agenti sia omogenea. Per ogni faglia ricadente in ogni settore sono stati calcolati i valori di slip tendency al variare del
PROGETTO SPOT - Report integrato di fine progetto
16
rapporto dei valori tra i tre assi degli sforzi. Nella Tabella L2-1 vengono elencati i risultati di sintesi in cui è riportato il valore medio di slip tendency per ogni settore.
Tabella L2-1 – Sintesi dei risultati dei valori di slip tendency media.
Area Adriatico Settentrionale Adriatico Centrale Calabria Ionica Canale di Sicilia
Il calcolo dei tassi di movimento di lungo termine è stato fatto per una singola struttura nell’ambito del progetto SPOT, per la quale le elaborazioni congiunte dei gruppi di lavoro INGV e CNR-ISMAR lo rendevano possibile. La struttura scelta come esempio si colloca all’interno dell’Adriatico Centrale, a largo di San Benedetto del Tronto; si tratta di una piega anticlinale ben evidente negli orizzonti di sottosuolo mappati sul dato sismico. Sono state studiate le deformazioni registrate da un orizzonte stratigrafico datato con metodi stratigrafici e radiometrici a 800.000 fa (equivalente a quello mostrato in Figura L1-4B), e si è preventivamente tenuto conto dell’effetto della compattazione dei depositi sottostanti. Una volta eliminato l’effetto di compattazione, infatti, la deformazione tettonica residua può essere imputata alle strutture tettoniche presenti nell’area. La procedura di retro-deformazione permette di calcolare l’entità di movimento delle faglie per produrre la deformazione tettonica osservata. Tenendo conto che tale deformazione è successiva alla deposizione dell’orizzonte stesso, si può calcolare il tasso medio di movimento della faglia (o delle faglie) studiate. Nel caso d’esempio sono state analizzate separatamente due strutture, una più superficiale ed una più profonda, che potrebbero aver causato la deformazione osservata, ottenendo un valore di slip rate medio di 0.3 mm/a riferito alla struttura più profonda e 0.4 mm/a per la struttura superficiale negli ultimi 800.000 anni. Tuttavia, non è possibile escludere che entrambe le strutture abbiano contribuito, in parte e con tassi di movimento inferiori, a produrre le deformazioni osservate.
Il terzo metodo utilizzato per la caratterizzazione delle strutture tettoniche è la valutazione dei tassi di sismicità nelle varie aree di studio. A tale scopo è stato utilizzato come dato di partenza il Catalogo Parametrico dei Terremoti Italiani (CPTI15; Rovida et al., 2016), con copertura temporale dal 1000 al 2014 e comprendente un totale di 4584 terremoti con stima della magnitudo momento (Figura L4-2, mappa di sinistra). Il tasso di sismicità è rappresentato da una distribuzione frequenza-magnitudo secondo due modelli, quello Gutenberg–Richter caratteristico e quello Gutenberg–Richter con decadimento graduale, detto tapered (Kagan, 2002). Per ottenere delle stime robuste con tale metodo è necessario avere un elevato numero di eventi distribuito su un ampio intervallo di magnitudo nel catalogo dei terremoti. A questo scopo abbiamo definito delle macro-regioni, alla scala di tutta la penisola italiana, per raggruppare la sismicità relativa al medesimo contesto tettonico. Assumendo che la sismicità di lungo termine possa avere una distribuzione circa uniforme in tutta l’area su cui si è basato il calcolo, si può ottenere il tasso annuale di sismicità ancorando le distribuzioni modellate alla distribuzione di sismicità osservata nella singola area di studio. Le aree di studio che si trovano a cavallo tra regioni caratterizzate da strutture tettoniche diverse sono state suddivise e trattate separatamente. Le distribuzioni frequenza-magnitudo così ottenute sono mostrate nei grafici di destra della Figura L2-2.
PROGETTO SPOT - Report integrato di fine progetto
17
Figura L2-2 – A sinistra, mappa della sismicità dal catalogo CPTI15 (Rovida et al., 2016) con l’indicazione delle aree di indagine. Le aree di studio sono identificate con le seguenti sigle: AS, Adriatico Settentrionale; AC, Adriatico Centrale; AM, Adriatico Meridionale; CI, Calabria Ionica, CS, Canale di Sicilia. Le suddivisioni ulteriori, mostrate in mappa in grigio, sono identificate con lettere minuscole e numeri. A destra, diagrammi della distribuzione frequenza-magnitudo cumulativa, per magnitudo uguale o maggiore a 4.5, nelle varie aree di studio, o loro suddivisioni, nelle sue varianti caratteristica (CHR) e tapered (TGR), ottenute per l’analisi di completezza storica (Hist) e statistica (Stat). La distribuzione ASb2 viene mostrata solo per la versione con completezza statistica. Nei grafici, le distribuzioni CSa e CSb sono quasi completamente coincidenti.
Scenari e loro caratteristiche
Nelle analisi di pericolosità, l’utilizzo degli scenari è molto diffuso come strumento complementare alle analisi probabilistiche. Queste ultime cercano di esplorare completamente il campo di variabilità di un fenomeno, mentre gli scenari esplorano alcuni aspetti particolari di un fenomeno in condizioni specificate.
In questo progetto si è convenuto di realizzare undici scenari (Figura L2-3) che si ritiene possano rappresentare eventi “tipo” della sismicità naturale in prossimità dei pozzi produttivi oggetto di questo studio (Figura 1). Tali scenari sono distribuiti su tutte le aree di studio, comprendendo sia casi nel campo vicino, sia casi nel campo lontano rispetto ai siti di interesse posti sulla costa. Essi sono stati inoltre individuati per rappresentare la variabilità dei contesti tettonici interessati e dei possibili meccanismi di rottura, coerentemente con i meccanismi delle faglie studiate nel progetto. Viceversa, la magnitudo momento (Mw) del terremoto di scenario è stata prefissata e posta pari a 6.5 per tutti gli scenari, al fine di poter effettuare dei confronti tra scenari a partire dallo stesso quantitativo di energia rilasciata. Le dimensioni delle rotture (lunghezza, larghezza e spostamento), sono state fissate sulla base dei valori attesi da relazioni di scala in ragione della magnitudo prefissata, del contesto tettonico e del meccanismo di rottura (Leonard, 2014). I parametri così ottenuti costituiscono il dato di ingresso per stimare lo scuotimento utilizzando delle leggi di attenuazione e per effettuare le simulazioni numeriche per gli tsunami.
PROGETTO SPOT - Report integrato di fine progetto
18
Figura L2-3 – Mappa generale delle rotture (faglie a geometria planare rettangolare) per gli scenari di scuotimento e di tsunami in relazione alle aree d’indagine. I parametri dei singoli scenari sono riportati in Tabella L2-2.
La Tabella L2-2 mostra i valori dei parametri dei piani di rottura cosismica (Figura L2-3) che costituiscono il dato d’ingresso per gli scenari. Per ogni piano di rottura sono anche mostrate le stime di slip tendency medie, ottenute campionando i valori dell’area di rottura sulla faglia a cui essa è associata. Questi valori non determinano in assoluto la possibilità che uno scenario possa verificarsi, ma misurano la propensione di ogni singola faglia a riattivarsi. Essi quindi rappresentano la stima di una delle cause che concorrono alla riattivazione. Valori bassi di slip tendency non significano che il verificarsi dello scenario sia impossibile, ma solo che il contributo di tale fattore è più basso rispetto agli altri casi analizzati. Viceversa, un valore di slip tendency alto consente ragionevolmente di ritenere che la riattivazione possa avvenire con un contributo inferiore di altri fattori.
La scelta di una singola magnitudo prefissata comporta che ci si possa attendere una diversa frequenza di accadimento degli scenari nelle diverse aree. Il tasso di sismicità di lungo termine, assunto come stazionario, è rappresentato dalla distribuzione frequenza-magnitudo nelle varie aree d’indagine (Figura L2-2) ed è la base su cui abbiamo stimato la frequenza relativa degli scenari secondo due alternative (Tabella L2-2): (i) che la magnitudo dello scenario venga superata; (ii) che la magnitudo dello scenario sia compresa tra quella prefissata e una appena superiore. Questi valori consentono di valutare la possibilità di accadimento di uno scenario, o gruppi di scenari, rispetto agli altri.
Tabella L2-2 – Parametri dei piani di rottura per gli scenari di scuotimento e tsunami (le sigle sono le medesime di
Figura L2-3). I valori di latitudine (Lat), longitudine (Lon) e profondità (depth) sono riferiti al baricentro del piano di
rottura. L: lunghezza, W: larghezza, D: spostamento. Sl. T.: slip tendency. Stress R e T: campo di sforzo rispettivamente
per meccanismo inverso e trascorrente. I valori riportati nelle sei colonne di destra rappresentano la frequenza
annuale normalizzata degli scenari considerati per due casi: magnitudo momento uguale o maggiore a 6.5, e per
magnitudo momento compresa tra 6.5 e 6.6, ricavate sulla base di due modelli di distribuzione frequenza-magnitudo,
CHR: Gutenberg-Richter caratteristica, TGR: Gutenberg-Richter tapered, adottando due stime di completezza, Hist:
PROGETTO SPOT - Report integrato di fine progetto
19
storica, Stat: statistica. Si tratta quindi di una misura relativa di quanto frequentemente potrebbe verificarsi uno
scenario in un’area di studio rispetto a uno scenario in un’altra.
Linea di ricerca 3. Identificazione di frane sottomarine e instabilità di sedimenti lungo i margini continentali, che possano essere attivati da scuotimento sismico generato dalle faglie sismogeniche
Attraverso l’analisi dei dati morfobatimetrici e di sismica a riflessione ad alta risoluzione, sono state riconosciute diverse evidenze di trasporto di massa lungo le scarpate continentali, ma solo su due aree si è concentrato lo studio della pericolosità legata a frane sottomarine capaci di generare maremoti, e cioè il Bacino di Gela (area di studio del Canale di Sicilia) e la zona a mare di Capo Rizzuto (area di studio nell’immediato offshore della Calabria ionica; Figura L3-1). Queste due aree coincidono con quelle dove si è riscontrata la presenza sia di frane sottomarine sia di faglie potenzialmente sismogeniche. Per ciascuno dei casi in esame è stata infatti calcolata la soglia di accelerazione di picco sufficiente ad attivare la frana ed è stata confrontata con l’accelerazione calcolata sulla base degli scenari di rottura sismica identificati nel progetto.
Figura L3-1 – Mappe di ubicazione e parametri dimensionali di ciascuna delle tre frane studiate e modellate per la loro capacità di generare onde di maremoto.
Nel Bacino di Gela, in parte sul fronte della Falda di Gela e circa a sud del pozzo Panda, è presente una grande frana messa in posto nel tardo Pleistocene e non più attiva, con una nicchia di distacco lunga più di 50 km e uno spessore di sedimenti di circa 700 m (Trincardi e Argnani, 1990; Figura L1-1A). La scarpata continentale del Bacino di Gela, invece, presenta numerosi e più recenti corpi di frana. La Northern Twin Slide (NTS) è la più grande di due frane adiacenti, verificatesi nelle ultime migliaia di anni (Minisini et al., 2007; Minisini e Trincardi, 2009), ultimi episodi di una sequenza il cui inizio viene datato circa 87.000 anni fa (Kuhlmann et al., 2017). Nel settore sud-orientale del Bacino di Gela, sono stati identificati 16 depositi di trasporto di massa esposti al fondo mare e 5 sepolti, con una superficie media di 30 km2. I 16 depositi esposti a fondo mare (con superficie variabile da 2 a 94 km2) sono stati mappati insieme alle loro nicchie di frana, che variano da 1 a 13 km di lunghezza, con una lunghezza media di 5 km. La più grande di queste frane è la Southern Gela Basin Slide, SGBS (FiguraL3-2). Le cause di instabilità di pendio non sono necessariamente riconducibili alla sismicità dell’area; si ritiene, infatti, che l'instabilità sia favorita dalla presenza di livelli vulcano-clastici con resistenza al taglio diversa rispetto ai sedimenti circostanti (Kuhlmann et., 2017) e di corpi sedimentari depositati dalle correnti di fondo, che aumentano il gradiente al ciglio della piattaforma continentale (Gauchery et al., 2019).
Nell'offshore di Capo Rizzuto, Mar Ionio, è stata individuata una possibile frana di grandi dimensioni (Capo Rizzuto Slide-CRS). Sebbene la documentazione della frana resti incompleta a causa della lacunosa copertura dei dati sismici, le sue grandi dimensioni (25 km3) la rendono di interesse per la modellazione tsunamigenica.
Per queste due aree, il Dipartimento di Fisica e Astronomia dell’Università di Bologna ha elaborato complessivamente tre scenari di maremoto, due per la prima area (-NTS e -SGBS) e uno per la seconda (CRS).
PROGETTO SPOT - Report integrato di fine progetto
21
Figura L3-2 – A: area della scarpata del Bacino di Gela interessata da frane sottomarine individuabili grazie ai dati batimetrici multibeam ad alta risoluzione. B: mappatura delle frane nel settore SE della scarpata. C: profilo sismico a riflessione CHIRP attraverso la frana SGBS.
Per l’analisi di stabilità del pendio, senza carico sismico, è stato usato il metodo dell’equilibrio limite,
modificato utilizzando il criterio della Minima Deviazione Litostatica (MDL) (Tinti e Manucci, 2006; Paparo e
Tinti, 2017). Per l’analisi della stabilità di un pendio sotto carico sismico, le accelerazioni sismiche sono
inserite nel modello in modo statico mediante i valori di PGA (Peak Ground Acceleration, Accelerazione di
Picco al Suolo). Per quanto riguarda la scelta dei valori di PGA, in questo studio si è tenuto conto sia delle
PROGETTO SPOT - Report integrato di fine progetto
22
accelerazioni definite nelle Mappe di Pericolosità Sismica italiana (MPS; Meletti e Montaldo, 2007) che
delle accelerazioni fornite da EUCENTRE e generate da 4 scenari di rottura identificati da INGV (CS_01,
CS_02, KR_01, KR_02), che si trovano in prossimità dei corpi di frana: NTS, SGBS e CRS. Per ognuno degli
scenari è stata quindi analizzata la stabilità con e senza carico sismico (Tabella L3-1).
Tabella L3-1 – Riepilogo delle PGA usate per l’analisi di stabilità delle frane individuate in prossimità della scarpata del
Bacino di Gela e al largo di Capo Rizzuto, per gli scenari di rottura associati ad un terremoto di Mw = 6.5. I valori medi
di PGA sono calcolati sul baricentro del corpo di frana prima del crollo.
PGA
NTS SGBS CRS
Carico sismico dedotto dalle Mappe di Pericolosità Sisimica
Probabilità di eccedenza del 10% in 50 anni (periodo di ritorno di 475 anni)
0.05 0.05 0.15
Probabilità di eccedenza del 2% in 50 anni (periodo di ritorno di 2475 anni)
0.15 0.10 0.25
Carico sismico definito nel progetto
Faglia CS_01 con Mw=6.5 0.05 0.08
Faglia CS_02 con Mw=6.5 0.34 0.04
Faglia KR_01 con Mw=6.5
0.17
Faglia KR_02 con Mw=6.5
0.06
Riassumendo, nel caso della frana NTS si ha instabilità solo nel caso di rottura della faglia CS_02, ma la frana
ha possibilità di innesco da carico sismico, anche in base ai valori riportati nella MPS. Nel caso SGBS, le PGA
definite nelle MPS e quelle calcolate da EUCENTRE dall’ipotesi di terremoto di Mw = 6.5 lungo la faglia
CS_01, non sono sufficienti a destabilizzare il corpo di frana. Ciò significa che l’innesco della frana potrebbe
risiedere in un’altra causa, che non sia quella sismica. Per la frana CRS invece non si ha instabilità per le PGA
identificate nel progetto, la sua mobilizzazione potrebbe avvenire per un carico sismico maggiore di 6.5.
Linea di ricerca 4. Modellazione dell’onda di maremoto potenzialmente prodotta dalle frane sottomarine
L’approccio utilizzato per sviluppare gli scenari di maremoto prevede differenti fasi, che comprendono la generazione della frana, la sua evoluzione, la generazione dell’onda di tsunami o maremoto e la propagazione del maremoto stesso. Nei casi in esame, in primo luogo si è proceduto alla ricostruzione del corpo di frana e della morfologia pre-evento, utilizzando le evidenze batimetriche (nicchie di frana) e i depositi lungo scarpata. Si è poi studiata la stabilità del pendio mediante una variante del metodo dell’Equilibrio Limite (Paparo e Tinti, 2017) e si è ricavata la soglia di PGA sufficiente ad attivare la frana. Nella fase successiva è stato simulato il moto della frana mediante il modello numerico UBO-BLOCK1. Il codice UBO-BLOCK1 calcola il moto della frana, vista come un insieme di blocchi interagenti che si possono deformare conservando tuttavia il volume complessivo (Tinti et al., 1997). Infine, la propagazione del maremoto è stata calcolata con il modello numerico UBO-TSUFD (metodo alle differenze finite per la risoluzione delle equazioni di Navier-Stokes non lineari in approssimazione di acqua bassa (Tinti e Tonini, 2013), che consente la produzione di mappe di pericolosità sulla base delle altezze massime raggiunte dalle onde lungo la costa.
PROGETTO SPOT - Report integrato di fine progetto
23
La modellazione indica che un maremoto ipoteticamente generato dalla frana NTS si propaga e colpisce prima di tutto la costa, in prossimità della città di Santa Croce Camerina, dopo circa 12 minuti dall’innesco della frana; e che, dopo 30 minuti, tutta la costa del Bacino di Gela viene raggiunta dallo tsunami. I valori massimi di elevazione dell’onda a costa indicano che il Golfo di Gela, tra le località di Licata e Santa Croce Camerina, è caratterizzato da valori sempre superiori a 1 m, che si riducono velocemente fino a meno di 0.5 m nel tratto di costa a ovest di Licata, mentre rimangono costanti tra 0.5 m ed 1 m fino a Capo Passero (Figura L4-1A).
Figura L4-1 – A: massimi valori di elevazione del mare lungo la costa causati dalla messa in posto della frana NTS; e B: della frana SGBS. Le località sono proiettate nel grafico sottostante, dove i valori massimi sono mostrati in funzione della longitudine in coordinate UTM.
La frana SGBS genera un maremoto che colpisce dopo circa 10 minuti l’isola di Gozo nell’arcipelago di Malta, e arriva dopo 20 minuti sulle coste della Sicilia. Dopo circa 30 minuti, tutte le coste dell’arcipelago di Malta sono state colpite mentre le coste siciliane lo sono solo in minima parte. La costa tra Santa Croce Camerina e Pozzallo è quella maggiormente colpita, con valori superiori ai 4 m presso il primo dei due insediamenti (Figura L4-1B). Il Golfo di Gela presenta valori generalmente compresi nella fascia tra 1 e 2 m, con solo alcuni valori superiori, in particolare vicino a Licata.
Il maremoto generato dalla frana CRS (Figura L4-2) colpisce prima di tutto la costa di Capo Rizzuto con un primo fronte negativo e, dopo circa 12 minuti dall’innesco della frana, raggiunge Botricello e Crotone. Il grafico a destra nella FiguraL4-2 mostra i valori massimi di elevazione delle onde di maremoto lungo la costa, nel tratto di costa tra le città di Capo Rizzuto e Capo Colonna, i cui picchi raggiungono valori superiori a 25 m. In generale, un tratto di costa calabrese di oltre 250 km è mediamente interessato da onde con elevazione massima che si pone nell’intervallo tra 1 e 3 m.
A B
PROGETTO SPOT - Report integrato di fine progetto
24
Figura L4-2 – Massimi valori di elevazione del mare lungo la costa causati dalla messa in posto della frana CRS. Il grafico a destra mostra la curva con i valori delle massime elevazioni raggiunte lungo la costa.
Linea di ricerca 5. Modellazione dell’onda di maremoto potenzialmente prodotta dall’attivazione delle faglie di scenario
Modellazione dell’onda di maremoto
L’INGV ha effettuato simulazioni numeriche relative a 11 scenari di tsunami sismicamente indotti. La modellazione è stata eseguita utilizzando il codice multi-GPU Tsunami-HySEA sviluppato dal gruppo EDANYA (Università di Malaga), che risolve le equazioni di acque basse nella loro approssimazione non-lineare.
I dati di input richiesti dalla simulazione sono la griglia batimetrica e la condizione iniziale per la propagazione dello tsunami, cioè la componente verticale in mare della deformazione cosismica dello scenario.
I dati batimetrici provengono da EMODnet (EMODnet Bathymetry Consortium, 2018) e hanno una risoluzione spaziale di ~ 90 m. Le simulazioni sono state eseguite utilizzando la riflessione pura come condizione al contorno per la parte di griglia corrispondente alla terra; questa condizione viene implementata fissando tutti i valori topografici positivi a un valore significativamente più grande. Le simulazioni degli scenari sono state eseguite considerando due diverse condizioni per la profondità minima del mare (cioè 5 e 10 m).
Per ogni simulazione sono stati estratti i valori delle altezze massime lungo le coste italiane applicando la legge di Green:
𝜂1
𝜂𝐷= (
𝐻𝐷
𝐻1)
1 4⁄
dove 𝜂1 e 𝜂𝐷 sono l’altezza dell’onda di tsunami alle profondità 𝐻1 = 1 𝑚 e 𝐻𝐷 (5 o 10 m), rispettivamente.
I risultati ottenuti relativamente a questi scenari sono stati utilizzati dal gruppo di ricerca del Consorzio ReLUIS per la stima dell’impatto dell’onda sul costruito. Una stima accurata dell’impatto di uno tsunami necessiterebbe di modellare numericamente la fase di inondazione ad una scala spaziale dell’ordine dei
PROGETTO SPOT - Report integrato di fine progetto
25
metri (o poche decine di metri). In assenza di dati topo-batimetrici accurati e con risoluzione sufficiente, di concerto con i colleghi dell’Università di Bologna, si è deciso di applicare il metodo dell’Energy Grade Line (EGL; Kriebel et al., 2017) per la trasformazione dell’altezza dello tsunami sulla linea di costa (𝜂1) in altezza (approssimata) della colonna d’acqua alle varie posizioni necessarie per la stima dell’impatto dell’onda sul costruito.
Il metodo dell’EGL prevede l’estrapolazione dell’altezza della colonna d’acqua lungo una serie di profili topografici 1D (da Trieste a San Vito Lo Capo, con una spaziatura di circa 20 m), senza tenere conto, tra le altre cose, di eventuali effetti di amplificazione che potrebbero derivare dalla propagazione 2D. Per considerare in modo speditivo i possibili effetti di focalizzazione dell’energia dello tsunami a terra, i valori di 𝜂1 per questa analisi sono stati quindi ulteriormente amplificati, moltiplicandoli per un fattore 2 (es.: Figura L5-1).
Figura L5-1 – Deformazione cosismica del fondo del mare e altezze massime dello tsunami stimate lungo le coste italiane associate allo scenario CS_2 (Fig. L4-1, Tabella L2-2). La profondità minima del mare adottata in questo scenario è 5 m.
Linea di ricerca 6. Modellazione di scenari di danno su strutture e infrastrutture costiere, e conseguenti
perdite, dovute ai terremoti generati dalle faglie di scenario (incluso webGIS dedicato)
La fondazione EUCENTRE si è occupata di definire gli scenari di danno sulle strutture e infrastrutture della
costa derivanti dall’innesco di faglie nelle aree di estrazione. Gli scenari sono stati generati per aree della
costa che ricadono entro un raggio di 100 km dalle piattaforme produttive, così come riportato in Figura 1.
Il raggio di 100 km è considerato come il massimo raggio di impatto di terremoti anche di magnitudo
importante. Una volta identificate le strutture e infrastrutture per le quali calcolare lo scenario di danno
ricadenti all’interno dell’area indicata, ne è stata quantificata la vulnerabilità sismica, cioè la loro
propensione a danneggiarsi se soggette a scuotimento sismico. Tale obiettivo è stato raggiunto assegnando
ad ogni struttura/infrastruttura un set di curve di fragilità. Nota la vulnerabilità, è stato necessario
quantificare lo scenario di scuotimento rappresentativo per gli eventi innescati. Sono state poi selezionate
le relazioni di attenuazione del moto sismico da implementare ed è stato sviluppato un tool per il calcolo
dell’accelerazione al sito risentita dalle strutture/infrastrutture in seguito all’attivazione delle faglie
identificate da INGV. Infine, è stato creato un WebGIS che consente di georeferenziare i dati delle
strutture/infrastrutture esposte, le sorgenti degli eventi e visualizzare su mappa e in tabelle gli scenari di
PROGETTO SPOT - Report integrato di fine progetto
26
danno calcolati. In Figura L6-1 è mostrato lo schema seguito per il calcolo dello scenario di danno e per la
pubblicazione dei dati di input e di output sul WebGIS.
Figura L6-1 – Schema per il calcolo dello scenario di danno.
Descrizione della base dati
La base dati è composta dalle strutture e infrastrutture identificate all’interno dell’area di studio. Sono stati
presi in considerazione innanzi tutto gli edifici residenziali presenti in 1563 comuni. Le informazioni su
numero di edifici, numero di abitazioni, popolazione e superficie abitativa per ogni comune sono fornite dal
database ISTAT 2001. In secondo luogo, sono state prese in considerazione le scuole. Sono stati identificati
11167 edifici scolastici dal Censimento Anagrafe dell’Edilizia Scolastica reperito tramite il MIUR; per 418
scuole sono disponibili anche le schede di verifica sismica di livello 1 e 2. Inoltre, per quanto riguarda le
infrastrutture della rete viabilistica, sono stati presi in considerazione 4837 ponti, 11514 opere di sostegno
e 286 gallerie. Le informazioni relative ai ponti derivano da ANAS e dalle Schede di verifica sismica di livello
1 e 2 disponibili presso il DPC. Il database sulle opere di sostegno contiene dati parziali di geometria
derivanti dalla scheda di accatastamento ANAS. Per le gallerie è nota da ANAS e Autostrade per l’Italia la
georeferenziazione e alcune caratteristiche geometriche, mentre la profondità è stata stabilita utilizzando il
DEM (Digital Elevation Model). Infine le informazioni relative alle opere portuali (banchine, gru, ecc.)
derivano dal MIT e dalle autorità portuali. I porti presi in esame sono 55, di cui 47 con livello di conoscenza
limitato e 8 per i quali sono disponibili maggiori informazioni.
Curve di fragilità adottate
Le curve di fragilità sono relazioni che indicano il danno atteso in una costruzione in funzione di un
parametro di intensità sismica (accelerazione di picco, intensità spettrale, intensità macrosismiche), ovvero
sono rappresentative della probabilità di superamento di un certo livello di danno al variare dell’input
sismico (Whitman et al., 1973). In ascissa esse riportano l’accelerazione di picco al suolo (PGA) al sito in cui
si trova la struttura, e in ordinata la probabilità di raggiungere le condizioni che identificano il livello di
danno, ad esempio il collasso. Sono state pertanto definite le curve di fragilità da assegnare alle opere
prese in considerazione; per alcune opere le curve di fragilità sono state selezionate tra quelle presenti in
letteratura, per altre sono state sviluppate dai ricercatori di EUCENTRE. In Tabella L6-1 sono riportate le
curve adottate e i relativi livelli di danno.
PROGETTO SPOT - Report integrato di fine progetto
27
Tabella L6-1 – Curve di fragilità implementate.
Opera Livelli di danno Origine delle curve
Edilizia residenziale D1, D2, D3, D4, D5 Curve sviluppate da EUCENTRE con la metodologia SP-BELA (Simplified Pushover-Based Earthquake Loss Assessment; Borzi et al., 2008a, 2008b)
Scuole D1, D2, D3, D4, D5 Curve sviluppate da EUCENTRE con la metodologia SP-BELA (Simplified Pushover-Based Earthquake Loss Assessment; Borzi et al., 2008a, 2008b)
Ponti Danno, collasso Curve sviluppate da EUCENTRE in funzione del livello di dettaglio delle informazioni disponibili