李鵬

摘 要 空域容量評估是空中交通安全、高效實施的關(guān)鍵，其中管制員的工作負荷是空域容量大小的主要影響因素。本文以塔臺管制員語音負荷為研究對象，利用MATLAB對管制話音進行分析，研究了端點檢測方法對管制員話音通話數(shù)據(jù)的獲取，提取出特征參數(shù)，利用SPSS軟件進行因子分析，利用提取的特征參數(shù)表征出管制話音負荷大小。最后應(yīng)用收集到的某塔臺管制通話為例，求出多段通話的語音負荷，統(tǒng)計出對應(yīng)航空器架次，擬合出語音負荷與航空器架次之間的函數(shù)關(guān)系，根據(jù)國際民航組織9426號文件中的規(guī)定進行了容量評估，為空域容量評估提供技術(shù)支撐。

關(guān)鍵詞 空中交通管理 容量評估 管制負荷 語音分析

中圖分類號：X951文獻標識碼：A

0引言

管制部門是保證民航安全高效運行的基礎(chǔ)。為保證航空安全績效水平，對管制員的工作量和工作質(zhì)量進行有效的評估具有重要意義，通過對管制負荷的有效評估，進而可以得出對應(yīng)扇區(qū)的容量值，當扇區(qū)容量過大時，對一定時間內(nèi)扇區(qū)的航空器數(shù)量進行有效的控制，從而確保民航在空管這一環(huán)的安全性和高效性。

國內(nèi)外用于空域容量負荷評估的方法較多，包括有建立數(shù)學模型、快速模擬、動態(tài)模擬等等。目前最常用管制負荷評估方法包括英國研究的DORATASK法，DORATASK將管制員的工作負荷分為語音通話負荷，操作負荷和思考負荷，將這幾個負荷相加即得管制負荷。斯坦福研究中心研發(fā)的相對容量評估過程（The Relative Capacity Estimating Process，RECEP），它將管制員所有可以考慮的管制行為都進行量化，包括通話時間、飛行進程單時間、集中于雷達屏幕上監(jiān)視的時間以及沖突解決的時間，從而得出一個小時內(nèi)總的工作時間，RECEP模型將空域容量定義為每小時內(nèi)空域內(nèi)的航空器架次，并且認為管制負荷的極限值極大決定了空域容量的大小，他們研究得出工作時間與每小時的航空器架次之間呈平方關(guān)系。2006年，美國的Arthur P.Smith，Anand D.Mundra等人對加裝ADS-B（Automatic Dependent Surveillance-Broadcast，廣播式自動相關(guān)監(jiān)視）設(shè)備對管制員工作負荷的影響。Peng Wei，Chittayong Surakitbanbarn等人利用動態(tài)密度對管制負荷進行評估，動態(tài)密度是指影響空中交通管制復(fù)雜度和難度的各種因素或者說各種變量。本文主要研究通過對管制員的語音負荷進行定量評估，利用得到的管制員負荷進一步得出該空域的容量，得到更加科學的空域容量評估后，有利于空域的安全和管制員的身心健康。

1管制語音信號分析

語音信號分析，是為了從收集到的語音信號中提取出研究所需要的特征參數(shù)，從而可以找出各變量與指標之間的關(guān)系，進行后續(xù)的研究。

1.1語音信號預(yù)處理

首先要對原始語音信號進行讀入，把音頻信號轉(zhuǎn)化為計算機可以處理的數(shù)字信號，在MATLAB中通常用wavread函數(shù)完成，得到的數(shù)字化音頻信號通常存貯在列向量中。對語音信號的預(yù)處理，通常包含兩步：分幀和加窗。

語音信號進行分幀處理就是基于語音信號的“準穩(wěn)態(tài)”特征，將其在分割成極短時間的多個語音片段，在穩(wěn)定狀態(tài)下逐幀提取出特征參數(shù)。同時為了使各幀之間可以平穩(wěn)過渡，防止由于超過能保持穩(wěn)定狀態(tài)的極短時間，導(dǎo)致參數(shù)發(fā)生變化，通常分幀技術(shù)會引入幀移的概念（），使相鄰兩幀之間部分重疊，幀移圖如圖1所示。

幀移是指后一幀對前一幀的位移量（inc），而相鄰兩幀之間的重疊部分overlap=wlen-inc。對于長為N的語音信號，分幀公式為：

其中，N為信號長度，wlen為幀長，inc為幀移，overlap為相鄰幀重疊部分，該長為N的語音信號將被分為fn幀。本文中，分幀的實現(xiàn)主要是：首先調(diào)用wavread函數(shù)，將格式為.wav的音頻文件讀入轉(zhuǎn)化為數(shù)字信號。分幀主要利用enframe函數(shù)，每一幀的數(shù)據(jù)順序儲存在數(shù)組的每一列中，列數(shù)即為幀數(shù)。加窗則是加在設(shè)定的幀長wlen 上，hamming（wlen）。

1.2端點檢測法

特征參數(shù)提取主要運用語音信號單參數(shù)雙門限端點檢測法檢測出是否有話音部分，因此，需要提前設(shè)定區(qū)分背景噪音和實際話音的閾值，以及判斷靜默段是否為短暫停頓以及話音段長度是否足以構(gòu)成一句話的參數(shù)。

閾值設(shè)定：根據(jù)前導(dǎo)無話段的短時平均能量也就是背景噪音的能量設(shè)置出兩個門限值，一個較大為T2，一個較小為T1。噪聲的估計通常依據(jù)前導(dǎo)無話段來實現(xiàn)。根據(jù)波形圖估計出前導(dǎo)無話段的時長IS（單位：s），則前導(dǎo)無話段的幀數(shù)NIS為公式2所示：

其中，IS為前導(dǎo)無話段時長。由于管制員工作環(huán)境的特殊性，本文將其閾值看為定值。本文閾值的確定是先取前導(dǎo)無話段IS長度為0.25秒，根據(jù)式（2）得出前導(dǎo)無話段的幀數(shù)NIS，對無話段逐幀求短時平均能量，存入矩陣eis中，顯示eis中各幀的能量，根據(jù)這些能量值取T1，T2。本文選取塔臺通話錄音的前導(dǎo)無話段短時平均能量大小設(shè)定T1和T2分別為0.0001和0.05。然后對語音信號特征參數(shù)進行提?。ㄌ崛〉恼Z音信號特征參數(shù)為短時能量、短時平均過零率、通話次數(shù)和時長以及飽和度）。

1.3語言信號分析

分幀后的語音信號，按列存儲在數(shù)組中，幀數(shù)為，按照求過零率的思路，求每幀的過零率。首先對分幀后的矩陣逐列調(diào)用，讀取每列的數(shù)據(jù)，即得到每幀的所有采樣點信息，按幀長循環(huán)，判斷相鄰取樣點之間的符號變化情況，當相鄰取樣點乘積為負，則說明信號過零一次，計數(shù)加一，繼續(xù)循環(huán)，一幀循環(huán)完后將該幀過零情況存入行向量中。最終求和可得該段語音信號的過零率。對數(shù)字化音頻信號求列向量長度，得到采樣點數(shù)N，結(jié)合采樣頻率可求得信號的時間坐標time，在繪制語音波形圖和后續(xù)求飽和度時會用上。對分幀后的信號逐幀點平方求和，分別得到每幀的短時平均能量etemp，求和可得語音信號能量setemp。

最后基于短時平均能量的單參數(shù)雙門限端點檢測法，檢測出該段語音中的通話次數(shù)并得到voiceseg結(jié)構(gòu)數(shù)據(jù)?？傻玫酵ㄔ捳Z音次數(shù)，通話時長和飽和度。流程圖2為本文語音信號分析程序?qū)崿F(xiàn)的整體思路：（見圖2）。

2管制語音特征參數(shù)分析

本文選擇使用SPSS軟件進行因子分析，找出各參數(shù)之間的內(nèi)在聯(lián)系，最終實現(xiàn)用語音參數(shù)表示出管制負荷的大小，用少數(shù)幾個因子去描述原資料的大部分信息。

2.1因子分析

因子分析是將有內(nèi)在聯(lián)系的幾個變量，找到支撐它們相互關(guān)系的公共因子，使變量能用公共因子和特殊因子這兩部分來表示的多變量統(tǒng)計分析方法。本文的思路是將語音分析得到的五項特征參數(shù)（通話語音的次數(shù)、能量、時長、飽和度、過零率）用SPSS進行因子分析，設(shè)它們分別為1，2，3，4，5作為五個隨機變量。

關(guān)系式3就是因子模型，表示各變量之間的關(guān)系，fi為公共因子，ei為特殊因子，bij（i=1，2，…，m）為因子的載荷，表示變量xi對公共因子fi的影響。其中，公共因子和特殊因子各自內(nèi)部之間以及兩者之間都是互不相關(guān)，相互獨立的。特殊因子與各自對應(yīng)的變量有關(guān)，用于補齊公共因子無法完全表示出來的部分。

2.2語音特征參數(shù)因子分析運用

提取9段不同繁忙時刻的管制錄音，每段3分鐘，分別提取特征參數(shù)，進行因子分析。在SPSS軟件上的操作步驟具體如下：

（1）數(shù)據(jù)的輸入，在軟件中分別輸入9段管制語音的特征參數(shù);

（2）按順序點擊分析→降維→因子分析→選擇要分析的變量→在“描述”中選擇KMO和Bartlett球度檢驗→在“抽取”中選擇“主要成分提取法”→在“旋轉(zhuǎn)”中選擇“最大方差法（varimax）”→開始分析;

（3）得到分析結(jié)果。

分析結(jié)果如表1，2，3所示：

根據(jù)結(jié)構(gòu)效度分析結(jié)果看出，KMO值大于0.7，根據(jù)Kaiser的觀點，屬于中等，KMO值越大表明越適合進行因子分析。

結(jié)合式3中描述的因子模型，可以因子載荷矩陣表示出因子模型，如式3-2所示：

其中，為公共因子，～為特殊因子。利用主要成分提取法抽取出因子得分系數(shù)矩陣，將公共因子表示為各變量的線性組合。

其中1～5 ，為選取的變量。語音信號的通話次數(shù)、時長等特征的變化是受語音負荷的變化而相應(yīng)變化影響的，因此語音負荷就是這五項特征參數(shù)的公共因子。語音負荷因此依據(jù)表3可得語音負荷的表達式為：

表4為提取的不同通話錄音提取參數(shù)，并計算負荷所得結(jié)果。

由表4可以看出，不同特征參數(shù)大小，語音負荷的大小有明顯差異。

3空域容量評估

根據(jù)《空中交通服務(wù)規(guī)劃手冊》（DOC.9426）DORATASK法所指出的：當航班量達到該扇區(qū)的容量值時，管制員的平均工作負荷不得超過工作負荷峰值的80%，并且工作負荷超過最大值的90%的時間不得超過總工作時間的2.5%。那么說明當負荷達到峰值的80%時，此時扇區(qū)內(nèi)的航班量即為該扇區(qū)的容量。

3.1最小二乘法回歸分析

為了得到負荷峰值80%時的航空器數(shù)量，需要得到管制負荷與航空器數(shù)量之間的關(guān)系式。因此，選用最小二乘法回歸分析，將多個管制負荷與航空器數(shù)量一一對應(yīng)的散點擬合成方程。如式7所示：

其中的a，b使F（a，b）最小，即函數(shù)F分別對a，b求偏導(dǎo)，令偏導(dǎo)等于零，得到方程組可求出a，b，得出的該條擬合函數(shù)就是最接近散點圖的直線，得到最小二乘法擬合的函數(shù)。

3.2基于話音負荷的容量評估

本文空域容量評估主要的依據(jù)是國際民航組織《空中交通服務(wù)規(guī)劃手冊》（DOC.9426）中的DORATASK法所指出的：當航班量達到該扇區(qū)的容量值時，管制員的工作負荷不得超過負荷最大值的80%，并且工作負荷超過最大值的90%的時間不得超過總工作時間的2.5%。因此可以基本確定空域容量的評估方法。

基于話音負荷的空域容量評估過程基本如下：

計算出多個時間片管制員的話音負荷;統(tǒng)計每個時間片對應(yīng)的航空器架次;根據(jù)第（1）、（2）步的結(jié)果，做出散點圖，判斷負荷與航空器架次的大致關(guān)系;利用MATLAB進行函數(shù)擬合，得到話音負荷與航空器架次之間的對應(yīng)函數(shù)關(guān)系;由于目前大多使用雷達管制，管制員的通話相比程序管制有所減少，因此當平均工作負荷達到峰值的70%時對應(yīng)的航空器架次即是空域容量。

根據(jù)上述原理，收集某塔臺管制單位話音數(shù)據(jù)，首先計算管制員的22個10分鐘時間片的話音負荷，并統(tǒng)計相應(yīng)航空器架次，得到飛機架次與話音負荷一一對應(yīng)的多個散點，對應(yīng)關(guān)系如表5所示：

利用MATLAB的polyfit函數(shù)對散點圖進行回歸分析，由此得到擬合結(jié)果：

其中，N為航空器架次。

由于本文選取的是塔臺席通話錄音，因此仍以峰值的80%進行計算，即工作負荷達到峰值的80%時，此時航空器的流量N為扇區(qū)的容量值。根據(jù)函數(shù)可以求得管制負荷的峰值和對應(yīng)的航空器數(shù)量。可以得到采集語音數(shù)據(jù)的該塔臺管制單位塔臺席的語音負荷峰值Wmax為810107.6871，那么峰值80%時對應(yīng)的航空器架數(shù)約為3架，說明該管制單位塔臺席較合理的安排是10分鐘的時間片內(nèi)同時掌控的航空器架次不超過3架。

4結(jié)束語

本文從管制話音特征參數(shù)入手，對一線通話錄音進行語音分析，提取出特征參數(shù)。為了使特征參數(shù)具體量化出語音負荷的大小，利用因子分析的方法對數(shù)據(jù)進行分析，求出本文所確立的語音負荷。最后利用MATLAB進行擬合可得到語音負荷與航空器架次之間關(guān)系的方程式，根據(jù)國際民航組織9426號文件要求，確定對應(yīng)的航空器架次即為空域容量，能夠有效提高空域容量評估的效率。

參考文獻

[1] Tuan，P.L.&H.S.Procter，etal.Advanced Productivity Analysis Method for Air Traffic Control Operations[J]. FAA-RD-76-164， 1976.

[2] Smith，A.P.&A.D.Mundra.Impact of ADS-B on Controller Workload： Results from Alaska's Capstone Program[C]. ieee/aiaa 25TH Digital Avionics Systems Conference，Portland，OR，USA，2006.

[3] Peng Wei， C.Surakitbanbarn，et al. Workload evaluation of sectorized air traffic control and stream management[C].? Integrated Communications， Navigation and Surveillance Conference，Herndon，VA，USA，2012.

[4] 宋知用.MATLAB在語音信號分析與合成中的應(yīng)用[M].北京：北京航空航天大學出版社，2013.

[5] 榮立元.運用SPSS軟件因子分析學生綜合成績的管理[J].山東紡織經(jīng)濟，2014（02）.

[6] Arnab，M.&P.John.Estimating the capacity of europes airspace using a simulation model of air traffic controller workload[C].Proceedings of Transportation Research Board 80th Annual Meeting，2001：30-43.

[7] International Civil Aviation Organization. Manual air transport service plan[S]. No.9426-AN/924-1999 Montreal：ICAO，1999.