1.一种基于IBBO的复杂环网方向保护MBPS求取方法，其特征在于，包括以下步骤：S1：输入电网参数，包括网络拓扑、方向保护配置；设置IBBO算法的初始参数，包括种群规模N、最大迭代次数、最大迁入概率、最大迁出概率、最大突变率、精英保留参数t％；S2：形成保护关联矩阵R，构造适应度评价函数；S3：初始化IBBO算法的种群H i ，i＝1，2，…，L；S4：初始化迭代次数T＝1；S5：处理每个种群每个个体的约束条件，对于不满足约束的个体执行约束修复策略；S6：评价种群每个个体的适应度HSI，即Habitat Suitability Index；S7：按照适应度将种群个体从优到劣进行排序；S8：判断是否达到最大迭代次数T max ，若满足转到S15，否则进入S9；S9：根据余弦迁移模型计算每个种群对应的个体数量、迁入和迁出概率；S10：对种群执行自适应迁移算子和差分变异算子操作；S11：评价新一代种群的个体适应度，按照适应度将种群从优到劣排序；S12：执行精英策略：将上一代种群L×t％个最好个体覆盖新一代种群L×t％个最差个体；S13：再次按照适应度将种群从优到劣进行排序；S14：迭代次数加1，T＝T+1，转到S5；S15：算法结束，输出最优个体对应的二进制编码，根据二进制编码与系统保护的对应关系，得到最小断点集；所述的S2中，保护关联矩阵R，定义如下：r ij 表示R中第i行第j列元素；采用如下所示适应度评价函数：其中，x i 依次对应于系统中第i个方向保护，若该保护设为断点，则x i 为1，否则为0。
2.根据权利要求1所述的一种基于IBBO的复杂环网方向保护MBPS求取方法，其特征在于：所述的S1中，种群规模N＝300，最大迭代次数＝500，最大迁入概率＝1.0，最大迁出概率＝1.0，最大突变率＝0.01，精英保留参数＝10％。
3.根据权利要求1所述的一种基于IBBO的复杂环网方向保护MBPS求取方法，其特征在于：所述的S5中，处理每个种群个体的约束条件、对于不满足约束的个体执行约束修复策略的具体步骤如下：S5.1：输入保护关联矩阵R和种群个体X；S5.2：对于X中值为1的元素，其对应的保护是断点，删除保护关联矩阵中这些保护对应的行和列；S5.3：判断保护关联矩阵R中是否存在全为零的行，若存在，则说明该行对应的保护已经从复杂环网中解环，可以计算定值，删除该保护对应的行和列；S5.4：重复S5.2，直到保护关联矩阵R中不存在元素全为零的行；若上述步骤执行结束后，R为空，则X为断点集，满足约束条件；若不为空，则不满足约束条件；对不满足约束的X执行约束修复策略：将仍无法从环网中解开的保护全部设为断点，将这些保护在X中对应的元素修正为1。
4.根据权利要求1所述的一种基于IBBO的复杂环网方向保护MBPS求取方法，其特征在于：所述的S9中，根据余弦迁移模型计算每个种群对应的迁入概率λ i 和迁出概率μ i 的公式如下：其中，I和E分别表示最大迁入、迁出概率；k i 为第i个种群的个体数量，N＝S max 为种群能够容纳的最大个体数量。
5.根据权利要求1所述的一种基于IBBO的复杂环网方向保护MBPS求取方法，其特征在于：所述的S10中，对种群执行自适应迁移算子操作的公式如下：H i (SIV)←αH i (SIV)+(1-α)H j (SIV)其中，H i (SIV)和H j (SIV)表示种群中的两个个体H i 和H j 的特征适应度向量，HSI i 和HSI j 表示这两个个体的适应度，T为迭代次数，T max 为最大迭代次数，α为(0,1)间的实数，ε为极小值，保证分母不为零；←表示将右侧根据个体H i 和H j 的特征向量计算得到的新特征向量替代个体H i 原来的特征向量，实现个体的迁移；所述的S10中，对种群执行差分变异算子操作的公式如下：H i (SIV)←H r1 (SIV)+F·(H r2 (SIV)+H r3 (SIV))其中，H r1 (SIV)、H r2 (SIV)、H r3 (SIV)为互不相同且与H i (SIV)不同的个体；F为大于0的尺度因子。
