ricostruzione tridimensionale delle caratteristiche idrogeologiche ...

6 Ricostruzione 3D delle caratteristiche idrogeologiche nel bacino tra Ticino e Oglio
6.2 Modellazione 3D del sistema idrogeologico
In questa fase del lavoro verrà simulato il volume dell’intero
sistema idrogeologico in esame attraverso una griglia di calcolo 3D:
questo è un passaggio fondamentale nell’ottica della realizzazione di
un modello di flusso successivo e in relazione al processo di
feedback iterativo descritto nel capitolo 3. Infatti, nel mondo di
Modflow (software per il calcolo del flusso nel saturo, paragrafo 9.4)
l’estensione e la discretizzazione della griglia, la stratificazione del
sistema e delle proprietà entro e attraverso gli strati sono elementi
fondamentali per la costruzione del modello di flusso. In particolare la
stratificazione del sistema è importante perché all’interno di un
singolo strato il movimento dell’acqua è solamente orizzontale; il
movimento verticale, invece, viene simulato attraverso l’inserimento
di molteplici strati ed è così possibile distinguere i diversi gradienti
verticali (Feinstein, 2010).
Viene brevemente descritto il procedimento generale di
costruzione di una griglia tramite Gocad.
La griglia (l’oggetto modellistico tridimensionale è denominato
SGrid) viene creata in base alle dimensioni (estensione e quote) delle
superfici limite del sistema che deve riprodurre: la griglia è
inizialmente regolare, ossia tutte le celle in cui essa è discretizzata
hanno le medesime dimensioni. E’ di tipo cell-centered (Figura 6.13),
ovvero il valore delle proprietà viene assegnato a centro cella.
L’origine dell’SGrid (X0, Y0, Z0) si trova in basso a sinistra; gli assi X,
Y e Z sono denominati rispettivamente u, v e k e sono definiti da tre
corrispondenti vettori step_v, step_u e step_w. Ogni asse è diviso in
nodi, corrispondenti ai vertici di ciascuna cella, i cui indici sono i lungo
l’asse X, j lungo l’asse Y e k lungo l’asse Z. La somma dei nodi lungo
ogni asse è definita dai valori nv, nu e nw.
61