---
title: Q/CR 824-2021 轨道短波指标深度解析：从标准到改进
author: 齐钒羽
date: 2026-03-23
summary: 逐条拆解 Q/CR 824-2021 附录A/B 中 TII 和 RCI 指标的计算流程，详解改进处理流程（粗定位+细分析），分析速度依赖问题与空间域解决方案
canonical: https://blog.sysu-sam.com/@qifanyu/qcr824-indicator-analysis
---

# Q/CR 824-2021 轨道短波指标深度解析：从标准到改进

## 背景

Q/CR 824-2021《铁路基础设施动态检测 轨道短波状态车辆动态响应检测及评价方法》定义了两个核心指标：

- **TII（轨道冲击指数）**：评价单点冲击型轨道短波病害（焊缝、接头、掉块）
- **RCI（钢轨波磨指数）**：评价区段周期类轨道短波病害（钢轨波磨）

本文基于实测数据（8车编组 67 km/h、6车编组 88.5 km/h，同一 80km 区段），逐条拆解标准附录中的计算步骤，找出实现中的速度依赖问题，并介绍一套完整的改进处理流程。

---

## 一、数据结构：等时间采样 + 1m 里程打包

原始数据的关键特征：

| 属性 | 值 |
|------|-----|
| 采样方式 | 等时间采样，5000 Hz |
| 里程分辨率 | 1m（0.001 km），同一米内所有点标记相同里程值 |
| 速度分辨率 | 1m，同一米内速度恒定 |
| 每米采样点数 | ~250点（67 km/h）或 ~200点（88 km/h） |

原始里程是阶梯状的——每跨过 1m 边界才跳变一次：

```
点号:    0      1      2    ... 251    252    253  ...
里程: 225.501 225.501 225.501 ... 225.501 225.502 225.502 ...
      |<-- 252点，同一个里程值 -->|<-- 下一米 -->|
```

**里程重建**：使用逐 1m 步长内的瞬时速度积分，将阶梯里程平滑为连续单调里程，每个 1m 边界精确锚定。

---

## 二、TII 计算流程逐步拆解（附录A）

### A.2 数据预处理：带通滤波

> 标准原文：在 20Hz ~ 500Hz 范围内对轴箱垂向加速度数据进行带通滤波。

TII 的滤波截止频率直接以 **Hz** 给出，**不涉及速度转换**。这是 TII 和 RCI 的关键区别。

### A.3 移动有效值计算

> 标准原文：采用 **3m 的窗长**移动计算带通滤波后的轴箱垂向加速度的有效值。

$$a_{i,\text{RMS}} = \sqrt{\frac{1}{K} \sum_{j=i}^{i+K-1} a_j^2}$$

这里的"3m"是**轨道上的空间距离**，不是时间。但数据是等时间采样的（5000Hz），所以需要把 3m 转换成采样点数：

$$K = \frac{3\text{m}}{v\text{(m/s)}} \times 5000\text{Hz}$$

| 速度 | 3m 对应点数 |
|------|------------|
| 60 km/h | 900 点 |
| 67 km/h | 806 点 |
| 80 km/h | 675 点 |
| 90 km/h | 600 点 |

**问题**：原始实现使用全段平均速度（一个标量）计算窗长，整段 80km 用同一个固定点数窗口。当速度变化时，实际覆盖的空间距离会偏离 3m。

### A.4 单元有效值

> 将轨道按照 50m 长度划分统计单元，将该单元内的移动有效值最大值记为单元有效值。

按里程分段，每段取 max。这一步不涉及速度，无需修改。

### A.5 归一化

> 轴箱垂向加速度移动有效值除以**该速度等级对应的**单元有效值的平均值，得到轨道冲击指数。

$$\lambda = \frac{a_{\text{RMS}}}{\overline{a}_{\text{RMS}}}$$

其中 $\overline{a}_{\text{RMS}}$ 是**全段所有 50m 单元最大值的算术平均**——这个"全局"是标准要求的，不是简化。目的是建立"正常水平"参考基准，使 TII 表示"该点相对于全段平均异常多少倍"。

"该速度等级"指的是铁路线路的设计速度等级分类（60/80/120/160 km/h 等），不是瞬时速度。同一次检测运行属于同一速度等级。

### TII 速度依赖汇总

| 步骤 | 是否涉及速度 | 原实现 |
|------|------------|--------|
| A.2 带通 20-500Hz | **否** | 正确 |
| A.3 移动窗 3m → 点数 | **是** | 用了全局均速 |
| A.4 50m 分段取 max | **否** | 正确 |
| A.5 归一化（全段均值） | **否**（标准要求全局） | 正确 |

---

## 三、RCI 计算流程逐步拆解（附录B）

### B.2 数据预处理：波长域带通滤波

> 标准原文：分别在 40mm~150mm 和 150mm~300mm 范围内对轴箱垂向加速度数据进行带通滤波。

RCI 的滤波范围以**波长（mm）**给出，但 Butterworth 滤波器接受的是频率（Hz）。转换关系：

$$f(\text{Hz}) = \frac{v(\text{m/s})}{\lambda(\text{m})}$$

| 波长 | 60 km/h | 67 km/h | 80 km/h | 90 km/h |
|------|---------|---------|---------|---------|
| 300mm → Hz | 55.6 | 62.0 | 74.1 | 83.3 |
| 40mm → Hz | 416.7 | 465.2 | 555.6 | 625.0 |

**这是 RCI 最核心的速度依赖点**：同一个物理波长（如 100mm 波磨），在不同速度下对应不同的时域频率。如果用全局均速设计滤波器，那么在实际速度偏离均值的位置，滤波器的通带对目标波长的覆盖就会偏移。

### B.3 移动有效值

> 采用 **30m 的窗长**移动计算。

与 TII 的 3m 窗同理，30m 也是空间距离，需要用速度换算成时间点数。RCI 用 30m 是因为波磨是周期性病害，需要覆盖多个完整波长周期。

### RCI 速度依赖汇总

| 步骤 | 是否涉及速度 | 原实现 |
|------|------------|--------|
| B.2 波长 → Hz 转换 | **是**（核心） | 用了全局均速 |
| B.3 移动窗 30m → 点数 | **是** | 用了全局均速 |
| B.4 50m 分段取 max | **否** | 正确 |
| B.5 归一化 | **否** | 正确 |

---

## 四、TII vs RCI 本质区别

两者计算链路相同（滤波 → 移动RMS → 50m归一化），但有两个关键差异：

| | TII | RCI |
|--|-----|-----|
| **滤波参数** | 20-500 **Hz**（固定） | 40-300 **mm**（需用 $f=v/\lambda$ 转换） |
| **移动窗长** | **3m**（捕捉单点冲击） | **30m**（覆盖周期性波磨） |
| **速度依赖程度** | 仅窗长换算 | 滤波截止频率 + 窗长换算（双重依赖） |

为什么窗长差 10 倍？
- TII 检测的是焊缝/接头等**单点缺陷**，冲击能量集中在几厘米内，3m 窗口足够
- RCI 检测的是**波磨**（波长 40-300mm 的连续波纹），30m 内至少包含 100-750 个完整波长周期

---

## 五、基准法的修正方案

### 修正思路

标准本身说的是"3m 窗长"和"波长 40-300mm"，**没有说用全局平均速度**——那是实现时做的简化。标准 4.2 节也要求"宜进行等速检测"，说明标准的预期工况是速度恒定的，此时简化与精确等价。

但当速度变化时，需要更严格的实现：

### 修正 1：空间移动RMS（消除 A.3/B.3 窗长的速度依赖）

原实现：用全局均速把 3m/30m 换算成固定点数，然后做等点数滑动平均。

```python
# 原: 固定点数窗口
v_ms = np.nanmean(d['spd']) / 3.6     # 全局平均速度
win = int(round(3.0 / v_ms * 5000))   # 3m -> 固定点数
rm = moving_rms(af, win)              # 等点数滑窗
```

新实现：直接用里程定义窗口边界——对每个点 $i$，找 $[\text{mil}[i]-1.5\text{m},\ \text{mil}[i]+1.5\text{m}]$ 范围内的所有点算 RMS。

```python
# 新: 按里程距离定义窗口, 速度从计算链路中消失
def moving_rms_spatial(x, mil_km, window_m):
    half_w_km = window_m / 2000.0
    cs = np.cumsum(x**2)                     # 前缀和
    lo = np.searchsorted(mil, mil - half_w)  # 窗口左边界
    hi = np.searchsorted(mil, mil + half_w)  # 窗口右边界
    return sqrt((cs[hi] - cs[lo]) / (hi - lo))
```

这样无论速度如何变化，窗口始终精确覆盖 3m（或 30m）的轨道距离。

### 修正 2：按局部速度分段滤波（消除 B.2 的速度依赖）

原实现：全段用一个 $v_{\text{mean}}$ 算截止频率，一次性滤波。

新实现：按速度变化点将数据分段，每段用段内均速计算 $f = v_{\text{local}}/\lambda$，独立滤波。段边界使用填充-截取消除滤波器瞬态效应。

```python
def bp_filter_by_speed(acc, spd, fs, wl_lo, wl_hi):
    # 按速度变化点分段
    for each segment [s, e] with constant speed:
        v_local = mean(spd[s:e])
        lo_hz = v_local / wl_hi   # 300mm -> Hz
        hi_hz = v_local / wl_lo   # 40mm  -> Hz
        # 扩展边界消除瞬态, 滤波后截取
        result[s:e] = filter(acc[s-pad:e+pad])[pad:-pad]
```

---

## 六、改进处理流程详解

除了对基准法的逐点修正，还有一套更完整的改进处理流程。该流程分为**粗定位**和**细分析**两条并行链路，从根本上解决速度依赖问题，并提供更丰富的诊断信息。

![改进处理流程图](https://blog.sysu-sam.com/images/qifanyu/qcr824-flowchart.png)

### 整体架构

输入为轴箱垂向加速度时域信号（等时间采样 5000Hz）。根据速度等级 $v_{\min}, v_{\max}$ 计算各频段的滤波截止频率，然后分两条路径处理：

- **左侧——粗定位**：快速找出异常位置，给出指标值和等级
- **右侧——细分析**：对异常区段做深入频谱诊断，给出波长、成因、置信度

---

### 粗定位路径（4步）

#### Step 1：多方式带通滤波

Q/CR 824 只规定了带通滤波，没有限定滤波器类型。改进流程提供三种可选实现：

| 滤波器 | 特点 | 适用场景 |
|--------|------|---------|
| **IIR**（Butterworth） | 实时性好，计算快 | 在线监测 |
| **FFT + IFFT** | 频域截断，过渡带最陡 | 离线精确分析 |
| **FIR** | 线性相位，无相位畸变 | 需要精确时间/空间定位 |

滤波截止频率由**速度等级**和**目标波长**共同决定：

$$H_p = \left[\frac{v_{\min}}{\lambda_{\max}},\ \frac{v_{\max}}{\lambda_{\min}}\right]$$

例如，对于短波粗糙度频段（30-200mm），在速度范围 60-90 km/h 下：

$$H_p = \left[\frac{16.7}{0.200},\ \frac{25.0}{0.030}\right] = [83.3,\ 833.3]\ \text{Hz}$$

> **改进点**：Q/CR 824 使用单一速度值设计滤波器；改进流程使用速度范围 $[v_{\min}, v_{\max}]$ 确保在速度波动时目标波长始终落在通带内。

#### Step 2：时域到空间域转换（核心步骤）

这是整个改进流程的**关键创新**。

Q/CR 824 的全部计算都在时域（等时间采样）上进行，速度变化会导致同一波长在不同位置对应不同频率。改进流程在此步完成**空间重采样**：

```
等时间信号 (5000Hz, 间距随速度变化)
      |
      | 速度 + 里程 → 先积分后空间采样
      | "积分时间准确"
      ↓
等空间信号 (固定步长, 间距恒定)
      |
      ↓
空间曲线分波长度输出 → 短/中/长三个频段
```

具体操作：
1. 利用速度和里程信息，将时域信号逐点映射到里程轴
2. 在等间距的里程网格上插值，得到等空间采样信号
3. 此后所有操作在空间域进行，**速度退出计算链路**

> 图中标注"每5000个点完成一次采样"指的是：原始 5000Hz 数据先在时域积分处理，然后按空间位置重新采样。

#### Step 3：三频段分离 + 冲击抑制

空间域信号被分为三个独立的波长频段：

| 频段 | 波长范围 | 对应病害类型 |
|------|---------|------------|
| **短波粗糙度** | 30-150mm | 钢轨波磨（短波磨） |
| **中波** | 150-300mm | 钢轨波磨（长波磨）、轨面不平顺 |
| **长波** | 300-1000mm | 焊缝凹陷、道床沉降等 |

Q/CR 824 将 40-300mm 合并为一个频段；改进流程将其细分为三段，各自独立处理。

每个频段还经过**冲击信号抑制**：

- **短波频段**：接头长度 3m 冲击抑制——检测到接头位置后，将接头区段（±1.5m）的短波响应置零或平滑，避免接头冲击被误判为波磨
- **中/长波频段**：Hampel 滤波器抑制——利用 Hampel 中值检测器识别并替换离群脉冲

> **改进点**：Q/CR 824 没有专门处理"冲击干扰波磨指标"的问题。改进流程在每个频段独立抑制冲击，使 RCI 不被 TII 型病害污染。

左侧还有一条专门的**冲击检测支路**：

```
带通滤波 FIR (5mm-10mm 物理平滑)
    ↓
轮轨接触斑 → 冲击信号 → 峰峰值
    ↓
冲击指数 / 波质量指数
    ↓
Hampel 去除冲击 → 连续性 / 多波判断
    ↓
车轮周长 3m → 轨道冲击指数
```

该支路使用更窄的带通（5-10mm 物理平滑对应轮轨接触斑尺度），专门提取冲击信号的峰峰值和冲击指数，并利用车轮周长（约 3m）判断是否为车轮扁疤引起的周期性冲击。

#### Step 4：移动有效值 + 双阈值判定

三个频段各自计算移动有效值后，使用**双阈值**进行判定：

| 阈值类型 | 含义 |
|---------|------|
| **幅值阈值** | 移动有效值超过管理值（对应 Q/CR 824 的 I/II 级） |
| **长度阈值** | 超限区段的持续长度（波磨不是一个点，而是一段） |

> **改进点**：Q/CR 824 仅用幅值阈值（TII/RCI 的 I/II 级管理值）。改进流程增加长度阈值——即使幅值刚超 I 级，如果持续数百米，也应关注；反之，如果仅在一两米内超限，可能是瞬态干扰而非真正的波磨。

最终输出：**超过阈值** → 确认病害；**疑似波磨** → 需细分析确认。

---

### 细分析路径

粗定位找到可疑位置后，细分析路径对这些区段进行深入诊断：

#### 时间域分析

```
轴箱加速度 → 带通滤波 → STFT 快速谱检查 → 特征频率提取
```

- **STFT（短时傅里叶变换）**：生成时频图，观察特征频率随里程的变化
- 如果某频率在 STFT 中表现为**水平亮带**（频率不随速度变化的空间域表现），说明是固定波长的周期性病害
- 如果表现为**竖直亮线**（宽频瞬态），说明是单点冲击

#### 空间域分析

```
小波去噪 / EMD 分解 → 反演及重采样 → 频域积分
```

- **小波去噪**：多尺度分解，分离不同波长成分
- **EMD（经验模态分解）**：自适应分解为 IMF 分量，无需预设频段
- **反演及重采样**：从加速度反演出轨面不平顺（位移），需要频域二次积分

#### 输出

| 输出项 | 说明 |
|--------|------|
| **波磨位置** | 精确到米级的波磨起止里程 |
| **波磨波长** | 从空间频谱主峰反推的波长（mm） |
| **钢轨表面剥离/高拨块** | 区分波磨与剥离等不同病害类型 |
| **置信度** | 综合幅值、长度、频谱一致性给出的可信程度评分 |

---

### 粗定位 vs 细分析 对比

| | 粗定位 | 细分析 |
|--|--------|--------|
| **目的** | 快速筛查"哪里有问题" | 深入诊断"问题是什么" |
| **速度** | 快，可实时 | 慢，离线处理 |
| **输出** | 指标值 + 等级 | 波长 + 成因 + 置信度 |
| **适用** | 日常检测、养护决策 | 专项分析、维修方案制定 |

---

## 七、三种方法完整对比

| 维度 | Q/CR 824 基准法 | 基准法修正版 | 改进处理流程 |
|------|----------------|------------|------------|
| **滤波设计** | 单一均速算截止频率 | 分段用局部速度 | 速度范围 $[v_{\min}, v_{\max}]$ + 空间域固定波长 |
| **移动窗长** | 均速换算固定点数 | 里程定义空间窗口 | 空间域固定点数（等效固定米数） |
| **频段划分** | TII: 20-500Hz; RCI: 40-300mm 合并 | 同基准法 | 三频段：30-150 / 150-300 / 300-1000mm |
| **冲击处理** | 无 | 无 | Hampel + 接头 3m 抑制 |
| **判定标准** | 仅幅值阈值 | 仅幅值阈值 | 幅值 + 长度双阈值 |
| **诊断能力** | 只给指标值和等级 | 同基准法 | 粗定位 + 细分析（波长、成因、置信度） |
| **速度依赖** | 强（滤波 + 窗长） | 弱（仅 RCI 滤波分段近似） | 无（空间域后速度退出） |

---

## 八、实测数据对比

对同一段 80km 轨道，使用修正后的基准法和空间域改进法分别计算指标。

### 8车编组（均速 67 km/h，速度极为恒定）

| 指标 | 修正前基准法 | 修正后基准法 | 改进法 |
|------|------------|------------|--------|
| TII左 归一化基数 | 1.7115 | **1.7150** | 1.7096 |
| TII左 最大值 | 6.36 | **6.24** | 6.24 |
| TII左 速度相关r | 0.088 | **0.094** | 0.094 |
| RCI左 归一化基数 | 0.8884 | **0.8881** | 0.8813 |
| RCI左 最大值 | 8.38 | **8.40** | 8.42 |

**关键发现**：修正后 TII 基准法最大值与改进法**完全一致**（6.24 vs 6.24）。说明在恒速条件下，用里程定义窗口的基准法与空间域改进法在 TII 上已等价。

### 6车编组（均速 88.5 km/h）

| 指标 | 修正前基准法 | 修正后基准法 | 改进法 |
|------|------------|------------|--------|
| TII左 最大值 | 7.48 | **7.59** | 7.48 |
| TII右 最大值 | 6.98 | **7.22** | 7.48 |
| RCI右 最大值 | 5.87 | **5.90** | 5.94 |

6car 数据速度略有波动，修正后基准法的 TII 检测灵敏度有所提升（右轴从 6.98 提高到 7.22），更接近改进法的 7.48。

---

## 九、总结

### 速度依赖性全景图

```
原始加速度 (等时间采样 5000Hz)
    |
    +--- Q/CR 824 基准法 ------+--- 改进处理流程 --------+
    |                          |                         |
    |  A.2 带通 20-500Hz       |  Step1: 多方式带通滤波   |
    |    [不用速度]            |    (IIR/FFT/FIR可选)    |
    |                          |                         |
    |  B.2 波长->Hz            |  Step2: 空间域转换       |
    |    [需要速度] ←修正:分段  |    速度彻底退出          |
    |                          |                         |
    |  A.3/B.3 窗长->点数      |  Step3: 三频段分离       |
    |    [需要速度] ←修正:里程窗 |    + Hampel冲击抑制      |
    |                          |                         |
    |  A.5/B.5 归一化          |  Step4: 移动有效值       |
    |    [不用速度,标准要求全局] |    + 幅值/长度双阈值     |
    |                          |                         |
    |  输出: 指标值 + I/II级    |  输出: 指标 + 波长 +     |
    |                          |        成因 + 置信度     |
    +--------------------------+-------------------------+
```

### 核心结论

1. **标准本身没有要求用全局平均速度**，那是实现时的简化。标准的"3m 窗长"就是 3m 空间距离，"波长 40-300mm"就是固定的空间波长。

2. **TII 的速度依赖较轻**：仅窗长换算依赖速度，滤波本身（20-500Hz）不依赖。

3. **RCI 的速度依赖较重**：滤波截止频率和窗长都依赖速度，是双重依赖。

4. **改进处理流程的核心是空间域转换**：等时间信号重采样为等空间信号后，所有参数都以空间单位表达（cpm、m、mm），速度彻底退出计算链路。

5. **改进流程额外提供**：三频段细分、Hampel 冲击抑制、双阈值判定、STFT/小波/EMD 细分析，以及置信度评估——这些都是 Q/CR 824 标准中未涵盖的实用增强。

6. 对于恒速检测数据（如本文的 67/88 km/h），基准法修正版与改进法等价；**改进法的优势在变速工况下才充分体现**。