CN114169178B - Transient process pressure pulsation numerical processing system and method based on EMD - Google Patents
Transient process pressure pulsation numerical processing system and method based on EMD Download PDFInfo
- Publication number
- CN114169178B CN114169178B CN202111529103.3A CN202111529103A CN114169178B CN 114169178 B CN114169178 B CN 114169178B CN 202111529103 A CN202111529103 A CN 202111529103A CN 114169178 B CN114169178 B CN 114169178B
- Authority
- CN
- China
- Prior art keywords
- pressure
- calculation
- transient
- unit
- pulsation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000010349 pulsation Effects 0.000 title claims abstract description 42
- 230000001052 transient effect Effects 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012545 processing Methods 0.000 title claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 65
- 238000004458 analytical method Methods 0.000 claims abstract description 32
- 238000012360 testing method Methods 0.000 claims abstract description 28
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 17
- 238000000611 regression analysis Methods 0.000 claims abstract description 10
- 238000013461 design Methods 0.000 claims abstract description 9
- 238000013178 mathematical model Methods 0.000 claims abstract description 5
- 238000005457 optimization Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000011144 upstream manufacturing Methods 0.000 claims description 5
- 238000003672 processing method Methods 0.000 claims description 4
- 230000001105 regulatory effect Effects 0.000 claims description 4
- 230000000630 rising effect Effects 0.000 claims description 4
- 230000002411 adverse Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 abstract 1
- 230000002349 favourable effect Effects 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control Of Non-Positive-Displacement Pumps (AREA)
Abstract
The invention provides a transient process pressure pulsation value processing system and a method thereof based on EMD, wherein a regression analysis module collects pressure pulsation signal values in transient debugging tests of a unit, extracts IMF components, takes the sum of the IMF components as pulsation pressure and takes the final allowance as mean pressure; performing constant flow calculation analysis and unsteady flow calculation analysis to obtain a mathematical model closest to a debugging test; obtaining a mean pressure calculation process, a unit rotating speed and a surge chamber water level fluctuation process and an extremum; and the presumption contrast calculation analysis module is used for carrying out one-dimensional transient process calculation to obtain a calculation mean pressure and a calculation rotating speed value, so that parameters under the least favorable transient working condition can be obtained, if the extreme values of the parameters can not fully meet the design requirements, the parameters are optimized, and the next round of presumption contrast calculation is carried out until all the parameters meet the design requirements. The invention forms a complete set of pressure pulsation regression analysis and presumption contrast calculation method which is suitable for transient test of the unit.
Description
Technical Field
The invention relates to the field of hydropower station unit pressure pulsation analysis, in particular to a transient process pressure pulsation numerical processing system and method based on empirical mode analysis (EMD).
Background
In the transient process of the unit, the instability of the unit is increased, the safety problem exists due to the fact that the pressure pulsation is too high, and the instability of the pressure pulsation and the like can also occur when the unit is put into operation for the first time. Therefore, a unit debugging test is needed, the pressure pulsation value is predicted, whether the pressure pulsation value is in a safe range or not is judged, and potential safety hazards are reduced. The invention reduces the error of analyzing the pressure pulsation value and analyzes the running environment of the unit in a safe range.
Disclosure of Invention
The invention provides a transient process pressure pulsation numerical processing system based on empirical mode analysis (EMD), which comprises the following components: the regression analysis module and the presumption contrast calculation analysis module; wherein,
the regression analysis module is used for collecting pressure pulsation signal values in a transient debugging test of the unit, extracting IMF components from the volute pressure pulsation signal and the draft tube inlet pressure pulsation signal by using an empirical mode analysis method, taking the sum of the IMF components as pulsation pressure and taking the final allowance as average pressure; and adopting the same calculation conditions as transient test operation, including water level, guide vane opening, unit output and the like, and carrying out constant flow calculation analysis and non-constant flow calculation analysis. Parameter correction is carried out through constant flow calculation and analysis, and a mathematical model closest to a debugging test is obtained; and obtaining the mean pressure calculation, the unit rotating speed and the surge chamber water level fluctuation process and the extremum through unsteady flow calculation and analysis.
The speculation contrast calculation analysis module adopts the same calculation conditions as the transient debugging test working condition of the unit to be predicted, including water level, guide vane opening, unit output and the like, performs one-dimensional transient process calculation to obtain calculated mean pressure and calculated rotating speed values, namely the maximum and minimum pressures of an upstream side pressure pipeline, a volute inlet and a draft tube inlet of the water inlet ball valve under the most unfavorable transient working condition, the maximum rising value of the rotating speed, the maximum and minimum surge value of a pressure regulating chamber, the distribution of the maximum pressure and the minimum pressure of the whole water delivery system along a pipeline and the like, and if the extreme value of the parameters can not fully meet the design requirement, guide vane closing rule optimization or appropriate operation condition optimization measures are required to be provided on the basis, and the next round of speculation contrast calculation is performed until all the parameters meet the design requirement.
Correspondingly, by adopting the system, the invention also provides an EMD-based transient process pressure pulsation numerical processing method,
the method is characterized by comprising the following steps of:
the regression analysis module is used for collecting pressure pulsation signal values in transient debugging test of the unit, extracting IMF components from volute pressure pulsation signals and draft tube inlet pressure pulsation signals by using an empirical mode analysis method, taking the sum of the IMF components as pulsation pressure and taking the final allowance as average pressure; adopting the same calculation conditions as the transient test operation to perform constant flow calculation analysis and unsteady flow calculation analysis; parameter correction is carried out through constant flow calculation and analysis, and a mathematical model closest to a debugging test is obtained; calculating and analyzing the unsteady flow to obtain a mean pressure calculation, a unit rotating speed and a surge chamber water level fluctuation process and an extremum;
the prediction contrast calculation analysis module adopts the same calculation conditions as the transient debugging test working condition of the unit to be predicted to perform one-dimensional transient process calculation to obtain a calculation mean pressure and a calculation rotating speed value, namely the maximum and minimum pressures of an upstream side pressure pipeline, a volute inlet and a draft tube inlet of the water inlet ball valve under the least adverse transient working condition, the maximum rising value of the rotating speed, the maximum and minimum surge value of the pressure regulating chamber and the distribution of the maximum pressure and the minimum pressure of the whole water delivery system along the pipeline are obtained, if the extreme value of the parameters can not fully meet the design requirements, guide vane closing rule optimization is required to be performed or proper operation condition optimization measures are provided on the basis, and the next round of prediction contrast calculation is performed until all the parameters meet the design requirements.
Wherein, the calculation formula of the middle-dimension transient process is as follows:
(1) predictive mean pressure calculation formula
Predicted mean pressure = calculated mean pressure at the corresponding time + calculated error at the corresponding time.
(2) Calculation formula for predicting instantaneous total pressure
Predicted instantaneous total pressure = predicted mean pressure at the respective moment + pulsating pressure at the respective moment.
(3) Calculation formula for predicting rotating speed of unit
Predicted unit speed = corresponding time computer unit speed + corresponding time unit speed calculation error.
The invention provides a transient process pressure pulsation numerical value processing system and a transient process pressure pulsation numerical value processing method based on empirical mode analysis, which form a complete set of pressure pulsation regression analysis and presumption contrast calculation methods applicable to transient tests of units, and correct and perfect pressure pulsation numerical value calculation errors; the invention also provides a new thought and method for inversion and predictive analysis of transient tests of similar units, has reference significance for design selection of units, calculation of transient processes and determination of control standards, and is beneficial to promoting construction and perfection of domestic unit localization processes and related standard specification systems.
Drawings
Fig. 1 is a schematic structural view of the present invention.
Detailed Description
The specific technical scheme of the invention is described by combining the embodiments.
According to the steps shown in fig. 1:
test environment: when the double-motor is normally operated under the load of 225MW, 75% of rated load is thrown at the same time, the upstream water level 397.28, the downstream water level 76.88, the guide vanes are normally closed, and the ball valves are closed. Arranging a pressure sensor, collecting a pressure pulsation value, decomposing test data by using an empirical mode analysis (EMD) method through 9 times of EMD, and obtaining the frequency range and the maximum amplitude of each IMF. The sum of the IMF components of each order is taken as pressure pulsation, and the final margin is taken as mean pressure. According to EMD decomposition, it can be seen that for ball valve and volute inlet measuring points, the maximum energy of the first machine appears in the middle and low frequency bands with the frequency of 30-50 Hz; and for the draft tube measuring point, the frequency band of the maximum energy of the first machine is lower between 20 and 80 Hz. And carrying out one-dimensional hydraulic transition rechecking calculation according to the to-be-predicted double-engine load dump 100% test working condition by adopting a field actually-measured guide vane closing rule, and then carrying out calculation error correction and pressure pulsation superposition on the basis of a calculated value according to the completed single-engine load dump 100% load and double-engine load dump 75% load test regression analysis result to respectively obtain the average pressure of the to-be-predicted test working condition and the instantaneous pressure extreme value containing the pressure pulsation. According to the test result of the double-machine load shedding 100% of the hydraulic unit of the first machine, the actual measurement results and the predicted results of the ball valve inlet, the volute inlet and the draft tube inlet measuring points of the double-machine load shedding of the first machine and the second machine are respectively compared, and the predicted value and the actual measurement value have good consistency no matter the pressure change process or the average pressure and the instantaneous pressure value.
Claims (3)
1. An EMD-based transient process pressure pulsation numerical processing system, comprising: the regression analysis module and the presumption contrast calculation analysis module;
the regression analysis module is used for: collecting pressure pulsation signal values in transient debugging test of the unit, respectively extracting IMF components from volute pressure pulsation signals and draft tube inlet pressure pulsation signals by using an empirical mode analysis method, taking the sum of the IMF components as pulsation pressure, and taking the final allowance as average pressure; adopting the same calculation conditions as the transient test operation to perform constant flow calculation analysis and unsteady flow calculation analysis; parameter correction is carried out through constant flow calculation and analysis, and a mathematical model closest to a debugging test is obtained; calculating and analyzing the unsteady flow to obtain a mean pressure calculation, a unit rotating speed and a surge chamber water level fluctuation process and an extremum;
the speculation contrast calculation analysis module is used for: and calculating a one-dimensional transient process by adopting calculation conditions which are the same as the transient debugging test working conditions of the unit to be predicted, obtaining a calculated mean pressure and a calculated rotating speed value, namely obtaining the maximum and minimum pressures of an upstream side pressure pipeline, a volute inlet and a draft tube inlet of the water inlet ball valve under the least adverse transient working conditions, and the maximum rising value of the rotating speed, the maximum and minimum surge value of a pressure regulating chamber and the distribution of the maximum pressure and the minimum pressure of the whole water delivery system along a pipeline.
2. An EMD-based transient process pressure pulsation value processing method, characterized in that the processing system of claim 1 is adopted, comprising the following processes:
the regression analysis module is used for collecting pressure pulsation signal values in transient debugging test of the unit, extracting IMF components from volute pressure pulsation signals and draft tube inlet pressure pulsation signals by using an empirical mode analysis method, taking the sum of the IMF components as pulsation pressure and taking the final allowance as average pressure; adopting the same calculation conditions as the transient test operation to perform constant flow calculation analysis and unsteady flow calculation analysis; parameter correction is carried out through constant flow calculation and analysis, and a mathematical model closest to a debugging test is obtained; calculating and analyzing the unsteady flow to obtain a mean pressure calculation, a unit rotating speed and a surge chamber water level fluctuation process and an extremum;
the prediction contrast calculation analysis module adopts the same calculation conditions as the transient debugging test working condition of the unit to be predicted to perform one-dimensional transient process calculation to obtain a calculation mean pressure and a calculation rotating speed value, namely the maximum and minimum pressures of an upstream side pressure pipeline, a volute inlet and a draft tube inlet of the water inlet ball valve under the least adverse transient working condition, the maximum rising value of the rotating speed, the maximum and minimum surge value of the pressure regulating chamber and the distribution of the maximum pressure and the minimum pressure of the whole water delivery system along the pipeline are obtained, if the extremum of the parameters can not fully meet the design requirements, guide vane closing rule optimization or operation condition optimization measures are required to be provided on the basis, and the next round of prediction contrast calculation is performed until all the parameters meet the design requirements.
3. The EMD-based transient pressure pulsation value processing method according to claim 2, wherein the one-dimensional transient calculation formula is:
(1) predictive mean pressure calculation formula
Predicted mean pressure = calculated mean pressure at the respective time + calculated error at the respective time;
(2) calculation formula for predicting instantaneous total pressure
Predicted instantaneous total pressure = predicted mean pressure at the respective moment + pulsating pressure at the respective moment;
(3) calculation formula for predicting rotating speed of unit
Predicted unit speed = corresponding time computer unit speed + corresponding time unit speed calculation error.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111529103.3A CN114169178B (en) | 2021-12-14 | 2021-12-14 | Transient process pressure pulsation numerical processing system and method based on EMD |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111529103.3A CN114169178B (en) | 2021-12-14 | 2021-12-14 | Transient process pressure pulsation numerical processing system and method based on EMD |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114169178A CN114169178A (en) | 2022-03-11 |
| CN114169178B true CN114169178B (en) | 2024-03-22 |
Family
ID=80486529
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111529103.3A Active CN114169178B (en) | 2021-12-14 | 2021-12-14 | Transient process pressure pulsation numerical processing system and method based on EMD |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114169178B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009041851A1 (en) * | 2007-09-24 | 2009-04-02 | Central Institute Of Aviation Motors (Ciam) | Method for monitoring the operating modes of a compressor and a device for carrying out said method |
| CN111104734A (en) * | 2019-12-12 | 2020-05-05 | 中国电建集团中南勘测设计研究院有限公司 | Inversion prediction method for pumped storage power station unit load shedding test |
| CN112528805A (en) * | 2020-12-04 | 2021-03-19 | 湖南五凌电力科技有限公司 | Water turbine pressure pulsation analysis method, device, equipment and storage medium |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10983507B2 (en) * | 2016-05-09 | 2021-04-20 | Strong Force Iot Portfolio 2016, Llc | Method for data collection and frequency analysis with self-organization functionality |
-
2021
- 2021-12-14 CN CN202111529103.3A patent/CN114169178B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009041851A1 (en) * | 2007-09-24 | 2009-04-02 | Central Institute Of Aviation Motors (Ciam) | Method for monitoring the operating modes of a compressor and a device for carrying out said method |
| CN111104734A (en) * | 2019-12-12 | 2020-05-05 | 中国电建集团中南勘测设计研究院有限公司 | Inversion prediction method for pumped storage power station unit load shedding test |
| CN112528805A (en) * | 2020-12-04 | 2021-03-19 | 湖南五凌电力科技有限公司 | Water turbine pressure pulsation analysis method, device, equipment and storage medium |
Non-Patent Citations (2)
| Title |
|---|
| 水电站水力过渡过程数值仿真的理论与实践;巨江;诸亮;刘菁;徐春燕;;西北水电;20090630(第03期);全文 * |
| 黑麇峰抽水蓄能电站机组甩负荷试验反演预测及主要特性分析;郑建兴;刘平;杨晖;曾艳梅;;水电站机电技术;20161215(第S2期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114169178A (en) | 2022-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105240187B (en) | The method that turbine state monitoring and fault diagnosis are realized based on chaology | |
| CN109190166B (en) | Cavitation judgment and state evaluation method and system for vane pump | |
| CN113236581A (en) | Intelligent parallel pump system and optimal adjustment method | |
| CN111159844A (en) | Abnormity detection method for exhaust temperature of gas turbine of power station | |
| CN110702437B (en) | Evaluation method for health condition of hydraulic system of hydroelectric generating set | |
| CN112052871B (en) | Rocket engine health diagnosis method and system based on support vector machine | |
| CN114169178B (en) | Transient process pressure pulsation numerical processing system and method based on EMD | |
| CN110439825A (en) | A method of on-line monitoring sealing ring gap | |
| CN108959739B (en) | Method and device for analyzing pressure pulsation of pumped storage power station in transition process | |
| Guo et al. | Application of qualitative trend analysis in fault diagnosis of entrained-flow coal-water slurry gasifier | |
| CN108681614B (en) | Turbofan engine mutation fault diagnosis method based on improved Gaussian particle filtering | |
| CN113705120A (en) | Method for formulating optimal regulation and control scheme of rear valve of water transfer engineering pump | |
| CN113446146A (en) | Online water turbine efficiency test method | |
| CN103499452A (en) | Automatic interpretation method for engine testing data | |
| CN109344522B (en) | Method and system for calculating grinding quantity of stationary blade of axial flow compressor | |
| CN115342101A (en) | Hydraulic system internal leakage detection method and system for hydraulic generator speed regulator | |
| CN111238582A (en) | Method for measuring flow of draft tube of pumped storage unit | |
| CN114997309A (en) | Water feed pump fault early warning method and device | |
| CN112761740B (en) | Turbine blade fault testing system and intelligent fault diagnosis method thereof | |
| CN211777784U (en) | Unit efficiency monitoring device | |
| CN110705786B (en) | Evaluation method for optimal operation of cascade hydropower station | |
| CN114165382A (en) | Method and system for testing absolute efficiency of hydroelectric generating set | |
| CN208621312U (en) | A kind of pressure regulator Performance Test System | |
| CN117539168B (en) | Hydraulic turbine cavitation diagnosis system and method based on semi-physical simulation | |
| CN104458307B (en) | A kind of power canal development hydraulic turbine energy response tests method for controlling water level |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |