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

rock_physics_open/fluid_models/gas_model/gas_properties.py

@@ -4,7 +4,7 @@ from scipy.constants import gas_constant
 
 from rock_physics_open.equinor_utilities.conversions import celsius_to_kelvin
 
-AIR_WEIGHT = 28.8  # g/mol
+AIR_WEIGHT = 28.8 * 1.0e-3  # kg/mol
 
 
 def gas_properties(
@@ -20,10 +20,11 @@ def gas_properties(
     :param model: for future use
     :return: vel_gas [m/s], den_gas [kg/m^3], k_gas [Pa], eta_gas [cP]
     """
-    den_gas = gas_density(celsius_to_kelvin(temperature), pressure * 1e-6, gas_gravity)
-    k_gas = gas_bulk_modulus(
-        celsius_to_kelvin(temperature), pressure * 1e-6, gas_gravity
-    )
+
+    den_gas = gas_density(celsius_to_kelvin(temperature), pressure, gas_gravity)
+
+    k_gas = gas_bulk_modulus(celsius_to_kelvin(temperature), pressure, gas_gravity)
+
     vel_gas = (k_gas / den_gas) ** 0.5
 
     eta_gas = lee_gas_viscosity(celsius_to_kelvin(temperature), pressure, gas_gravity)
@@ -33,9 +34,9 @@ def gas_properties(
 
 def molecular_weight(gas_gravity: np.ndarray | float) -> np.ndarray | float:
     """
-    calculates molecluar weight of a gas from gas gravity.
+    calculates molecular weight of a gas from gas gravity.
     :param gas_gravity: molar mass of gas relative to air molar mas.
-    :return: The volume of the gas in g/mol.
+    :return: The volume of the gas in kg/mol.
     """
     return gas_gravity * AIR_WEIGHT
 
@@ -47,9 +48,10 @@ def molar_volume(
     """
     calculates molar volume using the ideal gas law.
     :param absolute_temperature: The absolute temperature of the gas in kelvin.
-    :param pressure: Confining pressure in MPa.
-    :return: The volume of the gas in cc/mol.
+    :param pressure: Confining pressure in Pa.
+    :return: The volume of the gas in m^3/mol.
     """
+
     return gas_constant * absolute_temperature / pressure
 
 
@@ -62,8 +64,8 @@ def ideal_gas_density(
     calculates molar volume using the ideal gas law.
     :param gas_gravity: molar mass of gas relative to air molar mas.
     :param absolute_temperature: The absolute temperature of the gas in kelvin.
-    :param pressure: Confining pressure in MPa.
-    :return: The density of the gas in g/cc
+    :param pressure: Confining pressure in Pa.
+    :return: The density of the gas in kg/m^3
     """
     return molecular_weight(gas_gravity) / molar_volume(absolute_temperature, pressure)
 
@@ -89,13 +91,11 @@ def ideal_gas(
     :param gas_gravity: molar mass of gas relative to air molar mas.
     :param absolute_temperature: The absolute temperature of the gas in kelvin.
     :param pressure: Confining pressure in Pa.
-    :return: ideal_gas_density, ideal_gas_velocity
+    :return: ideal_gas_velocity [m/s], ideal_gas_density [kg/m^3],
     """
-    ideal_gas_den = 1000 * ideal_gas_density(
-        absolute_temperature, pressure * 1e6, gas_gravity
-    )
+    ideal_gas_den = ideal_gas_density(absolute_temperature, pressure, gas_gravity)
     ideal_gas_vel = ideal_gas_primary_velocity(absolute_temperature, gas_gravity)
-    return ideal_gas_den, ideal_gas_vel
+    return ideal_gas_vel, ideal_gas_den
 
 
 def pseudoreduced_temperature(
@@ -132,13 +132,13 @@ def pseudoreduced_pressure(
     Tech., 22, 889-892.
 
     :param gas_gravity: molar mass of gas relative to air molar mas.
-    :param pressure: Confining pressure in MPa.
-    :return: Pseudoreduced pressure in MPa.
+    :param pressure: Confining pressure in Pa.
+    :return: Pseudoreduced pressure in Pa.
     """
     return pressure / (4.892 - 0.4048 * gas_gravity)
 
 
-def compressability_factor(
+def compressibility_factor(
     absolute_temperature: np.ndarray | float,
     pressure: np.ndarray | float,
     gas_gravity: np.ndarray | float,
@@ -149,11 +149,13 @@ def compressability_factor(
 
     :param gas_gravity: molar mass of gas relative to air molar mas.
     :param absolute_temperature: The absolute temperature of the gas in kelvin.
-    :param pressure: Confining pressure in MPa.
-    :return: The density of the gas in g/cc
+    :param pressure: Confining pressure in Pa.
+    :return: Gas compressibility - unitless
     """
     tpr = pseudoreduced_temperature(absolute_temperature, gas_gravity)
-    ppr = pseudoreduced_pressure(pressure, gas_gravity)
+
+    # Pseudoreduced pressure has unit MPa in equation
+    ppr = pseudoreduced_pressure(pressure, gas_gravity) * 1.0e-6
 
     return (
         (0.03 + 0.00527 * (3.5 - tpr) ** 3) * ppr
@@ -162,7 +164,7 @@ def compressability_factor(
         - 0.52
         + 0.109
         * (3.85 - tpr) ** 2
-        / exp((0.45 + 8 * (0.56 - 1 / tpr) ** 2) * ppr**1.2 / tpr)
+        / exp((0.45 + 8.0 * (0.56 - 1 / tpr) ** 2) * ppr**1.2 / tpr)
     )
 
 
@@ -176,18 +178,17 @@ def gas_density(
 
     :param gas_gravity: molar mass of gas relative to air molar mas.
     :param absolute_temperature: The absolute temperature of the gas in kelvin.
-    :param pressure: Confining pressure in MPa.
-    :return: The density of the gas in g/cc
+    :param pressure: Confining pressure in Pa.
+    :return: The density of the gas in kg/m^3
     """
-    ideal_gas_den, ideal_gas_vel = ideal_gas(
-        absolute_temperature, pressure * 1e-6, gas_gravity
-    )
-    return ideal_gas_den / compressability_factor(
+
+    _, ideal_gas_den = ideal_gas(absolute_temperature, pressure, gas_gravity)
+    return ideal_gas_den / compressibility_factor(
         absolute_temperature, pressure, gas_gravity
     )
 
 
-def compressability_rate_per_pseudoreduced_pressure(
+def compressibility_rate_per_pseudoreduced_pressure(
     absolute_temperature: np.ndarray | float,
     pressure: np.ndarray | float,
     gas_gravity: np.ndarray | float,
@@ -198,10 +199,12 @@ def compressability_rate_per_pseudoreduced_pressure(
     :param gas_gravity: molar mass of gas relative to air molar mas.
     :param absolute_temperature: The absolute temperature of the gas in kelvin.
     :param pressure: Confining pressure in MPa.
-    :return: The density of the gas in g/cc
+    :return: Derivative of the compressibility factor (unitless) with respect to pseudoreduced pressure
     """
     tpr = pseudoreduced_temperature(absolute_temperature, gas_gravity)
-    ppr = pseudoreduced_pressure(pressure, gas_gravity)
+
+    # Pseudoreduced pressure is expected to be in MPa in the expression
+    ppr = pseudoreduced_pressure(pressure, gas_gravity) * 1.0e-6
 
     return (
         0.03
@@ -226,15 +229,16 @@ def gas_bulk_modulus(
 
     :param gas_gravity: molar mass of gas relative to air molar mas.
     :param absolute_temperature: The absolute temperature of the gas in kelvin.
-    :param pressure: Confining pressure in MPa.
-    :return: The bulk modulus of the gas in MPa.
+    :param pressure: Confining pressure in Pa.
+    :return: The bulk modulus of the gas in Pa.
     """
-    z = compressability_factor(absolute_temperature, pressure, gas_gravity)
-    dz_dppr = compressability_rate_per_pseudoreduced_pressure(
+    z = compressibility_factor(absolute_temperature, pressure, gas_gravity)
+    dz_dppr = compressibility_rate_per_pseudoreduced_pressure(
         absolute_temperature, pressure, gas_gravity
     )
 
-    ppr = pseudoreduced_pressure(pressure, gas_gravity)
+    # Set ppr in unit MPa in order to use it in calculation of gamma_0
+    ppr = pseudoreduced_pressure(pressure, gas_gravity) * 1.0e-6
 
     # Equation 11b
     gamma_0 = (
@@ -247,6 +251,44 @@ def gas_bulk_modulus(
     return gamma_0 * pressure / (1 - dz_dppr * ppr / z)
 
 
+def gas_viscosity(
+    absolute_temperature: np.ndarray | float,
+    pressure: np.ndarray | float,
+    gas_gravity: np.ndarray | float,
+) -> np.ndarray | float:
+    """
+    The gas viscosity of hydrocarbon gas, using equations 12 and 13 of Batzle & Wang [1].
+
+    :param absolute_temperature: The absolute temperature of the gas in kelvin.
+    :param pressure: Confining pressure in Pa.
+    :param gas_gravity: molar mass of gas relative to air mas.
+    :return: The gas viscosity of the gas in cP.
+    """
+    temp_pr = pseudoreduced_temperature(absolute_temperature, gas_gravity)
+
+    # Pseudoreduced pressure should be in unit MPa
+    pres_pr = pseudoreduced_pressure(pressure, gas_gravity) * 1.0e-6
+
+    eta_1 = 0.0001 * (
+        temp_pr * (28.0 + 48.0 * gas_gravity - 5.0 * gas_gravity**2)
+        - 6.47 * gas_gravity**-2
+        + 35.0 * gas_gravity**-1
+        + 1.14 * gas_gravity
+        - 15.55
+    )
+    return eta_1 * (
+        0.001
+        * pres_pr
+        * (
+            (1057.0 - 8.08 * temp_pr) / pres_pr
+            + (796.0 * pres_pr**0.5 - 704.0)
+            / (((temp_pr - 1.0) ** 0.7) * (pres_pr + 1.0))
+            - 3.24 * temp_pr
+            - 38.0
+        )
+    )
+
+
 def lee_gas_viscosity(
     absolute_temperature: np.ndarray | float,
     pressure: np.ndarray | float,

rock_physics_open/fluid_models/oil_model/dead_oil_density.py

@@ -10,17 +10,19 @@ def pressure_adjusted_dead_oil_density(
 
     Uses equation 18 from Batzle & Wang [1].
 
-    :param reference_density: The density (g/cc) of the dead oil at 15.6 degrees Celsius
+    :param reference_density: The density [kg/m^3] of the dead oil at 15.6 degrees Celsius
         and atmospheric pressure.
-    :param pressure: Pressure (MPa) to adjust to.
+    :param pressure: Pressure [Pa] to adjust to.
     :return: Density of oil at given pressure and 21 degrees Celsius (~70 degrees
-    Farenheit).
+    Farenheit). [kg/m^3]
     """
-    return (
-        reference_density
-        + (0.00277 * pressure - 1.71 * 10**-7 * pressure**3)
-        * (reference_density - 1.15) ** 2
-        + 3.49 * 10**-4 * pressure
+    pressure_mpa = pressure / 1e6
+    density_gcc = reference_density / 1000.0
+    return 1000.0 * (
+        density_gcc
+        + (0.00277 * pressure_mpa - 1.71e-7 * pressure_mpa**3)
+        * (density_gcc - 1.15) ** 2
+        + 3.49e-4 * pressure_mpa
     )
 
 
@@ -33,11 +35,14 @@ def temperature_adjusted_dead_oil_density(
 
     Uses equation 19 from Batzle & Wang [1].
 
-    :param density_at_21c: The density (g/cc) of the dead oil at 21 degrees Celsius
-    :param temperature: Temperature (Celsius) of oil.
-    :return: Density of oil at given temperature.
+    :param density_at_21c: The density [kg/m^3] of the dead oil at 21 degrees Celsius
+    :param temperature: Temperature [°C] of oil.
+    :return: Density of oil at given temperature. [kg/m^3]
     """
-    return density_at_21c / (0.972 + 3.81 * (10**-4) * (temperature + 17.78) ** 1.175)
+    density_at_21c_gcc = density_at_21c / 1000.0
+    return (
+        1000.0 * density_at_21c_gcc / (0.972 + 3.81e-4 * (temperature + 17.78) ** 1.175)
+    )
 
 
 def dead_oil_density(
@@ -51,10 +56,10 @@ def dead_oil_density(
     Uses equation 18 & 19 from Batzle & Wang [1].
 
     :param reference_density: Density of oil at 15.6 degrees Celsius and atmospheric
-        pressure. (g/cc)
-    :param pressure: Pressure (MPa) of oil
-    :param temperature: Temperature (Celsius) of oil.
-    :return: density of dead oil at given conditions.
+        pressure [kg/m^3]
+    :param pressure: Pressure [Pa] of oil
+    :param temperature: Temperature [°C] of oil.
+    :return: density of dead oil at given conditions (kg/m^3).
     """
     density_p = pressure_adjusted_dead_oil_density(pressure, reference_density)
     return temperature_adjusted_dead_oil_density(temperature, density_p)

rock_physics_open/fluid_models/oil_model/dead_oil_velocity.py

@@ -12,17 +12,19 @@ def dead_oil_velocity(
     Uses equation 20a from Batzle & Wang [1].
 
     :param reference_density: Density of oil at 15.6 degrees Celsius and atmospheric
-        pressure (g/cc)
-    :param pressure: Pressure (MPa) of oil
-    :param temperature: Temperature (Celsius) of oil.
+        pressure [kg/m^3]
+    :param pressure: Pressure [Pa] of oil
+    :param temperature: Temperature [°C] of oil.
     :return: primary velocity of dead oil in m/s.
     """
+    pressure_mpa = pressure * 1e-6
+    density_gcc = reference_density / 1000.0
     return (
-        2096 * np.sqrt(reference_density / (2.6 - reference_density))
+        2096 * np.sqrt(density_gcc / (2.6 - density_gcc))
         - 3.7 * temperature
-        + 4.64 * pressure
+        + 4.64 * pressure_mpa
         + 0.0115
-        * (4.12 * np.sqrt(1.08 * reference_density**-1 - 1) - 1)
+        * (4.12 * np.sqrt(1.08 * density_gcc**-1 - 1) - 1)
         * temperature
-        * pressure
+        * pressure_mpa
     )

rock_physics_open/fluid_models/oil_model/live_oil_density.py

@@ -14,17 +14,18 @@ def live_oil_density(
     Equation 24 in Batzle & Wang [1].
 
     :param reference_density: Density of the oil without dissolved gas
-        at 15.6 degrees Celsius and atmospheric pressure. (g/cc)
-    :param pressure: Pressure (MPa) of oil (for future implementation only)
+        at 15.6 degrees Celsius and atmospheric pressure. [kg/m^3]
+    :param pressure: Pressure [Pa] of oil (for future implementation only)
     :param gas_oil_ratio: The volume ratio of gas to oil [l/l]
-    :param temperature: Temperature (Celsius) of oil.
+    :param temperature: Temperature [°C] of oil.
     :param gas_gravity: molar mass of gas relative to air molar mas.
-    :return: Density of live oil [g/cc].
+    :return: Density of live oil [kg/m^3].
     """
+    density_gcc = reference_density / 1000.0
     b0 = live_oil_volume_factor(
         temperature, reference_density, gas_oil_ratio, gas_gravity
     )
-    return (reference_density + 0.0012 * gas_gravity * gas_oil_ratio) / b0
+    return 1000.0 * (density_gcc + 0.0012 * gas_gravity * gas_oil_ratio) / b0
 
 
 def live_oil_pseudo_density(
@@ -40,16 +41,17 @@ def live_oil_pseudo_density(
     Equation 22 in Batzle & Wang [1].
 
     :param reference_density: Density of the oil without dissolved gas
-        at 15.6 degrees Celsius and atmospheric pressure. (g/cc)
+        at 15.6 degrees Celsius and atmospheric pressure. [kg/m^3]
     :param gas_oil_ratio: The volume ratio of gas to oil [l/l]
-    :param temperature: Temperature (Celsius) of oil.
+    :param temperature: Temperature [°C] of oil.
     :param gas_gravity: molar mass of gas relative to air molar mas.
-    :return: Pseudo-density of live oil.
+    :return: Pseudo-density of live oil [kg/m^3].
     """
+    density_gcc = reference_density / 1000.0
     b0 = live_oil_volume_factor(
         temperature, reference_density, gas_oil_ratio, gas_gravity
     )
-    return (reference_density / b0) / (1 + 0.001 * gas_oil_ratio)
+    return 1000.0 * (density_gcc / b0) / (1 + 0.001 * gas_oil_ratio)
 
 
 def live_oil_volume_factor(
@@ -61,17 +63,18 @@ def live_oil_volume_factor(
     """
     Volume factor derived by Standing (1962), equation 23 in Batzle & Wang [1].
     :param reference_density: Density of the oil without dissolved gas
-        at 15.6 degrees Celsius and atmospheric pressure. (g/cc)
+        at 15.6 degrees Celsius and atmospheric pressure. [kg/m^3]
     :param gas_oil_ratio: The volume ratio of gas to oil [l/l]
-    :param temperature: Temperature (Celsius) of oil.
+    :param temperature: Temperature [°C] of oil.
     :param gas_gravity: molar mass of gas relative to air molar mas.
-    :return: A volume factor in calculating pseudo-density of live oil.
+    :return: A volume factor in calculating pseudo-density of live oil [unitless].
     """
+    density_gcc = reference_density / 1000.0
     return (
         0.972
         + 0.00038
         * (
-            2.4 * gas_oil_ratio * np.sqrt(gas_gravity / reference_density)
+            2.4 * gas_oil_ratio * np.sqrt(gas_gravity / density_gcc)
             + temperature
             + 17.8
         )

rock_physics_open/fluid_models/oil_model/live_oil_velocity.py

@@ -11,10 +11,10 @@ def live_oil_velocity(
     Substitute Equation 22 in Equation 20 of Batzle & Wang [1].
 
     :param reference_density: Density of the oil without dissolved gas
-        at 15.6 degrees Celsius and atmospheric pressure. (g/cc)
-    :param pressure: Pressure (MPa) of oil
+        at 15.6 degrees Celsius and atmospheric pressure. [kg/m^3]
+    :param pressure: Pressure [Pa] of oil
     :param gas_oil_ratio: The volume ratio of gas to oil [l/l]
-    :param temperature: Temperature (Celsius) of oil.
+    :param temperature: Temperature [°C] of oil.
     :param gas_gravity: molar mass of gas relative to air molar mas.
     :return: Primary wave velocity of live oil [m/s].
     """

rock_physics_open/fluid_models/oil_model/oil_properties.py

@@ -17,13 +17,13 @@ def oil_properties(
     gas_gravity: np.ndarray | float,
 ) -> np.ndarray | float:
     """
-    :param temperature: Temperature (Celsius) of oil.
-    :param pressure: Pressure (Pa) of oil
+    :param temperature: Temperature [°C] of oil.
+    :param pressure: Pressure [Pa] of oil
     :param rho0: Density of the oil without dissolved gas at 15.6 degrees Celsius and
-                 atmospheric pressure. (kg/m^3)
+                 atmospheric pressure. [kg/m^3]
     :param gas_oil_ratio: The volume ratio of gas to oil [l/l]
     :param gas_gravity: Gas Gravity, molar mass of gas relative to air molar mas.
-    :return: vel_oil, den_oil, k_oil
+    :return: vel_oil [m/s], den_oil [kg/m^3], k_oil [Pa]
     """
     # Since live_oil with gas_oil_ratio=0.0 is not equal to dead oil
     # we use an apodization function to interpolate between the two
@@ -41,10 +41,11 @@ def oil_properties(
         window = np.clip((np.abs(x) - length / 2) / (length / 2), 0, 1)
         return 1 - window
 
-    loil_den, loil_vel = live_oil(
-        temperature, pressure, rho0, gas_oil_ratio, gas_gravity
-    )
-    doil_den, doil_vel = dead_oil(temperature, pressure, rho0)
+    (
+        loil_vel,
+        loil_den,
+    ) = live_oil(temperature, pressure, rho0, gas_oil_ratio, gas_gravity)
+    doil_vel, doil_den = dead_oil(temperature, pressure, rho0)
     window = triangular_window(gas_oil_ratio)
     den_oil = doil_den * window + (1 - window) * loil_den
     vel_oil = doil_vel * window + (1 - window) * loil_vel
@@ -59,20 +60,16 @@ def dead_oil(
 ) -> tuple[np.ndarray | float, np.ndarray | float]:
     """
     :param reference_density: Density of the oil without dissolved gas
-        at 15.6 degrees Celsius and atmospheric pressure. kg/m3
+        at 15.6 degrees Celsius and atmospheric pressure. [kg/m^3]
     :param gas_oil_ratio: The volume ratio of gas to oil [l/l]
     :param gas_gravity: molar mass of gas relative to air molar mas.
-    :param pressure: Pressure (Pa) of oil
-    :param temperature: Temperature (Celsius) of oil.
-    :return: dead_oil_density, dead_oil_velocity
+    :param pressure: Pressure [Pa] of oil
+    :param temperature: Temperature [°C] of oil.
+    :return: dead_oil_density [kg/m^3], dead_oil_velocity [m/s]
     """
-    dead_oil_den = 1000 * dead_oil_density(
-        temperature, pressure * 1e-6, reference_density / 1000
-    )
-    dead_oil_vel = dead_oil_velocity(
-        temperature, pressure * 1e-6, reference_density / 1000
-    )
-    return dead_oil_den, dead_oil_vel
+    dead_oil_den = dead_oil_density(temperature, pressure, reference_density)
+    dead_oil_vel = dead_oil_velocity(temperature, pressure, reference_density)
+    return dead_oil_vel, dead_oil_den
 
 
 def live_oil(
@@ -84,12 +81,12 @@ def live_oil(
 ) -> tuple[np.ndarray | float, np.ndarray | float]:
     """
     :param reference_density: Density of the oil without dissolved gas
-        at 15.6 degrees Celsius and atmospheric pressure. (kg/m^3)
+        at 15.6 degrees Celsius and atmospheric pressure. [kg/m^3]
     :param gas_oil_ratio: The volume ratio of gas to oil [l/l]
     :param gas_gravity: molar mass of gas relative to air molar mas.
-    :param pressure: Pressure (Pa) of oil
-    :param temperature: Temperature (Celsius) of oil.
-    :return: live_oil_density, live_oil_velocity
+    :param pressure: Pressure [Pa] of oil
+    :param temperature: Temperature [°C] of oil.
+    :return: live_oil_density , live_oil_velocity
     """
     if np.any(
         pressure
@@ -99,18 +96,51 @@ def live_oil(
             "Pressure is below bubble point of oil, estimated elastic properties can be inaccurate",
             stacklevel=1,
         )
-    live_oil_den = 1000 * live_oil_density(
+    live_oil_den = live_oil_density(
         temperature,
-        pressure * 1e-6,
-        reference_density / 1000,
+        pressure,
+        reference_density,
         gas_oil_ratio,
         gas_gravity,
     )
     live_oil_vel = live_oil_velocity(
         temperature,
-        pressure * 1e-6,
-        reference_density / 1000,
+        pressure,
+        reference_density,
         gas_oil_ratio,
         gas_gravity,
     )
-    return live_oil_den, live_oil_vel
+    return (
+        live_oil_vel,
+        live_oil_den,
+    )
+
+
+def oil_viscosity(
+    temperature: np.ndarray | float,
+    pressure: np.ndarray | float,
+    reference_density: np.ndarray | float,
+) -> np.ndarray | float:
+    """
+    Calculate dead oil viscosity. If dissolved gas is present in the oil, the reference density
+    should be substituted by live oil density.
+
+    Equations 25a, 25b, 26a & 26b in Batzle and Wang 1992
+
+    Based on Beggs and Robinson 1975
+
+    :param temperature: Temperature [°C] of oil
+    :param pressure: Pressure [Pa] of oil
+    :param reference_density: Density of the oil without dissolved gas
+    """
+    # Change unit in pressure to MPa
+    pressure_mpa = pressure / 1.0e6
+    # Change unit in density to g/cc
+    density_gcc = reference_density / 1000.0
+
+    y_factor = 10 ** (5.693 - 2.863 / density_gcc)
+    eta_t = -1.0 + 10 ** (0.505 * y_factor * (17.8 + temperature) ** -1.163)
+    i_factor = 10 ** (
+        18.6 * (0.1 * np.log10(eta_t) + (np.log10(eta_t) + 2) ** -0.1 - 0.985)
+    )
+    return eta_t + 0.145 * pressure_mpa * i_factor

rock_physics_open/sandstone_models/__init__.py

@@ -22,6 +22,7 @@ from .patchy_cement_fluid_substitution_model import (
 from .patchy_cement_model import (
     constant_cement_model_pcm,
     patchy_cement_model_cem_frac,
+    patchy_cement_model_dry,
     patchy_cement_model_weight,
 )
 from .patchy_cement_optimisation import (
@@ -48,6 +49,7 @@ __all__ = [
     "friable_shaly_sand_sandy_shale_model",
     "patchy_cement_pressure_fluid_substitution",
     "constant_cement_model_pcm",
+    "patchy_cement_model_dry",
     "patchy_cement_model_cem_frac",
     "patchy_cement_model_weight",
     "patchy_cement_model_optimisation",

rock_physics_open/sandstone_models/constant_cement_optimisation.py

@@ -1,7 +1,10 @@
 import numpy as np
 
 from rock_physics_open.equinor_utilities import gen_utilities
-from rock_physics_open.t_matrix_models import gen_opt_routine, save_opt_params
+from rock_physics_open.equinor_utilities.optimisation_utilities import (
+    gen_opt_routine,
+    save_opt_params,
+)
 
 from .curvefit_sandstone_models import curvefit_constant_cement
 

rock_physics_open/sandstone_models/friable_optimisation.py

@@ -1,7 +1,10 @@
 import numpy as np
 
 from rock_physics_open.equinor_utilities import gen_utilities
-from rock_physics_open.t_matrix_models import gen_opt_routine, save_opt_params
+from rock_physics_open.equinor_utilities.optimisation_utilities import (
+    gen_opt_routine,
+    save_opt_params,
+)
 
 from .curvefit_sandstone_models import curvefit_friable