A method and system for evaluating the risk of loss of transmission path section caused by strong convective gales

Abstract

A strong convective gale caused transmission channel cross section loss risk assessment method and system, including, based on the radar echo characteristic determination monitored object is in the strong convective gale type;From the strong convective gale type corresponding wind profile and wind speed monitoring data, the wind speed at the height where the monitored object is located is calculated;Then from two aspects, risk assessment, 1) with minimum electrical clearance distance as the index of wind deviation caused channel cross section loss risk judgment, that is, when the distance between the monitored object and the surrounding object is greater than the minimum electrical clearance distance, there is a channel cross section loss risk, 2) with the wind speed of the monitored object as the index of wind deviation caused channel cross section loss risk judgment, that is, the minimum wind speed corresponding to the maximum allowed wind deflection angle of the monitored object in the channel is inversely deduced with the minimum electrical clearance distance.The wind speed at the place where the monitored object is located is accurately calculated, and the maximum inversion wind speed leading to the wind deflection of the monitored object is taken as the basis for judgment, so that the strong convective gale risk is accurately evaluated.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of power transmission risk assessment, and particularly relates to a strong convective gale caused power transmission passage section loss risk assessment method and system. BACKGROUND

[0002] The existing strong convective gale risk assessment mostly adopts the 10-minute average wind speed at 10 meters above the ground to carry out wind load calculation according to the 4 basic wind profiles specified in the specification, and then assesses the wind deviation risk of the power transmission line. In fact, the strong convective gale such as tornado and downburst has a unique wind field structure, and its wind profile is different from the basic wind profile given in the specification. Therefore, the existing technology is not accurate in strong convective gale risk assessment. SUMMARY

[0003] In order to solve the problem of inaccurate strong convective gale risk assessment in the prior art, the present application provides a strong convective gale caused power transmission passage section loss risk assessment method, which comprises the following steps:

[0004] Determine the type of strong convective gale at the location of the monitored object based on the radar echo characteristics;

[0005] Calculate the wind speed at the height of the monitored object based on the wind profile corresponding to the type of strong convective gale and the wind speed monitoring data at the location of the monitored object;

[0006] Determine whether there is a wind deviation risk in the power transmission passage section where the monitored object is located based on the wind speed at the height of the monitored object by using any one of the following methods:

[0007] Method one: based on the wind speed at the height of the monitored object, calculate the wind deviation of the monitored object;

[0008] Method two: based on the wind speed at the height of the monitored object, and based on the minimum electrical clearance distance set between the monitored object and the surrounding objects, carry out back calculation.

[0009] Optionally, the calculation of the wind speed at the height of the monitored object based on the wind profile corresponding to the type of strong convective gale and the wind speed monitoring data at the location of the monitored object comprises:

[0010] According to the wind speed monitoring data at the location of the monitored object, interpolate the average wind speed at the set distance and set time from the ground at the location of the monitored object according to the inverse distance weighting method;

[0011] Based on the average wind speed at the set distance and set time from the ground at the location of the monitored object and the wind profile corresponding to the type of strong convective gale, obtain the wind speed at the height of the monitored object.

[0012] Optionally, the windage calculation is performed on the monitored object based on the wind speed at the height of the monitored object to determine whether there is a windage risk on the section of the power transmission channel where the monitored object is located, comprising:

[0013] The wind load of the monitored object is calculated based on the wind speed at the height of the monitored object and the parameters of the monitored object.

[0014] The windage risk assessment index is calculated based on the wind load of the monitored object and the parameters of the monitored object.

[0015] Based on the windage risk assessment index and the set minimum electrical clearance distance, it is determined whether there is a windage tripping risk.

[0016] Optionally, the monitored object includes an insulator string or an arbitrary point of a conductor.

[0017] Optionally, the wind load of the monitored object is calculated based on the wind speed at the height of the monitored object and the parameters of the monitored object, comprising:

[0018] When the monitored object is an insulator string, the wind load of the insulator string is calculated based on the conductor parameters, the wind speed at the height of the insulator string, the horizontal span and the angle between the wind direction and the line direction.

[0019] The conductor parameters include the wind pressure unevenness coefficient, the wind load body shape coefficient, the calculated outer diameter of the conductor and the cross-sectional area of the conductor.

[0020] Optionally, the windage risk assessment index is calculated based on the wind load of the monitored object and the parameters of the monitored object, comprising:

[0021] When the monitored object is an insulator string, the windage angle of the insulator string is calculated based on the wind load of the insulator string, the self-weight of the insulator string, the wind load of the arbitrary point of the conductor, the self-weight of the conductor and the cross-sectional area of the conductor.

[0022] The horizontal windage of the insulator string is calculated based on the windage angle of the insulator string and the length of the insulator string.

[0023] The shortest distance from the bottom vertex of the insulator string to the tower body is obtained by subtracting the horizontal windage of the insulator string from the distance between the top of the insulator string and the tower body, and the shortest distance from the bottom vertex of the insulator string to the tower body is taken as the windage risk assessment index of the insulator string.

[0024] Optionally, the wind load of the monitored object is calculated based on the wind speed at the height of the monitored object and the parameters of the monitored object, comprising:

[0025] When the monitored object is an arbitrary point on a conductor, the wind load at the arbitrary point on the conductor is calculated based on the conductor parameters and the wind speed at the height of the arbitrary point on the conductor;

[0026] The conductor parameters include: conductor wind pressure unevenness coefficient, wind load body shape coefficient, calculated outer diameter of the conductor and conductor cross-sectional area.

[0027] Optionally, the calculating of the windage risk assessment index based on the wind load of the monitored object in combination with the parameters of the monitored object includes:

[0028] When the monitored object is any point on the conductor, the wind deflection angle at any point on the conductor is calculated based on the conductor's own weight and the wind load at any point on the conductor;

[0029] Calculate the sag at any point on the conductor based on the span, the distance from any point on the conductor to the hanging point, and the maximum sag;

[0030] Calculate the horizontal windage and vertical windage at any point on the conductor based on the windage angle at any point on the conductor and the sag at any point on the conductor;

[0031] Based on the horizontal windage and vertical windage at any point on the conductor, the minimum distance between any point on the conductor and surrounding objects after windage occurs is calculated using three-dimensional spatial geometric distance, and the minimum distance is used as a windage risk assessment indicator for any point on the conductor.

[0032] Optionally, the wind load at any point on the conductor is calculated as follows:

[0033]

[0034] Where g4 is the wind load at any point on the conductor, α is the conductor wind pressure unevenness coefficient, C is the wind load shape coefficient, D is the calculated outer diameter of the conductor, v is the wind speed, A is the conductor cross-sectional area, and ρ is the air density.

[0035] Optionally, the wind load of the insulator string is calculated as follows:

[0036]

[0037] Where, P j is the wind load of the insulator string, j is the serial number of the insulator string, L s is the horizontal span, θ is the angle between wind direction and line direction, α is the conductor wind pressure unevenness coefficient, C is the wind load body shape coefficient, D is the calculated outer diameter of the conductor, v is the wind speed, and ρ is the air density.

[0038] Optionally, the wind deflection angle of the insulator string is calculated as follows:

[0039]

[0040] Where, is the wind deflection angle of the insulator string, L s , L c are horizontal and vertical span respectively; g1, g4 are the gravity of the arbitrary point of the conductor and the wind load of the arbitrary point of the conductor respectively; G j , P j are the weight of the insulator string and the wind load of the insulator string, j is the serial number of the insulator string, A is the cross-sectional area of the conductor, tg -1 is arctanx, d is the split distance of the conductor, and λ is the length of the insulator string.

[0041] Optionally, the wind deflection angle of the arbitrary point of the conductor is calculated according to the following formula:

[0042]

[0043] In the formula, ξ is the wind deflection angle of the arbitrary point of the conductor, and g1 and g4 are the gravity of the conductor and the wind load of the arbitrary point of the conductor respectively.

[0044] Optionally, the sag of the arbitrary point of the conductor is calculated according to the following formula:

[0045]

[0046] In the formula, f x is the sag of the arbitrary point, f m is the maximum sag, l is the span, l x is the distance between the arbitrary point of the conductor and the hanging point, and x is the arbitrary point.

[0047] Optionally, the horizontal wind deflection and the vertical wind deflection of the arbitrary point of the conductor are calculated according to the following formula:

[0048]

[0049] In the formula, d x is the horizontal wind deflection of the arbitrary point of the conductor, f x is the sag of the arbitrary point, ξ is the wind deflection angle of the arbitrary point of the conductor, λ is the length of the insulator string, is the wind deflection angle of the insulator string, d y is the vertical wind deflection of the arbitrary point of the conductor, and x is the arbitrary point.

[0050] Optionally, based on the wind speed at the height of the monitored object, the minimum electrical clearance distance set between the monitored object and the surrounding object is inversely calculated to determine whether there is a wind deflection risk in the cross section of the power transmission channel where the monitored object is located, comprising:

[0051] For each line in the channel, the maximum allowed wind deflection angle of the monitored object is inversely calculated based on the set minimum electrical clearance distance;

[0052] calculate the minimum wind speed of each line in the channel to cause the monitored object to yaw when the maximum allowable yaw angle of the monitored object is reached;

[0053] take the maximum value of the minimum wind speed of each line in the channel to cause the monitored object to yaw as the minimum wind speed corresponding to the maximum allowable yaw angle of the monitored object in the channel;

[0054] when the wind speed at the height of the monitored object is greater than the minimum wind speed corresponding to the maximum allowable yaw angle of the monitored object in the channel, the cross section of the power transmission channel where the monitored object is located has a risk of yawing.

[0055] In still another aspect, the application further provides a system for evaluating the risk of losing the cross section of a power transmission channel caused by strong convective winds, comprising:

[0056] a strong convective wind type determination module configured to determine the type of strong convective wind at the location of the monitored object based on radar echo characteristics;

[0057] a wind speed calculation module configured to calculate the wind speed at the height of the monitored object based on the wind profile corresponding to the type of strong convective wind and the wind speed monitoring data at the location of the monitored object;

[0058] an evaluation module configured to perform yaw calculation on the monitored object based on the wind speed at the height of the monitored object to determine whether the cross section of the power transmission channel where the monitored object is located has a risk of yawing, or to perform back calculation based on the wind speed at the height of the monitored object and the minimum electrical clearance distance set for the monitored object and the surrounding objects to determine whether the cross section of the power transmission channel where the monitored object is located has a risk of yawing.

[0059] Optionally, the evaluation module comprises:

[0060] a yaw calculation sub-module configured to perform yaw calculation on the monitored object based on the wind speed at the height of the monitored object to determine whether the cross section of the power transmission channel where the monitored object is located has a risk of yawing;

[0061] a back calculation sub-module configured to perform back calculation on the minimum electrical clearance distance set for the monitored object and the surrounding objects based on the wind speed at the height of the monitored object to determine whether the cross section of the power transmission channel where the monitored object is located has a risk of yawing.

[0062] Optionally, the yaw calculation sub-module comprises:

[0063] a wind load calculation unit configured to calculate the wind load of the monitored object based on the wind speed at the height of the monitored object and the parameters of the monitored object;

[0064] The index calculation unit is configured to calculate a wind deviation risk assessment index based on the wind load of the monitored object and the parameter of the monitored object.

[0065] The risk assessment unit is configured to determine whether there is a wind deviation trip-out risk based on the wind deviation risk assessment index and a set minimum electrical clearance distance.

[0066] Optionally, the wind load calculation unit comprises:

[0067] The insulator string wind load calculation subunit is configured to calculate the wind load of the insulator string based on the conductor parameter, the wind speed at the height where the insulator string is located, the horizontal span and the angle between the wind direction and the line direction when the monitored object is an insulator string.

[0068] The conductor arbitrary point wind load calculation subunit is configured to calculate the wind load of the conductor arbitrary point based on the conductor parameter and the wind speed at the height where the conductor arbitrary point is located when the monitored object is a conductor arbitrary point.

[0069] The conductor parameter comprises a conductor wind pressure unevenness coefficient, a wind load body shape coefficient, a calculated outer diameter of the conductor and a conductor cross-sectional area.

[0070] Optionally, the conductor arbitrary point wind load calculation subunit calculates the wind load of the conductor arbitrary point by the following formula:

[0071]

[0072] In the formula, g4 is the wind load of the conductor arbitrary point, α is the conductor wind pressure unevenness coefficient, C is the wind load body shape coefficient, D is the calculated outer diameter of the conductor, v is the wind speed, A is the conductor cross-sectional area, and ρ is the air density.

[0073] Optionally, the insulator string wind load calculation subunit calculates the wind load of the insulator string by the following formula:

[0074]

[0075] In the formula, P is the wind load of the insulator string, j is the serial number of the insulator string, L is the horizontal span, θ is the angle between the wind direction and the line direction, α is the conductor wind pressure unevenness coefficient, C is the wind load body shape coefficient, D is the calculated outer diameter of the conductor, v is the wind speed, and ρ is the air density. j s

[0076] Optionally, the index calculation unit comprises:

[0077] ​​The insulator string index calculation subunit is configured to, when the monitored object is an insulator string, calculate an insulator string wind deflection angle based on an insulator string wind load, an insulator string self weight, a conductor arbitrary point wind load, a conductor self weight and a conductor cross-sectional area; calculate an insulator string horizontal wind deflection based on the insulator string wind deflection angle and a length of the insulator string; obtain a shortest distance from a bottom vertex of the insulator string to a tower body by subtracting the insulator string horizontal wind deflection from a distance from a top of the insulator string to the tower body; and take the shortest distance from the bottom vertex of the insulator string to the tower body as a wind deflection risk assessment index of the insulator string.

[0078] The conductor arbitrary point index calculation subunit is configured to, when the monitored object is a conductor arbitrary point, calculate a conductor arbitrary point wind deflection angle based on a conductor self weight and a conductor arbitrary point wind load; calculate a conductor arbitrary point sag based on a span, a distance between the conductor arbitrary point and a hanging point and a maximum sag; calculate a conductor arbitrary point horizontal wind deflection and a conductor arbitrary point vertical wind deflection based on the conductor arbitrary point wind deflection angle and the conductor arbitrary point sag; and calculate a minimum distance between the conductor arbitrary point and surrounding objects after wind deflection by using three-dimensional space geometric distance based on the conductor arbitrary point horizontal wind deflection and the conductor arbitrary point vertical wind deflection, and take the minimum distance as a wind deflection risk assessment index of the conductor arbitrary point.

[0079] Optionally, the insulator string index calculation subunit calculates the insulator string wind deflection angle by the following formula:

[0080]

[0081] In the formula, ξ is the insulator string wind deflection angle, g1 and g4 are respectively the conductor self weight and the conductor arbitrary point wind load, L is the length of the insulator string, j is the insulator string serial number, A is the conductor cross-sectional area, and tg is arctanx. L is the length of the insulator string. s L is the length of the insulator string. c L is the length of the insulator string. j L is the length of the insulator string. j L is the length of the insulator string. -1 L is the length of the insulator string.

[0082] Optionally, the conductor arbitrary point index calculation subunit calculates the conductor arbitrary point wind deflection angle by the following formula:

[0083]

[0084] In the formula, ξ is the conductor arbitrary point wind deflection angle, g1 and g4 are respectively the conductor self weight and the conductor arbitrary point wind load.

[0085] Optionally, the conductor arbitrary point index calculation subunit calculates the conductor arbitrary point sag by the following formula:

[0086] L is the length of the insulator string.

[0087] wherein f x is the sag at an arbitrary point, f m is the maximum sag, l is the span, l x is the distance between an arbitrary point of the conductor and the hanging point, and x is the arbitrary point.

[0088] Optionally, the arbitrary point of the conductor index calculation subunit calculates the horizontal wind deflection and the vertical wind deflection at the arbitrary point of the conductor by the following formula:

[0089]

[0090] wherein d x is the horizontal wind deflection at an arbitrary point of the conductor, f x is the sag at an arbitrary point, ξ is the wind deflection angle at the arbitrary point of the conductor, λ is the length of the insulator string, is the wind deflection angle of the insulator string, d y is the vertical wind deflection at an arbitrary point of the conductor.

[0091] Optionally, the back-stepping sub-module is specifically configured to:

[0092] back-step the maximum allowable wind deflection angle of the monitored object based on the set minimum electrical clearance distance for each line in the channel;

[0093] calculate the minimum wind speed at which each line in the channel causes the wind deflection of the monitored object to reach the maximum allowable wind deflection angle of the monitored object;

[0094] take the maximum value in the minimum wind speed at which each line in the channel causes the wind deflection of the monitored object as the minimum wind speed corresponding to the maximum allowable wind deflection angle of the monitored object in the channel;

[0095] when the wind speed at the height of the monitored object is greater than the minimum wind speed corresponding to the maximum allowable wind deflection angle of the monitored object in the channel, the cross section of the power transmission channel in which the monitored object is located has a wind deflection risk.

[0096] In still another aspect, the present application also provides a computing device, comprising: at least one processor and a memory;

[0097] the memory is configured to store one or more programs;

[0098] when the one or more programs are executed by the at least one processor, a strong convective gale induced power transmission channel cross section loss risk assessment method as described above is implemented.

[0099] In still another aspect, the present application also provides a computer readable storage medium having a computer program stored thereon, wherein the computer program is executed to implement a strong convective gale induced power transmission channel cross section loss risk assessment method as described above.

[0100] Compared with the prior art, the application has the beneficial effects that:

[0101] The application provides a strong convective gale caused transmission channel section loss risk assessment method, which comprises the following steps: determining the type of strong convective gale at the place where a monitored object is located based on radar echo characteristics; calculating the wind speed at the height where the monitored object is located based on the wind profile corresponding to the type of strong convective gale and the wind speed monitoring data at the place where the monitored object is located; performing wind deviation calculation on the monitored object based on the wind speed at the height where the monitored object is located to determine whether there is a wind deviation risk of the transmission channel section where the monitored object is located; or performing back calculation based on the wind speed at the height where the monitored object is located and the minimum electrical clearance distance set between the monitored object and surrounding objects to determine whether there is a wind deviation risk of the transmission channel section where the monitored object is located. The application can more accurately calculate the wind speed at the place where the monitored object is located, and uses the maximum inversion wind speed causing the wind deviation of the monitored object in the channel as the basis for judging the wind deviation caused channel section loss risk, so that the strong convective gale risk assessment is more accurate. BRIEF DESCRIPTION OF DRAWINGS

[0102] Figure 1 The application provides a strong convective gale caused transmission channel section loss risk assessment method flow chart;

[0103] Figure 2 The application provides a downburst wind profile schematic diagram;

[0104] Figure 3 The application provides a wind deviation angle schematic diagram of a suspension insulator string and a conductor;

[0105] Figure 4 The application provides a suspension insulator string and tower body wind deviation distance schematic diagram;

[0106] Figure 5 The application provides a conductor to tree wind deviation distance schematic diagram;

[0107] Figure 6 The application provides a squall line wind radar echo characteristic;

[0108] Figure 7 The application provides a multi-cell thunderstorm radar echo characteristic;

[0109] Figure 8 The application provides a tornado or downburst echo characteristic. DETAILED DESCRIPTION

[0110] The application provides a strong convective gale caused transmission channel section loss risk assessment method, which accurately calculates the applied wind load, uses a wind deviation generalization model to quickly calculate the wind deviation angle, and further realizes the assessment of the risk of the transmission channel caused by the strong convective gale leading to the wind deviation tripping of all lines and causing the loss of the section.

[0111] Embodiment 1:

[0112] A strong convective gale caused transmission channel cross section loss risk assessment method, as shown in Figure 1 , comprising:

[0113] Step 1: Determine the type of strong convective gale at the location of the monitored object based on radar echo characteristics;

[0114] Step 2: Based on the wind profile corresponding to the type of strong convective gale and the wind speed monitoring data at the location of the monitored object, calculate the wind speed at the height of the monitored object;

[0115] Step 3: Based on the wind speed at the height of the monitored object, calculate the wind deflection of the monitored object to determine whether there is a wind deflection risk at the cross section of the transmission channel where the monitored object is located; or based on the wind speed at the height of the monitored object, based on the minimum electrical clearance distance set between the monitored object and the surrounding objects, perform back calculation to determine whether there is a wind deflection risk at the cross section of the transmission channel where the monitored object is located.

[0116] The following further introduces each step:

[0117] A strong convective gale caused transmission channel cross section loss risk assessment method, containing a strong convective gale type research and judgment method based on radar echo characteristics, a preferred wind field model method based on the type of strong convective gale, a conductor hanging point height wind speed calculation method based on monitoring data, and a transmission channel wind deflection caused cross section loss risk assessment method based on digital twinning.

[0118] Step 1: Determine the type of strong convective gale at the location of the monitored object based on radar echo characteristics, specifically including:

[0119] Based on radar echo characteristics to judge the type of strong wind, mainly including four types of strong convective gale: squall line wind, thunderstorm wind, downburst, and tornado.

[0120] The radar echo characteristics and the type of strong convective gale have the following corresponding relationship:

[0121] If the radar echo is a strip structure, the reflectivity factor is greater than 45 dBz, and the horizontal scale is about tens to hundreds of kilometers long and about tens of kilometers to two hundred kilometers wide, it is a squall line wind;

[0122] If there are multiple discrete cylindrical or spherical echo clusters in the radar echo, the reflectivity factor is greater than 45 dBz, and the range is at least 2500 km 2 , it is a multi-cell thunderstorm wind system;

[0123] If there is a bow echo in the radar echo, it is a tornado or a downburst; the bow echo in the northern hemisphere is counterclockwise, with a near-ground diameter of 25 meters to several hundred meters, an air diameter of several kilometers, and a radar echo reflectivity factor greater than 55 dBz, which is a tornado; the horizontal scale of the bow echo is generally 1 km to 10 km, and the reflectivity factor is between 45 dBz and 55 dBz, which is a downburst.

[0124] Before introducing step 2, the wind profile of different types of strong wind is further introduced:

[0125] The wind speed profile of squall line or multi-cell thunderstorm gale is:

[0126]

[0127] In the formula: z w is the height of the wind measurement point from the ground (m), usually 10 m; z is the height of the calculation point (m); is the 10 min average wind speed at the sampling point height (m / s), usually the 10 min average wind speed at 10 m above the ground; v z is the 10 min average wind speed at the calculation point (m / s); a is the roughness index, and 0.15 is taken for class B topography (according to GB 50009).

[0128] The wind speed profile of tornado is preferably the kuo-wen three-dimensional tornado model considering the boundary layer effect, and the parameters include the distance from the tornado center, the maximum tangential wind speed corresponding radius, the maximum tangential wind speed and the height:

[0129] T(η,r)=f(r)[1-e -πη cos(2bπη)]

[0130] R(η,r)=f(r){0.672e -πη sin[(b+1)πη]}

[0131] W(η,r)=93r 3 exp(-5r)V max [1-e -πη cos(2bπη)]

[0132]

[0133] In the formula: T(η,r), R(η,r), and W(η,r) are tangential velocity, radial velocity, and vertical velocity, respectively; V max is the maximum tangential velocity; b=1.2e -0.8r^4 ; f(r) is the height; η is the relative height.

[0134] Relative height: η = z / δ(r'), where z is the height at which the wind speed is to be calculated, and δ(r') is the thickness of the tornado boundary layer.

[0135] Tornado boundary layer thickness: δ(r')=δ0[1-exp(-0.5r 2 )]

[0136] Where: r = r' / r max , r is the ratio of the distance between the tower location and the tornado center to the radius corresponding to the maximum tangential wind speed; r' is the distance between the tower location and the tornado center; r max is the radius corresponding to the maximum tangential wind speed; δ0 is the boundary layer thickness of a tornado when r>>1, which is generally taken as 457m.

[0137] like Figure 2 The wind speed profile of the downburst shown is preferably the Wood & Kwok model which mainly considers the radial wind speed:

[0138]

[0139] Where δ is 0.5V max The corresponding height (m) is generally 6 times the height of the maximum wind speed; erf is the error function.

[0140] Step 2: Calculate the wind speed at the height of the monitored object based on the wind profile corresponding to the severe convective wind type and the wind speed monitoring data at the location of the monitored object, specifically including:

[0141] Based on the wind speed monitoring data at the location of the monitored object, the average wind speed at the location of the monitored object at a set distance from the ground and a set time is interpolated using the inverse distance weighted method;

[0142] The wind speed at the height of the monitored object is obtained based on the average wind speed at the set distance from the ground for the set time and the wind profile corresponding to the strong convective wind type.

[0143] Specifically, the monitoring objects include: insulator strings and arbitrary points on conductors. This embodiment takes the monitoring objects of insulator strings and arbitrary points on conductors as an example to further introduce this solution:

[0144] Based on the real-time monitoring data of the average wind speed at a height of 10m above the ground for 10 minutes with a spatial resolution of 1km×1km, the average wind speed at a height of 10m above the ground for 10 minutes at the tower location is interpolated using the inverse distance weighted method. Combined with the wind type and the wind profile corresponding to the severe convective wind type, the wind speed is converted to the height of the conductor hanging point.

[0145] Includes windage risk assessment of any point on insulator strings and conductors.

[0146] Before calculating the wind speed at the height of the monitored object, a three-dimensional digital geometric model of the important power transmission channel is established, that is:

[0147] The three-dimensional digital geometric model of the important power transmission channel is established by using power transmission line design data and operation and maintenance records. It includes a refined model of the tower-insulator string-conductor system, terrain, and surrounding object model.

[0148] Step 3: Based on the wind speed at the height of the monitored object, wind deflection calculation is performed on the monitored object to determine whether there is a wind deflection risk in the section of the power transmission channel where the monitored object is located; or based on the wind speed at the height of the monitored object, the minimum electrical clearance distance between the monitored object and the surrounding objects is calculated to determine whether there is a wind deflection risk in the section of the power transmission channel where the monitored object is located. Specifically, there are two ways to determine whether the section of the power transmission channel where the monitored object is located has a wind deflection risk:

[0149] Method one: the wind load of the monitored object is calculated based on the wind speed at the height of the monitored object and the parameters of the monitored object, including:

[0150] When the monitored object is an insulator string, the wind load of the insulator string is calculated based on the conductor parameters, the wind speed at the height of the insulator string, the horizontal span, and the angle between the wind direction and the line direction.

[0151] Wherein, the conductor parameters include: conductor wind pressure uneven coefficient, wind load body type coefficient, conductor calculation outer diameter and conductor cross-sectional area.

[0152] Further, the wind deflection risk assessment index is calculated based on the wind load of the monitored object and the parameters of the monitored object, including:

[0153] When the monitored object is an insulator string, the insulator string wind deflection angle is calculated based on the insulator string wind load, the insulator string self weight, the conductor arbitrary point wind load, the conductor self weight and the conductor cross-sectional area.

[0154] Based on the insulator string wind deflection angle and the length of the insulator string, the horizontal wind deflection of the insulator string is calculated.

[0155] The shortest distance from the bottom vertex of the insulator string to the tower body is obtained by subtracting the horizontal wind deflection of the insulator string from the distance between the top of the insulator string and the tower body, and the shortest distance from the bottom vertex of the insulator string to the tower body is taken as the wind deflection risk assessment index of the insulator string.

[0156] Further, the wind load of the monitored object is calculated based on the wind speed at the height of the monitored object and the parameters of the monitored object, including:

[0157] When the monitored object is an arbitrary point of the conductor, the wind load of the arbitrary point of the conductor is calculated based on a conductor parameter and a wind speed at a height where the arbitrary point of the conductor is located.

[0158] The conductor parameter includes a conductor wind pressure unevenness coefficient, a wind load shape coefficient, a calculated outer diameter of the conductor, and a conductor cross-sectional area.

[0159] Further, the wind deflection risk assessment index is calculated based on the wind load of the monitored object in combination with the parameter of the monitored object, and includes:

[0160] When the monitored object is an arbitrary point of the conductor, the wind deflection angle of the arbitrary point of the conductor is calculated based on a conductor self-weight and the wind load of the arbitrary point of the conductor.

[0161] The sag of the arbitrary point of the conductor is calculated based on a span, a distance between the arbitrary point of the conductor and a hanging point, and a maximum sag.

[0162] The horizontal wind deflection and the vertical wind deflection of the arbitrary point of the conductor are calculated based on the wind deflection angle of the arbitrary point of the conductor and the sag of the arbitrary point of the conductor.

[0163] Based on the horizontal wind deflection and the vertical wind deflection of the arbitrary point of the conductor, the minimum distance between the arbitrary point of the conductor and surrounding objects after wind deflection is calculated by using three-dimensional space geometry distance, and the minimum distance is taken as the wind deflection risk assessment index of the arbitrary point of the conductor.

[0164] The wind deflection angle of the insulator string and the arbitrary point of the conductor is as shown in Figure 3 , and the specific calculation formula is as follows:

[0165] The wind deflection angle of the insulator string is calculated as:

[0166]

[0167] In the formula, L s and L c are horizontal and vertical spans (m) respectively; g1 and g4 are the conductor self-weight and the wind load of the arbitrary point of the conductor (kg / m·mm 2 ); j is the serial number of the insulator string; G j and P j are the weight of the insulator string and the wind load of the insulator string (kg); and A is the cross-sectional area of the conductor (mm 2 ). Wherein:

[0168]

[0169] In the formula, g4 is the wind load of the arbitrary point of the conductor, j is the serial number of the insulator string, and a is the conductor wind pressure unevenness coefficient; C is the wind load shape coefficient; D is the calculated outer diameter of the conductor (mm); and L sm is the horizontal span, v is the wind speed, θ is the angle between the wind direction and the line direction, and ρ is the air density.

[0170] Table 1 Wind pressure unevenness coefficient

[0171] Design wind speed (m / s) v＜20 20≤v＜30 30≤v＜35 v≥35 Wind pressure unevenness coefficient 1 0.85 0.75 0.7

[0172] Insulator string horizontal wind deviation d λ = 0.5 m

[0173]

[0174] where λ is the insulator string length (m), is the insulator string wind deviation angle.

[0175] Conductor arbitrary point wind deviation angle calculation:

[0176]

[0177] where g1 and g4 are the conductor arbitrary point self-weight and conductor arbitrary point wind load (kg / m·mm 2 ) respectively.

[0178] Conductor arbitrary point sag f x = 0.5 m

[0179]

[0180] where x is the horizontal distance from the tower (m), l is the span (m), f m is the maximum sag (m), and l x is the distance between the conductor arbitrary point and the hanging point.

[0181] Conductor arbitrary point horizontal wind deviation d x , vertical wind deviation d y , respectively:

[0182]

[0183] Insulator string and tower wind deviation risk assessment:

[0184] When the shortest distance d1 from the bottom vertex of the suspension insulator string to the tower is less than the minimum electrical clearance distance [d] specified in the specification, there is a risk of insulator string wind deviation tripping the tower. Wherein d1 is obtained by subtracting the insulator string horizontal wind deviation from the distance between the insulator string top and the tower, as shown in Figure 4 .

[0185] Conductor risk assessment of wind deviation on buildings, trees, railways, roads, etc.:

[0186] According to the above steps, the position space coordinate information of the conductor under the wind deviation state is calculated, and on the basis of the three-dimensional coordinate system, the minimum distance d between the conductor and the surrounding trees after the wind deviation is calculated by using the three-dimensional space geometry r As shown in Figure 5 When d r is less than the minimum electrical clearance distance [d] specified in the specification, there is a risk of conductor wind deviation discharge triggering tripping of surrounding objects.

[0187] The wind deviation tripping risk of all lines in the power transmission channel is evaluated, and when there is a possibility of tripping, it is considered that there is a risk of loss of channel section, and an alarm is issued.

[0188] Method two: based on the wind speed at the height of the monitored object, the minimum electrical clearance distance between the monitored object and the surrounding object is calculated to determine whether the power transmission channel section where the monitored object is located has a wind deviation risk, comprising:

[0189] For each line in the channel, the maximum allowed wind deviation angle of the monitored object is calculated based on the set minimum electrical clearance distance;

[0190] The minimum wind speed of each line in the channel to make the monitored object wind deviation when reaching the maximum allowed wind deviation angle of the monitored object is calculated;

[0191] The maximum value of the minimum wind speed of each line in the channel to make the monitored object wind deviation is taken as the minimum wind speed corresponding to the maximum allowed wind deviation angle of the monitored object in the channel;

[0192] When the wind speed at the height of the monitored object is greater than the minimum wind speed corresponding to the maximum allowed wind deviation angle of the monitored object in the channel, the power transmission channel section where the monitored object is located has a wind deviation risk.

[0193] Embodiment 2:

[0194] In the following, taking an insulator string as an example, the above-mentioned embodiment two is used to judge whether the power transmission channel section where the monitored object is located has a wind deviation risk, and the present application is further introduced:

[0195] 1. Identify the wind field type based on the radar echo characteristics

[0196] (a) Squall line wind radar echo characteristics, as shown in Figure 6 ;

[0197] (b) Multi-cell thunderstorm radar echo characteristics, as shown in Figure 7 ;

[0198] (c) Tornado or downburst echo characteristics, as shown in Figure 8 .

[0199] 2. According to the measured 10m height 10min average wind speed, according to the judgment of the wind field type, select the wind profile, calculate the insulator string height wind speed V1.

[0200] 3. Calculate the minimum wind speed V2 of each line in the channel that causes the insulator string wind deviation; that is, the minimum wind speed V2 corresponding to the maximum allowed wind deviation angle in the channel where the insulator is located.

[0201] 4. When V1 is greater than V2, it is judged that the channel section is lost.

[0202] Compared with the boundary layer wind profile specified in the specification, the wind profile model based on strong convective wind type identification and differentiation can calculate the wind speed at the insulator string height more accurately. At the same time, the maximum inversion wind speed that causes the wind deviation of the insulator string of the transmission line in the channel is used as the basis for judging the risk of wind deviation causing the channel section loss, which is more accurate.

[0203] Embodiment 3:

[0204] A strong convective wind induced transmission channel section loss risk assessment method, including a strong convective wind type identification method based on radar echo characteristics, an optimal wind field model method based on strong convective wind type, a conductor hanging point height wind speed calculation method based on monitoring data, and a transmission channel wind deviation induced section loss risk assessment method based on digital twinning.

[0205] The strong convective wind type is identified based on radar echo characteristics, mainly including four types of strong convective wind: squall line wind, thunderstorm wind, downburst, and tornado.

[0206] The radar echo characteristics and the strong convective wind type have the following corresponding relationship:

[0207] If the radar echo is a strip structure, the reflectivity factor is greater than 45dBz, and the horizontal scale is about tens to hundreds of kilometers long and tens of kilometers to two hundred kilometers wide, it is a squall line wind;

[0208] If the radar echo exists multiple discrete cylindrical or spherical echo clusters, the reflectivity factor is greater than 45dBz, and the range is at least 2500km 2 It is a multi-cell thunderstorm wind system;

[0209] If the radar echo exists an arc-shaped echo, it is a tornado or a downburst; the arc-shaped echo in the northern hemisphere is counterclockwise, the near-ground diameter is 25 meters to several hundred meters, the air diameter is several kilometers, and the radar echo reflectivity factor is greater than 55dBz, which is a tornado; the horizontal scale of the arc-shaped echo is generally 1km to 10km, and the reflectivity factor is between 45dBz and 55dBz, which is a downburst.

[0210] The wind profile of different strong wind types is optimized as follows:

[0211] The wind speed profile of squall line or multi-cell thunderstorm gale is:

[0212]

[0213] wherein z w is the height of the measuring point from the ground (m), usually 10 m; z is the height of the calculating point (m); is the 10 min average wind speed of the sampling point (m / s), usually the 10 min average wind speed at the height of 10 m from the ground; v z is the 10 min average wind speed of the calculating point (m / s); a is the roughness index, 0.15 for B type landform (according to GB 50009).

[0214] The wind speed profile of tornado is preferably the Kuo-Wen three-dimensional tornado model considering the boundary layer effect, and the parameters include the distance from the tornado center, the radius corresponding to the maximum tangential wind speed, the maximum tangential wind speed and the height:

[0215] T(η,r)=f(r)[1-e -πη cos(2bπη)]

[0216] R(η,r)=f(r){0.672e -πη sin[(b+1)πη]}

[0217] W(η,r)=93r 3 exp(-5r)V max [1-e -πη cos(2bπη)]

[0218]

[0219] wherein T(η,r), R(η,r) and W(η,r) are respectively the tangential velocity, the radial velocity and the vertical velocity; V max is the maximum tangential velocity; b=1.2e -0.8r^4 ; f(r) is the tangential wind speed outside the boundary layer height; η is the relative height.

[0220] Relative height: η=z / δ(r'), z is the height at which the wind speed is calculated, and δ(r') is the tornado boundary layer thickness.

[0221] Tornado boundary layer thickness: δ(r')=δ0[1-exp(-0.5r 2 )]

[0222] wherein r=r' / r max , r is the ratio of the distance of the tower from the tornado center to the radius corresponding to the maximum tangential wind speed; r' is the distance of the tower from the tornado center; r maxis the maximum tangential wind speed corresponding to the radius; δ0 is the boundary layer thickness of the tornado when r >> 1, generally taken as 457 m.

[0223] As shown in the downburst, the wind speed profile of the downburst is preferably a Wood & Kwok model mainly considering the radial wind speed as shown in the following formula: Figure 2

[0224]

[0225] wherein, δ is 0.5V max corresponding height (m), generally taken as 6 times the height where the maximum wind speed is located; erf is an error function.

[0226] According to the real-time monitoring data of 10-minute average wind speed at a height of 10 m from the ground with a spatial resolution of 1 km x 1 km, the 10-minute average wind speed at a height of 10 m from the ground at the location of the tower is interpolated according to the inverse distance weighting method, and the wind speed at the height of the conductor hanging point is converted according to the above-mentioned selected wind profile in combination with the type of strong wind.

[0227] The wind deflection risk assessment of any point of the conductor and insulator is shown in the following formula: Figure 3

[0228] ① The three-dimensional digital geometric model of the important power transmission channel is established by using the power transmission line design data and operation and maintenance records. The three-dimensional digital geometric model includes a fine model of the power transmission tower-insulator-conductor system, terrain, and surrounding object model of the channel.

[0229] ② Insulator string wind deflection angle calculation:

[0230]

[0231] wherein, L s , L c are the horizontal and vertical span (m) respectively; g1, g4 are the self-weight and wind load of the conductor (kg / m·mm 2 ); G j , P j are the weight and wind load of the insulator string (kg); A is the cross-sectional area of the conductor (mm 2 ).

[0232] wherein:

[0233] The wind load at any point of the conductor is calculated according to the following formula:

[0234]

[0235] wherein, g4 is the wind load at any point of the conductor, α is the wind pressure unevenness coefficient of the conductor, C is the wind load shape coefficient, D is the calculation outer diameter of the conductor, v is the wind speed, A is the cross-sectional area of the conductor, and ρ is the air density. ​​

[0236] The wind load of the insulator string is calculated by the following formula:

[0237]

[0238] In the formula, P j is the wind load of the insulator string, j is the serial number of the insulator string, L s is the horizontal span, θ is the angle between the wind direction and the line direction, α is the wind pressure unevenness coefficient of the conductor, C is the body shape coefficient of the wind load, D is the calculated outer diameter of the conductor, v is the wind speed, and ρ is the air density.

[0239] Table 1 Wind pressure unevenness coefficient

[0240]

[0241] ③ Horizontal wind deviation d λ of the insulator string is:

[0242]

[0243] In the formula, λ is the length of the insulator (m).

[0244] ④ Wind deviation angle at any point in the span of the conductor is calculated by:

[0245]

[0246] In the formula, g1 and g4 are the self-weight and wind load of the conductor (kg / m·mm 2 ).

[0247] The sag f x at any point is:

[0248]

[0249] In the formula, x is the horizontal distance from the tower (m); L is the span (m); f m is the maximum sag (m).

[0250] ⑤ Horizontal wind deviation d x and vertical wind deviation d y at any point in the span of the conductor are respectively:

[0251]

[0252] In the formula, d x is the horizontal wind deviation at any point of the conductor, f x is the sag at any point, ξ is the wind deviation angle at any point of the conductor, λ is the length of the insulator string, is the wind deviation angle of the insulator string, d y is the vertical wind deviation at any point of the conductor, and x is the arbitrary point.

[0253] (6) Insulator and tower body windage risk assessment:

[0254] As shown in Figure 4 , when the shortest distance d1 from the bottom vertex of the overhanging insulator string to the tower body is less than the minimum electrical clearance distance [d] specified in the specification, there is a risk of insulator windage tripping of the tower body.

[0255] (7) Conductor against buildings, trees, railways, roads, etc. Windage risk assessment:

[0256] As shown in Figure 5 , the position spatial coordinate information of the conductor under windage is calculated according to the above steps, and the minimum distance d r between the conductor and the surrounding trees after windage is calculated based on the three-dimensional coordinate system. r When d r is less than the minimum electrical clearance distance [d] specified in the specification, there is a risk of conductor windage discharge causing tripping of the surrounding objects.

[0257] (8) Assess the windage tripping risk of all lines in the power transmission channel. When all of them have the possibility of tripping, it is considered that there is a risk of loss of channel section, and an alarm is issued.

[0258] Example 4:

[0259] In another aspect, the application also provides a strong convective wind induced power transmission channel section loss risk assessment system, comprising:

[0260] A strong convective wind type determination module is configured to determine the type of strong convective wind at the location of the monitored object based on radar echo characteristics.

[0261] A wind speed calculation module is configured to calculate the wind speed at the height of the monitored object based on the wind profile corresponding to the type of strong convective wind and the wind speed monitoring data at the location of the monitored object.

[0262] An assessment module is configured to perform windage calculation on the monitored object based on the wind speed at the height of the monitored object to determine whether there is a windage risk in the power transmission channel section where the monitored object is located, or to perform reverse calculation on the set minimum electrical clearance distance based on the wind speed at the height of the monitored object to determine whether there is a windage risk in the power transmission channel section where the monitored object is located.

[0263] Optionally, the assessment module comprises:

[0264] A windage calculation sub-module is configured to perform windage calculation on the monitored object based on the wind speed at the height of the monitored object to determine whether there is a windage risk in the power transmission channel section where the monitored object is located.

[0265] An inverse calculation module is configured to inversely calculate the set minimum electrical clearance distance based on the wind speed at the height of the monitored object to determine whether there is a windage yaw risk in the cross section of the power transmission channel where the monitored object is located.

[0266] Optionally, the windage yaw calculation submodule comprises:

[0267] A wind load calculation unit is configured to calculate the wind load of the monitored object based on the wind speed at the height of the monitored object and the parameters of the monitored object.

[0268] An index calculation unit is configured to calculate a windage yaw risk assessment index based on the wind load of the monitored object in combination with the parameters of the monitored object.

[0269] A risk assessment unit is configured to determine whether there is a windage yaw tripping risk based on the windage yaw risk assessment index in combination with the set minimum electrical clearance distance.

[0270] Optionally, the wind load calculation unit comprises:

[0271] An insulator string wind load calculation subunit is configured to calculate the wind load of the insulator string based on the conductor parameters, the wind speed at the height of the insulator string, the horizontal span and the angle between the wind direction and the line direction when the monitored object is an insulator string.

[0272] A conductor arbitrary point wind load calculation subunit is configured to calculate the wind load of the conductor arbitrary point based on the conductor parameters and the wind speed at the height of the conductor arbitrary point when the monitored object is a conductor arbitrary point.

[0273] The conductor parameters comprise a conductor wind pressure unevenness coefficient, a wind load shape coefficient, a calculated outer diameter of the conductor and a conductor cross-sectional area.

[0274] Optionally, the conductor arbitrary point wind load calculation subunit calculates the wind load of the conductor arbitrary point by the following formula:

[0275]

[0276] In the formula, g4 is the wind load of the conductor arbitrary point, α is the conductor wind pressure unevenness coefficient, C is the wind load shape coefficient, D is the calculated outer diameter of the conductor, v is the wind speed, A is the conductor cross-sectional area, and ρ is the air density.

[0277] Optionally, the insulator string wind load calculation subunit calculates the wind load of the insulator string by the following formula:

[0278]

[0279] In the formula, P is the wind load of the insulator string, j is the serial number of the insulator string, L is the horizontal span, v is the wind speed, C is the wind load shape coefficient, D is the calculated outer diameter of the conductor, and ρ is the air density. j s ​is the horizontal span, θ is the angle between wind direction and line direction, α is the conductor wind pressure unevenness coefficient, C is the wind load body shape coefficient, D is the calculated outer diameter of the conductor, v is the wind speed, and ρ is the air density.

[0280] Optionally, the indicator calculation unit includes:

[0281] An insulator string index calculation subunit is configured to, when the monitored object is an insulator string, calculate the wind deflection angle of the insulator string based on the wind load of the insulator string, the deadweight of the insulator string, the wind load at any point on the conductor, the deadweight of the conductor, and the cross-sectional area of ​​the conductor; calculate the horizontal wind deflection of the insulator string based on the wind deflection angle of the insulator string and the length of the insulator string; and subtract the horizontal wind deflection of the insulator string from the distance between the top of the insulator string and the tower body to obtain the shortest distance from the bottom vertex of the insulator string to the tower body, and use the shortest distance from the bottom vertex of the insulator string to the tower body as the wind deflection risk assessment index of the insulator string;

[0282] The conductor arbitrary point index calculation subunit is used to calculate the wind deflection angle of any point on the conductor based on the conductor's own weight and the wind load at any point on the conductor when the monitored object is any point on the conductor; calculate the sag at any point on the conductor based on the span, the distance between any point on the conductor and the hanging point, and the maximum sag; calculate the horizontal wind deflection and vertical wind deflection at any point on the conductor based on the wind deflection angle and the sag at any point on the conductor; and calculate the minimum distance between any point on the conductor and the surrounding objects after wind deflection occurs using three-dimensional spatial geometric distance based on the horizontal wind deflection and vertical wind deflection at any point on the conductor, and use the minimum distance as the wind deflection risk assessment index for any point on the conductor.

[0283] Optionally, the insulator string index calculation subunit calculates the insulator string wind deflection angle using the following formula:

[0284]

[0285] Where, is the wind deflection angle of the insulator string, L s , L c are the horizontal and vertical spans respectively; g1 and g4 are the conductor deadweight and wind load respectively; G j 、P j is the weight of the insulator string and its wind load, j is the serial number of the insulator string, A is the cross-sectional area of ​​the conductor, tg -1 is arctanx, d is the conductor splitting spacing, and λ is the insulator string length.

[0286] Optionally, the conductor arbitrary point index calculation subunit calculates the wind deflection angle of the conductor arbitrary point using the following formula:

[0287]

[0288] In the formula, ξ is the wind deflection angle at any point of the conductor, g1 and g4 are the self-weight of the conductor and the wind load at any point of the conductor respectively.

[0289] Optionally, the conductor arbitrary point index calculation subunit calculates the sag at any point of the conductor by the following formula:

[0290]

[0291] In the formula, f x is the sag at any point, f m is the maximum sag, l is the span, l x is the distance between the arbitrary point of the conductor and the hanging point, and x is the arbitrary point.

[0292] Optionally, the conductor arbitrary point index calculation subunit calculates the horizontal wind deflection and the vertical wind deflection at any point of the conductor by the following formula:

[0293]

[0294] In the formula, d x is the horizontal wind deflection at any point of the conductor, f x is the sag at any point, ξ is the wind deflection angle at any point of the conductor, λ is the length of the insulator string, is the wind deflection angle of the insulator string, d y is the vertical wind deflection at any point of the conductor.

[0295] Optionally, the backstepping calculation subunit is specifically configured to:

[0296] backstepping the maximum allowable wind deflection angle of the monitored object based on the set minimum electrical clearance distance for each line in the channel;

[0297] calculating the minimum wind speed at which each line in the channel causes the wind deflection of the monitored object when reaching the maximum allowable wind deflection angle of the monitored object;

[0298] taking the maximum value in the minimum wind speed at which each line in the channel causes the wind deflection of the monitored object as the minimum wind speed corresponding to the maximum allowable wind deflection angle of the monitored object in the channel;

[0299] when the wind speed at the height of the monitored object is greater than the minimum wind speed corresponding to the maximum allowable wind deflection angle of the monitored object in the channel, the cross section of the power transmission channel where the monitored object is located has a wind deflection risk.

[0300] Embodiment 5:

[0301] Based on the same inventive concept, the present application further provides a computer device, which comprises a processor and a memory, the memory is used to store a computer program, the computer program comprises program instructions, and the processor is used to execute the program instructions stored in the computer storage medium. The processor can be a central processing unit (CPU), and can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc., which are the computing core and control core of the terminal, and are suitable for implementing one or more instructions, and are specifically suitable for loading and executing one or more instructions in the computer storage medium to implement a corresponding method flow or a corresponding function, so as to implement the steps of the strong convection gale induced transmission path section loss risk assessment method in the above embodiment.

[0302] Embodiment 6:

[0303] Based on the same inventive concept, the present application further provides a storage medium, specifically a computer readable storage medium (Memory), which is a memory device in the computer device, and is used to store programs and data. It can be understood that the computer readable storage medium herein can include the built-in storage medium in the computer device, and of course can also include the expansion storage medium supported by the computer device. The computer readable storage medium provides a storage space, and the storage space stores the operating system of the terminal. Moreover, one or more instructions suitable for being loaded and executed by the processor are also stored in the storage space, and the instructions can be one or more computer programs (including program codes). It should be noted that the computer readable storage medium herein can be a high-speed RAM memory, or a non-volatile memory such as at least one disk memory. One or more instructions stored in the computer readable storage medium can be loaded and executed by the processor to implement the steps of the strong convection gale induced transmission path section loss risk assessment method in the above embodiment.

[0304] Those skilled in the art will appreciate that embodiments of the present application can be readily used as a method, a system, or a computer program product. Accordingly, the present application can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application can take the form of a computer program product on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROMs, optical storage devices, and the like) embodying computer readable program code.

[0305] The present application is described in reference to the flowchart illustrations and / or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processing system or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks.

[0306] These computer program instructions can also be stored in a computer- readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instructions which implement the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks.

[0307] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks.

[0308] The foregoing is merely illustrative of the principles of this application and various modifications can be made by those skilled in the art without departing from the scope and spirit of the application. The above specification, examples and data provide a complete description of the manufacture and use of the composition of the application. The compositions herein have structural requirements as set forth in the description, figures, and claims. Since many modifications, variations, and changes in detail can be made to the described embodiments, it is intended that only such modifications and variations are within the scope of the application as disclosed in the foregoing description, examples, and the appended claims.

Claims

1. A method for evaluating the risk of loss of section of a transmission path due to severe convective gales, characterized in that, The method comprises the following steps: determining a strong convective gale type at the location of the monitored object based on radar echo characteristics; calculating the wind speed at the height of the monitored object based on the wind profile corresponding to the strong convective gale type and the wind speed monitoring data at the location of the monitored object; performing wind deviation calculation on the monitored object based on the wind speed at the height of the monitored object to determine whether there is a wind deviation risk in the transmission channel section where the monitored object is located; or performing back calculation based on the minimum electrical clearance distance set between the monitored object and the surrounding objects based on the wind speed at the height of the monitored object to determine whether there is a wind deviation risk in the transmission channel section where the monitored object is located; the calculation of the wind speed at the height of the monitored object based on the wind profile corresponding to the strong convective gale type and the wind speed monitoring data at the location of the monitored object comprises: interpolating the wind speed monitoring data at the location of the monitored object to the average wind speed at the set distance from the ground and the set time at the location of the monitored object according to the inverse distance weighting method; obtaining the wind speed at the height of the monitored object based on the average wind speed at the set distance from the ground and the set time at the location of the monitored object and the wind profile corresponding to the strong convective gale type; the wind deviation calculation on the monitored object based on the wind speed at the height of the monitored object to determine whether there is a wind deviation risk in the transmission channel section where the monitored object is located comprises: calculating the wind load of the monitored object based on the wind speed at the height of the monitored object and the parameters of the monitored object; calculating the wind deviation risk evaluation index based on the wind load of the monitored object and the parameters of the monitored object; determining whether there is a wind deviation tripping risk based on the minimum electrical clearance distance set and the wind deviation risk evaluation index; the monitored object comprises an insulator string and an arbitrary point of a conductor; the back calculation based on the minimum electrical clearance distance set between the monitored object and the surrounding objects based on the wind speed at the height of the monitored object to determine whether there is a wind deviation risk in the transmission channel section where the monitored object is located comprises: back calculating the maximum allowable wind deviation angle of the monitored object based on the minimum electrical clearance distance set for each line in the channel; calculating the minimum wind speed at which each line in the channel causes the wind deviation of the monitored object when the maximum allowable wind deviation angle of the monitored object is reached; taking the maximum value in the minimum wind speed at which each line in the channel causes the wind deviation of the monitored object as the minimum wind speed corresponding to the maximum allowable wind deviation angle of the monitored object in the channel; when the wind speed at the height of the monitored object is greater than the minimum wind speed corresponding to the maximum allowable wind deviation angle of the monitored object in the channel, the transmission channel section where the monitored object is located has a wind deviation risk.

2. The method of claim 1, wherein, the calculation of the wind load of the monitored object based on the wind speed at the height of the monitored object and the parameters of the monitored object comprises: when the monitored object is an insulator string, calculating the wind load of the insulator string based on the conductor parameters, the wind speed at the height of the insulator string, the horizontal span and the angle between the wind direction and the line direction. The conductor parameters include the wind pressure unevenness coefficient of the conductor, the wind load body shape coefficient, the calculated outer diameter of the conductor and the cross-sectional area of the conductor.

3. The method of claim 1, wherein, The calculating of the windage risk assessment index based on the wind load of the monitored object and the parameters of the monitored object includes: When the monitored object is an insulator string, the wind deflection angle of the insulator string is calculated based on the wind load of the insulator string, the dead weight of the insulator string, the wind load at any point on the conductor, the dead weight of the conductor, and the cross-sectional area of ​​the conductor; Calculating the horizontal windage of the insulator string based on the windage angle of the insulator string and the length of the insulator string; The distance between the top of the insulator string and the tower body is subtracted from the horizontal wind deviation of the insulator string to obtain the shortest distance from the bottom vertex of the insulator string to the tower body. The shortest distance from the bottom vertex of the insulator string to the tower body is used as the wind deviation risk assessment index of the insulator string.

4. The method of claim 1, wherein, The calculating of the wind load of the monitored object based on the wind speed at the height of the monitored object and the parameters of the monitored object includes: When the monitored object is an arbitrary point on a conductor, the wind load at the arbitrary point on the conductor is calculated based on the conductor parameters and the wind speed at the height of the arbitrary point on the conductor; The conductor parameters include: conductor wind pressure unevenness coefficient, wind load body shape coefficient, calculated outer diameter of the conductor and conductor cross-sectional area.

5. The method of claim 1, wherein, The calculating of the windage risk assessment index based on the wind load of the monitored object and the parameters of the monitored object includes: When the monitored object is any point on the conductor, the wind deflection angle at any point on the conductor is calculated based on the conductor's own weight and the wind load at any point on the conductor; Calculate the sag at any point on the conductor based on the span, the distance from any point on the conductor to the hanging point, and the maximum sag; Calculate the horizontal windage and vertical windage at any point on the conductor based on the windage angle at any point on the conductor and the sag at any point on the conductor; Based on the horizontal windage and vertical windage at any point on the conductor, the minimum distance between any point on the conductor and surrounding objects after windage occurs is calculated using three-dimensional spatial geometric distance, and the minimum distance is used as a windage risk assessment indicator for any point on the conductor.

6. The method of claim 4, wherein, The wind load at any point of the conductor is calculated as follows: Where g4 is the wind load at any point on the conductor, α is the conductor wind pressure unevenness coefficient, C is the wind load shape coefficient, D is the calculated outer diameter of the conductor, v is the wind speed, A is the conductor cross-sectional area, and ρ is the air density.

7. The method of claim 2, wherein, The wind load of the insulator string is calculated as follows: In the formula, P j is the wind load of the insulator string, j is the serial number of the insulator string, L s is the horizontal span, θ is the angle between the wind direction and the line direction, α is the wind pressure unevenness coefficient of the conductor, C is the body shape coefficient of the wind load, D is the calculated outer diameter of the conductor, and v is the wind speed, and ρ is the air density.

8. The method of claim 3, wherein, The wind deflection angle of the insulator string is calculated as follows: wherein L is the windage angle of the insulator string s L c are the horizontal and vertical span, respectively; g1, g4 are the self-weight of the conductor and the wind load at any point of the conductor, respectively; G j P j is the weight of the insulator string and the wind load of the insulator string, j is the serial number of the insulator string, A is the cross-sectional area of the conductor, d is the split distance of the conductor, and λ is the length of the insulator string.

9. The method of claim 5, wherein, The wind deflection angle at any point of the conductor is calculated as follows: Where ξ is the wind deflection angle at any point on the conductor, g1 and g4 are the conductor weight and wind load at any point on the conductor, respectively.

10. The method of claim 5, wherein, The sag at any point of the conductor is calculated as follows: where f x is the sag at any point, f m is the maximum sag, l is the span, l x is the distance from any point on the conductor to the point of suspension, x is any point.

11. The method of claim 5, wherein, The horizontal windage and vertical windage at any point on the conductor are calculated as follows: where d x is the horizontal wind deflection at an arbitrary point of the conductor, f x is the sag at an arbitrary point, ξ is the wind deflection angle at an arbitrary point of the conductor, λ is the length of the insulator string, is the wind deflection angle of the insulator string, d y is the vertical wind deflection at an arbitrary point of the conductor, x is an arbitrary point.

12. A system for implementing the method for assessing the risk of loss of cross section of the convection and windage-induced transmission path according to any one of claims 1 to 11, characterized in that, include: A strong convective wind type determination module is used to determine the strong convective wind type at the location of the monitored object based on radar echo characteristics; A wind speed calculation module, configured to calculate the wind speed at the height of the monitored object based on the wind profile corresponding to the severe convective wind type and the wind speed monitoring data at the location of the monitored object; The evaluation module is configured to: based on the wind speed at the height of the monitored object, perform wind deviation calculation on the monitored object to determine whether there is a wind deviation risk on the section of the power transmission channel where the monitored object is located; or based on the wind speed at the height of the monitored object, perform back-stepping calculation based on the minimum electrical clearance distance set between the monitored object and the surrounding object to determine whether there is a wind deviation risk on the section of the power transmission channel where the monitored object is located.

13. The system of claim 12, wherein, The evaluation module comprises: a wind deviation calculation submodule configured to: based on the wind speed at the height of the monitored object, perform wind deviation calculation on the monitored object to determine whether there is a wind deviation risk on the section of the power transmission channel where the monitored object is located; a back-stepping calculation submodule configured to: based on the wind speed at the height of the monitored object, perform back-stepping calculation on the minimum electrical clearance distance set to determine whether there is a wind deviation risk on the section of the power transmission channel where the monitored object is located.

14. The system of claim 13, wherein, The wind deviation calculation submodule comprises: a wind load calculation unit configured to: based on the wind speed at the height of the monitored object and the parameters of the monitored object, calculate the wind load of the monitored object; an index calculation unit configured to: based on the wind load of the monitored object and the parameters of the monitored object, calculate a wind deviation risk evaluation index; a risk evaluation unit configured to: based on the wind deviation risk evaluation index and the minimum electrical clearance distance set, determine whether there is a wind deviation tripping risk.

15. A computer device, comprising: comprises: at least one processor and a memory; the memory is configured to store one or more programs; when the one or more programs are executed by the at least one processor, a method for evaluating the loss risk of a power transmission channel section caused by strong convective gales is implemented, as claimed in any one of claims 1 to 11.

16. A computer-readable storage medium, characterized in that, a computer program is stored thereon, and when the computer program is executed, a method for evaluating the loss risk of a power transmission channel section caused by strong convective gales is implemented, as claimed in any one of claims 1 to 11.

Images