吐哈油田紅連線配電網(wǎng)理論線損計(jì)算

1、.電力系統(tǒng)分析課程設(shè)計(jì) 任務(wù)書題 目吐哈油田紅連線配電網(wǎng)理論線損計(jì)算 學(xué)生姓名學(xué)號專業(yè)班級設(shè)計(jì)（論文）內(nèi)容及基本要求1）設(shè)計(jì)要求 掌握電力系統(tǒng)潮流計(jì)算的基本概念和計(jì)算方法，并用于實(shí)際電網(wǎng)的理論線損計(jì)算當(dāng)中；針對吐哈油田紅連線配電網(wǎng)，進(jìn)行線損計(jì)算程序設(shè)計(jì)。 其目的在于加深學(xué)生對電力系統(tǒng)穩(wěn)態(tài)分析課程中基本概念和計(jì)算方法的掌握和理解，培養(yǎng)學(xué)生運(yùn)用所學(xué)知識分析和解決問題的能力。2）內(nèi)容學(xué)習(xí)并熟悉VC+編程語言的環(huán)境。掌握變壓器、線路等元件的等值電路和等值參數(shù)的求解方法，掌握二叉樹理論，并根據(jù)網(wǎng)絡(luò)接線圖寫出網(wǎng)絡(luò)的數(shù)據(jù)結(jié)構(gòu)。掌握基本的、常用的理論線損的計(jì)算方法。針對吐哈油田紅連線配電網(wǎng)，進(jìn)行理論線損計(jì)算。

2、利用VC+語言編程并完成調(diào)試，并對計(jì)算結(jié)果進(jìn)行分析。整理課程設(shè)計(jì)論文。起止時(shí)間指導(dǎo)教師簽名2010年6 月 日系（教研室）主任簽名2010年6 月 日學(xué)生簽名2010年6 月 日*;目錄 一、設(shè)計(jì)目的意義1二、課程設(shè)計(jì)原理1三、課程設(shè)計(jì)內(nèi)容1四、結(jié)果分析3五、設(shè)計(jì)總結(jié)5六、同組同學(xué)名單.5七、參考文獻(xiàn)5附錄.6一、設(shè)計(jì)目的意義近年來，能源危機(jī)遍及世界各國，開發(fā)和節(jié)約能源已成為當(dāng)務(wù)之急。能源危機(jī)促使人們重視能源的開發(fā)和利用，節(jié)能本身就是一種能源，是一種不可忽視的能源，是一種必須加以管理的能源。電能是油田生產(chǎn)的重要動力，隨著油田的發(fā)展、油氣勘探開發(fā)的深入，用電量將不斷增大,同時(shí)電能在傳輸、使用中的

3、損耗隨之增大。因此，做好油田配電網(wǎng)的優(yōu)化，切實(shí)有效地抓好電力管理工作，減少電能在生產(chǎn)輸送、分配及利用中的損耗，提高電能的利用效率，提高決策的科學(xué)性，對于確保油田的正常生產(chǎn)，提高油田的經(jīng)濟(jì)效益具有十分重要的意義。長期以來我國油田的6kV、10kV配電網(wǎng)建設(shè)未得到應(yīng)有的重視，處于資金短缺、網(wǎng)絡(luò)結(jié)構(gòu)陳舊、技術(shù)裝備落后、自動化水平低、事故發(fā)生時(shí)恢復(fù)供電所需時(shí)間長的狀況，造成了我國油田配電網(wǎng)網(wǎng)損長期居高不下。油田的單位電耗逐年上升，采油成本逐年增加。油田配電網(wǎng)的優(yōu)化運(yùn)行可以達(dá)到減低網(wǎng)損，具有巨大的經(jīng)濟(jì)意義。近年來人們采用各種控制方案解決石油開采中的普通存在的電能浪費(fèi)現(xiàn)象，降低網(wǎng)損也是配電網(wǎng)絡(luò)重構(gòu)和網(wǎng)絡(luò)規(guī)

4、劃的重要目標(biāo)之一。油田配電網(wǎng)的優(yōu)化運(yùn)行和優(yōu)化規(guī)劃的基礎(chǔ)是油田配電網(wǎng)的理論線損計(jì)算。二、課程設(shè)計(jì)原理前推回代潮流計(jì)算法： 前推回代潮流計(jì)算法針對10(6)kV配電網(wǎng)的單電源樹狀結(jié)構(gòu)特點(diǎn)，在利用節(jié)點(diǎn)雙親孩子兄弟鏈表對節(jié)點(diǎn)與支路進(jìn)行順層次和逆層次遍歷的基礎(chǔ)上，很容易實(shí)現(xiàn)配電網(wǎng)的“前推回代”潮流算法。潮流計(jì)算的范圍是從饋線出口(樹根)至各配變高壓側(cè)(樹葉)。節(jié)點(diǎn)和支路順層次遍歷主要用于“前推”過程，即由根節(jié)點(diǎn)(饋線首端)的電壓和支路阻抗、支路電流計(jì)算各節(jié)點(diǎn)電壓。節(jié)點(diǎn)和支路逆層次遍歷則主要用于“回代”過程，即由各葉子節(jié)點(diǎn)(末端節(jié)點(diǎn))的功率(電流)向上層累加求得其雙親節(jié)點(diǎn)和相關(guān)支路的功率(電流)。在某一層

5、所有節(jié)點(diǎn)和支路的功率(電流)都求得后，即可按同樣方法求得其雙親層所有節(jié)點(diǎn)和支路的功率(電流)。這樣就避免了求解潮流方程過程中大量的矩陣運(yùn)算23，利用前推回代潮流計(jì)算的方法，直接求出各點(diǎn)電壓進(jìn)而求出各段的功率損耗。 三、課程設(shè)計(jì)內(nèi)容1學(xué)習(xí)并熟悉VC+編程語言的環(huán)境。2掌握變壓器、線路等元件的等值電路和等值參數(shù)的求解方法，掌握二叉樹理論，并根據(jù)網(wǎng)絡(luò)接線圖寫出網(wǎng)絡(luò)的數(shù)據(jù)結(jié)構(gòu)。3掌握基本的、常用的理論線損的計(jì)算方法。4針對吐哈油田紅連線配電網(wǎng)，進(jìn)行理論線損計(jì)算。5利用VC+語言編程并完成調(diào)試，并對計(jì)算結(jié)果進(jìn)行分析。6整理課程設(shè)計(jì)論文。吐哈油田紅連線（10個(gè)節(jié)點(diǎn) ） 表1 各節(jié)點(diǎn)采用的變壓器型號和負(fù)荷電

6、量負(fù)荷節(jié)點(diǎn)編號變壓器型號有功電量（kW.h）無功電量（kvar.h）4S7-80/624.86029.0655S7-80/68.610.057S7-80/68.610.059S7-80/68.610.05表2 各支路導(dǎo)線的型號和長度起始節(jié)點(diǎn)終止節(jié)點(diǎn)導(dǎo)線型號饋線長度（km）12LJ-1201.606723LJ-950.252134LJ-1200.074935LJ-950.161626LJ-1200.930267LJ-1200.122668LJ-1200.293589LJ-1200.1135810LJ-1200.4206 四、結(jié)果分析10,總有功損耗：2.20686,線路有功損耗：0.139919

7、,變壓器有功損耗：2.06694,網(wǎng)損率：3.50814,線損率：0.222423,系統(tǒng)總有功：62.9069,系統(tǒng)總有功負(fù)荷：60.7,總無功損耗：11.6853,線路無功損耗：0.171208,變壓器無功損耗：11.5141,無功網(wǎng)損率：14.8318,無功線損率：0.217309,系統(tǒng)總無功：78.7853,系統(tǒng)總無功負(fù)荷：67.1,系統(tǒng)pp：0.489759,子節(jié)點(diǎn)：1,10,9,7,5,4,.2,10,9,7,5,4,.3,5,4,.4,.5,.6,10,9,7,.7,.8,10,9,.9,.10,.理論分析結(jié)果仍與實(shí)際測量結(jié)果之間存在偏差，是由于互感器和測量儀器誤差、導(dǎo)線參數(shù)不精確

8、、未考慮沿線電壓降落、抽油機(jī)負(fù)荷的正弦近似以及負(fù)荷的非平穩(wěn)性等因素造成的通過程序的運(yùn)行，得出了線損的數(shù)據(jù)結(jié)果，其中系統(tǒng)總有功：62.9069KVA，系統(tǒng)總無功負(fù)荷：67.1KVA,而總有功損耗：2.20686KVA, 總無功損耗：11.6853KVA,所以理論線損比較小。五、設(shè)計(jì)總結(jié)本次課程設(shè)計(jì)主要是對吐哈油田紅連線配電網(wǎng)理論線損計(jì)算 。在這次課程設(shè)計(jì)的過程中，自己能夠查閱了相關(guān)資料，對課程設(shè)計(jì)的題目、要求和具體內(nèi)容深刻的分析，并完成了此次設(shè)計(jì)。通過本次設(shè)計(jì)，我能夠鞏固所學(xué)的基本理論、專業(yè)知識，并綜合運(yùn)用所學(xué)知識來解決實(shí)際的工程問題，學(xué)習(xí)工程設(shè)計(jì)的基本技能和基本方法?？傊?，此次課程設(shè)計(jì)，使我能

9、把在課堂上學(xué)習(xí)的理論知識應(yīng)用到實(shí)踐中，更好的發(fā)現(xiàn)了自己在學(xué)習(xí)中的不足之處。在設(shè)計(jì)中，通過查閱資料，解決了在設(shè)計(jì)中所遇到的一些問題。通過此次課程設(shè)計(jì)，我受益匪淺，學(xué)到了很多東西。六、同組同學(xué)的名單蘇錄峰杜敏娜梁鵬偉朱海峰白瑞濤彭興德七、參考文獻(xiàn)附錄程序說明由于程序中定義了多種變量，所以先對每種變量的含義作以解釋： intiNodeNo 節(jié)點(diǎn)編號； intiNodeType 節(jié)點(diǎn)類型；float NodeU 源點(diǎn)電壓； float NodeP 節(jié)點(diǎn)負(fù)荷有功；float NodeQ 節(jié)點(diǎn)負(fù)荷無功； int iFatherNodeNo 父節(jié)點(diǎn)編號；int iChildNodeNo 子節(jié)點(diǎn)編號；floa

10、t InputArcUnitR 輸出弧電阻；float InputArcUnitX 輸出弧電抗； float fInputArcLength 輸出弧長度； float NodeInputI 節(jié)點(diǎn)輸入電流；float NodeOutputI 節(jié)點(diǎn)輸出電流； CComplex NodeLoadI 負(fù)荷電流； float NodeTP0 節(jié)點(diǎn)變壓器空載損耗功率； float NodeTPk 節(jié)點(diǎn)變壓器短路損耗功率；float NodeTSn 節(jié)點(diǎn)變壓器的額定功率程序分析其中第一列為節(jié)點(diǎn)編號，第二列為節(jié)點(diǎn)類型（0表示電源節(jié)點(diǎn),2表示中間節(jié)點(diǎn)，5表示負(fù)荷節(jié)點(diǎn)），第三列為節(jié)點(diǎn)電壓，第四列和第五列分別為幾點(diǎn)

11、的有功P和無功Q，第六列為父節(jié)點(diǎn)編號，第七列和第八列為子結(jié)點(diǎn)編號，第九列至第12列為父節(jié)點(diǎn)到兩個(gè)子節(jié)點(diǎn)之間的單位長度電阻和電線桿個(gè)數(shù)，第13至17為變壓器的五個(gè)參數(shù)，依次為P0，PK，I0，UK，和SN。10, 1, 0, 6.3, 64.11, 92.79, -1, 2, -1, -1, -1, 0.27, 1.6067, 0, 0, 0, 0, 0, 0, 0, 0, 80, 2, 2, 6.3, 0, 0, 1, 3, 6, -1, -1, 0.33, 0.2521, 0.27, 0.9302, 0, 0, 0, 0, 0, 0, 80, 3, 2, 6.3, 0, 0, 2, 4, 5

12、, -1, -1, 0.27, 0.0749, 0.33, 0.1616, 0, 0, 0, 0, 0, 0, 80, 4, 5, 6.3, 24.86, 29.065, 3, -1, -1, -1, -1, 0.27, 0, 0.27, 0, 0, 0, 0, 0, 0.27, 1.65, 80, 5, 5, 6.3, 8.6, 10.05, 3, -1, -1, -1, -1, 0.27, 0, 0.27, 0, 0, 0, 0, 0, 0.27, 1.65, 80, 6, 2, 6.3, 0, 0, 2, 7, 8, -1, -1, 0.27, 0.1226, 0.27, 0.2935,

13、 0, 0, 0, 0, 0, 0, 80, 7, 5, 6.3, 8.6, 10.05, 6, -1, -1, -1, -1, 0.27, 0, 0.27, 0, 0, 0, 0, 0, 0.27, 1.65, 80, 8, 2, 6.3, 0, 0, 6, 9, 10, -1, -1, 0.27, 0.1135, 0.27, 0.4206, 0, 0, 0, 0, 0, 0, 80, 9, 5, 6.3, 8.6, 10.05, 8, -1, -1, -1, -1, 0.27, 0, 0.27, 0, 0, 0, 0, 0, 0.27, 1.65, 80,10, 5, 6.3, 0, 0,

14、 8, -1, -1, -1, -1, 0.33, 0, 0.33, 0, 0, 0, 0, 0, 0, 0, 80,1, 2, 3, 4, 5, 6, 7, 8, 9, 10,程序設(shè)計(jì)：/ (1) TideCalculateDlg.cpp : implementation file/#include stdafx.h#include TideCalculate.h#include TideCalculateDlg.h#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endif/

15、CAboutDlg dialog used for App Aboutclass CAboutDlg : public CDialogpublic:CAboutDlg();/ Dialog Data/AFX_DATA(CAboutDlg)enum IDD = IDD_ABOUTBOX ;/AFX_DATA/ ClassWizard generated virtual function overrides/AFX_VIRTUAL(CAboutDlg)protected:virtual void DoDataExchange(CDataExchange* pDX); / DDX/DDV suppo

16、rt/AFX_VIRTUAL/ Implementationprotected:/AFX_MSG(CAboutDlg)/AFX_MSGDECLARE_MESSAGE_MAP();CAboutDlg:CAboutDlg() : CDialog(CAboutDlg:IDD)/AFX_DATA_INIT(CAboutDlg)/AFX_DATA_INITvoid CAboutDlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(CAboutDlg)/AFX_DATA_MAPBEGIN_MESSAG

17、E_MAP(CAboutDlg, CDialog)/AFX_MSG_MAP(CAboutDlg)/ No message handlers/AFX_MSG_MAPEND_MESSAGE_MAP()/ CTideCalculateDlg dialogCTideCalculateDlg:CTideCalculateDlg(CWnd* pParent /*=NULL*/): CDialog(CTideCalculateDlg:IDD, pParent)/AFX_DATA_INIT(CTideCalculateDlg)/ NOTE: the ClassWizard will add member in

18、itialization here/AFX_DATA_INIT/ Note that LoadIcon does not require a subsequent DestroyIcon in Win32m_hIcon = AfxGetApp()-LoadIcon(IDR_MAINFRAME);void CTideCalculateDlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(CTideCalculateDlg)/ NOTE: the ClassWizard will add DD

19、X and DDV calls here/AFX_DATA_MAPBEGIN_MESSAGE_MAP(CTideCalculateDlg, CDialog)/AFX_MSG_MAP(CTideCalculateDlg)ON_WM_SYSCOMMAND()ON_WM_PAINT()ON_WM_QUERYDRAGICON()ON_BN_CLICKED(IDC_BUTTON, OnButton)/AFX_MSG_MAPEND_MESSAGE_MAP()/ CTideCalculateDlg message handlersBOOL CTideCalculateDlg:OnInitDialog()CD

20、ialog:OnInitDialog();/ Add About. menu item to system menu./ IDM_ABOUTBOX must be in the system command range.ASSERT(IDM_ABOUTBOX & 0xFFF0) = IDM_ABOUTBOX);ASSERT(IDM_ABOUTBOX AppendMenu(MF_SEPARATOR);pSysMenu-AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);/ Set the icon for this dialog. The fram

21、ework does this automatically/ when the applications main window is not a dialogSetIcon(m_hIcon, TRUE);/ Set big iconSetIcon(m_hIcon, FALSE);/ Set small icon/ TODO: Add extra initialization herereturn TRUE; / return TRUE unless you set the focus to a controlvoid CTideCalculateDlg:OnSysCommand(UINT n

22、ID, LPARAM lParam)if (nID & 0xFFF0) = IDM_ABOUTBOX)CAboutDlg dlgAbout;dlgAbout.DoModal();elseCDialog:OnSysCommand(nID, lParam);/ If you add a minimize button to your dialog, you will need the code below/ to draw the icon. For MFC applications using the document/view model,/ this is automatically don

23、e for you by the framework.void CTideCalculateDlg:OnPaint() if (IsIconic()CPaintDC dc(this); / device context for paintingSendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);/ Center icon in client rectangleint cxIcon = GetSystemMetrics(SM_CXICON);int cyIcon = GetSystemMetrics(SM_CYICON);CRe

24、ct rect;GetClientRect(&rect);int x = (rect.Width() - cxIcon + 1) / 2;int y = (rect.Height() - cyIcon + 1) / 2;/ Draw the icondc.DrawIcon(x, y, m_hIcon);elseCDialog:OnPaint();/ The system calls this to obtain the cursor to display while the user drags/ the minimized window.HCURSOR CTideCalculateDlg:O

25、nQueryDragIcon()return (HCURSOR) m_hIcon;void CTideCalculateDlg:OnButton() / TODO: Add your control notification handler code hereinti(0),j(0),k(0);intiNodeNumber(0);charcComma( ); CDistributionNodeData*pData= NULL; float Los=0.0,LosQ=0.0; FILE *fp;CString FilePath；FilePath=吐哈油田紅連線數(shù)據(jù).txt；fp=fopen(Fi

26、lePath,r)；char ctem；fscanf(fp,%d,&iNodeNumber)；fscanf(fp,%c,&ctem)；pData= new CDistributionNodeDataiNodeNumber；for(i=0;iiNodeNumber;i+)fscanf(fp,%d,&pDatai.iNodeNo)；CString FilePath;FilePath=吐哈油田紅連線數(shù)據(jù).txt ;fp=fopen(FilePath,r);char ctem;fscanf(fp,%d,&iNodeNumber);fscanf(fp,%c,&ctem);pData= new CDist

27、ributionNodeDataiNodeNumber;for(i=0;iiNodeNumber;i+)fscanf(fp,%d,&pDatai.iNodeNo);fscanf(fp,%c,&ctem);fscanf(fp,%d,&pDatai.iNodeType);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.NodeU);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.NodeP); pDatai.NodeP=pDatai.NodeP/30/24;fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDa

28、tai.NodeQ); pDatai.NodeQ=pDatai.NodeQ/30/24;fscanf(fp,%c,&ctem);fscanf(fp,%d,&pDatai.iFatherNodeNo);fscanf(fp,%c,&ctem);fscanf(fp,%d,&pDatai.iChildNodeNo0);fscanf(fp,%c,&ctem);fscanf(fp,%d,&pDatai.iChildNodeNo1);fscanf(fp,%c,&ctem);/*fscanf(fp,%d,&pDatai.iChildNodeNo2);fscanf(fp,%c,&ctem);fscanf(fp,

29、%d,&pDatai.iChildNodeNo3);fscanf(fp,%c,&ctem);*/fscanf(fp,%f,&pDatai.InputArcUnitR0);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.fInputArcLength0); pDatai.fInputArcLength0=pDatai.fInputArcLength0*55/1000;fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.InputArcUnitR1);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.f

30、InputArcLength1); pDatai.fInputArcLength1=pDatai.fInputArcLength1*55/1000;fscanf(fp,%c,&ctem);/*fscanf(fp,%f,&pDatai.InputArcUnitR2);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.fInputArcLength2);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.InputArcUnitR3);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.fInputArcL

31、ength3);fscanf(fp,%c,&ctem);*/fscanf(fp,%f,&pDatai.NodeTP0);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.NodeTPk);fscanf(fp,%c,&ctem); fscanf(fp,%f,&pDatai.NodeTI0);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.NodeTUk);fscanf(fp,%c,&ctem);fscanf(fp,%f,&pDatai.NodeTSn);fscanf(fp,%c,&ctem);/(2) 補(bǔ)各個(gè)支路的電抗值在第一次編

32、程時(shí)，數(shù)據(jù)格式中只列出了支路的電阻值，而計(jì)算過程中用到支路電抗值，因?yàn)?txt文件修改比較麻煩，所以在程序中補(bǔ)了各支路的電抗值。這也是本文數(shù)據(jù)寫入的一個(gè)缺點(diǎn)。 float Sort;for(i=0;iiNodeNumber;i+)Sort=pDatai.InputArcUnitR0；if(fabs(Sort- 0.27)=0.00001)pDatai.InputArcUnitX0=0.342；if(fabs(Sort-0.33=0.00001)pDatai.InputArcUnitX0=0.335；Sort=pDatai.InputArcUnitR1；if(fabs(Sort-0.27)=0.

33、00001)pDatai.InputArcUnitX1=0.342；if(fabs(Sort-0.33)=0.00001)pDatai.InputArcUnitX1=0.335；Sort=pDatai.InputArcUnitR2；if(fabs(Sort-0.27)=0.00001)pDatai.InputArcUnitX2=0.342；if(fabs(Sort-0.33)=0.00001)pDatai.InputArcUnitX2=0.335；Sort=pDatai.InputArcUnitR3；if(fabs(Sort-0.27)=0.00001)pDatai.InputArcUnitX

34、3=0.342；if(fabs(Sort-0.33)0;i-)tem=-1；if(pDatai.iNodeType=5)tem=i;while(tem=0)Code=pDatatem.iFatherNodeNo；if(Code=0)pDataCode-1.TwigNodepDataCode-1.TwigNum=I；pDataCode-1.TwigNum+;tem=Code-1；Loss=0；LossQ=0；/(4) 找本支路參數(shù)for(i=iNodeNumber-1;i0;i-)Code=pDatai.iFatherNodeNo；Rtem=0； Xtem=0； Ltem=0；if(Code=-

35、1)continue；if(Code!=-1)tem=pDataCode-1.iChildNodeNo0；if(tem=pDatai.iNodeNo)Rtem=pDataCode-1.InputArcUnitR0；Xtem=pDataCode-1.InputArcUnitX0；Ltem=pDataCode-1.fInputArcLength0；tem=pDataCode-1.iChildNodeNo1；if(tem=pDatai.iNodeNo)Rtem=pDataCode-1.InputArcUnitR1；Xtem=pDataCode-1.InputArcUnitX1；Ltem=pDataC

36、ode-1.fInputArcLength1；tem=pDataCode-1.iChildNodeNo2；if(tem=pDatai.iNodeNo)Rtem=pDataCode-1.InputArcUnitR2；Xtem=pDataCode-1.InputArcUnitX2；Ltem=pDataCode-1.fInputArcLength2；tem=pDataCode-1.iChildNodeNo3；if(tem=pDatai.iNodeNo) Rtem=pDataCode-1.InputArcUnitR3；Xtem=pDataCode-1.InputArcUnitX3；Ltem=pData

37、Code-1.fInputArcLength3；/(5) 計(jì)算本支路的有功損耗和無功損耗float Tatalp=0；float Power；int tem;float ftem4,ftemQ4,Utem,Utem0,ftem0；int Limit=50；Power=Tatalp=0；float LossP=0,LossQ=0,Loss2(0),Load(0),LossQ2(0),LoadQ(0)；float loadpp=0；for(i=iNodeNumber-1;i=0;i-)if(pDatai.iNodeType=5) LossP=(float)(pDatai.NodeTPk*(1.5*

38、pDatai.NodeP*pDatai.NodeP+pDatai.NodeQ*pDatai.NodeQ)/(pDatai.NodeTSn*pDatai.NodeTSn)+pDatai.NodeTP0； LossQ=0.04*(1.5*pDatai.NodeP*pDatai.NodeP+pDatai.NodeQ*pDatai.NodeQ)/(pDatai.NodeTSn)+0.024*pDatai.NodeTSn； pDatai.DltaNodeP=LossP； pDatai.DltaNodeQ=LossQ； Loss2+=LossP； LossQ2+=LossQ； Load+=pDatai.N

39、odeP； LoadQ+=pDatai.NodeQ； loadpp+=pDatai.NodeP*pDatai.NodeP； if( pDatai.NodeTP0=0 & pDatai.NodeTPk=0)pDatai.SampleType=1；else pDatai.NodeP=0； pDatai.NodeQ=0； pDatai.DltaNodeP=0； pDatai.DltaNodeQ=0； float Loss=0；float Ieff(0),AvS(0),DeltaPL(0),TemP,TemQ,TemPij；int Code(0),TemCode；float Rtem(0)；float

40、 Xtem(0)；float Ltem(0)；/(6) 計(jì)算干線的的有功損耗和無功損耗float temLossP(0),temLossQ(0);float TemP=TemQ=TemPij=0;float Temppij=0;int CodeEnd;for(j=0;jpDatai.TwigNum;+j)CodeEnd=pDatai.TwigNodej;/1.5*pi2if(pDataCodeEnd.SampleType)TemP+=pDataCodeEnd.NodeP*pDataCodeEnd.NodeP;elseTemP+=1.5*pDataCodeEnd.NodeP*pDataCodeE

41、nd.NodeP;Temppij+=pDataCodeEnd.NodeP*pDatai.DltaNodeP;TemQ+=pDataCodeEnd.NodeQ;/pijfor(k=j+1;kpDatai.TwigNum;+k)TemCode=pDatai.TwigNodek;TemPij+=2*pDataTemCode.NodeP*pDataCodeEnd.NodeP;if(pDatai.iNodeType=5)TemP=1.5*pDatai.NodeP*pDatai.NodeP;TemQ=pDatai.NodeQ;Temppij+=pDatai.NodeP*pDatai.DltaNodeP;pDataCode-1.DltaNodeP+=pDatai.DltaNodeP;pDataCode-1.DltaNodeQ+=pDatai.DltaNodeQ;TemQ=(TemQ+pDatai.DltaNodeQ)*(TemQ+pDatai.DltaNodeQ);temLossP=(TemP+TemPij+TemQ+2*Temppij+pDatai.DltaNodeP*pDatai.DltaN