利用ARIMA-BP预测中国2000年至2035年各省碳排放量数据

介绍

中国在减少碳排放方面发挥着至关重要的作用。中国政府表示,“中国将积极稳妥推进二氧化碳排放达峰碳中和”。从图1可以看出,中国碳排放总量自2000年以来一直快速上升,2023年碳排放总量将达到126亿吨,占全球碳排放量374亿吨的三分之一,成为全球碳排放量最大的国家。世界上最大的碳排放国。

湖南财政经济学院SCI:利用ARIMA-BP预测中国2000年至2035年各省碳排放量,附预测数据

图片

其次,碳排放存在显着的地区差异。2023年,碳排放最高的省份是山西,最低的是青海。华东地区是碳排放量最大的地区,约占全国比例的30%。2各省碳排放量存在显着差异,是实现碳排放市场均衡协调的重要途径。

控制碳排放总量、优化碳排放结构是实现“碳达峰和碳中和目标”的必要手段。因此,基于这些问题,本文结合ARIMA时间序列分析和BP神经网络算法对全国和省份的碳排放进行测算和预测,分析我国碳排放的空间演化特征,同时分解驱动因素。LMDI(对数平均Divisia指数)因素。研究结论为我国不同地区因地制宜制定科学合理的减排政策提供决策参考。

1、ARIMA-BP神经网络模型构建

首先利用ARIMA模型对数据进行预测,得到预测误差序列。其次,建立了一个神经元和三个神经元的两层隐含层BP神经网络。将误差数据按时间顺序输入到神经网络中,输入节点设置为4,输出节点设置为1。通过滚动窗口,上一周期的误差会继续传入神经网络,如下部分输入,模型将不断修正和预测。最后利用组合模型预测碳排放强度变化(图2、图3)。具体步骤如下:

图片

图2.组合预测模型示意图。

图片

图3.BP神经网络结构图。

2、碳排放预测与分析

图片

图4.2000-2021年各省平均碳排放量(万吨)。

图片

图5.2000年至2035年能源消费强度预测图。

图片

图6.2000-2021年各省平均碳排放量(万吨)。

3、中国碳排放空间分布特征

图片

图7.2000年至2035年能源消费强度预测图。

图片

图8.2000-2035年能源比重预测图。

图片

图9.2000-2035年中国碳排放总量变化趋势预测

图片

图10.2000-2035年人均GDP预测(单位:元/人)。

4、研究结论

(1)采用ARIMA-BP神经网络模型预测2000—2035年各省碳排放总量。总体来看,我国碳排放总量逐年增加,但随着产业结构的不断优化,增速碳排放量逐渐减少。从碳排放的具体结构变化来看,碳排放总量为“第二产业>居民生活>第三产业>第一产业”,其中第二产业和居民生活的碳增长速度较快,而第一产业和第三产业变化幅度较小。

(2)利用ArcGIS描述中国各省份碳排放的空间分布特征。各省份碳排放空间分布呈现出典型的“东部>中部>西部”和“北部>南部”的不平衡分布格局,反映了不同经济发展阶段碳排放的特征。根据标准椭圆差分析,2000年至2035年碳排放中心将向西北方向移动。

(3)基于新发展理念,将各省份碳排放划分为不同群体,分析群体效应的异质性。数字经济水平高、产业结构先进、新型优质生产力地区碳排放量相对较小,群体差异效应显着,说明产业结构和新业态在减排中的重要性。

(4)基于LMDI分解方法,分析了2000-2035年能源碳排放的驱动因素。我国能源消费强度为“第二产业>第三产业>第一产业”,碳排放强度的降低主要集中在第二产业,能源结构逐渐由煤炭为主向多元化能源结构转变;能源消费强度效应是推动碳排放持续增长的主要因素,人均碳排放量持续增长GDP和能源消费结构效应是抑制碳排放的主要因素。人口规模效应在第一年对碳排放增长有促进作用,但后期大多具有抑制作用,产业结构效应影响相对较小。

数据信息

中国2000年至2035年各省碳排放量数据

图片

数据格式:xlsx

数据容量:22.05KB

数据时间范围:2000-2035

### ARIMA-BP联合预测模型的实现与应用 #### 1. ARIMA-BP联合预测模型简介 ARIMA-BP联合预测模型结合了自回归积分滑动平均模型(ARIMA)和反向传播神经网络(BP Neural Network)。这种组合旨在发挥两种方法的优势:ARIMA擅长捕捉线性关系并处理非平稳时间序列,而BP神经网络则能更好地拟合复杂的非线性模式。通过先使用ARIMA对原始数据进行预处理,去除趋势性和季节成分,再利用BP神经网络进一步学习残差中的潜在规律,可以提高整体预测性能。 #### 2. 数据准备阶段 对于给定的时间序列数据集,在构建ARIMA-BP混合模型之前,通常需要完成如下几个步骤: - **数据清洗**:删除异常值、填补缺失值等操作。 - **探索性数据分析**:绘制时序图观察是否存在明显趋势或周期特性;计算ACF/PACF函数辅助判断可能存在的滞后效应。 - **检验平稳性**:采用ADF单位根测试验证原序列是否满足弱平稳条件[^3]。 ```matlab % ADF Test for Stationarity Check [h,pValue,stat,critValue,reg]=adftest(data,'model','ARD'); disp(['Is the series stationary? ', num2str(h)]); ``` #### 3. 构建ARIMA部分 一旦确认输入信号具备一定的稳定性,则可以根据样本统计特征选择合适的p,d,q参数来定义ARIMA(p,d,q),并通过极大似然估计法求解最优系数矩阵θ̂ 。值得注意的是,当面对具有较强记忆特性的长期依赖结构时,还可以考虑引入分数阶差分项f以增强描述能力。 ```matlab % Fit an ARIMA Model to Data Mdl = arima('Constant',0,'D',1,'Seasonality',12,... 'MALags',[1 12],'SMALags',[],'Variance',NaN); EstMdl = estimate(Mdl,data); residuals = infer(EstMdl,data); % Extract Residual Series ``` #### 4. 训练BP神经网络模块 接下来基于上述得到的残差序列作为新的训练目标y* ,设计一个多层感知机架构用于逼近剩余误差分布情况。具体来说,就是选取适当数隐藏节点数h以及激活函数形式σ() 来搭建前馈型全连接层,并借助梯度下降算法迭代更新权重w直至收敛至全局极小点附近位置停止优化过程。 ```python from sklearn.neural_network import MLPRegressor nn_model = MLPRegressor(hidden_layer_sizes=(5,), activation='relu', solver='adam', max_iter=2000) nn_model.fit(X_train, y_residuals) # X_train should be lagged residuals or other features ``` #### 5. 综合输出最终结果 最后一步便是将两个子系统的贡献相加起来形成总的预期值ŷ :一方面来自ARIMA部分所产生的基础估计值μ̂ ;另一方面则是由经过调校后的ANN所给出调整增Δ̂ 。如此这般便完成了整个流程的设计思路介绍。 ```matlab forecast_arima = forecast(EstMdl,numPeriods); final_forecast = forecast_arima + nn_model.predict(lagged_features); ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值