CN104715605A - 一种基于vsp分布的交通运行数据与排放数据耦合方法与系统 - Google Patents

Abstract

本发明提出了一种基于VSP分布的交通运行数据与机动车排放数据间的耦合方法与系统，通过对交通运行数据与机动车排放数据的特征进行挖掘与分析，以VSP分布为中介变量，建立了交通流量、平均速度、机动车行驶工况等交通运行参数与排放率、排放因子等机动车排放参数之间的动态映射关系。该发明简化了机动车尾气排放的物理意义，提高了交通排放量化评估的效率与准确性。

Description

一种基于VSP分布的交通运行数据与排放数据耦合方法与系统

【技术领域】

本发明涉及一种基于VSP分布的交通运行与排放数据耦合方法与系统。

【背景技术】

交通部门与环保部门经过长期的建设与发展，已经形成各自较为完备的数据采集与分析系统，并且已经获得了大量数据。

研究一个地区的交通运行情况，首先需要经过交通调查，了解该地区的交通特征及其发展趋势。通过对交通运行数据进行定量分析，掌握各种特征参数的具体数据，为针对具体情况进行科学的交通规划、线性设计和交通管理提供理论依据。常见的交通调查方法有定点测量(Point Method)、沿路段测量(Along a Length Method)、浮动车测量(Moving Observe)、短距离测量法(Short Section Method)等。感应线圈检测技术、超声波检测技术、红外线检测技术、移动信息采集技术等随着交通控制或交通信息服务的需求与20世纪60年代开始出现，如今获得了普遍应用。远程微波检测技术、视频图像式检测技术和光信标式检测技术等于20世纪80年代末研制成功，90年代获得应用。21世纪初，人们又开始利用GPS检测技术检测交通流参数。

在制定交通策略以降低机动车尾气排放的过程中，利用实验手段进行尾气排放检测，收集机动车实际的油耗排放数据，是一切有关机动车油耗排放研究的基础，也是建立油耗排放模型和评价交通策略的基本依据。目前国内外主要有四种机动车油耗和尾气排放检测方法，分别是台架测试法、隧道测试法、遥感检测法和车载测试法。

但是受限于采集方式与应用目的的不同，上述两类数据的相关特征未进行深入挖掘，无法实现准确、有效地关联。

【发明内容】

针对现有技术存在的问题，本发明提出了一种交通运行数据与机动车排放数据间的耦合方法与系统，通过对交通运行数据与机动车排放数据的特征进行挖掘与分析，以VSP分布为中介变量，建立了交通流量、平均速度、机动车行驶工况等交通运行参数与排放率、排放因子等机动车排放参数之间的动态映射关系。考虑到机动车在实际道路上的行驶速度对排放的影响，结合既有交通运行参数数据与机动车排放测试结果建立交通运行与机动车排放数据库，利用数据集群和并行计算技术，实现对交通运行与排放数据库的快速反应。同时集成GPS、GIS和数据挖掘技术，实现交通运行数据与机动车排放数据的耦合，对路段排放率进行测算。该发明简化了机动车尾气排放的物理意义，提高了交通排放量化评估的效率与准确性。

【附图说明】

图1交通运行数据与排放数据耦合方法设计思路

图2轻型车行驶工况数据库示例

图3轻型车在快速路VSP分布

图4机动车油耗和尾气排放实测数据示例

图5机动车排放率数据库示例

图60-4万公里国III标准轻型汽油车排放因子示例

图7机动车排放因子数据库示例

图8交通运行数据与排放数据耦合系统

【具体实施方式】

下面通过具体的实施例并结合附图对本发明作进一步详细的描述。

本发明基于VSP建立油耗和排放微观统计模型，能够弥补传统尾气预测模型对于交通特征考虑的不足，通过大量实测数据，利用VSP对北京市各类型道路上的驾驶行为进行描述。综合考虑北京市路段平均长度及路段平均行驶时间，以3分钟为间隔计算平均速度，建立不同等级道路上各平均速度区间的VSP分布，将车辆的平均行驶速度与微观瞬时行驶状态相结合，以充分反映不同行驶状态下的速度对于油耗和排放的影响。建立以基本交通参数为输入变量的油耗和排放的速度修正模型，实现交通运行数据与机动车尾气排放数据的有效衔接，从而为交通排放量化评估提供准确有效的手段。发明的设计思路如图1所示。

本发明的基于VSP分布的交通运行数据与排放数据耦合方法，包括以下步 骤：

第一步，采集交通运行数据，进一步包括：

步骤1.1，数据采集，利用车载GPS设备，采集机动车逐秒行驶工况数据，数据主要字段包括采集日期、时间、经度、纬度、速度、方向角、高程等，图2为实际采集到的轻型车逐秒行驶工况数据示例，将数据按照统一的格式进行整理，通过初步质量控制，去除不符合要求的记录，建立机动车行驶工况数据库。

步骤1.2，进行车型划分，将机动车按车重进行分类，分为微型车、轻型车、中型车和重型车。

步骤1.3，进行道路类型匹配。我国城市道路依据通行能力与交通功能划分为快速路、主干路、次干路、支路四种类型。道路类型的匹配工作在GIS系统中完成。每一条工况数据都包含经纬度字段，据此可以在已有的北京路网上创建点，通过判断该点所处的路段名得到其道路类型。经过道路类型匹配之后，所有机动车行驶工况数据依据其所属道路类型被划分为四个部分。

步骤1.4，进行短行程划分。对于每一种车型与道路类型，所有数据被划分到时长为180秒的短行程。这样处理是因为较长的行程时间，比如固定的行驶周期，由于包含多种行驶状态而无法刻画机动车在现实道路上的行驶特征，选择180秒作为瞬时速度的集成粒度既可以包含机动车在基本路段的行驶特征，又可以包含其在交叉口的行驶特征。每一个短行程都需要计算其平均速度，以便进行下一步的速度聚类。

步骤1.5，进行平均速度聚类，将上述过程得到的所有短行程依据平均速度进行聚类，得到其所属的速度区间，方法如下式所示：

步骤1.6，进行VSP聚类，VSP由瞬时速度和加速度计算得到，因为采集得到的行驶工况数据包含连续的逐秒速度，瞬时加速度由相邻两秒的速度差得出。VSP的计算如下式所示：

$\mathrm{VSP}=\frac{\mathrm{Av}+{\mathrm{Bv}}^2+{\mathrm{Cv}}^3+\mathrm{mva}}{f}$

其中，v为机动车速度，单位m/s；a为加速度，单位m/s²；A、B、C、m、f为常数。

为了可以更好地反映VSP分布与平均速度之间的关系，以1kW/t为间隔对 VSP进行划分并聚类，得到VSP bin，具体计算见下式：

$\∀:\mathrm{VSP}\∈[n-0.5,n+0.5),\mathrm{VSPbin}=n$

步骤1.7，建立VSP分布数据库，以大量不同车型逐秒的行驶工况数据为基础，经过上述处理过程，最后对所有车型、道路类型、速度区间相同的VSP进行聚类，统计每一个速度区间对应的VSP个数占该类VSP总数的百分比，即VSP分布，如下式所示：

$R_{i,k}=\frac{N_{i,k}}{N_k}$

其中，R_i,k为第k个平均速度区间第i个VSP bin的分布率；N_k为第k个平均速度区间的VSP总数；N_i,k为第k个平均速度区间VSP bin为i的VSP个数。图3为轻型车在快速路VSP分布示例。

第二步，采集机动车排放数据，进一步包括：

步骤2.1，采集排放测试数据，利用PEMS车载尾气测试设备采集机动车油耗和尾气排放数据，数据主要字段包括采集日期、时间、速度、NO_X、HC、CO逐秒排放量、逐秒油耗量等。图4为收集到的北京市实际路网中按照正常交通流行驶车辆的油耗和尾气排放数据示例。

步骤2.2，划分燃料类型，燃料类型包括：汽油、柴油、液化石油气(Liquefied Petroleum Gas，LPG)、液化天然气(Liquefied Natural Gas，LNG)、压缩天然气(Compressed Natural Gas，CNG)、混合动力、双燃料。

步骤2.3，划分排放标准，为有效的控制机动车污染和排放，按照车辆所执行的排放标准对车辆进行分类，根据测试车辆的生产年份，确定车辆的排放标准类型，划分车辆排放标准。表1是汽油车各阶段排放标准实施日期，对于在国I标准实施年份之前生产的车辆，将其定义为国0标准。

表1汽油车各阶段排放标准实施日期

排放标准	实施时间
		国I	1999.1.1
国II	2002.8.1
		国III	2005.12.30
国IV	2008.3.1
		国V	2013.2.1

步骤2.4，划分行驶里程，对同类车型的行驶里程进行划分，表2为轻型车 行驶里程划分方式。

表2机动车行驶里程划分方式

类别	行驶里程(万公里)
		1	[0，4)
2	[4，8)
		3	[8，16)
4	[16，24)
		5	[24，+∞)

步骤2.5，计算平均排放率，计算每种车型、燃料类型、排放标准、行驶里程相应分类的每个VSP bin的平均排放率，如下式所示：

${\mathrm{ER}}_i=\frac{1}{m}\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{er}}_j$

其中，ER_i为第i个VSP bin的平均排放率，单位g/s；er_j是VSP bin为i的第j个VSP对应的排放率，单位为g/s；m是VSP bin为i的VSP个数。

步骤2.6，建立排放率数据库，将不同车型、排放标准、行驶里程的上述计算结果整理、分类，建立各VSP bin排放率数据库。图5为排放率数据库示例。

第三步，对交通运行数据和排放数据进行耦合，进一步包括：

步骤3.1，计算排放因子。排放因子是机动车尾气排放测算中的重要参数，能够直观反映机动车的排放水平。排放因子定义为机动车行驶单位距离排放的污染物质量，各平均速度区间排放因子的计算方法如下式所示：

${\mathrm{EF}}_k=(\underset{i}{\mathrm{\Σ}}{\mathrm{ER}}_i\×\mathrm{VSP}{\mathrm{bin}}_i)/v\×3600$

其中，EF_k为第k平均速度区间的排放因子，单位为g/km；ER_i是第i个VSP bin的平均排放率，单位为g/s；VSP bin_i是第k平均速度区间第i个VSP bin的分布值；v是第k平均速度区间的中值，单位为km/h。

通过VSP分布与VSP bin平均排放率的映射关系，可以得到不同车型、道路类型、燃料类型、排放标准、行驶里程、平均速度区间的排放因子。图6为行驶里程为0-4万公里的国III标准轻型汽油车在四种道路类型、不同速度区间的NO_X、CO、HC排放因子示例。

步骤3.2，计算油耗因子。油耗因子定义为机动车行驶单位距离消耗的燃油 质量，计算方法如下式所示：

${\mathrm{EF}}_{\mathrm{Fuel}}=\frac{1}{\%C}(\frac{12}{44}{\mathrm{EF}}_{{\mathrm{CO}}_2}+\frac{12}{28}{\mathrm{EF}}_{\mathrm{CO}}+\frac{12}{13}{\mathrm{EF}}_{\mathrm{HC}})$

其中，EF_Fuel为油耗因子；EF_CO、EF_HC分别为CO₂、CO、HC排放因子；％C为燃料中C的质量比例。

对于实测数据无法覆盖到的分类情况，可以利用排放预测模型得到其相应的油耗因子。整理上述计算结果，建立车型、道路类型、平均速度、行驶里程、燃料类型、排放标准机动车数据库CO2、CO、HC、NO_X油耗因子。图7为机动车油耗因子数据库示例。

第四步，路段排放量计算，对城市道路网或重要交通走廊机动车尾气排放进行测算，评估其尾气排放强度，为制定合理的减排策略与减排目标提供科学依据。路段排放量计算方法如下式所示：

${\mathrm{ER}}_{\sec \mathrm{tion}}=l\·\underset{i,j,k,m}{\mathrm{\Σ}}{\mathrm{EF}}_{\mathrm{ijkm}}^v\·q({\mathrm{Vehicle}}_i,{\mathrm{Fuel}}_j,{S\tan \mathrm{dard}}_k,{\mathrm{Mileage}}_m)$

其中，ER_section为测算路段排放率，l为测算路段长度，q(Vehicle_i,Fuel_j,Standard_k,Mileage_m)为第i种车型、第j种燃料类型、第k种排放标准、第m类行驶里程的交通量，为上述分类对应的速度区间为v的机动车排放因子。

本发明的基于VSP分布的交通运行数据与排放数据耦合的系统，可以对交通运行数据、机动车排放数据格式化处理，将交通运行数据与机动车排放数据进行耦合，对路段排放率进行测算。系统结构如图8所示。该系统包括：交通运行数据采集单元，机动车排放数据采集单元，数据耦合单元，路段排放量计算单元。其中：

交通运行数据采集单元对交通运行数据进行采集、分析和处理，最终建立VSP分布数据库。进一步包括：

工况数据采集模块，利用车载GPS设备，采集机动车逐秒行驶工况数据，数据主要字段包括采集日期、时间、经度、纬度、速度、方向角、高程等，将数据按照统一的格式进行整理，通过初步质量控制，去除不符合要求的记录，建立机动车行驶工况数据库。

车型划分模块，将机动车按车重进行分类，分为微型车、轻型车、中型车和重型车。

道路类型匹配模块，我国城市道路依据通行能力与交通功能划分为快速路、主干路、次干路、支路四种类型。道路类型的匹配工作在GIS系统中完成。每 一条工况数据都包含经纬度字段，据此可以在已有的北京路网上创建点，通过判断该点所处的路段名得到其道路类型。经过道路类型匹配之后，所有机动车行驶工况数据依据其所属道路类型被划分为四个部分。

短行程划分模块，对于每一种车型与道路类型，所有数据被划分到时长为180秒的短行程。这样处理是因为较长的行程时间，比如固定的行驶周期，由于包含多种行驶状态而无法刻画机动车在现实道路上的行驶特征，选择180秒作为瞬时速度的集成粒度既可以包含机动车在基本路段的行驶特征，又可以包含其在交叉口的行驶特征。每一个短行程都需要计算其平均速度，以便进行下一步的速度聚类。

平均速度聚类模块，将上述过程得到的所有短行程依据平均速度进行聚类，得到其所属的速度区间，方法如下式所示：

VSP聚类模块，VSP由瞬时速度和加速度计算得到，因为采集得到的行驶工况数据包含连续的逐秒速度，瞬时加速度由相邻两秒的速度差得出。VSP的计算如下式所示：

$\mathrm{VSP}=\frac{\mathrm{Av}+{\mathrm{Bv}}^2+{\mathrm{Cv}}^3+\mathrm{mva}}{f}$

其中，v为机动车速度，单位m/s；a为加速度，单位m/s²；A、B、C、m、f为常数。

为了可以更好地反映VSP分布与平均速度之间的关系，以1kW/t为间隔对VSP进行划分并聚类，得到VSP bin，具体计算见下式：

$\∀:\mathrm{VSP}\∈[n-0.5,n+0.5),\mathrm{VSPbin}=n$

VSP分布数据库建立模块，以大量不同车型逐秒的行驶工况数据为基础，经过上述处理过程，最后对所有车型、道路类型、速度区间相同的VSP进行聚类，统计每一个速度区间对应的VSP个数占该类VSP总数的百分比，即VSP分布，如下式所示：

$R_{i,k}=\frac{N_{i,k}}{N_k}$

其中，R_i,k为第k个平均速度区间第i个VSP bin的分布率；N_k为第k个平均速度区间的VSP总数；N_i,k为第k个平均速度区间VSP bin为i的VSP个数。图3为轻型车在快速路VSP分布示例。

机动车排放数据采集单元对机动车排放数据进行采集、分析和处理，最终生成排放率数据库。进一步包括：

采集排放数据模块，利用PEMS车载尾气测试设备采集机动车油耗和尾气排放数据，数据主要字段包括采集日期、时间、速度、NO_X、HC、CO逐秒排放量、逐秒油耗量等。

燃料类型划分模块，燃料类型包括：汽油、柴油、液化石油气(Liquefied Petroleum Gas，LPG)、液化天然气(Liquefied Natural Gas，LNG)、压缩天然气(Compressed Natural Gas，CNG)、混合动力、双燃料。

排放标准划分模块，为有效的控制机动车污染和排放，按照车辆所执行的排放标准对车辆进行分类，根据测试车辆的生产年份，确定车辆的排放标准类型，划分车辆排放标准。

行驶里程划分模块，对同类车型的行驶里程进行划分。

平均排放率计算模块，计算每种车型、燃料类型、排放标准、行驶里程相应分类的每个VSP bin的平均排放率，如下式所示：

${\mathrm{ER}}_i=\frac{1}{m}\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{er}}_j$

其中，ER_i为第i个VSP bin的平均排放率，单位g/s；er_j是VSP bin为i的第j个VSP对应的排放率，单位为g/s；m是VSP bin为i的VSP个数。

排放率数据库建立模块，将不同车型、排放标准、行驶里程的上述计算结果整理、分类，建立各VSP bin排放率数据库。

数据耦合单元对交通运行数据和机动车排放数据进行耦合。进一步包括：排放因子计算模块，各平均速度区间排放因子的计算方法如下式所示：

${\mathrm{EF}}_k=(\underset{i}{\mathrm{\Σ}}{\mathrm{ER}}_i\×\mathrm{VSP}{\mathrm{bin}}_i)/v\×3600$

其中，EF_k为第k平均速度区间的排放因子，单位为g/km；ER_i是第i个VSP bin的平均排放率，单位为g/s；VSP bin_i是第k平均速度区间第i个VSP bin的分布值；v是第k平均速度区间的中值，单位为km/h。

通过VSP分布与VSP bin平均排放率的映射关系，可以得到不同车型、道路类型、燃料类型、排放标准、行驶里程、平均速度区间的排放因子。

油耗因子计算模块，计算方法如下式所示：

${\mathrm{EF}}_{\mathrm{Fuel}}=\frac{1}{\%C}(\frac{12}{44}{\mathrm{EF}}_{{\mathrm{CO}}_2}+\frac{12}{28}{\mathrm{EF}}_{\mathrm{CO}}+\frac{12}{13}{\mathrm{EF}}_{\mathrm{HC}})$

其中，EF_Fuel为油耗因子；EF_CO、EF_HC分别为CO₂、CO、HC排放因子；％C为燃料中C的质量比例。

对于实测数据无法覆盖到的分类情况，可以利用排放预测模型得到其相应的油耗因子。整理上述计算结果，建立车型、道路类型、平均速度、行驶里程、燃料类型、排放标准机动车数据库CO2、CO、HC、NO_X油耗因子。

路段排放量计算单元计算得到路段排放量。路段排放量计算方法如下式所示：

${\mathrm{ER}}_{\sec \mathrm{tion}}=l\·\underset{i,j,k,m}{\mathrm{\Σ}}{\mathrm{EF}}_{\mathrm{ijkm}}^v\·q({\mathrm{Vehicle}}_i,{\mathrm{Fuel}}_j,{S\tan \mathrm{dard}}_k,{\mathrm{Mileage}}_m)$

其中，ER_section为测算路段排放率，l为测算路段长度，q(Vehicle_i,Fuel_j,Standard_k,Mileage_m)为第i种车型、第j种燃料类型、第k种排放标准、第m类行驶里程的交通量，为上述分类对应的速度区间为v的机动车排放因子。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (10)

1.一种基于VSP分布的交通运行数据与排放数据耦合的方法，

具体包括：

第一步，采集交通运行数据；

第二步，采集机动车排放数据；

第三步，对交通运行数据和排放数据进行耦合；

第四步，进行路段排放量估算。

2.如权利要求1所述的方法，其中第一步进一步包括：

步骤1.1，工况数据采集，利用车载GPS设备，采集机动车逐秒行驶工况数据，包括采集日期、时间、经度、纬度、速度、方向角、高程等，将数据按照统一的格式进行整理，通过初步质量控制，标记不符合要求的记录，建立机动车行驶工况数据库；

步骤1.2，进行车型划分，将机动车按车重进行分类，分为微型车、轻型车、中型车和重型车；

步骤1.3，进行道路类型匹配，道路类型分为快速路、主干路、次干路、支路四种，道路类型的匹配工作在GIS系统中完成，每一条工况数据都包含经纬度字段，可以在已有的路网上创建点，通过判断该点所处的路段名得到其道路类型；

步骤1.4，进行短行程划分，对于每一种车型与道路类型，所有数据被划分到时长为180秒的短行程，每一个短行程都需要计算其平均速度，以便进行下一步的速度聚类；

步骤1.5，进行平均速度聚类，将上述过程得到的所有短行程依据平均速度进行聚类，得到其所属的速度区间，方法如下式所示：

步骤1.6，进行VSP聚类，VSP由瞬时速度和加速度计算得到，VSP的计算如下式所示：

$\mathrm{VSP}=\frac{\mathrm{Av}+{\mathrm{Bv}}^2+{\mathrm{Cv}}^2+\mathrm{mva}}{f}$

其中，v为机动车速度，单位m/s；a为加速度，单位m/s²；A、B、C、m、f为常数；

为了可以更好地反映VSP分布与平均速度之间的关系，以1kW/t为间隔对VSP进行划分并聚类，得到VSP bin，具体计算见下式：

$\∀:\mathrm{VSP}\∈[n-0.5,n+0.5),\mathrm{VSPbin}=n$

步骤1.7，建立VSP分布数据库，以大量不同车型逐秒的行驶工况数据为基础，经过上述处理过程，最后对所有车型、道路类型、速度区间相同的VSP进行聚类，统计每一个速度区间对应的VSP个数占该类VSP总数的百分比，即VSP分布，如下式所示：

$R_{i,k}=\frac{N_{i,k}}{N_k}$

其中，R_i,k为第k个平均速度区间第i个VSP bin的分布率；N_k为第k个平均速度区间的VSP总数；N_i,k为第k个平均速度区间VSPbin为i的VSP个数。

3.如权利要求1所述的方法，其中第二步进一步包括：

步骤2.1，采集排放测试数据，利用PEMS车载尾气测试设备采集机动车油耗和尾气排放数据，包括采集日期、时间、速度、CO₂、NO_X、HC、CO逐秒排放量、逐秒油耗量等；

步骤2.2，划分燃料类型，燃料类型包括：汽油、柴油、液化石油气(Liquefied Petroleum Gas，LPG)、液化天然气(Liquefied NaturalGas，LNG)、压缩天然气(Compressed Natural Gas，CNG)、混合动力、双燃料；

步骤2.3，划分排放标准，按照车辆所执行的排放标准对车辆进行分类，根据测试车辆的生产年份，确定车辆的排放标准类型；

步骤2.4，划分行驶里程，对同类车型的行驶里程进行划分；

步骤2.5，计算平均排放率，计算每种车型、燃料类型、排放标准、行驶里程相应分类的每个VSP bin的平均排放率，如下式所示：

${\mathrm{ER}}_i=\frac{1}{m}\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{er}}_j$

其中，ER_i为第i个VSP bin的平均排放率，单位g/s；er_j是VSPbin为i的第j个VSP对应的排放率，单位为g/s；m是VSP bin为i的VSP个数；

步骤2.6，建立排放率数据库，将不同车型、排放标准、行驶里程的上述计算结果整理、分类，建立各VSP bin排放率数据库。

4.如权利要求1所述的方法，其中第三步进一步包括：

步骤3.1，计算排放因子，各平均速度区间排放因子的计算方法如下式所示：

${\mathrm{EF}}_k=(\underset{i}{\mathrm{\Σ}}{\mathrm{ER}}_i\×\mathrm{VSP}{\mathrm{bin}}_i)/v\×3600$

其中，EF_k为第k平均速度区间的排放因子，单位为g/km；ER_i是第i个VSP bin的平均排放率，单位为g/s；VSP bin_i是第k平均速度区间第i个VSP bin的分布值；v是第k个速度区间的均值，单位为km/h；

步骤3.2，计算油耗因子，油耗因子定义为机动车行驶单位距离消耗的燃油质量，计算方法如下式所示：

${\mathrm{EF}}_{\mathrm{Fuel}}=\frac{1}{\%C}(\frac{12}{44}{\mathrm{EF}}_{{\mathrm{CO}}_2}+\frac{12}{28}{\mathrm{EF}}_{\mathrm{CO}}+\frac{12}{13}{\mathrm{EF}}_{\mathrm{HC}})$

其中，EF_Fuel为油耗因子；EF_CO、EF_HC分别为CO₂、CO、HC排放因子；％C为燃料中C的质量比例。

5.如权利要求1所述的方法，其中第四步路段排放量计算方法如下式所示：

${\mathrm{ER}}_{\sec \mathrm{tion}}=l\·\underset{i,j,k,m}{\mathrm{\Σ}}{\mathrm{EF}}_{\mathrm{ijkm}}^v\·q({\mathrm{Vehicle}}_i,{\mathrm{Fuel}}_j,{S\tan \mathrm{dard}}_k,{\mathrm{Mileage}}_m)$

其中，ER_section为测算路段排放率，l为测算路段长度，q(Vehicle_i,Fuel_j,Standard_k,Mileage_m)为第i种车型、第j种燃料类型、第k种排放标准、第m类行驶里程的交通量，为上述分类对应的速度区间为v的机动车排放因子。

6.一种基于VSP分布的交通运行数据与排放数据耦合的系统，包括交通运行数据采集单元，用于采集交通运行数据；机动车排放数据采集单元，用于采集机动车排放数据；数据耦合单元，用于对交通运行数据和排放数据进行耦合；路段排放量计算单元，进行路段排放量估算。

7.如权利要求6所述的系统，其中交通运行数据处理单元进一步包括：

工况数据采集模块，利用车载GPS设备，采集机动车逐秒行驶工况数据，包括采集日期、时间、经度、纬度、速度、方向角、高程等，将数据按照统一的格式进行整理，通过初步质量控制，标记不符合要求的记录，建立机动车行驶工况数据库；

车型划分模块，将机动车按车重进行分类，分为微型车、轻型车、中型车和重型车；

道路类型匹配模块，道路类型分为快速路、主干路、次干路、支路四种，道路类型的匹配工作在GIS系统中完成，每一条工况数据都包含经纬度字段，可以在已有的路网上创建点，通过判断该点所处的路段名得到其道路类型；

短行程划分模块，对于每一种车型与道路类型，所有数据被划分到时长为180秒的短行程，每一个短行程都需要计算其平均速度；

平均速度聚类模块，将上述过程得到的所有短行程依据平均速度进行聚类，得到其所属的速度区间，方法如下式所示：

VSP聚类模块，VSP由瞬时速度和加速度计算得到，VSP的计算如下式所示：

$\mathrm{VSP}=\frac{\mathrm{Av}+{\mathrm{Bv}}^2+{\mathrm{Cv}}^2+\mathrm{mva}}{f}$

其中，v为机动车速度，单位m/s；a为加速度，单位m/s²；A、B、C、m、f为常数。

为了可以更好地反映VSP分布与平均速度之间的关系，以1kW/t为间隔对VSP进行划分并聚类，得到VSP bin，具体计算见下式：

$\∀:\mathrm{VSP}\∈[n-0.5,n+0.5),\mathrm{VSPbin}=n$

VSP分布数据库建立模块，以大量不同车型逐秒的行驶工况数据为基础，经过上述处理过程，最后对所有车型、道路类型、速度区间相同的VSP进行聚类，统计每一个速度区间对应的VSP个数占该类VSP总数的百分比，即VSP分布，如下式所示：

$R_{i,k}=\frac{N_{i,k}}{N_k}$

其中，R_i,k为第k个平均速度区间第i个VSP bin的分布率；N_k为第k个平均速度区间的VSP总数；N_i,k为第k个平均速度区间VSPbin为i的VSP个数。

8.如权利要求6所述的系统，其中机动车排放数据处理单元进一步包括：

排放测试数据采集模块，利用PEMS车载尾气测试设备采集机动车油耗和尾气排放数据，包括采集日期、时间、速度、CO₂、NO_X、HC、CO逐秒排放量、逐秒油耗量等；

燃料类型划分模块，燃料类型包括：汽油、柴油、液化石油气(Liquefied Petroleum Gas，LPG)、液化天然气(Liquefied NaturalGas，LNG)、压缩天然气(Compressed Natural Gas，CNG)、混合动力、双燃料；

排放标准划分模块，按照车辆所执行的排放标准对车辆进行分类，根据测试车辆的生产年份，确定车辆的排放标准类型；

行驶里程划分模块，对同类车型的行驶里程进行划分；

平均排放率计算模块，计算每种车型、燃料类型、排放标准、行驶里程相应分类的每个VSP bin的平均排放率，如下式所示：

${\mathrm{ER}}_i=\frac{1}{m}\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{er}}_j$

其中，ER_i为第i个VSP bin的平均排放率，单位g/s；er_j是VSPbin为i的第j个VSP对应的排放率，单位为g/s；m是VSP bin为i的VSP个数；

排放率数据库建立模块，将不同车型、排放标准、行驶里程的上述计算结果整理、分类，建立各VSP bin排放率数据库。

9.如权利要求6所述的系统，其中数据耦合单元进一步包括：

排放因子计算模块，各平均速度区间排放因子的计算方法如下式所示：

${\mathrm{EF}}_k=(\underset{i}{\mathrm{\Σ}}{\mathrm{ER}}_i\×\mathrm{VSP}{\mathrm{bin}}_i)/v\×3600$

其中，EF_k为第k平均速度区间的排放因子，单位为g/km；ER_i是第i个VSP bin的平均排放率，单位为g/s；VSP bin_i是第k平均速度区间第i个VSP bin的分布值；v是第k平均速度区间的中值，单位为km/h。

油耗因子计算模块，计算方法如下式所示：

${\mathrm{EF}}_{\mathrm{Fuel}}=\frac{1}{\%C}(\frac{12}{44}{\mathrm{EF}}_{{\mathrm{CO}}_2}+\frac{12}{28}{\mathrm{EF}}_{\mathrm{CO}}+\frac{12}{13}{\mathrm{EF}}_{\mathrm{HC}})$

其中，EF_Fuel为油耗因子；EF_CO、EF_HC分别为CO₂、CO、HC排放因子；％C为燃料中C的质量比例。

10.如权利要求6所述的系统，其中路段排放量计算单元计算路段排放量的方法如下式所示：

${\mathrm{ER}}_{\sec \mathrm{tion}}=l\·\underset{i,j,k,m}{\mathrm{\Σ}}{\mathrm{EF}}_{\mathrm{ijkm}}^v\·q({\mathrm{Vehicle}}_i,{\mathrm{Fuel}}_j,{S\tan \mathrm{dard}}_k,{\mathrm{Mileage}}_m)$

其中，ER_section为测算路段排放率，l为测算路段长度，q(Vehicle_i,Fuel_j,Standard_k,Mileage_m)为第i种车型、第j种燃料类型、第k种排放标准、第m类行驶里程的交通量，为上述分类对应的速度区间为v的机动车排放因子。