rock-physics-open 0.3.0__py3-none-any.whl → 0.3.1__py3-none-any.whl

rock_physics_open/equinor_utilities/various_utilities/gassmann_mod.py

@@ -1,7 +1,26 @@
+import numpy as np
+import numpy.typing as npt
+
 from rock_physics_open.equinor_utilities import std_functions
 
 
-def gassmann_model(k_min, k_fl, rho_fl, k_dry, mu, rho_dry, por):
+def gassmann_model(
+    k_min: npt.NDArray[np.float64],
+    k_fl: npt.NDArray[np.float64],
+    rho_fl: npt.NDArray[np.float64],
+    k_dry: npt.NDArray[np.float64],
+    mu: npt.NDArray[np.float64],
+    rho_dry: npt.NDArray[np.float64],
+    por: npt.NDArray[np.float64],
+) -> tuple[
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+]:
     """
     Gassmann model to go from dry rock to saturated state.
 
@@ -25,7 +44,7 @@ def gassmann_model(k_min, k_fl, rho_fl, k_dry, mu, rho_dry, por):
     Returns
     -------
     tuple
-        vp_sat, vs_sat, rho_sat, ai_sat, vpvs_sat, k_sat, mu : (np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray).
+        vp_sat, vs_sat, rho_sat, ai_sat, vpvs_sat, k_sat, mu : np.ndarray
         vp_sat, vs_sat:  saturated velocities [m/s], rho_sat: saturated density [kg/m^3], ai_sat: saturated acoustic
         impedance [kg/m^3 x m/s], vpvs_sat: saturated velocity ratio [unitless], k_sat, mu: saturated bulk modulus and
         shear modulus (the latter unchanged from dry state) [Pa].

rock_physics_open/equinor_utilities/various_utilities/gassmann_sub_mod.py

@@ -1,17 +1,28 @@
+import numpy as np
+import numpy.typing as npt
+
 from rock_physics_open.equinor_utilities import std_functions
 
 
 def gassmann_sub_model(
-    k_min,
-    k_fl_orig,
-    rho_fl_orig,
-    k_fl_sub,
-    rho_fl_sub,
-    k_sat_orig,
-    mu,
-    rho_sat_orig,
-    por,
-):
+    k_min: npt.NDArray[np.float64],
+    k_fl_orig: npt.NDArray[np.float64],
+    rho_fl_orig: npt.NDArray[np.float64],
+    k_fl_sub: npt.NDArray[np.float64],
+    rho_fl_sub: npt.NDArray[np.float64],
+    k_sat_orig: npt.NDArray[np.float64],
+    mu: npt.NDArray[np.float64],
+    rho_sat_orig: npt.NDArray[np.float64],
+    por: npt.NDArray[np.float64],
+) -> tuple[
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+]:
     """
     Gassmann model to go from one saturated state to another.
 
@@ -24,7 +35,7 @@ def gassmann_sub_model(
     rho_fl_orig : np.ndarray
         Original fluid density [lg/m^3].
     k_fl_sub : np.ndarray
-        Sunstituted fluid  bulk modulus [Pa].
+        Substituted fluid  bulk modulus [Pa].
     rho_fl_sub : np.ndarray
         Substituted fluid density [lg/m^3].
     k_sat_orig : np.ndarray
@@ -39,7 +50,7 @@ def gassmann_sub_model(
     Returns
     -------
     tuple
-        vp_sat, vs_sat, rho_sat, ai_sat, vpvs_sat, k_sat, mu : (np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray).
+        vp_sat, vs_sat, rho_sat, ai_sat, vpvs_sat, k_sat, mu : np.ndarray
         vp_sat, vs_sat:  saturated velocities [m/s], rho_sat: saturated density [kg/m^3], ai_sat: saturated acoustic
         impedance [kg/m^3 x m/s], vpvs_sat: saturated velocity ratio [unitless], k_sat, mu: saturated bulk modulus and
         shear modulus (the latter unchanged from dry state) [Pa].

rock_physics_open/equinor_utilities/various_utilities/hs_average.py

@@ -1,7 +1,26 @@
+import numpy as np
+import numpy.typing as npt
+
 from rock_physics_open.equinor_utilities import std_functions
 
 
-def hs_average(k1, mu1, rhob1, k2, mu2, rhob2, f):
+def hs_average(
+    k1: npt.NDArray[np.float64],
+    mu1: npt.NDArray[np.float64],
+    rhob1: npt.NDArray[np.float64],
+    k2: npt.NDArray[np.float64],
+    mu2: npt.NDArray[np.float64],
+    rhob2: npt.NDArray[np.float64],
+    f: npt.NDArray[np.float64],
+) -> tuple[
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+]:
     """
     BMix of two phases by Hashin-Shtrikman model. Derived properties are also returned.
 

rock_physics_open/equinor_utilities/various_utilities/pressure.py

@@ -1,7 +1,15 @@
 import numpy as np
+import numpy.typing as npt
 
 
-def pressure(rho, tvd_msl, water_depth, p_form, tvd_p_form, n):
+def pressure(
+    rho: npt.NDArray[np.float64],
+    tvd_msl: npt.NDArray[np.float64],
+    water_depth: float,
+    p_form: float,
+    tvd_p_form: float,
+    n: float,
+) -> tuple[npt.NDArray[np.float64], npt.NDArray[np.float64]]:
     """
     Function to estimate overburden pressure and vertical effective stress (lithostatic pressure)
     based on density.

rock_physics_open/equinor_utilities/various_utilities/reflectivity.py

@@ -1,9 +1,19 @@
+from typing import Literal, cast
+
 import numpy as np
+import numpy.typing as npt
 
 from rock_physics_open.equinor_utilities import gen_utilities, std_functions
 
 
-def reflectivity(vp_inp, vs_inp, rho_inp, theta=0.0, k=2.0, model="AkiRichards"):
+def reflectivity(
+    vp_inp: npt.NDArray[np.float64],
+    vs_inp: npt.NDArray[np.float64],
+    rho_inp: npt.NDArray[np.float64],
+    theta: float = 0.0,
+    k: float = 2.0,
+    model: Literal["AkiRichards", "SmithGidlow"] = "AkiRichards",
+) -> tuple[npt.NDArray[np.float64], npt.NDArray[np.bool_]]:
     """
     Reflectivity model according to Aki and Richards or Smith and Gidlow for weak contrasts
     and angles less than critical angle.
@@ -35,12 +45,17 @@ def reflectivity(vp_inp, vs_inp, rho_inp, theta=0.0, k=2.0, model="AkiRichards")
         idx_inp: index to accepted part of the input arrays [bool].
     """
 
-    vp, vs, rho, theta, k = gen_utilities.dim_check_vector(
-        (vp_inp, vs_inp, rho_inp, theta, k)
+    vp, vs, rho, theta_, k_ = cast(
+        list[npt.NDArray[np.float64]],
+        gen_utilities.dim_check_vector((vp_inp, vs_inp, rho_inp, theta, k)),
     )
 
-    idx_inp, (vp, vs, rho, theta, k) = gen_utilities.filter_input_log(
-        [vp, vs, rho, theta, k], positive=True
+    idx_inp, (vp, vs, rho, theta_, k_) = cast(
+        tuple[npt.NDArray[np.bool_], list[npt.NDArray[np.float64]]],
+        gen_utilities.filter_input_log(
+            [vp, vs, rho, theta_, k_],
+            positive=True,
+        ),
     )
 
     if np.any(~idx_inp):
@@ -69,16 +84,17 @@ def reflectivity(vp_inp, vs_inp, rho_inp, theta=0.0, k=2.0, model="AkiRichards")
                 > 1
             ) * "s"
             raise ValueError(
-                "{0:} reflectivity: Missing or illegal values in input log{1:}: {2:}interpolation of input log{1:} "
-                "is needed\n".format(model, log_str, pl_str)
+                "{0:} reflectivity: Missing or illegal values in input log{1:}: {2:}interpolation of input log{1:} is needed\n".format(
+                    model, log_str, pl_str
+                )
             )
 
     if model == "AkiRichards":
-        refl_coef = std_functions.aki_richards(vp, vs, rho, theta, k)
+        refl_coef = std_functions.aki_richards(vp, vs, rho, theta_, k_)
     elif model == "SmithGidlow":
-        refl_coef = std_functions.smith_gidlow(vp, vs, rho, theta, k)
+        refl_coef = std_functions.smith_gidlow(vp, vs, rho, theta_, k_)
     else:
-        raise ValueError(
+        raise ValueError(  # pyright: ignore[reportUnreachable] | Kept for backward compatibility
             f'{__file__}: unknown model: {model}, should be one of "AkiRichards", "SmithGidlow"'
         )
 

rock_physics_open/equinor_utilities/various_utilities/timeshift.py

@@ -1,7 +1,13 @@
 import numpy as np
+import numpy.typing as npt
 
 
-def time_shift_pp(tvd, vp_base, vp_mon, multiplier):
+def time_shift_pp(
+    tvd: npt.NDArray[np.float64],
+    vp_base: npt.NDArray[np.float64],
+    vp_mon: npt.NDArray[np.float64],
+    multiplier: int,
+) -> tuple[npt.NDArray[np.float64], npt.NDArray[np.float64]]:
     """
     Cumulative time shift calculation for 4D case. According to Equinor standard
     the time shift is negative for an increase in velocity from base to monitor
@@ -43,7 +49,14 @@ def time_shift_pp(tvd, vp_base, vp_mon, multiplier):
     return owt_pp_shift, twt_pp_shift
 
 
-def time_shift_ps(tvd, vp_base, vp_mon, vs_base, vs_mon, multiplier):
+def time_shift_ps(
+    tvd: npt.NDArray[np.float64],
+    vp_base: npt.NDArray[np.float64],
+    vp_mon: npt.NDArray[np.float64],
+    vs_base: npt.NDArray[np.float64],
+    vs_mon: npt.NDArray[np.float64],
+    multiplier: int,
+) -> npt.NDArray[np.float64]:
     """
     Cumulative time shift calculation for 4D case. According to Equinor standard
     the time shift is negative for an increase in velocity from base to monitor

rock_physics_open/equinor_utilities/various_utilities/vp_vs_rho_set_statistics.py

@@ -1,6 +1,8 @@
 import os
+from typing import Literal, cast
 
 import numpy as np
+import numpy.typing as npt
 import pandas as pd
 from sklearn.metrics import r2_score
 
@@ -8,18 +10,18 @@ from .display_result_statistics import disp_result_stats
 
 
 def vp_vs_rho_stats(
-    vp_observed,
-    vs_observed,
-    rho_observed,
-    vp_estimated,
-    vs_estimated,
-    rho_estimated,
-    fname,
-    estimated_set_names,
-    well_names,
-    file_mode="a",
-    disp_results=True,
-):
+    vp_observed: npt.NDArray[np.float64] | list[npt.NDArray[np.float64]],
+    vs_observed: npt.NDArray[np.float64] | list[npt.NDArray[np.float64]],
+    rho_observed: npt.NDArray[np.float64] | list[npt.NDArray[np.float64]],
+    vp_estimated: npt.NDArray[np.float64] | list[npt.NDArray[np.float64]],
+    vs_estimated: npt.NDArray[np.float64] | list[npt.NDArray[np.float64]],
+    rho_estimated: npt.NDArray[np.float64] | list[npt.NDArray[np.float64]],
+    fname: str,
+    estimated_set_names: str | list[str],
+    well_names: str | list[str],
+    file_mode: Literal["a", "w"] = "a",
+    disp_results: bool = True,
+) -> None:
     """
     Utility to estimate statistics between vp-vs-rho sets - observed and estimated values. The results are displayed
     on screen (optional) and saved to a .csv file. If the file exists, the results will be appended.
@@ -83,33 +85,42 @@ def vp_vs_rho_stats(
     ]
     est_frame = pd.DataFrame(columns=est_frame_columns, index=estimated_set_names)
     est_frame.index.name = "Estimated set name"
-    est_frame.iloc[:, 0] = well_names
+    est_frame.iloc[:, 0] = well_names  # pyright: ignore[reportArgumentType]
 
     # If inputs are found to satisfy expectations in _verify, and they are numpy arrays, cast to lists, and run through
     if isinstance(vp_observed, np.ndarray):
         vp_observed = [vp_observed]
+    if isinstance(vs_observed, np.ndarray):
         vs_observed = [vs_observed]
+    if isinstance(rho_observed, np.ndarray):
         rho_observed = [rho_observed]
+    if isinstance(vp_estimated, np.ndarray):
         vp_estimated = [vp_estimated]
+    if isinstance(vs_estimated, np.ndarray):
         vs_estimated = [vs_estimated]
+    if isinstance(rho_estimated, np.ndarray):
         rho_estimated = [rho_estimated]
 
     for i in range(len(vp_observed)):
-        res = []
+        res: list[float] = []
         for obs, est in zip(
             [vp_observed[i], vs_observed[i], rho_observed[i]],
             [vp_estimated[i], vs_estimated[i], rho_estimated[i]],
         ):
-            res.append(np.mean(np.abs((est.flatten() - obs.flatten()) / obs.flatten())))
-            res.append(
-                np.sqrt(
-                    np.mean(np.square((est.flatten() - obs.flatten()) / obs.flatten()))
-                )
+            rmae = float(
+                np.mean(np.abs((est.flatten() - obs.flatten()) / obs.flatten()))
+            )
+            rrmse: float = np.sqrt(
+                np.mean(np.square((est.flatten() - obs.flatten()) / obs.flatten()))
             )
-            res.append(r2_score(obs.flatten(), est.flatten()))
+            r2 = cast(float, r2_score(obs.flatten(), est.flatten()))
+
+            res.append(rmae)
+            res.append(rrmse)
+            res.append(r2)
 
         res_dict = dict(zip(est_frame_columns[1:], res))
-        est_frame.iloc[i, 1:] = res_dict
+        est_frame.iloc[i, 1:] = res_dict  # pyright: ignore[reportArgumentType]
         if disp_results:
             disp_result_stats(
                 estimated_set_names[i], res, est_frame_columns[1:], values_only=True
@@ -122,7 +133,12 @@ def vp_vs_rho_stats(
         est_frame.to_csv(fname, mode=file_mode)
 
 
-def _verify(*args, set_names=None, well_names=None, file_mode=None):
+def _verify(
+    *args: npt.NDArray[np.float64] | list[npt.NDArray[np.float64]],
+    set_names: list[str],
+    well_names: list[str],
+    file_mode: Literal["a", "w"],
+):
     """Verify that arguments are either numpy arrays or lists of numpy arrays.
     Raises
     ------
@@ -142,7 +158,7 @@ def _verify(*args, set_names=None, well_names=None, file_mode=None):
         if isinstance(arg, np.ndarray):
             arg = [arg]
         for this_arg in arg:
-            if not isinstance(this_arg, np.ndarray):
+            if not isinstance(this_arg, np.ndarray):  # pyright: ignore[reportUnnecessaryIsInstance] | For backward compatibility
                 raise ValueError(f"{__file__}: input not numpy array: {type(arg)}")
             if np.any(np.isnan(this_arg)):
                 raise ValueError(f"{__file__}: input contains NaNs")
@@ -151,4 +167,4 @@ def _verify(*args, set_names=None, well_names=None, file_mode=None):
         if not len(arg) == len(set_names) == len(well_names):
             raise ValueError(f"{__file__}: mismatch in argument lengths")
     if not (file_mode == "a" or file_mode == "w"):
-        raise ValueError(f'{__file__}: file_mode must be one of ["a", "w"]')
+        raise ValueError(f'{__file__}: file_mode must be one of ["a", "w"]')  # pyright: ignore[reportUnreachable] | For backward compatibility

rock_physics_open/equinor_utilities/various_utilities/vrh_3_min.py

@@ -1,9 +1,31 @@
+import numpy as np
+import numpy.typing as npt
+
 from rock_physics_open.equinor_utilities import std_functions
 
 
 def min_3_voigt_reuss_hill(
-    vp1, vs1, rhob1, f1, vp2, vs2, rhob2, f2, vp3, vs3, rhob3, f3
-):
+    vp1: npt.NDArray[np.float64],
+    vs1: npt.NDArray[np.float64],
+    rhob1: npt.NDArray[np.float64],
+    f1: npt.NDArray[np.float64],
+    vp2: npt.NDArray[np.float64],
+    vs2: npt.NDArray[np.float64],
+    rhob2: npt.NDArray[np.float64],
+    f2: npt.NDArray[np.float64],
+    vp3: npt.NDArray[np.float64],
+    vs3: npt.NDArray[np.float64],
+    rhob3: npt.NDArray[np.float64],
+    f3: npt.NDArray[np.float64],
+) -> tuple[
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+    npt.NDArray[np.float64],
+]:
     """
     Mix of three phases by Voigt-Reuss-Hill model. The fractions should add up to 1 with input of vp, vs and rho.
 

rock_physics_open/sandstone_models/friable_models.py

@@ -124,7 +124,7 @@ def friable_model_dry(k_min, mu_min, phi, p_eff, phi_c, coord_num_func, n, shear
         Critical porosity [fraction].
     coord_num_func : str
         Indication if coordination number should be calculated from porosity or kept constant.
-    n : float
+    n : float | None
         Coordination number [unitless].
     shear_red : float
         Shear reduction factor [fraction].
@@ -135,18 +135,17 @@ def friable_model_dry(k_min, mu_min, phi, p_eff, phi_c, coord_num_func, n, shear
         k, mu : (np.ndarray, np.ndarray).
         Bulk modulus k [Pa], shear modulus mu [Pa] of dry rock.
     """
-    # Expand floats to arrays
-    phi, phi_c, n, shear_red = gen_utilities.dim_check_vector(
-        (phi, phi_c, n, shear_red)
+    # Expand floats to arrays, check for equal length
+    phi, phi_c, shear_red, k_min, mu_min, p_eff = gen_utilities.dim_check_vector(
+        (phi, phi_c, shear_red, k_min, mu_min, p_eff)
     )
-
     # Valid porosity values are less or equal to the critical porosity
     # Use filter_input_log to remove values that do not comply with this
     (
         idx_phi,
-        (k_min, mu_min, phi, p_eff, phi_c, n, shear_red, _),
+        (phi, phi_c, shear_red, k_min, mu_min, p_eff, _),
     ) = gen_utilities.filter_input_log(
-        (k_min, mu_min, phi, p_eff, phi_c, n, shear_red, phi_c - phi)
+        (phi, phi_c, shear_red, k_min, mu_min, p_eff, phi_c - phi)
     )
 
     # Dry rock properties of high-porosity end member calculated with

rock_physics_open/sandstone_models/patchy_cement_model.py

@@ -1,6 +1,11 @@
 import numpy as np
 
 from rock_physics_open.equinor_utilities import std_functions
+from rock_physics_open.equinor_utilities.gen_utilities import (
+    dim_check_vector,
+    filter_input_log,
+    filter_output,
+)
 
 from .constant_cement_models import constant_cement_model_dry
 from .friable_models import friable_model_dry
@@ -299,6 +304,14 @@ def patchy_cement_model_dry(
         mu : dry rock shear modulus [Pa],
         rho_dry : dry rock density [kg/m3].,
     """
+    # There are cases which suffer from a lack of consistency check at this stage,
+    # add dim_check_vector and filter input/output
+    phi, k_min, mu_min, rho_min, k_cem, mu_cem, rho_cem, p_eff = dim_check_vector(
+        (phi, k_min, mu_min, rho_min, k_cem, mu_cem, rho_cem, p_eff)
+    )
+    (idx, (phi, k_min, mu_min, rho_min, k_cem, mu_cem, rho_cem, p_eff)) = (
+        filter_input_log((phi, k_min, mu_min, rho_min, k_cem, mu_cem, rho_cem, p_eff))
+    )
 
     k_zero, mu_zero = std_functions.hashin_shtrikman_walpole(
         k_cem, mu_cem, k_min, mu_min, FRAC_CEM_UP, bound="lower"
@@ -366,4 +379,6 @@ def patchy_cement_model_dry(
 
     rho_dry = (1 - phi - frac_cem) * rho_min + frac_cem * rho_cem
 
+    k_dry, mu, rho_dry = filter_output(idx, (k_dry, mu, rho_dry))
+
     return k_dry, mu, rho_dry

rock_physics_open/ternary_plots/ternary_plot_utilities.py

@@ -183,15 +183,15 @@ def _make_mesh(ax):
     # Tick mark text
     xy1 = _ternary_coord_trans(1 - i, i, np.zeros(len(i)))
     for j in range(len(i)):
-        text_handles.append(ax.text(xy1[j, 0], xy1[j, 1] - 0.025, "%.1f" % (i[j])))
+        text_handles.append(ax.text(xy1[j, 0], xy1[j, 1] - 0.025, "%.1f" % (i[j, 0])))
     xy1 = _ternary_coord_trans(i, np.zeros(len(i)), 1 - i)
     for j in range(len(i)):
         text_handles.append(
-            ax.text(xy1[j, 0] - 0.05, xy1[j, 1] + 0.025, "%.1f" % (i[j]))
+            ax.text(xy1[j, 0] - 0.05, xy1[j, 1] + 0.025, "%.1f" % (i[j, 0]))
         )
     xy1 = _ternary_coord_trans(np.zeros(len(i)), 1 - i, i)
     for j in range(len(i)):
-        text_handles.append(ax.text(xy1[j, 0], xy1[j, 1] + 0.025, "%.1f" % (i[j])))
+        text_handles.append(ax.text(xy1[j, 0], xy1[j, 1] + 0.025, "%.1f" % (i[j, 0])))
     plt.setp(text_handles, fontsize=10)
 
     return text_handles, line_handles

rock_physics_open/version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '0.3.0'
-__version_tuple__ = version_tuple = (0, 3, 0)
+__version__ = version = '0.3.1'
+__version_tuple__ = version_tuple = (0, 3, 1)
 
 __commit_id__ = commit_id = None

{rock_physics_open-0.3.0.dist-info → rock_physics_open-0.3.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: rock_physics_open
-Version: 0.3.0
+Version: 0.3.1
 Summary: Equinor Rock Physics Module
 Author-email: Harald Flesche <hfle@equinor.com>, Eivind Jahren <ejah@equinor.com>, Jimmy Zurcher <jiz@equinor.com>
 Maintainer-email: Harald Flesche <hfle@equinor.com>, Eirik Ola Aksnes <eoaksnes@equinor.com>, Sivert Utne <sutn@equinor.com>, Einar Wigum Arbo <earb@equinor.com>
@@ -76,7 +76,7 @@ Alternatively, you can update the dependencies in your `pyproject.toml` file:
 <!-- x-release-please-start-version -->
 ```toml
 dependencies = [
-    "rock-physics-open == 0.3.0",
+    "rock-physics-open == 0.3.1",
 ]
 ```
 <!-- x-release-please-end-version -->

{rock_physics_open-0.3.0.dist-info → rock_physics_open-0.3.1.dist-info}/RECORD

@@ -1,5 +1,5 @@
 rock_physics_open/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-rock_physics_open/version.py,sha256=5zTqm8rgXsWYBpB2M3Zw_K1D-aV8wP7NsBLrmMKkrAQ,704
+rock_physics_open/version.py,sha256=gGLpQUQx-ty9SEy9PYw9OgJWWzJLBnCpfJOfzL7SjlI,704
 rock_physics_open/equinor_utilities/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 rock_physics_open/equinor_utilities/anisotropy.py,sha256=toEsuIW_gHr34gZQpam3sQTl0ap6nAq1pW5hcX9GOHI,4469
 rock_physics_open/equinor_utilities/conversions.py,sha256=d6SGPbzYsmYBvhDqnohg1lPjAeYGOc0z_vdsuFeNlb4,194
@@ -16,15 +16,13 @@ rock_physics_open/equinor_utilities/gen_utilities/dict_to_float.py,sha256=-fiOsN
 rock_physics_open/equinor_utilities/gen_utilities/dim_check_vector.py,sha256=2ebRxWL4XFpnLsbi0j43KQff9ydRyR8hOB1KSqQY8JI,3964
 rock_physics_open/equinor_utilities/gen_utilities/filter_input.py,sha256=K-DUsS8nVOWick6CsnsWrUWnOtzPq5ofc_Z8kPBU_hs,5419
 rock_physics_open/equinor_utilities/gen_utilities/filter_output.py,sha256=iwoHuyzJtksEJKjuLwScd96Bov6SyifU-AqvosZ9VD0,3528
-rock_physics_open/equinor_utilities/machine_learning_utilities/__init__.py,sha256=5QwVh4dbH9sTaj33ok_7ZjGCxOjEKnU6IFMBKl4VZJQ,909
+rock_physics_open/equinor_utilities/machine_learning_utilities/__init__.py,sha256=tXpPTpkviwEbok3Ff8R7m4Og8d97fjdsZaaM599jjEA,469
 rock_physics_open/equinor_utilities/machine_learning_utilities/base_pressure_model.py,sha256=UHBdYODHbXey4IwAhCbedcb5SwGCLkVJK98zKIOaIfY,4722
 rock_physics_open/equinor_utilities/machine_learning_utilities/dummy_vars.py,sha256=y2hWACUX9sk-kRH7FNTaO2-9ZNDHR-HvZzSpIz9ILxM,1479
 rock_physics_open/equinor_utilities/machine_learning_utilities/exponential_model.py,sha256=ZWtA8aNKVDkf8vAyH8IzJx1l8feipxQOGhITo5Ytv0g,3738
-rock_physics_open/equinor_utilities/machine_learning_utilities/friable_pressure_models.py,sha256=cetCzp8QzJZmfp0Mxga0-VNid9y2jObO5q9HxgK7hiE,7017
-rock_physics_open/equinor_utilities/machine_learning_utilities/import_ml_models.py,sha256=BFGnvnOfLHlBEnEy2BxbLteb4kUnURCNm2OS2m7ejUY,3195
-rock_physics_open/equinor_utilities/machine_learning_utilities/patchy_cement_pressure_models.py,sha256=7F60cS_TrffiMozrVdNkhl2lKcyeSQgQoFsGCDXokEc,8361
+rock_physics_open/equinor_utilities/machine_learning_utilities/import_ml_models.py,sha256=a8KJTfeUCNz0zuDKwS4M5ZFIiWdIA8oHvEPxG4FzYvM,2391
 rock_physics_open/equinor_utilities/machine_learning_utilities/polynomial_model.py,sha256=WO6yOoygTD8B73SVWMoS5W-oMTAYm3yo8Ztj-LPR3s4,3673
-rock_physics_open/equinor_utilities/machine_learning_utilities/run_regression.py,sha256=L0bCEsMCvcOn6Ji9oDLztY-TaEEYvp1b7bYZQSJ_-II,4799
+rock_physics_open/equinor_utilities/machine_learning_utilities/run_regression.py,sha256=7Ocw-tJtedJeuwcmCBQM_pHDmQilvxc2tO_oa-8WIdE,4826
 rock_physics_open/equinor_utilities/machine_learning_utilities/sigmoidal_model.py,sha256=dIPtjPcjCOBJtQwuUfRRci8rgBfFwwcPcbuKPqqL7VM,7315
 rock_physics_open/equinor_utilities/optimisation_utilities/__init__.py,sha256=8IZ_5Yq83YUB7_cxTkC8WeMkQGUHt83IHzA6GJ4rxDc,364
 rock_physics_open/equinor_utilities/optimisation_utilities/opt_subst_utilities.py,sha256=X77V2hi6PKpKwrf9ygPXeIjJxr-LkztucVXlyAjMBHU,14164
@@ -32,30 +30,30 @@ rock_physics_open/equinor_utilities/snapshot_test_utilities/__init__.py,sha256=l
 rock_physics_open/equinor_utilities/snapshot_test_utilities/compare_snapshots.py,sha256=PcbCM29iKtVUyRDB9kEFIbboX4fLT_Co2aaztCyqE20,5729
 rock_physics_open/equinor_utilities/snapshot_test_utilities/snapshots.py,sha256=_UvX78IoPDPTCAbblHD9fH_Yn-LwkDOA9MuU69bPdbM,2691
 rock_physics_open/equinor_utilities/std_functions/__init__.py,sha256=MLOfMxRazS9tGOVvuLkUp-E4b6jRWOP4YCLf8si624c,1121
-rock_physics_open/equinor_utilities/std_functions/backus_ave.py,sha256=sYaXSfqjR69U_UFaMmrD5aEbR_0FrL7ZIWGg-Yi8QNo,1607
-rock_physics_open/equinor_utilities/std_functions/dvorkin_nur.py,sha256=TerlLcuP92ihf5vHaJES5KLyleG3CR7NWSExDI1jLIY,2138
-rock_physics_open/equinor_utilities/std_functions/gassmann.py,sha256=sxZ_d-mJXAxC_qfYuYwmpOLVwh3jQvQLCH7ezQx02nM,3641
-rock_physics_open/equinor_utilities/std_functions/hashin_shtrikman.py,sha256=dz_-n-nObXQDStveJYWofSV1fMJQoDDPLg6sUNXOVFA,5902
-rock_physics_open/equinor_utilities/std_functions/hertz_mindlin.py,sha256=KcDx9JRn45aaQjcx6n1GCMqiNGf_PoqGiM4QkFC68-E,1403
-rock_physics_open/equinor_utilities/std_functions/moduli_velocity.py,sha256=31M1qHDuuMx_iGi8POCnPmmTM07MlR7Hj3eI2Ejm2jk,1251
-rock_physics_open/equinor_utilities/std_functions/reflection_eq.py,sha256=DA-Av9fudwomZb9cFsA8vfhjKIu_QwK50tRqgId1WM8,2544
-rock_physics_open/equinor_utilities/std_functions/rho.py,sha256=3yzfJ4felR2DQHSGu2sLDP2-OvplEtogHbJEXb8U-lk,1406
-rock_physics_open/equinor_utilities/std_functions/voigt_reuss_hill.py,sha256=M5qRiDxWndeWdXqBuV1HxwdyTifUGu_OvwfVcoNRvEI,3064
-rock_physics_open/equinor_utilities/std_functions/walton.py,sha256=I2mQuSTu4OOO9KHQvmgRHr1AuzTYIWSby9NofYMmq9s,1250
-rock_physics_open/equinor_utilities/std_functions/wood_brie.py,sha256=Esb19YhjNsZLAq85fjDdzBI32RLjAZcUQeV-lA1VYlY,1929
+rock_physics_open/equinor_utilities/std_functions/backus_ave.py,sha256=KPkQrDcUL9DU4BnzcIMvgI61sJcU5aVoT62ZMnbfEpw,1997
+rock_physics_open/equinor_utilities/std_functions/dvorkin_nur.py,sha256=KgKe_xHHFRq-eWWBn_506ubGD0O_RmbF2NApWYXNqUM,2435
+rock_physics_open/equinor_utilities/std_functions/gassmann.py,sha256=oLWnGThD55z3nm4wyRDySeIdpTwRtRR-_kGmb5_cI4Q,4325
+rock_physics_open/equinor_utilities/std_functions/hashin_shtrikman.py,sha256=il6UbxmIskK_cObkTLSma2VPFqswPTLGI1kO8e6x-l4,6856
+rock_physics_open/equinor_utilities/std_functions/hertz_mindlin.py,sha256=Ekn4ig6YpkFM0m_Nlx5SwbP7SfT3rOVW2zaxiNQy7CI,1665
+rock_physics_open/equinor_utilities/std_functions/moduli_velocity.py,sha256=Tbmw_OSQaefQoRqOCW1MIjWRlZZVMqSi9UhkYgmwDx8,1521
+rock_physics_open/equinor_utilities/std_functions/reflection_eq.py,sha256=iXuWMGghwhUhyCQgGRm7xS_rhN4k6N042o3nDfgQzl8,3237
+rock_physics_open/equinor_utilities/std_functions/rho.py,sha256=xDIa2x82w2XtkxrgwYYeRnqzT8zWlsfWlp5ntVQpQDs,1696
+rock_physics_open/equinor_utilities/std_functions/voigt_reuss_hill.py,sha256=v1tRVYsH102rE504qJOdkjC2HCrL4FVlU4TI0X0yqKU,3803
+rock_physics_open/equinor_utilities/std_functions/walton.py,sha256=p97TGeHGgEQrrc-c21cb7__0rmlKYgBb82GE1thuwDM,1501
+rock_physics_open/equinor_utilities/std_functions/wood_brie.py,sha256=r4OVEYCVn-v7G2aQ-Yo79T4Zk1xa-Kh9Q8Ll6MJSuyI,2485
 rock_physics_open/equinor_utilities/various_utilities/Equinor_logo.gif,sha256=E9pXxAQig-6chS8Uz3rhFYcNQC4Nk0vbsja0pkZ-lEk,892
 rock_physics_open/equinor_utilities/various_utilities/Equinor_logo.ico,sha256=vXHlDI4DTCgxmm84_FB3BEMK-NORhcSR6wq01sb6QD0,19598
 rock_physics_open/equinor_utilities/various_utilities/__init__.py,sha256=UfxgTybUY_O48gou1DaL2bJYZKTViajN556J4EOjXLs,716
-rock_physics_open/equinor_utilities/various_utilities/display_result_statistics.py,sha256=sn06BZ3is2Yf519xMJc9vEHXq2DW5_ynWF2QKpK0Usk,3001
-rock_physics_open/equinor_utilities/various_utilities/gassmann_dry_mod.py,sha256=5BG6ozz2P69oBVejd8ZOekdN3T5q_XMgK_FcnvBHODc,1364
-rock_physics_open/equinor_utilities/various_utilities/gassmann_mod.py,sha256=JiQRjuFfJCIY6FuGaKs5XXeCnkhqiFyBu7_wzGYGfZg,1362
-rock_physics_open/equinor_utilities/various_utilities/gassmann_sub_mod.py,sha256=6IqLwY6gXd7aLpbke7oDjCNqsfjLAbgRErjq7FE_sbQ,1797
-rock_physics_open/equinor_utilities/various_utilities/hs_average.py,sha256=h6lSX66q5iWBjhPOR6qxJ0St0oYU5jYtqcrpa9OUKlk,1054
-rock_physics_open/equinor_utilities/various_utilities/pressure.py,sha256=_FckX63cGH1Le_ogfHfXQO2v8xWn-rBFSd_lfMEicfI,3300
-rock_physics_open/equinor_utilities/various_utilities/reflectivity.py,sha256=SoWjbUs6Xjii9OExiRDB70DiprPbfsl-sDQyI3v2GVw,3121
-rock_physics_open/equinor_utilities/various_utilities/timeshift.py,sha256=R_n6mlFFohlGqsi5-2geA4I0xaEw1IHBXzuNlT22v-g,2340
-rock_physics_open/equinor_utilities/various_utilities/vp_vs_rho_set_statistics.py,sha256=zre9bchR-x_l_5f7XRArrEMtzS5aIf6vAVXV9xUTY98,5105
-rock_physics_open/equinor_utilities/various_utilities/vrh_3_min.py,sha256=otMmJuiKWgMTNXYRmgjJMYH-6ND9-YQo9UzsuNVsP8w,1985
+rock_physics_open/equinor_utilities/various_utilities/display_result_statistics.py,sha256=8SfCrLhC9nmCTO2NWlEzg2EFJ76lNbA5OXtK15TYs0w,3245
+rock_physics_open/equinor_utilities/various_utilities/gassmann_dry_mod.py,sha256=l-7Pfw3wgDkVSWAcMNu_CPocRoUHMIghRVwIXL6qCAs,1757
+rock_physics_open/equinor_utilities/various_utilities/gassmann_mod.py,sha256=ZgummAFq9gxVwVKDMfqdG-iUvUnEHdfK6b7iI52BZ3g,1755
+rock_physics_open/equinor_utilities/various_utilities/gassmann_sub_mod.py,sha256=JW9KKrLW3qa8TBqDqeFJMjgqLom1pw86qMV-6oUnShw,2209
+rock_physics_open/equinor_utilities/various_utilities/hs_average.py,sha256=2lArjOpwAmZKNuBIwQXppZ7zUo2JlrwzyMdUkQPhI7A,1522
+rock_physics_open/equinor_utilities/various_utilities/pressure.py,sha256=XYw6zu2Lf9VrO4h4i8xRy2guMQSZGR0N1UnIfyhkvDI,3491
+rock_physics_open/equinor_utilities/various_utilities/reflectivity.py,sha256=kxcmooLicw9pGv_TF-O7_uKC44yAqqc0WbcYA4mvpx4,3658
+rock_physics_open/equinor_utilities/various_utilities/timeshift.py,sha256=kX7E4hnusHh5W3aCkaUAezZHn7BRQzVqPZzv2uBpVvU,2709
+rock_physics_open/equinor_utilities/various_utilities/vp_vs_rho_set_statistics.py,sha256=6M3V57icSKS2yKt5Wi8c6EIuGOoxyfSVt0c9Z5ny14c,6235
+rock_physics_open/equinor_utilities/various_utilities/vrh_3_min.py,sha256=c-hnf1FLgpLjFcys_9WRkVHelQ1VYz1Gcf0b8-2Gxt4,2592
 rock_physics_open/fluid_models/__init__.py,sha256=H3jeo69HuZCiuZ_meQe27NXIKTMOq2qiZ_UXKYyFmQs,201
 rock_physics_open/fluid_models/brine_model/__init__.py,sha256=PqVXs5rI421j7xL-Z5pFfNIdC0cUameBJHLOUqEO7wo,86
 rock_physics_open/fluid_models/brine_model/brine_properties.py,sha256=lBnai7l0gr_XGqVIaIPT-WWSudiMfBJAO-_v1FH99uo,5756
@@ -74,11 +72,11 @@ rock_physics_open/sandstone_models/constant_cement_models.py,sha256=kHRcPEC114uo
 rock_physics_open/sandstone_models/constant_cement_optimisation.py,sha256=D9ej8PaQ5iVMaCf0LKDQGgJ6CSl8JEFIaxobwddzs_Y,3745
 rock_physics_open/sandstone_models/contact_cement_model.py,sha256=EVCicRGnW13QBWuvyscFYLy8JpgANPla7ECDeZLfhdA,4488
 rock_physics_open/sandstone_models/curvefit_sandstone_models.py,sha256=pZ7u7bVNYG254fwtyyrD1SMIFn2onHBu_XZCfVfajzI,3988
-rock_physics_open/sandstone_models/friable_models.py,sha256=GmcdgT_yYhajIMEszZbuqeGb2JyQChi8o7NoctJc3aI,6091
+rock_physics_open/sandstone_models/friable_models.py,sha256=cvfLt7GQHn8EKKUk2Z86l-2o-_nXwDXDkQ1ekPQ6L6s,6153
 rock_physics_open/sandstone_models/friable_optimisation.py,sha256=YCg20EJQaQevEmilPWCbuqZ7GycKX84nmlca5W-jfRM,3467
 rock_physics_open/sandstone_models/friable_shalysand_sandyshale_models.py,sha256=eMXRCrGdCC9VyYwjg498M0GTAMlWar18iz5gS7Rz16w,6884
 rock_physics_open/sandstone_models/patchy_cement_fluid_substitution_model.py,sha256=21ujg7K-f7vykKxYwFynPenZLVZkdBAiAEwt4-K4r5I,16123
-rock_physics_open/sandstone_models/patchy_cement_model.py,sha256=EOF5kRuH9QXeVvrFps07CJWMPfWz4N2UoZTutPUPK5Y,10504
+rock_physics_open/sandstone_models/patchy_cement_model.py,sha256=zKdAAFJXgblCpmihRiEz_8DPd40nniJW97-0daCZoq0,11161
 rock_physics_open/sandstone_models/patchy_cement_optimisation.py,sha256=Tf3jRT0c1r3IgVPUFWTDqh4-XB5381dqeH7CPE27QYg,7017
 rock_physics_open/sandstone_models/unresolved_cemented_sandshale_models.py,sha256=xFF_FrQw_3DoYADK7-3Lecr_KnurNQLsXy-22eN-92c,3901
 rock_physics_open/sandstone_models/unresolved_friable_sandshale_models.py,sha256=HpRBIQw2POoQbJ2DlwoyJSDvKvhJQyhgmkjoGQFo1vE,3671
@@ -138,10 +136,10 @@ rock_physics_open/ternary_plots/__init__.py,sha256=mrmn0T4G91UPdCTEOr4HSEd1d8ZCR
 rock_physics_open/ternary_plots/gen_ternary_plot.py,sha256=uLwoeV19Tfj6Cy18NWcSoOV18DcYBJUYofb8zEYmSzM,1836
 rock_physics_open/ternary_plots/shale_prop_ternary.py,sha256=5-NB4LoXRP7JPQYaJBZwD5mnBLMjisCTEjj2u6u_boE,8228
 rock_physics_open/ternary_plots/ternary_patches.py,sha256=oYaIkCbzC4L19VWBz5iphVdQnpHX9kdy0p2NK65aIIM,7256
-rock_physics_open/ternary_plots/ternary_plot_utilities.py,sha256=weFDWyzd7wgr3xB63Q-xtT-KoJA6FNIDWuTkA7jr94Q,5669
+rock_physics_open/ternary_plots/ternary_plot_utilities.py,sha256=7CM4BuxwUjfPAGLDROBLFs1oiZ48xIZYIrgG46jlM-c,5678
 rock_physics_open/ternary_plots/unconventionals_ternary.py,sha256=KWjcBrjI8Tgr2FYpz1ey-DPuhUeSC7b5mpN8OzwNa4I,1823
-rock_physics_open-0.3.0.dist-info/licenses/LICENSE,sha256=46mU2C5kSwOnkqkw9XQAJlhBL2JAf1_uCD8lVcXyMRg,7652
-rock_physics_open-0.3.0.dist-info/METADATA,sha256=31O7aQyYTnFWoMjWrUpnxurqmsEVz16zwepO9_APppI,3901
-rock_physics_open-0.3.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-rock_physics_open-0.3.0.dist-info/top_level.txt,sha256=zFBbfbv2pshFJvcvcLLigDmlJjjg08iDdb6ABgncaP8,18
-rock_physics_open-0.3.0.dist-info/RECORD,,
+rock_physics_open-0.3.1.dist-info/licenses/LICENSE,sha256=46mU2C5kSwOnkqkw9XQAJlhBL2JAf1_uCD8lVcXyMRg,7652
+rock_physics_open-0.3.1.dist-info/METADATA,sha256=nT1qZefrmIhvMwB4iAF-_yLXV1u9P7OJrivLHxVmzlg,3901
+rock_physics_open-0.3.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+rock_physics_open-0.3.1.dist-info/top_level.txt,sha256=zFBbfbv2pshFJvcvcLLigDmlJjjg08iDdb6ABgncaP8,18
+rock_physics_open-0.3.1.dist-info/RECORD,,