gea-program 0.8.1__tar.gz → 0.10.0__tar.gz

{gea_program-0.8.1 → gea_program-0.10.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gea-program
-Version: 0.8.1
+Version: 0.10.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
@@ -125,7 +125,7 @@ gea dashboard --catalog-well volve_f12_f14_production_excerpt:15/9-F-12:10000 --
 gea client-report --report accuracy --out client_report
 gea model-cards --out model_cards
 gea sbom --out sbom
-gea accept                     # the product gate (302 checks)
+gea accept                     # the product gate (314 checks)
 gea help drift                 # the help library, by the job (16 pages; the same text is on every dashboard page)
 gea guide                      # the click-by-click tester guide (docs/TESTER_GUIDE.md)
 gea gui                        # the desktop window (pip install "gea-program[desktop]")
@@ -528,6 +528,29 @@ so one geophone can still say when each rig was working - and the Seismic
 Station Report and the station page carry the section, with the tracks drawn
 against time.
 
+### The Seismicity Response Area packet
+
+The Railroad Commission of Texas declares a Seismicity Response Area around a
+cluster of felt earthquakes and expects the operators inside it to run a
+response plan. The plans on file have one shape - a circle of about 9 km, a
+plan written against M 3.5, an 18-month goal, a response group meeting
+within 48 hours of a threshold event, quarterly checkpoints with Commission
+staff, disposal wells in two tiers by the base of the Wolfcamp - and the
+Commission's December 2023 Notice to Operators names the data: four daily
+parameters (maximum and average surface injection pressure, injection
+volume, maximum injection rate) and three bottomhole-pressure methods, one
+of which is a permanent downhole probe. `gea/sra.py` builds the packet to
+that shape from what the site holds. The operator declares what the program
+cannot measure (`set-well`: surface position and datum, API and UIC, depth
+tier, which channels are the pressure, the rate and the downhole gauge, the
+BHP method); `sra-define` puts the area on the site; `sra-report --catalog
+texnet.csv` writes the packet - membership by geodesic distance, the daily
+record per well rolled up by month, the catalogue's events against the plan
+with the trigger and the goal clock, the schedule, and every gap by name -
+and a daily export under the Notice's own parameter names. It will not
+invent a volume for a well with no rate channel, a tier from depth alone, or
+an aftershock: those are declarations, and the packet says so.
+
 ### Audit / Update
 
 The Audit / Update tab holds three things. The whole audit log

{gea_program-0.8.1 → gea_program-0.10.0}/README.md

@@ -86,7 +86,7 @@ gea dashboard --catalog-well volve_f12_f14_production_excerpt:15/9-F-12:10000 --
 gea client-report --report accuracy --out client_report
 gea model-cards --out model_cards
 gea sbom --out sbom
-gea accept                     # the product gate (302 checks)
+gea accept                     # the product gate (314 checks)
 gea help drift                 # the help library, by the job (16 pages; the same text is on every dashboard page)
 gea guide                      # the click-by-click tester guide (docs/TESTER_GUIDE.md)
 gea gui                        # the desktop window (pip install "gea-program[desktop]")
@@ -489,6 +489,29 @@ so one geophone can still say when each rig was working - and the Seismic
 Station Report and the station page carry the section, with the tracks drawn
 against time.
 
+### The Seismicity Response Area packet
+
+The Railroad Commission of Texas declares a Seismicity Response Area around a
+cluster of felt earthquakes and expects the operators inside it to run a
+response plan. The plans on file have one shape - a circle of about 9 km, a
+plan written against M 3.5, an 18-month goal, a response group meeting
+within 48 hours of a threshold event, quarterly checkpoints with Commission
+staff, disposal wells in two tiers by the base of the Wolfcamp - and the
+Commission's December 2023 Notice to Operators names the data: four daily
+parameters (maximum and average surface injection pressure, injection
+volume, maximum injection rate) and three bottomhole-pressure methods, one
+of which is a permanent downhole probe. `gea/sra.py` builds the packet to
+that shape from what the site holds. The operator declares what the program
+cannot measure (`set-well`: surface position and datum, API and UIC, depth
+tier, which channels are the pressure, the rate and the downhole gauge, the
+BHP method); `sra-define` puts the area on the site; `sra-report --catalog
+texnet.csv` writes the packet - membership by geodesic distance, the daily
+record per well rolled up by month, the catalogue's events against the plan
+with the trigger and the goal clock, the schedule, and every gap by name -
+and a daily export under the Notice's own parameter names. It will not
+invent a volume for a well with no rate channel, a tier from depth alone, or
+an aftershock: those are declarations, and the packet says so.
+
 ### Audit / Update
 
 The Audit / Update tab holds three things. The whole audit log

{gea_program-0.8.1 → gea_program-0.10.0}/gea/__init__.py

@@ -122,7 +122,7 @@ from .profile_catalog import (
     read_pangaea_txt,
 )
 
-__version__ = "0.8.1"
+__version__ = "0.10.0"
 __all__ = [
     "aging_rate", "rate_psi_yr", "accuracy_psi",
     "Sensor", "SimulatorConfig", "DownholeEngine",

{gea_program-0.8.1 → gea_program-0.10.0}/gea/__main__.py

@@ -281,9 +281,36 @@ def main(argv=None) -> int:
     p_ws.add_argument("--action", type=str, default="list",
                       choices=["init", "add-file", "add-catalog", "add-live", "remove", "list", "migrate", "refresh", "audit",
                                "add-seismic", "refresh-seismic", "remove-seismic", "sar-film", "refresh-all", "add-track", "refresh-track", "remove-track",
-                               "seismic-dataset", "seismic-field", "add-site", "sites", "site-report", "remove-site", "rename"])
+                               "seismic-dataset", "seismic-field", "add-site", "sites", "site-report", "remove-site", "rename",
+                               "set-well", "sra-define", "sra-report", "associate"])
     p_ws.add_argument("--stations", type=str, nargs="+", default=None, help="add-track: two or more array station ids of this workspace")
     p_ws.add_argument("--site", type=str, default=None, help="site-report/remove-site: the site id")
+    p_ws.add_argument("--api", type=str, default=None, help="set-well: the API number")
+    p_ws.add_argument("--uic", type=str, default=None, help="set-well: the UIC permit number")
+    p_ws.add_argument("--depth-tier", type=str, choices=["shallow", "deep"], default=None, help="set-well: above or below the named formation base")
+    p_ws.add_argument("--formation", type=str, default=None, help="set-well: the basis of the tier, e.g. 'completed below the base of the Wolfcamp'")
+    p_ws.add_argument("--channel-pressure", type=str, default=None, help="set-well: the channel that is the surface injection pressure (psi)")
+    p_ws.add_argument("--channel-rate", type=str, default=None, help="set-well: the channel that is the injection rate")
+    p_ws.add_argument("--rate-unit", type=str, default=None, help="set-well: the rate channel's unit: bbl/min or bbl/d")
+    p_ws.add_argument("--bhp-method", type=str, choices=["calculated", "dip_in", "probe"], default=None, help="set-well: how bottomhole pressure is known")
+    p_ws.add_argument("--channel-bhp", type=str, default=None, help="set-well: the downhole gauge channel (probe method)")
+    p_ws.add_argument("--top-interval-ft", type=float, default=None, help="set-well: top of the disposal interval, ft (calculated method)")
+    p_ws.add_argument("--gradient-psi-ft", type=float, default=None, help="set-well: fluid gradient, psi/ft (calculated method)")
+    p_ws.add_argument("--radius-km", type=float, default=None, help="sra-define: the area's radius in km (Gardendale: 9.08)")
+    p_ws.add_argument("--plan-date", type=str, default=None, help="sra-define: the plan date, YYYY-MM-DD")
+    p_ws.add_argument("--threshold-m", type=float, default=3.5, help="sra-define: the magnitude the plan is written against (3.5)")
+    p_ws.add_argument("--goal-months", type=int, default=18, help="sra-define: months without a threshold event that meet the goal (18)")
+    p_ws.add_argument("--response-hours", type=int, default=48, help="sra-define: hours within which the response group meets (48)")
+    p_ws.add_argument("--checkpoint-months", type=int, default=3, help="sra-define: months between checkpoints with Commission staff (3)")
+    p_ws.add_argument("--formation-boundary", type=str, default="base of the Wolfcamp", help="sra-define: the stratigraphic boundary between the two tiers")
+    p_ws.add_argument("--triggering-event", type=str, default="", help="sra-define: the catalogue id of the event the plan answers")
+    p_ws.add_argument("--catalog", type=str, default=None, help="sra-report: a TexNet earthquake-catalogue export (CSV)")
+    p_ws.add_argument("--aftershocks", type=str, nargs="*", default=None, help="sra-report: catalogue ids the operator declares as aftershocks (exempt)")
+    p_ws.add_argument("--start", type=str, default=None, help="sra-report: first day of the record, YYYY-MM-DD")
+    p_ws.add_argument("--end", type=str, default=None, help="sra-report: last day of the record, YYYY-MM-DD")
+    p_ws.add_argument("--vp", type=float, default=5.8, help="associate: the P velocity of the model, km/s (declared)")
+    p_ws.add_argument("--depth-km", type=float, default=6.0, help="associate: the hypocentral depth the location is made under, km (declared)")
+    p_ws.add_argument("--rms-tol", type=float, default=0.15, help="associate: the RMS residual an event must fit within, s")
     p_ws.add_argument("--wells", type=str, nargs="+", default=None, help="add-site: well ids of this workspace")
     p_ws.add_argument("--seismic", type=str, nargs="+", default=None, help="add-site: seismic station ids of this workspace")
     p_ws.add_argument("--tracks", type=str, nargs="+", default=None, help="add-site: track ids of this workspace")
@@ -362,6 +389,9 @@ def main(argv=None) -> int:
     p_w0.add_argument("--seed", type=int, default=1)
     p_w0.add_argument("--connect", type=str, default=None, help="host:port - instead of listening, connect to a listening tap and push frames")
 
+    p_sra = sub.add_parser("sra", help="the Seismicity Response Area packet: the self-test on a labelled scene")
+    p_sra.add_argument("--action", choices=["selftest"], default="selftest")
+    p_sra.add_argument("--json", action="store_true")
     p_dr = sub.add_parser("doctor", help="which code is running and can it serve: Python, the package, duplicates, the page, PyPI; with --workspace also the site folder and the port")
     p_dr.add_argument("--workspace", type=str, default=None)
     p_dr.add_argument("--host", type=str, default="127.0.0.1")
@@ -395,7 +425,8 @@ def main(argv=None) -> int:
     p_se.add_argument("--action", choices=["info", "spectrum", "lines", "fetch", "stations", "detect", "selftest", "convert", "response", "remove-response",
                                            "beam", "array-detect", "locate", "array-selftest", "sar-film", "bearings", "track", "track-selftest",
                                            "signatures", "multi-beam", "multi-bearings", "multi-track", "signature-selftest", "field-selftest",
-                                           "harmonics", "harmonic-selftest", "datum", "geodesy-selftest", "array-qc", "qc-selftest", "unlisted", "unlisted-selftest", "bearing-sigma", "coverage", "uncertainty-selftest"], default="info")
+                                           "harmonics", "harmonic-selftest", "datum", "geodesy-selftest", "array-qc", "qc-selftest", "unlisted", "unlisted-selftest", "bearing-sigma", "coverage", "uncertainty-selftest",
+                                           "assoc-selftest"], default="info")
     p_se.add_argument("--signatures", type=str, default=None, help="multi-beam/multi-bearings: the signatures JSON from --action signatures")
     p_se.add_argument("--histories", type=str, nargs="+", default=None, help="track: two or more bearing-history JSONs (from --action bearings), one per array")
     p_se.add_argument("--truth", type=str, default=None, help="track: ground truth CSV (point_id, lat, lon, utc[, note]) - the lateral's points in time")
@@ -851,6 +882,18 @@ def main(argv=None) -> int:
                           + (f"; verdict {v['verdict']} ({v['hit']} of {v['of']} inside their ellipse)" if v else "")
                           + (f" - {row.get('detail', '')}" if row.get("detail") else ""))
             return 0 if r["ok"] else 1
+        if a.action == "assoc-selftest":
+            from . import seismic_assoc as _SA
+            r = _SA.selftest()
+            print(f"association selftest [{r['label']}]: {r['status']} - {r['n_picks']} picks, {len(r['events'])} events, errors {r['errors_km']} km")
+            for k, v in r["checks"].items():
+                print(f"  {'ok ' if v else 'BAD'} {k}")
+            for e in r["events"]:
+                print(f"  {e['id']}: {e['origin_utc']}  {e['lat']:.4f}, {e['lon']:.4f}  {e['n_stations']} stations  RMS {e['rms_s']} s  "
+                      f"misfit region {e['region_km']['extent_east_km']} x {e['region_km']['extent_north_km']} km  magnitude: not estimated")
+            for x in r["catalogue"]["rows"]:
+                print(f"  vs catalogue: {x['event']} {x['standing']}")
+            return 0 if r["status"] == "OK" else 1
         if a.action == "uncertainty-selftest":
             from . import seismic_uncertainty as UQ
             r = UQ.selftest()
@@ -1443,6 +1486,18 @@ def main(argv=None) -> int:
             return 0
         print(f"notifications: {'configured' if n.cfg else 'not configured'}" + (f" - {len(n.cfg['channels'])} channel(s), {len(n.cfg['rules'])} rule(s), quiet {n.cfg['quiet_s']} s" if n.cfg else " (gea notify --example)"))
         return 0
+    elif a.cmd == "sra":
+        from . import sra as _SRA
+        r = _SRA.selftest()
+        if a.json:
+            print(json.dumps({k: v for k, v in r.items() if k != "packet"}, indent=1))
+        else:
+            print(f"sra selftest [{r['label']}]: {r['status']}")
+            for k, v in r["checks"].items():
+                print(f"  {'ok ' if v else 'BAD'} {k}")
+            print(f"  deep well volume {r['deep_mean_volume_bbl_day']} bbl/day; probe offset {r['probe_offset_psi']} psi")
+            print(_SRA.report_text(r["packet"]))
+        return 0 if r["status"] == "OK" else 1
     elif a.cmd == "doctor":
         from .doctor import run as _doctor
         return _doctor(a.workspace, a.host, a.port, a.json)
@@ -1630,6 +1685,44 @@ def main(argv=None) -> int:
                 if not a.site:
                     raise SystemExit("remove-site needs --site <id>")
                 ws.remove_site(a.site, actor=a.actor); print("removed (folder kept):", a.site); return 0
+            elif a.action == "set-well":
+                if not a.well:
+                    raise SystemExit("set-well needs --well <id>")
+                surface = {"lat": a.lat, "lon": a.lon, "datum": a.datum} if (a.lat is not None or a.lon is not None) else None
+                disposal = {k: v for k, v in {"api_number": a.api, "uic_number": a.uic, "depth_tier": a.depth_tier, "formation_basis": a.formation,
+                                              "channel_pressure": a.channel_pressure, "channel_rate": a.channel_rate, "rate_unit": a.rate_unit,
+                                              "bhp_method": a.bhp_method, "channel_bhp": a.channel_bhp, "top_interval_ft": a.top_interval_ft,
+                                              "gradient_psi_ft": a.gradient_psi_ft}.items() if v is not None}
+                if not surface and not disposal:
+                    raise SystemExit("set-well: nothing to declare (give --lat/--lon, or any of --api --uic --depth-tier --channel-pressure ...)")
+                w = ws.set_well(a.well, surface=surface, disposal=disposal or None, actor=a.actor)
+                print(f"well {a.well}: surface {w.get('surface')}; disposal {w.get('disposal')}"); return 0
+            elif a.action == "sra-define":
+                if not a.site or not a.name or a.lat is None or a.lon is None or a.radius_km is None or not a.plan_date:
+                    raise SystemExit("sra-define needs --site, --name, --lat, --lon, --radius-km and --plan-date")
+                sra = ws.sra_define(a.site, actor=a.actor, name=a.name, lat=a.lat, lon=a.lon, radius_km=a.radius_km, plan_date=a.plan_date,
+                                    datum=a.datum or "WGS84", threshold_m=a.threshold_m, goal_months=a.goal_months, response_hours=a.response_hours,
+                                    checkpoint_months=a.checkpoint_months, formation_boundary=a.formation_boundary,
+                                    triggering_event=a.triggering_event, note=a.note or "")
+                print(f"sra: {sra['name']} on site {a.site} - {sra['radius_km']} km around {sra['centre']['lat']}, {sra['centre']['lon']} "
+                      f"({sra['centre']['datum']}); plan {sra['plan_date']}; M {sra['threshold_m']}; goal {sra['goal_months']} months"); return 0
+            elif a.action == "sra-report":
+                from .sra import report_text as _sra_text
+                if not a.site:
+                    raise SystemExit("sra-report needs --site <id>")
+                r = ws.write_sra_packet(a.site, catalog_csv=a.catalog, aftershocks=a.aftershocks, start=a.start, end=a.end, actor=a.actor)
+                print(_sra_text(r["packet"]))
+                print("  report:", r["paths"]["html"]); print("  daily export:", r["paths"]["daily_csv"], f"({r['export_rows']} rows)"); return 0
+            elif a.action == "associate":
+                from .seismic_assoc import report_text as _assoc_text
+                if not a.site:
+                    raise SystemExit("associate needs --site <id>")
+                r = ws.associate_site(a.site, vp_km_s=a.vp, depth_km=a.depth_km, band=a.band, start=a.start, end=a.end, catalog_csv=a.catalog,
+                                      rms_tol_s=a.rms_tol, actor=a.actor)
+                for sid, x in r["stations"].items():
+                    print(f"  {sid}: {x['status']}, {x['picks']} pick(s)")
+                print(_assoc_text(r["association"], r["catalogue"]))
+                print("  written:", os.path.join(ws.path, "reports", "sites", a.site, "association.json")); return 0
             elif a.action == "seismic-dataset":
                 r = ws.write_seismic_dataset_report(actor=a.actor)
                 inv = r["inventory"]

{gea_program-0.8.1 → gea_program-0.10.0}/gea/acceptance_tests.py

@@ -2914,7 +2914,8 @@ def section_am_report_samples(tmp: str) -> None:
             "model_card_well_baseline.html", "model_card_well_test_detector.html", "model_card_rock_density_inventory.html",
             "model_card_strata_property_estimator.html", "gauge_drift_report_monitored_SYNTHETIC.html", "alarm_event_report_SYNTHETIC.html",
             "data_resilience_report_SYNTHETIC.html", "sla_report_SYNTHETIC.html", "sat_protocol.html", "seismic_station_report_SYNTHETIC.html", "seismic_track_report_SYNTHETIC.html",
-            "seismic_dataset_report_SYNTHETIC.html", "seismic_field_report_SYNTHETIC.html", "site_report_SYNTHETIC.html"}
+            "seismic_dataset_report_SYNTHETIC.html", "seismic_field_report_SYNTHETIC.html", "site_report_SYNTHETIC.html",
+            "sra_packet_SYNTHETIC.html"}
     have = {p.name for p in d.glob("*.html")}
     md = (d / "SAMPLES.md").read_text(encoding="utf-8") if (d / "SAMPLES.md").exists() else ""
     texts = {n: (d / n).read_text(encoding="utf-8", errors="replace") for n in sorted(have & want)}
@@ -4593,9 +4594,248 @@ def section_bb_supervisor(tmp: str) -> None:
        "panel's own page carries the Stop and Restart controls, the authorisation step behind a popup, and the auto-restart switch")
 
 
+def section_bc_sra(tmp: str) -> None:
+    """Section BC - the Seismicity Response Area packet. The Railroad Commission's operator-led response plans
+    and its Notice to Operators on disposal-well monitoring name a shape; the packet is built to it from what
+    the site holds, and names what it is missing rather than leaving a column blank."""
+    import csv as _csv
+    import subprocess as _sp
+    from datetime import datetime as _dt, timedelta as _td, timezone as _tz
+    from . import sra as SR
+    from .client_reports import forbidden_terms, sra_packet_report, write
+    from .workspace import Workspace, WorkspaceError
+    pkg = Path(__file__).parent
+    d = Path(tmp, "bc"); d.mkdir()
+
+    # -- the self-test: the labelled scene ----------------------------------------------
+    r = SR.selftest()
+    ok(r["status"] == "OK" and all(r["checks"].values()) and r["label"] == "SIMULATION_SELF_TEST"
+       and abs(r["deep_mean_volume_bbl_day"] - 8640.0) / 8640.0 < 0.03 and abs(r["probe_offset_psi"] - 0.465 * 9800) < 15,
+       "BC1 the self-test on a labelled scene: membership by geodesic distance with a NAD27 well's shift named, the four daily parameters from "
+       "the well's own channels with the volume integrated from the rate to within 3 %, two shut-in days at zero, the probe's offset from "
+       "surface pressure recovered, the shallow well refusing a volume it has no rate channel for, the catalogue's trigger and exempt "
+       "aftershock and goal clock, the gaps named, the checkpoints from the plan date")
+
+    # -- the area: declarations, refused when not an area --------------------------------
+    sra = SR.define("Test SRA", 31.95, -102.25, 9.08, "2026-02-25", "NAD27", triggering_event="x")
+    bad = 0
+    for kw in ({"datum": "UNKNOWN"}, {"radius_km": 0}, {"plan_date": "not a date"}, {"name": "  "}):
+        try:
+            SR.define(**{**{"name": "T", "lat": 31.95, "lon": -102.25, "radius_km": 9.08, "plan_date": "2026-02-25", "datum": "WGS84"}, **kw}); 
+        except ValueError:
+            bad += 1
+    cps = SR.checkpoints(sra)
+    ok(sra["centre"]["datum"] == "NAD27" and sra["centre_wgs84"]["shift_m"] > 10 and sra["area_km2"] == 259.0 and sra["threshold_m"] == 3.5
+       and sra["goal_months"] == 18 and sra["response_hours"] == 48 and sra["checkpoint_months"] == 3 and bad == 4
+       and [c["date"] for c in cps][:2] == ["2026-05-25", "2026-08-25"] and cps[-1]["n"] == "goal" and cps[-1]["date"] == "2027-08-25"
+       and "not a measurement" in " ".join(sra["not_a_measurement"]).lower().replace("never derived here", "not a measurement"),
+       "BC2 an area is a centre on a stated datum (carried to WGS84 with its shift named), a radius, a plan date and the plan's numbers - "
+       "M 3.5, 18 months, 48 hours, quarterly - recorded as declarations; an unknown datum, a zero radius, a non-date or a blank name is refused; "
+       "the checkpoints fall every three months from the plan date and the goal date closes the list")
+
+    # -- the catalogue reader: a TexNet-shaped export ------------------------------------
+    cat = d / "texnet.csv"
+    with open(cat, "w", newline="") as f:
+        w = _csv.writer(f); w.writerow(["EventID", "Origin Date", "Origin Time", "Local Magnitude", "Latitude", "Longitude", "Depth of Hypocenter (Km)"])
+        w.writerow(["texnet2026aaaa", "2026-02-20", "04:15:10", "4.4", "31.950", "-102.250", "7.2"])      # before the plan date
+        w.writerow(["texnet2026aaab", "2026-03-05", "11:02:33", "3.7", "31.960", "-102.240", "6.8"])      # counts
+        w.writerow(["texnet2026aaac", "2026-03-06", "01:00:00", "3.6", "31.955", "-102.245", "6.5"])      # declared aftershock
+        w.writerow(["texnet2026aaad", "2026-03-09", "21:40:00", "2.0", "31.940", "-102.260", "6.1"])
+        w.writerow(["texnet2026aaae", "2026-03-12", "03:00:00", "3.9", "32.500", "-102.700", "8.0"])      # outside
+        w.writerow(["", "", "", "", "", "", ""])                                                           # a blank line is skipped, not an event
+    c = SR.read_catalog(str(cat))
+    se = SR.seismicity(sra, c, as_of=_dt(2026, 4, 1, tzinfo=_tz.utc), aftershocks=["texnet2026aaac"])
+    ok(c["n_events"] == 5 and c["n_skipped"] == 1 and c["columns"]["origin_date"] == "Origin Date" and c["columns"]["origin_time"] == "Origin Time"
+       and c["events"][0]["time"] == "2026-02-20T04:15:10Z" and c["events"][0]["depth_km"] == 7.2
+       and se["n_inside"] == 4 and se["n_at_or_above"] == 3 and se["n_counted"] == 1 and se["n_exempt"] == 1
+       and se["largest"]["event_id"] == "texnet2026aaaa" and se["trigger"]["status"] == "TRIGGERED" and se["trigger"]["event_id"] == "texnet2026aaab"
+       and se["trigger"]["respond_by"] == "2026-03-07T11:02:33Z" and se["goal"]["clock_started"] == "2026-03-05" and se["goal"]["status"] == "RUNNING"
+       and se["goal"]["goal_date"] == "2027-09-05",
+       "BC3 a TexNet-shaped export is read by its column names (date and time in separate columns, the depth under its long name), a blank line "
+       "is skipped and counted; against the plan: the M 4.4 before the plan date is the largest but does not count, the declared aftershock is "
+       "exempt, the M 3.9 outside the circle is outside, the M 3.7 triggers the 48-hour response with its deadline, and the 18-month clock "
+       "starts from it")
+
+    # -- the daily record from a stream: refused without a calendar --------------------------------
+    from .ports import LiveStream, StreamChannel
+    import numpy as _np
+    st = LiveStream(name="x", source_format="synthetic", index_kind="time_s", index=_np.arange(100) * 600.0,
+                    channels={"P": StreamChannel("P", "psi", _np.full(100, 1500.0))}, meta={})
+    rr = SR.daily_records(st, {"channel_pressure": "P"})
+    sc = SR.synthetic_scene(days=3)
+    dr = SR.daily_records(sc["streams"]["deep-1"], {**sc["wells"][0]["disposal"], "rate_unit": "furlongs"})
+    ok(rr["status"] == "REFUSED" and "calendar" in rr["detail"]
+       and dr["parameters"]["injection_volume_bbl"]["status"] == "NOT RECORDED" and "furlongs" in dr["parameters"]["injection_volume_bbl"]["why"]
+       and dr["parameters"]["max_surface_injection_pressure_psi"]["status"] == "RECORDED",
+       "BC4 a daily record needs a calendar: a stream with no start time is refused, not dated from zero; a rate in a unit the program does not "
+       "convert leaves the volume NOT RECORDED with the unit named, while the pressures are still recorded")
+
+    # -- through the workspace and the command line ----------------------------------------
+    t0 = _dt(2026, 3, 1, tzinfo=_tz.utc); n = 6 * 144
+    hist = d / "swd.csv"
+    with open(hist, "w", newline="") as f:
+        w = _csv.writer(f); w.writerow(["timestamp", "P_surf_psi", "Q_bpm", "P_bh_psi"])
+        for i in range(n):
+            w.writerow([(t0 + _td(seconds=600 * i)).strftime("%Y-%m-%dT%H:%M:%SZ"), "1500.0", "5.0", f"{1500.0 + 0.465 * 9000:.1f}"])
+    ws = Workspace.create(str(d / "ws"), "Pad 3", actor="tester")
+    w1 = ws.add_well_file(str(hist), display="SWD 1", actor="tester")
+    try:
+        ws.set_well(w1["id"], disposal={"depth_tier": "middle"}, actor="tester"); tier_bad = False
+    except WorkspaceError:
+        tier_bad = True
+    run = lambda *a: _sp.run([sys.executable, "-m", "gea", "workspace", "--path", ws.path, "--actor", "tester", *a], capture_output=True, text=True)
+    r1 = run("--action", "set-well", "--well", w1["id"], "--lat", "31.96", "--lon", "-102.24", "--datum", "NAD27", "--api", "42-329-12345", "--uic", "123456",
+             "--depth-tier", "deep", "--formation", "completed below the base of the Wolfcamp", "--channel-pressure", "P_surf_psi", "--channel-rate", "Q_bpm",
+             "--rate-unit", "bbl/min", "--bhp-method", "probe", "--channel-bhp", "P_bh_psi")
+    r2 = run("--action", "add-site", "--name", "SRA block", "--wells", w1["id"], "--client", "Acme Operating")
+    r3 = run("--action", "sra-report", "--site", "SRA-block")                                        # no area yet: refused
+    r4 = run("--action", "sra-define", "--site", "SRA-block", "--name", "Test SRA", "--lat", "31.95", "--lon", "-102.25", "--datum", "WGS84",
+             "--radius-km", "9.08", "--plan-date", "2026-02-25")
+    r5 = run("--action", "sra-report", "--site", "SRA-block", "--catalog", str(cat), "--aftershocks", "texnet2026aaac")
+    ws2 = Workspace(ws.path)
+    out = Path(ws.path, "reports", "sites", "SRA-block")
+    rows = list(_csv.DictReader(open(out / "sra_daily_export.csv", encoding="utf-8")))
+    md = (out / "sra_packet.md").read_text(encoding="utf-8")
+    audit = Path(ws.path, "records", "audit.jsonl").read_text(encoding="utf-8")
+    pkj = json.loads((out / "sra_packet.json").read_text(encoding="utf-8"))["packet"]
+    ok(tier_bad and r1.returncode == 0 and r2.returncode == 0 and r3.returncode != 0 and "sra-define" in (r3.stderr + r3.stdout)
+       and r4.returncode == 0 and r5.returncode == 0 and "COMPLETE" in r5.stdout and "TRIGGERED" in r5.stdout
+       and ws2.site("SRA-block")["sra"]["name"] == "Test SRA" and pkj["status"] == "COMPLETE" and pkj["wells"][0]["inside"] and pkj["wells"][0]["depth_tier"] == "deep"
+       and len(rows) == 6 and rows[0]["API number"] == "42-329-12345" and rows[0]["UIC permit number"] == "123456"
+       and abs(float(rows[1]["Injection volume (barrels per day)"]) - 7200.0) < 10 and rows[1]["Maximum surface injection pressure (pounds per square inch)"] == "1500.0"
+       and abs(float(rows[1]["Bottomhole pressure (pounds per square inch)"]) - (1500 + 0.465 * 9000)) < 1 and rows[1]["BHP method"] == "probe"
+       and "Maximum surface injection pressure (pounds per square inch)" in rows[0]
+       and not forbidden_terms(md) and "well.declare" in audit and "site.sra_define" in audit and "site.sra_packet" in audit,
+       "BC5 `set-well` records what the operator declares - surface position and datum, API and UIC, depth tier by the named formation, which "
+       "channels are the surface pressure, the rate and the downhole gauge, the BHP method - and refuses a tier that is not shallow or deep; "
+       "`sra-define` puts the area on the site; `sra-report` is refused until it does; then the packet is COMPLETE for a well with everything "
+       "declared, the daily export carries the Notice's parameter names with the API and UIC on every line, the volume integrates to 7,200 bbl "
+       "a day from 5 bbl/min, the probe is the bottomhole pressure, the report passes the vocabulary gate, and every step is audited")
+
+    # -- the service offers it and the page carries it ------------------------------------------
+    svc_src = (pkg / "service.py").read_text(encoding="utf-8"); page_src = (pkg / "web" / "app.html").read_text(encoding="utf-8")
+    r6 = _sp.run([sys.executable, "-m", "gea", "sra"], capture_output=True, text=True)
+    ok("'sra_report'" in svc_src and "sra_packet.html" in svc_src and "Seismicity Response Area</th>" in page_src and "sra-define" in page_src
+       and r6.returncode == 0 and "SIMULATION_SELF_TEST" in r6.stdout and "OK" in r6.stdout
+       and (pkg / "help" / "sra.md").is_file(),
+       "BC6 the Sites page carries the area and the packet's status beside each site; `gea sra` runs the self-test from the command line and "
+       "labels it SIMULATION_SELF_TEST; the help page exists")
+
+
+def section_bd_association(tmp: str) -> None:
+    """Section BD - association: several stations heard the same thing, or they did not. The onset picker,
+    the associator with its grid search and its refusals, the location with its misfit region, the array as
+    one station, the catalogue comparison, and the workspace run that feeds the SRA packet."""
+    import csv as _csv
+    import subprocess as _sp
+    import numpy as _np
+    from . import seismic as S
+    from . import seismic_assoc as SA
+    from .workspace import Workspace, WorkspaceError
+    d = Path(tmp, "bd"); d.mkdir()
+
+    # -- the self-test -----------------------------------------------------------------------
+    r = SA.selftest()
+    ok(r["status"] == "OK" and all(r["checks"].values()) and r["label"] == "SIMULATION_SELF_TEST"
+       and r["errors_km"]["ev-A"] < 1.0 and r["errors_km"]["ev-B"] < 1.0,
+       "BD1 the self-test on a labelled scene: picks at every station, two events located to within a kilometre and a third of a second with all "
+       "five stations and RMS under 0.05 s, a two-station event refused rather than located, a spurious burst left unused, finite misfit "
+       "regions, no magnitude, and against the catalogue one AGREES, one DIFFERS by its separation and one catalogue event the stations did "
+       "not associate")
+
+    # -- the picker: the onset is the first energy, not half a window early ------------------
+    sc = SA.synthetic_scene()
+    fr = sc["frame"]
+    from .seismic_array import local_xy
+    devs = []
+    for sid, s in sc["stations"].items():
+        ex, ey = local_xy(sc["truth"][0]["lat"], sc["truth"][0]["lon"], fr["lat0"], fr["lon0"])
+        tt = SA.travel_time(ex, ey, s, sc["vp"], sc["depth_km"])
+        pk = [p for p in SA.picks(sc["traces"][sid], band=(2.0, 20.0)) if abs(p.time - (sc["truth"][0]["t0"] + tt)) < 2.0]
+        devs.append(pk[0].time - sc["truth"][0]["t0"] - tt if pk else 9.0)
+    quiet = S.Trace("XX", "Q", "", "HHZ", 0.0, 100.0, _np.random.default_rng(1).normal(0, 1, 60000).astype(_np.float32), "synthetic", "FLOAT32")
+    ok(all(abs(x) < 0.08 for x in devs) and not SA.picks(quiet, band=(2.0, 20.0)) and all(p.snr > 10 for p in SA.picks(sc["traces"]["STA5"], band=(2.0, 20.0))),
+       "BD2 the picker puts every onset within 80 ms of the true arrival (the first refinement walked back into the noise and every pick came "
+       "out 150 ms early, which the location absorbed into the origin time and nobody would have seen), picks nothing on a minute of plain "
+       "noise, and carries a signal-to-noise ratio on every pick")
+
+    # -- the locator: refuses two stations, refines about the coarse best, names the model ----
+    two = SA.locate({"STA1": 10.0, "STA2": 10.5}, sc["stations"], sc["vp"], sc["depth_km"])
+    ex, ey = local_xy(sc["truth"][1]["lat"], sc["truth"][1]["lon"], fr["lat0"], fr["lon0"])
+    pset = {sid: sc["truth"][1]["t0"] + SA.travel_time(ex, ey, s, sc["vp"], sc["depth_km"]) for sid, s in sc["stations"].items()}
+    exact = SA.locate(pset, sc["stations"], sc["vp"], sc["depth_km"])
+    off = SA.locate({**pset, "STA2": pset["STA2"] + 0.6}, sc["stations"], sc["vp"], sc["depth_km"])
+    ok(two["status"] == "NOT ASSOCIATED" and "three" in two["detail"]
+       and exact["status"] == "LOCATED" and abs(exact["best"]["x_km"] - ex) < 0.1 and abs(exact["best"]["y_km"] - ey) < 0.1 and exact["best"]["rms_s"] < 0.01
+       and abs(exact["best"]["t0"] - sc["truth"][1]["t0"]) < 0.02 and not exact["region_km"]["reaches_grid_edge"]
+       and "flat earth" in exact["model"]["basis"] and exact["model"]["depth_km"] == sc["depth_km"]
+       and off["best"]["rms_s"] > 0.15 and abs(off["residuals_s"]["STA2"]) == max(abs(v) for v in off["residuals_s"].values()),
+       "BD3 the locator refuses two stations, finds exact picks to within 100 m of the true epicentre and 20 ms of the origin time (the first "
+       "version refined about the stations' centroid instead of the coarse best node and lost every event more than a kilometre from the "
+       "centre), names the straight-ray flat-earth model and the declared depth on the result, and a pick 0.6 s off shows as the largest residual")
+
+    # -- the associator: the outlier is dropped, the event stands on the rest ----------------
+    picks_all = []
+    for sid, tr in sc["traces"].items():
+        picks_all += SA.picks(tr, band=(2.0, 20.0))
+    bad = [SA.Pick("STA4", sc["truth"][0]["t0"] + 3.2, 50.0, 1.0)]       # a wrong pick at STA4 near event A
+    picks_mix = [p for p in picks_all if not (p.station == "STA4" and abs(p.time - sc["truth"][0]["t0"]) < 3.0)] + bad
+    res = SA.associate(picks_mix, sc["stations"], sc["vp"], sc["depth_km"], frame_info=fr)
+    evA = sorted(res["events"], key=lambda e: e["origin_utc"])[0]
+    ok(res["n_events"] == 2 and evA["n_stations"] == 4 and evA["dropped"] == ["STA4"] and evA["rms_s"] < 0.05
+       and evA["magnitude"] is None and "not estimated" in evA["magnitude_basis"],
+       "BD4 a wrong pick at one station is dropped as the largest residual and the event stands on the other four; an event carries no "
+       "magnitude and says why")
+
+    # -- an array is one station ----------------------------------------------------------------
+    raw = [SA.Pick("ARR", 100.0 + k * 0.03, 20.0 + k, 1.0, channel=f"S{k}") for k in range(6)] + [SA.Pick("ARR", 130.0, 9.0, 0.5, channel="S3")]
+    col = SA.collapse_array_picks(raw, "ARR")
+    ok(len(col) == 1 and abs(col[0].time - 100.075) < 0.02 and "6 sensor" in col[0].channel,
+       "BD5 an array's sensors' picks within half a second become one pick at the median time at the array's position, and a burst one sensor "
+       "alone saw is not an arrival at the array")
+
+    # -- through the workspace and the command line, into the packet ---------------------------
+    for sid, tr in sc["traces"].items():
+        S.write_mseed([tr], str(d / f"{sid}.mseed"))
+    cat = d / "cat.csv"
+    with open(cat, "w", newline="") as f:
+        w = _csv.writer(f); w.writerow(["EventID", "Origin Time", "Local Magnitude", "Latitude", "Longitude", "Depth"])
+        for c in sc["catalog"]:
+            w.writerow([c["event_id"], c["time"], c["magnitude"], c["lat"], c["lon"], 6.0])
+    ws = Workspace.create(str(d / "ws"), "Pad 3", actor="tester")
+    ids = [ws.add_seismic_station(sid, [str(d / f"{sid}.mseed")], st.lat, st.lon, "tester", datum=st.datum, band=[2.0, 20.0])["id"]
+           for sid, st in sc["stations"].items()]
+    ws.add_site("Assoc block", seismic=ids[:2], client="Acme", actor="tester")
+    try:
+        ws.associate_site("Assoc-block"); few = False
+    except WorkspaceError as e:
+        few = "three" in str(e)
+    ws.add_site("Assoc block 5", seismic=ids, client="Acme", actor="tester")
+    run = lambda *a: _sp.run([sys.executable, "-m", "gea", "workspace", "--path", ws.path, "--actor", "tester", *a], capture_output=True, text=True)
+    r1 = run("--action", "associate", "--site", "Assoc-block-5", "--catalog", str(cat))
+    r2 = run("--action", "sra-define", "--site", "Assoc-block-5", "--name", "Test SRA", "--lat", "31.95", "--lon", "-102.25", "--radius-km", "9.08", "--plan-date", "2026-02-25")
+    r3 = run("--action", "sra-report", "--site", "Assoc-block-5", "--catalog", str(cat))
+    aj = json.loads(Path(ws.path, "reports", "sites", "Assoc-block-5", "association.json").read_text(encoding="utf-8"))
+    md = Path(ws.path, "reports", "sites", "Assoc-block-5", "sra_packet.md").read_text(encoding="utf-8")
+    pk = json.loads(Path(ws.path, "reports", "sites", "Assoc-block-5", "sra_packet.json").read_text(encoding="utf-8"))["packet"]
+    audit = Path(ws.path, "records", "audit.jsonl").read_text(encoding="utf-8")
+    r4 = _sp.run([sys.executable, "-m", "gea", "seismic", "--action", "assoc-selftest"], capture_output=True, text=True)
+    ok(few and r1.returncode == 0 and "2 event(s)" in r1.stdout and "AGREES" in r1.stdout and "DIFFERS" in r1.stdout and "NOT ASSOCIATED BY THE STATIONS" in r1.stdout
+       and aj["association"]["n_events"] == 2 and all(v["status"] == "PICKED" for v in aj["stations"].values())
+       and r2.returncode == 0 and r3.returncode == 0 and pk["association"]["n_inside"] == 2
+       and "Events located by the site's own stations" in md and "AGREES" in md and "no magnitude is estimated here" in md
+       and "site.associate" in audit and r4.returncode == 0 and "SIMULATION_SELF_TEST" in r4.stdout,
+       "BD6 `--action associate` refuses a site with fewer than three stations, picks five miniSEED stations, associates and locates two events, "
+       "sets them against the catalogue, writes association.json and audits it; the SRA packet then carries the events the site's own stations "
+       "located beside the catalogue's, inside or outside the area, with the model named and no magnitude; `gea seismic --action assoc-selftest` "
+       "runs the labelled scene from the command line")
+
+
 def main() -> int:
     print("GEA-Program - ACCEPTANCE SUITE (the product gate)")
-    with tempfile.TemporaryDirectory() as tmp:
+    # a handle a job subprocess still holds at the end must not turn a finished gate into a traceback on Windows:
+    # the checks have already passed or failed on their own, and the folder is the temp folder
+    with tempfile.TemporaryDirectory(ignore_cleanup_errors=True) as tmp:
         section_a_cli(tmp)
         section_b_las(tmp)
         section_c_reconciler(tmp)
@@ -4648,6 +4888,8 @@ def main() -> int:
         section_az_uncertainty(tmp)
         section_ba_sites(tmp)
         section_bb_supervisor(tmp)
+        section_bc_sra(tmp)
+        section_bd_association(tmp)
     if _FAILS:
         print(f"[ACCEPTANCE] {len(_FAILS)} FAILURES ({_PASS} passed):")
         for f in _FAILS: