CN114725911B - 一种风电场汇集站系统变压器中性点接地电阻的估算方法 - Google Patents
一种风电场汇集站系统变压器中性点接地电阻的估算方法 Download PDFInfo
- Publication number
- CN114725911B CN114725911B CN202210357880.2A CN202210357880A CN114725911B CN 114725911 B CN114725911 B CN 114725911B CN 202210357880 A CN202210357880 A CN 202210357880A CN 114725911 B CN114725911 B CN 114725911B
- Authority
- CN
- China
- Prior art keywords
- model
- transformer
- point
- phase
- wind turbine
- 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
- 230000007935 neutral effect Effects 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 40
- 230000001360 synchronised effect Effects 0.000 claims abstract description 27
- 238000013178 mathematical model Methods 0.000 claims abstract description 14
- 239000004020 conductor Substances 0.000 claims description 25
- 238000004088 simulation Methods 0.000 claims description 10
- 239000003990 capacitor Substances 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 5
- 238000004804 winding Methods 0.000 claims description 5
- 230000004907 flux Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000005294 ferromagnetic effect Effects 0.000 claims description 3
- 230000005291 magnetic effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
- H02H9/00—Emergency protective circuit arrangements for limiting excess current or voltage without disconnection
- H02H9/04—Emergency protective circuit arrangements for limiting excess current or voltage without disconnection responsive to excess voltage
- H02H9/041—Emergency protective circuit arrangements for limiting excess current or voltage without disconnection responsive to excess voltage using a short-circuiting device
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
- H02H9/00—Emergency protective circuit arrangements for limiting excess current or voltage without disconnection
- H02H9/02—Emergency protective circuit arrangements for limiting excess current or voltage without disconnection responsive to excess current
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/16—Estimation of constants, e.g. the rotor time constant
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/24—Vector control not involving the use of rotor position or rotor speed sensors
- H02P21/28—Stator flux based control
- H02P21/30—Direct torque control [DTC] or field acceleration method [FAM]
-
- 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/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种风电场汇集站系统变压器中性点接地电阻的估算方法,属于风电场汇集站过电压防护领域。所述方法首先根据风电场内的风力机组的一台风力机和一台永磁同步发电机建立起风力机组的数学模型;然后结合模型1‑5,建立起风电场汇集站等效电路模型。在系统变压器不同的接地电阻情况下,进行接地故障试验。将试验结果进行描点绘图,最终选取合适的接地电阻阻值。最后给出风电场汇集站35kV侧系统变压器接地电阻阻值选取的建议。
Description
技术领域
本发明属于风电场汇集站过电压防护领域,具体涉及一种风电场汇集站系统变压器中性点接地电阻的估算方法。
背景技术
随着风电场规模的不断扩大,累计装机容量增加以及电缆网的大量采用,接地电容电流成线性比例增长,风电场的安全稳定运行的问题日益突出。当风电场出现单相接地故障时,由于电容电流较大,接地点形成的故障电弧难以熄灭,间歇性电弧接地造成的过电压又将使设备受到的危害进一步加大。然而,中性点接地电阻阻值的不合理导致多台风机同时脱网,带来的电量损失日趋严重。
为尽可能避免风电场由于中性点接地电阻阻值的不合理带来的日趋严重的电量损失,有必要对变压器中性点接地电阻阻值的选择进行研究。目前,风电场35/220kV系统变压器的35kV侧中性点常用的经电阻接地方式分为三种:中性点经高电阻接地(接地电阻值>500Ω)、中性点经中电阻接地(10Ω<接地电阻值<500Ω)、中性点经小电阻接地(接地电阻值<10Ω)。中性点接地电阻阻值的选取问题的一直没有一个通用的方法,在风电场汇集站过电压防护领域还存在以下问题:
1)在风电场汇集线路连接多个风电机组,发生故障情况极为复杂。
2)汇集线路35kV/220系统变压器的35kV侧多采用中性点经电阻接地方式,但是接地电阻阻值选择比较困难。
因此,需要提出一种风电场汇集站35/220kV系统变压器的35kV侧中性点接地电阻的估算方法,为风电场汇集站过电压防护及评估提供科学分析手段。
发明内容
本发明要解决的技术问题是针对风电场汇集站连接的风电机组数量多,35kV侧系统变压器中性点接地电阻阻值选择不合理,导致风电场安全稳定运行问题。
本发明的目的是这样实现的,本发明提供了一种风电场汇集站系统变压器中性点接地电阻的估算方法,所述系统变压器为35/220kV系统变压器,所述风电场汇集站包括35/220kV系统变压器、系统集电线路、220kV的三相电源、母线、风电机组、升压变压器和风电集电线路;所述35/220kV系统变压器的35kV侧采用中性点经电阻接地方式;
所述估算方法,建立风电场汇集站集中参数等效电路模型,并通过该模型得到中性点接地电阻阻值的取值范围,具体的,包括以下步骤:
步骤1,对风电机组中的风力机和永磁同步发电机建立数学模型;
步骤1.1,使用简化的空气动力模型来描述风力机的机械功率P、风力机的风速γ和风力机的机械转矩T1的关系,并将该简化的空气动力模型记为模型1,其表达式如下:
T1=ρπγr3C/(2λ)
P=ωT1=ρπr2γCγ3/2
其中,
ρ为空气密度,r为风力机的叶轮半径,入为风力机的叶尖速比,ω为风力机的机械转速;
C为风力机的功率系数,式中,β为风力机的浆距角,λ0为浆距系数,/>
步骤1.2,忽略永磁同步发电机的铁磁饱和、涡流损耗和磁滞损耗,令永磁铁在定子绕组中建立的磁通为正弦分布,建立永磁同步发电机在d-q同步旋转坐标系下的数学模型,并记为模型2,其表达式如下:
其中,L为定子的电感,rs为定子的电阻,id为定子d轴的电流,iq为定子q轴的电流,ud为定子d轴的电压,uq为定子q轴的电压,ωe为坐标系的旋转角速度,ωg为永磁同步发电机的转子速度,n为转子极对数,Ψ1为永磁体产生的磁链;
根据模型2,永磁同步发电机的有功功率P1、永磁同步发电机的无功功率Q1、永磁同步发电机的电磁转矩T1的表达式如下:
步骤2,将35/220kV系统变压器和升压变压器统称为变压器,并设变压器中有N个支路,N为正整数,建立变压器的等值电路模型,记为模型3,模型3的列写支路电压方程如下:
其中,V1,V2,…,VN分别为N个支路的绕组电压;I1,I2,…,IN分别为N个支路的绕组电流;ч=1,2,…,N;6=1,2,…,N;在列写支路电压方程中,当ч、б相等时,Zчб是支路б的自阻抗,当ч、б不相等时,Zчб是支路б对支路ч的互阻抗;
步骤3,将系统集电线路和风电集电线路统称为集电线路,建立集电线路集中参数模型,并记为模型4,模型4包括集电线路的正序电抗X1、集电线路的正序电容C1和集电线路零序电容C0的表达式;
不考虑架空地线,按单根导线来计算,并令集电线路的零序电阻等于正序电阻,集电线路的正序电抗X1的表达式为:
其中,d1为三相导线间的平均几何距离,式中,dab为a相导线与b相导线间的距离,dbc为b相导线与c相导线间的距离,dca为c相导线与a相导线间的距离,a相导线、b相导线和c相导线半径相等均记为r1,f1为频率;
令集电线路的零序电抗为3.5倍的正序电抗,集电线路的正序电容C1、集电线路零序电容C0的表达式如下:
其中,d2为三相导线对其镜像的互几何均距,式中,H1是a相导线到镜像的距离,H2是b相导线到镜像的距离,H3是c相导线到镜像的距离;
步骤4,对220kV的三相电源采用电压源模型,并记为模型5,模型5的输出电压Us的表达式为:
Us=E-IsRs
其中,Rs为模型5的内部电阻,E为模型5的电动势,Is为模型5的输出电流值;
步骤5,根据步骤1,步骤2,步骤3,步骤4建立的数学模型,在仿真软件ATP Draw环境下,建立风电场汇集站集中参数等效电路模型;
设风电场汇集站模型包括一个220kV的三相电源、一条35kV母线、一台35/220kV系统变压器、两条相同的系统集电线路、M个相同的风电机组、M个与风电机组连接的升压变压器、M条与升压变压器连接的风电集电线路;记220kV的三相电源为三相电源S220、35kV母线为母线M35,35/220kV系统变压器为变压器B11,两条相同的系统集电线路分别记为集电线路LA和集电线路LB,M个相同的风电机组记为机组Aj,M个与风电机组连接的升压变压器记为变压器Bj,M条与升压变压器连接的风电集电线路记为集电线路Lj,j=1,2…M;
其中,M个机组Aj依次与变压器Bj、集电线路Lj串联,集电线路Lj接母线M35,三相电源S220依次与集电线路LA、变压器B11、集电线路LB串联,集电线路LB接M35;
步骤6,在步骤5得到的风电场汇集站集中参数等效电路模型中,改变变压器B11的35kV侧接地电阻阻值,进行单相接地故障下的仿真测试;
设定仿真测试中35kV侧接地电阻阻值按照每次增加10Ω等值变化,给定100个测试35kV侧接地电阻阻值,并记为测试电阻值Ri,i=1,2,…,100;
启动仿真测试,检测得到与测试电阻值Ri对应的故障点残流值Ii和中性点电压值Ui,i=1,2,…,100;
步骤7,将步骤6记录的数据以测试电阻值Ri为横坐标,故障点残流值Ii和中性点电压值Ui为纵坐标进行描点绘图,得到两条曲线,分别记为R-I曲线和R-U曲线,记两曲线的交点为点0;
做点0对横轴的平行线,交纵轴为点Y0,再以点Y0为起点,向上垂直平移0.1pu,在纵轴上得到点Y1,做点Y1与横轴的平行线J1,并将J1与R-I曲线相交点的测试电阻值记为RI1、J1与R-U曲线相交点的测试电阻值记为RU1,进行如下判断:
若RU1-Ri1≤1Ω,则中性点接地电阻R的取值范围为{R|Ri1≤R≤RU1,R为实数},估算结束;
若RU1-Ri1>1Ω,则进入步骤8;
所述pu为标幺值的单位;
步骤8,设在点Y0和点Y1之间在以下估算过程中共在纵轴上得到F-1个点,将该F-1个点中的任意一个点记为点YD,D为估算序号,D=2,3…F,点YD为点Y0向上平移0.1/2D-1pu;
按照D=2,3…F的顺序,依次做点YD与横轴的平行线JD,记录JD与R-I曲线相交点的测试电阻值RiD、JD与R-U曲线相交点的测试电阻值RUD,并进行以下判断;
若RUD-RiD≤1Ω,则中性点接地电阻R的取值范围为{R|RiD≤R≤RUD,R为实数},估算结束;
若RUD-RiD>1Ω,进入下一个点YD+1的判断,直到两条曲线上对应的两个测试电阻值之差≤1Ω,即得到中性点接地电阻R的取值范围。
与现有的估算方法相比,本发明的有益效果如下:
1、考虑到风电场35/220kV系统汇集线路风电机组数量多、易发生接地故障问题,选取了N个风电机组。基于风力机的详细模型,进行简化风力机和永磁同步电机发电机的模型,有效的减少了仿真规模和时间,但可以有效仿真汇集站故障情况下对风电场的影响。
2、中性点接地电阻阻值选取范围是0-1000Ω,涵盖了小电阻,中电阻和高电阻接地,给中性点接地电阻的选取提供有效建议。
3、将故障情况下的过电压和过电流试验数据绘制曲线,有效给出了接地电阻的取值范围。
附图说明
图1为本发明估算方法的流程图。
图2为本发明实施例中风电场汇集站集中参数等效电路模型示意图。
图3为本发明模型2中三相导线及其镜像的位置图。
图4为本发明步骤7-步骤8所述方法的示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的说明:
本发明为一种风电场汇集站系统变压器中性点接地电阻的估算方法,所述系统变压器为35/220kV系统变压器,所述风电场汇集站包括35/220kV系统变压器、系统集电线路、220kV的三相电源、母线、风电机组、升压变压器和风电集电线路。所述35/220kV系统变压器的35kV侧采用中性点经电阻接地方式。
图1是本发明估算方法的流程图。由该图可见,本发明的评估方法为建立风电场汇集站集中参数等效电路模型,并通过该模型得到中性点接地电阻阻值的取值范围,具体的,包括以下步骤:
步骤1,对风电机组中的风力机和永磁同步发电机建立数学模型;
步骤1.1,使用简化的空气动力模型来描述风力机的机械功率P、风力机的风速γ和风力机的机械转矩T1的关系,并将该简化的空气动力模型记为模型1,其表达式如下:
T1=ρπγr3C/(2λ)
P=ωT1=ρπr2γCγ3/2
其中,
ρ为空气密度,r为风力机的叶轮半径,λ为风力机的叶尖速比,ω为风力机的机械转速;
C为风力机的功率系数,式中,β为风力机的浆距角,λ0为浆距系数,/>
步骤1.2,忽略永磁同步发电机的铁磁饱和、涡流损耗和磁滞损耗,令永磁铁在定子绕组中建立的磁通为正弦分布,建立永磁同步发电机在d-q同步旋转坐标系下的数学模型,并记为模型2,其表达式如下:
其中,L为定子的电感,rs为定子的电阻,id为定子d轴的电流,iq为定子q轴的电流,ud为定子d轴的电压,uq为定子q轴的电压,ωe为坐标系的旋转角速度,ωg为永磁同步发电机的转子速度,n为转子极对数,Ψ1为永磁体产生的磁链。
根据模型2,永磁同步发电机的有功功率P1、永磁同步发电机的无功功率Q1、永磁同步发电机的电磁转矩T1的表达式如下:
步骤2,将35/220kV系统变压器和升压变压器统称为变压器,并设变压器中有N个支路,N为正整数,建立变压器的等值电路模型,记为模型3,模型3的列写支路电压方程如下:
其中,V1,V2,…,VN分别为N个支路的绕组电压;I1,I2,…,IN分别为N个支路的绕组电流;ч=1,2,…,N;6=1,2,…,N;在列写支路电压方程中,当ч、б相等时,Zчб是支路6的自阻抗,当ч、б不相等时,Zч6是支路6对支路ч的互阻抗。
步骤3,将系统集电线路和风电集电线路统称为集电线路,建立集电线路集中参数模型,并记为模型4,模型4包括集电线路的正序电抗X1、集电线路的正序电容C1和集电线路零序电容C0的表达式;
不考虑架空地线,按单根导线来计算,并令集电线路的零序电阻等于正序电阻,集电线路的正序电抗X1的表达式为:
其中,d1为三相导线间的平均几何距离,式中,dab为a相导线与b相导线间的距离,dbc为b相导线与c相导线间的距离,dca为c相导线与a相导线间的距离,a相导线、b相导线和c相导线半径相等均记为r1,f1为频率;
令集电线路的零序电抗为3.5倍的正序电抗,集电线路的正序电容C1、集电线路零序电容C0的表达式如下:
其中,d2为三相导线对其镜像的互几何均距,式中,H1是a相导线到镜像的距离,H2是b相导线到镜像的距离,H3是c相导线到镜像的距离。
图3给出了三相导线及其镜像的位置图。
步骤4,对220kV的三相电源采用电压源模型,并记为模型5,模型5的输出电压Us的表达式为:
Us=E-IsRs
其中,Rs为模型5的内部电阻,E为模型5的电动势,Is为模型5的输出电流值;
步骤5,根据步骤1,步骤2,步骤3,步骤4建立的数学模型,在仿真软件ATP Draw环境下,建立风电场汇集站集中参数等效电路模型。
图2给出了本实施例中风电场汇集站集中参数等效电路模型示意图。
设风电场汇集站模型包括一个220kV的三相电源、一条35kV母线、一台35/220kV系统变压器、两条相同的系统集电线路、M个相同的风电机组、M个与风电机组连接的升压变压器、M条与升压变压器连接的风电集电线路;记220kV的三相电源为三相电源S220、35kV母线为母线M35,35/220kV系统变压器为变压器B11,两条相同的系统集电线路分别记为集电线路LA和集电线路LB,M个相同的风电机组记为机组Aj,M个与风电机组连接的升压变压器记为变压器Bj,M条与升压变压器连接的风电集电线路记为集电线路Lj,j=1,2…M。
其中,M个机组Aj依次与变压器Bj、集电线路Lj串联,集电线路Lj接母线M35,三相电源S220依次与集电线路LA、变压器B11、集电线路LB串联,集电线路LB接M35。
步骤6,在步骤5得到的风电场汇集站集中参数等效电路模型中,改变变压器B11的35kV侧接地电阻阻值,进行单相接地故障下的仿真测试;
设定仿真测试中35kV侧接地电阻阻值按照每次增加10Ω等值变化,给定100个测试35kV侧接地电阻阻值,并记为测试电阻值Ri,i=1,2,…,100;
启动仿真测试,检测得到与测试电阻值Ri对应的故障点残流值Ii和中性点电压值Ui,i=1,2,…,100;
步骤7,将步骤6记录的数据以测试电阻值Ri为横坐标,故障点残流值Ii和中性点电压值Ui为纵坐标进行描点绘图,得到两条曲线,分别记为R-I曲线和R-U曲线,记两曲线的交点为点0;
做点0对横轴的平行线,交纵轴为点Y0,再以点Y0为起点,向上垂直平移0.1pu,在纵轴上得到点Y1,做点Y1与横轴的平行线J1,并将J1与R-I曲线相交点的测试电阻值记为RI1、J1与R-U曲线相交点的测试电阻值记为RU1,进行如下判断:
若RU1-RI1≤1Ω,则中性点接地电阻R的取值范围为{R|RI1≤R≤RU1,R为实数},估算结束;
若RU1-RI1>1Ω,则进入步骤8;
所述pu为标幺值的单位;
步骤8,设在点Y0和点Y1之间在以下估算过程中共在纵轴上得到F-1个点,将该F-1个点中的任意一个点记为点YD,D为估算序号,D=2,3…F,点YD为点Y0向上平移0.1/2D-1pu;
按照D=2,3…F的顺序,依次做点YD与横轴的平行线JD,记录JD与R-I曲线相交点的测试电阻值RID、JD与R-U曲线相交点的测试电阻值RUD,并进行以下判断;
若RUD-RID≤1Ω,则中性点接地电阻R的取值范围为{R|RID≤R≤RUD,R为实数},估算结束;
若RUD-RID>1Ω,进入下一个点YD+1的判断,直到两条曲线上对应的两个测试电阻值之差≤1Ω,即得到中性点接地电阻R的取值范围。
图4为本发明步骤7-步骤8所述方法的示意图。由该图可见,在点Y1和点Y0作为上下边界的区间,R-I曲线和R-U曲线呈V字状,即在该区间内一定存在两个对应测试电阻值之差≤1Ω的区间,利用该估算方法能够得到中性点接地电阻R的取值范围。另外调整Y0和Y1的平移距离、调整点YD向上平移距离及调整对应测试电阻值之差可以进行各种电阻要求下中性点接地电阻值的估算。
由以上可见,本发明对风力机组中的风力机和永磁同步发电机建立数学模型,然后结合模型2-模型5,建立起风电场汇集站等效电路模型。在系统变压器不同的接地电阻情况下,进行接地故障试验。将试验结果进行描点绘图,最终给出风电场汇集站35/220kV系统变压器35kV侧接地电阻阻值选取的建议。
Claims (1)
1.一种风电场汇集站系统变压器中性点接地电阻的估算方法,所述系统变压器为35/220kV系统变压器,所述风电场汇集站包括35/220kV系统变压器、系统集电线路、220kV的三相电源、母线、风电机组、升压变压器和风电集电线路;所述35/220kV系统变压器的35kV侧采用中性点经电阻接地方式;
所述估算方法,其特征在于,建立风电场汇集站集中参数等效电路模型,并通过该模型得到中性点接地电阻阻值的取值范围,具体的,包括以下步骤:
步骤1,对风电机组中的风力机和永磁同步发电机建立数学模型;
步骤1.1,使用简化的空气动力模型来描述风力机的机械功率P、风力机的风速Υ和风力机的机械转矩T1的关系,并将该简化的空气动力模型记为模型1,其表达式如下:
T1=ρπΥr3C/(2λ)
P=ωT1=ρπr2ΥCΥ3/2
其中,
ρ为空气密度,r为风力机的叶轮半径,λ为风力机的叶尖速比,ω为风力机的机械转速;
C为风力机的功率系数,式中,β为风力机的浆距角,λ0为浆距系数,/>
步骤1.2,忽略永磁同步发电机的铁磁饱和、涡流损耗和磁滞损耗,令永磁铁在定子绕组中建立的磁通为正弦分布,建立永磁同步发电机在d-q同步旋转坐标系下的数学模型,并记为模型2,其表达式如下:
其中,L为定子的电感,rs为定子的电阻,id为定子d轴的电流,iq为定子q轴的电流,ud为定子d轴的电压,uq为定子q轴的电压,ωe为坐标系的旋转角速度,ωg为永磁同步发电机的转子速度,n为转子极对数,Ψ1为永磁体产生的磁链;
根据模型2,永磁同步发电机的有功功率P1、永磁同步发电机的无功功率Q1、永磁同步发电机的电磁转矩T1的表达式如下:
步骤2,将35/220kV系统变压器和升压变压器统称为变压器,并设变压器中有N个支路,N为正整数,建立变压器的等值电路模型,记为模型3,模型3的列写支路电压方程如下:
其中,V1,V2,...,VN分别为N个支路的绕组电压;I1,I2,...,IN分别为N个支路的绕组电流;ч=1,2,…,在列写支路电压方程中,当ч、/>相等时,/>是支路/>的自阻抗,当ч、/>不相等时,/>是支路/>对支路ч的互阻抗;
步骤3,将系统集电线路和风电集电线路统称为集电线路,建立集电线路集中参数模型,并记为模型4,模型4包括集电线路的正序电抗X1、集电线路的正序电容C1和集电线路零序电容C0的表达式;
不考虑架空地线,按单根导线来计算,并令集电线路的零序电阻等于正序电阻,集电线路的正序电抗X1的表达式为:
其中,d1为三相导线间的平均几何距离,
式中,dab为a相导线与b相导线间的距离,dbc为b相导线与c相导线间的距离,dca为c相导线与a相导线间的距离,a相导线、b相导线和c相导线半径相等均记为r1,f1为频率;
令集电线路的零序电抗为3.5倍的正序电抗,集电线路的正序电容C1、集电线路零序电容C0的表达式如下:
其中,d2为三相导线对其镜像的互几何均距,式中,H1是a相导线到镜像的距离,H2是b相导线到镜像的距离,H3是c相导线到镜像的距离;
步骤4,对220kV的三相电源采用电压源模型,并记为模型5,模型5的输出电压Us的表达式为:
Us=E-IsRs
其中,Rs为模型5的内部电阻,E为模型5的电动势,Is为模型5的输出电流值;
步骤5,根据步骤1,步骤2,步骤3,步骤4建立的数学模型,在仿真软件ATP Draw环境下,建立风电场汇集站集中参数等效电路模型;
设风电场汇集站模型包括一个220kV的三相电源、一条35kV母线、一台35/220kV系统变压器、两条相同的系统集电线路、M个相同的风电机组、M个与风电机组连接的升压变压器、M条与升压变压器连接的风电集电线路;记220kV的三相电源为三相电源S220、35kV母线为母线M35,35/220kV系统变压器为变压器B11,两条相同的系统集电线路分别记为集电线路LA和集电线路LB,M个相同的风电机组记为机组Aj,M个与风电机组连接的升压变压器记为变压器Bj,M条与升压变压器连接的风电集电线路记为集电线路Lj,j=1,2...M;
其中,M个机组Aj依次与变压器Bj、集电线路Lj串联,集电线路Lj接母线M35,三相电源S220依次与集电线路LA、变压器B11、集电线路LB串联,集电线路LB接M35;
步骤6,在步骤5得到的风电场汇集站集中参数等效电路模型中,改变变压器B11的35kV侧接地电阻阻值,进行单相接地故障下的仿真测试;
设定仿真测试中35kV侧接地电阻阻值按照每次增加10Ω等值变化,给定100个测试35kV侧接地电阻阻值,并记为测试电阻值Ri,i=1,2,...,100;
启动仿真测试,检测得到与测试电阻值Ri对应的故障点残流值Ii和中性点电压值Ui,i=1,2,...,100;
步骤7,将步骤6记录的数据以测试电阻值Ri为横坐标,故障点残流值Ii和中性点电压值Ui为纵坐标进行描点绘图,得到两条曲线,分别记为R-I曲线和R-U曲线,记两曲线的交点为点0;
做点0对横轴的平行线,交纵轴为点Y0,再以点Y0为起点,向上垂直平移0.1pu,在纵轴上得到点Y1,做点Y1与横轴的平行线J1,并将J1与R-I曲线相交点的测试电阻值记为RI1、J1与R-U曲线相交点的测试电阻值记为RU1,进行如下判断:
若RU1-RI1≤1Ω,则中性点接地电阻R的取值范围为{R|RI1≤R≤RU1,R为实数},估算结束;
若RU1-RI1>1Ω,则进入步骤8;
所述pu为标幺值的单位;
步骤8,设在点Y0和点Y1之间在以下估算过程中共在纵轴上得到F-1个点,将该F-1个点中的任意一个点记为点YD,D为估算序号,D=2,3...F,点YD为点Y0向上平移0.1/2D-1pu;
按照D=2,3...F的顺序,依次做点YD与横轴的平行线JD,记录JD与R-I曲线相交点的测试电阻值RID、JD与R-U曲线相交点的测试电阻值RUD,并进行以下判断;
若RUD-RID≤1Ω,则中性点接地电阻R的取值范围为{R|RID≤R≤RUD,R为实数},估算结束;
若RUD-RID>1Ω,进入下一个点YD+1的判断,直到两条曲线上对应的两个测试电阻值之差≤1Ω,即得到中性点接地电阻R的取值范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210357880.2A CN114725911B (zh) | 2022-04-06 | 2022-04-06 | 一种风电场汇集站系统变压器中性点接地电阻的估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210357880.2A CN114725911B (zh) | 2022-04-06 | 2022-04-06 | 一种风电场汇集站系统变压器中性点接地电阻的估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114725911A CN114725911A (zh) | 2022-07-08 |
CN114725911B true CN114725911B (zh) | 2024-03-26 |
Family
ID=82240925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210357880.2A Active CN114725911B (zh) | 2022-04-06 | 2022-04-06 | 一种风电场汇集站系统变压器中性点接地电阻的估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114725911B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243230A (zh) * | 2015-11-02 | 2016-01-13 | 合肥工业大学 | 一种计及卸荷电路的直驱风电场的等值建模方法 |
CN105576654A (zh) * | 2016-01-08 | 2016-05-11 | 南方电网科学研究院有限责任公司 | 直驱风电场等值方法和系统 |
CN107272645A (zh) * | 2017-06-28 | 2017-10-20 | 中国能源建设集团甘肃省电力设计院有限公司 | 中性点经电阻接地的光伏电站并网故障模型与分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107944142B (zh) * | 2017-11-24 | 2019-12-03 | 国电联合动力技术有限公司 | 高电压穿越能力仿真评估模型及基于其的仿真评估方法 |
-
2022
- 2022-04-06 CN CN202210357880.2A patent/CN114725911B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243230A (zh) * | 2015-11-02 | 2016-01-13 | 合肥工业大学 | 一种计及卸荷电路的直驱风电场的等值建模方法 |
CN105576654A (zh) * | 2016-01-08 | 2016-05-11 | 南方电网科学研究院有限责任公司 | 直驱风电场等值方法和系统 |
CN107272645A (zh) * | 2017-06-28 | 2017-10-20 | 中国能源建设集团甘肃省电力设计院有限公司 | 中性点经电阻接地的光伏电站并网故障模型与分析方法 |
Non-Patent Citations (2)
Title |
---|
大规模风电集中接入对系统短路电流的影响分析;周专;吕盼;宋新甫;张增强;刘建亮;;电工电气;20151115(第11期);全文 * |
规模化风电接入对系统短路电流的影响分析;周专;赵志强;王新刚;张增强;宋新甫;关洪浩;;四川电力技术;20151020(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114725911A (zh) | 2022-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108490304B (zh) | 发电机定子绕组单相接地故障位置定位方法 | |
CN105375532B (zh) | 一种适用于双馈风电场联络线的故障选相方法 | |
Chandrasekar et al. | Dynamic phasor modeling of type 3 DFIG wind generators (including SSCI phenomenon) for short-circuit calculations | |
CN105470950B (zh) | 故障分析用永磁直驱风电场电磁暂态等值模型建立方法 | |
CN112615393B (zh) | 基于矢量拟合的直驱风电机组控制器参数辨识方法及装置 | |
CN110880743B (zh) | 一种基于肯德尔秩相关性的风电场站送出线路纵联保护方法 | |
CN104215904A (zh) | 一种全功率变流器风电机组低电压穿越测试系统及方法 | |
CN104820741B (zh) | 兼顾风场分散性与机组差异性的风电场动态等值方法 | |
CN112907075A (zh) | 一种电力系统综合负荷模型参数辨识方法 | |
CN107167703A (zh) | 一种风电场集电线路单相接地故障测距方法 | |
CN110165644B (zh) | 一种基于暂态电流时频特性的新能源场站纵联保护方法 | |
CN108614189B (zh) | 双馈风电场送出线三相短路故障方向判断方法及方向元件 | |
Elizondo et al. | Model reduction, validation, and calibration of wind power plants for dynamic studies | |
CN106897514B (zh) | 一种全功率变换型新能源场站的短路电流计算模型的建立方法 | |
CN114725911B (zh) | 一种风电场汇集站系统变压器中性点接地电阻的估算方法 | |
Patel et al. | On-line load test for induction machine stator inter-turn fault detection under stator electrical asymmetries | |
Howard et al. | Calculation of fault current contribution of Type I wind turbine-generators | |
CN112649756A (zh) | 风电场集电线单相接地故障测距方法、系统、介质及设备 | |
Su et al. | Power quality assessment of wind turbines by Matlab/Simulink | |
CN116187082A (zh) | 风电场单机等值建模方法 | |
Qiao et al. | Faulty Phase Selection and Location Method of Rotor Ground Fault for Variable Speed Pumped Storage Generator-Motor | |
CN105550496B (zh) | 一种双馈风力发电机定子绕组匝间短路数学建模方法 | |
Kocar et al. | Nonlinear network equivalents of systems with inverter based resources to study unbalanced faults in steady state | |
CN113341223A (zh) | 一种适用于电网谐波分析及谐波源定位的方法 | |
Fang et al. | Modeling of Doubly-Fed Induction Generator with strong slot harmonic |
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 |