πŸ’§ IGW-NET Β· Quick Tutorial 24 of 31

Tutorial 24: T-PROGS 3D Geologic Model

Build 3D geologic property fields using T-PROGS (transition probability geostatistics) for realistic heterogeneity representation.

IGW-NET Tutorial 24 Prereq: MAGNET4WATER account 6 sections

This tutorial covers

  1. Extract Borehole Data
  2. Define the 3D Domain
  3. Fit the Transition Probability Model
  4. Generate the 3D Geologic Model
  5. Customize the 3D Visualization
  6. What's Next

1Extract Borehole Data

Step 1 β€” Zoom to the Area of Interest

Navigate to the Hastings greater area in southwest Michigan. This region has abundant water well borehole records in the MAGNET Data Center β€” the raw material for T-PROGS modeling.

Step 1 β€” Zoomed to the Hastings greater area in southwest Michigan
Step 1 β€” Zoomed to the Hastings greater area in southwest Michigan

Step 2 β€” Extract Categorized Borehole Data

Click '3D TP' under Analysis Tools. This opens the TP Model Data Options interface. There are multiple ways to add categorized borehole data: upload a shapefile or text file, enter data manually, or extract directly from the MAGNET Data Center.

Click 'Draw TP Area to extract borehole data from data center' under Borehole Well Data. Draw the extraction area on the map, then click 'SaveShape'. After a few moments, the data table populates with borehole records β€” each containing location, well ID, elevation, interval thickness, and material type.

Step 2 β€” TP Model Data Options interface with borehole data extracted from the MAGNET Data Center
Step 2 β€” TP Model Data Options interface with borehole data extracted from the MAGNET Data Center

Step 3 β€” Preview Borehole Data on Map

Click 'ShowOnMap' to display borehole markers. Click any marker to open its Borehole Information Table β€” showing the full lithology log with depth intervals and categorized material types. This is your raw data: the observations that T-PROGS will use to build the 3D geological model.

Step 3 β€” Borehole markers on map with Borehole Information Table showing lithology intervals
Step 3 β€” Borehole markers on map with Borehole Information Table showing lithology intervals

Data from the Data Center β€” Not from Scratch

The MAGNET difference: In traditional T-PROGS workflows, preparing borehole data is the most labor-intensive step β€” finding well logs, digitizing lithology intervals, categorizing materials, formatting files. MAGNET's Data Center has already processed millions of water well records across the USA and Canada. Categorized lithology is ready to extract for any area β€” click, draw, extract. The hours of data preparation become seconds.

Multiple data sources: You're not limited to the Data Center. You can upload your own borehole data (shapefiles or text files) β€” from site investigations, geotechnical reports, or proprietary databases. You can also combine Data Center records with your own field data for maximum coverage. The platform handles all formats through a unified interface.

2Define the 3D Domain

Step 4 β€” Set Top and Bottom Boundaries

Open Domain Attributes ('DomainAttr'). In the Aquifer Attributes tab:

Top Elevation: Select 'DEM Resolution' β€” the land surface from the Digital Elevation Model
Bottom Elevation: Select 'Data Center' β€” the bedrock top surface (interface between unconsolidated and consolidated materials)

This defines the 3D volume that T-PROGS will fill: from ground surface down to bedrock. The unconsolidated materials between these surfaces β€” sand, gravel, clay, till β€” are what the boreholes describe and what T-PROGS will model.

Step 4 β€” Domain Attributes: DEM for top elevation, Data Center bedrock surface for bottom
Step 4 β€” Domain Attributes: DEM for top elevation, Data Center bedrock surface for bottom

3Fit the Transition Probability Model

Step 5 β€” Configure and Run Pre-Simulation

In the TP Model Data Options interface, click 'Create TP Model'. In the TP Model Options interface:

Grid NX: 75 (Y adjusts automatically to maintain roughly square cells)
Zmax: 'MAX DEM' (maximum terrain elevation)
Zmin: 'Min RockTop Elevation' (minimum bedrock surface)

Click 'Pre-Simulation Only'. This calculates the observed transition probabilities from the borehole data β€” counting how often each lithology transition occurs at each lag distance β€” then fits exponential Markov Chain models to the observed data. The processing runs on the cloud; a prompt indicates when it's complete.

Step 5 β€” TP Model Options interface with grid settings and Pre-Simulation parameters
Step 5 β€” TP Model Options interface with grid settings and Pre-Simulation parameters

Step 6 β€” Review the Markov Chain Model

Click 'Proportion/Lens Length Results >>' to view the fitted transition probability models. For each category combination (AQ→AQ, AQ→MAQ, AQ→CM, etc.), you see the observed data (dots) and the fitted exponential model (curve). The fit quality tells you how well the mathematical model captures the geological structure. Key parameters include material proportions (how much of the subsurface is AQ vs. CM) and mean lens lengths (average thickness of each material type).

Step 6 β€” Markov Chain transition probability model: observed data (dots) and fitted exponential models (curves)
Step 6 β€” Markov Chain transition probability model: observed data (dots) and fitted exponential models (curves)

4Generate the 3D Geologic Model

Step 7 β€” Run the Conditional Simulation

Click 'Simulation' in the TP Model Data Options interface. The conditional sequential simulation process begins β€” filling every cell in the 3D grid with a lithology category. Cells containing borehole data receive the observed lithology. Cells without data are assigned categories that honor both the nearby observations and the fitted transition probability structure. Processing runs on the cloud.

Step 7 β€” TP simulation in progress: conditional sequential simulation filling the 3D grid
Step 7 β€” TP simulation in progress: conditional sequential simulation filling the 3D grid

Step 8 β€” View the Four-Panel Output

When simulation completes, four visualization panels appear automatically:

N-S Cross-Section: Centerline slice showing material variability with boreholes
W-E Cross-Section: Perpendicular centerline slice
3D Volume Chart: Interactive full or cropped 3D view with optional borehole display
3D Surface Plot: Interactive 3D surface with fence diagrams, terrain overlay, and layer import

Each material type is color-coded: AQ (bright blue) = aquifer, MAQ = marginal aquifer, PCM = partially confining, CM (red) = confining material. The boreholes are visible in the cross-sections, confirming that the simulation honors the observed data.

Step 8 β€” Four-panel output: N-S and W-E cross-sections, 3D volume chart, and 3D surface plot
Step 8 β€” Four-panel output: N-S and W-E cross-sections, 3D volume chart, and 3D surface plot

5Customize the 3D Visualization

Step 9 β€” Add Terrain Overlay to 3D Surface

Double-click in the 3D Surface Plot to enlarge it. From the first dropdown next to 'TPSurface', select 'Surface with Edges' to show colored cells with black grid lines. Then select 'TPbmp_top.png' from the top option bar to drape a terrain rendering on top of the geologic model β€” the DEM as a translucent overlay. Adjust transparency and vertical offset via the '>>' Other Options menu.

Step 9 β€” Customized 3D surface plot with terrain overlay draped on top of the geologic model
Step 9 β€” Customized 3D surface plot with terrain overlay draped on top of the geologic model

Step 10 β€” Display Wells and Bedrock Surface

Hide the TP model cells (select 'Hidden' from the TPSurface dropdown). Add the bedrock surface by selecting 'TPbmp_bot.png' from the second dropdown β€” the 3D bedrock topography appears. Add color-coded wells by selecting 'TP Well Only' from the last dropdown. Now you see the structural framework: terrain on top, bedrock below, and the borehole evidence in between. The geology that T-PROGS modeled fills this space.

Step 10 β€” 3D display with wells, bedrock surface, and terrain overlay β€” geology hidden to show structure
Step 10 β€” 3D display with wells, bedrock surface, and terrain overlay β€” geology hidden to show structure

Key Concepts

T-PROGS belongs in the conceptual model: T-PROGS generates geology β€” it belongs in Part I of the modeling workflow (Conceptual Model), not in the numerical simulation. The geological model defines the material distribution; the numerical model solves flow through that distribution. Getting the geology right is the foundation β€” everything downstream (flow, transport, calibration, uncertainty) depends on it.

One realization vs. many: This tutorial generates a single T-PROGS realization β€” one possible arrangement of materials that honors the data and statistics. But there are infinitely many plausible arrangements. Tutorials 15–18 showed how Monte Carlo simulation explores this uncertainty for random K fields. The same principle applies to T-PROGS: generate many realizations, simulate flow through each, and accumulate statistics. The uncertainty in geology propagates to uncertainty in predictions.

Data density matters differently: In a data-rich area (many boreholes), T-PROGS is tightly constrained β€” most cells are close to an observation, and the conditional simulation has little freedom. In a data-sparse area, the simulation has more freedom, and different realizations will differ more. The transition probability model provides the statistical "glue" that fills gaps consistently β€” but more data always reduces uncertainty.

From T-PROGS to flow model: The categorical output (AQ, MAQ, PCM, CM) converts directly to a 3D hydraulic conductivity field. Each category maps to a K value (e.g., AQ = 30 m/d, CM = 0.001 m/d). This K field feeds into IGW-NET's flow solver. The geological realism of T-PROGS β€” lenses, pinch-outs, connectivity patterns β€” produces flow behavior that uniform or simple layered models cannot capture.

6What's Next

Continue to 3D visualization and data analysis:

Tutorial 25: 3D Flow Visualization β€” immersive 3D flow fields, water table, and layers
Tutorial 26: 3D Point Data Analysis β€” borehole and well data in 3D
Tutorial 27: DataNET-based Model β€” build a model from federated data services