Complex resistivity and spectral induced polarization techniques for environmental applications

The potential and limits of the complex resistivity and the Spectral Induced Polarization techniques for environmental applications are investigated starting from the controlled laboratory conditions. The solution of the 2D forward problem for the complex resistivity is achieved on cylindrical-shaped models by a transformed formulation and using the Complete Electrode Model to define boundary conditions. The finite element algorithm, implemented in Matlab to solve numerically these equations, has proven to be accurate and stable through a comparison with the 3D full formulation and synthetic simulations. Accordingly, this formulation can be used to produce 2D reconstruction on cylindrical laboratory devices provided with one horizontal mid-section of electrodes (e.g. lysimeters). The whole geophysical process (data processing, forward modelling and inversion) has been included within EIDORS through an interactive Matlab algorithm with numerous new function added. Simulations with Gaussian noise, performed on heterogeneous models with anomalies included, provide low-error inverted models whereas anomalies are well-detected both for shape and position. EIDORS software could act as a base code for tomographic inversion of time-domain and frequency-domain data, both for laboratory and for in field investigations, because of its high-flexibility and reliability now reached by forward and inversion routines, with the primary aim of to share data and promote collaboration between groups working these fields. A new experimental device was presented to perform 2D and 3D measurement on cylindrical physical models in order to simulate typical subsoil configurations. The cylindrical symmetry avoids to have boundary effects of prismatic tanks and allows us to provide horizontal cross-section within the different layers. Preliminary water and sand models have highlighted the main problems due to saturation, electrodes, array and frequency range. The experimental AC instrument, after the calibration, provides reliable amplitude and phase dataset with stable and repeatable measurements. Measurements are not affected by the inductive coupling for frequency up to 5-10 Hz. Physical models built up in the laboratory have concerned the preparation of a typical sand-clay-sand multilayer model, simulating a shallow aquifer by saturation of the top layer. The HFE 7100, having physical properties similar to the chlorinated solvents, was spilled from the top to simulate a point source spilling, that is for example a failure of a little storage tank. Laboratory experiments have confirmed the reliability of the algorithm proposed. After HFE spilling, impedance amplitude increases progressively as expected, taking into account the percentage ratio between the global volume of HFE spilled and the sand volume. On the other hand the phase increases less evidently, since they are probably more sensible to the particular working conditions. Evidences from SIP phase spectra seems to be also convincing as we detect a single peak for frequency of about 1-2 Hz. These results confirm that complex resistivity tomography can play an important for contaminant detection, stressing once more the importance to define codified criteria for acquiring CR ad SIP in field data.

La tesi concerne la valutazione delle potenzialità e dei limiti delle tecniche di Resistività Complessa e di Polarizzazione Indotta Spettrale per l'individuazione di contaminazione del sottosuolo, a partire dalle condizioni controllate di laboratorio. La soluzione del problema diretto 2D per indagini in laboratorio è stata realizzata tramite una formulazione trasformata e utilizzando il Complete Electrode Model per definire le condizioni al contorno. L'algoritmo agli elementi finiti, implementato in Matlab per risolvere numericamente queste equazioni, ha dimostrato di essere stabile e accurato attraverso un confronto con la formulazione 3D inclusa nel software EIDORS e per mezzo di simulazioni sintetiche con rumore gaussiano aggiunto. Di conseguenza, questa formulazione può essere utilizzata per produrre modelli 2D su dispositivi cilindrici di laboratorio provvisti di una sezione orizzontale in mezzeria (p.es. lisimetri). L'intero processing geofisico (elaborazione dei dati, modellazione e inversione) è stato incluso all'interno di EIDORS tramite un algoritmo in Matlab interattivo, per mezzo dell'aggiunta di numerose nuove funzioni. Simulazioni con rumore gaussiano, eseguite su modelli eterogenei con anomalie in ampiezza e fase, hanno fornito dei modelli finali con valori limitati di errore e una buona ricostruzione delle anomalie sia in forma che in posizione. Pertanto il software EIDORS può costituire un codice di base per l'inversione tomografica di dati nel dominio del tempo e nel dominio della frequenza, sia per indagini in laboratorio che in sito, a causa della sua flessibilità e dell'affidabilità raggiunta dalla modellazione e dall'inversione, con l'obiettivo primario di condividere i dati e promuovere la collaborazione tra i diversi gruppi di ricerca. Per i fini suddetti è stato progettato un nuovo dispositivo sperimentale per eseguire misure 2D e 3D su modelli fisici di laboratorio, capaci di simulare le configurazioni tipiche del sottosuolo. La simmetria cilindrica permette di minimizzare gli effetti di bordo e di ricavare delle sezioni orizzontali all'interno dei singoli strati. Modelli preliminari di acqua e sabbia hanno messo in luce i principali problemi nelle campagne di misura provocati dalla saturazione, dal materiale costituente gli elettrodi, dagli array e dal range di frequenza. La strumentazione sperimentale in corrente alterenata, dopo la calibrazione e la minimizzazione dei suddetti errori, fornisce misure di ampiezza e fase affidabili, stabili e ripetibili. Le misurazioni non sono affette dall'accoppiamento induttivo per frequenze fino a 5-10 Hz. I modelli fisici costruiti in laboratorio hanno riguardato la preparazione di una tipica alternanza sabbia-argilla-sabbia, simulando un acquifero superficiale per saturazione del livello superiore. Lo sversamento di HFE 7100, avente caratteristiche fisiche simili ai solventi clorurati, è stato condotto dall'alto, ipotizzando la rottura di un piccolo serbatoio di stoccaggio. Gli esperimenti di laboratorio hanno confermato l'affidabilità dell'algoritmo proposto. Dopo lo sversamento di HFE, l'ampiezza dell'impedenza aumenta progressivamente come previsto, tenendo conto del rapporto percentuale tra il volume globale di HFE sversato e il volume di sabbia. D'altra parte la fase evidenzia un aumento meno marcato ed è probabilmente più sensibile alle particolari condizioni di lavoro e al rapporto volumetrico tra contaminante e materiale contaminato. Inoltre negli spettri di fase si individua un unico picco per frequenze di circa 1-2 Hz. Questi risultati confermano l'affidabilità della tecnica tomografica a resistività complessa per l'individuazione di contaminanti, sottolineando ancora una volta l'importanza di definire dei criteri codificati per l'acquisizione in sito.