基于改进的Big Bang-Big Crunch的分布式光伏并网接纳能力计算方法
【技术领域】
[0001] 本发明属于配电系统规划技术领域，特别是涉及一种基于改进的BigBang-Big Crunch的分布式光伏并网接纳能力计算方法。
【背景技术】
[0002] 在全球性化石能源逐渐枯竭和环境问题日益恶化的背景下，光伏发电被认为是一 种具有巨大潜力的可再生能源利用形式。据全球权威能源机构预测，太阳能在本世纪中期 将成为人类能源构成中的重要组成部分，而到本世纪末将成为人类能源构成中的主要部 分。由于不可避免地会受到天气等外界环境因素的影响，光伏发电具有间歇性和波动性，因 此光伏发电大量并网将对所接入的配电系统的安全运行产生多方面的不利影响。在当前国 内光伏产业大力发展和分布式光伏持续大规模接入的新形势下，如何分析和评估配电系统 接纳光伏发电的最大潜力以及如何控制和规划分布式光伏的接入成为最近几年困扰配电 系统规划工作者的实际难题，开展对该问题的深入研究具有重要的实际意义，将为未来新 型配电系统的规划、建设和运行提供较好的理论依据和技术支撑。
[0003]目前，考虑分布式发电接入配电网所带来的有功损耗、电压分布、短路电流、电能 质量等方面的影响，针对配电网对分布式电源的接纳能力，国内外已有较多的研究成果。相 关的分析方法主要包括带约束的最优化方法、仿真模拟方法以及基于拟合的解析方法等。 最优化方法以系统最大可接入容量为目标函数，以配电系统潮流等为约束，将光伏接纳能 力问题建模为最优化模型，并进行求解；仿真模拟方法通过模拟仿真技术逐渐增加系统各 点接入容量最终得到光伏极限接纳能力；解析方法力图通过数学解析推导的方式得到并网 接纳能力与电压和电流的关系，从而获得最大接入能力的解析解。由于计算规模的限制以 及实际系统与理想系统间存在的差异性，目前仿真模拟方法和解析方法应用于实际配电系 统还存在一定的局限性。
【发明内容】
[0004] 为了解决上述问题，本发明的目的在于提供一种基于改进的BigBang-Big Crunch的分布式光伏并网接纳能力计算方法。
[0005] 为了达到上述目的，本发明提供的基于改进的BigBang-BigCrunch的分布式光 伏并网接纳能力计算方法包括按顺序执行的下列步骤：
[0006] 步骤一、基础数据获取：获取待研究的配电系统中包括网架结构、负荷水平、电气 参数在内的基础数据；
[0007] 步骤二、模型构建：利用步骤一获取的基础数据构建分布式光伏并网最大接入能 力问题的最优化模型，并确定该最优化模型的目标函数和约束条件；
[0008] 步骤三、模型转换：以上述约束条件作为罚函数，将上述最优化模型转换成求解最 小值的无约束适应度函数；
[0009] 步骤四、算法参数初始化：初始化改进的BigBang-BigCrunch算法的参数，设定 迭代次数t= 0,改进的BigBang-BigCrunch算法的碎片解个数为n，碎片解的编码长度 为配电系统中可接入分布式光伏的节点总数，并确定收缩因子a以及迭代终止条件；
[0010] 步骤五、初始碎片解生成：在无约束适应度函数的可行域内随机产生n个候选解， 即初始碎片解，形成初始碎片解集合；
[0011] 步骤六、适应度评价：根据上述碎片解对应的分布式光伏并网接入情况进行交流 潮流计算、谐波潮流计算和短路计算，并将节点电压、支路电流、电压总谐波畸变率和短路 电流结果带入无约束适应度函数的公式中，求解出对应的适应度；
[0012] 步骤七、引力解生成：根据上述碎片解的适应度，应用坍缩公式坍缩成一个候选 解，即引力解，并计算引力解的适应度：
[0013] 步骤八、爆炸迭代循环：迭代次数增加1，复制上一代最优碎片解到下一代碎片解 集合中，将上一代引力解作为参考，根据爆炸公式生成其余碎片解，形成新的碎片解集合；
[0014] 步骤九、循环终止判别：判断迭代是否完成，若未完成，则重复步骤六；若迭代结 束，则得到最终的最优解或近似最优解，同时输出相关结果，包括分布式光伏的接入位置和 容量以及最大接入能力。
[0015] 在步骤二中，最优化模型的目标函数为配电网所有可接入节点接纳的分布式光伏 的总容量，表示为：
[0016]
(1)
[0017] 其中，npv为系统中可接入分布式光伏的节点总数；PPVil为第i个接入点接入分布 式光伏的装机容量；
[0018] 约束条件包括系统潮流约束、运行电压水平约束、支路电流限制、短路电流限制、 电压谐波畸变率限制，分别表示为：
[0019]
(2):
[0020] K V A V工，随，i = 1，2,…，N (3)
[0021] K 1^_，i = 1，2,…，N，j G N(i) (4)
[0022] VTHD^;VTHD.max ⑶
[0023] IsCL^SISCL.max (6)
[0024] 其中，PvQsjlU为节点i注入的有功与无功功率；V 别为节点i和j的电 压幅值别为系统导纳矩阵的实部和虚部；0y为节点i和j的电压相角差；vliniax、 \_分别为节点电压V 上限和下限；N为节点总数，IJPI_x分别为第i条支路通过的 电流和其最大允许电流；L为支路总数；VTHD和V__分别为系统谐波电压畸变率和其最大 允许值；1^和ISa._分别为支路最大短路电流和断路器最大开断电流。
[0025] 在步骤三中，所述的求解最小值的无约束适应度函数表示为；
[0026]
C7)
[0027]其中，x;G{〇,1}，x;= 0表示该项约束未越限，x;= 1表示该项约束越限；1;为 惩罚因子；1:、12、13、14分别为运行电压水平约束、支路电流限制约束、电压谐波畸变率限制 约束、短路电流限制约束的惩罚因子。
[0028] 在步骤七中，所述的坍缩公式为：
[0029]
(8 >
[0030]其中，Xlk是第i个碎片解的第k维分量，XA是算法坍缩阶段得到的引力解的第k维分量；匕是碎片解i的适应度。
[0031] 在步骤八中，所述的爆炸公式为：
[0032]
C9)-
[0033] 其中，t是迭代次数；r是随机数，rG(-1，1) ^_和1_分别是无约束适应度函 数可行域的上限和下限。
[0034] 本发明提供的基于改进的BigBang-BigCrunch的分布式光伏并网接纳能力计算 方法从配电网规划角度对配电系统中分布式光伏的最大可接入容量问题进行建模，建立以 分布式光伏接入容量最大化为目标函数，以系统潮流约束、运行电压水平约束、支路电流限 制、电压谐波畸变率限制、短路电流限制等为约束条件的并网光伏最大可接入能力模型，通 过改进的BigBang-BigCrunch优化搜索方法来进行求解，本发明能够为光伏接入系统的 规划、建设和运营提供指导。
[0035] 本发明提供的基于改进的BigBang-BigCrunch的分布式光伏并网接纳能力计算 方法的有益效果：
[0036] (1)模型描述更全面：对分布式光伏的接入个数没有限制，可以实现一个、多个以 及全网分布式光伏的配电网接入能力确定。
[0037] (2)求解方法更有效：需要调整的参数少，只需对收缩因子a进行调整即可，对于 不同规模的配电网也不需要进行复杂的参数调整。同时，在避免局部最优解方面，改进的 BigBang-BigCrunch算法随机性很强，可以扩大搜索范围、增加种群多样性，具有更好的 抗早熟特性，更不容易陷入局部最优解，可以保证以较大的几率找到全局最优解。
【附图说明】
[0038] 图1为本发明提供的基于改进的BigBang-BigCrunch的分布式光伏并网接纳能 力计算方法流程图；
[0039] 图2为10kV33节点算例拓扑图及节点编号图；
【具体实施方式】
[0040] 下面结合附图及实施例本发明提供的基于改进的BigBang-BigCrunch的分布式 光伏并网接纳能力计算方法进行详细说明。
[0041] 下面以图2所示的33节点算例为实施例，结合图1所示的流程图对本发明提供的 基于改进的BigBang-BigCrunch的分布式光伏并网接纳能力计算方法进行详细说明。
[0042] 如图1所示，本发明提供的基于改进的BigBang-BigCrunch的分布式光伏并网 接纳能力计算方法包括顺序执行的下列步骤：
[0043] 步骤一、基础数据获取：获取待研究的配电系统中包括网架结构、负荷水平、电气 参数在内的基础数据；
[0044] 在本实施例中，获取33节点配电系统的网架结构、负荷水平、电气参数等基础数 据。其中，网架结构如图2所示，节点名称如图中所示，负荷数据如表1所示，支路电气参数 如表2所示，为便于计算，设定分布式光伏的接入位置为全部节点，功率因数为1. 0,并网分 布式光伏各次谐波电压幅值如表3所示；
[0045] 表1 33节点算例负荷参数
[0046]
[0047] 表2 33节点算例支路电气参数
[0048]
[0052] 步骤二、模型构建：利用步骤一获取的基础数据构建分布式光伏并网最大接入能 力问题的最优化模型，并确定该最优化模型的目标函数和约束条件；
[0053] 其中最优化模型的目标函数为配电网所有可接入节点接纳的分布式光伏的总容 量，表示为：
[0054]
( 1 )
[0055] 其中，npv为系统中可接入分布式光伏的节点总数；PPVil为第i个接入点接入分布 式光伏的装机容量；
[0056] 约束条件包括系统潮流约束、运行电压水平约束、支路电流限制、短路电流限制、 电压谐波畸变率限制，分别表示为：
[0057]
、2 )
[0058] 彡VAVwi= 1，2,…，N(3