Output distribution calculation method, system and equipment for reducing startup and shutdown times of mixed-flow water turbine set based on binary arithmetic progression algorithm and medium
By optimizing the output distribution of the mixed-flow turbine unit based on the binary arithmetic series algorithm, the problem of equipment aging caused by frequent start-ups and shutdowns was solved, and the safe and stable operation and lifespan extension of the equipment were achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-10-10
- Publication Date
- 2026-03-10
AI Technical Summary
Mixed-flow turbine units frequently start and stop during peak-shaving operations, leading to equipment aging and increased failure rates. Traditional calculation methods are inefficient and cannot effectively avoid vibration zones, resulting in equipment wear and failure.
By adopting a binary arithmetic series algorithm, a calculation formula is established by defining the stable zone and the restricted operating zone of a single unit, optimizing the output distribution scheme, adjusting the load curve in real time, avoiding the unit from entering the vibration zone, and reducing the number of start-ups and shutdowns.
It effectively reduces the number of start-ups and shutdowns of mixed-flow turbine units, reduces equipment wear and failure rate, improves equipment lifespan and operational safety, and enhances calculation efficiency and accuracy.
Smart Images
Figure CN121643091A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of hydropower operation optimization, and in particular to a power distribution calculation method, system, device and medium for reducing the start-stop times of a Francis turbine unit based on a binary arithmetic progression algorithm. BACKGROUND
[0002] With the construction of new power systems, the functions of large hydropower stations are undergoing an important transformation. Traditional large hydropower stations mainly bear basic output and provide stable power output to meet basic electricity demand. However, with the large-scale integration of new energy such as wind and solar energy, the power system has significantly increased requirements for flexibility and regulation capacity. New power systems are dominated by new energy, which is intermittent and volatile, so they need power sources with fast response and regulation capacity to balance power supply and demand and ensure stable system operation.
[0003] In this context, the functions of large hydropower stations are gradually shifting from bearing basic output to bearing peak shaving tasks. Peak shaving tasks require power sources to quickly adjust power generation according to the needs of the power grid to respond to changes in output. Large hydropower units using Francis turbines need to constantly change the number of operating units to avoid running in the vibration zone in response to changes in output, resulting in an increase in start-stop times. Frequent start-stop processes can cause component aging and equipment failure. SUMMARY
[0004] In view of the above existing problems, the present application is proposed.
[0005] To solve the above technical problems, the present application provides the following technical solutions: In a first aspect, a power distribution calculation method for reducing the start-stop times of a Francis turbine unit based on a binary arithmetic progression algorithm, comprising: Defining the single-machine stable zone and the limited operating zone of the hydropower unit through actual experiments on the power station; Inputting the water head value, and performing a first calculation to obtain the unit output range under the input water head value according to the single-machine stable zone and the limited operating zone of the hydropower unit; Determining the total output, and outputting a unit output distribution scheme in combination with the number of units according to the unit output range; Performing a second calculation according to the unit output distribution scheme, and outputting the latest output distribution scheme; Optimizing the latest output distribution scheme, performing a third calculation, and outputting the optimized distribution scheme; Applying the optimized distribution scheme to real-time calculation of the combined vibration zone, checking whether the issued load curve falls into the combined vibration zone, and making corresponding adjustments.
[0006] As a preferred scheme of the output distribution calculation method for reducing the start-stop times of a mixed-flow hydraulic turbine unit based on a binary arithmetic progression algorithm, the input water head value is calculated according to the single-unit stable zone and the limited operation zone of the hydraulic turbine unit to obtain the unit output range under the input water head value, including: The input water head value is established as a first calculation formula. The output range is calculated according to the first calculation formula and in combination with the single-unit stable zone and the limited operation zone of the hydraulic turbine unit. The unit output range under the current water head value is output.
[0007] The preferred technical scheme has the beneficial effects that the first calculation formula is established, the single-unit stable zone and the limited operation zone of the hydraulic turbine unit are combined, the unit output range under the corresponding water head value can be calculated according to any input water head value, the limitation of only relying on fixed water head data is avoided, the system can correctly arrange the unit to operate in the stable zone or the short-time limited operation zone after accurately calculating the unit output range, the unit is prevented from entering the prohibited operation zone due to improper output distribution, and thus the equipment wear and tear and the fault risk are reduced.
[0008] As a preferred scheme of the output distribution calculation method for reducing the start-stop times of a mixed-flow hydraulic turbine unit based on a binary arithmetic progression algorithm, the second calculation is performed according to the unit output distribution scheme to output the latest output distribution scheme, including: The second calculation formula is constructed according to the unit output distribution scheme to calculate the adjustable output value of the current distribution scheme. Wherein, F(n, NO) represents the minimum output value of various combinations of the units operating in the stable zone and the limited operation zone when the number of operating units NO is given, NO is the total number of operating units, represents the adjustable output value of various combinations of the units operating in the stable zone and the limited operation zone when the number of operating units NO is given, the minimum output value of various combinations of the units operating in the stable zone and the limited operation zone, is the lower limit value of the limited operation zone, n is the number of units operating in the limited operation zone, is the upper limit value of the limited operation zone, is the upper limit value of the single-unit stable zone, is the lower limit value of the single-unit stable zone. The second judgment formula is constructed according to the adjustable output value of the current distribution scheme to judge whether the unit output distribution scheme is feasible. The adjustable capacity is proportionally distributed according to the currently feasible unit output distribution scheme, and the latest unit output distribution scheme is output.
[0009] The beneficial effects of the preferred technical solution are that the application converts the complex unit combination problem into a calculable mathematical model by constructing the second calculation formula, and the traditional method needs to enumerate all possible unit combinations to find a feasible output distribution scheme, which is time-consuming and has a large amount of calculation; the patent directly calculates the total output lower limit and total adjustable range corresponding to each combination by using the second calculation formula, and then combines the second judgment formula to judge the feasibility, and directly outputs the feasible unit output distribution scheme.
[0010] As a preferred scheme of the output distribution calculation method for reducing the start-stop times of the mixed-flow hydraulic turbine unit based on the binary arithmetic progression algorithm, the latest output distribution scheme is optimized, the third calculation is performed, and the optimized distribution scheme is output, including: The latest output distribution scheme is optimized again, the third calculation is performed, and the optimal output point under the current distribution scheme limit operating area is obtained: Wherein, PZ_XZQ_E represents the optimal output point under the current distribution scheme limit operating area, H is the current water head value, Hmax and Hmin are the highest water head value and the lowest water head value respectively; According to the optimal output point under the current distribution scheme limit operating area, a quadratic optimization distribution formula is constructed; According to the quadratic optimization distribution formula, the optimized distribution scheme is output.
[0011] The beneficial effects of the preferred technical solution are that the vibration levels of different output points within the limit operating area are different, the application outputs the optimal output point of the limit operating area by constructing the third calculation formula, finds a theoretically minimum vibration point, and then combines the third judgment formula to approximate the actual output direction of the unit to the optimal point, so that the system actively and purposefully adjusts the actual output direction of the unit to the known operating point with the best vibration performance under the current water head.
[0012] In the second aspect, the application provides an output distribution calculation system for reducing the start-stop times of the mixed-flow hydraulic turbine unit based on the binary arithmetic progression algorithm, which comprises an initialization module, a first calculation module, a scheme output module, a second calculation module, a third calculation module and an adjustment module. The initialization module is used to define the single-machine stable area and the limit operating area of the power station hydraulic turbine unit through actual experiments of the power station. The first calculation module is used to input the water head value, and perform the first calculation to obtain the unit output range under the input water head value according to the single-machine stable area and the limit operating area of the hydraulic turbine unit. The scheme output module is used to determine the total output, output the unit output distribution scheme according to the unit output range and the number of units. The second calculation module is used to perform a second calculation based on the unit output allocation scheme and output the latest output allocation scheme; The third calculation module is used to optimize the latest power allocation scheme, perform a third calculation, and output the optimized allocation scheme. The adjustment module is used to apply the optimized allocation scheme, calculate the combined vibration zone in real time, check whether the issued load curve falls into the combined vibration zone, and make corresponding adjustments.
[0013] Thirdly, a computer device includes a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of a method for calculating the output distribution of a mixed-flow turbine unit based on a binary arithmetic series algorithm to reduce the number of start-ups and shutdowns.
[0014] Fourthly, a computer-readable storage medium storing a computer program thereon, characterized in that, when the computer program is executed by a processor, it implements the steps of a method for calculating the output distribution of a mixed-flow turbine unit based on a binary arithmetic series algorithm to reduce the number of start-ups and shutdowns.
[0015] Compared with the prior art, the beneficial effects of the present invention are as follows: During peak-shaving operations, the most direct method to avoid generating units falling into the vibration zone is frequent start-up and shutdown to match load changes. This leads to accelerated aging of unit components, increased failure rate, and soaring maintenance costs. This invention calculates and actively utilizes the restricted operating range for short periods, allowing for adjustments to the output of online units to respond to load fluctuations, rather than necessitating start-up and shutdown. This fundamentally reduces unnecessary start-up and shutdown operations, lowers mechanical wear and fatigue damage, and extends overhaul cycles and equipment lifespan. Attached Figure Description
[0016] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the following description of the embodiments will be briefly introduced. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0017] Figure 1 This is a schematic diagram of the overall process of a method for calculating the output distribution of a mixed-flow turbine unit based on a binary arithmetic series algorithm to reduce the number of start-ups and shutdowns, as described in one embodiment of the present invention. Detailed Implementation
[0018] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the protection scope of the present invention.
[0019] Example 1, referring to Figure 1 As an embodiment of the present invention, a method for calculating output distribution to reduce the number of start-ups and shutdowns of a mixed-flow turbine unit based on a binary arithmetic series algorithm is provided, comprising: It should be noted that with the construction of new power systems, the function of large hydropower stations is gradually shifting from providing basic power output to undertaking peak-shaving tasks. Peak-shaving tasks require power sources to quickly adjust their power generation according to the grid's demand to cope with changes in output. After large hydropower units using mixed-flow turbines are shifted to peak-shaving tasks, a certain vibration zone exists. To cope with output changes and avoid operating in the vibration zone, the number of operating units needs to be constantly changed, resulting in an increase in the number of start-ups and shutdowns. Frequent start-ups and shutdowns can easily cause component aging and equipment failure.
[0020] Therefore, to address the aforementioned issues such as frequent unit start-ups and shutdowns leading to significant equipment wear, the following steps (S1-S6) are implemented: defining the single-unit stable zone and restricted operating zone of the power station turbine unit; performing a first calculation to obtain the unit output range under the input head value; and outputting a unit output allocation scheme, thus solving the problem of low allocation scheme generation efficiency. A second calculation is performed to determine the feasibility of the unit output allocation scheme. A third calculation is performed to output an optimized allocation scheme, ensuring that the unit's output within the restricted operating zone is at the point of minimum vibration, avoiding frequent start-ups and shutdowns, and guaranteeing equipment operating safety and lifespan. Finally, the issued load curve is checked to ensure it falls within the combined vibration zone, and appropriate adjustments are made.
[0021] S1: Define the single-unit stable zone and restricted operating zone of the power station turbine unit through actual power station experiments; S2: Input the head value, and perform the first calculation based on the single unit's stable zone and restricted operating zone to obtain the unit's output range under the input head value; S3: Determine the total output, and output the unit output allocation scheme based on the unit output range and the number of units; S4: Based on the unit output allocation scheme, perform the second calculation and output the latest output allocation scheme; S5: Optimize the latest power allocation scheme, perform the third calculation, and output the optimized allocation scheme; S6: Apply the optimized allocation scheme, calculate the combined vibration zone in real time, check whether the issued load curve falls into the combined vibration zone, and make corresponding adjustments.
[0022] Example 2, refer to Figure 1 As an embodiment of the present invention, based on the above embodiment, a method for calculating the output distribution of a mixed-flow turbine unit by reducing the number of start-ups and shutdowns is provided based on a binary arithmetic series algorithm.
[0023] In step S1, through actual power plant experiments, the single-unit stable zone and restricted operating zone of the power plant turbine are defined, including the following steps A1-A3: A1: Assign values to the maximum and minimum values of the single unit's stable zone under different head conditions obtained through actual tests at the power station, and output a table of operating parameters for the turbine unit under different head values.
[0024] Table 1 shows the operating parameters of the turbine unit under different head values: Table 1. Operating parameters of the turbine unit at different head values.
[0025] For example, Table 2 shows a specific example of a power plant, where the range from 0 to minimum output is defined as the single-unit vibration zone: Table 2. Operating parameters of a power station turbine unit under different head values.
[0026] A2: Based on the operating parameter table for different head values of the turbine unit, the stable zone of a single unit is obtained.
[0027] Based on the operating parameter table for different head values of the turbine unit, the stable zone of a single unit is determined to be [Pmin, Pmix].
[0028] A3: The maximum and minimum values of the restricted operating area under full head are obtained through experiments at the power station, and the restricted operating area is obtained by assigning these values.
[0029] Assign values to the restricted operating range [PLMin, PLMax] obtained from the power station's tests under full head conditions.
[0030] It should be noted that when the mixed-flow turbine operates in this region, some data may slightly exceed the standard, allowing for short-term dwell times. According to relevant standards, the annual cumulative operating hours in this range must not exceed a certain standard, and this is called the restricted operating zone. This range is the minimum intersection of all allowable short-term dwell time ranges obtained from tests under all head conditions, where [PLMin,PLMax]∈(0,Pmin).
[0031] In step S2, an input water head value is entered. According to the single-unit stable region and the restricted operation region of the water turbine unit, a first calculation is performed to obtain the unit output range at the input water head value, including the following steps B1 - B3: B1: Enter the water head value and establish a first calculation formula.
[0032] B2: According to the first calculation formula, in combination with the single-unit stable region and the restricted operation region of the water turbine unit, perform an output range calculation.
[0033] B3: Output the unit output range at the current water head value.
[0034] In the embodiment of the present application, the interpolation calculation method is adopted for the first calculation formula in B2. The specific steps are as follows: Assume that the input water head value is H. Query Table 1. If H = Hn (Hn corresponds to H1 - H11 in Table 1), then directly assign values: Among them, is the minimum output at the current water head, is the maximum output at the current water head, corresponds to Pmin1 - Pmin11 in Table 1, corresponds to Pmax1 - Pmax11 in Table 1.
[0035] If Hn < H <= Hn + 1, then calculate according to the interpolation calculation formula: In an optional embodiment, the first calculation formula in step B2 can also adopt a method of fitting the output characteristic surface based on a neural network. Collect a large amount of historical operation data of the hydropower station to form a sample set, construct a neural network model suitable for regression tasks, input the prepared sample set into the network for training. When it is necessary to query the output value at a specific calculated water head, input the water head value into the trained neural network, and output the corresponding minimum and maximum output prediction values.
[0036] In another optional embodiment, the first calculation formula in step B2 can also adopt an analytical calculation method based on a parametric formula. According to the physical working principle of the water turbine generator, deduce the theoretical formula for the water head power generation capacity, and use the characteristic curve of the water turbine (usually provided by the manufacturer) or a large amount of historical operation data to determine the key parameters in the above theoretical formula. Incorporate the fitted parametric formula into the calculation program. When the output at a certain water head is required, the system calls these analytical formulas, substitutes the current water head value into the corresponding formula, and outputs the calculated result.
[0037] It should be noted that the head of the hydropower station changes in real time, and the stable operation area of the unit strictly depends on the current head; the traditional method uses the parameters corresponding to a fixed head or manual estimation, which has a large error. In this step, through interpolation calculation, it is ensured that no matter how the head changes, the system can obtain the upper and lower limits of the stable area that best matches the current working condition.
[0038] In step S3, determine the total output. According to the output range of the unit and combined with the number of units, output the unit output distribution plan, including the following steps C1 - C3: C1: Determine the total output value according to the production demand.
[0039] C2: Make a first judgment according to the total output value, combined with the output range of the unit and the number of units under the current head value.
[0040] C3: Output the unit output distribution plan according to the result of the first judgment.
[0041] In the implementation manner of this application, the specific steps of making the first judgment in step C2 are as follows: Calculate that when a given head H is fixed, according to the number of operating units NO, determine the specific given output PZ and how to distribute it to NO units, where the maximum number of units in the power station is NO_MAX.
[0042] If (PZ / Pmax) > NO_MAX, it means that the output PZ has exceeded the maximum total output of the units under the given head of this power station, and an alarm needs to be given, indicating that this output is illegal.
[0043] If (PZ / NO) > Pmax, it means that the output PZ has exceeded the maximum output value of the given number of operating units NO, and an alarm reminder is needed, indicating that it is lower than the minimum number of starting units, and units need to be added.
[0044] If (PZ / NO) < PLMin, it means that the output PZ is lower than the lowest value of all restricted operation areas of the given number of operating units NO, and an alarm reminder is needed, indicating that units need to be arranged to stop or arranged to run no-load for a short time. And when NO = 1 and PZ < Pmin, directly remind that there is only 1 unit and it is lower than the lower limit of the stable operation area.
[0045] If Pmin <= (PZ / NO) <= Pmax, it means that when the output PZ is distributed according to the number of operating units NO, it can be completely distributed to the stable operation area. Set PSET_AVE to represent the average output of the units in the stable operation area, and directly calculate PSET_AVE = PZ / NO, and output the unit output distribution plan.
[0046] In an optional implementation, the first judgment output unit allocation scheme in step C2 can also adopt a sequential load allocation method based on dynamic programming, which defines stages, states and decisions, regards each unit as a stage, and the state of each stage is the total allocated load. In the kth stage, how much load is allocated to the kth unit? The Bellman equation is established and recursively derived from back to front, and the decision is made from front to back to output the allocation scheme.
[0047] In another optional implementation, the first judgment output unit allocation scheme in step C2 can also adopt the classic allocation method based on the principle of constant incremental rate. For each unit, its input-output characteristic curve under the fixed head is established through test or historical operating data. The constant incremental rate equation is constructed based on the input-output characteristic curve, and the constant incremental rate equation is solved to output the allocation scheme.
[0048] It should be noted that the present invention, through the first judgment, determines whether the total output exceeds the maximum capacity of the power station or the maximum carrying capacity of the current number of units in operation, thereby avoiding the issuance of impossible tasks. It also determines whether the output of a single unit is too low after average distribution, which may lead to entering the restricted area or even failure to operate, and provides early warning that the number of operating units needs to be reduced, thus preventing equipment from being forced to operate in dangerous conditions or system control failure due to unreasonable instructions.
[0049] In step S4, a second calculation is performed based on the unit output allocation scheme, and the latest output allocation scheme is output, including the following steps D1-D3: D1: Based on the unit output allocation scheme, construct a second calculation formula to calculate the adjustable output value of the current allocation scheme.
[0050] D2: Based on the adjustable output value of the current allocation scheme, construct a second judgment formula to determine whether the unit output allocation scheme is feasible.
[0051] The second judgment formula is expressed as follows: Where PZ is the given total output, and F(n,NO) represents the minimum output value of various combinations of units operating in the stable and restricted operating regions when the number of operating units NO is given. This represents the adjustable output value of various combinations of units in the stable and restricted operating regions when the number of operating units NO is constant.
[0052] Calculate the total output and the number of units in operation. When F(n,NO)<=PZ<=F(n+1,NO), return the value of n to obtain the unit operation combination.
[0053] If total output PZ> + This indicates that the given output PZ cannot be fully allocated to the stable zone and the restricted operating zone, prompting an alarm. It means that with a given total output and a given number of units in operation, it is impossible to completely avoid the prohibited operating zone, and the number of units in operation needs to be adjusted.
[0054] If total output PZ <= + This indicates that the given output PZ can be allocated to a combination of the stable region and the restricted operating region, and the allocation scheme is feasible.
[0055] D3: Based on the currently feasible unit output allocation scheme, perform adjustable capacity proportional allocation and output the latest output allocation scheme.
[0056] Based on the currently feasible unit output allocation scheme, calculate the adjustable output of units in the stable zone: P_WDQ_T = (Pmax - Pmin) * N, and the adjustable output of units in the restricted operating zone: P_XZQ_T = (PLMax - PLMin) * (NO - N). Allocate the power output proportionally to the adjustable capacity. Average power output of units in the stable operating range: PZ_WDQ =(PZ-F(N,NO))*P_WDQ_T / F_T(N,NO) / N+Pmin Average power output distribution of units in restricted operating areas: PZ_XZQ =(PZ-F(N,NO))*P_XZQ_T / F_T(N,NO) / (N0-N)+PLMin In this embodiment of the application, the second calculation formula in step D1 uses a binary arithmetic series algorithm, and the formula is expressed as follows: Where F(n,NO) represents the minimum output value of various combinations of units operating in the stable and restricted operating regions when the number of operating units NO is given, and NO is the total number of operating units. This represents the adjustable output value of various combinations of units operating in the stable and restricted operating regions, given a fixed number of operating units (NO). It also represents the minimum output value of various combinations of units operating in the stable and restricted operating regions. This represents the lower limit of the restricted operating area, where n is the number of units operating within the restricted operating area. To limit the upper limit of the operating area, This represents the upper limit of the stable region for a single machine. This is the lower limit of the stable region for a single machine.
[0057] It should be noted that when When F = 1, F(1, NO) indicates that one unit is operating in the stable region, and NO-1 units are operating in the restricted operating region. The value of F(1, NO) is the minimum output value under this combination. Because the entire plant's AGC (Automatic Generation Control) must not disengage, at least one unit is required to operate in the stable region. The range of values is 1 to NO. For a given operating unit NO, the range of F(n,NO) and F(n+1,NO) is a constant. Therefore, F(n,NO) is actually a binary arithmetic series.
[0058] In an alternative implementation, the second calculation formula in step D1 can also be a real-time standby capacity assessment method based on sensitivity analysis, which treats the current load allocation scheme as the optimal solution to an optimization problem, checks which units' output is exactly at their upper or lower limits, calculates the adjustable range, and outputs the adjustable output value under the current optimal allocation scheme.
[0059] In another alternative implementation, the second calculation formula in step D1 can also be based on the feasible region traversal method of linear programming, reconstructing the optimization model, taking the output of each unit as the decision variable, establishing a new objective function based on the maximum upper reserve and the maximum lower reserve of the system, and solving the function to output the adjustable range.
[0060] It should be noted that the present invention calculates the adjustable output value of the allocation scheme using a binary arithmetic series, clarifying the range within which the total output can be adjusted without loss under a selected unit combination. By systematically traversing the allocation combinations and using a binary arithmetic series for calculation and judgment, the algorithm can automatically locate a feasible unit combination that can bear the current total output.
[0061] In step S5, the latest power allocation scheme is optimized by performing a third calculation and outputting the optimized allocation scheme, including the following steps E1-E3: E1: Perform a second optimization on the latest power allocation scheme, and a third calculation to obtain the optimal power output point under the current allocation scheme's restricted operating area.
[0062] E2: Construct a secondary optimization allocation formula based on the optimal output point in the operating area under the current allocation scheme.
[0063] It should be noted that while the adjustable capacity proportional allocation algorithm guarantees that power output will be allocated to the corresponding range, the restricted operating zone is the minimum intersection of the lower limit operating zone and the full head operating zone. The vibration and sway characteristics of the units operating at different positions within this zone also vary significantly. The optimal point should be determined based on the input head value. Therefore, the power output allocated to the units in the restricted operating zone using the above algorithm may not be the optimal solution. Thus, a secondary optimization allocation of the corresponding power output is required.
[0064] Construct a quadratic optimization allocation formula: Where Hmax and Hmin are the highest and lowest head values, respectively, H is the input head value, and PZ_XZQ_E is the optimal output point in the limit operation zone.
[0065] E3: Output the optimized allocation scheme based on the quadratic optimization allocation formula.
[0066] Based on the quadratic optimization allocation formula, the calculation of the quadratic optimization allocation amount can be divided into four cases: (1) If (PZ_XZQ -PZ_XZQ_E) * (NO - N)<= (Pmax - PZ_WDQ) * N, The average output value of the secondary distribution of the unit in the restricted operating area is then: PZ_XZQ_S=PZ_XZQ_E Average output value of the unit under secondary distribution in the stable operating range: PZ_WDQ_S=PZ_WDQ+(PZ_XZQ -PZ_XZQ_E) * (NO - N) / N PZ_WDQ_S represents the average output value of the unit under secondary distribution in the stable operating area.
[0067] (2) If (PZ_XZQ_E-PZ_XZQ ) * (NO - N)<= ( PZ_WDQ-Pmin) * N, The average output value of the secondary distribution of the unit in the restricted operating area is then: PZ_XZQ_S=PZ_XZQ_E Average output value of the unit under secondary distribution in the stable operating range: PZ_WDQ_S=PZ_WDQ -(PZ_XZQ_E-PZ_XZQ ) * (NO - N) / N (3) If (PZ_XZQ -PZ_XZQ_E) * (NO - N)>= (Pmax - PZ_WDQ) * N, The average output value of the secondary distribution of the unit in the restricted operating area is then: PZ_XZQ_S=PZ_XZQ - (Pmax - PZ_WDQ) *N / (NO -N) Average output value of the unit under secondary distribution in the stable operating range: PZ_WDQ_S=Pmax (4) If (PZ_XZQ_E-PZ_XZQ ) * (NO - N)>= (PZ_WDQ-Pmin ) * N, The average output value of the secondary distribution of the unit in the restricted operating area is then: PZ_XZQ_S=PZ_XZQ + (PZ_WDQ-Pmin ) *N / (NO -N) Average output value of the unit under secondary distribution in the stable operating range: PZ_WDQ_S=Pmin In this embodiment of the application, the third calculation formula in step E1 is expressed as follows: Wherein, PZ_XZQ_E represents the optimal output point under the current allocation scheme's restricted operating area, H is the current head value, and Hmax and Hmin are the highest and lowest head values, respectively.
[0068] In an optional implementation, the third calculation in step E1 can also employ a swarm intelligence search method based on particle swarm optimization, which initializes and encodes particles, defines fitness functions and handles constraints, and iteratively updates and converges to obtain the optimal output point under the current allocation scheme's restricted operating area.
[0069] In another optional implementation, the third calculation in step E1 can also adopt the optimal output point calculation scheme based on function transformation and generalized inverse. This algorithm eliminates inequality constraints (such as limiting the operating range and oscillation safety threshold) through function transformation, and then uses generalized inverse to solve the indeterminate equation system to output the optimal output point under the limited operating range.
[0070] It should be noted that this invention calculates a theoretically optimal output point for the units within the restricted operating zone by constructing a third calculation formula. It introduces four judgment scenarios: when the units in the stable zone have sufficient adjustment margin, the units in the restricted zone are adjusted by increasing or decreasing the adjustment margin to make them reach the optimal output point; when the units in the stable zone have insufficient adjustment margin, they are allowed to operate at the boundary, thereby making the units in the restricted zone as close to the optimal output point as possible. Under the premise that the total output and the number of units in operation remain unchanged, the optimal allocation scheme for the equipment is found, reducing the number of start-ups and shutdowns.
[0071] In step S6, the optimized allocation scheme is applied, the combined vibration zone is calculated in real time, and the issued load curve is checked to see if it falls within the combined vibration zone. Appropriate adjustments are made, including the following steps F1-F3: F1: Calculates the number of combined vibration zones in real time based on the input head value and the single-machine stability zone.
[0072] The formula for calculating the number of combined vibration zones is as follows: Wherein, N represents the presence of N joint vibration zones.
[0073] F2: Calculate the range of each joint vibration zone based on the number of joint vibration zones.
[0074] The formula for calculating the range of each joint vibration zone is expressed as follows: Where ZDQ(k).PMIN is the lower limit of the joint vibration zone of the kth segment, ZDQ(k).PMAX is the upper limit of the joint vibration zone of the kth segment, and k is a positive integer from 1 to N.
[0075] F3: Check whether the issued load curve falls within the combined vibration zone and make corresponding adjustments.
[0076] Check whether the load curve values issued by the dispatcher fall within the joint vibration zone, so as to apply to the dispatcher for modification in advance.
[0077] In summary, this invention constructs an intelligent load allocation decision system based on operating range definition, head adaptive calculation, joint vibration zone identification, binary arithmetic series optimization, and secondary optimization allocation. This system enables full-process control of mixed-flow turbine units, from operating condition perception, safety boundary calculation, output allocation to vibration control. It solves the problems of low calculation accuracy and efficiency, frequent unit start-ups and shutdowns, and high equipment wear in traditional manual allocation methods, thereby improving the automation level of hydropower station operation.
[0078] Example 3 illustrates a schematic scheme for a method of calculating output distribution to reduce the number of start-ups and shutdowns of a mixed-flow turbine unit based on a binary arithmetic series algorithm. It should be noted that the technical solution of this system for calculating output distribution to reduce the number of start-ups and shutdowns of a mixed-flow turbine unit based on a binary arithmetic series algorithm is based on the same concept as the aforementioned method for calculating output distribution to reduce the number of start-ups and shutdowns of a mixed-flow turbine unit based on a binary arithmetic series algorithm. Details not described in detail in this embodiment can be found in the description of the aforementioned method for calculating output distribution to reduce the number of start-ups and shutdowns of a mixed-flow turbine unit based on a binary arithmetic series algorithm.
[0079] This embodiment also provides a power distribution calculation system for reducing the number of start-ups and shutdowns of mixed-flow turbine units based on a binary arithmetic series algorithm, including an initialization module, a first calculation module, a scheme output module, a second calculation module, a third calculation module, and an adjustment module; The initialization module is used to define the single-unit stable zone and restricted operating zone of the power station turbine unit through actual power station experiments; The first calculation module is used to input the head value and perform a first calculation based on the single-unit stable zone and restricted operating zone of the turbine unit to obtain the unit output range under the input head value; The scheme output module is used to determine the total output and output the unit output allocation scheme based on the unit output range and the number of units. The second calculation module is used to perform a second calculation based on the unit output allocation scheme and output the latest output allocation scheme; The third calculation module is used to optimize the latest power allocation scheme, perform a third calculation, and output the optimized allocation scheme. The adjustment module is used to apply the optimized allocation scheme, calculate the combined vibration zone in real time, check whether the issued load curve falls into the combined vibration zone, and make corresponding adjustments.
[0080] This embodiment also provides a computer device, including a memory and a processor. The memory stores a computer program. The processor executes the computer program to implement the steps of a method for calculating the output distribution of a mixed-flow turbine unit based on a binary arithmetic series algorithm to reduce the number of start-ups and shutdowns.
[0081] This embodiment proposes a computer-readable storage medium storing a computer program, characterized in that, when the computer program is executed by a processor, it implements the steps of a method for calculating the output distribution of a mixed-flow turbine unit based on a binary arithmetic series algorithm to reduce the number of start-ups and shutdowns.
[0082] Based on the above description of the implementation methods, those skilled in the art can clearly understand that the present invention can be implemented using software and necessary general-purpose hardware, and of course, it can also be implemented using hardware. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as a computer floppy disk, read-only memory (ROM), random access memory (RAM), flash memory, hard disk, or optical disk, etc., including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods of the various embodiments of the present invention.
[0083] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.
Claims
1. A power distribution calculation method for reducing the number of start-stop operations of a Francis turbine unit based on a binary arithmetic progression algorithm, characterized by, The application relates to a method for calculating a power distribution scheme of a power station. The method comprises the following steps: defining a single-machine stable region and a limited operation region of a water turbine unit of the power station through actual experiments of the power station; inputting a water head value, and performing first calculation to obtain a unit output range under the input water head value according to the single-machine stable region and the limited operation region of the water turbine unit; determining a total output, and outputting a unit output distribution scheme according to the unit output range and the number of units; performing second calculation according to the unit output distribution scheme, and outputting a latest output distribution scheme; optimizing the latest output distribution scheme, performing third calculation, and outputting an optimized distribution scheme; 2. The output distribution calculation method for reducing the number of start and stop of a Francis turbine unit based on a binary arithmetic progression algorithm according to claim 1, characterized in that, applying the optimized distribution scheme to perform real-time calculation on a joint vibration region, checking whether a load curve issued falls into the joint vibration region, and performing corresponding adjustment. The method for defining the single-machine stable region and the limited operation region of the water turbine unit of the power station through actual experiments of the power station comprises the following steps: assigning maximum and minimum values of the single-machine stable region of a single unit under different water heads obtained through actual experiments of the power station, and outputting an operation parameter table of the water turbine unit under different water head values; obtaining the single-machine stable region according to the operation parameter table of the water turbine unit under different water head values; 3. The output distribution calculation method for reducing the number of start-up and shut-down of a Francis turbine unit based on a binary arithmetic progression algorithm according to claim 2, characterized in that, assigning maximum and minimum values of the limited operation region under full water head obtained through experiments, and obtaining the limited operation region. The method for inputting the water head value, performing first calculation to obtain the unit output range under the input water head value according to the single-machine stable region and the limited operation region of the water turbine unit comprises the following steps: inputting the water head value, and establishing a first calculation formula; performing output range calculation according to the first calculation formula and the single-machine stable region and the limited operation region of the water turbine unit; 4. The output distribution calculation method for reducing the number of start-up and shut-down of a Francis turbine unit based on a binary arithmetic progression algorithm according to claim 3, characterized in that, outputting the unit output range under the current water head value. The method for determining the total output, outputting the unit output distribution scheme according to the unit output range and the number of units comprises the following steps: determining a total output value according to production demand; performing first judgment according to the total output value, the unit output range under the current water head value and the number of units; 5. The output distribution calculation method for reducing the number of start-up and shut-down of a Francis turbine unit based on a binary arithmetic progression algorithm according to claim 4, characterized in that, outputting the unit output distribution scheme according to the first judgment result. The method for performing second calculation and judging whether the unit output distribution scheme is feasible comprises the following steps: Wherein, F(n, NO) represents the minimum output value of various combinations of units running in the stable region and the limit operation region when the given number of operating units is NO, and NO is the total number of operating units, represents the adjustable output value of various combinations of units running in the stable region and the limit operation region when the given number of operating units is NO, and the minimum output value of various combinations of units running in the stable region and the limit operation region, is the lower limit value of the limit operation region, n is the number of units running in the limit operation region, is the upper limit value of the limit operation region, is the upper limit value of the single-unit stable region, is the lower limit value of the single-unit stable region; constructing a second calculation formula to calculate an adjustable output value of the current distribution scheme according to the unit output distribution scheme; constructing a second judgment formula to judge whether the unit output distribution scheme is feasible according to the adjustable output value of the current distribution scheme; 6. The output distribution calculation method for reducing the number of start-up and shut-down of a Francis turbine unit based on a binary arithmetic progression algorithm according to claim 5, characterized in that, performing adjustable capacity proportional distribution according to the currently feasible unit output distribution scheme, and outputting a latest unit output distribution scheme. The method for optimizing the latest output distribution scheme, performing third calculation and outputting the optimized distribution scheme comprises the following steps: performing second optimization on the latest output distribution scheme, performing third calculation and obtaining a best output point under the limited operation region of the current distribution scheme: wherein PZ_XZQ_E represents the best output point under the limited operation region of the current distribution scheme, H is the current water head value, Hmax and Hmin are respectively the highest water head value and the lowest water head value; constructing a second optimization distribution formula according to the best output point under the limited operation region of the current distribution scheme; 7. The output distribution calculation method for reducing the number of start-up and shut-down of a Francis turbine unit based on a binary arithmetic progression algorithm according to claim 6, characterized in that, outputting the optimized distribution scheme according to the second optimization distribution formula. The method for performing real-time calculation on the joint vibration region comprises the following steps: performing real-time calculation on the number of joint vibration regions according to the current water head value and the single-machine stable region: Wherein, N represents N joint vibration zones; According to the number of joint vibration zones, the range of each joint vibration zone is calculated: Wherein, ZDQ(k).PMIN is the lower limit value of the kth joint vibration zone, ZDQ(k).PMAX is the upper limit value of the kth joint vibration zone, and k is a positive integer from 1 to N.
8. A power distribution calculation system for reducing the number of start and stop of a Francis turbine unit based on a binary arithmetic progression algorithm, applying the method according to any one of claims 1 to 7, characterized in that, The method comprises an initialization module, a first calculation module, a scheme output module, a second calculation module, a third calculation module and an adjustment module. The initialization module is used to define the single-machine stable zone and the limited operation zone of the hydroelectric generating set through actual experiments of the power station. The first calculation module is used to input the water head value, and perform first calculation to obtain the unit output range under the input water head value according to the single-machine stable zone and the limited operation zone of the hydroelectric generating set. The scheme output module is used to determine the total output, and output the unit output distribution scheme according to the unit output range and the number of units. The second calculation module is used to perform second calculation according to the unit output distribution scheme, and output the latest output distribution scheme. The third calculation module is used to optimize the latest output distribution scheme, perform third calculation, and output the optimized distribution scheme. The adjustment module is used to apply the optimized distribution scheme, calculate the joint vibration zone in real time, check whether the issued load curve falls into the joint vibration zone, and make corresponding adjustment. 9.A computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer device is configured to perform the method according to any one of claims 1-8 when the computer program is executed by the processor. The processor executes the computer program to realize the steps of the power distribution calculation method for reducing the start-stop times of the Francis turbine set based on the binary arithmetic progression algorithm in any one of claims 1 to 7.
10. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to realize the steps of the power distribution calculation method for reducing the start-stop times of the Francis turbine set based on the binary arithmetic progression algorithm in any one of claims 1 to 7.