【正文】 1)在電網(wǎng)規(guī)劃階段,通過(guò)潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無(wú)功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。 (2)在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長(zhǎng)及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對(duì)規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。 (3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開(kāi)機(jī)方式,有功、無(wú)功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。 　　(4)預(yù)想事故、設(shè)備退出運(yùn)行對(duì)靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案。 總結(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算。 —拉夫遜法潮流計(jì)算基本步驟牛頓—拉夫遜法潮流計(jì)算的基本步驟不外乎有如下幾步。 形成節(jié)點(diǎn)導(dǎo)納矩陣。 設(shè)各節(jié)點(diǎn)電壓的初始值 和 求出功率偏差根據(jù)公式算出雅可比矩陣。解如下線性方程組，得到依據(jù)如下公式計(jì)算修正量，并形成迭代公式：這樣反復(fù)迭代計(jì)算，直至所有節(jié)點(diǎn)|ΔU|ε和|Δδ|ε為止。計(jì)算線路功率分布，PV節(jié)點(diǎn)無(wú)功功率和平衡節(jié)點(diǎn)注入功率。第四章程序流程圖及程序代碼根據(jù)牛頓—拉夫遜法潮流計(jì)算的基本步驟可以得到如下程序流程圖%*******************主程序開(kāi)始********************clc。clear allformat long。%********************定義全局變量*****************global nSW nPQ nPV。 %平衡節(jié)點(diǎn)，PQ節(jié)點(diǎn)，PV節(jié)點(diǎn)個(gè)數(shù)global nb。 %節(jié)點(diǎn)個(gè)數(shù)global nl。 %支路數(shù) global bus。 %bus:節(jié)點(diǎn)數(shù)據(jù)global line。 %line:支路數(shù)據(jù)global Y。 %Y:節(jié)點(diǎn)導(dǎo)納矩陣global nodenum。 %節(jié)點(diǎn)編號(hào)對(duì)應(yīng)表global lPQ。 %線路損耗矩陣global myf。 %輸出文件句柄%*******************子程序調(diào)用***************openfile。 %打開(kāi)bus line數(shù)據(jù)文件 change。 %節(jié)點(diǎn)重新編號(hào) ybus。 %建立節(jié)點(diǎn)導(dǎo)納矩陣 NR。 %牛頓拉夫遜法解非線性方程序組 PQ。 %計(jì)算節(jié)點(diǎn)注入功率 flow。 %計(jì)算線路功率及損耗 ret。 %恢復(fù)原節(jié)點(diǎn)編號(hào) output。 %計(jì)算結(jié)果輸出%*************************主程序結(jié)束*****************%******************打開(kāi)bus line的數(shù)據(jù)文件*****************function openfile global nb nl。global line。global bus。global myf。 [dfile,pathname]=uigetfile(39。*.m39。,39。Select Data File39。)。 %打開(kāi)數(shù)據(jù)文件if pathname == 0 error(39。 you must select a valid data file39。)else lfile =length(dfile)。 % strip off .m eval(dfile(1:lfile2))。 %執(zhí)行數(shù)據(jù)文件中語(yǔ)句end [nl,ml]=size(line)。 %求線路數(shù)nl[nb,mb]=size(bus)。 %求節(jié)點(diǎn)數(shù)nb outfile=39。39。for I=1:length(dfile) if dfile(I)==39。.39。 break。 end outfile=strcat(outfile,dfile(I))。 endoutfile=strcat(outfile,39。.txt39。)。 %輸出文件名形成myf=fopen(outfile,39。w39。)。 %打開(kāi)輸出文件%**********************給節(jié)點(diǎn)編號(hào)*********************function changeglobal nb。global nl。global nPQ。global bus。global line。global nodenum。 nSW = 0。 nPV = 0。 nPQ = 0。 for I = 1:nb, % nb為總節(jié)點(diǎn)數(shù) type= bus(I,6)。 if type == 3, %求平衡節(jié)點(diǎn)數(shù) nSW = nSW + 1。 SW(nSW,:)=bus(I,:)。 elseif type == 2, %求PV節(jié)點(diǎn)數(shù) nPV = nPV +1。 PV(nPV,:)=bus(I,:)。else %求PQ節(jié)點(diǎn)數(shù) nPQ = nPQ + 1。 PQ(nPQ,:)=bus(I,:)。 endend bus=[PQ。PV。SW]。 %按PQ PV SW順序排列的節(jié)點(diǎn)數(shù)據(jù)矩陣newbus=[1:nb]39。nodenum=[newbus bus(:,1)]。 %形成節(jié)點(diǎn)號(hào)對(duì)應(yīng)表bus(:,1)=newbus。 for I=1:nl %根據(jù)節(jié)點(diǎn)號(hào)對(duì)應(yīng)表重新排列線路矩陣節(jié)點(diǎn)編號(hào) for J=1:2 for k=1:nb if line(I,J)==nodenum(k,2) line(I,J)=nodenum(k,1)。 break end end endend %******************建立節(jié)點(diǎn)導(dǎo)納矩陣*****************function Y = ybus(bus,line)global nl。global nb。global bus。global line。global Y。global myf。 Y=zeros(nb,nb)。 %對(duì)導(dǎo)納矩陣賦初值0 for k=1:nl I=line(k,1)。 %讀入線路參數(shù) J=line(k,2)。 Zt=line(k,3)+j*line(k,4)。 if J~=0 %接地支路分母為零 Yt=1/Zt。 %非接地支路為阻抗的倒數(shù) end Ym=line(k,5)+j*line(k,6)。 %求Ym K=line(k,7)。 %求K if (K==0)amp。(J~=0) % 普通線路: K=0。 Y(I,I)=Y(I,I)+Yt+Ym。 Y(J,J)=Y(J,J)+Yt+Ym。 Y(I,J)=Y(I,J)Yt。 Y(J,I)=Y(I,J)。 end if (K==0)amp。(J==0) % 對(duì)地支路: K=0,J=0,R=X=0。 Y(I,I)=Y(I,I)+Ym。 end if K0 % 變壓器線路:Zt和Ym為折算到i側(cè)的值,K在j側(cè) Y(I,I)=Y(I,I)+Yt+Ym。 Y(J,J)=Y(J,J)+Yt/K/K。 Y(I,J)=Y(I,J)Yt/K。 Y(J,I)=Y(I,J)。 end if K0 % 變壓器線路:Zt和Ym為