Modellazione numerica del comportamento di una struttura in muratura

In questi ultimi anni il tema della sicurezza sismica degli edifici storici in muratura ha assunto particolare rilievo in quanto a partire soprattutto dall’ordinanza 3274 del 2003, emanata in seguito al sisma che colpì il Molise nel 2002, la normativa ha imposto un monitoraggio ed una classificazione degli edifici storici sotto tutela per quanto riguarda la vulnerabilità sismica (nel 2008, quest’anno, scade il termine per attuare quest’opera di classificazione). Si è posto per questo in modo più urgente il problema dello studio del comportamento degli edifici storici (non solo quelli che costituiscono monumento, ma anche e soprattutto quelli minori) e della loro sicurezza. Le Linee Guida di applicazione dell’Ordinanza 3274 nascono con l’intento di fornire strumenti e metodologie semplici ed efficaci per affrontare questo studio nei tempi previsti. Il problema si pone in modo particolare per le chiese, presenti in grande quantità sul territorio italiano e di cui costituiscono gran parte del patrimonio culturale; questi edifici, composti di solito da grandi elementi murari, non presentano comportamento scatolare, mancando orizzontamenti, elementi di collegamento efficace e muri di spina interni e sono particolarmente vulnerabili ad azioni sismiche; presentano inoltre un comportamento strutturale a sollecitazioni orizzontali che non può essere colto con un approccio globale basato, ad esempio, su un’analisi modale lineare: non ci sono modi di vibrare che coinvolgano una sufficiente parte di massa della struttura; si hanno valori dei coefficienti di partecipazione dei varii modi di vibrare minori del 10% (in generale molto più bassi). Per questo motivo l’esperienza e l’osservazione di casi reali suggeriscono un approccio di studio degli edifici storici sacri in muratura attraverso l’analisi della sicurezza sismica dei cosiddetti “macroelementi” in cui si può suddividere un edificio murario, i quali sono elementi che presentano un comportamento strutturale autonomo. Questo lavoro si inserisce in uno studio più ampio iniziato con una tesi di laurea dal titolo “Analisi Limite di Strutture in Muratura. Teoria e Applicazione all'Arco Trionfale” (M. Temprati), che ha studiato il comportamento dell’arco trionfale della chiesa collegiata di Santa Maria del Borgo a San Nicandro Garganico (FG). Suddividere un edificio in muratura in più elementi è il metodo proposto nelle Linee Guida, di cui si parla nel primo capitolo del presente lavoro: la vulnerabilità delle strutture può essere studiata tramite il moltiplicatore di collasso quale parametro in grado di esprimere il livello di sicurezza sismica. Nel secondo capitolo si illustra il calcolo degli indici di vulnerabilità e delle accelerazioni di danno per la chiesa di Santa Maria del Borgo, attraverso la compilazione delle schede dette “di II livello”, secondo quanto indicato nelle Linee Guida. Nel terzo capitolo viene riportato il calcolo del moltiplicatore di collasso a ribaltamento della facciata della chiesa. Su questo elemento si è incentrata l’attenzione nel presente lavoro. A causa della complessità dello schema strutturale della facciata connessa ad altri elementi dell’edificio, si è fatto uso del codice di calcolo agli elementi finiti ABAQUS. Della modellazione del materiale e del settaggio dei parametri del software si è discusso nel quarto capitolo. Nel quinto capitolo si illustra l’analisi condotta tramite ABAQUS sullo stesso schema della facciata utilizzato per il calcolo manuale nel capitolo tre: l’utilizzo combinato dell’analisi cinematica e del metodo agli elementi finiti permette per esempi semplici di convalidare i risultati ottenibili con un’analisi non-lineare agli elementi finiti e di estenderne la validità a schemi più completi e più complessi. Nel sesto capitolo infatti si riportano i risultati delle analisi condotte con ABAQUS su schemi strutturali in cui si considerano anche gli elementi connessi alla facciata. Si riesce in questo modo ad individuare con chiarezza il meccanismo di collasso di più facile attivazione per la facciata e a trarre importanti informazioni sul comportamento strutturale delle varie parti, anche in vista di un intervento di ristrutturazione e miglioramento sismico.

Monitoraggio e analisi di segnali ORP, pH, DO finalizzato al controllo del processo di predenitrificazione eseguito su impianto pilota a F.A.

Il lavoro sperimentale condotto nell’ambito della presente tesi è stato svolto su un impianto pilota a flusso continuo con schema predenitro/nitro in scala di laboratorio, alimentato con refluo sintetico ed installato presso i laboratori della Sezione ACS PROT IDR - Gestione Risorse Idriche dell’ENEA (sede di Bologna). Le finalità della sperimentazione erano mirate ad una verifica preliminare sull’utilizzo di segnali indiretti, quali pH, ORP, DO, nel monitoraggio in tempo reale dei processi biologici, mirata alla possibilità di implementare un sistema di controllo automatico di un impianto di depurazione, capace di minimizzare i costi di gestione e in grado di garantire, in ogni caso, elevate efficienze depurative. Nel presente studio sono stati esaminati gli aspetti metodologici e tecnologici alla base dell’automazione di un impianto di trattamento di acque reflue, avendo il duplice obiettivo di garantire una buona qualità dell’effluente e un contestuale risparmio energetico. In particolare si è studiato un problema molto diffuso: il controllo della rimozione dell’azoto in un sistema a fanghi attivi con predenitrificazione per l’ossidazione congiunta di carbonio ed azoto, ponendolo in relazione con la variazione dei tempi di commutazione in un processo SBR. L’obiettivo primario della sperimentazione effettuata è stato quello di portare il sistema in un regime di equilibrio stazionario, in modo da poter definire delle condizioni di regime per le quali, ad un ingresso costante e noto, corrispondesse un’uscita altrettanto costante che mantenesse tali condizioni fino a quando non fosse cambiato l’ingresso e/o le condizioni al contorno. A questo scopo si è provveduto, in primo luogo, allo studio di un influente sintetico con cui alimentare l’impianto, di composizione tale da simulare in tutte le sue caratteristiche un refluo civile reale e, successivamente, alla valutazione di una configurazione impiantistica idonea per il raggiungimento dell’equilibrio dei processi. La prima parte della sperimentazione è stata dedicata all’avviamento dell’impianto, effettuando un monitoraggio continuo dell’intero sistema allo scopo di impostare i parametri operativi e in modo da ottenere condizioni stabili per effettuare le successive sperimentazioni.

Analisi termo-strutturale dell'impianto di scarico di un motoveicolo ad elevate prestazioni

Il lavoro svolto, in collaborazione con “Ducati Motor Holding”, si propone l’analisi, mediante l’applicazione di un software FEM (HyperWorks 9.0), della struttura rappresentata dai condotti di un girotubi di scarico motociclistico. Il primo obiettivo consiste nell’ottenimento di un modello agli elementi finiti della struttura. In questo senso si è anche compiuta un’indagine sperimentale per la caratterizzazione di un componente chiave del sistema (un supporto antivibrante in gomma). Le simulazioni di calcolo strutturale vogliono coinvolgere gli effetti che i carichi termici hanno nell’alterare le proprietà meccaniche dei materiali e nell’indurre sollecitazioni derivate dalle dilatazioni termiche. Si svolgono inoltre analisi di tipo dinamico per determinare i modi di vibrare della struttura e come le frequenze naturali vengano influenzate dagli effetti termici. Sulla base di ulteriori simulazioni sulla risposta in frequenza del sistema e di un riscontro sperimentale (prova accelerometrica), si è inoltre tentata una prima validazione del modello dal punto di vista dinamico. La finalità del lavoro risiede dunque nella determinazione delle criticità strutturali che il girotubi può manifestare durante il funzionamento, ma anche nell’esplorazione e applicazione delle potenzialità del software applicato, per l’analisi degli aspetti sopra esposti.

Sviluppo di un modello quantitativo di valutazione e valorizzazione del servizio di taglio materiali: il caso G.D

La necessità di sviluppare questo progetto è nata dal fatto che quando due aziende desiderano collaborare in maniera costruttiva con rapporti di fornitura di materiali o lavorazioni, bisogna evitare che si verifichino scostamenti nei rapporti a discapito di uno dei due contraenti. La coscienza dell'inefficienza di un rapporto squilibrato ha portato le due aziende a intavolare un progetto migliorativo basato sulla massima trasparenza e propositività. Quando G.D ha deciso di avvalersi della fornitura da parte di Transmec per quanto riguarda il servizio di taglio del materiale, l'ha fatto per andare a cercare un miglior livello di efficacia ed efficienza. Il raggiungimento degli attuali livelli di performance non sarebbe stato possibile se entrambe le aziende non avessero collaborato attivamente nell'ottenimento degli standard desiderati. Questo studio è nato dal fatto che il modello interno adottato precedentemente da G.D per la ripartizione dei costi delle Business Unit non era più efficace se utilizzato come strumento di pagamento. Si rendeva quindi necessaria un'analisi accurata della situazione produttiva utilizzando gli approcci analitici resi disponibili nel corso degli anni come lo strumento per l'allocazione dei costi secondo le attività (Activity Based Costing) e la mappatura dei processi utilizzata nei casi di re-ingegnerizzazione dei reparti.

Amenagement du territoire autour des sites Seveso: comparaison entre les approches en Italie et en region Wallonne

Il continuo verificarsi di gravi incidenti nei grandi impianti industriali ha spinto gli Stati membri dell’Unione Europea a dotarsi di una politica comune in materia di prevenzione dei grandi rischi industriali. Anche a seguito della pressione esercitata dall’opinione pubblica sono state implementate, nel corso degli ultimi quarant’anni, misure legislative sempre più efficaci per la prevenzione e la mitigazione dei rischi legati ad attività industriali particolarmente pericolose. A partire dagli ultimi anni dello scorso secolo, l’Unione Europea ha emanato una serie di direttive che obbligano gli Stati membri ad essere garanti della sicurezza per l’uomo e per l’ambiente nelle zone circostanti a stabilimenti a rischio di incidente rilevante. In quest’ottica è stata pubblicata nel 1982 la Direttiva Seveso I [82/501/EEC], che è stata ampliata nel 1996 dalla Direttiva Seveso II [96/82/CE] ed infine emendata nel dicembre 2003 dalla Direttiva Seveso III [2003/105/CE]. Le Direttive Seveso prevedono la realizzazione negli Stati membri di una valutazione dei rischi per gli stabilimenti industriali che sono suscettibili a incendi, esplosioni o rilasci di gas tossici (quali, ad esempio, le industrie chimiche, le raffinerie, i depositi di sostanze pericolose). La Direttiva Seveso II è stata trasposta in legge belga attraverso “l’Accord de Coopération” del 21 giugno 1999. Una legge federale nel giugno del 2001 [M.B. 16/06/2001] mette in vigore “l’Accord de Coopération”, che è stato in seguito emendato e pubblicato il 26 aprile del 2007 [M.B. 26/04/2007]. A livello della Regione Vallona (in Belgio), la tematica del rischio di incidente rilevante è stata inclusa nelle disposizioni decretali del Codice Vallone della Pianificazione Territoriale, dell’Urbanismo e del Patrimonio [CWATUP]. In questo quadro la Regione Vallona ha elaborato in collaborazione con la FPMs (Faculté Polytechnique de Mons) una dettagliata metodologia di analisi del rischio per gli stabilimenti a rischio di incidente rilevante. In Italia la Direttiva Seveso II è stata recepita dal Decreto Legislativo n°334 emanato nell’agosto del 1999 [D. Lgs. 334/99], che ha introdotto per la prima volta nel quadro normativo italiano i concetti fondamentali di “controllo dell’urbanizzazione” e “requisiti minimi di sicurezza per la pianificazione territoriale”. Il Decreto Legislativo 334/99 è attualmente in vigore, modificato ed integrato dal Decreto Legislativo n°238 del 21 settembre 2005 [D. Lgs. 238/05], recepimento italiano della Direttiva Seveso III. Tra i decreti attuativi del Decreto Legislativo 334/99 occorre citare il Decreto Ministeriale n°151 del 2001 [D. M. 151/01] relativo alla pianificazione territoriale nell’intorno degli stabilimenti a rischio di incidente rilevante. L’obiettivo di questo lavoro di tesi, che è stato sviluppato presso la Faculté Polytechnique di Mons, è quello di analizzare la metodologia di quantificazione del rischio adottata nella Regione Vallona, con riferimento alla pianificazione territoriale intorno agli stabilimenti a rischio di incidente rilevante, e di confrontarla con quella applicata in Italia. La metodologia applicata in Vallonia è di tipo “probabilistico” ovvero basata sul rischio quale funzione delle frequenze di accadimento e delle conseguenze degli scenari incidentali. Il metodo utilizzato in Italia è “ibrido”, ovvero considera sia le frequenze che le conseguenze degli scenari incidentali, ma non la loro ricomposizione all’interno di un indice di rischio. In seguito al confronto teorico delle due metodologie, se ne è effettuato anche una comparazione pratica tramite la loro applicazione ad un deposito di GPL. Il confronto ha messo in luce come manchino, nella legislazione italiana relativa agli stabilimenti a rischio di incidente rilevante, indicazioni di dettaglio per la quantificazione del rischio, a differenza di quanto accade nella legislazione belga. Ciò lascia all’analista di rischio italiano una notevole arbitrarietà nell’effettuare ipotesi ed assunzioni che rendono poi difficile la comparazione del rischio di stabilimenti differenti. L’auspicio è che tale lacuna possa essere rapidamente superata.

Approccio modellistico al processo di selezione e separazione dimensionale dei RSU ed applicazione ad un impianto TMB

Nel presente lavoro sono stati illustrati i risultati preliminari di uno studio mirato alla modellazione del processo di vagliatura di RSU. Il modello, sviluppato sulla base di ipotesi sia cinematiche che probabilistiche, è stato implementato ed applicato per la verifica dell’influenza delle condizioni operative sull’efficienza del processo di separazione. La modellazione è stata sviluppata studiando la cinematica di una particella all’interno di un vaglio rotante, prima studiando il comportamento della singola particella e poi studiando il comportamento di più particelle. Lo sviluppo di tale modello, consente di determinare l’efficienza di rimozione per le diverse frazioni merceologiche di un rifiuto, ciascuna caratterizzata da una propria distribuzione dimensionale,. La validazione è stata svolta effettuando una campagna di misurazioni utilizzando un vaglio in scala (1:10) e un vaglio di dimensioni reali . Da queste misurazioni è emerso che il modello è ben rappresentativo della realtà, in alcuni casi, si sono evidenziate delle discrepanze dal comportamento reale del rifiuto, ciò è senz’altro da attribuite alla forte eterogeneità dei materiali costituenti il rifiuto solido urbano e dall’interferenza che le une causano sulle altre. Si è appurato che, variando i parametri quali: l’angolo di inclinazione longitudinale del vaglio, la velocità di rotazione si ha un diverso comportamento del rifiuto all’interno del vaglio e quindi una variazione dell’efficienza di separazione di quest’ultimo. In particolare è stata mostrata una possibile applicazione del modello legata alla valutazione delle caratteristiche del sopravaglio nel caso questo sia destinato alla produzione di CDR. Attraverso la formula di Dulong si è calcolato il potere calorifico superiore, (sia allo stato umido che allo stato secco), ottenuto effettuando delle variazioni dei parametri citati precedentemente, calcolato sulla massa uscente dal sopravaglio e quindi quella utilizzata come combustibile per gli impianti di termovalorizzazione.( C.D.R.) Si è osservato il legame esistente tra le condizioni operative inclinazione e velocità di rotazione del tamburo ed il potere calorifico del flusso selezionato con l’obiettivo di massimizzare quest’ultimo. Dalle è prove è emerso come il fattore umidità giochi un ruolo fondamentale per l’ottimizzazione del CDR, l’umidità infatti incide negativamente sul potere calorifico del rifiuto invertendo il rapporto normalmente direttamente proporzionale tra percentuale di sopravaglio e potere calorifico, dunque a causa dell’umidità all’aumentare del flusso di sopravaglio si ha una diminuzione del potere calorifico. Lo sviluppo e la validazione di un modello teorico, ci permette di affrontare il problema della selezione e separazione dei RSU,non più solo tramite un procedimento induttivo basato unicamente sulle osservazioni, ma tramite un procedimento deduttivo grazie al quale "conoscendo l’input si prevede quale sarà l’output". Allo stato attuale la scelta e il dimensionamento delle unità di separazione dimensionale nei processi di selezione è basata unicamente su conoscenze di natura empirica L’obiettivo principale è stato dunque quello di costruire uno strumento in grado di determinare le caratteristiche dei flussi uscenti dall’unità di vagliatura, note che siano le caratteristiche del flusso in ingresso ed i parametri operativi Questo approccio permette di eseguire valutazioni oggettive circa la reale rispondenza delle caratteristiche dei flussi in uscita a standard prefissati,peraltro consente di eseguire valutazioni sulle uscite di processo rispetto alla variabilità delle caratteristiche dei rifiuti in ingresso .

Modellazioni FEM dinamiche per analisi di avarie dovute a vibrazioni su componenti di motociclo

Il lavoro svolto, in collaborazione con “Ducati Motor Holding”, si propone l’analisi, mediante l’applicazione di un software FEM (HyperWorks 9.0), delle strutture rappresentate dall’assieme paracoppa-parasassi e dal silenziatore. Il primo obiettivo consiste nell’ottenimento di un modello agli elementi finiti della struttura. Le simulazioni di calcolo consistono in analisi di tipo dinamico per determinare i modi di vibrare delle strutture. Sulla base di ulteriori simulazioni sulla risposta al transitorio del sistema si è studiato il loro comportamento nelle condizioni di funzionamento. La finalità del lavoro risiede dunque nella determinazione delle criticità strutturali che i componenti hanno manifestato durante prove sperimentali verificandole attraverso calcolo numerico e applicazione delle potenzialità del software applicato, per l’analisi degli aspetti sopra esposti.

La trasformata di Fourier nella valutazione d'opzioni

Lo scopo di questa tesi è di fornire una analisi sistematica delle condizioni necessarie per l’esistenza delle formule di valutazione di opzioni che impiegano la trasformata di Fourier. Perchè la trasformata di Fourier? Se assumiamo che il processo del prezzo del sottostante sia un processo di Lévy noi non conosciamo solitamente la sua Legge di probabilità ma sempre la sua funzione caratteristica che è proprio la trasformata di Fourier della Legge.

Valutazione di impatto ambientale di un gruppo elettrogeno mediante Life Cycle Assessment: il caso COGEM.

Il presente lavoro ha come oggetto l’analisi di impatto ambientale, svolta mediante la metodologia Life Cycle Assessment (LCA) di un gruppo elettrogeno prodotto da COGEM s.r.l., azienda italiana situata a Castel d’Argile, nel bolognese, con l’obiettivo di supportare eventuali scelte di riprogettazione del prodotto anche in termini di Design for Disassembly. Dopo una prima analisi del contesto attuale in cui si colloca, la metodologia LCA è stata studiata nel dettaglio per poterla poi applicare al prodotto in oggetto. Esso è stato individuato mediante un’analisi delle vendite di COGEM, in seguito si è svolta una fase di raccolta dati e si sviluppata l’analisi LCA usando il software SimaPro 7.1. I risultati ottenuti hanno consentito di individuare le possibili aree di miglioramento dell’impatto ambientale dell’intero ciclo di vita del gruppo elettrogeno. In particolare si sono valutate due soluzioni innovative: un gruppo elettrogeno alimentato a olio vegetale e uno progettato in ottica DFD per consentire un corretto smaltimento dei rifiuti.

La sensitività delle curve a fatica oligociclica rispetto alle colate di derivazione per i provini: applicazione alla caratterizzazione di un materiale per cappe di turboalternatore

In questo studio vengono riportati i risultati di prove di fatica oligociclica eseguiti su provini dello stesso materiale ottenuti con uguali processi tecnologici ma provenienti da differenti colate di metallo. Il materiale in questione è un acciaio di elevata qualità frequentemente utilizzato per la realizzazione di cappe per turboalternatori. Obiettivo dello studio è stato ricavare i coefficienti necessari per tracciare le curve di fatica del materiale, non ancora presenti in letteratura, ed infine indagare la bontà del risultato ottenuto con un’analisi statistica delle curve e dei risultati ottenuti. Nella prima parte è descritto l’attuale stato dell’arte e la situazione in cui si colloca il presente studio. Nella seconda parte viene fornita una descrizione dettagliata del materiale studiato, delle condizioni nelle quali sono state eseguite le prove e delle attrezzature utilizzate a tale scopo. Si conclude esponendo i risultati ottenuti, comprensivi dei confronti e delle considerazioni derivate dalle analisi statistiche eseguite.

Analisi dei processi di Supply Chain. Il problema dei mancanti: studio di due casi reali.

L'elaborato di tesi sviluppato ha come obiettivo l'ottimizzazione della Supply Chian di contesti manufatturieri con particolare attenzione alle procedure di approvvigionamento dei materiali e al problema dei materiali "mancanti" in produzione.Viene esposto il risultato di un’attività di analisi di due casi reali per l’ottimizzazione dell'impatto logistico del problema dell'alimentazione di sistemi di fabbricazione e montaggio. Lo studio è stato affrontato definendo un approccio basato su una duplice analisi del processo, top-down e bottom-up. Le due analisi condotte hanno permesso di procedere alla mappatura dell’intero processo end to end (che inizia con le richieste del cliente e termina con la consegna dei prodotti allo stesso) e di individuare le criticità in esso presenti, sulla base delle quali sono state fornite una serie di linee guida per la risoluzione del problema dei mancanti. I risultati ottenuti hanno evidenziato che l’ottimizzazione dell’area di pianificazione permette di migliorare la qualità dei dati in ingresso all’intero processo con conseguente diminuzione del numero di materiali mancanti. Tuttavia, la trasversalità dell’argomento ai processi aziendali impedisce la definizione di un approccio al problema focalizzato su un’area specifica. Infine, è stata eseguita una bozza di valutazione economica per la stima dei tempi di recupero degli investimenti necessari.

Valutazione delle proprietà di diverse misure d'intensità dello scuotimento sismico

La progettazione sismica negli ultimi anni ha subito una forte trasformazione, infatti con l’introduzione delle nuove normative tecniche si è passati dallo svolgere una verifica delle capacità locali dei singoli elementi ad una progettazione con un approccio di tipo probabilistico, il quale richiede il soddisfacimento di una serie di stati limite ai quali viene attribuita una certa probabilità di superamento. La valutazione dell’affidabilità sismica di una struttura viene condotta di solito attraverso metodologie che prendono il nome di Probabilistic Seismic Design Analysis (PSDA) in accordo con la procedura del Performance Based Earthquake Engineering (PBEE). In questa procedura di tipo probabilistico risulta di notevole importanza la definizione della misura d’intensità sismica, la quale può essere utilizzata sia come predittore della risposta strutturale a fronte di un evento sismico, sia come parametro per definire la pericolosità di un sito. Queste misure d’intensità possono essere definite direttamente dalla registrazione dell’evento sismico, come ad esempio l’accelerazione di picco del terreno, oppure sulla base della risposta, sia lineare che non, della struttura soggetta a tale evento, ovvero quelle che vengono chiamate misure d’intensità spettrali. Come vedremo è preferibile l’utilizzo di misure d’intensità che soddisfino certe proprietà, in modo da far risultare più efficace possibile le risoluzione del problema con l’approccio probabilistico PBEE. Obbiettivo principale di questa dissertazione è quello di valutare alcune di queste proprietà per un gran numero di misure d’intensità sismiche a partire dai risultati di risposta strutturale ottenuti mediante analisi dinamiche non lineari nel tempo, condotte per diverse tipologie di strutture con differenti proprietà meccaniche.

Valutazione del rischio di incidente rilevante con contaminazione del suolo: il caso di un oleodotto.

Il presente lavoro di tesi si colloca nell’ambito della valutazione del rischio di incidente rilevante. Ai sensi della normativa europea (direttive Seveso) e del loro recepimento nella legislazione nazionale (D.Lgs. 334/99 e s.m.i.) un incidente rilevante è costituito da un evento incidentale connesso al rilascio di sostanze pericolose in grado di causare rilevanti danni all’uomo e/o all’ambiente. Ora, se da un lato esistono indici di rischio quantitativi per il bersaglio ”uomo” da tempo definiti e universalmente adottati nonché metodologie standardizzate e condivise per il loro calcolo, dall’altro non vi sono analoghi indici di rischio per il bersaglio “ambiente” comunemente accettati né, conseguentemente, procedure per il loro calcolo. Mancano pertanto anche definizioni e metodologie di calcolo di indici di rischio complessivo, che tengano conto di entrambi i bersagli citati dalla normativa. Al fine di colmare questa lacuna metodologica, che di fatto non consente di dare pieno adempimento alle stesse disposizioni legislative, è stata sviluppata all’interno del Dipartimento di Ingegneria Chimica, Mineraria e delle Tecnologie Ambientali dell’Università degli Studi di Bologna una ricerca che ha portato alla definizione di indici di rischio per il bersaglio “ambiente” e alla messa a punto di una procedura per la loro stima. L’attenzione è stata rivolta in modo specifico al comparto ambientale del suolo e delle acque sotterranee (falda freatica) ed ai rilasci accidentali da condotte di sostanze idrocarburiche immiscibili e più leggere dell’acqua, ovvero alle sostanze cosiddette NAPL – Non Acqueous Phase Liquid, con proprietà di infiammabilità e tossicità. Nello specifico si sono definiti per il bersaglio “ambiente” un indice di rischio ambientale locale rappresentato, punto per punto lungo il percorso della condotta, dai volumi di suolo e di acqua contaminata, nonché indici di rischio ambientale sociale rappresentati da curve F/Vsuolo e F/Sacque, essendo F la frequenza con cui si hanno incidenti in grado di provocare contaminazioni di volumi di suolo e di superfici di falda uguali o superiori a Vsuolo e Sacque. Tramite i costi unitari di decontaminazione del suolo e delle acque gli indici di rischio ambientale sociale possono essere trasformati in indici di rischio ambientale sociale monetizzato, ovvero in curve F/Msuolo e F/Macque, essendo F la frequenza con cui si hanno incidenti in grado di provocare inquinamento di suolo e di acque la cui decontaminazione ha costi uguali o superiori a Msuolo ed Macque. Dalla combinazione delle curve F/Msuolo e F/Macque è possibile ottenere la curva F/Mambiente, che esprime la frequenza degli eventi incidentali in grado di causare un danno ambientale di costo uguale o superiore a Mambiente. Dalla curva di rischio sociale per l’uomo ovvero dalla curva F/Nmorti, essendo F la frequenza con cui si verificano incidenti in grado di determinare un numero di morti maggiore o uguale ad Nmorti, tramite il costo unitario della vita umana VSL (Value of a Statistical Life), è possibile ottenete la curva F/Mmorti, essendo F la frequenza con cui si verificano incidenti in grado di determinare un danno monetizzato all’uomo uguale o superiore ad Mmorti. Dalla combinazione delle curve F/Mambiente ed F/Mmorti è possibile ottenere un indice di rischio sociale complessivo F/Mtotale, essendo F la frequenza con cui si verifica un danno economico complessivo uguale o superiore ad Mtotale. La procedura ora descritta è stata implementata in un apposito software ad interfaccia GIS denominato TRAT-GIS 4.1, al fine di facilitare gli onerosi calcoli richiesti nella fase di ricomposizione degli indici di rischio. La metodologia è stata fino ad ora applicata ad alcuni semplici casi di studio fittizi di modeste dimensioni e, limitatamente al calcolo del rischio per il bersaglio “ambiente”, ad un solo caso reale comunque descritto in modo semplificato. Il presente lavoro di tesi rappresenta la sua prima applicazione ad un caso di studio reale, per il quale sono stati calcolati gli indici di rischio per l’uomo, per l’ambiente e complessivi. Tale caso di studio è costituito dalla condotta che si estende, su un tracciato di 124 km, da Porto Marghera (VE) a Mantova e che trasporta greggi petroliferi. La prima parte del lavoro di tesi è consistita nella raccolta e sistematizzazione dei dati necessari alla stima delle frequenze di accadimento e delle conseguenze per l’uomo e per l’ambiente degli eventi dannosi che dalla condotta possono avere origine. In una seconda fase si è proceduto al calcolo di tali frequenze e conseguenze. I dati reperiti hanno riguardato innanzitutto il sistema “condotta”, del quale sono stati reperiti da un lato dati costruttivi (quali il diametro, la profondità di interramento, la posizione delle valvole sezionamento) e operativi (quali la portata, il profilo di pressione, le caratteristiche del greggio), dall’altro informazioni relative alle misure di emergenza automatiche e procedurali in caso di rilascio, al fine di stimare le frequenze di accadimento ed i termini “sorgente” (ovvero le portate di rilascio) in caso di rotture accidentali per ogni punto della condotta. In considerazione delle particolarità della condotta in esame è stata sviluppata una procedura specifica per il calcolo dei termini sorgente, fortemente dipendenti dai tempi degli interventi di emergenza in caso di rilascio. Una ulteriore fase di raccolta e sistematizzazione dei dati ha riguardato le informazioni relative all’ambiente nel quale è posta la condotta. Ai fini del calcolo del rischio per il bersaglio “uomo” si sono elaborati i dati di densità abitativa nei 41 comuni attraversati dall’oleodotto. Il calcolo dell’estensione degli scenari incidentali dannosi per l’uomo è stato poi effettuato tramite il software commerciale PHAST. Allo scopo della stima del rischio per il bersaglio “ambiente” è stata invece effettuata la caratterizzazione tessiturale dei suoli sui quali corre l’oleodotto (tramite l’individuazione di 5 categorie di terreno caratterizzate da diversi parametri podologici) e la determinazione della profondità della falda freatica, al fine di poter calcolare l’estensione della contaminazione punto per punto lungo la condotta, effettuando in tal modo l’analisi delle conseguenze per gli scenari incidentali dannosi per l’ambiente. Tale calcolo è stato effettuato con il software HSSM - Hydrocarbon Spill Screening Model gratuitamente distribuito da US-EPA. La ricomposizione del rischio, basata sui risultati ottenuti con i software PHAST e HSSM, ha occupato la terza ed ultima fase del lavoro di tesi; essa è stata effettuata tramite il software TRAT-GIS 4.1, ottenendo in forma sia grafica che alfanumerica gli indici di rischio precedentemente definiti. L’applicazione della procedura di valutazione del rischio al caso dell’oleodotto ha dimostrato come sia possibile un’analisi quantificata del rischio per l’uomo, per l’ambiente e complessivo anche per complessi casi reali di grandi dimensioni. Gli indici rischio ottenuti consentono infatti di individuare i punti più critici della condotta e la procedura messa a punto per il loro calcolo permette di testare l’efficacia di misure preventive e protettive adottabili per la riduzione del rischio stesso, fornendo al tempo gli elementi per un’analisi costi/benefici connessa all’implementazione di tali misure. Lo studio effettuato per la condotta esaminata ha inoltre fornito suggerimenti per introdurre in alcuni punti della metodologia delle modifiche migliorative, nonché per facilitare l’analisi tramite il software TRAT-GIS 4.1 di casi di studio di grandi dimensioni.