Modeling method and device for one-dimensional water power model of irrigation area in case of insufficient cross-section data

Abstract

The application provides a method and device for modeling a one-dimensional water power model of an irrigation area under the condition of insufficient cross-section data, which comprises the following steps: obtaining a hydraulic property interpolation table of an upstream cross-section and a hydraulic property interpolation table of a downstream cross-section corresponding to a target cross-section to be interpolated in a ditch of the irrigation area; determining an interpolation coefficient based on the relative position relationship between the target cross-section and the upstream cross-section and the downstream cross-section; in the hydraulic parameter space, taking the hydraulic property interpolation table of the upstream cross-section and the hydraulic property interpolation table of the downstream cross-section as operation objects, and performing interpolation by using the interpolation coefficient to obtain a hydraulic property interpolation table of the target cross-section; and constructing a one-dimensional water power model of the irrigation area based on the hydraulic property interpolation table of the target cross-section. In the application, the hydraulic property interpolation table is directly taken as an operation object under the condition of insufficient cross-section data, and interpolation is performed in the hydraulic parameter space, so that the accuracy of the obtained hydraulic property interpolation table of the target cross-section is improved.

Description

Technical Field

[0001] This invention relates to the field of hydrodynamic analysis technology, and in particular to a method and apparatus for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient. Background Technology

[0002] One-dimensional unsteady flow dynamic models based on the Saint-Venant equations are important tools for simulating and analyzing water conveyance and distribution processes in irrigation districts. The Saint-Venant equations include continuity and momentum equations. When using finite difference or finite volume methods for numerical solutions, the functional relationships between hydraulic characteristics such as cross-sectional area, water surface width, and hydraulic radius and water level values ​​at each cross-section of the spatial computational grid are required to support the calculation of discrete terms. In practical modeling, if the spacing between adjacent computational cross-sections is too large, the discretization error of the spatial derivative term will increase, potentially leading to distorted simulation results. Generally, the spacing between computational cross-sections is required to meet numerical stability conditions. This means that the number of computational cross-sections is often much larger than the number of measured cross-sections, necessitating some method to supplement the cross-sectional data at missing locations.

[0003] Currently, the automatic cross-section interpolation method is as follows: between two known cross-sections, a chord (including primary and secondary chords) connecting the corresponding geometric features is constructed, and linear interpolation based on distance weights is performed at both ends of the chord to generate the geometric shape of the intermediate cross-section, and the hydraulic characteristics at each water level are calculated accordingly.

[0004] This method is inaccurate in determining the hydraulic characteristics of each section of the irrigation canal when cross-sectional data is insufficient. Summary of the Invention

[0005] This invention provides a method and apparatus for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, in order to solve the problem of inaccurate hydraulic characteristics of each cross-section in the irrigation area ditch determined under the condition of insufficient cross-sectional data in the prior art, and to improve the accuracy of the hydraulic characteristics of the obtained cross-section.

[0006] This invention provides a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, comprising the following steps: Obtain the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section corresponding to the target section to be interpolated in the irrigation area's ditches; Based on the relative positional relationship between the target section and the upstream and downstream sections, the interpolation coefficients are determined; In the hydraulic parameter space, the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section are used as the operation objects, and the interpolation coefficients are used to perform interpolation to obtain the hydraulic characteristic interpolation table of the target section. Based on the hydraulic characteristic interpolation table of the target section, a one-dimensional hydrodynamic model of the irrigation area is constructed.

[0007] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient includes, in the hydraulic parameter space, using the hydraulic characteristic interpolation tables of the upstream and downstream cross-sections as the operation objects, and performing interpolation using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target cross-section, including: The water level values ​​at each interpolation point in the hydraulic characteristic interpolation tables of the upstream and downstream sections are converted into relative water depths to obtain reference hydraulic characteristic interpolation tables for the upstream and downstream sections respectively. Interpolation is performed using the interpolation coefficients and the reference hydraulic characteristic interpolation tables for the upstream and downstream sections to obtain the reference hydraulic characteristic interpolation table for the target section. The relative water depths in the reference hydraulic characteristic interpolation table of the target cross section are converted to obtain the hydraulic characteristic interpolation table of the target cross section.

[0008] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient is provided. The reference hydraulic characteristic interpolation table of the target cross-section includes a relative water depth-flow cross-sectional area interpolation table, a relative water depth-water surface width interpolation table, a relative water depth-hydraulic radius interpolation table, and a flow cross-sectional area-relative water depth interpolation table. Interpolation is performed using the interpolation coefficients and the reference hydraulic characteristic interpolation tables for the upstream and downstream sections to obtain the reference hydraulic characteristic interpolation table for the target section, including: Using the interpolation coefficients and the relative water depth-flow cross-sectional area interpolation table between the upstream and downstream sections, the relative water depth-flow cross-sectional area interpolation table of the target section is determined; and using the interpolation coefficients and the relative water depth-hydraulic radius interpolation table between the upstream and downstream sections, the relative water depth-hydraulic radius interpolation table of the target section is determined. Numerical differentiation is performed on each flow cross-sectional area in the relative water depth-flow cross-sectional area interpolation table of the target cross-section to obtain the relative water depth-water surface width interpolation table of the target cross-section. The relative water depth-flow cross-sectional area interpolation table is converted to obtain the flow cross-sectional area-relative water depth interpolation table.

[0009] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient is provided, which converts the water level values ​​at each interpolation point in the hydraulic characteristic interpolation table of the upstream and downstream cross-sections into relative water depths, including: For the hydraulic characteristic interpolation table of the upstream section or the downstream section, perform the following steps: Obtain the lowest water level value from the hydraulic characteristic interpolation table; The difference between the water level value at each interpolation point in the hydraulic characteristic interpolation table and the lowest water level value is taken as the relative water depth at each interpolation point in the hydraulic characteristic interpolation table.

[0010] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient is provided, wherein the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section are within the effective water depth range of the target cross-section, and the method further includes: Determine the first maximum relative water depth in the reference hydraulic characteristic interpolation table of the upstream section and the second maximum relative water depth in the reference hydraulic characteristic interpolation table of the downstream section; The effective water depth range of the target cross section is obtained based on the interpolation coefficient, the first maximum relative water depth, and the second maximum relative water depth.

[0011] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient involves converting the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section to obtain the hydraulic characteristic interpolation table of the target cross-section, including: Obtain the bottom elevation of the target section, which is based on the interpolation coefficient, the bottom elevation of the upstream section, and the bottom elevation of the downstream section; Based on the bottom elevation, the relative water depth of each interpolation point in the reference hydraulic characteristic interpolation table of the target section is converted to obtain the basic hydraulic characteristic interpolation table of the target section; The basic hydraulic characteristic interpolation table is sparsified to obtain the hydraulic characteristic interpolation table of the target section.

[0012] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient includes sparse processing of the basic hydraulic characteristic interpolation table to obtain the hydraulic characteristic interpolation table of the target cross-section, comprising: The basic hydraulic characteristic interpolation table of the target section is divided into segments to obtain several interpolation segments; For each interpolation segment, the interpolation points to be removed are determined based on the linearity of the interpolation segment, and the interpolation points to be removed are processed to obtain the hydraulic characteristic interpolation table of the target section.

[0013] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, before converting the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section to obtain the hydraulic characteristic interpolation table of the target cross-section, the method further includes: The physical rationality of the reference hydraulic characteristic interpolation table is verified, including monotonicity verification, non-negativity verification and / or endpoint continuity verification. If the result of the physical rationality verification is that the verification is passed, the step of converting the relative water depths in the reference hydraulic characteristic interpolation table of the target section to obtain the hydraulic characteristic interpolation table of the target section is executed. If the physical rationality check fails, an error message is generated and the current interpolation process for the target section is terminated.

[0014] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient is provided. This method obtains hydraulic characteristic interpolation tables for the upstream and downstream sections corresponding to the target cross-section to be interpolated in the irrigation area's ditches, including: The topographic data and hydrodynamic interpolation parameters of the upstream section and the downstream section are obtained respectively. Based on the topographic data of the upstream section and the hydrodynamic interpolation parameters, a hydraulic characteristic interpolation table for the upstream section is constructed; and based on the topographic data of the downstream section and the hydrodynamic interpolation parameters, a hydraulic characteristic interpolation table for the downstream section is constructed.

[0015] According to the present invention, a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient is provided, the method further comprising: The data conditions for each ditch in the irrigation area are determined. The data conditions include scattered data with multiple cross-sections, scattered data with two cross-sections (inlet and outlet), design parameters of the ditch, or partial shape parameters of the ditch. For each ditch, if the data condition of the ditch is scattered data with multiple cross-sections or scattered data with two cross-sections (inlet and outlet), the cross-section to be interpolated is taken as the target cross-section, and the two cross-sections with scattered data are taken as the upstream and downstream cross-sections of the target cross-section, respectively.

[0016] This invention also provides a modeling device for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, comprising the following modules: The hydraulic characteristic interpolation table acquisition module is used to acquire the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section corresponding to the target section to be interpolated in the irrigation area's ditches. The interpolation coefficient determination module is used to determine the interpolation coefficients based on the relative positional relationship between the target section and the upstream section and the downstream section. The interpolation module is used to perform interpolation in the hydraulic parameter space, using the hydraulic characteristic interpolation tables of the upstream section and the downstream section as operation objects, and using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target section. The one-dimensional hydrodynamic model construction module is used to construct a one-dimensional hydrodynamic model of the irrigation area based on the hydraulic characteristic interpolation table of the target cross section.

[0017] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the one-dimensional hydrodynamic modeling method for irrigation areas when cross-sectional data is insufficient, as described above.

[0018] The present invention also provides a non-transitory computer-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the one-dimensional hydrodynamic modeling method for irrigation areas when cross-sectional data is insufficient, as described in any of the above embodiments.

[0019] The present invention also provides a computer program product, including a computer program that, when executed by a processor, implements the one-dimensional hydrodynamic modeling method for irrigation areas when cross-sectional data is insufficient, as described above.

[0020] This invention provides a method and apparatus for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient. Considering the significant differences in cross-sectional morphology between irrigation ditches and natural rivers, and the varying cross-sectional characteristics of different water system types, constructing chords (including primary and secondary chords) connecting the corresponding geometric features of the cross-section and performing linear interpolation based on distance weights at both ends of the chords to generate the geometric shape of the intermediate cross-section and determine its hydraulic characteristics would result in inaccurate hydraulic characteristics due to significant differences between the constructed intermediate cross-section's combined shape and the actual shape. Therefore, this application, when cross-sectional data is insufficient, uses the hydraulic characteristic interpolation table of the upstream and downstream cross-sections corresponding to the target cross-section as the direct operation object, performing interpolation in the hydraulic parameter space without considering the actual shape of the target cross-section, thus improving the accuracy of the determined actual hydraulic characteristics of the target cross-section and consequently improving the accuracy of the one-dimensional hydrodynamic model. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0022] Figure 1 This is one of the flowcharts illustrating the method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, as provided by the present invention.

[0023] Figure 2This is the second flowchart illustrating the method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, as provided by this invention.

[0024] Figure 3 This is the third flowchart illustrating the method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, as provided by this invention.

[0025] Figure 4 This is a schematic diagram of the structure of the one-dimensional hydrodynamic modeling device for irrigation areas when cross-sectional data is insufficient, as provided by the present invention.

[0026] Figure 5 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

[0027] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0028] The following is combined with Figures 1 to 3 This invention describes a method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient. Figure 1 This is one of the flowcharts illustrating the method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, as provided in this invention. Figure 1 As shown, the method includes the following: Step 101: Obtain the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section corresponding to the target section to be interpolated in the irrigation area's ditches.

[0029] An irrigation district can be understood as a specific area requiring agricultural irrigation. Its internal water network is complex, including various levels of canals for water distribution and various levels of ditches for drainage. Ditches are a general term for the water system within an irrigation district; they can be natural rivers or artificially excavated channels or ditches. When performing hydrodynamic simulations of an irrigation district, ditches are typically discretized along their direction into a series of computational cross-sections. When measured cross-sectional data is insufficient, the number of computational cross-sections far exceeds the number of measured cross-sections. Therefore, there are a large number of cross-sections without measured data; these cross-sections are the target cross-sections. Target cross-sections can be cross-sections in ditches where cross-sectional data is insufficient; the cross-sectional data can include cross-sectional shape and dimensions, etc. The upstream and downstream cross sections are two reference cross sections located upstream and downstream of the target cross section along the water flow direction, possessing known hydraulic characteristics. For example, the upstream and downstream cross sections can be cross sections with measured topographic data, or they can be the inlet and outlet design cross sections of a ditch, or they can be cross sections with hydraulic characteristics determined by the one-dimensional hydrodynamic modeling method for irrigation areas when the cross section data provided by this invention is insufficient. That is, after the hydraulic characteristics of the target cross section are determined, the target cross section can serve as the upstream or downstream cross section for other cross sections whose hydraulic characteristics need to be determined.

[0030] Hydraulic characteristic interpolation tables are tabular data describing the variation of cross-sectional hydraulic parameters with water level or depth. For example, a hydraulic characteristic interpolation table can be one or more tables indicating the functional relationship between water level (or water depth) and key hydraulic parameters such as cross-sectional area, water surface width, and hydraulic radius. For instance, a hydraulic characteristic interpolation table may include one or more of the following: water level-cross-sectional area interpolation table, water level-water surface width interpolation table, water level-hydraulic radius interpolation table, and cross-sectional area-water level interpolation table. These tables are essential foundational data for numerically solving one-dimensional hydrodynamic models (such as models based on the Saint-Venant equations).

[0031] The process of obtaining the hydraulic characteristic interpolation tables for the upstream and downstream sections can be done by directly calling existing data or by pre-generating them based on the measured topographic data or design parameters of these two sections.

[0032] Step 102: Determine the interpolation coefficients based on the relative positional relationship between the target section and the upstream and downstream sections.

[0033] The relative positional relationship mainly refers to the spatial distribution of the three cross-sections along the centerline of the ditch. For example, the relative positional relationship can be the distance between the target cross-section and the upstream and / or downstream cross-sections. For instance, the mileage coordinates of the upstream cross-section can be obtained. Mileage coordinates of the downstream section and the mileage coordinates of the target section Based on these mileage coordinates, the interpolation coefficients can be calculated. interpolation coefficients The calculation method can be = The interpolation coefficients here This reflects the spatial proximity of the target cross-section to the upstream cross-section. Of course, nonlinear or other more complex algorithms can also be used to determine the interpolation coefficients as needed.

[0034] Step 103: In the hydraulic parameter space, using the hydraulic characteristic interpolation tables of the upstream section and the downstream section as the operation objects, interpolation is performed using interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target section.

[0035] Unlike traditional interpolation of cross-sectional shapes in geometric space, this scheme performs interpolation in hydraulic parameter space, using hydraulic characteristic interpolation tables as the operating objects. Specifically, for corresponding entries in the hydraulic characteristic interpolation tables of upstream and downstream sections, weighted calculations are performed using previously determined interpolation coefficients to obtain the hydraulic characteristics of the target section. For example, the flow cross-sectional area of ​​the target section at a certain water level can be obtained by weighting the flow cross-sectional areas of the upstream and downstream sections at the same water level using interpolation coefficients. By performing such interpolation calculations on a series of water points, a complete hydraulic characteristic interpolation table for the target section can be generated.

[0036] Step 104: Based on the hydraulic characteristic interpolation table of the target section, construct a one-dimensional hydrodynamic model of the irrigation area.

[0037] After obtaining the interpolation tables of hydraulic characteristics for all target cross-sections, the hydraulic data for all cross-sections on the entire canal computational grid were completed. These newly generated interpolation tables were then combined with the existing interpolation tables of known cross-sections in the irrigation area to form a complete and continuous cross-section database. This database can be used as input to any one-dimensional hydrodynamic model based on the Saint-Venant equations or similar theories for subsequent unsteady flow simulations.

[0038] This embodiment directly performs spatial interpolation on the hydraulic characteristic interpolation table in the hydraulic parameter space to generate hydraulic characteristic data for the target section without actual measured data. This bypasses the difficult problem of matching the geometric shape of the section in traditional methods, thus making it applicable to sections of any shape. It realizes automated and batch modeling under the condition of insufficient section data, and provides efficient and reliable data support for building large-scale and refined one-dimensional hydrodynamic models of irrigation areas.

[0039] In one embodiment, the hydraulic characteristic interpolation table for the target cross section includes a water level-flow cross-sectional area interpolation table, a water level-water surface width interpolation table, a water level-hydraulic radius interpolation table, and a flow cross-sectional area-water level interpolation table. Similarly, the hydraulic characteristic interpolation tables for the upstream and downstream cross sections may also include a water level-flow cross-sectional area interpolation table, a water level-water surface width interpolation table, a water level-hydraulic radius interpolation table, and a flow cross-sectional area-water level interpolation table.

[0040] The water level-flow cross-sectional area interpolation table can represent the correspondence between different water levels and flow cross-sectional areas; the water level-water surface width interpolation table can represent the correspondence between different water levels and water surface widths; the water level-hydraulic radius interpolation table can represent the correspondence between different water levels and hydraulic radii; and the flow cross-sectional area-water level interpolation table can represent the correspondence between different flow cross-sectional areas and water level interpolation tables.

[0041] In one embodiment, step 103 above may include, for example: Figure 2 The following steps are shown: Step 201: Convert the water level values ​​at each interpolation point in the hydraulic characteristic interpolation tables of the upstream and downstream sections into relative water depths to obtain reference hydraulic characteristic interpolation tables for the upstream and downstream sections respectively.

[0042] Considering that irrigation ditches typically have longitudinal slopes, the bottom elevations (i.e., lowest point elevations) of upstream and downstream sections often differ. Direct interpolation using absolute water level coordinates can confuse the effects of bottom elevation changes and cross-sectional morphology changes, leading to unclear physical meanings in the interpolation results. To address this issue, this embodiment introduces water depth normalization. Specifically, for each known cross-section (upstream and downstream), its hydraulic characteristic interpolation table is converted from water level based on absolute elevation to relative water depth based on the bottom elevation of each cross-section. The resulting table is the reference hydraulic characteristic interpolation table. Specifically, the four interpolation tables for the upstream and downstream cross-sections can be converted to the corresponding reference hydraulic characteristic interpolation tables.

[0043] In one embodiment, converting the water level values ​​at each interpolation point in the hydraulic characteristic interpolation tables of the upstream and downstream sections into relative water depths may include: performing the following steps on the hydraulic characteristic interpolation tables of the upstream or downstream sections: obtaining the lowest water level value in the hydraulic characteristic interpolation table; and using the difference between the water level value at each interpolation point in the hydraulic characteristic interpolation table and the lowest water level value as the relative water depth at each interpolation point in the hydraulic characteristic interpolation table.

[0044] This lowest water level corresponds to the elevation of the bottom of the ditch at that cross-section, and can be denoted as... When normalizing water depth, the lowest water level value from each hydraulic characteristic interpolation table is used. Subtract the water level value from the interpolation table This yields the water depth value. For example, the lowest water level in the hydraulic characteristic interpolation table for the upstream section is... The water level value at each interpolation point in the hydraulic characteristic interpolation table of the upstream section is reduced by the lowest water level value. The relative water depth at the interpolation point is obtained. The lowest water level in the hydraulic characteristic interpolation table for the downstream section is... Similarly, the water level value at each interpolation point in the hydraulic characteristic interpolation table for the downstream section is reduced by the lowest water level value. The relative water depth at the interpolation point is obtained.

[0045] In one embodiment, the relative water depths in the reference hydraulic characteristic interpolation table of the target section are within the effective water depth range of the target section, and the method further includes: First, determine the first maximum relative water depth in the reference hydraulic characteristic interpolation table for the upstream section and the second maximum relative water depth in the reference hydraulic characteristic interpolation table for the downstream section. These are also the bottom elevations of these two sections. In the converted reference hydraulic characteristic interpolation table, the maximum water depth is taken as the effective water depth range for that section; that is, the first maximum relative water depth of the upstream section is denoted as... The second maximum relative water depth at the downstream section is denoted as .

[0046] Then, based on the interpolation coefficients, the first maximum relative water depth, and the second maximum relative water depth, the effective water depth range of the target section is obtained.

[0047] The effective water depth range of the target cross-section can be from 0 to the maximum relative water depth. The maximum relative water depth of the target cross-section... The determination method can be: ,in These are the interpolation coefficients.

[0048] This embodiment ensures that the interpolated target section data covers a reasonable and necessary water depth range by limiting the effective water depth range of the target section, reducing meaningless extrapolation and guaranteeing the validity of the data.

[0049] Step 202: Interpolate using interpolation coefficients and reference hydraulic characteristic interpolation tables for upstream and downstream sections to obtain the reference hydraulic characteristic interpolation table for the target section.

[0050] Under a unified relative water depth coordinate system, the differences in various hydraulic parameters (such as cross-sectional area and hydraulic radius) between upstream and downstream sections at the same relative water depth purely reflect the changes in cross-sectional morphology along the course of the flow. In this case, the interpolation coefficients determined in Example 1 can be used to interpolate the reference hydraulic characteristic interpolation tables for the upstream and downstream sections.

[0051] In one embodiment, the reference hydraulic characteristic interpolation table for the target cross section includes a relative water depth-flow cross section area interpolation table, a relative water depth-water surface width interpolation table, a relative water depth-hydraulic radius interpolation table, and a flow cross section area-relative water depth interpolation table.

[0052] The relative water depth-flow cross-sectional area interpolation table represents the correspondence between different relative water depths and flow cross-sectional areas; the relative water depth-water surface width interpolation table represents the correspondence between different relative water depths and water surface widths; and the relative water depth-hydraulic radius interpolation table represents the correspondence between different relative water depths and hydraulic radii. The flow cross-sectional area-relative water depth interpolation table represents the correspondence between different flow cross-sectional areas and relative water depths, and it can be derived from the relative water depth-flow cross-sectional area interpolation table.

[0053] Considering the continuity equations in Saint-Venant's equations, the flow area... With water surface width There are strict differential constraints between them: That is, the width of the water surface equals the derivative of the flow area with respect to the relative water depth. This relationship is directly derived from the geometric definition of the cross-section and is a strictly valid physical constraint. If for... and Perform independent linear interpolation separately, and the interpolated results are as follows. and May no longer be satisfied This constraint introduces errors due to mass non-conservation in subsequent hydrodynamic model calculations. To address this issue, this embodiment employs a strategy that uses the cross-sectional area of ​​the flow path as the primary interpolation parameter, while the water surface width is derived from the primary interpolation parameter. This eliminates the inherent inconsistency introduced by the independent interpolation of the two parameters without increasing the input data requirements, ensuring the physical self-consistency of the hydraulic characteristic parameters generated by the interpolation.

[0054] That is, step 202 above may include the following steps: determining the relative water depth-flow cross-sectional area interpolation table of the target cross-section using interpolation coefficients and the relative water depth-flow cross-sectional area interpolation table of the upstream and downstream cross-sections; and determining the relative water depth-hydraulic radius interpolation table of the target cross-section using interpolation coefficients and the relative water depth-hydraulic radius interpolation table of the upstream and downstream cross-sections. Then, numerical differentiation is performed on each flow cross-sectional area in the relative water depth-flow cross-sectional area interpolation table of the target cross-section to obtain the relative water depth-water surface width interpolation table of the target cross-section. Finally, the relative water depth-flow cross-sectional area interpolation table is transformed to obtain the flow cross-sectional area-relative water depth interpolation table.

[0055] For example, for each relative water depth obtained by sampling within the effective water depth of the target cross section The cross-sectional area of ​​the target section The cross-sectional area of ​​the upstream section can be used to determine the flow rate. and the cross-sectional area of ​​the downstream section We obtain this through linear weighting: ,in These are the interpolation coefficients. By performing this operation on a series of sampled relative water depths, a relative water depth-flow cross-sectional area interpolation table for the target cross-section can be obtained.

[0056] Similarly, for each relative water depth obtained by sampling within the effective water depth of the target cross section Hydraulic radius of the target section It can be determined by the hydraulic radius of the upstream section. Hydraulic radius of downstream section We obtain this through linear weighting: ,in These are the interpolation coefficients. By performing this operation on a series of sampled relative water depths, a relative water depth-hydraulic radius interpolation table for the target cross-section can be obtained.

[0057] For example, considering that the relative water depth-flow cross-sectional area interpolation table is discrete data, a central difference scheme can be used for calculation at the center interpolation point. For instance, for the relative water depth point... Its water surface is wide can be = The calculation yielded the result. For the end interpolation points, a one-sided difference method can be used.

[0058] This embodiment uses a strategy that prioritizes the cross-sectional area of ​​the flow path as the primary interpolation parameter, while the water surface width is derived from the primary interpolation parameter. This strategy ensures the internal physical consistency of the interpolated hydraulic characteristic parameters. This constraint-consistent interpolation strategy guarantees that the generated target cross-sectional data satisfies the inherent constraints of the Saint-Venant equations, reduces errors introduced by independent interpolation, and improves the stability and accuracy of the final one-dimensional hydrodynamic model calculation.

[0059] In addition, the water surface width and cross-sectional area at each relative water depth determined by numerical differentiation strictly satisfy the differential constraint relationship, making the constructed one-dimensional hydrodynamic model more accurate.

[0060] Step 203: Convert the relative water depths in the reference hydraulic characteristic interpolation table of the target section to obtain the hydraulic characteristic interpolation table of the target section.

[0061] This step can be considered as the inverse transformation of the water depth normalized coordinates. After obtaining the reference hydraulic characteristic interpolation table of the target section, it needs to be transformed from relative water depth coordinates back to absolute water level coordinates for use by the one-dimensional hydrodynamic model.

[0062] This embodiment effectively decouples the longitudinal slope effect of the ditch from the cross-sectional shape variation effect by introducing water depth normalization, making the interpolation process more physically reasonable in hydraulics and significantly improving the accuracy of the hydraulic characteristic data of the target cross-section generated by interpolation.

[0063] In one embodiment, step 203 above may include the following steps: First, obtain the bottom elevation of the target section. The bottom elevation is based on interpolation coefficients. , bottom elevation of the upstream section and the bottom elevation of the downstream section Obtained. For example, the bottom elevation of the target section. The determination method can be: .

[0064] Then, based on the bottom elevation, the relative water depth of each interpolation point in the reference hydraulic characteristic interpolation table of the target section is converted to obtain the basic hydraulic characteristic interpolation table of the target section.

[0065] For example, by adding the bottom elevation value to the relative water depth value of each row of the reference hydraulic characteristic interpolation table, the final hydraulic characteristic interpolation table of the target section can be obtained.

[0066] Then, the basic hydraulic characteristic interpolation table is sparsified to obtain the hydraulic characteristic interpolation table of the target section.

[0067] Among them, the basic hydraulic characteristic interpolation table can be sparsed according to the linearity of the basic hydraulic characteristic interpolation table to obtain the hydraulic characteristic interpolation table of the target section.

[0068] Alternatively, in one embodiment, sparse processing of the basic hydraulic characteristic interpolation table to obtain the hydraulic characteristic interpolation table of the target section may include the following steps: first, dividing the basic hydraulic characteristic interpolation table of the target section into segments to obtain several interpolation segments; then, for each interpolation segment, determining the interpolation points to be removed within the interpolation segment based on the linearity of the interpolation segment, and removing the interpolation points to be removed to obtain the hydraulic characteristic interpolation table of the target section.

[0069] One method for dividing the basic hydraulic characteristic interpolation table of the target section into segments to obtain several interpolation segments is as follows: The interpolation points in the basic hydraulic interpolation table are sequentially numbered; a numbering range is selected based on the sequential numbering of each interpolation point; and interpolation segments are constructed based on the interpolation points corresponding to the numbering range. For example, interpolation segments can be constructed by selecting interpolation points with a numbering range from 1 to n.

[0070] The above method of determining whether an interpolation segment contains interpolation points to be removed based on the linearity of the interpolation segment, and removing these interpolation points to obtain the hydraulic characteristic interpolation table of the target section can be as follows: Based on the two endpoints of the interpolation segment, construct a linear equation for the interpolation segment, detect the deviation between all interpolation points in the interpolation segment and the linear equation, and if there are at least three consecutive interpolation points whose deviations are less than a threshold, then it is determined that there are interpolation points to be extracted among these at least three consecutive interpolation points. These interpolation points to be extracted can be the middle interpolation point among the at least three interpolation points, rather than the end interpolation point.

[0071] Optionally, if the error between any interpolation point in the interpolation segment and the linear equation is not less than a threshold, the interpolation point with the largest deviation from the linear equation in the interpolation segment can be detected as the target interpolation point. Then, according to the sequential numbering corresponding to the target interpolation point, a new numbering range is selected, and the interpolation segment is reconstructed. That is, an interpolation segment can be constructed first, and after the processing of the interpolation segment is completed, the next interpolation segment can be constructed, and so on.

[0072] In one embodiment, before converting the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section to obtain the hydraulic characteristic interpolation table of the target cross-section, the method further includes: Perform physical rationality checks on the reference hydraulic characteristic interpolation table. Physical rationality checks include monotonicity checks, non-negativity checks, and / or endpoint continuity checks.

[0073] If the physical rationality verification result is passed, the step of converting the relative water depths in the reference hydraulic characteristic interpolation table of the target section to obtain the hydraulic characteristic interpolation table of the target section is executed.

[0074] That is, the physical rationality of the interpolation tables of relative water depth-flow cross-sectional area, relative water depth-water surface width, relative water depth-hydraulic radius, and flow cross-sectional area-relative water depth for the target cross-section is verified respectively.

[0075] If the relative water depth-flow cross-sectional area interpolation table, the relative water depth-water surface width interpolation table, the relative water depth-hydraulic radius interpolation table, and the flow cross-sectional area-relative water depth interpolation table of the target cross-section all pass the physical rationality verification, then the step of converting each relative water depth in the reference hydraulic characteristic interpolation table of the target cross-section is performed to obtain the hydraulic characteristic interpolation table of the target cross-section.

[0076] If the rationality check fails, an error message is generated and the current interpolation process for the target section is terminated. In some application scenarios, terminating the current interpolation process means pausing the step of converting the relative water depths in the reference hydraulic characteristic interpolation table of the target section to obtain the hydraulic characteristic interpolation table of the target section. After terminating the current interpolation process, the hydraulic characteristic interpolation table can be automatically corrected first. Once the automatically corrected reference hydraulic characteristic interpolation table passes the physical rationality check, the current interpolation process can be restarted, i.e., the step of converting the relative water depths in the reference hydraulic characteristic interpolation table of the target section to obtain the hydraulic characteristic interpolation table of the target section will be initiated.

[0077] In this embodiment, after verifying the physical rationality of the reference hydraulic characteristic interpolation table of the target section, the basic hydraulic characteristic interpolation table of the target section is obtained, which makes the determined basic hydraulic characteristic interpolation table more accurate.

[0078] In one embodiment, obtaining the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section corresponding to the target section to be interpolated in the irrigation area's ditches includes: Topographic data and hydrodynamic interpolation parameters for the upstream and downstream cross-sections are acquired separately. The topographic data includes a series of cross-section scatter coordinates (distance from starting point, elevation), arranged in ascending order of distance from starting point. The hydrodynamic interpolation parameters include water level information and accuracy information, specifically the minimum water level (default is the minimum elevation in the cross-section scatter coordinates), the maximum water level (default is the maximum elevation in the cross-section scatter coordinates), the minimum interpolation interval (default is 0.01m), and the relative interpolation accuracy (default is 1×10⁻⁶). -5 ).

[0079] Then, based on the topographic data and hydrodynamic interpolation parameters of the upstream section, a hydraulic characteristic interpolation table for the upstream section is constructed; and based on the topographic data and hydrodynamic interpolation parameters of the downstream section, a hydraulic characteristic interpolation table for the downstream section is constructed.

[0080] For example, the method for constructing the hydraulic characteristic interpolation table of the upstream section based on the topographic data and hydrodynamic interpolation parameters of the upstream section can be as follows: Construct a basic hydraulic characteristic interpolation table (e.g., a basic water level-water surface width interpolation table, a water level-flow cross-sectional area interpolation table, a water level-hydraulic radius interpolation table, and a flow cross-sectional area-water level interpolation table) based on the topographic data and hydrodynamic interpolation parameters. For example, the correspondence between different water levels and water surface width, flow cross-sectional area, and hydraulic radius can be determined based on the obtained topographic data and hydrodynamic interpolation parameters, and then the basic water level-water surface width interpolation table, water level-flow cross-sectional area interpolation table, and water level-hydraulic radius interpolation table can be constructed according to this correspondence. Furthermore, since the flow cross-sectional area and water level have a monotonically increasing relationship, the basic flow cross-sectional area-water level interpolation table can be obtained by transforming the basic water level-flow cross-sectional area interpolation table.

[0081] Then, the basic hydraulic characteristic interpolation table is sparsified to obtain the hydraulic characteristic interpolation table. The method for sparsifying the basic hydraulic characteristic interpolation table to obtain the hydraulic characteristic interpolation table can be as follows: the basic hydraulic characteristic interpolation table is divided into segments to obtain several interpolation segments; for each interpolation segment, the number of interpolation points to be removed is determined based on the linearity of the interpolation segment, and these points are then removed to obtain the hydraulic characteristic interpolation table. The specific process is similar to the method described above for sparsifying the basic hydraulic characteristic interpolation table of the target section, and will not be elaborated here. The method for constructing the hydraulic characteristic interpolation table of the downstream section based on the topographic data and hydrodynamic interpolation parameters is similar and will not be elaborated here.

[0082] In this embodiment, by sparse processing of the basic hydraulic characteristic interpolation table, it is possible to reduce memory overhead and improve computational efficiency while preserving the hydraulic characteristics at different water levels.

[0083] Considering that irrigation systems consist of canals (water distribution), ditches (drainage), and rivers (natural, serving both water diversion and drainage functions), the data conditions and cross-sectional characteristics of each level of the water system exhibit systematic differences, making them unsuitable for single-method processing. This invention addresses this reality by proposing a hierarchical adaptive processing strategy. Within a unified hydraulic characteristic interpolation table framework, it provides appropriate processing paths for canal sections with different data conditions. In one embodiment, the method further includes: First, determine the data conditions for each ditch in the irrigation area. Data conditions include scattered data with multiple cross-sections, scattered data with both inlet and outlet cross-sections, and data containing the ditch's design parameters or partial shape parameters. Scattered data with multiple cross-sections can be considered for important main canals or natural rivers with a large number of measured cross-sections. Scattered data with both inlet and outlet cross-sections can be considered for some branch canals, with measurement data only at their beginning and end. Data containing the ditch's design parameters can be considered for standard lined canals, with known bottom width, slope coefficient, design water depth, etc. Data containing partial shape parameters can be considered for some small ditches, with only the approximate opening width and depth known.

[0084] Then, for each ditch, if the ditch's data conditions are scattered data with multiple cross-sections or scattered data with two cross-sections (inlet and outlet), the cross-section to be interpolated is taken as the target cross-section, and the two cross-sections with scattered data are taken as the upstream and downstream cross-sections of the target cross-section, respectively.

[0085] In other words, for each ditch, a corresponding processing path is adopted based on its data conditions. Specifically, if the data conditions of the ditch are "scattered data with multiple cross-sections" or "scattered data with two cross-sections, namely inlet and outlet," then the method in the aforementioned embodiments of the present invention can be applied. In this case, the cross-section that needs to have its data supplemented can be taken as the target cross-section, and two spatially adjacent cross-sections with scattered data can be taken as the upstream and downstream cross-sections of the target cross-section, respectively, for interpolation calculation. For ditches with other data conditions, hydraulic characteristic interpolation tables can be directly generated using analytical formulas or other methods.

[0086] All four processing paths output hydraulic characteristic interpolation tables in the same format, completing the cross-sectional data of the entire irrigation district's water system within a unified table framework. Finally, these data are combined into a complete cross-sectional dataset that can be directly used in a one-dimensional hydrodynamic model. This strategy covers the most common data conditions encountered in actual irrigation district modeling, supports automated batch modeling of the entire irrigation district's water system, helps reduce manual modeling workload, and adapts to the engineering application needs of digital twin irrigation districts.

[0087] This embodiment proposes a systematic modeling strategy for complex water systems in irrigation districts. By classifying ditches under different data conditions, the interpolation method proposed in this invention can work in conjunction with other methods, greatly enhancing the universality and automation capabilities of this invention in practical engineering applications, and adapting to the needs of large-scale modeling such as digital twin irrigation districts.

[0088] In one embodiment, a method for modeling a one-dimensional hydrodynamic model of an irrigation district when cross-sectional data is insufficient may include, for example: Figure 3 The following steps are shown: Step 301: Calculate the interpolation coefficients based on the target section, the upstream section of the target section, and the mileage coordinates of the downstream section.

[0089] For example, the mileage coordinates of the upstream section. Mileage coordinates of the downstream section and the mileage coordinates of the target section Based on these mileage coordinates, the interpolation coefficients can be calculated. For example, interpolation coefficients The calculation method can be = .

[0090] Step 302: Convert the hydraulic characteristic interpolation tables of the upstream and downstream sections from absolute water level to relative water depth, and determine the bottom elevation and effective water depth range of the target section.

[0091] For example, when normalizing water depth, the lowest water level value in each hydraulic characteristic interpolation table is taken, and this lowest water level value is subtracted from each water level value in the interpolation table to obtain the relative water depth at each interpolation point in the table. For instance, the lowest water level value in the hydraulic characteristic interpolation table for the upstream section is... The lowest water level in the hydraulic characteristic interpolation table for the downstream section is [value missing]. Among them, the lowest water level value It is also the bottom elevation of the upstream section; similarly, the lowest water level value. This is the bottom elevation of the downstream section. In the converted reference hydraulic characteristic interpolation table, the maximum water depth is taken as the effective water depth range for this section; that is, the maximum relative water depth of the upstream section is denoted as... The maximum relative water depth at the downstream section is denoted as .

[0092] Using the interpolation coefficients determined in step 301, the bottom elevation and maximum relative water depth of the interpolation section are determined using linear interpolation. For example, the maximum relative water depth of the target section... The determination method can be: ,in These are the interpolation coefficients. The bottom elevation of the target section. The determination method can be: .

[0093] Step 303: Generate a relative water depth sampling sequence within the effective water depth range of the target section and determine the reference hydraulic characteristic interpolation table of the target section.

[0094] The relative water depth sampling sequence includes several relative water depths obtained from sampling. For example, from 0 to... Sampling points are generated at 1cm intervals. For any relative water depth obtained from the sampling, the corresponding cross-sectional area is calculated. The calculation method can be: Similarly, its hydraulic radius The calculation method can be: This yields the relative water depth-flow cross-sectional area interpolation table and the relative water depth-hydraulic radius interpolation table for the target cross-section. The relative water depth-water surface width interpolation table for the target cross-section is obtained by numerical differentiation based on the relative water depth-flow cross-sectional area interpolation table. For example, for a relative water depth point... Its water surface is wide can be = The calculations yielded the target cross-sectional area-relative water depth interpolation table. Additionally, the cross-sectional area-relative water depth interpolation table for the target cross-section was obtained by reverse mapping the relative water depth-flow cross-sectional area interpolation table.

[0095] Step 304: Perform a physical rationality check on the reference hydraulic characteristic interpolation table of the target section.

[0096] The reference hydraulic characteristic interpolation tables include relative water depth-water surface width interpolation tables, relative water depth-flow cross-sectional area interpolation tables, relative water depth-hydraulic radius interpolation tables, and flow cross-sectional area-relative water depth interpolation tables. The physical rationality of the relative water depth-water surface width interpolation tables, relative water depth-flow cross-sectional area interpolation tables, relative water depth-hydraulic radius interpolation tables, and flow cross-sectional area-relative water depth interpolation tables is verified.

[0097] Step 305: If the physical rationality verification is passed, convert the relative water depth of each interpolation point in the reference hydraulic characteristic interpolation table of the target section according to the bottom elevation to obtain the basic hydraulic characteristic interpolation table of the target section.

[0098] The foundation hydraulic characteristic interpolation tables include the foundation water level-water surface width interpolation table, the water level-flow cross-sectional area interpolation table, the water level-hydraulic radius interpolation table, and the flow cross-sectional area-water level interpolation table.

[0099] Optionally, if the physical rationality check fails, an error message is generated and the current interpolation process for the target section is terminated.

[0100] Step 306: Sparsify the basic hydraulic characteristic interpolation table to obtain the hydraulic characteristic interpolation table of the target section.

[0101] The basic water level-water surface width interpolation table, water level-flow cross-sectional area interpolation table, water level-hydraulic radius interpolation table, and flow cross-sectional area-water level interpolation table are sparsed to obtain the target cross-section's water level-water surface width interpolation table, water level-flow cross-sectional area interpolation table, water level-hydraulic radius interpolation table, and flow cross-sectional area-water level interpolation table.

[0102] Step 307: Based on the hydraulic characteristic interpolation table of the target section, construct a one-dimensional hydrodynamic model of the irrigation area.

[0103] That is, a one-dimensional hydrodynamic model of the irrigation area is constructed based on the water level-water surface width interpolation table, the water level-flow cross-sectional area interpolation table, the water level-hydraulic radius interpolation table, and the flow cross-sectional area-water level interpolation table of the target cross-section.

[0104] The following describes the modeling device for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, as provided by the present invention. The modeling device for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, as described below, can be referred to in correspondence with the modeling method for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, as described above. Figure 4 As shown, the one-dimensional hydrodynamic modeling device 400 for irrigation areas when cross-sectional data is insufficient may include the following modules: The hydraulic characteristic interpolation table acquisition module 401 is used to acquire the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section corresponding to the target section to be interpolated in the irrigation area's ditches. The interpolation coefficient determination module 402 is used to determine the interpolation coefficients based on the relative positional relationship between the target section and the upstream section and the downstream section. Interpolation module 403 is used to perform interpolation in the hydraulic parameter space, using the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section as operation objects, and using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target section. The one-dimensional hydrodynamic model construction module 404 is used to construct a one-dimensional hydrodynamic model of the irrigation area based on the hydraulic characteristic interpolation table of the target section.

[0105] According to the present invention, a one-dimensional hydrodynamic modeling device 400 for irrigation areas with insufficient cross-sectional data includes an interpolation module 403 that, in the hydraulic parameter space, uses the hydraulic characteristic interpolation tables of the upstream and downstream cross-sections as operation objects, and performs interpolation using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target cross-section. The interpolation module 403 includes: The water level values ​​at each interpolation point in the hydraulic characteristic interpolation tables of the upstream and downstream sections are converted into relative water depths to obtain reference hydraulic characteristic interpolation tables for the upstream and downstream sections respectively. Interpolation is performed using the interpolation coefficients and the reference hydraulic characteristic interpolation tables for the upstream and downstream sections to obtain the reference hydraulic characteristic interpolation table for the target section. The relative water depths in the reference hydraulic characteristic interpolation table of the target cross section are converted to obtain the hydraulic characteristic interpolation table of the target cross section.

[0106] According to the present invention, a modeling device 400 for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient is provided. The reference hydraulic characteristic interpolation table of the target cross-section includes a relative water depth-flow cross-sectional area interpolation table, a relative water depth-water surface width interpolation table, a relative water depth-hydraulic radius interpolation table, and a flow cross-sectional area-relative water depth interpolation table. The interpolation module 403 uses the interpolation coefficients and the reference hydraulic characteristic interpolation table of the upstream section and the downstream section to perform interpolation, thereby obtaining the reference hydraulic characteristic interpolation table of the target section, including: Using the interpolation coefficients and the relative water depth-flow cross-sectional area interpolation table between the upstream and downstream sections, the relative water depth-flow cross-sectional area interpolation table of the target section is determined; and using the interpolation coefficients and the relative water depth-hydraulic radius interpolation table between the upstream and downstream sections, the relative water depth-hydraulic radius interpolation table of the target section is determined. Numerical differentiation is performed on each flow cross-sectional area in the relative water depth-flow cross-sectional area interpolation table of the target cross-section to obtain the relative water depth-water surface width interpolation table of the target cross-section. The relative water depth-flow cross-sectional area interpolation table is converted to obtain the flow cross-sectional area-relative water depth interpolation table.

[0107] According to the present invention, a one-dimensional hydrodynamic modeling device 400 for irrigation areas with insufficient cross-sectional data includes an interpolation module 403 that converts the water level values ​​at each interpolation point in the hydraulic characteristic interpolation table of the upstream and downstream cross-sections into relative water depths, including: For the hydraulic characteristic interpolation table of the upstream section or the downstream section, perform the following steps: Obtain the lowest water level value from the hydraulic characteristic interpolation table; The difference between the water level value at each interpolation point in the hydraulic characteristic interpolation table and the lowest water level value is taken as the relative water depth at each interpolation point in the hydraulic characteristic interpolation table.

[0108] According to the present invention, a one-dimensional hydrodynamic modeling device 400 for irrigation areas with insufficient cross-sectional data is provided, wherein the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section are within the effective water depth range of the target cross-section, and the device 400 further includes: Determine the first maximum relative water depth in the reference hydraulic characteristic interpolation table of the upstream section and the second maximum relative water depth in the reference hydraulic characteristic interpolation table of the downstream section; The effective water depth range of the target cross section is obtained based on the interpolation coefficient, the first maximum relative water depth, and the second maximum relative water depth.

[0109] According to the present invention, a one-dimensional hydrodynamic modeling device 400 for irrigation areas with insufficient cross-sectional data includes an interpolation module 403 that converts the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section to obtain the hydraulic characteristic interpolation table of the target cross-section, comprising: Obtain the bottom elevation of the target section, which is based on the interpolation coefficient, the bottom elevation of the upstream section, and the bottom elevation of the downstream section; Based on the bottom elevation, the relative water depth of each interpolation point in the reference hydraulic characteristic interpolation table of the target section is converted to obtain the basic hydraulic characteristic interpolation table of the target section; The basic hydraulic characteristic interpolation table is sparsified to obtain the hydraulic characteristic interpolation table of the target section.

[0110] According to the present invention, a modeling device 400 for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient includes an interpolation module 403 that performs sparsification processing on the basic hydraulic characteristic interpolation table to obtain the hydraulic characteristic interpolation table of the target cross-section, comprising: The basic hydraulic characteristic interpolation table of the target section is divided into segments to obtain several interpolation segments; For each interpolation segment, the interpolation points to be removed are determined based on the linearity of the interpolation segment, and the interpolation points to be removed are processed to obtain the hydraulic characteristic interpolation table of the target section.

[0111] According to the present invention, a one-dimensional hydrodynamic modeling device 400 for irrigation areas with insufficient cross-sectional data, before converting the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section to obtain the hydraulic characteristic interpolation table of the target cross-section, the interpolation module 403 is further configured to: The physical rationality of the reference hydraulic characteristic interpolation table is verified, including monotonicity verification, non-negativity verification and / or endpoint continuity verification. If the result of the physical rationality verification is that the verification is passed, the step of converting the relative water depths in the reference hydraulic characteristic interpolation table of the target section to obtain the hydraulic characteristic interpolation table of the target section is executed. If the physical rationality check fails, an error message is generated and the current interpolation process for the target section is terminated.

[0112] According to the present invention, a one-dimensional hydrodynamic modeling device 400 for irrigation areas when cross-sectional data is insufficient includes a hydraulic characteristic interpolation table acquisition module 401 that acquires the hydraulic characteristic interpolation tables of the upstream and downstream sections corresponding to the target cross-section to be interpolated in the irrigation area's ditches, comprising: The topographic data and hydrodynamic interpolation parameters of the upstream section and the downstream section are obtained respectively. Based on the topographic data of the upstream section and the hydrodynamic interpolation parameters, a hydraulic characteristic interpolation table for the upstream section is constructed; and based on the topographic data of the downstream section and the hydrodynamic interpolation parameters, a hydraulic characteristic interpolation table for the downstream section is constructed.

[0113] According to the present invention, a modeling device 400 for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, wherein the hydraulic characteristic interpolation table acquisition module 401 is further used for: The data conditions for each ditch in the irrigation area are determined. The data conditions include scattered data with multiple cross-sections, scattered data with two cross-sections (inlet and outlet), design parameters of the ditch, or partial shape parameters of the ditch. For each ditch, if the data condition of the ditch is scattered data with multiple cross-sections or scattered data with two cross-sections (inlet and outlet), the cross-section to be interpolated is taken as the target cross-section, and the two cross-sections with scattered data are taken as the upstream and downstream cross-sections of the target cross-section, respectively.

[0114] Figure 5 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 5 As shown, the electronic device may include: a processor 510, a communication interface 520, a memory 530, and a communication bus 540, wherein the processor 510, the communication interface 520, and the memory 530 communicate with each other through the communication bus 540. The processor 510 can call logical instructions in the memory 530 to execute a one-dimensional hydrodynamic modeling method for irrigation districts when cross-sectional data is insufficient. This method includes: obtaining the hydraulic characteristic interpolation tables of the upstream and downstream cross-sections corresponding to the target cross-section to be interpolated in the irrigation district's ditches; determining interpolation coefficients based on the relative positional relationship between the target cross-section and the upstream and downstream cross-sections; in the hydraulic parameter space, using the hydraulic characteristic interpolation tables of the upstream and downstream cross-sections as operation objects, performing interpolation using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target cross-section; and constructing a one-dimensional hydrodynamic model of the irrigation district based on the hydraulic characteristic interpolation table of the target cross-section.

[0115] Furthermore, the logical instructions in the aforementioned memory 530 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0116] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute a method for modeling a one-dimensional hydrodynamic model of an irrigation area when the cross-sectional data provided by the above methods is insufficient. This method includes: obtaining hydraulic characteristic interpolation tables of the upstream and downstream sections corresponding to the target section to be interpolated in the irrigation area's ditches; determining interpolation coefficients based on the relative positional relationship between the target section and the upstream and downstream sections; in the hydraulic parameter space, using the hydraulic characteristic interpolation tables of the upstream and downstream sections as operation objects, performing interpolation using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target section; and constructing a one-dimensional hydrodynamic model of the irrigation area based on the hydraulic characteristic interpolation table of the target section.

[0117] In another aspect, the present invention also provides a non-transitory computer-readable storage medium storing a computer program thereon. When executed by a processor, this computer program implements a method for modeling a one-dimensional hydrodynamic model of an irrigation district when the cross-sectional data provided by the methods described above is insufficient. The method includes: obtaining hydraulic characteristic interpolation tables for upstream and downstream sections of the target section to be interpolated in the irrigation district's ditches; determining interpolation coefficients based on the relative positional relationship between the target section and the upstream and downstream sections; performing interpolation using the interpolation coefficients in the hydraulic parameter space, with the hydraulic characteristic interpolation tables of the upstream and downstream sections as the operation objects, to obtain the hydraulic characteristic interpolation table of the target section; and constructing a one-dimensional hydrodynamic model of the irrigation district based on the hydraulic characteristic interpolation table of the target section.

[0118] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0119] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0120] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for modeling a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, characterized in that, include: Obtain the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section corresponding to the target section to be interpolated in the irrigation area's ditches; Based on the relative positional relationship between the target section and the upstream and downstream sections, the interpolation coefficients are determined; In the hydraulic parameter space, the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section are used as the operation objects, and the interpolation coefficients are used to perform interpolation to obtain the hydraulic characteristic interpolation table of the target section. Based on the hydraulic characteristic interpolation table of the target section, a one-dimensional hydrodynamic model of the irrigation area is constructed.

2. The method according to claim 1, characterized in that, In the hydraulic parameter space, using the hydraulic characteristic interpolation tables of the upstream section and the downstream section as the operation objects, interpolation is performed using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target section, including: The water level values ​​at each interpolation point in the hydraulic characteristic interpolation tables of the upstream and downstream sections are converted into relative water depths to obtain reference hydraulic characteristic interpolation tables for the upstream and downstream sections respectively. Interpolation is performed using the interpolation coefficients and the reference hydraulic characteristic interpolation tables for the upstream and downstream sections to obtain the reference hydraulic characteristic interpolation table for the target section. The relative water depths in the reference hydraulic characteristic interpolation table of the target cross section are converted to obtain the hydraulic characteristic interpolation table of the target cross section.

3. The method according to claim 2, characterized in that, The reference hydraulic characteristic interpolation tables for the target cross section include a relative water depth-flow cross section area interpolation table, a relative water depth-water surface width interpolation table, a relative water depth-hydraulic radius interpolation table, and a flow cross section area-relative water depth interpolation table. Interpolation is performed using the interpolation coefficients and the reference hydraulic characteristic interpolation tables for the upstream and downstream sections to obtain the reference hydraulic characteristic interpolation table for the target section, including: Using the interpolation coefficients and the relative water depth-flow cross-sectional area interpolation table between the upstream and downstream sections, the relative water depth-flow cross-sectional area interpolation table of the target section is determined; and using the interpolation coefficients and the relative water depth-hydraulic radius interpolation table between the upstream and downstream sections, the relative water depth-hydraulic radius interpolation table of the target section is determined. Numerical differentiation is performed on each flow cross-sectional area in the relative water depth-flow cross-sectional area interpolation table of the target cross-section to obtain the relative water depth-water surface width interpolation table of the target cross-section. The relative water depth-flow cross-sectional area interpolation table is converted to obtain the flow cross-sectional area-relative water depth interpolation table.

4. The method according to claim 2, characterized in that, Converting the water level values ​​at each interpolation point in the hydraulic characteristic interpolation table for the upstream and downstream sections into relative water depths includes: For the hydraulic characteristic interpolation table of the upstream section or the downstream section, perform the following steps: Obtain the lowest water level value from the hydraulic characteristic interpolation table; The difference between the water level value at each interpolation point in the hydraulic characteristic interpolation table and the lowest water level value is taken as the relative water depth at each interpolation point in the hydraulic characteristic interpolation table.

5. The method according to claim 4, characterized in that, The relative water depths in the reference hydraulic characteristic interpolation table of the target section are within the effective water depth range of the target section, and the method further includes: Determine the first maximum relative water depth in the reference hydraulic characteristic interpolation table of the upstream section and the second maximum relative water depth in the reference hydraulic characteristic interpolation table of the downstream section; The effective water depth range of the target cross section is obtained based on the interpolation coefficient, the first maximum relative water depth, and the second maximum relative water depth.

6. The method according to any one of claims 2 to 5, characterized in that, The relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section are converted to obtain the hydraulic characteristic interpolation table of the target cross-section, including: Obtain the bottom elevation of the target section, which is based on the interpolation coefficient, the bottom elevation of the upstream section, and the bottom elevation of the downstream section; Based on the bottom elevation, the relative water depth of each interpolation point in the reference hydraulic characteristic interpolation table of the target section is converted to obtain the basic hydraulic characteristic interpolation table of the target section; The basic hydraulic characteristic interpolation table is sparsified to obtain the hydraulic characteristic interpolation table of the target section.

7. The method according to claim 6, characterized in that, The basic hydraulic characteristic interpolation table is sparsified to obtain the hydraulic characteristic interpolation table of the target section, including: The basic hydraulic characteristic interpolation table of the target section is divided into segments to obtain several interpolation segments; For each interpolation segment, the interpolation points to be removed are determined based on the linearity of the interpolation segment, and the interpolation points to be removed are processed to obtain the hydraulic characteristic interpolation table of the target section.

8. The method according to any one of claims 2 to 5, characterized in that, Before converting the relative water depths in the reference hydraulic characteristic interpolation table of the target cross-section to obtain the hydraulic characteristic interpolation table of the target cross-section, the method further includes: The physical rationality of the reference hydraulic characteristic interpolation table is verified, including monotonicity verification, non-negativity verification and / or endpoint continuity verification. If the result of the physical rationality verification is that the verification is passed, the step of converting the relative water depths in the reference hydraulic characteristic interpolation table of the target section to obtain the hydraulic characteristic interpolation table of the target section is executed. If the physical rationality check fails, an error message is generated and the current interpolation process for the target section is terminated.

9. The method according to any one of claims 1 to 5, characterized in that, Obtain the hydraulic characteristic interpolation tables for the upstream and downstream sections of the target cross-section in the irrigation area's ditches, including: The topographic data and hydrodynamic interpolation parameters of the upstream section and the downstream section are obtained respectively. Based on the topographic data of the upstream section and the hydrodynamic interpolation parameters, a hydraulic characteristic interpolation table for the upstream section is constructed; and based on the topographic data of the downstream section and the hydrodynamic interpolation parameters, a hydraulic characteristic interpolation table for the downstream section is constructed.

10. The method according to any one of claims 1 to 5, characterized in that, The method further includes: The data conditions for each ditch in the irrigation area are determined. The data conditions include scattered data with multiple cross-sections, scattered data with two cross-sections (inlet and outlet), design parameters of the ditch, or partial shape parameters of the ditch. For each ditch, if the data condition of the ditch is scattered data with multiple cross-sections or scattered data with two cross-sections (inlet and outlet), the cross-section to be interpolated is taken as the target cross-section, and the two cross-sections with scattered data are taken as the upstream and downstream cross-sections of the target cross-section, respectively.

11. A modeling device for a one-dimensional hydrodynamic model of an irrigation area when cross-sectional data is insufficient, characterized in that, include: The hydraulic characteristic interpolation table acquisition module is used to acquire the hydraulic characteristic interpolation table of the upstream section and the hydraulic characteristic interpolation table of the downstream section corresponding to the target section to be interpolated in the irrigation area's ditches. The interpolation coefficient determination module is used to determine the interpolation coefficients based on the relative positional relationship between the target section and the upstream section and the downstream section. The interpolation module is used to perform interpolation in the hydraulic parameter space, using the hydraulic characteristic interpolation tables of the upstream section and the downstream section as operation objects, and using the interpolation coefficients to obtain the hydraulic characteristic interpolation table of the target section. The one-dimensional hydrodynamic model construction module is used to construct a one-dimensional hydrodynamic model of the irrigation area based on the hydraulic characteristic interpolation table of the target cross section.

12. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that, When the processor executes the computer program, it implements the one-dimensional hydrodynamic modeling method for irrigation areas when cross-sectional data is insufficient, as described in any one of claims 1 to 10.

13. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the one-dimensional hydrodynamic modeling method for irrigation areas when cross-sectional data is insufficient, as described in any one of claims 1 to 10.

Images