1.一种空间负荷预测中确定元胞负荷合理最大值的CEEMD方法，其特征在于，它包括以下步骤：

1)基础资料整合

将待预测区域的网架结构、市政信息、用地信息和历史负荷数据相关的基础资料进行归纳和整合；

2)建立电力地理信息系统

首先配准基图，然后在待预测区域内创建用地信息图层、I类元胞图层和II类元胞图层，并建立与之对应的数据库，进而形成电力地理信息系统；

3)合理最大值的确定

采用CEEMD对I类元胞负荷时间序列进行分解，将分解得到的各本征模态分量，本征模态分量的英文全称是Intrinsic Mode Function，英文缩写是IMF，通过判据筛选出其中的高频分量并予以剔除，将余下IMF分量重构为主体分量，所述主体分量中的最大值即为I类元胞负荷合理最大值；

①CEEMD分解各元胞负荷时间序列

首先向原始信号中加入g组辅助白噪声，每组噪声均由正噪声和负噪声组成，从而生成

2组混合的时间序列M1和M2；

其中，S为原始数据；N为添加的白噪声；M1为原始信号与正噪声的叠加；M2为原始信号与负噪声的叠加，从而得到2g个信号；

然后对M1和M2进行EMD分解，从而得到2组IMF的集合，分别包含正噪声和负噪声；

最后将包含正噪声和负噪声的每对IMF分量整合为最终的IMF分量：其中，Cij为第j个信号的第i个IMF分量；Ci为待处理信号经过CEEMD后得到最终的第i个IMF分量；j＝1,2,…,2g，g为信号的个数；i＝1,2,…,v，v为元胞负荷时间序列经CEEMD后得到的IMF分量的个数；

②建立基于游程检验的高频分量识别判据

一方面，电力负荷中的高频分量具有随机波动性，即高频分量的时间序列中的数据与顺序无关；另一方面，游程检验技术能够判断观察值的顺序是否随机；所以采用游程检验技术识别电力负荷中的高频分量；

采用所述游程检验技术对各I类元胞历史负荷数据进行CEEMD，分别将每个IMF分量的所有数据与各自的中位数进行比较，用1替换该IMF分量中大于中位数的数据，用0替换该IMF分量中小于中位数的数据，等于中位数的数据用其前一个数据的替换值来替换；若第一个数据等于中位数，则用其后第一个不等于中位数的数据的替换值来替换，得到对应IMF分量的只含有0和1的时间序列，一个游程是指一段没有间断的相同数序列，该时间序列中0的总个数记为n1，1的总个数记为n2，当n1>20或n2>20时，游程数R的分布近似为正态分布，则：其中，n1为0的个数；n2为1的个数；R为游程数；E(R)、D(R)分别为游程数的数学期望、方差；Rα为临界值；Zα为正态分布值；α为显著性水平；

α是发生“弃真错误”的最大概率，不宜取得过大，不超过0.1；但也不是α越小越好，因为在样本容量取定的情况下，α的减小，尽管保证发生“弃真错误”的可能性变小，然而同时使得发生“取伪错误”的概率增大；所以，α的值为0.05，对应有Zα＝±1.96；

建立反映I类元胞负荷随机波动的高频分量识别判据，计算公式为式(6)；

其中，Rα.min为游程数的临界值下限；RIMFi为第i个IMF分量的游程数，i＝1,2,…,v，v为元胞负荷时间序列经CEEMD后得到的IMF分量的个数；Rα.max为游程数的临界值上限，若某个IMF分量的游程数满足式(6)，则说明该IMF分量反映了I类元胞负荷的随机波动性，为高频分量；

③重构生成主体分量

首先采用CEEMD技术将各I类元胞负荷序列分解得到一系列的IMF分量，然后对每个IMF分量求出其各自的游程数临界值及游程数，利用式(6)的高频分量识别判据，判断出高频分量并予以剔除，剩余的IMF分量能够刻画元胞负荷的规律性和趋势性，将这些IMF分量按式(7)进行重构得到主体分量，将其中的最大值作为I类元胞负荷合理最大值；

其中，X为主体分量；i＝1,2,…,v，v为元胞负荷时间序列经CEEMD后得到的IMF分量的个数；h为高频本征模态分量的个数；

4)基于I类元胞的SLF

利用步骤3)得到历史年的I类元胞负荷合理最大值，结合灰色理论法、线性回归法、指数平滑法中的任意一种或多种传统负荷预测方法对各I类元胞进行负荷预测，均能够得到基于I类元胞的SLF结果；

5)空间电力负荷的网格化

①生成II类元胞

在根据步骤2)所建立的电力地理信息系统中，以等大小正方形网格生成II类元胞，建立II类元胞图层；

②确定负荷密度均衡系数

首先，求出每个I类元胞的负荷密度，在每个单位时间段内的I类元胞的负荷密度由式(8)求得，dik＝Pik/Si (8)

其中，dik为第i个I类元胞在第k个单位时间段内的负荷密度，i＝1,2,…,n，n为I类元胞的个数；k＝1,2,…,f，f为单位时间段的个数；Pik为第i个I类元胞在第k个单位时间段内的负荷最大值；Si为第i个I类元胞的面积；

考虑到同类负荷在发展程度不同的I类元胞中负荷密度也有差异，在此，引入负荷密度均衡系数，记作β，βik＝(dik‑dk.min)/(dk.max‑dk.min) (9)其中，βik为第i个I类元胞在第k个单位时间段内的负荷密度均衡系数；i＝1,2,…,n，n为I类元胞的个数；k＝1,2,…,f，f为单位时间段的个数；dik为第i个I类元胞在第k个单位时间段内的负荷密度；dk.max为第k个单位时间段内的所有I类元胞中负荷密度的最大值；dk.min为第k个单位时间段内的所有I类元胞中负荷密度的最小值；

③求取各分类负荷密度指标

各I类元胞内第j类用地负荷密度的最大值乘以对应的负荷密度均衡系数，即能够得到各I类元胞内第j类用地负荷密度，计算公式为式(10)，Dij＝βikDj (10)

其中，Dij为第i个I类元胞内第j类用地负荷密度，i＝1,2,…,n，n为I类元胞的个数；j＝

1,2,…,m，m为用地类型的类数；βik为第i个I类元胞在第k个单位时间段内的负荷密度均衡系数；Dj为各I类元胞内第j类用地负荷密度的最大值；

对于任意一个I类元胞，其负荷等于该I类元胞内各类用地面积与对应分类负荷密度乘积的和，计算公式为式(11)，其中，Pik为第i个I类元胞在第k个单位时间段内的负荷最大值，i＝1,2,…,n，n为I类元胞的个数；k＝1,2,…,f，f为单位时间段的个数；Sij为第i个I类元胞内第j类用地面积；j＝

1,2,…,m，m为用地类型的类数；Dij为第i个I类元胞内第j类用地负荷密度；

结合式(10)和式(11)，得式(12)，

其中，Pik为第i个I类元胞在第k个单位时间段内的负荷最大值，i＝1,2,…,n，n为I类元胞的个数；k＝1,2,…,f，f为单位时间段的个数；βik为第i个I类元胞在第k个单位时间段内的负荷密度均衡系数；Sij为第i个I类元胞内第j类用地面积，j＝1,2,…,m，m为用地类型的类数；Dj为各I类元胞内第j类用地负荷密度的最大值；

式(12)的矩阵表示形式为式(13)，

P＝BSD＝CD (13)

其中，P为I类元胞的单位时间段内负荷最大值矩阵；B为各负荷密度的均衡系数矩阵；S为I类元胞内用地面积矩阵；D为各I类元胞内分类负荷密度最大值矩阵；C为I类元胞内各类用地的面积与其对应的负荷密度均衡系数乘积的矩阵；

I类元胞的单位时间段内负荷最大估计值与分类负荷密度最大估计值之间的关系表示为式(14)，其中， 为I类元胞的单位时间段内负荷最大估计值矩阵；C为I类元胞内各类用地的面积与其对应的负荷密度均衡系数乘积的矩阵； 为各I类元胞内分类负荷密度最大估计值矩阵；

考虑到量测误差有正有负，把所有I类元胞的单位时间段内负荷最大值与估计值之差的平方和的总和记作Q，计算公式为式(15)，其中，Pi为第i个I类元胞的单位时间段内负荷最大值矩阵，i＝1,2,…,n，n为I类元胞的个数； 为第i个I类元胞的单位时间段内负荷最大估计值矩阵，i＝1,2,…,n，n为I类元胞的个数；Q为所有I类元胞的单位时间段内负荷最大值与估计值之差的平方和的总和；

利用最小二乘法原理对式(15)分类负荷密度最大估计值矩阵 进行求解，得到分类负荷密度最大估计值，计算公式为式(16)和式(17)，其中，P为I类元胞的单位时间段内负荷最大值矩阵；C为I类元胞内各类用地的面积与其对应的负荷密度均衡系数乘积的矩阵； 为各I类元胞内分类负荷密度最大估计值矩阵；

同种用地类型的平均负荷密度值获取方法式(18)，

其中， 表示第j类用地类型的平均负荷密度值； 为各I类元胞内第j类用地负荷密度最大估计值，j＝1,2,…,m,m为用地类型的类数；βik为第i个I类元胞在第k个单位时间段内的负荷密度均衡系数；Sij为第i个I类元胞内第j类用地面积；

计算得到的各种用地类型的平均负荷密度值即为各分类负荷密度指标；

在根据步骤2)所建立的电力地理信息系统中，结合根据步骤2)所建立的用地信息图层，确定出每个根据步骤5)子步骤①所生成的II类元胞内各种用地类型的面积，并将其乘以根据式(18)所求出的对应的各种用地类型的平均负荷密度值，再乘以根据步骤5)子步骤②所确定的对应的负荷密度均衡系数，计算出各II类元胞目标时间段的负荷预测值，从而实现对空间电力负荷预测结果的网格化。