市域铁路对埋地管道感性耦合干扰建模及仿真

解绍锋; 李卫兰; 黄大锐; 钟帆

doi:10.3969/j.issn.0258-2724.20220652

市域铁路对埋地管道感性耦合干扰建模及仿真

doi: 10.3969/j.issn.0258-2724.20220652

西南交通大学电气工程学院，四川成都 611756

基金项目: 国家自然科学基金项目（51877182）；四川省科技厅科技项目（2021YJ0028）

详细信息

作者简介:
解绍锋（1976—），男，教授，博士生导师，研究方向为牵引供电系统和电能质量，E-mail：sfxie@home.swjtu.edu.cn

中图分类号: U228
计量
- 文章访问数: 161
- HTML全文浏览量: 60
- PDF下载量: 50
- 被引次数: 4
出版历程
- 收稿日期: 2022-09-28
- 修回日期: 2023-05-04
- 网络出版日期: 2024-09-14
- 刊出日期: 2023-05-08

Modeling and Simulation of Inductive Coupling Interference from Suburban Railways to Buried Pipelines

School of Electrical Engineering, Southwest Jiaotong University, Chengdu 611756, China

摘要

摘要:
随着市域铁路和油气管道的快速发展，铁路与埋地管道平行或交叉铺设的情况已不可避免. 为分析市域铁路牵引供电系统对油气管道安全运行的影响，首先，建立铁路交流牵引供电系统与邻近埋地管道的数学模型，并利用CDEGS软件进行仿真验证；然后，采用等效距离方法将斜接近和平行接近2种情况进行统一，探究牵引供电系统对邻近埋地管道的影响机理，分析土壤电阻率、导线与埋地管道间接近距离、机车负荷电流、铁路与埋地管道间并行长度、管道涂层电阻率、动车组电流谐波次数等因素对管道沿线感应耦合电压分布的影响；最后，结合贯通地线的设置，提出4种抑制感应耦合电压的方案进行对比分析. 研究结果表明：当市域铁路牵引供电系统与埋地管道的最大间距与最小间距比值小于4.5时，采用等效距离方法的偏差在5%以内；管道上的感应耦合电压极大值随着土壤电阻率的增加而增大；接近距离由50 m增大到250 m时，感应耦合电压降幅为50.6%；机车负荷电流在200～1000 A变化时，感应耦合电压增幅较大；并行长度在2～10 km变化时，感应耦合电压从22.6 V增加到170.7 V；谐波含有率和谐波次数对感应耦合电压的影响不可忽视；在贯通地线基础上增设一根回流线的方案对抑制感应耦合电压的效果最佳.
- 市域铁路 /
- 埋地管道 /
- 感应耦合电压 /
- 斜接近 /
- 谐波
Abstract:
With the rapid development of suburban railways and oil and gas pipelines, parallel laying or cross-laying is inevitable for suburban railways and buried pipelines. In order to assess the impact of the traction power supply system of suburban railways on the safe operation of oil and gas pipelines, firstly, the mathematical model of the alternating current (AC) traction power supply system of the suburban railways and the adjacent buried pipelines was established, and the simulation based on CDEGS software was conducted. Next, the two cases of the oblique approach and parallel approach were unified by the equivalent distance method. The influence mechanism of the traction power supply system of the suburban railways on the adjacent buried pipeline was investigated, and the influence of factors including soil resistivity, distance between the conductor and the buried pipeline, locomotive load current, parallel length of railway and buried pipeline, pipeline coating resistivity, and number of current harmonics of electric multiple units on the inductive coupling voltage distribution along the pipeline was investigated. Finally, combined with the setting of a through ground line, four schemes to suppress the inductive coupling voltage were put forward for comparative analysis. The results show that the error of the equivalent distance method is within 5% when the ratio of the maximum distance to the minimum distance between the traction power supply system of the suburban railway and the buried pipeline is less than 4.5. The maximum value of the inductive coupling voltage of the pipeline increases with increasing soil resistivity. The decrease in the inductive coupling voltage is 50.6% when the distance between the suburban railway and the buried pipeline varies from 50 m to 250 m. The increase in the inductive coupling voltage rises significantly when the locomotive load current varies from 200 A to 1 000 A. The inductive coupling voltage increases from 22.6 V to 170.7 V when the parallel length varies from 2 km to 10 km. The harmonic content and the number of harmonics have a significant influence on the inductive coupling voltage. The best inductive coupling voltage suppression effect is achieved by adding a return line on the basis of the through ground line.
- suburban railway /
- buried pipelines /
- inductive coupling voltage /
- oblique approach /
- harmonics

HTML全文

汽车声学包是指具有降低车内声压尤其是空气声声压的声学部件组合，对车内声学特征分布具有重要的影响^[1]. 前围声学包作为隔离动力总成噪声的主要结构件，是整个声学包的重要部件之一，其吸隔声性能的优劣直接影响到整车NVH （noise, vibration and harshness）性能水平^[2]. 好的声学包设计不仅能够降低车内噪声水平，提升车内声学品质，而且能够促进车辆轻量化与降低成本. 因此，对汽车声学包性能提升进行研究具有重要的现实与工程意义.

早期由于技术手段限制（有限元等方法在中高频噪声计算准确性较低^[3]），汽车声学包设计和开发主要基于实际样车与声学包部件，并通过试验与对标车进行比较. 随着统计能量分析（SEA）理论的发展并达到工程化^[4]，加上大型商业软件如VAONE的推广，使得车内中高频噪声的分析得到了一定的改善. 唐荣江等^[5]针对车辆前围结构特点，依据有限元-统计能量分析（FE-SEA）模型建模原则，通过在表面创建声腔来确保能量在模型中的正确传递路径，其仿真结果与测试结果误差小于1.6 dB（A），验证了FE-SEA方法的有效性. Lee等^[6]基于SEA方法建立了某纯电动汽车整车模型，分析了对车内噪声贡献较大的多个声学包子系统，并通过正交试验设计方法对其进行吸隔声优化，优化之后的驾驶员头部声腔声压级降低了1.1 dB（A）. SEA方法在声学包工程设计中的有效性毋庸置疑，然而为构建SEA模型所需要的声学包材料参数（如孔隙率、黏性特征长度、热力特征长度等）和特性参数（如模态密度、内损耗因子、耦合损耗因子等）的准确值获取较难^[7]，导致基于SEA方法的仿真模型在简化与调参后仅能在部分特定工况与试验结果存在可比性，而在其他工况下模型置信度欠佳^[8].

虽然准确的声学包材料参数和特性参数较难获得，但是声学包吸隔声性能及其组成部件的厚度和面积等参数就相对容易获取^[9]. 在大数据与智能化的背景下，数据驱动为声学包的开发提供了更多的技术手段（如近似模型），可利用易获取的试验数据之间的模糊映射关系来代替声振动力学之间的精确推导关系^[10]，从而指导声学包的正向开发. 贺岩松等^[11]以车辆声学包材料和厚度作为设计参数，驾驶员耳旁声压级、声学包总重量与价格作为设计目标，基于Kring方法建立了设计参数与目标的近似模型，其拟合准确度达到0.95，能够较好地指导声学包优化设计. 另外，Schaefer等^[12]采用神经网络方法构建了前围声学包材料类型参数与成本和重量映射关系，在不降低原始吸隔声性能的条件下，通过粒子群算法对声学包进行了优化，优化后的声学包总重量降低了15.21%，同时，成本下降了9.7%. 然而，对于近似模型，其拟合与预测准确度可类比于仿真模型的建模与求解准确度，传统的浅层网络面对复杂问题拟合能力有限^[13]，因此，为了获得准确的声学包吸隔声性能预测结果，需要提升近似模型的预测性能. 深度信念网络（DBNs）作为深度学习的经典算法之一，以其优异的特征提取能力，在信息检索、图像识别等问题上获得很好的应用^[14-15]，因此，该方法可为声学包吸隔声性能预测提供借鉴.

由于汽车声学包匹配的非规则性以及噪声传递特性的复杂性，直接构建底层设计参数与声学包吸隔声性能的近似模型可能存在欠拟合的风险^[16]. 声学包设计参数匹配与其整体隔声性能属于弱相关关系，而对与之相关的声包部件或路径的传递损失或插入损失关联度较高，若绕过声包部件或路径的吸隔声特性直接建立近似模型则会导致设计参数与目标间关联的复杂度剧增，从而影响预测精度. 因此，直接建立声学包设计参数与吸隔声性能的映射模型欠佳，究其原因在于没有考虑车辆声振传递特性与声学包正向开发流程^[1]. 实际上，声学包性能对于驾乘人员来说是显性的，但是其成因却是隐性的，即声学包并不是一个完全意义上的独立结构，而是由众多子系统和零部件所组成^[17]. 因此，对于前围声学包系统，其性能、重量与成本的多目标优化不能仅局限于声学包系统性能（上层）与设计参数（底层）两个层面，而需对隐藏于系统之下的子系统性能进行匹配分析. 因此，本文将从噪声传递关系与层级目标分解角度出发，提出前围声学包吸隔声性能、重量与成本多目标预测与优化方法，并将其应用于具体车型的前围声学包进行验证，可为今后声学包的分析与优化研究提供参考.

1. 研究方法

1.1 基于层级分解的前围声学包多目标优化模型

根据声学包噪声传递关系与正向开发流程，前围声学包系统可划分为系统级、子系统级与零部件级，据此可建立前围声学包多层级目标分解架构如图1所示. 其中，顶层设计目标为前围系统的吸隔声性能、重量和成本，而与之相连的相邻下层级设计变量为各子系统的性能指标（前围钣金传递损失、内前围隔音垫插入损失、外前围隔音垫插入损失、前风挡传递损失）与工艺指标（前围覆盖率、前围泄漏量），其中工艺指标与声学包性能、重量和成本都有关系^[18]，与子系统级相连的底层设计变量为零部件的材料组成、分块厚度及分块面积. 图1所示的多层级目标分解架构具备如下两个优点：1）按照层级进行设计目标与设计变量的梯度分解有利于降低设计难度，并且符合声学包的开发流程；2）中间层级（子系统层）的存在降低了整个优化模型的非线性度与复杂度，从而能够更容易地避开局部最优解. 另外，图1所示的多层级目标分解架构的层级之间并非全连接关系，这主要取决于层级之间的整体相关性与层级之内的局部协调性，反映在前围声学包正向设计流程与声振传递关系中，从而可实现各层级设计变量间的解耦或部分解耦，进一步降低了复杂问题的非线性，有利于提升分析与计算的准确度.

图 1 前围声学包多层级目标分解架构

Figure 1. Multi-level goals decomposition architecture of firewall acoustic package

下载: 全尺寸图片幻灯片

依据图1所示的分解架构，将其进行数值化建模，建立声学包多目标的预测模型与优化模型如下：

$\left\{ { \begin{aligned} \text{min}&\quad\left\{{y}_{k}^{[L]}\right\}\text{，}\quad k=1,2,\cdots,K,\;L=1,2,\cdots,{L}_{0},\\ &\quad {y}_{k}^{[L]}={f}_{k}^{[L]}({x}_{1}^{[L\text{ + 1]}},{x}_{2}^{[L\text{ + 1]}},\cdots,{x}_{R}^{[L + 1]}),\\ {\rm{s.t.}}\;\; &\quad{c}_{m}{}^{[L]}({x}^{[L]})=0,\quad m=1,2,\cdots,M,\\ &\quad {g}_{n}{}^{[L]}({x}^{[L]})\leqslant 0,\quad n=1,2,\cdots,N,\\ \mathrm{var}.&\quad{x}_{r}^{[L\text{+1]}}\in {X}_{r}^{[L\text{+1]}},\quad r=1,2,\cdots,\\ &\quad R_1,\cdots, R_2,\cdots,R_3,\cdots，R,\\ &\quad {y}_{k}^{[L]}\in {Y}_{k}^{[L]}\text{，} \end{aligned} } \right.$

(1)

式中：L₀为层级分解架构总层数；L为层级分解架构中的拟计算层级； ${y}_{k}^{[L]}$ 为相邻上层级设计目标；K为该层级设计目标个数； ${x}_{r}^{[L+1]}$ 为相邻下层级的设计变量；R为该层级设计变量个数； ${X}_{r}^{[L + 1]}$ 为第r个设计变量在第L + 1层的可行域区间； ${Y}_{k}^{\left[L\right]}$ 为第k个设计变量在第L层的可行域区间； ${f}_{k}^{[L]}$ 为预测模型，用以描述设计变量与设计目标间的映射关系，因此， ${f}_{k}^{[L]}$ 模型的精准度对预测结果的准确性非常重要，对此本文提出SVR-DBNs （support vector regression- deep belief networks）方法作为前围声学包的预测模型，其具体计算过程将在下一小节进行介绍； ${c}_{m}^{[L]}$ 和 ${g}_{n}^{[L]}$ 为约束条件.

由于层级目标的分解，整个预测模型可看作前向传递的计算过程，而优化模型则可看作反向传递的求解过程，其中的多层级目标是与传统优化模型的不同之处. 基于该模型，即可对前围声学包的性能、重量和成本进行多目标优化设计.

1.2 SVR-DBNs方法提出

1.2.1 DBNs模型

DBNs是由多个受限玻尔兹曼机（RBM）堆叠而成的网络结构^[19]，其核心是采用逐层贪婪学习算法去优化网络的连接权重与偏置，从而提取数据特征同时进行模式识别，DBNs的拓扑结构如图2所示. RBM为DBNs的重要组成单元^[20]，是一种特殊类型的马尔可夫随机场，由一层二元随机可见单元（通常是Bernoulli或Gauss）v∈{0,1}^v和一层二元随机隐藏单元（通常是Bernoulli）h∈{0,1}^h组成. RBM不同层内单元之间的连接受到限制，即所有可见节点都连接到隐藏节点，而可见-可见或隐藏-隐藏节点内没有连接.

图 2 DBNs拓扑结构

Figure 2. DBNs topology

下载: 全尺寸图片幻灯片

在RBM模型中，不同层中各单元间的权重和偏置概率分布通过能量函数来定义. 对于Bernoulli-Bernoulli RBM，能量函数为

$E({\boldsymbol{v}},{{\boldsymbol{h}}};\theta ) = - \sum\limits_{i = 1}^V {\sum\limits_{j = 1}^H {{w_{ij}}{v_i}{h_j} - } } \sum\limits_{i = 1}^V {{a_i}{v_i} - \sum\limits_{j = 1}^H {{b_j}{h_j}} } \text{，}$

(2)

式中：θ=(w_ij, a_i, b_j)为模型参数；v_i和h_j分别为可见单元i和隐藏单元j输入参数；w_ij为v_i和h_j之间的权重；a_i 和 b_j 为偏置；V和H分别为可见和隐藏单元的数量.

给定能量函数和模型参数，可见单元和隐藏单元的联合分布为

$p({\boldsymbol{v}},{\boldsymbol{h}};\theta ) = \frac{{\exp ( - E({\boldsymbol{v}},{\boldsymbol{h}};\theta ))}}{{\displaystyle\sum\limits_v {\displaystyle\sum\limits_h {\exp ( - E({\boldsymbol{v}},{\boldsymbol{h}};\theta ))} } }} .$

(3)

因为可见-可见或隐藏-隐藏单元连接不存在，可见单元v_i和隐藏单元h_j的概率是独立的. 因此，条件分布p(h|v)和p(v|h)可以通过式（4）和式（5）进行计算.

$p({h_j} = 1|{\boldsymbol{v}};\theta ) = \sigma \left(\sum\limits_{i = 1}^V {{w_{ij}}{v_i} + {b_j}} \right),$

(4)

$p({v_i} = 1|{\boldsymbol{h}};\theta ) = \sigma \left(\sum\limits_{j = 1}^H {{w_{ij}}{h_j} + {a_i}} \right),$

(5)

式中： $\sigma \left(x\right)={1 + \mathrm{e}\mathrm{x}\mathrm{p}\;x}^{-1}$ .

RBM 模型参数θ可以通过最速梯度下降或共轭梯度下降来估计，以最大化训练数据的似然概率. w_ij、a_i和b_j的更新规则为

$\left\{ { \begin{gathered} \vartriangle {w_{ij}} = \varepsilon ({\langle {v_i}{h_j}\rangle _{{\rm{data}}}} - {\langle {v_i}{h_j}\rangle _{{\rm{model}}}}), \\ \vartriangle {a_i} = \varepsilon ({\langle {v_i}\rangle _{{\rm{data}}}} - {\langle {v_i}\rangle _{{\rm{model}}}}), \\ \vartriangle {b_j} = \varepsilon ({\langle {h_j}\rangle _{{\rm{data}}}} - {\langle {h_j}\rangle _{{\rm{model}}}}), \\ \end{gathered} } \right.$

(6)

式中：ε为学习率，范围为0～1； $\langle {{\text{•}}}\rangle_{{\rm{data}}}$ 为对训练数据的期望； $\langle {{\text{•}}}\rangle_{{\rm{model}}}$ 为对模型定义分布的相同期望，实际中通常利用比对散度的方法来求解^[21].

此外，在模型训练中为了提升计算效率，还可采用分批训练的方式^[22].

DBNs由多个RBM堆积而成，其底层将接收实值数据，第一个 RBM 模型是Gauss-Bernoulli RBM，其余的 RBM 模型是Bernoulli-Bernoulli RBM. DBNs的学习过程包括预训练和微调两个阶段. 在预训练阶段，DBNs将自下而上逐层贪婪地训练. 每个RBM 都是单独训练的，层之间的权重和层的偏差是固定的并保存以供进一步分析. 此过程可以表示为

$p({h_j}^{(l - 1)} = 1|{{\boldsymbol{h}}^{(l)}};{\theta ^{(l)}}) = \sigma \left(\sum\limits_{i = 1}^V {w_{ij}^{(l)}v_i^{(l)} + {b_j}^{(l - 1)}} \right),$

(7)

式中： l 为层数，底层单元的状态为h⁽⁰⁾.

根据式（6），Gauss-Bernoulli RBM也有类似的形式. 此学习过程将一直持续至最后一个隐藏层. 预训练阶段简化了深度架构的优化过程，并允许使用多参数来提高特征融合能力防止过度拟合. 在微调阶段，通过反向传播算法使用标记的输入数据对DBNs模型的权重和偏差将进行更新.

1.2.2 SVR-DBNs改进模型

常规的DBNs 模型是将线性回归作为顶层结构，所提取的特征与目标非线性较强的情况下可能导致预测效果欠佳. 因此，在本文中，提出了一种新的DBN回归模型SVR-DBNs，该模型使用支持向量回归（SVR）层代替DBN结构顶部的线性回归层，提升了模型应对强非线性问题的处理能力. SVR-DBNs以与一般 DBNs相同的方式执行预训练阶段和微调阶段. 不同之处在于，较低层的权重和偏差为通过反向传播来自顶部SVR层的梯度来学习，因此，应该根据倒数第二层的激活来区分SVR目标，故SVR-DBNs可适应维度更高、非线性更复杂的问题计算. 在实践中，L2-SVR 被广泛用于解决回归问题，若令目标为 f (w)并用倒数第二个激活h替换输入x；那么，L2-SVR的反向传播误差为

${\left. {\frac{{\partial f({\boldsymbol{w}})}}{{\partial {{\boldsymbol{h}}_n}}}} \right|_{k = 2}} = {\boldsymbol{w}} - 2C{\boldsymbol{w}}{t_n}(\max (1 - {{\boldsymbol{w}}^{\text{T}}}{{\boldsymbol{h}}_n}{t_n},0)),$

(8)

式中：w为SVR超平面的法向量； ${{\boldsymbol{h}}}_{n}$ 为第n个隐含层单元； $C > 0\mathrm{为}$ 常数，控制样本满足约束的程度；t_n∈{−1, + 1}.

反向传播算法的学习过程与一般的DBNs相同. 图3为SVR-DBNs的训练过程，其中 $，{e}_{\mathrm{m}\mathrm{a}\mathrm{x}}\mathrm{为}$ 人为设置的误差迭代最大值，例如10⁻⁶；l₀为设置的迭代停止数目.

图 3 SVR-DBN训练流程

Figure 3. SVR-DBN training flowchart

下载: 全尺寸图片幻灯片

1.3 基础数据获取

测试结果能够真实反映前围声学包的吸隔声效果，本文所用数据在长安汽车研究院的混响室-半消声室内进行收集，如图4所示，测试方法参考国际汽车工程师学会标准SAE J1400—2017^[23]，激励声源为总能量120 dB的白噪声. 为收集构建近似模型所需样本数，采用拉丁超立方试验设计^[11]构造样本组合，其中设计因素有前围钣金、外前围隔音垫和内前围隔音垫的厚度与对应的厚度面积占比. 根据工程经验，厚度设计水平上（下）限为原始方案数值分别上（下）浮动10%（20%），不同厚度对应面积占比设计水平为原始方案百分比分别增加（减少）10%（20%），同时有不同厚度面积占比总和为100%约束的限制，覆盖率设计水平为99.3%～99.7%，间隔0.1%，泄漏量设计水平为0.10%～0.30%，间隔0.05%，表1所示为汽车前围声学包原始方案. 通过拉丁超立方试验设计，生成了45组测试样本，外加1组原始状态，总共有46组样本. 测试中，依据标准SAE J1400—2017分别获取前围系统传递损失、前围钣金传递损失、外前围插入损失和内前围插入损失，其中部件孔洞由密封件进行密封处理.

图 4 前围声学包系统及零部件测试环境

Figure 4. Firewall acoustic package system and parts test environment

下载: 全尺寸图片幻灯片

表 1 前围声学包系统及零部件原始状态

Table 1. Original state of firewall acoustic package system

部件名称	材料名称	厚度/ mm	面积占比/%	密度/ （kg·m⁻³）	单价/ （元·m⁻³）
前围钣金	高强度钢	① 1.0	② 30	7820	39100
前围钣金	高强度钢	③ 1.2	④ 70	7820	39100
外前围隔音垫	玻璃纤维	⑤ 10.0	⑥ 20	100	5800
外前围隔音垫	玻璃纤维	⑦ 20.0	⑧ 80	100	5800
内前围隔音垫	PU 泡沫	⑨ 5.0	⑩ 35	64	7200
	PU 泡沫	⑪ 10.0	⑫ 65	64	7200
	棉毡	⑬ 3.0	⑭ 20	113	4000
	棉毡	⑮ 5.0	⑯ 80	113	4000
前围覆盖率		⑰ 99.5%			±0.1%, ±25 元
前围泄漏量		⑱ 0.20%			±0.05%, ±35 元
注：①～⑱表示构建预测模型的设计变量.

下载: 导出CSV

| 显示表格

2. 研究结果与讨论

2.1 基于SVR-DBNs的声学包预测结果

通过试验测试获取了准确的前围声学包系统与零部件传递损失与插入损失结果，在此基础上引入提出的SVR-DBNs方法，并结合多层级目标分解架构，建立有效适应于工程应用的声学包正向开发模型. 对于吸隔声性能目标，在构建第三层级至第二层级模型时，将零部件厚度与厚度对应的面积占比作为输入，将子系统传递损失或插入损失作为输出，为方便计算，选择传递损失和插入损失的1/3倍频点综合总值MTL （mean transmission loss）作为目标. 构建第二层级至第一层级模型时，将零部件的MTL、覆盖率和泄漏量作为输入，前围声学包系统的MTL作为输出. 对于重量和成本目标的建模与预测也是类似的流程. 在测试结果中随机选取35组样本作训练集，剩余11组做测试集进行学习，并取10次运行结果的平均值作为评价指标以降低随机误差. 对于多层级目标分解架构中的相邻层级，选择4层的SVR-DBNs网络（即含有2层隐含层）作为子模型进行建模，其中各子模型的输入层节点数为输入数据维度，隐含层节点数目通过网格寻优的方法在2与输入层维度范围内确定，输出层为基于高斯核函数的SVR. 建模所用软件为MATLAB R2020a，配置Intel 1165G7 CPU与16G内存，GPU为Intel Iris Xe Graphics with 96 EUs.

经过网络学习，图5所示为模型 ${f}_{1}^{\left[1\right]}$ ～ ${f}_{3}^{\left[1\right]}$ 在训练集与测试集中的计算结果，其中， ${f}_{1}^{\left[1\right]}$ ～ ${f}_{3}^{\left[1\right]}$ 模型相对误差最大值在各自的训练集与测试集中分别为 ${f}_{1}^{\left[1\right]}$ 的4.52%和 ${f}_{1}^{\left[1\right]}$ 的4.74%.

图 5 模型

${f}_{1}^{\left[1\right]}$ ～

${f}_{3}^{\left[1\right]}$ 的训练与测试结果

Figure 5. Training and test results of model

${f}_{1}^{\left[1\right]}$ −

${f}_{3}^{\left[1\right]}$

下载: 全尺寸图片幻灯片

表2统计了多层级模型中各子模型在测试集中的准确度，采用决定系数R²来评价. 由表2可知：第一层级模型的R²均高于0.975，其中 ${f}_{3}^{\left[1\right]}$ 的决定系数值达到0.986，同时第二层级模型的R²均不低于0.970，表明所提的SVR-DBNs方法在前围声学包预测中的有效性；由于重量与成本目标相对于MTL性能目标具有较好的线性度，所以重量和成本模型的拟合精度要高于MTL性能模型. 实际上，成本和重量在不考虑复杂因素如人工成本、加工成本等的条件下还可以直接使用线性累加模型，但为了所提方法的泛化性，本文仍选取SVR-DBNs建立重量和成本的预测模型.

表 2 多层级分解模型在测试集上的准确度

Table 2. Accuracy of multi-level model on test set

目标	第一～二层级（R²）	第二～三层级（R²）
MTL	0.975	0.970、0.974、0.968
重量	0.983	0.979、0.984、0.981
成本	0.986	0.988、0.978、0.980

下载: 导出CSV

| 显示表格

为验证SVR-DBNs的优越性，将常用的DBNs、SVR与BPNN模型应用于多层级目标分解架构并进行建模预测，其中DBNs的结构参数与SVR-DBNs一致但顶层采用线性回归模型，SVR采用高斯核函数并通过网格参数寻优，反向传播神经网络（BPNN）经过交叉验证择优选择模型结构，并通过遗传算法进行参数初始化. 表3所示为4种模型对前围声学包MTL、重量和成本的预测结果比较，可见：BPNN在测试集中的准确性相对最低，其R²在MTL、重量和成本上分别为0.917、0.924和0.931，同时均方误差（MSE）分别为0.771、1.267和13.457；SVR和DBNs的准确性较高，其中，DBNs的R²在MTL、重量和成本上分别为0.971、0.979和0.974，同时MSE分别为0.341，0.611和7.491. 对比可知，SVR-DBNs方法拥有最高的R²和最低的MSE，验证了所提方法的优越性.

表 3 声学包多目标优化结果

Table 3. Multi-goals optimization results of acoustic package

预测模型	MTL		重量		成本
预测模型	R²	MSE/dB	R²	MSE/kg	R²	MSE/元
BPNN	0.917	0.771	0.924	1.267	0.931	13.457
SVR	0.955	0.364	0.959	0.846	0.961	7.950
DBNs	0.971	0.341	0.979	0.611	0.974	7.491
SVR-DBNs	0.975	0.296	0.983	0.579	0.986	4.360

下载: 导出CSV

| 显示表格

2.2 优化结果与验证

基于前围声学包多层级目标分解架构与SVR-DBNs模型，对前围声学包的MTL、重量和成本进行多目标优化. 为了方便对多目标进行综合评价，引入“加权目标”，考虑到各个目标的重要程度不同，需设置不同的权重系数. 由于前围声学包主要目的是吸隔声，所以认为其MTL相较成本与重量更为重要，因此，拟定MTL、重量与成本的权重系数分别为0.50、0.25和0.25. 设计变量为前围系统各零部件的分块厚度与不同厚度的面积占比，根据工程可行性，设计空间为变量在原始值基础上±20%的区间范围，约束条件有两项：前围声学包最大厚度不大于其在原始状态的最大厚度；各零部件不同厚度对应的面积占比之和为100%，所建立的优化模型如式（9）所示. 遗传算法有收敛速度快、计算准确等优点，因此采用遗传算法求解式（9），其中遗传算法参数设置为种群大小200，最大遗传代数500，交叉概率0.7，变异概率0.01，代沟0.95. 优化的参数结果如表4所示，根据优化参数得到的多目标预测结果与实测结果如表5所示.

表 4 前围声学包多目标优化设计参数

Table 4. Multi-goals optimization design parameters of firewall acoustic package

部件名称	材料名称	厚度/mm	面积占比/%
前围钣金	高强度钢	1.0	32
前围钣金	高强度钢	1.2	68
外前围隔音垫	玻璃纤维	14.0	28
外前围隔音垫	玻璃纤维	20.0	72
内前围隔音垫	PU 泡沫	6.0	37
	PU 泡沫	10.0	63
	棉毡	4.0	28
	棉毡	5.0	72
前围覆盖率/%		99.6
前围泄露量/%		0.20

下载: 导出CSV

| 显示表格

表 5 前围声学包多目标优化结果与实测结果

Table 5. Multi-goals optimization results and measured results of firewall acoustic package

状态		MTL/dB	重量/kg	成本/元	加权目标
原始状态		45.4	21.4	281.6	1.0247
优化后	预测值	48.6	19.3	267.3	0.9588
优化后	实测值	47.9	19.5	269.2	0.9693

下载: 导出CSV

| 显示表格

$\qquad\qquad\qquad\left\{ {\begin{aligned} {\rm{min}}&\;\;\;\;\; {{Y}} = 0.4(1/y_1^{[1]})/{\xi _1}\; + \;0.3 y_2^{[1]}/{\xi _2}\; + \;0.3 y_3^{[1]}/{\xi _3},\\ & \quad y_{{k_1}}^{[1]} = f_{{k_1}}^{[1]}(x_1^{[2{\rm{]}}},x_2^{[2{\rm{]}}},\cdots,x_{{R_2}}^{[2]})\;, \quad{k_1} = 1,2,3,\;{R_{\rm{2}}} = 1,2,\cdots,{\rm{11}},\\ & \quad y_{{k_2}}^{[2]} = f_{{k_2}}^{[2]}(x_1^{[3{\rm{]}}},x_2^{[3{\rm{]}}},\cdots,x_{{R_3}}^{[3]})\;, \;y_{{k_2}}^{[2]} = x_{{k_2}}^{[2]},\quad{k_2} = 1,2,\cdots,9,\;{R_{\rm{3}}} = 1,2,\cdots,{\rm{16}},\\ {\rm{s}}{\rm{.}}\;{\rm{t}}{\rm{.}}&\;\;\;\;\;\;H = \;\displaystyle\sum\limits_{{R_3}{\rm{ = 1}},3,5,7} {\max {\rm{\{ }}x_{{R_{\rm{3}}}}^{[{\rm{3]}}},x_{{R_{\rm{3}}}{\rm{ + 1}}}^{[{\rm{3]}}}{\rm{\} }}} \leqslant {\rm{36}}{\rm{.2}},\\ &\quad {s_j} = x_{{R_{\rm{3}}}}^{[{\rm{3]}}}\;{\rm{ + }}x_{{R_{\rm{3}}}{\rm{ + 1}}}^{[{\rm{3]}}}{\rm{ = 100{\text{%}} }},\quad{R_3} = 9,11,13,15,\;j = ({R_3} - 7)/2,\\ {\mathop{\rm var}} .&\;\;\;\; x_{{R_{\rm{2}}}}^{[{\rm{2]}}} \in [0.8 x_{{R_1}}^{[{\rm{2]}}},\;\;1.2 x_{{R_1}}^{[{\rm{2]}}}],\;\;\\ &\quad x_{{R_{\rm{3}}}}^{[{\rm{3]}}} \in [0.8 x_{{R_{\rm{3}}}}^{[{\rm{3]}}},\;\;1.2 x_{{R_{\rm{3}}}}^{[{\rm{3]}}}],\; \end{aligned} } \right.$

(9)

式中： ${\xi }_{1}、{\xi }_{2}$ 与 ${\xi }_{3}$ 分别为设计目标 ${y}_{1}^{\left[1\right]}、$ ${y}_{2}^{\left[1\right]}$ 与 ${y}_{3}^{\left[1\right]}$ 的归一化常数，一般取数据集目标的统计均值； ${x}_{{R}_{3}}^{\left[3\right]}$ , ${R}_{3}\in \{1, 2,\cdots ,8\}$ 时为厚度参数， ${R}_{3}\in \left\{9, 10,\cdots ,16\right\}$ 时为面积占比参数，由于MTL需要越大越好，因此优化时取 ${y}_{1}^{\left[1\right]}$ 的倒数，其余参数同式（1）.

根据表5可知：多层级架构模型优化的前围声学包加权目标为0.9588（MTL、重量和成本分别为48.6 dB、19.3 kg与267.3元），较原始状态的加权目标1.0247（MTL、重量和成本分别为45.4 dB、21.4 kg与281.6元）提升了6.44%（MTL、重量和成本分别改进了7.05%、9.94%与5.07%）.

为了验证优化实际效果，对前围声学包按照优化状态进行了调整，并在混响室-半消声室内进行了验证，实测加权目标为0.9693（MTL、重量和成本分别为47.9 dB、19.5 kg与268.2元），预测结果与实测结果接近，两者相对误差为1.09%（MTL、重量和成本相对误差绝对值分别为1.44%、1.04%与0.71%），验证了所提方法的准确性与有效性.

3. 结　论

本文提出了SVR-DBNs方法对汽车前围声学包吸隔声性能、重量和成本多目标进行研究，并从噪声传递关系与层级目标分解角度出发，提出声学包多层级目标预测与优化方法，并在实际车型上进行了验证. 所提的SVR-DBNs方法对前围声学包性能、重量与成本目标预测准确度均在0.975以上，优于传统的BPNN、SVR与DBNs模型. 此外，研究中构建了前围声学包多层级架构模型，降低了整个优化模型的非线性度与复杂度，能够更容易地避开局部最优解. 基于该模型的优化结果与实测结果接近，两者加权目标相对误差为1.09%（MTL、重量和成本相对误差绝对值分别为1.44%、1.04%与0.71%），优化后的实测结果较前围声学包原始状态性能、重量和成本分别提升了5.51%、9.01%与4.40%，验证了所提方法的有效性与准确性.

致谢：感谢西南交通大学基础研究培育支持计划资助项目（2682021ZTPY068）、科技创新项目（2682022CX008，XJ2021KJZK054）的支持.

图 1 市域铁路对邻近埋地管道电磁干扰示意

Figure 1. Electromagnetic interference from suburban railways to adjacent buried pipelines

下载: 全尺寸图片幻灯片

图 2 埋地管道对地分布参数模型

Figure 2. Parametric model of buried pipeline during ground distribution

下载: 全尺寸图片幻灯片

图 3 市域铁路牵引供电系统线路与管道空间位置示意

Figure 3. Spatial location of lines of traction power supply system of suburban railways and pipelines

下载: 全尺寸图片幻灯片

图 4 市域铁路牵引供电系统与埋地管道位置示意（单位：mm）

Figure 4. Location of traction power supply system of suburban railways and buried pipelines (unit: mm)

下载: 全尺寸图片幻灯片

图 5 CDEGS仿真分析模型示意

Figure 5. CDEGS simulation model

下载: 全尺寸图片幻灯片

图 6 平行接近时数学模型与仿真模型结果对比

Figure 6. Comparison of values calculated by mathematical model with those calculated by simulation model during parallel approach

下载: 全尺寸图片幻灯片

图 7 斜接近时数学模型与仿真模型结果对比

Figure 7. Comparison of values calculated by mathematical model with those calculated by simulation model during oblique approach

下载: 全尺寸图片幻灯片

图 8 直接供电方式下市域铁路与管道空间位置

Figure 8. Spatial location of suburban railway and pipeline under direct power supply method

下载: 全尺寸图片幻灯片

图 9 土壤电阻率对管道感应耦合电压的影响

Figure 9. Influence of soil resistivity on inductive coupling voltage of pipeline

下载: 全尺寸图片幻灯片

图 10 接近距离对管道感应耦合电压的影响

Figure 10. Influence of distance on inductive coupling voltage of pipeline

下载: 全尺寸图片幻灯片

图 11 负荷电流对埋地管道感应耦合电压的影响

Figure 11. Influence of load currents on inductive coupling voltage of buried pipeline

下载: 全尺寸图片幻灯片

图 12 并行长度对埋地管道感应耦合电压的影响

Figure 12. Influence of parallel length on inductive coupling voltage of buried pipeline

下载: 全尺寸图片幻灯片

图 13 管道涂层电阻率对埋地管道感应耦合电压的影响

Figure 13. Influence of pipeline coating resistivity on inductive coupling voltage of buried pipeline

下载: 全尺寸图片幻灯片

图 14 交叉角度对埋地管道感应耦合电压的影响

Figure 14. Influence of crossover angle on inductive coupling voltage of buried pipeline

下载: 全尺寸图片幻灯片

图 15 单次谐波分量对埋地管道感应耦合电压的影响

Figure 15. Influence of single harmonic component on inductive coupling voltage of buried pipeline

下载: 全尺寸图片幻灯片

图 16 合成谐波对埋地管道感应耦合电压的影响

Figure 16. Influence of synthetic harmonics on inductive coupling voltage of buried pipeline

下载: 全尺寸图片幻灯片

图 17 基本型仿真模型横截面示意

Figure 17. Cross-section of basic simulation model

下载: 全尺寸图片幻灯片

图 18 4种抑制方案下仿真模型横截面示意

Figure 18. Cross-section of simulation model under four suppression schemes

下载: 全尺寸图片幻灯片

图 19 不同抑制方案下埋地管道感应耦合电压

Figure 19. Inductive coupling voltage of buried pipeline under different suppression schemes

下载: 全尺寸图片幻灯片

表 1 直供加回流供电方式牵引网导线规格参数

Table 1. Specifications of traction network conductors for direct power supply with return wire

导线类型	导线规格	计算半径/mm	电阻/ （Ω·km⁻¹）	等效半径/mm
接触线	CTA-120	6.45	0.148	5.03
承力索	JTMH-95	5.90	0.303	5.30
回流线	LBGLJ-185	7.85	0.146	6.11

下载: 导出CSV

表 2 钢轨规格参数

Table 2. Specifications of rail

参数	规格/（kg·m⁻¹）	截面面积/cm²	计算半径/mm	电阻/（Ω·km⁻¹）	等效半径/mm
取值	60	77.08	109.10	0.135	12.79

下载: 导出CSV

表 3 埋地管道参数

Table 3. Parameters of buried pipeline

参数	规格
管道金属相对磁导率	300
管道金属电阻率/（× 10⁻⁷ Ω·m）	1.7
管道涂层相对介电常数	5
管道涂层电阻率/（× 10⁵ Ω·m²）	1.0
管道半径/m	0.3
管道涂层厚度/m	0.003
管道埋深/m	1.5

下载: 导出CSV

表 4 负荷电流变化时埋地管道感应耦合电压最大值

Table 4. Maximum values of inductive coupling voltage of buried pipeline during load current variations

负荷电流幅值/A	数学模型/V	仿真模型/V	偏差/%
200	32.5	30.7	5.5
400	64.9	61.5	5.2
600	97.1	92.2	5.0
800	129.5	123.2	4.9
1000	162.2	154.4	4.8

下载: 导出CSV

表 5 间距比值变化时埋地管道感应耦合电压最大值

Table 5. Maximum values of inductive coupling voltage of buried pipeline with varying distance ratios

R_i	数学模型/V	仿真模型/V	偏差/%
1.5	75.9	73.5	3.2
2.0	72.2	70.2	2.8
2.5	69.7	67.8	2.7
3.0	67.3	65.2	3.1
3.5	65.6	63.2	3.7
4.0	63.8	61.1	4.2
4.5	62.4	59.4	4.8

下载: 导出CSV

表 6 某动车组谐波参数

Table 6. Harmonic parameters of electric multiple unit

主要谐波次数/次	含有率/%
3	1.38
5	2.84
21	3.00
23	4.07
25	2.40

下载: 导出CSV

参考文献(25)

[1]	曹建猷. 电气化铁道供电系统[M]. 北京:中国铁道出版社,1983.
[2]	李群湛,贺建闽. 牵引供电系统分析[M]. 3版. 成都:西南交通大学出版社,2012.
[3]	国家发展改革委,交通运输部,国家铁路局,等. 关于推动都市圈市域(郊)铁路加快发展的意见[Z]. 2020.
[4]	万红霞,李婷婷,宋东东,等. 杂散电流对埋地管道的腐蚀及排流方式的研究进展[J]. 表面技术,2021,50(4): 125-134. WAN Hongxia, LI Tingting, SONG Dongdong, et al. Research progress of stray current on corrosion and drainage method of buried pipeline[J]. Surface Technology, 2021, 50(4): 125-134.
[5]	梁毅,杜艳霞. 交流干扰和阴极保护协同作用下的腐蚀评判标准与机理研究进展[J]. 中国腐蚀与防护学报,2020,40(3): 215-222. doi: 10.11902/1005.4537.2019.053 LIANG Yi, DU Yanxia. Research progress on evaluation criteria and mechanism of corrosion under cathodic protection and AC interference[J]. Journal of Chinese Society for Corrosion and Protection, 2020, 40(3): 215-222. doi: 10.11902/1005.4537.2019.053
[6]	朱久国. 交流电气化铁路对埋地油气管道电磁干扰特性研究[D]. 成都: 西南交通大学,2018.
[7]	胡士信,路民旭,杜艳霞,等. 管道交流腐蚀的新观点[J]. 腐蚀与防护,2010,31(6): 419-424.
[8]	吴荫顺,曹备. 阴极保护和阳极保护:原理、技术及工程应用[M]. 北京: 中国石化出版社,2007: 223.
[9]	WAKELIN R G, SHELDON C. Investigation and mitigation of AC corrosion on a 300 mm diameter natural gas pipeline[C]//The 59th NACE Annual Conference. New Orleans: Corrosion, 2004: 04205.1-04205.9.
[10]	FLOYD R. Testing and mitigation of AC corrosion on 8 line: a field study[C]//The 59th NACE Annual Conference. New Orleans: Corrosion, 2004: 04210.1-04210.9.
[11]	LINHARDT P, BALL G. AC corrosion: results from laboratory investigations and from a failure analysis[C]//The 61th NACE Annual Conference. Houston: Corrosion, 2006: 06160.1-06160.9.
[12]	任增珺. 日东管道杂散电流干扰检测与防护[J]. 油气储运,2015,34(1): 111-114. doi: 10.6047/j.issn.1000-8241.2015.01.024 REN Zengjun. Detection and protection of stray current interference of Rizhao–Dongming oil pipeline[J]. Oil & Gas Storage and Transportation, 2015, 34(1): 111-114. doi: 10.6047/j.issn.1000-8241.2015.01.024
[13]	CARSON J R. Wave propagation in overhead wires with ground return[J]. Bell System Technical Journal, 1926, 5(4): 539-554 doi: 10.1002/j.1538-7305.1926.tb00122.x
[14]	高攸纲,沈远茂,石丹. 交流电气化铁道对周围电气及电子系统的阻性耦合影响[J]. 邮电设计技术,2007(3): 57-60. doi: 10.3969/j.issn.1007-3043.2007.03.013 GAO Youzang, SHEN Yuanmao, SHI Dan. Resistive coupling effects of AC electric railway to surrounding electric and electronic system[J]. Designing Techniques of Posts and Telecommunications, 2007(3): 57-60. doi: 10.3969/j.issn.1007-3043.2007.03.013
[15]	蒋俊. 交流线路正常运行时对平行敷设油气管道的电磁影响[J]. 电网技术,2008,32(2): 78-80,92. JIANG Jun. Electromagnetic influence of normally operating AC power transmission line on gas/oil pipeline parallel to transmission line[J]. Power System Technology, 2008, 32(2): 78-80,92.
[16]	汪可. 电气化铁路对油气管道的影响及防护措施[D]. 成都: 西南交通大学,2013.
[17]	YONG J, XIA B, YONG H, et al. Harmonic voltage induction on pipelines: measurement results and methods of assessment[C]//IEEE Transactions on Power Delivery. [S.l.]: IEEE, 2018: 2170-2179.
[18]	陈民武,朱久国,解绍锋,等. 牵引供电系统对埋地管道阻性耦合交流干扰建模及仿真[J]. 中国铁道科学,2018,39(2): 80-86. doi: 10.3969/j.issn.1001-4632.2018.02.10 CHEN Minwu, ZHU Jiuguo, XIE Shaofeng, et al. Modeling and simulation of resistive coupling AC interference of traction power supply system to buried pipeline[J]. China Railway Science, 2018, 39(2): 80-86. doi: 10.3969/j.issn.1001-4632.2018.02.10
[19]	盛望群. 基于CDEGS的交流电气化铁路对沿线油气管道电磁干扰影响研究[J]. 铁道科学与工程学报,2020,17(8): 2101-2108. SHENG Wangqun. Study on influence of AC electrified railway on electromagnetic interference of oil and gas pipelines along the line based on CDEGS[J]. Journal of Railway Science and Engineering, 2020, 17(8): 2101-2108.
[20]	CHARALAMBOUS C A, DEMETRIOU A, LAZARI A L, et al. Effects of electromagnetic interference on underground pipelines caused by the operation of high voltage AC traction systems: the impact of harmonics[J]. IEEE Transactions on Power Delivery, 2018, 33(6): 2664-2672. doi: 10.1109/TPWRD.2018.2803080
[21]	齐磊,崔翔,郭剑,等. 特高压交流输电线路正常运行时对输油输气管道的感性耦合计算模型[J]. 中国电机工程学报,2010,30(21): 121-126. QI Lei, CUI Xiang, GUO Jian, et al. Inductive coupling modelling of normally operating UHV AC transmission line to adjacent oil/gas pipeline[J]. Proceedings of the CSEE, 2010, 30(21): 121-126.
[22]	国家质量监督检验检疫总局. 接地系统的土壤电阻率、接地阻抗和地面电位测量导则第1部分:常规测量:GB/T 17949.1—2000[S]. 北京: 中国标准出版社,2000.
[23]	解绍锋,孙镜堤,骆冰祥,等. 高速铁路对邻近普速铁路电力电缆的干扰机理[J]. 西南交通大学学报,2021,56(1): 206-213. XIE Shaofeng, SUN Jingdi, LUO Bingxiang, et al. Mechanism of high-speed railway interference on power cables of adjacent normal-speed railway[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 206-213.
[24]	周胜军,谈萌. 基于监测数据的高铁动车组谐波特性分析[J]. 电力科学与技术学报,2018,33(3): 128-133. doi: 10.3969/j.issn.1673-9140.2018.03.020 ZHOU Shengjun, TAN Meng. Harmonic characteristics analysis of electric multiple units in high-speed railway based on the monitoring data[J]. Journal of Electric Power Science and Technology, 2018, 33(3): 128-133. doi: 10.3969/j.issn.1673-9140.2018.03.020
[25]	周秀荣. 电气化铁路对通信线路的干扰影响及防护措施[J]. 中国铁路,2007(6): 66-68. doi: 10.3969/j.issn.1001-683X.2007.06.017

施引文献

期刊类型引用(0)
其他类型引用(4)

附加材料(0)

访问统计

点击查看大图

图(19) / 表(6)

计量

文章访问数: 161
HTML全文浏览量: 60
PDF下载量: 50
被引次数: 4

1. 研究方法
1.1 基于层级分解的前围声学包多目标优化模型
1.2 SVR-DBNs方法提出
1.3 基础数据获取
2. 研究结果与讨论
2.1 基于SVR-DBNs的声学包预测结果
2.2 优化结果与验证
3. 结　论

市域铁路对埋地管道感性耦合干扰建模及仿真

doi: 10.3969/j.issn.0258-2724.20220652

作者简介: 解绍锋（1976—），男，教授，博士生导师，研究方向为牵引供电系统和电能质量，E-mail：sfxie@home.swjtu.edu.cn

计量

出版历程

Modeling and Simulation of Inductive Coupling Interference from Suburban Railways to Buried Pipelines

1. 研究方法

1.1 基于层级分解的前围声学包多目标优化模型

1.2 SVR-DBNs方法提出

1.2.1 DBNs模型

1.2.2 SVR-DBNs改进模型

1.3 基础数据获取

2. 研究结果与讨论

2.1 基于SVR-DBNs的声学包预测结果

2.2 优化结果与验证

3. 结 论

期刊类型引用(0)

其他类型引用(4)

计量

出版历程

目录

1. 研究方法

1.1 基于层级分解的前围声学包多目标优化模型

1.2 SVR-DBNs方法提出

1.3 基础数据获取

2. 研究结果与讨论

2.1 基于SVR-DBNs的声学包预测结果

2.2 优化结果与验证

3. 结 论

作者简介:
解绍锋（1976—），男，教授，博士生导师，研究方向为牵引供电系统和电能质量，E-mail：sfxie@home.swjtu.edu.cn

3. 结　论

3. 结　论