pyrestoolbox 2.1.4__tar.gz → 2.2.1__tar.gz

{pyrestoolbox-2.1.4/pyrestoolbox.egg-info → pyrestoolbox-2.2.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyrestoolbox
-Version: 2.1.4
+Version: 2.2.1
 Summary: pyResToolbox - A collection of Reservoir Engineering Utilities
 Home-page: https://github.com/mwburgoyne/pyResToolbox
 Author: Mark W. Burgoyne
@@ -12,7 +12,6 @@ Classifier: License :: OSI Approved :: GNU General Public License v3 or later (G
 Classifier: Operating System :: OS Independent
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: requests
 Requires-Dist: numpy
 Requires-Dist: scipy
 Requires-Dist: pandas
@@ -20,7 +19,6 @@ Requires-Dist: tabulate
 Requires-Dist: gwr_inversion
 Requires-Dist: mpmath
 Requires-Dist: openpyxl
-Requires-Dist: setuptools
 Dynamic: author
 Dynamic: author-email
 Dynamic: description

{pyrestoolbox-2.1.4 → pyrestoolbox-2.2.1}/pyrestoolbox/brine/brine.py

@@ -21,25 +21,39 @@
           Contact author at mark.w.burgoyne@gmail.com
 """
 
-import sys
-from collections import Counter
-import glob
-from enum import Enum
-import pkg_resources
-
 import numpy as np
 import numpy.typing as npt
-
 import pandas as pd
+
+from typing import Tuple
 from tabulate import tabulate
-from typing import Union, List, Tuple
 
 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
 
+def _Eq41(t, input_array):
+    """Eq 4.1 from McCain Petroleum Reservoir Fluid Properties"""
+    t2 = t / 100
+    return (
+        input_array[1] * t2 ** 2 + input_array[2] * t2 + input_array[3]
+    ) / (input_array[4] * t2 ** 2 + input_array[5] * t2 + 1)
+
+# Spivey coefficient tables (shared by brine_props and brine_props_co2)
+_RHOW_T70_ARR = [0, -0.127213, 0.645486, 1.03265, -0.070291, 0.639589]
+_EWT_ARR = [0, 4.221, -3.478, 6.221, 0.5182, -0.4405]
+_FWT_ARR = [0, -11.403, 29.932, 27.952, 0.20684, 0.3768]
+_DM2T_ARR = [0, -0.00011149, 0.000175105, -0.00043766, 0, 0]
+_DM32T_ARR = [0, -0.0008878, -0.0001388, -0.00296318, 0, 0.51103]
+_DM1T_ARR = [0, 0.0021466, 0.012427, 0.042648, -0.081009, 0.525417]
+_DM12T_ARR = [0, 0.0002356, -0.0003636, -0.0002278, 0, 0]
+_EMT_ARR = [0, 0, 0, 0.1249, 0, 0]
+_FM32T_ARR = [0, -0.617, -0.747, -0.4339, 0, 10.26]
+_FM1T_ARR = [0, 0, 9.917, 5.1128, 0, 3.892]
+_FM12T_ARR = [0, 0.0365, -0.0369, 0, 0, 0]
+
 def brine_props(p: float, degf: float, wt: float=0, ch4_sat: float=0) -> Tuple:
     """ Calculates Brine properties from modified Spivey Correlation per McCain Petroleum Reservoir Fluid Properties pg 160
         Returns Tuple of (Bw (rb/stb), Density (sg), viscosity (cP), Compressibility (1/psi), Rw GOR (scf/stb))
@@ -48,12 +62,7 @@ def brine_props(p: float, degf: float, wt: float=0, ch4_sat: float=0) -> Tuple:
         wt: Salt wt% (0-100)
         ch4_sat: Degree of methane saturation (0 - 1)
     """
-
-    def Eq41(t, input_array):  # From McCain Petroleum Reservoir Fluid Properties
-        t2 = t / 100
-        return (
-            input_array[1] * t2 ** 2 + input_array[2] * t2 + input_array[3]
-        ) / (input_array[4] * t2 ** 2 + input_array[5] * t2 + 1)
+    Eq41 = _Eq41
 
     Mpa = p * 0.00689476  # Pressure in mPa
     degc = (degf - 32) / 1.8  # Temperature in deg C
@@ -62,17 +71,17 @@ def brine_props(p: float, degf: float, wt: float=0, ch4_sat: float=0) -> Tuple:
         1000 * (wt / 100) / (58.4428 * (1 - (wt / 100)))
     )  # Molar concentration of NaCl from wt % in gram mol/kg water
 
-    rhow_t70_arr = [0, -0.127213, 0.645486, 1.03265, -0.070291, 0.639589]
-    Ewt_arr = [0, 4.221, -3.478, 6.221, 0.5182, -0.4405]
-    Fwt_arr = [0, -11.403, 29.932, 27.952, 0.20684, 0.3768]
-    Dm2t_arr = [0, -0.00011149, 0.000175105, -0.00043766, 0, 0]
-    Dm32t_arr = [0, -0.0008878, -0.0001388, -0.00296318, 0, 0.51103]
-    Dm1t_arr = [0, 0.0021466, 0.012427, 0.042648, -0.081009, 0.525417]
-    Dm12t_arr = [0, 0.0002356, -0.0003636, -0.0002278, 0, 0]
-    Emt_arr = [0, 0, 0, 0.1249, 0, 0]
-    Fm32t_arr = [0, -0.617, -0.747, -0.4339, 0, 10.26]
-    Fm1t_arr = [0, 0, 9.917, 5.1128, 0, 3.892]
-    Fm12t_arr = [0, 0.0365, -0.0369, 0, 0, 0]
+    rhow_t70_arr = _RHOW_T70_ARR
+    Ewt_arr = _EWT_ARR
+    Fwt_arr = _FWT_ARR
+    Dm2t_arr = _DM2T_ARR
+    Dm32t_arr = _DM32T_ARR
+    Dm1t_arr = _DM1T_ARR
+    Dm12t_arr = _DM12T_ARR
+    Emt_arr = _EMT_ARR
+    Fm32t_arr = _FM32T_ARR
+    Fm1t_arr = _FM1T_ARR
+    Fm12t_arr = _FM12T_ARR
 
     rhow_t70 = Eq41(degc, rhow_t70_arr)
     Ewt = Eq41(degc, Ewt_arr)
@@ -171,7 +180,7 @@ def brine_props(p: float, degf: float, wt: float=0, ch4_sat: float=0) -> Tuple:
 
     try:
         mch4w = np.exp(A_t * np.power(np.log(Mpa - vap_pressure), 2) + B_t * np.log(Mpa - vap_pressure) + C_t)  # Eq 4.15
-    except:
+    except (ValueError, FloatingPointError):
         mch4w = 0
     
     u_arr = [
@@ -352,9 +361,9 @@ BBL2CUFT = 5.614583333 # cuft in a bbl
 
 #============================================================================
 #    ***  Mutual solubilities between CO2 and Brine - Calculated with  ***
-#  A Phase-Partitioning Model for CO2�Brine Mixtures at Elevated Temperatures 
+#  A Phase-Partitioning Model for CO2�Brine Mixtures at Elevated Temperatures 
 #  and Pressures: Application to CO2-Enhanced Geothermal Systems
-#  Nicolas Spycher & Karsten Pruess, Transp Porous Med (2010) 82:173�196
+#  Nicolas Spycher & Karsten Pruess, Transp Porous Med (2010) 82:173�196
 #  DOI 10.1007/s11242-009-9425-y
 #============================================================================
 
@@ -712,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):
@@ -852,9 +870,9 @@ class CO2_Brine_Mixture():
         self.Bprime = B
     
     #============================================================================
-    #  A Phase-Partitioning Model for CO2�Brine Mixtures at Elevated Temperatures 
+    #  A Phase-Partitioning Model for CO2�Brine Mixtures at Elevated Temperatures 
     #  and Pressures: Application to CO2-Enhanced Geothermal Systems
-    #  Nicolas Spycher & Karsten Pruess, Transp Porous Med (2010) 82:173�196
+    #  Nicolas Spycher & Karsten Pruess, Transp Porous Med (2010) 82:173�196
     #  DOI 10.1007/s11242-009-9425-y
     #============================================================================
     
@@ -1069,26 +1087,18 @@ class CO2_Brine_Mixture():
         Mpa = pBar * 0.1                    # Pressure in mPa
         tKel = degc + CEL2KEL               # Temperature in deg K
         
-        def Eq41(t, input_array):  # From McCain Petroleum Reservoir Fluid Properties
-            t2 = t / 100
-            return (input_array[1] * t2 ** 2 + input_array[2] * t2 + input_array[3]) / (input_array[4] * t2 ** 2 + input_array[5] * t2 + 1)
-    
-        # Table 4-6 Coefficients
-        rhow_t70_arr = [0, -0.127213, 0.645486, 1.03265, -0.070291, 0.639589]
-        Ewt_arr = [0, 4.221, -3.478, 6.221, 0.5182, -0.4405]
-        Fwt_arr = [0, -11.403, 29.932, 27.952, 0.20684, 0.3768]
-        
-        # Table 4-7 Coefficients
-        Dm2t_arr = [0, -0.00011149, 0.000175105, -0.00043766, 0, 0]
-        Dm32t_arr = [0, -0.0008878, -0.0001388, -0.00296318, 0, 0.51103]
-        Dm1t_arr = [0, 0.0021466, 0.012427, 0.042648, -0.081009, 0.525417]
-        Dm12t_arr = [0, 0.0002356, -0.0003636, -0.0002278, 0, 0]
-        
-        # Table 4-8 Coefficients
-        Emt_arr = [0, 0, 0, 0.1249, 0, 0]
-        Fm32t_arr = [0, -0.617, -0.747, -0.4339, 0, 10.26]
-        Fm1t_arr = [0, 0, 9.917, 5.1128, 0, 3.892]
-        Fm12t_arr = [0, 0.0365, -0.0369, 0, 0, 0]
+        Eq41 = _Eq41
+        rhow_t70_arr = _RHOW_T70_ARR
+        Ewt_arr = _EWT_ARR
+        Fwt_arr = _FWT_ARR
+        Dm2t_arr = _DM2T_ARR
+        Dm32t_arr = _DM32T_ARR
+        Dm1t_arr = _DM1T_ARR
+        Dm12t_arr = _DM12T_ARR
+        Emt_arr = _EMT_ARR
+        Fm32t_arr = _FM32T_ARR
+        Fm1t_arr = _FM1T_ARR
+        Fm12t_arr = _FM12T_ARR
         
         # Table 4-14 Mao-Duan Coefficients
         d = [0, 2885310, -11072.577, -9.0834095, 0.030925651, -0.0000274071, -1928385.1, 5621.6046, 13.82725, -0.047609523, 0.000035545041]
@@ -1115,7 +1125,7 @@ class CO2_Brine_Mixture():
         
         # -- CO2-Free Brine Density (gm/cm3)
         def brine_denw(Mpa):
-            # cw(T, p), in MPa�1, of pure water at temperature T and pressure p,
+            # cw(T, p), in MPa�1, of pure water at temperature T and pressure p,
             cwtp = (1 / 70) * (1 / (Ewt * (Mpa / 70) + Fwt))        # Eq 4.2
     
             #Density of pure water at temperature T and pressure p.
@@ -1168,7 +1178,7 @@ class CO2_Brine_Mixture():
             return (1.0 + mRat * xRat) / (vPhi * xRat / MwB + 1.0 / rhoBRnoCO2) # --Equation 18 of Garcia paper
         
         # Correct CO2 free brine viscosity for dissolved CO2
-        # Using approach from "Viscosity Models and Effects of Dissolved CO2", Akand W. Islam and Eric S. Carlson (Jul 2012), Energy Fuels 2012, 26, 8, 5330�5336, https://doi.org/10.1021/ef3006228
+        # Using approach from "Viscosity Models and Effects of Dissolved CO2", Akand W. Islam and Eric S. Carlson (Jul 2012), Energy Fuels 2012, 26, 8, 5330�5336, https://doi.org/10.1021/ef3006228
         def co2_vis_brine(cP_brine, xCO2):
             # Uses CO2 free brine viscosity (cP) and mole fraction CO2 in brine (xCO2), and returns cP
             return cP_brine * (1 + 4.65 * xCO2**1.0134)
@@ -1237,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.1.4 → pyrestoolbox-2.2.1}/pyrestoolbox/classes/classes.py

@@ -28,11 +28,13 @@ class z_method(Enum):  # Gas Z-Factor calculation model
     HY = 1
     WYW = 2
     BUR = 3
+    BNS = 3
 
 class c_method(Enum):  # Gas critical properties calculation method
     PMC = 0
     SUT = 1
     BUR = 2
+    BNS = 3
 
 class pb_method(Enum):  # Bubble point calculation method
     STAN = 0

{pyrestoolbox-2.1.4 → pyrestoolbox-2.2.1}/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.1.4 → pyrestoolbox-2.2.1}/pyrestoolbox/docs/changelist.rst

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