A water quality response simulation method under different inflow and outflow scheduling scenarios based on an IWIND-LR model

Abstract

The application relates to a quantitative simulation method for water quality response of a target water body under different inflow and outflow scheduling scenarios based on an IWIND-LR model, and belongs to the technical field of water resource management. The method comprises the following steps: establishing a hydrodynamic-water quality-ecological model of the target water body based on an IWIND-LR model; adjusting a momentum equation, a water quality equation and an algae equation in the model based on the hydrodynamic and water quality mechanism processes corresponding to the inflow and outflow scheduling, specifically comprising adding an inflow flow term and an outflow flow term in the momentum equation, adding an inflow water quality term and an outflow water quality term in the water quality equation, and adding an inflow algae term and an outflow algae term in the algae equation, so as to quantitatively express the influence of the inflow and outflow on the hydrodynamic and water quality processes of the target water body; coupling the adjusted momentum equation, water quality equation and algae equation with a hydrodynamic-water quality-ecological simulation framework of the IWIND-LR model, generating an inflow and outflow scheduling model of the target water body through boundary condition setting, specifically comprising rebuilding a gridded hydrodynamic simulation module, a gridded water quality simulation module and a gridded ecological simulation module based on the adjusted equations, and coupling based on the simulation framework of the IWIND-LR model; formulating an inflow and outflow scheduling scenario scheme; simulating the different inflow and outflow scheduling scenarios through the inflow and outflow scheduling model, and generating water quality response simulation results of the target water body under each scheduling scenario. The method described in the application can realize quantitative simulation of water quality response of the target water body under different inflow and outflow scheduling scenarios, and provide scientific support for managers to formulate scheduling strategies.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of water environment numerical simulation, in particular to a quantitative simulation method for water quality response of a target water body under different inflow and outflow scheduling scenarios based on an IWIND-LR model, which is suitable for water quality simulation analysis of inflow and outflow scheduling scenarios of rivers, reservoirs, lakes and other water areas. BACKGROUND

[0002] In the daily management of lakes, reservoirs and the like, water body inflow and outflow scheduling is a control means that comprehensively considers hydrology, water power, water quality and other factors, and through water diversion facilities, clear water is introduced into the lake to improve the water environment of the lake, so as to realize the sustainable development of regional social economy and lake ecological environment, which is the most commonly used and most convenient management method. Through scientific scheduling of inflow and outflow, water resources can be reasonably utilized, the water quality environment of the water body can be improved, and the problems of water resource shortage and pollution can be alleviated. Different inflow and outflow scheduling will have different effects on the water quality of the water body, and the effects are not only related to the inflow and outflow flow, but also related to meteorological factors and other comprehensive factors. For a long time, subjective qualitative analysis has been mainly used, and there is a lack of scientific and systematic quantitative basis, which greatly reduces the reliability and feasibility of scheduling decision-making. The lake is a highly complex nonlinear dynamic system, so it is impossible to rely only on experience and qualitative estimation to simulate and evaluate the effects of inflow and outflow scheduling on the water quality of the water body.

[0003] In order to accurately evaluate the effects of different inflow and outflow scheduling schemes on the water quality of the water body, the present application provides a three-dimensional hydrodynamic-water quality model capable of simulating the migration and transformation process of pollutants in the lake body as a scientific and technological means to numerically simulate different inflow and outflow scheduling scenarios, thereby forming a quantitative simulation method for water quality response of a target water body under different inflow and outflow scheduling scenarios based on an IWIND-LR model. SUMMARY

[0004] In order to make the scheduling management of the water body more scientific and reasonable, the present application provides a quantitative simulation method for water quality response of a target water body under different inflow and outflow scheduling scenarios based on an IWIND-LR model.

[0005] In order to achieve the above-mentioned object, the present application adopts the following technical solutions:

[0006] A quantitative simulation method for water quality response of a target water body under different inflow and outflow scheduling scenarios based on an IWIND-LR model, comprising the following steps:

[0007] Step 1, building a hydrodynamic-water quality-ecological model of the target water body based on the IWIND-LR model;

[0008] Step 2, based on the water dynamics and water quality mechanism processes corresponding to the inflow and outflow scheduling, adjusting the momentum equation, water quality equation and algae equation in the above model, specifically including adding inflow and outflow flow terms in the momentum equation, adding inflow and outflow water quality terms in the water quality equation, and adding inflow and outflow algae terms in the algae equation, quantitatively expressing the influence of inflow and outflow on the water dynamics and water quality processes of the target water body;

[0009] Step 3, coupling the adjusted momentum equation, water quality equation and algae equation with the water dynamics-water quality-water ecology simulation framework of the IWIND-LR model, generating the inflow and outflow scheduling model of the target water body through boundary condition setting, specifically including rebuilding the grid-based water dynamics simulation module, grid-based water quality simulation module and grid-based water ecology simulation module based on the adjusted equation, and coupling based on the simulation framework of the IWIND-LR model;

[0010] Step 4, developing inflow and outflow scheduling scenario schemes;

[0011] Step 5, simulating the different inflow and outflow scheduling scenarios through the inflow and outflow scheduling model, generating the water quality response simulation results of the target water body corresponding to each scheduling scenario.

[0012] Through the above technical solution, combined with the water dynamics-water quality-water ecology model built based on the IWIND-LR model, the grid-discretized inflow and outflow scheduling simulation model of the target water body is established, which can quantitatively analyze the water quality impact of the scheduling scenario on the target water body by setting the location, operation mode and other scenario parameters of the inflow and outflow scheduling scenario in the inflow and outflow scheduling scenario evaluation process.

[0013] As a preferred, the above step 1, the water dynamics-water quality-water ecology model of the target water body built based on the IWIND-LR model, comprises the following sub-steps, characterized in that:

[0014] Step 1-1, discretizing the target water body into a plurality of grid cells based on the IWIND-LR model, considering that using curved grid can better match the boundary shape of the water body, and improving the calculation efficiency of the model under the condition of ensuring spatial accuracy, therefore, in the horizontal direction, the target water body is divided into a plurality of orthogonal curved grids; for each orthogonal curved grid, cutting into a plurality of grid cells in the vertical direction according to the sigma coordinate, and performing terrain interpolation on the grid cells according to the edge terrain data of the target water body.

[0015] Step 1-2, building the water dynamics model, water quality model and water ecology model of the target water body based on the IWIND-LR model and the grid cells;

[0016] Step 1-3, coupling the hydrodynamic model, the water quality model and the water ecological model based on the IWIND-LR model to generate a hydrodynamic-water quality-water ecological model of the target water body.

[0017] Based on the above technical solution, the dynamic processes of the hydrodynamics, water quality and water ecology inside the water body are numerically reproduced based on the IWIND-LR model, the dynamic changes inside the water body are tracked in all directions, so as to establish a quantitative relationship between various scheduling scenarios and the responses of the hydrodynamics, water quality and water ecology inside the water body, and provide a scientific decision basis for water body management.

[0018] As preferred, in the step 2, the hydrodynamic and water quality mechanism processes corresponding to the inflow and outflow scheduling are used to adjust the momentum equation, the water quality equation and the algae equation in the above model, specifically including adding the inflow flow term and the outflow flow term in the momentum equation, adding the inflow water quality term and the outflow water quality term in the water quality equation, and adding the inflow algae term and the outflow algae term in the algae equation, so as to quantitatively express the influence of the inflow and outflow on the hydrodynamic and water quality processes of the target water body, and the following sub-steps are included, characterized in that:

[0019] Step 2-1, in the momentum equation in the hydrodynamic-water quality-water ecological model, the inflow flow term and the outflow flow term quantitatively expressing the influence of the inflow and outflow scheduling processes on the hydrodynamics of the target water body are added;

[0020] In implementation, when adjusting the momentum equation in the hydrodynamic-water quality-water ecological model, the inflow flow term and the outflow flow term of the inflow and outflow scheduling scenarios can be added in the momentum equation, for example, the original momentum equation can be:

[0021]

[0022] Wherein, mx and my respectively represent the square root values of the components of the metric tensor along the diagonal direction; m = mx * my constitutes the Jacobian determinant or the determinant formed by the square root values of the metric tensor; u, v and w respectively represent the velocity components along the x, y and z coordinate axis directions; H represents the total depth (H = ζ + h); f refers to the Coriolis force parameter; p represents the pressure; Av refers to the longitudinal turbulent velocity or the eddy viscosity, and Qu and Qv are other momentum source and sink terms except for the inflow and outflow scheduled by the inflow and outflow scheduling scenarios.

[0023] The adjusted momentum equation can be:

[0024]

[0025] Wherein, Qu and Qv are other momentum source, sink items except for the inflow and outflow scheduling scenario scheduling inflow and outflow; Qwdu and Qwdv respectively represent the inflow and outflow scheduling scenario scheduling inflow flow items along the x coordinate axis direction and along the y coordinate axis direction in the model; Qreu and Qrev respectively represent the inflow and outflow scheduling scenario scheduling outflow flow items along the x coordinate axis direction and along the y coordinate axis direction in the model.

[0026] Step 2-2, in the water quality equation in the hydrodynamic-water quality-water ecological model, add inflow water quality items and outflow water quality items quantitatively expressing the influence of the inflow and outflow scheduling process on the water quality of the target water body.

[0027] In the water quality equation in the hydrodynamic-water quality-water ecological model, add the inflow and outflow water quality items of the inflow and outflow scheduling scenario, including chemical oxygen demand, total phosphorus, total nitrogen and dissolved oxygen indicators.

[0028] (1) Chemical oxygen demand

[0029] The original chemical oxygen demand equation can be:

[0030] COD = OC RPOC + OC LPOC + OC DOC + ∑ x＝c,d,g OC Bx + SOD …… (5); wherein: COD is the chemical oxygen demand, converted into the amount of oxygen (mg / L). OC RPOC is the oxygen consumption of inert particulate organic carbon; OC LPOC is the oxygen consumption of unstable particulate organic carbon; OC DOC is the oxygen consumption of dissolved organic carbon; OC Bx is the oxygen consumption of algae x (subscript x is c, d, g, corresponding to different algae); SOD is the oxygen consumption of endogenous organic carbon (all in O2, g / m 3 ).

[0031] The control equation of the oxygen consumption of inert and unstable particulate organic carbon is:

[0032]

[0033] The control equation of the oxygen consumption of dissolved organic carbon is:

[0034]

[0035] Wherein, OC RPOC , OC LPOC , OC DOC are the oxygen consumptions of inert particulate organic carbon, unstable particulate organic carbon and dissolved organic carbon (in O2, g / m 3); FCRP, FCLP, and FCDP represent the fractions of preyed inert particulate organic carbon, unstable particulate organic carbon, and dissolved organic carbon, respectively; K RPOC K LPOC The dissolution rates (d) of inert particulate organic carbon and unstable particulate organic carbon are respectively. -1 ); WS RP With WS LP , representing the settling velocities (m / d) of inert particulate organic carbon and unstable particulate organic carbon, respectively; WRPOC, WLPOC, and WDOC represent the external loads (gC / d) of inert particulate organic carbon, unstable particulate organic carbon, and dissolved organic carbon, respectively, excluding the inflow and outflow scheduling scenarios; FCD x KHR represents the fraction of dissolved organic carbon secreted by algae during metabolism when dissolved oxygen concentration is infinitely high. x DO represents the half-saturated concentration of dissolved oxygen corresponding to dissolved organic carbon excreted by algae; DO is the dissolved oxygen concentration (gO2 / m³). 3 ); K HR The heterotrophic respiration rate of dissolved organic carbon (d) -1 Denit is the rate of denitrification (d). -1 ); R OC The carbon-to-oxygen ratio is 2.67 g O2 / g C.

[0036] The governing equation for oxygen consumption of algae x (where subscripts x represent c, d, and g, corresponding to different algae) is:

[0037] Where: subscript x represents c, d, g, corresponding to different species of algae; OC Bx Oxygen consumption of algae x (in O2, g / m³) 3 ); t is time (d); P x Productivity (d) -1 ); BM x Metabolic rate (d) -1 ); PR x Predation rate (d) -1 ); WS x Settlement velocity (m / d); WB x To schedule other algal exogenous loads (gC / d) besides the inflow / outflow scheduling scenario; V is the computational grid volume (m³). 3 ), R OC The carbon-to-oxygen ratio is 2.67 g O2 / g C.

[0038] The governing equation for endogenous organic carbon oxygen consumption is:

[0039]

[0040] In the formula, SOD is the oxygen consumption of endogenous organic carbon (expressed as O2, g / m³). 3 ), DO is the dissolved oxygen concentration (g / m³) 3 ), KH SOD The half-saturation constant of dissolved oxygen required for the oxidation of endogenous organic carbon (in O2, g / m³). 3 KSOD is the oxidation rate of endogenous organic carbon (d). -1 BFSOD is the deposition flux of endogenous organic carbon (in O2 terms, g / (m³)). 2 ·d). WSOD is the external load of other endogenous organic carbon (as O2, g / d) other than the inflow and outflow scheduling scenarios.

[0041] The adjusted water quality equation can be:

[0042]

[0043] Among them, WRPOC, WLPOC, and WDOC represent inert particulate organic carbon, unstable particulate organic carbon, and dissolved organic carbon, respectively, representing other external loads (gC / d) excluding inflow and outflow scheduling scenarios; WB x The external algal load (gC / d) is the amount of algae generated outside of the inflow / outflow scheduling scenario; WSOD is the external load of endogenous organic carbon (g / d as O2) outside of the inflow / outflow scheduling scenario. WRPOC 调度入流 WLPOC 调度入流 WDOC 调度入流 WB x调度入流 The inflow load (gC / d) for inert particulate organic carbon, unstable particulate organic carbon, dissolved organic carbon, and different species of algae inflow and outflow scheduling scenarios are WSOD. 调度入流 Schedule inflow load (in O2, g / d) for inflow and outflow scheduling scenarios of endogenous organic carbon; WRPOC 调度出流 WLPOC 调度出流 WDOC 调度出流 WB x调度出流 The scenarios represent inert particulate organic carbon, unstable particulate organic carbon, dissolved organic carbon, and different species of algae, respectively, for scheduling outflow loads (gC / d) and WSOD. 调度出流 The outflow load (in terms of O2, g / d) is scheduled for the inflow and outflow scheduling scenario of endogenous organic carbon.

[0044] (2) Total phosphorus

[0045] The original total phosphorus equation can be:

[0046] TP=RPOP+LPOP+DOP+PO4……………………(16);

[0047] Where: TP is total phosphorus (mg / L), and RPOP is the concentration of inert particulate organic phosphorus (g·P / m³). 3 ); LPOP is the concentration of unstable particulate organic phosphorus (g·P / m³). 3 DOP is the concentration of soluble organic phosphorus (gP / m³). -3 ); PO4t represents total phosphate (gP / m -3 ).

[0048] The kinetic equations for inert and unstable particulate organophosphorus compounds are as follows:

[0049]

[0050]

[0051] Wherein: RPOP is the concentration of inert particulate organic phosphorus (g·P / m³). 3 ); LPOP is the concentration of unstable particulate organic phosphorus (g·P / m³). 3 FPRx represents the fraction of inert particulate organic phosphorus produced by algal metabolism; FPLx represents the fraction of unstable particulate organic phosphorus produced by algal metabolism; FPRP represents the fraction of predatory inert particulate organic phosphorus produced; FPLP represents the fraction of predatory unstable particulate organic phosphorus produced; APC represents the average phosphorus-to-carbon ratio of all algae (gP / gC); KRPOP represents the hydrolysis rate of inert particulate organic phosphorus (day). -1 KLPOP represents the hydrolysis rate (days) of unstable particulate organophosphorus compounds. -1 WRPOP represents the external load (P / day) of inert particulate organophosphorus compounds other than those scheduled inflow / outflow under the inflow / outflow scheduling scenario; WLPOP represents the external load (P / day) of unstable particulate organophosphorus compounds other than those scheduled inflow / outflow under the inflow / outflow scheduling scenario.

[0052] The kinetic equation for soluble organophosphorus compounds is:

[0053]

[0054] Where: DOP is the concentration of soluble organophosphorus compounds (gP / m³) -3 FPDx represents the fraction of soluble organic phosphorus produced by algal metabolism; FPDPx represents the fraction of soluble organic phosphorus produced by predation; KDOP represents the rate of soluble organic phosphorus mineralization (day). -1 ); WDOP refers to the external load of dissolved organic phosphorus (gP / m³) other than the inflow / outflow scheduling scenario. -3 ).

[0055] The kinetic equation for total phosphate is:

[0056]

[0057] Wherein: PO4t is total phosphate (gP / m -3 The composition is PO4d + PO4p; PO4d is soluble phosphate (gP / m -3 ); PO4p is particulate (adsorbed) phosphate (gP / m -3 FPIx represents the fraction of inorganic phosphorus produced by algal metabolism; FPIP represents the fraction of inorganic phosphorus produced by predation; WSTSS represents the suspended solids settling velocity (m / day), provided by the hydraulic model; BFPO4d represents the phosphate sediment-water exchange coefficient (gP / m³). 2 / day), exists only at the bottom layer; WPO4t(WPO4p+WPO4d) is the total phosphate external load (gP / day) other than the inflow and outflow scheduling scenario.

[0058] The adjusted water quality equation can be:

[0059]

[0060]

[0061] Wherein, WRPOP represents the external load of inert particulate organophosphorus compounds (P / day) other than those under the inflow / outflow scheduling scenario; WLPOP represents the external load of unstable particulate organophosphorus compounds (P / day) other than those under the inflow / outflow scheduling scenario; and WDOP represents the external load of soluble organophosphorus compounds (gP / m³) other than those under the inflow / outflow scheduling scenario. -3 WPO4t(WPO4p+WPO4d) represents the total external phosphate load (gP / day) other than the inflow and outflow scheduling scenarios. 调度入流 WLPOP 调度入流 The inflow load (gP / day) is scheduled under inflow and outflow scheduling scenarios for inert particulate organic phosphorus and unstable particulate organic phosphorus, respectively. 调度入流 Adjusting the inflow load (gP / m³) for soluble organic phosphorus -3 ), WPO4 调度入流 Scheduling inflow load (gP / day) for total phosphate inflow and outflow scheduling scenarios; WRPOP 调度出流 WLPOP 调度出流 The outflow load (gP / day) is scheduled under inflow and outflow scheduling scenarios for inert particulate organic phosphorus and unstable particulate organic phosphorus, respectively. 调度出流 Adjusting the outflow load (gP / m³) for soluble organic phosphorus inflow / outflow scheduling scenarios -3 ), WPO4 调度出流The outflow load (gP / day) is scheduled for the total phosphate inflow and outflow scheduling scenario.

[0062] (3) Total nitrogen

[0063] The original total nitrogen equation can be:

[0064] TN=RPON+LPON+DON+NH4+NO3……………………(25);

[0065] Wherein: TN is total nitrogen (mg / L), and RPON is the concentration of inert particulate organic nitrogen (g·N / m³). 3 ); LPON is the concentration of unstable particulate organic nitrogen (g·N / m³). 3 DON represents the concentration of dissolved organic nitrogen (gN / m³). -3 ); NH4 represents the ammonia nitrogen concentration (gN / m³). -3 NO3 represents the nitrate nitrogen concentration (gN / m³). -3 ).

[0066] The kinetic equations for inert and unstable organic nitrogen are as follows:

[0067]

[0068]

[0069] Wherein: RPON is the concentration of inert particulate organic nitrogen (g·N / m³). 3 ); LPON is the concentration of unstable particulate organic nitrogen (g·N / m³). 3 FNRx represents the fraction of inert particulate organic nitrogen produced by algal metabolism; FNLx represents the fraction of unstable particulate organic nitrogen produced by algal metabolism; FNRP represents the fraction of organic nitrogen produced by algae predating inert particulate organic nitrogen; FNLP represents the fraction of organic nitrogen produced by algae predating unstable particulate organic nitrogen; ANCx represents the N / C ratio (gN / gC) of a certain algae (x); KRPON represents the hydrolysis rate of inert particulate organic nitrogen (day). -1 KLPON represents the hydrolysis rate of unstable particulate organic nitrogen (days). -1 WRPON represents inert particulate organic nitrogen external loads (gN / day) other than those scheduled inflow / outflow in the inflow / outflow scheduling scenario; WLPON represents unstable particulate organic nitrogen external loads (gN / day) other than those scheduled inflow / outflow in the inflow / outflow scheduling scenario.

[0070] The kinetic equation for dissolved organic nitrogen is:

[0071]

[0072] Where: DON is the concentration of dissolved organic nitrogen (gN / m³) 3);FNDx is the fraction of dissolved organic nitrogen produced by the metabolism of algae x; FNDP is the fraction produced by predation of dissolved organic nitrogen; KDON is the mineralization rate of dissolved organic nitrogen (day -1 );BFDON is the exchange of dissolved organic nitrogen that only occurs at the deep bottom layer (gC / m 2 / day); WDON is the external load of other dissolved organic nitrogen except for the inflow and outflow scheduled in the inflow and outflow scheduling scenario (gN / day).

[0073] The kinetic equation for ammonia nitrogen is as follows:

[0074]

[0075] Where: FNIx is the fraction of inorganic nitrogen produced by the metabolism of algae; FNIP is the fraction produced by the predation of inorganic nitrogen by algae; PNx is the preference coefficient of algae for ammonia nitrogen uptake (0 < PNx < 1); KNit is the nitrification rate (day -1 );BFNH4 is the sediment-water ammonia nitrogen exchange that only occurs at the deep bottom layer (gN / m 2 / day); WNH4 is the external load of other ammonia nitrogen except for the inflow and outflow scheduled in the inflow and outflow scheduling scenario (gN / day).

[0076] The kinetic equation for nitrate nitrogen is as follows:

[0077]

[0078] Where: ANDC is the mass of nitrate nitrogen consumed per unit mass of oxidized dissolved organic carbon (0.933 gN / gC); BFNO3 is the sediment-water nitrate exchange that only occurs at the deep bottom layer (gN / m 2 / day); WNO3 is the external load of other nitrates except for the inflow and outflow scheduled in the inflow and outflow scheduling scenario (gN / day).

[0079] The adjusted water quality equation can be:

[0080]

[0081]

[0082] Among them, WRPON represents the external source load of inert particulate organic nitrogen (gN / day) other than the inflow / outflow scheduling scenario; WLPON represents the external source load of unstable particulate organic nitrogen (gN / day) other than the inflow / outflow scheduling scenario; WDON represents the external source load of soluble organic nitrogen (gN / day) other than the inflow / outflow scheduling scenario; WNH4 represents the external source load of ammonia nitrogen (gN / day) other than the inflow / outflow scheduling scenario; and WNO3 represents the external source load of nitrate (gN / day) other than the inflow / outflow scheduling scenario. 调度入流 WLPON 调度入流 WDON 调度入流 WNH3 调度入流 WNO3 调度入流 The inflow load (gN / day) is scheduled for inflow and outflow scenarios for inert particulate organic nitrogen, unstable particulate organic nitrogen, dissolved organic nitrogen, ammonia nitrogen, and nitrate nitrogen, respectively; WRPON 调度出流 WLPON 调度出流 WDON 调度出流 WNH3 调度出流 WNO3 调度出流 The outflow load (gN / day) is scheduled for inflow and outflow scheduling scenarios for inert particulate organic nitrogen, unstable particulate organic nitrogen, dissolved organic nitrogen, ammonia nitrogen, and nitrate nitrogen, respectively.

[0083] (4) Dissolved oxygen

[0084] The original dissolved oxygen equation can be:

[0085]

[0086] Where: AONT is the amount of dissolved oxygen consumed per unit quantity of ammonia nitrogen nitrification (4.33 g O2 / g N); AOCR is the dissolved oxygen-carbon ratio during respiration (2.67 g O2 / g C); Kr is the reoxygenation coefficient (day). -1 The reoxygenation term applies only to the surface layer; DOs is the dissolved oxygen saturation concentration (gO2 / m³). 3 SOD is the oxygen demand of sediments (gO2 / m³). 2 / day), only applicable to the bottom layer; WDO is the external dissolved oxygen load (gO2 / day) other than the inflow / outflow scheduling scenario; PNx is the preference of algae for ammonia nitrogen uptake (0 <PNx<1)。

[0087] The adjusted water quality equation can be:

[0088]

[0089] Wherein, WDO represents the external dissolved oxygen load (gO2 / day) other than the inflow / outflow scheduling scenario. 调度入流 Adjusting the inflow load (gO2 / day) for dissolved oxygen inflow / outflow scheduling scenarios; WDO 调度出流 Adjust the outflow load (gO2 / day) for the dissolved oxygen inflow and outflow scheduling scenario.

[0090] Steps 2-3: In the algae equation of the gridded hydrodynamic-water quality-aquatic ecosystem model, add inflow algae terms and outflow algae terms that quantitatively express the impact of the inflow scheduling process on the algae of the target water body.

[0091] (5) Algae

[0092] The original algal equation can be:

[0093]

[0094] Where the subscript x corresponds to different types of algae: Bx is the biomass of algae x (g·C / m³). 3 ); t represents time (days); Px represents productivity (days). -1 ); BMx is the metabolic rate (days) -1 ); PRx is the predation rate (days) -1 ); WSx is the settling velocity (m / day); WBx is the external algal load (g·C / day) other than the inflow / outflow scheduling scenario; V is the calculation unit volume of the model.

[0095] The adjusted water quality equation can be:

[0096]

[0097] Wherein, WBx represents the exogenous algal load (g·C / day) other than the inflow / outflow scheduling scenario. 调度入流 Adjusting inflow load (g·C / day) for different algal inflow / outflow scheduling scenarios; WBx 调度出流 Schedule outflow load (g·C / day) for different algae inflow and outflow scheduling scenarios.

[0098] Based on the above technical solution, the inflow scheduling and outflow scheduling processes in the inflow scheduling scenario can be parametrically expressed in IWIND-LR, which can facilitate quantitative analysis of the water quality impact of the scheduling scenario on the target water body.

[0099] Preferably, step 3 above, which involves coupling the adjusted momentum equation, water quality equation, and algae equation with the hydrodynamic-water quality-aquatic ecosystem simulation framework of the IWIND-LR model to generate an inflow-outflow scheduling model for the target water body, includes the following sub-steps, characterized in that:

[0100] Step 3-1: Rebuild the hydrodynamic simulation module of the target water body based on the adjusted momentum equation;

[0101] Step 3-2: Based on the adjusted water quality equation, rebuild the water quality simulation module of the target water body to dynamically simulate the changes in the concentrations of chemical oxygen demand, total nitrogen, total phosphorus and dissolved oxygen in the water body;

[0102] Step 3-3: Rebuild the aquatic ecosystem simulation module for the target water body based on the adjusted algae equation;

[0103] Steps 3-4: Using the newly rebuilt hydrodynamic simulation module, water quality simulation module, and aquatic ecosystem simulation module, coupled based on the simulation framework of the IWIND-LR model, and establishing the inflow and outflow scheduling model of the target water body through boundary condition settings.

[0104] Based on the above technical solution, the hydrodynamic-water quality-water ecology model can be adjusted to combine with the actual scheduling scenario and couple with the three-dimensional water quality model of the water body to construct the inflow and outflow scheduling scenario assessment simulation calculation kernel, forming the inflow and outflow scheduling scenario assessment calculation engine to analyze the impact on water quality.

[0105] Preferably, step 4 above, which involves formulating an inflow / outflow scheduling scenario scheme, includes the following sub-steps, characterized in that:

[0106] Step 4-1: In terms of time scale, the inflow and outflow scheduling scenario takes the day as the smallest time unit. Outflow and inflow can be carried out simultaneously or in stages, and can be divided into 4 categories: (1) Only inflow is scheduled, without active outflow control; (2) Active water release is carried out first, and then inflow is scheduled, with the two alternating, but water is not released when inflow occurs, and inflow is not scheduled when water is released; (3) Inflow is scheduled first, and then active water release is carried out, with the two alternating, and water is not released when inflow is scheduled, and inflow is not scheduled when water is released; (4) Both inflow and water release are scheduled, and the two can be carried out simultaneously or alternately.

[0107] Step 4-2: On a spatial scale, the inflow and outflow scheduling scenarios can be flexibly combined according to the number of inflow outlets. Specifically, they can be divided into three categories: (1) single-point inflow; (2) multi-point inflow; (3) simultaneous inflow from all inflow outlets.

[0108] Step 4-3: Based on the inflow and outflow capacity of the target water body, and combined with the decision-maker's scheduling needs, set different inflow and outflow rates to form different water scheduling scenarios.

[0109] Based on the above technical solution, considering the temporal and spatial distribution of inflow and outflow in the inflow and outflow scheduling, the inflow and outflow of water are set as variable boundary conditions. Different inflow and outflow rates are set according to the scheduling needs of decision-makers, forming different water scheduling scenarios.

[0110] Preferably, step 5 above, which involves simulating different scheduling scenarios using the inflow / outflow scheduling model and generating simulation results of the target water body's water quality response for each scheduling scenario, includes the following sub-steps, characterized in that:

[0111] Step 5-1: Transmit inflow / outflow scheduling scenario data to the inflow / outflow scheduling simulation module;

[0112] Step 5-2: Output the model simulation results, including time series plots of the target water body's water quality response, and comparison plots of two-dimensional and three-dimensional visualization layers.

[0113] Based on the above technical solutions, the time series diagram of water quality results at space stations and the comparison diagram of two-dimensional and three-dimensional visualization layers of the target water body can be used to intuitively view the impact of different inflow and outflow scheduling scenarios on water quality.

[0114] The present invention has the following advantages over the prior art:

[0115] 1. The water quality response simulation method based on the IWIND-LR model under different inflow and outflow scheduling scenarios disclosed in this application is adopted. The hydrodynamic-water quality-water ecology model of the target water body is built through the IWIND-LR model, and relevant factors of inflow and outflow scheduling scenarios are introduced into the hydrodynamic-water quality-water ecology model to construct an inflow and outflow scheduling simulation model suitable for the target water body.

[0116] 2. Through dynamic simulation of the model, taking into account the time and spatial distribution of the inflow and outflow scheduling process, the inflow and outflow of water are set as variable boundary conditions. Different inflow and outflow rates are set according to the scheduling needs of decision-makers, thereby generating simulated scheduling schemes. The impact of the scheduling schemes on water quality is verified through model simulation, thereby improving the scientificity and accuracy of water quality management.

[0117] 3. By accurately expressing the scheduling plan and quantitatively and dynamically simulating the hydrodynamic-water quality-water ecology process, the limitations of insufficient accuracy when relying solely on expert experience for water quality scheduling and management can be effectively overcome. This is of great significance for water bodies with unstable water quality that require scheduling and water replenishment to achieve water quality standards. Attached Figure Description

[0118] Figure 1 This is a flowchart illustrating a water quality response simulation method based on the IWIND-LR model under different inflow and outflow scheduling scenarios, as described in this application embodiment.

[0119] Figure 2 Flowchart for building a grid-discrete hydrodynamic-water-quality-water-ecology model for a target water body;

[0120] Figure 3 A schematic diagram illustrating the specific development process of the inflow / outflow scheduling module;

[0121] Figure 4 A schematic diagram illustrating the algorithm coupling process between the inflow / outflow scheduling simulation module and the grid-discrete hydrodynamic-water-quality-water-ecology model;

[0122] Figure 5 This is a schematic diagram of the inflow / outflow scheduling scenario for the implementation case of this application;

[0123] Figure 6 This is a schematic diagram showing the water quality simulation results under different inflow and outflow scheduling scenarios in the implementation case of this application. Detailed Implementation

[0124] To make the purpose, technical solution, and advantages of this application clearer, the following description is provided in conjunction with the appendix. Figures 1-6 The present application will be further described in detail below with reference to embodiments. It should be understood that the specific embodiments described herein are for illustrative purposes only and are not intended to limit the scope of the application.

[0125] The following will describe the specific implementation methods. Figure 1 The processing flow shown is explained in detail below:

[0126] Step 101: Construct a hydrodynamic-water quality-aquatic ecosystem model for the target lake based on the IWIND-LR model.

[0127] The IWIND-LR model is a water environment model upgraded from the internationally widely used Environmental Fluid Dynamics Code (EFDC) as its computational kernel. It includes comprehensive water quality and ecological environment indicators, pollutant migration and transformation, dissolved oxygen, organic nutrients (carbon, nitrogen, phosphorus, silicon), inorganic nutrients (carbon, nitrogen, phosphorus, silicon), algae (cyanobacteria, green algae, and diatoms), aquatic vegetation (emergent and submerged vegetation), water age, and endogenous pollutants and eutrophication dynamics. It has been widely applied in many surface waters in China, such as Taihu Lake, Chaohu Lake, Erhai Lake, Poyang Lake, Dianchi Lake, Fuxian Lake, the Yangtze River, and the Wusong River.

[0128] Based on the bottom topography and boundaries of the target water body, the grid of the study area is discretized and topographic interpolated. Water body boundary grids, inflow river flow boundaries, meteorological factor-driven boundaries, and initial water level fields are set to construct a grid-discrete hydrodynamic model. On the basis of the hydrodynamic model, a water quality process module for the target water body is coupled, considering inflow water quality conditions, pollution source concentration boundaries, and initial conditions, to achieve the construction of a grid-discrete water quality model. Based on the hydrodynamic and water quality model, a water ecological process module is coupled, considering inflow algae concentration conditions and the initial algae concentration of the water body, to achieve the establishment of a grid-discrete water ecological model. A standardized data framework format for the hydrodynamic, water quality, and water ecological models of the target water body is established to achieve automatic data transfer between the three sub-models. The process of building the grid-discrete hydrodynamic-water quality-water ecological model is described in [link to documentation]. Figure 2 As shown.

[0129] Step 102: Based on the hydrodynamic and water quality mechanisms corresponding to the inflow and outflow scheduling, adjust the momentum equation, water quality equation, and algae equation in the above model. Specifically, add inflow and outflow terms to the momentum equation, inflow and outflow water quality terms to the water quality equation, and inflow and outflow algae terms to the algae equation, quantitatively expressing the impact of inflow and outflow on the hydrodynamic and water quality processes of the target water body.

[0130] The adjustments to the momentum equation, water quality equation, and algae equation in step 102 can be as follows: In the momentum equation of the hydrodynamic-water quality-aquatic ecosystem model, add inflow and outflow terms that quantitatively express the impact of the inflow scheduling process on the hydrodynamics of the target water body; in the water quality equation of the hydrodynamic-water quality-aquatic ecosystem model, add inflow and outflow water quality terms that quantitatively express the impact of the inflow scheduling process on the water quality of the target water body; in the algae equation of the hydrodynamic-water quality-aquatic ecosystem model, add inflow and outflow algae terms that quantitatively express the impact of the inflow scheduling process on the algae of the target water body.

[0131] The specific considerations for developing the inflow / outflow scheduling module are as follows: Figure 3 As shown.

[0132] Step 103: Couple the adjusted momentum equation, water quality equation, and algae equation with the hydrodynamic-water quality-aquatic ecosystem simulation framework of the IWIND-LR model to generate the inflow and outflow scheduling model for the target water body:

[0133] Using the newly rebuilt hydrodynamic simulation module, water quality simulation module, and aquatic ecosystem simulation module, coupled with the simulation framework based on the IWIND-LR model, and by setting boundary conditions, an inflow and outflow scheduling model for the target water body is established.

[0134] The inflow and outflow boundary conditions, based on the original grid-discrete hydrodynamic-water quality model boundary conditions (including inflow and outflow water bodies), include new inflow and outflow boundaries from the inflow and outflow scheduling module. These new boundaries are primarily obtained through the inflow and outflow scheduling scenario scheme defined in step 104, and then converted into inflow and outflow flow rate time series, inflow and outflow water quality concentration time series, and inflow and outflow algae concentration time series formats that can be directly used by the inflow and outflow scheduling simulation module. The algorithm coupling process between the inflow and outflow scheduling simulation module and the grid-discrete hydrodynamic-water quality-aquatic ecosystem model is as follows: Figure 4 As shown.

[0135] Step 104, Develop inflow / outflow scheduling scenario plan:

[0136] By combining the inflow and outflow capacity of the target water body and taking into account the temporal and spatial distribution of inflow and outflow, the inflow and outflow of the water body are set as variable boundary conditions. Different inflow and outflow rates are set according to the decision-maker's scheduling needs to form different water volume scheduling scenarios, which serve as input data for the inflow and outflow scheduling simulation module.

[0137] In the implementation case, in terms of time, both inflow and outflow were scheduled simultaneously, and the scheduling was carried out from March to May; in terms of space, the two scheduling outlets were used for single-point inflow and simultaneous inflow from both inflow outlets, while the only outflow outlet was used for outflow; in terms of water volume, a constant inflow and outflow scheduling of 1, 2, 3, and 4 cubic meters per second was adopted.

[0138] Eight scheduling modes were set up: Mode 1: From March to May, Inlet #1 would receive water at a rate of 1 cubic meter per second, with a constant outflow of 1 cubic meter per second at the outlet; Mode 2: From March to May, Inlet #1 would receive water at a rate of 2 cubic meters per second, with a constant outflow of 2 cubic meters per second at the outlet; Mode 3: From March to May, Inlet #2 would receive water at a rate of 1 cubic meter per second, with a constant outflow of 1 cubic meter per second at the outlet; Mode 4: From March to May, Inlet #2 would receive water at a rate of 2 cubic meters per second, with a constant outflow of 2 cubic meters per second at the outlet; Mode 5: From March to May, Inlet #1... Simultaneous inflow at inlet #2, with inflow rates of 1 and 1 cubic meters per second respectively, and a constant outflow rate of 2 cubic meters per second at the outlet; sixth, simultaneous inflow at inlet #1 and #2 from March to May, with inflow rates of 1 and 2 cubic meters per second respectively, and a constant outflow rate of 3 cubic meters per second at the outlet; seventh, simultaneous inflow at inlet #1 and #2 from March to May, with inflow rates of 2 and 1 cubic meters per second respectively, and a constant outflow rate of 3 cubic meters per second at the outlet; eighth, simultaneous inflow at inlet #1 and #2 from March to May, with inflow rates of 2 and 2 cubic meters per second respectively, and a constant outflow rate of 4 cubic meters per second at the outlet.

[0139] For detailed procedures, please see [link / document / documentation].Figure 5 As shown.

[0140] Step 105: Using the inflow / outflow scheduling model, simulate the different scheduling scenarios and generate simulation results of the target water body's water quality response for each scheduling scenario:

[0141] After completing the research and development and construction of the target water body inflow and outflow scheduling simulation module, the model was run using the same time step as the three-dimensional hydrodynamic and water quality model of the target water body. Model stability was achieved through continuous adjustments to model parameters and initial states. Water quality simulation results for different inflow and outflow scheduling scenarios were generated using the 104 steps generated. The simulation results for four scheduling scenarios are shown below. Figure 6 As shown in the figure, the different shapes of the lines represent the simulation results of different inflow and outflow scheduling schemes. The figure effectively evaluates the impact of inflow and outflow scheduling on water quality.

[0142] In summary, this invention details a water quality simulation method for inflow and outflow scheduling scenarios based on the IWIND-LR model. From standardized data conversion channels, efficient model construction, complete inflow and outflow scheduling module development, and flexible result output, it directly establishes a quantitative relationship between different scheduling scenarios and water quality, forming a water quality simulation method for inflow and outflow scheduling scenarios based on the IWIND-LR model. This method can provide important theoretical extensions and corresponding decision support for the daily scheduling and management of water bodies such as lakes and reservoirs.

[0143] This invention is not limited to the content described in the above examples, but is defined by the scope of the claims. Furthermore, any modifications, additions, or equivalent substitutions made by those skilled in the art based on these examples are within the scope of protection of the claims. It should be understood that the specific embodiments described herein are for illustrative purposes only and are not intended to limit this application.

Claims

1. A method for quantitatively simulating water quality response of a target water body under different scenarios of inflow and outflow scheduling based on an IWIND-LR model, characterized in that, The method comprises the following steps: Step 1, building a hydrodynamic-water quality-ecological model of a target water body based on an IWIND-LR model; Step 2, adjusting a momentum equation, a water quality equation and an algae equation in the hydrodynamic-water quality-ecological model based on a water dynamic and water quality mechanism process corresponding to flow-in and flow-out scheduling, specifically including adding a flow-in flow rate term and a flow-out flow rate term in the momentum equation, adding a flow-in water quality term and a flow-out water quality term in the water quality equation, and adding a flow-in algae term and a flow-out algae term in the algae equation, so as to quantitatively express the influence of flow-in and flow-out on the hydrodynamic and water quality process of the target water body; Step 3, coupling the adjusted momentum equation, water quality equation and algae equation with a hydrodynamic-water quality-ecological simulation framework of the IWIND-LR model to generate a flow-in and flow-out scheduling model of the target water body, specifically including rebuilding a hydrodynamic simulation module based on the adjusted momentum equation, rebuilding a water quality simulation module based on the adjusted water quality equation, and rebuilding an ecological simulation module based on the adjusted algae equation, and coupling based on the simulation framework of the IWIND-LR model to generate the flow-in and flow-out scheduling model of the target water body through boundary condition setting; Step 4, formulating a flow-in and flow-out scheduling scenario scheme, specifically including setting the water body inflow and outflow as variable boundary conditions by combining the flow-in and flow-out capacity of the target water body, taking into account the time and spatial distribution of the flow-in and flow-out scheduling inflow and outflow, setting different inflow and outflow according to the scheduling requirements of the decision maker to form different water quantity scheduling scenarios as input data of the flow-in and flow-out scheduling model; Step 5, simulating the different flow-in and flow-out scheduling scenarios through the flow-in and flow-out scheduling model to generate a simulation result of the water quality response of the target water body corresponding to each scheduling scenario.

2. The method of claim 1, wherein, The adjustment of the momentum equation, water quality equation and algae equation in the hydrodynamic-water quality-ecological model based on the water dynamic and water quality mechanism process corresponding to the flow-in and flow-out scheduling comprises: In the momentum equation in the hydrodynamic-water quality-ecological model, the flow-in flow rate term and the flow-out flow rate term are added to quantitatively express the influence of the flow-in and flow-out scheduling process on the hydrodynamic process of the target water body; In the water quality equation in the hydrodynamic-water quality-ecological model, the flow-in water quality term and the flow-out water quality term are added to quantitatively express the influence of the flow-in and flow-out scheduling process on the water quality of the target water body; In the algae equation in the hydrodynamic-water quality-ecological model, the flow-in algae term and the flow-out algae term are added to quantitatively express the influence of the flow-in and flow-out scheduling process on the algae of the target water body.

3. The method of claim 1, wherein, The coupling of the adjusted momentum equation, water quality equation and algae equation with the hydrodynamic-water quality-ecological simulation framework of the IWIND-LR model to generate the flow-in and flow-out scheduling model of the target water body comprises: Rebuilding a hydrodynamic simulation module of the target water body based on the adjusted momentum equation; Rebuilding a water quality simulation module of the target water body based on the adjusted water quality equation to dynamically simulate the concentration changes of chemical oxygen demand, total nitrogen, total phosphorus and dissolved oxygen in the water body; Rebuilding an ecological simulation module of the target water body based on the adjusted algae equation; The water inflow and outflow scheduling model of the target water body is established based on the simulation framework of the IWIND-LR model by coupling the reestablished water power simulation module, water quality simulation module and water ecological simulation module and setting boundary conditions.

4. The method of claim 1, wherein, The scheduling scenario scheme of the inflow and outflow is formulated, and the scheme includes: In the time scale, the inflow and outflow scheduling scenario takes day as the minimum time unit, and the outflow and inflow can be performed synchronously or in steps; In the spatial scale, the inflow and outflow scheduling scenario can be flexibly combined according to the number of inflow ports, and the single-point inflow, multi-point inflow and simultaneous inflow of all inflow ports can be adopted; Based on the inflow and outflow capacity of the target water body, different inflow and outflow are set in combination with the scheduling demand of the decision maker to form different water quantity scheduling scenarios.

5. The method of claim 1, wherein, The scheduling scenarios are simulated through the inflow and outflow scheduling model to generate water quality response simulation results of the target water body under each scheduling scenario; the water quality response simulation results include a time series graph of the water quality response of the target water body, and a two-dimensional and three-dimensional visualization layer comparison graph.

Images