CN107436162A - A kind of power network line waves Occurrence forecast method and system - Google Patents

A kind of power network line waves Occurrence forecast method and system Download PDF

Info

Publication number
CN107436162A
CN107436162A CN201710640834.2A CN201710640834A CN107436162A CN 107436162 A CN107436162 A CN 107436162A CN 201710640834 A CN201710640834 A CN 201710640834A CN 107436162 A CN107436162 A CN 107436162A
Authority
CN
China
Prior art keywords
msub
mrow
temperature
analyzed
waving
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
CN201710640834.2A
Other languages
Chinese (zh)
Other versions
CN107436162B (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.)
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd
Hunan Xiangdian Test Research Institute Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd
Hunan Xiangdian Test Research Institute Co Ltd
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 State Grid Corp of China SGCC, State Grid Hunan Electric Power Co Ltd, Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd, Hunan Xiangdian Test Research Institute Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201710640834.2A priority Critical patent/CN107436162B/en
Publication of CN107436162A publication Critical patent/CN107436162A/en
Application granted granted Critical
Publication of CN107436162B publication Critical patent/CN107436162B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass

Abstract

The present invention relates to power transmission line plant maintenance field, disclose a kind of power network line and wave Occurrence forecast method and system, think power network icing fast prediction, carry out Anti-galloping emergency measure in advance to prepare, strong technical support is provided, the present invention collects the history to be analyzed waved a little and waves situation first, to the history situation of waving analyze and screen drawing and to be analyzed wave a temperature and relative humidity for waving initial stage;Warm and humid factor set is obtained according to temperature and relative humidity calculation, then using temperature as abscissa, relative humidity is that warm and humid factor set is plotted in rectangular coordinate system by ordinate, and the angle for calculating each warm and humid factor pair and abscissa obtains angle threshold to carry out sentencing knowledge to waving situation about waving a little to be analyzed;The temperature and relative humidity to be analyzed waved a little is monitored in real time, builds new warm and humid factor set, and waving situation about waving a little to be analyzed is carried out sentencing knowledge with reference to angle threshold.

Description

A kind of power network line waves Occurrence forecast method and system
Technical field
The present invention relates to power transmission line plant maintenance field, more particularly to a kind of power network line to wave Occurrence forecast method and be System.
Background technology
In recent years, as social and economic stability is fast-developing, power load increases, power network coverage rate augmentation, by easily waving The transmission line of electricity in region also increases therewith, simultaneously as Extreme Weather-climate Events occur again and again so that easily waves regional extent not Disconnected increase, in the winter time, powerline ice-covering is waved problem and highlighted increasingly, and serious threat the stability and safety of operation of power networks Property.
At present, for the forecast for transmission line galloping mainly by establishing model, calculating future waves index to carry out, should Method has the preferable value of forecasting, but waving for zonule shaft tower for large-scale transmission line galloping, due to existing The local meteorological change in field comparatively fast, is easily reported by mistake, it is necessary to consider the real time humiture change at scene, progress fast scroll forecast.So And using traditional forecasting procedure, computationally intensive due to needing the numerical computations to become more meticulous on a large scale, calculating time length, it is difficult to Meet to the quick determination requirement of line oscillation.
Therefore, now need to provide a kind of temperature and relative humidity progress real-time monitored that can be waved based on scene a little, so as to real The power network line of existing quick recognition waves Occurrence forecast method and system.
The content of the invention
Present invention aims at provide a kind of power network line to wave Occurrence forecast method and system, this method and system by temperature Degree and relative humidity synthesis carry out threshold calculations into a variable, and introduce coordinate space and enter line range judgement, make result more Intuitively, and using relative quantity as reference quantity, the scope of application is broad, introduces standard deviation as deviation range, makes result of calculation more Accurately.
To achieve the above object, the invention provides a kind of power network line to wave Occurrence forecast method, including:
Collect the history to be analyzed waved a little and wave situation, analysis is carried out to the history situation of waving and draws institute with screening State and to be analyzed wave a temperature and relative humidity for waving initial stage;
The temperature of initial stage is waved according to history and relative humidity calculation obtains temperature-wet factor set, then using temperature as Abscissa, relative humidity is that the temperature-wet factor set is plotted in rectangular coordinate system by ordinate, and it is warm-wet to calculate each The angle of factor pair and abscissa obtains angle threshold to carry out sentencing knowledge to the waving situation about waving a little to be analyzed;
Monitor in real time it is described it is to be analyzed wave a little current temperature and relative humidity, build new temperature-wet factor set, and tie Close the angle threshold the waving situation about waving a little to be analyzed is carried out sentencing knowledge.
To achieve the above object, the present invention provides a kind of power network line and waves Occurrence forecast system, including:
First module:Situation is waved for collecting the history to be analyzed waved a little, the history situation of waving is divided Analysis and screening, which are drawn, described to be analyzed waves a temperature and relative humidity for waving initial stage;
Second unit:Temperature-wet factor set is obtained for waving the temperature of initial stage and relative humidity calculation according to history, Then using temperature as abscissa, relative humidity is that the temperature-wet factor set is plotted in rectangular coordinate system, and calculated by ordinate The angle of each temperature-wet factor pair and abscissa obtains angle threshold to be carried out to the situation of waving to be analyzed waved a little Sentence knowledge;
Third unit:For monitor in real time it is described it is to be analyzed wave a little current temperature and relative humidity, build new temperature- Wet factor set, and the waving situation about waving a little to be analyzed is carried out sentencing knowledge with reference to the angle threshold.
The invention has the advantages that:
The present invention provides a kind of power network line and waves Occurrence forecast method and system, first to the history to be analyzed waved a little Wave situation and screened according to the duration is waved, the temperature and relative humidity waving event and early stage occurs after Analysis and Screening Relation and temperature-wet factor set is calculated accordingly, further obtain angle threshold, and monitor described to be analyzed wave a little in real time Temperature and relative humidity, then sentence to know by angle threshold and described to be analyzed wave whether a future waves;This method Temperature is subjected to threshold calculations into a variable with relative humidity synthesis with system, and introduces coordinate space and enters line range judgement, Make result more directly perceived, and using relative quantity as reference quantity, the scope of application is broad, introduces standard deviation as deviation range, makes meter It is more accurate to calculate result, in addition, during analysis, can also be in time according to temperature and relative humidity the modification angle monitored Threshold value, accuracy and the flexibility of prediction result are improved, be power network icing fast prediction, carry out Anti-galloping emergency measure in advance Prepare, there is provided strong technical support.
Below with reference to accompanying drawings, the present invention is further detailed explanation.
Brief description of the drawings
The accompanying drawing for forming the part of the application is used for providing a further understanding of the present invention, schematic reality of the invention Apply example and its illustrate to be used to explain the present invention, do not form inappropriate limitation of the present invention.In the accompanying drawings:
Fig. 1 is the process schematic of the preferred embodiment of the present invention;
Fig. 2 is that to be analyzed wave of the preferred embodiment of the present invention will not a little wave schematic diagram;
Fig. 3 is that to be analyzed wave of the preferred embodiment of the present invention can wave schematic diagram.
Embodiment
Embodiments of the invention are described in detail below in conjunction with accompanying drawing, but the present invention can be defined by the claims Implement with the multitude of different ways of covering.
Embodiment 1
The present embodiment provides a kind of power network line and waves Occurrence forecast method, reference picture 1, including:
Collect the history to be analyzed waved a little and wave situation, analysis is carried out to the history situation of waving and draws institute with screening State and to be analyzed wave a temperature and relative humidity for waving initial stage;
The temperature of initial stage is waved according to history and relative humidity calculation obtains temperature-wet factor set, then using temperature as Abscissa, relative humidity is that the temperature-wet factor set is plotted in rectangular coordinate system by ordinate, and it is warm-wet to calculate each The angle of factor pair and abscissa obtains angle threshold to carry out sentencing knowledge to the waving situation about waving a little to be analyzed;
Monitor in real time it is described it is to be analyzed wave a little current temperature and relative humidity, build new temperature-wet factor set, and tie Close the angle threshold the waving situation about waving a little to be analyzed is carried out sentencing knowledge.
Specifically, by taking certain 500kV transmission line gallopings point of certain province's power network as an example, the transmission line galloping point is designated as Q, collect Q points and wave event 10 times in history, represented to wave event each time with i, wherein, 0<I≤10, collection are waved each time (it is determined that to hour) S between at the beginning of eventi, end time (it is determined that to hour) EiAnd early stage takes place in the event of waving Temperature TiWith relative humidity Hi, what deserves to be explained is, the event of waving take place refer to early stage wave event occur before it is one small When, and waving event each time is occurred in the case of less than 0 °.
Record the temperature for waving event collected and relative humidity data tentatively establishes temperature-wet factor set (Ti, Hi), Q points Temperature and relative humidity value it is as shown in table 1 below:
Table 1Q points 10 temperature and rh value for waving event in history
Sequence number Temperature Ti(℃) Relative humidity Hi(%)
1 -1.99 95
2 -1.84 98
3 -1.74 93
4 -2.06 95
5 -1.87 100
6 -1.68 98
7 -1.9 93
8 -1.52 96
9 -1.68 95
10 -1.46 98
According to S between the waving each time at the beginning of event of collectioniWith end time EiHolding for event is waved in calculating each time Continuous duration:
Li=Ei-Si
In formula, LiRepresent that Q points wave the duration of event each time, history is waved into duration and is more than a certain setting Value wave event reservation, and by history wave duration less than or equal to a certain setting value wave event reject, it is necessary to Illustrate, because the line oscillation duration is more long, the circuit within usual 1 hour, which is swung, not to be regarded as waving, Therefore, the setting value is set as 2 hours, screened according to above-mentioned screening conditions 10 events of waving historical to Q points, The lasting event for waving event is all higher than 2 hours due to above-mentioned 10 times, therefore according to waving the new temperature of event establishment-wet for above-mentioned 10 times Factor set (Tj, Hj), wherein 0<j≤10.
Choose above-mentioned 10 times and wave the temperature maximum T that early stage occurs for eventmax, temperature minimum value TminAnd relative humidity Maximum Hmax, relative humidity minimum value HminObtain Tmax=-1.46, Tmin=-2.06, Hmax=100, Hmin=93, and will sieve The data of choosing substitute into formula:
With by temperature-wet factor set (Tj, Hj) enter line translation, obtain the temperature being easy to be drawn in two-dimensional coordinate of specification-wet because Subset (tj, hj), wherein 0<J≤10, its result of calculation are as shown in table 2 below:
Temperature and rh value after the conversion of table 2
Sequence number Temperature ti(℃) Relative humidity hi(%)
1 0.12 0.29
2 0.37 0.71
3 0.53 0.00
4 0.00 0.29
5 0.32 1.00
6 0.63 0.71
7 0.27 0.00
8 0.90 0.43
9 0.63 0.29
10 1.00 0.71
Further, using temperature as abscissa, relative humidity is ordinate by the temperature obtained after conversion-wet factor set (tj, hj) be plotted in two-dimensional coordinate space, and above-mentioned related data is substituted into formulaCalculate temperature-wet factor set (tj, hj) in the angle theta of each temperature-wet factor pair and abscissaj, and according to formulaWithθ is calculatedjAverage value be θmean=0.74, standard deviation θstd=0.54, and according to θj's Average value and standard deviation pass through formula θmaxmeanstdAnd θminmeanstdThe maximum θ of angle threshold is calculatedmax =1.28, the minimum value θ of angle thresholdmin=0.20, then the scope of angle threshold is [0.20,1.28].
Then, the current temperature T of Q points is monitored in real time by measurement apparatusobsWith relative humidity Hobs, built based on above-mentioned thought Current temperature-wet the factor set of vertical Q points, therefrom chooses a data Tobs=-1.5 DEG C and Hobs=94%, and by the data group Cheng Wen-wet factor pair (tobs, hobs), and temperature-wet factor pair (t is obtained after calculating and conversionobs, hobs) value for (0.93, 0.14) its angle value θ, is then calculatedobs=0.15, as shown in Fig. 2 due toIt can then obtain Go out conclusion, the point Q to be analyzed that waves will not be waved for 1 hour in future.It should be noted that in the current number to monitoring in real time According to, it is necessary to judge data, that is, the temperature data and relative humidity data in the data being monitored to need before being analyzed To meet the maximum less than or equal to historical data simultaneously, more than or equal to the minimum value of historical data, to ensure prison in real time The data of survey can meet temperature that the environmental condition of icing and ensureing obtains during the data to monitoring in real time are analyzed- Wet factor pair (tobs, hobs) in tobsAnd hobsValue be between 0 and 1.In addition, what deserves to be explained is, if the number of monitoring Temperature data and relative humidity data in are more than or equal to without the maximum met less than or equal to historical data simultaneously Any of temperature data and relative humidity data in the minimum value of historical data, or the data of monitoring data do not meet Less than or equal to the maximum of historical data, more than or equal to the minimum value of historical data, still, the data of the monitoring are in reality Waved in situation, then it is right in historical data to be added to the temperature data in the data of the monitoring and relative humidity data Historical data is adjusted, be then based on above-mentioned thought recalculate angle threshold with improve the accuracy of prediction result with flexibly Property.
Embodiment 2
From the temperature T of the Q points real-time monitoredobsWith relative humidity HobsIt is middle to choose other a data Tobs=-1.6 DEG C And Hobs=98%, and the data is formed into temperature-wet factor pair, and obtain temperature-wet factor pair (t by calculating and conversionobs, hobs) it is (0.76,0.71), its angle value θ is then calculatedobs=0.75, as shown in figure 3, due to 0.75 ∈ [0.20, 1.28], then it may be concluded that the point Q to be analyzed that waves can be waved for 1 hour in future.
Embodiment 3
Corresponding with above method embodiment, the present embodiment discloses a kind of power network line and waves Occurrence forecast system, bag Include:
First module:Situation is waved for collecting the history to be analyzed waved a little, the history situation of waving is divided Analysis and screening, which are drawn, described to be analyzed waves a temperature and relative humidity for waving initial stage;
Second unit:Temperature-wet factor set is obtained for waving the temperature of initial stage and relative humidity calculation according to history, Then using temperature as abscissa, relative humidity is that the temperature-wet factor set is plotted in rectangular coordinate system, and calculated by ordinate The angle of each temperature-wet factor pair and abscissa obtains angle threshold to be carried out to the situation of waving to be analyzed waved a little Sentence knowledge;
Third unit:For monitor in real time it is described it is to be analyzed wave a little current temperature and relative humidity, build new temperature- Wet factor set, and the waving situation about waving a little to be analyzed is carried out sentencing knowledge with reference to the angle threshold.
The concrete processing procedure of above-mentioned each unit can refer to above method embodiment, repeat no more.
As described above, the present invention, which provides a kind of power network line, waves Occurrence forecast method and system, first to dance to be analyzed The history of dynamic point is waved situation and screened according to the duration is waved, the temperature waved event and early stage occurs after Analysis and Screening With the relation of relative humidity and temperature-wet factor set is calculated accordingly, further obtain angle threshold, and treated in real time described in monitoring Temperature and relative humidity a little are waved in analysis, then sentence that knowledge is described to be analyzed to wave whether a future waves by angle threshold It is dynamic;Temperature is carried out threshold calculations by this method and system with relative humidity synthesis into a variable, and introduces coordinate space progress Scope judges, makes result more directly perceived, and using relative quantity as reference quantity, the scope of application is broad, introduces standard deviation as deviation Scope, make result of calculation more accurate, be power network icing fast prediction, carry out Anti-galloping emergency measure in advance and prepare, there is provided Strong technical support.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, for the skill of this area For art personnel, the present invention can have various modifications and variations.Within the spirit and principles of the invention, that is made any repaiies Change, equivalent substitution, improvement etc., should be included in the scope of the protection.

Claims (4)

1. a kind of power network line waves Occurrence forecast method, it is characterised in that including:
Collect the history to be analyzed waved a little and wave situation, analysis is carried out to the history situation of waving and draws described treat with screening A temperature and relative humidity for waving initial stage are waved in analysis;
The temperature of initial stage is waved according to history and relative humidity calculation obtains temperature-wet factor set, then using temperature as horizontal seat Mark, relative humidity are that the temperature-wet factor set is plotted in rectangular coordinate system by ordinate, and are calculated every in temperature-wet factor set The angle of one temperature-wet factor pair and abscissa obtains angle threshold to sentence to the waving situation about waving a little to be analyzed Know;
Monitor in real time it is described it is to be analyzed wave a little current temperature and relative humidity, build new temperature-wet factor set, and combine institute State angle threshold the waving situation about waving a little to be analyzed is carried out sentencing knowledge.
2. power network line according to claim 1 waves Occurrence forecast method, it is characterised in that
Calculate temperature-wet factor set and specifically include following steps:
(1) select it is to be analyzed wave point Q, collect Q points and wave in history in event for n times, between waving each time at the beginning of event i Si, end time Ei, and wave the temperature T of initial stageiWith relative humidity Hi, it is preliminary to establish temperature-wet factor set (Ti, Hi), wherein 0<i≤n;
(2) duration to be analyzed waved point Q and wave event each time is calculated:
Li=Ei-Si
In formula, LiRepresent the duration to be analyzed waved point Q and wave event each time;
By history wave duration more than a certain setting value wave event retain, and by history wave duration be less than or Event of waving equal to a certain setting value is rejected, and is obtained m times by screening and is waved event, m≤n, and establishes new temperature-wet accordingly Factor set (Tj, Hj), wherein 0<j≤m;
(3) the temperature maximum T obtained by screening is chosenmax, temperature minimum value TminAnd relative humidity maximum Hmax, phase To humidity minimum value Hmin,
Pass through below equation:
<mrow> <msub> <mi>t</mi> <mi>j</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>min</mi> </msub> </mrow> <mrow> <msub> <mi>T</mi> <mi>max</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>min</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>
<mrow> <msub> <mi>h</mi> <mi>j</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>H</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>H</mi> <mi>max</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>
Simultaneously temperature described in shift step (2)-wet factor set is calculated, obtains temperature-wet factor set (t of specificationj, hj), wherein 0<j≤ m;
Calculate angle threshold and specifically include following steps:
(1) using temperature as abscissa, relative humidity is ordinate, by temperature-wet factor set (tj, hj) it is plotted in two-dimensional coordinate space In, calculate temperature-wet factor set (tj, hj) in the angle theta of each temperature-wet factor pair and abscissaj
<mrow> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>h</mi> <mi>j</mi> </msub> <msub> <mi>t</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
Calculate θjAverage value θmeanWith standard deviation θstd
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>m</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>;</mo> </mrow> 1
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>s</mi> <mi>t</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mi>m</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>;</mo> </mrow>
According to θjAverage value θmeanWith standard deviation θstdCalculate the maximum of angle threshold and the minimum value of angle threshold:
θmaxmeanstd
θminmeanstd
It is [θ to obtain angle threshold scopemin, θmax];
Monitor in real time and to be analyzed wave the current temperature T of point QobsWith relative humidity Hobs, calculated by above-mentioned formula and change got in return To it is to be analyzed wave the current temperature of point Q-wet factor set and therefrom choose to be analyzed wave the current temperature of point Q-wet factor pair (tobs, hobs), angle, θ is calculatedobs, by θobsIt is compared with above-mentioned angle threshold scope, works as θobs∈[θmin, θmax] when, judge The point Q to be analyzed that waves can be waved in following certain time;WhenWhen, judge following one The interior point Q to be analyzed that waves that fixes time will not be waved.
3. a kind of power network line waves Occurrence forecast system, it is characterised in that including:
First module:Wave situation for collecting the history to be analyzed waved a little, to the history situation of waving carry out analysis with Screening, which is drawn, described to be analyzed waves a temperature and relative humidity for waving initial stage;
Second unit:Temperature-wet factor set is obtained for waving the temperature of initial stage and relative humidity calculation according to history, then Using temperature as abscissa, relative humidity is that the temperature-wet factor set is plotted in rectangular coordinate system by ordinate, and is calculated each The angle of individual temperature-wet factor pair and abscissa obtains angle threshold to sentence to the waving situation about waving a little to be analyzed Know;
Third unit:For monitor in real time it is described it is to be analyzed wave a little current temperature and relative humidity, build new temperature-wet because Subset, and the waving situation about waving a little to be analyzed is carried out sentencing knowledge with reference to the angle threshold.
4. a kind of power network line according to claim 3 waves Occurrence forecast system, it is characterised in that
Temperature-wet factor set is calculated in the second unit and specifically includes following steps:
(1) select it is to be analyzed wave point Q, collect Q points and wave in history in event for n times, between waving each time at the beginning of event i Si, end time Ei, and wave the temperature T of initial stageiWith relative humidity Hi, it is preliminary to establish temperature-wet factor set (Ti, Hi), wherein 0<i≤n;
(2) duration to be analyzed waved point Q and wave event each time is calculated:
Li=Ei-Si
In formula, LiRepresent the duration to be analyzed waved point Q and wave event each time;
By history wave duration more than a certain setting value wave event retain, and by history wave duration be less than or Event of waving equal to a certain setting value is rejected, and is obtained m times by screening and is waved event, m≤n, and establishes new temperature-wet accordingly Factor set (Tj, Hj), wherein 0<j≤m;
(3) the temperature maximum T obtained by screening is chosenmax, temperature minimum value TminAnd relative humidity maximum Hmax, phase To humidity minimum value Hmin,
Pass through below equation:
<mrow> <msub> <mi>t</mi> <mi>j</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>min</mi> </msub> </mrow> <mrow> <msub> <mi>T</mi> <mi>max</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>min</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>
<mrow> <msub> <mi>h</mi> <mi>j</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>H</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>min</mi> </msub> </mrow> <mrow> <msub> <mi>H</mi> <mi>max</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>min</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow> 2
Simultaneously temperature described in shift step (2)-wet factor set is calculated, obtains temperature-wet factor set (t of specificationj, hj), wherein 0<j≤ m;
Angle threshold is calculated in the second unit and specifically includes following steps:
(1) using temperature as abscissa, relative humidity is ordinate, by temperature-wet factor set (tj, hj) it is plotted in two-dimensional coordinate space In, calculate temperature-wet factor set (tj, hj) in the angle theta of each temperature-wet factor pair and abscissaj
<mrow> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>h</mi> <mi>j</mi> </msub> <msub> <mi>t</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
Calculate θjAverage value θmeanWith standard deviation θstd
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>m</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>;</mo> </mrow>
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>s</mi> <mi>t</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mi>m</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>;</mo> </mrow>
According to θjAverage value θmeanWith standard deviation θstdCalculate the maximum of angle threshold and the minimum value of angle threshold:
θmaxmeanstd
θminmeanstd
It is [θ to obtain angle threshold scopemin, θmax];
The third unit monitors in real time to be analyzed waves the current temperature T of point QobsWith relative humidity Hobs, by above-mentioned formula Calculate and conversion obtain it is to be analyzed wave the current temperature of point Q-wet factor set and therefrom choose it is to be analyzed wave the current temperature of point Q- Wet factor pair (tobs, hobs), angle, θ is calculatedobs, by θobsIt is compared with above-mentioned angle threshold scope, works as θobs∈ [θmin, θmax] when, judge that the point Q to be analyzed that waves can be waved in following certain time;When When, judge that the point Q to be analyzed that waves will not be waved in following certain time.
CN201710640834.2A 2017-07-31 2017-07-31 A kind of power network line waves Occurrence forecast method and system Active CN107436162B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710640834.2A CN107436162B (en) 2017-07-31 2017-07-31 A kind of power network line waves Occurrence forecast method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710640834.2A CN107436162B (en) 2017-07-31 2017-07-31 A kind of power network line waves Occurrence forecast method and system

Publications (2)

Publication Number Publication Date
CN107436162A true CN107436162A (en) 2017-12-05
CN107436162B CN107436162B (en) 2019-07-12

Family

ID=60460274

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710640834.2A Active CN107436162B (en) 2017-07-31 2017-07-31 A kind of power network line waves Occurrence forecast method and system

Country Status (1)

Country Link
CN (1) CN107436162B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107464024A (en) * 2017-08-17 2017-12-12 国网湖南省电力公司 Overhead transmission line galloping Forecasting Methodology and system based on the experiment of icing shape
CN107958312A (en) * 2017-12-12 2018-04-24 国网湖南省电力有限公司 Transmission line galloping Forecasting Methodology, system and storage medium based on inversion algorithm
CN109446474A (en) * 2018-09-05 2019-03-08 中国电力科学研究院有限公司 A kind of methods of exhibiting and system of waving property of data
CN109948839A (en) * 2019-03-06 2019-06-28 国网河南省电力公司电力科学研究院 The prediction of overhead transmission line galloping risk and method for early warning and system
CN110929808A (en) * 2019-12-11 2020-03-27 国网湖南省电力有限公司 Multi-element intelligent correction method and system for waving temperature
CN112365057A (en) * 2020-11-12 2021-02-12 国网湖南省电力有限公司 Regional power transmission line galloping event model construction method and galloping event prediction method
CN113536568A (en) * 2021-07-15 2021-10-22 国网湖南省电力有限公司 Galloping prediction method and system based on analysis of sensitivity of key influence factors of galloping of north-south power transmission line

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070145828A1 (en) * 2005-12-23 2007-06-28 Taylor George F Smart power switch for broadband communications network
KR101102496B1 (en) * 2010-10-18 2012-01-05 (주)인디스디앤아이 Method for deriving state determination equation, apparatus and method for diagnosis of power equipment using the same
CN103559557A (en) * 2013-11-01 2014-02-05 国家电网公司 Gallop warning method and system based on electric transmission line of electrical power system
CN104239970A (en) * 2014-09-04 2014-12-24 国家电网公司 Power transmission line gallop risk early-warning method based on Adaboost
CN104422480A (en) * 2013-08-28 2015-03-18 国家电网公司 Power transmission line online monitoring circuit and intelligent insulator suspension ring

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070145828A1 (en) * 2005-12-23 2007-06-28 Taylor George F Smart power switch for broadband communications network
KR101102496B1 (en) * 2010-10-18 2012-01-05 (주)인디스디앤아이 Method for deriving state determination equation, apparatus and method for diagnosis of power equipment using the same
CN104422480A (en) * 2013-08-28 2015-03-18 国家电网公司 Power transmission line online monitoring circuit and intelligent insulator suspension ring
CN103559557A (en) * 2013-11-01 2014-02-05 国家电网公司 Gallop warning method and system based on electric transmission line of electrical power system
CN104239970A (en) * 2014-09-04 2014-12-24 国家电网公司 Power transmission line gallop risk early-warning method based on Adaboost

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107464024A (en) * 2017-08-17 2017-12-12 国网湖南省电力公司 Overhead transmission line galloping Forecasting Methodology and system based on the experiment of icing shape
CN107958312A (en) * 2017-12-12 2018-04-24 国网湖南省电力有限公司 Transmission line galloping Forecasting Methodology, system and storage medium based on inversion algorithm
CN107958312B (en) * 2017-12-12 2021-08-31 国网湖南省电力有限公司 Electric transmission line galloping prediction method and system based on inversion algorithm and storage medium
CN109446474A (en) * 2018-09-05 2019-03-08 中国电力科学研究院有限公司 A kind of methods of exhibiting and system of waving property of data
CN109948839A (en) * 2019-03-06 2019-06-28 国网河南省电力公司电力科学研究院 The prediction of overhead transmission line galloping risk and method for early warning and system
CN109948839B (en) * 2019-03-06 2023-02-03 国网河南省电力公司电力科学研究院 Method and system for predicting and early warning galloping risk of overhead transmission line
CN110929808A (en) * 2019-12-11 2020-03-27 国网湖南省电力有限公司 Multi-element intelligent correction method and system for waving temperature
CN112365057A (en) * 2020-11-12 2021-02-12 国网湖南省电力有限公司 Regional power transmission line galloping event model construction method and galloping event prediction method
CN113536568A (en) * 2021-07-15 2021-10-22 国网湖南省电力有限公司 Galloping prediction method and system based on analysis of sensitivity of key influence factors of galloping of north-south power transmission line
CN113536568B (en) * 2021-07-15 2022-06-24 国网湖南省电力有限公司 Galloping prediction method and system based on key influence factors of power transmission line galloping

Also Published As

Publication number Publication date
CN107436162B (en) 2019-07-12

Similar Documents

Publication Publication Date Title
CN107436162A (en) A kind of power network line waves Occurrence forecast method and system
CN110035090B (en) False data injection attack detection method for smart grid
CN103713237B (en) A kind of power system transmission line short trouble diagnostic method
CN105353268B (en) One kind is used for the judgement of transmission line of electricity distribution traveling wave fault and localization method
CN107153150A (en) A kind of power distribution network over-voltage fault type recognition method and device
CN105425107B (en) A kind of method and its system of active power distribution network fault diagnosis and location
CN104155519B (en) Harmonic source positioning method with harmonic wave sequence relation analysis and electric power calculation being combined
CN103198184A (en) Low-frequency oscillation character noise-like identification method in electric power system
CN102931728B (en) Online identification and visualization method for power grid disturbances based on multi-resolution wavelet analysis
CN105004939A (en) Composite electric energy quality disturbance signal quantitative analysis method
CN110542879B (en) Method and system for predicting operation performance variation trend of capacitor voltage transformer
CN109858140A (en) One kind being based on comentropy discrete type Bayesian network water cooler method for diagnosing faults
CN101741641A (en) Method for reliability test of communication network services based on link circuits
CN105354768A (en) Assessment method and system for electric transmission and transformation equipment states
CN105764137A (en) Indoor location method and system
CN111209535B (en) Power equipment successive fault risk identification method and system
CN116050599A (en) Line icing fault prediction method, system, storage medium and equipment
CN113644741B (en) Massive fault recording data processing method based on edge calculation
CN106651206A (en) Method for evaluating testability evaluation index system of relay protection
CN113189513A (en) Ripple-based redundant power supply current sharing state identification method
CN113285471A (en) Method, device and equipment for sensing and positioning sub-supersynchronous oscillation source of offshore wind power plant
CN108646145B (en) Method and system for positioning flashover tower of power transmission line
CN115829168B (en) New energy equipment remote coordination operation and maintenance management system based on data analysis
CN109946640B (en) Special transformer terminal voltage and current loss judgment method based on metering automation system data
CN115912625A (en) Online monitoring system and method for power transmission line

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
CB02 Change of applicant information

Address after: 410000 Shaoshan North Road, Yuhua District, Hunan, No. 388, No.

Applicant after: State Grid Hunan Electric Power Co., Ltd.

Applicant after: National network Hunan Electric Power Co., Ltd. disaster prevention and Reduction Center

Applicant after: Hunan Xiangdian Electric Power Text & Research Co., Ltd.

Applicant after: State Grid Corporation of China

Address before: 410000 Shaoshan North Road, Yuhua District, Hunan, No. 388, No.

Applicant before: State Grid Hunan Electric Power Company

Applicant before: State Grid Hunan Electric Power Company prevents and reduces natural disasters center

Applicant before: Hunan Xiangdian Electric Power Text & Research Co., Ltd.

Applicant before: State Grid Corporation of China

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant