一种多环芳烃排放的投入产出分析和结构分解方法及系统技术方案

本发明专利技术公开的一种多环芳烃排放的投入产出分析和结构分解方法，包括以下步骤：获取若干年份的投入产出表并设定行业分类；根据价格指数，将投入产出表转换为指定年份价格；确定各个年份各个行业的各多环芳烃排放系数；进而获取多环芳烃总量；依据多环芳烃总量得到一列排放强度矩阵；通过投入产出表获取消费和投入系数矩阵，进行排放连锁分析，得到基于消费端和基于投入端的连锁矩阵；获取后向联系值和前向联系值，确定关键排放行业；本发明专利技术用于确定消费端和投入端的多环芳烃排放量，能明确在产业链中占排放主导的关键行业；通过年份结构分解分析，确定各社会经济因素对不同年份间的多环芳烃排放变化量的贡献，可用于制定各行业的多环芳烃减排策略。多环芳烃减排策略。多环芳烃减排策略。

【技术实现步骤摘要】
一种多环芳烃排放的投入产出分析和结构分解方法及系统


[0001]本专利技术涉及环境经济学的研究领域，特别涉及一种多环芳烃排放的投入产出分析和结构分解方法及系统。

技术介绍

[0002]现有的投入产出分析方法和结构分解方法都是针对二氧化碳，二氧化硫，PM2.5等污染物的排放，且缺少同时基于消费端和投入端的针对多环芳烃的分析方法。作为燃烧的副产物，多环芳烃具有致癌性，但是却没有针对多环芳烃的基于消费端和投入端的投入产出分析方法，而且也缺少了分析社会经济因素对多环芳烃变化量的结构分解方法。因此，针对多环芳烃的关键排放行业，行业排放特征及流动模式都不明确。此外，社会经济因素对多环芳烃排放的影响也有待明确。

技术实现思路

[0003]本专利技术的目的在于克服现有技术的缺点与不足，提供一种多环芳烃排放的投入产出分析和结构分解方法及系统。
[0004]本专利技术的第一目的在于提供一种多环芳烃排放的投入产出分析和结构分解方法。
[0005]本专利技术的第二目的在于提供一种多环芳烃排放的投入产出分析和结构分解系统。
[0006]本专利技术的第一目的通过以下的技术方案实现：
[0007]一种多环芳烃排放的投入产出分析和结构分解方法，包括以下步骤：
[0008]通过数据获取和分类模块获取若干年份的投入产出表，并设定行业分类；
[0009]根据各年份的各个行业价格指数，将投入产出表通过双重通缩法转换为指定年份的价格；
[0010]根据技术分割法确定各个年份各个行业的各种燃料燃烧的多环芳烃排放系数、工业过程的多环芳烃排放系数以及农业残留物燃烧的多环芳烃排放系数；
[0011]通过各多环芳烃排放系数获取指定年份各个行业生产过程所排放的多环芳烃总量；
[0012]依据多环芳烃总量得到一列排放强度矩阵，将该一列排放强度矩阵转换为对角线矩阵，得到该指定年份多环芳烃的排放强度矩阵E；
[0013]通过投入产出表获取消费系数矩阵A和投入系数矩阵B，并通过消费系数矩阵A和单位矩阵I得到Leontief逆矩阵L,通过投入系数矩阵B和单位矩阵I得到Ghosh逆矩阵G；
[0014]将各个最终使用转换为对角矩阵，得到各个最终使用矩阵F；将各个初始投入转换为对角线矩阵，得到各个初始投入矩阵V；
[0015]根据排放强度矩阵E、Leontief逆矩阵L、Ghosh逆矩阵G，进行排放连锁分析，得到基于消费端的连锁矩阵和基于投入端的连锁矩阵；通过基于消费端的连锁矩阵获取各个行业的后向联系值，通过基于投入端的连锁矩阵获取各个行业的前向联系值，根据后向联系值和前向联系值确定关键排放行业；
[0016]进行矩阵相乘运算，得到基于消费端的多环芳烃排放矩阵和投入端的多环芳烃排放矩阵；
[0017]利用结构分解方法分析社会经济因素对于多环芳烃排放变化量的影响；
[0018]通过两级分解法得出各个社会经济因素对不同年份间的多环芳烃排放变化量的贡献。
[0019]进一步地，所述设定行业分类，具体为：分别设定农业、采掘业、轻工业、重工业、能源生产业、建筑业、交通业、其他行业。
[0020]进一步地，所述通过各多环芳烃排放系数获取指定年份各个行业生产过程所排放的多环芳烃总量，具体为：确定各个行业在指定年份的各种燃料消耗量、特定工业产品生产量、农业残留物的燃烧量，将其与对应的多环芳烃排放系数相乘，得到各个行业在生产过程中所排放的多环芳烃总量；燃料消耗量从统计年鉴中得到，由于燃料消耗的行业分类少于投入产出表的行业分类，所以需要根据生产燃料的行业对其他行业的投入比例，将燃料消耗的行业分类进行拆分。
[0021]进一步地，所述依据多环芳烃总量得到一列排放强度矩阵，将该一列排放强度矩阵转换为对角线矩阵，得到该指定年份多环芳烃的排放强度矩阵E，具体为：得到各个行业在生产过程中所排放的多环芳烃总量之后，将其除以投入产出表中各个行业的总产出值，得到一列排放强度矩阵；然后将其转换为对角线矩阵，即为该年份的多环芳烃排放强度矩阵E。
[0022]进一步地，所述通过投入产出表获取消费系数矩阵A和投入系数矩阵B，并通过消费系数矩阵A和单位矩阵I得到Leontief逆矩阵L,通过投入系数矩阵B和单位矩阵I得到Ghosh逆矩阵G，具体为：将每个投入产出表的中间投入除以每列的行业的总产出，即可得到消费系数矩阵A，将单位矩阵I减去消费系数矩阵，再进行逆矩阵化，即可得到Leontief逆矩阵L，L＝(I
‑
A)
‑1；将每个投入产出表的中间投入除以每行的行业的总产出，即可得到投入系数矩阵B，将单位矩阵I减去投入系数矩阵B，再进行逆矩阵化，即可得到Ghosh逆矩阵G，G＝(I
‑
B)
‑1。
[0023]进一步地，所述各个最终使用包括：农村居民消费，城镇居民消费，政府消费，固定资本形成，存货变动，出口；所述初始投入矩阵包括：劳动者报酬，生产税净额，固定资产折旧，营业盈余，进口。
[0024]进一步地，所述根据排放强度矩阵E、Leontief逆矩阵L、Ghosh逆矩阵G，进行排放连锁分析，得到基于消费端的连锁矩阵和基于投入端的连锁矩阵；通过基于消费端的连锁矩阵获取各个行业的后向联系值，通过基于投入端的连锁矩阵获取各个行业的前向联系值，根据后向联系值和前向联系值确定关键排放行业，具体为：根据排放强度矩阵，Leontief逆矩阵和Ghosh逆矩阵，进行排放连锁分析；在Matlab中进行矩阵相乘运算，基于消费端的排放连锁矩阵l为排放强度矩阵E乘以Leontief逆矩阵L，即l＝E*L，在消费端排放连锁矩阵中，求每列之和，将总值与列数相除得到平均值，再将每列之和与平均值相除，得到各个行业的后向联系值；在投入端排放连锁矩阵中，求每行之和，将总值与行数相除得到平均值，再将每行之和与平均值相除，得到各个行业的前向联系值；若某行业的后向与前向联系值都大于1，则该行业被确定为关键排放行业。
[0025]进一步地，所述进行矩阵相乘运算，得到基于消费端的多环芳烃排放矩阵和投入
端的多环芳烃排放矩阵；具体为：在Matlab中进行矩阵相乘运算，得到基于消费端和投入端的多环芳烃排放矩阵；其中，基于消费端的排放矩阵为Qc，是排放强度矩阵E乘以Leontief逆矩阵L乘以最终使用矩阵F得到，即Qc＝E*L*F；基于投入端的排放矩阵为Qi，是初始投入矩阵V乘以排放强度矩阵E乘以Ghosh逆矩阵G得到，即Qi＝V*G*E；在消费端的排放矩阵中，数值代表了每行的行业因每列的最终使用所排放的多环芳烃质量，则列之和代表了此行业的多环芳烃消费端排放量；而在投入端的排放矩阵中，数值代表了每列的行业因每行的初始投入所排放的多环芳烃质量，行之和代表了此行业的多环芳烃投入端排放量。
[0026]进一步地，所述利用结构分解方法分析社会经济因素对于多环芳烃排放变化量的影响；具体为：利用结构分解方法分析社会经济因素对于多环芳烃排放变化量的影响；将最终使用矩阵F分解成三个因素：最终使本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种多环芳烃排放的投入产出分析和结构分解方法，其特征在于，包括以下步骤：通过数据获取和分类模块获取若干年份的投入产出表，并设定行业分类；根据各年份的各个行业价格指数，将投入产出表通过双重通缩法转换为指定年份的价格；根据技术分割法确定各个年份各个行业的各种燃料燃烧的多环芳烃排放系数、工业过程的多环芳烃排放系数以及农业残留物燃烧的多环芳烃排放系数；通过各多环芳烃排放系数获取指定年份各个行业生产过程所排放的多环芳烃总量；依据多环芳烃总量得到一列排放强度矩阵，将该一列排放强度矩阵转换为对角线矩阵，得到该指定年份多环芳烃的排放强度矩阵E；通过投入产出表获取消费系数矩阵A和投入系数矩阵B，并通过消费系数矩阵A和单位矩阵I得到Leontief逆矩阵L,通过投入系数矩阵B和单位矩阵I得到Ghosh逆矩阵G；将各个最终使用转换为对角矩阵，得到各个最终使用矩阵F；将各个初始投入转换为对角线矩阵，得到各个初始投入矩阵V；根据排放强度矩阵E、Leontief逆矩阵L、Ghosh逆矩阵G，进行排放连锁分析，得到基于消费端的连锁矩阵和基于投入端的连锁矩阵；通过基于消费端的连锁矩阵获取各个行业的后向联系值，通过基于投入端的连锁矩阵获取各个行业的前向联系值，根据后向联系值和前向联系值确定关键排放行业；进行矩阵相乘运算，得到基于消费端的多环芳烃排放矩阵和投入端的多环芳烃排放矩阵；利用结构分解方法分析社会经济因素对于多环芳烃排放变化量的影响；通过两级分解法得出各个社会经济因素对不同年份间的多环芳烃排放变化量的贡献。2.根据权利要求1所述的一种多环芳烃排放的投入产出分析和结构分解方法，其特征在于，所述设定行业分类，具体为：分别设定农业、采掘业、轻工业、重工业、能源生产业、建筑业、交通业、其他行业。3.根据权利要求1所述的一种多环芳烃排放的投入产出分析和结构分解方法，其特征在于，所述通过各多环芳烃排放系数获取指定年份各个行业生产过程所排放的多环芳烃总量，具体为：确定各个行业在指定年份的各种燃料消耗量、特定工业产品生产量、农业残留物的燃烧量，将其与对应的多环芳烃排放系数相乘，得到各个行业在生产过程中所排放的多环芳烃总量；燃料消耗量从统计年鉴中得到，由于燃料消耗的行业分类少于投入产出表的行业分类，所以需要根据生产燃料的行业对其他行业的投入比例，将燃料消耗的行业分类进行拆分。4.根据权利要求1所述的一种多环芳烃排放的投入产出分析和结构分解方法，其特征在于，所述依据多环芳烃总量得到一列排放强度矩阵，将该一列排放强度矩阵转换为对角线矩阵，得到该指定年份多环芳烃的排放强度矩阵E，具体为：得到各个行业在生产过程中所排放的多环芳烃总量之后，将其除以投入产出表中各个行业的总产出值，得到一列排放强度矩阵；然后将其转换为对角线矩阵，即为该年份的多环芳烃排放强度矩阵E。5.根据权利要求1所述的一种多环芳烃排放的投入产出分析和结构分解方法，其特征在于，所述通过投入产出表获取消费系数矩阵A和投入系数矩阵B，并通过消费系数矩阵A和单位矩阵I得到Leontief逆矩阵L,通过投入系数矩阵B和单位矩阵I得到Ghosh逆矩阵G，具
体为：将每个投入产出表的中间投入除以每列的行业的总产出，即可得到消费系数矩阵A，将单位矩阵I减去消费系数矩阵，再进行逆矩阵化，即可得到Leontief逆矩阵L，L＝(I
‑
A)
‑1；将每个投入产出表的中间投入除以每行的行业的总产出，即可得到投入系数矩阵B，将单位矩阵I减去投入系数矩阵B，再进行逆矩阵化，即可得到Ghosh逆矩阵G，G＝(I
‑
B)
‑1。6.根据权利要求1所述的一种多环芳烃排放的投入产出分析和结构分解方法，其特征在于，所述各个最终使用包括：农村居民消费，城镇居民消费，政府消费，固定资本形成，存货变动，出口；所述初始投入矩阵包括：劳动者报酬，生产税净额，固定资产折旧，营业盈余，进口。7.根据权利要求1所述的一种多环芳烃排放的投入产出分析和结构分解方法，其特征在于，所述根据排放强度矩阵E、Leontief逆矩阵L、Ghosh逆矩阵G，进行排放连锁分析，得到基于消费端的连锁矩阵和基于投入端的连锁矩阵；通过基于消费端的连锁矩阵获取各个行业的后向联系值，通过基于投入端的连锁矩阵获取各个行业...

【专利技术属性】
技术研发人员：李瑞飞，张晋，董良，
申请(专利权)人：暨南大学，
类型：发明
国别省市：