Design face storm surge calculation method and device

Abstract

The application discloses a design area rainstorm calculation method and device, and the method comprises the following steps: obtaining an area rainfall Cv value according to characteristic parameters of annual maximum period rainfall of each rainfall station in a basin F ; obtaining a mean value EX F and a skewness coefficient Cs F of the area rainfall according to a mean value EX0 and a skewness coefficient CS0 of the annual maximum period rainfall of the rainfall station; and determining annual maximum design rainstorm of each period according to the area rainfall Cv F value, the mean value EX F and the skewness coefficient Cs F of the area rainfall. The application can obtain the design area rainfall by the rainfall data of the representative station.

Description

Technical Field

[0001] This invention relates to the field of computer technology, and in particular to a method and apparatus for calculating design surface storm surge. Background Technology

[0002] In water conservancy engineering design, using rainfall data to deduce the design rainfall, and then using the design rainfall to deduce the design flood, is a commonly used method for calculating the design flood. This method is particularly effective in watersheds where flow data is scarce or where human activity is significant. However, calculating the design flood requires the design rainfall process of the watershed, and there are generally two methods for calculating the design areal rainfall: the "direct method" and the "indirect method." The "direct method" uses concurrent observation data from various rain gauge stations, employing the arithmetic mean method or area weighting method to calculate the daily areal rainfall for each year, and then selects the areal rainfall with the maximum for each time period Δt each year. Then, the frequency calculation method is used to deduce the design areal rainfall XFP for each corresponding time period. The "indirect method" first calculates the design point rainfall Xop of the maximum for each time period Δt at the central rain gauge station of the watershed, and then converts the design point rainfall into the design areal rainfall XFP using the established point-area relationship. The former method has high data requirements, demanding a long period of concurrent observation data from rain gauge stations in the calculated watershed, which is generally difficult to meet. Furthermore, this method involves a large computational workload, but its results are highly accurate. The latter method has weaker point-area relationships in large and medium-sized watersheds, resulting in significant errors when converting design point rainfall to design areal rainfall. In addition, the latter method cannot fully utilize a large amount of measured heavy rainfall data. Summary of the Invention

[0003] The purpose of this invention is to provide a method and apparatus for calculating design surface storm surge, aiming to solve the above-mentioned problems in the prior art.

[0004] This invention provides a method for calculating design surface rainfall, comprising:

[0005] The areal rainfall Cv was estimated based on the characteristic parameters of the annual maximum rainfall during the period at each rain gauge station within the watershed. F value;

[0006] Based on the mean EX0 of the annual maximum rainfall during the period represented by the rain gauge station and the skewness coefficient CS0, the mean of the areal rainfall EX is estimated. F and skewness coefficient Cs F ;

[0007] Based on the areal rainfall Cv F Value, mean of areal rainfall EX F and skewness coefficient Cs F Determine the maximum annual design rainfall for each period.

[0008] This invention provides a design surface stormfall estimation calculation device, comprising:

[0009] The first estimation module is configured to estimate the areal rainfall Cv according to the characteristic parameters of the annual maximum period rainfall of each rainfall station in the basin F ;

[0010] The second estimation module is configured to estimate the mean value EX of the areal rainfall according to the mean value EX0 and the skewness coefficient CS0 of the annual maximum period rainfall of the rainfall station F ; F ;

[0011] The calculation module is configured to determine the annual maximum design storm of each period according to the Cv value of the areal rainfall, the mean value EX of the areal rainfall and the skewness coefficient CS F ; F ; F .

[0012] The embodiment of the present application also provides a design areal storm derivation calculation device, which comprises a memory, a processor and a computer program stored in the memory and executable on the processor, and the computer program is used for realizing the steps of the design areal storm derivation calculation method.

[0013] The embodiment of the present application also provides a computer readable storage medium, which stores an implementation program of information transmission, and the program is used for realizing the steps of the design areal storm derivation calculation method when executed by a processor.

[0014] By using the embodiment of the present application, the statistical relationship between point rainfall and areal rainfall is used to find a method for calculating the statistical parameters of areal rainfall by using the statistical parameters of point rainfall, so as to derive the design areal rainfall from the rainfall data of representative stations.

[0015] The above description is only a summary of the technical scheme of the present application, in order to more clearly understand the technical means of the present application, the content of the specification can be implemented, and in order to make the above and other purposes, characteristics and advantages of the present application more obvious and easy to understand, the specific embodiments of the present application are described as follows. BRIEF DESCRIPTION OF DRAWINGS

[0016] In order to more clearly illustrate the specific embodiments of the present application or the technical scheme in the prior art, the following will briefly introduce the drawings needed to be used in the specific embodiments or the prior art description. Obviously, the drawings in the following description are some embodiments of the present application, and those skilled in the art can also obtain other drawings according to these drawings without creating any creative labor.

[0017] Figure 1 is the flow chart of the design areal storm derivation calculation method of the embodiment of the present application;

[0018] Figure 2 is a schematic diagram of a design surface storm rainfall calculation device of the first embodiment of the device of the present application;

[0019] Figure 3 is a schematic diagram of a design surface storm rainfall calculation device of the second embodiment of the device of the present application. DETAILED DESCRIPTION

[0020] The technical solutions of the present application will be described below in conjunction with the embodiments, obviously, the described embodiments are part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0021] In the description of the present application, it should be understood that the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise" and the like indicate the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation of the present application.

[0022] In addition, the terms "first", "second" are only for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, the features defined with "first", "second" can explicitly or implicitly include one or more of the features. In the description of the present application, the meaning of "multiple" is two or more, unless otherwise specifically limited. In addition, the terms "mounting", "connecting", "connecting" should be broadly understood, for example, it can be fixedly connected, or it can be detachably connected, or integrally connected; it can be mechanically connected, or it can be electrically connected; it can be directly connected, or it can be indirectly connected through an intermediate medium, or it can be the communication between two elements. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0023] Method embodiment

[0024] According to the embodiment of the present application, a design surface storm rainfall calculation method is provided, Figure 1 is a flowchart of the design surface storm rainfall calculation method of the embodiment of the present application, as Figure 1 shown, the design surface storm rainfall calculation method according to the embodiment of the present application specifically includes:

[0025] Step 101: Estimate the areal rainfall Cv based on the characteristic parameters of the annual maximum rainfall period at each rain gauge station within the watershed. F Value; Step 101 specifically includes:

[0026] Assuming that the rainfall at each rain gauge station within the watershed is independent, the areal rainfall Cv can be estimated using Formula 1. F value:

[0027]

[0028] in, EX0 is the deviation coefficient representing the maximum annual rainfall at a rain gauge station. 2 / (EX0) 2 The maximum annual periodic isal rainfall is calculated based on the rainfall from the corresponding representative rain gauge station, where m is the number of rain gauge stations in the watershed.

[0029] When the rainfall data at each rain gauge station within the watershed are not independent of each other, the areal rainfall Cv can be estimated using Formula 2. F value:

[0030]

[0031] Where i and j represent different rain gauge stations, i ≠ j. The annual maximum hourly areal rainfall is calculated based on the annual maximum hourly areal rainfall at each corresponding station.

[0032] Step 102: Based on the mean EX0 of the annual maximum rainfall during the period represented by the rain gauge station and the skewness coefficient CS0, estimate the mean of the areal rainfall EX. F and skewness coefficient Cs F ;

[0033] Step 103, based on the areal rainfall Cv F Value, mean of areal rainfall EX F and skewness coefficient Cs F Determine the maximum annual design rainfall for each time period. Step 103 specifically includes:

[0034] According to Cs F Determine Φ p , where Φ p The deviation coefficient of a type P-III distribution with frequency P is represented by the deviation coefficient from the mean.

[0035] The annual maximum design rainfall X for each time period is determined according to Formula 3. p :

[0036] X p =EX F (Cv F Φ p +1) Formula 3.

[0037] The technical solutions of the embodiments of the present application are described in detail below.

[0038] The annual maximum storm of each period can be considered to obey P-III distribution, and the P-III distribution has three parameters, mean EX, coefficient of dispersion C v and coefficient of skewness C s . After the three parameters are determined, the annual maximum design storm of each period can be obtained by the following formula:

[0039] X p = EX(C v Φ p +1) (1)

[0040] In the formula, EX represents the mean, C v represents the coefficient of dispersion, and Φ p represents the coefficient of dispersion of P-III distribution with frequency P, which is only related to P and Cs.

[0041] It can be seen from formula (1) that X p is uniquely determined as long as EX, C v and C s are determined. In the prior art, it is considered in practice that the mean EX F and the coefficient of skewness Cs F of the annual maximum period rainfall are independent of the area change, and the mean EX0 and the coefficient of skewness Cs0 of the representative point storm annual maximum period rainfall can be selected to replace them. The embodiment of the present application focuses on the estimation of the coefficient of dispersion of the area rainfall.

[0042] (2) Estimation method of the coefficient of dispersion Cv F of the area rainfall:

[0043] The coefficient of dispersion Cv F of the area rainfall can be represented by the following formula:

[0044]

[0045] In the formula, σ F represents the mean square deviation of the annual maximum period area rainfall.

[0046]

[0047]

[0048] Wherein: i≠j.

[0049] It can be considered that the second-order origin moment of the corresponding annual maximum period area rainfall and the rainfall of the representative station in the same period is equal to the mean of the second-order origin moments of the rainfall of all stations; and the square of the mean of the rainfall of the representative station is equal to the mean of the squares of the means of the rainfall of all stations in the basin. Namely:

[0050]

[0051]

[0052] where m is the number of rainfall stations in the basin; EX0 is the average of the rainfall of the representative station in the corresponding year; EXi is the average of the rainfall of the i station in the corresponding year; EX0i is the second-order origin moment of the rainfall of the representative station in the corresponding year; EXii is the second-order origin moment of the rainfall of the i station in the corresponding year. i 2 i 2

[0053] If the rainfall stations in the basin are independent of each other, formula (4) can be expressed as:

[0054]

[0055] If the rainfall stations in the basin are not independent of each other, formula (4) can be expressed as:

[0056]

[0057] where i≠j, r ij is the correlation coefficient of the rainfall of the i station and the j station; σ i is the mean square deviation of the rainfall of the i station.

[0058] After formulae (2) to (8) are combined and arranged, the following two cases can be obtained.

[0059] If the rainfall stations in the basin are independent of each other, then:

[0060]

[0061] where CV0 is the coefficient of variation of the annual maximum period rainfall of the representative station in the corresponding year; EX0 2 / (EX0) 2 is calculated based on the rainfall of the representative station in the corresponding year, and the annual maximum period rainfall of the representative station can be used to approximate the calculation.

[0062] If the rainfall stations in the basin are not independent of each other, then:

[0063]

[0064] where i≠j.

[0065] where: is calculated based on the rainfall of each station in the corresponding year, and the annual maximum period rainfall of each station in the basin can be used to approximate the calculation. ​​​

[0066] In the embodiment of the present application, the surface rainfall Cv can be seen from the formula (9) and (10) F The value can be estimated by the characteristic parameters of the annual maximum period rainfall of each rainfall station in the basin, and the mean EX and the skewness coefficient Cs of the surface rainfall F and the skewness coefficient Cs F can be estimated by the mean EX0 and the skewness coefficient CS0 of the annual maximum period rainfall of the representative station. The parameters and the design value of the surface rainfall estimated according to the formula derived in this paper make the calculation more convenient and fast, and the data can be fully applied.

[0067] Example 1

[0068] The Nanshui Reservoir in Guangdong is located on the Nanshui River, a tributary of the Beijiang River in Ruyuan Yao Autonomous County, Guangdong. The control basin area is 608 km2. The reservoir was impounded in February 1967, and the total storage capacity is 12.805×108m3, belonging to a multi-year regulation reservoir. There are four rainfall stations with long observation data in the reservoir, which are Dam station, Tixia station, Baitu station and Pingxi station. Due to the need of flood review of the reservoir, two methods are used to calculate the design storm of the Nanshui Reservoir: one method is a direct method to calculate the annual surface average daily rainfall of each station from 1961 to 1998 by using the arithmetic mean method, and then to calculate the surface rainfall of the annual maximum period Δt (Δt = 1, 3, 5, 7d), and then to calculate the design surface rainfall; the second method is to select a representative station, select Tixia station as the representative station, calculate the rainfall parameters and the correlation coefficient of each rainfall station in the annual maximum period Δt (Δt = 1, 3, 5, 7d), and then substitute into formula (10) to calculate the surface rainfall of each period, and then to calculate the design surface rainfall of the corresponding period.

[0069] The calculation results are shown in Table 1.

[0070] Table 1 Calculation results of design surface rainfall in annual maximum period

[0071]

[0072] Table 2 Correlation coefficient table of annual maximum rainfall of each station

[0073]

[0074] Among them: the CV of Tixia station is respectively: 1 day 0.6, 3 days 0.5, 5 days 0.52, 7 days 0.47. As can be seen from Table 1, the surface rainfall estimated by the representative station adopts the technical scheme of the embodiment of the present application, which can accurately estimate the surface rainfall. The technical scheme of the embodiment of the present application provides a simple, convenient and high-precision method for deriving the design surface rainfall.

[0075] In summary, the embodiment of the present application seeks a method for calculating the statistical parameter of the surface rainfall by the statistical relationship between the point rainfall and the surface rainfall, so as to achieve the purpose of calculating the design surface rainfall by the rainfall data of the representative station.

[0076] Device embodiment one

[0077] According to the embodiment of the present application, a design surface storm calculation device is provided. Figure 2 is a design surface storm calculation device of the device embodiment one of the present application, as shown in Figure 2 According to the embodiment of the present application, the design surface storm calculation device specifically comprises:

[0078] The first estimation module 20 is used for estimating the surface rainfall Cv value according to the annual maximum period rainfall characteristic parameters of each rainfall station in the basin. F The first estimation module 20 is specifically used for:

[0079] In the case that the rainfall of each rainfall station in the basin is independent, the surface rainfall Cv value is estimated according to the formula 1: F

[0080]

[0081] Among them, is the dispersion coefficient of the annual maximum period rainfall of the representative rainfall station, EX0 2 / (EX0) 2 is the annual maximum period surface rainfall calculated according to the rainfall of the corresponding representative rainfall station, and m is the number of rainfall stations in the basin.

[0082] In the case that the rainfall of each rainfall station in the basin is not independent, the surface rainfall Cv value is estimated according to the formula 2: F

[0083]

[0084] Among them, i and j represent different rainfall stations, i≠j, is the annual maximum period surface rainfall calculated according to the annual maximum period surface rainfall of each corresponding station.

[0085] The second estimation module 22 is used for estimating the mean value EX of the surface rainfall and the skewness coefficient Cs according to the mean value EX0 of the annual maximum period rainfall of the representative rainfall station and the skewness coefficient CS0. F F

[0086] The calculation module 24 is used for calculating the mean value EX of the surface rainfall and the skewness coefficient Cs according to the surface rainfall Cv value, the mean value EX of the surface rainfall and the skewness coefficient Cs F F F ​​​​​​Determine the annual maximum design storm of each period. The calculation module 24 is specifically used for:

[0087] According to Cs F Determine Φ p , wherein Φ p represents the coefficient of variation of P-III type distribution with frequency P;

[0088] Determine the annual maximum design storm X of each period according to formula 3 p :

[0089] X p = EX F (Cv F Φ p +1) Formula 3.

[0090] The embodiment of the application is a device corresponding to the above-mentioned method embodiment, and the specific operation of each module can be understood with reference to the description of the method embodiment, which will not be repeated here.

[0091] Device embodiment two

[0092] The embodiment of the application provides a design surface storm derivation calculation device, as shown in the figure, comprising: a memory 30, a processor 32 and a computer program stored on the memory 30 and executable on the processor 32, the computer program is executed by the processor 32 to realize the steps as described in the method embodiment. Figure 3

[0093] Device embodiment two

[0094] The embodiment of the application provides a computer readable storage medium, and the computer readable storage medium stores an implementation program of information transmission, and the program is executed by the processor 32 to realize the steps as described in the method embodiment.

[0095] The computer readable storage medium described in the embodiment includes but is not limited to ROM, RAM, magnetic disk or optical disk and the like.

[0096] The above describes specific embodiments of the present specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be executed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multi-task processing and parallel processing are possible or can be advantageous.

[0097] ​In the 1930s, it was clear to distinguish whether an improvement in a technology was in hardware (e.g., improvement in circuit structure of diodes, transistors, switches, etc.) or in software (e.g., improvement in method flow). However, as technology has evolved, many improvements in method flow today can be considered as direct improvements in hardware circuit structure. Designers almost always obtain the corresponding hardware circuit structure by programming the improved method flow into the hardware circuit. Therefore, it cannot be said that an improvement in a method flow cannot be implemented by a hardware entity module. For example, a programmable logic device (PLD) (e.g., a field programmable gate array (FPGA)) is an integrated circuit whose logic function is determined by user programming of the device. A digital system is "integrated" on a PLD by the designer programming it, rather than by asking a chip manufacturer to design and fabricate a custom integrated circuit chip. Moreover, instead of manually fabricating an integrated circuit chip, this programming is now mostly implemented by "logic compiler" software, which is similar to software compilers used in program development, and the original code to be compiled is written in a specific programming language, called a hardware description language (HDL), of which there are many, such as ABEL (Advanced Boolean Expression Language), AHDL (Altera Hardware Description Language), Confluence, CUPL (Cornell University Programming Language), HDCal, JHDL (Java Hardware Description Language), Lava, Lola, MyHDL, PALASM, RHDL (Ruby Hardware Description Language), etc., the most commonly used being VHDL (Very-High-Speed Integrated Circuit Hardware Description Language) and Verilog. It should be clear to those skilled in the art that, by simply logically programming a method flow in one of the above hardware description languages and programming it into an integrated circuit, a hardware circuit implementing the logical method flow can be easily obtained.

[0098] The controller can be implemented in any suitable way, for example, the controller can take the form of a microprocessor or processor and a computer readable medium storing computer readable program code, e.g. software or firmware, executable by the (micro)processor, logic gates, switches, an application specific integrated circuit (ASIC), a programmable logic controller and an embedded microcontroller, examples of which include but are not limited to the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320, the memory controller can also be implemented as part of the control logic of the memory. Those skilled in the art will also know that, in addition to being implemented in pure computer readable program code, the controller can equally well be implemented to perform the same functions using logic gates, switches, an application specific integrated circuit, a programmable logic controller and an embedded microcontroller, etc. by means of a logical programming of the method steps. The controller can thus be considered as a hardware component, and the means comprised therein for performing the various functions can be considered as structures within the hardware component. Alternatively, the means for performing the various functions can even be considered as both a software module implementing the method and a structure within the hardware component.

[0099] The systems, apparatuses, modules or units illustrated by the above embodiments can be implemented by computer chips or entities, or by products with certain functions. A typical implementation device is a computer. Specifically, the computer can be a personal computer, a laptop computer, a cellular phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

[0100] For the sake of description, the above apparatuses are described in various units with functions respectively. Of course, the functions of the units can be implemented in one or more software and / or hardware in implementing the embodiments of the present specification.

[0101] Those skilled in the art will appreciate that one or more embodiments of the present specification can be provided as a method, a system or a computer program product. Therefore, one or more embodiments of the present specification can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present specification can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROMs, optical storage devices, etc.) containing computer usable program code.

[0102] The specification is presented with reference to flow diagrams and / or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the specification. It will be understood that each block of the flow diagrams and / or block diagrams, and combinations of blocks in the flow diagrams and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processing element or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flow diagrams and / or block diagrams block or blocks. Figure 1 The flow diagram and / or block diagram in the flow diagrams and / or block diagrams can represent one or more of any appropriate circuitry configured to perform the specified functions. In this regard, one or more flow diagrams and / or block diagrams in the flow diagrams and / or block diagrams can represent a device or devices configured to perform one or more of the functions described herein. Figure 1 The flow diagram and / or block diagram in the flow diagrams and / or block diagrams can represent one or more of any appropriate circuitry configured to perform the specified functions. In this regard, one or more flow diagrams and / or block diagrams in the flow diagrams and / or block diagrams can represent a device or devices configured to perform one or more of the functions described herein.

[0103] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flow diagrams and / or block diagrams block or blocks. Figure 1 The flow diagram and / or block diagram in the flow diagrams and / or block diagrams can represent one or more of any appropriate circuitry configured to perform the specified functions. In this regard, one or more flow diagrams and / or block diagrams in the flow diagrams and / or block diagrams can represent a device or devices configured to perform one or more of the functions described herein. Figure 1 The flow diagram and / or block diagram in the flow diagrams and / or block diagrams can represent one or more of any appropriate circuitry configured to perform the specified functions. In this regard, one or more flow diagrams and / or block diagrams in the flow diagrams and / or block diagrams can represent a device or devices configured to perform one or more of the functions described herein.

[0104] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flow diagrams and / or block diagrams block or blocks. Figure 1 The flow diagram and / or block diagram in the flow diagrams and / or block diagrams can represent one or more of any appropriate circuitry configured to perform the specified functions. In this regard, one or more flow diagrams and / or block diagrams in the flow diagrams and / or block diagrams can represent a device or devices configured to perform one or more of the functions described herein. Figure 1 The flow diagram and / or block diagram in the flow diagrams and / or block diagrams can represent one or more of any appropriate circuitry configured to perform the specified functions. In this regard, one or more flow diagrams and / or block diagrams in the flow diagrams and / or block diagrams can represent a device or devices configured to perform one or more of the functions described herein.

[0105] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0106] The memory can include non-persistent memory and / or persistent memory, such as flash memory, read-only memory (ROM), and / or volatile or non-volatile random access memory (RAM), among others. The memory is an example of computer-readable media.

[0107] Computer-readable media includes permanent and non-permanent, movable and non-movable media that can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible to a computing device. According to the definition herein, computer-readable media does not include transitory media such as modulated data signals and carriers.

[0108] It should also be noted that the terms "comprising", "containing", or any other variant thereof are intended to cover non-exclusive inclusion, such that processes, methods, articles or devices that include a series of elements not only include those elements, but also include other elements not explicitly listed, or inherent to such processes, methods, articles or devices. Without more limitations, the element defined by the statement "comprising a" does not exclude the presence of additional identical elements in the process, method, article or device that includes the element.

[0109] One or more embodiments of the present specification can be described in the general context of computer-executable instructions being executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, and the like that perform particular tasks or implement particular abstract data types. One or more embodiments of the present specification can also be practiced in a distributed computing environment, in which tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules can be located in both local and remote computer storage media, including storage devices.

[0110] Each embodiment in the present specification is described in a progressive manner, and the same or similar parts between each embodiment can be referred to each other, and each embodiment focuses on the difference from other embodiments. In particular, for system embodiments, since they are basically similar to method embodiments, the description is relatively simple, and the relevant parts can be referred to the part of the method embodiment.

[0111] The above merely provides the examples of the present document and is not intended to limit the present document. For those skilled in the art, the present document can have various modifications and changes. Any modification, equivalent replacement, improvement, etc. within the spirit and principle of the present document shall be included in the scope of claims of the present document.

Claims

1. A design storm hydrograph calculation method, characterized by, comprising: The areal rainfall is estimated according to the characteristic parameters of the annual maximum period rainfall of each rainfall station in the basin values; specifically including: In the case of independent rainfall at each rainfall station in a drainage basin, the areal rainfall is estimated according to Formula 1 Value: Formula 1 ; wherein, is the exponent coefficient representing the annual maximum hourly rainfall of the rainfall station, is the annual maximum hourly areal rainfall calculated according to the rainfall of the corresponding representative rainfall station, and m is the number of rainfall stations in the basin. According to the mean and the coefficient of skewness representing the annual maximum period rainfall of the rainfall station and the coefficient of skewness , the mean and the coefficient of skewness of the areal rainfall are estimated and the coefficient of skewness ; specifically comprising: In the case of the rainfall at each rainfall station in a river basin is not independent, the areal rainfall is estimated according to Formula 2 Value: Formula 2; where i and j represent different rain gauges, , is the annual maximum period area rainfall calculated from the annual maximum period area rainfall at each respective station. According to the planar rainfall intensity The mean value of the planar rainfall intensity And the coefficient of skewness Determine the annual maximum design storm of each period; specifically including: According to determined Φ p wherein Φ p represents the coefficient of variation of the P-III distribution at frequency P; The annual maximum design storm X for each period is determined according to Formula 3 p : X p = ( Φ p +1) Formula 3.

2. A design storm hydrograph calculation device, characterized by, comprising: The first estimation module is configured to estimate the areal rainfall according to the annual maximum time period rainfall characteristic parameters of each rainfall station in the basin values; in particular for: In the case of independent rainfall at each rainfall station in a drainage basin, the areal rainfall is estimated according to Formula 1 Value: Formula 1 ; wherein, is the representative coefficient of the annual maximum hourly rainfall of the rainfall station, is the annual maximum hourly areal rainfall calculated from the rainfall of the corresponding representative rainfall station, and m is the number of rainfall stations in the basin. a second estimation module configured to estimate the mean and the coefficient of skewness of the areal rainfall according to the mean of the annual maximum period rainfall of the rain gauges and the coefficient of skewness of the areal rainfall and the coefficient of skewness ; in particular for: In the case of the rainfall at each rainfall station in a river basin is not independent, the areal rainfall is estimated according to Formula 2 Value: Formula 2: where i and j represent different rain gauges, , is the annual maximum period area rainfall calculated from the annual maximum period area rainfall of each respective station. A computing module is configured to determine the annual maximum design storm of each period according to the face rain amount , the mean value of the face rain amount , and the skewness coefficient ; in particular for: According to determined Φ p wherein Φ p represents the coefficient of variation of the P-III distribution at frequency P; The annual maximum design storm X for each period is determined according to Formula 3 p : X p = ( Φ p +1) Formula 3.

3. A design storm hydrograph calculation device, characterized by, a memory, a processor, and a computer program stored on the memory and executable on the processor, which, when executed by the processor, implements the steps of the design surface storm surge calculation method according to claim 1. The computer readable storage medium stores an information transmission implementation program, and the program, when executed by a processor, implements the steps of the design surface storm surge calculation method according to claim 1.

4. A computer-readable storage medium, characterized in that, ​

Images