CN114118509B - 一种道路交通二氧化碳排放的预测方法及预测系统 - Google Patents
一种道路交通二氧化碳排放的预测方法及预测系统 Download PDFInfo
- Publication number
- CN114118509B CN114118509B CN202110923903.7A CN202110923903A CN114118509B CN 114118509 B CN114118509 B CN 114118509B CN 202110923903 A CN202110923903 A CN 202110923903A CN 114118509 B CN114118509 B CN 114118509B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- new energy
- vehicles
- target year
- year
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 title claims abstract description 192
- 229910002092 carbon dioxide Inorganic materials 0.000 title claims abstract description 96
- 239000001569 carbon dioxide Substances 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000004083 survival effect Effects 0.000 claims abstract description 97
- 239000000446 fuel Substances 0.000 claims abstract description 80
- 238000010276 construction Methods 0.000 claims description 36
- 230000035699 permeability Effects 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 16
- 230000014759 maintenance of location Effects 0.000 claims description 12
- 239000000126 substance Substances 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000007423 decrease Effects 0.000 claims description 4
- 238000013077 scoring method Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 description 15
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 6
- 229910052799 carbon Inorganic materials 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000002567 autonomic effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
- Y02P90/845—Inventory and reporting systems for greenhouse gases [GHG]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Geometry (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Primary Health Care (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Automation & Control Theory (AREA)
- Educational Administration (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于交通运输行业二氧化碳排放量(碳排放)核算和二氧化碳排放空间格局预测技术领域,具体地说,涉及一种道路交通二氧化碳排放的预测方法,该方法包括:针对未来目标年份的待测区域,获取未来目标年份的加权路网密度空间分布密度数据集;得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率;计算未来目标年份的机动车的保有量,同时,得到未来目标年份的机动车的燃油车保有量;将来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量,以及燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排放空间网格数据。
Description
技术领域
本发明属于交通运输行业二氧化碳排放量(碳排放)核算和二氧化碳排放空间格局预测技术领域,具体地说,涉及一种道路交通二氧化碳排放的预测方法及预测系 统。
背景技术
全球二氧化碳排放量约有25%来自交通运输,我国道路交通的碳排放在交通运输中的比重超过90%,随着交通快速机动化发展阶的到来,道路交通对区域碳排放 的贡献将进一步增加。目前,道路交通二氧化碳排放预测主要以指数分解等自上而下 的分析方法为主,通过研究社会、经济、能源、排放等宏观总体因子之间的关系,来 预测未来道路交通二氧化碳排放量的变化趋势,缺乏自下而上基于演化思维的排放 路径建立方法。因此,如何将自上而下的宏观思维和自下而上的微观自主性思维相结 合,对适用于道路交通二氧化碳进行路径分析,有利于精确化制定区域减排政策,助 力城市治理体系和治理能力现代化。
道路交通二氧化碳排放路径建立方法具有重要意义:
(1)路径建立方法可以最大可能分析不同社会经济发展、政策设置情景下的减 排潜力,识别中长期深度减排的关键驱动技术与政策方案,有利于加快推广成熟减排 技术促进高质量达峰实现,并提前部署深度减排需要的各项技术和政策工具;
(2)路径分析不仅能够大幅减少温室气体排放,还能有效分析减排行为的效益,评估减排措施及其效益,保证稳定性,助力气候与经济发展的共赢;
(3)融合自上而下的宏观思维和自下而上的微观自主性思维路径分析可以有效建立一个全区域的宏观-中观-微观的减排分析机制,能实现多尺度、跨区域的减排潜 力分析,有利于多尺度减排目标的合理、公平、高效设置。
为通过节能减排,传统汽油车向新能源汽车的转型是一项重要的举措,其中,电动汽车因其节能减排的优势将在这次转型中发挥重要作用。电动汽车的能源利用效 率比传统燃油汽车高出46%以上,并具有13%至68%的二氧化碳减排潜力。尽管 交通电气化被认为是交通部门的重要组成部分,但仍然缺乏详细的中期和长期目标, 特别是缺少电动汽车推广战略及其减排作用方式。因此,需要对道路交通碳排放未来 碳排放的空间格局的进行更高精度的研究,以此实现区域道路交通二氧化碳排放达 峰目标的合理制定。
对于道路交通二氧化碳排放预测,通常是通过构架不同社会发展情景下的机动车数量变化,使用相关统计模型,进行碳排放在时间序列上的预测;但是,现有的预 测方法存在如下问题:
(1)未充分考虑机动车存活率和新能源车销售渗透率的发展趋势和特点差异, 无法反应传统燃油车逐步报废和道路交通电气化推广下的道路交通二氧化碳排放的 时空差异;
(2)不能同时识别道路交通二氧化碳排放在时间上和空间上的异质性,无法准 确判断不同区域的达峰时间,无法针对预测结果,提出更具针对性的减排目标。
发明内容
为解决现有技术存在的上述缺陷,本发明提出了一种道路交通二氧化碳排放的预测方法,具体地,涉及一种耦合自上而下和自下而上相结合,基于空间演化理论的 道路交通二氧化碳排放预测方法,通过城市建设用地扩张和人口变化,预测加权路网 空间变化,根据车辆存活曲线和新能源车销售渗透率预测道路交通电气化率,根据加 权路网密度,分配燃油车数量至空间网格,通过建立燃油车数量与二氧化碳排放之间 的函数关系,预测未来年的道路交通二氧化碳排放量。
本发明提供了一种道路交通二氧化碳排放的预测方法,该方法包括:
针对未来目标年份的待测区域,将待测区域在空间上划分为多个基准网格,每个基准网格为10km×10km的空间基准网格,结合人口空间分布预测结果和城市建设用 地空间扩张预测结果,获取每一个基准网格的加权路网空间分布密度数据,进而获取 未来目标年份的加权路网密度空间分布密度数据集;
根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能源车的车龄存活概率;
根据Gompertz预测模型和增长率预测模型,计算未来目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并结合新能源车的车龄存活 概率,得到目标年份的新能源车保有量,利用目标年份的机动车总量减去目标年份的 新能源车保有量,得到未来目标年份的机动车的燃油车保有量;
将来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量,以及燃油车数 量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排放空间网格数 据,实现对道路交通二氧化碳排放的预测。
作为上述技术方案的改进之一,所述目标年份的加权路网密度空间分布密度数据集的具体获取过程为:
采用基于共享社会经济路径开发的全球2010-2100年间每十年的人口网格预测数据,按照预先设定的行政级别,提取不同年份的人口数据,再结合待测区域的实际 情况的逐年总人口预测数据,得到符合实际的人口分布数据集D1,作为人口空间分 布预测结果;
根据历史建设用地扩张空间数据,以预先设定的基准年2020年60%的土地利用现状数据为基础,基于热点阈值打分的方法,得到目标年份城市建设用地空间分布数 据集D2,作为城市建设用地空间扩张预测结果;
从城市建设用地空间扩张预测结果中,提取目标年份城市建设用地扩张面积;并结合对应的人口空间分布预测结果,获取第i个基准网格的加权路网空间分布密度数 据Roadni:
其中,Roadni代表第n目标年份的第i个基准网格的加权路网空间分布密度数 据;其中,在本实施例中,n为2025,2030,2035,2040,2050和2060年;
ULni为第n目标年份的第i个基准网格的城市建设用地扩张面积;UL2020i为2020 年的第i个基准网格的建设用地扩张面积;Popni为第n目标年份的第i个基准网格 的人口数据;Pop2020i为2020年的第i个基准网格的人口数据;Road2020i为2020年 第i个基准网格的加权路网空间分布密度数据;
将通过上式,将获取的每一个基准网格的加权路网空间分布密度数据进行汇总,得到目标年份加权路网空间分布密度数据集D3。
作为上述技术方案的改进之一,所述根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新 能源车的车龄存活概率;其具体过程为:
假设机动车分为新能源车和燃油车两大类;其中,新能源车包括:乘用车和商用车;燃油车包括:乘用车和商用车;
根据新能源车的历史年份2000-2020每年的注册量和保有量,分别拟合得到新 能源车乘用车和商用车存活曲线,该曲线反映新车群体随着车龄增加其存活比例逐 渐下降的规律:
新能源车中的乘用车的存活概率:
其中,为第k目标年份,车型j为乘用车,车龄为i的新能源车中的乘 用车的存活概率;n0,j(k)1为第k目标年份车型j为乘用车的注册量;ni,j(k+1)1为 车型j为乘用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新乘用车保有 量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度;Tj表示新能源车中的车型j为乘用车的寿命;
新能源车中的商用车的存活概率:
其中,为第k目标年份,车型j为商用车,车龄为i的新能源车中的商用 车的存活概率;n0,j(k)2为第k目标年份车型j为商用车的注册量;ni,j(k+1)2为车 型j为商用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新商用车保有量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度,Tj1表示车型j为新能源车中的商用车的寿命。
作为上述技术方案的改进之一,所述根据Gompertz预测模型和增长率预测模型,计算未来目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗 透率,并结合新能源车的车龄存活概率,得到目标年份的新能源车保有量,利用目标 年份的机动车总量减去目标年份的新能源车保有量,得到未来目标年份的机动车的 燃油车保有量;其具体过程为:
获取未来目标年份的机动车保有量的具体过程如下:
假设机动车的车型分为乘用车和商用车;根据Gompertz预测模型,计算未来目 标年份的乘用车保有量:
其中,Vp为未来目标年份的乘用车的保有量,其单位为辆/千人;K为饱和保有量;a为第一常数值,b为第二常数值,二者均为已知的常数值;g为人均GDP,其单位 为元/人;
根据增长率预测模型,计算未来目标年份的商用车的保有量:
V=a′+b′×H (7)
其中,V为未来目标年份的商用车的保有量,其单位为万辆;a′、b′为常数值,根 据2002-2010年商用货车保有量和货运量的线性关系拟合得到;H为货运量,万吨;
则未来目标年份的机动车保有量A:
A=Vp+V; (8)
未来目标年份的新能源车保有量的具体计算过程如下:
假设新能源车包括:乘用车和商用车;
基于基准年的新能源车车龄分布,根据未来每年的新能源车的新车渗透率,预测未来年份每年的新能源车新车注册量,结合车辆存活概率曲线,进而得到目标年份的 乘用车与商用车的新能源车的保有量:
未来目标年份新能源车中的乘用车保有量NVPk:
NVPk=∑j∑iNVPi,j,k (10)
其中,NVPi,j,k为k年的乘用车新能源车车龄分布,表示k年车龄为i的乘用车新 能源车保有量,NVPk为k年新能源车中的乘用车保有量,NSPj,k-i为第k-i年的新能 源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年(2020年) 的车辆车龄分布,即NVP2020,为k年、车型j为新能源乘用车、车龄为i的 存活概率;
其中,
NSPj,k-i=SPk-i×MSPk-i (11)
其中,SPk-i为k-i年乘用车的总销售量,MSPk-i为第k-i年的新能源车的新车 渗透率;
未来目标年份新能源车中的商用车保有量NVk:
NVk=∑j∑iNVi,j,k (13)
其中,NVi,j,k为k年的商用车新能源车车龄分布,表示k年车龄为i的商用车新 能源车保有量,NVk为k年新能源车中的商用车保有量,NSj,k-i为第k-i年的商用车 新能源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年(2020 年)的车辆车龄分布,即NV2020,为k年j车型为商用车新能源车、车龄为i 的存活概率;
NSj,k-i=Sk-i×MSk-i (14)
其中,Sk-i为k-i年商用车的总销售量,MSk-i为第k-i年的商用车新能源车的 新车渗透率;
目标年份的新能源车保有量等于未来目标年份新能源车中的乘用车保有量NVPk加上未来目标年份新能源车中的商用车保有量NVk;
将未来目标年份的机动车保有量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量:
具体地,未来目标年份的机动车保有量A和k年新能源车乘用车和商用的保有量NVPk、NVk,得到未来目标年份的机动车的燃油车保有量FVk;
FVk=A-NVPk-NVk (15)。
作为上述技术方案的改进之一,所述将未来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基 准网格的燃油车保有量,以及燃油车数量与二氧化碳排放的函数关系,预测目标年份 的道路交通二氧化碳排放空间网格数据,实现对道路交通二氧化碳排放的预测;其具 体过程为:
针对待测区域在空间上划分的多个10km×10km的基准网格,将未来目标年份的机动车的燃油车保有量分配至每个基准网格,得到未来目标年份的各个基准网格的 机动车的燃油车保有量Vnj:
其中,Vnj为目标年份的第n个基准网格的机动车的燃油车保有量;′RLS2020j为 基准年份的第n个基准网格的加权路网空间分布密度数据,RLS2020j为基准年份的n 个基准网格的加权路网空间分布密度数据的总和;Vn为目标年份的机动车的燃油车 的保有量,Vn=FVk;
基于燃油车数量与二氧化碳排放的函数关系:
其中,Cn为预测的第n个基准网格的目标年份的道路交通二氧化碳排放空间网 格数据;a1为第一权重值;b1为第二权重值;c为常数;RLSn为目标年份的n个基 准网格的加权路网空间分布密度数据的总和;RLS2020为目标年份的第n个基准网格 的加权路网空间分布密度数据,即目标年份加权路网空间分布密度数据集D3; V2020为基准年份的第n个基准网格的机动车的燃油车的保有量;C2020为基准年份 的二氧化碳排放量,记为已知值;
根据得到的Cn,获取目标年份的道路交通二氧化碳排放空间网格数据C:
C=n×Cn (18)
其中,n为针对待测区域在空间上所划分的基准网格的数量;
并将上述计算结果作为预测结果,再根据二氧化碳排放总量与预测结果的差异进行倍数调整,得到横坐标为时间、纵坐标为二氧化碳排放量的变化曲线,根据该变 化曲线的变化率平滑得到逐年的二氧化碳排放量。
本发明还提供了一种道路交通二氧化碳排放的预测系统,该系统包括:
密度数据获取模块,用于针对未来目标年份的待测区域,将待测区域在空间上划分为多个基准网格,每个基准网格为10km×10km的空间基准网格,结合人口空间分 布预测结果和城市建设用地空间扩张预测结果,获取每一个基准网格的加权路网空 间分布密度数据,进而获取未来目标年份的加权路网密度空间分布密度数据集;
新能源车存活概率获取模块,用于根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能 源车的车龄存活概率;
燃油车保有量获取模块,根据Gompertz预测模型和增长率预测模型,计算未来 目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并 结合新能源车的车龄存活概率,得到目标年份的新能源车保有量,利用目标年份的机 动车总量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保 有量;和
预测模块,用于将来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量, 以及燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排 放空间网格数据,实现对道路交通二氧化碳排放的预测。
本发明与现有技术相比的有益效果是:
1、本发明的方法充分考虑了车辆存活曲线与新能源车销售渗透率,结合车辆 存活曲线和新能源车销售渗透率的变化测算道路交通二氧化碳排放时,充分考虑了 传统燃油车逐步淘汰和新能源车的推广,能够更有效地揭示道路交通电气化推广下 的道路交通二氧化碳排放的时空差异;
2、空间模型的构建能够同时识别道路交通二氧化碳排放在时间上和空间上的异质性,有助于筛选出该排放的典型区域,提出更具针对性的减排目标;
3、在空间网格层面上,实现任意网格区域内的道路交通二氧化碳排放的预测, 从而完成对不同区域的道路交通二氧化碳排放达峰的快速预测。
附图说明
图1是本发明的一种道路交通二氧化碳排放的预测方法的流程图;
图2是图1的本发明的一种道路交通二氧化碳排放的预测方法的步骤1)的流程图。
具体实施方式
现结合附图和实例对本发明作进一步的描述。
如图1所示,本发明提供了一种道路交通二氧化碳排放的预测方法,该方法具 体包括:
步骤1)针对未来目标年份的待测区域,将待测区域在空间上划分为多个基准网格,每个基准网格为10km×10km的空间基准网格,结合人口空间分布预测结果和城 市建设用地空间扩张预测结果,获取每一个基准网格的加权路网空间分布密度数据, 进而获取未来目标年份的加权路网密度空间分布密度数据集;
具体地,如图2所示,采用基于共享社会经济路径(SSP)开发的全球2010-2100 年间每十年的人口网格预测数据,按照预先设定的行政级别,提取不同年份的人口数 据,再结合待测区域的实际情况的逐年总人口预测数据,得到符合实际的人口分布数 据集D1,作为人口空间分布预测结果;
根据历史建设用地扩张空间数据,以预先设定的基准年(在本实施例中为2020年)60%的土地利用现状数据为基础,基于热点阈值打分的方法(2025年全国65% 城镇化水平、2030年全国70%城镇化水平、2035年全国75%城镇化水平),得到目 标年份城市建设用地空间分布数据集D2,作为城市建设用地空间扩张预测结果;
从城市建设用地空间扩张预测结果中,提取目标年份城市建设用地扩张面积;并结合对应的人口空间分布预测结果,获取第i个基准网格的加权路网空间分布密度数 据Roadni:
其中,Roadni代表第n目标年份的第i个基准网格的加权路网空间分布密度数 据;其中,在本实施例中,n为2025,2030,2035,2040,2050和2060年;
ULni为第n目标年份的第i个基准网格的城市建设用地扩张面积;UL2020i为2020 年的第i个基准网格的建设用地扩张面积;Popni为第n目标年份的第i个基准网格 的人口数据;Pop2020i为2020年的第i个基准网格的人口数据;Road2020i为2020年 第i个基准网格的加权路网空间分布密度数据;
将通过上式将获取的每一个基准网格的加权路网空间分布密度数据进行汇总,得到目标年份加权路网空间分布密度数据集D3。
步骤2)根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能源车的车龄存活概率;
具体地,假设机动车分为新能源车和燃油车两大类;其中,新能源车包括:乘用 车和商用车;燃油车包括:乘用车和商用车;
根据新能源车的历史年份2000-2020每年的注册量和保有量,分别拟合得到新 能源车乘用车和商用车存活曲线,该曲线反映新车群体随着车龄增加其存活比例逐 渐下降的规律:
新能源车中的乘用车的存活概率:
其中,为第k目标年份,车型j为乘用车,车龄为i的新能源车中的乘 用车的存活概率;n0,j(k)1为第k目标年份车型j为乘用车的注册量;ni,j(k+1)1为 车型j为乘用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新乘用车保有 量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度;Tj表示新能源车中的车型j为乘用车的寿命;
新能源车中的商用车的存活概率:
其中,为第k目标年份,车型j为商用车,车龄为i的新能源车中的商用 车的存活概率;n0,j(k)2为第k目标年份车型j为商用车的注册量;ni,j(k+1)2为车 型j为商用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新商用车保有量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度,Tj1表示车型j为新能源车中的商用车的寿命。
步骤3)根据Gompertz预测模型和增长率预测模型,计算未来目标年份的机动 车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并结合新能源车的 车龄存活概率,得到目标年份的新能源车保有量,利用目标年份的机动车总量减去目 标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量;
具体地,获取未来目标年份的机动车保有量的具体过程如下:
假设机动车的车型分为乘用车和商用车;根据Gompertz预测模型,计算未来目 标年份的乘用车保有量:
其中,Vp为未来目标年份的乘用车的保有量,其单位为辆/千人;K为饱和保有量;a为第一常数值,b为第二常数值,二者均为已知的常数值;g为人均GDP,其单位 为元/人;
根据增长率预测模型,计算未来目标年份的商用车的保有量:
V=a′+b′×H (7)
其中,V为未来目标年份的商用车的保有量,其单位为万辆;a′、b′为常数值,根 据2002-2010年商用货车保有量和货运量的线性关系拟合得到;H为货运量,万吨;
则未来目标年份的机动车保有量A:
A=Vp+V; (8)
未来目标年份的新能源车保有量的具体计算过程如下:
假设新能源车包括:乘用车和商用车;
基于基准年的新能源车车龄分布,根据未来每年的新能源车的新车渗透率,预测未来年份每年的新能源车新车注册量,结合车辆存活概率曲线,进而得到目标年份的 乘用车与商用车的新能源车的保有量:
未来目标年份新能源车中的乘用车保有量NVPk:
NVPk=∑j∑iNVPi,j,k (10)
其中,NVPi,j,k为k年的乘用车新能源车车龄分布,表示k年车龄为i的乘用车新 能源车保有量,NVPk为k年新能源车中的乘用车保有量,NSPj,k-i为第k-i年的新能 源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年(2020年) 的车辆车龄分布,即NVP2020,为k年、车型j为新能源乘用车、车龄为i的 存活概率;
其中,
NSPj,k-i=SPk-i×MSPk-i (11)
其中,SPk-i为k-i年乘用车的总销售量,MSPk-i为第k-i年的新能源车的新车 渗透率;
未来目标年份新能源车中的商用车保有量NVk:
NVk=∑j∑iNVi,j,k (13)
其中,NVi,j,k为k年的商用车新能源车车龄分布,表示k年车龄为i的商用车新 能源车保有量,NVk为k年新能源车中的商用车保有量,NSj,k-i为第k-i年的商用车 新能源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年(2020 年)的车辆车龄分布,即NV2020,为k年j车型为商用车新能源车、车龄为i 的存活概率;
NSj,k-i=Sk-i×MSk-i (14)
其中,Sk-i为k-i年商用车的总销售量,MSk-i为第k-i年的商用车新能源车的 新车渗透率;
目标年份的新能源车保有量等于未来目标年份新能源车中的乘用车保有量NVPk加上未来目标年份新能源车中的商用车保有量NVk;
将未来目标年份的机动车保有量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量:
具体地,未来目标年份的机动车保有量A和k年新能源车乘用车和商用的保有量NVPk、NVk,得到未来目标年份的机动车的燃油车保有量FVk;
FVk=A-NVPk-NVk (15)。
步骤4)将未来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量,以及 燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排放空 间网格数据,实现对道路交通二氧化碳排放的预测。
具体地,
针对待测区域在空间上划分的多个10km×10km的基准网格,将未来目标年份的机动车的燃油车保有量分配至每个基准网格,得到未来目标年份的各个基准网格的 机动车的燃油车保有量Vnj:
其中,Vnj为目标年份的第n个基准网格的机动车的燃油车保有量;′RLS2020j为 基准年份的第n个基准网格的加权路网空间分布密度数据,RLS2020j为基准年份的n 个基准网格的加权路网空间分布密度数据的总和;Vn为目标年份的机动车的燃油车 的保有量,Vn=FVk;
基于燃油车数量与二氧化碳排放的函数关系:
其中,Cn为预测的第n个基准网格的目标年份的道路交通二氧化碳排放空间网 格数据;a1为第一权重值;b1为第二权重值;c为常数;RLSn为目标年份的n个基 准网格的加权路网空间分布密度数据的总和;RLS2020为目标年份的第n个基准网格 的加权路网空间分布密度数据,即目标年份加权路网空间分布密度数据集D3; V2020为基准年份的第n个基准网格的机动车的燃油车的保有量;C2020为基准年份 的二氧化碳排放量,记为已知值;
根据得到的Cn,获取目标年份的道路交通二氧化碳排放空间网格数据C:
C=n×Cn (18)
其中,n为针对待测区域在空间上所划分的基准网格的数量;
并将上述计算结果作为预测结果,再根据二氧化碳排放总量与预测结果的差异进行倍数调整,得到横坐标为时间、纵坐标为二氧化碳排放量的变化曲线,根据该变 化曲线的变化率平滑得到逐年的二氧化碳排放量。
本发明还提供了一种道路交通二氧化碳排放的预测系统,该系统包括:
密度数据获取模块,用于针对未来目标年份的待测区域,将待测区域在空间上划分为多个基准网格,每个基准网格为10km×10km的空间基准网格,结合人口空间分 布预测结果和城市建设用地空间扩张预测结果,获取每一个基准网格的加权路网空 间分布密度数据,进而获取未来目标年份的加权路网密度空间分布密度数据集;
新能源车存活概率获取模块,用于根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能 源车的车龄存活概率;
燃油车保有量获取模块,根据Gompertz预测模型和增长率预测模型,计算未来 目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并 结合新能源车的车龄存活概率,得到目标年份的新能源车保有量,利用目标年份的机 动车总量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保 有量;和
预测模块,用于将来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量, 以及燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排 放空间网格数据,实现对道路交通二氧化碳排放的预测。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管 参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明 的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均 应涵盖在本发明的权利要求范围当中。
Claims (2)
1.一种道路交通二氧化碳排放的预测方法,该方法包括:
针对未来目标年份的待测区域,将待测区域在空间上划分为多个基准网格,每个基准网格为10km×10km的空间基准网格,结合人口空间分布预测结果和城市建设用地空间扩张预测结果,获取每一个基准网格的加权路网空间分布密度数据,进而获取未来目标年份的加权路网密度空间分布密度数据集;
根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能源车的车龄存活概率;
根据Gompertz预测模型和增长率预测模型,计算未来目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并结合新能源车的车龄存活概率,得到目标年份的新能源车保有量,利用目标年份的机动车总量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量;
将未来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量,以及燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排放空间网格数据,实现对道路交通二氧化碳排放的预测;
所述目标年份的加权路网密度空间分布密度数据集的具体获取过程为:
采用基于共享社会经济路径开发的全球2010-2100年间每十年的人口网格预测数据,按照预先设定的行政级别,提取不同年份的人口数据,再结合待测区域的实际情况的逐年总人口预测数据,得到符合实际的人口分布数据集D1,作为人口空间分布预测结果;
根据历史建设用地扩张空间数据,以预先设定的基准年2020年60%的土地利用现状数据为基础,基于热点阈值打分的方法,得到目标年份城市建设用地空间分布数据集D2,作为城市建设用地空间扩张预测结果;
从城市建设用地空间扩张预测结果中,提取目标年份城市建设用地扩张面积;并结合对应的人口空间分布预测结果,获取第i个基准网格的加权路网空间分布密度数据Roadni:
其中,Roadni代表第n目标年份的第i个基准网格的加权路网空间分布密度数据;其中,n为2025,2030,2035,2040,2050和2060年;
ULni为第n目标年份的第i个基准网格的城市建设用地扩张面积;UL2020i为2020年的第i个基准网格的建设用地扩张面积;Popni为第n目标年份的第i个基准网格的人口数据;Pop2020i为2020年的第i个基准网格的人口数据;Road2020i为2020年第i个基准网格的加权路网空间分布密度数据;
将通过上式,将获取的每一个基准网格的加权路网空间分布密度数据进行汇总,得到目标年份加权路网空间分布密度数据集D3;
所述根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能源车的车龄存活概率,其具体过程为:
假设机动车分为新能源车和燃油车两大类;其中,新能源车包括:乘用车和商用车;燃油车包括:乘用车和商用车;
根据新能源车的历史年份2000-2020每年的注册量和保有量,分别拟合得到新能源车乘用车和商用车存活曲线,该曲线反映新能源车群体随着车龄增加其存活比例逐渐下降的规律:
新能源车中的乘用车的存活概率:
其中,为第k目标年份,车型j为乘用车,车龄为i的新能源车中的乘用车的存活概率;n0,j(k)1为第k目标年份车型j为乘用车的注册量;ni,j(k+1)1为车型j为乘用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新乘用车保有量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度;Tj表示新能源车中的车型j为乘用车的寿命;
新能源车中的商用车的存活概率:
其中,为第k目标年份,车型j为商用车,车龄为i的新能源车中的商用车的存活概率;n0,j(k)2为第k目标年份车型j为商用车的注册量;ni,j(k+1)2为车型j为商用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新商用车保有量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度,Tj1表示车型j为新能源车中的商用车的寿命;
所述根据Gompertz预测模型和增长率预测模型,计算未来目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并结合新能源车的车龄存活概率,得到目标年份的新能源车保有量,利用目标年份的机动车总量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量,其具体过程为:
获取未来目标年份的机动车保有量的具体过程如下:
假设机动车的车型分为乘用车和商用车;根据Gompertz预测模型,计算未来目标年份的乘用车保有量:
其中,Vp为未来目标年份的乘用车的保有量,其单位为辆;K为饱和保有量;a为第一常数值,b为第二常数值,二者均为已知的常数值;g为人均GDP,其单位为元/人;
根据增长率预测模型,计算未来目标年份的商用车的保有量:
V=a′+b′×H (7)
其中,V为未来目标年份的商用车的保有量,其单位为万辆;a′、b′为常数值,根据2002-2010年商用货车保有量和货运量的线性关系拟合得到;H为货运量,其单位为万吨;
则未来目标年份的机动车保有量A:
A=Vp+V; (8)
未来目标年份的新能源车保有量的具体计算过程如下:
假设新能源车包括:乘用车和商用车;
基于基准年的新能源车车龄分布,根据未来每年的新能源车的新车渗透率,预测未来年份每年的新能源车新车注册量,结合车辆存活概率曲线,进而得到目标年份的乘用车与商用车的新能源车的保有量:
未来目标年份新能源车中的乘用车保有量NVPk:
NVPk=∑j∑iNVPi,j,k (10)
其中,NVPi,j,k为k年的乘用车新能源车车龄分布,表示k年车龄为i的乘用车新能源车保有量,NVPk为k年新能源车中的乘用车保有量,NSPj,k-i为第k-i年的新能源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年2020年的车辆车龄分布,即NVP2020,为k年、车型j为新能源乘用车、车龄为i的存活概率;
其中,
NSPj,k-i=SPk-i×MSPk-i (11)
其中,SPk-i为k-i年乘用车的总销售量,MSPk-i为第k-i年的新能源车的新车渗透率;
未来目标年份新能源车中的商用车保有量NVk:
NVk=∑j∑iNVi,j,k (13)
其中,NVi,j,k为k年的商用车新能源车车龄分布,表示k年车龄为i的商用车新能源车保有量,NVk为k年新能源车中的商用车保有量,NSj,k-i为第k-i年的商用车新能源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年2020年的车辆车龄分布,即NV2020,为k年j车型为商用车新能源车、车龄为i的存活概率;
NSj,k-i=Sk-i×MSk-i (14)
其中,Sk-i为k-i年商用车的总销售量,MSk-i为第k-i年的商用车新能源车的新车渗透率;
目标年份的新能源车保有量等于未来目标年份新能源车中的乘用车保有量NVPk加上未来目标年份新能源车中的商用车保有量NVk;
将未来目标年份的机动车保有量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量:
具体地,未来目标年份的机动车保有量A和k年新能源车乘用车和商用车的保有量NVPk、NVk,得到未来目标年份的机动车的燃油车保有量FVk;
FVk=A-NVPk-NVk (15);
所述将未来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量,以及燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排放空间网格数据,实现对道路交通二氧化碳排放的预测,其具体过程为:
针对待测区域在空间上划分的多个10km×10km的基准网格,将未来目标年份的机动车的燃油车保有量分配至每个基准网格,得到未来目标年份的各个基准网格的机动车的燃油车保有量Vnj:
其中,Vnj为目标年份的第n个基准网格的机动车的燃油车保有量;′RLS2020j为基准年份的第n个基准网格的加权路网空间分布密度数据,RLS2020j为基准年份的n个基准网格的加权路网空间分布密度数据的总和;Vn为目标年份的机动车的燃油车的保有量,Vn=FVk;
基于燃油车数量与二氧化碳排放的函数关系:
其中,Cn为预测的第n个基准网格的目标年份的道路交通二氧化碳排放空间网格数据;a1为第一权重值;b1为第二权重值;c为常数;RLSn为目标年份的n个基准网格的加权路网空间分布密度数据的总和;RLS2020为目标年份的第n个基准网格的加权路网空间分布密度数据,即目标年份加权路网空间分布密度数据集D3;V2020为基准年份的第n个基准网格的机动车的燃油车的保有量;C2020为基准年份的二氧化碳排放量,记为已知值;
根据得到的Cn,获取目标年份的道路交通二氧化碳排放空间网格数据C:
C=n×Cn (18)
其中,n为针对待测区域在空间上所划分的基准网格的数量;
并将所述目标年份的道路交通二氧化碳排放空间网格数据C作为预测结果,再根据二氧化碳排放总量与预测结果的差异进行倍数调整,得到横坐标为时间、纵坐标为二氧化碳排放量的变化曲线,根据该变化曲线的变化率平滑得到逐年的二氧化碳排放量。
2.一种道路交通二氧化碳排放的预测系统,其特征在于,该系统基于上述权利要求1所述的道路交通二氧化碳排放的预测方法实现,该系统包括:
密度数据获取模块,用于针对未来目标年份的待测区域,将待测区域在空间上划分为多个基准网格,每个基准网格为10km×10km的空间基准网格,结合人口空间分布预测结果和城市建设用地空间扩张预测结果,获取每一个基准网格的加权路网空间分布密度数据,进而获取未来目标年份的加权路网密度空间分布密度数据集;
新能源车存活概率获取模块,用于根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能源车的车龄存活概率;
燃油车保有量获取模块,根据Gompertz预测模型和增长率预测模型,计算未来目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并结合新能源车的车龄存活概率,得到目标年份的新能源车保有量,利用目标年份的机动车总量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量;和
预测模块,用于将未来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量,以及燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排放空间网格数据,实现对道路交通二氧化碳排放的预测;
所述密度数据获取模块中,对于目标年份的加权路网密度空间分布密度数据集的具体获取过程为:
采用基于共享社会经济路径开发的全球2010-2100年间每十年的人口网格预测数据,按照预先设定的行政级别,提取不同年份的人口数据,再结合待测区域的实际情况的逐年总人口预测数据,得到符合实际的人口分布数据集D1,作为人口空间分布预测结果;
根据历史建设用地扩张空间数据,以预先设定的基准年2020年60%的土地利用现状数据为基础,基于热点阈值打分的方法,得到目标年份城市建设用地空间分布数据集D2,作为城市建设用地空间扩张预测结果;
从城市建设用地空间扩张预测结果中,提取目标年份城市建设用地扩张面积;并结合对应的人口空间分布预测结果,获取第i个基准网格的加权路网空间分布密度数据Roadni:
其中,Roadni代表第n目标年份的第i个基准网格的加权路网空间分布密度数据;其中,n为2025,2030,2035,2040,2050和2060年;
ULni为第n目标年份的第i个基准网格的城市建设用地扩张面积;UL2020i为2020年的第i个基准网格的建设用地扩张面积;Popni为第n目标年份的第i个基准网格的人口数据;Pop2020i为2020年的第i个基准网格的人口数据;Road2020i为2020年第i个基准网格的加权路网空间分布密度数据;
将通过上式,将获取的每一个基准网格的加权路网空间分布密度数据进行汇总,得到目标年份加权路网空间分布密度数据集D3;
所述新能源车存活概率获取模块中,根据历史年份的新能源车的注册量和保有量,得到新能源车中的乘用车的存活概率和新能源车中的商用车的存活概率,作为新能源车的车龄存活概率的具体过程为:
假设机动车分为新能源车和燃油车两大类;其中,新能源车包括:乘用车和商用车;燃油车包括:乘用车和商用车;
根据新能源车的历史年份2000-2020每年的注册量和保有量,分别拟合得到新能源车乘用车和商用车存活曲线,该曲线反映新能源车群体随着车龄增加其存活比例逐渐下降的规律:
新能源车中的乘用车的存活概率:
其中,为第k目标年份,车型j为乘用车,车龄为i的新能源车中的乘用车的存活概率;n0,j(k)1为第k目标年份车型j为乘用车的注册量;ni,j(k+1)1为车型j为乘用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新乘用车保有量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度;Tj表示新能源车中的车型j为乘用车的寿命;
新能源车中的商用车的存活概率:
其中,为第k目标年份,车型j为商用车,车龄为i的新能源车中的商用车的存活概率;n0,j(k)2为第k目标年份车型j为商用车的注册量;ni,j(k+1)2为车型j为商用车的车辆在第k+i目标年份仍然正常使用的车辆数,即新商用车保有量;
该存活概率符合两参数的Weibull分布:
其中,bj表示失效陡度,Tj1表示车型j为新能源车中的商用车的寿命;
所述燃油车保有量获取模块中,根据Gompertz预测模型和增长率预测模型,计算未来目标年份的机动车的保有量,同时,基于未来目标年份的新能源车的新车渗透率,并结合新能源车的车龄存活概率,得到目标年份的新能源车保有量,利用目标年份的机动车总量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量的具体过程为:
获取未来目标年份的机动车保有量的具体过程如下:
假设机动车的车型分为乘用车和商用车;根据Gompertz预测模型,计算未来目标年份的乘用车保有量:
其中,Vp为未来目标年份的乘用车的保有量,其单位为辆;K为饱和保有量;a为第一常数值,b为第二常数值,二者均为已知的常数值;g为人均GDP,其单位为元/人;
根据增长率预测模型,计算未来目标年份的商用车的保有量:
V=a′+b′×H (7)
其中,V为未来目标年份的商用车的保有量,其单位为万辆;a′、b′为常数值,根据2002-2010年商用货车保有量和货运量的线性关系拟合得到;H为货运量,其单位为万吨;
则未来目标年份的机动车保有量A:
A=Vp+V; (8)
未来目标年份的新能源车保有量的具体计算过程如下:
假设新能源车包括:乘用车和商用车;
基于基准年的新能源车车龄分布,根据未来每年的新能源车的新车渗透率,预测未来年份每年的新能源车新车注册量,结合车辆存活概率曲线,进而得到目标年份的乘用车与商用车的新能源车的保有量:
未来目标年份新能源车中的乘用车保有量NVPk:
NVPk=∑j∑iNVPi,j,k (10)
其中,NVPi,j,k为k年的乘用车新能源车车龄分布,表示k年车龄为i的乘用车新能源车保有量,NVPk为k年新能源车中的乘用车保有量,NSPj,k-i为第k-i年的新能源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年2020年的车辆车龄分布,即NVP2020,为k年、车型j为新能源乘用车、车龄为i的存活概率;
其中,
NSPj,k-i=SPk-i×MSPk-i (11)
其中,SPk-i为k-i年乘用车的总销售量,MSPk-i为第k-i年的新能源车的新车渗透率;
未来目标年份新能源车中的商用车保有量NVk:
NVk=∑j∑iNVi,j,k (13)
其中,NVi,j,k为k年的商用车新能源车车龄分布,表示k年车龄为i的商用车新能源车保有量,NVk为k年新能源车中的商用车保有量,NSj,k-i为第k-i年的商用车新能源车注册量,这是一个动态迭代的公式,当k-i等于2020时代表基准年2020年的车辆车龄分布,即NV2020,为k年j车型为商用车新能源车、车龄为i的存活概率;
NSj,k-i=Sk-i×MSk-i (14)
其中,Sk-i为k-i年商用车的总销售量,MSk-i为第k-i年的商用车新能源车的新车渗透率;
目标年份的新能源车保有量等于未来目标年份新能源车中的乘用车保有量NVPk加上未来目标年份新能源车中的商用车保有量NVk;
将未来目标年份的机动车保有量减去目标年份的新能源车保有量,得到未来目标年份的机动车的燃油车保有量:
具体地,未来目标年份的机动车保有量A和k年新能源车乘用车和商用车的保有量NVPk、NVk,得到未来目标年份的机动车的燃油车保有量FVk;
FVk=A-NVPk-NVk (15);
所述预测模块中,将未来目标年份的机动车的燃油车保有量分配至每个基准网格,基于每一个基准网格的加权路网空间分布密度数据和每个基准网格的燃油车保有量,以及燃油车数量与二氧化碳排放的函数关系,预测目标年份的道路交通二氧化碳排放空间网格数据,实现对道路交通二氧化碳排放的预测的具体过程为:
针对待测区域在空间上划分的多个10km×10km的基准网格,将未来目标年份的机动车的燃油车保有量分配至每个基准网格,得到未来目标年份的各个基准网格的机动车的燃油车保有量Vnj:
其中,Vnj为目标年份的第n个基准网格的机动车的燃油车保有量;′RLS2020j为基准年份的第n个基准网格的加权路网空间分布密度数据,RLS2020j为基准年份的n个基准网格的加权路网空间分布密度数据的总和;Vn为目标年份的机动车的燃油车的保有量,Vn=FVk;
基于燃油车数量与二氧化碳排放的函数关系:
其中,Cn为预测的第n个基准网格的目标年份的道路交通二氧化碳排放空间网格数据;a1为第一权重值;b1为第二权重值;c为常数;RLSn为目标年份的n个基准网格的加权路网空间分布密度数据的总和;RLS2020为目标年份的第n个基准网格的加权路网空间分布密度数据,即目标年份加权路网空间分布密度数据集D3;V2020为基准年份的第n个基准网格的机动车的燃油车的保有量;C2020为基准年份的二氧化碳排放量,记为已知值;
根据得到的Cn,获取目标年份的道路交通二氧化碳排放空间网格数据C:
C=n×Cn (18)
其中,n为针对待测区域在空间上所划分的基准网格的数量;
并将所述目标年份的道路交通二氧化碳排放空间网格数据C作为预测结果,再根据二氧化碳排放总量与预测结果的差异进行倍数调整,得到横坐标为时间、纵坐标为二氧化碳排放量的变化曲线,根据该变化曲线的变化率平滑得到逐年的二氧化碳排放量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110923903.7A CN114118509B (zh) | 2021-08-12 | 2021-08-12 | 一种道路交通二氧化碳排放的预测方法及预测系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110923903.7A CN114118509B (zh) | 2021-08-12 | 2021-08-12 | 一种道路交通二氧化碳排放的预测方法及预测系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114118509A CN114118509A (zh) | 2022-03-01 |
CN114118509B true CN114118509B (zh) | 2022-08-23 |
Family
ID=80440898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110923903.7A Active CN114118509B (zh) | 2021-08-12 | 2021-08-12 | 一种道路交通二氧化碳排放的预测方法及预测系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114118509B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115455681B (zh) * | 2022-09-02 | 2023-08-25 | 江苏城乡建设职业学院 | 一种面向多载运工具的通勤交通碳排放空间分布估算方法 |
CN117271992B (zh) * | 2023-09-19 | 2024-03-26 | 东莞市东莞通股份有限公司 | 一种基于大数据的城市道路碳排放监测管理系统 |
CN117094868B (zh) * | 2023-10-19 | 2024-02-02 | 交通运输部规划研究院 | 公路碳排放量评估方法、装置及可读存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140129789A (ko) * | 2013-04-30 | 2014-11-07 | 전자부품연구원 | 운송 물품의 이산화탄소 배출량 예측 방법 및 시스템 |
CN112365041A (zh) * | 2020-11-06 | 2021-02-12 | 国网冀北电力有限公司张家口供电公司 | 基于系统动力学的城市电动汽车保有量的分类预测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105868865A (zh) * | 2016-04-13 | 2016-08-17 | 华北电力大学 | 基于多元线性回归与比例替代法的电动汽车保有量预测方法 |
JP7064356B2 (ja) * | 2018-03-14 | 2022-05-10 | 株式会社日立製作所 | 将来状態推定装置および将来状態推定方法 |
CN110956343A (zh) * | 2020-01-02 | 2020-04-03 | 国家电网有限公司 | 一种基于碳排放约束的电动汽车保有量预测方法及系统 |
-
2021
- 2021-08-12 CN CN202110923903.7A patent/CN114118509B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140129789A (ko) * | 2013-04-30 | 2014-11-07 | 전자부품연구원 | 운송 물품의 이산화탄소 배출량 예측 방법 및 시스템 |
CN112365041A (zh) * | 2020-11-06 | 2021-02-12 | 国网冀北电力有限公司张家口供电公司 | 基于系统动力学的城市电动汽车保有量的分类预测方法 |
Non-Patent Citations (4)
Title |
---|
城镇景观格局对区域碳排放影响及其差别化管控研究;郭杰等;《中国人口·资源与环境》;20181011(第10期);全文 * |
基于Richards-BP神经网络的交通运输碳排放量组合预测模型;高洁等;《长安大学学报(自然科学版)》;20130715(第04期);全文 * |
基于存活概率的动态车龄分布模型;杨方等;《中国安全科学学报》;20050630(第06期);全文 * |
面向城市群交通规划的道路碳排测算与分析方法;陈娇杨等;《公路交通科技》;20180215(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114118509A (zh) | 2022-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114118509B (zh) | 一种道路交通二氧化碳排放的预测方法及预测系统 | |
CN111400662B (zh) | 一种考虑电动汽车充电需求的空间负荷预测方法 | |
CN110390349A (zh) | 基于XGBoost模型的公交车客流量预测建模方法 | |
CN107180274B (zh) | 一种电动汽车充电设施规划典型场景选取和优化方法 | |
CN105160428A (zh) | 高速公路电动汽车快速充电站的规划方法 | |
CN108805335B (zh) | 一种公共自行车调度方法 | |
CN114417614B (zh) | 一种中心城区机动车空间管控措施下碳排放减排测算方法 | |
CN113962450B (zh) | 一种核算城市道路交通逐日二氧化碳排放量的方法及其系统 | |
CN112446609B (zh) | 一种基于大数据的充电场站建站方法及系统 | |
CN110019569B (zh) | 一种获取城市轨道交通运营状态信息的方法 | |
CN109035768B (zh) | 一种出租车绕路行为的识别方法 | |
CN110956343A (zh) | 一种基于碳排放约束的电动汽车保有量预测方法及系统 | |
CN110889086B (zh) | 一种评估共享电动租赁汽车对城市汽车尾气排放影响的方法 | |
CN115660501A (zh) | 电动汽车充电负荷可调裕度评估方法 | |
CN112036624B (zh) | 一种基于区域内电动汽车充电负荷预测的电网调度方法 | |
CN112070259B (zh) | 一种预测城市空载出租车数量的方法 | |
CN111126878A (zh) | 一种基于生态指数的城市交通运行评价方法 | |
CN112632374B (zh) | 一种考虑定制公交的居民出行方式选择分析方法 | |
CN112949939B (zh) | 基于随机森林模型的出租车载客热点预测方法 | |
CN113962599A (zh) | 一种城市轨道交通网络化运营管理方法及系统 | |
CN112419707B (zh) | 基于gps数据空间匹配的车辆运行效率评价方法及系统 | |
CN114240180A (zh) | 一种基于混合蚁群算法的多企业共享单车联合调度方法 | |
CN109558978B (zh) | 基于出行距离的区域交通方式划分方法 | |
CN112348341A (zh) | 一种地面公交实时竞争指数的计算方法 | |
CN109658000B (zh) | 一种计算静态交通指数的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |