DOI: 10.5281/zenodo.15674141
title: Numerical investigation of different petrophysical joint inversion schemes

description: Dataset Description

This dataset supports the study on petrophysical joint inversion (PJI) with temperature- and frequency-dependent modeling of subsurface electrical properties. It contains the scripts and resources necessary to reproduce all synthetic model scenarios, inversion experiments, and visualizations presented in the manuscript.

The dataset is organized into three main components:

1. Synthetic Model Generation

Directory: 01_synthetic_modelsThis directory includes the model setup and synthetic data generation scripts for three different synthetic model scenarios:





Model_IAssumes purely electrolytic conduction with no temperature effects (Archie conditions).





01_create_model.py: Defines hydrogeological properties (porosity, water, ice, air content).




02_create_syn_data.py: Generates synthetic seismic and electrical data under temperature-independent conditions.






Model_IIBuilds on Model I but incorporates temperature-dependent electrical conductivity.





01_create_model.py: Identical hydrogeological setup to Model I.




02_create_syn_data.py: Includes temperature effects in electrical forward modeling.






Model_IIIExtends Model II by including cation exchange capacity (CEC) for modeling frequency-dependent conductivity via the dynamic Stern layer model.





01_create_model.py: Defines spatially variable CEC and hydrogeological parameters.




02_create_syn_data.py: Produces synthetic data reflecting both surface conduction and temperature effects.





2. Joint Inversion Execution

Directory: 02_joint_inversionContains inversion scripts for each of the six evaluated PJI schemes (labeled PJI A to F in the manuscript):





PJI_A: Archie implementation applied to synthetic model I




PJI_B: Archie implementation with synthetic model II (temperature-dependent conductivity)




PJI_C: Archie implementation with synthetic model III (surface conduction effects)




PJI_D: Same as PJI C, with porosity constraint




PJI_E: Waxman–Smits implementation with dual-frequency data on model III




PJI_F: Same as PJI E, with porosity constraint



Each subdirectory contains:





03_joint_inversion.py: Performs the inversion for the corresponding PJI configuration using the appropriate synthetic data.



3. Visualization and Post-Processing





04_plot_truemodel.py: Visualizes the synthetic hydrogeological and geophysical parameter distributions across models I–III.




05_plot_results.py: Plots the inversion results for all PJI schemes, including hydrogeological and geophysical properties.



Code Availability

All scripts for running the petrophysical joint inversion, including the core inversion engine and petrophysical models, are available on GitHub:





https://github.com/geophilik/petrojoint




https://github.com/florian-wagner/four-phase-inversion



These repositories provide further documentation and source code for extending or adapting the inversion schemes.

Citation

If you use this dataset or scripts, please cite the associated manuscript and include appropriate attribution to the original software repositories.