pyrestoolbox 3.8.0__tar.gz → 3.8.1__tar.gz

{pyrestoolbox-3.8.0 → pyrestoolbox-3.8.1}/Cargo.lock

@@ -152,7 +152,7 @@ dependencies = [
 
 [[package]]
 name = "pyrestoolbox-native"
-version = "3.8.0"
+version = "3.8.1"
 dependencies = [
  "pyo3",
  "rug",

{pyrestoolbox-3.8.0 → pyrestoolbox-3.8.1}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "pyrestoolbox-native"
-version = "3.8.0"
+version = "3.8.1"
 edition = "2021"
 license = "GPL-3.0-or-later"
 description = "Rust acceleration extensions for pyResToolbox"

{pyrestoolbox-3.8.0 → pyrestoolbox-3.8.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyrestoolbox
-Version: 3.8.0
+Version: 3.8.1
 Classifier: Programming Language :: Python :: 3
 Classifier: Programming Language :: Rust
 Classifier: License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)

{pyrestoolbox-3.8.0 → pyrestoolbox-3.8.1}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "maturin"
 
 [project]
 name = "pyrestoolbox"
-version = "3.8.0"
+version = "3.8.1"
 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.8.0 → pyrestoolbox-3.8.1}/pyrestoolbox/docs/changelist.rst

@@ -1,3 +1,7 @@
+Changelist in 3.8.1:
+
+- **NODAL: condensate that drops out in the tubing is treated as a saturated light liquid, not a dead oil.** The gas VLP marches gave dropped-out condensate the Beggs-Robinson dead-oil viscosity (Rs = 0), which is 2-5 times too viscous for a liquid in equilibrium with gas (up to about 16 cP at low temperature) and made the Woldesemayat-Ghajar BHP with dropout exceed the no-dropout value. The dropped liquid now takes the Beggs-Robinson live-oil viscosity at the Standing (1947) solution GOR at local pressure and temperature, with the bubble point at the local pressure (0.25-1.6 cP for 45-55 API condensate), in Python and Rust. Only gas wells with ``pr`` supplied (condensate dropout on) are affected; for 5 MMscf/d with CGR 10, 10 bbl/d water and 3000 psia reservoir pressure (6000 ft tubing), BHP moves from 877.6 to 850.4 psia for WG and from 635.6 to 652.9 psia for HB (its holdup correlation runs through the liquid viscosity number), with Gray and BB within 0.4%. With dropout on, WG BHP is now below the no-dropout value, as it should be with less liquid in the tubing.
+
 Changelist in 3.8.0:
 
 - **NODAL: water viscosity in the VLP marches comes from the library's brine model (documented values move).** The marches used a private water-viscosity formula that evaluated a natural-log fit in degC as a base-10 fit in degF, so water viscosity was 12-30x too low (0.054 cP at 100 degF against 0.68; floored at 0.01 cP above about 200 degF), understating liquid friction and holdup in every method. Both the Python and Rust marches now use the same gas-free brine viscosity ``brine_props`` reports (IAPWS-2008 water with IF97 density, Jones-Dole NaCl ratio, Kestin pressure factor), with the salinity taken from ``wsg`` via McCain's gamma_w = 1 + 0.695e-6 ppm; Rust carries a port of that chain matching Python to 1e-14. Wells without water are unchanged. With water the change is modest for small water rates (the documented gas examples with 10 bbl/d move by -0.1 to -45.6 psi; WG moves most because realistic viscosity puts the liquid Reynolds number below Chisholm's laminar threshold) and larger with more water (5 MMscf/d with 500 bbl/d: HB +9%, WG +5%). The pure-Python march is slower (about 15 ms per ``fbhp`` against 3); the Rust path is about 0.2-0.3 ms.

{pyrestoolbox-3.8.0 → pyrestoolbox-3.8.1}/pyrestoolbox/nodal/nodal.py

@@ -1048,8 +1048,17 @@ def _condensate_dropout(cgr, qg_mmscfd, p_avg, pr, osg, qw_bwpd, wsg):
 
 
 def _condensate_vis(pr, cgr_local, gsg, api, temp_f, p_avg, oil_vis):
+    """Viscosity (cP) of condensate dropped out in the tubing.
+
+    The dropped liquid is in equilibrium with the gas at the local pressure, so
+    it is a saturated live liquid: Beggs-Robinson live-oil viscosity at the
+    Standing (1947) solution GOR at local p and T, with Pb = p. Before 3.8.1 it
+    was given the dead-oil viscosity (Rs = 0), 2-5x too viscous for a light
+    condensate (up to ~16 cP at low temperature).
+    """
     if pr > 14.7 and cgr_local > 0.01:
-        return _oil_viscosity_full(gsg, api, temp_f, rsb=0.0, pb=14.7,
+        rs_sat = _standing_rs(gsg, p_avg, temp_f, api)
+        return _oil_viscosity_full(gsg, api, temp_f, rsb=rs_sat, pb=p_avg,
                                    press_psia=p_avg)
     return oil_vis
 

{pyrestoolbox-3.8.0 → pyrestoolbox-3.8.1}/src/vlp/march.rs

@@ -86,8 +86,11 @@ fn condensate_vis(
     pr: f64, cgr_local: f64, gsg: f64, api: f64, temp_f: f64,
     p_avg: f64, oil_vis: f64,
 ) -> f64 {
+    // Dropped condensate is a saturated live liquid: Beggs-Robinson at the
+    // Standing solution GOR at local p, T with Pb = p (mirrors Python _condensate_vis)
     if pr > 14.7 && cgr_local > 0.01 {
-        oil_viscosity_full(gsg, api, temp_f, 0.0, 14.7, p_avg, 1.0, 1.0)
+        let rs_sat = standing_rs(gsg, p_avg, temp_f, api);
+        oil_viscosity_full(gsg, api, temp_f, rs_sat, p_avg, p_avg, 1.0, 1.0)
     } else {
         oil_vis
     }