serve.py 175 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798991001011021031041051061071081091101111121131141151161171181191201211221231241251261271281291301311321331341351361371381391401411421431441451461471481491501511521531541551561571581591601611621631641651661671681691701711721731741751761771781791801811821831841851861871881891901911921931941951961971981992002012022032042052062072082092102112122132142152162172182192202212222232242252262272282292302312322332342352362372382392402412422432442452462472482492502512522532542552562572582592602612622632642652662672682692702712722732742752762772782792802812822832842852862872882892902912922932942952962972982993003013023033043053063073083093103113123133143153163173183193203213223233243253263273283293303313323333343353363373383393403413423433443453463473483493503513523533543553563573583593603613623633643653663673683693703713723733743753763773783793803813823833843853863873883893903913923933943953963973983994004014024034044054064074084094104114124134144154164174184194204214224234244254264274284294304314324334344354364374384394404414424434444454464474484494504514524534544554564574584594604614624634644654664674684694704714724734744754764774784794804814824834844854864874884894904914924934944954964974984995005015025035045055065075085095105115125135145155165175185195205215225235245255265275285295305315325335345355365375385395405415425435445455465475485495505515525535545555565575585595605615625635645655665675685695705715725735745755765775785795805815825835845855865875885895905915925935945955965975985996006016026036046056066076086096106116126136146156166176186196206216226236246256266276286296306316326336346356366376386396406416426436446456466476486496506516526536546556566576586596606616626636646656666676686696706716726736746756766776786796806816826836846856866876886896906916926936946956966976986997007017027037047057067077087097107117127137147157167177187197207217227237247257267277287297307317327337347357367377387397407417427437447457467477487497507517527537547557567577587597607617627637647657667677687697707717727737747757767777787797807817827837847857867877887897907917927937947957967977987998008018028038048058068078088098108118128138148158168178188198208218228238248258268278288298308318328338348358368378388398408418428438448458468478488498508518528538548558568578588598608618628638648658668678688698708718728738748758768778788798808818828838848858868878888898908918928938948958968978988999009019029039049059069079089099109119129139149159169179189199209219229239249259269279289299309319329339349359369379389399409419429439449459469479489499509519529539549559569579589599609619629639649659669679689699709719729739749759769779789799809819829839849859869879889899909919929939949959969979989991000100110021003100410051006100710081009101010111012101310141015101610171018101910201021102210231024102510261027102810291030103110321033103410351036103710381039104010411042104310441045104610471048104910501051105210531054105510561057105810591060106110621063106410651066106710681069107010711072107310741075107610771078107910801081108210831084108510861087108810891090109110921093109410951096109710981099110011011102110311041105110611071108110911101111111211131114111511161117111811191120112111221123112411251126112711281129113011311132113311341135113611371138113911401141114211431144114511461147114811491150115111521153115411551156115711581159116011611162116311641165116611671168116911701171117211731174117511761177117811791180118111821183118411851186118711881189119011911192119311941195119611971198119912001201120212031204120512061207120812091210121112121213121412151216121712181219122012211222122312241225122612271228122912301231123212331234123512361237123812391240124112421243124412451246124712481249125012511252125312541255125612571258125912601261126212631264126512661267126812691270127112721273127412751276127712781279128012811282128312841285128612871288128912901291129212931294129512961297129812991300130113021303130413051306130713081309131013111312131313141315131613171318131913201321132213231324132513261327132813291330133113321333133413351336133713381339134013411342134313441345134613471348134913501351135213531354135513561357135813591360136113621363136413651366136713681369137013711372137313741375137613771378137913801381138213831384138513861387138813891390139113921393139413951396139713981399140014011402140314041405140614071408140914101411141214131414141514161417141814191420142114221423142414251426142714281429143014311432143314341435143614371438143914401441144214431444144514461447144814491450145114521453145414551456145714581459146014611462146314641465146614671468146914701471147214731474147514761477147814791480148114821483148414851486148714881489149014911492149314941495149614971498149915001501150215031504150515061507150815091510151115121513151415151516151715181519152015211522152315241525152615271528152915301531153215331534153515361537153815391540154115421543154415451546154715481549155015511552155315541555155615571558155915601561156215631564156515661567156815691570157115721573157415751576157715781579158015811582158315841585158615871588158915901591159215931594159515961597159815991600160116021603160416051606160716081609161016111612161316141615161616171618161916201621162216231624162516261627162816291630163116321633163416351636163716381639164016411642164316441645164616471648164916501651165216531654165516561657165816591660166116621663166416651666166716681669167016711672167316741675167616771678167916801681168216831684168516861687168816891690169116921693169416951696169716981699170017011702170317041705170617071708170917101711171217131714171517161717171817191720172117221723172417251726172717281729173017311732173317341735173617371738173917401741174217431744174517461747174817491750175117521753175417551756175717581759176017611762176317641765176617671768176917701771177217731774177517761777177817791780178117821783178417851786178717881789179017911792179317941795179617971798179918001801180218031804180518061807180818091810181118121813181418151816181718181819182018211822182318241825182618271828182918301831183218331834183518361837183818391840184118421843184418451846184718481849185018511852185318541855185618571858185918601861186218631864186518661867186818691870187118721873187418751876187718781879188018811882188318841885188618871888188918901891189218931894189518961897189818991900190119021903190419051906190719081909191019111912191319141915191619171918191919201921192219231924192519261927192819291930193119321933193419351936193719381939194019411942194319441945194619471948194919501951195219531954195519561957195819591960196119621963196419651966196719681969197019711972197319741975197619771978197919801981198219831984198519861987198819891990199119921993199419951996199719981999200020012002200320042005200620072008200920102011201220132014201520162017201820192020202120222023202420252026202720282029203020312032203320342035203620372038203920402041204220432044204520462047204820492050205120522053205420552056205720582059206020612062206320642065206620672068206920702071207220732074207520762077207820792080208120822083208420852086208720882089209020912092209320942095209620972098209921002101210221032104210521062107210821092110211121122113211421152116211721182119212021212122212321242125212621272128212921302131213221332134213521362137213821392140214121422143214421452146214721482149215021512152215321542155215621572158215921602161216221632164216521662167216821692170217121722173217421752176217721782179218021812182218321842185218621872188218921902191219221932194219521962197219821992200220122022203220422052206220722082209221022112212221322142215221622172218221922202221222222232224222522262227222822292230223122322233223422352236223722382239224022412242224322442245224622472248224922502251225222532254225522562257225822592260226122622263226422652266226722682269227022712272227322742275227622772278227922802281228222832284228522862287228822892290229122922293229422952296229722982299230023012302230323042305230623072308230923102311231223132314231523162317231823192320232123222323232423252326232723282329233023312332233323342335233623372338233923402341234223432344234523462347234823492350235123522353235423552356235723582359236023612362236323642365236623672368236923702371237223732374237523762377237823792380238123822383238423852386238723882389239023912392239323942395239623972398239924002401240224032404240524062407240824092410241124122413241424152416241724182419242024212422242324242425242624272428242924302431243224332434243524362437243824392440244124422443244424452446244724482449245024512452245324542455245624572458245924602461246224632464246524662467246824692470247124722473247424752476247724782479248024812482248324842485248624872488248924902491249224932494249524962497249824992500250125022503250425052506250725082509251025112512251325142515251625172518251925202521252225232524252525262527252825292530253125322533253425352536253725382539254025412542254325442545254625472548254925502551255225532554255525562557255825592560256125622563256425652566256725682569257025712572257325742575257625772578257925802581258225832584258525862587258825892590259125922593259425952596259725982599260026012602260326042605260626072608260926102611261226132614261526162617261826192620262126222623262426252626262726282629263026312632263326342635263626372638263926402641264226432644264526462647264826492650265126522653265426552656265726582659266026612662266326642665266626672668266926702671267226732674267526762677267826792680268126822683268426852686268726882689269026912692269326942695269626972698269927002701270227032704270527062707
  1. #!/usr/bin/env python3
  2. # -*- coding: utf-8 -*-
  3. """windscada serve (M7c): 交互后端 — 工作台(全场)+分台页, 全部带时间窗.
  4. 设计: ①台页只渲染"有问题的系统"的证据图(分册两型: A=判据线+摆幅标注, B=本台vs全场四分位带), 无问题不出图;
  5. ②台内交互查询: 任一契约通道×时间窗×聚合 即时出图; ③工作台: 系统分类问题(全系统)+故障统计梳理+SOP控制策略;
  6. ④时间窗: 月度类数据真窗; 判级类(变桨/偏航/蓄能=近90日, 曲线=2025H2判别窗)固定时间窗如实标注.
  7. 端口 8033. 判级=筛查级/候选, P级归综合报告线."""
  8. try:
  9. from app_common.app_common_guanlan.api import install_root as _install_root
  10. except ImportError: # 极端兜底(包结构异常时按位置上跳)
  11. from pathlib import Path as _P
  12. def _install_root(_f): return _P(_f).resolve().parents[3]
  13. import sys as _sys, pathlib as _plb
  14. _sys.path.insert(0, str(_install_root(__file__)))
  15. from src import paths as P
  16. import sys, json, pathlib, re, threading, urllib.parse, collections
  17. import hashlib, time
  18. sys.path.insert(0, str(_install_root(__file__))) # 安装根入 sys.path (先于 src.* 导入)
  19. from src import paths as _P # 路径唯一真源 (与 cwd 无关)
  20. import numpy as np, pandas as pd
  21. import math
  22. def _jclean(o):
  23. """NaN/Inf → null (非法 JSON 防线; pd.DataFrame 会把 None 变 NaN)."""
  24. if isinstance(o, float) and (math.isnan(o) or math.isinf(o)): return None
  25. if isinstance(o, dict): return {k: _jclean(v) for k, v in o.items()}
  26. if isinstance(o, (list, tuple)): return [_jclean(x) for x in o]
  27. return o
  28. # ── 客户交付脱敏 (2026-08-28 demo 门禁逮出 7 处特征频率 + 21 处判据门槛上屏) ──
  29. # 规则来自 memory client-deliverable-hide-core-parameters: 判据名与实测值可见, 数值门槛与特征频率隐去。
  30. # 分级而非一刀切 — 一刀切会把业主必须知道的口径也糊掉:
  31. # 脱 特征频率 (134.9Hz 是我们解出来的轴承几何, 属核心参数)
  32. # 脱 判据门槛 (≥40kW / |z|≥3 / ≥2.5× — 我们的判别阈值)
  33. # 留 时间口径 (≥1h 段 / ≥3 天 — 业主要据此理解统计范围)
  34. # 留 物理常数 (Cp ≤ 0.593 贝兹极限 — 公开知识, 不是我们的)
  35. # 留 台数与占比 (≥2 个系统报警 — 业主自己数得出来)
  36. # 内部模式: WINDSCADA_INTERNAL=1
  37. import os as _os
  38. INTERNAL = _os.environ.get('WINDSCADA_INTERNAL', '') == '1'
  39. # ★ Hz 后不能用 \b: 中文属于 \w, "597Hz线" 里 z 与 线 之间没有词边界 ⇒ \b 不成立, 实测漏脱 WTG05 的 597Hz。
  40. # ★ 前瞻也不能用 \w: 中文同属 \w, "HS_gen测点274Hz" 里 274 前是"点" ⇒ 断言失败漏脱。
  41. # 只需排除数字与小数点 (防止把 1274 截成 274)。
  42. _RE_HZ = re.compile(r'(?<![\d.])\d+(?:\.\d+)?\s*(?=Hz)')
  43. _KEEP_UNIT = ('h', '小时', '天', '日', '月', '年', '台', '次', '条', '起', '个')
  44. # ★ 判据符号前缀: |z|≥3 / σ≥2 / ×≥2.5 这类无单位阈值必须脱 —
  45. # 否则会被"无单位小整数=计数口径"那条规则放行 (实测 |z|≥3 漏脱)。
  46. _RE_THR = re.compile(r'([≥≤])\s*(\d+(?:\.\d+)?)\s*([a-zA-Z%℃×/·²]*)')
  47. _CRIT_MARK = ('z', 'σ', '×', '倍', 'dev', 'resid', '选择性', '峰', '结构度', 'ratio')
  48. def _redact_text(t):
  49. if INTERNAL or not isinstance(t, str) or not t:
  50. return t
  51. if '0.593' in t or '贝兹' in t or 'Betz' in t: # 物理常数不动
  52. return _RE_HZ.sub('▪▪', t)
  53. t = _RE_HZ.sub('▪▪', t)
  54. def _th(m):
  55. sign, num, unit = m.group(1), m.group(2), m.group(3) or ''
  56. pre = t[max(0, m.start() - 8):m.start()]
  57. if any(c in pre for c in _CRIT_MARK): # |z|≥3 / ×≥2.5 = 判据阈值, 必脱
  58. return f'{sign}▪{unit}'
  59. if unit in _KEEP_UNIT or (not unit and float(num) <= 10 and float(num) == int(float(num))):
  60. return m.group(0) # 时间/计数口径保留 (≥1h段 / ≥2 个系统)
  61. return f'{sign}▪{unit}'
  62. return _RE_THR.sub(_th, t)
  63. def _redact(o):
  64. if INTERNAL:
  65. return o
  66. if isinstance(o, str):
  67. return _redact_text(o)
  68. if isinstance(o, dict):
  69. return {k: _redact(v) for k, v in o.items()}
  70. if isinstance(o, list):
  71. return [_redact(v) for v in o]
  72. return o
  73. # ── 语言投影 (2026-08-31) ────────────────────────────────────────────
  74. # 诊断: 骨架页 HTML 里**一个中文字都没有** (实测 0 字符), 内容全部由 JS 从 API 拉取,
  75. # 而 API 不认 lang —— 所以"英文版"拿到的一直是中文数据。此前三轮翻译 HTML 源码,
  76. # 改的是一个没有内容的文件, 这是英文残留反复修不好的真正原因。
  77. # 修法: 在 API 唯一出口 _jdump 上做文本层投影。**不带 lang 时原样返回**, 中文版行为零变化。
  78. _REQ = threading.local()
  79. try:
  80. from src.ontology import codes as _codes
  81. except Exception: # 码表缺失不应让服务起不来
  82. _codes = None
  83. def _lang_of():
  84. return getattr(_REQ, 'lang', '') or ''
  85. def _enum_fallback(field, value, lang):
  86. """枚举值映射不出码时的回落 —— **走翻译, 不留中文**。
  87. 根因 (2026-08-31 实测): `part` 这类字段被当成枚举, 但它在工单里其实是**现场手写的
  88. 自由文本** ('齿轮箱高速轴' / 'IGBT、直流保险' / '(振动线未收录)' / '无')。码表只有 18 条,
  89. 覆盖不了; 而 normalize 失败时原实现**原样返回中文**, 于是英文页上直接显示中文。
  90. 把这些硬塞进码表是错的 —— 它们不是枚举语义, 是描述。正确做法是让它们走翻译层。
  91. 查不到才留原文 (不猜译)。
  92. """
  93. if lang == 'zh' or not isinstance(value, str) or not value.strip():
  94. return value
  95. try:
  96. from windscada import en_text as _et
  97. except Exception:
  98. return value
  99. for scope in (field, 'part', None):
  100. try:
  101. en, _src = _et.lookup(value, scope) if scope else _et.lookup(value)
  102. except Exception:
  103. en = None
  104. if en:
  105. return en
  106. return value
  107. # ★这些字段 **不翻译** —— 它们是 JS 的判别值与查表键, 不是给人读的文案。
  108. # 翻了会静默劈叉: 页面上看不出错, 但整段内容消失或整片颜色变灰。
  109. # 2026-09-02 线上实测 (用户逮「有问题的看不见? 只有 good」):
  110. # · st 被翻成 Alarm → `rows.filter(r=>r.st==='报警')` 零命中 → **报警 12 台整表不渲染**,
  111. # 而柱状条读的是 sysdist (**键**没被翻) 仍显示 "Alarm 12" ⇒ 数字与明细自相矛盾。
  112. # · ostate 被翻成 Good/Excellent → `C[c.ostate]` (C 是中文键) 全部 MISS
  113. # → 融合矩阵**每一格颜色都退成灰**。
  114. # 判据: 凡是 JS 里 `X[obj.字段]` 查中文键表、或 `obj.字段==='中文'` 比较的字段, 一律进本表。
  115. # 显示层 (ui_en.value_script) 负责把渲染出来的中文文本换成英文, 所以人看到的仍是英文。
  116. # tests/test_keep_zh_discriminators.py 扫源码断言本表齐全, 防以后加字段漏登记。
  117. _KEEP_ZH = frozenset((
  118. 'st', 'ostate', 'wlevel', '状态', '设备状态', '证据状态', '部件', 'liveness',
  119. 'loc', # 模型本地/云端档位: `v.loc==='本地'` 决定"数据不出本地"这句能不能说
  120. # ★第三类判别用法, 本闸的扫描器 (tests/test_keep_zh_discriminators.py) 抓不到:
  121. # 不是比较也不是查表, 而是拼进 URL 路径当 ID (`/problem/${t}/${encodeURIComponent(x)}`)。
  122. # syss 被翻成 'Yaw' 后, 浮层请求 sys=Yaw 永远查不到烤的是 sys=偏航 的快照键,
  123. # 点开必是"本时间窗无记录"——页面不报错, 数据静默消失, 和 st/ostate 那两类同一后果不同机制。
  124. # 2026-09-03 用户实测逮出 (Watch list → 系统名链接)。显示文字仍由 sysName(x) 客户端译。
  125. 'syss',
  126. ))
  127. # ★按**值**豁免 —— 给 `k` 这种一名多用的字段用。
  128. # `k` 在告警排行里是显示串 (`"7122 偏航液压油位低"`, 要翻),
  129. # 在控制状态与部件配色里是枚举键 (`'停机'` / `'齿轮箱'`, 不能翻)。
  130. # 键名分不开两者, 只能按值。2026-09-02 实测: 我先把 'k' 整个放进 _KEEP_ZH,
  131. # 结果**报警名全部不再翻译** —— 修一个洞挖一个洞, 按值豁免才收得住边界。
  132. # 来源: JS 里 CO / LOSSPLAIN 两张中文键表的键 + `x.k==='停机'` 比较。
  133. _KEEP_ZH_VALS = frozenset((
  134. '齿轮箱', '发电机', '主轴承', '变桨叶片', '偏航', '测风', '其他', # CO 配色键
  135. '限电·绑定', '限电命令·未绑定', '停机(调度令)', '停机', '正常发电', '低风待机', # LOSSPLAIN 键
  136. ))
  137. def _proj(o, lang):
  138. """把 `key_zh`/`key_en` 投影成裸 `key`。
  139. 英文侧取不到 `_en` 时写 `[translation pending]` 而**不回落中文** ——
  140. 回落会让缺口永远不被发现, 正是英文版残留三轮没被逮住的机制。
  141. """
  142. if isinstance(o, list):
  143. return [_proj(x, lang) for x in o]
  144. if not isinstance(o, dict):
  145. return o
  146. suf = '_' + lang
  147. other = '_zh' if lang == 'en' else '_en'
  148. # ★裸 `zh` / `en` 成对 (framework.py 的 SYSTEMS/RECORDS/stages 就是这个写法)。
  149. # 它们不带下划线, 走不到下面的后缀规则; 而末尾"残留中文一律写 [translation pending]"
  150. # 的兜底会把 zh 打成 pending —— **旁边的 en 明明就是对的**。实测 145 处 pending
  151. # 里的大多数是这一种。消费端 JS 读的是 `.zh`, 所以英文版把 en 的值投影到 zh 上。
  152. if lang != 'zh' and isinstance(o, dict) and 'zh' in o and 'en' in o:
  153. o = dict(o, zh=o['en'])
  154. # ★只有 k_en **有值**时才让裸 k 让位 (2026-09-01 parity 闸逮)。
  155. # 原来无条件让位: 上游给了 action_en=None, 裸 action 就被顶成 None ——
  156. # 页面上显示 "None", 比留中文糟得多。空的语言副本等于没给, 不该覆盖原值。
  157. out = {}
  158. source_fields = {'name', '步骤', '故障', '排查', '原因', '处置', 'check'}
  159. originals = {key: o.get(key + '_zh', o.get(key)) for key in source_fields} if lang == 'en' else {}
  160. originals = {key: value for key, value in originals.items() if isinstance(value, str) and _CJK.search(value)}
  161. langed = {k[:-len(suf)] for k in o if k.endswith(suf) and o[k] not in (None, '', [])}
  162. for k, v in o.items():
  163. if k == '_original':
  164. out[k] = dict(v) if isinstance(v, dict) else v
  165. continue
  166. if k.endswith(other): # 另一语言的副本不出 API
  167. continue
  168. if k.endswith(suf):
  169. if v not in (None, '', []):
  170. out[k[:-len(suf)]] = _proj(v, lang)
  171. continue
  172. if k in langed: # 已有 key_<lang>, 裸 key 让位
  173. continue
  174. v2 = _proj(v, lang)
  175. if k in _KEEP_ZH or (isinstance(v2, str) and v2 in _KEEP_ZH_VALS):
  176. out[k] = v2 # 判别值/查表键: 原样留中文 (见 _KEEP_ZH 注释)
  177. continue
  178. # 枚举层先查码表 — 这些值有确定译名, 不该退化成 [translation pending]。
  179. # 列表型枚举 (如 facets) 逐项查。
  180. if lang != 'zh' and _codes is not None and k in _codes.CODES:
  181. if isinstance(v2, str):
  182. c = _codes.normalize(k, v2)
  183. if c:
  184. v2 = _codes.label(k, c, lang)
  185. else:
  186. v2 = _enum_fallback(k, v2, lang)
  187. elif isinstance(v2, list):
  188. v2 = [(_codes.label(k, _codes.normalize(k, x), lang)
  189. if isinstance(x, str) and _codes.normalize(k, x) else x) for x in v2]
  190. # 单数字段名对应复数码表 (facet ← facets) 也走一遍
  191. elif lang != 'zh' and _codes is not None and k.rstrip('s') in _codes.CODES:
  192. f = k.rstrip('s')
  193. if isinstance(v2, str):
  194. c = _codes.normalize(f, v2)
  195. if c:
  196. v2 = _codes.label(f, c, lang)
  197. elif isinstance(v2, list):
  198. v2 = [(_codes.label(f, _codes.normalize(f, x), lang)
  199. if isinstance(x, str) and _codes.normalize(f, x) else x) for x in v2]
  200. # ★这里原来只判 isinstance(v2, str) —— **字符串列表整个绕过了英文投影**。
  201. # _proj 递归进 list, 但 _proj(str) 原样返回, 所以 ["偏航","变桨"] 这种
  202. # 一路穿到 API。实测漏掉 5 条路径 (parts[] / syss[] / srcs[] /
  203. # open_items[] / blind[]), 是英文页数据层残留中文的唯一来源。
  204. if lang == 'en':
  205. # Action counts use display labels as dictionary keys, so _en_deep
  206. # cannot reach them. Only this display map uses the established
  207. # work-order action glossary; all other dictionary/URL IDs stay put.
  208. if k == '动作统计' and isinstance(v2, dict):
  209. from src.windscada.maint.workorder_en import ACTION
  210. v2 = {ACTION.get(action) or action: count for action, count in v2.items()}
  211. v2 = _en_deep(v2, k)
  212. out[k] = v2
  213. for key, original in originals.items():
  214. if key in out and out[key] != original:
  215. out.setdefault('_original', {})[key] = original
  216. return out
  217. def _en_deep(v, field='', _depth=0):
  218. """对字符串与**字符串列表/嵌套列表**一律走 _data_text_en。
  219. 只对 str 生效的旧写法会让 list-of-str 静默穿透 —— 这类漏法的特征是
  220. "抽查单值全对, 页面上却成片中文", 因为漏的是容器不是值。
  221. """
  222. if isinstance(v, str):
  223. return _data_text_en(v, field) if _CJK.search(v) else v
  224. if isinstance(v, list) and _depth < 4:
  225. return [_en_deep(x, field, _depth + 1) for x in v]
  226. return v
  227. _DATA_EN_CACHE = {}
  228. _DATA_EN_MISS = set() # 查不到的原文集合 (与命中缓存分开: 未命中的返回值随请求模式而异)
  229. def _data_text_en(v, field=''):
  230. """数据层中文文本 → 英文。取不到才写 [translation pending]。
  231. ★为什么要有这一层: 原实现是"残留中文一律 [translation pending]" —— 干净但过于绝对。
  232. 实测 /api/fleet 一次就打出 **2199 处** pending, 而其中大量是厂商报警名、工单原文、
  233. 部件名、判级词, 这些在本项目里**已经有译源**: 厂商中英对照表 3834 条、IPC 部件表
  234. 2113 条、枚举码表、工单名构词表。不查就直接标 pending, 等于把已有的资产扔掉,
  235. 还让英文页看起来像没做完。
  236. 查不到仍写 pending —— **不猜译**。这条不能松: 猜译出来的厂商报警名会误导现场。
  237. """
  238. if v in _DATA_EN_CACHE:
  239. return _DATA_EN_CACHE[v]
  240. if v in _DATA_EN_MISS:
  241. return v if getattr(_REQ, 'keep_orig', False) else '[translation pending]'
  242. out = None
  243. # ⓪⁻ 数据层**整串**直译表 (融合面盲区/未决项长文)。放在最前: 这些是整段人写的
  244. # 分析结论, 走模板反解或分层词表都只会碎成片段。
  245. try:
  246. from src.windscada import data_tpl_en as _tplN
  247. out = getattr(_tplN, 'NOTES', {}).get(v)
  248. if out is None:
  249. from src.windscada.en_text import standard_term
  250. out = standard_term(v)
  251. if out is None:
  252. out = getattr(_tplN, 'RECORD_EN', {}).get(v)
  253. except Exception:
  254. out = None
  255. if out is not None:
  256. _DATA_EN_CACHE[v] = out
  257. return out
  258. # ⓪ 「码 + 中文描述」→ 本体 AlarmCode 的权威 name_en (查表, 不翻译)
  259. try:
  260. from src.windscada import en_text as _et0
  261. out = _et0.alarm_en(v)
  262. except Exception:
  263. out = None
  264. if out is not None:
  265. _DATA_EN_CACHE[v] = out
  266. return out
  267. # ① 服务端拼串走**模板反解** (精确, 非翻译): 抽数字得形态 → 查英文模板 → 数字按序填回
  268. try:
  269. from src.windscada import data_tpl_en as _tpl
  270. out = _tpl.render(v)
  271. except ValueError as e:
  272. print('[_data_text_en] 模板占位符对不上, 已跳过不填: %s' % e, file=sys.stderr)
  273. out = None
  274. except Exception:
  275. out = None
  276. if out is not None:
  277. _DATA_EN_CACHE[v] = out
  278. return out
  279. try:
  280. from src.windscada import en_text as _et
  281. # ① 分层查表 (厂商报警 → IPC 部件 → 枚举/文案)
  282. for sc in ('fault', 'part', 'ui', None):
  283. hit, _src = _et.lookup(v, sc) if sc else _et.lookup(v)
  284. if hit:
  285. out = hit
  286. break
  287. # ② 部件自由文本解析 (复合名/带注解/多部件连写)
  288. if out is None:
  289. out = _et.part_en(v)
  290. except Exception:
  291. pass
  292. # ②b 词表直查 —— 部件名/判级词/系统名作为**独立字段**出现时 (fus.表[].部件 / .振动结论 /
  293. # rel.表[].部件), 既不是模板串也不在 en_text 分层表里, 而 data_tpl_en 的
  294. # CHAN/GRADE2/SYSTEM_EN 早就登记了同一批词 (供模板内部用)。不查就整片留中文
  295. # (2026-09-03 实逮 fus.表 328 字)。**同一批词只因出现位置不同就译不出, 是接线漏不是词条漏。**
  296. if out is None:
  297. try:
  298. from src.windscada import data_tpl_en as _t2
  299. out = _t2.CHAN.get(v) or _t2.SYSTEM_EN.get(v)
  300. if out is None:
  301. g = _t2.GRADE2.get(v) or _t2.GRADE.get(v)
  302. if g: # 判级词在表里是小写短语 (模板内嵌用), 独立成字段要首字母大写
  303. out = g[0].upper() + g[1:]
  304. except Exception:
  305. pass
  306. if out is None:
  307. try:
  308. from src.windscada.maint import workname_en as _wn
  309. # ★只用 work_name_en (逐段都命中才拼, 命不中返回 None), **不用 category_en**。
  310. # category_en 是正则归类兜底: 它对看不懂的文本也会给出一个"像样"的答案。
  311. # 实测它把化验样注记译成 work item:
  312. # '末样2025-01-09 (距振动现窗>12月, stale: 只作检测链回测)'
  313. # → 'Condition monitoring — work item' ← 错, 那不是工单
  314. # '末样2026-07-13: 华标判正常(NB/T宽限); 重点值…'
  315. # → 'Generator — work item' ← 错且误导
  316. # 我在本函数的说明里写了"不猜译"然后把猜译器接了进来。
  317. # **看得懂就译, 看不懂就标 pending —— pending 是诚实的, 猜译不是。**
  318. out = _wn.work_name_en(v)
  319. except Exception:
  320. pass
  321. if out:
  322. _DATA_EN_CACHE[v] = out
  323. return out
  324. # 查不到: 经典页写 [translation pending] (显示层再藏); v2 (keep=1) **留原文** —— 前端按原文
  325. # 标注为"手册原文", 不丢信息也不冒充英文。命中与未命中分开缓存, 两种模式互不污染。
  326. _DATA_EN_MISS.add(v)
  327. return v if getattr(_REQ, 'keep_orig', False) else '[translation pending]'
  328. def maint_framework_view():
  329. """运维管理体系三维骨架 —— **唯一实现** (在线 API 与快照预烤共用)。
  330. ★同 dq_findings_view 的教训: 这段原来也是两份, 在线那份我接了 unmapped_top 的
  331. 英文化, 预烤那份没有 ⇒ 在线干净、快照里还是中文, 而上云的是快照。
  332. 抽出来一次, 消掉这类"改了一半"的可能。
  333. """
  334. import json as _j, pathlib as _p
  335. from src.windscada.maint import framework as _fw
  336. db = _j.loads(_P.objects_json().read_text(encoding='utf-8'))
  337. wo = [(v.get('props') or {}) for k, v in db.items() if k.startswith('workorder/')]
  338. r = dict(systems=_fw.system_load(wo)['rows'], tiers=_fw.systems_by_tier(),
  339. stages=_fw.stage_matrix(), records=_fw.RECORDS,
  340. gaps=_fw.capability_gaps(wo), coverage=_fw.coverage(wo))
  341. # unmapped_top = **未映射到系统码的现场手写部件名**, 不是枚举, 码表投影对它无效。
  342. # 走 part_en 自由文本解析; 空值在这个语境下写 "(not recorded)" 而非破折号。
  343. lang = _lang_of()
  344. if lang and lang != 'zh' and (r.get('gaps') or {}).get('unmapped_top'):
  345. try:
  346. from src.windscada import en_text as _et
  347. out = []
  348. for it in r['gaps']['unmapped_top']:
  349. if isinstance(it, (list, tuple)) and len(it) == 2:
  350. e = _et.part_en(it[0]) or it[0]
  351. out.append(['(not recorded)' if e == '\u2014' else e, it[1]])
  352. else:
  353. out.append(it)
  354. r['gaps'] = dict(r['gaps'], unmapped_top=out)
  355. except Exception as e:
  356. print('[maint_framework_view] unmapped_top 英文化失败: %s' % e, file=sys.stderr)
  357. return r
  358. def dq_findings_view():
  359. """数据可信度 finding 行 —— **唯一实现**。
  360. ★2026-08-31: 这段原来存在两份 —— 一份在 /api/dq_findings 处理分支, 一份在快照的
  361. `_dq()` 预烤闭包里。我给 API 那份补了 sensitivity / age_basis 等字段, 预烤那份没动
  362. ⇒ **在线看得到, 快照里没有**。而上云交付的恰恰是快照。
  363. 同型教训已在 memory `library-gate-not-wired-to-production`: 落库不等于接上生产路径。
  364. 这里直接消掉重复: 两边都调本函数。
  365. """
  366. import json as _j, pathlib as _p
  367. db = _j.loads(_P.objects_json().read_text(encoding='utf-8'))
  368. KEEP = ('finding/data_quality/', 'finding/caliber/',
  369. 'finding/system_coverage/', 'finding/tier_trend/')
  370. rows = []
  371. for k, o in db.items():
  372. if not k.startswith(KEEP) or not isinstance(o, dict):
  373. continue
  374. p2 = o.get('props') or {}
  375. rows.append(dict(
  376. id=k, kind=k.split('/')[1],
  377. name=p2.get('name_zh') or p2.get('name'),
  378. name_zh=p2.get('name_zh') or p2.get('name'),
  379. name_en=p2.get('name_en'),
  380. grade=p2.get('grade'), claim_window=p2.get('claim_window'),
  381. evidence=p2.get('evidence_zh') or p2.get('evidence') or [],
  382. evidence_zh=p2.get('evidence_zh') or [],
  383. evidence_en=p2.get('evidence_en') or [],
  384. caveat=p2.get('caveat'), caveat_en=p2.get('caveat_en'),
  385. impact=p2.get('discriminator_zh') or p2.get('disposition_zh')
  386. or p2.get('implication_zh') or p2.get('corrected_comparison_zh'),
  387. impact_en=p2.get('discriminator_en') or p2.get('disposition_en')
  388. or p2.get('implication_en') or p2.get('corrected_comparison_en'),
  389. got_wrong=p2.get('what_i_got_wrong_zh'),
  390. got_wrong_en=p2.get('what_i_got_wrong_en'),
  391. quality_flag=p2.get('quality_flag_zh'), quality_flag_en=p2.get('quality_flag_en'),
  392. robust=p2.get('robust_zh'), robust_en=p2.get('robust_en'),
  393. basis=p2.get('basis_zh'), basis_en=p2.get('basis_en'),
  394. sensitivity=(p2.get('sensitivity') or {}).get('note_zh'),
  395. sensitivity_en=(p2.get('sensitivity') or {}).get('note_en'),
  396. age_basis=p2.get('age_basis_zh'), age_basis_en=p2.get('age_basis_en'),
  397. ))
  398. rows.sort(key=lambda x: (x['kind'] != 'data_quality', x['id']))
  399. r = dict(rows=rows, n=len(rows))
  400. if not rows:
  401. # ★2026-09-18 用户令("页面只能基于输入数据重算")之后实逮: 清过产物再重算的机器上
  402. # `objects.json` 里 **0 条 finding/* 对象** ⇒ 本页只是"空", 不说原因就等于让人以为
  403. # "系统自查没问题"。这四族的生成端确实不在重算链上(见族表 guanlan/ontology 说明),
  404. # 所以如实写清, 并明确"不用旧版产出补"。
  405. r['note'] = ('本体对象库里没有 finding/* 对象 (四族: data_quality / caliber / system_coverage / '
  406. 'tier_trend) ⇒ 本页**如实为空**, 不是"自查无发现"。这些 finding 的生成端不在重算链上'
  407. '(唯一写方 scripts/ingest_ops_2025.py 的输入 data/raw/工作库 不在现场数据里), '
  408. '运行期不从交付包补齐(用户令 2026-09-17); 要出件须研发补生成端。')
  409. return r
  410. def _en_page(page):
  411. """lang=en 时把页面文案换成英文。**中文版不经过这里** (lang 空 → 直接返回原页)。
  412. 用的是 ui_en 的**整条短语表**, 不是 i18n_en 的逐词词表 —— 后者会产出
  413. "Verification 状态机" 这类混杂串 (已实测作废)。缺条目由 ui_en.audit 报出,
  414. 不静默留在页面上。
  415. """
  416. lang = _lang_of()
  417. if not lang or lang == 'zh':
  418. return page
  419. try:
  420. from src.windscada import ui_en
  421. page = ui_en.apply(page)
  422. # ★译后语法闸: 保护规则一松一紧都能把 JS 译坏, 而症状是**空白骨架页** ——
  423. # 此时"可见中文=0", 自查全绿。必须用外部解析器判, 不能靠数引号。
  424. for _m in re.finditer(r'<script[^>]*>(.*?)</script>', page, re.S):
  425. _code = _m.group(1)
  426. if _code.strip() and 'application/json' not in _m.group(0):
  427. ui_en.check_js(_code, label='英文页内联脚本')
  428. # 数据值 (报警/定论/工单…) 既做比较又要显示。比较照旧用中文, 显示层在 render
  429. # 之后替换文本节点 —— 只碰显示不碰逻辑, 所以判断不会失效。
  430. return page.replace('</body>', ui_en.value_script() + '</body>') \
  431. if '</body>' in page else page + ui_en.value_script()
  432. except Exception as e:
  433. # ★退回中文页比给个坏页安全, 但**必须响亮** —— 静默退回会让"英文版没做完"
  434. # 看起来像"英文版做完了只是有中文残留"。守护步骤失败必须可见。
  435. print('[_en_page] 🔴 英文化失败, 已退回中文页: %s: %s' % (type(e).__name__, e),
  436. file=sys.stderr)
  437. return page
  438. def _jdump(obj, **kw):
  439. o = _redact(_jclean(obj))
  440. lang = _lang_of()
  441. if lang and lang != 'zh': # 中文版走这行以外的路径, 一字不动
  442. o = _proj(o, lang)
  443. return json.dumps(o, allow_nan=False, **kw)
  444. def _gap_page(path: str, missing: 'pathlib.Path | None' = None, what: str = '这一块') -> str:
  445. r"""产物缺失时的兜底页 (2026-09-17 现场问题驱动)。
  446. 原来 `/static/…`、`/turbines/…` 取不到文件就回一个裸 `404` 文本 —— 浏览器里表现为**一块空白**,
  447. 用户看到的现象是"这个板块没有信息", 而不是"缺件"。最典型的就是总览页那个
  448. `<iframe src="/static/index.html">`(全场状态总览页)。
  449. ★2026-09-19 更新 (用户令「所有的计算均要形成观澜的源代码」): 该页**已有生成端** ——
  450. `scripts/windscada_overview_build.py`(用组件同一个快照烘焙器从重算产物生成), 重算链第 ⑦b 步会跑它,
  451. 所以现在"不在位"只意味着**还没重算**, 不再是"包内没有生成端"。本页说明文字已据此改写。
  452. 现在改成结构化说明: 缺哪一件、它归哪条产物线、**怎么补**、去哪儿看全局状态。
  453. 口径与 CMS 的"无产物"页一致 (docs §11.4: 无产物时必须还能打开页面, 且要写明缺件与补齐办法)。
  454. """
  455. import html as _h
  456. rel = ''
  457. try:
  458. rel = pathlib.Path(missing).relative_to(pathlib.Path(_P.ROOT)).as_posix() if missing else ''
  459. except Exception:
  460. rel = str(missing or '')
  461. return ('<!doctype html><meta charset="utf-8"><title>产物不在位</title>'
  462. '<style>body{font:14px/1.7 system-ui,"Microsoft YaHei",sans-serif;color:#183548;background:#f2f6f8;'
  463. 'margin:0;display:flex;align-items:center;justify-content:center;min-height:100vh}'
  464. '.b{background:#fff;border:1px solid #d4e3e7;border-radius:10px;padding:22px 26px;max-width:800px;'
  465. 'box-shadow:0 6px 20px rgba(14,74,87,.08)}h1{font-size:17px;margin:0 0 10px}'
  466. 'code{background:#eef4f6;border-radius:4px;padding:1px 5px}.s{color:#667f8e;font-size:13px}'
  467. 'a{color:#087f89}li{margin:3px 0}</style>'
  468. f'<div class="b"><h1>{_h.escape(what)}的数据产物不在位</h1>'
  469. f'<p class="s">请求路径 <code>{_h.escape(path)}</code></p>'
  470. + (f'<p>缺的是 <code>{_h.escape(rel)}</code>(属 <code>outputs/&lt;场站&gt;/</code> 下的产物)。</p>'
  471. if rel else '')
  472. + '<p>两种情形分开补:</p><ul>'
  473. '<li><b>有生成端的产物</b>(alarms / workorders / temp_monthly / 振动窗索引谱 / '
  474. '<code>windscada/index.html</code> …): 放好 '
  475. '<code>data/raw/&lt;场站&gt;/</code> 后跑 <code>python scripts/rebuild_all.py</code>, '
  476. '或点运维控制台的「数据重算」。</li>'
  477. '<li><b>仍缺生成端的那几族</b>(如 <code>sop/</code>、<code>paradigm_r1/</code> 里的人工底稿、'
  478. '<code>windcms/</code> 的厂家报告转录): 按用户令 2026-09-19「所有的计算均要形成观澜的源代码」'
  479. '正在逐族补(清单见 <code>docs/源代码化清单_v0.1.md</code>); 在此之前页面如实显示缺件, '
  480. '且按用户令 2026-09-17 <b>运行期不从交付包补齐</b>。'
  481. '(2026-09-19 已补上生成端的: 变桨面 <code>pitch/**</code> ← <code>scripts/pitch_face_build.py</code>; '
  482. '振动在升/闭环 <code>m5_cms_tcm/component_history.json</code>、三层基线 <code>baseline_38.json</code> '
  483. '← <code>scripts/component_history_build.py</code> / <code>scripts/baseline_38_build.py</code>; '
  484. '事实契约 <code>guanlan/**</code> ← 重算台账生成 claim; '
  485. '发布层 <code>ontology/release_r1|r2</code> ← <code>scripts/ontology_release_build.py</code>。)</li></ul>'
  486. '<p class="s">逐件来源见 <code>outputs/&lt;场站&gt;/_provenance.json</code>;'
  487. '缺口清单见 <code>docs/系统设计说明.md</code> §7 与 docs/重算缺口与补件清单_v0.1.md。</p>'
  488. '<p><a href="/v2" target="_top">← 返回工作台</a> · '
  489. '<a href="http://127.0.0.1:28084/ops" target="_top">运维控制台 · 数据重算</a> · '
  490. '<a href="http://127.0.0.1:28084/healthz" target="_top">全部模块状态</a></p></div>')
  491. from http.server import ThreadingHTTPServer, BaseHTTPRequestHandler
  492. from src.windscada.config import farm
  493. from src.windscada import taxonomy
  494. from src.windscada.subsys import temp_nbm, hydraulic, yaw as yawmod, pitch as pitchmod
  495. from src.windscada.perf import availability
  496. from src.windscada import i18n
  497. CFG = farm(); ST = pathlib.Path(CFG['store'])
  498. _CJK = re.compile(r'[\u4e00-\u9fff]')
  499. def vibcms_results():
  500. """windcms 评估报告 → 结果层转录 (2026-08-28 用户令: 生接改融合·只显示结果·模型与计算隐藏).
  501. 取最新一期报告md; '融合级(模型)/CMS红黄'两列=模型输出, 不出结果层; 分析功能留独立cms.
  502. ★2026-09-18 用户令"页面只能基于输入数据重算、不许用旧版产出补"之后实逮: 清过产物再重算的机器上
  503. `windcms/` 里没有 `报告_CMS振动状态评估报告_*.md`, 于是 `sorted(glob)[-1]` 抛
  504. **IndexError: list index out of range**, 页面看到的是"windcms 报告解析失败: list index out of range"
  505. —— 与 `ProductsMissing` 的教训同一个病: **缺产物被报成了程序坏了**。现在按缺产物如实回结构化的
  506. `no_report`, 并写清它由谁生成。★2026-09-19 再修: 报告步的生成端(`scripts/windcms.py report`)与它的
  507. 上游(`rudong_model_run.py`/`rudong_fusion_run.py`)都已随包并上链, 缺件只剩"还没重算到那一步"这一个原因,
  508. 故这里再附一句"重算是否正在跑"(读 run/ops_job.json), 免得把"正在跑"读成"系统坏了"。
  509. """
  510. reps = sorted(ST.parent.glob('windcms/报告_CMS振动状态评估报告_*.md'))
  511. if not reps:
  512. # ★2026-09-19 修: 原话还写着"本包六层链的 model_run/fusion 两步脚本未随包" —— 那两步已在
  513. # 2026-09-19 按口径重建并上链(scripts/rudong_model_run.py / rudong_fusion_run.py), 话说反了。
  514. # 真实原因只有两种: ① 重算还没走到 ④b 的报告步(清过产物后这一步要跑很久); ② windcms 产物被清且未重算。
  515. try:
  516. from src import opsjob as _oj
  517. _run = (' ' + _oj.running_text()) if _oj.running_text() else ''
  518. except Exception:
  519. _run = ''
  520. return dict(date=None, overall='', action=[], grades=[], diff=dict(only_cms=[], only_handoff=[]),
  521. no_report=True,
  522. note_zh='本机没有 CMS 振动状态评估报告 (由重算链 ④b 的 report 步生成: '
  523. 'python scripts/windcms.py report)。如实为空, 不用旧版产出补。' + _run,
  524. err='无产物: windcms/报告_CMS振动状态评估报告_*.md 不在位')
  525. try:
  526. rp = reps[-1]
  527. md = rp.read_text(encoding='utf-8') # 不写 encoding 会按系统 locale(cp936) 读 → 中文报告乱码
  528. date = rp.stem.rsplit('_', 1)[-1]
  529. overall = ''
  530. if '### 3.1' in md:
  531. seg = md.split('### 3.1')[1].split('###')[0]
  532. overall = next((l.strip() for l in seg.splitlines() if l.strip().startswith('全场')), '')
  533. def _tab(sec):
  534. return [[c.strip() for c in l.strip('|').split('|')] for l in sec.splitlines() if l.startswith('| WTG')]
  535. action = _tab(md.split('### 3.2')[1].split('## ')[0]) if '### 3.2' in md else []
  536. grades = []
  537. if '## 附录 A' in md:
  538. for r in _tab(md.split('## 附录 A')[1]):
  539. # 列: 机组/前/后/齿/发/综合/融合级(模型·隐)/CMS红黄(隐)/行动
  540. grades.append(dict(t=r[0], 主轴承前=r[1], 主轴承后=r[2], 齿轮箱=r[3], 发电机=r[4],
  541. 综合=r[5], 行动=r[8] if len(r) > 8 else ''))
  542. # 与融合矩阵(handoff 定谳链)口径差异台
  543. from src.windscada.subsys import fusion as _f
  544. ft, _ = _f.fusion_table()
  545. hset = {r.turbine for _, r in ft.iterrows() if _f.verdict_class(str(r.振动结论))[1] in ('bad', 'warn')}
  546. cset = {r[0] for r in action}
  547. diff = dict(only_cms=sorted(cset - hset), only_handoff=sorted(hset - cset))
  548. # ★2026-09-19: 把**这个窗真实覆盖的日期区间**一并回给页面。原来只回 date(出件日) 与表格,
  549. # 页面上看不到"数据从哪天到哪天" ⇒ 用户拿 2026-03~04 的导出却只看到 4 月的谱时,
  550. # 无法判断是"数据没呈现"还是"源件本身就那么点"(实测源件最早只到 2026-03-16, 逐台更晚)。
  551. _span, _win = '', ''
  552. try:
  553. from src.windcms import config as _wcfg, data as _wdata
  554. _w2 = _wcfg.farm(str(CFG.get('key') or 'rudong'))
  555. _span = _wdata.span_text(_w2) or ''
  556. _win = ','.join(_wdata.windows(_w2))
  557. except Exception:
  558. _span = _win = ''
  559. return dict(date=date, overall=overall, action=action, grades=grades, diff=diff,
  560. 时间范围=_span, 窗=_win)
  561. except Exception as e:
  562. return dict(error=f'windcms 报告解析失败: {e}') # 响亮, 不静默空表
  563. LOCK = threading.Lock()
  564. _CACHE = {}
  565. # ── 产物指纹 → **运行中自动重载** (2026-09-16 用户令: "重算后的产物可即时用于页面呈现") ──────
  566. # 背景 (实测): _load() 把产物**一次性**读进内存, 之后只在启动时读一次 ⇒ 重算完成后页面仍显示旧数,
  567. # 必须重启组件服务才更新。判定依据: 改产物后 /detail/api/fleet 等仍回旧值, 重启组件后才变。
  568. # 做法: 取"页面取数依赖的产物文件"的 (路径, mtime_ns, size) 摘要当指纹; _load() 前比一次
  569. # (带 TTL, 免得每请求都 stat 一遍), 指纹变了就重载 ⇒ 重算产物**不再需要重启服务**。
  570. # 也提供 /api/reload 显式重载 (运维控制台重算结束后可点一下)。
  571. _STAMP_TTL = 2.0 # 秒: 指纹有效期 (stat ~150 个文件 ≈ 1~3 ms, 不值得每请求都做)
  572. _STAMP = {'val': None, 'at': 0.0}
  573. def _product_files():
  574. """页面取数依赖的产物文件清单 (顺序稳定: 供指纹; 只列**产物**, 不含 reference/ 随包契约)。
  575. 覆盖 _load() 直读的件 + 它经 taxonomy/temp_nbm/hydraulic/fusion 间接读的件:
  576. windscada/*.parquet|csv|json · ontology/*.json · pitch/*.parquet
  577. m5_cms_tcm/{handoff_vibration_v2,component_history,baseline_38}.json · windcms/报告_CMS*.md
  578. """
  579. try:
  580. from src.windcms.config import cms_out as _cms_out # 与写侧同一解析口 (WINDCMS_OUT)
  581. _cms_dir = _cms_out()
  582. except Exception:
  583. _cms_dir = _P.cms()
  584. pats = ((ST, ('*.parquet', '*.csv', '*.json')),
  585. (_P.ont(), ('*.json',)),
  586. (_P.pitch(), ('*.parquet',)),
  587. (_P.m5(), ('handoff_vibration_v2.json', 'component_history.json', 'baseline_38.json')),
  588. (_cms_dir, ('报告_CMS振动状态评估报告_*.md',)))
  589. files = []
  590. for d, ps in pats:
  591. for pat in ps:
  592. files.extend(sorted(d.glob(pat)))
  593. return files
  594. def products_stamp(force=False):
  595. """产物指纹 (sha1 前 16 位)。force=True 时忽略 TTL 立即重算。"""
  596. now = time.time()
  597. if not force and _STAMP['val'] is not None and (now - _STAMP['at']) < _STAMP_TTL:
  598. return _STAMP['val']
  599. h = hashlib.sha1()
  600. for p in _product_files():
  601. try:
  602. st = p.stat()
  603. h.update(f'{_P.rel(p)}|{st.st_mtime_ns}|{st.st_size}\n'.encode('utf-8'))
  604. except OSError:
  605. h.update(f'{_P.rel(p)}|MISSING\n'.encode('utf-8'))
  606. _STAMP.update(val=h.hexdigest()[:16], at=now)
  607. return _STAMP['val']
  608. class ProductsMissing(RuntimeError):
  609. """产物缺失 (outputs/<场>/… 被清空或尚未生成)。
  610. 2026-09-11 人工检查空状态时逮: 原来 `temp_monthly.parquet` 一缺, `_load()` 直接抛 FileNotFoundError,
  611. 外层处理器把它当**未知异常**回 500 + traceback —— 人看到的是"服务器错误", 分不清"没算产物"还是"程序坏了"。
  612. 现在这类缺失单独成类, 由请求层回**结构化的无数据**(HTTP 200 + err=no_products), 与 `adhoc_query` 的
  613. `no_source` 一致。真 bug 仍然照旧 500 并带 traceback (不静默)。"""
  614. def __init__(self, name, path):
  615. super().__init__(f'产物缺失: {name} (期望 {path})')
  616. self.what, self.path = name, path
  617. def reload_products(reason=''):
  618. """显式清缓存 (下一次请求重载)。返回清前的指纹, 供日志/接口回显。"""
  619. with LOCK:
  620. old = _CACHE.get('__stamp')
  621. _CACHE.clear()
  622. print(f'[reload] 清产物缓存 ({reason or "手动"}) 旧指纹={old}', flush=True)
  623. return old
  624. def _load():
  625. with LOCK:
  626. stamp = products_stamp()
  627. if _CACHE.get('__loaded'):
  628. if _CACHE.get('__stamp') == stamp:
  629. return
  630. # 产物变了 (典型: 刚跑完重算) → 重载。**先建后换**: 下面任何一步抛异常都不会破坏旧缓存,
  631. # 请求照旧能用旧数 (降级但不空白), 同时日志留痕。
  632. print(f'[reload] 产物指纹变化 {_CACHE.get("__stamp")} → {stamp}, 重载', flush=True)
  633. tmp = {}
  634. try:
  635. for key, name in (('tm', 'temp_monthly.parquet'), ('al', 'alarms.parquet'),
  636. ('lm', 'loss_monthly.parquet'), ('bins', 'powercurve_bins.parquet'),
  637. ('pcd', 'powercurve_dev.parquet')):
  638. f = ST / name
  639. if not f.exists():
  640. raise ProductsMissing(name, f)
  641. tmp[key] = pd.read_parquet(f)
  642. tmp['al']['month'] = tmp['al']['t_on'].dt.to_period('M').astype(str)
  643. tmp['pcd'] = tmp['pcd'].set_index('turbine')
  644. for key, name in (('duty', 'duty_monthly.parquet'),):
  645. f = ST / name
  646. tmp[key] = pd.read_parquet(f) if f.exists() else None
  647. try:
  648. tmp['sysmx'] = taxonomy.system_matrix()
  649. tmp['treg'] = temp_nbm.registry()
  650. tmp['treg'] = tmp['treg'][0] if isinstance(tmp['treg'], tuple) else tmp['treg']
  651. tmp['hyd'], _ = hydraulic.registry()
  652. tmp['hyd'] = tmp['hyd'].set_index('turbine')
  653. except FileNotFoundError as e: # 这些派生件同样在产物仓里; 缺了就按"无产物"处理
  654. raise ProductsMissing(getattr(e, 'filename', '派生产物'), getattr(e, 'filename', ST))
  655. zp = _P.pitch() / 'pitch_zero_monthly.parquet'
  656. tmp['zero'] = pd.read_parquet(zp) if zp.exists() else None
  657. tmp['__stamp'] = stamp
  658. tmp['__loaded'] = True
  659. _CACHE.update(tmp) # 原子提交: 失败时不留下半截缓存
  660. except Exception:
  661. if _CACHE.get('__loaded'):
  662. print('[reload] 重载失败 → 继续用上一份缓存 (页面不会空白, 但数是旧的; 看上面的异常)', flush=True)
  663. raise
  664. WINDOWS = ['近30日', '近90日', '2026年', '2026H1', '2025H2', '全程']
  665. def _month_end(m):
  666. """'YYYY-MM' → 该月最后一天 'YYYY-MM-DD'。"""
  667. import calendar
  668. y, mo = int(m[:4]), int(m[-2:])
  669. return f'{m}-{calendar.monthrange(y, mo)[1]:02d}'
  670. def months_of(win):
  671. """窗 → 覆盖到的月份列表(月度类数据按此过滤)。"""
  672. _load()
  673. all_m = sorted(_CACHE['tm'].month.unique())
  674. if win == '全程': return all_m
  675. if win == '2025H2': return [m for m in all_m if '2025-07' <= m <= '2025-12']
  676. if win == '2026H1': return [m for m in all_m if '2026-01' <= m <= '2026-06']
  677. if win == '2026年': return [m for m in all_m if m >= '2026-01']
  678. import re as _re
  679. if _re.fullmatch(r'\d{4}-\d{2}', win): # 单月窗 (2026-08-27 用户令: 按月份过滤和选择)
  680. return [m for m in all_m if m == win]
  681. m2 = _re.fullmatch(r'(\d{4}-\d{2})~(\d{4}-\d{2})', win)
  682. if m2: # 月区间窗 "2025-04~2025-10"
  683. return [m for m in all_m if m2.group(1) <= m <= m2.group(2)]
  684. m3 = _re.fullmatch(r'(\d{4}-\d{2}-\d{2})~(\d{4}-\d{2}-\d{2})', win)
  685. if m3: # ★自定义起止(含)→ 相交的月 (用户令 2026-09-21)
  686. a, b = sorted((m3.group(1), m3.group(2)))
  687. return [m for m in all_m if _month_end(m) >= a and f'{m}-01' <= b]
  688. n = 1 if win == '近30日' else 3
  689. return all_m[-n:]
  690. def win_range(win):
  691. """窗 → (起, 止) 日期串,**含两端**(用户令 2026-09-21:时间窗支持自定义起止日期)。
  692. 预设/逐月/月区间窗一律落到"该时间窗覆盖月的月首/月末",于是与 `months_of()` 同口径;
  693. 日区间窗原样返回。日粒度件(停机事件/报警/日粒度判据)用它做**含端**过滤。
  694. """
  695. import re as _re
  696. if win:
  697. if _re.fullmatch(r'\d{4}-\d{2}-\d{2}~\d{4}-\d{2}-\d{2}', win):
  698. a, b = win.split('~')
  699. return (a, b) if a <= b else (b, a)
  700. if _re.fullmatch(r'\d{4}-\d{2}', win):
  701. return (f'{win}-01', _month_end(win))
  702. m2 = _re.fullmatch(r'(\d{4}-\d{2})~(\d{4}-\d{2})', win)
  703. if m2:
  704. return (f'{m2.group(1)}-01', _month_end(m2.group(2)))
  705. if win == '2025H2':
  706. return ('2025-07-01', '2025-12-31')
  707. if win == '2026H1':
  708. return ('2026-01-01', '2026-06-30')
  709. ms = months_of(win)
  710. if not ms:
  711. return ('0001-01-01', '9999-12-31')
  712. return (f'{ms[0]}-01', _month_end(ms[-1]))
  713. def in_win(ts, win):
  714. """按窗过滤**有日粒度**的时间列(含两端)。月度列请用 `months_of()`(闰月/跨月取整不同)。"""
  715. import pandas as pd
  716. a, b = win_range(win)
  717. t = pd.to_datetime(ts, errors='coerce')
  718. return (t >= pd.Timestamp(a)) & (t <= pd.Timestamp(b) + pd.Timedelta(days=1) - pd.Timedelta(seconds=1))
  719. def span_of(win):
  720. """窗 → 计算层用的 (起, 止) 元组(含);判级/曲线等按窗重算的口子都吃这个。"""
  721. return win_range(win)
  722. # 九系统的 slug(前端 hash 用)↔ 中文名(服务端取数用)。**必须双向可认**:
  723. # 2026-09-22 用户测试 3/5 实逮 —— v2 的「该问题页」「命中系统」链接过去传的是 slug(`pitch`),
  724. # 而 `/api/problem` 与页面按中文名取数 ⇒ 返回 0 条后页面还写"该系统在判级时间窗内无非常态"(对报警机组是**假陈述**)。
  725. SYS_SLUG = {'变桨': 'pitch', '偏航': 'yaw', '齿轮箱': 'gearbox', '发电机': 'generator',
  726. '主轴承': 'main_bearing', '变流器': 'converter', '塔架与基础': 'tower',
  727. '主控与传感网': 'control', '叶片与叶根': 'blade'}
  728. SYS_FROM_SLUG = {v: k for k, v in SYS_SLUG.items()}
  729. def sys_norm(x):
  730. """系统标识归一:中文名 / 前端 slug / 大小写变体 → 规范中文名;认不出返回 None。"""
  731. if not x:
  732. return None
  733. s = str(x).strip()
  734. if s in SYS_SLUG:
  735. return s
  736. low = s.lower().replace('-', '_')
  737. if low in SYS_FROM_SLUG:
  738. return SYS_FROM_SLUG[low]
  739. for zh, slug in SYS_SLUG.items(): # 兜底: 去空格/全半角差异后再比
  740. if s.replace(' ', '') == zh.replace(' ', '') or low == slug.replace('-', '_'):
  741. return zh
  742. return None
  743. # ── 按窗重算:后台算 + 缓存(用户令 2026-09-21)─────────────────────────────────
  744. # 判级矩阵 ≈25 s/窗、七镜头 ≈9~44 s/窗 —— 直接同步算会把 HTTP 请求挂几十秒(远端还要过网关,
  745. # 长请求有被掐的风险)。所以统一口径:**缓存命中即回;未命中就起后台线程,本次先回 pending,
  746. # 页面轮询再取**。同一窗只算一次;算完常驻内存(进程内)。
  747. import threading as _thr
  748. _WIN_CACHE: dict = {'sysmx': {}, 'curves': {}, 'm9': {}}
  749. _WIN_BUSY: dict = {'sysmx': set(), 'curves': set(), 'm9': set()}
  750. _WIN_ERR: dict = {'sysmx': {}, 'curves': {}, 'm9': {}}
  751. _WIN_LOCK = _thr.Lock()
  752. # ★重活**串行**(2026-09-21 服务器实逮): 预热与请求各自起线程 ⇒ 同时跑两个按窗重算(各几百 MB 级
  753. # 中间帧),小内存机器上进程被系统杀掉(日志停在"按窗重算启动"两行、随后网关 503)。
  754. # 现在全部重活排队执行:峰值 = 1 份中间帧。
  755. _HEAVY_LOCK = _thr.Lock()
  756. def _mem_mb():
  757. """本进程可用内存(MB);拿不到就返回 None(不因为这些诊断信息把重算搞挂)。"""
  758. try:
  759. import ctypes
  760. class _MS(ctypes.Structure):
  761. _fields_ = [('dwLength', ctypes.c_ulong), ('dwMemoryLoad', ctypes.c_ulong),
  762. ('ullTotalPhys', ctypes.c_ulonglong), ('ullAvailPhys', ctypes.c_ulonglong),
  763. ('ullTotalPageFile', ctypes.c_ulonglong), ('ullAvailPageFile', ctypes.c_ulonglong),
  764. ('ullTotalVirtual', ctypes.c_ulonglong), ('ullAvailVirtual', ctypes.c_ulonglong),
  765. ('ullAvailExtendedVirtual', ctypes.c_ulonglong)]
  766. st = _MS()
  767. st.dwLength = ctypes.sizeof(_MS)
  768. ctypes.windll.kernel32.GlobalMemoryStatusEx(ctypes.byref(st))
  769. return int(st.ullAvailPhys // (1024 * 1024))
  770. except Exception:
  771. return None
  772. def _win_key(win):
  773. a, b = span_of(win)
  774. return f'{a}~{b}'
  775. def _win_get(kind, win, fn):
  776. """→ 值 或 None(None = 正在算/刚起算)。失败**记名**(`_WIN_ERR`)而不是静默当"没数据"。"""
  777. key = _win_key(win)
  778. with _WIN_LOCK:
  779. if key in _WIN_CACHE[kind]:
  780. return _WIN_CACHE[kind][key]
  781. if key not in _WIN_BUSY[kind]:
  782. _WIN_BUSY[kind].add(key)
  783. def _run():
  784. try:
  785. with _HEAVY_LOCK: # 重活串行: 同一时刻只算一份
  786. m0 = _mem_mb()
  787. print(f'[win] {kind} 按窗重算开算 {key}(可用内存 {m0} MB)', flush=True)
  788. v = fn()
  789. print(f'[win] {kind} 按窗重算完成 {key}(可用内存 {_mem_mb()} MB)', flush=True)
  790. except Exception as e: # 守护失败必响亮
  791. _WIN_ERR[kind][key] = f'{type(e).__name__}: {e}'[:200]
  792. print(f'[win] {kind} 按窗重算失败 {key}: {_WIN_ERR[kind][key]}', flush=True)
  793. v = None
  794. with _WIN_LOCK:
  795. if v is not None:
  796. _WIN_CACHE[kind][key] = v
  797. _WIN_BUSY[kind].discard(key)
  798. _thr.Thread(target=_run, name=f'win-{kind}-{key}', daemon=True).start()
  799. print(f'[win] {kind} 按窗重算启动 {key}(首次约数十秒,页面会先出 pending)', flush=True)
  800. return None
  801. def _win_wait(kind, win, timeout=120.0):
  802. """等某窗的按窗重算落地(预热用)。超时/失败返回 None,绝不抛。"""
  803. import time as _t
  804. key = _win_key(win)
  805. t0 = _t.time()
  806. while _t.time() - t0 < timeout:
  807. with _WIN_LOCK:
  808. if key in _WIN_CACHE[kind]:
  809. return _WIN_CACHE[kind][key]
  810. if key in _WIN_ERR[kind]:
  811. return None
  812. _t.sleep(1.0)
  813. return None
  814. def sysmx_of(win):
  815. """→ (判级矩阵, pending)。`pending=True` 时返回的是"另一口径"的旧矩阵,页面必须如实标注。"""
  816. from src.windscada import taxonomy
  817. m = _win_get('sysmx', win, lambda: taxonomy.system_matrix(CFG, span=span_of(win)))
  818. if m is None:
  819. return _CACHE['sysmx'], True
  820. return m, False
  821. def curves_of(win):
  822. """→ (七镜头分箱件 或 None, pending)。按窗重算(干净窗就是正式产物本身)。"""
  823. from src.windscada.perf import curves as _cv
  824. a, b = span_of(win)
  825. if (a, b) == (_cv.WIN[0], '2025-12-31'): # 与正式产物的判别窗一致 ⇒ 直接用产物, 不重算
  826. import pandas as _pd
  827. return _pd.read_parquet(ST / 'curve_lenses.parquet'), False
  828. st = _win_get('curves', win, lambda: _cv.build_store(CFG, span=(a, b), write=False))
  829. return st, st is None
  830. def m9_of(win):
  831. """→ (控制参数一致性, pending)。按窗重算(走窄仓 ≈5 s/窗);干净窗直接用正式产物。"""
  832. from src.windscada.perf import control as _cm
  833. a, b = span_of(win)
  834. if (a, b) == (_cm.WIN[0], '2025-12-31'):
  835. return _cm.registry(CFG), False
  836. v = _win_get('m9', win, lambda: _cm.registry(CFG, span=(a, b)))
  837. if v is None:
  838. return _cm.registry(CFG), True # 未就位: 先给上一口径, 页面按 pending 标注
  839. return v, False
  840. TS = {'报警': 0, '危险': 0, '良好': 1, '不可判': 2, '优秀': 3}
  841. # ---------- 图规格 (分册两型) ----------
  842. def fig_temp(t, ch, win):
  843. tm = _CACHE['tm']; ms = months_of(win)
  844. g = tm[tm.channel == ch]
  845. piv = g.pivot_table(index='month', columns='turbine', values='med').reindex(ms)
  846. if t not in piv.columns or piv[t].dropna().empty: return None
  847. nm = i18n.ch_cn(ch, with_tag=False)
  848. return dict(kind='band', title=f'{t} {nm} 月度 (发电态P>500kW)', unit='℃', months=ms,
  849. self=[None if v != v else round(v, 2) for v in piv[t]],
  850. med=[None if v != v else round(v, 2) for v in piv.median(axis=1)],
  851. q1=[None if v != v else round(v, 2) for v in piv.quantile(.25, axis=1)],
  852. q3=[None if v != v else round(v, 2) for v in piv.quantile(.75, axis=1)],
  853. note='灰带=全场四分位区间, 深线=全场中位, 红线=本台')
  854. def fig_duty(t, ch, title, unit, win, thresholds=None):
  855. d = _CACHE['duty']
  856. if d is None: return None
  857. ms = months_of(win)
  858. g = d[d.channel == ch]
  859. piv = g.pivot_table(index='month', columns='turbine', values='val').reindex(ms)
  860. if t not in piv.columns or piv[t].dropna().empty: return None
  861. return dict(kind='band', title=f'{t} {title}', unit=unit, months=ms,
  862. self=[None if v != v else round(v, 2) for v in piv[t]],
  863. med=[None if v != v else round(v, 2) for v in piv.median(axis=1)],
  864. q1=[None if v != v else round(v, 2) for v in piv.quantile(.25, axis=1)],
  865. q3=[None if v != v else round(v, 2) for v in piv.quantile(.75, axis=1)],
  866. thresholds=thresholds or [], note='灰带=全场四分位区间')
  867. def fig_alarm_bars(t, codes, title, win):
  868. al = _CACHE['al']; ms = months_of(win)
  869. g = al[(al.turbine == t) & (al.code.isin(codes)) & (al.month.isin(ms))]
  870. if not len(g): return None
  871. cnt = g.groupby('month').size().reindex(ms).fillna(0)
  872. return dict(kind='bars', title=f'{t} {title} 月度条数', unit='条', months=ms,
  873. vals=[int(v) for v in cnt], note='时间窗内合计 %d 条' % int(cnt.sum()))
  874. def fig_curve(t):
  875. b = _CACHE['bins']; g = b[b.turbine == t].sort_values('bin')
  876. if not len(g): return None
  877. dv = _CACHE['pcd'].loc[t]
  878. return dict(kind='xy', title=f'{t} 功率曲线 vs fleet中位 [固定判别窗 2025-H2, 机舱风口径]', unit='kW',
  879. x=[float(x) for x in g['bin']], self=[round(float(v), 1) for v in g.p],
  880. med=[round(float(v), 1) for v in g.fleet_p], xlab='风速 m/s',
  881. note=f"偏差 {dv.dev_w:+.1%} · 风速计偏置 {dv.ws_bias:+.2f} m/s · {dv['判别']}")
  882. def fig_zero(t, win):
  883. z = _CACHE['zero']
  884. if z is None: return None
  885. tt = t.replace('WTG', '')
  886. g = z[(z.turbine.astype(str).str.contains(tt)) | (z.turbine == t)] if 'turbine' in z.columns else None
  887. if g is None or not len(g): return None
  888. ms = months_of(win)
  889. mcol = [c for c in g.columns if 'month' in c.lower() or '月' in c]
  890. if not mcol: return None
  891. g = g[g[mcol[0]].astype(str).isin(ms)]
  892. # 只画零位偏差本身 (样本数量纲不同, 混画会把判据线压扁 → 移入图注)
  893. n_med = int(g['n_rows'].median()) if 'n_rows' in g.columns and len(g) else 0
  894. return dict(kind='multiline', title=f'{t} 桨叶零位偏差 月度', unit='°', months=[str(x) for x in g[mcol[0]]],
  895. series=[dict(name='零位偏差', vals=[round(float(v), 3) if v == v else None for v in g['zero_dev']])],
  896. thresholds=[dict(v=0.6, label='报警 0.6°'), dict(v=-0.6, label='')],
  897. note=f'判据: 变桨分册 §3.2 (报警门 ±0.6°); 月样本中位 {n_med} 行')
  898. SYS_ALARM = {'偏航': yawmod.FAM, '变桨': {'蓄能缸检': ['3411', '3412', '3413'], '泵/过热': ['7106', '7100'], '油位/滤芯': ['7101', '7103', '7104']}}
  899. SYS_CODES = {
  900. '变桨': ['3130', '3145', '3173', '3224', '3225', '3411', '3412', '3413', '3418', '7100', '7101', '7102', '7103', '7104', '7106', '9306', '9309'],
  901. '偏航': ['10100', '10104', '10105', '10113', '10125', '10201', '10400', '10401', '7122', '1012'],
  902. '齿轮箱': ['4101', '4102', '4103'],
  903. '变流器': ['13902', '13905', '19024', '19170', '19285', '20000', '63036'],
  904. '主控与传感网': ['8108', '8110', '8171', '8184', '8185', '8186', '8187', '63209', '63004'],
  905. }
  906. def fig_sys_alarms(t, sysname, win):
  907. """该台该系统报警族 月度条数 (与工单处理史并读)."""
  908. codes = SYS_CODES.get(sysname)
  909. if not codes: return None
  910. ms = months_of(win)
  911. al = _CACHE['al']
  912. g = al[(al.turbine == t) & (al.code.isin(codes)) & (al.month.isin(ms))]
  913. if not len(g): return None
  914. v = g.groupby('month').size().reindex(ms).fillna(0)
  915. top = g.groupby('code').size().sort_values(ascending=False).head(3)
  916. note = '时间窗内合计 %d 条; top: ' % len(g) + '; '.join(f'{c} {i18n.alarm_label(c)[:12]}×{n}' for c, n in top.items())
  917. return dict(kind='bars', title=f'{t} {sysname}系统 报警月度条数', unit='条', months=ms,
  918. vals=[int(x) for x in v], note=note)
  919. def turbine_problems(t, win):
  920. """有问题的系统 → 证据图 (分册式); 无问题不出."""
  921. _load()
  922. sysd = _CACHE['sysmx'][t]
  923. treg = _CACHE['treg']
  924. out = []
  925. for s in taxonomy.SYSTEMS:
  926. st = sysd[s]['状态']
  927. if st == '优秀': continue
  928. figs = []
  929. # 温度证据
  930. for _, r in treg[(treg.turbine == t)].iterrows():
  931. ch = r['channel']
  932. nm = ch.replace('tmp_wtc_', '').replace('_mean', '').replace('_max', '')
  933. if taxonomy.TEMP2SYS.get(nm) == s and str(r['判']) != '—':
  934. f = fig_temp(t, ch, win)
  935. if f: f['title'] = f"{t} {i18n.ch_cn(ch, with_tag=False)} — {r['判']} ({r['dev_K']:+.1f}K)"; figs.append(f)
  936. # 系统专属证据
  937. if s == '变桨':
  938. hv = _CACHE['hyd'].loc[t]
  939. if hv['判'] != '—':
  940. f = fig_duty(t, 'HydPress_band', '变桨液压10min摆幅 月度', 'bar', win,
  941. thresholds=[dict(v=52, label='正常充压滞环≈52bar')])
  942. if f: f['title'] = f"{t} 蓄能行为 — {hv['判']}"; figs.append(f)
  943. for nm, codes in SYS_ALARM['变桨'].items():
  944. f = fig_alarm_bars(t, codes, f'变桨{nm}报警', win)
  945. if f: figs.append(f)
  946. f = fig_zero(t, win)
  947. if f: figs.append(f)
  948. if s == '偏航':
  949. yrow_d = dict(zip(['D1润滑不可能', 'D2过热', 'D2变流器警告', 'D2解缆', 'D3制动液压错误', 'D3油位低', 'D4偏航失败', 'D4风传感'], [0]*8))
  950. f = fig_duty(t, 'dot_wtc_YawLubPu_timeon', '偏航润滑泵运行时长 月度', '秒/日', win)
  951. if f: figs.append(f)
  952. for nm, codes in yawmod.FAM.items():
  953. f = fig_alarm_bars(t, codes, '偏航·' + nm.lstrip('D1234'), win)
  954. if f: figs.append(f)
  955. if s in ('主轴承', '齿轮箱', '发电机'):
  956. vf = vib_scalar_fig(t, s)
  957. if vf: figs.append(vf)
  958. # 振动驱动的问题 (taxonomy 振动映射) → 配套分析曲线: 部件双端温度 + 主轴承柱塞 duty
  959. if '振动:' in sysd[s]['依据']:
  960. vib_chs = {'主轴承': ['tmp_wtc_MBearGTm_mean', 'tmp_wtc_MBearHTm_mean'],
  961. '发电机': ['tmp_wtc_GenBeRTm_mean', 'tmp_wtc_GenBeGTm_mean'],
  962. '齿轮箱': ['tmp_wtc_IMSGenTm_mean', 'tmp_wtc_HSGenTmp_mean']}[s]
  963. for ch in vib_chs:
  964. f = fig_temp(t, ch, win)
  965. if f: f['note'] = (f.get('note') or '') + ' | 振动判级交叉轴: 温度=同机制果侧观测(观测面≠源数)'; figs.append(f)
  966. if s == '主轴承':
  967. f = fig_duty(t, 'din_wtc_MBLubPis_timeon', '主轴承润滑柱塞运行时长 月度 (因侧执行链)', '秒/日', win)
  968. if f: f['note'] = (f.get('note') or '') + ' | 29#签名: 单次时长塌→回升→2026-07单月突变(泵打不满)'; figs.append(f)
  969. # 振动摘要配传动链图形 (2026-08-28 用户令: 关联图形填入) — comp=依据里的部件短语, sev=结论级
  970. _why = sysd[s]['依据']
  971. _m = re.search(r'振动: ([^ ]+) ([^ (]+)', _why)
  972. _comp, _lv = (_m.group(1), _m.group(2)) if _m else (s, '')
  973. _sev = 'var(--bad)' if any(k in _lv for k in ('定论', '报警')) else ('var(--warn)' if '预警' in _lv else 'var(--ok)')
  974. figs.append(dict(kind='vib', title=f'{t} {s} 振动分析摘要',
  975. comp=_comp, sev=_sev,
  976. text=_why + ' | 谱级深潜与逐窗趋势 → 振动融合板块·振动分析入口 (windcms)'))
  977. if '配对差' in sysd[s]['依据']:
  978. pair_chs = {'主轴承': ['tmp_wtc_MBearHTm_mean', 'tmp_wtc_MBearGTm_mean'],
  979. '发电机': ['tmp_wtc_GenBeGTm_mean', 'tmp_wtc_GenBeRTm_mean']}.get(s, [])
  980. for ch in pair_chs:
  981. f = fig_temp(t, ch, win)
  982. if f: figs.append(f)
  983. if '润滑' in sysd[s]['依据']:
  984. ch = {'主轴承': 'din_wtc_MBLubPis_timeon', '齿轮箱': 'din_wtc_GBLubPis_timeon', '发电机': 'dot_wtc_GenLubPu_timeon'}[s]
  985. f = fig_duty(t, ch, f'{s}润滑 月度', '秒/日', win)
  986. if f: figs.append(f)
  987. if s == '主控与传感网':
  988. if '测风' in sysd[s]['依据'] or '风速计' in sysd[s]['依据']:
  989. f = fig_alarm_bars(t, yawmod.FAM['D4风传感'], '风传感故障报警', win)
  990. if f: figs.append(f)
  991. if '曲线判风速计' in sysd[s]['依据'] or '风速计偏置' in sysd[s]['依据']:
  992. f = fig_curve(t)
  993. if f: figs.append(f)
  994. # ★每个问题配: 该系统报警月度图 + 工单处理史 (用户令)
  995. try:
  996. from src.windscada.subsys import workorder as womod
  997. hist = womod.history(t, s)
  998. except Exception as e:
  999. hist = dict(n=0, rows=[], note=f'工单不可用 {e}'[:60], 动作统计={}, 窗='')
  1000. af = fig_sys_alarms(t, s, win)
  1001. if af: figs.append(af)
  1002. # 排版 (2026-08-28 用户令"排版更合理"): ①同标题图去重 (配对差/温度证据两处会重复挂同一通道)
  1003. # ②按阅读顺序排 — 结论定位(传动链) → 证据强度(振动标量) → 时序(温度/润滑/报警) → 曲线
  1004. seen, uniq = set(), []
  1005. for _f in figs:
  1006. k = (_f.get('kind'), _f.get('title'))
  1007. if k in seen: continue
  1008. seen.add(k); uniq.append(_f)
  1009. _ORD = {'vib': 0, 'vibbar': 1, 'band': 2, 'line': 3, 'bar': 4, 'xy': 5, 'text': 6}
  1010. figs = sorted(uniq, key=lambda x: _ORD.get(x.get('kind'), 9))
  1011. out.append(dict(sys=s, 状态=st, 依据=i18n.humanize(sysd[s]['依据']), figs=figs, orders=hist))
  1012. dv = _CACHE['pcd'].loc[t] if t in _CACHE['pcd'].index else None
  1013. if dv is not None and abs(dv.dev_w) > 0.03:
  1014. f = fig_curve(t)
  1015. out.append(dict(sys='功率曲线', 状态='报警' if 'A类' in str(dv['判别']) else '良好',
  1016. 依据=f"偏差 {dv.dev_w:+.1%}, 风速计偏置 {dv.ws_bias:+.2f} m/s → {dv['判别']}", figs=[f] if f else []))
  1017. out.sort(key=lambda x: TS.get(x['状态'], 9))
  1018. return out
  1019. LENSES = [('wsb', 'grd_wtc_ActPower_mean', 'L1 功率曲线 (风速-功率)', '风速 m/s (机舱风)', 'kW'),
  1020. ('wsb', 'cp', 'L6 风能利用系数 Cp (风速-Cp)', '风速 m/s', 'Cp'),
  1021. ('wsb', 'lam', 'L5 叶尖速比 λ (风速-λ)', '风速 m/s', 'λ'),
  1022. ('grb', 'tq', 'L4 转矩-转速 (发电机侧)', '发电机转速 rpm', 'kNm'),
  1023. ('pwb', 'tur_wtc_PitcPosA_mean', 'L2 功率-桨距 (控制律/削峰)', '功率 kW', '°'),
  1024. ('pwb', 'tur_wtc_GenRpm_mean', 'L3 功率-发电机转速 (饱和点)', '功率 kW', 'rpm'),
  1025. ('pwb', 'p3', 'L2b 三叶桨距极差 (集距/不平衡)', '功率 kW', '°'),
  1026. ('pwb', 'ratio', 'L7 转速比 gen/rot (传动链)', '功率 kW', '—')]
  1027. def curve_view(win='2025H2'):
  1028. """特性曲线多镜头 (SOP §4.6c): fleet 中位+四分位带 + 过绝对锚的离群台曲线.
  1029. ★2026-09-21 用户令「发电性能随时间窗变化」: 七镜头**按所选时间窗重算**(首次后台算,先回 pending),
  1030. 7 张月度时序图按所选窗的月份过滤。原口径(2025H2 干净判别窗)仍可作为窗之一被选中,
  1031. 选到它时直接用正式产物,不重算。
  1032. """
  1033. import pandas as pd
  1034. from src.windscada.perf import curves as cv
  1035. a, b = span_of(win)
  1036. ms = months_of(win)
  1037. store, pending = curves_of(win)
  1038. if pending:
  1039. return dict(building=True, win=win, span=[a, b], months=ms,
  1040. note=f'正在按所选时间窗 {a} ~ {b} 重算七镜头曲线(首次约 10~45 秒),完成后自动刷新。'
  1041. f'期间下方月度时序图已按所选窗过滤,可直接看。')
  1042. wlabel = f'{a} ~ {b}'
  1043. figs = []
  1044. for xc, yc, title, xlab, unit in LENSES:
  1045. r = cv.lens(xc, yc, store_df=store, win_label=wlabel)
  1046. if not r: continue
  1047. outs = sorted(r['离群'].items(), key=lambda kv: -abs(kv[1]['z']))[:3]
  1048. a = r['anchor']
  1049. series = [dict(name='全场中位', vals=r['fleet'])]
  1050. series += [dict(name=f"{t.replace('WTG','')}# {v['resid']:+.4g}{a['单位']}", vals=r['per_t'][t]) for t, v in outs]
  1051. figs.append(dict(kind='multiline', title=title, unit=unit, months=[f"{x:g}" for x in r['x']],
  1052. series=series, xlab=xlab,
  1053. note=(f"显著门 |残差|≥{a['门']}{a['单位']} ({a['说明']}) ∧ |z|≥3 → 离群 "
  1054. f"{'/'.join(t.replace('WTG','')+'#' for t,_ in outs) if outs else '无'}"
  1055. f" | 灰域=全场四分位, 逐档中位聚合, 已剥限电(只正常发电态)"
  1056. f" | 样本 {r['n_total']:,} 个10min点 / {r['n_turbines']} 台 / 最小档 {r['n_min_bin']:,}"
  1057. f" | 时间窗 {r['win']}(随所选时间窗)"),
  1058. band=dict(q1=r['q1'], q3=r['q3'])))
  1059. # 残差视图 (§4.6c④ per机必算同型残差): 绝对量尺度上四分位带只有几个像素, 残差面才看得见
  1060. if outs:
  1061. rs = [dict(name=f"{t.replace('WTG','')}# {v['resid']:+.4g}{a['单位']}",
  1062. vals=[None if (r['per_t'][t][i] is None or r['fleet'][i] is None) else round(r['per_t'][t][i] - r['fleet'][i], 4)
  1063. for i in range(len(r['x']))]) for t, v in outs]
  1064. q1d = [round(r['q1'][i] - r['fleet'][i], 4) for i in range(len(r['x']))]
  1065. q3d = [round(r['q3'][i] - r['fleet'][i], 4) for i in range(len(r['x']))]
  1066. figs.append(dict(kind='multiline', title=title.split(' (')[0] + ' · 同型机群残差 (本台−全场中位)', unit=unit,
  1067. months=[f"{x:g}" for x in r['x']], series=rs, xlab=xlab,
  1068. thresholds=[dict(v=a['门'], label=f"显著门 +{a['门']}{a['单位']}"), dict(v=-a['门'], label='')],
  1069. band=dict(q1=q1d, q3=q3d),
  1070. note=(f"虚线=物理绝对锚 ±{a['门']}{a['单位']} ({a['说明']}); 灰域=全场四分位残差带; 出带且过锚才算离群"
  1071. f" | 样本 {r['n_total']:,} 点 / {r['n_turbines']} 台 | 时间窗 {r['win']}")))
  1072. # ---- 时序件: 控制参数月度 (M9b) + 偏航动态月度 (M12) ----
  1073. def _ts(path, val, title, unit, xlab, picks=None, note='', thr=None):
  1074. f = ST / path
  1075. if not f.exists(): return None
  1076. d = pd.read_parquet(f)
  1077. if val not in d.columns: return None
  1078. piv = d.pivot_table(index='month', columns='turbine', values=val)
  1079. # ★2026-09-21 用户令: 月度时序图按**所选时间窗**的月份过滤(月度件无法按日切,如实写"按月取整")
  1080. if ms:
  1081. piv = piv.reindex([m for m in [str(x) for x in piv.index] if m in set(ms)])
  1082. if not len(piv): return None
  1083. msx = [str(x) for x in piv.index]
  1084. med = piv.median(axis=1)
  1085. dev = (piv.sub(med, axis=0)).abs().mean()
  1086. picks = picks or list(dev.sort_values(ascending=False).head(3).index)
  1087. series = [dict(name='全场中位', vals=[None if v != v else round(float(v), 3) for v in med])]
  1088. series += [dict(name=t.replace('WTG', '') + '#', vals=[None if v != v else round(float(v), 3) for v in piv[t]])
  1089. for t in picks if t in piv.columns]
  1090. # 样本量与窗随图走 (2026-08-28 门禁 F1): 7 张时序图原来只有结论没有依据规模,
  1091. # 违反第一性原理②"Sample size declared"。逐月×台的非空格数即样本量。
  1092. _n = int(piv.notna().to_numpy().sum())
  1093. _nt = int(piv.shape[1])
  1094. note2 = (note + f" | 样本 {_n:,} 个月×台 / {_nt} 台 / {len(msx)} 个月"
  1095. f" | 时间窗 {msx[0]}~{msx[-1]}(月度聚合,随所选时间窗·按月取整)") if msx else note
  1096. return dict(kind='multiline', title=title, unit=unit, months=msx, series=series, xlab=xlab,
  1097. note=note2, thresholds=thr or [],
  1098. band=dict(q1=[None if v != v else round(float(v), 3) for v in piv.quantile(0.25, axis=1)],
  1099. q3=[None if v != v else round(float(v), 3) for v in piv.quantile(0.75, axis=1)]))
  1100. for args in (
  1101. ('control_monthly.parquet', 'p_cap', '时序① 满发功率封顶 月度 (取每月高位稳定值; 阶跃=配置变更)', 'kW', '月', None,
  1102. '灰域=全场四分位; 参数级分组见发电性能页 (4175/4200/4225 三组); 阶跃=配置变更, 平移=工况'),
  1103. ('control_monthly.parquet', 'w_cap', '时序② 转速封顶 月度 (取每月高位稳定值)', 'rpm', '月', None, '两轴独立分组 (1672/1680)'),
  1104. ('control_monthly.parquet', 'pitch_rated', '时序③ 额定段桨距角 月度 (标定漂移)', '°', '月',
  1105. ['WTG03', 'WTG34', 'WTG19'], '离群台=桨距调度同档差 top3'),
  1106. ('yaw_dynamic_monthly.parquet', 'err_sd', '时序④ 对风散布 σ 月度 (原始, 未清洗)', '°', '月',
  1107. ['WTG08', 'WTG16', 'WTG04'], '三台19个月全程高 → 非发作型; 清洗后降至5-7°, 判为机舱位置通道拖偏(A类数据质量)非对风故障'),
  1108. ('yaw_dynamic_monthly.parquet', 'travel_day', '时序⑤ 偏航行程 月度 (抗采样物理量)', '°/日', '月', None,
  1109. '行程=活动量代理, 非真磨损量; 与状态位次数互核'),
  1110. ('yaw_dynamic_monthly.parquet', 'twist_span', '时序⑥ 扭缆角月跨度 (ScYawPos 累计位置)', '°', '月', None, '解缆动作后回零'),
  1111. ('yaw_dynamic_monthly.parquet', 'pump_med', '时序⑦ 偏航泵压 月度 (小站泵压)', 'bar', '月', None, '全场极齐 → 无泵压异常'),
  1112. ):
  1113. f = _ts(*args)
  1114. if f: figs.append(f)
  1115. return dict(figs=figs, physics=cv.physics_check(),
  1116. note=('镜头口径 (SOP §4.6c): X轴优先功率(直测干净量); 机舱风 self-ref → 只判形状/相对, '
  1117. '绝对达成率须现场测风; per机异常必算同型机群残差(禁眼估); 工况段筛选=剔除0kW与低转速混合态档'))
  1118. _VIBSENS = {'主轴承': ['Main_bearing_front', 'Main_bearing_rear'],
  1119. '齿轮箱': ['Gear_planet', 'Gear_IMS', 'Gear_HS_rotor_side', 'Gear_HS_generator_side'],
  1120. '发电机': ['Generator_DE', 'Generator_NDE']}
  1121. _SENS_CN = {'Main_bearing_front': '主轴承前', 'Main_bearing_rear': '主轴承后', 'Gear_planet': '行星级',
  1122. 'Gear_IMS': '中间轴', 'Gear_HS_rotor_side': '高速轴转子侧', 'Gear_HS_generator_side': '高速轴发电机侧',
  1123. 'Generator_DE': '发电机DE', 'Generator_NDE': '发电机NDE'}
  1124. _VIBMET = [('Peak', '冲击峰值'), ('Kurtosis', '峭度'), ('Rms_Vel', '速度有效值'), ('iso_rms', 'ISO烈度')]
  1125. _VIBC = {}
  1126. def vib_scalar_fig(t, sysname):
  1127. """TCM 标量: 本台 vs 全场中位 (按测点×指标) + 六窗趋势. 数据=振动线 m5_cms_tcm 产物, 本页只转录不判级."""
  1128. import pandas as _pd
  1129. if 'z' not in _VIBC:
  1130. base = _P.m5()
  1131. try:
  1132. _VIBC['z'] = _pd.read_parquet(base / 'fleet_scalar_z.parquet')
  1133. _VIBC['tr'] = _pd.read_parquet(base / 'scalar_trend_6win.parquet')
  1134. except Exception:
  1135. _VIBC['z'] = _VIBC['tr'] = None
  1136. z, tr = _VIBC['z'], _VIBC['tr']
  1137. if z is None: return None
  1138. sens = _VIBSENS.get(sysname, [])
  1139. if not sens: return None
  1140. g = z[(z.turbine == t) & (z.sensor.isin(sens)) & (z.meas.isin([m for m, _ in _VIBMET]))]
  1141. if not len(g): return None
  1142. # 取最高功率档 (载荷高处信噪比最好); 该档缺则退最近档
  1143. order = sorted(g.bin.unique())
  1144. gb = g[g.bin == order[-1]]
  1145. rows = []
  1146. for sn in sens:
  1147. for mk, mcn in _VIBMET:
  1148. r = gb[(gb.sensor == sn) & (gb.meas == mk)]
  1149. if not len(r): continue
  1150. r = r.iloc[0]
  1151. tv = None
  1152. if tr is not None:
  1153. q = tr[(tr.turbine == t) & (tr.sensor == sn) & (tr.metric == mk)]
  1154. if len(q): tv = float(q.iloc[0]['ratio'])
  1155. rows.append(dict(测点=_SENS_CN.get(sn, sn), 指标=mcn, 本台=round(float(r.val), 3),
  1156. 全场中位=round(float(r.fleet_med), 3), z=round(float(r.z), 2),
  1157. 倍中位=round(float(r.val) / float(r.fleet_med), 2) if r.fleet_med else None,
  1158. 六窗比=round(tv, 2) if tv else None))
  1159. if not rows: return None
  1160. return dict(kind='vibbar', size='L', title=f'{t} {sysname} 振动标量 · 本台 vs 全场 ({str(order[-1]).replace("WPS-ActivePower ", "功率 ").rstrip(",")} kW 档)',
  1161. rows=rows, note='数据=TCM 六窗标量 (振动线 m5_cms_tcm 产物转录); z=fleet 稳健离群度, 六窗比=末窗/首窗 (>1 上升)')
  1162. def single_problem(t, sys, win):
  1163. """只回该台该系统这一个问题 (用户令: 点开就是一个问题, 其他不带出来)."""
  1164. return [p for p in turbine_problems(t, win) if p['sys'] == sys]
  1165. def fleet_view(win):
  1166. _load()
  1167. ms = months_of(win)
  1168. # ★2026-09-21 用户令「判级也按所选时间窗重算」: 判级矩阵按窗算(后台+缓存,首次 pending)。
  1169. sysmx, sysmx_pending = sysmx_of(win)
  1170. # ① 系统分类问题 (全系统)
  1171. systems, sysdist = {}, {}
  1172. for s in taxonomy.SYSTEMS:
  1173. rows = [dict(t=t, st=sysmx[t][s]['状态'], why=i18n.humanize(sysmx[t][s]['依据']))
  1174. for t in CFG['turbines'] if sysmx[t][s]['状态'] in ('报警', '不可判')] # 良好不出 (用户令: 只显示有问题的)
  1175. rows.sort(key=lambda r: TS.get(r['st'], 9))
  1176. systems[s] = rows
  1177. # 全场分布 (系统入口卡用: rows 只含问题台, 不能当全场分母 — 会显示"共12台")
  1178. from collections import Counter as _C
  1179. sysdist[s] = dict(_C(sysmx[t][s]['状态'] for t in CFG['turbines']))
  1180. # 良好台清单单列 (系统详情页第三档)
  1181. systems[s + '·良好'] = [dict(t=t, st='良好', why=i18n.humanize(sysmx[t][s]['依据']))
  1182. for t in CFG['turbines'] if sysmx[t][s]['状态'] == '良好']
  1183. # ①b 关注清单: ≥2 系统报警 (逐台独立台页)
  1184. watch = []
  1185. for t in CFG['turbines']:
  1186. al_sys = [x for x in taxonomy.SYSTEMS if sysmx[t][x]['状态'] == '报警']
  1187. if len(al_sys) >= 2:
  1188. watch.append(dict(t=t, n=len(al_sys), syss=al_sys,
  1189. # ★不截断: 截到 70 字会把依据切成半句, 而模板反解要完整形态
  1190. # ⇒ 英文侧整条落回中文 (2026-09-03 实逮, 总览 7 卡 4 张中文)。
  1191. # 显示长度归前端 CSS 管。
  1192. why=';'.join(f"{x}: {i18n.humanize(sysmx[t][x]['依据'])}" for x in al_sys)))
  1193. watch.sort(key=lambda r: -r['n'])
  1194. n_alarm_t = sum(1 for t in CFG['turbines'] if any(sysmx[t][x]['状态'] == '报警' for x in taxonomy.SYSTEMS))
  1195. # ② 故障统计梳理 (真时间窗)
  1196. al = _CACHE['al']; a = al[al.month.isin(ms)]
  1197. top_txt = a.groupby(['code', 'text']).agg(n=('code', 'size'), dur_h=('dur_s', lambda x: x.sum() / 3600)).reset_index()
  1198. pareto_n = top_txt.sort_values('n', ascending=False).head(12)
  1199. pareto_d = top_txt.sort_values('dur_h', ascending=False).head(12)
  1200. monthly = a.groupby('month').size().reindex(ms).fillna(0)
  1201. per_t = a.groupby('turbine').size().sort_values(ascending=False).head(10)
  1202. # 经典故障分析: 条数与时长必须配对看 (2026-08-28)。两个独立排行榜分不出
  1203. # "高频短时(信号抖动或重复触发)" 与 "低频长时(硬故障)" — 处置方向相反, 混在一起会误派工。
  1204. # 第三维=影响台数, 分离"单台刷屏"与"全场批次共性"。
  1205. # 残月识别 (2026-08-28 审核逮): 2026-07 只有 6 天数据(覆盖 19.4%), 却被当整月
  1206. # 并进"2026年"的占比与月度趋势 → 末柱"下降"是假象。占比类分母也因此不完整。
  1207. import calendar as _cal
  1208. _lmc = _CACHE['lm']; _lmc = _lmc[_lmc.month.astype(str).isin(ms)]
  1209. _cov = {}
  1210. for _mm in ms:
  1211. _sub = _lmc[_lmc.month.astype(str) == _mm]
  1212. if not len(_sub):
  1213. _cov[_mm] = 0.0; continue
  1214. _y, _mo = int(_mm[:4]), int(_mm[-2:])
  1215. _cal_h = _cal.monthrange(_y, _mo)[1] * 24
  1216. _cov[_mm] = round(float(_sub['rows_'].sum() / 6 / max(_sub['turbine'].nunique(), 1) / _cal_h), 3)
  1217. _nt = a.groupby(['code', 'text'])['turbine'].nunique().rename('nt')
  1218. _med = a.groupby(['code', 'text'])['dur_s'].median().rename('med_s')
  1219. qd = top_txt.set_index(['code', 'text']).join([_nt, _med]).reset_index()
  1220. qd = qd.sort_values('n', ascending=False).head(30)
  1221. faults = dict(
  1222. pareto_n=[dict(k=f"{r.code} {i18n.alarm_label(r.code, r.text)[:18]}", v=int(r.n)) for _, r in pareto_n.iterrows()],
  1223. pareto_d=[dict(k=f"{r.code} {i18n.alarm_label(r.code, r.text)[:18]}", v=round(float(r.dur_h), 1)) for _, r in pareto_d.iterrows()],
  1224. monthly=dict(months=ms, vals=[int(v) for v in monthly],
  1225. cov=[_cov.get(x, 0.0) for x in ms],
  1226. per_day=[round(float(monthly[x]) / max(_cov.get(x, 0.0) * _cal.monthrange(int(x[:4]), int(x[-2:]))[1], 1e-9), 1)
  1227. if _cov.get(x, 0) > 0.02 else None for x in ms]),
  1228. cov=_cov, cov_min=min(_cov.values()) if _cov else 1.0,
  1229. partial=[x for x in ms if _cov.get(x, 1) < 0.5],
  1230. quad=[dict(code=str(r.code), k=i18n.alarm_label(r.code, r.text)[:20], n=int(r.n),
  1231. h=round(float(r.dur_h), 1), nt=int(r.nt), med=round(float(r.med_s), 1))
  1232. for _, r in qd.iterrows()],
  1233. n_codes=int(len(top_txt)),
  1234. per_t=[dict(k=k, v=int(v)) for k, v in per_t.items()], total=int(len(a)))
  1235. # ③ SOP 控制策略
  1236. lm = _CACHE['lm']; l = lm[lm.month.astype(str).isin(ms)]
  1237. hrs = l.groupby('state')['rows_'].sum() / 6
  1238. loss = l.groupby('state')['loss'].sum() / 1000
  1239. hrs_tot = max(hrs.sum(), 1)
  1240. stop_h = hrs.get('停机', 0) + hrs.get('停机(调度令)', 0)
  1241. disp_h = hrs.get('停机(调度令)', 0)
  1242. lw_h = hrs.get('低风待机', 0)
  1243. # 口径与 availability.summary 单源一致: 1 − 停机/(总 − 调度令 − 低风待机)
  1244. avail = 1 - (stop_h - disp_h) / max(hrs_tot - disp_h - lw_h, 1)
  1245. pcd = _CACHE['pcd']
  1246. # 能量账闭合 (2026-08-28 经典分析): 四行工况表看不出"电量流向"。
  1247. # 理论可发 = 实际上网 + Σ各态损失; 过物理上限核 (不可超 38台×4.0MW×窗时长)。
  1248. _act = l.groupby('state')['act'].sum() / 1000
  1249. _wf = [dict(k=str(k), loss=round(float(loss.get(k, 0)), 0), act=round(float(_act.get(k, 0)), 0),
  1250. h=round(float(hrs.get(k, 0)), 0)) for k in hrs.index]
  1251. _A = float(_act.sum()); _L = float(loss.sum())
  1252. _cap = len(CFG['turbines']) * 4.0 * float(hrs_tot) / max(len(CFG['turbines']), 1)
  1253. # 逐月能量 (2026-08-28 报表端口): 汇报要环比与趋势, 全时间窗聚合给不出。
  1254. # 残月覆盖率一并带出 — 报表里拿残月和整月比环比会读反 (故障月度图已踩过一次)。
  1255. _em = []
  1256. for _mm in ms:
  1257. _sub = l[l.month.astype(str) == _mm]
  1258. if not len(_sub):
  1259. continue
  1260. _a1 = float(_sub['act'].sum()) / 1000
  1261. _l1 = float(_sub['loss'].sum()) / 1000
  1262. _em.append(dict(m=_mm, act=round(_a1, 0), loss=round(_l1, 0), theo=round(_a1 + _l1, 0),
  1263. loss_pct=round(_l1 / max(_a1 + _l1, 1) * 100, 1),
  1264. cov=_cov.get(_mm, 1.0),
  1265. eflh=round(_a1 / max(len(CFG['turbines']) * 4.0, 1), 0)))
  1266. energy = dict(monthly=_em, act=round(_A, 0), loss=round(_L, 0), theo=round(_A + _L, 0),
  1267. loss_pct=round(_L / max(_A + _L, 1) * 100, 1),
  1268. eflh=round(_A / max(len(CFG['turbines']) * 4.0, 1), 0),
  1269. cap=round(_cap, 0), cap_ok=bool(_A + _L < _cap),
  1270. items=sorted(_wf, key=lambda r: -r['loss']))
  1271. control = dict(
  1272. energy=energy,
  1273. states=[dict(k=k, h=round(float(v), 0), pct=round(float(v / hrs_tot * 100), 1), loss=round(float(loss.get(k, 0)), 0)) for k, v in hrs.items()],
  1274. avail=round(float(avail * 100), 1),
  1275. curve=dict(sigma=round(float(pcd.dev_w.std() * 100), 2),
  1276. cands=[dict(t=i, dev=round(float(r.dev_w * 100), 2), 判=r['判别']) for i, r in pcd.iterrows() if r['判别'] != '—'],
  1277. note='固定判别窗 2025-H2'),
  1278. anchors='限电命令面占比 2025≈0.197 → 2026≈0.769 (专项锚)')
  1279. m8 = None
  1280. try:
  1281. from src.windscada.perf import faults as fmod
  1282. if (ST / 'stop_events.parquet').exists():
  1283. m8 = dict(mtbf=fmod.mtbf_summary(ms), stop_pareto=fmod.stop_pareto(ms), seasonal=fmod.seasonal())
  1284. except Exception as e:
  1285. m8 = dict(err=str(e)[:120])
  1286. m9 = None
  1287. try:
  1288. from src.windscada.perf import control as cmod
  1289. if (ST / 'control_profile.parquet').exists():
  1290. m9, m9_pending = m9_of(win) # ★按窗重算(用户令 2026-09-21)
  1291. except Exception as e:
  1292. m9 = dict(err=str(e)[:120])
  1293. m9_pending = False
  1294. try:
  1295. from src.windscada.perf import reliability as relmod
  1296. rel = relmod.overview(CFG, span=span_of(win))
  1297. except Exception as e:
  1298. rel = dict(err=str(e)[:120])
  1299. try:
  1300. from src.windscada.subsys import fusion as fusmod
  1301. ftab, fmeta = fusmod.fusion_table(CFG, all_turbines=True)
  1302. mrows, mheads = fusmod.matrix(CFG, ms) # 整体时间窗: 矩阵现窗/记忆轨迹随全局窗选择器
  1303. cap = fusmod.capability()
  1304. ev, ev_bounds = fusmod.drivetrain_events(ms, CFG)
  1305. lvcnt = {}
  1306. for r in mrows:
  1307. lv = r['振动']['level']; lvcnt[lv] = lvcnt.get(lv, 0) + 1
  1308. fkpi = dict(定论=lvcnt.get('bad', 0), 预警=lvcnt.get('warn', 0), 候选观察=lvcnt.get('note', 0),
  1309. 销案正常=lvcnt.get('ok', 0), 未列=lvcnt.get('unlisted', 0),
  1310. 收录台=len(CFG['turbines']) - lvcnt.get('unlisted', 0), 全场=len(CFG['turbines']),
  1311. 机制链=sum(1 for r in mrows if '机制链' in r['润滑'].get('note', '')),
  1312. 系数=dict(wear=cap['classes'].get('wear_progressive', {}).get('coef'),
  1313. thermal=cap['classes'].get('thermal_acute', {}).get('coef'),
  1314. unknown=cap['classes'].get('unknown', {}).get('coef')))
  1315. # 矩阵补 报警/工单 两列 (2026-08-28 用户令) — 报警=时间窗内真窗计数+族分布; 工单=历史台账(窗见列头)
  1316. # 报警/工单两列只统计传动链三部件 (2026-08-28 用户令: 只说发电机/齿轮箱/主轴) —
  1317. # 本矩阵每行是"该台传动链状态", 全场报警计数会把变桨/偏航/主控的量混进来 (316 条里绝大多数与传动链无关)
  1318. _DT_PAT = '齿轮|齿箱|润滑油|油冷|滤芯|发电机|定子|绕组|滑环|主轴承|主轴|轴承'
  1319. _a = _CACHE['al']; _aw = _a[_a.month.isin(ms)]
  1320. _aw = _aw[_aw.text.str.contains(_DT_PAT, na=False)]
  1321. _acnt = _aw.groupby('turbine').size()
  1322. _arank = _acnt.rank(ascending=False, method='min')
  1323. _FAM = [('齿轮箱', '齿轮|齿箱|润滑油|油冷|滤芯'), ('发电机', '发电机|定子|绕组|滑环'),
  1324. ('主轴承', '主轴承|主轴')]
  1325. try:
  1326. from src.windscada.subsys import workorder as _wo
  1327. _wod = _wo.load(CFG)
  1328. except Exception:
  1329. _wod = None
  1330. for _r in mrows:
  1331. _t = _r['turbine']; _g = _aw[_aw.turbine == _t]
  1332. _fam, _taken = [], set()
  1333. _CODEFAM = {'3225': '变桨叶片'} # 文本为空的码走码号归族 (3225=变桨液压, 现场西门子资料对译表)
  1334. for _nm, _pat in _FAM:
  1335. _byc = _g.code.astype(str).map(_CODEFAM) == _nm
  1336. _sub = _g[~_g.index.isin(_taken) & (_g.text.str.contains(_pat, na=False) | _byc)]
  1337. _taken |= set(_sub.index)
  1338. if len(_sub): _fam.append(dict(k=_nm, n=int(len(_sub))))
  1339. _oth = len(_g) - sum(f['n'] for f in _fam)
  1340. if _oth > 0: _fam.append(dict(k='其他', n=int(_oth)))
  1341. # 刷屏闸: 单码占比过半 → 条数不代表"问题多"而是一个码在抖 (08# 4609/4918=93.7% 码3225 集中2026-01)
  1342. _burst = None
  1343. if len(_g):
  1344. _bc = _g.groupby('code').size().sort_values(ascending=False)
  1345. _share = float(_bc.iloc[0]) / len(_g)
  1346. if _share >= 0.5:
  1347. _sub2 = _g[_g.code == _bc.index[0]]
  1348. _burst = dict(code=str(_bc.index[0]), share=round(_share, 3), n=int(_bc.iloc[0]),
  1349. months=sorted({str(m)[:7] for m in _sub2.t_on.dt.to_period('M').astype(str)}),
  1350. med_s=float(_sub2.dur_s.median()) if 'dur_s' in _sub2 else None)
  1351. # 明细 (用户令"浮窗说明之前发生问题"): top 码 + 首末时间
  1352. _top = []
  1353. if len(_g):
  1354. for (_c2, _tx), _sub3 in _g.groupby(['code', 'text']):
  1355. _top.append(dict(code=str(_c2), text=str(_tx)[:34], n=int(len(_sub3)),
  1356. first=str(_sub3.t_on.min())[:10], last=str(_sub3.t_on.max())[:10]))
  1357. _top.sort(key=lambda x: -x['n'])
  1358. _r['报警'] = dict(n=int(_acnt.get(_t, 0)), rank=int(_arank.get(_t, 0)) if _t in _arank else None,
  1359. tot=len(CFG['turbines']), fam=_fam, burst=_burst, items=_top[:6])
  1360. if _wod is not None and len(_wod):
  1361. _w = _wod[_wod.turbine == _t]
  1362. _wtxt = (_w.get('故障名称', '').astype(str) + ' ' + _w.get('故障位置二级', '').astype(str)
  1363. + ' ' + _w.get('维修对象', '').astype(str) + ' ' + _w.get('元器件名称', '').astype(str))
  1364. _w = _w[_wtxt.str.contains(_DT_PAT, na=False)]
  1365. _acts = _w['维修动作'].replace('', pd.NA).dropna().value_counts() if len(_w) else None
  1366. _wi = []
  1367. for _, _wr in _w.sort_values('t_report', ascending=False).head(6).iterrows():
  1368. _wi.append(dict(date=str(_wr.get('t_report'))[:10], name=str(_wr.get('故障名称') or '')[:30],
  1369. act=str(_wr.get('维修动作') or ''), part=str(_wr.get('维修对象') or _wr.get('元器件名称') or '')[:14]))
  1370. _r['工单'] = dict(n=int(len(_w)),
  1371. acts=[dict(k=str(k), n=int(v)) for k, v in (_acts.head(3).items() if _acts is not None else [])],
  1372. items=_wi,
  1373. last=(lambda v: str(v.max())[:10] if len(v) else None)(_w.t_report.dropna()[_w.t_report.dropna().dt.year > 2000]) if len(_w) else None)
  1374. else:
  1375. _r['工单'] = dict(n=0, acts=[], last=None)
  1376. _wyr = ''
  1377. if _wod is not None and len(_wod):
  1378. _tv = _wod.t_report.dropna()
  1379. _tv = _tv[_tv.dt.year > 2000] # 空日期落 1970 epoch, 剔除后再报窗 (否则列头写"1970~")
  1380. if len(_tv): _wyr = f"{_tv.min():%Y}~{_tv.max():%Y}"
  1381. mheads = dict(mheads, 报警=f"传动链三部件 · {ms[0]}~{ms[-1]}" if ms else '传动链三部件',
  1382. 工单=f"传动链三部件 · 台账 {_wyr}" if _wyr else '台账不可用')
  1383. # 决策链进度盘 (2026-08-28 用户问"这个到底能怎么用"): 原来只有 29# 一台的静态六格,
  1384. # 是展示牌不是工作面 — 打开它做不了任何决定。改为对全部需跟踪台报"走到第几步、卡在哪、下一步做什么"。
  1385. # 判定规则单源在此, 前端只渲染不判断。
  1386. _CLOOP = {}
  1387. try:
  1388. _ch = json.loads((ST.parent / 'm5_cms_tcm' / 'component_history.json').read_text(encoding='utf-8'))
  1389. for _k, _v in _ch.get('summary', {}).items():
  1390. if _k.startswith('★换件闭环') and isinstance(_v, list):
  1391. for _r in _v:
  1392. _CLOOP[_r['turbine']] = _r
  1393. if _k.startswith('★当前在升') and isinstance(_v, list):
  1394. for _r in _v:
  1395. _CLOOP.setdefault(_r['turbine'], {}).update(rising=_r)
  1396. except Exception:
  1397. pass
  1398. _HIT = {'bad', 'warn', 'note'}
  1399. # ★一台可能有多行 (29#/17# 各 2 行: 齿轮箱 + 主轴承)。用 {turbine: row} 直接建字典会被后一行覆盖,
  1400. # 实测 29# 因此取到"齿轮箱二级行星内齿圈·监视"那行, 而它的决策链讲的是"主轴承+集中润滑泵·定论"
  1401. # ⇒ 浮窗里部件、温度通道、油样全是另一个部件的 (memory: 画图须确认各维度来自同一个体)。
  1402. _SEV = {'bad': 0, 'warn': 1, 'note': 2, 'stale': 3, 'ok': 4, 'unlisted': 5}
  1403. _tbi = {}
  1404. for _row in ftab.to_dict('records'):
  1405. _k = _row['turbine']
  1406. if _k not in _tbi or _SEV.get(_row.get('级'), 9) < _SEV.get(_tbi[_k].get('级'), 9):
  1407. _tbi[_k] = _row
  1408. board = []
  1409. for _r in mrows:
  1410. t = _r['turbine']
  1411. srcs = [k for k in ('振动', '温度', '润滑', '油液') if (_r.get(k) or {}).get('level') in _HIT]
  1412. lvl = (_r.get('振动') or {}).get('level')
  1413. tb = _tbi.get(t, {})
  1414. mech = ('机制链' in ((_r.get('润滑') or {}).get('note') or '')
  1415. or '机制链' in (tb.get('振动结论') or ''))
  1416. n_wo = (_r.get('工单') or {}).get('n') or 0
  1417. cl = _CLOOP.get(t) or {}
  1418. closed = bool(cl.get('verdict') == '恢复')
  1419. rising = cl.get('rising')
  1420. if not (lvl in _HIT or (srcs and lvl != 'ok')):
  1421. continue
  1422. # 六步: done=已达 / open=未达 / unknown=数据不足不可判
  1423. steps = [
  1424. dict(k='证据', st='done' if srcs else 'open', v=f"{len(srcs)} 源", d='/'.join(srcs)),
  1425. dict(k='机制', st='done' if mech else 'open',
  1426. v='已定性' if mech else '未定性', d='根因定到部件与失效模式'),
  1427. dict(k='判级', st='done' if lvl in _HIT else 'open',
  1428. v=tb.get('证据状态') or lvl or '—', d=f"设备状态 {tb.get('设备状态') or '—'}"),
  1429. dict(k='排期', st='done' if t == 'WTG29' else 'open',
  1430. v='情景A 2026-09' if t == 'WTG29' else '未排',
  1431. d='进检修排程沙盘 (风险×损失双轴)'),
  1432. # ★这一列一律"不可判", 不是偷懒: 工单台账止 2024-11 而判级时间窗在 2026,
  1433. # 台账里的历史工单不可能是针对本次问题的 ⇒ 用 n>0 判"已派工"会造出假的完成态。
  1434. dict(k='动作', st='unknown',
  1435. v=(f"台账 {n_wo} 单" if n_wo else '台账无'),
  1436. d='工单台账 2020~2024-11, 25-26 现场未提供 — 历史单不对应本次问题, 须向现场核实'),
  1437. # ★历史闭环不占本格: 05#/20# 2024~2025 换过齿轮箱且已验证恢复, 但本轮又有新的候选级证据,
  1438. # 把历史闭环填进第六格会出现"验收已达而机制未定"的自相矛盾行。历史闭环走 hist_loop 单独标记。
  1439. dict(k='验收', st='open', v='未验',
  1440. d='本轮动作执行后复测同一判据是否回落' + (
  1441. f" (该台 {cl.get('replace_date')} 有过一次闭环: {cl.get('component')} 降至 {cl.get('ratio')}×)" if closed else '')),
  1442. ]
  1443. _first = next((x for x in steps if x['st'] != 'done'), None)
  1444. stuck = _first['k'] if _first else None
  1445. stuck_kind = _first['st'] if _first else None
  1446. NEXT = {'证据': '补测: 四源全为常态, 先确认监测面有效',
  1447. '机制': '定性: 现场/取样把根因定到部件与失效模式 — 机制不清则排期与动作都无依据',
  1448. '判级': '送审: 按 RV-1 触发相应审级',
  1449. '排期': '排程: 进沙盘比情景, 出建议窗口',
  1450. '动作': '核实: 向现场调取 2025~2026 工单 — 台账止于 2024-11, 系统无法判定是否已派工',
  1451. '验收': '复测: 动作已执行, 复测同一判据是否回落'}
  1452. board.append(dict(t=t, lvl=lvl, ostate=tb.get('设备状态'), estate=tb.get('证据状态'),
  1453. part=tb.get('部件'), steps=steps, stuck=stuck, stuck_kind=stuck_kind,
  1454. next=NEXT.get(stuck, ''), done=sum(1 for x in steps if x['st'] == 'done'),
  1455. rising=rising, hist_loop=(cl if closed else None), srcs=srcs))
  1456. _ORD = {'bad': 0, 'warn': 1, 'note': 2}
  1457. board.sort(key=lambda r: (_ORD.get(r['lvl'], 3), -r['done']))
  1458. _cl_rows = [v for v in _CLOOP.values() if v.get('verdict') == '恢复']
  1459. fus = dict(链盘=dict(rows=board, closed=_cl_rows,
  1460. stuck=dict(collections.Counter(r['stuck'] for r in board if r['stuck'])),
  1461. wo_bound='工单台账 2020~2024-11 (现场未提供 25-26)'),
  1462. 表=ftab.to_dict('records'),
  1463. 矩阵=mrows, 列窗=mheads, kpi=fkpi,
  1464. 事件=ev, 事件边界=ev_bounds,
  1465. 能力=dict(classes=cap['classes'], blind=[str(b) for b in cap['blind']], loop=cap.get('loop', {})),
  1466. 窗=dict(振动=fmeta.get('振动窗', '—'), SCADA=fmeta.get('SCADA窗', ''), 油样=fmeta.get('油样窗', '')),
  1467. 色标=fmeta.get('色标', {}), 数据时点=fmeta.get('数据时点', ''), handoff日期=fmeta.get('handoff日期', ''),
  1468. 证据窗末=fmeta.get('证据窗末', ''),
  1469. 纪律=str(fmeta.get('纪律', '')), 盲区=[str(b) for b in fmeta.get('盲区', [])],
  1470. open_items=[str(o) for o in fmeta.get('open_items', [])], gap=fmeta.get('gap'))
  1471. except Exception as e:
  1472. fus = dict(err=str(e)[:160])
  1473. kpi = dict(报警台=n_alarm_t, 关注台=len(watch), 全场=len(CFG['turbines']))
  1474. return dict(win=win, months=ms, all_months=sorted(_CACHE['tm'].month.unique()),
  1475. systems=systems, sysdist=sysdist, watch=watch, kpi=kpi, rel=rel, fus=fus, faults=faults, control=control, m8=m8, m9=m9,
  1476. # ★2026-09-21: 判级轴现在**按所选窗重算**;首次是后台算,这一份仍是旧口径 ⇒ 如实标 pending
  1477. sysmx_pending=bool(sysmx_pending),
  1478. m9_pending=bool(m9_pending),
  1479. win_pending=bool(sysmx_pending or m9_pending),
  1480. sysmx_span=list(span_of(win)),
  1481. note=('判级轴(变桨/偏航/蓄能/温度)按所选时间窗重算;曲线按所选时间窗重算;'
  1482. '故障统计/五态/温度月轨迹/停机台账=所选窗真窗'
  1483. + ('(⚠ 判级矩阵正在按所选时间窗重算,下面系统卡暂为上一份口径,稍后自动刷新)' if sysmx_pending else '')))
  1484. # ---------- 本体网页面 (ontology 的网页体现: 机制链/决策台/机制库/对象浏览器; 只读, 复用 MCP 纯函数) ----------
  1485. def ontology_view():
  1486. from src.ontology import mcp_server as om
  1487. d = om.db()
  1488. from collections import Counter
  1489. stats = dict(Counter(v['type'] for v in d.values()))
  1490. mechs = [dict(id=o['id'], **{k: o['props'].get(k) for k in ('name', 'status', 'facets', 'source')},
  1491. affects=o.get('links', {}).get('affects', []))
  1492. for o in d.values() if o['type'] == 'FailureMechanism']
  1493. # 29# 机制链 (§2 规划补齐): 机制为中心的分层链
  1494. ev = [dict(id=i, **{k: d[i]['props'].get(k) for k in ('facet', 'claim_window', 'summary', 'source')})
  1495. for i in d['mechanism/WTG29/lub_starvation']['links']['evidencedBy']]
  1496. chain29 = dict(
  1497. mechanism=dict(id='mechanism/WTG29/lub_starvation', **d['mechanism/WTG29/lub_starvation']['props']),
  1498. evidence=ev,
  1499. verdict=dict(id='verdict/WTG29/TCM-16', **{k: d['verdict/WTG29/TCM-16']['props'].get(k)
  1500. for k in ('level', 'title', 'review_tier', 'claim_window')}),
  1501. scenario=dict(id='scenario/WTG29/repair_scheduling_v1',
  1502. **{k: d['scenario/WTG29/repair_scheduling_v1']['props'].get(k)
  1503. for k in ('key_insight', 'recommendation', 'scenarios')}),
  1504. action=dict(id='action/WTG29/repair_pump', **d['action/WTG29/repair_pump']['props']),
  1505. acceptance=dict(id='acceptance/WTG29/state', **d.get('acceptance/WTG29/state', {}).get('props', {})),
  1506. tree=dict(id='decision/WTG29/acceptance_tree', **d['decision/WTG29/acceptance_tree']['props']))
  1507. decisions = [dict(id=o['id'], **{k: o['props'].get(k) for k in ('decided_by', 'when', 'decision', 'status', 'outcome')})
  1508. for o in d.values() if o['id'].startswith('decision/user') or o['id'].startswith('decision/pending')]
  1509. roll = d.get('meta/roll_state', {}).get('props', {})
  1510. constraints = ('C1 证据必带claim_window+观测面+源身份 · C2 判级必带审级+证据链接 · C3 独立源数按观测面去重禁按条数 · '
  1511. 'C4 判读三分(支持|排除|不可判, 不可判≠阴性) · C5 状态措辞窗相关必带窗 · C6 动作必带可逆级(L3须审级)')
  1512. weather = d.get('mwindow/rudong/climatology', {}).get('props', {}).get('monthly', {})
  1513. return dict(stats=stats, constraints=constraints, mechanisms=mechs, chain29=chain29,
  1514. decisions=decisions, roll=roll, weather=weather)
  1515. def ontology_list(otype, turbine, limit=60):
  1516. from src.ontology import mcp_server as om
  1517. return om.ont_query(type=otype or '', turbine=turbine or '', limit=int(limit)).get('data', [])
  1518. def ontology_obj(oid):
  1519. from src.ontology import mcp_server as om
  1520. r = om.ont_get(id=oid)
  1521. lk = om.ont_links(id=oid).get('data', {})
  1522. return dict(obj=r.get('data'), links=lk, err=r.get('error'))
  1523. # ---------- AI 问答 (P2 离线编排的网页入口; dsh + ontology MCP; 串行锁=一次一问, windcms五类错教训) ----------
  1524. CUR_WIN = ['2026年'] # 最近一次 fleet 查询的窗, 供报告助手复用
  1525. ASK = {'state': 'idle', 'q': '', 'model': '', 'answer': '', 'verify': None, 't0': 0, 'err': '',
  1526. 'step': '', 'steps': [], 'review': None, 'xrev': False, 'lang': 'zh'}
  1527. # step/steps: 推理进度 (229s 全程无反馈会被读成"给不出结果"); review: 另一模型的交叉审核结果
  1528. _ASK_STEP_KEY_EN = {
  1529. '启动': 'Starting',
  1530. '预取': 'Prefetch',
  1531. '推理': 'Reasoning',
  1532. '取证': 'Evidence retrieval',
  1533. '成文': 'Drafting',
  1534. '校闸': 'Grounding gate',
  1535. '交叉审核': 'Cross-review',
  1536. '完成': 'Done',
  1537. '失败': 'Failed',
  1538. }
  1539. def _ask_step_text(k, d='', lang='zh'):
  1540. k = '' if k is None else str(k)
  1541. d = '' if d is None else str(d)
  1542. if lang != 'en':
  1543. return f'{k}· {d}' if d else k
  1544. key = _ASK_STEP_KEY_EN.get(k, k)
  1545. detail = d.strip()
  1546. m = re.search(r'第\s*(\d+\s*/\s*\d+)\s*轮', detail)
  1547. if m:
  1548. detail = 'round ' + m.group(1).replace(' ', '')
  1549. else:
  1550. m = re.search(r'(\d{1,2})\s*一站式证据包', detail)
  1551. if m:
  1552. alias = 'focus unit'
  1553. try:
  1554. from src.windscada.deid import unitize_answer
  1555. alias = unitize_answer(f"WTG{int(m.group(1)):02d}")
  1556. except Exception:
  1557. pass
  1558. detail = f'{alias} evidence pack'
  1559. elif '闭环' in detail and '证据' in detail:
  1560. n = ''.join(re.findall(r'\d+', detail))
  1561. detail = f'closure and trend evidence: {n} records' if n else 'closure and trend evidence'
  1562. elif '全场' in detail and '判级' in detail:
  1563. n = ''.join(re.findall(r'\d+', detail))
  1564. detail = f'fleet ratings: {n} units' if n else 'fleet ratings'
  1565. elif '换模型' in detail or '复核' in detail:
  1566. detail = 'reviewing with another model'
  1567. elif re.search(r'[\u3400-\u9fff\u3000-\u303f\uff00-\uffef]', detail):
  1568. detail = ''
  1569. if re.search(r'[\u3400-\u9fff\u3000-\u303f\uff00-\uffef]', key):
  1570. key = 'Working'
  1571. if re.search(r'[\u3400-\u9fff\u3000-\u303f\uff00-\uffef]', detail):
  1572. detail = ''
  1573. return f'{key} · {detail}' if detail else key
  1574. # 模型档位 (2026-08-28 用户裁: 去掉 qwen3:30b-a3b, 用 27B)。
  1575. # 每档标注适用场景与实测耗时 — 让人按"这个问题值不值得等"来选, 而不是凭型号大小猜。
  1576. ASK_MODELS = {'qwen3:8b': 'deploy/dsh/model.qwen8b.patch.yml',
  1577. 'qwen3.8:27b': 'deploy/dsh/model.qwen27b.patch.yml',
  1578. # 2026-09-06 用户裁 "加档": 部署包写的是 "本机 Ollama 的 DeepSeek 模型", 机上实有 r1 14b/32b 却选不到.
  1579. # 本地档一律走 fast_agent 快路径 (Ollama 直连), 不走 dsh, 故 patch 为 None; r1 不支持 tools 时快路径自动退化为预取直答 (见 fast_agent._chat).
  1580. 'deepseek-r1:14b': None, 'deepseek-r1:32b': None,
  1581. 'qwen3:30b-a3b': None, # 2026-09-07: 8B 校闸未过的自动升档目标 (30B 参数 / 3B 激活, 单轮 ~20 s)
  1582. 'deepseek-v4-pro': None}
  1583. # ★这张表原来只有本地模型时代的三档, **不含实际部署的 deepseek-chat / kimi-fast /
  1584. # qwen-plus** —— 于是 MODEL_SPEC[sel.value] 落空, v.warn 为空, 徽章就显示「零外发」。
  1585. # 云端模型对用户说"零外发"是**假保证**, 比时间写错严重得多 (2026-09-01 用户说"太慢"
  1586. # 时顺带查出来的)。补全 + 下面 mdlNote 的兜底改为"按外发处理", 宁可说重不可说轻。
  1587. MODEL_SPEC = {
  1588. 'qwen3:8b': dict(nm='轻量 8B', loc='本地', t='约 30 到 60 秒', use='默认档: 快速简答问答 (用户裁 2026-09-07)'),
  1589. 'qwen3.8:27b': dict(nm='标准 27B', loc='本地', t='约 2 到 4 分钟', use='复杂问题档: 多台横向比对与机制推断 (用户裁 2026-09-07)'),
  1590. 'qwen3:30b-a3b': dict(nm='快 30B (a3b)', loc='本地', t='约 20 到 60 秒', use='8B 校闸未过时的自动升档目标; 也可直接选'),
  1591. 'deepseek-r1:14b': dict(nm='本机 DeepSeek 14B', loc='本地', t='约 1 到 3 分钟', use='交叉审核审核档 (作答请优先 qwen 27B; 用户裁 2026-09-07)'),
  1592. 'deepseek-r1:32b': dict(nm='本机 DeepSeek 32B', loc='本地', t='约 3 到 8 分钟', use='备用: 更强更慢, 引用格式不稳 (约三成概率不出文)'),
  1593. 'deepseek-chat': dict(nm='DeepSeek', loc='云端', t='约 3 到 10 秒', use='默认档, 状态查询与机制推断',
  1594. warn='内容会外发到模型厂商, 场名与台号已去标识化'),
  1595. 'kimi-fast': dict(nm='Kimi', loc='云端', t='约 5 到 20 秒', use='长文归纳',
  1596. warn='内容会外发到模型厂商, 场名与台号已去标识化'),
  1597. 'qwen-plus': dict(nm='通义千问', loc='云端', t='约 5 到 20 秒', use='备用云端档',
  1598. warn='内容会外发到模型厂商, 场名与台号已去标识化'),
  1599. 'deepseek-v4-pro': dict(nm='云端 DS', loc='云端', t='约 2 到 6 分钟', use='复杂推理与交叉审核',
  1600. warn='内容会外发到模型厂商, 场名与台号已去标识化'),
  1601. }
  1602. def _warm(model):
  1603. # 预热+常驻: 空载请求把模型拉进内存并 keep_alive 30m (与 dsh 冷启并行, 消掉 18GB 级冷加载)
  1604. try:
  1605. import urllib.request
  1606. urllib.request.urlopen(urllib.request.Request(
  1607. 'http://127.0.0.1:11434/api/generate', method='POST',
  1608. data=json.dumps(dict(model=model, prompt='', keep_alive='30m')).encode(),
  1609. headers={'Content-Type': 'application/json'}), timeout=120).read()
  1610. except Exception:
  1611. pass
  1612. def _ask_worker(q, model, lang='zh'):
  1613. import subprocess, threading, time as _t
  1614. out = _P.ont() / 'ask_last.out'
  1615. try:
  1616. if MODEL_SPEC.get(model, {}).get('loc') == '本地': # 本地档 (qwen3 / deepseek-r1) 走直连快路径; 原按 'qwen' 前缀硬判, 加 r1 档后改按档表 loc
  1617. # 本地模型走直连快路径 (fast_agent: ollama chat+tools 同进程, 砍 dsh/node/MCP 三层;
  1618. # 实测 8b 116s→19s; 同一套工具与闸, 服务端强制双闸). deepseek 云端仍走 dsh 完整编排.
  1619. threading.Thread(target=_warm, args=(model,), daemon=True).start()
  1620. sys.path.insert(0, str(pathlib.Path('.').resolve()))
  1621. sys.path.insert(0, str(pathlib.Path('scripts').resolve()))
  1622. from src.ontology.fast_agent import answer_escalate as _fast # 2026-09-07: 8B 拦了自动升 a3b 重答 (用户令)
  1623. def _on_step(k, d=''):
  1624. ASK['step'] = _ask_step_text(k, d, lang)
  1625. ASK.setdefault('steps', []).append(ASK['step'])
  1626. def _on_partial(txt):
  1627. # 先快后准 (用户裁 2026-09-07): 30 s 内先出初答 (契约直出), 状态 partial, 前端继续轮询; 27B 出来后整段替换
  1628. ASK.update(answer=txt, state='partial', partial=True)
  1629. r = _fast(q, model, on_step=_on_step, lang=lang, on_partial=_on_partial)
  1630. _esc = f" · 升档 {r['escalated_from']}→{r['model']}" if r.get('escalated_from') else ''
  1631. ans = r['answer'] + f"\n\n({r['gate']} · {r['rounds']}轮 · 直连快路径{_esc})"
  1632. out.write_text(ans, encoding='utf-8') # 答案是中文: 缺 encoding 在中文 Windows 上按 GBK 写, 可能崩
  1633. from ontology_p2_verify import verify as _v0
  1634. ASK.update(answer=ans, verify=_v0(ans, lang=lang), partial=False, refined_from=r.get('escalated_from'))
  1635. # 交叉审核: 换一个本地模型复核同一份证据 (同模型自审 = 自己给自己判卷)
  1636. if ASK.get('xrev'):
  1637. _on_step('交叉审核', '换模型复核中')
  1638. from src.ontology.fast_agent import cross_review as _xr
  1639. # 审核档 = DeepSeek r1:14b (用户裁 2026-09-07); 作答本身是 DeepSeek 时换 qwen 27B 审 —— 始终跨家族 (同家族自审 = 自己判卷)
  1640. _rm = 'deepseek-r1:14b' if not model.startswith('deepseek') else 'qwen3.8:27b'
  1641. ASK['review'] = _xr(q, r['answer'], r.get('facts') or [], model=_rm)
  1642. ASK.update(state='done')
  1643. return
  1644. if False:
  1645. q = '/no_think ' + q # (快路径内已带 /no_think; dsh 路径仅 deepseek)
  1646. cmd = ['deploy/dsh/dsh_run.sh', '--profile', 'headless',
  1647. '--patch', 'deploy/dsh/headless.patch.ontology.yml']
  1648. mp = ASK_MODELS.get(model)
  1649. if mp:
  1650. cmd += ['--patch', mp]
  1651. if lang == 'en':
  1652. cmd.append(q + ' Answer in English only. End every conclusion sentence with [object/id]. Finally run claim_check.')
  1653. else:
  1654. cmd.append(q + ' 每个结论句末尾用[对象id]标注来源。最后调用claim_check。')
  1655. # 2026-09-16: 走 src.proc.run —— 本函数跑在 detail 服务进程里 (无可见控制台), 裸 spawn
  1656. # 会让 Windows 给这条"云端档问答"新建可见控制台窗口 (本地档不走这里)。
  1657. from src import proc as _proc
  1658. r = _proc.run(cmd, capture_output=True, text=True, timeout=600)
  1659. ans = r.stdout.strip() or ('(空输出) stderr: ' + r.stderr[-500:])
  1660. out.write_text(ans, encoding='utf-8')
  1661. sys.path.insert(0, str(pathlib.Path('scripts').resolve()))
  1662. from ontology_p2_verify import verify as _v
  1663. ASK.update(answer=ans, verify=_v(ans, lang=lang), state='done')
  1664. except Exception as e:
  1665. _msg = str(e)
  1666. if 'Connection refused' in _msg or 'Errno 61' in _msg or '11434' in _msg:
  1667. _msg = '本机模型未启动 (Ollama 127.0.0.1:11434 拒连); 离线版不转云, 请先启动 Ollama' # 2026-09-07 真停实测: 原样抛 urlopen 错误给用户
  1668. ASK.update(state='error', err=_msg[:300])
  1669. # ── 报告定制助手 (2026-08-28 用户令: 留一个互动交互窗, 报告按定制, 离线模型支持) ──
  1670. # 两件事: ①听懂"给集团看的月报, 重点讲设备问题"这类话 → 选出该放哪些块
  1671. # ②基于**真实数字**写一段总结叙述 — 为图表补充文字结论, 供汇报使用。
  1672. # 防编造: 数字只许用下面这份摘要里的, 生成后过 grounding 闸复核, 越界即标注。
  1673. RPT_BLOCK_KEYS = ['kpi', 'water', 'emon', 'blame', 'ram', 'parts', 'attrib',
  1674. 'pareto', 'quad', 'frate', 'units', 'chain', 'loop', 'plan', 'caliber']
  1675. RPT_SYS = ('你是风电场报告编排助手。根据用户要求做两件事, 严格按格式输出, 不要解释:\n'
  1676. 'BLOCKS: 用逗号分隔的块名 (只能从给定清单里选, 按报告顺序排)\n'
  1677. 'TEXT: 一段中文总结 (150~300 字, 可分 2~3 句一组)\n'
  1678. '铁律: TEXT 里出现的每个数字都必须来自"可用数据"里, 不许自己算不许估计; '
  1679. '没有依据的话就不写。不要写套话, 不要写"综上所述"。')
  1680. def rpt_compose(need, canvas, model='qwen3:8b'):
  1681. import time as _t
  1682. fl = fleet_view(CUR_WIN[0] if CUR_WIN else '2026年')
  1683. E = (fl.get('control') or {}).get('energy') or {}
  1684. R = fl.get('rel') or {}
  1685. M = ((fl.get('m8') or {}).get('mtbf')) or {}
  1686. F = fl.get('faults') or {}
  1687. BD = ((fl.get('fus') or {}).get('链盘')) or {}
  1688. facts = [
  1689. f"统计期 {fl.get('win')}; 38 台 4.0MW",
  1690. f"上网电量 {E.get('act')} MWh; 理论可发 {E.get('theo')}; 损失 {E.get('loss')} MWh; 损失率 {E.get('loss_pct')}%",
  1691. f"等效满发 {E.get('eflh')} h; 时间可用率 {(fl.get('control') or {}).get('avail')}%",
  1692. f"总停机 {R.get('总停机h')} h = 设备类 {R.get('设备类停机h')} + 外部类 {R.get('外部类停机h')}",
  1693. f"平均停机间隔(MTBO){M.get('MTBF_h')} h; 停机 {M.get('停机事件数')} 起; 单次停机时长(MDT){M.get('MDT_h')} h",
  1694. f"报警 {F.get('total')} 条 / {F.get('n_codes')} 个码",
  1695. f"需行动机组 {len(BD.get('rows') or [])} 台; 卡点 {BD.get('stuck')}",
  1696. # 占比直接给出: 不给模型就得自己算, 算对了也会被接地闸判为"未接地数字"(实测 31.4/46.6)
  1697. f"损失构成: " + '; '.join(
  1698. f"{x['k']} {x['loss']} MWh (占损失 {round(x['loss'] / max(E.get('loss') or 1, 1) * 100, 1)}%)"
  1699. for x in (E.get('items') or [])),
  1700. f"停机 {R.get('总停机h')} h 中设备类占 {round((R.get('设备类停机h') or 0) / max(R.get('总停机h') or 1, 1) * 100, 1)}%",
  1701. f"端到端闭环 {len(BD.get('closed') or [])} 例",
  1702. ]
  1703. prompt = (f"可用数据:\n" + '\n'.join('- ' + x for x in facts) +
  1704. f"\n\n可选块 (只能从这里选): {', '.join(RPT_BLOCK_KEYS)}"
  1705. f"\n块的含义: kpi=一级指标 water=能量瀑布 emon=逐月发电 blame=损失分账 ram=可靠性象限 "
  1706. f"parts=部件停机 attrib=归因覆盖 pareto=故障帕累托 quad=故障四象限 frate=月度故障率 "
  1707. f"units=重点机组 chain=决策链 loop=闭环实证 plan=下期建议 caliber=口径边界"
  1708. f"\n\n当前画布: {', '.join(canvas) if canvas else '(空)'}"
  1709. f"\n用户要求: {need}")
  1710. import urllib.request as _ur # 顶部只 import urllib.parse, urllib.request 需显式导入
  1711. body = dict(model=model, messages=[dict(role='system', content=RPT_SYS),
  1712. dict(role='user', content='/no_think ' + prompt)],
  1713. stream=False, keep_alive='30m', options=dict(num_ctx=8192, num_predict=800))
  1714. t0 = _t.time()
  1715. try:
  1716. req = _ur.Request('http://127.0.0.1:11434/api/chat', data=json.dumps(body).encode(),
  1717. headers={'Content-Type': 'application/json'})
  1718. with _ur.urlopen(req, timeout=300) as r:
  1719. txt = (json.loads(r.read())['message'].get('content') or '')
  1720. except Exception as e:
  1721. return dict(err=f'本地模型不可达: {str(e)[:120]}')
  1722. txt = re.sub(r'<think>.*?</think>', '', txt, flags=re.S).strip()
  1723. mb = re.search(r'BLOCKS\s*[::]\s*(.+)', txt)
  1724. mt = re.search(r'TEXT\s*[::]\s*(.+)', txt, re.S)
  1725. blocks = [b.strip() for b in re.split(r'[,,\s]+', mb.group(1))] if mb else []
  1726. blocks = [b for b in blocks if b in RPT_BLOCK_KEYS]
  1727. narrative = (mt.group(1).strip() if mt else '').strip()
  1728. # 接地闸: 叙述里的数字必须在 facts 里出现 (同 fast_agent 的 grounding 口径)
  1729. bad = []
  1730. try:
  1731. from src.windscada.audit import grounding as _gr
  1732. ok, b = _gr(narrative, '\n'.join(facts))
  1733. bad = (b or {}).get('numbers') or []
  1734. except Exception:
  1735. pass
  1736. return dict(blocks=blocks, narrative=narrative, model=model,
  1737. secs=round(_t.time() - t0, 1), ungrounded=bad, raw=txt[:900])
  1738. ASK_CAT_HINTS = {
  1739. '设备类': '(设备状态类: 从本体的机组档案、对象查询与机制检索三条路取证)',
  1740. '检修排程决策': '(检修排程类: 用 ont_scenario/ont_get 取排期沙盘与验收树; 建议级动作注明待人裁)',
  1741. '损失分析': '(损失分析类: 用 ont_loss/ont_scenario 取证; 电量kWh为主, 金额必带电价区间0.391~0.85口径; 限电类损失≠设备责任)',
  1742. '安全风险预测': '(风险预测类: 用 ont_risk/ont_mechanism_search 取证; 区分升级中/稳定高位/盲区, 未列≠正常, 不可判≠阴性)',
  1743. }
  1744. ASK_CAT_HINTS_EN = {
  1745. '设备类': '(Equipment state: use unit profiles, object queries and mechanism search for evidence)',
  1746. '检修排程决策': '(Maintenance scheduling: use scenarios and acceptance trees; mark advisory actions as pending human decision)',
  1747. '损失分析': '(Loss analysis: use loss accounts and scenarios; energy in kWh/MWh, monetary value must carry the tariff-range caveat)',
  1748. '安全风险预测': '(Risk prediction: distinguish escalating, stable-high and blind-spot cases; no listing is not proof of normality)',
  1749. }
  1750. def ask_start(q, model, cat='', xrev=False, lang='zh', display_q=''):
  1751. import threading, time as _t
  1752. if not q.strip():
  1753. return dict(err='问题为空' if lang != 'en' else 'The question is empty')
  1754. if ASK['state'] in ('running', 'partial'): # partial = 初答已出、27B 复核中, 仍占槽
  1755. return dict(err='已有问题在推理中 (串行锁, 一次一问)', state=ASK['state'])
  1756. if model not in ASK_MODELS:
  1757. return dict(err=f'未知模型 {model}; 可选 {list(ASK_MODELS)}')
  1758. hints = ASK_CAT_HINTS_EN if lang == 'en' else ASK_CAT_HINTS
  1759. q2 = (hints.get(cat, '') + ' ' + q).strip()
  1760. if lang == 'en':
  1761. try:
  1762. from src.windscada.deid import ununitize_question
  1763. q2 = ununitize_question(q2)
  1764. except Exception as e:
  1765. print('[ask] 代号反解失败: %s' % e, file=sys.stderr)
  1766. ASK.update(state='running', q=display_q or q, model=model, answer='', verify=None, err='', t0=_t.time(),
  1767. step='启动' if lang != 'en' else 'Starting', steps=[], review=None, lang=lang)
  1768. ASK['xrev'] = bool(xrev)
  1769. threading.Thread(target=_ask_worker, args=(q2, model, lang), daemon=True).start()
  1770. return dict(state='running', model=model)
  1771. def ask_status():
  1772. import time as _t
  1773. d = dict(ASK)
  1774. d['elapsed'] = int(_t.time() - ASK['t0']) if ASK['t0'] else 0
  1775. return d
  1776. # ---------- 台内即时查询 ----------
  1777. _RAW_CACHE = {}
  1778. def adhoc_query(t, ch, win, agg):
  1779. from src.windscada.data import load_10min, contracted_cols
  1780. _load()
  1781. key = t
  1782. with LOCK:
  1783. if key not in _RAW_CACHE:
  1784. if len(_RAW_CACHE) > 6: _RAW_CACHE.clear()
  1785. try:
  1786. _RAW_CACHE[key] = load_10min(t, CFG, groups=['A.功率', 'A.风况', 'A.转速', 'A.状态', 'B.变桨', 'B.温度NBM', 'B.偏航', 'B.润滑液压'])
  1787. except Exception as e:
  1788. # 2026-09-11 场景① 逮: 原始件不在时这里原来直接把异常抛给 HTTP 处理器 → 前端拿到 500,
  1789. # 面板只显示"服务器错误", 人看不出是"没放数据"还是"程序坏了"。改为**结构化的无数据**回答:
  1790. # 与成功路径同一个形状 (kind/title/unit/months/series/note), 前端按正常渲染即可显示原因。
  1791. src = CFG.get('src_10min')
  1792. return dict(kind='multiline', title=f'{t} {i18n.ch_cn(ch)}', unit='', months=[], series=[],
  1793. err='no_source',
  1794. note=f'无原始件: 本机应在 {src}\\{t}.csv (即 data\\raw\\<场站名称>\\scada_10min\\)'
  1795. f' — 放入后重跑 scripts\\rebuild_from_raw.py --scada')
  1796. d = _RAW_CACHE[key]
  1797. if ch not in d.columns: return dict(err=f'通道 {ch} 不在契约内')
  1798. ms = months_of(win)
  1799. dd = d[d.ts.dt.to_period('M').astype(str).isin(ms)]
  1800. if agg == 'daily':
  1801. g = dd.groupby(dd.ts.dt.date)[ch].median()
  1802. xs = [str(x) for x in g.index]
  1803. else:
  1804. g = dd.groupby(dd.ts.dt.to_period('M').astype(str))[ch].median() if agg == 'med' else \
  1805. dd.groupby([dd.ts.dt.to_period('M').astype(str), dd.ts.dt.date])[ch].sum().groupby(level=0).mean()
  1806. xs = list(g.index)
  1807. agg_cn = {'med': '月中位', 'sum': '月duty(秒/日)', 'daily': '逐日中位'}.get(agg, agg)
  1808. fig = dict(kind='multiline', title=f'{t} {i18n.ch_cn(ch)} ({agg_cn}, {win})', unit='', months=xs,
  1809. series=[dict(name=i18n.ch_cn(ch, with_tag=False), vals=[None if v != v else round(float(v), 3) for v in g])], note='台内即时查询')
  1810. # fleet 带 (若月仓有)
  1811. for store, col in ((_CACHE['tm'], 'med'), (_CACHE['duty'], 'val')):
  1812. if store is None or agg == 'daily': continue
  1813. gg = store[store.channel == ch]
  1814. if len(gg):
  1815. piv = gg.pivot_table(index='month', columns='turbine', values=col).reindex(xs)
  1816. fig = dict(kind='band', title=fig['title'], unit='', months=xs,
  1817. self=[None if v != v else round(float(v), 3) for v in piv.get(t, pd.Series(index=xs, dtype=float))],
  1818. med=[None if v != v else round(float(v), 3) for v in piv.median(axis=1)],
  1819. q1=[None if v != v else round(float(v), 3) for v in piv.quantile(.25, axis=1)],
  1820. q3=[None if v != v else round(float(v), 3) for v in piv.quantile(.75, axis=1)], note='灰带=全场四分位 (月仓)')
  1821. break
  1822. return fig
  1823. def channel_list():
  1824. from src.windscada.data import contract
  1825. c = contract(CFG)
  1826. out = {}
  1827. for g, cols in c['groups'].items():
  1828. out[g] = [dict(code=k, cn=i18n.ch_cn(k, with_tag=False) + ' · ' + k) for k, v in cols.items() if v['liveness'] == '活']
  1829. return out
  1830. # ---------- HTML ----------
  1831. from src.windscada.design import CSS as BASE_CSS
  1832. # 经典页模板(CHART_JS / PAGE_FLEET / PAGE_PROBLEM / PAGE_TURBINE)已随 P5 迁入前端模块;
  1833. # 模块间调用只许走 api.py(见 configs/modules.yaml 的 allow)——这里取一次,名字与迁移前不变。
  1834. from app_frontEnd.app_frontEnd_guanlan.api import classic_pages as _classic
  1835. CHART_JS = _classic.CHART_JS
  1836. PAGE_FLEET = _classic.PAGE_FLEET
  1837. PAGE_PROBLEM = _classic.PAGE_PROBLEM
  1838. PAGE_TURBINE = _classic.PAGE_TURBINE
  1839. class H(BaseHTTPRequestHandler):
  1840. def log_message(self, *a): pass
  1841. def _send(self, body, ctype='text/html; charset=utf-8', code=200):
  1842. if 'text/html' in ctype and not self.path.startswith('/v2'): # 2026-09-06 用户批: /v2 = detail_style_baseline, 不注入 unified-review 样式/导航条 (基线令牌须等于门户)
  1843. from src.windscada.design import GUANLAN_REVIEW_CSS
  1844. text = body.decode('utf-8') if isinstance(body, bytes) else body
  1845. # ★2026-09-22 用户测试 4 实逮: 这里原来写死 `http://127.0.0.1:18084/观澜_如东样板_门户_单文件.html`
  1846. # —— 用户从别的机器点「返回观澜门户」等于访问**自己机器的 loopback**,且那文件在服务端并不存在
  1847. # ⇒ 网关回 {"err":"not found","path":"/观澜_…html"}。
  1848. # 修法: 用根相对 `/`(门户就在网关根)+ `data-abs="1"` 让**网关不要按组件前缀改写它**
  1849. # (不加标记会被改写成 /detail/,永远回不到门户 —— 见 guanlan_gateway.rewrite 的说明)。
  1850. text = text.replace('</header>', '</header><nav class="action-path" aria-label="页面操作">'
  1851. '<a href="/" data-abs="1" target="_top">← 返回观澜门户</a>'
  1852. '<a href="/detail/v2" data-abs="1">分析总览 / 重新选择</a></nav>', 1)
  1853. if '</head>' in text:
  1854. body = text.replace('</head>', '<style id="guanlan-unified-review-v1">' + GUANLAN_REVIEW_CSS + '</style></head>', 1)
  1855. b = body.encode('utf-8') if isinstance(body, str) else body
  1856. self.send_response(code)
  1857. self.send_header('Content-Type', ctype)
  1858. self.send_header('Content-Length', str(len(b)))
  1859. self.end_headers()
  1860. self.wfile.write(b)
  1861. def do_POST(self):
  1862. try:
  1863. u = urllib.parse.urlparse(self.path)
  1864. n = int(self.headers.get('Content-Length', 0))
  1865. body = json.loads(self.rfile.read(n) or b'{}')
  1866. if u.path == '/api/ask':
  1867. self._send(_jdump(ask_start(body.get('q', ''), body.get('model', 'qwen3:8b'), body.get('cat', ''),
  1868. body.get('xrev', False), body.get('lang', 'zh'), body.get('display_q', '')),
  1869. ensure_ascii=False), 'application/json; charset=utf-8')
  1870. else:
  1871. self._send('404', code=404)
  1872. except Exception as e:
  1873. self._send(_jdump(dict(err=str(e)[:300]), ensure_ascii=False), 'application/json; charset=utf-8', 500)
  1874. def do_GET(self):
  1875. try:
  1876. u = urllib.parse.urlparse(self.path)
  1877. # encoding='utf-8' 显式给定 (默认 latin-1 会把中文参数解成乱码 — "齿轮箱"→"齿轮箱",
  1878. # 浏览器 encodeURIComponent 路径侥幸正常, curl/直传中文即崩; 2026-08-28 四链接口实逮)
  1879. q = dict(urllib.parse.parse_qsl(u.query, encoding='utf-8', errors='replace'))
  1880. _REQ.lang = q.get('lang', '') # 本请求的语言, _jdump 出口据此投影文本层
  1881. _REQ.keep_orig = q.get('keep') == '1' # v2: 译不出留原文, 不写 [translation pending]
  1882. win = q.get('win', '2026年')
  1883. if u.path == '/':
  1884. _pg = PAGE_FLEET.replace('%CSS%', BASE_CSS).replace('%JS%', CHART_JS)
  1885. self._send(_en_page(_pg))
  1886. elif u.path == '/v2':
  1887. # 新前端 (2026-09-03 用户裁"完全重新做一套, 在原来的基础上"): 同一套 /api/*,
  1888. # 表现层双语原生 —— 不经 _en_page/ui_en 那三层事后翻译。迁移期与 / 并存,
  1889. # 页签逐个搬完后再切默认。
  1890. from src.windscada.ui import build as _v2
  1891. self._send(_v2.render('en' if q.get('lang') == 'en' else 'zh'))
  1892. elif u.path == '/v2/snapshot':
  1893. # 新版单文件快照 (静态站用): 接口响应烤进页面, 拦截器实时优先快照兜底。
  1894. # 取数函数与在线 API 同一实现; 英文走同一套 _redact → _proj 投影 (keep=1 留原文)。
  1895. from src.windscada.ui import snapshot as _snap
  1896. from src.ontology import mcp_server as _omc
  1897. import json as _j3, pathlib as _p3
  1898. def _codes():
  1899. db = _j3.loads(_p3.Path(_omc.OBJ_PATH).read_text(encoding='utf-8'))
  1900. return sorted(k.split('/')[-1] for k in db if k.startswith('alarmcode/'))
  1901. _ctx = dict(fleet_view=fleet_view, curve_view=curve_view, vibcms_results=vibcms_results,
  1902. survey=__import__('src.ontology.maintenance', fromlist=['survey']).survey,
  1903. ont=_omc, alarm_codes=_codes, turbine_problems=turbine_problems,
  1904. single_problem=single_problem, _redact=_redact, _jclean=_jclean, _proj=_proj, _REQ=_REQ)
  1905. _v2_lang = 'en' if q.get('lang') == 'en' else 'zh'
  1906. _html, _rep = _snap.bake(_v2_lang, _ctx,
  1907. peer=q.get('peer', ''), noindex=q.get('noindex') == '1',
  1908. slim=(q.get('slim') != '0' if _v2_lang == 'en' else False))
  1909. print('[v2 snapshot] %s' % _j3.dumps(_rep, ensure_ascii=False)[:600], flush=True)
  1910. self._send(_html)
  1911. elif u.path.startswith('/turbine/'):
  1912. t = u.path.rsplit('/', 1)[1]
  1913. _pg = PAGE_TURBINE.replace('%CSS%', BASE_CSS).replace('%JS%', CHART_JS).replace('%T%', t)
  1914. self._send(_en_page(_pg))
  1915. elif u.path == '/api/maint_survey':
  1916. try:
  1917. from src.ontology.maintenance import survey as _sv
  1918. r = dict(rows=_sv())
  1919. except Exception as e:
  1920. r = dict(rows=[], error=f'维护面盘点失败: {e}'[:200])
  1921. self._send(_jdump(r, ensure_ascii=False), 'application/json; charset=utf-8')
  1922. elif u.path == '/api/ont_chain':
  1923. try:
  1924. from src.ontology import mcp_server as _om
  1925. kind = q.get('kind', 'fault')
  1926. fn = {'fault': lambda: _om.ont_chain_fault(code=q.get('code', '')),
  1927. 'prev': lambda: _om.ont_chain_preventive(system=q.get('system', '')),
  1928. 'tree': lambda: _om.ont_fault_tree(system=q.get('system', '')),
  1929. 'plan': lambda: _om.ont_maint_plan(top=q.get('top', 10)),
  1930. 'doc': lambda: _om.ont_doc_search(keyword=q.get('keyword', ''), kind=q.get('dkind', ''))}
  1931. r = fn[kind]() if kind in fn else dict(text=f'未知链 {kind}', error=True)
  1932. except Exception as e:
  1933. r = dict(text=f'本体查询失败: {e}'[:200], error=True)
  1934. self._send(_jdump(r, ensure_ascii=False), 'application/json; charset=utf-8')
  1935. elif u.path == '/api/maint_std':
  1936. # 运维标准四判据 + 成熟度评分 (docs/如东运维标准化设计.md)
  1937. try:
  1938. import json as _j, pathlib as _p
  1939. from src.windscada.maint import standards as _st
  1940. _db = _j.loads(_P.objects_json()
  1941. .read_text(encoding='utf-8'))
  1942. _wo = [(v.get('props') or {}) for k, v in _db.items()
  1943. if k.startswith('workorder/')]
  1944. _r = _st.run_all(_wo)
  1945. _r['score'] = _st.score(_wo, _r)
  1946. except Exception as e:
  1947. _r = dict(err=f'运维标准判据失败: {e}'[:200])
  1948. self._send(_jdump(_r, ensure_ascii=False), 'application/json; charset=utf-8')
  1949. elif u.path == '/api/maint_framework':
  1950. # 运维管理体系 —— 与快照预烤**同一实现** (maint_framework_view)
  1951. try:
  1952. _r = maint_framework_view()
  1953. except Exception as e:
  1954. _r = dict(err=f'运维体系读取失败: {e}'[:200])
  1955. self._send(_jdump(_r, ensure_ascii=False), 'application/json; charset=utf-8')
  1956. elif u.path == '/api/dq_findings':
  1957. # 数据可信度发现 —— 与快照预烤**同一实现** (dq_findings_view)
  1958. try:
  1959. _r = dq_findings_view()
  1960. except Exception as e:
  1961. _r = dict(rows=[], err=f'数据可信度读取失败: {e}'[:200])
  1962. self._send(_jdump(_r, ensure_ascii=False), 'application/json; charset=utf-8')
  1963. elif u.path == '/api/vibcms':
  1964. self._send(_jdump(vibcms_results(), ensure_ascii=False), 'application/json; charset=utf-8')
  1965. elif u.path == '/api/curves':
  1966. # ★2026-09-21: 接 win(按所选时间窗重算;未命中就起后台线程 + 回 building,页面轮询)
  1967. self._send(_jdump(curve_view(q.get('win') or '2026年'), ensure_ascii=False),
  1968. 'application/json; charset=utf-8')
  1969. elif u.path == '/api/fleet':
  1970. CUR_WIN[0] = win
  1971. self._send(_jdump(fleet_view(win), ensure_ascii=False), 'application/json; charset=utf-8')
  1972. elif u.path.startswith('/problem/'):
  1973. import urllib.parse as _up
  1974. _, _, t, sysn = u.path.split('/', 3)
  1975. _raw = _up.unquote(sysn)
  1976. _zh = sys_norm(_raw)
  1977. # ★2026-09-22 用户测试 3/5: 认得出就换成规范中文名(页面/接口都按中文名取数),
  1978. # 认不出就把原文交给页面 —— 页面据此**如实说"系统标识无法识别"**,
  1979. # 绝不写成"该系统在判级时间窗内无非常态"(那对报警机组是假陈述)。
  1980. self._send(PAGE_PROBLEM.replace('%CSS%', BASE_CSS).replace('%JS%', CHART_JS)
  1981. .replace('%T%', t).replace('%SYS%', _zh or _raw)
  1982. .replace('%SYSRAW%', _raw).replace('%SYSOK%', '1' if _zh else '0'))
  1983. elif u.path == '/api/problem':
  1984. _zh = sys_norm(q['sys']) or q['sys']
  1985. self._send(_jdump(dict(problems=single_problem(q['t'], _zh, win), sys=_zh,
  1986. sys_ok=bool(sys_norm(q['sys']))), ensure_ascii=False),
  1987. 'application/json; charset=utf-8')
  1988. elif u.path == '/api/turbine':
  1989. self._send(_jdump(dict(problems=turbine_problems(q['t'], win)), ensure_ascii=False), 'application/json; charset=utf-8')
  1990. elif u.path == '/api/query':
  1991. agg = q.get('agg', 'med')
  1992. self._send(json.dumps(adhoc_query(q['t'], q['ch'], win, 'sum' if agg == 'dutyd' else agg), ensure_ascii=False), 'application/json; charset=utf-8')
  1993. elif u.path == '/api/snapshot':
  1994. _lang = urllib.parse.parse_qs(u.query).get('lang', [''])[0]
  1995. # 单文件快照 (2026-08-28 用户: 我怎么发给别人)。
  1996. # 把各窗 fleet / curves / ontology / 维护盘点全部内联进 HTML, 并在页内改写
  1997. # window.fetch 直接返回内联数据 —— 前端代码一行不改, 收件人双击即开, 不需要
  1998. # 装 Python、不需要连内网、断网也能看。数据按当前受众设置脱敏后再内联。
  1999. try:
  2000. _wins = ['2026年', '2026H1', '2025H2', '全程']
  2001. snap, _skip = {}, []
  2002. def _nk(path, **kw):
  2003. """快照 key 归一化: 参数**解码后**排序拼接。
  2004. 必须与页内 window.__NK 同语义 —— 直接用 URL 字符串做 key 会被
  2005. encodeURIComponent 与 urllib quote 的编码差异 (! * ' ( ) 一编一不编)
  2006. 和参数顺序打掉, 而这种失效是**静默**的: 命中不了就退到
  2007. "快照版不含此功能", 看起来像功能没做, 不像 key 错了。
  2008. """
  2009. a = sorted('%s=%s' % (k, v) for k, v in kw.items() if v not in (None, ''))
  2010. return path + ('?' + '&'.join(a) if a else '')
  2011. def _bake(key, fn):
  2012. """烤一个 API 响应。失败**记名**不静默 —— 守护步骤失败必须响亮。
  2013. ★2026-08-31: 此前直接 _redact(_jclean(...)), **绕过 _jdump 里的语言投影**,
  2014. 于是 ?lang=en 的快照内联的仍是中文数据 (实测 57 万→11 万中文只是
  2015. 页面文案变了, 数据没变)。这里补上同一套投影。
  2016. """
  2017. try:
  2018. # ★顺序是 **_redact 先, _proj 后**, 不要"为了与在线同序"去调它:
  2019. # 曲线图注等模板的 key 是照**脱敏后**的文本登记的 (含 ▪ 掩码),
  2020. # 先投影的话字符串里还是真实数字, 模板全部对不上, 整批退回
  2021. # [translation pending]。我按"同序"改过一次, 当场被 parity 闸逮住。
  2022. # 在线路径没有这个问题是因为脱敏在页面层统一做, 不在数据层。
  2023. _v = _redact(_jclean(fn()))
  2024. if _lang and _lang != 'zh':
  2025. _v = _proj(_v, _lang)
  2026. snap[key] = _v
  2027. return True
  2028. except Exception as e:
  2029. _skip.append('%s: %s' % (key, str(e)[:60]))
  2030. return False
  2031. for _w in _wins:
  2032. _bake(_nk('/api/fleet', win=_w), lambda w=_w: fleet_view(w))
  2033. _bake('/api/curves', curve_view)
  2034. _bake('/api/ontology', ontology_view)
  2035. # ★后加的面板接口也要预烤 —— 漏了会在离线快照里空掉。
  2036. # 新增接口时必须回到这里补一行, 否则快照静默少一块。
  2037. def _maint_std():
  2038. import json as _j2, pathlib as _p2
  2039. from src.windscada.maint import standards as _st2
  2040. _db2 = _j2.loads(_P.objects_json()
  2041. .read_text(encoding='utf-8'))
  2042. _wo2 = [(v.get('props') or {}) for k, v in _db2.items()
  2043. if k.startswith('workorder/')]
  2044. _r2 = _st2.run_all(_wo2)
  2045. _r2['score'] = _st2.score(_wo2, _r2)
  2046. return _r2
  2047. _maint_fw = maint_framework_view # ★与在线 API 同一实现
  2048. _dq = dq_findings_view # ★与在线 API 同一实现
  2049. _bake('/api/maint_std', _maint_std)
  2050. _bake('/api/maint_framework', _maint_fw)
  2051. _bake('/api/dq_findings', _dq)
  2052. _bake('/api/maint_survey',
  2053. lambda: dict(rows=__import__('src.ontology.maintenance',
  2054. fromlist=['survey']).survey()))
  2055. # ── 扩展预烤 (2026-08-30 用户"完整展示") ───────────────────────
  2056. # 原先只烤上面 4 个, 另外 7 个接口在离线版里是死的: 部件问题下钻、
  2057. # 本体对象/清单/链、振动详情点开全是"快照版不含此功能" —— 页面能渲染
  2058. # 掩盖了功能不可用, 我此前按"十页可交互"汇报是说过头了。
  2059. # ask / rpt_compose 要现跑模型, 不烤 (页面另给在线提示), 其余全量枚举。
  2060. _bake('/api/vibcms', vibcms_results)
  2061. # 离线桩: 快照里没有服务端跑模型, 所以这条**有意**与在线不同。
  2062. # note 双写 —— 靠翻译层会得到 [translation pending], 而这是给读者看的一句话。
  2063. _bake('/api/ask_status', lambda: dict(
  2064. offline=True, models=[],
  2065. note_zh='离线版: 模型问答需连接系统',
  2066. note_en='Offline copy: the model Q&A needs a live connection to the system'))
  2067. # 部件问题: (机组, 系统) 从各窗 systems 展开, 每个窗各烤一份
  2068. _pairs = set()
  2069. for _w in _wins:
  2070. _fv = snap.get(_nk('/api/fleet', win=_w)) or {}
  2071. for _sysn, _rows in (_fv.get('systems') or {}).items():
  2072. for _r in (_rows or []):
  2073. if isinstance(_r, dict) and _r.get('t'):
  2074. _pairs.add((_r['t'], _sysn))
  2075. # ★系统名两种形态都要烤 (2026-09-01 实逮"点开是空浮层")。
  2076. # systems 的分组键带**等级后缀** (`塔架与基础·良好`), 而页面上
  2077. # `/problem/<t>/<sys>` 链接发出的是**裸系统名** (`塔架与基础`) ——
  2078. # 于是 key 差一截, 取不到数, 浮层开了但里面是"本时间窗无记录"。
  2079. # **这种失效最难发现: 浮层确实弹出来了, 看着像"这台这系统没数据",
  2080. # 而不像"key 错了"。** 两种形态都烤成同一份数据即可。
  2081. for _t, _sy in sorted(_pairs):
  2082. _forms = {_sy}
  2083. if '·' in _sy:
  2084. _forms.add(_sy.split('·')[0])
  2085. for _w in _wins:
  2086. for _f in sorted(_forms):
  2087. # ★每种形态烤**它自己**的数据。上一版这里捕获的是带等级的 _sy,
  2088. # 于是裸名 key 拿到的是带等级那份 —— 实测差出 0 vs 1 条,
  2089. # **被对拍闸逮到**。lambda 默认参数必须捕获循环变量 _f 本身。
  2090. _bake(_nk('/api/problem', t=_t, sys=_f, win=_w),
  2091. lambda t=_t, sy=_f, w=_w: dict(problems=single_problem(t, sy, w)))
  2092. # 整机视图: 页面上机组号链到 /turbine/<t>, 弹层要用
  2093. for _t in sorted({t for t, _ in _pairs}):
  2094. for _w in _wins:
  2095. _bake(_nk('/api/turbine', t=_t, win=_w),
  2096. lambda t=_t, w=_w: dict(problems=turbine_problems(t, w)))
  2097. # 本体清单 + 每个对象详情 (清单里出现的 id 才可能被点开)
  2098. _oids = set()
  2099. for _ty in sorted((ontology_view().get('stats') or {})):
  2100. _k = _nk('/api/ont_list', type=_ty, limit=60)
  2101. if _bake(_k, lambda ty=_ty: dict(rows=ontology_list(ty, '', 60))):
  2102. for _row in (snap[_k].get('rows') or []):
  2103. if isinstance(_row, dict) and _row.get('id'):
  2104. _oids.add(_row['id'])
  2105. for _oid in sorted(_oids):
  2106. _bake(_nk('/api/ont_obj', id=_oid), lambda i=_oid: ontology_obj(i))
  2107. # 本体链: kind 五种, 参数取页面下拉的实际取值域
  2108. try:
  2109. from src.ontology import mcp_server as _omc
  2110. _bake(_nk('/api/ont_chain', kind='plan', top='12'),
  2111. lambda: _omc.ont_maint_plan(top=12))
  2112. _syslist = sorted({sy for _, sy in _pairs}) or ['齿轮箱']
  2113. for _sy in _syslist:
  2114. _bake(_nk('/api/ont_chain', kind='prev', system=_sy),
  2115. lambda sy=_sy: _omc.ont_chain_preventive(system=sy))
  2116. _bake(_nk('/api/ont_chain', kind='tree', system=_sy),
  2117. lambda sy=_sy: _omc.ont_fault_tree(system=sy))
  2118. # ★原来码源是**另一处烤好的** /api/ont_list?limit=60 的行, 于是全场 555 个码
  2119. # 只进了按码号排序的前 60 个; 而页面上"点一下就查"的示例码
  2120. # (3225/64000/64001/4110/13121) 恰好都不在那 60 个里
  2121. # ⇒ **页面在推荐自己答不上的码**, 点下去整条抢修链没反应。
  2122. # "排序后取前 N"的裁剪最容易把真正要用的挡在外面。
  2123. # ★改为**直接从本体读**, 不再依赖别处烤没烤、烤了多少 ——
  2124. # 跨模块借数据是这次截断的根因, 借来的上限管不着。
  2125. # (中途踩坑: 我先写成 _omc.ont_list(...), 而 mcp_server 根本没有这个函数,
  2126. # AttributeError 被外层 try 吞掉 ⇒ 整段 fault 一条都没烤, 还静默。)
  2127. import json as _js
  2128. _db = _js.loads(pathlib.Path(_omc.OBJ_PATH).read_text(encoding='utf-8'))
  2129. _codes = sorted(k.split('/')[-1] for k in _db
  2130. if k.startswith('alarmcode/'))
  2131. for _c in [c for c in _codes if c]:
  2132. _bake(_nk('/api/ont_chain', kind='fault', code=_c),
  2133. lambda c=_c: _omc.ont_chain_fault(code=c))
  2134. except Exception as _e:
  2135. _skip.append('ont_chain: %s' % str(_e)[:60])
  2136. if _skip:
  2137. print('[snapshot] 未烤成 %d 项: %s' % (len(_skip), '; '.join(_skip[:6])),
  2138. file=sys.stderr)
  2139. print('[snapshot] 预烤 %d 键 (%s)' % (
  2140. len(snap), ', '.join(sorted({k.split('?')[0] for k in snap}))), file=sys.stderr)
  2141. if _lang and _lang != 'zh':
  2142. # ★key 统一在**全部烤完之后**加 lang: 英文页的 fetch 会自动追加
  2143. # ?lang=en, 预烤 key 不带就永远对不上 (miss 是静默的, 退到 err 提示,
  2144. # 看着像功能没做)。
  2145. # ★不能在 _bake 内部改 —— 后续代码还要用原 key 从 snap 取值,
  2146. # 改了就 KeyError (本次实测踩到)。
  2147. # ★存双 key: 带 lang 与不带。页面首帧的 fetch 可能早于 fetch patch
  2148. # 生效, 那一发不带 lang —— 只存带 lang 的会 miss 一次 (实测首屏 fleet)。
  2149. #
  2150. # ★★ 加 lang 必须走**同一个归一化器**, 不能字符串拼接 (2026-09-01 实逮)。
  2151. # __NK 把参数**按字母排序**, 所以正确形态是 `?lang=en&sys=…&t=…&win=…`;
  2152. # 而拼接得到的是 `?sys=…&t=…&win=…&lang=en` —— 差一个顺序, 全部 miss。
  2153. # 后果极隐蔽: 十个标签页照常显示 (它们命中的是**不带 lang** 的那份),
  2154. # 只有**点开下钻**才走带 lang 的 key ⇒ 浮层弹出来却是"本时间窗无记录",
  2155. # 看着像"这台这系统没数据", 不像 key 错了。用户报的"有些点不开"就是它。
  2156. def _with_lang(k):
  2157. base, _, q = k.partition('?')
  2158. kv = dict(x.split('=', 1) for x in q.split('&') if '=' in x) if q else {}
  2159. kv['lang'] = _lang
  2160. return base + '?' + '&'.join(sorted('%s=%s' % (a, b) for a, b in kv.items()))
  2161. snap = dict(snap, **{_with_lang(k): v for k, v in snap.items()})
  2162. page = PAGE_FLEET.replace('%CSS%', BASE_CSS).replace('%JS%', CHART_JS)
  2163. if _lang and _lang != 'zh':
  2164. # 快照页面本身也走英文化 (文案表 + 显示层脚本), 并改 lang 属性 ——
  2165. # 此前只投影了 API 路径, 快照页漏了, 导出的 index.html 一直是 lang="zh"
  2166. page = _en_page(page).replace('<html lang="zh"', '<html lang="%s"' % _lang, 1)
  2167. inject = ('<script>window.__SNAP__='
  2168. + json.dumps(snap, ensure_ascii=False, allow_nan=False) + ';\n'
  2169. '/* 快照模式: 拦截 fetch 从内联数据取。key 走归一化(参数解码后排序)而非 */'
  2170. '/* 原样 URL —— 编码差异(encodeURIComponent 不编 !*\'() 而 quote 编)与参数 */'
  2171. '/* 顺序都会让字符串比对落空, 且落空是静默的: 退到 err 提示, 看着像功能没做。 */'
  2172. 'window.__NK=function(u){var s=String(u),i=s.indexOf("?");if(i<0)return s;'
  2173. 'var p=s.slice(0,i),a=s.slice(i+1).split("&").filter(Boolean).map(function(kv){'
  2174. 'var j=kv.indexOf("=");if(j<0)return kv;var k=kv.slice(0,j),v=kv.slice(j+1);'
  2175. 'try{v=decodeURIComponent(v.replace(/\\+/g," "))}catch(e){}'
  2176. 'return v===""?null:k+"="+v}).filter(function(x){return x!==null}).sort();'
  2177. 'return p+(a.length?"?"+a.join("&"):"")};'
  2178. 'window.__HIT={ok:0,miss:[]};window.__RF=window.fetch;'
  2179. # ★没烤的窗按钮直接摘掉: 让用户点一个必然空白的选项, 比少给几个窗更糟。
  2180. # 默认窗取**已烤窗里最新的那个**, 供上面的回落候选用。
  2181. 'window.__WINS=(function(){var s={};Object.keys(window.__SNAP__).forEach('
  2182. 'function(k){var m=/[?&]win=([^&]*)/.exec(k);if(m)s[decodeURIComponent(m[1])]=1});'
  2183. 'return Object.keys(s)})();'
  2184. 'window.__DEFWIN=window.__WINS[0]||"";'
  2185. 'window.__WINLABEL=' + json.dumps(
  2186. 'Latest completed build' if _lang == 'en'
  2187. else '最新一次完整分析') + ';'
  2188. # ★快照里**去掉窗选择器**, 换成固定标签 (2026-09-01 用户令
  2189. # "设定不随时间窗动, 就显示最新的")。
  2190. # 理由: 选择器给约 20 个窗而只烤了 4 个, 选中其余任何一个都空白;
  2191. # 而这 4 个窗的内容**完全一样** (fus 40 行 / 需跟踪 76 台次 逐窗相同),
  2192. # 所以留着选择器只会让人以为能切、切了还空。
  2193. # ★先前那版是"摘掉没烤的 option" —— 摘过头了: 选项显示文本已英文化
  2194. # ("2026 full year"), 而 __WINS 是中文键, 匹配不上就把有效的也摘了。
  2195. # **文本匹配在译过的 UI 上天然不可靠**, 所以改为整体替换而非逐项筛。
  2196. 'window.__PINWIN=function(){try{'
  2197. 'var sel=document.querySelector("select#win");if(!sel)return;'
  2198. 'var lab=document.createElement("span");lab.id="winpin";'
  2199. 'lab.style.cssText="font-size:12.5px;color:var(--muted);padding:5px 9px;'
  2200. 'border:1px solid var(--line);border-radius:7px;background:var(--panel)";'
  2201. 'lab.textContent=window.__WINLABEL||"Latest build";'
  2202. 'lab.title="This snapshot is fixed to the latest completed build; '
  2203. 'the panels below do not vary by window.";'
  2204. 'sel.parentNode.replaceChild(lab,sel)}catch(e){}};'
  2205. 'try{new MutationObserver(function(){window.__PINWIN()})'
  2206. '.observe(document.documentElement,{childList:true,subtree:true})}catch(e){}'
  2207. 'document.addEventListener("DOMContentLoaded",function(){'
  2208. 'setTimeout(window.__PINWIN,300)});'
  2209. # 问答接口不能被离线拦截: 云上 nginx /api/ 后面有真服务在跑,
  2210. # 其余接口才走预烤。放行名单写死, 免得日后新增接口误穿。
  2211. 'window.fetch=function(u,o){var k=window.__NK(u);'
  2212. 'if(/^\\/api\\/(ask|ask_status|ask_models|rpt_compose|rpt_export)\\b/.test(String(u)))return window.__RF(u,o);'
  2213. # ★第 4 个候选: 把 win 换成默认窗 (2026-09-01 用户令"不随时间窗动, 就显示最新的")。
  2214. # 窗选择器给了约 20 个选项而快照只烤了 4 个 —— 选中其余任何一个
  2215. # 都会整页空白 (表现为"该系统全场无异常"、振动页空), 看着像没数据,
  2216. # 实则是没烤。这些面板的内容本就不随窗变 (四个窗 fus 完全一样),
  2217. # 所以回落到默认窗是**正确答案**而不是遮丑。
  2218. # ★同 src/windscada/ui/snapshot.py 的坑: 这里不能 encodeURIComponent —— k 与
  2219. # __SNAP__ 键都是原文(未编码)中文, 编码后拼出的键谁都不匹配, 兜底静默失效
  2220. # (2026-09-03 实逮: 换窗点浮层报"本时间窗无记录", 默认窗其实有数据)。
  2221. 'var c=[k,k.replace(/&?win=[^&]*/,"").replace(/\\?$/,""),String(u).split("?")[0],'
  2222. 'k.replace(/win=[^&]*/,"win="+(window.__DEFWIN||""))];'
  2223. 'for(var i=0;i<c.length;i++){var d=window.__SNAP__[c[i]];'
  2224. 'if(d!==undefined){window.__HIT.ok++;return Promise.resolve({ok:true,'
  2225. 'json:function(){return Promise.resolve(d)}})}}'
  2226. 'if(window.__HIT.miss.indexOf(k)<0)window.__HIT.miss.push(k);'
  2227. 'return Promise.resolve({ok:true,json:function(){return Promise.resolve('
  2228. '{err:"离线版不含此功能 (模型问答与报告生成需连接系统)"})}})};'
  2229. '</script>')
  2230. # ★英文版改走**源码级生成**, 不再注入运行时翻译 hook
  2231. # (2026-08-31 用户令: 不要打补丁式翻译, 整个系统重新生成一遍)。
  2232. # hook 的本质是给中文系统套一层英文 —— 只译写入 DOM 的那一刻,
  2233. # JS 逻辑与 JSON 键仍是中文, 长尾补不完且源码里满是中文。
  2234. # 改为在页面拼好后整页替换 (见本函数末尾与 i18n_en.to_english_source)。
  2235. # ★中文版完全不受影响: 不带 lang=en 时这条分支不执行, 国内照常用。
  2236. # ★ 单文件版的下钻: 页面上 74 个链接指向 /problem/<t>/<sys> 与 /turbine/<t>,
  2237. # 那是服务端独立页面 —— 快照里点了就是死链 (404/白页)。此前"十页可交互"
  2238. # 的说法没覆盖这一层。改法: 拦截这两类链接, 用已预烤的同一份数据在页内
  2239. # 开浮层, 图表复用共享的 drawFig, 与在线版看到的是同一批图。
  2240. # ★注入块里的文案要**自己**按语言选: 它加在 _en_page 之后,
  2241. # 英文化那一遍已经跑完了, 这里的中文字面量永远不会被翻译。
  2242. # 实逮: 英文页浮层里弹出"本时间窗无记录"。
  2243. # ('状态'/'依据' 是 JSON 键, '年' 是窗名数据值 —— 都不能动。)
  2244. inject += ('<script>(function(){var M=null;'
  2245. 'function esc(s){return String(s==null?"":s).replace(/[&<>]/g,function(c){'
  2246. 'return {"&":"&amp;","<":"&lt;",">":"&gt;"}[c]})}'
  2247. 'window.__CLOSEM=function(){if(M)M.style.display="none"};'
  2248. 'function box(){if(M)return M;M=document.createElement("div");'
  2249. 'M.style.cssText="position:fixed;left:0;top:0;right:0;bottom:0;z-index:9999;'
  2250. 'background:rgba(0,0,0,.6);overflow:auto";'
  2251. 'M.onclick=function(e){if(e.target===M)window.__CLOSEM()};'
  2252. 'document.body.appendChild(M);return M}'
  2253. 'async function open2(t,sys){'
  2254. 'var w=(window.D&&window.D.win)||"2026年";'
  2255. 'var u=sys?("/api/problem?t="+t+"&sys="+encodeURIComponent(sys)+"&win="+encodeURIComponent(w))'
  2256. ':("/api/turbine?t="+t+"&win="+encodeURIComponent(w));'
  2257. 'var d={};try{d=await (await fetch(u)).json()}catch(e){}'
  2258. 'var ps=(d&&d.problems)||[];'
  2259. 'var h=\'<div style="max-width:1120px;margin:26px auto;background:var(--card,#111a24);\''
  2260. '+\'border:1px solid var(--line,#22303d);border-radius:10px;padding:18px 20px">\''
  2261. '+\'<div style="display:flex;justify-content:space-between;align-items:flex-start;gap:12px">\''
  2262. '+"<h3 style=\'margin:0 0 4px\'>"+esc(t)+(sys?" · "+esc(sys):"")+"</h3>"'
  2263. '+\'<button onclick="window.__CLOSEM()" style="background:none;border:1px solid \''
  2264. '+\'var(--line,#22303d);color:inherit;border-radius:6px;padding:2px 9px;cursor:pointer">×</button></div>\';'
  2265. 'if(!ps.length)h+=\'<p class="small">' + ('No records in this window' if _lang == 'en' else '本时间窗无记录') + '</p>\';'
  2266. 'ps.forEach(function(p){h+=\'<div class="card mb"><b>\'+esc(p.sys)+\'</b> · \'+esc(p["状态"])'
  2267. '+\'<p class="small">\'+esc(p["依据"])+\'</p>\';'
  2268. 'if(p.figs&&p.figs.length&&typeof drawFig==="function"){h+=\'<div class="grid">\';'
  2269. 'p.figs.forEach(function(f){try{h+=drawFig(f)}catch(e){}});h+="</div>"}'
  2270. 'h+="</div>"});'
  2271. 'h+="</div>";var b=box();b.innerHTML=h;b.style.display="block";b.scrollTop=0}'
  2272. 'document.addEventListener("click",function(e){'
  2273. 'var a=e.target&&e.target.closest?e.target.closest(\'a[href^="/problem/"],a[href^="/turbine/"]\'):null;'
  2274. 'if(!a)return;e.preventDefault();e.stopPropagation();'
  2275. 'var q=a.getAttribute("href").split("/").filter(Boolean);'
  2276. 'open2(decodeURIComponent(q[1]||""),q[0]==="problem"?decodeURIComponent(q[2]||""):"")},true);'
  2277. 'document.addEventListener("keydown",function(e){if(e.key==="Escape")window.__CLOSEM()});'
  2278. '})();</script>')
  2279. page = page.replace('<script>%JS%', inject + '<script>%JS%') if '%JS%' in page else page.replace('</head>', inject + '</head>')
  2280. if inject not in page:
  2281. page = page.replace('</head>', inject + '</head>')
  2282. # ★这里**曾经**调 i18n_en.to_english_source(page) 做逐词替换。
  2283. # 它跑在 inject 之后, 作用在含 __SNAP__ JSON 的整页上, 把预烤 key 里的
  2284. # 中文吃掉: win=2026年 → win=2026, win=全程 → win=full window,
  2285. # 于是页面按原 key 取数全部 miss —— **而 miss 是静默的**, 退到 err 提示,
  2286. # 看着像功能没做。_PROT 保护不了 JSON 里的 key。
  2287. # 页面英文化已在前面由 _en_page (整段短语表 + NEVER 数据值保护) 完成,
  2288. # 这里不再需要第二遍, 已删。
  2289. # ★成品页语法闸 —— 校验**最终字节**, 不是中间态 (2026-08-31)。
  2290. # _en_page 里也有一道, 但它跑在 value_script 与 __SNAP__ 注入**之前**,
  2291. # 护不住后注入的脚本。而快照是要上云的成品, 坏了就是白屏,
  2292. # 且白屏时"可见中文=0"会让所有自查指标全绿 —— 最坏的一种失败。
  2293. try:
  2294. from src.windscada import ui_en as _uie
  2295. _nchk = 0
  2296. for _m in re.finditer(r'<script([^>]*)>(.*?)</script>', page, re.S):
  2297. if 'application/json' in _m.group(1) or not _m.group(2).strip():
  2298. continue
  2299. _uie.check_js(_m.group(2), label='快照脚本#%d' % (_nchk + 1))
  2300. _nchk += 1
  2301. print('[snapshot] 成品页语法闸: %d 个脚本全部通过 node --check' % _nchk,
  2302. file=sys.stderr)
  2303. if _nchk == 0:
  2304. print('[snapshot] 🔴 语法闸切出 0 个脚本 —— **本次未校验**, 非通过',
  2305. file=sys.stderr)
  2306. except Exception as _e:
  2307. print('[snapshot] 🔴 成品页语法闸失败: %s' % _e, file=sys.stderr)
  2308. raise
  2309. if _lang == 'en':
  2310. _n1 = len(re.findall(r'[\u4e00-\u9fff]', page))
  2311. print('[snapshot] 英文快照: 页面残留中文 %d 字 (含内联数据)' % _n1,
  2312. file=sys.stderr)
  2313. # ★脱敏必须是**最后一步**, 且只对上云那一版做 (2026-09-01 用户令)。
  2314. # 中文快照是内部工作副本, 要保留真实场名/机型, 不脱敏。
  2315. # 英文快照是要发到公网的, 而它此前带着 rudong×366 / Siemens×194 /
  2316. # SWT×89 —— **同时页面上写着 "site, owner and OEM removed"**,
  2317. # 产物与自己的声明不符, 这比不声明更糟。
  2318. # 放最后: 对象 id 与引用它的地方必须换成同一个词, 只能整体做。
  2319. if _lang == 'en':
  2320. from src.windscada import deid_public as _dp
  2321. page = _dp.scrub_or_die(page, label='英文快照 (上云件)')
  2322. print('[snapshot] 脱敏回扫: 干净 (rudong/OEM/机型/省市/业主/坐标 全零)',
  2323. file=sys.stderr)
  2324. fn = (f'Guanlan_Snapshot_{CUR_WIN[0]}_EN.html' if _lang == 'en'
  2325. else f'观澜快照_{CUR_WIN[0]}_{"上报版" if not INTERNAL else "内部版"}.html')
  2326. data = page.encode('utf-8')
  2327. self.send_response(200)
  2328. self.send_header('Content-Type', 'text/html; charset=utf-8')
  2329. self.send_header('Content-Length', str(len(data)))
  2330. self.send_header('Content-Disposition',
  2331. "attachment; filename*=UTF-8''" + urllib.parse.quote(fn))
  2332. self.end_headers()
  2333. self.wfile.write(data)
  2334. except Exception as e:
  2335. import traceback
  2336. self._send(_jdump(dict(err=str(e)[:200], tb=traceback.format_exc()[-500:]),
  2337. ensure_ascii=False), 'application/json; charset=utf-8', 500)
  2338. elif u.path == '/api/rpt_export':
  2339. # 导出 Word/Excel (2026-08-28 用户令: 导出 demo)。图表块转数据表, 不放位图 —
  2340. # Word 里的位图既不能编辑也无法核对, 上报集团要的是能核的数字。
  2341. q = urllib.parse.parse_qs(u.query)
  2342. fmt = q.get('fmt', ['docx'])[0]
  2343. blocks = [x for x in q.get('blocks', [''])[0].split(',') if x]
  2344. aud = q.get('aud', ['in'])[0]
  2345. w = q.get('win', ['2026年'])[0]
  2346. narr = q.get('narr', [''])[0] or None
  2347. try:
  2348. from src.windscada import report_export as _rx
  2349. fl = fleet_view(w)
  2350. if aud == 'up': # 上报版走同一套脱敏规则
  2351. fl = _redact(_jclean(fl))
  2352. data = (_rx.to_xlsx if fmt == 'xlsx' else _rx.to_docx)(fl, blocks, w, aud, narr)
  2353. ext = 'xlsx' if fmt == 'xlsx' else 'docx'
  2354. mime = ('application/vnd.openxmlformats-officedocument.spreadsheetml.sheet' if fmt == 'xlsx'
  2355. else 'application/vnd.openxmlformats-officedocument.wordprocessingml.document')
  2356. fn = f'运行分析报告_{w}_{"上报版" if aud == "up" else "内部版"}.{ext}'
  2357. self.send_response(200)
  2358. self.send_header('Content-Type', mime)
  2359. self.send_header('Content-Length', str(len(data)))
  2360. self.send_header('Content-Disposition',
  2361. "attachment; filename*=UTF-8''" + urllib.parse.quote(fn))
  2362. self.end_headers()
  2363. self.wfile.write(data)
  2364. except Exception as e:
  2365. import traceback
  2366. self._send(_jdump(dict(err=str(e)[:200], tb=traceback.format_exc()[-400:]),
  2367. ensure_ascii=False), 'application/json; charset=utf-8', 500)
  2368. elif u.path in ('/api/facts', '/api/facts/claim'):
  2369. # 结构化事实契约 (用户批 2026-09-06 "8033 /api 切到契约生成"): 只读 derived/ 产物, 不算数不改写;
  2370. # 生成器 = scripts/guanlan_facts_contract.py render, 契约 sha 随体返回供前端/CI 对拍.
  2371. import json as _json
  2372. _dd = _P.guanlan() / 'derived'
  2373. try:
  2374. if u.path == '/api/facts':
  2375. _body = _dd.joinpath('detail_cards.json').read_text(encoding='utf-8')
  2376. else:
  2377. _cid = urllib.parse.parse_qs(u.query).get('id', [''])[0]
  2378. _qa = _json.loads(_dd.joinpath('qa_refs.json').read_text(encoding='utf-8'))
  2379. _hit = next((r for r in _qa.get('refs', []) if r.get('claim_id') == _cid), None)
  2380. _body = _jdump(dict(contract_sha256=_qa.get('contract_sha256'), claim=_hit, err=None if _hit else 'unknown claim_id'), ensure_ascii=False)
  2381. self._send(_body, 'application/json; charset=utf-8')
  2382. except FileNotFoundError as e:
  2383. self._send(_jdump(dict(err='facts derived 缺失: %s' % e), ensure_ascii=False), 'application/json; charset=utf-8', 503)
  2384. elif u.path == '/api/rpt_compose':
  2385. q = urllib.parse.parse_qs(u.query)
  2386. if q.get('probe', [''])[0]:
  2387. # 前端只探"编排接口是否存在" (v2 bindReport probe=1); 原先照跑 rpt_compose('', [], model) 走模型,
  2388. # 汇报定制页首屏空等 ~23 s (2026-09-06 实测) —— 探针只答存在性, 不算数.
  2389. self._send(_jdump(dict(ok=True, probe=True), ensure_ascii=False), 'application/json; charset=utf-8'); return
  2390. r = rpt_compose(q.get('need', [''])[0], [x for x in q.get('canvas', [''])[0].split(',') if x],
  2391. q.get('model', ['qwen3:8b'])[0])
  2392. self._send(_jdump(r, ensure_ascii=False), 'application/json; charset=utf-8')
  2393. elif u.path == '/api/ask_status':
  2394. self._send(_jdump(ask_status(), ensure_ascii=False), 'application/json; charset=utf-8')
  2395. elif u.path == '/api/ask_models':
  2396. # 与云端问答微服务 (windscada_ask_serve.py) 同形: {models: {id: {...}}}。此前本地没有
  2397. # 这条路由, 经典页的模型下拉在本地永远是"无可用模型" —— v2 与经典页同源取清单。
  2398. self._send(_jdump(dict(models={k: {} for k in ASK_MODELS}), ensure_ascii=False),
  2399. 'application/json; charset=utf-8')
  2400. elif u.path == '/api/ontology':
  2401. self._send(_jdump(ontology_view(), ensure_ascii=False), 'application/json; charset=utf-8')
  2402. elif u.path == '/api/ont_list':
  2403. self._send(_jdump(dict(rows=ontology_list(q.get('type', ''), q.get('t', ''), q.get('limit', 60))),
  2404. ensure_ascii=False), 'application/json; charset=utf-8')
  2405. elif u.path == '/api/ont_obj':
  2406. self._send(_jdump(ontology_obj(q.get('id', '')), ensure_ascii=False), 'application/json; charset=utf-8')
  2407. elif u.path == '/api/channels':
  2408. self._send(_jdump(channel_list(), ensure_ascii=False), 'application/json; charset=utf-8')
  2409. elif u.path.startswith('/static/'):
  2410. f = ST / u.path[len('/static/'):]
  2411. if f.is_file():
  2412. self._send(f.read_bytes(), 'text/html; charset=utf-8')
  2413. else:
  2414. # 裸 404 会让总览 iframe 变成"一块空白"(用户 2026-09-17 报的正是这个现象) ⇒ 结构化缺件页
  2415. self._send(_gap_page(u.path, f, '总览(全场状态)'), 'text/html; charset=utf-8')
  2416. elif u.path.startswith('/turbines/'):
  2417. f = ST / u.path.lstrip('/')
  2418. if f.is_file():
  2419. raw = f.read_text(encoding='utf-8', errors='ignore')
  2420. light = '<style>:root[data-theme="light"]{--ground:#f3f6f9;--panel:#fff;--panel-2:#f8fbfc;--ink:#0c1622;--muted:#5a6b7d;--line:#d8e1e9;--line-2:#eaf0f5;--accent:#0d7183;--accent-soft:#e6f3f5;--accent-ink:#095764;--ok:#137a68;--warn:#8a5a12;--bad:#a32d38;--unknown:#4d5d6e;--info:#0d7183}</style><script>(function(){try{document.documentElement.setAttribute("data-theme","light")}catch(e){}})();</script>'
  2421. raw = raw.replace('</head>', light + '</head>', 1) if '</head>' in raw else light + raw
  2422. self._send(raw)
  2423. else:
  2424. self._send(_gap_page(u.path, f, '机组详情页'), 'text/html; charset=utf-8')
  2425. elif u.path.startswith('/windcms/'):
  2426. # windcms 振动分析本体嵌入 (用户令: 原振动分析板块进入振动融合板块; 静态挂载非复制,
  2427. # 正本仍在 outputs/rudong/windcms 由振动线维护, 本路由只读)
  2428. base = (ST.parent / 'windcms').resolve()
  2429. f = (base / u.path[len('/windcms/'):]).resolve()
  2430. if base in f.parents or f == base:
  2431. if f.is_dir(): f = f / 'index.html'
  2432. if f.is_file():
  2433. ext = f.suffix.lower()
  2434. mime = {'.html': 'text/html; charset=utf-8', '.json': 'application/json; charset=utf-8',
  2435. '.png': 'image/png', '.svg': 'image/svg+xml', '.css': 'text/css',
  2436. '.js': 'application/javascript', '.md': 'text/plain; charset=utf-8'}.get(ext, 'application/octet-stream')
  2437. self._send(f.read_bytes(), mime)
  2438. else:
  2439. self._send('404', code=404)
  2440. else:
  2441. self._send('403', code=403)
  2442. elif u.path == '/api/reload':
  2443. # 产物重载 (2026-09-16 用户令: 重算产物即时可用于页面呈现)。
  2444. # 平时**不需要**调它 —— _load() 每次请求都会比产物指纹, 变了自动重载;
  2445. # 这个端点给"重算结束想立刻生效"的显式动作, 也便于运维界面/脚本确认当前指纹。
  2446. _old = reload_products(urllib.parse.parse_qs(u.query).get('why', ['manual'])[0])
  2447. self._send(_jdump(dict(ok=True, old_stamp=_old, stamp=products_stamp(force=True),
  2448. files=len(_product_files())), ensure_ascii=False),
  2449. 'application/json; charset=utf-8')
  2450. else:
  2451. self._send('404', code=404)
  2452. except ProductsMissing as e:
  2453. # 产物缺失不是"程序坏了": 回 200 + 结构化无数据, 前端按正常渲染显示原因 (2026-09-11 空状态检查)
  2454. self._send(_jdump(dict(kind='multiline', title='无产物', unit='', months=[], series=[],
  2455. err='no_products',
  2456. note=f'{e.what} 不存在 → {_P.rel(e.path)}; 先放原始件再跑 '
  2457. f'scripts/rebuild_from_raw.py (SCADA 侧加 --scada) 或 '
  2458. f'scripts/rebuild_all.py 生成产物。'
  2459. f'"包内没有生成端"的件按用户令 2026-09-17 **运行期不由交付包补齐**: '
  2460. f'要么等研发补生成端(逐族可逆性见 docs/系统设计说明.md §13.6), '
  2461. f'要么由人离线补救(products_restore_missing.py --stash <…>, 不在运行链里); '
  2462. f'2026-09-16 起清除产物不留备份, 故没有 --on 还原'), ensure_ascii=False),
  2463. 'application/json; charset=utf-8')
  2464. except Exception as e:
  2465. import traceback
  2466. self._send(_jdump(dict(err=str(e), tb=traceback.format_exc()[-500:]), ensure_ascii=False), 'application/json; charset=utf-8', 500)
  2467. if __name__ == '__main__':
  2468. # 统一日志口径 (用户令 2): 之后每一行都带 时间戳/级别/组件(详情见 src/logfile.py)
  2469. import pathlib as _pl, sys as _sys0
  2470. _sys0.path.insert(0, str(_install_root(__file__)))
  2471. from src import logfile as _lf; _lf.prefix_stdout('detail')
  2472. # 控制台可能是 GBK(中文 Windows 代码页 936): 正文里的 ✔ ✗ ✅ ⚠ 这类字符编不出来会抛
  2473. # UnicodeEncodeError, 脚本干成了事却以退出码 1 结束(同类坑见 src/console.py)。降级为 '?' 而不是崩;
  2474. # 不用 import 是为了兼顾 python -m 与直接当脚本跑两种启动方式。
  2475. import sys as _sys
  2476. for _s in (_sys.stdout, _sys.stderr):
  2477. try: _s.reconfigure(errors='replace')
  2478. except Exception: pass
  2479. import argparse
  2480. ap = argparse.ArgumentParser()
  2481. import os
  2482. ap.add_argument('--port', type=int, default=int(os.environ.get('PORT', 8033))) # PORT 环境变量优先, 供 preview 分配端口
  2483. ap.add_argument('--host', default='127.0.0.1', help='0.0.0.0 = 对局域网开放 (演示用, 结束请关)')
  2484. ap.add_argument('--farm', default=None, help='场名 (默认 rudong; 可用场见 app_ETL/configs/farms/)')
  2485. a = ap.parse_args()
  2486. print(f'windscada serve :{a.port}', flush=True)
  2487. # --host 0.0.0.0 供局域网演示 (2026-08-28 用户: 怎么发给别人)。
  2488. # 默认仍只绑本机 — 内部数据不该因为一次演示就对整个网段开放, 要开是显式动作。
  2489. if getattr(a, 'farm', None):
  2490. # ★CFG/ST 在模块加载时就绑定了 (第 81 行), 这里必须显式重绑而不是只切 current —
  2491. # 否则 --farm 只改了配置模块里的指针, 服务用的仍是启动时那个场 (静默不生效)。
  2492. # 用模块字典赋值而非 global: main 里更早处已引用过 CFG, global 声明会触发
  2493. # "used prior to global declaration"。改模块属性等效且不受声明位置约束。
  2494. from src.windscada import config as _cfg
  2495. import sys as _sysmod
  2496. _me = _sysmod.modules[__name__]
  2497. _cfg.set_current(a.farm)
  2498. _me.CFG = _cfg.farm()
  2499. _me.ST = pathlib.Path(_me.CFG['store'])
  2500. _CACHE.clear()
  2501. _STAMP['val'] = None # 换场后产物指纹必须重算, 否则新场的缓存判为"未变化"
  2502. print(f'场: {_me.CFG["name"]} ({a.farm}, {_me.CFG["n_turbines"]} 台, 仓 {_me.ST})')
  2503. else:
  2504. print(f'场: {CFG["name"]} (默认; 可用 {list(__import__("src.windscada.config", fromlist=["available"]).available())})')
  2505. _host = getattr(a, 'host', '127.0.0.1')
  2506. if _host != '127.0.0.1':
  2507. import socket as _sk
  2508. try:
  2509. _ip = _sk.gethostbyname(_sk.gethostname())
  2510. except Exception:
  2511. _ip = _host
  2512. print(f'⚠ 已对局域网开放: http://{_ip}:{a.port} (同网段任何人可访问, 演示结束请关闭)')
  2513. # ★按窗重算的**预热**(用户令 2026-09-21:判级/曲线随所选时间窗变化)。
  2514. # 判级矩阵 ≈25 s/窗、镜头 ≈10~45 s/窗 —— 不预热的话重启后第一个访客要等 pending 轮询。
  2515. # 这里在**后台守护线程**里顺序把预设窗算进进程缓存(页面照常可用,未就位的窗如实标 pending)。
  2516. def _warm():
  2517. if _os.environ.get('GUANLAN_WIN_WARM', '1') == '0':
  2518. print('[warm] 按窗预热被 GUANLAN_WIN_WARM=0 关掉(按需算,首次看窗会 pending 一轮)', flush=True)
  2519. return
  2520. import time as _t
  2521. n = 0
  2522. for w in list(WINDOWS):
  2523. try:
  2524. sysmx_of(w) # 未命中即起后台线程
  2525. if _win_wait('sysmx', w, timeout=180) is not None:
  2526. n += 1
  2527. _t.sleep(2) # 留口气: 小内存机器上连续重算更容易被系统盯上
  2528. except Exception as e:
  2529. print(f'[warm] 判级 {w} 失败: {type(e).__name__}: {e}'[:160], flush=True)
  2530. for w in ('2026年',):
  2531. try:
  2532. curves_of(w) # 未命中即起后台线程
  2533. _win_wait('curves', w, timeout=240)
  2534. except Exception as e:
  2535. print(f'[warm] 镜头 {w} 失败: {type(e).__name__}: {e}'[:160], flush=True)
  2536. print(f'[warm] 按窗缓存预热完成: 判级 {n}/{len(WINDOWS)} 窗(可用内存 {_mem_mb()} MB)', flush=True)
  2537. _thr.Thread(target=_warm, name='win-warm', daemon=True).start()
  2538. ThreadingHTTPServer((_host, a.port), H).serve_forever()