CN115906410A - 基于通用型辅助变量法的DET与cosSYST程序耦合的动态特性分析方法 - Google Patents

Abstract

本发明公开了基于通用型辅助变量法的DET与cosSYST程序耦合的动态特性分析方法，包括以下步骤：构建离散动态事件树的DET仿真模型，根据DET仿真对象状态转移类型，通过增加通用型辅助TRIP变量的方式，修改输入文件TRIP卡；通用型辅助TRIP变量为变量型或逻辑型；通过设定cosSYST程序的仿真时间，控制仿真时间步长，解析cosSYST程序每个仿真步长的输出结果文件；根据DET仿真对象状态转移类型，回溯执行cosSYST程序。本发明具有的优点是实现了DET与cosSYST程序动态耦合，简化了cosSYST程序输入文件中DET状态转移对象TRIP的设定流程和方法，减小了建模复杂度，提高了建模效率。

Description

基于通用型辅助变量法的DET与cosSYST程序耦合的动态特性分析方法

技术领域

本发明涉及一种核电厂安全分析领域，尤其涉及基于通用型辅助变量法的DET与cosSYST耦合的动态特性分析方法。

背景技术

传统的核电厂安全分析方法包括确定论安全分析方法和概率安全评价方法两种。使用确定论安全分析方法进行安全评价时，我国科研工作者自主开发了相关软件系统，即COSINE热工水力系统分析软件，该软件中的一个模块单元cosSYST程序是根据确定论实现安全分析评估的核心模块单元。概率安全评价方法基于传统事件树/故障树方法。单一的确定论或概率论方法在分析核电厂事故动态特性上存在局限性，例如法国电力集团开展的核电厂案例分析结果表明，传统的分析方法在估计由于核电厂严重事故现象导致的放射性释放频率上相当保守；因此，为了实现对复杂系统演化过程中的硬件/软件/流程/人类行为之间可能的复杂交互模型进行显式建模，进而发展了动态事件树方法；而且动态事件树方法与传统分析方法的主要不同之处在于核电厂事故物理进程由核电厂系统动态仿真演化结果驱动，减少了对分析人员经验的依赖，同时减少了不必要的保守的工程判断。

目前，离散动态事件树(DDET)是动态事件树方法在实际工程应用中应用较为广泛的一种方法(以下除特殊说明，DDET简称DET)，DET 方法根据其有效的分支规则生成一系列含时间动态变化的事件序列，具体通过确定分支条件的方式确定系统演化路径，即当条件满足时，就会导致复杂系统在不同分支路径上的演化，进而得到一系列事件序列集用以生成DET模型；DET模型生成后，所有由分支表示的事件序列都可以利用确定论安全分析方法进行仿真模拟，也即可开展确定论分析；通过DET与确定论分析软件的动态耦合，除了可以获得传统的确定论分析计算结果之外，还可获得计算分析结果的可能性也就是发生概率，进而为核电厂的设计、运行及管理提供风险指引型的决策支持信息。

传统地，基于核电厂确定论安全分析，离散动态事件树方法需要大量的计算资源并产生大量的数据，基于传统的分析人员手动设置分析边界条件的方法不能实现上述需求，因此，发展一种可自动实现DET方法与确定论分析软件cosSYST耦合的核电厂安全分析方法、技术方案和分析工具，进而实现在概率的环境下对核电厂的事故演化过程进行动态模拟，以获得核电厂的安全裕度特性以及安全性能。

发明内容

为了解决上述问题，克服传统的核电厂单一的确定论安全分析和概率安全评价的不足，本发明提供了基于通用型辅助变量法的DET与cosSYST 耦合的动态特性分析方法，通用性更强，适用范围更广，只需要区分DET 仿真对象状态转移属于参数型状态转移还是逻辑性状态转移两种即可；不需要在对仿真对象区分运行型还是需求型，一次转移还是多次转移，以及转移条件的逻辑判断关系。

本发明是通过以下技术方案实现的：基于通用型辅助变量法的DET 与cosSYST耦合的动态特性分析方法，包括以下步骤：

S1、构建离散动态事件树的DET仿真模型，根据DET仿真对象状态转移类型，通过增加通用型辅助TRIP变量的方式，修改输入文件的TRIP 卡；其中通用型辅助TRIP变量为变量型或逻辑型；

S2、确定DET仿真对象应用于cosSYST程序的仿真时间和仿真时间步长，并运行cosSYST程序；

S3、根据cosSYST程序输出仿真时间步长的结果文件，解析结果文件获取DET仿真对象状态转移控制TRIP变量变化的时间信息；

S4、识别DET仿真对象状态转移发生的DET分支时刻，并依据DET 仿真模型的分支规则获取cosSYST程序回溯再启动的cosSYST再启动号；

S5、确定DET仿真对象状态转移类型，更新cosSYST程序计算的再启动输入文件，并生成用于DET成功分支和失败分支的cosSYST程序计算的再启动输入文件；

S6、依据DET分支时刻确定cosSYST程序回溯再启动时间，并以更新后的成功和失败再启动文件为输入，结合当前仿真步长cosSYST计算结果文件，回溯执行cosSYST程序，完成DET成功分支和失败分支的下一个仿真步长的模拟；

S7、循环步骤S2-步骤S6的过程，一直到仿真时间达到用户指定的仿真任务时间，DET与cosSYST程序的动态耦合结束。

通过上述技术方案，通过增加辅助TRIP变量的方法，简化了cosSYST 程序输入文件中，DET仿真对象状态转移TRIP的设定流程和方法，降低了复杂性，提高了效率。

进一步地，在步骤S1中，以核电厂系统、设备或人员操作作为DET 仿真对象，DET仿真对象的状态转移通过cosSYST中的TRIP变量控制； TRIP变量逻辑值分为真(True)，假(False)；DET仿真对象的TRIP变量空间向量

其中n为DET仿真对象数量，s_i代表第i 个DET仿真对象对应的TRIP的逻辑值，即

进一步地，在步骤S2中，确定核电厂始发事故的分析对象，按照确定论分析方法编程模拟核电厂事故的初始cosSYST程序输入卡文件 input.dat和初始再启动计算文件restart.in；在第一个仿真步长结束之后， input.dat文件不再需要，其数据会记录在输出文件中的存储数据文件 restart.bin里。在后续的每次仿真中输入文件为再启动计算文件restart.in，其中cosSYST程序要求初始输入卡文件名称“input”和再启动计算文件名称“restart”不能更改，“.dat”和“.in”代表文件类型。

通过上述技术方案，该部分是cosSYST与DET耦合所必需的部分， restart.in文件中的DET仿真对象的TRIP变量需要按通用型辅助变量方式建模，使其内容上简化。

进一步地，在步骤S2中，所述cosSYST输入卡文件input.dat中除按照确定论分析方法编程模拟核电厂事故的数据之外，还需设置将所有DET 仿真对象的TRIP变量信息输出，初始再启动计算文件restart.in文件内容包括①、cosSYST程序processCtr数据块，用于设置再启动计算；②、cosSYST程序calcCtr数据块，用于设置计算时间步长，输出结果频率等信息；③、cosSYST程序trip数据块，用于设置逻辑与控制信息。

通过上述技术方案，该步骤是重启调用cosSYST程序必须的，该技术方案作为整体的一部分。

进一步地，在步骤S3和步骤S4中，识别DET仿真对象状态转移发生的DET分支时刻的具体方法如下：

第1步：根据步骤S3获得DET仿真对象状态转移控制TRIP变量变化的时间信息；

第2步：判断仿真时间步长是否存在分支；若不存在分支，则更新 cosSYST的再启动文件restart.in中的仿真时间，cosSYST程序输入卡中的仿真时间为结束时仿真时间与cosSYST程序的仿真步长之和，然后以更新后的再启动计算文件restart.in和结果文件restart.bin为输入，再次运行 cosSYST程序；

第3步：若存在分支，则通过解析结果文件获得DET仿真对象状态转移TRIP触发最小时间，然后根据DET仿真模型的分支规则获得与其 cosSYST回溯再启动信息块编号对应的cosSYST再启动号。

cosSYST程序作为封装的工具软件，在一个仿真时间步长内是不对外交互的，因此只能通过解析每一个步长的仿真输出结果，利用上述编号确定DET分支节点；该方法可以自动准确识别需要修改的DET分支点。

进一步地，在步骤S4中，所述DET回溯再启动号识别方法如下：

第(1)步：根据DET分支时刻获取的TRIP触发时间集合，识别其中最小触发时间TRIPTime_min以及TRIPTime_min对应TRIP变量比较并选择 TRIPTime_min对应的TRIP变量中最小TRIPNum_min的DET仿真对象状态转移进行分支；

第(2)步：判断TRIPTime_min对应的DET仿真对象是否为单次状态转移类型，若是则执行第(4)步；否则执行第(3)步；

第(3)步：本仿真步长DET仿真对象可能发生多次状态转移，解析 trip.plt文件信息，比较获得其中本步长DET分支仿真对象最小的状态转移时刻，并赋值给TRIPTime_min，然后执行第(4)步；

第(4)步：将TRIPTime_min替换{TRIPT1，TRIPT2，...，TRIPTN}中对应的元素，然后执行第(5)步。

第(5)步：计算DET回溯再启动信息块编号StepNo。

通过上述技术方案，可以自动、准确识别DET仿真对象状态转移发生的时刻。

更进一步地，所述StepNo是利用下述公式计算：

Δt＝CpuTime_max×rstStepFreq

其中：

StepNo：根据DET仿真对象状态转移的最小TRIP触发时间 TRIPTime_min计算的cosSYST回溯再启动信息块编号；

CpuTime_max：cosSYST程序输入卡中时间控制卡中的CPU最大仿真步长时间；

INT[]：取整函数；

rstStepFreq：cosSYST程序输入卡中时间控制卡的再启动频率；

TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最小TRIP触发时间；

Δt：cosSYST程序输入卡中时间控制卡中两个再启动信息块的时间间隔。

通过上述技术方案，可以自动、准确的根据DET仿真对象状态转移发生的时刻确定cosSYST再启动时刻。

进一步地，在步骤S3中，所述cosSYST再启动文件的更新方法是在每个cosSYST仿真步长结束，需要按照以下两类情况更新再启动文件：①、DET分支不存在DET仿真对象状态转移；②、DET分支存在DET仿真对象状态转移。

通过上述技术方案，其属于实际物理过程，按DET状态是否转移划分，结果就存在这两种可能，即成功和失效，不考虑部分成功的中间状态。

更进一步地，针对DET分支不存在DET仿真对象状态转移，cosSYST 再启动文件更新方法如下：

第1步：将当前节点cosSYST再启动文件中calcCtr数据块中的rstNo 再启动号更新为本次仿真步长仿真结果output.dat文件中最后一个再启动编号rstNoLast；

第2步：将当前节点cosSYST再启动文件中calcCtr数据块中的 endTime仿真结束时间按下列公式更新；；

endTime＝endTime+ΔT

其中：endTime：cosSYST程序输入卡中时间控制卡的仿真时间，也即当前节点结束时的仿真时间；ΔT：每次cosSYST程序的仿真步长。

通过上述技术方案，通过用户设置ΔT，可以方便控制cosSYST的仿真步长，并具备动态调整步长，进一步提高计算效率的可能。方法中是在每一个ΔT步长内识别DET状态对象转移，若客观上系统运行平稳，预计没有状态转移，则ΔT可以动态设定长一些，减少cosSYST停止、调用以及解析输出文件的时间，整体上提高运行效率。

更进一步地，针对DET分支存在DET仿真对象状态转移，cosSYST 再启动文件更新方法如下：

第1步：将当前节点cosSYST再启动文件restart.in中rstNo的再启动号更新为ResNum；

第2步：将当前节点cosSYST再启动文件restart.in中的endTime仿真结束时间按下列公式进行更新：

其中：INT[]：取整函数；TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最小TRIP触发时间；Δt：cosSYST程序输入卡中时间控制卡中两个再启动信息块的时间间隔；△T：每次cosSYST 程序的仿真步长；

第3步：分类更新DET仿真对象状态转移控制TRIP信息。

通过上述技术方案，自动根据再启动时间自动动态确定下一步 cosSYST仿真步长。

与现有技术相比，本发明具有以下优点：

(1)本发明基于离散动态事件树方法和cosSYST程序，给出了一种确定论与概率论动态耦合的核电厂事故实现分析方法，弥补传统安全分析方法在处理核电厂事故过程中系统、设备和操纵人员随机失效与核电厂动态过程耦合影响的不足，提高了耦合分析效率，减少对分析人员经验和判断的依赖，分析结果趋于现实；

(2)本发明可自动获得大量的核电厂事故仿真结果，从而获得核电厂事故工况下的安全特性，为基于风险指引的核电厂设计、运行和管理决策支持提供了支撑，优化核电厂设计、运行和管理，最终提高核电厂的安全性和经济性；

(3)本发明中的离散动态事件树方法在动态概率安全分析中普遍使用，cosSYST程序核电厂核安全分析领域广泛使用的商业软件，便于被核电工程和应用人员接受，有利于现实推广应用；

(4)本方法适用性广，可模拟核电厂事故分析过程中绝大多数核电厂系统、设备、人员操作状态转移、系统过程参数变化过程；

(5)本方法通用型强，可用于核电厂事故及瞬态的动态性能分析，也适用于可以利用cosSYST进行仿真模拟的其他复杂系统的性能分析。

附图说明

图1是本发明的结构示意图；

图2是核电厂状态演化的DET模型结构示意图；

图3是DET与cosSYST程序动态耦合示意图；

图4是DET仿真对象单次状态转移示意图；

图5是DET仿真对象多次状态转移示意图；

具体实施方式

下面结合附图对本发明做进一步的说明。

基于通用型辅助变量法的DET与cosSYST耦合的动态特性分析方法，如图1所示，包括以下步骤：

S1、构建离散动态事件树的DET仿真模型，根据DET仿真对象状态转移类型，通过增加通用型辅助TRIP变量的方式，修改输入文件的TRIP 变量；其中通用型辅助TRIP变量为变量型或逻辑型；

S2、确定DET仿真对象应用于cosSYST程序的仿真时间和仿真时间步长，并运行cosSYST程序；

S3、根据cosSYST程序输出仿真时间步长的结果文件，解析结果文件获取DET仿真对象状态转移控制TRIP变量变化的时间信息；

S4、识别DET仿真对象状态转移发生的DET分支时刻，并依据DET 仿真模型的分支规则获取cosSYST程序回溯再启动的cosSYST再启动号；

S5、确定DET仿真对象状态转移类型，更新cosSYST程序计算的再启动输入文件，并生成用于DET成功分支和失败分支的cosSYST程序计算的再启动输入文件；

S6、依据DET分支时刻确定cosSYST程序回溯再启动时间，并以更新后的成功和失败再启动文件为输入，结合当前仿真步长cosSYST计算结果文件，回溯执行cosSYST程序，完成DET成功分支和失败分支的下一个仿真步长的模拟；

S7、循环步骤S2-步骤S6的过程，一直到仿真时间达到用户指定的仿真任务时间，DET与cosSYST程序的动态耦合结束。

基于上述技术方案，步骤S1的具体实现原理如下：

首先以需要仿真模拟的核电厂系统、设备或人员操作为DET仿真对象，然后根据DET分析方法进行下列分析：

如图2所示，图中每个节点代表由DET仿真对象状态组合构成的一个核电厂特定状态S_i，其中根节点代表DET与cosSYST动态耦合仿真时核电厂的初始状态S₀，中间节点代表核电厂系统、设备不同状态的组合，终节点代表预先定义的核电厂某一个吸收终止状态，如达到核电厂堆芯损毁或稳定状态；

DET与cosSYST动态耦合时将从根节点开始，随着时间变化，DET 仿真对象状态将发生随机性转移；在达到预设的分支规则时，从该时刻产生的“分支”代表仿真对象的一次状态转移；且从父节点到子节点代表仿真对象的一次状态转移；

当DET仿真对象状态转移触发条件满足时，DET仿真对象状态转移仅存在状态转移成功和状态转移失败的两种状态，父节点对应DET模型产生的新的分支将演化为成功节点和失败节点，其中上节点代表成功节点，节点编码方式是“父节点名-1；下节点代表转移失败节点，节点编码方式是“父节点名称-2”；终节点代表DET与cosSYST耦合的终止，节点名称编码方式是“父节点名称-f”，并且从父节点到子节点之间的系统演化过程由cosSYST程序的1个或多个仿真步长计算。

在上述方案的基础上，如图2所示，当DET与cosSYST程序耦合时，就是在每个节点之间进行计算；与传统的事件树方法不同，DET中节点分支发生的时间是由cosSYST仿真的物理过程决定的，减小了传统方法中对分析人员经验依赖或保守的工程判断。

实施例

基于上述技术方案，步骤S2-步骤S7的具体实现过程如下：

如图3所示，针对图中所示符号定义如下：

StepNo：根据DET仿真对象状态转移的最小TRIP触发时间 TRIPTimemin计算的cosSYST回溯再启动信息块编号；

CpuTime_max：cosSYST程序输入卡中时间控制卡中的CPU最大仿真步长时间；

ErrNum：cosSYST异常终止计数；

MT：运行型状态转移DET仿真对象的状态转移任务时间；

NorNum：cosSYST正常时间终止计数；

rstStepFreq：cosSYST程序输入卡中时间控制卡的再启动频率；

ResNum：与cosSYST回溯重启信息块编号StepNo对应的cosSYST 再启动号；

SimTime：cosSYST程序输入卡中时间控制卡的仿真时间，也即当前节点结束时的仿真时间；

TaskTime：DET与cosSYST耦合仿真的任务时间，也是cosSYST仿真最终终止的仿真终止时间；

TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最小TRIP触发时间；

rstNoLast：当前cosSYST仿真步长计算结果文件中最后一个再启动编号；

△t：cosSYST程序输入卡中时间控制卡中的重启频率；

△T：每次cosSYST程序的仿真步长。

结合图3说明DET与cosSYST程序动态耦合的总体过程如下：

第1步：编写DET与cosSYST耦合仿真事故分析输入文件input.da 和初始再启动文件restart.in文件；

第2步：根据DET仿真对象状态转移类型，通过增加辅助变量的方法，修改restart.in中的TRIP信息，形成新的restart.in输入文件；

第3步：在DET根节点，以cosSYST事故分析输入文件input.dat为初始输入文件，以SimTime＝△T作为仿真时间，运行cosSYST程序；

第4步：在cosSYST仿真步长结束以后，生成结果文件分别为trip.plt，output.dat和restart.bin；解析trip.plt文件中的DET仿真对象状态控制TRIP 变量变化的时间信息；

第5步：判断本次仿真步长是否存在分支；若不存在分支，则更新 cosSYST初始的再启动文件restart.in中的仿真时间，SimTime＝SimTime+ △T，然后以restart.bin和restart.in为输入，再次运行cosSYST分析；

第6步：若存在分支，则通过解析output.dat文件获得DET仿真对象状态转移TRIP触发最小时间TRIPTime_min，然后根据规则获得cosSYST程序回溯再启动编号ResNum；

第7步：根据DET仿真对象状态转移类别，更新restart.in，分别生成 DET成功分支cosSYST再启动输入文件R0-1和DET失败分支cosSYST 再启动输入文件R0-2；

第8步：分别以成功分支再启动文件和失败分支再启动文件，以及当前节点cosSYST计算输出的restart.bin文件作为输入，分别运行cosSYST， cosSYST的节点输出文件名称根据规则更新；

第9步：循环第2步-第8步的过程，一直到节点的cosSYST仿真时间SimTime达到仿真任务时间TaskTime或者达到cosSYST其他终止条件。

在上述方案的基础上，在第1步中，确定核电厂始发事故的分析对象，按照确定论分析方法编程模拟核电厂事故的初始cosSYST程序输入卡文件，其中输入卡文件包括初始cosSYST输入文件input.dat和初始再启动计算文件restart.in；并且始发事故cosSYST仿真模型需调试稳定，典型的事故序列分析结果符合预期；input.dat调试稳定后，将restart.in输入卡中的processCtr数据块中的问题类型设置为restartProblem；calcCtr数据块中的结束时间初始化为△T。

在第1步中，restart.in文件内容包括①、cosSYST程序processCtr数据块，用于设置再启动计算；②、cosSYST程序calcCtr数据块，用于设置计算时间步长，输出结果频率等信息；③、cosSYST程序trip数据块，用于设置逻辑与控制信息。具体地，以图2结合图3中根节点第1次DET 分支为例：

第①步：根据第3步获得DET仿真对象状态转移控制TRIP变量变化的时间信息；

第②步：判断仿真时间步长是否存在分支；若不存在分支，则更新 cosSYST初始的再启动文件restart.in中的仿真时间，cosSYST程序输入卡中时间控制卡的仿真时间为结束时仿真时间与cosSYST程序的仿真步长之和，然后以初始再启动计算文件restart.in和结果文件restart.bin为输入，再次运行cosSYST程序；

第③步：若存在分支，则通过解析结果文件获得DET仿真对象状态转移TRIP触发最小时间，然后根据DET仿真模型的分支规则获得与其 cosSYST回溯再启动信息块编号对应的cosSYST再启动号。

上述判断仿真时间步长是否存在DET分支的判断和识别方法如下：

第(1)步：初始化DET仿真对象状态转移TRIP触发时间集合变量 {TRIPT1，TRIPT2，…,TRIPT N}和当前节点TRIP触发事件集合变量 {NowTRIPT1，NowTRIPT 2，…,NowTRIPTN}，每个触发时间为“-1”。其中N代表与DET仿真对象状态转移相关的TRIP逻辑数量；

第(2)步：打开cosSYST输出的TRIP信息文件，提取输出文件中的所有与DET仿真对象状态转移相关的TRIP变量的触发时间，获得当前仿真步长DET仿真对象状态转移TRIP触发时间集合{NowTRIPT1， NowTRIPT 2，…,NowTRIPT N}；

第(3)步：逐一对应比较{TRIPT1，TRIPT 2，…,TRIPT N}和 {NowTRIPT1，NowTRIPT2，…,NowTRIPTN}的元素，获得TRIP触发时间变化的信息集合{NowTRIPTi}，其中i∈(1，N)；若{NowTRIPTi}为空，则本步长不发生DET仿真对象状态转移，不存在DET分支，本步长仿真DET分支识别结束。若{NowTRIPTi}不为空，则在本步长仿真中存在DET 仿真对象状态转移过程，存在DET分支，本步长DET分支识别结束。

在上述方案的基础上，在第2步中，通过如下方式设置通用型辅助 TRIP变量；具体地，通用型辅助TRIP变量分参数型和逻辑型状态转移对象的DET与cosSYST耦合方法，并且其分别按单次状态转移和多次状态转移两种类别的具体设置方法如下：

(1)参数型状态转移TRIP辅助变量设置方法

cosSYST程序中，需求型单次状态转移辅助TRIP修改方式如下： TRIP信息：

"CC1"控制DET仿真对象状态转移的判断逻辑，即初始TRIP逻辑

"TF-T"：＝time>-1.0

辅助TRIP变量，其逻辑值是一个“真”的常量

"TF-F"：＝time<-1.0

辅助TRIP变量，其逻辑值是一个“假”的常量。

1、检测到CC1为True且状态将发生改变，则其分支逻辑如下：

状态不变："CC1"＝time>-1.0

状态转移："CC1"继承初始TRIP逻辑

2、检测到CC1为False且状态将发生改变，则其分支逻辑如下：

状态不变："CC1"＝time<-1.0

状态转移："CC1"继承初始TRIP逻辑

其中：

CC1：初始cosSYST的输入文件中满足唯一性、独立性要求的字符串均可，根据核电厂事故特性设置的控制DET仿真对象正常需求状态转移的TRIP逻辑。本发明中讨论了CC1初始逻辑值为“假”和“真”两种状态下，在不同的初始状态条件下满足上述所示的不同的状态转移规则。

当DET分支过程中检测到状态转移恢复信号后，将状态转移后的CC1 逻辑重置为转移前的控制逻辑，等待后续的参数型分支信号触发并按上述的规则进行分支，即可实现CC1的多次状态转移。

(2)逻辑状态转移辅助TRIP变量设置方法

cosSYST程序中，逻辑状态转移辅助TRIP卡修改按如下设置。

辅助变量永真TRIP："TF-T"＝time>-1.0

辅助变量永假TRIP："TF-F"＝time<-1.0

CC1单次需求状态转移DET仿真对象的TRIP逻辑

辅助变量"TFN"TRIP："TFN"＝原"CC1"TRIP逻辑

1、检测到CC1为True且状态将发生改变，则其分支逻辑如下：

状态不变："CC1"＝"TFN"&&"TF-F"

状态转移："CC1"继承初始TRIP逻辑

2、检测到CC1为False且状态将发生改变，则其分支逻辑如下：

状态不变："CC1"＝"TFN"&&"TF-T"

状态转移："CC1"继承初始TRIP逻辑

其中：

CC1：上述所控制DET仿真对象可进行1次状态转移的TRIP逻辑，当检测到TRIP状态转移信号时，当初始状态为True时执行上述的第一个状态转移规则，当初始状态为False时执行上述中第二个状态转移规则，状态发生转移后其值便不再变化。

"TF-T"：辅助TRIP变量，其逻辑值是一个“真”的常量；

"TF-F"：辅助TRIP变量，其逻辑值是一个“假”的常量；

"TFN"：辅助TRIP变量，其逻辑值是一个与CC1逻辑值完全相同的常量；

当DET分支过程中检测到状态转移恢复信号后，将状态转移后的CC1 逻辑重置为转移前的控制逻辑，等待后续分支信号继续执行上述的分支规则即可实现CC1的多次状态转移。

在上述方案的基础上，以图2结合图3中cosSYST输出文件第2个仿真步长，也就是2△T结束时的分支，说明DET回溯再启动号识别方法如下：

第(1)步：根据DET分支时刻获取的TRIP触发时间集合，识别其中最小触发时间TRIPTime_min以及TRIPTime_min对应的TRIP变量，比较并选择TRIPTime_min对应的TRIP变量中最小的DET仿真对象状态转移进行分支；

第(2)步：判断TRIPTime_min对应的DET仿真对象是否为单次状态转移类型，若是则执行第(4)步；否则执行第(3)步；

第(3)步：本仿真步长DET仿真对象可能发生多次状态转移，解析 trip.plt文件中的信息，比较获得其中本步长DET分支仿真对象最小的状态转移时刻，并赋值给TRIPTime_min，然后执行第(4)步；

第(4)步：将TRIPTimemin替换{TRIPT1，TRIPT 2，…,TRIPT N} 中对应的元素，然后执行第(5)步。

第(5)步：在TRIP list信息块信息中获取DET回溯再启动信息块编号StepNo。

在上述方案的基础上，StepNo是利用下述公式计算：

Δt＝CpuTime_max×rstStepFreq

其中：

StepNo：根据DET仿真对象状态转移的最开始TRIP触发时间 TRIPTime_min计算的cosSYST回溯再启动编号；

CpuTime_max：cosSYST程序输入卡中时间控制卡中的CPU最大仿真步长时间；

INT[]：取整函数；

rstStepFreq：cosSYST程序输入卡中时间控制卡的再启动频率；

TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最开始TRIP触发时间；

Δt：cosSYST程序输入卡中时间控制卡中两个再启动信息块的时间间隔。

通过上述技术方案，可以自动、准确的根据DET仿真对象状态转移发生的时刻确定cosSYST再启动时刻。

通过上述公式，可以确定一次仿真步长过程中，比DET仿真对象状态转移时刻小，并且最近的cosSYST再启动编号StepNo；通过本仿真步长cosSYST输出文件中的trip.pit解析，即可查找到与StepNo对应的 cosSYST回溯计算需要的再启动号ResNum；利用ResNum更新cosSYST 再启动卡中的再启动编号，结合当前节点cosSYST的restart.bin文件，如图3中的R0节点的restart.bin文件，就可以在ResNum对应的时间点继续运行cosSYST程序进行仿真。

在上述方案的基础上，在第4步中，cosSYST再启动文件的更新方法是在每个cosSYST仿真步长结束后，需要按照以下两类情况更新再启动文件：①、DET分支不存在DET仿真对象状态转移；②、DET分支存在 DET仿真对象状态转移。

具体地，针对DET分支不存在DET仿真对象状态转移，cosSYST再启动文件更新方法如下：

第1步：将当前节点cosSYST再启动文件中的calcCtr数据块中的 rstNo再启动号更新为本次仿真步长仿真结果output.dat文件中最后一个再启动编号rstNoLast；

第2步：将当前节点cosSYST再启动文件中的calcCtr数据块中的 endTime仿真结束时间按下列公式更新；

endTime＝endTime+ΔT

其中：endTime：cosSYST程序输入卡中时间控制卡的仿真时间，也即当前节点结束时的仿真时间；ΔT：每次cosSYST程序的仿真步长。

具体地，针对DET分支存在DET仿真对象状态转移，cosSYST再启动文件更新方法如下：

第1步：将当前节点cosSYST再启动文件restart.in中rstNo再启动号更新为ResNum；

第2步：将当前节点cosSYST再启动文件restart.in中endTime仿真结束时间按下列公式进行更新：

其中：INT[]：取整函数；TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最开始的trip触发时间；Δt：cosSYST程序输入卡中时间控制卡中两个再启动信息块的时间间隔；ΔT：每次 cosSYST程序的仿真步长；

第3步：分类更新DET仿真对象状态转移控制TRIP信息。

在上述方案的基础上，通用型TRIP变量为变量型或逻辑性变量设置，具体设置内容如下：

若DET仿真对象状态转移控制TRIP属于参数型，则按如下方式更新TRIP信息：

1、检测到CC１为True且状态将发生改变，则其分支逻辑如下：

状态不变：″CC1″＝time＞-1.0

状态转移："CC1"继承初始TRIP逻辑

2、检测到CC1为False且状态将发生改变，则其分支逻辑如下：

状态不变："CC1"＝time<-1.0

状态转移："CC1"继承初始TRIP逻辑

若识别DET仿真对象初始状态为true(真)状态，则上下分支再启动文件按如下方式更新：

上分支(成功分支、状态转移分支)：再启动文件restart.in不更新。也即实际过程中程序将按照既定的转移规则发生状态转移；

下分支(失败分支、状态不转移分支)：在再启动文件restart.in中按 cosSYST规则在TRIP数据块中修改原TRIP信息如下：

"CC1"＝time>-1.0

其中CC1是初始DET仿真对象状态转移控制的TRIP名。

若识别DET仿真对象初始状态为false(假)状态，则上下分支再启动文件按如下方式更新：

上分支(成功分支、状态转移分支)：再启动文件restart.in不更新。也即实际过程中程序将按照既定的转移规则发生状态转移；

下分支(失败分支、状态不转移分支)：在再启动文件restart.in中按cosSYST规则在TRIP数据块中修改原TRIP信息如下：

"CC1"＝time<-1.0

其中CC1是初始DET仿真对象状态转移控制的TRIP名。

若DET仿真对象状态转移控制TRIP属于逻辑性型，则按如下方式更新TRIP信息：

辅助变量永真TRIP："TF-T"＝time>-1.0

辅助变量永假TRIP："TF-F"＝time<-1.0

CC1单次需求状态转移DET仿真对象的TRIP逻辑

辅助变量"TFN"TRIP："TFN"原"CC1"TRIP逻辑

1、检测到CC1为True且状态将发生改变，则其分支逻辑如下：

状态不变："CC1"＝"TFN"&&"TF-F"

状态转移："CC1"继承初始TRIP逻辑

2、检测到CC1为False且状态将发生改变，则其分支逻辑如下：

状态不变："CC1"＝"TFN"&&"TF-T"

状态转移："CC1"继承初始TRIP逻辑

若识别DET仿真对象初始状态为true(真)状态，则上下分支再启动文件按如下方式更新：

上分支(成功分支、状态转移分支)：再启动文件restart.in不更新。也即实际过程中程序将按照既定的转移规则发生状态转移；

下分支(失败分支、状态不转移分支)：在再启动文件restart.in中按 cosSYST规则在TRIP数据块中修改原TRIP信息如下：

"CC1"＝"TFN"&&"TF-T"

其中

CC1：是初始DET仿真对象状态转移控制的TRIP名；

TFN：任意辅助变量TRIP名，该TRIP名符合cosSYST要求，并满足独立性和唯一性即可。

TF-T：辅助变量TRIP名，逻辑值为永“真”

若识别DET仿真对象初始状态为false(假)状态，则上下分支再启动文件按如下方式更新：

上分支(成功分支、状态转移分支)：再启动文件restart.in不更新。也即实际过程中程序将按照既定的转移规则发生状态转移；

下分支(失败分支、状态不转移分支)：在再启动文件restart.in中按 cosSYST规则在TRIP数据块中修改原TRIP信息。

"CC1"＝"TFN"&&"TF-F"

CC1：是初始DET仿真对象状态转移控制的TRIP名；

TFN：任意辅助变量TRIP名，该TRIP名符合cosSYST要求，并满足独立性和唯一性即可。

TF-F：辅助变量TRIP名，逻辑值为永“假”

其中：TF-T是cosSYST软件中用户可以自定义一个TRIP变量名，是为便于实现DET与cosSYST程序动态耦合而额外增加的一个辅助变量，其逻辑值是永“真”。在实际应用过程中可根据需要修改，符合唯一性、独立性要求的字符串均可；

TF-F是cosSYST软件中用户可以自定义一个TRIP变量名，是为便于实现DET与cosSYST程序动态耦合而额外增加的一个辅助变量，其逻辑值是永“假”。在实际应用过程中可根据需要修改，符合唯一性、独立性要求的字符串均可；

效果实施例

根据事故特征和仿真计算资源条件，初始化DET与cosSYST动态耦合相关参数，包括：

StepNo：根据DET仿真对象状态转移的最小TRIP触发时间 TRIPTime_min计算的cosSYST回溯再启动编号，初始值为0；

CpuTime_max：cosSYST程序输入卡中时间控制卡中的CPU最大仿真步长时间，初始值用户自定义；

ErrNo：cosSYST异常终止计数，初始值为0；

MT：运行型状态转移DET仿真对象的状态转移任务时间，根据涉及到的运行型状态转移对象的数量确定需要的任务时间随机变量的数量，并根据用户自定义的离散方式对每个任务时间随机变量进行离散，获得离散的任务时间序列；

NorNo：cosSYST正常时间终止计数；

rstStepFreq：cosSYST程序输入卡中时间控制卡的再启动频率，用户自定义；

ResNum：与cosSYST回溯重启信息块编号StepNo对应的cosSYST 再启动号，初始值为0；

SimTime：cosSYST程序输入卡中时间控制卡的仿真时间，也即当前节点结束时的仿真时间，初始值为△T；

TaskTime：DET与cosSYST耦合仿真的任务时间，也是cosSYST仿真最终终止的仿真终止时间，初始值用户自定义；

TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最小TRIP触发时间，初始值为0；

rstNoLast：当前cosSYST仿真步长计算结果文件中最后一个再启动编号，初始值为0；

△t：cosSYST程序输入卡中时间控制卡中两个再启动信息块的时间间隔；

△T：每次cosSYST程序的仿真步长；

如图3所示，DET与cosSYST程序动态具体耦合过程：

第1步：将DET根节点R0的初始cosSYST输入文件input.dat和初始再启动文件restart.in作为一对文件放入计算队列中；

第2步：判断计算队列中是否存在cosSYST输入文件。若存在 cosSYST输入文件，则按照“后进先出”的原则，选择一个cosSYST输入文件作为cosSYST程序的输入，假定为R0的输入文件input.dat，运行 cosSYST仿真，获得R0的输出结果output文件夹；若计算队列中不存在 cosSYST输入文件，则执行第5步；

第3步：以output文件夹中的trip.plt为输入，利用DET分支判别方法，判断是否存在DET分支。若不存在DET分支，则按照没有DET分支的类别，更新本次计算结果对应的上一个cosSYST仿真步长的再启动文件restart.in。将更新后的再启动文件restart.in和对应的cosSYST的 restart.bin文件作为一个计算案例放入计算队列中；然后执行第2步；

若存在分支，则执行第4步；

第4步：根据DET仿真对象类别，生成DET成功分支和失败分支，分别命名为R0-1和R0-2，按照再启动文件更新方法生成DET成功分支 R0-1的再启动文件restart.in和失败分支R0-2的再启动文件restart.in，并分别与cosSYST的restart.bin文件组成两个计算案例放入计算队列中，然后执行第2步；

第5步：每个DET分支都达到用户指定的仿真任务时间或cosSYST 仿真终止条件，DET与cosSYST程序的动态耦合结束；

综上所述，DET与cosSYST程序动态耦合以后，可以通过解析每个分支节点及终节点的输出文件，提取不同分支序列的核电厂的物理特性；可以通过统计分析、聚类分析等，获得核电厂的安全特性等。

以上所述实施方式仅表达了本发明的一种或多种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。

Claims (2)

1.基于通用型辅助变量法的DET与cosSYST程序耦合的动态特性分析方法，其特征在于，包括以下步骤：

S1、构建离散动态事件树的DET仿真模型，根据DET仿真对象的状态转移类型，通过增加通用型辅助TRIP变量方式修改cosSYST程序输入文件的TRIP数据块内容；其中通用型辅助TRIP变量为变量型或逻辑型；这里的DET仿真对象为核电厂系统、设备或人员操作，DET仿真对象的状态转移通过cosSYST中的TRIP变量控制；TRIP变量的逻辑值分为真True，假False；DET仿真对象的TRIP变量的空间向量

其中n为DET仿真对象数量，s_i代表第i个DET仿真对象所对应的TRIP逻辑值；即

S2、确定DET仿真对象应用于cosSYST程序的仿真时间和仿真时间步长，并运行cosSYST程序；此时的DET仿真对象为核电厂始发事故的分析对象，按照确定论分析方法编程模拟核电厂事故的初始cosSYST程序输入卡文件input.dat和初始再启动计算文件restart.in，在第一个仿真步长结束之后，input.dat文件不再需要，其数据会记录在输出文件restart.bin里，在后续的每次仿真中输入文件为再启动计算文件restart.in，其中cosSYST程序要求初始输入卡文件名称“input”和再启动计算文件名称“restart”不能更改，“.dat”和“.in”代表文件类型；

S3、根据cosSYST程序输出仿真时间步长的结果文件，解析结果文件获取DET仿真对象状态转移控制TRIP变量变化的时间信息；

S4、识别DET仿真对象状态转移发生的DET分支时刻，并依据DET仿真模型的分支规则获取cosSYST程序回溯再启动的cosSYST再启动号；上述识别DET仿真对象状态转移发生的DET分支时刻的具体方法如下：第①步：根据步骤S3获得DET仿真对象状态转移控制TRIP变量变化的时间信息；第②步：判断仿真时间步长是否存在分支；若不存在分支，则更新cosSYST程序初始的再启动文件restart.in中的仿真时间，cosSYST程序输入卡中时间控制卡的仿真时间为结束时仿真时间与cosSYST程序的仿真步长之和，然后以更新后的再启动计算文件restart.in和结果文件restart.bin为输入，再次运行cosSYST程序；第③步：若存在分支，则通过解析结果文件获得DET仿真对象状态转移TRIP触发最小时间，然后根据DET仿真模型的分支规则获得与其cosSYST回溯重启信息块编号对应的cosSYST再启动号；这里的DET回溯再启动号识别方法如下：第(1)步：根据DET分支时刻获取的TRIP触发时间集合，识别其中最小触发时间TRIPTime_min以及TRIPTime_min对应TRIP变量，比较并选择TRIPTime_min对应的TRIP变量中最小TRIPNum_min的DET仿真对象状态转移进行分支；第(2)步：判断TRIPTime_min对应的DET仿真对象是否为单次状态转移类型，若是则执行第(4)步；否则执行第(3)步；第(3)步：本仿真步长DET仿真对象发生多次状态转移，解析trip.plt文件中所有的TRIP信息块信息，比较获得其中本步长DET分支仿真对象最小的状态转移时刻，并赋值给TRIPTime_min，然后执行第(4)步，这里的trip.plt文件是cosSYST输出结果output文件夹中的一个文件，里面包含了在input.dat文件中设置输出的所有TRIP变量的信息，该文件名“trip”可以在input.dat文件中设置任意名称符合cosSYST程序要求即可，“.plt”代表文件类型；第(4)步：将TRIPTime_min替换{TRIPT 1，TRIPT 2，…,TRIPT N}中对应的元素，然后执行第(5)步；第(5)步：计算DET回溯重启信息块编号StepNo，这里的StepNo是利用下述公式计算：

Δt＝CpuTime_max×rstStepFreq

其中：StepNo：根据DET仿真对象状态转移的最小TRIP触发时间TRIPTime_min计算的cosSYST回溯重启信息块编号；CpuTime_max：cosSYST程序输入卡中时间控制卡中的CPU最大仿真步长时间；INT[]：取整函数；rstStepFreq：cosSYST程序输入卡中时间控制卡的再启动频率；TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最小TRIP触发时间；△t：cosSYST程序输入卡中时间控制卡中两个再启动信息块的时间间隔；

S5、确定DET仿真对象状态转移类型，更新cosSYST程序计算的再启动输入文件，并生成用于DET成功分支和失败分支的cosSYST程序计算的再启动输入文件；这里的cosSYST再启动文件的更新方法是在每个cosSYST仿真步长结束，需要按照以下两类情况更新再启动文件：①DET分支不存在DET仿真对象状态转移，此时的cosSYST再启动文件更新方法为：第1步：将当前节点cosSYST再启动文件中的calcCtr数据块的rstNo再启动号更新为本次仿真步长仿真结果output.dat文件中最后一个再启动编号rstNoLast；第2步：将当前节点cosSYST再启动文件的calcCtr数据块中的endTime仿真结束时间按下列公式更新；SendTime＝endTime+ΔT；其中：endTime：cosSYST程序输入卡中时间控制卡的仿真时间，也即当前节点结束时的仿真时间；△T：cosSYST程序的仿真步长；②DET分支存在DET仿真对象状态转移；此时的cosSYST再启动文件更新方法如下：第1步：将当前节点cosSYST再启动文件restart.in中rstNo再启动号更新为ResNum；第2步：将当前节点cosSYST再启动文件restart.in中endTime仿真结束时间按下列公式进行更新：

其中：INT[]：取整函数；TRIPTime_min：当前cosSYST仿真步长计算结果中DET仿真对象状态转移最小TRIP触发时间；△t：cosSYST程序输入卡中时间控制卡中两个再启动信息块的时间间隔；△T：每次cosSYST程序的仿真步长；第3步：分类更新DET仿真对象状态转移控制TRIP信息；

S6、依据DET分支时刻确定cosSYST程序回溯再启动时间，并以更新后的成功和失败再启动文件为输入，结合当前仿真步长cosSYST计算结果文件，回溯执行cosSYST程序，完成DET成功分支和失败分支的下一个仿真步长的模拟；

S7、循环步骤S2-步骤S6的过程，一直到仿真时间达到用户指定的仿真任务时间，DET与cosSYST程序的动态耦合结束。

2.根据权利要求1所述的基于通用型辅助变量法的DET与cosSYST耦合的动态特性分析方法，其特征在于，在步骤S2中，所述cosSYST输入卡文件input.dat中除按照确定论分析方法编程模拟核电厂事故的数据之外，还需设置将所有TRIP变量信息输出，初始再启动计算文件restart.in文件内容包括①cosSYST程序processCtr数据块，用于设置再启动计算；②cosSYST程序calcCtr数据块，用于设置计算时间步长，输出结果频率等信息；③cosSYST程序TRIP数据块，用于设置TRIP逻辑与控制信息。

Images