Method and device for determining height of surge chamber of hydropower station and electronic equipment

Abstract

The application provides a method and device for determining the height of a surge chamber of a hydropower station and an electronic device. The method is based on the head loss of the surge chamber of the target hydropower station, the head loss of the diversion tunnel of the target hydropower station, and the flow data in the main pipe of the pressure pipeline corresponding to the current working condition of the target hydropower station. The method is processed by a surge chamber surge fusion algorithm, comprehensively considers the influence of the water flow inertia in the surge chamber connecting pipe and the head loss, and can make the obtained surge chamber surge determination algorithm more accurate. Then, according to the derivative of the surge extreme value with respect to time being a preset derivative value, the surge chamber surge determination algorithm is processed to obtain the surge extreme value occurrence time. Therefore, by using the surge extreme value occurrence time and the surge initial parameter corresponding to the current working condition, the surge chamber surge determination algorithm is processed, a more accurate surge extreme value can be obtained, and then the height of the surge chamber can be reasonably determined while ensuring safety according to the surge extreme value.

Description

Technical Field

[0001] This application relates to the field of water conservancy and hydropower engineering technology, and in particular to a method, device and electronic equipment for determining the height of a surge tank in a hydropower station. Background Technology

[0002] Hydropower station turbines require a stable water flow during operation; excessively large or small flows can damage the hydraulic machinery. In particular, excessively large flows can cause catastrophic accidents such as turbine vibration, collisions, or contact failures, leading to overload or destructive oscillations. Therefore, surge tanks are needed to regulate the water flow and stabilize the pressure. However, if the height of the surge tank is improperly designed, the safe operation of the hydraulic machinery cannot be guaranteed.

[0003] Therefore, determining the appropriate height of the pressure regulating chamber while ensuring safety has become a pressing technical problem. Summary of the Invention

[0004] In view of this, the purpose of this application is to provide a method, device and electronic equipment for determining the height of the surge tank of a hydropower station, so as to solve or partially solve the above-mentioned technical problems.

[0005] To achieve the above objectives, the first aspect of this application provides a method for determining the height of a surge chamber in a hydropower station, comprising:

[0006] Determine the head loss of the surge tank of the target hydropower station, and determine the head loss of the water diversion tunnel of the target hydropower station;

[0007] Obtain the flow rate data in the main pressure pipeline corresponding to the current operating condition of the target hydropower station. Based on the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline, process them through the surge chamber surge fusion algorithm to obtain the surge chamber surge determination algorithm.

[0008] Based on the derivative of the surge extreme value with respect to time as a preset derivative value, the surge extreme value is processed according to the surge determination algorithm of the pressure regulating chamber to obtain the surge extreme value occurrence time;

[0009] Determine the initial surge parameters corresponding to the current operating condition, process the surge extreme value and the surge initial parameters through the surge determination algorithm of the surge chamber to obtain the surge extreme value, and determine the height of the surge chamber based on the surge extreme value.

[0010] Based on the same inventive concept, a second aspect of this application provides a device for determining the height of a surge tank in a hydropower station, comprising:

[0011] The head loss determination module is configured to determine the head loss of the surge tank of the target hydropower station, and to determine;

[0012] The fusion processing module is configured to acquire the flow data in the main pressure pipeline corresponding to the current operating condition of the target hydropower station, and process the head loss of the surge chamber, the head loss of the water diversion tunnel and the flow data in the main pressure pipeline through the surge chamber fusion algorithm to obtain the surge chamber determination algorithm.

[0013] The time determination module is configured to process the surge determination algorithm based on the derivative of the surge extreme value with respect to time as a preset derivative value to obtain the time of the surge extreme value occurrence.

[0014] The height determination module is configured to determine the initial parameters of the surge corresponding to the current operating condition, process the surge extreme value and the initial parameters of the surge using the surge extreme value and the surge determination algorithm to obtain the surge extreme value, and determine the height of the surge chamber based on the surge extreme value.

[0015] Based on the same inventive concept, a third aspect of this application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable by the processor, wherein the processor, when executing the computer program, implements the method described in the first aspect above.

[0016] As can be seen from the above, the method, device, and electronic equipment for determining the height of a surge chamber in a hydropower station provided in this application determine the head loss of the surge chamber and the head loss of the water diversion tunnel of the target hydropower station. Then, based on the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow data in the main pressure pipeline corresponding to the current operating conditions of the target hydropower station, the surge chamber surge fusion algorithm is used to comprehensively consider the influence of the inertia of the water flow in the surge chamber connecting pipe and the head loss, thereby making the surge chamber surge determination algorithm more accurate. Then, based on the derivative of the surge extreme value with respect to time as a preset derivative value, the surge extreme value occurrence time is obtained by processing the surge extreme value using the surge extreme value occurrence time and the surge initial parameters corresponding to the current operating conditions through the surge chamber surge determination algorithm. Therefore, based on the surge extreme value, the height of the surge chamber can be reasonably determined while ensuring safety. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in this application or related technologies, the drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the drawings described below are only embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0018] Figure 1 This is a flowchart of a method for determining the height of a surge chamber in a hydropower station, as described in an embodiment of this application.

[0019] Figure 2A This is a schematic diagram illustrating the process of determining the height of the surge tank in a hydropower station according to an embodiment of this application.

[0020] Figure 2B This is a schematic diagram of the layout of a hydropower station according to an embodiment of this application;

[0021] Figure 3 This is a structural block diagram of the surge tank height determination device of a hydropower station according to an embodiment of this application;

[0022] Figure 4 This is a schematic diagram of an electronic device according to an embodiment of this application. Detailed Implementation

[0023] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with specific embodiments and the accompanying drawings.

[0024] It should be noted that, unless otherwise defined, the technical or scientific terms used in the embodiments of this application should have the ordinary meaning understood by one of ordinary skill in the art to which this application pertains. The terms "first," "second," and similar terms used in the embodiments of this application do not indicate any order, quantity, or importance, but are merely used to distinguish different components. Terms such as "comprising" or "including" mean that the element or object preceding the word encompasses the elements or objects listed after the word and their equivalents, without excluding other elements or objects. Terms such as "connected" or "linked" are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. Terms such as "upper," "lower," "left," and "right" are only used to indicate relative positional relationships; when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0025] It is understood that before using the technical solutions of the various embodiments in this application, users will be informed of the type, scope of use, and usage scenarios of the personal information involved in an appropriate manner, and user authorization will be obtained.

[0026] For example, upon receiving a user's active request, a prompt message is sent to the user to explicitly inform them that the requested operation will require the acquisition and use of the user's personal information. This allows the user to independently choose, based on the prompt message, whether to provide personal information to the software or hardware such as electronic devices, applications, servers, or storage media performing the operations described in this application.

[0027] As an optional but not limited implementation, in response to a user's active request, sending a prompt message to the user can be done via a pop-up window, where the prompt message can be presented in text format. Furthermore, the pop-up window can also include a selection control allowing the user to choose "agree" or "disagree" to provide personal information to the electronic device.

[0028] It is understood that the above notification and user authorization process is merely illustrative and does not limit the implementation of this application. Other methods that comply with relevant laws and regulations may also be applied to the implementation of this application.

[0029] In hydropower stations with long water diversion systems, steep slopes, and low unit installation elevations, extending the short pipe of the impedance-type surge tank—that is, constructing a long connecting pipe between the surge tank shaft and the tunnel—reduces the amount of mountain excavation and thus lowers construction costs. When the connecting pipe is short, its impact can be ignored, and it is still treated as a conventional impedance-type surge tank in calculations and analyses. However, when the connecting pipe is long, significant water flow inertia and head loss will occur within the connecting pipe during the water flow into and out of the surge tank. If this impact is ignored, it will lead to certain calculation deviations, affecting the accuracy of engineering design and even causing power station safety accidents.

[0030] In addition, existing technologies estimate the extreme values ​​of surge waves in the surge chamber during the hydropower station design phase using the following method, with the calculation formula as follows:

[0031] When the operating condition is load shedding:

[0032]

[0033]

[0034]

[0035]

[0036] When the operating condition is increased load:

[0037]

[0038]

[0039]

[0040]

[0041]

[0042] in, This represents the head loss caused by the total flow into and out of the upstream surge tank. This indicates the head loss within the entire water diversion tunnel. Indicates the length of the water diversion tunnel. Indicates the area of ​​the water diversion tunnel. Indicates the flow velocity in the water diversion tunnel. This indicates the area of ​​the upstream surge tank well.

[0043] However, it does not consider the effects of water flow inertia and head loss within the connecting pipes, resulting in inaccurate calculations. Furthermore, this formula is not applicable to any single-tunnel, multi-machine hydropower station; calculations for different layouts require manual conversion by designers, leading to low computational efficiency.

[0044] In addition, existing technologies also provide a calculation model for surge waves in upstream surge chambers with connecting pipes based on the method of characteristics. The calculation formula is as follows:

[0045] Bottom branch pipe:

[0046]

[0047]

[0048] Connecting pipe:

[0049]

[0050]

[0051]

[0052] Connection point between the connecting pipe and the well:

[0053]

[0054] Pressure regulating chamber shaft:

[0055]

[0056]

[0057]

[0058] in, Indicates the flow rate of the water diversion tunnel. Indicates the flow rate within the main pressure pipeline. Indicates the flow rate inside the connecting pipe. , , This indicates the piezometric head at the bifurcation point of each pipeline. This indicates the inflow and outflow flow rates of the surge tank's main well. This indicates the water head in the piezometer at the bottom of the well. Indicates surge in the pressure regulating chamber. Represents the cross-sectional area of ​​the large well. Indicates the elevation of the pressure regulating chamber floor. This indicates the head loss of the pressure regulating chamber through the orifice. This represents the impedance loss coefficient.

[0059] However, this method requires complete technical parameters for both the water conveyance system and the electromechanical system when using the method of characteristics to solve the water conveyance system of a power station. Furthermore, the method is complex to program and involves a large amount of computation. In the initial design phase of a power station, data is often incomplete, and the method of characteristics cannot meet the speed requirements of the initial design phase.

[0060] This application provides a method for determining the height of a surge chamber in a hydropower station. The method determines the head loss of the surge chamber and the head loss of the water diversion tunnel of the target hydropower station. Then, based on the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline corresponding to the current operating conditions of the target hydropower station, a surge chamber surge fusion algorithm is used. This algorithm comprehensively considers the influence of water flow inertia and head loss in the surge chamber connecting pipe, thereby making the surge chamber surge determination algorithm more accurate. The surge extreme value is then processed using the surge extreme value with respect to time as a preset derivative value, and the surge extreme value occurrence time is obtained. Since this application does not ignore the influence of water flow inertia and head loss in the surge chamber connecting pipe, using the surge extreme value occurrence time and the initial surge parameters corresponding to the current operating conditions through the surge chamber surge determination algorithm, a more accurate surge extreme value can be obtained. Therefore, based on the surge extreme value, the height of the surge chamber can be reasonably determined while ensuring safety.

[0061] like Figure 1 As shown, the method in this embodiment includes:

[0062] Step 101: Determine the head loss of the surge tank of the target hydropower station and the head loss of the water diversion tunnel of the target hydropower station.

[0063] In this step, in order to better ensure the safety of the power station, the head loss of the surge tank of the target hydropower station and the head loss of the water diversion tunnel of the target hydropower station are determined to comprehensively consider the influence of the water flow inertia and head loss in the surge tank connecting pipe, thereby avoiding deviations that could affect the accuracy of the engineering design.

[0064] Step 102: Obtain the flow rate data in the main pressure pipeline corresponding to the current operating condition of the target hydropower station. Based on the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline, process them using the surge chamber surge fusion algorithm to obtain the surge chamber surge determination algorithm.

[0065] In this step, the surge wave fusion algorithm is used to process the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline corresponding to the current operating conditions of the target hydropower station. This comprehensively considers the influence of the water flow inertia and head loss in the connecting pipe of the surge chamber, thereby making the surge wave determination algorithm of the surge chamber more accurate.

[0066] Step 103: Based on the preset derivative value of the surge extreme value with respect to time, the surge extreme value is processed according to the surge determination algorithm of the pressure regulating chamber to obtain the surge extreme value occurrence time.

[0067] In this step, since the surge determination algorithm for the surge chamber is relatively accurate, the surge determination algorithm is used to process the surge extreme value based on the derivative of the surge extreme value with respect to time as a preset derivative value. This can obtain a more accurate surge extreme value occurrence time, reduce deviations, and ensure the safety of the power station.

[0068] Step 104: Determine the initial surge parameters corresponding to the current operating condition, process the surge extreme value and the initial surge parameters using the surge determination algorithm of the surge chamber to obtain the surge extreme value, and determine the height of the surge chamber based on the surge extreme value.

[0069] In this step, since this application does not ignore the influence of water flow inertia and head loss in the surge chamber connecting pipe, the surge extreme value is obtained by processing the surge extreme value through the surge determination algorithm with the surge extreme value and the surge initial parameters corresponding to the current working condition. Then, the height of the surge chamber can be reasonably determined based on the surge extreme value while ensuring safety.

[0070] The top elevation of the surge chamber is the maximum surge peak plus the safety freeboard, which is, for example, 2-3 meters. The bottom elevation of the surge chamber should be at least a safe distance from the minimum surge peak, which is, for example, 1 meter.

[0071] In addition, the minimum surge peak value may not be used as a reference value for the bottom elevation of the surge chamber, but it can be used as a reference value for the bottom elevation of the pressure steel pipe inlet of the target hydropower station. For example, the bottom elevation of the pressure steel pipe inlet should be at least 1 meter away from the minimum surge peak value.

[0072] The above scheme determines the head loss of the surge chamber and the head loss of the water diversion tunnel of the target hydropower station. Then, based on the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline corresponding to the current operating conditions of the target hydropower station, a surge fusion algorithm is used. This comprehensively considers the influence of water flow inertia and head loss in the surge chamber connecting pipe, thus making the surge determination algorithm more accurate. The surge extreme value is then processed using the surge determination algorithm based on a preset derivative value of the surge extreme value with respect to time to obtain the surge extreme value occurrence time. Since this application does not ignore the influence of water flow inertia and head loss in the surge chamber connecting pipe, using this surge extreme value occurrence time and the surge initial parameters corresponding to the current operating conditions through the surge determination algorithm yields a more accurate surge extreme value. Based on this surge extreme value, the height of the surge chamber can be reasonably determined while ensuring safety.

[0073] In some embodiments, step 101 includes:

[0074] Acquire surge data of the pressure regulating chamber The turbine flow rate of the target hydropower station is... The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is Head loss in the water diversion tunnel The water turbine's flow rate is The water flow velocity in the water diversion tunnel at that time The large well area of ​​the pressure regulating chamber The area inside the water diversion tunnel and the flow velocity inside the water diversion tunnel ;

[0075] Based on the surge data of the pressure regulating chamber The water turbine's flow rate is The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is Head loss in the water diversion tunnel The water turbine's flow rate is The water flow velocity in the water diversion tunnel at that time The large well area of ​​the pressure regulating chamber and the area inside the water diversion tunnel The head loss of the surge tank of the target hydropower station is determined by the following formula. :

[0076]

[0077] in, ;

[0078] Based on the preset head loss coefficient of the water diversion tunnel and the flow velocity inside the water diversion tunnel The head loss of the water diversion tunnel of the target hydropower station is determined by the following formula. :

[0079] .

[0080] In the above scheme, surge data from the pressure regulating chamber are comprehensively considered. The water turbine's flow rate is Head loss generated by the pressure regulating chamber The water turbine's flow rate is Head loss in the water diversion tunnel The water turbine's flow rate is Water flow velocity in the water diversion tunnel Large well area of ​​the pressure regulating chamber and the area inside the water diversion tunnel This can reduce the head loss in the surge tank of a standard hydroelectric power station. More accurate.

[0081] In addition, by comprehensively considering the pre-set head loss coefficient of the water diversion tunnel... and the flow velocity inside the water diversion tunnel This can reduce the head loss in the water diversion tunnel of the target hydropower station. More accurate.

[0082] In some embodiments, in step 102, the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline are processed using a surge chamber surge algorithm to obtain a surge chamber surge determination algorithm, including:

[0083] Step A1: Based on the head loss of the surge chamber and the head loss of the water diversion tunnel, the first nonlinear surge equation of the surge chamber is obtained by combining the pre-constructed continuity equation of the surge chamber with the pre-constructed momentum equation of the water diversion tunnel.

[0084] Step A2: Based on the flow rate data in the main pressure pipeline, obtain the nonlinear surge equation for the second pressure regulating chamber using the nonlinear surge equation of the first pressure regulating chamber.

[0085] Step A3: Using the pre-set surge parameters under the current operating condition, the surge equation of the surge chamber is obtained by combining the pre-constructed nonlinear vibration equation with the surge equation of the second surge chamber.

[0086] In the above scheme, the surge wave fusion algorithm is used to process the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow data in the main pressure pipeline corresponding to the current operating conditions of the target hydropower station. This comprehensively considers the influence of the inertia of the water flow in the connecting pipe of the surge chamber and the head loss, thereby making the surge wave determination algorithm of the surge chamber more accurate.

[0087] In some embodiments, step A1 includes:

[0088] The head loss in the pressure regulating chamber and the head loss of the water diversion tunnel Substituting the continuity equation of the surge tank and simultaneously solving the momentum equation of the water diversion tunnel, the flow velocity inside the water diversion tunnel is... After elimination processing, the nonlinear surge equation of the first pressure regulating chamber is obtained;

[0089] The continuity equation for the pressure regulating chamber is as follows:

[0090]

[0091] in, This indicates the length of the water diversion tunnel. This indicates the area of ​​the water diversion tunnel. Indicates the flow rate inside the water diversion tunnel. This indicates the head loss of the water diversion tunnel. This represents the surge data of the surge chamber. This indicates the length of the connecting pipe to the pressure regulating chamber. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the flow rate within the connecting pipe of the pressure regulating chamber. Represents gravitational acceleration. Indicates time, This indicates the head loss caused by the inflow and outflow from the pressure regulating chamber;

[0092] The momentum equation for the water diversion tunnel is as follows:

[0093]

[0094] in, This indicates the area of ​​the large well in the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, Indicates the flow rate inside the water diversion tunnel. This represents the flow rate data within the main pressure pipeline of the target hydropower station;

[0095] The nonlinear surge equation for the first pressure regulating chamber is as follows:

[0096]

[0097] in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, The turbine's flow rate is indicated as Water head loss in the water diversion tunnel This represents the flow rate data within the main pressure pipeline of the target hydropower station. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel mentioned above.

[0098] In the above scheme, based on the rigid water hammer theory and considering the inertia of the water flow and head loss in the surge tank connecting pipe, the continuity equation of the surge tank and the momentum equation of the water diversion tunnel can be obtained.

[0099] In the continuity equation of the pressure regulating chamber:

[0100]

[0101] The surge data for the upstream surge chamber are all defined as positive when upward, with the reference being the surge data from the upstream reservoir. This indicates the length of the connecting pipe to the upstream pressure regulating chamber. This indicates the area of ​​the connecting pipe in the upstream pressure regulating chamber. This indicates the flow rate in the connecting pipe of the upstream pressure regulating chamber. This indicates the head loss caused by water flowing into and out of the upstream surge tank. This indicates the length of the connecting pipe to the upstream pressure regulating chamber. This indicates the area of ​​the connecting pipe in the upstream pressure regulating chamber.

[0102] In the momentum equation of the water diversion tunnel:

[0103]

[0104] This indicates the area of ​​the large well in the upstream surge tank. This indicates the flow rate data within the main pressure pipeline upstream of the target hydropower station.

[0105] Taking into account the head loss of the surge tank and head loss of water diversion tunnel The head loss of the surge tank and head loss of water diversion tunnel Substituting the continuity equation of the surge tank and simultaneously solving the momentum equation of the water diversion tunnel, the flow velocity inside the water diversion tunnel is... Elimination processing can make the obtained nonlinear surge equation of the first pressure regulating chamber more accurate.

[0106] In some embodiments, the second sub-surge chamber nonlinear surge equation includes a first sub-surge chamber nonlinear surge equation and a second sub-surge chamber nonlinear surge equation.

[0107] Step A2 includes:

[0108] Step A21: In response to the current operating condition of the target hydropower station being load shedding, and the flow rate data in the main pressure pipeline corresponding to the load shedding being 0, the flow rate data in the main pressure pipeline of the target hydropower station is substituted into the nonlinear surge equation of the first surge chamber as the flow rate data in the main pressure pipeline of the target hydropower station. The nonlinear surge equation for the first sub-pressure regulating chamber is obtained.

[0109] The nonlinear surge equation for the first pressure regulating chamber is as follows:

[0110]

[0111] in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This represents the flow rate data within the main pressure pipeline of the target hydropower station. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: Head loss within the aforementioned water diversion tunnel;

[0112] The nonlinear surge equation for the first sub-pressure regulating chamber is as follows:

[0113]

[0114] in, This represents the surge data of the surge chamber. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity within the aforementioned water diversion tunnel. Or,

[0115] Step A22, in response to the current operating condition of the target hydropower station being an increased load, and the flow rate data in the main pressure pipeline corresponding to the increased load being... The flow rate data in the main pressure pipeline is... Substituting the flow data in the main pressure pipeline of the target hydropower station into the nonlinear surge equation of the first surge chamber... The nonlinear surge equation for the second sub-pressure regulating chamber is obtained.

[0116] The nonlinear surge equation for the first pressure regulating chamber is as follows:

[0117]

[0118] in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This represents the flow rate data within the main pressure pipeline of the target hydropower station. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: Head loss within the aforementioned water diversion tunnel;

[0119] The nonlinear surge equation for the second sub-pressure regulating chamber is as follows:

[0120]

[0121] in, This represents the simple harmonic solution of the pressure regulating chamber under increased load conditions. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time.

[0122] In the above scheme, when the target hydropower station is currently operating under load shedding conditions, it is assumed that the target hydropower station has a single-unit flow rate of [missing information]. A multi-machine hydropower station with one tunnel, when the same hydraulic unit The unit was operating normally. Due to a sudden load shedding caused by an accident, the flow rate in the main pressure pipeline (i.e., the flow rate data in the main pressure pipeline corresponding to the load shedding) changed from... The flow rate suddenly decreases to 0. Substituting the flow rate data in the main pressure pipeline of the target hydropower station into the nonlinear surge equation of the first surge chamber, the flow rate data in the main pressure pipeline of the target hydropower station is calculated as follows: This allows us to obtain a more accurate nonlinear surge equation for the first sub-pressure chamber.

[0123] Alternatively, if the target hydropower station is currently operating under increased load conditions, assuming the flow rate per unit is... A multi-machine hydropower station with one tunnel, when the same hydraulic unit The number of generating units operating normally has increased due to load dispatching by the power center. At that time, the flow rate in the main pressure pipeline was... Change to The flow rate data within the main pressure pipeline is... Substituting the flow data in the main pressure pipeline of the target hydropower station into the nonlinear surge equation of the first surge chamber This allows us to obtain a more accurate nonlinear surge equation for the second sub-pressure chamber.

[0124] In some embodiments, the surge equation of the pressure regulating chamber includes a first sub-equation and a second sub-equation.

[0125] Step A3 includes:

[0126] Step A31: Transform the pre-constructed nonlinear vibration equation to obtain the solution form equation corresponding to the nonlinear vibration equation;

[0127] The nonlinear vibration equation is specifically as follows:

[0128]

[0129] in, Indicates the independent variable. Indicates the independent variable The second derivative, Indicates frequency, Indicates time, Represents the nonlinear disturbance parameter. , Represents nonlinear forces;

[0130] The solution equation is specifically as follows:

[0131]

[0132] in, Indicates the independent variable. Represents the phase angle a periodic function, Represents the nonlinear disturbance parameter. , Indicates amplitude. Indicates the first A periodic function, The amplitude represents the total number of terms in a periodic function. and phase angle All times The function, specifically:

[0133]

[0134] in, Indicates amplitude. Indicates the phase angle. Indicates time, Indicates frequency, Represents the nonlinear disturbance parameter. , Indicates the first A periodic function, This represents the total number of terms in a periodic function. Indicates the equivalent damping ratio. Amplitude The function, Indicates the natural frequency of the equivalent period. Amplitude The function.

[0135] Step A32, determine the first-order approximate solution equation of the function, wherein the first-order approximate solution equation is specifically:

[0136]

[0137] in, Indicates the independent variable. Indicates amplitude, Indicates the phase angle. Indicates time, Represents the nonlinear disturbance parameter. , It represents pi (π).

[0138] Step A33, set the preset Substituting into the aforementioned nonlinear vibration equation, we obtain the nonlinear asymptotic equation;

[0139] The nonlinear asymptotic equation is specifically as follows:

[0140]

[0141] in, This represents the nonlinear disturbance parameter of the surge equation in the surge chamber under load shedding conditions. This represents the surge data of the surge chamber. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time.

[0142] Step A34, in response to the current operating condition being load shedding, and under the preset surge parameters... Substituting the nonlinear asymptotic equation into the first-order approximate solution equation yields the first substituted first-order approximate solution equation.

[0143] The first approximate solution equation after substitution is as follows:

[0144]

[0145] in, Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the nonlinear disturbance parameter of the surge equation in the surge chamber under load shedding conditions. This represents the nonlinear force in the surge equation of the surge chamber under load shedding conditions. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time.

[0146] Step A35: Solve the first approximate solution equation after substitution using ordinary differential equations to obtain the first ordinary differential equation, and use the first ordinary differential equation as the first sub-equation.

[0147] The first ordinary differential equation is as follows:

[0148]

[0149] in, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Indicates time, , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel at that time. Or,

[0150] Step A36, in response to the current operating condition being an increased load, and the preset surge parameters under the increased load being... and ,Will and Substituting the nonlinear asymptotic equation into the equation yields the nonlinear asymptotic equation after substitution.

[0151] The nonlinear asymptotic equation after substitution is specifically as follows:

[0152]

[0153] in, The nonlinear disturbance parameters represent the surge equation of the surge chamber under increased load conditions. This represents the simple harmonic solution of the pressure regulating chamber under increased load conditions. Indicates time, This indicates the number of operating units of the target hydropower station. This indicates the area of ​​the large well in the pressure regulating chamber. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time Indicates the surge phase, Indicates the amplitude of the surge wave. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel mentioned above.

[0154] Step A37: Substitute the nonlinear asymptotic equation into the first-order approximate solution equation to obtain the second first-order approximate solution equation after substitution.

[0155] The second, approximate first-order solution equation after substitution is as follows:

[0156]

[0157] in, The nonlinear disturbance parameters represent the surge equation of the surge chamber under increased load conditions. This represents the simple harmonic solution of the pressure regulating chamber under increased load conditions. Indicates time, This indicates the number of operating units of the target hydropower station. This indicates the area of ​​the large well in the pressure regulating chamber. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time Indicates the surge phase, Indicates the amplitude of the surge wave. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel mentioned above.

[0158] Step A38: Solve the second first-order approximate solution equation after substitution using ordinary differential equations to obtain the second ordinary differential equation, and use the second ordinary differential equation as the second sub-equation;

[0159] The second ordinary differential equation is as follows:

[0160]

[0161] in, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time Represents pi (π). , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. Indicates time, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the swell phase.

[0162] In the above scheme, for the nonlinear vibration equation, amplitude and phase angle All times The function, and the solution of the formal equation, when When the equation is a free vibration differential equation, its solution is a frequency of... The resonant solution; when At this time, nonlinear disturbances exist, and the solution to the equation contains higher harmonics. Since higher-order solutions are excessively complex, and their computational complexity increases geometrically with increasing order, and for practical engineering applications, the first-order approximate solution already possesses high computational accuracy, this application only uses a first-order approximate solution (i.e., the first-order approximate solution equation), which is as follows:

[0163] .

[0164] In some embodiments, the surge extreme value occurrence time includes the first surge extreme value occurrence time and the second surge extreme value occurrence time.

[0165] Step 103 includes:

[0166] Step 1031: In response to determining that the current operating condition is load shedding and that the derivative of the surge extreme value with respect to time is 0, the preset derivative value of 0 is substituted into the first sub-equation to obtain the first derivative value, which is then substituted into the equation.

[0167] The first sub-equation is as follows:

[0168]

[0169] in, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Indicates time, , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above;

[0170] Substituting the first derivative value into the equation is as follows:

[0171]

[0172] in, Indicates the extreme value of the surge. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Represents pi (π). , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: Head loss within the aforementioned water diversion tunnel;

[0173] Step 1032: Substitute the first derivative value into the equation for simplification to obtain the equation for the occurrence time of the surge extreme value in the first voltage regulating chamber, and solve the equation for the occurrence time of the surge extreme value in the first voltage regulating chamber to obtain the occurrence time of the first surge extreme value.

[0174] The specific time equation for the occurrence of the surge extreme value in the first pressure regulating chamber is as follows:

[0175]

[0176] in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. Indicates time, This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Represents pi (π). , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss within the water diversion tunnel at that time. Or,

[0177] Step 1033: In response to determining that the current operating condition is an increased load and that the derivative of the surge extreme value with respect to time is 0, the preset derivative value of 0 is substituted into the second sub-equation to obtain the second derivative value, which is then substituted into the equation.

[0178] The second sub-equation is as follows:

[0179]

[0180] in, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time Represents pi (π). , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. Indicates time, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase;

[0181] The second derivative value is substituted into the equation as follows:

[0182]

[0183] in, Indicates the extreme value of the surge. Indicates time, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load.

[0184] Step 1034: Substitute the second derivative value into the equation for simplification to obtain the equation for the occurrence time of the surge extreme value in the second surge chamber, and solve the equation for the occurrence time of the surge extreme value in the second surge chamber to obtain the occurrence time of the second surge extreme value.

[0185] The specific time equation for the occurrence of the surge extreme value in the second voltage regulating chamber is as follows:

[0186]

[0187] in, Indicates time, This indicates the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time.

[0188] In the above scheme, the current operating condition is determined to be load shedding. When the surge in the upstream pressure regulating chamber with the connecting pipe reaches an extreme value, there is... (That is, the derivative of the surge extreme with respect to time is 0), then the preset derivative value is 0. Substituting the preset derivative value of 0 into the first sub-equation yields the first derivative value. Substituting this first derivative value into the equation allows for a more accurate determination of the time when the first surge extreme occurs. Or,

[0189] The current operating condition is determined to be an increased load. When the surge in the upstream pressure regulating chamber with the connecting pipe reaches an extreme value, there is... (That is, the derivative of the surge extreme with respect to time is 0), then the preset derivative value is 0. Substituting the preset derivative value of 0 into the second sub-equation, we can obtain the second derivative value. Substituting this second derivative value into the equation can more accurately determine the time of occurrence of the second surge extreme.

[0190] In some embodiments, the initial parameters of the surge include the initial amplitude and the initial phase of the surge, and the extreme values ​​of the surge include a first extreme value and a second extreme value.

[0191] In step 104, determining the initial surge parameters corresponding to the current operating condition, and processing the surge extreme value using the surge extreme value occurrence time and the initial surge parameters through the surge determination algorithm in the surge chamber to obtain the surge extreme value, includes:

[0192] Step 1041: In response to the current operating condition being load shedding, substitute the initial amplitude of the surge corresponding to the load shedding, the initial phase of the surge corresponding to the load shedding, and the occurrence time of the first surge extreme value into the first sub-equation to obtain the first surge extreme value.

[0193] The initial amplitude of the surge corresponding to the load shedding Specifically:

[0194]

[0195] in, This indicates the number of operating units of the target hydropower station. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as Water head loss in the water diversion tunnel;

[0196] The initial phase of the surge corresponding to the load shedding Specifically:

[0197]

[0198] in, This indicates the number of operating units of the target hydropower station. The turbine's flow rate is indicated as Water head loss in the water diversion tunnel This represents the initial amplitude of the surge corresponding to the load shedding. So that An angle in the fourth quadrant. Represents pi;

[0199] The first sub-equation is as follows:

[0200]

[0201] in, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Indicates time, , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel at that time. Or,

[0202] Step 1042: In response to the current operating condition being an increased load, the initial amplitude of the surge corresponding to the increased load, the initial phase of the surge corresponding to the increased load, and the occurrence time of the second surge extreme value are substituted into the second sub-equation to obtain the second surge extreme value;

[0203] The initial amplitude of the surge corresponding to the increased load Specifically:

[0204]

[0205] in, This indicates the number of currently operating generating units at the target hydropower station. This indicates the original number of operating units of the target hydropower station. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This indicates the area of ​​the large well in the pressure regulating chamber. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe of the pressure regulating chamber;

[0206] The initial phase of the surge corresponding to the increased load Specifically:

[0207]

[0208] in, This indicates the number of currently operating generating units at the target hydropower station. This indicates the original number of operating units of the target hydropower station. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the initial amplitude of the surge corresponding to the increased load;

[0209] The second sub-equation is as follows:

[0210]

[0211] in, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as follows: The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as follows: The head loss in the water diversion tunnel at that time Represents pi (π). , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. Indicates time, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the swell phase.

[0212] In the above scheme, the current operating condition of the multi-unit hydropower station (i.e., the target hydropower station) is load shedding, and will Substituting these equations into the continuity equation of the surge tank, the momentum equation of the water diversion tunnel, and the head loss of the surge tank of the target hydropower station... In addition, the head loss of the water diversion tunnel of the target hydropower station is considered. The initial amplitude of the surge corresponding to the load shedding and the initial phase of the surge corresponding to the load shedding are obtained as follows:

[0213]

[0214] in, It is an angle in the fourth quadrant.

[0215] Substituting the initial surge amplitude corresponding to the load shedding, the initial surge phase corresponding to the load shedding, and the occurrence time of the first surge extreme value into the first sub-equation, a more accurate first surge extreme value can be obtained. Here, the occurrence time of the first surge extreme value represents the occurrence time of the highest surge extreme value, and the first surge extreme value represents the highest surge extreme value.

[0216] The current operating condition of the multi-unit hydropower station (i.e., the target hydropower station) is increased load, which will... Substituting these equations into the continuity equation of the surge tank, the momentum equation of the water diversion tunnel, and the head loss of the surge tank of the target hydropower station... In addition, the head loss of the water diversion tunnel of the target hydropower station is considered. The initial amplitude and initial phase of the surge corresponding to the increased load are obtained as follows:

[0217]

[0218] in, Angles in the first quadrant.

[0219] Substituting the initial surge amplitude corresponding to the load increase, the initial surge phase corresponding to the load increase, and the occurrence time of the second surge extreme value into the second sub-equation, a more accurate second surge extreme value can be obtained. Here, the occurrence time of the second surge extreme value represents the occurrence time of the minimum surge extreme value, and the second surge extreme value represents the minimum surge extreme value.

[0220] In some embodiments, this application relies on multiple single-tunnel multi-machine hydropower stations (i.e., target hydroelectric power stations) as actual engineering projects to verify the accuracy and effectiveness of the analytical formula for the extreme value of surge waves in the upstream connecting pipe surge chamber of a single-tunnel multi-machine hydropower station under load shedding / load increase conditions derived in this application.

[0221] Taking hydropower stations with one tunnel and two turbines, one tunnel and three turbines, and one tunnel and four turbines as examples, transient process calculation software was used as the reference for verification experiments to verify the analytical solution of surge waves in a surge-regulating chamber with an upstream connecting pipe under load shedding / load increase conditions. The specific operating conditions are as follows:

[0222] SOC-1 (Operating Condition 1): N generating units are operating normally, while suddenly shedding load;

[0223] SOC-2 (Operating Condition 2): N units are running under no-load conditions, and are simultaneously started up to increase to full load;

[0224] (1) As can be seen from Table 1, for a hydropower station with a single tunnel and two turbines, the analytical solution and software calculation results for the surge chamber of the connecting pipe are basically consistent. The calculation error of the surge extreme value is within 0.1%, and the error of the surge extreme value occurrence time is within 7.0 s, indicating high calculation accuracy. Table 1 is as follows:

[0225] Table 1. Comparison of analytical and software values ​​of surge waves in a typical operating condition of a single-tunnel, two-machine hydropower station with upstream connecting pipe surge chamber.

[0226]

[0227] (2) As can be seen from Table 2, for a hydropower station with a three-machine layout, the calculation results of the surge extreme value obtained by analytical formula and analytical model are basically consistent with the calculation value of software. The calculation error of the surge extreme value is within 0.1%, and the time error of the surge extreme value is about 8 s, which has high calculation accuracy.

[0228] Table 2 Comparison of analytical and software values ​​of surge waves in the upstream surge chamber with connecting pipe of the Yidong Sanji Hydropower Station under typical operating conditions.

[0229]

[0230] (3) As can be seen from Table 3, for a hydropower station with a one-tunnel-four-machine layout, the analytical calculation results of the surge extreme value are basically consistent with the software calculation values. The calculation error of the surge extreme value is within 0.5%, and the time error of the surge extreme value occurrence is within 6.0 s, which has high calculation accuracy.

[0231] Table 3. Comparison of analytical and software values ​​of surge waves in a surge chamber with an upstream connecting pipe under typical operating conditions of a one-tunnel, four-machine hydropower station.

[0232]

[0233] Based on the above calculation results, it can be seen that the method for solving the extreme values ​​of surge waves in a simple working condition of a multi-machine hydropower station with a surge chamber connected to a pipe in this application can accurately reflect the extreme value characteristics of surge waves in the surge chamber and has high calculation accuracy.

[0234] This application describes in detail the process of determining the height of the surge tank of a hydropower station using an example, as follows:

[0235] Upstream connection pipe surge chamber single-tunnel multi-unit hydropower station layout as follows: Figure 2A As shown, combined with Figure 2B Now proceed to step S1.

[0236] Step S1: Establish the unsteady flow equations for the tunnel-upstream surge chamber considering the inertia of the water flow in the connecting pipe and the head loss (i.e., the continuity equation of the surge chamber and the momentum equation of the water diversion tunnel), and derive the surge wave equation for the upstream surge chamber with connecting pipe (i.e., the surge wave equation of the surge chamber); give the definition of the nonlinear asymptotic method.

[0237] Step S2: Use the nonlinear asymptotic method to derive the analytical expressions for the surge extreme value (i.e., the first surge extreme value) and occurrence time (i.e., the occurrence time of the first surge extreme value) of the pressure regulating chamber of the upstream connecting pipe under the load shedding condition (simple condition).

[0238] Step S3: Use the nonlinear asymptotic method to derive the analytical expressions for the surge extreme value (i.e., the second surge extreme value) and occurrence time (i.e., the occurrence time of the second surge extreme value) of the pressure regulating chamber of the upstream connecting pipe under the increased load condition (simple condition).

[0239] Then, based on the surge peak value of the surge chamber of the upstream connecting pipe corresponding to step S2, and in conjunction with the safety freeboard, determine the top elevation of the upstream and downstream reservoirs. Here, the safety freeboard is, for example, 2-3 meters. Based on the surge peak value of the surge chamber of the upstream connecting pipe corresponding to step S3, determine the minimum safe distance that the upstream and downstream reservoirs should be from the lowest surge peak value. Here, the safe distance is, for example, 1 meter.

[0240] The process ends after step S3 is completed.

[0241] It should be noted that the method in this embodiment can be executed by a single device, such as a computer or server. The method can also be applied in a distributed scenario, where multiple devices cooperate to complete the task. In such a distributed scenario, one of these devices may execute only one or more steps of the method in this embodiment, and the multiple devices will interact with each other to complete the method described.

[0242] It should be noted that the above description describes some embodiments of this application. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in a different order than that shown in the above embodiments and still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require a specific or sequential order to achieve the desired result. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0243] Based on the same inventive concept, and corresponding to any of the above embodiments, this application also provides a device for determining the height of a surge tank in a hydropower station.

[0244] refer to Figure 3 The device for determining the height of the surge tank of the hydropower station includes:

[0245] The head loss determination module 301 is configured to determine the head loss of the surge tank of the target hydropower station, and to determine;

[0246] The fusion processing module 302 is configured to acquire the flow data in the main pressure pipeline corresponding to the current operating condition of the target hydropower station, and process the head loss of the surge chamber, the head loss of the water diversion tunnel and the flow data in the main pressure pipeline through the surge chamber fusion algorithm to obtain the surge chamber determination algorithm.

[0247] The time determination module 303 is configured to process the surge determination algorithm based on the derivative of the surge extreme value with respect to time as a preset derivative value to obtain the time of occurrence of the surge extreme value;

[0248] The height determination module 304 is configured to determine the initial parameters of the surge corresponding to the current operating condition, process the surge extreme value and the initial parameters of the surge using the surge extreme value and the surge determination algorithm to obtain the surge extreme value, and determine the height of the surge chamber based on the surge extreme value.

[0249] For ease of description, the above devices are described in terms of function, divided into various modules. Of course, in implementing this application, the functions of each module can be implemented in one or more software and / or hardware.

[0250] The apparatus described above is used to implement the corresponding method for determining the height of the surge tank of a hydropower station in any of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0251] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, this application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for determining the height of the surge tank of a hydropower station as described in any of the above embodiments.

[0252] Figure 4 This illustration shows a more specific hardware structure diagram of an electronic device provided in this embodiment. The device may include: a processor 401, a memory 402, an input / output interface 403, a communication interface 404, and a bus 405. The processor 401, memory 402, input / output interface 403, and communication interface 404 are interconnected internally via the bus 405.

[0253] The processor 401 can be implemented using a general-purpose CPU (Central Processing Unit), microprocessor, application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0254] The memory 402 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, etc. The memory 402 can store the operating system and other applications. When the technical solutions provided in the embodiments of this specification are implemented by software or firmware, the relevant program code is stored in the memory 402 and is called and executed by the processor 401.

[0255] Input / output interface 403 is used to connect input / output modules to realize information input and output. Input / output modules can be configured as components in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. Input devices may include keyboards, mice, touch screens, microphones, various sensors, etc., and output devices may include displays, speakers, vibrators, indicator lights, etc.

[0256] Communication interface 404 is used to connect a communication module (not shown in the figure) to enable communication between this device and other devices. The communication module can communicate via wired means (such as USB, Ethernet cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0257] Bus 405 includes a pathway for transmitting information between various components of the device, such as processor 401, memory 402, input / output interface 403, and communication interface 404.

[0258] It should be noted that although the above-described device only shows the processor 401, memory 402, input / output interface 403, communication interface 404, and bus 405, in specific implementations, the device may also include other components necessary for normal operation. Furthermore, those skilled in the art will understand that the above-described device may only include the components necessary for implementing the embodiments of this specification, and not necessarily all the components shown in the figures.

[0259] The electronic equipment described above is used to implement the corresponding method for determining the height of the surge tank of a hydropower station in any of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0260] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, this application also provides a non-transitory computer-readable storage medium storing computer instructions for causing the computer to execute the method for determining the height of the surge tank of a hydropower station as described in any of the above embodiments.

[0261] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. 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, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transfer medium that can be used to store information accessible by a computing device.

[0262] The computer instructions stored in the storage medium of the above embodiments are used to cause the computer to execute the method for determining the height of the surge tank of a hydropower station as described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0263] Those skilled in the art should understand that the discussion of any of the above embodiments is merely exemplary and is not intended to imply that the scope of this application (including the claims) is limited to these examples; within the framework of this application, the technical features of the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations of different aspects of the embodiments of this application as described above, which are not provided in the details for the sake of brevity.

[0264] Additionally, to simplify the description and discussion, and to avoid obscuring the embodiments of this application, the well-known power / ground connections to integrated circuit (IC) chips and other components may or may not be shown in the provided drawings. Furthermore, the apparatus may be shown in block diagram form to avoid obscuring the embodiments of this application, and this also takes into account the fact that the details of the implementation of these block diagram apparatuses are highly dependent on the platform on which the embodiments of this application will be implemented (i.e., these details should be fully understood by those skilled in the art). While specific details (e.g., circuits) have been set forth to describe exemplary embodiments of this application, it will be apparent to those skilled in the art that the embodiments of this application can be implemented without these specific details or with variations thereof. Therefore, these descriptions should be considered illustrative rather than restrictive.

[0265] Although this application has been described in conjunction with specific embodiments thereof, many substitutions, modifications, and variations of these embodiments will be apparent to those skilled in the art from the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may be used with the embodiments discussed.

[0266] The embodiments of this application are intended to cover all such substitutions, modifications, and variations that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the embodiments of this application should be included within the protection scope of this application.

Claims

1. A method for determining the height of a surge chamber in a hydropower station, characterized in that, include: Determine the head loss of the surge tank of the target hydropower station, and determine the head loss of the water diversion tunnel of the target hydropower station; Obtain the flow rate data in the main pressure pipeline corresponding to the current operating condition of the target hydropower station. Based on the head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline, process them through the surge chamber surge fusion algorithm to obtain the surge chamber surge determination algorithm. Based on the derivative of the surge extreme value with respect to time as a preset derivative value, the surge extreme value is processed according to the surge determination algorithm of the pressure regulating chamber to obtain the surge extreme value occurrence time; Determine the initial parameters of the surge corresponding to the current operating condition, process the surge extreme value and the initial parameters of the surge using the surge determination algorithm of the surge chamber to obtain the surge extreme value, and determine the height of the surge chamber based on the surge extreme value; The determination of the head loss of the surge tank of the target hydropower station and the determination of the head loss of the water diversion tunnel of the target hydropower station include: Acquire surge data of the pressure regulating chamber The turbine flow rate of the target hydropower station is... The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is Head loss in the water diversion tunnel The water turbine's flow rate is The water flow velocity in the water diversion tunnel at that time The large well area of ​​the pressure regulating chamber The area inside the water diversion tunnel and the flow velocity inside the water diversion tunnel ; Based on the surge data of the pressure regulating chamber The water turbine's flow rate is The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is Head loss in the water diversion tunnel The water turbine's flow rate is The water flow velocity in the water diversion tunnel at that time The large well area of ​​the pressure regulating chamber and the area inside the water diversion tunnel The head loss of the surge tank of the target hydropower station is determined by the following formula. : in, ; Based on the preset head loss coefficient of the water diversion tunnel and the flow velocity inside the water diversion tunnel The head loss of the water diversion tunnel of the target hydropower station is determined by the following formula. : ； The head loss of the surge chamber, the head loss of the water diversion tunnel, and the flow rate data in the main pressure pipeline are processed using a surge chamber surge algorithm to obtain a surge chamber surge determination algorithm, including: Based on the head loss of the surge chamber and the head loss of the water diversion tunnel, the first nonlinear surge equation of the surge chamber is obtained by combining the pre-constructed continuity equation of the surge chamber with the pre-constructed momentum equation of the water diversion tunnel. Based on the flow rate data in the main pressure pipeline, the nonlinear surge equation of the second pressure regulating chamber is obtained through the nonlinear surge equation of the first pressure regulating chamber; By using the pre-set surge parameters under the current operating condition, and combining them with the pre-constructed nonlinear vibration equation and the surge equation of the second surge chamber, the surge equation of the surge chamber is obtained.

2. The method according to claim 1, characterized in that, The head loss based on the surge chamber and the head loss of the water diversion tunnel are combined using a pre-constructed continuity equation for the surge chamber and a pre-constructed momentum equation for the water diversion tunnel to obtain the first nonlinear surge equation for the surge chamber, which includes: The head loss in the pressure regulating chamber and the head loss of the water diversion tunnel Substituting the continuity equation of the surge tank and simultaneously solving the momentum equation of the water diversion tunnel, the flow velocity inside the water diversion tunnel is... After elimination processing, the nonlinear surge equation of the first pressure regulating chamber is obtained; The continuity equation for the pressure regulating chamber is as follows: in, This indicates the length of the water diversion tunnel. This indicates the area of ​​the water diversion tunnel. Indicates the flow rate inside the water diversion tunnel. This indicates the head loss of the water diversion tunnel. This represents the surge data of the surge chamber. This indicates the length of the connecting pipe to the pressure regulating chamber. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the flow rate within the connecting pipe of the pressure regulating chamber. Represents gravitational acceleration. Indicates time, This indicates the head loss caused by the inflow and outflow from the pressure regulating chamber; The momentum equation for the water diversion tunnel is as follows: in, This indicates the area of ​​the large well in the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, Indicates the flow rate inside the water diversion tunnel. This represents the flow rate data within the main pressure pipeline of the target hydropower station; The nonlinear surge equation for the first pressure regulating chamber is as follows: in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, The turbine's flow rate is indicated as Water head loss in the water diversion tunnel This represents the flow rate data within the main pressure pipeline of the target hydropower station. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time.

3. The method according to claim 1, characterized in that, The second nonlinear surge equation for the pressure regulating chamber includes the first nonlinear surge equation for the pressure regulating chamber and the second nonlinear surge equation for the pressure regulating chamber. The process of obtaining the second nonlinear surge equation for the pressure regulating chamber based on the flow rate data in the main pressure pipeline and the nonlinear surge equation of the first pressure regulating chamber includes: In response to the target hydropower station's current operating condition being load shedding, and the corresponding flow rate data in the main pressure pipeline being 0, the flow rate data in the main pressure pipeline of the target hydropower station is substituted into the nonlinear surge equation of the first surge chamber as 0. The nonlinear surge equation for the first sub-pressure regulating chamber is obtained. The nonlinear surge equation for the first pressure regulating chamber is as follows: in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This represents the flow rate data within the main pressure pipeline of the target hydropower station. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as Head loss within the aforementioned water diversion tunnel; The nonlinear surge equation for the first sub-pressure regulating chamber is as follows: in, This represents the surge data of the surge chamber. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity within the water diversion tunnel at that time; or, In response to the target hydropower station's current operating condition being an increased load, and the corresponding flow rate data in the main pressure pipeline is as follows: The flow rate data in the main pressure pipeline is... Substituting the flow data in the main pressure pipeline of the target hydropower station into the nonlinear surge equation of the first surge chamber... The nonlinear surge equation for the second sub-pressure regulating chamber is obtained. The nonlinear surge equation for the first pressure regulating chamber is as follows: in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the surge data of the surge chamber. Indicates time, The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This represents the flow rate data within the main pressure pipeline of the target hydropower station. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as Head loss within the aforementioned water diversion tunnel; The nonlinear surge equation for the second sub-pressure regulating chamber is as follows: in, This represents the simple harmonic solution of the pressure regulating chamber under increased load conditions. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time.

4. The method according to claim 3, characterized in that, The surge equation for the pressure regulating chamber includes a first sub-equation and a second sub-equation; The process of obtaining the surge equation of the surge chamber by using the preset surge parameters under the current operating condition, through a pre-constructed nonlinear vibration equation, and combining it with the surge equation of the second surge chamber, includes: The pre-constructed nonlinear vibration equation is transformed to obtain the solution form equation corresponding to the nonlinear vibration equation; The nonlinear vibration equation is specifically as follows: in, Indicates the independent variable. Indicates the independent variable The second derivative, Indicates frequency, Indicates time, Represents the nonlinear disturbance parameter. , Represents nonlinear forces; The solution equation is specifically as follows: in, Indicates the independent variable. Represents the phase angle a periodic function, Represents the nonlinear disturbance parameter. , Indicates amplitude. Indicates the first A periodic function, This represents the total number of terms and amplitude of a periodic function. and phase angle All times The function, specifically: in, Indicates amplitude. Indicates the phase angle. Indicates time, Indicates frequency, Represents the nonlinear disturbance parameter. , Indicates the first A periodic function, This represents the total number of terms in a periodic function. Indicates the equivalent damping ratio. Amplitude The function, Indicates the natural frequency of the equivalent period. Amplitude The function; Determine the first-order approximate solution equation of the function, wherein the first-order approximate solution equation is specifically: in, Indicates the independent variable. Indicates amplitude. Indicates the phase angle. Indicates time, Represents the nonlinear disturbance parameter. , Represents pi; Preset Substituting into the aforementioned nonlinear vibration equation, we obtain the nonlinear asymptotic equation; The nonlinear asymptotic equation is specifically as follows: in, This represents the nonlinear disturbance parameter of the surge equation in the surge chamber under load shedding conditions. This represents the surge data of the surge chamber. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time; In response to the current operating condition being load shedding, and under the preset surge parameters... Substituting the nonlinear asymptotic equation into the first-order approximate solution equation yields the first substituted first-order approximate solution equation. The first approximate solution equation after substitution is as follows: in, Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the nonlinear disturbance parameter of the surge equation in the surge chamber under load shedding conditions. This represents the nonlinear force in the surge equation of the surge chamber under load shedding conditions. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time; Solve the first approximate solution equation by substituting it into an ordinary differential equation to obtain the first ordinary differential equation, and use the first ordinary differential equation as the first sub-equation. The first ordinary differential equation is as follows: in, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Indicates time, , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel at that time; or, In response to the current operating condition being an increased load, and under the increased load, the preset surge parameters are as follows: and ,Will and Substituting the nonlinear asymptotic equation into the equation yields the nonlinear asymptotic equation after substitution. The nonlinear asymptotic equation after substitution is specifically as follows: in, The nonlinear disturbance parameters represent the surge equation of the surge chamber under increased load conditions. This represents the simple harmonic solution of the pressure regulating chamber under increased load conditions. Indicates time, This indicates the number of operating units of the target hydropower station. This indicates the area of ​​the large well in the pressure regulating chamber. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time Indicates the surge phase, Indicates the amplitude of the surge wave. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as Head loss within the aforementioned water diversion tunnel; Substituting the nonlinear asymptotic equation into the first-order approximate solution equation yields the second first-order approximate solution equation after substitution. The second, approximate first-order solution equation after substitution is as follows: in, The nonlinear disturbance parameters represent the surge equation of the surge chamber under increased load conditions. This represents the simple harmonic solution of the pressure regulating chamber under increased load conditions. Indicates time, This indicates the number of operating units of the target hydropower station. This indicates the area of ​​the large well in the pressure regulating chamber. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time Indicates the surge phase, Indicates the amplitude of the surge wave. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as Head loss within the aforementioned water diversion tunnel; Solve the second first-order approximate solution equation by substituting it into an ordinary differential equation to obtain the second ordinary differential equation, and use the second ordinary differential equation as the second sub-equation; The second ordinary differential equation is as follows: in, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time Represents pi (π). , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. Indicates time, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase.

5. The method according to claim 4, characterized in that, The time of occurrence of the surge extreme value includes the time of occurrence of the first surge extreme value and the time of occurrence of the second surge extreme value; The process of determining the surge extreme value based on the surge extreme value with respect to time, using a preset derivative value, is performed according to the surge determination algorithm of the surge chamber to obtain the surge extreme value occurrence time, including: In response to determining that the current operating condition is load shedding and that the derivative of the surge extreme value with respect to time is 0, the preset derivative value of 0 is substituted into the first sub-equation to obtain the first derivative value, which is then substituted into the equation. The first sub-equation is as follows: in, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Indicates time, , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above; Substituting the first derivative value into the equation is as follows: in, Indicates the extreme value of the surge. Indicates time, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Represents pi (π). , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as Head loss within the aforementioned water diversion tunnel; Substituting the first derivative value into the equation for simplification, we obtain the equation for the occurrence time of the surge extreme value in the first voltage regulating chamber. Solving the equation for the occurrence time of the surge extreme value in the first voltage regulating chamber, we obtain the occurrence time of the first surge extreme value. The specific time equation for the occurrence of the surge extreme value in the first pressure regulating chamber is as follows: in, , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. Indicates time, This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Represents pi (π). , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss within the water diversion tunnel at that time; or, In response to determining that the current operating condition is an increased load and that the derivative of the surge extreme value with respect to time is 0, the preset derivative value of 0 is substituted into the second sub-equation to obtain the second derivative value, which is then substituted into the equation. The second sub-equation is as follows: in, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time Represents pi (π). , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. Indicates time, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase; The second derivative value is substituted into the equation as follows: in, Indicates the extreme value of the surge. Indicates time, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load; Substituting the second derivative value into the equation for simplification, the equation for the occurrence time of the surge extreme value in the second surge chamber is obtained. The equation for the occurrence time of the surge extreme value in the second surge chamber is then solved to obtain the occurrence time of the second surge extreme value. The specific time equation for the occurrence of the surge extreme value in the second voltage regulating chamber is as follows: in, Indicates time, This indicates the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time.

6. The method according to claim 5, characterized in that, The initial parameters of the surge include the initial amplitude and the initial phase of the surge, and the extreme values ​​of the surge include the first extreme value and the second extreme value of the surge. The process of determining the initial surge parameters corresponding to the current operating condition, and then processing the surge extreme values ​​using the surge extreme value occurrence time and the initial surge parameters through the surge determination algorithm in the surge chamber to obtain the surge extreme value, includes: In response to the current operating condition being load shedding, the initial amplitude of the surge corresponding to the load shedding, the initial phase of the surge corresponding to the load shedding, and the occurrence time of the first surge extreme value are substituted into the first sub-equation to obtain the first surge extreme value; The initial amplitude of the surge corresponding to the load shedding Specifically: in, This indicates the number of operating units of the target hydropower station. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as Water head loss in the water diversion tunnel; The initial phase of the surge corresponding to the load shedding Specifically: in, This indicates the number of operating units of the target hydropower station. The turbine's flow rate is indicated as Water head loss in the water diversion tunnel This represents the initial amplitude of the surge corresponding to the load shedding. So that An angle in the fourth quadrant. Represents pi; The first sub-equation is as follows: in, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase, Represents pi (π). , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the load shedding. This indicates the initial phase of the surge corresponding to the load shedding. Indicates time, , This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time The turbine's flow rate is indicated as The head loss of the water diversion tunnel at that time; or, In response to the current operating condition being an increased load, the initial amplitude of the surge corresponding to the increased load, the initial phase of the surge corresponding to the increased load, and the occurrence time of the second surge extreme value are substituted into the second sub-equation to obtain the second surge extreme value; The initial amplitude of the surge corresponding to the increased load Specifically: in, This indicates the number of currently operating generating units at the target hydropower station. This indicates the original number of operating units of the target hydropower station. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the area of ​​the water diversion tunnel. The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This indicates the area of ​​the large well in the pressure regulating chamber. , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe of the pressure regulating chamber; The initial phase of the surge corresponding to the increased load Specifically: in, This indicates the number of currently operating generating units at the target hydropower station. This indicates the original number of operating units of the target hydropower station. The turbine's flow rate is indicated as The head loss of the water diversion tunnel mentioned above, This indicates the initial amplitude of the surge corresponding to the increased load; The second sub-equation is as follows: in, , This indicates the flow rate referenced by the water turbine. This indicates the area of ​​the large well in the pressure regulating chamber. , , This indicates the natural frequency of the surge equivalent period in the pressure regulating chamber. Represents gravitational acceleration. This indicates the area of ​​the large well in the pressure regulating chamber. This indicates the area of ​​the water diversion tunnel. This indicates the area of ​​the connecting pipe of the pressure regulating chamber. This indicates the length of the water diversion tunnel. This indicates the length of the connecting pipe to the pressure regulating chamber. , The water turbine's flow rate is indicated as The head loss generated by the pressure regulating chamber at that time The water turbine's flow rate is indicated as The head loss in the water diversion tunnel at that time Represents pi (π). , The turbine's flow rate is indicated as The water flow velocity inside the water diversion tunnel at that time This represents the initial amplitude of the surge corresponding to the increased load. This indicates the initial phase of the surge corresponding to the increased load. Indicates time, Indicates the extreme value of the surge. Indicates the amplitude of the surge wave. Indicates the surge phase.

7. A device for determining the height of a surge tank in a hydropower station, characterized in that, The method for determining the height of a surge chamber in a hydropower station according to any one of claims 1 to 6 includes: The head loss determination module is configured to determine the head loss of the surge tank of the target hydropower station and the head loss of the water diversion tunnel of the target hydropower station. The fusion processing module is configured to acquire the flow data in the main pressure pipeline corresponding to the current operating condition of the target hydropower station, and process the head loss of the surge chamber, the head loss of the water diversion tunnel and the flow data in the main pressure pipeline through the surge chamber fusion algorithm to obtain the surge chamber determination algorithm. The time determination module is configured to process the surge determination algorithm based on the derivative of the surge extreme value with respect to time as a preset derivative value to obtain the time of the surge extreme value occurrence. The height determination module is configured to determine the initial parameters of the surge corresponding to the current operating condition, process the surge extreme value and the initial parameters of the surge using the surge extreme value and the surge determination algorithm to obtain the surge extreme value, and determine the height of the surge chamber based on the surge extreme value.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method as described in any one of claims 1 to 6.

Images