面向网络分配的货运方式选择与OD矩阵标定框架——基于欧洲多式联运模型的实证研究

《European Transport Research Review》：From modal-choice estimation to freight OD-matrix assignment: a calibrated framework for European multimodal transport models

　　

货运方式选择模型通常以其再现总体方式分担率的能力来评判。然而，对于网络分配而言，这并不足够：模型还必须再现吨位在方式起讫点（OD）矩阵中的分布位置。本文提出了一种将方式选择估计与货运矩阵分配相联系的标定框架。基于公路、内河航运（IWW）和铁路的汇总OD-商品流数据，估计了一个行为多项Logit（MNL）模型。该模型使用了相对于公路表达的变量、商品组常数、方式特定距离样条函数以及铁路市场结构修正项。由于稀疏的铁路和内河航运矩阵无法仅靠行为模型精确再现，因此在非公路效用中添加了有界基年校准层。该框架在三个欧洲数据集上进行了测试：欧洲NUTS2矩阵、“比荷卢+”NUTS3矩阵以及德国国内NUTS3矩阵。评估在总体方式分担率层面和OD单元层面进行，采用加权绝对百分比误差（WAPE）。详细的标定结果被解释为样本内基年重建，而转移性则通过重复的OD对留出实验单独检验，用以区分基年重建与未见OD对上的表现。结果表明，行为层提供了有用的结构，但基年方式OD矩阵的再现需要详细的有界标定。随后，所估计的模型被集成到专为货运设计的交通网络建模软件Nodus中，使得标定后的框架可用于网络分配。

论文解读：面向网络分配的货运方式选择与OD矩阵标定框架

一、研究背景与问题提出

在欧洲，货运政策日益依赖于能够比较不同定价、基础设施或网络性能假设下公路、铁路和内河航运方案的综合模型。然而，现有货运模型面临一个核心矛盾：一个方式选择模型可能很好地再现了总体的铁路或内河航运分担率，却将这些吨位分配到了错误的起讫点（OD）单元上。这样的模型完全不适合用于网络分配。这一问题对于非公路方式尤为尖锐，因为铁路和内河航运矩阵极其稀疏：许多OD-商品单元在技术上可行，但并未承载任何观测到的非公路流量，而少数走廊和商品却占据了绝大部分交通量。

因此，研究的核心问题在于：如何估计、标定和评估一个汇总货运方式选择模型，使其既能保持行为可解释性，又能精确再现方式OD矩阵，以满足网络分配的需求？传统的离散选择模型（如多项式Logit模型，MNL）虽然能提供行为结构，但在处理稀疏的非公路矩阵时，会为所有技术可行的模式赋予正概率，从而将少量非公路交通扩散到大量实际未使用的OD单元上。同时，成本和时间的强相关性也导致系数估计不稳定。为此，研究人员提出了一个结合行为模型与有界基年校准层的双层框架。

二、主要关键技术方法

1. 1.
   行为多项式Logit模型（MNL）：基于汇总OD-商品流数据，构建了一个以公路为参照的行为模型。该模型创新性地采用了相对成本和时间比率（取对数）、商品组常数、方式特定距离样条函数（阈值设为150km和300km）以及铁路市场结构修正项（以内河航运可用性作为代理变量）。
2. 2.
   有界基年校准层：为解决行为模型无法精确再现稀疏非公路矩阵的问题，在非公路方式的效用函数中添加了有界校准常数。这些常数通过迭代比例校正算法在效用空间中更新，并被限制在[-8, +8]区间内，以防止过拟合。
3. 3.
   模型评估指标：采用加权绝对百分比误差（WAPE）和偏差（Bias）作为OD单元层面的主要评估指标，克服了平均绝对百分比误差（MAPE）在稀疏矩阵中不稳定的缺陷。
4. 4.
   转移性诊断：通过重复的OD对留出实验（Holdout Experiment），将20%的OD对完全排除在训练集外，以评估模型对未见OD对的预测能力，区分基年重建与真正的预测性能。
5. 5.
   模型实现：所估计的模型被集成到专为货运设计的交通网络建模软件Nodus中，实现了标定后框架的网络分配功能。

三、研究结果

6.1 行为层结果

• •
  内河航运（IWW）成本效应最稳健：在所有三个数据集中，IWW成本比率的系数均为负值，表明当IWW成本相对于公路上升时，其效用下降。例如，欧洲数据集中，相对IWW成本每增加10%，IWW的几率比约降低30%。
• •
  移动时间系数不稳定：在欧洲NUTS2数据集中，IWW和铁路的移动时间系数为负，符合预期。但在两个NUTS3数据集中，这些系数达到零边界。这并非表明运输时间不重要，而是因为在更精细的空间分辨率下，移动时间与成本、距离和地理特征高度相关，无法提供独立的额外信息。
• •
  距离效应因数据集而异：距离样条系数显示出强烈的数据集特异性。在欧洲层面，短距离对IWW和铁路不利；而在比荷卢+和德国，长距离段系数更大，表明非公路方式在更长距离上更具竞争力。
• •
  铁路市场结构修正项有效：IWW可用性变量在铁路效用中的系数为正，表明位于河运走廊的OD对，铁路也更有可能被使用，起到了市场结构指示器的作用。

6.2 距离样条阈值的敏感性分析

• •
  将距离样条阈值从基准的150/300公里分别调整为100/250公里和200/400公里后，预校准WAPE的变化非常小（最大差异约为0.031）。这表明模型结果对合理的阈值选择不敏感，基准阈值被保留作为可解释的共同分段。

6.3 铁路成本项的敏感性分析

• •
  在行为模型中引入铁路成本项后，只有欧洲数据集表现出轻微的改善（预校准铁路WAPE从1.088降至1.073），而在两个NUTS3数据集中，铁路成本系数达到零边界。该变量未能提供跨数据集的稳定成本响应，且对基年校准后的矩阵再现几乎没有影响。因此，更简洁的不含铁路成本的模型被选为参考模型。

6.4 基年校准层结果

• •
  应用详细的“方式-OD-商品组”校准后，所有数据集和方式的WAPE值均大幅下降，表明该层有效地恢复了基年矩阵的地理分布和稀疏性。例如，欧洲数据集的内河航运WAPE从1.428降至0.087。这是样本内基年重建的成功，而非预测能力的证明。

6.5 校准常数分布分析

• •
  未加权的校准常数分布显示，大量常数处于负边界（-8），尤其对于铁路（80%-89%）。这对应了大量技术上可行但无观测流量的OD单元，校准层通过施加极负的效用修正来抑制这些单元的虚假正概率。
• •
  按观测吨位加权后，几乎没有任何重要流量与负边界的常数相关联。相反，承载非公路流量的单元，其中位数校准常数为正值（约2.4至3.2），表明行为模型在这些关键走廊上低估了非公路方式的吸引力。

6.6 OD对留出实验与转移性诊断

• •
  在留出实验中，使用较粗略的“方式-商品组-距离”校准层。结果显示，所有27次比较（3个数据集×3种方式×3次随机分割）中，校准后的WAPE均低于预校准WAPE，表明模型结构具有一定的转移能力。
• •
  然而，改善幅度有限。例如，欧洲数据集的IWW WAPE从1.451降至1.382。这表明粗略校准层能捕捉部分可转移的结构（如商品组和距离效应），但无法像详细的OD级校准那样精确再现基年地理特征。该实验应被视为转移性的诊断，而非对详细校准的替代。

四、总结与讨论

本研究提出了一个连接货运方式选择估计与矩阵分配的双层标定框架。结论如下：

1. 1.
   行为层提供必要结构：提出的相对公路MNL规范，特别是稳健的IWW成本效应和铁路市场结构修正项，为理解系统性的方式选择行为提供了有用基础。
2. 2.
   有界校准层是实现矩阵再现的关键：仅靠行为模型无法精确再现稀疏的非公路OD矩阵。详细的、有界的“方式-OD-商品组”校准是生成可用于网络分配的基年矩阵的必要环节。
3. 3.
   区分基年重建与预测能力：评估策略明确区分了样本内基年重建（通过详细校准实现低WAPE）与对未见数据的转移性（通过留出实验评估）。留出实验显示，模型具有一定但有限的转移能力。
4. 4.
   框架的应用模式：对于基年分配，应采用详细的“方式-OD-商品组”校准；对于政策情景分析，则应采用枢轴点（pivot-point）方法，即固定基年校准常数，仅更新行为属性，以此模拟围绕基年的变化。

研究结论翻译：本研究证明了将行为方式选择模型与有界基年校准相结合的有效性，为欧洲多式联运货运网络分配提供了一个实用且透明的框架。该框架成功解决了传统模型在稀疏矩阵分配中的不足，并为未来的政策敏感性分析和模型转移性研究奠定了基础。