rock-physics-open 0.0__py3-none-any.whl

rock_physics_open/t_matrix_models/t_matrix_opt_fluid_sub_petec.py

@@ -0,0 +1,163 @@
+import numpy as np
+
+from rock_physics_open.equinor_utilities import gen_utilities
+
+from .curvefit_t_matrix_min import curve_fit_2_inclusion_sets
+from .opt_subst_utilities import gen_mod_routine, gen_sub_routine, load_opt_params
+
+
+def run_t_matrix_with_opt_params_petec(
+    min_k,
+    min_mu,
+    min_rho,
+    fl_k_orig,
+    fl_rho_orig,
+    fl_k_sub,
+    fl_rho_sub,
+    vp,
+    vs,
+    rhob,
+    phi,
+    angle,
+    perm,
+    visco,
+    tau,
+    freq,
+    f_name,
+    fluid_sub=True,
+):
+    """
+        Based on the input file with parameters for the optimally fitted model, the correct modelling version is run.
+        Fluid substitution follows, in case it is selected. If not, the vp_sub and vs_sub will contain the same values as
+        the input logs.
+    def run_t_matrix_with_opt_params_exp(fl_k_orig, fl_rho_orig, fl_k_sub, fl_rho_sub,
+                                     vp, vs, rhob,
+                                     phi, vsh,
+                                     angle, perm, visco, tau, freq,
+                                     f_name, fluid_sub=True):
+        Parameters
+        ----------
+        min_k : np.ndarray
+            Effective mineral bulk modulus [Pa].
+        min_mu : np.ndarray
+            Effective mineral shear modulus [Pa].
+        min_rho : np.ndarray
+            Effective mineral density [kg/m^3].
+        fl_k_orig : np.ndarray
+            Effective in situ fluid bulk modulus [Pa].
+        fl_rho_orig : np.ndarray
+            Effective in situ fluid density [kg/m^3].
+        fl_k_sub : np.ndarray
+            Effective substituted fluid bulk modulus [Pa].
+        fl_rho_sub : Effective substituted density [kg/m^3].
+        vp : np.ndarray
+            Compressional velocity [m/s].
+        vs : np.ndarray
+            Shear velocity [m/s].
+        rhob : np.ndarray
+            Bulk density [kg/m^3].
+        phi : np.ndarray
+            Porosity [fraction].
+        angle : float
+            Angle of symmetry plane [degrees]
+        perm : float
+            Permeability [mD].
+        visco : float
+            Viscosity [cP].
+        tau : float
+            Relaxation time constant [s].
+        freq : float
+            Signal frequency [Hz].
+        f_name : str
+            File name for parameter file for optimal parameters.
+        fluid_sub : bool
+            Boolean parameter to perform fluid substitution.
+
+        Returns
+        -------
+        tuple
+            Tuple of np.ndarrays: vp and vs for pressure substituted case, vp, vs and density for fluid substituted case, vp and vs for
+            optimal fitted model, vp and vs residuals (observed logs minus modelled values).
+    """
+    opt_type, opt_params, opt_dict = load_opt_params(f_name)
+    y_data = np.stack([vp, vs], axis=1)
+    y_shape = y_data.shape
+    phi, angle, perm, visco, tau, freq, def_vpvs = gen_utilities.dim_check_vector(
+        (phi, angle, perm, visco, tau, freq, 1.0)
+    )
+
+    rho_sub = rhob + (fl_rho_sub - fl_rho_orig) * phi
+    # Set None values for inputs that will be defined in the different cases
+    x_data_new = None
+
+    opt_fcn = curve_fit_2_inclusion_sets
+    # Generate x_data according to method min
+    x_data = np.stack(
+        (
+            phi,
+            min_k,
+            min_mu,
+            min_rho,
+            fl_k_orig,
+            fl_rho_orig,
+            angle,
+            perm,
+            visco,
+            tau,
+            freq,
+            def_vpvs,
+        ),
+        axis=1,
+    )
+    if fluid_sub:
+        x_data_new = np.stack(
+            (
+                phi,
+                min_k,
+                min_mu,
+                min_rho,
+                fl_k_sub,
+                fl_rho_sub,
+                angle,
+                perm,
+                visco,
+                tau,
+                freq,
+                def_vpvs,
+            ),
+            axis=1,
+        )
+    rho_mod = min_rho * (1.0 - phi) + fl_rho_orig * phi
+
+    if fluid_sub:
+        v_sub, v_mod, v_res = gen_sub_routine(
+            opt_fcn, x_data, x_data_new, y_data, opt_params
+        )
+        vp_sub, vs_sub = [arr.flatten() for arr in np.split(v_sub, 2, axis=1)]
+        vp_mod, vs_mod = [arr.flatten() for arr in np.split(v_mod, 2, axis=1)]
+        vp_res, vs_res = [arr.flatten() for arr in np.split(v_res, 2, axis=1)]
+    else:
+        v_mod = gen_mod_routine(opt_fcn, x_data, y_shape, opt_params)
+        vp_mod, vs_mod = [arr.flatten() for arr in np.split(v_mod, 2, axis=1)]
+        vp_sub = vp
+        vs_sub = vs
+        vp_res = vp_mod - vp
+        vs_res = vs_mod - vs
+
+    rho_res = rho_mod - rhob
+    ai_sub = vp_sub * rho_sub
+    vpvs_sub = vp_sub / vs_sub
+
+    return (
+        vp_sub,
+        vs_sub,
+        rho_sub,
+        ai_sub,
+        vpvs_sub,
+        vp_mod,
+        vs_mod,
+        rho_mod,
+        vp_res,
+        vs_res,
+        rho_res,
+    )

rock_physics_open/t_matrix_models/t_matrix_opt_forward_model_exp.py

@@ -0,0 +1,72 @@
+import numpy as np
+
+from rock_physics_open.equinor_utilities import gen_utilities
+
+from .curvefit_t_matrix_exp import curvefit_t_matrix_exp
+from .opt_subst_utilities import gen_mod_routine, load_opt_params, opt_param_info
+
+
+def run_t_matrix_forward_model_with_opt_params_exp(
+    fl_k, fl_rho, phi, vsh, angle, perm, visco, tau, freq, f_name
+):
+    """Based on the input file with parameters for the optimally fitted model, a forward modelling is done
+    with inputs of mineral properties, fluid properties and porosity per sample. Other parameters (constants)
+    can also be varied from their setting when the optimal parameters were found.
+
+    Parameters
+    ----------
+    fl_k : np.ndarray.
+        Effective in situ fluid bulk modulus [Pa].
+    fl_rho : np.ndarray
+        Effective in situ fluid density [kg/m^3].
+    phi : np.ndarray
+        Porosity [fraction].
+    vsh : np.ndarray
+        Shale volume [fraction].
+    angle : float
+        Angle of symmetry plane
+    perm : float
+        Permeability [mD].
+    visco : float
+        Viscosity [cP].
+    tau : float
+        Relaxation time constant [s].
+    freq : float
+        Signal frequency [Hz].
+    f_name : str
+        File name for parameter file for optimal parameters.
+
+    Returns
+    -------
+    tuple
+        Tuple of np.ndarrays: vp [m/s], vs [m/s], rho [kg/m^3], ai [kg/m^3 x m/s], vp/vs [fraction] for forward model.
+    """
+    opt_type, opt_params, opt_dict = load_opt_params(f_name)
+    scale_val = opt_param_info()[1]
+    phi, angle, perm, visco, tau, freq, def_vpvs = gen_utilities.dim_check_vector(
+        (phi, angle, perm, visco, tau, freq, 1.0)
+    )
+    rho_mod = (
+        (1.0 - vsh) * opt_dict["rho_carb"] * scale_val["rho_carb"]
+        + vsh * opt_dict["rho_sh"] * scale_val["rho_sh"]
+    ) * (1.0 - phi) + phi * fl_rho
+    y_shape = (phi.shape[0], 2)
+
+    if opt_type != "exp":
+        raise ValueError(
+            f"{__file__}: incorrect type of optimal parameter input file, must come from EXP "
+            f"optimisation"
+        )
+    # No Need for preprocessing
+    opt_fcn = curvefit_t_matrix_exp
+    # Generate x_data according to method min
+    x_data = np.stack(
+        (phi, vsh, fl_k, fl_rho, angle, perm, visco, tau, freq, def_vpvs), axis=1
+    )
+    v_mod = gen_mod_routine(opt_fcn, x_data, y_shape, opt_params)
+    if not isinstance(v_mod, np.ndarray):
+        raise ValueError(f"{__file__}: no solution to forward model")
+    vp_mod, vs_mod = [arr.flatten() for arr in np.split(v_mod, 2, axis=1)]
+    vpvs_mod = vp_mod / vs_mod
+    ai_mod = vp_mod * rho_mod
+    return vp_mod, vs_mod, rho_mod, ai_mod, vpvs_mod

rock_physics_open/t_matrix_models/t_matrix_opt_forward_model_min.py

@@ -0,0 +1,86 @@
+import numpy as np
+
+from rock_physics_open.equinor_utilities import gen_utilities
+
+from .curvefit_t_matrix_min import curve_fit_2_inclusion_sets
+from .opt_subst_utilities import gen_mod_routine, load_opt_params
+
+
+def run_t_matrix_forward_model_with_opt_params_petec(
+    min_k, min_mu, min_rho, fl_k, fl_rho, phi, angle, perm, visco, tau, freq, f_name
+):
+    """Based on the input file with parameters for the optimally fitted model, a forward modelling is done
+    with inputs of mineral properties, fluid properties and porosity per sample. Other parameters (constants)
+    can also be varied from their setting when the optimal parameters were found.
+
+    Parameters
+    ----------
+    min_k : np.ndarray
+        Effective mineral bulk modulus [Pa].
+    min_mu : np.ndarray
+        Effective mineral shear modulus [Pa].
+    min_rho : np.ndarray
+        Effective mineral density [kg/m^3].
+    fl_k : np.ndarray
+        Effective in situ fluid bulk modulus [Pa].
+    fl_rho : np.ndarray
+        Effective in situ fluid density [kg/m^3].
+    phi : np.ndarray
+        Porosity [fraction].
+    angle : float
+        Angle of symmetry plane [degrees]
+    perm : float
+        Permeability [mD].
+    visco : float
+        Viscosity [cP].
+    tau : float
+        Relaxation time constant [s].
+    freq : float
+        Signal frequency [Hz].
+    f_name : str
+        File name for parameter file for optimal parameters.
+
+    Returns
+    -------
+    tuple
+        Tuple of np.ndarrays: vp [m/s], vs [m/s], rho [kg/m^3], ai [kg/m^3 x m/s], vp/vs [fraction] for forward model.
+    """
+    opt_type, opt_params, opt_dict = load_opt_params(f_name)
+    phi, angle, perm, visco, tau, freq, def_vpvs = gen_utilities.dim_check_vector(
+        (phi, angle, perm, visco, tau, freq, 1.0)
+    )
+    rho_mod = min_rho * (1.0 - phi) + fl_rho * phi
+    y_shape = (phi.shape[0], 2)
+
+    if opt_type != "min":
+        raise ValueError(
+            f"{__file__}: incorrect type of optimal parameter input file, must come from PETEC "
+            f"optimisation"
+        )
+    # No Need for preprocessing
+    opt_fcn = curve_fit_2_inclusion_sets
+    # Generate x_data according to method min
+    x_data = np.stack(
+        (
+            phi,
+            min_k,
+            min_mu,
+            min_rho,
+            fl_k,
+            fl_rho,
+            angle,
+            perm,
+            visco,
+            tau,
+            freq,
+            def_vpvs,
+        ),
+        axis=1,
+    )
+    v_mod = gen_mod_routine(opt_fcn, x_data, y_shape, opt_params)
+    if not isinstance(v_mod, np.ndarray):
+        raise ValueError(f"{__file__}: no solution to forward model")
+    vp_mod, vs_mod = [arr.flatten() for arr in np.split(v_mod, 2, axis=1)]
+    vpvs_mod = vp_mod / vs_mod
+    ai_mod = vp_mod * rho_mod
+    return vp_mod, vs_mod, rho_mod, ai_mod, vpvs_mod

rock_physics_open/t_matrix_models/t_matrix_parameter_optimisation_exp.py

@@ -0,0 +1,172 @@
+import numpy as np
+
+from rock_physics_open.equinor_utilities import gen_utilities
+
+from .curvefit_t_matrix_exp import curvefit_t_matrix_exp
+from .opt_subst_utilities import gen_opt_routine, opt_param_info, save_opt_params
+from .t_matrix_parameter_optimisation_min import DEF_VP_VS_RATIO
+
+
+def t_matrix_optimisation_exp(
+    k_fl: np.ndarray,
+    rho_fl: np.ndarray,
+    por: np.ndarray,
+    vsh: np.ndarray,
+    vp: np.ndarray,
+    vs: np.ndarray,
+    rhob: np.ndarray,
+    angle: float = 0.0,
+    k_r: float = 50.0,
+    eta_f: float = 1.0,
+    tau: float = 1.0e-7,
+    freq: float = 1.0e3,
+    file_out_str: str = "opt_params_exp.pkl",
+    display_results: bool = False,
+    well_name: str = "Unknown well",
+    **opt_kwargs,
+):
+    """T-Matrix optimisation adapted to an exploration setting, where detailed well information is generally not known.
+    Inclusion parameters are optimised for the whole well.
+
+    Parameters
+    ----------
+    por :
+        Inclusion porosity [ratio].
+    vsh :
+        Shale volume [ratio].
+    k_fl :
+        Fluid bulk modulus [Pa].
+    rho_fl :
+        Fluid density [kg/m^3].
+    vp :
+        Compressional velocity log [m/s].
+    vs :
+        Shear velocity log [m/s].
+    rhob :
+        Bulk density log [kg/m^3].
+    angle : float
+        Angle of symmetry plane [degrees]
+    k_r :
+        Permeability [mD].
+    eta_f :
+        Fluid viscosity [cP].
+    tau :
+        Relaxation time constant [s].
+    freq :
+        Signal frequency [Hz].
+    file_out_str :
+        Output file name (string) to store optimal parameters (pickle format).
+    display_results :
+        D isplay optimal parameters in a window after run.
+    well_name :
+        Name of well to be displayed in info box title.
+    opt_kwargs :
+        Additional keywords to be passed to optimisation function
+
+    Returns
+    -------
+    tuple
+        vp_mod, vs_mod - modelled logs, vp_res, vs_res - residual logs.
+    """
+    # 1. Preparation that is independent of search for minimum possible aspect ratio for second inclusion set
+
+    # Optimisation function for selected parameters
+    opt_fun = curvefit_t_matrix_exp
+    # rho_min = (rhob - por * rho_fl) / (1 - por)
+    # EXP adapted inputs: include fluid data and other params in x_data
+    # expand single value parameters to match logs length
+    por, angle, k_r, eta_f, tau, freq, def_vpvs = gen_utilities.dim_check_vector(
+        (por, angle, k_r, eta_f, tau, freq, DEF_VP_VS_RATIO)
+    )
+    x_data = np.stack(
+        (por, vsh, k_fl, rho_fl, angle, k_r, eta_f, tau, freq, def_vpvs), axis=1
+    )
+    # Set weight to vs to give vp and vs similar influence on optimisation
+    y_data = np.stack([vp, vs * DEF_VP_VS_RATIO], axis=1)
+
+    # 2. Search for minimum aspect ratio of inclusion set no. 2, given that inclusion set no. 1 will at least represent
+    # 50% of inclusions. Minimum aspect ratio is linked to the porosity, and the most conservative estimate is to use
+    # the maximum value. This is likely to deteriorate the optimisation results, so a search for a more appropriate
+    # value that still produces valid results is made
+    percentiles = [50, 75, 80, 85, 90, 95, 99, 100]
+    valid_result = False
+    vel_mod = None
+    vel_res = None
+    opt_params = None
+    scale_val = opt_param_info()[1]
+    while not valid_result and percentiles:
+        try:
+            # Make sure that parameters are not in conflict with T Matrix assumptions
+            min_v1 = 0.5
+            max_por = np.percentile(por, percentiles[0])
+            min_a2 = (1.0 - min_v1) * max_por
+            # Test with all parameters in the range 0.0 - 1.0
+            # Params:      f_ani f_con a1 a2 v1 k_carb mu_carb rho_carb k_sh mu_sh rho_sh
+            lower_bound = np.array(
+                [
+                    0.0,
+                    0.0,
+                    0.5,
+                    min_a2,
+                    min_v1,  # f_ani f_con a1 a2 v1
+                    35.0 / scale_val["k_carb"],
+                    30.0 / scale_val["mu_carb"],
+                    2650.0 / scale_val["rho_carb"],  # k_carb, mu_carb, rho_carb
+                    15.0 / scale_val["k_sh"],
+                    7.0 / scale_val["mu_sh"],
+                    2500.0 / scale_val["rho_sh"],
+                ],
+                dtype=float,
+            )  # k_sh, mu_sh, rho_sh
+            upper_bound = np.array(
+                [
+                    1.0,
+                    1.0,
+                    1.0,
+                    0.30,
+                    1.0,  # f_ani f_con a1 a2 v1
+                    1.0,
+                    1.0,
+                    1.0,  # k_carb, mu_carb, rho_carb
+                    1.0,
+                    1.0,
+                    1.0,
+                ],
+                dtype=float,
+            )  # k_sh, mu_sh, rho_sh
+            x0 = (upper_bound + lower_bound) / 2.0
+            # Optimisation step without fluid substitution
+            vel_mod, vel_res, opt_params = gen_opt_routine(
+                opt_fun, x_data, y_data, x0, lower_bound, upper_bound, **opt_kwargs
+            )
+            valid_result = True
+        except ValueError:
+            percentiles.pop(0)
+            valid_result = False
+
+    if not valid_result:
+        raise ValueError(
+            f"{__file__}: unable to find stable value for T Matrix optimisation, second inclusion"
+        )
+
+    # Reshape outputs and remove weight from vs
+    vp_mod, vs_mod = [arr.flatten() for arr in np.split(vel_mod, 2, axis=1)]
+    vp_res, vs_res = [arr.flatten() for arr in np.split(vel_res, 2, axis=1)]
+    vs_mod = vs_mod / DEF_VP_VS_RATIO
+    vs_res = vs_res / DEF_VP_VS_RATIO
+    vpvs_mod = vp_mod / vs_mod
+    # Calculate the modelled density
+    rhob_mod = (
+        (1.0 - vsh) * opt_params[7] * scale_val["rho_carb"]
+        + vsh * opt_params[10] * scale_val["rho_sh"]
+    ) * (1.0 - por) + por * rho_fl
+    ai_mod = vp_mod * rhob_mod
+    rhob_res = rhob_mod - rhob
+    # Save the optimal parameters
+    save_opt_params("exp", opt_params, file_out_str, well_name=well_name)
+    if display_results:
+        from .opt_subst_utilities import opt_param_to_ascii
+
+        opt_param_to_ascii(file_out_str, well_name=well_name)
+
+    return vp_mod, vs_mod, rhob_mod, ai_mod, vpvs_mod, vp_res, vs_res, rhob_res

rock_physics_open/t_matrix_models/t_matrix_parameter_optimisation_min.py

@@ -0,0 +1,159 @@
+import numpy as np
+
+from rock_physics_open.equinor_utilities import gen_utilities
+
+from .curvefit_t_matrix_min import curve_fit_2_inclusion_sets
+from .opt_subst_utilities import gen_opt_routine, save_opt_params
+
+# Trade-off between calcite, dolomite and quartz, vs is weighted by this in order to make it count as much as vp
+# in the optimisation
+DEF_VP_VS_RATIO = 1.8
+
+
+def t_matrix_optimisation_petec(
+    k_min: np.ndarray,
+    mu_min: np.ndarray,
+    rho_min: np.ndarray,
+    k_fl: np.ndarray,
+    rho_fl: np.ndarray,
+    por: np.ndarray,
+    vp: np.ndarray,
+    vs: np.ndarray,
+    rhob: np.ndarray,
+    angle: float = 0.0,
+    k_r: float = 50.0,
+    eta_f: float = 1.0,
+    tau: float = 1.0e-7,
+    freq: float = 1.0e3,
+    file_out_str: str = "opt_params_min.pkl",
+    display_results: bool = False,
+    well_name: str = "Unknown well",
+    **opt_kwargs,
+):
+    """T-Matrix optimisation adapted to a case with detailed information available, such as in a development or production
+    setting. Mineral and fluid composition should be known on a sample basis. Inclusion parameters are regarded as
+    unknown and they are optimised for.
+
+    Parameters
+    ----------
+    k_min :
+        Effective mineral bulk modulus [Pa].
+    mu_min :
+        Effective mineral shear modulus [Pa].
+    rho_min :
+        Effective mineral bulk density [kg/m^3].
+    k_fl :
+        Effective fluid bulk modulud [Pa].
+    rho_fl :
+        Effective fluid density [kg/m^3].
+    por :
+        Inclusion porosity [ratio].
+    vp :
+        Compressional velocity log [m/s].
+    vs :
+        Shear velocity log [m/s].
+    rhob :
+        Bulk density log [kg/m^3].
+    angle : float
+        Angle of symmetry plane [degrees]
+    k_r :
+        Permeability [mD].
+    eta_f :
+        Fluid viscosity [cP].
+    tau :
+        Relaxation time constant [s].
+    freq :
+        Signal frequency [Hz].
+    file_out_str :
+        Output file name (string) to store optimal parameters (pickle format).
+    display_results :
+        Display optimal parameters in a window after run.
+    well_name:
+        Name of well to be displayed in info box title.
+    opt_kwargs:
+        Additional keywords to be passed to optimisation function
+
+    Returns
+    -------
+    tuple
+        vp_mod, vs_mod, rho_mod, ai_mod, vpvs_mod - modelled logs, vp_res, vs_res, rho_res - residual logs.
+    """
+    # 1. Preparation that is independent of search for minimum possible aspect ratio for second inclusion set
+
+    # Optimisation function for selected parameters, with effective mineral properties known and 2 inclusion sets
+    opt_fun = curve_fit_2_inclusion_sets
+    rhob_mod = rho_min * (1 - por) + rho_fl * por
+    rhob_res = rhob - rhob_mod
+    # PETEC adapted inputs: include fluid data and other params in x_data
+    por, angle, k_r, eta_f, tau, freq, def_vp_vs_ratio = gen_utilities.dim_check_vector(
+        (por, angle, k_r, eta_f, tau, freq, DEF_VP_VS_RATIO)
+    )
+    x_data = np.stack(
+        (
+            por,
+            k_min,
+            mu_min,
+            rho_min,
+            k_fl,
+            rho_fl,
+            angle,
+            k_r,
+            eta_f,
+            tau,
+            freq,
+            def_vp_vs_ratio,
+        ),
+        axis=1,
+    )
+    # Set weight to vs to give vp and vs similar influence on optimisation
+    y_data = np.stack([vp, vs * DEF_VP_VS_RATIO], axis=1)
+
+    # 2. Search for minimum aspect ratio of inclusion set no. 2, given that inclusion set no. 1 will at least represent
+    # 50% of inclusions. Minimum aspect ratio is linked to the porosity, and the most conservative estimate is to use
+    # the maximum value. This is likely to deteriorate the optimisation results, so a search for a more appropriate
+    # value that still produces valid results is made
+    percentiles = [50, 75, 80, 85, 90, 95, 99, 100]
+    valid_result = False
+    vel_mod = None
+    vel_res = None
+    opt_params = None
+    while not valid_result and percentiles:
+        try:
+            # Make sure that parameters are not in conflict with T Matrix assumptions
+            min_v1 = 0.5
+            max_por = np.percentile(por, percentiles[0])
+            min_a2 = (1.0 - min_v1) * max_por
+            # Test with all parameters in the range 0.0 - 1.0 for best optimiser performance
+            # Params:               f_ani f_con a1   a2    v1
+            lower_bound = np.array([0.0, 0.0, 0.5, min_a2, min_v1], dtype=float)
+            upper_bound = np.array([1.0, 1.0, 1.0, 0.30, 1.0], dtype=float)
+            x0 = (upper_bound + lower_bound) / 2.0
+            # Optimisation step without fluid substitution
+            vel_mod, vel_res, opt_params = gen_opt_routine(
+                opt_fun, x_data, y_data, x0, lower_bound, upper_bound, **opt_kwargs
+            )
+            valid_result = True
+        except ValueError:
+            percentiles.pop(0)
+            valid_result = False
+
+    if not valid_result:
+        raise ValueError(
+            f"{__file__}: unable to find stable value for T Matrix optimisation, second inclusion"
+        )
+
+    # Reshape outputs and remove weight from vs
+    vp_mod, vs_mod = [arr.flatten() for arr in np.split(vel_mod, 2, axis=1)]
+    ai_mod = vp_mod * rhob_mod
+    vp_res, vs_res = [arr.flatten() for arr in np.split(vel_res, 2, axis=1)]
+    vs_mod = vs_mod / DEF_VP_VS_RATIO
+    vs_res = vs_res / DEF_VP_VS_RATIO
+    vpvs_mod = vp_mod / vs_mod
+    # Save the optimal parameters
+    save_opt_params("min", opt_params, file_out_str, well_name=well_name)
+    if display_results:
+        from .opt_subst_utilities import opt_param_to_ascii
+
+        opt_param_to_ascii(file_out_str, well_name=well_name)
+
+    return vp_mod, vs_mod, rhob_mod, ai_mod, vpvs_mod, vp_res, vs_res, rhob_res

rock_physics_open/t_matrix_models/t_matrix_vector/__init__.py

@@ -0,0 +1,12 @@
+from .array_functions import array_inverse, array_matrix_mult
+from .calc_pressure import calc_pressure_vec
+from .pressure_input import pressure_input_utility
+from .t_matrix_vec import t_matrix_porosity_vectorised
+
+__all__ = [
+    "array_inverse",
+    "array_matrix_mult",
+    "calc_pressure_vec",
+    "pressure_input_utility",
+    "t_matrix_porosity_vectorised",
+]