[发明专利]一种电网多重时间区间同步仿真方法、装置、介质及设备

权利要求书

1.一种电网多重时间区间同步仿真方法，其特征在于，所述方法包括：

根据电网的时间变量计算所述电网的微分代数方程组，并将待仿真的时间区间划分为若干子区间；

根据所述时间区间的初始值，采用多维插值预测方法和蜂群优化跟踪预测方法预测每个子区间的初始值；

根据所述微分代数方程组和每一子区间的初始值分别对每一子区间进行单独仿真，获取每一子区间的仿真结果；

通过每一子区间构成的环形网络传递每一子区间的仿真结果，得到每一子区间对应的边界值；

根据每一子区间的初始值和对应的边界值计算得到修正量，并对每一子区间的初始值进行修正；

根据修正后的初始值，对每一子区间重新仿真，直到仿真得到的每一区间对应的边界值与仿真的初始值相同，输出仿真结果。

2.如权利要求1所述的电网多重时间区间的仿真方法，其特征在于，所述根据电网的时间变量计算所述电网仿真的微分代数方程组，并将待仿真的时间区间划分为若干子区间，具体包括：

根据电网的时间变量，计算所述电网动态仿真的微分代数方程组，所述微分代数方程组包括微分方程和代数方程；

将所述电网需要仿真的时间区间按照预设时间长度划分为若干子区间；

在区间划分时，子区间为等间隔划分，子区间内部的排列点为不等间距插入；在区间边界处，状态变量不能突变，保持连续；控制变量和代数变量能够突变，在区间边界处发生阶跃；除第一个积分子区间初始值为已知，其他积分子区间初始值均未知；

其中，所述时间变量包括状态变量x(t)、代数变量y(t)和控制变量u(t)，所述微分方程为所述代数方程为：0＝g(x(t),y(t),u(t))，f和g分别代表微分方程和代数方程的代数表达式。

3.如权利要求1所述的电网多重时间区间同步仿真方法，其特征在于，所述根据所述时间区间的初始值，采用多维插值预测方法和蜂群优化跟踪预测方法预测每个子区间的初始值，具体包括：

根据预设的相对位置在每一子区间插入若干排列点；

根据蜂群优化的跟踪预测方法计算每次的预测因子W^j+1；

将所述时间区间的初始值作为第一个子区间的初始值，对于所述若干子区间中的任一子区间，根据第一个子区间的初始值和所述时间变量的分布情况，采用张量积计算每一子区间的初始值；

其中，每次的预测因子W^j+1＝W^j+L(r)·(W^j-W^j-1)+L′(r)·(W^j-W^j-1)，L(r)为勒让德多项式，L′(r)为勒让德多项式的一阶导数；每一子区间p的初始值预测为为所述第一子区间的初始值中的状态变量的初始值，为子区间p的初始值，h为子区间的区间长度，W为子区间p的预测因子，n为子区间的数量，p为子区间的积分维度，q为子区间p中排列点的个数，i_n为第n个状态变量x(n)的插值水平，为第n个状态变量x(n)的插值集合。

4.如权利要求1所述的电网多重时间区间同步仿真方法，其特征在于，所述根据所述微分代数方程组和每一子区间的初始值分别对每一子区间进行单独仿真，获取每一子区间的仿真结果，具体包括：

计算每一子区间上的排列点的欧米茄函数系数和性质条件；

根据所述性质条件计算得到所述欧米茄函数系数的时间积分；

根据所述欧米茄函数系数、每一子区间的状态变量多项式和代数变量多项式，计算得到每一子区间的状态变量和代数变量作为仿真结果；

其中，所述欧米茄函数系数为所述性质条件为：所述欧米茄函数系数的为拉格朗日多项式的积分：拉格朗日多项式为子区间p的状态变量多项式近似为子区间p的代数变量多项式近似为x_p,r、y_p,r,和u_p分别为子区间p的第r个排列点的状态变量、代数变量以及控制变量的数值，子区间p上的排列点数量为3，分别为τ_p,1,τ_p,2和τ_p,3，δ_r,q代表克罗内克尔delta，τ_p,q为数值在[0,1]区间范围内的排列点的相对位置，微分代数方程组的控制变量u(t)在每个子区间内保持恒定。