Transmission line galloping judgment and early warning method based on grey correlation model

Abstract

The present invention relates to a power transmission line galloping criterion and early warning method based on a grey correlation model, comprising the following steps: step 1: constructing a power transmission line galloping monitoring system based on an ultrasonic transducer to achieve more accurate and more relaxed monitoring of the dancing trajectory of the conductor; step 2: establishing a power transmission line galloping online monitoring and early warning criterion based on the monitored dancing trajectory; step 3: establishing an index system according to the main causes of transmission line galloping failures; step 4: using a hierarchical analysis method to assign weights to each index according to the influence of each index; step 5: using a grey correlation model to perform correlation analysis on current data and reference data, thereby calculating the risk index of the current line galloping failure and predicting the possibility of the failure; the present invention has the advantages of establishing an index system, constructing monitoring and early warning criterion, effectively preventing galloping failures and reducing the harm of galloping.

Description

Technical Field

[0001] The invention belongs to the technical field of power transmission line galloping warning, and in particular relates to a power transmission line galloping criterion warning method based on a grey correlation model. Background Art

[0002] At present, the online monitoring of power transmission line galloping mainly adopts three methods based on acceleration sensor technology, fiber optic sensor technology and image processing technology. The online monitoring technology of power transmission line galloping based on acceleration sensor mainly uses the acceleration sensor installed on the line to measure the line galloping signal. The measurement results are comprehensively processed by the control center and the corresponding processing results are issued. However, the installation and maintenance are very inconvenient, and the service life is limited. The maintenance cost is high and it is not suitable for large-scale application and promotion; fiber optic sensor technology uses multiple fiber optic sensors evenly distributed on the transmission line to form a quasi-distributed fiber optic sensor network. The load changes of the transmission line are transmitted to the fiber optic Bragg grating through the metal plate, and the collected stress and temperature information are sent to the upper computer control center, which is controlled by the control center. The information is processed comprehensively, but the measurement is susceptible to temperature interference, sensor layout and signal transmission. Image processing technology is used to install a camera near the line to obtain image information of line dancing in real time or at a fixed time, and use a wireless network to transmit the image information to the upper computer control center. The control center analyzes and processes the line dancing image information to obtain the warning level. However, it is difficult to make intelligent judgments on the dancing of the line with ordinary visible light, and it is also difficult to clearly and accurately monitor the dancing of the line at night and in rainy and foggy weather. Therefore, it is very necessary to provide a transmission line dancing criterion warning method based on a grey correlation model to establish an indicator system, construct monitoring and early warning criteria, effectively prevent dancing failures, and reduce the hazards of dancing. Summary of the invention

[0003] The purpose of the present invention is to overcome the shortcomings of the prior art and to provide a transmission line galloping criterion warning method based on a grey correlation model, which can establish an indicator system, construct monitoring and early warning criteria, effectively prevent galloping failures, and reduce the hazards of galloping.

[0004] The object of the present invention is achieved by: a power transmission line galloping criterion early warning method based on a grey correlation model, the method comprising the following steps:

[0005] Step 1: Build a transmission line galloping monitoring system based on ultrasonic transducers to achieve more accurate and more relaxed monitoring of the conductor's galloping trajectory;

[0006] Step 2: Based on the monitored dancing trajectory, establish online monitoring and early warning criteria for transmission line dancing, including line tension criteria for transmission line dancing and phase distance criteria based on the relationship between the dancing amplitude and the phase distance;

[0007] Step 3: Establish an indicator system based on the main causes of transmission line galloping failures;

[0008] Step 4: Use the analytic hierarchy process to assign weights to each indicator according to its influence;

[0009] Step 5: Based on this, the grey correlation model is used to analyze the correlation between the current data and the reference data, so as to calculate the risk index of the current line galloping fault and predict the possibility of the fault.

[0010] The power transmission line dancing monitoring system based on ultrasonic transducer constructed in step 1 includes: ultrasonic transducer transmitting system, ultrasonic transducer receiving system, signal control system, signal receiving system and power supply system; wherein the ultrasonic transducer transmitting system includes: ultrasonic transducer 1 and driving board installed on the sub-spacer rod of the transmission line conductor; the ultrasonic transducer receiving system includes: ultrasonic transducer 2, ultrasonic transducer 3 and signal amplification circuit respectively fixed under the line to be tested; the signal control system includes: main control board and signal isolation fence; the signal acquisition system includes: data acquisition card and host computer.

[0011] The construction of the power transmission line dancing monitoring system based on ultrasonic transducer in step 1 is specifically as follows: when the ultrasonic transducer installed on the transmission line sub-spacer is working, the mechanical wave it emits at time t0 is received by ultrasonic transducer 2 and ultrasonic transducer 3 at time t1 and time t2 respectively. At this time, the distance l2 between ultrasonic transducer 1 and ultrasonic transducer 2 and the distance l3 between ultrasonic transducer 1 and ultrasonic transducer 3 can be respectively expressed as follows: l2=v(t1-t0), l3=v(t2-t0), where v is the speed of sound wave transmission; after considering the influence of wind speed and temperature on the speed of sound wave transmission respectively, the distances l2 and l3 between ultrasonic transducer 1 and ultrasonic transducer 2 and ultrasonic transducer 3 at this moment can be obtained. When the distance l1 between ultrasonic transducer 2 and ultrasonic transducer 3 has been measured, the position (X, Y) of the power transmission line and ultrasonic transducer 1 at this moment can be expressed as follows: P=(l1+l2+l3) / 2, Y=2S / l1, Wherein, S is the area of ​​the triangle enclosed by ultrasonic transducers 1, 2 and 3; P is the intermediate variable for calculating the area of ​​the triangle using Heron's formula; after using the above formula to obtain the dancing state of the outdoor observation line at a certain moment, let the ultrasonic transducer 1 work in time-sharing mode, and emit mechanical waves again after a specific time interval, so as to obtain a series of discrete dancing displacement coordinates (X, Y) of the transmission line. After processing this displacement using the spline interpolation method, the dancing trajectory of the line can be obtained.

[0012] The line tension criterion of the transmission line galloping in step 2 specifically includes the following steps:

[0013] Step 2.1: Calculation of line galloping tension;

[0014] Step 2.2: Establishment of tension criterion: At any time, the tension on point B is the largest, so the suspension point B is selected as the establishment point of the transmission line dancing tension criterion, and the maximum tension that the overhead transmission line can withstand is set as F max When the line dances, the horizontal tension at point B is T'0, and the vertical tension is T' By To ensure that the line is not damaged, it is required to: Where T′0=T0+ΔT0; T′ By =T By +ΔT By ,but To ensure that there is no disconnection or hardware damage during the transmission line galloping, the maximum value of the two items on the left side of formula (12) must be smaller than the item on the right side, that is, Therefore, the relationship between the excessive line tension caused by the transmission line dancing is: Considering that in actual measurement, the ice coating of the conductor has a certain influence on the calculation of the line tension, and it is difficult to measure accurately. At the same time, due to the derivation of the line tension calculation formula, the influence of torsional galloping on the line is ignored. Therefore, the reliability coefficient K1 is introduced in formula (13). In the actual system, the measurement error of the system, the torsional factor of the conductor and the safety margin condition are mainly considered. The line tension hazard criterion when the transmission line gallops is obtained by comprehensive determination of these conditions: (T0+ΔT 0_max ) 2 +(T By +ΔT By_max ) 2 >(K1F max ) 2 (14).

[0015] The calculation of the line tension in step 2.1 is as follows: the length of the transmission line at rest before the dancing is set to L0; the length of the line after the dancing is set to L; E is the comprehensive elastic modulus of the line; A is the cross-sectional area of ​​the line at rest under normal circumstances; ΔT is the change in line tension before and after the dancing; according to Hooke's law: The shape of the overhead transmission line when it is stationary, with the span set to l and the height difference of the line suspension points to h AB ; The height difference angle is β; The tension at the line suspension point is T A , T B ; The horizontal tension of the conductor at rest is T0; The vertical tension component of the conductor suspension point is T Ay , T By; The static displacement of the transmission line before the galloping is y0(x); The weight per unit length of the line when it gallops is W; According to the oblique parabola method, the deflection curve equation of the transmission line before the galloping is: When the transmission line gallops, the galloping displacement equation is approximately: In the formula, A0 is the dancing amplitude; n is the dancing half-wave number; ω is the dancing angular frequency; t is the dancing time; therefore, when the transmission line is dancing, the line deflection curve equation at any time is: By integrating equations (2) and (4) by line length, the length of the transmission line between points A and B can be obtained. The length of the transmission line before galloping is: The length of the transmission line when it is dancing is: Substituting equations (5) and (6) into equation (1), we can obtain that the variation of horizontal tension due to the galloping of the transmission line is: Substitute equations (2) and (4) into equation (7), and let k c =EA / L0, we can get: When the transmission line gallops, the maximum value of the tension change is crucial to the line. The maximum value of the horizontal tension change can be obtained from the above formula: According to the relationship between the horizontal tension and vertical tension of the transmission line, calculate the change of the vertical tension component ΔT at point B By , before the transmission line gallops, the vertical tension component at point B is: When the transmission line gallops, the slope of point B is: Therefore, when the transmission line gallops, the change in the vertical tension component at point B is: Maximum value of vertical tension component ΔT By_max for:

[0016] The specific criterion for establishing the phase distance in step 2 is: set the adjacent two-phase transmission lines a and b, and the minimum distance between the lines in the static state before the galloping occurs is D ab , the safety distance between lines of this voltage level is D s When the transmission line dances, the system detects that the dancing amplitude of line a is A0 and the dancing amplitude of line b is B0. When the line dances, the relationship between the phases flickers due to insufficient phase spacing: D ab -(A0+B0)<D s (15), considering that there are certain errors in actual measurement, the reliability coefficient K2 is introduced in formula (15). In the actual system, the main consideration is the existence of system measurement errors and the guarantee of certain safety margins. Based on these conditions, the criterion for the hazard of the galloping phase spacing of the transmission line is obtained as follows: D ab -(A0+B0)<K2D s .

[0017] The step 4 uses the analytic hierarchy process to assign weights to each indicator according to the degree of influence of each indicator, which specifically includes the following steps:

[0018] Step 4.1: Establishment of structural hierarchical model: For each relevant factor, a multi-level hierarchical structure with mutual influence between each level can be established according to different attributes;

[0019] Step 4.2: Construct a pairwise comparison discriminant matrix: For the established multi-level hierarchical structure, assume that the elements of the next layer dominated by the previous layer criterion C are u1, u2, ..., u n , in order to describe the two elements u under the criterion C i 、u j The relative importance of the two comparison elements needs a scale to measure, so the pairwise comparison discriminant matrix for the criterion C composed of a total of n comparison elements is: J = (u ij ) n×m (16), where u ij For element u i and u j The relative importance of criterion C is proportional to the scale. For the galloping hazard index A as the overall target layer, the indicators of the criterion layer B are natural conditions B1, line parameters B2 and anti-galloping measures B3. The discriminant matrix A can be constructed as follows: For natural condition B1, the indicators of criterion layer C are relative humidity C1, temperature C2, lateral wind speed C3 and precipitation mode C4. Similarly, the judgment matrix B1 can be obtained as follows: For the line parameter B2, the indicators of the criterion layer C are the number of conductor splits C5, the conductor diameter C6 and the span C7. Similarly, the judgment matrix B2 can be obtained as follows: For the anti-dancing measure B3, the indicators of the criterion layer C are the number of rotating spacers C8, the number of double pendulum anti-dancing devices C9 and the number of interphase spacers C 10 , similarly, the judgment matrix B3 is:

[0020] Step 4.3: Calculate the index weight and consistency test: Use the characteristic root method to solve the characteristic root of the judgment matrix A, that is: Aω=λ max ω(17), where λ max is the maximum characteristic root of A; ω is the corresponding eigenvector. After normalizing the obtained ω, the weight under the single criterion can be obtained. At this point, it is necessary to perform a consistency test. When the order of the judgment matrix is ​​greater than 2, it is usually difficult to construct a matrix that satisfies the consistency. However, the judgment matrix should deviate from the consistency condition to a certain degree. Therefore, it is necessary to identify whether the consistency of the judgment matrix is ​​acceptable.

[0021] Step 4.4: Through the above calculation steps, we can get the result of hierarchical single sorting. Find the total ranking result of each indicator, where W (k) is the total order of the k-th level elements on the target; p (k) The order of the elements in the k-1th layer by the elements in the kth layer.

[0022] The consistency check in step 4.3 specifically includes the following steps:

[0023] Step 4.31: Calculate the consistency index CI:

[0024] Step 4.32: The random consistency index RI is related to the order of the judgment matrix. Generally speaking, the larger the matrix order, the greater the possibility of random deviation from consistency.

[0025] Step 4.33: Calculate the consistency ratio CR: CR = CI / RI (19). If CR < 0.1, the consistency of the judgment matrix is ​​considered acceptable; if CR ≥ 0.1, the judgment matrix should be appropriately modified.

[0026] The dancing warning algorithm based on the grey correlation model in step 5 specifically includes the following steps:

[0027] Step 5.1: Grey correlation analysis method: The discrete behavior observation values ​​of system factors are converted into piecewise continuous broken lines by linear interpolation, and then a model for measuring the degree of correlation is constructed based on the geometric characteristics of the broken lines. Suppose the system behavior sequence is: For ρ∈(0,1), let: Then ε(d0,d i ) satisfies the grey relational axiom, where ρ is the discrimination coefficient; ε(d0,d i ) is d0 and d i Grey relational degree; Grey relational degree ε(d0,d i ) is often abbreviated as ε 0i , k-point correlation coefficient ε(d0(k),d i (k)) is abbreviated as ε 0i (k); Grey relational degree R i It is the weighted average value of the grey correlation coefficient: Among them, w k is the weight of indicator k;

[0028] Step 5.2: Dancing warning algorithm: Analyze historical data to obtain three representative reference indicator sets for comparison: danger indicator set, intermediate indicator set and safety indicator set; when calculating, first use meteorological data and ledger information to obtain the 11 indicator values ​​required for warning to form the indicator set to be tested; then combine the weights obtained by the hierarchical analysis method to calculate the grey correlation between the indicator set to be tested and the three reference indicator sets; finally, according to the order of the correlation degree, the dancing warning result is obtained.

[0029] The beneficial effects of the present invention are as follows: the present invention is a transmission line galloping criterion and early warning method based on the grey correlation model. In use, the present invention constructs a transmission line galloping monitoring system based on an ultrasonic transducer. By directly measuring the line galloping displacement, the galloping trajectory of the conductor can be monitored more accurately and in a more relaxed environment, and the basic characteristics of the line galloping can be more realistically characterized. On the basis of the monitoring system, the Hooke's law of elastic materials is used to derive the horizontal tension equation of the galloping line and the vertical tension equation of the high suspension point, and establish the line galloping tension criterion. At the same time, according to the relationship between the galloping amplitude and the phase-to-phase distance, a phase-to-phase distance criterion is established. Furthermore, according to the main cause of the galloping fault, a An indicator system is established, and the hierarchical analysis method is used to assign weights to each indicator according to the degree of influence of each indicator. The three factors of meteorological conditions, line parameters and anti-dancing measures are comprehensively considered. Based on this, the grey correlation model is used to perform correlation analysis on the current data and the reference data, so as to calculate the risk index of the current line dancing failure. The early warning result divides the dancing risk index into 4 levels, which can effectively provide decision-making support for operation and maintenance personnel, effectively realize line dancing early warning, and accurately guide operation and maintenance personnel to formulate inspection strategies, so as to further reduce the hazards caused by the dancing of transmission lines. The present invention has the advantages of establishing an indicator system, constructing monitoring and early warning criteria, effectively preventing dancing failures, and reducing the hazards of dancing. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 This is a schematic diagram of installing the ultrasonic transducer of the present invention on a wire.

[0031] Figure 2 It is a schematic diagram of installing the ultrasonic transducer of the present invention on the span.

[0032] Figure 3 It is a schematic diagram of the force on the oblique parabola of the power transmission line of the present invention.

[0033] Figure 4 It is a schematic diagram of the dancing phase spacing of the transmission line of the present invention.

[0034] Figure 5 Schematic diagram of the index system constructed in the present invention.

[0035] Figure 6It is a flow chart of the present invention. DETAILED DESCRIPTION

[0036] The present invention will be further described below in conjunction with the accompanying drawings.

[0037] Example 1

[0038] like Figure 1-6 As shown, a transmission line galloping judgment and early warning method based on a grey correlation model comprises the following steps:

[0039] Step 1: Build a transmission line galloping monitoring system based on ultrasonic transducers to achieve more accurate and more relaxed monitoring of the conductor's galloping trajectory;

[0040] Step 2: Based on the monitored dancing trajectory, establish online monitoring and early warning criteria for transmission line dancing, including line tension criteria for transmission line dancing and phase distance criteria based on the relationship between the dancing amplitude and the phase distance;

[0041] Step 3: Establish an indicator system based on the main causes of transmission line galloping failures;

[0042] Step 4: Use the analytic hierarchy process to assign weights to each indicator according to its influence;

[0043] Step 5: Based on this, the grey correlation model is used to analyze the correlation between the current data and the reference data, so as to calculate the risk index of the current line galloping fault and predict the possibility of the fault.

[0044] The power transmission line dancing monitoring system based on ultrasonic transducer constructed in step 1 includes: ultrasonic transducer transmitting system, ultrasonic transducer receiving system, signal control system, signal receiving system and power supply system; wherein the ultrasonic transducer transmitting system includes: ultrasonic transducer 1 and driving board installed on the sub-spacer rod of the transmission line conductor; the ultrasonic transducer receiving system includes: ultrasonic transducer 2, ultrasonic transducer 3 and signal amplification circuit respectively fixed under the line to be tested; the signal control system includes: main control board and signal isolation fence; the signal acquisition system includes: data acquisition card and host computer.

[0045] In this embodiment, ① the input end of the ultrasonic transducer 1 is connected to the output end of the driving board; the input end of the driving board is connected to the output end of the main control board of the signal control system, and the output end is connected to the input end of the ultrasonic transducer 1; ② the spatial distances between the ultrasonic transducer 2 and the ultrasonic transducer 3 and the ultrasonic transducer 1 are l2 and l3 respectively, and the distance between the ultrasonic transducer 2 and the ultrasonic transducer 3 is l1; the output ends of the ultrasonic transducer 2 and the ultrasonic transducer 3 are respectively connected to the input end of the signal amplification circuit, and the output end of the signal amplification circuit is connected to the input end of the main control board of the signal control system; ③ the input end of the main control board is connected to the signal amplification circuit of the ultrasonic transducer receiving system. The output end of the main control board is connected to the driving board and the signal isolation fence of the ultrasonic transducer transmitting system, and the power supply end of the main control board is connected to the power supply system; the input end of the signal isolation fence is connected to the output end of the main control board, and the output end of the signal isolation fence is connected to the data acquisition card of the signal acquisition system; ④ the input end of the data acquisition card is connected to the signal isolation fence of the signal control system, and the output end is connected to the host computer system; the input end of the host computer system is connected to the data acquisition card through the PCI-E card slot of the host computer; ⑤ the power supply system is respectively connected to the driving board of the ultrasonic transducer transmitting system and the main control board of the signal control system, and supplies power to the driving board and the main control board respectively.

[0046] The construction of the power transmission line dancing monitoring system based on ultrasonic transducer in step 1 is specifically as follows: The installation diagram of the ultrasonic transducer is as follows Figure 1 and Figure 2 As shown; when the ultrasonic transducer installed on the transmission line sub-spacer is working, the mechanical wave it emits at time t0 is received by ultrasonic transducer 2 and ultrasonic transducer 3 at time t1 and time t2 respectively. At this time, the distance l2 between ultrasonic transducer 1 and ultrasonic transducer 2 and the distance l3 between ultrasonic transducer 1 and ultrasonic transducer 3 can be expressed as follows: l2 = v(t1-t0), l3 = v(t2-t0), where v is the speed of sound wave transmission; after considering the influence of wind speed and temperature on the speed of sound wave transmission respectively, the distances l2 and l3 between ultrasonic transducer 1 and ultrasonic transducer 2 and ultrasonic transducer 3 at this moment can be obtained. When the distance l1 between ultrasonic transducer 2 and ultrasonic transducer 3 has been measured, the position (X, Y) of the transmission line and ultrasonic transducer 1 at this moment can be expressed as follows: P = (l1+l2+l3) / 2, Y=2S / l1, Wherein, S is the area of ​​the triangle enclosed by ultrasonic transducers 1, 2 and 3; P is the intermediate variable for calculating the area of ​​the triangle using Heron's formula; after using the above formula to obtain the dancing state of the outdoor observation line at a certain moment, let the ultrasonic transducer 1 work in time-sharing mode, and emit mechanical waves again after a specific time interval, so as to obtain a series of discrete dancing displacement coordinates (X, Y) of the transmission line. After processing this displacement using the spline interpolation method, the dancing trajectory of the line can be obtained.

[0047] The present invention is a transmission line galloping criterion early warning method based on a grey correlation model. In use, the present invention constructs a transmission line galloping monitoring system based on an ultrasonic transducer. By directly measuring the line galloping displacement, the galloping trajectory of the conductor can be monitored more accurately and in a more relaxed environment, and the basic characteristics of the line galloping can be more realistically characterized. On the basis of the monitoring system, the Hooke's law of elastic materials is used to derive the horizontal tension equation of the galloping line and the vertical tension equation of the high suspension point, and establish the line galloping tension criterion; at the same time, the phase distance criterion is established according to the relationship between the galloping amplitude and the phase distance; furthermore, an indicator is established according to the main cause of the galloping fault. system, and assign weights to each indicator according to the degree of influence of each indicator by using the hierarchical analysis method, comprehensively consider the three factors of meteorological conditions, line parameters and anti-dancing measures, and then use the grey correlation model to perform correlation analysis on the current data and the reference data, so as to calculate the risk index of the current line dancing failure. The early warning result divides the dancing risk index into 4 levels, which can effectively provide decision-making support for operation and maintenance personnel, effectively realize line dancing early warning, and accurately guide operation and maintenance personnel to formulate inspection strategies, so as to further reduce the hazards caused by the dancing of transmission lines. The present invention has the advantages of establishing an indicator system, constructing monitoring and early warning criteria, effectively preventing dancing failures, and reducing the hazards of dancing.

[0048] Example 2

[0049] like Figure 1-6 As shown, a transmission line galloping judgment and early warning method based on a grey correlation model comprises the following steps:

[0050] Step 1: Build a transmission line galloping monitoring system based on ultrasonic transducers to achieve more accurate and more relaxed monitoring of the conductor's galloping trajectory;

[0051] Step 2: Based on the monitored dancing trajectory, establish online monitoring and early warning criteria for transmission line dancing, including line tension criteria for transmission line dancing and phase distance criteria based on the relationship between the dancing amplitude and the phase distance;

[0052] Step 3: Establish an indicator system based on the main causes of transmission line galloping failures;

[0053] Step 4: Use the analytic hierarchy process to assign weights to each indicator according to its influence;

[0054] Step 5: Based on this, the grey correlation model is used to analyze the correlation between the current data and the reference data, so as to calculate the risk index of the current line galloping fault and predict the possibility of the fault.

[0055] The line tension criterion of the transmission line galloping in step 2 specifically includes the following steps:

[0056] Step 2.1: Calculation of line galloping tension;

[0057] Step 2.2: Establishment of tension criterion: At any time, Figure 3 The tension at point B is the largest, so the suspension point B is selected as the establishment point of the transmission line dancing tension criterion, and the maximum tension that the overhead transmission line can withstand is set as F max When the line dances, the horizontal tension at point B is T'0, and the vertical tension is T' By To ensure that the line is not damaged, it is required to: Where T′0=T0+ΔT0; T′ By =T By +ΔT By ,but To ensure that there are no accidents such as line breakage and hardware damage when the transmission line is dancing, the maximum value of the two items on the left side of formula (12) must be smaller than the item on the right side, that is, Therefore, the relationship between the excessive line tension caused by the transmission line dancing is: Considering that in actual measurement, the ice coating of the conductor has a certain influence on the calculation of the line tension, and it is difficult to measure accurately. At the same time, due to the derivation of the line tension calculation formula, the influence of torsional galloping on the line is ignored. Therefore, the reliability coefficient K1 is introduced in formula (13), and the value of K1 is preferably 0.6-0.8. In the actual system, the measurement error of the system, the torsion factor of the conductor and the safety margin are mainly considered. The comprehensive determination of these conditions leads to the conclusion that the line tension hazard criterion when the transmission line gallops is: (T0+ΔT 0_max ) 2 +(T By +ΔT By_max ) 2 >(K1F max ) 2 (14).

[0058] The calculation of line tension in step 2.1 is as follows: Since the frequency of the transmission line dancing is low (0.1-3 Hz), the elastic deformation of the line is slow within one cycle of the dancing. Ignoring the inertial force of any node of the line moving along the length direction of the line, the line tension change caused by the dancing of the transmission line and the line length change caused by the dancing are within the linear elastic range of the material. Therefore, the line tension change can be calculated according to Hooke's law; the length of the transmission line at rest before the dancing is set to L0; the line length after the dancing is set to L; E is the comprehensive elastic modulus of the line; A is the cross-sectional area of ​​the line at rest under normal circumstances; ΔT is the change in line tension before and after the dancing; According to Hooke's law: like Figure 3As shown in Figure 2, the shape of the overhead transmission line when it is stationary ( Figure 3 The solid line part AB), the set span is l; the height difference of the line suspension point is h AB ; The height difference angle is β; The tension at the line suspension point is T A , T B ; The horizontal tension of the conductor at rest is T0; The vertical tension component of the conductor suspension point is T Ay , T By ; The static displacement of the transmission line before the galloping is y0(x); The weight per unit length of the line when it gallops is W; According to the oblique parabola method, the deflection curve equation of the transmission line before the galloping is: When the transmission line gallops, the galloping displacement equation is approximately: In the formula, A0 is the dancing amplitude; n is the dancing half-wave number; ω is the dancing angular frequency; t is the dancing time; therefore, when the transmission line is dancing, the line deflection curve equation at any time is: By integrating equations (2) and (4) by line length, the length of the transmission line between points A and B can be obtained. The length of the transmission line before galloping is: The length of the transmission line when it is dancing is: Substituting equations (5) and (6) into equation (1), we can obtain that the variation of horizontal tension due to the galloping of the transmission line is: Substitute equations (2) and (4) into equation (7), and let k c =EA / L0, we can get: When the transmission line gallops, the maximum value of the tension change is crucial to the line. The maximum value of the horizontal tension change can be obtained from the above formula: According to the relationship between the horizontal tension and vertical tension of the transmission line, calculate Figure 3 The change in the vertical tension component of point B (high suspension point) ΔT By , before the transmission line gallops, the vertical tension component at point B is: When the transmission line gallops, the slope of point B is: Therefore, when the transmission line gallops, the change in the vertical tension component at point B (high suspension point) is: Maximum value of vertical tension component ΔT By_max for:

[0059] The specific criterion for establishing the phase distance in step 2 is: Figure 4 As shown, assume that the adjacent two-phase transmission lines a and b, in the static state before the galloping occurs, the minimum distance between the lines is D ab , the safety distance between lines of this voltage level is D sWhen the transmission line dances, the system detects that the dancing amplitude of line a is A0 and the dancing amplitude of line b is B0. When the line dances, the relationship between the phases flickers due to insufficient phase spacing: D ab -(A0+B0)<D s (15), considering that there are certain errors in actual measurement, the reliability coefficient K2 is introduced in formula (15), and the value of K2 is preferably in the range of 1.2 to 1.3. In the actual system, the main consideration is the measurement error of the system and the guarantee of a certain safety margin. Based on these conditions, the criterion for the hazard of the inter-phase spacing of the transmission line is: D ab -(A0+B0)<K2D s .

[0060] In this embodiment, the establishment of the index system in step 3 is specifically as follows: The important factors causing the galloping are as follows: ① Conductor icing: The formation of icing must meet the three conditions of high relative humidity, suitable temperature and suitable wind speed that can make water droplets in the air move; ② Wind excitation: In addition to icing, the formation of galloping also requires stable laminar wind excitation, which is related to wind speed, wind direction and line axis; ③ Line structure parameters: mainly including the type of transmission line, conductor cross-sectional area, line span, etc.; The galloping of the transmission line depends on natural conditions, basic line parameters and anti-galloping measures. Considering the importance of the indicators and the difficulty of obtaining the indicator values, 10 indicators are selected as reference indicators, and the constructed index system is as follows Figure 5 shown.

[0061] The step 4 uses the analytic hierarchy process to assign weights to each indicator according to the degree of influence of each indicator, which specifically includes the following steps:

[0062] Step 4.1: Establishment of structural hierarchical model: For each relevant factor, a multi-level hierarchical structure with mutual influence between levels can be established according to different attributes. The structural hierarchical model constructed by 10 reference indicators is as follows: Figure 5 As shown;

[0063] Step 4.2: Construct a pairwise comparison discriminant matrix: For the established multi-level hierarchical structure, assume that the elements of the next layer dominated by the previous layer criterion C are u1, u2, ..., u n , in order to describe the two elements u under the criterion C i 、u j The relative importance of , requires a scale to measure, the scale assignment is shown in Table 1, then the total number of comparison elements n for the criterion C constitutes a pairwise comparison discriminant matrix: J = (u ij ) n×m (16), where u ij For element u iand u j The proportional scale of the importance of the relative criterion C is shown in Table 1. For the galloping hazard index A as the overall target layer, the indicators of the criterion layer B are natural conditions B1, line parameters B2 and anti-dancing measures B3. The discriminant matrix A can be constructed as follows: For natural condition B1, the indicators of criterion layer C are relative humidity C1, temperature C2, lateral wind speed C3 and precipitation mode C4. Similarly, the judgment matrix B1 can be obtained as follows: For the line parameter B2, the indicators of the criterion layer C are the number of conductor splits C5, the conductor diameter C6 and the span C7. Similarly, the judgment matrix B2 can be obtained as follows: For the anti-dancing measure B3, the indicators of the criterion layer C are the number of rotating spacers C8, the number of double pendulum anti-dancing devices C9 and the number of interphase spacers C 10 , similarly, the judgment matrix B3 is:

[0064] In this embodiment, Table 1 AHP ​​scale

[0065]

[0066] Step 4.3: Calculate the index weight and consistency test: Use the characteristic root method to solve the characteristic root of the judgment matrix A, that is: Aω=λ max ω(17), where λ max is the maximum characteristic root of A; ω is the corresponding eigenvector. After normalizing the obtained ω, the weight under the single criterion can be obtained. At this point, it is necessary to perform a consistency test. When the order of the judgment matrix is ​​greater than 2, it is usually difficult to construct a matrix that satisfies the consistency. However, the judgment matrix should deviate from the consistency condition to a certain degree. Therefore, it is necessary to identify whether the consistency of the judgment matrix is ​​acceptable.

[0067] Step 4.4: Through the above calculation steps, we can get the result of hierarchical single sorting. The total ranking results of each indicator are shown in Table 3, where W (k) is the total order of the k-th level elements on the target; p (k) The order of the elements in the k-1th layer by the elements in the kth layer.

[0068] In this embodiment, Table 3 Index weights

[0069]

[0070] The consistency check in step 4.3 specifically includes the following steps:

[0071] Step 4.31: Calculate the consistency index CI:

[0072] Step 4.32: The random consistency index RI is related to the order of the judgment matrix. Generally speaking, the larger the matrix order, the greater the possibility of random deviation from the consistency. The corresponding relationship between RI and n is shown in Table 2;

[0073] In this embodiment, Table 2 shows the relationship between the random consistency index and the matrix order.

[0074]

[0075] Step 4.33: Calculate the consistency ratio CR: CR = CI / RI (19). If CR < 0.1, the consistency of the judgment matrix is ​​considered acceptable; if CR ≥ 0.1, the judgment matrix should be appropriately modified.

[0076] The dancing warning algorithm based on the grey correlation model in step 5 specifically includes the following steps:

[0077] Step 5.1: Grey correlation analysis method: The discrete behavior observation values ​​of system factors are converted into piecewise continuous broken lines by linear interpolation, and then a model for measuring the degree of correlation is constructed based on the geometric characteristics of the broken lines. Suppose the system behavior sequence is: For ρ∈(0,1), let: Then ε(d0,d i ) satisfies the grey relational axiom, where ρ is the discrimination coefficient; ε(d0,d i ) is d0 and d i Grey relational degree; Grey relational degree ε(d0,d i ) is often abbreviated as ε 0i , k-point correlation coefficient ε(d0(k),d i (k)) is abbreviated as ε 0i (k); Grey relational degree R i It is the weighted average value of the grey correlation coefficient: Among them, w k is the weight of indicator k;

[0078] Step 5.2: Dancing warning algorithm: The calculation process of the dancing warning algorithm proposed based on the grey correlation theory is as follows: Figure 6 As shown in the figure, the historical data are analyzed to obtain three representative reference indicator sets for comparison: the danger indicator set, the intermediate indicator set and the safety indicator set; when calculating, the 11 indicator values ​​required for the early warning are first obtained through meteorological data and ledger information to form the indicator set to be tested; then, combined with the weights obtained by the hierarchical analysis method, the grey correlation between the indicator set to be tested and the three reference indicator sets is calculated respectively; finally, according to the arrangement order of the correlation degree, the dancing warning result is obtained.

[0079] In this embodiment, the specific process of the dancing warning algorithm is as follows:

[0080] ① Determination of reference index values: Based on the historical data of galloping and the occurrence mechanism of galloping, the three index sets are analyzed as reference index sets. The index of the safety index set is the value of the state with the minimum possibility of galloping, and each index of the danger index set is the value of the state with the maximum possibility of galloping. The index of the intermediate index set is obtained by referring to the critical value of galloping, as shown in Table 4. Since galloping is impossible when the line is not covered with ice, the relative humidity and temperature need to be preprocessed in the calculation process. If the relative humidity is less than 60%, it is calculated as 60%; if the temperature is higher than 10℃, it is calculated as 10℃;

[0081] Table 4 Reference indicator set

[0082]

[0083] Note: For the values ​​of double pendulum anti-dance device, rotary spacer and phase spacer, "1" means installed, "0" means not installed, and "-1" is only for reference and has no practical significance.

[0084] The precipitation modes and corresponding index values ​​during calculation are shown in Table 5.

[0085] Table 5 Precipitation mode index values

[0086]

[0087] ② Calculation of transverse wind speed: Except for transverse wind speed, other required index data can be directly read from the ledger information or meteorological data; the transverse wind speed is related to the wind speed, wind direction, line height and line direction. According to the usual calculation method, the law of average wind speed change along the height is: Among them, U′ and Z are the average wind speed and height above the ground at any point respectively; Z s , U s are the standard height (generally set at 10m) and the average wind speed at the standard height; α is the ground roughness coefficient, the rougher the ground, the greater its value; as shown in Table 6;

[0088] Table 6 Roughness coefficient

[0089]

[0090] Considering that winter winds are generally northerly and that in the collected data, all wind directions during dancing are northerly, we only consider the angle θ between the line direction and the north-south direction. Among them, W x W is the latitude of the tower x; y is the latitude of the tower y; J x is the longitude of the tower x; Jy is the longitude of the tower y; the calculation formula of the north wind lateral wind speed U is: U = U′sinθ(25);

[0091] ③ Normalization of index values: The index set to be measured is: D0 = [d1, d2, d3, d4, d5, d6, d7, d8, d9, d 10 ]; the safety index set is: D1 = [60, 10, 0, 0, 1, 23, 76, 120, 1, 1, 1]; the intermediate index set is: D2 = [70, 5, 4, 0.5, 4, 30.00, 400, 0, 0, 0]; the danger index set is: D3 = [90, -2, 15, 1, 6, 39.90, 600, -1, -1, -1]; D = [D0 D1D2 D3] T ; The evaluation indicators often have different dimensions and orders of magnitude, so they cannot be directly compared. It is necessary to use the extreme value processing method to normalize the matrix D:

[0092] For the air temperature C2, the number of double pendulum anti-dance devices C8, the number of rotary spacers C9 and the number of interphase spacers C 10 For these four indicators, the smaller the value, the greater the risk index, so the minimum value processing method is used: For the six indicators of relative humidity C1, lateral wind speed C3, precipitation mode C4, number of conductor splits C5, conductor diameter C6, and span C7, the larger the value, the greater the risk index, so the maximum value processing method is adopted: Among them, M j is the maximum value of the data taken by index j; m j is the minimum value of the data taken by index j; the resulting normalized matrix is:

[0093] ④ Calculation of correlation coefficient matrix: Let ρ = 0.5. According to the grey correlation analysis method, the correlation coefficient ε between the jth indicator in the indicator set to be tested and the ith typical indicator set can be obtained by formula (6): 0i (j), the correlation matrix of the correlation coefficient is obtained as follows:

[0094] ⑤ Calculation of grey correlation and dancing risk index: According to the correlation coefficient row vector E i and the indicator weight vector W i , the correlation R between the test data and the three typical indicator sets can be obtained i (i=1,2,3) is: The greater the correlation, the closer the index sequence D0 to the reference sequence D iThe closer, the greater the probability of the occurrence of dancing. Therefore, the dancing risk index of the indicator sequence to be tested and the probability of dancing can be obtained according to the order of the correlation degree. The results of the descending order of correlation degree and the corresponding warning results are shown in Table 7.

[0095] Table 7 Early warning results

[0096] Order of association Dancing Hazard Index Dancing Possibilities <![CDATA[R1>R2>R3]]> 0 impossible <![CDATA[R2>R1>R3]]> 1 Almost impossible <![CDATA[R2>R3>R1]]> 2 There is a certain possibility <![CDATA[R3>R2>R1]]> 3 And it is possible

[0097] The present invention is a transmission line galloping criterion early warning method based on a grey correlation model. In use, the present invention constructs a transmission line galloping monitoring system based on an ultrasonic transducer. By directly measuring the line galloping displacement, the galloping trajectory of the conductor can be monitored more accurately and in a more relaxed environment, and the basic characteristics of the line galloping can be more realistically characterized. On the basis of the monitoring system, the Hooke's law of elastic materials is used to derive the horizontal tension equation of the galloping line and the vertical tension equation of the high suspension point, and establish the line galloping tension criterion; at the same time, the phase distance criterion is established according to the relationship between the galloping amplitude and the phase distance; furthermore, an indicator is established according to the main cause of the galloping fault. system, and assign weights to each indicator according to the degree of influence of each indicator by using the hierarchical analysis method, comprehensively consider the three factors of meteorological conditions, line parameters and anti-dancing measures, and then use the grey correlation model to perform correlation analysis on the current data and the reference data, so as to calculate the risk index of the current line dancing failure. The early warning result divides the dancing risk index into 4 levels, which can effectively provide decision-making support for operation and maintenance personnel, effectively realize line dancing early warning, and accurately guide operation and maintenance personnel to formulate inspection strategies, so as to further reduce the hazards caused by the dancing of transmission lines. The present invention has the advantages of establishing an indicator system, constructing monitoring and early warning criteria, effectively preventing dancing failures, and reducing the hazards of dancing.

Claims

1. A transmission line galloping judgment and early warning method based on a grey correlation model, characterized by: The method comprises the following steps: Step 1: Build a transmission line galloping monitoring system based on ultrasonic transducers to achieve more accurate and more relaxed monitoring of the conductor's galloping trajectory; Step 2: Based on the monitored dancing trajectory, establish online monitoring and early warning criteria for transmission line dancing, including line tension criteria for transmission line dancing and phase distance criteria based on the relationship between the dancing amplitude and the phase distance; Step 3: Establish an indicator system based on the main causes of transmission line galloping failures; Step 4: Use the analytic hierarchy process to assign weights to each indicator according to its influence; Step 5: Based on this, the grey correlation model is used to analyze the correlation between the current data and the reference data, so as to calculate the risk index of the current line galloping fault and predict the possibility of the fault.

2. The power transmission line galloping judgment and early warning method based on the grey correlation model as claimed in claim 1 is characterized by: The power transmission line dancing monitoring system based on ultrasonic transducer constructed in step 1 includes: ultrasonic transducer transmitting system, ultrasonic transducer receiving system, signal control system, signal receiving system and power supply system; wherein the ultrasonic transducer transmitting system includes: ultrasonic transducer 1 and driving board installed on the sub-spacer rod of the transmission line conductor; the ultrasonic transducer receiving system includes: ultrasonic transducer 2, ultrasonic transducer 3 and signal amplification circuit respectively fixed under the line to be tested; the signal control system includes: main control board and signal isolation fence; the signal acquisition system includes: data acquisition card and host computer.

3. The power transmission line galloping judgment and early warning method based on the grey correlation model as claimed in claim 2 is characterized by: The construction of the power transmission line dancing monitoring system based on ultrasonic transducer in step 1 is specifically as follows: when the ultrasonic transducer installed on the transmission line sub-spacer is working, the mechanical wave it emits at time t0 is received by ultrasonic transducer 2 and ultrasonic transducer 3 at time t1 and time t2 respectively. At this time, the distance l2 between ultrasonic transducer 1 and ultrasonic transducer 2 and the distance l3 between ultrasonic transducer 1 and ultrasonic transducer 3 can be respectively expressed as follows: l2=v(t1-t0), l3=v(t2-t0), where v is the speed of sound wave transmission; after considering the influence of wind speed and temperature on the speed of sound wave transmission respectively, the distances l2 and l3 between ultrasonic transducer 1 and ultrasonic transducer 2 and ultrasonic transducer 3 at this moment can be obtained. When the distance l1 between ultrasonic transducer 2 and ultrasonic transducer 3 has been measured, the position (X, Y) of the power transmission line and ultrasonic transducer 1 at this moment can be expressed as follows: P=(l1+l2+l3) / 2, Y=2S / l1, Wherein, S is the area of ​​the triangle enclosed by ultrasonic transducers 1, 2 and 3; P is the intermediate variable for calculating the area of ​​the triangle using Heron's formula; after using the above formula to obtain the dancing state of the outdoor observation line at a certain moment, let the ultrasonic transducer 1 work in time-sharing mode, and emit mechanical waves again after a specific time interval, so as to obtain a series of discrete dancing displacement coordinates (X, Y) of the transmission line. After processing this displacement using the spline interpolation method, the dancing trajectory of the line can be obtained.

4. The power transmission line galloping judgment and early warning method based on grey correlation model as claimed in claim 1 is characterized by: The line tension criterion of the transmission line galloping in step 2 specifically includes the following steps: Step 2.1: Calculation of line galloping tension; Step 2.2: Establishment of tension criterion: At any time, the tension on point B is the largest, so the suspension point B is selected as the establishment point of the transmission line dancing tension criterion, and the maximum tension that the overhead transmission line can withstand is set as F max When the line dances, the horizontal tension at point B is T'0, and the vertical tension is T' By To ensure that the line is not damaged, it is required to: Where T′0=T0+ΔT0; T′ By =T By +ΔT By ,but To ensure that there is no disconnection or hardware damage during the transmission line galloping, the maximum value of the two items on the left side of formula (12) must be smaller than the item on the right side, that is, Therefore, the relationship between the excessive line tension caused by the transmission line dancing is: Considering that in actual measurement, the ice coating of the conductor has a certain influence on the calculation of the line tension, and it is difficult to measure accurately. At the same time, due to the derivation of the line tension calculation formula, the influence of torsional galloping on the line is ignored. Therefore, the reliability coefficient K1 is introduced in formula (13). In the actual system, the measurement error of the system, the torsional factor of the conductor and the safety margin condition are mainly considered. The line tension hazard criterion when the transmission line gallops is obtained by comprehensive determination of these conditions: (T0+ΔT 0_max ) 2 +(T By +ΔT By_max ) 2 >(K1F max ) 2 (14).

5. The power transmission line galloping judgment and early warning method based on grey correlation model as claimed in claim 4 is characterized by: The calculation of the line tension in step 2.1 is as follows: the length of the transmission line at rest before the dancing is set to L0; the length of the line after the dancing is set to L; E is the comprehensive elastic modulus of the line; A is the cross-sectional area of ​​the line at rest under normal circumstances; ΔT is the change in line tension before and after the dancing; according to Hooke's law: The shape of the overhead transmission line when it is stationary, with the span set to l and the height difference of the line suspension points to h AB ; The height difference angle is β; The tension at the line suspension point is T A 、T B ; The horizontal tension of the conductor at rest is T0; The vertical tension component of the wire suspension point is T Ay 、T By ; The static displacement of the transmission line before galloping is y0(x); The weight per unit length of the line when it gallops is W; According to the oblique parabola method, the deflection curve equation of the transmission line before galloping is: When the transmission line gallops, the galloping displacement equation is approximately: In the formula, A0 is the dancing amplitude; n is the dancing half-wave number; ω is the dancing angular frequency; t is the dancing time; therefore, when the transmission line is dancing, the line deflection curve equation at any time is: By integrating equations (2) and (4) by line length, the length of the transmission line between points A and B can be obtained. The length of the transmission line before galloping is: The length of the transmission line when it is dancing is: Substituting equations (5) and (6) into equation (1), we can obtain that the variation of horizontal tension due to the galloping of the transmission line is: Substitute equations (2) and (4) into equation (7), and let k c =EA / L0, we can get: When the transmission line gallops, the maximum value of the tension change is crucial to the line. The maximum value of the horizontal tension change can be obtained from the above formula: According to the relationship between the horizontal tension and vertical tension of the transmission line, calculate the change of the vertical tension component ΔT at point B By , before the transmission line gallops, the vertical tension component at point B is: When the transmission line gallops, the slope of point B is: Therefore, when the transmission line gallops, the change in the vertical tension component at point B is: Maximum value of vertical tension component ΔT By_max for:

6. The power transmission line galloping judgment and early warning method based on grey correlation model as claimed in claim 4 is characterized by: The specific criterion for establishing the phase distance in step 2 is: set the adjacent two-phase transmission lines a and b, and the minimum distance between the lines in the static state before the galloping occurs is D ab , the safety distance between lines of this voltage level is D s When the transmission line dances, the system detects that the dancing amplitude of line a is A0 and the dancing amplitude of line b is B0. When the line dances, the relationship between the phases flickers due to insufficient phase spacing: D ab -(A0+B0)<D s (15), considering that there are certain errors in actual measurement, the reliability coefficient K2 is introduced in formula (15). In the actual system, the main consideration is the existence of system measurement errors and the guarantee of certain safety margins. Based on these conditions, the criterion for the hazard of the galloping phase spacing of the transmission line is obtained as follows: D ab -(A0+B0)<K2D s .

7. The power transmission line galloping judgment and early warning method based on grey correlation model as claimed in claim 1 is characterized by: The step 4 uses the analytic hierarchy process to assign weights to each indicator according to the degree of influence of each indicator, which specifically includes the following steps: Step 4.1: Establishment of structural hierarchical model: For each relevant factor, a multi-level hierarchical structure with mutual influence between each level can be established according to different attributes; Step 4.2: Construct a pairwise comparison discriminant matrix: For the established multi-level hierarchical structure, assume that the elements of the next layer dominated by the previous layer criterion C are u1, u2, ..., u n , in order to describe the two elements u under the criterion C i 、u j The relative importance of the two comparison elements needs a scale to measure, so the pairwise comparison discriminant matrix for the criterion C composed of a total of n comparison elements is: J = (u ij ) n×m (16), where u ij For element u i and u j The relative importance of criterion C is proportional to the scale. For the galloping hazard index A as the overall target layer, the indicators of the criterion layer B are natural conditions B1, line parameters B2 and anti-galloping measures B3. The discriminant matrix A can be constructed as follows: For natural condition B1, the indicators of criterion layer C are relative humidity C1, temperature C2, lateral wind speed C3 and precipitation mode C4. Similarly, the judgment matrix B1 can be obtained as follows: For the line parameter B2, the indicators of the criterion layer C are the number of conductor splits C5, the conductor diameter C6 and the span C7. Similarly, the judgment matrix B2 can be obtained as follows: For the anti-dancing measure B3, the indicators of the criterion layer C are the number of rotating spacers C8, the number of double pendulum anti-dancing devices C9 and the number of interphase spacers C 10 , similarly, the judgment matrix B3 is: Step 4.3: Calculate the index weight and consistency test: Use the characteristic root method to solve the characteristic root of the judgment matrix A, that is: Aω=λ max ω(17), where λ max is the maximum characteristic root of A; ω is the corresponding eigenvector. After normalizing the obtained ω, the weight under the single criterion can be obtained. At this point, it is necessary to perform a consistency test. When the order of the judgment matrix is ​​greater than 2, it is usually difficult to construct a matrix that satisfies the consistency. However, the judgment matrix should deviate from the consistency condition to a certain degree. Therefore, it is necessary to identify whether the consistency of the judgment matrix is ​​acceptable. Step 4.4: Through the above calculation steps, we can get the result of hierarchical single sorting. Find the total ranking result of each indicator, where W (k) is the total order of the k-th level elements on the target; p (k) The order of the elements in the k-1th layer by the elements in the kth layer.

8. The power transmission line galloping judgment and early warning method based on grey correlation model as claimed in claim 7 is characterized by: The consistency check in step 4.3 specifically includes the following steps: Step 4.31: Calculate the consistency index CI: Step 4.32: The random consistency index RI is related to the order of the judgment matrix. Generally speaking, the larger the matrix order, the greater the possibility of random deviation from consistency. Step 4.33: Calculate the consistency ratio CR: CR = CI / RI (19). If CR < 0.1, the consistency of the judgment matrix is ​​considered acceptable; if CR ≥ 0.1, the judgment matrix should be appropriately modified.

9. The power transmission line galloping judgment and early warning method based on grey correlation model as claimed in claim 1, characterized in that: The dancing warning algorithm based on the grey correlation model in step 5 specifically includes the following steps: Step 5.1: Grey correlation analysis method: The discrete behavior observation values ​​of system factors are converted into piecewise continuous broken lines by linear interpolation, and then a model for measuring the degree of correlation is constructed based on the geometric characteristics of the broken lines. Suppose the system behavior sequence is: For ρ∈(0,1), let: Then ε(d0,d i ) satisfies the grey relational axiom, where ρ is the discrimination coefficient; ε(d0,d i ) is d0 and d i Grey relational degree; Grey relational degree ε(d0,d i ) is often abbreviated as ε 0i , k-point correlation coefficient ε(d0(k),d i (k)) is abbreviated as ε 0i (k); Grey relational degree R i It is the weighted average value of the grey correlation coefficient: Among them, w k is the weight of indicator k; Step 5.2: Dancing warning algorithm: Analyze historical data to obtain three representative reference indicator sets for comparison: danger indicator set, intermediate indicator set and safety indicator set; when calculating, first use meteorological data and ledger information to obtain the 11 indicator values ​​required for warning to form the indicator set to be tested; then combine the weights obtained by the hierarchical analysis method to calculate the grey correlation between the indicator set to be tested and the three reference indicator sets; finally, according to the order of the correlation degree, the dancing warning result is obtained.

Images