机械通风及地温对寒区隧道防冻长度的影响

陶亮亮; 曾艳华; 周小涵; 田啸宇

doi:10.3969/j.issn.0258-2724.20211002

机械通风及地温对寒区隧道防冻长度的影响

doi: 10.3969/j.issn.0258-2724.20211002

1.
西南交通大学交通隧道工程教育部重点实验室，四川成都 610031
2.
重庆大学土木工程学院，重庆 400045

基金项目: 国家自然科学基金（52104076）；四川省交通运输科技项目（2021-B-01，2021-B-03）

详细信息

作者简介:
陶亮亮（1995—），男，博士研究生，研究方向为隧道通风与防灾，E-mail：taolliang@126.com

通讯作者:
周小涵（1988—），男，副研究员，研究方向为隧道通风与防灾、隧道与地下结构稳定性，E-mail：zhouxh2008@126.com

中图分类号: U451.2
计量
- 文章访问数: 415
- HTML全文浏览量: 132
- PDF下载量: 43
- 被引次数: 0
出版历程
- 收稿日期: 2021-12-14
- 修回日期: 2022-04-19
- 网络出版日期: 2023-05-11
- 刊出日期: 2022-12-01

Effect of Mechanical Ventilation and Ground Temperature on Anti-Freezing Length of Tunnels in Cold Regions

1.
Key Laboratory of Transportation Tunnel Engineering, Ministry of Education, Southwest Jiaotong University, Chengdu 610031, China
2.
School of Civil Engineering, Chongqing University, Chongqing 400045, China

摘要

摘要:
为揭示通风参数及地温对寒区隧道冻害的影响，基于传热学理论推导了隧道围岩-衬砌-风流三维非稳态数值传热控制方程,分析了不同节点的数值传热差分方程，建立了高地温寒区隧道三维温度场数值计算模型；基于数值分析研究了机械通风速度、机械通风时间和地温对高地温寒区隧道防冻长度的影响. 结果表明：1）不考虑机械通风，孜拉山隧道的防冻长度超过1200 m；地温每升高5 ℃隧道防冻长度将减小约100 m. 2）考虑机械通风影响，当地温为10～30 ℃时，每天以2.5 m/s的速度通风2.0 h可减少防冻长度215 m；不同地温下机械通风速度每增大0.5 m/s，防冻长度减小约20 m，地温对防冻长度衰减速率影响很小；当机械通风时间小于2.0 h时，不同地温条件下增大机械通风速度对防冻长度的影响不大.
- 寒区隧道 /
- 高地温 /
- 数值分析 /
- 机械通风 /
- 防冻长度
Abstract:
In order to reveal the effect of ventilation parameters and ground temperature on freezing damage in tunnels in cold regions, the three-dimensional unsteady numerical heat transfer control equations for the tunnel surrounding rock, lining, and airflow are developed based on heat transfer theory. The numerical heat transfer difference equations of different nodes are analyzed, and a three-dimensional temperature field numerical calculation model is established for tunnels in cold regions with high ground temperatures. The effect of mechanical ventilation velocity, mechanical ventilation time, and ground temperature on the anti-freezing length of tunnels in cold regions with high ground temperatures is studied based on numerical analysis. The results show that 1) without considering mechanical ventilation, the anti-freezing length of Zilashan tunnel will exceed 1 200 m and decrease by about 100 m when the ground temperature increases by every 5 ℃. 2) By considering the effect of mechanical ventilation, when the ground temperature is 10–30 ℃, ventilation at a speed of 2.5 m/s for 2.0 h per day can reduce the anti-freezing length by 215 m; the anti-freezing length decreases by about 20 m as mechanical ventilation velocity increases every 0.5 m/s under different ground temperatures, and the effect of ground temperature on the decay rate of anti-freezing length is slight; when the mechanical ventilation time is less than 2.0 h, increasing mechanical ventilation velocity under different ground temperature conditions has little effect on anti-freezing length.
- tunnel in cold region /
- high ground temperature /
- numerical analysis /
- mechanical ventilation /
- anti-freezing length

HTML全文

数字式仪表具有读数准确方便、测量速度快、能提供数字信号输出、方便进行目视观察和数字记录等优点，已被广泛应用到变电站、工业、军事、航空等领域. 然而，在实际应用中，由于技术限制和各类其他原因，变电站中一些数字仪表的数据读取和记录始终依靠人工观察. 但是，人工观察在整个工作过程中存在手动工作量大、实时性差、效率低等问题.

对变电站数字仪表数据进行智能化监测和管理能有效增强实时性，减少安全隐患. 随着变电站智能化发展，人工采集已不能满足变电站对数据实时性的需求，对变电站数字仪表数据的实时处理已成为当前重要研究问题之一. 近年来，相关研究人员提出使用计算机视觉处理技术实现变电站数字仪表数据读取，这类方法大多先使用相机获取图像，再运用数字图像处理相关技术获取数字仪表读数，进而完成识别任务.

数字仪表读数识别采用的方法有支持向量机、神经网络、模板匹配、穿针引线法等. Cui等^[1]提出一种数字仪表数字识别方法，采用基于方向梯度直方图的线性支持向量机，其拥有较好的精确度以及鲁棒性，但不适用于大样本数据的处理；卜令正等^[2]提出一种改进的卷积神经网络（CNN）实现手写数字与数码管数字的识别，具有较快的检测速度，但是有一定精度损失；郭兰英等^[3]提出一种基于可变形CNN的数显仪表示数分割与识别方法，准确率高，但是检测速率相对较慢；陈刚等^[4]提出一种基于特征提取的数字仪表数码快速直接识别算法，能够快速识别字符，但其在复杂环境中鲁棒性较差；郭爽^[5]采用穿线法对数字仪表进行识别，有较好的识别效果，但是依赖简单的识别环境；高菊等^[6]在模板匹配的基础上，通过增加对字符特征的匹配，减少对相近字符的误判，适用于某一种类型仪表的数字识别，但对不同字体、不同倾斜程度的数字字符识别效果不佳，且模板匹配法容易遭受噪声影响.

近年来，随着计算机视觉技术发展迅速，深度学习的方法在图像目标检测领域中应用愈加广泛. 2014年提出的基于区域的CNN （RCNN）^[7]、2015年提出的全卷积网络（FCN）^[8-9]以及2015年提出的Faster RCNN^[10-11]等标志着计算机视觉领域目标检测的巨大进步. 与此同时，You Only Look Once （YOLO）系列算法也相继被提出，包括2015年YOLO-v1^[12]、2016年YOLO-v2^[13]、2017年YOLO-9000^[14]、2018年YOLO-v3^[15]、2019年YOLO-v4^[16]以及YOLO-v5^[17]. YOLO-v4作为检测速率与精确度并肩的检测算法，已经在各个行业被广泛应用，如煤矸石检测^[18]、人体检测^[19]、安全帽检测^[20]等. 相比于其他深度学习目标检测模型，YOLO-v4模型可兼顾检测性能和预测速度，能满足大部分实际检测任务的实时性需求.

尽管YOLO-v4模型的计算效率已大大提高，但对于变电站数字仪表的智能监测这一实际应用而言，其实时性和存储复杂度仍有待进一步改进. 针对深度学习模型的轻量化问题，常用的轻量级网络为MobileNet系列模型，包括2017年Howard提出的MobileNet-v1^[21]、2018年MobileNet-v2^[22]、2019年MobileNet-v3^[23]等. 该系列模型中构建的深度可分离卷积模块能在低参数量的情况下有效提取特征，但也存在非线性能力不足、模型压缩程度有限等问题.

针对上述问题，面向变电站数字仪表检测识别的实际应用需求，本文提出一种基于轻量级YOLO-v4模型的数字仪表检测识别方法，在几乎不损失模型准确度的同时，有效降低整个模型的计算和存储复杂度. 主要工作包括以下3个部分：

1）考虑到传统深度可分离卷积模块存在的非线性能力不足、压缩程度有限等问题，本文结合非线性拟合能力更强的Mish函数，提出轻量级网络模块bneck-m. 进一步，结合有效通道注意力（ECA）模块，构建ECA改进的bneck-m （ECA-bneck-m）模块，有效提高其非线性处理能力和计算效率.

2）考虑到传统YOLO-v4模型无法满足实际场景中变电站数字仪表检测识别的实时性和低存储复杂度需求，基于ECA-bneck-m模块，本文提出一种轻量级YOLO-v4模型，在几乎不损失模型检测准确度的前提下，提升模型的推理速度.

3）从鄂尔多斯变电站实际场景中采集数字仪表图像，并使用Albumentations框架对图像进行弹性变化、旋转、色彩空间变换等方式的数据增强，构建一个数据量充足、能较为全面反映变电站数字仪表可能情况的数字仪表检测识别数据集.

1. 深度学习数字仪表检测识别方法

1.1 总体研究思路

面向变电站实际场景中对数字仪表检测识别实时性和准确性的应用要求，以鄂尔多斯变电站采集的数字仪表图像作为研究对象，基于由特征提取网络以及路径聚合网络（PANet）组成的轻量级YOLO-v4网络，提出深度学习数字仪表检测识别方法，其总体框架如图1所示，主要包括数字仪表预处理、数字仪表数据集构建、用于数字仪表目标检测以及数字字符识别的轻量级YOLO-v4模型、数字仪表检测识别等部分.

图 1 深度学习数字仪表检测识别方法总体框架

Figure 1. General framework of deep learning-based method for detection and recognition of digital instruments

下载: 全尺寸图片幻灯片

1.2 基于轻量级YOLO-v4的数字仪表检测算法

类似于传统自然图像目标检测，最简单的变电站数字仪表检测方法为模板匹配.

1.2.1 模板匹配

模板匹配是在一幅图像中寻找与另一幅模板图像相似部分的技术^[24]. 该方法通过目标图像在待匹配图像中进行遍历，选择一定的匹配方式得到每个起始像素点的匹配度，匹配度最高的位置就是候选匹配位置，达到匹配查找的效果，如图2所示. 图中，T(m,n)为模板图像，A(u,v)为待识别图像，(m,n)为模板图像像素点坐标，u、v分别为待识别图像的宽、高.

图 2 模板匹配方式示意

Figure 2. Illustration of template matching

下载: 全尺寸图片幻灯片

通过将模板图像在待识别图像上进行平移，被模板覆盖的搜索区域为子图 $A_{ij}$ (m,n)，i、j为子图左上角坐标，i∈[1, u−M]，j∈[1, v−N]，M、N分别为模板图像的宽、高. T(m,n)和 $A_{ij}$ (m,n)间的误差值为D(i, j)，如式（1）所示. 误差值越小，则匹配度越高，相似性越高.

尽管这种方法计算量小，且检测速度快，但其鲁棒性较差，在复杂环境下并不适用.

$D(i,j) = {\sum\limits_{m = 1}^M {\sum\limits_{n = 1}^N {[A_{ij}(m,n) - T(m,n)]^2} } }.$

(1)

1.2.2 轻量级YOLO-v4模型

YOLO-v4的算法框架由主干特征提取网络（CSPDarknet-53）、空间金字塔池化（SPP）和PANet构成，如图3所示.

图 3 YOLO-v4模型框架

Figure 3. Framework of YOLO-v4 model

下载: 全尺寸图片幻灯片

面向变电站数字仪表检测识别，本文结合Mish函数和注意力机制，提出一种轻量级 YOLO-v4 模型，包含 4 个改进思路：

1） bneck是一种特殊的深度可分离卷积模块，能有效提取目标的特征信息，且拥有较少的参数量. 如图4所示，bneck由3 × 3深度可分离卷积（Dwise）、SE （squeeze-and-excitation）通道注意力机制和残差连接结构构成. 图中：NL代表非线性激活函数，Pool为池化处理，FC为全连接层，ReLU与Hard-σ （σ为Sigmoid函数）为激活函数，+ 表示残差连接，Dwise用于提取图像特征，SE通道注意力用于提高特征的表达能力，残差连接结构用于解决由于网络层数过深而产生过拟合的问题.

图 4 bneck模块示意

Figure 4. Structure of bneck module

下载: 全尺寸图片幻灯片

bneck模块采用H-Swish函数作为模型的激活函数，使用光滑的激活函数的模型拥有较好的性能^[25]. H-Swish、Mish函数曲线如图5所示，图中，x为函数自变量. 可以看出，H-Swish函数存在拐点，Mish函数呈现为光滑曲线. 因此，选用Mish函数作为激活函数，提出改进的深度可分离卷积模块bneck-m，提高整个模块的特征提取能力.

图 5 H-Swish函数和Mish函数曲线

Figure 5. H-Swish and Mish function curves

下载: 全尺寸图片幻灯片

基于上述分析，所提出bneck-m模块结构与图4相同，由逆残差结构、深度可分离卷积结构和SE注意力模块构成，在提高模块非线性能力的同时，自适应地从特征的通道维度学习有用特征信息，同时抑制干扰信息.

2）在1）基础上，为进一步降低bneck-m模块参数量，且提高其计算效率，结合ECA模块，将SE模块（图6）替换为ECA模块（图7），提出如图8所示的ECA-bneck-m模块. 图6中：e为输入特征，H × W × C为通道维度. SE模块对每个特征通道单独使用全局平均池化（GAP），然后使用2个FC层，最后，使用Hard-σ函数生成通道权值，将通道权重与原特征相乘，获得输出特征 $\tilde e$ . 其中，FC层的设计是为了捕捉非线性的跨通道交互，并通过降维控制其复杂性. 如所示，ECA模块先对每个通道单独全局平均池化，然后，通过考虑每个通道及其 $s$ 个邻居来捕获局部跨通道交互信息，即使用 $s$ 个1维卷积捕获局部跨通道交互信息，最后，用Sigmoid函数来生成通道权值. 相较SE模块，ECA模块的 $s$ 个1维卷积结构复杂度低，且推理速度快. 因此，ECA模块的整体性能更佳.

图 6 SE模块示意

Figure 6. Structure of SE module

下载: 全尺寸图片幻灯片

图 7 ECA模块示意

Figure 7. Structure of ECA module

下载: 全尺寸图片幻灯片

图 8 ECA-bneck-m模块示意

Figure 8. Structure of ECA-bneck-m module

下载: 全尺寸图片幻灯片

相比于bneck-m模块，ECA-bneck-m模块具有参数量小和特征提取能力强的特点，能有效降低模型参数量和存储大小.

3）传统YOLO-v4主干网络CSPDarkNet-53的参数量大，导致推理时间较长，无法满足变电站数字仪表检测识别要求，本文基于Mish函数提出ECA-bneck-m模块，构建轻量级特征提取网络，并将其作为YOLO-v4模型新的主干网络.

4）空间金字塔池化（SPP）层能提取多尺度空间特征，提升模型对于空间布局和物体变形的鲁棒性，如所示. SPP对特征分别采用4种最大池化操作处理并使用 ${k_l}$ 表示，池化得到图像特征q_l，如式（2）所示. 最后，通过级联得到输出特征q （式（3））.

图 9 SPP层结构

Figure 9. Structure of SPP layer

下载: 全尺寸图片幻灯片

$q_l = k_lo, \;\;\;l=1,2,3,4,$

(2)

$q = \sum\limits_{l = 1}^4 {q_l},$

(3)

式中： $o$ 为主干网络提取的图像特征.

由于变电站数字仪表目标不同于常规目标，本文结合实际采集的仪表图像，对SPP层进行优化，通过对不同池化组合的研究，获取最优尺度，使改进后的模型更好地适用于数字仪表的检测识别，提高模型精度.

轻量级YOLO-v4模型结构如图10所示，将其应用于变电站数字仪表的智能检测识别，从实际拍摄的图像中定位数字仪表的位置，为数字仪表读数提供更为精确的数字仪表数显区域，降低计算复杂度，同时提高整个模型的识别准确度.

图 10 轻量级YOLO-v4模型结构示意

Figure 10. Structure of lightweight YOLO-v4 model

下载: 全尺寸图片幻灯片

1.3 数字仪表检测数据集构建

由于鄂尔多斯变电站日常摄像机所拍摄的数字仪表图像受到实际环境因素的影响，图像会出现亮度较暗等情况. 因此，需要对实际采集的数字仪表图像进行预处理，本文采用伽马变换，提高原始数字仪表图像的质量.

考虑到深度学习模型需要大量的标签数据，而变电站中实际能够采集到的数字仪表图像数据有限，无法满足深度学习模型的训练需求. 因此，本文使用基于Albumentations框架的数据增强方式，实现数字仪表图像数据量级的有效提升. 该框架包括多种图像扩充的方式，例如弹性变化、旋转、高斯噪声以及色相和饱和度值偏移（HSV）等. 基于上述4种处理方法，对采集到的变电站数字仪表图像进行数据增强，并构建得到变电站数字仪表检测识别数据集，用于轻量级YOLO-v4模型的训练和学习.

1.4 基于轻量级YOLO-v4的数字仪表数显区域定位和读数识别算法

尽管1.2节已完成对变电站数字仪表的精确定位，但还需对数字仪表的数显区域和读数进行定位与识别，为工作人员智能监测变电站运行情况与异常情况提供及时的参考和辅助.

如图1所示，本文采用轻量级YOLO-v4模型，先对数字仪表的数显区域进行定位，然后定位与识别其中的数字字符，根据数字字符坐标顺序依次输出分类结果，得到最终的数字仪表读数识别结果.

2. 实验结果以及分析

2.1 数据集构建

考虑到深度学习模型需要大量标签样本，本文使用Albumentations框架对实际采集的数字仪表图像进行扩充，进而构建变电站数字仪表检测识别数据集（见表1）. 原始数字仪表图像数据集从1571幅扩充到5000幅，原始数字字符图像数据集从1201幅扩充到5000幅. 数据扩充效果如图11所示.

表 1 图像数据扩充结果

Table 1. Image data expansion results 幅

数据集	数字仪表	数字字符	总计
原数据集	1571	1201	2772
扩充数据集	5000	5000	10000

下载: 导出CSV

| 显示表格

图 11 Albumentations数据扩充效果

Figure 11. Data expansion results of Albumentations

下载: 全尺寸图片幻灯片

参考文献[26-29]，基于扩充数据集，按照9∶1的比例制作训练集及测试集. 其中，用于仪表检测的数据集总数为5000幅，4500幅为训练集，500幅为测试集；用于字符检测的数据集总数为5000幅，4500幅为训练集，500幅为测试集.

2.2 实验环境与参数设置

实验环境配置为：操作系统Windows 10，CPU AMD Ryzen 5 3600，内存16 GB，显卡NVIDIA Geforce RTX 2060 （6 GB）. 模型的训练平台为深度学习开发工具Pytorch 1.2.0，编程语言为Python.

轻量级YOLO-v4模型相关训练参数为：学习率为0.001，Batch_size为4，Epoch设置为180次.

2.3 性能评价指标

采用目标检测常用的性能评价指标衡量模型性能，包括平均精度（average precision, AP）、平均精确度均值（mean average precision, mAP）、精确率（P）、召回率（R）以及F₁值.

P是模型对某类正确预测为正样本的数量T_P与该类所有预测为正样本总数（即T_P与预测为正样本的负样本数量F_P之和）的比值，如式（4）所示.

$P = \frac{{{T_{{{\rm{P}}}}}}}{{{T_{{{\rm{P}}}}} + {F_{{{\rm{P}}}}}}}.$

(4)

R是T_P与实际为该类正样本总数（T_P与预测错误的正样本数量F_N之和）的比值，如式（5）所示.

$R = \frac{{{T_{{{\rm{P}}}}}}}{{{T_{{{\rm{P}}}}} + {F_{{{\rm{N}}}}}}}.$

(5)

F₁值（式（6））为精确率和召回率的加权调和平均，是分类问题中对性能综合评估的常用评价指标.

$F_1 = \frac{{2RP}}{{R + P}}.$

(6)

AP即精确率召回曲线的面积，AP值越高，则检测效果越好. mAP是所有类别AP的均值，作为综合评价指标.

2.4 实验结果与对比分析

2.4.1 伽马变换效果

对数字仪表图像采用1.3节基于伽马变换的预处理操作，当变换参数 $\gamma$ =0.7时，图像效果如图12所示. 数字仪表图像的亮度值得到有效提升，图像特征信息更加明显.

图 12 伽马变换预处理效果示意

Figure 12. Preprocessing results of Gamma transformation

下载: 全尺寸图片幻灯片

2.4.2 预选框的选取与设置

由于数字仪表和字符目标较为固定，使用YOLO-v4原有的预选框尺寸效果较差，需重新设定预选框尺寸. 数字仪表数据集类型和目标尺寸均较为单一，直接通过k-means聚类算法得到的聚类中心往往较为集中，进而影响整个模型训练效果.

本文采用改进的预选框尺寸聚类方法. 首先，使用k-means方法对人工标注的数字仪表数据集重新进行聚类，获取初始预选框尺寸. 在此基础上，通过线性尺度的缩放操作，对预选框尺寸进行调整，进而获得满足变电站数字仪表检测识别所需的预选框尺寸，如式（7）所示.

$\left\{ {\begin{array}{*{20}{l}} {({\tilde{\alpha}_1} ,{\tilde{\beta}_1} ) = 0.5({{\alpha}_1},{{\beta}_1})}, \\ {({\tilde{\alpha}_9} ,{\tilde{\beta}_9} ) = 3({{\alpha}_9},{{\beta}_9})}, \\ {({\tilde{\alpha}_\delta } ,{\tilde{\beta}_\delta} ) = \left( {\dfrac{{{{\alpha}_\delta} - {{\alpha}_1}}}{{{{\alpha}_9} - {{\alpha}_1}}}({\tilde{\alpha}_9} - {\tilde{\alpha}_1} ),{\tilde{\alpha}_\delta} \dfrac{{{{\beta}_\delta}}}{{{{\alpha}_\delta}}}} \right)}, \;\;\;\;2 \leqslant \delta \leqslant 8, \end{array}} \right.$

(9)

式中： $\alpha _\delta$ 和 $\beta _\delta$ 分别为k-means聚类预选框尺寸的宽和高； $\tilde \alpha_\delta$ 和 $\tilde \beta _\delta$ 分别为线性尺寸缩放后的预选框尺寸的宽和高； δ为聚类预选框的编号， $1 \leqslant \delta \leqslant 9$ .

因此，通过上述计算得到的变电站数字仪表检测识别对应的预选框尺寸如表2所示.

表 2 k-means预选框聚类结果

Table 2. k-means clustering results of prior box

模型	特征层
模型	13 × 13	26 × 26	52 × 52
仪表检测模型	（204, 149）	（84, 174）	（5, 16）
	（221, 448）	（128, 227）	（21, 36）
	（288, 144）	（174, 479）	（71, 131）
字符识别模型	（159, 191）	（94, 127）	（14, 24）
	（163, 270）	（127, 167）	（42, 62）
	（297, 876）	（131, 633）	（70, 308）

下载: 导出CSV

| 显示表格

2.4.3 SPP层的尺度优化

输入图像通过轻量级YOLO-v4模型的主干特征提取网络后，需经过SPP层，对空间尺度为13 × 13的特征进行不同尺寸的最大池化，进一步学习输入图像的多尺度属性. 先根据经验固定2组核大小为1 × 1和13 × 13，再从常用池化尺寸{3 × 3，5 × 5，7 × 7，9 × 9，11 × 11}中选取经典组合{3 × 3，5 × 5 }、{5 × 5，7 × 7}、{7 × 7，9 × 9 }、{7 × 7，11 × 11}、{5 × 5，9 × 9}. 通过5组池化尺度进行实验对比，结果如表3所示. 当尺度组合设置为{7 × 7, 9 × 9}时，整个模型的mAP值最高. 因此，用于数字仪表图像识别的轻量级YOLO-v4模型中剩余2个最大池化尺度组合设置为{7 × 7, 9 × 9}.

表 3 SPP层不同池化尺度性能对比结果

Table 3. Performance comparison of SPP layer at different pooling scales

池化尺度	mAP/%
{3 × 3, 5 × 5}	99.75
{5 × 5, 7 × 7}	99.69
{7 × 7, 9 × 9}	99.78
{7 × 7, 11 × 11}	99.74
{5 × 5, 9 × 9}	99.68

下载: 导出CSV

| 显示表格

2.4.4 模型训练参数设置

为直观展示轻量级YOLO-v4模型的训练过程，图13给出了模型的学习曲线. 当迭代次数超过180次后，整个模型趋于稳定收敛. 因此，实验中迭代次数Epoch设置为180次.

图 13 轻量级YOLO-v4模型的学习曲线

Figure 13. Learning curves of lightweight YOLO-v4 model

下载: 全尺寸图片幻灯片

2.4.5 不同深度学习模型性能对比分析

为突出YOLO-v4模型在改进过程中的轻量化变化程度，表4给出上述4种模型大小的对比结果. 由于激活函数不会影响模型大小及参数量，基于bneck和bneck-m模块的轻量级YOLO-v4模型相比原始模型（基于DarkNet-53模块），大小均从244.0 MB减少到53.8 MB，参数从63986151个缩减到14018719个. 而基于ECA-bneck-m模块轻量级YOLO-v4网络的模型大小仅为48.0 MB，参数量为12506463个，均小于原始模型，上述实验结果表明，所提出的模型在模型存储大小和参数量上更有优势.

表 4 不同网络模型大小对比结果

Table 4. Comparison results of different model sizes

网络模型	参数量/个	模型大小/MB
YOLO-v4 （DarkNet-53）	63986151	244.0
YOLO-v4 （bneck）	14018719	53.8
YOLO-v4 （bneck-m）	14018719	53.8
YOLO-v4 （ECA-bneck-m）	12506463	48.0

下载: 导出CSV

| 显示表格

基于2.1节训练集与测试集，为验证轻量级YOLO-v4 （ECA-bneck-m）模型在检测准确度和实时性上的优势，表5给出不同识别模型性能对比结果，选择Faster-RCNN、YOLO-v3、YOLO-v4、轻量级YOLO-v4 （bneck）和YOLO-v4 （bneck-m）5种模型与本文模型作对比分析. 表6给出轻量级YOLO-v4 （ECA-bneck-m）模型的每类字符识别、仪表检测与数显区域定位结果. 而图14给出数字仪表检测与读数识别结果.

表 5 不同深度学习目标检测模型对比结果

Table 5. Comparison of different deep learning detection models

网络	mAP/%	FPS/（帧·s⁻¹）
Faster-RCNN	83.88	6.0
YOLO-v3	99.64	30.0
YOLO-v4	99.80	24.0
轻量级YOLO-v4 （bneck）	99.58	33.7
轻量级YOLO-v4 （bneck-m）	99.75	35.6
轻量级YOLO-v4 （ECA-bneck-m）	99.78	36.9

下载: 导出CSV

| 显示表格

表 6 轻量级YOLO-v4（ECA-bneck-m）测试结果

Table 6. Lightweight YOLO-v4 （ECA-bneck-m） test results

类别	P/%	R/%	F₁
字符 0 识别	99.83	99.65	1.00
字符 1 识别	99.43	98.87	0.99
字符 2 识别	98.58	100.00	0.99
字符 3 识别	100.00	100.00	1.00
字符 4 识别	98.56	100.00	0.99
字符 5 识别	96.58	99.30	0.98
字符 6 识别	95.27	98.60	0.97
字符 7 识别	99.05	100.00	1.00
字符 8 识别	97.54	100.00	0.99
字符 9 识别	100.00	99.25	1.00
仪表检测	97.22	100.00	0.99
数显区域定位	98.25	100.00	0.99

下载: 导出CSV

| 显示表格

图 14 数字仪表检测与读数识别结果示意

Figure 14. Recognition results of digital instrument detection and reading

下载: 全尺寸图片幻灯片

根据表5、6所示的实验结果，YOLO-v4模型的mAP值比Faster-RCNN高15.92%，比YOLO-v3高0.16%. 从实时性角度，YOLO-v4的每秒处理帧数（frames per second，FPS）比Faster-RCNN高18.0帧/s，但略低于YOLO-v3模型（仅降低6.0帧/s）. 可以看出，相比于YOLO-v3与Faster-RCNN模型，YOLO-v4模型可以在获得较高检测准确度的同时，保证较快的推理速度.

与YOLO-v4模型相比，改进的轻量级YOLO-v4 （ECA-bneck-m）模型在仅降低0.02% mAP值的前提下，推理速度FPS从原始的24.0帧/s提升到36.9帧/s. 相比YOLO-v4 （bneck）和YOLO-v4 （bneck-m）这2种模型，轻量级YOLO-v4 （ECA-bneck-m）模型的mAP值分别提高0.20%和0.03%，同时推理速度FPS也分别提高3.2 帧/s和1.3 帧/s. 实验结果说明了ECA-bneck-m模块的有效性，同时也证明了轻量级YOLO-v4 （ECA-bneck-m）模型更符合变电站数字仪表检测识别的实际应用需求. 此外，从类别识别结果上看，使用改进ECA-bneck-m模块设计的轻量级YOLO-v4模型在每个类别的精确率均达到95%以上，而且在测试图像中存在噪声或其他复杂情况下，仍然拥有较好的准确度和鲁棒性.

3. 结　论

面向变电站数据仪表智能检测识别的实际应用需求，本文提出一种基于轻量级YOLO-v4模型的数字仪表检测识别方法. 该方法以YOLO-v4模型为基础，通过设计改进的ECA-bneck-m模块，对YOLO-v4主干特征提取网络进行轻量化设计，进而得到轻量级YOLO-v4 （ECA-bneck-m）模型，在几乎不损失准确度的前提下，有效降低模型的计算和存储复杂度. 此外，本文从鄂尔多斯变电站实际场景下采集了数字仪表图像数据，并使用Albumentations框架进行数据增强，构建了变电站数字仪表检测识别数据集. 实验结果表明，基于轻量级YOLO-v4模型的数字仪表检测识别方法能够实时、准确地实现数字仪表检测与读数识别，进而推动变电站数字仪表检测识别的智能化发展.

图 1 围岩径向节点划分示意

Figure 1. Schematic diagram of surrounding rock radial node division

下载: 全尺寸图片幻灯片

图 2 模型试验平台及测点布置（单位：cm)

Figure 2. Model test platform and measuring point layout (unit: cm)

下载: 全尺寸图片幻灯片

图 3 试验结果与数值结果对比（截面D）

Figure 3. Comparison between test results and numerical results (section D)

下载: 全尺寸图片幻灯片

图 4 孜拉山隧道纵横断面

Figure 4. Vertical section and cross section of Zilashan Tunnel

下载: 全尺寸图片幻灯片

图 5 孜拉山隧道原始地温分布

Figure 5. Original geotemperature distribution of Zilashan Tunnel

下载: 全尺寸图片幻灯片

图 6 不同机械通风速度下围岩纵向温度分布

Figure 6. Longitudinal temperature distribution of surrounding rock under different mechanical speeds and ground temperature

下载: 全尺寸图片幻灯片

图 7 不同机械通风速度下围岩温度及不同地温下防冻长度

Figure 7. Surrounding rock temperature and anti-freezing length under different mechanical ventilation speeds and geotemperature

下载: 全尺寸图片幻灯片

图 8 不同通风时间下通风速度对围岩温度及防冻长度影响

Figure 8. Influence of ventilation speed on surrounding rock temperature and anti-freezing length under different ventilation time

下载: 全尺寸图片幻灯片

表 1 昌都地区气温参数（近20年）

Table 1. Temperature parameters in Changdu area (nearly 20 years) ℃

月份	1月	2月	3月	4月	5月	6月	7月	8月	9月	10月	11月	12月
均温	−1.6	1.0	4.7	8.1	12.2	15.3	16.3	15.5	13.1	8.4	2.4	−1.5
最低均温	−9.6	−6.8	−2.4	1.4	5.3	9.2	12.6	10.0	7.3	1.8	−5.0	−9.3

下载: 导出CSV

表 2 数值计算工况

Table 2. Numerical calculation conditions

工况	机械通风风速/（m·s⁻¹）	每天机械通风时间/h	地温/℃
1～5	0	1.0，1.5，2.0， 2.5，3.0	孜拉山实际地温
6～10	2.5
11～15	3.0
16～20	3.5
21～25	4.0
26～30	0	2.0	10，15，20， 25，30
31～35	2.5
36～40	3.0
41～45	3.5
46～50	4.0

下载: 导出CSV

参考文献(21)

[1]	周幼吾, 郭东信, 邱国庆, 等. 中国冻土[M]. 北京: 科学出版社, 2000.
[2]	ILKEN Z, GUNERHAN H. An investigation about the relations between the results of heat conduction problems with and without phase change[J]. International Communications in Heat and Mass Transfer, 1996, 23(6): 899-905. doi: 10.1016/0735-1933(96)00072-3
[3]	MOTTAGHY D, RATH V. Latent heat effects in subsurface heat transport modelling and their impact on palaeotemperature reconstructions[J]. Geophysical Journal International, 2006, 164(1): 236-245. doi: 10.1111/j.1365-246X.2005.02843.x
[4]	TAKUMI K, TAKASHI M, KOUICHI F. An estimation of inner temperatures at cold region tunnel for heat insulator design[J]. Journal of Structural Engineering, 2008, 54(A): 32-38.
[5]	BRONFENBRENER L. The modelling of the freezing process in fine-grained porous media: application to the frost heave estimation[J]. Cold Regions Science and Technology, 2009, 56(2/3): 120-134.
[6]	张国柱,夏才初,殷卓. 寒区隧道轴向及径向温度分布理论解[J]. 同济大学学报(自然科学版),2010,38(8): 1117-1122,1160. doi: 10.3969/j.issn.0253-374x.2010.08.003 ZHANG Guozhu, XIA Caichu, YIN Zhuo. Analytical solution to axial and radial temperature of tunnel in cold region[J]. Journal of Tongji University (Natural Science), 2010, 38(8): 1117-1122,1160. doi: 10.3969/j.issn.0253-374x.2010.08.003
[7]	夏才初,张国柱,肖素光. 考虑衬砌和隔热层的寒区隧道温度场解析解[J]. 岩石力学与工程学报,2010,29(9): 1767-1773. XIA Caichu, ZHANG Guozhu, XIAO Suguang. Analytical solution to temperature fields of tunnel in cold region considering lining and insulation layer[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(9): 1767-1773.
[8]	冯强,蒋斌松. 多层介质寒区公路隧道保温层厚度计算的一种解析方法[J]. 岩土工程学报,2014,36(10): 1879-1887. FENG Qiang, JIANG Binsong. Analytical method for insulation layer thickness of highway tunnels with multilayer dielectric in cold regions[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1879-1887.
[9]	ZHOU X H, ZENG Y H, FAN L. Temperature field analysis of a cold-region railway tunnel considering mechanical and train-induced ventilation effects[J]. Applied Thermal Engineering, 2016, 100: 114-124.
[10]	周小涵,曾艳华,白赟,等. 基于圆形断面的隧道温度场有限差分计算模型[J]. 隧道建设,2016,36(11): 1332-1336. ZHOU Xiaohan, ZENG Yanhua, BAI Yun, et al. Finite difference calculation model for tunnel temperature field based on circular cross-section[J]. Tunnel Construction, 2016, 36(11): 1332-1336.
[11]	高焱. 寒区高速铁路隧道温度场理论与保温技术研究[D]. 成都: 西南交通大学, 2017.
[12]	高焱,朱永全,赵东平,等. 列车活塞风影响下寒区隧道温度场的变化规律[J]. 西安建筑科技大学学报(自然科学版),2017,49(1): 118-124. doi: 10.15986/j.1006-7930.2017.01.019 GAO Yan, ZHU Yongquan, ZHAO Dongping, et al. Varying rule of temperature field in the tunnel in cold region under the influence of the train piston wind[J]. Journal of Xi’an University of Architecture & Technology (Natural Science Edition), 2017, 49(1): 118-124. doi: 10.15986/j.1006-7930.2017.01.019
[13]	LU T S, ZHANG G Z, LIU S Y, et al. Numerical investigation of the temperature field and thermal insulation design of cold-region tunnels considering airflow effect[J]. Applied Thermal Engineering, 2021, 191: 116923.1-116923.15.
[14]	WU H, ZHONG Y J, XU W, et al. Experimental investigation of ground and air temperature fields of a cold-region road tunnel in NW China[J]. Advances in Civil Engineering, 2020, 2020: 4732490.1-4732490.13.
[15]	ZHAO X, ZHANG H W, LAI H P, et al. Temperature field characteristics and influencing factors on frost depth of a highway tunnel in a cold region[J]. Cold Regions Science and Technology, 2020, 179: 103141.1-103141.17.
[16]	TAO L L, REN X C, ZHAO D X, et al. Numerical study on effect of natural wind and piston wind on anti-freezing length of tunnels with high geo-temperature in cold region[J]. International Journal of Thermal Sciences, 2022, 172: 107372.1-107372.11.
[17]	JIANG H Q, NIU F J, MA Q G, et al. Numerical analysis of heat transfer between air inside and outside the tunnel caused by piston action[J]. International Journal of Thermal Sciences, 2021, 170: 107164.1-107164.15.
[18]	KANG F C, LI Y C, TANG C A. Numerical study on airflow temperature field in a high-temperature tunnel with insulation layer[J]. Applied Thermal Engineering, 2020, 179: 115654.1-115654.13.
[19]	张建荣,刘照球. 混凝土对流换热系数的风洞实验研究[J]. 土木工程学报,2006,39(9): 39-42,61. doi: 10.3321/j.issn:1000-131X.2006.09.006 ZHANG Jianrong, LIU Zhaoqiu. A study on the convective heat transfer coefficient of concrete in wind tunnel experiment[J]. China Civil Engineering Journal, 2006, 39(9): 39-42,61. doi: 10.3321/j.issn:1000-131X.2006.09.006
[20]	ZENG Y H, TAO L L, YE X Q, et al. [J]. Tunnelling and Underground Space Technology, 2020, 99: 103381.1-103381.16. ZENG Y H,TAO L L,YE X Q,et al. Temperature reduction for extra-long railway tunnel with high geotemperature by longitudinal ventilation[J]. Tunnelling and Underground Space Technology,2020,99: 103381.1-103381.16.
[21]	中铁二院工程集团有限责任公司. 铁路隧道运营通风设计规范: TB 10068—2010[S]. 北京: 中国铁道出版社, 2010.