gea-program 0.35.0__tar.gz → 0.37.0__tar.gz

{gea_program-0.35.0/gea_program.egg-info → gea_program-0.37.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gea-program
-Version: 0.35.0
+Version: 0.37.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
@@ -648,6 +648,34 @@ against it; `gea seedlink --config ... --hello` prints what a server says of
 itself; `gea help seedlink`. A TexNet station is served by EarthScope at
 `rtserve.iris.washington.edu:18000`, network `TX`.
 
+### Economics
+
+The Economics tab puts each well's days on its drilling record, feet
+drilled and feet per day beside what it cost by the operator's own declared
+day rates, items and AFE - cost per foot, per day and by section, never
+divided while the rates leave time uncovered - and sets each well against
+the U.S. Energy Information Administration's daily WTI and Henry Hub spot
+prices while it was drilled. The prices ship with the program and refresh
+on request; any other series (a rig count) is imported with its source.
+`gea economics`, `gea help economics`.
+
+### The well in 3D
+
+The 3D tab draws each well's path from its directional survey by the
+minimum-curvature method, below sea level from its declared rig floor, among
+the mapped formations - interpreted horizon surfaces, or a formation-top
+map's contours at their elevations - with its picks on it. Colour the wells
+by formation, inclination, dogleg severity, or a drilling recorder's gamma,
+rate of penetration or weight on bit along the bit depth; replay the bit
+going down the hole from the recorder's bit depth or the survey stations'
+own times. Turn, pan and zoom by pointer, touch or keys. The Volve field
+(Equinor's open data: three original surveys, the field's formation picks,
+five mapped horizons) ships as the worked example, checked against itself.
+Until then the workspace's formation-top contours are drawn on their own.
+A workspace well is drawn once it has a position, a reference elevation and
+a survey: `gea trajectory` attaches one, verbatim, with its azimuth
+reference stated. `gea well3d`, `gea help well3d`.
+
 ### The control room
 
 The Control room tab opens the wall: a screen read from across the room, one

{gea_program-0.35.0 → gea_program-0.37.0}/README.md

@@ -609,6 +609,34 @@ against it; `gea seedlink --config ... --hello` prints what a server says of
 itself; `gea help seedlink`. A TexNet station is served by EarthScope at
 `rtserve.iris.washington.edu:18000`, network `TX`.
 
+### Economics
+
+The Economics tab puts each well's days on its drilling record, feet
+drilled and feet per day beside what it cost by the operator's own declared
+day rates, items and AFE - cost per foot, per day and by section, never
+divided while the rates leave time uncovered - and sets each well against
+the U.S. Energy Information Administration's daily WTI and Henry Hub spot
+prices while it was drilled. The prices ship with the program and refresh
+on request; any other series (a rig count) is imported with its source.
+`gea economics`, `gea help economics`.
+
+### The well in 3D
+
+The 3D tab draws each well's path from its directional survey by the
+minimum-curvature method, below sea level from its declared rig floor, among
+the mapped formations - interpreted horizon surfaces, or a formation-top
+map's contours at their elevations - with its picks on it. Colour the wells
+by formation, inclination, dogleg severity, or a drilling recorder's gamma,
+rate of penetration or weight on bit along the bit depth; replay the bit
+going down the hole from the recorder's bit depth or the survey stations'
+own times. Turn, pan and zoom by pointer, touch or keys. The Volve field
+(Equinor's open data: three original surveys, the field's formation picks,
+five mapped horizons) ships as the worked example, checked against itself.
+Until then the workspace's formation-top contours are drawn on their own.
+A workspace well is drawn once it has a position, a reference elevation and
+a survey: `gea trajectory` attaches one, verbatim, with its azimuth
+reference stated. `gea well3d`, `gea help well3d`.
+
 ### The control room
 
 The Control room tab opens the wall: a screen read from across the room, one

{gea_program-0.35.0 → gea_program-0.37.0}/gea/__init__.py

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

{gea_program-0.35.0 → gea_program-0.37.0}/gea/__main__.py

@@ -475,6 +475,41 @@ def main(argv=None) -> int:
     p_ck.add_argument("--note", type=str, default="", help="--confirm: what it was checked against (required)")
     p_ck.add_argument("--actor", type=str, default="cli")
 
+    p_ec = sub.add_parser("economics", help="each well's time, footage and declared cost against the market while it was drilled; declare rates, items and an AFE; refresh or import price series")
+    p_ec.add_argument("--workspace", type=str, required=True)
+    p_ec.add_argument("--json", action="store_true")
+    p_ec.add_argument("--well", type=str, default=None, help="the well a rate, item or AFE is for")
+    p_ec.add_argument("--rate", type=float, default=None, help="a day rate, USD/day, over --from/--to (open when left out)")
+    p_ec.add_argument("--item", type=float, default=None, help="a one-off cost, USD (a bit, a cement job), optionally --section N")
+    p_ec.add_argument("--afe", type=float, default=None, help="the well's AFE, USD")
+    p_ec.add_argument("--what", type=str, default="", help="what the rate or item pays for")
+    p_ec.add_argument("--section", type=int, default=None)
+    p_ec.add_argument("--from", dest="date_from", type=str, default=None, help="rate start, ISO date or date-time on the recorder's clock")
+    p_ec.add_argument("--to", dest="date_to", type=str, default=None, help="rate end")
+    p_ec.add_argument("--source", type=str, default="", help="where the number came from (contract, invoice, AFE number)")
+    p_ec.add_argument("--refresh-prices", action="store_true", help="fetch the latest EIA spot-price copies into the workspace")
+    p_ec.add_argument("--import-series", type=str, default=None, help="a CSV of dates and values to keep beside the prices (a rig count)")
+    p_ec.add_argument("--name", type=str, default="")
+    p_ec.add_argument("--unit", type=str, default="")
+    p_ec.add_argument("--actor", type=str, default="cli")
+
+    p_tj = sub.add_parser("trajectory", help="attach a well's directional survey (WITSML trajectory or a table of MD, Inc, Azi), kept verbatim beside the well, so the 3D tab can draw its path")
+    p_tj.add_argument("--workspace", type=str, required=True)
+    p_tj.add_argument("--well", type=str, required=True)
+    p_tj.add_argument("--file", type=str, required=True)
+    p_tj.add_argument("--azimuth-reference", type=str, required=True, choices=("grid", "true", "magnetic"))
+    p_tj.add_argument("--grid", type=str, default=None, help="grid azimuths: the grid, e.g. TX_CENTRAL or UTM14N")
+    p_tj.add_argument("--declination", type=float, default=None, help="magnetic azimuths: the declination used, degrees east positive")
+    p_tj.add_argument("--md-unit", type=str, default="ft", choices=("ft", "m"))
+    p_tj.add_argument("--depth-reference", type=str, default=None, help="KB, DF or GL: must match the well's declared reference elevation")
+    p_tj.add_argument("--source", type=str, required=True, help='where the survey came from, e.g. "RRC imaged records, directional survey"')
+    p_tj.add_argument("--actor", type=str, default="cli")
+
+    p_3d = sub.add_parser("well3d", help="the 3D scene in words: each well's path, how it was checked, the surfaces and the replay (the 3D tab draws it)")
+    p_3d.add_argument("--scene", type=str, default="volve", choices=("volve", "workspace"))
+    p_3d.add_argument("--workspace", type=str, default=None)
+    p_3d.add_argument("--json", action="store_true")
+
     p_wl = sub.add_parser("wall", help="the control-room board as text: a line per well, worst first, its state and why (the same board the wall screen shows)")
     p_wl.add_argument("--workspace", type=str, default=None, help="the workspace folder, e.g. C:\\site")
     p_wl.add_argument("--json", action="store_true")
@@ -1937,6 +1972,71 @@ def main(argv=None) -> int:
         else:
             print(_SC.text(r))
         return 0 if r["ready"] else 2
+    elif a.cmd == "economics":
+        import json as _json
+        from . import economics as _EC
+        from .cli import _safe_stdout
+        from .workspace import Workspace, WorkspaceError
+        _safe_stdout()
+        try:
+            ws = Workspace(a.workspace)
+        except WorkspaceError as e:
+            raise SystemExit(f"economics: {e}")
+        try:
+            if a.rate is not None:
+                r = _EC.set_rate(ws, a.well or "", a.actor, a.rate, a.what, a.source, a.date_from, a.date_to)
+                print(f"{a.well}: {r['usd_per_day']:,.2f} USD/day {r['what']} from {r['from'] or 'open'} to {r['to'] or 'open'}")
+                return 0
+            if a.item is not None:
+                r = _EC.add_item(ws, a.well or "", a.actor, a.item, a.what, a.source, a.section)
+                print(f"{a.well}: {r['usd']:,.2f} USD {r['what']}")
+                return 0
+            if a.afe is not None:
+                r = _EC.set_afe(ws, a.well or "", a.actor, a.afe, a.source)
+                print(f"{a.well}: AFE {r['usd']:,.2f} USD ({r['source']})")
+                return 0
+            if a.refresh_prices:
+                for r in _EC.refresh_prices(ws, a.actor):
+                    print(f"{r['id']}: {r['status']}" + (f" to {r['as_of']} ({r['rows']} days)" if r.get('as_of') else f" - {r.get('why', '')}"))
+                return 0
+            if a.import_series:
+                with open(a.import_series, "rb") as f:
+                    r = _EC.import_series(ws, a.actor, a.name, a.unit, a.source, f.read())
+                print(f"{r['name']}: {r['rows']} values, {r['first']} to {r['as_of']} ({r['unit']})")
+                return 0
+        except (ValueError, KeyError, OSError, WorkspaceError) as e:
+            raise SystemExit(f"economics: {e.args[0] if e.args else e}")
+        e = _EC.economics(ws)
+        print(_json.dumps(e, indent=1) if a.json else _EC.text(e))
+        return 0
+    elif a.cmd in ("trajectory", "well3d"):
+        import json as _json
+        from . import well3d as _W3
+        from .cli import _safe_stdout
+        from .workspace import Workspace, WorkspaceError
+        _safe_stdout()
+        ws = None
+        if a.cmd == "trajectory" or a.scene == "workspace":
+            if not a.workspace:
+                raise SystemExit(f"{a.cmd} needs --workspace, e.g. --workspace C:\\site")
+            try:
+                ws = Workspace(a.workspace)
+            except WorkspaceError as e:
+                raise SystemExit(f"{a.cmd}: {e}")
+        if a.cmd == "trajectory":
+            try:
+                with open(a.file, "rb") as f:
+                    content = f.read()
+                r = _W3.attach_trajectory(ws, a.well, a.file, content, a.actor, a.azimuth_reference, a.md_unit, a.depth_reference, a.grid,
+                                          a.declination, a.source)
+            except (ValueError, OSError, WorkspaceError) as e:
+                raise SystemExit(f"trajectory: {e}")
+            print(f"{a.well}: {r['stations']} stations, {r['md_range_m'][0]:,.1f} to {r['md_range_m'][1]:,.1f} m MD, azimuths {r['azimuth_reference']}"
+                  + (f" ({r['grid']})" if r.get('grid') else "") + f"; kept as wells/{a.well}/trajectory/{r['file']}")
+            return 0
+        sc = _W3.scene(ws, a.scene)
+        print(_json.dumps(sc, indent=1) if a.json else _W3.text(sc))
+        return 0
     elif a.cmd in ("wall", "handover"):
         import json as _json
         from . import control_room as _CR

{gea_program-0.35.0 → gea_program-0.37.0}/gea/acceptance_tests.py

@@ -8613,6 +8613,363 @@ def section_cd_control_room(tmp: str) -> None:
        + (f" [failed: {', '.join(failed3)}]" if failed3 else ""))
 
 
+
+def section_ce_well3d(tmp: str) -> None:
+    """Section CE - the well in 3D: the minimum-curvature method against a known arc and against the survey companies'
+    own positions, the shipped Volve excerpt checked against itself (frame, picks, surfaces, provenance), a workspace
+    well drawn from an attached survey among its formation-top contours with its recorder replay, the API and the tab."""
+    import base64 as _b64
+    import gzip as _gz
+    import hashlib as _hl
+    import http.cookiejar
+    import math as _m
+    import subprocess as _sp
+    import urllib.request
+    import urllib.error
+    from . import edr as E
+    from . import geodesy as G
+    from . import helplib as H
+    from . import shapefile as SF
+    from . import well3d as W
+    from .workspace import Workspace
+    from .service import Service, Users
+    pkg = Path(__file__).parent
+    d = Path(tmp, "ce"); d.mkdir()
+
+    R = 1000.0
+    q = W.min_curvature([0.0, _m.pi * R / 2], [0.0, 90.0], [0.0, 0.0])
+    p_arc = W.Path([0.0, _m.pi * R / 2], [0.0, 90.0], [0.0, 0.0], 0.0, 0.0, 0.0)
+    mid = p_arc.at(_m.pi * R / 4)
+    tr4 = W.read_witsml_trajectory(_gz.decompress((pkg / "reference" / "volve_traj_F-4.xml.gz").read_bytes()).decode())[0]
+    tr10 = W.read_witsml_trajectory(_gz.decompress((pkg / "reference" / "volve_traj_F-10.xml.gz").read_bytes()).decode())[0]
+    f5 = [t for n in ("36in", "17.5in", "12.25in", "8.5in") for t in W.read_witsml_trajectory(_gz.decompress((pkg / "reference" / f"volve_traj_F-5_{n}.xml.gz").read_bytes()).decode())]
+    rep5, _ = W.replay_from_stations(f5)
+    rep4, note4 = W.replay_from_stations([tr4])
+    csv_ok = "MD (ft),Inc (deg),Azi (deg)\n0,0,0\n1000,0,0\n2000,10,90\n3000,30,90\n4000,45,95\n5750,60,100\n"
+    refused = []
+    for bad in ("Depth,Angle\n0,0\n", "MD,Inc,Azi\n0,0,0\n0,1,2\n", "MD,Inc,Azi\n0,0,0\n100,x,2\n", "MD,Inc,Azi\n0,0,0\n100,200,2\n"):
+        try:
+            W.read_survey_csv(bad)
+        except ValueError:
+            refused.append(True)
+    st_csv = W.read_survey_csv(csv_ok, "ft")
+    c1 = {
+        "arc": abs(q["n"][1] - R) < 1e-9 and abs(q["v"][1] - R) < 1e-9 and abs(mid[1] - R * (1 - _m.cos(_m.pi / 4))) < 1e-6
+               and abs(mid[2] - R * _m.sin(_m.pi / 4)) < 1e-6 and abs(q["dls"][1] - 90.0 / (_m.pi * R / 2) * 30.0) < 1e-9,
+        "units": tr4["azi_ref"] == "true north" and abs(tr4["stations"][40]["incl"] - _m.degrees(0.446106147512686)) < 1e-9 and tr10["grid_cor_deg"] == -0.95
+                 and len(tr4["stations"]) == 87 and abs(st_csv[-1]["md"] - 5750 * 0.3048) < 1e-9,
+        "company": all(x is not None and x < 0.01 for x in (W.file_position_miss(tr4["stations"]), W.file_position_miss(tr10["stations"]),
+                                                            W.file_position_miss(f5[-1]["stations"]))),
+        "csv": len(refused) == 4,
+        "replay": rep5 is not None and len(rep5["t"]) == 100 and rep5["t"][0].startswith("2007-11-16") and rep5["t"][-1].startswith("2008-07-29")
+                  and rep4 is None and "1 distinct time" in note4 and all(m0 <= m1 or t0 < t1 for (t0, m0), (t1, m1) in zip(zip(rep5["t"], rep5["md"]), zip(rep5["t"][1:], rep5["md"][1:]))),
+    }
+    failed1 = [k for k, v in c1.items() if not v]
+    ok(not failed1,
+       "CE1 the minimum-curvature method: a quarter circle of 1,000 m radius ends 1,000 m out and 1,000 m down, and its middle lies on the arc "
+       "(292.9 m out, 707.1 m down), not on a chord; the WITSML reader takes radians and degrees alike; the paths computed from the Volve files "
+       "reproduce the TVD and displacements the survey companies wrote in them to under a centimetre; a survey table without its columns, with a "
+       "repeated depth or a bad row refused; the replay only from stations measured while drilling (F-5: 100 steps, Nov 2007 to Jul 2008; F-4's "
+       "one-time resurvey not replayed)"
+       + (f" [failed: {', '.join(failed1)}]" if failed1 else ""))
+
+    sc = W.volve_scene()
+    wl = {w["id"]: w for w in sc["wells"]}
+    picks = W.volve_picks()
+    dms = lambda dd, mm, ss: dd + mm / 60 + ss / 3600
+    frame = []
+    for wn, la, lo in (("NO 15/9-F-1 C", dms(58, 26, 29.907), dms(1, 53, 14.708)), ("NO 15/9-F-11 B", dms(58, 26, 29.960), dms(1, 53, 14.870))):
+        sb = next(p for p in picks if p["WELL"] == wn and p["PICKS"] == "Seabed")
+        a = G.unproject(float(sb["EASTING"]), float(sb["NORTHING"]), G.utm(31, True, "INTL1924"))
+        b = G.unproject(float(sb["EASTING"]), float(sb["NORTHING"]), G.utm(31, True, "WGS84"))
+        frame.append((G.geodesic_m(a[0], a[1], la, lo)["distance_m"], G.geodesic_m(b[0], b[1], la, lo)["distance_m"]))
+    prov = json.loads((pkg / "reference" / "volve_3d.provenance.json").read_text(encoding="utf-8"))
+    hashes_ok = all(_hl.sha256((pkg / "reference" / k).read_bytes()).hexdigest() == v["sha256"] for k, v in prov["outputs"].items())
+    pts = prov["outputs"]["volve_surfaces.json.gz"]["points"]
+    # the repository stores every text file with LF line ends (.gitattributes: text=auto eol=lf), so a package file holding
+    # CRLF would reach a checkout - and the published wheel - as different bytes from the ones hashed here
+    crlf = []
+    for dp, dn, fn in os.walk(pkg):
+        dn[:] = [x for x in dn if x not in ("__pycache__", "catalog_operator")]
+        for f_ in fn:
+            b_ = Path(dp, f_).read_bytes()
+            if b"\0" not in b_[:8000] and b"\r\n" in b_:
+                crlf.append(os.path.relpath(os.path.join(dp, f_), pkg))
+    c2 = {
+        "frame": all(a < 0.1 and 140 < b < 160 for a, b in frame) and "ED50" in sc["frame"]["crs"],
+        "picks": all(wl[k]["pick_miss_max_m"] < 0.3 and abs(wl[k]["ref_elev_m"] - 54.9) < 1e-6 for k in ("F-4", "F-5", "F-10"))
+                 and abs(wl["F-10"]["td_md"] - 5300.07) < 0.01 and len(sc["picks_only"]) == 31,
+        "azimuth": wl["F-4"]["azimuth"]["applied"] == "as grid" and "14 m" in wl["F-4"]["azimuth"]["note"] and wl["F-5"]["azimuth"]["file_reference"] == "grid north",
+        "surfaces": len(sc["surfaces"]) == 5 and all(s["affine_miss_m"] < 0.01 and s["tie_median_m"] < 7.0 and len(s["ties"]) >= 23 for s in sc["surfaces"])
+                    and all(pts[s["id"]] == [s["n_source_points"], s["n_points"]] for s in sc["surfaces"])
+                    and [s["id"] for s in sc["surfaces"]] == ["Ty", "SHETLAND", "BCU", "Hugin_Fm_Top", "Hugin_Base"],
+        "line_ends": not crlf and all(k.endswith(".gz") or b"\r\n" not in (pkg / "reference" / k).read_bytes() for k in prov["outputs"]),
+        "provenance": hashes_ok and "commercial purposes" in prov["license"] and "may not be sold" in prov["license"] and len(prov["sources"]) == 12
+                      and all(len(v["commit"]) == 40 for v in prov["sources"].values()) and "Equinor ASA" in sc["attribution"] and "not sold" in sc["attribution"],
+    }
+    failed2 = [k for k, v in c2.items() if not v]
+    ok(not failed2,
+       "CE2 the Volve excerpt checked against itself: its eastings and northings are ED50 zone 31N (the slots unproject within 0.1 m of the "
+       "operator's own log headers, about 150 m off if read as WGS84); the three paths land on Equinor's picks to under 0.3 m from a 54.9 m rig "
+       "floor; 15/9-F-4's azimuths read as grid, with the 14 m it would miss by otherwise stated; five surfaces drawn as delivered, every fourth "
+       "line, each within a few metres of its picks at the median; every shipped file hashed in its provenance with the pinned public copy it "
+       "came from (and no package file with CRLF line ends, which a checkout would change from the bytes hashed), the licence (commercial use, no sale) and the attribution on the page"
+       + (f" [failed: {', '.join(failed2)}]" if failed2 else ""))
+
+    ev = E.selftest_events(str(d))
+    wsp = str(d / "ws"); ws = Workspace.create(wsp, "Pad 3", actor="tester")
+    wb = ws.add_well_file(ev["path"], display="Well B", actor="tester"); wc = ws.add_well_file(ev["path"], display="Well C", actor="tester")
+    ws.set_well(wb["id"], surface={"lat": 31.42, "lon": -104.10, "datum": "NAD83", "elevation": {"ft": 3000.0, "datum": "KB", "source": "test plat"}}, actor="tester")
+    lat1, lon1 = 31.42, -104.10
+    kmlon = 1.0 / (111.195 * _m.cos(_m.radians(lat1)))
+    SF.write(str(d / "Top_A"), "PolyLine", [([[(lon1 + (i - 1.5) * kmlon + 0.0005 * k, lat1 - 0.05 + 0.001 * k) for k in range(100)]], {"ElevTop_ft": -1500.0 - 100.0 * i}) for i in range(4)],
+             [("ElevTop_ft", "F", 10, 1)], SF.NAD83_GEOGCS)
+    SF.write(str(d / "Faults"), "PolyLine", [([[(lon1 - 0.02 + 0.001 * j, lat1) for j in range(40)]], {"TYPE": "normal fault"})], [("TYPE", "C", 20, 0)], SF.NAD83_GEOGCS)
+    ws.add_layer(str(d / "Top_A.shp"), "Top A", actor="tester", formation="A", role="top", citation="test")
+    ws.add_layer(str(d / "Faults.shp"), "Faults", actor="tester", citation="test")
+    sc0 = W.workspace_scene(ws)                                   # no survey yet: the formation tops alone
+    att_refused = 0
+    for kw in ({"azimuth_reference": ""}, {"azimuth_reference": "grid"}, {"azimuth_reference": "magnetic"}, {"azimuth_reference": "true", "source": " "}):
+        k = {"azimuth_reference": "true", "source": "s", **kw}
+        try:
+            W.attach_trajectory(ws, wb["id"], "s.csv", csv_ok.encode(), "tester", k["azimuth_reference"], "ft", None, None, None, k["source"])
+        except ValueError:
+            att_refused += 1
+    ws.refresh_dashboard(actor="tester")
+    Users(str(Path(wsp, "users.json"))).add("v", "viewer-pass-1", "viewer"); Users(str(Path(wsp, "users.json"))).add("op", "op-pass-123", "operator")
+    svc = Service(wsp, port=0, scheduler=False).start()
+    base_url = svc.url.rstrip("/")
+
+    def opener(name, pw):
+        o = urllib.request.build_opener(urllib.request.HTTPCookieProcessor(http.cookiejar.CookieJar()))
+        rq = urllib.request.Request(base_url + "/api/login", method="POST", data=json.dumps({"name": name, "password": pw}).encode())
+        rq.add_header("Content-Type", "application/json"); rq.add_header("X-GEA-Action", "1"); o.open(rq, timeout=30)
+        return o
+
+    def call(o, method, path, body=None):
+        rq = urllib.request.Request(base_url + path, method=method, data=(json.dumps(body).encode() if body is not None else None))
+        if method == "POST":
+            rq.add_header("Content-Type", "application/json"); rq.add_header("X-GEA-Action", "1")
+        try:
+            with o.open(rq, timeout=120) as rr_:
+                return rr_.status, json.loads(rr_.read())
+        except urllib.error.HTTPError as e:
+            return e.code, json.loads(e.read())
+    try:
+        V, O = opener("v", "viewer-pass-1"), opener("op", "op-pass-123")
+        st_l, ls = call(V, "GET", "/api/3d/scenes")
+        st_v, vs = call(V, "GET", "/api/3d/scene/volve")
+        st_n, _ = call(V, "GET", "/api/3d/scene/nope")
+        body = {"filename": "survey.csv", "content_b64": _b64.b64encode(csv_ok.encode()).decode(), "azimuth_reference": "grid", "grid": "TX_CENTRAL",
+                "md_unit": "ft", "source": "test survey"}
+        st_pv, _ = call(V, "POST", f"/api/wells/{wb['id']}/trajectory", body)
+        st_b64, _ = call(O, "POST", f"/api/wells/{wb['id']}/trajectory", {**body, "content_b64": "%%%"})
+        st_at, at = call(O, "POST", f"/api/wells/{wb['id']}/trajectory", body)
+        st_w, wsc = call(V, "GET", "/api/3d/scene/workspace")
+    finally:
+        svc.stop()
+    kept = (Path(wsp, "wells", wb["id"], "trajectory", "survey.csv").read_bytes() == csv_ok.encode())
+    audit = [e["action"] for e in ws.audit_log(500)]
+    cli = _sp.run([sys.executable, "-m", "gea", "well3d", "--scene", "workspace", "--workspace", wsp], capture_output=True, text=True, cwd=str(pkg.parent))
+    page_src = (pkg / "web" / "app.html").read_text(encoding="utf-8")
+    help_md = (pkg / "help" / "well3d.md").read_text(encoding="utf-8")
+    w0 = (wsc.get("wells") or [{}])[0]
+    c3 = {
+        "attach": att_refused == 4 and st_pv == 403 and st_b64 == 400 and st_at == 200 and at["stations"] == 6 and at["grid"] == "TX_CENTRAL" and kept
+                  and at["sha256"] == _hl.sha256(csv_ok.encode()).hexdigest() and "well.trajectory" in audit,
+        "scene": st_l == 200 and [s["id"] for s in ls["scenes"]] == ["workspace", "volve"] and st_v == 200 and len(vs["wells"]) == 3 and st_n == 404
+                 and st_w == 200 and len(wsc["wells"]) == 1 and abs(w0["td_md"] - 1752.6) < 0.01 and "TX_CENTRAL grid north" in w0["azimuth"]["applied"]
+                 and w0["replay"]["kind"] == "drilling recorder" and "Gamma" in w0["channels"] and not w0["replay"]["t"][0].endswith("Z")
+                 and [n["id"] for n in wsc["not_placed"]] == [wc["id"]] and wsc["layers"][0]["lines"][0]["z"] == round(1500 * 0.3048, 1)
+                 and any(n["display"] == "Faults" and "no depth" in n["why"] for n in wsc["layers_not_drawn"]),
+        "tops_alone": sc0["wells"] == [] and [l["name"] for l in sc0["layers"]] == ["Top A"] and sc0["frame"]["area"] == "the placed wells"
+                      and {n["id"]: n["why"] for n in sc0["not_placed"]}[wb["id"]] == "no directional survey attached",
+        "cli": cli.returncode == 0 and "Well B: 6 stations" in cli.stdout and "contour lines at their elevations" in cli.stdout,
+        "page": all(x in page_src for x in ("VIEWS['3d']", "['#/3d', '3D']", "getContext('webgl'", "OES_element_index_uint", "w3Tube", "w3SurfMesh",
+                                             "/api/3d/scene/", "w3Play", "the recorder's clock", "ArrowLeft", "see-through", "/trajectory`", "w3Hit")),
+        "help": any(t["topic"] == "well3d" for t in H.topics()) and not H.check("well3d") and H.VIEW_TOPIC.get("3d") == "well3d"
+                and "a well with no survey" in help_md.split("will not call a measurement")[1],
+    }
+    failed3 = [k for k, v in c3.items() if not v]
+    ok(not failed3,
+       "CE3 a workspace well in 3D: a survey attached only with its azimuth reference stated (a grid named, a declination given) and its source, "
+       "kept byte for byte and audited; the path turned from TX_CENTRAL grid north by the convergence there, set below sea level from the "
+       "declared rig floor among the formation-top contours at their elevations (a fault trace listed as having no depth; the tops drawn alone "
+       "before any survey), the recorder's bit depth replayed on its own clock and its gamma along the path; a well without a position listed, not drawn; the API, the CLI and the 3D "
+       "tab (WebGL by hand, the replay, turn and zoom by pointer and keys); the help page"
+       + (f" [failed: {', '.join(failed3)}]" if failed3 else ""))
+
+
+
+def section_cf_economics(tmp: str) -> None:
+    """Section CF - economics: the shipped EIA spot prices checked against their provenance and against two days EIA
+    published, series refreshed and imported only when they read, a well's cost from declared rates done by hand and
+    left undivided when the rates leave time uncovered, the market while the well was drilled, the API, the CLI, the tab."""
+    import base64 as _b64
+    import gzip as _gz
+    import hashlib as _hl
+    import http.cookiejar
+    import subprocess as _sp
+    import urllib.request
+    import urllib.error
+    from datetime import datetime as _dt
+    from . import economics as EC
+    from . import edr as E
+    from . import helplib as H
+    from .workspace import Workspace
+    from .service import Service, Users
+    pkg = Path(__file__).parent
+    d = Path(tmp, "cf"); d.mkdir()
+
+    prov = json.loads((pkg / "reference" / "eia_prices.provenance.json").read_text(encoding="utf-8"))
+    hashes = all(_hl.sha256((pkg / "reference" / s["file"]).read_bytes()).hexdigest() == s["sha256"]
+                 and _hl.sha256(_gz.decompress((pkg / "reference" / s["file"]).read_bytes())).hexdigest() == s["sha256_uncompressed"]
+                 and len(s["source"]["commit"]) == 40 for s in prov["series"].values())
+    wd, wv, _ = EC.read_series(None, "wti")
+    hd, hv, _ = EC.read_series(None, "henry_hub")
+    refused = 0
+    for bad in ("Date,Price\n2020-01-01,1\n2020-01-01,2\n", "Date,Price\n2020-01-01,1\n2020-01-02,x\n", "Date,Price\n2020-13-01,1\n2020-01-02,2\n", "Date,Price\n"):
+        try:
+            EC._parse_series(bad)
+        except ValueError:
+            refused += 1
+    bd, bv, _ = EC._parse_series("Date,Price\n2020-01-03,3\n2020-01-01,1\n2020-01-02,\n")
+    ev = E.selftest_events(str(d))
+    wsp = str(d / "ws"); ws = Workspace.create(wsp, "Pad 3", actor="tester")
+    wb = ws.add_well_file(ev["path"], display="Well B", actor="tester"); wc = ws.add_well_file(ev["path"], display="Well C", actor="tester")
+    later = ("Date,Price\n" + "".join(f"{x},{y}\n" for x, y in zip(wd, wv)) + "2099-01-01,50.0\n").encode()
+    hh_short = "Date,Price\n2001-01-02,9.0\n2001-01-03,9.5\n".encode()
+    calls = {}
+
+    def fake(u):
+        calls[u] = True
+        if "oil-prices" in u:
+            return later
+        if "natural-gas" in u:
+            return hh_short
+        raise OSError("no network")
+    ref = EC.refresh_prices(ws, "tester", fetch=fake)
+    ref_status = {r["id"]: r["status"] for r in ref}
+    imp_refused = 0
+    for kw in ({"name": "", "unit": "rigs", "source": "s"}, {"name": "Rigs", "unit": "", "source": "s"}, {"name": "Rigs", "unit": "rigs", "source": ""}, {"name": "wti", "unit": "x", "source": "s"}):
+        try:
+            EC.import_series(ws, "tester", kw["name"], kw["unit"], kw["source"], b"Date,Count\n2026-06-26,300\n2026-07-03,301\n")
+        except ValueError:
+            imp_refused += 1
+    imp = EC.import_series(ws, "tester", "Permian rig count", "rigs", "test", b"Week,Rigs\n2026-06-26,300\n2026-07-03,301\n2026-07-10,299\n")
+    c1 = {
+        "shipped": hashes and prov["series"]["wti"]["eia_series"] == "RWTC" and prov["series"]["henry_hub"]["eia_series"] == "RNGWHHD"
+                   and "public domain" in prov["license"].lower() and wv[wd.index("2020-04-20")] == -36.98 and hv[hd.index("2021-02-17")] == 23.86
+                   and len(wd) == prov["series"]["wti"]["rows"] and wd == sorted(wd),
+        "parse": refused == 4 and bd == ["2020-01-01", "2020-01-03"] and bv == [1.0, 3.0],
+        "refresh": ref_status == {"wti": "REFRESHED", "henry_hub": "NOT REFRESHED"} and EC.read_series(ws, "wti")[0][-1] == "2099-01-01"
+                   and EC.read_series(ws, "henry_hub")[2]["kind"] == "shipped" and len(calls) == 2,
+        "import": imp_refused == 4 and imp["id"] == "permian_rig_count" and imp["rows"] == 3 and imp["date_column"] == "Week" and imp["value_column"] == "Rigs"
+                  and Path(wsp, "economics", "series", "permian_rig_count.csv").read_bytes().startswith(b"Week,Rigs"),
+    }
+    failed1 = [k for k, v in c1.items() if not v]
+    ok(not failed1,
+       "CF1 the price series: EIA's daily WTI (RWTC) and Henry Hub (RNGWHHD) spot prices shipped with their provenance and hashes, carrying the "
+       "-36.98 of 20 April 2020 and the 23.86 of 17 February 2021 as EIA published them; a series with a repeated or bad date or a bad number "
+       "refused, a blank day left out; a refreshed copy kept only when it reads and runs later than the one held; an imported series kept byte "
+       "for byte with its name, unit and source"
+       + (f" [failed: {', '.join(failed1)}]" if failed1 else ""))
+
+    ws.refresh_dashboard(actor="tester")
+    rec = json.loads(Path(ws.reports_dir, "wells", wb["id"], "drilling_record.json").read_text(encoding="utf-8"))
+    t0, t1 = _dt.fromisoformat(rec["start"]), _dt.fromisoformat(rec["end"])
+    days = (t1 - t0).total_seconds() / 86400.0
+    mid = "2026-07-02T12:00:00"
+    EC.set_rate(ws, wb["id"], "tester", 30000, "rig spread", "contract", rec["start"], mid)
+    part = EC.economics(ws)
+    pb = next(w for w in part["wells"] if w["well_id"] == wb["id"])
+    bad = 0
+    for args in ((30000, "rig", "c", "2026-07-02", "2026-07-05"), (30000, "", "c", None, None), (30000, "rig", "", None, None), (-5, "rig", "c", "2027-01-01", None),
+                 (30000, "rig", "c", "2027-01-02", "2027-01-01")):
+        try:
+            EC.set_rate(ws, wb["id"], "tester", *args)
+        except ValueError:
+            bad += 1
+    EC.set_rate(ws, wb["id"], "tester", 32000, "rig spread", "contract", mid, None)
+    EC.add_item(ws, wb["id"], "tester", 18000, "bit", "invoice", 1)
+    EC.set_afe(ws, wb["id"], "tester", 150000, "AFE 1")
+    full = EC.economics(ws)
+    fb = next(w for w in full["wells"] if w["well_id"] == wb["id"])
+    first = (_dt.fromisoformat(mid) - t0).total_seconds() / 86400.0
+    by_hand = round(30000 * first + 32000 * (days - first) + 18000, 2)
+    win = [v for x, v in zip(wd, wv) if rec["start"][:10] <= x <= rec["end"][:10]]
+    c2 = {
+        "partial": pb["cost"]["complete"] is False and pb["cost"]["per_ft"] is None and "cover" in pb["cost"]["why_incomplete"]
+                   and abs(pb["cost"]["usd"] - round(30000 * first, 2)) < 0.01,
+        "refused": bad == 5,
+        "by_hand": fb["cost"]["complete"] and abs(fb["cost"]["usd"] - by_hand) < 0.01 and abs(fb["cost"]["per_ft"] - round(by_hand / rec["feet_drilled"], 2)) < 0.01
+                   and abs(fb["cost"]["vs_afe_pct"] - round(100 * (by_hand - 150000) / 150000, 1)) < 0.05 and fb["days"] == round(days, 2),
+        "sections": all(s["complete"] for s in fb["sections"]) and fb["sections"][0]["cost_usd"] > 18000,
+        "market": fb["market"]["wti"]["n"] == len(win) and abs(fb["market"]["wti"]["mean"] - round(sum(win) / len(win), 3)) < 1e-9
+                  and full["chart"]["wells"][0]["from"] == rec["start"][:10] and "permian_rig_count" in full["chart"]["series"],
+        "correlation": "spearman_cost_per_ft_vs_wti" not in full["correlation"] and full["correlation"]["n_wells"] == 1
+                       and abs(EC._spearman([1, 2, 3, 4], [10, 20, 30, 40]) - 1) < 1e-12 and abs(EC._spearman([1, 2, 3, 4], [4, 3, 2, 1]) + 1) < 1e-12,
+        "audit": all(a in [x["action"] for x in ws.audit_log(500)] for a in ("economics.rate", "economics.item", "economics.afe", "economics.prices.refresh", "economics.series.import")),
+    }
+    failed2 = [k for k, v in c2.items() if not v]
+    ok(not failed2,
+       "CF2 a well's cost: two day rates and a bit over the drilling record's own time come to the dollar the hand sum gives, with cost per foot "
+       "and the difference from the AFE; while the rates leave time uncovered the cost is shown, marked incomplete and not divided by the footage; "
+       "overlapping periods, a rate with no source or no purpose, a negative rate and a period that ends before it starts refused; the WTI mean "
+       "while drilling equal to the days EIA published in that window; no correlation across fewer than ten wells; every declaration audited"
+       + (f" [failed: {', '.join(failed2)}]" if failed2 else ""))
+
+    Users(str(Path(wsp, "users.json"))).add("v", "viewer-pass-1", "viewer"); Users(str(Path(wsp, "users.json"))).add("op", "op-pass-123", "operator")
+    svc = Service(wsp, port=0, scheduler=False).start()
+    base_url = svc.url.rstrip("/")
+
+    def opener(name, pw):
+        o = urllib.request.build_opener(urllib.request.HTTPCookieProcessor(http.cookiejar.CookieJar()))
+        rq = urllib.request.Request(base_url + "/api/login", method="POST", data=json.dumps({"name": name, "password": pw}).encode())
+        rq.add_header("Content-Type", "application/json"); rq.add_header("X-GEA-Action", "1"); o.open(rq, timeout=30)
+        return o
+
+    def call(o, method, path, body=None):
+        rq = urllib.request.Request(base_url + path, method=method, data=(json.dumps(body).encode() if body is not None else None))
+        if method == "POST":
+            rq.add_header("Content-Type", "application/json"); rq.add_header("X-GEA-Action", "1")
+        try:
+            with o.open(rq, timeout=120) as rr_:
+                return rr_.status, json.loads(rr_.read())
+        except urllib.error.HTTPError as e:
+            return e.code, json.loads(e.read())
+    try:
+        V, O = opener("v", "viewer-pass-1"), opener("op", "op-pass-123")
+        st_g, eg = call(V, "GET", "/api/economics")
+        st_vr, _ = call(V, "POST", "/api/economics/rate", {"well_id": wc["id"], "usd_per_day": 1000, "what": "rig", "source": "s"})
+        st_or, _ = call(O, "POST", "/api/economics/rate", {"well_id": wc["id"], "usd_per_day": 25000, "what": "rig spread", "source": "contract"})
+        st_bad, _ = call(O, "POST", "/api/economics/rate", {"well_id": wc["id"], "usd_per_day": "x", "what": "rig", "source": "s"})
+        st_imp, imp2 = call(O, "POST", "/api/economics/import", {"name": "Midland WTI", "unit": "USD/bbl", "source": "test", "content_b64": _b64.b64encode(b"Date,Price\n2026-07-01,70\n2026-07-02,71\n").decode()})
+        st_rm, _ = call(O, "POST", "/api/economics/remove", {"well_id": wc["id"], "kind": "rate", "index": 0})
+        st_rm2, _ = call(O, "POST", "/api/economics/remove", {"well_id": wc["id"], "kind": "rate", "index": 0})
+        st_ac, _ = call(O, "POST", "/api/economics/nope", {})
+    finally:
+        svc.stop()
+    cli = _sp.run([sys.executable, "-m", "gea", "economics", "--workspace", wsp], capture_output=True, text=True, cwd=str(pkg.parent))
+    page_src = (pkg / "web" / "app.html").read_text(encoding="utf-8")
+    help_md = (pkg / "help" / "economics.md").read_text(encoding="utf-8")
+    c3 = {
+        "api": st_g == 200 and len(eg["wells"]) == 2 and st_vr == 403 and st_or == 200 and st_bad == 400 and st_imp == 200 and imp2["id"] == "midland_wti"
+               and st_rm == 200 and st_rm2 == 404 and st_ac == 404,
+        "cli": cli.returncode == 0 and "USD/ft" in cli.stdout and "incomplete" not in cli.stdout.split("Well B")[1].split("Well C")[0] and "U.S. Energy Information Administration" in cli.stdout,
+        "page": all(x in page_src for x in ("VIEWS.economics", "['#/economics', 'Economics']", "ecChart", "ec-hit", "ArrowLeft", "/api/economics/rate",
+                                             "/api/economics/import", "/api/economics/refresh", "incomplete", "a gap of more than six days is not bridged")),
+        "help": any(t["topic"] == "economics" for t in H.topics()) and not H.check("economics") and H.VIEW_TOPIC.get("economics") == "economics"
+                and "fewer than ten wells" in help_md.split("will not call a measurement")[1],
+    }
+    failed3 = [k for k, v in c3.items() if not v]
+    ok(not failed3,
+       "CF3 the API (a viewer reads, only an operator declares; a bad number, an unknown action and a second removal of the same rate named), "
+       "the CLI and the Economics tab: one chart per series on its own axis with the wells' drilling windows shaded and a crosshair by pointer or "
+       "arrows, the wells' cost table, the cost-basis forms, refresh and import; the help page"
+       + (f" [failed: {', '.join(failed3)}]" if failed3 else ""))
+
+
 def main() -> int:
     print("GEA-Program - ACCEPTANCE SUITE (the product gate)")
     # a handle a job subprocess still holds at the end must not turn a finished gate into a traceback on Windows:
@@ -8698,6 +9055,8 @@ def main() -> int:
         section_cb_charts(tmp)
         section_cc_site_map(tmp)
         section_cd_control_room(tmp)
+        section_ce_well3d(tmp)
+        section_cf_economics(tmp)
     if _FAILS:
         print(f"[ACCEPTANCE] {len(_FAILS)} FAILURES ({_PASS} passed):")
         for f in _FAILS: