“碳减排”视域下内河流域梯级枢纽联合通航调度优化

高攀; 方志伟; 赵旭

doi:10.3969/j.issn.0258-2724.20230002

“碳减排”视域下内河流域梯级枢纽联合通航调度优化

doi: 10.3969/j.issn.0258-2724.20230002

三峡大学经济与管理学院，湖北宜昌 443002

基金项目: 国家社会科学基金重大项目（19ZDA089）；湖北省人民政府智力成果采购项目（HBZW-2019-07）

详细信息

作者简介:
高攀（1987—），男，讲师，硕士生导师，研究方向为交通运输管理、物流与供应链管理，E-mail：gaopan1987@ctgu.edu.cn

通讯作者:
赵旭（1982—），男，教授，博士生导师，研究方向为内河航运与船闸调度，E-mail：zhaoxu@ctgu.edu.cn

中图分类号: U641.7；TP18
计量
- 文章访问数: 115
- HTML全文浏览量: 31
- PDF下载量: 29
- 被引次数: 0
出版历程
- 收稿日期: 2023-01-03
- 修回日期: 2023-06-06
- 网络出版日期: 2024-12-18

Optimization of Joint Navigation Scheduling of Cascade Hubs in Inland River Basin from Perspective of Carbon Emission Reduction

College of Economics & Management, China Three Gorges University, Yichang 443002, China

摘要

摘要:
为解决梯级枢纽联合通航调度中船闸运行规则不统一、船舶调度不同步等问题，以船舶综合通航效率、待闸成本和碳排放为决策目标，构建考虑船舶优先级的多维非线性规划（MDNP）模型；然后，拟用改进的回溯多目标模拟退火算法（IBMOSA）对MDNP进行求解，提出梯级枢纽联合调度的优化方案；最后，以“三峡-葛洲坝”梯级枢纽为例，验证MDNP模型与IBMOSA的有效性与可靠性. 结果表明：MDNP模型能够有效兼顾船舶通航效率和公平性，且IBMOSA具有较好的收敛性和全局性；通过编制协同排闸计划，对各船闸闸次安排进行合理布局，可规避梯级船闸的倒闸现象，减少船舶整体过坝时间；与原始调度方案相比，联合调度方案对3个决策目标的优化效率均接近40%，梯级枢纽通航拥堵缓解率约为16%，有效化解了梯级枢纽间的通航矛盾，提升了三峡水域的整体通航效益.
- 交通管理 /
- 联合调度 /
- 梯级枢纽 /
- 模拟退火
Abstract:
In view of non-uniform operation rules of ship locks and non-synchronous ship scheduling in the joint navigation scheduling of cascade hubs, a multi-dimensional nonlinear programming (MDNP) model considering ship priority was constructed with the comprehensive navigation efficiency of ships, the cost of waiting for lock opening, and the carbon emission as the decision-making objectives. Then, the improved backtracking multi-objective simulated annealing algorithm (IBMOSA) was used to solve the MDNP, and the optimization scheme for joint scheduling of cascade hubs was proposed. Finally, the effectiveness and reliability of the MDNP model and IBMOSA were verified by taking the “Three Gorges-Gezhouba” Cascade Hub as an example. The results show that the MDNP model can effectively take into account the navigation efficiency of ships and fairness, and IBMOSA has good convergence and global property. At the same time, through the coordinated drainage sluice plan formulation, the opening of each lock is reasonably arranged, which avoids the phenomenon of switching operation of cascade locks and reduces the overall dam crossing time of ships. Compared with that of the original scheduling scheme, the optimization efficiency of the three decision-making objectives of the joint scheduling scheme is nearly 40%, and the traffic congestion relief rate of cascade hubs is about 16%, which effectively solves the navigation contradiction among cascade hubs and improves the overall navigation efficiency of the Three Gorges water area.
- traffic management /
- joint scheduling /
- cascade hub /
- simulated annealing

HTML全文

有砟道床由大量散体道砟颗粒组成，可保持列车稳定行驶、扩散列车作用在铁轨上的载荷、减小振动吸收轮轨冲击等功能^[1]. 但由于道砟颗粒的尺寸形状具有很强的随机性，并且大多数道砟颗粒存在一定棱角，当其在道床中受到载荷时会产生局部破碎甚至整体破坏，会影响道床整体稳定性使得铁轨发生变形和破损，还会使得道床整体发生不同程度的沉降^[2]. 因此，深入探讨道砟颗粒破碎过程以及破碎机理对有砟道床的劣化评估与养护维修具有重要研究价值.

由于道砟颗粒的几何外形随机性高，导致其破碎机理较为复杂，但仍可以将道砟颗粒破碎分为棱角破碎和整体破坏. 据研究表明，道砟颗粒的破碎与其拉伸强度相关：Mcdowell等^[3]通过对道砟颗粒压碎实验数据的分析发现，其拉伸强度服从威布尔分布. 自离散单元法（DEM）提出以后，针对道砟颗粒破碎的数值模拟逐渐增多：Indraratna等^[4]在二维尺度使用球型颗粒构造几何外形不规则颗粒并进行三轴压缩实验，探讨其在循环载荷下的破碎机理；Lim^[5]在三维尺度使用8个球颗粒构建了形状规则的道砟颗粒并进行道砟箱实验，此方法可以更好地模拟道砟颗粒的几何特征但仍存在很大区别；严颖等^[6]使用离散元软件对道砟颗粒的压破实验进行数值模拟得到了道砟颗粒的尺寸越大其破碎率越大，并对道砟颗粒的离散元接触模型粘结参数进行了标定. 由于缺乏对道砟颗粒几何形态特征的统计分析，在DEM中对刚度、摩擦、阻尼以及颗粒间黏结强度等细观力学参数的标定可能存在偏差^[7]. 为了更准确地研究道砟颗粒的受载破碎行为，改进道砟颗粒形状以及道砟颗粒内部的性能参数成为当前的趋势^[8].

同时，在离散元软件中为了准确地分析道砟颗粒破碎情况以及其内部微裂纹萌生情况，需要选用更为精细的接触模型以及设置与真实道砟颗粒破碎情况相对应的黏结强度. 目前，国内外学者对道砟颗粒离散元接触参数的研究主要基于单个几何外形评价参数，在接触模型选用上多使用Potyondy等^[9]提出的平行粘结模型. 但由于道砟颗粒的几何外形具有随机性以及复杂性，使用单一参数无法准确地描述道砟颗粒的几何外形，所以需要设置多个参数对道砟颗粒的破碎进行定义以此确定其黏结强度.

反向传播（BP）神经网络通过对数据集内部包含的关系进行预测并按照误差反向传播训练，通过不断迭代训练出最优模型，其拥有很强非线性映射能力，很适合求解多元非线性问题^[11]. Ren等^[12]在离散元软件Particle Flow Code（PFC）中使用BP 神经网络对平行粘结模型进行了细观黏结键参数标定，验证了 BP 神经网络在离散元方法中应用的有效性；Wang等^[13]采用神经网络对参数标定以及校核了标定结果流程进行了研究. 但由于 BP 神经网络在训练过程中的拟合程度取决于样本的规律性和数量，规律性不强时预测结果不准确. 通过将遗传算法（GA）与BP神经网络结合，使用GA的迭代搜索寻优功能选取最优BP神经网络参数，增加BP神经网络模型预测精度，可以很好地提高预测精度^[14-15].

为此，本文通过激光光栅扫描的方法得到道砟颗粒的几何参数并对其进行几何外形尺寸定义. 将构造好的道砟颗粒导入离散元软件中，使用刚性块（Rigid-Block）进行填充，并与传统球颗粒（Ball）填充方法对比分析，寻找更能模拟道砟颗粒复杂几何外形以及棱角特性的建模方法，并通过Softbond模型建立黏结键. 依据道砟颗粒单轴压碎实验数据与离散元仿真加载数据对道砟颗粒黏结强度进行标定. 在此基础上，使用GA-BP神经网络模型建立起针对不同几何外形尺寸的道砟颗粒离散元接触参数预测模型，研究分析道砟颗粒几何外形尺寸与其黏结强度的关系，以此揭示道砟颗粒破碎的基本规律.

1. 道砟颗粒单轴破碎实验分析

1.1 实验设置

道砟颗粒的尺寸形状具有随机性，为了对其进行统一分析，首先对道砟颗粒的尺寸形状进行定义. 对于道砟颗粒此类不规则形状多面体，将其放置于一个长方体内，且该长方体的每一个面都与该多面体相切且得到的长方体体积最小时，此长方体可以定义为最小外接矩形.

在单颗粒压碎实验中，首先选取一定数量的不同尺寸与形状的道砟. 将道砟颗粒扫描后获取三维数据，并采用自编算法读取道砟颗粒形状参数. 基于最小包容长方体准则，得到道砟的长轴尺寸（L）、中轴尺寸（W）和短轴尺寸（T），并将长轴、中轴和短轴尺寸求和取均值，定义为等效粒径 $D_{\mathrm{ed}}$ ，如图1.

图 1 道砟颗粒最小外接矩形尺寸

Figure 1. Definition of minimum bounding rectangle dimensions for ballast particles

下载: 全尺寸图片幻灯片

为了进一步研究不同尺寸道砟的力学特性，选取材质为花岗岩的道砟颗粒进行研究. 在选取颗粒过程中，为了更好地研究颗粒几何外形对其破坏特性的影响，选取26组不同尺寸的扁平状、针状和类圆状道砟颗粒并对其进行单轴破碎实验.

采用的实验设备为量程为100 kN的电液伺服万能试验机，如图2所示. 将道砟颗粒的尺寸数据输入至计算机中，并将道砟放置于辅助垫块上，在实验过程中使用道砟颗粒面积最大的面作为底面放置在实验台上的. 将上加载装置下降至接近道砟的位置，此时记录道砟颗粒的摆放形式，以供后续仿真实验还原实验布置. 将加载装置的加载速度设置为0.6 mm/min，数据记录方式设置为力-位移，这样可以准确地记录道砟颗粒所受的加载力，避免了由于道砟颗粒外形复杂导致的共直线问题. 在加载过程中持续观察道砟的破裂情况，以便于停止压缩工作. 当道砟出现明显的贯穿裂纹时停止压缩工作.

图 2 单轴压缩实验装置

Figure 2. Experimental setup for uniaxial compression test

下载: 全尺寸图片幻灯片

1.2 加载力曲线分析

图3所示为道砟颗粒加载力曲线. 由图3可以看到：在压载过程中，力矩随轴向位移的增加存在较大的波动.

图 3 道砟压碎实验加载力曲线

Figure 3. Loading force curve of ballast breakage experiment

下载: 全尺寸图片幻灯片

根据实验过程中观察到的现象，可以将道砟颗粒破碎分为3个阶段，每一个阶段都对应存在一个区域内的峰值：1）在第一个阶段——转动阶段，由于道砟颗粒的形状具有棱角特性，在受到力的压载时，道砟颗粒产生了一定转动，使得道砟颗粒与加载板的距离变大，导致其受到的力减小. 由此，当加载力达到第一个峰值后出现了小范围下降，但由于道砟颗粒的形状不规则，有些较为平整道砟的转动阶段并不明显. 2）随着加载板的轴向位移逐渐增加，道砟颗粒已被压紧，颗粒不再发生转动，此时进入第二阶段——棱角破碎阶段. 由于加载力持续上升达到第二个峰值，并且道砟颗粒边缘棱角部分与加载板接触面积相对较小，使得棱角部分受到的应力超出了其强度，导致棱角破裂，并且道砟颗粒表面由于长期受到外界因素腐蚀导致强度下降，道砟外表面也会产生局部脱落. 此时，道砟颗粒所受加载力再次下降. 3）加载板轴向位移继续增加，道砟颗粒所受力也随之增加直到最高峰值道砟整体破坏，此时对应第三个阶段——整体破坏阶段. 当道砟破坏后，加载力急剧下降，同时压载实验结束.

1.3 等效应力及尺寸影响分析

为进一步验证实验的有效性，对进行加载的26组道砟颗粒载荷数据进行等效应力定义，如式（1），其与等效粒径的关系如图4.

图 4 等效粒径与等效应力关系

Figure 4. Relationship between equivalent stress and equivalent particle size

下载: 全尺寸图片幻灯片

$\sigma_{{\mathrm{ef}}}=\frac{F}{D_{{\mathrm{ed}}}^2},$

(1)

式中: $F$ 为法向最大荷载.

由图4可以看出：道砟颗粒的峰值等效应力随着道砟颗粒等效粒径的增大而减小，与已有研究结果一致^[6]，验证了本实验的有效性；同一等效粒径的道砟颗粒会得出大小不一的等效应力，这是由于道砟颗粒的几何外形具有很高的随机性，在受载过程中，内部应力分布不同，所能承载的载荷也不同.

为分析不同形状颗粒的承载能力，本文对颗粒的尺寸进行进一步定义，根据Zingg^[16]的分类方法，采用伸长率 $= W/L$ 和扁平率 $= T/W$ 定量表征道砟的形态，并得到不同形态参数道砟对应的峰值载荷均值，如图5所示. 由图5可知：道砟颗粒的峰值载荷随着伸长率的增加而增加，尤其是在伸长率为0.9～1.0时，此时的峰值载荷急剧增大，在伸长率为0.9以下时，道砟颗粒的峰值载荷随其尺寸变化较小；当扁平率处于0.8～0.9时，道砟颗粒的峰值载荷最高，在0.5～0.6区间，道砟颗粒抵抗载荷的能力很弱. 由此可见，当道砟颗粒的几何形状近似于类矩形时，其承载载荷的能力较强.

图 5 不同伸长率及扁平率对应的峰值载荷均值

Figure 5. Average peak load corresponding to different elongation and flatness ratios

下载: 全尺寸图片幻灯片

2. 离散元单颗粒破碎分析

2.1 离散元建模

本文对单轴加载试验的道砟颗粒进行激光光栅扫描，构建道砟颗粒真实形状的三维模型. 其中激光光栅扫描建模仪如图6所示，经过扫描后的形状信息将以三角形网格的形式构建出道砟颗粒三维外形.

图 6 激光光栅扫描系统和道砟颗粒构造

Figure 6. Laser grating scanning system and ballast particle structure

下载: 全尺寸图片幻灯片

将扫描后的道砟颗粒文件导入PFC6.0中构建离散元道砟颗粒. 在PFC中，使用球颗粒填充几何模型是近年来最常用的方式，但球颗粒在填充过程中会遇到内应力过大，无法准确模拟道砟颗粒外形等问题. 为解决以上问题，本文将采用Rigid-block进行道砟颗粒填充，并与传统球颗粒填充方式做对比.

如图7所示，对同一道砟颗粒三维模型建立ball填充样本以及Rigid-block填充样本，其中：σ₁为加载力，采用六面体网格法填充球颗粒，Rigid-block填充采用hexahedral方法并设定组成道砟颗粒的刚性块最小边界为9 mm，最大边界为12 mm，刚性块间的最大夹角设定为45°. 在此填充方法下，形状复杂的道砟颗粒也可以很好地被建立，但刚性块之间会存在一些重叠. 为减小刚性块之间的重叠量，通过rblock erode命令改变刚性块的圆度. 侵蚀会产生新的几何形状，通常会产生包含的圆角刚性块，其中面平面保持不变，并且此形状在尺寸上与先前的几何形状基本相同，从而剔除了许多有源触点，这种方法可以提高计算效率. 使用Rigid-block填充的道砟颗粒从填充效果来看，可以更好地还原道砟颗粒的几何形状，尤其是在棱角部分.

图 7 道砟颗粒建模方式以及单轴加载仿真构建

Figure 7. Ballast particle modeling method and uniaxial loading simulation

下载: 全尺寸图片幻灯片

2.2 道砟颗粒离散元接触模型

针对单颗粒破碎仿真，传统平行粘结模型可以很好地反映颗粒间的接触，但当黏结键所受应力超过抗拉强度后，黏结键直接断裂，并没有考虑到岩石存在的软化阶段，无法提供丰富的微观力学过程信息且无法反映岩石的抗压抗拉强度比. 本文将采用Softbond模型对颗粒间以及刚性面间进行粘结，此模型可以反映岩石的软化特性，并且可以更真实地体现道砟颗粒内微裂纹的产生情况.

其中Softbond模型黏结键法向应力 ${\sigma ^{}}$ （ ${\sigma ^{}} \gt 0$ ）和剪切应力 ${\tau ^{}}$ 更新模式如式（2）.

$\left\{\begin{gathered} {\sigma ^{}}{\text{ = }} - \frac{{{F_{\mathrm{n}}}}}{A} + \psi \frac{{\left\| {{M_{\mathrm{b}}}} \right\|R}}{I}, \\ {\tau } = \left\{ {\begin{array}{*{20}{l}} \dfrac{{\left\| {{F_{\mathrm{s}}}} \right\|}}{A} + {0,}\quad{2{\mathrm{D}}}, \\ \dfrac{{\left\| {{F_{\mathrm{s}}}} \right\|}}{A} + {\psi \dfrac{{\left| {{M_{\mathrm{t}}}} \right|R}}{J},}\quad{3{\mathrm{D}}} , \end{array}} \right. \\ \end{gathered}\right.$

(2)

式中： ${F_{\mathrm{n}}}$ 为法向力， ${F_{\mathrm{s}}}$ 为切向力， ${M_{\mathrm{b}}}$ 为扭转力矩， ${M_{\text{t}}}$ 为弯曲力矩， $A$ 、 $I$ 和 $J$ 分别为黏结键横截面的面积、转动惯量和极惯性矩， $\psi$ 为矩贡献因子， $R$ 为黏结键半径.

为SoftBond模型的力学响应规律，图中： ${\sigma _{\mathrm{c}}}$ 为抗拉强度； ${k_{\mathrm{n}}}$ 为黏结键法向刚度； $\zeta$ 为软化系数； ${l_{c}}$ 为最大法向应力时键的伸长量，如式（3）； $l$ 为软化阶段键的伸长量；γ为软化阈值.

图 8 SoftBond模型的力学响应规律

Figure 8. Mechanical response pattern of SoftBond model

下载: 全尺寸图片幻灯片

${l_{\text{c}}} = \frac{{{\sigma _{\text{c}}}}}{{{k_{\text{n}}}}},\quad\zeta \geqslant 0,0 \leqslant \gamma \leqslant 1.$

(3)

当 $\zeta \ne 1.0$ 且 $\gamma \ne 1.0$ 时，黏结键会出现软化行为. 如果黏结键外围 $\sigma \gt {\sigma _{\mathrm{c}}}$ ，则黏结键进入软化状态，软化过程黏结键受到的应力为

${\sigma ^*} = \frac{{{\sigma _{\mathrm{c}}}({l^*} - l)}}{{\zeta {l_{\mathrm{c}}}}},$

(4)

式中： ${l^*}$ 为黏结键最大伸长量.

在选择接触模型后，进行部分宏观接触参数的设置. 对于黏结模型，性能参数与强度参数的设置尤为重要. 首先是性能参数，主要包括杨氏模量 ${E^*}$ ，切向法向刚度比 $k_{\mathrm{r}}$ . 严颖等^[6]通过使用单轴压缩道砟材料样本计算得到了道砟材料的杨氏模量 ${E^*}$ =60 GPa，切向法向刚度比 $k_{\mathrm{r}}$ =0.18，本文以此数据以及花岗岩的宏观特性为参考进行性能参数设置. 强度参数包括软化因子，黏结键破坏阈值、黏结强度. 考虑到道砟颗粒属于脆性岩石，虽然存在软化阶段但处于软化阶段的时间较短，软化因子取值1.0，黏结键破坏阈值取0.8.

针对黏结强度的设置，本文提出峰值等效应力标定法，在仿真过程中计算不同黏结强度下的等效应力破坏峰值，与实际实验等效应力破坏峰值对比，当峰值大小相似时，即确定道砟颗粒的黏结强度. 通过大量仿真发现，不同尺寸的道砟颗粒通常对应不同大小的黏结强度，所以本文将对黏结强度的标定作详细研究.

2.3 道砟颗粒构造方法对比

由于要通过仿真的方式对黏结强度进行大量测试，所以要选择一种更能模拟实际道砟破坏的建模方法，上文分别对球颗粒以及刚性块进行填充构建可破碎道砟颗粒，并如表1所示获取了道砟颗粒的宏观接触参数，本节对2种建模方法的优劣势进行分析，并选择出最合理的一种建模方法.

表 1 道砟颗粒棱角特性仿真接触参数

Table 1. Simulation contact parameters for angular characteristics of ballast particle

微观参数	数值
杨氏模量/GPa	70
剪切刚度比	0.18
软化系数	1.0
软化阈值	0.8
摩擦系数	0.1

下载: 导出CSV

| 显示表格

按照等效应力峰值标定图7中建立的道砟颗粒黏结强度，并对球颗粒以及刚性块建立的道砟颗粒进行单轴压缩仿真，如图9所示.

图 9 不同构造方式的道砟颗粒单轴压缩加载力曲线

Figure 9. Uniaxial compression loading force curves of ballast particles in different construction modes

下载: 全尺寸图片幻灯片

从图9中可以看到：在相同接触参数下，球颗粒的加载力峰值更大且只有一个峰值，其破碎过程曲线更类似于圆柱形岩土样本单轴压缩的破坏曲线，并没有与实际实验的加载力曲线一样存在不同阶段及多个峰值；在使用球颗粒进行压载时，轴向应变在0.02 mm处发生整体破坏，与实际压载实验中的轴向应变不符，但使用刚性块建立的道砟颗粒的轴向应变与实际加载轴向应变非常相近.

由图9（b）可以看到明显的加载力波动：在第一个峰值时，对应道砟颗粒压碎的第一个阶段——转动阶段，道砟颗粒转动后，导致加载板与道砟间存在间隙，所以在第一个峰值后存在一段加载力为0的阶段；当加载板再次接触道砟颗粒后，加载力逐渐上升并且伴随着小幅度波动，此时反映了道砟颗粒压碎实验的第二阶段——棱角破碎阶段；使用刚性块填充的道砟与实际实验的数据趋势相近，并且峰值对应的加载力大小一致.

通过对比得出，使用刚性块方法构建的道砟颗粒其几何外形以及尖角特征与真实道砟颗粒更为相近，在受到加载力作用时，所受应力大小更贴近真实情况. 由此可以看出，使用刚性块填充道砟颗粒，更能真实反映道砟压碎过程，并且获得的数据更为准确.

2.4 道砟颗粒压缩破碎分析

上文已验证使用Rigid-block建立道砟颗粒更具有合理性，所以将26组进行单颗粒压碎实验的道砟颗粒在PFC6.0中使用刚性块建模并逐个进行单轴破碎仿真实验，依据等效应力峰值标定出每一个颗粒所对应的黏结强度. 并对其中一组道砟颗粒压碎过程的细观力学特性进行分析，如图10.

图 10 道砟颗粒受压过程中的破碎形式

Figure 10. Breakage patterns of ballast particles during compression process

下载: 全尺寸图片幻灯片

如图10（a）所示的道砟颗粒外轮廓类似于矩形，此类形状的道砟颗粒在压碎过程中不易产生转动. 由图10可以看出：在压碎过程中，由于应力集中，道砟颗粒首先产生棱角破碎，并且棱角破碎的位置与大小基本与实际实验中的情况一致；随着轴向位移增加，道砟颗粒产生一条明显贯穿道砟颗粒的宏观裂纹；当轴向位移进一步增加，道砟颗粒发生整体破坏.

图11为道砟颗粒所受加载力与裂纹萌生的关系曲线. 可以看出：道砟颗粒在轴向应变为0～1 mm发生小幅度转动，转动阶段并未产生微观裂纹，当轴向位移接近1 mm时，加载力曲线产生了一个小峰值，此时发生了道砟颗粒的棱角破碎；随着轴向位移增加到1～2 mm，加载力产生了多个峰值，此时微观裂纹开始萌生并且不断积累，道砟颗粒内部产生了一定的宏观裂纹，如图10（b）所示在道砟颗粒外形上，还未能看出明显地岩体破坏；当轴向位移达到2～3 mm，加载力峰值较小，但仍产生了大量的微观裂纹，这是由于道砟颗粒内部的宏观裂纹扩展延伸至道砟颗粒整体使其破碎为多个碎块，此时，道砟颗粒在外观上可以看到明显的贯穿裂纹，可以将道砟颗粒定义为整体破坏；当进一步增加轴向位移后，裂纹的产生数量变少，道砟颗粒破碎为多个碎块.

图 11 道砟颗粒破碎中所受加载力和微裂纹数量随轴向应变变化曲线

Figure 11. Variation of loading force and number of microcracks with axial strain during ballast particle breakage

下载: 全尺寸图片幻灯片

通过总结分析所有道砟颗粒的破碎过程以及破碎情况：在加载力达到峰值以前，产生的裂纹主要分布在道砟顶部或底部的接触载荷点附近，反映了不规则道砟颗粒的棱角特性；随着裂纹逐渐积累，裂纹生长速度提高并产生裂缝，从而引起道砟局部结构断裂；道砟的起裂点一般出现在棱角特性最明显的位置，而裂纹则由该点向道砟另一端的接触区进行扩张；当道砟表面越平整时，与墙体间的接触分布越均匀，轴向力引起的弯矩也越小，这种情况下，道砟被压碎的可能性很小，故而在室内实验中观察到的峰值强度很高. 综上所述，道砟的局部棱角结构和形状分别决定其起裂点位置和裂纹方向，并最终呈现出不同的破碎模式和峰值响应.

3. 机器学习参数标定

3.1 机器学习模型选用

本文在进行单颗粒破碎仿真时发现，不同尺寸形状的道砟颗粒需要设置不同的黏结强度参数. 通过对所有道砟颗粒进行离散元仿真，得出了每种道砟颗粒需要对应的黏结强度. 依据26组实验数据以及对应的仿真数据建立数据集，在分析数据集过程发现，道砟颗粒的黏结强度随尺寸的增大而增大，但相近尺寸道砟颗粒的黏结强度也会有较大区别，是由于道砟颗粒的粒径以及几何外形尺寸存在很大的随机性^[17]. 通过传统线性拟合等方式无法准确地获取对应的黏结强度参数，若逐一对不同道砟颗粒进行黏结强度参数标定，则效率较低. 为了更好地得到道砟颗粒黏结强度，本文基于获取的黏结强度数据集进行神经网络预测. 采用GA-BP神经网络的方式对道砟颗粒离散元黏结键接触参数进行标定.

使用神经网络模型不但可以从几何外形评价指标（道砟颗粒的等效粒径、最小外接矩形的长轴、中轴和短轴）对黏结强度进行分析，而且可以预测所需尺寸道砟颗粒的粘结参数，解决了实验数据组不足的问题.

由于道砟颗粒的几何尺寸参数对离散元黏结强度影响较大，并且等效粒径相似的道砟颗粒会因为其最小外接矩形的3个轴尺寸不同，导致最终标定的黏结强度不同. 为了更好地寻找到神经网络的参数阈值，提高训练准确度和速度，可通过遗传算法（GA）建立神经网络模型^[18-19].

3.2 GA-BP神经网络构建

本文26组仿真数据中17组数据拟合神经网络，9组数据作为验证集，该验证集专门用来对比BP 与 GA-BP 神经网络的预测误差，并采用随机种子进行数据分组. BP 神经网络相关参数设置如下：训练函数采用Levenberg-Marquardt，训练次数设置为1000，学习速率设置为0.01，最小误差设置为0.00001，归一化函数采用Mapminmax，如图12所示，通过自编算法选取隐含层节点数量. 通过计算每一隐含层数量对应的预测结果与原数据的均方误差（MSE），选取出MSE最小的隐含层数量为9层，输入层数量为4层，输出层数量为1层.

图 12 隐含层节点数对应训练集MSE变化

Figure 12. MSE change with number of hidden layer nodes in training set

下载: 全尺寸图片幻灯片

由于 BP 神经网络神经元间的参数阈值是随机选取的，神经网络对变化幅度较大的数据集拟合能力不强，利用GA的种群迭代对BP神经网络的参数阈值进行选优可以有效提高拟合精度. 本研究中GA的参数设置如下：初始种群规模为 60，最大进化代数为 210，交叉概率为 0.8，变异概率为 0.2，自变量上下限分别设置为3和−3.

图13为BP神经网络和GA-BP神经网络模型预测值与真实值的对比. 由图13可以明显看出：经过GA优化后的BP神经网络训练结果更准确，预测值与原数据拟合度达到了0.98160，而未经过优化的BP神经网络模型拟合度仅有0.80828. 图14为GA-BP与BP神经网络的预测误差对比，可以看出：在样本编号为6和9处，未优化的BP神经网络误差明显，而GA-BP神经网络的预测误差较小. 由此证明，经过此参数训练的GA-BP神经网络模型可以很好地预测道砟颗粒的黏结强度.

图 13 BP神经网络和GA-BP神经网络模型预测值与真实值的对比

Figure 13. Comparison of predicted values and actual values of BP and GA-BP neural network models

下载: 全尺寸图片幻灯片

3.3 道砟颗粒黏结强度预测

使用训练好的GA-BP神经网络模型对100个不同尺寸形状的道砟颗粒进行黏结强度预测. 其中道砟颗粒的选取按照铁道行业道砟级配铺设标准进行，过筛尺寸分别为16.0～25.0、25.0～35.5、35.5～45.0、45.0～56.0 mm和56.0～63.0 mm，每一级配范围对应20组道砟颗粒. 如图15所示，首先将100个道砟颗粒进行激光扫描，获取道砟颗粒最小外接矩形的长轴、中轴和短轴并计算其等效粒径. 通过分析100个道砟颗粒所属的级配尺寸范围与等效粒径尺寸，并使用等效粒径概念重新将上述5组道砟颗粒过筛尺寸进行分配为25～39、39～48、48～56、56～64 mm和64～80 mm. 使用长轴、中轴、短轴和等效粒径4个参数进行神经网络预测.

图 14 GA-BP与BP神经网络的预测误差对比

Figure 14. Comparison of prediction errors between GA-BP and BP neural networks

下载: 全尺寸图片幻灯片

图 15 5组不同等效粒径的道砟颗粒

Figure 15. Five groups of ballast particles with different equivalent particle sizes

下载: 全尺寸图片幻灯片

首先，将道砟颗粒等效粒径、长轴、中轴、短轴4个参数的数据进行归一化处理；将归一化后的数据导入训练好的GA-BP神经网络模型中进行黏结强度预测，100个道砟颗粒黏结强度预测结果如图16. 从图中可以看出：黏结强度随着道砟颗粒等效粒径的增加而增加，并且等效粒径处于[25，39）、[39，48）、[48，56） mm 3个区间的道砟颗粒黏结强度为120～180 MPa，56～64 mm和64～80 mm尺寸的道砟颗粒的黏结强度密集分布在170～200 MPa，因为道砟颗粒的尺寸越大，其黏结强度受到几何外形差异的影响越小. 由此可以看出，道砟颗粒的尺寸对其离散元接触参数有较大影响.

图 16 100个道砟颗粒的GA-BP神经网络模型黏结强度预测结果

Figure 16. GA-BP neural network model prediction results for bond strength of 100 ballast particles

下载: 全尺寸图片幻灯片

由于实际铁路中使用的道砟颗粒等效粒径大小不一，在离散元仿真中研究人员通常会选取几十种典型道砟颗粒进行道床建模. 为更好地模拟典型道砟颗粒的破碎情况，对于黏结强度参数的设置本文提供2种方案：1）提取道砟颗粒的长轴、中轴和短轴尺寸，并计算出与其对应的等效粒径，将以上4个参数导入训练好的GA-BP神经网络模型中进行黏结强度预测，此方法可以较为精确地获取每个道砟颗粒的黏结强度；2）当所需要标定的道砟颗粒数量较多，并且道砟颗粒的尺寸差异较大时，可以计算出道砟颗粒的等效粒径并依据表2获取黏结强度. 在离散元仿真中，每一组道砟颗粒设置一个对应的黏结强度，此方法可以提高配置模型参数的效率同时也可以获取较为精确的黏结强度参数.

表 2 依据等效粒径分组的道砟颗粒黏结强度

Table 2. Bond strength of ballast particles grouped according to equivalent particle size

等效粒径范围/mm	黏结强度范围/MPa	黏结强度均值/MPa
[25，39）	130.37～155.34	151.85
[39，48）	132.22～186.58	159.45
[48，56）	138.58～195.69	166.71
[56，64）	164.99～202.46	175.29
[64，80）	177.84～201.76	185.29

下载: 导出CSV

| 显示表格

为进一步验证神经网络标定结果，将等效粒径为54.15 mm，长、中、短轴分别为71.58、58.2、32.28 mm的道砟颗粒进行神经网络黏结强度预测，得到其预测黏结强度为165.18 MPa. 选用与表1相同的接触参数，并设定相应黏结强度，使用Rigid-Block建模进行单颗粒破碎仿真并与实验数据相对比，如图17所示. 由图17可知：神经网络预测结果所标定的黏结强度可以很好地反映道砟颗粒在单颗粒破碎实验中的峰值加载力，破碎裂纹也与实际岩石相似. 由此可见，使用GA-BP神经网络可以高效准确地预测道砟颗粒的黏结强度.

图 17 单颗粒破碎仿真与实际实验加载力对比

Figure 17. Comparison of single particle breakage simulation and actual experimental loading force

下载: 全尺寸图片幻灯片

有砟道床由多种材质的道砟石组成，不同材质岩石的宏观力学特性不同，导致其在颗粒尺度上，组成岩石样本的颗粒黏结键承受的应力水平不一致. 针对不同种类岩石，需要配置不同的接触参数重新进行标定，并构建相对应岩石的数据集以建立神经网络模型进行黏结强度预测.

4. 结　论

1）通过对单个道砟颗进行单轴压碎实验，确定了其破坏力以及等效应力，结合其破碎过程将道砟颗粒破碎分为转动阶段、棱角破碎阶段和整体破坏阶段. 此外，使用最小外接矩形法对道砟颗粒的长轴、中轴、短轴以及等效力矩进行规定.

2）使用PFC6.0对扫描道砟颗粒的构建方式进行对比分析，使用刚性块构建的道砟颗粒可以更精确地反映道砟颗粒在破碎过程中的特性. 对使用刚性块填充的道砟颗粒进行单轴压碎仿真得到，道砟的局部棱角结构和形状决定其起裂点位置和裂纹方向，以及不同几何外形道砟颗粒的破碎模式和峰值响应.

3）利用 GA-BP 神经网络构建其预测数学模型. 黏结强度随着道砟颗粒等效粒径的增加而增加，提供神经网络标定以及等效粒径分组标定2种方案对黏结强度参数进行设置，同时使用预测结果进行离散元仿真，进一步验证了神经网络模型的正确性，为道砟颗粒破碎的离散元数值分析提供了一定理论参考.

图 1 “五闸一机”联合通航系统运作流程

Figure 1. Operation flow of joint navigation system with five locks and one ship lift

下载: 全尺寸图片幻灯片

图 2 IBMOSA的流程框架

Figure 2. Framework of IBMOSA

下载: 全尺寸图片幻灯片

图 3 IBMOSA与NSGA的收敛速度

Figure 3. Convergence speeds of IBMOSA and NSGA

下载: 全尺寸图片幻灯片

图 4 IBMOSA与NSGA的开闸时间计划波动值对比

Figure 4. Comparison of schedule fluctuation values for lock opening of IBMOSA and NSGA

下载: 全尺寸图片幻灯片

图 5 优化模型的灵敏度与稳定性验证

Figure 5. Sensitivity and stability verification of optimization model

下载: 全尺寸图片幻灯片

图 6 各船闸的闸次时间安排

Figure 6. Schedule of ship lock opening

下载: 全尺寸图片幻灯片

图 7 船舶待闸时间比较

Figure 7. Comparison of time costs of ships waiting for lock opening

下载: 全尺寸图片幻灯片

表 1 变量说明

Table 1. Variable specification

变量类型	符号	说明
模型参数	$U_{{\mathrm{weight}}}$	静态权重与动态权重的相关系数
	$\alpha$	静态权重在综合权重中的比例系数
	${L_k}$	船舶 $k$ 的长度， ${\text{m}}$
	${P_1}$	燃油单位价格，元/L
	${C_1}$	船舶单位时间供电燃油消耗，L/min
	${A_i}$	闸室 $i$ 的闸室面积，m²
	${X_k}$	船舶 $k$ 的横坐标
	${F_i}$	由于倒闸（葛洲坝和升船机）或换向（三峡南线和北线）引起闸室 $i$ 的附加运行时间， ${\text{s}}$
	${a_k}$	船舶 $k$ 的船舶面积，m²
	${Q_i}$	船闸 $i$ 在1个调度期中开闸总次数
	${\omega _i}$	葛洲坝船闸 $i$ 的开闸次数在 3 个葛洲坝船闸总开闸次数中的分配比例
	${\gamma _{m}}$	船舶信息m的权重系数
	${K_1}$	葛洲坝运行状态相关的参数
	${X_i}$	闸室 $i$ 的长度， ${\text{m}}$
	${n_i}$	闸室 $i$ 的开闸次数，次
	${T_{i jk}}$	船舶 $k$ 在闸室 $i$ 第 $j$ 闸次的到闸时刻
	$\beta$	动态权重在综合权重中的比例系数
	${W_k}$	船舶 $k$ 的宽度， ${\text{m}}$
	${P_2}$	人工成本，元/min
	$\mathrm{E_{CO}}_{_2}$	船舶 ${{\mathrm{CO}}}_2$ 排放系数，kg/L
	${B_i}$	升船机的闸室面积，m²
	${Y_k}$	船舶 $k$ 的纵坐标
	${E_i}$	闸室 $i$ 运行一次的时间， ${\text{s}}$
	$N$	申请过闸的船舶总数，艘
	${P_{i jk}}$	船舶 $k$ 通过闸室i 第 $j$ 次升降时的重量， ${\text{t}}$
	${M_i}$	闸室i计划期内通过的船舶面积和，m²
	$\theta$	影响船舶信息优先级的量化值
	${K_2}$	葛洲坝运行状态相关的参数
	${Y_i}$	闸室 $i$ 的宽度， ${\text{m}}$
决策变量	${{{W}}_{\mathrm{S}}}$	静态权重值
	${S_{i jk}}$	船舶 $k$ 是否通过船闸（或升船机）的0-1 决策变量， ${S_{i jk}} = 1$ 表示船舶 $k$ 通过闸室 $i$ 第 $j$ 闸次，反之 ${S_{i jk}} = 0$
	${d_{i j}}$	闸室 $i$ 第 $j$ 闸次的运行方向，1 为上行，0 为下行
	${t_{i jk}}$	船舶 $k$ 在闸室 $i$ 第 $j$ 闸次的过闸时刻
	${W_{\mathrm{D}}}$	动态权重值
	$W$	综合优先级权重

下载: 导出CSV

表 2 通航拥堵优化效率

Table 2. Optimization efficiency of navigation congestion

目标函数值	联合调度方案	原始调度方案	优化效率
总目标值	76035.4	193925.0	0.392086631
${F_1}$ 平均值	74452.1	189937.3	0.391982512
${F_2}$ 平均值	1583.3	3987.7	0.397045916
${F_3}$ 平均值	527.8	1329.2	0.397080951
综合通航效率值	158830.0	189940.0	0.166417290

下载: 导出CSV

参考文献(15)

[1]	交通运输部. 交通运输部关于印发《内河航运发展纲要》的通知[EB/OL]. (2020-05-29)[2023-04-17]. https://www.gov.cn/zhengce/zhengceku/2020-06/04/content_5517185.htm.
[2]	赵旭,尹熙琛,高攀. 主导权异质性视角下的过坝方式货运量分担率研究[J]. 运筹与管理,2020,29(8): 192-201. ZHAO Xu, YIN Xichen, GAO Pan. Research on the share of freight volume of via-dam method under the perspective of dominance heterogeneity[J]. Operations Research and Management Science, 2020, 29(8): 192-201.
[3]	长江三峡通航管理局. 2022年三峡枢纽航运通过量近1.6亿吨[EB/OL]. (2023-01-01)[2023-04-17]. https://cjhy.mot.gov.cn/xw/slxw/202301/t20230104_287113.shtml.
[4]	长江三峡通航管理局. 宜昌“充电站”为绿色长江赋能 [EB/OL]. (2023-01-07)[2023-04-17]. https://news.hubeidaily.net/mobile/1040558.
[5]	JI B, YUAN X H, YUAN Y B, et al. An adaptive large neighborhood search for solving generalized lock scheduling problem: comparative study with exact methods[J]. IEEE Transactions on Intelligent Transportation Systems, 2020, 21(8): 3344-3356. doi: 10.1109/TITS.2019.2926405
[6]	张红颖,周子林,李彪. 基于多Agent的通航运力资源协同调度[J]. 交通运输系统工程与信息,2020,20(1): 214-221. ZHANG Hongying, ZHOU Zilin, LI Biao. Collaborative schedule of general aviation resource based on multi-agent[J]. Journal of Transportation Systems Engineering and Information Technology, 2020, 20(1): 214-221.
[7]	马亮,胡宸瀚,金福才,等. 铁水运输调度双层多目标约束优化模型研究[J]. 西南交通大学学报,2023,58(2): 357-366,397 MA Liang, HU Chenhan, JIN Fucai, et al. Research on double-layer multi-objective constrained optimization model for hot metal transportation scheduling[J]. Journal of Southwest Jiaotong University, 2023, 58(2): 357-366,397
[8]	高攀,刘顺,赵旭,等. 绿色通航视域下过坝船舶预约调度双目标优化研究[J]. 交通运输系统工程与信息,2022,22(5): 293-299,336. GAO Pan, LIU Shun, ZHAO Xu, et al. Bi-objective optimization of ship dam-passing appointment scheduling considering green navigation[J]. Journal of Transportation Systems Engineering and Information Technology, 2022, 22(5): 293-299,336.
[9]	GOLAK J A P, DEFRYN C, GRIGORIEV A. Optimizing fuel consumption on inland waterway networks: local search heuristic for lock scheduling[J]. Omega, 2022, 109: 1-12. doi: 10.1016/j.omega.2021.102580
[10]	JI B, ZHANG D Z, YU S S, et al. Optimally solving the generalized serial-lock scheduling problem from a graph-theory-based multi-commodity network perspective[J]. European Journal of Operational Research, 2021, 288(1): 47-62. doi: 10.1016/j.ejor.2020.05.035
[11]	JI B, ZHANG D Z, YU S S, et al. An exact approach to the generalized serial-lock scheduling problem from a flexible job-shop scheduling perspective[J]. Computers & Operations Research, 2021, 127: 1-17.
[12]	JI B, ZHANG D Z, YU S S, et al. Mathematical programming models for scheduling multiple cascaded waterway locks[J]. Computers & Industrial Engineering, 2021, 156: 1-13.
[13]	章建玲,王澎涛,袁晓辉,等. 基于进化算法和气象信息的三峡-葛洲坝通航调度研究[J]. 水电能源科学,2017,35(2): 91-94. ZHANG Jianling, WANG Pengtao, YUAN Xiaohui, et al. Co-navigation scheduling of Three Gorges-gezhouba Dam based on evolutionary algorithm and meteorological information[J]. Water Resources and Power, 2017, 35(2): 91-94.
[14]	鲁彦汝. 基于随机Petri网的三峡-葛洲坝通航系统联合调度研究[D]. 武汉:武汉理工大学,2017.
[15]	甘家华,毛新华,赵京. 考虑动态交通流的停车场车辆最优疏散模型[J]. 西南交通大学学报,2021,56(1): 123-130. GAN Jiahua, MAO Xinhua, ZHAO Jing. Optimal evacuation model of parking vehicles in dynamic traffic flows[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 123-130.