gea-program 0.85.0__tar.gz → 0.87.0__tar.gz

{gea_program-0.85.0/gea_program.egg-info → gea_program-0.87.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gea-program
-Version: 0.85.0
+Version: 0.87.0
 Summary: GEA-Program: downhole gauge monitoring - canonical measurement records, quality rules, drift evaluation with a re-fit change log, well-test validation, alarm management, model cards, store-and-forward, monthly SLA measurement, FAT/SAT protocols and a dashboard, generated from live objects in the client's reporting outline.
 Author: Daniel T. Murphy
 License: MPL-2.0
@@ -847,7 +847,7 @@ unsigned, the key hand-overs - and shows the reports as they were signed. A
 reader made from the workspace knows the site's key fingerprints, so it
 says whether a dossier is the site's; it allows no network request, and it
 leaves the trusted time-stamp to `verify_dossier.py`. For a client who wants
-it installed, every release also carries `GEA-Reader-Setup-0.85.0.exe`: GEA
+it installed, every release also carries `GEA-Reader-Setup-0.87.0.exe`: GEA
 Reader as a Windows program of its own - its own icon and taskbar button,
 Check with GEA Reader on a dossier's right-click menu, fingerprints the client
 saves kept on that computer, nothing of GEA-Program. `gea help reader`.
@@ -981,6 +981,47 @@ against it. The same contours drawn at the declared top in the 3D tab and
 carried in the dossier. Every figure a SIMULATION, said on every line. The
 **Pressure plume** tab, `gea plume`, `gea help plume`.
 
+### Fluid front
+
+Where the injected water itself has gone. The pressure plume is miles wide;
+the water, which has to push the formation's own out of the pore space, is
+hundreds of feet from the well. On a fine grid around each injection well
+the wells' steady flow is solved once per well (conjugate gradients, the
+sealing faults closed) and the injected water moved with it - the fraction
+of the thickness it fills, spreading under its own weight against the
+formation water's by the sharp-interface equations, so a denser water runs
+out along the base and a lighter one along the top. The tip, the half and
+the whole thickness against a piston of the same volume; the day the water
+reaches each producer, its share of that producer's stream and the TDS that
+predicts, beside the TDS the producer's own report measures and the share
+that implies; continued and shut in for five years; every injected barrel
+accounted - in place, produced back, across a fault, beyond the grid. Drawn
+in plan, in section and, at its depths, in the 3D tab. Every figure a
+SIMULATION, said on every line. The **Fluid front** tab, `gea front`,
+`gea help front`.
+
+### Calibration
+
+The plume and the front, run many times over the numbers nobody knows
+exactly and ranked against what is measured. Each uncertain number - the
+permeability, the thickness, the porosity, the compressibility - gets a
+range, declared with its source or a default spread named as one; each
+fault layer is open in half the runs and sealing in the other half. Every
+run is scored against the pressures in the ledger (a gauge in an
+observation well, an injector's or a producer's own bottomhole pressure,
+each compared in its shape through time) and the producers' measured TDS;
+the runs the measurements do not contradict are kept, and when none comes
+close the page says NO RUN FITS and to look outside the ranges. Every
+forecast the other pages give as one number - the area of review, each
+fault's rise and the share of kept runs that take it past the friction
+line, the front's radius, each producer's share and arrival - is given as
+a range: what the declarations alone allow beside what the measurements
+leave. And what to measure next, ranked: a falloff test, a gauge at a
+neighbour, a water analysis, each with how much it would narrow the area
+of review and the front. Run by a person, kept with the hash of its inputs,
+marked stale when they change. Every figure a SIMULATION. The
+**Calibration** tab, `gea calibrate`, `gea help calibration`.
+
 ### Plan vs actual and the survey tools
 
 Each well against its directional plan: at every survey station, how far
@@ -1645,7 +1686,7 @@ attaches the zips to the GitHub release of that tag.
 
 ### The Windows installer
 
-Every release also carries `GEA-Program-Setup-0.85.0.exe` (for 0.85.0): the
+Every release also carries `GEA-Program-Setup-0.87.0.exe` (for 0.87.0): the
 same Windows kit installed like any Windows program, for the signed-in user
 alone and with no administrator rights (`%LOCALAPPDATA%\Programs\GEA-Program`).
 It asks for the site workspace (an existing one such as `C:\site` is used as
@@ -1664,7 +1705,7 @@ it, removes the program and leaves the workspace. Silent, for a site's
 software deployment:
 
 ```
-GEA-Program-Setup-0.85.0.exe /VERYSILENT /SUPPRESSMSGBOXES /NORESTART /WORKSPACE="C:\site" /SITENAME="Pad 3" /TASKS="desktopicon"
+GEA-Program-Setup-0.87.0.exe /VERYSILENT /SUPPRESSMSGBOXES /NORESTART /WORKSPACE="C:\site" /SITENAME="Pad 3" /TASKS="desktopicon"
 ```
 
 The setup is built by Inno Setup 6 from `tools/installer.iss` in the same

{gea_program-0.85.0 → gea_program-0.87.0}/README.md

@@ -808,7 +808,7 @@ unsigned, the key hand-overs - and shows the reports as they were signed. A
 reader made from the workspace knows the site's key fingerprints, so it
 says whether a dossier is the site's; it allows no network request, and it
 leaves the trusted time-stamp to `verify_dossier.py`. For a client who wants
-it installed, every release also carries `GEA-Reader-Setup-0.85.0.exe`: GEA
+it installed, every release also carries `GEA-Reader-Setup-0.87.0.exe`: GEA
 Reader as a Windows program of its own - its own icon and taskbar button,
 Check with GEA Reader on a dossier's right-click menu, fingerprints the client
 saves kept on that computer, nothing of GEA-Program. `gea help reader`.
@@ -942,6 +942,47 @@ against it. The same contours drawn at the declared top in the 3D tab and
 carried in the dossier. Every figure a SIMULATION, said on every line. The
 **Pressure plume** tab, `gea plume`, `gea help plume`.
 
+### Fluid front
+
+Where the injected water itself has gone. The pressure plume is miles wide;
+the water, which has to push the formation's own out of the pore space, is
+hundreds of feet from the well. On a fine grid around each injection well
+the wells' steady flow is solved once per well (conjugate gradients, the
+sealing faults closed) and the injected water moved with it - the fraction
+of the thickness it fills, spreading under its own weight against the
+formation water's by the sharp-interface equations, so a denser water runs
+out along the base and a lighter one along the top. The tip, the half and
+the whole thickness against a piston of the same volume; the day the water
+reaches each producer, its share of that producer's stream and the TDS that
+predicts, beside the TDS the producer's own report measures and the share
+that implies; continued and shut in for five years; every injected barrel
+accounted - in place, produced back, across a fault, beyond the grid. Drawn
+in plan, in section and, at its depths, in the 3D tab. Every figure a
+SIMULATION, said on every line. The **Fluid front** tab, `gea front`,
+`gea help front`.
+
+### Calibration
+
+The plume and the front, run many times over the numbers nobody knows
+exactly and ranked against what is measured. Each uncertain number - the
+permeability, the thickness, the porosity, the compressibility - gets a
+range, declared with its source or a default spread named as one; each
+fault layer is open in half the runs and sealing in the other half. Every
+run is scored against the pressures in the ledger (a gauge in an
+observation well, an injector's or a producer's own bottomhole pressure,
+each compared in its shape through time) and the producers' measured TDS;
+the runs the measurements do not contradict are kept, and when none comes
+close the page says NO RUN FITS and to look outside the ranges. Every
+forecast the other pages give as one number - the area of review, each
+fault's rise and the share of kept runs that take it past the friction
+line, the front's radius, each producer's share and arrival - is given as
+a range: what the declarations alone allow beside what the measurements
+leave. And what to measure next, ranked: a falloff test, a gauge at a
+neighbour, a water analysis, each with how much it would narrow the area
+of review and the front. Run by a person, kept with the hash of its inputs,
+marked stale when they change. Every figure a SIMULATION. The
+**Calibration** tab, `gea calibrate`, `gea help calibration`.
+
 ### Plan vs actual and the survey tools
 
 Each well against its directional plan: at every survey station, how far
@@ -1606,7 +1647,7 @@ attaches the zips to the GitHub release of that tag.
 
 ### The Windows installer
 
-Every release also carries `GEA-Program-Setup-0.85.0.exe` (for 0.85.0): the
+Every release also carries `GEA-Program-Setup-0.87.0.exe` (for 0.87.0): the
 same Windows kit installed like any Windows program, for the signed-in user
 alone and with no administrator rights (`%LOCALAPPDATA%\Programs\GEA-Program`).
 It asks for the site workspace (an existing one such as `C:\site` is used as
@@ -1625,7 +1666,7 @@ it, removes the program and leaves the workspace. Silent, for a site's
 software deployment:
 
 ```
-GEA-Program-Setup-0.85.0.exe /VERYSILENT /SUPPRESSMSGBOXES /NORESTART /WORKSPACE="C:\site" /SITENAME="Pad 3" /TASKS="desktopicon"
+GEA-Program-Setup-0.87.0.exe /VERYSILENT /SUPPRESSMSGBOXES /NORESTART /WORKSPACE="C:\site" /SITENAME="Pad 3" /TASKS="desktopicon"
 ```
 
 The setup is built by Inno Setup 6 from `tools/installer.iss` in the same

{gea_program-0.85.0 → gea_program-0.87.0}/gea/__init__.py

@@ -123,7 +123,7 @@ from .profile_catalog import (
     read_pangaea_txt,
 )
 
-__version__ = "0.85.0"
+__version__ = "0.87.0"
 __all__ = [
     "aging_rate", "rate_psi_yr", "accuracy_psi",
     "Sensor", "SimulatorConfig", "DownholeEngine",

{gea_program-0.85.0 → gea_program-0.87.0}/gea/__main__.py

@@ -565,7 +565,7 @@ def main(argv=None) -> int:
     p_fl.add_argument("--workspace", type=str, default=None, help="the workspace folder; without one the self-test runs on labelled scenes")
     p_fl.add_argument("--well", type=str, default=None, help="the well a declaration, a report or a removal is for")
     p_fl.add_argument("--declare", action="store_true", help="declare the well's role, formation or density (with --source)")
-    p_fl.add_argument("--role", type=str, default=None, choices=("injection", "production"))
+    p_fl.add_argument("--role", type=str, default=None, choices=("injection", "production", "observation"), help="--declare: injection, production, or observation (a gauge or a water analysis, no flow)")
     p_fl.add_argument("--formation", type=str, default=None, help="--declare: the formation the well is completed in; alone with --porosity etc.: the formation to declare")
     p_fl.add_argument("--density-sg", type=float, default=None, help="--declare: the injected water's specific gravity")
     p_fl.add_argument("--density-ppg", type=float, default=None, help="--declare: the injected water's density in lb/gal")
@@ -603,6 +603,31 @@ def main(argv=None) -> int:
     p_pm.add_argument("--write", action="store_true", help="write reports/plume/plume.json and .md")
     p_pm.add_argument("--actor", type=str, default="cli")
     p_pm.add_argument("--json", action="store_true")
+    p_ff = sub.add_parser("front", help="the fluid front (SIMULATION): where the injected water itself has gone - the fraction of the thickness it fills on a fine grid around each injection well, moved by the wells' steady flow and spread by its own weight against the formation water's; the day it reaches each producer, its share of the stream and the TDS that predicts beside the measured one; every injected barrel accounted")
+    p_ff.add_argument("--workspace", type=str, default=None, help="the workspace folder; without one the self-test runs (the piston radius, the doublet's breakthrough, the tongue, a sealing fault)")
+    p_ff.add_argument("--well", type=str, default=None, help="one injection well's front (default: every injection well in the ledger)")
+    p_ff.add_argument("--cells", type=int, default=None, help="settings: cells per side of the grid, odd, 41 to 141 (default 81)")
+    p_ff.add_argument("--projection-days", type=int, default=None, help="settings: how far beyond the record to project, 365 to 7300 (default 1825)")
+    p_ff.add_argument("--write", action="store_true", help="write reports/front/front.json and .md")
+    p_ff.add_argument("--actor", type=str, default="cli")
+    p_ff.add_argument("--json", action="store_true")
+    p_cb = sub.add_parser("calibrate", help="the calibration (SIMULATION): the plume and the front run many times over the uncertain numbers - permeability, thickness, porosity, compressibility, each fault open or sealing - ranked against the measured pressures and the producers' water; the ranges that leaves, and which measurement would narrow them most")
+    p_cb.add_argument("--workspace", type=str, default=None, help="the workspace folder; without one the self-test runs (a formation whose true numbers are not the declared ones, found from its measurements)")
+    p_cb.add_argument("--run", action="store_true", help="run the ensemble now and keep the result")
+    p_cb.add_argument("--formation", type=str, default=None, help="--run: one formation; --range: the formation the range is for")
+    p_cb.add_argument("--runs", type=int, default=None, help="settings: the runs, 8 to 200 (default 48)")
+    p_cb.add_argument("--sigma-psi", type=float, default=None, help="settings: the pressure measurements' error, psi (default 25)")
+    p_cb.add_argument("--sigma-share", type=float, default=None, help="settings: the error of a share of injected water read from a TDS (default 0.05)")
+    p_cb.add_argument("--well-radius-ft", type=float, default=None, help="settings: the wellbore radius for a flowing well's own draw (default 0.354)")
+    p_cb.add_argument("--sample-sealing", type=str, default=None, choices=("yes", "no"), help="settings: whether each fault layer is open in half the runs and sealing in the other half (default yes)")
+    p_cb.add_argument("--range", dest="range_param", type=str, default=None, choices=("k_md", "thickness_ft", "porosity", "c_t_per_psi"), help="declare an uncertain number's range for --formation: --low, --high, --source (or --clear)")
+    p_cb.add_argument("--low", type=float, default=None)
+    p_cb.add_argument("--high", type=float, default=None)
+    p_cb.add_argument("--clear", action="store_true", help="--range: back to the default spread")
+    p_cb.add_argument("--source", type=str, default="", help="where a declared range comes from")
+    p_cb.add_argument("--write", action="store_true", help="write reports/calibration/calibration.json and .md from the kept result")
+    p_cb.add_argument("--actor", type=str, default="cli")
+    p_cb.add_argument("--json", action="store_true")
     p_os = sub.add_parser("seismicity", help="operations and seismicity: each well's drilling, completion and disposal windows against the earthquake catalogue around it - the events during and after each, the levels a regulator's protocol would reach, and the rate during operations against the rest of the record; the Mahoning Valley worked example")
     p_os.add_argument("--workspace", type=str, default=None)
     p_os.add_argument("--scene", type=str, default=None, choices=("workspace", "mahoning"), help="your wells (the default with --workspace) or the public worked example")
@@ -2738,6 +2763,75 @@ def main(argv=None) -> int:
         else:
             print(_PM.text(p))
         return 0
+    elif a.cmd == "front":
+        import json as _json
+        from . import front as _FF
+        from .cli import _safe_stdout
+        from .workspace import Workspace, WorkspaceError
+        _safe_stdout()
+        if not a.workspace:
+            r = _FF.selftest()
+            print(_json.dumps(r, default=str) if a.json else "\n".join(f"{'PASS' if c['pass'] else 'FAIL'}  {c['check']}: {c['value']}" for c in r["checks"]) + f"\n{'OK' if r['ok'] else 'FAILED'} - {r['label']}")
+            return 0 if r["ok"] else 1
+        try:
+            ws = Workspace(a.workspace)
+        except WorkspaceError as e:
+            raise SystemExit(f"front: {e}")
+        try:
+            if a.cells is not None or a.projection_days is not None:
+                st_ = _FF.set_settings(ws, a.actor, a.cells, a.projection_days)
+                print("settings: " + ", ".join(f"{k} {v}" for k, v in st_.items() if k in ("cells", "projection_days")))
+                return 0
+            if a.write:
+                from . import __version__ as _ver86
+                r = _FF.write_report(ws, _ver86, a.actor)
+                print(f"written: {', '.join(r['files'])} under {r['folder']}")
+                p = r["page"]
+            else:
+                p = _FF.page(ws, a.well)
+        except (_FF.FrontError, WorkspaceError, ValueError, OSError) as e:
+            raise SystemExit(f"front: {e.args[0] if e.args else e}")
+        print(_json.dumps(p, default=str) if a.json else _FF.text(p))
+        return 0
+    elif a.cmd == "calibrate":
+        import json as _json
+        from . import calibration as _CB
+        from .cli import _safe_stdout
+        from .workspace import Workspace, WorkspaceError
+        _safe_stdout()
+        if not a.workspace:
+            r = _CB.selftest()
+            print(_json.dumps(r, default=str) if a.json else "\n".join(f"{'PASS' if c['pass'] else 'FAIL'}  {c['check']}: {c['value']}" for c in r["checks"]) + f"\n{'OK' if r['ok'] else 'FAILED'} - {r['label']} ({r['runs']} runs)")
+            return 0 if r["ok"] else 1
+        try:
+            ws = Workspace(a.workspace)
+        except WorkspaceError as e:
+            raise SystemExit(f"calibrate: {e}")
+        try:
+            if a.range_param:
+                if not a.formation:
+                    raise SystemExit("calibrate: --range needs --formation")
+                _CB.set_range(ws, a.actor, a.formation, a.range_param, None if a.clear else a.low, None if a.clear else a.high, a.source)
+                rg_ = _CB.page(ws)["ranges"].get(a.formation, {}).get(a.range_param)
+                print(f"{a.formation} {a.range_param}: " + (f"{rg_['low']:g} to {rg_['high']:g} ({rg_['kind']}: {rg_['source']})" if rg_ else "kept; the formation has no permeability, thickness and porosity declared yet"))
+                return 0
+            if any(v is not None for v in (a.runs, a.sigma_psi, a.sigma_share, a.well_radius_ft, a.sample_sealing)):
+                st_ = _CB.set_settings(ws, a.actor, a.runs, a.sigma_psi, a.sigma_share, (None if a.sample_sealing is None else a.sample_sealing == "yes"), a.well_radius_ft)
+                print("settings: " + ", ".join(f"{k} {v}" for k, v in st_.items() if k in ("runs", "sigma_psi", "sigma_share", "well_radius_ft", "sample_sealing")))
+                if not a.run:
+                    return 0
+            if a.run:
+                r = _CB.run(ws, a.actor, a.formation)
+                print(f"run: {r['seconds']} s; " + "; ".join(f"{f['formation']} {f['status']}" + (f" ({f['n_kept']} of {f['n_ok']} kept)" if f.get('n_kept') is not None else "") for f in r["formations"]))
+            if a.write:
+                from . import __version__ as _ver87
+                r = _CB.write_report(ws, _ver87, a.actor)
+                print(f"written: {', '.join(r['files'])} under {r['folder']}")
+            p = _CB.page(ws)
+        except (_CB.CalibrationError, WorkspaceError, ValueError, OSError) as e:
+            raise SystemExit(f"calibrate: {e.args[0] if e.args else e}")
+        print(_json.dumps(p, default=str) if a.json else _CB.text(p))
+        return 0
     elif a.cmd == "seismicity":
         import json as _json
         from . import op_seismicity as _OS