Modeling method of Saint-Venant equation set containing gate based on model order reduction technology

Abstract

The invention discloses a gate-containing Saint-Venant equation set modeling method based on a model order reduction technology, and belongs to the technical field of hydraulic engineering, and the gate-containing Saint-Venant equation set modeling method comprises the following steps: S1, constructing a gate-containing original Saint-Venant equation set; s2, discretizing the constructed original Saint-Venant equation set to obtain a nonlinear function in a high-dimensional system; s3, based on a simulation purpose, performing system simulation by using the original Saint-Venant equation set, and collecting state variables to construct a snapshot matrix; and S4, based on the nonlinear function in the high-dimensional system, simplifying the snapshot matrix by using a model order reduction technology, and further constructing a low-dimensional order reduction model corresponding to the original Saint-Venant equation set, thereby completing the construction of the Saint-Venant equation set containing the gate. According to the method, the model is subjected to order reduction processing, so that the simulation precision is ensured, the calculation efficiency is greatly improved, and the requirements of real-time simulation and control are met; the method has a wide application prospect and can be used in the fields of river scheduling, water conservancy project optimization, flood control early warning and the like.

Description

Technical Field

[0001] The invention belongs to the technical field of water conservancy engineering, and in particular relates to a modeling method of a Saint-Venant equation group including a gate based on a model order reduction technology. Background Art

[0002] The Saint-Venant equations are a set of partial differential equations used to describe the motion of one-dimensional non-steady open channel water flow. They are widely used in water flow simulation and prediction in water conservancy projects such as rivers, reservoirs and irrigation canals. In practical applications, multiple gates are often arranged on the river to control the water level and flow. The presence of gates has a significant impact on the dynamic characteristics of the water flow, which increases the complexity of the mathematical model. However, due to the introduction of nonlinear and complex boundary conditions, the numerical solution process of the Saint-Venant equations containing gates has a large amount of calculation and low simulation efficiency, which makes it difficult to meet the needs of real-time simulation and control. Therefore, the existing numerical simulation methods of the Saint-Venant equations containing gates have the problems of low computational efficiency and difficulty in real-time application. Summary of the invention

[0003] In view of the above-mentioned deficiencies in the prior art, the modeling method of the Saint-Venant equations group including gates based on the model reduction technology provided by the present invention solves the problem that in the traditional method, when the Saint-Venant equations group including gates are used for simulation, the numerical solution process has a large amount of calculation, the simulation efficiency is low, and it is difficult to meet the requirements of real-time applications.

[0004] In order to achieve the above-mentioned invention object, the technical solution adopted by the present invention is: a modeling method of the Saint-Venant equations group including gates based on model reduction technology, which is used to realize water flow simulation and prediction simulation of water conservancy projects, including the following steps: S1. Construct the original Saint-Venant equations including the gate; S2, discretize the constructed original Saint-Venant equations to obtain nonlinear functions in the high-dimensional system; S3. Based on the simulation purpose, the original Saint-Venant equations are used to simulate the system, collect state variables and construct a snapshot matrix; S4. Based on the nonlinear functions in the high-dimensional system, the snapshot matrix is ​​simplified using model reduction technology, and then a low-dimensional reduced-order model corresponding to the original Saint-Venant equations is constructed to complete the construction of the Saint-Venant equations including the gate.

[0005] Furthermore, in step S1, the original Saint-Venant equations including the gate include the continuity equation and the momentum equation, and according to the gate position, the flow relationship of the gate is used as the boundary condition of the original Saint-Venant equations.

[0006] Furthermore, in the original Saint-Venant equations, the continuity equation is:

[0007] The momentum equation is:

[0008] The boundary conditions are:

[0009] In the formula, A represents the cross-sectional area of ​​water flow, Q Indicates flow rate, t Indicates time, represents the length of the river, h Indicates water depth. g represents the acceleration due to gravity, Indicates the bottom slope, represents the resistance slope, represents the flow rate through the gate, represents the flow coefficient, Indicates the gate opening, Indicates the gate width, and Respectively represent the water level before and after the gate.

[0010] Furthermore, in step S2, the discretization processing of the Saint-Venant equations includes spatial discretization and temporal discretization; In the control volume i At , the continuity equation is discretized as:

[0011] The momentum equation is discretized as:

[0012] In the formula, Indicates i Units at time step n The depth of water, Indicates i The flow flux at the right boundary of the unit, Indicates i The flow flux at the left boundary of the unit, Indicates i Units at time step n of traffic, Indicates i The momentum flux at the right boundary of the unit, Indicates i The momentum flux at the left boundary of the unit, Indicates i Units at time step n The friction slope, represents the time step of time discretization, Represents the step size of discretization of river channel length.

[0013] Further, i The momentum flux at the right boundary of the unit The calculation formula is:

[0014] In the formula, Indicates Units at time step of traffic, Indicates Units at time step of traffic, Indicates Units at time step The depth of water, Indicates Units at time step The depth of water, represents the step size of discretization of river channel length, Represents the acceleration due to gravity.

[0015] Further, i Units at time step n Friction slope The calculation formula is:

[0016] In the formula, represents the Manning roughness coefficient, Indicates i Units at time step n of traffic, represents the cross-sectional area of ​​water flow, Represents the hydraulic parameter wetted perimeter.

[0017] Furthermore, in step S2, when the Saint-Venant equations are discretized, when the gate position is When the gate boundary conditions are met, the momentum flux Replaced with the corresponding gate flow rate, which is expressed as:

[0018] In the formula, Indicates j The gate flow rate of each gate is represents the flow coefficient, Indicates j The gates are in the time step n The opening degree, Indicates jThe width of the gate, Indicates i Unit j The water depth of the gate at time step n, Indicates the control volume of the gate, the subscript j Indicates the number of gates, .

[0019] Furthermore, in step S2, the system state vector after discretization is It is expressed as:

[0020] In the formula, M Quantity represents the water depth of each control volume, and Quantity represents the flow rate of each control volume; Dynamical equations of high-dimensional systems:

[0021] In the formula, Represents a high-dimensional system dimensional nonlinear function, which can be expressed as:

[0022] In the formula, the symbol Indicates that it is located at M Units and M -1 momentum flux between cells, Respectively represent the first unit... M The water depth of the unit, Respectively represent the first unit to the M The friction slope of each unit, represents the acceleration due to gravity, Indicates the bottom slope; Among them, M The component is the water depth change rate of the high-dimensional system, which is expressed as:

[0023] back M The component is the flow rate change rate of the high-dimensional system, which is expressed as:

[0024] The momentum flux is expressed as:

[0025] In the formula, Indicates gate overflow.

[0026] Furthermore, the step S3 includes the following sub-steps: S31. Based on the simulation purpose, at the selected sampling time or working condition, the original Saint-Venant equation is used to simulate the system and obtain the high-dimensional state vector of the system ;in, represents the sampling time; S32, the obtained high-dimensional state vector is used as the column vector of the snapshot matrix S in turn, and then the snapshot matrix S is constructed. : .

[0027] Furthermore, the step S4 includes the following sub-steps: S41. Snapshot Matrix Perform singular value decomposition to obtain the left singular vector matrix , singular value matrix Σ and right singular vector matrix V, take the largest k Singular values ​​to construct the reduced-order basis matrix ;

[0028] In the formula, represents the reduced-order basis matrix , represents a set of real matrices consisting of matrices of dimension (2M,k), Indicates the number of cross-sectional elements; S42, according to the reduced-order basis matrix , construct the reduced-order state variables ;

[0029] In the formula, represents the system state vector after discretization of the original Saint-Venant equations, represents the transpose operator, Indicated by k The vector set consisting of dimensional real vectors; S43. Based on the constructed reduced-order state variables, the high-dimensional system obtained by discretizing the original Saint-Venant equations is projected onto the reduced-order basis matrix to construct a low-dimensional reduced-order model;

[0030] In the formula, represents the state variable value of the low-dimensional reduced-order model, represents the nonlinear function in the high-dimensional system obtained by discretizing the original Saint-Venant equations, Represents the transpose operator.

[0031] The beneficial effects of the present invention are: (1) In the method of the present invention, by reasonably selecting the reduced-order basis when constructing the reduced-order model, the simulation results maintain good consistency with the high-precision model, and the error is within an acceptable range.

[0032] (3) By reducing the order of the model, the degree of freedom of the reduced-order model is much smaller than that of the original model. While ensuring the simulation accuracy, the computational efficiency is greatly improved, meeting the needs of real-time simulation and control.

[0033] (3) The method of the present invention has broad application prospects and can be used in the fields of river channel scheduling, water conservancy project optimization and flood prevention and early warning. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] Figure 1 A flow chart of the modeling method of the Saint-Venant equations group including gates based on the model reduction technology provided by the present invention. DETAILED DESCRIPTION

[0035] The specific implementation modes of the present invention are described below so that those skilled in the art can understand the present invention. However, it should be clear that the present invention is not limited to the scope of the specific implementation modes. For those of ordinary skill in the art, as long as various changes are within the spirit and scope of the present invention as defined and determined by the attached claims, these changes are obvious, and all inventions and creations utilizing the concept of the present invention are protected.

[0036] like Figure 1 As shown, an embodiment of the present invention provides a modeling method of Saint-Venant equations including gates based on model reduction technology, which is used to realize water flow simulation and prediction simulation of water conservancy projects, including the following steps: S1. Construct the original Saint-Venant equations including the gate; S2, discretize the constructed original Saint-Venant equations to obtain nonlinear functions in the high-dimensional system; S3. Based on the simulation purpose, the original Saint-Venant equations are used to simulate the system, collect state variables and construct a snapshot matrix; S4. Based on the nonlinear functions in the high-dimensional system, the snapshot matrix is ​​simplified using model reduction technology, and then a low-dimensional reduced-order model corresponding to the original Saint-Venant equations is constructed to complete the construction of the Saint-Venant equations including the gate.

[0037] In step S1 of the embodiment of the present invention, considering that N gates are arranged on the river channel, the traditional Saint-Venant equations need to be expanded; therefore, in this embodiment, the original Saint-Venant equations containing the gates include the continuity equation and the momentum equation, and according to the gate position, the flow relationship of the gate is used as the boundary condition of the original Saint-Venant equations; The continuity equation is:

[0038] The momentum equation is:

[0039] The boundary conditions are:

[0040] In the formula, A represents the cross-sectional area of ​​water flow, Q Indicates flow rate, t Indicates time, represents the length of the river, h Indicates water depth. g represents the acceleration due to gravity, Indicates the bottom slope, represents the resistance slope, represents the flow rate through the gate, represents the flow coefficient, Indicates the gate opening, Indicates the gate width, and Respectively represent the water level before and after the gate.

[0041] In step S2 of the embodiment of the present invention, the discretization processing of the Saint-Venant equations includes spatial discretization and temporal discretization; Among them, spatial discretization means that the river length L Divide into M equidistant space units, each unit length is:

[0042] The spatial node positions are:

[0043] Each control volume i lie in The range is ,in:

[0044] Time discretization refers to setting the time step size , the time nodes are: .

[0045] Based on this, in this embodiment, in the control volume i At , the continuity equation is discretized as:

[0046] The momentum equation is discretized as:

[0047] In the formula, Indicates i Units at time step n The depth of water, Indicates i The flow flux at the right boundary of the unit, Indicates i The flow flux at the left boundary of the unit, Indicates i Units at time step n of traffic, Indicates i The momentum flux at the right boundary of the unit, Indicates i The momentum flux at the left boundary of the unit, Indicates i Units at time step n The friction slope, represents the time step of time discretization, Represents the step size of discretization of river channel length.

[0048] The central difference method is used to calculate the i The momentum flux at the right boundary of the unit for:

[0049] In the formula, Indicates Units at time step of traffic, Indicates Units at time step of traffic, Indicates Units at time step The depth of water, Indicates Units at time step The depth of water, represents the step size of discretization of river channel length, Represents the acceleration due to gravity.

[0050] Among them, the Manning formula is used to calculate the i Units at time step n Friction slope The is:

[0051] In the formula, represents the Manning roughness coefficient, Indicates i Units at time stepn of traffic, represents the cross-sectional area of ​​water flow, It represents the hydraulic parameter wetted perimeter, which is obtained by measuring the specific cross section.

[0052] In step S2 of the embodiment of the present invention, when the Saint-Venant equations are discretized, when the gate position is When the gate boundary conditions are met, the momentum flux Replaced with the corresponding gate flow rate, which is expressed as:

[0053] In the formula, Indicates j The gate flow rate of each gate is represents the flow coefficient, Indicates j The gates are in the time step n The opening degree, Indicates j The width of the gate, Indicates i Unit j The water depth of the gate at time step n, Indicates the control volume of the gate, the subscript j Indicates the number of gates, .

[0054] In step S2 of the embodiment of the present invention, the system state vector after the discretization of the original Saint-Venant equations is It is expressed as:

[0055] In the formula, M Quantity represents the water depth of each control volume, and Quantity represents the flow rate of each control volume; Dynamical equations of high-dimensional systems:

[0056] In the formula, Represents a high-dimensional system dimensional nonlinear function, which can be expressed as:

[0057] In the formula, the symbol Indicates that it is located at M Units and M -1 momentum flux between cells, Respectively represent the first unit... M The water depth of the unit, Respectively represent the first unit to the M The friction slope of each unit, represents the acceleration due to gravity, Indicates the bottom slope; Among them, M The component is the water depth change rate of the high-dimensional system, which is expressed as:

[0058] back M The component is the flow rate change rate of the high-dimensional system, which is expressed as:

[0059] The momentum flux is expressed as:

[0060] In the formula, Indicates gate overflow.

[0061] Step S3 of the embodiment of the present invention includes the following sub-steps: S31. Based on the simulation purpose, at the selected sampling time or working condition, the original Saint-Venant equation is used to simulate the system and obtain the high-dimensional state vector of the system ;in, represents the sampling time; S32, the obtained high-dimensional state vector is used as the column vector of the snapshot matrix S in turn, and then the snapshot matrix S is constructed. :

[0062] In step S4 of the embodiment of the present invention, in order to improve the simulation efficiency, the model reduction technology is used in this embodiment to simplify the expanded Saint-Venant equations. Based on this, step S4 of this embodiment includes the following sub-steps: S41. Snapshot Matrix Perform singular value decomposition to obtain the left singular vector matrix , singular value matrix Σ and right singular vector matrix V, take the largest k Singular values ​​to construct the reduced-order basis matrix ;

[0063] In the formula, represents the reduced-order basis matrix , represents a set of real matrices consisting of matrices of dimension (2M,k), Indicates the number of cross-sectional elements; S42, according to the reduced-order basis matrix , construct the reduced-order state variables ;

[0064] In the formula, represents the system state vector after discretization of the original Saint-Venant equations, represents the transpose operator, Indicated by k The vector set consisting of dimensional real vectors; S43. Based on the constructed reduced-order state variables, the high-dimensional system obtained by discretizing the original Saint-Venant equations is projected onto the reduced-order basis matrix to construct a low-dimensional reduced-order model;

[0065] In the formula, represents the state variable value of the low-dimensional reduced-order model, represents the nonlinear function in the high-dimensional system obtained by discretizing the original Saint-Venant equations, Represents the transpose operator.

[0066] In this embodiment, the state variables in the reduced-order system are a It is the core variable of the low-dimensional reduced-order model, representing the projection of the original high-dimensional system state on the reduced-order basis. is the state vector of the low-dimensional system, representing the high-dimensional state In the reduced base The projection on Each component ( ) represents a high-dimensional state In the In order to calculate the reduced-order state variables, , we need to convert the high-dimensional state Projection to a reduced-order basis Then, the above-mentioned low-dimensional reduced-order model was established.

[0067] The present invention uses specific embodiments to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present invention.

[0068] Those skilled in the art will appreciate that the embodiments described herein are intended to help readers understand the principles of the present invention, and should be understood that the protection scope of the present invention is not limited to such specific statements and embodiments. Those skilled in the art can make various other specific variations and combinations that do not deviate from the essence of the present invention based on the technical revelations disclosed by the present invention, and these variations and combinations are still within the protection scope of the present invention.

Claims

1. A modeling method for the Saint-Venant equations including gates based on model reduction technology, characterized in that: The method is used to realize water flow simulation and prediction simulation of water conservancy projects, including the following steps: S1. Construct the original Saint-Venant equations including the gate; S2, discretize the constructed original Saint-Venant equations to obtain nonlinear functions in the high-dimensional system; S3. Based on the simulation purpose, the original Saint-Venant equations are used to simulate the system, collect state variables and construct a snapshot matrix; S4. Based on the nonlinear functions in the high-dimensional system, the snapshot matrix is ​​simplified using model reduction technology, and then a low-dimensional reduced-order model corresponding to the original Saint-Venant equations is constructed to complete the construction of the Saint-Venant equations including the gate.

2. The method for modeling the gated Saint-Venant equations based on model order reduction technology according to claim 1, characterized in that: In step S1, the original Saint-Venant equations including the gate include the continuity equation and the momentum equation, and according to the gate position, the flow relationship of the gate is used as the boundary condition of the original Saint-Venant equations.

3. The Saint-Venant equations modeling method including gates based on model order reduction technology according to claim 2, characterized in that: In the original Saint-Venant equations, the continuity equation is: The momentum equation is: The boundary conditions are: In the formula, A represents the cross-sectional area of ​​water flow, Q Indicates flow rate, t Indicates time, represents the length of the river, h Indicates water depth. g represents the acceleration due to gravity, Indicates the bottom slope, represents the resistance slope, represents the flow rate through the gate, represents the flow coefficient, Indicates the gate opening, Indicates the gate width, and Respectively represent the water level before and after the gate.

4. The Saint-Venant equations modeling method including gates based on model order reduction technology according to claim 3, characterized in that: In step S2, the discretization process of the Saint-Venant equations includes spatial discretization and temporal discretization; In the control volume i At , the continuity equation is discretized as: The momentum equation is discretized as: In the formula, Indicates i Units at time step n The depth of water, Indicates i The flow flux at the right boundary of the unit, Indicates i The flow flux at the left boundary of the unit, Indicates i Units at time step n of traffic, Indicates i The momentum flux at the right boundary of the unit, Indicates i The momentum flux at the left boundary of the unit, Indicates i Units at time step n The friction slope, represents the time step of time discretization, Represents the step size of discretization of river channel length.

5. The Saint-Venant equations modeling method including gates based on model order reduction technology according to claim 4, characterized in that: No. i The momentum flux at the right boundary of the unit The calculation formula is: In the formula, Indicates Units at time step of traffic, Indicates Units at time step of traffic, Indicates Units at time step The depth of water, Indicates Units at time step The depth of water, represents the step size of discretization of river channel length, Represents the acceleration due to gravity.

6. The Saint-Venant equations modeling method including gates based on model order reduction technology according to claim 4, characterized in that: No. i Units at time step n Friction slope The calculation formula is: In the formula, represents the Manning roughness coefficient, Indicates i Units at time step n of traffic, represents the cross-sectional area of ​​water flow, Represents the hydraulic parameter wetted perimeter.

7. The method for modeling Saint-Venant equations including gates based on model order reduction technology according to claim 4, characterized in that: In step S2, when the Saint-Venant equations are discretized, when the gate position is When the gate boundary conditions are met, the momentum flux Replaced with the corresponding gate flow rate, which is expressed as: In the formula, Indicates j The gate flow rate of each gate is represents the flow coefficient, Indicates j The gates are in the time step n The opening degree, Indicates j The width of the gate, Indicates i Unit j The water depth of the gate at time step n, Indicates the control volume of the gate, the subscript j Indicates the number of gates, .

8. The method for modeling Saint-Venant equations including gates based on model order reduction technology according to claim 4, characterized in that: In step S2, the system state vector after discretization processing It is expressed as: In the formula, M Quantity represents the water depth of each control volume, and Quantity represents the flow rate of each control volume; Dynamical equations of high-dimensional systems: In the formula, Represents a high-dimensional system dimensional nonlinear function, which can be expressed as: In the formula, the symbol Indicates that it is located at M Unit and M -1 momentum flux between cells, Respectively represent the first unit... M The water depth of the unit, Respectively represent the first unit to the M The friction slope of each unit, represents the acceleration due to gravity, Indicates the bottom slope; Among them, M The component is the water depth change rate of the high-dimensional system, which is expressed as: back M The component is the flow rate change rate of the high-dimensional system, which is expressed as: The momentum flux is expressed as: In the formula, Indicates gate overflow.

9. The method for modeling Saint-Venant equations including gates based on model order reduction technology according to claim 1, characterized in that: The step S3 comprises the following sub-steps: S31. Based on the simulation purpose, at the selected sampling time or working condition, the original Saint-Venant equation is used to simulate the system and obtain the high-dimensional state vector of the system ;in, represents the sampling time; S32, the obtained high-dimensional state vector is used as the column vector of the snapshot matrix S in turn, and then the snapshot matrix S is constructed. : .

10. The method for modeling Saint-Venant equations including gates based on model order reduction technology according to claim 6, characterized in that: The step S4 comprises the following sub-steps: S41. Snapshot Matrix Perform singular value decomposition to obtain the left singular vector matrix , singular value matrix Σ and right singular vector matrix V, take the largest k Singular values ​​to construct the reduced-order basis matrix ; In the formula, represents the reduced-order basis matrix , represents a set of real matrices consisting of matrices of dimension (2M,k), Indicates the number of cross-sectional elements; S42, according to the reduced-order basis matrix , construct the reduced-order state variables ; In the formula, represents the system state vector after discretization of the original Saint-Venant equations, represents the transpose operator, Indicated by k The vector set consisting of dimensional real vectors; S43. Based on the constructed reduced-order state variables, the high-dimensional system obtained by discretizing the original Saint-Venant equations is projected onto the reduced-order basis matrix to construct a low-dimensional reduced-order model; In the formula, represents the state variable value of the low-dimensional reduced-order model, represents the nonlinear function in the high-dimensional system obtained by discretizing the original Saint-Venant equations, Represents the transpose operator.

Images