November 04, 2024

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

November 04, 2024/ Saul Montoya

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.

October 09, 2024

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

October 09, 2024/ Saul Montoya

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.

September 25, 2024

Geological Modeling of Sedimentary Layers from Geospatial Data with Python and Aquifer App - Tutorial

September 25, 2024/ Saul Montoya

A practical example for constructing a geological model using Python and Aquifer App, based on geospatial data in ESRI Shapefile format and raster data in Tiff format. The tutorial demonstrates the complete procedure of spatial data processing using libraries such as Geopandas, Rasterio, and Pyvista to create surface and orientation files, which are then input into Aquifer App. Finally, the project is run locally in Gempy, generating the geometry of lithology and contacts.

September 04, 2024

Open, modify, run, and read output from HEC-RAS models with Python - Tutorial

September 04, 2024/ Saul Montoya

Python can communicate with other softwares through standards in Windows (COMs) and thus provides a higher level of interoperability for river flow models made in HEC-RAS. The amount of tools provided by HEC-RAS through the COM is huge; in this tutorial we will cover some examples of HEC-RAS model interaction for two applied cases from Python scripts in Jupyter Lab.

August 15, 2024

How to create a Geological Model from a Point Shapefile with Gempy and Aquifer App - Tutorial

August 15, 2024/ Saul Montoya

This is an applied case were we build a geological model only with lithology information stored as a point shapefile. The tutorial covers all steps from raster (array) generation for all surfaces together with the orientation sampling and format of surfaces/orientations as Gempy input files. The generated data was inserted in the Aquifer App that implements an interface to create Gempy models. Finally the lithology and layer surface geometry was exported as Vtk to be represented on Paraview with the initial data to evaluate the accuracy of the simulation.

August 12, 2024

Geological modeling of a petroleum system with Gempy, Gemgis and Aquifer App - Tutorial

August 12, 2024/ Saul Montoya

Aquifer App offers a friendly, clean and powerful way to create Gempy scripts for geological modeling. We have developed an applied case of a 3D petroleum system model for a part of the Williston Basin, USA that contains the Nesson Anticline. Top layer information was provided as raster format where elevation and orientation for random points were extracted with Python codes using the Gemgis package and exported in Gempy input file format. On Aquifer App the processed CSV files were inserted and the geological model was set up with the corresponding geological sequence.

August 01, 2024

NetCDF for water resources with Python for dummies (CHIRPS dataset) - Tutorial

August 01, 2024/ Saul Montoya

NetCDF has become a popular choice for storing and delivering precipitation and water resources related data. Its capacities to store multiple geospatial raster layers over time allow another level of abstraction on data analysis, however the format is of limited use on normal desktop applications and most times we are required to use a programming language such as Python or R.

July 17, 2024

Export 3d Geological Units as Vtk format with Gempy and Aquifer App - Tutorial

July 17, 2024/ Saul Montoya

Aquifer App offers a friendly, clean and powerful way to create Gempy scripts for geological modeling. The platform also allows users to run the scripts and download the whole modeling project that has the 3D geometry of the geological units and faults in Vtk format. This tutorial shows an applied case of geological data visualization, model creation, model run, export and representation of the Vtk files in Paraview.

July 12, 2024

Regional geological modeling with surface topography with Gempy and Aquifer App - Tutorial

July 12, 2024/ Saul Montoya

The coupling of Python scripts with the Aquifer App platform allows us to generate geological models with a versatility not seen before. This is a basic example of regional geological modeling in Aquifer App with surface topography from a digital elevation model in TIF format. The tutorial also makes a 3D representation of the contacts and orientations of the one layer geology in a Jupyter notebook for better analysis of the geological sequence.

July 03, 2024

Geological modeling of a recumbent fold with Gempy and Aquifer App - Tutorial

July 03, 2024/ Saul Montoya

The coupling of Python scripts with the Aquifer App platform allows us to generate geological models with a versatility not seen before. This is a basic example of geological modeling of an recumbent fold in Aquifer App that can create models using Gempy. The tutorial also makes a 3D representation of the contacts and orientations of the layers in a Jupyter notebook for better analysis of the geological sequence.

July 02, 2024

Flood simulation from direct rainfall with Python and Landlab - Tutorial

July 02, 2024/ Saul Montoya

Simulation of direct precipitation events over a surface is important to assess the impacted areas and manage disaster situations. Response to precipitation is evaluated over cells in a 2D model and runoff is routed to the stream networks. Landlab is a Python library for a variety of surface processes and can model shallow water flow over topography with a variety of algorithms. We have developed an applied case of direct rainfall for flood modeling over a geospatial raster where an initial water height is applied and then flows over the surface. Resulting water heights are exported as geospatial rasters.

June 28, 2024

Geological modeling of an anticline structure with Gempy and Aquifer App - Tutorial

June 28, 2024/ Saul Montoya

The coupling of Python scripts with the Aquifer App platform allows us to generate geological models with a versatility not seen before. This is a basic example of geological modeling of an anticline structure in Aquifer App that can create models using Gempy. The tutorial also makes a 3D representation of the contacts and orientations of the layers in a Jupyter notebook for better analysis of the geological sequence.

June 21, 2024

Missing Crop Vegetation Areas Detection with Python and Scikit Learn

June 21, 2024/ Saul Montoya

An applied case for the recognition of missing crop vegetation areas based on a drone orthophoto. Contours have been identified from an enhanced combination of raster bands with a marching squares method to find constant valued contours and then exported as geospatial polygons.

June 18, 2024

Simplify / densify stream network online with Hatari Utils - Tutorial

June 18, 2024/ Saul Montoya

Watershed delineation and stream network determination are common tasks on the hydrological analysis for any area of interest. There is a parameter on the stream network determination that allows us to have simple stream networks that resemble the main water courses or dense stream networks that map all the permanent / temporary streams. This is an applied example of stream network determination on our online platform Hatari Utils where you can specify the number of upstream cells and interactively review the result.

June 11, 2024

Geological modeling of horizontal sedimentary layers with Gempy and Aquifer App - Tutorial

June 11, 2024/ Saul Montoya

The coupling of Python scripts with the Aquifer App platform allows us to generate geological models with a versatility not seen before. This is a basic example of geological modeling of two horizontal layers in Aquifer App that can create models using Gempy. The tutorial also makes a 3D representation of the contacts and orientations of the layers in a Jupyter notebook for better analysis of the geological sequence.

June 07, 2024

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

June 07, 2024/ Saul Montoya

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.

June 04, 2024

Watershed delineation for basins containing internal lakes / reservoirs with Hatari Utils - Tutorial

June 04, 2024/ Saul Montoya

Hydrological basins could have a diversity of shapes due to the morphology and it´s a challenge to have a process that can delineate basins with different characteristics in a short amount of time. Hatari Utils is an online platform for hydrogeological analysis that has a tool for basin delineation, in this tutorial we have tested the capabilities of Hatari Utils to delineate a basin that has an internal lake.

June 03, 2024

Modeling Brine Density vs Concentration Regression Lines with Phreeqc and Aquifer App - Tutorial

June 03, 2024/ Saul Montoya

Complex geochemical simulations are entirely possible to be performed with Phreeqc coupled with Aquifer App and Python. Brines can be simulated at different concentrations to obtain relations that are input of other variable density flow models. In this case we have model one brine, sodium bicarbonate, with the REACTION keyword with moles values that range from 0.5 to 12 moles. Values of mass, volume, concentration, and density were processed in Python from the dataframes generated from Aquifer App.

May 29, 2024

Simulation of water evaporation with Phreeqc and Aquifer App - Tutorial

May 29, 2024/ Saul Montoya

We define evaporation as the process of removing water from the chemical system in Phreeqc and this can be achieved with the REACTION keyword and a negative reaction coefficient. We have developed an applied case of water evaporation from precipitation water and the mix of the resulting water to restore the original volume. The whole Phreeqc setup, simulation and result analysis process has been in Aquifer App.

May 28, 2024

Maps of Days without Rain generation with Python and Rasterio - Tutorial

May 28, 2024/ Saul Montoya

Water resurces management requires not only some direct values as precipitation amounts but also more elaborated data as days without rain. Based on the CHIRPS dataset for year 2022 stored on NetCDF format we have elaborated a map of days without rain for a given country and exported the results as a fully geospatial raster. There is a particular discussion about the data type and the data value to calculate and store value just for the selected location.

Blog