1 背景
自18世纪中叶以来，由于快速城市化、人口增长和技术发展，导致一氧化二氮（N2O）、 甲烷（CH4）和二氧化碳（CO 2）等温室气体浓度急剧上升，引发了全球变暖、海平面上 升、极端天气以及环境污染等一系列问题，严重制约了社会、经济、生态的可持续发展， 威胁人类生存与健康[1]。
由于温室气体排放增加引起的全球变暖、极 端高温和热浪、龙卷风、飓风、干旱和洪水等自然灾害成为新常态，已经成为世界各国 政府和学术界关注的焦点[ 2]，因此，估算大气中CO2浓度是我们研究全球变暖等问题的 最可靠的方法，探讨CO2循环和碳源汇收支的变化规律是应对全球气候变化的关键所在。故本文针对巴西1960年-2014年固体燃料消耗产生的CO2排放量来进行分析和预测，针对特定的数据进行建模分析，最终得出相应的结论。
2 数据和方法说明
本文所运用到的数据是全球暖化数据集中的全球国家CO2排放情况表（分燃料状态）（年）其中的巴西的数据，得到数据后，对数据进行了相应的筛选，其数据展示如下：
表1 1960年-2014年巴西固体燃料消耗产生的CO2排放量原始数据
| SgnYear | Cntrnm | Region | IncomeGroup | Solid_CO2m | Liquid_CO2m |
| 1960 | 巴西 | 拉丁美洲 | 中等偏上 | 4968.79 | 39049.88 |
| 1961 | 巴西 | 拉丁美洲 | 中等偏上 | 4682.76 | 41503.11 |
| ... | ... | ... | ... | ... | ... |
| 2014 | 巴西 | 拉丁美洲 | 中等偏上 | 73666.36 | 339028 |
3理论
略
4 实证分析
巴西固体燃料消耗产生的CO2排放量描述性统计分析
首先展示原始数据（前6行），如下图，随后进行整体数据的描述性统计分析：
表3 整体数据描述性统计
| Solid_CO2m | Liquid_CO2m | |||||
| min | 4683 | 39050 | ||||
| 1st Qu | 9487 | 118811 | ||||
| median | 35750 | 150336 | ||||
| mean | 32212 | 161237 | ||||
| 3st Qu | 48522 | 229956 | ||||
| max | 73666 | 339029 | ||||
| SgnYear | Cntrnm | Region | IncomeGroup | |||
| Length | 55 | 55 | 55 | 55 | ||
| calss | character | character | character | character | ||
| mode | character | character | character | character | ||
从表3可以看出，对巴西固体和液体燃料消耗产生的CO2排放量以及其他数据进行了描述性统计，得到了最大最小值，均值以及1/4分位数和3/4分位数，其中前四个变量为非数值型变量。且下图4画出了1960年-2014年巴西固体燃料消耗产生的CO2排放量的时序图。
ARIMA模型的构建
进行ARIMA模型构建之前，要对时间序列数据纯随机性和平稳性检验。可以判断数据是否具有建模的价值以及是否适合ARIMA模型。下面对巴西固体燃料消耗产生的CO2排放量数据进行纯随机性检验和平稳性检验结果如下表4和表5：
表4 纯随机检验
| 滞后期数 | 卡方统计量 | P值 |
| 滞后6期P值 | 234.39 | 0.000 |
| 滞后12期P值 | 350.1 | 0.000 |
下面进行自动定阶的函数，计算得到模型应该采用ARIMA(2,1,2)，拟合得到模型系数：
表 7 模型定阶系数
| Coefficients: | |||||
| s.e. | ar1 | ar2 | ma1 | ma2 | drift |
| -0.1213 | -0.8560 | -0.1862 | 0.9513 | 1236.9922 | |
| 0.1035 | 0.1002 | 0.0863 | 0.1673 | 330.2231 | |
| Sigma^2=8135344: likelihood=-505.11 | |||||
| Aic=1022.23 AICc=1024.01 BIC=1034.16 | |||||
随后进行模型判断和误差的计算：
最后进行预测，预测3期，即未来3年巴西的巴西固体燃料消耗产生的CO2排放量，
5 结论
巴西1960年-2014年固体燃料消耗产生的CO2排放量来进行分析和预测，针对特定的数据进行建模分析，最终得出相应的结论。ARIMA模型的预测方面的还可行性，针对预测的结果，可以对政策调整和其他方面的策略判断做出相应的参考，在理论上具有一定的参考价值。
本文代码
dataset1<- read.xlsx("巴西不同燃料的排放量.xlsx", sheet = 1)
dataset1###首先展示数据前6行
head(dataset1,6)###随后对整体数据进行描述性添加分析
summary(dataset1)###画出1960年-2014年巴西固体燃料消耗产生的CO2排放量的时间序列图形dataset1$Solid_CO2Emission
HG_I<-ts(dataset1$Solid_CO2Emission,start=c(1960),frequency=1)
HG_I
plot(HG_I,type="o",pch=20,main="1960年-2014年巴西固体燃料消耗产生的CO2排放量时间序列图",xlab = "年份/Y",ylab="排放量",col = "green")#白噪声检验
for(i in 1:2) print(Box.test(HG_I,type = "Ljung-Box",lag=6*i))
###P值很少，很明显为非白噪声，可继续建模library(stats)
ndiffs(HG_I)
###结果显示为需要1阶差分
##但是个人看 2 阶才能平稳
diff.HG_I<-diff(HG_I,2)
plot(diff.HG_I,main='2阶差分图')ADF2<-adf.test(diff.HG_I) #1阶差分单位根检验
ADF2# 确定ARIMA模型中的p,q
# 这里有两种方法，一种是凭对知识点的理解通过ACF函数图和PACF函数图自行判断
# p,q的值另一种是通过软件的算法自动预测。
acf(diff.HG_I,main='差分后acf',lag.max = 12)
pacf(diff.HG_I,main='差分后pacf',lag.max = 12)###模型拟合
HG_I.fit<-auto.arima(HG_I)
HG_I.fit #模型预测
per_HG_I<-forecast(HG_I.fit,h=3)
per_HG_I
plot(per_HG_I)
巴西co2数据和BG
