CN114725911A - 一种风电场汇集站系统变压器中性点接地电阻的估算方法 - Google Patents

一种风电场汇集站系统变压器中性点接地电阻的估算方法 Download PDF

Info

Publication number
CN114725911A
CN114725911A CN202210357880.2A CN202210357880A CN114725911A CN 114725911 A CN114725911 A CN 114725911A CN 202210357880 A CN202210357880 A CN 202210357880A CN 114725911 A CN114725911 A CN 114725911A
Authority
CN
China
Prior art keywords
model
transformer
wind power
point
phase
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.)
Granted
Application number
CN202210357880.2A
Other languages
English (en)
Other versions
CN114725911B (zh
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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202210357880.2A priority Critical patent/CN114725911B/zh
Publication of CN114725911A publication Critical patent/CN114725911A/zh
Application granted granted Critical
Publication of CN114725911B publication Critical patent/CN114725911B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H9/00Emergency protective circuit arrangements for limiting excess current or voltage without disconnection
    • H02H9/04Emergency protective circuit arrangements for limiting excess current or voltage without disconnection responsive to excess voltage
    • H02H9/041Emergency protective circuit arrangements for limiting excess current or voltage without disconnection responsive to excess voltage using a short-circuiting device
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H9/00Emergency protective circuit arrangements for limiting excess current or voltage without disconnection
    • H02H9/02Emergency protective circuit arrangements for limiting excess current or voltage without disconnection responsive to excess current
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/16Estimation of constants, e.g. the rotor time constant
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24Vector control not involving the use of rotor position or rotor speed sensors
    • H02P21/28Stator flux based control
    • H02P21/30Direct torque control [DTC] or field acceleration method [FAM]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power 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为风力机的功率系数,
Figure BDA0003582182050000021
式中,β为风力机的浆距角,λ0为浆距系数,
Figure BDA0003582182050000022
步骤1.2,忽略永磁同步发电机的铁磁饱和、涡流损耗和磁滞损耗,令永磁铁在定子绕组中建立的磁通为正弦分布,建立永磁同步发电机在d-q同步旋转坐标系下的数学模型,并记为模型2,其表达式如下:
Figure BDA0003582182050000023
Figure BDA0003582182050000031
其中,L为定子的电感,rs为定子的电阻,id为定子d轴的电流,iq为定子q轴的电流,ud为定子d轴的电压,uq为定子q轴的电压,ωe为坐标系的旋转角速度,ωg为永磁同步发电机的转子速度,n为转子极对数,Ψ1为永磁体产生的磁链;
根据模型2,永磁同步发电机的有功功率P1、永磁同步发电机的无功功率Q1、永磁同步发电机的电磁转矩T1的表达式如下:
Figure BDA0003582182050000032
Figure BDA0003582182050000033
Figure BDA0003582182050000034
步骤2,将35/220kV系统变压器和升压变压器统称为变压器,并设变压器中有N个支路,N为正整数,建立变压器的等值电路模型,记为模型3,模型3的列写支路电压方程如下:
Figure BDA0003582182050000035
其中,V1,V2,…,VN分别为N个支路的绕组电压;I1,I2,…,IN分别为N个支路的绕组电流;ч=1,2,…,N;6=1,2,…,N;在列写支路电压方程中,当ч、б相等时,Zчб是支路б的自阻抗,当ч、б不相等时,Zчб是支路б对支路ч的互阻抗;
步骤3,将系统集电线路和风电集电线路统称为集电线路,建立集电线路集中参数模型,并记为模型4,模型4包括集电线路的正序电抗X1、集电线路的正序电容C1和集电线路零序电容C0的表达式;
不考虑架空地线,按单根导线来计算,并令集电线路的零序电阻等于正序电阻,集电线路的正序电抗X1的表达式为:
Figure BDA0003582182050000041
其中,d1为三相导线间的平均几何距离,
Figure BDA0003582182050000042
式中,dab为a相导线与b相导线间的距离,dbc为b相导线与c相导线间的距离,dca为c相导线与a相导线间的距离,a相导线、b相导线和c相导线半径相等均记为r1,f1为频率;
令集电线路的零序电抗为3.5倍的正序电抗,集电线路的正序电容C1、集电线路零序电容C0的表达式如下:
Figure BDA0003582182050000043
Figure BDA0003582182050000044
其中,d2为三相导线对其镜像的互几何均距,
Figure BDA0003582182050000045
式中,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为风力机的功率系数,
Figure BDA0003582182050000071
式中,β为风力机的浆距角,λ0为浆距系数,
Figure BDA0003582182050000072
步骤1.2,忽略永磁同步发电机的铁磁饱和、涡流损耗和磁滞损耗,令永磁铁在定子绕组中建立的磁通为正弦分布,建立永磁同步发电机在d-q同步旋转坐标系下的数学模型,并记为模型2,其表达式如下:
Figure BDA0003582182050000073
Figure BDA0003582182050000074
其中,L为定子的电感,rs为定子的电阻,id为定子d轴的电流,iq为定子q轴的电流,ud为定子d轴的电压,uq为定子q轴的电压,ωe为坐标系的旋转角速度,ωg为永磁同步发电机的转子速度,n为转子极对数,Ψ1为永磁体产生的磁链。
根据模型2,永磁同步发电机的有功功率P1、永磁同步发电机的无功功率Q1、永磁同步发电机的电磁转矩T1的表达式如下:
Figure BDA0003582182050000075
Figure BDA0003582182050000081
Figure BDA0003582182050000082
步骤2,将35/220kV系统变压器和升压变压器统称为变压器,并设变压器中有N个支路,N为正整数,建立变压器的等值电路模型,记为模型3,模型3的列写支路电压方程如下:
Figure BDA0003582182050000083
其中,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的表达式为:
Figure BDA0003582182050000084
其中,d1为三相导线间的平均几何距离,
Figure BDA0003582182050000085
式中,dab为a相导线与b相导线间的距离,dbc为b相导线与c相导线间的距离,dca为c相导线与a相导线间的距离,a相导线、b相导线和c相导线半径相等均记为r1,f1为频率;
令集电线路的零序电抗为3.5倍的正序电抗,集电线路的正序电容C1、集电线路零序电容C0的表达式如下:
Figure BDA0003582182050000091
Figure BDA0003582182050000092
其中,d2为三相导线对其镜像的互几何均距,
Figure BDA0003582182050000093
式中,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为风力机的功率系数,
Figure FDA0003582182040000011
式中,β为风力机的浆距角,λ0为浆距系数,
Figure FDA0003582182040000012
步骤1.2,忽略永磁同步发电机的铁磁饱和、涡流损耗和磁滞损耗,令永磁铁在定子绕组中建立的磁通为正弦分布,建立永磁同步发电机在d-q同步旋转坐标系下的数学模型,并记为模型2,其表达式如下:
Figure FDA0003582182040000013
Figure FDA0003582182040000014
其中,L为定子的电感,rs为定子的电阻,id为定子d轴的电流,iq为定子q轴的电流,ud为定子d轴的电压,uq为定子q轴的电压,ωe为坐标系的旋转角速度,ωg为永磁同步发电机的转子速度,n为转子极对数,Ψ1为永磁体产生的磁链;
根据模型2,永磁同步发电机的有功功率P1、永磁同步发电机的无功功率Q1、永磁同步发电机的电磁转矩T1的表达式如下:
Figure FDA0003582182040000021
Figure FDA0003582182040000022
Figure FDA0003582182040000023
步骤2,将35/220kV系统变压器和升压变压器统称为变压器,并设变压器中有N个支路,N为正整数,建立变压器的等值电路模型,记为模型3,模型3的列写支路电压方程如下:
Figure FDA0003582182040000024
其中,V1,V2,...,VN分别为N个支路的绕组电压;I1,I2,...,IN分别为N个支路的绕组电流;ч=1,2,…,
Figure FDA0003582182040000027
在列写支路电压方程中,当ч、
Figure FDA00035821820400000212
相等时,
Figure FDA0003582182040000029
是支路
Figure FDA00035821820400000213
的自阻抗,当ч、
Figure FDA0003582182040000028
不相等时,
Figure FDA00035821820400000210
是支路
Figure FDA00035821820400000211
对支路ч的互阻抗;
步骤3,将系统集电线路和风电集电线路统称为集电线路,建立集电线路集中参数模型,并记为模型4,模型4包括集电线路的正序电抗X1、集电线路的正序电容C1和集电线路零序电容C0的表达式;
不考虑架空地线,按单根导线来计算,并令集电线路的零序电阻等于正序电阻,集电线路的正序电抗X1的表达式为:
Figure FDA0003582182040000025
其中,d1为三相导线间的平均几何距离,
Figure FDA0003582182040000026
式中,dab为a相导线与b相导线间的距离,dbc为b相导线与c相导线间的距离,dca为c相导线与a相导线间的距离,a相导线、b相导线和c相导线半径相等均记为r1,f1为频率;
令集电线路的零序电抗为3.5倍的正序电抗,集电线路的正序电容C1、集电线路零序电容C0的表达式如下:
Figure FDA0003582182040000031
Figure FDA0003582182040000032
其中,d2为三相导线对其镜像的互几何均距,
Figure FDA0003582182040000033
式中,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的取值范围。
CN202210357880.2A 2022-04-06 2022-04-06 一种风电场汇集站系统变压器中性点接地电阻的估算方法 Active CN114725911B (zh)

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 true CN114725911A (zh) 2022-07-08
CN114725911B 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 (4)

* Cited by examiner, † Cited by third party
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 中国能源建设集团甘肃省电力设计院有限公司 中性点经电阻接地的光伏电站并网故障模型与分析方法
US20200082043A1 (en) * 2017-11-24 2020-03-12 Guodian United Power Technology Company Ltd. Simulation evaluation model of high voltage ride through capability, simulation evaluation method based on the same and storage medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
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 中国能源建设集团甘肃省电力设计院有限公司 中性点经电阻接地的光伏电站并网故障模型与分析方法
US20200082043A1 (en) * 2017-11-24 2020-03-12 Guodian United Power Technology Company Ltd. Simulation evaluation model of high voltage ride through capability, simulation evaluation method based on the same and storage medium

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周专;吕盼;宋新甫;张增强;刘建亮;: "大规模风电集中接入对系统短路电流的影响分析", 电工电气, no. 11, 15 November 2015 (2015-11-15) *
周专;赵志强;王新刚;张增强;宋新甫;关洪浩;: "规模化风电接入对系统短路电流的影响分析", 四川电力技术, no. 05, 20 October 2015 (2015-10-20) *

Also Published As

Publication number Publication date
CN114725911B (zh) 2024-03-26

Similar Documents

Publication Publication Date Title
CN108490304B (zh) 发电机定子绕组单相接地故障位置定位方法
Chandrasekar et al. Dynamic phasor modeling of type 3 DFIG wind generators (including SSCI phenomenon) for short-circuit calculations
US20160202308A1 (en) Method and system for identifying full parameters of element by fault recording and fault locating method
Schulz et al. Dynamic models of turbine generators derived from solid rotor equivalent circuits
CN104991161A (zh) 双馈风力发电机定子绕组匝间短路建模及特征提取方法
CN111308264B (zh) 一种基于余弦相似度的配电网单相接地故障区段定位方法
CN111123105B (zh) 基于高频信号注入的电机匝间短路故障诊断方法
CN113326673B (zh) 一种同步电机的vbr模型电磁暂态仿真方法及系统
CN104820741B (zh) 兼顾风场分散性与机组差异性的风电场动态等值方法
CN109239513B (zh) 一种发电机定子绕组单相接地故障的选相方法
CN107167703A (zh) 一种风电场集电线路单相接地故障测距方法
Fan et al. Effect of damper winding and stator slot skewing structure on no-load voltage waveform distortion and damper bar heat in large tubular hydro generator
CN112816831A (zh) 风电场集电线单相接地故障定位方法
Howard et al. Calculation of fault current contribution of Type I wind turbine-generators
CN114725911B (zh) 一种风电场汇集站系统变压器中性点接地电阻的估算方法
CN115184718B (zh) 一种基于时域分析的三芯电缆故障定位方法
CN107862167B (zh) 一种考虑瞬时铁耗的开关磁阻电机建模方法
CN113341223B (zh) 一种适用于电网谐波分析及谐波源定位的方法
CN111475929B (zh) 基于配电网实境试验平台监测数据的反演验证方法及系统
CN105550496B (zh) 一种双馈风力发电机定子绕组匝间短路数学建模方法
Rehman et al. Detection of rotor inter-turn short circuit fault in doubly-fed induction generator using FEM simulation
Elez et al. Detection of inter-coil short circuits in synchronous machine armature winding on the basis of analysis of machine magnetic field
Kocar et al. Nonlinear network equivalents of systems with inverter based resources to study unbalanced faults in steady state
Singh et al. Validation and analysis of wind power plant models using short-circuit field measurement data
CN109782172A (zh) 一种同步电机参数试验测量方法

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