CN109241693B - 测算风力发电机组防雷等级的方法 - Google Patents
测算风力发电机组防雷等级的方法 Download PDFInfo
- Publication number
- CN109241693B CN109241693B CN201811361780.7A CN201811361780A CN109241693B CN 109241693 B CN109241693 B CN 109241693B CN 201811361780 A CN201811361780 A CN 201811361780A CN 109241693 B CN109241693 B CN 109241693B
- Authority
- CN
- China
- Prior art keywords
- calculating
- lightning
- wind
- lightning protection
- generating set
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D80/00—Details, components or accessories not provided for in groups F03D1/00 - F03D17/00
- F03D80/30—Lightning protection
-
- 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Abstract
本发明提供一种测算风力发电机组防雷等级的方法,包括以下步骤:S1、根据气象部门的数据计算风场所在地年平均雷击密度;S2、根据回击电流强度对等效截面半径的影响,计算与风机接收相同雷击次数的等效截收面积;S3、根据风场所在地形选择合适的环境因素,并依据风力涡轮发电机系统防雷规范中的表述,计算风力发电机组的年平均直击雷数量;S4、根据步骤S3得到的风力发电机组的年平均直击雷数量,确定风力发电机组防雷等级。本发明提供的测算风力发电机组防雷等级的方法,便于得到更加科学的雷击等效截收面积,从而能够更加有效地为不同机型的风电机组选择合适的雷电防护等级,提高风电机组防雷保护设计施工过程中的经济型与有效性。
Description
技术领域
本发明涉及风力发电机组领域,具体涉及一种测算风力发电机组防雷等级的方法。
背景技术
对于一般建筑物而言,为了对避雷性能进行分析首先需要得到避雷设施的闪击距离或者吸引半径,随后使用数学方法计算得到等效截收面积,通过考虑当地的年雷击密度得到遭受雷击概率。因此,通过计算得到准确的闪击距离以及吸引半径是对建筑物雷击评估的关键数据,对合理规避由雷电造成的风险就显得尤为重要。
根据IEC 61400标准中的定义,在计算风力发电设备等效截收面积时可将其视为静止的高塔看待,高度为底座杆塔与叶片长度之和。在地面以高度值的3倍为长度形成的范围为等效截收面积。可以将其视为是折线法的变形,然而这种方法并不包含任何物理上的含义,对于高度较高的物体来说过分夸大了其等效截收面积,仅仅是为了满足标准中雷击风险评估而提出的折中方案。无法得到准确的等效截收面积。而准确的等效截收面积对于为不同机型的风电机组选择合适的雷电防护等级尤为重要。
发明内容
针对现有技术中的缺陷,本发明提供的测算风力发电机组防雷等级的方法,便于得到更加科学合理的雷击等效截收面积,从而能够更加有效地为不同机型的风电机组选择合适的雷电防护等级,提高风电机组防雷保护设计施工过程中的经济型与有效性。
本发明提供的一种测算风力发电机组防雷等级的方法,包括以下步骤:
S1、根据气象部门的数据计算风场所在地年平均雷击密度Ng;
S2、根据回击电流强度对等效截面半径的影响,计算与风机接收相同雷击次数的等效截收面积Ae;
S3、根据风场所在地形选择合适的环境因素Cd,并依据风力涡轮发电机系统防雷规范中的表述,计算风力发电机组的年平均直击雷数量Nd;
S4、根据步骤S3得到的风力发电机组的年平均直击雷数量Nd,确定风力发电机组防雷等级。
可选地,步骤S2包括以下步骤:
S21、根据IEC62305规范计算闪击距离S;
S22、根据闪击距离S和机组高度h,计算雷击区域在地面上投影的半径r与机组高度h之间的函数r(h);
S23、根据发电机组中塔杆的高度H、叶片的长度L以及旋转角度θ,计算出发电机组在旋转过程中的高度h(θ)和宽度W(θ),所述旋转角度θ为位于塔杆最顶端的叶片与竖直轴线之间的夹角;
S24、将步骤S23中得到的h(θ)代入步骤S22所得到的函数r(h)中,计算雷击区域在地面上投影的半径r;
S25、根据规则建筑物的等效截收面积计算方法,计算等效截收面积与旋转角度θ之间的函数Ae(θ);并将步骤S24得到半径r和步骤S23得到的宽度W(θ)代入Ae(θ)中,得到等效截收面积Ae。
由上述技术方案可知,本发明的有益效果:本发明提供的一种测算风力发电机组防雷等级的方法,包括以下步骤:S1、根据气象部门的数据计算风场所在地年平均雷击密度Ng;S2、根据回击电流强度对等效截面半径的影响,计算与风机接收相同雷击次数的等效截收面积Ae;S3、根据风场所在地形选择合适的环境因素Cd,并依据风力涡轮发电机系统防雷规范中的表述,计算风力发电机组的年平均直击雷数量Nd;S4、根据步骤S3得到的风力发电机组的年平均直击雷数量Nd,确定风力发电机组防雷等级。本发明提供的测算风力发电机组防雷等级的方法,便于得到更加科学合理的雷击等效截收面积,从而能够更加有效地为不同机型的风电机组选择合适的雷电防护等级,提高风电机组防雷保护设计施工过程中的经济型与有效性。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中,类似的元件或部分一般由类似的附图标记标识。附图中,各元件或部分并不一定按照实际的比例绘制。
图1为本发明的计算流程图;
图2为风力发电机组雷击区域的示意图;
图3为风力发电机组的结构示意图;
图4为风力发电机组等效截收面积的示意图。
附图标记:
1-杆塔、2-叶片。
具体实施方式
下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案,因此只作为示例,而不能以此来限制本发明的保护范围。
请参阅图1-4,本实施例提供的一种测算风力发电机组防雷等级的方法,包括以下步骤:
S1、根据气象部门的数据计算风场所在地的年平均雷击密度Ng,年平均雷击密度Ng可以通过闪电定位仪探测到的单位面积上的闪电总数计算得到。在缺乏闪电定位仪的地区,可以根据公式Ng=0.1×Td计算年平均雷击密度Ng;其中Ng的单位为km2year-1,Td表示年平均雷暴日,年平均雷暴日根据气象部门的统计数据计算。
S2、根据回击电流强度对等效截面半径的影响,计算与风机接收相同雷击次数的等效截收面积Ae,具体步骤如下:
S21、根据IEC62305规范计算闪击距离其中Ip是雷击电流幅值强度峰值;Ip可以通过气象部门的数据查得;
S22、根据闪击距离S和机组高度h,计算雷击区域在地面上投影的半径r与机组高度h之间的函数r(h);如图2所示,
当h≥βS时:r(h)=S ①
当h<βS时:
式②中β为防雷经验值;计算时β的取值范围为0.6-0.8。为了便于计算,β的近似值通常取为0.735。
S23、根据发电机组中塔杆的高度H、叶片的长度L以及旋转角度θ,计算出发电机组在旋转过程中的高度h(θ)和宽度W(θ),所述旋转角度θ为位于塔杆最顶端的叶片与竖直轴线之间的夹角;如图3所示,
h(θ)=H+Lcos(θ) ④
S24、将步骤S23中得到的h(θ)代入步骤S22所得到的函数r(h)中,计算雷击区域在地面上投影的半径r;
当I≤IB时,r=10×I0.65 ⑤
当I>IB时,
其中I表示雷击的电流值,IB是一个临界电流强度,与β的取值对应,表示当h=βS时的电流强度。因此,IB也可用θ表示:
由于因此/>
即
将式④代入式⑦中,得到
S25、如图4所示,根据规则建筑物的等效截收面积计算方法,计算等效截收面积与旋转角度θ之间的函数Ae(θ);即等效截收面积Ae(θ)等于以雷击区域在地面上投影的直径2r为长度,以W(θ)为宽度围成的矩形的面积加上雷击区域在地面上投影的面积。
Ae(θ)=πr2+2rW(θ) ⑨
并将步骤S24得到半径r和步骤S23得到的宽度W(θ)代入Ae(θ)中,得到等效截收面积Ae。
由于有些标准和文献资料中提到的等效截收面积与风力发电机高度3倍和平方之间存在正比关系,并未能包含任何物理含义而只是简单假设了一个较大的区域,使其能够笼统的涵盖任何情况下的等效截收面积。采用本文描述的方法,计算能够更加准确。
S3、根据风场所在地形选择合适的环境因素Cd,并依据风力涡轮发电机系统防雷规范中的表述,计算风力发电机组的年平均直击雷数量
Nd=Ng·Ae·Cd ⑩
将式⑨代入式⑩中,并对角度和电流大小做分段积分
其中Cd为环境因素,当风力发电设备位于平原地区时,Cd可取值为1,当风力发电设备位于山丘地区时,Cd可取值为2。
S4、根据步骤S3得到的风力发电机组的年平均直击雷数量,确定风力发电机组防雷等级。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (1)
1.一种测算风力发电机组防雷等级的方法,其特征在于,包括以下步骤:
S1、根据气象部门的数据计算风场所在地年平均雷击密度Ng;
S2、根据雷击电流强度对等效截面半径的影响,计算与风机接收相同雷击次数的等效截收面积Ae;具体步骤如下:
S21、根据IEC62305规范计算闪击距离其中Ip是雷击电流幅值强度峰值;
S22、根据闪击距离S和机组高度h,计算雷击区域在地面上投影的半径r与机组高度h之间的函数r(h);
当h≥βS时:r(h)=S;
当h<βS时:
式中β为防雷经验值;计算时β的取值范围为0.6-0.8;
S23、根据发电机组中塔杆的高度H、叶片的长度L以及旋转角度θ,计算出发电机组在旋转过程中的高度h(θ)和宽度W(θ),所述旋转角度θ为位于塔杆最顶端的叶片与竖直轴线之间的夹角;
h(θ)=H+Lcos(θ)
S24、将步骤S23中得到的h(θ)代入步骤S22所得到的函数r(h)中,计算雷击区域在地面上投影的半径r;
当I≤IB时,r=10×I0.65;
当I>IB时,
其中I表示雷击的电流值,IB是一个临界电流强度,与β的取值对应,表示当h=βS时的电流强度;因此,IB也可用θ表示:
由于因此/>
即则
S25、根据规则建筑物的等效截收面积计算方法,计算等效截收面积与旋转角度θ之间的函数Ae(θ);
Ae(θ)=πr2+2rW(θ)
并将步骤S24得到半径r和步骤S23得到的宽度W(θ)代入Ae(θ)中,得到等效截收面积Ae;
S3、根据风场所在地形选择合适的环境因素Cd,并依据风力涡轮发电机系统防雷规范中的表述,计算风力发电机组的年平均直击雷数量Nd;
Nd=Ng·Ae·Cd
将步骤S25得到的函数Ae(θ)代入上式,并对角度和电流大小做分段积分
其中Cd为环境因素,当风力发电设备位于平原地区时,Cd取值为1,当风力发电设备位于山丘地区时,Cd取值为2;
S4、根据步骤S3得到的风力发电机组的年平均直击雷数量Nd,确定风力发电机组防雷等级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811361780.7A CN109241693B (zh) | 2018-11-15 | 2018-11-15 | 测算风力发电机组防雷等级的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811361780.7A CN109241693B (zh) | 2018-11-15 | 2018-11-15 | 测算风力发电机组防雷等级的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109241693A CN109241693A (zh) | 2019-01-18 |
CN109241693B true CN109241693B (zh) | 2023-07-21 |
Family
ID=65074787
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811361780.7A Active CN109241693B (zh) | 2018-11-15 | 2018-11-15 | 测算风力发电机组防雷等级的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109241693B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111594398B (zh) * | 2019-02-20 | 2022-03-29 | 中车株洲电力机车研究所有限公司 | 一种风力发电机组主动智能防雷方法及系统 |
CN111752990B (zh) * | 2020-06-09 | 2024-02-02 | 贵州省气象灾害防御技术中心 | 一种雷电易发性等级划分方法 |
CN115276062B (zh) * | 2022-07-22 | 2024-04-26 | 国网江苏省电力有限公司苏州供电分公司 | 一种基于雷击概率的风储联合系统功率两层调控方法与系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175936A (zh) * | 2011-01-19 | 2011-09-07 | 广东电网公司电力科学研究院 | 一种配网避雷器给定置信水平下非限制性预期运行寿命评估方法 |
CN104133936A (zh) * | 2013-07-11 | 2014-11-05 | 高磊 | 考虑周边物体时计算建构筑物各部位雷击概率的方法 |
CN105760694A (zh) * | 2016-03-10 | 2016-07-13 | 吉林建筑大学 | 一种建筑物截收相同雷击次数的等效面积计算方法 |
JP2017082689A (ja) * | 2015-10-28 | 2017-05-18 | 正昭 茆原 | 風力発電設備への落雷の検出方法及び検出装置、並びに、風力発電設備への落雷状況の分析方法及び分析装置 |
-
2018
- 2018-11-15 CN CN201811361780.7A patent/CN109241693B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175936A (zh) * | 2011-01-19 | 2011-09-07 | 广东电网公司电力科学研究院 | 一种配网避雷器给定置信水平下非限制性预期运行寿命评估方法 |
CN104133936A (zh) * | 2013-07-11 | 2014-11-05 | 高磊 | 考虑周边物体时计算建构筑物各部位雷击概率的方法 |
JP2017082689A (ja) * | 2015-10-28 | 2017-05-18 | 正昭 茆原 | 風力発電設備への落雷の検出方法及び検出装置、並びに、風力発電設備への落雷状況の分析方法及び分析装置 |
CN105760694A (zh) * | 2016-03-10 | 2016-07-13 | 吉林建筑大学 | 一种建筑物截收相同雷击次数的等效面积计算方法 |
Non-Patent Citations (3)
Title |
---|
呼包鄂地区闪电定位资料分析及误差订正;王汉堃;《中国优秀硕士学位论文全文数据库(基础科学辑)》;20180715(第7期);第2章 * |
风力发电机雷电放电全过程电流测量系统设计;姒天军 等;《仪表技术与传感器》;20171031(第10期);全文 * |
风机雷电监测及雷击暂态特性研究;王瑞峰;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20180315(第3期);第5章 * |
Also Published As
Publication number | Publication date |
---|---|
CN109241693A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109241693B (zh) | 测算风力发电机组防雷等级的方法 | |
El‐Asha et al. | Quantification of power losses due to wind turbine wake interactions through SCADA, meteorological and wind LiDAR data | |
Boudia et al. | Investigation of wind power potential at Oran, northwest of Algeria | |
Zhang | Wind resource assessment and micro-siting: science and engineering | |
Khahro et al. | Evaluation of wind power production prospective and Weibull parameter estimation methods for Babaurband, Sindh Pakistan | |
Gualtieri et al. | Comparing methods to calculate atmospheric stability-dependent wind speed profiles: A case study on coastal location | |
Lee et al. | Assessment of offshore wind energy at Younggwang in Korea | |
WO2018019278A1 (zh) | 基于小型雷达的输电线路雷电监测预警方法及存储介质 | |
Fredianelli et al. | Analytical assessment of wind turbine noise impact at receiver by means of residual noise determination without the wind farm shutdown | |
Ahmad et al. | Field implementation and trial of coordinated control of wind farms | |
US20210312101A1 (en) | Method for determining a wind turbine layout | |
Martin et al. | Impact of near-future turbine technology on the wind power potential of low wind regions | |
Sanusi et al. | Estimation of wind energy potential in Southwestern Nigeria | |
Benmedjahed et al. | Wind Potential Assessment of Three Coastal Sites in Algeria: Calculation and Modeling of Wind Turbine Noise using Matlab | |
Long et al. | On the lightning incidence to wind farms | |
Kheder-Haddouche et al. | Feasibility study of a wind farm in el golea region in the Algerian sahara | |
WO2023035708A1 (zh) | 风机叶片声音信号采集最优位置及其选取方法 | |
Boudia et al. | On the use of wind energy at Tlemcen, North-western region of Algeria | |
US11243238B2 (en) | Method for warning about lightning activity in wind farms | |
Benmedjahed et al. | Temporal Assessment of Wind Energy Resource in Algerian Desert Sites: Calculation and Modelling of Wind Noise | |
Sharma et al. | Estimation of Weibull parameters and wind power density at different heights for Akal site at Jaisalmer in Rajasthan | |
Malcolm et al. | Estimation of the failure rate of wind turbine electrical systems exposed to lightning strikes | |
CN111709610A (zh) | 雷击风险的评估方法、装置、存储介质和电子设备 | |
Mohammed | Analysis of wind energy potential in Kano Nigeria using WindPRO: a case study of Kwankwasiyya city | |
Arroyo et al. | Comparative Analysis of the Weibull Model and Observed Wind Data in the City of Floridablanca, Colombia. |
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 |