基于工業(yè)攪拌器選型的層次分析法權(quán)重計算方法①

吳長松 吳夢楚 梁佳赟 雷雨田 楊 斌

(西北大學(xué)化工學(xué)院)

基于工業(yè)攪拌器選型的層次分析法權(quán)重計算方法①

吳長松 吳夢楚 梁佳赟 雷雨田 楊 斌

(西北大學(xué)化工學(xué)院)

為了應(yīng)用層次分析法進(jìn)行工業(yè)攪拌器選型，采用兩種常用標(biāo)度針對攪拌器選型實例建立比較矩陣，用4種矩陣計算方法求解權(quán)重向量，最終得到了綜合權(quán)重排序，并確定了每種計算方法求解矩陣的計算機耗時。綜合權(quán)重排序結(jié)果顯示90/9～99/9標(biāo)度法相比于1～9標(biāo)度法的誤差較小，較為符合實際，更適于工業(yè)攪拌器的層次分析法選型。在90/9～99/9標(biāo)度下比較綜合權(quán)重結(jié)果與耗時量可知，方根法(RMS)相比其他3種矩陣計算方法更適于攪拌器選型的層次分析法應(yīng)用，在同樣滿足一致性和計算精度的基礎(chǔ)上，其計算時間最少(效率最高)，耗時分別是和積法(SPM)的2.50%、特征根法(EM)的32.36%、行和歸一化法(NRA)的76.34%。

攪拌器選型 層次分析法 標(biāo)度 權(quán)重向量

選擇合適的攪拌器是攪拌設(shè)備設(shè)計過程中相當(dāng)關(guān)鍵的一步，攪拌器的型式對攪拌效果以及攪拌設(shè)備制造和操作成本有很大影響[1，2]。長期以來，攪拌器的設(shè)計參數(shù)繁多復(fù)雜且不一致，設(shè)計過程缺乏系統(tǒng)成熟的理論，選型過程過多依靠設(shè)計人員的專業(yè)知識和經(jīng)驗。近年來，人們開始應(yīng)用層次分析法來改善攪拌器選型效率，提高選型合理性[3，4]。層次分析法是Saaty T L提出的一種多準(zhǔn)則決策方法[5]，它將復(fù)雜的決策問題分解成多個組成因素，按目標(biāo)層、準(zhǔn)則層及方案層等多個層次排列，按照一定標(biāo)度將人的主觀判斷進(jìn)行客觀的量化，利用量化信息將決策的方式數(shù)量化，為復(fù)雜的決策問題提供一種具有實用性和有效性的決策方法。運用層次分析法進(jìn)行決策時，關(guān)鍵在于單一比較矩陣權(quán)重向量和綜合權(quán)重系數(shù)的求解。筆者基于工業(yè)攪拌器選型實例，采用1～9標(biāo)度和90/9～99/9標(biāo)度建立比較矩陣，運用特征根法(EM)、方根法(RMS)、和積法(SPM)、行和歸一化法(NRA)4種常用的算法求解矩陣的權(quán)重向量，得到總的權(quán)重系數(shù)[6，7]，根據(jù)結(jié)果比較得到更為適宜工業(yè)攪拌器選型的標(biāo)度和計算方法。

1 標(biāo)度方法介紹

1.1 1～9標(biāo)度

根據(jù)遞階層次模型，將相同層次的因素之間進(jìn)行兩兩比較。為了數(shù)量化判斷思維來簡化問題分析過程，層次分析法定義了標(biāo)度的概念，表1給出了由Satty T L最早提出的1～9標(biāo)度方法，因其矩陣元素簡單，可以直觀地看出元素之間的對比關(guān)系，因而使用較為廣泛。

表1 1～9標(biāo)度及其含義因素關(guān)系對比

1.2 90/9～99/9標(biāo)度

90/9～99/9標(biāo)度對標(biāo)度元素進(jìn)行了完善，增大了標(biāo)度的區(qū)分度，提高了計算精度，使之更符合實際應(yīng)用(表2)[8]。

表2 90/9～99/9標(biāo)度及其含義因素關(guān)系對比

1.3 比值與語言程度的對應(yīng)關(guān)系

標(biāo)度選定后，分析各準(zhǔn)則層對目標(biāo)層以及各方案層對準(zhǔn)則層的影響，確定兩兩比較矩陣的比值關(guān)系，用比值去對應(yīng)語言程度，然后通過語言程度去對應(yīng)標(biāo)度值，從而建立比較矩陣(表3)[9]。

表3 比值與語言程度的對應(yīng)關(guān)系

2 權(quán)重向量W的計算方法

假設(shè)比較矩陣A=(aij)n×n(i，j=1，2，…，n)，比較矩陣A的權(quán)重向量W=(ω1，ω2，…，ωn)T，比較矩陣的最大特征值為λmax。

特征值法(EM)：

AW=λmaxW

(1)

在特征值計算中，若A為比較矩陣且滿足一致性要求，則向量W為矩陣A的特征向量。在矩陣左側(cè)乘以λmax，運用線性代數(shù)的方法求解特征值，求得的向量W即為特征向量，對特征向量W進(jìn)行歸一化處理得到權(quán)重向量W=(ω1，ω2，…，ωn)T。

方根法(RMS)：

(2)

和積法(SPM)[10]：

(3)

行和歸一化(NRA)：

(4)

最大特征值λmax的計算公式為：

(5)

其中，(AW)i表示向量AW中的第i個因素。

3 實例計算與結(jié)果分析

3.1 建立比較矩陣

工業(yè)攪拌器實例來源于Scargiali F等的生物攪拌器研究[11]，表4為該文獻(xiàn)中列出的直葉圓盤渦輪攪拌器、開啟渦輪攪拌器、A310攪拌器和螺旋槳攪拌器4種攪拌器要達(dá)到自由表面渦流的臨界轉(zhuǎn)速Nc、臨界單位體積功率P/Vc、臨界功率準(zhǔn)數(shù)NPc、臨界單位體積傳質(zhì)系數(shù)KLac。

表4 攪拌器參數(shù)

根據(jù)層次結(jié)構(gòu)模型和攪拌器參數(shù)，分別用1～9標(biāo)度和90/9～99/9標(biāo)度建立準(zhǔn)則層對目標(biāo)層比較矩陣A，方案層對準(zhǔn)則層比較矩陣B1、B2、B3、B4[12]。

1～9標(biāo)度：

90/9～99/9標(biāo)度：

3.2 求解綜合權(quán)重

基于上述建立的比較矩陣，運用EM法、RMS法、SPM法和NRA法解矩陣，得到總的權(quán)重系數(shù)[13，14]見表5、6。

表5 攪拌器選擇方案的總排序(1～9標(biāo)度)

表6 攪拌器選擇方案的總排序(90/9～99/9標(biāo)度)

由表5、6可知，采用1～9標(biāo)度和90/9～99/9標(biāo)度計算得到的排序結(jié)果皆為A310型、斜葉開啟渦輪、直葉圓盤渦輪和螺旋槳型攪拌器。結(jié)果表明A310型攪拌器效果最好，綜合性能優(yōu)于其余3種攪拌器。在1～9標(biāo)度下，由于比值與標(biāo)度之間數(shù)值差距較大，造成不同方案的綜合權(quán)重系數(shù)結(jié)果彼此偏差較大，存在較大的誤差[15]。90/9～99/9標(biāo)度比1～9標(biāo)度比值與標(biāo)度對應(yīng)效果較好，最終結(jié)果更符合實際。因此，應(yīng)用層次分析法進(jìn)行工業(yè)攪拌器選型可選擇90/9～99/9標(biāo)度建立比較矩陣。同樣由表5、6可知，4種計算方法求解權(quán)重向量的結(jié)果極為接近，精度基本一致。

3.3 分析計算耗時

基于MATLAB編寫程序分別計算出在兩種標(biāo)度下用4種算法求解100個單一隨機比較矩陣時所耗時間，得到求解矩陣的計算耗時，結(jié)果見表7、8。

表7 不同計算方法求解矩陣耗時(1～9標(biāo)度)

表8 不同計算方法求解矩陣耗時(90/9～99/9標(biāo)度)

由表7、8可知， RMS法求解矩陣時耗時最小，在90/9～99/9標(biāo)度下所用時間僅為SPM法的2.50%、EM法的32.36%、NRA法的76.34%。即使采用1～9標(biāo)度，RMS法耗時也僅為NRA法的88.00%。在層次分析法選型過程中，攪拌器備選方案越多，RMS法計算效率的優(yōu)勢越明顯。

4 結(jié)束語

針對實際工業(yè)攪拌器選型的層次分析法應(yīng)用，所選取的兩個標(biāo)度中，90/9～99/9標(biāo)度更符合實際，更適于層次分析法的攪拌器選型。在90/9～99/9標(biāo)度下比較綜合權(quán)重結(jié)果和耗時量可知，RMS法在同樣符合計算精度和排序一致性的前提下，相比其余3種矩陣計算方法耗時更少，更適于攪拌器選型的層次分析法應(yīng)用。綜上，在工業(yè)攪拌器選型中選用RMS法求解矩陣的權(quán)重向量較為適宜。

[1] Tatterson G B，Kyser E A.Mixing in Fed-Batch and Continuous Flow Processes in Nonstandard Geometries[J].AIChE Journal,1991,37(2):269～273.

[2] 宮磊,周麗,崔文璟,等.攪拌槳組合數(shù)值模擬優(yōu)化及在谷氨酰胺轉(zhuǎn)胺酶發(fā)酵中的應(yīng)用[J].化工進(jìn)展,2014,33(10):2570～2575.

[3] 代麗,朱愛華,趙勻.應(yīng)用層次分析法計算分插機構(gòu)優(yōu)化目標(biāo)的權(quán)重[J].農(nóng)業(yè)工程學(xué)報,2013,29(2):60～65.

[4] 王嘉駿,馮連芳,顧雪萍,等.層次分析法在攪拌器選型中的應(yīng)用[J].計算機與應(yīng)用化學(xué),1999,(3):207～210.

[5] Saaty T L.Axiomatic Foundation of the Analytic Hierarchy Progress[J].Management Science,1986,32(7):841～855.

[6] 鞏春領(lǐng).大跨度斜拉橋施工風(fēng)險分析與對策研究[D].上海：同濟(jì)大學(xué),2006.

[7] 王應(yīng)明.判斷矩陣排序方法綜述[J].管理科學(xué)學(xué)報,1995,(3):101～114.

[8] 侯岳衡,沈德家.指數(shù)標(biāo)度及其與幾種標(biāo)度的比較[J].系統(tǒng)工程理論與實踐,1995,15(10):43～46.

[9] 何堃.層次分析法的標(biāo)度研究[J].系統(tǒng)工程理論與實踐,1997,17(6):59～62,104.

[10] 魏翠萍.層次分析法中和積法的最優(yōu)化理論基礎(chǔ)及性質(zhì)[J].系統(tǒng)工程理論與實踐,1999,19(9):113～115,119.

[11] Scargiali F,Busciglio A,Grisafi F,et al.Mass Transfer and Hydrodynamic Characteristics of Unbaffled Stirred Bio-Reactors:Influence of Impeller Design[J].Biochemical Engineering Journal,2014,82(1):42～47.

[12] 吳祈宗,李有文.層次分析法中矩陣的判斷一致性研究[J].北京理工大學(xué)學(xué)報,1999,19(4):502～505.

[13] 薛勇勇,劉陽,劉琳琳，等.基于層次分析法優(yōu)化天然氣凈化工藝[J].化工進(jìn)展,2016,35(5):1298～1302.

[14] 王蓮芬.層次分析法中排序權(quán)數(shù)的計算方法[J].系統(tǒng)工程理論與實踐,1987,7(2):31～37.

[15] 駱正清,楊善林.層次分析法中幾種標(biāo)度的比較[J].系統(tǒng)工程理論與實踐,2004,24(9):51～60.

WeightCalculationMethodofAnalyticHierarchyProcessBasedonIndustrialAgitatorSelection

WU Chang-song,WU Meng-chu,LIANG Jia-yun,LEI Yu-tian,YANG Bin

(CollegeofChemicalEngineering,NorthwestUniversity)

In order to apply the analytic hierarchy process (AHP) to select industrial agitators,two kinds of commonly-used scales were adopted to establish the comparison matrix for the agitator selection and the weight vectors were solved by four matrix calculation methods to obtain the weight sequencing and the computational cost of each method.The sequencing result of comprehensive weights shows that,compared with the 1～9 scale methods,the 90/9～99/9scale method has less error and it’s more realistic and more suitable for AHP selection of the industrial agitators.Comparing the results of comprehensive weights and time consumption in the 90/9～99/9scales shows that,compared to the other three matrix calculation methods,the root mean square (RMS) method is more suitable for application of the agitator selection’s APH.On the basis of the same consistency and calculation accuracy,the calculation time of RMS is the least(i.e.the highest efficiency) and the time-consumption is 2.50% of the sum-product method (SPM),32.36% of the eigenvector method (EM) and 76.34% of the row and normalization (NRA) which shows that the RMS has the highest calculation efficiency.

agitator selection,AHP,scale,weight vector

國家自然科學(xué)基金項目(11102153)；陜西省自然科學(xué)基礎(chǔ)研究計劃青年人才資助項目(2013JQ2001)。

吳長松(1996-)，本科生，從事過程流體機械的研究。

聯(lián)系人楊斌(1978-)，副教授，從事過程流體機械、環(huán)境流體力學(xué)的研究，binyang@nwu.edu.cn。

TQ051.9

A

0254-6094(2017)04-0397-05

2017-03-18)