基准排放因子更新对隧道需风量的影响研究

王旭; 王明年; 严涛; 于丽

doi:10.3969/j.issn.0258-2724.20210585

基准排放因子更新对隧道需风量的影响研究

doi: 10.3969/j.issn.0258-2724.20210585

王旭^{1, 2, 3,},
王明年^4, ,,
严涛⁴,
于丽⁴

1.
中国铁道科学研究院集团有限公司，北京 100081
2.
国铁川藏科创中心（成都）有限公司，四川成都 610213
3.
川藏铁路技术创新中心有限公司，四川成都 610213
4.
西南交通大学土木工程学院，四川成都 610031

基金项目: 中国国家铁路集团科研开发重点项目（N2020G045）；西藏自治区重点研发与转化计划（XZ201801-GB-07）

详细信息

作者简介:
王旭（1993—），男，助理研究员，博士，研究方向为隧道与地下工程，E-mail：soar1618@163.com

通讯作者:
王明年（1963—），男，教授，博士，研究方向为隧道与地下工程，E-mail：19910622@163.com

中图分类号: U459.2
计量
- 文章访问数: 330
- HTML全文浏览量: 137
- PDF下载量: 25
- 被引次数: 9
出版历程
- 收稿日期: 2021-07-20
- 修回日期: 2022-04-26
- 网络出版日期: 2023-08-07
- 刊出日期: 2022-07-06

Influence of Base Emission Factor Update on Tunnel Fresh-Air Demand

WANG Xu^{1, 2, 3
,},
WANG Mingnian^{4
, ,},
YAN Tao⁴,
YU Li⁴

1.
China Academy of Railway Sciences Corporation Limited, Beijing 100081, China
2.
Sichuan—Tibet Technology Innovation Center (Chengdu) of National Railway Co., Ltd., Chengdu 610213, China
3.
Sichuan—Tibet Railway Technology Innovation Center Co., Ltd., Chengdu 610213, China
4.
School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China

摘要

摘要:
为解决交通车辆持续更新导致的城市隧道内通风系统设计浪费和运营闲置，研究基准排放因子更新对隧道需风量的影响，并建立2种应对车辆更新的基准排放因子计算方法. 首先，通过理论分析明确需风量计算式中随时间更新的关键参数为基准排放因子；随后，根据定量分析得到基准排放因子和纵坡-车速系数更新对需风量的影响；最后，结合国外理念和实际设计经验，提出2种考虑时间更新的基准排放因子计算方法. 研究结果表明：我国道路隧道通风设计规范在制定基准排放因子时参考了世界道路协会（PIARC）通风报告，并分析得到了参考理由；与2000年相比，2021年的CO、NO_x和PM基准排放限值分别降低了81.6%、76.7%和97.9%；建立了使用机动车污染物排放限值作为基准排放因子的最不利设计方法，计算出2018年各污染物基准排放因子（CO为0.0011 m³/（辆·km），PM为0.4610 m²/（辆·km）），相比于2014年《细则》分别降低84.3%和77.0%. 建立并验证了基准排放因子计算方法，为城市道路隧道通风系统设计提供参考.
- 城市道路隧道 /
- 通风设计 /
- 需风量 /
- 基准排放因子 /
- 排放清单
Abstract:
To solve the design waste and idle operation of ventilation systems in urban tunnels caused by the continuous update of traffic vehicles, the influence of base emission factor update on tunnel fresh-air demand was studied, and two calculation methods of base emission factor in response to vehicle update were established. To begin with, the key parameters updating over time in the fresh-air demand formula were determined as base emission factors through theoretical analysis. Then, the influence of base emission factors and gradient-speed factor update on fresh-air demand was obtained based on quantitative analysis. Finally, according to foreign concepts and actual design experience, two calculation methods for base emission factors considering time were proposed. The results show that Chinese road tunnel ventilation design standards refer to the PIARC ventilation report when formulating the base emission factors, with the reasons explained. Compared with those in 2000, the base CO, NO_x, and PM emission limits in 2021 are reduced by 81.6%, 76.7%, and 97.9% respectively. The most unfavorable design method using vehicle pollutant emission limits as base emission factors is established, and the base emission factors of each pollutant in 2018 are calculated as follows: the base CO emission factor is 0.0011 m³/(veh·km), and the base PM emission factor is 0.4610 m²/(veh·km), which are reduced by 84.3% and 77.0% compared with China Standard in 2014. The research results establish and verify the calculation method of base emission factors, which provides a reference for the design of urban road tunnel ventilation systems.
- urban road tunnel /
- ventilation design /
- fresh-air demand /
- base emission factor /
- emission inventory

HTML全文

道岔是机车车辆从一股轨道转入或越过另一股轨道的线路设备，是铁路轨道的重要组成部分，也是线路上的薄弱环节^[1]. 轨距和轨底坡是高速道岔重要的设计技术条件，其参数的设置直接影响行车安全和行车品质. 60N钢轨在区间线路上的使用取得了成功，为提升列车过岔的平稳性和轮轨接触力学特性，在高速道岔区同样需要应用60N钢轨. 列车运行过程中，车轮始终处于动态磨损状态，对轮轨接触行为影响较大.

致谢：京沪高速铁路股份有限公司科技研究项目（京沪科研-2020-11）.

针对线路上的轨距和轨底坡参数，国内外学者做了大量的研究. 杜星等^[2]建立了LMD车轮和CHN60钢轨匹配的动力学模型，分析同一轨道在不同轨底坡条件下的动力学行为发现，轨底坡变化对列车直线运行时的平稳性、舒适性影响很大. 钱瑶等^[3]对比分析了不同轨底坡下60N钢轨和高速车轮LMA、XP55、S1002G匹配时的轮轨接触行为，匹配较优的轨底坡是1/20和1/30. 陈嵘等^[4-5]研究了我国地铁线路常用的LM型面与CHN60钢轨在不同轨距和非对称轨底坡下的轮轨接触特性，发现增大轨距和改变轨底坡可改善轮轨匹配关系. Cui等^[6]提出一种优化轮轨廓形的“正向求解法”，车辆动力学行为结果表明，轨底坡1/30、轨距1435 mm是优化后廓形的最优轨道参数. Sánchez等^[7]开发了一个严谨的测量轨距和轨底坡的程序确保测量过程的准确性. 李超等^[8]分析了转辙器区采用动态轨距加宽技术的轮轨接触作用，该技术可减轻尖轨磨耗和滚动接触疲劳. Ye等^[9]通过建立多体动力学模型，分析轨道参数对车轮磨耗的影响，发现轨道参数对车轮磨耗和脱轨安全性有较大影响. 闫正等^[10]分析了高速动车组车轮踏面和高速60N钢轨道岔断面的静态接触特性，发现适当地增大轨距和轨底坡有利于改善轮轨接触状态. 上述研究表明，改变轨距和轨底坡参数对轮轨接触行为影响较大，且现行的1/40轨底坡和1435 mm轨距往往不是最优的轨道参数，而国内的研究主要集中在区间线路或是针对CHN60钢轨道岔，因此，有必要对新型350 km/h 60N钢轨18号高速道岔的合理轨距和轨底坡展开研究.

本文根据实测LMA磨耗车轮型面和60N钢轨高速道岔关键断面，基于迹线法原理和三维非赫兹滚动接触理论，建立道岔区轮轨滚动接触模型，计算不同轨距和轨底坡下的接触几何参数和静力学指标，并与CHN60钢轨高速道岔进行对比.

1. 道岔关键断面和磨耗车轮

1.1 60N钢轨关键断面

利用CAD导出各关键断面离散数据，将尖轨和基本轨分开，并线性插值，得到如图1所示60N钢轨18号高速道岔转辙器区钢轨模型，X为沿钢轨的纵向里程坐标，Y为钢轨横坐标，Z为钢轨竖向坐标，尖轨顶宽35 mm钢轨断面位于轮载过渡段，其接触行为较为复杂，轮轨相互作用剧烈，对研究道岔区的轮轨接触行为具有代表性，因此，本文选取尖轨顶宽35 mm关键断面进行计算. 60N与CHN60钢轨高速道岔在尖轨顶宽35 mm处钢轨廓形见图2.

图 1 60N钢轨高速道岔直尖轨模型

Figure 1. 60N rail model of high-speed turnout straight switch rail

下载: 全尺寸图片幻灯片

图 2 尖轨顶宽35 mm处钢轨廓形

Figure 2. Rail profiles at 35 mm width of top surface of switch rail

下载: 全尺寸图片幻灯片

1.2 LMA磨耗车轮

跟踪记录某线路上运行的CRH2型动车组上LMA车轮型面演变^[11]，从标准LMA磨耗车轮到运营里程15万km，再到25万km，磨耗量增大较为明显，本文取这三种磨耗车轮进行分析. 车轮凹型磨耗对轮对接触具有较大影响^[12]，随着列车运营里程的增加，凹型磨耗车轮的磨耗量大致呈线性增大，如图3所示. 对磨耗车轮廓形通过三次样条函数进行插值、拟合.

图 3 不同运营里程LMA车轮型面

Figure 3. LMA wheel profile with different operating mileages

下载: 全尺寸图片幻灯片

2. 计算原理

2.1 轮轨接触几何计算原理

利用二维迹线法原理，不考虑车轮摇头角，由最小距离法搜索轮轨接触点. 接触点位置的求解有两个等价几何条件：1）轮轨接触点处轮轨垂直距离为0，非接触点轮轨垂直距离大于0；2）轮轨接触点处轮轨的轮廓线具有相同的斜率. 本文采用条件1求解轮轨接触点，再用条件2对结果加以验证.

2.2 UIC 519标准计算等效锥度

EN 15302标准^[13]中，等效锥度的计算采用轮对周期运动的假设，这个方法称之为UIC 519标准^[14]等效锥度. 自由轮对在轨道上的运动用微分方程表示为

$\ddot y + \frac{{v}^{2}}{e{r}_{0}}\text{Δ}r=0,$

(1)

式中： $y$ 为轮对在钢轨上的横向位移； $v$ 为车辆前进速度； $e$ 为接触点跨距； ${r_0}$ 为名义滚动圆半径； $\Delta r$ 为右侧和左侧车轮滚动半径差.

最后应用Klingel公式计算等效锥度 ${\gamma _e}$ ：

${\gamma _e} = 2e{r_0}{\left(\frac{{\text{π }}}{\lambda }\right)^2}$

(2)

式中： $\lambda$ 为轮对在钢轨上做正弦运动的波长，其求解见文献[14].

2.3 轮轨滚动接触力学原理

本文在接触力学部分采用Kalker的三维非赫兹滚动接触理论及其数值程序CONTACT，这是目前为止最为完善的滚动接触理论^[15]. 该理论将轮轨接触问题转化为数学规划问题，利用Bossinesq-Cerruti公式可得轮轨滚动接触离散模型为：

$\left\{ {\begin{array}{l} \underset{{p}_{Jj}}{\mathrm{min}\;C}=\dfrac{1}{2}{p}_{Ii}{A}_{IiJj}{p}_{Jj} +\\\\ \quad [({g}_{0J}-q){p}_{J{\textit{z}}}\; + ({W}_{J\tau }-{u}_{J\tau }){p}_{J\tau }]{A}_{0},\\ {p}_{J{\textit{z}}}\geqslant 0,\;\left|{p}_{J\tau }\right|\leqslant {b}_{J},　　\forall x\in {A}_{{\rm{c}}},\end{array} } \right.$

(3)

式中： $C$ 为接触问题的余能； $I、J$ 为接触斑上离散单元的编号； $i、j、{\textit{z}}$ 为 ${X}、{Y}、{Z}$ 坐标轴的方向； $\tau$ 为切向，对应 ${{XOY}}$ 平面； ${p_{Ii}}$ 为单元 $I$ 沿方向 $i$ 的接触应力； ${p}_{Jj}、{p}_{J{\textit{z}}}、{p}_{J\tau }$ 分别为单元J上沿Y轴、Z轴及方向 $\tau$ 的接触应力； ${A_{IiJj}}$ 为影响函数，表示单元J方向j上单位力引起的单元I方向i上位移； ${g_{0J}}$ 为单元J中心处轮轨变形前接触面间的法向间隙；q为渗透量； ${W_{J\tau }}$ 为单元J中心处从前时刻 $t'$ 到当前时刻t总的刚性滑动量； ${u_{J\tau }}$ 为轮轨滚动接触过程中时刻t时单元J中心处的弹性位移差； ${A_0}$ 为单元J的面积； ${b_J}$ 为是单元J中心处Coulomb极限摩擦力；x为轮轨在X轴上坐标； ${A_{\rm{c}}}$ 为接触斑内坐标集合；

车轮和尖轨的接触为异型接触，容易产生接触疲劳现象. 本文采用基于安定图模型的表面滚动接触疲劳因子预测不同工况下的轮轨表面滚动接触疲劳伤损. 轮轨接触斑中任一点的轮轨表面滚动接触疲劳系数 $f_{\rm{I}}({\text{•}})$ 定义为

$\left\{ \begin{gathered} f_{\rm{I}}(x,y) = {f_{\rm{t}}} - \frac{k}{{{p_{\textit{z}}}(x,y)}}, \\ {f_{\rm{t}}} = \frac{{\sqrt {{{({p_x}(x,y))}^2} + {{({p_y}(x,y))}^2}} }}{{{p_{\textit{z}}}(x,y)}}, \\ \end{gathered} \right.$

(4)

式中： ${f_{\rm{t}}}$ 为牵引系数； $y、{\textit{z}}$ 分别为接触斑在 ${Y}、{Z}$ 轴上坐标值； ${p_x}({\text{•}})$ 、 ${p_y}({\text{•}})$ 、 ${p_{\textit{z}}}({\text{•}})$ 分别为轮轨纵向、横向和法向接触应力； $k$ 为材料自身的剪切屈服强度，取350 MPa.

式（4）是根据赫兹接触理论得到，对于非赫兹接触问题，依据条带法，认为条带中间位置符合赫兹假设条件，从而将式（4）的应用扩展至非赫兹接触的范围，求得接触斑内任一单元的轮轨表面滚动接触疲劳系数，轮轨表面滚动接触疲劳因子定义为接触斑内滚动接触疲劳系数的最大值.

3. 计算结果与分析

3.1 轮轨接触点

轮轨几何接触是解释轮轨接触关系的基础^[4]. 利用迹线法原理计算轮轨接触点，计算参数有：轮背距1353 mm，名义滚动圆半径460 mm，轮背到名义滚动圆处水平距离为70 mm，以Y轴正向为正，轮对横移量取−12 ~ 12 mm，横移量步长取0.5 mm. 以往研究表明，加宽轨距和改变轨底坡往往能够改善轮轨匹配关系^[3-5]，因此，对轨距1433、1435、1437、1439 mm，轨底坡1/10、1/20、1/30、1/40、1/50进行计算分析.

不同轨距和轨底坡参数下的接触点分布如图4和图5所示. 由图4可见：随着轨距的变化，轮轨接触点分布存在明显的差异；随着轨距的增加，轮轨接触点更多地分布在尖轨顶部，不易发生轮缘接触；随着轨距的增大，轮轨接触点由基本轨转移到尖轨所需要的横移量增大，也即轮载过渡延后，有利于减小尖轨受力，但会增大轮载过渡时轮轨接触点的跳跃，横向不平顺增大；在相同轨距下，随着车轮的磨耗，轮轨接触点更多地分布在基本轨，轮载过渡位置延后，横向不平顺增大.

图 4 不同轨距下轮轨接触点分布

Figure 4. Distribution of wheel-rail contact points under different rail gauges

下载: 全尺寸图片幻灯片

图 5 不同轨底坡下轮轨接触点分布

Figure 5. Distribution of wheel-rail contact points under different rail cants

下载: 全尺寸图片幻灯片

由图5可见：不同轨底坡条件下，轮轨接触点的分布存在较大差异. 当轨底坡为1/10、1/20时，发生轮缘接触需要的轮对横移量最大，但和其余轨底坡相差较小；当轨底坡为1/10、1/20时，轮载过渡延后，横向不平顺增大；不同轨底坡下，发生轮载过渡时所需要的横移量最大为1/10，其次为1/20和1/30，最后为1/40和1/50.

3.2 等效锥度

等效锥度作为轮轨接触线性化指标，被广泛用于表征轮轨接触几何特征. 计算等效锥度的方法有简化法、谐波法和UIC 519法，UIC 519法采用轮对周期运动的假设，计算更准确^[16]，本文通过UIC 519标准^[14]计算等效锥度.

在不同轨距下，磨耗车轮和60N关键断面匹配时的等效锥度见图6. 增大轨距有利于减小车轮踏面的等效锥度，从而提升列车过岔的平稳性；当轨距为1439 mm时，等效锥度基本上小于0.05，横移量相同条件下，其等效锥度是轨距1435 mm时的20%左右，极大减小了车轮等效锥度，提升了列车过岔平稳性；轨距变化对运营里程为25万km车轮的等效锥度影响较小.

图 6 不同轨距下等效锥度

Figure 6. Equivalent conicity of different rail gauges

下载: 全尺寸图片幻灯片

图7为不同轨底坡和不同磨耗车轮条件下的等效锥度. 对于标准车轮和运营里程25万km车轮，1/30、1/40、1/50轨底坡条件下的等效锥度相差较小；当轨底坡为1/10、1/20，在横移量小于6 mm时，和轨底坡为1/30、1/40、1/50相比，等效锥度普遍较大，车辆过岔平稳性较差，横移量大于8 mm时，结果相反；在车轮运营里程为15万km时，1/30轨底坡条件下的等效锥度较小，列车过岔平稳性较好.

图 7 不同轨底坡下等效锥度

Figure 7. Equivalent conicity of different rail cants

下载: 全尺寸图片幻灯片

3.3 接触斑面积

相同法向轮轨力作用下，轮轨接触斑面积越大，其接触应力越小. 轮轨接触应力是影响轮轨磨耗和接触疲劳的重要因素. 利用非赫兹接触理论，计算轮轨接触斑面积和滚动接触疲劳因子，分析轨道参数取值对轮轨静力学接触行为的影响. 不同轨道参数工况下的参数取值：轴重14 t，单侧车轮轮心施加一半轴重，剪切模量82 GPa，泊松比0.28，摩擦系数0.3，划分网格单元0.2 mm × 0.2 mm，轮对横移量取0 ~ 12 mm. 选取右轮轨作为分析对象.

不同轨距和不同磨耗车轮条件下的接触斑面积如图8所示. 当横移量小于9 mm时，轨距和接触斑面积参数呈负相关，而当轮对横移量大于9 mm时，轨距和接触斑面积参数大致呈正相关；轮对横移量小于8 mm时，不同轨距条件下的轮轨接触斑面积相差较小，但当轮对横移量大于8 mm时，轨距越大，轮轨接触斑面积普遍越大；在轮对横移量大于8 mm时，轨距加宽有利于增大轮轨接触斑面积，减小轮轨接触应力.

图 8 不同轨距下接触斑面积

Figure 8. Contact spot area under different rail gauges

下载: 全尺寸图片幻灯片

不同轨底坡和不同磨耗车轮条件下的接触斑面积如图9所示. 1/10轨底坡下的轮轨接触斑面积普遍较小，且随车轮磨耗量的增大，其接触斑面积远小于其余轨底坡条件下的，说明1/10轨底坡下的轮轨接触力学性能较差，且会随着车轮磨耗变得更差. 由标准轮轨条件下，轨底坡和接触斑面积大致呈负相关. 随车轮磨耗量增大，1/30 ~ 1/50轨底坡的接触斑面积相差较小.

图 9 不同轨底坡下接触斑面积

Figure 9. Contact spot area under different rail cants

下载: 全尺寸图片幻灯片

3.4 滚动接触疲劳因子

利用式（4）计算轮轨滚动接触疲劳因子，分析轮轨在不同轨道参数下的接触疲劳现象. 不同轨道参数下的表面滚动接触疲劳因子分布如图10、11所示. 由图10可知：在轮对横移量小于7 mm时，不同轨距下的表面滚动接触疲劳因子相差较小，在横移量大于 7 mm时，增大轨距可延缓轮轨表面进入滚动接触疲劳区；车轮运营里程达到15万km，在横移量为3 ~ 7 mm时，易发生轮轨间的两点接触，导致轮轨表面材料易进入疲劳区，从而产生疲劳破坏；增大轨距有利于减少轮轨材料出现接触疲劳现象，延长轮轨服役寿命.

图 10 不同轨距下表面滚动接触疲劳因子

Figure 10. Surface rolling contact fatigue factors under different rail gauges

下载: 全尺寸图片幻灯片

图 11 不同轨底坡下表面接触疲劳因子

Figure 11. Surface contact fatigue factors under different rail cants

下载: 全尺寸图片幻灯片

由图11可知：标准轮轨条件下，轨底坡和滚动接触疲劳因子大致呈正相关；车轮磨耗导致轮轨滚动接触疲劳因子减小，原因是轮轨过渡位置延后，有利于减少尖轨磨耗，在横移量大于8 mm时，1/10和1/20轨底坡能延缓轮轨材料进入滚动接触疲劳区，但在横移量较小时，其滚动接触疲劳因子明显较大；1/10和1/20轨底坡下的轮轨滚动接触疲劳因子普遍较大，易引起轮轨材料进入滚动接触疲劳区，降低轮轨材料的使用寿命，轨底坡为1/30、1/40、1/50时，滚动接触疲劳因子相差较小.

随着车轮磨耗加深，道岔区的轮轨接触行为变得更为复杂. 运营里程15万km车轮和60N钢轨高速道岔在轨距1439 mm条件下发生轮轨两点接触时的接触斑分布如图12所示，接触斑中箭头指向表示切向应力合力的方向，箭头长短代表合力大小，滑动区轮轨表面间发生了相对滑移，黏着区轮轨间存在滑动趋势，但没有相对滑移. 由图12（a）可见，轨距1439 mm，轮对横移量5 ~ 7 mm时，发生了轮载过渡行为，随着横移量的增大，轮载更多由尖轨承载. 由图12（b）可见：在轮对横移量5 ~ 6 mm时，尖轨上存在较大的滑动区，在轮对横移量7 mm时，基本轨上存在较大的滑动区，两点接触导致了较大滑动区的存在，易导致车轮在钢轨上空转，引起道岔钢轨的磨损，从而缩短道岔区钢轨服役寿命. 需要指出的是，法向接触应力较小时同样可能存在较大滑动区，如图12中轮对横移量为5、7 mm时所示.

图 12 轮轨接触斑分布

Figure 12. Wheel-rail contact spot distribution

下载: 全尺寸图片幻灯片

3.5 60N和CHN60钢轨高速道岔轮轨接触几何对比

为对比60N和CHN60钢轨高速道岔的轮轨接触几何行为，在轨距1439 mm和轨底坡1/30、1/40、1/50条件下，计算60N、CHN60钢轨高速道岔35 mm顶宽关键断面与三种磨耗车轮的等效锥度，如图13.

图 13 等效锥度对比

Figure 13. Comparison of equivalent conicities

下载: 全尺寸图片幻灯片

由图13可知：相同轨底坡和轮对横移量条件下，相比于CHN60钢轨，60N钢轨高速道岔区的等效锥度更小，列车过岔平稳性更优；车轮运营里程达到25万km时，当车轮横移量小于5 mm时，60N钢轨高速道岔区等效锥度普遍更小，具有较好的过岔平稳性，车轮横移量大于5 mm时，60N钢轨高速道岔区等效锥度普遍更大，轮对的对中性能更好；随车轮磨耗量增加，等效锥度整体呈增加趋势.

4. 结　论

1）轨距对轮轨接触行为影响较大. 轨距加宽有利于减少轮缘接触，较大程度减小等效锥度，提升列车过岔的平稳性；轨距加宽可减小轮对横移量大于8 mm时的轮轨接触应力和表面滚动接触疲劳因子，减少轮轨材料发生接触疲劳，延长尖轨使用寿命. 轮对横移量小于8 mm时，轨距加宽对运营里程25万km车轮的接触性能影响较小.

2）轨底坡对轮轨接触行为影响较大. 标准轮轨条件下，轨底坡和接触斑面积大致呈负相关，与滚动接触疲劳因子大致呈正相关. 轨底坡为1/10和1/20时，轮载过渡位置延后，横向不平顺增大，车轮横移量小于6 mm时，等效锥度普遍较大；1/10和1/20轨底坡下的接触斑面积普遍较小，轮轨滚动接触疲劳因子普遍较大，较易引起轮轨材料进入滚动接触疲劳区，降低轮轨材料的使用寿命，且1/10轨底坡对车轮磨耗的适应性较差. 轨底坡为1/30、1/40、1/50时，轮轨接触参数相差较小，匹配性能较优.

3）和CHN60钢轨高速道岔相比，60N钢轨的等效锥度普遍更小，列车过岔平稳性更优；车轮运营里程为25万km时，当轮对横移量小于5 mm时，60N钢轨的等效锥度普遍更小，当轮对横移量大于5 mm时，结果相反.

4）车轮磨耗易引起道岔区轮轨间的两点接触，在较小轮轨法向接触应力下，接触斑上也易出现较大滑动区，导致车轮空转，引起钢轨伤损. 随车轮磨耗量增加，轮轨间等效锥度整体呈增加趋势.

致谢：京沪高速铁路股份有限公司科技研究项目（京沪科研-2020-11）的支持.

图 1 柴油各车型不同坡度不同车速下PM基准排放因子

Figure 1. Comparison of base PM emission factors for different diesel vehicles under different gradients and speeds

下载: 全尺寸图片幻灯片

图 2 算例隧道需风量计算结果对比

Figure 2. Comparison of calculation results for tunnel fresh air demand

下载: 全尺寸图片幻灯片

表 1 3个计算式中排放产生率计算参数对比

Table 1. Comparison of calculation parameters for emission generation rate in three calculation formulas

计算式	基准排放因子	区域因子	交通量	海拔因子	纵坡-车速系数	时间因子	质量因子	非排放颗粒
式（1）	q_co=0.007/0.015	f_a	f_dLN_m	f_h	f_iv	每年 2% 衰减	f_m,j
式（2）	q_ex,j(v,i)	f_e	n_veh,j=DLα_cat,j	f_h	v,i	f_t,j	f_m,j	4.7vq_ne,j（v）/1000
式（3）	g_cat,j(v,i)	f_e	n_veh,j=DLα_cat,j	f_h	v,i	f_t,j	f_m,j	g_non-ex,j

下载: 导出CSV

表 2 我国道路隧道需风量计算规范中随时间变化参数取值对比

Table 2. Comparison of parameter values updating over time in Chinese road tunnel ventilation design standards

参数	2014年《细则》	2017年上海《标准》
基准排放因子	q_co=0.007/0.015；q_pm=2. 来源：综合参考 2004 年PIARC 通风报告和国内科研单位成果	不同车型、燃油类型、坡度和车速下 NO₂、CO 和 PM 的基准排放因子. 来源：参考 2012年 PIARC 通风报告
基准排放年	2000 年. 来源：我国车辆保有量中绝大部分机动车为 2000 年后生产的	2010 年. 来源：2012 年PIARC 通风报告排放因子基准排放年为 2010 年
时间因子	2%，最高 30 年. 来源：尾气污染物年递减率一般超过 10%，但安全考虑按 2% 取值	不同车型、燃油类型、污染物类别每 5 年一个衰减系数，最高至 2030 年. 来源：参考 2012 年 PIARC 通风报告
纵坡-车速系数	给出了车速为 10～100 km/h、坡度为−4%～4% 的修正因子. 来源：参照了1987 年和1991年 PIARC 通风报告	已集成在基准排放因子中，车速范围为 0～100 km/h，坡度为−6%～6%. 来源：参考 2012 年 PIARC通风报告

下载: 导出CSV

表 3 汽车排放标准对比（以轻型汽车为例）

Table 3. Comparison of vehicle emission standards（light-duty vehicles） g/km

排放标准	执行年	CO	NO_x	PM
欧 1/国 Ⅰ	1992 年/2001 年	2.72/2.72		0.140 0/0.140 0
欧 2/国 Ⅱ	1997 年/2004 年	2.20/2.20		0.140 0/0.140 0
欧 3/国 Ⅲ	2000 年/2007 年	2.30/2.30	0.150/0.150	0.140 0/0.140 0
欧 4/国 Ⅳ	2005 年/2010 年	1.00/1.00	0.080/0.080	0.140 0/0.140 0
欧 5/国 Ⅴ	2008 年/2018 年	1.00/1.00	0.060/0.060	0.005 0/0.004 5
欧 6/国 Ⅵ a	2014 年/2019 年	1.00/ 0.70	0.060/0.060	0.005 0/0.004 5
国 Ⅵ b	2021 年	0.50	0.035	0.003 0

下载: 导出CSV

表 4 汽油小客车基准排放因子

Table 4. Base emission factors for gasoline passenger cars g/km

排放标准	CO	NO_x	PM 2.5
国Ⅰ前	25.72	1.971	0.028
国Ⅰ	6.71	0.409	0.026
国 Ⅱ	2.52	0.324	0.011
国 Ⅲ	1.18	0.100	0.007
国Ⅳ	0.68	0.032	0.003
国Ⅴ	0.46	0.017	0.003

下载: 导出CSV

表 5 各排放标准下汽车组成比例

Table 5. Proportion of vehicles under various emission standards %

排放标准	2012PIARC 通风报告	2019PIARC 通风报告	2018 年汽车环境管理年报
欧 1 前/ 国Ⅰ前	12.37	2.12	0.10
欧 1/国Ⅰ	3.73	1.84	2.90
欧 2/国 Ⅱ	11.83	2.36	4.50
欧 3/国 Ⅲ	28.15	5.10	19.10
欧 4/国Ⅳ	35.05	14.66	42.50
欧 5/国Ⅴ	8.87	32.84	30.90
欧 6/国Ⅵ	0	41.10	0

下载: 导出CSV

表 6 2018年各污染物基准排放因子（轻型汽车）

Table 6. Base emission factors of various pollutants in 2018 （light-duty vehicles） g/km

排放阶段	保有量比例/%	CO	NO_x	PM
国Ⅰ前	0.1	1.350	0.092	0.098
国Ⅰ	2.9
国 Ⅱ	4.5
国 Ⅲ	19.1
国Ⅳ	42.5
国Ⅴ	30.9

下载: 导出CSV

参考文献(26)

[1]	张宏,姚延钢,杨晓勤. 城市道路轻型汽车行驶工况构建[J]. 西南交通大学学报,2019,54(6): 1139-1146,1154. ZHANG Hong, YAO Yangang, YANG Xiaoqin. Light-duty vehicles driving cycle construction based on urban roads[J]. Journal of Southwest Jiaotong University, 2019, 54(6): 1139-1146,1154.
[2]	中华人民共和国交通运输部. 公路隧道通风设计细则: JTG/T D70/2-02—2014[S]. 北京: 人民交通出版社, 2014.
[3]	WANG M N, WANG X, YU L, et al. Field measurements of the environmental parameter and pollutant dispersion in urban undersea road tunnel[J]. Building and Environment, 2019, 149: 100-108. doi: 10.1016/j.buildenv.2018.11.036
[4]	WANG X, WANG M N, CHEN J, et al. Analysis of calculation of fresh-air demand for road tunnel ventilation design in China[J]. Tunnelling and Underground Space Technology, 2020, 103: 103469.1-103469.8. doi: 10.1016/j.tust.2020.103469
[5]	吴德兴,任小峰,郑国平. 公路特长隧道通风与照明系统协同集约化设计探索[J]. 交通运输研究,2017,3(6): 40-45. doi: 10.16503/j.cnki.2095-9931.2017.06.006 WU Dexing, REN Xiaofeng, ZHENG Guoping. Collaborative intensive design of ventilation and lighting system of highway extra-long tunnel[J]. Transportation Standardization, 2017, 3(6): 40-45. doi: 10.16503/j.cnki.2095-9931.2017.06.006
[6]	朱凌,刘杨青. 高速公路隧道运营管理费用的分析研究[J]. 公路交通科技(应用技术版),2013(9): 219-220.
[7]	DENG Y W, CHEN C, LI Q, et al. Measurements of real-world vehicle CO and NO_x fleet average emissions in urban tunnels of two cities in China[J]. Atmospheric Environment, 2015, 122: 417-426. doi: 10.1016/j.atmosenv.2015.08.036
[8]	王明年,王旭,于丽,等. 翔安海底隧道运营环境及污染物分布规律[J]. 西南交通大学学报,2020,55(4): 695-703. WANG Mingnian, WANG Xu, YU Li, et al. Operating environment and pollutant distribution in Xiang’an undersea tunnel[J]. Journal of Southwest Jiaotong University, 2020, 55(4): 695-703.
[9]	刘洋. 西安城市公路隧道空气污染物浓度分布及通风方式研究[D]. 西安: 西安工程大学, 2018.
[10]	宁艳涛,贺肖杰. 青岛胶州湾海底公路隧道污染物浓度水平测试[J]. 建筑热能通风空调,2020,39(10): 72-75,71. doi: 10.3969/j.issn.1003-0344.2020.10.018 NING Yantao, HE Xiaojie. Field measurement of contaminant concentration in underwater road tunnel[J]. Building Energy & Environment, 2020, 39(10): 72-75,71. doi: 10.3969/j.issn.1003-0344.2020.10.018
[11]	王东伟,苟红松,戴新. 南昌红谷隧道通风及洞口污染物排放分析[J]. 暖通空调,2021,51(3): 59-65. WANG Dongwei, GOU Hongsong, DAI Xin. Analysis on ventilation and pollutant emission at entrance of Honggu Tunnel in Nanchang[J]. Heating Ventilating & Air Conditioning, 2021, 51(3): 59-65.
[12]	刘靖晗. 城市道路隧道空气污染状况实验研究与数值模拟[D]. 济南: 山东建筑大学, 2020.
[13]	张继业,郑伟范. 交通流随机行为的研究进展[J]. 西南交通大学学报,2016,51(3): 534-545. ZHANG Jiye, ZHENG Weifan. Research on stochastic behavior of traffic flow[J]. Journal of Southwest Jiaotong University, 2016, 51(3): 534-545.
[14]	上海市隧道工程轨道交通设计研究院. 道路隧道设计标准: DG/TJ 08-2033—2017 [S]. 上海: 同济大学出版社, 2017.
[15]	STURM P, BRANDT R, BURKHART J F, et al. Road tunnels: vehicle emissions and air demand for ventilation[M]. Paris: PIARC World Road Association, 2019.
[16]	Technical Committee on Road Tunnels. Road tunnels: emissions, environment, ventilation[M]. Paris: PIARC World Road Association, 1995.
[17]	STURM P, ZUMSTEG F, ALARCON A E, et al. Road tunnels: vehicle emissions and air demand for ventilation[M]. Paris: PIARC World Road Association, 2004.
[18]	STURM P, BRANDT R, BORCHIELLINI R, et al. Road tunnels: vehicle emissions and air demand for ventilation[M]. Paris: PIARC World Road Association, 2012.
[19]	交通部重庆公路科学研究院. 公路隧道通风照明设计规范: JTJ 026.1—1999[S]. 北京: 人民交通出版社, 1999.
[20]	牛国华. 机动车排放因子模型数据库研究[D]. 武汉: 武汉理工大学, 2011.
[21]	刘恒宇,汝宜红. 考虑交通拥堵及工作量平衡性的一致性车辆路径问题[J]. 西南交通大学学报,2016,51(5): 931-937. doi: 10.3969/j.issn.0258-2724.2016.05.016 LIU Hengyu, RU Yihong. Consistent vehicle routing problem considering traffic congestion and workload balance[J]. Journal of Southwest Jiaotong University, 2016, 51(5): 931-937. doi: 10.3969/j.issn.0258-2724.2016.05.016
[22]	NTZIACHRISTOS L, GKATZOFLIAS D, KOURIDIS C, et al. COPERT: a European road transport emission inventory model[C]//Information Technologies in Environmental Engineering: Proceedings of the 4th International ICSC Symposium Thessaloniki. Greece: Springer Berlin Heidelberg, 2009: 491-504.
[23]	中华人民共和国生态环境部. 道路机动车大气污染物排放清单编制技术指南（试行）[R/OL].（2015-01-13）[2023-09-09]. https://www.mee.gov.cn/ywdt/hjnews/201501/ t20150113_294091.shtml.
[24]	徐倩. 基于城市交通特征的中国机动车排放清单建立[D]. 广州: 暨南大学, 2020.
[25]	车轮飞. 公路隧道通风系统设计工程实录[M]. 北京: 中国建筑工业出版社, 2014.
[26]	大气环境司. 中国机动车环境管理年报（2018）[R]. 北京: 生态环境部, 2018.

施引文献

期刊类型引用(2)

1.	陈永，常婷，张冰旺. 混沌映射与中国剩余定理增强的切换认证方案. 西安电子科技大学学报. 2024(04): 192-205 . 百度学术
2.	陈永，刘雯，詹芝贤. 基于混合密钥增强的LTE-R车地认证密钥协商方案. 铁道学报. 2023(06): 69-79 . 百度学术