CN105095670A - Drawing method of electric transmission line icing-thickness growth rate distribution diagram based on number of continuous freezing rain days - Google Patents

Drawing method of electric transmission line icing-thickness growth rate distribution diagram based on number of continuous freezing rain days Download PDF

Info

Publication number
CN105095670A
CN105095670A CN201510519026.1A CN201510519026A CN105095670A CN 105095670 A CN105095670 A CN 105095670A CN 201510519026 A CN201510519026 A CN 201510519026A CN 105095670 A CN105095670 A CN 105095670A
Authority
CN
China
Prior art keywords
mrow
msub
icing
growth rate
ice coating
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
CN201510519026.1A
Other languages
Chinese (zh)
Other versions
CN105095670B (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
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
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 filed Critical State Grid Corp of China SGCC
Priority to CN201510519026.1A priority Critical patent/CN105095670B/en
Publication of CN105095670A publication Critical patent/CN105095670A/en
Application granted granted Critical
Publication of CN105095670B publication Critical patent/CN105095670B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a drawing method of an electric transmission line icing-thickness growth rate distribution diagram based on number of continuous freezing rain days. The method comprises the steps that 1, an icing database is established; 2, the icing-thickness is preliminarily estimated; 3, an icing-thickness optimization calibration model is established, and an icing-thickness calibration value is calculated when historical continuous freezing rain happens; 4, the icing growth rate is calculated; 5, the recurrence interval of the icing growth rate is calculated; 6, interpolation, smoothness and electric system GIS platform display are conducted on calculated results. The drawing method of the electric transmission line icing-thickness growth rate distribution diagram based on the number of the continuous freezing rain days has the advantages that the operability is strong, and the practicability is strong; the applicability is strong, the portability is good, and the icing-thickness growth rate distribution diagram of any region can be drawn; the icing-thickness growth situation of the historical icing process can be reflected scientifically and reasonably, and an important breakthrough point is provided for studying the historical icing growth rule.

Description

Method for drawing distribution diagram of ice coating thickness growth rate of power transmission line based on continuous rime days
Technical Field
The invention belongs to the technical field of power transmission and distribution, and particularly relates to a method for drawing a distribution map of the ice coating thickness growth rate of a power transmission line based on the number of continuous rimes days.
Background
Icing of a power transmission line is one of main disasters for safe and stable operation of a large power grid. When the line is seriously coated with ice, especially when the ice coating of the power transmission line is rapidly increased in a short time, the risk probability of power grid accidents such as tower falling, line breaking and the like is greatly improved.
In the existing technology for drawing the icing distribution map, the icing thickness is obtained mainly by using a test formula or an empirical formula, then the icing thickness reappearing period is calculated, the icing thickness under different reappearing periods is obtained, and the spatial distribution of the icing thickness of a specified area is drawn. However, in the actual operation process of the power grid, the ice coating growth speed of the areas with the same ice thickness on the ice area map is greatly different, and the damage degree to the power grid is obviously different. For example, in 2014, the icing thickness of the power transmission line in the area of the Yueyang micro-terrain in Hunan in 8-10 days in 2 months, and the average icing thickness growth rate reaches 5 mm/day, so that the power transmission line is tripped. In addition, most of the icing thickness calculation formulas adopted by the ice area map are experimental formulas or empirical formulas, the application range of the formulas is limited, the situation that the actual icing thickness of the power grid is increased is considered to be less, and the applicability of the calculation result is not strong.
Therefore, a distribution diagram for the ice coating thickness growth rate of the power transmission line needs to be drawn, so that scientific and effective guidance is provided for deploying ice-resisting measures in advance for a power grid and making ice-resisting decisions, and the safe and stable operation of the power transmission line is guaranteed.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: aiming at the defects of the prior art, the method for drawing the distribution map of the ice coating thickness growth rate of the power transmission line based on the number of continuous rime days is provided, scientific, reliable and comprehensive guidance can be provided for ice-resistant deployment of a power grid by using the method, and the method is novel in thought, clear in flow, high in accuracy and strong in practicability.
The solution of the invention is:
a method for drawing a distribution map of the ice coating thickness growth rate of a power transmission line based on the number of continuous rimes days comprises the following steps:
(1) establishing an icing database:
establishing a weather icing database according to weather element data of each weather observation station of a weather department, wherein the weather element data comprises the ice coating thickness of an electric wire, observation records of the daily rime phenomenon in winter (12 months-2 months next year), wind speed, daily total precipitation and continuous precipitation duration;
establishing an electric power icing observation database according to meteorological record data of each manual observation whistle of an electric power department, wherein the meteorological record data comprise the icing thickness of a power transmission line, the rime phenomenon observation record, the wind speed, the total daily precipitation and the continuous precipitation duration;
wherein not every site has corresponding line icing thickness data;
(2) preliminarily estimating the accumulated ice coating thickness of each continuous rime occurrence time period (namely the time period of the rime phenomenon occurring for continuous multiple days) of each station, wherein the station comprises a weather observation station and an artificial observation whistle;
counting out through the meteorological icing database and the electric power icing observation database established in the step (1)The period of continuous rime occurrence in winter (12 months-2 months of the next year) and the duration of each period of continuous rime occurrence in the history each year; calculating a preliminary estimate T of the accumulated ice coating thickness of the jth continuous rime occurrence period of the ith year according to an ice coating thickness growth formula by using data in a weather ice coating database and a power ice coating observation databaseij(i 1, 2.. times, n), (j 1, 2.. times, m), wherein n is the total historical years and m is the number of continuous rime occurrence periods in the ith year;
(3) determining an ice coating thickness regression correction model, and calculating a correction value of accumulated ice coating thickness of each continuous rime occurrence period;
(3.1) establishing an icing thickness regression correction modelWherein,a correction value representing the accumulated ice coating thickness of the jth continuous rime occurrence period in the ith year, wherein a and b are regression parameters;
(3.2) taking the site recorded with the icing thickness of the power transmission line as a sample, and calculating regression parameters a and b by using the initial estimation values of the icing thickness of the power transmission line of the sample in the step (1) and the accumulated icing thickness of the sample in the step (2) by adopting a least square method;
(3.3) calculating a correction value of the accumulated ice coating thickness of the jth continuous rime occurrence period of the ith year y ^ i j , ( i = 1 , 2 , ... , n ) , ( j = 1 , 2 , ... , m ) ;
(4) Calculating the average ice coating growth rate;
(4.1) calculated according to the step (3.3)And the duration of the continuous rime occurrence period, and calculating the average ice coating growth rate x of the jth continuous rime occurrence period in the ith yearij(i 1, 2.. times, n), (j 1, 2.. times, m), the calculation formula is:
(4.2) calculating the maximum average icing thickness growth rate x for year ii(i ═ 1, 2.. times, n), the calculation formula is:
xi=max(xi1,xi2,…,xij,…,xim)
(5) calculating the icing growth rate reappearing period;
calculating the ice coating thickness growth rate values of the single site in different reappearance periods by utilizing a probability distribution function according to the annual maximum average ice coating thickness growth rate of the single site calculated in the step (4);
sequentially calculating the average ice coating thickness growth rate of all the sites in different reappearance periods;
(6) interpolating and smoothing the calculation result obtained in the step (5) and displaying the calculation result in the power grid GIS;
determining the latitude and longitude ranges of the power grid GIS display space according to the requirements of the display area, and determining the spatial interpolation resolution; according to the resolution, equally dividing the space range into a plurality of grids;
and (5) interpolating the average ice coating growth rate of all the sites in the display area obtained by calculation in the step (5) in different reappearance periods to the divided grid points by using an interpolation method, smoothing interpolation results, displaying images on a power grid GIS platform, and outputting a distribution diagram drawing result.
Further, in the step (2), a JONES model is adopted as an ice coating thickness increasing formula, and the formula is as follows:
<math> <mrow> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>&pi;</mi> </mrow> </mfrac> <munderover> <mo>&Integral;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>h</mi> </munderover> <msup> <mrow> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>p&rho;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mn>3.6</mn> <mi>v</mi> <mi>w</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&rsqb;</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mi>d</mi> <mi>t</mi> </mrow> </math>
where ρ is1The density of ice is 0.8g/cm3,ρ0The density of water is 1g/cm3P is the total daily precipitation, w is the liquid water content, and is a function of the total daily precipitation p, w is 0.067p0.846(ii) a v is wind speed, h isContinuous precipitation days.
Further, in the step (3.2), the stations are classified according to the altitude in consideration of differences of icing in different altitude areas, regression parameters a and b are respectively calculated, and an icing thickness regression correction model of the stations in a corresponding altitude range is established.
Further, in the step (5), the recurrence period is calculated by using a pearson-type III probability distribution function as follows:
<math> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>&beta;</mi> <mi>&alpha;</mi> </msup> <mrow> <mi>&Gamma;</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mrow> <mi>&alpha;</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>&beta;</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> </msup> </mrow> </math>
<math> <mrow> <mi>&alpha;</mi> <mo>=</mo> <mfrac> <mn>4</mn> <msubsup> <mi>C</mi> <mi>s</mi> <mn>2</mn> </msubsup> </mfrac> <mo>,</mo> <mi>&beta;</mi> <mo>=</mo> <mfrac> <mn>2</mn> <mrow> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <msub> <mi>C</mi> <mi>s</mi> </msub> <msub> <mi>C</mi> <mi>v</mi> </msub> </mrow> </mfrac> <mo>,</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <msub> <mi>C</mi> <mi>v</mi> </msub> <mo>/</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mfrac> </mrow> </math>
<math> <mrow> <msub> <mi>C</mi> <mi>v</mi> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>/</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mrow> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>3</mn> <mo>)</mo> <msubsup> <mi>C</mi> <mi>v</mi> <mn>3</mn> </msubsup> </mrow> </mfrac> </mrow> </math>
wherein, p (x)i) The maximum average icing thickness growth rate x is expressed as a function of probability densityiIs (a) a gamma function,is xiAverage value of (d);
the maximum ice coating growth rate of the single station per year calculated in the step (4) is brought into a formula to obtain ice coating thickness growth rate curves under different probability densities, and the ice coating thickness growth rate curves are obtained according to the probability p (x)i) Relation to the recurrence period TThe average ice accretion rate of the ice coating at different reoccurrence periods of the station was obtained.
Has the advantages that:
the method considers the growth condition of the icing thickness of the power transmission line under the condition of continuous rime, constructs data of the historical icing thickness growth rate, calculates the recurrence period of the icing thickness growth rate, interpolates the calculation results of different recurrence periods onto grid points by using an interpolation method, and locally adjusts the distribution diagram according to actual operation experience to obtain the distribution diagram of the icing thickness growth rate of the power transmission line with strong actual applicability. Has the following advantages:
1. the method has the advantages of clear conditioning, strong operability and strong practicability, and fills the blank that no ice coating growth rate distribution diagram drawing technology exists at present.
2. The method has wide applicability and good portability, and can be used for drawing the icing thickness growth rate distribution map of any region.
3. The calculation method and the correction method of the historical icing growth rate can scientifically and reasonably reflect the growth condition of the icing thickness in the historical icing process, and provide an important breakthrough point for researching the historical icing growth rule;
4. the distribution diagram of the icing growth rate of the power transmission line drawn by the invention provides an important reference for scientifically, reasonably, economically and effectively deploying anti-icing measures and making ice-melting decisions, and reduces the loss of the power grid caused by icing.
5. The method is based on meteorological element data of each meteorological observation station of the meteorological department and meteorological record data of each artificial observation post of the power department, and the data of the historical icing thickness increase rate is constructed, so that the observation stations are multiple, and the precision is higher.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The present invention will be described in more detail with reference to the accompanying drawings and embodiments.
Example 1:
taking Hunan as an example, as shown in FIG. 1, the specific implementation process drawn by the Hunan power grid icing growth rate distribution diagram based on the number of days of rime is as follows:
(1) and establishing an icing database.
Establishing a Hunan meteorological icing database by observing and recording the ice coating thickness of the historical electric wire in Hunan, the rime phenomenon in the past years (63 years in 1951 and 2013), precipitation at meteorological stations and wind speed data; the icing occurrence time period, the icing thickness, the rime phenomenon observation record, the daily total precipitation amount, the continuous precipitation duration and the wind speed recorded by the manual icing observation whistle of the power grid in Hunan province in the process of icing every time in 2011-plus 2013 are established as an icing observation database of the power grid in Hunan province.
(2) Preliminarily estimating the accumulated ice coating thickness of each continuous rime occurrence period of each station, wherein the stations comprise a weather observation station and an artificial observation whistle; according to an icing thickness increasing formula based on a physical process, the JONES model is adopted to preliminarily estimate the icing thickness by utilizing the Hunan power grid icing observation database and the meteorological icing database established in the step (1), and the formula is as follows:
<math> <mrow> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>&pi;</mi> </mrow> </mfrac> <munderover> <mo>&Integral;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>h</mi> </munderover> <msup> <mrow> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>p&rho;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mn>3.6</mn> <mi>v</mi> <mi>w</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&rsqb;</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mi>d</mi> <mi>t</mi> </mrow> </math>
wherein, Tij(i 1, 2.. times, n), (j 1, 2.. times, m) is a preliminary estimation value of the accumulated ice coating thickness of the j-th continuous rime occurrence period in the ith year, n is the total historical year, and m is the number of the i-th continuous rime occurrence period; rho1The density of ice is 0.8g/cm3,ρ0The density of water is 1g/cm3P is total daily precipitation, w is 0.067p0.846V is wind speed, and h is continuous precipitation day.
(3) Determining an ice coating thickness regression correction model, and calculating a correction value of accumulated ice coating thickness of each continuous rime occurrence period;
(3.1) establishing an icing thickness regression correction modelWherein,a correction value representing the accumulated ice coating thickness of the jth continuous rime occurrence period in the ith year, wherein a and b are regression parameters;
(3.2) taking the site recorded with the icing thickness of the power transmission line as a sample, and calculating regression parameters a and b by using the initial estimation values of the icing thickness of the power transmission line of the sample in the step (1) and the accumulated icing thickness of the sample in the step (2) by adopting a least square method;
considering the difference of ice coating in different altitudes of Hunan, the icing point in Hunan area is divided into the following parts according to the altitude: 0-400 m, 400-800 m, 800-1200 m, and more than 1200 m. And establishing a regression correction model of the station icing thickness in the range of the corresponding altitude. The model is as follows:
y ^ i j ( 0 - 400 ) = 0.91 T i j + 0.06
y ^ i j ( 400 - 800 ) = 0.94 T i j + 0.33
y ^ i j ( 800 - 1200 ) = 1.11 T i j + 0.65
y ^ i j ( 1200 ) = 1.38 T i j + 0.24
(3.3) calculating a correction value of the accumulated ice coating thickness of the jth continuous rime occurrence period in the ith year according to the model in (3.2)(4) Calculating the average ice coating growth rate;
(4.1) calculated according to the step (3.3)And the duration of the continuous rime occurrence period, and calculating the average ice coating growth rate x of the jth continuous rime occurrence period in the ith yearij(i 1, 2.. times, n), (j 1, 2.. times, m), the calculation formula is:
(4.2) calculating the maximum average icing thickness growth rate x for year ii(i ═ 1, 2.. times, n), the calculation formula is:
xi=max(xi1,xi2,…,xij,…,xim)
(5) and calculating the ice coating growth rate reappearing period.
Calculating the ice coating thickness growth rate values of the single site in different reappearance periods by utilizing a probability distribution function according to the annual maximum average ice coating thickness growth rate of the single site calculated in the step (4);
the pearson-type III probability distribution function is as follows:
<math> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>&beta;</mi> <mi>&alpha;</mi> </msup> <mrow> <mi>&Gamma;</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mrow> <mi>&alpha;</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>&beta;</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> </msup> </mrow> </math>
<math> <mrow> <mi>&alpha;</mi> <mo>=</mo> <mfrac> <mn>4</mn> <msubsup> <mi>C</mi> <mi>s</mi> <mn>2</mn> </msubsup> </mfrac> <mo>,</mo> <mi>&beta;</mi> <mo>=</mo> <mfrac> <mn>2</mn> <mrow> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <msub> <mi>C</mi> <mi>s</mi> </msub> <msub> <mi>C</mi> <mi>v</mi> </msub> </mrow> </mfrac> <mo>,</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <msub> <mi>C</mi> <mi>v</mi> </msub> <mo>/</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mfrac> </mrow> </math>
<math> <mrow> <msub> <mi>C</mi> <mi>v</mi> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>/</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mrow> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>3</mn> <mo>)</mo> <msubsup> <mi>C</mi> <mi>v</mi> <mn>3</mn> </msubsup> </mrow> </mfrac> </mrow> </math>
wherein, p (x)i) The maximum average icing thickness growth rate x is expressed as a function of probability densityiIs (a) a gamma function,is xiAverage value of (d);
and (4) substituting the historical maximum ice coating growth rate per year calculated in the step (4) into a formula to obtain ice coating thickness growth rate curves under different probability densities, and obtaining average ice coating thickness growth rates of the site with the recurrence periods of 15 years, 30 years, 50 years and 100 years according to the curves. Taking the Yueyang station in Hunan as an example, the ice coating thickness growth rate in different reappearance periods is as follows:
TABLE 1 growth rate of icing thickness at different reproduction periods of Yueyang station
Frequency of Meet once in 15 years Meet once in 30 years Meet one meeting in 50 years Meet every 100 years
Icing growth rate (mm/day) 3.5 4.2 4.8 5.5
And performing corresponding reappearance period calculation on all the sites in Hunan.
(6) And (5) interpolation and smoothing of the calculation result and display of a power grid GIS platform.
And the Hunan power grid GIS shows the spatial latitude range: 24 ° N-31 ° N, longitude space range: 108E-115E, and setting the spatial interpolation resolution to be 1 km. The spatial range is equally divided into a grid of 1km x 1km, depending on the resolution. And (4) interpolating the icing thickness growth rate of all observation sites in Hunan province in different reappearance periods obtained by the calculation in the step (5) to a grid point by using a creatan interpolation method, displaying images on a power grid GIS platform, and outputting a distribution diagram drawing result.

Claims (4)

1. A method for drawing a distribution map of the ice coating thickness growth rate of a power transmission line based on the number of continuous rimes days is characterized by comprising the following steps:
(1) establishing an icing database:
establishing a weather icing database according to weather element data of each weather observation station of a weather department, wherein the weather element data comprises the ice coating thickness of an electric wire, observation records of the daily rime phenomenon in winter (12 months-2 months next year), wind speed, daily total precipitation and continuous precipitation duration;
establishing an electric power icing observation database according to meteorological record data of each manual observation whistle of an electric power department, wherein the meteorological record data comprise the icing thickness of a power transmission line, the rime phenomenon observation record, the wind speed, the total daily precipitation and the continuous precipitation duration;
(2) preliminarily estimating the accumulated ice coating thickness of each continuous rime occurrence period of each station, wherein the stations comprise a weather observation station and an artificial observation whistle;
counting the historical continuous rime occurrence time period in winter (12 months-next 2 months) every year and the time length of each continuous rime occurrence time period through the weather icing database and the electric power icing observation database established in the step (1); calculating a preliminary estimate T of the accumulated ice coating thickness of the jth continuous rime occurrence period of the ith year according to an ice coating thickness growth formula by using data in a weather ice coating database and a power ice coating observation databaseij(i 1, 2.. times, n), (j 1, 2.. times, m), wherein n is the total historical years and m is the number of continuous rime occurrence periods in the ith year;
(3) determining an ice coating thickness regression correction model, and calculating a correction value of accumulated ice coating thickness of each continuous rime occurrence period;
(3.1) establishing an icing thickness regression correction modelWherein,a correction value representing the accumulated ice coating thickness of the jth continuous rime occurrence period in the ith year, wherein a and b are regression parameters;
(3.2) taking the site recorded with the icing thickness of the power transmission line as a sample, and calculating regression parameters a and b by using the initial estimation values of the icing thickness of the power transmission line of the sample in the step (1) and the accumulated icing thickness of the sample in the step (2) by adopting a least square method;
(3.3) calculating a correction value of the accumulated ice coating thickness of the jth continuous rime occurrence period of the ith year y ^ i j , ( i = 1 , 2 , ... , n ) , ( j = 1 , 2 , ... , m ) ;
(4) Calculating the average ice coating growth rate;
(4.1) calculated according to the step (3.3)And the duration of the continuous rime occurrence period, and calculating the average ice coating growth rate x of the jth continuous rime occurrence period in the ith yearij(i 1, 2.. times, n), (j 1, 2.. times, m), the calculation formula is:
(4.2) calculating the maximum average icing thickness growth rate x for year ii(i ═ 1, 2.. times, n), the calculation formula is:
xi=max(xi1,xi2,…,xij,…,xim)
(5) calculating the icing growth rate reappearing period;
calculating the ice coating thickness growth rate values of the single site in different reappearance periods by utilizing a probability distribution function according to the annual maximum average ice coating thickness growth rate of the single site calculated in the step (4);
sequentially calculating the average ice coating thickness growth rate of all the sites in different reappearance periods;
(6) interpolating and smoothing the calculation result obtained in the step (5) and displaying the calculation result in the power grid GIS;
determining the latitude and longitude ranges of the power grid GIS display space according to the requirements of the display area, and determining the spatial interpolation resolution; according to the resolution, equally dividing the space range into a plurality of grids;
and (5) interpolating the average ice coating growth rate of all the sites in the display area obtained by calculation in the step (5) in different reappearance periods to the divided grid points by using an interpolation method, smoothing interpolation results, displaying images on a power grid GIS platform, and outputting a distribution diagram drawing result.
2. The method for drawing the distribution map of the ice coating thickness growth rate of the power transmission line based on the number of continuous rimes days according to claim 1, wherein in the step (2), a JONES model is adopted as an ice coating thickness growth formula, and the formula is as follows:
<math> <mrow> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>&pi;</mi> </mrow> </mfrac> <munderover> <mo>&Integral;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>h</mi> </munderover> <msup> <mrow> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>p&rho;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mn>3.6</mn> <mi>v</mi> <mi>w</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&rsqb;</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mi>d</mi> <mi>t</mi> </mrow> </math>
where ρ is1Is iceDensity of 0.8g/cm3,ρ0The density of water is 1g/cm3P is the total daily precipitation, w is the liquid water content, and is a function of the total daily precipitation p, w is 0.067p0.846(ii) a v is wind speed and h is continuous precipitation day.
3. The method for drawing the distribution map of the ice coating thickness growth rate of the power transmission line based on the number of continuous rime days according to claim 2, wherein in the step (3.2), the stations are classified according to the altitude, regression parameters a and b are respectively calculated, and a regression correction model of the ice coating thickness of the stations in the corresponding altitude range is established.
4. The method for drawing the distribution map of the growth rate of the ice coating thickness of the power transmission line based on the number of days of continuous rime according to any one of claims 1 to 3, wherein in the step (5), the calculation of the recurrence period is performed by using a Pearson-III type probability distribution function, which is as follows:
<math> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>&beta;</mi> <mi>&alpha;</mi> </msup> <mrow> <mi>&Gamma;</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mrow> <mi>&alpha;</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>&beta;</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mover> <mn>1</mn> <mo>&CenterDot;</mo> </mover> </msub> <mo>-</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> </msup> </mrow> </math>
<math> <mrow> <mi>&alpha;</mi> <mo>=</mo> <mfrac> <mn>4</mn> <msubsup> <mi>C</mi> <mi>s</mi> <mn>2</mn> </msubsup> </mfrac> <mo>,</mo> <mi>&beta;</mi> <mo>=</mo> <mfrac> <mn>2</mn> <mrow> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <msub> <mi>C</mi> <mi>s</mi> </msub> <msub> <mi>C</mi> <mi>v</mi> </msub> </mrow> </mfrac> <mo>,</mo> <msub> <mi>&alpha;</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <msub> <mi>C</mi> <mi>v</mi> </msub> <mo>/</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mfrac> </mrow> </math>
<math> <mrow> <msub> <mi>C</mi> <mi>v</mi> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>/</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mrow> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>3</mn> <mo>)</mo> <msubsup> <mi>C</mi> <mi>v</mi> <mn>3</mn> </msubsup> </mrow> </mfrac> </mrow> </math>
wherein, p (x)i) The maximum average icing thickness growth rate x is expressed as a function of probability densityiIs (a) a gamma function,is xiAverage value of (d);
the maximum ice coating growth rate of the single station per year calculated in the step (4) is brought into a formula to obtain ice coating thickness growth rate curves under different probability densities, and the ice coating thickness growth rate curves are obtained according to the probability p (x)i) Relation to the recurrence period TThe average ice accretion rate of the ice coating at different reoccurrence periods of the station was obtained.
CN201510519026.1A 2015-08-21 2015-08-21 Drawing method of electric transmission line icing-thickness growth rate distribution diagram based on continuous freezing rain Active CN105095670B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510519026.1A CN105095670B (en) 2015-08-21 2015-08-21 Drawing method of electric transmission line icing-thickness growth rate distribution diagram based on continuous freezing rain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510519026.1A CN105095670B (en) 2015-08-21 2015-08-21 Drawing method of electric transmission line icing-thickness growth rate distribution diagram based on continuous freezing rain

Publications (2)

Publication Number Publication Date
CN105095670A true CN105095670A (en) 2015-11-25
CN105095670B CN105095670B (en) 2017-03-22

Family

ID=54576088

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510519026.1A Active CN105095670B (en) 2015-08-21 2015-08-21 Drawing method of electric transmission line icing-thickness growth rate distribution diagram based on continuous freezing rain

Country Status (1)

Country Link
CN (1) CN105095670B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105808819A (en) * 2016-02-17 2016-07-27 南方电网科学研究院有限责任公司 Calculation method for real-time distribution model of icing of power transmission line
CN106595551A (en) * 2016-12-13 2017-04-26 武汉大学 Icing thickness detection method for power transmission line icing image based on deep learning
CN107436972A (en) * 2017-07-07 2017-12-05 国网湖南省电力公司 Powerline ice-covering speedup variation characteristic computational methods and system under special climate
CN107464024A (en) * 2017-08-17 2017-12-12 国网湖南省电力公司 Overhead transmission line galloping Forecasting Methodology and system based on the experiment of icing shape
CN112116132A (en) * 2020-08-27 2020-12-22 中国南方电网有限责任公司超高压输电公司检修试验中心 Method for predicting conductor icing thickness in complex mountain environment

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110649505B (en) * 2019-09-09 2020-12-11 国网北京市电力公司 Three-dimensional monitoring method and device for icing of power transmission line

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU556358A1 (en) * 1975-12-01 1977-04-30 Московский институт инженеров сельскохозяйственного производства им.В.П.Горячкина Device for measuring and recording ice loads on power lines
CN101672667A (en) * 2008-09-10 2010-03-17 华东电力试验研究院有限公司 Method for forecasting wire icing of transmission line
CN104462660A (en) * 2014-11-14 2015-03-25 贵州电力试验研究院 Drawing method for winter icing thickness distribution of field electric transmission line
CN104636601A (en) * 2015-01-05 2015-05-20 中国南方电网有限责任公司 Electric transmission line icing thickness forecasting method for introducing dynamic de-icing process

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU556358A1 (en) * 1975-12-01 1977-04-30 Московский институт инженеров сельскохозяйственного производства им.В.П.Горячкина Device for measuring and recording ice loads on power lines
CN101672667A (en) * 2008-09-10 2010-03-17 华东电力试验研究院有限公司 Method for forecasting wire icing of transmission line
CN104462660A (en) * 2014-11-14 2015-03-25 贵州电力试验研究院 Drawing method for winter icing thickness distribution of field electric transmission line
CN104636601A (en) * 2015-01-05 2015-05-20 中国南方电网有限责任公司 Electric transmission line icing thickness forecasting method for introducing dynamic de-icing process

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘浩梁等: "微气象条件下架空地线的覆冰增长", 《电力建设》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105808819A (en) * 2016-02-17 2016-07-27 南方电网科学研究院有限责任公司 Calculation method for real-time distribution model of icing of power transmission line
CN105808819B (en) * 2016-02-17 2017-12-26 南方电网科学研究院有限责任公司 Calculation method for real-time distribution model of icing of power transmission line
CN106595551A (en) * 2016-12-13 2017-04-26 武汉大学 Icing thickness detection method for power transmission line icing image based on deep learning
CN106595551B (en) * 2016-12-13 2019-01-04 武汉大学 Ice covering thickness detection method in powerline ice-covering image based on deep learning
CN107436972A (en) * 2017-07-07 2017-12-05 国网湖南省电力公司 Powerline ice-covering speedup variation characteristic computational methods and system under special climate
CN107436972B (en) * 2017-07-07 2020-11-06 国网湖南省电力有限公司 Method and system for calculating icing and acceleration change characteristics of power transmission line in special climate
CN107464024A (en) * 2017-08-17 2017-12-12 国网湖南省电力公司 Overhead transmission line galloping Forecasting Methodology and system based on the experiment of icing shape
CN112116132A (en) * 2020-08-27 2020-12-22 中国南方电网有限责任公司超高压输电公司检修试验中心 Method for predicting conductor icing thickness in complex mountain environment
CN112116132B (en) * 2020-08-27 2023-08-01 中国南方电网有限责任公司超高压输电公司检修试验中心 Prediction method for wire icing thickness in complex mountain area

Also Published As

Publication number Publication date
CN105095670B (en) 2017-03-22

Similar Documents

Publication Publication Date Title
CN105095670B (en) Drawing method of electric transmission line icing-thickness growth rate distribution diagram based on continuous freezing rain
CN107316095B (en) Regional weather drought level prediction method coupled with multi-source data
Noriah et al. Probable maximum precipitation for 24 h duration over southeast Asian monsoon region—Selangor, Malaysia
Bellaire et al. Forecasting the formation of critical snow layers using a coupled snow cover and weather model
Hitchens et al. Spatial and temporal characteristics of heavy hourly rainfall in the United States
CN114019579B (en) High space-time resolution near-surface air temperature reconstruction method, system and equipment
CN105631115A (en) Refined covering ice model establishment method for power transmission line
Zhu et al. A new procedure to estimate the rainfall erosivity factor based on Tropical Rainfall Measuring Mission (TRMM) data
CN111242404B (en) Extreme evaluation method and system for heavy rainfall induced flood incident
CN105975791B (en) A kind of sparse region Calculation of precipitation fusion method based on double exponential smoothings
Chavan et al. Probable maximum precipitation estimation for catchments in Mahanadi river basin
CN103090831A (en) Judgment method of icing thickness of icing area electric transmission line
CN112686426B (en) Incoming water quantity early warning method and system based on hydropower station basin key points
CN110097223B (en) Early warning method for damage of power transmission line under typhoon disaster
Kolendowicz et al. Atmospheric circulation and sounding-derived parameters associated with thunderstorm occurrence in Central Europe
Chubb et al. Estimation of wind-induced losses from a precipitation gauge network in the Australian snowy mountains
CN107748933B (en) Meteorological element message data error correction method and fog, sunrise, cloud sea and rime prediction method
CN107765348B (en) Precipitation weather message data error correction method and rime landscape prediction method
Sokol et al. Nowcasting of precipitation–advective statistical forecast model (SAM) for the Czech Republic
CN108510098B (en) Power transmission line corridor snow depth estimation and early warning method and system based on satellite remote sensing
Neilsen et al. Development and verification of daily gridded climate surfaces in the Okanagan Basin of British Columbia
KR101870743B1 (en) Tidal flat experience index forecasting apparatus and the method thereof
Feng A new method for ionospheric short-term forecast using similar-day modeling
CN112561140A (en) China four-season rainfall prediction method based on cooperative change of east Asia subtropical zone torrent and extreme torrent
Jones et al. The development of a US climatology of extreme ice loads

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant