pyrestoolbox 2.2__tar.gz → 2.2.2__tar.gz

{pyrestoolbox-2.2/pyrestoolbox.egg-info → pyrestoolbox-2.2.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyrestoolbox
-Version: 2.2
+Version: 2.2.2
 Summary: pyResToolbox - A collection of Reservoir Engineering Utilities
 Home-page: https://github.com/mwburgoyne/pyResToolbox
 Author: Mark W. Burgoyne

{pyrestoolbox-2.2 → pyrestoolbox-2.2.2}/pyrestoolbox/brine/brine.py

@@ -23,12 +23,14 @@
 
 import numpy as np
 import numpy.typing as npt
+import pandas as pd
 
 from typing import Tuple
+from tabulate import tabulate
 
 import pyrestoolbox.gas as gas # Needed for Z-Factor
 from pyrestoolbox.classes import z_method, c_method, pb_method, rs_method, bo_method, uo_method, deno_method, co_method, kr_family, kr_table, class_dic
-from pyrestoolbox.shared_fns import convert_to_numpy, process_output
+from pyrestoolbox.shared_fns import convert_to_numpy, process_output, halley_solve_cubic
 from pyrestoolbox.validate import validate_methods
 from pyrestoolbox.constants import R, psc, tsc, degF2R, tscr, scf_per_mol, CUFTperBBL, WDEN, MW_CO2, MW_H2S, MW_N2, MW_AIR, MW_H2
 
@@ -719,31 +721,40 @@ class CO2_Brine_Mixture():
     #  Cubic Polynomial Solver: f(Z) = Z**3 + E2*Z**2 + E1*Z + E1 = 0
     #=======================================================================
         self.repeat = False
-        Z = np.roots(np.array([1.0, e2, e1, e0]))
-        Z = np.array([x for x in Z if np.isreal(x)]) # Keep only real results
-        if len(Z) > 1: # Evaluate which root to use per Eqs 25 and 26 in Spycher & Pruess (2003)
-            vgas, vliq = max(Z), min(Z)
-            
+
+        # Try Halley for all roots
+        roots = halley_solve_cubic(e2, e1, e0, flag=0)
+
+        # Fallback to np.roots if Halley returned None
+        if roots is None:
+            Z = np.roots(np.array([1.0, e2, e1, e0]))
+            Z = np.array([x for x in Z if np.isreal(x)])  # Keep only real results
+            roots = np.real(Z)
+
+        if len(roots) > 1: # Evaluate which root to use per Eqs 25 and 26 in Spycher & Pruess (2003)
+            vgas, vliq = max(roots), min(roots)
+
             w1 = self.pBar*(vgas - vliq)
             w2 = RGASCON * self.tKel * np.log((vgas - self.bMix)/(vliq - self.bMix)) + self.aMix/(self.tKel**0.5 * self.bMix) * np.log((vgas + self.bMix) * vliq / ((vliq + self.bMix) * vgas))
-            
+
             if w2 - w1 > 0:
-                Z[0] = max(Z)
+                result = max(roots)
                 if self.CO2_sat:          # CO2 was saturated in previous iteration, but now its not
                     self.CO2_sat = False
                     self.repeat = True
             else:
-                Z[0] = min(Z)
+                result = min(roots)
                 if not self.CO2_sat:
                     self.CO2_sat = True
                     self.repeat = True
         else:
+            result = roots[0]
             if self.CO2_sat:          # CO2 was saturated in previous iteration, but now its not
                 self.CO2_sat = False
                 self.repeat = True
-                
-            
-        return np.real(Z[0])
+
+
+        return np.real(result)
 
         
     def MolarVolume(self):
@@ -1236,4 +1247,116 @@ class CO2_Brine_Mixture():
         # Undersaturated compressibility = 1/V dV/dP
         c_usat = 1 - sg_CO2_Brine / sg_CO2_Brine_ # 1/Bar    
         
-        return ([sg_CO2_Brine, sg_brine, rhowtp], [cP_CO2_brine, cP_brine, cP_freshwater], viscosblty, [bw, brine_res_vol, bw_freshwater], rs, c_usat)
+        return ([sg_CO2_Brine, sg_brine, rhowtp], [cP_CO2_brine, cP_brine, cP_freshwater], viscosblty, [bw, brine_res_vol, bw_freshwater], rs, c_usat)
+
+
+def make_pvtw_table(
+    pi: float,
+    degf: float,
+    wt: float = 0,
+    ch4_sat: float = 0,
+    pmin: float = 500,
+    pmax: float = 10000,
+    nrows: int = 20,
+    export: bool = False,
+) -> dict:
+    """ Generates a PVTW (water PVT) table over a pressure range using brine_props (Spivey correlation).
+        Follows the pattern of make_bot_og from the oil module.
+
+        pi: Initial (reference) pressure (psia)
+        degf: Temperature (deg F)
+        wt: Salt wt% (0-100), default 0
+        ch4_sat: Degree of methane saturation (0 - 1), default 0
+        pmin: Minimum pressure for table (psia), default 500
+        pmax: Maximum pressure for table (psia), default 10000
+        nrows: Number of rows in table, default 20
+        export: If True, writes PVTW.INC (ECLIPSE keyword) and pvtw_table.xlsx
+
+        Returns dict with keys:
+            table: pandas DataFrame with columns Pressure, Bw, Density, Viscosity, Cw, Rsw
+            pref: Reference pressure (psia)
+            bw_ref: Bw at reference pressure (rb/stb)
+            cw_ref: Compressibility at reference pressure (1/psi)
+            visw_ref: Viscosity at reference pressure (cP)
+            rsw_ref: Rsw at reference pressure (scf/stb)
+            den_ref: Density (sg) at reference pressure
+    """
+    # Build pressure grid, ensuring pi is included
+    pressures = list(np.linspace(pmin, pmax, nrows))
+    if pi not in pressures:
+        pressures.append(pi)
+    pressures = sorted(set(pressures))
+
+    bws, dens, visws, cws, rsws = [], [], [], [], []
+    for p in pressures:
+        bw, lden, visw, cw, rsw = brine_props(p=p, degf=degf, wt=wt, ch4_sat=ch4_sat)
+        bws.append(bw)
+        dens.append(lden)
+        visws.append(visw)
+        cws.append(cw)
+        rsws.append(rsw)
+
+    df = pd.DataFrame()
+    df["Pressure (psia)"] = pressures
+    df["Bw (rb/stb)"] = bws
+    df["Density (sg)"] = dens
+    df["Viscosity (cP)"] = visws
+    df["Cw (1/psi)"] = cws
+    df["Rsw (scf/stb)"] = rsws
+
+    # Reference properties at pi
+    bw_ref, den_ref, visw_ref, cw_ref, rsw_ref = brine_props(
+        p=pi, degf=degf, wt=wt, ch4_sat=ch4_sat
+    )
+
+    if export:
+        # Write ECLIPSE PVTW keyword
+        # PVTW format: Pref  Bw  Cw  Visw  Viscosibility
+        # Viscosibility set to 0 (constant viscosity assumption)
+        pvtw_line = f"  {pi:.1f}  {bw_ref:.6f}  {cw_ref:.6e}  {visw_ref:.4f}  0.0"
+        fileout = f"-- Generated by pyResToolbox make_pvtw_table\n"
+        fileout += f"-- Temperature: {degf:.1f} deg F, Salt: {wt:.1f} wt%, CH4 sat: {ch4_sat:.2f}\n"
+        fileout += f"PVTW\n{pvtw_line} /\n"
+        with open("PVTW.INC", "w") as f:
+            f.write(fileout)
+
+        # Write full table to Excel
+        df.to_excel("pvtw_table.xlsx", index=False)
+
+    return {
+        "table": df,
+        "pref": pi,
+        "bw_ref": bw_ref,
+        "cw_ref": cw_ref,
+        "visw_ref": visw_ref,
+        "rsw_ref": rsw_ref,
+        "den_ref": den_ref,
+    }
+
+
+class SoreideWhitson:
+    """ Soreide-Whitson (1992) model for gas solubility in water/brine.
+
+        Planned support for multicomponent gas mixtures containing:
+        C1, C2, C3, nC4, CO2, H2S, N2, H2
+
+        Will calculate:
+        - Mole fraction of dissolved gas components in aqueous phase
+        - Mole fraction of vaporised water in gas phase
+        - Water content of gas (stb/mmscf)
+        - Gas solubility in water (scf/stb)
+
+        Supports fresh and saline water (NaCl equivalent).
+
+        Reference:
+            Soreide, I. and Whitson, C.H., "Peng-Robinson Predictions for Hydrocarbons,
+            CO2, N2, and H2S with Pure Water and NaCl Brine", Fluid Phase Equilibria,
+            77, 217-240, 1992.
+    """
+
+    def __init__(self, **kwargs):
+        raise NotImplementedError(
+            "SoreideWhitson model is not yet implemented. "
+            "Planned support includes multicomponent gas (C1, C2, C3, nC4, CO2, H2S, N2, H2) "
+            "solubility in fresh/saline water using the Soreide-Whitson (1992) PR-EOS approach."
+        )

{pyrestoolbox-2.2 → pyrestoolbox-2.2.2}/pyrestoolbox/docs/brine.rst

@@ -161,5 +161,106 @@ Usage example for 175 Bara x 85 degC and 0% NaCl brine:
 
     >>> mix = brine.CO2_Brine_Mixture(pres = 175, temp = 85)
     >>> mix.Rs  # Returns sm3 dissolved CO2 / sm3 Brine
-    24.742923469934272   
+    24.742923469934272
+
+pyrestoolbox.brine.make_pvtw_table
+======================
+
+.. code-block:: python
+
+    make_pvtw_table(pi, degf, wt=0, ch4_sat=0, pmin=500, pmax=10000, nrows=20, export=False) -> dict
+
+Generates a PVTW (water PVT) table over a pressure range using the Spivey correlation (brine_props).
+Optionally exports ECLIPSE PVTW keyword file and Excel spreadsheet.
+
+.. list-table:: Inputs
+   :widths: 10 15 40
+   :header-rows: 1
+
+   * - Parameter
+     - Type
+     - Description
+   * - pi
+     - float
+     - Initial (reference) pressure (psia)
+   * - degf
+     - float
+     - Temperature (deg F)
+   * - wt
+     - float
+     - Salt weight% in the brine (0 - 100). Default 0
+   * - ch4_sat
+     - float
+     - Degree of methane saturation (0 - 1). Default 0
+   * - pmin
+     - float
+     - Minimum table pressure (psia). Default 500
+   * - pmax
+     - float
+     - Maximum table pressure (psia). Default 10000
+   * - nrows
+     - int
+     - Number of table rows. Default 20
+   * - export
+     - bool
+     - If True, writes PVTW.INC and pvtw_table.xlsx. Default False
+
+.. list-table:: Return dict keys
+   :widths: 10 15 40
+   :header-rows: 1
+
+   * - Key
+     - Type
+     - Description
+   * - table
+     - DataFrame
+     - Pressure, Bw, Density, Viscosity, Cw, Rsw
+   * - pref
+     - float
+     - Reference pressure (psia)
+   * - bw_ref
+     - float
+     - Bw at reference pressure (rb/stb)
+   * - cw_ref
+     - float
+     - Compressibility at reference pressure (1/psi)
+   * - visw_ref
+     - float
+     - Viscosity at reference pressure (cP)
+   * - rsw_ref
+     - float
+     - Rsw at reference pressure (scf/stb)
+   * - den_ref
+     - float
+     - Density (sg) at reference pressure
+
+Examples:
+
+.. code-block:: python
+
+    >>> from pyrestoolbox import brine
+    >>> result = brine.make_pvtw_table(pi=3000, degf=200, wt=0, ch4_sat=0)
+    >>> print(result['bw_ref'])
+    >>> print(result['table'].head())
+
+pyrestoolbox.brine.SoreideWhitson
+======================
+
+.. code-block:: python
+
+    SoreideWhitson(**kwargs) -> raises NotImplementedError
+
+Placeholder for the Soreide-Whitson (1992) PR-EOS model for gas solubility in water/brine.
+
+Planned support includes multicomponent gas mixtures (C1, C2, C3, nC4, CO2, H2S, N2, H2)
+solubility in fresh and saline water.
+
+Planned return values:
+- Mole fraction of dissolved gas components in aqueous phase
+- Mole fraction of vaporised water in gas phase
+- Water content of gas (stb/mmscf)
+- Gas solubility in water (scf/stb)
+
+Reference: Soreide, I. and Whitson, C.H., "Peng-Robinson Predictions for Hydrocarbons,
+CO2, N2, and H2S with Pure Water and NaCl Brine", Fluid Phase Equilibria, 77, 217-240, 1992.
 

{pyrestoolbox-2.2 → pyrestoolbox-2.2.2}/pyrestoolbox/docs/changelist.rst

@@ -1,3 +1,8 @@
+Changelist in 2.2:
+
+- Bugfixes.
+
+
 Changelist in 2.1.3:
 
 - Updated viscosity parameters for BUR method.

{pyrestoolbox-2.2 → pyrestoolbox-2.2.2}/pyrestoolbox/docs/gas.rst

@@ -210,22 +210,7 @@ Examples:
     
     >>> gas.gas_z(p=[1000, 2000], sg=0.75, degf=160, cmethod='SUT', n2 = 0.02, co2 = 0.17)
     array([0.91920553, 0.87196032])
- 
- def gas_ug(
-    p: npt.ArrayLike,
-    sg: float,
-    degf: float,
-    zmethod: z_method = z_method.DAK,
-    cmethod: c_method = c_method.PMC,
-    co2: float = 0,
-    h2s: float = 0,
-    n2: float = 0,
-    h2: float = 0,
-    tc: float = 0,
-    pc: float = 0,
-    zee: float = 0,
-    ugz = False
-    
+
 pyrestoolbox.gas.gas_ug
 ===================
 
@@ -741,7 +726,7 @@ pyrestoolbox.gas.gas_rate_radial
     gas_rate_radial(k, h, pr, pwf, r_w, r_ext, degf, zmethod='DAK, cmethod='PMC', S = 0, D = 0, sg = 0.75, n2 = 0, co2 = 0, h2s = 0, tc  = 0, pc = 0) -> float or np.array
 
 Returns gas rate (mscf/day) for radial flow using Darcy pseudo steady state equation & gas pseudopressure. 
-Arrays can be used for any one of k, h, pr or pwf, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
+Arrays can be used for any one of k, h, pr or pwf, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
 
 .. list-table:: Inputs
    :widths: 10 15 40
@@ -806,7 +791,6 @@ Examples:
 
 .. code-block:: python
 
-    >>> from pyrestoolbox import pyrestoolbox as rtb
     >>> gas.gas_rate_radial(k=5, h=50, pr=2000, pwf=750, r_w=0.3, r_ext=1500, degf=180, sg = 0.75, D = 0.01, S=5)
     10269.669190157822
     
@@ -822,7 +806,7 @@ pyrestoolbox.gas.gas_rate_linear
     gas_rate_linear(k, pr, pwf, area, length, degf, zmethod='DAK, cmethod='PMC', sg = 0.75, n2 = 0, co2 = 0, h2s = 0, tc  = 0, pc = 0) -> float or np.array
 
 Returns gas rate (mscf/day) for linear flow using Darcy steady state equation & gas pseudopressure. 
-Arrays can be used for any one of k, pr, pwf or area, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
+Arrays can be used for any one of k, pr, pwf or area, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
 
 
 .. list-table:: Inputs

{pyrestoolbox-2.2 → pyrestoolbox-2.2.2}/pyrestoolbox/docs/library.rst

@@ -65,7 +65,7 @@ Examples:
 
 .. code-block:: python
 
-    >>> library.prop(comp = 'C3', prop = 'VTran'), rtb.comp_library.prop(comp = 'C3', prop = 'VTran', model='SRK')
+    >>> library.prop(comp = 'C3', prop = 'VTran'), library.prop(comp = 'C3', prop = 'VTran', model='SRK')
     (-0.06381, 0.09075)
 
 
@@ -158,5 +158,5 @@ Example:
 
     >>> library.df
 
-.. image:: https://github.com/mwburgoyne/pyResToolbox/blob/main/docs/img/properties_df.png
+.. image:: img/properties_df.png
     :alt: DataFrame of Component Library data

{pyrestoolbox-2.2 → pyrestoolbox-2.2.2}/pyrestoolbox/docs/oil.rst

@@ -57,7 +57,7 @@ pyResToolBox uses class objects to track calculation options through the functio
      - Method for calculating oil formation volume factor. Defaults to 'MCAIN'
        Options are:
         + 'STAN': Standing Correlation
-        + 'MCAIN': McCain approach, calculating from densities � Default
+        + 'MCAIN': McCain approach, calculating from densities � Default
 
 Users can specify which calculation method to use either by passing an option string, or a class object to any given function. The implementation of class objects should make it easier to program in an IDE that supports type hinting
 
@@ -148,7 +148,7 @@ Examples:
 
 .. code-block:: python
 
-    >>> oil.ja_sg(mw=150, ja=0.5)
+    >>> oil.oil_ja_sg(mw=150, ja=0.5)
     0.8583666666666667
 
 pyrestoolbox.oil.oil_twu_props
@@ -239,7 +239,7 @@ pyrestoolbox.oil.oil_pbub
 
 .. code-block:: python
 
-    pbub(api, degf, rsb, sg_g =0, sg_sp =0, pbmethod ='VALMC') -> float
+    oil_pbub(api, degf, rsb, sg_g =0, sg_sp =0, pbmethod ='VALMC') -> float
 
 Returns bubble point pressure (psia) calculated with different correlations. 
 At least one of sg_g and sg_sp must be supplied. This function will make simple assumption to estimate missing gas sg if only one is provided.
@@ -344,15 +344,15 @@ Returns solution gas oil ratio (scf/stb) calculated with different correlations.
    * - degf
      - float
      - Oil Temperature (deg F)
-   * - sg_g
+   * - sg_sp
      - float
-     - Weighted average specific gravity of surface gas, inclusive of gas evolved after separation (relative to air).   
+     - Separator gas gravity (relative to air).
    * - p
      - float
-     - Pressure (psia). 
+     - Pressure (psia).
    * - pb
      - float
-     - Original bubble point pressure (psia) 
+     - Original bubble point pressure (psia)
    * - rsb
      - float
      - Original solution GOR at original bubble point pressure (scf/stb)
@@ -367,7 +367,7 @@ Examples:
 
 .. code-block:: python
 
-    >>> oil.oil_rs(api = 43, degf = 185, sg_sp =0 .72, p = 3000, pb = 5179.5, rsb = 2370)
+    >>> oil.oil_rs(api = 43, degf = 185, sg_sp=0.72, p = 3000, pb = 5179.5, rsb = 2370)
     1017.9424240354475
     
     >>> oil.oil_rs(api=43, degf=185, sg_sp=0.72, p=3000, rsb =2370)
@@ -610,7 +610,7 @@ pyrestoolbox.oil.make_bot_og
 
 .. code-block:: python
 
-    make_bot_og(pi, api, degf, sg_g, pmax, pb =0, rsb =0, pmin =14.7, nrows = 20, wt =0, ch4_sat =0, comethod='EXPLT', zmethod='DAK', rsmethod='VELAR', cmethod'PMC', denomethod='SWMH', bomethod='MCAIN', pbmethod='VALMC', export=False) -> tuple
+    make_bot_og(pi, api, degf, sg_g, pmax, pb =0, rsb =0, pmin =14.7, nrows = 20, wt =0, ch4_sat =0, comethod='EXPLT', zmethod='DAK', rsmethod='VELAR', cmethod='PMC', denomethod='SWMH', bomethod='MCAIN', pbmethod='VALMC', export=False) -> tuple
 
 Creates data required for Oil-Gas-Water black oil tables. Returns dictionary of results, with index:
  - bot: Pandas table of blackoil data (for PVTO == False), or Saturated properties to pmax (if PVTO == True)
@@ -709,7 +709,7 @@ Examples:
     >>> results = oil.make_bot_og(pvto=False, pi=4000, api=38, degf=175, sg_g=0.68, pmax=5500, pb=4500, nrows=10, export=True)
     >>> df, st_deno, st_deng, res_denw, res_cw, visw, pb, rsb, rsb_frac, usat = results['bot'], results['deno'], results['deng'], results['denw'], results['cw'], results['uw'], results['pb'], results['rsb'], results['rsb_scale'], results['usat']
     >>> df
-.. image:: https://github.com/mwburgoyne/pyResToolbox/blob/main/docs/img/bot_img.png
+.. image:: img/bot_img.png
     :alt: Black Oil Table DataFrame
 
 pyrestoolbox.oil.sg_evolved_gas
@@ -749,9 +749,8 @@ Examples:
 
 .. code-block:: python
 
-    >>> oil.sg_st_gas(114.7, rsp=1500, api=42, sg_sp=0.72, degf_sp=80)
-    1.1923932340625523 
-    
+    >>> oil.sg_evolved_gas(p=2000, degf=185, rsb=2370, api=43, sg_sp=0.72)
+    0.7872810977386344
 
 pyrestoolbox.oil.sg_st_gas
 =======================
@@ -891,7 +890,7 @@ pyrestoolbox.oil.oil_rate_radial
     oil_rate_radial(k, h, pr, pwf, r_w, r_ext, uo, bo, S = 0, vogel = False, pb = 0) -> float or np.array
 
 Returns liquid rate (stb/day) for radial flow using Darcy pseudo steady state equation with optional Vogel correction.
-Arrays can be used for any one of k, h, pr or pwf, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
+Arrays can be used for any one of k, h, pr or pwf, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
 
 .. list-table:: Inputs
    :widths: 10 15 40
@@ -952,7 +951,7 @@ pyrestoolbox.oil.oil_rate_linear
     oil_rate_linear(k, pr, pwf, area, length, uo, bo, vogel = False, pb = 0) -> float or np.array
 
 Returns liquid rate (stb/day) for linear flow using Darcy steady state equation with optional Vogel correction.
-Arrays can be used for any one of k, pr, pwf or area, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
+Arrays can be used for any one of k, pr, pwf or area, returning corresponding 1-D array of rates. Using more than one input array � while not prohibited - will not return expected results 
 
 .. list-table:: Inputs
    :widths: 10 15 40

{pyrestoolbox-2.2 → pyrestoolbox-2.2.2}/pyrestoolbox/docs/simtools.rst

@@ -52,7 +52,8 @@ Examples:
 
 .. code-block:: python
 
-    >>> results = rtb.simtools.ix_extract_problem_cells()
+    >>> from pyrestoolbox import simtools
+    >>> results = simtools.ix_extract_problem_cells()
     >>> wells, grid_pres, grid_sat, grid_comp = results
     >>> grid_sat
     Processing TEST.PRT
@@ -64,7 +65,7 @@ Examples:
     Grid Saturation Change                 310  32,121,24                       13
     Grid Composition Change               1627  35,212,25                      544 
     
-.. image:: https://github.com/mwburgoyne/pyResToolbox/blob/main/docs/img/grid_sat_df.png
+.. image:: img/grid_sat_df.png
     :alt: DSorted ataFrame of grid blocks with saturation related convergence issues    
     
     
@@ -112,10 +113,10 @@ Examples:
 
 .. code-block:: python
 
-    >>> from pyrestoolbox import pyrestoolbox as rtb
+    >>> from pyrestoolbox import simtools
     >>> import matplotlib.pyplot as plt
     >>> ReDs = [1.5, 2, 3, 5, 10, 25, 1000]
-    >>> tds, pds = rtb.simtools.influence_tables(ReDs=ReDs, export=True)
+    >>> tds, pds = simtools.influence_tables(ReDs=ReDs, export=True)
     >>> 
     >>> for p, pd in enumerate(pds):
     >>>     plt.plot(tds, pd, label = str(ReDs[p]))
@@ -129,7 +130,7 @@ Examples:
     >>> plt.title('Constant Terminal Rate Solution')
     >>> plt.show()
 
-.. image:: https://github.com/mwburgoyne/pyResToolbox/blob/main/docs/img/influence.png
+.. image:: img/influence.png
     :alt: Constant Terminal Rate influence tables
 
 
@@ -173,10 +174,10 @@ Examples:
 
 .. code-block:: python
 
-    >>> rtb.simtools.zip_check_sim_deck(['FIELD_A.DATA', 'FIELD_B.afi'], console_summary=False)
+    >>> simtools.zip_check_sim_deck(['FIELD_A.DATA', 'FIELD_B.afi'], console_summary=False)
     ['INCLUDE/GridOpts.inc', 'INCLUDE/ZCORN_COORD.GRDECL', 'EPS.ixf']
     
-    >>> rtb.simtools.zip_check_sim_deck()
+    >>> simtools.zip_check_sim_deck()
       Index  File Name
     -------  ------------
           0  FIELD_A.DATA
@@ -247,7 +248,7 @@ Examples:
 
 .. code-block:: python
 
-    >>> rtb.simtools.rr_solver(zi =np.array([0.7, 0.15, 0.1, 0.05]), ki = np.array([50, 5, 0.5, 0.01]))
+    >>> simtools.rr_solver(zi =np.array([0.7, 0.15, 0.1, 0.05]), ki = np.array([50, 5, 0.5, 0.01]))
     (6,
     array([0.7406252 , 0.1570315 , 0.09469948, 0.00764382]),
     array([0.0148125 , 0.0314063 , 0.18939896, 0.76438224]),
@@ -342,9 +343,12 @@ pyrestoolbox.simtools.rel_perm_table
           
           
 Examples:
-    >>> from pyrestoolbox import pyrestoolbox as rtb
+
+.. code-block:: python
+
+    >>> from pyrestoolbox import simtools
     >>> import matplotlib.pyplot as plt
-    >>> df = rtb.simtools.rel_perm_table(rows=25, krtable='SGOF', krfamily='LET', kromax =1, krgmax =1, swc =0.2, sorg =0.15, Lo=2.5, Eo = 1.25, To = 1.75, Lg = 1.2, Eg = 1.5, Tg = 2.0)
+    >>> df = simtools.rel_perm_table(rows=25, krtable='SGOF', krfamily='LET', kromax =1, krgmax =1, swc =0.2, sorg =0.15, Lo=2.5, Eo = 1.25, To = 1.75, Lg = 1.2, Eg = 1.5, Tg = 2.0)
     >>> plt.plot(df['Sg'], df['Krgo'], c = 'r', label='Gas')
     >>> plt.plot(df['Sg'], df['Krog'], c = 'g', label='Oil')
     >>> plt.title('SGOF Gas Oil LET Relative Permeability Curves')
@@ -354,12 +358,12 @@ Examples:
     >>> plt.grid('both')
     >>> plt.plot()
 
-.. image:: https://github.com/mwburgoyne/pyResToolbox/blob/main/docs/img/sgof.png
+.. image:: img/sgof.png
     :alt: SGOF Relative Permeability Curves
 
 .. code-block:: python
 
-    >>> df = rtb.simtools.rel_perm_table(rows=25, krtable='SWOF', kromax =1, krwmax =0.25, swc =0.15, swcr = 0.2, sorw =0.15, no=2.5, nw=1.5)
+    >>> df = simtools.rel_perm_table(rows=25, krtable='SWOF', kromax =1, krwmax =0.25, swc =0.15, swcr = 0.2, sorw =0.15, no=2.5, nw=1.5)
     >>> plt.plot(df['Sw'], df['Krow'], c = 'g', label='Oil')
     >>> plt.plot(df['Sw'], df['Krwo'], c = 'b', label='Water')
     >>> plt.title('SWOF Water Oil Corey Relative Permeability Curves')
@@ -369,7 +373,7 @@ Examples:
     >>> plt.grid('both')
     >>> plt.plot()
     
-.. image:: https://github.com/mwburgoyne/pyResToolbox/blob/main/docs/img/swof.png
+.. image:: img/swof.png
     :alt: SWOF Relative Permeability Curves