北京交通大学校徽 北京交通大学

面向堤防管涌监测的
DAS 识别方法研究

专业硕士学位论文答辩

答辩人:奚洋 专业:软件工程 学校:北京交通大学 时间:2026 年 5 月
目录
北京交通大学校徽 北京交通大学

研究背景与问题提出

管涌风险、DAS 连续监测价值、多通道时序问题、研究任务边界。

数据来源与样本可信性

现场数据、两个采集日、通道映射、标签依据、时间块防泄漏。

数据处理与模型方法

Diff 主输入、小波去噪、三频段拼接、多尺度卷积、Transformer。

实验结果与空间解释

整体识别性能、加噪鲁棒性、管涌空间定位、通道敏感性边界。

结论、不足与展望

总结研究链条、说明结果边界,提出多点位、多日期和在线验证方向。

面向堤防管涌监测的 DAS 识别方法研究 02 / 20
一、研究背景与问题提出
北京交通大学校徽 北京交通大学

本文研究问题由真实堤防 DAS 数据驱动

论文技术路线图
论文技术路线:从真实 DAS 数据到识别与管涌空间定位
边界:管涌存在来自现场先验,模型任务是 DAS 响应识别。 03 / 20
一、研究背景与问题提出
北京交通大学校徽 北京交通大学

管涌风险隐蔽,监测目标是及时识别异常响应

堤防管涌形成机制示意图
堤防管涌形成机制示意
问题本质 渗流带动细颗粒迁移并形成集中通道,异常往往早于明显破坏。
监测目标 关注异常响应阶段的识别,减少只依赖事后人工确认的滞后。
方法入口 DAS 提供沿线连续通道记录,为多通道时序识别提供数据基础。

问题到方法的对应关系

现场困难
数据表现
方法需求
早期响应弱
异常幅值不稳定
使用 Diff 与去噪增强变化信息
背景扰动多
频段差异明显
引入多频段输入而非单频段判断
响应具有空间连续性
相邻通道相关但不完全同形
利用多通道卷积和局部通道融合
研究入口:将现场异常监测问题转化为多通道时序识别问题。 04 / 20
一、研究背景与问题提出
北京交通大学校徽 北京交通大学

DAS 数据是空间通道 × 时间采样的多通道时序信号

Φ-OTDR 监测原理与堤防管涌场景映射
Φ-OTDR 监测原理与堤防管涌场景映射
Xi ∈ RC × L,三频段拼接后为 R3C × L,正式输入为 33 × 5000

C 为空间通道数,L 为时间窗口长度。该表示同时保留沿线空间位置和连续振动过程。

空间通道 响应沿相邻通道扩散,对应多通道输入与局部通道融合。
时间采样 异常存在持续过程,对应一维卷积与 Transformer。
频段差异 不同频段保留不同响应信息,对应三频段拼接输入。
输入规格:三频段拼接后为 33 × 5000。 05 / 20
一、研究背景与问题提出
北京交通大学校徽 北京交通大学

本文任务限定为已知管涌响应识别与识别后定位

现场先验

水下光缆 a 段约 195 m 处已知管涌存在。

DAS 采集

沿光缆获得多通道振动记录,形成空间通道 × 时间采样数据。

样本构建

根据通道范围、时间窗口和标签规则构造二分类样本。

阶段一:管涌响应识别 判断样本是否属于已知管涌的 DAS 振动响应,输出管涌概率。
阶段二:管涌空间定位 对识别后的正类事件估计通道级响应中心,并进行辅助解释。

模型识别

MS 1-D CNN-Transformer 建模局部变化、通道关联和时序依赖。

定位估计

基于空间能量与多频段融合,给出通道级中心估计。

辅助解释

基于输入梯度的通道敏感性分析用于观察模型关注区域。

边界:本文是“先识别、后定位”的两阶段任务。 06 / 20
二、数据来源与样本可信性
北京交通大学校徽 北京交通大学

数据来自水下光缆 a 段两个采集日的真实 DAS 记录

团洲垸现场位置与水下光缆布设示意
现场位置与水下光缆 a 段布设
2024-07-07 水下光缆 a 段现场测量图
水下光缆 a 段现场测量
采集日目标区响应特征样本用途
2024-07-07ch125--135强响应数据训练、验证与辅助参考测试
2024-07-08ch55--65弱响应数据训练、验证与主盲测集
统一采集口径水下光缆 a 段原始采样率 5 000 Hz单文件时长 1 s
数据来源:水下光缆 a 段两个采集日真实 DAS 记录。 07 / 20
二、数据来源与样本可信性
北京交通大学校徽 北京交通大学

标签依据来自现场点位、通道映射和目标区赋标规则

现场点位与光缆作业图
现场点位与水下光缆作业证据
通道响应分布与目标区示意
目标区与背景区通道响应分布
通道映射与目标区域示意
通道映射与目标区赋标依据
证据层解决的问题进入样本构建的作用
现场点位管涌空间基准从哪里来限定目标物理区域
通道映射两个采集日通道号为何不同把不同绝对通道对齐到同一物理点位
目标区赋标正负类样本如何界定目标区为正类,背景区与其他光纤区为负类
标签依据:现场点位、通道映射和目标区赋标规则共同约束。 08 / 20
二、数据来源与样本可信性
北京交通大学校徽 北京交通大学

时域、频域和相邻通道证据支撑样本可分性

Raw 与 Diff 频域特征对比
Raw 与 Diff 频域特征对比
目标区频域响应特征
目标区频域响应特征
时频响应特征
时频响应特征
相邻通道响应特征
相邻通道响应特征

从现象到方法

特征差异不是只来自单一图形。时域异常、频段能量差异、时频局部结构和相邻通道关联共同支撑后续 Diff、三频段输入与局部通道融合设计。

证据链:特征分析解释为什么需要 Diff、频段选择和局部通道融合。 09 / 20
二、数据来源与样本可信性
北京交通大学校徽 北京交通大学

时间块隔离用于控制滑动窗口重叠造成的数据泄漏

滑动窗口切片与时间块隔离示意图
滑动窗口切片与时间块隔离示意
边界:主盲测集固定来自第二采集日弱响应隔离时间块。 10 / 20
三、数据处理与模型方法
北京交通大学校徽 北京交通大学

数据处理将连续 DAS 记录转为 33 × 5000 模型输入

DAS 信号处理流程图
Raw / Diff 到多通道样本矩阵的处理链条
去噪前后信号示例 A
去噪前后时域响应示例
去噪前后信号示例 B
弱响应条件下的去噪效果
输入口径:正式主输入为去噪后 Diff,标准化只作为对照说明。 11 / 20
三、数据处理与模型方法
北京交通大学校徽 北京交通大学

Diff、小波去噪和三频段输入分别处理弱异常、噪声和频带差异

Raw 和 Diff 时域响应对比
Raw 与 Diff 时域响应对比
Raw 和 Diff 标准化响应对比
标准化后 Raw 与 Diff 对比
小波去噪前后功率谱对比
小波去噪前后功率谱密度对比
ASNR = 10 log10[(1 / |Ωe| Σ x²(t)) / (1 / |Ωb| Σ x²(t))];ΔASNR = ASNRout - ASNRin

Ωe 表示目标事件段,Ωb 表示背景段;ASNR 衡量目标段相对背景段的能量增强。

方法对应:Diff 处理弱异常,小波处理噪声,三频段处理频带差异。 12 / 20
三、数据处理与模型方法
北京交通大学校徽 北京交通大学

主模型通过多尺度卷积、局部通道融合和 Transformer 完成识别

MS 1-D CNN-Transformer 模型整体架构图
MS 1-D CNN-Transformer 模型整体架构
输入统一为 Diff,33 × 5000。 13 / 20
三、数据处理与模型方法
北京交通大学校徽 北京交通大学

多尺度卷积和局部通道融合处理持续时间变化与相邻通道关联

多尺度卷积特征提取示意图
多尺度卷积提取不同时间尺度异常
局部通道融合示意图
局部通道融合示意
yk,t = ΣcΣτ wk,c,τxc,t+τ + bk

多通道一维卷积在时间方向提取局部模式,同时在通道方向利用目标通道及相邻通道响应。

卷积核 9 侧重短时突变;stride = 2,padding = 4,pool = 4。
卷积核 17 侧重中等时间尺度波动;padding = 8。
卷积核 33 侧重较长时间跨度趋势;padding = 16。
方法对应:多尺度处理持续时间变化,局部通道融合处理相邻通道关联。 14 / 20
三、数据处理与模型方法
北京交通大学校徽 北京交通大学

四种模型采用统一输入和训练口径,差异集中在结构设计

Transformer 编码器分类流程
Transformer 编码器后的分类流程
Attention(Q, K, V) = softmax(QKT / √dk)V

自注意力用于在时序特征序列中建立长程依赖,与 BiLSTM 的递归式建模形成对照。

1-D CNN单尺度卷积前端 + 全局平均池化 + 线性分类头。
1-D CNN-BiLSTM单尺度卷积前端 + BiLSTM + 线性分类头。
MS 1-D CNN-BiLSTM多尺度卷积分支 + 局部通道融合 + BiLSTM + 分类头。
MS 1-D CNN-Transformer多尺度卷积分支 + 局部通道融合 + Transformer 编码器 + 分类头。
控制变量:统一输入与训练口径,只比较结构差异。 15 / 20
四、实验结果与空间解释
北京交通大学校徽 北京交通大学

主盲测集和辅助参考测试集验证模型识别性能

主盲测集 PR 曲线
主盲测集 PR 曲线
辅助参考测试集 PR 曲线
辅助参考测试集 PR 曲线
测试集模型RecallF1评价口径
主盲测集MS 1-D CNN-Transformer0.99000.9940用于检验隔离时段下的识别能力
辅助参考测试集MS 1-D CNN-Transformer1.00000.9980用于观察另一类跨日参考样本表现
主盲测集四模型 F10.9723 / 0.9879 / 0.9920 / 0.9940结构增强后指标逐步提高
识别结果口径:高 F1 建立在当前数据划分和样本构建规则之上。 16 / 20
四、实验结果与空间解释
北京交通大学校徽 北京交通大学

低信噪比条件下模型鲁棒性差异扩大

加噪条件下鲁棒性对比
不同加噪条件下四种模型识别性能对比
条件MS 1-D CNN-TransformerMS 1-D CNN-BiLSTM1-D CNN说明
5 dBF1 = 0.9817F1 = 0.9734中等噪声下差距开始出现
0 dBF1 = 0.9777F1 = 0.9282信噪比下降后差距扩大
-5 dBF1 = 0.9382F1 = 0.4894F1 = 0.2876低信噪比下最能体现结构稳定性
鲁棒性结论:低信噪比下主模型 F1 下降较小,但仍需更多现场噪声场景验证。 17 / 20
四、实验结果与空间解释
北京交通大学校徽 北京交通大学

管涌空间定位给出两个采集日的通道级中心估计

管涌空间定位结果
两个采集日的多频段融合管涌空间定位结果
采集日基准中心多频段融合中心最大波动判定口径
第一采集日ch130.20ch130.731.24 通道小于 2 通道阈值
第二采集日ch66.65ch67.710.58 通道小于 2 通道阈值
定位口径:先识别、后定位;结果为通道级中心估计。 18 / 20
四、实验结果与空间解释
北京交通大学校徽 北京交通大学

通道敏感性与推理代价共同说明模型边界

基于输入梯度的通道敏感性分析
基于输入梯度的通道敏感性分析
模型复杂度与推理代价对比
模型复杂度与推理代价对比
解释口径:通道敏感性解释模型关注区域,不作为正式管涌空间定位结论。 19 / 20
五、结论、不足与展望
北京交通大学校徽 北京交通大学

本文完成真实 DAS 数据到识别与空间解释的离线验证

主要结论 1:数据链条

完成真实堤防 DAS 数据整理、两个采集日样本对齐、防泄漏样本构建和四类数据集划分。

数据链条支撑后续模型训练、主盲测评价和辅助参考测试。

主要结论 2:识别模型

MS 1-D CNN-Transformer 在主盲测集上 F1 = 0.9940,辅助参考测试集 F1 = 0.9980

-5 dB 条件下主模型 F1 = 0.9382,体现低信噪比下的稳定性。

主要结论 3:空间解释

识别后管涌空间定位得到 ch130.73 和 ch67.71 两个通道级中心估计。

基于输入梯度的通道敏感性分析用于辅助解释,不替代定位结果。

不足 1:数据范围

当前验证集中于同一物理管涌点的跨日样本重构场景。

不能直接推广为完全陌生场景泛化

不足 2:定位真值

定位结果仍是通道级中心估计,缺少独立物理测量或多源传感真值校核。

后续需要与现场巡检、水位和渗流观测联合验证。

展望

扩展多点位、多日期和不同水文条件下的 DAS 数据。

进一步开展在线滑窗、告警阈值和长期连续运行验证。

研究总结

真实堤防 DAS 数据、可解释的数据集构建、深度学习识别和识别后的空间定位构成可核对的研究链条。

结束语

请各位老师批评指正。