Industrial nitrogen metabolism greenhouse gas statistical method, statistical system, statistical device and storage medium

Abstract

The invention discloses an industrial nitrogen metabolism greenhouse gas statistical method, a statistical system, a device and a storage medium, which are used for counting the temperature increasing potential total value of the influence of greenhouse gas emission of an industrial cluster, and the method comprises the following steps: dividing the industrial cluster into a plurality of sub-emission systems; establishing a sub-system emission model corresponding to each sub-emission system; each sub-emission system calculates a temperature increasing potential value based on the sub-system emission model, and the step of calculating the temperature increasing potential value comprises the substeps of obtaining input material data, and calculating the gas emission amount of various greenhouse gases in the sub-system emission model according to the input material data; carrying out normalization processing on the gas emissions of the various greenhouse gases according to the temperature increasing potential coefficients of the corresponding greenhouse gases and carrying out summation so as to obtain temperature increasing potential values of the sub-emission systems; the temperature increasing potential values of all the sub-emission systems are integrated to obtain the temperature increasing potential total value, regional greenhouse gas emission is considered in a unified mode, and the influence on the temperature increasing trend is visually and rapidly obtained.

Description

Technical Field

[0001] The present invention relates to the field of greenhouse gas computer accounting technology, and in particular to an industrial nitrogen metabolism greenhouse gas statistics method and a statistics system, device, and storage medium. Background Art

[0002] Against the backdrop of rapid urbanization, human activities such as farmland expansion, industrial agglomeration, and fossil fuel combustion have disrupted the natural nitrogen cycle, causing problems such as non-point source nitrogen pollution in watersheds. At the same time, it has also caused a significant increase in greenhouse gas emissions such as nitrous oxide (N2O), carbon dioxide (CO2), and methane (CH4).

[0003] Currently, major international environmental organizations are paying more attention to CO2 emissions. CH4 and N2O are also receiving more and more attention because of their higher global warming potential values, which are 28 and 265 times that of CO2 respectively, and their significant impact on climate change.

[0004] Therefore, it is necessary to detect greenhouse gas emissions from the river basin where the industrial cluster is located, so as to evaluate the impact of the greenhouse effect. From the perspective of existing carbon accounting methods, they can be summarized into three categories: actual measurement method, emission factor method, and mass balance method. The actual measurement method usually directly obtains the greenhouse gas emission flux in the river basin through field sampling. However, the actual measurement method can only provide specific values, which makes it difficult to trace the source and identify the key emission areas and sources in the river basin from a macro level.

[0005] The emission factor method and mass balance method are mostly limited to a single industry (agriculture, animal husbandry, aquaculture, etc.) and a single process (agricultural planting process, livestock breeding process, sewage discharge process, etc.). They lack unified consideration of the multi-process emissions of multiple greenhouse gases such as N2O, CO2, and CH4 from an industrial perspective. The emission of a single greenhouse gas cannot accurately consider the impact of the region on the greenhouse effect. Summary of the invention

[0006] The present invention aims to solve at least one of the technical problems existing in the prior art. To this end, the present invention proposes an industrial nitrogen metabolism greenhouse gas statistical method and a statistical system, device, and storage medium, which uniformly consider the greenhouse gas emissions of industrial clusters in a region, and can intuitively and quickly derive the impact of industrial clusters in a region on the warming trend.

[0007] The industrial nitrogen metabolism greenhouse gas statistical method according to the first aspect of the present invention is used to calculate the total warming potential of greenhouse gases emitted in industrial clusters, wherein there are multiple greenhouse gases, and the industrial nitrogen metabolism greenhouse gas statistical method includes: dividing the industrial cluster into multiple sub-emission systems; establishing a sub-system emission model corresponding to each sub-emission system according to the differences of each sub-emission system, wherein the establishment of the sub-system emission model includes establishing a material flow framework based on the material flow analysis method according to the flow process of the input material in the sub-emission system, and establishing a material flow framework based on the emission factor method according to the flow process of different greenhouse gases in the material. The emission situation in the flow framework derives a relationship model between the input materials and the emissions of various greenhouse gases in the sub-emission system; each sub-emission system calculates the warming potential value based on its own sub-system emission model, and the calculation of the warming potential value includes: obtaining the input material data of the corresponding sub-emission system, and calculating the gas emissions of various greenhouse gases in the sub-system emission model according to the input material data; normalizing the gas emissions of various greenhouse gases according to the warming potential coefficients of the corresponding greenhouse gases and summing them up to obtain the warming potential value of the sub-emission system; integrating the warming potential values ​​of each sub-emission system to obtain the total warming potential value.

[0008] The industrial nitrogen metabolism greenhouse gas statistics method according to the embodiment of the present invention has at least the following beneficial effects:

[0009] The greenhouse impact statistical method of the present invention first divides the industrial cluster into multiple sub-emission systems. In each sub-emission system, a material flow framework is established based on the material flow analysis method. The material flow framework can characterize the existence of the input material in the sub-emission system, the flow direction after the chemical reaction and the generated greenhouse gas. Then, based on the emission factor method and the emission factors of various greenhouse gases, the emission of various greenhouse gases emitted by the input material is considered, so as to establish a subsystem emission model for each sub-emission system. Then, according to the acquired input material data, the gas emission of various greenhouse gases is calculated. Different greenhouse gases have different effects on the warming trend. Various greenhouse gases are normalized and summed to obtain the warming potential value of the sub-emission system, and the total warming potential value is integrated. The impact of the regional industrial cluster on the warming trend can be known based on the total warming potential value. This design uniformly considers the greenhouse gas emissions of the industrial clusters in the region, and can intuitively and quickly derive the impact of the industrial clusters in the region on the warming trend.

[0010] According to some embodiments of the present invention, establishing a subsystem emission model corresponding to each sub-emission system includes: dividing the subsystem emission model into a plurality of unit models according to different types of activities producing biochemical reactions in the sub-emission system based on a material flow analysis method, wherein the unit model characterizes the relationship between an input substance and at least one greenhouse gas emission;

[0011] The calculation of the gas emission of various greenhouse gases in the subsystem emission model based on the input material data includes: inputting the input material data into each unit model to calculate to obtain the unit emission of various greenhouse gases; adding the unit emission of the same type of greenhouse gas to obtain the gas emission of the greenhouse gas.

[0012] According to some embodiments of the present invention, the gas emissions of various greenhouse gases are normalized according to the warming potential coefficients of the corresponding greenhouse gases and summed up to obtain the warming potential value of the sub-emission system, which includes: taking the product of the gas emissions of the greenhouse gas and the warming potential coefficient of the greenhouse gas as the gas potential value of the greenhouse gas; and summing the gas potential values ​​of multiple greenhouse gases to obtain the warming potential value of the sub-emission system.

[0013] According to some embodiments of the present invention, there is a farmland system in the plurality of sub-emission systems, and the warming potential value of the farmland system is:

[0014] E 农田系统 =E 化肥施用 +E 粪便还田 +E 秸秆还田 +E 稻田种植 +E 能源消费1 +E 淋溶径流 +E 大气氮沉降 ;

[0015] Among them, E 农田系统 is the warming potential of the farmland system, E 化肥施用 The greenhouse gas potential value of nitrogen input into the soil for agricultural fertilizer application, E 粪便还田 E is the greenhouse gas potential value generated by the nitrogen input of livestock and poultry manure into the soil. 秸秆还田 E is the greenhouse gas potential value of nitrogen input into the soil when straw is returned to the field; 稻田种植 is the greenhouse gas potential value produced during rice cultivation, E 能源消费1 The greenhouse gas potential value of agricultural energy consumption; E 淋溶径流 is the gas potential value of greenhouse gases caused by leaching runoff, E 大气氮沉降 is the greenhouse gas potential of soils due to atmospheric nitrogen deposition;

[0016] E 化肥施用 =N 化肥 ×f cropland ×GWP1;

[0017] Among them, N 化肥 is the fertilizer nitrogen input, f croplad is the emission factor of N2O greenhouse gas in agricultural soil, GWP1 is the warming potential coefficient of N2O greenhouse gas;

[0018] E 粪便还田 =N 粪便 ×f cropland ×GWP1;

[0019]

[0020] Among them, i is different types of livestock and poultry, N livestock is the nitrogen excretion of livestock and poultry feces and urine, is the rate of livestock and poultry manure returned to the field, A i is the annual output or end-year inventory of the i-th type of livestock and poultry, N i is the annual average nitrogen excretion of feces and urine of the i-th livestock and poultry, T i is the i-th livestock and poultry breeding cycle, with 365 days per year;

[0021] E 秸秆还田 =N 秸秆 ×f cropland ×GWP1;

[0022]

[0023] Among them, j is different crop types, Str j is the straw yield of the jth crop, is the nitrogen content of the jth crop straw, is the rate of returning straw to the field, j is the crop yield, R str / see is the straw-grain ratio;

[0024] E 大气氮沉降 =∑N 化肥,粪便,秸秆 ×Rate dep ×f dep ×GWP1;

[0025] Among them, Rate dep is the observed volatilization coefficient of nitrogen input in the region, f dep Emission factors of N2O greenhouse gases from soils caused by atmospheric nitrogen deposition;

[0026] E 淋溶径流 =∑N 化肥,粪便,秸秆 ×f leaching ×GWP1;

[0027] Among them, f leaching N2O greenhouse gas emission factors from soils due to nitrogen leaching and runoff losses;

[0028] E 稻田种植 =∑EF d ×AD d ×GWP2;

[0029] Among them, EFd is the CH4 emission factor for rice fields of different types, AD d is the rice planting area corresponding to the emission factor; d is the rice field type, and GWP2 is the greenhouse warming potential coefficient of CH4.

[0030] According to some embodiments of the present invention, there is a livestock and poultry breeding system in the plurality of sub-emission systems, and the warming potential value of the livestock and poultry breeding system is:

[0031] E 畜禽养殖系统 =E 畜禽粪便 +E 牲畜肠道发酵 +E 能源消费2 ;

[0032] Among them, E 畜禽养殖系统 is the warming potential value of livestock and poultry breeding system, E 畜禽粪便 E is the warming potential of greenhouse gases emitted by livestock and poultry manure during livestock and poultry breeding. 牲畜肠道发酵 is the gas potential value of greenhouse gases produced by livestock intestinal fermentation during livestock and poultry breeding, E 能源消费2 The greenhouse gas potential value of energy consumption in livestock and poultry farming;

[0033] E 畜禽粪便 =∑EF i1 ×AD i ×GWP1+∑EF i2 ×AD i ×GWP2;

[0034] E 牲畜肠道发酵 =∑EF i3 ×AD i ×GWP2;

[0035] Among them, EF i1 EF is the emission factor of CH4 greenhouse gas from livestock and poultry manure; i2 EF is the emission factor of N2O greenhouse gas from livestock and poultry manure; i3 AD is the greenhouse gas emission factor of CH4 from enteric fermentation of livestock and poultry; i The number of livestock and poultry slaughtered for different types;

[0036] Among the plurality of said sub-discharge systems there is an aquaculture system, and the warming potential of the aquaculture system is:

[0037] E 水产养殖系统 =E 水产养殖反硝化 +E 能源消费3 ;

[0038] Among them, E 水产养殖系统 is the warming potential of the aquaculture system, E 水产养殖反硝化E is the greenhouse gas potential produced by denitrification in water bodies during nitrogen fertilizer application in aquaculture. 能源消费3 Greenhouse gas potentials for energy consumption in freshwater aquaculture;

[0039] E 水产养殖反硝化 =N Surplus ×f AQ ×GWP1;

[0040] Among them, N Surplus represents the surplus nitrogen in the aquaculture system of the water plant, f AQ is the N2O emission factor for the residual nitrogen in the aquaculture system;

[0041] N Surplus =AQIN bait +AQIN fer +AQIN wf +AQIN dep +AQIN sprout -AQOUT prod ;

[0042] Among them, AQIN bait Amount of nitrogen input into feed for aquaculture systems, AQIN fer AQIN is the nitrogen fertilizer input for feed ratio in aquaculture system. dep Indicates atmospheric nitrogen deposition; AQIN sprout AQOUT is the nitrogen input to the aquaculture system. prod Nitrogen exports to aquaculture systems, capture fisheries and wild fisheries;

[0043] AQIN bait =∑W s ×N s ÷f s,bait ;

[0044] AQIN fer =AQIN bait ×f fer ;

[0045] AQIN dep =Rate dep ×Area AQ ;

[0046] AQIN wf =∑WF s ×N s ;

[0047] AQIN sprout =Spruot s ×N s,sprQut ;

[0048] AQOUT prod =∑W s ×N s ;

[0049] Among them, s is different types of aquatic products, W s is the output of the sth aquatic product, N s is the nitrogen content of the sth aquatic product, f s,bait is the feed conversion coefficient of the sth aquatic product, f fer Indicates the ratio parameter of nitrogen fertilizer in aquatic product feed, Rate dep is the observed volatilization coefficient of nitrogen input in the area, Area AQ is the freshwater aquaculture area in the study area, WF i is the output of wild aquatic products of category s, N i,F is the nitrogen content of the sth aquatic product, N i,sprout is the nitrogen content of the sth aquatic product seedlings;

[0050] Among the multiple emission sub-systems, there are urban green space and forest land systems, and the warming potential value of the urban green space and forest land system is:

[0051] E 城市绿地与林地系统 =E 森林脱氮 +E 城市草地化肥施用 +E 能源消费4 ;

[0052] Among them, E 城市绿地与林地系统 is the warming potential of the urban green space and forest system, E 森林脱氮 Gas potential of greenhouse gases produced by forest denitrification, E 城市草地化肥施用 Greenhouse gas potential from soil fertilization for urban grasslands, E 能源消费4 Gas potentials for greenhouse gases produced for energy consumption in the forest sector;

[0053] E 城市草地化肥施用 =UG fer ×UG area ×f UG ×GWP1;

[0054] E 森林脱氮 =∑FR area ×f forest ×GWP1;

[0055] Among them, UG fer is the amount of nitrogen fertilizer applied per unit area of ​​urban green lawn, UG area is the area of ​​urban green lawn, f UG is the emission factor of soil N2O greenhouse gas, FR area is the area corresponding to different forest types, fforest is the forest N2O greenhouse gas emission factor corresponding to different forest types;

[0056] Among the plurality of said sub-discharge systems there is a sewage treatment system, and the warming potential value of the sewage treatment system is:

[0057] E 污水处理系统 =E 工业废水处理过程 +E 生活废水处理过程 ;

[0058] Among them, E 污水处理系统 is the warming potential of the sewage treatment system.

[0059] E 工业废水处理过程 is the greenhouse gas potential value generated during the treatment of industrial wastewater, E 生活废水处理过程 The greenhouse gas potential value generated during the treatment of domestic wastewater;

[0060] E 生活废水处理过程 =TN1×f TN,removed,1 ×EF f1 ×GWP1+COD1×f COD,removed,1 ×EF f2 ×GWP2;

[0061] E 工业废水处理过程 =TN2×f TN,removed,2 ×EF f1 ×GWP1+COD2×f COD,removed,2 ×EF f2 ×GWP2;

[0062]

[0063] Among them, P is the permanent population of the city, X TN is the daily per capita total nitrogen pollution intensity, X COD is the pollution intensity of chemical oxygen demand per capita per day, f collect is the collection rate of urban domestic sewage in each region, TN1 is the total nitrogen content of domestic wastewater collected through the urban domestic sewage pipe network and entering the domestic sewage treatment plant, COD1 is the chemical oxygen demand collected through the urban domestic sewage pipe network and entering the domestic sewage treatment plant, IP is the output of industrial products, Y TN Y is the pollution coefficient of producing unit industrial product TN, COD is the pollution coefficient of COD per unit of industrial product production, TN2 is the total nitrogen produced by industrial wastewater product production, COD2 is the chemical oxygen demand produced by industrial wastewater product production, and f TN,remobed,1 is the average TN removal rate of sewage treatment plants in each region of the basin, f COD,removed,1 is the average COD treatment rate of sewage treatment plants in various regions within the basin; TN,removed,2is the average treatment rate of TN in the terminal treatment of industrial wastewater, f COD,removed,2 The average treatment rate of COD in the terminal treatment of industrial wastewater, EF f1 is the emission factor of N2O greenhouse gas in the sewage treatment process, EF f2 is the emission factor of CH4 greenhouse gas in the sewage treatment process;

[0064] Among the plurality of said sub-emission systems, there is a solid waste disposal system, and the warming potential value of the solid waste disposal system is:

[0065] E 固废处置系统 =E 垃圾填埋 +E 垃圾焚烧 +E 能源消费5 ;

[0066] Among them, E 固废处置系统 is the greenhouse gas potential of the solid waste disposal system, E 垃圾填埋 is the greenhouse gas potential value generated during the landfill of domestic waste, E 垃圾焚烧 is the potential value of greenhouse gases produced during the incineration of domestic waste, E 能源消费5 The gas potential of greenhouse gases generated by the energy input of waste incineration;

[0067] E 垃圾填埋 =∑ g MSW g，填埋 ×DOC g ×MCF×DOC f ×F×a1×(1-R1)×(1-OX)×GWP2;

[0068] Among them, g is the type of different waste components, MSW g，填埋 is the landfill volume of different types of waste components, DOC i is the degradable organic carbon fraction of different garbage components g, MCF is the correction factor for CH4 greenhouse gas, DOC f is the proportion of decomposed degradable organic carbon, F is the content of CH4 greenhouse gas in landfill gas, R1 is the recovery rate of CH4 greenhouse gas, OX is the oxidation factor, and a1 is the relative molecular mass ratio of CH4 and C;

[0069] E 垃圾焚烧 =∑ g MSW g，焚烧 ×dm g ×FC g ×FCF I ×OF×a2+MSW g，焚烧 ×EF b ×GWP1;

[0070] Among them, MSW g，焚烧is the incineration amount of different types of garbage components, dm g is the dry matter content of different garbage components g, FC g is the fossil carbon component of different garbage components g, OF is the oxidation factor, EF is b is the emission factor of N2O greenhouse gas in the waste incineration process, a2 is the relative molecular mass ratio of CO2 and C;

[0071] Among the plurality of said sub-discharge systems, there is a surface water system, and the warming potential value of the surface water system is:

[0072] E 地表水系统 =E 工业废水-地表水排放过程 +E 生活废水-地表水排放过程 ; E 地表水系统 is the greenhouse gas potential of the surface water system,

[0073] E 工业废水-地表水排放过程 It is the greenhouse gas potential value of pollutants produced by industrial wastewater after direct discharge or discharge into surface water after treatment.

[0074] E 生活废水-地表水排放过程 It is the greenhouse gas potential value of pollutants produced by direct discharge of domestic wastewater or after treatment into surface water;

[0075] E 生活废水-地表水排放过程 =(TN3+TN4)×EF f3 ×GWP1+(COD3+COD4)×EF f4 ×GWP2;

[0076] E 工业废水-地表水排放过程 =TN5×EF f3 ×GWP1+COD5×EF f4 ×GWP2;

[0077]

[0078] Among them, TN3 is the total nitrogen content of domestic wastewater entering the surface water after being treated in the sewage treatment plant, COD3 is the chemical oxygen demand of domestic wastewater entering the surface water after being treated in the sewage treatment plant, TN4 is the total nitrogen content of domestic wastewater directly discharged into the surface water system without being collected by the sewage pipe network, COD4 is the chemical oxygen demand of domestic wastewater directly discharged into the surface water system without being collected by the sewage pipe network, TN5 is the total nitrogen content of industrial wastewater entering the surface water system after treatment, COD5 is the chemical oxygen demand of industrial wastewater entering the surface water system after treatment, EF f3 is the emission factor of N2O greenhouse gas from wastewater discharged into rivers and lakes, EF f4 It is the emission factor of CH4 greenhouse gas from off-site discharge of wastewater into rivers and lakes.

[0079] According to some embodiments of the present invention, E 能源消费1 、E 能源消费2 、E 能源消费3 、E 能源消费4 and E 能源消费5 Both are:

[0080]

[0081] in, is the gas potential value of greenhouse gas caused by the combustion of fossil fuels in any one of energy consumption 1, energy consumption 2, energy consumption 3, energy consumption 4, and energy consumption 5, k is the type of energy consumed, is the k-th energy consumption, NCV k is the net calorific value of the k-th energy source, is the emission factor corresponding to each of the greenhouse gases produced by the k-th energy consumption, The carbon oxidation rate factor corresponding to the various greenhouse gases produced by the k-th energy consumption, GWP (1，2，3) is the corresponding warming potential coefficient of N2O, CH4 and CO2 greenhouse gases.

[0082] According to some embodiments of the present invention, there is a human consumption system in the plurality of emission sub-systems, and the warming potential value of human consumption is:

[0083] E 人类消费系统 =E 居民家庭能源消费 +E 交通运输能源消费 ;

[0084] Among them, E 人类消费系统 is the greenhouse gas potential of human consumption system, E 居民家庭能源消费 The gas potential value of greenhouse gases generated by household energy consumption, E 交通运输能源消费 Gas potentials for greenhouse gases produced by transport energy consumption;

[0085] Among the plurality of emission sub-systems, there are industrial and energy systems, and the warming potential of the industrial and energy systems is:

[0086] E 工业与能源系统 =∑E 工业分行业能源消费 ;

[0087] E 工业与能源系统 is the total greenhouse gas emissions from industry and energy systems,

[0088] E 工业分行业能源消费 Greenhouse gas emissions from energy consumption by industrial sub-sector;

[0089] E 居民家庭能源消费 、E 交通运输能源消费 and E 工业分行业能源消费Both are:

[0090]

[0091] in, is the gas potential value of greenhouse gases caused by the combustion of any energy fossil fuel in household energy consumption, transportation energy consumption, and industrial sub-sector energy consumption, k is the type of energy consumed, is the k-th energy consumption, NCV k is the net calorific value of the k-th energy source, is the emission factor corresponding to each of the greenhouse gases produced by the k-th energy consumption, The carbon oxidation rate factor corresponding to the various greenhouse gases produced by the k-th energy consumption, GWP (1，2，3) is the corresponding warming potential coefficient of N2O, CH4 and CO2 greenhouse gases.

[0092] The statistical system according to the second aspect of the present invention is used to execute the industrial nitrogen metabolism greenhouse gas statistical method disclosed in any of the above embodiments.

[0093] According to the control device of the third aspect of the present invention, the control device includes a memory and a processor, the memory stores a computer program, and the processor implements an industrial nitrogen metabolism greenhouse gas statistics method disclosed in any of the above embodiments when executing the computer program.

[0094] According to the computer-readable storage medium of the fourth aspect of the present invention, the computer-readable storage medium stores a computer program, characterized in that when the computer program is executed by a processor, an industrial nitrogen metabolism greenhouse gas statistics method disclosed in any of the above embodiments is implemented.

[0095] Additional aspects and advantages of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0096] The above and / or additional aspects and advantages of the present invention will become apparent and easily understood from the description of the embodiments in conjunction with the following drawings, in which:

[0097] Figure 1 This is a flow chart of one embodiment of the industrial nitrogen metabolism greenhouse gas statistics method of the present invention. DETAILED DESCRIPTION

[0098] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0099] It should be noted that, although the functional modules are divided in the device schematic diagram and the logical order is shown in the flowchart, in some cases, the steps shown or described may be performed in a different order than the module division in the device or the order in the flowchart. The terms "first", "second", etc. in the specification, claims and the above drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence.

[0100] In the description of the present invention, “several” means one or more, “more” means more than two, “greater than”, “less than”, “exceed”, etc. are understood as not including the number itself, and “above”, “below”, “within”, etc. are understood as including the number itself.

[0101] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which this application belongs. The terms used herein are only for the purpose of describing the embodiments of this application and are not intended to limit this application.

[0102] like Figure 1 As shown, the industrial nitrogen metabolism greenhouse gas statistical method according to the first aspect of the present invention is used to calculate the total warming potential of greenhouse gases emitted in industrial clusters, wherein there are multiple greenhouse gases, and the industrial nitrogen metabolism greenhouse gas statistical method includes:

[0103] S110, divide the industrial cluster into multiple emission sub-systems;

[0104] S120. According to the differences among the sub-emission systems, a sub-system emission model corresponding to each sub-emission system is established, wherein the establishment of the sub-system emission model includes establishing a material flow framework based on the material flow analysis method according to the flow process of the input material in the sub-emission system, and deriving a relationship model between the input material in the sub-emission system and the emission of various greenhouse gases in the material flow framework based on the emission factor method according to the emission conditions of different greenhouse gases;

[0105] S130. Each emission subsystem calculates a warming potential value based on its own emission subsystem model. The calculation of the warming potential value includes:

[0106] S140, obtaining input material data of the corresponding sub-emission system, and calculating the gas emissions of various greenhouse gases in the sub-system emission model according to the input material data;

[0107] S150, normalizing the gas emissions of various greenhouse gases according to the warming potential coefficients of the corresponding greenhouse gases and summing them up to obtain the warming potential value of the sub-emission system;

[0108] S160. Integrate the warming potential values ​​of each emission sub-system to obtain a total warming potential value.

[0109] The greenhouse impact statistical method of the present invention first divides the industrial cluster into multiple sub-emission systems. In each sub-emission system, a material flow framework is established based on the material flow analysis method. The material flow framework can characterize the existence of the input material in the sub-emission system, the flow direction after the chemical reaction and the generated greenhouse gas. Then, based on the emission factor method and the emission factors of various greenhouse gases, the emission of various greenhouse gases emitted by the input material is considered, so as to establish a subsystem emission model for each sub-emission system. Then, according to the acquired input material data, the gas emission of various greenhouse gases is calculated. Different greenhouse gases have different effects on the warming trend. Various greenhouse gases are normalized and summed to obtain the warming potential value of the sub-emission system, and the total warming potential value is integrated. The impact of the regional industrial cluster on the warming trend can be known based on the total warming potential value. This design uniformly considers the greenhouse gas emissions of the industrial clusters in the region, and can intuitively and quickly derive the impact of the industrial clusters in the region on the warming trend.

[0110] In some embodiments of the present invention, establishing a subsystem emission model corresponding to each sub-emission system includes:

[0111] Based on the material flow analysis method, the subsystem emission model is divided into a plurality of unit models according to different types of activities producing biochemical reactions in the sub-emission system, wherein the unit model represents the relationship between the input material and at least one greenhouse gas emission;

[0112] Calculating the gas emissions of various greenhouse gases in the subsystem emission model according to the input material data includes:

[0113] Input material data are fed into each unit model calculation to obtain unit emissions of various greenhouse gases;

[0114] The unit emissions of the same type of greenhouse gas are added together to obtain the gas emissions of that greenhouse gas.

[0115] In each subsystem emission model, according to different emission activities, it can be divided into multiple unit models. Take the farmland system as an example. In the farmland system, the greenhouse gas emissions of different activities such as fertilizer application, manure return to the field, rice field planting, and even users burning energy during the planting process are different. Therefore, the activities that may emit greenhouse gases in the sub-emission system can be further refined. It should be noted that each unit model may only emit one greenhouse gas or multiple greenhouse gases. When the unit model only emits one greenhouse gas, the unit model only needs to calculate the emission of the greenhouse gas in combination with the emission factor of the greenhouse gas, and then combine the warming potential of the greenhouse gas to obtain the warming potential value of the sub-emission system. When the unit model emits multiple greenhouse gases, it is necessary to use the emission factors of different greenhouse gases to calculate the emissions of the multiple greenhouse gases respectively, and then use the emissions of different greenhouse gases to correct the warming potential values ​​of the corresponding greenhouse gases and then combine them to calculate the warming potential value of the sub-emission system.

[0116] Specifically, the emission factors of different greenhouse gases in the subsystem emission model vary based on the different sub-emission systems and can be found in the Greenhouse Gas Inventory Guidelines and Provincial Greenhouse Gas Inventory Compilation Guidelines issued by the IPCC (Intergovernmental Panel on Climate Change).

[0117] In some embodiments of the present invention, the gas emissions of various greenhouse gases are normalized according to the warming potential coefficients of the corresponding greenhouse gases and summed to obtain the warming potential value of the sub-emission system, including:

[0118] The product of the gas emission of the greenhouse gas and the warming potential coefficient of the greenhouse gas is taken as the gas potential value of the greenhouse gas;

[0119] The gas potential values ​​of multiple greenhouse gases are summed to obtain the warming potential value of the sub-emission system.

[0120] There are many types of greenhouse gases, usually N2O, CO2, and CH4. The emission of CO2 is generally used as a reference to determine the impact of the greenhouse effect. The global warming potential values ​​of CH4 and N2O are 28 and 265 times that of CO2, respectively. Therefore, the global warming potential value of CO2 can be used as a reference value. If the warming potential coefficient of CO2 is set to 1, the warming potential coefficient of CH4 is 28, and the warming potential coefficient of N2O is 265. These warming potential coefficients can also be obtained by referring to the assessment report of the United Nations Framework Convention on Climate Change (UFFCCC, 2013).

[0121] It should be noted that this design mainly analyzes the metabolic process of nitrogen-containing input materials in industrial clusters. N2O emissions have the same roots and sources as CO2 and CH4 emissions. In addition, for the establishment of subsystem emission models, the material flow analysis method establishes a material flow framework based on the flow process of input materials in the sub-emission system. The data here can be obtained from the regional provincial, municipal, and county statistical yearbooks. Different unit scale data can be selected based on the selected area. The amount of input materials can be obtained from the regional provincial or municipal comprehensive energy balance table and the energy consumption data of industrial sub-sectors above designated size.

[0122] A subsystem emission model corresponding to each sub-emission system is established, wherein the establishment of the subsystem emission model includes establishing a material flow framework based on the material flow analysis method according to the flow process of the input material in the sub-emission system, and deriving a relationship model between the input material in the sub-emission system and the emission of various greenhouse gases in the material flow framework based on the emission factor method.

[0123] In some embodiments of the present invention, there is a farmland system among the plurality of sub-emission systems, and the warming potential value of the farmland system is:

[0124] E 农田系统 =E 化肥施用 +E 粪便还田 +E 秸秆还田 +E 稻田种植 +E 能源消费1 +E 淋溶径流 +E 大气氮沉降 ;

[0125] Among them, E 农田系统 is the warming potential of the farmland system, E 化肥施用 The greenhouse gas potential value of nitrogen input into the soil for agricultural fertilizer application, E 粪便还田 E is the greenhouse gas potential value generated by the nitrogen input of livestock and poultry manure into the soil. 秸秆还田 E is the greenhouse gas potential value of nitrogen input into the soil when straw is returned to the field; 稻田种植 is the greenhouse gas potential value produced during rice cultivation, E 能源消费1 The greenhouse gas potential value of agricultural energy consumption; E 淋溶径流 is the gas potential value of greenhouse gases caused by leaching runoff, E 大气氮沉降 is the greenhouse gas potential of soils due to atmospheric nitrogen deposition;

[0126] E 化肥施用 =N 化肥 ×f cropland ×GWP1;

[0127] Among them, N化肥 is the nitrogen input of chemical fertilizer (nitrogen fertilizer, compound fertilizer), croplad is the emission factor of N2O greenhouse gas in agricultural soil, GWP1 is the warming potential coefficient of N2O greenhouse gas;

[0128] E 粪便还田 =N 粪便 ×f cropland ×GWP1;

[0129]

[0130] Among them, i is different types of livestock and poultry, N livestock is the nitrogen excretion of livestock and poultry feces and urine, is the rate of livestock and poultry manure returned to the field, A i is the annual output or end-year inventory of the i-th type of livestock and poultry, N i is the annual average nitrogen excretion of feces and urine of the i-th livestock and poultry, T i is the i-th livestock and poultry breeding cycle, the number of days per year is 365, it can be understood that, for example, the livestock and poultry types include dairy cows, non-dairy cows, buffaloes, sheep, pigs, poultry and are arranged in sequence, when i is 1, it is a dairy cow, N1 is the annual average excretion of feces and urine nitrogen by dairy cows, the following j, d, s, k, g and other category parameters can also be understood in this way, and will not be repeated later;

[0131] E 秸秆还田 =N 秸秆 ×f cropland ×GWP1;

[0132]

[0133] Where j is a different crop type (specifically, it may include rice, potatoes, sugarcane, peanuts, soybeans, dry grains, etc.), Str j is the straw yield of the jth crop, is the nitrogen content of the jth crop straw, is the rate of returning straw to the field, j is the crop yield, R str / see is the straw-grain ratio;

[0134] E 大气氮沉降 =∑N 化肥,粪便,秸秆 ×Rate dep ×f dep ×GWP1;

[0135] Among them, Rate dep is the observed volatilization coefficient of nitrogen input in the region, f dep Emission factors of N2O greenhouse gases from soils caused by atmospheric nitrogen deposition;

[0136] E 淋溶径流=∑N 化肥,粪便,秸秆 ×f leaching ×GWP1;

[0137] Among them, f leaching N2O greenhouse gas emission factors from soils due to nitrogen leaching and runoff losses;

[0138] E 稻田种植 =∑EF d ×AD d ×GWP2;

[0139] Among them, EF d is the CH4 emission factor for each type of paddy field, AD d is the rice planting area corresponding to the emission factor; d is the rice field type (specifically, it can include double-season early rice, double-season late rice, and single-season rice); GWP2 is the warming potential coefficient of CH4 greenhouse gas.

[0140] In some embodiments of the present invention, there is a livestock and poultry breeding system among the plurality of sub-emission systems, and the warming potential value of the livestock and poultry breeding system is:

[0141] E 畜禽养殖系统 =E 畜禽粪便 +E 牲畜肠道发酵 +E 能源消费2 ;

[0142] Among them, E 畜禽养殖系统 is the warming potential value of livestock and poultry breeding system, E 畜禽粪便 E is the warming potential of greenhouse gases emitted by livestock and poultry manure during livestock and poultry breeding. 牲畜肠道发酵 is the gas potential value of greenhouse gases produced by livestock intestinal fermentation during livestock and poultry breeding, E 能源消费2 The greenhouse gas potential value of energy consumption in livestock and poultry farming;

[0143] E 畜禽粪便 =∑EF i1 ×AD i ×GWP1+∑EF i2 ×AD i ×GWP2;

[0144] E 牲畜肠道发酵 =∑EF i3 ×AD i ×GWP2;

[0145] Among them, EF i1 EF is the emission factor of CH4 greenhouse gas from livestock and poultry manure; i2 EF is the emission factor of N2O greenhouse gas from livestock and poultry manure; i3 AD is the greenhouse gas emission factor of CH4 from enteric fermentation of livestock and poultry; iThe number of livestock and poultry slaughtered for different types;

[0146] Among the plurality of said sub-discharge systems there is an aquaculture system, and the warming potential of the aquaculture system is:

[0147] E 水产养殖系统 =E 水产养殖反硝化 +E 能源消费3 ;

[0148] Among them, E 水产养殖系统 is the warming potential of the aquaculture system, E 水产养殖反硝化 E is the greenhouse gas potential produced by denitrification in water bodies during nitrogen fertilizer application in aquaculture. 能源消费3 Greenhouse gas potentials for energy consumption in freshwater aquaculture;

[0149] E 水产养殖反硝化 =N Surplus ×f AQ ×GWP1;

[0150] Among them, N Sirplus represents the surplus nitrogen in the aquaculture system of the water plant, f AQ is the N2O emission factor for the residual nitrogen in the aquaculture system;

[0151] N Surplus =AQIN bait +AQIN fer +AQIN wf +AQIN dep +AQIN sprout -AQOUT prod ;

[0152] Among them, AQIN bait Amount of nitrogen input into feed for aquaculture systems, AQIN fer AQIN is the nitrogen fertilizer input for feed ratio in aquaculture system. dep Indicates atmospheric nitrogen deposition; AQIN sprout AQOUT is the nitrogen input to the aquaculture system. prod Nitrogen exports to aquaculture systems, capture fisheries and wild fisheries;

[0153] AQIN bait =∑W s ×N s ÷f s,bait ;

[0154] AQIN fer =AQIN bait ×f fer ;

[0155] AQINdep =Rate dep ×Area AQ ;

[0156] AQIN wf =∑WF s ×N s ;

[0157] AQIN sprout =Spruot s ×N s,sprout ;

[0158] AQOUT prod =∑W s ×N s ;

[0159] Among them, s is different types of aquatic products, W s is the output of the sth aquatic product, N s is the nitrogen content of the sth aquatic product, f s,bait is the feed conversion coefficient of the sth aquatic product, f fer Indicates the ratio parameter of nitrogen fertilizer in aquatic product feed, Rate dep is the observed volatilization coefficient of nitrogen input in the area, Area AQ is the freshwater aquaculture area in the study area, WF i is the output of wild aquatic products of category s, N i,F is the nitrogen content of the sth aquatic product, N i,sprout is the nitrogen content of the sth aquatic product seedlings;

[0160] Among the multiple emission sub-systems, there are urban green space and forest land systems, and the warming potential value of the urban green space and forest land system is:

[0161] E 城市绿地与林地系统 =E 森林脱氮 +E 城市草地化肥施用 +E 能源消费4 ;

[0162] Among them, E 城市绿地与林地系统 is the warming potential of the urban green space and forest system, E 森林脱氮 Gas potential of greenhouse gases produced by forest denitrification, E 城市草地化肥施用 Greenhouse gas potential from soil fertilization for urban grasslands, E 能源消费4 Gas potentials for greenhouse gases produced for energy consumption in the forest sector;

[0163] E 城市草地化肥施用 =UG fer ×UG area ×f UG ×GWP1;

[0164] E 森林脱氮 =∑FR area ×f forest ×GWP1;

[0165] Among them, UG fer is the amount of nitrogen fertilizer applied per unit area of ​​urban green lawn, UG area is the area of ​​urban green lawn, f UG is the emission factor of soil N2O greenhouse gas, FR area is the area corresponding to different forest types, f forest is the forest N2O greenhouse gas emission factor corresponding to different forest types;

[0166] Among the plurality of said sub-discharge systems there is a sewage treatment system, and the warming potential value of the sewage treatment system is:

[0167] E 污水处理系统 =E 工业废水处理过程 +E 生活废水处理过程 ;

[0168] Among them, E 污水处理系统 is the warming potential of the sewage treatment system, E 工业废水处理过程 is the greenhouse gas potential value generated during the treatment of industrial wastewater, E 生活废水处理过程 The greenhouse gas potential value generated during the treatment of domestic wastewater;

[0169] E 生活废水处理过程 =TN1×f TN,removed,1 ×EF f1 ×GWP1+COD1×f COD,removed,1 ×EF f2 ×GWP2;

[0170] E 工业废水处理过程 =TN2×f TN,removed,2 ×EF f1 ×GWP1+COD2×f COD,removed,2 ×EF f2 ×GWP2;

[0171]

[0172] Among them, P is the permanent urban population (10,000 people), X TN is the per capita daily total nitrogen pollution intensity (g / person / day), X COD is the per capita daily chemical oxygen demand pollution intensity (g / person / day), f collectis the collection rate of urban domestic sewage in each region, TN1 is the total nitrogen content of domestic wastewater collected through the urban domestic sewage pipe network and entering the domestic sewage treatment plant, COD1 is the chemical oxygen demand collected through the urban domestic sewage pipe network and entering the domestic sewage treatment plant, IP is the output of industrial products, Y TN Y is the pollution coefficient of producing unit industrial product TN, COD is the pollution coefficient of COD per unit of industrial product production, TN2 is the total nitrogen produced by industrial wastewater product production, COD2 is the chemical oxygen demand produced by industrial wastewater product production, and f TN,removed,1 is the average TN removal rate of sewage treatment plants in each region of the basin, f COD,removed,1 is the average COD treatment rate of sewage treatment plants in various regions within the basin; TN,removed,2 is the average treatment rate of TN in the terminal treatment of industrial wastewater, f COD,removed,2 The average treatment rate of COD in the terminal treatment of industrial wastewater, EF f1 is the emission factor of N2O greenhouse gas in the sewage treatment process, EF f2 is the emission factor of CH4 greenhouse gas in the sewage treatment process;

[0173] Among the plurality of said sub-emission systems, there is a solid waste disposal system, and the warming potential value of the solid waste disposal system is:

[0174] E 固废处置系统 =E 垃圾填埋 +E 垃圾焚烧 +E 能源消费5 ;

[0175] Among them, E 固废处置系统 is the greenhouse gas potential of the solid waste disposal system, E 垃圾填埋 is the greenhouse gas potential value generated during the landfill of domestic waste, E 垃圾焚烧 is the potential value of greenhouse gases produced during the incineration of domestic waste, E 能源消费5 The gas potential of greenhouse gases generated by the energy input of waste incineration;

[0176] E 垃圾填埋 =∑ g MSW g，填埋 ×DOC g ×MCF×DOC f ×F×a1×(1-R1)×(1-OX)×GWP2;

[0177] Among them, g is the type of different waste components, MSW g，填埋 is the landfill volume of different types of waste components, DOC i is the degradable organic carbon fraction of different garbage components g, MCF is the correction factor for CH4 greenhouse gas, DOCf is the ratio of decomposed degradable organic carbon, F is the content of CH4 greenhouse gas in landfill gas, R1 is the recovery rate of CH4 greenhouse gas, OX is the oxidation factor, a1 is the relative molecular mass ratio of CH4 and C, specifically 16 / 12;

[0178] E 垃圾焚烧 =∑ g MSW g，焚烧 ×dm g ×FC g ×FCF I ×OF×a2+MSW g，焚烧 ×EF b ×GWP1;

[0179] Among them, MSW g，焚烧 is the incineration amount of different types of garbage components, dm g is the dry matter content of different garbage components g, FC g is the fossil carbon component of different garbage components g, OF is the oxidation factor, EF is b is the emission factor of N2O greenhouse gas in the waste incineration process, a2 is the relative molecular mass ratio of CO2 and C, specifically 44 / 12;

[0180] Among the plurality of said sub-discharge systems, there is a surface water system, and the warming potential value of the surface water system is:

[0181] E 地表水系统 =E 工业废水-地表水排放过程 +E 生活废水-地表水排放过程 ; E 地表水系统 is the greenhouse gas potential of the surface water system,

[0182] E 工业废水-地表水排放过程 It is the greenhouse gas potential value of pollutants produced by industrial wastewater after direct discharge or discharge into surface water after treatment.

[0183] E 生活废水-地表水排放过程 It is the greenhouse gas potential value of pollutants produced by direct discharge of domestic wastewater or after treatment into surface water;

[0184] E 生活废水-地表水排放过程 =(TN3+TN4)×EF f3 ×GWP1+(COD3+COD4)×EF f4 ×GWP2;

[0185] E 工业废水-地表水排放过程 =TN5×EF f3 ×GWP1+COD5×EF f4 ×GWP2;

[0186]

[0187] Among them, TN3 is the total nitrogen content of domestic wastewater entering the surface water after being treated in the sewage treatment plant, COD3 is the chemical oxygen demand of domestic wastewater entering the surface water after being treated in the sewage treatment plant, TN4 is the total nitrogen content of domestic wastewater directly discharged into the surface water system without being collected by the sewage pipe network, COD4 is the chemical oxygen demand of domestic wastewater directly discharged into the surface water system without being collected by the sewage pipe network, TN5 is the total nitrogen content of industrial wastewater entering the surface water system after treatment, COD5 is the chemical oxygen demand of industrial wastewater entering the surface water system after treatment, EF f3 is the emission factor of N2O greenhouse gas from wastewater discharged into rivers and lakes, EF f4 It is the emission factor of CH4 greenhouse gas from off-site discharge of wastewater into rivers and lakes.

[0188] In some embodiments of the present invention, the E in the above sub-emission system 能源消费1 、E 能源消费2 、E 能源消费3 、E 能源消费4 and E 能源消费5 The following unit models can be used for calculation:

[0189]

[0190] in, is the gas potential value of greenhouse gas caused by the combustion of any fossil fuel of energy consumption 1, energy consumption 2, energy consumption 3, energy consumption 4, and energy consumption 5, k is the type of energy consumed (specifically, it can be coal, natural gas, oil, etc.), is the k-th energy consumption, NCV k is the net calorific value of the k-th energy source, is the emission factor corresponding to each of the greenhouse gases produced by the k-th energy consumption, The carbon oxidation rate factor corresponding to the various greenhouse gases produced by the k-th energy consumption, GWP (1，2，3) is the corresponding warming potential coefficient of N2O, CH4 and CO2 greenhouse gases.

[0191] The CO2, CH 4, The warming potential values ​​of N2O are calculated separately, and then the warming potential values ​​of multiple greenhouse gases are added together to obtain the warming potential value of the energy combustion. Finally, the warming potential value of each energy combustion is calculated.

[0192] In some embodiments of the present invention, there is a human consumption system among the plurality of emission sub-systems, and the warming potential value of human consumption can be calculated using the following unit model:

[0193] E 人类消费系统 =E 居民家庭能源消费 +E 交通运输能源消费 ;

[0194] Among them, E 人类消费系统 is the greenhouse gas potential of human consumption system, E 居民家庭能源消费 The gas potential value of greenhouse gases generated by household energy consumption, E 交通运输能源消费 Gas potentials for greenhouse gases produced by transport energy consumption;

[0195] Among the plurality of emission sub-systems, there are industrial and energy systems, and the warming potential of the industrial and energy systems is:

[0196] E 工业与能源系统 =∑E 工业分行业能源消费 ;

[0197] E 工业与能源系统 is the total greenhouse gas emissions from industry and energy systems,

[0198] E 工业分行业能源消费 Greenhouse gas emissions from energy consumption by industrial sub-sector;

[0199] E 居民家庭能源消费 、E 交通运输能源消费 and E 工业分行业能源消费 Both are:

[0200]

[0201] in, is the gas potential value of greenhouse gases caused by the combustion of any energy fossil fuel in household energy consumption, transportation energy consumption, and industrial sub-sector energy consumption, k is the type of energy consumed, is the k-th energy consumption, NCV k is the net calorific value of the k-th energy source, is the emission factor corresponding to each of the greenhouse gases produced by the k-th energy consumption, The carbon oxidation rate factor corresponding to the various greenhouse gases produced by the k-th energy consumption, GWP (1，2，3) is the corresponding warming potential coefficient of N2O, CH4 and CO2 greenhouse gases.

[0202] Specifically, E 居民家庭能源消费 、E 交通运输能源消费 and E 工业分行业能源消费 The calculation of can adopt the calculation method of energy consumption 1, energy consumption 2, energy consumption 3, energy consumption 4, and energy consumption 5 mentioned above, which will not be described in detail here.

[0203] The statistical system according to the second aspect of the present invention is used to execute the industrial nitrogen metabolism greenhouse gas statistical method disclosed in any of the above embodiments.

[0204] According to the control device of the third aspect of the present invention, the control device includes a memory and a processor, the memory stores a computer program, and the processor implements an industrial nitrogen metabolism greenhouse gas statistics method disclosed in any of the above embodiments when executing the computer program.

[0205] The control device can be any intelligent terminal including a central computer, a remote device terminal computer, etc.

[0206] According to the computer-readable storage medium of the fourth aspect of the present invention, the computer-readable storage medium stores a computer program, characterized in that when the computer program is executed by a processor, an industrial nitrogen metabolism greenhouse gas statistics method disclosed in any of the above embodiments is implemented.

[0207] The memory, as a non-transient computer-readable storage medium, can be used to store non-transient software programs and non-transient computer executable programs. In addition, the memory may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some embodiments, the memory may optionally include a memory remotely disposed relative to the processor, and these remote memories may be connected to the processor via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0208] The embodiments described in the embodiments of the present application are intended to more clearly illustrate the technical solutions of the embodiments of the present application and do not constitute a limitation on the technical solutions provided in the embodiments of the present application. Those skilled in the art will appreciate that with the evolution of technology and the emergence of new application scenarios, the technical solutions provided in the embodiments of the present application are also applicable to similar technical problems.

[0209] Those skilled in the art will appreciate that the technical solutions shown in the figures do not constitute a limitation on the embodiments of the present application, and may include more or fewer steps than shown in the figures, or a combination of certain steps, or different steps.

[0210] The device embodiments described above are merely illustrative, and the units described as separate components may or may not be physically separated, that is, they may be located in one place or distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0211] Those skilled in the art will appreciate that all or some of the steps in the methods disclosed above, and the functional modules / units in the systems and devices may be implemented as software, firmware, hardware, or a suitable combination thereof.

[0212] The terms "first", "second", "third", "fourth", etc. (if any) in the specification of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0213] The preferred embodiments of the present invention are described above with reference to the accompanying drawings, but the scope of the rights of the present invention is not limited thereto. Any modification, equivalent substitution and improvement made by a person skilled in the art without departing from the scope and essence of the present invention should be within the scope of the rights of the present invention.

[0214] The technical features of the above-described embodiments may be arbitrarily combined. To make the description concise, not all possible combinations of the technical features in the above-described embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0215] Although the embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the claims and their equivalents.

Claims

1. The statistical method of greenhouse gases from industrial nitrogen metabolism is used to calculate the total warming potential of greenhouse gases emitted in industrial clusters, including: There are many types of greenhouse gases, and the characteristic is that the statistical method of industrial nitrogen metabolism greenhouse gases includes: Divide the industrial cluster into multiple emission sub-systems; According to the differences of each sub-emission system, a sub-system emission model corresponding to each sub-emission system is established, wherein the establishment of the sub-system emission model includes establishing a material flow framework based on the material flow analysis method according to the flow process of the input material in the sub-emission system, and obtaining a relationship model between the input material in the sub-emission system and the emission of various greenhouse gases in the material flow framework based on the emission factor method; Each sub-emission system calculates the warming potential value based on its own sub-system emission model, which includes: Obtain input material data of the corresponding sub-emission system, and calculate the gas emissions of various greenhouse gases in the sub-system emission model according to the input material data; The gas emissions of various greenhouse gases are normalized according to the warming potential coefficients of the corresponding greenhouse gases and summed to obtain the warming potential value of the sub-emission system; The warming potential values ​​of the various emission subsystems are integrated to obtain a total warming potential value.

2. The method for calculating greenhouse gas statistics of industrial nitrogen metabolism according to claim 1, characterized in that: The subsystem emission model corresponding to each sub-emission system includes: Based on the material flow analysis method, the subsystem emission model is divided into a plurality of unit models according to different types of activities producing biochemical reactions in the sub-emission system, wherein the unit model represents the relationship between the input material and at least one greenhouse gas emission; Calculating the gas emissions of various greenhouse gases in the subsystem emission model according to the input material data includes: Input material data are fed into each unit model calculation to obtain unit emissions of various greenhouse gases; The unit emissions of the same type of greenhouse gas are added together to obtain the gas emissions of that greenhouse gas.

3. The method for calculating greenhouse gas statistics of industrial nitrogen metabolism according to claim 2, characterized in that: The gas emissions of various greenhouse gases are normalized according to the warming potential coefficients of the corresponding greenhouse gases and summed to obtain the warming potential value of the sub-emission system, including: The product of the gas emission of the greenhouse gas and the warming potential coefficient of the greenhouse gas is taken as the gas potential value of the greenhouse gas; The gas potential values ​​of multiple greenhouse gases are summed to obtain the warming potential value of the sub-emission system.

4. The industrial nitrogen metabolism greenhouse gas statistical method according to claim 2, characterized in that: Among the plurality of sub-emission systems, there is a farmland system, and the warming potential value of the farmland system is: AND 农田系统 =And 化肥施用 +E 粪便还田 +E 秸秆还田 +E 稻田种植 + AND 能源消费1 +E 淋溶径流 +E 大气氮沉降 ; Among them, E 农田系统 is the warming potential of the farmland system, E 化肥施用 The greenhouse gas potential value of nitrogen input into the soil for agricultural fertilizer application, E 粪便还田 E is the greenhouse gas potential value generated by the nitrogen input of livestock and poultry manure into the soil. 秸秆还田 E is the greenhouse gas potential value of nitrogen input into the soil when straw is returned to the field; 稻田种植 is the greenhouse gas potential value produced during rice cultivation, E 能源消费1 The greenhouse gas potential value of agricultural energy consumption; E 淋溶径流 is the gas potential value of greenhouse gases caused by leaching runoff, E 大气氮沉降 is the greenhouse gas potential of soils due to atmospheric nitrogen deposition; E 化肥施用 =N 化肥 ×f cropland ×GWP1; Among them, N 化肥 is the fertilizer nitrogen input, f croplad is the emission factor of N2O greenhouse gas in agricultural soil, GWP1 is the warming potential coefficient of N2O greenhouse gas; E 粪便还田 =N 粪便 ×f cropland ×GWP1; Among them, i is different types of livestock and poultry, N licestock is the nitrogen excretion of livestock and poultry feces and urine, is the rate of livestock and poultry manure returned to the field, A i is the annual output or end-year inventory of the i-th type of livestock and poultry, N i is the annual average nitrogen excretion of feces and urine of the i-th livestock and poultry, T i is the i-th livestock and poultry breeding cycle, with 365 days per year; E 秸秆还田 =N 秸秆 ×f cropland ×GWP1; Among them, j is different crop types, Str j is the straw yield of the jth crop, is the nitrogen content of the jth crop straw, is the rate of returning straw to the field, j is the crop yield, R str / see is the straw-grain ratio; AND 大气氮沉降 =∑N 化肥，粪便，秸秆 ×Rate dep ×f dep ×GWP1; Among them, Rate dep is the observed volatilization coefficient of nitrogen input in the region, f dep Emission factors of N2O greenhouse gases from soils caused by atmospheric nitrogen deposition; AND 淋溶径流 =∑N 化肥，粪便，秸秆 ×f leaching ×GWP1; Among them, f leaching N2O greenhouse gas emission factors from soils due to nitrogen leaching and runoff losses; AND 稻田种植 =∑EF d ×AD d ×GWP2; Among them, EF d is the CH4 emission factor for each type of paddy field, AD d is the rice planting area corresponding to the emission factor; d is the rice field type, and GWP2 is the greenhouse warming potential coefficient of CH4.

5. The industrial nitrogen metabolism greenhouse gas statistical method according to claim 4, characterized in that: Among the plurality of sub-emission systems, there is a livestock and poultry breeding system, and the warming potential value of the livestock and poultry breeding system is: AND 畜禽养殖系统 =And 畜禽粪便 +E 牲畜肠道发酵 +E 能源消费2 ; Among them, E 畜禽养殖系统 is the warming potential value of livestock and poultry breeding system, E 畜禽粪便 E is the warming potential of greenhouse gases emitted by livestock and poultry manure during livestock and poultry breeding. 牲畜肠道发酵 is the gas potential value of greenhouse gases produced by livestock intestinal fermentation during livestock and poultry breeding, E 能源消费2 The greenhouse gas potential value of energy consumption in livestock and poultry farming; AND 畜禽粪便 =∑EF i1 ×AD i ×GWP2+∑EF i2 ×AD i ×GWP2; AND 牲畜肠道发酵 =∑EF i3 ×AD i ×GWP2; Among them, EF i1 EF is the emission factor of CH4 greenhouse gas from livestock and poultry manure; i2 EF is the emission factor of N2O greenhouse gas from livestock and poultry manure; i3 AD is the greenhouse gas emission factor of CH4 from enteric fermentation of livestock and poultry; i The number of livestock and poultry slaughtered for different types; Among the plurality of said sub-discharge systems there is an aquaculture system, and the warming potential of the aquaculture system is: AND 水产养殖系统 =And 水产养殖反硝化 +E 能源消费3 ; Among them, E 水产养殖系统 is the warming potential of the aquaculture system, E 水产养殖反硝化 E is the greenhouse gas potential produced by denitrification in water bodies during nitrogen fertilizer application in aquaculture. 能源消费3 Greenhouse gas potentials for energy consumption in freshwater aquaculture; E 水产养殖反硝化 =N Surplus ×f AQ ×GWP1; Among them, N Surplus represents the surplus nitrogen in the aquaculture system of the water plant, f AQ is the N2O emission factor for the residual nitrogen in the aquaculture system; N Surplus =AQIN bait +AQIN fer +AQIN wf +AQIN dep +AQIN sprout -AQOUT prod ; Among them, AQIN bait Amount of nitrogen input into feed for aquaculture systems, AQIN fer AQIN is the nitrogen fertilizer input for feed ratio in aquaculture system. dep Indicates atmospheric nitrogen deposition; AQIN sprout AQOUT is the nitrogen input to the aquaculture system. prod Nitrogen exports to aquaculture systems, capture fisheries and wild fisheries; AQIN bait =∑W s ×N s ÷f s，bait ; AQIN fer =AQIN bait ×f fer ; AQIN dep =Rate dep ×Area AQ ; AQIN wf =∑WF s ×N s ; AQIN sprout =Sprues s ×N s,sprout ; AQOUT prod =∑W s ×N s ; Among them, s is different types of aquatic products, W s is the output of the sth aquatic product, N s is the nitrogen content of the sth aquatic product, f s，bait is the feed conversion coefficient of the sth aquatic product, f fer Indicates the ratio parameter of nitrogen fertilizer in aquatic product feed, Rate dep is the observed volatilization coefficient of nitrogen input in the area, Area AQ is the freshwater aquaculture area in the study area, WF i is the output of wild aquatic products of category s, N i，F is the nitrogen content of the sth aquatic product, N i，sprout is the nitrogen content of the sth aquatic product seedlings; Among the multiple emission sub-systems, there are urban green space and forest land systems, and the warming potential value of the urban green space and forest land system is: AND 城市绿地与林地系统 =And 森林脱氮 +E 城市草地化肥施用 +E 能源消费4 ; Among them, E 城市绿地与林地系统 is the warming potential of the urban green space and forest system, E 森林脱氮 Gas potential of greenhouse gases produced by forest denitrification, E 城市草地化肥施用 Greenhouse gas potential from soil fertilization for urban grasslands, E 能源消费4 Gas potentials for greenhouse gases produced for energy consumption in the forest sector; Yes 城市草地化肥施用 =AND fer ×UG area ×f UG ×GWP1; AND 森林脱氮 =∑FR area ×f forest ×GWP1; Among them, UG fer is the amount of nitrogen fertilizer applied per unit area of ​​urban green lawn, UG area is the area of ​​urban green lawn, f UG is the emission factor of soil N2O greenhouse gas, FR area is the area corresponding to different forest types, f forest is the forest N2O greenhouse gas emission factor corresponding to different forest types; Among the plurality of said sub-discharge systems there is a sewage treatment system, and the warming potential value of the sewage treatment system is: AND 污水处理系统 =And 工业废水处理过程 +E 生活废水处理过程 ; Among them, E 污水处理系统 is the warming potential of the sewage treatment system. E 工业废水处理过程 is the greenhouse gas potential value generated during the treatment of industrial wastewater, E 生活废水处理过程 The greenhouse gas potential value generated during the treatment of domestic wastewater; AND 生活废水处理过程 =TN1×f TN，removed,1 ×EF f1 ×GWP1+COD1×f COD,removed，1 ×EF f2 ×GWP2; AND 工业废水处理过程 =TN2×f TN，removed,2 ×EF f1 ×GWP1+COD2×f COD,removed，2 ×EF f2 ×GWP2; Among them, P is the permanent population of the city, X TN is the daily per capita total nitrogen pollution intensity, X COD is the pollution intensity of chemical oxygen demand per capita per day, fc ollect is the collection rate of urban domestic sewage in each region, TN1 is the total nitrogen content of domestic wastewater collected through the urban domestic sewage pipe network and entering the domestic sewage treatment plant, COD1 is the chemical oxygen demand collected through the urban domestic sewage pipe network and entering the domestic sewage treatment plant, IP is the output of industrial products, Y TN Y is the pollution coefficient of producing unit industrial product TN, COD is the pollution coefficient of COD per unit of industrial product production, TN2 is the total nitrogen produced by industrial wastewater product production, COD2 is the chemical oxygen demand produced by industrial wastewater product production, and f TN，removed，1 is the average TN removal rate of sewage treatment plants in each region of the basin, f COD,removed，1 is the average COD treatment rate of sewage treatment plants in various regions within the basin; TN，removed，2 is the average treatment rate of TN in the terminal treatment of industrial wastewater, f COD,removed,2 The average treatment rate of COD in the terminal treatment of industrial wastewater, EF f1 is the emission factor of N2O greenhouse gas in the sewage treatment process, EF f2 is the emission factor of CH4 greenhouse gas in the sewage treatment process; Among the plurality of said sub-emission systems, there is a solid waste disposal system, and the warming potential value of the solid waste disposal system is: AND 固废处置系统 =And 垃圾填埋 +E 垃圾焚烧 +E 能源消费5 ; Among them, E 固废处置系统 is the greenhouse gas potential of the solid waste disposal system, E 垃圾填埋 is the greenhouse gas potential value generated during the landfill of domestic waste, E 垃圾焚烧 is the potential value of greenhouse gases produced during the incineration of domestic waste, E 能源消费5 The gas potential of greenhouse gases generated by the energy input of waste incineration; E 垃圾填埋 =∑ g MSW g，填埋 ×DOC g ×MCF×DOC f ×F×a1×(1-R1)×(1-OX)×GWP2; Among them, g is the type of different waste components, MSW g，填埋 is the landfill volume of different types of waste components, DOC i is the degradable organic carbon fraction of different garbage components g, MCF is the correction factor for CH4 greenhouse gas, DOC f is the proportion of decomposed degradable organic carbon, F is the content of CH4 greenhouse gas in landfill gas, R1 is the recovery rate of CH4 greenhouse gas, OX is the oxidation factor, and a1 is the relative molecular mass ratio of CH4 and C; E 垃圾焚烧 ＝ <h2 style=";text-align:left;direction:ltr">∑<h2 style=";text-align:left;direction:ltr"> g <h2 style=";text-align:left;direction:ltr"> MSW<h2 style=";text-align:left;direction:ltr"> g，焚烧 <h2 style=";text-align:left;direction:ltr"> ×dm<h2 style=";text-align:left;direction:ltr"> g <h2 style=";text-align:left;direction:ltr"> ×FC<h2 style=";text-align:left;direction:ltr"> g <h2 style=";text-align:left;direction:ltr"> ×FCF<h2 style=";text-align:left;direction:ltr"> I <h2 style=";text-align:left;direction:ltr"> ×OF×a2+MSW<h2 style=";text-align:left;direction:ltr"> g，焚烧 <h2 style=";text-align:left;direction:ltr"> ×EF<h2 style=";text-align:left;direction:ltr"> b <h2 style=";text-align:left;direction:ltr"> ×GWP1; Among them, MSW g，焚烧 is the incineration amount of different types of garbage components, dm g is the dry matter content of different garbage components g, FC g is the fossil carbon component of different garbage components g, OF is the oxidation factor, EF is b is the emission factor of N2O greenhouse gas in the waste incineration process, a2 is the relative molecular mass ratio of CO2 and C; Among the plurality of said sub-discharge systems, there is a surface water system, and the warming potential value of the surface water system is: AND 地表水系统 =And 工业废水-地表水排放过程 +E 生活废水-地表水排放过程 ; E 地表水系统 is the greenhouse gas potential of the surface water system, E 工业废水-地表水排放过程 It is the greenhouse gas potential value of pollutants produced by industrial wastewater after direct discharge or discharge into surface water after treatment. E 生活废水-地表水排放过程 It is the greenhouse gas potential value of pollutants produced by direct discharge of domestic wastewater or after treatment into surface water; E 生活废水-地表水排放过程 =(TN3+TN4)×EF f3 ×GWP1+(COD3+COD4)×EF f4 ×GWP2; AND 工业废水-地表水排放过程 =TN5×EF f3 ×GWP1+COD5×EF f4 ×GWP2; Among them, TN3 is the total nitrogen content of domestic wastewater entering the surface water after being treated in the sewage treatment plant, COD3 is the chemical oxygen demand of domestic wastewater entering the surface water after being treated in the sewage treatment plant, TN4 is the total nitrogen content of domestic wastewater directly discharged into the surface water system without being collected by the sewage pipe network, COD4 is the chemical oxygen demand of domestic wastewater directly discharged into the surface water system without being collected by the sewage pipe network, TN5 is the total nitrogen content of industrial wastewater entering the surface water system after treatment, COD5 is the chemical oxygen demand of industrial wastewater entering the surface water system after treatment, EF f3 is the emission factor of N2O greenhouse gas from wastewater discharged into rivers and lakes, EF f4 It is the emission factor of CH4 greenhouse gas from off-site discharge of wastewater into rivers and lakes.

6. The industrial nitrogen metabolism greenhouse gas statistics method according to claim 5, characterized in that: E 能源消费1 、E 能源消费2 、E 能源消费3 、E 能源消费4 and E 能源消费5 Both are: in, is the gas potential value of greenhouse gas caused by the combustion of fossil fuels in any one of energy consumption 1, energy consumption 2, energy consumption 3, energy consumption 4, and energy consumption 5, k is the type of energy consumed, is the k-th energy consumption, NCVk is the net calorific value of the k-th energy, is the emission factor corresponding to each of the greenhouse gases produced by the k-th energy consumption, The carbon oxidation rate factor corresponding to the various greenhouse gases produced by the k-th energy consumption, GWP (1，2，3) is the corresponding warming potential coefficient of N2O, CH4 and CO2 greenhouse gases.

7. The industrial nitrogen metabolism greenhouse gas statistics method according to claim 5, characterized in that: Among the plurality of emission sub-systems, there is a human consumption system, and the warming potential value of human consumption is: AND 人类消费系统 =And 居民家庭能源消费 +E 交通运输能源消费 ; Among them, E 人类消费系统 is the greenhouse gas potential of human consumption system, E 居民家庭能源消费 The gas potential value of greenhouse gases generated by household energy consumption, E 交通运输能源消费 Gas potentials for greenhouse gases produced by transport energy consumption; Among the plurality of emission sub-systems, there are industrial and energy systems, and the warming potential of the industrial and energy systems is: AND 工业与能源系统 =∑E 工业分行业能源消费 ; E 工业与能源系统 is the total greenhouse gas emissions from industry and energy systems, E 工业分行业能源消费 Greenhouse gas emissions from energy consumption by industrial sub-sector; E 居民家庭能源消费 、E 交通运输能源消费 and E 工业分行业能源消费 Both are: in, is the gas potential value of greenhouse gases caused by the combustion of any energy fossil fuel in household energy consumption, transportation energy consumption, and industrial sub-sector energy consumption, k is the type of energy consumed, is the k-th energy consumption, NCVk is the net calorific value of the k-th energy, is the emission factor corresponding to each of the greenhouse gases produced by the k-th energy consumption, The carbon oxidation rate factor corresponding to the various greenhouse gases produced by the k-th energy consumption, GWP (1，2，3) is the corresponding warming potential coefficient of N2O, CH4 and CO2 greenhouse gases.

8. A statistical system, characterized in that: Used to implement the industrial nitrogen metabolism greenhouse gas statistics method as described in any one of claims 1 to 7.

9. A control device, characterized in that: The control device includes a memory and a processor, the memory stores a computer program, and the processor implements the industrial nitrogen metabolism greenhouse gas statistics method according to any one of claims 1 to 7 when executing the computer program.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the industrial nitrogen metabolism greenhouse gas statistics method according to any one of claims 1 to 7 is implemented.

Images