一种考虑相位跳变的电力系统电压暂降水平评估方法
【技术领域】
[0001] 本发明涉及电力系统电能质量评估领域，具体涉及一种考虑相位跳变的电力系统 电压暂降水平评估方法。
【背景技术】
[0002] 电压暂降是电力系统不可避免的电能质量扰动事件，被认为是最严重的电能质量 问题。随着越来越多对电压暂降敏感的用电设备，以及风电、太阳能发电等间歇性电源接入 电网，电压暂降所造成的经济损失和社会影响日益增长，引起了供电企业、电力用户以及学 术界的广泛关注。对电力系统或关心节点进行电压暂降水平评估，有助于供电企业和电力 用户预知潜在的暂降危害，采取必要的措施规避或缓解由电压暂降造成的风险。
[0003] 电压暂降评估从层次上可分为单一暂降事件评估、节点评估以及系统评估。其中， 系统评估以节点评估为基础，节点评估又建立在对单一暂降事件评估的基础上，换言之，节 点评估和系统评估是以单一暂降事件为基础的统计描述。因此，提取单一暂降事件的重要 特征是关键。
[0004] 已被提出的单一暂降事件特征有暂降幅值、持续时间、相位跳变、波形点等。在以 往的研究中，暂降幅值和持续时间是最受关注的两个特征量。为描述某个暂降事件的严重 程度，提出了暂降能量指标、严重度指标;为描述敏感设备与暂降的兼容程度，制定了电压 耐受曲线;在节点和系统层面，提出了电压暂降幅值-持续时间统计表格。然而，越来越多的 研究表明，许多设备对相位跳变同样是十分敏感的，如可调速设备、交流接触器、整流器、风 力发电机等。现有的暂降指标、统计表格和参考曲线，均只考虑了暂降的幅值和持续时间， 对于需要考虑相位跳变的情况无法适用。
[0005]造成暂降的主要原因包括电力系统短路故障、大电机启动、变压器励磁涌流等。其 中，系统短路故障是最主要原因。发生短路故障时，相位跳变伴随幅值的降低出现，因此从 电气联系上来说，相位跳变与幅值的关系，比持续时间和幅值的关系更为紧密。现有的暂降 评估方法无法体现幅值和相位跳变的关系，无法揭示二者范围以及变化规律。
[0006] 与本发明有关的现有技术及其缺点如下，现有的暂降水平评估方法主要包括实测 统计法、仿真模拟法以及状态估计法。
[0007] 实测统计法:实测统计法通过在现场安装监测装置，监测安装节点一个周期内的 暂降情况。该方法的好处在于监测结果能可靠反映监测周期的情况。存在的缺点为：1)电力 系统复杂庞大，在全网安装需耗费巨大的购置设备、搭建通信通道、后期维护等费用；2) - 个不长的监测周期，无法消除该监测周期内情况的个异性，不能客观、稳定反映该处的暂降 平均水平;3)对于还未安装监测装置的节点，实测统计将耗费较大的时间成本。
[0008] 仿真模拟法:仿真模拟法通过仿真模拟软件，参考实际电网的故障率及故障条件 进行短路计算。主要方法包括临界距离法、故障点法、解析法。临界距离法适用于简单的辐 射型网络，对复杂网络不适用;故障点法中故障点的设置对结果有较大影响;解析法物理意 义明确，结果准确，但解析曲线交点的求解较为困难。其次，这三种方法均只仅能求出对应 幅值区间的暂降频次，无法反映相位跳变情况。
[0009] 状态估计法:状态估计法可以通过状态估计方程，基于监测节点的暂降频次，估计 未监测节点的暂降频次。该方法同样无法排除特定监测周期的个异性，且结果仅有幅值和 频次信息，无法反映相位跳变情况。
【发明内容】
[0010] 有鉴于此，本发明的目的是提供一种考虑相位跳变的电力系统电压暂降水平评估 方法，考虑故障线路阻抗与系统阻抗之间的阻抗角、故障距离、故障类型、变压器对暂降类 型的变换作用等不确定因素，基于蒙特卡罗模拟方法进行多次短路计算，然后计算关心节 点的暂降幅值和相位跳变。引入考虑暂降幅值和相位跳变的统计表格及统计图，将蒙特卡 罗模拟结果表示在图表中，便于直接观察和细致比较。
[0011] 本发明采用以下方案实现：一种考虑相位跳变的电力系统电压暂降水平评估方 法，具体包括以下步骤：
[0012] 步骤S1:确定蒙特卡罗模拟的故障条件；
[0013] 步骤S2:根据步骤S1中确定的故障条件，基于电压暂降分压模型，计算暂降幅值和 相位跳变角；
[0014] 步骤S3:将步骤S2中的计算结果制成暂降幅值和相位跳变初步统计表格以及细分 统计表格；
[0015] 步骤S4:将计算结果形成考虑幅值和相位跳变的图形化表示。
[0016] 进一步地，所述步骤S1通过产生随机数，模拟在不同阻抗角、故障距离、故障类型、 变压器对暂降的变换作用等情况，确定蒙特卡罗模拟的故障条件，具体包括以下步骤：
[0017] 步骤S11:生成随机数确定阻抗角：阻抗角α为故障阻抗ZF与系统阻抗Zs之间的夹 角，如式（1)所示，其中 Zs = Rs+jXs，Zf=Rf+jXF，
[0019] 由于阻抗角绝大多数位于-60°~10°之间，故生成随机数x模拟每次故障可能的阻 抗角RIAi，如式(2)所示：
[0020] RIAi = xX70-60 (2)
[0021] 其中，X服从(〇，1)均匀分布。
[0022] 步骤S12:确定短路故障类型：电力系统中的短路故障分为三相对称短路、单相接 地短路、两相短路、两相接地短路，且四种短路故障类型发生的概率依次为5%、65%、10%、 20%，生成随机数y模拟每次发生的短路类型RFh，如式(3)所示：
[0024]其中，1^1^ = 1,2,3,4分别表示发生三相对称短路、单相接地短路、两相短路、两相 接地短路。
[0025]步骤S13:生成随机数确定故障位置:故障位置用故障点到关心节点的距离标幺值 P表示，P服从(〇，1)均匀分布;对于长度为L的线路，故障点到PCC点的距离可表示为l=pL。 [0026]步骤S14:生成随机数确定变压器对暂降类型的变换作用:对称和不对称短路故障 造成的电压暂降分为7种类型，去除含有零序分量的种类型，设备侧能感受到的暂降类型有 5种。其中，三相短路故障导致A类暂降，单相接地短路故障导致C、D类暂降，两相短路故障导 致C、D类暂降，以及两相接地短路故障导致F、G类暂降，生成随机数z服从(0，1)分布。
[0027]在暂降类型的模拟中，若所述步骤S12中确定短路类型为三相短路，不论z为何值， 变换结果都为A类暂降；若确定短路类型为单相接地短路或两相短路，当0〈z < 0.5时为C类 暂降，当0.5 < z〈l时为D类暂降；若确定短路类型为两相接地短路，当0〈z < 0.5时为F类暂 降，当0.5<z〈l时为G类暂降。
[0028] 进一步地，所述步骤S2具体包括以下步骤：
[0029] 步骤S21:根据生成的故障条件进行短路计算:三相对称故障可直接应用分压公式 求解;不对称故障需应用对称分量法，首先利用各边界条件求出正、负、零序电压，再叠加旋 转因子作用合成为A、B、C三相电压相量，分别计算各相故障时电压幅值Uaf、Ubf、Ucf和相角 Φα、Φβ、Φ。；
[0030] 步骤S22:计算暂降幅值和相位跳变角：以故障前Α相电压为0°、Β相电压为-120°、C 相电压为120°、三相幅值均为l.Opu为基准，将计算所得A、B、C三相电压相量与故障前基准 电压相量进行比较，求出暂降幅值和相位跳变角；其中暂降幅值为故障时电压幅值与故障 前电压幅值之比，相位跳变角为故障时电压相角和故障前电压相角之差，如式(4)所示：
[0032] 进一步地，所述步骤S3具体包括以下步骤：
[0033]步骤S31:建立初步统计表格:根据IEEE对电压暂降幅值的定义，将暂降幅值划分 为小于0 · lpu、0 · 1至0 · 2pu、0 · 3至0 · 4pu、0 · 5至0 · 6pu、0 · 6至0 · 7pu、0 · 7至0 · 8pu以及0 · 8至 0.9pu这9个区间；根据实际可能出现的相位跳变范围，将相位跳变划分为大于-70°、-70° 至-50°、-50° 至-30°、-30° 至-10°、-10° 至 10°、10° 至30°、30° 至50°、50° 至70° 以及大于70° 这9个区间。
[0034] 步骤S32:建立细分统计表格：将幅值划分为0.6至0.7pu、0.7至0.8pu以及0.8至 0.9pu这3个区间，将相位跳变划分为-30°至-25°、-25°至-20°、-20°至-15°、-15°至-10°、-10°至-5°、-5°至0°、0°至5°以及5°至10°这8个区间。
[0035]进一步地，所述步骤S4建立考虑暂降幅值和相位跳变的三维频次图及二维等高线 图，用以直观体现暂降幅值和相位跳变之间的关系，以及各区域暂降出现的频率规律。
[0036]与现有技术相比，本发明的有益效果是:本发明提供了一种同时考虑幅值和相位 跳变的电压暂降评估方法。基于蒙特卡罗仿真方法，模拟了实际系统中阻抗角、短路故障类 型、故障发生位置、变压器对暂降类型变换等诸多不确定因素，评估关心节点的暂降幅值和 相位跳变水平。首次提出考虑幅值和相位跳变的统计表格，并将可能的幅值和相位跳变范 围进行了初步划分和细致划分，能从整体和局部两个层面认识和比较关心节点的暂降水 平。幅值和相位跳变三维频次图和二维等高线图能给人更加直观的观察和判断。因此，本发 明能为供电企业、电力用户等提供较完整的考虑幅值和相位跳变的暂降评估方法和呈现方 法，为治理、投资提供决策依据。
【附图说明】
[0037] 图1为本发明的电压暂降分压模型；
[0038] 图2为本发明的设备终端感受的电压暂降类型；
[0039] 图3为本发明的基于蒙特卡罗模拟的考虑相位跳变的电压暂降水平评估流程图；
[0040] 图4为本发明的某实际电网电气接线图；
[00411图5为本发明的STU-11节点三维频次图；
[0042]图6为本发明的STU-11节点二维等高线图。
【具体实施方式】
[0043] 下面结合附图及实施例对本发明做进一步说明。
[0044] 本实施提供一种考虑相位跳变的电力系统电压暂降水平评估方法，具体包括以下 步骤：
[0045] 步骤S1:确定蒙特卡罗模拟的故障条件；
[0046] 步骤S2:根据步骤S1中确定的故障条件，基于电压暂降分压模型，计算暂降幅值和 相位跳变角；
[0047] 步骤S3:将步骤S2中的计算结果制成暂降幅值和相位跳变初步统计表格以及细分 统计表格；
[0048]步骤S4:将计算结果形成考虑幅值和相位跳变的图形化表示。
[0049] 在本实施例中，所述步骤S1通过产生随机数，模拟在不同阻抗角、故障距离、故障 类型、变压器对暂降的变换作用等情况，确定蒙特卡罗模拟的故障条件，具体包括以下步 骤：
[0050] 步骤SI 1:生成随机数确定阻抗角：电压暂降分压模型如图1所示。阻抗角α为故障 阻抗ZF与系统阻抗Zs之间的夹角，如式（1)所示，其中Zs = Rs+jXs，ZF = RF+jXF，
[0052] 由于阻抗角绝大多数位于-60°~10°之间，故生成随机数x模拟