CN112541211A - Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel - Google Patents

Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel Download PDF

Info

Publication number
CN112541211A
CN112541211A CN202011380660.9A CN202011380660A CN112541211A CN 112541211 A CN112541211 A CN 112541211A CN 202011380660 A CN202011380660 A CN 202011380660A CN 112541211 A CN112541211 A CN 112541211A
Authority
CN
China
Prior art keywords
lightning
current
distance
scanning
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011380660.9A
Other languages
Chinese (zh)
Other versions
CN112541211B (en
Inventor
邓冶强
王羽
高磊
陈小月
文习山
蓝磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202011380660.9A priority Critical patent/CN112541211B/en
Publication of CN112541211A publication Critical patent/CN112541211A/en
Application granted granted Critical
Publication of CN112541211B publication Critical patent/CN112541211B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Structural Engineering (AREA)
  • Locating Faults (AREA)

Abstract

The invention relates to the field of lightning hazard protection of a power transmission line, in particular to a lightning trip-out rate scanning calculation method considering shielding of a plurality of lines in a power transmission channel. By adopting a new computer scanning method, based on the translation of the lightning leading linear equation, the specific exposed arc section of the strike distance circle of each wire of the power transmission line does not need to be solved, the lightning triggering width of the projection of the corresponding exposed arc section on the ground is directly obtained, and the processing of the position relation of different strike distance circles is skipped, so that the method is suitable for the complicated multi-circuit tower type shielding strike calculation.

Description

Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel
Technical Field
The invention belongs to the field of lightning hazard protection of power transmission lines, and particularly relates to a lightning trip-out rate scanning calculation method considering shielding of multiple lines in a power transmission channel.
Background
With the rapid construction of ultra-high and extra-high voltage lines in China and the shortage of land resources, available line corridors are reduced day by day, so that some power transmission channels with large transmission capacity and compact line arrangement gradually appear. The channels are mostly composed of two or more than two loops of important power transmission lines, the distance between the lines in the power transmission channels is small, and when the lightning damage risk of the dense power transmission channel lines is evaluated by the existing line lightning protection evaluation system, the influence of the shielding effect of the power transmission lines in the same channel is not considered, so that the trip-out rate of the evaluated dense channel lines is possibly high.
At present, the lightning shielding effect of a power transmission line is mainly calculated by adopting an electrical geometric model method (EGM) in engineering, and the method mainly has the idea that a stroke distance circle of a lead and a ground wire is drawn according to the change of a lightning current amplitude and a pilot incidence angle, and an exposure arc and the projection distance of the exposure arc on the ground are calculated according to the position relation (intersection, phase separation and tangency) of the stroke distance circle. However, in the classic EGM algorithm, the position relationship of the impact distance circle is easy to process because only the condition of a single-base-tower transmission line is considered; under the condition that a plurality of base tower lines of dense transmission channels are arranged in parallel, the number of the wire and the ground wire impact distance circle is rapidly increased, the position relation of the circle becomes very complex, repeated or omitted conditions are avoided when the position relation of the impact distance circle is discussed in a classification mode, and the traditional shielding failure calculation method based on the geometric relation becomes extremely complex.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a lightning trip-out rate scanning calculation method considering shielding of a plurality of lines in a power transmission channel.
In order to solve the technical problems, the invention adopts the following technical scheme: the lightning trip-out rate scanning calculation method considering shielding of a plurality of lines in a power transmission channel comprises the following steps:
step 1, obtaining tower parameters, terrain parameters and lightning activity parameters of left and right lines of a power transmission channel;
step 2, determining a attack distance and attack distance coefficients, numbering all ground wire attack distance circles of the left side tower and the right side tower, and determining a lightning current pilot starting equation and a current lightning current pilot amplitude and a scanning range under the current lightning current pilot amplitude and the current lightning current amplitude according to the current lightning current incident angle and the current lightning current amplitude;
step 3, scanning is carried out by changing the cross-sectional distance of the lightning current pilot equation on the x axis, the intersection point coordinates of the current scanning equation and each stroke distance circle are obtained, wherein the intersection point with the highest vertical coordinate is determined as the current scanning point, and the number of the stroke distance circle to which the current scanning point belongs and whether the stroke distance circle is shielded by the ground are judged;
step 4, when the number of the pitch circle to which the scanning point belongs changes, the scanning of the previous pitch circle under the leading equation is represented to be completed, the lightning triggering width of the previous pitch circle under the current lightning leading equation is obtained, and the lightning triggering widths of all pitch circles are sequentially determined in the scanning process of the lightning leading according to the method;
and 5, circularly considering the lightning triggering widths of all the stroke distance circles under different lightning incident angles and lightning current amplitudes, and calculating the tripping rate of the corresponding lead of the different stroke distance circles based on the lightning leading incident angle distribution function and the lightning current amplitude distribution function.
In the above-mentioned lightning trip-out rate scanning calculation method considering shielding of multiple lines in a power transmission channel, the implementation of step 1 includes:
inquiring a transmission line operation database to obtain left and right line parameters of the transmission line, wherein the left and right line parameters comprise the intermediate distance, the height and the voltage grade of each phase conducting ground wire of a left tower and a right tower, and the ground inclination angles of the left tower and the right tower in a channel, the altitude of the ground where the left tower and the right tower are located and the distance between the left tower and the right tower are obtained according to a GIS system; acquiring the lightning activity condition of a line in a power transmission channel according to the recorded data of the lightning positioning system; or selecting according to the following formula;
lightning amplitude probability distribution, lightning incidence angle probability distribution and striking rod rate:
Figure BDA0002808399390000021
P(φ)＝0.75×cos3(φ)
g＝1:4
wherein I is the lightning current amplitude; phi is the lightning incident angle (-90 degrees); g is the striking rod rate.
In the above-mentioned lightning trip-out rate scanning calculation method considering shielding of multiple lines in a power transmission channel, the implementation of step 2 includes:
step 2.1, determining the hit distance and the hit distance coefficient by adopting the following formula:
rc＝10×I0.65，
rgw＝kgw×rc，
rg＝kg×rc，
wherein r isc,rgw,rgRespectively corresponding to the conducting wire, the lightning conductor and the earth strike distance; k is a radical ofgwAnd kgRespectively a lightning conductor strike distance coefficient and an earth strike distance coefficient, wherein kgw＝1,kg0.9; i is the lightning current amplitude;
step 2.2, numbering all ground wire strike distance circles of the left and right side towers, and determining a lightning current pilot initial equation, a current lightning current pilot amplitude and a scanning range under the incident angle, the strike distance circle number and the scanning range according to the current lightning current incident angle and the current lightning current amplitude;
step 2.2.1, calculating the hitting distances of all hitting distance circles according to the lightning current amplitude, and drawing all hitting distance circles;
step 2.2.2, obtaining the minimum value x of the abscissa of all the stroke distance circlesminOrdinate minimum value yminMaximum value x of abscissamaxMaximum value y of ordinatemaxAnd with (x)min,ymin),(xmin,ymax),(xmax,ymin)，(xmax,ymax) Making four points into a rectangle;
step 2.2.3, when the pilot incident angle is a positive angle, the tangent value of the pilot incident angle is taken as a slope, and the positive angle and the slope respectively pass through the upper left point (x) of the rectanglemin,ymax) And lower right point (x)max,ymin) Making two parallel lines which intersect (a) at the x-axismin0) and (a)max0) two points, wherein the part between the two points is the scanning range under the current lightning current incidence angle and the current lightning current amplitude;
step 2.2.4, when the pilot incident angle is a negative angle, the tangent value of the pilot incident angle is taken as a slope, and the tangent value passes through the lower left point (x) of the rectanglemin,ymin) And the upper right point (x)max,ymax) Making two parallel lines which intersect (a) at the x-axismin0) and (a)max0) two points, wherein the part between the two points is the scanning range under the current lightning current incidence angle and the current lightning current amplitude;
step 2.2.5, the lightning current lead equation is expressed by a linear equation, and the expression is as follows:
x＝ky+a
where k is the slope of the linear equation, k is tan ψ, ψ is the pilot incident angle, and a is the intersection of the linear equation and the x-axis, and its value is at aminAnd amaxIn the meantime.
In the above-mentioned lightning trip-out rate scanning calculation method considering shielding of multiple lines in a power transmission channel, the implementation of step 3 includes:
step 3.1, scanning is carried out by changing the cross-sectional distance of the lightning current pilot equation on the x axis, intersection point coordinates of the current scanning equation and each stroke distance circle are obtained, wherein the intersection point with the highest vertical coordinate is determined as the current scanning point, and the stroke distance circle number of the current scanning point is judged;
step 3.2, according to the stroke distance circle number of the current scanning point, judging which base tower the point belongs to on the left side and the right side, and performing coordinate transformation according to intersection point coordinate information and the ground inclination angle of the corresponding tower, wherein the coordinate transformation formula is as follows:
Figure BDA0002808399390000041
wherein theta is the ground inclination angle of the tower to which the tower belongs, (x, y) is the coordinate of the current scanning point, and (x ', y') is the coordinate after the coordinate transformation of the scanning point;
step 3.3, the vertical distance between the current scanning point and the ground is taken as the transformed vertical coordinate y ', and the vertical coordinate y ' is compared with the ground hit distance, if y ' is more than or equal to rgIndicating that the current scanning point is not shielded by the ground, and the projection point of the point on the x axis belongs to the lightning triggering width of the current scanning circle; if when y' < rgAnd when the current scanning point is shielded by the ground, the projection point of the point on the x axis is not within the lightning guiding width of the current scanning circle, and the influence of the point is not considered when the lightning guiding width of the current scanning circle is calculated.
In the above-mentioned lightning trip-out rate scanning calculation method considering shielding of multiple lines in a power transmission channel, the implementation of step 5 includes:
step 5.1, the lightning strike width of the lead A is W under the lightning current amplitude value I and the lightning leading incident angle psi1(I, ψ); considering the pilot incidence angle distribution function p (ψ) and the lightning current amplitude distribution function p (i), the calculation formula of the wire a shielding failure rate W is:
Figure BDA0002808399390000042
in the formula IcFor lightning-resistance of wire around the wire, Wall(I, psi) is the sum of lightning-induced widths of all lightning conductors and conducting wires of the tower where the conducting wire is located, and the calculation expression is as follows:
Figure BDA0002808399390000043
calculating the tripping rate of the corresponding lead of different attack distance circles based on the lightning pilot incident angle distribution function;
step 5.2, according to the number N of lightning falling times of the lineLAnd obtaining the wire shielding failure trip rate SFFOR as follows:
SFFOR＝NLηW
wherein SFFOR is the shielding failure times, times/100 km/a; n is a radical ofgIn order of density of ground lightning, times/100 km.a, wherein the number of landings of the line is NLCan be calculated by a rule recommendation formula, and the expression is as follows
Figure BDA0002808399390000044
In the formula, hTIs the height of the tower, m; b is the distance between two ground wires, m; n is a radical ofgIs the density of the ground flash, times/km2·a。
Compared with the prior art, the invention has the beneficial effects that: (1) in the classic EGM algorithm, the position relation of the impact distance circle is easy to process because only the condition of a single-base-tower power transmission line is considered; under the condition that a plurality of base tower lines of a dense power transmission channel are erected in parallel, the number of wire and ground wire impact distance circles is rapidly increased, the position relation of the circles becomes very complex, repeated or omitted conditions are avoided when the position relation of the impact distance circles is discussed in a classification mode, the traditional shielding impact calculation method based on the geometric relation becomes extremely complex, and the calculation problem of the shielding condition of the plurality of lines can be well processed by the scanning calculation method adopted at present.
(2) By adopting a new computer scanning method, based on the translation of the lightning leading linear equation, the specific exposed arc section of the strike distance circle of each wire of the power transmission line does not need to be solved, the lightning triggering width of the projection of the corresponding exposed arc section on the ground is directly obtained, and the processing of the position relation of different strike distance circles is skipped, so that the method is suitable for the complicated multi-circuit tower type shielding strike calculation.
(3) In the invention, only the number of the hit distance circles and the positions of the hit distance circles are considered in the calculation, so that the tower types (such as a wine glass tower and a cat head tower) with the same total number of lines and different position relations can be realized by the same calculation program, and the calculation program has universal applicability.
The method is suitable for the shielding effect shielding around attack calculation considering that a plurality of lines are erected in the power transmission channel in parallel.
Drawings
FIG. 1 is a flow chart of a lightning trip-out rate scan calculation method according to an embodiment of the invention;
FIG. 2 is a schematic diagram of the number of pitch circle and the scanning range according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a scanning point selection according to an embodiment of the present invention;
fig. 4 is a schematic diagram of determining the width of the pitch circle lightning guide based on a scanning method according to an embodiment of the invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the following embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
The present invention is further illustrated by the following examples, which are not to be construed as limiting the invention.
In the aspect of research on shielding effects of multiple lines erected in parallel in a power transmission channel, the embodiment provides a shielding failure calculation method based on an EGM (earth generalized model) calculation principle and through computer scanning traversal solution, so that the shielding failure trip rate of each base tower line with the shielding effects of the multiple lines erected in parallel in the power transmission channel is obtained, and guidance is provided for lightning damage risk assessment and lightning protection measure transformation work of dense channel lines.
As shown in fig. 1, a lightning trip rate scanning calculation method considering shielding effects of a plurality of lines in a power transmission channel comprises the following steps:
and S1, acquiring tower parameters, terrain parameters and lightning activity parameters of the left and right lines of the power transmission channel.
And collecting power transmission channel line parameters, lightning activity and topographic data to prepare for subsequent modeling and calculation.
And inquiring a transmission line operation database to obtain left and right line parameters of the transmission line, wherein the left and right line parameters comprise the intermediate distance, the height and the voltage grade of each phase ground wire of the left and right towers, the ground inclination angles of the left and right towers in the channel, the altitude of the ground where the left and right towers are located and the distance between the left and right towers are obtained according to the GIS system. According to the thunder and lightning positioning system, the thunder and lightning activity condition of the line in the power transmission channel is obtained, if the actual thunder and lightning activity information cannot be obtained, the following recommended formula can be selected:
the lightning amplitude probability distribution, the lightning incident angle probability distribution and the striking rod rate are as follows: (actual can be selected according to the actual thunder and lightning data of the local area)
Figure BDA0002808399390000061
P(φ)＝0.75×cos3(φ)
g＝1:4
Wherein I is the lightning current amplitude; phi is the lightning incident angle (-90 degrees); g is the striking rod rate.
S2, determining relevant parameters of the electrical geometry models such as a strike distance and a strike distance coefficient, numbering all ground wire strike distance circles of the left and right side towers, and determining a lightning current pilot initiation equation and a current lightning current pilot amplitude and a scanning range under the current lightning current pilot amplitude and the current lightning current amplitude according to the current lightning current incident angle and the current lightning current amplitude.
And determining relevant parameters of the electric geometric model such as a range and a range coefficient. In the method of the embodiment, the impact distance and the impact distance coefficient adopt the following formulas:
rc＝10×I0.65，
rgw＝kgw×rc，
rg＝kg×rc，
wherein r isc,rgw,rgRespectively corresponding to the conducting wire, the lightning conductor and the earth strike distance; k is a radical ofgwAnd kgRespectively, the lightning conductor strike distance coefficient and the earth strike distance coefficient, wherein kgw＝1,kg0.9; i is a thunderThe magnitude of the current.
All ground wire strike distance circles of the left side tower and the right side tower are numbered, a lightning current pilot starting equation and a current lightning current pilot amplitude and a current scanning range under the incident angle are determined according to the current lightning current incident angle and the current lightning current amplitude, and the strike distance circle number and the scanning range are shown in figure 2.
S3, scanning is carried out by changing the cross-sectional distance of the lightning current pilot equation on the x axis, the intersection point coordinates of the current scanning equation and each stroke distance circle are obtained, wherein the intersection point with the highest vertical coordinate is determined as the current scanning point, and the number of the stroke distance circle to which the current scanning point belongs and whether the stroke distance circle is shielded by the ground are judged.
Scanning is carried out by changing the cross-sectional distance of the lightning current pilot equation on the x axis, the intersection point coordinates of the current scanning equation and each stroke distance circle are obtained, the intersection point with the highest vertical coordinate is determined as the current scanning point, the number of the stroke distance circle to which the current scanning point belongs is judged, and the scanning point is selected as shown in figure 3.
And judging which base tower the point belongs to on the left side and the right side according to the stroke distance circle number of the current scanning point, carrying out coordinate transformation according to intersection point coordinate information and the ground inclination angle of the corresponding tower, acquiring the vertical distance between the current scanning point and the ground, comparing the vertical distance with the ground stroke distance, and judging whether the current scanning point is shielded by the ground.
S4, when the number of the pitch circle to which the scanning point belongs changes, the scanning of the previous pitch circle under the leading equation is completed, the lightning strike width of the previous pitch circle under the current lightning leading equation is obtained, and the lightning strike widths of all pitch circles are sequentially determined in the scanning process of the lightning leading according to the method.
The schematic diagram of the scanning method for determining the width of the pitch circle lightning strike is shown in figure 4.
And S5, circularly considering lightning strike widths of all the strike distance circles under different lightning incident angles and lightning current amplitudes, and calculating the tripping rate of the corresponding lead of the different strike distance circles based on the lightning leading incident angle distribution function and the lightning current amplitude distribution function.
In specific implementation, as shown in fig. 1, a lightning trip-out rate scanning calculation method considering shielding effects of multiple lines in a power transmission channel includes the following steps:
the method comprises the steps of firstly, inquiring an electric transmission line operation database to obtain left and right line parameters of the electric transmission line, wherein the left and right line parameters comprise the intermediate distance, the height and the voltage grade of each phase ground wire of a left pole tower and a right pole tower, and the ground inclination angles of the left and right pole towers in a channel, the altitude of the ground where the left and right pole towers are located and the distance between the left and right pole towers are obtained according to a. According to the thunder and lightning positioning system, the thunder and lightning activity condition of the line in the power transmission channel is obtained, if the actual thunder and lightning activity information cannot be obtained, the actual thunder and lightning activity information can be selected according to the following recommendation formula
The lightning amplitude probability distribution, the lightning incident angle probability distribution and the striking rod rate are as follows: (actual can be selected according to the actual thunder and lightning data of the local area)
Figure BDA0002808399390000071
P(φ)＝0.75×cos3(φ)
g＝1:4
Wherein I is the lightning current amplitude; phi is the lightning incident angle (-90 degrees); g is the striking rod rate;
and step two, determining relevant parameters of the electric geometric model such as a range and a range coefficient. In the method, the impact distance and the impact distance coefficient adopt the following formulas:
rc＝10×I0.65，
rgw＝kgw×rc，
rg＝kg×rc，
wherein r isc,rgw,rgRespectively corresponding to the conducting wire, the lightning conductor and the earth strike distance; k is a radical ofgwAnd kgRespectively, the lightning conductor strike distance coefficient and the earth strike distance coefficient, wherein kgw＝1,kg0.9; i is the lightning current amplitude.
Numbering all ground wire strike distance circles of the left and right side towers, and determining a scanning range and a lightning current pilot equation according to the current lightning current incident angle and the current lightning current amplitude, wherein the scanning range is obtained by the following method:
A. and calculating the stroke distances of all stroke distance circles according to the amplitude of the lightning current, and drawing all stroke distance circles.
B. Acquiring the minimum value x of the abscissa of all the hit distance circlesminOrdinate minimum value yminMaximum value x of abscissamaxMaximum value y of ordinatemaxAnd with (x)min,ymin),(xmin,ymax),(xmax,ymin)，(xmax,ymax) Four points are made into a rectangle.
C. When the pilot incident angle is positive, the pilot incident angle tangent value is used as a slope, and the pilot incident angle tangent value passes through the upper left point (x) of the rectanglemin,ymax) And lower right point (x)max,ymin) Making two parallel lines which intersect (a) at the x-axismin0) and (a)max0) two points, and the portion between the two points is the scanning range under the current lightning current incidence angle and the current amplitude.
D. When the pilot incidence angle is negative, the tangent value of the pilot incidence angle is used as a slope, and the slope passes through the lower left point (x) of the rectanglemin,ymin) And the upper right point (x)max,ymax) Making two parallel lines which intersect (a) at the x-axismin0) and (a)max0) two points, and the portion between the two points is the scanning range under the current lightning current incidence angle and the current amplitude.
The lightning leading equation is expressed by a linear equation as follows
x＝ky+a
Where k is the slope of the linear equation, k is tan ψ, ψ is the pilot incident angle, and a is the intersection of the linear equation and the x-axis, and its value is at aminAnd amaxIn the meantime.
And step three, scanning by changing the cross-sectional distance of the lightning current pilot equation on the x axis, obtaining the coordinates of intersection points of the current scanning equation and each stroke distance circle, wherein the intersection point with the highest vertical coordinate is determined as the current scanning point, and judging the number of the stroke distance circle to which the current scanning point belongs, wherein the selection schematic diagram of the scanning point is shown in figure 3.
Judging which base tower the point belongs to on the left side and the right side according to the stroke distance circle number of the current scanning point, and carrying out coordinate transformation according to intersection point coordinate information and the ground inclination angle of the corresponding tower, wherein the coordinate transformation formula is as follows:
Figure BDA0002808399390000091
wherein theta is the ground inclination angle of the tower to which the tower belongs, (x, y) is the coordinate of the current scanning point, and (x ', y') is the coordinate after the coordinate transformation of the scanning point.
Taking the vertical distance between the current scanning point and the ground as the transformed vertical coordinate y ', comparing the vertical distance with the ground hit distance, and if y' is more than or equal to rgIndicating that the current scanning point is not shielded by the ground, and the projection point of the point on the x axis belongs to the lightning triggering width of the current scanning circle; if when y' < rgAnd when the current scanning point is shielded by the ground, the projection point of the point on the x axis is not within the lightning guiding width of the current scanning circle, and the influence of the point is not considered when the lightning guiding width of the current scanning circle is calculated.
Step four, when the number of the pitch circle to which the scanning point belongs changes, the scanning of the previous pitch circle under the pilot equation is completed, the lightning triggering width of the previous pitch circle under the current lightning pilot equation is obtained, the lightning triggering widths of all pitch circles are sequentially determined in the scanning process of the lightning pilot according to the method, and a schematic diagram of the lightning triggering width of the pitch circle determined by the scanning method is shown in fig. 4.
And step five, circularly considering the lightning strike widths of all the stroke distance circles under different lightning incident angles and lightning current amplitudes, taking a lead A as an example, wherein the lightning strike width is W under the lightning current amplitude I and the lightning leading incident angle psi1(I, ψ). Considering the pilot incidence angle distribution function p (ψ) and the lightning current amplitude distribution function p (i), the calculation formula of the wire a shielding failure rate W is:
Figure BDA0002808399390000092
in the formula IcFor lightning-resistance of wire around the wire, Wall(I, ψ) ofThe sum of lightning-induced widths of all lightning conductors and the conducting wire of the tower where the conducting wire is located is calculated as follows:
Figure BDA0002808399390000093
and calculating the tripping rate of the corresponding wire of different attack distance circles based on the lightning pilot incident angle distribution function.
Finally, according to the number N of lightning strike of the lineLAnd obtaining the wire shielding failure trip rate SFFOR as follows:
SFFOR＝NLηW
wherein SFFOR is the shielding failure times, times/100 km/a; n is a radical ofgThe number of lightning strike on the line is N/(100 km.a)LCan be calculated by a rule recommendation formula, and the expression is as follows
Figure BDA0002808399390000101
In the formula, hTIs the height of the tower, m; b is the distance between two ground wires, m; n is a radical ofgIs the density of the lightning, times/(km)2·a)。
While the invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention.

Claims (5)

1. The lightning trip-out rate scanning and calculating method considering shielding of a plurality of lines in a power transmission channel is characterized by comprising the following steps of:
step 1, obtaining tower parameters, terrain parameters and lightning activity parameters of left and right lines of a power transmission channel;
step 2, determining a attack distance and attack distance coefficients, numbering all ground wire attack distance circles of the left side tower and the right side tower, and determining a lightning current pilot starting equation and a current lightning current pilot amplitude and a scanning range under the current lightning current pilot amplitude and the current lightning current amplitude according to the current lightning current incident angle and the current lightning current amplitude;
step 3, scanning is carried out by changing the cross-sectional distance of the lightning current pilot equation on the x axis, the intersection point coordinates of the current scanning equation and each stroke distance circle are obtained, wherein the intersection point with the highest vertical coordinate is determined as the current scanning point, and the number of the stroke distance circle to which the current scanning point belongs and whether the stroke distance circle is shielded by the ground are judged;
step 4, when the number of the pitch circle to which the scanning point belongs changes, the scanning of the previous pitch circle under the leading equation is represented to be completed, the lightning triggering width of the previous pitch circle under the current lightning leading equation is obtained, and the lightning triggering widths of all pitch circles are sequentially determined in the scanning process of the lightning leading according to the method;
and 5, circularly considering the lightning triggering widths of all the stroke distance circles under different lightning incident angles and lightning current amplitudes, and calculating the tripping rate of the corresponding lead of the different stroke distance circles based on the lightning leading incident angle distribution function and the lightning current amplitude distribution function.
2. The method of claim 1 for calculating a lightning trip rate scan taking into account multiple line shielding in a power transmission channel, wherein the step 1 is implemented by:
inquiring a transmission line operation database to obtain left and right line parameters of the transmission line, wherein the left and right line parameters comprise the intermediate distance, the height and the voltage grade of each phase conducting ground wire of a left tower and a right tower, and the ground inclination angles of the left tower and the right tower in a channel, the altitude of the ground where the left tower and the right tower are located and the distance between the left tower and the right tower are obtained according to a GIS system; acquiring the lightning activity condition of a line in a power transmission channel according to the recorded data of the lightning positioning system; or selecting according to the following formula;
lightning amplitude probability distribution, lightning incidence angle probability distribution and striking rod rate:
Figure FDA0002808399380000011
P(φ)＝0.75×cos3(φ)
g＝1:4
wherein I is the lightning current amplitude; phi is the lightning incident angle (-90 degrees); g is the striking rod rate.
3. The method of claim 1 for calculating a lightning trip rate scan taking into account multiple line shielding in a power transmission channel, wherein the step 2 is implemented by:
step 2.1, determining the hit distance and the hit distance coefficient by adopting the following formula:
rc＝10×I0.65，
rgw＝kgw×rc，
rg＝kg×rc，
wherein r isc,rgw,rgRespectively corresponding to the conducting wire, the lightning conductor and the earth strike distance; k is a radical ofgwAnd kgRespectively a lightning conductor strike distance coefficient and an earth strike distance coefficient, wherein kgw＝1,kg0.9; i is the lightning current amplitude;
step 2.2, numbering all ground wire strike distance circles of the left and right side towers, and determining a lightning current pilot initial equation, a current lightning current pilot amplitude and a scanning range under the incident angle, the strike distance circle number and the scanning range according to the current lightning current incident angle and the current lightning current amplitude;
step 2.2.1, calculating the hitting distances of all hitting distance circles according to the lightning current amplitude, and drawing all hitting distance circles;
step 2.2.2, obtaining the minimum value x of the abscissa of all the stroke distance circlesminOrdinate minimum value yminMaximum value x of abscissamaxMaximum value y of ordinatemaxAnd with (x)min,ymin),(xmin,ymax),(xmax,ymin)，(xmax,ymax) Making four points into a rectangle;
step 2.2.3, when the pilot incident angle is a positive angle, the tangent value of the pilot incident angle is taken as a slope, and the positive angle and the slope respectively pass through the upper left point (x) of the rectanglemin,ymax) And lower right point (x)max,ymin) Making two parallel lines which intersect (a) at the x-axismin0) and (a)max0) two points, wherein the part between the two points is the scanning range under the current lightning current incidence angle and the current lightning current amplitude;
step 2.2.4, when the pilot incident angle is a negative angle, the tangent value of the pilot incident angle is taken as a slope, and the tangent value passes through the lower left point (x) of the rectanglemin,ymin) And the upper right point (x)max,ymax) Making two parallel lines which intersect (a) at the x-axismin0) and (a)max0) two points, wherein the part between the two points is the scanning range under the current lightning current incidence angle and the current lightning current amplitude;
step 2.2.5, the lightning current lead equation is expressed by a linear equation, and the expression is as follows:
x＝ky+a
where k is the slope of the linear equation, k is tan ψ, ψ is the pilot incident angle, and a is the intersection of the linear equation and the x-axis, and its value is at aminAnd amaxIn the meantime.
4. The method of claim 1 for calculating a lightning trip rate scan taking into account multiple line shielding in a power transmission channel, wherein the step 3 is implemented by:
step 3.1, scanning is carried out by changing the cross-sectional distance of the lightning current pilot equation on the x axis, intersection point coordinates of the current scanning equation and each stroke distance circle are obtained, wherein the intersection point with the highest vertical coordinate is determined as the current scanning point, and the stroke distance circle number of the current scanning point is judged;
step 3.2, according to the stroke distance circle number of the current scanning point, judging which base tower the point belongs to on the left side and the right side, and performing coordinate transformation according to intersection point coordinate information and the ground inclination angle of the corresponding tower, wherein the coordinate transformation formula is as follows:
Figure FDA0002808399380000031
wherein theta is the ground inclination angle of the tower to which the tower belongs, (x, y) is the coordinate of the current scanning point, and (x ', y') is the coordinate after the coordinate transformation of the scanning point;
step 3.3, the vertical distance between the current scanning point and the ground is taken as the transformed vertical coordinate y ', and the vertical coordinate y ' is compared with the ground hit distance, if y ' is more than or equal to rgIndicating that the current scanning point is not shielded by the ground, and the projection point of the point on the x axis belongs to the lightning triggering width of the current scanning circle; if when y' < rgAnd when the current scanning point is shielded by the ground, the projection point of the point on the x axis is not within the lightning guiding width of the current scanning circle, and the influence of the point is not considered when the lightning guiding width of the current scanning circle is calculated.
5. The method of claim 1 for calculating a lightning trip rate scan taking into account multiple line shielding in a power transmission channel, wherein the step 5 is implemented by:
step 5.1, the lightning strike width of the lead A is W under the lightning current amplitude value I and the lightning leading incident angle psi1(I, ψ); considering the pilot incidence angle distribution function p (ψ) and the lightning current amplitude distribution function p (i), the calculation formula of the wire a shielding failure rate W is:
Figure FDA0002808399380000032
in the formula IcFor lightning-resistance of wire around the wire, Wall(I, psi) is the sum of lightning-induced widths of all lightning conductors and conducting wires of the tower where the conducting wire is located, and the calculation expression is as follows:
Figure FDA0002808399380000033
calculating the tripping rate of the corresponding lead of different attack distance circles based on the lightning pilot incident angle distribution function;
step 5.2, according to the number N of lightning falling times of the lineLAnd obtaining the wire shielding failure trip rate SFFOR as follows:
SFFOR＝NLηW
wherein SFFOR is the shielding failure times, times/100 km/a; n is a radical ofgTo density of the lightning, sub/100 kma, wherein
Number of landmine times N of lineLCan be calculated by a rule recommendation formula, and the expression is as follows
Figure FDA0002808399380000041
In the formula, hTIs the height of the tower, m; b is the distance between two ground wires, m; n is a radical ofgIs the density of the ground flash, times/km2·a。
CN202011380660.9A 2020-11-30 2020-11-30 Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel Active CN112541211B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011380660.9A CN112541211B (en) 2020-11-30 2020-11-30 Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011380660.9A CN112541211B (en) 2020-11-30 2020-11-30 Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel

Publications (2)

Publication Number Publication Date
CN112541211A true CN112541211A (en) 2021-03-23
CN112541211B CN112541211B (en) 2022-06-07

Family

ID=75016754

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011380660.9A Active CN112541211B (en) 2020-11-30 2020-11-30 Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel

Country Status (1)

Country Link
CN (1) CN112541211B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072992A (en) * 2010-10-26 2011-05-25 国网电力科学研究院武汉南瑞有限责任公司 Detailed terrain data-based transmission line failure-shielding and lightning-protection performance evaluation method
JP2013164651A (en) * 2012-02-09 2013-08-22 Chugoku Electric Power Co Inc:The Failure prediction device for lightning arrester for power transmission, failure prediction program and failure prediction method
CN103837769A (en) * 2014-02-27 2014-06-04 广西电网公司电力科学研究院 Lightening damage early-warning method and system for electric transmission line
CN105242133A (en) * 2015-09-18 2016-01-13 南京信息工程大学 Improved method for calculating lightning trip-out rate of distribution line
CN106919799A (en) * 2017-03-07 2017-07-04 四川大学 Phase conductor back flash-over rate computational methods on ultra high-tension transmission line staggered tower upward slope position
CN107703381A (en) * 2017-09-19 2018-02-16 国网浙江省电力公司 A kind of computational methods for screen effect of being struck by lightning between neighbouring transmission line of electricity
JP2019009951A (en) * 2017-06-28 2019-01-17 東京電力ホールディングス株式会社 Thunder accident risk determination system, thunder accident risk determination method, and program
CN111814355A (en) * 2020-07-24 2020-10-23 南方电网科学研究院有限责任公司 Method and device for protecting high-voltage transmission line from lightning shielding failure and terminal equipment

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072992A (en) * 2010-10-26 2011-05-25 国网电力科学研究院武汉南瑞有限责任公司 Detailed terrain data-based transmission line failure-shielding and lightning-protection performance evaluation method
JP2013164651A (en) * 2012-02-09 2013-08-22 Chugoku Electric Power Co Inc:The Failure prediction device for lightning arrester for power transmission, failure prediction program and failure prediction method
CN103837769A (en) * 2014-02-27 2014-06-04 广西电网公司电力科学研究院 Lightening damage early-warning method and system for electric transmission line
CN105242133A (en) * 2015-09-18 2016-01-13 南京信息工程大学 Improved method for calculating lightning trip-out rate of distribution line
CN106919799A (en) * 2017-03-07 2017-07-04 四川大学 Phase conductor back flash-over rate computational methods on ultra high-tension transmission line staggered tower upward slope position
JP2019009951A (en) * 2017-06-28 2019-01-17 東京電力ホールディングス株式会社 Thunder accident risk determination system, thunder accident risk determination method, and program
CN107703381A (en) * 2017-09-19 2018-02-16 国网浙江省电力公司 A kind of computational methods for screen effect of being struck by lightning between neighbouring transmission line of electricity
CN111814355A (en) * 2020-07-24 2020-10-23 南方电网科学研究院有限责任公司 Method and device for protecting high-voltage transmission line from lightning shielding failure and terminal equipment

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SAKAE TANIGUCHI等: "Improved Method of Calculating Lightning Stroke Rate to Large-sized Transmission Lines Based on Electric Geometry Model", 《IEEE TRANSACTIONS ON DIELECTRICS AND ELECTRICAL INSULATION》 *
樊春雷等: "高压输电线路绕击跳闸率的研究", 《电瓷避雷器》 *
耿厚花: "特高压双回输电线路耐雷计算方法", 《中国电力》 *
苏杰等: "高速铁路牵引网暴露弧面法直击雷跳闸率计算", 《中国电机工程学报》 *

Also Published As

Publication number Publication date
CN112541211B (en) 2022-06-07

Similar Documents

Publication Publication Date Title
CN110309527B (en) Electrical geometric model-based lightning damage risk assessment method for overhead distribution line
CN112883878B (en) Automatic point cloud classification method under transformer substation scene based on three-dimensional grid
CN110908014B (en) Galloping refined correction forecasting method and system
CN107958486A (en) A kind of generation method and device of conducting wire vector model
CN105242133B (en) A kind of improvement distribution line lightning trip-out rate calculation method
CN109064056B (en) Power transmission line lightning protection measure selection method based on grey correlation analysis method
CN111751638A (en) Analysis method, device and equipment for lightning trip-out rate of power transmission line and storage medium
CN112785117A (en) Multidimensional analysis method for risk of full-line lightning shielding failure tripping of power transmission line
CN107703381B (en) Method for calculating lightning stroke shielding effect between adjacent power transmission lines
CN113011099A (en) Method for calculating and correcting lightning trip-out rate of power transmission line
CN116109133A (en) Risk early warning method for thunderstorm weather power transmission line
CN115825984A (en) Power transmission line inspection system and inspection method based on laser point cloud
CN112541211B (en) Lightning trip-out rate scanning calculation method considering shielding of multiple lines in power transmission channel
CN106919799B (en) Phase conductor back flash-over rate calculation method on ultra high-tension transmission line staggered tower upward slope position
CN110687371B (en) Method and system for determining lightning shielding failure performance of same-tower multi-circuit line
CN104076250A (en) Method for analyzing shielding failure and shielding failure trip-out rate of upper-phase conductors of double-circuit lines on same tower
CN113267667B (en) Three-dimensional lightning trip-out rate calculation method considering cross-crossing power transmission line shielding
CN110543610A (en) Real-time lightning transmission line risk assessment method
Zhuang et al. Improve the electrogeometric model by the analysis results of leader propagation model for transmission lines
CN109521333A (en) The multiple-circuit on same tower Characteristic of Lightning Shielding Failure analysis method of shielding action between meter and conducting wire
CN113487601A (en) Method for calculating lightning protection angle based on high-precision three-dimensional point cloud data
Wen et al. Evaluation of whole span lightning shielding flashover risk based on 3-D laser scanning technology
CN112199841A (en) Simulation calculation method for differential evaluation of lightning damage
CN112053365A (en) Method for extracting laser point cloud data of power transmission line conductor based on linear characteristics
CN117078014B (en) Power transmission line lightning disaster risk level assessment method and device based on laser point cloud data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant