基于統(tǒng)計過程控制的QAR數(shù)據(jù)風(fēng)險預(yù)測方法

孫瑞山，張思遠，汪 磊

（中國民航大學(xué) 民航安全科學(xué)研究所，天津 300300）

0 引言

飛行記錄數(shù)據(jù)（FDR/QAR/ACARS）包括大量飛行、操作和環(huán)境信息，目前主要應(yīng)用于飛機故障檢測和簡單的運行超限管理，大量數(shù)據(jù)的利用缺乏系統(tǒng)性，沒有得到有效、系統(tǒng)地利用，造成信息浪費。

在飛行品質(zhì)監(jiān)控方面，民航業(yè)對飛行數(shù)據(jù)的使用大多數(shù)還是以基于觸發(fā)邏輯的“超限管理”為主，依靠FDR和QAR數(shù)據(jù)中出現(xiàn)的超限事件的等級與次數(shù)來監(jiān)控航班及飛行員，根據(jù)各個超限事件的“軟、硬”監(jiān)控標(biāo)準(zhǔn)來考核飛行安全。

傳統(tǒng)的“超限管理”存在以下不足：

（1）由于超限監(jiān)測軟件只是監(jiān)測觸發(fā)邏輯表達式的事件，因此大量的沒有觸及限值的飛行數(shù)據(jù)信息被浪費，特別是處在超限范圍內(nèi)接近觸發(fā)而沒有觸發(fā)超限事件的信息，因為被過濾掉而沒有得到科學(xué)合理的整合、分析和運用。

（2）缺乏對噪音的管理。影響飛行品質(zhì)的因素可以大致分為兩類因素：偶然因素和異常因素。偶然因素是始終存在的并且不能被徹底消除的，例如飛機的起飛滑跑路徑不可能是一條直線；飛機飛行的航跡不可能與航線嚴(yán)絲合縫等；異常因素時而存在，對飛行品質(zhì)影響很大，但是異常因素是可以被消除的，例如由于飛行員過于疲勞導(dǎo)致飛機飛離飛行包線等。

倘若在飛行中由于較大的偶然因素和較小的異常因素引起了超限，那么在傳統(tǒng)的超限管理中責(zé)任的界定將是困難的。

（3）虛發(fā)警報和漏發(fā)警報。在實際飛行中，單純地由偶然因素引起超限的概率雖然很小，但不是絕對不可能發(fā)生的。倘若發(fā)生，傳統(tǒng)的超限管理根據(jù)“數(shù)據(jù)超限即告警”的原則發(fā)出“虛假的”警報。這種錯誤將會造成相關(guān)人員尋找并不存在的異常因素，對航空安全管理資源造成了浪費。如圖1，數(shù)據(jù)在很長的一段時間內(nèi)都處于正常，突然有一個點出現(xiàn)異常，之后又迅速回歸正常。這種異常很有可能是由于偶然因素引起的。

圖1 虛發(fā)警報

同樣地，在實際飛行中，盡管出現(xiàn)了異常因素，但是尚未引起超限。在這種情況中，傳統(tǒng)的超限管理就會漏發(fā)警報。如圖2，假設(shè)數(shù)據(jù)大于等于1.6為超限，那么所有的數(shù)據(jù)均未超限，在當(dāng)前的判斷邏輯下不會發(fā)出告警，但此時飛機已經(jīng)處于危險的邊緣，應(yīng)當(dāng)發(fā)出告警，此種情況稱為“漏發(fā)警報”[1]。

圖2 漏發(fā)警報

1 改進的QAR數(shù)據(jù)分析方法

過程控制的概念與實施過程監(jiān)控的方法在20世紀(jì)20年代就由美國的休哈特（W.A.Shewhart）提出，它能科學(xué)地區(qū)分出過程中的偶然波動與異常波動，從而對過程的異常及時告警?？刂茍D是對過程質(zhì)量特性值進行測定、記錄、評估和監(jiān)察過程是否處于控制狀態(tài)的一種用統(tǒng)計方法設(shè)計的圖。

傳統(tǒng)的QAR數(shù)據(jù)超限分析沒有對飛行中的偶然因素和異常因素進行有效的界定。為了區(qū)分異常因素，我們首先可以假設(shè)過程中只存在偶然因素，不存在異常因素或者異常因素已經(jīng)被消除。根據(jù)當(dāng)前圖中的最小波動，應(yīng)用統(tǒng)計學(xué)原理我們可以設(shè)計出控制圖的控制限[2]。

1.1 QAR數(shù)據(jù)風(fēng)險告警的關(guān)口前移

實現(xiàn)QAR數(shù)據(jù)風(fēng)險告警的關(guān)口前移就必須對若干個特征值進行系統(tǒng)的分析，盡可能多地發(fā)現(xiàn)尚未引起超限的異常因素，如下圖3：

圖3 告警關(guān)口前移

這樣就做到了風(fēng)險告警關(guān)口前移，起到了預(yù)防作用。

在某些情況下，沒有任何先兆，控制圖會突然顯示異常，這時需要查出異常因素，采取相關(guān)措施，保證消除異常因素并不再出現(xiàn)，并將其納入標(biāo)準(zhǔn)。每進行一次這樣循環(huán)，就會消除一個影響飛行品質(zhì)的異常因素，令其永不再出現(xiàn)，也起到了預(yù)防作用[3～5]。

1.2 控制限的確定

1.2.1 ±3σ原理

我們將虛發(fā)警報這類錯誤稱為第I類錯誤，發(fā)生第I類錯誤的概率記為α；將漏發(fā)警報的錯誤稱為第II類錯誤，發(fā)生第II類錯誤的概率記為β。

孤立地看，這兩類錯誤都可以減少甚至消除，但是同時避免兩類錯誤是不可能的。如果放寬控制界限可以減少第I類錯誤例如將控制界限從μ±3σ放寬至μ±5σ，則

此時α=0.0001%，即每百萬次約有一次第I類錯誤，但是由于控制界限的放寬，使得第II類錯誤更容易出現(xiàn)，即β增大。如果壓縮控制界限，可以減小第II類錯誤的發(fā)生概率β，但會增加第I類錯誤的發(fā)生概率α。經(jīng)驗證明，休哈特所提出的±3σ控制限在很多情況下都接近最優(yōu)。這個原理也被稱為±3σ原理。

需要特別指出的是，±3σ控制限是一個用來參考的控制限，在實際應(yīng)用中，由于應(yīng)用背景的差異，如果我們認(rèn)為α過大而β過小，我們可以根據(jù)分析需要隨時更改控制限[6～8]。

1.2.2 改進后的風(fēng)險預(yù)測告警準(zhǔn)則

若統(tǒng)計量服從正態(tài)分布，根據(jù)±3σ原理，α=1-p(x-μ)≤0.27%，這是一個相當(dāng)小的概率值，我們幾乎可以肯定，只要有一個點超限就意味著飛行中出現(xiàn)了異常。但正如上文中的分析，α很小就意味著β很大，如果我們只看到一個點未超限是遠遠不能判斷整個飛行過程是合乎規(guī)范的。若取得連續(xù)的若干個飛行數(shù)據(jù)，設(shè)為n個，根據(jù)乘法原理，總的βn會比單個點的β小很多。由此，我們可以得到結(jié)論：在控制限內(nèi)連續(xù)的若干點，即使有個別的點超限，整個飛行過程仍然可以被看做是合格的。

由以上結(jié)論我們得到如下三個判斷飛行過程合格的情況：

（1）連續(xù)25個點都未超過控制限；

（2）連續(xù)35個點最多有1個點超過控制限；

（3）連續(xù)100個點最多有2個點超過控制限；

這是一個與同數(shù)量級的小概率。同理，計算得到α1=0.0654，α3=0.0026 。

有了三種判斷飛行過程合格的判斷準(zhǔn)則，同理，我們不難得到判斷飛行過程異常的準(zhǔn)則。ISO8258：1991即國標(biāo)GB/T4091-2001引用了8種判斷異常的準(zhǔn)則，在沒有積累相關(guān)QAR數(shù)據(jù)分析經(jīng)驗的情況下，我們不妨直接引用，待積累足夠的相關(guān)經(jīng)驗后可以做出適應(yīng)飛行品質(zhì)管理的相應(yīng)的調(diào)整。

如下是8種判斷飛行過程出現(xiàn)異常因素的判斷準(zhǔn)則：

（1）一點落在控制限外；

（2）連續(xù)9點落在控制限一側(cè)；

（3）連續(xù)6點遞增或遞減；

（4）連續(xù)14點中相鄰的點上下交替；

（5）連續(xù)3點中有2點落在中心線同一側(cè)的2σ區(qū)以外。如圖5。

圖4 準(zhǔn)則5圖示

（6）連續(xù)5點中有4點落在中心線同一側(cè)的1σ區(qū)以外，如圖5。

圖5 準(zhǔn)則6圖示

（7）連續(xù)15點在1σ區(qū)中心線上下。

（8）連續(xù)8點在中心線兩側(cè)，且無一個在1σ區(qū)中。

需要特別指出的是，控制限和超限警告限不是一個概念，超限警告限是一種規(guī)范，而控制限是用來區(qū)分飛行中的偶然因素和異常因素的，二者不能混淆[9,10]。

2 案例演示

2.1 QAR原始數(shù)據(jù)

目前飛機重著陸的診斷判定主要是依賴QAR數(shù)據(jù)，關(guān)鍵指標(biāo)是著陸接地瞬間主起落架的垂直加速度，在空地電門轉(zhuǎn)換瞬間，機載QAR設(shè)備可有效記錄飛機垂直加速度值的變化VERT ACCEL為Vertical Acceleration的縮寫，代表飛機垂直加速度，QAR原始數(shù)據(jù)表中縮寫為VRTG，對應(yīng)各數(shù)值意為重力加速度g的倍數(shù)。限于篇幅原始數(shù)據(jù)表略。

2.2 建立控制圖

分析“飛機重著陸風(fēng)險”，依據(jù)慣例，截取QAR數(shù)據(jù)中飛機著陸前2秒和著陸后5秒的數(shù)據(jù)進行分析，截取的數(shù)據(jù)如下表1。

表1 截取后的QAR數(shù)據(jù)

QAR以1/16秒為間隔記錄飛機垂直加速度的數(shù)據(jù)，按照時間順序?qū)?shù)據(jù)輸入MINITAB軟件運行得到單值-移動極差控制圖。

圖6 QAR數(shù)據(jù)均值-移動極差控制圖

依據(jù)判斷準(zhǔn)則，可以得到QAR著陸數(shù)據(jù)的單值控制圖檢驗結(jié)果如表2所示。

控制圖共發(fā)出8個告警，他們可能由于異常因素引起的。事實上，對于這組未超限的QAR數(shù)據(jù)來說，這次檢驗的虛發(fā)警報太多了，故我們可以放寬控制限，也可以減少判斷異常的準(zhǔn)則，例如，我們不妨將控制限放寬到±4σ，將判斷為異常的準(zhǔn)則只留下“1.一點落在控制限外”。重新得到分析結(jié)果如下圖。

圖7 修正后的QAR數(shù)據(jù)均值-移動極差控制圖

表2 QAR著陸數(shù)據(jù)分析結(jié)果

檢驗結(jié)果如表3所示。

表3 修正后的QAR著陸數(shù)據(jù)分析結(jié)果

取兩個圖告警結(jié)果的交集，得到點43已經(jīng)非常接近超限，查閱原數(shù)據(jù)得知點43為1.441，距離B738著陸垂直載荷輕度超限1.68已經(jīng)非常接近，超限風(fēng)險較大，應(yīng)當(dāng)提前引起注意。

3 結(jié)束語

本文中對于偶然因素和一場因素的判斷準(zhǔn)則引用了ISO8258：1991的規(guī)定，實際上，該判斷準(zhǔn)則在生產(chǎn)領(lǐng)域應(yīng)用廣泛，在民航業(yè)的數(shù)據(jù)分析中適用性可能降低。

在運用本文方法進行QAR數(shù)據(jù)分析時需要根據(jù)經(jīng)驗重新設(shè)定判斷準(zhǔn)則，但控制限的設(shè)定原則是不變的，應(yīng)當(dāng)基于ISO8258：1991推薦的判斷準(zhǔn)則，在此基礎(chǔ)上科學(xué)合理地增加或者刪除。

[1]張公緒,孫靜.新編質(zhì)量管理學(xué)[M].北京:高等教育出版社，2003.

[2]Vender Wiel S A.Algorithmic Statistical Process Control:Concepts and Application[J].Rechnometrics,1992,34(3).

[3]美國質(zhì)量協(xié)會.注冊質(zhì)量工程師手冊克勞士比中國學(xué)院譯[M].北京:中國城市出版社,2003.

[4]孫瑞山,張思遠.人的可靠性綜合分析模式及應(yīng)用[J].工業(yè)工程，2014,17(2).

[5]汪磊,孫瑞山,吳昌旭,崔振新,陸正.基于飛行QAR數(shù)據(jù)的重著陸風(fēng)險定量評價模型[J].中國安全科學(xué)學(xué)報,2014,24(2).

[6]趙會珍,張聚濤.理解SPC控制圖的關(guān)鍵——小概率原理[J].統(tǒng)計與決策,2009,（22）.

[7]康軍.關(guān)于統(tǒng)計過程控制中樣本組數(shù)量的研究.統(tǒng)計與決策,2012,（19）.

[8]George-Box,Luceno A.Discrete Proportional-IntegralAdjustment andStatistical Process Control[J].Journal of Quality Technology,1997,29(3).

[9]Sun R S,Chen Y F,Liu XY,et al.A Method of Analysis Integrating HCR and ETA Modeling for Determining Risks Associated With Inadequate Flight Separation Events[J].Journal of Aviation Technology and Engineering,2011,1(1).

[10]Wang L,Wu C X,Sun R S.An Analysis of Flight Quick Access Recorder(QAR)Data and Its Applications in Preventing Landing Incidents[J].Reliability Engineering&System Safety,2014,(127).