1.城市峡谷街道机动车污染物排放与扩散的实时估计方法，其特征在于：城市道路交通网络建模部分A，机动车污染物排放部分和扩散部分B、C，以及污染物浓度显示部分D；
交通网路建模部分A，基于CTM模型对城区路网进行建模、仿真，将道路划分成不同链路，每条链路划分成不同元胞，基于信号灯配时建立元胞之间传输关系，估计计算建模街道的交通状况；
城市道路机动车污染物排放估计部分B，利用基于智能手机的交通数据实时采集系统，实时获取机动车实际工况下的行驶参数，结合VSP模型实时估计机动车污染物排放因子；
城市道路机动车污染物扩散估计部分C，针对城市道路的峡谷特性，选取OSPM模型，模拟道路机动车污染物的扩散过程，计算街道处不同污染物的时空浓度；
污染物浓度显示部分D，主要包括地图界面部分和选项栏部分，将污染物浓度分为6个等级，每个等级分配一个颜色，利用模型计算的数据，将污染物等级显示在地图上。
2.根据权利要求1所述的城市峡谷街道机动车污染物排放与扩散的实时估计方法，其特征在于：基于CTM模型对城区路网进行元胞划分的方法；划分步骤包括：
(1)将目标区域以交叉口和区域边界为节点划分成n条链路，每条链路沿车辆行驶方向划分成两部分，称之为方向链路；
(2)每条方向链路按l米长度划分成多个元胞，将方向链路与交叉口连接的元胞分为分离元胞和融合元胞，分离元胞组成交叉口处不同行驶方向的车道，融合元胞接收来自交叉口的车辆；
(3)对于丁字路口，分离元胞分为两种：直行元胞和左转元胞，或直行元胞和右转元胞，或左转元胞和右转元胞；对于十字路口，分离元胞分为三种：直行元胞、左转元胞和右转元胞；
(4)若交叉口处分离车道的长度为ld米，相同行驶方向的分离元胞个数则为round(ld/l)个，其中round(·)为四舍五入函数。
3.根据权利要求1所述的城市峡谷街道机动车污染物排放与扩散的实时估计方法，其特征在于：对已划分元胞和元胞之间截面流量的编号方法；编号方法如下：
首先对不同方向链路进行编号，具体编号方法为：链路编号+方向缩写，链路编号采用罗马数字；其次对方向链路内的元胞进行编号；
元胞编号方法：首先，按照直行、左转、右转将所有分离元胞分别编号为S、L、R，其次，自融合元胞或边界元胞起，沿车辆行驶方向除分离元胞外依次编为正整数1、2、…、n，分离元胞则在此顺序下继续编号，所以分离元胞的最终编号为S[n+1],S[n+2]；最后结合方向链路编号和元胞编号即可区分区域内所有的元胞；最后对截面流量进行编号；
截面编号方法：在元胞编号基础上编号，首先将每个元胞的输出流量编号为y元胞编号，其余的流量根据指向的元胞编号记为yU+元胞编号。
4.根据权利要求1所述的城市峡谷街道机动车污染物排放与扩散的实时估计方法，其特征在于：为得到路段内的车辆数和路段截面的车流量，对并列的元胞内车辆数和并列截面的总车流量的相加方法：
元胞内车辆数相加原则为：若链路内有n列元胞，相加的多个元胞的编号中的数字相加为n+1；
截面车流量相加原则为：若链路内有n列元胞，除链路两端的截面流量外，相加的两个截面流量的编号中的数字相加为n，且无下标“U”，两端的截面流量，相加的多个截面流量的编号中的数字相加为n+1。
5.根据权利要求1所述的城市峡谷街道机动车污染物排放与扩散的实时估计方法，其特征在于：利用基于智能手机的交通数据采集系统和VSP模型对实际行驶工况下机动车尾气排放因子进行估计：
通过基于智能手机的实时交通数据采集系统采集机动车不同行驶状况下的动态交通数据，包括实时的经纬度、瞬时速度﹑瞬时加速度以及行驶轨迹信息；
交通数据采集系统基于android智能手机操作系统编写API应用程序，实时调用GPS模块和加速度传感器模块数据并记录；应用程序以1Hz的固定采样频率实时获取数据信息，包括系统时间、GPS模块以及加速度模块信息，并以数据库文件形式保存在手机外部存储器内，由GPS模块能够获得位置信息和速度信息，由加速度传感器能够获得加速度信息；
对采集的交通动态数据进行预处理方法，得到瞬时速度与瞬时加速度的秒数据，结合VSP模型逐秒计算可获得各秒的排放水平，逐次查询获得对应的排放量清单，累加后可得到数据对应的出行过程的总排放量。
6.根据权利要求1所述的城市峡谷街道机动车污染物排放与扩散的实时估计方法，其特征在于：污染物浓度显示部分以OSPM模型计算所得结果为基础，为将模型的实时估计结果直观可视化，本方法将污染物浓度分为6个等级，并分配不同的颜色，根据OSPM模型计算得出的污染物浓度值，选取代表路段污染物浓度的颜色，将其实时显示在对应街道处，用户通过查看路段处的颜色来查询街道处污染物的浓度值。
7.根据权利要求1所述的城市峡谷街道机动车污染物排放与扩散的实时估计方法，其特征在于：实现城市峡谷街道机动车污染物实时估计的整个实施过程；
第一部分是交通建模与仿真；本方法首先基于CTM模型对城区路网进行建模、仿真；
根据CTM模型将空间和时间离散化；空间离散化即元胞划分方法，将街道路段划分成多个具有相同长度的元胞，同时将时间划分为等时间间隔，作为时间步长，且满足在一个时间步长内，车辆以最大速度行驶的距离小于元胞长度；本方法的元胞划分方法为：选取元胞长度为l米，时间间隔为t秒，l取20，t取1；具体划分步骤如下：
(1)将目标区域以交叉口和区域边界为节点划分成n条链路，每条链路沿车辆行驶方向划分成两部分，称之为方向链路，代表同一条道路不同的行驶方向；
(2)每条方向链路按l米长度划分成多个元胞，将方向链路与交叉口连接的元胞分为分离元胞和融合元胞，分离元胞组成交叉口处不同行驶方向的车道，分离元胞内的车辆将驶入交叉口，融合元胞接收来自交叉口的车辆；
(3)分离元胞根据交叉口类型，元胞的划分方法不同，对于丁字路口，分离元胞分为两种：直行元胞和左转元胞，或直行元胞和右转元胞，或左转元胞和右转元胞；对于十字路口，分离元胞分为三种：直行元胞、左转元胞和右转元胞；
(4)若交叉口处分离车道的长度为ld米，相同行驶方向的分离元胞个数则为i＝round(ld/l)个，其中round(·)为四舍五入函数；
按以上元胞划分步骤，将路网划分为元胞，并在此基础上建立交通网络拓扑图；
为区分不同元胞、不同截面流量以及方便编程，需要对交通网络拓扑图中的元胞和截面流量进行编号；本方法以方向链路为单位进行编号，不同方向链路内元胞和截面流量的编号方法相同；
首先对不同方向链路进行编号，具体编号方法为：链路编号+方向缩写，链路编号采用罗马数字，方向缩写采用N、S、W、E，分别代表北、南、西、东；其次对方向链路内的元胞进行编号；
元胞编号方法：首先，按照直行、左转、右转将所有分离元胞分别编号为S、L、R，其次，自融合元胞或边界元胞起，沿车辆行驶方向除分离元胞外依次编为正整数1、2、…、n，分离元胞则在此顺序下继续编号，所以分离元胞的最终编号为S[n+1],S[n+2],…,S[n+i]，其余方向以此类推；最后结合方向链路编号和元胞编号即可区分区域内所有的元胞；最后对截面流量进行编号；
截面流量编号方法：在元胞编号基础上编号，首先将每个元胞的输出流量记为y元胞编号，其余的流量根据指向的元胞编号记为yU+元胞编号；
根据每个元胞对应的道路几何信息和通行能力，确定以下元胞参数：
R：一个时间步长T内能流入元胞的最大车辆数；
S：一个时间步长T内能流出元胞的最大车辆数；
Q：一个时间步长T内元胞之间截面能流过的最大车辆数；
N：元胞中能容纳的最大车辆数；
Nc：路段为临界密度时元胞中的车辆数；
建立元胞之间交通流传输关系；元胞之间的连接包括三种方式：简单连接﹑融合和分离，各个连接的计算方法如下：
(1)简单连接
交通拓扑图中，两元胞之间只有单个箭头的连接即为简单连接；设元胞1为上游元胞，元胞2为下游元胞，简单连接的元胞1和元胞2之间实际的传输流量为：
其中：y12为T时间内元胞1实际的流出车辆数，也为T时间内元胞2实际的流入车辆数；n1和n2分别为元胞1和元胞2内现有的车辆数； 分别为T时间内元胞1能流出的和元胞2能流入的最大车辆数；N2元胞2内能容纳的最大车辆数；δ2为元胞2的自由流速度和反向波速度的比，即δ＝Nc/(N-Nc)；
(2)融合连接
交通拓扑图中，由多个箭头汇集到一个箭头的连接即为融合连接；设上游元胞为元胞S和元胞L，它们同时流入下游的元胞1，元胞S、元胞R和元胞1之间的实际传输流量公式为：
YU1≥YS+YR时：
YU1＜YS+YR时：
yU1＝yS+yR
其中
式中：αS和αR为融合比，即元胞1分配给元胞S和元胞R的空间比例，而且αS+αR＝1；mid函数表示取三个值的中间值；YS和YR表示T时间内元胞S和元胞R能流出的最大车辆数；YU1表示T时间内元胞1能够接收的最大车辆数；yS和yR表示T时间内元胞S和元胞R实际流出的车辆数；yU1表示T时间内元胞1实际接收的车辆数；
融合连接因不同交叉路口和不同红绿灯相位实际为以下几种情况：无连接、简单连接和融合连接；实际为融合连接的连接用实际传输流量公式计算；仿真时，根据信号灯配时选择元胞之间不同的连接方式；
对于三相位信号灯的丁字路口，有3个融合元胞，没有接收来自右转分离元胞内车辆的融合元胞，与其他两个元胞做融合连接时，实际为无连接或者简单连接：当直行信号灯为绿时，该融合元胞只接收到来自对应直行元胞的车辆，实际只与一个元胞连接，所以此时为简单连接；当无直行的方向左转信号灯为绿时，该融合元胞只接收到来自左转元胞的车辆，亦然；当有直行的方向左转信号灯为绿时，此时无直行的方向左转信号灯为红，没有车辆流入该融合元胞，此时为无连接；其余两个融合元胞情况一致，实际为简单连接或融合连接；
对于十字路口，有4个融合元胞，它们的情况一致，实际都为简单连接或融合连接，且最多为两个元胞融合到一个元胞；
(3)分离连接
元胞之间由一个箭头分离为多个箭头的连接即为分离连接；对于丁字路口，只有两个行车方向选择，所以路口前一个元胞分离为两个元胞，对于十字路口，由三个行车方向选择，所以路口前一个元胞分离为三个元胞；两种情况计算方法类似，具体说明十字路口，丁字路口只需去掉一个元胞的相关参数即可；设上游元胞为元胞15，下游元胞为元胞L、元胞S和元胞R，所以元胞15与元胞L、元胞S、元胞R之间的实际传输流量为：
yUL＝βLy15；yUS＝βSy15；yUR＝βRy15
其中βL、βS和βR为分离比，即车流从元胞15进入元胞L、元胞S和元胞R的比例，且βL+βS+βR＝1；
利用以上元胞连接方法，并分别为这两个路口配置交通信号灯配时方案，用实时的交通流量作为路段的边界值进行仿真，即可得到该路段不同元胞内的车辆数、元胞之间的传输流量交通数据；
为得到路段内的车辆数和路段截面的车流量，即链路内两个方向链路中并列元胞的总车辆数和并列截面的总车流量，还需要对并列的元胞内车辆数和并列截面的总车流量进行相加，为满足并列条件，仿真时元胞内车辆数相加原则为：若链路内有n列元胞，相加的多个元胞的编号中的数字相加为n+1；截面车流量相加原则为：若链路内有n列元胞，除链路两端的截面流量外，相加的两个截面流量的编号中的数字相加为n，且无下标“U”，两端的截面流量，相加的多个截面流量的编号中的数字相加为n+1；
至此，路网内每个路段的车辆数和每个截面的车流量即可仿真得到；
第二部分是对机动车污染物排放量的计算；本方法计算机动车主要的排放物，包括氮化物NOX、碳氢化物HC和一氧化碳CO；机动车污染物排放量可由污染物排放源强计算得出，每种污染物的排放强源可用下式计算：
式中，Q为强源[mg·(m·s)-1]；Ni为某类车的车流量(辆·h-1)；Ei为某类车的排放因子[g·(km·辆)-1]；n为机动车分类数；
由排放强源表达式能够看出，机动车污染物排放量与交通流量和机动车污染物排放因子有着直接的关系；交通流量是指在单位时间内，通过道路某横断面的车辆数，单位是辆·-1
h ；机动车排放因子是指单辆机动车运行单位里程或消耗单位燃料排放污染物的量，单位是g/km，它是反映机动车排放状况的最基本的参数，也是确定机动车污染物排放总量及其环境影响的重要依据；本方法中，交通流量由第一部分的CTM模型仿真得出，机动车污染物排放因子根据机动车实际工况下的行驶参数，结合排放因子模型实时模拟估计得出；
针对机动车排放因子在实际应用过程中存在的不同车辆技术类型、不同工况差异性问题，本方法利用基于智能手机的交通数据采集系统实时采集城区街道机动车实际道路行驶参数，基于VSP模型对实际行驶工况下机动车尾气排放因子进行估计，结合第一部分中计算得到的交通流量确定交通源的排放因子；
通过基于智能手机的交通数据实时采集系统采集机动车不同行驶状况下的动态交通数据，获得车辆实时的经纬度、瞬时速度﹑瞬时加速度以及行驶轨迹信息；对采集的机动车行驶数据进行分析，估计机动车实际工况下污染物排放因子；
交通动态采集系统基于android智能手机操作系统编写API应用程序，实时调用GPS模块和加速度传感器模块数据并记录；应用程序以1Hz的固定采样频率实时获取数据信息，包括系统时间、GPS模块以及加速度模块信息，并以数据库文件形式保存在手机外部存储器内，供后续分析使用；由GPS模块能够获得位置信息和速度信息，由加速度传感器能够获得加速度信息，包括三个方向上的加速度，垂直屏幕向上的方向为Z轴，屏幕由左至右方向为X轴，由下至上方向为Y轴；
采集的交通动态数据经过数据预处理方法后，得到瞬时速度与瞬时加速度的秒数据；
根据处理后的速度与加速度信息秒数据，结合VSP模型逐秒计算可获得各秒的排放水平，依据所得排放水平逐次查询表获得对应的排放量清单，由于清单对应的即为当前一秒，所以排放清单即这一秒的排放量，累加后可得到数据对应的出行过程的总排放量；VSP定义为单位质量机动车的瞬时功率，是发送机克服滚动阻力和空气阻力，增加机动车的动能时能所输出的功率；VSP计算公式如下：
其中：KE为机动车动能；PE为机动车势能；Frolling为车轮的阻力；FAerodynamic为空气阻力；
v为机动车的运行速度；m为机动车的总质量；ξi为滚动质量系数；g为重力加速度；h为机动车距海平面的距离；CR为轮胎转动阻力系数；ρα为空气密度；CD为风阻系数；A为机动车前沿面积；vm为机动车顶风风速；a为机动车加速度；grade为道路坡度；
本方法中假定路段内机动车的车型为：轻型车和重型车，并根据交通时段和GPS信息分配不同的比例；根据已标定的特征选取相关参数简化得到相关类型车辆VSP计算公式；其中，轻型车取值如下，ξi取0.1，CR取0.0135，CD×A/m取0.0005，vm和grade忽略不计：
3
VSP＝v×(1.1×a+0.132)+0.000302×v
通过对上式相关参数进行标定，得到重型车的VSP表达式为：
VSP＝v×(a+0.09199)+0.000169×v3
估算出行距离；GPS接收模块获得的经纬度位置信息是十进制小数形式的数据，单位为度；因此只需要依次计算一段数据中相邻两个经纬度坐标间的距离，并累加即可得到整个出行过程的行驶距离；假设地球是一个完美的球体，即半径为平均半径，记做R；设一点A的经纬度坐标为(LngA,LatA)，另一点B的经纬度坐标为(LngB,LatB)；以0度经线和赤道为基准，东经经度取正值，西经经度取负值，北纬取负值并加上90度，南纬取正值加上90度，由此能够得到转换后的A、B两点的坐标，分别为(MLngA,MLatA)和(MLngB,MLatB)，则根据三角推导得出计算两点距离Distance的如下公式所示：
C＝sin(MLatA)×sin(MLatB)×cos(MLngA-MLngB)+cos(MLatA)×cos(MLatB)
Distance(A,B)＝R×Arccos(C)×π/180
其中，Distance与R的单位保持一致；
基于VSP排放模型，并结合时间段内所行驶的距离即可估计出不同行驶工况下污染物排放因子
其中： 的单位为g/km，Distance(A,B)为通过经纬度信息所计算得到的两点间距离；
至此，即可计算出不同类型车辆的排放因子，再结合车辆类型查表获得各类型车辆每种污染物的排放因子，将其带入源强计算公式，结合对应类型车辆的车流量，即可计算出不同类型车辆各种污染物的排放量，公式将不同类型车辆排放的同种污染物相加，最终即得到路段不同污染物的总排放量；
第三部分基于OSPM扩散模型模拟计算道路机动车尾气污染物扩散浓度；城市街道两旁都有建筑物，考虑街道内的空气流动对污染物传输和扩散研究，这就不仅需要考虑排放源﹑气象条件和环境因素的影响，而且要考虑街道两旁建筑物形成的小尺度污染扩散过程，称为街道峡谷污染扩散；针对峡谷类型路段，本方法选取OSPM模型，应用模型模拟的方式计算道路机动车尾气污染物扩散浓度；
OSPM模型认为街道峡谷中受体点的污染物浓度由两个主要的贡献源组成，一个是由街道中交通排放导致的浓度部分，另一个是区域其他所有源贡献所致的背景浓度部分；前者又主要由交通排放直接引起的扩散浓度和街区内因气流涡旋引起污染物循环形成的浓度构成，并分别用高斯烟羽模型和箱模型来计算；于是街区内受体点污染物的总浓度为：
Ctot＝Cd+Cr+Cbk
式中：Cd是为直接扩散浓度，Cr为街道峡谷中循环浓度，Cbk为城市背景浓度；
利用高斯烟羽模型计算背风面直接贡献浓度Cd的公式如下：
d1＝min(Lmax,Lr)
d6＝min[max(Lmax,Lr),x1]
d7＝max(Lmax,x1)-x1
迎风面直接贡献浓度Cd的公式如下：
d4＝min[(Lmax-Lr),x1]
d5＝max[(Lmax-Lr),x1]-x1
利用箱模型计算循环区域内的污染物浓度Cr的公式如下：
d2＝min(Lmax,0.5Lr)
d3＝max[0,d1Ls/(0.5Lr-1)]
其中：Q为街道机动车排放源强，由第二部分获得；ub为街道处风速，由ub＝β·ut得出，式中ut为峡谷顶部的风速，风速转换系数β的值取0.34； 为垂直湍流速度
的标准方差，a为经验常数，根据风洞实测确定，取为0.1，σwo为交通引起的湍流，在模式中取为定值σwo＝0.1m/s；L为街道宽度；Lmax＝L/sin(θ)，θ是风向与街道的夹角；Lr＝2Hr为循环区下边长度，r是低风速修正因子，ub＞2m/s时r＝1.0，ub＜2m/s时r＝0.5；H为建筑物高度；
h0为汽车尾气初始混合高度，取h0＝2m；x1＝ub(H-h0)/σw为积分上限；R＝Max[0,cos(2rθ)]为权重因子；σwt为峡谷顶部湍流速度的方差；ud为对流速度；
气象参数采用中国气象局网站的常规监测数据，包括实时的风速、风向和背景浓度值；
街道特征通过实地考量获得，包括建筑物平均高度、街道宽度和建筑物位置；
将OSPM模型和相关参数输入计算机，结合第一部分得出的交通流数据和第二部分得出的排放源强数据，就能够实时仿真计算出路网内街道不同位置各种污染物的浓度；
第四部分为污染物浓度显示部分；本方法污染物浓度显示界面主要包括污染物浓度地图显示界面以及选项栏；
界面右上角为选项栏，本方法将污染物浓度分为6个等级，并分配不同的颜色，依次为绿、黄、橙、红、酒红、深红；对于不同的污染物，根据其浓度的数量级，每种颜色对应的浓度阶梯值不同，氮氧化物和碳氢化物的浓度阶梯值分为0～100、100～200、200～700、700～
1200、1200～2340、＞2340，单位为μg/m3，一氧化碳浓度阶梯值分为0～5、5～10、10～35、35～60、60～90、＞90，单位为mg/m3；
根据第三部分OSPM模型计算得出的污染物浓度值，以及以上污染物浓度等级划分方法，选取代表路段污染物浓度的颜色，将其实时显示在对应街道处，用户通过查看路段处的颜色来查询街道处污染物的浓度值；
显示界面只显示一种污染物的浓度等级，用户能够通过切换污染物选项栏查看不同污染物的浓度等级；
同时将风力、风向天气信息显示到选项栏供用户查看。