一种换流站交流侧阻抗在线测试方法与流程

本发明涉及交流电网阻抗测试技术领域，具体涉及一种换流站交流侧阻抗在线测试方法。

背景技术：

高压直流输电(highvoltagedirectcurrenttransmission，hvdc)具有迅速灵活的功率调节能力和实现区域电网异步互联的功能，是新能源的集约化开发和资源优化配置的关键技术之一。随着电力系统“源-网-荷”各环节中电力电子变流器渗透率不断提高，尤其输配电线路较长、区域电网中隔离变压器较多时，在hvdc换流站交流侧的连接点，电网呈现以低短路比为显著特征的弱电网特性。换流站的变流器与弱电网的动态交互作用易引发谐振和宽频谐振，甚至在稳定裕度不足时导致变流器-弱电网互联系统发生振荡失稳，危害系统稳定运行。

近年来，国内外学者针对电力电子化系统的宽频带振荡问题展开了广泛的研究，阻抗法已成为分析该系统稳定性问题的有效工具。实时且准确的阻抗特性数据是应用阻抗法分析稳定性的基础，故阻抗在线测试是稳定性在线分析的关键环节，阻抗在线测试方法的研究具有重要意义。国内外已有阻抗测试方法主要面向于低电压、小容量新能源发电装备及微电网系统，而hvdc换流站具有高压大容量的特点，内部电路结构复杂多样，且存在耦合的多时间尺度控制机制，导致了外部阻抗特性的复杂化。已有研究缺乏对于hvdc换流站交流侧阻抗在线测试方法的讨论和设计，如何兼顾换流站阻抗测试的快速性和准确性是换流站-电网互联系统稳定性分析及应用中亟待解决的问题。

技术实现要素：

本发明提出的一种换流站交流侧阻抗在线测试方法，可解决现有阻抗测试方法复杂且准确率低的技术问题。

为实现上述目的，本发明采用了以下技术方案：

一种换流站交流侧阻抗在线测试方法，包括以下步骤：

s1、生成测试频率序列；

s2、在换流站交流侧串联第一组三相扰动电压；

s3、提取换流站的第一组三相响应电流；

s4、在换流站交流侧串联第二组三相扰动电压；

s5、提取换流站的第二组三相响应电流；

s6、计算换流站交流侧在测试频率点的阻抗值；

进一步的，

步骤s1、生成测试频率序列的具体实现过程包括：

s11-测试频率序列的初步确定：根据测试频率密度需求，令同步旋转坐标系下的测试频率值在水平对数坐标轴等间距分布；

s12-并网电流频谱分析：记录时长1秒的换流站入网电流波形，并变换至与电网电压同步的旋转坐标系下，通过离散傅里叶变换输出电流频谱。首先，记录频谱中直流分量幅值，即abc坐标系下并网电流基波分量幅值if0；其次，记录系统谐波谐振频率fr，即并网电流频谱中除去直流分量和基波频率整数倍分量后幅值尖峰所在频率。

s13-测试频率调整：将初步确定的测试频率值均调整为整数，若频率值为基波频率整数倍或系统谐波谐振频率，则将该频率值替换为相近整数，最终生成由n个测试点组成的测试频率序列(f＝[f1,f2…fm…fn])；

进一步的，

步骤s2、在换流站交流侧串联第一组三相扰动电压，该电压为等效d通道扰动电压，其形式为三个周期性二进制序列，具体实现过程包括：

s21-给定扰动量幅值，基于静止-旋转坐标系间转换关系，根据步骤s1所得的频率序列依次计算用于合成各频率d通道扰动umd(1≤m≤n)的两组三相量(um1a，um1b，um1c)、(um2a，um2b，um2c)；其中，扰动量幅值可设置为基波电压幅值的3％至10％。

s22-叠加合成三相扰动电压，具体表达式分别为：

s23-应用软件算法拟合，生成与三相扰动电压(uda，udb，udc)频谱相近的三个二进制序列(vda，vdb，vdc)，该序列中频率fm(1≤m≤n)扰动电压分量在同步旋转坐标系下的复数表示为(um1d，um1q)。

进一步的，

步骤s3、提取换流站的第一组三相响应电流的具体实现过程包括：

s31-提取响应电流的基波幅值if1，并计算基波电流幅值偏移率l1：

若扰动电压引入的基波电流幅值偏移率l1超过允许值α％，则说明选取的扰动量幅值过大，需将扰动量的幅值降低比例β％重新进行扰动注入与响应提取，α和β的具体值可结合实际系统需要灵活设置。

s32-依次提取频率fm(1≤m≤n)的响应电流分量，根据同步旋转坐标系下的谐波幅值和相位以复数形式记录为(im1d，im1q)；

进一步的，

步骤s4、在换流站交流侧串联第二组三相扰动电压，该电压为等效q通道扰动电压，形式为三个周期性二进制序列，具体实现过程包括：

s41-给定扰动量幅值，基于静止-旋转坐标系间转换关系，根据步骤s1所得的频率序列依次计算用于合成各频率q通道扰动umq(1≤m≤n)的两组三相量(um3a，um3b，um3c)、(um4a，um4b，um4c)；其中，扰动量幅值可设置为基波电压幅值的3％至10％。

s42-叠加合成三相扰动电压，具体表达式分别为：

s43-通过软件算法拟合，生成与三相扰动电压(uqa，uqb，uqc)频谱相近的三个二进制序列(vqa，vqb，vqc)，该序列中频率fm(1≤m≤n)扰动电压分量在同步旋转坐标系下的复数表示为(um2d，um2q)。

进一步的，

步骤s5、提取换流站的第二组三相响应电流具体实现过程包括：

s51-提取响应电流的基波幅值if2，并计算基波电流幅值偏移率l2：

若扰动电压引入的基波电流幅值偏移率l2超过允许值α％，则将扰动量的幅值降低比例β％重新进行扰动注入与响应提取。

s52-依次提取频率fm(1≤m≤n)的响应电流分量，根据同步旋转坐标系下的谐波幅值和相位以复数形式记录为(im2d，im2q)；

进一步的，

步骤s6、计算换流站交流侧在测试频率点阻抗值的具体实现过程中采用不同通道阻抗独立计算的形式；首先，d通道阻抗包括zdd(m)和zdq(m)，计算方法为联立以下两个等式：

d通道阻抗计算结果为：

q通道阻抗包括zqd(m)和zqq(m)，计算方法为联立以下两个等式：

q通道阻抗计算结果为：

由上述技术方案可知，本发明的换流站交流侧阻抗在线测试方法可以高效获取具有实时性的hvdc换流站交流侧阻抗特性数据，该数据可应用于换流站-电网互联系统的稳定性在线分析，为高压直流输电系统的运行提供指导。

附图说明

图1为本发明提出的阻抗测试方法流程示意图；

图2为本发明实施例特高压直流输电系统结构图；

图3为本发明实施例扰动合成示意图；

图4为本发明实施例静止-旋转坐标系频率转换关系示意图；

图5为本发明实施例扰动注入主电路示意图；

图6为本发明实施例换流站交流侧阻抗测试结果。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。

如图1所示，本实施例所述的换流站交流侧阻抗在线测试方法，主要包括以下步骤：s1、生成测试频率序列；s2、在换流站交流侧串联第一组三相扰动电压；s3、提取换流站的第一组三相响应电流；s4、在换流站交流侧串联第二组三相扰动电压；s5、提取换流站的第二组三相响应电流；s6、计算换流站交流侧在测试频率点的阻抗值；

以下结合具体应用具体说明：

昌吉-古泉±1100kv特高压直流输电工程系统结构如图2所示，该工程线路全长3284km，额定直流电流5455a，额定容量12000mw，采用双极输电的方式，整流站和逆变站中每极含有2个串联的12脉波换流器，每个12脉波换流器两端的直流电压为550kv。新疆昌吉换流站(整流站)接入750kv交流电网，安徽古泉换流站(逆变站)采用分层接入的方式，高端换流变接入500kv交流电网，低端换流变接入1000kv交流电网，换流站交流母线均配置交流滤波器组。本实施例根据该高压直流输电系统的工程参数，在matlab/simulink中搭建仿真模型，并应用提出的一种换流站交流侧阻抗测试方法在500kv交流电网的接入点进行逆变站交流侧阻抗测试。

步骤s1、生成测试频率序列的具体实现过程包括：

s11-测试频率序列的初步确定：本实施例将测试频率范围设置为1hz-2000hz，对数坐标系下每十倍频程等间距分布10个频率点。

s12-并网电流频谱分析：记录时长1秒的换流站入网电流波形，并变换至与电网电压同步的旋转坐标系下，通过离散傅里叶变换输出电流频谱。首先，记录频谱中直流分量幅值，即abc坐标系下并网电流基波分量幅值if0，该值为9240a。其次，记录谐波谐振频率fr，即并网电流频谱中除去直流分量和基波频率整数倍分量后幅值尖峰所在频率；当电网阻抗等效电感分量和等效电阻分量分别为5mh和0.06ω时，记录的系统谐波谐振频率fr＝12hz。

s13-测试频率调整：将初步确定的测试频率值均调整为整数，若频率值为基波整数倍频或系统谐波谐振频率，则将该频率值替换为相近的整数频率点，最终生成由33个测试点组成的测试频率序列(f＝[12345678910131620253240526379102126159198253316398502631794995126015851995])；

步骤s2、在换流站交流侧串联第一组三相扰动电压；该电压为等效d通道扰动电压，合成方法如图3第一个子图所示，三相谐波电压经坐标变换在同步旋转坐标系下为旋转矢量，若同时注入关于d轴对称的正向旋转矢量vp1和反向旋转矢量vn1，由于矢量旋转速度相同且方向相反，q轴分量保持为零，系统中仅存在d通道的扰动电压。

根据静止-旋转坐标系坐标变换公式，以电网基波频率f0＝50hz为例，频率大于50hz的正序谐波在旋转坐标系下为正向旋转矢量，频率小于50hz的正序谐波和所有频率的负序谐波在旋转坐标系下为反向旋转矢量，同步旋转坐标系下正负序谐波对应的旋转矢量与abc坐标系下正负序量的频率转换关系如图4所示。

扰动电压注入时形式为三个周期性二进制序列，具体实现过程包括：

s21-将扰动量幅值设置为交流侧电网电压的5％，结合图3所示的扰动合成方法和图4所示的静止-旋转坐标系频率转换关系，依次计算用于合成各频率d通道扰动umd(1≤m≤33)的两组三相量(um1a，um1b，um1c)、(um2a，um2b，um2c)；

s22-叠加合成三相扰动电压，具体表达式分别为：

s23-应用软件算法拟合，生成与三相扰动电压(uda，udb，udc)频谱相近的三个二进制序列(vda，vdb，vdc)，该序列中频率fm(1≤m≤33)扰动电压分量在同步旋转坐标系下的复数表示为(um1d，um1q)。二进制序列扰动电压注入方式如图5所示，令扰动电压在换流站的并网点与电网电压串联；进一步提取换流站并网点的响应电流可实现换流站交流侧阻抗值的测量。

步骤s3、提取换流站的第一组三相响应电流的具体实现过程包括：

s31-提取响应电流的基波幅值if1，并计算基波电流幅值偏移率l1：

若扰动电压引入的基波电流幅值偏移率l1超过允许值％，则说明选取的扰动量幅值过大，需将扰动量的幅值降低比例％重新进行扰动注入与响应提取。本实施例将值设置为3，值设置为30，即：若当扰动注入引起的基波电流幅值偏移率l1超过3％，则将扰动电压的幅值降低30％重新进行扰动注入与响应提取。

s32-依次提取频率fm(1≤m≤n)的响应电流分量，根据同步旋转坐标系下的谐波幅值和相位以复数形式记录为(im1d，im1q)；

步骤s4、在换流站交流侧串联第二组三相扰动电压，该电压为等效q通道扰动电压，形式为三个周期性二进制序列，具体实现过程包括：

s41-将扰动量幅值设置为交流侧电网电压的5％，结合图3第二个子图所示的扰动合成方法和图4所示的静止-旋转坐标系频率转换关系，依次计算用于合成各频率q通道扰动umq(1≤m≤33)的两组三相量(um3a，um3b，um3c)、(um4a，um4b，um4c)；

s42-叠加合成三相扰动电压，具体表达式分别为：

s43-通过软件算法拟合，生成与三相扰动电压(uqa，uqb，uqc)频谱相近的三个二进制序列(vqa，vqb，vqc)，该序列中频率fm(1≤m≤33)扰动电压分量在同步旋转坐标系下的复数表示为(um2d，um2q)。

步骤s5、提取换流站的第二组三相响应电流具体实现过程包括：

s51-提取响应电流的基波幅值if2，并计算基波电流幅值偏移率l2：

若扰动电压引入的基波电流幅值偏移率l2超过允许值％，则将扰动量的幅值降低％重新进行扰动注入与响应提取。

s52-依次提取频率fm(1≤m≤33)的响应电流分量，根据同步旋转坐标系下的谐波幅值和相位以复数形式记录为(im2d，im2q)；

步骤s6、计算换流站交流侧在测试频率点阻抗值的具体实现过程中采用不同通道阻抗独立计算的形式；首先，d通道阻抗包括zdd(m)和zdq(m)，计算方法为联立以下两个等式：

d通道阻抗计算结果为：

q通道阻抗包括zqd(m)和zqq(m)，计算方法为联立以下两个等式：

q通道阻抗计算结果为：

根据测试的33个频率点阻抗计算结果绘制的阻抗曲线如图6所示，阻抗zdd在低频段存在负阻特性区间，而弱电网存在不可忽略的阻感性电网阻抗，两种阻抗特性的交互作用是换流站-电网互联系统主要的致振诱因；此外，由于换流站交流滤波器组的存在，随着频率的升高阻抗幅值整体呈现下降趋势。

综上所述，本发明实施例的换流站交流侧阻抗在线测试方法可以高效获取具有实时性的hvdc换流站交流侧阻抗特性数据，该数据可应用于换流站-电网互联系统的稳定性在线分析，为高压直流输电系统的运行提供指导。

以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。