CN104537186B - 一种模拟特高压直流输电线路雷击暂态过程的方法 - Google Patents
一种模拟特高压直流输电线路雷击暂态过程的方法 Download PDFInfo
- Publication number
- CN104537186B CN104537186B CN201510012592.3A CN201510012592A CN104537186B CN 104537186 B CN104537186 B CN 104537186B CN 201510012592 A CN201510012592 A CN 201510012592A CN 104537186 B CN104537186 B CN 104537186B
- Authority
- CN
- China
- Prior art keywords
- switch
- lightning
- thunderbolt
- transmission line
- current
- 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
Landscapes
- Suspension Of Electric Lines Or Cables (AREA)
Abstract
本发明公开了一种模拟特高压直流输电线路雷击暂态过程的方法。该方法仿真人员根据雷电录波、防雷准则和/或解析研究需要定给定雷电流参数,包括幅值I0,波头参数T1,波尾参数T2;根据式1选定所用雷电流数学模型i0(t)作为电流源的控制输入;信号处理控制模块运行控制第三开关S3数值,若S3为1,为造成闪络的绕击模式;若第三开关S3数值为0,未发生绝缘子闪络,开关始终为打开状态,此时为不造成闪络的绕击模式;第二开关S2先闭合,表示雷击作用于杆塔上;第一开关S1打开,则S1=0;信号处理控制模块运行控制第三开关S3值的变化,若第三开关S3值为1,发生绝缘子闪络,开关已经闭合,此时为反击模式。本发明在一种模式中可以实现多种雷击类型的仿真,应用范围广。
Description
技术领域
本发明涉及输电线路雷击暂态过程仿真,特别是涉及一种模拟特高压直流输电线路雷击暂态过程的方法。
背景技术
研究特高压直流输电线路雷击下暂态电气量的变化特征是实现雷击故障分析、耐雷水平评估等工作的重要前提。由于求解特高压直流输电线路雷击暂态电气量的数学模型复杂,故只能通过数值仿真的方法获得以上数据。目前用于特高压直流输电线路雷击暂态过程仿真的模型均采取用双指数函数模型表示的雷电流i0作为恒流源,再将雷电通道等效阻抗Z0与之并联后直接接于特高压直流输电线路上进行仿真。该方法对于其他雷电流数学模型的应用,各类雷击事件的模型整合尚未有提及。
当前关于特高压直流输电线路雷击暂态过程仿真的雷击模型主要实现途径是在EMTDC等电磁暂态仿真软件中已有仿特高压直流输电系统的情况下,建立用双指数函数描述的雷电流恒流源,根据彼德逊等效定则建立雷击等值电路,然后通过开关将雷击电路接入直流输电线路中。具体情况如下图1所示,其中虚线部分雷电流i0受双指数函数(式1)模型控制的恒流源,Z0为雷电通道波阻抗,根据国际标准取取300Ω。在仿真平台上进行雷击仿真时,直接通过对控制开关S1实行外部输入控制其闭合。图中的实线部分为基于电磁暂态仿真软件平台的特高于直流输电系统。
目前用于特高压直流输电线路雷击暂态过程仿真的雷击模型主要存在以下两方面的问题:
1)雷电流模拟方式单一,目前仿真只采用适合短波头波尾参数的双指数函数作为数学函数模型,对于适用于长波头波尾参数的其他雷电流数学模型并未提及。
2)在描述复杂多样的雷击过程时仅采用电流源通过开关直接作用于输电线路的方式进行处理,实际上只描述了非故障性绕击一种情况,对于反击以及绕击后闪络等情况并没有包含在内。
发明内容
本发明针对现有技术缺点,提供一种可根据雷电流情况自动模拟不造成闪络的绕击,造成闪络的绕击和反击多种情形的模拟特高压直流输电线路雷击暂态过程的方法,使参数与模型更匹配。
本发明基于电磁暂态仿真软件平台,提出一种用于特高压直流输电线路雷击暂态过程仿真的综合雷击模型。本发明可应用多种数学函数模型描述雷电流,并给定根据雷电波头波尾参数选择模型时的判据,而后提供包含三种典型雷击事件的综合雷击模型建立方法。
本发明目的通过如下技术方案实现:
一种模拟特高压直流输电线路雷击暂态过程的方法,包括如下步骤:
1)仿真人员根据雷电录波、防雷准则和/或解析研究需要定给定雷电流参数,包括幅值I0,波头参数T1,波尾参数T2;雷电通道波阻抗Z0取300Ω;根据式1选定所用雷电流数学模型i0(t)作为电流源的控制输入;
式1
T1<5μs时选取双指数函数模型;当T1>5μs,且T2<150μs时采用脉冲函数模型;当T1>5μs,且T2>150μs时采用海德勒函数模型;t表示时间;η为修正参数,由仿真人员确定;
2)在各监测点安装示波器记录雷击暂态过程中的电气暂态量,控制第一开关S1、第二开关S2、第三开关S3的输出情况进行仿真;第一开关S1、第二开关S2、第三开关S3用0和1数字量表示其开闭情况;开关打开时,用0表示;开关闭合时,用1表示;第一开关S1、第二开关S2分别代表雷击发生于线路上、雷击发生于杆塔上;第三开关S3初值为0;根据雷击点电位U0、杆塔电位Ut之间实时电压差的绝对值与绝缘子伏秒特性Uf是否有交点决定第三开关S3取值是否发生变化;绝缘子的伏秒特性Uf如式2,其中t为仿真时的当前时刻,该值通过电磁暂态仿真软件的时钟元件获取,tp表示仿真时雷击发生的时刻,在仿真进行前人为提前给定;则S3的取值如式3所示:
式2
式3
第一开关S1闭合,则S1=1,表示雷击作用于输电线路上;开关S2打开,则S2=0;信号处理控制模块运行控制第三开关S3数值,若S3为1,发生绝缘子闪络,为造成闪络的绕击模式;若第三开关S3数值为0,未发生绝缘子闪络,为不造成闪络的绕击模式;
第二开关S2先闭合,表示雷击作用于杆塔上;第一开关S1打开,则S1=0;信号处理控制模块运行控制第三开关S3值的变化,若第三开关S3值为1,发生绝缘子闪络,为反击模式。
相对于现有技术,本发明具有如下优点:
1)本发明通过输入参数控制第一开关S1、第二开关S2通断以及第三开关S3初始状态,再通过数学逻辑运算实现雷击点电位U0、杆塔电位Ut之间实时电压差的绝对值与绝缘子伏秒特性Uf的大小比较,并根据比较控制S3通断,进而确定雷击模型,在一种模式中可以实现多种雷击类型的仿真,应用范围广。
2)本发明用于特高压直流输电线路雷击暂态过程仿真的综合雷击模型是基于电磁暂态仿真平台提出的,按照三种典型模型的性质特点,将雷电流底层数学模型整合为一分段函数(式1),以波头波尾参数T1/T2的大小范围作为各函数段的划分条件。模型中点线部分表示利用电磁暂态软件实现仿真模型时,雷击点电位U0、杆塔电位Ut与信号处理控制模块之间的虚拟链接关系(参数传递关系),控制模块于电磁暂态仿真软件下根据数学运算与逻辑元件实现,方法简单实用。
附图说明
图1现有特高压直流输电线路雷击暂态过程仿真电路原理图;
图2特高压直流输电线路雷击暂态过程模拟的等效电气模型图;
图3特高压直流输电线路雷击暂态过程模拟模型图。
具体实施方式
为更好地理解本发明,下面结合附图对本发明作进一步的说明,但本发明的实施方式不限如此。
如图2、3所示,特高压直流输电线路雷击暂态过程模拟等效电气中其中虚线部分为综合雷击模型,包括雷电流输入端口、雷电通道波阻抗Z0、第一开关S1、第二开关S2、第三开关S3、基于电磁暂态仿真软件平台的特高于直流输电系统、信号处理控制模块和塔杆T;其中雷电流输入端口i0和雷电通道波阻抗Z0并联连接,一并联端接地,另一并联端分别与第一开关S1和第二开关S2的一端连接,第二开关S2的另一端与塔杆T的塔顶连接,两个第三开关S3的一端分别连接塔杆T的横担,另一端连接直流输电线路,与直流输电线路的连接点分别为C2和B2点;第三开关S3与塔杆T的横担的连接点分别为C1和B1点;信号处理控制模块分别与第一开关S1、第二开关S2、第三开关S3以及特高于直流输电系统连接。信号处理控制模块对第一开关S1、第二开关S2、第三开关S3实现断开闭合控制,以及用雷击点电位U0与杆塔电位Ut之间实时电压差的绝对值与绝缘子伏秒特性Uf比较,判断是否发生闪络;信号处理控制模块可选用EMTDC型信号处理控制模块,该模块现有技术应用于电磁暂态仿真软件平台。
雷击线路一般分为不造成闪络的绕击(非故障性绕击),造成闪络的绕击和反击三种典型情况。当雷电流直接作用于线路后,绝缘子两端由杆塔电位与线路之间造成的电压差过大时会造成闪络,此时发生故障性绕击,若不闪络则发生的是非故障性绕击。而反击则为雷电流先作用于杆塔再造成绝缘子发生闪络导致线路通过杆塔接地。本发明模拟特高压直流输电线路雷击暂态过程的方法可以不改变电路情况下自动模拟不造成闪络的绕击,造成闪络的绕击和反击三种典型情况。其中电流源i0(t)以及雷电通道波阻抗Z0构成彼德逊等效电路,表示雷电流底层模型。Z0表示雷电通道波阻抗,根据国际标准取300Ω,T表示用多波阻抗模型描述的杆塔。i0(t)由式1的雷电流模型中的数学函数决定:
式1
T1<5μs时选取双指数函数模型;当T1>5μs,且T2<150μs时采用脉冲函数模型;当T1>5μs,且T2>150μs时采用海德勒函数模型;
雷电流模型中的数学函数受控电流源作为雷电流发生器,是实现雷击仿真模型的前提;雷电流底层数学模型有以下三种:式1双指数函数模型、式2海德勒函数模型、式3脉冲函数模型:
(1)双指数函数模型
式4
(2)海德勒函数模型
式5
(3)脉冲函数模型
式6
其中T1表示波头参数,T2表示波尾参数,I0表示雷电流幅值大小,在海德勒函数模型和脉冲函数模型中还须考虑修正参数η;以上参数在进行仿真前,由仿真人员根据雷电录波、防雷准则或解析研究需要自行给定。
一种模拟特高压直流输电线路雷击暂态过程的方法,包括如下步骤:
1)仿真人员根据雷电录波、防雷准则和/或解析研究需要定给定雷电流参数,包括幅值I0,波头参数T1,波尾参数T2;根据式1选定所用雷电流数学模型i0(t)作为电流源的控制输入;
式1
T1<5μs时选取双指数函数模型;当T1>5μs,且T2<150μs时采用脉冲函数模型;当T1>5μs,且T2>150μs时采用海德勒函数模型;t表示时间;η为修正参数,由仿真人员确定;雷电通道波阻抗Z0取300Ω
2)在各监测点安装示波器记录雷击暂态过程中的电气暂态量,控制第一开关S1、第二开关S2、第三开关S3的输出情况进行仿真;第一开关S1、第二开关S2、第三开关S3用0和1数字量表示其开闭情况;开关打开时,用0表示;开关闭合时,用1表示;第一开关S1、第二开关S2分别代表雷击发生于线路上、雷击发生于杆塔上,故其初始值根据仿真人员选定雷击模型后确定,而状态在后续仿真中均不改变;第三开关S3初值为0;根据雷击点电位U0、杆塔电位Ut之间实时电压差的绝对值与绝缘子伏秒特性Uf是否有交点决定第三开关S3取值是否发生变化;绝缘子的伏秒特性Uf如式2,其中t为仿真时的当前时刻,该值通过电磁暂态仿真软件的时钟元件获取,tp表示仿真时雷击发生的时刻,在仿真进行前人为提前给定;则S3的取值如式3所示。
式2
式3
第一开关S1应先闭合,则S1=1,表示雷击作用于输电线路上;开关S2打开,则S2=0;信号处理控制模块运行控制第三开关S3数值,若S3为1,发生绝缘子闪络,开关已经闭合,为造成闪络的绕击模式;若第三开关S3数值为0,未发生绝缘子闪络,开关始终为打开状态,此时为不造成闪络的绕击模式;
第二开关S2应先闭合,表示雷击作用于杆塔上;第一开关S1打开,则S1=0;信号处理控制模块运行控制第三开关S3值的变化,若第三开关S3值为1,发生绝缘子闪络,开关已经闭合,此时为反击模式。
本发明通过输入参数控制第一开关S1、第二开关S2通断以及第三开关S3初始状态,再通过数学逻辑运算实现雷击点电位U0、杆塔电位Ut之间实时电压差的绝对值与绝缘子伏秒特性Uf的大小比较,并根据比较控制S3通断,进而确定雷击模型,在一种模式中可以实现多种雷击类型的仿真,应用范围广。
本发明用于特高压直流输电线路雷击暂态过程仿真的综合雷击模型是基于电磁暂态仿真平台提出的,按照三种典型模型的性质特点,将雷电流底层数学模型整合为一分段函数(式1),以波头波尾参数T1/T2的大小范围作为各函数段的划分条件。模型中点线部分表示利用电磁暂态软件实现仿真模型时,雷击点电位U0、杆塔电位Ut与信号处理控制模块之间的虚拟链接关系(参数传递关系),控制模块于电磁暂态仿真软件下根据数学运算与逻辑元件实现,方法简单实用。
Claims (2)
1.一种模拟特高压直流输电线路雷击暂态过程的方法,其特征在于包括如下步骤:
1)仿真人员根据雷电录波、防雷准则和/或解析研究需要定给定雷电流参数,包括幅值I0,波头参数T1,波尾参数T2;雷电通道波阻抗Z0取300Ω;根据式1选定所用雷电流数学模型i0(t)作为电流源的控制输入;
T1<5μs时选取双指数函数模型;当T1>5μs,且T2<150μs时采用脉冲函数模型;当T1>5μs,且T2>150μs时采用海德勒函数模型;t表示时间;η为修正参数,由仿真人员确定;
2)在各监测点安装示波器记录雷击暂态过程中的电气暂态量,控制第一开关S1、第二开关S2、第三开关S3的输出情况进行仿真;第一开关S1、第二开关S2、第三开关S3用0和1数字量表示其开闭情况;开关打开时,用0表示;开关闭合时,用1表示;第一开关S1、第二开关S2分别代表雷击发生于线路上、雷击发生于杆塔上;第三开关S3初值为0;根据雷击点电位U0、杆塔电位Ut之间实时电压差的绝对值与绝缘子伏秒特性Uf是否有交点决定第三开关S3取值是否发生变化;绝缘子的伏秒特性Uf如式2,其中t为仿真时的当前时刻,该值通过电磁暂态仿真软件的时钟元件获取,tp表示仿真时雷击发生的时刻,在仿真进行前人为提前给定;则S3的取值如式3所示:
第一开关S1闭合,则S1=1,表示雷击作用于输电线路上;开关S2打开,则S2=0;信号处理控制模块运行控制第三开关S3数值,若S3为1,发生绝缘子闪络,为造成闪络的绕击模式;若第三开关S3数值为0,未发生绝缘子闪络,为不造成闪络的绕击模式;
第二开关S2先闭合,表示雷击作用于杆塔上;第一开关S1打开,则S1=0;信号处理控制模块运行控制第三开关S3值的变化,若第三开关S3值为1,发生绝缘子闪络,为反击模式。
2.根据权利要求1所述的模拟特高压直流输电线路雷击暂态过程的方法,其特征在于,所述信号处理控制模块选用EMTDC型信号处理控制模块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510012592.3A CN104537186B (zh) | 2015-01-09 | 2015-01-09 | 一种模拟特高压直流输电线路雷击暂态过程的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510012592.3A CN104537186B (zh) | 2015-01-09 | 2015-01-09 | 一种模拟特高压直流输电线路雷击暂态过程的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104537186A CN104537186A (zh) | 2015-04-22 |
CN104537186B true CN104537186B (zh) | 2017-07-18 |
Family
ID=52852713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510012592.3A Active CN104537186B (zh) | 2015-01-09 | 2015-01-09 | 一种模拟特高压直流输电线路雷击暂态过程的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104537186B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241664B (zh) * | 2018-09-29 | 2021-05-25 | 贵州电网有限责任公司 | 一种关于风力发电机雷电电磁暂态特性的分析计算方法 |
CN110598348B (zh) * | 2019-09-23 | 2022-08-02 | 四川大学 | 一种基于Matlab的核信号产生及处理的仿真方法 |
CN112699576B (zh) * | 2021-01-27 | 2024-02-06 | 沈阳工业大学 | 一种特高压直流雷电电磁暂态工程分析方法 |
CN113191103A (zh) * | 2021-03-05 | 2021-07-30 | 南方电网科学研究院有限责任公司 | 一种计算配电线路雷电电磁暂态的方法和装置 |
CN113919183B (zh) * | 2021-12-13 | 2022-04-12 | 国网江西省电力有限公司电力科学研究院 | 基于rtds的直流控保系统雷击故障响应特性模拟方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530453A (zh) * | 2013-09-30 | 2014-01-22 | 南方电网科学研究院有限责任公司 | 雷击输电线路引起特高压直流系统发生闭锁的分析方法 |
CN103529332A (zh) * | 2013-10-28 | 2014-01-22 | 昆明理工大学 | 一种基于电压相关度和小波变换暂态能量分布特性的特高压直流输电线路雷击干扰识别方法 |
CN103646148A (zh) * | 2013-12-20 | 2014-03-19 | 国家电网公司 | 一种特高压输电线路雷电反击性能仿真计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6826730B2 (en) * | 1998-12-15 | 2004-11-30 | Texas Instruments Incorporated | System and method for controlling current in an integrated circuit |
-
2015
- 2015-01-09 CN CN201510012592.3A patent/CN104537186B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530453A (zh) * | 2013-09-30 | 2014-01-22 | 南方电网科学研究院有限责任公司 | 雷击输电线路引起特高压直流系统发生闭锁的分析方法 |
CN103529332A (zh) * | 2013-10-28 | 2014-01-22 | 昆明理工大学 | 一种基于电压相关度和小波变换暂态能量分布特性的特高压直流输电线路雷击干扰识别方法 |
CN103646148A (zh) * | 2013-12-20 | 2014-03-19 | 国家电网公司 | 一种特高压输电线路雷电反击性能仿真计算方法 |
Non-Patent Citations (1)
Title |
---|
结合矩量法和边界元法的综合接地算法研究;陈凡 等;《高电压技术》;20120331;第38卷(第3期);第729-736页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104537186A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104537186B (zh) | 一种模拟特高压直流输电线路雷击暂态过程的方法 | |
CN101425681B (zh) | 电网系统继电保护在线校核方法 | |
CN104734364B (zh) | 智能变电站继电保护装置sv、goose输入测试的方法及系统 | |
CN103728514B (zh) | 一种低压变频器故障诊断方法 | |
CN108415269B (zh) | 基于fpga的mmc系统的仿真平台 | |
CN104866665A (zh) | 基于接口等值与交互的含电力电子设备的混合仿真方法 | |
Jalili-Marandi et al. | A real-time dynamic simulation tool for transmission and distribution power systems | |
CN102831792B (zh) | 配网故障仿真系统及其仿真方法 | |
CN105045952A (zh) | 模块化多电平换流器多维度建模方法与仿真方法 | |
CN205608420U (zh) | 交直流并联电网线路继电保护实时仿真系统 | |
CN108647479A (zh) | 一种避雷器故障暂态波形诊断方法及装置 | |
CN104505941A (zh) | 智能变电站网络的监控系统及模拟运行系统 | |
CN103678811B (zh) | 电力牵引交流传动变流器中间直流回路故障建模仿真方法 | |
CN110399671A (zh) | 电网的继电保护装置的整定计算方法、装置及终端 | |
CN201876517U (zh) | 大气过压侵入波监测系统 | |
CN110414118A (zh) | 一种基于分离式建模的Boost变换器建模方法及应用 | |
CN111953657B (zh) | 一种序列-数据联合驱动的配电网cps网络攻击辨识方法 | |
CN104732453A (zh) | 基于互联互通影响因子的scd文件比对方法 | |
CN104410159A (zh) | 一种基于实时数字仿真的电网黑启动全过程校核方法 | |
CN111884190A (zh) | 一种限流电抗器的拓扑及绝缘设计方法及装置 | |
CN104820152B (zh) | 便携式直流开关级差仿真测试系统及其实现方法 | |
CN104699903B (zh) | 一种智能变电站二次信号传输系统评估装置及方法 | |
CN108169610A (zh) | 一种基于实时数字仿真器的配电网单相接地故障测试方法 | |
CN108647427B (zh) | 一种基于rtds主变压器空载冲击的仿真方法 | |
Kezunovic et al. | Optimized fault location |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |