[发明专利]数据驱动的直流受端电网动态无功储备需求评估方法

权利要求书

1.一种数据驱动的直流受端电网动态无功储备需求评估方法，其特征在于该方法包括以下步骤：

(1)建立直流受端电网的离线样本库，离线样本库中的每个样本包括直流受端电网的特征量和动态无功储备需求量，定义直流受端电网为电网中的一个从直流输电线接收电能的区域电网，具体包括以下步骤：

(1-1)提取直流受端电网中直流逆变换流站近区的暂态电压稳定特征量，过程如下：

(1-1-1)确定直流逆变换流站近区的范围：

定义电网中的PQ节点为有功功率注入和无功功率注入给定，PQ节点为电压幅值和电压相角待求，设定电网中500kV母线和逆变换流站母线为PQ节点，根据电网潮流方程的快速分解法，简化的修正方程的PQ节点的无功功率为：

-B”ΔV＝ΔQ/V

其中，B”是电网中不包括PV节点时节点的导纳矩阵的虚部，ΔV是电网中不包括PV节点时节点的电压幅值的修正量，ΔQ是电网中不包括PV节点时的节点无功功率注入的修正量，修正量的取值为任意实数，V是电网中不包括PV节点时的点电压幅值；

把上式增广到包含PV节点和PQ节点，则增广后的修正方程PQ节点的无功功率为：

-B”⁰ΔV＝ΔQ⁰/V

其中，B”⁰是电网中包含PV节点和PQ节点时节点导纳矩阵的虚部，ΔQ⁰是电网中包含PV节点和PQ节点时节点无功功率注入的修正量，修正量的取值为任意实数；

计算PQ节点与PV节点之间的无功功率-电压灵敏度矩阵S_VQ如下：

S_VQ＝ΔV/ΔQ⁰＝-(B”⁰)^-1/V

根据上述无功功率-电压灵敏度矩阵S_VQ，得到逆变换流站母线与电网中500kV母线之间的无功功率-电压灵敏度向量S'_VQ：

其中，S_Bus500是电网中500kV母线的集合，下标Bus500_i是集合S_Bus500中的第i个元素，下标n₅₀₀是S_Bus500中的元素数量，是逆变换流站母线与集合S_Bus500中的第i个元素Bus500_i之间的无功功率-电压灵敏度；

定义|S'_VQ|为S'_VQ里各个元素的绝对值的集合|S'_VQ|：

定义|S'_VQ|_0.90大于|S'_VQ|中90％的元素，但小于|S'_VQ|中10％的元素，并定义S_NZ500为直流逆变换流站近区的母线集合，即：

定义电网中电压等级低于500kV的母线为低压母线，则电网中低压母线对S_Bus500中母线的无功功率-电压灵敏度矩阵S″_VQ为：

S_Bus'＝{Bus₁,...,Bus_n'}

其中，S_Bus'电网中低压母线的集合，下标n'是S_Bus'中的元素数量，Bus_i是S_Bus'的第i个元素，为S_Bus'的第j个元素对S_Bus500的第i个元素的无功功率-电压灵敏度；

定义|S″_VQ|为S″_VQ中各个元素的绝对值的集合|S″_VQ|：

定义为S_Bus500的第i个元素所代表的S_Bus'中的元素的集合，并定义S_Bus500,Bus'为全部组成的集合，即：

其中，是|S″_VQ|中每一列的最大值，

定义S_{BusNZ500,Bus'}为S_Bus500,Bus'中与直流逆变换流站近区的母线集合S_NZ500相对应的部分，即：

根据上述步骤中的定义和表达式，得到直流逆变换流站近区的范围S_NZ为：

S_NZ＝S_NZ500∪S_{BusNZ500,Bus'}∪S_Inv

其中，S_Inv是电网中逆变换流站母线的集合，Bus_Inv,i是S_Inv中的第i个元素，下标n_Inv是S_Inv中的元素数量；

(1-1-2)提取上述直流逆变换流站近区的暂态电压稳定特征量：

定义S_NZ,F为直流逆变换流站近区的特征量的集合S_NZ,F为：

其中，P_g是电网中发电机母线的有功出力，Q_g是电网中发电机母线的无功出力，P_l是电网中负荷母线的有功负荷量，Q_l是电网中负荷母线的无功负荷量，V是电网中母线的电压幅值，P_d是电网中逆变换流站母线接收的有功功率；

(1-2)定义量化评估上述直流逆变换流站近区暂态电压稳定性指标I_STVS如下：

I_STVS＝T_span,max/T_th

其中，T_th为v(t)低电压阈值，T_end为v(t)结束的时刻，V_th为v(t)低电压持续时间阈值，v(t)为母线的电压轨迹，T_span,max是v(t)中电压持续低于V_th的最长时间段，I_STVS是T_span,max与T_th的比值，I_STVS的取值范围是[0,+∞)；

(1-3)利用数值仿真迭代方法，计算上述直流逆变换流站近区的动态无功功率储备需求，过程如下：

迭代计算的步骤如下：

第1步，设定直流受端电网的运行方式，用数值仿真方法计算直流受端电网的潮流方程，得到潮流计算结果；

第2步，定义S_Q为直流受端电网中直流逆变换流站近区当前的动态无功功率源的容量，记录直流逆变换流站近区当前S_Q，并根据上述(1-3)中的潮流计算结果，计算上述步骤(1-1-2)中的直流逆变换流站近区的特征量的集合S_NZ,F，记录计算得到的当前S_NZ,F；

第3步，用数值仿真方法，计算得到直流逆变换流站近区S_NZ中母线的电压轨迹v(t)，并根据电压轨迹v(t)，计算上述步骤(1-2)的直流逆变换流站近区暂态电压稳定性指标I_STVS，记录计算得到的当前I_STVS；

第4步，取I_STVS中的最大值，记为max(I_STVS)，计算1与max(I_STVS)的差的绝对值，对该绝对值进行判断，若|(max(I_STVS)-1)|0.01，则判定直流受端电网中直流逆变换流站近区当前的动态无功功率源的需求量为S_Q，停止迭代计算，若|(max(I_STVS)-1)|≥0.01，则进行第五步；

第5步，设定一个修正系数α，α＞0，根据修正系数，对上述(1-3)中的当前动态无功功率源的容量S_Q进行修正，使动态无功功率源中的每个动态无功源容量的修正量与相应动态无功功率源的容量成正比，所有动态无功源容量的修正量的总和为α(max(I_STVS)-1)，返回上述第2步；

(1-4)建立直流受端电网的离线样本库，步骤如下：

第1步，将一天分隔为96个相等的时间间隔，记为96个时刻，设定96个时刻直流受端电网的运行方式为直流受端电网的设定运行方式；

第2步，以96个设定运行方式为中心，分别在96个设定运行方式的邻域内随机选择100个运行方式，分别计算9600个运行方式下的直流逆变换流站近区的特征量的集合S_NZ,F和直流逆变换流站近区当前的动态无功功率源容量S_Q；

设定直流受端电网的运行方式为Oper_i，将Oper_i的邻域记为U(Oper_i)：

其中，S_Bus是直流受端电网中母线的集合，是S_Bus中的第m个元素在Oper_i下的有功功率，是S_Bus中的第m个元素在Oper_i下的无功功率，是S_Bus中的第m个元素在Oper_i下的发电机有功出力，是S_Bus中的第m个元素在Oper_i下的发电机电压幅值，是全部的绝对值平均值，是全部的绝对值平均值，是全部的绝对值平均值，是全部的绝对值平均值；

第3步，将上述(1-4)中的直流逆变换流站近区的特征量的集合S_NZ,F和直流逆变换流站近区当前的动态无功功率源容量S_Q存储到直流受端电网的离线样本库，并为每个样本附加一个Mark：

Mark＝{Oper_i,Case_j}

其中，Mark是该样本的标识，Oper_i是该样本所对应的设定运行方式编号，Case_j是该样本在相应的设定运行方式所生成的样本集合中的编号；

第4步，将直流受端电网的离线样本库记为S_spl，构建S_spl，S_spl中每个样本的表述如下：

(2)根据上述步骤(1-4)的离线样本库，建立数据挖掘模型，并训练数据挖掘模型，具体包括以下步骤：

(2-1)建立以上述步骤(1-1)中的直流逆变换流站近区暂态电压稳定特征量为输入，以上述步骤(1-1)中的直流逆变换流站近区的动态无功功率容量为输出的支持向量回归模型，将上述步骤(1)中直流受端电网的离线样本库S_spl中的样本按照相应的运行方式分成96类，96类样本分别对应一个子支持向量回归模型，总支持向量回归模型是96个子支持向量回归模型的并联，总支持向量回归模型、子支持向量回归模型的输入为S_NZ,F，总支持向量回归模型、子支持向量回归模型的输出为S_Q；

(2-2)对上述支持向量回归模型进行训练和调优，包括以下步骤：

(2-2-1)对上述步骤(1)的直流受端电网的离线样本库S_spl的样本数据进行预处理，包括以下步骤：

第1步，将S_spl中的样本按照相应运行方式分成96类：

S_spl＝{s_spl,1,...,s_spl,96}

其中，s_spl,i是Oper_i对应的样本子集合；

第2步，计算s_spl,i中的样本的平均值μ(s_spl,i)，计算s_spl,i中的样本的标准差σ(s_spl,i)；

第3步，得到预处理后的样本z：

x∈s_spl,i

其中，x是s_spl,i中的一个样本，z的输入为S_NZ,F'，输出为S_Q'；

(2-2-2)训练支持向量回归模型：

定义一个支持向量回归模型的评价函数Cost：

Cost＝|(S_Q'-(S_Q')_predict)/S_Q'|，

其中，(S_Q')_predict是支持向量回归模型的输出值，S_Q'是直流逆变换流站近区的动态无功储备需求的实际值；

利用评价函数Cost，计算得到支持向量回归模型的评价值，根据上述评价值，将步骤(2-1)中样本集S_spl分成训练集、交叉检验集和测试集3个部分，其中，训练集占70％，交叉检验集占15％，测试集占15％；

(2-2-3)调优支持向量回归模型，包括以下步骤：

第1步，设定c的松弛系数C，C的取值范围为C_min≤C≤C_max，设定支持向量回归模型的核函数参数为γ，γ的取值范围γ_min≤γ≤γ_max；

第2步，定义为C的取值上限和下限的平均值，为γ的取值上限和下限的平均值，对松弛系数C和核函数参数γ进行判断，若或者则结束调优，并将当前的(C,γ)作为支持向量回归模型最优的松弛系数C和核函数参数γ；若且则进行第3步；

第3步，定义Cost_(C,γ)为一组与上述(2-2-3)中的(C,γ)相对应的支持向量回归模型的Cost，对分别训练支持向量回归模型，然后分别计算Cost_(C,γ)；

第4步，分别计算和找出其中的最大者，将(C,γ)的取值范围进行缩减，若最大者是则减小(C,γ)的取值范围，使返回上述第1步；

(3)周期性地更新上述步骤(1)的直流受端电网的离线样本库S_spl，并周期性地更新支持向量回归模型，具体包括：

(3-1)直流受端电网离线样本库的周期性更新方法：

设定离线样本库的更新周期为15分钟，周期性离线样本库的具体步骤如下：

第1步，记当前直流受端电网的运行方式为Oper₀，当前时刻为T₀，记录Oper₀和T₀；

第2步，记运行方式Oper₀邻域中的运行方式构成集合为在Oper₀邻域中随机选取100个运行方式，构成集合

第3步，根据步骤(1-1-2)和步骤(1-3)，计算的各运行方式的直流逆变换流站近区特征量S_NZ,F和直流逆变换流站近区的动态无功功率源容量S_Q，并将的各运行方式的直流逆变换流站近区特征量S_NZ,F和直流逆变换流站近区的动态无功功率源容量S_Q添加到当前直流受端电网样本库中；

第4步，从直流受端电网离线样本库中，删去与当前时刻T₀之前24小时的运行方式相对应的样本；

第5步，等待到T₀+15分钟时刻，返回上述(3-1)中的第1步；

(3-2)支持向量回归模型的周期性更新方法：

设定支持向量回归模型的更新周期为15分钟；

更新支持向量回归模型的步骤如下：

第1步，记当前时刻为T₀，读取离线样本库S_spl,0；

第2步，利用S_spl,0，训练并更新T₀对应的子支持向量回归模型；

第3步，更新总支持向量回归模型；

第4步，等待到T₀+15分钟时刻，返回上述(3-2)中的第1步；

(4)利用步骤(2)和步骤(3)训练得到的支持向量回归模型，在线计算直流受端电网动态无功功率源容量，具体包括：

(4-1)在线获取直流受端电网当前运行状态：

从电网的能量管理系统中获取S_NZ中的母线的电气量S_StateVar，表述如下：

S_StateVar＝{StateVar_i|i∈S_NZ}

StateVar_i＝{P_g,i,Q_g,i,P_l,i,Q_l,i,V_i,θ_i,P_d},i∈S_NZ

其中，P_g,i是直流受端电网中母线i的发电机有功功率，Q_g,i是直流受端电网中母线i的发电机无功功率，P_l,i是直流受端电网中母线i的有功功率，Q_l,i是直流受端电网中母线i的无功功率，V_i是直流受端电网中母线i的电压幅值，θ_i是直流受端电网中母线i的电压相角，P_d是直流受端电网中逆变换流站母线接收的有功功率；

S_StateVar将用于后续步骤中提取直流受端电网中直流逆变换流站近区的暂态电压稳定特征量；

(4-2)在线提取直流受端电网中直流逆变换流站近区的暂态电压稳定特征量：

定义S_NZ,F为直流逆变换流站近区的特征量的集合S_NZ,F为：

其中，P_g是电网中发电机母线的有功出力，Q_g是电网中发电机母线的无功出力，P_l是电网中负荷母线的有功负荷量，Q_l是电网中负荷母线的无功负荷量，V是电网中母线的电压幅值，P_d是电网中逆变换流站母线接收的有功功率；

(4-3)在线计算直流逆变换流站近区的动态无功功率源容量，包括以下步骤：

(4-3-1)计算直流受端电网运行方式之间的距离

其中，是发电机有功功率的权重，是发电机无功功率的权重，是直流受端电网有功功率的权重，是直流受端电网无功功率的权重，ω_V是电压幅值的权重，ω_θ是电压相角的权重，是逆变换流站接收的有功功率的权重，下标Oper_i代表运行方式i，下标Oper_j代表运行方式j，P_g,k是S_NZ500中母线的发电机有功出力，Q_g,k是S_NZ500中母线的发电机无功出力，P_l,k是S_NZ500中母线的负荷有功量，Q_l,k是S_NZ500中母线的负荷无功量，V_k是S_NZ500中母线的电压幅值，θ_k是S_NZ500中母线的电压相角，P_d是S_NZ500中的直流输电有功功率；

(4-3-2)从步骤(2)和步骤(3)得到的总支持向量回归模型中选择适用于当前运行方式的子支持向量回归模型，分别计算当前运行方式与96个设定运行方式之间的距离从96个中选择与的最小值相对应的子支持向量回归模型，作为适用于当前运行方式的子支持向量回归模型；

(4-3-3)根据上述步骤(4-3-2)选择的子支持向量回归模型，以直流逆变换流站近区的暂态电压稳定特征量S_NZ,F作为模型输入，计算得到直流逆变换流站近区的动态无功储备需求S_Q，实现直流受端电网动态无功储备需求的评估。