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 PDF

Info

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
Application number
CN201711121497.2A
Other languages
Chinese (zh)
Other versions
CN107844680B (en
Inventor
张军
寇攀高
田海平
肖剑
谢捷敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd, State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201711121497.2A priority Critical patent/CN107844680B/en
Publication of CN107844680A publication Critical patent/CN107844680A/en
Application granted granted Critical
Publication of CN107844680B publication Critical patent/CN107844680B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas 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

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

A kind of Hydropower Unit synthesis oscillation area acquisition methods based on gradation of image superposition
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)

  1. 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. 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. 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. 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. 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. 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. 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. 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. 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.
CN201711121497.2A 2017-11-14 2017-11-14 Hydroelectric generating set comprehensive vibration region acquisition method based on image gray level superposition Active CN107844680B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
陈喜阳 等；: "《基于CWT 灰度矩的水电机组振动征兆提取》", 《电力系统自动化》 *

Cited By (7)

* Cited by examiner, † Cited by third party
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
Carballo et al. A methodology to determine the power performance of wave energy converters at a particular coastal location
KR101411420B1 (en) System and method for controlling wind power generator in wind farm through simulation modeling
González et al. A new and efficient method for optimal design of large offshore wind power plants
Mycek et al. Experimental study of the turbulence intensity effects on marine current turbines behaviour. Part I: One single turbine
Saulnier et al. Wave groupiness and spectral bandwidth as relevant parameters for the performance assessment of wave energy converters
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
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
Amaro et al. Big data in power systems leveraging grid optimization and wave energy integration
CN202599970U (en) Wind speed forecasting device of wind power plant and power forecasting system of wind power plant
Germain et al. Protocols for testing marine current energy converters in controlled conditions. Where are we in 2018
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
Wyman et al. A workflow and estimate for the Economic viability of offshore wind projects
Bilal et al. Economic analysis of wind electricity generation in the northern coast of Senegal
Cottura et al. Dynamic Modeling of an Offshore Floating Wind Turbine for Application in the Mediterranean Sea. Energies 2021, 14, 248
Sørensen et al. Towards the North Sea wind power revolution
Al Shibli et al. Artificial intelligent nonlinear auto-regressive external input neural network modeling, design and control of a sea wave electro-mechanical power generating system

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