本发明提出了一种基于暂态电压灵敏度的配电网调压控制方法。对分布式光伏密集接入下的直流受端近区电网，采样直流闭锁后各节点的电压直至趋于稳态，得到多组节点电压向量并分组；根据分组计算各节点的暂态电压裕度，通过与裕度阈值进行比较得到电压薄弱点；依次扰动无功补偿设备的无功，计算各无功补偿设备的平均暂态电压灵敏度并排序，根据排序结果得到无功补偿设备调节顺序；根据调节顺序，结合无功补偿设备的约束条件依次调节各无功补偿设备，直至各电压薄弱节点的暂态电压稳定裕度均大于裕度阈值。本发明协调各无功补偿设备参与特高压直流受端近区电网调压，从而提高电压薄弱点暂态电压稳定性。
1.一种基于暂态电压灵敏度的配电网调压控制方法，其特征在于，包括以下步骤：步骤1，将多个分布式光伏电源、多个无功补偿设备接入特高压直流受端近区交流配电网中，模拟直流闭锁后，对特高压直流受端近区交流配电网的节点电压进行采样直至特高压直流受端近区交流配电网趋于稳态，得到多组特高压直流受端近区交流配电网中节点电压向量，进一步构建特高压直流受端近区交流配电网中节点电压分组；
步骤2：根据特高压直流受端近区交流配电网中节点电压分组，计算特高压直流受端近区交流配电网中节点的暂态电压裕度，进一步通过与裕度阈值进行比较，得到特高压直流受端近区交流配电网中电压薄弱节点；
步骤3：依次对无功补偿设备进行无功扰动，计算特高压直流受端近区交流配电网中电压薄弱节点的无功补偿设备无功扰动后暂态电压灵敏度，计算无功补偿设备无功扰动后平均暂态电压灵敏度，将无功补偿设备无功扰动后平均暂态电压灵敏度进行从大至小排序得到排序后无功补偿设备无功扰动后暂态电压灵敏度，根据排序后无功补偿设备无功扰动后暂态电压灵敏度得到无功补偿设备调节顺序；
步骤4：根据无功补偿设备的条件顺序，结合无功补偿设备的约束条件依次调节各无功补偿设备；
步骤5：重复执行步骤4直至特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压裕度均大于裕度阈值；
步骤1所述特高压直流受端近区交流配电网中节点的数量为：N；
步骤1所述分布式光伏电源的数量为：Mpv，即特高压直流受端近区交流配电网中所述分布式光伏电源接入节点的数量为Mpv；
步骤1所述无功补偿设备的数量为：Msvg，即特高压直流受端近区交流配电网中所述无功补偿设备接入节点的数量为Msvg；
步骤1所述特高压直流受端近区交流配电网中节点电压向量为：(u1(m),u2(m),...,uN(m))，m∈[1,M]其中，(u1(m),u2(m),...,uN(m))表示第m次采样的特高压直流受端近区交流配电网中节点电压向量，ui(m)表示第m次采样的特高压直流受端近区交流配电网中第i个节点的电压，i∈[1,N]，N为特高压直流受端近区交流配电网中节点的数量，M为特高压直流受端近区交流配电网的节点电压直至特高压直流受端近区交流配电网趋于稳态的采样次数；
步骤1所述进一步构建特高压直流受端近区交流配电网中节点电压分组：其中，特高压直流受端近区交流配电网中第i个节点的电压的第n个分组为：(ui(n),ui(n+1),...,ui(n+τcr/T‑1))i∈[1,N],n∈[1,M‑τcr/T+1]其中，特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中第k个电压为：
其中，τcr为电压二元表的时间参数，T为采样的间隔时间，N为特高压直流受端近区交流配电网中节点的数量，M为特高压直流受端近区交流配电网的节点电压直至特高压直流受端近区交流配电网趋于稳态的采样次数，M‑τcr/T+1为特高压直流受端近区交流配电网中第i个节点的电压的分组数量，τcr/T特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中电压的数量；
步骤2所述计算特高压直流受端近区交流配电网中节点的暂态电压裕度为：ξi＝min{ξi,1,ξi,2,...,ξi,τcr/T}i∈[1,N],n∈[1,M‑τcr/T+1],k∈[1,τcr/T]其中，ξi,n为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组的暂态电压裕度，ξi为特高压直流受端近区交流配电网中第i个节点的暂态电压裕度，ucr为电压二元表的电压参数，τcr为电压二元表的时间参数，uN为额定电压，T为采样的间隔时间，ui(n+k‑1)为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中第k个电压，ui(n+k)为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中第k+1个电压，N为特高压直流受端近区交流配电网中节点的数量，M为特高压直流受端近区交流配电网的节点电压直至特高压直流受端近区交流配电网趋于稳态的采样次数，M‑τcr/T+1为特高压直流受端近区交流配电网中第i个节点的电压的分组数量，τcr/T特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中电压的数量；
步骤2所述进一步通过与裕度阈值进行比较为：ξ1,ξ2,...,ξN依次与裕度阈值进行比较，若小于等于裕度阈值则被判定为电压薄弱节点；
步骤2所述的特高压直流受端近区交流配电网中电压薄弱节点为：(Weak1,Weak2,...,WeakK)其中，K为特高压直流受端近区交流配电网中的电压薄弱节点的数量，Weakj为特高压直流受端近区交流配电网中第j个薄弱节点，即特高压直流受端近区交流配电网中序号为Weakj的节点，j∈[1,K]；
特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压裕度为：
2.根据权利要求1所述的基于暂态电压灵敏度的配电网调压控制方法，其特征在于：步骤3所述无功补偿设备包括：
步骤1所述的多个分布式光伏电源为Mpv台分布式光伏电源、步骤1所述的多个无功补偿设备为Msvg台无功补偿设备；
第l台无功补偿设备的无功由Ql注入ΔQl后，扰动后第l台无功补偿设备的无功为Ql+ΔQl；
步骤3中所述计算特高压直流受端近区交流配电网中电压薄弱节点的无功补偿设备无功扰动后暂态电压灵敏度为：
其中，λl,j为特高压直流受端近区交流配电网中第j个电压薄弱节点的第l台无功补偿设备无功扰动后暂态电压灵敏度，Ql+ΔQl为扰动后第l台无功补偿设备的无功，Ql为第l台无功补偿设备的无功，ξj(Ql+ΔQl)为步骤2所述特高压直流受端近区交流配电网中第j个薄弱节点即第Weakj节点的第l台无功补偿设备无功扰动后暂态电压裕度，ξj(Ql)为步骤2所述特高压直流受端近区交流配电网中第j个薄弱节点即第Weakj节点的暂态电压裕度，Mpv+Msvg为无功补偿设备的数量，K为特高压直流受端近区交流配电网中的电压薄弱节点的数量；
步骤3所述计算无功补偿设备无功扰动后平均暂态电压灵敏度为：l∈[1,Mpv+Msvg]
其中，λl为第l台无功补偿设备无功扰动后平均暂态电压灵敏度，Mpv+Msvg为无功补偿设备的数量；
步骤3所述将无功补偿设备无功扰动后平均暂态电压灵敏度进行从大至小排序为：将 从大至小排序，得到步骤3所述排序后无功补偿设备无功扰动后暂态电压灵敏度为：
其中， 为排序后第z个位置上无功补偿设备无功扰动后暂态电压灵敏度， 等效于第dgz台无功补偿设备无功扰动后平均暂态电压灵敏度，步骤3所述无功补偿设备调节顺序为：由于灵敏度计算量大，需要计算各无功补偿设备无功摄动时的暂态响应计算，在工程实践中，可通过并行运算来加快优化计算速度。
3.根据权利要求2所述的基于暂态电压灵敏度的配电网调压控制方法，其特征在于：步骤4所述根据无功补偿设备的调节顺序为：根据步骤3所述无功补偿设备的调节顺序步骤4所述结合无功补偿设备的约束条件依次调节各无功补偿设备为：若第dgz台无功补偿设备为分布式光伏系统，其约束条件为：z∈[1,Mpv+Msvg]
dgz∈[1,Mpv+Msvg]其中， 为第dgz台无功补偿设备的分布式光伏系统无功输出， 为第dgz台无功补偿设备的分布式光伏系统无功调节最大值；
根据第dgz台无功补偿设备的约束条件，每次调节增加量均为Q ，并通过步骤2得到特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压裕度即若 均大于裕度阈值则调节停止，否则继续参照step
若 达到 中仍存在小于等于裕度阈值的情况，则通过第dgz+1台无功补偿设备继续调节；
若第dgz台无功补偿设备为无功补偿设备，其约束条件为：z∈[1,Mpv+Msvg]
dgz∈[1,Mpv+Msvg]其中， 为第dgz台无功补偿设备的无功补偿设备无功输出， 为第dgz台无功补偿设备的无功补偿设备无功调节最大值；
根据第dgz台无功补偿设备的约束条件，每次调节增加量均为Q ，并通过步骤2得到特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压裕度即若 均大于裕度阈值则调节停止，否则继续参照step
若 达到 中仍存在小于等于裕度阈值的情况，则通过第dgz+1台无功补偿设备继续调节。
一种基于暂态电压灵敏度的配电网调压控制方法
技术领域
[0001] 本发明属于特高压直流受端近区电网调压控制技术领域，尤其涉及一种基于暂态电压灵敏度的配电网调压控制方法。
背景技术
[0002] 我国的能源中心和负荷中心呈现严重的逆向分布。西北地区的优质清洁能源和西南地区的水电通过多条特高压直流输电线路输送到位于华东及华中地区的负荷中心，这既
产生了巨大的社会、经济和环境效益，同时也深刻改变了现代电网的结构，使其运行与控制
具有新的特征。而在作为负荷中心和直流落点的安徽省，新能源装机呈爆发式增长，近三年
内光伏规模增长8倍，风电规模增长近1倍。但新能源并网电压等级低，分布式特点显著，全
省110kV及以下电压等级并网的容量超过全省总容量的95％，且分散在35kV、10kV及380V等
各个电压等级电网中。2018年底，吉泉特高压直流落地安徽古泉换流站。针对特高压直流受
端近区电网的暂态电压稳定问题，常用的解决办法是加装调相、SVG、STATCOM等动态无功补
偿装置，而依靠分布式光伏提供无功支撑也是一种极为有效的解决途径。图1为分布式光伏
密集接入的特高压直流受端近区电网系统结构图。对于类似分布式光伏密集接入的特高压
直流受端近区电网而言，直流受端系统外来输电比例大幅度提高，且换流站需要消耗大量
的无功功率。由于电力电子装置时间常数极小，动态响应速度非常快，故系统的等效惯性时
间常数减小，直流受端电压调节能力减弱。并且由于分布式光伏的耐压能力低，特高压直流
闭锁后的系统电压波动，可能造成分布式光伏大规模无序脱网，进一步恶化系统的电压，受
端近区电网电压稳定面临更大的故障冲击威胁。综上所述，对于类似分布式光伏密集接入
的特高压直流受端近区电网，暂态电压问题值得得到更多关注，需要提出更加有效的控制
方法来保证特高压直流受端近区电网中的分布式光伏不会因为电压问题而脱网。
[0003] 传统的光伏有功控制系统功率因数通常不具备常规机组的无功支撑能力。为了能够充分挖掘分布式光伏密集接入的特高压直流受端近区电网的调压能力，需要研究分布式
光伏参与特高压直流受端近区电网调压的协调优化方法。
[0004] 每个光伏电站和无功补偿设备的当前输出无功均能通过WAMS实施监测。在固定的故障校验水平(直流闭锁)，便能够找到当前运行点下的电压薄弱点，匹配出电压薄弱点附
近各光伏电站和无功补偿设备的暂态电压灵敏度与分析出当前的约束条件，从而优化计算
出各光伏电站和无功补偿设备的无功补偿量。通过各光伏电站和无功补偿设备来提供无功
支撑，从而提高系统电压稳定性。
发明内容
[0005] 本发明的目的是提供一种基于暂态电压灵敏度的配电网调压控制方法。
[0006] 本发明通过建立分布式光伏密集接入下的特高压直流受端近区电网等效仿真模型，获取直流闭锁后特高压直流受端近区电网的电压波动曲线；计算各节点的暂态电压稳
定裕度，将每个节点的暂态电压稳定裕度ξu,i与裕度阈值比较，识别配电网模型中的电压薄
弱节点，再依次扰动包含分布式光伏在内的各无功补偿设备，计算其对各电压薄弱节点的
暂态电压灵敏度，然后将第k个无功补偿设备对各电压薄弱点的暂态电压灵敏度进行加权
求和得到平均暂态电压灵敏度，根据平均暂态电压灵敏度对各无功补偿设备进行排序，结
合构建的约束条件模型，按照灵敏度大小依次调节各无功补偿设备，直至各电压薄弱节点
的暂态电压稳定裕度均大于裕度阈值。本发明实现了协调包括分布式光伏在内的各无功补
偿设备参与特高压直流受端近区电网调压，改善了特高压直流受端近区电网直流闭锁扰动
下的暂态电压稳定性。
[0007] 本发明采用如下的技术方案：
[0008] 一种基于暂态电压灵敏度的配电网调压控制方法，包括步骤：
[0009] 步骤1，将多个分布式光伏电源、多个静态无功补偿设备特高压直流受端近区交流配电网中，模拟直流闭锁后，对特高压直流受端近区交流配电网的节点电压进行采样直至
特高压直流受端近区交流配电网趋于稳态，得到多组特高压直流受端近区交流配电网中节
点电压向量，进一步构建特高压直流受端近区交流配电网中节点电压分组；
[0010] 步骤2：根据特高压直流受端近区交流配电网中节点电压分组，计算特高压直流受端近区交流配电网中节点的暂态电压裕度，进一步通过与裕度阈值进行比较，得到特高压
直流受端近区交流配电网中电压薄弱节点；
[0011] 步骤3：依次对无功补偿设备进行无功扰动，计算特高压直流受端近区交流配电网中电压薄弱节点的无功补偿设备无功扰动后暂态电压灵敏度，计算无功补偿设备无功扰动
后平均暂态电压灵敏度，将无功补偿设备无功扰动后平均暂态电压灵敏度进行从大至小排
序得到排序后无功补偿设备无功扰动后暂态电压灵敏度，根据排序后无功补偿设备无功扰
动后暂态电压灵敏度得到无功补偿设备调节顺序；
[0012] 步骤4：根据无功补偿设备的条件顺序，结合无功补偿设备的约束条件依次调节各无功补偿设备；
[0013] 步骤5：重复执行步骤4直至特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度均大于裕度阈值。
[0014] 作为优选，步骤1所述特高压直流受端近区交流配电网中节点的数量为：N步骤1所述分布式光伏电源的数量为：Mpv，即特高压直流受端近区交流配电网中所述分布式光伏电
源接入节点的数量为Mpv；
[0015] 步骤1所述静态无功补偿设备的数量为：Msvg，即特高压直流受端近区交流配电网中所述静态无功补偿设备接入节点的数量为Msvg；
[0016] 步骤1所述特高压直流受端近区交流配电网中节点电压向量为：
[0017] (u1(m),u2(m),...,uN(m))m∈[1,M]
[0018] 其中，(u1(m),u2(m),...,uN(m))表示第m次采样的特高压直流受端近区交流配电网中节点电压向量，ui(m)表示第m次采样的特高压直流受端近区交流配电网中第i个节点
的电压，i∈[1,N]，N为特高压直流受端近区交流配电网中节点的数量，M为特高压直流受端
近区交流配电网的节点电压直至特高压直流受端近区交流配电网趋于稳态的采样次数；
[0019] 步骤1所述进一步构建特高压直流受端近区交流配电网中节点电压分组：
[0020]
[0021] 其中，特高压直流受端近区交流配电网中第i个节点的电压的第n个分组为：
[0022] (ui(n),ui(n+1),...,ui(n+τcr/T‑1))
[0023] i∈[1,N],n∈[1,M‑τcr/T+1]
[0024] 其中，特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中第k个电压为：
[0025] ui(n+k‑1)
[0026] k∈[1,τcr/T]
[0027] 其中，τcr为电压二元表的时间参数，T为采样的间隔时间，N为特高压直流受端近区交流配电网中节点的数量，M为特高压直流受端近区交流配电网的节点电压直至特高压直
流受端近区交流配电网趋于稳态的采样次数，M‑τcr/T+1为特高压直流受端近区交流配电网
中第i个节点的电压的分组数量，τcr/T特高压直流受端近区交流配电网中第i个节点的电压
的第n个分组中电压的数量；
[0028] 作为优选，步骤2所述计算特高压直流受端近区交流配电网中节点的暂态电压裕度为：
[0029]
[0030]
[0031] i∈[1,N],n∈[1,M‑τcr/T+1],k∈[1,τcr/T]
[0032] 其中，ξi,n为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组的暂态电压裕度，ξi为特高压直流受端近区交流配电网中第i个节点的暂态电压裕度，ucr为电
压二元表的电压参数，τcr为电压二元表的时间参数，uN为额定电压，T为采样的间隔时间，ui
(n+k‑1)为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中第k个电压，
ui(n+k)为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中第k+1个电
压，N为特高压直流受端近区交流配电网中节点的数量，M为特高压直流受端近区交流配电
网的节点电压直至特高压直流受端近区交流配电网趋于稳态的采样次数，M‑τcr/T+1为特高
压直流受端近区交流配电网中第i个节点的电压的分组数量，τcr/T特高压直流受端近区交
流配电网中第i个节点的电压的第n个分组中电压的数量；
[0033] 步骤2所述进一步通过与裕度阈值进行比较为：
[0034] ξ1,ξ2,...,ξN依次与裕度阈值进行比较，若小于等于裕度阈值则被判定为电压薄弱节点；
[0035] 步骤2所述的特高压直流受端近区交流配电网中电压薄弱节点为：
[0036] (Weak1,Weak2,...,WeakK)
[0037] 其中，K为特高压直流受端近区交流配电网中的电压薄弱节点的数量，Weakj为特高压直流受端近区交流配电网中第j个薄弱节点，即特高压直流受端近区交流配电网中序
号为Weakj的节点，j∈[1,K]；
[0038] 特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度为：
[0039]
[0040] 作为优选，步骤3所述无功补偿设备包括：
[0041] 步骤1所述的Mpv台分布式光伏电源、步骤1所述的Msvg台静态无功补偿设备；
[0042] 无功补偿设备的数量为：Mpv+Msvg；
[0043] 第l台无功补偿设备的无功为：Ql；
[0044] 第l台无功补偿设备的扰动量为：ΔQl；
[0045] 步骤3所述无功扰动为：
[0046] 第l台无功补偿设备的无功由Ql注入ΔQl后，扰动后第l台无功补偿设备的无功为Ql+ΔQl；
[0047] 步骤3中所述计算特高压直流受端近区交流配电网中电压薄弱节点的无功补偿设备无功扰动后暂态电压灵敏度为：
[0048]
[0049] 其中，λl,j为特高压直流受端近区交流配电网中第j个电压薄弱节点的第l台无功补偿设备无功扰动后暂态电压灵敏度，Ql+ΔQl为扰动后第l台无功补偿设备的无功，Ql为第
l台无功补偿设备的无功，ξj(Ql+ΔQl)为步骤2所述特高压直流受端近区交流配电网中第j
个薄弱节点即第Weakj节点的第l台无功补偿设备无功扰动后暂态电压裕度，
[0050] ξj(Ql)为步骤2所述特高压直流受端近区交流配电网中第j个薄弱节点即第Weakj节点的暂态电压裕度，Mpv+Msvg为无功补偿设备的数量，K为特高压直流受端近区交流配电网
中的电压薄弱节点的数量；
[0051] 步骤3所述计算无功补偿设备无功扰动后平均暂态电压灵敏度为：
[0052]
[0053]
[0054] l∈[1,Mpv+Msvg]
[0055] 其中，αl为第l台无功补偿设备无功扰动后平均暂态电压灵敏度，λl为第l台无功补偿设备无功扰动后平均暂态电压灵敏度，Mpv+Msvg为无功补偿设备的数量；
[0056] 步骤3所述将无功补偿设备无功扰动后平均暂态电压灵敏度进行从大至小排序为：
[0057] 将 从大至小排序，得到步骤3所述排序后无功补偿设备无功扰动后暂态电压灵敏度为：
[0058]
[0059] 其中， 为排序后第z个位置上无功补偿设备无功扰动后暂态电压灵敏度，等效于第dgz台无功补偿设备无功扰动后平均暂态电压灵敏度，
[0060]
[0061] 步骤3所述无功补偿设备调节顺序为：
[0062]
[0063] 需要进一步说明，灵敏度计算量大，此时需要计算各无功补偿设备无功摄动时的暂态响应计算，但各无功补偿设备的摄动计算本质上可并行进行，因此，在工程实践中，可
通过并行运算来加快优化计算速度。
[0064] 作为优选，步骤4所述根据无功补偿设备的调节顺序为：
[0065] 根据步骤3所述无功补偿设备的调节顺序
[0066] 步骤4所述结合无功补偿设备的约束条件依次调节各无功补偿设备为：
[0067] 若第dgz台无功补偿设备为分布式光伏系统，其约束条件为：
[0068]
[0069] z∈[1,Mpv+Msvg]
[0070] dgz∈[1,Mpv+Msvg]
[0071] 其中， 为第dgz台无功补偿设备的分布式光伏系统无功输出， 为第dgz台无功补偿设备的分布式光伏系统无功调节最大值；
[0072] 根据第dgz台无功补偿设备的约束条件，每次调节增加量均为Qstep，并通过步骤2得到特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度即
[0073] 若 均大于裕度阈值则调节停止，否则继续step
参照Q 调节；
[0074] 若 达到 中仍存在小于等于裕度阈值的情况，则通过第dgz+1台无功补偿设备继续调节；
[0075] 若第dgz台无功补偿设备为静态无功补偿设备，其约束条件为：
[0076]
[0077] z∈[1,Mpv+Msvg]
[0078] dgz∈[1,Mpv+Msvg]
[0079] 其中， 为第dgz台无功补偿设备的静态无功补偿设备无功输出， 为第dgz台无功补偿设备的静态无功补偿设备无功调节最大值；
[0080] 根据第dgz台无功补偿设备的约束条件，每次调节增加量均为Qstep，并通过步骤2得到特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度即
[0081] 若 均大于裕度阈值则调节停止，否则继step
续参照Q 调节；
[0082] 若 达到 中仍存在小于等于裕度阈值的情况，则通过第dgz+1台无功补偿设备继续调节；
[0083] 和现有技术相比，本发明具有如下优点和有益效果：
[0084] 可实时监测各机组输出无功和出力水平，匹配当下运行状态的暂态电压灵敏度，快速协调各无功补偿装置暂态无功支撑量，从而提高暂态电压稳定裕度。
[0085] 引入分布式光伏做无功补偿设备，减少了SVG等动态无功补偿设备的投资，实现基于暂态电压灵敏度的协调包含分布式光伏在内的各无功补偿设备参与特高压直流受端近
区电网调压，为提高电压薄弱点暂态电压裕度、增强含高比例可再生能源的特高压直流受
端近区电网中无功调节能力提供了新的思路。
附图说明
[0086] 图1：本发明仿真场景示意图；
[0087] 图2：本发明方法流程图；
[0088] 图3：本发明计算暂态电压裕度的示意图；
[0089] 图4：本发明直流闭锁故障时分布式光伏不参与调压模式下故障时特高压直流受端近区电网的电压响应曲线；
[0090] 图5：本发明分布式光伏系统提供无功支撑时不断迭代求解下电压薄弱点的电压响应曲线；
[0091] 图6：是直流闭锁后分布式光伏系统的无功出力曲线。
具体实施方式
[0092] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于
本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他
实施例，都属于本发明保护的范围。
[0093] 本发明具体实施设置直流双极闭锁故障场景，对某特高压直流受端近区电网进行DigSILENT仿真，考察分布式光伏系统在特高压直流受端近区电网中直流双极闭锁参与调
压的效果。DigSILENT模型中的发电机、变压器、输电线路、励磁、调速器等原件的模型和参
数均为现场实际调研参数。
[0094] 本发明特高压直流受端近区电网仿真场景示意图如图1所示，
[0095] 如图2所示为本发明的方法流程图。
[0096] 下面结合图1至图6介绍本发明的具体实施方式为一种基于暂态电压灵敏度的配电网调压控制方法，包括以下步骤：
[0097] 步骤1，将多个分布式光伏电源、多个静态无功补偿设备特高压直流受端近区交流配电网中，模拟直流闭锁后，对特高压直流受端近区交流配电网的节点电压进行采样直至
特高压直流受端近区交流配电网趋于稳态，得到多组特高压直流受端近区交流配电网中节
点电压向量，进一步构建特高压直流受端近区交流配电网中节点电压分组；
[0098] 步骤1所述特高压直流受端近区交流配电网中节点的数量为：N＝147
[0099] 步骤1所述分布式光伏电源的数量为：Mpv＝10，即特高压直流受端近区交流配电网中所述分布式光伏电源接入节点的数量为Mpv＝10；
[0100] 步骤1所述静态无功补偿设备的数量为：Msvg＝2，即特高压直流受端近区交流配电网中所述静态无功补偿设备接入节点的数量为Msvg＝2；
[0101] 步骤1所述特高压直流受端近区交流配电网中节点电压向量为：
[0102] (u1(m),u2(m),...,uN(m)) m∈[1,M]
[0103] 其中，(u1(m),u2(m),...,uN(m))表示第m次采样的特高压直流受端近区交流配电网中节点电压向量，ui(m)表示第m次采样的特高压直流受端近区交流配电网中第i个节点
的电压，i∈[1,N]，N＝147为特高压直流受端近区交流配电网中节点的数量，M＝900为特高
压直流受端近区交流配电网的节点电压直至特高压直流受端近区交流配电网趋于稳态的
采样次数；
[0104] 步骤1所述进一步构建特高压直流受端近区交流配电网中节点电压分组：
[0105]
[0106] 其中，特高压直流受端近区交流配电网中第i个节点的电压的第n个分组为：
[0107] (ui(n),ui(n+1),...,ui(n+τcr/T‑1))
[0108] i∈[1,N],n∈[1,M‑τcr/T+1]
[0109] 其中，特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中第k个电压为：
[0110] ui(n+k‑1)
[0111] k∈[1,τcr/T]
[0112] 其中，τcr＝1为电压二元表的时间参数，T＝0.01为采样的间隔时间，N＝147为特高压直流受端近区交流配电网中节点的数量，M＝900为特高压直流受端近区交流配电网的节
点电压直至特高压直流受端近区交流配电网趋于稳态的采样次数，M‑τcr/T+1＝801为特高
压直流受端近区交流配电网中第i个节点的电压的分组数量，τcr/T＝100特高压直流受端近
区交流配电网中第i个节点的电压的第n个分组中电压的数量；
[0113] 步骤2：根据特高压直流受端近区交流配电网中节点电压分组，计算特高压直流受端近区交流配电网中节点的暂态电压裕度，进一步通过与裕度阈值进行比较，得到特高压
直流受端近区交流配电网中电压薄弱节点；
[0114] 如图3所示，步骤2所述计算特高压直流受端近区交流配电网中节点的暂态电压裕度为：
[0115]
[0116]
[0117] i∈[1,N],n∈[1,M‑τcr/T+1],k∈[1,τcr/T]
[0118] 其中，ξi,n为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组的暂态电压裕度，ξi为特高压直流受端近区交流配电网中第i个节点的暂态电压裕度，ucr为电
压二元表的电压参数，τcr＝1为电压二元表的时间参数，uN为额定电压，T＝0.01为采样的间
隔时间，ui(n+k‑1)为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中
第k个电压，ui(n+k)为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中
第k+1个电压，N为特高压直流受端近区交流配电网中节点的数量，M＝900为特高压直流受
端近区交流配电网的节点电压直至特高压直流受端近区交流配电网趋于稳态的采样次数，
M‑τcr/T+1＝801为特高压直流受端近区交流配电网中第i个节点的电压的分组数量，τcr/T
＝100为特高压直流受端近区交流配电网中第i个节点的电压的第n个分组中电压的数量；
[0119] 步骤2所述进一步通过与裕度阈值进行比较为：
[0120] ξ1,ξ2,...,ξN依次与裕度阈值进行比较，若小于等于裕度阈值则被判定为电压薄弱节点；
[0121] 步骤2所述的特高压直流受端近区交流配电网中电压薄弱节点为：
[0122] (Weak1,Weak2,...,WeakK)
[0123] 其中，K为特高压直流受端近区交流配电网中的电压薄弱节点的数量，Weakj为特高压直流受端近区交流配电网中第j个薄弱节点，即特高压直流受端近区交流配电网中序
号为Weakj的节点，j∈[1,K]；
[0124] 特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度为：
[0125]
[0126] 步骤3：依次对无功补偿设备进行无功扰动，计算特高压直流受端近区交流配电网中电压薄弱节点的无功补偿设备无功扰动后暂态电压灵敏度，计算无功补偿设备无功扰动
后平均暂态电压灵敏度，将无功补偿设备无功扰动后平均暂态电压灵敏度进行从大至小排
序得到排序后无功补偿设备无功扰动后暂态电压灵敏度，根据排序后无功补偿设备无功扰
动后暂态电压灵敏度得到无功补偿设备调节顺序；
[0127] 步骤3所述无功补偿设备包括：
[0128] 步骤1所述的Mpv＝10台分布式光伏电源、步骤1所述的Msvg＝2台静态无功补偿设备；
[0129] 无功补偿设备的数量为：Mpv+Msvg＝12；
[0130] 第l台无功补偿设备的无功为：Ql；
[0131] 第l台无功补偿设备的扰动量为：ΔQl；
[0132] 步骤3所述无功扰动为：
[0133] 第l台无功补偿设备的无功由Ql注入ΔQl后，扰动后第l台无功补偿设备的无功为Ql+ΔQl；
[0134] 步骤3中所述计算特高压直流受端近区交流配电网中电压薄弱节点的无功补偿设备无功扰动后暂态电压灵敏度为：
[0135]
[0136] 其中，λl,j为特高压直流受端近区交流配电网中第j个电压薄弱节点的第l台无功补偿设备无功扰动后暂态电压灵敏度，Ql+ΔQl为扰动后第l台无功补偿设备的无功，Ql为第
l台无功补偿设备的无功，ξj(Ql+ΔQl)为步骤2所述特高压直流受端近区交流配电网中第j
个薄弱节点即第Weakj节点的第l台无功补偿设备无功扰动后暂态电压裕度，
[0137] ξj(Ql)为步骤2所述特高压直流受端近区交流配电网中第j个薄弱节点即第Weakj节点的暂态电压裕度，Mpv＝60+Msvg为无功补偿设备的数量，K为特高压直流受端近区交流
配电网中的电压薄弱节点的数量；
[0138] 步骤3所述计算无功补偿设备无功扰动后平均暂态电压灵敏度为：
[0139]
[0140]
[0141] l∈[1,Mpv+Msvg]
[0142] 其中，αl为第l台无功补偿设备无功扰动后平均暂态电压灵敏度，λl为第l台无功补偿设备无功扰动后平均暂态电压灵敏度，Mpv+Msvg＝12为无功补偿设备的数量；
[0143] 步骤3所述将无功补偿设备无功扰动后平均暂态电压灵敏度进行从大至小排序为：
[0144] 将 从大至小排序，得到步骤3所述排序后无功补偿设备无功扰动后暂态电压灵敏度为：
[0145]
[0146] 其中， 为排序后第z个位置上无功补偿设备无功扰动后暂态电压灵敏度，等效于第dgz台无功补偿设备无功扰动后平均暂态电压灵敏度，
[0147]
[0148] 步骤3所述无功补偿设备调节顺序为：
[0149]
[0150] 需要进一步说明，灵敏度计算量大，此时需要计算各无功补偿设备无功摄动时的暂态响应计算，但各无功补偿设备的摄动计算本质上可并行进行，因此，在工程实践中，可
通过并行运算来加快优化计算速度。
[0151] 步骤4：根据无功补偿设备的条件顺序，结合无功补偿设备的约束条件依次调节各无功补偿设备；
[0152] 步骤4所述根据无功补偿设备的调节顺序为：
[0153] 根据步骤3所述无功补偿设备的调节顺序
[0154] 步骤4所述结合无功补偿设备的约束条件依次调节各无功补偿设备为：
[0155] 若第dgz台无功补偿设备为分布式光伏系统，其约束条件为：
[0156]
[0157] z∈[1,Mpv+Msvg]
[0158] dgz∈[1,Mpv+Msvg]
[0159] 其中， 为第dgz台无功补偿设备的分布式光伏系统无功输出， 为第dgz台无功补偿设备的分布式光伏系统无功调节最大值；
[0160] 根据第dgz台无功补偿设备的约束条件，每次调节增加量均为Qstep，并通过步骤2得到特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度即
[0161] 若 均大于裕度阈值则调节停止，否则继step
续参照Q 调节；
[0162] 若 达到 中仍存在小于等于裕度阈值的情况，则通过第dgz+1台无功补偿设备继续调节；
[0163] 若第dgz台无功补偿设备为静态无功补偿设备，其约束条件为：
[0164]
[0165] z∈[1,Mpv+Msvg]
[0166] dgz∈[1,Mpv+Msvg]
[0167] 其中， 为第dgz台无功补偿设备的静态无功补偿设备无功输出， 为第dgz台无功补偿设备的静态无功补偿设备无功调节最大值；
[0168] 根据第dgz台无功补偿设备的约束条件，每次调节增加量均为Qstep，并通过步骤2得到特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度即
[0169] 若 均大于裕度阈值则调节停止，否则继续step
参照Q 调节；
[0170] 若 达到 中仍存在小于等于裕度阈值的情况，则通过第dgz+l台无功补偿设备继续调节；
[0171] 步骤5：重复执行步骤4直至特高压直流受端近区交流配电网中各电压薄弱节点的暂态电压稳定裕度均大于裕度阈值。
[0172] 图4为模拟直流闭锁后，采样得到的直流受端近区交流配电网各节点的暂态电压曲线，可以看到各节点能保持暂态电压稳定，但是陶楼变220kV母线的电压最低。
[0173] 图5为陶楼变220kV母线电压在迭代计算过程中的电压取线变化情况，可以看到暂态电压明显的上升了，陶楼变220kV母线的暂态电压稳定性得到了改善。
[0174] 图6为无功补偿设备提供的无功支撑，可以看到灵敏度最大的义井光伏电站首先达到了可调节的最大容量，然后是灵敏度次大的林庄光伏电站，最后皖能垃圾焚烧电站。
[0175] 表1为分布式光伏不参与调压时的仿真运行工况。故障发生前，各无功补偿设备的无功出力、无功可调节容量见表2。
[0176] 表1 仿真运行工况
[0177]
[0178] 表2 以陶楼一区为例、分布式发电不参与调压下各无功补偿设备出力及备用
[0179]
[0180] 直流闭锁后，特高压直流受端近区电网受到1000MW的功率扰动。设置在t＝1.0s发生直流闭锁故障，在t＝1.1s切除换流站无功补偿设备，特高压直流受端近区电网各地区代
表节点的电压动态响应曲线如图4所示，设系统需满足的电压二元表为(0.75pu，1s)，计算
表明，系统总体上能够保持暂态电压稳定，但陶楼一区暂态电压稳定裕度偏低，其中，220kV
母线暂态电压稳定裕度最低，为0.574。
[0181] 按照公式(2)计算出各分布式光伏的暂态电压灵敏度，其暂态灵敏度如表3所示
[0182] 表3 各无功补偿设备的暂态电压灵敏度
[0183]
[0184] 首先选取电压灵敏度最大和可调容量最大的义井光伏电站作为调节节点，设暂态电压稳定裕度下限ξu＝0.67，于是计算出义井应该补偿的容量为18.8009Mvar。使义井光伏
电站的无功输出在1.2s增加18.8009MVar，再次仿真得到暂态电压曲线，由于电力系统的非
线性，所得到的结果并不一定满足要求，所以将得到的曲线再次计算暂态电压裕度来检验
是否满足要求。经过3次迭代后暂态稳定裕度满足要求，其迭代计算过程如表4所示。
[0185] 表4 基于灵敏度分析的无功优化迭代过程
[0186]
[0187] 由表4可见，随着分布式光伏提供更多的无功支撑，而暂态稳定裕度也随之增大，陶楼一区的220kV母线的暂态电压裕度由0.5740提升到了0.6710，满足了预设要求。
[0188] 以上所述实施例仅表达了本发明的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，
在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范
围。因此，本发明专利的保护范围应以所附权利要求为准。
