[發明專利]一種用交替投影法修正配電網相關性模型的隨機潮流算法

在審

• 申請號: 201711296766.9
• 申請日: 2017-12-08
• 公開（公告）號: CN108123437A
• 公開（公告）日: 2018-06-05
• 發明（設計）人: 路學剛;朱杰;蔣亞坤;王珍意;趙川;趙瑩
• 申請（專利權）人: 云南電網有限責任公司
• 主分類號: H02J3/00
• 分類號: H02J3/00;H02J3/06
• 代理公司: 昆明大百科專利事務所 53106
• 代理人: 李云

1.一種用交替投影法修正配電網相關性模型的隨機潮流算法，其特征在于，所述方法包括：

步驟S1，獲得輸入變量的累積分布函數及各變量間的相關系數矩陣；

步驟S2，采用重要抽樣法抽取樣本，得到樣本矩陣；

步驟S3，基于交替投影法，通過不斷迭代找到最接近的相關系數矩陣；

步驟S4，修正最接近的相關系數矩陣，保證相關系數矩陣的對稱性與正定性；

步驟S5，對修正后的相關系數矩陣，進行Cholesky分解、相關變換和排序，確定最終的樣本矩陣；

步驟S6，將樣本矩陣帶入節點，進行潮流計算，得到節點電壓與支路功率，計算相對誤差指標。

2.如權利要求1所述的一種用交替投影法修正配電網相關性模型的隨機潮流算法，其特征在于，所述步驟S1的具體方法為：確定各節點接入的新能源類型，建立輸出功率的模型，確定輸入變量的累積分布函數，計算各變量間的相關系數矩陣。

3.如權利要求1所述的一種用交替投影法修正配電網相關性模型的隨機潮流算法，其特征在于，步驟S2所述的采用重要抽樣法抽取樣本，抽樣公式為：Z_k為輸入隨機變量的概率密度函數，A為抽樣的分層的總層數，x_k,a和y_k,a是第a層左邊界和右邊界，表示求反函數。

4.如權利要求1所述的一種用交替投影法修正配電網相關性模型的隨機潮流算法，其特征在于，所述步驟S3具體為：基于交替投影法，通過不斷迭代找到最接近的相關系數矩陣；修正原理為：

X＝P_U(P_S(P_U(…P_S(P))))→P

P_U(P)＝P-diag(diag(P-I))

P_S(P)＝Z×diag(max(v_i,0))×Z^T

P＝ZAZ^T,A＝diag(v_i)

當迭代次數接近無窮大時，輸出X趨向于P；

具體的修正步驟為：

ΔP＝0

Y＝P

X是最終修正的相關系數矩陣，Y是中間變量，I是單位矩陣；‘eig’,‘max’,‘·’,‘×’,‘diag’,‘sqrt’,‘·/’是MATLAB中的函數；E是特征向量矩陣，v是特征值的列向量；w是滿足:v＞τ₁·max(v)的元素，[E]_w是E中與w相對應的列數；τ₁,τ₂取值不同將影響誤差大小，一般τ₁＝10^-7,τ₂＝10^-4；‘||X||’是Frobenius范數，計算公式為：

X是n階矩陣。

5.如權利要求1所述的一種用交替投影法修正配電網相關性模型的隨機潮流算法，其特征在于，所述步驟S4具體為：修正最接近的相關系數矩陣，保證相關系數矩陣的對稱性與正定性；修正步驟為：

[E,v]＝eig(X)

δ＝τ₃max(v)

X＝E×diag(v)×E^T

X_V＝diag(X)

V＝sqrt(max(X_V,δ)./diag(X))

X＝(X+X^T)/2

X是最終修正的相關系數矩陣。‘eig’,‘max’,‘·’,‘×’,‘diag’,‘sqrt’,‘·/’是MATLAB中的函數；E是特征向量矩陣，v是特征值的列向量；τ₃取值不同將影響誤差大，一般τ₃＝10^-7。