一、反推场景
OD 矩阵是交通规划、路网设计、项目评估的核心基础数据，传统大规模 OD 调查成本高、周期长、更新难。OD 反推（ODME）是利用路段实测流量，结合路网与初始 OD，通过 “分配 — 校核 — 迭代” 反推全路网 OD 矩阵，属于交通分配的逆过程。
- 高速公路 / 干线公路网：缺乏全域 OD，仅有收费站、断面流量或 ETC 数据；
- 城市宏观模型更新：旧 OD 过时、新调查不足，需用现状路段流量校准；
- 项目可研 / 工可：无历史 OD，仅路网 + 观测流量，快速生成基线 OD；
- 路网改扩建 / 瓶颈分析：反推现状 OD，用于 “有无项目” 对比及流量再分配。
二、数据准备
1、路网（线层）必备字段（双向）
AB-V/BA-V：观测流量（pcu/d，与 OD 口径一致）AB-C/BA-C：基准通行能力（pcu/d）
- 比如：120km/h 高速、单向 2 车道 → 单车道 2200 pcu/h → 2×2200/0.125=35200 pcu/d
AB-T/BA-T：路段行驶时间（min），由长度 / 速度计算AB-S/BA-S：设计速度（km/h）- 拓扑：节点连通、无悬挂路段、方向一致。
注：若采用自定义阻抗函数，可以在路段属性表里
- 新建一个字段，比如
AB_FY、BA_FY，用公式算出每条路段的初始广义费用（自由流状态下）。 - 把
AB_FY、BA_FY当成 “伪时间”，填入Time字段，然后继续用 BPR 函数，把它当作延误函数的基础值。 - 本质上是用 BPR 来描述 “广义费用随流量增加而上升” 的关系。
2、交通小区（面层 / 点层）
- 按行政区 / 功能区划分，质心与路网节点连接（Centroid Connector）；
- 小区编号连续，与 OD 矩阵行列号一一对应。
注：高速公路OD反推，可以将区域内每个收费站作为单独小区。
3、初始 OD 矩阵（种子矩阵）
- 行列数 = 小区数，对角线 = 0（无区内出行）；
- 三种常用来源：
- 1）历史 OD 或小样本调查扩算；
- 2）全 1 矩阵（无先验信息时）；
- 3）重力模型生成（基于小区阻抗）
- 如果没有历史OD，但有各收费站的进出口量，可以利用重力模型（后文有对重力模型的介绍），用P/A约束获得OD分布矩阵作为初始OD，将该矩阵作为反推OD的输入。
- 该操作能够在保证PA总量的同时，使反推OD精准贴合实际流量。
三、反推过程（P、A约束的反推）
- 出行生成：计算各小区P_i、A_j；
- 阻抗矩阵：用路网生成小区间时间阻抗；
- 重力模型：双约束生成Gravity_OD.mtx（带 P/A 结构）；
- ODME 反推：以 Gravity_OD 为种子，用路段观测流量反推，得到ODME_Final.mtx；
- 矩阵平衡：强制行和 = P、列和 = A，得到ODME_Balanced.mtx；
- 分配校验：误差达标 → 最终 OD 可用。
四、实例分析
1、首先，获得路网自由流阻抗：
2、其次选择应用重力模型，用P、A约束初始OD矩阵
3、获得初始OD之后，对初始OD进行索引转换后，作为OD反推的输入
4、最后进行OD反推
5、精度检验
对比反推结果和现状数据，大概10%以内的误差（这个误差可接受，因为现状数据之间可能存在匹配性误差，比如路段和收费站数据统计口径的问题）
如果存在较大误差，需要分析流量流向，手动调整OD单元→再分配→再修正，直至误差控制在10%以内。
五、问题解答
对以上提到的重力模型进行补充介绍，在交通分布阶段，此方法的应用尤为重要。
1、重力模型介绍
在介绍初始OD矩阵之间，先讲一下公路/交通规划主流的重力模型（双约束逆幂重力模型）,对应TransCAD中的inverse power重力模型。
模型主要考虑两个因素：距离、区域空间位置的影响
1）距离：Fij——参数b（阻抗敏感系数）
① 物理含义
表征出行量对出行阻抗（时间、费用、距离）的敏感程度： 阻抗越大（路越远、耗时越久、费用越高），出行意愿下降越快，b 就是控制衰减快慢的核心参数。
②取值规律与直观理解
b 越大出行对阻抗越敏感：
- 距离稍远、时间略增，出行量就大幅衰减；
- 短距离出行占比极高，长距离跨区出行很少。 适用：城市内部短途出行、普通国省道、集散道路。
b 越小出行对阻抗越不敏感：
- 即便距离远、耗时久，依然有大量跨区出行；
- 长距离、跨区域出行活跃。 适用：高速公路、干线公路、城际 / 省际长途出行（你当前项目）。
③工程常用取值（公路体系）
- 城市内部出行：b取1.8~2.5
- 普通国省干线：b取1.2~1.8
- 高速公路 / 城际长途：b取0.8~1.5（你高速项目优先取 1.0～1.5）
④和阻抗的联动
用 ** 广义费用（元）** 代替出行时间做阻抗，规则不变： b 越大，出行对 “通行费 + 时间 + 里程” 综合成本越敏感，大家更倾向选择低成本近路。
2）摩擦因子K（区域调整系数/平衡因子）
① 物理含义
②典型使用场景举例
③TransCAD 里的表现
④区分：Fij 和 Kij
很多人混淆两个 “因子”，一句话分清：
2、为什么进行OD反推，需要加P、A约束
3、如果没有现状P、A数据，应该怎么获得
4、运用重力模型获取初始OD时为什么b取值没有强制要求
反推的流量拟合与矩阵平衡的P、A约束会显著削弱b的敏感性，导致最终反推OD对b取值不敏感（因为路段、P、A的约束，反推过程把b强制修正了）
b的取值只改变迭代起点，不改变终点。但b的取值会影响收敛速度，只要不是极端取值，在以上约束条件下，最终都会走向收敛。
5、没有现状OD，之后预测过程（交通分布、诱增交通等）中重力模型参数该怎么确定
本次首先采用经验衰减系数生成重力模型初始 OD 矩阵，该参数仅作为 OD 反推的种子初值；由于 OD 反推以路段实测流量为核心约束，并辅以出行产生、吸引总量平衡，初始参数带来的分布偏差可被完全修正。最终利用校准后的现状 OD 矩阵与项目自定义阻抗矩阵，重新标定得到本地化衰减系数b，该标定值作为唯一参数应用于远期交通分布预测。
