CN-121093747-B - Method and device for determining drainage basin nitrogen output path and electronic equipment

CN121093747BCN 121093747 BCN121093747 BCN 121093747BCN-121093747-B

Abstract

The invention relates to the field of water quality management, in particular to a method and a device for determining a nitrogen output path of a river basin and electronic equipment. The method comprises the steps of obtaining initial continuous daily flow data and initial discrete nitrogen concentration data corresponding to an outlet of a target river basin in preset duration, expanding the initial discrete nitrogen concentration data to obtain target continuous nitrogen concentration data, separating the initial continuous daily flow data, determining target surface flow data, target underground flow data and target soil flow data in the initial continuous daily flow data, determining surface nitrogen concentration data, underground nitrogen concentration data and soil nitrogen concentration data from the target continuous nitrogen concentration data, and determining a nitrogen output path corresponding to the target river basin based on the surface nitrogen concentration data, the underground nitrogen concentration data and the soil nitrogen concentration data. The method for accurately determining the nitrogen output path under the sparse water quality data condition is provided.

Inventors

HUANG XUAN
XIE YUTING
XIA YONGQIU
ZHANG YUHANG
LI YUANHANG
REN YANYU

Assignees

河海大学
中国科学院南京土壤研究所

Dates

Publication Date: 20260512
Application Date: 20250805

Claims (9)

1. A method for determining a basin nitrogen output path, the method comprising: acquiring initial continuous daily flow data corresponding to an outlet of a target river basin in a preset time period and initial discrete nitrogen concentration data corresponding to the outlet of the target river basin; Expanding the initial discrete nitrogen concentration data based on the relation between the initial continuous daily flow data and the initial discrete nitrogen concentration data to obtain target continuous nitrogen concentration data; Separating the initial continuous daily flow data, and determining target surface flow data, target underground flow data and target soil flow data in the initial continuous daily flow data; Determining surface nitrogen concentration data, underground nitrogen concentration data and in-soil nitrogen concentration data from the target continuous nitrogen concentration data based on the target surface flow data, the target underground flow data and a relationship between the target in-soil flow data and the target continuous nitrogen concentration data; Determining a nitrogen output path corresponding to the target river basin based on the surface nitrogen concentration data, the underground nitrogen concentration data and the in-soil flowing nitrogen concentration data; The step of separating the initial continuous daily flow data and determining the target surface flow data, the target underground flow data and the target soil flow data in the initial continuous daily flow data comprises the following steps: acquiring a drainage basin type corresponding to the target drainage basin; Extracting a drainage basin characteristic parameter corresponding to the target drainage basin according to the drainage basin type, wherein the drainage basin characteristic parameter comprises a landform parameter, a land utilization parameter and a hydrological parameter; According to the drainage basin characteristic parameters, calculating a first initial filtering core parameter and a second initial filtering core parameter through a mapping model; Correcting the first initial filter core parameter and the second initial filter core parameter based on a preset parameter range to obtain a first target initial filter core parameter and a second target filter core parameter; And separating the initial continuous daily flow data based on the first target initial filtering core parameter and the second target filtering core parameter, and determining target surface flow data, target underground flow data and target soil flow data in the initial continuous daily flow data.
2. The method of claim 1, wherein expanding the initial discrete nitrogen concentration data based on a relationship between the initial continuous daily flow data and the initial discrete nitrogen concentration data to obtain target continuous nitrogen concentration data comprises: acquiring a drainage basin type corresponding to the target drainage basin; Acquiring a target parameter corresponding to the target river basin based on the river basin type, and calculating a correction coefficient corresponding to the target river basin based on the target parameter; Inputting the initial continuous daily flow data and the initial discrete nitrogen concentration data into an expansion model based on preset data, and expanding the initial discrete nitrogen concentration data to obtain standby continuous nitrogen concentration data; determining sampling time characteristics corresponding to each sub-data in the standby continuous nitrogen concentration data according to the time corresponding relation between the standby continuous nitrogen concentration data and the initial continuous daily flow data; Determining a time layering weight corresponding to each sub-data based on each sampling time characteristic; And correcting the standby continuous nitrogen concentration data based on the correction coefficient and the time layering weight to obtain the target continuous nitrogen concentration data.
3. The method according to claim 2, wherein the inputting the initial continuous daily flow data and the initial discrete nitrogen concentration data into an expansion model based on preset data expands the initial discrete nitrogen concentration data to obtain spare continuous nitrogen concentration data, comprises: Dividing a time axis according to the intensity of hydrologic events to obtain a stormwater flushing time interval, a conventional rainfall time interval and a rainless stable time interval; binding a concentration variation expected tag for each time interval; Determining an initial sub-data expansion model matched with the concentration change expected label corresponding to each time interval from the preset data expansion model aiming at the concentration change expected label corresponding to each time interval; Grading the initial discrete nitrogen concentration data according to data representativeness to obtain grade A nitrogen concentration data, grade B nitrogen concentration data and grade C nitrogen concentration data; Determining the A-level nitrogen concentration data and the B-level nitrogen concentration data as anchor point data, and determining weight information of each anchor point data according to a time interval corresponding to each anchor point data; according to the weight information of each anchor point data, calculating the weight duty ratio corresponding to each time interval; Calculating the sensitivity coefficient corresponding to each time interval according to the weight duty ratio corresponding to each time interval; Correcting initial parameters in the initial sub-data expansion model corresponding to each time interval based on the sensitivity coefficient to obtain a target sub-data expansion model corresponding to each time interval; and expanding the initial discrete nitrogen concentration data based on each target sub-data expansion model to obtain standby continuous nitrogen concentration data.
4. The method of claim 3, wherein the correcting the initial parameters in the initial sub-data expansion model corresponding to each time interval based on the sensitivity coefficient to obtain the target sub-data expansion model corresponding to each time interval includes: acquiring a reference time point corresponding to each anchor point data in each time interval; predicting predicted nitrogen concentration data corresponding to each reference time point in each time interval based on the initial sub-data expansion model corresponding to each time interval; calculating a first relation between the predicted nitrogen concentration data corresponding to each reference time point and the anchor point data corresponding to the reference time point in each time interval; And correcting initial parameters in the initial sub-data expansion model corresponding to each time interval based on the first relation and the sensitivity coefficient to obtain the target sub-data expansion model corresponding to each time interval.
5. The method of claim 1, wherein the separating the initial continuous daily flow data based on the first time target initial filter core parameters and the second time target filter core parameters, determining target surface flow data, target subsurface flow data, and target in-soil flow data in the initial continuous daily flow data, comprises: calculating a first daily acceleration of the initial continuous daily flow data; if the first daily acceleration is larger than a preset acceleration threshold, reducing the first target initial filter core parameters; Calculating initial surface flow data based on the reduced first time target initial filter core parameters; Calculating the second daily acceleration corresponding to the initial surface flow data; If the second daily acceleration rate is larger than the first daily acceleration rate, calculating initial underground flow data and initial in-soil flow data based on the second target initial filter core parameters; adding the calculated initial surface flow data, the initial underground flow data and the initial soil flow data to obtain predicted continuous flow data; If the ratio of the absolute value of the difference between the predicted continuous daily flow data and the initial continuous daily flow data to the initial continuous daily flow data is smaller than or equal to a preset ratio threshold value, determining the initial surface flow data, the initial underground flow data and the initial in-soil flow data as the target surface flow data, the target underground flow data and the target in-soil flow data; And if the ratio of the absolute value of the difference between the predicted continuous daily flow data and the initial continuous daily flow data to the initial continuous daily flow data is larger than the preset ratio threshold, correcting the initial surface flow data, the initial underground flow data and the initial soil flow data to obtain the target surface flow data, the target underground flow data and the target soil flow data.
6. The method of claim 1, wherein the determining surface nitrogen concentration data, subsurface nitrogen concentration data, and in-soil flow nitrogen concentration data from the target continuous nitrogen concentration data based on the target surface flow data, the target subsurface flow data, and a relationship between the target in-soil flow data and the target continuous nitrogen concentration data comprises: according to the corresponding relation between the target continuous nitrogen concentration data and the initial continuous daily flow data, selecting an underground flow candidate end member with a first ratio between the target underground flow data and the initial continuous daily flow data larger than a first preset ratio from the target continuous nitrogen concentration data; If the daily fluctuation of the nitrogen concentration corresponding to the candidate end members of the underground flow for a continuous preset number of days is smaller than a preset fluctuation threshold value, determining each candidate end member of the underground flow as an effective end member of the underground flow; According to the corresponding relation between the target continuous nitrogen concentration data and the initial continuous daily flow data, selecting an in-soil flow candidate end member with a second ratio smaller than a second preset ratio between the target surface flow data and the initial continuous daily flow data from first residual nitrogen concentration data, wherein the first residual nitrogen concentration data represents other nitrogen concentration data except the effective end member of the underground flow in the target continuous nitrogen concentration data; calculating the correlation between each in-soil flow candidate end member and target in-soil flow data corresponding to the in-soil flow candidate end members; Screening out the effective end members of the in-soil flow from the candidate end members of the in-soil flow according to the correlation; Determining each second residual nitrogen concentration data as a ground surface stream end member candidate set, wherein the second residual nitrogen concentration data characterizes other nitrogen concentration data except the effective end member of the underground stream and the effective end member of the soil middle stream in the target continuous nitrogen concentration data; Screening the ground surface flow data and the initial continuous daily flow data from the ground surface flow end member candidate set according to the corresponding relation between the target continuous nitrogen concentration data and the initial continuous daily flow data, wherein a third ratio between the target ground surface flow data and the initial continuous daily flow data is larger than a third preset ratio, and the flow is larger than a ground surface flow effective end member with a preset quantile; And determining the surface nitrogen concentration data, the underground nitrogen concentration data and the in-soil flow nitrogen concentration data according to the surface flow effective end member, the underground flow effective end member and the in-soil flow effective end member.
7. The method of claim 6, wherein said determining said surface nitrogen concentration data, said subsurface nitrogen concentration data, and said in-soil flow nitrogen concentration data from said surface flow effective end member, said subsurface flow effective end member, and said in-soil flow effective end member comprises: extracting characteristics of the effective end members of the underground flow to obtain an underground flow nitrogen concentration reference value and an underground flow nitrogen concentration stability coefficient; based on the characteristic extraction of the effective end member of the in-soil flow, obtaining an average value of the concentration of the in-soil flow nitrogen and a response slope of the in-soil flow; Extracting characteristics based on the effective end members of the surface flow to obtain a surface flow nitrogen concentration peak value and a nonlinear response index; calculating the underground nitrogen concentration data according to the underground flow nitrogen concentration reference value and the underground flow nitrogen concentration stability coefficient; calculating the data of the concentration of the in-soil flowing nitrogen according to the average value of the concentration of the in-soil flowing nitrogen and the response slope of the in-soil flowing nitrogen; and calculating the surface nitrogen concentration data according to the surface nitrogen concentration peak value and the nonlinear response index.
8. A basin nitrogen output path determination device, the device comprising: the acquisition module is used for acquiring initial continuous daily flow data corresponding to the outlet of the target river basin in a preset time period and initial discrete nitrogen concentration data corresponding to the outlet of the target river basin; The expansion module is used for expanding the initial discrete nitrogen concentration data based on the relation between the initial continuous daily flow data and the initial discrete nitrogen concentration data to obtain target continuous nitrogen concentration data; The separation module is used for separating the initial continuous daily flow data and determining target surface flow data, target underground flow data and target soil flow data in the initial continuous daily flow data; the method comprises the steps of separating initial continuous daily flow data, determining target surface flow data, target underground flow data and target in-soil flow data in the initial continuous daily flow data, acquiring a drainage basin type corresponding to the target drainage basin, extracting drainage basin characteristic parameters corresponding to the target drainage basin according to the drainage basin type, wherein the drainage basin characteristic parameters comprise landform parameters, land utilization parameters and hydrological parameters, calculating first initial filtering core parameters and second initial filtering core parameters through a mapping model according to the drainage basin characteristic parameters, correcting the first initial filtering core parameters and the second initial filtering core parameters based on a preset parameter range to obtain first target initial filtering core parameters and second target filtering core parameters, separating the initial continuous daily flow data based on the first target initial filtering core parameters and the second target filtering core parameters, and determining target surface flow data, target underground flow data and target in-soil flow data in the initial continuous daily flow data; The first determining module is used for determining surface nitrogen concentration data, underground nitrogen concentration data and in-soil nitrogen concentration data from the target continuous nitrogen concentration data based on the target surface flow data, the target underground flow data and the relation between the target in-soil flow data and the target continuous nitrogen concentration data; The second determining module is used for determining a nitrogen output path corresponding to the target river basin based on the surface nitrogen concentration data, the underground nitrogen concentration data and the in-soil flow nitrogen concentration data.
9. An electronic device, comprising: A memory and a processor in communication with each other, the memory having stored therein computer instructions, the processor executing the computer instructions to perform the basin nitrogen output path determination method of any one of claims 1 to 7.

Description

Method and device for determining drainage basin nitrogen output path and electronic equipment Technical Field The invention relates to the field of water quality management, in particular to a method and a device for determining a nitrogen output path of a river basin and electronic equipment. Background With the advancement of industrialization and urbanization, human activities have been aggravated by intervention of natural nitrogen circulation, and water eutrophication has become a global prominent environmental problem. The nitrogen is used as a key limiting nutrient element, the ecological health of the receiving water body is directly affected in the river basin conveying process, and the nitrogen pollution management and the water environment protection are based on the accurate cognition of a nitrogen conveying hydrologic path. Currently, sparsity of water quality monitoring data is a major technical challenge, particularly in less developed economic areas, remote mountainous areas and small and medium-sized watercourses. The traditional nitrogen output path analysis technology has obvious limitations that the hydrologic process separation technology adopts a two-component mode to ignore key intermediate paths such as soil shallow lateral flow and the like, the chemical tracing method has high requirements on sampling density and data quality, the accuracy is reduced when the data is sparse, and the load estimation method confuses hydrologic path separation and load attribution and relies on continuous complete water quality data. The core problems of the prior art include strong data dependence, incomplete path identification, poor sparse data processing capability and limited applicability, and difficulty in meeting the nitrogen pollution treatment requirement of a limited river basin of a monitoring resource. Therefore, there is a need to develop a technical method for precisely determining the nitrogen output path under the condition of sparse water quality data. Disclosure of Invention In view of the above, the invention provides a method and a device for determining a nitrogen output path of a river basin and electronic equipment, so as to provide a method for accurately determining the nitrogen output path under the condition of sparse water quality data. In a first aspect, the present invention provides a method for determining a drainage basin nitrogen output path, the method comprising: Acquiring initial continuous daily flow data corresponding to the outlet of a target river basin in a preset time period and initial discrete nitrogen concentration data corresponding to the outlet of the target river basin; Expanding the initial discrete nitrogen concentration data based on the relation between the initial continuous daily flow data and the initial discrete nitrogen concentration data to obtain target continuous nitrogen concentration data; separating the initial continuous daily flow data, and determining target surface flow data, target underground flow data and target soil flow data in the initial continuous daily flow data; Determining surface nitrogen concentration data, underground nitrogen concentration data and in-soil nitrogen concentration data from the target continuous nitrogen concentration data based on the target surface flow data, the target underground flow data and the relationship between the target in-soil flow data and the target continuous nitrogen concentration data; And determining a nitrogen output path corresponding to the target river basin based on the surface nitrogen concentration data, the underground nitrogen concentration data and the in-soil flowing nitrogen concentration data. In an alternative implementation mode, the method comprises the steps of expanding initial discrete nitrogen concentration data based on a relation between the initial continuous daily flow data and the initial discrete nitrogen concentration data to obtain target continuous nitrogen concentration data, obtaining a basin type corresponding to the target basin, obtaining target parameters corresponding to the target basin based on the basin type, calculating correction coefficients corresponding to the target basin based on the target parameters, inputting the initial continuous daily flow data and the initial discrete nitrogen concentration data into an expansion model based on preset data to expand the initial discrete nitrogen concentration data to obtain standby continuous nitrogen concentration data, determining sampling time characteristics corresponding to each sub-data in the standby continuous nitrogen concentration data according to a time corresponding relation between the standby continuous nitrogen concentration data and the initial continuous daily flow data, determining time layering weights corresponding to each sub-data based on the sampling time characteristics, and correcting the standby continuous nitrogen concentration data based on the correction coefficients