CN117236526A

CN117236526A - River methane annual emission determining method, river methane annual emission determining device, electronic equipment and storage medium

Info

Publication number: CN117236526A
Application number: CN202311500838.2A
Authority: CN
Inventors: 夏星辉; 王君峰
Original assignee: Beijing Normal University
Current assignee: Beijing Normal University
Priority date: 2023-11-13
Filing date: 2023-11-13
Publication date: 2023-12-15
Anticipated expiration: 2043-11-13
Also published as: CN117236526B

Abstract

The invention provides a method, a device, electronic equipment and a storage medium for determining methane annual emission of rivers, relates to the technical field of methane accounting, and aims to solve the technical problem that methane emission conditions of rivers in various local areas cannot be known. The determining method comprises the following steps: acquiring a database; determining at least one modeling river, and extracting river information of the modeling river; constructing an accounting model based on river information of the modeling river; dividing the region to be counted into at least two sub-regions; determining at least one river to be analyzed from the subareas; calculating the dissolved concentration and bubbling flux of methane based on river information and an accounting model of a river to be analyzed; calculating the discharge flux of river methane to be analyzed; determining methane annual emission of all rivers in the subarea; and finally, calculating the annual emission of river methane in the area to be counted. The method provided by the invention can obtain the total emission amount of river methane and the spatial distribution of methane emission, and can intuitively know the emission conditions of river methane in the whole and partial areas.

Description

River methane annual emission determining method, river methane annual emission determining device, electronic equipment and storage medium

Technical Field

The invention relates to the technical field of methane accounting, in particular to a river methane annual emission determining method, a river methane annual emission determining device, electronic equipment and a storage medium.

Background

Methane is a powerful greenhouse gas whose global warming potential (GWP, global warming potential) is 27 times that of carbon dioxide (time scale 100 years), and is one of the important factors leading to global climate change. The river receives a large amount of organic carbon from the land ecosystem and is transformed in the river by methanogenic bacteria to produce methane which is discharged to the atmosphere. It is reported in literature that rivers produce large amounts of methane each year, a hot zone of global methane emissions, contributing significantly to global methane balance accounting.

The accurate nuclear calculation of river methane emission is the basis of greenhouse gas emission reduction work, and has very important significance for maintaining carbon balance and relieving global climate change. The current method for accounting river methane emission in the literature mainly comprises the steps of collecting published paper methane emission data, and obtaining the total emission of river methane by using Monte Carlo data simulation upscaling method. The Monte Carlo simulation based on collecting methane emission data accounts for the following disadvantages in terms of river methane emissions:

(1) This method relies heavily on the reliability and spatial distribution of the raw data of the collected literature.

(2) The main influencing factors and mechanisms of the methane emission of the river cannot be revealed, and the space-time difference of the methane emission cannot be reasonably explained.

(3) The calculated methane emission is a total amount, the spatial distribution of river methane emission cannot be obtained, and the emission reduction work of regional pertinence cannot be carried out according to the emission.

Therefore, the Monte Carlo simulation method based on collecting methane emission data has great uncertainty in accounting for the emission of large-scale river methane, which is disadvantageous to the establishment of greenhouse gas emission inventory and the development of greenhouse gas emission reduction work.

Therefore, the invention provides a method for determining the annual methane emission of a river, which is necessary to solve the technical problem that the methane emission condition of the river in each local area cannot be known.

Disclosure of Invention

The invention aims to solve or improve the technical problem that the methane emission condition of rivers in various local areas cannot be known in the related technology.

The invention provides a method for determining annual methane emission of a river.

The second aspect of the invention provides a river methane annual emission determining device.

A third aspect of the invention provides an electronic device.

A fourth aspect of the present invention provides a storage medium.

The invention provides a method for determining annual emission of methane in a river, which comprises the following steps: acquiring a database, wherein the database comprises river information; determining at least one modeling river for modeling, and extracting river information of the modeling river from the river information; constructing an accounting model based on river information of the modeling river; dividing the region to be counted into at least two sub-regions; determining at least one river to be analyzed from the subareas; extracting river information of a river to be analyzed from the river information; transmitting river information of the river to be analyzed to an accounting model, and calculating the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model; calculating the methane emission flux of the river to be analyzed based on the dissolved concentration and bubbling flux of the methane of the river to be analyzed; determining the annual emission of methane in the river in each sub-area based on the emission flux of methane in all the rivers to be analyzed in the corresponding sub-area; and determining the annual emission of river methane in the area to be counted based on the annual emission of river methane in all the subareas.

According to the river methane annual emission determining method provided by the invention, a database is firstly obtained, and the database comprises river information; the database is constructed in advance, each river can be sampled through sampling points to know the river information of each river, data of related documents published in the global scope can be collected to know the river information of each river, modeling is started after the river information is obtained, specifically, at least one modeling river for modeling is determined, then the river information of the modeling river is extracted from the river information, and an accounting model is constructed based on the river information of the modeling river; the method mainly comprises the steps of extracting representative rivers, and then constructing the association relation between the dissolved concentration and bubbling flux of methane in the rivers and the parameters such as the water temperature of the rivers, the concentration of soluble organic carbon in the water body of the rivers, the flow rate of the rivers, the water depth of the rivers, the concentration of ammonium in the water body of the rivers, the total phosphorus concentration in the water body of the rivers and the like, so that the dissolved concentration and the bubbling flux of the methane can be calculated according to the corresponding relation in the accounting model for other rivers in the later stage. After modeling, the calculation of the dissolved concentration and bubble flux for each river was started. Specifically, the area to be counted is divided into at least two sub-areas, for example, the area to be counted is divided into grid shapes according to longitude and latitude, each grid area is determined to be one sub-area, then at least one river to be analyzed is determined from the sub-areas, and the database comprises river information of each river collected in advance and river information of the river to be analyzed, so that the dissolving concentration and bubbling flux of methane of the river to be analyzed can be calculated based on the river information and the accounting model of the river to be analyzed; the method comprises the steps of calculating the methane emission flux of a river to be analyzed based on the methane concentration and the bubbling flux of the river to be analyzed, wherein the methane concentration is the amount of methane dissolved in a unit volume of water body, and the bubbling flux is the amount of methane emitted to the atmosphere in a unit time unit area in a bubbling mode; the emission flux is equal to bubbling flux and diffusion flux, the diffusion flux is obtained through concentration calculation, the annual emission of river methane in each subarea is determined based on the emission flux of methane in all rivers to be analyzed in each subarea, and finally the annual emission of river methane in the area to be counted is determined based on the annual emission of river methane in all subareas. The invention establishes a set of method for calculating the global river methane emission, can predict the dissolved concentration and the emission flux of the river methane in any area, further estimate the total emission amount of the river methane, and accumulate the total emission amount of the methane in all subareas to obtain the global river methane total emission amount. In addition, when the method is used for calculating the methane emission of the river, the methane emission of each river segment, river basin and global scale can be calculated, the spatial distribution of the methane emission of the river can be obtained, namely, the methane emission conditions of the river in the whole and local areas can be intuitively known, and the traditional Monte Carlo simulation method can only calculate the total methane emission amount and cannot intuitively know the methane emission conditions of the river in each local area.

In some technical solutions, optionally, the database further includes a frozen soil area and a non-frozen soil area, the river information of the river to be analyzed is transmitted to an accounting model, and the step of calculating the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model specifically includes: determining whether a river to be analyzed is in a frozen soil area or a non-frozen soil area; under the condition that a river to be analyzed is in a frozen soil area, river information of the river to be analyzed is conveyed to a first accounting model, and the dissolved concentration and bubbling flux of methane of the river to be analyzed are calculated through the first accounting model; under the condition that the river to be analyzed is in the frozen soil area, river information of the river to be analyzed comprises water temperature of the river to be analyzed, concentration of soluble organic carbon in water body of the river to be analyzed, flow rate of the river to be analyzed and water depth of the river to be analyzed; under the condition that the river to be analyzed is in a non-frozen soil area, river information of the river to be analyzed is conveyed to a second accounting model, and the dissolved concentration and bubbling flux of methane of the river to be analyzed are calculated through the second accounting model; under the condition that the river to be analyzed is in a non-frozen soil area, the river information of the river to be analyzed comprises the water temperature of the river to be analyzed, the ammonium concentration in the water body of the river to be analyzed, the total phosphorus concentration in the water body of the river to be analyzed, the flow velocity of the river to be analyzed and the water depth of the river to be analyzed.

In the technical scheme, the database also comprises a frozen soil area and a non-frozen soil area, wherein the non-frozen soil area and the frozen soil area can be divided based on the distributed data products of the frozen soil and the non-frozen soil in the global space, or the whole area to be counted can be analyzed to belong to the frozen soil area and the whole area to be counted can be analyzed to belong to the non-frozen soil area by referring to data, basic data and the like, and the dominant factors of the dissolution concentration of methane and bubbling flux are different in the non-frozen soil area and the frozen soil area, so that the river to be analyzed is also required to be determined in the frozen soil area or the non-frozen soil area; in the case of a river under analysis in a frozen soil region, river information of the river under analysis includes water temperature of the river under analysis, concentration of soluble organic carbon in a water body of the river under analysis, flow rate of the river under analysis, and water depth of the river under analysis.

The water temperature of the river to be analyzed, the concentration of the soluble organic carbon in the water body of the river to be analyzed, the flow rate of the river to be analyzed, the water depth of the river to be analyzed, the dissolved concentration of methane in the river to be analyzed and the bubbling flux of the methane in the river to be analyzed satisfy the following relational expression:

；

wherein,Logas a logarithmic function, CFor the dissolved concentration of methane in the river to be analyzed,Tfor the water temperature of the river to be analysed [DOC]Is the concentration of the soluble organic carbon in the water body,vis the flow rate of the river to be analyzed,F _E for the bubbling flux of methane in the river to be analyzed,his the water depth of the river to be analyzed;

namely, under the condition that the river to be analyzed is in the frozen soil area, when the water temperature of the river to be analyzed, the concentration of the soluble organic carbon in the water body of the river to be analyzed, the flow rate of the river to be analyzed, the water depth of the river to be analyzed and other information are extracted from the database, the dissolved concentration of methane in the river to be analyzed and the bubbling flux of methane in the river to be analyzed can be calculated according to the formula.

Under the condition that the river to be analyzed is in a non-frozen soil area, river information of the river to be analyzed comprises water temperature of the river to be analyzed, ammonium concentration in a water body of the river to be analyzed, total phosphorus concentration in the water body of the river to be analyzed, flow velocity of the river to be analyzed and water depth of the river to be analyzed;

the water temperature of the river to be analyzed, the ammonium concentration in the water body of the river to be analyzed, the total phosphorus concentration in the water body of the river to be analyzed, the flow rate of the river to be analyzed, the water depth of the river to be analyzed, the dissolved concentration of methane in the river to be analyzed and the bubbling flux of methane in the river to be analyzed satisfy the following relational expression:

；

Wherein,Logas a logarithmic function,Cfor the dissolved concentration of methane in the river to be analyzed,Tfor the water temperature of the river to be analysed [NH ₄ ]Is the concentration of ammonium radical in water bodyTP]Is the concentration of total phosphorus in the water body,vis the flow rate of the river to be analyzed,F _E for the bubbling flux of methane in the river to be analyzed,his the water depth of the river to be analyzed.

Namely, under the condition that the river to be analyzed is in a non-frozen soil area, when the water temperature of the river to be analyzed, the ammonium concentration in the water body of the river to be analyzed, the total phosphorus concentration in the water body of the river to be analyzed, the flow rate of the river to be analyzed, the water depth of the river to be analyzed and other information are extracted from the database, the dissolved concentration of methane in the river to be analyzed and the bubbling flux of methane in the river to be analyzed can be calculated according to the formula.

Wherein the accounting model comprises a first accounting model and a second accounting model.

According to the invention, the dissolved concentration and bubbling flux of methane can be calculated according to the region (frozen soil region or non-frozen soil region) where the river to be analyzed is located and then according to different calculation formulas, so that the measurement accuracy is improved.

In some embodiments, optionally, in the step of calculating the methane emission flux of the river to be analyzed based on the dissolved concentration of methane and the bubbling flux of the river to be analyzed, the dissolved concentration of methane, the bubbling flux of methane, and the methane emission flux satisfy the following relationship:

；

Wherein,F _D for the diffusion flux of methane,kas a function of the gas transmission coefficient,Cis the dissolved concentration of methane,C _eq is the theoretical equilibrium concentration of methane when the water body is in equilibrium with the atmosphere,Ffor the discharge flux of methane,F _E is the bubbling flux of methane,vfor the flow rate of the river to be analyzed,hfor the depth of the river to be analyzed,u ₁₀ for wind speeds at a height of 10m from the water surface,Sc _CH4 is the schmitt number of methane,nfor Schmidt coefficient, when wind speed is less than 3.6 m.s ^-1 Time of dayn2/3, when the wind speed is more than or equal to 3.6ms ^-1 Time of daynIs 1/2 of the total number of the components,βfor the bensen coefficient corrected for water temperature and salinity,V _m is the gas molar volume of the local atmosphere,Pfor the partial pressure of methane in the local atmosphere,Tto be separated intoThe water temperature of the river is analyzed,lnin order to be a natural logarithm,S‰is the salinity of the water body.

In some embodiments, optionally, in the step of determining the annual emission of methane from the river in each sub-zone based on the emission fluxes of methane from all the rivers to be analyzed in the corresponding sub-zone, the emission fluxes of methane and the annual emission of methane from the sub-zone satisfy the following formula:

；

wherein,Emissionis the annual emission of methane from rivers in a sub-area,the symbols are accumulated and the symbols are accumulated,ifor the number of river sequences to be analyzed,nfor the total number of rivers to be analyzed within the sub-zone, F _i For rivers to be analysediIs used for the methane emission flux of the (a),A _i for rivers to be analysediIs a water surface area of the water tank.

In some embodiments, optionally, after the step of determining the annual emission of methane in the river in the corresponding sub-area, the method further comprises: comparing the annual emission of river methane in the subarea with a preset emission threshold; and when the annual emission amount of river methane in the subarea is larger than a preset emission amount threshold value, determining the corresponding subarea as an emission exceeding area.

In the technical scheme, after the step of determining the annual emission of the river methane in the corresponding subarea, the annual emission of the river methane in the subarea is compared with a preset emission threshold, and if the annual emission of the river methane in the subarea is greater than the preset emission threshold, the corresponding subarea is indicated as an emission exceeding area. The invention can judge whether the annual emission of the river methane in each subarea exceeds the standard, can intuitively know the emission condition of the river methane in the local area, and solves the technical problem that the traditional Monte Carlo simulation method can only calculate the total emission of the methane and can not know the local emission.

In some embodiments, optionally, the method for determining annual emission of methane in a river further comprises: calculating population density of each emission exceeding area; comparing the population density of the emission exceeding area with a preset population density; and when the population density of the emission exceeding area is greater than or equal to the preset population density, determining the emission exceeding area as an emission reduction potential area.

In the technical scheme, after the corresponding subarea is determined to be the emission exceeding area, the population density of each emission exceeding area is calculated, if the population density of the emission exceeding area is greater than or equal to the preset population density, the area is proved to experience serious interference of artificial activities and has emission reduction potential, and the emission exceeding area is determined to be the emission reduction potential area. The invention can calculate and obtain the emission reduction potential of methane, namely, the population density of the emission exceeding area is larger, which indicates that the emission reduction potential is larger, so that targeted emission reduction measures can be formulated according to the emission conditions of different areas, and the invention has important value for greenhouse gas emission reduction work.

The second aspect of the present invention provides a river methane annual emission determining apparatus comprising: the acquisition module is used for acquiring a database, wherein the database comprises river information; a determination module for determining at least one modeled river for modeling; the acquisition module is also used for extracting river information of the modeling river from the river information; the determining module is also used for constructing an accounting model based on river information of the modeling river; the region dividing module is used for dividing the region to be counted into at least two sub-regions; a determining module for determining at least one river to be analyzed from the sub-region; the acquisition module is also used for extracting river information of the river to be analyzed from the river information; the data transmission module is used for transmitting river information of the river to be analyzed to the accounting model so as to calculate the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model; the determining module is also used for calculating the methane emission flux of the river to be analyzed based on the dissolved concentration and bubbling flux of the methane of the river to be analyzed; the determining module is further used for determining the annual emission of methane in the river in each subarea based on the emission flux of methane in all the rivers to be analyzed in the corresponding subarea; the determining module is further used for determining the annual emission of river methane in the area to be counted based on the annual emission of river methane in all the subareas.

A third aspect of the present application provides an electronic device, comprising a memory and a processor, the memory having stored thereon a computer program or instructions which when executed by the processor implement a method for determining annual emission of methane in a river as provided in any of the first aspects of the present application.

A fourth aspect of the present application provides a storage medium having a program or instructions stored thereon, which when executed by a processor implements a method for determining annual emission of methane in a river as provided in any one of the first aspects of the present application.

Drawings

FIG. 1 shows one of the flow charts of the river methane annual emission determination method provided by the embodiment of the application;

FIG. 2 shows a second flow chart of a river methane annual emission determination method provided by an embodiment of the present application;

FIG. 3 shows a third flow chart of a river methane annual emission determination method provided by an embodiment of the present application;

FIG. 4 shows a block diagram of a river methane annual emission determining apparatus provided by an embodiment of the present application;

FIG. 5 shows a block diagram of an electronic device provided by an embodiment of the application;

fig. 6 is a schematic diagram of a hardware structure of an electronic device implementing an embodiment of the present application.

Detailed Description

As shown in fig. 1, the embodiment provides a method for determining annual emission of methane in a river, which includes the following steps:

s102: acquiring a database, wherein the database comprises river information;

s104: determining at least one modeling river for modeling, and extracting river information of the modeling river from the river information;

s106: constructing an accounting model based on river information of the modeling river;

s108: dividing the area to be counted into at least two sub-areas, and determining at least one river to be analyzed from the sub-areas;

s110: extracting river information of a river to be analyzed from the river information;

s112: transmitting river information of the river to be analyzed to an accounting model so as to calculate the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model;

s114: calculating the methane emission flux of the river to be analyzed based on the dissolved concentration and bubbling flux of the methane of the river to be analyzed;

s116: determining the annual emission of methane in the river in each sub-area based on the emission flux of methane in all the rivers to be analyzed in the corresponding sub-area;

s118: and determining the annual emission of river methane in the area to be counted based on the annual emission of river methane in all the subareas.

According to the river methane annual emission determining method provided by the embodiment, a database is firstly obtained, and the database comprises river information; the database is constructed in advance, each river can be sampled through sampling points to know the river information of each river, data of related documents published in the global scope can be collected to know the river information of each river, modeling is started after the river information is obtained, specifically, at least one modeling river for modeling is determined, then the river information of the modeling river is extracted from the river information, and an accounting model is constructed based on the river information of the modeling river; the method mainly comprises the steps of extracting representative rivers, and then constructing the association relation between the dissolved concentration and bubbling flux of methane in the rivers and the parameters such as the water temperature of the rivers, the concentration of soluble organic carbon in the water body of the rivers, the flow rate of the rivers, the water depth of the rivers, the concentration of ammonium in the water body of the rivers, the total phosphorus concentration in the water body of the rivers and the like, so that the dissolved concentration and the bubbling flux of the methane can be calculated according to the corresponding relation in the accounting model for other rivers in the later stage. After modeling, the calculation of the dissolved concentration and bubble flux for each river was started. Specifically, dividing the area to be counted into at least two sub-areas, for example, dividing the world into grid shapes according to longitude and latitude, determining each grid area as one sub-area, and then determining at least one river to be analyzed from the sub-areas, wherein the database comprises river information of each river acquired in advance and river information of the river to be analyzed, so that the river information of the river to be analyzed can be conveyed to an accounting model to calculate the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model; the method comprises the steps of calculating the methane emission flux of a river to be analyzed based on the methane concentration and the bubbling flux of the river to be analyzed, wherein the methane concentration is the amount of methane dissolved in a unit volume of water body, and the bubbling flux is the amount of methane emitted to the atmosphere in a unit time unit volume in a bubbling mode; the emission flux is equal to bubbling flux and diffusion flux, the annual emission amount of the river methane in each subarea is determined based on the emission flux of the methane of all the rivers to be analyzed in each subarea, and the annual emission amount of the river methane in the area to be counted is determined based on the annual emission amount of the river methane in all the subareas. The invention establishes a set of method for calculating the global river methane emission, can predict the dissolved concentration and the emission flux of the river methane in any area, further estimate the total emission amount of the river methane, and accumulate the total emission amount of the methane in all subareas to obtain the global river methane total emission amount. In addition, when the method is used for calculating the methane emission of the river, the methane emission of each river segment, river basin and global scale can be calculated, the spatial distribution of the methane emission of the river can be obtained, namely, the methane emission conditions of the river in the whole and local areas can be intuitively known, and the traditional Montechol simulation method can only calculate the total methane emission, so that the methane emission conditions of the river in each local area can not be intuitively known.

In some embodiments, optionally, the database further includes a frozen soil area and a non-frozen soil area, and the step of delivering river information of the river to be analyzed to the accounting model to calculate the dissolved concentration and bubbling flux of methane of the river to be analyzed by the accounting model specifically includes: determining whether a river to be analyzed is in a frozen soil area or a non-frozen soil area; under the condition that a river to be analyzed is in a frozen soil area, river information of the river to be analyzed is conveyed to a first accounting model, and the dissolved concentration and bubbling flux of methane of the river to be analyzed are calculated through the first accounting model; under the condition that the river to be analyzed is in the frozen soil area, river information of the river to be analyzed comprises water temperature of the river to be analyzed, concentration of soluble organic carbon in water body of the river to be analyzed, flow rate of the river to be analyzed and water depth of the river to be analyzed; under the condition that the river to be analyzed is in a non-frozen soil area, river information of the river to be analyzed is conveyed to a second accounting model, and the dissolved concentration and bubbling flux of methane of the river to be analyzed are calculated through the second accounting model; under the condition that the river to be analyzed is in a non-frozen soil area, the river information of the river to be analyzed comprises the water temperature of the river to be analyzed, the ammonium concentration in the water body of the river to be analyzed, the total phosphorus concentration in the water body of the river to be analyzed, the flow velocity of the river to be analyzed and the water depth of the river to be analyzed.

In this embodiment, the database further includes a frozen soil area and a non-frozen soil area, wherein the non-frozen soil area and the frozen soil area can be divided based on the distributed data of the spatial distribution of the frozen soil and the non-frozen soil in the whole world, or the areas of the whole area to be counted belong to the frozen soil area and the areas of the non-frozen soil area can be analyzed by referring to data and basic data, and the like, and the factors leading in the dissolution concentration and bubbling flux of methane are different in the non-frozen soil area and the frozen soil area, so the invention also needs to determine whether the river to be analyzed is in the frozen soil area or the non-frozen soil area; under the condition that the river to be analyzed is in the frozen soil area, river information of the river to be analyzed comprises water temperature of the river to be analyzed, concentration of soluble organic carbon in water body of the river to be analyzed, flow rate of the river to be analyzed and water depth of the river to be analyzed;

；

Wherein,F _D for the diffusion flux of methane,kas a function of the gas transmission coefficient,Cis the dissolved concentration of methane,C _eq is the theoretical equilibrium concentration of methane when the water body is in equilibrium with the atmosphere,Ffor the discharge flux of methane,F _E is the bubbling flux of methane,vfor the flow rate of the river to be analyzed,hfor the depth of the river to be analyzed,u ₁₀ for wind speeds at a height of 10m from the water surface,Sc _CH4 is the schmitt number of methane,nfor Schmidt coefficient, when wind speed is less than 3.6 m.s ^-1 Time of dayn2/3, when the wind speed is more than or equal to 3.6ms ^-1 Time of daynIs 1/2 of the total number of the components,βfor the bensen coefficient corrected for water temperature and salinity,V _m is the gas molar volume of the local atmosphere,Pfor the partial pressure of methane in the local atmosphere,Tfor the water temperature of the river to be analyzed,lnin order to be a natural logarithm,S‰is the salinity of the water body.

In some embodiments, optionally, in the step of determining the annual emission of methane from the river in each sub-zone based on the emission flux of methane from all the rivers to be analyzed in the corresponding sub-zone, the emission flux of methane and the annual emission of methane from the sub-zone satisfy the following formula:

；

In some embodiments, optionally, after the step of determining the annual emission of river methane in the corresponding sub-area, further comprising: comparing the annual emission of river methane in the subarea with a preset emission threshold; and when the annual emission amount of river methane in the subarea is larger than a preset emission amount threshold value, determining the corresponding subarea as an emission exceeding area.

In this embodiment, after the step of determining the annual emission of river methane in the corresponding sub-zone, the annual emission of river methane in the sub-zone is compared with a preset emission threshold, and if the annual emission of river methane in the sub-zone is greater than the preset emission threshold, the corresponding sub-zone is indicated as an out-of-standard emission zone. The invention can judge whether the annual emission of the river methane in each subarea exceeds the standard, can intuitively know the emission condition of the river methane in the local area, and solves the technical problem that the traditional Monte Carlo simulation method can only calculate the total emission of the methane and can not know the local emission.

In some embodiments, optionally, the methane year emission determination method further comprises: calculating population density of each emission exceeding area; comparing the population density of the emission exceeding area with a preset population density; and when the population density of the emission exceeding area is greater than or equal to the preset population density, determining the emission exceeding area as an emission reduction potential area.

In this embodiment, after determining that the corresponding sub-region is an emission exceeding region, calculating population density of each emission exceeding region, if population density of the emission exceeding region is greater than or equal to a preset population density, indicating that the region is subject to serious interference of human activities, and has emission reduction potential, determining the emission exceeding region as an emission reduction potential region. The invention can calculate and obtain the emission reduction potential of methane, namely, the population density of the emission exceeding area is larger, which indicates that the emission reduction potential is larger, so that targeted emission reduction measures can be formulated according to the emission conditions of different areas, and the invention has important value for greenhouse gas emission reduction work.

As shown in fig. 2, another embodiment of the present invention provides a methane annual emission determining method, including the steps of:

s202: and establishing a database.

Specifically, the method is characterized in that high-density in-situ sampling monitoring is carried out in the whole country, the sampling year is 2017-2022, the water system involved is provided with 90 sampling points such as Zhujiang, changjiang, huai river, huang river, sea river and Liaohe, etc., seasonal sampling is carried out, the sampling time is spring (4 months), summer (8 months) and autumn and winter (10 months to 12 months), 258 monitoring data are provided in total, in-situ sampling monitoring of a river in a certain city is carried out, the sampling year is 2018-2020, the river is a river around a certain city, 19 sampling points are provided, and the sampling is carried out for the week (3 months to 11 months), and 196 monitoring data are provided in total; data of relevant literature published worldwide was collected, published year 1990 to 2023, 295 was involved, total 16769 pieces of monitored data, and a global river methane emission database was constructed by integrating all the above data.

S204: key influencing factors are identified.

Specifically, the method comprises the steps of dividing a non-frozen soil area and a frozen soil area by published data products of the spatial distribution of the frozen soil and the non-frozen soil in the world, identifying that the factors influencing the dissolved concentration and bubbling flux of methane in the non-frozen soil area and the frozen soil area are different by utilizing various data analysis means such as correlation analysis, regression analysis, random forest analysis and the like, wherein the factors influencing the dissolved concentration and bubbling flux of methane in the non-frozen soil area are nutrient salt levels including ammonia nitrogen (NH ₄ ⁺ ) And Total Phosphorus (TP), the dominant factor in the frozen soil zone affecting methane sequestration concentration and bubble flux is the dissolved organic carbon concentration (DOC); wherein NH is ₄ ⁺ TP and DOC are all conventional water chemistry parameters that are easy to monitor.

S206: and respectively establishing prediction models of methane dissolved concentration and bubbling flux in the frozen soil and non-frozen soil areas.

Specifically, a linear regression model of the methane dissolved concentration and bubbling flux of frozen soil and non-frozen soil areas and key influence factors thereof is established by utilizing database data, and the specific results are as follows:

non-frozen soil region:

；

frozen soil area:

；

wherein,Logas a logarithmic function,CandF _E represents the concentration of methane dissolved in the water (. Mu. Mol.L) ^-1 ) And bubble flux (mmol.m) ^-2 ·d ^-1 ），TIs the water temperature (DEG C) of the river to be analyzed [NH ₄ ]Is the concentration of ammonium radical (mg.L) in water body ^-1 ），[TP]Is the concentration (mg.L) of total phosphorus in the water body ^-1 ），vIs the flow rate (m.s) of the river to be analyzed ^-1 ），hIs the water depth (m) of the river to be analyzed [DOC]Is the concentration (mg.L) of the soluble organic carbon in the water body ^-1 ）。

S208: the diffusion flux and the total flux of river methane were calculated.

Specifically, the diffusion flux of river methane is calculated according to a water-gas interface diffusion model method, and the specific formula is as follows:

；

wherein,F _D for the diffusion flux of methane,Fis the emission of methaneFlux, i.e. total flux (mmol.m) ^-2 ·d ^-1 ），kIs the gas transmission coefficient (m.d ^-1 ），C _eq Is the theoretical equilibrium concentration (mu mol.L) of methane when water body and atmosphere are balanced ^-1 ）。

Wherein,kcalculation is performed according to hydrologic and meteorological parameters:

；

wherein,vfor the average flow rate (m.s) of the river to be analyzed ^-1 ），hFor the water depth (m) of the river to be analyzed,u ₁₀ is the wind speed (m.s) at the position 10m from the water surface ^-1 ），Sc _CH4 Is the schmitt number of methane,nfor Schmidt coefficient, when wind speed is less than 3.6 m.s ^-1 When the wind speed is greater than or equal to 3.6ms, n is 2/3 ^-1 And n is 1/2.

Wherein,Sc _CH4 using water temperatureTC) is calculated by the following formula:

；

wherein,Tfor the river water temperature (c) to be analyzed,C _eq calculation was performed according to henry's law:

；

Wherein,βto the Benson coefficient (L.L) corrected for water temperature and salinity ^-1 ·atm ^-1 ），V _m Is the gas molar volume (L.mol) ^-1 ），PIs the partial pressure (atm) of methane in the local atmosphere.

Wherein the method comprises the steps ofβThe calculation is performed by the following formula:

；

wherein,Tfor the river water temperature (K) to be analyzed,lnis natural logarithm, S per mill is salinity of water body.

S210: accounting for annual emissions of river methane.

Specifically, the emission flux of methane is calculated by using the nutrient salt or DOC data of each grid or river reach, the emission amount of the current grid or river reach methane is obtained according to the emission flux and the water surface area, and the annual emission amount of the river methane is accumulated and calculated finally, the specific formula is as follows:

；

wherein,Emissionannual emission of methane (Tg. Yr) for river ^-1 ），The symbols are accumulated and the symbols are accumulated,iis the number of river segments of a river,F _i for the discharge flux (mmol.m) at the grid or river reach ^-2 ·d ^-1 ），A _i Is the water surface area (m) ² ) 16 is the molar mass of methane, 365 is the number of days per year, 10 ^-15 Is a conversion coefficient.

S212: the emission hot zone of river methane is identified and emission reduction potential is determined.

Specifically, according to the spatial distribution of river methane emissions, a river segment with methane emissions higher than the global average value is determined as the emission hot zone of methane. The emission of river methane is divided and accumulated spatially according to land utilization type or population density, and river or high population density area (more than 20 people/km) ² ) The river of the river means that the river is subjected to serious artificial activity, has emission reduction potential, and can improve water quality by reducing artificial management strategies such as nutrient salt input and the like, thereby reducing methane emission of the river. Theoretically, when the concentration of the nutrient salt is reduced to the level of the natural state, the reduction amount of the river methane emission is the maximum emission reduction potential of the river methane emission.

As shown in fig. 3, another embodiment of the present invention provides a methane annual emission determining method including the steps of:

s302: collecting data and establishing a global river methane database;

s304: identifying control factors influencing the methane concentration flux of the river, and establishing a prediction model of the methane concentration flux;

s306: calculating the emission flux of methane on the river reach according to the prediction model;

s308: calculating the discharge amount of methane on the river reach by utilizing the discharge flux of methane and the corresponding water surface area, and accumulating to obtain the total discharge amount of the river;

s310: and identifying hot spot areas of methane emission according to the spatial distribution of the methane emission of the river, and determining the emission reduction potential of the river methane in the artificially disturbed area.

The invention can realize the following beneficial effects:

(1) The invention establishes a set of optimized method for calculating the global river methane emission, predicts the concentration and flux of the river methane by using a prediction model of the concentration and flux of the methane on a global scale, further estimates the emission of the river methane and accumulates the emission to obtain the global river methane total emission.

(2) The parameters in the prediction model of the methane concentration and flux are conventional water chemistry parameters which are easy to measure, can realize high-frequency monitoring, are not limited by data acquisition, and are easy to popularize and apply in a large scale.

(3) When the method provided by the invention is used for calculating the methane emission of the river, the resolution is higher, the methane emission of the river reach, the river basin and the global scale can be calculated, the spatial distribution of the methane emission of the river can be obtained, and the traditional Monte Carllo upscaling method can only calculate the total methane emission.

(4) By adopting the river methane emission accounting method, the river methane emission under the global scale is visualized, the whole and partial area river methane emission conditions can be intuitively known, the methane emission hot zone is identified, the methane emission reduction potential can be calculated, and the targeted emission reduction measures are formulated according to the emission conditions of different areas, so that the river methane emission accounting method has important value for greenhouse gas emission reduction work.

(5) Although the method and the model provided by the invention are aimed at the calculation of the methane emission of the river, the method and the model are not limited to the river, and the method and the model can be applied to the calculation of the methane emission of other inland water bodies, such as lakes, ponds, reservoirs and the like.

As shown in fig. 4, a second aspect of the present invention provides a river methane annual emission determining apparatus 1 including an acquisition module 12, a region dividing module 14, a determining module 16, and a data delivery module 18, the acquisition module 12 being configured to acquire a database including river information; the determination module 16 is used to determine at least one modeled river for modeling; the acquisition module 12 is further configured to extract river information modeling a river from the river information; the determination module 16 is also configured to construct an accounting model based on river information modeling the river; the region dividing module 14 is configured to divide the region to be counted into at least two sub-regions; the determining module 16 is configured to determine at least one river to be analyzed from the sub-area; the acquisition module 12 is further used for extracting river information of the river to be analyzed from the river information; the data transmission module 18 is used for transmitting river information of the river to be analyzed to the accounting model so as to calculate the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model; the determination module 16 is further configured to calculate a methane emission flux of the river to be analyzed based on the dissolved concentration and the bubbling flux of methane of the river to be analyzed; the determining module 16 is further configured to determine annual emission of methane from the river in each sub-area based on the emission flux of methane from all the rivers to be analyzed in the corresponding sub-area; the determining module 16 is further configured to determine annual emissions of river methane for the area to be counted based on the annual emissions of river methane for all the sub-areas.

The river methane annual emission determining device in the embodiment of the application can be a device, and can also be a component, an integrated circuit or a chip in a terminal. The device may be a mobile electronic device or a non-mobile electronic device. The mobile electronic device may be a mobile phone, a tablet computer, a notebook computer, a palm computer, a vehicle-mounted electronic device, a wearable device, an ultra mobile personal computer, a netbook or a personal digital assistant, and the non-mobile electronic device may be a server, a network attached memory, a personal computer, a television, a teller machine or a self-service machine, which is not particularly limited in the embodiments of the present application.

The methane year emission determination device in the embodiment of the application may be a device having an operating system. The operating system may be an Android operating system, an iOS operating system, or other possible operating systems, and the embodiment of the present application is not limited specifically.

The river methane annual emission determining device provided by the embodiment of the application can realize each process realized by the method embodiment, and in order to avoid repetition, the description is omitted.

As shown in fig. 5, a third aspect of the present application provides an electronic device 700, including a processor 701, a memory 702, and a program or an instruction stored in the memory 702 and capable of running on the processor 701, where the program or the instruction is executed by the processor 701 to implement the above-mentioned methane annual emission determining method and achieve the same technical effects, and is not repeated herein.

It should be noted that, the electronic device in the embodiment of the present application includes a mobile electronic device and a non-mobile electronic device.

It should be noted that, in the embodiment of the present application, the electronic device includes a mobile electronic device, such as a mobile phone, and may also include a non-mobile electronic device, such as a computer.

Fig. 6 is a schematic diagram of a hardware structure of another electronic device 2000 implementing an embodiment of the present application.

The electronic device 2000 includes, but is not limited to: radio frequency unit 2001, network module 2002, audio output unit 2003, input unit 2004, sensor 2005, display unit 2006, user input unit 2007, interface unit 2008, memory 2009, and processor 2010.

Those skilled in the art will appreciate that the electronic device 2000 may further include a power source 2011 (such as a battery) for powering the various components, where the power source 2011 may be logically connected to the processor 2010 through a power management system to perform functions such as managing charging, discharging, and power consumption. The electronic device structure shown in fig. 4 does not constitute a limitation of the electronic device, and the electronic device may include more or less components than shown, or may combine certain components, or may be arranged in different components, which are not described in detail herein.

Wherein the user input unit 2007 receives a first input.

Processor 2010 generates and stores a corresponding original operation record according to the first input, wherein the original operation record includes at least one original operation node;

the user input unit 2007 receives a second input to a target operation node among the operation nodes.

Processor 2010 generates adjusted simulated operational records in response to the second input.

And controlling the electronic equipment to run corresponding programs or functions according to the simulated operation records.

Optionally, the first input comprises at least one input step, and each original operation node comprises one input step and a corresponding operation result.

Wherein, the operation result is: after receiving the input step, the program or function of the electronic device outputs a feedback result according to the input step.

The input unit 2004 acquires a program or function corresponding to the first input.

The memory 2009 records each input step and the corresponding operation result according to the input sequence of the input steps.

The processor 2010 correspondingly stores the program or function corresponding to the first input, the input steps and the operation result in the input order, and forms an original operation record.

Optionally, the display unit 2006 displays an identification associated with the original operation record.

The user input unit 2007 receives a third input of the identification.

The display unit 2006 displays the original operation nodes in the original operation record in the input order in response to the third input.

Optionally, the processor 2010 adjusts the target input step corresponding to the target operation node according to the second input, to obtain an adjusted analog input step;

the processor 2010 controls the electronic equipment to run a program or a function corresponding to the target input step according to the analog input step so as to obtain an analog operation result corresponding to the analog input step;

the processor 2010 generates corresponding simulation operation nodes according to the simulation input steps and the simulation operation results, and generates simulation operation records according to the simulation operation nodes;

the input sequence corresponding to the analog operation node is the same as the input sequence corresponding to the target operation node.

Optionally, the user input unit 2007 receives a running input.

Processor 2010, in response to the run input, controls the electronic device to run a corresponding program or function according to the simulated operation record.

Optionally, the processor 2010 separately determines simulated operation results for each simulated operation node in each of the plurality of simulated operation records.

The display unit 2006 displays corresponding prompt information when any two simulated operation records exist and the simulated operation results of corresponding simulated operation nodes in any two simulated operation records are different.

According to the embodiment of the application, the first input of the user is saved, the first input is formed into the original operation node according to each operation step, so that the user can trace back the operation node with the error after the operation error occurs and carry out targeted correction, after the correction, the operation is formed into a complete operation record and is carried out according to the saved correct node and the corrected node, the user is prevented from manually operating from beginning, on one hand, the quick correction of misoperation is realized, on the other hand, the possibility of misoperation again is fundamentally avoided without the need of the user, and the interactive experience of the user is improved.

It should be appreciated that in embodiments of the present application, the input unit 2004 may include a graphics processor 5082 and a microphone 5084, the graphics processor 5082 processing image data of still pictures or video obtained by an image capturing device (e.g., a camera) in a video capturing mode or an image capturing mode.

The display unit 2006 may include a display panel 5122, and the display panel 5122 may be configured in the form of a liquid crystal display, an organic light emitting diode, or the like. The user input unit 2007 includes a touch panel 5142 and other input devices 5144. The touch panel 5142 is also referred to as a touch screen. The touch panel 5142 may include two parts of a touch detection device and a touch controller. Other input devices 5144 can include, but are not limited to, a physical keyboard, function keys (e.g., volume control keys, switch keys, etc.), a trackball, a mouse, a joystick, and so forth, which are not described in detail herein. Memory 2009 may be used to store software programs as well as various data including, but not limited to, application programs and an operating system. Processor 2010 may integrate an application processor with a modem processor, wherein the application processor primarily handles operating systems, user interfaces, applications, etc., and the modem processor primarily handles wireless communications. It will be appreciated that the modem processor described above may not be integrated into the processor 2010.

A fourth aspect of the present application provides a storage medium having stored thereon a program or instructions which when executed by a processor implements a method for determining annual emission of methane in a river as provided by any of the embodiments of the first aspect of the present application.

The processor is a processor in the electronic device in the above embodiment. Storage media include computer readable storage media such as computer read only memory, random access memory, magnetic or optical disks, and the like.

The embodiment of the application further provides a chip, the chip comprises a processor and a communication interface, the communication interface is coupled with the processor, the processor is used for running programs or instructions, the processes of the control method embodiment of the electronic equipment can be realized, the same technical effects can be achieved, and the repetition is avoided, and the description is omitted here.

It should be understood that the chips referred to in the embodiments of the present application may also be referred to as system-on-chip chips, chip systems, or system-on-chip chips, etc.

The application establishes a set of methods for accounting global methane emission, identifying emission hot spot areas and determining methane emission reduction potential, and further establishes a prediction model of the river methane emission by collecting data of the river methane emission in the global range, identifying influence factors capable of controlling the river methane emission by using related analysis, random forest and regression analysis, accounting the emission amount of the global river methane, analyzing the spatial distribution of the global methane emission and identifying the hot spot areas of the methane emission, calculating the emission reduction potential of the river methane, and providing theoretical basis and data support for the greenhouse gas emission reduction work targeted by the subsequent development area. The application is used for solving the problems that the current uncertainty of the discharge amount of methane in a river is large and the hot spot area of methane discharge cannot be identified, and the main application is as follows:

(1) The method has the advantages that the in-situ monitoring of the national river and the collection of the literature data published in the global scope are used for establishing a database, the main factors affecting the emission of the river methane in the global scope are identified, the prediction model of the emission of the river methane is established for estimating the emission of the river methane, the problem that the influence factors and the mechanism of the emission of the river methane in the global scale are unclear is solved, and the problem that the large-scale emission of the river methane is calculated by the Monte Carlo simulation method based on the collection of the emission data of the methane to cause large uncertainty is improved.

(2) The global raster data of the influencing factors are input into the prediction model, so that the spatial distribution of global methane emission is obtained, the hot spot area of methane emission is identified, the current greenhouse gas emission reduction work can be served, the emission quantity of river methane after manual management in the future can be predicted, the emission reduction potential of methane is clear, and the problem that the spatial distribution of methane emission is unclear is solved.

(3) The parameters in the methane emission prediction model are conventional water chemistry parameters (water temperature, ammonia nitrogen, total phosphorus, organic carbon concentration and the like) which are easy to measure, and can be obtained through daily observation of large scale and high space-time resolution, so that the methane emission rate with high space-time resolution is calculated, and the problem of large accounting error caused by large space-time variability of methane emission is solved.

Claims

1. A method for determining annual methane emission in a river, comprising:

acquiring a database, wherein the database comprises river information;

determining at least one modeling river for modeling, and extracting river information of the modeling river from the river information;

constructing an accounting model based on river information of the modeled river;

dividing the region to be counted into at least two sub-regions;

determining at least one river to be analyzed from the sub-area;

extracting river information of the river to be analyzed from the river information;

transmitting river information of the river to be analyzed to the accounting model, and calculating the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model;

calculating the methane emission flux of the river to be analyzed based on the dissolved concentration and bubbling flux of the methane of the river to be analyzed;

determining the annual emission of methane in the river in each sub-area based on the emission flux of methane in all the rivers to be analyzed in the corresponding sub-area;

and determining the annual emission of river methane in the area to be counted based on the annual emission of river methane in all the subareas.

2. The method for determining annual methane emission in a river according to claim 1, wherein the database further comprises a frozen soil area and a non-frozen soil area, the step of transmitting the river information of the river to be analyzed to the accounting model, and calculating the dissolved concentration and bubbling flux of methane in the river to be analyzed by the accounting model specifically comprises:

Determining whether the river to be analyzed is in the frozen soil area or the non-frozen soil area;

the river information of the river to be analyzed is transmitted to a first accounting model under the condition of the frozen soil area, and the dissolved concentration and bubbling flux of methane of the river to be analyzed are calculated through the first accounting model; wherein, in the case that the river to be analyzed is in the frozen soil region, the river information of the river to be analyzed includes the water temperature of the river to be analyzed, the concentration of the soluble organic carbon in the water body of the river to be analyzed, the flow rate of the river to be analyzed and the water depth of the river to be analyzed;

the river information of the river to be analyzed is transmitted to a second accounting model under the condition that the river to be analyzed is in the non-frozen soil area, and the dissolved concentration and bubbling flux of methane of the river to be analyzed are calculated through the second accounting model; wherein, under the condition that the river to be analyzed is in the non-frozen soil area, the river information of the river to be analyzed comprises the water temperature of the river to be analyzed, the ammonium concentration in the water body of the river to be analyzed, the total phosphorus concentration in the water body of the river to be analyzed, the flow velocity of the river to be analyzed and the water depth of the river to be analyzed.

3. The method for determining annual emission of river methane according to claim 2, wherein,

under the condition that the river to be analyzed is in the frozen soil area, the water temperature of the river to be analyzed, the concentration of the soluble organic carbon in the water body of the river to be analyzed, the flow rate of the river to be analyzed, the water depth of the river to be analyzed, the dissolved concentration of the methane of the river to be analyzed and the bubbling flux of the methane of the river to be analyzed satisfy the following relation:

；

wherein,Logas a logarithmic function,Cfor the dissolved concentration of methane in the river to be analyzed,Tfor the water temperature of the river to be analyzed [DOC]Is the concentration of the soluble organic carbon in the water body,vis the flow rate of the river to be analyzed,F _E for the bubbling flux of methane of the river to be analyzed,hfor the water depth of the river to be analyzed;

under the condition that the river to be analyzed is in the non-frozen soil area, the water temperature of the river to be analyzed, the ammonium concentration in the water body of the river to be analyzed, the total phosphorus concentration in the water body of the river to be analyzed, the flow rate of the river to be analyzed, the water depth of the river to be analyzed, the dissolved concentration of methane in the river to be analyzed and the bubbling flux of methane in the river to be analyzed satisfy the following relational expression:

；

Wherein,Logas a logarithmic function,Cfor the dissolved concentration of methane in the river to be analyzed,Tfor the water temperature of the river to be analyzed [NH ₄ ]Is the concentration of ammonium radical in water bodyTP]Is the concentration of total phosphorus in the water body,vis the flow rate of the river to be analyzed,F _E for the bubbling flux of methane of the river to be analyzed,hfor the water depth of the river to be analyzed.

4. The method for determining annual methane emission from a river according to claim 1, wherein in the step of calculating the methane emission flux from the river to be analyzed based on the dissolved concentration of methane and the bubbling flux of methane, the dissolved concentration of methane, the bubbling flux of methane, and the methane emission flux satisfy the following relationship:

；

wherein,F _D for the diffusion flux of methane,kas a function of the gas transmission coefficient,Cis the dissolved concentration of methane,C _eq is the theoretical equilibrium concentration of methane when the water body is in equilibrium with the atmosphere,Ffor the discharge flux of methane,F _E is the bubbling flux of methane,vfor the flow rate of the river to be analyzed,hfor the water depth of the river to be analyzed,u ₁₀ for wind speeds at a height of 10m from the water surface,Sc _CH4 is the schmitt number of methane,nfor Schmidt coefficient, when wind speed is less than 3.6 m.s ^-1 Time of dayn2/3, when the wind speed is more than or equal to 3.6ms ^-1 Time of daynIs 1/2 of the total number of the components,βfor the bensen coefficient corrected for water temperature and salinity,V _m is the gas molar volume of the local atmosphere,Pfor the partial pressure of methane in the local atmosphere,Tfor the water temperature of the river to be analyzed,lnis a natural logarithm of the number of the pairs,S‰is the salinity of the water body.

5. The annual emission determining method of river methane in accordance with claim 1, wherein in the step of determining the annual emission of river methane in the corresponding sub-zone based on the emission fluxes of methane in all the rivers to be analyzed in each sub-zone, the emission fluxes of methane and the annual emission of methane in the sub-zone satisfy the following formula:

；

wherein,Emissionis the annual emission of river methane in a sub-area,the symbols are accumulated and the symbols are accumulated,ifor the number of river sequences to be analyzed,nfor the total number of said rivers to be analysed within a sub-zone,F _i for rivers to be analysediIs used for the methane emission flux of the (a),A _i for rivers to be analysediIs a water surface area of the water tank.

6. The method for determining annual emission of river methane in accordance with claim 1, further comprising, after said step of determining annual emission of river methane in the corresponding sub-area:

comparing the annual emission of river methane in the subarea with a preset emission threshold;

And when the annual emission amount of river methane in the subarea is larger than a preset emission amount threshold value, determining the corresponding subarea as an emission exceeding area.

7. The river methane annual emission determination method of claim 6, further comprising:

calculating population density of each emission exceeding area;

comparing the population density of the emission exceeding area with a preset population density;

and when the population density of the emission exceeding area is greater than or equal to the preset population density, determining the emission exceeding area as an emission reduction potential area.

8. A river methane annual emission amount calculation apparatus, comprising:

the acquisition module is used for acquiring a database, wherein the database comprises river information;

a determination module for determining at least one modeled river for modeling;

the acquisition module is also used for extracting river information of the modeling river from the river information;

the determining module is further used for constructing an accounting model based on river information of the modeling river;

the region dividing module is used for dividing the region to be counted into at least two sub-regions;

a determining module for determining at least one river to be analyzed from the sub-area;

The acquisition module is also used for extracting river information of the river to be analyzed from the river information;

the data transmission module is used for transmitting river information of the river to be analyzed to the accounting model so as to calculate the dissolved concentration and bubbling flux of methane of the river to be analyzed through the accounting model;

the determining module is further used for calculating the emission flux of the methane of the river to be analyzed based on the dissolved concentration and bubbling flux of the methane of the river to be analyzed;

the determining module is further used for determining methane annual emission of the river in each subarea based on the emission flux of methane of all the rivers to be analyzed in the corresponding subarea;

the determining module is further used for determining the methane year emission of the river in the area to be counted based on the methane year emission of the river in all the subareas.

9. An electronic device, comprising:

a memory and a processor, the memory having stored thereon a computer program or instructions which when executed implement the river methane annual emission determination method of any one of claims 1 to 7.

10. A storage medium having stored thereon a program or instructions which when executed by a processor implements the river methane annual emission determination method according to any one of claims 1 to 7.