KR20170045354A - 희소 코딩 방법론을 활용한 가스 터빈 센서 고장 검출 - Google Patents

희소 코딩 방법론을 활용한 가스 터빈 센서 고장 검출 Download PDF

Info

Publication number
KR20170045354A
KR20170045354A KR1020177009440A KR20177009440A KR20170045354A KR 20170045354 A KR20170045354 A KR 20170045354A KR 1020177009440 A KR1020177009440 A KR 1020177009440A KR 20177009440 A KR20177009440 A KR 20177009440A KR 20170045354 A KR20170045354 A KR 20170045354A
Authority
KR
South Korea
Prior art keywords
sensor
readings
gas turbine
data
gambling
Prior art date
Application number
KR1020177009440A
Other languages
English (en)
Other versions
KR101955305B1 (ko
Inventor
시옹 테 고
차오 유안
아미트 챠크라보르티
매튜 에반스
Original Assignee
지멘스 에너지, 인코포레이티드
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 지멘스 에너지, 인코포레이티드 filed Critical 지멘스 에너지, 인코포레이티드
Publication of KR20170045354A publication Critical patent/KR20170045354A/ko
Application granted granted Critical
Publication of KR101955305B1 publication Critical patent/KR101955305B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D3/00Indicating or recording apparatus with provision for the special purposes referred to in the subgroups
    • G01D3/08Indicating or recording apparatus with provision for the special purposes referred to in the subgroups with provision for safeguarding the apparatus, e.g. against abnormal operation, against breakdown
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/14Testing gas-turbine engines or jet-propulsion engines
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0224Process history based detection method, e.g. whereby history implies the availability of large amounts of data
    • G05B23/024Quantitative history assessment, e.g. mathematical relationships between available data; Functions therefor; Principal component analysis [PCA]; Partial least square [PLS]; Statistical classifiers, e.g. Bayesian networks, linear regression or correlation analysis; Neural networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Automation & Control Theory (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Biology (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Combustion & Propulsion (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Testing And Monitoring For Control Systems (AREA)
  • Testing Of Engines (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

가스 터빈들을 모니터링하는 데 사용되는 센서들의 고장들을 인식(및/또는 예측)하기 위한 방법 및 시스템은, 희소 코딩 프로세스를 수집된 센서 판독들에 적용하고 그리고 희소 코딩 프로세스로부터의 L-1 노름 잔차들이 잠재적 센서 문제를 표시하는 것으로 정의한다. 잔차 센서 판독들의 그룹의 추가의 평가는 그룹을 카테고리화하기 위해 그리고 잡음있는 데이터보다는 결함있는 센서와 더 연관될 가능성이 있는 것으로 고려될 상당한 이상점들(비정상 데이터)이 존재하는지를 결정하기 위해 수행된다. 현재의 비정상 결과를 이전 결과들의 세트와 비교하는 평가에 시간 컴포넌트가 도입되며, 상당한 수의 이전 판독들이 또한 비정상 값을 갖는 경우 결함있는 센서 결정을 한다. 시간 컴포넌트를 고려함으로써, 허위 포지티브들의 수가 감소된다.

Description

희소 코딩 방법론을 활용한 가스 터빈 센서 고장 검출{GAS TURBINE SENSOR FAILURE DETECTION UTILIZING A SPARSE CODING METHODOLOGY}
관련 출원들에 대한 상호-참조
본 출원은, 2014년 9월 10일 출원된 미국 가출원 번호 제 62/048,588호를 우선권으로 주장하며, 이 미국 가출원은 인용에 의해 본원에 포함된다.
기술 분야
본 발명은, 전력 생성 네트워크(power generation network)들에서 활용되는 가스 터빈(gas turbine)들의 관리감독(supervision), 및 특정 센서(sensor)의 고장(failure)을 예측하는 능력을 포함하여 터빈의 성능을 모니터링(monitor)하기 위해 사용되는 센서들 중 임의의 센서의 고장을 검출하는 능력을 개선하기 위한 방법론에 관한 것이다.
화석 연료들의 효율적인 사용은 안정적인 전력 네트워크를 유지하는 데 중요하다. 이러한 유형의 연료를 전기 에너지(electrical energy)로 변환하는 특히 효율적인 수단은 가스 터빈이다. 가스 터빈 컴포넌트(component)들은 매우 고온의 환경에서 그리고 다양한 부하 상태들 하에서 동작한다. 열 피로 및 마모로 인한 부분들의 열화는 현실적인 관심사이다. 효율적인 동작을 계속해서 보장하기 위하여 마모를 검출 및 제어하기 위해서뿐만 아니라 필요에 따라서는 마모된 부분들을 수리 또는 교체하기 위해서 유지보수가 수행된다.
가스 터빈의 성능은 통상적으로, 가스 터빈의 동작의 다양한 양상들을 평가하기 위해 다양한 상이한 센서들(즉, 전력 센서들, 온도 센서들, 압력 센서들 등)을 사용함으로써 모니터링된다(monitored). 불행하게도, 센서 판독(sensor reading)들 그 자체는 상대적으로 잡음이 있는 경향이 있어서, 때때로 센서가 적절하게 동작하는지를 알기 어렵다.
결함있는 센서(faulty sensor)들을 검출하는 문제를 연구하기 위해 사용된 다양한 접근법들이 종래기술에 존재한다. 예컨대, 일 접근법은 센서 판독들에서 3개의 특정 유형들의 에러(error)들, 즉, 짧은 결함(short fault)들, 지속적인 결함(constant fault)들 및 잡음 결함(noise fault)들을 모니터링한다. 센서 판독들에서 "결함있는 센서" 시그니처(signature)들을 탐색하기 위해, 주요 컴포넌트 분석(PCA; principal component analysis)을 활용하는 전통적인 통계적 방법론들이 수행되었다.
가스 터빈 센서들의 결함 검출을 수행하기 위한 이들 접근법들 및 다른 접근법들은 임계치-기반 접근법을 사용한다. 즉, 특정 센서 판독이 주어진 임계치를 초과하는 경우, 센서는 "결함있음(faulty)"으로 선언된다. 작업가능한 동안, 이 접근법은 센서의 작업 환경의 시간 컴포넌트를 고려하지 않는다. 부가적으로, 일부 센서 판독들에는 큰 잡음 컴포넌트가 존재할 수 있기 때문에, 임계치 접근법은 잡음있는 데이터(noisy data)와 결함있는 센서 간의 차이를 판별하는 데 어려움을 겪을 수 있다. 결과적으로, 이 임계치-기반 접근법은 많은 수의 허위 포지티브(false positive)들을 산출할 수 있는데, 이들은 센서들의 불필요한 검사들로 옮겨져서, 가스 터빈의 성능을 방해할 수 있다.
종래 기술의 남아있는 필요성들은, 전력 생성 네트워크들에서 활용되는 가스 터빈들의 관리감독 및 특정 센서의 고장을 예측하는 능력을 포함하여 센서 고장을 정확하게 검출하는 능력을 개선하기 위한 방법론에 관한 본 발명에 의해 처리된다.
본 발명에 따르면, 센서 판독들을 분석하기 위해 희소 코딩 프로세스(sparse coding process)가 사용된다. 희소 코딩 접근법은 이 태스크(task)에 아주 적합한 것으로 고려되는데, 그 이유는 정상 센서 판독(normal sensor reading)들의 방대한 집합 가운데서 최소 개수의 "결함있는 센서" 판독들을 식별하는 것이 그 목적이기 때문이다. 희소 코딩 프로세스로부터의 잔차(residual)들은 잠재적으로 비정상 센서 판독(abnormal sensor reading)들로서 정의되며, 시간 팩터(time factor)를 분석에 도입하기 위해 그리고 잔차를 "잡음있는 데이터" 또는 "결함있는 센서"로서 정의하기 위한 결정을 하기 위해 추가로 프로세싱된다(processed).
본 발명의 프로세스의 희소 코딩 부분은 2개의 메인 컴포넌트(main component)들, 즉, (1) 딕셔너리 학습(dictionary learning) 및 (2) 딕셔너리 학습에 의해 정의된 좌표 공간 내의 수집된 데이터의 표현들을 컴퓨팅(computing)하고 그리고 딕셔너리의 벡터 공간(vector space) 내에 피팅(fit)되지 않는 데이터를 비정상성(anomaly)들로서 정의하는 것에 의한 비정상성 검출(anomaly detection)로 이루어진다. 본 발명에 따르면, 비정상성들(잔차들)은 정상 동작 범위를 초과하는 센서 판독들의 형태를 취한다. 본 발명의 바람직한 실시예에서, L-1 노름(L-1 norm)은 (종래의 L-2 노름 대신에) 잔차를 정의하기 위해 사용되는데, 그 이유는 L-1 노름("최소 절대 편차(least absolute deviation)" 또는 "최소 절대 에러(least absolute error)" 값으로 또한 알려져 있음)이 허위 포지티브 결과들을 산출할 가능성이 더 낮기 때문인데, 허위 포지티브 결과들은, 통상적으로 잡음있는 가스 터빈 센서 판독들을 연구할 때 문제가 된다.
희소 코딩에 의해 식별된 잔차들을 평가하기 위해 그리고 잔차들이 ― 결함있는 센서가 존재하는 경우 ― 결함있는 센서를 표시하는지를 결정하기 위해, 사후-프로세싱 기능(post-processing function)들이 사용된다. 일 예시적 사후-프로세싱 접근법은, 잔차들의 세트(set)의 모델(model)을 생성하기 위해 일-클래스(one-class) 지원 벡터 머신(SVM; Support Vector Machine)을 이용한다. 다른 사후-프로세싱 접근법은 잔차들의 그룹(group)의 특정 값들에 기반하여 임계값을 개발한다. 어느 경우이든, 시간-기반 판독들의 세트는 사용되고 있는 모델에 대해 비교되며, 시간에 따라 수집된 충분한 양의 데이터가 기준들을 충족하는 경우, "결함있는 센서" 검출이 선언된다. 그 후에, 수리 또는 교체되어야 할 특정 센서를 식별하는 "센서 고장" 메시지(message)가 가스 터빈 유지보수를 담당하는 직원에게 전송된다.
일 실시예에서, 본 발명은 적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법의 형태를 취하며, 방법은, (1) 정상 동작 상태들에 대해 알려진 센서 판독들과 연관된 값들을 정의하는 기저 벡터(basis vector)들의 딕셔너리를 생성하는 단계; (2) L-1 노름 잔차 데이터(L-1 norm residual data)의 세트를 식별하기 위해, 생성된 딕셔너리를 사용하여, 희소 코딩 프로세스를 센서 판독 데이터의 세트에 적용하는 단계; (3) 최대-값의 L-1 노름 잔차 데이터(largest-valued L-1 norm residual data)의 미리 결정된 서브세트(subset)를 비정상 센서 판독(abnormal sensor reading)들로 카테고리화(categorize)하기 위해 L-1 노름 잔차 데이터를 평가하는 단계; (4) 비정상 센서 판독들을 복수의 이전(prior) 센서 판독들과 비교하고 그리고 이전 센서 판독들 중 미리 정의된 개수의 이전 센서 판독들이 또한 비정상 센서 판독들로서 지정되는 경우, 지정된 비정상 센서 판독들을 센서 고장과 연관된 것으로 정의하는 단계; (5) 특정 센서를 식별하는 센서 고장 신호를 가스 터빈 직원(gas turbine personnel)에게 송신하는 단계; 및 (6) 고장난 센서를 서비스(service)로부터 제거하고 그리고 고장난 센서를 수리 또는 교체하는 단계를 포함한다.
다른 실시예에서, 본 발명은 가스 터빈 센서의 고장을 검출하기 위한 시스템(system)으로서 설명되며, 시스템은 이러한 컴포넌트들, 즉, (1) 센서 판독들의 데이터베이스(database); (2) 센서 판독들의 데이터베이스와 통신하는 센서 모니터링 시스템 컴포넌트(sensor monitoring system component) ― 센서 모니터링 시스템 컴포넌트는 프로그램 저장 디바이스(program storage device) 및 프로세서(processor)를 포함하고, 프로그램 저장 디바이스는 방법 단계들을 수행하기 위해 프로세서에 의해 실행가능한 프로그램 명령들의 세트를 고정형의 유형적 매체(fixed tangible medium)에 포함하고, 방법 단계들은, (a) 정상 동작 상태들에 대해 알려진 센서 판독들과 연관된 값들을 정의하는 기저 벡터들의 딕셔너리를 생성하는 단계; (b) L-1 노름 잔차 데이터의 세트를 식별하기 위해, 생성된 딕셔너리를 사용하여, 희소 코딩 프로세스를 센서 판독 데이터의 세트에 적용하는 단계; (c) 최대-값의 L-1 노름 잔차 데이터의 미리 결정된 서브세트를 비정상 센서 판독들로 카테고리화하기 위해 L-1 노름 잔차 데이터를 평가하는 단계; (d) 비정상 센서 판독들을 복수의 이전 센서 판독들과 비교하고 그리고 이전 센서 판독들 중 미리 정의된 개수의 이전 센서 판독들이 또한 비정상 센서 판독들로서 지정되는 경우, 지정된 비정상 센서 판독들을 센서 고장과 연관된 것으로 정의하는 단계 및 (e) 수리 또는 교체되도록 하기 위해, 고장난 센서를 식별하는 정보를 포함하는 센서 고장 신호를 가스 터빈 직원에게 송신하는 단계임 ―; 및 (3) 비정상(abnormal)으로서 분류된 센서 판독들을 저장하기 위한 모니터링 시스템 데이터베이스를 포함한다.
본 발명의 또 다른 실시예에서, 희소 코딩 프로세스에 의해 평가된 센서 판독들은 (고장이 이미 발생한 상태와는 대조적으로) 임박한 센서 고장을 예측하기 위해 사용될 수 있다.
본 발명의 다른 양상들 및 실시예들 그리고 추가의 양상들 및 실시예들은 다음의 논의의 과정 동안 그리고 첨부 도면들을 참조하여 명백해질 것이다.
이제 도면들을 참조하면, 여러 도면들에서 동일한 숫자들은 동일한 부분들을 표현한다:
도 1은 결함 센서들을 검출하는 본 발명의 시스템 및 방법이 사용될 수 있는 예시적 가스 터빈 환경의 간략화된 도면이고;
도 2는 수집된 센서 판독들을 평가하기 위해 그리고 하나 또는 그 초과의 센서들이 결함있는지 및 수리 또는 교체될 필요가 있는지를 결정하기 위해 활용될 수 있는 본 발명의 예시적 프로세스의 흐름도이고;
도 3은, 센서 고장들을 검출하기 위해 L-2 노름(L-2 norm)을 사용하는 것과 연관된 문제들을 예시하는, L-2 노름을 희소 코딩 프로세스에 적용함으로써 컴퓨팅된(computed) 잔차들의 다이어그램(diagram)이고;
도 4는 L-1 노름을 희소 코딩 프로세스에 적용함으로써 본 발명의 교시들에 따라 컴퓨팅된 잔차들의 다이어그램이고;
도 5는, 어느 잔차들이 적절하게 "비정상"인지를 결정하고 센서의 고장을 표시하기 위해, 일-클래스 지원 벡터 머신(SVM; Support Vector Machine) 모델 접근법을 사용하여 L-1 노름 잔차들을 평가하기 위한 방법의 흐름도이고; 그리고
도 6은 L-1 노름 잔차들을 평가하고 센서 고장의 존재를 검출하기 위한 대안적 방법의 흐름도이다.
도 1은, 플랜트 전기 부하(plant electric load)(3)를 공급하는 발생기(2)를 갖는 통상의 가스 터빈 전력 플랜트(1)의 간략화된 도면이다. 발생기(2)는 가스 터빈 엔진(gas turbine engine)(5)에 의해 동력이 공급되는 축(4)에 의해 구동된다. 가스 터빈 엔진(5) 그 자체는, 압축기(5.1), 연소 섹션(combustion section)(5.2), 터빈(5.3) 및 아마도 조정가능한 입구 베인(inlet vane)들(5.4)의 세트를 포함한 많은 수의 개별 컴포넌트들로 이루어진다.
연료는 밸브(valve)(6)를 통해 연소 섹션(5.2)에 공급된다. 가스 터빈 전력 플랜트(1)의 수용가능한 동작을 유지하기 위해, 다수의 센서들(7)이 사용되어 다양한 컴포넌트들의 동작을 모니터링(monitor)하고, 측정된 센서 판독들을 개별 제어 모듈(control module)(8)에 전달한다. 제어 모듈(8)은 가스 터빈 전력 플랜트(1)와 병치될(co-located) 수 있거나 또는 터빈 그 자체로부터 떨어져 있을 수 있다. 도 1의 다이어그램에서, 센서들(7)은 연소기 유입 공기 센서(combustor inlet air sensor)(7.1), 연소 온도 센서(7.2), 및 블레이드 경로 온도 센서(blade path temperature sensor)(7.3)를 포함한다. 가스 터빈의 성능을 모니터링하기 위해 사용되는 더 많은 센서들이 존재하여서, 온도, 압력, 회전, 진동 등과 같은 상태들을 측정한다는 것이 이해되어야 한다. 실제로, 주어진 가스 터빈 전력 플랜트에 따라 200개에 가까운 상이한 센서들이 활용될 수 있는 것이 가능하다.
제어 모듈(8)은 센서들(7)로부터 입력을 수신하고, 제어 신호들을 당해 기술 분야에 알려진 바와 같은 밸브들, 모터(motor)들, 및 액추에이터(actuator)들에 송신한다. 제어기는, 당해 기술 분야에 알려진 바와 같은 연관된 하드웨어(hardware) 및 소프트웨어(software)를 갖는 하나 또는 그 초과의 프로세서들을 하나 또는 그 초과의 위치들에 포함할 수 있다.
가스 터빈의 성능이 정확한 센서 판독들에 종속되는 것을 고려하면, 센서가 고장나기 시작하여 부정확한 판독들을 전달할 때를 인지할 수 있는 것이 중요하다. 위에서 언급된 바와 같이, 센서 거동을 이해하는 중요한 양상은, 일회성의 결함있는 판독(일회성의 결함있는 판독은 다양한 환경 팩터들로 인한 것일 수 있고 따라서 무시될 수 있음)과 일련의 계속되는 결함있는 판독들(일련의 계속되는 결함있는 판독들은 프리젠테이션(presentation)에서 더 미묘할 수 있음) 간을 판별하는 능력을 포함하는 것인데, 그 이유는 일련의 계속되는 결함있는 판독들의 상태가, 고장나서 교체될 필요가 있는 센서로부터 초래될 가능성이 더 많기 때문이다. 더욱이, 위에서 언급된 바와 같이, 가스 터빈 센서들에 의해 수집된 통상의 아날로그 센서 판독(analog sensor reading)들은 상대적으로 큰 잡음 컴포넌트를 포함하는 경향이 있어서, "잡음있는 데이터"와 "결함있는 센서" 간을 구별하는 능력을 제한한다.
기본 임계치 센서 판독을 "알람(alarm)"으로서 사용하는 다양한 종래 기술의 기법들과 대조적으로, 본 발명의 방법론은 센서 판독들의 시계열 데이터 세트(time series data set)에 대한 분석을 수행한다. 수용가능한 센서 판독들의 유니버스(universe)를 적절하게 특징짓는 딕셔너리를 유도하기 위해 희소 코딩 기법이 사용되며, 여기서 기저 벡터들의 정의된 경계들 외측에 있는 임의의 판독들("잔차(residual)"로 지칭됨)은 잠재적 문제로 분류된다. 상황을 "결함있는 센서"로서 정의하기 위한 시간의 함수로써 충분한 양의 비정상 판독들이 존재하는지를 결정하기 위해, 잔차들의 모델이 생성되어 사용된다.
그러므로, 본 발명에 따르면, 제어 모듈(8)의 성능들은 센서 모니터링 시스템(10)을 포함하도록 확장되며, 센서 모니터링 시스템(10)은 아래에서 상세하게 논의되는 바와 같이, 수집된 센서 판독들을 분석하고 센서 판독들에서 특정 센서의 고장의 표시자인 임의의 텔-테일 변화(tell-tale change)들을 찾기 위해 활용될 수 있거나 또는 다음 수 일 내에 센서 고장이 발생할 가능성이 있다는 것을 합리적으로 예측할 수 있다. 아래에서 상세하게 논의될 바와 같이, 센서 모니터링 시스템(10)은 (수년 전으로 되돌아갈 수 있는) 모든 이력 센서 판독들의 데이터베이스(12)를 포함하는데; 실제로 이 데이터베이스는 종래의 제어 모듈(8)의 부분을 형성할 수 있다. 센서 모니터링 시스템(10)은 메모리(memory)(16) 및 프로세서(18)를 갖는 프로세서 컴포넌트(14)를 포함하며, 이들 엘리먼트(element)들은 센서 판독 데이터의 실제 평가를 수행하고 센서들 그 자체의 성능을 모니터링하는 데 사용된다. 그 후에, "비정상" 센서 판독의 경보의 형태를 취하는 센서 모니터링 시스템(10)으로부터의 출력은, 센서가 교체 또는 수리될 필요가 있는지의 여부를 결정하는 데 있어서 플랜트 직원에 의해 사용된다. 도 1에 도시된 바와 같은 센서 모니터링 시스템(10)의 특정 아키텍처(architecture)는 단지 예시적인 것으로 고려되며, 아래에서 상세하게 설명되는 바와 같은 본 발명의 방법을 수행하도록 다양한 다른 구성들의 하드웨어 및 소프트웨어 컴포넌트들이 구성될 수 있다.
도 2는 수집된 센서 판독들을 평가하고 하나 또는 그 초과의 센서들이 결함있는지 및 수리 또는 교체될 필요가 있는지를 결정하기 위해 센서 모니터링 시스템(10)에 의해 활용될 수 있는 본 발명의 예시적 프로세스의 흐름도이다.
프로세스는, 단계(200)에 도시된 바와 같이, 실제 센서 판독 데이터를 수집하고 데이터를, 데이터가 나중에 리트리빙되고(retrieved) 분석될 수 있는 엘리먼트에 저장함으로써 시작된다. 예컨대, 도 1의 제어기(8)는 이러한 기능들을 수행하기 위해 메모리 엘리먼트들 및 프로세싱 성능들을 포함하는 것으로 가정된다. 일단 데이터가 수집되면, 예비 단계는, 제거되지 않으면 분석의 결과들을 왜곡(skew)할 명백한 이상점(outlier)들 및 리던던트(redundant) 데이터를 제거하기 위해 센서 판독들을 "사전프로세싱(preprocessing)"하는 단계(단계(210))를 수반할 수 있다. 이 사전프로세싱 단계가 의무적인 것이라고 고려되지 않지만(그리고 흐름도에서 점선 박스(box)로서 도시됨), 본 발명의 바람직한 실시예에서, 사전프로세싱 단계는 효율적인 방식으로 정확한 해결책(solution)에 도달하기 위해 활용된다.
단계(220)는 수집된 데이터 세트에 대해 "희소 코딩"을 수행하는 실제 프로세스와 연관된다. 아래에서 더 상세하게 논의될 바와 같이, 희소 코딩은, 데이터 세트를 정확하게 표현하기 위해(이 경우, 센서 판독들을 정확하게 표현하기 위해) 사용될 수 있는 기저 벡터들의 세트를 정의하는 것을 포함한다. "잔차들"로 지칭되는, 희소 코딩 벡터 공간 내에 피팅(fit)되지 않은 소수의 데이터 포인트(few data point)들은, 잡음이 있을 수 있거나 또는 그렇지 않으면 결함있는 센서들과 연관될 수 있는 그러한 데이터이다. 데이터가 벡터 공간 외측에 있는 레벨(level)을 결정하기 위해 다양한 상이한 벡터 노름들이 사용될 수 있지만, 센서 판독들에 기인하는 잡음있는 데이터는 L-1 노름(즉, 최소 절대 편차 노름(least absolute deviation norm))을 사용하는 것에 제일 적합하다.
그 후에, 단계(220)에서 식별된 바와 같은 잔차들은, 데이터 포인트들이 단지 잡음있는 데이터인지 또는 정말로 결함있는 센서와 연관되었는지를 궁극적으로 결정하기 위해 추가로 평가된다. 다음 단계(도 2의 흐름도에서 단계(230)로 도시됨)는 잔차들을 2개의 그룹들, 즉, 트레이닝 세트(training set)로서 사용되는 제 1 그룹 및 테스트 세트(test set)로서 사용되는 제 2 그룹으로 분할한다. 그 후에, 트레이닝 잔차들은 아래에서 논의될 바와 같이, 잔차들의 그룹을 평가하고 그리고 소수의 가장 큰 잔차들을 센서의 "비정상" 동작과 연관된 것으로 정의하기 위한 한계치를 결정하는 데 적합한 모델을 개발하기 위해 사용된다(단계(240)).
일단 모델의 파라미터(parameter)들이 식별되면, 테스트 데이터에 대해 희소 코드 프로세싱(sparse code processing)을 수행함으로써 생성된 잔차들이 이 모델을 통과하며(단계(250)), 여기서 비교가 수행되고, 테스트되고(tested) 있는 현재의 "잔차"가 모델 내에 있는지 또는 모델 외측에 있는지의 여부의 결정이 이루어진다(단계(260)). 그 후에, 모델 외측에 있는 그러한 잔차들은 "비정상"으로 식별되고, 그와 같이 라벨링되고(labeled), 추가의 사용을 위해 메모리에 저장된다.
프로세스의 이 시점에서 분석의 시간 컴포넌트가 도입된다. "비정상" 데이터 포인트의 디스커버리(discovery)에 따라, 프로세스의 다음 단계(단계(270))는 다수의 이전의 시간 포인트들에 걸쳐 데이터 포인트들을 조사(review)하고 이들 데이터 포인트들 중 상당한 수가 또한 "비정상"으로서 식별되었는지를 확인하는 것이다. 또한, 이 시점에서, 주어진 시간 기간에 걸쳐 발견되도록 요구되는 "비정상" 판독들의 수를 평가하는 데 있어서 상이한 측정들이 사용될 수 있다.
시간의 함수로써 센서 판독 데이터의 상당한 퍼센티지(percentage)가 계속해서 "비정상"으로 플래그되는(flagged) 경우(단계(280)), 연관된 센서가 결함이 있을 수 있고 실제 디바이스를 평가할 필요가 있다는 결정이 이루어진다(단계(290)). 그렇지 않으면, 현재의 "비정상" 판독은 메모리 엘리먼트에 유지되고, 센서 판독들을 모니터링하는 프로세스가 계속된다(단계(300)). 센서가 고장나기 시작함에 따라 센서 판독 데이터가 저하되기 시작하여 더 많은 수의 잔차들을 생성할 수 있기 때문에, 도 2의 흐름도에서 간략히 서술된 바와 같은 프로세스는 또한, 센서 고장들의 발생 전에 센서 고장들을 예측하는 데 유용하다는 것이 이해되어야 한다.
다음의 논의는, 특히, 희소 코딩 기법이 본 발명의 가스 터빈/센서 환경에 적용됨에 따른 희소 코딩 기법의 간략한 개요를 포함하며, 그 후에 잔차들을 분석하고 시간 팩터를 분석에 도입하기 위한 모델을 개발하는 데 사용될 수 있는 2개의 상이한 사후-프로세싱 기법들의 논의가 뒤따른다. 그 후에, 실제 가스 터빈 센서 판독들에 대한 본 발명의 방법론의 적용의 예가 뒤따른다.
"희소 코딩" 프로세스는 일반적으로, 입력 신호의 간결한 표현을 발견하기 위한 알고리즘(algorithm)들의 클래스(class)와 관련된다(이 경우, 입력 신호는 가스 터빈 센서 판독들의 형태를 취함). 특히, 신호 분석 기법으로서의 희소 코딩은, 효율적인 방식으로 데이터를 표현하는 오버-컴플리트 베이스(over-complete base)들의 세트들(때때로 기저 벡터들의 세트로 또한 지칭됨)을 학습하기 위한 방법론의 부분이다. "희소(sparsity)"는, (소수의 큰 계수들을 제외한) 변환 계수들의 대부분이 제로(zero)로 근사화될 때 존재하는 제한된 수의 비-제로 컴포넌트들을 나타낸다. 본 발명에 따르면, 연구되고 있는 데이터는 시계열 센서 판독들(각각의 센서에 대해 수집된 데이터의 상이한 세트)이다. 가스 터빈 센서 데이터를 연구하기 위한 종래의 기법들(이를테면, 위에서 언급된 PCA 접근법)과는 대조적으로, 희소 코딩은 기저 벡터들이 입력 데이터에서 구조들 및 패턴(pattern)들(뿐만 아니라 비정상성들)을 더 양호하게 캡처(capture)하는 능력을 제공한다(희소 코딩이 기저 벡터들의 오버-컴플리트 세트이기 때문임). 희소 코딩은 주어진 입력 데이터 세트 X가 딕셔너리
Figure pct00001
에서 하나 또는 그 초과의 기저 벡터들의 선형 조합으로서 표현되도록 허용한다.
본 발명에 따르면, 결함있는 센서 판독들은 "비정상성들"로서 정의되며, 비정상성 검출과 함께 희소 코딩 프로세스를 활용함으로써 컴퓨팅될(computed) 수 있다. 희소 코딩은 2개의 메인 컴포넌트들, 즉, 딕셔너리 학습 컴포넌트(dictionary learning component) 및 "딕셔너리 공간(dictionary space)"을 형성하는 기저 벡터들로 입력 데이터를 표현하기 위한 좌표 컴포넌트(coordinate component)로 이루어진다. 딕셔너리 학습 부분의 경우, 아래의 최적화 문제가 해결된다:
Figure pct00002
X의 작은 노름에 도달하기 위해 딕셔너리 D가 임의적으로 크게 성장하지 않는다는 것을 보장하기 위하여 제약
Figure pct00003
이 부과된다. 딕셔너리 학습 문제의 제 1 항은 "재구성" 항(즉, 기저 벡터들의 세트로 데이터를 표현함)이고, 제 2 항은 (α의 선택된 값의 함수로써) 표현에서 희소성(sparseness)을 촉진하기 위해 포함된다.
최적화 문제는 D와 X에서 양면볼록(biconvex)인데; 즉, 하나의 변수가 고정될 때, 나머지 최적화 문제는 볼록 최적화 문제이다. 그러므로, 최적의 해결책으로의 수렴이 보장되도록 그리고 승수들의 교번 방향 방법(ADDM; Alternating Direction Method of Multipliers)과 같은 접근법들(또는 상대적으로 또한 강건하고 데이터의 임의적-스케일(arbitrary-scale)의 최적화를 제공할 수 있는 다른 접근법들)이 바람직한 것으로 간주되도록, 각각의 변수에 대해 교번적인 방식으로 해결함으로써 수용가능한 결과가 제공되는 것이 뒤따른다. 궁극적으로, 센서 판독 데이터를 기저 벡터들의 선형 조합으로서 정의하기 위해 사용될 수 있는 기저 벡터들의 딕셔너리 D가 획득된다. 실제 데이터에 기반하여 벡터들이 실시간으로 구축되는(built) 온라인-학습 딕셔너리(online-learning dictionary)를 포함한 임의의 적절한 딕셔너리 학습 알고리즘이 활용될 수 있다. 부가적으로, 딕셔너리를 초기화하기 위해 임의적인 값들을 사용하기보다는, 바람직한 접근법은 k-평균 알고리즘(k-means algorithm)을 사용하여 데이터 클러스터(data cluster)들의 적절한 센터(center)들을 발견하여, 그 센터들을 프로세스를 초기화하는 데 사용하는 것이다. 이 기법은 학습 프로세스에 대한 수렴의 속도를 개선하는 것으로 고려된다.
이 딕셔너리는, "정상 값"의 임의의 착신 판독(incoming reading)이 기저 벡터들의 선형 조합으로서 정확하게 표현될 수 있는 수용가능한 센서 판독들의 유니버스(universe)로 생각될 수 있다. 본 발명에 따르면, 이 딕셔너리는, 연관된 센서의 성능을 모니터링하고 딕셔너리 공간 외측에 있는 센서 판독들(이들 판독들은 "잔차들"로 지칭되며, 하기의 논의에서 E로 표기됨)의 존재를 검출하기 위해 희소 코딩 프로세스와 함께 사용된다. 즉, 주어진 딕셔너리 D에 대해, 본 발명의 프로세스는 아래의 관계식을 사용하여 착신 센서 판독들을 딕셔너리 D에 대해 평가함으로써 진행된다:
Figure pct00004
여기서
Figure pct00005
는 행렬 X의 각각의 행의 L-2 노름의 합이고,
Figure pct00006
는 행렬 X의 각각의 열의 L-1 노름 잔차들의 합이며, α 및 β는 작은 포지티브 스칼라(small positive scalar)들이다. 이 접근법은 또한 "다중 측정 벡터(multiple measurement vector)"(MMV) 접근법으로 알려져 있다. 항
Figure pct00007
은 "잔차"이며, 착신 센서 데이터가 수용 값들 너머로 벗어나서 센서가 아마도 "결함있음"으로서 평가될 필요가 있는지를 결정하는 데 사용된다. β 팩터는, 선택된 센서 판독들을 잔차로서 식별하는 감도의 레벨을 조정하기 위해 모니터링 직원에 의해 사용될 수 있다.
도시된 바와 같이, 비정상성 검출을 위해 사용되는 종래의 L-2 노름보다는, L-1 노름이 이러한 "잔차"를 정의하는 데 사용된다. 이는, L-2 노름이 상당한 수의 큰 잔차들을 (그리고 결과적으로, 결함있는 센서들의 관점들에서 상당한 수의 허위 포지티브들을) 생성하기 때문이다. 아래에서 정의되는 바와 같은 데이터포인트 값(datapoint value)들의 합인 L-1 노름은, 큰 잡음 팩터를 포함할 수 있는 데이터를 평가할 때 바람직한 대책인 것으로 간주되며, 목표는 잡음 에러와 비정상 센서 판독들 간의 차이를 판별할 수 있는 것이다. 논의된 바와 같은 L-1 노름의 사용은 잔차에서 희소성을 촉진할 뿐만 아니라 데이터에 대한 더 양호한 피팅(fit)을 제공한다. 그에 반해, 잡음있는 데이터(이 경우, 여기서는 가스 터빈 센서 데이터)를 처리할 때 L-2 노름이 사용될 경우, 결과적인 잔차들은 (선호되는 더 작은 잔차 값들 대신에) 매우 많으며 상대적으로 큰 경향이 있을 것이다.
도 3 및 도 4는 잔차들을 생성하는 데 있어서 L-2 노름을 사용하는 것(도 3에 도시되는 바와 같음)과 L-1 노름을 사용하는 것(도 4에 도시되는 바와 같음) 간의 결과들의 차이를 예시한다. L-2 노름이 상당히 많은 수의 데이터 포인트들을 "잔차들"로서 식별하는 것을 초래하고, 이는 바람직하지 않은 수의 허위 포지티브 "결함있는 센서" 알람들을 초래할 수 있다는 것이 명백하다. 반면에, L-1 노름을 사용함으로써 생성된 잔차들은 더 적고 훨씬 더 뚜렷하다. 이들 그래프(graph)들과 관련된 "실험" 설명 세부사항들은 아래에서 상세하게 설명될 것이다.
본 발명의 방법론의 논의를 다시 참조하면, 수학식(2)으로 표현되는 바와 같은 문제는 ADMM 최적화 기법을 적용함으로써 해결된다. 이 경우, 잔차들
Figure pct00008
의 세트의 형태를 취하는 해결책은 노름으로부터의 경계들 외측에 있는 센서 판독들의 세트이며, 센서의 결함있는 동작과 연관될 수 있다. 아래의 목록에서 보여지는 바와 같이, 반복적 방식(아래의 단계 3)으로 잔차 E가 업데이트될(updated) 때 종래의 수축 연산자(shrinkage operator)(아래의 단계 2에서 사용되는 바와 같음)가 엘리먼트-방식 "소프트 임계처리(soft thresholding)" 연산자로 대체되도록, MMV 알고리즘(위에서 양면볼록 변수들 D 및 X를 논의할 때 논의된 바와 같음)이 변경된다. 소프트 임계처리는, 여전히 데이터를 가능한 한 깨끗하게 분할하지만, 최대 마진(maximum margin)을 수정하고 그리고 일부 잘못 라벨링된(mislabeled) 데이터가 생성되게 하는 "느슨한(slack)" 변수를 도입하는 수정된 절차이다. 아래의 목록은 희소 코딩의 L-1 잔차들을 컴퓨팅(compute)하기 위해 사용되는 알고리즘을 상세하게 보여주며, 여기서 바로 아래에 포함된 테이블(table)은 알고리즘에 포함된 파라미터들을 정의한다:
Figure pct00009
희소 코딩의 L-1 잔차들을 컴퓨팅하기 위한 알고리즘 ( 잔차들은 가능하게는 잠재적인 센서 고장을 표시하는 것으로 분석됨)
해결:
Figure pct00010
입력: 신호들
Figure pct00011
, 딕셔너리
Figure pct00012
초기화:
Figure pct00013
로 설정
수렴까지 반복:
1.
Figure pct00014
,
여기서 L은 증강된-라그랑지안(Augmented-Lagrangian)이고 아래와 같이 정의됨:
Figure pct00015
2.
Figure pct00016
,
여기서
Figure pct00017
임. 이 최적화 문제는 행 수축 연산자(row shrinkage operator)를 사용함으로써 해결될 수 있음:
Figure pct00018
3.
Figure pct00019
,
여기서
Figure pct00020
임. 이 최적화 문제는 엘리먼트-방식 소프트 임계처리 연산자
Figure pct00021
를 사용함으로써 해결될 수 있으며, 여기서 아래와 같음.
Figure pct00022
4.
Figure pct00023
5.
Figure pct00024
6.
Figure pct00025
7.
Figure pct00026
인 경우, 중지됨.
출력:
Figure pct00027
이 알고리즘에 의해 생성되고, "잔차들"로서 정의된
Figure pct00028
의 값들은 가스 터빈 센서의 결함있는 동작과 연관된 잠재적 데이터 포인트들이다. 도 2의 흐름도를 다시 참조하면, 생성된 잔차들은 단계(220)로부터 출력된다. 이들 잔차들을 평가하고 그리고 시간 컴포넌트를 고려하기 위해, 본 발명에 따라 활용될 수 있는 적어도 2개의 상이한 접근법들이 존재한다.
제 1 접근법에서, 방사형 기저 함수(RBF; radial basis function) 커널(kernel)을 사용하여 일-클래스 지원 벡터 머신(SVM; support vector machine) 모델을 학습하기 위해 트레이닝 잔차가 사용된다. 일반적으로 말하면, SVM 모델은 데이터(이 경우, 잔차들)를 분석하고 패턴들을 인식하기 위해 학습 알고리즘들과 관련하여 사용된다. SVM은 개별 카테고리(category)들의 모델을 구축하고 그 후에 새로-도착하는 데이터를 적절한 카테고리에 할당한다. 여기서 사용되는 바와 같은 "일-클래스" SVM은 단지 단일 카테고리만을 생성하고, 새로-도착하는 데이터를 평가하도록 기능하고, 데이터가 카테고리 내에 있는지의 여부를 결정한다. SVM 알고리즘의 효율성을 증가시키기 위해, 전체 세트 대신에 트레이닝 잔차들의 샘플링(sampling)을 사용하는 것이 가능하다. 그 후에, SVM "일 클래스" 모델은 새로-도착하는 잔차들에 적용된다. 이 모델에서, 잔차 데이터 대부분은 생성된 모델의 경계들 내에 속할 것이고, 상대적으로 소수의 값들만이 "이상점들"로 남겨진다. 그 후에, 이들 이상점들은, 센서 그 자체의 추가의 평가를 프롬프트(prompt)할 수 있는 "비정상" 센서 판독들로서 정의된다. 실제로, 비정상 판독들은, 센서 고장이 이미 발생했다고 표시할 수 있거나, 또는 대안적으로, 센서가 고장나기 시작하고 주의가 필요하다는 예측으로서 사용될 수 있다.
허위 포지티브들의 수를 감소시키기 위해, 이전의 T 시간 단계들에 대한 센서 판독 데이터를 조사하는 것이 제안된다. 조사되는 데이터의 높은 퍼센티지(p%로 표기됨)가 또한 (생성된 일-클래스 SVM 모델에 따라) 비정상으로서 정의되는 경우, 현재의 데이터 포인트는 개연성있는 "결함있는 센서"를 경보하는 것으로 정의된다.
도 5의 흐름도는 센서 결함 검출의 이러한 예시적인 프로세스를 간략히 서술한다. 흐름도는, 특정 센서에 대한 데이터가 조사되고 있고 최대 센서 판독 값의 20% 미만인 판독들 및 이상점들을 제거하기 위해 데이터가 클리닝된다고(cleaned) 가정한다.
프로세스는 단계(500)에서, 다양한 가스 터빈 센서들로부터 원시 데이터(raw data)를 수집함으로써 시작된다(위에서 논의된 바와 같이, 이 데이터는 제어기(8)의 메모리 컴포넌트 내에 저장됨). 논의를 위해, 터빈의 생성된 출력 전력(MW)과 연관된 센서가 평가되는 것으로 가정된다. 단계(510)에서 도시되는 바와 같이, 데이터는 바람직하게, 이 파라미터와 연관되지 않은 센서들뿐만 아니라, 데이터가 전력 센서 출력과 강하게 상관되지 않은 임의의 센서들을 제거하기 위해 사전-프로세싱된다(pre-processed).
일단 사전프로세싱되면(preprocessed), 나머지 데이터는 위에서 간략히 서술된 바와 같이 희소 코딩 프로세스를 겪는다. 도 5에 도시된 바와 같이, 이 프로세스는 먼저 이 데이터에 대한 기저 벡터들의 딕셔너리를 구축함으로써(단계(520)), 그리고 그 후에 수학식(2)에 도시된 바와 같은 비정상성 검출과 함께 희소 코딩 알고리즘을 센서 판독 데이터에 적용하여, 이 데이터로부터 L-1 노름 잔차들을 컴퓨팅함으로써(단계(530)) 시작된다. 딕셔너리를 구축하는 다양한 방법들이 단계(520)에서 사용될 수 있지만, 효율적인 방식으로 수렴되는 바람직한 프로세스는, 데이터를 클러스터링(clustering)하는 단계(단계(522)), 클러스터들의 센터들을 발견하기 위해 k-평균 알고리즘을 사용하는 단계(단계(524)), 및 이들 "센터" 값들을 이용하여 딕셔너리를 초기화하는 단계(526)를 포함한다. 일단 구축되면, 딕셔너리는 L-1 노름 잔차들을 식별하기 위해 단계(530)의 희소 코딩 프로세스에서 사용된다. 일단 식별되면, L-1 노름 잔차들은 그 후에, 이들 잔차들 중 어느 것이, 결함있는 센서를 표시하는 정말로 "비정상" 센서 판독들인지를 결정하기 위해 데이터의 시간-기반 조사를 포함하여 추가로 평가된다.
도 5를 참조하면, L-1 노름 잔차들의 세트는 먼저, 분석의 나머지에서 수반되는 데이터포인트들의 수를 감소시키기 위해 샘플링될(sampled) 수 있다(단계(540)). 이 단계는 선택적이지만, 다음의 단계(550)에서 사용되는 SVM 알고리즘의 효율성을 개선하는 것으로 고려된다. 실제로, SVM 알고리즘은 그 후에, 잔차들의 샘플링된 그룹을 "트레이닝 데이터"로서 사용하고, 적절한 일-클래스 모델을 구축하며, 적절한 일-클래스 모델은 그 후에, 새로-도착한 잔차들을 클래스 내에 있는 것으로 또는 클래스 외측에 존재하는 것(이상점 데이터는 본 목적들을 위해 "비정상"으로서 정의됨)으로 분류하는 데 사용된다.
개발된 일-클래스 모델의 정확도를 검증하기 위해(그리고 도 2의 일반적 흐름도에서 간략히 서술된 바와 같이), 테스트 데이터의 세트는 일-클래스 SVM 모델을 통과한다. 도 5를 참조하면, 테스트 데이터는 (단계(520)에서 구축된 딕셔너리를 사용하여) L-1 잔차들의 테스트 데이터의 세트를 컴퓨팅하기 위해 프로세싱된다(단계(600)). 그 후에, L-1 잔차들의 테스트 데이터 세트는 단계(550)에서 생성된 SVM 모델을 통과한다(단계(610)로서 도시됨). (즉, 단계(620)에서의 평가 동안에 SVM 모델의 경계들 외측에 있는) "비정상"으로서 식별되는 임의의 데이터 포인트들이 수집 및 저장된다(단계(630)).
본 발명에 따르면, 분석과 연관된 시간 팩터가 이제 도입되며, 여기서 단계(640)에서 이전의 T 시간 단계들의 세트에 대한 데이터의 조사가 수행된다("T"는 분석을 수행하는 개인에 의해 결정된 값임). 그 후에, 단계(650)에서 결정이 이루어지는데, 여기서 이전의 T 시간 단계들 중 p% 초과가 또한 "비정상" 데이터 포인트들을 포함한 경우(또한, p는 사용자에 의해 선택된 값임), 연관된 센서는 아마도 "결함있음"으로 정의되며, 실제 센서의 평가가 적절하다(단계(660)). 그렇지 않고, 이전의 T 시간 단계들 중 p% 미만이 이들 비정상 판독들을 갖지 않는 경우, 현재의 비정상 판독은 잡음있는 판독, 가짜 에러 등에 기인할 수 있다고 결론내려질 수 있으며; 어느 경우이든, p% 임계치가 초과되지 않으면, 센서가 적절하게 동작하고 있고 어떠한 추가의 평가도 요구되지 않는 것으로 추정될 수 있다(단계(670)).
일-클래스 SVM 모델을 사용하는 대신에, 잔차 데이터를 추가로 연구하고 데이터 중 임의의 데이터가 관련 센서의 결함있는 성능과 연관될 수 있는지를 확인하기 위해, (잔차들을 지정하기 위해 L-1 노름을 활용하는) 희소 코딩 방법론을 다른 종류들의 분석들과 조합하는 것이 가능하다. 도 6은 하나의 이러한 대안적 분석 프로세스의 흐름도이며, 이러한 경우 여기서, 착신 센서 판독들이 결함있는 센서를 표시하는지를 결정하기 위해 일-클래스 SVM 모델 대신에, 임계치 분석이 사용된다. 도 6의 흐름도에서 도시되는 바와 같이, 프로세스의 초기 단계들은 도 5에 도시된 것들과 동일한데; 즉, 트레이닝 데이터를 사전프로세싱하고(단계(510)), 기저 벡터들의 딕셔너리를 구축하고(단계(520)), 희소 코드 프로세싱을 활용하여(단계(530)), L-1 잔차들의 세트를 컴퓨팅한다.
일단 트레이닝 데이터와 연관된 잔차들이 컴퓨팅되면, 도 6에 도시된 바와 같은 프로세스는, 단계(700)에서, 트레이닝 데이터에서 식별된 잔차들의 세트를 분석하고, 임계치를 초과하는 임의의 값이 "비정상" 데이터 포인트로서 지정되도록 그룹에 대한 적절한 임계값을 결정함으로써 계속된다. 이 임계치는 모든 센서들에 대해 균일할 수 있거나 또는 센서 단위 기반으로 정의될 수 있다. 예컨대, 5%의 가장 큰 잔차 값들이 "비정상"으로서 정의되고 나머지 잔차들은 무시되도록, 95%의 임계값이 사용될 수 있다.
프로세스의 이 시점에서, 단계(520)에서 구축된 것과 동일한 딕셔너리를 사용하여 그리고 단계(530)에 의해 예시된 동일한 희소 코딩 프로세스를 수행하여, 테스팅 데이터(testing data)에 대한 L-1 잔차들이 다음에 결정된다(단계(710)). 일단 테스팅 데이터에 대한 L-1 잔차들이 컴퓨팅되면, 이들은 단계(700)에서 정의된 임계치에 대해 (단계(720)에서) 비교된다. 주어진 테스트 데이터 잔차가 (결정(725)에서) 임계치 미만인 경우, 데이터는 단계(730)에서 "잡음있는" 데이터(또는 어떠한 추가의 연구도 요구되지 않는다는 것을 표시하는 일부 유사한 지정 유형)로 지정된다.
단계(725)에서의 임계치 테스트의 결과가 포지티브라고 가정하면, 임계치-초과 잔차는 "비정상"으로 지정되고 저장된다(단계(740)에서 도시된 바와 같음). 분석의 시간 양상이 이 시점에서 도입되고, 이전의 T2 시간 단계들에 대한 데이터가 조사된다(단계(750)). 그 후에, 다음 단계(760)에서 결정이 이루어지는데, 여기서 T2 이전의 결과들 중 q%가 모두 임계치를 초과한 경우, 현재의 데이터 포인트가 센서 결함을 표시한다는 결정이 이루어지고, 센서 그 자체의 추가의 연구가 정당화된다(단계(770)). 그렇지 않고, 결정 단계(760)에서 임계치가 초과되지 않는 경우, 현재의 데이터 포인트는 "잡음있음"(또는 유사한 지정)으로 라벨링되고, 조사 프로세스가 계속된다.
도 6에서 간략히 서술된 바와 같은 프로세스의 바람직한 실시예에서, q=100의 값이 사용되어, 특정 센서가 결함이 있을 가능성이 있는 것으로 평가되는 것을 요청하기 전에 과거 T2 값들 모두가 "비정상"으로서 지정되는 것이 요구된다는 것을 주목해야 한다.
결함있는 가스 터빈 센서를 검출하는 데 활용되는 희소 코딩의 상세한 예들
센서 결함을 검출하기 위한 본 발명의 희소 코딩 접근법의 정확도를 결정하기 위한 실험적 접근법에서, 주어진 가스 터빈에 대한 이력 데이터의 세트가 연구되었다. 센서 데이터(MW)는, 머신이 계속 작동되고 있는지의 여부 및 판독들이 연속적인 실제 값들인지 여부를 표시하는 센서이다. 데이터는 다음과 같이 사전-프로세싱되었는데: (1) MW 출력의 최대 값 아래의 20% 미만의 값들의 명백한 이상점들 및 다른 데이터 포인트들이 제거되었고; (2) (예컨대, "ROC 미만 영역" 기준들을 사용하여) MW 판독들과 관련이 없는 리던던트 센서들(및 그들의 데이터)을 제거하기 위해 피처 선택(feature selection)이 수행되었고; 그리고 (3) MW와의 약한 선형 상관 관계를 나타낸 모든 나머지 센서들이 제거되었다. 궁극적으로, 15개의 센서들의 세트가 추가의 평가를 위해 남았다.
나머지 15개의 센서들과 연관된 데이터의 세트는 2개의 그룹들로 분할되었으며, 제 1 그룹은 트레이닝 데이터로서 사용되고, 제 2 그룹은 테스팅 데이터로서 사용되었다. 트레이닝 데이터는 에러가 없는 것으로 가정되었다. 트레이닝 데이터의 랜덤하게(randomly) 선택된 위치들에서 랜덤 에러(random error)(이 경우, 증가하는 램프(increasing ramp)의 형태로)가 도입되었다.
도 2, 도 5 및 도 6의 흐름도들에서 간략히 서술된 바와 같은 프로세스들이 데이터를 프로세싱(process)하기 위해 사용되었으며, 희소 코딩 알고리즘으로서 수학식(2)에서 사용된 α = 0.001 및 β = 0.01의 값들이 데이터에 적용되고, 잔차들이 계산되었다. "비정상" 데이터(즉, 잠재적으로 결함있는 센서)를 정의하기 위한 임계치 방법들 및 일-클래스 SVM 모델 둘 모두가 평가되었다. 일-클래스 SVM 모델의 경우, 3000개의 트레이닝 잔차들의 샘플링된 세트가 사용되어 모델을 형성하였다. SVM 모델의 시간-기반 컴포넌트의 경우, p=50의 값이 선택되었고, T=50의 값이 사용되었다(즉, 결함있는 센서의 가능성을 선언하기 위해서는 과거 50 시간 단계들에 대한 데이터의 적어도 50%가 "비정상"이도록 요구되었음). 임계치 평가 프로세스의 경우, q=100의 값이 선택되었고, T2=24와 함께 사용되었다(즉, 현재의 데이터 포인트를 "비정상"으로서 지정하고 결함있는 센서의 가능성을 표시하기 위해서는 과거 24 시간 단계들에 대한 잔차들의 100%가 임계치를 초과해야 함).
위에서 언급된 바와 같이, 본 발명의 현저한 양상은, 희소 코딩 프로세스들과 함께 통상적으로 사용되는 L-2 노름 대신에, L-1 노름을 사용하여 잔차들을 계산하는 것이다. 일반적으로, (행렬에서 벡터들로서 디스플레이되는(displayed)) 데이터 세트의 "노름"은 아래와 같이 정의될 수 있다:
Figure pct00029
L-1 노름의 경우(여기서 j=1), 이 관계식은 아래와 같이 단순화되고:
Figure pct00030
그리고 종종 "최소 절대 편차(least absolute deviation)"(LAD) 또는 "최소 절대 에러(least absolute error)"(LAE)로 지칭된다. 이와 비교하여, L-2는 아래와 같이 종래의 유클리드 "거리" 측정(Euclidean "distance" measure)으로서 표현될 수 있다:
Figure pct00031
가스 터빈 센서들에 의해 수집된 데이터가 상대적으로 큰 잡음 컴포넌트를 포함할 수 있기 때문에, L-2 노름보다 L-1 노름이 이러한 원하지 않는 정보에 덜 민감할 것이라고 결정되었다. 비교를 위해, L-2 노름을 사용하여 계산들의 세트가 또한 준비되었다. 위에서 언급된 도 3 및 도 4에 대한 참조가 이제 이루어지는데, 여기서 이들 2개의 플롯(plot)들은 L-1 노름을 사용할 때 획득된 결과들의 개선을 예시한다. 도 3은 L-2 잔차들로서 상대적으로 많은 수의 스파이크(spike)들을 도시하는, 희소 코딩 프로세스의 결과들의 플롯이다. 이에 반해, (본 발명에 따라 L-1 노름을 사용하는 것에 기반한) 도 4의 플롯은 훨씬 더 적은 수의 잔차들을 포함한다. 명백하게, L-2 노름은 더 많은 스파이크들을 생성하는데, 이는 실제 센서 판독들에서 작은 편차들을 검출하는 능력이 검출을 하는 것을 더 어렵게 만들며, 여기서 이러한 미묘한 데이터는 센서 고장을 더 많이 표시할 수 있고 본 발명의 목적들에 있어서 더 중요할 수 있다.
또한 비교를 위해, "비정상" 데이터의 최종 결정을 하는 시간 컴포넌트를 무시한 결과들의 세트가 생성되었다. 아래의 테이블에서, 이 결과들의 세트는 "베이스라인 알고리즘(baseline algorithm)"으로 식별된다. 테이블에서, 잔차들(비정상성들)은 "+"로 표기되고, 정상 데이터는 "-"로 표기된다. 특정 태블러 형태(tabular form)는, 정상 데이터가 비정상성으로서 잘못 분류된 경우뿐만 아니라 비정상성이 정상 데이터로 잘못 분류된 경우의 인스턴스(instance)들의 형태로 "허위 포지티브들"을 (헤더(header)들 "예측" 및 "실제"를 사용하여) 보여준다.
Figure pct00032
위의 테이블의 조사로부터, 본 발명의 접근법들(SVM 모델 및 임계치) 둘 모두가, 베이스라인 접근법과 비교할 때 (더 낮은 수들의 허위-포지티브들의 관점들에서) 개선된 결과들을 제공한다는 것이 명백하다. (L-2 및 L-1 노름 둘 모두를 사용한 각각의 인스턴스에서 계산된 값들을 포함하여) 총 10개의 랜덤 시간 포인트들에 대한 더 완전한 세트의 결과들이 본 명세서의 부록(appendix)에서 보여진다. 부록에 포함된 완전한 세트의 결과들은 추가로, 각각의 인스턴스에서 L-2 노름 대신에 L-1 노름의 사용이 더 정확한 결과를 제공한다는 것을 보여주며, 그에 따라 L-1 노름의 선택이 더 양호한 연산임을 입증한다.
위에서 설명된 방법은 위에서 설명된 바와 같이 컴퓨터(computer)에 의해 실행되는 프로그램 모듈(program module)들에 의해 구현될 수 있다. 일반적으로, 프로그램 모듈들은, 특정 태스크(task)들을 수행하거나 또는 특정한 추상적 데이터 유형들을 구현하는 루틴(routine)들, 오브젝트(object)들, 컴포넌트들, 데이터 구조들 등을 포함한다. 본원에서 사용되는 바와 같은 "프로그램"이라는 용어는 단일 프로그램 모듈 또는 협력하여 동작하는 다수의 프로그램 모듈들을 함축할 수 있다. 본 개시내용은, 개인 컴퓨터(PC; personal computer)들, 핸드-헬드 디바이스(hand-held device)들, 멀티-프로세서 시스템(multi-processor system)들, 마이크로프로세서-기반 프로그램가능 가전제품(microprocessor-based programmable consumer electronics), 네트워크 PC(network PC)들, 미니-컴퓨터(mini-computer)들, 메인프레임 컴퓨터(mainframe computer)들 등을 포함한 다양한 유형들의 컴퓨터들 상에서 구현될 수 있다. 본 개시내용은 또한, 통신 네트워크를 통해 링크된(linked) 원격 프로세싱 디바이스(remote processing devices)들에 의해 태스크들이 수행되는 분산형 컴퓨팅 환경들에서 이용될 수 있다. 분산형 컴퓨팅 환경에서, 모듈들은 로컬(local) 및 원격 메모리 저장 디바이스들 둘 모두에 위치될 수 있다.
위의 방법론을 구현하기 위한 예시적인 프로세싱 모듈은 하드와이어되거나(hardwired) 또는 개별 메모리에 저장될 수 있으며, 개별 메모리는, ROM 또는 다른 유형의 하드 자기 드라이브(hard magnetic drive), 광학 저장소, 테이프(tape) 또는 플래시 메모리(flash memory)와 같은 컴퓨터 판독가능 매체로부터 프로세서 또는 복수의 프로세서들의 메인 메모리(main memory)로 판독된다. 메모리 매체들에 저장된 프로그램의 경우, 모듈에서의 명령들의 시퀀스(sequence)들의 실행은 프로세서로 하여금 본원에서 설명된 프로세스 단계들을 수행하게 한다. 본 개시내용의 실시예들은 하드웨어와 소프트웨어의 임의의 특정 조합으로 제한되지 않으며, 전술한 내용을 구현하기 위해 요구되는 컴퓨터 프로그램 코드는 당업자에 의해 개발될 수 있다.
본원에서 이용되는 바와 같은 "컴퓨터-판독가능 매체"라는 용어는, 명령들을 하나 또는 그 초과의 프로세서들에 제공하는 또는 제공하는 데 참여하는 임의의 유형적 머신-인코딩된 매체(tangible machine-encoded medium)를 의미한다. 예컨대, 컴퓨터-판독가능 매체는 하나 또는 그 초과의 광학 또는 자기 메모리 디스크(disk)들, 플래시 드라이브(flash drive)들 및 카드(card)들, 판독-전용 메모리, 또는 통상적으로 메인 메모리를 구성하는 DRAM과 같은 랜덤 액세스 메모리(random access memory)일 수 있다. 이러한 매체들은 유형적이지 않은 전파된 신호들을 제외한다. 캐시된(cached) 정보는 컴퓨터-판독가능 매체 상에 저장되는 것으로 고려된다. 컴퓨터-판독가능 매체들의 통상의 수단들은 당해 기술 분야에 잘 알려져 있으며, 여기서 상세하게 설명될 필요가 없다.
전술한 상세한 설명은 모든 관점에서 제한적인 것이 아니라 예시적이고 대표적인 것으로서 이해되어야 하며, 본원의 본 개시내용의 범위는 상세한 설명으로부터 결정되는 것이 아니라, 특허법에 의해 허용되는 최대 범위에 따라 해석되는 청구항들로부터 결정될 것이다. 본 개시내용의 범위 및 사상으로부터 벗어남이 없이 다양한 수정들이 당업자들에 의해 구현될 수 있다는 것이 이해되어야 한다.
부록 - 잔차에 대한 영향; L-2 노름을 L-1 노름과 비교
Figure pct00033
Figure pct00034
Figure pct00035
Figure pct00036
Figure pct00037
Figure pct00038
Figure pct00039
Figure pct00040
Figure pct00041
Figure pct00042

Claims (20)

  1. 적어도 하나의 센서(sensor)로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장(gas turbine sensor failure)을 검출하기 위한 방법으로서,
    정상 동작 상태(normal operating condition)들에 대해 알려진 센서 판독들과 연관된 값들을 정의하는 기저 벡터(basis vector)들의 딕셔너리(dictionary)를 생성하는 단계;
    L-1 노름 잔차 데이터(L-1 norm residual data)의 세트(set)를 식별하기 위해, 생성된 딕셔너리를 사용하여, 희소 코딩 프로세스(sparse coding process)를 센서 판독 데이터의 세트에 적용하는 단계;
    최대-값의 L-1 노름 잔차 데이터(largest-valued L-1 norm residual data)의 미리 결정된 서브세트(subset)를 비정상 센서 판독(abnormal sensor reading)들로 카테고리화(categorize)하기 위해 상기 L-1 노름 잔차 데이터를 평가하는 단계;
    상기 비정상 센서 판독들을 복수의 이전(prior) 센서 판독들과 비교하고 그리고 상기 이전 센서 판독들 중 미리 정의된 개수의 이전 센서 판독들이 또한 비정상 센서 판독들로서 지정되는 경우, 지정된 비정상 센서 판독들을 센서 고장과 연관된 것으로 정의하는 단계;
    특정 센서를 식별하는 센서 고장 신호를 가스 터빈 직원(gas turbine personnel)에게 송신하는 단계; 및
    고장난 센서를 서비스(service)로부터 제거하고 그리고 상기 고장난 센서를 수리 또는 교체하는 단계를 포함하는,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  2. 제 1 항에 있어서,
    상기 딕셔너리를 생성하는 단계에서,
    상기 센서 판독 데이터를 복수의 클러스터(cluster)들로 조직화하는 단계,
    k-평균 알고리즘(k-means algorithm)을 사용하여 각각의 클러스터의 센터(center)를 결정하는 단계, 및
    각각의 클러스터의 결정된 센터를 이용하여 상기 딕셔너리를 초기화하는 단계가 수행되는,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  3. 제 1 항에 있어서,
    상기 희소 코딩 프로세스는 아래의 관계식,
    Figure pct00043
    을 수행하고,
    여기서 D는 생성된 딕셔너리이고, X는 복수의 센서 판독들이고, E는 L-1 노름 잔차들이고, 그리고 α 및 β는 작은 포지티브 스칼라(small positive scalar)들인,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  4. 제 3 항에 있어서,
    최소화 함수(minimization function)는 아래의 반복적 단계들,
    a.
    Figure pct00044
    ― 여기서 L은 증강된-라그랑지안(Augmented-Lagrangian)이고 그리고
    Figure pct00045
    로서 정의됨 ―,
    b.
    Figure pct00046
    ― 여기서
    Figure pct00047
    이고, 행 수축 연산자(row shrinkage operator), 즉,
    Figure pct00048
    를 사용하는 경우 ―,
    c.
    Figure pct00049
    ― 여기서
    Figure pct00050
    Figure pct00051
    이고, 여기서
    Figure pct00052
    임 ―,
    d.
    Figure pct00053

    e.
    Figure pct00054

    f.
    Figure pct00055

    g.
    Figure pct00056
    인 경우 중지됨 ― 여기서
    Figure pct00057
    는 L-1 정상 잔차 데이터의 세트로서 정의되고, Z는 보조 변수이고, M은 라그랑주 승수(Lagrange multiplier)이고,
    Figure pct00058
    는 페널티 계수(penalty coefficient)이고,
    Figure pct00059
    Figure pct00060
    를 업데이트(updating)하기 위한 파라미터(parameter)이고, 그리고
    Figure pct00061
    는 알고리즘을 종결시키기 위해 Z 내지 X의 근사(closeness)를 측정하는 데 사용되는 에러 허용한계(error tolerance)임 ―
    을 수행하는,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  5. 제 1 항에 있어서,
    상기 평가를 수행하는 단계에서,
    일-클래스 지원 벡터 머신 모델(one-class Support Vector Machine model)이 생성되고, 그리고 각각의 잔차는 상기 클래스 내에 있는 것 또는 상기 클래스 외측에 있는 것으로 분류되고, 상기 클래스 외측의 L-1 노름 잔차들은 비정상 센서 판독들로서 정의되는,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  6. 제 5 항에 있어서,
    상기 일-클래스 SVM 모델을 생성하기 전에 상기 데이터 세트의 크기를 감소시키기 위해, 상기 L-1 노름 잔차 데이터의 세트가 샘플링되는(sampled),
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  7. 제 5 항에 있어서,
    상기 비교를 수행하는 단계에서,
    복수의 이전 시간 단계 데이터(prior time step data)가 조사되고(reviewed), 그리고 상기 복수의 이전 시간 단계 데이터의 p%가 또한 비정상 센서 판독으로서 지정되는 경우, 센서 고장 결정이 이루어지는,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  8. 제 7 항에 있어서,
    상기 p의 값은 적어도 50과 동일한,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  9. 제 1 항에 있어서,
    상기 L-1 노름 잔차 데이터는, 임계치 퍼센티지(threshold percentage)를 정의하고 그리고 정의된 임계치를 초과하는 모든 값들을 비정상 센서 판독 데이터로서 정의함으로써 평가되는,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  10. 제 9 항에 있어서,
    상기 비교를 수행하는 단계에서,
    복수의 이전 시간 단계 데이터가 조사되고, 그리고 상기 복수의 이전 시간 단계 데이터의 q%가 또한 비정상 센서 판독으로서 지정되는 경우, 센서 고장 결정이 이루어지는,
    적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서 고장을 검출하기 위한 방법.
  11. 가스 터빈 센서의 고장을 검출하기 위한 시스템(system)으로서,
    센서 판독들의 데이터베이스(database);
    상기 센서 판독들의 데이터베이스와 통신하는 센서 모니터링 시스템 컴포넌트(sensor monitoring system component) ― 상기 센서 모니터링 시스템 컴포넌트는 프로그램 저장 디바이스(program storage device) 및 프로세서(processor)를 포함하고, 상기 프로그램 저장 디바이스는 방법 단계들을 수행하기 위해 상기 프로세서에 의해 실행가능한 프로그램 명령들의 세트를 고정형의 유형적 매체(fixed tangible medium)로 포함하고, 상기 방법 단계들은,
    정상 동작 상태들에 대해 알려진 센서 판독들과 연관된 값들을 정의하는 기저 벡터들의 딕셔너리를 생성하는 단계,
    L-1 노름 잔차 데이터의 세트를 식별하기 위해, 생성된 딕셔너리를 사용하여, 희소 코딩 프로세스를 센서 판독 데이터의 세트에 적용하는 단계,
    최대-값의 L-1 노름 잔차 데이터의 미리 결정된 서브세트를 비정상 센서 판독들로 카테고리화하기 위해 상기 L-1 노름 잔차 데이터를 평가하는 단계,
    상기 비정상 센서 판독들을 복수의 이전 센서 판독들과 비교하고 그리고 상기 이전 센서 판독들 중 미리 정의된 개수의 이전 센서 판독들이 또한 비정상 센서 판독들로서 지정되는 경우, 지정된 비정상 센서 판독들을 센서 고장과 연관된 것으로 정의하는 단계, 및
    수리 또는 교체되도록 하기 위해, 고장난 센서를 식별하는 정보를 포함하는 센서 고장 신호를 가스 터빈 직원에게 송신하는 단계임 ―;
    비정상(abnormal)으로서 분류된 센서 판독들을 저장하기 위한 모니터링 시스템 데이터베이스를 포함하는,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  12. 제 11 항에 있어서,
    상기 실행가능 프로그램 명령들은, 상기 프로세서에 의해 실행될 때,
    상기 센서 판독 데이터를 복수의 클러스터들로 조직화하고,
    k-평균 알고리즘을 사용하여 각각의 클러스터의 센터를 결정하고, 그리고
    각각의 클러스터의 결정된 센터를 이용하여 상기 딕셔너리를 초기화함으로써, 기저 벡터들의 딕셔너리를 생성하는 단계를 수행하도록 상기 프로세서에 추가로 명령하는,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  13. 제 10 항에 있어서,
    상기 실행가능 프로그램 명령들은, 상기 프로세서에 의해 실행될 때,
    Figure pct00062
    에 의해, 희소 코딩 프로세스를 적용하는 단계를 수행하도록 상기 프로세서에 추가로 명령하고,
    여기서 D는 생성된 딕셔너리이고, X는 복수의 센서 판독들이고, E는 L-1 노름 잔차들이고, 그리고 α 및 β는 작은 포지티브 스칼라들인,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  14. 제 13 항에 있어서,
    상기 실행가능 프로그램 명령들은, 상기 프로세서에 의해 실행될 때,
    아래의 반복식,
    a.
    Figure pct00063
    ― 여기서 L은 증강된-라그랑지안이고 그리고
    Figure pct00064
    로서 정의됨 ―,
    b.
    Figure pct00065
    ― 여기서
    Figure pct00066
    이고, 행 수축 연산자, 즉,
    Figure pct00067
    를 사용하는 경우 ―,
    c.
    Figure pct00068
    ― 여기서
    Figure pct00069
    Figure pct00070
    이고, 여기서
    Figure pct00071
    임 ―,
    d.
    Figure pct00072

    e.
    Figure pct00073

    f.
    Figure pct00074

    g.
    Figure pct00075
    인 경우 중지됨 ― 여기서
    Figure pct00076
    는 L-1 정상 잔차 데이터의 세트로서 정의되고, Z는 보조 변수이고, M은 라그랑주 승수이고,
    Figure pct00077
    는 페널티 계수이고,
    Figure pct00078
    Figure pct00079
    를 업데이트하기 위한 파라미터이고, 그리고
    Figure pct00080
    는 알고리즘을 종결시키기 위해 Z 내지 X의 근사를 측정하는 데 사용되는 에러 허용한계임 ― 을 수행함으로써, 희소 코딩 잔차 검출 최소화(sparse coding residual detection minimization)를 수행하도록 상기 프로세서에 추가로 명령하는,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  15. 제 10 항에 있어서,
    상기 실행가능 프로그램 명령들은, 상기 프로세서에 의해 실행될 때,
    상기 평가하는 단계에서 일-클래스 지원 벡터 머신 모델을 생성하도록 상기 프로세서에 추가로 명령하며,
    각각의 L-1 노름 잔차는 상기 클래스 내에 또는 상기 클래스 외측에 있는 것으로 분류되고, 생성된 클래스 외측의 L-1 노름 잔차들은 비정상 센서 판독들로서 정의되는,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  16. 제 15 항에 있어서,
    상기 실행가능 프로그램 명령들은, 상기 프로세서에 의해 실행될 때,
    복수의 이전 시간 단계 데이터를 조사함으로써 상기 비교하는 단계를 수행하도록 상기 프로세서에 추가로 명령하고, 그리고
    상기 복수의 이전 시간 단계 데이터의 p%가 또한 비정상 센서 판독으로서 지정되는 경우, 센서 고장 결정이 이루어지는,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  17. 제 10 항에 있어서,
    상기 실행가능 프로그램 명령들은, 상기 프로세서에 의해 실행될 때,
    임계치 퍼센티지를 정의하고 그리고 정의된 임계치를 초과하는 모든 값들을 비정상 센서 판독 데이터로서 정의함으로써 상기 L-1 노름 잔차들을 평가하는 단계를 수행하도록 상기 프로세서에 추가로 명령하는,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  18. 제 17 항에 있어서,
    상기 실행가능 프로그램 명령들은, 상기 프로세서에 의해 실행될 때,
    복수의 이전 시간 단계 데이터가 조사되는 것에 의해 상기 비교하는 단계를 수행하도록 상기 프로세서에 추가로 명령하고, 그리고
    상기 복수의 이전 시간 단계 데이터의 q%가 또한 비정상 센서 판독으로서 지정되는 경우, 센서 고장 결정이 이루어지는,
    가스 터빈 센서의 고장을 검출하기 위한 시스템.
  19. 적어도 하나의 센서로부터 수집된 판독들에 기반하여 가스 터빈 센서의 고장을 예측하기 위한 방법으로서,
    정상 동작 상태들에 대해 알려진 센서 판독들과 연관된 값들을 정의하는 기저 벡터들의 딕셔너리를 생성하는 단계;
    L-1 노름 잔차 데이터의 세트를 식별하기 위해, 생성된 딕셔너리를 사용하여, 희소 코딩 프로세스를 센서 판독 데이터의 세트에 적용하는 단계;
    최대-값의 L-1 노름 잔차 데이터의 미리 결정된 서브세트를 비정상 센서 판독들로 카테고리화하기 위해 상기 L-1 노름 잔차 데이터를 평가하는 단계;
    상기 비정상 센서 판독들을 복수의 이전 센서 판독들과 비교하고 그리고 상기 이전 센서 판독들 중 미리 정의된 개수의 이전 센서 판독들이 또한 비정상 센서 판독들로서 지정되는 경우, 지정된 비정상 센서 판독들을 예측된 센서 고장과 연관된 것으로 정의하는 단계; 및
    잠재적 수리 또는 교체를 위해 센서를 식별하는 메시지(message)를 가스 터빈 직원에게 송신하는 단계를 포함하는,
    가스 터빈 센서의 고장을 예측하기 위한 방법.
  20. 제 19 항에 있어서,
    상기 희소 코딩 프로세스는 아래의 관계식,
    Figure pct00081
    을 수행하고,
    여기서 D는 생성된 딕셔너리이고, X는 복수의 센서 판독들이고, E는 상기 L-1 노름 잔차들이고, 그리고 α 및 β는 작은 포지티브 스칼라들인,
    가스 터빈 센서의 고장을 예측하기 위한 방법.
KR1020177009440A 2014-09-10 2015-09-03 희소 코딩 방법론을 활용한 가스 터빈 센서 고장 검출 KR101955305B1 (ko)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201462048588P 2014-09-10 2014-09-10
US62/048,588 2014-09-10
PCT/US2015/048265 WO2016040082A1 (en) 2014-09-10 2015-09-03 Gas turbine sensor failure detection utilizing a sparse coding methodology

Publications (2)

Publication Number Publication Date
KR20170045354A true KR20170045354A (ko) 2017-04-26
KR101955305B1 KR101955305B1 (ko) 2019-03-08

Family

ID=54186277

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020177009440A KR101955305B1 (ko) 2014-09-10 2015-09-03 희소 코딩 방법론을 활용한 가스 터빈 센서 고장 검출

Country Status (5)

Country Link
US (1) US10557719B2 (ko)
EP (1) EP3191797B1 (ko)
KR (1) KR101955305B1 (ko)
CN (1) CN107111311B (ko)
WO (1) WO2016040082A1 (ko)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10378936B2 (en) * 2015-09-04 2019-08-13 International Business Machines Corporation Identification of failed sensors in a system of interconnected devices
KR101765235B1 (ko) * 2016-11-28 2017-08-04 한국건설기술연구원 사물인터넷 기반 센서와 무인비행체를 이용한 시설물 유지관리 시스템 및 그 방법
US9996078B1 (en) 2017-04-07 2018-06-12 Pratt & Whitney Canada Corp. Pre-emptive fault detection through advanced signal analysis
US20180355797A1 (en) * 2017-06-12 2018-12-13 United Technologies Corporation Nacelle anti-icing troubleshooting for a two valve system
WO2019194698A1 (en) * 2018-04-04 2019-10-10 Siemens Aktiengesellschaft Validation device and method for validating sensor data
CN108828402A (zh) * 2018-04-26 2018-11-16 广东电网有限责任公司 一种高压电缆故障诊断的方法
US10967514B2 (en) * 2018-08-22 2021-04-06 Abb Schweiz Ag Apparatus and method to monitor robot mechanical condition
US11327476B2 (en) * 2018-11-13 2022-05-10 Microsoft Technology Licensing, Llc Sensor fall curve identification
EP3722998A1 (en) * 2019-04-11 2020-10-14 Teraki GmbH Data analytics on pre-processed signals
US11448570B2 (en) * 2019-06-04 2022-09-20 Palo Alto Research Center Incorporated Method and system for unsupervised anomaly detection and accountability with majority voting for high-dimensional sensor data
CN112131907A (zh) * 2019-06-24 2020-12-25 华为技术有限公司 一种训练分类模型的方法及装置
EP3779620A1 (en) * 2019-08-13 2021-02-17 Siemens Aktiengesellschaft Automatic calculation of measurement confidence in flexi-ble modular plants and machines
DE102019219820A1 (de) * 2019-12-17 2021-06-17 Robert Bosch Gmbh Verfahren, Computerprogramm, elektronisches Speichermedium und Vorrichtung zur Überwachung eines Sensorsystems.
GB201919198D0 (en) * 2019-12-23 2020-02-05 Univ Surrey Sensor fault prediction method and apparatus
CN111091159B (zh) * 2019-12-27 2023-03-17 哈尔滨工程大学 一种自主式水下机器人推进器弱故障特征提取方法
US11657122B2 (en) * 2020-07-16 2023-05-23 Applied Materials, Inc. Anomaly detection from aggregate statistics using neural networks
CN112098105B (zh) * 2020-08-10 2022-09-20 北京化工大学 基于机匣宽频振动信号的燃气轮机叶片故障监测预警方法
CN111950505B (zh) * 2020-08-24 2023-08-29 湖南科技大学 一种ssa-aann的风力发电机传感器状态评估方法
CN112098850B (zh) * 2020-09-21 2024-03-08 山东工商学院 基于sdo算法的锂离子电池电压故障诊断方法及系统
CN112257773B (zh) * 2020-10-19 2022-06-10 重庆邮电大学 基于无线传感器网络多测量向量的机械装备故障诊断方法
CN115524150B (zh) * 2022-09-13 2024-04-05 西安交通大学 基于稀疏时域同步平均的旋转机械故障检测方法
CN117389202A (zh) * 2023-12-13 2024-01-12 南京德克威尔自动化有限公司 一种用于智能工厂总线控制io模块

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070124113A1 (en) * 2005-11-28 2007-05-31 Honeywell International, Inc. Fault detection system and method using multiway principal component analysis

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59149513A (ja) 1983-02-04 1984-08-27 Toshiba Corp 制御系の異常判別方法
US5465321A (en) * 1993-04-07 1995-11-07 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Hidden markov models for fault detection in dynamic systems
US7082359B2 (en) * 1995-06-07 2006-07-25 Automotive Technologies International, Inc. Vehicular information and monitoring system and methods
JP3606446B2 (ja) 2000-06-06 2005-01-05 日本ビクター株式会社 音声符号化装置及びその方法
US6609036B1 (en) * 2000-06-09 2003-08-19 Randall L. Bickford Surveillance system and method having parameter estimation and operating mode partitioning
US6741919B1 (en) * 2003-02-26 2004-05-25 General Electric Company Methods and apparatus for detecting impending sensor failure
US7447609B2 (en) 2003-12-31 2008-11-04 Honeywell International Inc. Principal component analysis based fault classification
US7477960B2 (en) 2005-02-16 2009-01-13 Tokyo Electron Limited Fault detection and classification (FDC) using a run-to-run controller
US7233932B2 (en) 2005-05-31 2007-06-19 Honeywell International, Inc. Fault detection system and method using approximate null space base fault signature classification
US7677075B2 (en) 2006-09-29 2010-03-16 General Electric Company Methods and apparatus for evaluating sensors and/or for controlling operation of an apparatus that includes a sensor
GB2477324A (en) * 2010-02-01 2011-08-03 Rolls Royce Plc Device monitoring
CN101937207B (zh) * 2010-08-27 2012-09-05 上海交通大学 机械设备状态的智能可视化监测与诊断方法
US8290745B2 (en) 2010-12-17 2012-10-16 General Electric Company Systems and methods for identifying faulty sensors within a power generation system
CN103020935B (zh) * 2012-12-10 2015-09-02 宁波大学 一种自适应在线字典学习的图像超分辨率方法
CN103295242B (zh) 2013-06-18 2015-09-23 南京信息工程大学 一种多特征联合稀疏表示的目标跟踪方法
CN103592014B (zh) * 2013-11-06 2016-03-23 重庆工商大学 一种车载称重系统的传感器标定方法
US9568378B2 (en) * 2013-12-18 2017-02-14 Siemens Energy, Inc. Multi functional sensor system for gas turbine combustion monitoring and control

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070124113A1 (en) * 2005-11-28 2007-05-31 Honeywell International, Inc. Fault detection system and method using multiway principal component analysis

Also Published As

Publication number Publication date
CN107111311B (zh) 2020-01-14
US20180231394A1 (en) 2018-08-16
CN107111311A (zh) 2017-08-29
WO2016040082A1 (en) 2016-03-17
US10557719B2 (en) 2020-02-11
KR101955305B1 (ko) 2019-03-08
EP3191797B1 (en) 2018-08-01
EP3191797A1 (en) 2017-07-19

Similar Documents

Publication Publication Date Title
KR101955305B1 (ko) 희소 코딩 방법론을 활용한 가스 터빈 센서 고장 검출
CN107111309B (zh) 利用监督式学习方法的燃气涡轮机故障预测
KR101903283B1 (ko) 발전 설비의 자동 진단 시스템 및 자동 진단 방법
Que et al. A data-driven health prognostics approach for steam turbines based on xgboost and dtw
US9483049B2 (en) Anomaly detection and diagnosis/prognosis method, anomaly detection and diagnosis/prognosis system, and anomaly detection and diagnosis/prognosis program
WO2020049087A1 (en) Computer-implemented method, computer program product and system for anomaly detection and/or predictive maintenance
CN106404403B (zh) 用于涡轮机的分析的方法和系统
CN115034248A (zh) 用于设备的自动诊断方法、系统和存储介质
CN107710089B (zh) 工厂设备诊断装置以及工厂设备诊断方法
EP3105644A1 (en) Method of identifying anomalies
JP2010181188A (ja) プラント運転状態監視方法
JP2011070334A (ja) 異常診断装置及び異常診断方法
Calabrese et al. An event based machine learning framework for predictive maintenance in industry 4.0
Hong et al. Remaining useful life prediction using time-frequency feature and multiple recurrent neural networks
Palmé et al. Early warning of gas turbine failure by nonlinear feature extraction using an auto-associative neural network approach
US11339763B2 (en) Method for windmill farm monitoring
KR102350636B1 (ko) 플랜트 고장 예지 장치 및 방법
Aksan et al. Review of the application of deep learning for fault detection in wind turbine
KR102212022B1 (ko) 양수 수차의 건전성 자동 판정 방법 및 이를 위한 시스템
KR102350635B1 (ko) 플랜트 고장 예지 장치 및 방법
Qosim et al. Fault Detection System Using Machine Learning on Synthesis Loop Ammonia Plant
WO2006107295A1 (en) Model-based detection, diagnosis of turbine engine faults
KR20210101914A (ko) 고장 예지 모델 생성 장치 및 방법
MOHAMED development and application of an importance measures for predictive maintenance in oil and gas installations
Loyer Empirical Evaluation of the Impact of Data Pre-Processing on the Performance of Predictive SHM of Jet Engines

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right