CN107844680A - A kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition - Google Patents
A kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition Download PDFInfo
- Publication number
- CN107844680A CN107844680A CN201711121497.2A CN201711121497A CN107844680A CN 107844680 A CN107844680 A CN 107844680A CN 201711121497 A CN201711121497 A CN 201711121497A CN 107844680 A CN107844680 A CN 107844680A
- Authority
- CN
- China
- Prior art keywords
- measuring point
- measured value
- vibration
- throw
- area
- 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.)
- Granted
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 77
- 238000003786 synthesis reaction Methods 0.000 title claims abstract description 77
- 230000010355 oscillation Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 92
- 238000013316 zoning Methods 0.000 claims abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 241000239290 Araneae Species 0.000 claims description 9
- DNHVXYDGZKWYNU-UHFFFAOYSA-N lead;hydrate Chemical compound O.[Pb] DNHVXYDGZKWYNU-UHFFFAOYSA-N 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/45—Analysis of texture based on statistical description of texture using co-occurrence matrix computation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30164—Workpiece; Machine component
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Probability & Statistics with Applications (AREA)
- Marketing (AREA)
- Computing Systems (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Physics (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Hydraulic Turbines (AREA)
Abstract
The invention discloses a kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition, including:Step 1:Gather the measured value of each measuring point in Hydropower Unit;Step 2:Isopleth cloud atlas according to each measuring point in the measured value plot step 1 in step 1;Isopleth cloud atlas is divided into non-vibration area and vibrating area;Step 3:The gray matrix E of each isopleth cloud atlas in obtaining step 2i;Step 4:According to synthesis oscillation zoning partial objectives for, each gray matrix E is setiWeight coefficient wiAnd threshold value b;Step 5:Summation is weighted to gray matrix and obtains synthetical matrix, and color range is carried out to the element in synthetical matrix and adjusts to form Gray matrix E;Step 6:Gray Output matrix in step 5 is obtained into the cloud atlas in target synthesis oscillation area into image.The present invention uses the above method, can automatically generate the cloud atlas in synthesis oscillation area, intuitively reacts synthesis oscillation area, and can meet the acquisition demand in the synthesis oscillation area of different target.
Description
Technical field
The invention belongs to hydrogenerator quality inspection, and in particular to a kind of Hydropower Unit synthesis based on gradation of image superposition is shaken
Dynamic area's acquisition methods.
Background technology
Each hydroelectric power plant enters transformation, enlarging peak period, as new machine progressively puts into operation, the aging year by year of old unit, each Water of Power Plant
Group of motors gradually exposes various problems, and unplanned outage happens occasionally, and wherein unit vibration is to influence the stable fortune of unit safety
An important factor for row, the unit being chronically under nuisance vibration, easily produce fatigue rupture and shorten the unit uptime,
Therefore, unit vibration area is divided, for instructing the operation of unit anti-vibration, ensureing that unit safety stable operation is significant.Mesh
The mode that preceding vibrating area obtains has the following disadvantages:1) synthesis oscillation area obtains the quantitative analysis for lacking reflection measuring point characteristic, keeps away
The operational decisions that shake tend to rely on synthesis oscillation area, when using the vibrating area of single measuring point as Main Analysis foundation, to each survey
Point data exceedes the often unrealistic with regard to carrying out anti-vibration measure of vibrating area, need to possess heightened awareness to each measuring point operation characteristic
Synthetic determination can be carried out to vibrating area;2) due to odjective causes such as time and heads, can only often obtain under some heads, load
Vibration data, it is known that vibrating area under data obtain the vibrating area being difficult under other unknown data obtain provide it is effective enough
Foundation, when further carrying out synthesis oscillation area to the vibrating area of each measuring point and obtaining, the synthesis oscillation area under unknown data
Acquisition will be more difficult.3) lack synthesis oscillation area figure and automatically generate means, the figure in conventional synthesis oscillation area is mostly based on
Each measuring point vibrating area, then using expertise, final vibrating area synthesis Man Graphics are carried out, the synthesis different to target is shaken
Dynamic area is obtained (such as comprehensive throw vibrating area, comprehensive fixed component vibrating area, comprehensive hydraulic turbine vibrating area, comprehensive electric generating machine vibration
Area etc.), Transformation of Unit or longtime running cause the situation that vibrating area changes therewith, expertise need to be relied on repeatedly, waste time and energy,
It is less efficient.
When automatically generating means for shortage synthesis oscillation area figure and causing the synthesis oscillation area acquisitions different to target
Waste time and energy, the problem of efficiency is low, it is necessary to which a kind of Hydropower Unit synthesis oscillation acquisition side based on gradation of image superposition is provided
Method solves this problem.
The content of the invention
It is real it is an object of the invention to provide a kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition
The existing fast and effective automatic need quickly generated for obtaining unit synthesis oscillation area figure, meeting the synthesis oscillation area of different target
Ask, meanwhile, the synthesis oscillation area figure of generation can effectively reflect Power Plant nuisance vibration situation.
The present invention provides a kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition, including following step
Suddenly:
Step 1:Gather the measured value of each measuring point in Hydropower Unit;
Wherein, the measuring point of Hydropower Unit comprise at least above lead X to throw, on lead Y-direction throw, under lead X to throw, under lead Y
To throw, water lead X to throw, water lead Y-direction throw, upper spider radial vibration, lower bearing bracket radial vibration, lower bearing bracket vertical vibration,
Any two after out frame radial vibration, top cover radial vibration, top cover vertical vibration, stator in pressure fluctuation, it is each to survey
The measured value of point is that measuring point data of the corresponding measuring point under known hydroelectric generating set load and head is calculated;
Step 2:Isopleth cloud atlas according to each measuring point in the measured value plot step 1 in step 1;
Wherein, the vibration Division according to the measured value of each measuring point is obtained using load, head as two-dimensional coordinate, and measured value is mesh
Target isopleth cloud atlas;
When the measured value of measuring point exceedes the preset limit of corresponding measuring point, the measured value of the measuring point belongs to vibrating area, not less than
During the preset limit of corresponding measuring point, the measured value of the measuring point belongs to non-vibration area;
The isopleth cloud atlas is divided into non-vibration area corresponding with the vibrating area classification of measured value and vibrating area;
Step 3:The gray matrix E of each isopleth cloud atlas in obtaining step 2i;
Wherein, the gray scale in the non-vibration area in each isopleth cloud atlas in step 2 is arranged to 0, by vibrating area
Gray scale be arranged to 1;
EiRepresent gray matrix corresponding to i-th of measuring point;
Step 4:According to each gray matrix E in synthesis oscillation zoning partial objectives for setting steps 3iWeight coefficient wiAnd
Threshold value b;
Wherein, preset comprehensive vibration Division target and the weight coefficient w of each gray matrixiAnd threshold value b correspondence
Relation, the corresponding relation obtain according to expertise;
0 < b < Σ wi, 0≤i≤n;
Wherein, wiThe weight coefficient of i-th of gray matrix is represented, n is the quantity of gray matrix;
Step 5:The weight coefficient of each gray matrix set according to step 4 is weighted summation to gray matrix and obtained
Synthetical matrix, and color range is carried out to the element in synthetical matrix and adjusts to form Gray matrix E;
Wherein, the element in synthetical matrix weighted sum obtained more than threshold value b is arranged to same gray scale, remaining
Element is arranged to another gray scale and obtains Gray matrix E;
Step 6:Gray Output matrix in step 5 is obtained into the cloud atlas in target synthesis oscillation area into image.
The synthesis oscillation zoning partial objectives for is to obtain the synthesis oscillation area of specific function, and wherein specific function can be
Solve all measuring points to the combined influence of vibration, part measuring point to the combined influence of vibration, fixed component to the combined influence of vibration,
Throw to the combined influence of vibration, the combined influence of hydraulic turbine vibrating area, Generator Vibration area combined influence, but do not limit
In this.Therefore, synthesis oscillation zoning partial objectives for can be any one following synthesis oscillation area, but be also not necessarily limited to this, such as:
The larger measuring point vibration of comprehensive throw vibrating area, comprehensive fixed component vibrating area, the vibrating area comprising specific measuring point, weight coefficient
Area, the vibrating area for containing all measuring points, comprehensive hydraulic turbine vibrating area, magnet synthetic electricity generator vibrating area.In order to realize specific function,
The weight coefficient of gray matrix corresponding to each measuring point is set according to expertise, then is combined to obtain Gray matrix
And then obtain representing the cloud atlas in the target synthesis oscillation area of specific function, the scope to synthesis oscillation area is would know that from cloud atlas,
And then carry out anti-vibration operation according to cloud atlas.
Preferably, the isopleth cloud atlas that each measuring point is drawn in the step 2 comprises the following steps:
Step 21:Build load, head, the three-dimensional system of coordinate of measured value;
Step 22:Measured value according to each measuring point in step 1 builds grid corresponding with load, head coordinate in step 21
Matrix;
Step 23:Measured value according to each measuring point in step 1 is generated to measured value corresponding to grid matrix using quick convex closure
Algorithm, triangulation enter row interpolation;
Step 24:Generated using the measured value for exceeding limits value in step 23 interpolation using load, head as coordinate, measured value is mesh
The isopleth cloud atlas of each measuring point of target;
Wherein, step 22 specifically comprises the following steps:
First, the measured value according to each measuring point in step 1 gets the load of the Hydropower Unit and the maximum of head
And minimum value, and load setting and head section are got respectively according to the maximum and minimum value of load and head;
Then, the load setting is divided into N number of equally spaced interval and by described M equally spaced interval of head interval division
To form M × N grid matrixs;
Wherein, N, M are positive integer.
The measured value of the measuring point gathered using step 1 enters row interpolation to measured value corresponding to grid matrix, then to being obtained after interpolation
Measured value divided, select the measured value generation isopleth cloud atlas more than corresponding measuring point limits value, because measured value exceedes limits value,
Then belong to vibrating area depending on it, otherwise belong to non-vibration area, therefore obtained isopleth cloud atlas medium cloud graph region is shaking of matching
Dynamic area.Load, head, the three-dimensional system of coordinate of measured value of structure are really that load, head are two-dimensional coordinate, and measured value is target ginseng
Number.Because the size of data such as the vibration of unit, throw, hydraulic pressure is different under different load, head, give with machine under load, head
The vibration of group, throw, the data such as hydraulic pressure are substantially close, and change is smaller, and therefore, using load, head as coordinate, measured value is target
Isopleth cloud atlas can intuitively get unit vibration area and non-vibration area.The process of drawing isoline cloud atlas is born known to utilizing
Measured value under lotus, head enters row interpolation to the back end of the grid matrix of unknown measured value can improve the reliable of isopleth cloud atlas
Property.
Preferably, the process for obtaining the measured value of measuring point is as follows:
The amplitude of the measuring point data in each cycle in several continuous sampling periods is obtained, the amplitude is each cycle
The maximum of interior measuring point data and the difference of minimum value;
The average value of measuring point amplitude within several described continuous sampling periods is calculated, the average value is one of measuring point
Measured value.
Preferably, the synthesis oscillation zoning partial objectives for is to obtain the synthesis oscillation area of specific function, is influenceed described specific
The weight coefficient of gray matrix corresponding to the measuring point of function is more than the weight coefficient of gray matrix corresponding to remaining measuring point.
Such as specific function is combined influence of the throw to vibration, then the power of the gray matrix of each measuring point of throw is represented
Weight coefficient is greater than the weight coefficient of the gray matrix of remaining measuring point.
Preferably, gray scale weight coefficient all same corresponding to all measuring points of the specific function, remaining measuring point pair are influenceed
The weight coefficient all same for the gray matrix answered.
Preferably, the measuring point of the Hydropower Unit include above lead X to throw, on lead Y-direction throw, under lead X to throw, under lead
Y-direction throw, water lead X to throw, water lead Y-direction throw, upper spider radial vibration, lower bearing bracket radial vibration, lower bearing bracket vertical vibration,
Pressure fluctuation after out frame radial vibration, top cover radial vibration, top cover vertical vibration, stator.
Preferably, when the synthesis oscillation zoning partial objectives for is comprehensive throw vibrating area, above lead X to throw, on lead Y-direction pendulum
Degree, under lead X to throw, under lead Y-direction throw, water leads the weight system that X leads to throw, water gray matrix corresponding to Y-direction throw measuring point
Number is 1, and the weight coefficient of gray matrix is 0 corresponding to remaining measuring point.
Preferably, when the synthesis oscillation zoning partial objectives for is comprehensive fixed component vibrating area, upper spider radial vibration, under
Frame in radial vibration, lower bearing bracket vertical vibration, out frame radial vibration, top cover radial vibration, top cover vertical vibration measuring point pair
The weight coefficient for the gray matrix answered is 1, and the weight coefficient of gray matrix is 0 corresponding to remaining measuring point.
Preferably, when the synthesis oscillation zoning partial objectives for is the vibrating area of all measuring points, all measuring points in Hydropower Unit
The weight coefficient of corresponding gray matrix is 1.
Beneficial effect
Compared with existing synthesis oscillation area acquisition methods, advantages of the present invention has:It is each in Hydropower Unit according to collecting
The measured value of individual measuring point and the synthesis oscillation area cloud atlas that Hydropower Unit can be got after being handled automatically, synthesis oscillation area cloud atlas
Intuitively react Hydropower Unit nuisance vibration situation;Secondly, the present invention is directed to the synthesis oscillation area of different target, need to only set pair
Answer each gray matrix EiWeight coefficient wiAnd threshold value b, carry out the synthesis oscillation area that gradation of image superposition can obtain target
Cloud atlas, meet that the synthesis oscillation area of different target quickly generates needs;It is in addition, comprehensive during synthesis oscillation area obtained by the present invention
The data of multiple measuring points are closed, the reliability in obtained synthesis oscillation area is higher, more presses close to the actual conditions of unit.
Brief description of the drawings
Fig. 1 is that a kind of Hydropower Unit synthesis oscillation area acquisition based on gradation of image superposition provided in an embodiment of the present invention is shown
It is intended to;
Fig. 2 is the isopleth cloud atlas of single measuring point provided in an embodiment of the present invention.
Embodiment
Below in conjunction with embodiment, the present invention is described further.
As shown in figure 1, the present invention provides a kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition,
Comprise the following steps:
Step 1:Gather the measured value of each measuring point in Hydropower Unit;
Wherein, in the present embodiment Hydropower Unit measuring point include above lead X to throw (micron), on lead Y-direction throw (micron),
Under lead X to throw (micron), under lead Y-direction throw (micron), water leads X and leads Y-direction throw (micron), upper machine to throw (micron), water
Frame radial vibration (micron), lower bearing bracket radial vibration (micron), lower bearing bracket vertical vibration (micron), out frame radial vibration
Pressure fluctuation (KPa) after (micron), top cover radial vibration (micron), top cover vertical vibration (micron), stator.Other feasible realities
Apply in example, measuring point comprises at least above-mentioned two, can also include other kinds of.The measured value for obtaining measuring point is in Hydropower Unit pair
Location settings sensor is answered to gather measured value, sensor is frequently located at Turbines, generator, runner key position, example
As vertical turbine spiral case, draft tube and top cover at, surveyed with the sensor of the patterns such as displacement, speed, current vortex, pressure
Measure and obtain data.
Wherein, the process for obtaining the measured value of measuring point is as follows:
The amplitude of the measuring point data in each cycle in several continuous sampling periods is obtained, the amplitude is each cycle
The maximum of interior measuring point data and the difference of minimum value;And calculate measuring point within several described continuous sampling periods amplitude it is flat
Average, the average value are a measured value of measuring point.
Step 2:Isopleth cloud atlas according to each measuring point in the measured value plot step 1 in step 1;
Wherein, the vibration Division according to the measured value of each measuring point is obtained using load, head as two-dimensional coordinate, and measured value is mesh
Target isopleth cloud atlas.
Wherein, when the measured value of measuring point exceedes the preset limit of corresponding measuring point, the measured value of the measuring point belongs to vibrating area, not
During more than the preset limit for corresponding to measuring point, the measured value of the measuring point belongs to non-vibration area.The limits value of measuring point is according to country
Standard《Hydraulic turbine basic technology condition》With《Hydraulic turbine generating fundamental technical specifications》And the safe and stable operation requirement of unit
Set.The isopleth cloud atlas is divided into non-vibration area corresponding with the vibrating area classification of measured value and vibrating area, such as Fig. 2 institutes
Show, be the isopleth cloud atlas of single measuring point, it has intuitively reacted the vibrating area of measuring point and non-vibration area, it is illustrated that grey area
For vibrating area, clear area is non-vibration area.
As shown in figure 1, XiThe measured value of i-th of measuring point is represented, wherein n measuring point altogether, then correspond to and obtain n width isopleth clouds
Scheme, Hydropower Unit includes 13 groups of measuring points in the present embodiment, then corresponds to and obtain 13 secondary isopleth cloud atlas.
Step 3:The gray matrix E of each isopleth cloud atlas in obtaining step 2i;
Wherein, the gray scale in the non-vibration area in each isopleth cloud atlas in step 2 is arranged to 0, by vibrating area
Gray scale be arranged to 1, and then form the gray matrix E of imagei。
Step 4:According to each gray matrix E in synthesis oscillation zoning partial objectives for setting steps 3iWeight coefficient wiAnd
Threshold value b;
Wherein, 0 < b < Σ wi, 0≤i≤n;
Wherein, preset comprehensive vibration Division target and the weight coefficient w of each gray matrixiAnd threshold value b correspondence
Relation, the corresponding relation obtain according to expertise.Synthesis oscillation zoning partial objectives for is to obtain the synthesis of specific function
Vibrating area, wherein specific function can understand all measuring points to the combined influence of vibration, part measuring point to the synthesis shadow of vibration
Ring, fixed component to the combined influence of vibration, throw to the combined influence of vibration, hydraulic turbine vibrating area combined influence, generate electricity
The combined influence in machine vibration area, but it is not restricted to this.Therefore, synthesis oscillation zoning partial objectives for can be that any one is comprehensive as follows
Vibrating area is closed, but is also not necessarily limited to this, such as:Comprehensive throw vibrating area, comprehensive fixed component vibrating area, include specific measuring point
The larger measuring point vibrating area of vibrating area, weight coefficient, the vibrating area for containing all measuring points, comprehensive hydraulic turbine vibrating area, synthesis
Generator Vibration area.Specifically, to obtain different synthesis oscillation areas, then change the weight coefficient and threshold of each gray matrix
Value.
Because corresponding relation obtains according to expertise, therefore the weight coefficient and threshold value according to setting obtain most
The reliability of termination fruit is stronger, and, synthesis oscillation zoning partial objectives for is to obtain the synthesis oscillation of specific function in the present embodiment preferably
Area, influence the measuring point of the specific function gray matrix weight coefficient be more than remaining measuring point gray matrix weight system
Number, influence the gray scale weight coefficient all same of all measuring points of the specific function, the weight system of the gray matrix of remaining measuring point
Number all same.
Step 5:The weight coefficient of each gray matrix set according to step 4 is weighted summation to gray matrix and obtained
Synthetical matrix, and color range is carried out to the element in synthetical matrix and adjusts to form Gray matrix E;
Wherein, summation is weighted to gray matrix according to equation below and obtains synthetical matrix E`, E`=Σ wi·Ei, then
Element in synthetical matrix E` that weighted sum obtains more than threshold value b is arranged to same gray scale, such as 50, surplus element
It is arranged to another gray scale and obtains Gray matrix E.The span of gray scale is 0 to 255, what surplus element was set
Another gray scale is preferably comparative strong color range, such as 255.
Step 6:Gray Output matrix in step 5 is obtained into the cloud atlas in target synthesis oscillation area into image.
As shown in figure 1, the cloud atlas that combination 1 or combination 2 or combination 3 or combination 4 obtain is the cloud atlas in target synthesis oscillation area,
Combination 1, combination 2, combination 3 and combination 4 have corresponded to the cloud atlas in the synthesis oscillation area of different target respectively.
In the present embodiment, combination 1 is that synthesis oscillation zoning partial objectives for is to obtain integrating throw vibrating area, corresponding each measuring point
The weight coefficient of corresponding gray matrix and being set as threshold value:On lead X to throw, on lead Y-direction throw, under lead X to throw,
Under lead Y-direction throw, to lead the weight coefficient that X leads gray matrix corresponding to Y-direction throw measuring point to throw, water be 1 to water, remaining measuring point
The weight coefficient of corresponding gray matrix is 0, and threshold value b is 4.
Combination 2 is that synthesis oscillation zoning partial objectives for is comprehensive fixed component vibrating area, corresponds to gray scale corresponding to each measuring point
The weight coefficient of matrix and being set as threshold value:Upper spider radial vibration, lower bearing bracket radial vibration, lower bearing bracket vertical vibration,
Out frame radial vibration, top cover radial vibration, the weight coefficient of gray matrix corresponding to top cover vertical vibration measuring point are 1,
The part measuring point is fixed component measuring point, and the weight coefficient of gray matrix is 0 corresponding to remaining measuring point, and threshold value b is 2.
Combination 3 is that synthesis oscillation zoning partial objectives for is corresponding for the larger measuring point vibrating area of weight coefficient, corresponding each measuring point
The weight coefficient of gray matrix and being set as threshold value:On lead X to throw, on lead Y-direction throw, under lead X to throw, under lead
Y-direction throw, water lead X to throw, water lead Y-direction throw, the weight coefficient of out frame radial vibration is 0.5, belong to waterpower exciting
The upper spider radial vibration in region, lower bearing bracket radial vibration, lower bearing bracket vertical vibration, top cover radial vibration, top cover vertical vibration,
The weight coefficient of gray matrix corresponding to pressure fluctuation is 1 after stator, and threshold value b is 3.
Combination 4 is the vibrating area that synthesis oscillation zoning partial objectives for is all measuring points, Gray Moment corresponding to corresponding each measuring point
The weight coefficient of battle array and being set as threshold value:On lead X to throw, on lead Y-direction throw, under lead X to throw, under lead Y-direction throw,
Water leads X and leads Y-direction throw, upper spider radial vibration, lower bearing bracket radial vibration, lower bearing bracket vertical vibration, out frame to throw, water
The weight coefficient of gray matrix corresponding to pressure fluctuation is after radial vibration, top cover radial vibration, top cover vertical vibration, stator
1, threshold value b are 0.5.
It should be appreciated that combinations thereof 1- combinations 4 are that the present invention will be described with specific embodiment, combination side of the invention
Formula not limited to this, if such as want to stress to understand to the combined influence of vibration at a certain or certain two measuring point, can be by one or two
The weight of gray matrix corresponding to individual measuring point is arranged to 1, and weight coefficient corresponding to remaining measuring point is arranged to 0.5.
In addition, the isopleth cloud atlas that each measuring point is drawn in above-mentioned steps 2 comprises the following steps:
Step 21:Build load, head, the three-dimensional system of coordinate of measured value;
Step 22:Measured value according to each measuring point in step 1 builds grid corresponding with load, head coordinate in step 21
Matrix;
Step 23:Measured value according to each measuring point in step 1 is generated to measured value corresponding to grid matrix using quick convex closure
Algorithm, triangulation enter row interpolation;
Step 24:Generated using the measured value for exceeding limits value in step 23 interpolation using load, head as coordinate, measured value is mesh
The isopleth cloud atlas of each measuring point of target;
Wherein, step 22 specifically comprises the following steps:
First, the measured value according to each measuring point in step 1 gets the load of the Hydropower Unit and the maximum of head
And minimum value, and load setting and head section are got respectively according to the maximum and minimum value of load and head;
Then, the load setting is divided into N number of equally spaced interval and by described M equally spaced interval of head interval division
To form M × N grid matrixs;
Wherein, N, M are positive integer.Such as M=100, N=100.
A kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition provided by the invention, can be automatic
Synthesis oscillation area cloud atlas is generated, intuitively reacts synthesis oscillation area, further, it is also possible to according to the actual requirements, it is different for target
Synthesis oscillation area the gray matrix of measuring point is set corresponding to weight coefficient and threshold value, and then obtain target synthesis oscillation area
Cloud atlas, meet various demands, have a wide range of application.
It is emphasized that example of the present invention is illustrative, rather than it is limited, therefore the present invention is unlimited
Example described in embodiment, it is every drawn by those skilled in the art's technique according to the invention scheme other
Embodiment, present inventive concept and scope are not departed from, whether modification or replace, also belong to protection model of the invention
Enclose.
Claims (9)
- A kind of 1. Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition, it is characterised in that:Including following step Suddenly:Step 1:Gather the measured value of each measuring point in Hydropower Unit;Wherein, the measuring point of Hydropower Unit comprise at least above lead X to throw, on lead Y-direction throw, under lead X to throw, under lead Y-direction pendulum Degree, water lead X and lead Y-direction throw, upper spider radial vibration, lower bearing bracket radial vibration, lower bearing bracket vertical vibration, stator to throw, water Any two after support radial vibration, top cover radial vibration, top cover vertical vibration, stator in pressure fluctuation, each measuring point Measured value is that measuring point data of the corresponding measuring point under known hydroelectric generating set load and head is calculated;Step 2:Isopleth cloud atlas according to each measuring point in the measured value plot step 1 in step 1;Wherein, the vibration Division according to the measured value of each measuring point is obtained using load, head as two-dimensional coordinate, and measured value is target Isopleth cloud atlas;When the measured value of measuring point exceedes the preset limit of corresponding measuring point, the measured value of the measuring point belongs to vibrating area, not less than corresponding During the preset limit of measuring point, the measured value of the measuring point belongs to non-vibration area;The isopleth cloud atlas is divided into non-vibration area corresponding with the vibrating area classification of measured value and vibrating area;Step 3:The gray matrix E of each isopleth cloud atlas in obtaining step 2i;Wherein, the gray scale in the non-vibration area in each isopleth cloud atlas in step 2 is arranged to 0, by the ash of vibrating area Degree color range is arranged to 1;EiRepresent gray matrix corresponding to i-th of measuring point;Step 4:According to each gray matrix E in synthesis oscillation zoning partial objectives for setting steps 3iWeight coefficient wiAnd threshold value b;Wherein, preset comprehensive vibration Division target and the weight coefficient w of each gray matrixiAnd threshold value b corresponding relation, The corresponding relation obtains according to expertise;0 < b < Σ wi, 0≤i≤n;Wherein, wiThe weight coefficient of i-th of gray matrix is represented, n is the quantity of gray matrix;Step 5:The weight coefficient of each gray matrix set according to step 4 is weighted summation to gray matrix and integrated Matrix, and color range is carried out to the element in synthetical matrix and adjusts to form Gray matrix E;Wherein, the element in synthetical matrix weighted sum obtained more than threshold value b is arranged to same gray scale, surplus element It is arranged to another gray scale and obtains Gray matrix E;Step 6:Gray Output matrix in step 5 is obtained into the cloud atlas in target synthesis oscillation area into image.
- 2. according to the method for claim 1, it is characterised in that:The isopleth cloud atlas of each measuring point is drawn in the step 2 Comprise the following steps:Step 21:Build load, head, the three-dimensional system of coordinate of measured value;Step 22:Measured value according to each measuring point in step 1 builds grid square corresponding with load, head coordinate in step 21 Battle array;Step 23:According to each measuring point in step 1 measured value to measured value corresponding to grid matrix using quick convex closure generating algorithm, Triangulation enters row interpolation;Step 24:Using the measured value generation for the preset limit for exceeding corresponding measuring point in step 23 interpolation using load, head as seat Mark, measured value are the isopleth cloud atlas of each measuring point of target;Wherein, step 22 specifically comprises the following steps:First, the measured value according to each measuring point in step 1 gets the load of the Hydropower Unit and the maximum of head and most Small value, and load setting and head section are got respectively according to the maximum and minimum value of load and head;Then, the load setting is divided into N number of equally spaced interval and by described M equally spaced interval of head interval division with shape Into M × N grid matrixs;Wherein, N, M are positive integer.
- 3. according to the method for claim 1, it is characterised in that:The process for obtaining the measured value of measuring point is as follows:The amplitude of the measuring point data in each cycle in several continuous sampling periods is obtained, the amplitude is to be surveyed in each cycle The maximum of point data and the difference of minimum value;The average value of measuring point amplitude within several described continuous sampling periods is calculated, the average value is a survey of measuring point Value.
- 4. according to the method for claim 1, it is characterised in that:The synthesis oscillation zoning partial objectives for is to obtain specific function Synthesis oscillation area, influence gray matrix corresponding to the measuring point of the specific function weight coefficient be more than remaining measuring point corresponding to The weight coefficient of gray matrix.
- 5. according to the method for claim 4, it is characterised in that:Influence gray scale corresponding to all measuring points of the specific function Weight coefficient all same, the weight coefficient all same of gray matrix corresponding to remaining measuring point.
- 6. according to any described methods of claim 1-5, it is characterised in that:The measuring point of the Hydropower Unit include above lead X to Throw, on lead Y-direction throw, under lead X to throw, under lead Y-direction throw, water lead X to throw, water lead Y-direction throw, upper spider radially shake Dynamic, lower bearing bracket radial vibration, lower bearing bracket vertical vibration, out frame radial vibration, top cover radial vibration, top cover vertical vibration, Pressure fluctuation after stator.
- 7. according to the method for claim 6, it is characterised in that:The synthesis oscillation zoning partial objectives for is that comprehensive throw is vibrated Qu Shi, above lead X to throw, on lead Y-direction throw, under lead X to throw, under lead Y-direction throw, water leads X and leads Y-direction throw to throw, water The weight coefficient of gray matrix corresponding to measuring point is 1, and the weight coefficient of gray matrix is 0 corresponding to remaining measuring point.
- 8. the method according to right wants 6, it is characterised in that:The synthesis oscillation zoning partial objectives for is that comprehensive fixed component is shaken During dynamic area, upper spider radial vibration, lower bearing bracket radial vibration, lower bearing bracket vertical vibration, out frame radial vibration, top cover are radially The weight coefficient of gray matrix is 1 corresponding to vibration, top cover vertical vibration measuring point, the power of gray matrix corresponding to remaining measuring point Weight coefficient is 0.
- 9. according to the method for claim 6, it is characterised in that:The synthesis oscillation zoning partial objectives for is shaken for all measuring points During dynamic area, the weight coefficient of gray matrix corresponding to all measuring points is 1 in Hydropower Unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711121497.2A CN107844680B (en) | 2017-11-14 | 2017-11-14 | Hydroelectric generating set comprehensive vibration region acquisition method based on image gray level superposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711121497.2A CN107844680B (en) | 2017-11-14 | 2017-11-14 | Hydroelectric generating set comprehensive vibration region acquisition method based on image gray level superposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107844680A true CN107844680A (en) | 2018-03-27 |
CN107844680B CN107844680B (en) | 2020-04-07 |
Family
ID=61678780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711121497.2A Active CN107844680B (en) | 2017-11-14 | 2017-11-14 | Hydroelectric generating set comprehensive vibration region acquisition method based on image gray level superposition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107844680B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109057273A (en) * | 2018-09-11 | 2018-12-21 | 洛阳理工学院 | A kind of automatic Mei Feng robot and its application method |
CN109632081A (en) * | 2018-11-23 | 2019-04-16 | 积成电子股份有限公司 | Vibration of wind generating set feature extraction and otherness sentence method for distinguishing |
CN111092442A (en) * | 2019-12-19 | 2020-05-01 | 国网浙江省电力有限公司紧水滩水力发电厂 | Hydroelectric generating set multi-dimensional vibration region fine division method based on decision tree model |
CN112131958A (en) * | 2020-08-28 | 2020-12-25 | 成都信息工程大学 | Method for automatically identifying southwest low vortex |
CN113505509A (en) * | 2021-07-08 | 2021-10-15 | 河北工业大学 | High-precision motor magnetic field prediction method based on improved U-net |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103500425A (en) * | 2013-10-12 | 2014-01-08 | 国家电网公司 | Method for intelligently evaluating running state of hydroelectric generating set |
CN203414278U (en) * | 2013-08-23 | 2014-01-29 | 中国水利水电科学研究院 | System for detecting abnormal condition of vibration of hydroelectric generating set in real time |
CN104463435A (en) * | 2014-11-20 | 2015-03-25 | 云南电网公司电力科学研究院 | Online analysis and evaluation method for stable running areas of hydroelectric generating set |
CN106022965A (en) * | 2016-05-30 | 2016-10-12 | 云南电力试验研究院(集团)有限公司 | Method for establishing power generating vibrating area of hydroelectric generating set |
CN107218997A (en) * | 2017-05-16 | 2017-09-29 | 南京南瑞继保电气有限公司 | A kind of Hydropower Unit vibration monitoring channel abnormal detection method based on operating mode's switch |
-
2017
- 2017-11-14 CN CN201711121497.2A patent/CN107844680B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN203414278U (en) * | 2013-08-23 | 2014-01-29 | 中国水利水电科学研究院 | System for detecting abnormal condition of vibration of hydroelectric generating set in real time |
CN103500425A (en) * | 2013-10-12 | 2014-01-08 | 国家电网公司 | Method for intelligently evaluating running state of hydroelectric generating set |
CN104463435A (en) * | 2014-11-20 | 2015-03-25 | 云南电网公司电力科学研究院 | Online analysis and evaluation method for stable running areas of hydroelectric generating set |
CN106022965A (en) * | 2016-05-30 | 2016-10-12 | 云南电力试验研究院(集团)有限公司 | Method for establishing power generating vibrating area of hydroelectric generating set |
CN107218997A (en) * | 2017-05-16 | 2017-09-29 | 南京南瑞继保电气有限公司 | A kind of Hydropower Unit vibration monitoring channel abnormal detection method based on operating mode's switch |
Non-Patent Citations (1)
Title |
---|
陈喜阳 等;: "《基于CWT 灰度矩的水电机组振动征兆提取》", 《电力***自动化》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109057273A (en) * | 2018-09-11 | 2018-12-21 | 洛阳理工学院 | A kind of automatic Mei Feng robot and its application method |
CN109057273B (en) * | 2018-09-11 | 2024-04-02 | 洛阳理工学院 | Automatic seam beautifying robot and using method thereof |
CN109632081A (en) * | 2018-11-23 | 2019-04-16 | 积成电子股份有限公司 | Vibration of wind generating set feature extraction and otherness sentence method for distinguishing |
CN111092442A (en) * | 2019-12-19 | 2020-05-01 | 国网浙江省电力有限公司紧水滩水力发电厂 | Hydroelectric generating set multi-dimensional vibration region fine division method based on decision tree model |
CN112131958A (en) * | 2020-08-28 | 2020-12-25 | 成都信息工程大学 | Method for automatically identifying southwest low vortex |
CN112131958B (en) * | 2020-08-28 | 2022-05-20 | 成都信息工程大学 | Method for automatically identifying southwest low vortex |
CN113505509A (en) * | 2021-07-08 | 2021-10-15 | 河北工业大学 | High-precision motor magnetic field prediction method based on improved U-net |
Also Published As
Publication number | Publication date |
---|---|
CN107844680B (en) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107844680A (en) | A kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition | |
KR101411420B1 (en) | System and method for controlling wind power generator in wind farm through simulation modeling | |
Carballo et al. | A methodology to determine the power performance of wave energy converters at a particular coastal location | |
Mycek et al. | Experimental study of the turbulence intensity effects on marine current turbines behaviour. Part I: One single turbine | |
González et al. | A new and efficient method for optimal design of large offshore wind power plants | |
Betti et al. | Development of a control-oriented model of floating wind turbines | |
Ovando et al. | Emulation of a low power wind turbine with a DC motor in Matlab/Simulink | |
Taner | Economic analysis of a wind power plant: A case study for the Cappadocia region | |
Lande-Sudall et al. | Co-located offshore wind and tidal stream turbines: Assessment of energy yield and loading | |
CN115758942A (en) | Offshore wind turbine pile foundation load calculation method and system based on coupling effect of atmosphere-sea wave-ocean-structure | |
Rasool et al. | Quantifying the reduction in power variability of co-located offshore wind-wave farms | |
Ahwide et al. | Estimation of electricity generation in libya using processing technology of wind available data: The case study in derna | |
Messineo et al. | Evaluating the performances of small wind turbines: A case study in the south of Italy | |
Pecher | Experimental testing and evaluation of WECs | |
Mishra et al. | Design of non-linear controller for ocean wave energy plant | |
CN202599970U (en) | Wind speed forecasting device of wind power plant and power forecasting system of wind power plant | |
Folley et al. | The adequacy of phase-averaged models for modelling wave farms | |
Commin et al. | Close-proximity tidal phasing for ‘firm’electricity supply | |
Blanco et al. | Development of a laboratory test bench for the emulation of wave energy converters | |
KR101400338B1 (en) | System and method for controlling wind power generator in wind farm | |
Sørensen et al. | Towards the North Sea wind power revolution | |
Saporu et al. | Wind Energy Potential of Maiduguri, Borno State, Nigeria | |
Khalid et al. | Reliability assessment of tidal stream energy: significance for large-scale deployment in the UK | |
Cottura et al. | Dynamic Modeling of an Offshore Floating Wind Turbine for Application in the Mediterranean Sea. Energies 2021, 14, 248 | |
Bilal et al. | Economic analysis of wind electricity generation in the northern coast of Senegal |
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 |