A Method, System, Terminal and Medium for Monitoring and Analyzing Icing Galloping of High-Voltage Lines

Abstract

The present invention discloses a method, system, terminal and medium for monitoring and analyzing the galloping of ice-covered high-voltage lines, which relates to the technical field of power safety. The key points of its technical solution are as follows: The present invention collects actual vibration signals at relatively stable monitoring positions near the suspension points in high-voltage lines, which are less affected by the torque action. Based on the pre-constructed amplitude fitting function and wavelength fitting function, the vibration situation at the lowest point of the sag is inversely deduced. At the same time, in combination with the drag and lift determined by the angle of attack and wind speed values in the wind speed and wind direction data, with the goal that the resultant force direction of the drag, lift and gravity at the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment, the estimated ice-covered thickness at the lowest point of the sag corresponding to the minimum gravity error is solved. This method does not require collecting a large amount of real-time data, and can accurately and reliably determine the estimated ice-covered thickness through simple mechanical analysis, with low implementation difficulty.

Description

Technical Field

[0001] The present invention relates to the technical field of power safety, and more specifically, it relates to a method, system, terminal and medium for monitoring and analyzing icing galloping of high-voltage lines. Background Art

[0002] Icing galloping of high-voltage lines refers to a large-amplitude and low-frequency self-excited vibration phenomenon caused by the formation of a non-circular cross-section on the surface of the line after icing under specific meteorological conditions. This non-circular cross-section leads to unstable aerodynamic parameters. Icing can change the physical characteristics of the line, such as increasing the weight of the line and changing the aerodynamic characteristics of the line. Therefore, icing galloping of high-voltage lines may cause damage to the power supply line, such as serious accidents like wire breakage and tower collapse, affecting the stable operation of the power grid and power supply safety. Through monitoring, accidents can be detected and measures can be taken in a timely manner to reduce the occurrence of accidents. Monitoring of icing galloping of high-voltage lines is crucial for ensuring the safe, stable and reliable operation of the power system and is an indispensable part of the power industry.

[0003] Currently, the methods for monitoring icing galloping of high-voltage lines mainly include the monitoring method based on image recognition technology and the monitoring method based on supervised learning algorithms. For the monitoring method based on image recognition technology, in some areas with relatively complex terrain, there is often a large amount of fog in winter, and the video monitoring lens is prone to icing, making it difficult for the camera to accurately and clearly capture the icing image of the line, resulting in poor accuracy and reliability of the monitoring of icing galloping of high-voltage lines. For the monitoring method based on supervised learning algorithms, generally a large amount of meteorological data needs to be collected. In actual working conditions, high-voltage lines are often far from meteorological stations, and the prediction accuracy of the icing thickness in the adjacent area of the line channel is poor.

[0004] Therefore, how to research and design a method, system, terminal and medium for monitoring and analyzing icing galloping of high-voltage lines that can overcome the above defects is an urgent problem for us to solve currently. Summary of the Invention

[0005] To solve the deficiencies in the prior art, the purpose of the present invention is to provide a method, system, terminal and medium for monitoring and analyzing icing galloping of high-voltage lines. This method does not require collecting a large amount of real-time data, and can accurately and reliably determine the estimated icing thickness through simple mechanical analysis, with low implementation difficulty.

[0006] The above technical purpose of the present invention is achieved through the following technical solutions:

[0007] In the first aspect, a method for monitoring and analyzing icing galloping of high-voltage lines is provided, including the following steps:

[0008] Simulate and construct an amplitude fitting function that characterizes the amplitude change caused by vibration transmission and a wavelength fitting function that characterizes the wavelength change caused by vibration transmission in the ice-covered and dancing state of a high-voltage line;

[0009] Collect the actual vibration signals at the monitoring positions in the high-voltage line and the wind speed and wind direction data at the lowest point of the sag. The actual vibration signals include vertical vibration signals, lateral vibration signals, and horizontal vibration signals;

[0010] Combine the actual vibration signals, the horizontal length of the line between the monitoring position and the lowest point of the sag, the amplitude fitting function, and the wavelength fitting function to solve for the estimated vibration signal at the lowest point of the sag;

[0011] Divide the displacement directions of the lowest point of the sag in each vibration displacement segment from the estimated vibration signal;

[0012] Determine the angle of attack and the wind speed value based on the wind speed and wind direction data, and determine the drag and lift of the lowest point of the sag based on the angle of attack and the wind speed value;

[0013] With the goal that the resultant force direction of the drag, lift, and gravity at the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment, solve for the estimated ice-covering thickness at the lowest point of the sag corresponding to the minimum gravity error.

[0014] Furthermore, the construction process of the amplitude fitting function is specifically as follows:

[0015] Collect the amplitude values at each line position during the dancing of the simulated high-voltage line under different preset ice-covering thicknesses, and convert the amplitude values at each line position into a first relative parameter, where the first relative parameter is the ratio of the amplitude value at each line position to the amplitude value corresponding to the monitoring position;

[0016] Using the line position as the abscissa and the first relative parameter corresponding to the amplitude value of the line position as the ordinate, fit to obtain the amplitude fitting function corresponding to the preset ice-covering thickness;

[0017] And / or, the construction process of the wavelength fitting function is specifically as follows:

[0018] Collect the wavelengths of each changing vibration wave during the dancing of the simulated high-voltage line under different preset ice-covering thicknesses, and convert the wavelengths of each changing vibration wave into a second relative parameter, where the second relative parameter is the ratio of the wavelength of each changing vibration wave to the wavelength of the changing vibration wave corresponding to the monitoring position;

[0019] Using the serial number of the changing vibration wave as the abscissa and the second relative parameter corresponding to the wavelength of the changing vibration wave as the ordinate, fit to obtain the wavelength fitting function corresponding to the preset ice-covering thickness.

[0020] Furthermore, the expression of the amplitude fitting function is specifically as follows:

[0021] ;

[0022] Among them, represents the amplitude change coefficient at the lateral line position of the high-voltage line ; represents the amplitude fitting function simulated and constructed laterally when the preset ice coating thickness of the high-voltage line is ; And / or, the expression of the wavelength fitting function is specifically:

[0023] ;

[0024] ;

[0025] Among them, represents the wavelength change coefficient of the th changing vibration wave of the high-voltage line laterally ; represents the wavelength fitting function simulated and constructed laterally when the preset ice coating thickness of the high-voltage line is ; on the lateral direction

[0026] Furthermore, the solution process of the estimated vibration signal at the lowest point of the sag is specifically:

[0027] Input the horizontal length of the line between the monitoring position and the lowest point of the sag into the amplitude fitting function to calculate the amplitude change coefficient at the lowest point of the sag;

[0028] Calculate the estimated amplitude of the estimated vibration signal by multiplying the actual amplitude in the actual vibration signal by the amplitude change coefficient;

[0029] Solve for the wavelength of the changing vibration wave corresponding to the lowest point of the sag based on the wavelength fitting function, the actual wavelength at the monitoring position, and the horizontal length of the line between the monitoring position and the lowest point of the sag;

[0030] Determine the phase difference based on the position of the lowest point of the sag in the corresponding changing vibration wave, and calculate the estimated phase of the estimated vibration signal in combination with the initial phase in the actual vibration signal;

[0031] Keep the vibration frequency unchanged, and construct the estimated vibration signal at the lowest point of the sag by combining the estimated amplitude and the estimated phase.

[0032] Furthermore, the solution expression of the estimated vibration signal at the lowest point of the sag is:

[0033] ;

[0034] Among them, represents the actual wavelength at the monitoring position; Indicates the wavelength change coefficient of the th changing vibration wave in the lateral direction of the high-voltage line; The th; Indicates the horizontal length of the line between the monitoring position and the lowest point of the sag; Indicates the phase difference of the high-voltage line in the lateral direction; Indicates the sequence number of the changing vibration wave at the lowest point of the sag on the high-voltage line; Indicates the wavelength change coefficient of the th changing vibration wave in the lateral direction of the high-voltage line; The th; Indicates the actual vibration signal at the monitoring position; Indicates the actual amplitude in the actual vibration signal; Indicates the angular frequency; Indicates the time; Indicates the initial phase; Indicates the vibration signal at the lowest point of the sag in the lateral direction; Indicates the amplitude change coefficient of the lowest point of the sag on the high-voltage line in the lateral direction; Indicates the predicted vibration signal at the lowest point of the sag; Indicates the vibration signal at the lowest point of the sag in the horizontal direction; Indicates the vibration signal at the lowest point of the sag in the vertical direction.

[0035] Furthermore, the process of determining the displacement direction of the vibration displacement segment is specifically as follows:

[0036] Project the predicted vibration signal onto the horizontal-vertical plane to obtain a displacement signal, where the horizontal-vertical plane is a reference plane constructed perpendicular to the lateral direction;

[0037] Screen out the peak points and valley points in the displacement signal, and take the displacement segment between adjacent peak points and valley points as the vibration displacement segment. The direction from the endpoint with an earlier time to the endpoint with a later time in the vibration displacement segment is taken as the displacement direction.

[0038] Furthermore, the process of solving the predicted ice coating thickness at the lowest point of the sag corresponding to the minimum gravity error is specifically as follows:

[0039] Match the preset gravity corresponding to the preset ice coating thickness according to the amplitude fitting function and wavelength fitting function selected for solving the predicted vibration signal;

[0040] Solve for the predicted gravity at the lowest point of the sag by ensuring that the resultant direction of the resistance, lift, and gravity at the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment;

[0041] Calculate the corresponding gravity error based on the absolute value of the difference between the preset gravity and the estimated gravity, and select the preset gravity corresponding to the minimum gravity error among all gravity errors.

[0042] Use the preset ice-covered thickness corresponding to the selected preset gravity as the estimated ice-covered thickness at the lowest point of the sag.

[0043] In a second aspect, a monitoring and analysis system for ice galloping of high-voltage lines is provided. This system is used to implement a monitoring and analysis method for ice galloping of high-voltage lines as described in any one of the first aspects, including:

[0044] A simulation analysis module, used to simulate and construct an amplitude fitting function characterizing the amplitude change caused by vibration transmission and a wavelength fitting function characterizing the wavelength change in the ice galloping state of the high-voltage line.

[0045] A data acquisition module, used to collect the actual vibration signals at the monitoring positions in the high-voltage line and the wind speed and direction data at the lowest point of the sag. The actual vibration signals include vertical vibration signals, lateral vibration signals, and horizontal vibration signals.

[0046] A vibration tracing module, used to solve for the estimated vibration signal at the lowest point of the sag by combining the actual vibration signals, the horizontal length of the line between the monitoring position and the lowest point of the sag, the amplitude fitting function, and the wavelength fitting function.

[0047] A vibration segmentation module, used to divide the displacement directions of the lowest point of the sag in each vibration displacement segment from the estimated vibration signal.

[0048] An aerodynamic analysis module, used to determine the angle of attack and wind speed value based on the wind speed and direction data, and determine the drag and lift at the lowest point of the sag according to the angle of attack and wind speed value.

[0049] An ice-coverage estimation module, used to solve for the estimated ice-covered thickness at the lowest point of the sag corresponding to the minimum gravity error with the goal that the resultant force direction of the drag, lift, and gravity at the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment.

[0050] In a third aspect, a computer terminal is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements a monitoring and analysis method for ice galloping of high-voltage lines as described in any one of the first aspects.

[0051] In a fourth aspect, a computer-readable medium is provided, on which a computer program is stored. When the computer program is executed by a processor, it can implement a monitoring and analysis method for ice galloping of high-voltage lines as described in any one of the first aspects.

[0052] Compared with the prior art, the present invention has the following beneficial effects:

[0053] 1. A method for monitoring and analyzing the galloping of high-voltage lines provided by the present invention collects actual vibration signals at relatively stable monitoring positions near the suspension points in high-voltage lines, which are less affected by torque. Based on the pre-constructed amplitude fitting function and wavelength fitting function, it inversely deduces the vibration situation at the lowest point of the sag. At the same time, combining the resistance and lift determined by the angle of attack and wind speed values in the wind speed and direction data, with the goal that the resultant force direction of the resistance, lift, and gravity at the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment, it solves for the estimated ice coating thickness at the lowest point of the sag corresponding to the minimum gravity error. This method does not require collecting a large amount of real-time data, and can accurately and reliably determine the estimated ice coating thickness through simple mechanical analysis, with low implementation difficulty;

[0054] 2. When establishing the amplitude fitting function and wavelength fitting function in the present invention, the amplitude values at each line position and the wavelengths of each varying vibration wave are all converted into relative parameters, which is convenient for direct deformation processing based on the expression of the actual vibration signal, so as to quickly obtain the expression of the estimated vibration signal;

[0055] 3. When dividing the vibration displacement segments in the present invention, considering that the vertical direction is the main vibration direction of the high-voltage line and the vibration degree in the horizontal direction is weak, the estimated vibration signal is projected onto the transverse vertical plane to obtain the displacement signal, and the displacement direction of each vibration displacement segment is determined based on the peak points and valley points in the displacement signal, which can reduce the error influence in the mechanical analysis process and effectively improve the accuracy of solving the estimated ice coating thickness;

[0056] 4. When solving the estimated ice coating thickness at the lowest point of the sag in the present invention, by traversing and analyzing the amplitude fitting function and wavelength fitting function under different preset ice coating thicknesses, the preset gravity corresponding to the minimum gravity error is obtained, and the preset ice coating thickness corresponding to the selected preset gravity is used as the estimated ice coating thickness at the lowest point of the sag, effectively ensuring the rationality of the selection of the amplitude fitting function and wavelength fitting function. Description of the Drawings

[0057] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, form a part of this application, and do not constitute a limitation to the embodiments of the present invention. In the drawings:

[0058] Figure 1 is the flowchart in Embodiment 1 of the present invention;

[0059] Figure 2 is the system block diagram in Embodiment 2 of the present invention. Detailed Embodiments

[0060] To make the objectives, technical solutions, and advantages of the present invention more clear and understandable, the present invention will be further described in detail below with reference to the embodiments and the accompanying drawings. The illustrative embodiments of the present invention and their descriptions are only used to explain the present invention and are not intended to limit the present invention.

[0061] Embodiment 1: A method for monitoring and analyzing the galloping of a high-voltage line under ice coating, as Figure 1 shown, includes the following steps:

[0062] S1: Simulate and construct an amplitude fitting function that characterizes the amplitude change caused by vibration transmission and a wavelength fitting function that characterizes the wavelength change in the state of galloping of the high-voltage line under ice coating;

[0063] S2: Collect the actual vibration signals at the monitoring positions in the high-voltage line and the wind speed and wind direction data at the lowest point of the sag. The actual vibration signals include vertical vibration signals, lateral vibration signals, and horizontal vibration signals;

[0064] S3: Combine the actual vibration signals, the horizontal length of the line between the monitoring position and the lowest point of the sag, the amplitude fitting function, and the wavelength fitting function to solve for the estimated vibration signal at the lowest point of the sag;

[0065] S4: Divide the displacement directions of the lowest point of the sag in each vibration displacement segment from the estimated vibration signal;

[0066] S5: Determine the angle of attack and the wind speed value according to the wind speed and wind direction data, and determine the resistance and lift of the lowest point of the sag based on the angle of attack and the wind speed value;

[0067] S6: With the objective that the resultant force direction of the resistance, lift, and gravity of the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment, solve for the estimated ice coating thickness at the lowest point of the sag corresponding to the minimum gravity error.

[0068] In step S1, the simulation and construction of the amplitude fitting function and the wavelength fitting function can collect data through a wind tunnel test or can be simulated and analyzed using ANSYS software.

[0069] Taking ANSYS software as an example, the structural parameters of the high-voltage line are established through the finite element analysis method, including but not limited to data such as the height of the transmission tower, the line span, and the line sag.

[0070] The specific process of constructing the amplitude fitting function is as follows: Collect the amplitude values of each line position during the galloping of the simulated high-voltage line under different preset ice coating thicknesses, and convert the amplitude values of each line position into a first relative parameter. The first relative parameter is the ratio of the amplitude value of each line position to the amplitude value corresponding to the monitoring position; Fit the amplitude fitting function corresponding to the preset ice coating thickness with the line position as the abscissa and the first relative parameter corresponding to the amplitude value of the line position as the ordinate.

[0071] For example, if the amplitude value at line position 1 is 2m and the amplitude value at the monitoring position is 0.4m, then the first relative parameter of line position 1 is 5.

[0072] The specific expression of the amplitude fitting function is:

[0073] ;

[0074] where represents the amplitude change coefficient at the line position on the transverse of the high-voltage line; is the amplitude fitting function simulated and constructed on the transverse when the preset ice coating thickness of the high-voltage line is of the high-voltage line. It should be noted that the amplitude value of the line position may not be converted into the first relative parameter, and the least squares method can be directly sampled to perform curve fitting on the amplitude value of the line position, which is not restricted here.

[0075] The specific process of constructing the wavelength fitting function is as follows: collect the wavelengths of the varying vibration waves when the simulated high-voltage line dances under different preset ice coating thicknesses, and convert the wavelengths of the varying vibration waves into the second relative parameter, where the second relative parameter is the ratio of the wavelength of each varying vibration wave to the wavelength of the corresponding varying vibration wave at the monitoring position; use the serial number of the varying vibration wave as the abscissa and the second relative parameter corresponding to the wavelength of the varying vibration wave as the ordinate to fit the wavelength fitting function corresponding to the preset ice coating thickness.

[0076] The specific expression of the wavelength fitting function is:

[0077] The specific expression of the wavelength fitting function is:

[0078] ;

[0079] where represents the wavelength change coefficient of the th varying vibration wave on the transverse of the high-voltage line; is the wavelength fitting function simulated and constructed on the transverse when the preset ice coating thickness of the high-voltage line is of the high-voltage line.

[0080] Similarly, the wavelength of the varying vibration wave may not be converted into the second relative parameter, and the least squares method can be directly sampled to perform curve fitting on the wavelength of the varying vibration wave, which is not restricted here. In addition, the corresponding amplitude fitting functions and wavelength fitting functions are constructed for each of the vertical, transverse, and horizontal directions.

[0081] When establishing the amplitude fitting function and wavelength fitting function of the present invention, the amplitude values at each line position and the wavelengths of each changing vibration wave are converted into relative parameters, which facilitates direct deformation processing based on the expression of the actual vibration signal, so as to quickly obtain the expression of the predicted vibration signal.

[0082] In step S2, the monitoring position needs to ensure that effective vibration signals can be obtained, generally 5 - 20 m away from the suspension point. Vibration traceability processing is required for vertical vibration signals, lateral vibration signals, and horizontal vibration signals.

[0083] In step S3, the solution process of the predicted vibration signal at the lowest point of the sag is specifically as follows: The horizontal line length between the monitoring position and the lowest point of the sag is input into the amplitude fitting function to calculate the amplitude change coefficient at the lowest point of the sag; The predicted amplitude of the predicted vibration signal is calculated by multiplying the actual amplitude in the actual vibration signal by the amplitude change coefficient; The wavelength of the changing vibration wave corresponding to the lowest point of the sag is solved based on the wavelength fitting function, the actual wavelength at the monitoring position, and the horizontal line length between the monitoring position and the lowest point of the sag; The phase difference is determined according to the position of the lowest point of the sag in the corresponding changing vibration wave, and the predicted phase of the predicted vibration signal is calculated in combination with the initial phase in the actual vibration signal; Keeping the vibration frequency unchanged, the predicted vibration signal at the lowest point of the sag is constructed by combining the predicted amplitude and the predicted phase.

[0084] Taking the amplitude fitting function and wavelength fitting function constructed in step S1 as an example, the solution expression of the predicted vibration signal at the lowest point of the sag is:

[0085] ;

[0086] Among them, represents the actual wavelength at the monitoring position; represents the wavelength change coefficient of the th changing vibration wave of the high-voltage line in the lateral direction; represents the horizontal line length between the monitoring position and the lowest point of the sag; represents the phase difference of the high-voltage line in the lateral direction; represents the serial number of the changing vibration wave where the lowest point of the sag on the high-voltage line is located; represents the wavelength change coefficient of the th changing vibration wave of the high-voltage line in the lateral direction; represents the actual vibration signal at the monitoring position; represents the actual amplitude in the actual vibration signal; represents the angular frequency; represents the time; represents the initial phase; Indicates the vibration signal at the lowest point of the sag in the transverse direction ; Indicates the amplitude change coefficient of the lowest point of the sag of the high-voltage line in the transverse direction ; Indicates the predicted vibration signal at the lowest point of the sag Indicates the vibration signal at the lowest point of the sag in the horizontal direction ; Indicates the vibration signal at the lowest point of the sag in the vertical direction ;

[0087] In step S4, the process of determining the displacement direction of the vibration displacement segment is as follows: project the predicted vibration signal onto the transverse-vertical plane to obtain a displacement signal, where the transverse-vertical plane is a reference plane constructed perpendicular to the transverse direction; screen out the peak points and valley points in the displacement signal, and the displacement segment between adjacent peak points and valley points is used as the vibration displacement segment, and the direction from the endpoint with an earlier time to the endpoint with a later time in the vibration displacement segment is used as the displacement direction.

[0088] When dividing the vibration displacement segment in the present invention, considering that the vertical direction is the main vibration direction of the high-voltage line and the vibration degree in the horizontal direction is weak, the predicted vibration signal is projected onto the transverse-vertical plane to obtain a displacement signal, and the displacement direction of each vibration displacement segment is determined based on the peak points and valley points in the displacement signal, which can reduce the error influence in the mechanical analysis process and effectively improve the accuracy of solving the predicted ice coating thickness.

[0089] In step S5, the corresponding drag coefficient and lift coefficient can be matched from the database according to the angle of attack, and the drag coefficient and lift coefficient can be obtained through simulation analysis or wind tunnel tests.

[0090] In step S6, the process of solving the predicted ice coating thickness at the lowest point of the sag corresponding to the minimum gravity error is as follows: match the corresponding preset gravity according to the preset ice coating thickness corresponding to the amplitude fitting function and wavelength fitting function selected for solving the predicted vibration signal; solve the predicted gravity at the lowest point of the sag with the resultant force direction of the drag, lift, and gravity at the lowest point of the sag being consistent with the displacement direction of the corresponding vibration displacement segment; calculate the corresponding gravity error with the absolute value of the difference between the preset gravity and the predicted gravity, and select the preset gravity corresponding to the minimum gravity error among all gravity errors; use the preset ice coating thickness corresponding to the selected preset gravity as the predicted ice coating thickness at the lowest point of the sag.

[0091] When the present invention solves the estimated ice coating thickness at the lowest point of the sag, by traversing and analyzing the amplitude fitting function and the wavelength fitting function under different preset ice coating thicknesses, the preset gravity corresponding to the minimum gravity error is obtained, and the preset ice coating thickness corresponding to the selected preset gravity is used as the estimated ice coating thickness at the lowest point of the sag, effectively ensuring the rationality of the selection of the amplitude fitting function and the wavelength fitting function.

[0092] It should be noted that the lift is perpendicular to the application direction of the wind load, the drag is in the same direction as the application direction of the wind load, and the angle of attack is the angle between the application direction of the wind load and the axis direction of the high-voltage line.

[0093] Embodiment 2: A monitoring and analysis system for ice-induced galloping of high-voltage lines, which is used to implement a monitoring and analysis method for ice-induced galloping of high-voltage lines as described in Embodiment 1, as Figure 2 shown, including a simulation analysis module, a data acquisition module, a vibration source tracing module, a vibration segmentation module, an aerodynamic analysis module, and an ice coating estimation module.

[0094] Among them, the simulation analysis module is used to simulate and construct an amplitude fitting function characterizing the amplitude change caused by vibration transmission and a wavelength fitting function characterizing the wavelength change in the ice-induced galloping state of the high-voltage line; the data acquisition module is used to collect the actual vibration signals at the monitoring positions in the high-voltage line and the wind speed and wind direction data at the lowest point of the sag, and the actual vibration signals include vertical vibration signals, lateral vibration signals, and horizontal vibration signals; the vibration source tracing module is used to solve the estimated vibration signal at the lowest point of the sag by combining the actual vibration signals, the horizontal length of the line between the monitoring position and the lowest point of the sag, the amplitude fitting function, and the wavelength fitting function; the vibration segmentation module is used to divide the displacement direction of the lowest point of the sag in each vibration displacement segment from the estimated vibration signal; the aerodynamic analysis module is used to determine the angle of attack and the wind speed value according to the wind speed and wind direction data, and determine the drag and lift at the lowest point of the sag based on the angle of attack and the wind speed value; the ice coating estimation module is used to solve the estimated ice coating thickness at the lowest point of the sag corresponding to the minimum gravity error with the goal that the resultant force direction of the drag, lift, and gravity at the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment.

[0095] The present invention also records a computer terminal, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements a monitoring and analysis method for ice-induced galloping of high-voltage lines as described in Embodiment 1.

[0096] The present invention also records a computer-readable medium, on which a computer program is stored. When the computer program is executed by the processor, it can implement a monitoring and analysis method for ice-induced galloping of high-voltage lines as described in Embodiment 1.

[0097] Working principle: The present invention collects actual vibration signals at a relatively stable monitoring position near the suspension point in a high-voltage line, which is less affected by the torque effect. Based on the pre-constructed amplitude fitting function and wavelength fitting function, it inversely deduces the vibration situation at the lowest point of the sag. At the same time, combining the resistance and lift determined by the angle of attack and wind speed values in the wind speed and direction data, with the goal that the resultant force direction of the resistance, lift, and gravity at the lowest point of the sag is consistent with the displacement direction of the corresponding vibration displacement segment, it solves for the estimated ice coating thickness at the lowest point of the sag corresponding to the minimum gravity error. This method does not require collecting a large amount of real-time data, and can accurately and reliably determine the estimated ice coating thickness through simple mechanical analysis, with low implementation difficulty.

[0098] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0099] The present application is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to the embodiments of the present application. It should be understood that each flow and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, and the combination of the flows and / or blocks in the flowchart and / or block diagram can also be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate means for implementing the functions specified in Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.

[0100] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured article including instruction means, and the instruction means implements the functions specified in Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.

[0101] These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process. Thus, the instructions executed on the computer or other programmable device provide for implementing the functions in the processFigure 1 One process or multiple processes and / or boxes Figure 1 Steps of the functions specified in one box or multiple boxes.

[0102] The specific embodiments described above further elaborate on the purpose, technical solutions, and beneficial effects of the present invention. It should be understood that the above description is only the specific embodiments of the present invention and is not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.

Claims

1. A method for monitoring and analyzing ice dancing of high-voltage lines, characterized in that: The following steps are involved: Simulate and construct the amplitude fitting function to characterize the amplitude change caused by vibration transmission of high-voltage lines in the state of ice dancing, and the wavelength fitting function to characterize the wavelength change caused by vibration transmission; Collect the actual vibration signal of the monitoring position in the high-voltage line and the wind speed and direction data at the lowest point of the sag. The actual vibration signal includes vertical vibration signal, lateral vibration signal and horizontal vibration signal; The estimated vibration signal at the lowest point of the sag is obtained by combining the actual vibration signal, the horizontal length of the line between the monitoring position and the lowest point of the sag, the amplitude fitting function and the wavelength fitting function; Divide the displacement direction of the lowest point of the sag in each vibration displacement segment from the estimated vibration signal; Determine the wind attack angle and wind speed value according to the wind speed and direction data, and determine the drag and lift at the lowest point of the sag according to the wind attack angle and wind speed value; With the goal of keeping the combined force direction of drag, lift and gravity at the lowest point of the sag consistent with the displacement direction of the corresponding vibration displacement segment, the estimated ice thickness at the lowest point of the sag corresponding to the minimum gravity error is solved.

2. A method for monitoring and analyzing ice dancing of high-voltage lines according to claim 1, characterized in that: The construction process of the amplitude fitting function is specifically as follows: Collect the amplitude value of each line position when simulating the high-voltage line dancing under different preset ice thicknesses, and convert the amplitude value of each line position into a first relative parameter, which is the ratio of the amplitude value of each line position to the amplitude value corresponding to the monitoring position; The amplitude fitting function corresponding to the preset ice thickness is obtained by fitting the line position as the horizontal coordinate and the first relative parameter corresponding to the amplitude value of the line position as the vertical coordinate; And / or, the construction process of the wavelength fitting function is specifically as follows: The wavelength of each changing vibration wave when the simulated high-voltage line dances under different preset ice thicknesses is collected, and the wavelength of each changing vibration wave is converted into a second relative parameter, which is the ratio of the wavelength of each changing vibration wave to the wavelength of the changing vibration wave corresponding to the monitoring position; The wavelength fitting function corresponding to the preset ice thickness is obtained by fitting with the serial number of the changing vibration wave as the horizontal coordinate and the second relative parameter corresponding to the wavelength of the changing vibration wave as the vertical coordinate.

3. A method for monitoring and analyzing ice dancing of high-voltage lines according to claim 2, characterized in that: The expression of the amplitude fitting function is specifically: ； in, Indicates that the high voltage line is horizontal Upper line position The amplitude variation coefficient at ; Indicates that the preset ice thickness of the high-voltage line is When in horizontal The amplitude fitting function constructed by the above simulation; And / or, the expression of the wavelength fitting function is specifically: ； in, Indicates that the high voltage line is horizontal Previous The wavelength variation coefficient of the changing vibration wave; Indicates that the preset ice thickness of the high-voltage line is When in horizontal The wavelength fitting function constructed by the simulation above.

4. The method for monitoring and analyzing ice dancing of high-voltage lines according to claim 1 is characterized in that: The specific solution process of the estimated vibration signal at the lowest point of the sag is as follows: The horizontal length of the line between the monitoring position and the lowest point of the sag is input into the amplitude fitting function to calculate the amplitude variation coefficient of the lowest point of the sag; The estimated amplitude of the estimated vibration signal is calculated by the product of the actual amplitude in the actual vibration signal and the amplitude variation coefficient; The wavelength of the changing vibration wave corresponding to the lowest point of the sag is obtained based on the wavelength fitting function, the actual wavelength at the monitoring position, and the horizontal length of the line between the monitoring position and the lowest point of the sag; The phase difference is determined by the position of the lowest point of the sag in the corresponding changing vibration wave, and the estimated phase of the estimated vibration signal is calculated in combination with the initial phase in the actual vibration signal; Keeping the vibration frequency unchanged, the estimated vibration signal of the lowest point of the sag is constructed by combining the estimated amplitude and the estimated phase.

5. A method for monitoring and analyzing ice dancing of high-voltage lines according to claim 4, characterized in that: The estimated vibration signal solution expression at the lowest point of the sag is: ； in, Indicates the actual wavelength at the monitoring location; Indicates that the high voltage line is horizontal Previous The wavelength variation coefficient of the changing vibration wave; Indicates the horizontal length of the line between the monitoring position and the lowest point of the sag; Indicates that the high voltage line is horizontal The phase difference on Indicates the changing vibration wave number of the lowest point of sag on the high-voltage line; Indicates that the high voltage line is horizontal Previous The wavelength variation coefficient of the changing vibration wave; Indicates the actual vibration signal at the monitoring location; Represents the actual amplitude in the actual vibration signal; represents the angular frequency; Indicates time; represents the initial phase; Indicates that the lowest point of the sag is in the horizontal direction. Vibration signal on Indicates that the high voltage line is horizontal Amplitude variation coefficient of the lowest point of the upper sag; The estimated vibration signal indicating the lowest point of the sag; Indicates that the lowest point of the sag is horizontal Vibration signal on Indicates that the lowest point of the sag is vertical The vibration signal on.

6. The method for monitoring and analyzing ice dancing of high-voltage lines according to claim 1 is characterized in that: The displacement direction determination process of the vibration displacement section is specifically as follows: The displacement signal is obtained by projecting the estimated vibration signal onto the transverse plane, where the transverse plane is a reference plane constructed by the vertical and transverse directions; The peak points and trough points in the displacement signal are screened out, and the displacement segment between adjacent peak points and trough points is used as the vibration displacement segment, and the direction from the endpoint in front in time to the endpoint in back in time in the vibration displacement segment is used as the displacement direction.

7. The method for monitoring and analyzing ice dancing of high-voltage lines according to claim 1 is characterized in that: The process of solving the estimated ice thickness at the lowest sag point corresponding to the minimum gravity error is specifically as follows: The preset ice thickness corresponding to the amplitude fitting function and the wavelength fitting function selected by solving the estimated vibration signal matches the corresponding preset gravity; The estimated gravity at the lowest point of the sag is obtained by keeping the direction of the combined force of the drag, lift and gravity at the lowest point of the sag consistent with the displacement direction of the corresponding vibration displacement section; The corresponding gravity error is calculated by the absolute value of the difference between the preset gravity and the estimated gravity, and the preset gravity corresponding to the minimum gravity error among all gravity errors is selected; The preset ice thickness corresponding to the selected preset gravity is used as the estimated ice thickness at the lowest point of the sag.

8. A high-voltage line ice dancing monitoring and analysis system, characterized in that: The system is used to implement a high-voltage line ice dancing monitoring and analysis method as described in any one of claims 1 to 7, comprising: A simulation analysis module is used to simulate and construct an amplitude fitting function to characterize the amplitude change caused by vibration transmission of a high-voltage line in an ice-covered dancing state, as well as a wavelength fitting function to characterize the wavelength change caused by the vibration transmission; A data acquisition module is used to collect the actual vibration signal of the monitoring position in the high-voltage line and the wind speed and direction data at the lowest point of the sag. The actual vibration signal includes a vertical vibration signal, a lateral vibration signal and a horizontal vibration signal; The vibration tracing module is used to combine the actual vibration signal, the horizontal length of the line between the monitoring position and the lowest point of the sag, the amplitude fitting function and the wavelength fitting function to obtain the estimated vibration signal of the lowest point of the sag; A vibration segmentation module is used to divide the displacement direction of the lowest point of the sag in each vibration displacement segment from the estimated vibration signal; The aerodynamic analysis module is used to determine the wind attack angle and wind speed value according to the wind speed and direction data, and determine the drag and lift at the lowest point of the sag according to the wind attack angle and wind speed value; The ice estimation module is used to solve the estimated ice thickness at the lowest point of the sag corresponding to the minimum gravity error, with the goal of keeping the direction of the combined force of resistance, lift and gravity at the lowest point of the sag consistent with the displacement direction of the corresponding vibration displacement segment.

9. A computer terminal comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, a high-voltage line ice dancing monitoring and analysis method as described in any one of claims 1-7 is implemented.

10. A computer readable medium having a computer program stored thereon, characterized in that: The computer program is executed by a processor to implement a high-voltage line ice dancing monitoring and analysis method as described in any one of claims 1 to 7.

Images