一种电网电压暂降薄弱区域的辨识方法及装置与流程

本发明涉及电力系统电能质量评估领域，特别涉及一种电网电压暂降薄弱区域的辨识方法及装置。

背景技术：

随着现代社会高新技术的发展，电力系统中敏感设备逐渐增多，电网的电压暂降会导致敏感性设备跳闸或者误操作，影响整个终端用户进程，导致巨大的生产和经济损失。电压暂降按照电气与电子工程师协会ieee的定义，是指工频条件下电压均方根减小到0.1～0.9倍额定电压之间、持续时间为0.5周波(以我国工频50hz算，1周波是20ms)至1min的短时间电压变动现象。电压暂降包含两个特征：电压暂降幅值和电压暂降持续时间。因为系统故障是不可预测的，所以本质上讲电压暂降是随机发生的。因此，对电力系统中电压暂降问题进行分析、评估具有十分重要的现实意义。

针对电压暂降问题提出的有效对策和缓解计划需要对系统电压暂降性能评估。估计系统线路上的预期下垂频率(esf)对于了解系统电压暂降性能非常有效。通常，可以使用薄弱区域(aov)和系统故障统计的概念来进行esf估计。aov可以表示导致敏感负荷点的电压幅值跌落至低于给电压阈值的故障范围。通过将系统故障统计量乘以线路的总长度和aov中包含的线路数来计算esf。因此，准确地确定aov对于评估系统电压暂降性能非常重要。

现有技术中，对于薄弱区域的辨识方法，如临界距离法和故障定位法，往往对大规模电力系统的临界点和薄弱区域的计算有很多缺陷或过于耗时，并且在某些情况下会导致在计算临界点时产生分歧，并不能对电网电压暂降薄弱区域进行准确且快速的辨识。因此，如何提供一种运算较少且具有较高的参考可信度的电网电压暂降薄弱区域的辨识方法，以便于对电力系统中电压暂降问题进行分析，是现今急需解决的问题。

技术实现要素：

本发明的目的是提供一种电网电压暂降薄弱区域的辨识方法及装置，以利用黄金分割法和割线迭代法的混合方法，确定电网电压暂降薄弱区域，运算较少且具有较高的参考可信度。

为解决上述技术问题，本发明提供一种电网电压暂降薄弱区域的辨识方法，包括：

根据获取的每条故障线路的故障类型，选择对应的电压方程；其中，所述电压方程为用于计算对应的故障线路上任意一点故障引起的敏感性负荷接入点处的跌落幅值；

根据所述电压方程求解每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值，并计算每条故障线路对应的薄弱性指标向量；其中，所述薄弱性指标向量的数值为对应的故障线路与薄弱区域的关系；

根据所述薄弱性指标向量的数值，利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域。

可选的，所述根据所述电压方程求解每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值，并计算每条故障线路对应的薄弱性指标向量之前，还包括：

通过潮流分析计算每条故障线路的故障前电压；

利用序导纳矩阵获取每条故障线路的序阻抗矩阵。

可选的，所述根据获取的每条故障线路的故障类型，选择对应的电压方程，包括：

若获取的当前故障线路的故障类型为a相单相接地故障，则选择的电压方程为：和其中，s为所述当前故障线路的敏感性负荷接入点，和分别为所述当前故障线路的s点处的a相、b相和c相的相电压幅值，为所述当前故障线路的s点处的故障前电压；为复数运算符，为所述当前故障线路的故障位置k点处的故障前的序电压；和分别为所述当前故障线路的s点与k点之间的正、负和零序互阻抗；和分别为所述当前故障线路的k点处的正、负和零序自阻抗；

若所述当前故障线路的故障类型为b相和c相的两相短路故障，则选择的电压方程为：和

若所述当前故障线路的故障类型为b相和c相的两相短路接地故障，则选择的电压方程为：和

若所述当前故障线路的故障类型为三相短路故障，则选择的电压方程为：其中，为所述当前故障线路的s点处的相电压幅值。

可选的，所述根据所述电压方程求解每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值，并计算每条故障线路对应的薄弱性指标向量，包括：

根据所述当前故障线路对应的电压方程，求解所述当前故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值；其中，若所述当前故障线路对应的电压方程为三相短路故障的电压方程，则所述当前故障线路的首端点故障引起的敏感性负荷接入点处的跌落幅值|f(fi)|和所述当前故障线路的尾端点故障引起的敏感性负荷接入点处的跌落幅值|f(ti)|分别为和若所述当前故障线路对应的电压方程不为三相短路故障的电压方程，则|f(fi)|为和中的最小值，|f(ti)|为和中的最小值，i为所述当前故障线路，fi为所述当前故障线路的首端点，ti为所述当前故障线路的尾端点；

利用所述当前故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值分别减去电压阈值的差与0的比较，确定所述当前故障线路的首端点和尾端点的薄弱性指标矢量为0或1；其中，所述当前故障线路的首端点的薄弱性指标矢量所述当前故障线路的尾端点的薄弱性指标矢量vth为所述电压阈值；

将所述当前故障线路的首端点和尾端点的薄弱性指标矢量相加，确定所述当前故障线路的薄弱性指标向量；其中，所述当前故障线路的薄弱性指标向量

可选的，所述根据所述薄弱性指标向量的数值，利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域，包括：

若所述当前故障线路的薄弱性指标向量的数值为0，则所述当前故障线路中不包含所述薄弱区域；

若所述当前故障线路的薄弱性指标向量的数值为1，则利用所述当前故障线路的首端点和尾端点的跌落幅值和所述当前故障线路的中间点的跌落幅值构造第一二次插值方程，并根据所述第一二次插值方程的根，利用所述割线迭代法计算第一临界点，确定所述当前故障线路中包含的所述薄弱区域；其中，所述第一二次插值方程的根的数量为1，所述第一临界点的数量为1；

若所述当前故障线路的薄弱性指标向量的数值为1，则利用所述黄金分割法求解所述当前故障线路对应的电压方程，获取所述当前故障线路的对应的敏感性负荷接入点处的最大跌落幅值，并判断所述最大跌落幅值是否大于所述电压阈值；

若是，则所述当前故障线路为所述薄弱区域；

若否，利用所述当前故障线路的首端点和尾端点的跌落幅值和所述最大跌落幅值构造第二二次插值方程，并根据所述第二二次插值方程的根，利用所述割线迭代法计算第二临界点，确定所述当前故障线路中包含的所述薄弱区域；其中，所述第二二次插值方程的根的数量为2，所述第二临界点的数量为2。

可选的，所述利用所述黄金分割法求解所述当前故障线路对应的电压方程，获取所述当前故障线路的对应的敏感性负荷接入点处的最大跌落幅值，包括：

根据黄金比例公式确定所述当前故障线路的第一预设点和第二预设点；其中，所述黄金比例公式为所述第一预设点为所述第二预设点为为所述当前故障线路的第三预设点，为所述当前故障线路的第四预设点；

通过所述当前故障线路对应的电压方程，计算所述第一预设点和所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值；其中，若所述当前故障线路对应的电压方程为三相短路故障的电压方程，则所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值和所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值分别为和若所述当前故障线路对应的电压方程不为三相短路故障的电压方程，则为和中的最小值，为和中的最小值；

判断所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值是否大于所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值；

若所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值大于所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值，则令和

若所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值不大于所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值，则令和

判断所述第四预设点减去所述第三预设点的差的绝对值是否小于预设公差；

若所述第四预设点减去所述第三预设点的差的绝对值小于预设公差，则迭代终止，获取所述当前故障线路对应的敏感性负荷接入点处的最大跌落幅值；其中，所述当前故障线路对应的敏感性负荷接入点处的最大跌落幅值为

若所述第四预设点减去所述第三预设点的差的绝对值不小于预设公差，则执行所述通过所述当前故障线路对应的电压方程，计算所述第一预设点和所述第二预设点的跌落幅值的步骤。

可选的，所述根据所述第一二次插值方程的根，利用所述割线迭代法计算第一临界点，包括：

根据所述第一二次插值方程的根，确定第一初始点和第二初始点；其中，所述第一初始点和所述第二初始点分别为所述当前故障线路上的不同的点，且与所述第一二次插值方程的根的距离小于距离阈值；

通过确定所述当前故障线路的估计点；其中，为所述估计点，和分别为所述第一初始点和所述第二初始点，和分别为所述第一初始点和所述第二初始点故障引起的敏感性负荷接入点处的跌落幅值；

判断所述估计点故障引起的敏感性负荷接入点处的跌落幅值减去所述电压阈值的差是否小于预定公差；

若否，则令执行所述通过确定估计点的步骤；

若是，则所述估计点为所述第一临界点。

可选的，所述根据所述薄弱性指标向量的数值，利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域之前，还包括：

存储每条故障线路与对应的电压方程的对应关系和每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值。

此外，本发明还提供了一种电网电压暂降薄弱区域的辨识装置，包括：

选择模块，用于根据获取的每条故障线路的故障类型，选择对应的电压方程；其中，所述电压方程为用于计算对应的故障线路上任意一点故障引起的敏感性负荷接入点处的跌落幅值；

计算模块，用于根据所述电压方程求解每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值，并计算每条故障线路对应的薄弱性指标向量；其中，所述薄弱性指标向量的数值为对应的故障线路与薄弱区域的关系；

确定模块，用于根据所述薄弱性指标向量的数值，利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域。

可选的，该装置还包括：

潮流分析模块，用于通过潮流分析计算每条故障线路的故障前电压；

获取模块，用于利用序导纳矩阵获取每条故障线路的序阻抗矩阵。

本发明所提供的一种电网电压暂降薄弱区域的辨识方法，包括：根据获取的每条故障线路的故障类型，选择对应的电压方程；其中，电压方程为用于计算对应的故障线路上任意一点故障引起的敏感性负荷接入点处的跌落幅值；根据电压方程求解每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值，并计算每条故障线路对应的薄弱性指标向量；其中，薄弱性指标向量的数值为对应的故障线路与薄弱区域的关系；根据薄弱性指标向量的数值，利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域；

可见，本发明通过根据获取的每条故障线路的故障类型，选择对应的电压方程，可以计算每条故障线路的上任意一点故障引起的敏感性负荷接入点处的跌落幅值；通过计算每条故障线路对应的薄弱性指标向量，可以获取每条故障线路的临界点的数量；通过利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域，可以利用黄金分割法和割线迭代法的混合方法，计算临界点，从而确定电网电压暂降薄弱区域，使得运算较少且具有较高的参考可信度，能够有效的用于电力系统电压暂降性能评估分析，对广大电力客户的高质量的生产经营具有重要意义。此外，本发明还提供了一种电网电压暂降薄弱区域的辨识装置，同样具有上述有益效果。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

图1为本发明实施例所提供的一种电网电压暂降薄弱区域的辨识方法的流程图；

图2为本发明实施例所提供的一种电网电压暂降薄弱区域的辨识装置的结构图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

请参考图1，图1为本发明实施例所提供的一种电网电压暂降薄弱区域的辨识方法的流程图。该方法可以包括：

步骤101：根据获取的每条故障线路的故障类型，选择对应的电压方程；其中，电压方程为用于计算对应的故障线路上任意一点故障引起的敏感性负荷接入点处的跌落幅值。

其中，本步骤中的故障线路可以为存在对应的故障类型的故障的线路，故障线路可以为实际使用过程中发生故障的线路，如实施例所提供的方法在线路实际使用过程中发生故障后，将该线路确定为故障线路；也可以为使用本实施例所提供的方法对电网电压暂降薄弱区域进行辨识的过程中，设置发生故障的线路，如实施例所提供的方法还包括对每条线路设置对应故障类型的故障的步骤，使每条线路均成为对应的故障线路。

可以理解的是，本步骤的目的可以为通过每条故障线路的故障类型，获取用于计算对应的故障线路上任意一点故障引起的敏感性负荷接入点处的跌落幅值的电压方程。对于故障类型和电压方程的设置具体，可以由设计人员自行设置，只要可以计算对应的故障类型的故障线路上任意一点故障引起的敏感性负荷接入点处的跌落幅值，也就是故障后敏感性负荷接入点处的相电压幅值，本实施例对此不受任何限制。

具体的，可以使用正序阻抗(z1)，负序阻抗(z2)和零序阻抗(z0)的三个对称分量的系统阻抗模型来计算故障线路上由任意一点故障引起的敏感性负荷接入点处的跌落幅值，如本步骤可以为：若获取的当前故障线路的故障类型为a相单相接地故障，则选择的电压方程为：和其中，s为当前故障线路的敏感性负荷接入点，和分别为当前故障线路的s点处的a相、b相和c相的相电压幅值，为当前故障线路的s点处的故障前电压；为复数运算符，为当前故障线路的故障位置k点处的故障前的序电压；和分别为所述当前故障线路的s点与k点之间的正、负和零序互阻抗；和分别为所述当前故障线路的k点处的正、负和零序自阻抗；

若当前故障线路的故障类型为b相和c相的两相短路故障，则选择的电压方程为：和

若当前故障线路的故障类型为b相和c相的两相短路接地故障，则选择的电压方程为：和

若当前故障线路的故障类型为三相短路故障，则选择的电压方程为：其中，为当前故障线路的s点处的相电压幅值。

其中，当前故障线路可以为本步骤中全部故障线路中的任意一条故障线路。

可以理解的是，根据单相接地故障，两相短路故障，两相短路接地故障和三相短路故障这四种故障类型，可以上述选择对应的电压方程求解故障类型对应的故障线路上的任意一点的相电压幅值。若故障类型为三相短路故障，则故障线路上的任意一点的跌落幅值为相电压幅值若故障类型不为三相短路故障，则故障线路上的任意一点(k)故障引起的敏感性负荷接入点处的跌落幅值为三个相电压幅值和中的最小值。

需要说明的是，对于上述电压方程中的和可以通过如和这三个公式对应的得到。其中，是故障位置k点处的策动点序阻抗，和分别是故障线路的首端点f和尾端点t的策动点序阻抗，是故障线路的首端点f和尾端点t之间的传输序阻抗，是故障线路的首端点f和尾端点t之间的线路序阻抗，是故障线路的任意一点s和故障位置k之间的传输序阻抗，是故障线路的敏感性负荷接入点s和首端点f之间的传输序阻抗，是故障线路的敏感性负荷接入点s和尾端点t之间的传输序阻抗，是故障位置k的故障前的序电压，和分别是故障线路的首端点f和尾端点t的故障前的序电压，并且p是故障位置k与故障线路f-t(0≤p≤1)的整个长度的比例长度。

对应的，为了对上述三个公式进行求解从而获取上述电压方程中的和使用对应的电压方程之前，也就是步骤102之前，还可以包括通过潮流分析计算每条故障线路的故障前电压；利用序导纳矩阵获取每条故障线路的序阻抗矩阵z⁰¹²的步骤，以获取每条故障线路的任意一点的故障前电压和序阻抗矩阵z⁰¹²。

步骤102：根据电压方程求解每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处跌落幅值，并计算每条故障线路对应的薄弱性指标向量；其中，薄弱性指标向量的数值为对应的故障线路与薄弱区域的关系。

可以理解的是，每条故障线路对应的薄弱性指标向量的数值可以为对应的故障线路中的临界点的数量，也就是对应的故障线路与薄弱区域的关系。如若故障线路中的不存在临界点，则故障线路中不包含薄弱区域；若故障线路中的存在一个临界点，则故障线路中包含一个薄弱区域，也就是故障线路的一端在薄弱区域中；若故障线路中的存在两个临界点，则故障线路中包含两个薄弱区域或均在薄弱区域，也就是故障线路的两端均在薄弱区域中。

具体的，本步骤可以为：根据当前故障线路对应的电压方程，求解当前故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值；其中，若当前故障线路对应的电压方程为三相短路故障的电压方程，则当前故障线路的首端点故障引起的敏感性负荷接入点处的跌落幅值|f(fi)|和当前故障线路的尾端点故障引起的敏感性负荷接入点处的跌落幅值|f(ti)|分别为和若当前故障线路对应的电压方程不为三相短路故障的电压方程，则|f(fi)|为和中的最小值，|f(ti)|为和中的最小值，i为当前故障线路，fi为当前故障线路的首端点，ti为当前故障线路的尾端点；

利用当前故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值分别减去电压阈值的差与0的比较，确定当前故障线路的首端点和尾端点的薄弱性指标矢量为0或1；其中，当前故障线路的首端点的薄弱性指标矢量当前故障线路的尾端点的薄弱性指标矢量vth为电压阈值；

将当前故障线路的首端点和尾端点的薄弱性指标矢量相加，确定当前故障线路的薄弱性指标向量；其中，当前故障线路的薄弱性指标向量

可以理解的是，由于每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值在步骤103中同样需要使用，因此，可以存储求解得到的每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值，为了方便对每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值的存储和本步骤中的计算过程，可以通过对应的两个跌落幅值向量获取矩阵形式的每条故障线路的首端点和尾端点的跌落幅值，如，

两个跌落幅值为和其中，vmag,f为每条故障线路的首端点故障引起的敏感性负荷接入点处的跌落幅值向量，vmag,t为每条故障线路的尾端点故障引起的敏感性负荷接入点处的跌落幅值向量，n为全部故障线路对应的数量，i表示当前故障线路为大于等于1小于等于n的正整数。

对应的，每条故障线路的首端点的薄弱性指标矢量、尾端点的薄弱性指标矢量和薄弱性指标向量分别为和其中，

需要说明的是，通过上述方式求解出的当前故障线路的首端点和尾端点的薄弱性指标矢量和分别可以为1或0。对应的当前故障线路的薄弱性指标向量lvii可以为0、1或2，代表当前故障线路的中的临界点的数量。

步骤103：根据薄弱性指标向量的数值，利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域。

具体的，本步骤可以为若当前故障线路的薄弱性指标向量的数值为0，则当前故障线路中不包含薄弱区域；

若当前故障线路的薄弱性指标向量的数值为1，则利用当前故障线路的首端点和尾端点的跌落幅值和当前故障线路的中间点的跌落幅值构造第一二次插值方程，并根据第一二次插值方程的根，利用割线迭代法计算第一临界点，确定当前故障线路中包含的薄弱区域；其中，第一二次插值方程的根的数量为1，第一临界点的数量为1；

若当前故障线路的薄弱性指标向量的数值为1，则利用黄金分割法求解当前故障线路对应的电压方程，获取当前故障线路的对应的敏感性负荷接入点处的最大跌落幅值，并判断最大跌落幅值是否大于电压阈值；

若是，则当前故障线路为薄弱区域；

若否，利用当前故障线路的首端点和尾端点的跌落幅值和最大跌落幅值构造第二二次插值方程，并根据第二二次插值方程的根，利用割线迭代法计算第二临界点，确定当前故障线路中包含的薄弱区域；其中，第二二次插值方程的根的数量为2，第二临界点的数量为2。

其中，利用黄金分割法求解当前故障线路对应的电压方程，获取当前故障线路的对应的敏感性负荷接入点处的最大跌落幅值的步骤，可以包括：根据黄金比例公式确定所述当前故障线路的第一预设点和第二预设点；其中，所述黄金比例公式为所述第一预设点为所述第二预设点为为所述当前故障线路的第三预设点，为所述当前故障线路的第四预设点；

通过所述当前故障线路对应的电压方程，计算所述第一预设点和所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值；其中，若所述当前故障线路对应的电压方程为三相短路故障的电压方程，则所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值和所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值分别为和若所述当前故障线路对应的电压方程不为三相短路故障的电压方程，则为和中的最小值，为和中的最小值；

判断所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值是否大于所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值；

若所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值大于所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值，则令和

若所述第一预设点故障引起的敏感性负荷接入点处的跌落幅值不大于所述第二预设点故障引起的敏感性负荷接入点处的跌落幅值，则令和

判断所述第四预设点减去所述第三预设点的差的绝对值是否小于预设公差；

若所述第四预设点减去所述第三预设点的差的绝对值小于预设公差，则迭代终止，获取所述当前故障线路对应的敏感性负荷接入点处的最大跌落幅值；其中，所述当前故障线路对应的敏感性负荷接入点处的最大跌落幅值为

若所述第四预设点减去所述第三预设点的差的绝对值不小于预设公差，则执行所述通过所述当前故障线路对应的电压方程，计算所述第一预设点和所述第二预设点的跌落幅值的步骤。

其中，根据第一二次插值方程的根，利用割线迭代法计算第一临界点的步骤，可以包括：

根据第一二次插值方程的根，确定第一初始点和第二初始点；其中，第一初始点和第二初始点分别为当前故障线路上的不同的点，且与第一二次插值方程的根的距离小于距离阈值；

通过确定当前故障线路的估计点；其中，为估计点，和分别为第一初始点和第二初始点，和分别为第一初始点和第二初始点故障引起的敏感性负荷接入点处的跌落幅值；

判断估计点故障引起的敏感性负荷接入点处的跌落幅值减去电压阈值的差是否小于预定公差；

若否，则令执行通过确定估计点的步骤；

若是，则估计点为第一临界点。

可以理解的是，对于利用当前故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值和当前故障线路的中间点故障引起的敏感性负荷接入点处的跌落幅值(|f(0.5)|)构造的第一二次插值方程和利用当前故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值和最大跌落幅值构造的第二二次插值方程的具体构造过程以及求解过程，可以采用与现有技术相似的方式，本实施例对此不受任何限制。

需要说明的是，对于根据第二二次插值方程的根，利用割线迭代法计算第二临界点的过程，可以采用与根据第一二次插值方程的根，利用割线迭代法计算第一临界点的步骤相似的方式，获取第二临界点对应的两个点。

需要说明的是，本实施例的目的是为了获取每条故障线路对应的薄弱区域，对于当前故障线路来说，在确定当前故障线路对应的薄弱区域后，可以通过判断是否已确定全部故障线路对应的薄弱区域的步骤，确保可以达到本实施例的目的。对应的，若未确定全部故障线路对应的薄弱区域，则选择下一故障线路作为当前故障线路，确定该当前故障线路对应的薄弱区域。本实施例对此不做任何限制。

本实施例中，本发明实施例通过根据获取的每条故障线路的故障类型，选择对应的电压方程，可以计算每条故障线路的上任意一点故障引起的敏感性负荷接入点处的跌落幅值；通过计算每条故障线路对应的薄弱性指标向量，可以获取每条故障线路的临界点的数量；通过利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域，可以利用黄金分割法和割线迭代法的混合方法，计算临界点，从而确定电网电压暂降薄弱区域，使得运算较少且具有较高的参考可信度，能够有效的用于电力系统电压暂降性能评估分析，对广大电力客户的高质量的生产经营具有重要意义。

请参考图2，图2为本发明实施例所提供的一种电网电压暂降薄弱区域的辨识装置的结构图。该装置可以包括：

选择模块100，用于根据获取的每条故障线路的故障类型，选择对应的电压方程；其中，电压方程为用于计算对应的故障线路上任意一点故障引起的敏感性负荷接入点处的跌落幅值；

计算模块200，用于根据所述电压方程求解每条故障线路的首端点和尾端点故障引起的敏感性负荷接入点处的跌落幅值，并计算每条故障线路对应的薄弱性指标向量；其中，所述薄弱性指标向量的数值为对应的故障线路与薄弱区域的关系；

确定模块300，用于根据薄弱性指标向量的数值，利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域。

可选的，该装置还可以包括：

潮流分析模块，用于通过潮流分析计算每条故障线路的故障前电压；

获取模块，用于利用序导纳矩阵获取每条故障线路的序阻抗矩阵。

本实施例中，本发明实施例通过选择模块100根据获取的每条故障线路的故障类型，选择对应的电压方程，可以计算每条故障线路的上任意一点故障引起的敏感性负荷接入点处的跌落幅值；通过计算模块200计算每条故障线路对应的薄弱性指标向量，可以获取每条故障线路的临界点的数量；通过确定模块300利用黄金分割法和/或割线迭代法确定每条故障线路对应的薄弱区域，可以利用黄金分割法和割线迭代法的混合方法，计算临界点，从而确定电网电压暂降薄弱区域，使得运算较少且具有较高的参考可信度，能够有效的用于电力系统电压暂降性能评估分析，对广大电力客户的高质量的生产经营具有重要意义。

说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。

结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(ram)、内存、只读存储器(rom)、电可编程rom、电可擦除可编程rom、寄存器、硬盘、可移动磁盘、cd-rom、或技术领域内所公知的任意其它形式的存储介质中。

以上对本发明所提供的电网电压暂降薄弱区域的辨识方法及装置进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。