CN110502859B - 用于电气耦合园区综合能源系统的多速率动态仿真方法 - Google Patents

Abstract

本发明涉及电气耦合技术领域，具体公开一种用于电气耦合园区综合能源系统的多速率动态仿真方法，所述方法包括：包括：建立电气耦合园区综合能源系统及其对应的电气耦合数学模型；采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解。本发明提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法，既具有较高的计算精度又有较高计算效率。

Description

用于电气耦合园区综合能源系统的多速率动态仿真方法

技术领域

本发明涉及电气耦合技术领域，尤其涉及一种用于电气耦合园区综合能源系统的多速率动态仿真方法。

背景技术

随着能源与环境问题的日益突出，能够实现多种能源形式融合并高效利用的电气耦合综合能源系统已经成为当今能源技术领域的研究热点。其中，电力能源清洁高效、易于转化利用、可控性强，是综合能源系统的核心能源形式；而天然气因其经济环保、储量丰富等特征，同样被认为是未来综合能源体系中的重要组成部分。由于电力和燃气的动态特征和行为模式差异巨大，在综合能源系统中，对这两种能源形式进行准确、可靠的协同建模与仿真分析成为成为研究重点。

目前，电气系统由配电网、微型燃气轮机机组及P2G装置(P2G,即Power to Gas,概念源于德国可再生能源导入量的扩大，利用风力发电、太阳能发电等的剩余电力电解水生成氢，然后提供给现有的燃气管道网络；或者利用电力、水及大气中的CO₂，通过甲烷化反应制造甲烷，提供燃气)组成，配电网模型多种多样，常用的一种是将配电网中各元件及负荷等效为恒定的阻抗，利用节点分析法求解整个系统如下所示：

I＝Y

式中：I表示阻抗网络中各节点注入电流矩阵，Y表示网络的节点导纳矩阵，U表示网络各点的电压矩阵。

天然气管网模型主要以下式表示：

式中：A为管道截面积，p为管道节点压力，M为管道的流量，D为管道的直径，c为管道内的声速，λ为摩擦系数。

配电网和天然气系统采用微型燃气轮机机组和P2G装置进行双向耦合，微型燃气轮机机组由原动机和同步发电机组成，具体模型如下：

同步发电机模型：

式中：δ为同步电机的功角，ω_*为电机的转速标幺值，ω_N为电机的额定转速，P_m为输入的机械功率的标幺值，P_e为电磁功率的标幺值，D为同步电机阻尼系数，H为惯性时间常数的一半。

原动机模型：

P_m＝f(p,M)

式中：P_m表示燃气轮机的输出功率，p表示电磁功率，M为管道流量信号。

将P2G装置看作配电网中的用电器，P2G模型如下所示：

式中：M表示P2G装置产生的气体流量，P为P2G装置需要的电功率，μ为转换效率，H为天然气燃烧的热值。

对电气耦合园区综合能源系统进行仿真,就是求解电气耦合园区综合能源系统各个节点处的电流或者电压等电力状态参数以及流量或者压力等气体状态参数的过程。

通过仿真，可以预知整个电气耦合园区综合能源系统各个节点处的状态参数，从而指导实际生产和运行。

现有的仿真方法大多是采用较小的步长，对天然气系统和电气系统在相同的时间尺度上依次进行求解，虽然求解的精度较高，但是计算所耗费的时间长，占用的空间和资源较多。

因此，急需针对研究出一种既能满足计算精度又有较高计算效率的仿真方法。

发明内容

本发明的一个目的在于，提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法，既具有较高的计算精度又有较高计算效率。

为达以上目的，本发明提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法，包括：

建立电气耦合园区综合能源系统及其对应的电气耦合数学模型；

采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解。

优选的，所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型包括：

确定配电网中线路的阻抗及变压器参数，然后建立电力能源系统拓扑连接模型；

确定天然气管网中各支路的长度、直径及摩擦系数，然后建立天然气能源系统拓扑连接模型；

分别确定微型燃气轮机机组和P2G装置二者的位置、容量及动态参数，通过所述微型燃气轮机机组和P2G装置耦合所述天然气能源系统拓扑连接模型和电力能源系统拓扑连接模型，得到所述电气耦合园区综合能源系统对应的电气耦合拓扑连接模型；

针对所述电气耦合拓扑连接模型建立电气耦合数学模型。

优选的，所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型还包括：

指定所述电气耦合数学模型的仿真参数；所述仿真参数包括电气耦合数学模型的仿真终止时间T、配电网模型的仿真步长Δt₁、天然气管网模型的仿真步长Δt₂、天然气管网模型的空间步长Δx，且满足Δt₂＝NΔt₁，N为正整数，并设置仿真的扰动事件信息。

优选的，电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型，所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型之后，还包括：

将所述天然气管网模型和同步发电机模型转化为差分化形式。

优选的，所述采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解包括：

设定外部循环仿真时刻初值t_n＝0以及设定电气耦合数学模型中若干状态变量的初值；

设置配电网模型的求解次数s＝1且s为正整数；

设置当前时刻t＝t_n+sΔt₁，以t_n时刻天然气管网末端出口流量M_k(t_n)作为微型燃气轮机机组燃料系统的注入流量，以仿真步长Δt₁求解配电网模型，得到微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率、功角、天然气管网与微型燃气轮机连接点的压力p(t)、电网中各点的电压及P2G流向天然气管网中的流量；使s＝s+1；

判断t时刻配电网模型是否有扰动发生，若有则修改电力能源系统拓扑连接模型的拓扑结构和动态参数，否则直接进入下一步；

判断s是否大于或等于N+1，若是，进入下一步；

令外部循环仿真时刻t_n＝t_n+Δt₂，将求得的天然气管网与微型燃气轮机连接点的压力p(t_n)作为天然气管网的末端压力，P2G注入天然气管网中的流量作为边界条件，以仿真步长Δt₂求解天然气管网模型，得到天然气管网系统各节点压力及各支路流量，并更新天然气管网流向微型燃气轮机的流量M_k(t_n)；

判断在t_n+Δt₂时刻天然气管网模型是否发生扰动，若有扰动，则根据扰动信息修改天然气能源系统拓扑连接模型的拓扑结构和动态参数，否则直接进行下一步骤；

判断仿真时间t_n是否大于仿真终止时刻T，若是则进入下一步；

输出仿真结果，仿真任务结束。

本发明的有益效果在于：提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法，可以用于电气耦合技术领域中的电气耦合数学模型求解情形，通过适当减少天然气管网模型的求解次数，可以以较小的计算精度损失换取较大的计算速度提升，能有效提高仿真效率。

附图说明

为了更清楚地说明本实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。

图1是本发明实施例提供的用于电气耦合园区综合能源系统的多速率动态仿真方法的流程图；

图2是本发明实施例提供的电气耦合园区综合能源系统的拓扑结构图；

图3～图10是采用本实施例所提供的多速率动态仿真方法和采用传统的单速率求解方法求解所得的各状态参数的对比示意图。

具体实施方式

为使得本发明的目的、特征、优点能够更加的明显和易懂，下面将结合本实施例中的附图，对本实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

本实施例提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法，适用于电气耦合领域中的应用场景，可以提高电气耦合园区综合能源系统模型的求解效率，所述用于电气耦合园区综合能源系统的多速率动态仿真方法由一种监控装置来执行，通过软件和/或硬件实现。

图1是本实施例一提供的用于电气耦合园区综合能源系统的多速率动态仿真方法的流程图。

参见图1，所述用于电气耦合园区综合能源系统的多速率动态仿真方法包括如下步骤：

S10：参见图1，建立电气耦合园区综合能源系统及其对应的电气耦合数学模型；

具体地，S10包括：

S101：确定配电网中线路的阻抗及变压器参数，然后建立电力能源系统拓扑连接模型1；

S102：确定天然气管网中各支路的长度、直径及摩擦系数，然后建立天然气能源系统拓扑连接模型2；

本实施例中，天然气管网的各个参数可参照下表1进行确定。

表1传输管道参数表

S103：分别确定微型燃气轮机机组3和P2G装置4二者的位置、容量及动态参数，通过所述微型燃气轮机机组和P2G装置耦合所述天然气能源系统拓扑连接模型2和电力能源系统拓扑连接模型1，得到所述电气耦合园区综合能源系统对应的电气耦合拓扑连接模型；

本实施例中，同步电机的各个参数可参照下表2进行确定。

表2同步电机参数表

本实施例中，微型燃气轮机的各个参数可参照下表3进行确定。

表3微型燃气轮机参数表

参数名称	参数值
		微型燃气轮机基准功率	250kW
有功功率参考值	0.85p.u.
		阻尼系数	0.03
燃料系统延迟时间常数	0.1s
		燃料系统延迟时间常数	1.0s
负荷限定时间常数	3.0s
		负荷限制	1.2
温度循环控制增益	1.0
		功率控制比例增益	0.01
功率控制积分增益	5

本实施例中，P2G装置的各个参数可参照下表4进行确定。

表4 P2G装置参数表

参数名称	参数值
		天然气热值	<![CDATA[35588kJ/m<sup>3</sup>]]>
P2G转化效率	0.62
		天然气密度	<![CDATA[0.7174kg/m<sup>3</sup>]]>

S104：针对所述电气耦合拓扑连接模型建立电气耦合数学模型；

一般地，电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型。可以理解的是，不同的电气耦合拓扑连接模型可以对应不同的电气耦合数学模型，为了方便说明，本实施例以背景技术中所提及的常用的电气耦合数学模型进行说明。即本实施例中：

配电网模型为：

I＝YU

式中：I表示阻抗网络中各节点注入电流矩阵，Y表示网络的节点导纳矩阵，U表示网络各点的电压矩阵。

天然气管网模型为：

式中：A为管道截面积，p为管道节点压力，M为管道的流量，D为管道的直径，c为管道内的声速，λ为摩擦系数。

同步发电机模型：

式中：δ为同步电机的功角，ω_*为电机的转速标幺值，ω_N为电机的额定转速，P_m为输入的机械功率的标幺值，P_e为电磁功率的标幺值，D为同步电机阻尼系数，H为惯性时间常数的一半。

原动机模型为：

P_m＝f(p,M)

式中：P_m表示燃气轮机的输出功率，p表示电磁功率，M为管道流量信号。

P2G模型为：

式中：M表示P2G装置产生的气体流量，P为P2G装置需要的电功率，μ为转换效率，H为天然气燃烧的热值。

S105：指定所述电气耦合数学模型的仿真参数；所述仿真参数包括电气耦合数学模型的仿真终止时间T、配电网模型的仿真步长Δt₁、天然气管网模型的仿真步长Δt₂、天然气管网模型的空间步长Δx，且满足Δt₂＝NΔt₁，N为正整数，并设置仿真的扰动事件信息。

S20：电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型，将所述天然气管网模型和同步发电机模型转化为差分化形式；

具体地，

1)同步电机差分化模型为：

式中，δ为同步电机的功角，ω_*为电机的转速标幺值，ω_N为电机的额定转速，P_m为输入的机械功率的标幺值，P_e为电磁功率的标幺值，D为同步电机阻尼系数，H为惯性时间常数的一半。

2)天然气管网差分化模型为：

式中，p_m表示每一微元管道的末端压力，p_k表示每一微元管道的始端压力，M_k表示每一微元管道的流入流量，M_m表示每一微元管道的流出流量，A为管道截面积，D为管道的直径，c为管道内的声速，λ为摩擦系数。

S30：采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解。

具体地，S30包括：

S301：设定外部循环仿真时刻初值t_n＝0以及设定电气耦合数学模型中若干状态变量的初值。

具体地，若干状态变量可以包括天然气管道的支路流量和节点压力以及微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率和功角等。

S302：设置配电网模型的求解次数s＝1且s为正整数；

S303：设置当前时刻t＝t_n+sΔt₁，以t_n时刻天然气管网末端出口流量M_k(t_n)作为微型燃气轮机机组燃料系统的注入流量，以仿真步长Δt₁求解配电网模型，得到微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率、功角、天然气管网与微型燃气轮机连接点的压力p(t)、电网中各点的电压及P2G流向天然气管网中的流量；使s＝s+1；

具体地，S303包括：

S3031：求解同步发电机模型，根据t时刻微型燃气轮机输出的机械功率P_m(t)和电网电磁功率P_e(t)，求解同步电机差分化模型，得到t+Δt₁时刻电机的功角δ(t+Δt₁)以及微型燃气轮机机组的转速ω_r*(t+Δt₁)；

S3032：求解注入电网侧的电流，根据同步电机功角δ(t+Δt₁)和端口电压U(t)，求得端口有功功率P_e(t+Δt₁)、无功功率Q_e(t+Δt₁)及电流I(t+Δt₁)，如下式所示：

式中，P_e和Q_e为微型燃气轮机机组注入电网的有功和无功功率，U为注入点的电压，δ为功角，

为注入电流的共轭，S_e为注入电网的视在功率，X_d为同步电机的d轴同步电抗，E_q为发电机空载电动势；

S3033：求解配电网模型，根据同步电机端口电流I(t+Δt₁)，更新电网节点注入电流矩阵，求解配电网模型得到电网中各节点的电压；

S3034：求解原动机模型，根据天然气管网末端出口流向微型燃气轮机的流量M_k(t_n)，求解差分化的微型燃气轮机模型，得到微型燃气轮机输出的机械功率P_m(t+Δt₁)以及天然气管网与微型燃气轮机连接点的压力p(t+Δt₁)；

S3035：求解P2G模型，根据P2G装置在电力能源系统拓扑连接模型接入点的电压，求解P2G模型得到其注入天然气管网中的流量。

S304：判断t时刻配电网模型是否有扰动发生，若有则修改电力能源系统拓扑连接模型的拓扑结构和动态参数，否则直接进入下一步；

S305：判断s是否大于或等于N+1，若是，进入下一步，否则跳转至步骤S303；

S306：令外部循环仿真时刻t_n＝t_n+Δt₂，将求得的天然气管网与微型燃气轮机连接点的压力p(t_n)作为天然气管网的末端压力，P2G注入天然气管网中的流量作为边界条件，以仿真步长Δt₂求解天然气管网模型，得到天然气管网系统各节点压力及各支路流量，并更新天然气管网流向微型燃气轮机的流量M_k(t_n)；

具体地，以仿真步长Δt_2求解所得的天然气管网模型为联立天然气管网系统的元件约束和拓扑约束后的系统方程：

式中，A为系数矩阵，b为常数向量，p表示各点的压力向量，M表示各支路的流量向量，其中，M又可以写为M＝[M_k；M_m]^T

式中，M_k表示管道的流入流量向量，M_m表示管道的流出流量向量。

S307：判断在t_n+Δt₂时刻天然气管网模型是否发生扰动，若有扰动，则根据扰动信息修改天然气能源系统拓扑连接模型的拓扑结构和动态参数，否则直接进行下一步骤；

S308：判断仿真时间t_n是否大于仿真终止时刻T，若是则进入下一步,否则跳转至步骤S302；

S309：输出仿真结果，仿真任务结束。

可以理解的是，所述仿真结果包括电气耦合拓扑连接模型以及求解所得的各状态参数。

结合上述步骤描述，下面给出具体算例：

本例以Matlab编程语言实现了用于电气耦合园区综合能源系统的多速率动态仿真方法，对本实施例提供的多速率动态仿真方法进行验证与分析，并与传统的单速率求解仿真方法仿真结果进行比较。仿真测试的硬件平台为Intel Core(TM)i5-6300HQ CPU@2.30GHz，8GB RAM的4核PC机；软件环境为64位Windows 10操作系统。

本算例中天然气经过50m管道后一部分再经100m管道流向一个恒定的天然气负荷，另一部分经过150m管道流向微型燃气轮机机组发电，配电网模型采用IEEE 33节点算例。

仿真场景设置如下：

设置仿真时间T＝1000s，仿真步长Δt₁＝0.001s，Δt₂＝0.05s，Δx＝25m，天然气流量的基值为50kg/s，600s时分别发生天然气负荷流量由50kg/s减少到45kg/s的扰动和电网25号节点阻抗标幺值增加0.1+j0.1的扰动，并将结果与单速率求解的结果进行对比，对比结果参见图3～图10。

需要说明的是，在传统的单速率求解步骤中，总是对配电网模型和天然气管网模型进行同步求解，即每次对配电网模型进行求解后，立刻使用求解所得的配电网模型对天然气管网模型进行求解，这样的做法可以提高求解精度，但是会导致运算次数较多，极大地降低了运算效率。本实施例采用的求解方法是，每对配电网模型进行N次求解才对天然气管网模型进行一次求解，这样虽然会损失一定的计算精度，但是却可以极大的减少运算次数，从而提高运算效率。针对上述参数进行对比，采用本实施例提供的多速率动态仿真方法所使用的运算时间为556.58s，用传统的单速率求解方法所使用的运算时间为827.22s，运算时间减少了大概30％。图3～图10为采用本实施例所提供的多速率动态仿真方法和采用传统的单速率求解方法求解所得的各状态参数的对比示意图，由图3～图10可知，采用本实施例所提供的多速率动态仿真方法后，各状态参数与采用传统的单速率求解方法求解所得的数值虽有差异，但是差异极小，均未达1％，但是运算时间却可以减少超过30％，其实际应用价值较高。

本实施例提供的用于电气耦合园区综合能源系统的多速率动态仿真方法，可以用于电气耦合技术领域中的电气耦合数学模型求解情形，减少天然气管网模型的求解次数，可以以较小的计算精度损失换取较大的计算速度提升，能有效提高仿真效率。

在本申请所提供的实施例中，应该理解到，所揭露的系统，单元，装置和方法，可以通过其它的方式实现。例如，以上所描述的所有实施例仅仅是示意性的，例如，上述单元或者模块等的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元、模块以及组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个计算机可读存储介质中，包括若干指令用以使得一台终端设备(可以是手机，笔记本，或者其他电子设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (1)

1.一种用于电气耦合园区综合能源系统的多速率动态仿真方法，其特征在于，包括：

建立电气耦合园区综合能源系统及其对应的电气耦合数学模型；

采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解；

所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型包括：

确定配电网中线路的阻抗及变压器参数，然后建立电力能源系统拓扑连接模型；

确定天然气管网中各支路的长度、直径及摩擦系数，然后建立天然气能源系统拓扑连接模型；

分别确定微型燃气轮机机组和P2G装置二者的位置、容量及动态参数，通过所述微型燃气轮机机组和P2G装置耦合所述天然气能源系统拓扑连接模型和电力能源系统拓扑连接模型，得到所述电气耦合园区综合能源系统对应的电气耦合拓扑连接模型；

针对所述电气耦合拓扑连接模型建立电气耦合数学模型；

所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型还包括：

指定所述电气耦合数学模型的仿真参数；所述仿真参数包括电气耦合数学模型的仿真终止时间T、配电网模型的仿真步长Δt₁、天然气管网模型的仿真步长Δt₂、天然气管网模型的空间步长Δx，且满足Δt₂＝NΔt₁，N为正整数，并设置仿真的扰动事件信息；

电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型，所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型之后，还包括：

将所述天然气管网模型和同步发电机模型转化为差分化形式；

所述采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解包括：

设定外部循环仿真时刻初值t_n＝0以及设定电气耦合数学模型中若干状态变量的初值；

设置配电网模型的求解次数s＝1且s为正整数；

设置当前时刻t＝t_n+sΔt₁，以t_n时刻天然气管网末端出口流量M_k(t_n)作为微型燃气轮机机组燃料系统的注入流量，以仿真步长Δt₁求解配电网模型，得到微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率、功角、天然气管网与微型燃气轮机连接点的压力p(t)、电网中各点的电压及P2G流向天然气管网中的流量；使s＝s+1；

判断t时刻配电网模型是否有扰动发生，若有则修改电力能源系统拓扑连接模型的拓扑结构和动态参数，否则直接进入下一步；

判断s是否大于或等于N+1，若是，进入下一步；

令外部循环仿真时刻t_n＝t_n+Δt₂，将求得的天然气管网与微型燃气轮机连接点的压力p(t_n)作为天然气管网的末端压力，P2G注入天然气管网中的流量作为边界条件，以仿真步长Δt₂求解天然气管网模型，得到天然气管网系统各节点压力及各支路流量，并更新天然气管网流向微型燃气轮机的流量M_k(t_n)；

判断在t_n+Δt₂时刻天然气管网模型是否发生扰动，若有扰动，则根据扰动信息修改天然气能源系统拓扑连接模型的拓扑结构和动态参数，否则直接进行下一步骤；

判断仿真时间t_n是否大于仿真终止时刻T，若是则进入下一步；

输出仿真结果，仿真任务结束。

Images