railwise-ai 1.2.32 → 1.2.34

package/agent-pack/assets/tool/cross_section.ts

@@ -0,0 +1,356 @@
+/// <reference path="../env.d.ts" />
+import { tool } from "nb-railwise/tool"
+
+// ============================================================
+// Tunnel / excavation cross-section analysis tools
+// ============================================================
+
+// ============================================================
+// Tool: Convergence calculation from monitoring pairs
+// ============================================================
+
+export const convergence_calc = tool({
+  description:
+    "收敛量计算：根据隧道断面上对称测点对的坐标，计算净空收敛量（水平收敛、拱顶沉降等）。地铁隧道健康监测和施工监测中的核心分析工具。data_analyst 在处理隧道断面监测数据时必须调用此工具。",
+  args: {
+    sectionId: tool.schema.string().describe("断面编号"),
+    chainage: tool.schema.string().optional().describe("里程桩号，如 K12+345.678"),
+    initialPairs: tool.schema
+      .array(
+        tool.schema.object({
+          pairId: tool.schema.string().describe("测点对编号，如 H1（水平收敛对1）"),
+          pointA: tool.schema.object({
+            x: tool.schema.number().describe("A点X坐标(m)"),
+            y: tool.schema.number().describe("A点Y坐标(m)"),
+          }),
+          pointB: tool.schema.object({
+            x: tool.schema.number().describe("B点X坐标(m)"),
+            y: tool.schema.number().describe("B点Y坐标(m)"),
+          }),
+        }),
+      )
+      .min(1)
+      .describe("初始观测（基准）的测点对坐标"),
+    currentPairs: tool.schema
+      .array(
+        tool.schema.object({
+          pairId: tool.schema.string().describe("测点对编号（与初始观测对应）"),
+          pointA: tool.schema.object({
+            x: tool.schema.number().describe("A点X坐标(m)"),
+            y: tool.schema.number().describe("A点Y坐标(m)"),
+          }),
+          pointB: tool.schema.object({
+            x: tool.schema.number().describe("B点X坐标(m)"),
+            y: tool.schema.number().describe("B点Y坐标(m)"),
+          }),
+        }),
+      )
+      .min(1)
+      .describe("本期观测的测点对坐标"),
+    crownInitial: tool.schema
+      .object({
+        x: tool.schema.number().describe("拱顶X坐标(m)"),
+        y: tool.schema.number().describe("拱顶Y坐标(m)"),
+      })
+      .optional()
+      .describe("拱顶初始坐标（若提供则计算拱顶沉降）"),
+    crownCurrent: tool.schema
+      .object({
+        x: tool.schema.number().describe("拱顶X坐标(m)"),
+        y: tool.schema.number().describe("拱顶Y坐标(m)"),
+      })
+      .optional()
+      .describe("拱顶本期坐标"),
+    alertThreshold: tool.schema.number().positive().optional().describe("收敛报警值(mm)"),
+  },
+  async execute(args) {
+    // Build initial pair map
+    const initialMap = new Map(args.initialPairs.map((p) => [p.pairId, p]))
+
+    const convergences = args.currentPairs.map((curr) => {
+      const init = initialMap.get(curr.pairId)
+      if (!init)
+        return {
+          pair_id: curr.pairId,
+          error: `未找到对应的初始观测数据`,
+        }
+
+      const initDist = Math.sqrt((init.pointB.x - init.pointA.x) ** 2 + (init.pointB.y - init.pointA.y) ** 2)
+      const currDist = Math.sqrt((curr.pointB.x - curr.pointA.x) ** 2 + (curr.pointB.y - curr.pointA.y) ** 2)
+      const convergence = (initDist - currDist) * 1000 // mm, positive = converging
+      const rate = initDist > 0 ? (convergence / (initDist * 1000)) * 100 : 0
+
+      let status = "🟢 正常"
+      if (args.alertThreshold) {
+        const ratio = Math.abs(convergence) / args.alertThreshold
+        if (ratio >= 1.0) status = "🔴 超限"
+        else if (ratio >= 0.85) status = "🟠 接近阈值"
+        else if (ratio >= 0.7) status = "🟡 关注"
+      }
+
+      return {
+        pair_id: curr.pairId,
+        initial_distance_m: Number(initDist.toFixed(4)),
+        current_distance_m: Number(currDist.toFixed(4)),
+        convergence_mm: Number(convergence.toFixed(3)),
+        convergence_rate_pct: Number(rate.toFixed(4)),
+        status,
+      }
+    })
+
+    // Crown settlement
+    let crownSettlement: Record<string, unknown> | undefined
+    if (args.crownInitial && args.crownCurrent) {
+      const dy = (args.crownCurrent.y - args.crownInitial.y) * 1000 // mm
+      const dx = (args.crownCurrent.x - args.crownInitial.x) * 1000 // mm
+      const total = Math.sqrt(dx * dx + dy * dy)
+
+      let status = "🟢 正常"
+      if (args.alertThreshold) {
+        const ratio = Math.abs(dy) / args.alertThreshold
+        if (ratio >= 1.0) status = "🔴 超限"
+        else if (ratio >= 0.85) status = "🟠 接近阈值"
+        else if (ratio >= 0.7) status = "🟡 关注"
+      }
+
+      crownSettlement = {
+        vertical_mm: Number(dy.toFixed(3)),
+        horizontal_mm: Number(dx.toFixed(3)),
+        total_mm: Number(total.toFixed(3)),
+        direction: dy < 0 ? "沉降" : "上抬",
+        status,
+      }
+    }
+
+    const maxConvergence = convergences
+      .filter((c) => "convergence_mm" in c)
+      .reduce(
+        (max, c) => {
+          const val = Math.abs((c as { convergence_mm: number }).convergence_mm)
+          return val > max.value ? { pair_id: c.pair_id, value: val } : max
+        },
+        { pair_id: "", value: 0 },
+      )
+
+    return JSON.stringify({
+      section_id: args.sectionId,
+      chainage: args.chainage ?? "未指定",
+      convergences,
+      crown_settlement: crownSettlement,
+      max_convergence: {
+        pair_id: maxConvergence.pair_id,
+        value_mm: maxConvergence.value,
+      },
+      message: `✅ 断面 ${args.sectionId} 收敛分析完成：最大收敛 ${maxConvergence.pair_id}(${maxConvergence.value.toFixed(3)}mm)${crownSettlement ? `，拱顶${(crownSettlement as { direction: string }).direction} ${Math.abs(crownSettlement.vertical_mm as number).toFixed(3)}mm` : ""}`,
+    })
+  },
+})
+
+// ============================================================
+// Tool: Cross-section profile comparison
+// ============================================================
+
+export const profile_comparison = tool({
+  description:
+    "断面轮廓对比分析：将隧道/基坑实测断面轮廓与设计轮廓进行对比，计算各点的超欠挖量和净空余量。适用于隧道开挖断面检测、盾构管片姿态分析。",
+  args: {
+    sectionId: tool.schema.string().describe("断面编号"),
+    designProfile: tool.schema
+      .array(
+        tool.schema.object({
+          angle: tool.schema.number().describe("从拱顶顺时针方向的角度(度)，0=拱顶，180=仰拱"),
+          radius: tool.schema.number().positive().describe("设计半径/距中心距离(m)"),
+        }),
+      )
+      .min(3)
+      .describe("设计断面轮廓（极坐标表示：角度+半径）"),
+    measuredProfile: tool.schema
+      .array(
+        tool.schema.object({
+          angle: tool.schema.number().describe("测量方向角(度)"),
+          radius: tool.schema.number().positive().describe("实测距中心距离(m)"),
+        }),
+      )
+      .min(3)
+      .describe("实测断面轮廓"),
+    overbreakLimit: tool.schema.number().positive().default(150).describe("允许超挖限值(mm)，默认150mm"),
+    underbreakLimit: tool.schema.number().positive().default(0).describe("允许欠挖限值(mm)，默认0mm（不允许欠挖）"),
+  },
+  async execute(args) {
+    // Sort design profile by angle for interpolation
+    const design = [...args.designProfile].sort((a, b) => a.angle - b.angle)
+
+    // Interpolate design radius at given angle
+    function designRadius(angle: number): number {
+      const normAngle = ((angle % 360) + 360) % 360
+      for (let i = 0; i < design.length - 1; i++) {
+        if (normAngle >= design[i]!.angle && normAngle <= design[i + 1]!.angle) {
+          const ratio = (normAngle - design[i]!.angle) / (design[i + 1]!.angle - design[i]!.angle)
+          return design[i]!.radius + ratio * (design[i + 1]!.radius - design[i]!.radius)
+        }
+      }
+      // Wrap around or nearest
+      return design[0]!.radius
+    }
+
+    const comparisons = args.measuredProfile.map((m) => {
+      const dRadius = designRadius(m.angle)
+      const diff = (m.radius - dRadius) * 1000 // mm, positive = overbreak
+      const isOverbreak = diff > 0
+      const isUnderbreak = diff < 0
+
+      let status = "🟢 合格"
+      if (isOverbreak && diff > args.overbreakLimit) {
+        status = "🔴 超挖超限"
+      } else if (isUnderbreak && Math.abs(diff) > args.underbreakLimit) {
+        status = "🔴 欠挖"
+      } else if (isOverbreak) {
+        status = "🟡 超挖（限内）"
+      }
+
+      return {
+        angle_deg: m.angle,
+        design_radius_m: Number(dRadius.toFixed(4)),
+        measured_radius_m: Number(m.radius.toFixed(4)),
+        deviation_mm: Number(diff.toFixed(1)),
+        type: isOverbreak ? "超挖" : isUnderbreak ? "欠挖" : "吻合",
+        status,
+      }
+    })
+
+    const overbreakPoints = comparisons.filter((c) => c.type === "超挖")
+    const underbreakPoints = comparisons.filter((c) => c.type === "欠挖")
+    const maxOverbreak =
+      overbreakPoints.length > 0
+        ? overbreakPoints.reduce((max, c) => (c.deviation_mm > max.deviation_mm ? c : max), overbreakPoints[0]!)
+        : null
+    const maxUnderbreak =
+      underbreakPoints.length > 0
+        ? underbreakPoints.reduce(
+            (max, c) => (Math.abs(c.deviation_mm) > Math.abs(max.deviation_mm) ? c : max),
+            underbreakPoints[0]!,
+          )
+        : null
+
+    const avgDeviation = comparisons.reduce((s, c) => s + c.deviation_mm, 0) / comparisons.length
+    const failCount = comparisons.filter((c) => c.status.includes("超限") || c.status.includes("欠挖")).length
+
+    return JSON.stringify({
+      section_id: args.sectionId,
+      measured_points: comparisons.length,
+      overbreak_count: overbreakPoints.length,
+      underbreak_count: underbreakPoints.length,
+      fail_count: failCount,
+      average_deviation_mm: Number(avgDeviation.toFixed(1)),
+      max_overbreak: maxOverbreak
+        ? {
+            angle_deg: maxOverbreak.angle_deg,
+            deviation_mm: maxOverbreak.deviation_mm,
+          }
+        : null,
+      max_underbreak: maxUnderbreak
+        ? {
+            angle_deg: maxUnderbreak.angle_deg,
+            deviation_mm: maxUnderbreak.deviation_mm,
+          }
+        : null,
+      details: comparisons,
+      assessment:
+        failCount === 0
+          ? `✅ 断面 ${args.sectionId} 开挖轮廓合格：${overbreakPoints.length}处超挖均在限值内，无欠挖`
+          : `❌ 断面 ${args.sectionId} 有 ${failCount} 处不合格（超挖超限或欠挖），需处理`,
+      message: `断面 ${args.sectionId}：${comparisons.length}个测点，超挖 ${overbreakPoints.length}处，欠挖 ${underbreakPoints.length}处，${failCount === 0 ? "全部合格" : `${failCount}处不合格`}`,
+    })
+  },
+})
+
+// ============================================================
+// Tool: Section clearance check
+// ============================================================
+
+export const clearance_check = tool({
+  description:
+    "限界（建筑净空）检查：对比隧道实测断面与限界要求，判断是否侵入建筑限界。运营隧道健康监测中用于结构变形是否影响行车安全的关键判定。",
+  args: {
+    sectionId: tool.schema.string().describe("断面编号"),
+    gaugeProfile: tool.schema
+      .array(
+        tool.schema.object({
+          angle: tool.schema.number().describe("方向角(度)"),
+          minRadius: tool.schema.number().positive().describe("该方向上的最小净空要求(m)"),
+        }),
+      )
+      .min(3)
+      .describe("限界轮廓（极坐标，定义各方向最小净空）"),
+    measuredProfile: tool.schema
+      .array(
+        tool.schema.object({
+          angle: tool.schema.number().describe("测量方向角(度)"),
+          radius: tool.schema.number().positive().describe("实测距中心距离(m)"),
+        }),
+      )
+      .min(3)
+      .describe("实测断面轮廓"),
+    safetyMargin: tool.schema.number().default(50).describe("安全余量(mm)，低于此值发出预警，默认50mm"),
+  },
+  async execute(args) {
+    const gauge = [...args.gaugeProfile].sort((a, b) => a.angle - b.angle)
+
+    function gaugeRadius(angle: number): number {
+      const norm = ((angle % 360) + 360) % 360
+      for (let i = 0; i < gauge.length - 1; i++) {
+        if (norm >= gauge[i]!.angle && norm <= gauge[i + 1]!.angle) {
+          const ratio = (norm - gauge[i]!.angle) / (gauge[i + 1]!.angle - gauge[i]!.angle)
+          return gauge[i]!.minRadius + ratio * (gauge[i + 1]!.minRadius - gauge[i]!.minRadius)
+        }
+      }
+      return gauge[0]!.minRadius
+    }
+
+    const checks = args.measuredProfile.map((m) => {
+      const minR = gaugeRadius(m.angle)
+      const margin = (m.radius - minR) * 1000 // mm
+
+      let status: string
+      if (margin < 0) {
+        status = "🔴 侵入限界"
+      } else if (margin < args.safetyMargin) {
+        status = "🟠 余量不足"
+      } else {
+        status = "🟢 合格"
+      }
+
+      return {
+        angle_deg: m.angle,
+        gauge_min_radius_m: Number(minR.toFixed(4)),
+        measured_radius_m: Number(m.radius.toFixed(4)),
+        margin_mm: Number(margin.toFixed(1)),
+        status,
+      }
+    })
+
+    const violations = checks.filter((c) => c.margin_mm < 0)
+    const warnings = checks.filter((c) => c.margin_mm >= 0 && c.margin_mm < args.safetyMargin)
+    const minMargin = checks.reduce((min, c) => (c.margin_mm < min.margin_mm ? c : min), checks[0]!)
+
+    return JSON.stringify({
+      section_id: args.sectionId,
+      total_points: checks.length,
+      violations: violations.length,
+      warnings: warnings.length,
+      min_margin: {
+        angle_deg: minMargin.angle_deg,
+        margin_mm: minMargin.margin_mm,
+        status: minMargin.status,
+      },
+      details: checks,
+      assessment:
+        violations.length > 0
+          ? `🔴 断面 ${args.sectionId} 有 ${violations.length} 处侵入建筑限界！最小余量 ${minMargin.margin_mm.toFixed(1)}mm（${minMargin.angle_deg}°方向），必须立即处理！`
+          : warnings.length > 0
+            ? `🟠 断面 ${args.sectionId} 有 ${warnings.length} 处净空余量不足 ${args.safetyMargin}mm，需密切关注`
+            : `✅ 断面 ${args.sectionId} 建筑限界检查通过，最小余量 ${minMargin.margin_mm.toFixed(1)}mm`,
+      message: `限界检查 ${args.sectionId}：${violations.length > 0 ? `🔴 ${violations.length}处侵限` : warnings.length > 0 ? `🟠 ${warnings.length}处余量不足` : "✅ 全部通过"}，最小余量 ${minMargin.margin_mm.toFixed(1)}mm`,
+    })
+  },
+})

package/agent-pack/assets/tool/deformation_rate.ts

@@ -0,0 +1,234 @@
+/// <reference path="../env.d.ts" />
+import { tool } from "nb-railwise/tool"
+
+// ============================================================
+// Deformation rate analysis & trend prediction
+// ============================================================
+
+// ============================================================
+// Tool: Deformation rate calculation
+// ============================================================
+
+export const deformation_rate = tool({
+  description:
+    "变形速率计算与趋势分析。根据监测点的时间-变形量序列，计算各期变形速率、累计变形量，并用线性回归进行趋势预测。城市轨道监测中判断变形是否收敛的核心分析工具。data_analyst 在分析自动化监测数据趋势时必须调用此工具。",
+  args: {
+    pointId: tool.schema.string().describe("监测点编号"),
+    data: tool.schema
+      .array(
+        tool.schema.object({
+          date: tool.schema.string().describe("观测日期，格式 YYYY-MM-DD 或 YYYY-MM-DD HH:mm"),
+          value: tool.schema.number().describe("累计变形量(mm)，正值表示沉降/收敛方向"),
+        }),
+      )
+      .min(2)
+      .describe("按时间顺序排列的监测数据序列"),
+    alertThreshold: tool.schema.number().positive().optional().describe("报警控制值(mm)，若提供则输出预警分析"),
+    rateThreshold: tool.schema.number().positive().optional().describe("速率控制值(mm/d)，若提供则判断速率是否超限"),
+    predictionDays: tool.schema.number().int().positive().default(7).describe("向前预测天数，默认7天"),
+  },
+  async execute(args) {
+    const n = args.data.length
+
+    // Parse dates to day offsets
+    const t0 = new Date(args.data[0]!.date).getTime()
+    const days = args.data.map((d) => (new Date(d.date).getTime() - t0) / 86400000)
+    const values = args.data.map((d) => d.value)
+
+    // Period-by-period rates
+    const rates: Array<{ period: string; rate_mm_per_day: number; increment_mm: number; days: number }> = []
+    for (let i = 1; i < n; i++) {
+      const dt = days[i]! - days[i - 1]!
+      const dv = values[i]! - values[i - 1]!
+      const rate = dt > 0 ? dv / dt : 0
+      rates.push({
+        period: `${args.data[i - 1]!.date} → ${args.data[i]!.date}`,
+        rate_mm_per_day: Number(rate.toFixed(4)),
+        increment_mm: Number(dv.toFixed(4)),
+        days: Number(dt.toFixed(2)),
+      })
+    }
+
+    // Overall statistics
+    const totalDays = days[n - 1]! - days[0]!
+    const totalDeformation = values[n - 1]! - values[0]!
+    const avgRate = totalDays > 0 ? totalDeformation / totalDays : 0
+    const latestValue = values[n - 1]!
+    const latestRate = rates.length > 0 ? rates[rates.length - 1]!.rate_mm_per_day : 0
+
+    // Linear regression: y = a + b*x
+    const meanX = days.reduce((s, v) => s + v, 0) / n
+    const meanY = values.reduce((s, v) => s + v, 0) / n
+    const ssxy = days.reduce((s, x, i) => s + (x - meanX) * (values[i]! - meanY), 0)
+    const ssxx = days.reduce((s, x) => s + (x - meanX) * (x - meanX), 0)
+    const b = ssxx > 0 ? ssxy / ssxx : 0
+    const a = meanY - b * meanX
+
+    // R² (coefficient of determination)
+    const ssyy = values.reduce((s, y) => s + (y - meanY) * (y - meanY), 0)
+    const r2 = ssyy > 0 ? (ssxy * ssxy) / (ssxx * ssyy) : 0
+
+    // Residual standard error
+    const residuals = values.map((y, i) => y - (a + b * days[i]!))
+    const rmse = Math.sqrt(residuals.reduce((s, r) => s + r * r, 0) / Math.max(n - 2, 1))
+
+    // Prediction
+    const lastDay = days[n - 1]!
+    const predictions = Array.from({ length: args.predictionDays }, (_, i) => {
+      const predDay = lastDay + i + 1
+      const predDate = new Date(t0 + predDay * 86400000)
+      const predValue = a + b * predDay
+      return {
+        date: predDate.toISOString().slice(0, 10),
+        predicted_mm: Number(predValue.toFixed(4)),
+        day_offset: Number(predDay.toFixed(1)),
+      }
+    })
+
+    // Stability assessment
+    const last3Rates = rates.slice(-3).map((r) => Math.abs(r.rate_mm_per_day))
+    const avgLast3Rate = last3Rates.length > 0 ? last3Rates.reduce((s, v) => s + v, 0) / last3Rates.length : 0
+
+    let stability: string
+    if (avgLast3Rate < 0.01) {
+      stability = "✅ 已收敛：近期速率趋近于零，变形基本稳定"
+    } else if (avgLast3Rate < 0.05) {
+      stability = "🟢 趋于收敛：变形速率逐渐减小"
+    } else if (Math.abs(latestRate) > Math.abs(avgRate) * 1.5) {
+      stability = "🔴 加速变形：最新速率明显大于平均速率，需加密监测"
+    } else if (avgLast3Rate < Math.abs(avgRate)) {
+      stability = "🟡 减速变形：速率有所减小但尚未收敛，继续监测"
+    } else {
+      stability = "🟠 等速变形：速率基本稳定，关注发展趋势"
+    }
+
+    // Alert analysis
+    let alertAnalysis: Record<string, unknown> | undefined
+    if (args.alertThreshold) {
+      const ratio = Math.abs(latestValue) / args.alertThreshold
+      const predMax = Math.max(...predictions.map((p) => Math.abs(p.predicted_mm)))
+      const predRatio = predMax / args.alertThreshold
+      const daysToThreshold = b !== 0 ? (args.alertThreshold * Math.sign(b) - a) / b - lastDay : Infinity
+
+      alertAnalysis = {
+        current_ratio_pct: Number((ratio * 100).toFixed(1)),
+        predicted_max_ratio_pct: Number((predRatio * 100).toFixed(1)),
+        estimated_days_to_threshold:
+          daysToThreshold > 0 && isFinite(daysToThreshold)
+            ? Number(daysToThreshold.toFixed(1))
+            : "不会达到（趋势方向相反或速率为零）",
+        alert_level:
+          ratio >= 1.0 ? "🔴 已超阈值" : ratio >= 0.85 ? "🟠 接近阈值" : ratio >= 0.7 ? "🟡 需关注" : "🟢 正常",
+      }
+    }
+
+    // Rate threshold check
+    let rateAlert: string | undefined
+    if (args.rateThreshold) {
+      rateAlert =
+        Math.abs(latestRate) > args.rateThreshold
+          ? `🔴 最新速率 ${Math.abs(latestRate).toFixed(4)} mm/d 超过限值 ${args.rateThreshold} mm/d`
+          : `🟢 最新速率 ${Math.abs(latestRate).toFixed(4)} mm/d 在限值 ${args.rateThreshold} mm/d 内`
+    }
+
+    return JSON.stringify({
+      point_id: args.pointId,
+      data_count: n,
+      monitoring_duration_days: Number(totalDays.toFixed(1)),
+      latest_value_mm: latestValue,
+      total_deformation_mm: Number(totalDeformation.toFixed(4)),
+      average_rate_mm_per_day: Number(avgRate.toFixed(4)),
+      latest_rate_mm_per_day: latestRate,
+      rates,
+      regression: {
+        equation: `y = ${a.toFixed(4)} + ${b.toFixed(4)} × t`,
+        slope_mm_per_day: Number(b.toFixed(4)),
+        intercept_mm: Number(a.toFixed(4)),
+        r_squared: Number(r2.toFixed(4)),
+        rmse_mm: Number(rmse.toFixed(4)),
+      },
+      predictions,
+      stability_assessment: stability,
+      alert_analysis: alertAnalysis,
+      rate_alert: rateAlert,
+      message: `✅ ${args.pointId} 变形分析：累计 ${latestValue}mm，最新速率 ${latestRate}mm/d，${stability}`,
+    })
+  },
+})
+
+// ============================================================
+// Tool: Multi-point deformation comparison
+// ============================================================
+
+export const deformation_comparison = tool({
+  description:
+    "多测点变形对比分析。同时对比多个监测点的变形量和速率，找出最大变形点、异常点。用于编制监测日报/周报中的断面对比分析。",
+  args: {
+    points: tool.schema
+      .array(
+        tool.schema.object({
+          id: tool.schema.string().describe("测点编号"),
+          latestValue: tool.schema.number().describe("最新累计变形量(mm)"),
+          previousValue: tool.schema.number().describe("上期累计变形量(mm)"),
+          daysBetween: tool.schema.number().positive().describe("两期间隔天数"),
+        }),
+      )
+      .min(1)
+      .describe("各监测点数据"),
+    alertThreshold: tool.schema.number().positive().optional().describe("统一报警控制值(mm)"),
+    rateThreshold: tool.schema.number().positive().optional().describe("速率控制值(mm/d)"),
+  },
+  async execute(args) {
+    const analyzed = args.points.map((p) => {
+      const increment = p.latestValue - p.previousValue
+      const rate = increment / p.daysBetween
+      const absVal = Math.abs(p.latestValue)
+
+      let status = "🟢 正常"
+      if (args.alertThreshold) {
+        const ratio = absVal / args.alertThreshold
+        if (ratio >= 1.0) status = "🔴 超限"
+        else if (ratio >= 0.85) status = "🟠 接近阈值"
+        else if (ratio >= 0.7) status = "🟡 关注"
+      }
+      if (args.rateThreshold && Math.abs(rate) > args.rateThreshold) {
+        status = "🔴 速率超限"
+      }
+
+      return {
+        point_id: p.id,
+        latest_mm: p.latestValue,
+        increment_mm: Number(increment.toFixed(4)),
+        rate_mm_per_day: Number(rate.toFixed(4)),
+        status,
+      }
+    })
+
+    // Sort by absolute latest value descending
+    const sorted = [...analyzed].sort((a, b) => Math.abs(b.latest_mm) - Math.abs(a.latest_mm))
+    const maxPoint = sorted[0]!
+    const alertCount = analyzed.filter((a) => a.status.includes("超限") || a.status.includes("接近")).length
+
+    const avgDeformation = analyzed.reduce((s, a) => s + Math.abs(a.latest_mm), 0) / analyzed.length
+    const maxRate = analyzed.reduce(
+      (max, a) => (Math.abs(a.rate_mm_per_day) > Math.abs(max.rate_mm_per_day) ? a : max),
+      analyzed[0]!,
+    )
+
+    return JSON.stringify({
+      total_points: analyzed.length,
+      alert_count: alertCount,
+      max_deformation: {
+        point_id: maxPoint.point_id,
+        value_mm: maxPoint.latest_mm,
+      },
+      max_rate: {
+        point_id: maxRate.point_id,
+        rate_mm_per_day: maxRate.rate_mm_per_day,
+      },
+      average_deformation_mm: Number(avgDeformation.toFixed(4)),
+      details: sorted,
+      message: `✅ ${analyzed.length}个测点对比：最大变形 ${maxPoint.point_id}(${maxPoint.latest_mm}mm)，最大速率 ${maxRate.point_id}(${maxRate.rate_mm_per_day}mm/d)，${alertCount}个测点预警`,
+    })
+  },
+})