电力系统潮流计算方法分析
电力系统潮流分析
电力系统潮流分析潮流分析是电力系统中一种重要的计算方法，用于分析电力系统中各节点电压、功率和电流的分布情况。
通过潮流分析可以评估电力系统的稳定性和可靠性，为电力系统的规划、运行和控制提供参考依据。
本文将介绍电力系统潮流分析的基本原理、计算方法以及应用范围。
一、潮流分析的基本原理在电力系统中，各节点以母线表示，节点之间通过线路连接。
潮流分析基于以下几个基本原理：1. 电压平衡原理：电力系统中的节点电压必须满足节点处功率平衡方程，即节点出注入电流之和为零。
2. 潮流方程：潮流方程描述了电力系统中各节点之间电压、功率和电流之间的关系。
潮流方程是通过母线注入导纳矩阵、支路导纳和节点注入功率来表达。
3. 网络拓扑：电力系统中的节点和线路之间形成了复杂的拓扑结构，潮流分析需要考虑节点之间的相互连接关系。
二、潮流分析的计算方法潮流分析通常采用迭代法来计算各节点的电压、功率和电流。
常用的迭代法包括高斯-赛德尔迭代法和牛顿-拉夫逊迭代法。
1. 高斯-赛德尔迭代法：该方法是最简单的潮流计算方法之一。
它通过假设电力系统中所有节点电压的初始值，逐步迭代更新节点电压，直到满足收敛条件为止。
2. 牛顿-拉夫逊迭代法：该方法通过建立功率不平衡方程的雅可比矩阵，采用牛顿迭代和拉夫逊补偿的方法来求解节点电压。
牛顿-拉夫逊迭代法具有更快的收敛速度和更高的计算精度。
三、潮流分析的应用范围潮流分析在电力系统中有广泛的应用，包括但不限于以下几个方面：1. 系统规划：潮流分析可以用于电力系统的规划和设计，评估系统瓶颈、优化系统结构和参数配置。
2. 运行控制：潮流分析可以用于电力系统的运行控制，评估节点电压的合理范围、分析负荷变化对系统的影响。
3. 网络优化：潮流分析可以用于电力系统的网络优化，寻找最优输电线路和改善电力系统的供电可靠性。
4. 风电并网：潮流分析可以用于风电并网系统的规划和运行，评估并网系统的可靠性和电力系统与风电场的相互影响。
电力系统中的潮流计算与优化方法
电力系统中的潮流计算与优化方法潮流计算是电力系统运行和规划中的重要环节，它用于计算电力系统中各节点的电压、相角、有功、无功功率以及线路、变压器等的潮流分布情况。
对电力系统进行潮流计算可以帮助电力系统运行人员了解系统的稳定性、可靠性以及容载能力，也可以为电力系统规划提供数据支持。
本文将介绍电力系统潮流计算的基本方法与优化技术。
一、潮流计算的基本方法1.1 普通潮流计算方法潮流计算的基本方法是牛顿-拉夫逊迭代法（Newton-Raphson Iteration Method）和高尔顿法（Gauss-Seidel Method）。
牛顿-拉夫逊迭代法主要是通过不断迭代求解雅可比矩阵的逆，直到迭代误差小于给定阀值时停止迭代；高尔顿法则是逐一更新所有节点的电压与相角，直至所有节点的迭代误差都小于给定阀值。
1.2 快速潮流计算方法在大型电力系统中，普通的潮流计算方法计算速度较慢。
因此，研究人员提出了一些针对快速潮流计算的方法，如快速牛顿-拉夫逊法（Fast Newton-Raphson Method）和DC潮流计算方法。
快速牛顿-拉夫逊法通过简化牛顿-拉夫逊法的迭代公式，减少计算量，提高计算速度；DC潮流计算方法则是将潮流计算问题转化为一个线性方程组的求解问题，进一步提升计算效率。
二、潮流计算的优化技术2.1 改进的潮流计算算法为了提高潮流计算的准确性和收敛速度，研究人员提出了一些改进的潮流计算算法。
其中，改进的牛顿-拉夫逊法（Improved Newton-Raphson Method）是一种结合牛顿-拉夫逊法和割线法的算法，通过混合使用这两种方法，实现在减小迭代误差的同时加快计算速度。
此外，基于粒子群优化算法（Particle Swarm Optimization）和遗传算法（Genetic Algorithm）的潮流计算算法也得到了广泛研究和应用。
2.2 潮流优化潮流计算不仅可以用于分析电力系统的工作状态，还可以作为优化问题的约束条件。
有关电力系统三种潮流计算方法的比较.docx
电 力 系 统 三 种 潮 流 计 算 方 法 的 比 较一、高斯 -赛德尔迭代法：以导纳矩阵为基础， 并应用高斯 -- 塞德尔迭代的算法是在电力系统中最早得到应用的潮流计算方法，目前高斯一塞德尔法已很少使用。
将所求方程 f ( x ) 0 改写为 x( x )不能直接得出方程的根，给一个猜测值x 0 得 x 1( x 0 )又可取 x1 为猜测值，进一步得：x 2 ( x 1 )反复猜测x k 1 迭代则方程的根( x k )优点：1. 原理简单，程序设计十分容易。
2. 导纳矩阵是一个对称且高度稀疏的矩阵，因此占用内存非常节省。
3. 就每次迭代所需的计算量而言，是各种潮流算法中最小的，并且和网络所包含的节点数成正比关系。
缺点：1. 收敛速度很慢。
2. 对病态条件系统，计算往往会发生收敛困难：如节点间相位角差很大的重负荷系统、包含有负电抗支路 (如某些三绕组变压器或线路串联电容等 )的系统、具有较长的辐射形线路的系统、长线路与短线路接在同一节点上，而且长短 线路的长度比值又很大的系统。
3. 平衡节点所在位置的不同选择，也会影响到收敛性能。
二、牛顿 -拉夫逊法： 求解 f ( x ) 0设 x x 0 x ，则 按牛顿二项式展开：当 △x 不大，则取线性化（仅取一次项） 则可得修正量对 得：作变量修正：x k 1xk x k ，求解修正方程 20 世纪 牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。
自从60 年代中期采用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、计算速度方面都超过了其他方法，成为直到目前仍被广泛采用的方法。
优点：1. 收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代 4—5 次便可以收敛到一个非常精确的解。
而且其迭代次数与所计算网络的规模基本无关。
2. 具有良好的收敛可靠性， 对于前面提到的对以节点导纳矩阵为基础的高斯一塞德尔法呈病态的系统，牛顿法均能可靠地收敛。
电力系统三种潮流计算方法的比较
电力系统三种潮流计算方法的比较电力系统潮流计算是电力系统分析和运行控制中最重要的问题之一、它通过计算各节点电压和各支路电流的数值来确定电力系统各个节点和支路上的电力变量。
常见的潮流计算方法有直流潮流计算方法、高斯-赛德尔迭代法和牛顿-拉夫逊迭代法。
以下将对这三种方法进行比较。
首先，直流潮流计算方法是最简单和最快速的计算方法之一、它假设整个系统中的负载功率都是直流的，忽略了交流电力系统中的复杂性。
直流潮流计算方法非常适用于传输和配电系统，尤其是对于稳定的系统，其结果比较准确。
然而，该方法忽略了交流电力系统中的变压器的磁耦合和饱和效应，可能会导致对系统状态误判。
因此，直流潮流计算方法的适用范围有限。
其次，高斯-赛德尔迭代法是一种迭代方法，通过反复迭代计算来逼近系统的潮流分布。
该方法首先进行高斯潮流计算，然后根据计算结果更新节点电压，并再次进行计算，直到收敛为止。
高斯-赛德尔迭代法考虑了变压器的复杂性，计算结果比直流潮流计算方法更准确。
然而，该方法可能发生收敛问题，尤其是在系统变压器的串联较多或系统中存在不良条件时。
此外，该方法的计算速度较慢，尤其是对于大型电力系统而言。
最后，牛顿-拉夫逊迭代法是一种基于牛顿法的迭代方法，用于解决非线性潮流计算问题。
该方法通过线性化系统等式并迭代求解来逼近系统的潮流分布。
与高斯-赛德尔迭代法相比，牛顿-拉夫逊迭代法收敛速度更快，所需迭代次数更少。
此外，该方法可以处理系统中的不平衡和非线性元件，计算结果更准确。
然而，牛顿-拉夫逊迭代法需要建立和解算雅可比矩阵，计算量相对较大。
综上所述，电力系统潮流计算方法根据应用需求和系统特点选择合适的方法。
直流潮流计算方法适用于稳定的系统，计算简单、快速，但适用范围有限。
高斯-赛德尔迭代法适用于一般的交流电力系统，考虑了变压器复杂性，但可能存在收敛问题和计算速度较慢的缺点。
牛顿-拉夫逊迭代法适用于复杂的非线性系统，收敛速度快且计算结果准确，但需要较大的计算量。
电力系统潮流计算
电力系统潮流计算电力系统潮流计算是电力系统运行分析中的重要环节。
它通过对电力系统中各节点的电压、相角以及功率等参数进行计算和分析，从而得出电力系统的稳态运行状态。
本文将从潮流计算的基本原理、计算方法、应用及其发展等方面进行阐述。
一、潮流计算的基本原理电力系统潮流计算的基本原理是基于潮流方程建立的。
潮流方程是一组非线性的方程，描述了电力系统中各节点的电压、相角以及功率之间的关系。
潮流计算的目的就是求解这组非线性方程，以确定电力系统的电压幅值、相角及有功、无功功率的分布情况。
二、潮流计算的基本方法潮流计算的基本方法主要有直接法、迭代法以及牛顿-拉夫逊法。
直接法是通过直接求解潮流方程得到电力系统的潮流状况，但对于大规模复杂的电力系统来说，直接法计算复杂度高。
迭代法是通过对电力系统的节点逐个进行迭代计算，直到满足预设的收敛条件。
牛顿-拉夫逊法是一种较为高效的迭代法，它通过近似潮流方程的雅可比矩阵，实现了计算的高效和稳定。
三、潮流计算的应用潮流计算在电力系统运行与规划中起着重要作用。
首先，潮流计算可以用于电力系统的稳态分析，确定电力系统在各种工况下的电压、相角等参数，以判断电力系统是否存在潮流拥挤、电压失调等问题。
其次，潮流计算还可以用于电力系统的优化调度，通过调整电力系统的发电机出力、负荷组织等参数，以改善电力系统的经济性和可靠性。
此外，潮流计算还可以用于电力系统规划，通过对电力系统进行潮流计算，可以为新建电源、输电线路以及变电站等设备的规划和选择提供科学依据。
四、潮流计算的发展随着电力系统的规模不断扩大和复杂度的提高，潮流计算技术也得到了迅速的发展。
传统的潮流计算方法在计算效率和计算精度上存在一定的局限性。
因此，近年来研究者提出了基于改进的迭代方法、高精度的求解算法以及并行计算等技术，以提高潮流计算的速度和准确性。
此外，随着可再生能源的不断融入电力系统，潮流计算还需要考虑多种能源的互联互通问题，这对潮流计算提出了新的挑战，需要进一步的研究和改进。
电力系统潮流计算方法分析
电力系统潮流计算方法分析1.黎曼法是最简单和最直接的计算方法。
该方法直接利用电力系统的基本方程式，即功率平衡方程式和节点电压方程式来计算潮流分布。
然而，黎曼法需要利用复杂的矩阵方程式来解决系统中节点电压的计算，计算量大且计算速度较慢，对大型复杂系统不适用。
2.高斯-赛德尔法是一种迭代法，将电网中的节电清设置为未知数，并采用全局迭代求解。
该方法通过迭代计算不断逼近潮流分布，直到满足系统中所有节点的电压和功率平衡方程为止。
高斯-赛德尔法具有迭代次数多、耗时较长的缺点，但计算稳定可靠，对于小型系统具有较好的适用性。
3.牛顿-拉夫逊法是一种基于牛顿迭代思想的高效潮流计算方法。
该方法通过利用电力系统中的雅可比矩阵，将潮流计算问题转化为解非线性方程组的问题。
牛顿-拉夫逊法的迭代速度和稳定性较高，适用于大型复杂系统的潮流计算。
综上所述，电力系统潮流计算方法可以选择黎曼法、高斯-赛德尔法和牛顿-拉夫逊法等不同的算法进行计算。
选择合适的计算方法应根据系统的规模、复杂度以及计算时间要求来综合考虑。
实际应用中，通常会根据具体情况采用不同的方法进行潮流计算，以获得准确和高效的结果。
同时，随着电力系统的发展和智能化技术的应用，也出现了一些基于机器学习和深度学习的潮流计算方法。
这些方法利用大数据和智能算法，通过学习和分析系统历史数据，能够更好地预测和计算系统潮流分布，提高计算效率和准确性。
这些方法在未来的电力系统潮流计算中具有潜力和广阔的应用前景。
总结起来，电力系统潮流计算是电力系统分析和规划的重要工作，不同的计算方法有不同的优劣势，合理选择计算方法对于准确评估系统稳定性和可靠性至关重要。
随着技术的进步和应用的发展，电力系统潮流计算方法也在不断演化和改进，以满足电力系统智能化和可持续发展的需求。
电力系统潮流计算方法分析
电力系统潮流计算方法分析电力系统潮流计算是电力系统运行中的基础性分析方法之一，它用于求解电力系统中各个节点的电压、相角以及线路的功率、电流等变量。
潮流计算是电力系统规划、运行和控制等方面的重要工具。
本文将对电力系统潮流计算方法进行分析。
电力系统潮流计算方法主要有两种，即直接法和迭代法。
直接法又分为解析法和数值法，迭代法包括高斯赛德尔迭代法、牛顿-拉夫逊迭代法等。
解析法是通过电力系统各个节点之间的网络拓扑关系和节点电压平衡条件的方程式，直接求解节点电压和线路功率等参数。
解析法的优点是计算速度快，但其适用范围较窄，主要适用于小型简单电力系统，对于大型复杂电力系统的潮流计算会出现计算量庞大的问题。
数值法是通过将连续变量离散化，将微分方程转化为差分方程，并利用数值解法求解离散的方程组来得到电力系统潮流计算结果。
数值法的优点是适用范围广，能够处理大型复杂电力系统的潮流计算，但其缺点是计算速度相对较慢。
在迭代法中，高斯赛德尔迭代法是一种经典的迭代法，它通过先假设节点电压的初值，然后利用节点注入功率与节点电压之间的关系不断迭代计算，最终达到收敛条件为止。
高斯赛德尔迭代法的优点是收敛速度快，计算精度高，但其缺点是收敛性有时不易保证，并且计算速度会随着系统规模的增大而变慢。
牛顿-拉夫逊迭代法是一种基于牛顿迭代法的改进方法，它引入雅可比矩阵，通过牛顿迭代法的迭代过程来求解节点电压和线路功率等参数。
牛顿-拉夫逊迭代法的优点是收敛性好，计算速度快，但其缺点是在实际应用中需要预先计算雅可比矩阵，会增加计算的复杂度。
综上所述，电力系统潮流计算方法有直接法和迭代法两种，其中直接法包括解析法和数值法，迭代法包括高斯赛德尔迭代法和牛顿-拉夫逊迭代法。
在实际应用中，根据电力系统的规模和复杂程度选择合适的方法进行潮流计算，以得到准确可靠的计算结果。
此外，随着计算机技术的不断发展，还可以利用并行计算和分布式计算等方法来提高潮流计算的效率。
电力系统中的潮流计算方法及精度评估研究
电力系统中的潮流计算方法及精度评估研究概述电力系统潮流计算是电力系统运行和规划的关键技术之一。
它用于计算电力系统中各节点的电压和功率流向，以评估系统的稳定性、安全性和经济性。
本文将介绍电力系统中常用的潮流计算方法，并探讨潮流计算结果的精度评估方法。
一、潮流计算方法1. 高斯-赛德尔迭代法高斯-赛德尔迭代法是最早应用于电力系统潮流计算的方法之一。
该方法通过迭代计算每个节点的电压值，直到满足潮流平衡方程。
然而，由于其收敛速度较慢，只适用于较小规模的电力系统。
2. 牛顿-拉夫逊迭代法牛顿-拉夫逊迭代法是目前应用较广的潮流计算方法。
该方法通过建立潮流计算的牛顿方程组，并迭代求解节点电压值。
相比高斯-赛德尔迭代法，牛顿-拉夫逊迭代法具有更快的收敛速度和更好的稳定性。
3. 直流潮流计算法直流潮流计算法是一种快速计算潮流的方法，主要用于大规模电力系统的运行和规划。
该方法基于直流潮流模型，忽略了交流系统中的谐波和动态特性，降低了计算的复杂性。
然而，由于其模型简化，直流潮流计算法在评估系统安全性和稳定性方面的准确性较低。
二、潮流计算结果的精度评估1. 误差分析法误差分析法是一种常用的潮流计算结果的精度评估方法。
它通过比较潮流计算结果与实际测量值之间的差异来评估计算结果的准确性。
误差分析法通常涉及计算误差、输入误差和观测误差等方面的考虑。
2. 灵敏度分析法灵敏度分析法是一种用于评估潮流计算结果的精度和稳定性的方法。
通过计算各个输入参数对潮流计算结果的影响程度，可以评估计算结果对输入参数变化的敏感度，并识别不确定性因素。
3. 置信区间分析法置信区间分析法是一种用于评估潮流计算结果的不确定性的方法。
它通过构建置信区间，表示潮流计算结果的可信程度。
置信区间分析法可以在统计学框架下对潮流计算结果进行准确的可信度评估。
三、研究展望1. 基于深度学习的潮流计算方法近年来，深度学习在电力系统领域取得了显著的应用成果。
基于深度学习的潮流计算方法能够利用大量的数据和高级模型进行潮流计算，提高计算效率和准确性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
电力系统潮流分析—基于牛拉法和保留非线性的随机潮流，姓名：***学号：***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 ijV e jf Y G jB =+=+ （1-7）假设系统的网络中一共设有n 个节点，平衡节点的电压是已知的，平衡节点表示如下。
n n n V e jf =+（1-8）}除了平衡节点以外的所有2(1)n -个节点是需要求解的量。
每个节点可列出两个方程式。
假定系统中前m 个节点为P-Q 节点,第1m +到1n -个节点为P-V 节点。
对于PQ 节点，i P 和i Q 的值是固定的，对于PV 节点，i P 和i V 的值是固定的。
()()01,2,,()()0i is ij ij i ij j ij j j j jj i j i ij ij ij j j ij j i is i j j j i j i i m f f f e G e G e P P B B Q Q f f f G e e G e B B ∈∈∈∈⎧∆=---+=⎪=⋅⋅⋅⎨∆=--++=⎪⎩∑∑∑∑ （1-9）2222()()01,2,,1()0i is ij ij i ij jijj ji jj i j ii is ii i m m n ff fe G e G e P P B Bf V V e ∈∈⎧∆=---+=⎪=++⋅⋅⋅-⎨⎪∆=-+=⎩∑∑（1-10）选定电压初始值，按泰勒级数展开,忽略,i i e f ∆∆二次方程及以后各项,得到修正方程如下:W J U ∆=-∆（1-11）其中：22111111T mmm m n n W P Q P Q P UP U++--⎡⎤∆=∆∆∆∆∆∆∆∆⎣⎦，[]11111Tm mm m n n U e fe f e f e f ++--∆=∆∆∆∆∆∆∆∆，11111111111111111111111111111111m m m m n n m m m m n n m m m m m m m m m m m n P P P P P P P P ef e f e f e f Q Q Q Q Q Q Q Q ef e f e f e f P P P P P P P e f e f e f e J ++--++--++-∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∂∂∂∂∂∂∂∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∂∂∂∂∂∂∂∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∂∂∂∂∂∂=1111111111111111111111222221111111m n m m m m m m m m m m m m n n m m m m m m m m mm m m n n m m m m m m m P f Q Q Q Q Q Q Q Q e f e f e f e f P P P P P P P P ef e f e f e f U U U U U ef e f e -++--++++++++++--+++++∂∆∂∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∆∂∆∂∆∂∆∂∆∂∂∂∂∂2221111111111111111111112222221111111111m m m m m n n n n n n n n n n m m m m n n n n n n n n mmm m U U U f e f P P P P P P P P e f e f e f e f U U U U U U ef e f e f +++++----------++--------++∂∆∂∆∂∆∂∂∂∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∂∂∂∂∂∂∂∂∆∂∆∂∆∂∆∂∆∂∆∂∆∂∂∂∂∂∂221111n n n n U U e f ----⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥∂∆⎢⎥⎢⎥∂∂⎣⎦雅克比矩阵J 各元素的计算公式如下：22()0ii ij i ij i jj i iij i ij i j j j j P Q G e B f e f P Q B e G f j i f e U U e f ⎧∂∆∂∆=-=-+⎪∂∆∂∆⎪⎪∂∆∂∆⎪==-≠⎨∂∆∂∆⎪⎪∂∆∂∆⎪==∂∂⎪⎩（1-12）【111122()()()()22n iij j ij j ii i ii j in iij j ij j ii i ii ii j jn iij j ij j ii i ii ij i ni ij j ij j ii i ii i j j iiji i iP G e B f G e B fe P Gf B e G f B e f Q G f B e G f B e e Q G e B f G e B f f U e e U f f ====∂∆⎧=----⎪∂⎪⎪∂∆=-+-+⎪∂⎪⎪∂∆⎪=+-+∂⎪⎨∂∆⎪=-∆-++⎪∂⎪∂∆⎪=-⎪∂⎪⎪∂∆=-∂⎩∑∑∑∑j i =⎪（1-13）一般雅克比矩阵表示为：()()()()()()()()()()ij i ij i i ij ij j ij j ii i ii i j j iij i ij i iij ij j ij j ii i ii i j j iij i ij i i ijij j ij j ii i ii i j j i G e B f j i P H G e B f G e B f j i eB e G f j i P N G f B e B e G f j i f B e G f Q M G f B e B e G f e ∈∈∈-+⎧≠∂∆⎪==⎨----=∂⎪⎩-⎧≠∂∆⎪==⎨-++-=∂⎪⎩-⎧∂∆⎪==⎨++-∂⎪⎩∑∑∑22()()()()()()0()2()0()2()ij i ij i i ijij j ij j ii i ii i j i j ii ij i j i iji j j i j i j iG e B f j i Q L G e B f G e B f j i f j i U R e j i e j i U S f j i f ∈⎧⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪≠⎪⎪==⎨⎪⎪+⎧≠⎪∂∆⎪==⎨⎪--++=∂⎪⎪⎩⎪⎪≠⎧∂∆==⎪⎨-=∂⎩⎪⎪≠⎧∂∆⎪==⎨⎪-=∂⎩⎩∑ （1-14）牛顿拉夫逊方法求解框图如下：保留非线性法求解过程/图 牛顿拉夫逊潮流计算法求解框图与牛顿法的不同之处在于，第一是假设雅克比矩阵在迭代过程中不变，即取初值θ和U形成的雅克比矩阵来迭代；第二是计算出来的修正量一直是初始值的修正量。
由于保留非线性只对直角坐标形式的公式不存在截断误差，因此为了减小计算误差，本文以直角坐标形式的牛拉法为基础编写了保留非线性潮流计算方法的程序。
迭代公式为：∆x(k+1)＝－J-1[y(x(0))－y s＋y(∆x(k))] （1-14）迭代过程和牛拉法相类似，流程图如下所示：图保留非线性法求解框图蒙特卡罗模拟法蒙特卡罗模拟原理《蒙特卡罗模拟方法的思想是，是当求解问题是一不确定事件的平均值时，我们通过构建模型并采用某特定的“实验”，就可以实验中此事件发生的频率去估算概率。
蒙特卡罗模拟步骤1）根据不同新能源的特点建立新能源输出功率的样本，规模为N ；2）将得到的N 个样本值带入对应接入新能源的各节点，得到接入光伏后的各节点的值。
3）按照所述的牛顿拉夫逊法进行确定性潮流计算，得到N 组关于节点的电压，支路功率与网损的数据等。
4）运用数学上的统计原理，可以求出输出变量的分布情况。
拉丁超立方采样法拉丁超立方采样原理拉丁超立方采样由M. 、和在1979年提出，它通过分层采样使采样点能够覆盖到整个随机变量的分布范围。
该方法分成两步:1）采样：所有的输入变量可以通过分层采样，使得样本点更加准确均匀的分布；—2）排列：改变初次采样得到的样本数据的顺序，令变量数据之间的关联程度最小，或者通过排序达到指定的相关系数。
拉丁超立方采样优点1)可以使采样得到的数据较为全面地覆盖变量所分布的范围，同时分层使得采样时不会再采到一样或相似的数据，更准确地体现变量的总体情况，同时减小了样本规模。
一些文献证明了拉丁超立方采样与简单随机采样在采样规模同是M 时，两种方法抽取到的变量假设是独立的，那么它们的联合覆盖空间百分比平均值表示如下：221100%1100%1l m M P M M P M ⎧-⎡⎤=⨯⎪⎢⎥⎪⎣⎦⎨-⎪⎡⎤=⨯⎪⎢⎥+⎣⎦⎩（1-16）可以看出,当M 大于等于2时，一式大于二式，表明拉丁超立方采样比随机采样覆盖的范围大。
比如当M=20时，按式（1-16）计算得：90.25%l P =,81.86%m P =.2)拉丁超立方采样的稳健性好。
假设一输出随机变量Y 满足下式：1ni i i Y c X ==∑（1-17）i c 是常数，Y 是输入随机变量i X 的线性函数。
在相同采样规模下，进行一定次数的蒙特卡罗模拟，每一次都能获得一个关于Y 的分布情况。
由每个Y 的分布的期望值可以得到一个新的分布。
用方差Z σ表示这个分布的离散程度。
若Z σ越大，表明不同仿真间的差异越大，算法的稳健性越不好。