CN105512812A - 基于动态仿真模型的核电厂设备故障预警分析方法及系统 - Google Patents
基于动态仿真模型的核电厂设备故障预警分析方法及系统 Download PDFInfo
- Publication number
- CN105512812A CN105512812A CN201510872454.2A CN201510872454A CN105512812A CN 105512812 A CN105512812 A CN 105512812A CN 201510872454 A CN201510872454 A CN 201510872454A CN 105512812 A CN105512812 A CN 105512812A
- Authority
- CN
- China
- Prior art keywords
- early warning
- monitoring equipment
- warning analysis
- nuclear power
- power plant
- 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.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 164
- 238000005094 computer simulation Methods 0.000 title claims abstract description 47
- 238000012544 monitoring process Methods 0.000 claims abstract description 127
- 238000000034 method Methods 0.000 claims abstract description 60
- 230000008569 process Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 16
- 230000007246 mechanism Effects 0.000 claims description 12
- 230000010354 integration Effects 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000007812 deficiency Effects 0.000 description 4
- 230000036541 health Effects 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000006378 damage Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 208000033999 Device damage Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G07—CHECKING-DEVICES
- G07C—TIME OR ATTENDANCE REGISTERS; REGISTERING OR INDICATING THE WORKING OF MACHINES; GENERATING RANDOM NUMBERS; VOTING OR LOTTERY APPARATUS; ARRANGEMENTS, SYSTEMS OR APPARATUS FOR CHECKING NOT PROVIDED FOR ELSEWHERE
- G07C3/00—Registering or indicating the condition or the working of machines or other apparatus, other than vehicles
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种基于动态仿真模型的核电厂设备故障分析方法及系统。该系统包括:核电站实时/历史数据库,用于提供动态仿真模型正常工作所需基本参数;基于核电站全范围模拟机构建的动态仿真模型,用于提供各监测设备的基准参数值;实时数据采集中心,用于测量各监测设备的实时运行参数值;动态预警分析平台,用于将各监测设备的实时运行参数值及相应的基准参数值输入趋势预警模型,在趋势预警模型中应用综合预警分析策略对各监测设备进行动态预警分析,并对经综合预警分析策略判定为故障设备的监测设备启动预警。本发明可提早识别核电厂运行中的故障设备并启动预警,以防止故障设备的进一步损害,避免核电厂事故的发生,提高核电厂的运营效率。
Description
技术领域
本发明属于核电站设备监控领域,更具体地说,本发明涉及一种基于动态仿真模型的核电厂设备故障预警方法及系统。
背景技术
核电厂设备的运行状况对整个核电厂的正常运行生产都具有十分重要的意义。一方面,当关键设备发生故障时,可能对整个系统及其他相关设备产生重大影响,造成非计划停堆,降低生产效率,甚至可能造成相关排放物排放超标。另一方面,由于生产环境的限制,核电厂关键设备一般要求具备高温高压、耐腐蚀、抗辐射等特性,这就决定了设备制造成本高,维修成本高,生产周期长的特点。
目前,核电厂主要采用的设备报警方法为高低限预警方法、高低线预警方法及基于机理模型的预警方法。
其中,高低线报警方法的技术缺陷在于:只有当故障测点达到预警线时才触发报警,此时,故障信息已经比较明显,设备损坏已比较严重,其不能提早发现异常设备,因而无法保障核电厂的安全稳定运行。
高低线预警方法在高低线报警方法的基础上增加了一定的裕度。但是,高低线预警方法并不能提早发现故障设备,不能从根本上解决设备故障预警的问题。
基于机理模型的预警方法通过对设备的运行机理进行分析,建立设备精确的数学模型,将模型产生的结果与实际测量值进行比较,产生残差值,通过状态估计等方法对残差值进行分析,以此评估设备是否发生故障。该基于机理模型的预警方法的技术缺陷在于:一方面,该方法需要对系统进行线性化处理,随着现代大生产的发展,现代设备,尤其核电关键设备的结构趋于复杂,耦合程度越来越高,采用机理模型的方法会造成模型精度不高,误差增大。另一方面,建模成本高,易受现场环境噪声的干扰。
发明内容
本发明的目的在于:针对现有技术的上述缺陷,提供一种能提早且精准地识别核电厂故障设备,提早启动预警,以确保核电厂安全稳定运行的基于动态仿真模型的核电厂设备故障预警方法及系统。
为了实现上述发明目的,本发明提供了一种基于动态仿真模型的核电厂设备故障预警方法,该方法包括如下步骤:
基于核电站全范围模拟机构建动态仿真模型;
通过动态仿真模型获取表征核电厂监测设备健康运行的基准参数值;
将采集的各个监测设备的实时运行参数值及其相应的基准参数值输入趋势预警模型,根据趋势预警模型并应用综合预警分析策略对核电厂各台监测设备进行动态预警分析,并对经综合预警分析策略判定为故障设备的监测设备启动预警。
本发明另提供有一种基于动态仿真模型的核电厂设备故障预警系统,该系统包括:
用作数据接口的核电站实时/历史数据库,用于为动态仿真模型提供其正常工作所需的基本参数;
基于核电站全范围模拟机构建的动态仿真模型,该动态仿真模型根据输入的确保其正常工作所需的基本参数计算及输出表征核电厂各监测设备健康运行的基准参数值;
实时数据采集中心,用于测量各监测设备的实时运行参数值;
动态预警分析平台,用于将该实时数据采集中心所测监测设备的实时运行参数值及相应的基准参数值输入趋势预警模型,根据趋势预警模型并应用综合预警分析策略对各监测设备进行动态预警分析,并对经综合预警分析策略判定为故障设备的监测设备启动预警。
应用本发明基于动态仿真模型的核电厂设备故障预警方法及系统,一方面,可提早且精准地识别核电厂监测设备中运行异常的监测设备,在该类故障设备的损坏程度尚不足以触发报警之前提前启动预警,以便指导并及时安排检修人员对故障设备进行点修或关停处理,可防止故障设备的进一步损坏,避免核电厂事故的发生,确保核电厂安全稳定运行,提高核电厂的运营效率;另一方面,本发明还可将反映监测设备是否出现故障及故障严重程度的预警因子以动态的图形化形式在趋势预警模型中进行实时显示,以便核电厂设备管理人员据此直观地了解现场监测设备的运行状态,合理安排及指导巡检人员的巡检工作。
附图说明
下面结合附图和具体实施方式,对本发明基于动态仿真模型的核电厂设备故障预警方法及系统进行详细说明,其中:
图1为本发明较佳实施例提供的基于动态仿真模型的核电厂设备故障预警方法的流程图;
图2为图1的故障预警方法中包括的综合动态预警分析方法的流程图;
图3为图2的综合动态预警分析方法中包括的残差预警分析方法的流程图;
图4为图2的综合动态预警分析方法中包括的曲线面积积分预警分析方法的流程图;
图5为图2的综合动态预警分析方法包括的波动值预警分析方法的流程图;
图6为本发明另一较佳实施例提供的基于动态仿真模型的核电厂设备故障预警系统的结构框图;
图7为图6所示的核电厂设备故障预警系统的动态预警分析平台的结构框图;
图8为图7所示的动态预警分析平台预警分析模块的结构框图。
具体实施方式
为了使本发明的发明目的、技术方案及其技术效果更加清晰,以下结合附图和具体实施方式,对本发明进一步详细说明。应当理解的是,本说明书中描述的具体实施方式仅仅是为了解释本发明,并非为了限定本发明。
如图1所示,针对现有技术中高低限预警法不能提早发现设备异常的技术缺陷,以及基于机理模型的预警法会造成模型精度不高,预警误差增大,建模成本较高,且易受环境噪声干扰的技术缺陷,本发明提出了一种基于动态仿真模型200的核电厂设备故障预警方法,该方法包括以下步骤:
S100、基于核电站全范围模拟机构建动态仿真模型200。核电厂仿真模型的建立,需要借助核电厂全范围仿真模拟机平台,借助仿真平台设计的动态仿真模型200具有很好的可靠性及实用性。在核电厂全范围模拟机平台上,利用图形化的搭建方式,建立所需的设备模型,进行编译后,将生成的设备源码转译为所需的语言,即可脱离核电厂全范围模拟机平台建立动态仿真模型200。
S200、通过所建立的动态仿真模型200模拟计算及获取表征核电厂各监测设备健康运行的基准参数值。在本发明中,为获取各监测设备的基准参数值,将核电厂实时/历史数据库作为接口,从中获取动态仿真模型200正常运转所需的基本参数,将上述基本参数输入动态仿真模型200,即可从动态仿真模型200输出获取核电厂各监测设备的基准参数值。
S300、将采集的各监测设备的实时运行参数及其基准参数值输入趋势预警模型,基于趋势预警模型并应用综合性预警分析策略对监测设备进行动态预警分析,并对经综合预警分析策略判定为故障设备的监测设备启动预警。
相比在核电厂中广泛应用的高低限报警方案或基于机理模型的故障预警方案而言,本发明核电厂设备故障预警方法通过对核电厂各监测设备的历史、实时数据进行建模分析,并应用综合预警分析策略对监测设备当前工作状态进行综合预警分析,提早且准确无误地识别故障监测设备,以在该监测设备的故障严重程度尚未达到触发高低限报警操作之前提早启动预警,提早排查核电厂中异常工作的监测设备,变事后处理为事前控制,以便及时通知维修人员对故障设备进行点修、维修,避免故障设备的彻底损坏,防止核电厂事故的发生,提高了核电厂运行的安全性及经济效益。
可选地,本发明还可包括位于步骤S200及步骤S300之间的步骤S230,该步骤S230如下:
在对核电厂监测设备执行动态预警分析之前,对所测的监测设备的实时运行参数值进行滤波处理。为减少现场干扰因素对传感器工作及其采集数据准确性的影响,减少误报率,优选地,通过二次函数平滑均值滤波算法对采集数据进行滤波处理,用以有效抑制由于环境噪声引起的干扰或毛刺现象,提高预警准确率。
请结合图2,在本发明核电厂设备故障预警方法中,上述步骤S300进一步包括:
步骤S301、应用残差预警分析、曲线面积积分预警分析及波动值预警分析相结合的综合预警分析策略对核电厂监测设备运行异常与否进行综合分析及判定,根据判定结果决定是否启动预警机制。本发明应用残差预警分析、曲线面积积分预警分析及波动值预警分析相结合的综合预警分析策略,上述三种预警方式相互补充,对基于其中一种预警分析策略所作设备故障预判作进一步核实,从而大大提高了预警准确率,避免了应用单一预警分析策略可能造成的错误预警。
步骤S302、以动态形式将监测设备相对于预设残差预警线、预设偏移预警线及预设波动预警线的残差、偏移及波动程度在趋势预警模型中进行实时显示。在本发明中,通过将反映监测设备故障大小及设备损坏程度的若干预警因子(例如相对于设定参考基准的残差、偏移率及波动率)以动态的图形化的形式进行实时显示,核电厂设备管理人员可直观地了解到各监测设备是否出现故障及故障大小,以便在核电厂设备故障预警系统的自动预警机制失灵时人工辅助介入,启动预警,进一步提高核电厂设备运行的安全性。
下面将对本发明综合预警分析策略中采用的残差预警分析方法作详细介绍。
请结合图3,该残差预警分析方法具体包括:
步骤S301a、对于待执行综合预警分析的监测设备,在设定的采样时间段内采集该台监测设备的实时运行参数;
步骤S302a、计算该台监测设备的实际运行参数值及该台监测设备的基准参数值之间的残差值;
步骤S303a、判断该残差值中的最大值是否超过设定的残差阈值;如是,则执行下一步骤S304a,否则,执行另一步骤S305a。
步骤S304a、做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
步骤S305a、判定该监测设备无故障运行,返回步骤S301a。
下面将对本发明综合预警分析策略中采用的曲线面积积分法作详细介绍。
请结合图4,该残差预警分析方法具体包括:
步骤S301b、基于该台监测设备的基准参数值及实时运行参数,在趋势预警模型中分别绘制表征该台监测设备的理想/健康运行状态及实际运行状态的基准曲线及实时运行曲线;
步骤S302b、将该条基准曲线作为参考线,在设定的采样时间段内以该参考曲线为基准对该条实时运行曲线作曲线面积积分运算,将运算所得曲面面积转换成直角三角形面积,根据该直角三角形面积计算该采样时间段内该台监测设备实际运行曲线相对于该基准线的偏移率;
步骤S303b、判断该偏移率是否超过设定的偏移率阈值(系统默认设定,或由管理人员凭经验自行设定)。如是,则执行步骤S304b,否则,执行步骤S305b。
步骤S304b、做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
步骤S305b、判定该台监测设备无故障运行,返回步骤S301b。
下面将对本发明综合预警分析策略中采用的波动值预警分析方法作详细介绍。
请结合图5,该波动值预警分析方法具体包括:
步骤S301c、针对待执行动态预警分析的监测设备,基于该监测设备的基准参数值在趋势预警模型中绘制表征该台监测设备在安全/理想状态下的基准曲线;
将该条基准曲线作为参考线,利用均方差原理计算该监测设备在指定的采样时间段内相对于该参考线的波动率;
步骤S303c、判断该波动率是否超过设定的波动率阈值;如是,则执行步骤S304c,否则,执行步骤S305c。
步骤S304c、做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
步骤S305c、判定该台监测设备无故障运行,返回步骤S301c。
本发明基于核电厂监测设备的历史数据及实测数据构建趋势预警模型,在趋势预警模型中应用上述综合预警分析策略对监测设备运行异常与否进行综合判定,不仅能够提早发现故障设备,还达到了极高的预警成功率。核电厂管理人员基于本发明综合预警分析结果提早安排故障设备检修或关停处理,能有效避免核电厂事故的发生,因此,本发明可为核电厂的安全运行提供了强有力保障。
另外,本发明还可将反映监测设备是否出现故障及故障严重程度的若干预警因子(例如监测设备实际运行曲线相对于设定参考基准的残差、偏移率及波动率)以动态的图形化的形式在趋势预警模型中进行实时显示,核电厂设备管理人员可直观地了解到各监测设备是否出现故障及故障严重程度,以便在核电厂设备故障预警系统的自动预警机制失灵时人工辅助介入,启动预警,进一步地提高了核电厂设备运行的安全性及稳定性。核电厂设备管理人员也能据此直观地了解现场监测设备的运行状态,合理安排巡检人员的巡检工作。
基于图1至图5所述的基于动态仿真模型200的核电厂故障预警方法,本发明还提出了一种基于动态仿真模型200的核电厂故障预警系统。
请结合图6,该核电厂故障预警系统包括:
作为数据接口的核电站实时/历史数据库100,用于为动态仿真模型200提供其正常工作所需的基本参数;
基于核电站全范围模拟机构建的动态仿真模型200;核电站全范围模拟机模拟真实核电站运行,基于核电站全范围模拟机构建的动态仿真模型200可根据取自核电站实时/历史数据库100的,保证其正常工作所需的基本参数模拟计算及输出表征核电厂各监测设备健康运行的基准参数值;
实时数据采集中心400,用于测量各监测设备的实时运行参数值;较佳地,该实时数据采集中心400可以包括用于测量核电厂各类监测设备的各种传感器。
动态预警分析平台300,用于将该实时数据采集中心400所测监测设备的实时运行参数值及相应的基准参数值输入趋势预警模型,在趋势预警模型中应用综合预警分析策略对核电厂监测设备进行动态预警分析,并对经综合预警分析策略判定为故障设备的监测设备启动预警。
可选地,该核电厂设备故障预警系统还包括滤波模块,用于在动态预警分析平台300执行综合性预警分析之前,对实时数据采集中心400获取的各种监测设备的实时运行参数进行滤波处理,以有效抑制由于环境噪声引起的干扰或毛刺现象,提高预警准确率。
请结合图7,该动态预警分析平台300进一步包括:
预警分析模块301,用于应用残差预警分析、曲面面积积分预警分析及波动值预警分析相结合的综合预警分析策略对监测设备运行异常与否进行综合判定,并对经综合预警分析策略判定为故障设备的监测设备启动预警机制;
预警结果显示模块302,用于将反映监测设备是否出现故障及故障严重程度的若干预警因子以动态的图形化的形式在趋势预警模型中进行实时显示,该预警因子包括该监测设备相对于相应参考基准(比如,预设的残差阈值、偏移率阈值及波动率阈值)的残差、偏移率及波动率。
请结合图8,该预警分析模块301进一步包括:
残差预警分析子模块3011,其进一步包括残差值计算单元及第一预警分析单元;其中:该残差值计算单元用于针对待执行预警分析的每台监测设备,将该实时数据采集中心400在设定的采样时间范围内采集的该台监测设备的实时运行参数值与经动态仿真模型200计算得到的该台监测设备的基准参数值进行比较,计算该实时运行参数值与该基准参数值之间的残差值;该第一预警分析单元用于判断该残差值中的最大值是否超过设定的残差阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实;
曲面面积积分预警分析子模块3012,其进一步包括偏移率计算单元及第二预警分析单元;其中,该偏移率计算单元用于针对待执行预警分析的每台监测设备,在趋势预警模型中分别绘制表征该台监测设备的理想/健康运行状态及实际运行状态的基准曲线及实时运行曲线,将该条基准曲线作为参考线,并在设定的采样时间范围内将该台监测设备的实时运行曲线以该参考线为基准作曲线面积积分运算,将曲面面积换算成直角三角形面积,根据该直角三角形面积分析该段采样时间内该设备运行相对于该参考线的偏移率;该第二预警分析单元用于判断该偏移值是否超过设定的偏移阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实;
波动值预警分析子模块3013,其进一步包括波动率计算单元及第三预警分析单元;其中:该波动率计算单元用于针对待执行预警分析的每台核电站设备,绘制表征该台监测设备在理想/健康状态下的基准曲线,将该条基准曲线作为参考线,利用均方差原理计算该设备在指定的采样时间范围内相对于该参考线的波动率;该第三预警分析单元用于判断该台监测设备实际运行过程中相对于该参考线的波动率是否超过设定的波动率阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
综上所述,相比现有的核电厂设备预警方案,在核电厂部署本发明基于动态仿真模型200的核电厂设备故障预警方法及系统,可提早且精准地甄别核电厂监测设备中运行存在异常的监测设备,在故障设备的损坏程度尚不足以触发报警之前提前启动预警,以便及时指导并安排检修人员对故障设备进行点修或关停处理,可防止故障设备的进一步损坏,避免核电厂事故的发生,由此确保了核电厂安全稳定地运行。本发明还可将反映监测设备是否出现故障及故障严重程度的预警因子以动态的图形化形式在趋势预警模型中进行实时显示,以便核电厂设备管理人员据此直观地了解现场监测设备的运行状态,合理安排及指导巡检人员的巡检工作。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(ReadOnlyMemory,ROM)或随机存储记忆体(RandomAccessMemory,RAM)等。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。此外,尽管本说明书中使用了一些特定的术语,但这些术语只是为了方便说明,并不对本发明构成任何限制。
Claims (11)
1.一种基于动态仿真模型的核电厂设备故障预警方法,其特征在于,该方法包括如下步骤:
基于核电站全范围模拟机构建动态仿真模型;
通过动态仿真模型模拟计算获得各监测设备的基准参数值;
将采集的各监测设备的实时运行参数值及其相应的基准参数值输入趋势预警模型,根据趋势预警模型并应用综合预警分析策略对核电厂各台监测设备进行动态预警分析,并对经综合预警分析策略判定为故障设备的监测设备启动预警。
2.根据权利要求1所述的核电厂设备故障预警方法,其特征在于,该方法还包括:
在对监测设备执行动态预警趋势分析之前,对所测监测设备的实时运行参数值进行滤波处理。
3.根据权利要求1所述的核电厂设备故障预警方法,其特征在于,该基于趋势预警模型并应用综合性预警分析策略对监测设备进行动态预警分析的步骤具体包括:
应用残差预警分析、曲面面积积分预警分析及波动值预警分析相结合的综合性预警分析策略对核电厂各监测设备运行异常与否进行综合判定。
4.根据权利要求3所述的核电厂设备故障预警方法,其特征在于:
该残差预警分析包括如下步骤:
在设定的采样时间段内采集监测设备的实时运行参数值,计算该监测设备的实际运行参数值及其基准参数值之间的残差值;
判断该残差值中的最大值是否超过设定的残差阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实;
该曲面面积积分预警分析包括如下步骤:
在趋势预警模型中分别绘制反映监测设备的基准参数值及实时运行参数值的基准曲线及实时运行曲线;
将该条基准曲线作为参考线,并在设定的采样时间段内,以该参考线为基准对该条实时运行曲线作曲线面积积分运算;
将曲面面积换算成直角三角形面积,根据该直角三角形面积分析该段采样时间段内该台监测设备实际运行相对于该参考线的偏移率;
判断该偏移率是否超过设定的偏移率阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实;
该波动值预警分析包括如下步骤:
在趋势预警模型中绘制反映监测设备的基准参数值的基准曲线,将该条基准曲线作为参考线,利用均方差原理计算该台监测设备在指定的采样时间段内相对于该参考线的波动率;
判断该监测设备相对于该参考线的波动率是否超过设定的波动率阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
5.根据权利要求1所述的核电厂设备故障预警方法,其特征在于,该方法还包括:
将反映监测设备是否出现故障及故障严重程度的若干预警因子以动态的图形化的形式在趋势预警模型中进行实时显示,该预警因子包括该监测设备相对于相应参考基准的残差、偏移率及波动率。
6.一种基于动态仿真模型的核电厂设备故障预警系统,其特征在于,该系统包括:
用作数据接口的核电站实时/历史数据库,用于为动态仿真模型提供其正常工作所需的基本参数;
基于核电站全范围模拟机构建的动态仿真模型,该动态仿真模型根据输入的确保其正常工作所需的基本参数模拟计算及输出各监测设备的基准参数值;
实时数据采集中心,用于测量各监测设备的实时运行参数值;
动态预警分析平台,用于将该实时数据采集中心所测监测设备的实时运行参数值及相应的基准参数值输入趋势预警模型,根据趋势预警模型并应用综合预警分析策略对各监测设备进行动态预警分析,并对经综合预警分析策略判定为故障设备的监测设备启动预警。
7.根据权利要求6所述的基于动态仿真模型的核电厂设备故障预警系统,其特征在于:
该系统还包括:
滤波模块,用于该动态预警分析平台执行监测设备预警分析之前,对该实时数据采集中心获取的监测设备的实时运行参数值进行滤波处理。
8.根据权利要求7所述的基于动态仿真模型的核电厂设备故障预警系统,其特征在于:
该动态预警分析平台包括:
预警分析模块,用于应用残差预警分析、曲面面积积分预警分析及波动值预警分析相结合的综合预警分析策略对监测设备运行异常与否进行综合判定,并对经综合预警分析策略判定为故障设备的监测设备启动预警机制;
预警结果显示模块,用于将反映监测设备是否出现故障及故障严重程度的若干预警因子以动态的图形化的形式在趋势预警模型中进行实时显示,该预警因子包括该监测设备相对于相应参考基准的残差、偏移率及波动率。
9.根据权利要求8所述的基于动态仿真模型的核电厂设备故障预警系统,其特征在于:
所述预警分析模块包括:
残差预警分析子模块,其进一步包括残差值计算单元及第一预警分析单元;
其中:该残差值计算单元用于针对待执行预警分析的每台监测设备,将所述实时数据采集中心在设定的采样时间范围内采集的该台监测设备的实时运行参数值与经动态仿真模型模拟计算得到的该台监测设备的基准参数值进行比较,计算该实时运行参数值与该基准参数值之间的残差值;
该第一预警分析单元用于判断该残差值中的最大值是否超过设定的残差阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
10.根据权利要求9所述的基于动态仿真模型的核电厂设备故障预警系统,其特征在于:
所述预警分析模块还包括:
曲面面积积分预警分析子模块,其进一步包括偏移率计算单元及第二预警分析单元;
其中,该偏移率计算单元用于针对待执行预警分析的每台监测设备,在趋势预警模型中分别绘制表征该台监测设备的理想/健康运行状态及实际运行状态的基准曲线及实时运行曲线,将该条基准曲线作为参考线,并在设定的采样时间范围内将该台监测设备的实时运行曲线以该参考线为基准作曲线面积积分运算,将曲面面积换算成直角三角形面积,根据该直角三角形面积分析该段采样时间内该设备运行相对于该参考线的偏移率;
该第二预警分析单元用于判断该偏移值是否超过设定的偏移阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
11.根据权利要求9所述的基于动态仿真模型的核电厂设备故障预警系统,其特征在于:
所述预警分析模块还包括:波动值预警分析子模块,其进一步包括波动率计算单元及第三预警分析单元;
其中:该波动率计算单元用于针对待执行预警分析的核电站设备,在趋势预警模型中绘制表征该台监测设备处于理想/健康运行状态的基准曲线,将该条基准曲线作为参考线,利用均方差原理计算该台监测设备在指定的采样时间范围内相对于该参考线的波动率;
该第三预警分析单元用于判断该台监测设备实际运行过程中相对于该参考线的波动率是否超过设定的波动率阈值;如是,则做出该台监测设备可能出现故障的初步预警分析结果,并请求通过综合预警分析策略中包含的另外一种或两种预警分析方法对该预警分析结果作进一步核实。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510872454.2A CN105512812A (zh) | 2015-12-02 | 2015-12-02 | 基于动态仿真模型的核电厂设备故障预警分析方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510872454.2A CN105512812A (zh) | 2015-12-02 | 2015-12-02 | 基于动态仿真模型的核电厂设备故障预警分析方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105512812A true CN105512812A (zh) | 2016-04-20 |
Family
ID=55720774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510872454.2A Pending CN105512812A (zh) | 2015-12-02 | 2015-12-02 | 基于动态仿真模型的核电厂设备故障预警分析方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105512812A (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107204208A (zh) * | 2017-04-26 | 2017-09-26 | 岭东核电有限公司 | 一种便携式的中子探测器信号动态仿真装置 |
CN107369121A (zh) * | 2017-07-24 | 2017-11-21 | 中煤电气有限公司 | 一种煤矿信息化系统的实时数据模拟方法 |
CN107844619A (zh) * | 2017-06-21 | 2018-03-27 | 国家电网公司 | 500kV充油式电力设备运行数据分析与故障预警平台及方法 |
CN108052019A (zh) * | 2017-12-15 | 2018-05-18 | 中广核工程有限公司 | 核电站全范围模拟机流体模型计算溢出的抑制方法及系统 |
CN108172314A (zh) * | 2017-11-23 | 2018-06-15 | 中国核电工程有限公司 | 数字化压水堆核电厂主控室不可用故障处理策略分析方法 |
CN108490919A (zh) * | 2018-04-17 | 2018-09-04 | 安徽华电六安电厂有限公司 | 基于在线状态监测的分散控制系统设备可靠性评估方法 |
CN108875964A (zh) * | 2017-05-12 | 2018-11-23 | 深圳烟草工业有限责任公司 | 一种卷烟制造设备点检分析方法及系统 |
CN110175935A (zh) * | 2019-05-30 | 2019-08-27 | 中广核核电运营有限公司 | 核电数据监测方法、装置、计算机设备和存储介质 |
CN110658798A (zh) * | 2018-06-29 | 2020-01-07 | 株洲中车时代电气股份有限公司 | 轨道交通车辆传动控制单元的同步测试系统及方法 |
CN110866616A (zh) * | 2019-11-01 | 2020-03-06 | 许继集团有限公司 | 一种变电站二次设备故障预警方法及装置 |
CN111105886A (zh) * | 2019-12-18 | 2020-05-05 | 岭东核电有限公司 | 强辐射环境下的免校准温度测量方法及装置 |
CN111400850A (zh) * | 2018-12-28 | 2020-07-10 | 阿里巴巴集团控股有限公司 | 设备故障分析方法、装置、设备和存储介质 |
CN111524336A (zh) * | 2020-04-01 | 2020-08-11 | 广州盛康动力设备有限公司 | 一种发电机组预警方法及系统 |
CN111597637A (zh) * | 2020-05-19 | 2020-08-28 | 北京铁道工程机电技术研究所股份有限公司 | 一种动车现场部件的故障检测方法、装置、服务器及系统 |
CN111816338A (zh) * | 2020-06-08 | 2020-10-23 | 核动力运行研究所 | 一种核电厂信息系统健康监测与故障定位系统及方法 |
CN112069606A (zh) * | 2019-05-22 | 2020-12-11 | 赵英田 | 用于制氢储氢加氢站点的实时安全监测系统及监测方法 |
CN112447308A (zh) * | 2019-08-30 | 2021-03-05 | 海南核电有限公司 | 避免机组停机停堆的预警系统 |
CN113219939A (zh) * | 2021-04-07 | 2021-08-06 | 山东润一智能科技有限公司 | 基于残差自回归的设备故障预测方法及系统 |
CN113238530A (zh) * | 2021-04-25 | 2021-08-10 | 中国核电工程有限公司 | 监测系统显示功能设计方法及系统、验证方法、显示方法 |
WO2021174980A1 (zh) * | 2020-03-04 | 2021-09-10 | 广州明珞装备股份有限公司 | 一种维修任务分配方法、系统、装置及存储介质 |
CN113420954A (zh) * | 2021-05-08 | 2021-09-21 | 中国电建集团华东勘测设计研究院有限公司 | 基于机理模型的工程管理业务信息化方法 |
CN113554526A (zh) * | 2020-11-30 | 2021-10-26 | 国网北京市电力公司 | 电力设备的故障预警方法、装置、存储介质及处理器 |
CN113991855A (zh) * | 2021-10-29 | 2022-01-28 | 国网上海市电力公司 | 综合能源系统运行初期的性能监测及故障预警方法、系统 |
CN114881335A (zh) * | 2022-05-12 | 2022-08-09 | 华能澜沧江水电股份有限公司 | 基于大数据的水电机组运行数据趋势分析预警方法及系统 |
CN114964367A (zh) * | 2022-05-20 | 2022-08-30 | 中国船级社武汉规范研究所 | 一种基于时变参数的lng罐箱故障预报方法及系统 |
CN118118319A (zh) * | 2024-04-30 | 2024-05-31 | 南京风船云聚信息技术有限公司 | 一种基于大数据的网络设备智能诊断方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003042792A3 (en) * | 2001-11-14 | 2003-07-31 | Sap Ag | Inventory early warning agent with correction by error correlation calculation |
CN102705078A (zh) * | 2012-04-19 | 2012-10-03 | 哈尔滨工程大学 | 基于灰色模型的柴油机故障预测方法 |
CN103196682A (zh) * | 2013-03-21 | 2013-07-10 | 北京交通大学 | 基于d-s证据理论的信息融合的列车悬挂系统故障分离方法 |
CN104766139A (zh) * | 2015-03-27 | 2015-07-08 | 大唐淮南洛河发电厂 | 基于工业互联网的火电厂设备故障诊断与检修优化方法及系统 |
-
2015
- 2015-12-02 CN CN201510872454.2A patent/CN105512812A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003042792A3 (en) * | 2001-11-14 | 2003-07-31 | Sap Ag | Inventory early warning agent with correction by error correlation calculation |
CN102705078A (zh) * | 2012-04-19 | 2012-10-03 | 哈尔滨工程大学 | 基于灰色模型的柴油机故障预测方法 |
CN103196682A (zh) * | 2013-03-21 | 2013-07-10 | 北京交通大学 | 基于d-s证据理论的信息融合的列车悬挂系统故障分离方法 |
CN104766139A (zh) * | 2015-03-27 | 2015-07-08 | 大唐淮南洛河发电厂 | 基于工业互联网的火电厂设备故障诊断与检修优化方法及系统 |
Non-Patent Citations (2)
Title |
---|
朱纪洪 等: "舵机特征模型及其故障检测方法", 《航空学报》 * |
郭鹏 等: "风电机组齿轮箱温度趋势状态检测及分析方法", 《中国电机工程学报》 * |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107204208A (zh) * | 2017-04-26 | 2017-09-26 | 岭东核电有限公司 | 一种便携式的中子探测器信号动态仿真装置 |
CN108875964A (zh) * | 2017-05-12 | 2018-11-23 | 深圳烟草工业有限责任公司 | 一种卷烟制造设备点检分析方法及系统 |
CN107844619A (zh) * | 2017-06-21 | 2018-03-27 | 国家电网公司 | 500kV充油式电力设备运行数据分析与故障预警平台及方法 |
CN107369121A (zh) * | 2017-07-24 | 2017-11-21 | 中煤电气有限公司 | 一种煤矿信息化系统的实时数据模拟方法 |
CN108172314A (zh) * | 2017-11-23 | 2018-06-15 | 中国核电工程有限公司 | 数字化压水堆核电厂主控室不可用故障处理策略分析方法 |
CN108052019A (zh) * | 2017-12-15 | 2018-05-18 | 中广核工程有限公司 | 核电站全范围模拟机流体模型计算溢出的抑制方法及系统 |
CN108052019B (zh) * | 2017-12-15 | 2021-01-19 | 中广核工程有限公司 | 核电站全范围模拟机流体模型计算溢出的抑制方法及系统 |
CN108490919A (zh) * | 2018-04-17 | 2018-09-04 | 安徽华电六安电厂有限公司 | 基于在线状态监测的分散控制系统设备可靠性评估方法 |
CN110658798A (zh) * | 2018-06-29 | 2020-01-07 | 株洲中车时代电气股份有限公司 | 轨道交通车辆传动控制单元的同步测试系统及方法 |
CN111400850A (zh) * | 2018-12-28 | 2020-07-10 | 阿里巴巴集团控股有限公司 | 设备故障分析方法、装置、设备和存储介质 |
CN111400850B (zh) * | 2018-12-28 | 2024-04-09 | 阿里巴巴集团控股有限公司 | 设备故障分析方法、装置、设备和存储介质 |
CN112069606A (zh) * | 2019-05-22 | 2020-12-11 | 赵英田 | 用于制氢储氢加氢站点的实时安全监测系统及监测方法 |
CN112069606B (zh) * | 2019-05-22 | 2022-04-05 | 赵英田 | 用于制氢储氢加氢站点的实时安全监测系统及监测方法 |
CN110175935A (zh) * | 2019-05-30 | 2019-08-27 | 中广核核电运营有限公司 | 核电数据监测方法、装置、计算机设备和存储介质 |
CN112447308A (zh) * | 2019-08-30 | 2021-03-05 | 海南核电有限公司 | 避免机组停机停堆的预警系统 |
CN110866616A (zh) * | 2019-11-01 | 2020-03-06 | 许继集团有限公司 | 一种变电站二次设备故障预警方法及装置 |
CN111105886A (zh) * | 2019-12-18 | 2020-05-05 | 岭东核电有限公司 | 强辐射环境下的免校准温度测量方法及装置 |
WO2021174980A1 (zh) * | 2020-03-04 | 2021-09-10 | 广州明珞装备股份有限公司 | 一种维修任务分配方法、系统、装置及存储介质 |
CN111524336A (zh) * | 2020-04-01 | 2020-08-11 | 广州盛康动力设备有限公司 | 一种发电机组预警方法及系统 |
CN111597637A (zh) * | 2020-05-19 | 2020-08-28 | 北京铁道工程机电技术研究所股份有限公司 | 一种动车现场部件的故障检测方法、装置、服务器及系统 |
CN111816338A (zh) * | 2020-06-08 | 2020-10-23 | 核动力运行研究所 | 一种核电厂信息系统健康监测与故障定位系统及方法 |
CN113554526A (zh) * | 2020-11-30 | 2021-10-26 | 国网北京市电力公司 | 电力设备的故障预警方法、装置、存储介质及处理器 |
CN113219939A (zh) * | 2021-04-07 | 2021-08-06 | 山东润一智能科技有限公司 | 基于残差自回归的设备故障预测方法及系统 |
CN113238530A (zh) * | 2021-04-25 | 2021-08-10 | 中国核电工程有限公司 | 监测系统显示功能设计方法及系统、验证方法、显示方法 |
CN113420954A (zh) * | 2021-05-08 | 2021-09-21 | 中国电建集团华东勘测设计研究院有限公司 | 基于机理模型的工程管理业务信息化方法 |
CN113991855A (zh) * | 2021-10-29 | 2022-01-28 | 国网上海市电力公司 | 综合能源系统运行初期的性能监测及故障预警方法、系统 |
CN114881335A (zh) * | 2022-05-12 | 2022-08-09 | 华能澜沧江水电股份有限公司 | 基于大数据的水电机组运行数据趋势分析预警方法及系统 |
CN114881335B (zh) * | 2022-05-12 | 2024-06-28 | 华能澜沧江水电股份有限公司 | 基于大数据的水电机组运行数据趋势分析预警方法及系统 |
CN114964367A (zh) * | 2022-05-20 | 2022-08-30 | 中国船级社武汉规范研究所 | 一种基于时变参数的lng罐箱故障预报方法及系统 |
CN114964367B (zh) * | 2022-05-20 | 2024-02-23 | 中国船级社武汉规范研究所 | 一种基于时变参数的lng罐箱故障预报方法及系统 |
CN118118319A (zh) * | 2024-04-30 | 2024-05-31 | 南京风船云聚信息技术有限公司 | 一种基于大数据的网络设备智能诊断方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105512812A (zh) | 基于动态仿真模型的核电厂设备故障预警分析方法及系统 | |
CN105551549B (zh) | 一种核电设备运行状况在线监测方法及系统 | |
US8099672B2 (en) | System and method for continuous online safety and reliability monitoring | |
CN108803569A (zh) | 电站锅炉诊断专家系统及其故障诊断方法 | |
CN104298225A (zh) | 化工过程异常工况因果关系推理模型建模与图形化展示方法 | |
KR20120046821A (ko) | 범용 센서 자가 진단 장치 및 그 진단 방법 | |
WO2015007040A1 (zh) | 核电站运行状态监控方法、装置及系统 | |
CN113420465B (zh) | 一种基于数字孪生模型的液压支架全寿命周期管理方法 | |
KR101853480B1 (ko) | 해양플랜트 예지보전 시스템 | |
CN106228472B (zh) | 一种变电站设备智能巡检故障分析方法 | |
CN110688774A (zh) | 一种电网故障仿真推演系统及方法 | |
CN103760895A (zh) | 新能源汽车车载充电机控制装置硬件在环测试设备及方法 | |
CN114239864A (zh) | 一种基于数字孪生的换流阀故障检修方法及终端 | |
WO2023240902A1 (zh) | 核电机组风险预测方法、预测系统及评估系统 | |
CN102184130A (zh) | 一种芯片流水线的测试方法和装置 | |
CN108051637A (zh) | 一种智能电能表时钟电池故障诊断方法 | |
CN115076452A (zh) | 一种基于数字孪生的高压隔膜泵单向阀健康管理系统及方法 | |
CN110210722A (zh) | 一种核电厂系统监督方案的开发方法 | |
CN117038048B (zh) | 一种用于医疗器械的远程故障处理方法及系统 | |
CN117375246A (zh) | 电力设备管理方法、系统及设备 | |
CN112363432A (zh) | 一种水电站辅助设备的监测系统及监测方法 | |
CN116859870A (zh) | 一种火电厂机组运行装置 | |
CN113238530B (zh) | 监测系统显示功能设计方法及系统、验证方法、显示方法 | |
CN114065990A (zh) | 一种火电厂中给水泵的监测方法、存储介质和电子装置 | |
CN113386976B (zh) | 一种大型飞机燃油系统全模试验方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160420 |