Methods, systems, and media for assessing lightning strike risk of transmission lines based on multi-source data

Abstract

This invention relates to the field of power equipment risk assessment technology, and discloses a method, system, and medium for assessing the lightning strike risk of transmission lines based on multi-source data. The method includes: acquiring three-dimensional laser point cloud data, lightning statistics data, and meteorological data of the transmission line; the three-dimensional laser point cloud data includes terrain parameters of the tower locations; calculating the conductor spatial position of the electrical geometric model of the transmission line under the influence of dynamic wind deflection based on the meteorological data, as well as the dynamic wind deflection correction coefficient of the electrical geometric model; and accurately calculating the tripping rate of the transmission line, especially the backflashover tripping rate, based on the terrain correction coefficient corresponding to the terrain parameters, the conductor spatial position, the dynamic wind deflection correction coefficient, and the lightning statistics data, thereby accurately assessing the lightning strike risk of the line.

Description

Technical Field

[0001] This invention relates to the field of power equipment risk assessment technology, and in particular to a method, system and medium for assessing the lightning strike risk of transmission lines based on multi-source data. Background Technology

[0002] Lightning strikes are one of the main causes of power transmission line tripping, seriously affecting the safe and stable operation of the power grid. Statistics show that lightning strikes account for 40%-70% of high-voltage transmission line faults. Lightning strikes mainly take two forms: backflashover and lightning strikes. For 500kV and above lines, due to their higher insulation levels, the risk of backflashover is relatively low, making backflashover the primary threat. Lightning strikes not only cause insulator flashovers and equipment damage, but can also trigger large-scale power outages, impacting social and economic development and residential electricity use.

[0003] Currently, the lightning strike risk assessment of transmission lines mainly relies on the following methods: (1) Procedure method: The lightning strike rate is calculated based on empirical formulas, but it cannot reflect the influence of micro-topography and dynamic weather. (2) Electrical geometry model: The lightning attraction range is analyzed through geometric relationships, but traditional EGM does not consider factors such as conductor wind deflection and dynamic terrain. (3) Improved electrical geometry model: A height correction coefficient is introduced, but it still relies on two-dimensional slice calculation, which is difficult to handle complex terrain and multiple circuits on the same tower. (4) Data-driven method: The lightning location system and meteorological data are combined, but the real-time performance is insufficient and there is a lack of multi-source data fusion mechanism.

[0004] Existing technologies generally suffer from problems such as a single data source and static models, mostly relying on lightning location system data and 3D laser point cloud data, without considering dynamic factors such as weather and wind deviation. Summary of the Invention

[0005] Therefore, it is necessary to propose a method, system, and medium for assessing the lightning strike risk of transmission lines based on multi-source data to address the above problems.

[0006] A method for assessing the lightning strike risk of transmission lines based on multi-source data, the method comprising:

[0007] Acquire three-dimensional laser point cloud data, lightning statistics data, and meteorological data of the transmission line. The three-dimensional laser point cloud data includes terrain parameters of the tower location.

[0008] The spatial position of the conductor in the electrical geometric model of the transmission line under the influence of dynamic wind deflection is calculated based on meteorological data, as well as the dynamic wind deflection correction coefficient of the electrical geometric model.

[0009] The transmission line tripping rate is calculated based on the terrain correction coefficient corresponding to the terrain parameters, the spatial position of the conductor, the dynamic wind deflection correction coefficient, and the lightning statistics.

[0010] The lightning strike risk level of the transmission line is determined based on the tripping rate of the transmission line.

[0011] The meteorological data includes wind deflection amplitude and wind speed. The calculation of the conductor spatial position of the transmission line's electrical geometric model under the influence of dynamic wind deflection based on the meteorological data, and the dynamic wind deflection correction coefficient of the electrical geometric model, specifically includes:

[0012] according to Calculate the spatial position of the conductors in the electrical geometry model of the transmission line under the influence of dynamic wind deflection, where y( (t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, t is time, the spatial position of the conductor changes with time, H is the horizontal tension, w is the weight per unit length, A is the wind deflection amplitude, and x is the horizontal tension. Let be the projected coordinates of the conductor on the horizontal plane.

[0013] according to Calculate the dynamic wind deflection correction coefficient for the electrical geometry model, where, This represents the dynamic wind deflection correction coefficient for the electrical geometry model. For real-time wind speed, As the reference wind speed, The angle between the wind direction and the line.

[0014] The lightning statistics include annual average ground flash density and lightning current amplitude. The calculation of transmission line tripping rate based on the terrain correction coefficient corresponding to the terrain parameters, the conductor spatial location, the dynamic wind deflection correction coefficient, and the lightning statistics specifically includes:

[0015] The projected area of ​​the exposed arc in the electrical geometry model is determined based on the spatial position of the conductor and the equivalent radius of the exposed arc in the rotating coordinate system.

[0016] The probability distribution function of the lightning current amplitude is determined using the maximum likelihood estimation method.

[0017] The transmission line tripping rate is determined based on the terrain correction coefficient, exposed arc projection area, dynamic wind deflection correction coefficient, probability distribution function of lightning current amplitude, and annual average ground flash density corresponding to the terrain parameters. The transmission line tripping rate includes the backflashover tripping rate.

[0018] Specifically, determining the projected area of ​​the exposed arc in the electrical geometry model based on the spatial position of the conductor and the equivalent radius of the exposed arc in the rotating coordinate system includes:

[0019] according to Determine the projected area of ​​the exposed arc in the electrical geometry model, where, The starting angle of the exposed arc. The termination angle of the exposed arc, Let r0 be the projected area of ​​the exposed arc in the electrical geometry model, and r0 be the equivalent radius of the conductor in the rotating coordinate system. ,t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, h g (I) is the leading development level, K curve This is the sag correction factor.

[0020] Specifically, determining the probability distribution function of the lightning current amplitude using the maximum likelihood estimation method includes:

[0021] use Determine the probability distribution function of the lightning current amplitude, where, Let be the probability distribution function of the lightning current amplitude, α be the overall amplitude level control parameter, β be the lightning current amplitude distribution concentration control parameter, and I be the lightning current amplitude.

[0022] The method of determining the transmission line tripping rate based on the terrain correction coefficient, exposed arc projected area, dynamic wind deflection correction coefficient, probability distribution function of lightning current amplitude, and annual average ground flash density corresponding to the terrain parameters, wherein the transmission line tripping rate includes the backflashover tripping rate, specifically including:

[0023] according to Determine the trip rate due to the circuit breaker, where I max I is the maximum lightning current amplitude. min This is the minimum lightning current amplitude. For the tripping rate due to the circuit breaker, This represents the annual average ground flash density. For the projected area of ​​the exposed arc, This is the dynamic wind deflection correction factor. is the terrain correction coefficient, and P(I) is the corrected probability distribution function of lightning current amplitude.

[0024] Terrain correction factor The slope calculation is expressed as follows:

[0025]

[0026] J represents the slope angle.

[0027] Specifically, determining the lightning risk level of a transmission line based on its tripping rate includes:

[0028] The weighting coefficient for the tripping rate of the transmission line is determined based on the current voltage level of the transmission line.

[0029] The overall tripping rate is determined by combining the weighting coefficients and the transmission line tripping rate.

[0030] The lightning risk level of the transmission line is determined based on the matching of the comprehensive tripping rate with the preset risk level range.

[0031] Specifically, determining the weighting coefficient for the tripping rate of the transmission line based on its current voltage level includes:

[0032] according to The weighting coefficients for determining the transmission line tripping rate are determined, wherein the transmission line tripping rate includes the backlash tripping rate, and... is the weighting coefficient for the tripping rate caused by the inrush, b0 is the benchmark weight for the tripping rate caused by the inrush, and V is the line voltage level.

[0033] Specifically, determining the overall tripping rate by combining the weighting coefficients and the transmission line tripping rate includes:

[0034] according to Determine the overall trip rate, where, The overall tripping rate is given by b(V), where b(V) is the weighting coefficient for the tripping rate caused by the stunt. For the tripping rate due to the circuit breaker, For environmental weighting factors, f(X) env () is a meteorological environment correction factor.

[0035] A transmission line lightning strike risk assessment system based on multi-source data, the system comprising:

[0036] The data acquisition module acquires three-dimensional laser point cloud data, lightning statistics data, and meteorological data of the transmission line. The three-dimensional laser point cloud data includes terrain parameters of the tower's location.

[0037] The wind deflection correction module calculates the spatial position of the conductors in the electrical geometric model of the transmission line under the influence of dynamic wind deflection based on meteorological data, as well as the dynamic wind deflection correction coefficient of the electrical geometric model.

[0038] The transmission line tripping rate determination module calculates the transmission line tripping rate based on the terrain correction coefficient corresponding to the terrain parameters, the spatial position of the conductor, the dynamic wind deflection correction coefficient, and the lightning statistics data.

[0039] The risk level determination module determines the lightning risk level of the transmission line based on the transmission line tripping rate.

[0040] A computer-readable storage medium storing a computer program that, when executed by a processor, causes the processor to perform the steps of the method described above.

[0041] The embodiments of the present invention have the following beneficial effects:

[0042] This invention provides a detailed assessment of the lightning strike risk of transmission lines by combining 3D laser point cloud data, lightning statistics, and meteorological data. The 3D laser point cloud data provides detailed information about the terrain where the towers are located, facilitating accurate terrain correction. Based on meteorological data, the method can dynamically calculate the spatial position of the conductors under wind deflection and its correction coefficient, improving the accuracy of the electrical geometry model. By combining terrain correction coefficients, conductor spatial position, wind deflection correction coefficients, and lightning statistics, the tripping rate of the transmission line, especially the backflashover tripping rate, can be accurately calculated, thereby accurately assessing the line risk. Attached Figure Description

[0043] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0044] in:

[0045] Figure 1 This is a flowchart illustrating an embodiment of a method for assessing the lightning strike risk of transmission lines based on multi-source data provided by the present invention.

[0046] Figure 2 This is a flowchart illustrating another embodiment of the transmission line lightning strike risk assessment method based on multi-source data provided by the present invention.

[0047] Figure 3 This is a schematic diagram of an embodiment of a transmission line lightning strike risk assessment system based on multi-source data provided by the present invention;

[0048] Figure 4 A schematic diagram of the structure of an embodiment of the medium provided by the present invention. Detailed Implementation

[0049] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0050] like Figure 1 As shown, Figure 1 This is a flowchart illustrating an embodiment of a method for assessing lightning strike risk of transmission lines based on multi-source data provided by the present invention. The method includes:

[0051] S101: Acquire 3D laser point cloud data, lightning statistics and meteorological data of the transmission line. The 3D laser point cloud data includes terrain parameters of the tower location.

[0052] In one implementation scenario, three-dimensional laser point cloud data, lightning statistics, and meteorological data of the transmission line are acquired. The three-dimensional laser point cloud data includes the geometric parameters of the tower and the terrain parameters of the tower's location.

[0053] S102: Based on meteorological data, calculate the spatial position of the conductor in the electrical geometric model of the transmission line under the influence of dynamic wind deflection, and the dynamic wind deflection correction coefficient of the electrical geometric model.

[0054] In one implementation scenario, meteorological data includes wind deflection amplitude and wind speed. The spatial position of the conductors in the electrical geometry model of the transmission line under the influence of dynamic wind deflection is calculated using the formula shown below:

[0055] ;

[0056] Where y( (t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, t is time, the spatial position of the conductor changes with time, H is the horizontal tension, w is the weight per unit length, A is the wind deflection amplitude, and x is the horizontal tension. Let be the projected coordinates of the conductor on the horizontal plane.

[0057] The dynamic wind deflection correction factor for the electrical geometry model is calculated using the formula shown below:

[0058]

[0059] in, This represents the dynamic wind deflection correction coefficient for the electrical geometry model. For real-time wind speed, As the reference wind speed, The angle between the wind direction and the line.

[0060] S103: Calculate the transmission line tripping rate based on the terrain correction coefficient corresponding to the terrain parameters, the spatial position of the conductor, the dynamic wind deflection correction coefficient, and lightning statistics.

[0061] In one implementation scenario, the transmission line tripping rate includes the backflashover tripping rate. Lightning statistics include the annual average ground flash density and lightning current amplitude. Based on the conductor's spatial location and the equivalent radius of the exposed arc in the rotating coordinate system, the projected area of ​​the exposed arc in the electrical geometry model is determined. The probability distribution function of the lightning current amplitude is determined using the maximum likelihood estimation method. The transmission line tripping rate is determined based on the terrain correction coefficients corresponding to the terrain parameters, the projected area of ​​the exposed arc, the dynamic wind deflection correction coefficient, the probability distribution function of the lightning current amplitude, and the annual average ground flash density. The transmission line tripping rate includes the backflashover tripping rate, as shown in the following formula:

[0062] ;

[0063] Among them, I max I is the maximum lightning current amplitude. min This is the minimum lightning current amplitude. For the tripping rate due to the circuit breaker, This represents the annual average ground flash density. For the projected area of ​​the exposed arc, This is the dynamic wind deflection correction factor. Let P(I) be the terrain correction coefficient, and P(I) be the probability distribution function of the corrected lightning current amplitude, where I is the lightning current amplitude. Terrain Correction Coefficient The slope calculation is expressed as follows:

[0064]

[0065] J represents the slope angle.

[0066] In another implementation scenario, the transmission line tripping rate includes the backflashover tripping rate and the backflashover tripping rate. The spatial position of the conductors in the electrical geometry model of the transmission line under the influence of dynamic wind deflection can be calculated based on meteorological data, along with the dynamic wind deflection correction coefficient for the electrical geometry model. The backflashover tripping rate can be calculated based on the conductor spatial position, the dynamic wind deflection correction coefficient, and lightning statistics. Furthermore, the tower attraction range can be corrected based on the terrain correction coefficients corresponding to the geometric and terrain parameters. The backflashover tripping rate can then be calculated based on the corrected tower attraction range and lightning statistics.

[0067] Specifically, when calculating the backflashover trip rate, the tower attraction range is corrected based on the terrain correction coefficients corresponding to the impulse grounding resistance, geometric parameters, and terrain parameters to determine the corrected tower attraction range. The tower attraction range is determined according to the following formula:

[0068] ;

[0069] in, For the attraction range of the tower, The tower height is a geometric parameter. To reduce the grounding resistance, This is the terrain correction factor.

[0070] Furthermore, lightning overvoltage simulations were performed on the transmission line to determine the critical lightning current amplitude; the probability of backflashover critical lightning current when the lightning current amplitude exceeds the critical lightning current amplitude was also determined. The probability of backflashover critical lightning current is shown in the following formula:

[0071] ;

[0072] Among them, I crit This is the critical lightning current amplitude. To counteract the critical lightning current probability of flashover.

[0073] Furthermore, based on the annual average ground flashover density, the critical lightning current probability of backflashover, and the tower attraction range, the backflashover trip rate is determined. The backflashover trip rate is determined according to the following formula:

[0074] ;

[0075] in, To counter the tripping rate, The annual average lightning density (times / km²·year) is given. To account for the tower attraction range of micro-topography correction, P(I>I) crit ) represents the critical lightning current probability for counter-flashover.

[0076] When calculating the backflashover trip rate, meteorological data for the transmission line is first obtained, including wind amplitude and wind speed. Based on the meteorological data, the spatial position of the conductors in the electrical geometry model of the transmission line under the influence of dynamic wind deflection is calculated, as well as the dynamic wind deflection correction coefficient of the electrical geometry model. The backflashover trip rate is then calculated based on the terrain correction coefficient corresponding to the terrain parameters, the conductor spatial position, the dynamic wind deflection correction coefficient, and lightning statistics.

[0077] S105: Determine the lightning risk level of a transmission line based on its tripping rate.

[0078] In one implementation scenario, the weighting coefficient for the transmission line tripping rate is determined based on the current voltage level of the transmission line. Specifically, the voltage level adjustment range for the transmission line is shown in the table below:

[0079]

[0080] Furthermore, the comprehensive tripping rate is determined by combining weighting coefficients and transmission line tripping rates; the lightning risk level of the transmission line is determined based on the matching between the comprehensive tripping rate and the preset risk level range. Specifically, the preset risk level range is constructed using the comprehensive tripping rate, in units of trips / (100km·year).

[0081] The first preset risk level range is [0, 0.05).

[0082] The second preset risk level range is [0.05, 0.15).

[0083] The third preset risk level range is: [0.15, 0.30).

[0084] The fourth preset risk level range is [0.30, +∞).

[0085] If the overall tripping rate is within the first preset risk level range, the lightning strike risk level of the transmission line is Level I (low risk); if the overall tripping rate is within the second preset risk level range, the lightning strike risk level of the transmission line is Level II (moderate risk); if the overall tripping rate is within the third preset risk level range, the lightning strike risk level of the transmission line is Level III (higher risk); if the overall tripping rate is within the fourth preset risk level range, the lightning strike risk level of the transmission line is Level IV (high risk).

[0086] For areas with severe lightning activity (more than 90 thunderstorm days per year), the threshold range between Level III and Level IV is narrowed by 50%.

[0087] For important power transmission corridors (lines that fail the N-1 check), the risk level is automatically upgraded by 1 level.

[0088] Furthermore, based on a pre-defined rule base for matching measures, differentiated lightning protection upgrade strategies corresponding to risk levels I-IV are determined. The rule base for matching measures is shown below:

[0089] The generation of lightning protection upgrade recommendations includes a rule base for matching measures:

[0090]

[0091] As described above, this invention comprehensively assesses the lightning strike risk of transmission lines by combining three-dimensional laser point cloud data, lightning statistics, and meteorological data. The three-dimensional laser point cloud data provides detailed information about the terrain where the towers are located, facilitating accurate terrain correction. Based on meteorological data, the method can dynamically calculate the spatial position of the conductor under wind deflection and its correction coefficient, improving the accuracy of the electrical geometry model. By combining terrain correction coefficients, conductor spatial position, wind deflection correction coefficients, and lightning statistics, the tripping rate of the transmission line, especially the backflashover tripping rate, can be accurately calculated, thereby accurately assessing the line risk.

[0092] like Figure 2 As shown, Figure 2This is a flowchart illustrating another embodiment of the transmission line lightning strike risk assessment method based on multi-source data provided by the present invention. The method includes:

[0093] S201: Acquire 3D laser point cloud data, lightning statistics and meteorological data of transmission lines. The 3D laser point cloud data includes terrain parameters of the tower location.

[0094] S202: Based on meteorological data, calculate the spatial position of the conductor in the electrical geometric model of the transmission line under the influence of dynamic wind deflection, and the dynamic wind deflection correction coefficient of the electrical geometric model.

[0095] In one implementation scenario, meteorological data includes wind deflection amplitude and wind speed. The spatial position of the conductors in the electrical geometry model of the transmission line under the influence of dynamic wind deflection is calculated using the formula shown below:

[0096] ;

[0097] Where y( (t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, t is time, H is horizontal tension, w is weight per unit length, A is wind deflection amplitude, and x is the horizontal tension. The projected coordinates of the conductor on the horizontal plane:

[0098] ;

[0099] Where R(I) is the leader attraction radius corresponding to the lightning current amplitude I. Let K be the polar coordinate angle centered on the tower. proj ( Terrain projection distortion correction coefficient.

[0100] The dynamic wind deflection correction factor for the electrical geometry model is calculated using the formula shown below:

[0101] ;

[0102] in, This represents the dynamic wind deflection correction coefficient for the electrical geometry model. For real-time wind speed, As the reference wind speed, The angle between the wind direction and the line.

[0103] S203: Determine the projected area of ​​the exposed arc in the electrical geometry model based on the spatial position of the conductor and the equivalent radius of the exposed arc in the rotating coordinate system.

[0104] In one implementation scenario, the projected area of ​​the exposed arc in the electrical geometry model is determined based on the spatial location of the conductor and the equivalent radius of the exposed arc in the rotating coordinate system, as shown in the following formula:

[0105] ;

[0106] in, The starting angle of the exposed arc. The termination angle of the exposed arc, Let r0 be the projected area of ​​the exposed arc in the electrical geometry model, and r0 be the equivalent radius of the conductor in the rotating coordinate system. ,t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, h g (I) is the leading development level, K curve This is the sag correction factor.

[0107] S204: Determine the probability distribution function of lightning current amplitude using the maximum likelihood estimation method.

[0108] In one implementation scenario, lightning statistics include the annual average ground flash density and lightning current amplitude. The probability distribution function of the lightning current amplitude is determined using the maximum likelihood estimation method. The corrected probability distribution function is shown in the following equation:

[0109] ;

[0110] in, Let α be the probability distribution function of the lightning current amplitude, β be the overall amplitude level control parameter, β be the lightning current amplitude distribution concentration control parameter, and I be the lightning current amplitude. Parameters α and β are updated quarterly.

[0111] S205: The transmission line tripping rate is determined based on the terrain correction coefficient, exposed arc projection area, dynamic wind deflection correction coefficient, probability distribution function of lightning current amplitude, and annual average ground flash density corresponding to the terrain parameters. The transmission line tripping rate includes the backflashover tripping rate.

[0112] In one implementation scenario, the bypass trip rate is determined based on the terrain correction coefficient, the projected area of ​​the exposed arc, the dynamic wind deflection correction coefficient, the probability distribution function of the lightning current amplitude, and the annual average ground flash density, as shown in the following formula:

[0113] ;

[0114] Among them, I max I is the maximum lightning current amplitude. min This is the minimum lightning current amplitude. For the tripping rate due to the circuit breaker, This represents the annual average ground flash density. For the projected area of ​​the exposed arc, This is the dynamic wind deflection correction factor. Let P(I) be the terrain correction factor, and P(I) be the corrected probability distribution function of the lightning current amplitude. Terrain Correction Factor The slope calculation is expressed as follows:

[0115]

[0116] J represents the slope angle.

[0117] In another implementation scenario, the transmission line tripping rate includes the backflashover tripping rate and the backflashover tripping rate. When calculating the backflashover tripping rate, it is determined based on the terrain correction coefficient, the projected area of ​​the exposed arc, the dynamic wind deflection correction coefficient, the probability distribution function of the lightning current amplitude, and the annual average ground flash density.

[0118] When calculating the backflashover trip rate, the tower attraction range is corrected according to the terrain correction coefficients corresponding to the geometric parameters and terrain parameters. The transmission line trip rate is determined based on the corrected tower attraction range and lightning statistics. The transmission line trip rate includes the backflashover trip rate.

[0119] Specifically, the tower attraction range is corrected based on the terrain correction coefficients corresponding to the impulse grounding resistance, geometric parameters, and terrain parameters to determine the corrected tower attraction range. The geometric parameters include tower height, conductor heights, phase-to-phase distances, ground wire height, and insulator string length; the terrain parameters include slope, aspect, and relative elevation. The tower attraction range is determined using the following formula:

[0120] ;

[0121] in, For the attraction range of the tower, The tower height is a geometric parameter. To reduce the grounding resistance, This is the terrain correction factor.

[0122] Furthermore, lightning overvoltage simulations were performed on the transmission line to determine the critical lightning current amplitude. The probability of a backflashover critical lightning current when the lightning current amplitude exceeds the critical lightning current amplitude was determined, as shown in the following formula:

[0123] ;

[0124] Among them, I crit This is the critical lightning current amplitude. To counteract the critical lightning current probability of flashover.

[0125] Furthermore, based on the annual average ground flashover density, the critical lightning current probability of backflashover, and the tower attraction range, the transmission line tripping rate is determined. The transmission line tripping rate includes the backflashover tripping rate, as shown in the following formula:

[0126] ;

[0127] in, To counter the tripping rate, The annual average lightning density (times / km²·year) is given. To account for the tower attraction range of micro-topography correction, P(I>I) crit ) represents the critical lightning current probability for counter-flashover.

[0128] S206: Determine the weighting coefficient of the transmission line tripping rate based on the current voltage level of the transmission line.

[0129] In one implementation scenario, the transmission line tripping rate includes the backlash tripping rate, and the weighting coefficient of the transmission line tripping rate is determined according to the following formula:

[0130]

[0131] in, is the weighting coefficient for the tripping rate caused by the inrush, b0 is the benchmark weight for the tripping rate caused by the inrush, and V is the line voltage level.

[0132] In another implementation scenario, the transmission line tripping rate includes backflashover tripping rate and backflashover tripping rate.

[0133] The weighting coefficient for the back-off trip rate is determined according to the following formula:

[0134] ;

[0135] in, , where a0 is the weighting coefficient for the backflashover trip rate, V is the base weight for the backflashover trip rate, and V is the line voltage level.

[0136] The weighting coefficient for the tripping rate due to shunting is determined according to the following formula:

[0137]

[0138] in, is the weighting coefficient for the tripping rate caused by the inrush, b0 is the benchmark weight for the tripping rate caused by the inrush, and V is the line voltage level.

[0139] S207: Determine the overall tripping rate by combining the weighting coefficient and the transmission line tripping rate.

[0140] In one implementation scenario, the transmission line tripping rate includes the tripping rate caused by backlash, and the overall tripping rate is determined according to the following formula:

[0141] ;

[0142] in, The overall tripping rate is given by b(V), where b(V) is the weighting coefficient for the tripping rate caused by the stunt. For the tripping rate due to the circuit breaker, For environmental weighting factors, f(X) env () is a meteorological environment correction factor.

[0143] In another implementation scenario, the transmission line tripping rate includes backflashover tripping rate and backflashover tripping rate.

[0144] The overall trip rate is determined according to the following formula:

[0145] ;

[0146] Where a(V) is the weighting coefficient for the backflashover tripping rate, and b(V) is the weighting coefficient for the bypass tripping rate. For environmental weighting factors, f(X) env) It is a meteorological environment correction factor.

[0147] S208: Determine the lightning risk level of the transmission line based on the matching of the comprehensive tripping rate and the preset risk level range.

[0148] It should be noted that step S208 is in Figure 1 The implementation scenarios shown have been discussed in detail and will not be repeated here.

[0149] As described above, this invention extracts tower geometry and micro-topographical features from collected 3D laser point cloud data, and accurately calculates backflashover and bypass tripping rates by combining lightning statistical distribution data and meteorological data. A dynamic weighting algorithm is used to comprehensively assess tripping risk and match it with a preset risk level threshold matrix, outputting differentiated lightning protection upgrade strategies. By dynamically adjusting the weighting coefficients of backflashover and bypass tripping rates to adapt to the characteristics of different voltage levels and matching them with a preset risk level threshold matrix, accurate risk assessment and standardized upgrades are achieved. Overall, this invention effectively integrates multi-source data, improves the accuracy and economy of transmission line lightning risk assessment, and significantly enhances the safe and stable operation capability of the power grid.

[0150] like Figure 3 As shown, Figure 3 This is a schematic diagram of an embodiment of a transmission line lightning strike risk assessment system based on multi-source data provided by the present invention. A transmission line lightning strike risk assessment system 10 based on multi-source data includes:

[0151] The data acquisition module 11 is used to acquire three-dimensional laser point cloud data and lightning statistics data of the transmission line. The three-dimensional laser point cloud data includes the geometric parameters of the tower and the terrain parameters of the tower's location.

[0152] In one implementation scenario, the data acquisition module 11 acquires three-dimensional laser point cloud data and lightning statistics of the transmission line. The three-dimensional laser point cloud data includes the geometric parameters of the tower and the terrain parameters of the tower's location.

[0153] The wind deflection correction module 12 is used to calculate the spatial position of the conductor in the electrical geometric model of the transmission line under the influence of dynamic wind deflection based on meteorological data, as well as the dynamic wind deflection correction coefficient of the electrical geometric model.

[0154] In one implementation scenario, in the wind deflection correction module 12, the spatial position of the conductor in the electrical geometric model of the transmission line under the influence of dynamic wind deflection is calculated based on meteorological data, as well as the dynamic wind deflection correction coefficient of the electrical geometric model.

[0155] The transmission line tripping rate determination module 13 is used to calculate the transmission line tripping rate based on the terrain correction coefficient corresponding to the terrain parameters, the spatial position of the conductor, the dynamic wind deflection correction coefficient, and lightning statistics.

[0156] In one implementation scenario, in the transmission line tripping rate determination module 13, lightning statistics include the annual average ground flash density and lightning current amplitude. Based on the spatial location of the conductor and the equivalent radius of the exposed arc in the rotating coordinate system, the projected area of ​​the exposed arc in the electrical geometric model is determined. The probability distribution function of the lightning current amplitude is determined using the maximum likelihood estimation method. The transmission line tripping rate is determined based on the terrain correction coefficient corresponding to the terrain parameters, the projected area of ​​the exposed arc, the dynamic wind deflection correction coefficient, the probability distribution function of the lightning current amplitude, and the annual average ground flash density. The transmission line tripping rate includes the backflashover tripping rate.

[0157] In another implementation scenario, the transmission line tripping rate includes both backflashover tripping rate and backflashover tripping rate. The tower attraction range can be corrected based on the terrain correction coefficients corresponding to the impulse grounding resistance, geometric parameters, and terrain parameters, thus determining the corrected tower attraction range. Lightning overvoltage simulation of the transmission line is performed to determine the critical lightning current amplitude. The probability of backflashover critical lightning current when the lightning current amplitude exceeds the critical lightning current amplitude is determined. Based on the annual average ground flashover density, the probability of backflashover critical lightning current, and the tower attraction range, the transmission line tripping rate is determined, including the backflashover tripping rate. Furthermore, the conductor spatial position of the transmission line's electrical geometric model under the influence of dynamic wind deflection can be calculated based on meteorological data, along with the dynamic wind deflection correction coefficient for the electrical geometric model. The backflashover tripping rate is then calculated based on the conductor spatial position, the dynamic wind deflection correction coefficient, and lightning statistics.

[0158] The risk level determination module 14 is used to determine the lightning risk level of a transmission line based on the transmission line tripping rate.

[0159] In one implementation scenario, the weighting coefficient of the transmission line tripping rate is determined based on the current voltage level of the transmission line; the comprehensive tripping rate is determined by combining the weighting coefficient and the transmission line tripping rate; and the lightning risk level of the transmission line is determined based on the matching of the comprehensive tripping rate with the preset risk level range.

[0160] like Figure 4 As shown, Figure 4 This is a schematic diagram of the structure of an embodiment of the medium provided by the present invention. The medium 20 stores at least one computer program 21, which is executed by a processor to perform the following... Figure 1 and Figure 2 The method shown is detailed above and will not be repeated here. In one embodiment, the medium 20 can be a storage chip, hard disk, portable hard disk, USB flash drive, optical disk, or other read / write storage device, or even a server, etc.

[0161] Furthermore, the processes depicted in the accompanying drawings do not necessarily have to be performed in the specific or sequential order shown to achieve the desired result. In some implementations, multitasking and parallel processing are possible or may be advantageous.

[0162] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, the embodiments of apparatus, devices, and non-volatile computer-readable storage media are basically similar to the method embodiments, and therefore described more simply; relevant parts can be referred to the descriptions of the method embodiments.

[0163] The apparatus, device, non-volatile computer-readable storage medium and method provided in the embodiments of this specification are corresponding. Therefore, the apparatus, device and non-volatile computer storage medium also have similar beneficial technical effects as the corresponding method. Since the beneficial technical effects of the method have been described in detail above, the beneficial technical effects of the corresponding apparatus, device and non-volatile computer storage medium will not be repeated here.

[0164] The systems, devices, modules, or units described in the above embodiments can be implemented by computer chips or entities, or by products with certain functions. A typical implementation device is a computer. Specifically, a computer can be, for example, a personal computer, laptop computer, cellular phone, camera phone, smartphone, personal digital assistant, media player, navigation device, email device, game console, tablet computer, wearable device, or any combination of these devices.

[0165] For ease of description, the above apparatus is described by dividing it into various functional units. Of course, in implementing this specification, the functions of each unit can be implemented in one or more software and / or hardware components. Those skilled in the art will understand that the embodiments of this specification can be provided as methods, systems, or computer program products. Therefore, the embodiments of this specification can take the form of entirely hardware embodiments, entirely software embodiments, or embodiments combining software and hardware aspects. Furthermore, the embodiments of this specification can take the form of computer program products 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.

[0166] This specification is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this specification. 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 processor, 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 a machine for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0167] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0168] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0169] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0170] Memory may include non-persistent storage in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0171] Computer-readable media includes both permanent and non-permanent, removable and non-removable media that can store information by any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.

[0172] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0173] This specification can be described in the general context of computer-executable instructions that are executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform a specific task or implement a specific abstract data type. This specification can also be practiced in distributed computing environments, where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0174] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, the system embodiments are basically similar to the method embodiments, so the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments.

[0175] The above description discloses only preferred embodiments of the present invention and should not be construed as limiting the scope of the present invention. Therefore, equivalent variations made in accordance with the claims of the present invention are still within the scope of the present invention.

Claims

1. A method for assessing the lightning strike risk of transmission lines based on multi-source data, characterized in that, The method includes: Acquire three-dimensional laser point cloud data, lightning statistics data, and meteorological data of the transmission line, wherein the three-dimensional laser point cloud data includes terrain parameters of the tower location; The calculation of the conductor spatial position of the transmission line's electrical geometric model under the influence of dynamic wind deflection, based on meteorological data, and the dynamic wind deflection correction coefficient of the electrical geometric model, specifically includes: according to Calculate the spatial position of the conductors in the electrical geometry model of the transmission line under the influence of dynamic wind deflection, where t is time and y( ,t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, H represents the horizontal tension, w represents the weight per unit length, A represents the wind deflection amplitude, and x represents the horizontal tension. Let be the projected coordinates of the conductor on the horizontal plane; according to Calculate the dynamic wind deflection correction coefficient for the electrical geometry model, where, This represents the dynamic wind deflection correction coefficient for the electrical geometry model. For real-time wind speed, As the reference wind speed, The angle between the wind direction and the line; The transmission line tripping rate is calculated based on the terrain correction coefficient corresponding to the terrain parameters, the spatial position of the conductor, the dynamic wind deflection correction coefficient, and the lightning statistics data. Determining the lightning risk level of a transmission line based on its tripping rate, specifically includes: determining a weighting coefficient for the tripping rate based on the current voltage level of the transmission line; determining a comprehensive tripping rate by combining the weighting coefficient and the tripping rate; and determining the lightning risk level of the transmission line based on the matching between the comprehensive tripping rate and a preset risk level range. The determination of the weighting coefficient for the tripping rate based on the current voltage level of the transmission line specifically includes: according to The weighting coefficients for determining the transmission line tripping rate are determined, wherein the transmission line tripping rate includes the backlash tripping rate, and... is the weighting coefficient for the tripping rate caused by the inrush, b0 is the benchmark weight for the tripping rate caused by the inrush, and V is the line voltage level.

2. The method for assessing the lightning strike risk of transmission lines based on multi-source data according to claim 1, characterized in that, The lightning statistics include the annual average ground flash density and lightning current amplitude. The calculation of the transmission line tripping rate based on the terrain correction coefficient corresponding to the terrain parameters, the conductor spatial location, the dynamic wind deflection correction coefficient, and the lightning statistics specifically includes: Based on the spatial position of the conductor and the equivalent radius of the exposed arc in the rotating coordinate system, determine the projected area of ​​the exposed arc in the electrical geometry model; The probability distribution function of the lightning current amplitude is determined using the maximum likelihood estimation method. The transmission line tripping rate is determined based on the terrain correction coefficient, exposed arc projection area, dynamic wind deflection correction coefficient, probability distribution function of lightning current amplitude, and annual average ground flash density corresponding to the terrain parameters. The transmission line tripping rate includes the backflashover tripping rate.

3. The method for assessing the lightning strike risk of transmission lines based on multi-source data according to claim 2, characterized in that, The step of determining the projected area of ​​the exposed arc in the electrical geometry model based on the spatial position of the conductor and the equivalent radius of the exposed arc in the rotating coordinate system specifically includes: according to Determine the projected area of ​​the exposed arc in the electrical geometry model, where, The starting angle of the exposed arc. The termination angle of the exposed arc, Let r0 be the projected area of ​​the exposed arc in the electrical geometry model, and r0 be the equivalent radius of the conductor in the rotating coordinate system. ,t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, h g (I) is the leading development level, K curve This is the sag correction factor.

4. The method for assessing the lightning strike risk of transmission lines based on multi-source data according to claim 2, characterized in that, The method of determining the probability distribution function of lightning current amplitude using the maximum likelihood estimation method specifically includes: use Determine the probability distribution function of the lightning current amplitude, where, Let be the probability distribution function of the lightning current amplitude, α be the overall amplitude level control parameter, β be the lightning current amplitude distribution concentration control parameter, and I be the lightning current amplitude.

5. The method for assessing the lightning strike risk of transmission lines based on multi-source data according to claim 3 or 4, characterized in that, The transmission line tripping rate is determined based on the terrain correction coefficient, exposed arc projected area, dynamic wind deflection correction coefficient, probability distribution function of lightning current amplitude, and annual average ground flash density corresponding to the terrain parameters. The transmission line tripping rate includes the backflashover tripping rate, specifically including: according to Determine the trip rate due to the circuit breaker, where I max I is the maximum lightning current amplitude. min This is the minimum lightning current amplitude. For the bypass trip rate, This represents the annual average ground flash density. For the projected area of ​​the exposed arc, This is the dynamic wind deflection correction factor. Let P(I) be the terrain correction coefficient, and P(I) be the probability distribution function of the lightning current amplitude, where the terrain correction coefficient is... The slope calculation formula is as follows: J represents the slope angle.

6. The method for assessing the lightning strike risk of transmission lines based on multi-source data according to claim 1, characterized in that, The determination of the overall tripping rate by combining the weighting coefficients and the transmission line tripping rate specifically includes: according to Determine the overall trip rate, where, The overall tripping rate is given by b(V), where b(V) is the weighting coefficient for the tripping rate caused by the stunt. For the bypass trip rate, For environmental weighting factors, f(X) env () is a meteorological environment correction factor.

7. A transmission line lightning strike risk assessment system based on multi-source data, characterized in that, The system includes: The data acquisition module is used to acquire three-dimensional laser point cloud data, lightning statistics data and meteorological data of the transmission line. The three-dimensional laser point cloud data includes the terrain parameters of the tower location. The wind deflection correction module is used to calculate the spatial position of the conductors in the electrical geometric model of the transmission line under the influence of dynamic wind deflection, and the dynamic wind deflection correction coefficient of the electrical geometric model, based on meteorological data. The meteorological data includes wind deflection amplitude and wind speed. The calculation of the spatial position of the conductors in the electrical geometric model of the transmission line under the influence of dynamic wind deflection, and the dynamic wind deflection correction coefficient of the electrical geometric model, based on meteorological data, specifically includes: according to Calculate the spatial position of the conductors in the electrical geometry model of the transmission line under the influence of dynamic wind deflection, where t is time and y( ,t) represents the spatial position of the conductor in the electrical geometry model under dynamic wind deflection, H represents the horizontal tension, w represents the weight per unit length, A represents the wind deflection amplitude, and x represents the horizontal tension. Let be the projected coordinates of the conductor on the horizontal plane; according to Calculate the dynamic wind deflection correction coefficient for the electrical geometry model, where, This represents the dynamic wind deflection correction coefficient for the electrical geometry model. For real-time wind speed, As the reference wind speed, The angle between the wind direction and the line; The transmission line tripping rate determination module is used to calculate the transmission line tripping rate based on the terrain correction coefficient corresponding to the terrain parameters, the spatial position of the conductor, the dynamic wind deflection correction coefficient, and the lightning statistics data. The risk level determination module is used to determine the lightning risk level of a transmission line based on its tripping rate. Specifically, determining the lightning risk level based on the tripping rate includes: determining a weighting coefficient for the tripping rate based on the current voltage level of the transmission line; determining a comprehensive tripping rate by combining the weighting coefficient and the tripping rate; and determining the lightning risk level of the transmission line based on the matching between the comprehensive tripping rate and a preset risk level range. The determination of the weighting coefficient for the tripping rate based on the current voltage level of the transmission line specifically includes: according to The weighting coefficients for determining the transmission line tripping rate are determined, wherein the transmission line tripping rate includes the backlash tripping rate, and... is the weighting coefficient for the tripping rate caused by the inrush, b0 is the benchmark weight for the tripping rate caused by the inrush, and V is the line voltage level.

8. A computer-readable storage medium, characterized in that, The system stores a computer program that, when executed by a processor, causes the processor to perform the steps of the method as described in any one of claims 1 to 6.

Images