CN116562051B - Land sea nitrogen and phosphorus load trend estimation method - Google Patents
Land sea nitrogen and phosphorus load trend estimation method Download PDFInfo
- Publication number
- CN116562051B CN116562051B CN202310768295.6A CN202310768295A CN116562051B CN 116562051 B CN116562051 B CN 116562051B CN 202310768295 A CN202310768295 A CN 202310768295A CN 116562051 B CN116562051 B CN 116562051B
- Authority
- CN
- China
- Prior art keywords
- load
- phosphorus
- nitrogen
- sea
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 title claims abstract description 236
- 229910052698 phosphorus Inorganic materials 0.000 title claims abstract description 146
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 title claims abstract description 145
- 239000011574 phosphorus Substances 0.000 title claims abstract description 145
- 229910052757 nitrogen Inorganic materials 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004088 simulation Methods 0.000 claims abstract description 33
- 238000006243 chemical reaction Methods 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 16
- 241001465754 Metazoa Species 0.000 claims abstract description 14
- 235000015097 nutrients Nutrition 0.000 claims description 70
- 239000000126 substance Substances 0.000 claims description 37
- 239000002352 surface water Substances 0.000 claims description 34
- YUWBVKYVJWNVLE-UHFFFAOYSA-N [N].[P] Chemical compound [N].[P] YUWBVKYVJWNVLE-UHFFFAOYSA-N 0.000 claims description 25
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 23
- 238000009792 diffusion process Methods 0.000 claims description 20
- 239000002689 soil Substances 0.000 claims description 8
- 238000009304 pastoral farming Methods 0.000 claims description 3
- NNJPGOLRFBJNIW-HNNXBMFYSA-N (-)-demecolcine Chemical compound C1=C(OC)C(=O)C=C2[C@@H](NC)CCC3=CC(OC)=C(OC)C(OC)=C3C2=C1 NNJPGOLRFBJNIW-HNNXBMFYSA-N 0.000 claims description 2
- 238000003306 harvesting Methods 0.000 claims description 2
- 230000014759 maintenance of location Effects 0.000 claims description 2
- 238000011002 quantification Methods 0.000 claims description 2
- 230000000717 retained effect Effects 0.000 claims description 2
- 239000003344 environmental pollutant Substances 0.000 abstract description 14
- 231100000719 pollutant Toxicity 0.000 abstract description 14
- 239000003653 coastal water Substances 0.000 abstract description 7
- 238000005094 computer simulation Methods 0.000 abstract description 5
- 238000011160 research Methods 0.000 description 8
- 239000010865 sewage Substances 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 210000003608 fece Anatomy 0.000 description 6
- 238000013508 migration Methods 0.000 description 6
- 230000005012 migration Effects 0.000 description 6
- MMDJDBSEMBIJBB-UHFFFAOYSA-N [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] Chemical compound [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] MMDJDBSEMBIJBB-UHFFFAOYSA-N 0.000 description 5
- 229930002868 chlorophyll a Natural products 0.000 description 5
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 239000013049 sediment Substances 0.000 description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000010871 livestock manure Substances 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 238000011144 upstream manufacturing Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000011324 bead Substances 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 244000144972 livestock Species 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 239000002028 Biomass Substances 0.000 description 1
- 241000283690 Bos taurus Species 0.000 description 1
- 244000068988 Glycine max Species 0.000 description 1
- 235000010469 Glycine max Nutrition 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 241001494479 Pecora Species 0.000 description 1
- 241000282887 Suidae Species 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 description 1
- 238000009360 aquaculture Methods 0.000 description 1
- 244000144974 aquaculture Species 0.000 description 1
- 238000004178 biological nitrogen fixation Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012851 eutrophication Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 238000011010 flushing procedure Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 239000010842 industrial wastewater Substances 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000813 microbial effect Effects 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 244000144977 poultry Species 0.000 description 1
- 238000009374 poultry farming Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 235000013311 vegetables Nutrition 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a land sea nitrogen and phosphorus load trend estimation method, which relates to the field of pollutant management, and comprises the steps of obtaining statistical data of different data types; the data types include land use portion, economic portion, population portion, agricultural portion, animal husbandry portion, river basin information, and reservoir information; determining total nitrogen and phosphorus load into the sea by utilizing a nitrogen and phosphorus load flexible estimation model according to the statistical data; performing time scale conversion on the total nitrogen and phosphorus load entering the sea to obtain the nitrogen and phosphorus concentration of each sea entrance; and determining the space-time distribution and biochemical trend of Liu Yuandan phosphorus output load after sea entry by using a estuary coast numerical simulation model according to the nitrogen and phosphorus concentration of each sea entry. The invention conveniently and efficiently realizes the whole process load estimation and dynamic simulation of Liu Yuandan phosphorus load output-sea entrance confluence-estuary coastal water environment trend.
Description
Technical Field
The invention relates to the field of pollutant management, in particular to a land sea nitrogen and phosphorus load trend estimation method.
Background
Dissolved inorganic nitrogen (Dissolved Inorganic Nitrogen, DIN) and dissolved inorganic phosphorus (Dissolved Inorganic Phosphorus, DIP) are the main nutrients supporting phytoplankton growth and reproduction and maintaining stable water ecosystem, but the input of excessive nitrogen and phosphorus load easily causes estuary coastal water eutrophication. Along with the rapid development of social economy of coastal cities, high-strength human activities output a large amount of nitrogen, phosphorus and other nutrient substances for downstream estuary water areas, and Liu Yuandan phosphorus output load becomes the most main external contribution source of nutrient substance circulation of the estuary coastal water areas, so that the estimation of Liu Yuandan phosphorus output load is particularly important for downstream water ecological management and protection.
The method for estimating Liu Yuandan phosphorus output load is generally divided into two types: firstly, a measuring method is adopted, reference points and control points are arranged at different positions of a river, even monitoring points are arranged at a sewage outlet, and the output load of nutrient substances of the river is calculated according to monitoring data; and secondly, the method can carry out mathematical modeling according to statistical data or estimate the nitrogen and phosphorus output load by means of the existing model, and has low time and economic cost and wide application. Depending on the model complexity, empirical models such as ECM (Export Coefficient Model), conceptual models such as Global NEWS 2 (Global Nutrient Export from WaterSheds 2)、MARINA (Model to Assess River Inputs of Nutrients to seAs) , mechanism models such as SWAT (oil AND WATER ASSESSMENT Tool) and the like can be classified.
The biochemical trend of Liu Yuandan phosphorus loading into the sea generally comprises migration and diffusion in the horizontal direction, movement in the vertical direction, deposit burial and release at the water body-sediment interface, microbial denitrification, participation in primary production and the like, and has important significance in promoting the circulation of bio-geochemical substances and maintaining the stability of a water ecological system. Because the process is numerous and the calculation is complex, the method is usually realized by means of a numerical simulation model, the numerical simulation model is formed by formulating physical phenomena of a real environment, and the method is more suitable for large-scale area research and is more widely applied. The representative numerical simulation models include FVCOM, EFDC, MIKE and the like, the calculation principles are approximately the same, structured or unstructured grids are used as calculation units, and the hydrodynamic water quality equation is solved.
The invention discloses a quantitative identification method for river basin nitrate nitrogen source river entering load and river entering coefficient, which comprises the following steps: the method comprises the following steps of (1) collecting and analyzing isotope traceable samples; (2) Calculating the contribution proportion of the river nitrate nitrogen source in the flow field; (3) Basin outlet runoff monitoring and nitrate nitrogen sample collection and analysis; (4) calculating the total river load of the nitrate nitrogen at the outlet of the river basin; (5) And (5) calculating river basin nitrate nitrogen source river entering load and river entering coefficient.
The invention discloses an agricultural non-point source phosphorus pollution estimation method based on soil attribute space distribution, which comprises the following three steps: (1) The establishment of the response relationship between the soil attribute space distribution and the non-point source pollution load comprises the steps of typical small research area selection, soil attribute space distribution, typical area non-point source pollution SWAT model simulation and the establishment of the response relationship between the soil attribute and the non-point source pollution load space distribution; (2) collecting and detecting soil samples in the area to be estimated; and (3) estimating the pollution load of the agricultural non-point source.
The invention discloses a method for estimating pollutant load of a coastal zone, which mainly comprises the following two steps: (1) Determining a pollution source and main pollutants, and determining a point source pollution source and a non-point source pollution source according to the pollutant formation characteristics and the data acquisition conditions of a research area; determining main pollutants according to statistics of pollutant content; point source pollution including industrial wastewater; the non-point source pollution refers to water pollution caused by various pollutants which are converged into the receiving water body under the action of large-area precipitation and runoff flushing; (2) Pollution load estimation, including point source pollution load estimation and non-point source pollution load estimation; the main pollutants are total nitrogen, total phosphorus, chemical oxygen demand and ammonia nitrogen.
CN104951986A discloses a river basin agriculture non-point source pollutant lake inlet load estimation method, which comprises the following steps: acquiring typical small watershed data in a target watershed, constructing a typical small watershed SWAT, and acquiring Load Lr and Load Lo in the typical small watershed according to the SWAT; acquiring the pollution discharge coefficients of the planting industry source, the livestock and poultry farming industry source, the aquaculture industry source and the rural living source of each sub-basin of the typical small basin, and acquiring Load Ls in the typical small basin by combining basic information survey data; obtaining a target river basin agricultural non-point source pollutant ditch reduction coefficient Factor cr according to Load Ls and Load Lr; acquiring Length Lr according to the water system data; acquiring a river channel reduction coefficient Factor rr of the agricultural non-point source pollutant in the target river basin according to Load Lr、LoadLo and Length Lr; acquiring the agricultural non-point source nitrogen and phosphorus emission load of a target river basin and the total Length Br of a river channel in the target river basin; and obtaining the load of the agricultural non-point source pollutant entering the lake of the target river basin according to the nitrogen and phosphorus discharge load of the agricultural non-point source of the target river basin, the factors cr and rr and the Length Br.
Mayorga et al set up a basin nitrogen-phosphorus output load model (Global NEWS 2) with a Global scale (resolution of 1 degree x1 degree), the pollutants are divided into point sources and diffusion sources, basin boundaries are divided according to a geographic information model, then hydrologic information of the basin is obtained through hydrologic model simulation, nutrient input and related parameters are obtained from different databases and documents, and the dissolved inorganic and organic nitrogen-phosphorus output load of the Global main basin is effectively estimated.
Strokal et al further perfects the Global NEWS 2 model according to the characteristics of the Chinese river basin, proposes a MARINA model, corrects related parameters and updates data sources aiming at the characteristics of the Chinese river basin, supplements the nitrogen and phosphorus output load contributed by human excreta which is not connected with a sewage system, and successfully realizes the estimation of the nitrogen and phosphorus sea load of the six river basins in the sub-river basin scale and the prediction of the nitrogen and phosphorus sea load in future situations.
Qin the content and degradation coefficient of TOC, TN, TP in the Liaohua wetland sediment are obtained through field sampling and laboratory analysis, a MIKE model is coupled with a built-in ECO Lab module to simulate hydrodynamic force-water quality of the Liaohua wetland, and the result is compared with the spatial-temporal distribution and migration trend of the Liaohua wetland in different hydrologic periods TOC, TN, TP on the sediment surface layer and the sediment bottom layer.
The general nitrogen and phosphorus load estimation experience model is a black box model, so that mechanism interpretation of results is difficult, the data demand of the mechanism model is large, when a weather and soil database is incomplete, a plurality of parameters are difficult to acquire and determine, more time and economic cost are required to be input, and the cost performance is low; the biochemical trend simulation of the traditional sea-entering nitrogen-phosphorus load often does not incorporate the process of generating and converging Liu Yuandan phosphorus load into calculation, the close connection between an upstream river basin and a downstream estuary coast is cut, and the migration and transformation process and potential ecological influence after Liu Yuandan phosphorus load is output to the estuary coast water area cannot be estimated; in addition, similarly, the source of sea nitrogen and phosphorus load and the main contribution source of simulation results cannot be resolved by adopting in-situ sampling or literature data as numerical simulation model input.
Therefore, a method for conveniently and efficiently realizing the whole process load estimation and dynamic simulation of Liu Yuandan phosphorus load output-inlet seaport confluence-estuary coastal water environment trend is needed.
Disclosure of Invention
The invention aims to provide a land sea nitrogen and phosphorus load returning estimation method which can conveniently and efficiently realize the whole process load estimation and dynamic simulation of Liu Yuandan phosphorus load output-sea port confluence-estuary coastal water environment returning.
In order to achieve the above object, the present invention provides the following solutions:
a land sea nitrogen and phosphorus load trend estimation method comprises the following steps:
Acquiring statistical data of different data types; the data types include land use portion, economic portion, population portion, agricultural portion, animal husbandry portion, river basin information, and reservoir information;
determining total nitrogen and phosphorus load into the sea by utilizing a nitrogen and phosphorus load flexible estimation model according to the statistical data;
performing time scale conversion on the total nitrogen and phosphorus load entering the sea to obtain the nitrogen and phosphorus concentration of each sea entrance;
And determining the space-time distribution and biochemical trend of Liu Yuandan phosphorus output load after sea entry by using a estuary coast numerical simulation model according to the nitrogen and phosphorus concentration of each sea entry.
Optionally, determining the total nitrogen-phosphorus load into the sea by using a nitrogen-phosphorus load flexible estimation model according to the statistical data specifically comprises the following steps:
According to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of nutrient substances output from a diffusion source and a point source to the surface water body;
according to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of the nutrient substances output from the surface water body to the outlet of the sub-watershed;
According to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of the nutrient substances output from the outlet of the sub-river basin to the river mouth;
And determining the total nitrogen and phosphorus load entering the sea according to the load of the nutrient substances output from the diffusion source and the point source to the surface water body, the load of the nutrient substances output from the surface water body to the outlet of the sub-river basin and the load of the nutrient substances output from the outlet of the sub-river basin to the river mouth.
Optionally, the time scale conversion is performed on the total nitrogen and phosphorus load entering the sea to obtain the nitrogen and phosphorus concentration of each sea entrance, which specifically comprises the following steps:
Determining the sea-entering nitrogen and phosphorus concentration according to the sea-entering total nitrogen and phosphorus load and the daily sea-entering flow data of the river basin;
And distributing the sea-entering nitrogen and phosphorus concentration according to the runoff ratio to obtain the nitrogen and phosphorus concentration of each sea-entering mouth.
Optionally, the estuary coast numerical simulation model is MIKE a 21 FM model.
Optionally, the spatial-temporal distribution is the concentration of the dissolved inorganic nitrogen and the dissolved inorganic phosphorus at different positions at different moments; the biochemical trend is the growth and propagation process of phytoplankton participated by the dissolved inorganic nitrogen and the dissolved inorganic phosphorus.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
The invention acquires the statistical data of different data types; the data types include land use portion, economic portion, population portion, agricultural portion, animal husbandry portion, river basin information, and reservoir information; determining total nitrogen and phosphorus load into the sea by utilizing a nitrogen and phosphorus load flexible estimation model according to the statistical data; performing time scale conversion on the total nitrogen and phosphorus load entering the sea to obtain the nitrogen and phosphorus concentration of each sea entrance; and determining the space-time distribution and biochemical trend of Liu Yuandan phosphorus output load after sea entry by using a estuary coast numerical simulation model according to the nitrogen and phosphorus concentration of each sea entry. The overall process load estimation and dynamic simulation of Liu Yuandan phosphorus load output-inlet coastal confluence-estuary coastal water environment trend can be conveniently and efficiently realized through the nitrogen and phosphorus load flexible estimation model and the estuary coastal numerical simulation model.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of a Liu Haidan phosphorus load trend estimation method provided by the invention;
FIG. 2 is a flow chart of a Liu Haidan phosphorus load trend estimation method provided by the invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention aims to provide a land sea nitrogen and phosphorus load returning estimation method which can conveniently and efficiently realize the whole process load estimation and dynamic simulation of Liu Yuandan phosphorus load output-sea port confluence-estuary coastal water environment returning.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
As shown in FIG. 1 and FIG. 2, the land sea nitrogen and phosphorus load trend estimation method provided by the invention can be divided into three parts, namely Liu Yuandan phosphorus output load estimation, sea entrance nitrogen and phosphorus load distribution and nitrogen and phosphorus load trend simulation after sea entrance, and comprises the following steps:
Step 101: acquiring statistical data of different data types; the data types include land use, economic, population, agricultural, animal husbandry, watershed, and reservoir information.
Step 102: and determining the total nitrogen and phosphorus load entering the sea by utilizing a nitrogen and phosphorus load flexible estimation model according to the statistical data.
Step 102, specifically includes: according to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of nutrient substances output from a diffusion source and a point source to the surface water body; according to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of the nutrient substances output from the surface water body to the outlet of the sub-watershed; according to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of the nutrient substances output from the outlet of the sub-river basin to the river mouth; and determining the total nitrogen and phosphorus load entering the sea according to the load of the nutrient substances output from the diffusion source and the point source to the surface water body, the load of the nutrient substances output from the surface water body to the outlet of the sub-river basin and the load of the nutrient substances output from the outlet of the sub-river basin to the river mouth.
Liu Yuandan phosphorus output load estimation is realized by means of a nitrogen-phosphorus load flexible estimation model (Flexible Model to Access Nutrient Loads from Basins, F-MANLB) provided by the invention. The model is improved based on MARINA 2.0, the characteristics of a conceptual model (combination of a mechanism formula and an empirical coefficient, moderate data demand and accuracy) are reserved, and three parts of contents are mainly improved based on MARINA 2.0: 1) The drainage basin hydrologic information is extracted and generalized, and pretreatment such as grading, trend, confluence and the like of the river in the dividing area can be omitted; 2) Correcting hierarchical relationships among partial sub-watershed and partial formula parameters according to literature and related data to make them more realistic, e.g. correcting the downstream outlet of east river as a bead triangle watershed rather than a bead river mouth and formulating the formulaThe correction of 0.8 in (2) is 0.95 (the maximum nitrogen and phosphorus removal rate achieved by the four-stage treatment of the existing sewage treatment plant is 95%); wherein,The ratio of inorganic nitrogen in the urban population excreta, to which the sewage system is connected for the point source, to surface water, wherein FE represents the coefficient, pnt represents the point source, DIN represents the inorganic nitrogen in dissolved form, hum represents the human, con represents the connection (sewage system), urb represents the town; /(I)Is the removal rate of N in the sewage system treatment process in the sub-basin j. 3) The computer programming language packaging model formula is adopted, the whole 4-component design (comprising an input file, a model main body, an output file and a configuration file) is adopted, wherein the input file, the output file and the configuration file are all in a well-known Excel workbook format, and the data can be conveniently and efficiently arranged and edited.
The required data types include the following major classes: 1) The land utilization part is required to collect a river basin land utilization type data set, and uses geographic information system software to carry out mask extraction, reclassifying and statistics on each sub-river basin/city to obtain an agricultural land occupation ratio Agr and a vegetation occupation ratio Nat; 2) The economic part is the domestic total production value (GDP); 3) The population portion includes town resident population Pop urb and country resident population Pop rur; 4) The agricultural part comprises various crop planting areas (rice, corn, wheat, soybean, vegetables and the like are respectively expressed by Area 1、Area2、Area3; 5) The animal husbandry part comprises the feeding amount of various livestock (pigs, cattle, sheep, poultry and the like, which are respectively expressed by Ani 1、Ani2、Ani3. Cndot.); 6) The basin information includes natural flow (surface runoff) Q nat, water consumption Q con, and actual flow at basin exit Q act of the basin; 7) The reservoir information includes the reservoir flow rate Q res, the volume V, and the depth H of all the tributaries or main flows within the flow area.
The F-MANLB model quantitatively estimates the nitrogen and phosphorus loads discharged by different sources in a flow field by defining the flow field to form sub-flow fields, wherein the specific forms of nitrogen and phosphorus are divided into DIN and DIP, the discharge forms are divided into a diffusion source and a point source, and the diffusion source comprises animal manure (ani) of animal husbandry, chemical fertilizer application (fe), resident excrement (hum. Unconon) without a sewage system, atmospheric nitrogen sedimentation (dep) and biological nitrogen fixation (bio) of crops and natural vegetation; point source animal husbandry animal manure (ani), resident excreta (hum. Uncon. Url), unconnected and connected to a sewage system.
The F-MANLB model estimates the nitrogen and phosphorus output load of a target river basin through three main progressive steps, wherein the first step is to quantify the load of nutrient substances output from a pollution source to an surface water body, the second step is to quantify the load of nutrient substances output from the surface water body to a sub-river basin outlet, and the third step is to quantify the load of nutrient substances output from the sub-river basin outlet to a river mouth. And outputting the calculated results of the three steps and the nitrogen and phosphorus load amounts discharged by different sources in different forms. The general formula is shown as equation (1), and the main steps are shown as formulas (2) to (6).
The general formula:
（1）
Wherein M F.y.j is the load (kg) of nutrient substances (F: DIN, DIP) output to the estuary; RS F.y.j is the load (kg/a) of nutrient F from the pollution source y to the surface water body; FE riv.F.outlet.j is the ratio (0-1) of RS F.y.j output to sub-basin j outlet; FE riv.F.mouth.j is The ratio of the output to the estuary (0-1). F is a nutrient substance, and is herein classified into DIN and DIP according to the form. In addition, y represents a source such as animal manure, human excrement, etc., and different sources may be discharged in two forms: the diffusion source and the point source are denoted by dif and pnt, respectively. j represents a sub-basin (a single basin may be represented if the sub-basin is not distinguished).
Step 1: quantification of the load of nutrients from diffuse and point sources to surface water RS F.y.j:
（2）
（3）
Wherein RSdif F.y.j is nutrient F is the load (kg) output from sub-basin j diffusion source y to surface water; WSdif E.y.j is nutrient element (E: N, P) is load (kg/a) output from sub-basin j diffusion source y into the ground; g F.j is the proportion (0-1) of nutrient (F) remaining in the soil after application to the agricultural land by animal grazing and crop harvesting; FE ws.F.j is the ratio (0-1) of nutrient (F) to surface water entering basin (j). Similarly RSpnt F.y.j is nutrient F is the load (kg) output from sub-basin j point source y to surface water; RSpnt E.y.j is nutrient E is the load (kg/a) of surface water output from the point source y of the sub-basin j; FEpnt F.j is the ratio of nutrient F to surface water in sub-basin j (0-1). RS F.y.j in equation (1) is the sum of RSdif F.y.j and RSpnt F.y.j.
Step 2: quantifying the load FE riv.F.oulet.j of nutrients from the surface water output to the outlet of the sub-basin:
（4）
Wherein D F.j is the ratio (0-1) of nutrient substances (DIN, DIP) retained in the reservoir of the sub-basin j; l F.j is the ratio (0-1) of nutrient (DIN, DIP) retention or loss in the subflow j water system; FQrem j is the ratio (0-1) of nutrients (DIN, DIP) in the sub-basin j water system to be removed by the various types of water.
Step 3: quantifying the load of nutrient output from the outlet of the sub-basin to the estuary (FE riv.F.mouth.j):
（5）
（6）
Wherein, FE riv.F.mouth.jmT、FEriv.F.mouth.jmC、FEriv.F.mouth.jdC is the proportion of nutrient substance F output to the river mouth from the middle stream branch jmT, the middle stream main stream jmC and the downstream main stream jdC respectively; jmTFEriv.F.outlet.jmC、jmTFEriv.F.outlet.jdC The ratio of nutrient F output from the downstream tributary jmT outlet to the downstream main stream jmC and downstream main stream jdC, respectively; jmCFEriv.F.outlet.jdC Is the ratio of nutrient F output from the midstream main stream jmC to the downstream main stream jdC.
Step 103: and performing time scale conversion on the total nitrogen and phosphorus load entering the sea to obtain the nitrogen and phosphorus concentration of each sea entrance.
Step 103, specifically includes: determining the sea-entering nitrogen and phosphorus concentration according to the sea-entering total nitrogen and phosphorus load and the daily sea-entering flow data of the river basin; and distributing the sea-entering nitrogen and phosphorus concentration according to the runoff ratio to obtain the nitrogen and phosphorus concentration of each sea-entering mouth.
Load distribution of nitrogen and phosphorus at sea entrance: because the estuary coastal numerical simulation model uses the grid as a calculation unit, the requirements on the terrain precision are higher, the coastal contour and the subdivision degree of the estuary have a certain influence on the simulation precision, and it is necessary to distinguish the main estuary and assign an input time sequence file for each estuary. In addition, the output result of the F-MANLB model is total nitrogen and phosphorus load (kg/a) of the sea in the time scale of year, and the estuary coast numerical simulation model generally needs a time sequence file with a finer scale, in order to simulate the sea entry trend of the land source nitrogen and phosphorus load in the year, the output result of the F-MANLB model needs to be converted to match the input precision requirement of the estuary coast numerical simulation model, and the following provides a conversion method of day scale data:
1) Since the human social production activities are less affected by weather, the daily nitrogen and phosphorus output load is considered to be the same, i.e., the daily nitrogen and phosphorus output load is the total nitrogen and phosphorus load M F divided by the number of days (365 days, for example). The daily sea intake flow data D of the river basin is collected, and the daily nitrogen and phosphorus output load is divided by the flow to obtain the sea intake nitrogen and phosphorus concentration C F.
2) And determining the sea inlet runoff ratio Fr i of each sea inlet according to literature and other data, and distributing the sea inlet nitrogen and phosphorus concentration according to the runoff ratio to obtain the nitrogen and phosphorus concentration C F.i of each sea inlet. As shown in the formulas (7) to (8):
（7）
（8）
Wherein, C F is the concentration (mg/l) of nitrogen and phosphorus entering the sea; m F is total nitrogen-phosphorus load (kg/a) into the sea; d is basin inflow (10 8 m3）;10-5 represents concentration unit conversion; subscript F represents nitrogen and phosphorus (DIN, DIP) of different forms; C F.i is nitrogen and phosphorus concentration (mg/l) of each sea inlet; fr i is sea inflow runoff ratio (0-1) of each sea inlet; i is sea inlet label (i=1, 2, 3).
Step 104: and determining the space-time distribution and biochemical trend of Liu Yuandan phosphorus output load after sea entry by using a estuary coast numerical simulation model according to the nitrogen and phosphorus concentration of each sea entry. The estuary coast numerical simulation model is constructed according to actual requirements by a two-dimensional water simulation platform MIKE and FM which are developed by Denmark DHI company and take unstructured grids as calculation units. The space-time distribution is the concentration of the dissolved inorganic nitrogen and the dissolved inorganic phosphorus at different positions at different moments; the biochemical trend is the growth and propagation process of phytoplankton participated by the dissolved inorganic nitrogen and the dissolved inorganic phosphorus.
Nitrogen-phosphorus loading sea-going simulation: the result M F.y.j obtained by the F-MANLB model is the DIN and DIP total output load (kg) of the current year of the target river basin, the daily DIN and DIP concentrations C F.i at different sea inlets at the downstream of the river basin can be obtained after the conversion and distribution of the unit in the step 2, and the DIN and DIP concentrations C F.i are used as data input of a MIKE numerical simulation model after being manufactured into a time series file.
The trend of nitrogen and phosphorus into sea backlashes means that dissolved inorganic nitrogen and dissolved inorganic phosphorus participate in the biological geochemical processes of phytoplankton growth and reproduction (phytoplankton biomass is expressed as chlorophyll a concentration, mug/l), diffusion, sedimentation and the like, and the biological geochemical processes are specifically characterized by the concentration space-time distribution of chlorophyll a (Chlorophyll-a, chl-a), DIN and DIP. The model is mainly divided into 5 parts of data collection, data introduction, model parameter setting, model debugging and final model operation, and the specific method comprises the following steps:
1) Data collection and time series file production. And importing a flow and water quality time series file as an offshore boundary. Determining the boundary of the research area, and collecting the runoff amount of each sea entrance at the upstream boundary, the tide level of the downstream boundary, the wind speed, the wind direction and the topography or water depth data of the research area. In addition, water quality data such as temperature, salinity, dissolved inorganic nitrogen, dissolved inorganic phosphorus, chlorophyll a, and dissolved oxygen (Dissolved Oxygen, DO) are also required to be collected. The different types of data will be produced as time series files, respectively, as data input at the model boundaries. After finishing, all types of data are divided into two parts, one part is used for preliminary simulation (data set 1), the other part is used for calibration and verification (data set 2), and time series files are manufactured for each sea entrance and open sea boundary according to the model format requirement. In addition, the daily nitrogen and phosphorus concentrations C F.i for each sea entrance obtained in step 103 were prepared as DIN and DIP time series files, respectively (dataset 3). Each type of data is a time series file, such as a DIN time series file, containing DIN concentrations (mg/l) at the same time of day at different boundaries. And drawing a boundary contour, importing a water depth topographic dataset and performing spatial interpolation to obtain the integral topography of the research area. Generating grids, optimizing and adjusting the number and the size of the grids, and deriving grid files.
2) Setting a model time range and a step length in a MIKE 21 FM model master setting file, and importing a grid file, a wind field time sequence file (the wind speed and the wind direction of monitoring points in a research area are respectively m/s and DEG) and a roughness file (the Manning roughness coefficient of each grid is calculated according to the terrain grid file, and the unit is m 1/3/s).
3) The open sea boundary is imported into a tidal time series file, and relevant parameters are set. Since MIKE, model 21 and FM are modular in design, only the hydrodynamic module and ECO Lab ecological simulation module are used here, and data set 1 is employed. In the hydrodynamic module, a flow time series file (in m n/s) is respectively imported for each sea entrance at an upstream boundary, and a tide level time series file (in m) is imported for the open sea at a downstream boundary. In order to ensure the accuracy of the simulation result, a salt temperature submodule in the hydrodynamic module is started, and temperature and salinity time series files (with the units of the temperature and the PSU respectively) are imported for all sea inlets and the open sea boundaries. Different templates can be selected according to requirements in the ECO Lab module or related formulas can be written by the ECO Lab module, wherein a 'EutrophicationModel 1' template is taken as an example: after the template file is imported, DIN, DIP, DO and Chl-a time series files are imported for each sea entrance, and the data of the open sea boundary and other required inputs such as non-actual measurement data and referents are set as constants. Finally, setting nitrogen-phosphorus migration conversion coefficients such as nitrogen-phosphorus diffusion rate, deposition rate, phytoplankton absorptivity and the like.
4) Running the model, model rating and verification is performed using another portion of the data (dataset 2). The evaluation method uses the correlation coefficient (R 2) and Nash coefficient (NSE), where R 2 is close to 1 and NSE is between 0-1 (the closer to 1 the better), the model results are considered acceptable, otherwise the correlation parameters are readjusted. The nitrogen-phosphorus migration conversion coefficient refers to parameters such as nitrogen-phosphorus diffusion rate, deposition rate, phytoplankton absorptivity, etc., and may be a globally uniform constant, or a function of other parameters (flow rate, temperature, salinity, etc.). The specific usage is that in step 3), after the ECO Lab module is imported into the MIKE 21: 21 FM model, the parameters are set in the model master profile.
5) After model calibration is completed, DIN and DIP input at each sea entrance in the ECO Lab module are changed into corresponding F-MANLB Liu Yuandan phosphorus output load model estimation results (data set 3), and the estuary coast numerical simulation model is operated again to obtain space-time distribution and biochemical trend of Liu Yuandan phosphorus output load after sea entrance, wherein the space-time distribution of nitrogen and phosphorus after sea entrance is reflected as concentration (mg/l) of DIN and DIP at different positions at different moments, the biochemical trend mainly means that DIN and DIP participate in the growth and propagation process of phytoplankton, and the concentration space-time distribution of chlorophyll a (Chl-a) is expressed as (mug/l), namely the concentration of DIN, DIP and Chl-a at different positions at different moments in a research area.
The invention provides a land sea nitrogen and phosphorus load trend estimation system, which comprises:
The acquisition module is used for acquiring statistical data of different data types; the data types include land use, economic, population, agricultural, animal husbandry, watershed, and reservoir information.
And the sea-entering total nitrogen and phosphorus load determining module is used for determining the sea-entering total nitrogen and phosphorus load by utilizing a nitrogen and phosphorus load flexible estimation model according to the statistical data.
And the time scale conversion module is used for performing time scale conversion on the total nitrogen and phosphorus load entering the sea to obtain the nitrogen and phosphorus concentration of each sea entrance.
And the simulation module is used for determining the space-time distribution and biochemical trend of Liu Yuandan phosphorus output load after the sea is input according to the nitrogen and phosphorus concentration of each sea entrance by using a river mouth coast numerical simulation model.
As an optional implementation manner, the sea-going total nitrogen-phosphorus load determining module specifically includes:
And the first load determining unit is used for quantifying the load of the nutrient substances output from the diffusion source and the point source to the surface water body by using the nitrogen-phosphorus load flexible estimation model according to the statistical data.
And the second load determining unit is used for quantifying the load of the nutrient substances output from the surface water body to the outlet of the sub-watershed by using the nitrogen-phosphorus load flexible estimation model according to the statistical data.
And the third load determining unit is used for quantifying the load of the nutrient substances output from the outlet of the sub-river basin to the river mouth by using the nitrogen-phosphorus load flexible estimation model according to the statistical data.
And the sea-entering total nitrogen and phosphorus load determining unit is used for determining the sea-entering total nitrogen and phosphorus load according to the load of the nutrient substances output from the diffusion source and the point source to the surface water body, the load of the nutrient substances output from the surface water body to the sub-basin outlet and the load of the nutrient substances output from the sub-basin outlet to the estuary.
As an alternative embodiment, the time scale conversion module specifically includes:
And the sea-entering nitrogen and phosphorus concentration determining unit is used for determining sea-entering nitrogen and phosphorus concentration according to the sea-entering total nitrogen and phosphorus load and the daily sea-entering flow data of the river basin.
And the distribution unit is used for distributing the sea-entering nitrogen and phosphorus concentration according to the runoff ratio to obtain the nitrogen and phosphorus concentration of each sea-entering mouth.
According to the invention, a Liu Yuandan phosphorus output load estimation conceptual model is coupled on the basis of a estuary coast numerical simulation model, good balance is achieved in terms of data demand, accuracy and mechanism interpretation performance, consistency of Liu Yuandan phosphorus output load and sea nitrogen and phosphorus input load is ensured, the method has good applicability to plain river network areas, accuracy and universality of estimation simulation of Liu Haidan phosphorus load migration and conversion processes in different areas are effectively improved, and more scientific and comprehensive references can be provided for land and sea overall planning management such as coastal urban social and economic development and estuary coast ecological environment protection.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In view of the foregoing, this description should not be construed as limiting the invention.
Claims (4)
1. A land sea nitrogen and phosphorus load trend estimation method is characterized by comprising the following steps:
Acquiring statistical data of different data types; the data types include land use portion, economic portion, population portion, agricultural portion, animal husbandry portion, river basin information, and reservoir information;
determining total nitrogen and phosphorus load into the sea by utilizing a nitrogen and phosphorus load flexible estimation model according to the statistical data;
performing time scale conversion on the total nitrogen and phosphorus load entering the sea to obtain the nitrogen and phosphorus concentration of each sea entrance;
determining the space-time distribution and biochemical trend of Liu Yuandan phosphorus output load after sea entry by using a estuary coast numerical simulation model according to the nitrogen and phosphorus concentration of each sea entry;
The method for determining the total nitrogen and phosphorus load into the sea by utilizing the nitrogen and phosphorus load flexible estimation model according to the statistical data comprises the following steps:
According to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of nutrient substances output from a diffusion source and a point source to the surface water body;
according to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of the nutrient substances output from the surface water body to the outlet of the sub-watershed;
According to the statistical data, the nitrogen-phosphorus load flexible estimation model is utilized to quantify the load quantity of the nutrient substances output from the outlet of the sub-river basin to the river mouth;
determining total nitrogen and phosphorus load entering sea according to the load of the nutrient substances output from the diffusion source and the point source to the surface water body, the load of the nutrient substances output from the surface water body to the outlet of the sub-river basin and the load of the nutrient substances output from the outlet of the sub-river basin to the river mouth;
quantification of the load of nutrients output from diffuse and point sources to surface water rsf.y.j:
Wherein rsdiff.y.j is the load of nutrient F output from sub-basin j diffusion source y to surface water; wsdife.y.j is the load that the nutrient element is output into the land from sub-basin j diffusion source y; gf.j is the proportion of nutrients remaining in the soil after application to the agricultural land, through animal grazing and crop harvesting; fews.f.j is the proportion of nutrients entering the surface water of the basin; rspntf.y.j is the load of nutrient F output from sub-basin j point source y to surface water; rspnte.y.j is the load of nutrient element E to surface water output from sub-basin j point source y; fepntf.j is the ratio of nutrient F to sub-basin j surface water; the RSF.y.j is the sum of RSdifF.y.j and RSpntF.y.j;
Quantifying the amount of load FEriv.F.oulet.j of nutrients output from the surface water body to the outlet of the sub-basin:
wherein DF.j is the ratio of nutrient substances retained in the reservoir of the sub-basin j; LF.j is the ratio of nutrient retention or loss in the subflow j water system; FQremj is the ratio of nutrients removed by various types of water in the sub-basin j water system;
quantifying the load quantity FEriv.F.mouth.j of nutrient output from the sub-basin outlet to the estuary:
Wherein, FE riv.F.mouth.jmT、FEriv.F.mouth.jmC is the proportion of nutrient substance F output from midstream branch flow jmT and midstream trunk flow jmC to river mouth; jmTFEriv.F.outlet.jmC、jmTFEriv.F.outlet.jdC The ratio of nutrient F output from the downstream tributary jmT outlet to the downstream main stream jmC and downstream main stream jdC, respectively; jmCFEriv.F.outlet.jdC Is the ratio of nutrient F output from the midstream main stream jmC to the downstream main stream jdC.
2. The Liu Haidan phosphorus load trend estimation method according to claim 1, wherein the time scale conversion is performed on the total nitrogen and phosphorus load of the sea to obtain the nitrogen and phosphorus concentration of each sea entrance, specifically including:
Determining the sea-entering nitrogen and phosphorus concentration according to the sea-entering total nitrogen and phosphorus load and the daily sea-entering flow data of the river basin;
And distributing the sea-entering nitrogen and phosphorus concentration according to the runoff ratio to obtain the nitrogen and phosphorus concentration of each sea-entering mouth.
3. The Liu Haidan phosphorus load trend estimation method of claim 1, wherein the estuary coastal numerical simulation model is a MIKE model FM model.
4. The Liu Haidan phosphorus load trend estimation method according to claim 1, wherein the spatial-temporal distribution is the concentration of dissolved inorganic nitrogen and dissolved inorganic phosphorus at different positions at different moments; the biochemical trend is the growth and propagation process of phytoplankton participated by the dissolved inorganic nitrogen and the dissolved inorganic phosphorus.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310768295.6A CN116562051B (en) | 2023-06-28 | 2023-06-28 | Land sea nitrogen and phosphorus load trend estimation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310768295.6A CN116562051B (en) | 2023-06-28 | 2023-06-28 | Land sea nitrogen and phosphorus load trend estimation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116562051A CN116562051A (en) | 2023-08-08 |
| CN116562051B true CN116562051B (en) | 2024-04-26 |
Family
ID=87496647
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310768295.6A Active CN116562051B (en) | 2023-06-28 | 2023-06-28 | Land sea nitrogen and phosphorus load trend estimation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116562051B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117892983B (en) * | 2024-03-14 | 2024-05-28 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | Method and system for determining offshore area land pollution load distribution |
| CN118225995B (en) * | 2024-05-22 | 2024-08-02 | 宁波瀚陆海洋科技有限公司 | Bay nitrogen and phosphorus pollutant monitoring method and system |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102663267A (en) * | 2012-05-15 | 2012-09-12 | 南京大学 | Method for determining drainage basin surface source pollution load of semi-humid region |
| CN103449656A (en) * | 2012-05-30 | 2013-12-18 | 东莞理工学院 | Sequence fed-batch reactor and biological aerated filter (SFBR-BAF) biological nitrogen and phosphorus removal technology |
| CN103473463A (en) * | 2013-09-17 | 2013-12-25 | 中国环境科学研究院 | Method of quantitative determination of background concentration of nitrogen and phosphorus of water body of lake basins |
| CN104951986A (en) * | 2015-06-03 | 2015-09-30 | 中国农业科学院农业资源与农业区划研究所 | Estimation method for load of watershed agricultural non-point source pollutants into lake |
| CN106250666A (en) * | 2016-06-22 | 2016-12-21 | 中国水利水电科学研究院 | A kind of irrigation system ecumene lake nitrogen and phosphorus loading analogy method |
| CN108647401A (en) * | 2018-04-19 | 2018-10-12 | 北京师范大学 | A kind of basin nitrogen and phosphorus pollution appraisal procedure based on space remote sensing technology |
| CN209442809U (en) * | 2018-09-06 | 2019-09-27 | 岭南生态文旅股份有限公司 | A kind of riverway water quality purification system for coping with point-source pollution |
| CN114021829A (en) * | 2021-11-11 | 2022-02-08 | 东莞理工学院 | Land use pattern prediction and optimization method considering non-point source pollution control |
| CN114580167A (en) * | 2022-03-02 | 2022-06-03 | 广西大学 | Method for analyzing river-entering load of nitrogen and phosphorus from sediment source by real-time monitoring and CSSI technology |
| CN115641028A (en) * | 2022-10-31 | 2023-01-24 | 同济大学 | Plain river network surface source pollution load estimation method based on revenue and expenditure model |
| CN115809605A (en) * | 2022-12-19 | 2023-03-17 | 浙江大学 | Prediction method for farmland runoff nitrogen and phosphorus loss concentration and runoff |
-
2023
- 2023-06-28 CN CN202310768295.6A patent/CN116562051B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102663267A (en) * | 2012-05-15 | 2012-09-12 | 南京大学 | Method for determining drainage basin surface source pollution load of semi-humid region |
| CN103449656A (en) * | 2012-05-30 | 2013-12-18 | 东莞理工学院 | Sequence fed-batch reactor and biological aerated filter (SFBR-BAF) biological nitrogen and phosphorus removal technology |
| CN103473463A (en) * | 2013-09-17 | 2013-12-25 | 中国环境科学研究院 | Method of quantitative determination of background concentration of nitrogen and phosphorus of water body of lake basins |
| CN104951986A (en) * | 2015-06-03 | 2015-09-30 | 中国农业科学院农业资源与农业区划研究所 | Estimation method for load of watershed agricultural non-point source pollutants into lake |
| CN106250666A (en) * | 2016-06-22 | 2016-12-21 | 中国水利水电科学研究院 | A kind of irrigation system ecumene lake nitrogen and phosphorus loading analogy method |
| CN108647401A (en) * | 2018-04-19 | 2018-10-12 | 北京师范大学 | A kind of basin nitrogen and phosphorus pollution appraisal procedure based on space remote sensing technology |
| CN209442809U (en) * | 2018-09-06 | 2019-09-27 | 岭南生态文旅股份有限公司 | A kind of riverway water quality purification system for coping with point-source pollution |
| CN114021829A (en) * | 2021-11-11 | 2022-02-08 | 东莞理工学院 | Land use pattern prediction and optimization method considering non-point source pollution control |
| CN114580167A (en) * | 2022-03-02 | 2022-06-03 | 广西大学 | Method for analyzing river-entering load of nitrogen and phosphorus from sediment source by real-time monitoring and CSSI technology |
| CN115641028A (en) * | 2022-10-31 | 2023-01-24 | 同济大学 | Plain river network surface source pollution load estimation method based on revenue and expenditure model |
| CN115809605A (en) * | 2022-12-19 | 2023-03-17 | 浙江大学 | Prediction method for farmland runoff nitrogen and phosphorus loss concentration and runoff |
Non-Patent Citations (1)
| Title |
|---|
| Qiangqiang Rong, Meirong Su,et al..Simulation of dissolved nutrient export from the Dongjiang river basin with a grid-based NEWS model."ICEMEE 2018".2018,第1-4页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116562051A (en) | 2023-08-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116562051B (en) | Land sea nitrogen and phosphorus load trend estimation method | |
| Aazami et al. | Monitoring and assessment of water health quality in the Tajan River, Iran using physicochemical, fish and macroinvertebrates indices | |
| Wagner et al. | Landscape drivers of regional variation in the relationship between total phosphorus and chlorophyll in lakes | |
| Pantus et al. | Quantifying and evaluating ecosystem health: A case study from Moreton Bay, Australia | |
| Kim et al. | Simulation of algal bloom dynamics in a river with the ensemble Kalman filter | |
| Arheimer et al. | Watershed modelling of nonpoint nitrogen losses from arable land to the Swedish coast in 1985 and 1994 | |
| Alexander et al. | Estimates of diffuse phosphorus sources in surface waters of the United States using a spatially referenced watershed model | |
| Costa et al. | Choosing an appropriate water quality model—A review | |
| Pai et al. | Identifying priority subwatersheds in the Illinois River drainage area in Arkansas watershed using a distributed modeling approach | |
| Mouri et al. | Using land-use management policies to reduce the environmental impacts of livestock farming | |
| Pedde et al. | Modeling sources of nutrients in rivers draining into the Bay of Bengal—a scenario analysis | |
| Fatema et al. | Spatial variation of water quality parameters in a mangrove estuary | |
| Petihakis et al. | Ecosystem modelling: Towards the development of a management tool for a marine coastal system part-II, ecosystem processes and biogeochemical fluxes | |
| Xu et al. | Estimate of nutrient sources and transport into Bohai Bay in China from a lower plain urban watershed using a SPARROW model | |
| CN114462698B (en) | Phosphorus emission pollution load prediction method for drainage basin catchment area | |
| Yuan et al. | Defining the sources and the fate of nitrate by using dual isotopes and a Bayesian isotope mixing model: Water–nitrate management in cascade dams of Lancang river | |
| Srinivasan et al. | Development of a national-scale framework to characterise transfers of N, P and Escherichia coli from land to water | |
| Wei et al. | Nutrient transport and transformation in macrotidal estuaries of the French Atlantic coast: a modeling approach using the Carbon-Generic Estuarine Model | |
| Shibata et al. | Comparison of spatial and temporal changes in riverine nitrate concentration from terrestrial basins to the sea between the 1980s and the 2000s in Japan: Impact of recent demographic shifts | |
| Forbes et al. | Physical factors control phytoplankton production and nitrogen fixation in eight Texas reservoirs | |
| Donyaii | Evaluation of climate change impacts on the optimal operation of multipurpose reservoir systems using cuckoo search algorithm | |
| Hasanyar et al. | How much do bacterial growth properties and biodegradable dissolved organic matter control water quality at low flow? | |
| Povilaitis | Source apportionment and retention of nutrients and organic matter in the Merkys river basin in southern Lithuania | |
| Mouri et al. | Estimation of total nitrogen transport and retention during flow in a catchment using a mass balance model incorporating the effects of land cover distribution and human activity information | |
| Pearce et al. | Metabolic regimes of three mid-order streams in southern Ontario, Canada exposed to contrasting sources of nutrients |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |