A method for predicting and analyzing ecological water supplement quantity of inland lakes

Abstract

The application discloses a kind of inland lake ecological water supplement amount prediction analysis method, comprising: obtaining the underwater topographic map of target lake, constructs the plane two-dimensional water flow movement basic equation of lake wind lake current and the plane two-dimensional water quality equation describing salt transport in lake conversion;Obtain the ecological water supplement flow and salt concentration of target lake, fit the plane two-dimensional water quality equation;To the finite difference discrete of continuity equation, momentum equation, the plane two-dimensional water flow movement basic equation is converted into linear algebra equation group;Linear algebra equation group is solved by ADI row-by-row method in time period, direction, and the water level, water depth and salt concentration of output prediction are exported;The water depth of each grid ecological water supplement is calculated, and then the comprehensive water depth of target lake ecological water supplement in future time is calculated, and the ecological water supplement amount of target lake prediction is obtained.

Description

Technical Field

[0001] This invention relates to the field of lake ecology research, specifically to a method for predicting and analyzing the ecological water replenishment of inland lakes. Background Technology

[0002] Inland lakes, as important ecological environments, often serve as ideal breeding grounds for birds. For example, the Relict Gull is a typical bird species in desert-semi-desert ecosystems, using wetland habitats as its habitat. According to long-term observations and research in the Relict Gull Nature Reserve, the Relict Gulls stay at Tao-Ahaizi for about 5 months, arriving in early April each year and flying away by the end of September. The Relict Gulls only use the island in the middle of the lake as their breeding ground; they feel insecure outside of the island. All adult Relict Gulls participate in breeding. The nests of the Relict Gulls are simply constructed using pebbles from the island in the middle of the lake. During the breeding season, each nest contains 2-3 eggs, and the incubation of the eggs is completed by both parents.

[0003] As a bird species that exclusively breeds on lake islands in arid to semi-arid environments with considerable instability, the Relict Gull's population fluctuations are closely related to ecological and environmental influences. Climate and lake water level changes are the main factors affecting its breeding. The Relict Gull's breeding habitat is significantly affected by the water surface area. If the lake area is too small, the lake island will be connected to the outside land, and various small mammals will threaten the incubation of the eggs; if the lake area is too high, the lake island will be submerged, making it impossible for the bird to lay eggs and breed.

[0004] It is evident that maintaining a reasonable water level in inland lakes is crucial for protecting the surrounding ecological environment. However, inland areas suffer from water scarcity and limited available ecological water transfer volumes. Therefore, rationally planning ecological water transfer is essential for maintaining inland lake water levels. Consequently, there is an urgent need to develop a predictive analysis method for ecological water replenishment to inland lakes. Summary of the Invention

[0005] To address the aforementioned shortcomings of existing technologies, this invention provides a method for predicting and analyzing the ecological water replenishment of inland lakes. By performing hydrodynamic analysis and fitting on the target lake, the method enables reasonable calculation and prediction of the future ecological water replenishment of the target lake.

[0006] To achieve the above-mentioned objectives, the technical solution adopted by this invention is as follows: A method for predicting and analyzing the ecological water replenishment of inland lakes is provided, which includes: S1: Obtain the underwater topographic map of the target lake and divide it into... N A grid is used to construct the basic equations of two-dimensional planar water flow motion of lake wind-generated lake currents and the two-dimensional planar water quality equations describing the transport and transformation of salt in lakes; S2: Collect the salinity concentration at several points in the target lake, and assign a value to the salinity concentration of the grid according to the salinity concentration at each point; obtain the ecological water replenishment flow and salinity concentration of the water replenishment of the target lake, fit the two-dimensional planar water quality equation, and obtain the converged two-dimensional planar water quality equation for each grid. S3: Discretize the fundamental equations of two-dimensional planar water flow motion into each grid and set the time step for the fitting process. By performing finite difference discretization on the continuity equation and momentum equation, the basic equations of two-dimensional planar water flow motion are transformed into a system of linear algebraic equations. S4: Adjust the time step Split into two half-steps The Analytic Difference (ADI) method is used to implicitly solve the linear algebraic equations in different time periods and directions, outputting the predicted water level and depth. These values ​​are then input into a convergent two-dimensional planar water quality equation to obtain the predicted future time. Salt concentration at the location of the grid ; S5: Based on the ideal salinity concentration of the target lake Calculate the predicted ecological replenishment water depth for each grid cell, and then calculate the future time horizon of the target lake. The comprehensive water depth for ecological water replenishment is used to obtain the predicted ecological water replenishment volume for the target lake.

[0007] Further, step S1 includes: S11: Obtain underwater topographic maps of the target lake and construct a two-dimensional grid system on the underwater topographic maps, dividing the plane of each underwater topographic map into... N There are 1 grid, and the center coordinates of each grid are 1. ; S12: Construct the basic two-dimensional planar water flow motion equations for lake wind-driven currents and the two-dimensional planar water quality equations describing the transport and transformation of salt in lakes; The fundamental equations of two-dimensional planar water flow motion include: Continuity equation: ; x The momentum equation for direction: ; y The momentum equation for direction: ; in, For grid Location in time t water level, They are respectively x direction, y Flux in the direction, and , They are respectively x direction, y Water flow velocity in the direction, For instantaneous water depth, For the water depth that changes over time, g It is the acceleration due to gravity. C For the Xie Cai coefficient, The Coriolis force coefficient for Earth's rotation. For wind resistance, They are respectively x direction, y directional wind speed component, For the density of water, Atmospheric pressure, E This is the eddy viscosity coefficient; The two-dimensional planar water quality equation describing the transport and transformation of salt in lakes is as follows: ; in, c This refers to the salt concentration in the lake. Salt content x direction, y Diffusion coefficient in the direction, k The attenuation coefficient is... S For source and sink items, , For source and sink flow, This represents the salt concentration in the source and sink terms.

[0008] Further, step S2 includes: S21: Collect salinity concentrations at several points in the target lake, and assign salinity concentration values ​​to the grid based on the salinity concentration at each point. Specifically, this includes: If there are points within the grid, the salinity concentration collected at those points is assigned as the salinity concentration of the grid, and the grid with points is used as the reference grid. If there are no points within the grid, it is used as the interpolation grid, and the salinity concentration of the interpolation grid is linearly interpolated using the salinity concentration of the reference grid. ; in, The salt concentration is obtained by linear interpolation of the interpolation grid. For interpolation grid in x The salt concentration is obtained by linear interpolation in the direction. For interpolation grid in y The salt concentration is obtained by linear interpolation of the direction. S22: Obtain the ecological water replenishment flow and salinity concentration of the target lake, input them into the two-dimensional planar water quality equation, and obtain the water depth data of the corresponding grid positions at each historical time based on the underwater topographic map. Fit the two-dimensional planar water quality equation to calculate the salinity. x direction, y Diffusion coefficient in the direction and attenuation coefficient k This yields the convergent two-dimensional planar water quality equation for each grid cell; Further, step S4 includes: S41: Adjust the time step Split into two half-steps The ADI line-by-line method is used to implicitly solve the linear algebraic equations in different time periods and directions. x Direction is implicitly calculated, for x The unknowns of the linear algebraic equations in the direction are solved implicitly. y The linear algebraic equations of direction and continuity are solved explicitly to obtain the half-step size. Corresponding flux median ; S42: Based on flux median along y Direction is implicitly calculated, for y The unknowns of the linear algebraic equations in the direction are solved implicitly. x The linear algebraic equations for direction and continuity are solved explicitly to obtain the final flux value with an integer step length. ; S43: Final flux value Input the linear algebraic equation of the continuity equation and solve it explicitly to obtain the water level of integer step length. and water depth ; S44: Utilizing water levels in whole steps and water depth Verify the set time step Whether it converges; If satisfied and , Let the water level residual and water depth residual be used to determine the time step. If convergence is achieved, proceed to step S45; otherwise, determine the time step. If the circuit does not converge, return to step S3 to adjust the time step. And execute steps S3-S44 until the time step is reached. Once convergence is achieved, proceed to step S45. S45: Adjust the time step The water level calculated after convergence , water depth Final flux value Inputting a convergent two-dimensional planar water quality equation, the predicted future time is calculated. Grid Salt concentration at location .

[0009] Further, step S5 includes: S51: Based on the ideal salinity concentration of the target lake Calculate the predicted difference in salt concentration Differences in salt concentration Input a convergent two-dimensional planar water quality equation to calculate the predicted future time. Grid Location of ecological water replenishment water depth ; S52: Based on the ecological water replenishment depth corresponding to each grid. Calculate the target lake in the future time Comprehensive water depth for ecological water replenishment Utilizing comprehensive water depth Calculate the predicted future time based on the area of ​​the target lake. The amount of ecological water replenishment for the target lake.

[0010] .

[0011] The beneficial effects of this invention are as follows: This invention achieves the analysis of hydrodynamics and dynamic changes in salinity concentration fields within inland lakes by constructing a two-dimensional planar equation for the fundamental motion of lake currents and a two-dimensional planar water quality equation describing the transport and transformation of salt in lakes. Based on the coupled calculation of salinity concentration field and hydrodynamics, it correlates water level changes with salinity concentration distribution, and uses the salinity concentration change process to predict future water level changes, thereby obtaining the required ecological water replenishment. This provides effective data support for ensuring the ideal water level of target lakes and protecting the ecological environment, and provides accurate data support for the scheduling and planning of future ecological water replenishment. Attached Figure Description

[0012] Figure 1 This is an underwater topographic map of the Tao-Ahai lake, as measured in actual measurements. Detailed Implementation

[0013] The specific embodiments of the present invention are described below to enable those skilled in the art to understand the present invention. However, it should be understood that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, various changes are obvious as long as they are within the spirit and scope of the present invention as defined and determined by the appended claims. All inventions utilizing the concept of the present invention are protected.

[0014] A method for predicting and analyzing the ecological water replenishment of inland lakes includes the following steps: S1: Obtain the underwater topographic map of the target lake and divide it into... N Using a grid, we construct a two-dimensional planar equation for the fundamental motion of lake currents and a two-dimensional planar water quality equation describing the transport and transformation of salt in the lake. Step S1 specifically includes: S11: Obtain underwater topographic maps of the target lake and construct a two-dimensional grid system on the underwater topographic maps, dividing the plane of each underwater topographic map into... N There are 1 grid, and the center coordinates of each grid are 1. ; like Figure 1 As shown, this embodiment uses the underwater topographic map of Tao-Ahaizi Lake measured in actual conditions as an example. Figure 1 In this context, "Bathymetry" represents water depth, "grid spacing 40 meter" indicates a grid size of 40 meters, the measurement grid system uses the 2000 National Geodetic Grid System with a Gauss-Kruger projection 3º zone, and the elevation system uses the 1985 National Elevation System. The total east-west length of Tao-Ahaizi is approximately 6.52 km, and the total north-south width is approximately 3.56 km. Considering the complex geometric boundaries of the lake, a rectangular grid is used for mesh generation, and a relatively dense computational grid is adopted to ensure computational accuracy. The longitudinal and transverse computational grid lengths of Tao-Ahaizi are determined to be Δ. x =40m、△ y =40m, the total number of grids involved in the calculation is 89×163=14507.

[0015] S12: Construct the basic two-dimensional planar water flow motion equations for lake wind-driven currents and the two-dimensional planar water quality equations describing the transport and transformation of salt in lakes; The fundamental equations of two-dimensional planar water flow motion include: Continuity equation: ; x The momentum equation for direction: ; y The momentum equation for direction: ; in, For grid Location in time t water level, They are respectively x direction, y Flux in the direction, and , They are respectively x direction, y Water flow velocity in the direction, For instantaneous water depth, For the water depth that changes over time, g It is the acceleration due to gravity. C For the Xie Cai coefficient, The Coriolis force coefficient for Earth's rotation. For wind resistance, They are respectively x direction, y directional wind speed component, For the density of water, Atmospheric pressure, E This is the eddy viscosity coefficient; The two-dimensional planar water quality equation describing the transport and transformation of salt in lakes is as follows: ; in, c This refers to the salt concentration in the lake. Salt content x direction, y Diffusion coefficient in the direction, k The attenuation coefficient is... S For source and sink items, , For source and sink flow, The salt concentration in the source and sink terms; in this embodiment, salt is an inert substance, neither decreasing nor increasing, but only migrating and diffusing with the water flow, so the decrease coefficient is taken. k =0.

[0016] S2: Collect salinity concentrations at several points in the target lake, assign salinity values ​​to the grid based on the salinity concentration at each point, obtain the ecological water replenishment flow and salinity concentration of the target lake, fit the two-dimensional planar water quality equation, and obtain the converged two-dimensional planar water quality equation for each grid. Step S2 specifically includes: S21: Collect salinity concentrations at several points in the target lake, and assign salinity concentration values ​​to the grid based on the salinity concentration at each point. Specifically, this includes: If there are points within the grid, the salinity concentration collected at those points is assigned as the salinity concentration of the grid, and the grid with points is used as the reference grid. If there are no points within the grid, it is used as the interpolation grid, and the salinity concentration of the interpolation grid is linearly interpolated using the salinity concentration of the reference grid. ; in, The salt concentration is obtained by linear interpolation of the interpolation grid. For interpolation grid in x The salt concentration is obtained by linear interpolation in the direction. For interpolation grid in y The salt concentration is obtained by linear interpolation of the direction.

[0017] In this embodiment, the salt concentration of the interpolation grid is at x direction and y When performing linear interpolation in the direction, the salt concentration corresponding to the nearest reference grid at both ends of the interpolation grid is taken as a reference. The difference between the known salt concentrations of the two reference grids is divided by the number of interpolation grids in the middle to obtain the slope of the salt concentration change. Then, taking one of the known salt concentrations as a reference, the salt concentration of each interpolation grid in the middle is calculated to achieve linear interpolation of the salt concentration of the interpolation grid.

[0018] S22: Obtain the ecological replenishment flow and salinity concentration of the target lake. In this embodiment, the inflow of water into Tao-Ahaizi Lake consists of runoff (combining surface and groundwater runoff), precipitation, and mine water replenishment. The annual average values ​​of runoff (combining surface and groundwater runoff) and precipitation are taken, with mine water replenishment being the primary consideration. Input the data into the two-dimensional planar water quality equation, and obtain the water depth data corresponding to the grid positions at each historical time based on the underwater topographic map. Fit the two-dimensional planar water quality equation to calculate the salinity concentration. x direction, y Diffusion coefficient in the direction and attenuation coefficient k The convergent two-dimensional planar water quality equations for each grid are obtained. The ecological replenishment flow rate is determined based on the ecological replenishment capacity that can be provided to the target lake. Generally, the ecological replenishment flow rate corresponding to the maximum ecological replenishment capacity is taken. Both the ecological replenishment flow rate and the salinity concentration are fixed values. The ecological replenishment flow rate is used as the source and sink flow rate.

[0019] S3: Discretize the fundamental equations of two-dimensional planar water flow motion into each grid and set the time step for the fitting process. The continuity equation and momentum equation are discretized by finite difference, including the use of second-order central difference for linear terms and first-order upwind scheme for nonlinear terms, thus transforming the basic equations of two-dimensional planar water flow motion into a system of linear algebraic equations. In this embodiment, the linear terms in the basic equations of planar two-dimensional water flow motion include gravity terms, Coriolis force terms, atmospheric pressure terms, and eddy viscosity terms, while the nonlinear terms include inertial terms, wind resistance terms, and water density terms. The time partial derivatives are expressed using a forward difference scheme.

[0020] The time-difference and space-difference methods for the continuity equation are as follows: ; ; The linear algebraic equation obtained after the continuity equation is discretized by finite difference is: ; in, For grid Location in time t +1 water level, For time t +1 water depth, They are respectively x Directional grid Location flux, y Directional grid Flux at a location.

[0021] In this embodiment, x Inertial term of direction For example, the spatial difference method for performing a first-order upwind scheme is as follows: ; in, For grid The instantaneous water depth at the location; by x Water level gradient term in the direction For example, the method for performing central difference is as follows: ; in, For grid Location in time t water level, For grid Location in time t The water level.

[0022] The finite difference discretization method for the remaining linear and nonlinear terms is similar. x The momentum equations of direction are discretized in terms of time. t +1 x Directional flux Linear algebraic equations with unknowns y The momentum equations of direction are discretized in terms of time. t +1 y Directional flux A system of linear algebraic equations is formed by combining linear algebraic equations with unknowns and transforming them into continuous equations.

[0023] S4: Adjust the time step Split into two half-steps The Analytic Difference (ADI) method is used to implicitly solve the linear algebraic equations in different time periods and directions, outputting the predicted water level, water depth, and flux. These values ​​are then input into a convergent two-dimensional planar water quality equation to obtain the predicted future time. Salt concentration at the location of the grid Step S4 specifically includes: S41: Adjust the time step Split into two half-steps The ADI line-by-line method is used to implicitly solve the linear algebraic equations in different time periods and directions. x Direction is implicitly calculated, for x The unknowns of the linear algebraic equations in the direction are solved implicitly. y The linear algebraic equations of direction and continuity are solved explicitly to obtain the half-step size. Corresponding flux median ; S42: Based on flux median along y Direction is implicitly calculated, for y The unknowns of the linear algebraic equations in the direction are solved implicitly. x The linear algebraic equations for direction and continuity are solved explicitly to obtain the final flux value with an integer step length. ; S43: Final flux value Input the linear algebraic equation of the continuity equation and solve it explicitly to obtain the water level of integer step length. and water depth ; S44: Utilizing water levels in whole steps and water depth Verify the set time step Whether it converges; If satisfied and , Let the water level residual and water depth residual be used to determine the time step. If convergence is achieved, proceed to step S45; otherwise, determine the time step. If the circuit does not converge, return to step S3 to adjust the time step. And execute steps S3-S44 until the time step is reached. Once convergence is achieved, proceed to step S45. S45: Adjust the time step The water level calculated after convergence , water depth Final flux value Inputting a convergent two-dimensional planar water quality equation, the predicted future time is calculated. Grid Salt concentration at location .

[0024] S5: Based on the ideal salinity concentration of the target lake The difference in salinity concentration is calculated and input into a convergent two-dimensional planar water quality equation. The predicted ecological replenishment depth for each grid cell is then calculated, and the future time horizon of the target lake is calculated. The comprehensive water depth for ecological water replenishment is used to obtain the predicted ecological water replenishment volume for the target lake. Step S5 specifically includes: S51: Based on the ideal salinity concentration of the target lake Calculate the predicted difference in salt concentration Differences in salt concentration Input a convergent two-dimensional planar water quality equation to calculate the predicted future time. Grid Location of ecological water replenishment water depth ; S52: Based on the ecological water replenishment depth corresponding to each grid. Calculate the target lake in the future time Comprehensive water depth for ecological water replenishment Utilizing comprehensive water depth Calculate the predicted future time based on the area of ​​the target lake. The amount of ecological water replenishment for the target lake.

[0025] .

[0026] In this embodiment, during the long-term cumulative ecological water replenishment process, salt will accumulate. Even when the area of ​​the target lake reaches the maximum constrained area, the lake salinity will still be higher than the ideal salt concentration. In such cases, timely warnings and desalination measures are implemented to reduce lake salinity.

[0027] This invention analyzes the hydrodynamics and dynamic changes of salinity concentration fields in inland lakes by constructing a two-dimensional planar equation for the fundamental motion of lake currents and a two-dimensional planar water quality equation describing the transport and transformation of salt in the lake. Based on the coupled calculation of salinity concentration field and hydrodynamics, it correlates water level changes with salinity concentration distribution, predicts future water level and salinity changes, and thus obtains the required ecological water replenishment. This provides effective data support for ensuring ideal water levels in target lakes and protecting the ecological environment, and provides accurate data support for the scheduling and planning of future ecological water replenishment.

Claims

1. A method for predicting and analyzing the ecological water replenishment of inland lakes, characterized in that, include: S1: Obtain the underwater topographic map of the target lake and divide it into... N A grid is used to construct the basic equations of two-dimensional planar water flow motion of lake wind-generated lake currents and the two-dimensional planar water quality equations describing the transport and transformation of salt in lakes; S2: Collect the salinity concentration at several points in the target lake, and assign a value to the salinity concentration of the grid according to the salinity concentration at each point; obtain the ecological water replenishment flow and salinity concentration of the water replenishment of the target lake, fit the two-dimensional planar water quality equation, and obtain the converged two-dimensional planar water quality equation for each grid. S3: Discretize the fundamental equations of two-dimensional planar water flow motion into each grid and set the time step for the fitting process. By performing finite difference discretization on the continuity equation and momentum equation, the basic equations of two-dimensional planar water flow motion are transformed into a system of linear algebraic equations. S4: Adjust the time step Split into two half-steps The Analytic Difference (ADI) method is used to implicitly solve the linear algebraic equations in different time periods and directions, outputting the predicted water level and depth. These values ​​are then input into a convergent two-dimensional planar water quality equation to obtain the predicted future time. Salt concentration at the location of the grid ; S5: Based on the ideal salinity concentration of the target lake Calculate the predicted ecological replenishment water depth for each grid cell, and then calculate the future time horizon of the target lake. The comprehensive water depth for ecological water replenishment is used to obtain the predicted ecological water replenishment volume for the target lake.

2. The method for predicting and analyzing the ecological water replenishment of inland lakes according to claim 1, characterized in that, Step S1 includes: S11: Obtain underwater topographic maps of the target lake and construct a two-dimensional grid system on the underwater topographic maps, dividing the plane of each underwater topographic map into... N There are 1 grid, and the center coordinates of each grid are 1. ; S12: Construct the basic two-dimensional planar water flow motion equations for lake wind-driven currents and the two-dimensional planar water quality equations describing the transport and transformation of salt in lakes; The fundamental equations of two-dimensional planar water flow motion include: Continuity equation: ； x The momentum equation for direction: ； y The momentum equation for direction: ； in, For grid Location in time t water level, They are respectively x direction, y Flux in the direction, and , They are respectively x direction, y Water flow velocity in the direction, For instantaneous water depth, For the water depth that changes over time, g It is the acceleration due to gravity. C For the Xie Cai coefficient, The Coriolis force coefficient for Earth's rotation. For wind resistance, They are respectively x direction, y directional wind speed component, For the density of water, Atmospheric pressure, E This is the eddy viscosity coefficient; The two-dimensional planar water quality equation describing the transport and transformation of salt in lakes is as follows: ； in, c This refers to the salt concentration in the lake. Salt content x direction, y Diffusion coefficient in the direction, k The attenuation coefficient is... S For source and sink items, , For source and sink flow, This represents the salt concentration in the source and sink terms.

3. The method for predicting and analyzing the ecological water replenishment of inland lakes according to claim 2, characterized in that, Step S2 includes: S21: Collect salinity concentrations at several points in the target lake, and assign salinity concentration values ​​to the grid based on the salinity concentration at each point. Specifically, this includes: If there are points within the grid, the salinity concentration collected at those points is assigned as the salinity concentration of the grid, and the grid with points is used as the reference grid. If there are no points within the grid, it is used as the interpolation grid, and the salinity concentration of the interpolation grid is linearly interpolated using the salinity concentration of the reference grid. ； in, The salt concentration is obtained by linear interpolation of the interpolation grid. For interpolation grid in x The salt concentration is obtained by linear interpolation in the direction. For interpolation grid in y The salt concentration is obtained by linear interpolation of the direction. S22: Obtain the ecological water replenishment flow and salinity concentration of the target lake, input them into the two-dimensional planar water quality equation, and obtain the water depth data of the corresponding grid positions at each historical time based on the underwater topographic map. Fit the two-dimensional planar water quality equation to calculate the salinity. x direction, y Diffusion coefficient in the direction and attenuation coefficient k This yields the convergent two-dimensional planar water quality equation for each grid cell.

4. The method for predicting and analyzing the ecological water replenishment of inland lakes according to claim 3, characterized in that, Step S4 includes: S41: Adjust the time step Split into two half-steps The ADI line-by-line method is used to implicitly solve the linear algebraic equations in different time periods and directions. x Direction is implicitly calculated, for x The unknowns of the linear algebraic equations in the direction are solved implicitly. y The linear algebraic equations of direction and continuity are solved explicitly to obtain the half-step size. Corresponding flux median ; S42: Based on flux median along y Direction is implicitly calculated, for y The unknowns of the linear algebraic equations in the direction are solved implicitly. x The linear algebraic equations for direction and continuity are solved explicitly to obtain the final flux value with an integer step length. ; S43: Final flux value Input the linear algebraic equation of the continuity equation and solve it explicitly to obtain the water level of integer step length. and water depth ; S44: Utilizing water levels in whole steps and water depth Verify the set time step Whether it converges; If satisfied and , Let the water level residual and water depth residual be used to determine the time step. If convergence is achieved, proceed to step S45; otherwise, determine the time step. If the circuit does not converge, return to step S3 to adjust the time step. And execute steps S3-S44 until the time step is reached. Once convergence is achieved, proceed to step S45. S45: Adjust the time step The water level calculated after convergence , water depth Final flux value Inputting a convergent two-dimensional planar water quality equation, the predicted future time is calculated. Grid Salt concentration at location .

5. The method for predicting and analyzing the ecological water replenishment of inland lakes according to claim 4, characterized in that, Step S5 includes: S51: Based on the ideal salinity concentration of the target lake Calculate the predicted difference in salt concentration Differences in salt concentration Input a convergent two-dimensional planar water quality equation to calculate the predicted future time. Grid Location of ecological water replenishment water depth ; S52: Based on the ecological water replenishment depth corresponding to each grid. Calculate the target lake in the future time Comprehensive water depth for ecological water replenishment Utilizing comprehensive water depth Calculate the predicted future time based on the area of ​​the target lake. The ecological water replenishment volume of the target lake; 。

Images