第二章 PSASP 潮流计算
一、实验目的
理解电力系统分析中潮流计算的相关概念，掌握用PSASP 软件对系统潮流进行计算的过程。学会在文本方式下和图形方式下的对潮流计算结果进行分析。
二、预习要求
复习《电力系统分析》中有关潮流计算的内容，了解有关潮流计算的功能，理解常用潮流计算方法，了解PQ 、PV 和V θ（平衡节点，在PSASP 中称为Slack 节点）的设置。
三、实验内容
（一） PSASP 潮流计算概述
潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。
PSASP 潮流计算的流程和结构如下图所示：
潮流计算
各种计算公共部分
图形方式
文本方式
以一个图所示9节点系统为例，计算其在常规、规划两种运行方式下的潮流。规划运行方式即在常规运行方式下，其中接于一条母线（STNB-230）处的负荷增加，对原有电网进行改造后的运行方式，具体方法为：在母线GEN3-230和STNB-230之间增加一回输电线，增加发电机3的出力及其出口变压器的容量，新增或改造的元件如下图虚线所示。
（二）数据准备
1. 指定数据目录及基准容量
双击PSASP图标，弹出PSASP封面后，按任意键，即进入PSASP主画面：
在该画面中，要完成的工作如下：
(1)指定数据目录
第一次可通过“创建数据目录”按钮，建立新目录，如：C:\CLJS。以后可通过“选择数据目录”按钮，选择该目录。
(2)给定系统基准容量
系统基准容量项中，键入该系统基准容量，如100MVA。建立了数据之后，该数不要轻易改动。
2. 录入系统潮流计算数据
基础数据(系统参数)如下：
(6)区域数据
在PSASP主画面中，点击“图形支持环境”按钮，进入图形支持环境，再点击“编辑模式”按钮，即进入系统单线图编辑窗口，分别录入系统母线、交流线、变压器、发电机和负荷的数据，以下以变压器为例：
注意：变压器是发I侧为标准侧， I、J侧互换后，变压器的等效π型等值电路不同，故其I、J侧不能互换。
（三）潮流计算基础方案的定义
在做电力系统分析计算之前，PSASP要求首先确定潮流计算的基础方案，即定义待计算电网的规模、结构和运行方式，以便从已建立的电网基础数据库中抽
方案名数据组方案解释
常规方式常规常规运行方式
规划方式常规＋新建规划运行方式在文本环境窗口中点击“计算”下拉菜单中的“方案定义”或在图形支持的运行模式窗口中点击“作业”下拉菜单中的“方案定义”，弹出方案定义窗口。
当键入新的方案名时，“已选项”框中没有内容，可从窗口下部左侧候选的数据组中加以选择。在“已选项”框中可以删除已选的数据组和改变其次序。“已选项”框中的最终内容和次序，即是新定义的方案内容。
（四）潮流计算作业的定义和执行
潮流计算的基础方案给出了待计算电网的网络结构、参数和各节点发电、负荷等基本数据，再配以不同的计算控制信息(包括发电、负荷的按比例修改等)，
作业
号
方案名计算方法
允许
误差
迭代次数
上限
电压
上限
电压
下限
1 常规方式牛拉法1e-550 1.10 0.90
在文本和图形环境下均可以定义不同的潮流计算作业，每一作业需给定一个作业号。
潮流作业的定义时要给定潮流作业号，选择该作业所对应的方案、给定电压的上下限、选择潮流计算的方法等。
点击作业页上的“刷新”按钮，按作业定义从基础数据库中重新抽取数据，以刷新该潮流作业计算数据；点击“删除”按钮，可删除该潮流作业；点击“编辑...”按钮，激活以下的数据栏，可进一步给出该潮流作业所包含的其它内容。
计算方法数据栏要选择计算方法及给定迭代时允许误差及迭代次数上限，在PSASP时，提供了5种计算方法：① PQ分解法；②牛顿法(功率式)；③最佳乘子法；④牛顿法(电流式)；⑤ PQ-牛顿法(功率式)，在计算时可根据实际系统特性选择合适的算法。
数据修改：当选择或定义了一个作业之后，便生成了该作业的计算数据。可以仅对计算数据进行调整和修改，而不影响基础数据。潮流计算的作业计算数据包括交流线、变压器、发电机、负荷和直流线的数据(参见作业数据的调整和修改)。数据回存：使用“数据回存”功能，可将计算数据所做的修改存回到基础
数据中。
定义好潮流作业后，点击“计算”按钮，执行定义好的作业的潮流计算，在计算过程窗口中显示迭代过程。关闭该窗口，返回潮流计算信息窗口。
（五）结果输出
一个潮流作业成功执行(指收敛)后，结果保存在潮流计算结果库中，查看时，可根据需要做进一步的编辑，在输出形式上有报表、图示和曲线等。
1．报表输出
在文本支持窗口中，点击“结果｜潮流”，弹出：
再点击“报表输出”按钮，进入：
设置输出基本信息，选择输出范围，输出方式，选中对象，选择变量，点击“输出”按钮，所得报表如下：
在文本支持窗口中，点击“结果｜潮流｜图示化输出”，弹出：
(1) 母线支路图示化
在图示化选择窗口中，选中“母线支路图示化”，点击“确定”按钮，弹出
母线支路选择窗口，在这窗口中设置基本信息，在对象中选择“母线”，再在母
线信息栏中选定一个母线：1，点击“确定”按钮退出后，即显示该母线结果图
如下：
(2) 区域图示化
在图示化选择窗口中，选中“区域图示化”，点击“确定”按钮，弹出区域选择窗口，设置基本信息，选择区域１，点击“确定”按钮退出后，即显示该区域结果图如下：
3．图形支持环境结果显示
在图形支持环境的运行模式下，选取菜单“视图｜潮流结果”或都点击工具条“浏览结果”按钮，弹出作业选择窗口，选择某一曾经计算过的作业号后，点
击“确定”按钮，即显示该作业的原计算结果如下：
四、实验报告
实验完成后要求完成有关潮流计算的实验报告，内容如下：
1、在文本方式下的潮流计算结果的报表（Excel格式并将其内容粘贴到Word 文档中）。
2、在图形方式下潮流计算结果的图形显示（用屏幕粘贴的方式粘贴到Word 文档中）。
3、比较作业1、2潮流计算结果，文字说明潮流结果的变化和原因。
注：实验报告交Word电子文档，文档命名为“学号＋姓名”。
简介几种潮流计算 电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算，下面简单介绍三种潮流计算方法。 一、基于多口逆向矩阵的并行潮流计算方法 多口逆向矩阵方法是求解线性方程组的普通并行方法，它只是修改了串行方法的几个部分，并且非常适用于从串行到并行的编程。该方法已用于一些电力系统并行分析方法，比如说机电暂态稳定分析和小信号稳定性，并且并行效率高。基于多口逆向矩阵方法，本文提出了一种并行牛顿潮流算法。对一个划分几个网络的大型互联系统模型的仿真结果表明这种并行算法是正确的并且效率很高。 关键词:并行潮流计算，串行潮流计算，多口逆向矩阵方法，线性方程组，电力系统分析 随着电力系统规模的扩大，尤其是区域互联网络，人们要求速度更快效率更高的功率计算，传统的串行计算越来越难满足要求，特别是对实时控制。作为电力系统的基本计算，它的效率的提高会使其他为基础的计算速度都得到提高。因为传统串行计算变的越来越难满足要求，并行计算成为提高潮流计算效率的需要。潮流计算的主要步骤是求解稀疏线性方程组，因此对并行方法的研究主要集中在线性方程组的并行求解。根据不同的实现方案，并行算法分为多因子方法、稀疏向量方法等等。多口逆向矩阵方法在各种问题中是一种求解线性方程组的通用方法。在这篇论文中，通过最常见的电力系统中的节点电压方程来说明这种方法。多口逆向矩阵法不需要在矩阵中集中调整边界点，我们根据子网的密度把矩阵分裂并且把边界节点集中在顶部，整个网络的节点电压方程组如下： 消去上矩阵中对应子网的部分，只保留边界部分。经过网络分割，边界矩阵TT Y 注入电流向量T I 被分为主控制网和各个子网。设定主控制网矩阵为 TT Y ,子网i 的为 TTi Y 。注入电流矩阵分割为子网i 为 Ti I ，即 () 31 0∑=+=k i TTi TT TT Y Y Y () 41 ∑== k i Ti T I I
潮流计算系统是根据电网的结构和参数以及系统的运行状态，运用负荷点吸收和电源点发出的有功和无功功率(PQ节点)、电压控制点的电压幅值和有功功率(PV节点)、平衡点的电压幅值和相角(Vq节点)值，来计算电网中的功率、电压、电流的分布以及各元件的功率损耗、电压损耗的过程。 潮流计算主要用于研究系统运行中：负荷变化和网络结构改变会不会危及系统安全，系统各母线电压是否在允许范围之内，系统各元件(线路、变压器等)是否过负荷，对可能出现的过负荷，事先应采取哪些预防措施，检验所提的规划方案是否满足各种运行方式的要求。另外潮流计算是系统分析计算的基础，为稳定计算提供初始运行方式。短路计算也需要以不同的潮流方式做支持。 本系统可以很直观的将电网中的潮流流向，各元件的损耗在图形中显示出来，还可以输出全网损耗表、潮流计算结果表、各元件潮流损失明细表以及各种分析报表，为我们的日常工作提供可靠的理论依据。
软件功能及特点 1. 软件功能实现智能化。 图形分层显示。只需输入变电站名称，变压器及母线名称自动生成，如有多台变压器，系统将自动给出其编号。 导线名称根据变电站名称自动生成；变电站内最高电压等级母线的颜色自动赋予该变电站；母线颜色自动赋予与其相连的导线。 图形中各设备的参数自动和图形保存在一起，方便在多台计算机上操作。图形从一台计算机复制到另外一台计算机时，如果本系统数据库中没有某型号设备的名牌参数，系统自动从图形中提取该型号设备的名牌参数，并保存到本系统数据库中，保证计算能顺利进行。 变压器能够自动分辨与其三侧相连母线的电压等级。 根据导线和变压器开关、刀闸的开、闭状态，系统自动分析此设备是否带电并给出相应的状态，并用相应的颜色标识。 1.智能动态显示有功和无功的真实潮流流向。 2. 可实现大图的分页打印，可以把全网图、变电站结构图打印成册，打印效果精确细腻。 3. 元件库和程序分离，操作者可以根据习惯自行定义。 4. 增加单独绘制图形模板功能，可绘制出变电站、母线、间隔等结构图形模板，再次增加此类图形时，可直接从模板中提取，大大提高工作效率。 5. 采用热点吸附技术。绘图时能够自动寻找热点并与其连接。 6．方便高效的图形输入平台，自动拓扑分析，图模自动转换。 7．可以输入名牌值、有名值或标幺值，也可以混合输入。 8．可以分片、分组进行计算，分片、分组显示计算结果潮流图，可通过打印机、绘图仪输出。 9.数据可以通过图形直接录入。也可以通过电子表格到入数据库。 10. 计算结果可以进行分压、分片、分组统计，可以得到变压器、导线、等元器件的详细损失情况。 11.各种输出报表均可导出Excel格式文件，方便进行数据的二次处理。 12．图形导航功能，便于整图的查找和定位。
潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新
的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图
电力系统的潮流计算 西安交通大学自动化学院 2012.10 3．1 电网结构如图3—11所示，其额定电压为10KV 。已知各节点的负荷功率及参数： MVA j S )2.03.0(2 +=， MVA j S )3.05.0(3+=， MVA j S )15.02.0(4+= Ω+=)4.22.1(12j Z ，Ω+=)0.20.1(23j Z ，Ω+=)0.35.1(24j Z 试求电压和功率分布。 解：（1）先假设各节点电压均为额定电压，求线路始端功率。 0068.00034.0)21(103.05.0)(2 2223232232323j j jX R V Q P S N +=++=++=?0019.00009.0)35.1(10 15.02.0)(2 2 224242242424j j jX R V Q P S N +=++=++=?
则： 3068.05034.023323j S S S +=?+= 1519.02009.024424j S S S +=?+= 6587.00043.122423' 12 j S S S S +=++= 又 0346 .00173.0)4.22.1(106587.00043.1)(2 2 212122'12'1212j j jX R V Q P S N +=++=++=? 故： 6933.00216.112'1212 j S S S +=?+= （2） 再用已知的线路始端电压kV V 5.101 =及上述求得的线路始端功率 12 S ，求出线 路 各 点 电 压 。
kV V X Q R P V 2752.05 .104.26933.02.10216.1)(11212121212=?+?=+=? kV V V V 2248.101212=?-≈ kV V V V kV V X Q R P V 1508.100740.0) (24242 2424242424=?-≈?=+=? kV V V V kV V X Q R P V 1156.101092.0) (23232 2323232323=?-≈?=+=? （3）根据上述求得的线路各点电压，重新计算各线路的功率损耗和线路始端功率。 0066.00033.0)21(12.103.05.02 2 223j j S +=++=? 0018.00009.0)35.1(15 .1015.02.02 2 224j j S +=++=? 故 3066.05033.023323j S S S +=?+= 1518.02009.024424j S S S +=?+= 则 6584.00042.122423' 12 j S S S S +=++= 又 0331.00166.0)4.22.1(22 .106584.00042.12 2 212j j S +=++=? 从而可得线路始端功率 6915.00208.112 j S +=
§4.5牛顿-拉夫逊法计算潮流有关问题 一、初值、收敛性和多值解 1.初值：初值选择不好，比较大，破坏了牛顿 法的基础，不收敛。选择的原则。 2.收敛性：牛顿-拉夫逊法具有平方收敛特性，高斯-塞德尔法、PQ 分解法为一阶收敛特性。 X Δ
3.多值解 对于非线性方程组，解的可能性有： ?有实际意义的解 ?有解，但在实际中无意义 （PV节点或平衡节点的无功功率超过允许值，平衡节点 的有功功率超过允许值；节点的电压过高或过低） 对策：调整运行参数，PV节点、PQ节点相互转化 ?无解，或无实数解 给定的网络结构和运行方式不合理；PV节点数目过少 对策：调整运行方式，增加PV节点 z问题很复杂，至今尚未很好解决
二、稀疏矩阵技术 1.稀疏矩阵表示法 ?节点导纳矩阵：高度稀疏的N阶复数对称方阵。因此记录矩阵的下三角。 用数组表示 数组1：记录矩阵对角元素的数值； 数组2：记录矩阵非对角元素的数值（按列存储）； 数组3：记录矩阵非对角元素的行号； 数组4：记录矩阵非对角元素的按行排的位置数；
?雅可比矩阵：高度稀疏的2N阶实数方阵，其形式对称但数值不对称。其稀疏程度与节点导纳矩阵相同，可根据节点导纳矩阵形成。
2.高斯消去法 求解牛顿-拉夫逊法潮流计算的修正方程，可以采用矩阵求逆的方法。但是由于潮流计算的雅可比矩阵通常是一个高度稀疏的矩阵，其逆阵则是一个满矩阵，因此用求逆的方法会增加额外的存储单元和计算工作量。而用高斯消去法则可以保持方程组原有的稀疏性，可以大大减少计算所需的内存和时间。
3.节点的优化编号 ?静态优化法：按静态联结支路数的多少编号。 统计好网络中各节点联结的支路数后，按联结支路数的多少，由少到多，顺序编号。 ?半动态优化法：按动态联结支路数的多少编号。 先只编一个联结支路数最小的节点号，并立即将其消去；再编消去第一个节点后联结支路数最小的节点号，再立即将其消去……依此类推。 ?动态优化法：按动态增加支路数的多少编号。 不首先进行节点编号，而是寻找消去后出现的新支路数最少的节点，并为其编号，且立即将其消去； 然后再寻找第二个消去后出现的新支路数最少的节 点并为其编号，再立即将其消去……依此类推。
电力系统分析课程设计 设计题目9节点电力网络潮流计算 指导教师 院（系、部）电气与控制工程学院 专业班级 学号 姓名 日期
电气工程系课程设计标准评分模板
目录 1 PSASP软件简介 (1) 1.1 PSASP平台的主要功能和特点 (1) 1.2 PSASP的平台组成 (2) 2 牛顿拉夫逊潮流计算简介 (3) 2.1 牛顿—拉夫逊法概要 (3) 2.2 直角坐标下的牛顿—拉夫逊潮流计算 (5) 2.3 牛顿—拉夫逊潮流计算的方法 (6) 3 九节点系统单线图及元件数据 (7) 3.1 九节点系统单线图 (7) 3.2 系统各项元件的数据 (8) 4 潮流计算的结果 (10) 4.1 潮流计算后的单线图 (10) 4.2 潮流计算结果输出表格 (10) 5 结论 (14)
电力系统分析课程设计任务书9节点系统单线图如下： 基本数据如下：
表3 两绕组变压器数据 负荷数据
1 PSASP软件简介 “电力系统分析综合程序”（Power System Analysis Software Package,PSASP）是一套历史悠久、功能强大、使用方便的电力系统分析程序，是高度集成和开发具有我国自主知识产权的大型软件包。 基于电网基础数据库、固定模型库以及用户自定义模型库的支持，PSASP可进行电力系统（输电、供电和配电系统）的各种计算分析，目前包括十多个计算机模块，PSASP的计算功能还在不断发展、完善和扩充。 为了便于用户使用以及程序功能扩充，在PSASP7.0中设计和开发了图模一体化支持平台，应用该平台可以方便地建立电网分析的各种数据，绘制所需要的各种电网图形（单线图、地理位置接线图、厂站主接线图等）；该平台服务于PSASP 的各种计算，在此之外可以进行各种分析计算，并输出各种计算结果。 1.1PSASP平台的主要功能和特点 PSASP图模一体化支持平台的主要功能和特点可概括为： 1. 图模支持平台具备MDI多文档操作界面，是一个单线图图形绘制、元件数据录入编辑、各种计算功能、结果显示、报表和曲线输出的集成环境。用户可以方便地建立电网数据、绘制电网图形、惊醒各种分析计算。人机交互界面全部汉化，界面良好，操作方便。 2. 真正的实现了图模一体化。可边绘图边建数据，也可以在数据已知的情况下进行图形自动快速绘制；图形、数据自动对应，所见即所得。 3. 应用该平台可以绘制各种电网图形，包括单线图、地理位置接线图、厂站主接线图等。 ●所有图形独立于各种分析计算，并为各计算模块所共享； ●可在图形上进行各种计算操作，并在图上显示各种计算结果； ●同一系统可对应多套单线图，多层子图嵌套； ●单线图上可细化到厂站主接线结构；
第一章 简单电力系统的分析和计算 一、 基本要求 掌握电力线路中的电压降落和功率损耗的计算、变压器中的电压降落和功率损耗的计 算；掌握辐射形网络的潮流分布计算；掌握简单环形网络的潮流分布计算；了解电力网络的简化。 二、 重点内容 1、电力线路中的电压降落和功率损耗 图3-1中，设线路末端电压为2U 、末端功率为222~jQ P S +=，则 （1）计算电力线路中的功率损耗 ① 线路末端导纳支路的功率损耗： 222 2* 222~U B j U Y S Y -=?? ? ??=? ……………（3-1） 则阻抗支路末端的功率为： 222~~~Y S S S ?+=' ② 线路阻抗支路中的功率损耗： ()jX R U Q P Z I S Z +'+'==?2 2 22222 ~ ……（3-2） 则阻抗支路始端的功率为： Z S S S ~ ~~21?+'=' ③ 线路始端导纳支路的功率损耗： 2121* 122~U B j U Y S Y -=?? ? ??=? …………（3-3） 则线路始端的功率为： 111~ ~~Y S S S ?+'= ~~~图3-3 变压器的电压和功率 ~2 ? U （2）计算电力线路中的电压降落 选取2U 为参考向量，如图3-2。线路始端电压 U j U U U δ+?+=2 1 其中 2 2 2U X Q R P U '+'= ? ； 222U R Q X P U '-'=δ ……………（3-4） 则线路始端电压的大小： ()()2 221U U U U δ+?+= ………………（3-5） 一般可采用近似计算： 2 2 2221U X Q R P U U U U '+'+ =?+≈ ………………（3-6）
家里蹲大学 电力系统稳态课程设计 题目名称：电力系统潮流计算 系别：物理与电气工程系 专业：电气工程及其自动化 学号： 姓名： qq 2316670882 指导老师：要仿真文件联系我 日期： 2014年6月3日
电力系统稳态课程设计任务书 主要内容： 一、课程设计目的 1.掌握电力系统潮流计算的基本原理； 2.掌握并能熟练运用PWS 仿真软件； 3.采用PWS 软件，做出系统接线图的潮流计算仿真结果； 二、课程设计任务 110KV 系统结线如图所示，图1中，发电厂A 装有额定功率为25+j18的发电机一台，满载运行，除供应发电机电压负荷12+j10MV A ，余下均通过两台7SF -1000/110型变压器输入系统。变压器变比为121/6.3KV 。 图1 系统结线图 变电所I 装设有两台7SF -16000/110型变压器，变比为115.5/11KV ,有如下试验数据： K P =86KW; K U %=10.5 0P =23.5KW; 0 I %=0.9 变电所II 装设有一台7SF -10000/11型变压器，变比为110/10KV ，有如下试验数据： K P =59KW; K U %=10.5 0P =16.5KW; 0I %=1.0 发电厂A 装设的两台7SF -10000/11型变压器的试验数据与变电所II 的变压器相同。 各变电所负荷、线路长度和所选导线均已示于图1。设图中与等值系统S 连接处母线电压为116KV ，试求各变电所和发电厂低压母线线电压。
基本要求： 1、按学校规定的格式编写设计论文。 2、论文主要内容有： ①课题名称。 ②设计任务和要求。 ③手算潮流和PWS的应用以及仿真结果。 ④收获体会、存在问题和进一步的改进意见等。 参考资料: [1] 何仰赞、温增银.电力系统分析[M]. 华中科技大学出版社2010.3 [2] 西安交通大学等.电力系统计算[M].北京：水利电力出版社，1993.12 [3] 陈衍.电力系统稳态分析[M].北京：水利电力出版社，2004.1 [4] 李光琦.电力系统暂态分析[M].北京：水利电力出版社，2002.5 [5] 于永源，杨绮雯. 电力系统分析（第二版）[M]. 北京：中国电力出版社，2004.3
电力系统潮流分析 —基于牛拉法和保留非线性的随机潮流 ， 姓名：*** 学号：***
1 潮流算法简介 常规潮流计算 常规的潮流计算是在确定的状态下。即：通过已知运行条件（比如节点功率或网络结构等）得到系统的运行状态（比如所有节点的电压值与相角、所有支路上的功率分布和损耗等）。 常规潮流算法中的一种普遍采用的方法是牛顿-拉夫逊法。当初始值和方程的精确解足够接近时，该方法可以在很短时间内收敛。下面简要介绍该方法。 牛顿拉夫逊方法原理 对于非线性代数方程组式（1-1），在待求量x 初次的估计值(0)x 附近，用泰勒级数（忽略二阶和以上的高阶项）表示它，可获得如式（1-2）的线性化变换后的方程组，该方程组被称为修正方程组。'()f x 是()f x 对于x 的一阶偏导数矩阵，这个矩阵便是重要的雅可比矩阵J 。 12(,,,)01,2, ,i n f x x x i n == （1-1） (0)'(0)(0)()()0f x f x x +?= （1-2） ' 由修正方程式可求出经过第一次迭代之后的修正量(0)x ?，并用修正量(0)x ?与估计值(0) x 之和，表示修正后的估计值(1)x ，表示如下（1-4）。 (0)'(0)1(0)[()]()x f x f x -?=- （1-3） (1)(0)(0)x x x =+? （1-4） 重复上述步骤。第k 次的迭代公式为： '()()()()()k k k f x x f x ?=- （1-5） (1)()()k k k x x x +=+? （1-6） 当采用直角坐标系解决潮流方程，此时待解电压和导纳如下式： i i i ij ij ij V e jf Y G jB =+=+ （1-7） 假设系统的网络中一共设有n 个节点，平衡节点的电压是已知的，平衡节点表示如下。 n n n V e jf =+ （1-8） }
电力系统潮流计算 Company number：【0089WT-8898YT-W8CCB-BUUT-202108】
电力系统 课程设计题目: 电力系统潮流计算 院系名称：电气工程学院 专业班级：电气F1206班 学生姓名： 学号： 指导教师：张孝远 1 2 节点的分类 (5) 3 计算方法简介 (6) 牛顿—拉夫逊法原理 (6) 牛顿—拉夫逊法概要 (6) 牛顿法的框图及求解过程 (8) MATLAB简介 (9) 4 潮流分布计算 (10)
系统的一次接线图 (10) 参数计算 (10) 丰大及枯大下地潮流分布情况 (14) 该地区变压器的有功潮流分布数据 (15) 重、过载负荷元件统计表 (17) 5 设计心得 (17) 参考文献 (18) 附录:程序 (19) 原始资料 一、系统接线图见附件1。 二、系统中包含发电厂、变电站、及其间的联络线路。500kV变电站以外的系统以一个等值发电机代替。各元件的参数见附件2。 设计任务 1、手动画出该系统的电气一次接线图，建立实际网络和模拟网络之间的联系。 2、根据已有资料，先手算出各元件的参数，后再用Matlab表格核算出各元件的参数。 3、潮流计算 1）对两种不同运行方式进行潮流计算，注意110kV电网开环运行。 2）注意将电压调整到合理的范围 110kV母线电压控制在106kV～117kV之间; 220kV母线电压控制在220 kV～242kV之间。 附件一：
72 水电站2 水电站1 30 3x40 C 20+8 B 2x8 A 2x31.5 D 4x7.5 水电站5 E 2x10 90+120 H 12.5+31.5 F G 1x31.5 水电站3 24 L 2x150 火电厂 1x50 M 110kV线路220kV线路课程设计地理接线示意图 110kV变电站220kV变电站牵引站火电厂水电站500kV变电站
简单环网潮流计算算例 图中所示为110kV 闭式电力网，A 为某发电厂的高压母线，其运行电压为117kV 。网络各元件的参数如下： 线路每公里的参数为 线路I 、II r 0 = 0.27 Ω， x 0 = 0.423 Ω， b 0 = 2.69×10-6S Ω III r 0 = 0.45 Ω， x 0 = 0.44 Ω， b 0 = 2.58×10-6S Ω 线路I 的长度为60km ，线路II 为50km ，线路III 为40km 。 各变电所每台变压器的额定容量、励磁功率和归算到110kV 电压级的阻抗分别为 变电所b S N =20MV A ，ΔS 0=0.05+j0.6MV A ，R T =4.84Ω，X T =63.5Ω 变电所c S N =10MV A ，ΔS 0=0.03+j0.35MV A ，R T =11.4Ω，X T =127Ω 负荷功率 S LDb =24+j18MV A ，S LDc =12+j9MV A 试求电力网的功率分布及最大电压损耗。 解： 1． 计算网络参数及制订等值电路 线路I ： var M 975.0var M 1101061.12 1Q S 1061.1S 60102.6938.25j 2.1660)423.0j 27.0(24BI 46I I -=???-=??=??=Ω +=Ω?+=---B Z 线路II ： var M 815.0var M 1101035.12 1Q S 1035.1S 50102.6915.21j 5.1350)423.0j 27.0(24BII 46II II -=???-=??=??=Ω +=Ω?+=---B Z 图（a ）
潮流计算，顾名思义是用来计算电力系统中各节点以及线路的注入功率和流动功率的。 1、对于节点处来说，要想求得节点处的注入功率，根据功率的计算公式： * ??==+I U S jQ P 必须知道各节点的电压和电流值来求得有功和无功功率；对于线路来说，要想求得流动功率，只需要用线路两端节点处的功率相加即可。 2、从上面的分析中可以看出，要想求出功率，必须知道节点电压值。这就是潮流计算的首要工作。 求节点电压的一半方法从电路中就知道可以用节点电压方程来解，应用到电力系统中，一样可以对网络列解节点电压方程。节点电压方程是通过系统节点导纳矩阵形成的电压与电流之间关系的方程。如下： B B B U Y I =其中I B 为各节点注入电流的列向量。 对于在电路中接触到的问题，应用节点电压方程求电压值是因为各节点电流量已知，而在电力系统中我们却无法知道各节点电流，所以要想利用这个方法来求电压，必须用已知量来替代电流，电力系统为我们提供了的就是各节点处的注入功率。功率与电流的关系在1中已经写明，所以以第i 个节点为例，电压方程就变成了这样一个形式： i i j n j ij i jQ P U Y U +=* =*∑ 1，电压可以表示成直角坐标形式：i i i jf e U +=，或者极坐标形式： θ∠i U 3、由于在电压方程中出现了电压相乘的情况，所以节点电压方程变成了非线性的。所以要想解出对应的电压或者功率只能采取迭代的方式。 4、电力系统中各个节点都有四个变量： 节点注入有功功率和无功功率以及节点电压的大小和相位角。 根据节点类型的不同这些变量可以分为三类：可控变量（主要指电源发出的有功、无功功率），不可控变量或者称为扰动变量（指负荷消耗的有功、无功功率），状态变量（母线或节点电压的大小和相位角，表征系统的状态）。 对于不同的节点上述四个变量的已知情况不尽相同，但是整体来说有意义的划分将节点分成了三类：已知注入功率的PQ 节点；已知注入有功和电压幅值的PV 节点；已知电压的平衡节点。 下面以直角坐标形式为例进行分析： 潮流计算的第一步要先对各节点电压设定初值，电力系统的实际情况可以允许假设所有节点电压的标幺值都为 1.0。因为各节点处的功率和电压都有一个等式，我们假设为),(;;);,(f e Q f e P i i ；所以根据假设的电压初值可以计算出节点功率值，进一步可以得出 )],()[,(][][f e Q f e P Q P Q P i i s i s i i i -=?。要想进行迭代就必须有修正方程式以获得不平衡量。 电力系统潮流计算的修正方程式的一侧是功率变化量（PQ 节点为有功和无功的变化量，PV 节点的为有功变化量和电压变化量的平方），另一侧则是电压变化量（例如直角坐标就是f e ??，），通过上段中所说方法先求得功率变化量，通过修正方程就可以得出电压变化
电力系统的潮流计算 —MATLAB的运用 20121093班：陈煜珏 指导老师：李咸善 摘要：潮流计算，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。 关键词：潮流计算、牛顿—拉夫逊算法 一、潮流计算的牛顿—拉夫逊算法的原理 牛顿—拉夫逊算法是求解非线性方程的一种有效且收敛速度快的迭代计算方法。 对于一维非线性代数方程 f(x)=0 设其准确解为x（*），而x(0)为其近似解，它与准确解之间的差为△x(0) 则x(1)= x(0)+△x(0)，此时有： f(x(1))=f(x(0)+△x(0))=0 将f(x(1))用泰勒级数展开得： f(x(1))=f(x(0)+△x(0))=f(x(0))+f'(x(0))△x(0)+[f''(x(0))(△x(0))2]/2+…
若所取的|△x(0)|足够小，则△x(0)二阶及以上阶次的各项均可略去得: f(x(0))+f'(x(0))△x(0) =0 这是对于修正量的修正方程，利用它可以解出： △x(0)=-[ f'(x(0))]-1 f(x(0)) 此时若x(0)与准确值x（*）还有差距则再次进行迭代直到满足精度要求为止。 整个迭代过程满足以下的迭代格式： △x(k)=-[ f'(x(k))]-1 f(x(k)); x(k+1)= x(k)+△x(k) k=0,1,2,3…… 牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代4~5次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值。 二、极坐标表示的牛顿—拉夫逊潮流算法 在电力系统潮流计算中，需要将全部节点分为PQ节点、PV节点和平衡节点三类。设系统有n个节点，其中有m个PQ节点，而除了
9节点电力系统潮流计算 课程设计设计题目 指导教师 院（系、部） 专业班级 学号 姓名 期9节点电力网络潮流计算 电气与控制工程学院 日 电气工程系课程设计标准评分模板 目录 1 PSASP软件简介 (1) 1.1 PSASP平台的主要功能和特点 (1) 1.2 PSASP的平台组成 (2) 2 牛顿拉夫逊潮流计算简介 (3) 2.1 牛顿—拉夫逊法概要.................................................................
3 2.2 直角坐标下的牛顿—拉夫逊潮流计算 (5) 2.3 牛顿—拉夫逊潮流计算的方法 (6) 3 九节点系统单线图及元件数据 (7) 3.1 九节点系统单线图 (7) 3.2 系统各项元件的数据 (8) 4 潮流计算的结果 (10) 4.1 潮流计算后的单线图............................................................... 10 4.2 潮流计算结果输出表格........................................................... 10 5 结论 (14) 电力系统分析课程设计任务书 9节点系统单线图如下： 基本数据如下：
表3 两绕组变压器数据 负荷数据 1 PSASP软件简介 “电力系统分析综合程序”（Power System Analysis Software Package,PSASP）是一套历史悠久、功能强大、使用方便的电力系统分析程序，是高度集成和开发具有我国自主知识产权的大型软件包。 基于电网基础数据库、固定模型库以及用户自定义模型库的支持，PSASP可进行电力系统（输电、供电和配电系统）的各种计算分析，目前包括十多个计算机模块，PSASP的计算功能还在不断发展、完善和扩充。 为了便于用户使用以及程序功能扩充，在PSASP7.0中设计和开发了图模一体化支持平台，应用该平台可以方便地建立电网分析的各种数据，绘制所需要的各种电网图形（单线图、地理位置接线图、厂站主接线图等）；该平台服务于PSASP的各种计算，在此之外可以进行各种分析计算，并输出各种计算结果。 1.1 PSASP平台的主要功能和特点 PSASP图模一体化支持平台的主要功能和特点可概括为： 1. 图模支持平台具备MDI多文档操作界面，是一个单线图图形绘制、元件数据录入编辑、各种计算功能、结果显示、报表和曲线输出的集成环境。用户可以方便地建立电网数据、绘制电网图形、惊醒各种分析计算。人机交互界面全部汉化，界面良好，操作方便。 2. 真正的实现了图模一体化。可边绘图边建数据，也可以在数
第三章 电力系统的潮流计算 3-1 电力系统潮流计算就是对给定的系统运行条件确定系统的运行状态。系 统运行条件是指发电机组发出的有功功率和无功功率（或极端电压），负荷的有 功功率和无功功率等。运行状态是指系统中所有母线（或称节点）电压的幅值和 相位，所有线路的功率分布和功率损耗等。 3-2 电压降落是指元件首末端两点电压的相量差。 电压损耗是两点间电压绝对值之差。当两点电压之间的相角差不大时， 可以近似地认为电压损耗等于电压降落的纵分量。 电压偏移是指网络中某点的实际电压同网络该处的额定电压之差。电压 偏移可以用kV 表示，也可以用额定电压的百分数表示。 电压偏移= %100?-N N V V V 功率损耗包括电流通过元件的电阻和等值电抗时产生的功率损耗和电压 施加于元件的对地等值导纳时产生的损耗。 输电效率是是线路末端输出的有功功率2P 与线路首端输入的有功功率 1P 之比。 输电效率= %1001 2 ?P P 3-3 网络元件的电压降落可以表示为 ()? ? ? ? ? +=+=-2221V V I jX R V V δ? 式中，?2V ?和? 2V δ分别称为电压降落的纵分量和横分量。 从电压降落的公式可见，不论从元件的哪一端计算，电压降落的纵、横分量计算公式的结构都是一样的，元件两端的电压幅值差主要有电压降落的纵分量决定，电压的相角差则由横分量决定。在高压输电线路中，电抗要远远大于电阻，即R X ??，作为极端的情况，令0=R ，便得 V QX V /=?，V PX V /=δ 上式说明，在纯电抗元件中，电压降落的纵分量是因传送无功功率而产生的，而电压降落的横分量则是因为传送有功功率产生的。换句话说，元件两端存在电压幅值差是传送无功功率的条件，存在电压相角差则是传送有功功率的条件。 3-4 求解已知首端电压和末端功率潮流计算问题的思路是，将该问题转化成 已知同侧电压和功率的潮流计算问题。
潮流计算（matlab）实例计算 潮流例题：根据给定的参数或工程具体要求（如图） 本设计选择Matlab 进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿- 拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到 更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f（x） = 0 的单根 附 近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不 ，收集和查阅资料；学习相关软件（软件自选：
潮流计算（matlab）实例计算 平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y 。 （2）设个节点电压的初始值U 和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也 ２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化（增加，减少或保持不变）。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“ Matrix Laboratory ”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图
电力系统 课程设计题目: 电力系统潮流计算 院系名称：电气工程学院 专业班级：电气F1206班 学生姓名： 学号： 指导教师：张孝远 1 2 2.2 节点的分类 (5) 3 计算方法简介 (6) 3.1 牛顿—拉夫逊法原理 (6) 3.1.1 牛顿—拉夫逊法概要 (6) 3.1.2 牛顿法的框图及求解过程 (8) 3.2 MATLAB简介 (9) 4 潮流分布计算 (10) 4.1 系统的一次接线图 (10) 4.2 参数计算 (10) 4.3 丰大及枯大下地潮流分布情况 (14) 4.3.1 该地区变压器的有功潮流分布数据 (15) 4.3.2 重、过载负荷元件统计表 (17)
5 设计心得 (17) 参考文献 (18) 附录:程序 (19) 原始资料 一、系统接线图见附件1。 二、系统中包含发电厂、变电站、及其间的联络线路。500kV变电站以外的系统以一个等值发电机代替。各元件的参数见附件2。 设计任务 1、手动画出该系统的电气一次接线图，建立实际网络和模拟网络之间的联系。 2、根据已有资料，先手算出各元件的参数，后再用Matlab表格核算出各元件的参数。 3、潮流计算 1）对两种不同运行方式进行潮流计算，注意110kV电网开环运行。 2）注意将电压调整到合理的范围 110kV母线电压控制在106kV～117kV之间; 220kV母线电压控制在220 kV～242kV之间。 附件一：
72 水电站2 水电站1 30 3x40 C 20+8 B 2x8 A 2x31.5 D 4x7.5 水电站5 E 2x10 90+120 H 12.5+31.5 F G 1x31.5 水电站3 24 L 2x150 火电厂 1x50 M 110kV线路220kV线路课程设计地理接线示意图 110kV变电站220kV变电站牵引站火电厂水电站500kV变电站
银川能源学院 课程设计 课程名称：电力系统分析 设计题目：简单环形网络的潮流计算 学院：电力学院 专业：电气工程及其自动化 班级：电气（本）1202班 姓名：罗通 学号：1210240073 成绩： 指导教师：李莉、张彦迪 日期：2014年12月8日—2014年12月19日
潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法，是电力系统规划和运营中不可缺少的一个重要组成部分。可以说，它是电力系统分析中最基本、最重要的计算，是系统安全、经济分析和实时控制与调度的基础。常规潮流计算的任务是根据给定的运行条件和网络结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。是电力系统研究人员长期研究的一个课题。它既是对电力系统规划设计和运行方式的合理性、可靠性及经济性进行定量分析的依据，又是电力系统静态和暂态稳定计算的基础。
前言------------------------------------------------------------------------------------------2 第一章：简单环形网络的潮流计算原理--------------------------------------4 1.1 电力线路和变压器上的功率损耗、电压降落及电能损耗--------------- 4 1.2电压降落、电压损耗、电压偏移及电压调整的概念---------------------- 5 1.3闭环网的潮流计算步骤---------------------------------------------------------- 6第二章：简单环形网络的潮流计算过程-------------------------------------- 7 2.1参数整理---------------------------------------------------------------------------- 7 2.2计算网络参数及等效电路------------------------------------------------------- 8 2.3电力系统潮流计算的运用------------------------------------------------------- 10 2.4注意事项---------------------------------------------------------------------------- 10 第三章：P-Q分解法的基本潮流算法-------------------------------------------11 3.1 P-Q分解法的原理----------------------------------------------------------------11 3.2 P-Q分解法的特点 ------------------------------------------ 13 3.3 P-Q分解法的潮流计算步骤 --------------------------------- 14 总结-------------------------------------------------------------------------------------------16谢辞-------------------------------------------------------------------------------------------17参考文献------------------------------------------------------------------------------------18