An estuary and coast gate-controlled water system water body exchange simulation method

Abstract

This invention discloses a method for simulating water exchange in estuarine and coastal gate-controlled water systems. This method decomposes the water exchange period into multiple time periods, simulating them sequentially by time. By extracting the model's execution results from the previous time period and analyzing the internal and external water level processes and the relationship between safe high and low water levels for water exchange, the opening and closing of the water system gates are determined, enabling simulation of water exchange during complex scheduling of multi-gate water systems. This method achieves automated simulation of water exchange during normal scheduling of estuarine and coastal gate-controlled water systems through repetitive programmable operations such as data extraction, intersection point finding, and modification of model configuration files. It avoids human intervention and interference during the simulation process, offering advantages such as high efficiency and wide applicability, and has broad application prospects in coastal city construction, tidal flat and island development and utilization.

Description

Technical Field

[0001] This invention relates to a water exchange simulation method for estuarine and coastal gate control systems, specifically a water exchange simulation method under normal operation of estuarine and coastal gate control systems under given high and low water exchange levels, belonging to the field of coastal dynamics technology. Background Technology

[0002] In the development of many new coastal cities, reclaimed tidal flats, wetland parks, and (artificial) islands along river estuaries, a series of enclosed waterways and waterways are often formed within them. To prevent the deterioration of water quality and the aquatic environment in these enclosed water systems, it is necessary to connect them to open waterways for water exchange. Affected by factors such as astronomical tides and storm surges, the tidal levels in estuarine and coastal waters are constantly changing. For water safety considerations such as preventing flooding due to tides, these water systems often employ gate control to utilize tidal dynamics to exchange water within a certain safe water exchange level range. When the water level fluctuates outside this safe range, the gates are closed to stop the water exchange. This design involves the frequent opening and closing of the water system gates during the water exchange process, and the timing of these opening and closing does not depend solely on a single water level or level difference, but rather on the tidal processes and the safe water exchange level range both inside and outside the water system.

[0003] Numerical simulation is an important research tool for water exchange. However, existing mathematical models for water exchange have relatively simple gate control, making it difficult to automate the frequent opening and closing of gates in the aforementioned water system under specific logical conditions. In the MIKE software, widely used in estuary and coastal areas, sluice gates have three control methods:

[0004] (1) The gate is controlled to open and close by setting the opening degree k (the value ranges between 0 and 1, where 0 represents closed, 1 represents fully open, and a fraction between 0 and 1 represents the gate being partially open);

[0005] (2) Control the opening and closing of the gate according to the water level at the given location;

[0006] (3) Control the opening and closing of the gate according to the water level difference between the two given locations.

[0007] None of these three control methods can achieve the automatic opening and closing of the gates in the aforementioned estuary-coastal water system under the condition of safe water exchange level, based on the tidal process inside and outside the system and the safe water exchange level range, which greatly affects the efficiency of water exchange numerical simulation. Summary of the Invention

[0008] Objective: In view of at least one of the above technical problems, the present invention provides a method for simulating water exchange in estuary and coastal gate control systems under given high water exchange level (upper limit of safe water exchange range) and low water exchange level (lower limit of safe water exchange range).

[0009] This invention decomposes the water exchange period into multiple time periods and simulates them sequentially by time. By extracting and analyzing the model operation results of the previous time period, the opening and closing of the water system gates are determined, realizing the simulation of water exchange under complex scheduling of multi-gate water systems under given high and low water exchange levels. This method achieves automated simulation of water exchange under complex scheduling of sluice gates in estuarine and coastal sluice-controlled water systems through repetitive programmable operations such as data extraction, intersection point finding, and modification of model configuration files. It avoids human intervention and interference during the simulation process, and has advantages such as high efficiency and wide applicability, showing broad application prospects in coastal city construction, tidal flat and island development and utilization, and other fields.

[0010] The technical solution adopted in this invention is as follows:

[0011] Firstly, this application provides a method for simulating water exchange in estuarine and coastal gate-controlled water systems, including:

[0012] S1. Determine the high and low water exchange levels of the estuary and coastal water system based on the upper and lower limits of the safe water exchange level range.

[0013] S2. Set up n water level monitoring points inside the water system and m water level monitoring points outside the water system, respectively;

[0014] S3. Establish a mathematical model of water exchange encompassing the water system, and set the model configuration file, including: setting the opening degree of all gates to 0, setting the initial water level of the model between the high and low water exchange levels, and setting the initial concentration of transported substances within and outside the water system to 0; run the model until the hydrodynamic field stabilizes, and record the running time as T. w ;

[0015] S4. Modify the configuration file of the water exchange mathematical model established in S3, and run the modified model. The running time is T1, T / 2≤T1≤3T / 4, where T is the tidal period. The modifications include: A. Extracting T from S3. w A. The hydrodynamic field and the water exchange and transport material field at any given time are used as the initial conditions for this model operation; B. The initial concentration of water exchange and transport material in the water system is set to 1 (initial concentration of water body in the water system, dimensionless), and the initial concentration outside the water system is set to 0 (initial concentration of water body outside the water system, dimensionless).

[0016] S5. Extract the water level time series from n water level monitoring points within the water system and m water level monitoring points outside the water system from the model execution results in S4. Average these time series to obtain the expected water level values ​​for the n water level monitoring points within the water system and the m water level monitoring points outside the water system. Based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection of the water level values ​​of the monitoring points within and outside the water system.int0 Modify the configuration file of the water exchange mathematical model established in S4, and run the modified model. The running time is T2, T / 2≤T2≤3T / 4. The modifications include: A. Extracting T from S4. int0 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening degree of the inlet gate to 1 and the opening degree of the outlet gate to 0;

[0017] Repeat the cyclic steps in sequence until the cumulative simulation time reaches the test requirements.

[0018] Furthermore, the cyclical steps include:

[0019] S6. Extract the water level time series of n water level monitoring points in the water system from the model running results of the previous step, and take the average to obtain the time series of the expected water level values ​​of the n water level monitoring points in the water system; based on the time series of the expected water level values ​​of the n water level monitoring points in the water system, determine the time T corresponding to the first intersection point of the water level of the water level monitoring points and the high water level of the water exchange. int1 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T3, where T / 2 ≤ T3 ≤ 3T / 4. The modifications include: A. Extracting the T value from the model running results of the previous step. int1 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening of the inlet gate to 0 and the opening of the outlet gate to 0.

[0020] S7. Extract the water level time series from the n water level monitoring points within the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the n water level monitoring points within the water system; extract the water level time series from the m water level monitoring points outside the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the m water level monitoring points outside the water system; based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the time series of the expected water level values ​​for the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system. int2 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T4, where T / 2 ≤ T4 ≤ 3T / 4. The modifications include: A. Extracting the T from the previous step. int2 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening of the inlet gate to 0 and the opening of the outlet gate to 1.

[0021] S8. Extract the water level time series of n water level monitoring points in the water system from the model running results of the previous step, and take the average to obtain the time series of the expected water level values ​​of the n water level monitoring points in the water system; based on the time series of the expected water level values ​​of the n water level monitoring points in the water system, determine the time T corresponding to the first intersection point between the water level of the water level monitoring points and the low water level of the water exchange. int3 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T5, where T / 2 ≤ T5 ≤ 3T / 4. The modifications include: A. Extracting the T value from the model running results of the previous step. int3 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening of the inlet gate to 0 and the opening of the outlet gate to 0.

[0022] S9. Extract the water level time series from the n water level monitoring points within the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the n water level monitoring points within the water system; extract the water level time series from the m water level monitoring points outside the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the m water level monitoring points outside the water system; based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the time series of the expected water level values ​​for the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system. int4 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T6, where T / 2 ≤ T6 ≤ 3T / 4. The modifications include: A. Extracting the T value from the model running results of the previous step. int4 A. The hydrodynamic field and the water exchange and transport material field at any given time are used as the initial conditions for this model operation; B. Set the opening degree of the inlet gate to 1 and the opening degree of the outlet gate to 0.

[0023] Secondly, this application provides a water exchange simulation device for a river estuary and coastal gate control system, including a processor and a storage medium;

[0024] The storage medium is used to store instructions;

[0025] The processor is configured to operate according to the instructions to execute the method according to the first aspect.

[0026] Thirdly, this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method described in the first aspect.

[0027] Fourthly, this application provides a computer device including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the method described in the first aspect.

[0028] Fifthly, this application provides a computer program product, including a computer program that, when executed by a processor, implements the method described in the first aspect.

[0029] Beneficial effects: The water exchange simulation method for estuary and coastal gate control water systems provided by this invention automatically controls the opening and closing of inlet and outlet gates through logical judgment. It realizes the automated simulation of water exchange during complex scheduling of multi-gate water systems under given safe high water exchange and low water exchange conditions. It avoids human intervention and interference in the simulation process and has the advantages of high efficiency and wide applicability. It has broad application prospects in the fields of coastal city construction and (artificial) island development and utilization. Attached Figure Description

[0030] Figure 1 This is a schematic diagram of the layout of a gate control system in a coastal new city according to an embodiment of the present invention;

[0031] Figure 2 This is a schematic diagram illustrating the relative relationship between the high and low water levels during water exchange, the water exchange process, the non-water exchange process, and the external tidal level process in an embodiment of the present invention.

[0032] Figure 3 This is a schematic diagram illustrating the change of water level over time at an external water level monitoring point in the water system within the mathematical model of water exchange in this embodiment of the invention.

[0033] Figure 4 This is a schematic diagram of the water level process at the S5 water system monitoring point in an embodiment of the present invention;

[0034] Figure 5 This is a schematic diagram of the water level process at the monitoring point of the S6 water system in an embodiment of the present invention;

[0035] Figure 6 This is a schematic diagram of the water level process at the S7 water system monitoring point in an embodiment of the present invention;

[0036] Figure 7 This is a schematic diagram of the water level process at the monitoring point of the S8 water system in an embodiment of the present invention;

[0037] Figure 8 This is a schematic diagram of the water level process at the S9 water system monitoring point in an embodiment of the present invention;

[0038] Figure 9 This is a schematic diagram of the water change effect on the 5th day in this embodiment of the invention;

[0039] Figure 10 This is a schematic diagram of the water change effect after the 10th day in this embodiment of the invention. Detailed Implementation

[0040] The present invention will be further described below with reference to the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and should not be used to limit the scope of protection of the present invention.

[0041] In the description of this application, "several" means one or more, "multiple" means two or more, "greater than," "less than," and "exceeding" are understood to exclude the stated number, while "above," "below," and "within" are understood to include the stated number. The use of "first" and "second" in the description is merely for distinguishing technical features and should not be construed as indicating or implying relative importance, or implicitly indicating the number of indicated technical features, or implicitly indicating the order of the indicated technical features.

[0042] In the description of this application, the terms "one embodiment," "some embodiments," "illustrative embodiment," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples.

[0043] The term "and / or" simply describes the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, or B alone. Additionally, the character " / " generally indicates that the preceding and following related objects have an "or" relationship.

[0044] Example 1: This example provides a method for simulating water exchange in estuarine and coastal sluice gate control systems, such as... Figure 1 As shown, it includes:

[0045] S1. Determine the high water level h for water exchange in the estuary-coastal water system based on the upper and lower limits of the safe water exchange water level range. h With water exchange at low water level h l ;

[0046] S2. Set up n water level monitoring points inside the water system and m water level monitoring points outside the water system, where n≥1 and m≥1 respectively;

[0047] S3. Establish a mathematical model of water exchange encompassing the water system. Configure the model configuration file as follows: set all gate openings to 0 (closed); set the initial water level between the high and low water exchange levels; and set the initial concentrations of transported substances within and outside the water system to 0. Run the model until the hydrodynamic field stabilizes, and record the running time as T. w T w The water level outside the water system at the corresponding time should be lower than T.w The water level in the water system at that moment;

[0048] S4. Modify the configuration file of the water exchange mathematical model established in S3, and run the modified model. The running time is T1, T / 2≤T1≤3T / 4, where T is the tidal period. The modifications include the following two points:

[0049] A. Extract T from S3 w The hydrodynamic field and the water exchange and transport mass field at any given time are used as the initial conditions for this model operation; the hydrodynamic field includes the total water depth, the velocity component in the X direction, and the velocity component in the Y direction, and the water exchange and transport mass field includes the concentration.

[0050] B. Set the initial concentration of the water exchange and transport substances within the water system to 1 (initial concentration of the water body within the water system, dimensionless), and set the initial concentration outside the water system to 0 (initial concentration of the water body outside the water system, dimensionless).

[0051] S5. Extract the water level time series from n water level monitoring points within the water system from the model execution results in S4, and average them to obtain the time series of the expected water level values ​​for the n water level monitoring points within the water system; extract the water level time series from m water level monitoring points outside the water system from the model execution results in S4, and average them to obtain the time series of the expected water level values ​​for the m water level monitoring points outside the water system; based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the time series of the expected water level values ​​for the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection point of the water levels of the water level monitoring points inside and outside the water system. int0 ;

[0052] Modify the configuration file for the water exchange mathematical model established in S4, and run the modified model. The running time is T2, where T / 2 ≤ T2 ≤ 3T / 4. The modifications include the following two points:

[0053] A. Extract T from S4 int0 The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation;

[0054] B. Set the inlet gate opening to 1, indicating a fully open state, and the outlet gate opening to 0;

[0055] Repeat steps S6, S7, S8, and S9 in sequence until the cumulative simulation time meets the test requirements:

[0056] S6. Extract the water level time series of n water level monitoring points in the water system from the model running results of the previous step, and take the average to obtain the time series of the expected water level values ​​of the n water level monitoring points in the water system; based on the time series of the expected water level values ​​of the n water level monitoring points in the water system, determine the time T corresponding to the first intersection point of the water level of the water level monitoring points and the high water level of the water exchange.int1 ;

[0057] Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T3, where T / 2 ≤ T3 ≤ 3T / 4. The modifications include the following two points:

[0058] A. Extract T from the model execution results of the previous step. int1 The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation;

[0059] B. Set the inlet gate opening to 0 and the outlet gate opening to 0;

[0060] S7. Extract the water level time series from the n water level monitoring points within the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the n water level monitoring points within the water system; extract the water level time series from the m water level monitoring points outside the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the m water level monitoring points outside the water system; based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the time series of the expected water level values ​​for the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system. int2 ;

[0061] Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T4, where T / 2 ≤ T4 ≤ 3T / 4. The modifications include the following two points:

[0062] A. Extract T from the previous step int2 The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation;

[0063] B. Set the inlet gate opening to 0 and the outlet gate opening to 1;

[0064] S8. Extract the water level time series of n water level monitoring points in the water system from the model running results of the previous step, and take the average to obtain the time series of the expected water level values ​​of the n water level monitoring points in the water system; based on the time series of the expected water level values ​​of the n water level monitoring points in the water system, determine the time T corresponding to the first intersection point between the water level of the water level monitoring points and the low water level of the water exchange. int3 ;

[0065] Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T5, where T / 2 ≤ T5 ≤ 3T / 4. The modifications include the following two points:

[0066] A. Extract T from the model execution results of the previous step. int3 The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation;

[0067] B. Set the inlet gate opening to 0 and the outlet gate opening to 0;

[0068] S9. Extract the water level time series from the n water level monitoring points within the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the n water level monitoring points within the water system; extract the water level time series from the m water level monitoring points outside the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the m water level monitoring points outside the water system; based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the time series of the expected water level values ​​for the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system. int4 ;

[0069] Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T6, where T / 2 ≤ T6 ≤ 3T / 4. The modifications include the following two points:

[0070] A. Extract T from the model execution results of the previous step. int4 The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation;

[0071] B. Set the inlet gate opening to 1 and the outlet gate opening to 0.

[0072] Furthermore, in S1, the high water level during water exchange is the upper limit of the safe water exchange water level range, and the low water level during water exchange is the lower limit of the safe water exchange water level range; the high water level during water exchange h h The low tide level should be lower than the lowest tide level measured, predicted, or subsequently reported outside the water system, while the low tide level for water exchange should be higher than the highest low tide level measured, predicted, or subsequently reported outside the water system.

[0073] Furthermore, in S6, the cumulative simulation time... The calculation starts from the beginning of the S4 water exchange and is performed as follows:

[0074] ;

[0075] In the formula, x represents the number of times the model is run from the start of the S4 water exchange to the end of the simulation, which is the total number of model runs minus 1.

[0076] Furthermore, in some embodiments, if S5 is not a continuation of S4 calculation, but a continuation calculation that is interrupted after several calculations, then it is necessary to determine the water system state calculated in the previous step. The water system state includes: (1) the external tide level of the water system is in the process of rising tide and the water system is in the process of opening the gate to receive water; (2) the external tide level of the water system is in the process of rising tide and the water system is in the process of closing the gate; (3) the external tide level of the water system is in the process of falling tide and the water system is in the process of opening the gate to receive water; (4) the external tide level of the water system is in the process of falling tide and the water system is in the process of closing the gate.

[0077] If the previous calculation was performed when the external tide level of the water system was in the process of rising tide and the water system was in the process of opening the gate to receive water, then the current calculation is performed when the water system is in the process of closing the gate. Repeat steps S6, S7, S8 and S9 in sequence until the simulated cumulative time reaches the test requirements.

[0078] If the previous calculation was performed when the external tide level of the water system was in the process of rising tide and the water system was in the process of shutting off the gate, then the current calculation is performed when the water system is in the process of opening the gate and releasing water. Repeat steps S7, S8, S9 and S6 in sequence until the simulated cumulative time reaches the test requirements.

[0079] If the previous calculation was that the external tide level of the water system was in the process of ebb tide and the internal water system was in the process of opening the gate to release water, then the current calculation is the process of the water system closing the gate. Repeat S8, S9, S6 and S7 in sequence until the simulated cumulative time reaches the test requirements.

[0080] If the previous calculation was performed when the external tide level of the water system was in the process of ebb tide and the internal water system was in the process of gate closure, then the current calculation is performed when the water system is in the process of gate opening and water intake. Repeat steps S9, S6, S7 and S8 in sequence until the simulated cumulative time reaches the test requirements.

[0081] Simulation verification examples: such as Figure 1As shown, a coastal new city consists of nine artificial islands with a crisscrossing network of rivers and lakes. Narrow channels form between islands 1 and 5, 5 and 6, 1 to 4, and 6 to 9. For moisture control, four water exchange gates (Gate 1, Gate 2, Gate 3, and Gate 4) are installed at these four locations. When these gates are closed, the waterways and channels between the islands effectively form a closed system, negatively impacting the internal water quality and environment. To improve the water quality and environment within the system, the management has stipulated that water exchange can be conducted within a safe water exchange zone by opening and closing these four gates using tidal forces. Specifically, Gates 1 and 2 are inlet gates, and Gates 3 and 4 are outlet gates. The safe water exchange range is -0.25m to 0.25m (all water levels in this case are calculated from the local mean sea level, the same below), that is, the safe high water level for water exchange is 0.25m, and the safe low water level for water exchange is -0.25m. The water exchange process is described as follows: During high tide, when the water level inside the water system is higher than the safe low water level for water exchange, the inlet gate is opened and the outlet gate is closed, allowing water from outside the system to enter the water system through tidal forces; when the water level inside the water system is equal to the safe high water level for water exchange, for safety reasons, all gates are closed; during low tide, when the water level outside the system is equal to the water level inside the system, the outlet gate is opened and the inlet gate is closed, allowing water from inside the water system to flow out of the system under tidal action; when the water level inside the system is equal to the safe low water level for water exchange, all gates are closed, and water exchange stops. The relative relationships between the high water level for water exchange, the low water level for water exchange, and the water exchange and non-water exchange processes and the tidal process outside the system are shown in [reference needed]. Figure 2 .

[0082] To simulate the water exchange process described above, the simulation should be conducted according to the following steps:

[0083] S1: Determine the safe high water level h for water exchange in the estuary-coastal water system h With safe water exchange at low water level h l The values ​​are 0.25m and -0.25m respectively. During the water exchange simulation, the lowest high tide level outside the water system is higher than 0.25m, and the highest low tide level is lower than -0.25m. The safe water exchange high water level is h. h With safe water exchange at low water level h l Between the lowest high tide level and the highest low tide level outside the system;

[0084] S2: Nine water level monitoring points are set up within the water system, and one water level monitoring point is set up outside the water system. The locations of the monitoring points are as follows: Figure 1 As shown;

[0085] S3: Establish a mathematical model configuration file (S02TJ-0.m21fm) for water exchange, encompassing the water system. Start time: 2013-10-17 00:00:00, time step: 30s, calculation steps: 2820, end time: 2013-10-17 23:30:00; gates 1-4 are fully closed (i.e., gate opening is set to 0), initial water level is 0m, run the model until the hydrodynamic field stabilizes, and observe the water level changes over time at external monitoring points. Figure 3 As shown, the water level change curve tends to be smooth, indicating that the hydrodynamic field of the model is stable.

[0086] S4:

[0087] (1) Using the Data Extraction FM function of Mike Zero software, extract the hydrodynamic field of the model running result at 23:00:00 on October 17, 2013 in S3, including three parameters: total water depth, X-direction velocity component, and Y-direction velocity component, and save it as HYD_ini1.dfsu;

[0088] (2) Using the Data Extraction FM function of Mike Zero software, extract the water exchange transport material field at 23:00:00 on October 17, 2013 in S3, including the concentration, which is a total of 1 parameter. Set the initial concentration of the transport material inside the water system to 1 and the initial concentration outside the water system to 0, and save it as TSS_ini1.dfsu.

[0089] (3) Copy the water exchange mathematical model configuration file established in S3 (rename the file to S02TJ-1.m21fm) and make the following modifications: start time 2013-10-17 23:00:00, time step 30s, number of calculation steps 250, end time 2013-10-18 01:05:00; import the hot start field files of the hydrodynamic and material transport calculation modules respectively, namely HYD_ini1.dfsu and TSS_ini1.dfsu; the water exchange gates 1~4 are still in a completely closed state (i.e., the gate opening control factor is set to 0). Start the calculation with S02TJ-1.m21fm.

[0090] S5: Extract the water level time series of 9 water level monitoring points within the water system from the model execution results in S4, take the average, and obtain the time series of the expected water level values ​​of the 9 water level monitoring points within the water system; extract the water level time series of 1 water level monitoring point outside the water system from the model execution results in S4, and obtain the time series of the expected water level value of the 1 water level monitoring point outside the water system; write a program to automatically determine the time corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system, which is 2013-10-17 23:37:30. Figure 4 As shown.

[0091] (1) Using the Data Extraction FM function of Mike Zero software, the hydrodynamic field and water exchange and transport material field at 23:37:30 on October 17, 2013 were extracted from the model running results in S4 and saved as HYD_ini2.dfsu and TSS_ini2.dfsu, respectively;

[0092] (2) Copy the water exchange mathematical model configuration file established in S4 (rename the file to S02TJ-2.m21fm) and make the following modifications: start time 2013-10-17 23:37:30, time step 30s, number of calculation steps 500, end time 2013-10-18 03:47:30; import the hot start field files of the hydrodynamic and material transport calculation modules respectively, namely HYD_ini2.dfsu and TSS_ini2.dfsu; water exchange gates 1 and 2 are in the open state (i.e., the gate opening is set to 1), and water exchange gates 3 and 4 are still in the completely closed state (i.e., the gate opening is set to 0). Start the calculation with S02TJ-2.m21fm.

[0093] S6:

[0094] (1) Extract the water level time series of 9 water level monitoring points in the water system from the model operation results of the previous step, take the average, and obtain the time series of the expected water level values ​​of n water level monitoring points in the water system; determine whether the time series of the expected water level values ​​of the 9 water level monitoring points in the water system reaches the safe water exchange high water level h. h The time 0.25m corresponds to 2013-10-18 00:37:30. (For example...) Figure 5 As shown.

[0095] (2) Using the Data Extraction FM function of Mike Zero software, extract the hydrodynamic field and water exchange and transport material field at 00:37:30 on October 18, 2013, respectively, and save them as HYD_ini3.dfsu and TSS_ini3.dfsu respectively;

[0096] (2) Copy the water exchange mathematical model configuration file established in S4 (rename the file to S02TJ-3.m21fm) and make the following modifications: start time 2013-10-18 00:37:30, time step 30s, number of calculation steps 750, end time 2013-10-18 06:52:30; import the hot start field files of the hydrodynamic and material transport calculation modules respectively, namely HYD_ini3.dfsu and TSS_ini3.dfsu; water exchange gates 1 to 4 are all in a completely closed state (i.e., the gate opening is set to 0). Start the calculation with S02TJ-3.m21fm.

[0097] S7:

[0098] (1) Extract the water level time series of 9 water level monitoring points within the water system from the previous model operation results, take the average, and obtain the time series of the expected water level values ​​of the 9 water level monitoring points within the water system; extract the water level time series of 1 water level monitoring point outside the water system from the previous model operation results, and obtain the time series of the expected water level value of the 1 water level monitoring point outside the water system; compile a program to automatically determine the time corresponding to the first intersection of the time series of the expected water level values ​​of the water level monitoring points inside and outside the water system, that is, the time when the water level of the water level monitoring point outside the water system reaches the high water level of 0.25m during the ebb tide, which is 2013-10-18 05:49:30. Figure 6 As shown.

[0099] (2) Using the Data Extraction FM function of Mike Zero software, extract the hydrodynamic field and water exchange and transport material field at 05:49:30 on October 18, 2013, respectively, and save them as HYD_ini4.dfsu and TSS_ini4.dfsu respectively;

[0100] (3) Copy the water exchange mathematical model configuration file established in S4 (rename the file to S02TJ-4.m21fm) and make the following modifications: start time 2013-10-18 05:49:30, time step 30s, number of calculation steps 500, end time 2013-10-18 09:59:30; import the hot start field files of the hydrodynamic and material transport calculation modules respectively, namely HYD_ini4.dfsu and TSS_ini4.dfsu; water exchange gates 1 and 2 are both in a completely closed state (i.e., the gate opening is set to 0), and water exchange gates 3 and 4 are in an open state (i.e., the gate opening is set to 1). Start the calculation with S02TJ-4.m21fm.

[0101] S8:

[0102] (1) Extract the water level time series of 9 water level monitoring points in the water system from the model operation results of the previous step, take the average, and obtain the time series of the expected water level values ​​of the 9 water level monitoring points in the water system. Compile a program to automatically determine the time corresponding to the first intersection of the time series of the expected water level values ​​of the water level monitoring points in the water system with the low water level of the water exchange, that is, the time when the average water level of the 9 water level monitoring points in the water system reaches the safe low water level of the water exchange -0.25m, which is 07:34:00 on 2013-10-18. Figure 7 As shown.

[0103] (2) Using the Data Extraction FM function of Mike Zero software, extract the hydrodynamic field and water exchange and transport material field at 07:34:00 on October 18, 2013, respectively, and save them as HYD_ini5.dfsu and TSS_ini5.dfsu respectively;

[0104] (3) Copy the water exchange mathematical model configuration file established in S4 (rename the file to S02TJ-5.m21fm) and make the following modifications: start time 2013-10-18 07:34:00, time step 30s, number of calculation steps 750, end time 2013-10-18 13:49:00; import the hot start field files of the hydrodynamic and material transport calculation modules respectively, namely HYD_ini5.dfsu and TSS_ini5.dfsu; water exchange gates 1 to 4 are all in a completely closed state (i.e., the gate opening is set to 0). Start the calculation with S02TJ-5.m21fm.

[0105] S9:

[0106] (1) Extract the water level time series of 9 water level monitoring points within the water system from the model operation results in the previous step, take the average, and obtain the time series of the expected water level values ​​of the 9 water level monitoring points within the water system; extract the water level time series of 1 water level monitoring point outside the water system from the model operation results in the previous step, and obtain the time series of the expected water level value of the 1 water level monitoring point outside the water system; compile a program to automatically determine the time corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system, which is 2013-10-18 11:45:00. Figure 8 As shown.

[0107] (2) Using the Data Extraction FM function of Mike Zero software, extract the hydrodynamic field and water exchange and transport material field at 11:45:00 on October 18, 2013, respectively, and save them as HYD_ini6.dfsu and TSS_ini6.dfsu respectively;

[0108] (3) Copy the water exchange mathematical model configuration file established in S4 (rename the file to S02TJ-6.m21fm) and make the following modifications: start time 2013-10-18 11:45:00, time step 30s, number of calculation steps 750, end time 2013-10-18 18:00:00; import the hot start field files of the hydrodynamic and material transport calculation modules respectively, namely HYD_ini6.dfsu and TSS_ini6.dfsu; water exchange gates 1 and 2 are fully open (i.e., the gate opening is set to 1), and water exchange gates 3 and 4 are fully closed (i.e., the gate opening is set to 0). Start the calculation with S02TJ-6.m21fm.

[0109] Repeat steps S6, S7, S8, and S9 sequentially several times until the water system achieves the expected water exchange effect, then stop the calculation. It should be noted that the file copying, parameter settings, water level determination, time determination, and calculation initiation involved in steps S6 through S9 are all automated through computer programming, avoiding human intervention and interference during the simulation process.

[0110] The water change results on day 5 and day 10 in this case are shown below. Figure 9 , Figure 10 .

[0111] Example 2: Based on Example 1, this example provides a water exchange simulation device for a estuary and coastal gate control system, including a processor and a storage medium;

[0112] The storage medium is used to store instructions;

[0113] The processor is configured to operate according to the instructions to execute the method according to Embodiment 1.

[0114] Example 3: Based on Example 1, this example provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method described in Example 1.

[0115] Example 4: Based on Example 1, this example provides a computer device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the method described in Example 1.

[0116] Example 5: Based on Example 1, this example provides a computer program product, including a computer program that, when executed by a processor, implements the method described in Example 1.

[0117] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0118] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, as well as combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0119] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0120] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0121] The above description is only a preferred embodiment of this application. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principle of this application, and these improvements and modifications should also be considered within the scope of protection of this application.

Claims

1. A method for simulating water exchange in estuarine and coastal sluice gate control systems, characterized in that, include: S1. Determine the high and low water exchange levels of the estuary and coastal water system based on the upper and lower limits of the safe water exchange level range. S2. Set up n water level monitoring points inside the water system and m water level monitoring points outside the water system, respectively; S3. Establish a mathematical model of water exchange encompassing the water system, and set the model configuration file, including: setting the opening degree of all gates to 0, setting the initial water level of the model between the high and low water exchange levels, and setting the initial concentration of transported substances within and outside the water system to 0; run the model until the hydrodynamic field stabilizes, and record the running time as T. w ; S4. Modify the configuration file of the water exchange mathematical model established in S3, and run the modified model. The running time is T1, T / 2≤T1≤3T / 4, where T is the tidal period. The modifications include: A. Extracting T from S3. w A. The hydrodynamic field and the water exchange and transport material field at any given time are used as the initial conditions for this model operation; B. The initial concentration of water exchange and transport material in the water system is set to 1, and the initial concentration outside the water system is set to 0; S5. Extract the water level time series from n water level monitoring points within the water system and m water level monitoring points outside the water system from the model execution results in S4. Take the average of these two time series to obtain the expected water level values ​​for the n water level monitoring points within the water system and the m water level monitoring points outside the water system. Determine the time T corresponding to the first intersection of the water level values ​​of the monitoring points within and outside the water system. int0 Modify the configuration file of the water exchange mathematical model established in S4, and run the modified model. The running time is T2, T / 2≤T2≤3T / 4. The modifications include: A. Extracting T from S4. int0 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening degree of the inlet gate to 1 and the opening degree of the outlet gate to 0; Repeat the cyclic steps in sequence until the simulated cumulative time reaches the test requirements. The cyclic steps include: S6. Extract the water level time series of n water level monitoring points in the water system from the model running results of the previous step, and take the average to obtain the time series of the expected water level values ​​of the n water level monitoring points in the water system; based on the time series of the expected water level values ​​of the n water level monitoring points in the water system, determine the time T corresponding to the first intersection point of the water level of the water level monitoring points and the high water level of the water exchange. int1 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T3, where T / 2 ≤ T3 ≤ 3T / 4. The modifications include: A. Extracting the T value from the model running results of the previous step. int1 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening of the inlet gate to 0 and the opening of the outlet gate to 0. S7. Extract the water level time series from the n water level monitoring points within the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the n water level monitoring points within the water system; extract the water level time series from the m water level monitoring points outside the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the m water level monitoring points outside the water system; based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the time series of the expected water level values ​​for the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system. int2 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T4, where T / 2 ≤ T4 ≤ 3T / 4. The modifications include: A. Extracting the T from the previous step. int2 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening of the inlet gate to 0 and the opening of the outlet gate to 1. S8. Extract the water level time series of n water level monitoring points in the water system from the model running results of the previous step, and take the average to obtain the time series of the expected water level values ​​of the n water level monitoring points in the water system; based on the time series of the expected water level values ​​of the n water level monitoring points in the water system, determine the time T corresponding to the first intersection point between the water level of the water level monitoring points and the low water level of the water exchange. int3 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T5, where T / 2 ≤ T5 ≤ 3T / 4. The modifications include: A. Extracting the T value from the model running results of the previous step. int3 A. The hydrodynamic field and the water exchange and transport field at any given time are used as the initial conditions for this model operation; B. Set the opening of the inlet gate to 0 and the opening of the outlet gate to 0. S9. Extract the water level time series from the n water level monitoring points within the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the n water level monitoring points within the water system; extract the water level time series from the m water level monitoring points outside the water system from the model running results of the previous step, and average them to obtain the time series of the expected water level values ​​for the m water level monitoring points outside the water system; based on the time series of the expected water level values ​​for the n water level monitoring points within the water system and the time series of the expected water level values ​​for the m water level monitoring points outside the water system, determine the time T corresponding to the first intersection of the water levels of the water level monitoring points inside and outside the water system. int4 Modify the configuration file for the water exchange mathematical model established in S4, and run the model. The running time is T6, where T / 2 ≤ T6 ≤ 3T / 4. The modifications include: A. Extracting the T value from the model running results of the previous step. int4 A. The hydrodynamic field and the water exchange and transport material field at any given time are used as the initial conditions for this model operation; B. Set the opening degree of the inlet gate to 1 and the opening degree of the outlet gate to 0.

2. The water exchange simulation method for estuary and coastal gate control systems according to claim 1, characterized in that, In S1, the high water level during water exchange is the upper limit of the safe water exchange range, and the low water level during water exchange is the lower limit of the safe water exchange range; the high water level during water exchange h h The water level should be lower than the lowest high tide level measured, predicted, or subsequently reported outside the water system, and the water level should be changed to a lower level h. l It should be higher than the highest low tide level measured, predicted, or subsequently reported outside the water system.

3. The water exchange simulation method for estuary and coastal gate control systems according to claim 1, characterized in that, In S3, T w The water level outside the water system at the corresponding time should be lower than T. w The water level in the water system at that moment.

4. The water exchange simulation method for estuary and coastal gate control systems according to claim 1, characterized in that, A gate opening of 0 indicates a closed state; a gate opening of 1 indicates a fully open state.

5. The method for simulating water exchange in an estuary-coastal gate control system according to claim 1, characterized in that, The hydrodynamic field includes the total water depth, the velocity component in the X direction, and the velocity component in the Y direction. The water exchange and transport of substances field includes the concentration.

6. The water exchange simulation method for estuary and coastal gate control systems according to claim 1, characterized in that, In S6, the cumulative time is simulated. The calculation starts from the beginning of the S4 water exchange and is performed as follows: ； In the formula, x represents the number of times the model is run from the start of the S4 water exchange to the end of the simulation, which is the total number of model runs minus 1.

7. The water exchange simulation method for estuary and coastal gate control systems according to claim 1, characterized in that, If S5 is not a continuation of S4, but a continuation of the calculation after being interrupted several times, then it is necessary to determine the state of the water system calculated in the previous step. The state of the water system includes: (1) the external tide level of the water system is in the process of rising tide and the water system is in the process of opening the gate to receive water; (2) the external tide level of the water system is in the process of rising tide and the water system is in the process of closing the gate; (3) the external tide level of the water system is in the process of falling tide and the water system is in the process of opening the gate to receive water; (4) the external tide level of the water system is in the process of falling tide and the water system is in the process of closing the gate. If the previous calculation was performed when the external tide level of the water system was in the process of rising tide and the water system was in the process of opening the gate to receive water, then the current calculation is performed when the water system is in the process of closing the gate. Repeat steps S6, S7, S8 and S9 in sequence until the simulated cumulative time reaches the test requirements. If the previous calculation was performed when the external tide level of the water system was in the process of rising tide and the water system was in the process of shutting off the gate, then the current calculation is performed when the water system is in the process of opening the gate and releasing water. Repeat steps S7, S8, S9 and S6 in sequence until the simulated cumulative time reaches the test requirements. If the previous calculation was that the external tide level of the water system was in the process of ebb tide and the internal water system was in the process of opening the gate to release water, then the current calculation is the process of the water system closing the gate. Repeat S8, S9, S6 and S7 in sequence until the simulated cumulative time reaches the test requirements. If the previous calculation was performed when the external tide level of the water system was in the process of ebb tide and the internal water system was in the process of gate closure, then the current calculation is performed when the water system is in the process of gate opening and water intake. Repeat steps S9, S6, S7 and S8 in sequence until the simulated cumulative time reaches the test requirements.

8. A water exchange simulation device for a river estuary and coastal sluice gate control system, characterized in that, Including processor and storage media; The storage medium is used to store instructions; The processor is configured to operate according to the instructions to perform the method according to any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the method according to any one of claims 1 to 7.

Images