CN111881176A - Anomaly detection method for marine nuclear power based on logical distance characterization of safe operation domain - Google Patents
Anomaly detection method for marine nuclear power based on logical distance characterization of safe operation domain Download PDFInfo
- Publication number
- CN111881176A CN111881176A CN202010646535.1A CN202010646535A CN111881176A CN 111881176 A CN111881176 A CN 111881176A CN 202010646535 A CN202010646535 A CN 202010646535A CN 111881176 A CN111881176 A CN 111881176A
- Authority
- CN
- China
- Prior art keywords
- operating
- logical distance
- nuclear power
- data
- state
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2465—Query processing support for facilitating data mining operations in structured databases
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/284—Relational databases
- G06F16/285—Clustering or classification
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Linguistics (AREA)
- Probability & Statistics with Applications (AREA)
- Fuzzy Systems (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明属于核动力系统异常检测技术领域,公开了一种基于逻辑距离表征安全运行域的船用核动力异常检测方法,通过分析核动力系统历史运行数据,检索系统典型运行工况数据样本;构建标准运行工况样本库,遴选表征运行状态的特征参数;构建逻辑距离计算函数;通过历史运行数据进行训练学习、仿真计算或收集的异常数据样本的检验验证确定判定系统异常的逻辑距离阈值即安全域归属阈值,在线检测系统异常或检索系统历史运行数据中的异常情况。本发明直接从对象核动力装置运行状态数据出发构造异常检测模型,贴近实际系统,从机理角度比其他算法具有更强的针对性和适用性,具有可解释性强,检测力度和调节等优势。
The invention belongs to the technical field of nuclear power system abnormality detection, and discloses a marine nuclear power abnormality detection method based on a logical distance to represent a safe operation domain. The operating condition sample library is used to select characteristic parameters that characterize the operating state; build a logical distance calculation function; conduct training learning, simulation calculation, or test and verification of collected abnormal data samples through historical operating data to determine the logical distance threshold for judging system abnormality, that is, the safety domain Attribution threshold, online detection of system anomalies or retrieval of anomalies in the historical operating data of the system. The invention constructs an abnormality detection model directly from the operating state data of the target nuclear power plant, which is close to the actual system, has stronger pertinence and applicability than other algorithms from the perspective of mechanism, and has the advantages of strong interpretability, detection strength and adjustment.
Description
技术领域technical field
本发明属于核动力系统异常检测技术领域,尤其涉及一种基于逻辑距离表征安全运行域的船用核动力异常检测方法。The invention belongs to the technical field of abnormal detection of nuclear power systems, and in particular relates to a method for abnormal detection of marine nuclear power based on logical distance representation of a safe operation domain.
背景技术Background technique
目前,核动力系统异常状态自动检测技术是系统故障诊断技术的前端入口程序模块,也是辅助操纵员值班的伺服技术模块。现有国内外关于核动力系统异常检测和故障诊断技术的研究大致分为两大类:机理分析法和数据驱动法。At present, the automatic detection technology of abnormal state of nuclear power system is the front-end entry program module of system fault diagnosis technology, and it is also the servo technology module that assists operators on duty. The existing domestic and foreign research on abnormal detection and fault diagnosis technology of nuclear power system can be roughly divided into two categories: mechanism analysis method and data-driven method.
1)机理分析异常检测法。机理分析法发展较早也比较受青睐,是当前核动力系统异常检测和故障诊断的传统主流方法,其又分为模拟仿真法和专家规则集导向法。模拟仿真法是指依据核动力装置结构特点,利用物理、热工等方面的守恒方程和经验公式,构建与核动力系统基本相同的数值计算全范围仿真模型系统,在模型系统上进行异常、故障和事故演示,获得规律和经验,以指导实际系统的异常检测和故障诊断,其应用被拓展到了核动力系统设计验证和核动力操纵员培训中去,成为了目前核动力领域主流的工程设计验证和人员培训方法,然而应用到异常检测和故障诊断中,该方法存在与实际系统同步困难的问题;专家规则集法是指综合利用物理、热工方程和专家经验,设定事故辨识程序,制定事故处置规程,用于伺服检测系统异常、诊断系统故障、指导事故处置的技术方法,是直接作用在核动力装置上的异常检测和故障诊断方法,也是目前国内外核动力装置普遍采用的一种方法,然而其异常检测和事故辨识程序流程复杂,需要大量人员干预,很难保证其完备性,且容易诱发人因失误,而其诊断和处置过程要么固定死板无法权变,要么含糊不清不易理解,所以近几年来许多国内外专家提出了基于征兆或状态导向的处置规程研究,将事故的判断和处置分多个步骤进行,提高了诊断的适应性,也降低了误判的风险,但在系统异常自动检测方面还存在的短板和缺陷,并且核动力系统异擦汗那个和故障种类繁多,通过枚举的方式很难覆盖所有故障类型,容易漏检。1) Mechanism analysis anomaly detection method. The mechanism analysis method developed earlier and is more popular. It is the traditional mainstream method of abnormal detection and fault diagnosis of nuclear power system. It is divided into simulation method and expert rule set-oriented method. The simulation method refers to the construction of a full-scale simulation model system of numerical calculation that is basically the same as that of the nuclear power system based on the structural characteristics of the nuclear power plant, using conservation equations and empirical formulas in physics and thermal engineering, and performing abnormal and fault analysis on the model system. and accident demonstrations to obtain rules and experience to guide the abnormal detection and fault diagnosis of the actual system. Its application has been extended to nuclear power system design verification and nuclear power operator training, and has become the mainstream engineering design verification in the nuclear power field. However, when applied to anomaly detection and fault diagnosis, this method has the problem of being difficult to synchronize with the actual system; the expert rule set method refers to the comprehensive use of physics, thermal equations and expert experience to set accident identification procedures, formulate Accident handling procedures are used to detect abnormality of the servo system, diagnose system failures, and guide the technical methods of accident handling. It is an abnormality detection and fault diagnosis method directly acting on nuclear power plants. It is also a commonly used nuclear power plant at home and abroad. However, its anomaly detection and accident identification procedures are complex, require a large number of human intervention, it is difficult to ensure its completeness, and it is easy to induce human errors, and its diagnosis and disposal process is either fixed and rigid and cannot be changed, or it is ambiguous and difficult to understand. Therefore, in recent years, many domestic and foreign experts have put forward the research on symptom-based or state-oriented disposal procedures, which divides the judgment and disposal of accidents into multiple steps, which improves the adaptability of diagnosis and reduces the risk of misjudgment. There are still shortcomings and defects in the automatic detection of system anomalies, and the nuclear power system has a wide variety of faults and faults. It is difficult to cover all fault types through enumeration, and it is easy to miss detection.
2)数据驱动异常检测法,近年来随着信息技术的不断进步,核动力系统信息化程度不断提高,大量核动力运行数据被完整详细的记录下来,数据驱动的异常检测方法也逐渐兴起,受到青睐,并体现了其独特的优势,比如在线能力强,贴近实装,准确性高等。2) Data-driven anomaly detection method. In recent years, with the continuous progress of information technology, the degree of informatization of nuclear power system has been continuously improved. A large number of nuclear power operation data have been recorded in complete and detailed, and data-driven anomaly detection methods have gradually emerged. It is favored and reflects its unique advantages, such as strong online ability, close to actual installation, and high accuracy.
数据驱动的异常检测和故障诊断方法主要分为两大类:一类是基于核动力装置全范围仿真模型系统计算数据的故障诊断方法,利用核动力装置全范围仿真模型系统可以计算不同工况,不同程度的系统典型故障和事故数据,构建核动力装置异常、事故数据谱库,然后利用机器学习方法训练事故辨识诊断模型,用于异常检测和故障诊断,主要通过故障特征来甄别系统异常,可以称为“黑名单法”,该方法利用数据库技术和人工智能技术有效地结合,解决了机理分析法中模拟仿真法在线支持困难的问题,然而,该方法过分依赖于全范围仿真模型的准确性,而仿真模型与实际核动力装置之间存在误差,这个误差在系统运行过程中并未被定量标定,也未被确定置信区间,误差在异常和事故工况下很可能会被扩大,甚至不收敛,因此,该方法存在不确定因素较多,还需进一步研究,该方向吸引了大量国内学者研究。The data-driven anomaly detection and fault diagnosis methods are mainly divided into two categories: one is the fault diagnosis method based on the calculation data of the nuclear power plant full-scale simulation model system. The nuclear power plant full-scale simulation model system can be used to calculate different working conditions. Different levels of typical system faults and accident data, build a nuclear power plant abnormality and accident data spectrum database, and then use machine learning methods to train an accident identification and diagnosis model for abnormality detection and fault diagnosis. Known as the "blacklist method", this method utilizes the effective combination of database technology and artificial intelligence technology to solve the difficult problem of online support for the simulation method in the mechanism analysis method. However, this method relies too much on the accuracy of the full-scale simulation model. , and there is an error between the simulation model and the actual nuclear power plant. This error has not been quantitatively calibrated during the operation of the system, and the confidence interval has not been determined. The error is likely to be enlarged under abnormal and accident conditions, or even not Therefore, there are many uncertain factors in this method, and further research is needed. This direction has attracted a large number of domestic scholars to study.
第二类数据驱动方法是直接挖掘核动力系统实际运行数据,获得规律来进行异常检测和诊断,本案按该思路研究系统异常自动检测技术,并提出一种异常检测算法。The second type of data-driven method is to directly mine the actual operating data of the nuclear power system and obtain the rules for abnormal detection and diagnosis. This case studies the automatic detection technology of system abnormality according to this idea, and proposes an abnormality detection algorithm.
基于核动力系统实际运行数据分析,研究系统异常检测技术,通常的研究思路是梳理系统发生过的异常、事件、事故数据,进行特征提取,然后依据故障特征构建故障辨识算法或模型,用于在线故障诊断,然而,核动系统是一种高可靠性的系统,其发生异常、事件和事故的概率极低,历史运行数据中异常工况数据极少,类型分布零散,出现随机,而事故数据更是几乎没有,样本少、分布广和特征不明显使得基于故障特征的系统异常自动检测技术研究思路陷入瓶颈。Based on the analysis of the actual operation data of the nuclear power system, the system abnormality detection technology is studied. The usual research idea is to sort out the abnormality, event, and accident data that have occurred in the system, perform feature extraction, and then construct a fault identification algorithm or model according to the fault characteristics. Fault diagnosis, however, the nuclear power system is a highly reliable system, and its probability of occurrence of anomalies, events and accidents is extremely low, there are very few abnormal conditions in historical operation data, and the type distribution is scattered and random, while accident data There are almost no samples. The small number of samples, wide distribution and inconspicuous features make the research idea of automatic detection technology of system anomalies based on fault characteristics into a bottleneck.
既然历史运行数据大部分为常用工况的正常运行数据,因此开拓思路想到借鉴网络安全领域的“白名单”理念,尝试利用系统正常运行数据构建系统的安全运行域以检测系统异常,实际是构造了一个由运行状态参数组成的多维向量数据空间,当表征系统运行状态的参数模型超出该数据空间时就认为系统异常,反之认为系统运行正常,该思路构建的异常检测算法通常具有较低的漏检率。Since most of the historical operation data is the normal operation data of common working conditions, the idea is to learn from the "whitelist" concept in the field of network security, and try to use the normal operation data of the system to construct the safe operation domain of the system to detect system anomalies. A multi-dimensional vector data space composed of operating state parameters is established. When the parameter model representing the operating state of the system exceeds the data space, the system is considered abnormal. Otherwise, the system is considered to be operating normally. detection rate.
利用正常运行数据研究构建安全运行域方法来检测核动力系统异常,思路较新颖,且需要持续海量的核动力系统实际运行数据作为支撑,而往往在核电领域,研究机构与核动力运行单位相互独立,因此,目前在核动力行业内公开发表相关著作较少。Using normal operation data to study and construct a safe operation domain method to detect nuclear power system anomalies is a novel idea, and requires continuous and massive actual operation data of the nuclear power system as support. In the field of nuclear power, research institutions and nuclear power operation units are often independent of each other. , therefore, few relevant works have been published in the nuclear power industry.
通过上述分析,现有技术存在的问题及缺陷为:现有船用核动力系统在线检测技术时,故障样本稀少且不完备。Through the above analysis, the existing problems and defects of the existing technology are: when the existing marine nuclear power system online detection technology is used, the fault samples are sparse and incomplete.
解决以上问题及缺陷的难度为:The difficulty of solving the above problems and defects is as follows:
与陆上核设施不同,船用核动力系统具有任务剖面复杂、运行工况多变、舰船核动力系统空间狭小,安全设施相对简单等特点,事故的处理能力相对薄弱,其异常和故障的处置对操纵员的经验判断与临场处置依赖程度也更大,这些因素增大了船用核动力系统的运行安全风险,同时也给运行操纵人员增添了额外的精神压力和工作负担。Different from onshore nuclear facilities, marine nuclear power system has the characteristics of complex mission profile, changeable operating conditions, small space for ship nuclear power system, relatively simple safety facilities, etc., and its ability to handle accidents is relatively weak. There is also a greater reliance on the operator's experience judgment and on-site disposal. These factors increase the operational safety risk of the marine nuclear power system, and also add additional mental pressure and workload to the operator.
解决以上问题及缺陷的意义为:The significance of solving the above problems and defects is:
需要安全、智能、人性化的运行支持技术手段,分担操纵人员压力,减少因疲劳引起的人因失误。系统异常检测技术是核动力运行支持技术体系中的一项关键技术,如其技术成熟智能,可以更有效的辅助操纵员值班,降低操纵员工作压力,减少人因失误,规避安全风险。Safe, intelligent and humanized operation support technical means are needed to share the pressure of operators and reduce human errors caused by fatigue. System anomaly detection technology is a key technology in the nuclear power operation support technology system. If its technology is mature and intelligent, it can more effectively assist operators on duty, reduce operator work pressure, reduce human errors, and avoid safety risks.
发明内容SUMMARY OF THE INVENTION
针对现有技术存在的问题,本发明提供了一种基于逻辑距离表征安全运行域的船用核动力异常检测方法。In view of the problems existing in the prior art, the present invention provides a marine nuclear power anomaly detection method based on the logical distance characterizing the safe operation domain.
本发明是这样实现的,一种基于逻辑距离表征安全运行域的船用核动力异常检测方法,所述基于逻辑距离表征安全运行域的船用核动力异常检测方法包括:The present invention is implemented in the following way: a method for detecting anomalies in marine nuclear power based on a logical distance representing a safe operating domain, the method for detecting anomalies in marine nuclear power based on a logical distance representing a safe operating domain includes:
步骤一,通过分析核动力系统历史运行数据,检索系统典型运行工况数据样本;Step 1: Retrieve data samples of typical operating conditions of the system by analyzing the historical operating data of the nuclear power system;
步骤二,构建标准运行工况样本库,遴选表征运行状态的特征参数;Step 2: Build a sample library of standard operating conditions, and select characteristic parameters that characterize the operating state;
步骤三,在标准运行工况样本与待检测的系统运行状态数据之间,构建逻辑距离计算函数;Step 3: Build a logical distance calculation function between the standard operating condition sample and the system operating state data to be detected;
步骤四,通过历史运行数据进行训练学习、仿真计算或收集的异常数据样本的检验验证确定判定系统异常的逻辑距离阈值即安全域归属阈值,检测系统异常或检索系统历史运行数据中的异常情况。Step 4: Determine the logical distance threshold for judging system anomalies, that is, the security domain attribution threshold, to detect system anomalies or retrieve anomalies in the system historical operation data.
进一步,步骤二中,所述标准运行工况体系包括:稳态运行类工况、切换类工况以及启/停堆类工况;Further, in
所述每一个运行工况即为一个运行模态,对应一个运行域和一个标准运行工况;Each operating condition is an operating mode, corresponding to an operating domain and a standard operating condition;
所述稳态运行类工况包括行进一、行进二、行进三…行进N工况;The steady-state operating conditions include
所述切换类工况包括行进M切换至行进L,M≠L且M,L∈(1,2,┅,N;The switching conditions include switching from traveling M to traveling L, where M≠L and M, L∈(1, 2, ┅, N;
所述启/停堆类工况包括冷启动即物理启动、热启动、正常停堆、紧急停堆工况。The start/stop conditions include cold start, namely physical start, hot start, normal shutdown, and emergency shutdown.
进一步,步骤二中,所述标准运行工况体系构建方法还包括:采取模拟试验与数据分析两种策略相结合的方式遴选单模态标准工况;Further, in
所述模拟试验包括利用核动力系统配套的全范围模拟器依据典型工况运行参数系统运行状态作为标准状态;The simulation test includes the use of a full-scale simulator provided with the nuclear power system to take the operating state of the system as a standard state according to the operating parameters of typical operating conditions;
所述数据分析包括基于系统历史运行数据,计算获取样本中某一运行工况下所有运行数据样本之间的逻辑距离,并找出其中与其他所有运行数据样本距离之和最小的一个运行数据样本实例,作为表征该运行工况即运行模态的标准运行工况样本,或计算运行区域内数据样本集的逻辑中心作为标准运行工况。The data analysis includes calculating the logical distance between all the operating data samples under a certain operating condition in the acquired samples based on the historical operating data of the system, and finding out the one operating data sample with the smallest sum of distances from all other operating data samples. For example, as a standard operating condition sample representing the operating condition, that is, the operating mode, or calculating the logical center of the data sample set in the operating area as the standard operating condition.
进一步,步骤二中,所述表征运行状态的特征参数包括:堆右入口温度、堆左入口温度、堆右出口温度、堆左出口温度、右回路流量、左回路流量、稳压器压力、稳压器温度、稳压器水位、净化离子交换器进口温度、二回路负荷功率、核功率、净化水温度、净化水流量、1#主蒸汽压力、2#主蒸汽压力、1#蒸汽发生器水位、2#蒸汽发生器水位、1#蒸汽发生器给水流量、2#蒸汽发生器给水流量、1#蒸汽发生器蒸汽流量、2#蒸汽发生器蒸汽流量以及艉轴转速共计23个表征参数。Further, in
进一步,步骤三中,所述逻辑距离计算函数构建方法包括:Further, in
所述逻辑距离计算函数用于定量分析系统不同稳定运行状态之间逻辑关系远近的度量函数;包括稳态功率运行工况的逻辑距离计算函数以及启、停堆和工况切换瞬态过程的逻辑距离计算函数;The logical distance calculation function is used to quantitatively analyze the metric function of the logical relationship between different stable operating states of the system; including the logical distance calculation function of the steady-state power operating condition and the logic of the startup, shutdown and operating conditions switching transient process distance calculation function;
(1)稳态功率运行工况的逻辑距离计算函数:采用加权的欧氏距离算法用于计算系统稳态间的逻辑距离,具体如下:(1) Calculation function of logical distance for steady-state power operating conditions: The weighted Euclidean distance algorithm is used to calculate the logical distance between the steady states of the system, as follows:
稳态功率运行时,计算逻辑距离的两个多维向量分别为:During steady-state power operation, the two multidimensional vectors for calculating the logical distance are:
A=[a1,a2,a3,a4,…,a23]B=[b1,b2,b3,…,b23]A=[a1,a2,a3,a4,...,a23]B=[b1,b2,b3,...,b23]
对应的权值向量为:The corresponding weight vector is:
C=[c1,c2,c3,c4,…,c23]C=[c1,c2,c3,c4,...,c23]
c1+c2+,…,+c23=100c1+c2+,…,+c23=100
定义的加权欧氏距离计算函数为:The defined weighted Euclidean distance calculation function is:
(2)启、停堆和工况切换瞬态过程的逻辑距离计算函数:(2) Calculation function of logical distance for the transient process of startup, shutdown and switching of operating conditions:
用于计算运行状态逻辑距离的特征向量表示为:The eigenvectors used to calculate the operating state logical distance are expressed as:
A=[a1,a2,a3,a4,…,a23,a1′,a2′,a3′,…,ak′],k<23A=[a1, a2, a3, a4, ..., a23, a1', a2', a3', ..., ak'], k<23
B=[b1,b2,b3,b4,…,b23,b1′,b2′,b3′,…,bk′],k<23B=[b1, b2, b3, b4,..., b23, b1', b2', b3',..., bk'], k<23
对应权值向量为:The corresponding weight vector is:
C=[c1,c2,c3,c4,…,c23,c1′,c2′,c3′,…,ck′],k<23C=[c1, c2, c3, c4, ..., c23, c1', c2', c3', ..., ck'], k<23
c1+c2+…,c23=100c1+c2+…, c23=100
c1′+c2′+c3′,…,+ck′=100c1'+c2'+c3',...,+ck'=100
对于瞬态运行过程逻辑距离,定义为特征参数逻辑距离和特征参数变化速率即差分逻辑距离两个计算函数,具体如下所示:For the logical distance of the transient operation process, it is defined as two calculation functions of the characteristic parameter logical distance and the characteristic parameter change rate, that is, the differential logical distance, as shown below:
本发明的另一目的在于提供一种实施所述基于逻辑距离表征安全运行域的船用核动力异常检测方法的基于逻辑距离表征安全运行域的船用核动力异常检测系统,所述基于逻辑距离表征安全运行域的船用核动力异常检测系统包括:Another object of the present invention is to provide a marine nuclear power anomaly detection system based on a logical distance to characterize a safe operation domain by implementing the method for detecting a marine nuclear power anomaly based on a logical distance to characterize a safe operating domain. Marine nuclear power anomaly detection systems in the operational domain include:
构建模块,包括运行状态表征向量构建单元、核动力系统标准运行工况体系构建单元、逻辑距离函数构建单元、阈值确定单元、数值试验校验验证单元;用于进行运行状态表征向量构建、核动力系统标准运行工况体系构建、逻辑距离函数构建单元、阈值确定以及数值试验校验验证;Building modules, including operating state characterization vector building unit, nuclear power system standard operating condition system building unit, logical distance function building unit, threshold determination unit, numerical test verification unit; used for operating state characterization vector building, nuclear power system System construction of standard operating conditions, logical distance function construction unit, threshold determination and numerical test verification;
测试与试验模块,用于在线检测系统异常或检索系统历史运行数据中的异常情况。The test and experiment module is used for online detection of system anomalies or retrieval of anomalies in the historical operating data of the system.
进一步,所述构建模块包括:Further, the building blocks include:
运行状态表征向量构建单元,用于表征系统运行状态的特征参数,并构建系统运行状态特征向量;The operating state characterization vector construction unit is used to characterize the characteristic parameters of the operating state of the system, and construct the feature vector of the operating state of the system;
核动力系统标准运行工况体系构建单元,用于构建系统标准运行工况体系,并确定单模态标准运行工况;The nuclear power system standard operating condition system construction unit is used to construct the system standard operating condition system and determine the single-modal standard operating condition;
逻辑距离函数构建单元,用于分别构建稳态功率运行工况的逻辑距离计算函数以及瞬态运行过程的逻辑距离计算函数;The logical distance function construction unit is used to respectively construct the logical distance calculation function of the steady-state power operating condition and the logical distance calculation function of the transient operation process;
阈值确定单元,用于通过系统正常运行数据样本训练确定安全域归属阈值,并通过仿真计算或者收集的异常数据样本的检验验证安全域归属阈值;The threshold determination unit is used to determine the security domain attribution threshold through the training of data samples of normal operation of the system, and verify the security domain attribution threshold through simulation calculation or inspection of the collected abnormal data samples;
数值试验校验验证单元,用于进行稳态工况异常检测数值试验验证以及变换工况状态异常检测数值试验验证。Numerical test verification and verification unit is used to perform numerical test verification of abnormality detection in steady state conditions and numerical test verification of abnormality detection in changing conditions.
本发明另一目的在于提供一种计算机设备,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行如下步骤:Another object of the present invention is to provide a computer device, the computer device includes a memory and a processor, the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the following steps :
通过分析核动力系统历史运行数据,检索系统典型运行工况数据样本;By analyzing the historical operating data of the nuclear power system, the data samples of the typical operating conditions of the system are retrieved;
构建标准运行工况样本库,遴选表征运行状态的特征参数;Build a sample library of standard operating conditions, and select characteristic parameters that characterize the operating state;
在标准运行工况样本与待检测的系统运行状态数据之间,构建逻辑距离计算函数;Build a logical distance calculation function between the standard operating condition sample and the system operating state data to be detected;
通过历史运行数据进行训练学习、仿真计算或收集的异常数据样本的检验验证确定判定系统异常的逻辑距离阈值即安全域归属阈值,检测系统异常或检索系统历史运行数据中的异常情况。The logical distance threshold for judging system anomalies, that is, the security domain attribution threshold, is determined through training learning, simulation calculation, or the inspection and verification of abnormal data samples collected through historical operation data, to detect system anomalies or retrieve anomalies in system historical operation data.
本发明另一目的在于提供一种计算机可读存储介质,存储有计算机程序,所述计算机程序被处理器执行时,使得所述处理器执行如下步骤:Another object of the present invention is to provide a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, causes the processor to perform the following steps:
通过分析核动力系统历史运行数据,检索系统典型运行工况数据样本;By analyzing the historical operating data of the nuclear power system, the data samples of the typical operating conditions of the system are retrieved;
构建标准运行工况样本库,遴选表征运行状态的特征参数;Build a sample library of standard operating conditions, and select characteristic parameters that characterize the operating state;
在标准运行工况样本与待检测的系统运行状态数据之间,构建逻辑距离计算函数;Build a logical distance calculation function between the standard operating condition sample and the system operating state data to be detected;
通过历史运行数据进行训练学习、仿真计算或收集的异常数据样本的检验验证确定判定系统异常的逻辑距离阈值即安全域归属阈值,检测系统异常或检索系统历史运行数据中的异常情况Through historical operation data for training and learning, simulation calculation or the inspection and verification of collected abnormal data samples, the logical distance threshold for judging system abnormality, that is, the security domain attribution threshold, is used to detect system abnormality or retrieve abnormal conditions in system historical operation data.
结合上述的所有技术方案,本发明所具备的优点及积极效果为:Combined with all the above-mentioned technical solutions, the advantages and positive effects possessed by the present invention are:
本发明从核动力系统历史运行数据出发,基于网络安全领域的“白名单”理念,提出并构造了一种数据驱动的核动力系统安全运行域表征算法,用于辅助操纵员在线检测系统异常。Starting from the historical operation data of the nuclear power system and based on the "white list" concept in the network security field, the invention proposes and constructs a data-driven nuclear power system safe operation domain characterization algorithm for assisting the operator to detect system abnormalities online.
本发明基于数据驱动的理念,从系统正常运行数据出发,结合船用堆特点,构建了一种基于系统安全运行域的异常检测方法,并基于逻辑距离计算设计了系统安全运行域表征函数,用于在线监测船用核动力系统异常,最后,以系统常用运行工况为对象,对算法进行了数值试验验证。结果表明设计算法能有效检测系统异常和故障,具有良好可靠性和可解释性,并且其检测力度可调节。Based on the data-driven concept, the invention starts from the normal operation data of the system and combines the characteristics of the marine reactor to construct an anomaly detection method based on the system safe operation domain, and designs the system safe operation domain characterization function based on the logical distance calculation, which is used for The abnormality of the marine nuclear power system is monitored online. Finally, the algorithm is verified by numerical experiments with the common operating conditions of the system as the object. The results show that the designed algorithm can effectively detect system anomalies and faults, with good reliability and interpretability, and its detection strength can be adjusted.
本发明直接从对象核动力装置运行状态数据出发构造异常检测模型,贴近实际系统,是为系统个性化定制的异常检测方法,从机理角度比其他算法具有更强的针对性和适用性。The invention constructs an anomaly detection model directly from the operating state data of the target nuclear power plant, which is close to the actual system, is an anomaly detection method customized for the system, and has stronger pertinence and applicability than other algorithms from the perspective of mechanism.
本发明从正常运行数据出发构建系统安全运行域检测系统异常,有效解决了核动力系统可靠性高,异常和故障数据样本少,无法提取故障特征的瓶颈问题。The invention starts from normal operation data to construct a system safe operation domain to detect system abnormality, and effectively solves the bottleneck problem of high reliability of nuclear power system, few abnormality and fault data samples, and inability to extract fault features.
本发明通过专家经验和数值分析遴选出的表征系统运行状态的23个参数,能够良好的体现不同运行工况下系统特征;通过数值验证划分。The present invention selects 23 parameters to characterize the operating state of the system through expert experience and numerical analysis, which can well reflect the system characteristics under different operating conditions, and is divided by numerical verification.
本发明基于逻辑距离的安全运行域表征函数,通过显性的公式表征不同运行状态间的细微差异,与神经网络类数据驱动算法的异常检测技术相比,可解释性强,并且通过阈值调节,可线性调节算法检测力度,提高或降低对异常和故障检测的灵敏度。The present invention is based on the logical distance-based safe operating domain characterization function, and expresses the subtle differences between different operating states through explicit formulas. The detection strength of the algorithm can be adjusted linearly to increase or decrease the sensitivity to anomaly and fault detection.
本发明通过数值试验验证,设计异常检测算法,可以检测出系统的微小异常,漏检率极低,检测速度快,比如即使对征兆不明显且显现缓慢的一回路系统4mm尺寸的微小破口,在1分钟以内就可以发现明显故障征兆,比其他算法敏感性更强。The invention is verified by numerical experiments, and an abnormality detection algorithm is designed, which can detect the small abnormality of the system, the missed detection rate is extremely low, and the detection speed is fast. Obvious signs of failure can be found within 1 minute, which is more sensitive than other algorithms.
本发明检测算法针对不同异常和故障,检测表征结果呈现出不同特征,为进一步辨识异常和故障类型提供了信息空间,优于神经网络类数据驱动异常检测算法The detection algorithm of the present invention presents different characteristics in the detection and representation results for different anomalies and faults, provides an information space for further identifying the types of anomalies and faults, and is superior to the neural network data-driven anomaly detection algorithm
本发明算法可推广应用到核电厂机组的运行状态异常检测中。The algorithm of the invention can be extended and applied to the abnormal detection of the operating state of the nuclear power plant unit.
对比的技术效果或者实验效果。本发明图10、11体现了不同类型故障,检测结果呈现出不同特征。The technical effect or experimental effect of the comparison. Figures 10 and 11 of the present invention show different types of faults, and the detection results show different characteristics.
本发明遴选23个状态表征技术参数的数值分析依据。The invention selects the numerical analysis basis of 23 state characterization technical parameters.
本发明对比图10、11、9与图12(基于神经网络类的故障诊断结果)可看出神经网络类的异常检测及故障诊断算法仅能输出一个匹配最相似的分类号,无法体现故障的渐变过程,可解释性差;而本发明可以体现故障程度和不同故障的相近性,以及故障渐变过程,或者偶发时机,可解释性强。Comparing Figures 10, 11, 9 and Figure 12 (fault diagnosis results based on neural network) in the present invention, it can be seen that the abnormal detection and fault diagnosis algorithm of neural network can only output a classification number that matches the most similar, and cannot reflect the fault. The gradual change process has poor interpretability; while the present invention can reflect the degree of failure and the similarity of different faults, as well as the gradual change process of the fault, or the occasional time, and the interpretability is strong.
与基于故障特征分析的传统异常检测方法相比较,本案检测范围广,且不受故障样本稀少的限制。如图12,传统异常检测只能针对已知的经典的单一类型故障实施检测,或者针对单一参数进行系统异常检测,且通常还需要使用仿真模拟器计算获得或者补充故障训练样本,而模拟系统与实际系统存在固有计算误差,会影响诊断方法的有效性;对照图7、8、9,本案基于实际系统构建安全运行域,贴近实际,且可同时检测各种类型的异常,可检测未知的异常和故障,拓展了检测范围,降低了漏检率,同时可根据需求通过阈值调节检测力度,还有效规避了故障样本稀少且分布不均的现实问题。Compared with traditional anomaly detection methods based on fault feature analysis, this case has a wide detection range and is not limited by the scarcity of fault samples. As shown in Figure 12, traditional anomaly detection can only detect a known classic single type of fault, or perform system anomaly detection for a single parameter, and usually needs to use a simulation simulator to calculate or supplement the fault training samples, while the simulation system is similar to There are inherent calculation errors in the actual system, which will affect the effectiveness of the diagnosis method; compared with Figures 7, 8, and 9, this case builds a safe operation domain based on the actual system, which is close to reality, and can detect various types of anomalies at the same time, and can detect unknown anomalies. The detection range is expanded and the missed detection rate is reduced. At the same time, the detection intensity can be adjusted through the threshold value according to the requirements, and the practical problem of rare and uneven distribution of fault samples can be effectively avoided.
附图说明Description of drawings
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例中所需要使用的附图做简单的介绍,显而易见地,下面所描述的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下还可以根据这些附图获得其他的附图。In order to explain the technical solutions of the embodiments of the present application more clearly, the following will briefly introduce the drawings that need to be used in the embodiments of the present application. Obviously, the drawings described below are only some embodiments of the present application. For those of ordinary skill in the art, other drawings can also be obtained from these drawings without creative effort.
图1是本发明实施例提供的基于逻辑距离表征安全运行域的船用核动力异常检测方法流程图。FIG. 1 is a flowchart of a method for detecting anomalies in marine nuclear power based on a logical distance representing a safe operating domain provided by an embodiment of the present invention.
图2是本发明实施例提供的基于逻辑距离表征安全运行域的船用核动力异常检测系统结构示意图;2 is a schematic structural diagram of a marine nuclear power anomaly detection system based on a logical distance to characterize a safe operating domain provided by an embodiment of the present invention;
图中:1、构建模块;2、测试与试验模块;11、运行状态表征向量构建单元;12、核动力系统标准运行工况体系构建单元;13、逻辑距离函数构建单元;14、阈值确定单元;15、数值试验校验验证单元。In the figure: 1. Building module; 2. Test and experiment module; 11. Operating state representation vector building unit; 12. Nuclear power system standard operating condition system building unit; 13. Logical distance function building unit; 14. Threshold determination unit ; 15. Numerical test verification unit.
图3是本发明实施例提供的基于逻辑距离表征安全运行域的船用核动力异常检测系统架构图。FIG. 3 is an architecture diagram of a marine nuclear power anomaly detection system based on a logical distance representation of a safe operating domain provided by an embodiment of the present invention.
图4是本发明实施例提供的安全运行域与运行状态数据逻辑关系示意图。FIG. 4 is a schematic diagram of a logical relationship between a safe operating domain and operating state data provided by an embodiment of the present invention.
图5是本发明实施例提供的行进一样本数据部分关键参数分布示意图。FIG. 5 is a schematic diagram of the distribution of key parameters of a traveling sample data part according to an embodiment of the present invention.
图6是本发明实施例提供的行进一偶发性异常检测验证示意图。FIG. 6 is a schematic diagram of the detection and verification of an occasional anomaly in progress provided by an embodiment of the present invention.
图7是本发明实施例提供的行进一破坏性损伤故障检测验证示意图。FIG. 7 is a schematic diagram of a traveling-destructive damage fault detection verification provided by an embodiment of the present invention.
图8是本发明实施例提供的快速升负荷工况状态参数逻辑距离分布示意图。FIG. 8 is a schematic diagram of the distribution of logical distances of state parameters in a fast-rising load condition provided by an embodiment of the present invention.
图9是本发明实施例提供的快速升负荷工况系统异常检测验证示意图。FIG. 9 is a schematic diagram of anomaly detection and verification of a system under a fast-rising load condition provided by an embodiment of the present invention.
图10是本发明实施例提供的不同类型故障,检测结果呈现出不同特征图一。FIG. 10 shows different types of faults provided by the embodiment of the present invention, and the detection results show different characteristic diagram 1.
图11是本发明实施例提供的不同类型故障,检测结果呈现出不同特征图二。FIG. 11 shows different types of faults provided by an embodiment of the present invention, and the detection results show different characteristics, as shown in FIG. 2 .
图12是本发明实施例提供的传统采用神经网络模型对单一参数的跟踪与异常检测图。FIG. 12 is a diagram of tracking and anomaly detection of a single parameter using a traditional neural network model provided by an embodiment of the present invention.
具体实施方式Detailed ways
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, but not to limit the present invention.
针对现有技术存在的问题,本发明提供了一种基于逻辑距离表征安全运行域的船用核动力异常检测方法,下面结合附图对本发明作详细的描述。In view of the problems existing in the prior art, the present invention provides a method for detecting anomalies in marine nuclear power based on the logical distance characterizing the safe operation domain. The present invention is described in detail below with reference to the accompanying drawings.
如图1所示,本发明实施例提供的基于逻辑距离表征安全运行域的船用核动力异常检测方法包括:As shown in FIG. 1 , the method for detecting anomalies in marine nuclear power based on a logical distance characterizing a safe operating domain provided by an embodiment of the present invention includes:
S101,通过分析核动力系统历史运行数据,检索系统典型运行工况数据样本;S101, by analyzing the historical operating data of the nuclear power system, retrieve data samples of typical operating conditions of the system;
S102,构建标准运行工况样本库,遴选表征运行状态的特征参数;S102, build a standard operating condition sample library, and select characteristic parameters that characterize the operating state;
S103,在标准运行工况样本与待检测的系统运行状态数据之间,构建逻辑距离计算函数;S103, constructing a logical distance calculation function between the standard operating condition sample and the system operating state data to be detected;
S104,通过历史运行数据进行训练学习、仿真计算或收集的异常数据样本的检验验证确定判定系统异常的逻辑距离阈值即安全域归属阈值,检测系统异常或检索系统历史运行数据中的异常情况。S104: Determine the logical distance threshold for judging system abnormality, namely the security domain attribution threshold, by performing training learning, simulation calculation, or testing and verification of abnormal data samples collected through historical operating data, and detecting system abnormality or retrieving abnormal conditions in the historical operating data of the system.
步骤S102中,本发明实施例提供的标准运行工况体系包括:稳态运行类工况、切换类工况以及启/停堆类工况;In step S102, the standard operating condition system provided by the embodiment of the present invention includes: steady-state operating conditions, switching operating conditions, and start/stop operating conditions;
所述每一个运行工况即为一个运行模态,对应一个运行域和一个标准运行工况;Each operating condition is an operating mode, corresponding to an operating domain and a standard operating condition;
所述稳态运行类工况包括行进一、行进二、行进三…行进N工况;The steady-state operating conditions include
所述切换类工况包括行进M切换至行进L,M≠L且M,L∈(1,2,┅,N;The switching conditions include switching from traveling M to traveling L, where M≠L and M, L∈(1, 2, ┅, N;
所述启/停堆类工况包括冷启动即物理启动、热启动、正常停堆、紧急停堆工况。The start/stop conditions include cold start, physical start, hot start, normal shutdown, and emergency shutdown.
步骤S102中,本发明实施例提供的标准运行工况体系构建方法还包括:采取模拟试验与数据分析两种策略相结合的方式遴选单模态标准工况;In step S102, the method for constructing a standard operating condition system provided by the embodiment of the present invention further includes: selecting a single-modal standard operating condition by combining two strategies of simulation test and data analysis;
所述模拟试验包括利用核动力系统配套的全范围模拟器依据典型工况运行参数系统运行状态作为标准状态;The simulation test includes using a full-scale simulator matched with the nuclear power system to take the operating state of the system as a standard state according to the operating parameters of typical operating conditions;
所述数据分析包括基于系统历史运行数据,计算获取样本中某一运行工况下所有运行数据样本之间的逻辑距离,并找出其中与其他所有运行数据样本距离之和最小的一个运行数据样本实例,作为表征该运行工况即运行模态的标准运行工况样本,或计算运行区域内数据样本集的逻辑中心作为标准运行工况。The data analysis includes calculating the logical distance between all the operating data samples under a certain operating condition in the acquired samples based on the historical operating data of the system, and finding out the one operating data sample with the smallest sum of distances from all other operating data samples. For example, as a standard operating condition sample representing the operating condition, that is, the operating mode, or calculating the logical center of the data sample set in the operating area as the standard operating condition.
步骤S102中,本发明实施例提供的表征运行状态的特征参数包括:堆右入口温度、堆左入口温度、堆右出口温度、堆左出口温度、右回路流量、左回路流量、稳压器压力、稳压器温度、稳压器水位、净化离子交换器进口温度、二回路负荷功率、核功率、净化水温度、净化水流量、1#主蒸汽压力、2#主蒸汽压力、1#蒸汽发生器水位、2#蒸汽发生器水位、1#蒸汽发生器给水流量、2#蒸汽发生器给水流量、1#蒸汽发生器蒸汽流量、2#蒸汽发生器蒸汽流量以及艉轴转速共计23个表征参数。In step S102, the characteristic parameters that characterize the operating state provided by the embodiment of the present invention include: stack right inlet temperature, stack left inlet temperature, stack right outlet temperature, stack left outlet temperature, right loop flow, left loop flow, and pressure of the regulator. , pressure regulator temperature, pressure regulator water level, purification ion exchanger inlet temperature, secondary circuit load power, nuclear power, purified water temperature, purified water flow, 1# main steam pressure, 2# main steam pressure, 1# steam generation There are a total of 23 characteristic parameters including the water level of the steam generator, the water level of the 2# steam generator, the feed water flow of the 1# steam generator, the feed water flow of the 2# steam generator, the steam flow of the 1# steam generator, the steam flow of the 2# steam generator and the speed of the stern shaft. .
步骤S103中,本发明实施例提供的逻辑距离计算函数构建方法包括:In step S103, the method for constructing a logical distance calculation function provided by the embodiment of the present invention includes:
所述逻辑距离计算函数用于定量分析系统不同稳定运行状态之间逻辑关系远近的度量函数;包括稳态功率运行工况的逻辑距离计算函数以及启、停堆和工况切换瞬态过程的逻辑距离计算函数;The logical distance calculation function is used to quantitatively analyze the metric function of the logical relationship between different stable operating states of the system; including the logical distance calculation function of the steady-state power operating condition and the logic of the startup, shutdown and operating conditions switching transient process distance calculation function;
(1)稳态功率运行工况的逻辑距离计算函数:采用加权的欧氏距离算法用于计算系统稳态间的逻辑距离,具体如下:(1) Calculation function of logical distance for steady-state power operating conditions: The weighted Euclidean distance algorithm is used to calculate the logical distance between the steady states of the system, as follows:
稳态功率运行时,计算逻辑距离的两个多维向量分别为:During steady-state power operation, the two multidimensional vectors for calculating the logical distance are:
A=[a1,a2,a3,a4,…,a23]B=[b1,b2,b3,…,b23]A=[a1,a2,a3,a4,...,a23]B=[b1,b2,b3,...,b23]
对应的权值向量为:The corresponding weight vector is:
C=[c1,c2,c3,c4,…,c23]C=[c1,c2,c3,c4,...,c23]
c1+c2+,…,+c23=100c1+c2+,…,+c23=100
定义的加权欧氏距离计算函数为:The defined weighted Euclidean distance calculation function is:
(2)启、停堆和工况切换瞬态过程的逻辑距离计算函数:(2) Calculation function of logical distance for the transient process of startup, shutdown and switching of operating conditions:
用于计算运行状态逻辑距离的特征向量表示为:The eigenvectors used to calculate the operating state logical distance are expressed as:
A=[a1,a2,a3,a4,…,a23,a1′,a2′,a3′,…,ak′],k<23A=[a1, a2, a3, a4, ..., a23, a1', a2', a3', ..., ak'], k<23
B=[b1,b2,b3,b4,…,b23,b1′,b2′,b3′,…,bk′],k<23B=[b1, b2, b3, b4,..., b23, b1', b2', b3',..., bk'], k<23
对应权值向量为:The corresponding weight vector is:
C=[c1,c2,c3,c4,…,c23,c1′,c2′,c3′,…,ck′],k<23C=[c1, c2, c3, c4, ..., c23, c1', c2', c3', ..., ck'], k<23
c1+c2+…,c23=100c1+c2+…, c23=100
c1′+c2′+c3′,…,+ck′=100c1'+c2'+c3',...,+ck'=100
对于瞬态运行过程逻辑距离,定义为特征参数逻辑距离和特征参数变化速率即差分逻辑距离两个计算函数,具体如下所示:For the logical distance of the transient operation process, it is defined as two calculation functions of the characteristic parameter logical distance and the characteristic parameter change rate, that is, the differential logical distance, as shown below:
如图2-图3所示,本发明实施例提供的基于逻辑距离表征安全运行域的船用核动力异常检测系统包括:As shown in FIG. 2-FIG. 3, the marine nuclear power anomaly detection system based on the logical distance characterizing the safe operation domain provided by the embodiment of the present invention includes:
构建模块1,包括运行状态表征向量构建单元11、核动力系统标准运行工况体系构建单元12、逻辑距离函数构建单元13、阈值确定单元14、数值试验校验验证单元15;用于进行运行状态表征向量构建、核动力系统标准运行工况体系构建、逻辑距离函数构建单元、阈值确定以及数值试验校验验证;The
测试与试验模块2,用于在线检测系统异常或检索系统历史运行数据中的异常情况。Test and
本发明实施例提供的构建模块1包括:The
运行状态表征向量构建单元11,用于表征系统运行状态的特征参数,并构建系统运行状态特征向量;The operating state characterizing
核动力系统标准运行工况体系构建单元12,用于构建系统标准运行工况体系,并确定单模态标准运行工况;The nuclear power system standard operating condition
逻辑距离函数构建单元13,用于分别构建稳态功率运行工况的逻辑距离计算函数以及瞬态运行过程的逻辑距离计算函数;The logical distance
阈值确定单元14,用于通过系统正常运行数据样本训练确定安全域归属阈值,并通过仿真计算或者收集的异常数据样本的检验验证安全域归属阈值;
数值试验校验验证单元15,用于进行稳态工况异常检测数值试验验证以及变换工况状态异常检测数值试验验证。The numerical test verification and
下面结合具体实施例对本发明的技术方案作进一步说明。The technical solutions of the present invention will be further described below in conjunction with specific embodiments.
实施例:Example:
针对船用核动力系统在线检测技术时,故障样本稀少且不完备的问题,基于数据驱动的理念,从系统正常运行数据出发,结合船用堆特点,设计了一种基于系统安全运行域的异常检测算法,并基于逻辑距离计算设计了系统安全运行域表征函数,用于在线监测船用核动力系统异常,最后,以系统常用运行工况为对象,对算法进行了数值试验验证。结果表明设计算法能有效检测系统异常和故障,具有良好可靠性和可解释性,并且其检测力度可调节。Aiming at the problem of rare and incomplete fault samples in the online detection technology of marine nuclear power system, based on the data-driven concept, starting from the normal operation data of the system, and combining the characteristics of marine reactors, an anomaly detection algorithm based on the system's safe operation domain is designed. , and based on the logical distance calculation, the system safety operating domain characterization function is designed to monitor the abnormality of the marine nuclear power system on-line. The results show that the designed algorithm can effectively detect system anomalies and faults, with good reliability and interpretability, and its detection strength can be adjusted.
算法体系设计主要包括:运行状态表征向量设计、核动力系统标准运行工况体系设计、逻辑距离函数设计、阈值确定算法设计、数值试验校验验证五项主要内容。The algorithm system design mainly includes five main contents: operating state representation vector design, nuclear power system standard operating condition system design, logical distance function design, threshold determination algorithm design, and numerical test verification and verification.
基本步骤是通过分析核动力系统历史运行数据,检索系统典型运行工况数据样本,构建标准运行工况样本库,遴选表征运行状态的特征参数,在标准运行工况样本与待检测的系统运行状态数据之间,构建逻辑距离计算函数,在通过历史运行数据进行训练学习,最终确定判定系统异常的逻辑距离阈值,用于在线检测系统异常或检索系统历史运行数据中的异常情况。The basic steps are to analyze the historical operating data of the nuclear power system, retrieve data samples of typical operating conditions of the system, build a sample library of standard operating conditions, and select characteristic parameters that characterize the operating state. Between the data, a logical distance calculation function is constructed, and the historical operating data is used for training and learning, and finally the logical distance threshold for determining the abnormality of the system is determined, which is used to detect the abnormality of the system online or retrieve the abnormality in the historical operating data of the system.
1、表征核动力系统运行状态的特征向量的构建1. Construction of eigenvectors that characterize the operating state of the nuclear power system
核动力系统作为一个存在安全风险的能源系统,对其运维的可靠性要求都很高,因此其运行过程中会在线记录大量监测参数,将全部参数都作为表征系统运行状态的特征参数,固然可以保留最多的运行状态特征信息,但是一方面,监测参数数量庞大,使得设计检测算法计算量巨大,难以保证异常检测的实时性;另一方面,大量弱相关参数的重要度权值不容易确定,会影响检测结果的稳定性。因此需要对参数进行遴选,从而获得表征系统运行状态的强相关参数。考虑到可靠性冗余设计,参数间的耦合关系,参数与系统运行状态间的相关性强弱,数据可获取性和完备性,再经机理分析和数据试验筛选,遴选出如下23个参数作为表征系统运行状态的特征参数,构成系统运行状态特征向量。As an energy system with security risks, the nuclear power system has high reliability requirements for its operation and maintenance. Therefore, a large number of monitoring parameters are recorded online during its operation, and all parameters are used as characteristic parameters to characterize the operating state of the system. It can retain the most feature information of the running state, but on the one hand, the large number of monitoring parameters makes the design detection algorithm a huge amount of calculation, and it is difficult to ensure the real-time performance of anomaly detection; on the other hand, the importance weights of a large number of weakly correlated parameters are not easy to determine. , which will affect the stability of the test results. Therefore, it is necessary to select the parameters so as to obtain strongly correlated parameters that characterize the operating state of the system. Taking into account the reliability redundancy design, the coupling relationship between parameters, the correlation between parameters and the system operating state, the availability and completeness of data, and then through mechanism analysis and data test screening, the following 23 parameters were selected as The characteristic parameters that characterize the operating state of the system constitute the characteristic vector of the operating state of the system.
表1核动力系统运行状态表征参数Table 1 Characterization parameters of the operating state of the nuclear power system
2、系统标准运行工况的构建2. Construction of system standard operating conditions
1)标准运行工况体系构建1) Construction of standard operating condition system
船用核动力系统因船舶机动性需求,需航行在不同航速,因此需要设计多个稳态运行工况;同时还会时常在各工况间切换,增加了许多正常转换工况;核动力装置因其工作原理和设计特点,其启、停堆过程持续时间很长,是异常和故障易发的运行工况,且与功率运行差异较大,也需要单列为专门的运行工况,上述实际情况导致系统标准运行工况呈现多模态,因此需要构建一个系统标准运行工况体系,来支撑算法设计,船用核动力系统运行工况体系如下表所示,每一个运行工况即为一个运行模态,对应一个运行域和一个标准运行工况。Due to the mobility requirements of ships, the marine nuclear power system needs to sail at different speeds, so it is necessary to design multiple steady-state operating conditions; at the same time, it will often switch between various operating conditions, adding many normal conversion conditions; Its working principle and design characteristics, its start-up and shutdown process lasts a long time, it is an abnormal and fault-prone operating condition, and it is quite different from power operation, and it also needs to be listed as a special operating condition. The above actual situation As a result, the standard operating conditions of the system are multi-modal, so it is necessary to build a system of standard operating conditions to support the algorithm design. The operating conditions of the marine nuclear power system are shown in the following table, and each operating condition is an operating model. state, corresponding to an operating domain and a standard operating condition.
表2核动力系统运行工况体系Table 2 Nuclear power system operating condition system
2)单模态标准运行工况确定策略2) Determination strategy for single-modal standard operating conditions
核动力系统是一个由核物理、传热、热工水力、机械运动几种物理过程相互耦合组成的动态平衡系统,加之船用核动力系统出作为主推进动力外,还有各种其他负载,导致即使在同一运行工况下,其状态特征参数也不可能相同,但是单一运行工况下核动力系统运行状态之间逻辑距离比较接近,通常会构成一个运行区域,选取的理想标准运行工况应该在这个运行区域的逻辑中心位置上。在遴选单模态标准工况时主要采取两种策略相结合的方式,第一种方式是模拟试验,是利用核动力系统配套的全范围模拟器依据典型工况运行参数系统运行状态作为标准状态;第二种方式是数据分析,是基于系统历史运行数据,计算获取样本中某一运行工况下所有运行数据样本之间的逻辑距离,并找出其中与其他所有运行数据样本距离之和最小的一个运行数据样本实例,作为表征该运行工况(运行模态)的标准运行工况样本,也可以直接计算运行区域内数据样本集的逻辑中心作为标准运行工况,核动系统运行工况的运行域与运行状态数据样本关系示意图如图4所示。The nuclear power system is a dynamic balance system composed of several physical processes of nuclear physics, heat transfer, thermal hydraulics, and mechanical motion coupled with each other. Even under the same operating condition, the state characteristic parameters cannot be the same, but the logical distance between the operating states of the nuclear power system under a single operating condition is relatively close, which usually constitutes an operating area, and the ideal standard operating condition selected should be at the logical center of this operating area. In the selection of single-modal standard operating conditions, a combination of two strategies is mainly adopted. The first method is the simulation test, which uses the full-scale simulator of the nuclear power system to take the operating state of the system according to the operating parameters of typical operating conditions as the standard state. ; The second method is data analysis, which is based on the historical operating data of the system, calculates the logical distance between all operating data samples under a certain operating condition in the sample, and finds out the smallest sum of distances from all other operating data samples. An instance of an operating data sample of , as a standard operating condition sample characterizing the operating condition (operating mode), or it can directly calculate the logical center of the data sample set in the operating area as a standard operating condition, and check the system operating condition The schematic diagram of the relationship between the running domain and the running state data samples is shown in Figure 4.
3、基于逻辑距离的表征函数设计3. Design of characterization function based on logical distance
上述过程中,采用数据分析方法确定标准运行工况,需要计算不同运行数据样本之间的逻辑距离,在后面步骤中确定安全运行域边界或者辨识样本是否安全时,需要计算当前样本或者支持向量样本与标准运行工况之间的逻辑距离,可见,逻辑距离计算函数设计是算法设计的关键环节,需要设计一个具备良好问题适用性和鲁棒性的逻辑距离计算函数,以实现异常检测算法。In the above process, the data analysis method is used to determine the standard operating conditions, and the logical distance between different operating data samples needs to be calculated. When determining the boundary of the safe operating domain or identifying whether the sample is safe in the following steps, the current sample or support vector sample needs to be calculated. From the logical distance between the standard operating conditions, it can be seen that the design of the logical distance calculation function is the key link in the algorithm design. It is necessary to design a logical distance calculation function with good problem applicability and robustness to realize the anomaly detection algorithm.
1)稳态功率运行工况的逻辑距离计算函数。设计的逻辑距离计算函数是用来定量分析系统不同稳定运行状态之间逻辑关系远近的度量函数。表征运行状态的物理量是遴选出来的23个状态特征参数组成的多维向量,需要构建一个计算两个23维向量之间的逻辑距离的计算函数。因为23个特征参数,在辨识系统运行状态时,重要度并不相同,所以在计算时每个参数都需要进行加权处理,权值的选择需要依据运行经验、数值试验综合确定。基于上述分析,以欧氏距离为原始设计,定义了一种加权的欧氏距离算法用于计算系统稳态间的逻辑距离,具体如下。1) The logical distance calculation function of the steady-state power operating condition. The designed logical distance calculation function is a measure function used to quantitatively analyze the logical relationship between different stable operating states of the system. The physical quantity characterizing the running state is a multi-dimensional vector composed of 23 selected state characteristic parameters, and a calculation function to calculate the logical distance between two 23-dimensional vectors needs to be constructed. Because the importance of 23 characteristic parameters is not the same when identifying the operating state of the system, each parameter needs to be weighted during calculation, and the selection of weights needs to be comprehensively determined based on operating experience and numerical experiments. Based on the above analysis, with Euclidean distance as the original design, a weighted Euclidean distance algorithm is defined to calculate the logical distance between system steady states, as follows.
稳态功率运行时,设需要计算逻辑距离的两个多维向量分别为:When the steady-state power is running, the two multi-dimensional vectors that need to calculate the logical distance are as follows:
A=[a1,a2,a3,a4,…,a23]B=[b1,b2,b3,…,b23]A=[a1,a2,a3,a4,...,a23]B=[b1,b2,b3,...,b23]
对应的权值向量为:The corresponding weight vector is:
C=[c1,c2,c3,c4,…,c23]C=[c1,c2,c3,c4,...,c23]
c1+c2+,…,+c23=100c1+c2+,…,+c23=100
定义的加权欧氏距离计算函数为:The defined weighted Euclidean distance calculation function is:
2)对于启、停堆和工况切换瞬态过程,因特征参数处于单向性变化中,其运行状态是一个面向过程的时空区域,单纯通过特征参数在单一时刻的表征值难以体现系统的动态变化过程,因此,在该类工况下计算逻辑距离,不仅要考虑其表征值,还要考虑其时间序列的差分值,为降低问题的复杂度,因参数变化对系统运行状态变化贡献不一,所以并不需要将所有特征参数的差分值都引入到特征向量,可仅选取部分关键参数提取差分值即可,具体需根据实际问题进行分析。2) For the transient process of startup, shutdown and operating condition switching, because the characteristic parameters are in a one-way change, the operating state is a process-oriented space-time region, and it is difficult to reflect the system’s behavior simply by the characteristic values of the characteristic parameters at a single moment. Dynamic change process, therefore, to calculate the logical distance under this kind of working conditions, not only its characterization value, but also the difference value of its time series should be considered. First, it is not necessary to introduce the difference values of all the feature parameters into the feature vector, but only some key parameters can be selected to extract the difference values, which need to be analyzed according to the actual problem.
因此,在瞬态运行过程中,用于计算运行状态逻辑距离的特征向量表示为:Therefore, during transient operation, the eigenvector used to calculate the logical distance of the operating state is expressed as:
A=[a1,a2,a3,a4,…,a23,a1′,a2′,a3′,…,ak′],k<23A=[a1, a2, a3, a4, ..., a23, a1', a2', a3', ..., ak'], k<23
B=[b1,b2,b3,b4,…,b23,b1′,b2′,b3′,…,bk′],k<23B=[b1, b2, b3, b4,..., b23, b1', b2', b3',..., bk'], k<23
对应权值向量为:The corresponding weight vector is:
C=[c1,c2,c3,c4,…,c23,c1′,c2′,c3′,…,ck′],k<23C=[c1, c2, c3, c4, ..., c23, c1', c2', c3', ..., ck'], k<23
c1+c2+…,c23=100c1+c2+…, c23=100
c1′+c2′+c3′,…,+ck′=100c1'+c2'+c3',...,+ck'=100
关于瞬态运行过程逻辑距离,需定义为特征参数逻辑距离和特征参数变化速率(差分)逻辑距离两个计算函数,具体如下所示:Regarding the logical distance of the transient operation process, it needs to be defined as two calculation functions: the logical distance of the characteristic parameter and the logical distance of the rate of change of the characteristic parameter (difference), as shown below:
k<23k<23
4、安全运行域表征算法体系与归属阈值确定4. Determining the algorithm system for the characterization of the safe operation domain and the attribution threshold
对于系统单个运行工况的安全运行域表征算法,包括三个要素:标准运行工况及其特征向量,逻辑距离计算函数和安全域归属阈值,安全域归属阈值是在前两个要素设计完成后,通过系统正常运行数据样本训练确定的,并需要通过仿真计算或者收集的异常数据样本的检验验证,才能最终确定。For the safe operating domain characterization algorithm of a single operating condition of the system, it includes three elements: standard operating conditions and their eigenvectors, a logical distance calculation function and a safe domain belonging threshold. The safe domain belonging threshold is after the design of the first two elements is completed. , which is determined through the training of data samples of normal operation of the system, and needs to be verified through simulation calculation or the inspection and verification of abnormal data samples collected before it can be finally determined.
与核动力系统标准运行工况体系相对应,核动力系统安全运行域表征算法也是由不同模态构成的算法体系,具体如下表所示。Corresponding to the nuclear power system standard operating condition system, the nuclear power system safe operating domain characterization algorithm is also an algorithm system composed of different modes, as shown in the following table.
表3标准运行工况体系及逻辑距离计算函数Table 3 Standard operating condition system and logical distance calculation function
5、算法数值试验验证5. Algorithm numerical test verification
为验证算法的有效性,选取某船用核动力装置历史运行数据作为实验对象,分别针对常用的行进一稳态运行工况和运行瞬变节奏较快的快速升负荷工况进行了数值试验验证。In order to verify the effectiveness of the algorithm, the historical operating data of a marine nuclear power plant was selected as the experimental object, and numerical tests were carried out for the commonly used traveling-steady-state operating conditions and the fast-rising load conditions with fast operating transient rhythm.
(1)稳态工况异常检测数值试验验证(1) Numerical test verification of abnormal detection under steady state conditions
1)数值试验设定1) Numerical test setting
a、数据样本构造a. Data sample structure
首先,使用实际系统运行数据对设计算法进行检验,主要验证算法对稳态工况下随机偶发性偏离、扰动、异常的检测能力,测试稳态工况异常检测算法的训练样本和测试样本,通过抽取对象核动力装置在不同燃耗阶段的行进一状态数据构建,每个阶段的采样窗口宽度为5000秒,共选取了10个阶段的数据样本,选取的数据样本的部分关键特征参数时域曲线如图5所示。First, use the actual system operation data to test the design algorithm, mainly to verify the algorithm's ability to detect random and occasional deviations, disturbances, and anomalies under steady-state conditions, and to test the training samples and test samples of the anomaly detection algorithm under steady-state conditions. The travel-state data of the object nuclear power plant in different burnup stages are extracted to construct, the sampling window width of each stage is 5000 seconds, a total of 10 stages of data samples are selected, and the time-domain curves of some key characteristic parameters of the selected data samples are As shown in Figure 5.
从上述特征参数分布可以看出,特定工况的运行状态参数呈现区域聚集现象明显,可构造典型的安全运行区域。From the above characteristic parameter distribution, it can be seen that the operating state parameters of a specific working condition exhibit obvious regional aggregation phenomenon, and a typical safe operating area can be constructed.
其次,针对设备损伤引起的系统偏移和故障,因实际系统无法获得完整相关检测样本,因此通过与实际系统相对应的仿真模拟机获取测试样本,设定的设备损伤性异常和故障测试样本类型包括蒸汽发生器传热管微小破口、堆左回路冷端微小破口,安全阀卡开,二回路蒸汽管道破口,单根控制棒失控抽出等5类,涵盖堆芯,一、二回路,主辅系统和薄弱环节。Secondly, for the system offset and fault caused by equipment damage, because the actual system cannot obtain complete relevant test samples, the test samples are obtained through the simulation machine corresponding to the actual system, and the set equipment damage abnormality and fault test sample types Including micro-breaks in the heat transfer tubes of the steam generator, micro-breaks at the cold end of the left loop of the reactor, stuck open safety valves, breaks in the secondary-circuit steam pipes, and uncontrolled extraction of a single control rod, covering the reactor core, primary and secondary loops , main and auxiliary systems and weak links.
b、标准运行工况确定与参数权值分配b. Determination of standard operating conditions and allocation of parameter weights
取表征系统运行性状态的各特征参数所有样本数据值的平均值作为标准值,以之为分量得到行进一的标准运行工况特征向量;并依据遴选参数的重要程度不同,以及重复数值验证反馈,确定计算逻辑距离时各分量参数分配权值,具体如下表所示,因设计算法主要检测核动力一回路系统的异常和故障,故一回路及与一回路直接相关的参数权值较高。Take the average value of all sample data values of each characteristic parameter that characterizes the operational state of the system as the standard value, and use it as the component to obtain the characteristic vector of the standard operating condition of march one; , determine the distribution weights of each component parameter when calculating the logical distance, as shown in the following table, because the design algorithm mainly detects the abnormality and failure of the nuclear power primary circuit system, so the primary circuit and the parameters directly related to the primary circuit have higher weights.
表4行进一标准运行工况及参数权值分配Table 4. Standard operating conditions and parameter weight assignment
2)数值试验结果分析2) Analysis of numerical test results
a、使用实际系统数据对算法在稳态工况下随机偶发性偏离、扰动、异常的检测能力的验证结果如下图6所示。a. The verification results of the algorithm's ability to detect random and occasional deviations, disturbances, and anomalies under steady-state conditions using actual system data are shown in Figure 6 below.
从图中可以看出,正常运行时,实际系统运行状态与标准运行工况间的逻辑距离时域分布均匀,并呈现典型带状聚集现象。发生偶发性异常、故障后,其取值呈现阶跃性上扬,异常征兆明显,只要依据实际需求选取合适的阈值,便可以获得敏感性强且可线性调节的在线异常检测算法。通过对检测出的异常检测样本进行人工勘验,均确定了异常或偏离的真实存在性;对剩余未检测出异常的正常检测样本进行人工勘验,结果均为正常,具体如图7所示。因此,无论是从误检率角度,还是漏检率角度分析,该方法都具有优异的表现。It can be seen from the figure that during normal operation, the logical distance between the actual system operating state and the standard operating condition is uniformly distributed in the time domain, and presents a typical band-like aggregation phenomenon. After an occasional abnormality or failure occurs, its value shows a step increase, and the abnormal symptoms are obvious. As long as an appropriate threshold is selected according to the actual demand, an online abnormality detection algorithm with strong sensitivity and linear adjustment can be obtained. Through manual inspection of the detected abnormal detection samples, the true existence of the abnormality or deviation is confirmed; the manual inspection of the remaining normal detection samples with no abnormality detected, the results are all normal, as shown in Figure 7. . Therefore, the method has excellent performance both from the perspective of false detection rate and missed detection rate.
b、使用全范围模拟机计算获得的破坏性损伤数据样本,对异常检测算法的验证效果如下图7所示。从全范围模拟及计算运行数据数据分析,计算其与标准运行工况的逻辑距离,不难发现其值在时域上呈现典型周期性的高频干扰,其稳定性反而低于实际运行数据,因此,从稳定性角度判断,有理由相信基于安全运行域的异常检测方法对于实际系统的应用效果会优于模拟器。b. Using the destructive damage data samples obtained by the full-scale simulation computer calculation, the verification effect of the anomaly detection algorithm is shown in Figure 7 below. From the data analysis of full-scale simulation and calculation operation data, and calculating the logical distance between it and the standard operating conditions, it is not difficult to find that its value presents typical periodic high-frequency interference in the time domain, and its stability is lower than the actual operation data. Therefore, from the perspective of stability, it is reasonable to believe that the application effect of the anomaly detection method based on the safe operation domain for the actual system is better than that of the simulator.
为验证算法对破坏性损伤故障的检测能力,设置了5类破坏性故障检测数值试验,从测试结果分析,针对5类故障,设计的安全运行域算法以及定义的逻辑距离表征技术,均能够获得明显的故障征兆。In order to verify the detection ability of the algorithm for destructive damage faults, five types of destructive fault detection numerical tests are set up. From the analysis of the test results, the designed safe operation domain algorithm and the defined logical distance characterization technology can be obtained for the five types of faults. Obvious signs of failure.
检测算法针对不同类型的故障,其在逻辑距离上表征出的故障特征存在差异,比如一回路承压边界破损故障,其经高频滤波后的逻辑距离是逐渐偏离正常工况,并线性递增的,其时域值呈现为一条直线,可以用直线的斜率来表征故障程度(破口尺寸),即可以使用其斜率来定义异常检测阈值;反应性事故发生时,逻辑距离呈现出阶跃性扩大现象;二回路蒸汽管道破口时,系统热阱异常,逻辑距离会呈现整体性偏离正常区域,同时随着二回路工质的流失偏离会逐渐增大,后两者可以直接在时域内确定检测阈值。The detection algorithm is aimed at different types of faults, and there are differences in the fault characteristics represented by the logical distance. For example, the pressure boundary of the primary circuit is damaged. The logical distance after high-frequency filtering gradually deviates from the normal working condition and increases linearly. , its time domain value is presented as a straight line, and the slope of the straight line can be used to characterize the degree of failure (crack size), that is, its slope can be used to define the abnormal detection threshold; when a reactive accident occurs, the logical distance shows a step-by-step expansion Phenomenon; when the secondary circuit steam pipeline breaks, the system heat sink is abnormal, the logical distance will deviate from the normal area as a whole, and the deviation will gradually increase with the loss of the secondary circuit working medium. The latter two can be directly detected in the time domain. threshold.
试验结果同时表明,算法在检测破坏者性故障时,对不同类型故障在逻辑距离上会呈现出的不同表征现象,这可以作为异常检测后面一个环节故障类型的辨识与诊断的依据,在后续研究工作中将会进行针对性深入分析。The test results also show that when the algorithm detects destructive faults, different types of faults will show different characterization phenomena in the logical distance, which can be used as the basis for the identification and diagnosis of fault types in the next link of abnormal detection. In the follow-up research Targeted in-depth analysis will be carried out in the work.
(2)变换工况状态异常检测数值试验验证(2) Numerical test verification of abnormal state detection in changing conditions
1)数值试验设定1) Numerical test setting
针对变工况情况下算法的验证,系统选取了参数变化剧烈的快速升负荷过程作为对象工况,状态参数及其分配权值仍然选取与稳态工况相同方案。在计算状态参数逻辑距离时,遴选的标准工况是升功率前的初始工况对应的标准运行工况,即行进一标准运行工况;在计算参数变化速率逻辑距离时,设定计算参数变化速率(差分)的时间间隔是1秒钟,选取的标准工况分量是训练样本中对应特征参数所有样本差分值的中位数,如下式所示。For the verification of the algorithm under variable working conditions, the system selects the rapid load-rising process with drastic parameter changes as the object working condition, and the state parameters and their distribution weights still select the same scheme as the steady-state working condition. When calculating the logical distance of the state parameter, the selected standard operating condition is the standard operating condition corresponding to the initial operating condition before the power increase, that is, a standard operating condition; when calculating the logical distance of the parameter change rate, set the calculation parameter change The time interval of the rate (difference) is 1 second, and the selected standard condition component is the median of all sample difference values of the corresponding feature parameters in the training sample, as shown in the following formula.
a′i标准=(MAX(b′i)+Min(b′i))/2,i∈(1,2,…,23)a'i criterion=(MAX(b'i)+Min(b'i))/2,i∈(1,2,...,23)
2)数值试验结果分析2) Analysis of numerical test results
从实际系统历史运行数据库中提取典型从行进一快速升负荷的数据样本作为训练样本,首先按照设计的状态参数逻辑距离计算算法,计算训练样本与行进一标准运行工况之间的逻辑距离dAB,得到各训练样本的状态参数表征值逻辑距离如图8所示,图中在(300至400)*0.125秒时刻的阶跃变化是由于主冷却剂泵低速切为高速造成的。从图中可以看出,快速升负荷工况,随着系统负荷的不断提升,与标准运行工况的逻辑距离不断增加,并且明显具有区域性聚集表现,能够形成有效的安全运行域,表明基于安全运行域的异常检测算法可行,但是直接由状态参数表征值逻辑距离一个参数判定系统异常,信息并不充足,需要一些辅助参数才能得到结论。From the historical operation database of the actual system, the typical data samples of the running-to-fast-rising load are extracted as the training samples. First, according to the designed state parameter logical distance calculation algorithm, the logical distance dAB between the training sample and the running-standard operating condition is calculated, The logical distance of the state parameter characterization value of each training sample is obtained as shown in Figure 8. The step change at (300 to 400)*0.125 seconds in the figure is caused by the main coolant pump being switched from low speed to high speed. It can be seen from the figure that the rapid load increase condition, with the continuous increase of the system load, the logical distance from the standard operating condition continues to increase, and it has obvious regional aggregation performance, which can form an effective safe operation domain, indicating that based on The anomaly detection algorithm in the safe operation domain is feasible, but the system anomaly is determined directly by a parameter of the state parameter characterization value logical distance, and the information is not sufficient, and some auxiliary parameters are needed to reach the conclusion.
快速升负荷过程系统的运行状态不仅与特征参数的表征值有关,更多的信息其实存在于特征参数的变化速率中,因此从特征参数的差分数据中可能更容易检索到故障数据,图9展示了特征参数查分值与标准查分值之间的逻辑距离d’AB分布情况,及故障状态特征,从图中可以看出,在差分逻辑距离空间里,只要设置合适阈值,可以很容易检索出系统异常。验证表明选择合适的表征参数,设计算法可以有效检测出系统运行瞬态工况下发生的异常和故障。The operating state of the system during the rapid load escalation process is not only related to the characteristic values of the characteristic parameters, but more information actually exists in the rate of change of the characteristic parameters, so it may be easier to retrieve the fault data from the differential data of the characteristic parameters, as shown in Figure 9. The distribution of the logical distance d'AB between the characteristic parameter score value and the standard score value, and the fault state characteristics, it can be seen from the figure that in the differential logical distance space, as long as the appropriate threshold is set, it can be easily retrieved A system exception occurred. The verification shows that by selecting appropriate characterization parameters, the design algorithm can effectively detect the abnormality and faults that occur under the transient operating conditions of the system.
本发明图10、11体现了不同类型故障,检测结果呈现出不同特征。Figures 10 and 11 of the present invention show different types of faults, and the detection results show different characteristics.
本发明遴选23个状态表征技术参数的数值分析依据。The invention selects the numerical analysis basis of 23 state characterization technical parameters.
本发明对比图10、11、9与图12(基于神经网络类的故障诊断结果)可看出神经网络类的异常检测及故障诊断算法仅能输出一个匹配最相似的分类号,无法体现故障的渐变过程,可解释性差;而本发明可以体现故障程度和不同故障的相近性,以及故障渐变过程,或者偶发时机,可解释性强。Comparing Figures 10, 11, 9 and Figure 12 (fault diagnosis results based on neural network) in the present invention, it can be seen that the abnormal detection and fault diagnosis algorithm of neural network can only output a classification number that matches the most similar, and cannot reflect the fault. The gradual change process has poor interpretability; while the present invention can reflect the degree of failure and the similarity of different faults, as well as the gradual change process of the fault, or the occasional time, and the interpretability is strong.
与基于故障特征分析的传统异常检测方法相比较,本案检测范围广,且不受故障样本稀少的限制。如图12,传统异常检测只能针对已知的经典的单一类型故障实施检测,或者针对单一参数进行系统异常检测,且通常还需要使用仿真模拟器计算获得或者补充故障训练样本,而模拟系统与实际系统存在固有计算误差,会影响诊断方法的有效性;对照图7、8、9,本案基于实际系统构建安全运行域,贴近实际,且可同时检测各种类型的异常,可检测未知的异常和故障,拓展了检测范围,降低了漏检率,同时可根据需求通过阈值调节检测力度,还有效规避了故障样本稀少且分布不均的现实问题。Compared with traditional anomaly detection methods based on fault feature analysis, this case has a wide detection range and is not limited by the scarcity of fault samples. As shown in Figure 12, traditional anomaly detection can only detect a known classic single type of fault, or perform system anomaly detection for a single parameter, and usually needs to use a simulation simulator to calculate or supplement the fault training samples, while the simulation system is similar to There are inherent calculation errors in the actual system, which will affect the effectiveness of the diagnosis method; compared with Figures 7, 8, and 9, this case builds a safe operation domain based on the actual system, which is close to reality, and can detect various types of anomalies at the same time, and can detect unknown anomalies. The detection range is expanded and the missed detection rate is reduced. At the same time, the detection intensity can be adjusted through the threshold value according to the requirements, and the practical problem of rare and uneven distribution of fault samples can be effectively avoided.
在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上;术语“上”、“下”、“左”、“右”、“内”、“外”、“前端”、“后端”、“头部”、“尾部”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”等仅用于描述目的,而不能理解为指示或暗示相对重要性。In the description of the present invention, unless otherwise stated, "plurality" means two or more; the terms "upper", "lower", "left", "right", "inner", "outer" The orientation or positional relationship indicated by , "front end", "rear end", "head", "tail", etc. are based on the orientation or positional relationship shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, not An indication or implication that the referred device or element must have a particular orientation, be constructed and operate in a particular orientation, is not to be construed as a limitation of the invention. Furthermore, the terms "first," "second," "third," etc. are used for descriptive purposes only and should not be construed to indicate or imply relative importance.
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,都应涵盖在本发明的保护范围之内。The above are only specific embodiments of the present invention, but the protection scope of the present invention is not limited to this. Any person skilled in the art is within the technical scope disclosed by the present invention, and all within the spirit and principle of the present invention Any modifications, equivalent replacements and improvements made within the scope of the present invention should be included within the protection scope of the present invention.
Claims (9)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010646535.1A CN111881176B (en) | 2020-07-07 | 2020-07-07 | Anomaly detection method for marine nuclear power based on logical distance characterization of safe operation domain |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010646535.1A CN111881176B (en) | 2020-07-07 | 2020-07-07 | Anomaly detection method for marine nuclear power based on logical distance characterization of safe operation domain |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111881176A true CN111881176A (en) | 2020-11-03 |
| CN111881176B CN111881176B (en) | 2021-09-07 |
Family
ID=73150329
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010646535.1A Active CN111881176B (en) | 2020-07-07 | 2020-07-07 | Anomaly detection method for marine nuclear power based on logical distance characterization of safe operation domain |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111881176B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112684778A (en) * | 2020-12-24 | 2021-04-20 | 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) | Steam generator water supply system diagnosis method based on multi-source information reinforcement learning |
| CN114519382A (en) * | 2022-01-05 | 2022-05-20 | 哈尔滨工程大学 | Nuclear power plant key operation parameter extraction and abnormity monitoring method |
| WO2024216890A1 (en) * | 2023-04-19 | 2024-10-24 | 中广核研究院有限公司 | Nuclear reactor fault diagnosis method and apparatus, computer device, and storage medium |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000242328A (en) * | 1999-02-23 | 2000-09-08 | Toshiba Corp | Surveillance test apparatus and storage medium storing program |
| CN102298978A (en) * | 2011-05-17 | 2011-12-28 | 哈尔滨工程大学 | MFM (multilevel flow model)-based indeterminate fault diagnosis method for nuclear power plant for ship |
| CN106368816A (en) * | 2016-10-27 | 2017-02-01 | 中国船舶工业系统工程研究院 | Method for online abnormity detection of low-speed diesel engine of ship based on baseline deviation |
| CN107609313A (en) * | 2017-10-18 | 2018-01-19 | 哈尔滨工程大学 | A kind of passive safety system analysis method for reliability peculiar to vessel |
| CN107767975A (en) * | 2017-10-17 | 2018-03-06 | 中北大学 | A monitoring and fault diagnosis method for the quality and performance of critical components of nuclear power plants |
| CN108961696A (en) * | 2018-06-20 | 2018-12-07 | 中国船舶重工集团公司第七〇九研究所 | A kind of early warning system and method for early warning of ocean nuclear power platform |
| CN110569990A (en) * | 2019-08-02 | 2019-12-13 | 中国船舶重工集团公司第七一九研究所 | Operation and maintenance system and operation and maintenance method suitable for marine nuclear power platform |
| CN110738274A (en) * | 2019-10-26 | 2020-01-31 | 哈尔滨工程大学 | A data-driven fault diagnosis method for nuclear power plant |
-
2020
- 2020-07-07 CN CN202010646535.1A patent/CN111881176B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000242328A (en) * | 1999-02-23 | 2000-09-08 | Toshiba Corp | Surveillance test apparatus and storage medium storing program |
| CN102298978A (en) * | 2011-05-17 | 2011-12-28 | 哈尔滨工程大学 | MFM (multilevel flow model)-based indeterminate fault diagnosis method for nuclear power plant for ship |
| CN106368816A (en) * | 2016-10-27 | 2017-02-01 | 中国船舶工业系统工程研究院 | Method for online abnormity detection of low-speed diesel engine of ship based on baseline deviation |
| CN107767975A (en) * | 2017-10-17 | 2018-03-06 | 中北大学 | A monitoring and fault diagnosis method for the quality and performance of critical components of nuclear power plants |
| CN107609313A (en) * | 2017-10-18 | 2018-01-19 | 哈尔滨工程大学 | A kind of passive safety system analysis method for reliability peculiar to vessel |
| CN108961696A (en) * | 2018-06-20 | 2018-12-07 | 中国船舶重工集团公司第七〇九研究所 | A kind of early warning system and method for early warning of ocean nuclear power platform |
| CN110569990A (en) * | 2019-08-02 | 2019-12-13 | 中国船舶重工集团公司第七一九研究所 | Operation and maintenance system and operation and maintenance method suitable for marine nuclear power platform |
| CN110738274A (en) * | 2019-10-26 | 2020-01-31 | 哈尔滨工程大学 | A data-driven fault diagnosis method for nuclear power plant |
Non-Patent Citations (3)
| Title |
|---|
| 余刃等: "基于小波独立成分分析的核动力装置冗余传感器故障在线诊断方法研究 ", 《核动力工程》 * |
| 宋梅村等: "基于动态PCA的核动力装置传感器故障检测 ", 《武汉理工大学学报(交通科学与工程版)》 * |
| 陈进军等: "基于RBF神经网络的核动力装置故障诊断方法研究 ", 《热科学与技术》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112684778A (en) * | 2020-12-24 | 2021-04-20 | 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) | Steam generator water supply system diagnosis method based on multi-source information reinforcement learning |
| CN114519382A (en) * | 2022-01-05 | 2022-05-20 | 哈尔滨工程大学 | Nuclear power plant key operation parameter extraction and abnormity monitoring method |
| CN114519382B (en) * | 2022-01-05 | 2024-10-18 | 哈尔滨工程大学 | Nuclear power plant key operation parameter extraction and abnormality monitoring method |
| WO2024216890A1 (en) * | 2023-04-19 | 2024-10-24 | 中广核研究院有限公司 | Nuclear reactor fault diagnosis method and apparatus, computer device, and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111881176B (en) | 2021-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113792762B (en) | Chiller fault diagnosis method, system and medium based on Bayesian optimized LightGBM | |
| Li et al. | Diagnosis for multiple faults of chiller using ELM-KNN model enhanced by multi-label learning and specific feature combinations | |
| CN108376298A (en) | A kind of Wind turbines generator-temperature detection fault pre-alarming diagnostic method | |
| Cao et al. | PCA-SVM method with sliding window for online fault diagnosis of a small pressurized water reactor | |
| CN107092242B (en) | A kind of Industrial Process Monitoring method based on missing variable pca model | |
| CN115187832A (en) | Energy system fault diagnosis method based on deep learning and gram angular field image | |
| Peng et al. | An intelligent hybrid methodology of on-line system-level fault diagnosis for nuclear power plant | |
| CN111899905B (en) | A kind of fault diagnosis method and system based on nuclear power plant | |
| CN108051211A (en) | A kind of wind generator set main shaft holds temperature pre-warning diagnostic method | |
| CN114943281B (en) | Intelligent decision-making method and system for heat pipe cooling reactor | |
| CN111881176A (en) | Anomaly detection method for marine nuclear power based on logical distance characterization of safe operation domain | |
| CN112016251A (en) | Nuclear power device fault diagnosis method and system | |
| CN112036087A (en) | A multi-strategy fusion method and system for fault diagnosis of nuclear power key equipment | |
| Zhou et al. | Structural health monitoring of offshore wind power structures based on genetic algorithm optimization and uncertain analytic hierarchy process | |
| CN108446529A (en) | Organic rankine cycle system fault detection method based on broad sense cross-entropy-DPCA algorithms | |
| CN108388234A (en) | A kind of fault monitoring method dividing changeable gauge block pca model based on correlation | |
| CN111648992A (en) | Gas turbine compressor fault identification and early warning method | |
| CN111797533B (en) | A method and system for abnormal detection of operating parameters of nuclear power plant | |
| Yu et al. | A continuous learning monitoring strategy for multi-condition of nuclear power plant | |
| CN110880024A (en) | Fault identification method and system for nonlinear process based on discriminant kernel slow feature analysis | |
| CN105388884A (en) | Alarm system for detecting leakage fault of heat supply network based on identification algorithm driven by data and method | |
| Li et al. | Diagnosis for the refrigerant undercharge fault of chiller using deep belief network enhanced extreme learning machine | |
| CN103617105A (en) | Self-adaptation multilevel flow model equipment diagnosis method based on data driving | |
| CN112036496A (en) | Nuclear power device fault diagnosis method and system | |
| CN110289112B (en) | Nuclear power plant health state diagnosis method based on hierarchical analysis and fuzzy evaluation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |
















