CN106952005A - 一种考虑降雨型滑坡地质灾害的电力系统风险评估方法 - Google Patents

Abstract

本发明公开了一种考虑降雨型滑坡地质灾害的电力系统风险评估方法，所述方法依据降雨气象数据信息和输电杆塔周围地质数据信息估算降雨气象诱发滑坡灾害的累积降雨量和滑坡体导致输电杆塔变形的挠曲度；采用曲线拟合法获得输电杆塔损坏联合概率分布；基于网格化技术形成“气象-电网-网格”一体化图层，统计降雨覆盖的输电线路对应的有效输电杆塔数量并计算输电线路失效概率；然后采用非序贯蒙特卡洛法对输电线路进行抽样获得系统失负荷状态；基于交流潮流的最优削减负荷模型计算系统失负荷状态下系统负荷切除量；最后统计系统风险指标。

Description

一种考虑降雨型滑坡地质灾害的电力系统风险评估方法

技术领域

本发明涉及一种电力系统风险评估方法，具体涉及一种考虑降雨型滑坡地质灾害的电力系统风险评估方法。

背景技术

输电杆塔作为电网高压输电线路的支撑体，输电杆塔的完整性直接关系着电网的安全稳定运行和供电可靠性。随着“西电东送”、、“川藏电力联网工程”等战略的实施，越来越多的高压输电网络建设在丘陵、山区以及河流等危险地势，降雨灾害气象极易诱发地质灾害，不仅威胁着降雨灾害区域内人员生命财产安全，而且也威胁着输电杆塔的完整性，对输电线路的正常运行造成严重的影响。根据我国西南地区某电力公司统计资料，暴雨和持续性强降雨天气下已经发生多起输电杆塔损毁事故，包括滑坡地质灾害滑坡体冲撞输电杆塔而变形、输电杆塔地基下沉后倾斜等事故。因此，降雨型滑坡地质灾害引起输电线路断线事故已经受到了电力部门的重点关注，提前预防此类事故的发生显得极其重要。

一般而言，输电杆塔的选址是经过严格规划和帅选后制定的，杆塔具有较高的抗干扰能力。但是，降雨气象下雨水不断入渗到土壤中，土壤趋于饱和状态，最终导致边坡失稳进而引发滑坡地质灾害，威胁着输电杆塔的安全。当输电杆塔处于滑坡体上方时，杆塔下方土体被慢慢掏空易导致杆塔地基下沉；当输电杆塔处于滑坡体下方时，滑坡体从高势能下滑冲撞杆塔易导致杆塔变形甚至倾覆。因此，对输电杆塔周围降雨气象和地质环境的监测是预防此类事故发生的基础。

考虑降雨型滑坡地质灾害的电力系统风险评估方法是潜在滑坡地质灾害威胁输电系统带来系统风险的量化方法，为电网运行人员保证系统稳定运行和供电可靠性提供决策依据。为科学实施输电线路因杆塔损坏而断线的预防措施，首先需建立输电杆塔损坏概率模型，进而建立输电线路断线失效概率模型，然后建立系统风险评估模型，最后根据评估结果制定应急决策，从而避免或减轻降雨型滑坡地质灾害对电网造成的损失。

目前，针对降雨灾害气象诱发地质灾害对输电系统影响的风险评估方法鲜有研究。因此，结合降雨气象模型和滑坡致灾模型进行电网风险评估显得意义重大。

发明内容

本发明目的是提出一种考虑降雨型滑坡地质灾害的电力系统风险评估方法，针对降雨气象诱发地质灾害的特点，分析处于丘陵、山区和河流等位置的输电杆塔损坏情况，为科学合理地防范降雨型地质灾害对电力系统安全运行造成的威胁提供决策依据。

本发明提出的一种考虑降雨型滑坡地质灾害的电力系统风险评估方法，包括以下步骤：

Step1获取输电杆塔所处位置的降雨气象数据信息、地理环境数据信息、滑坡体下滑关联参数和电网数据信息；

Step2估算累积降雨量、滑坡体冲撞输电杆塔等效冲击力和冲撞后输电杆塔挠曲度；

Step3采用曲线拟合法获得降雨诱发岩土滑坡灾害的概率密度函数和滑坡体下滑导致输电杆塔变形的概率密度函数；

Step4建立降雨型滑坡灾害导致输电杆塔损坏的联合概率分布模型，基于网格化技术统计降雨区域内输电杆塔的数量并计算输电线路失效概率；

Step5采用蒙特卡洛法对输电线路状态进行抽样获得系统失负荷概率；

Step6计算系统失负荷状态下的最小切负荷量；

Step7统计系统风险指标。

其中，Step1中所述输电杆塔所处位置的降雨气象数据信息包括滑坡灾害发生时前期降雨的日降雨量和滑坡灾害当日实时降雨量；所述输电杆塔所处位置的地理环境数据信息包括岩体组成、坡体倾斜角度和滑坡体密度；所述滑坡体下滑关联参数包括滑坡体下滑时滑体平均厚度、滑面摩擦系数和滑体冲撞输电杆塔时距离地面的平均高度；所述电网数据信息包括输电线路位置、输电线路长度、输电线路档距、输电杆塔高度、输电杆塔弹性模量、输电杆塔平均截面积和输电杆塔平均宽度。

其中，Step2中所述累积降雨量r(t)的计算公式如式(1)所示：

r(t)＝qr₁+q²r₂+…+q^tr_t (1)

式(1)中r_t表示滑坡地质灾害发生前第t天的降雨量，mm；q为衰减系数，由杆塔所处区域内滑坡体组成性质和降雨蒸腾作用等因素决定。

其中，Step2中所述滑坡体冲撞杆塔等效冲击力F(x)的计算公式如式(2)所示：

式(2)中x表示等效冲击力F作用于杆塔时距离地面的高度，m；m表示滑坡体总质量，kg；g表示重力加速度，为常数；h表示滑坡体质心距离杆塔塔基的垂直高度，m；E表示输电杆塔主材的弹性模量，kN/mm²；S表示杆塔水平平均截面积，m²；a表示滑坡体作用与输电杆塔塔身时对应的宽度，m；μ表示滑面摩擦系数，根据历史统计取值；θ表示坡体倾斜角度，°。

其中，Step2中所述输电杆塔挠曲度w(x)的计算公式如式(3)所示：

其中，Step3所述降雨诱发岩土滑坡灾害的概率密度函数是根据降雨引发滑坡的历史数据，采用皮尔逊Ⅲ型概率密度曲线拟合得到，计算公式如式(4)所示：

式(4)中r表示年降雨极值；r_init表示变量r所能取得的最小值；α表示形状参数；β表示尺度参数；Γ(α)表示α的伽马函数。

其中，Step3中所述滑坡体下滑导致输电杆塔变形的概率密度函数是利用有限元分析软件ANSYS建立滑坡致灾模型进行仿真分析，采用威布尔概率密度曲线拟合得到，计算公式如式(5)所示：

式(5)中λ、γ分别表示尺度参数和形状参数。

其中，Step4中所述降雨型滑坡灾害导致输电杆塔损坏联合概率分布模型的计算公式如式(6)所示：

F(r,w)＝∫∫f(r)f(w)drdw (6)

其中，Step4中所述网格化技术是将降雨气象图层叠加到电网地理信息图层，并形成一个叠加总图层，然后将按比例缩小的网格图层叠加到叠加总图层上，最终形成“气象-电网-网格”一体化图层，网格图层中最小网格单元为：1km×1km。

其中，Step4中所述降雨区域内有效输电杆塔数量n的计算公式如式(7)所示：

式(7)中n_r表示一体化图层中降雨覆盖输电线路的最小网格单元个数；n_t表示一体化图层中输电线路的最小网格单元总数；L_t表示输电线路的实际长度，km；L_d表示输电线路的单位档距，km。

其中，Step4中所述输电线路失效概率的计算公式如式(8)所示：

式(8)中F_l(r,w)表示第l座杆塔损坏的概率。

其中，Step5中所述蒙特卡洛抽样法获得系统失负荷状态是采用非序贯蒙特卡洛法对输电线路进行抽样分析获得系统失负荷状态，系统失负荷概率的计算公式如式(9)所示：

式(9)中N表示系统失负荷状态数；M表示抽样总次数；m(z)表示系统在状态z出现的次数。

其中，Step6中所述系统失负荷状态下最小切负荷量是采用基于交流潮流的最优削减负荷模型，该模型不仅考虑了负荷有功功率和无功功率的切除，而且采用恒定功率因数法，从而得到系统失负荷状态下最小切负荷量，计算公式如式(10)～(17)所示：

0≤P_ci≤P_Li (13)

P_Gimin≤P_Gi≤P_Gimax (14)

Q_Gimin≤Q_Gi≤Q_Gimax (15)

P_Limin≤P_Li≤P_Limax (16)

U_imin≤U_i≤U_imax (17)

式(10)～(17)中I和IL分别表示系统节点和负荷节点数；表示节点i的功率因数角；P_Li表示节点i负荷的有功功率；P_Gi和Q_Gi分别表示节点i的发电机有功功率和无功功率；P_ci表示节点i负荷被切除的有功功率；U_i和U_j分别表示节点i和节点j的电压幅值；G_ij和B_ij分别表示节点i和节点j之间线路的电导和电纳；θ_ij表示节点i和节点j之间的电压相角差；P_Gimax和Q_Gimax分别表示节点i上发电机有功功率和无功功率上限值；P_Gimin和Q_Gimin分别表示节点i上发电机有功功率和无功功率下限值；P_Lmax和P_Limin分别节点i上负荷有功功率上限值和下限值；U_imax和U_imin分别表示节点i的电压上限值和下限值。

其中，Step7中所述系统供电负荷损失期望值的计算公式如式(18)所示：

式(18)中E_loss表示系统供电负荷损失期望值，MW；P_c(z)表示系统在状态z的削减负荷量，MW。

其中，Step7中所述系统供电负荷损失风险的计算公式如式(19)所示：

式(19)中R_risk表示系统供电负荷损失风险；P表示系统总负荷，MW。

与现有技术相比，本发明提出的考虑降雨型滑坡地质灾害的电力系统风险评估方法的优点是：本发明充分考虑了前期降雨对滑坡地质灾害的影响，经过反复仿真校验并对滑坡体下滑冲撞输电杆塔导致变形的过程进行了精细化建模，求得的输电线路失效概率更加符合实际情况，综合降雨气象诱发滑坡地质灾害过程和滑坡体导致输电杆塔变形过程，针对输电杆塔损坏引发输电线路失效导致系统损失负荷过程建立了风险评估模型。该风险评估模型评估结果在一定程度上科学且合理地量化了降雨气象诱发的地质灾害对系统造成风险大小，为电网运行人员在灾害天气下进行调度操作提供理论依据，为提高电力系统抗灾性能而采取更有针对性的预防控制措施。

附图说明

图1为本发明的风险评估流程图；

图2为滑坡体下滑冲撞输电杆塔示意图；

图3为“气象-电网-网格”一体化图层示意图。

具体实施方式

下面结合附图对本发明的具体实施方式进行详细说明。

如图1所示，本发明提出的一种考虑降雨型滑坡地质灾害的电力系统风险评估方法，具体包括以下步骤：

(1)获取输电杆塔所处位置的降雨气象数据信息、地理环境数据信息、滑坡体下滑关联参数和电网数据信息；

从国家气象中心或地方气象中心获取输电杆塔所处位置的降雨气象信息、从国家地质调查局或地方地质调查局获取输电杆塔所处位置的地理环境信息、从国家电网公司或地方电网公司获取电网信息是进行降雨型滑坡地质灾害的电力系统风险评估的基础。

(2)估算累积降雨量、滑坡体冲撞输电杆塔等效冲击力和冲撞后输电杆塔挠曲度；

结合前期降雨指数方法(API)估算累积降雨量，累积降雨量r(t)可按如下方程计算：

r(t)＝qr₁+q²r₂+…+q^tr_t

方程中r_t表示滑坡地质灾害发生前第t天的降雨量，mm；q为衰减系数，由杆塔所处区域内滑坡体组成性质和降雨蒸腾作用等因素决定。

根据滑坡体下滑能量守恒定律估算滑坡体冲撞输电杆塔等效冲击力，等效冲击力F(x)可按如下方程计算：

方程中x表示等效冲击力F作用于杆塔时距离地面的高度，m；m表示滑坡体总质量，kg；g表示重力加速度，为常数；h表示滑坡体在最初位置时其质心距离杆塔塔基的垂直高度，m；E表示输电杆塔主材的弹性模量，kN/mm²；S表示杆塔水平平均截面积，m²；a表示滑坡体作用与输电杆塔塔身时对应的宽度，m；μ表示滑面摩擦系数，根据历史统计取值；θ表示坡体倾斜角度，°。

根据钢体悬臂梁简化方法估算滑坡体冲撞输电杆塔的挠曲度，挠曲度w(x)用来衡量输电杆塔变形程度，可按如下方程计算：

(3)采用曲线拟合法获得降雨诱发岩土滑坡灾害的概率密度函数和滑坡体下滑导致输电杆塔变形的概率密度函数；

根据历史降雨气象统计数据，采用曲线拟合法获得降雨气象引发滑坡地质灾害概率密度函数，f(r)可按如下方程计算：

方程中r表示年降雨极值；r_init表示变量r所能取得的最小值；α表示形状参数；β表示尺度参数；Γ(α)表示α的伽马函数。

在有限元分析软件ANSYS中建立滑坡致灾模型，并采用曲线拟合法获得滑坡体致使输电杆塔变形的概率密度函数，f(w)可按如下方程计算：

方程中λ、γ分别表示尺度参数和形状参数。

(4)建立降雨型滑坡灾害导致输电杆塔损坏的联合概率分布模型，基于网格化技术统计降雨区域内输电杆塔的数量并计算输电线路失效概率；

降雨灾害气象条件下输电线路失效是降雨气象先诱发滑坡地质灾害，然后滑坡体冲撞输电杆塔而导致的。因此，输电线路失效概率模型是条件概率模型，输电线路失效概率F(r,w)可按如下方程计算：

F(r,w)＝∫∫f(r)f(w)drdw

一条完整的输电线路是由多座输电杆塔组成，一座输电杆塔损坏即可造成输电线路断线失效。因此，搜寻降雨气象区域内输电杆塔损坏情况是分析线路失效的基础。为此，采用网格化技术将降雨气象图、电网拓扑图和网格图叠加为“气象-电网-网格”一体化图层，一体化图层中最小网格单元代表1km×1km。则降雨区域内有效输电杆塔数量n可按如下方程计算：

方程中n_r表示一体化图层中降雨覆盖输电线路的最小网格单元个数；n_t表示一体化图层中输电线路的最小网格单元总数；L_t表示输电线路的实际长度，km；L_d表示输电线路的单位档距，km。

输电线路失效的概率是任一输电杆塔损坏的概率，失效概率p_k可按如下方程计算：

方程中F_l(r,w)表示第l座杆塔损坏的概率。

(5)采用蒙特卡洛法对输电线路状态进行抽样获得系统失负荷概率；

首先计算输电线路失效概率，然后采用非序贯蒙特卡洛法对输电线路状态进行抽样，最后获得输电系统失负荷概率，失负荷概率P_loss可按如下方程计算：

方程中N表示系统失负荷状态数；M表示抽样总数；m(z)表示系统在状态z出现的次数。

(6)计算系统失负荷状态下的最小切负荷量；

为了量化降雨型滑坡地质灾害对电力系统造成的风险水平，采用基于交流潮流的最优削减负荷模型计算系统失负荷状态下负荷削减量，为风险指标的计算提供数据基础。系统失负荷状态下负荷削减大小可按如下方程计算：

0≤P_ci≤P_Li

P_Gimin≤P_Gi≤P_Gimax

Q_Gimin≤Q_Gi≤Q_Gimax

P_Limin≤P_Li≤P_Limax

U_imin≤U_i≤U_imax

方程中I和IL分别表示系统节点和负荷节点数；表示节点i的功率因数角；P_Li表示节点i负荷的有功功率；P_Gi和Q_Gi分别表示节点i的发电机有功功率和无功功率；P_ci表示节点i负荷被切除的有功功率；U_i和U_j分别表示节点i和节点j的电压幅值；G_ij和B_ij分别表示节点i和节点j之间线路的电导和电纳；θ_ij表示节点i和节点j之间的电压相角差；P_Gimax和Q_Gimax分别表示节点i上发电机有功功率和无功功率上限值；P_Gimin和Q_Gimin分别表示节点i上发电机有功功率和无功功率下限值；P_Lmax和P_Limin分别节点i上负荷有功功率上限值和下限值；U_imax和U_imin分别表示节点i的电压上限值和下限值。

(7)统计系统风险指标。

根据风险理论可知，降雨型地质灾害对电力系统造成的风险指标就是系统失负荷概率与系统失负荷量的乘积，则系统供电负荷损失期望值可按如下方程计算：

方程中P_c(z)表示表示系统在状态z的削减负荷量，MW。

系统供电负荷损失风险值可按如下方程计算：

方程中P表示系统总负荷，MW。

Claims (9)

1.一种考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，包括以下步骤：

Step1获取输电杆塔所处位置的降雨气象数据信息、地理环境数据信息、滑坡体下滑关联参数和电网数据信息；

Step2估算累积降雨量、滑坡体冲撞输电杆塔等效冲击力和冲撞后输电杆塔挠曲度；

Step3采用曲线拟合法获得降雨诱发岩土滑坡灾害的概率密度函数和滑坡体下滑导致输电杆塔变形的概率密度函数；

Step4建立降雨型滑坡灾害导致输电杆塔损坏的联合概率分布模型，基于网格化技术统计降雨区域内输电杆塔的数量并计算输电线路失效概率；

Step5采用蒙特卡洛法对输电线路状态进行抽样获得系统失负荷概率；

Step6计算系统失负荷状态下的最小切负荷量；

Step7统计系统风险指标。

2.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step1中，所述输电杆塔所处位置的降雨气象数据信息包括滑坡灾害发生时前期降雨的日降雨量和滑坡灾害当日实时降雨量；

所述输电杆塔所处位置的地理环境数据信息包括岩体组成、坡体倾斜角度和滑坡体密度；

所述滑坡体下滑关联参数包括滑坡体下滑时滑体平均厚度、滑面摩擦系数和滑体冲撞输电杆塔时距离地面的平均高度；

所述电网数据信息包括输电线路位置、输电线路长度、输电线路档距、输电杆塔高度、输电杆塔弹性模量、输电杆塔平均截面积和输电杆塔平均宽度。

3.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step2中，所述累积降雨量r(t)可用式(1)计算：

r(t)＝qr₁+q²r₂+…+q^tr_t (1)

式(1)中r_t表示滑坡地质灾害发生前第t天的降雨量，mm；q为衰减系数，由杆塔所处区域内滑坡体组成性质和降雨蒸腾作用等因素决定，

滑坡体冲撞杆塔等效冲击力F(x)可用式(2)计算：

$F\left(x\right)=\sqrt{\frac{6{\mathrm{mghESa}}^2(1-\mathrm{\μ}\·\cot \mathrm{\θ})}{x^3}}---\left(2\right)$

式(2)中x表示等效冲击力F作用于杆塔时距离地面的高度，m；m表示滑坡体总质量，kg；g表示重力加速度，为常数；h表示滑坡体质心距离杆塔塔基的垂直高度，m；E表示输电杆塔主材的弹性模量，kN/mm²；S表示杆塔水平平均截面积，m²；a表示滑坡体作用与杆塔塔身时对应的宽度，m；μ表示滑面摩擦系数，根据历史统计取值；θ表示坡体倾斜角度，°

输电杆塔挠曲度w(x)可将式(2)带入式(3)计算：

$w\left(x\right)=-\frac{F\left(x\right){\mathrm{Hx}}^2}{2{\mathrm{ESa}}^2}+\frac{F\left(x\right)x^3}{4{\mathrm{ESa}}^2}---\left(3\right).$

4.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step3中，所述降雨诱发岩土滑坡灾害的概率密度函数是根据降雨引发滑坡的历史数据，并采用皮尔逊Ⅲ型概率密度曲线计算：

$f\left(r\right)=\frac{{\mathrm{\β}}^{\mathrm{\α}}}{\mathrm{\Γ}\left(\mathrm{\α}\right)}{(r-r_{init})}^{\mathrm{\α}-1}\exp (-\mathrm{\β}(r-r_{init}\left)\right)---\left(4\right)$

式(4)中r表示年降雨极值；r_init表示变量r所能取得的最小值；α表示形状参数；β表示尺度参数；Γ(α)表示α的伽马函数；

滑坡体下滑导致输电杆塔变形的概率密度函数是利用有限元分析软件ANSYS建立滑坡致灾模型进行仿真分析，并采用威布尔概率密度曲线计算：

$f\left(w\right)=\frac{\mathrm{\γ}}{\mathrm{\λ}}{\left(\frac{w}{\mathrm{\λ}}\right)}^{\mathrm{\β}-1}\exp (-{\left(\frac{w}{\mathrm{\λ}}\right)}^{\mathrm{\γ}})---\left(5\right)$

式(5)中λ、γ分别表示尺度参数和形状参数；

降雨型滑坡灾害导致输电杆塔损坏联合概率分布模型是条件概率模型，杆塔损坏的联合概率分布函数F(r,w)可用式(6)计算：

F(r,w)＝∫∫f(r)f(w)drdw (6)。

5.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step4中，所述网格化技术是将降雨气象图层叠加到电网地理信息图层，形成一个叠加图层，然后将按比例缩小的网格图层叠加到叠加图层上，最终形成“气象-电网-网格”一体化图层，网格图层中最小网格单元为：1km×1km，则降雨区域内有效输电杆塔数量n可用式(7)计算：

$n=\frac{n_r}{n_t}\·\frac{L_t}{L_d}---\left(7\right)$

式(7)中n_r表示一体化图层中降雨覆盖输电线路的最小网格单元个数；n_t表示一体化图层中输电线路的最小网格单元总数；L_t表示输电线路的实际长度，km；L_d表示输电线路的单位档距，km。

6.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step4中，所述输电线路失效概率可用式(8)计算：

$p_k=1-\left(\underset{l=1}{\overset{n}{\Π}}\right(1-F_l\left(r,w\right)\left)\right)---\left(8\right)$

式(8)中F_l(r,w)表示第l座杆塔损坏的概率。

7.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step5中，所述采用蒙特卡洛抽样法获得系统失负荷状态是采用非序贯蒙特卡洛法对输电线路进行抽样分析获得系统失负荷状态，系统失负荷概率p_loss可用式(9)计算：

$p_{loss}=\underset{z=1}{\overset{N}{\Σ}}\frac{m\left(z\right)}{M}---\left(9\right)$

式(9)中N表示系统失负荷状态数；M表示抽样总数；m(z)表示系统在状态z出现的次数。

8.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step6中，所述系统失负荷状态下最小切负荷量是采用基于交流潮流的最优削减负荷模型计算输电线路失效后系统削减负荷量，最优削减负荷模型同时考虑了负荷有功功率和无功功率的切除，并采用恒定功率因数法来计算输电线路失效后系统削减负荷量，具体公式为：

$\begin{array}{cc}s.t.& \underset{j\∈i}{\Σ}{P}_{Gi}=U_i\underset{j=1}{\overset{I}{\Σ}}{U}_j(G_{ij}{\mathrm{cos\θ}}_{ij}+B_{ij}{\mathrm{sin\θ}}_{ij})+\underset{i=1}{\overset{IL}{\Σ}}(P_{Li}-P_{ci})\end{array}---\left(11\right)$

0≤P_ci≤P_Li (13)

P_Gimin≤P_Gi≤P_Gimax (14)

Q_Gimin≤Q_Gi≤Q_Gimax (15)

P_Limin≤P_Li≤P_Limax (16)

U_imin≤U_i≤U_imax (17)

式中I和IL分别表示系统节点总数和负荷节点数；表示节点i的功率因数角；P_Li表示节点i负荷的有功功率；P_Gi和Q_Gi分别表示节点i的发电机有功功率和无功功率；P_ci表示节点i负荷被切除的有功功率；U_i和U_j分别表示节点i和节点j的电压幅值；G_ij和B_ij分别表示节点i和节点j之间线路的电导和电纳；θ_ij表示节点i和节点j之间的电压相角差；P_Gimax和Q_Gimax分别表示节点i上发电机有功功率和无功功率上限值；P_Gimin和Q_Gimin分别表示节点i上发电机有功功率和无功功率下限值；P_Lmax和P_Limin分别表示节点i上负荷有功功率上限值和下限值；U_imax和U_imin分别表示节点i的电压上限值和下限值。

9.根据权利要求1所述的考虑降雨型滑坡地质灾害的电力系统风险评估方法，其特征在于，Step7中，所述系统风险指标包括系统供电负荷损失期望值和供电负荷损失风险

系统供电负荷损失期望值E_loss可用式(18)计算：

$E_{loss}=\underset{z=1}{\overset{N}{\Σ}}\frac{m\left(z\right)}{M}{P}_c\left(z\right)---\left(18\right)$

式(18)中P_c(z)表示系统在状态z的削减负荷量，MW

系统供电负荷损失风险R_risk可用式(19)计算：

$R_{risk}=\frac{\underset{z=1}{\overset{N}{\Σ}}\frac{m\left(z\right)}{M}{P}_c\left(z\right)}{P}---\left(19\right)$

式(19)中P表示系统总负荷，MW。