Hydrology to ERT#

This workflow converts hydrologic model states into resistivity and forward ERT responses.

Steps#

  1. Load MODFLOW/ParFlow outputs.

  2. Convert water content to resistivity.

  3. Run forward modeling or inversion.

import numpy as np

from PyHydroGeophysX.model_output.water_content import MODFLOWWaterContent
from PyHydroGeophysX.petrophysics import water_content_to_resistivity

# idomain marks the active cells of the MODFLOW grid. The loader needs it to
# unpack the WaterContent records into (nlay, nrow, ncol), so it cannot be None.
idomain = np.loadtxt("examples/data/modflow/id.txt")

wc = MODFLOWWaterContent("examples/data/modflow", idomain=idomain).load_timestep(0, nlay=3)
rho = water_content_to_resistivity(wc, rhos=100.0, n=2.0, porosity=0.3)