技术领域
[0001] 本
发明涉及一种电
力技术，具体涉及一种基于区块链技术的配电网分布式状态感知与优化方法。
背景技术
[0002] 配电网可靠性是指直接向用户供电或分配
电能的配电网本身及其向用户供电的可靠性。在配电网可靠性评估中，其评估的分析过程可以由以下几个步骤组成：(1)对系统的运行状态选择；(2)对系统的状态进行估计，主要是进行潮流计算，看是否有越限、过负荷等情况；(3)计算各可靠性指标。
[0003] 当前，随着分布式发电技术的发展以及环境污染问题的日益严重，在靠近负荷侧引入分布式电源(以下简称DG)实现
可再生能源的就地消纳越来越受到全社会的重视。由于DG既可以独立地对负荷进行供电，也可以接入配电网，与电网一起共同对负荷供电，因此大量的DG接入后，会对配电网的网络结构和运行方式产生一定的影响：在传统的配电网中，都是由单一的电源点对负荷进行供电，是典型的
辐射状供电网络，网络中任何一条
馈线发生故障，都有可能导致接在馈线后的用户全部停电；而DG接入后，配电网变成一个由多个电源与负荷及其其它设备相联的复杂网络，当网络中某一条馈线发生故障，可能会出现由DG与主网络断开、单独对部分负荷进行供电的
孤岛运行方式。
[0004] DG的接入能有效降低配电网对上级电网的功率依赖，从而减少传统发电方式对化石
燃料的消耗。DG对配电网可靠性的影响与接入DG的
位置和容量都有关系，合理配置DG的并网位置及容量，可以有效的降低配电网的网损、改善电能
质量和提高配电网的供电可靠性。但随着配电网中DG接入数量的增加，传统的集中式优化方法对计算、存储资源的消耗很大，且过度依赖于可靠的通信网络，无法适应未来配电网的发展。
[0005] 区块链作为一种新颖的数据结构组织形式，由于去中心化、信息公开透明、数据安全可靠、历史信息可追溯等特点，已在互联网、金融等领域得到了初步的应用，近年来也逐渐被应用于电力市场领域的研究当中，不仅能完成现有的电力交易、需求响应等功能，还能在通信故障的情况下实现对配电网运行状态的分布式优化。为了进一步提高配电网的运行可靠性和经济性，需要对区域链技术加以更多的研究。
发明内容
[0006] 发明目的：针对
现有技术中传统的集中式优化方法对计算、存储资源的消耗很大，且过度依赖于可靠的通信网络的问题，本发明公开了一种基于区块链技术的配电网分布式状态感知与优化方法。
[0007] 技术方案：为实现上述技术目的，本发明采用了如下技术方案：
[0008] 一种基于区块链技术的配电网分布式状态感知与优化方法，其特征在于，包括以下步骤：
[0009] 步骤(1)、配电网优化问题模型：建立以
可再生能源的最大化利用率为目标函数的有功优化模型和以配电网有功损耗最低为目标函数的无功优化模型；
[0010] 步骤(2)、利用branch flow模型对优化问题分散化处理：忽略相
角变化对潮流计算的影响，建立单馈线配电网模型，利用分支流模型对所述步骤(1)中的优化模型进行分散化处理；
[0011] 步骤(3)、利用二阶锥规划凸松弛方法将模型凸化松弛为可求解形式：利用二阶锥规划凸松弛方法，将所述步骤(2)中经过分散化处理的非线性优化模型凸化松弛为线性优化模型；
[0012] 步骤(4)、基于交叉方向乘子法的全分布式求解：对所述步骤(3)中的线性优化模型进行全分布式求解，首先进行有功优化的分布式计算，在计算完成得到各电源有功出力后将出力值固定，再进行无功优化的分布式计算，得出最优调度方案；
[0013] 步骤(5)、基于区块链
节点的配电网分布式状态感知：在配电网的电源节点处设置包括测量和计算单元的区块链节点，构造区块链系统，测量节点之间
母线上的
电压和功率流动，并通过区块链节点之间的信息交互，估算出整个配电网的电压分布。
[0014] 作为优选，所述步骤(2)中的有功优化阶段中建立的有功模型如下：
[0015] 目标函数：以配电网内的可再生能源为最大化利用率作为优化目标，建立目标函数J1：
[0016]
[0017] 其中，Pwt、Ppv、Ppg、Pmt1、Pgt分别表示
风机、光伏、沼气发电、储能装置和微型
燃气轮机在时段t内的发电功率，Pmt2表示储能装置在时段t内从电网吸收的功率；
[0018] 约束条件包括：
[0019] a.功率平衡约束：
[0020]
[0021] 其中，PDGi、PMGi、Pgrd、Pload、Ploss分别表示分布式电源出力值、微网出力值、配电网从上级电网的购电量、负荷功率和系统有功损耗，NDG、NMG分别表示分布式电源和微网的数量；
[0022] b.各类电源有功出力约束：
[0023] 配电网中光伏、风机、沼气发电、储能装置和微型燃气轮机在各调度时段上的有功出力值都应当保持在允许的输出范围内：
[0024] Ppv，min≤Ppv≤Ppv，max
[0025] Pwt，min≤Pwt≤Pwt，max
[0026] Ppg，min≤Ppg≤Ppg，max
[0027] Pmt1，min≤Pmt1≤Pmt1，max
[0028] Pgt，min≤Pgt≤Pgt，max
[0029] c.微型燃气轮机爬坡率约束：
[0030] Pm，t-Pm，t-1≤URm
[0031] Pm，t-1-Pm，t≤DRm
[0032] 其中，URm、DRm分别为第m台微型燃气轮机的最大上升和最大下降速率；
[0033] d.储能装置约束
[0034] 为了保证优化调度的
可持续性，储能装置在完成一次日前调度后，应当回到初始的储能状态，也就是在调度周期内保持充放电的平衡：
[0035] ΔE＝Ein-Eout＝0
[0036] 其中，Ein、Eout分别表示调度周期内储能装置吸收和放出的电量；
[0037] 此外，储能装置的过充过放会影响其使用寿命，因此需要将储能的
荷电状态限制在合理的区间内：
[0038] Smin≤St≤Smax
[0039] 其中，St表示储能装置的荷电状态，Smin、Smax分别表示储能装置荷电状态的上下限；
[0040] e.系统安全约束
[0041] 分布式电源的接入改变了配电网中的潮流分布，在调整分布式电源出力值的同时，线路上流过的
电流也随之变化，因此必须将电流值限制在允许的范围内：
[0042] Iij≤Imax
[0043] 其中，Iij表示线路i、j上流过的电流，Imax表示线路上允许流过电流的最大值。
[0044] 作为优选，所述步骤(2)中的无功优化阶段中建立的无功模型如下：
[0045] 目标函数：以配电网有功损耗最低作为优化目标，建立目标函数J2：
[0046]
[0047] 其中，Pi为系统中各节点的注入的有功功率；
[0048] 约束条件包括：
[0049] a.潮流平衡约束
[0050]
[0051]
[0052] 其中，Pi、Qi、Vi表示节点i处输入的有功功率、
无功功率和节点电压，Gij、Bij、δij表示节点i、j之间的电导、电纳、电压与电流之间的相角差；
[0053] b.节点电压约束
[0054] Vi，min＜Vi＜Vi，max
[0055] 式中，Vi，min、Vi，max分别表示节点i电压幅值允许的上下限；
[0056] c.电源无功补偿约束
[0057] Qi，min≤Qi≤Qi，max
[0058] 其中，Qi表示电源i输出的无功功率，Qi，min、Qi，max分别表示电源i输出的无功功率的上下限。
[0059] 作为优选，所述步骤(2)中分散化处理后的优化模型中的潮流平衡约束为：
[0060]
[0061] 其中，lij＝Iij2，vi＝Vi2，vj＝Vj2，ri、xi表示线路i、j上的
电阻和电抗；pi、qi分别表示i节点的有功和无功负荷；Pi、Qi表示线路i、j上流过的有功和无功功率。
[0062] 作为优选，所述步骤(3)中凸化松弛后的线性优化模型为：
[0063]
[0064] 其中，lij＝Iij2， Pij表示线路i～j上流过的有功功率，并固定vi＝Vi2。
[0065] 作为优选，所述步骤(3)中，若凸松弛后的优化模型与原模型具有相同的最优解，该凸优化模型为适用于交叉方向乘子法求解的分散式配电网优化模型。
[0066] 作为优选，所述电源节点为
变压器节点、分布式电源节点和微网节点。
[0067] 作为优选，所述步骤(5)中，各区块链节点将配电网分段，每个节点负责其下游区段的状态感知，节点装设测量装置用于测量本地母线上的电压和功率流动，节点两两之间可以通信；在一条馈线系统中，选取节点1、节点3为关键节点，当P1＞0、P2＞0、Q1＞0、Q2＞0时，电压最低值为U3；当P1＜0、P2＜0、Q1＜0、Q2＜0时，电压最低值为U1；在除去上述两种情况后，区域内的最
低电压点才有可能位于中间的某个位置，此时通过节点1、节点3分别估算区域的最低电压值；
[0068] 由节点1估算区域的最低电压值：
[0069]
[0070] 由节点3估算区域的最低电压值：
[0071]
[0072] 取二者的平均值作为区域最低电压的估算值：
[0073]
[0074] 关键节点将遥测数据传递给具备计算能力的区块链节点，区块链节点利用上述方法计算出区域内电压的最高值和最低值，最后通过区块链节点之间的信息交互，估算出整个配电网大致的电压分布。
[0075] 有益效果：由于采用了上述技术方案，本发明具有如下技术效果：
[0076] 1.本发明提出了一种基于区块链技术的配电网分布式状态感知与优化方法，通过在所有的电源节点装设测量、计算、通信装置来建立区块链系统，利用区块链节点的测量数据，以估算的方式实现对配电网运行状态的分布式感知，提高了弱通信情况下的优化运行能力，使得区块链系统不仅能完成现有的电力交易、需求响应等功能，还能在通信故障的情况下实现对配电网运行状态的分布式优化，有效提高了配电网的运行可靠性和经济性。
[0077] 2.利用二阶锥优化凸松弛方法实现了集中式优化问题的分布式求解，使得问题适用于分散式优化方法；
[0078] 3.引入交叉方向乘子法后，解决了分散式优化模型的求解问题。
附图说明
[0080] 图2是本发明的基于区块链节点的配电网状态感知的原理图；
[0081] 图3是本发明的单馈线配电网示意图；
[0082] 图4是本发明的配电网分支流模型；
[0083] 图5是本发明的区块链
框架下的配电网分布式优化流程。
具体实施方式
[0084] 下面结合附图对本发明作更进一步的说明。
[0085] 一种基于区块链技术的配电网分布式状态感知与优化方法，如图1所示，包括以下步骤：
[0086] 步骤(1)、在配电网上建立优化模型，包括以可再生能源的最大化利用率为目标函数的有功优化模型和以配电网有功损耗最低为目标函数的无功优化模型。
[0087] 1)、有功优化模型如下：
[0088] 目标函数：以配电网内可再生能源风机、光伏、沼气发电、储能的最大化利用率作为优化目标
[0089]
[0090] 其中，Pwt、Ppv、Ppg、Pmt1、Pgt分别表示风机、光伏、沼气发电、储能装置和微型燃气轮机在时段t内的发电功率，Pmt2表示储能装置在时段t内从电网吸收的功率；
[0091] 约束条件：
[0092] a.功率平衡约束
[0093]
[0094] 其中，PDGi、PMGi、Pgrd、Pload、Ploss分别表示分布式电源出力值、微网出力值、配电网从上级电网的购电量、负荷功率和系统有功损耗，NDG、NMG分别表示分布式电源和微网的数量；
[0095] b.各类电源有功出力约束
[0096] 配电网中光伏、风机、沼气发电、储能装置和微型燃气轮机在各调度时段上的有功出力值都应当保持在允许的输出范围内：
[0097]
[0098] c.微型燃气轮机爬坡率约束
[0099] Pm，t-Pm，t-1≤URm
[0100] Pm，t-1-Pm，t≤DRm (7)
[0101] 其中，URm、DRm分别为第m台微型燃气轮机的最大上升和最大下降速率；
[0102] d.储能装置约束
[0103] 为了保证优化调度的可持续性，储能装置在完成一次日前调度后，应当回到初始的储能状态，也就是在调度周期内保持充放电的平衡：
[0104] ΔE＝Ein-Eout＝0 (8)
[0105] 其中，Ein、Eout分别表示调度周期内储能装置吸收和放出的电量；
[0106] 此外，储能装置的过充过放会影响其使用寿命，因此需要将储能的荷电状态限制在合理的区间内：
[0107] Smin≤St≤Smax (9)
[0108] 其中，St表示储能装置的荷电状态，Smin、Smax分别表示储能装置荷电状态的上下限；
[0109] e.系统安全约束
[0110] 分布式电源的接入改变了配电网中的潮流分布，在调整分布式电源出力值的同时，线路上流过的电流也随之变化，因此必须将电流值限制在允许的范围内：
[0111] Iij≤Imax (10)
[0112] 其中，Iij表示线路i、j上流过的电流，Imax表示线路上允许流过电流的最大值。
[0113] 2)、无功优化模型如下：
[0114] 目标函数：以配电网有功损耗最低作为优化目标
[0115]
[0116] 其中，Pi为系统中各节点的注入的有功功率；
[0117] 约束条件：
[0118] a.潮流平衡约束
[0119]
[0120] 其中，Pi、Qi、Vi表示节点i处输入的有功功率、无功功率和节点电压，Gij、Bij、δij表示节点i、j之间的电导、电纳、电压与电流之间的相角差；
[0121] b.节点电压约束
[0122] Vi，min＜Vi＜Vi，max (13)
[0123] 式中，Vi，min、Vi，max分别表示节点i电压幅值允许的上下限；
[0124] c.电源无功补偿约束
[0125] Qi，min≤Qi≤Qi，max (14)
[0126] 其中，Qi表示电源i输出的无功功率，Qi，min、Qi，max分别表示电源i输出的无功功率的上下限。
[0127] 步骤(2)、忽略相角变化对潮流计算的影响，建立单馈线配电网模型，利用分支流模型对所述步骤(1)中的优化模型进行分散化处理。
[0128] 由于配电网通常为辐射网，馈线上电压、电流的相角变化是很小的，忽略相角变化对于潮流计算而言是足够精确的，因此采用图3所示的模型用于描述单条馈线的潮流情况。利用分支流模型表示的馈线潮流平衡方程：
[0129]
[0130] 其中，lij＝Iij2，vi＝Vi2，vj＝Vj2，ri、xi表示线路i、j上的电阻和电抗；pi、qi分别表示i节点的有功和无功负荷；Pi、Qi表示线路i、j上流过的有功和无功功率。
[0131] 可见在上式表示的潮流平衡方程已经解除了节点之间的耦合关系，同时忽略了电压、电流的相角。
[0132] 步骤(3)、利用二阶锥规划凸松弛方法，将所述步骤C中经过分散化处理的非线性优化模型凸化松弛为线性优化模型。
[0133] 相比于传统的集中式模型，分散式电气模型解除了设备间的耦合关系，极大地提高了计算速度，能够实现优化问题的快速求解。由于潮流平衡约束的非线性特点，模型为非凸模型，仍不具备利用交叉方向乘子法求解的条件。利用凸松弛方法可将非线性约束松弛成为线性约束，从而将模型转
化成凸模型进行求解。众多的凸松弛方法中二阶锥凸松弛法是较为常见的凸松弛方法，其所建模型相对其它方法复杂度较低，非常适合应用与辐射型网络。根据二阶锥凸松弛方法，可将模型转换为凸优化模型：
[0134]
[0135] 其中，lij＝Iij2， Pij表示线路i～j上流过的有功功率，并固定vi＝Vi2。
[0136] 若凸松弛后的优化模型与原模型具有相同的最优解，称其该凸优化模型为精确的。对于辐射型配电网，只需满足微弱的条件，凸化后的模型便具有良好的精确性，因此不妨认为模型是精确的。至此便得到了适用于交叉方向乘子法求解的分散式配电网优化模型。
[0137] 步骤(4)、基于交叉方向乘子法对所述步骤(3)中的线性优化模型进行全分布式求解。本发明采用的优化运行模型分为有功优化和无功优化两个阶段，两阶段的优化任务均在区块链节点上利用交叉方向乘子法完成求解。各节点首先以配电网日运行成本最低为目标函数进行有功优化的分布式计算，在计算完成得到各电源有功出力后将出力值固定，进行以系统网损最低为目标函数的无功优化得出最优调度方案。在调度方案的
基础上计算交易信息并进行记录用于后期清算，具体流程如图5所示。
[0138] 由于在配电网中，普通负荷节点不具备对本地母线电气量的测量能力，只能依赖节点间的信息传递，通过逐点计算来获得系统的运行状态，这种方法在节点通讯故障时将无法发挥作用。因此若能够利用区块链节点，以及少数具有遥测能力的关键节点，独立感知区域内负荷节点的运行状态，对保证配电网在通信故障情况下的安全运行具有重要意义。本发明基于区块链技术的配电网分布式状态感知与优化方法，在配电网的变压器节点、分布式电源节点和微网(以下简称MG)节点处设置区块链节点，即通过在配电网的电源节点(变压器、DG、MG节点)配置计算单元和通信、测量模块构建区块链系统，测量节点之间母线上的电压和功率流动，并通过区块链节点之间的信息交互，估算出整个配电网的电压分布。
[0139] 如图2所示为本发明的基于区块链节点的配电网状态感知的原理图。
[0140] 根据区块链节点将配电网分段，每个节点负责其下游区段(图2虚线框)的状态感知，节点装设测量装置用于测量本地母线上的电压和功率流动，节点两两之间可以通信。将配电网中的某一段馈线简化表示成图3单馈线配电网示意图。图中所示的简单馈线系统中，节点1、节点3为关键节点。显然，当P1＞0、P2＞0、Q1＞0、Q2＞0时，电压最低值为V3；当P1＜0、P2＜0、Q1＜0、Q2＜0时，电压最低值为V1；在除去上述两种情况后，区域内的最低电压点才有可能位于中间的某个位置，此时通过节点1、节点3分别估算区域的最低电压值。
[0141] 由节点1估算区域的最低电压值：
[0142]
[0143] 由节点3估算区域的最低电压值：
[0144]
[0145] 为了保留一定的裕度，取二者的平均值作为区域最低电压的估算值：
[0146]
[0147] 关键节点将遥测数据传递给具备计算能力的区块链节点，后者利用上述方法计算出区域内电压的最高值和最低值。最后通过区块链节点之间的信息交互，估算出整个配电网大致的电压分布。
[0148] 本发明并不限于上文描述的实施方式，以上所述只仅是本发明的较佳
实施例，并不用以限制本发明。凡是依据本发明的技术实质所作的任何
修改、等同变换、改进等，均属于本发明所要求的保护范围。
