Modeling Variable Viscosity Groundwater Flow with MODFLOW 6 and Model Muse - Tutorial

In order to improve the accuracy of a groundwater flow simulation we need a strong conceptual model, a high quality observation dataset and a numerical code that can handle specific characteristics of the groundwater flow we want to evaluate. MODFLOW 6 can deal with variable density flow with the BUY package and also with variable viscosity with the VSC package that is implemented on the latest versions of Model Muse. This tutorial covers an applied case of variable viscosity modeling due to the injection of salinity and heat to a model with regional groundwater flow. The main features for the implementation of the VSC package are explained and the resulting concentrations on two observation points are evaluated with Flopy codes.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Tutorial in Groundwater Modeling of Faults and Fractures with Modflow6 - Disv on Model Muse

With Model Muse and MODFLOW6 DISV you can model faults and fractures having small cells close to the fracture alignment. In order to achieve small cell sizes close to the fracture aligment the normal MODFLO6 discretization schema (columns and rows) creates a series of unused/unwanted cells, but MODFLOW6 came with two new discretization options that allow us to have local refinements close to areas of interest decreasing the total amount of cells and thus the decreasing computational time.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Coupling local geological data into Model Muse with Python and Aquifer App - Tutorial

Spatial discretization for a mine related groundwater model has to come from geospatial data that has a system of reference (crs). The latest version of Modflow that is Modflow 6 implements the discretized by vertices option (DISV) that allows the creation of triangular, quadtree, voronoi meshes among other options. We have developed some Python scripts to create voronoi meshes from shapefiles where the user has to define a limit polygon, layers (point, line or polygon) and refinement levels. The code generates meshes with adequate performance and gives geospatial output for the final voronoi mesh as well as the intermediate steps.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Geospatial processing for a groundwater model in Model Muse in 14 minutes - Tutorial

Groundwater model creation requires a complete set of spatial data for the different hydraulic parameters, boundary conditions and other model items. Vector and raster data need to be preprocessed, converted, reprojected to fit the requirements of Model Muse. 

This tutorial covers an applied case of raster and vector data processing for a basin.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Basic example of seawater intrusion with variable density MODFLOW6 - BUY and Model Muse - Tutorial

One of the greatest new features of the latest Model Muse version is the implementation of the BUY package for variable density and seawater intrusion modeling. The BUY package was already available in Modflow 6 but the resources to build a model was limited to scripts in Flopy or to build the packages for Modflow by hand, now with Model Muse much stress and pain could be relieved to simulate this type of groundwater model much relevant on the current climate change scenarios.

This tutorial shows the complete procedure to implement variable density / seawater intrusion modeling on a simple coastal aquifer with regional flow. The model runs over two steady / transient stress periods of 50 year each and has a cell with of 20 on the intrusion and 15 layers.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

ModelMuse version 5.2 has been released, explore what is new version.

Even though that the version 5.2 was released on March 22 2024 we just found time to make a video explaning the new features of this version and where they are implemented on the graphical user interfase. This new version of Model Muse allows to simulate variable density and variable viscosity flow among other features of Modflow 6. Have a look on the video and wait for our comming applied tutorials related to these new features.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

A Python class to export Modflow 2005 model 3D geometry and parameters to Vtk - Tutorial

Generating 3D visualizations of groundwater models is essential to analyze the flow regime, perform quality checks and see the interaction of the groundwater body with external factors / boundary conditions. The Flopy library has tools to export the parameters, boundary conditions and results that we have modified and compiled within a Python class. The use of this class allows the generation of Vtk files on a friendly way and in few steps. The tutorial also includes a representation of the parameters generated in ParaView.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Automatic calibration of transient pumping test with MODFLOW 6, Python, FloPy and Scikit Learn - Tutorial

Python has awesome packages for machine learning that can be coupled to groundwater models and perform automatic calibration of hydraulic parameters. This tutorial covers the procedure to implement a neural network based on a set of parameters set with corresponding head values; the study case is on a transient pumping test model with an observation point located around 10 meters away from the pump well. The analysis predicts the calibration parameters from the spatially interpolated observed head values. 

Read More
2 Comments

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

How to export a MODFLOW 6 model grid to shapefile with Python and Flopy - Tutorial

We did this tutorial since the documentation or examples on the topic were not available. Flopy can export the model grid and model attributes to shapefile with a coordinate system of reference for the three types of discretizations of MODFLOW 6. We have done an applied example to export the grid of a discretized by vertices model (DISV) to the ESRI Shapefile model. The tutorial also show options for the final grid representation in geopandas.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Sensibility analysis of transient pumping test with MODFLOW-6, Flopy and SALib - Tutorial

This tutorial covers the whole procedure to perform a sensitivity analysis over a 72 hour pumping test plus recovery on the hydraulic response in an observation piezometer located at 11 meters from the well. the Since the study case it’s a transient model the sensitivities will vary over time and stage of pumping/recovery.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

How to import a Leapfrog geological model to Model Muse with Python - Tutorial

Having a geological model can enhance numerical models since it allows to represent higher accuracy the potential distribution of hydraulic parameters in the horizontal and vertical direction. The process to implement/insert a geological model into a Modflow model is a challenge due to restrictions on proprietary software and spatial tools; we have done the whole procedure to insert a geological model into a MODFLOW-NWT groundwater model with scripts in Python, Pyvista and others.

The code extracts the cell centroids of the modflow model and then compares its position with the different geological units exported from Leapfrog to Vtk. Once the corresponding lithology of the cell is identified a hydraulic parameter is assigned. The tutorial works from Vtks of a geological model done in Leapfrog, but it can work with any Vtk.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Convert geological units from Leapfrog (*.msh) to VTK using Python and GemGIS - Tutorial

This tutorial shows the complete procedure to convert a geological unit as a Leapfrog mesh (*.msh) to the VTK (Visualization Toolkit) format using Python and GemGIS. Follow these steps to seamlessly transfer your geological data for advanced visualization, analysis and comparison with other model results.

Read More
1 Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Geological modeling of sedimentary layers from geospatial rasters with Python and Gempy - Tutorial

Based on a coupled workflow on QGIS and Python it is possible to extract the required information for a Gempy model and run it for defined voxel sizes. This tutorial covers the whole procedure of spatial data preparation, data preprocessing in defined formats and geological modeling with Python and Gempy.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Groundwater modeling of a pumping test over a confined aquifer with MODFLOW-6 and FloPy - Tutorial

We have done an applied groundwater model for a pumping test on a confined aquifer. The wellsite is located on the Clairborne aquifer (Georgia, USA) and has two observations on different layers besides the well itself. The aquifer test has two periods of pumping and recovery that last 3 days each one and the numerical model was constructed with MODFLOW-6 based on a Python script with the FloPy package. Comparison among observed and simulated wells were done with plots in Matplotlib.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Export MODFLOW6 3D model geometry (as Vtk) with Python - Tutorial

Generating 3D visualizations of groundwater models is essential to analyze the flow regime, perform quality checks and see the interaction of the groundwater body with external factors / boundary conditions. The Flopy library has tools to export the parameters, boundary conditions and results that we have modified and compiled within a Python class. The use of this class allows the generation of Vtk files on a friendly way and in few steps. The tutorial also includes a representation of the parameters generated in ParaView.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Docker image containing the parallel version of MODFLOW 6 with application examples

Modflow on parallel processing seemed to lie in the far-off future, but now it’s possible to take advantage of most of your cores to run parallelized versions of MODFLOW 6. Since most MODFLOW modelers work on Windows we wanted to bring something more than a “packed” solution; our goal was to bring an enriched “packed” experience.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Lithium speciation calculation on salt lake brines with Phreeqc and Python - Tutorial

Lithium interaction with the aqueous phase on brines phase is not well understood or at least not well covered by the scientific publications. To evaluate brine composition, brine dynamics and the interaction of the different components is really a challenge with the available tools for geochemical modeling as Phreeqc since it doesn’t appear in all the databases and in most it is only associated with sulfates.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

Online tools for groundwater modeling preprocessing with Model Muse - Tutorial

If you want to speed up some steps in groundwater modeling you might be interested in the tools presented on this tutorial. We have developed some online tools that replace intense and time consuming tasks on desktop software and web servers and provide vector and raster data on the formats required by Model Muse. The tutorial covers the steps of DEM generation as Surfer Grid File (*.grd) with the definition of watersheds and river network based on two online services provided by Hatarilabs.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

How to split Modflow 6 groundwater model with Flopy and MF6Splitter - Tutorial

Modflow 6 has changed in many ways the traditional concepts of groundwater modeling; now a group of models can run under a "simulation" and these models can interchange flow and transport. Under this structure Flopy has implemented a tool to split a groundwater model into several models where their results can be reconstructed and compared to the original entire models. Capabilities and features like this make Flopy and Modflow 6 a great tool for developing advanced groundwater models that simulate complicated tasks or requirements to the groundwater flow regime.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.

 

How to direct download / clip elevations and import them to Model Muse with Python - Tutorial

If you work often with numerical models this process will save you much time. The traditional method of downloading elevation models (Aster DEM) from websites, and reprojecting / clipping them with GIS software can be time consuming and pulls you out of other critical tasks on model conceptualization and calibrations. We have done a complete process in Python under WSL that runs Ubuntu in Windows, this platform was chosen due to the complexities of installing GDAL in Windows.

There are two scripts, one to download the images and the second to translate and clip the image. Finally the elevation are imported as a Surfer Grid File (*.grd) into a model in Model Muse.

Read More
Comment

 

Suscribe to our online newsletter

Subscribe for free newsletter, receive news, interesting facts and dates of our courses in water resources.