CN117114428B - 一种电力设备气象灾害分析预警方法 - Google Patents
一种电力设备气象灾害分析预警方法 Download PDFInfo
- Publication number
- CN117114428B CN117114428B CN202311392865.2A CN202311392865A CN117114428B CN 117114428 B CN117114428 B CN 117114428B CN 202311392865 A CN202311392865 A CN 202311392865A CN 117114428 B CN117114428 B CN 117114428B
- Authority
- CN
- China
- Prior art keywords
- model
- power equipment
- steps
- constructing
- weather
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 89
- 238000004458 analytical method Methods 0.000 title claims abstract description 34
- 238000012876 topography Methods 0.000 claims abstract description 25
- 238000011156 evaluation Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 43
- 230000005540 biological transmission Effects 0.000 claims description 29
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 23
- 238000005299 abrasion Methods 0.000 claims description 22
- 238000004088 simulation Methods 0.000 claims description 18
- 239000012212 insulator Substances 0.000 claims description 17
- 230000015556 catabolic process Effects 0.000 claims description 16
- 238000013210 evaluation model Methods 0.000 claims description 15
- 239000012530 fluid Substances 0.000 claims description 15
- 238000010276 construction Methods 0.000 claims description 14
- 230000004907 flux Effects 0.000 claims description 13
- 238000009413 insulation Methods 0.000 claims description 13
- 239000004020 conductor Substances 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 11
- 238000005859 coupling reaction Methods 0.000 claims description 11
- 230000005684 electric field Effects 0.000 claims description 8
- 239000000725 suspension Substances 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 230000017525 heat dissipation Effects 0.000 claims description 5
- 101000867849 Rattus norvegicus Voltage-dependent R-type calcium channel subunit alpha-1E Proteins 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 239000011248 coating agent Substances 0.000 claims description 2
- 238000000576 coating method Methods 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000003746 surface roughness Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- SDIXRDNYIMOKSG-UHFFFAOYSA-L disodium methyl arsenate Chemical compound [Na+].[Na+].C[As]([O-])([O-])=O SDIXRDNYIMOKSG-UHFFFAOYSA-L 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 239000013013 elastic material Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000007716 flux method Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明提供一种电力设备气象灾害分析预警方法,属于电力设备气象灾害分析预警技术领域;所要解决的技术问题为:提供一种电力设备气象灾害分析预警方法;解决该技术问题采用的技术方案为:针对待预警区域地形特征,采集该区域气象及地形数据,构建典型微地形区电力设备逐站逐杆塔模型;基于降雨、风速、温度、湿度、气压的气象要素分类,构建气象要素对电力设备影响评估模型,定义电力设备灾害风险评估指数,构建电力设备气象灾害风险评估指标的模型;基于指标模型对当前区域内的电力设备气象灾害风险进行分析,并判断是否发出预警;本发明应用于电力设备气象灾害分析预警。
Description
技术领域
本发明提供一种电力设备气象灾害分析预警方法,属于电力设备气象灾害分析预警技术领域。
背景技术
目前随着新能源电力系统的不断投入建设,风电和光伏发电系统在未来会逐步成为供电系统的主体,但采用的风电、光伏等发电方式受气象条件的影响更大,如遇狂风暴雨、长时间阴雨天气、以及水灾内涝灾害都将造成新能源电力系统出现运行故障,因此对自然灾害因素的评估分析及预警,成为新能源电力系统能否正常运行的关键。
目前电网监控自然气象灾害提供有专用的监测预警平台,但在使用过程中仍然存在以下缺点与不足:气象预报数据主要由气象部门给出,预报主要针对固定的城镇区域,缺乏对不同区域、不同电力设备(不同输电线路、变电设备、配电设备)的危害影响程度的量化分析预报,预报信息针对性、指导性弱,将会导致设备所处位置的灾害预警失当,造成不必要停电;此外,由于外部气象数据与电网气象监测数据融合度低,不能充分利用电网安装的气象监测数据开展针对小空间尺度、小时间尺度的预报,对于特定区域的气象模拟与分析缺失,将导致预报精度较低。
发明内容
本发明为了克服现有技术中存在的不足,所要解决的技术问题为:提供一种电力设备气象灾害分析预警方法。
为了解决上述技术问题,本发明采用的技术方案为:一种电力设备气象灾害分析预警方法,包括如下的分析预警步骤:
步骤一:针对待预警区域地形特征,采集该区域气象及地形数据,构建典型微地形区电力设备逐站逐杆塔模型;
步骤二:基于降雨、风速、温度、湿度、气压的气象要素分类,构建气象要素对电力设备影响的评估模型,定义电力设备灾害风险评估指数,构建电力设备气象灾害风险评估指标的模型,其中:
针对降雨要素的评估模型构建内容包括:对流体动力模型和二维水动力模型的构建;
针对风速要素的评估模型构建内容包括:对输电塔进行有限元建模、对绝缘子串进行有限元建模、对导线进行有限元建模,并对金具磨损模型进行构建;
针对温度要素的评估模型构建内容包括:采用流体温度场耦合计算方法对配电变压器内外热交换进行仿真分析;
针对湿度、气压要素的评估模型构建内容包括:对空气电击穿模型进行构建;
步骤三:基于指标模型对当前区域内的电力设备气象灾害风险进行分析,并判断是否发出预警。
所述步骤二中构建流体动力模型的具体方法为:
将三维空间内整体连续的流体计算域分割为计算单元,以有限单元为最小计算尺度,求解微元上的流体力学控制方程,最终获取流场特征,具体采用有限体积法对数值进行处理:
首先对计算区域离散化并划分网格,并使每个网格点周围有一个互不重复的控制体积;
再对每一个控制体积积分,得出一组离散方程;
区域离散化过程结束之后获得几何要素:节点、控制体积、界面、网格线,构建流体动力模型。
所述步骤二中构建二维水动力模型的具体方法为:
针对雨水地表漫流及水流演进采用二维浅水方程SWEs进行模拟,矢量形式的表达式为:
;
其中:
;
式中,h为水深,q x 和q y 分别为x、y方向的单宽流量,g为重力加速度,u、v分别为x、y方向的流速,F和G分别为x、y方向的通量矢量,S为源项矢量,z b 为河床底高程,C f 为床面摩擦系数,C f =gn 2/h 1/3,其中n为曼宁系数;
该模型采用动力波法对水流演进过程进行模拟计算,先采用非结构网格剖分计算域,随后应用有限体积法,在单元i内,控制方程SWEs的积分表达式为:
;
式中:Ω为控制体i的体积,应用高斯散度定理,则式中通量项的面积分可以用线积分表示为:
;
式中:Г为控制体i的边界,n Г为边界Г所对应的外法线方向的单位向量,S b 为底坡源项,S f 为摩阻源项;
相应界面的通量向量F(Q)·n Г表示为如下的三角网格单元:
;
式中:n x 、n y 分别为n Г在x和y方向上的单位向量分量;
在该三角网格单元内,通量F(Q)·n Г线积分表达式为:
;
式中:k为单元边的编号,l k 为第i个单元格第k个边的边长,通量项采用基于Godunov格式的黎曼求解器求解。
所述步骤二中对输电塔的塔架进行有限元建模的具体方法为:
针对空间刚架模型进行建模,具体对塔架的所有杆件均采用空间梁单元模拟,首先在CAD建模环境中建立塔架的简化模型,然后导入有限元软件中对单元网格进行划分,并定义各个组成构件的截面特性和相关材料信息,最后输出建模结果。
所述步骤二中对绝缘子串进行有限元建模的具体方法为:
采用RBE2运动耦合刚性连接来模拟绝缘子串,由RBE2运动耦合刚性连接模拟安装绝缘子串受到风荷载或者线路不均匀覆冰情况被拉偏的位移参考点,以及在线路中产生的不平衡张力。
所述步骤二中对导线进行有限元建模的具体方法为:
将输电塔之间悬挂的导线定义为悬链线,定义点A、B分别为导线挂点,点O为悬链线弧垂最低点;
对整档导线中OC段进行受力分析,该段导线的垂直荷载为,C点沿导线切线方向的应力为/>,水平应力为/>,应力方向与水平夹角成θ倾角,该段导线的受力平衡方程为:
;
;
上述两式比值为该段导线任一点的切线斜率,即:
;
将上式对x微分得到:
;
将上式积分后整理得:
;
再将上式分离变量后对两端积分,并代入相应参数,求得坐标原点位于左侧悬挂点的导线悬链线方程:
;
式中:h为悬挂点高差;
为比载,架空导线单位长度、单位截面上承受的荷载;
为档距,相邻两悬挂点间垂直于荷载方向的投影距离;
为两侧悬挂点等高下的档内悬链线线长,计算公式为:
;
可得到等高悬点架空导线悬链线的计算公式为:
;
基于架空导线档距,垂直比载和水平应力,即可绘制出悬链线曲线,构建导线模型。
所述步骤二中对金具磨损模型进行构建的具体方法为:
采用Archard磨损模型,根据磨损率是局部接触压力和滑移距离的线性函数进行计算,计算公式为:
;
式中,V为磨损体积;K为磨损系数;H为材料硬度;P为接触压力;L为接触距离;
针对不同风载荷,结合输电线路舞动模型得到不同风载荷情况下金具的应力应变分布情况,然后基于Archard磨损模型构建金具磨损模型,模拟其磨损过程,进而提取金具磨损指数,得到金具的磨损区域及磨损程度。
所述步骤二中对配电变压器内外热交换进行仿真分析的具体方法为:
通过有限体积法模拟变电器的能量产生与传递过程,计算得到不同环境温度及风速条件下,变压器内部与外部的热交换速率,提取设备过热指数,得到不同温度及风速对变压器散热效率的影响程度。
所述步骤二中对空气电击穿模型进行构建的具体方法为:
采用玻尔兹曼解析法求解不同湿度不同压力空气的电击穿特性,分析湿度及压力对空气的电击穿特性的影响,具体针对不同湿度及不同气压计算出对应的临界电场强度,评估湿度及气压对空气绝缘特性的影响,提取空气绝缘指数,评估湿度及气压对空气绝缘及击穿特性的影响程度。
本发明相对于现有技术具备的有益效果为:本发明提供的电力设备气象灾害分析预警方法,具体通过选取典型微地形区,搭建气象要素实时同步监测系统,实现气象要素的实时准确监测,同时,基于计算流体动力学技术,实现典型微地形区电力设备逐站逐杆塔气象精准预测,并基于精准的气象预测结果,构建气象要素对电力设备影响评估模型,形成对电力设备气象灾害风险的判断规则,提取电力设备灾害风险关键指标,并结合自动化技术,实现对电网不同运维需求的定制化电力气象服务,进而实现电网气象灾害预警的精准化与智能化。
附图说明
下面结合附图对本发明做进一步说明:
图1为本发明的步骤流程图;
图2为本发明采用有限体积法对二维结构网格区域进行离散的状态图;
图3为本发明绝缘子串连接结构示意图;
图4为本发明不等高挂点架空导线结构示意图;
图5为本发明对变压器流体温度场分布计算的步骤流程图。
具体实施方式
如图1所示,本发明在充分采集自然环境与电网设备数据特征的基础上,构建了电网电力设备气象灾害风险关键指标体系,针对电网不同环节提供定制化、精准化的电力气象服务,变被动监测为主动预警,解决现有气象预报对电网针对性不强,缺乏量化指标的问题,全面提高电力设备应对气象灾害的应急反应和处置能力,最大限度预防和减轻气象灾害对电网安全运行的影响。
为实现上述目的,本发明采用的分析预警方法主要包括:
针对待预警区域地形特征,采集该区域气象及地形数据,构建典型微地形区电力设备逐站逐杆塔模型;
基于降雨、风速、温度、湿度、气压的气象要素分类,构建气象要素对电力设备影响的评估模型,定义电力设备灾害风险评估指数,构建电力设备气象灾害风险评估指标的模型;
将气象灾害与电力设备状态信息进行融合,基于指标模型对当前区域内的电力设备气象灾害风险进行分析,并判断是否发出预警。
其中,构建微地形区气象要素分析模型的原因在于:受限于计算性能及模式参数化方案的限制,气象数值模式大多对地形进行了网格剖分和平滑处理,只考虑了宏观地形对气象要素的作用,忽略了对微地形对气象参数的影响,模拟的云微物理特征和垂直运动特征均与实际情况存在一定差别,这种差别对定点、定量降水和天空状况、地面气温、湿度等要素的准确预报都会产生影响,因此预报的气象要素与实际环境下的气象要素存在较大误差,使得基于气象预报得到的电力气象灾害风险指标存在较大的不确定性,获得高精度、精细化的微地形区气象场是实现覆冰厚度精确分析的关键;
构建电力设备气象要素影响评估模型的原因在于:电力设备气象灾害风险指标体系的建立需以电力设备气象要素影响评估模型作为基础,但各类气象要素对电力设备的影响分析模型众多,且涉及大量不确定因素;如变电站和配电设备积水模型的建立中需考虑水流演进过程中干湿交替模拟精度不足,地形因素影响,地表漫流模拟精度不足等各类问题,因此需建立统一的评估模型对不确定因素进行分析处理,得到评估结果。
构建流体动力模型的方法为:
计算流体动力学(Computational Fluid Dynamic,CFD),是基于计算机平台求解流体流动、传热效应以及流场相关参数的分析方法和手段。CFD的基本原理为:将三维空间内整体连续的流体计算域,分割为极小的计算单元,以有限单元为最小计算尺度,求解微元上的流体力学控制方程,最终获取流场特征。
CFD数值计算方法主要是对计算域内单元的离散化,进而对控制方程进行离散化,通过有限差分法(FDM)、有限元法(FEM)和有限体积法(FVM)三种离散方法对控制方程进行离散求解;
本发明采用有限体积法FVM对数值进行处理,首先将计算区域离散化并划分网格,并使每个网格点周围有一个互不重复的控制体积;再对每一个控制体积积分,从而得出一组离散方程;如图2所示为直角坐标系下有限体积法的二维网格离散例子,从其二维结构网格区域离散示意图中可以看出,区域离散化过程结束之后可以获得如下四种几何要素:
节点(Node):待求解未知量的位置;
控制体积(Control volume):遵循控制方程或守恒方程的最小单位;
界面(Face):相邻节点间控制容积的分界面;
网格线(Grid line):连接相邻节点所形成的曲线簇。
有限体积法相比于有限差分法,对离散网格的正交性没有要求,因此可以大量使用非结构化网格求解复杂的流场模型问题,在实际计算时因网格计算量小,而对内存空间的占用较小,有条件更精细地划分流场细节;利用计算流体动力学仿真,进行微地形区空气流动特征、水汽分布特征的仿真分析,从而实现微气象数据与大范围气象数据的交互分析,为开展输电线路覆冰、舞动特征分析奠定基础。
构建二维水动力模型的方法为:
水动力学模型是以平面二维水动力学模型为基础,将计算区域上边界产生的洪水过程与区间的产、汇流过程,分别按上开边界条件和面源,以沿程旁侧入汇形式结合起来融入二维水动力学模型。通过计算区域内水域动边界的自动跟踪、调整、合理分配,解决各子区间内的产、汇流问题,并通过适宜的穿渠建筑物泄流曲线或泄流公式控制中边界过水问题。全面、准确地模拟计算区域内在不同标准、不同工程规模情况下洪水的纵、横向传播及串流状况。
针对雨水地表漫流及水流演进采用二维浅水方程SWEs进行模拟,矢量形式的表达式为:
;
其中:
;
式中,h为水深,q x 和q y 分别为x、y方向的单宽流量,g为重力加速度,u、v分别为x、y方向的流速,F和G分别为x、y方向的通量矢量,S为源项矢量,z b 为河床底高程,C f 为床面摩擦系数,C f =gn 2/h 1/3,其中n为曼宁系数。
该模型采用动力波法对水流演进过程进行模拟计算,先采用非结构网格剖分计算域,随后应用有限体积法,在单元i内,控制方程SWEs的积分表达式为:
;
式中:Ω为控制体i的体积,应用高斯散度定理,则式中通量项的面积分可以用线积分表示为(其他项在网格单元内可认为恒定):
;
式中:Г为控制体i的边界,n Г为边界Г所对应的外法线方向的单位向量,S b 为底坡源项,S f 为摩阻源项。
相应界面的通量向量F(Q)·n Г可以表示为如下的三角网格单元:
;
式中:n x 、n y 分别为n Г在x和y方向上的单位向量分量。
在该三角网格单元内,通量F(Q)·n Г线积分表达式为:
;
式中:k为单元边的编号,l k 为第i个单元格第k个边的边长。通量项采用基于Godunov格式的黎曼求解器求解。
在网格单元内,界面通量通过HLLC格式的近似黎曼求解器计算,通过静水重构来修正干湿边界处负水深问题,底坡源项使用底坡通量法进行处理,摩阻源项使用半隐式法计算以此提高模型稳定性。
进一步的,本发明在进行多源数据融合的气象灾害预测时,首先建立基于三维变分同化的数值天气预报模型,实现对监测区域气象要素(温度、湿度、风速、风向、降水等)1×1km区域的精细化预报。针对典型微地形区,抽取微地形定量表征参数,构建典型微地形3D模型,利用1×1km区域的数值预报同化资料作为输入场,结合电力设备附近气象站历史公共气象数据、三维地形等数据,基于CFD对不同空气路径影响下,微地形区空气流动特征、气温及水汽分布特征进行敏感性分析,建立微地形区温度、湿度、降水量及相态、风速、风向、光照等气象要素预报模型,实现典型微地形区电力设备逐站逐杆塔气象精准预测。
针对微地形地理特征指标抽取和建模时,通过提取电力设备所在区域的微地形特征参数(如地表粗糙度、坡度、坡向、湖泊距离及下垫面类型),制定符合局地的微地形参数化方案,为其数字化提供数据支撑,可采用激光雷达点云扫描、倾斜摄影等手段,准确获取微地形定量表征参数,构建典型微地形三维模型,准确掌握局地的微地形定量表征参数。
针对电网典型微地形区域建立观测站,本发明根据现场情况确定部署自动气象站或X波段雷达,并在迎风方向上,同时安装声雷达,并结合辐射、湍流观测系统、土壤温湿度和热流观测设备等,进行全方位的微地形条件下的微气象观测试验,全面采集数据;围绕观测点进行长时间的微气象要素监测,获取完整的时间序列数据,包括春秋季大风及冬季冰冻雨雪天气等数据。通过分析观测数据,获取微气象要素的时空变化规律及微地形对气象要素的影响。
针对电网1×1km数值预报模拟时,结合气象卫星资料、雷达资料、区域地形数据及水体数据,利用三维变分同化、资料预处理、观测资料同化、气象卫星资料同化、雷达资料同化和气溶胶资料等手段,研究气象预报数据与微地形地理信息、微气象监测数据、卫星资料、多普勒雷达资料等的同化方法,得到预报结果更准确的气象数值预测模型,实现对特定区域1×1km气象要素的短期定点预报。
同时,本发明为研究微地形对气象的影响,结合微地形3D模型,利用1×1km数值区域预报同化资料作为输入场,使用CFD软件对不同冷空气路径影响下,观测点塔杆附近的空气流动特征、气温及水汽分布特征进行敏感性试验并建模,研究基于1×1km气象同化资料的气象要素降尺度方法,使其能准确预报设备附近气象要素的时空变化特征,并确保与气象监测数据尽可能吻合;并通过研究不同微地形指标(如地表粗糙度、坡度、坡向、湖泊距离及下垫面类型)对微气象要素的影响,在此基础上对微气象要素预测模型进行订正,使其适应不同微地形条件下微气象要素的预报。此外,将该预测模型与数值同化预报系统相耦合,利用系统输出的气象要素对同一时间段内不同微地形条件下气象要素的变化特征进行模拟,并与观测值进行对比,利用人工智能方法对误差进行校正。
进一步的,本发明针对风速要素的评估模型构建需要基于对输电塔、绝缘子串、导线进行有限元建模,并对因风力导致金具磨损的模型进行构建,输电塔线体系是主要由输电塔、输电线和绝缘子串等构件组成的耦合体系,因此建立一个合理有效的塔线耦合体系有限元模型对后续各种工况的动力响应分析和动力特性研究非常重要。
(1)输电塔有限元建模:
输电塔的塔架由各种不同型号的等边角钢组成,本发明具体针对空间刚架模型进行建模,具体对塔架的所有杆件均采用空间梁单元模拟,首先在CAD建模环境中建立塔架的简化模型,然后导入有限元软件中对单元网格进行划分,并定义各个组成构件的截面特性和相关材料信息,最后输出建模结果。
(2)绝缘子串有限元建模:
如图3所示,绝缘子串作为输电线和输电塔之间的连接构件,在二者之间起着绝缘的作用,绝缘子串两端分别铰接于输电塔横担和导线上,在自重荷载作用下一般是垂直向下的,遇到风荷载或者线路不均匀覆冰等情况时会被拉偏,并在线路中产生不平衡张力,本发明采用RBE2运动耦合刚性连接来模拟绝缘子串,由RBE2运动耦合刚性连接模拟安装绝缘子串受到风力影响的位移参考点。
(3)导线有限元建模:
由于架空输电线路的档距很大,输电塔挂点间的导线长度远大于其截面尺寸,因此可忽略导线刚性对导线悬挂于空中的几何形状的影响,将导线看作一根处处铰接的柔性链条,在进行输电导线力学性能分析前,需要假设导线只承受轴向拉力,不能承受压力与弯矩,是一种理想的柔性结构,并且作用于导线上的荷载沿导线均匀分布并指向同一方向。
通过上述假设,可将输电塔之间悬挂的导线看作一条理想的悬链线,如图4所示为不等高挂点架空导线示意图及受力图,图中点A、B分别为导线挂点,点O为悬链线弧垂最低点;
对图中整档导线中OC段进行单独的受力分析,该段导线的垂直荷载为,C点沿导线切线方向的应力为/>,水平应力为/>,应力方向与水平夹角成θ倾角,该段导线的受力平衡方程如下:
;
;
上述两式比值为该段导线任一点的切线斜率,即:
;
将上式对x微分得到:
;
将上式积分后整理得:
;
再将上式分离变量后对两端积分,并代入相应参数即可求得坐标原点位于左侧悬挂点的导线悬链线方程:
;
式中:h为悬挂点高差;
为比载,架空导线单位长度、单位截面上承受的荷载;
为档距,相邻两悬挂点间垂直于荷载方向的投影距离;
为两侧悬挂点等高下的档内悬链线线长,计算公式为:
;
可得到等高悬点架空导线悬链线的计算公式为:
;
故只需要知道架空导线档距,垂直比载和水平应力,即可绘制出悬链线曲线;
本发明采用梁单元来模拟导线,针对导线抗拉不抗压的特性,利用软件中的亚弹性材料来定义其拉压方向的应力应变关系曲线来实现;当应变ε大于0时,弹性模量E取导线的实际弹性模量;应变ε小于0时,弹性模量E取为0。
最后将输电塔、绝缘子串和导线等有限元模型组合起来,就可得到输电塔线体系的耦合有限元模型。
上述杆塔、导线、绝缘子串建模计算采用了非线性有限元瞬态动力学分析算法。瞬态动力学分析又称为时间历程分析,是用于确定结构在任何随时间变化的载荷作用下的动力学响应的分析方法,可以确定结构在静载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力等。
对于空间有限元离散,采用,/>,/>(/>)描述物体内各点在时间0、时间t、时间/>的位形坐标,/>,/>(/>)表示各质点在时间t和时间/>的位移,采用平衡方程,基于虚位移原理进行计算,计算公式为:
;
式中:
为/>时刻应力;
为/>时刻应变;
为/>时刻的外力虚功。
由于导线脱冰跳跃过程是大位移、小变形的非线性运动过程,因此需要采用相应的空间运动描述和定义,如采用TL或者UL坐标等以及相应的度量,可得瞬态动力学分析求解的基本的运动方程为:
;
式中:
表示质量矩阵;
表示阻尼矩阵;
表示刚度矩阵;
表示节点加速度向量;
表示节点速度向量;
表示节点位移向量;
表示t时刻的载荷向量。
在任意给定的时间t,上式可以看做是一系列考虑了惯性力和阻尼力的静力学平衡方程,对时间积分,具体采用Newmark方法,设定:
;
;
要求在间隔内,满足线性假设:
;
间隔内的加速度为:
;
时间点处的位移解u,可以通过求解:
;
其中加速度可由下式:
;
代入,则得到从、/>、/>计算/>的公式:
;
重复迭代以上计算过程,即可确定给定时间内的导线各点脱冰跳跃的运动轨迹。
(4)对金具磨损的模型进行构建:
该步骤主要涉及金具磨损模型及磨损指数提取技术研究,由于金具磨损一般发生在接触应力大,并发生较大相对滑动的区域,为了模拟金具之间的磨损情况,采用Archard磨损模型,磨损率是局部接触压力和滑移距离的线性函数,计算公式为:
;
式中,V为磨损体积;K为磨损系数;H为材料硬度;P为接触压力;L为接触距离。
针对不同风载荷,结合输电线路舞动模型得到不同风载荷情况下金具的应力应变分布情况,基于 Archard 磨损模型针对应力大且接触并容易发生滑移的不稳构建金具磨损模型,模拟其磨损过程,进而提取金具磨损指数,得到金具的磨损区域及磨损程度,对相关部位重点关注与加强防护。
进一步的,本发明采用流体温度场耦合计算方法对配电变压器内外热交换进行仿真分析,主要涉及对散热模型构建及散热指数提取技术研究;
由于配电变压器季节性短时过载问题突出,在夏季用电高峰期,变压器的负载率出现短时激增,严重时超过额定值,内部温度迅速上升,可能引起变压器烧毁,威胁供电安全与可靠性。
本发明采用流体-温度场耦合计算方法对配电变压器内外热交换进行仿真分析,通过有限体积法模拟变电器的能量产生与传递过程,计算过程如图5所示,计算得到不同环境温度及风速条件下,变压器内部与外部的热交换速率,提取设备过热指数,得到不同温度及风速对变压器散热效率的影响程度,结合气象预报,对变压器过热风险做出动态预警。
进一步的,本发明针对湿度、气压要素的评估模型构建主要涉及对空气电击穿模型的构建;
衡量气体电击穿特性的一个重要指标是临界电场强度,临界电场强度指的是在气体中电子崩处于继续发展和逐渐消失的临界状态时对应的电场强度。当外界电场强度高于临界电场强度时,电子崩能够迅速发展,气体有被击穿的可能,因此在实际的绝缘设计中,应使得气体间隙所能承受的最大场强低于临界电场强度。
本发明主要采用玻尔兹曼解析法求解不同湿度不同压力空气的电击穿特性,进而分析湿度及压力对空气的电击穿特性的影响,具体针对不同湿度及不同气压计算出对应的临界电场强度,评估湿度及气压对空气绝缘特性的影响,提取空气绝缘指数,评估湿度及气压对空气绝缘及击穿特性的影响程度。
由于空气中的水分附着在绝缘材料表面,使电气设备的绝缘电阻降低,特别是使用年限较长的设备,由于内部有积尘吸附水分,潮湿程度将更严重,绝缘电阻更低。设备的泄露电流大大增加,甚至造成绝缘击穿,产生事故。针对这一问题,本发明通过提取设备闪络指数,评估湿度对重污秽区设备发生污闪的影响程度,并结合气象预报,对未来设备发生闪络风险做出动态预警。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (2)
1.一种电力设备气象灾害分析预警方法,其特征在于:包括如下的分析预警步骤:
步骤一:针对待预警区域地形特征,采集该区域气象及地形数据,构建典型微地形区电力设备逐站逐杆塔模型;
步骤二:基于降雨、风速、温度、湿度和气压的气象要素分类,构建气象要素对电力设备影响的评估模型,定义电力设备灾害风险评估指数,构建电力设备气象灾害风险评估指标的模型,其中:
针对降雨要素的评估模型构建内容包括:对流体动力模型和二维水动力模型的构建;
针对风速要素的评估模型构建内容包括:对输电塔进行有限元建模、对绝缘子串进行有限元建模和对导线进行有限元建模,并对金具磨损模型进行构建;
针对温度要素的评估模型构建内容包括:采用流体温度场耦合计算方法对配电变压器内外热交换进行仿真分析;
针对湿度和气压要素的评估模型构建内容包括:对空气电击穿模型进行构建;
所述步骤二中构建流体动力模型的具体方法为:
将三维空间内整体连续的流体计算域分割为计算单元,以有限单元为最小计算尺度,求解微元上的流体力学控制方程,最终获取流场特征,具体采用有限体积法对数值进行处理:
首先对计算区域离散化并划分网格,并使每个网格点周围有一个互不重复的控制体积;
再对每一个控制体积积分,得出一组离散方程;
区域离散化过程结束之后获得几何要素:节点、控制体积、界面和网格线,构建流体动力模型;
所述步骤二中构建二维水动力模型的具体方法为:
针对雨水地表漫流及水流演进采用二维浅水方程SWEs进行模拟,矢量形式的表达式为:
;
其中:
;
式中,h为水深,q x 和q y 分别为x、y方向的单宽流量,g为重力加速度,u、v分别为x、y方向的流速,F和G分别为x、y方向的通量矢量,S为源项矢量,z b 为河床底高程,C f 为床面摩擦系数,C f = gn 2/h 1/3,其中n为曼宁系数;
该模型采用动力波法对水流演进过程进行模拟计算,先采用非结构网格剖分计算域,随后应用有限体积法,在单元i内,控制方程SWEs的积分表达式为:
;
式中:Ω为控制体i的体积,应用高斯散度定理,则式中通量项的面积分可以用线积分表示为:
;
式中:Г为控制体i的边界,n Г为边界Г所对应的外法线方向的单位向量,S b 为底坡源项,S f 为摩阻源项;
相应界面的通量向量F(Q)·n Г表示为如下的三角网格单元:
;
式中:n x 、n y 分别为n Г在x和y方向上的单位向量分量;
在该三角网格单元内,通量F(Q)·n Г线积分表达式为:
;
式中:k为单元边的编号,l k 为第i个单元格第k个边的边长,通量项采用基于Godunov格式的黎曼求解器求解;
所述步骤二中对输电塔的塔架进行有限元建模的具体方法为:
针对空间刚架模型进行建模,具体对塔架的所有杆件均采用空间梁单元模拟,首先在CAD建模环境中建立塔架的简化模型,然后导入有限元软件中对单元网格进行划分,并定义各个组成构件的截面特性和相关材料信息,最后输出建模结果;
所述步骤二中对绝缘子串进行有限元建模的具体方法为:
采用RBE2运动耦合刚性连接来模拟绝缘子串,由RBE2运动耦合刚性连接模拟安装绝缘子串受到风荷载或者线路不均匀覆冰情况被拉偏的位移参考点,以及在线路中产生的不平衡张力;
所述步骤二中对导线进行有限元建模的具体方法为:
将输电塔之间悬挂的导线定义为悬链线,定义点A、B分别为导线挂点,点O为悬链线弧垂最低点;
对整档导线中OC段进行受力分析,该段导线的垂直荷载为,C点沿导线切线方向的应力为/>,水平应力为/>,应力方向与水平夹角成θ倾角,该段导线的受力平衡方程为:
;
;
上述两式比值为该段导线任一点的切线斜率,即:
;
将上式对x微分得到:
;
将上式积分后整理得:
;
再将上式分离变量后对两端积分,并代入相应参数,求得坐标原点位于左侧悬挂点的导线悬链线方程:
;
式中:h为悬挂点高差;
为比载,架空导线单位长度、单位截面上承受的荷载;
为档距,相邻两悬挂点间垂直于荷载方向的投影距离;
为两侧悬挂点等高下的档内悬链线线长,计算公式为:
;
可得到等高悬点架空导线悬链线的计算公式为:
;
基于架空导线档距,垂直比载和水平应力,即可绘制出悬链线曲线,构建导线模型;
所述步骤二中对金具磨损模型进行构建的具体方法为:
采用Archard磨损模型,根据磨损率是局部接触压力和滑移距离的线性函数进行计算,计算公式为:
;
式中,V为磨损体积;K为磨损系数;H为材料硬度;P为接触压力;L为接触距离;
针对不同风载荷,结合输电线路舞动模型得到不同风载荷情况下金具的应力应变分布情况,然后基于Archard磨损模型构建金具磨损模型,模拟其磨损过程,进而提取金具磨损指数,得到金具的磨损区域及磨损程度;
所述步骤二中对空气电击穿模型进行构建的具体方法为:
采用玻尔兹曼解析法求解不同湿度不同压力空气的电击穿特性,分析湿度及压力对空气的电击穿特性的影响,具体针对不同湿度及不同气压计算出对应的临界电场强度,评估湿度及气压对空气绝缘特性的影响,提取空气绝缘指数,评估湿度及气压对空气绝缘及击穿特性的影响程度;
步骤三:基于指标模型对当前区域内的电力设备气象灾害风险进行分析,并判断是否发出预警。
2.根据权利要求1所述的一种电力设备气象灾害分析预警方法,其特征在于:所述步骤二中对配电变压器内外热交换进行仿真分析的具体方法为:
通过有限体积法模拟变电器的能量产生与传递过程,计算得到不同环境温度及风速条件下,变压器内部与外部的热交换速率,提取设备过热指数,得到不同温度及风速对变压器散热效率的影响程度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311392865.2A CN117114428B (zh) | 2023-10-25 | 2023-10-25 | 一种电力设备气象灾害分析预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311392865.2A CN117114428B (zh) | 2023-10-25 | 2023-10-25 | 一种电力设备气象灾害分析预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117114428A CN117114428A (zh) | 2023-11-24 |
CN117114428B true CN117114428B (zh) | 2024-01-30 |
Family
ID=88807815
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311392865.2A Active CN117114428B (zh) | 2023-10-25 | 2023-10-25 | 一种电力设备气象灾害分析预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117114428B (zh) |
Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105741045A (zh) * | 2016-02-05 | 2016-07-06 | 中国水利水电科学研究院 | 一种洪水风险动态分析与展示系统及方法 |
WO2016129715A1 (ko) * | 2015-02-10 | 2016-08-18 | 주식회사 주빅스 | 환경 재난, 재해 조기감지를 위한 대기질 예측 및 관리 시스템 |
CN106099756A (zh) * | 2016-08-15 | 2016-11-09 | 国网山西省电力公司电力科学研究院 | 基于微气象预测的输电线路覆冰灾害分析装置 |
CN107292497A (zh) * | 2017-06-05 | 2017-10-24 | 国网陕西省电力公司电力科学研究院 | 基于层次分析‑熵权组合的输电线路风偏闪络风险评估方法 |
CN109214604A (zh) * | 2018-11-02 | 2019-01-15 | 国网湖南省电力有限公司 | 基于Adaboost的输变电设备的雾霾灾害风险预测方法及系统 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109919409A (zh) * | 2019-01-08 | 2019-06-21 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种变电站或换流站耦连回路地震风险评估方法 |
CN111898303A (zh) * | 2020-08-05 | 2020-11-06 | 苏州大圜科技有限公司 | 基于气象预报及水动力模拟的流域水位及内涝预报方法 |
CN112541636A (zh) * | 2020-12-17 | 2021-03-23 | 国网湖南省电力有限公司 | 输电线路覆冰风险预警方法、装置、介质和电子设备 |
CN113642206A (zh) * | 2021-07-06 | 2021-11-12 | 岭澳核电有限公司 | 核电站用boss头焊缝缺陷评估方法、装置、设备及介质 |
CN113987810A (zh) * | 2021-10-29 | 2022-01-28 | 中铁二院工程集团有限责任公司 | 一种流域尺度下暴雨洪水的数值预报方法 |
CN115657028A (zh) * | 2022-11-09 | 2023-01-31 | 上海无线电设备研究所 | 一种基于GEO-UAV BiSAR的海面晃动舰船复合散射建模方法 |
CN115983093A (zh) * | 2022-11-29 | 2023-04-18 | 国网新疆电力有限公司经济技术研究院 | 一种输电线路塔基地质灾害预测方法和装置 |
CN116090677A (zh) * | 2023-04-10 | 2023-05-09 | 湖南大学 | 考虑电力-通信-交通网络耦合的空地应急资源规划方法 |
CN116205342A (zh) * | 2022-12-30 | 2023-06-02 | 国网河南省电力公司电力科学研究院 | 基于精细化预测的电力气象预警方法 |
CN116384733A (zh) * | 2023-03-10 | 2023-07-04 | 国网河南省电力公司电力科学研究院 | 一种强对流天气下基于气象雷达的输电线路风险预警方法 |
CN116564045A (zh) * | 2023-03-27 | 2023-08-08 | 国网河南省电力公司社旗县供电公司 | 一种基于电力微气象的输电线路风偏灾害预警装置 |
CN116799946A (zh) * | 2023-04-26 | 2023-09-22 | 国网山东省电力公司电力科学研究院 | 一种基于风险调度的电网全息智能管控方法及平台 |
CN116911135A (zh) * | 2023-07-24 | 2023-10-20 | 武汉理工大学 | 计及波浪下砰击载荷的非线性水弹性时域计算评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11402362B2 (en) * | 2017-06-19 | 2022-08-02 | Nanjing Institute Of Geography | Stereoscopic monitoring and data mining system and method for harmful lake cyanobacteria bloom |
-
2023
- 2023-10-25 CN CN202311392865.2A patent/CN117114428B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016129715A1 (ko) * | 2015-02-10 | 2016-08-18 | 주식회사 주빅스 | 환경 재난, 재해 조기감지를 위한 대기질 예측 및 관리 시스템 |
CN105741045A (zh) * | 2016-02-05 | 2016-07-06 | 中国水利水电科学研究院 | 一种洪水风险动态分析与展示系统及方法 |
CN106099756A (zh) * | 2016-08-15 | 2016-11-09 | 国网山西省电力公司电力科学研究院 | 基于微气象预测的输电线路覆冰灾害分析装置 |
CN107292497A (zh) * | 2017-06-05 | 2017-10-24 | 国网陕西省电力公司电力科学研究院 | 基于层次分析‑熵权组合的输电线路风偏闪络风险评估方法 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109214604A (zh) * | 2018-11-02 | 2019-01-15 | 国网湖南省电力有限公司 | 基于Adaboost的输变电设备的雾霾灾害风险预测方法及系统 |
CN109919409A (zh) * | 2019-01-08 | 2019-06-21 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种变电站或换流站耦连回路地震风险评估方法 |
WO2019238140A1 (zh) * | 2019-01-08 | 2019-12-19 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种变电站或换流站耦连回路地震风险评估方法 |
CN111898303A (zh) * | 2020-08-05 | 2020-11-06 | 苏州大圜科技有限公司 | 基于气象预报及水动力模拟的流域水位及内涝预报方法 |
CN112541636A (zh) * | 2020-12-17 | 2021-03-23 | 国网湖南省电力有限公司 | 输电线路覆冰风险预警方法、装置、介质和电子设备 |
CN113642206A (zh) * | 2021-07-06 | 2021-11-12 | 岭澳核电有限公司 | 核电站用boss头焊缝缺陷评估方法、装置、设备及介质 |
CN113987810A (zh) * | 2021-10-29 | 2022-01-28 | 中铁二院工程集团有限责任公司 | 一种流域尺度下暴雨洪水的数值预报方法 |
CN115657028A (zh) * | 2022-11-09 | 2023-01-31 | 上海无线电设备研究所 | 一种基于GEO-UAV BiSAR的海面晃动舰船复合散射建模方法 |
CN115983093A (zh) * | 2022-11-29 | 2023-04-18 | 国网新疆电力有限公司经济技术研究院 | 一种输电线路塔基地质灾害预测方法和装置 |
CN116205342A (zh) * | 2022-12-30 | 2023-06-02 | 国网河南省电力公司电力科学研究院 | 基于精细化预测的电力气象预警方法 |
CN116384733A (zh) * | 2023-03-10 | 2023-07-04 | 国网河南省电力公司电力科学研究院 | 一种强对流天气下基于气象雷达的输电线路风险预警方法 |
CN116564045A (zh) * | 2023-03-27 | 2023-08-08 | 国网河南省电力公司社旗县供电公司 | 一种基于电力微气象的输电线路风偏灾害预警装置 |
CN116090677A (zh) * | 2023-04-10 | 2023-05-09 | 湖南大学 | 考虑电力-通信-交通网络耦合的空地应急资源规划方法 |
CN116799946A (zh) * | 2023-04-26 | 2023-09-22 | 国网山东省电力公司电力科学研究院 | 一种基于风险调度的电网全息智能管控方法及平台 |
CN116911135A (zh) * | 2023-07-24 | 2023-10-20 | 武汉理工大学 | 计及波浪下砰击载荷的非线性水弹性时域计算评估方法 |
Non-Patent Citations (2)
Title |
---|
基于知识图谱的电力变压器故障预测方法研究;董理科等;《高压电器》;第58卷(第11期);第151-157页 * |
极端龙卷风环境下超高层施工附属设施安全性态分析与评定;鲁逸;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》(第01期);第C038-1591页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117114428A (zh) | 2023-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Elshaer et al. | Variations in wind load on tall buildings due to urban development | |
CN111651885A (zh) | 一种智慧型海绵城市洪涝预报方法 | |
Al-Yahyai et al. | Review of the use of numerical weather prediction (NWP) models for wind energy assessment | |
Boudia et al. | Wind resource assessment in Algeria | |
Tabrizi et al. | Performance and safety of rooftop wind turbines: Use of CFD to gain insight into inflow conditions | |
CN103673960B (zh) | 用于预测输电线路上的覆冰状态的方法和装置 | |
Xue et al. | Impact of transmission tower-line interaction to the bulk power system during hurricane | |
Berge et al. | Wind in complex terrain. A comparison of WAsP and two CFD-models | |
Mochida et al. | Up-scaling CWE models to include mesoscale meteorological influences | |
CN112685965B (zh) | 一种台风中输电杆塔风险监测方法及系统 | |
CN102968554B (zh) | 基于安全裕度的杆塔覆冰灾害风险预测方法 | |
Zaïdi et al. | Numerical simulations of the microscale heterogeneities of turbulence observed on a complex site | |
CN108090285A (zh) | 一种适用于复杂地形输电线路风灾监测的微气象观测布点方法 | |
CN105447770A (zh) | 电网气象监测数据应用于精细化气象预报的评估方法 | |
Sharma et al. | Application of a new method to develop a CFD model to analyze wind characteristics for a complex terrain | |
Gentle | Concurrent wind cooling in power transmission lines | |
Li et al. | Wind tunnel study on the morphological parameterization of building non-uniformity | |
Uchida et al. | New evaluation technique for WTG design wind speed using a CFD-model-based unsteady flow simulation with wind direction changes | |
CN103018793B (zh) | 一种极端气象条件下的架线的云预警系统和方法 | |
CN117114428B (zh) | 一种电力设备气象灾害分析预警方法 | |
Phillips et al. | Dynamic rating of overhead transmission lines over complex terrain using a large-eddy simulation paradigm | |
JP2006317460A (ja) | 物理量の等値線図の作成方法及び物理量の推定方法 | |
Wang et al. | A study of function-based wind profiles based on least squares method: A case in the suburbs of Hohhot | |
Jin et al. | Integration of atmospheric stability in wind resource assessment through multi-scale coupling method | |
Cheggaga et al. | A neural network solution for extrapolation of wind speeds at heights ranging for improving the estimation of wind producible |
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 |