CN115204729A - Water and soil resource carbon emission accounting method, device, equipment and storage medium - Google Patents

Water and soil resource carbon emission accounting method, device, equipment and storage medium Download PDF

Info

Publication number
CN115204729A
CN115204729A CN202210909827.9A CN202210909827A CN115204729A CN 115204729 A CN115204729 A CN 115204729A CN 202210909827 A CN202210909827 A CN 202210909827A CN 115204729 A CN115204729 A CN 115204729A
Authority
CN
China
Prior art keywords
water
land
carbon emission
utilization
soil
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210909827.9A
Other languages
Chinese (zh)
Inventor
江文渊
张征云
陈启华
赵翌晨
闫佩
张彦敏
罗彦鹤
罗航
郭健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin Huanke Environmental Planning Technology Development Co ltd
Tianjin Academy of Ecological and Environmental Sciences
Original Assignee
Tianjin Huanke Environmental Planning Technology Development Co ltd
Tianjin Academy of Ecological and Environmental Sciences
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin Huanke Environmental Planning Technology Development Co ltd, Tianjin Academy of Ecological and Environmental Sciences filed Critical Tianjin Huanke Environmental Planning Technology Development Co ltd
Priority to CN202210909827.9A priority Critical patent/CN115204729A/en
Publication of CN115204729A publication Critical patent/CN115204729A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067—Enterprise or organisation modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
    • G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/80—Management or planning
    • Y02P90/84—Greenhouse gas [GHG] management systems

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Multimedia (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method, a device, equipment and a storage medium for accounting water and soil resource carbon emission. The accounting method comprises the following steps: the total carbon emission amount in the water resource development and utilization is obtained by accounting the sum of the carbon emission amounts of different links in the water resource development and utilization; obtaining the total carbon emission amount in the land resource development and utilization by calculating the difference between the total carbon emission amount of different carbon source lands in the land resource development and utilization and the total carbon emission amount of different carbon sink lands; and obtaining the total carbon emission in the development and utilization of the water and soil resources according to the sum of the total carbon emission in the development and utilization of the water resources and the total carbon emission in the development and utilization of the land resources. The invention can obtain accurate accounting of carbon emission in water and soil resource development and utilization under the condition of considering the coupling effect of the water and soil resources based on multi-source data.

Description

Water and soil resource carbon emission accounting method, device, equipment and storage medium
Technical Field
The invention relates to the technical field of carbon emission accounting methods, in particular to the technical field of carbon emission accounting methods generated in development and utilization of water and soil resources.
Background
Human life and production activities require investment in water and land resources. A large amount of energy is consumed in the development and utilization process of water resources and land resources, carbon emission is generated, the method is an indispensable important content in carbon emission accounting, and the method has important significance for realizing low-carbon development and intensive utilization of water and soil resources.
In the prior art, a carbon emission accounting method for water and soil resource development and utilization mainly focuses on single factors of water and soil resources, two basic resources are split, and the influence of the coupling effect of the two resources on resource energy allocation and the carbon emission effect are not considered. The accounting method for the water resource related carbon emission mainly comprises the steps of calculating the carbon water footprint of the industry and evaluating the carbon emission of a water system of a city, wherein the accounting is carried out on the basis of statistical data of industrial input and output, and used data are difficult to obtain on time by year and poor in timeliness; the latter relies on real-time detailed operation data of each water system, and has high data acquisition difficulty, poor operability in practical application and incapability of being applied to large-space-scale accounting; the accounting method for the land resource related carbon emission is mainly based on land utilization classification data, energy consumption data are coupled, and accounting is performed on a large space scale.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a method for accurately accounting carbon emission in development and utilization of water and soil resources based on multi-source data such as monitoring data, remote sensing images, survey data and the like, comprehensively considering the coupling effect of the water and soil resources and having excellent timeliness.
The invention also aims to provide a storage medium, a device and equipment which can carry out automation and intelligent implementation on the carbon emission accounting method.
The technical scheme of the invention is as follows:
a method for accounting carbon emission in development and utilization of water and soil resources comprises the following steps:
s1, accounting is carried out on carbon emission of different links in water resource development and utilization, and the total carbon emission in the water resource development and utilization is obtained through the sum of the carbon emission of the different links;
s2, the carbon emission amount of different carbon source lands and the carbon sink amount of different carbon sink lands in the land resource development and utilization are calculated, and the total carbon emission amount in the land resource development and utilization is obtained through the difference between the sum of the carbon emission amount of the different carbon source lands and the sum of the carbon sink amount of the different carbon sink lands;
s3, carrying out total accounting on carbon emission in water and soil resource development and utilization, wherein the total accounting model is as follows:
C net ＝C W +C L (22)
wherein, C net Total carbon emission for water and soil resource development and utilization, C w Represents the total carbon emission in the development and utilization of water resources, C L Representing the total carbon emission in the development and utilization of the land resourceAn amount;
wherein, the different links of water resource development and utilization comprise water taking, water making, water delivery, water treatment and water using links, and the water using link comprises a first to a third industrial water using links and a resident domestic water using link; the carbon source land comprises cultivated land, other agricultural land except the cultivated land, residential site, industrial and mining land and transportation land; the carbon sink lands include cultivated lands, woodlands, gardens, grasslands, water areas and water conservancy facilities, and unused lands.
According to some preferred embodiments of the invention, S1 comprises:
s11, carrying out carbon emission accounting on water taking, water making, water conveying and water treatment links in water resource development and utilization through the following calculation model:
Figure BDA0003773638510000021
wherein, W h CO generated by energy consumption in the water taking, water making, water delivery and water treatment links 2 Total amount of emissions; α =44/12, for conversion of carbon to CO 2 The conversion coefficient of (2); e i Terminal consumption of the ith fossil energy, NCV i Is the low calorific value, delta, of the ith fossil energy i Is the carbon content of the ith fossil energy, OR i Combustion oxidation rate for the ith fossil energy source; EH j For electric or thermal end consumption, beta j As electric or thermal CO 2 Discharge coefficient, where j =1 represents thermal power, and j =2 represents electrical power;
s12, carrying out carbon emission accounting on the water consumption ring section in water resource development and utilization through the following calculation model:
the first industrial water link carbon emission calculation model comprises:
W a ＝αG i h(2)
wherein, W a CO for irrigation water 2 Discharge amount, G i Is the irrigation area; h is carbon emission coefficient of unit irrigation area, and h =266.48kgC/hm can be selected 2 ；
A water link carbon emission calculation model for the second industry and the third industry:
W e ＝W g +W c ＝E g β heat generation +E c β Heat generation (9)
W s ＝E s β Heat generation (10)
Wherein, W e 、W g 、W c 、W s CO of water consumption links of second industry, building industry and third industry respectively 2 Discharge amount of, wherein E g 、E c 、E s Thermal terminal consumption, beta, for industry, construction, and third industries, respectively Heat generation Is thermal CO 2 The discharge coefficient;
a model for calculating carbon emission in a resident domestic water link:
Figure BDA0003773638510000031
wherein, W b CO as domestic water for residents 2 Emission amount, beta Heat generation Is thermal CO 2 Discharge coefficient, p is the density of water, Q Drinking water The water consumption for human body, P is the population number of the nuclear region, T is the average annual temperature, E r Terminal consumption of heat for residents' life, Q Washing machine Water consumption for bathing T i ' is the average air temperature in each season, i =1,2,3,4, which respectively represents spring, summer, autumn and winter seasons;
s13, developing and utilizing the total carbon emission C in the water resource by the following calculation model w And (4) performing accounting:
C W ＝W h +W a +W e +W s +W b (11-2)。
according to some preferred embodiments of the present invention, the terminal consumption of fossil energy and the terminal consumption of electricity and heat in the water intake, production, delivery and treatment links are obtained from energy consumption real-time monitoring data and/or monthly energy consumption survey data of water plants, desalination water plants, water delivery pumping stations, booster pumping stations and sewage treatment plants in the nuclear region.
According to some preferred embodiments of the present invention, the irrigation area is obtained by remote sensing image inversion and/or measured data of soil moisture sites.
According to some preferred embodiments of the present invention, the thermal terminal consumption of the industry, the building industry and the third industry is obtained through energy consumption real-time monitoring data and/or monthly energy consumption survey data of each energy enterprise and public institution in the checked area.
According to some preferred embodiments of the present invention, the data of the nuclear region population is obtained by a regional annual population dynamic update data survey.
According to some preferred embodiments of the present invention, the consumption of the residential heating power terminal is obtained by real-time monitoring data of a heating department.
According to some preferred embodiments of the invention, the average air temperature is obtained from daily monitoring data at a weather monitoring site.
According to some preferred embodiments of the invention, the per-person drinking water amount and the per-person bathing water amount are taken according to living habits of residents in the checked area.
According to some preferred embodiments of the invention, the obtaining of the irrigation area comprises:
s121, based on the obtained remote sensing image data, obtaining improved soil vertical drought indexes of pixels in different periods in a nucleated region through the following calculation model:
Figure BDA0003773638510000041
Figure BDA0003773638510000042
Figure BDA0003773638510000043
Figure BDA0003773638510000044
wherein MPDI is the improved soil vertical drought index, PDI is the soil vertical drought index, R and NIR respectively represent the spectral reflectivity of pixel red light and near red light in remote sensing image data, m represents the soil line slope, fvc represents the vegetation coverage, NDVI is the normalized vegetation index, NDVI max 、NDVI min Respectively representing the maximum value and the minimum value of the NDVI in the nucleated region;
s122, based on the improved soil vertical drought index, obtaining the soil water content of each pixel of the image of the nucleated region through the following calculation model to obtain different soil water content pixels:
SMC＝λMPDI η +e ε (7)
wherein SMC is the water content of the soil, lambda and eta are regression fitting parameters respectively, e is an exponential function, e ε Representing a random error; the regression fitting is obtained by establishing a regression model between the MPDI in the rainfall-free period in the remote sensing image data and the soil moisture content data actually measured by the soil moisture content site;
s123, judging whether the soil water content pixel is an irrigation pixel:
SMC i,start -SMC i,t ＞T
wherein T represents a soil moisture content threshold, SMC i,t The soil water content of the soil water content pixel i at the time t can be determined according to the soil water content and the soil water content at the time when the soil water content is obviously increased SMC i,start Correspondingly, the obtained pixel with the obviously increased soil water content is an irrigation pixel;
the soil water content threshold value can be determined according to the change condition of the soil water content after the perennial irrigation event recorded by the soil moisture content site of the nucleated area, and the minimum value of the change quantity can be taken;
s124, based on the obtained soil water content pixel and the irrigation pixel, obtaining an irrigation area through the following calculation model:
Area start,t ＝S p ×Num start,t (8)
in the formula, area start,t Indicates the irrigation area from the time t to the time when the water content of the soil is obviously increased, S p Indicates the total area, num, of the soil water content pixel start,t And the number of irrigation pixels from the time t to the time when the soil water content is obviously increased is shown.
According to some preferred embodiments of the invention, said S2 comprises:
s21, carrying out agricultural land carbon emission accounting through the following calculation model:
C d ＝C f +C p +C m +C e (12)
C f ＝αG f d (13)
C p ＝αG p f (14)
C m ＝αρ m s m hpg (15)
C e ＝αbS e +αcW e (16)
Figure BDA0003773638510000051
wherein, C d For ploughing CO 2 Discharge amount, C f Application of CO to fertilizers 2 Discharge amount, C p Using CO for pesticides 2 Discharge amount, C m Use of CO for agricultural films 2 Discharge amount, C e Use of CO for tillage and agricultural machinery 2 Discharging amount; c n For CO of other agricultural land 2 Discharge capacity; g f D =0.85754kgC/kg, representing the carbon emission coefficient of unit fertilizer application amount, in terms of the applied amount of the fertilizer; g p F =4.9341kgC/kg for the amount of pesticide used, representing the carbon emission coefficient per unit amount of pesticide used; ρ is a unit of a gradient m Is the agricultural film specific gravity s m The agricultural film coverage area is h, the agricultural film thickness is h, p is the theoretical coverage of the agricultural film, g =0.00384tC/mg, and the carbon emission coefficient of unit agricultural film usage is expressed; s. the e To the area of cultivated land, W e B =16.47kgC/hm for the total power of the agricultural machine 2 C =0.18kgC/kW, respectivelyCarbon emission coefficient of agricultural tillage and carbon emission coefficient of agricultural machinery; g k Denotes the number of livestock and poultry raised k, lambda k The carbon emission coefficient of the kth livestock and poultry is shown, k =1,2,3,4,5 represents 5 main livestock and poultry such as pigs, cattle, sheep, donkeys and poultry respectively;
s22, carrying out construction land carbon emission accounting through the following calculation model:
Figure BDA0003773638510000061
Figure BDA0003773638510000062
wherein, C h 、C g CO respectively used as traffic land, residential site and industrial and mining land 2 Discharge capacity; eg m 、Ed m The gasoline and diesel oil terminal consumption quantity is m =1,2, which respectively represents the life of residents or a third industry; NCV g 、NCV d The low heating value is gasoline and diesel oil; delta g 、δ d The carbon content of gasoline and diesel oil; OR (OR) g 、OR d The combustion oxidation rate of gasoline and diesel oil; e i The terminal consumption amount of the ith fossil energy is the terminal consumption amount of other fossil fuel varieties except gasoline and diesel consumed by the third industry and residents; NCV i Is the low-level calorific value, delta, of the ith fossil energy i Carbon content of the ith fossil energy, OR i Combustion oxidation rate of i-th fossil energy, EE j J =1,2,3,4 for power terminal consumption, respectively representing residential life, industry, construction industry, third industry, β Electric power To electric power CO 2 The discharge coefficient;
s23, obtaining the carbon sink amount of the carbon sink through the following model:
T i ＝αS i β i ，i＝1,2...6 (20)
wherein, T i 、S i 、β i Carbon sink amount, area and carbon sink coefficient of the land which are the ith land utilization mode, i6, namely 6 land utilization types of cultivated land, garden land, forest land, grassland, water area and water conservancy facility land and unused land respectively;
s24, obtaining the total amount of the land utilization carbon emission through the following calculation model:
Figure BDA0003773638510000063
according to some preferred embodiments of the present invention, the fertilizer application amount, the pesticide usage amount are obtained by real-time monitoring data of pesticide fertilizer usage monitoring points and/or field investigation data of monthly degrees.
According to some preferred embodiments of the present invention, the agricultural film coverage area and the cultivated land area are obtained by interpreting land use classification data through remote sensing images.
According to some preferred embodiments of the invention, the total power of the agricultural machine is obtained by monthly field survey data.
According to some preferred embodiments of the present invention, the consumption amount of the fossil energy terminal and the consumption amount of the power terminal are obtained by energy consumption real-time monitoring data of each enterprise and public institution and each residential user and/or monthly energy consumption survey data.
According to some preferred embodiments of the present invention, the area of the ith land utilization mode land is obtained by interpreting land utilization classification data from remote sensing images.
According to some preferred embodiments of the invention, the obtaining of the land use classification data comprises:
s211, collecting remote sensing image data in the checked area, and carrying out preprocessing of radiometric calibration, atmospheric correction, wave band synthesis, image mosaic and image cutting on the image data to obtain a plurality of image samples;
s212, respectively carrying out land use type labeling on the image samples, and forming a training set by the labeled image samples;
s213, carrying out image supervision and classification training through the training set to obtain a trained classification model, carrying out land utilization type classification on the preprocessed remote sensing image to be interpreted through the classification model, and counting the area of land corresponding to the utilization type according to the classification result to obtain the land utilization classification data.
According to the above method for accounting for carbon emissions in development and utilization of water and soil resources, there can be further obtained an accounting apparatus comprising: the system comprises a data acquisition module for acquiring original data required by the carbon emission accounting method in the water and soil resource development and utilization, a data processing module for processing the acquired original data to acquire calculation data which can be directly utilized by the carbon emission accounting method in the water and soil resource development and utilization, a water resource carbon emission accounting module for accounting the total carbon emission in the water resource development and utilization, a land resource carbon emission accounting module for accounting the total carbon emission in the land resource development and utilization, and a water and soil resource carbon emission total accounting module for performing the total accounting on the total carbon emission in the water and soil resource development and utilization.
According to the above method for accounting for carbon emission in water and soil resource development and utilization, there can be further obtained an accounting apparatus including a memory storing a readable program for implementing the above method for accounting for carbon emission in water and soil resource development and utilization, and a processor operable to execute the program.
According to the method for accounting for carbon emission in development and utilization of water and soil resources, a computer-readable storage medium storing a program and/or an algorithm for implementing the method for accounting for carbon emission in development and utilization of water and soil resources can be further obtained.
The invention has the following beneficial effects:
(1) The accounting method of the invention considers the coupling effect of water and soil resources in addition to the carbon accounting of each link of development and utilization of water resources and land resources, and the accounting result is more accurate and complete.
(2) The accounting method can rely on the prior multi-source data such as monitoring data, survey data, remote sensing images and the like which can be obtained in real time, thereby realizing the monthly, quarterly and annual accounting of carbon emission, overcoming the defect that the traditional carbon emission accounting mostly depends on annual statistical data, and leading the accounting result to be more time-efficient.
(3) The accounting method of the invention provides clear and reasonable accounting steps and an accurate and reliable data acquisition mode, and has strong practicability and operability.
(4) The accounting method or the accounting device, the terminal equipment and the like can automatically acquire and calculate the classification accounting result or the total accounting result, thereby realizing the automation and the intellectualization of the carbon emission accounting for the water and soil resource utilization.
Drawings
FIG. 1 is a flowchart of an embodiment of the accounting method of the present invention.
Fig. 2 is a structure of an embodiment of the accounting apparatus of the present invention.
Detailed Description
The present invention is described in detail below with reference to specific embodiments and examples, but it should be understood that the embodiments and examples are only for illustrative purposes and are not intended to limit the scope of the present invention. All reasonable variations and combinations that fall within the spirit of the invention are intended to be within the scope of the invention.
Referring to fig. 1, according to the technical solution of the present invention, some embodiments of the accounting method for carbon emission from water and soil resources include the following steps:
s1, accounting is carried out on carbon emission in water resource development and utilization
In some embodiments, S1 further comprises:
s11, carbon emission accounting is carried out on water taking, water production, water delivery and water treatment links in water resource development and utilization
Carbon emission in water resource development and utilization mainly comes from energy consumption in links of water taking, water making, water conveying, water using, water treatment and the like, such as energy required by operation of water pumps, pressure pumps, thermal equipment or sewage treatment equipment and the like used in all links.
In consideration of the relationship between energy consumption and carbon emission, in a specific embodiment, the present invention utilizes energy consumption, an energy combustion calorific value, an energy carbon content, and an oxidation rate to perform accounting of carbon emission associated with energy consumption.
Preferably, in order to avoid repeated calculation such as primary energy conversion, the carbon emission amount is calculated by using the terminal consumption amount of energy.
Further, in some embodiments, the carbon emission accounting in the water intake, water production, water delivery and water treatment stages is performed by the following calculation model:
Figure BDA0003773638510000091
wherein, W h CO generated for energy consumption in water taking, water making, water delivery and water treatment links 2 Total amount of emissions; α =44/12 for conversion of carbon to CO 2 The conversion coefficient of (2); e i The terminal consumption amount of the ith fossil energy source is the terminal consumption variety of the unit energy source in water taking, water making, water conveying and water treatment, and the terminal consumption variety comprises coal, gasoline, diesel oil and the like; NCV i Refers to the low calorific value, delta, of the ith fossil energy i Refers to the carbon content, OR, of the ith fossil energy i Refers to the combustion oxidation rate of the ith fossil energy; EH j Beta is the terminal consumption per unit of electric power or heat in water taking, water making, water delivery and water treatment j As electric or thermal CO 2 Discharge coefficient, where j =1 represents thermal power and j =2 represents electrical power.
The energy consumption of fossil energy, electric power and heating power terminals in water taking, water production, water delivery and water treatment can be energy consumption real-time monitoring data acquired from water affair departments and statistics, modification and power departments respectively or energy consumption investigation data counted periodically according to months.
S12, carrying out carbon emission accounting on water consumption ring joints in water resource development and utilization
The water consumption link is divided into first industrial water, second industrial water, third industrial water and resident domestic water according to water consumption departments.
Wherein the first industrial water mainly refers to agricultural water. Agricultural water is mainly used for agricultural production and is irrigated, and irrigation process diversion needs to use equipment such as water pump, and the equipment is driped irrigation to sprinkling irrigation, driping irrigation etc. need use and spout, can consume the energy during equipment use, and then produces the carbon and discharge.
In some embodiments, wherein the carbon emissions in the irrigation water are obtained by the following computational model:
W a ＝αG i h(2)
in the formula, W a CO for irrigation water 2 Discharge amount, G i Is the irrigation area; h =266.48kgC/hm 2 And represents the carbon emission coefficient per unit irrigation area.
The irrigation area can be obtained through remote sensing image inversion or actual measurement data of regional soil moisture content sites, and if the real-time irrigation area data are accumulated according to months, the monthly irrigation area is obtained.
The second industrial water includes industrial and construction water. The industrial water process comprises the steps of water taking, water making, water conveying, water using and water treatment. The building industrial water is mainly construction water and domestic water for constructors in the building construction process.
The third-generation industrial water includes water used in wholesale and retail industries, lodging and catering industries, warehousing and postal industries, and other service industries. The water using process comprises the steps of water taking, water delivery, water using, water treatment and the like, equipment which is put into the water using process comprises various water pumps, pressure pumps, heating power and sewage treatment facilities and the like, and the equipment and the facilities use the energy which is put into the equipment, such as fossil energy, electric power, heating power and the like.
The carbon emission in the links of water taking, water making, water conveying and water treatment of the second industry and the third industry is calculated in the S1 part of the invention.
In the second industrial water and the water consumption link in the industrial water process, the heat takes water as a carrier, or directly participates in the reaction, or serves as energy supply of the reaction and the process, and plays an important role in industrial production, and the part of energy is very huge, so that the carbon emission of industrial heat consumption is preferably used as the carbon emission of the industrial water link, and similarly, the building industrial water in the second industrial water and the water consumption link in the third industry are both accounted by the carbon emission of heat consumption.
Preferably, the carbon emission of the water consumption links of the second industry and the third industry is obtained by the following calculation model:
W e ＝W g +W c ＝E g β heat generation +E c β Heat generation (9)
W s ＝E s β Heat generation (10)
Wherein, W e 、W g 、W c 、W s CO of water consumption links of second, industrial, construction and third industries respectively 2 Discharge amount of, wherein E g 、E c 、E s Thermal terminal consumption, beta, for industry, construction, and third industries, respectively Heat generation Is thermal CO 2 The discharge coefficient.
The consumption of the heating power terminals in the second and third industries can be obtained through energy consumption real-time monitoring data of each energy-using enterprise and public institution or energy consumption investigation data counted regularly according to months.
Besides the first to third industrial water, residents also consume a large amount of water and soil resources and energy resources in life, and correspondingly generate carbon emission.
The domestic water is mainly used for drinking, bathing, washing, heating power service and the like in daily life, and the energy input in the water is mainly electric power and heating power and is used for supporting the domestic water supply of residents for pressurization, heating and the like.
Preferably, the invention uses the electricity and heat consumed by the residents using hot water and heat to calculate the carbon emission of the resident domestic water, and the calculation model is as follows:
Figure BDA0003773638510000111
wherein, W b CO as domestic water for residents 2 Emission amount, beta Heat generation Is thermal CO 2 Discharge coefficient, p is the density of water, Q Drinking water Water consumption for human body, P is regional personThe number of mouths, T being the annual average temperature, E r Heat terminal consumption for resident life, Q Washing machine Water consumption for bathing T i ' is the average air temperature in each season, i =1,2,3,4, which respectively represents spring, summer, autumn and winter seasons.
The embodiment obtains monthly heat supply data from a certain city development and reform department, and the average temperature can be obtained through daily monitoring data of a meteorological monitoring site; the water consumption for drinking and bathing can be calculated according to the living habits of residents in the checked area.
S13 total carbon emission C for water resource utilization through the following calculation model w And (4) carrying out accounting:
C W ＝W h +W a +W e +W s +W b (11-2)
s2, accounting is carried out on carbon emission in land resource development and utilization
The land utilization types comprise agricultural land, construction land and unused land, wherein the agricultural land comprises cultivated land, garden land, forest land and grassland; the construction land is divided into residential areas, industrial and mining land, traffic land, water areas and water conservancy facilities land. The carbon source land comprises cultivated land, residential site, industrial and mining land and transportation land, and the carbon sink land comprises cultivated land, forest land, garden land, grassland, water area, water conservancy facility land and unused land.
In some embodiments, S2 further comprises:
s21 agricultural carbon emission accounting
The agricultural land resource development and utilization activities mainly comprise activities such as land leveling, farming, harvesting and the like, and the carbon source is cultivated land. The cultivation of the farmland needs to put in chemical fertilizers, pesticides, agricultural films and other substances, and the production and the use of the substances generate carbon emission; agricultural tillage destroys organic carbon reservoirs in soil to cause carbon emission; agricultural machinery is needed during land leveling, ploughing and harvesting, and the agricultural machinery consumes fuels such as diesel oil, gasoline and electric power to generate carbon emission.
Preferably, the method performs the agricultural land carbon emission accounting through the following calculation model:
C d ＝C f +C p +C m +C e (12)
C f ＝αG f d (13)
C p ＝αG p f (14)
C m ＝αρ m s m hpg (15)
C e ＝αbS e +αcW e (16)
Figure BDA0003773638510000121
wherein, C d For ploughing CO 2 Discharge amount, C f Application of CO to fertilizers 2 Discharge amount, C p Using CO for pesticides 2 Discharge amount, C m Using CO for agricultural films 2 Discharge amount, C e Use of CO for tillage and agricultural machinery 2 Discharge amount, C n For CO from other agricultural land 2 Discharge capacity; g f D =0.85754kgC/kg, representing the carbon emission coefficient of unit fertilizer application amount; g p F =4.9341kgC/kg for pesticide usage, representing the carbon emission coefficient per pesticide usage; ρ is a unit of a gradient m Is the agricultural film specific gravity s m The agricultural film coverage area is h, the agricultural film thickness is h, p is the theoretical coverage of the agricultural film, g =0.00384tC/mg, and the carbon emission coefficient of unit agricultural film usage is expressed; s. the e To the area of cultivated land, W e B =16.47kgC/hm for the total power of the agricultural machine 2 C =0.18kgC/kW, which represents the carbon emission coefficient for agricultural tillage and the carbon emission coefficient for agricultural machinery, respectively, G k Denotes the number of livestock and poultry raised k, lambda k The carbon emission coefficient of the kth livestock and poultry is shown, and k =1,2,3,4,5 represents 5 main livestock and poultry such as pigs, cattle, sheep, donkeys and poultry.
The application amount of the fertilizer and the usage amount of the pesticide can be obtained through real-time monitoring data of pesticide and fertilizer usage monitoring points or field investigation data counted at regular intervals according to months; the agricultural film coverage area and the cultivated land area can be obtained through land utilization classification data obtained through remote sensing image interpretation; the total power of the agricultural machine can be obtained through field investigation data counted according to a monthly period.
S22 construction land carbon emission accounting
The carbon emissions from construction sites are mainly those generated by industrial and human activities carried thereon. The carbon emission of the transportation land mainly comprises carbon emission generated by the use of third industry and resident living transportation means, and the energy consumption varieties mainly comprise gasoline and diesel oil consumed by various vehicles; the carbon emission of residential sites and industrial and mining sites comprises carbon emission generated by terminal energy consumption of residential life, industry, building industry and third industry, and the energy consumption varieties comprise various fossil fuels, electric power and the like.
Preferably, the invention uses the energy consumption of various land-bearing industrial and residential life terminals to calculate the carbon emission of the construction land, and the calculation model is as follows:
Figure BDA0003773638510000131
Figure BDA0003773638510000132
wherein, C h 、C g CO respectively used as traffic land, residential site and industrial and mining land 2 Emission amount, eg m 、Ed m For gasoline and diesel end consumption, m =1,2, respectively for residential life or third industry, NCV g 、NCV d Is gasoline and diesel oil low heating value delta g 、δ d Is the carbon content of gasoline and diesel oil, OR g 、OR d The combustion oxidation rate of gasoline and diesel oil, E i For the ith fossil energy terminal consumption, the calculated energy variety is gasoline and diesel oil which are consumed by the residents and consumed by the third industry and the residents in the resident life, industry, construction industry and the third industrial fossil energy terminal energy varietyOther fossil energy varieties; NCV i Is the low calorific value delta of the ith fossil energy i Carbon content of the ith fossil energy, OR i Combustion oxidation rate of i-th fossil energy, EE j J =1,2,3,4 represents the residential life, industry, construction industry, third industry, respectively, for the consumption of the power terminal, β Electric power To electric power CO 2 The discharge coefficient.
The consumption of the second industry, the third industry and the resident bioactivation stone fuel and the electric power terminal can be monitored in real time through energy consumption of each enterprise and public institution and each resident user or energy consumption survey data counted regularly according to months.
S23 carbon sink amount accounting for carbon sink
The carbon sink land comprises cultivated land, forest land, garden land, grassland, water area and water conservancy facilities land and unused land.
Preferably, the carbon sequestration capacity of the carbon sequestration site is calculated by adopting an emission coefficient method, and a calculation model is as follows:
T i ＝αS i β i ，i＝1,2...6 (20)
wherein, T i 、S i 、β i The carbon sequestration amount, the area and the carbon sequestration coefficient of the ith land utilization mode are respectively, i =1,2.. 6 is respectively 6 land utilization types of cultivated land, forest land, garden land, grassland, water area, water conservancy facility land and unused land, and the carbon sequestration coefficient of each land utilization is as follows: cultivated land value is 0.042tC/hm 2 A, forest land and garden take values of 0.58tC/hm 2 A, meadow value 0.021tC/hm 2 A, the land area value of the water area and the water conservancy facilities is 0.257tC/hm 2 A, unused value 0.005tC/hm 2 And a, converting the annual carbon sink coefficient into a monthly value in calculating the monthly data to participate in calculation.
The required cultivated land, forest land, garden land, grassland, water area and land area not used can be obtained by the land use classification data obtained by remote sensing image interpretation.
S24 total carbon emission accounting for land utilization
The total carbon emission amount of the land for land use is obtained by the following calculation model, wherein the total carbon emission amount of the land for land use is obtained by subtracting the difference of the carbon sink amount of the carbon sink land from the sum of the carbon emission amount of the land for land use and the carbon emission amount of the construction land:
Figure BDA0003773638510000141
s3, carrying out total accounting on carbon emission in water and soil resource development and utilization
Preferably, the invention obtains the carbon emission in the development and utilization of the water and soil resources through the following calculation model:
C net ＝C W +C L (22)
wherein, C net The total carbon emission is utilized for water and soil resources.
Referring to fig. 2, according to the method for accounting carbon emission in the development and utilization of water and soil resources, an accounting device can be further obtained, which comprises:
a data acquisition module: and acquiring various data required in the calculation of the carbon emission of the water and soil resources of the checked area, wherein the various data comprise monitoring data, survey data, remote sensing image data and the like.
A data processing module: and processing the acquired data to obtain numerical data which can be directly utilized by accounting the carbon emission of the water and soil resources.
The water resource carbon emission accounting module: and carrying out accounting on the carbon emission in the water resource development and utilization, wherein the accounting comprises the carbon emission accounting in the water taking, water making, water conveying and water treatment links, the carbon emission accounting in the water using link and the total carbon emission accounting in the water resource utilization.
The land resource carbon emission accounting module comprises: and carrying out accounting on the carbon emission in the land resource development and utilization, including the carbon emission of agricultural land, the carbon emission of construction land, the carbon sink amount and the total carbon emission of land resources.
And (3) accounting the total carbon emission of water and soil resources: and the module for carrying out the total accounting on the carbon emission in the development and utilization of the water and soil resources comprises the accounting on the total carbon emission in the utilization of the water and soil resources.
According to the above method for accounting for carbon emission in water and soil resource development and utilization, there can be further obtained an accounting apparatus including a memory in which a program for implementing the above method for accounting for carbon emission in water and soil resource development and utilization is stored and a processor which can run the program.
Further, the apparatus may include: input/output device, network access device, bus, etc
In some embodiments, the device may be a desktop computer, a laptop computer, a cloud server, or the like.
According to the method for accounting for carbon emission in development and utilization of water and soil resources, a computer-readable storage medium storing a program and/or an algorithm for implementing the method for accounting for carbon emission in development and utilization of water and soil resources can be further obtained.
In some embodiments, the storage medium may be, for example, ROM/RAM, a magnetic disk, an optical disk, or the like.
Example 1
The carbon emission of water and soil resources at one month of a certain city is accounted through the specific implementation mode, wherein in S11:
the main energy varieties in the collected terminal consumption comprise coal, gasoline, diesel oil, electric power and heat, and the carbon emission is calculated according to the terminal consumption of the energy varieties.
And carrying out values as shown in table 1 on the combustion calorific value, the carbon content and the carbon oxidation rate of the fossil energy.
To power CO 2 The emission coefficient is estimated by using the standard coal consumption of power supply in a certain year, and the value is 736.53gCO 2 /KWh。
Thermal CO 2 The emission coefficient is 850.725tCO 2 /10TJ。
TABLE 1 fossil energy calorific value, carbon content and oxidation rate table
Figure BDA0003773638510000151
Figure BDA0003773638510000161
Energy consumption real-time monitoring data of various water plants, seawater desalination plants, water delivery pump stations, booster pump stations and sewage treatment plants in a certain market or energy consumption survey data counted periodically according to months are respectively obtained from a water affair part and a statistics, modification and power department.
In S12, an irrigation area selection improved vertical drought index (MPDI) is obtained through a remote sensing image inversion area, and soil water content inversion is carried out by taking optical data as a data source, wherein the inversion comprises the following steps:
firstly, calculating the vertical drought index (MPDI) of each pixel in each month of a certain city:
Figure BDA0003773638510000162
Figure BDA0003773638510000163
Figure BDA0003773638510000164
Figure BDA0003773638510000165
wherein MPDI is an improved soil vertical drought index for representing the water content of soil, PDI is a soil vertical drought index, R and NIR respectively represent the spectral reflectivity of pixel red light and near red light, m represents the slope of a soil line, fvc represents the vegetation coverage, NDVI is a normalized vegetation index, NDVI max 、NDVI min Respectively representing the maximum and minimum values of the region NDVI.
Secondly, in order to remove the influence of rainfall on the soil water content, selecting data of a rainfall-free period, calculating an MPDI value of the period, establishing a regression model with data of the soil water content actually measured by 19 soil moisture content sites in a certain market at the same time, and obtaining the soil water content of each pixel of each period of the region through inversion as follows:
SMC＝λMPDI η +e ε (7)
wherein SMC is the water content of the soil, lambda and eta are regression fitting parameters respectively, e is an exponential function, e ε Representing a random error; the regression fitting is obtained by establishing a regression model between the MPDI in the rainfall-free period in the remote sensing image data and the soil moisture content data actually measured by the soil moisture content site;
the method for judging the rainfall-free condition comprises the following steps:
P＞X
wherein, P is interpolation data of rainfall capacity of the regional meteorological station, X is a threshold value for judging no precipitation event, and the data is set after comprehensive judgment according to the historical record of precipitation event in the local soil moisture content data, and is set to be 2mm in an embodiment of a certain city.
Thirdly, judging whether the soil water content pixel is an irrigation pixel:
SMC i,start -SMC i,t ＞T
wherein the SMC i,start 、SMC i,t And the soil water content of the soil water content pixel i at the moment T and the moment when the soil water content is obviously increased by the start is shown as T, and the T represents the soil water content threshold value.
Finally, the area of the city irrigated every month in a certain year is calculated according to the following formula:
Area start,t ＝S p ×Num start,t (8)
Area start,t showing that the irrigation area from the time t to the time when the soil water content is obviously increased by the start time S p Indicates the area of the pixel element of the water content of the soil, num start,t And showing that the quantity of irrigation pixels from the time t to the time when the soil water content is obviously increased by the start time.
The remote sensing data required in the inversion calculation of the irrigation area can be obtained through medium-high resolution remote sensing image data meeting the requirements of spatial resolution, temporal resolution and spectral resolution, such as GF-1WFV, landsat8OLI data and HJ-1A/1BCCD data, and Landsat8OLI data of each month in a certain city is selected in the embodiment; the soil moisture content data can be obtained through actually measured data of soil moisture content sites, in the embodiment, the data of the soil moisture content sites in a certain market 19 is selected and comes from soil moisture content monitoring reports issued by the urban agricultural center; the rainfall data can be obtained through actually measured data of a weather station, and data of a weather station in a city 13 is selected in the embodiment.
The required population data is obtained by regional annual population dynamic update data investigation, the embodiment is obtained from a statistical department in a certain city, the consumption of the resident heating power terminal can be obtained by real-time monitoring data of a heating department, the embodiment obtains monthly heating data from a development and reform department in the certain city, the average temperature can be obtained by daily monitoring data of a meteorological monitoring site, and the embodiment obtains the average value from a meteorological site at a meteorological 13 position in the certain city; the water consumption for daily drinking, water consumption for daily bathing and daily number of daily bathing can be selected according to the living habits of residents in the nuclear area, embodiment Q Drinking water Value 2L/person.day, Q Washing machine The value is 32L/person-day, ts i The bathing days in spring and autumn are set to be 39 days, 91 days in summer and 26 days in winter. S21, the fertilizer application amount and the pesticide usage amount can be obtained through real-time monitoring data of pesticide and fertilizer usage monitoring points or field investigation data counted periodically monthly, and the embodiment is obtained through a soil moisture content monitoring report of a soil moisture content site at 19 places of a certain city; the agricultural film coverage area and the cultivated land area can be obtained through land utilization classification data obtained through remote sensing image interpretation; the total power of agricultural machinery can be obtained by field investigation data counted periodically according to months, and the embodiment is obtained from agricultural parts and statistical departments in a certain market.
The manner of acquiring land use classification data through remote sensing influence interpretation is as follows:
firstly, remote sensing image data of a certain city month by month is collected, and preprocessing such as radiometric calibration, atmospheric correction, wave band synthesis, image mosaic, image cutting and the like is carried out on the image to obtain an image sample of each month in a certain year.
Secondly, according to the obtained image samples, defining training samples according to the classification standard of the current land utilization state (GB/T21010-2007), supervising and classifying the images, dividing land utilization types of a certain city into cultivated land (containing greenhouse films), garden land, forest land, grassland, residential points, industrial and mining land, transportation land, water area, water conservancy facility land and unused land, performing classification precision verification by using a Kappa coefficient, wherein the Kappa coefficient of the embodiment is 0.92, and passing precision verification.
And finally, fusing the classified land utilization type data according to the types, and counting various land areas to obtain various land areas such as agricultural film coverage area, cultivated land area and the like in each month of a year.
The remote sensing image data can select medium-high resolution remote sensing image data which can simultaneously meet the requirements of spatial resolution, temporal resolution and spectral resolution, such as GF-1WFV/PMS and Landsat8OLI data, and the embodiment selects Landsat8OLI image monthly-by-month data in a certain market.
In S22, the energy varieties of the embodiment are raw coal, cleaned coal, other cleaned coal, molded coal, coke, crude oil, gasoline (only consumption of industry and construction industry is calculated), diesel oil (only consumption of industry and construction industry is calculated), kerosene, fuel oil, liquefied petroleum gas and natural gas
The consumption of the second industry, the third industry and the resident activated stone fuel and the power terminal can be monitored in real time through energy consumption of each enterprise and public institution and each resident user or energy consumption survey data counted regularly according to months, and the embodiment acquires the energy consumption data from a certain department of marketing and modification and a statistical department.
In S23, the required cultivated land, garden land, forest land, grassland, water area, and unused land area may be obtained by the land use classification data obtained by remote sensing image interpretation, and the embodiment is obtained by statistics of the land use classification data of the monthly land use in a certain city of land use year.
In S24, the embodiment obtains the total amount of the carbon emission of the land utilization in each month of a certain city in the embodiment, and obtains the total amount of the carbon emission of the water and soil resource utilization in each month of the certain city in the embodiment through S3. Firstly, the obtained result realizes the carbon emission accounting in the municipal water and soil resource utilization process, defines the carbon emission in each municipal water and soil resource utilization process, and visually embodies the carbon emission effect of the municipal water and soil resource and energy coupling. And secondly, monthly data of the carbon emission of the municipal water and soil resources are obtained, compared with other regional carbon emission accounting methods, only annual data can be calculated, and the shortage of annual dynamic data is difficult to obtain.
The above examples are merely preferred embodiments of the present invention, and the scope of the present invention is not limited to the above examples. All technical schemes belonging to the idea of the invention belong to the protection scope of the invention. It should be noted that modifications and adaptations to those skilled in the art without departing from the principles of the present invention should also be considered as within the scope of the present invention.

Claims (10)

1. A method for accounting carbon emission in development and utilization of water and soil resources is characterized by comprising the following steps:
s1, accounting is carried out on carbon emission of different links in water resource development and utilization, and the total carbon emission in the water resource development and utilization is obtained through the sum of the carbon emission of the different links;
s2, the carbon emission amount of different carbon source land and the carbon sink amount of different carbon sink land in the land resource development and utilization are calculated, and the total carbon emission amount in the land resource development and utilization is obtained through the difference between the sum of the carbon emission amount of the different carbon source land and the sum of the carbon sink amount of the different carbon sink land;
s3, carrying out total accounting on carbon emission in water and soil resource development and utilization, wherein the total accounting model is as follows:
C net ＝C W +C L (22)
wherein, C net Total carbon emission for water and soil resource development and utilization, C w Represents the total carbon emission, C, in the development and utilization of water resources L Representing a total amount of carbon emissions in the land resource development and utilization;
wherein, different links in the development and utilization of water resources comprise water taking, water production, water delivery, water treatment links and water utilization links, and the water utilization links comprise first to third industrial water utilization links and a resident domestic water utilization link; the carbon source land comprises cultivated land, other agricultural land except the cultivated land, residential site, industrial and mining land and transportation land; the carbon sink lands include cultivated lands, woodlands, gardens, grasslands, water areas and water conservancy facilities, and unused lands.
2. The carbon emission amount accounting method according to claim 1, wherein the S1 includes:
s11, carrying out carbon emission accounting on water taking, water production, water delivery and water treatment links in water resource development and utilization through the following calculation model:
Figure FDA0003773638500000011
wherein, W h CO generated by energy consumption in the water taking, water making, water conveying and water treatment links 2 Discharge capacity; α =44/12 for conversion of carbon to CO 2 The conversion coefficient of (2); e i Terminal consumption of the ith fossil energy, NCV i Is the low calorific value, delta, of the ith fossil energy i Is the carbon content of the ith fossil energy, OR i Combustion oxidation rate for the ith fossil energy source; EH j Beta is the terminal consumption per unit of electric power or heat in water taking, water making, water delivery and water treatment j As electric or thermal CO 2 A discharge coefficient, where j =1 represents thermal power, and j =2 represents electric power;
s12, carrying out carbon emission accounting on the water consumption ring section in water resource development and utilization through the following calculation model:
the first industrial water link carbon emission calculation model comprises:
W a ＝αG i h (2)
wherein, W a CO for irrigation water 2 Discharge amount, G i Is the irrigation area; h is the carbon emission coefficient of unit irrigation area;
a second industry and a third industry water link carbon emission calculation model:
W e ＝W g +W c ＝E g β heat generation +E c β Heat generation (9)
W s ＝E s β Heat generation (10)
Wherein, W e 、W g 、W c 、W s CO of water consumption links of second, industrial, construction and third industries respectively 2 Discharge amount of, wherein E g 、E c 、E s Thermal terminal consumption, beta, of industry, building, and third industries, respectively Heat generation Is thermal CO 2 The discharge coefficient;
a carbon emission calculation model in a resident domestic water link:
Figure FDA0003773638500000021
wherein, W b CO as domestic water for residents 2 Emission amount, beta Heat generation Is thermal CO 2 Discharge coefficient, p is the density of water, Q Drinking water The water consumption for human body, P is the population number of the nuclear region, T is the average annual temperature, E r Heat terminal consumption for resident life, Q Washing machine Water consumption for bathing T i ' is the average air temperature in each season, i =1,2,3,4, which respectively represents spring, summer, autumn, winter;
s13, developing and utilizing the total carbon emission C in the water resource by the following calculation model w And (4) performing accounting:
C W ＝W h +W a +W e +W s +W b (11-2)。
3. the method for accounting for carbon emission according to claim 2, wherein the terminal consumption of fossil energy and the terminal consumption of electricity and heat in the water intake, production, delivery and treatment links are obtained from energy consumption real-time monitoring data and/or monthly energy consumption survey data of water plants, desalination plants, water delivery pumping stations, booster pumping stations and sewage treatment plants in the nuclear region; and/or the irrigation area is obtained by remote sensing image inversion and/or soil moisture site actual measurement data; and/or the heat power terminal consumption of the industry, the building industry and the third industry is obtained through energy consumption real-time monitoring data of each energy enterprise and public institution in a checked area and/or energy consumption survey data counted monthly; and/or the data of the checked regional population is obtained by dynamic update data survey of regional annual population; and/or the consumption of the resident life heating power terminal is obtained through real-time monitoring data of a heating department, and/or the average temperature is obtained through daily monitoring data of a weather monitoring site; and/or the water consumption for per-person drinking and per-person bathing is taken according to the living habits of residents in the checked area.
4. A carbon emission amount accounting method according to claim 3, wherein the obtaining of the irrigation area comprises:
s121, based on the obtained remote sensing image data, obtaining the improved soil vertical drought index of each pixel in different periods in the nucleated area through the following calculation model:
Figure FDA0003773638500000031
Figure FDA0003773638500000032
Figure FDA0003773638500000033
Figure FDA0003773638500000034
wherein MPDI is the improved soil vertical drought index, PDI is the soil vertical drought index, R and NIR respectively represent the spectral reflectivity of red light and near red light of a pixel in remote sensing image data, m represents the slope of a soil line, fvc represents the vegetation coverage,NDVI is the normalized vegetation index, NDVI max 、NDVI min Respectively representing the maximum value and the minimum value of the NDVI in the nucleated region;
s122, based on the improved soil vertical drought index, obtaining the soil water content of each pixel of the image of the nucleated region through the following calculation model to obtain different soil water content pixels:
SMC＝λMPDI η +e ε (7)
wherein SMC is the soil water content, lambda and eta are regression fitting parameters respectively, e is an exponential function, e ε Representing a random error; the regression fitting is obtained by establishing a regression model between the MPDI in the rainfall-free period in the remote sensing image data and the soil moisture content data actually measured by the soil moisture content site;
s123, judging whether different soil water content pixels are irrigation pixels or not according to the soil water content threshold value, as follows:
when SMC is finished i,start -SMC i,t If the value is more than T, the pixel i is an irrigation pixel;
wherein T represents a soil moisture content threshold, SMC i,t The soil water content of the soil water content pixel i at the time t can be determined according to the soil water content and the soil water content at the time when the soil water content is obviously increased SMC i,start Correspondingly, the obtained pixel with the obviously increased soil water content is an irrigation pixel;
s124, based on the obtained soil water content pixel and the irrigation pixel, obtaining an irrigation area through the following calculation model:
Area start,t ＝S p ×Num start,t (8)
wherein, area start,t Indicates the irrigation area from the time t to the time when the water content of the soil is obviously increased, S p Indicates the total area, num, of the soil water content pixel start,t And the number of irrigation pixels from the time t to the time when the soil water content is obviously increased is shown.
5. The carbon emission amount accounting method according to claim 1, wherein the S2 includes:
s21, carrying out agricultural land carbon emission accounting through the following calculation model:
C d ＝C f +C p +C m +C e (12)
C f ＝αG f d (13)
C p ＝αG p f (14)
C m ＝αρ m s m hpg (15)
C e ＝αbS e +αcW e (16)
Figure FDA0003773638500000041
wherein, C d For ploughing CO 2 Discharge amount, C f Application of CO to fertilizers 2 Discharge amount, C p Using CO for pesticides 2 Discharge amount, C m Use of CO for agricultural films 2 Discharge amount, C e Use of CO for tillage and agricultural machinery 2 Discharge capacity; c n For CO from other agricultural land 2 Discharging amount; g f D =0.85754kgC/kg, representing the carbon emission coefficient of unit fertilizer application amount; g p The dosage of the pesticide is the usage amount; f =4.9341kgC/kg, representing the carbon emission coefficient per unit pesticide usage amount; rho m Is the agricultural film specific gravity; s m The agricultural film coverage area; h is the thickness of the agricultural film; p is the theoretical coverage of the agricultural film; g =0.00384tC/mg, representing the carbon emission coefficient per agricultural film usage; s. the e To the area of cultivated land, W e B =16.47kgC/hm2 and c =0.18kgC/kW of the total power of the agricultural machinery, and respectively represent an agricultural plowing carbon emission coefficient and an agricultural machinery use carbon emission coefficient; g k Denotes the k-th livestock breeding number, lambda k The carbon emission coefficient of the kth livestock and poultry is shown, k =1,2,3,4,5 represents 5 main livestock and poultry such as pigs, cattle, sheep, donkeys and poultry respectively;
s22, carrying out construction land carbon emission accounting through the following calculation model:
Figure FDA0003773638500000051
Figure FDA0003773638500000052
wherein, C h 、C g CO respectively used as traffic land, residential site and industrial and mining land 2 Discharge capacity; eg m 、Ed m The gasoline and diesel oil terminal consumption quantity is m =1,2, which respectively represents the life of residents or a third industry; NCV g 、NCV d The low heating value is gasoline and diesel oil; delta g 、δ d The carbon content of gasoline and diesel oil; OR (OR) g 、OR d The combustion oxidation rate of gasoline and diesel oil; e i The terminal consumption amount of the ith fossil energy is the terminal consumption amount of other fossil fuel varieties except gasoline and diesel consumed by the third industry and residents; NCV i Is the low-level calorific value, delta, of the ith fossil energy i Carbon content of the ith fossil energy, OR i Combustion oxidation rate of i-th fossil energy, EE j J =1,2,3,4 for power terminal consumption, respectively representing residential life, industry, construction industry, third industry, β Electricity To electric power CO 2 A discharge coefficient;
s23, obtaining the carbon sink amount of the carbon sink through the following model:
T i ＝αS i β i ，i＝1,2...6(20)
wherein, T i 、S i 、β i The carbon sequestration amount, the area and the carbon sequestration coefficient of the land which are respectively the ith land utilization mode, i =1,2.. 6, which are respectively 6 land utilization types of cultivated land, garden land, forest land, grassland, water area, water conservancy facility land and unused land;
s24, obtaining the total carbon emission amount in the land resource development and utilization through the following calculation model:
Figure FDA0003773638500000053
wherein, C L Representing the total amount of carbon emission in the land resource development and utilization.
6. The method of claim 5, wherein the fertilizer application amount and the pesticide usage amount are obtained from real-time monitoring data of pesticide and fertilizer usage monitoring points and/or monthly field survey data; and/or the agricultural film coverage area and the cultivated land area are obtained through land utilization classification data obtained through remote sensing image interpretation; and/or, the agricultural machine total power is obtained through monthly field survey data; and/or the fossil energy terminal consumption and the electric power terminal consumption are obtained through energy consumption real-time monitoring data of each enterprise and public institution and each resident user and/or monthly energy consumption survey data; and/or the area of the land in the ith land utilization mode is obtained through land utilization classification data obtained through remote sensing image interpretation.
7. The carbon emission amount accounting method according to claim 5, wherein the obtaining of the land use classification data includes:
s211, collecting remote sensing image data in the checked area, and carrying out preprocessing of radiometric calibration, atmospheric correction, wave band synthesis, image mosaic and image cutting on the image data to obtain a plurality of image samples;
s212, respectively carrying out land use type labeling on the image samples, and forming a training set by the labeled image samples;
s213, carrying out image supervision and classification training through the training set to obtain a trained classification model, carrying out land utilization type classification on the preprocessed remote sensing image to be interpreted through the classification model, and counting the area of land corresponding to the utilization type according to the classification result to obtain the land utilization classification data.
8. A computer-readable storage medium storing a program and/or an algorithm for realizing the method for accounting for carbon emissions in development and utilization of water and soil resources according to any one of claims 1 to 7.
9. An accounting apparatus that implements the method for accounting carbon emissions in development and utilization of water and soil resources according to any one of claims 1 to 7, comprising: the system comprises a data acquisition module for acquiring original data required by the carbon emission accounting method in the water and soil resource development and utilization, a data processing module for processing the acquired original data to acquire calculation data which can be directly utilized by the carbon emission accounting method in the water and soil resource development and utilization, a water resource carbon emission accounting module for accounting the total carbon emission in the water resource development and utilization, a land resource carbon emission accounting module for accounting the total carbon emission in the land resource development and utilization, and a water and soil resource carbon emission total accounting module for performing the total accounting on the total carbon emission in the water and soil resource development and utilization.
10. The apparatus for accounting for carbon emissions in development and utilization of water and soil resources according to any one of claims 1 to 7, comprising a memory storing a readable program for implementing the method for accounting for carbon emissions in development and utilization of water and soil resources according to any one of claims 1 to 7, and a processor operable to execute the program.
CN202210909827.9A 2022-07-29 2022-07-29 Water and soil resource carbon emission accounting method, device, equipment and storage medium Pending CN115204729A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210909827.9A CN115204729A (en) 2022-07-29 2022-07-29 Water and soil resource carbon emission accounting method, device, equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210909827.9A CN115204729A (en) 2022-07-29 2022-07-29 Water and soil resource carbon emission accounting method, device, equipment and storage medium

Publications (1)

Publication Number Publication Date
CN115204729A true CN115204729A (en) 2022-10-18

Family

ID=83586748

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210909827.9A Pending CN115204729A (en) 2022-07-29 2022-07-29 Water and soil resource carbon emission accounting method, device, equipment and storage medium

Country Status (1)

Country Link
CN (1) CN115204729A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116542430A (en) * 2023-07-07 2023-08-04 红杉天枰科技集团有限公司 Multi-dimensional water carbon emission intelligent analysis method and system
CN116975784A (en) * 2023-09-19 2023-10-31 四川省水利科学研究院 High-space-time resolution MPDI data set construction method, system and storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130142401A (en) * 2012-06-19 2013-12-30 건국대학교 산학협력단 Calculation system for reducing carbon emission of u-city and method therefor
CN103870678A (en) * 2014-02-18 2014-06-18 上海零碳建筑科技有限公司 Carbon emission metering method based on city planning
CN109496797A (en) * 2018-12-20 2019-03-22 张正卓 A kind of agricultural water-saving irrigation system and control method
CN114154871A (en) * 2021-12-02 2022-03-08 特斯联科技集团有限公司 Building carbon emission analysis system based on big data processing
CN114169714A (en) * 2021-11-24 2022-03-11 上海市节能减排中心有限公司 Regional carbon emission accounting method, system and computer storage medium
CN114638497A (en) * 2022-03-14 2022-06-17 江南大学 Household carbon emission accounting management method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130142401A (en) * 2012-06-19 2013-12-30 건국대학교 산학협력단 Calculation system for reducing carbon emission of u-city and method therefor
CN103870678A (en) * 2014-02-18 2014-06-18 上海零碳建筑科技有限公司 Carbon emission metering method based on city planning
CN109496797A (en) * 2018-12-20 2019-03-22 张正卓 A kind of agricultural water-saving irrigation system and control method
CN114169714A (en) * 2021-11-24 2022-03-11 上海市节能减排中心有限公司 Regional carbon emission accounting method, system and computer storage medium
CN114154871A (en) * 2021-12-02 2022-03-08 特斯联科技集团有限公司 Building carbon emission analysis system based on big data processing
CN114638497A (en) * 2022-03-14 2022-06-17 江南大学 Household carbon emission accounting management method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
柯丽娟: "遥感反演土壤含水量在灌溉用水管理中的应用研究", 《中国优秀硕士学位论文全文数据库 农业科技辑》, no. 4, pages 21 *
沈静: "遥感技术在灌溉面积监测上的应用研究", 《中国优秀硕士学位论文全文数据库 农业科技辑》, no. 10, pages 13 *
肖宏伟: "中国碳排放测算方法研究", 《阅江学刊》, no. 5, pages 48 *
赵荣钦等: "基于"水—能—碳"关联的城市水系统碳排放研究", 《地理学报》, vol. 76, no. 12, pages 3119 *
陈军娟等: "基于系统动力学的民勤绿洲农业系统碳排放仿真模拟研究", 《西南农业学报》, vol. 35, no. 6, pages 1433 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116542430A (en) * 2023-07-07 2023-08-04 红杉天枰科技集团有限公司 Multi-dimensional water carbon emission intelligent analysis method and system
CN116542430B (en) * 2023-07-07 2024-01-26 红杉天枰科技集团有限公司 Multi-dimensional water carbon emission intelligent analysis method and system
CN116975784A (en) * 2023-09-19 2023-10-31 四川省水利科学研究院 High-space-time resolution MPDI data set construction method, system and storage medium
CN116975784B (en) * 2023-09-19 2023-12-29 四川省水利科学研究院 High-space-time resolution MPDI data set construction method, system and storage medium

Similar Documents

Publication Publication Date Title
Giwa et al. A comprehensive review on biomass and solar energy for sustainable energy generation in Nigeria
CN112765800B (en) Design method of distributed water resource configuration model
CN115204729A (en) Water and soil resource carbon emission accounting method, device, equipment and storage medium
De Girolamo et al. Land use scenario development as a tool for watershed management within the Rio Mannu Basin
Ramachandra et al. Domestic energy consumption patterns in Uttara Kannada district, Karnataka state, India
CN111899126B (en) Three red line control index dividing method based on water circulation simulation
Frisk Simulation and Optimization of a Hybrid Renewable Energy System for application on a Cuban farm
Hu et al. Measuring water footprint on a lake basin scale: a case study of lake dianchi, China
Richards Policy and research implications of recent carbon-sequestering analysis
Roy et al. Multiple impacts of droughts and assessment of drought policy in major drought prone states in India
CN113642842A (en) Department report based ecological system production value business accounting method
Ramachandra et al. Decision support system to assess regional biomass energy potential
Rutibabara Environmental and Economic Cost Analysis of a Solar PV, Diesel and hybrid PV-Diesel water Pumping Systems for Agricultural Irrigation in Rwanda: Case study of Bugesera district
CN114416692B (en) Method for constructing river basin water environment management data resource system
Bansouleh Development of a spatial planning support system for agricultural policy formulation related to land and water resources in Borkhar & Meymeh district, Iran
Kaoma et al. Sustainability performance of lignocellulosic biomass-to-bioenergy supply chains for Rural Growth Centres in Zambia
Painter Estimated use of water in Georgia for 2015 and water-use trends, 1985–2015
Tidwell et al. Biofuel impacts on water.
Zhao et al. Energy-saving techniques in urban aquaponics farms by optimizing equipment operating scheme
CN111079087A (en) Method for measuring and calculating greenhouse gas emission in live pig breeding industry
Jensen et al. Environmental impact of Danish organic tomatoes grown in greenhouses: Quantifying the reduction potential from changes in energy supply towards 2030
Melnikova Assessment of renewable energy potentials based on GIS. A case study in southwest region of Russia
Zhou et al. A review on development water-saving agriculture in Asia
Li et al. Environment‐Friendly Behavior of New Agricultural Business Main Body Based on the Internet of Things
Adusumilli Economic policy and resource implications of biofuel feedstock production

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20221018

RJ01 Rejection of invention patent application after publication