如何自己做加盟网站,比较出名的外贸公司有哪些,网站开发设计文档模板,咖啡店网站首页怎么做二氧化碳排放及时空分布测度
整体求解过程概述(摘要) 面临全球气候变化的巨大挑战#xff0c;我国积极响应《巴黎协定》的号召#xff0c;提出“2030年前碳达峰#xff0c;2060 年前实现碳中和”的碳排放发展目标#xff0c;并将碳中和相关工作作为 2021 年的重点任务之一…二氧化碳排放及时空分布测度
整体求解过程概述(摘要) 面临全球气候变化的巨大挑战我国积极响应《巴黎协定》的号召提出“2030年前碳达峰2060 年前实现碳中和”的碳排放发展目标并将碳中和相关工作作为 2021 年的重点任务之一。目前我国主流的碳排放核算方法是基于产品、企业和项目的自下而上的方法成本较大且效率较低而且县级统计数据严重缺失难以进行精细化的管理和政策指导。 夜光灯数据被广泛应用在社会经济指标的测度中并表现出良好的解释能力。本文结合夜光灯和社会经济数据利用提升回归树(BRT)建立其与碳排放的关系并采用贝叶斯法优化(BO)模型超参数。本文利用 2000-2017 年中国省级数据进行有监督学习并通过 K 折交叉验证法对模型的表现结果进行评价结果表明BO-BRT 具有较高的拟合优度和泛化能力。 基于 BO-BRT 模型我们对 2000-2017 年中国县级碳排放数据进行测度并结合标准差椭圆对中国碳排放空间分布进行分析结果发现碳排放主轴逐渐向西旋转表明中国发展过程内蒙古、新疆等西部地区碳排放量逐渐增加。 本文结合夜光灯和社会经济数据建立的我国县域碳排放量的估计模型提供了一种成本低、精度高的小尺度碳排放核算方法对我国县域碳排放测度、碳减排措施的落实乃至其经济发展都有着重要指导意义同时针对中国全域碳排放分析的结果对国家加强顶层设计、统筹协调各个地区碳排放具有重大实践意义。
问题分析 在国内产业向低碳转型的过程中宏观政府部门如何制定减排措施以及微观企业如何执行落实、我国区域经济结构如何协同发展都是各部门在未来计划里值得深入探讨的问题而科学客观的碳排放数据则是政策实践强有力的理论支撑。科学测度碳排放不仅可以对人类文明产业发展起到警示告诫作用还能刺激各行业对于绿色和新清洁能源的技术创新与研究使用推动促进中国作为世界第二大经济体进行绿色转型的模范作用提早达到“碳高峰”对于人类社会发展与生态平衡二者之间的交锋也具有积极的示范意义。 在此上述背景之下寻找一种低成本又高精度的二氧化碳排放强度测算方法对于研究我国高速增长的经济和二氧化碳排放量之间的关系显得尤为重要。本文通过收集国家统计局数据库中 2000-2017 年各省的年度 GDP、人均 GDP、人口总数、城镇化率等面板数据以及中国碳核算数据库中估算的 1997 - 2017 年中国2735 个县域的 CO2 排放量数据结合并且从经过统一后的 DMSP/OLS 和 NPP/VIIRS卫星图像中得到了具有高质量并且时间跨度范围广度的稳定夜间灯光数据并据此建立中国小地域尺度下的二氧化碳排放量核算模型提出相应规划与建议。本文提出的结果有助于弥补现有的数据差距并有助于制定减少中国二氧化碳排放的战略。除了应用在碳排放领域具有借鉴价值而且对于今后人口分布、GDP 预测和污染物估计等社会指标测度也提供着重要的参考意义。
模型的建立与求解整体论文缩略图 全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可
(代码和文档not free)
import numpy as np
import pandas as pd
import statsmodels.api as sm
import matplotlib.pyplot as plt
from statsmodels.graphics.tsaplots import plot_acf,plot_pacf
import itertools
from sklearn.metrics import r2_score as rs
import warningswarnings.filterwarnings(ignore)#忽略输出警告
plt.rcParams[font.sans-serif][SimHei]#用来正常显示中文标签
plt.rcParams[axes.unicode_minus]False#用来显示负号dfpd.read_excel(E:\or.xlsx#指定Month列作为索引列#天然气CO2排放量
NGEdf1[Natural Gas Electric Power Sector CO2 Emissions]
NGE.head()
#折线图
fig, axplt.subplots(figsize(15,15))
NGE.plot(axax,fontsize15)
ax.set_title(天然气碳排放,fontsize25)
ax.set_xlabel(时间月,fontsize25)
ax.set_ylabel(碳排放量(百万总吨),fontsize25)
ax.legend(locbest,fontsize15)
ax.grid()
#分解时序
#STLSeasonal and Trend decomposition using Loess是一个非常通用和稳健强硬的分解时间序列的方法
import statsmodels.api as sm
#decompostiontsa.STL(NGE).fit()报错这里前面加上索引sm
decompostionsm.tsa.STL(NGE).fit()#statsmodels.tsa.api:时间序列模型和方法
decompostion.plot()
#趋势效益
trenddecompostion.trend
#季节效应
seasonaldecompostion.seasonal
#随机效应
residualdecompostion.resid
#平稳性检验
#自定义函数用于ADF检查平稳性
from statsmodels.tsa.stattools import adfuller as ADF
def test_stationarity(timeseries,alpha):#alpha为检验选取的显著性水平adfADF(timeseries)padf[1]#p值critical_valueadf[4][5%]#在95%置信区间下的临界的ADF检验值test_statisticadf[0]#ADF统计量if palpha and test_statisticcritical_value:print(ADF平稳性检验结果在显著性水平%s下数据经检验平稳%alpha)return Trueelse:print(ADF平稳性检验结果在显著性水平%s下数据经检验不平稳%alpha)return False
#将数据化为平稳数据
#一阶差分
NGE_diff1NGE.diff(1)
#十二步差分
NGE_seasonalNGE_diff1.diff(12)#非平稳序列经过d阶常差分和D阶季节差分变为平稳时间序列
print(NGE_seasonal)
#十二步季节差分平稳性检验结果
test_stationarity(NGE_seasonal.dropna(),1e-3)#使用dropna()去除NaN值全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可
