Scheduling method, system, device and storage medium for multi-line ship lock filling operation based on Fourier fast transform

Abstract

The invention discloses a scheduling method, system, equipment and storage medium for multi-line ship lock filling operation based on Fourier fast transform, belonging to the technical field of ship lock filling operation, the scheduling method comprises: respectively obtaining the main fluctuation curve of the pilot channel when each single-line ship lock is filled with water; superimposing the main fluctuation curves of all single-line ship locks to obtain the staggered operation time of multiple single-line ship locks when the total fluctuation of the pilot channel is minimum, and scheduling according to the staggered operation time; wherein, obtaining the main fluctuation curve of the pilot channel when the single-line ship lock is filled with water comprises: establishing a two-dimensional fluctuation mathematical model of the pilot channel upstream of the ship lock, obtaining the water surface fluctuation data of the pilot channel when the current single-line ship lock is filled with water under the initial water level of the pilot channel; establishing a Fourier fast change analysis model, obtaining the main fluctuation curve of the pilot channel when the current single-line ship lock is filled with water; the invention can control the fluctuation amplitude of the multi-line ship lock filling operation to ensure navigation safety in the pilot channel.

Description

Technical Field

[0001] The invention belongs to the technical field of ship lock water filling operation, and in particular relates to a scheduling method, system, equipment and storage medium for multi-line ship lock water filling operation based on Fourier fast transform. Background Art

[0002] With the growth of waterborne trade volume, the traffic density of the inland waterway network continues to rise; many single-line locks have been expanded and transformed into double-line locks or multi-line locks. Due to terrain restrictions, double-line locks or multi-line locks generally share pilot channels.

[0003] The multiple fluctuations in the pilot channel caused by the frequent operation of the ship locks chase and overlap each other, making the navigation flow conditions of the pilot channel extremely complex, and often facing the situation of reciprocating oscillation of the water level. In addition, the ship locks with double-line ship locks or multi-line ship locks sharing the pilot channel are also faced with more violent water surface fluctuations caused by the staggered operation or simultaneous operation of multiple-line ship locks. The violent water surface fluctuations have a huge impact on the navigation safety of the shared pilot channel. How to control the amplitude of the pilot channel water surface fluctuations when the multi-line ship locks or double-line ship locks are in operation is an important guarantee for ensuring the navigation safety of the double-line shared pilot channel. Summary of the invention

[0004] The purpose of the present invention is to overcome the deficiencies in the prior art and to provide a scheduling method, system, equipment and storage medium for multi-line ship lock filling operation based on Fourier fast transform, which can utilize the wave superposition law to control the fluctuation amplitude of the multi-line ship lock filling operation and ensure navigation safety in the pilot channel.

[0005] The present invention provides the following technical solutions:

[0006] In a first aspect, a scheduling method for water filling operation of a multi-line ship lock based on Fourier fast transform is provided, wherein the multi-line ship lock includes a plurality of single-line ship locks, and the scheduling method includes:

[0007] Obtain the main fluctuation curve of the navigation channel when each single-line ship lock is filled with water;

[0008] According to the superposition law of waveforms, the main fluctuation curves of all single-line ship locks are superimposed to obtain the staggered operation time of multiple single-line ship locks when the total fluctuation of the pilot channel is minimum, and the water filling operation of the multi-line ship locks is scheduled according to the obtained staggered operation time;

[0009] Among them, the main fluctuation curves of the pilot channel when the single-line ship lock is filled with water include:

[0010] Based on the current single-line ship lock water filling operation data, the relationship between the water level and water filling flow of the pilot channel, the pilot channel topographic data, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock water filling, a two-dimensional fluctuation mathematical model of the pilot channel is established;

[0011] The established two-dimensional wave mathematical model is used to obtain the water surface wave data of the pilot channel when the single-line ship lock is filled with water at the initial water level of the pilot channel.

[0012] Based on the water surface fluctuation data of the pilot channel, a Fourier rapid change analysis model is established to obtain the main fluctuations of the pilot channel when the single-line ship lock is filled with water at the initial water level of the pilot channel, as well as the frequency and fluctuation data of the main fluctuations;

[0013] According to the frequency and fluctuation data of the main fluctuations of the pilot channel, the main fluctuation curve of the pilot channel when the current single-line ship lock is filled with water is obtained.

[0014] Optionally, the current single-line ship lock filling operation data includes the current upstream and downstream water heads of the ship lock, the opening and operation time of the ship lock valve, the ship lock valve opening, the upstream water level of the ship lock and the lock chamber water level; the relationship between the pilot channel water level and the filling flow is designed based on the ship lock valve opening, the upstream water level of the ship lock and the lock chamber water level and verified according to measured data.

[0015] Optionally, the Fourier rapid change analysis model is established based on the water surface fluctuation data of the pilot channel, and the main fluctuations of the pilot channel when the current single-line ship lock is filled with water under the initial water level of the pilot channel, as well as the frequency and fluctuation data of the main fluctuations, are obtained, including:

[0016] Under the initial water level of the pilot channel, the water surface fluctuation data of the pilot channel when the current single-line ship lock is filled with water is standardized to obtain a water level fluctuation data set with equal intervals;

[0017] The water level fluctuation data set is transformed using Fourier fast transform to obtain a set of complex number arrays, each of which corresponds to a frequency and a modulus.

[0018] The fluctuations corresponding to the complex numbers whose modulus values ​​in the complex number array meet the requirements are taken as the main fluctuations of the pilot channel when the current single-line ship lock is filled with water, the fluctuation data corresponding to the complex numbers whose modulus values ​​meet the requirements are the fluctuation data of the main fluctuations of the pilot channel, and the corresponding frequency is the frequency of the main fluctuations of the pilot channel.

[0019] Optionally, in the method of using the fluctuations corresponding to the complex numbers whose modulus values ​​in the complex number array meet the requirements as the main fluctuations of the pilot channel when the current single-line ship lock is filled with water, the fluctuation data corresponding to the complex numbers whose modulus values ​​meet the requirements are the fluctuation data of the main fluctuations of the pilot channel, and the corresponding frequency is the frequency of the main fluctuations of the pilot channel,

[0020] The main fluctuation of the pilot channel when the current single-line ship lock is filled with water is the fluctuation corresponding to the complex number with the largest modulus value in the complex number array; the frequency of the main fluctuation of the pilot channel for:

[0021]

[0022] in, is the index corresponding to the complex number with the largest modulus value in the complex array, N is the number of data in the water level fluctuation dataset, and Δt is the sampling interval of the water level fluctuation dataset.

[0023] Optionally, the obtaining of the main fluctuation curve of the pilot channel when the single-line ship lock is currently filled with water according to the frequency and fluctuation data of the main fluctuation of the pilot channel comprises: determining the main fluctuation curve of the pilot channel when the single-line ship lock is currently filled with water according to the frequency and fluctuation data of the main fluctuation of the pilot channel by using a fluctuation formula;

[0024] The fluctuation formula is:

[0025]

[0026] Among them, Z(t) is the water level of the pilot channel, m is the number of major fluctuations in the pilot channel, is the modulus value of the main fluctuation; N is the number of data in the water level fluctuation data set, is the frequency of the main fluctuation, t is the time, and Z is the initial water level of the pilot channel.

[0027] Optionally, the two-dimensional fluctuation mathematical model of the pilot channel is established during the water filling process based on the operation data of the current single-line ship lock filling, the relationship between the water level and the filling flow of the pilot channel, the topographic data of the pilot channel, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock filling.

[0028] The two-dimensional wave mathematical model is constructed based on the river base flow and the non-steady flow; for the river base flow, the downstream flow of the pilot channel itself and its corresponding pilot channel water level are the initial boundaries; for the non-steady flow, when the ship lock is filled with water, the change process of the pilot channel flow is the boundary of the non-steady flow;

[0029] After the two-dimensional wave mathematical model is established, the calculation results of the two-dimensional wave mathematical model are verified using the water surface wave monitoring data of the navigation channel after the single-line ship lock is filled with water.

[0030] Optionally, in the two-dimensional fluctuation mathematical model established to obtain the water surface fluctuation data of the pilot channel when the current single-line ship lock is filled with water at the initial water level of the pilot channel, only the water surface fluctuation data of the current single-line ship lock within 4 hours of the initial water level of the pilot channel is obtained.

[0031] In a second aspect, a dispatching system for water filling operation of a multi-line ship lock based on Fourier fast transform is provided, wherein the multi-line ship lock includes a plurality of single-line ship locks, and the dispatching system includes:

[0032] The fluctuation curve acquisition module is used to respectively obtain the main fluctuation curve of the pilot channel when each single-line ship lock is filled with water;

[0033] The superposition scheduling module is used to superimpose the main fluctuation curves of all single-line ship locks according to the superposition law of waveforms, so as to obtain the staggered operation time of multiple single-line ship locks when the total fluctuation of the pilot channel is minimum, and to schedule the water filling operation of the multi-line ship locks according to the obtained staggered operation time;

[0034] Wherein, the fluctuation curve acquisition module includes a model building unit, a first fluctuation data acquisition unit, a second fluctuation data acquisition unit and a fluctuation curve generation unit;

[0035] A model building unit is used to establish a two-dimensional fluctuation mathematical model of the pilot channel based on the current single-line ship lock water filling operation data, the relationship between the pilot channel water level and the water filling flow, the pilot channel topographic data, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock is filled with water;

[0036] The first fluctuation data acquisition unit is used to obtain the water surface fluctuation data of the pilot channel when the single-line ship lock is filled with water at the initial water level of the pilot channel by using the established two-dimensional fluctuation mathematical model;

[0037] The second fluctuation data acquisition unit is used to establish a Fourier rapid change analysis model based on the water surface fluctuation data of the pilot channel, and obtain the main fluctuations of the pilot channel when the current single-line ship lock is filled with water under the initial water level of the pilot channel, as well as the frequency and fluctuation data of the main fluctuations;

[0038] The fluctuation curve generating unit is used to obtain the main fluctuation curve of the pilot channel when the single-line ship lock is currently filled with water according to the frequency and fluctuation data of the main fluctuation of the pilot channel.

[0039] In a third aspect, a computer device is provided, comprising a processor and a memory; wherein, when the processor executes a computer program stored in the memory, the steps of the scheduling method for multi-line ship lock filling operation based on Fourier fast transform as described in any one of the first aspects are implemented.

[0040] In a fourth aspect, a computer-readable storage medium is provided for storing a computer program; when the computer program is executed by a processor, the steps of the scheduling method for multi-line ship lock filling operation based on Fourier fast transform as described in any one of the first aspects are implemented.

[0041] Compared with the prior art, the present invention has the following beneficial effects:

[0042] The present invention combines the two-dimensional wave model of the upstream pilot channel of the lock with the Fourier rapid change analysis model, can fit the fluctuations caused by the filling of a single-line lock, and use the wave superposition law to determine the staggered operation time of multiple-line locks, thereby scheduling the filling operation of multiple-line locks to ensure the safety of ships when the multi-line shared pilot channel operates at staggered times. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 It is a schematic diagram of the steps of the scheduling method of the multi-line ship lock water filling operation based on Fourier fast transform of the present invention;

[0044] Figure 2 It is a schematic diagram of the steps of obtaining the main fluctuation curve when the single-line ship lock is filled with water according to the present invention;

[0045] Figure 3 is a graph showing the change of water level fluctuation over time in the berthing section of the upstream navigation channel after the ship lock is filled with water in the second embodiment of the present invention;

[0046] Figure 4 is a frequency and corresponding modulus distribution diagram of the second embodiment of the present invention;

[0047] Figure 5 The numerical simulation results of the second embodiment of the present invention are compared with the wave formula results obtained based on the Fourier fast transform;

[0048] Figure 6 is the maximum positive and negative wave heights under different staggered working conditions of the second embodiment of the present invention;

[0049] Figure 7 It is a structural block diagram of the dispatching system of the multi-line ship lock water filling operation based on Fourier fast transform of the present invention. DETAILED DESCRIPTION

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

[0051] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the term "comprising" and any of its variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0052] Example 1

[0053] like Figure 1 As shown, a scheduling method for water filling operation of a multi-line ship lock based on Fourier fast transform is provided, wherein the multi-line ship lock includes a plurality of single-line ship locks, and the scheduling method includes the following steps:

[0054] S100: Obtain the main fluctuation curve of the pilot channel when each single-line ship lock is filled with water.

[0055] Specifically, Figure 2 As shown in the figure, the main fluctuation curves of the navigation channel when the single-line ship lock is filled with water include:

[0056] S101: A two-dimensional fluctuation mathematical model of the pilot channel is established based on the current operation data of the single-line ship lock filling, the relationship between the water level and the filling flow of the pilot channel, the topographic data of the pilot channel, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock is filled with water.

[0057] Specifically, the current single-line ship lock water filling operation data includes the current upstream and downstream water heads h, the opening operation time t of the ship lock valve, and the o , lock valve opening, water level upstream of the lock and water level in the lock chamber; the relationship between the water level in the pilot channel and the filling flow rate is designed based on the lock valve opening, water level upstream of the lock and water level in the lock chamber and verified based on measured data.

[0058] The topographic data of the pilot channel usually includes the elevation and cross-sectional shape of the pilot channel.

[0059] According to the relationship between the water level of the pilot channel and the filling flow rate during the design filling process, the maximum flow rate Q during the single-line ship lock filling process is obtained. max ; Based on the measured topographic data of the upstream navigation channel of the ship lock, a two-dimensional wave mathematical model of the upstream navigation channel is established.

[0060] In the process of establishing the two-dimensional wave mathematical model of the upstream pilot channel, it is necessary to simulate the minimum navigable water level Z of the upstream pilot channel. min To the highest navigable water level Z max The fluctuation of the water level in the navigation channel after the ship lock is filled with water is calculated under different water levels between the two periods, and the variation process of the water level in the navigation channel after the ship lock is filled with water under different navigation water levels is obtained.

[0061] Since the waves generated by the filling of the ship lock are typical long gravity waves, and will be reflected at the upstream boundary and then transmitted back to the pilot channel; therefore, the position of the upstream boundary will affect the wave shape in the pilot channel; the position of the upstream boundary of the two-dimensional wave mathematical model is taken as 40 km upstream of the hub to reduce the impact of wave reflection on the wave of the pilot channel. The position of the upstream boundary makes the time for the waves caused by the filling of the ship lock to be reflected back to the pilot channel by the upstream boundary exceed 4 hours, so the calculation results of the two-dimensional wave mathematical model are more reliable within the first 4 hours.

[0062] The two-dimensional wave mathematical model is constructed based on the river base flow and non-steady flow. For the river base flow, the downstream flow of the pilot channel itself and its corresponding pilot channel water level are the initial boundaries. For the non-steady flow, the change process of the pilot channel flow when the lock is filled with water is the boundary of the non-steady flow.

[0063] Specifically, the calculation of unsteady flow in the two-dimensional wave mathematical model is based on the Reynolds mean stress equation for two-dimensional incompressible fluid and obeys the Boussinesq hypothesis and the hydrostatic pressure hypothesis.

[0064] More specifically, the governing equation for a two-dimensional unsteady flow is:

[0065] The basic equation is:

[0066] (1)

[0067] Where h is the total water head, is the river bottom elevation, and d is the still water depth.

[0068] The continuity equation describing the two-dimensional water flow in a plane is:

[0069] (2)

[0070] The momentum equation describing the two-dimensional water flow in a plane is:

[0071] (3)

[0072] (4)

[0073] ,

[0074] in, and It is the average velocity based on the water depth; and are the components of the vertical average velocity in the x and y directions, t is time, x, y, z are Cartesian coordinate systems, u and v are the velocity components in the x and y directions, g is the gravitational acceleration, is the relative density of water, , , and is the component of radiation stress; P a is the atmospheric pressure, is the relative density of water; S is the flow rate of the point source; and The flow rate of the source-sink flow is x and y Directional velocity component; and are the components of the shear stress of the water flow at the water surface in the x and y directions, and are the components of the riverbed shear stress in the x and y directions, and are the accelerations of the Earth's rotation in the x and y directions respectively.

[0075] Lateral stress It includes viscous friction, turbulent friction, and differential advection, and its value is estimated by the eddy viscosity formula based on the depth-averaged velocity gradient.

[0076] in, , ,

[0077] Where A is the flow cross-sectional area.

[0078] In the process of solving the base flow and non-steady flow of the river channel, the spatial discretization of the calculation area is to use the finite volume method to subdivide the continuous unity into non-overlapping units, which can be triangles or quadrilaterals.

[0079] The general form of the shallow water equations can generally be written as:

[0080] (5)

[0081] Among them, U is the conserved physical vector, F is the flux vector, and S is the source term.

[0082] In the Cartesian coordinate system, the two-dimensional shallow water equations can be written as:

[0083] (6)

[0084] Wherein, the superscripts I and V are the inviscid and viscous fluxes, respectively. and are the fluxes in the x and y directions respectively; integrating the i-th unit of the two-dimensional shallow water equations and rewriting them using the Gauss principle yields:

[0085] (7)

[0086] Where: n is the number of units, For unit The area of is the boundary of the unit; is the integral variable along the boundary. Here, single-point integration is used to calculate the area integral, and the integration point is located at the mass point of the unit. At the same time, the midpoint quadrature method is used to calculate the boundary integral. Formula (7) can be transformed into:

[0087] (8)

[0088] in, and are the i-th unit and The average value of and located in the unit center; is the number of boundaries of the cell; is the length of the jth unit.

[0089] Both first-order and second-order solutions can be used for spatial discretization. For two-dimensional cases, the approximate Riemann solution can be used to calculate the convective flow at the unit interface. When using the Roe method, the relevant variables on the left and right sides of the interface need to be estimated. In the second-order method, spatial accuracy can be obtained by using linear gradient techniques.

[0090] It is worth noting that the established two-dimensional wave mathematical model of the upstream pilot channel of the ship lock needs to be verified for effectiveness and used only after passing the verification. The specific verification method is: using the historical measured water level fluctuation data and comparing it with the water surface fluctuation data obtained by calculating the two-dimensional wave mathematical model, adjusting the calculation parameters such as the model roughness until the error between the calculated result and the measured data is less than 5%.

[0091] S102: using the established two-dimensional wave mathematical model, obtaining the water surface wave data of the pilot channel when the single-line ship lock is currently filled with water at the initial water level of the pilot channel.

[0092] Specifically, only the water surface fluctuation data of the current single-line ship lock under the initial water level of the pilot channel within 4 hours of water filling operation is obtained, so as to ensure the reliability of the obtained water surface fluctuation data; the water surface fluctuation data is the water level change data over time.

[0093] S103: Based on the water surface fluctuation data of the pilot channel, a Fourier rapid change analysis model is established, and the main fluctuations of the pilot channel under the initial water level of the pilot channel when the current single-line ship lock is filled with water, as well as the frequency and fluctuation data of the main fluctuations are obtained.

[0094] The specific steps include:

[0095] S1031: Standardize the water surface fluctuation data of the navigation channel when the front single-line ship lock is filled with water to obtain an equally spaced water level fluctuation data set.

[0096] Specifically, the time interval corresponding to each data in the water level fluctuation data set Z(t) is the same, and the number of data in the water level fluctuation data set N = 2 n . n is a constant.

[0097] S1032: Use Fourier fast transform to transform the water level fluctuation data set to obtain a set of complex number arrays, each complex number corresponds to a specific frequency, and each complex number corresponds to a modulus.

[0098] Specifically, the water level fluctuation dataset is time series data. The Fourier fast transform converts the time series data into frequency domain data. The water surface fluctuation dataset whose water level changes with time is transformed by the Fourier fast transform to obtain a set of complex number arrays Zi(t)=R(t)+I(t).

[0099] The absolute value of the complex array modulus obtained after Fourier rapid change analysis is the amplitude corresponding to the corresponding frequency.

[0100] The complex number obtained by Fourier fast transformation is converted to modulus using the following formula:

[0101] (9)

[0102] S1033: The fluctuations corresponding to the complex numbers whose modulus values ​​in the complex number array meet the requirements are taken as the main fluctuations of the pilot channel when the current single-line ship lock is filled with water. The fluctuation data corresponding to the complex numbers whose modulus values ​​meet the requirements are the fluctuation data of the main fluctuations of the pilot channel, and the corresponding frequency is the frequency of the main fluctuations of the pilot channel.

[0103] Specifically, the modulus values ​​of the complex arrays are compared. The larger the modulus value, the stronger the signal strength of this group of numbers in the frequency domain, that is, the fluctuation of the corresponding frequency is the main fluctuation of the lock discharge.

[0104] Specifically, the main fluctuation of the pilot channel when the current single-line ship lock is filled with water is the fluctuation corresponding to the complex number with the largest modulus value in the complex number array; the frequency of the main fluctuation of the pilot channel for:

[0105] (10)

[0106] in, is the index corresponding to the complex number with the largest modulus value in the complex array (starting from 0), N is the number of data in the water level fluctuation dataset, and Δt is the sampling interval of the water level fluctuation dataset.

[0107] Of course, in some other embodiments, multiple main fluctuations can be selected. The more main fluctuations are selected, the higher the accuracy of the simulation calculation will be. However, if too many main fluctuations are selected for each single-line lock, it will make it more difficult to superimpose multiple single-line locks.

[0108] S104: According to the frequency and fluctuation data of the main fluctuation of the pilot channel, a main fluctuation curve of the pilot channel when the single-line ship lock is currently filled with water is obtained.

[0109] Specifically, according to the frequency and fluctuation data of the main fluctuation of the pilot channel, the main fluctuation curve of the pilot channel when the current single-line ship lock is filled with water is determined using the fluctuation formula;

[0110] The fluctuation formula is:

[0111] (11)

[0112] Among them, Z(t) is the water level of the pilot channel, m is the number of major fluctuations, is the modulus value of the main fluctuation; N is the number of data in the water level fluctuation data set, is the frequency of the main fluctuation, t is the time, and Z is the initial water level of the pilot channel.

[0113] Specifically, the number and frequency of main fluctuations, the water level of main fluctuations and the initial water level of the approach channel are substituted into the fluctuation formula to obtain the main fluctuation curve when the current single-line ship lock is filled with water.

[0114] S200: According to the superposition law of waveforms, the main fluctuation curves of all single-line ship locks are superimposed to obtain the staggered operation time of multiple single-line ship locks when the total fluctuation of the pilot channel is minimum, and the water filling operation of the multi-line ship locks is scheduled according to the obtained staggered operation time.

[0115] Specifically, the wave height of the total fluctuation of the pilot channel is equal to the superposition of the wave heights of all one-way locks. The minimum total fluctuation can be the minimum amplitude. Of course, in some other embodiments, it can also be the minimum energy. Usually, in order to simplify the calculation, the minimum amplitude of the total fluctuation is selected. In order to achieve the minimum amplitude of the total fluctuation, it can be determined by traversing the staggered running time.

[0116] Specifically, as an example, the active fluctuation curves of two single-line ship locks are superimposed. The main fluctuation of the first single-line ship lock is Z1(t), and the main fluctuation of the second single-line ship lock is Z2(t). Then the total fluctuation Z(t) is:

[0117] Z(t)=Z1(t)+Z2(t+t0) (12)

[0118] Among them, t0 is the time of staggered operation.

[0119] Traverse the staggered operation time to obtain the optimal t0, and control the water filling operation time of the double-line ship lock according to the optimal staggered operation time.

[0120] For multi-line locks, the staggered operation times of multiple single-line locks are based on the operation time of the first single-line lock, and are scheduled according to the staggered operation time with the first single-line lock; it is worth noting that the staggered operation time can also be 0.

[0121] Example 2

[0122] Take the upstream water level of 145m and the working water head of the ship lock of 20.75m as an example. Figure 3As shown in the figure, the water level fluctuation of the upstream approach channel berthing section after the ship lock is filled with water is calculated using the two-dimensional wave mathematical model.

[0123] Perform equal time interpolation on the data of water level fluctuation over time so that the number of data satisfies N=2 n , in this example N=2 11 =2048.

[0124] like Figure 4 As shown in the figure, the water level fluctuation data was processed using the Fourier fast transform analysis model to obtain the frequency and corresponding modulus distribution diagram. According to the distribution data, the frequency corresponding to the maximum modulus value of 142 is the main frequency of the main fluctuation, and its corresponding wave height is 0.14m, and the corresponding frequency f i =0.0008, and its corresponding period is The calculated period T is about 20.83min. During the operation of the ship lock, the first fluctuation corresponding to the main frequency is the main factor affecting the water surface fluctuation of the approach channel. Therefore, during the operation of the ship lock, the first fluctuation is the main focus.

[0125] According to the main frequency of the main fluctuations and the corresponding wave height data, the wave process caused by the water filling of the ship lock is fitted with a formula. Substituting into the wave formula, the formula of the main fluctuation of the pilot channel is obtained:

[0126]

[0127] like Figure 5 As shown, the formula results of the main fluctuations of the pilot channel and the simulation results of the fluctuation values ​​show that the main fluctuation curve of the single-line ship lock obtained by the method of the present invention is highly consistent with the simulation results.

[0128] The operating conditions of the double-line ship locks with simultaneous operation and staggered operation of Δt=5min, 10min, 15min, 20min, 30min, and 40min were calculated, and the fluctuation amplitudes under different staggered operation modes were compared. Figure 6 It can be seen that when Δt=10min and 30min, the maximum wave height is significantly weakened; when Δt=20min and 40min, the maximum wave height is larger than that of other staggered working conditions, and the amplitude change of superimposed waves is closely related to the staggered time Δt.

[0129] Table 1 Fluctuation data statistics of the gate head location

[0130]

[0131] According to the superposition of sine functions, the lock is staggered. When , the wave height of the pilot channel is superimposed as:

[0132]

[0133] When 2 =2nT, , the fluctuation amplitude is superimposed; when 2 = hour, , fluctuation reduction (n is an integer).

[0134] From this we can see that satisfy When the time difference is too large, the fluctuation amplitudes generated by the time difference are completely superimposed, and the fluctuation amplitude is relatively large. satisfy , the fluctuation amplitude caused by staggered operation is completely offset, and the fluctuation amplitude is small. Therefore, in order to ensure that the fluctuation amplitude generated by the double-line ship lock filling process is small, try to choose staggered operation time when staggered operation is performed. , avoid choosing the wrong running time is an integer multiple of the fluctuation period.

[0135] Example 3

[0136] like Figure 7 As shown, a dispatching system for multi-line ship lock filling operation based on Fourier fast transform is provided, the multi-line ship lock includes multiple single-line ship locks, and the dispatching system includes:

[0137] The fluctuation curve acquisition module is used to respectively obtain the main fluctuation curve of the pilot channel when each single-line ship lock is filled with water;

[0138] The superposition scheduling module is used to superimpose the main fluctuation curves of all single-line ship locks according to the superposition law of waveforms, so as to obtain the staggered operation time of multiple single-line ship locks when the total fluctuation of the pilot channel is minimum, and to schedule the water filling operation of the multi-line ship locks according to the obtained staggered operation time;

[0139] Wherein, the fluctuation curve acquisition module includes a model building unit, a first fluctuation data acquisition unit, a second fluctuation data acquisition unit and a fluctuation curve generation unit;

[0140] A model building unit is used to establish a two-dimensional fluctuation mathematical model of the pilot channel based on the current single-line ship lock water filling operation data, the relationship between the pilot channel water level and the water filling flow, the pilot channel topographic data, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock is filled with water;

[0141] The first fluctuation data acquisition unit is used to obtain the water surface fluctuation data of the pilot channel when the single-line ship lock is filled with water at the initial water level of the pilot channel by using the established two-dimensional fluctuation mathematical model;

[0142] The second fluctuation data acquisition unit is used to establish a Fourier rapid change analysis model based on the water surface fluctuation data of the pilot channel, and obtain the main fluctuations of the pilot channel when the current single-line ship lock is filled with water under the initial water level of the pilot channel, as well as the frequency and fluctuation data of the main fluctuations;

[0143] The fluctuation curve generating unit is used to obtain the main fluctuation curve of the pilot channel when the single-line ship lock is currently filled with water according to the frequency and fluctuation data of the main fluctuation of the pilot channel.

[0144] For more specific processes of the above method, please refer to the corresponding contents disclosed in the aforementioned embodiments, which will not be repeated here.

[0145] Example 4

[0146] A computer device comprises a processor and a memory; wherein, when the processor executes a computer program stored in the memory, the steps of the scheduling method for multi-line ship lock filling operation based on Fourier fast transform of embodiment 1 are implemented.

[0147] For more specific processes of the above method, please refer to the corresponding contents disclosed in the aforementioned embodiments, which will not be repeated here.

[0148] Example 5

[0149] A computer-readable storage medium is used to store a computer program; when the computer program is executed by a processor, the steps of the scheduling method for multi-line ship lock filling operation based on Fourier fast transform in embodiment 1 are implemented.

[0150] For more specific processes of the above method, please refer to the corresponding contents disclosed in the aforementioned embodiments, which will not be repeated here.

[0151] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0152] Those skilled in the art can clearly understand that the technology in the embodiments of the present invention can be implemented by means of software plus a necessary general hardware platform. Based on this understanding, the technical solution in the embodiments of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product, which can be stored in a storage medium such as ROM / RAM, a disk, an optical disk, etc., and includes a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods of various embodiments of the present invention or certain parts of the embodiments.

[0153] The above are only preferred embodiments of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions under the concept of the present invention belong to the protection scope of the present invention. It should be pointed out that for ordinary technicians in this technical field, some improvements and modifications without departing from the principle of the present invention should be regarded as the protection scope of the present invention.

Claims

1. A dispatching method for multi-line ship lock filling operation based on Fourier fast transform, wherein the multi-line ship lock comprises a plurality of single-line ship locks, characterized in that: The scheduling method comprises: Obtain the main fluctuation curve of the navigation channel when each single-line ship lock is filled with water; According to the superposition law of waveforms, the main fluctuation curves of all single-line ship locks are superimposed to obtain the staggered operation time of multiple single-line ship locks when the total fluctuation of the pilot channel is minimum, and the water filling operation of the multi-line ship locks is scheduled according to the obtained staggered operation time; Among them, the main fluctuation curves of the pilot channel when the single-line ship lock is filled with water include: Based on the current single-line ship lock water filling operation data, the relationship between the pilot channel water level and the water filling flow, the pilot channel topographic data, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock water filling, a two-dimensional fluctuation mathematical model of the pilot channel is established; the current single-line ship lock water filling operation data includes the current ship lock upstream and downstream water heads, the ship lock valve opening operation time, the ship lock valve opening, the ship lock upstream water level and the lock chamber water level; the relationship between the pilot channel water level and the water filling flow is designed based on the ship lock valve opening, the ship lock upstream water level and the lock chamber water level and verified according to the measured data; the pilot channel topographic data includes the elevation and cross-sectional shape of the pilot channel; the two-dimensional fluctuation mathematical model is constructed based on the river base flow and non-constant flow; for the river base flow, the pilot channel's own discharge flow and its corresponding pilot channel water level are the initial boundaries; for the non-constant flow, when the ship lock is filled with water, the pilot channel flow change process is the boundary of the non-constant flow; The established two-dimensional wave mathematical model is used to obtain the water surface wave data of the pilot channel when the single-line ship lock is filled with water at the initial water level of the pilot channel. Based on the water surface fluctuation data of the approach channel when the current single-line ship lock is filled with water at the initial water level of the approach channel, a Fourier rapid change analysis model is established, and the main fluctuations of the approach channel when the current single-line ship lock is filled with water at the initial water level of the approach channel, as well as the frequency and fluctuation data of the main fluctuations are obtained, including: Under the initial water level of the pilot channel, the water surface fluctuation data of the pilot channel when the current single-line ship lock is filled with water is standardized to obtain a water level fluctuation data set with equal intervals; The water level fluctuation data set is transformed using Fourier fast transform to obtain a set of complex number arrays, each of which corresponds to a frequency and a modulus. The fluctuation corresponding to the complex number whose modulus value in the complex number array meets the requirement is taken as the main fluctuation of the pilot channel when the current single-line ship lock is filled with water, the fluctuation data corresponding to the complex number whose modulus value meets the requirement is the fluctuation data of the main fluctuation of the pilot channel, and the corresponding frequency is the frequency of the main fluctuation of the pilot channel; According to the frequency and fluctuation data of the main fluctuation of the pilot channel, the main fluctuation curve of the pilot channel when the single-line ship lock is filled with water is obtained by using the fluctuation formula; The fluctuation formula is: ； Where Z(t) is the water level of the pilot channel, m is the number of major fluctuations in the pilot channel, A i is the modulus value of the main fluctuation; N is the number of data in the water level fluctuation dataset, is the frequency of the main fluctuations, t For time, Z is the initial water level of the pilot channel.

2. The dispatching method for multi-line ship lock filling operation based on Fourier fast transform according to claim 1 is characterized in that: In the method of using the fluctuations corresponding to the complex numbers whose modulus values ​​in the complex number array meet the requirements as the main fluctuations of the pilot channel when the current single-line ship lock is filled with water, the fluctuation data corresponding to the complex numbers whose modulus values ​​meet the requirements are the fluctuation data of the main fluctuations of the pilot channel, and the corresponding frequency is the frequency of the main fluctuations of the pilot channel, The main fluctuation of the pilot channel when the current single-line ship lock is filled with water is the fluctuation corresponding to the complex number with the largest modulus value in the complex number array; the frequency of the main fluctuation of the pilot channel for: ； in, is the index corresponding to the complex number with the largest modulus value in the complex array, N is the number of data in the water level fluctuation dataset, Δt is the sampling interval of the water level fluctuation dataset.

3. The dispatching method for multi-line ship lock filling operation based on Fourier fast transform according to claim 1 is characterized in that: A two-dimensional fluctuation mathematical model of the pilot channel is established based on the operating data of the current single-line ship lock filling, the relationship between the water level and the filling flow of the pilot channel, the terrain data of the pilot channel, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock is filled with water. After the two-dimensional fluctuation mathematical model is established, the calculation results of the two-dimensional fluctuation mathematical model are verified using the water surface fluctuation monitoring data of the pilot channel after the single-line ship lock is filled with water.

4. The dispatching method for multi-line ship lock filling operation based on Fourier fast transform according to claim 1 is characterized in that: In the method of using the established two-dimensional fluctuation mathematical model to obtain the water surface fluctuation data of the pilot channel when the current single-line ship lock is filled with water at the initial water level of the pilot channel, only the water surface fluctuation data of the current single-line ship lock within 4 hours of filling operation at the initial water level of the pilot channel is obtained.

5. A dispatching system for multi-line ship lock filling operation based on Fourier fast transform, wherein the multi-line ship lock includes a plurality of single-line ship locks, characterized in that: The dispatching system comprises: The fluctuation curve acquisition module is used to respectively obtain the main fluctuation curve of the pilot channel when each single-line ship lock is filled with water; The superposition scheduling module is used to superimpose the main fluctuation curves of all single-line ship locks according to the superposition law of waveforms, so as to obtain the staggered operation time of multiple single-line ship locks when the total fluctuation of the pilot channel is minimum, and to schedule the water filling operation of the multi-line ship locks according to the obtained staggered operation time; Wherein, the fluctuation curve acquisition module includes a model building unit, a first fluctuation data acquisition unit, a second fluctuation data acquisition unit and a fluctuation curve generation unit; The model building unit is used to establish a two-dimensional fluctuation mathematical model of the pilot channel based on the current single-line ship lock water filling operation data, the relationship between the pilot channel water level and the water filling flow, the pilot channel topographic data, and the water surface fluctuation monitoring data of the pilot channel after the current single-line ship lock is filled with water; the current single-line ship lock water filling operation data includes the current ship lock upstream and downstream water heads, the ship lock valve opening operation time, the ship lock valve opening, the ship lock upstream water level and the lock chamber water level; the relationship between the pilot channel water level and the water filling flow is designed based on the ship lock valve opening, the ship lock upstream water level and the lock chamber water level and verified according to the measured data; the pilot channel topographic data includes the elevation and cross-sectional shape of the pilot channel; the two-dimensional fluctuation mathematical model is constructed based on the river base flow and non-constant flow; for the river base flow, the pilot channel's own discharge flow and its corresponding pilot channel water level are the initial boundaries; for the non-constant flow, when the ship lock is filled with water, the pilot channel flow change process is the boundary of the non-constant flow; The first fluctuation data acquisition unit is used to obtain the water surface fluctuation data of the pilot channel when the single-line ship lock is filled with water at the initial water level of the pilot channel by using the established two-dimensional fluctuation mathematical model; The second fluctuation data acquisition unit is used to establish a Fourier fast change analysis model based on the water surface fluctuation data of the pilot channel when the single-line ship lock is currently filled with water at the initial water level of the pilot channel, and obtain the main fluctuations of the pilot channel when the single-line ship lock is currently filled with water at the initial water level of the pilot channel, as well as the frequency and fluctuation data of the main fluctuations; specifically including: Under the initial water level of the pilot channel, the water surface fluctuation data of the pilot channel when the current single-line ship lock is filled with water is standardized to obtain a water level fluctuation data set with equal intervals; The water level fluctuation data set is transformed using Fourier fast transform to obtain a set of complex number arrays, each of which corresponds to a frequency and a modulus. The fluctuation corresponding to the complex number whose modulus value in the complex number array meets the requirement is taken as the main fluctuation of the pilot channel when the current single-line ship lock is filled with water, the fluctuation data corresponding to the complex number whose modulus value meets the requirement is the fluctuation data of the main fluctuation of the pilot channel, and the corresponding frequency is the frequency of the main fluctuation of the pilot channel; The fluctuation curve generating unit is used to obtain the main fluctuation curve of the pilot channel when the single-line ship lock is currently filled with water according to the frequency and fluctuation data of the main fluctuation of the pilot channel and the fluctuation formula; the fluctuation formula is: ； Where Z(t) is the water level of the pilot channel, m is the number of major fluctuations in the pilot channel, A i is the modulus value of the main fluctuation; N is the number of data in the water level fluctuation dataset, is the frequency of the main fluctuations, t For time, Z is the initial water level of the pilot channel.

6. A computer device, characterized in that: It comprises a processor and a memory; wherein, when the processor executes the computer program stored in the memory, it implements the steps of the scheduling method for multi-line ship lock filling operation based on Fourier fast transform as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: Used to store computer programs; when the computer programs are executed by the processor, the steps of the scheduling method for multi-line ship lock filling operation based on Fourier fast transform as described in any one of claims 1-4 are implemented.

Images