CN108321749B - Power transmission line windage yaw forecasting method and device based on numerical meteorological data - Google Patents

Power transmission line windage yaw forecasting method and device based on numerical meteorological data Download PDF

Info

Publication number
CN108321749B
CN108321749B CN201810063123.8A CN201810063123A CN108321749B CN 108321749 B CN108321749 B CN 108321749B CN 201810063123 A CN201810063123 A CN 201810063123A CN 108321749 B CN108321749 B CN 108321749B
Authority
CN
China
Prior art keywords
transmission line
power transmission
meteorological data
wind
wind speed
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
CN201810063123.8A
Other languages
Chinese (zh)
Other versions
CN108321749A (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
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
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, China Electric Power Research Institute Co Ltd CEPRI filed Critical State Grid Corp of China SGCC
Priority to CN201810063123.8A priority Critical patent/CN108321749B/en
Publication of CN108321749A publication Critical patent/CN108321749A/en
Application granted granted Critical
Publication of CN108321749B publication Critical patent/CN108321749B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02G—INSTALLATION OF ELECTRIC CABLES OR LINES, OR OF COMBINED OPTICAL AND ELECTRIC CABLES OR LINES
    • H02G7/00—Overhead installations of electric lines or cables
    • H02G7/14—Arrangements or devices for damping mechanical oscillations of lines, e.g. for reducing production of sound
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Wind Motors (AREA)

Abstract

The invention provides a method and a device for forecasting the windage yaw of a power transmission line based on numerical meteorological data. The technical scheme provided by the invention can forecast the windage yaw state of the target power transmission line in advance, provides technical support for timely and effectively preventing and controlling the windage yaw disaster of the power transmission line, and has wide application prospect in the field of power grid disaster prevention and reduction.

Description

一种基于数值气象数据的输电线路风偏预报方法和装置A method and device for forecasting wind deflection of transmission lines based on numerical meteorological data

技术领域technical field

本发明涉及输电线路防灾减灾技术领域，具体涉及一种基于数值气象数据的输电线路风偏预报方法和装置。The invention relates to the technical field of disaster prevention and mitigation of power transmission lines, in particular to a method and device for forecasting wind deflection of power transmission lines based on numerical meteorological data.

背景技术Background technique

输电线路风偏是指绝缘子串及其悬挂的输电线路在风荷载作用下产生的风偏摇摆现象。风偏可能造成严重的电气、机械故障，对电网安全稳定运行危害较大。为了提高运维阶段输电线路风偏防治工作的针对性和时效性，开展风偏预报十分必要。风偏角是描述输电线路风偏的关键指标，现有技术中的风偏研究主要通过刚体直棒法求解绝缘子串的最大静风偏角，解决了设计阶段静态风偏角及其模型参数的计算问题，却无法实现输电线路动态风偏角的实时预报和预警。准确预报输电线路的风偏角存在两方面困难：1)由于风场的不确定性及微地形因素，线路所在位置的实际风荷载与气象站点观测值有所差异，风偏角预报时应尽量采用实际输电线路处的风场数据；2)由于风场的脉动效应，基于平均风荷载计算得到的静态风偏角不能准确预报输电线路风偏状态。Wind deflection of transmission line refers to the phenomenon of wind deflection caused by the insulator string and its suspended transmission line under the action of wind load. The wind deviation may cause serious electrical and mechanical failures, which are more harmful to the safe and stable operation of the power grid. In order to improve the pertinence and timeliness of the prevention and control of wind deviation of transmission lines in the operation and maintenance stage, it is necessary to carry out wind deviation forecasting. The wind deflection angle is a key index to describe the wind deflection of the transmission line. The wind deflection research in the prior art mainly uses the rigid body straight rod method to solve the maximum static wind deflection angle of the insulator string, which solves the problem of the static wind deflection angle and its model parameters in the design stage. However, the real-time forecast and early warning of the dynamic wind deflection angle of the transmission line cannot be realized. There are two difficulties in accurately forecasting the wind deflection angle of transmission lines: 1) Due to the uncertainty of the wind field and micro-topographic factors, the actual wind load at the location of the line is different from the observation value of the meteorological station, and the wind deflection angle should be predicted as much as possible. Use the wind field data at the actual transmission line; 2) Due to the pulsation effect of the wind field, the static wind deflection angle calculated based on the average wind load cannot accurately predict the wind deflection state of the transmission line.

现有技术基于数值天气预报结果进行输电线路风偏闪络预警，但主要是针对静态风偏角和风偏闪络电压的预测，缺乏输电线路动态风偏角的预报，输电线路风偏状态预报结果准确度低。The prior art is based on numerical weather forecast results for transmission line wind deflection flashover warning, but it is mainly aimed at the prediction of static wind deflection angle and wind deflection flashover voltage. Accuracy is low.

发明内容SUMMARY OF THE INVENTION

为了克服上述现有技术中输电线路风偏状态预报结果准确度的不足，本发明提供一种基于数值气象数据的输电线路风偏预报方法和装置，先对接收的数值气象数据进行解析，并根据解析后的气象数据确定输电线路与风向的夹角和随机风速，然后根据输电线路与风向的夹角和随机风速计算随机风荷载，并根据随机风荷载计算输电线路的水平位移，最后根据输电线路的水平位移确定绝缘子动态风偏角，并根据绝缘子动态风偏角对输电线路风偏状态进行预报，实现基于数值气象数据和绝缘子动态风偏角的输电线路风偏预报，提高了输电线路风偏状态预报结果的准确度。In order to overcome the above-mentioned shortage of the accuracy of the wind deflection state forecasting result of the transmission line in the prior art, the present invention provides a method and device for forecasting the wind deflection of a transmission line based on numerical meteorological data. The analyzed meteorological data determines the angle between the transmission line and the wind direction and the random wind speed, and then calculates the random wind load according to the angle between the transmission line and the wind direction and the random wind speed, and calculates the horizontal displacement of the transmission line according to the random wind load. The horizontal displacement of the insulator determines the dynamic wind deflection angle of the insulator, and forecasts the wind deflection state of the transmission line according to the dynamic wind deflection angle of the insulator. The accuracy of the state forecast results.

为了实现上述发明目的，本发明采取如下技术方案：In order to realize the above-mentioned purpose of the invention, the present invention adopts the following technical solutions:

一方面，本发明提供一种基于数值气象数据的输电线路风偏预报方法，包括：In one aspect, the present invention provides a method for forecasting wind deflection of transmission lines based on numerical meteorological data, comprising:

对接收的数值气象数据进行解析，并根据解析后的气象数据确定输电线路与风向的夹角和随机风速；Analyze the received numerical meteorological data, and determine the angle between the transmission line and the wind direction and the random wind speed according to the analyzed meteorological data;

根据输电线路与风向的夹角和随机风速计算随机风荷载，并根据随机风荷载计算输电线路的水平位移；Calculate the random wind load according to the angle between the transmission line and the wind direction and the random wind speed, and calculate the horizontal displacement of the transmission line according to the random wind load;

根据输电线路的水平位移确定绝缘子动态风偏角，并根据绝缘子动态风偏角对输电线路风偏状态进行预报。The dynamic wind deflection angle of the insulator is determined according to the horizontal displacement of the transmission line, and the wind deflection state of the transmission line is predicted according to the dynamic wind deflection angle of the insulator.

所述对接收的数值气象数据进行解析，包括：The analysis of the received numerical meteorological data includes:

通过地理信息系统对数值气象数据进行坐标矫正，得到坐标矫正后的数值气象数据，并将坐标矫正后的数值气象数据进行保存；Perform coordinate correction on the numerical meteorological data through the geographic information system, obtain the numerical meteorological data after coordinate correction, and save the numerical meteorological data after coordinate correction;

将所有杆塔按照杆塔所在的地理坐标定位到与坐标矫正后的数值气象数据相同精度的网格中，按照杆塔所在的行和列，通过网格索引得到解析后的气象数据，所述解析后的气象数据包括数值气象数据对应的输电线路处的平均风速和风向。All towers are located in a grid with the same precision as the numerical meteorological data after coordinate correction according to the geographic coordinates where the towers are located, and the parsed meteorological data is obtained through the grid index according to the row and column where the towers are located. The meteorological data includes the average wind speed and wind direction at the transmission line corresponding to the numerical meteorological data.

根据解析后的气象数据，按下式确定输电线路与风向的夹角：According to the analyzed meteorological data, the angle between the transmission line and the wind direction is determined as follows:

θ＝arctan|(k2-k1)/(1+k1k2)|θ=arctan|(k 2 -k 1 )/(1+k 1 k 2 )|

其中，θ表示输电线路与风向的夹角，k1表示输电线路斜率，k2表示风向斜率。Among them, θ represents the angle between the transmission line and the wind direction, k 1 represents the slope of the transmission line, and k 2 represents the wind direction slope.

根据解析后的气象数据按下述过程确定随机风速：According to the analyzed meteorological data, the random wind speed is determined according to the following process:

根据标准高度处平均风速，通过

Figure BDA0001555915370000021
确定参考高度处的平均风速，根据
Figure BDA0001555915370000022
并通过谐波叠加法得到随机风速；其中，
Figure BDA0001555915370000023
表示参考高度处的平均风速，v表示标准高度处的平均风速，β表示地貌粗糙度指数。According to the average wind speed at the standard altitude, pass
Figure BDA0001555915370000021
Determine the average wind speed at the reference altitude, according to
Figure BDA0001555915370000022
And the random wind speed is obtained by the harmonic superposition method; among them,
Figure BDA0001555915370000023
Represents the average wind speed at the reference height, v represents the average wind speed at the standard height, and β represents the geomorphic roughness index.

根据输电线路与风向的夹角和随机风速，按下式计算随机风荷载：According to the angle between the transmission line and the wind direction and the random wind speed, the random wind load is calculated as follows:

F(t)＝0.625αμscdv(t)2Lsin2θF(t)=0.625αμ sc dv(t) 2 Lsin 2 θ

其中，F(t)表示t时刻的随机风荷载，v(t)表示t时刻的随机风速，α表示风压不均匀系数，μsc表示输电线路的体型系数，d表示输电线路的外径，L表示输电线路的跨度。Among them, F(t) is the random wind load at time t, v(t) is the random wind speed at time t, α is the wind pressure non-uniformity coefficient, μ sc is the body shape coefficient of the transmission line, d is the outer diameter of the transmission line, L represents the span of the transmission line.

所述根据随机风荷载计算输电线路的水平位移包括：The calculation of the horizontal displacement of the transmission line according to the random wind load includes:

按下式确定输电线路的非线性有限元动力方程：The nonlinear finite element dynamic equation of the transmission line is determined as follows:

Md″(t)+Cd′(t)+Kd(t)＝F(t)Md"(t)+Cd'(t)+Kd(t)=F(t)

其中，M表示绝缘子的质量矩阵，C表示绝缘子的刚度矩阵，K表示绝缘子的阻尼矩阵，d(t)表示输电线路的水平位移，d′(t)表示d(t)的一次导数，d″(t)表示d(t)的二次导数；Among them, M is the mass matrix of the insulator, C is the stiffness matrix of the insulator, K is the damping matrix of the insulator, d(t) is the horizontal displacement of the transmission line, d'(t) is the first derivative of d(t), d" (t) represents the second derivative of d(t);

求解输电线路的非线性有限元动力方程，得到d(t)。Solve the nonlinear finite element dynamic equation of the transmission line to obtain d(t).

根据输电线路的水平位移，按下式确定绝缘子动态风偏角：According to the horizontal displacement of the transmission line, the dynamic wind deflection angle of the insulator is determined as follows:

Figure BDA0001555915370000031
Figure BDA0001555915370000031

其中，

Figure BDA0001555915370000032
表示t时刻绝缘子动态风偏角，l表示绝缘子长度。in,
Figure BDA0001555915370000032
represents the dynamic wind deflection angle of the insulator at time t, and l represents the length of the insulator.

另一方面，本发明提供一种基于数值气象数据的输电线路风偏预报装置，包括：On the other hand, the present invention provides a transmission line wind deflection forecasting device based on numerical meteorological data, comprising:

确定模块，用于对接收的数值气象数据进行解析，并根据解析后的气象数据确定输电线路与风向的夹角和随机风速；The determination module is used to analyze the received numerical meteorological data, and determine the angle between the transmission line and the wind direction and the random wind speed according to the analyzed meteorological data;

计算模块，用于根据输电线路与风向的夹角和随机风速计算随机风荷载，并根据随机风荷载计算输电线路的水平位移；The calculation module is used to calculate the random wind load according to the angle between the transmission line and the wind direction and the random wind speed, and calculate the horizontal displacement of the transmission line according to the random wind load;

预报模块，用于根据输电线路的水平位移确定绝缘子动态风偏角，并根据绝缘子动态风偏角对输电线路风偏状态进行预报。The forecasting module is used to determine the dynamic wind deflection angle of the insulator according to the horizontal displacement of the transmission line, and forecast the wind deflection state of the transmission line according to the dynamic wind deflection angle of the insulator.

所述确定模块包括：The determining module includes:

接收单元，采用文件传输协议并按照预设间隔接收数值气象数据，所述数值气象数据包括标准高度处的平均风速和风向。The receiving unit adopts a file transfer protocol and receives numerical meteorological data at preset intervals, where the numerical meteorological data includes the average wind speed and wind direction at a standard altitude.

所述确定模块包括：The determining module includes:

解析单元，用于对接收的数值气象数据按下述过程进行解析：The parsing unit is used to parse the received numerical meteorological data according to the following process:

通过地理信息系统对数值气象数据进行坐标矫正，得到坐标矫正后的数值气象数据，并将坐标矫正后的数值气象数据进行保存；Perform coordinate correction on the numerical meteorological data through the geographic information system, obtain the numerical meteorological data after coordinate correction, and save the numerical meteorological data after coordinate correction;

将所有杆塔按照杆塔所在的地理坐标定位到与坐标矫正后的数值气象数据相同精度的网格中，按照杆塔所在的行和列，通过网格索引得到解析后的气象数据，所述解析后的气象数据包括数值气象数据对应的输电线路处的平均风速和风向。All towers are located in a grid with the same precision as the numerical meteorological data after coordinate correction according to the geographic coordinates where the towers are located, and the parsed meteorological data is obtained through the grid index according to the row and column where the towers are located. The meteorological data includes the average wind speed and wind direction at the transmission line corresponding to the numerical meteorological data.

所述确定模块包括：The determining module includes:

夹角确定单元，用于根据解析后的气象数据，按下式确定输电线路与风向的夹角：The angle determination unit is used to determine the angle between the transmission line and the wind direction according to the analyzed meteorological data as follows:

θ＝arctan|(k2-k1)/(1+k1k2)|θ=arctan|(k 2 -k 1 )/(1+k 1 k 2 )|

其中，θ表示输电线路与风向的夹角，k1表示输电线路斜率，k2表示风向斜率。Among them, θ represents the angle between the transmission line and the wind direction, k 1 represents the slope of the transmission line, and k 2 represents the wind direction slope.

随机风速确定单元，用于根据解析后的气象数据按下述过程确定随机风速：The random wind speed determination unit is used to determine the random wind speed according to the following process according to the analyzed meteorological data:

根据标准高度处平均风速，通过

Figure BDA0001555915370000033
确定参考高度处的平均风速，根据
Figure BDA0001555915370000034
并通过谐波叠加法得到随机风速；其中，
Figure BDA0001555915370000041
表示参考高度处的平均风速，v表示标准高度处的平均风速，β表示地貌粗糙度指数。According to the average wind speed at the standard altitude, pass
Figure BDA0001555915370000033
Determine the average wind speed at the reference altitude, according to
Figure BDA0001555915370000034
And the random wind speed is obtained by the harmonic superposition method; among them,
Figure BDA0001555915370000041
Represents the average wind speed at the reference height, v represents the average wind speed at the standard height, and β represents the geomorphic roughness index.

所述计算模块包括：The computing module includes:

随机风荷载计算单元，用于根据输电线路与风向的夹角和随机风速，按下式计算随机风荷载：The random wind load calculation unit is used to calculate the random wind load according to the angle between the transmission line and the wind direction and the random wind speed as follows:

F(t)＝0.625αμscdv(t)2Lsin2θF(t)=0.625αμ sc dv(t) 2 Lsin 2 θ

其中，F(t)表示t时刻的随机风荷载，v(t)表示t时刻的随机风速，α表示风压不均匀系数，μsc表示输电线路的体型系数，d表示输电线路的外径，L表示输电线路的跨度。Among them, F(t) is the random wind load at time t, v(t) is the random wind speed at time t, α is the wind pressure non-uniformity coefficient, μ sc is the body shape coefficient of the transmission line, d is the outer diameter of the transmission line, L represents the span of the transmission line.

所述计算模块还包括水平位移计算单元，所述水平位移计算单元包括：The calculation module also includes a horizontal displacement calculation unit, and the horizontal displacement calculation unit includes:

方程确定单元，用于按下式确定输电线路的非线性有限元动力方程：The equation determination unit is used to determine the nonlinear finite element dynamic equation of the transmission line as follows:

Md″(t)+Cd′(t)+Kd(t)＝F(t)Md"(t)+Cd'(t)+Kd(t)=F(t)

其中，M表示绝缘子的质量矩阵，C表示绝缘子的刚度矩阵，K表示绝缘子的阻尼矩阵，d(t)表示输电线路的水平位移，d′(t)表示d(t)的一次导数，d″(t)表示d(t)的二次导数；Among them, M is the mass matrix of the insulator, C is the stiffness matrix of the insulator, K is the damping matrix of the insulator, d(t) is the horizontal displacement of the transmission line, d'(t) is the first derivative of d(t), d" (t) represents the second derivative of d(t);

求解单元，用于求解输电线路的非线性有限元动力方程，得到d(t)。Solver element used to solve the nonlinear finite element dynamic equation of the transmission line to obtain d(t).

所述预报模块包括：The forecast module includes:

动态风偏角确定单元，用于根据输电线路的水平位移，按下式确定绝缘子动态风偏角：The dynamic wind deflection angle determination unit is used to determine the dynamic wind deflection angle of the insulator according to the horizontal displacement of the transmission line as follows:

Figure BDA0001555915370000042
Figure BDA0001555915370000042

其中，

Figure BDA0001555915370000043
表示t时刻绝缘子动态风偏角，l表示绝缘子长度；in,
Figure BDA0001555915370000043
represents the dynamic wind deflection angle of the insulator at time t, and l represents the length of the insulator;

预报单元，根据绝缘子动态风偏角对输电线路风偏状态进行预报。The forecasting unit forecasts the wind deflection state of the transmission line according to the dynamic wind deflection angle of the insulator.

与最接近的现有技术相比，本发明提供的技术方案具有以下有益效果：Compared with the closest prior art, the technical solution provided by the present invention has the following beneficial effects:

本发明提供的基于数值气象数据的输电线路风偏预报方法中，先对接收的数值气象数据进行解析，并根据解析后的气象数据确定输电线路与风向的夹角和随机风速，然后根据输电线路与风向的夹角和随机风速计算随机风荷载，并根据随机风荷载计算输电线路的水平位移，最后根据输电线路的水平位移确定绝缘子动态风偏角，并根据绝缘子动态风偏角对输电线路风偏状态进行预报，实现基于数值气象数据和绝缘子动态风偏角的输电线路风偏预报，提高了输电线路风偏状态预报结果的准确度；In the method for forecasting wind deflection of transmission lines based on numerical meteorological data provided by the present invention, the received numerical meteorological data is first analyzed, and the angle between the transmission line and the wind direction and the random wind speed are determined according to the analyzed meteorological data, and then the angle between the transmission line and the wind direction and the random wind speed are determined according to the analyzed meteorological data. Calculate the random wind load based on the angle between the wind direction and the random wind speed, and calculate the horizontal displacement of the transmission line according to the random wind load. It can forecast the wind deflection of transmission lines based on numerical meteorological data and the dynamic wind deflection angle of insulators, and improve the accuracy of the forecast results of wind deflection of transmission lines;

本发明实施例还提供的基于数值气象数据的输电线路风偏预报装置包括确定模块、计算模块和预报模块，确定模块，用于对接收的数值气象数据进行解析，并根据解析后的气象数据确定输电线路与风向的夹角和随机风速；计算模块，用于根据输电线路与风向的夹角和随机风速计算随机风荷载，并根据随机风荷载计算输电线路的水平位移；预报模块，用于根据输电线路的水平位移确定绝缘子动态风偏角，并根据绝缘子动态风偏角对输电线路风偏状态进行预报，实现基于数值气象数据和绝缘子动态风偏角的输电线路风偏预报，提高了输电线路风偏状态预报结果的准确度；The apparatus for forecasting wind deviation of transmission lines based on numerical meteorological data further provided by the embodiment of the present invention includes a determination module, a calculation module and a forecasting module. The determination module is used to analyze the received numerical meteorological data, and determine according to the analyzed meteorological data. The angle between the transmission line and the wind direction and random wind speed; the calculation module is used to calculate the random wind load according to the angle between the transmission line and the wind direction and the random wind speed, and the horizontal displacement of the transmission line is calculated according to the random wind load; the forecast module is used to calculate the horizontal displacement of the transmission line according to the random wind load. The horizontal displacement of the transmission line determines the dynamic wind deflection angle of the insulator, and forecasts the wind deflection state of the transmission line according to the dynamic wind deflection angle of the insulator, realizes the wind deflection forecast of the transmission line based on the numerical meteorological data and the dynamic wind deflection angle of the insulator, and improves the transmission line wind deflection. The accuracy of the forecast results of the wind deviation state;

本发明提供的技术方案中，通过地理信息系统对采用文件传输协议接收的数值气象数据进行坐标矫正和索引定位，准确获得解析后的数值气象数据；In the technical solution provided by the present invention, coordinate correction and index positioning are performed on the numerical meteorological data received by using the file transfer protocol through the geographic information system, and the parsed numerical meteorological data is accurately obtained;

本发明提供的技术方案采用目标输电线路处的平均风速和风向进行输电线路风偏预报，提高了随机风荷载的精度，避免了采用一般区域性气象预报数据引起的误差；The technical scheme provided by the invention adopts the average wind speed and wind direction at the target transmission line to forecast the wind deviation of the transmission line, which improves the accuracy of random wind loads and avoids errors caused by using general regional weather forecast data;

本发明提供的技术方案对目标输电线路风偏状态的提前预报，为及时、有效地开展输电线路风偏灾害防治提供了技术支撑，在电网防灾减灾领域具有广泛的应用前景。The technical solution provided by the invention provides technical support for the advance forecast of the wind deviation state of the target transmission line, provides technical support for timely and effective prevention and control of the wind deviation disaster of the transmission line, and has wide application prospects in the field of power grid disaster prevention and mitigation.

附图说明Description of drawings

图1是本发明实施例中基于数值气象数据的输电线路风偏预报方法流程图；Fig. 1 is a flow chart of a method for forecasting wind deflection of transmission lines based on numerical meteorological data in an embodiment of the present invention;

图2是本发明实施例中输电线路处的随机风速示意图；2 is a schematic diagram of random wind speed at the transmission line in the embodiment of the present invention;

图3是本发明实施例中输电线路的水平位移示意图；3 is a schematic diagram of the horizontal displacement of the transmission line in the embodiment of the present invention;

图4是本发明实施例中绝缘子动态风偏角示意图。4 is a schematic diagram of the dynamic wind deflection angle of an insulator in an embodiment of the present invention.

具体实施方式Detailed ways

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

本发明实施例提供一种基于数值气象数据的输电线路风偏预报方法，具体流程图如图1所示，具体过程如下：An embodiment of the present invention provides a method for forecasting wind deflection of a transmission line based on numerical meteorological data. The specific flowchart is shown in FIG. 1 , and the specific process is as follows:

S101：对接收的数值气象数据进行解析，并根据解析后的气象数据确定输电线路与风向的夹角和随机风速；S101: Analyze the received numerical meteorological data, and determine the angle between the transmission line and the wind direction and the random wind speed according to the analyzed meteorological data;

S102：根据S101确定的输电线路与风向的夹角和随机风速计算随机风荷载，并根据随机风荷载计算输电线路的水平位移；S102: Calculate the random wind load according to the angle between the transmission line and the wind direction and the random wind speed determined in S101, and calculate the horizontal displacement of the transmission line according to the random wind load;

S103：根据S102计算的输电线路的水平位移确定绝缘子动态风偏角，并根据绝缘子动态风偏角对输电线路风偏状态进行预报。S103: Determine the dynamic wind deflection angle of the insulator according to the horizontal displacement of the transmission line calculated in S102, and predict the wind deflection state of the transmission line according to the dynamic wind deflection angle of the insulator.

上述S101中，数值气象数据采用文件传输协议并按照预设间隔接收，数值气象数据包括标准高度处的平均风速和风向。In the above S101, the numerical meteorological data is received at preset intervals using a file transmission protocol, and the numerical meteorological data includes the average wind speed and wind direction at a standard altitude.

上述S101中，对接收的数值气象数据进行解析的具体过程如下：In the above S101, the specific process of parsing the received numerical meteorological data is as follows:

1)通过地理信息系统对数值气象数据进行坐标矫正，得到坐标矫正后的数值气象数据，并将坐标矫正后的数值气象数据进行保存；1) Perform coordinate correction on the numerical meteorological data through the geographic information system, obtain the numerical meteorological data after the coordinate correction, and save the numerical meteorological data after the coordinate correction;

2)将所有杆塔按照杆塔所在的地理坐标定位到与坐标矫正后的数值气象数据相同精度的网格中，按照杆塔所在的行和列，通过网格索引得到解析后的气象数据，所述解析后的气象数据包括数值气象数据对应的输电线路处的平均风速和风向。2) Position all towers in a grid with the same precision as the numerical meteorological data after coordinate correction according to the geographic coordinates where the towers are located, and obtain the parsed meteorological data through the grid index according to the row and column where the towers are located. The latter meteorological data includes the average wind speed and wind direction at the transmission line corresponding to the numerical meteorological data.

上述S101根据解析后的气象数据确定输电线路与风向的夹角和随机风速的具体过程如下：The specific process of determining the angle between the transmission line and the wind direction and the random wind speed according to the analyzed meteorological data in the above S101 is as follows:

1)根据解析后的气象数据，按下式确定输电线路与风向的夹角：1) According to the analyzed meteorological data, determine the angle between the transmission line and the wind direction as follows:

θ＝arctan|(k2-k1)/(1+k1k2)|θ=arctan|(k 2 -k 1 )/(1+k 1 k 2 )|

其中，θ表示输电线路与风向的夹角，k1表示输电线路斜率，k2表示风向斜率；Among them, θ represents the angle between the transmission line and the wind direction, k 1 represents the slope of the transmission line, and k 2 represents the slope of the wind direction;

2)根据解析后的气象数据按下述过程确定随机风速：2) According to the analyzed meteorological data, determine the random wind speed according to the following process:

根据标准高度处平均风速，通过

Figure BDA0001555915370000061
确定参考高度处的平均风速，根据
Figure BDA0001555915370000062
并通过谐波叠加法得到随机风速(如图2所示)；其中，
Figure BDA0001555915370000063
表示参考高度处的平均风速，v表示标准高度处的平均风速，β表示地貌粗糙度指数。According to the average wind speed at the standard altitude, pass
Figure BDA0001555915370000061
Determine the average wind speed at the reference altitude, according to
Figure BDA0001555915370000062
And the random wind speed is obtained by the harmonic superposition method (as shown in Figure 2); among them,
Figure BDA0001555915370000063
Represents the average wind speed at the reference height, v represents the average wind speed at the standard height, and β represents the geomorphic roughness index.

上述S102中，根据输电线路与风向的夹角和随机风速，按下式计算随机风荷载：In the above S102, according to the angle between the transmission line and the wind direction and the random wind speed, the random wind load is calculated as follows:

F(t)＝0.625αμscdv(t)2Lsin2θF(t)=0.625αμ sc dv(t) 2 Lsin 2 θ

其中，F(t)表示t时刻的随机风荷载，v(t)表示t时刻的随机风速，α表示风压不均匀系数，μsc表示输电线路的体型系数，d表示输电线路的外径，L表示输电线路的跨度。Among them, F(t) is the random wind load at time t, v(t) is the random wind speed at time t, α is the wind pressure non-uniformity coefficient, μ sc is the body shape coefficient of the transmission line, d is the outer diameter of the transmission line, L represents the span of the transmission line.

上述S102中，根据随机风荷载计算输电线路的水平位移包括：In the above S102, calculating the horizontal displacement of the transmission line according to the random wind load includes:

1)按下式确定输电线路的非线性有限元动力方程：1) Determine the nonlinear finite element dynamic equation of the transmission line as follows:

Md″(t)+Cd′(t)+Kd(t)＝F(t)Md"(t)+Cd'(t)+Kd(t)=F(t)

其中，M表示绝缘子的质量矩阵，C表示绝缘子的刚度矩阵，K表示绝缘子的阻尼矩阵，d(t)表示输电线路的水平位移，d′(t)表示d(t)的一次导数，d″(t)表示d(t)的二次导数；Among them, M is the mass matrix of the insulator, C is the stiffness matrix of the insulator, K is the damping matrix of the insulator, d(t) is the horizontal displacement of the transmission line, d'(t) is the first derivative of d(t), d" (t) represents the second derivative of d(t);

2)求解输电线路的非线性有限元动力方程，得到d(t)，如图3所示。2) Solve the nonlinear finite element dynamic equation of the transmission line to obtain d(t), as shown in Figure 3.

上述S103中，根据输电线路的水平位移，按下式确定绝缘子动态风偏角：In the above S103, according to the horizontal displacement of the transmission line, the dynamic wind deflection angle of the insulator is determined as follows:

Figure BDA0001555915370000071
Figure BDA0001555915370000071

其中，

Figure BDA0001555915370000072
表示t时刻绝缘子动态风偏角(如图4所示)，l表示绝缘子长度。in,
Figure BDA0001555915370000072
represents the dynamic wind deflection angle of the insulator at time t (as shown in Figure 4), and l represents the length of the insulator.

基于同一发明构思，本发明实施例还提供一种基于数值气象数据的输电线路风偏预报装置，包括确定模块、计算模块和预报模块，下面分别对上述几个模块的功能进行详细介绍：Based on the same inventive concept, an embodiment of the present invention also provides a wind deflection forecasting device for transmission lines based on numerical meteorological data, including a determination module, a calculation module and a forecasting module. The functions of the above modules are described in detail below:

其中的确定模块，用于对接收的数值气象数据进行解析，并根据解析后的气象数据确定输电线路与风向的夹角和随机风速；The determination module is used to analyze the received numerical meteorological data, and determine the angle between the transmission line and the wind direction and the random wind speed according to the analyzed meteorological data;

其中的计算模块，用于根据输电线路与风向的夹角和随机风速计算随机风荷载，并根据随机风荷载计算输电线路的水平位移；The calculation module is used to calculate the random wind load according to the angle between the transmission line and the wind direction and the random wind speed, and calculate the horizontal displacement of the transmission line according to the random wind load;

其中的预报模块，用于根据输电线路的水平位移确定绝缘子动态风偏角，并根据绝缘子动态风偏角对输电线路风偏状态进行预报。The forecast module is used to determine the dynamic wind deflection angle of the insulator according to the horizontal displacement of the transmission line, and forecast the wind deflection state of the transmission line according to the dynamic wind deflection angle of the insulator.

上述的确定模块包括：The above determination module includes:

1)接收单元，采用文件传输协议并按照预设间隔接收数值气象数据，所述数值气象数据包括标准高度处的平均风速和风向。1) A receiving unit, which adopts a file transfer protocol and receives numerical meteorological data at preset intervals, where the numerical meteorological data includes the average wind speed and wind direction at a standard altitude.

2)解析单元，用于对接收的数值气象数据按下述过程进行解析：2) Analysis unit, used to analyze the received numerical meteorological data according to the following process:

2-1)通过地理信息系统对数值气象数据进行坐标矫正，得到坐标矫正后的数值气象数据，并将坐标矫正后的数值气象数据进行保存；2-1) Perform coordinate correction on the numerical meteorological data through the geographic information system, obtain the numerical meteorological data after the coordinate correction, and save the numerical meteorological data after the coordinate correction;

2-2)将所有杆塔按照杆塔所在的地理坐标定位到与坐标矫正后的数值气象数据相同精度的网格中，按照杆塔所在的行和列，通过网格索引得到解析后的气象数据，所述解析后的气象数据包括数值气象数据对应的输电线路处的平均风速和风向。2-2) Position all towers in a grid with the same precision as the numerical meteorological data after coordinate correction according to the geographic coordinates where the towers are located, and obtain the parsed meteorological data through the grid index according to the row and column where the towers are located. The parsed meteorological data includes the average wind speed and wind direction at the transmission line corresponding to the numerical meteorological data.

3)夹角确定单元，用于根据解析后的气象数据，按下式确定输电线路与风向的夹角：3) The included angle determination unit is used to determine the included angle between the transmission line and the wind direction according to the analyzed meteorological data as follows:

θ＝arctan|(k2-k1)/(1+k1k2)|θ=arctan|(k 2 -k 1 )/(1+k 1 k 2 )|

其中，θ表示输电线路与风向的夹角，k1表示输电线路斜率，k2表示风向斜率。Among them, θ represents the angle between the transmission line and the wind direction, k 1 represents the slope of the transmission line, and k 2 represents the wind direction slope.

4)机风速确定单元，用于根据标准高度处平均风速，通过

Figure BDA0001555915370000073
确定参考高度处的平均风速，根据
Figure BDA0001555915370000074
并通过谐波叠加法得到随机风速(如图2所示)；其中，
Figure BDA0001555915370000075
表示参考高度处的平均风速，v表示标准高度处的平均风速，β表示地貌粗糙度指数。4) The wind speed determination unit is used to determine the average wind speed at the standard height by
Figure BDA0001555915370000073
Determine the average wind speed at the reference altitude, according to
Figure BDA0001555915370000074
And the random wind speed is obtained by the harmonic superposition method (as shown in Figure 2); among them,
Figure BDA0001555915370000075
Represents the average wind speed at the reference height, v represents the average wind speed at the standard height, and β represents the geomorphic roughness index.

上述的计算模块包括：The above calculation modules include:

1)随机风荷载计算单元，用于根据输电线路与风向的夹角和随机风速，按下式计算随机风荷载：1) The random wind load calculation unit is used to calculate the random wind load according to the angle between the transmission line and the wind direction and the random wind speed as follows:

F(t)＝0.625αμscdv(t)2Lsin2θF(t)=0.625αμ sc dv(t) 2 Lsin 2 θ

其中，F(t)表示t时刻的随机风荷载，v(t)表示t时刻的随机风速，α表示风压不均匀系数，μsc表示输电线路的体型系数，d表示输电线路的外径，L表示输电线路的跨度。Among them, F(t) is the random wind load at time t, v(t) is the random wind speed at time t, α is the wind pressure non-uniformity coefficient, μ sc is the body shape coefficient of the transmission line, d is the outer diameter of the transmission line, L represents the span of the transmission line.

2)水平位移计算单元，水平位移计算单元包括：2) Horizontal displacement calculation unit, the horizontal displacement calculation unit includes:

2-1)方程确定单元，用于按下式确定输电线路的非线性有限元动力方程：2-1) The equation determination unit is used to determine the nonlinear finite element dynamic equation of the transmission line as follows:

Md″(t)+Cd′(t)+Kd(t)＝F(t)Md"(t)+Cd'(t)+Kd(t)=F(t)

其中，M表示绝缘子的质量矩阵，C表示绝缘子的刚度矩阵，K表示绝缘子的阻尼矩阵，d(t)表示输电线路的水平位移，d′(t)表示d(t)的一次导数，d″(t)表示d(t)的二次导数；Among them, M is the mass matrix of the insulator, C is the stiffness matrix of the insulator, K is the damping matrix of the insulator, d(t) is the horizontal displacement of the transmission line, d'(t) is the first derivative of d(t), d" (t) represents the second derivative of d(t);

2-1)求解单元，用于求解输电线路的非线性有限元动力方程，得到d(t)。2-1) The solving unit is used to solve the nonlinear finite element dynamic equation of the transmission line to obtain d(t).

上述的预报模块包括：The aforementioned forecast modules include:

1)动态风偏角确定单元，用于根据输电线路的水平位移，按下式确定绝缘子动态风偏角：1) The dynamic wind deflection angle determination unit is used to determine the dynamic wind deflection angle of the insulator according to the horizontal displacement of the transmission line as follows:

Figure BDA0001555915370000081
Figure BDA0001555915370000081

其中，

Figure BDA0001555915370000082
表示t时刻绝缘子动态风偏角，l表示绝缘子长度；in,
Figure BDA0001555915370000082
represents the dynamic wind deflection angle of the insulator at time t, and l represents the length of the insulator;

2)预报单元，根据绝缘子动态风偏角对输电线路风偏状态进行预报。2) Forecasting unit, which forecasts the wind deflection state of the transmission line according to the dynamic wind deflection angle of the insulator.

为了描述的方便，以上所述装置的各部分以功能分为各种模块或单元分别描述。当然，在实施本申请时可以把各模块或单元的功能在同一个或多个软件或硬件中实现。For the convenience of description, each part of the device described above is divided into various modules or units by function and described respectively. Of course, when implementing the present application, the functions of each module or unit may be implemented in one or more software or hardware.

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。As will be appreciated by those skilled in the art, the embodiments of the present application may be provided as a method, a system, or a computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the present application. It will be understood that each flow and/or block in the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to the processor of a general purpose computer, special purpose computer, embedded processor or other programmable data processing device to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing device produce Means for implementing the functions specified in a flow or flow of a flowchart and/or a block or blocks of a block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory result in an article of manufacture comprising instruction means, the instructions The apparatus implements the functions specified in the flow or flow of the flowcharts and/or the block or blocks of the block diagrams.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded on a computer or other programmable data processing device to cause a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process such that The instructions provide steps for implementing the functions specified in the flow or blocks of the flowcharts and/or the block or blocks of the block diagrams.

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Those of ordinary skill in the art can still modify or equivalently replace the specific embodiments of the present invention with reference to the above embodiments. Any modifications or equivalent substitutions that depart from the spirit and scope of the present invention are all within the protection scope of the claims of the present invention for which the application is pending.

Claims (8)

1. A power transmission line windage yaw forecasting method based on numerical meteorological data is characterized by comprising the following steps:
analyzing the received numerical meteorological data, and determining an included angle between the power transmission line and the wind direction and a random wind speed according to the analyzed meteorological data;
calculating a random wind load according to an included angle between the power transmission line and the wind direction and a random wind speed, and calculating the horizontal displacement of the power transmission line according to the random wind load;
determining a dynamic windage yaw angle of the insulator according to the horizontal displacement of the power transmission line, and forecasting the windage yaw state of the power transmission line according to the dynamic windage yaw angle of the insulator;
the analyzing the received numerical meteorological data comprises the following steps:
coordinate correction is carried out on the numerical meteorological data through a geographic information system to obtain the numerical meteorological data after coordinate correction, and the numerical meteorological data after coordinate correction is stored;
positioning all towers into grids with the same precision as the numerical meteorological data after coordinate correction according to the geographic coordinates of the towers, and obtaining analyzed meteorological data through grid indexing according to the rows and columns of the towers, wherein the analyzed meteorological data comprise the average wind speed and the wind direction of the power transmission line corresponding to the numerical meteorological data;
according to the analyzed meteorological data, determining an included angle between the power transmission line and the wind direction according to the following formula:
θ＝arctan|(k2-k1)/(1+k1k2)|
wherein theta represents the included angle between the power transmission line and the wind direction, and k1Representing the slope, k, of the transmission line2Represents the wind direction slope; determining the random wind speed according to the analyzed meteorological data by the following process:
according to the average wind speed at standard altitude, by
Figure FDA0003363203900000011
Determining an average wind speed at a reference altitude, based on
Figure FDA0003363203900000012
Obtaining random wind speed by a harmonic superposition method; wherein,
Figure FDA0003363203900000013
denotes an average wind speed at a reference altitude, v denotes an average wind speed at a standard altitude, β denotes a terrain roughness index, and z denotes a reference altitude.
2. The method for forecasting the windage yaw of the power transmission line based on the numerical meteorological data as recited in claim 1, wherein the random wind load is calculated according to the included angle between the power transmission line and the wind direction and the random wind speed as follows:
F(t)＝0.625αμscdv(t)2Lsin2θ
wherein F (t) represents the random wind load at the time t, v (t) represents the random wind speed at the time t, alpha represents the wind pressure uneven coefficient, muscAnd the size coefficient of the power transmission line is represented, d represents the outer diameter of the power transmission line, and L represents the span of the power transmission line.
3. The method for forecasting the windage yaw of the power transmission line based on the numerical meteorological data as claimed in claim 2, wherein the calculating the horizontal displacement of the power transmission line according to the random wind load comprises:
determining a nonlinear finite element power equation of the power transmission line according to the following formula:
Md″(t)+Cd′(t)+Kd(t)＝F(t)
wherein M represents a mass matrix of the insulator, C represents a rigidity matrix of the insulator, K represents a damping matrix of the insulator, d (t) represents the horizontal displacement of the transmission line, d '(t) represents a first derivative of d (t), and d' (t) represents a second derivative of d (t);
solving the nonlinear finite element power equation of the transmission line to obtain d (t).
4. The method for forecasting the windage yaw of the power transmission line based on the numerical meteorological data as recited in claim 3, wherein the dynamic windage yaw angle of the insulator is determined according to the following formula according to the horizontal displacement of the power transmission line:
Figure FDA0003363203900000021
wherein,
Figure FDA0003363203900000022
and the dynamic windage yaw angle of the insulator at the time t is shown, and l represents the length of the insulator.
5. The utility model provides a transmission line windage yaw forecasting device based on numerical value meteorological data which characterized in that includes:
the determining module is used for analyzing the received numerical meteorological data and determining an included angle between the power transmission line and the wind direction and a random wind speed according to the analyzed meteorological data;
the calculation module is used for calculating random wind load according to the included angle between the power transmission line and the wind direction and the random wind speed and calculating the horizontal displacement of the power transmission line according to the random wind load;
the forecasting module is used for determining a dynamic windage yaw angle of the insulator according to the horizontal displacement of the power transmission line and forecasting the windage yaw state of the power transmission line according to the dynamic windage yaw angle of the insulator;
the determining module comprises:
the analysis unit is used for analyzing the received numerical meteorological data according to the following processes:
coordinate correction is carried out on the numerical meteorological data through a geographic information system to obtain the numerical meteorological data after coordinate correction, and the numerical meteorological data after coordinate correction is stored;
positioning all towers into grids with the same precision as the numerical meteorological data after coordinate correction according to the geographic coordinates of the towers, and obtaining analyzed meteorological data through grid indexing according to the rows and columns of the towers, wherein the analyzed meteorological data comprise the average wind speed and the wind direction of the power transmission line corresponding to the numerical meteorological data;
the determining module further comprises:
and the included angle determining unit is used for determining the included angle between the power transmission line and the wind direction according to the analyzed meteorological data and the following formula:
θ＝arctan|(k2-k1)/(1+k1k2)|
wherein theta represents the included angle between the power transmission line and the wind direction, and k1Representing the slope, k, of the transmission line2Represents the wind direction slope; the random wind speed determining unit is used for determining the random wind speed according to the analyzed meteorological data according to the following processes:
according to the average wind speed at standard altitude, by
Figure FDA0003363203900000032
Determining an average wind speed at a reference altitude, based on
Figure FDA0003363203900000033
Obtaining random wind speed by a harmonic superposition method; wherein,
Figure FDA0003363203900000034
denotes an average wind speed at a reference altitude, v denotes an average wind speed at a standard altitude, β denotes a terrain roughness index, and z denotes a reference altitude.
6. The device according to claim 5, wherein the computing module comprises:
the random wind load calculation unit is used for calculating the random wind load according to the included angle between the power transmission line and the wind direction and the random wind speed and according to the following formula:
F(t)＝0.625αμscdv(t)2Lsin2θ
wherein F (t) represents the random wind load at the time t, v (t) represents the random wind speed at the time t, alpha represents the wind pressure uneven coefficient, muscAnd the size coefficient of the power transmission line is represented, d represents the outer diameter of the power transmission line, and L represents the span of the power transmission line.
7. The device according to claim 6, wherein the calculation module further comprises a horizontal displacement calculation unit, and the horizontal displacement calculation unit comprises: the equation determining unit is used for determining a nonlinear finite element power equation of the power transmission line according to the following formula:
Md″(t)+Cd′(t)+Kd(t)＝F(t)
wherein M represents a mass matrix of the insulator, C represents a rigidity matrix of the insulator, K represents a damping matrix of the insulator, d (t) represents the horizontal displacement of the transmission line, d '(t) represents a first derivative of d (t), and d' (t) represents a second derivative of d (t);
and the solving unit is used for solving the nonlinear finite element power equation of the power transmission line to obtain d (t).
8. The device according to claim 7, wherein the forecasting module comprises:
the dynamic wind drift angle determining unit is used for determining the dynamic wind drift angle of the insulator according to the horizontal displacement of the power transmission line and the following formula:
Figure FDA0003363203900000031
wherein,
Figure FDA0003363203900000035
the dynamic wind deflection angle of the insulator at the moment t is shown, and l represents the length of the insulator;
and the forecasting unit is used for forecasting the windage yaw state of the power transmission line according to the dynamic windage yaw angle of the insulator.
CN201810063123.8A 2018-01-23 2018-01-23 Power transmission line windage yaw forecasting method and device based on numerical meteorological data Active CN108321749B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810063123.8A CN108321749B (en) 2018-01-23 2018-01-23 Power transmission line windage yaw forecasting method and device based on numerical meteorological data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810063123.8A CN108321749B (en) 2018-01-23 2018-01-23 Power transmission line windage yaw forecasting method and device based on numerical meteorological data

Publications (2)

Publication Number Publication Date
CN108321749A CN108321749A (en) 2018-07-24
CN108321749B true CN108321749B (en) 2022-03-18

Family

ID=62887376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810063123.8A Active CN108321749B (en) 2018-01-23 2018-01-23 Power transmission line windage yaw forecasting method and device based on numerical meteorological data

Country Status (1)

Country Link
CN (1) CN108321749B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271751B (en) * 2018-11-16 2022-08-26 重庆科技学院 Method for determining maximum dynamic wind deflection angle of suspension insulator string
CN109635444A (en) * 2018-12-13 2019-04-16 中国电力科学研究院有限公司 A kind of electric transmission line isolator dynamic windage yaw analysis method and device
CN113821897B (en) * 2021-11-25 2022-04-12 国网江西省电力有限公司电力科学研究院 Complex micro-terrain recognition and power transmission line windage yaw calculation method based on elevation image

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20100046908A (en) * 2008-10-28 2010-05-07 한국전력공사 Measuring system for power line tower on the sea
CN103683086A (en) * 2013-12-17 2014-03-26 国家电网公司 Method for electric transmission line windage yaw transient analysis under action of moving thunderstorm wind
CN105046581A (en) * 2015-07-09 2015-11-11 广西电网有限责任公司电力科学研究院 Multi-parameter wind load-based power transmission line strong wind tripping risk assessment method
CN105184492A (en) * 2015-09-10 2015-12-23 国网福建省电力有限公司 Electric power typhoon disaster resistance simulation analysis early warning system based on three-dimensional digital Earth
CN106096105A (en) * 2016-06-02 2016-11-09 浙江大学 Power transmission circuit caused by windage transient response computational methods
CN106096110A (en) * 2016-06-03 2016-11-09 国网电力科学研究院武汉南瑞有限责任公司 Flashover of power transmission circuit caused by windage yaw method for early warning based on numerical weather forecast
CN106500648A (en) * 2016-12-08 2017-03-15 北京国网富达科技发展有限责任公司 Power transmission circuit caused by windage monitoring method and device based on dynamic wind

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20100046908A (en) * 2008-10-28 2010-05-07 한국전력공사 Measuring system for power line tower on the sea
CN103683086A (en) * 2013-12-17 2014-03-26 国家电网公司 Method for electric transmission line windage yaw transient analysis under action of moving thunderstorm wind
CN105046581A (en) * 2015-07-09 2015-11-11 广西电网有限责任公司电力科学研究院 Multi-parameter wind load-based power transmission line strong wind tripping risk assessment method
CN105184492A (en) * 2015-09-10 2015-12-23 国网福建省电力有限公司 Electric power typhoon disaster resistance simulation analysis early warning system based on three-dimensional digital Earth
CN106096105A (en) * 2016-06-02 2016-11-09 浙江大学 Power transmission circuit caused by windage transient response computational methods
CN106096110A (en) * 2016-06-03 2016-11-09 国网电力科学研究院武汉南瑞有限责任公司 Flashover of power transmission circuit caused by windage yaw method for early warning based on numerical weather forecast
CN106500648A (en) * 2016-12-08 2017-03-15 北京国网富达科技发展有限责任公司 Power transmission circuit caused by windage monitoring method and device based on dynamic wind

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于数值天气预报的输电线路风偏闪络预警方法;王海涛，谷山强等;《电力系统保护与控制》;20170616;第121-127页 *
脉动风载下的输电线路风偏计算研究;杨悦;《中国优秀硕士学位论文全文数据库（工程科技Ⅱ辑）》;20170315;第C04-2702页第2-3章 *

Also Published As

Publication number Publication date
CN108321749A (en) 2018-07-24

Similar Documents

Publication Publication Date Title
CN102663251B (en) Physical prediction method for wind power station power based on computational fluid mechanics model
CN108321749B (en) Power transmission line windage yaw forecasting method and device based on numerical meteorological data
CN110908014B (en) A method and system for galloping refined correction forecasting
CN109059937B (en) On-satellite autonomous satellite orbit forecasting method
CN105260953A (en) Evaluation simulation system and evaluation method for theoretical output of integrated wind power plant
CN104794315B (en) Based on the Seam Floor Failure band depth determination method for improving fracture mechanics model
CN112132364A (en) Photovoltaic power station power prediction method, medium and equipment influenced by cloud layer
CN113553747A (en) Mechanical characteristic analysis and state evaluation method for icing power transmission tower
CN103942361A (en) Automatic calculation method and module of support hanger in three-dimensional process plant design
CN113313291B (en) Power grid wind disaster ensemble forecasting and verifying method based on initial disturbance
CN108808671A (en) A kind of short-term wind speed DATA PROCESSING IN ENSEMBLE PREDICTION SYSTEM method of wind power plant
AU2023202496A1 (en) An on-line monitoring method and system for fugitive carbon emission in cement industrial plant area
CN105241464B (en) A kind of crane rail slash figure fast matching method
CN108038333B (en) A Prediction Method of Flexible Machining Surface Error for Plane Milling with Large Disc Milling Cutter
CN110674571B (en) A method and system for downscaling calculation of wind speed at pass of transmission line
CN102495230B (en) Correction method used for prediction of wind speed
CN107330189B (en) Method and system for layered and refined prediction of adjacent grid temperature
CN103091551B (en) The wind power plant generated energy computational methods of unimodal Non-Gaussian Distribution
CN109033000A (en) A kind of the photovoltaic cloud computing cluster control system and method for spring high-efficient
CN112580870A (en) Disaster-affected pole tower determination method and device
CN118536799A (en) Short-coming early warning method and device for running risk of overhead transmission line under strong convection and strong wind
CN109115102B (en) Automatic layering and layer thickness data correction method in roadbed compaction process
CN110645939A (en) System and method for dynamically predicting relative vertical distance between two wires
CN117371581A (en) New energy generated power prediction method, device and storage medium
CN107423509B (en) A method and device for constructing aeroelastic model of transmission line tower system

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
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Patentee after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd.

Country or region after: China

Patentee after: STATE GRID CORPORATION OF CHINA

Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Patentee before: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd.

Country or region before: China

Patentee before: State Grid Corporation of China