[go: up one dir, main page]

CN111539842B - Overhead power transmission line icing prediction method based on meteorological and geographic environments - Google Patents

Overhead power transmission line icing prediction method based on meteorological and geographic environments Download PDF

Info

Publication number
CN111539842B
CN111539842B CN202010267424.XA CN202010267424A CN111539842B CN 111539842 B CN111539842 B CN 111539842B CN 202010267424 A CN202010267424 A CN 202010267424A CN 111539842 B CN111539842 B CN 111539842B
Authority
CN
China
Prior art keywords
data
icing
transmission lines
value
model
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.)
Active
Application number
CN202010267424.XA
Other languages
Chinese (zh)
Other versions
CN111539842A (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.)
Sichuan Shenzhou Lvchuang Smart Energy Co ltd
Original Assignee
Chengdu Sihan Technology 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 Chengdu Sihan Technology Co ltd filed Critical Chengdu Sihan Technology Co ltd
Priority to CN202010267424.XA priority Critical patent/CN111539842B/en
Publication of CN111539842A publication Critical patent/CN111539842A/en
Application granted granted Critical
Publication of CN111539842B publication Critical patent/CN111539842B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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
    • 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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Software Systems (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • General Engineering & Computer Science (AREA)
  • Game Theory and Decision Science (AREA)
  • Computing Systems (AREA)
  • Medical Informatics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明属于覆冰预测技术领域,具体涉及一种基于气象和地理环境的架空输电线路覆冰预测方法。本发明主要是基于气象数据和地理环境数据,利用机器学习方法对未来一段时间的架空输电线路覆冰情况进行预测和预警;根据目标线路的历史覆冰数据训练机器学习模型,并解决类别失衡,大数据训练等问题。整个过程包括数据处理,特征的抽取和覆冰预测模型建立。本发明在海量的数据上,通过数据处理技术,特征抽取技术和分类模型建立,实现架空输电线路的覆冰预测和预警,通过对输电线路所有历史的覆冰数据,结合覆冰的影响因子(气象,地理环境等)实现未来一段时间的覆冰预测。

Figure 202010267424

The invention belongs to the technical field of icing prediction, and in particular relates to a method for icing prediction of overhead transmission lines based on weather and geographical environment. The present invention is mainly based on meteorological data and geographic environment data, using machine learning methods to predict and warn the icing situation of overhead transmission lines in the future; train the machine learning model according to the historical icing data of the target line, and solve the category imbalance, Big data training and other issues. The whole process includes data processing, feature extraction and icing prediction model establishment. The present invention realizes icing prediction and early warning of overhead transmission lines through data processing technology, feature extraction technology and classification model establishment on massive data, and combines all historical icing data of transmission lines with the impact factor of icing ( Meteorology, geographical environment, etc.) to realize ice coverage prediction for a period of time in the future.

Figure 202010267424

Description

基于气象和地理环境的架空输电线路覆冰预测方法Icing Prediction Method for Overhead Transmission Lines Based on Meteorological and Geographical Environment

技术领域technical field

本发明属于覆冰预测技术领域,具体涉及一种基于气象和地理环境的架空输电线路覆冰预测方法。The invention belongs to the technical field of icing prediction, in particular to an icing prediction method for overhead transmission lines based on meteorological and geographical environments.

背景技术Background technique

输电线路覆冰是电力系统最重要的灾害之一,由于覆冰而使输电线路的荷重增加,造成线路断线、杆塔倒塌、冰闪跳闸、导线舞动、通信中断、绝缘子和金具等设备损毁等事故。覆冰灾害会严重威胁到电网的安全稳定运行,造成巨大的经济损失。我国很多地方都存在严重的覆冰灾害,预判覆冰的情况是电网运营的重要一环,因此,关于输电线路的覆冰预测是对电网有非常积极和重要意义的。Icing on transmission lines is one of the most important disasters in power systems. Due to icing, the load on transmission lines increases, causing line disconnection, tower collapse, ice flash tripping, wire galloping, communication interruption, damage to insulators and fittings, etc. ACCIDENT. Icing disasters will seriously threaten the safe and stable operation of the power grid and cause huge economic losses. There are serious icing disasters in many places in our country. Predicting the situation of icing is an important part of power grid operation. Therefore, the prediction of icing on transmission lines is very positive and important to the power grid.

输电线路的覆冰机理较为复杂,目前基于输电线路的覆冰较多的研究是基于物理方程的方法,也有通过动力学,热力学实验来进行研究。这些方法也存在很多问题,比如在应用时的难度较大,应用时的准确率较低。基于物理方程的复杂性和假设性在应用时存在很多问题,覆冰的机理复杂,很难通过一个物理方程进行描述,在通用性方面存在较大问题;另外,一般物理方程依赖很多经验系数,而这些系数很难被确定;而且物理方法无法对未来覆冰情况进行预测和判定。The mechanism of icing on transmission lines is relatively complicated. At present, most studies on icing on transmission lines are based on physical equations, and some studies are conducted through kinetic and thermodynamic experiments. There are also many problems in these methods, such as greater difficulty in application and lower accuracy in application. There are many problems in the application based on the complexity and assumptions of physical equations. The mechanism of icing is complex and difficult to describe through a physical equation, and there are big problems in generality. In addition, general physical equations rely on many empirical coefficients. However, these coefficients are difficult to determine; and physical methods cannot predict and judge future icing conditions.

随着大数据和机器学习的兴起,关于输电线路的覆冰数据越来越多,因此本发明采用基于大数据的机器学习方法来对输电线路的覆冰进行研究,并发明了基于气象和地理环境的架空输电线路覆冰预测方法。而传统的物理方法主要缺点如下:With the rise of big data and machine learning, there are more and more icing data about transmission lines, so the present invention adopts a machine learning method based on big data to study the icing of transmission lines, and invents a method based on weather and geography Environmental icing prediction methods for overhead transmission lines. The main disadvantages of traditional physical methods are as follows:

1、物理方法的主要是基于某一个局部进行观测给出较为合理的物理方程,当在大量应用时通用性会出线较大的问题,必须对每个局部(每个塔体或一段线路)建立物理方程,应用性非常低。1. The physical method is mainly based on the observation of a certain part to give a more reasonable physical equation. When it is used in a large number of applications, the versatility will be a big problem. It must be established for each part (each tower or a section of line) Physical equations, very low applicability.

2、物理方程本身是有很多影响因子和系数构建,但是系数的确定一般根据经验或多次观测过程给出,导致系数的确定较为复杂且难度较大,这些系数的合理性不易科学的验证。覆冰的形成机理极为复杂,在不同的时间点点,可能系数都可能不同。因此整个构建物理方程的过程难度较大。2. The physical equation itself is constructed with many influencing factors and coefficients, but the determination of the coefficients is generally given based on experience or multiple observations, which makes the determination of the coefficients more complicated and difficult, and the rationality of these coefficients is not easy to scientifically verify. The formation mechanism of icing is extremely complex, and the possible coefficients may be different at different points in time. Therefore, the whole process of constructing physical equations is more difficult.

3、物理方法主要是基于当前的情况来评估当前的覆冰状态,无法对未来一段时间进行预测和评估。3. The physical method is mainly based on the current situation to evaluate the current icing state, and cannot predict and evaluate for a period of time in the future.

发明内容Contents of the invention

本发明主要是基于气象数据和地理环境数据,利用机器学习方法对未来一段时间的架空输电线路覆冰情况进行预测和预警;根据目标线路的历史覆冰数据训练机器学习模型,并解决类别失衡,大数据训练等问题。整个过程包括数据处理,特征的抽取和覆冰预测模型建立。The present invention is mainly based on meteorological data and geographical environment data, and uses machine learning methods to predict and warn the icing situation of overhead transmission lines in the future; trains machine learning models according to historical icing data of target lines, and solves category imbalance, Big data training and other issues. The whole process includes data processing, feature extraction and icing prediction model establishment.

为实现上述目的,本发明采用如下技术方案:To achieve the above object, the present invention adopts the following technical solutions:

基于气象和地理环境的架空输电线路覆冰预测方法,如图1所示,包括以下步骤:The icing prediction method for overhead transmission lines based on meteorological and geographical environment, as shown in Figure 1, includes the following steps:

S1、数据收集:收集气象数据和地理环境高程数据,气象数据包括温度、湿度、降雨量、风速、风向、风力等级、日期;S1. Data collection: Collect meteorological data and geographic environment elevation data. Meteorological data include temperature, humidity, rainfall, wind speed, wind direction, wind force level, and date;

S2、数据处理:S2. Data processing:

将格式为hdf格式的高程数据通过hdf进行解析,并将缺失值标记为-1;Parse the elevation data in hdf format through hdf, and mark the missing value as -1;

将获得的温度数据中高于40摄氏度的统一为40摄氏度;Unify the obtained temperature data above 40 degrees Celsius as 40 degrees Celsius;

将获得的湿度小于0的值均设为0,湿度大于100的值均设置为100;Set the obtained values of humidity less than 0 to 0, and the values of humidity greater than 100 to 100;

S3、对收集的数据进行衍生计算,以提取更多的特征,具体为:S3. Perform derivative calculations on the collected data to extract more features, specifically:

通过高程数据进行衍生计算得到坡度和坡向,采用拟合曲面法,对于点e,用其最邻近的8个数据点,左上,上,右上,左,右,右下,下,右下这个8个点加上中心点e够成3*3的窗口(矩阵)(这里有个关键的问题,因为在权利要求里面也采用了这个描述,而权利要求是不能引用附图的,也就是附图不能用来解释权利要求,权利要求必须只能是文字描述,所以才想确定下是否可以这样描述这个窗口),如图2所示,第一列的数据依次为e5、e2、e6,第二列的数据依次为e1、e、e3,第三列的数据依次为e8、e4、e7,则中心点e的坡度Slope和坡向Aspect计算公式如下:The slope and aspect are obtained by deriving the elevation data, using the fitting surface method, for the point e, use its 8 nearest data points, upper left, upper, upper right, left, right, lower right, lower, lower right 8 points plus the center point e can form a 3*3 window (matrix) (here is a key problem, because this description is also used in the claims, and the claims cannot refer to the drawings, that is, the attached The picture cannot be used to explain the claim, the claim must only be described in words, so I want to determine whether this window can be described in this way), as shown in Figure 2, the data in the first column are e 5 , e 2 , e 6. The data in the second column are e 1 , e , and e 3 in sequence, and the data in the third column are e 8 , e 4 , and e 7 in sequence. The calculation formulas for the slope and aspect of the center point e are as follows:

Figure BDA0002441814160000021
Figure BDA0002441814160000021

Figure BDA0002441814160000022
Figure BDA0002441814160000022

其中,

Figure BDA0002441814160000031
cellsize是像元点的分辨率大小,如果点e周围存钱缺失值,则用邻近的值替代,若都缺少,则点e的坡度和坡向均为0;in,
Figure BDA0002441814160000031
cellsize is the resolution of the pixel point. If there is a missing value around point e, it will be replaced by the adjacent value. If both are missing, the slope and aspect of point e are both 0;

对获取的气象数据进行统计计算,包括:Perform statistical calculations on the acquired meteorological data, including:

最大值vmax=max(vt),t=t,t-1,t-2,t-3,t-4,...Maximum value v max =max(v t ), t=t, t-1, t-2, t-3, t-4,  …

最小值vmin=min(vt),t=t,t-1,t-2,t-3,t-4,...Minimum value v min =min(v t ), t=t, t-1, t-2, t-3, t-4,  …

均值

Figure BDA0002441814160000032
average
Figure BDA0002441814160000032

标准差

Figure BDA0002441814160000033
standard deviation
Figure BDA0002441814160000033

极值vp=vmax-vmin Extreme value v p =v max -v min

其中vt是指气象数据,t指最近的t个时间点,获得的气象数据中,根据当前时间获取的当前值,对温度进行最大、最小、均值、标准差、极值的计算;湿度进行最大、最小、均值、极值的计算;风速进行最大值计算;风力等级进行最大值计算;降雨量进行最大、最小、极值的计算;Among them, v t refers to meteorological data, and t refers to the latest t time points. Among the obtained meteorological data, the maximum, minimum, mean, standard deviation, and extreme value of the temperature are calculated according to the current value obtained at the current time; the humidity is calculated Calculation of maximum, minimum, average and extreme values; maximum calculation of wind speed; maximum calculation of wind power level; maximum, minimum and extreme calculation of rainfall;

对于获得的风向数据进行离散化,因风向数据为字符型,将其转化为数值型,离散化为分类变量,即将北,东东北,东北,东北北,东,东东南,东南,东南南,南,西南南,西南,西西南,西,西西北,西北,西北北16个方位转化为1-16的整数数值;Discretize the obtained wind direction data, because the wind direction data is a character type, convert it into a numerical type, and discretize it into a categorical variable, that is, north, east northeast, northeast, northeast north, east, east southeast, southeast, southeast south, South, Southwest South, Southwest, West Southwest, West, West Northwest, Northwest, Northwest North 16 orientations are converted into integer values from 1 to 16;

S4、判断步骤S3获得的数据量是否大于10万,若是,则进入步骤S5,否则进入步骤S6;S4, judging whether the amount of data obtained in step S3 is greater than 100,000, if so, proceed to step S5, otherwise proceed to step S6;

S5、基于架空输电线路覆冰中正样本和负样本之间是否相差5倍的数量级来进行类别失衡的判断,若相差5倍,则判定类别失衡,将数据输入bagging模式模型,否则采用机器学习的方法来获得预测结果;所述bagging模式模型为:将输入的数据集进行分割,获得N个子集,针对每个子集构建一个训练模型,即构建N个训练模型,每个子集经过对应的训练模型后,获得N个训练模型的结果,采用bagging算法对N个模型结果进行集成,获得预测结果;S5. Based on whether there is an order of magnitude difference of 5 times between the positive samples and negative samples in the icing of overhead transmission lines, the category imbalance is judged. If the difference is 5 times, the category imbalance is determined, and the data is input into the bagging mode model. Otherwise, the machine learning method is used. method to obtain prediction results; the bagging mode model is: the input data set is divided to obtain N subsets, and a training model is constructed for each subset, that is, N training models are constructed, and each subset undergoes a corresponding training model Finally, the results of N training models are obtained, and the bagging algorithm is used to integrate the results of N models to obtain prediction results;

S6、基于架空输电线路覆冰中正样本和负样本之间是否相差5倍的数量级来进行类别失衡的判断,若相差5倍,则判定类别失衡,将数据输入权重模型,否则采用机器学习的方法获得预测结果;所述权重模型为:在目标函数训练时取的权重不同,定义目标函数为S6. Based on whether there is an order of magnitude difference of 5 times between the positive samples and negative samples in the icing of overhead transmission lines, the category imbalance is judged. If the difference is 5 times, the category imbalance is determined, and the data is input into the weight model. Otherwise, the machine learning method is used. Obtain prediction result; Described weight model is: the weight that gets when objective function training is different, defines objective function as

Figure BDA0002441814160000041
Figure BDA0002441814160000041

其中k为类别变量,w为权重,L为loss函数,yk为观测样本中的目标变量,f(xk)为基础模型的输出结果,所述基础模型包括逻辑回归、支持向量机、决策树、boosting集成学习分类方法;Where k is the category variable, w is the weight, L is the loss function, y k is the target variable in the observation sample, f(x k ) is the output result of the basic model, and the basic model includes logistic regression, support vector machine, decision-making Tree, boosting integrated learning classification method;

权重w的计算公式为:The formula for calculating the weight w is:

Figure BDA0002441814160000042
Figure BDA0002441814160000042

其中,n为样本数量,mclass为class对应的样本数量,class取值为目标变量的类别;Among them, n is the number of samples, m class is the number of samples corresponding to the class, and the value of class is the category of the target variable;

根据获得的权重w进行训练后,利用训练好的模型获得预测结果。After training according to the obtained weight w, use the trained model to obtain the prediction result.

本发明的有益效果为,本发明在海量的数据上,通过数据处理技术,特征抽取技术和分类模型建立,实现架空输电线路的覆冰预测和预警,通过对输电线路所有历史的覆冰数据,结合覆冰的影响因子(气象,地理环境等)实现未来一段时间的覆冰预测;在大数据的基础上,结合了覆冰的形成机理,采用机器学习的方式对目标区域进行训练和预测。整个策略更加灵活,能有效的对每个目标区域进行预测,且在应用时更具有通用性和便捷性,不仅克服了物理方法的确点,而且本发使得输电线路的覆冰预测更加灵活,简易,易于应用。The beneficial effect of the present invention is that the present invention realizes icing prediction and early warning of overhead transmission lines through data processing technology, feature extraction technology and classification model establishment on massive data, and through all historical icing data of transmission lines, Combining the influencing factors of icing (meteorology, geographical environment, etc.) to realize icing prediction for a period of time in the future; on the basis of big data, combining the formation mechanism of icing, and using machine learning to train and predict the target area. The whole strategy is more flexible, can effectively predict each target area, and is more versatile and convenient in application. It not only overcomes the definite point of physical methods, but also makes the icing prediction of transmission lines more flexible and simple. , easy to apply.

附图说明Description of drawings

图1为本发明的逻辑结构示意图;Fig. 1 is a schematic diagram of the logical structure of the present invention;

图2为基于高程数据计算坡度和坡向建立的3*3数据窗口示意图;Figure 2 is a schematic diagram of a 3*3 data window based on the calculation of slope and aspect based on elevation data;

图3为本发明基于bagging模式构建的模型示意图;Fig. 3 is the model schematic diagram that the present invention builds based on bagging mode;

图4为本发明模型训练的示意图;Fig. 4 is the schematic diagram of model training of the present invention;

图5为本发明的AUC和ROC曲线指标示意图。Fig. 5 is a schematic diagram of AUC and ROC curve indicators of the present invention.

具体实施方式Detailed ways

下面结合附图对本发明进行进一步说明。The present invention will be further described below in conjunction with the accompanying drawings.

如图3所示,本发明采用多个子模型策略构建模型,子模型的基分类器选择了简单的逻辑回归(Logistic)模型,例如将总的数据集划分为41个子集,每个子集构建一个模型,总共构建41个模型。在应用时,数据都会经过41个模型,然后通过投票的方式将41个模型结果进行统一输出。As shown in Figure 3, the present invention adopts a plurality of sub-model strategies to build a model, and the base classifier of the sub-model selects a simple logistic regression (Logistic) model, for example, the total data set is divided into 41 subsets, and each subset constructs a models, a total of 41 models were constructed. During application, the data will pass through 41 models, and then the results of the 41 models will be output uniformly through voting.

由于使用的数据量较大,因此在训练时采用逐批次训练的方式进行。并在训练过程中引入参数优化,输出最优的模型。训练的流程如图4所示。Due to the large amount of data used, batch-by-batch training is adopted during training. And introduce parameter optimization in the training process to output the optimal model. The training process is shown in Figure 4.

模型构建后,基于测试集对训练后的模型进行效果评估,评估指标主要是通过AUC和ROC曲线,本发明的AUC值为0.9665,ROC曲线如图5所示。After the model is built, the effect of the trained model is evaluated based on the test set. The evaluation index is mainly through AUC and ROC curve. The AUC value of the present invention is 0.9665, and the ROC curve is shown in Figure 5.

Claims (1)

1.基于气象和地理环境的架空输电线路覆冰预测方法,其特征在于,包括以下步骤:1. The method for predicting icing of overhead transmission lines based on meteorological and geographical environment, is characterized in that, comprises the following steps: S1、数据收集:收集气象数据和地理环境高程数据,气象数据包括温度、湿度、降雨量、风速、风向、风力等级、日期;S1. Data collection: Collect meteorological data and geographic environment elevation data. Meteorological data include temperature, humidity, rainfall, wind speed, wind direction, wind force level, and date; S2、数据处理:S2. Data processing: 将格式为hdf格式的高程数据通过hdf进行解析,并将缺失值标记为-1;Parse the elevation data in hdf format through hdf, and mark the missing value as -1; 将获得的温度数据中高于40摄氏度的统一为40摄氏度;Unify the obtained temperature data above 40 degrees Celsius as 40 degrees Celsius; 将获得的湿度小于0的值均设为0,湿度大于100的值均设置为100;Set the obtained values of humidity less than 0 to 0, and the values of humidity greater than 100 to 100; S3、对收集的数据进行衍生计算,以提取更多的特征,具体为:S3. Perform derivative calculations on the collected data to extract more features, specifically: 通过高程数据进行衍生计算得到坡度和坡向,采用拟合曲面法,对于点e,用其邻近的8个数据将其包围形成3*3的窗口,第一列的数据依次为e5、e2、e6,第二列的数据依次为e1、e、e3,第三列的数据依次为e8、e4、e7,则中心点e的坡度Slope和坡向Aspect计算公式如下:The slope and aspect are obtained through derivative calculation of elevation data. Using the fitting surface method, for point e, it is surrounded by 8 adjacent data to form a 3*3 window. The data in the first column are e 5 , e 2 , e 6 , the data in the second column are e 1 , e, e 3 in sequence, and the data in the third column are e 8 , e 4 , e 7 in sequence, then the calculation formulas for the slope and aspect of the center point e are as follows :
Figure QLYQS_1
Figure QLYQS_1
Figure QLYQS_2
Figure QLYQS_2
其中,
Figure QLYQS_3
其中cellsize是像元点的分辨率大小,如果点e周围存钱缺失值,则用邻近的值替代,若都缺少,则点e的坡度和坡向均为0;
in,
Figure QLYQS_3
Among them, cellsize is the resolution of the pixel point. If there is a missing value around point e, it will be replaced by an adjacent value. If both are missing, the slope and aspect of point e are both 0;
对获取的气象数据进行统计计算,包括:Perform statistical calculations on the acquired meteorological data, including: 最大值vmax=max(vt),t=t,t-1,t-2,t-3,t-4,...Maximum value v max =max(v t ), t=t, t-1, t-2, t-3, t-4,  … 最小值vmin=min(vt),t=t,t-1,t-2,t-3,t-4,...Minimum value v min =min(v t ), t=t, t-1, t-2, t-3, t-4,  … 均值
Figure QLYQS_4
average
Figure QLYQS_4
标准差
Figure QLYQS_5
standard deviation
Figure QLYQS_5
极值vp=vmax-vmin Extreme value v p =v max -v min 其中vt是指气象数据,t指最近的t个时间点,获得的气象数据中,根据当前时间获得温度和湿度的当前值,对温度进行最大、最小、均值、标准差、极值的计算;湿度进行最大、最小、均值、极值的计算;风速进行最大值计算;风力等级进行最大值计算;降雨量进行最大、最小、极值的计算;Where v t refers to meteorological data, and t refers to the latest t time points. Among the obtained meteorological data, the current value of temperature and humidity is obtained according to the current time, and the maximum, minimum, mean, standard deviation, and extreme value of the temperature are calculated. ; Calculate the maximum, minimum, average, and extreme values of humidity; calculate the maximum value of wind speed; calculate the maximum value of wind power level; calculate the maximum, minimum, and extreme values of rainfall; 对于获得的风向数据进行离散化,因风向数据为字符型,将其转化为数值型,离散化为分类变量,即将北,东东北,东北,东北北,东,东东南,东南,东南南,南,西南南,西南,西西南,西,西西北,西北,西北北16个方位转化为1-16的整数数值;Discretize the obtained wind direction data, because the wind direction data is a character type, convert it into a numerical type, and discretize it into a categorical variable, that is, north, east northeast, northeast, northeast north, east, east southeast, southeast, southeast south, South, Southwest South, Southwest, West Southwest, West, West Northwest, Northwest, Northwest North 16 orientations are converted into integer values from 1 to 16; S4、判断步骤S3获得的数据量是否大于10万,若是,则进入步骤S5,否则进入步骤S6;S4, judging whether the amount of data obtained in step S3 is greater than 100,000, if so, proceed to step S5, otherwise proceed to step S6; S5、基于架空输电线路覆冰中正样本和负样本之间是否相差5倍的数量级来进行类别失衡的判断,若相差5倍,则判定类别失衡,将数据输入bagging模式模型,否则采用机器学习的方法来获得预测结果;所述bagging模式模型为:将输入的数据集进行分割,获得N个子集,针对每个子集构建一个训练模型,即构建N个训练模型,每个子集经过对应的训练模型后,获得N个训练模型的结果,采用bagging算法对N个模型结果进行集成,获得预测结果;S5. Based on whether there is an order of magnitude difference of 5 times between the positive samples and negative samples in the icing of overhead transmission lines, the category imbalance is judged. If the difference is 5 times, the category imbalance is determined, and the data is input into the bagging mode model. Otherwise, the machine learning method is used. method to obtain prediction results; the bagging mode model is: the input data set is divided to obtain N subsets, and a training model is constructed for each subset, that is, N training models are constructed, and each subset undergoes a corresponding training model Finally, the results of N training models are obtained, and the bagging algorithm is used to integrate the results of N models to obtain prediction results; S6、基于架空输电线路覆冰中正样本和负样本之间是否相差5倍的数量级来进行类别失衡的判断,若相差5倍,则判定类别失衡,将数据输入权重模型,否则采用机器学习的方法获得预测结果;所述权重模型为:在目标函数训练时取的权重不同,定义目标函数为S6. Based on whether there is an order of magnitude difference of 5 times between the positive samples and negative samples in the icing of overhead transmission lines, the category imbalance is judged. If the difference is 5 times, the category imbalance is determined, and the data is input into the weight model. Otherwise, the machine learning method is used. Obtain prediction result; Described weight model is: the weight that gets when objective function training is different, defines objective function as
Figure QLYQS_6
Figure QLYQS_6
其中k为类别变量,w为权重,L为loss函数,yk为观测样本中的目标变量,f(xk)为基础模型的输出结果,所述基础模型包括逻辑回归、支持向量机、决策树、boosting集成学习分类方法;Where k is the category variable, w is the weight, L is the loss function, y k is the target variable in the observation sample, f(x k ) is the output result of the basic model, and the basic model includes logistic regression, support vector machine, decision-making Tree, boosting integrated learning classification method; 权重w的计算公式为:The formula for calculating the weight w is:
Figure QLYQS_7
Figure QLYQS_7
其中,n为样本数量,mclass为class对应的样本数量,class取值为目标变量的类别;Among them, n is the number of samples, m class is the number of samples corresponding to the class, and the value of class is the category of the target variable; 根据获得的权重w进行训练后,利用训练好的模型获得预测结果。After training according to the obtained weight w, use the trained model to obtain the prediction result.
CN202010267424.XA 2020-04-08 2020-04-08 Overhead power transmission line icing prediction method based on meteorological and geographic environments Active CN111539842B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010267424.XA CN111539842B (en) 2020-04-08 2020-04-08 Overhead power transmission line icing prediction method based on meteorological and geographic environments

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010267424.XA CN111539842B (en) 2020-04-08 2020-04-08 Overhead power transmission line icing prediction method based on meteorological and geographic environments

Publications (2)

Publication Number Publication Date
CN111539842A CN111539842A (en) 2020-08-14
CN111539842B true CN111539842B (en) 2023-05-23

Family

ID=71974983

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010267424.XA Active CN111539842B (en) 2020-04-08 2020-04-08 Overhead power transmission line icing prediction method based on meteorological and geographic environments

Country Status (1)

Country Link
CN (1) CN111539842B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112183897B (en) * 2020-11-02 2024-06-28 成都卡普数据服务有限责任公司 Long-time prediction method for ice coating thickness of overhead transmission line based on deep learning
US20220171093A1 (en) * 2020-12-01 2022-06-02 International Business Machines Corporation Fractional ice cover predictions with machine learning, satellite, thermodynamics, and in-situ observations
CN112863114B (en) * 2021-01-18 2022-08-12 交通运输部科学研究院 A method and device for determining the time delay of road ice condensation
CN113065695A (en) * 2021-03-22 2021-07-02 广东海启星海洋科技有限公司 Sea ice grade prediction method and device based on multi-network fusion
CN114330061A (en) * 2021-12-28 2022-04-12 广西电网有限责任公司防城港供电局 An Analysis Method for Weak Links of Typical Towers of Transmission Lines

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156793A (en) * 2011-05-05 2011-08-17 中国电力工程顾问集团西南电力设计院 Icing degree classifying system for power transmission line
CN106682771A (en) * 2016-12-14 2017-05-17 云南电网有限责任公司电力科学研究院 Power transmission line coated ice thickness prediction method based on micro meteorological information
CN107092983A (en) * 2017-04-11 2017-08-25 北京国网富达科技发展有限责任公司 Transmission pressure ice covering thickness Forecasting Methodology and device
CN109800905A (en) * 2018-12-19 2019-05-24 国网重庆市电力公司检修分公司 The powerline ice-covering analysis method that mountain environment mima type microrelief microclimate influences
WO2019114160A1 (en) * 2017-12-14 2019-06-20 北京金风科创风电设备有限公司 Method and device for predicting ice formation, and model generation method and device
CN110136023A (en) * 2019-03-28 2019-08-16 清华大学 Icing Risk Prediction of Transmission Lines Based on Adaptive Reinforcement Learning
CN110427857A (en) * 2019-07-26 2019-11-08 国网湖北省电力有限公司检修公司 A kind of transmission line of electricity geological disasters analysis method based on Remote Sensing Data Fusion Algorithm

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103673960B (en) * 2012-08-30 2016-12-21 国际商业机器公司 For the method and apparatus predicting the ice coating state on transmission line of electricity

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156793A (en) * 2011-05-05 2011-08-17 中国电力工程顾问集团西南电力设计院 Icing degree classifying system for power transmission line
CN106682771A (en) * 2016-12-14 2017-05-17 云南电网有限责任公司电力科学研究院 Power transmission line coated ice thickness prediction method based on micro meteorological information
CN107092983A (en) * 2017-04-11 2017-08-25 北京国网富达科技发展有限责任公司 Transmission pressure ice covering thickness Forecasting Methodology and device
WO2019114160A1 (en) * 2017-12-14 2019-06-20 北京金风科创风电设备有限公司 Method and device for predicting ice formation, and model generation method and device
CN109800905A (en) * 2018-12-19 2019-05-24 国网重庆市电力公司检修分公司 The powerline ice-covering analysis method that mountain environment mima type microrelief microclimate influences
CN110136023A (en) * 2019-03-28 2019-08-16 清华大学 Icing Risk Prediction of Transmission Lines Based on Adaptive Reinforcement Learning
CN110427857A (en) * 2019-07-26 2019-11-08 国网湖北省电力有限公司检修公司 A kind of transmission line of electricity geological disasters analysis method based on Remote Sensing Data Fusion Algorithm

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HongYan Li等.Transmission Line Ice Coating Prediction Model Based on EEMD Feature Extraction.IEEE.2019,第7卷第40695-40706页. *
李军等.基于机载在线监测系统的输电线路覆冰预测研究.《现代信息科技》.2019,第3卷(第17期),第34-36+39页. *
郑雪莲.基于CB形态学的输电线路覆冰厚度检测.《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》.2014,(第3期),第C042-498页. *
陈勇等.基于PCA-GA-LSSVM的输电线路覆冰负荷在线预测模型.《电力系统保护与控制》.2019,第47卷(第10期),第110-119页. *

Also Published As

Publication number Publication date
CN111539842A (en) 2020-08-14

Similar Documents

Publication Publication Date Title
CN111539842B (en) Overhead power transmission line icing prediction method based on meteorological and geographic environments
CN110135630B (en) Short-term load demand forecasting method based on random forest regression and multi-step optimization
CN105426970B (en) A kind of meteorological intimidation estimating method based on discrete dynamic Bayesian network
CN104200484B (en) Distributed photovoltaic system ultra-short-term power output prediction method based on cloud cluster characteristic analysis
CN108227041B (en) Horizontal visibility forecast method based on site measured data and model results
CN102706389A (en) System and method for predicting ice accretion on surface of aircraft through rolling forecast within short time
CN112116160B (en) Important power transmission channel disaster monitoring method based on improved cellular automaton of optimized neural network
CN111428942B (en) Line icing thickness prediction method for extracting micro-terrain factors based on variable grid technology
CN112183897B (en) Long-time prediction method for ice coating thickness of overhead transmission line based on deep learning
CN109472396A (en) Mountain fire prediction technique based on depth e-learning
CN105701572A (en) Photovoltaic short-term output prediction method based on improved Gaussian process regression
CN112232591A (en) Intelligent early warning method for ice thickness based on meteorological factors
CN115267945A (en) Thunder and lightning early warning method and system based on graph neural network
CN115629160A (en) Air pollutant concentration prediction method and system based on space-time diagram
CN113689053B (en) Strong convection weather overhead line power failure prediction method based on random forest
CN113821895A (en) Construction method and device of power transmission line icing thickness prediction model and storage medium
CN109299208A (en) An intelligent visualization risk assessment method for transmission towers under typhoon disasters
CN114647985A (en) Training Method of Frontal Transit Prediction Model and Prediction Method of Frontal Transit
CN115935215B (en) Icing warning method and system for transmission lines based on deep learning in extreme weather
CN118536395A (en) A real-time correction method, device and storage medium for mesoscale numerical weather forecast
CN105005661A (en) Transmission tower material actual strength calculation method with practical result
CN104217260A (en) Combined filling system for measured wind speed loss values of multiple neighboring wind motors in wind field
Molinar et al. From data points to ampacity forecasting: Gated recurrent unit networks
CN113159358A (en) Method for predicting urban PM2.5 based on deep space-time residual error network
CN116128307A (en) Transmission line galloping early warning system based on multi-metadata transmission network

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
TR01 Transfer of patent right

Effective date of registration: 20241218

Address after: No. 1907, 19th Floor, Unit 1, Building 1, No. 722 Yizhou Avenue Middle Section, High tech Zone, Chengdu City, Sichuan Province 610094

Patentee after: Sichuan Shenzhou Lvchuang Smart Energy Co.,Ltd.

Country or region after: China

Address before: 610041 Fudi and Fucheng international t1f1907, No. 722, Yizhou Avenue, high tech Zone, Chengdu, Sichuan

Patentee before: CHENGDU SIHAN TECHNOLOGY Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right