MODFLOW – Discretizzazione verticale

Per i due problemi che poni puoi provare a seguire queste due strade
Il primo messaggio dice che 52190 celle del layer 1 (il layer superiore) hanno un carico idraulico di partenza inferiore alla quota del bottom del layer. In pratica quindi esistono 52190 celle che all’inizio della simulazione risultano insature. Questo fenomeno di solito capita quando si modificano il top e il bottom dei layer per adattarli alla topografia, operazione che tu hai dovuto fare per adattare il modello alla topografia reale del sito.
In sintesi succede questo: all’inizio, nella finestra di immissione che hai in fase di creazione di un nuovo modello, ti viene chiesto di definire il top e il bottom del modello e il numero del layer. Se ad esempio decidi che il top è 30 m, il bottom è 0 m e che il modello è suddiviso in 3 layer, il primo layer sarà di default compreso fra 20 e 30 m.
Il Groundwater Vistas assume sempre di default che il carico idraulico iniziale da attribuire alle celle del modello sia uguale al top che hai definito in fase di immissione, cioè nel caso dell’esempio 30 m. Se tu vai a modificare la topografia del modello è possibile che dopo aver introdotto il top e il bottom dei vari layer, il bottom del primo layer sia in aluni punti superiore ai 30 m e di conseguenza una parte delle celle all’inizio della simulazione risultano insature.
Il Modflow in generale si comporta bene quando una cella nel corso della simulazione da satura diviene insatura ma può avere difficoltà nel processo inverso, cioè nel rendere sature delle celle insature. Per questa ragione all’inizio di una simulazione è bene che non ci siano celle insature.
Per risolvere questo problema occorre

[attachment=2:17thf235]impostazione-carico-idraulico-iniziale.jpg[/attachment:17thf235]

Aprire la finestra Model -> Modflow -> Package Option (vedi immagine sovrastante)

e selezionare la linguetta Initial Head

[attachment=1:17thf235]Initial-Head.jpg[/attachment:17thf235]

Dovrebbe essere selezionata come opzione di defaul la voce

“Use Default Heads in Spreadsheet Below”

e se guardi nella tabella sottostante dovresti appunto trovare per i tre layer del modello il valore del top del modello che avevi inizialmente scelto in fase di creazione del tuo modello (se non hai modificato questo valore dovrebbe essere 100)

A questo punto è sufficiente scegliere una di queste due strade
1) modificare manualmente i valori della tabella mettendo la quota più elevata presente nella tua superficie topografica (in questo modo tutte le celle sono sicuramente sature)
2) selezionare dal menù a tendina la voce “top of layer 1” (in questo modo ad ogni cella viene associato un carico idraulico uguale alla superrficie topografica in quel punto)

[attachment=0:17thf235]top-of-layer-1.jpg[/attachment:17thf235]

Procedendo in questo modo il problema dovrebbe risolversi

Per quello che riguarda il secondo warning puoi provafe a cliccare sulla linguetta BCF-LPF (la trovi di fianco rispetto alla linguetta “initial head”) e spuntare la casella “Storage coefficient represent specific storage”

Questo warning potrebbe rimanere, ma non influenza i risultati della simulazione. Dovrebbe essere solamente un avvertimento per ricordare al modellatore questa caratteristica del software.

Se il modello continua a non girare magari prova a postare gli altri eventuali warning o messaggi di errore che compariranno.

Saluti!

Sempre in quella directory puoi trovare il documento GV5_commandreference che personalmente ho trovato abbastanza utile. Il manuale che hai riportato ti guida nella creazione del primo modello Modflow e ti illustra le varie funzionalità del programma mentre il commandreference ti indica con dettaglio tutte le caratteristiche dei pacchetti del programma.
Molto utile, sopratutto se hai un problema specifico, è l’help del programma, che con le funzioni di ricerca ti consente spesso di individuare rapidamente il punto che stai cercando.

Questi però sono testi che servono principalmente per usare il programma e non entrano nel merito delle tecniche per trasformare un problema reale in un modello numerico. Non sono cioè manuali di modellazione idrogeologica, ma sono esclusivamente documenti per guidare nell’uso del software.

Saluti!

Buongiorno a tutti,

Una lode particolare al Dott. Bernagozzi, per la puntualità e l’esattezza delle sue risposte.

Per ciò che concerne la tua domanda: hai mai pensato di usare qualche modello stocastico di ricostruzione del sottosuolo? la risposta definitiva (avverto: il prezzo da pagare, in termini di elaborazione e di acquisizione teorica della tecnica, è molto alto, ma le soddisfazioni saranno tante) ai problemi di discretizzazione della griglia potrebbe essere, ad esempio, una elaborazione in T-PROGS su interfaccia Aquaveo GMS. In questo caso servirebbero una buona distribuzione di sondaggi; occorre inoltre avere pazienza per ciò che concerne i problemi di convergenza.