Inverse rock physics modeling code for "Revealing the diversity of pore-filling fluids and cements in Martian crust"

This repository contains the Python codes for foward and INVERSE rock physics modeling.

The codes have been tested for Python version 3.8.10 (tags/v3.8.10:3d8993a) with WINDOWS 11 system (22631.3880).

## Potentail prerequisties for external libraries

This repository is designed to run with the standard libraries of Python version 3.8.10.

Upon that, some example programs with result plotting require the external library 'matplotlib'.

Dependancies on external libraries may arise only when comparing the results with other methods:

1. /libs/rpm_tools.py: To compare results from this repository with the external Python library 'rockphypy'.

2. /libs/Func_differential_solutions.py: To compare solutions from methods of this repository with the external Python library 'scipy'.

The 'import' statements for these external libraries are commented by default and should be uncommented when needed.

## Repository description

The repository is organized in the following:

libs/

    rpm_tools.py: Subroutines for essential tool functions in rock physics modeling.

    rpm_uncemented_sediment.py: Subroutines for rock physics modeling of uncemented seidments.

    rpm_cemented_sediment.py: Subroutines for rock physics modeling of cemented seidments.

    rpm_fractured_rock.py: Subroutines for rock physics modeling of fractured rocks.

    rpm_inv.py: Subroutines for INVERSE rock physics modeling.

    Func_differential_solutions.py:  Subroutines to solve differential equations.

    Func_property_PT.py: Subroutines to calcuate fliud physical properties across various pressure-temperature conditions.

    Func_binary_search.py: Subroutines to conduct the binary search for INVERSE rock physics modeling.

   

Settings.py: Settings essentail for rock physcis modeling. [Need to Be Set By User]

ex_rpm_uncem_sed.py: Example program and result plots for rock physics modeling of uncemented seidments.

ex_rpm_cem_sed.py: Example program and result plots for rock physics modeling of cemented seidments.

ex_rpm_rock_one_case.py: Example program for rock physics modeling of one rock case.

ex_rpm_rock.py: Example program and result plots for rock physics modeling of various rock models.

ex_inv_AR.py: Example program and result plots for INVERSE rock physics modeling of porosity and pore aspect ratio.

ex_inv_poreV.py: Example program for INVERSE rock physics modeling of cement/fluid volume in pore, with a given range of porosity and pore aspect ratio.

ex_inv_poreK.py: Example program and result plots for INVERSE rock physics modeling of pore-filling media, with a fixed porosity and pore aspect ratio.

ex_inv_pore_mat.py: Example program and result plots for INVERSE rock physics modeling of pore-filling media, for a given range of porosity and pore aspect ratio.

ex_brine.py: Example program and result plots for brine properties across various pressure-temperature conditions.

ex_gas.py: Example program and result plots for gas properties across various pressure-temperature conditions.