CN-121980128-A - Method for measuring and calculating carbon reserves of water body in subsidence ponding area

CN121980128ACN 121980128 ACN121980128 ACN 121980128ACN-121980128-A

Abstract

The invention relates to the technical field of mining area data measurement and calculation, and discloses a method for measuring and calculating carbon reserves of water bodies in a subsidence ponding area, which comprises the steps of obtaining geological mining condition data of a mining area where the subsidence ponding area is positioned, and calculating a subsidence data set of the subsidence ponding area; the method comprises the steps of acquiring a remote sensing image data sequence of a subsidence water accumulation area, preprocessing the remote sensing image data sequence to generate a water accumulation area data sequence, superposing the water accumulation area data sequence and subsidence data in a unified coordinate system to obtain a reservoir capacity time sequence, carrying out multipoint sampling and calculating water organic carbon density data, and calculating water carbon reserves at set moments.

Inventors

JIN XUELIANG
MA JUN
WEI XIANGPING
ZHENG LIUGEN
ZHU XIAOJUN
BI JIAYING

Assignees

淮北矿业股份有限公司
安徽大学

Dates

Publication Date: 20260505
Application Date: 20251217

Claims (10)

1. A method for measuring and calculating carbon reserves of a water body in a subsidence ponding area is characterized by comprising the following steps of: Step 100, obtaining geological mining condition data of a mining area where the subsidence water area is located, and calculating a subsidence data set of the subsidence water area based on the geological mining condition data; step 200, acquiring a remote sensing image data sequence of a subsidence water accumulation area, preprocessing, extracting a water accumulation boundary data sequence of the subsidence water accumulation area by adopting a water body index method, and generating a water accumulation area data sequence based on the water accumulation boundary data sequence; Step 300, overlapping the ponding area data sequence and the sinking data in a unified coordinate system, and obtaining a reservoir capacity time sequence; step 400, arranging a plurality of sampling points in a subsidence water accumulation area system by adopting a grid method, collecting a plurality of water sample data, and acquiring water body organic carbon density data of the subsidence water accumulation area at a set moment based on the plurality of water sample data; And 500, calculating the water body carbon reserves at the set moment based on the reservoir capacity time sequence and the water body organic carbon density data of the subsidence water accumulation area at the set moment.
2. The method for measuring and calculating the carbon reserves of the water body in the subsidence water accumulation area according to claim 1, further comprising the following steps: Step 600, constructing a carbon dynamic model, inputting the carbon reserves of the water body at the set moment in the step 500 into the carbon dynamic model, and outputting the change rate of the carbon reserves; step 700, generating dynamic water carbon reserves at set moments based on the carbon reserve change rate and the water carbon reserves at set moments.
3. The method for measuring and calculating carbon reserves of water bodies in a subsidence water accumulation area according to claim 2, wherein the geological mining condition data comprise average mining depth H, mining thickness m and coal seam inclination angle The expected running length L, the trend length L, the dip coefficient q, the horizontal movement coefficient b, the main influence radius r and the exploitation influence propagation angle of the working surface Inflection point offset s; the sinking data calculation formula is as follows: ; ; Wherein, the Representing the ground point A dip value at the point of the dip, Representing the maximum subsidence value of the earth surface, D representing the projected extent of the mined area on the coordinate plane, Are all coordinate variables differentiated by the production unit, Representing a double integral variable.
4. A method for measuring and calculating the carbon reserves of a water body in a subsidence water accumulation area according to claim 3, wherein the step 200 comprises the following specific steps: Step 201, performing radiometric calibration, atmospheric correction and geometric fine correction on a remote sensing image data sequence; Step 202, performing recognition processing on the corrected remote sensing image data sequence by adopting an improved normalized difference water body index method to generate an index image data sequence, wherein the water body index method has the following calculation formula: ; Wherein, the Representing an improved normalized difference water index value, Representing the earth's surface reflectivity of the green band of the satellite image, Representing the earth surface reflectivity of a short wave infrared band in a satellite image; step 203, binarizing the exponential image of the set image time phase by setting an empirical threshold, extracting water body pixels, and vectorizing to obtain accurate ponding boundary vector data And calculate the area thereof , wherein, The image phase is represented.
5. The method for measuring and calculating the carbon reserves of a water body in a subsidence water accumulation area according to claim 4, wherein the step 300 comprises the following specific steps: Step 301, each time phase to be acquired Ponding boundary vector data of (2) Carrying out space superposition analysis with the sinking data set under a unified coordinate system; step 302, for each time phase, determining the water surface elevation of the water accumulation area ; Calculating the water storage capacity of each time phase The calculation formula is as follows: ; Wherein, the Indicating time phase Is characterized by that its water-accumulating container is formed from several layers of water-accumulating container, Indicating time phase Is arranged in the water accumulation area of the water tank, Indicating time phase Is used for the water surface elevation of the water tank, Representation points The elevation of the ground surface to be estimated, Representing an area infinitesimal in a Cartesian coordinate system; Step 303, integrating the water reservoir capacity of each time phase Obtaining a water reservoir capacity data sequence 。
6. The method for measuring and calculating the carbon reserves of the water body in the subsidence water accumulation area according to claim 5, wherein the step 400 comprises the following specific steps: step 401, collecting a plurality of water samples at different water depths by adopting a grid point distribution method, and recording the coordinate position of each sampling point; Step 402, processing the water sample and obtaining total carbon and inorganic carbon of the water sample, and calculating the content of organic carbon based on the total carbon and the inorganic carbon, wherein the calculation formula is as follows: ; Wherein, the The total organic carbon content of the water sample is represented, The total carbon content of the water sample is represented, Representing the inorganic carbon content of a water sample; Step 403, arithmetic average is performed on TOC values of all sampling points to obtain Time-of-day average organic carbon density in water 。
7. The method for measuring and calculating the carbon reserves of the water body in the subsidence water accumulation area according to claim 6, wherein the calculation formula of the carbon reserves of the water body at the set time in the step 500 is as follows: ; Wherein, the And the carbon reserves of the water body at the set time are indicated.
8. The method for measuring and calculating the carbon reserves of the water body in the subsidence water accumulation area according to claim 7, wherein the calculation formula of the carbon dynamic model is as follows: ; Wherein, the Represents the organic carbon reserves of the water body in the water accumulation area at the moment t, The carbon reserve change rate at time t is indicated, The total external carbon input flux at time t is indicated, The total carbon output flux at time t is indicated, Represents the primary production carbon fixation flux in the water body at the moment t, And the respiration decomposition carbon emission flux in the water body at the moment t is represented.
9. The method for measuring and calculating the carbon reserves of the water body in the subsidence ponding region according to claim 8, wherein, ; ; ; ; ; Wherein, the Indicating the sedimentation of the atmosphere, Represents the precipitation amount of water to be produced, Represents the concentration of organic carbon in rainwater, Represents the dry sedimentation flux and the dry sedimentation rate, Representing an input of surface runoff, Indicating that the flow is continuously monitored, Represents the concentration of organic carbon measured by the input water sample of surface runoff, Indicating the input of the underground water, Representing the estimated replenishment quantity of the hydrogeologic model, Represents the organic carbon concentration of the underground water, Representing a sediment re-suspension input, For the wind speed of the wind, In order to achieve an average water depth, Is an experience coefficient; ; ; ; ; Wherein, the Indicating the outflow of the earth's surface, Represents the monitoring flow of the water outlet, Indicating the emission of greenhouse gases, Represents the actual measured CO 2 flux in the field, Indicating the in situ measured CH 4 flux, Is a global warming potential equivalent factor of methane, Indicating that the deposit is buried, Indicating the deposition rate of the particles, Indicating the organic carbon content of the deposit, Is the sediment volume weight.
10. The method for measuring and calculating the carbon reserves of the water body in the subsidence ponding region according to claim 9, wherein, ; ; Wherein, the Represents the concentration of chlorophyll a in the water body, Represents the maximum photosynthesis rate of the unit chlorophyll a under the optimal illumination, Represents the light attenuation coefficient of the water body, Representing the light utilization efficiency factor of the water column, Represents the daily average irradiance of the photosynthetic active radiation received by the water surface, Indicating the proportion of photosynthetically active radiation to the total solar radiation, Indicating the duration of the sunlight.

Description

Method for measuring and calculating carbon reserves of water body in subsidence ponding area Technical Field The invention relates to the technical field of mining area data measurement and calculation, in particular to a method for measuring and calculating carbon reserves of a water body in a subsidence ponding area. Background Coal resources are an indispensable component for the development of economic construction in China, but the problem of surface subsidence caused by large-scale exploitation is increasingly serious. In high-water-level mining areas, subsidence often results in groundwater exposure, forming large-area seasonal or permanent water-bearing areas. This not only causes land resource waste, ecological environment destruction and man-ground contradiction, but also significantly changes the carbon circulation process of the region. The subsidence water accumulation area is used as a staggered zone of the amphibious ecosystem, the carbon source sink function is active, and the subsidence water accumulation area is a key link for influencing the net carbon balance of the mining area. The method accurately calculates the carbon reserves of the water body in the subsidence water accumulation area, and is basic work for evaluating the carbon fixation capacity of the mining area, promoting the construction of green mines and realizing the 'double carbon' target. However, current measuring methods face three major bottlenecks: the underwater topography and the reservoir capacity are difficult to determine, the topography of a ponding area is complex and continuously changes along with sedimentation, the traditional measuring method has high cost and long period, underwater topography data with high space-time resolution are difficult to obtain, reservoir capacity calculation is inaccurate, and the method becomes a main error source for estimating the carbon reserves. The limitation of static calculation of the carbon reserves is that the existing method is mostly limited to 'snapshot' estimation of the carbon reserves at a certain moment, cannot reflect the dynamic change of the carbon reservoir along with seasons, years and mining activities, cannot predict future trends, and is difficult to meet the requirement of adaptive management. The mechanism of the carbon circulation process is unknown, and the key processes of controlling the change of carbon reserves (such as runoff and sedimentation), output (such as gas discharge and outflow), internal conversion (photosynthesis and respiration) and the like are lack of quantitative description, so that the key knowledge of the mechanism of carbon sink function formation and regulation is unclear. Disclosure of Invention The invention aims to provide a method for measuring and calculating the carbon reserves of the water body in the subsidence ponding area for solving the problems. The invention provides a method for measuring and calculating carbon reserves of a water body in a subsidence ponding area, which comprises the following steps: Step 100, obtaining geological mining condition data of a mining area where the subsidence water area is located, and calculating a subsidence data set of the subsidence water area based on the geological mining condition data; step 200, acquiring a remote sensing image data sequence of a subsidence water accumulation area, preprocessing, extracting a water accumulation boundary data sequence of the subsidence water accumulation area by adopting a water body index method, and generating a water accumulation area data sequence based on the water accumulation boundary data sequence; Step 300, overlapping the ponding area data sequence and the sinking data in a unified coordinate system, and obtaining a reservoir capacity time sequence; step 400, arranging a plurality of sampling points in a subsidence water accumulation area system by adopting a grid method, collecting a plurality of water sample data, and acquiring water body organic carbon density data of the subsidence water accumulation area at a set moment based on the plurality of water sample data; And 500, calculating the water body carbon reserves at the set moment based on the reservoir capacity time sequence and the water body organic carbon density data of the subsidence water accumulation area at the set moment. As a further optimization scheme of the invention, the method further comprises the following steps: Step 600, constructing a carbon dynamic model, inputting the carbon reserves of the water body at the set moment in the step 500 into the carbon dynamic model, and outputting the change rate of the carbon reserves; step 700, generating dynamic water carbon reserves at set moments based on the carbon reserve change rate and the water carbon reserves at set moments. As a further optimization scheme of the invention, the geological mining condition data comprise average mining depth H, mining thickness m and coal seam inclination angleThe expected ru