pyrestoolbox 3.7.3__tar.gz → 3.7.4__tar.gz

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyrestoolbox
-Version: 3.7.3
+Version: 3.7.4
 Classifier: Programming Language :: Python :: 3
 Classifier: Programming Language :: Rust
 Classifier: License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)
@@ -55,6 +55,7 @@ Includes functions to perform calculations including;
 - Generation of AQUTAB include file influence functions for use in ECLIPSE
 - Creation of Corey, LET and Jerauld relative permeability tables in Eclipse format, with curve fitting support
 - Calculation of Methane, CO2 and multicomponent gas saturated brine properties (Soreide-Whitson VLE)
+- Multi-salt gas-free brine viscosity: IAPWS-2008 water, an ion-additive salt ratio over any combination of Na, K, Mg, Ca, Sr, Ba, Cl, Br, sulfate and carbonate, and a measured pressure factor
 - Method recommendation engine for selecting appropriate correlations based on fluid composition
 - Sensitivity analysis with parameter sweeps and tornado charts
 
@@ -127,7 +128,7 @@ Module List
    * - `library <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/library.rst>`_
      - Return critical parameters for typical single components
    * - `brine <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/brine.rst>`_
-     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE)
+     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE), multi-salt gas-free brine viscosity
    * - `layer <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/layer.rst>`_
      - Lorenz coefficient from Beta value, Lorenz coefficient from flow fraction, Lorenz coefficient to flow fraction, Lorenz coefficient to permeability array
    * - `simtools <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/simtools.rst>`_

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/README.rst

@@ -29,6 +29,7 @@ Includes functions to perform calculations including;
 - Generation of AQUTAB include file influence functions for use in ECLIPSE
 - Creation of Corey, LET and Jerauld relative permeability tables in Eclipse format, with curve fitting support
 - Calculation of Methane, CO2 and multicomponent gas saturated brine properties (Soreide-Whitson VLE)
+- Multi-salt gas-free brine viscosity: IAPWS-2008 water, an ion-additive salt ratio over any combination of Na, K, Mg, Ca, Sr, Ba, Cl, Br, sulfate and carbonate, and a measured pressure factor
 - Method recommendation engine for selecting appropriate correlations based on fluid composition
 - Sensitivity analysis with parameter sweeps and tornado charts
 
@@ -101,7 +102,7 @@ Module List
    * - `library <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/library.rst>`_
      - Return critical parameters for typical single components
    * - `brine <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/brine.rst>`_
-     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE)
+     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE), multi-salt gas-free brine viscosity
    * - `layer <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/layer.rst>`_
      - Lorenz coefficient from Beta value, Lorenz coefficient from flow fraction, Lorenz coefficient to flow fraction, Lorenz coefficient to permeability array
    * - `simtools <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/simtools.rst>`_

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "maturin"
 
 [project]
 name = "pyrestoolbox"
-version = "3.7.3"
+version = "3.7.4"
 description = "pyResToolbox - A collection of Reservoir Engineering Utilities"
 license = {text = "GPL-3.0-or-later"}
 authors = [{name = "Mark W. Burgoyne", email = "mark.w.burgoyne@gmail.com"}]

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/pyrestoolbox/brine/__init__.py

@@ -1,4 +1,5 @@
 from .brine import *
+from .brine import recommended_method  # which entry point to use, and why
 from . import pr_vphi
 from .pr_vphi import V2_inf as V_phi_pr
 from . import vphi_route

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/pyrestoolbox/brine/_lib_vle_engine.py

@@ -37,10 +37,14 @@ RACHFORD-RICE SOLVER:
 
 References:
 - Soreide & Whitson, Fluid Phase Equilibria 77 (1992) 217-240 + Errata
+- Burgoyne & Nielsen, Fluid Phase Equilibria (2026) 114824,
+  doi:10.1016/j.fluid.2026.114824 - the refreshed framework implemented here;
+  "this work" in comments below refers to this paper. Includes the H2
+  correlations (A=-14.59, B=2.184, C=0.365, kij_NA=0.468), the CH4 rational
+  form (A=-2.1642, B=1.7325, C=0.2105, MC-3 alpha), the embedded salinity
+  BIPs and the H2S salting-out treatment.
 - Yan et al., Fluid Phase Equilibria 298 (2011) 180-189 (improved CO2)
 - Nielsen & Lia, Fluid Phase Equilibria (2022) - Robust RR solver
-- This work (H2 correlations): A=-14.59, B=2.184, C=0.365, kij_NA=0.468
-- This work (CH4 rational): A=-2.1642, B=1.7325, C=0.2105 (MC-3 alpha)
 
 Author: Mark Burgoyne, Markus H. Nielsen
 Date: 2025-2026
@@ -366,7 +370,7 @@ _SW_GASES_WITH_EMBEDDED_SALINITY = {
 GASES_WITH_EMBEDDED_SALINITY = _SW_GASES_WITH_EMBEDDED_SALINITY
 
 
-# --- Proposed freshwater kij_AQ (Paper 2, MC-3 alpha) ---
+# --- Proposed freshwater kij_AQ (Burgoyne & Nielsen 2026, MC-3 alpha) ---
 def kij_aq_co2_proposed(T_K: float, salinity_molal: float = 0.0) -> float:
     """CO2: cubic in T(K), n=614, MAE=0.0119 (MC-3 alpha)."""
     return -1.5127 + 9.4980e-3*T_K - 2.1680e-5*T_K**2 + 1.8500e-8*T_K**3
@@ -380,7 +384,7 @@ def kij_aq_n2_proposed(T_K: float, salinity_molal: float = 0.0) -> float:
     return -1.6669 + 3.447873e-3 * T_K
 
 def kij_aq_h2_proposed(T_K: float, salinity_molal: float = 0.0) -> float:
-    """H2: rational, n=162, MAE=0.0400 (Paper 2 MC-3 coefficients)."""
+    """H2: rational, n=162, MAE=0.0400 (Burgoyne & Nielsen 2026 MC-3 coefficients)."""
     Tr = T_K / 33.145
     return (-14.6157 + Tr) / (3.5494 + 0.2230 * Tr)
 
@@ -513,7 +517,7 @@ KIJ_AQ_FUNCTIONS: Dict[str, Callable] = KIJ_AQ_PROPOSED
 
 
 # =============================================================================
-# Embedded Salinity BIP Parameters (Paper 2, all 8 gases)
+# Embedded Salinity BIP Parameters (Burgoyne & Nielsen 2026, all 8 gases)
 # =============================================================================
 # Form: kij(T,m) = kij_fw(T) + (a0 + a1*Tr + a2*Tr^2)*m  [+ (b0 + b1*Tr)*m^2 for CO2]
 # Tr = T/Tc for each gas. Fitted with MC-3 alpha + proposed kij_fw.
@@ -1446,7 +1450,7 @@ class SWBinaryVLE:
                     params=EMBEDDED_SALINITY_PARAMS_DROPIN[self.gas])
                 kij_aq = kij_fw + delta
             else:
-                # Use Paper 2 embedded delta (proposed framework)
+                # Use Burgoyne & Nielsen 2026 embedded delta (proposed framework)
                 delta = calc_embedded_delta_kij(self.gas, T_K, self.salinity)
                 kij_aq = kij_fw + delta
             x_gas = self._calc_x_with_kij(T_K, P_Pa, kij_aq)
@@ -1702,7 +1706,7 @@ def get_gas_gas_bip(gas_a: str, gas_b: str) -> float:
 
 _SW_KVALUE_PARAMS = {
     # Light gases: Cross form [a, b, c, d, e, f]
-    # Fitted using proposed Paper 2 freshwater kij_AQ forms (Feb 2026)
+    # Fitted using proposed Burgoyne & Nielsen 2026 freshwater kij_AQ forms (Feb 2026)
     'H2':     [6.4295, 25.5844, -0.5985, -50.0000, -0.0007, -3.2598],
     'CO2':    [-5.9974, 26.3804, -0.9380, -16.3941, 0.0607, 0.3688],
     'N2':     [1.4998, 36.4929, -0.6419, -50.0000, 0.0110, -0.6328],

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/pyrestoolbox/brine/brine.py

@@ -367,7 +367,15 @@ def brine_props(p: float = None, degf: float = None, wt: float = None, ch4_sat:
     # Superseded the Mao-Duan chain in 3.7.3; NaCl viscosities move by 0.34%
     # mean and 1.36% max, which is the route change rather than a regression.
     # Mao-Duan remains available as route='mao_duan' for reproducibility.
-    ub_tpm = _brine_viscosity(degk, Mpa, m=m)  # cP
+    ub_tpm = _brine_viscosity(degk, Mpa, m=m)  # cP, gas-free
+
+    # Dissolved-CH4 viscosity correction, applied from 3.7.4. This function had
+    # returned the GAS-FREE viscosity for a brine it reports as methane
+    # saturated, which understated the measured effect (Ostermann 1985: +3 to 6%)
+    # by the full size of that effect. mch4 is the dissolved molality from
+    # Eqs 4.15-4.18; the correction is a function of mole fraction.
+    ub_tpm = ub_tpm * ch4_viscosity_factor(
+        degk, ch4_mole_fraction_from_molality(mch4))
 
     bw = Bw  # rb/stb (dimensionless ratio, same in both unit systems)
     lden = rhobtpbch4  # sg (g/cm3)
@@ -1322,11 +1330,14 @@ def make_pvtw_table(*args, **kwargs):
 # (1969) AIChE J 15:272. Murphy & Gaines report viscosity ratios but no mole
 # fractions; taking x_H2S from Soreide-Whitson instead put it 4.5 to 9.2% high
 # against Burgess & Germann, and since a scales as 1/x it biased the coefficient
-# low by the same amount. Implied coefficients on the correct basis are 2.31,
-# 1.54, 1.98, 1.31 and 0.27; the adopted value is the mean of the four below
-# 35 degC. Was 1.64 (S&W basis), and 1.50 before the primary source was read.
-_IC_A_H2S = 1.79       # Murphy & Gaines Table IV on Burgess & Germann solubility
-_IC_B = 1.0134         # Islam-Carlson exponent
+# low by the same amount. Implied coefficients on the correct basis are 2.20,
+# 1.47, 1.89, 1.25 and 0.26; the adopted value is the mean of the four below
+# 35 degC. Exponent set to UNITY 2026-07-31: the prior 1.0134 was borrowed
+# from Islam-Carlson's CO2 fit and never fitted to H2S, and the five points
+# cannot distinguish exponents (RMS 1.740 vs 1.738 pp). Was 1.79 with the
+# borrowed exponent, 1.64 (S&W basis), and 1.50 before the primary source
+# was read.
+_H2S_A = 1.70          # Murphy & Gaines Table IV on Burgess & Germann solubility
 
 # Calabrese et al. (2019) Eq. 25 CO2 viscosity increment (JCED 64:3831):
 # ln(mu/mu_brine) = e1 * exp(-e2*(T_K/T0 - 1)) * x_CO2
@@ -1363,6 +1374,31 @@ _CH4_A = 1.71739196e-03
 _CH4_B = 1239.77535        # K
 _CH4_K = 1.52860547e-03    # half-saturation mole fraction
 
+_MW_H2O = 18.015268
+
+
+def ch4_viscosity_factor(T_K, x_ch4):
+    """Dissolved-CH4 viscosity multiplier, Ostermann (1985) SPE 14211.
+
+    x_ch4 is the dissolved mole fraction on the SALT-FREE basis, which is what
+    the S&W flash returns and what the fit was performed against.
+
+    Shared by SoreideWhitson and brine_props so the two cannot diverge. Before
+    3.7.4 brine_props applied no CH4 viscosity correction at all: it returned the
+    gas-free viscosity for a brine it simultaneously reported as methane
+    saturated, understating the measured effect by 3 to 6 percent.
+    """
+    if x_ch4 <= 0.0:
+        return 1.0
+    return 1.0 + _CH4_A * np.exp(_CH4_B / T_K) * x_ch4 / (_CH4_K + x_ch4)
+
+
+def ch4_mole_fraction_from_molality(m_ch4):
+    """Dissolved CH4 molality (mol/kg water) -> salt-free mole fraction."""
+    if m_ch4 <= 0.0:
+        return 0.0
+    return m_ch4 / (m_ch4 + 1000.0 / _MW_H2O)
+
 # HC component molecular weights for SG-based splitting
 _MW_C1 = 16.043
 _MW_C2 = 30.069
@@ -1404,8 +1440,10 @@ from pyrestoolbox.brine._lib_vle_engine import calc_gas_brine_equilibrium as _ca
 
 
 class SoreideWhitson:
-    """ Soreide-Whitson (1992) VLE model for multicomponent gas solubility in water/brine,
-        with Garcia/Plyasunov density corrections and calibrated viscosity corrections.
+    """ Soreide-Whitson VLE model for multicomponent gas solubility in water/brine,
+        using by default the refreshed BIP relationships of Burgoyne & Nielsen (2026),
+        Fluid Phase Equilibria 114824, doi:10.1016/j.fluid.2026.114824, with
+        mass-balance density corrections and calibrated viscosity corrections.
 
         Uses the S&W VLE engine for gas-brine equilibrium calculations, supporting
         multicomponent gas mixtures containing: C1, C2, C3, nC4, CO2, H2S, N2, H2.
@@ -1439,7 +1477,8 @@ class SoreideWhitson:
             sg: Gas specific gravity — used to estimate HC split among C1-C4 (default 0.65)
             metric: Boolean for units (True=metric, False=oilfield). Default False.
             cw_sat: If True, also calculate saturated compressibility (default False)
-            framework: VLE framework. 'proposed' (default, Soreide-Whitson 1992 re-fit),
+            framework: VLE framework. 'proposed' (default, the Burgoyne & Nielsen 2026
+                refreshed Soreide-Whitson BIPs, doi:10.1016/j.fluid.2026.114824),
                 'sw_original' (original 1992 published), or 'dropin' (fitted to PR-EOS with
                 brine-aware water alpha). Affects kij and ks correlations.
             salinity_method: How salinity enters the flash. 'gamma_phi' (default, Sechenov
@@ -1897,10 +1936,10 @@ class SoreideWhitson:
                 T_K = (degf - 32.0) / 1.8 + 273.15
                 vis_factor *= np.exp(_CAL_E1 * np.exp(-_CAL_E2 * (T_K / _CAL_T0 - 1.0)) * x_i)
             elif gas_upper == 'H2S':
-                vis_factor *= (1.0 + _IC_A_H2S * x_i ** _IC_B)
+                vis_factor *= (1.0 + _H2S_A * x_i)
             elif gas_upper == 'CH4':
                 T_K = (degf - 32.0) / 1.8 + 273.15
-                vis_factor *= 1.0 + _CH4_A * np.exp(_CH4_B / T_K) * x_i / (_CH4_K + x_i)
+                vis_factor *= ch4_viscosity_factor(T_K, x_i)
             # C2H6, N2, H2, C3H8, nC4H10: no correction (factor *= 1.0)
 
         vis_gas_brine = vis_base_cP * vis_factor
@@ -1986,4 +2025,37 @@ class SoreideWhitson:
         else:
             rho_p1_gcc = rho_brine_p1_gcc
 
-        self.Cf_usat = 1.0 - rho_gas_brine_gcc / rho_p1_gcc  # 1/bar
+        self.Cf_usat = 1.0 - rho_gas_brine_gcc / rho_p1_gcc  # 1/bar
+
+# ============================================================================
+# Single recommended entry point
+# ============================================================================
+
+_BRINE_METHODS = {
+    'co2': ('CO2_Brine_Mixture', 'Spycher-Pruess mutual solubility'),
+    'sw': ('SoreideWhitson', 'Soreide-Whitson multicomponent PR flash'),
+    'ch4': ('brine_props', 'Duan-Mao/McCain Ch. 4 methane solubility'),
+}
+
+
+def recommended_method(y_CO2=0.0, y_H2S=0.0, y_N2=0.0, y_H2=0.0, ch4_only=False):
+    """Which brine entry point to use for a given equilibrium gas, and why.
+
+    The three entry points are not interchangeable wrappers: they carry three
+    different SOLUBILITY models, and each is the better choice somewhere. This
+    returns the recommendation rather than hiding it, because which model
+    produced a number is part of the answer.
+
+    Returns (function name, solubility model, reason).
+    """
+    others = float(y_H2S) + float(y_N2) + float(y_H2)
+    if float(y_CO2) >= 1.0 - 1e-12 and others <= 1e-12:
+        return _BRINE_METHODS['co2'] + (
+            'pure CO2: Spycher-Pruess scores 3.7% MARE against Yan (2011) '
+            'where the S&W flash scores 8.3%',)
+    if ch4_only and others <= 1e-12 and float(y_CO2) <= 1e-12:
+        return _BRINE_METHODS['ch4'] + (
+            'methane only, in oilfield units, following McCain Ch. 4; use '
+            'SoreideWhitson instead if you want the S&W flash or any other gas',)
+    return _BRINE_METHODS['sw'] + (
+        'any gas mixture, or any gas other than pure CO2',)

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/pyrestoolbox/brine/pr_vphi.py

@@ -110,7 +110,7 @@ VSHIFT = {
     'CH4':  -0.109632,
     'CO2':  -0.037913,
     'H2S':  -0.078975,
-    'N2':   -0.155288,
+    'N2':   -0.176510,
     'H2':   -0.177625,
     'C2H6': -0.073142,
     # C3H8 added 2026-07-25. NOT fitted to a densimetric data set, because none
@@ -208,15 +208,19 @@ def _check(gas, T, P):
 
 def V2_inf_raw(gas, T, P, m_nacl=0.0):
     """Unshifted V2inf from the exact EOS relation, cm3/mol. No range guard."""
-    V = v_water_liquid(T, P, m_nacl)
+    # The two derivatives of V2 = -(dP/dn2)/(dP/dV) are evaluated ANALYTICALLY
+    # at n1 = 1, n2 = 0, V = v_w (closed forms: manuscript Appendix B).
+    # Switched from finite differences 2026-07-31; the difference evaluation
+    # agrees to 4e-5 cm3/mol, its truncation error.
+    v = v_water_liquid(T, P, m_nacl)
     kij = _SW.get_kij_aq(_SW_NAME.get(gas, gas), T, m_nacl)
-    # dP/dn2 at fixed T, V, n1: one-sided, because n2 = 0 is the boundary
-    h = 1e-7
-    dPdn2 = (_P_mix(T, V, 1.0, h, gas, kij, m_nacl)
-             - _P_mix(T, V, 1.0, 0.0, gas, kij, m_nacl)) / h
-    dV = V * 1e-7
-    dPdV = (_P_mix(T, V + dV, 1.0, 0.0, gas, kij, m_nacl)
-            - _P_mix(T, V - dV, 1.0, 0.0, gas, kij, m_nacl)) / (2.0 * dV)
+    a1, b1 = _ab('H2O', T, m_nacl)
+    a2, b2 = _ab(gas, T)
+    a12 = np.sqrt(a1 * a2) * (1.0 - kij)
+    D = v * v + 2.0 * b1 * v - b1 * b1
+    dPdn2 = (R * T / (v - b1) + R * T * b2 / (v - b1) ** 2
+             - (2.0 * a12 * D - 2.0 * a1 * b2 * (v - b1)) / D ** 2)
+    dPdV = -R * T / (v - b1) ** 2 + 2.0 * a1 * (v + b1) / D ** 2
     return -dPdn2 / dPdV
 
 
@@ -228,9 +232,10 @@ def V2_inf(gas, T, P, m_nacl=0.0, s=None):
         T: K, within T_MIN..T_MAX and above the water saturation pressure
         P: MPa, up to P_MAX
         m_nacl: NaCl molality. Enters the S&W water alpha and kij only; there is
-            NO salinity term on the volume shift, consistent with the project's
-            decision that a salt effect on V_phi is real but not applied
-            (see salt_effect_vphi.py). Leave at 0 for the calibrated behaviour.
+            NO salinity term on the volume shift. The salt effect on V_phi
+            is applied one level up (vphi_route, as a relative fraction since
+            2026-07-30), never here, so there is no double counting. Leave at
+            0 for the calibrated behaviour.
         s: override the dimensionless volume shift (for refitting).
 
     Returns:

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.4}/pyrestoolbox/brine/vphi_route.py

@@ -42,26 +42,27 @@ DEFAULT_ROUTE = "auto"
 
 
 # ---------------------------------------------------------------- salt effect
-# ADOPTED 2026-07-25. V_phi falls with salinity. This term is fixed ENTIRELY
-# from direct measurement, with no parameter fitted to any brine-density data:
-# Tiepel & Gubbins (1972) dilatometry, 15 apparent molar volumes over five gases
-# and four electrolytes, every one negative, fitted as
-#     dV = -a*m/(1 + b*m),  a = 0.5914, b = 0.0416   [cm3/mol, m in mol/kg]
-# giving -0.57 at 1 molal and -2.45 at 5 molal. Enns, O'Sullivan & Smith and
-# Heusler & Gaiser agree on sign and magnitude. It is the conservative end of
-# the available range, so it under-corrects rather than over-corrects.
+# ADOPTED 2026-07-25 as an absolute additive term; FORM CHANGED to RELATIVE
+# 2026-07-30. V_phi falls with salinity. The dimensionless fraction is fixed
+# ENTIRELY from direct measurement, with no parameter fitted to any
+# brine-density data: the seven Tiepel & Gubbins (1972) KCl dilatometry points,
+# fitted as (generated by brine_props/code/relative_salt_shift.py)
+#     V_eff = V_phi(T,P) * (1 + g(m)),  g(m) = -c*m/(1 + d*m)
+# giving -1.56% at 1 molal and -5.85% at 5 molal. Enns and O'Sullivan & Smith
+# agree on sign and size; the magnitude is known to about a factor of two
+# (O'Sullivan NaCl shallower, Tiepel KCl at this value).
 #
 # The EOS's own salinity response (via the S&W water alpha and kij) is only
 # -0.03 to -0.07 cm3/mol per molal, an order of magnitude too small and linear
 # where the measurements saturate. It is therefore not used, and this term
 # carries the whole effect, so there is no double counting.
-SALT_A, SALT_B = 0.5914, 0.0416
+SALT_C, SALT_D = 0.017009386, 0.090684173   # 8 sig digits, generated - do not retype
 
 
-def salt_shift(m_nacl: float) -> float:
-    """Change in V_phi from dissolved salt, cm3/mol. Gas-generic, always <= 0."""
+def salt_fraction(m_nacl: float) -> float:
+    """Fractional change in V_phi from dissolved salt. Gas-generic, always <= 0."""
     m = max(float(m_nacl), 0.0)
-    return -SALT_A * m / (1.0 + SALT_B * m)
+    return -SALT_C * m / (1.0 + SALT_D * m)
 
 
 def pr_applicable(gas: str, T: float, P: float) -> bool:
@@ -90,17 +91,18 @@ def V_phi(gas: str, T: float, P: float, route: str = DEFAULT_ROUTE,
         P: pressure, MPa
         route: 'auto' uses PR where calibrated and Plyasunov elsewhere; 'pr'
             forces PR and raises outside its box; 'plyasunov' forces A12inf.
-        m_nacl: NaCl molality. Applies the literature-anchored salt shift; it is
-            NOT passed to the EOS, whose own salinity response is far too small.
-        salt_effect: set False to recover the pre-2026-07-25 freshwater V_phi.
+        m_nacl: NaCl molality. Applies the literature-anchored relative salt
+            fraction; it is NOT passed to the EOS, whose own salinity response
+            is far too small.
+        salt_effect: set False to recover the freshwater V_phi.
     """
     if route not in ROUTES:
         raise ValueError(f"route must be one of {ROUTES}, got {route!r}")
-    dV = salt_shift(m_nacl) if salt_effect else 0.0
+    f = 1.0 + (salt_fraction(m_nacl) if salt_effect else 0.0)
     if route == "plyasunov":
-        return float(_plyasunov_V_phi(gas, T, P)) + dV
+        return float(_plyasunov_V_phi(gas, T, P)) * f
     if route == "pr":
-        return float(_pr.V2_inf(gas, T, P)) + dV
+        return float(_pr.V2_inf(gas, T, P)) * f
     base = (float(_pr.V2_inf(gas, T, P)) if pr_applicable(gas, T, P)
             else float(_plyasunov_V_phi(gas, T, P)))
-    return base + dV
+    return base * f