基于SVR的区域交通碳排放预测研究

0 引言

随着工业化发展速度加快和新兴经济体规模的扩大，我国已成为世界上最大的能源生产和消费国，同时也是世界上CO2排放量最大的国家.研究表明，交通运输业是对全国碳排放强度下降和碳排放增长贡献最大的行业之一，且现阶段能源消耗量和碳排放量均呈逐年上升趋势[1].因此，交通运输业应承担起降低碳排放的重要责任.

目前，国内外专家学者关于交通碳排放的研究主要包括两大类：其一，为交通碳排放影响因素分析，多数学者站在整个国家或某具体省市的角度，宏观分析交通碳排放驱动力作用[2-3]；或研究不同交通领域和不同运输结构诸多碳排放影响因素[4]，同时对某一具体指标(如碳排放效率、能源利用率等)变化影响机理的研究也日益增加[5-6].但是我国区域经济发展不均衡，不同省市区域之间的交通碳排放水平存在明显差异，因此单纯从国家或某具体区域角度分析交通碳排放影响因素问题会存在一定局限性.其二，为交通碳排放预测及节能减排策略研究，大部分学者通过构建Kaya恒等式、IPAT或STIRPAT模型及其扩展模型、GM(1，1)模型等来直接预测碳排放总量[7-9]；或是结合线性回归法和情景预测法预测不同情景下的碳排放情况并提出针对性建议[10].但是，传统预测方法存在偏最小二乘回归和岭回归稳定性低、可解释性较弱、关键参数难确定等不足.此外，区域交通碳排放系统为复杂非线性系统，传统预测方法因缺乏数据样本学习过程而导致非线性关系模糊，计算过程复杂，预测结果精度降低.

英、汉语言富具文化意象。这些意象有时通用，有时却只是英、汉语言各自独有。如何翻译这些富含语言、文化特色的意象，译者翻译的主观能动性至关重要。是尽量保留这些意象，还是转换意象，或者干脆舍弃意象，翻译的目的决定原文信息在译文中要不要体现，或者以何种方式体现。试看一例：

综上所述，本研究综合考虑我国区域发展特征和交通碳排放研究现状，借鉴STIRPAT模型中各影响因素的选取方法并对其进行扩展，选择具有区域特征的碳排放影响因素分析我国区域交通碳排放变化规律；同时，为避免传统预测方法的诸多弊端，构建支持向量回归机(Support Vector Regression,SVR)模型对我国区域交通碳排放进行预测.该模型可有效解决预测模型中小样本、非线性等问题，且在高维模式识别中表现出一定优势，因此在行业碳排放预测方面具有良好的实用性[11]，可为交通碳排放预测提供新的科学方法.

该次数据证实，实验组患者的 TC（5.68±1.26）mmol/L、TG（2.69±1.04）mmol/L、LDL-C（4.15±0.92）mmol/L 以及ApoB（1.19±0.15）g/L 指标均高于对照组，对照组的HDL-C（1.40±0.46）mmol/L 以及 ApoA（1.20±0.21）g/L 高于实验组，组间对比差异有统计学意义（P<0.05）。对照组血浆黏度（PV）、全血黏度、HCT（血细胞比容）以及 ESR（血沉）等血液流变学指标高于实验组，两组对比差异有统计学意义（P<0.05）。

1 影响因素选取和模型建立

1.1 碳排放影响因素选取

STIRPAT模型是由York等提出的关于人口、富裕和技术对环境的随机影响模型，为提高模型的分析解释能力，模型各因素均可进行扩展改进[3].交通碳排放主要影响因素包括交通自身影响、经济发展和科学技术3个方面[3]，因此选取旅客周转量、货物周转量、人均GDP、机动车保有量、碳排放强度、能源结构和城市化率能源效率7项影响因素作为交通碳排放预测自变量.其中，货物周转量和旅客周转量反映交通自身影响，人均GDP和机动车保有量反映社会富裕程度影响，碳排放强度和能源结构反映交通运输业科技发展水平影响，城市化率作为扩展项用来反映区域特征.

人工栽培的黄芩品质受栽培区域生态因子及田间管理影响较大，如山东莒县所产黄芩中黄芩苷含量显著高于甘肃产黄芩[4]。生产过程中，由于各地生态因子及田间管理差异，出现了施肥种类差异大、中耕除草难以规范以及农药不合理使用等现象[5]，部分产区连作加重了黄芩根腐病的发生几率[6]，农残重金属超标越发严重，以上因素单一或综合限制了黄芩药材的产量及品质，造成市售黄芩药材质量不一，影响药材药效，同时也影响临床用药安全和黄芩的市场竞争力。中药无公害栽培体系尚不完善[7]，现阶段开展黄芩无公害栽培技术研究具有迫切性。

碳排放强度是单位GDP所产生的碳排放量，用来反映经济增长与碳排放之间的关系.计算方法为

 

2016年北京市客、货物周转量分别为1 889.31亿人公里和671.33亿吨公里，样本数据年均增长率分别为11.2%和3.6%.预测2017—2020年北京市客、货运周转量分别为2 100.92、2 336.22、2 597.88、2 888.84 亿 人公里和 695.50、720.53、746.47、773.35亿吨公里.

[1]LIN B,XIE C.Reduction potential of CO2,emissions in China׳s transport industry[J].Renewable&Sustainable Energy Reviews,2014(33):689-700.

为了假设 H1，参考 Harford et al.（2008）的研究，本文构造了如下所示的模型（1）来检验独立董事的海外背景对企业现金持有水平的影响。

 

式中：R为不同能源消耗所占比例；Ei为第i种能源消耗量；ci为第i种能源转换标准煤系数[12].

1.2 预测模型建立

1.2.1 支持向量回归机

支持向量回归机是一种基于SVM思想，利用统计学习理论进行回归计算的机器学习法[13].该方法适用于有限样本研究，可在理论上得到全局最优解，且计算过程复杂程度与样本维数无关，在函数逼近、回归预测等方面可达到最优效果.本研究选择ε-SVR模型进行区域交通碳排放预测，具体步骤如下：

给定线性训练集为

 

在得到的最优解中，若，则称训练集中对应输入的xi为支持向量.若存在，则 由 计 算可得

4.2.3 在学校建设中，体育特色的传承十分重要。在学校办学理念变更时，如何传承原有特色是管理层要深入思考的问题。

设Rn上的一个线性函数表达式为

 

式中：ω为权值向量；b为偏移常量.

依据数理统计学理论，函数估计问题可转化为式(5)和式(6)的最优化问题.

 

式中：ε为不敏感损失系数.

Study on formation and spatial distribution of oyster reefs in paleo-coast belt, Tianjin WANG Ya-ming(108)

由于个别数据在ε精度下不能完成估计，故引入松弛变量ξi和，将优化过程式(5)转化为式(7).

 

式中：C为惩罚系数.

为便于求解，通常引入拉格朗日乘子α和α*，将模型转换为其对偶问题.

 

式中：xi为输入向量；y为输出值.

为了对混合动力汽车做进一步的阐述，笔者参考了国际能源组织(IEA)的有关文献对混合动力车辆的定义，能量与功率传递路线具有如下特点的车辆称为混合动力车辆。

 

对应的线性回归函数表达式为

 

引入核函数K(xi,yi)，将线性问题转化为Hilbert空间中的非线性回归问题，构建ε-SVR模型为

 

最终为所求最优解，对应函数表达式为

 

1.2.2 基于ε-SVR的交通碳排放预测模型

获取n年某区域交通碳排放及各项交通碳排放影响因素数据组成样本集，i=1,2,…,n，构建区域交通碳排放SVR模型.具体步骤如下：

(1)对自变量和因变量进行归一化处理，保证所有数据在[0,1]之间.

 

碳排放强度是衡量交通运输业科技发展水平的指标之一，由历史数据可知其年均降低比例为6.0%，预测2017—2020年交通碳排放强度分别为0.112 6、0.105 9、0.099 5和0.093 5 t/万元.

(3)对训练样本集数据进行仿真，得到模型最优解和回归函数f(x)，将训练样本集和测试样本集的全部数据分别代入函数输出拟合值，并对拟合结果与真实值进行线性回归，依据相关系数结果判断模型的学习推广能力.若模型学习推广能力差，需返回步骤(2)直到获得最优解.

(4)将预测年限的交通碳排放影响因素值依照式(16)进行归一化处理，代入回归函数f(x)后输出结果并对其依照式(17)进行反归一化处理，即求得预测年限的交通碳排放数据.

 

2 数据来源与实例分析

2.1 数据来源

以北京市为例进行实例分析，所研究的交通运输业为交通运输、仓储及邮政业3个行业的综合.所有数据均来源于《中国能源统计年鉴》《北京市统计年鉴》，样本区间为1990—2016年；所涉及的人均GDP、碳排放强度指标均为参照1990年可比价重新计算的结果.

 

表1 交通碳排放及其影响因素数据Table 1 Data of carbon emissions and influence factors of transportation

  

年份1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 CO2排放量/万t 302.395 1 330.359 8 358.434 2 388.131 9 334.719 1 390.165 8 449.685 8 466.910 5 504.404 2 591.139 4 655.182 8 776.177 6 833.709 5 884.682 2 1 115.462 3 1 132.729 9 1 495.601 5 1 839.959 3 2 232.665 1 2 360.193 6 2 510.692 3 2 670.853 1 2 902.455 9 2 599.092 7 2 724.994 3 2 802.442 0 2 982.217 5旅客周转量/亿人公里119.80 138.35 166.68 165.61 190.67 207.69 221.92 227.22 241.98 270.15 313.99 346.26 396.16 393.31 758.01 838.08 825.45 960.35 1 042.00 1 146.48 1 399.54 1 528.65 1 595.79 1 498.77 1 602.72 1 747.68 1 889.31货物周转量/亿吨公里268.82 279.97 293.87 314.35 310.77 323.10 317.62 312.91 284.67 283.89 299.64 315.98 340.55 362.08 402.28 457.74 423.12 449.04 454.22 441.23 513.66 616.93 638.31 686.91 672.82 623.69 671.33人均GDP/万元0.463 5 0.549 4 0.645 8 0.800 6 1.024 0 1.269 0 1.425 4 1.662 1 1.912 8 2.140 7 2.412 7 2.698 0 3.073 0 3.477 7 4.091 6 4.599 3 5.172 2 6.009 6 6.449 1 6.694 0 7.385 6 8.165 8 8.747 5 9.464 8 9.999 5 10.649 7 11.500 0机动车保有量/万辆29.9 33.2 39.0 48.3 56.6 71.7 79.6 108.2 130.5 139.8 157.8 169.9 189.9 212.4 229.6 258.3 287.6 312.8 350.4 401.9 480.9 498.3 520.0 543.7 559.1 561.9 571.8碳排放强度/（t/万元）0.603 8 0.551 6 0.505 5 0.438 0 0.292 3 0.258 8 0.251 3 0.224 8 0.212 2 0.220 7 0.207 2 0.209 3 0.193 2 0.176 7 0.184 9 0.162 5 0.184 2 0.186 9 0.200 9 0.194 2 0.177 9 0.164 3 0.162 3 0.131 3 0.127 7 0.121 8 0.119 8能源结构/%10.98 12.16 11.50 7.25 11.13 15.67 15.07 5.95 13.12 12.87 13.87 10.49 9.97 13.67 15.09 11.23 13.85 17.47 23.23 24.72 23.57 20.43 27.88 15.98 16.21 15.12 16.73城市化率/%73.48 73.86 74.32 74.73 75.20 75.63 76.06 76.48 76.89 77.29 77.54 78.06 78.56 79.05 79.53 83.62 84.32 84.49 84.92 85.00 85.93 86.18 86.21 86.29 86.31 86.46 86.50

2.2 MATLAB运行结果及回归分析

依照ε-SVR预测模型构建步骤，对27年样本数据进行归一化处理.随机选取13年样本数据组成训练样本集，剩余14年数据样本组成测试样本集.运用MATLAB 7.0处理相关数据，设置ε的初始值为0.01，C和γ的取值范围均是[2 - 8，28]，网格宽度为0.8，采用5折交叉验证法处理训练样本集数据.系统对C的取值范围进行离散化处理，取不同的离散值使得支持向量回归机的训练误差最小，图1为参数选择效果3D图.最终得到C∗为48.502 9，ε∗为0.006 801，交叉验证均方误差(MSE)为0.008 040.

[2]DARIDO G,TORRES-MONTOYA M,MEHNDIRATTA S.Urban transport and CO2,emissions:some evidence from Chinese cities[J].Wiley Interdisciplinary Reviews Energy&Environment,2014,3(2):122-155.

  

图1 参数选择效果图(3D)Fig.1 The effect drawing of the paramater selection(3D)

  

图2 训练样本集回归结果Fig.2 Linear regression of the training set

  

图3 测试样本集回归结果Fig.3 Linear regression of the testing set

  

图4 训练预测值与实际值对比结果Fig.4 Forecast and actual values comparison of training set

2.3 未来交通碳排放预测

2.3.1 影响因素预测分析

根据北京市统计局相关数据、北京“十三五”发展规划纲要、北京“十三五”交通发展规划、北京“十三五”能源发展规划、北京市政府统计公报及政府年度工作报告等相关数据，对北京市2017—2020年交通碳排放7项影响因素值进行整理和预测分析.

  

图5 测试预测值与实际值对比Fig.5 Forecast and actual values comparison of testing set

(1)客、货物周转量.

式中：I为碳排放强度；Ct为历年交通运输行业碳排放总量；Gt为历年经济生产总值.

大力推行纳税评估智能化，是完善纳税评估、提升税收管理质量效率、充分挖掘数据与信息资源的需要。实现纳税评估智能化关键是纳税评估信息与相关数据自动分析比对，实现方法主要是开发相关软件和应用程序。研发纳税评估智能化软件，无缝对接大数据库的同时，还应实现行业模型与软件的匹配结合，通过预先设定的税种勾稽关系、数据模型等，实现计算机自行比对计算。评估人员只要进行简单操作，如输入纳税人代码，软件即可自动运用征管系统中信息资料和数据。通过针对评估对象纳税行为进行数据化测算，对所评估对象作出更科学、准确的风险分析报告，为评估人员制定下一步行动方案提供参考。

(2)人均GDP.

课堂互动不仅包括师生之间的互动，学生通过小组合作探究方式展开的交流也是一种课堂互动形式。高中数学教师可以将小组合作探究教学方法与课堂教学内容有机结合。首先，教师可以将班级学生分成几个小组；其次，教师可以将教学内容转换成学习任务布置给学生，并引导学生在小组中展开讨论和交流，共同探究完成教师布置的学习任务；最后，学生讨论完毕，每个小组可以选出一名学生将小组讨论的结果汇报给教师。教师根据学生的讨论结果讲解课程内容。

北京“十三五”发展规划纲要明确提出，到2020年北京市人口总数预测控制在2 300万人之内，年均增长率为1.16%；GDP年均增速保持在6.5%，不低于6%.因此，综合考虑北京市政治重要性、经济发展水平现状及中心产业转移政策等因素，设置预测年限GDP增长率分别为6.0%、6.5%、7.0%，对应年份人均GDP如表2所示.

 

表2 不同GDP增速下的人均GDP预测值Table 2 Forecast per capital GDP under different GDP growth rate

  

2020年13.576 8 13.834 7 14.096 4 GDP增速/%6.0 6.5 7.0人均GDP预测值/万元2017年11.955 3 12.011 7 12.068 1 2018年12.473 0 12.591 0 12.709 5 2019年13.013 2 13.198 2 13.385 0

(3)机动车保有量.

北京交通“十三五”规划中明确提出，预期到2020年北京市机动车保有量控制在630万辆以内，最大年均增长率为2.5%.样本数据增长率最小值为0.5%，因此本文将机动车保有量增长率设置为0.5%、1.5%、2.5%，对应年份机动车保有量如表3所示.

 

表3 不同增速下的机动车保有量预测值Table 3 Forecast vehicle population under different growth rate

  

量增速/%2017年2018年2019年2020年0.5574.7577.5580.4583.3 1.5580.4589.1597.9606.9 2.5586.1600.7615.8631.2机动车保有机动车保有量预测值/万辆

(4)碳排放强度.

海洋上有好几只无形的“手”能随时拨弄海面上孤零零的冰山——海风、洋流能直接推动冰山转向；日晒和气温升高能让冰山融化，改变冰山的形状。

(2)在样本集中选取m个样本组成训练样本集，剩下n-m个样本组成测试样本集.选择K(xi,x)=exp(-γ‖‖xi-x2)为径向基核函数处理训练样本集，构建ε-SVR模型.设定ε的初始值，应用格搜索和交叉验证法确定参数C和γ的取值，最终得到最优参数ε∗、C∗和γ∗.

(5)能源比例.

北京市交通“十三五”规划明确指出，到2020年北京市交通运输业能耗增长率将控制在6.2%之内，而清洁能源消费比例将继续提升.由样本数据可知，2016年北京市交通运输业标准煤消耗量为1 312.69万t，清洁能源转换为标准煤的消耗量为219.56万t且年均增长率为10.82%.预测2017—2020年清洁能源比例分别为17.45%、18.25%、19.09%和19.97%.

(6)城市化率.

2016年北京城市化率为86.50%，年均增长率为0.63%，预测2017—2020年北京城市化率分别为87.04%、87.59%、88.15%和88.70%.

2.3.2 交通碳排放预测分析

由不同情境下的预测结果(表4)可知，同一GDP增速下，机动车保有量增长率越低，交通碳排放总量越少，说明努力控制机动车保有量的增加，可有效缓解交通碳排放压力；相同机动车保有量增长率条件下，交通碳排放量随着GDP增长速度加快逐渐增大，即一味追求经济发展会带来更多交通运输业的减排压力；随时间推移，交通碳排放总量增长趋势渐缓，年均增长率由9.20%下降到6.12%，说明在相关部门的综合治理下，低碳节能是社会发展的必然趋势.

 

表4 不同条件下的碳排放预测值Table 4 Forecast value of carbon emissions under different situations

  

GDP增速/%量增长率/%2017年2018年2019年2020年0.53 115.467 73 295.690 63 489.246 33 695.690 4 1.53 122.100 43 309.116 73 509.373 43 722.630 1 2.53 128.729 83 322.528 43 529.922 93 750.293 4 0.53 117.970 73 300.879 13 497.290 63 706.755 6 1.53 124.602 73 314.302 53 517.410 93 733.682 1 2.53 131.231 63 327.711 43 537.953 43 761.331 5 0.53 120.472 83 306.085 43 505.401 73 717.959 0 1.53 127.104 33 319.506 13 525.515 13 744.871 9 2.53 133.732 53 332.912 23 546.050 53 772.506 9碳排放预测值/万t 6.0 6.5 7.0机动车保有

3 结论

基于IPCC清单指南的相关数据，确定区域交通碳排放相关指标的测算方法和预测模型，并以北京市为例进行实证分析，得出以下结论：

有关古代建筑的文字记录资料可分为3种：金石文字、题记与古代文献记载[4]4。有关旌义坊的文字记录，仅有金石文字和古代文献记载。

(1)建立支持向量回归机模型对区域交通碳排放进行预测研究，以北京市为例，交叉验证均方误差仅为0.008 040，训练和测试拟合相关系数分别为0.984 2和0.995 0，均方误差分别为0.002 369和0.001 163，可见模型具有很好的学习和推广能力.

(2)通过预测结果可知，区域交通碳排放总量整体呈上升趋势，社会仍然面临较大温室气体减排压力；但随时间推移，交通碳排放总量增长趋势逐渐变缓，这表明政府节能减排目标已初见成效.以北京市为例，到2020年碳排放总量平均增长到3 733.969 0万t，年均增长率则由9.20%下降到6.12%.

参考文献:

关于能源结构，本文借鉴产业结构指数构造法，将其定义为不同能源(本文重点研究电力和天然气两类清洁能源比例)消耗转换标准煤用量与行业消耗能源标准煤总量的比值，计算方法为

对训练样本集和测试样本集数据进行拟合分析，得到反归一化后的预测值与实际值线性回归结果如图2和图3所示.其中，训练结果线性回归方程为Q=0.963 0M+58.308 9，相关系数为0.984 2，均方误差为0.002 369；测试结果线性回归方程为Q=0.951 8M+66.107 4，相关系数为0.995 0，均方误差为0.001 164.对应预测值与实际值数据对比结果如图4和图5所示，可见模型具有良好的学习推广能力，可用作区域交通碳排放预测研究.

[3]陈亮,何涛,李巧茹,等.区域交通碳排放相关指标测算及影响因素分析[J].北京工业大学学报,2017,43(4):631-637.[CHEN L,HE T,LI Q R,et al.Index calculation and effect factor study of regional transport carbon emission[J].Journal of Beijing University of Technology,2017,43(4):631-637.]

[4]WANG T,LI H,ZHANG J,et al.Influencing factors of carbon emission in China's road freight transport[J].Procedia-Social and Behavioral Sciences,2012(43):54-64.

[5]卢建锋,傅惠,王小霞.区域交通运输业碳排放效率影响因素研究[J].交通运输系统工程与信息,2016,16(2):25-30.[LU J F,FU H,WANG X X.Research on the impact of regional transportation emissions efficiency factors[J].Journal of Transportation Systems Engineering and Information Technology,2016,16(2):25-30.]

[6]宋震,丛林.中国交通运输业能源效率及其影响因素研究[J].交通运输系统工程与信息,2016,16(1):19-25.[SONG Z,CONG L.Energyefficiencyand influencing factors of Chinese transportation industry[J].Journal of Transportation Systems Engineering and Information Technology,2016,16(1):19-25.]

[7]刘鸿远.城市交通碳排放预测与发展模式选择[D].天津:天津理工大学,2014.[LIU H Y.The forecast of carbon emission of urban transportation and choice of development mode take Tianjin as an example[D].Tianjin:Tianjin University of Technology,2014.]

[8]KAYA Y.Impacts of carbon dioxide emission control on GDP growth:interpretation of proposed scenarios[R].Paris:Presentation to the Energy and Industry Subgroup,Response Strategies Working Group,IPCC,1989:1-25.

[9]YORK R,ROSA E A,DIETZ T.STIRPAT,IPAT and ImPACT:analytic tools for unpacking the driving forces ofenvironmentalimpacts[J].EcologicalEconomics,2003,46(3):351-365.

[10]刘建翠.中国交通运输部门节能潜力和碳排放预测[J].资源科学,2011,33(4):640-646.[LIU J C.Energy saving potential and carbon emissions prediction for the transportation sector in China[J].Resources Science,2011,33(4):640-646.]

由于在此系统中，引入了电压负反馈，而负反馈电路要求反馈到输入端的输出信号与原输入信号同相。但因为电路中零极点的存在，不可避免地引入了相移，当相移达到180。时，会使系统产生自激振荡，系统变得不稳定。为了确定系统的稳定性，需要对其进行Bode图的分析。保证系统具有可靠的稳定性下，需要使Bode图中的幅值裕度大于10 dB，相位裕度大于45 deg。为了获得足够的系统稳定性，使系统运行稳定，还须满足下述条件[12-16]：

表4为上海与前10位主办城市网络中心度、中介度测算值的比对情况。从中可以发现，伦敦和墨尔本稳居第一等级，定为Alpha级；纽约、东京、巴黎位居第二等级，定为Alpha-级；温哥华、柏林、悉尼、伊斯坦布尔、多伦多位居第三等级，定为Beta级。上海与上述10个城市相比，各个指标上均存在较大的差距。上海国际著名体育城市的发展目标应该至少定位在Alpha-级，即向伦敦、墨尔本、纽约、东京、巴黎靠齐。

[11]宋杰鲲.基于支持向量回归机的中国碳排放预测模型[J].中国石油大学学报（自然科学版）,2012,36(1):182-187.[SONG J K.China’s carbon emissions prediction model based on support vector regression[J].Journal of China University of Petroleum(Edition of Natural Science),2012,36(1):182-187.]

[12]WOOLF D,AMONETTE J E,STREETPERROTT F A,et al.Sustainable biochar to mitigate global climate change[J].Nature Communications,2010,1(5):56.

[13]CORTES C,VAPNIK V.Support-vector networks[J].Machine Learning,1995,20(3):273-297.