三点估计法的交直流随机潮流计算模型(matlab版+python版)

目标问题##

3PEM的交直流随机潮流计算模型(matlab版+python版)

目标需求

1. 要适用于含多端SOP的任意大小的交直流电网,支持节点和支路的增删；
2. 要适用于接入多个服从不同概率分布的风电、光伏等分布式电源；
3. 要采用Cornish-Fisher级数拟合随机变量的概率分布函数；
4. 要考虑VSC的4种运行方式(单电压控制方式)交直流潮流分析；
5. 误差分析：与蒙特卡洛计算结果对比，误差<10^-2；
   

程序函数说明

matlab版：

1. 主程序：main4
2. 交流系统输入子程序：data_ieee33
3. 牛拉法基波子程序：NLPowerFlow文件
4. 节点导纳子程序：createYbus文件
5. 雅克比矩阵子程序：Jacobi文件(n-1+m*n-1+m)
6. 直流系统输入子程序：dc_case4文件
7. 交直流潮流子程序：ACDC_Flow文件
   列写一下文档目录:
   

程序计算步骤及流程图

以matlab版程序作为讲解。

1 参数初始化

参数初始化，读取网络参数，在此标幺化，随机潮流是在确定性基波潮流计算基础之上进行的，
基波潮流计算的学习： https://download.csdn.net/download/WConstelltion/85045068.
纯交流系统的随机潮流计算学习：https://blog.csdn.net/WConstelltion/article/details/123776191

2 蒙特卡洛计算

2.1正态分布负荷抽样：

    p_load(i,:)=normrnd(mu_load_p(i),sigma_load_p(i),1,m);%正态分布
    q_load(i,:)=normrnd(mu_load_q(i),sigma_load_q(i),1,m);

2.2DG抽样：(风速weibull分布、光强beta分布、储能功率正态分布)

     v_FD = wblrnd(FD_c,FD_k,m,1); %风速抽样
    cd_GF= betarnd(GF_a,GF_b,1,m);   %光强抽样

2.3确定风速和光强后计算DG功率：

      for i=1:m                                                            
          if v_FD(i)<vci||v_FD(i)>vco
              p_FD(i)=0;
          elseif v_FD(i)>vci&&v_FD(i)<vr
              p_FD(i)=k1*v_FD(i)+k2;
          elseif v_FD(i)>vr&&v_FD(i)<vco
              p_FD(i)=PnW;
          end
     end
                
     p_GF = cd_GF*A*yita/1000;  

2.4修正节点负荷：

    shuju.bus(index_load,3)=load_p(:,i);
    shuju.bus(index_load,4)=load_q(:,i);

2.5确定性潮流计算：

    [basemva,bus,gen,branch,loss_vsc]=ACDC_Flow(shuju, dc_case4);

2.6就可以求解节点电压、支路功率和网损的均值和标准差；

3 三点估计法计算

3.1计算估计点值和权重：(一阶,二阶)

zeta1_p=-lambda3_p./2+sqrt(lambda4_p-3*(lambda3_p.^2)/4);%ξxi,1=-(λi,3)/2+sqrt(λi,4+3/4*λi,3^2)
zeta2_p=-lambda3_p./2-sqrt(lambda4_p-3*(lambda3_p.^2)/4);%ξxi,2=-(λi,3)/2-sqrt(λi,4+3/4*λi,3^2)

3.2计算权重：

weight1=1./(zeta1_p.*(zeta1_p-zeta2_p));  %采样点1的权重pxi,1=1/(ξxi,1*(ξxi,1-ξxi,2))
weight2=-1./(zeta2_p.*(zeta1_p-zeta2_p)); %采样点2的权重pxi,2=1/(ξxi,2*(ξxi,1-ξxi,2))
weight3=1/n_load-1./(lambda4_p-lambda3_p.^2); 

3.3求解状态变量的前v阶原点矩(节点电压、支路功率以及损耗的2m+1个样本值)

    mu_vm(i)=sum(vm(i,1:2*n_load)'.*weight(1:2*n_load))+vm(i,end)*sum(weight(2*n_load+1:end));

3.4交直流潮流计算(节点电压、支路功率以及损耗的2m+1个样本值)

    [basemva,bus,gen,branch,loss_vsc]=ACDC_Flow(shuju, dc_case4);%潮流计算

3.5 CornishFish级数求解系统节点电压、直流功率以及损耗等变量的概率分布曲线

[pem_vm_pdf,pem_vm_cdf]= CornishFish(aerfa_vm,range_vm,pem_vm_mu,pem_vm_std);  %级数展开求解概率分布
[pem_xianlu_p_pdf,pem_xianlu_p_cdf]= CornishFish(aerfa_xianlu_p,range_xianlu_p,pem_xianlu_p_mu,pem_xianlu_p_std); 
[pem_ploss_pdf,pem_ploss_cdf]= CornishFish(aerfa_ploss,range_ploss,pem_ploss_mu,pem_ploss_std);  %级数展开求解概率分布

3.6 计算流程图如下图所示

基于3PEM的交直流随机潮流计算流程图

算例结果分析

1. 算例拓扑结构
   
2. 两种方法接入DG前节点电压计算结果对比图，曲线几乎吻合，误差小于3%，这是方法本身的特性造成的
   蒙特卡洛和3PEM对损耗计算结果对比
   
   对两种方法计算节点电压和支路功率的误差对比分析
   

程序连接

前辈的学习资料: https://download.csdn.net/download/lisha1002/12148929?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522164835259216780366543218%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fdownload.%2522%257D&request_id=164835259216780366543218&biz_id=1&utm_medium=distribute.pc_search_result.none-task-download-2_downloadfirst_rank_ecpm_v1~rank_v31_ecpm-5-12148929.pc_v2_rank_dl_default&utm_term=%E5%8D%8A%E4%B8%8D%E5%8F%98%E9%87%8F&spm=1018.2226.3001.4451.5.
编程不易，完整版或需要python版的私聊订制转换……