La valutazione del rischio per l’ambiente per i depositi di idrocarburi liquidi: un caso di studio

Il presente elaborato di tesi si pone come obiettivo quello di valutare il rischio per l’ambiente causato dal rilascio accidentale di idrocarburi liquidi; in particolare è stato preso a riferimento un deposito reale contenente tali sostanze e si è valutato il rischio ambientale da esso generato nelle matrici del suolo e sottosuolo e delle acque superficiali. Per far ciò sono stati applicati il metodo speditivo descritto nel Rapporto 57/2005 (Rapporto APAT) elaborato dal Gruppo misto APAT/ARPA/CNVVF, che ha provveduto all'individuazione di una metodologia speditiva per la valutazione del rischio per l’ambiente da incidenti rilevanti in depositi di idrocarburi liquidi con riferimento a suolo e sottosuolo, e il metodo speditivo per le acque superficiali descritto nel Rapporto ISPRA 92/2013, che riporta criteri e indirizzi tecnico-operativi per la valutazione delle analisi degli incidenti rilevanti con conseguenze per l’ambiente. Il Rapporto APAT propone una metodologia speditiva ad indici; per la sua applicazione è necessario calcolare 2 indici: l’Indice di Propensione al Rilascio, che si ottiene dalla conoscenza delle caratteristiche impiantistiche del deposito preso in esame e l’Indice di Propensione alla Propagazione, che si ottiene dalla conoscenza delle caratteristiche idrogeologiche del sito su cui insiste il deposito. Il Rapporto ISPRA presenta invece diverse metodologie per la valutazione del rischio ambientale maturate in ambito europeo e propone un nuovo metodo a indici per la valutazione del rischio ambientale per le acque superficiali. Per l’applicazione di tale metodo a indici si è calcolato dapprima l’Indice di Propensione al Rilascio, così come richiesto anche nel metodo speditivo per suolo e sottosuolo e si è successivamente classificato il corpo idrico su cui insiste il deposito. Attraverso l’applicazione dei 2 metodi è stato dunque possibile valutare il livello di criticità ambientale del deposito, determinando la necessità di procedere con una metodologia semplificata o una metodologia avanzata per la stima del rischio di contaminazione ambientale.

Studio con HEC-RAS di alcuni interventi per miglioramenti della funzionalità idraulica nel territorio del Consorzio di Bonifica della Romagna

Nel lavoro di tesi vengono esposti gli studi svolti su due aree situate nella giurisdizione riminese del Consorzio di Bonifica, riportanti forti problematiche idrauliche nelle reti di scolo: il territorio dei canali Tamagnino e Zonara Masere ed il Rio Gessi. La modellazione dei tratti è avvenuta attraverso i software Hec-Ras ed Hec-Hms, realizzati rispettivamente per l’analisi dei corsi d’acqua nelle pianure alluvionali attraverso la modellazione dell’evoluzione dei livelli e della corrente e per l’analisi dell’idrologia delle superfici di bacino. Il primo studio è dovuto a problematiche verificatesi fin dal 2014 causate da una variazione d’uso del suolo, il secondo studio è dovuto a dissesti idraulici iniziati nel 2013 causati alla presenza di un tombinamento che provoca l’ostruzione del tratto verso valle. Si esegue la modellazione dello stato di fatto per identificare le sezioni ed i manufatti problematici, in secondo luogo si realizzano le modellazioni delle ipotesi di progetto. Si definiscono infine le soluzioni praticabili e la spesa ipotizzata a seguito della realizzazione di un computo metrico estimativo.

Il terremoto del 24 agosto 2014 a Napa (California, USA): rilievo degli effetti di sito

Il lavoro di questa tesi è frutto dell’esperienza maturata in quattro mesi di tirocinio presso l’U.S. Geological Survey di Menlo Park, da ottobre 2015 a gennaio 2016. Durante questo periodo sono state eseguite quattro campagne di raccolta dati sismici e geotecnici presso la città di Napa, ove, il 24 agosto 2014, si è verificato un terremoto di magnitudo momento pari a 6.0. Quest’area è stata interessata in precedenza da terremoti minori, il più significativo dei quali risulta essere il terremoto di Yountville del 2000 di magnitudo momento pari a 5.1. I rilievi macrosismici, effettuati immediatamente dopo il terremoto di Napa del 2014, hanno fornito una classificazione dettagliata delle strutture danneggiate nel distretto colpito. Obiettivo di questa tesi è comprendere se la distribuzione dei danni osservata sia legata anche alle caratteristiche geologiche locali (effetti di sito), oltre che alla vicinanza all’epicentro. A questo scopo sono state acquisite 63 misure sismiche passive a stazione singola e 4 prove di sismica attiva multicanale, la cui inversione congiunta ha permesso di creare mappe dei riflettori sismici principali sotto le zone interessate dal terremoto. Allo stesso tempo tali mappe, interpretate alla luce della cartografia geologica, hanno permesso di costruire sezioni geologiche e di identificare le corrispondenze tra riflettori geologici e riflettori sismici. Si è così potuto osservare che le zone più pesantemente danneggiate dal terremoto sono quelle che ricadono ove si sono misurate amplificazioni per risonanza in medio-alta frequenza (> 3 Hz) mentre scarsi danni si sono registrati nelle aree caratterizzate da roccia affiorante (e curve di amplificazione per risonanza piatte) e danni minori nelle zone contraddistinte solo da bedrock sismico profondo. La distribuzione del danno, evidenziatosi prevalentemente dove le frequenze proprie di vibrazione del suolo sono mediamente sopra 3 Hz, risulta compatibile con gli intervalli di frequenza attesi per le strutture residenziali tipiche di Napa Valley, che sono costituite da 2-3 piani fuori terra, in legno e di altezza non superiore ai 10 metri.

Modellazione bidimensionale della propagazione nell'entroterra della inondazione da mare negli scenari dei previsti cambiamenti climatici

Le ricerche di carattere eustatico, mareografico, climatico, archeologico e geocronologico, sviluppatesi soprattutto nell’ultimo ventennio, hanno messo in evidenza che gran parte delle piane costiere italiane risulta soggetta al rischio di allagamento per ingressione marina dovuta alla risalita relativa del livello medio del mare. Tale rischio è la conseguenza dell’interazione tra la presenza di elementi antropici e fenomeni di diversa natura, spesso difficilmente discriminabili e quantificabili, caratterizzati da magnitudo e velocità molto diverse tra loro. Tra le cause preponderanti che determinano l’ingressione marina possono essere individuati alcuni fenomeni naturali, climatici e geologici, i quali risultano fortemente influenzati dalle attività umane soprattutto a partire dal XX secolo. Tra questi si individuano: - la risalita del livello del mare, principalmente come conseguenza del superamento dell’ultimo acme glaciale e dello scioglimento delle grandi calotte continentali; - la subsidenza. Vaste porzioni delle piane costiere italiane risultano soggette a fenomeni di subsidenza. In certe zone questa assume proporzioni notevoli: per la fascia costiera emiliano-romagnola si registrano ratei compresi tra 1 e 3 cm/anno. Tale subsidenza è spesso il risultato della sovrapposizione tra fenomeni naturali (neotettonica, costipamento di sedimenti, ecc.) e fenomeni indotti dall’uomo (emungimenti delle falde idriche, sfruttamento di giacimenti metaniferi, escavazione di materiali per l’edilizia, ecc.); - terreni ad elevato contenuto organico: la presenza di depositi fortemente costipabili può causare la depressione del piano di campagna come conseguenza di abbassamenti del livello della falda superficiale (per drenaggi, opere di bonifica, emungimenti), dello sviluppo dei processi di ossidazione e decomposizione nei terreni stessi, del costipamento di questi sotto il proprio peso, della carenza di nuovi apporti solidi conseguente alla diminuita frequenza delle esondazioni dei corsi d’acqua; - morfologia: tra i fattori di rischio rientra l’assetto morfologico della piana e, in particolare il tipo di costa (lidi, spiagge, cordoni dunari in smantellamento, ecc. ), la presenza di aree depresse o comunque vicine al livello del mare (fino a 1-2 m s.l.m.), le caratteristiche dei fondali antistanti (batimetria, profilo trasversale, granulometria dei sedimenti, barre sommerse, assenza di barriere biologiche, ecc.); - stato della linea di costa in termini di processi erosivi dovuti ad attività umane (urbanizzazione del litorale, prelievo inerti, costruzione di barriere, ecc.) o alle dinamiche idro-sedimentarie naturali cui risulta soggetta (correnti litoranee, apporti di materiale, ecc. ). Scopo del presente studio è quello di valutare la probabilità di ingressione del mare nel tratto costiero emiliano-romagnolo del Lido delle Nazioni, la velocità di propagazione del fronte d’onda, facendo riferimento allo schema idraulico del crollo di una diga su letto asciutto (problema di Riemann) basato sul metodo delle caratteristiche, e di modellare la propagazione dell’inondazione nell’entroterra, conseguente all’innalzamento del medio mare . Per simulare tale processo è stato utilizzato il complesso codice di calcolo bidimensionale Mike 21. La fase iniziale di tale lavoro ha comportato la raccolta ed elaborazione mediante sistema Arcgis dei dati LIDAR ed idrografici multibeam , grazie ai quali si è provveduto a ricostruire la topo-batimetria di dettaglio della zona esaminata. Nel primo capitolo è stato sviluppato il problema del cambiamento climatico globale in atto e della conseguente variazione del livello marino che, secondo quanto riportato dall’IPCC nel rapporto del 2007, dovrebbe aumentare al 2100 mediamente tra i 28 ed i 43 cm. Nel secondo e terzo capitolo è stata effettuata un’analisi bibliografica delle metodologie per la modellazione della propagazione delle onde a fronte ripido con particolare attenzione ai fenomeni di breaching delle difese rigide ed ambientali. Sono state studiate le fenomenologie che possono inficiare la stabilità dei rilevati arginali, realizzati sia in corrispondenza dei corsi d’acqua, sia in corrispondenza del mare, a discapito della protezione idraulica del territorio ovvero dell’incolumità fisica dell’uomo e dei territori in cui esso vive e produce. In un rilevato arginale, quale che sia la causa innescante la formazione di breccia, la generazione di un’onda di piena conseguente la rottura è sempre determinata da un’azione erosiva (seepage o overtopping) esercitata dall’acqua sui materiali sciolti costituenti il corpo del rilevato. Perciò gran parte dello studio in materia di brecce arginali è incentrato sulla ricostruzione di siffatti eventi di rottura. Nel quarto capitolo è stata calcolata la probabilità, in 5 anni, di avere un allagamento nella zona di interesse e la velocità di propagazione del fronte d’onda. Inoltre è stata effettuata un’analisi delle condizioni meteo marine attuali (clima ondoso, livelli del mare e correnti) al largo della costa emiliano-romagnola, le cui problematiche e linee di intervento per la difesa sono descritte nel quinto capitolo, con particolare riferimento alla costa ferrarese, oggetto negli ultimi anni di continui interventi antropici. Introdotto il sistema Gis e le sue caratteristiche, si è passati a descrivere le varie fasi che hanno permesso di avere in output il file delle coordinate x, y, z dei punti significativi della costa, indispensabili al fine della simulazione Mike 21, le cui proprietà sono sviluppate nel sesto capitolo.