logo资料库

三峡库区碳排放预测模型.pdf

发布时间：2022-05-29 发布人：admin 分类：说明书 资料大小：0.56M 资料格式：pdf 举报 版权申诉
第1页 / 共6页
第2页 / 共6页
第3页 / 共6页
第4页 / 共6页
第5页 / 共6页
第6页 / 共6页
资料共6页，全文预览结束
    第 44 卷第 6 期 2017 年 12 月 应 用 科 技 Applied Science and Technology Vol.44 No.6 Dec. 2017 DOI: 10.11991/yykj.201610006 网络出版地址：http://kns.cnki.net/kcms/detail/23.1191.U.20170428.0831.018.html 三峡库区碳排放预测模型 鲁祖亮1,2，李林1，曹龙舟1 1. 重庆三峡学院 非线性科学与系统结构重点实验室, 重庆 404100 2. 天津财经大学 数学与经济研究中心, 天津 300222 摘 要：为减小三峡工程对下游地区生态环境的影响，对三峡库区二氧化碳的排放进行预测。主要研究三峡库区的二氧化 碳排放预测问题，利用庞特里亚金最大值原理和马尔可夫过程建立三峡库区的二氧化碳排放模型，再将碳排放问题转化为 一阶偏微分方程，然后利用四阶变步长龙格-库塔方法计算出2015~2021年三峡库区的整体二氧化碳排放量。根据这些数 据，得出能源使用量与碳的排放量有着直接关系，二氧化碳的大量排放很可能导致全球的温室效应，然而一味减少能源使 用量会影响经济的增长，如何在两者之间达到平衡就显得很重要，同时为三峡库区环境治理提供数据支持。 关键词：三峡库区；碳排放量；环境治理；碳排放预测模型；哈密顿函数；庞特里亚金最大值原理；马尔科夫过程；龙格-库塔方法 中图分类号：O242.2，O175.2 文献标志码：A 文章编号：1009−671X(2017)06−0083−06 A Carbon emission prediction model in Three Gorges reservoir area LU Zuliang1,2, LI Lin1, CAO Longzhou1 1. Key Laboratory for Nonlinear Science and System Structure, Chongqing Three Gorges University, Chongqing 404100, China 2. Research Center for Mathematics and Economics, Tianjin University of Finance and Economics, Tianjin 300222, China Abstract: In order to reduce the ecological environment effects by the Three Gorges Project in the lower reaches of Three Gorges areas, a carbon emission prediction model was adopted. This paper mainly studies the prediction of carbon dioxide emissions in the Three Gorges Reservoir Area. By using Pontryagin's maximum principle and the Markov process, the model of carbon dioxide emissions in the Three Gorges Reservoir was established. The carbon emissions problem was transformed into a first-order differential equation. Then, the four-step variable step Runge Kutta method was used to calculate the total carbon dioxide emissions in the Three Gorges Reservoir Area during 2015~2021. According to these data, it can be derived that the amount of energy use and carbon emissions has a direct relationship, and a lot of carbon dioxide emissions is likely to lead to global greenhouse effect, however, the reduction of energy consumption will influence the growth of economy. So how to achieve a balance between the amount of energy and carbon dioxide emissions is very important, and we also can provide data for the environmental governance in the Three Gorges Reservoir area. Keywords: Three Gorges reservoir area; carbon emission; environmental governance; carbon emissions prediction model; Hamiltonian function; pontryagin maximum principle; Markov chain; Runge-Kutta method 三峡工程的全面完成，不仅标志着三峡大坝工 程的防洪、发电、航运、补水等综合效益开始全面发 挥，也标志着三峡移民任务的全面完成，三峡库区正 在由搬得出、稳得住逐步向实现安稳致富的战略目 收稿日期：2016−10−18. 网络出版日期：2017−04−28. 基金项目：国家自然科学基金项目(11201510, 11171251)；中国博 士后科学基金项目(2015M580197, 2017T100155)；重庆 市高校科研创新团队项目(CXTDX201601035)；重庆市 科委项目(cstc2015jcyjA20001)；教育部春晖计划项目 (Z2015139). 作者简介：鲁祖亮(1980−)，男，教授，博士. 通信作者：鲁祖亮，E-mail：zulianglux@126.com. 标转变，三峡库区整体进入了后三峡时代。 后三峡时代是指三峡工程建成后对下游生态、 生产、生活等方方面面的综合影响[1]。为了减小该 影响，中国政府计划投资1 700亿用于三峡工程后 期的规划工作。我们在关注三峡工程给我国带来巨 大的经济和社会效益的同时，也要重视三峡库区的 环境问题。环境问题影响人们的生活，给人们带来 诸多不便。如果不能有效解决这个问题，不但将危 害三峡库区及周边地区人们的生命财产安全，而且 将阻碍长江中下游地区甚至全国的可持续发展。
    ·84· 应 用 科 技 第 44 卷 二氧化碳过度排放会造成温室效应[2]，但是如 果全面减少二氧化碳的排放就会阻碍经济的发展。 就重庆三峡库区(包括万州区、涪陵区、巫溪县等 22个区、县，占整个三峡库区面积的85%)而言，碳 的排放主要来源于2个部分，其中一部分是来自人 们的日常活动，另一部分来源于工业生产。若是为 减少二氧化碳的排放而关停一些大型的工业企业， 将会阻碍经济的平稳快速增长。因此，寻求一种对 能源需求不高的经济增长路径是非常有必要的，需 要各行各业专家、学者的共同努力。 本文主要研究三峡库区的二氧化碳排放预测问 题，将利用庞特里亚金最大值原理和马尔可夫过程 建立三峡库区的二氧化碳排放模型，将碳排放问题 转化为一阶偏微分方程，然后利用四阶变步长龙格- 库塔方法计算出2015~2021年三峡库区的整体二氧 化碳排放量[3]，对三峡库区的碳排放情况进行预测， 为三峡库区环境治理提供数据支持。 1 数学模型 模型主要采用类似于文献[4]的处理方法，选取 效用函数： 式中：ct为每个社会个体在时间t的消费量， 为最 低消费水平，σ为风险厌恶系数， 称为Stone- Geary效用函数。这一效用函数表明人们的消费只 有达到了最低消费 ，即保证了生存下去的条件以 后，自己的需求、欲望才会得到满足，效用值才会 为正。 1.1 能源和经济之间动态关系的模型 在最优经济增长模型中[4-5]，全社会是由许许多 多的个体组成的，这些社会个体可以是单个成年人 或是家庭。每个社会个体在一段的时间内消费，都 会得到一定的效用值。全社会的效用值就是社会中 每一个社会个体效用值的总和。社会计划者的目 标，就是使未来无穷时间段内，全社会效用值贴现到 当前值的最大，即 式中：Ct为整个社会全体的消费，n为人口增长率， N0为初始年份的人口数。 每一个社会个体拥有一定量的资本以及一单位 的劳动力。本文的生产函数采用柯布−道格拉斯生 产函数的形式[6]，社会产出不仅与资本以及劳动力 有关，还与所消耗的能源有关。具体的生产函数形 式如下 式中：A0为总要素生产率，其增长速度为v，Kt为物 质资本投入，Et为能源投入，Lt为劳动力投入，α为 资本产出的弹性系数，γ为劳动力产出的弹性系 ， 数。其中各要素的弹性系数的和 这表明该生产函数的类型为递增报酬型，即按现有 技术水平，用扩大生产规模的方式来增加总产出，是 非常可行的。 考虑到能源投入来自两个方面：一是从国外进 口，另一个是本国生产。能源上的消费采用如下形 式[7-8]： 式中：θt为能源的平均价格，st为本国生产的能源的 份额， 为本国能源价格， 为外国能源的价格。 能源强度是指一个国家或地区、部门或行业单 位产值一定时间内消耗的能源量。根据定义可以给 出能源强度的表达式： 将式(4)带入式（3）中，得到 资本存量积累满足如下关系： 结合式(5)和式(2)的目标函数，得到如下最优 控制问题的模型： 相应的约束条件为 式中：t为时间，Nt为t时刻社会个体数量， 个社会个体的效用函数，ct为t时刻的每个个体的 消费，ρ为贴现率。将 带入式(1)中得 为每 其中消费Ct为控制变量，Kt为状态变量。 下面求解此最优控制问题，将采用传统的庞特 里亚金最大值原理。首先，定义哈密顿函数： u(ct)=(ct¯c)11¯cu(ct)¯cmaxU=w10u(ct)Ntexp(t)dt(1)u()u(ct)=(ct¯c)11maxU=w10N0exp((n)t)(ct¯c)11dt(2)Yt=A0exp(vt)KtE1tLt;0<<1(3)+(1)+>1t=pd;tst+pa;t(1st)Pd;tPd;t=EY(4)Yt=(A0exp(vt))1Ktt1Lt;0<<1Kt=YtKttEtCt=(YtttYt)KtCt=(1tt)(A0exp(vt))1Ktt1LtKtCt(5)maxU=w10N0exp((n)t)(ct¯c)11dt˙Kt=(1tt)(A0exp(vt))1Ktt1LtKtCtH=N0exp((n)t)(ct¯c)11+˙Kt
    第 6 期 鲁祖亮，等：三峡库区碳排放预测模型 ·85· 式中λ为资本的影子价格。其经济意义为求解资源 最优利用时，即在解决如何使有限资源的总产出最 大的过程中，得出相应的极小值，其解就是对偶解， 极小值作为对资源的经济评价，表现为影子价格。 然后，将哈密顿函数H对控制变量Ct求偏导，并令 其等于零： 最后，将哈密顿函数H对状态变量Kt求偏导，并令 它等于影子价格λ对时间的导数的相反数，得 型，n=1,2,3,4分别代表煤、石油、天然气、沼气等化 石燃料。 定义转移概率矩阵为 式中Pij从状态i到状态j的一步转移概率，且Pij满 足： ， 根据马尔可夫链的理论，在k时刻的状态向量 应满足： 在式(6)两边对t求偏导得 将式(6)和(7)分别代入式(8)，经过化简得 从而得到消费增长率的表达式： 2 模型中参数的估计 2.1 经济系统中的参数 对于生产函数 先在两边同时取对数，得到 在这里，为了更加符合实际，选取了形式复杂一 中含有 些的效用函数以后，导致了经济增长率 了未知函数Ct，这就使得不能根据能源强度τ的数 据来直接得出经济增长率的值。所以，采用数值方 法求解上述常微分方程(9)，得出未知函数Ct在具 体节点上的值，再根据实际经济意义，得到经济增长 率的值为 1.2 马尔可夫模型 马尔可夫过程是在给定当前知识或信息的情况 下，只有当前的状态用来预测将来，过去（即当前以 前的历史状态）对于预测将来（即当前以后的未来状 态）是无关的。即：未来的状态只与当前状态有关， 而与过去无关[9]。所以，马尔可夫模型中关键的 2个要素为初始状态向量与转移概率矩阵。 定 义 初 始 状 态 向 量 为 ，当 n=3时为产业结构预测模型，n=1,2,3分别代表第 一、第二、第三产业；当n=4时为能源混合预测模 式 中 ： ， ， ， , 为残差项。 再由《中国统计年鉴2010》以及《新中国 50年统计资料汇编》中的数据，通过线性回归的方 法得到相应参数的具体数据：总要素生产率的初始 值A0=0.017 4，总要素生产率增长速率v=0.001 44， 资本产出弹性系数α=0.781 5，劳动力产出的弹性系 数γ=0.390 2。 2.2 其他参数 效用函数中最低消费水平常数 ，本文采用重 庆民政局2014年颁布的最低生活保障作为数据来 源，具体数字为每人每月400元。初值 根据 《重庆统计年鉴2015》农场和城镇居民消费加权 平均得到[10]，数值为13 604.616元。 有关其他参数的估计详见参考文献[4]，具体结 果为资本折旧率δ=9.6%，贴现率ρ=0.065，风险厌恶 系数σ= 2.5，能源的不变平均价格θ=1 751元/吨油 当量。 人口增长率是根据联合国人口展望的数据以及 利用Lagrange插值方法，通过MATLAB软件计算 而来[11]，劳动力数据是用人口数据乘以劳动参与率 得到[12]，具体数据见表1。 @Ht@Ct=N0exp((n)t)(Ct¯C)=0(6)˙=@H@Kt=[(1tt)(A0exp(vt))1t1Lt](7)(n)N0exp((n)t)(Ct¯C)+N0exp((n)t)()(Ct¯C)1˙Ct˙=0(8)˙Ct=1[(n)+(1tt)(A0exp(vt))1t1Lt](Ct¯C)(9)gc=1[(n)+(1tt)(A0exp(vt))1t1Lt](Ct¯C)Ctg=˙CtCtgt=CtCt1Ct1s(0)=(s(0)1;s(0)2;;s(0)n)P=0BBBBBBBBBBBBBBBBBBBB@P11P12:::P1nP21P22:::P2n::::::::::::Pn1Pn2:::Pnn1CCCCCCCCCCCCCCCCCCCCA0⩽Pij⩽1(i;j=1;2;;n)n∑j=1Pij=1(i=1;2;;n)s(k)=s(0)PkYt=A0exp(vt)KtE1tLt;Y′=a0+vt+K′+L′+Y′=ln(Y=E)a0=ln(A0)K′=ln(K=E)L′=ln(L)¯CC(0)
    ·86· 应 用 科 技 第 44 卷 表1 2014~2021年人口增长率和劳动参与率 年份 人口增长率/% 劳动参与率 首先，将P2对角化得 2014-2015年 2015-2016年 2016-2017年 2017-2018年 2018-2019年 2019-2020年 2020-2021年 0.300 0.272 0.245 0.218 0.190 0.163 0.092 0.554 0.551 0.549 0.547 0.544 0.542 0.532 2.3 转移概率矩阵 根据国家统计局的数据，可以得出能源混合转 移概率矩阵为 于是 产业结构的转移概率矩阵为 3 三峡库区的碳排放预测 根 据 《 重 庆 统 计 年鉴2015》 的 数 据 显 示 ， 2014年，三峡库区第一、第二、第三产业的份额g分 别为9.2%、50.1%、40.7%。各产业能源强度的值 τ分别为27.658 7、165.924 6、51.768 3吨油当量/百 万元，再根据历史数据，各产业能源强度下降的速率 分别为2.1%、5.2%、4.1%。因此，计算出三峡库区 分产业的能源强度变化情况分别为 于是可以得到从2014年算起的第t年的能源 强度为 式中g1、g2、g3分别为第一产业、第二产业、第三产 业的份额。 这样，矩阵的乘法就转化为多项式运算。第 t年的能源强度τt可以写成 将式(11)以及前面得到的各参数代入式(9) 中，得初值问题： 由于人口增长率以及劳动参与率的数据不是常 数，所以在求解初值问题式(12)的时候，采用分段 的方法。由于表达式的形式复杂，不能给出问题的 解析解，因此考虑求该问题的数值近似解。本文利 用MATLAB软件编写四阶龙格-库塔方法的程序来 求解问题(12)的数值解。 龙格-库塔法是一种在工程上应用广泛的高精 度单步算法。理论基础来源于泰勒公式和使用斜率 P1=0BBBBBBBBBBBBBBBBBBBB@0:9870:0160:0000:0000:0000:8110:0290:1600:2530:0000:7250:0210:0000:4960:0000:5041CCCCCCCCCCCCCCCCCCCCAP2=0BBBBBBBBBBBB@0:9770:0080:0150:0000:9770:0230:0000:0230:9771CCCCCCCCCCCCA1=27:6587e0:021t106tYuan2=165:9246e0:052t106tYuan3=51:7683e0:041t106tYuan=g1t1+g2t2+g3t3=(g1tg2tg3t)0BBBBBBBBBBBB@1231CCCCCCCCCCCCA=(g10g20g30)P2t0BBBBBBBBBBBB@1231CCCCCCCCCCCCA(10)P2=VV1=0BBBBBBBBBBBB@10:57740:210400:57740:691300:57740:69131CCCCCCCCCCCCA0BBBBBBBBBBBB@0:97700010000:9541CCCCCCCCCCCCA0BBBBBBBBBBBB@10:57740:210400:57740:691300:57740:69131CCCCCCCCCCCCA1Pt2=0BBBBBBBBBBBB@10:57740:210400:57740:691300:57740:69131CCCCCCCCCCCCA0BBBBBBBBBBBB@0:97700010000:9541CCCCCCCCCCCCAt0BBBBBBBBBBBB@10:57740:210400:57740:691300:57740:69131CCCCCCCCCCCCA1t=0:0920:977t27:6587e0:021t+(0:060:977t+0:5+0:0610:954t)165:9246e0:052t+(0:0320:977t+0:50:0610:954t)51:7683e0:041t(11)8>>>>>>>>><>>>>>>>>>:˙Ct=1[(n)+(1tt)(A0exp(vt))1t1Lt](Ct¯C)C(0)=2112:7(12)
    第 6 期 鲁祖亮，等：三峡库区碳排放预测模型 ·87· 表4 能源强度及其下降速率 能源强度/（吨油当量·元−1） 能源强度下降速率 101.525 1×10−6 96.576 1×10−6 91.880 5×10−6 87.424 1×10−6 83.193 6×10−6 79.176 8×10−6 75.361 8×10−6 71.738 0×10−6 — 0.048 746 566 0.048 620 725 0.048 502 130 0.048 390 547 0.048 282 560 0.048 183 306 0.048 085 369 年份 2014 2015 2016 2017 2018 2019 2020 2021 再由式(4)得到 于是根据表3、4的数据，得出了2014~2021年 的能源使用量E，数据如表5。 表5 能源使用量 能源使用量/百万吨油当量 55.544 4 54.200 1 52.901 4 51.649 2 50.444 1 49.276 1 48.120 1 年份 2015 2016 2017 2018 2019 2020 2021 最后，根据能源使用量，求出碳排放的值。假定 碳排放系数均为恒定，其中煤、石油、天然气、非化 石能源的排放系数分别为1.005 2，0.753，0.617 3， 0。同时，根据《重庆统计年鉴2015》，煤、石油、天 然气、非化石能源的份额gi分别为0.586，0.134， 0.014 1，0.139。且份额的转化过程由转移概率矩阵 gi给出，于是，得到下面碳排放量P1的计算公式： 近似表达微分，它在积分区间多预计算出几个点的 斜率，然后进行加权平均，用做下一点计算的依据， 从而构造出精度更高的数值计算方法。 对于常微分方程初值问题： 其对应的四阶变步长龙格-库塔公式为[13] 为了求解人口增长率和劳动参与率，采用分段 的方法，其数值结果表2。 表2 消费及其增长率 年份 2015年 2016年 2017年 2018年 2019年 2020年 2021年 消费/亿元 2 166.1 2 221.7 2 279.0 2 338.2 2 399.5 2 462.6 2 526.3 消费增长率 0.025 275 715 0.025 668 252 0.025 791 061 0.025 976 305 0.026 216 748 0.026 297 145 0.025 866 970 根据表2中的消费增长率以及消费增长率和经 济 增 长 率 相 一 致 的 假 设 ， 结 合 《 重 庆 统 计 年鉴 2015》中三峡库区的GDP数据，推算出2014~2021 年三峡库区的社会总产出的数据如表3。 表3 社会产出 社会产出/亿元 5 609.57 5 751.36 5 898.98 6 051.12 6 208.31 6 371.07 6 538.61 6 707.75 年份 2014 2015 2016 2017 2018 2019 2020 2021 接下来考虑能源强度，根据式(11)同样算出 2014~2021年的能源强度数据及其下降速率，结果 如表4。 根据式(13)以及表5，得出碳排放量的结果如 表6所示。 8>><>>:y′=f(x;y)y(x0)=y08>>>>>>>>>>>>>>>>>>>>>><>>>>>>>>>>>>>>>>>>>>>>:yn+1=yn+h6(K1+2K2+2K3+K4)K1=f(xn;yn)K2=f(xn+h2;yn+h2K1)K3=f(xn+h2;yn+h2K2)K4=f(xn+h;yn+hK3)E=Ye=E(g1(t)g2(t)g3(t)g4(t))0BBBBBBBBBBBBBBBBBBBB@1:00520:7530:617301CCCCCCCCCCCCCCCCCCCCA=E(g1(0)g2(0)g3(0)g4(0))Pt10BBBBBBBBBBBBBBBBBBBB@1:00520:7530:617301CCCCCCCCCCCCCCCCCCCCA(13)
    ·88· 应 用 科 技 第 44 卷 表6 碳排放量 碳排放量/t 45.646 4×106 38.888 6×106 35.088 7×106 32.376 6×106 30.242 8×106 28.464 7×106 26.917 8×106 年份 2015 2016 2017 2018 2019 2020 2021 根据表5、6中的数据，可以绘制出能源使用量 和碳排放量曲线，如图1所示。 图1 能源使用量和碳排放量的趋势 从图1可以看到，三峡库区的能源使用量和碳 排放量的走势大体上一致，都呈现减少的趋势，但是 能源使用量减少趋势比较平缓，碳排放量在2016年 之前的减少相对比较剧烈，2016年之后比较平缓， 总体来说三峡库区能源使用量和碳排放量在今后几 年都会持续减少。 4 结论 1)应用庞特里亚金最大值原理和马尔可夫过程 研究三峡库区的碳排放问题，建立了污染问题的最 优控制问题模型。 2)利用经典的龙格-库塔方法求解实际的三峡 库区污染问题，将对三峡库区的环境保护研究发挥 越来越重要的作用。 3)根据数值结果，得出能源使用量与碳的排放 量有着直接关系，如果大量排放二氧化碳很可能导 致全球的温室效应，然而刻意减少能源使用量就会 本文引用格式： 影响三峡库区经济的增长。 本文计算出的2015～2021年的碳排放量，可以 用来预测三峡库区未来几年的碳排放的趋势，合理 调整各种能源的使用量，这将对三峡库区的环境治 理起到非常重要的作用。 参考文献： [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] empirical application[J]. 孙诗一. 四阶龙格−库塔法的原理及其应用[J]. 徐候君. 后三峡时代的库区移民新城规划建设探索——以 开县为例[J]. 西部人居环境学刊, 2016, 31(4): 79–83. 张申宁. 温室效应与生物相互作用[J]. 广东化工, 2016, 43(5): 117–118. 冯建强, 数学学习与研究, 2017(17): 3–5. 朱永彬, 王铮, 庞丽. 基于经济模拟的中国能源消费与碳 排放高峰预测[J]. 地理学报, 2009, 64(8): 935–944. HIGHFILL J, MCASEY M. An optimal control problem in economics[J]. International journal of mathematics & mathematical sciences, 1991, 14(3): 537–544. SEOK M Y, SONN Y H. Productive energy consumption and economic growth: An endogenous growth model and its Resource and energy economics, 1996, 18(2): 189–200. JOHN A, PECCHENINO R. An overlapping generation model of growth and the environment[J]. Economis journal, 1994, 104(427): 1393–1410. 舒元, 谢识予, 孔爱国. 现在经济增长模型[M]. 上海: 复 旦大学出版社, 1998. 林元烈. 应用随机过程[M]. 北京: 清华大学出版社, 2002. 张富民, 童泽圣. 重庆统计年鉴2015[M]. 北京: 中国统 计出版社, 2015. 刘卫国. MATLAB程序设计与应用[M]. 北京: 清华大 学出版社, 2002. 王金营, 蔺丽莉. 中国人口劳动参与率与未来劳动力供 给分析[J]. 人口学刊, 2006(4): 19–24. 李庆扬, 王能超, 易大义. 数值分析[M]. 北京: 清华大学 出版社, 2008. [11] [12] [13] 鲁祖亮, 李林, 曹龙舟. 三峡库区碳排放预测模型[J]. 应用科技, 2017, 44(6): 83–88. LU Zuliang, LI Lin, CAO Longzhou. A Carbon emission prediction model in Three Gorges reservoir area[J]. Applied science and technology, 2017, 44(6): 83–88. 6055504540353025202015201620172018201920202021E/
    分享到：
    赞 收藏

    相关推荐

    资料库

    开发技术

    共收录19306份资料，累计14个分类，关注成员有19位，主要包括：Actionscript,Javascript,PHP,VB,Perl,Delphi,Python,C#,C,Java,C++,Web开发,硬件开发,其它

    热门标签

    Actionscript Javascript PHP VB Perl Delphi Python C# C Java C++ Web开发 硬件开发 其它

    最新资料