一种居住建筑燃气碳排放量预测方法与流程

[0001]
本发明涉及二氧化碳预测技术领域，具体涉及一种居住建筑燃气碳排放量预测方法。


背景技术：

[0002]
气候变暖问题已成为全球关注的共同问题，低碳发展作为一种生态环保的发展模式成为社会普遍的共识。城市是碳排放的主体，我国约有85％碳排放来源于城市，根据相关统计数据，中国居住能耗占城市建筑部门能耗的16％-25％。住区作为居住的基本单元，其整体碳排放水平与住区结构、住宅形态密切相关。
[0003]
目前建筑燃气碳排放的核算往往依赖于能耗模拟软件进行计算，能耗模拟结果很难应用与城市规划管理，同时住区的实际情况千差万别，基于能耗软件的统一模拟往往很难准确核算不同住区的燃气碳排放。
[0004]
因此如何从规划领域入手实现居住建筑燃气碳排放的准确预测，成为亟待解决的问题。


技术实现要素：

[0005]
有鉴于此，本发明实施例提供了一种居住建筑燃气碳排放量预测方法，以解决现有技术中基于能耗软件的统一模拟对不同住区的燃气碳排放核算准确度不够高的问题。
[0006]
本发明实施例提供了一种居住建筑燃气碳排放量预测方法，包括：
[0007]
获取样本数据；样本数据包括居住建筑特征参数和与居住建筑特征参数对应的居住建筑单位面积碳排放值；
[0008]
根据特征重要性对样本数据进行评分排序，选取评分大于预设值对应的样本数据作为预测模型特征；
[0009]
根据预测模型特征建立预测模型，获得预测结果。
[0010]
可选地，居住建筑特征参数包括：居住建筑燃气排放数据、建筑面积数据和建筑形态特征；其中，建筑形态特征包括建筑面宽进深比、体型系数、四面窗墙比、平均窗墙比、梯户数、层高和层数。
[0011]
可选地，根据特征重要性对样本数据进行评分排序，选取评分大于预设值对应的样本数据作为预测模型特征包括：
[0012]
获取每个样本数据在决策树模型中每个子树中的平均贡献度；
[0013]
采用基尼指数评价平均贡献度，得到特征重要性评分。
[0014]
可选地，采用基尼指数评价所述平均贡献度，得到特征重要性评分包括：
[0015]
将特征重要性评分用f表示，将基尼指数用gi表示，假设有m个特征，计算每个样本数据中特征的基尼指数评分：
[0016][0017]
其中，k为类别总数，m、k、k均为自然数，p
mk
为第m个节点中第k个类别的特征所占比例；
[0018]
计算每个样本数据特征在第m个节点中的重要性：
[0019]
f
m
＝gi
m-gi
l-gi
r
[0020]
其中，gi
l
、gi
r
分别表示在第m个节点之后的两个新节点的基尼指数；
[0021]
计算每个样本数据特征在决策树的第i棵子树中的重要性：
[0022][0023]
其中，m为同一特征在决策树中出现的节点的集合，i为自然数；
[0024]
计算每个样本数据特征的得分：
[0025][0026]
其中，n为决策树的总数量。
[0027]
可选地，根据预测模型特征建立预测模型，获得预测结果包括：
[0028]
用极端梯度提升算法建立预测模型；
[0029]
通过预测模型特征训练预测模型；
[0030]
将待测建筑的特征参数代入预测模型，获得待测建筑的预测单位面积碳排放量；
[0031]
将预测单位面积碳排放量与待测建筑的居住建筑面积相乘，得到待测建筑的预测燃气碳排放总量。
[0032]
可选地，居住建筑单位面积碳排放值的计算如下：
[0033][0034]
其中，e为居住建筑单位面积碳排放量，p为统计一年居住建筑燃气排放量，d为燃气二氧化然排放因子，s为居住建筑的建筑面积。
[0035]
可选地，在根据特征重要性对样本数据进行评分排序之前，还包括：采用数据空值填充规则对样本数据进行预处理；
[0036]
对样本数据的各个特征进行归一化处理。
[0037]
可选地，对样本数据的各个特征进行归一化处理包括：
[0038][0039]
其中，x为样本数据的特征值，x
min
为特征值的最小值，x
max
为特征值的最大值，x
std
为特征值的归一化结果。
[0040]
可选地，通过预测模型特征训练预测模型包括：
[0041]
对数据集按照3:1的比例随机划分为训练集和测试集；
[0042]
按预设次数重复训练；
[0043]
选取均方误差小于等于样本数据平均值10％的预测模型为最优预测模型。
[0044]
可选地，选取均方误差小于等于样本数据平均值10％的预测模型为最优预测模型包括：
[0045]
获取每个预测模型的预测结果：
[0046][0047]
其中，为预测结果，j为自然数，f
j
(x
j
)表示第j棵树的模型函数，x
j
为预测模型特征；
[0048]
选取最优预测模型：
[0049][0050]
其中，y
t
表示经过训练的第t棵树的优化预测结果，t∈[1，n]，n为样本数据总数，ω(f
i
)表示正则项优化函数，为误差函数
[0051]
本发明实施例提供的一种居住建筑燃气碳排放量预测方法，根据居住建筑形态指标建立了居住建筑单位面积燃气碳排放预测模型，通过预测模型可以为规划中的居住小区、居住建筑提供燃气减碳建议，相较于传统线性回归预测模型、本模型具有更高的精度。本预测模型将规划居住小区建筑形态指标输入预测模型，得出规划小区居住建筑的单位面积燃气碳排放量，以碳排放量最小为目标，通过调整居住建筑形态参数，从而得出单位面积燃气碳排放量最小的居住建筑形态结果，从而确定小区实际建筑形态，实现了居住区规划时控制碳排放量的目标。
附图说明
[0052]
通过参考附图会更加清楚的理解本发明的特征和优点，附图是示意性的而不应理解为对本发明进行任何限制，在附图中：
[0053]
图1示出了本发明实施例中一种居住建筑燃气碳排放量预测方法的流程图；
[0054]
图2示出了本发明实施例中一种筛选后特征评分示意图；
[0055]
图3示出了本发明实施例中一种预测结果与实际结果的拟合曲线对比图；
[0056]
图4示出了本发明实施例中一种居住建筑燃气碳排放量预测终端的结构图。
具体实施方式
[0057]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0058]
本发明实施例提供了一种居住建筑燃气碳排放量预测方法，如图1所示，包括：
[0059]
步骤s10，获取样本数据。样本数据包括居住建筑特征参数和与居住建筑特征参数
对应的居住建筑单位面积碳排放值。
[0060]
在本实施例中，居住建筑特征参数包括：居住建筑燃气排放数据、建筑面积数据和建筑形态特征；其中，建筑形态特征包括建筑面宽进深比、体型系数、四面窗墙比、平均窗墙比、梯户数、层高、层数和经纬度坐标等。例如在南方城市的某小区与北方城市的某小区，两种小区的建筑形态存在明显差异，例如冬季供暖；而相同地域的不同小区之间则较为相似。通过居住建筑特征参数可以将不同地区的建筑进行分类，使居住建筑特征与年平均碳排放量互相关联，针对不同经纬度地区的不同小区，均能够取得较为精确的预测结果。
[0061]
在具体实施例中，居住建筑指每栋居住单元楼，获取居住建筑燃气排放数据、建筑面积数据、建筑形态特征参数以及与居住建筑特征参数对应的居住建筑单位面积碳排放值。
[0062]
居住建筑单位面积碳排放值的计算过程如下：
[0063][0064]
其中，e为居住建筑单位面积碳排放量，p为统计一年居住建筑燃气排放量，d为燃气二氧化然排放因子，s为居住建筑的建筑面积。
[0065]
步骤s20，根据特征重要性对样本数据进行评分排序，选取评分大于预设值对应的样本数据作为预测模型特征。
[0066]
在本实施例中，通过xgboost算法的feature_importance属性判定样本数据集的特征重要性，feature_importance属性会根据特征对模型的增益给出得分，根据特征得分对其进行排序。在具体实施例中，预设值为0.1，如图2所示，横坐标为特征评分值，纵坐标为特征名称，保留得分超过0.1的特征。
[0067]
步骤s30，根据预测模型特征建立预测模型，获得预测结果。
[0068]
在本实施例中，将步骤s20中经过筛选得出的特征按得分高低进行排序，选出最合适的一项或多项特征作为数据集和训练集，用xgboost算法建立对所述居住建筑单位面积碳排放预测模型，以利用所述模型进行居住建筑燃气碳排放预测。将预测结果与居住建筑面积相乘，得到居住建筑燃气碳排放总量。
[0069]
在具体实施例中，步骤s30包括：
[0070]
步骤s301，用极端梯度提升算法建立预测模型。
[0071]
步骤s302，通过预测模型特征训练预测模型。
[0072]
步骤s303，将待测建筑的特征参数代入预测模型，获得待测建筑的预测单位面积碳排放量。
[0073]
步骤s304，将预测单位面积碳排放量与待测建筑的居住建筑面积相乘，得到待测建筑的预测燃气碳排放总量。
[0074]
本发明实施例提供的一种居住建筑燃气碳排放量预测方法，根据居住建筑形态指标建立了居住建筑单位面积燃气碳排放预测模型，通过预测模型可以为规划中的居住小区、居住建筑提供燃气减碳建议，相较于传统线性回归预测模型、本模型具有更高的精度。本预测模型将规划居住小区建筑形态指标输入预测模型，得出规划小区居住建筑的单位面积燃气碳排放量，以碳排放量最小为目标，通过调整居住建筑形态参数，从而得出单位面积燃气碳排放量最小的居住建筑形态结果，从而确定小区实际建筑形态，实现了居住区规划
时控制碳排放量的目标。
[0075]
作为可选的实施方式，步骤s20包括：
[0076]
步骤s201，获取每个样本数据在决策树模型中每个子树中的平均贡献度。
[0077]
步骤s202，采用基尼指数评价平均贡献度，得到特征重要性评分。
[0078]
在本实施例中，特征重要性排序和选择的计算过程为xgboost以树模型作为基础学习器，决策树中每一个节点，通过特征条件分裂为两个部分，不断重复这一过程；判断每个特征在决策树中做了多大的贡献，特征越重要，对节点的纯度增加的效果越好，选出基尼指数大于0.1的特征作为最终训练模型的特征。
[0079]
作为可选的实施方式，步骤s202包括：
[0080]
步骤s2021，将特征重要性评分用f表示，将基尼指数用gi表示，假设有m个特征，计算每个样本数据中特征的基尼指数评分：
[0081][0082]
其中，k为类别总数，m、k、k均为自然数，p
mk
为第m个节点中第k个类别的特征所占比例。
[0083]
步骤s2022，计算每个样本数据特征在第m个节点中的重要性：
[0084]
f
m
＝gi
m-gi
l-gi
r
[0085]
其中，gi
l
、gi
r
分别表示在第m个节点之后的两个新节点的基尼指数。
[0086]
步骤s2023，计算每个样本数据特征在决策树的第i棵子树中的重要性：
[0087][0088]
其中，m为同一特征在决策树中出现的节点的集合，i为自然数。
[0089]
步骤s2024，计算每个样本数据特征的得分：
[0090][0091]
其中，n为决策树的总数量。
[0092]
作为可选的实施方式，在步骤s20之前，还包括：
[0093]
步骤s11，采用数据空值填充规则对样本数据进行预处理。
[0094]
在本实施例中，数据空值填充规则为，单位面积燃气碳排放数据如果为空值，首先用处于同一个小区相同层数住宅均值填充，若无小区信息，则用面宽进深比、体型系数、层数相同的住宅均值填充。因为周边小区的建筑具有相似性，因此用周边小区的特征参数均值对空值部分进行填充，能够提高预测结果的精度。
[0095]
步骤s12，对样本数据的各个特征进行归一化处理。
[0096]
在本实施例中，对样本数据的各个特征进行归一化处理通过以下公式计算：
[0097][0098]
其中，x为样本数据的特征值，x
min
为特征值的最小值，x
max
为特征值的最大值，x
std
为特征值的归一化结果。
[0099]
对样本数据的各个特征进行归一化处理的作用是剔除不同数据间不同量纲的影响。
[0100]
作为可选的实施方式，步骤s302包括：
[0101]
步骤s3021，对数据集按照3:1的比例随机划分为训练集和测试集。
[0102]
在本实施例中，xgboost预测模型参数如下：
[0103]
booster＝'gbtree'，计算模型选用树模型，解决模型预测中非线性问题；
[0104]
最大增量步长(max_delta_step)＝0(取值范围为：[0,∞])，设置每棵树的权重最大增量步长为0表示没有约束；
[0105]
最大深度(max_depth)＝1(取值范围为：[1,∞])树的最大深度为1；
[0106]
树的最小样本数(min_child_weight)＝4(取值范围为：[0,∞])，指节点中最小样本权重和，如果一个节点样本权重和小于4则拆分过程结束，在现行回归模型中，这个参数是指建立每个模型所需要的最小样本数。
[0107]
步骤s3022，按预设次数重复训练。
[0108]
在本实施例中，重复训练次数设置为十次。
[0109]
步骤s3023，选取均方误差小于等于样本数据平均值10％的预测模型为最优预测模型。
[0110]
在本实施例中，预测模型的计算过程如下：
[0111]
获取每个预测模型的预测结果：
[0112][0113]
其中，为预测结果，j为自然数，f
j
表示第j棵树的模型函数，x
j
为预测模型特征；
[0114]
选取最优预测模型：
[0115][0116]
其中，y
t
表示经过训练的第t棵树的优化预测结果，f
t
为第t轮的训练模型，t∈[1，n]，n为样本数据总数，ω(f
i
)表示正则项优化函数，为误差函数。
[0117]
选取最优预测模型的过程即为优化预测模型的过程，通过误差函数以及正则项优化函数优化预测模型，使预测模型的结果误差缩小到10％以内。
[0118]
计算训练集预测结果同测试集预测结果之间的均方误差(mse)，由于训练次数有限，因此设置了均方误差比样本数据平均值小于或等于10％的判断条件，从而选出最优预测模型。在具体实施例中，可以增加训练次数，同时降低均方误差的参考值，从而进一步提高模型的预测精度。
[0119]
以四川省金堂县小区燃气数据作为研究对象，具体步骤如下：
[0120]
步骤1：获取居住建筑燃气排放数据、建筑面积数据、建筑形态特征。
[0121]
步骤2：利用所述居住建筑燃气排放数据、建筑面积数据将所述居住建筑燃气排放数据转化为居住建筑单位面积碳排放数据。
[0122]
步骤3：将所述居住建筑单位面积碳排放数据、所述建筑形态特征进行预处理，得到样本数据集。
[0123]
步骤4：通过xgboost算法的feature_importance属性判定样本数据集的特征重要性，在jupyter notebook上搭建scikit-learn平台，使用xgboost算法中的model.feature_importances代码调用feature_importance函数输出特征重要性的分排序。
[0124]
步骤5：根据步骤4所得到的特征排序进行选择，保留得分大于0.1的特征作为应用于模型预测的特征，如图2所示。
[0125]
步骤6：将步骤5选择出的特征的数据集训练集，用xgboost算法建立对所述居住建筑单位面积碳排放预测模型，以利用所述模型进行居住建筑燃气碳排放预测。将预测结果与居住建筑面积相乘，得到居住建筑燃气碳排放总量。
[0126]
取1个实际样本作为预测模型，示例结果如下表：
[0127][0128][0129]
一个预测模型的结果拟合曲线如图3所示，实线y为实际结果，虚线y_pred为预测结果。根据表格结果显示，均方差mse最小为0.024，针对以上述一个实际样本的预测模型，选取均方差为0.024的预测模型为最优模型。当预测模型的均方差最小时，对应的最优预测模型的预测结果曲线与实际结果曲线的拟合度最高，预测模型的稳定性较好。
[0130]
模型参数根据数据样本进行自优化调整，本方法应用于其他数据时，应根据步骤6所述模型选择最优模型作为最终预测模型，其模型具体参数根据数据样本进行调整，即可应用于不同建筑进行碳排放量预测，并且预测拟合度较高，可靠性强。
[0131]
本发明实施例还提供了一种居住建筑燃气碳排放量预测终端，如图4所示，该居住建筑燃气碳排放量预测终端可以包括处理器41和存储器42，其中处理器41和存储器42可以通过总线或者其他方式连接，图4中以通过总线连接为例。
[0132]
处理器41可以为中央处理器(central processing unit，cpu)。处理器41还可以为其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路
(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。
[0133]
存储器42作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块。处理器41通过运行存储在存储器42中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及数据处理，即实现上述方法实施例中的居住建筑燃气碳排放量预测方法。
[0134]
存储器42可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器41所创建的数据等。此外，存储器42可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器42可选包括相对于处理器41远程设置的存储器，这些远程存储器可以通过网络连接至处理器41。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
[0135]
所述一个或者多个模块存储在所述存储器42中，当被所述处理器41执行时，执行如图1-3所示实施例中的居住建筑燃气碳排放量预测方法。
[0136]
上述居住建筑燃气碳排放量预测终端具体细节可以对应参阅图1至图3所示的实施例中对应的相关描述和效果进行理解，此处不再赘述。
[0137]
本领域技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述存储介质可为磁碟、光盘、只读存储记忆体(read-only memory，rom)、随机存储记忆体(random access memory，ram)、快闪存储器(flash memory)、硬盘(hard disk drive，缩写：hdd)或固态硬盘(solid-state drive，ssd)等；所述存储介质还可以包括上述种类的存储器的组合。
[0138]
虽然结合附图描述了本发明的实施例，但是本领域技术人员可以在不脱离本发明的精神和范围的情况下作出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。