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

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

【技术实现步骤摘要】
基于通用型辅助变量法的DET与cosSYST程序耦合的动态特性分析方法


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

技术介绍

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

【技术保护点】

【技术特征摘要】
1.基于通用型辅助变量法的DET与cosSYST程序耦合的动态特性分析方法，其特征在于，包括以下步骤：S1、构建离散动态事件树的DET仿真模型，根据DET仿真对象的状态转移类型，通过增加通用型辅助TRIP变量方式修改cosSYST程序输入文件的TRIP数据块内容；其中通用型辅助TRIP变量为变量型或逻辑型；这里的DET仿真对象为核电厂系统、设备或人员操作，DET仿真对象的状态转移通过cosSYST中的TRIP变量控制；TRIP变量的逻辑值分为真True，假False；DET仿真对象的TRIP变量的空间向量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，

【专利技术属性】
技术研发人员：陈浩尹，熊彬富，王贺，沙会娥，孟召灿，全国萍，余慧，赵强，李磊，
申请(专利权)人：国家电投集团科学技术研究院有限公司，
类型：发明
国别省市：