本发明涉及客流需求预测领域，具体涉及一种高铁客流需求的预测方法。
背景技术：
：在过去的十年间，中国的高铁事业取得突飞猛进地发展，高铁建设频频刷新世界纪录，取得举世瞩目的成就。现如今，高铁已然成为中国引以为豪的一项外交名片。高铁投资建设规模、列车运行方案等一系列问题的确定都需要客流数据的支持，良好的客流预测能够为高铁的建设及其运行效能的提升提供强有力的科学依据。不同于传统的铁路，高铁运量大、运速高，在其建成初期内，客流量会在短期内急剧增长，这是传统的客流预测方法所无法反映的，因此研究专门适用于高铁客流预测方法是十分有必要的。同一交通小区内，高铁线路与公路网两种交通运输方式既共同承担区域内的客运量又相互竞争。作为舒适方便、安全性好、耗时少、准点率高的交通方式，高铁新线的通车势必会吸引大量客流量，打破区域内的交通平衡，主要体现为周边公路网的旅客流量转移到高铁线路。因此，高铁客流需求的预测与周边公路网的旅客出行OD密切相关，利用公路网旅客出行OD推算高铁客流量是合理有效的。旅客出行OD矩阵的获取主要有两种方法，一种是人工调查，也可依靠客运发车班次及其客流量数据，推算出区域客流OD矩阵，但客运发车班次往往只涉及到大客车，数据不全面，难以推算出准确的旅客出行OD，而人工调查又费时费力，耗资较大。如果通过路段观测的机动车流量反推OD分布矩阵则相对简单，较容易实施，即根据实际路网情况，通过断面交通量调查确定一个出行分布预测模型，从而实现OD出行矩阵的反推。我国的高铁客流需求预测还不成熟，多以四阶段法为基础，在铁路客流预测方面的研究主要有：王芳采用基于时间序列分析的方法对铁路短期客运量进行预测，并对比季节调整方法、Box-Jenkins方法及灰色预测法，提出基于BP神经网络的预测模型。唐丽娜综合比较移动平均法、指数平滑法和基于Matlab的AR模型法的短时客流预测效果，发现AR模型的预测效果最优。郭文建立了遗传优化的BP神经网络模型预测铁路客流量，并仿真验证了该模型的合理性。陈剑龙采用灰色模型来预测铁路年客运总量，并利用后验差方法验证模型的精确度。在高速铁路方面，秦鸣利用噪声辅助经验模态分解法进行客流预测。在考虑高速铁路短期客流周期性和波动性的基础上，杨晓提出了反映短期客流时空变化特征的改进重力模型。技术实现要素：本发明的目的是为了解决目前基于传统四阶段法的客流预测数据收集较为困难，对各种运输方式之间的竞争合作关系考虑不足等缺陷，而提出一种基于公路网客流OD反推技术的高铁客流需求预测方法。一种基于公路网客流OD反推技术的高铁客流需求预测方法，包括以下步骤：步骤一、划分交通区并构建基础公路网；步骤二、筛选机动车交通量调查断面，获得客运机动车交通量数据，统计得到交通小区i各种类型车辆的产生指标和交通小区j各种类型车辆的吸引指标；步骤三、确定反推路网阻抗矩阵；步骤四、分析整理直接影响范围内的客运发车班次数据，结合车辆额定载客人数及实载率信息，推算大型客车相应的旅客流量，形成大型客车对应的旅客OD矩阵；步骤五、基于客运机动车交通量数据中的一组断面交通量数据，进行客车OD矩阵反推，分别获得小型客车OD矩阵、中型客车OD矩阵和大型客车OD矩阵；步骤六、基于客运机动车交通量数据中的另一组断面交通量数据，对小型客车OD分布矩阵、中型客车OD分布矩阵和大型客车OD分布矩阵进行验证和修正，并乘以各车型的平均载客人数最终加和形成客流OD矩阵。进一步地，步骤一所述的划分交通区并构建基础公路网的具体过程如下：步骤1.1、交通区划分：对研究区域进行交通区划分：研究区域包括直接影响区和间接影响区；直接影响区为拟建高铁经过的区域；以地级市及行政区划级别相等的行政区域建制为单位将直接影响区划分为交通中区；以县、县级市及行政区划级别相等的行政区域建制为单位将每个交通中区划分为交通小区；间接影响区为研究区域内除直接影响区以外的外围区，以省、地级市及省、地级市行政区划级别相等的行政区域建制为单位将间接影响区划分为交通大区；步骤1.2、公路网构建：绘制直接影响区内部的公路网络，作为OD反推路网；统计公路网中所有公路的属性，包括公路名称、公路等级、设计速度、车道数、里程。进一步地，所述步骤二具体实现过程如下：对于交通小区i到交通小区j内的各条公路，选取相互独立的国、省道断面，进行断面交通量调查，获得客运机动车交通量数据；每条调查路段共有两组断面交通量信息，分别针对每组断面进行交通量统计，第一组交通量作为客流OD推算的基础数据；第二组交通量不参与OD反推的过程，仅是用于OD反推结果的检验；不同类型的机动车流量主要有两个获取途径，一是通过各条国省道上的流量观测站提供的历史数据，通过历史数据获得断面流量；另一个是另行补充调查的断面流量，即针对没有流量观测站历史数据的道路通过调查获得断面流量；二者加和即获得客运机动车交通量；针对第一组交通量数据，根据各路段的交通流量及流向，可统计得到交通小区i各种类型车辆的产生指标Pi和交通小区j各种类型车辆的吸引指标Aj；对于不同类型的车辆，产生指标Pi表示从交通小区i出发的该类型车辆的数量，吸引指标Aj表示进入交通小区j的该类型车辆的数量进一步地，所述步骤三确定反推路网阻抗矩阵的过程如下：(3.1)如果路网阻抗矩阵选择距离阻抗矩阵，距离阻抗矩阵的确定过程如下：首先确定连接两个交通小区之间所有可通行的国省道，各条国省道的长度；然后，确定两交通小区之间的最短路径及其相对次短路径，计算各路径长度；最后，对最短路径及次短路径的长度进行加权平均，得到距离阻抗矩阵；(3.2)如果路网阻抗矩阵选择时间阻抗矩阵，时间阻抗矩阵的确定过程如下：首先确定两交通小区之间的最短路径及次短路径的长度，并通过调查得到各条公路的运行速度；然后，计算两交通小区利用各条路径的时间阻抗；最后，对各路径的时间阻抗进行加权平均，得到时间阻抗矩阵；将交通小区i到交通小区j的距离或行程时间记为tij。进一步地，所述步骤四的过程如下：步骤4.1、旅客流量推算：在已知发车班次、车型、额定载客量和载客率的条件下，即可确定出旅客流量，即旅客流量＝班车数量×额定载客量×载客率；步骤4.2、旅客OD矩阵确定：针对有客运发车班次数据的交通小区，每个发车班次起点站所对应的交通小区就是该车次旅客出行的起点，即O点，终点和途经站点所对应的交通小区就是该车次旅客出行的终点，即D点；对每个发车班次整理由起点至终点或途经站点的旅客人数，即旅客发送量，填入旅客OD表中；针对无客运发车班次数据的交通小区，依据出行平衡原理，将其它交通小区到达无客运发车班次数据的某个交通小区的旅客数量作为这些交通小区的旅客发送量，填入旅客OD表中；不考虑间接影响区中的交通大区之间，交通大区与交通小区之间的旅客出行量，将交通大区之间，交通大区与交通小区之间的旅客发送量设为零；根据上述交通小区旅客发送量对应的旅客OD表得到交通小区之间的旅客OD矩阵；对交通小区旅客出行量进行合并，得交通中区的旅客OD矩阵。进一步地，所述步骤4.2中得到交通小区之间的旅客OD矩阵后需要对旅客OD矩阵进行检验，具体过程如下：公路客运量是指在一定时期内各种公路客运交通工具实际运送的旅客人数；在统计时，无论出行远近或者票价多少均按一人一次客运量统计：各交通小区的客运量为公路运输企业报告期内实际运送的旅客人数之和；依据各交通小区发布的公路客运量数据对旅客OD矩阵进行检验，用以核查载客率取值是否合适，如果旅客OD矩阵与在公路客运量的数据相差小于10％，则表明载客率取值合适，推算的旅客OD矩阵准确；如果两者相差大于等于10％，则返回步骤4.1需重新调查相关地点客运发车班次的载客率，并计算旅客流量，直至旅客OD矩阵与在公路客运量的数据相差小于10％。进一步地，所述步骤五的过程如下：基于第一组断面交通量数据，进行客车OD矩阵反推：(5.1)针对小型客车和中型客车，采用无先验OD矩阵的OD反推，确定小型客车OD矩阵和中型客车OD矩阵，分别作为小型客车OD分布矩阵和中型客车OD分布矩阵；在无以往OD分布矩阵时，交通小区之间的出行交换量能用重力模型进行描述，即式中：qij为交通小区i到交通小区j的OD分布量；Pi为交通小区i的产生指标；Aj为交通小区j的吸引指标；tij为交通小区i到交通小区j的行程时间或距离；k，γ分别为重力模型参数；在进行OD分配时，采用容量限制迭代平衡的方法对OD分布量在路段上进行反复的分配计算，以求解路径选择概率及推算交通量；在迭代过程中，将OD分布量分配到某一路段上的公式为：式中：Qa为路段a上的分配交通量；为交通小区i到交通小区j的OD分布量选用路段a的概率；(5.2)针对大型客车，采用无先验OD与有先验OD相结合的OD反推方法：首先采用与上述小型客车和中型客车无先验OD反推确定小型客车OD矩阵和中型客车OD矩阵相同的方法，推算出大型客车的无先验OD矩阵，即大型客车OD矩阵；然后，依据步骤四中的大型客车对应的旅客OD矩阵对无先验OD反推的大型客车OD矩阵进行修正，修正过程为：大型客车OD矩阵乘以平均载客人数后，与大型客车的旅客OD矩阵进行比较，对于相差10％以上的数据，将大型客车OD矩阵中对应值修正为旅客OD矩阵除以平均载客人数；如果数据数据符合条件，不进行修正；得到最终的大型客车OD矩阵；将最终的大型客车OD矩阵作为先验OD矩阵；最后，采用最大熵值法进行有先验OD矩阵时的OD反推，确定大型客车OD分布矩阵。进一步地，所述步骤六的过程如下：基于两组断面进行机动车交通量统计中的另一组数据交通量进行大型客车、中型客车、小型客车的OD分布矩阵的精度验证：将步骤五得到的小型客车、中型客车、大型客车对应的OD分布矩阵重新分配到路网上，得到分配后的各路段交通量信息；以分配后的交通量与步骤二得到的第二组断面交通量间的相对误差作为评价客车OD反推精度验证指标；若相对误差大于误差阈值，则需对OD分布矩阵进行修正；当相对误差小于误差阈值，则表明反推的OD分布矩阵符合真实情况；依据小型客车OD分布矩阵、中型客车OD分布矩阵、大型客车OD分布矩阵及各车型的平均载客人数，可得到各种车型对应的客流OD矩阵，加和可得区域总客流OD矩阵。进一步地，所述步骤六中若相对误差大于误差阈值，则需对OD分布矩阵进行修正，修正过程如下：分别针对各种车型对应的OD分布矩阵命名为OD1，其交通分配后得到各条路段流量，其中第n个路段的流量为qn1，与第二组断面交通量qn2比较，相对误差大于误差阈值，则取路段n的交通量qn’＝(qn1+qn2)/2，按步骤五重新反推得到OD分布矩阵，记为OD3，再进行交通分配得qn3，验算其与qn1的相对误差，若仍不满足要求，则取交通量qn’＝(qn3+qn1)/2，按步骤五重新重新反推得到OD分布矩阵，记为OD4，与第二组断面交通量qn2进行相对误差比较，如此循环迭代下去直至误差满足精度要求为止。本发明的有益效果为：1、本发明解决现有的OD调查中客流预测数据收集困难的问题，通过断面流量调查实现OD矩阵的反推，不仅仅方便简单，而且大大减少了工作量，更节省了费用。以实施例中条件进行调查，现有的OD调查需要拦车问询，费时又费力，在两个月的时间内完成OD调查至少需要200人，而本发明仅需观测交通流量的大小，在同样的时间内完成数据收集仅需40人左右，不仅极大简化了数据收集，而且极大的减少了工作量，减少工作量至少60％以上。2、本发明充分考虑公路网与高铁线路的竞争关系，通过公路网旅客OD反推实现高铁客流预测，相比现有的方法，数据相对更加准确。附图说明图1为本发明的流程图；图2为交通中区和交通大区划分示意图；图3为交通小区划分示意图。具体实施方式具体实施方式一：结合图1说明本实施方式，一种基于公路网客流OD反推技术的高铁客流需求预测方法，包括以下步骤：步骤一、划分交通区并构建基础公路网；步骤二、筛选机动车交通量调查断面，获得客运机动车交通量数据，统计得到交通小区i各种类型车辆的产生指标和交通小区j各种类型车辆的吸引指标；步骤三、确定反推路网阻抗矩阵；步骤四、分析整理直接影响范围内的客运发车班次数据，结合车辆额定载客人数及实载率信息，推算大型客车相应的旅客流量，形成大型客车对应的旅客OD矩阵；步骤五、基于客运机动车交通量数据中的一组断面交通量数据，进行客车OD矩阵反推，分别获得小型客车OD矩阵、中型客车OD矩阵和大型客车OD矩阵；步骤六、基于客运机动车交通量数据中的另一组断面交通量数据，对小型客车OD分布矩阵、中型客车OD分布矩阵和大型客车OD分布矩阵进行验证和修正，并乘以平均载客人数最终加和形成客流OD矩阵。具体实施方式二：本实施方式中步骤一所述的划分交通区并构建基础公路网的具体过程如下：步骤1.1、交通区划分：对研究区域进行交通区划分：研究区域包括直接影响区和间接影响区；直接影响区为拟建高铁经过的区域；以地级市及行政区划级别相等的行政区域建制为单位将直接影响区划分为交通中区，如图2中的编号1-8所示的区域；以县、县级市及行政区划级别相等的行政区域建制为单位将每个交通中区划分为交通小区，如图3中的编号1至49所示的区域；间接影响区为研究区域内除直接影响区以外的外围区，根据实际情况以省、地级市及省、地级市行政区划级别相等的行政区域建制为单位将间接影响区划分为交通大区，如图2中的编号9-16所示的区域；步骤1.2、公路网构建：绘制直接影响区内部的公路网络，作为OD反推路网；由于高铁线路主要是提供长距离的客运出行，因此OD反推主要推算的是长距离的公路客运出行，故客流OD反推路网就是由直接影响区内部及周边国道和省道构成的干线公路网。统计公路网中所有公路的属性，包括公路名称、公路等级、设计速度、车道数、里程。其他步骤和参数与具体实施方式一相同。具体实施方式三：本实施方式所述步骤二具体实现过程如下：对于交通小区i到交通小区j内的各条公路，选取相互独立的国、省道断面，进行断面交通量调查，获得客运机动车交通量数据，这是OD反推时的依据交通量；在选择路段交通量调查断面时应当遵循以下四个基本原则：①OD覆盖原则：调查路段应包含所有OD点对的信息，即每个OD点对至少被一条路段检测到。②路段独立原则：保证所选断面是线性不相关的。③最大流截断原则：在调查断面数量相同的前提下，选择的断面应尽可能多地截断交通流。④最大流量比原则：该原则是针对单个OD点对而言的，即对某一特定的OD点对，应使得调查路段的交通量占该OD点对流量的比例尽可能大。每条调查路段共有两组断面交通量信息，分别针对每组断面进行交通量统计，第一组交通量作为客流OD推算的基础数据；第二组交通量不参与OD反推的过程，仅是用于OD反推结果的检验；后续计算涉及的交通量来源于第一组交通量；不同类型的机动车流量主要有两个获取途径，一是通过各条国省道上的流量观测站提供的历史数据，通过历史数据获得断面流量；另一个是另行补充调查的断面流量，即针对没有流量观测站历史数据的道路通过调查获得断面流量；二者加和即获得客运机动车交通量；流量观测站的历史数据一般含以下七种车型：小型货车、中型货车、大型货车、特大型货车、集装箱车、小型客车、中型客车和大型客车，非高速公路的观测断面还包括摩托车和拖拉机两种车型；交通量OD推算的主要对象为客运机动车，即大、中、小型客车，因此各类客运车辆的交通量及路段流量(标准车)为数据统计的重点。故实际调查，仅需调查小型客车、中型客车、大型客车的交通量即可。针对第一组交通量数据，根据各路段的交通流量及流向，可统计得到交通小区i各种类型车辆(小型客车、中型客车和大型客车)的产生指标Pi和交通小区j各种类型车辆的吸引指标Aj；对于不同类型的车辆，产生指标Pi表示从交通小区i出发的该类型车辆的数量，吸引指标Aj表示进入交通小区j的该类型车辆的数量。其他步骤和参数与具体实施方式一或二相同。具体实施方式四：本实施方式所述步骤三确定反推路网阻抗矩阵的过程如下：交通阻抗是反映交通区间交通便利程度的指标，是对交通区间交通设施状况和交通工具状况的综合反映。可用距离、时间或费用等表示，距离阻抗、时间阻抗及费用阻抗在一定条件下可相互转换。比较费用阻抗，距离阻抗和时间阻抗矩阵获取较易且精度可以保证，故推算路网的距离阻抗矩阵和时间阻抗矩阵。(3.1)如果路网阻抗矩阵选择距离阻抗矩阵，距离阻抗矩阵的确定过程如下：首先确定连接两个交通小区之间所有可通行的国省道，各条国省道的长度；然后，确定两交通小区之间的最短路径及其相对次短路径，计算各路径长度；最后，对最短路径及次短路径的长度进行加权平均，得到距离阻抗矩阵；(3.2)如果路网阻抗矩阵选择时间阻抗矩阵，时间阻抗矩阵的确定过程如下：首先确定两交通小区之间的最短路径及次短路径的长度，并通过调查得到各条公路的运行速度；然后，计算两交通小区利用各条路径的时间阻抗；最后，对各路径的时间阻抗进行加权平均，得到时间阻抗矩阵；将交通小区i到交通小区j的距离或行程时间记为tij。各条路径时间阻抗的计算采用美国联邦公路局的BPR函数，其表达式如下：式中的Va和Ca分别代表路段a的交通量和通行能力，ta(0)为路段自由行程时间；α、β分别为系数。本次OD反推过程中，为使获得的客运OD矩阵更加准确，将小型客车、中型客车及大型客车分别进行OD反推，即分别采用路段的小型客车、中型客车及大型客车数据进行OD反推。Va和Ca是路段a各种车型对应的交通量和通行能力，本发明采用如下方法来确定各类车型进行OD反推时的路段通行能力：将各类车型折算为标准小汽车的折算系数见表1。表1机动车折算系数假设路段a上各类型车的比例见表2。表2路段a的车型比例则路段a上各类客运车辆的通行能力分别为：式中：Cai′为路段a上i′类客运车辆的通行能力，i′＝1,2,3分别代表小型客车、中型客车和大型客车；γi′为路段a上i′类客运车辆的比例；PCEi′为i′类客运车辆的折算系数。其他步骤和参数与具体实施方式一至三之一相同。具体实施方式五：本实施方式所述步骤四的过程如下：步骤4.1、旅客流量推算：由于交通小区i到交通小区j之间的客运车辆均为大型客车，故这里推算的旅客流量为大型客车的旅客流量。在已知发车班次、车型、额定载客量和载客率的条件下，即可确定出旅客流量，即旅客流量＝班车数量×额定载客量×载客率；客运发车班次、车型和时刻表等信息的收集可通过现场走访各县市的公路长途客运站及互联网查询认证等方式。主要收集直接影响区内的客运数据，间接影响区由于涵盖范围较大，故不作考虑。额定载客量，即每辆车的额定载客人数或座位数，可通过客运发车班次的车型确定。载客率的确定可通过一个小规模的抽样调查完成。步骤4.2、旅客OD矩阵确定：针对有客运发车班次数据的交通小区，每个发车班次起点站所对应的交通小区就是该车次旅客出行的起点，即O点，终点和途经站点所对应的交通小区就是该车次旅客出行的终点，即D点；依据上述原理，对每个发车班次整理由起点至终点或途经站点的旅客人数，即旅客发送量，填入旅客OD表中；针对无客运发车班次数据的交通小区，依据出行平衡原理(旅客发送量与到达量大体相当的原理)，将除无客运发车班次数据的某个交通小区以外的其它交通小区到达无客运发车班次数据的某个交通小区的旅客数量(即D点对应的旅客数量)作为这些交通小区的旅客发送量(即O点对应的旅客数量)，填入旅客OD表中；不考虑间接影响区中的交通大区之间，交通大区与交通小区之间的旅客出行量，将交通大区之间，交通大区与交通小区之间的旅客发送量设为零；根据上述交通小区旅客发送量对应的旅客OD表得到交通小区之间的旅客OD矩阵；对交通小区旅客出行量进行合并，得交通中区的旅客OD矩阵。其他步骤和参数与具体实施方式一至四之一相同。具体实施方式六：本实施方式所述步骤4.2中得到交通小区之间的旅客OD矩阵后需要对旅客OD矩阵进行检验，具体过程如下：公路客运量是指在一定时期内(可取一年)各种公路客运交通工具实际运送的旅客人数；在统计时，无论出行远近或者票价多少均按一人一次客运量统计：各交通小区的客运量一般为公路运输企业报告期内实际运送的旅客人数之和；依据各交通小区相关部门发布的公路客运量数据对旅客OD矩阵进行检验，用以核查载客率取值是否合适，如果旅客OD矩阵与在公路客运量的数据相差小于10％，则表明载客率取值合适，推算的旅客OD矩阵准确；如果两者相差大于等于10％，则返回步骤4.1需重新调查相关地点客运发车班次的载客率，并计算旅客流量，直至旅客OD矩阵与在公路客运量的数据相差小于10％。其他步骤和参数与具体实施方式五相同。具体实施方式七：本实施方式所述步骤五的过程如下：基于第一组断面交通量数据，进行客车OD矩阵反推：OD反推是交通分配的逆过程。交通分配是已知交通分布(OD矩阵)的前提下，按照一定的规则将其分配道路网中的不同路径上，得到路段流量；而OD反推，就是根据现状路网的情况，通过观测路段上的机动车流量来确定一个出行分布模型(无先验OD的情况)或对已有的早期OD(或小样本调查得到的OD)分布矩阵进行修正(有先验OD的情况)。(5.1)针对小型客车和中型客车，采用无先验OD矩阵的OD反推，确定小型客车OD矩阵和中型客车OD矩阵，分别作为小型客车OD分布矩阵和中型客车OD分布矩阵；在无以往OD分布矩阵时，交通小区之间的出行交换量能用重力模型进行描述，即式中：qij为交通小区i到交通小区j的OD分布量(拟通过计算确定的未知量)；Pi为交通小区i的产生指标；Aj为交通小区j的吸引指标；tij为交通小区i到交通小区j的行程时间或距离(由路网确定，为已知量)；k，γ分别为重力模型参数；显然，在无以往OD分布矩阵时，OD反推的核心就是标定出一个重力模型。在进行OD分配时，采用容量限制迭代平衡的方法对OD分布量在路段上进行反复的分配计算，以求解路径选择概率及推算交通量；在迭代过程中，将OD分布量分配到某一路段上的公式为：式中：Qa为路段a上的分配交通量；为交通小区i到交通小区j的OD分布量选用路段a的概率；由于Qa为已知量，只需采用相应数学方法标定重力模型参数k及γ即可。(5.2)针对大型客车，采用无先验OD与有先验OD相结合的OD反推方法：首先采用与上述小型客车和中型客车无先验OD反推确定小型客车OD矩阵和中型客车OD矩阵相同的方法，推算出大型客车的无先验OD矩阵，即大型客车OD矩阵；然后，依据步骤四中的大型客车对应的旅客OD矩阵对无先验OD反推的大型客车OD矩阵进行修正，修正过程为：大型客车OD矩阵乘以平均载客人数后，与大型客车的旅客OD矩阵进行比较，对于相差10％以上的数据，将大型客车OD矩阵中对应值修正为旅客OD矩阵除以平均载客人数；如果数据数据符合条件，不进行修正；得到最终的大型客车OD矩阵；将最终的大型客车OD矩阵作为先验OD矩阵；最后，采用最大熵值法进行有先验OD矩阵时的OD反推，确定大型客车OD分布矩阵。由于是在先验OD历史数据上进行修正，其推算结果一般远远高于无先验OD的情况，因此，在条件允许的情况下，应尽量选用历史OD作为先验矩阵。其他步骤和参数与具体实施方式一至六之一相同。具体实施方式八：本实施方式所述步骤六的过程如下：基于两组断面进行机动车交通量统计中的另一组数据交通量进行大型客车、中型客车、小型客车的OD分布矩阵的精度验证：将步骤五得到的小型客车、中型客车、大型客车对应的OD分布矩阵重新分配到路网上，得到分配后的各路段交通量信息；以分配后的交通量与步骤二得到的第二组断面交通量间的相对误差作为评价客车OD反推精度验证指标；若相对误差大于误差阈值，则需对OD分布矩阵进行修正；当相对误差小于误差阈值，则表明反推的OD分布矩阵符合真实情况；依据小型客车OD分布矩阵、中型客车OD分布矩阵、大型客车OD分布矩阵及各车型的平均载客人数，可得到各种车型对应的客流OD矩阵，加和可得区域总客流OD矩阵。其他步骤和参数与具体实施方式一至七之一相同。具体实施方式九：本实施方式所述步骤六中若相对误差大于误差阈值，则需对OD分布矩阵进行修正，修正过程如下：分别针对各种车型(小型客车、中型客车、大型客车)对应的OD分布矩阵命名为OD1，其交通分配后得到各条路段流量，其中第n个路段的流量为qn1，与第二组断面交通量qn2比较，相对误差大于误差阈值，则取路段n的交通量qn’＝(qn1+qn2)/2，按步骤五重新反推得到OD分布矩阵，记为OD3，再进行交通分配得qn3，验算其与qn1的相对误差，若仍不满足要求，则取交通量qn’＝(qn3+qn1)/2，按步骤五重新重新反推得到OD分布矩阵，记为OD4，与第二组断面交通量qn2进行相对误差比较，如此循环迭代下去直至误差满足精度要求为止。其他步骤和参数与具体实施方式八相同。实施例：本实施例一种基于公路网客流OD反推技术的高铁客流需求预测方法具体是按照以下步骤实现的：具体实施方式以郑州至合肥区域公路网现状客流OD推算为例。依据本研究的目的并结合断面流量调查和客运发车班次基础数据，对公路网客流OD进行推算。研究范围为郑州至合肥间的2个省会城市和6个地级市，外围区、交通中区、交通小区划分结果分别见表2、表3、表4。表2外围区一览表表3交通中区一览表表4交通小区一览表小区编号小区名称小区编号小区名称小区编号小区名称1郑州市区18鄢陵县35颍上县2新郑市19襄城县36淮南市区3登封市20周口市区37凤台县4新密市21项城市38六安市区5巩义市22扶沟县39寿县6荥阳市23西华县40霍邱县7中牟县24商水县41舒城县8开封市区25太康县42金寨县9开封县26鹿邑县43霍山县10尉氏县27郸城县44合肥市区11通许县28淮阳县45长丰县12杞县29沈丘县46肥东县13兰考县30阜阳市区47肥西县14许昌市区31界首市48巢湖市15禹州市32临泉县49庐江县16长葛市33太和县17许昌县34阜南县由表2至表4可知，外围区基本上覆盖了项目直接影响区外的河南、安徽两省剩余地市及全国其他省市(不含港澳台地区)。交通中区划分至地级市，交通小区则细分到县级。OD反推由郑州至合肥间的81条国省干线公路构成，其中，国道30条(高速公路14条，普通国道16条)，省道51条(高速公路17条，普通省道34条)。反推路网中各条公路的名称、行政等级、技术等级、设计速度、车道数、里程等属性见表5。表5OD反推路网一览表主要依据最大流截断原则，即调查路段的选择以高等级公路(高速、一级公路)为主、低等级公路为辅(二级及以下等级公路)，最终选择了242个流量调查断面，其中，河南省210个断面，安徽省32个断面。调查断面数约占反推路网路段总数的40.5％，远高于最低10％的要求。各断面的机动车交通量调查结果见表6。表6OD机动车交通量一览表由表6可知，242个断面流量合计为428.7万辆标准小汽车，其中客运车辆164.2万辆标准小汽车。断面流量从高速公路、普通国道至普通省道呈递减趋势分布，平均日交通量分别为23333辆标准小汽车、14327辆标准小汽车和9260辆标准小汽车。高速公路、普通国道和普通省道小型客车、中型客车和大型客车的断面流量依次为8013辆、201辆、739辆，4228辆、108辆、400辆和3167辆、88辆、257辆。通过现场走访各县市的公路长途客运站并通过互联网查询认证后，收集了直接影响区内38个县市(即交通小区)的日公路客运发车班次、车型及时刻表等信息。剩余登封市、中牟县等11个县市的客运发车班次尚未获得，需作补遗处理。通过一个小规模的抽样调查确定了载客率的取值标准如下：到达站为终点站且日发车班次小于15次/日时，载客率为0.55～0.70；到达站为终点站且日发车班次大于15次/日时，载客率为0.70～0.75；到达站为途经站时，载客率为0.30。根据客运发车班次及载客率获得旅客流量，从而确定出旅客出行OD矩阵见表7。由于交通小区较多，这里仅列出前10个小区的出行OD作为示意。表7旅客出行OD矩阵OD1234567891011004913312618099581329175723063512731530557542212013598028520320107003618000196905807700481998219690011790003151424000042038380631932581179421553232130730635000032000082996107770383200080695305000381300001079970031000486800根据安徽省公路运输管理局发布的《公路客运调查技术方案》中公路客运量数据，对上述OD矩阵进行检验，核查客运车辆载客率取值是否合适。对比结果见表8。表8公路客运量与发车班次客运量的对比由表8可知，在9个核查市县中，6个市县发送量的误差在5％以内，3个在10％以内，说明客运发车班次的载客率取值是比较合适的。根据道路长度、运行速度等参数确定路网的距离阻抗矩阵和时间阻抗矩阵，见表9、表10，同样仅示意前10个交通小区。表9距离阻抗矩阵12345678910103559226823135176902350703612263236186563597002951648011214512942236290733138779695568122517304289129155171623636431420387810111971323803889380366073851611127712978360245397686145961551016024060109056129951711197353600表10时间阻抗矩阵大客车采用有先验OD矩阵的OD反推技术，中客车和小客车采用无先验OD矩阵的OD反推技术，分别得到三种车型的OD矩阵见表11至表13。表11大客车出行OD矩阵OD1234567891010158861017511112127563482158014116107130386104918113711410144901681461295751110141447161611118168140157117121611141028118275107374647572806542963111161165011048301291114210表12中客车出行OD矩阵OD123456789101026507910761621014446226012311231283491016115161146611601111617520111015209151621151115016117622112010611824366966039694211115113901106528171114010表13小客车出行OD矩阵在已获得交通量信息的242个断面交通量中，选取一组断面交通量用于客车OD矩阵反推，剩余的断面交通量用于反推OD的精度验证。小客车、中客车、大客车OD反推精度验证结果见表14。表14小客车、中客车、大客车OD反推精度验证结果注：表中流量均为道路双向流量，单位为veh/day。由表14可知，小客车与中客车分配流量与实际流量平均相对误差约为10％，其中70％的路段相对误差在10％以内；大客车OD反推过程中由于引入先验OD(小客车、中客车反推过程中不具备构造先验OD的条件)，其反推OD的精度进一步提高，路段分配流量与实际流量平均相对误差约为5％，其中80％的路段相对误差在10％以内。依据小客车、中客车、大客车现状出行OD矩阵及各车型的平均载客人数，可得到三种车型对应的客流OD矩阵，进而可得区域总客流OD矩阵见表15。表15区域客流出行OD矩阵本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。当前第1页1 2 3