1 引言

暂态稳定性分析计算是电力系统三大经典计算质疑，它是电力系统规划、设计、调度运行与控制必不可少的一项重要的计算任务。在正常的稳态运行情况下，各发电机输出的电磁转矩与原动机输入的机械转矩之间相互平衡，因而所有发电机的转子速度能一直维持恒定。但是，电力系统中不可避免地受到一些大扰动的冲击，例如，各类短路故障、大容量机组或大负荷的投入和切除等。电力系统的暂态稳定性分析，其主要目的是在大扰动的冲击时，各发电机维持同步运行的能力。

电力系统暂态稳定分析方法可以分为：时域仿真法、直接法、人工智能法、动态安全域方法。把以上方法当中的一种或者几种方法相互结合形成的混合法等。近年来，有关时域仿真算法许多的改进算法，主要目的是加快仿真速度或者提高计算精度。在数学上，暂态稳定时域仿真[1-2]可以表述为微分代数方程的初值问题。有联立求解和交替求解这两类方法。交替求解法是指微分方程组和代数方程组两者的求解彼此是独立但需要交替进行的。但会造成“交接误差”。由于显式数值积分方法的数值稳定性不太好，所以显式分离求解法的整体数值稳定性也不是很好。因此显式分离求解法在商业化的电力系统暂态稳定性计算程序中应用很少。联立求解法的基本过程为，先用隐式梯形积分公式把微分方程组进行差分化，它和代数方程组一起形成联立非线性方程组，然后求解此非线性方程组。所以，仿真算法的研究可以分为两种，一种是常微分数值积分方法，另一种是对非线性方程组的求解方法。常微分数值积分方法有改进欧拉法、隐式积分法、龙格-库塔法等，在此基础上也有改进方法。文献[3-5]研究了在电力系统中，中长期的暂态稳定性分析以及变步长变阶方法。文献[6]研究了不诚实牛顿法和变步长技术相结合的方法，提高了系统的仿真速度，验证算法的正确性和实用性。非线性方程组的求解方法一般用牛顿法迭代求解。文献[7]基于奇异摄动法，提出求解发电机状态变量和网络的节点电压的方法，计算量比牛顿迭代法要小，从而提高仿真速度。文献[8]通过引入新的状态变量，网络代数方程的形式发生了改变，计算过程中没有迭代和三角分解，大幅度的提高了计算速度。

本文应用摄动法进行电力系统暂态计算仿真研究。结合IEEE9及节点系统，从发电机的转速、机端电压、和仿真时间三个方面。验证该方法的有效性和实用性。

2 暂态稳定仿真的奇异摄动法

进行电力系统动态研究时，发电机和电力系统输电网络通常可以表达为：

(1)

0=g(x,U)

西方人饮食崇尚营养和天然，这与中国饮食讲究色香味俱全的文化有着天壤之别。区别的背后，是中国饮食为了追求色香味几近极致的“创举”，比如派生出擅长雕龙刻凤的厨师，对胡椒、生姜、八角、桂皮、大蒜等佐料推崇备至，对味精、酱油等调味品的使用习以为常……很难说谁对谁错，但可以肯定的是，饮食演变绝非一朝一夕，肯定锲入了太多的历史基因，比如餐桌上那仅供观瞻的雕刻，很大程度上来自于宫廷的幻象。

(2)

(3)忽略D绕组、Q绕组，通过补入阻尼项到转子运动方程中来近似处理。

综上可以得出：

(3)

x=xk+2εx1,U=Uk+2εU1

xk+1=xk+εf(xk,Uk)+εf(xk+1,Uk+1)

(4)

0=g(x,U)

矿井采用立斜井多水平开拓，生产集中在二、三水平，二水平标高-240m，三水平标高-517m。现主要生产采区为二水平戊三采区、三水平戊一上山采区、戊二下山采区、戊一下山采区、丁二采区。采掘工作面数目满足《煤矿安全规程》规定。存在采区接替紧张局面，二水平戊三采区、三水平戊一上山采区即将于近3年内结束，正在开发的有2个采区，即三水平戊一上山采区、三水平下延戊二上山采区。

(5)

式中：ε=h/2，h为步长。

为简单起见，式(3)、(4)可写为：

ε1:x1=f(x0,U0)

(6)

0=g(x,U)

(7)

其中，设x=xk+1,U=Uk+1，B=xk+εf(xk,Uk)很显然在式(4)中ε是一个很小的常数，可以应用奇异摄动原理进行电力系统仿真计算，把x和U分别展开为ε的幂级数，令

(8)

(9)

把式(8)、(9)代入(6)、(7)，展开和为泰勒级数，然后再对比两边ε的等次幂的系数，可得到相应的一系列方程，如果保留前面三项，则有：

ε0:x0=B

(10)

0=g(x0,U0)

(11)

x=B+εf(x,U)

关于这四首钢琴小曲的音乐分析，我们曾在《中国钢琴音乐的拓荒者—新发现的早期钢琴家李树化的照片史料》一文③中有简单的艺术分析，可供读者参考。

(12)

(13)

(14)

(15)

式中：xx1，xU1和UU1都是向量，每一个元素均可由向量x1,U1解出。

暗信号非均匀性的变化规律如图4所示，其随电子注量的变化趋势与暗电流的退化趋势类似。为了更加明显地给出暗信号非均匀性的退化趋势，图5绘出了辐照前后暗电流谱呈高斯曲线分布的情况。

因此，x0,U0可由(10)、(11)中求出，x1,U1可由(12)、(13)求出，x2,U2可由(14)、(15)求出。最后x和U可以表达为：

x≈x0+εx1+ε2x2

本文中我们采用哈希算法，主要的过程为：发送方采用单向哈希函数对消息进行计算，得到摘要并发送消息和摘要。接收方将接收到的消息，按同样方式进行哈希函数计算，并将新得出的结果与发送方的原摘要结果进行比对。如结果一致，说明消息完整。在本系统中，摘要信息的不可变，保证了需要存证信息的完整性和真实性。将需要存证的电子数据放在区块链中，避免数据被恶意篡改。

(16)

U≈U0+εU1+ε2U2

(17)

在计算中如果保留的次幂越高,精度也就会越高，但计算就会随之增大。实际计算表明，在一般的电力系统仿真计算中，保留前三项就能够满足精度要求。如果是在精度要求高时，则需要以保留更多的项来满足精度要求。

针对x0=xk+εf(xk,Uk)，x1=f(xk,Uk)，则可以得出：x0=xk+εx1

x=x0+εx1+ε2x2+o(ε3)=xk+2εx1+o(ε3)

(18)

U=U0+εU1+ε2U2+o(ε3)=Uk+2εU1+o(ε3)

(19)

采用隐式积分的方法求解微分方程时，一般先把微分方程换成一系列代数方程，再用求解代数方程的办法求出时段终值。隐式梯形积分法的差分方程可以表达为：

ε0:x0=B,0=g(x0,U0)

(20)

(21)

从而可以得出：

The comparative evaluation of the outcome after infrainguinal arterial disease (IAD) revascularization demonstrates controversial results[1].

因此，描述电力系统的动态方程(1)、(2)可以改写为：

(22)

因此，若想求解x、U，只需求解出x0、U0和x1、U1，进一步减少了计算量。

3 电力系统暂态计算的数学模型

在电力系统暂态稳定计算中，当选用的元件的数学模型不同，所得出的计算结果会受直接影响，稳定性分析计算的复杂程度也会有很大的影响。因此，选择合适的数学模型来表示元件本身，不仅计算过程简单还可以得出满足精度要求结果。数学模型是电力系统暂态计算中的一个重要的问题。在工程实际问题中，根据使用场合的不同对发电机的模型做不同程度的简化处理。本文选用发电机三阶模型，负荷为恒定阻抗模型。

科学技术的飞速前进，在弹性波理论与电磁波理论的基础上，涌现了大量的新型工程勘探技术，新型的与物理原理相匹配的数据采集设备不断引入。成本小而且勘探质量高、速度快。工程物探技术具有快速、信息量大等优点，与常规钻井方法相比，具有节省时间、节约成本、地形要求低、勘探精度高等优点。

3.1 发电机三阶实用模型

在实际电力系统，暂态计算仿真分析中，当要计及励磁系统的作用时，能选择的最简单的模型就是三阶模型。它较适用于凸极机。这种实用模型基于如下假定：

2)采用蜂窝夹层板理论计算得到四边简支类方形蜂窝夹层结构的固有频率与有限元模拟结果的一致性较好，进一步证明了采用改进Gibson公式得到的类方形蜂窝夹芯等效弹性参数的正确性，同时证明了将该振动理论运用到一般蜂窝夹层结构研究的可行性，为扩展研究其他类型蜂窝夹层结构振动特性奠定了基础。

(1)不考虑定子d绕组、q绕组的暂态；

(2)在定子电压方程中，设在速度变化不大的过渡过程中，其引起的误差很小；

● WISE-PaaS/SaaS Composer：流程可视化云端组态工具；支援客制化绘图元件，可将应用场景导入3D建模绘制与互动，并以毫秒等级的画面刷新速度，搭配WISE-PaaS/Dashboard将关键管理数据以视觉直观呈现，协助萃取数据价值与提升运营效率。

式中：x是同步发电机的状态变量;U是网络的电压向量。

(23)

在公式(23)中：为状态量，分别表示发电机功角、发电机转子角速度、发电机交轴暂态电势；ω0为发电机的同步转速；Pm为发电机机械输入功率；Pe为发电机输出有功功率；D为发电机阻尼系数；TJ为发电机转子惯性时间常数；为发电机定子开路时的励磁绕组时间常数；Eq为发电机的空载电势。图1为第i台发电机相关物理量向量图。

住院药房调剂工作看似简单，但其调剂效率与质量直接关系着临床用药的及时性与安全性。为此，本项目组应用PDCA循环法将药品调剂流程进行了优化改进，最终显著地减少了住院药房的领药次数和领药时间。经实践证明，优化后的药品调剂流程在新增设科室也可得到顺利实施。

3.2 电力网络数学模型

电力系统暂态计算和潮流、短路计算中一样，用相量把电力网络的节点电压方程可表示成：

YV=I

(24)

式中：分别为电力网络的节点电流和节点电压组成的列向量；Y表示节点导纳矩阵。上式所描述的网络方程在形式上为线性方程组，其中的导纳矩阵Y仅由电力网络的结构和参数所决定。

图1 第i台发电机的向量图

在电力系统中，为了求解电力网络方程，需要首先列写出每一个动态元件的本身的代数方程，并对其进行处理，从而才能和网络方程联立求解。

3.3 发电机与网络的连接

对于三阶同步发电机模型，在d-q坐标系下，定子电压方程都可统一表示为：

(25)

式中：分别表示同步电机q轴暂态电势和电抗以及d轴的暂态电抗。

我首先在一张白纸上画了一个锐角，拿着放大镜对着所画的角左照照，右照照，发现角的两条边变粗、变长了些，但角的度数看上去没什么变化。为了更直观些，我又画了一个直角及一个大的加号，透过放大镜看了看，发现直角都没有变化。

4 算例分析

本文的采用电科院3机9节点测试系统。仿真时间为5 s，设定故障为7号节点，0.1 s时发生三相短路故障，步长为0.001 s时，0.15 s时切出故障。

基于对浅阅读的偏见，许多人发出浅阅读是一种有害阅读的言论，似乎浅阅读就是一种很低下，有害无益，不值得提倡的阅读方式。事实是否如此呢？

4.1 仿真精度对比分析

如图2所示，表示利用牛顿法求解时，1、2号发电机电压仿真曲线。图3表示利用摄动法求解时，1、2号发电机电压仿真曲线。

由图2、3可以看出，用两种方法进行仿真时，摄动法和牛顿法的计算结果曲线基本相同。所以，从这个意义上来说，验证了它们的精度是相同的。

播种前应做好土壤消毒和种子处理这项基础工作，能够将棚内土壤病菌和种子带来的病菌有效消除，最大程度上避免病害的发生，从而减少用药和药害，是生产无公害蔬菜的基本措施之一。

图2 牛顿法求解时，1、2号发电机电压仿真曲线

图3 摄动法求解时，1、2号发电机电压仿真曲线

表1 不同算法的仿真时间对比

方法仿真时间/s牛顿法0.462摄动法0.075

4.2 仿真时间对比分析

由表1可以看出，当系统总仿真时间都为5s，牛顿法计算时间是0.462s，摄动法为0.075s，计算速度提高6倍。证明了摄动法的优越性。

5 结论

本文利用摄动法进行电力系统暂态计算，其计算速度比牛顿法要快，通过仿真，利用算例验证该方法在保持计算精度的情况下，而提高计算速度，取得较好的计算效果。

参考文献

[1] 倪以信,陈寿孙,张宝霖.动态电力系统的理论分析[M].北京：清华大学出版社,2002：147-156.

[2] 王锡凡,方万良,杜正春.现代电力系统分析[M].北京：科学出版社,2003:209-309.

[3] 汤涌,宋新立,刘文焯,等.电力系统全过程动态仿真的数值方法[J].电网技术,2002,26(9):7-12.

[4] 宋新立,汤涌,刘文焯,等.电力系统全过程动态仿真的组合数值积分算法研究[J].中国电机工程学报,2009,29(28):23-29.

[5] Astic J Y,Bihain A,Jerosolimski M.The mixed Adams BDF variable step size algorithm to simulate transient and long-term phenomena in power systems[J].IEEE Trans.on Power Systems,1994,9(2):929-935.

[6] 吴红斌,丁明.用于电力系统暂态稳定仿真的可变步长牛顿法[J].中国电机工程学报,2010,30(7):36-41.

[7] 毛承雄,樊俊,O.P.Malik,等.基于奇异摄动理论的电力系统动态仿真的一种新算法[J].电力系统及其自动化学报,1990,2(1):53-59.

[8] 毛承雄,路继明,樊俊.电力系统动态仿真新模型[J].电力系统及其自动化学报,1996,8(1):1-7.