潮流上机课程设计任务书_孙英云.doc

电力系统潮流上机课程设计任 务 书一、 目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、 主要内容1 掌握计算机潮流计算的原理：a) 复习电力系统分析基础中潮流的计算机算法一章，重点掌握节点分类、潮流算法介绍b) 详细阅读牛拉法部分，掌握潮流方程（极坐标、直角坐标）的写法，掌握雅可比矩阵的公式及排列顺序和潮流方程、变量顺序的关系，掌握迭代法收敛条件及迭代法的基本原理c) 设计程序框图，划分功能模块、并对每个模块的输入输出量进行细化。2 编写计算机潮流计算程序a) 学习了解IEEE标准格式数据，学习掌握C/C+读取数据的方法b) 设计计算机数据存储母线、支路数据的结构，并将所读取的数据存放于所设计的结构当中c) 学习节点排序、节点导纳阵计算方法，编写节点导纳阵生成模块d) 编写潮流方程不平衡量计算模块e) 编写雅可比矩阵生成子模块f) 利用给定的pfMatrix类，编写修正量计算模块g) 实现潮流计算主程序，并利用IEEE标准节点数据进行校验，要求能够输出计算结果、支路潮流等必要信息3 撰写报告4 思考题a) 潮流计算的方法有哪些？各有何特点？b) 如果交给你一个任务，请你用已有的潮流计算软件计算北京城市电网的潮流，你应该做哪些工作？（收集哪些数据，如何整理，计算结果如何分析）c) 设计中遇到的问题和解决的办法。三、 进度计划序号设计内容完成时间备注1用VC6建立工程，完成IEEE标准潮流数据的读入和输出模块一天2编写及调试形成导纳矩阵模块一天3计算不平衡功率一天4形成系数矩阵（牛顿-拉夫逊法形成雅可比矩阵）一天5牛顿法潮流计算一天6牛顿法潮流计算一天7计算线路潮流，网损，PV节点无功功率和平衡节点功率一天8写设计报告一天9验收一天四、 设计报告要求1 计算机计算潮流程序流程图、功能模块、并对每个模块的输入输出量；2 完整的潮流程序，并有程序注释。包括变量的定义，子程序的调用等，计算结果；3 回答思考题五、 考核方式总成绩=考试成绩+报告成绩 学生姓名：指导教师：孙英云 2016年6月27日参考资料1：IEEE 标准数据说明Partial Description of the IEEE Common Data Format for the Exchange of Solved Load Flow DataThe complete description can be found in the paper Common DataFormat for the Exchange of Solved Load Flow Data, Working Group on aCommon Format for the Exchange of Solved Load Flow Data, _IEEETransactions on Power Apparatus and Systems_, Vol. PAS-92, No. 6,November/December 1973, pp. 1916-1925.The data file has lines of up to 128 characters. The lines are groupedinto sections with section headers. Data items are entered in specificcolumns. No blank items are allowed, enter zeros instead. Floating pointitems should have explicit decimal point. No implicit decimal pointsare used.Data type codes: A - Alphanumeric (no special characters) I - Integer F - Floating point * - Mandatory itemTitle Data=First card in file.Columns 2- 9 Date, in format DD/MM/YY with leading zeros. If no date provided, use 0b/0b/0b where b is blank.Columns 11-30 Originators name (A)Columns 32-37 MVA Base (F*)Columns 39-42 Year (I)Column 44 Season (S - Summer, W - Winter)Column 46-73 Case identification (A)Bus Data *=Section start card *:-Columns 1-16 BUS DATA FOLLOWS (not clear that any more than BUS in 1-3 is significant) *Columns ?- ? NNNNN ITEMS (column not clear, I would not count on this)Bus data cards *:-Columns 1- 4 Bus number (I) *Columns 7-17 Name (A) (left justify) *Columns 19-20 Load flow area number (I) Dont use zero! *Columns 21-23 Loss zone number (I)Columns 25-26 Type (I) * 0 - Unregulated (load, PQ) 1 - Hold MVAR generation within voltage limits, (PQ) 2 - Hold voltage within VAR limits (gen, PV) 3 - Hold voltage and angle (swing, V-Theta) (must always have one)Columns 28-33 Final voltage, p.u. (F) *Columns 34-40 Final angle, degrees (F) *Columns 41-49 Load MW (F) *Columns 50-59 Load MVAR (F) *Columns 60-67 Generation MW (F) *Columns 68-75 Generation MVAR (F) *Columns 77-83 Base KV (F)Columns 85-90 Desired volts (pu) (F) (This is desired remote voltage if this bus is controlling another bus.Columns 91-98 Maximum MVAR or voltage limit (F)Columns 99-106 Minimum MVAR or voltage limit (F)Columns 107-114 Shunt conductance G (per unit) (F) *Columns 115-122 Shunt susceptance B (per unit) (F) *Columns 124-127 Remote controlled bus numberSection end card:-Columns 1- 4 -999Branch Data *=Section start card *:-Columns 1-16 BRANCH DATA FOLLOWS (not clear that any more than BRANCH is significant) *Columns 40?- ? NNNNN ITEMS (column not clear, I would not count on this)Branch data cards *:-Columns 1- 4 Tap bus number (I) * For transformers or phase shifters, the side of the model the non-unity tap is onColumns 6- 9 Z bus number (I) * For transformers and phase shifters, the side of the model the device impedance is on.Columns 11-12 Load flow area (I)Columns 13-14 Loss zone (I)Column 17 Circuit (I) * (Use 1 for single lines)Column 19 Type (I) * 0 - Transmission line 1 - Fixed tap 2 - Variable tap for voltage control (TCUL, LTC) 3 - Variable tap (turns ratio) for MVAR control 4 - Variable phase angle for MW control (phase shifter)Columns 20-29 Branch resistance R, per unit (F) *Columns 30-40 Branch reactance X, per unit (F) * No zero impedance linesColumns 41-50 Line charging B, per unit (F) * (total line charging, +B)Columns 51-55 Line MVA rating No 1 (I) Left justify!Columns 57-61 Line MVA rating No 2 (I) Left justify!Columns 63-67 Line MVA rating No 3 (I) Left justify!Columns 69-72 Control bus numberColumn 74 Side (I) 0 - Controlled bus is one of the terminals 1 - Controlled bus is near the tap side 2 - Controlled bus is near the impedance side (Z bus)Columns 77-82 Transformer final turns ratio (F)Columns 84-90 Transformer (phase shifter) final angle (F)Columns 91-97 Minimum tap or phase shift (F)Columns 98-104 Maximum tap or phase shift (F)Columns 106-111 Step size (F)Columns 113-119 Minimum voltage, MVAR or MW limit (F)Columns 120-126 Maximum voltage, MVAR or MW limit (F)Section end card:-Columns 1- 4 -999Loss Zone Data=Section start card-Columns 1-16 LOSS ZONES FOLLOWS (not clear that any more than LOSS is significant)Columns 40?- ? NNNNN ITEMS (column not clear, I would not count on this)Loss Zone Cards:-Columns 1- 3 Loss zone number (I)Columns 5-16 Loss zone name (A)Section end card:-Columns 1- 3 -99Interchange Data *=Section start card-Columns 1-16 INTERCHANGE DATA FOLLOWS (not clear that any more than first word is significant).Columns 40?- ? NNNNN ITEMS (column not clear, I would not count on this)Interchange Data Cards *:-Columns 1- 2 Area number (I) no zeros! *Columns 4- 7 Interchange slack bus number (I) *Columns 9-20 Alternate swing bus name (A)Columns 21-28 Area interch