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

吳長(zhǎng)松 吳夢(mèng)楚 梁佳赟 雷雨田 楊 斌

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

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

吳長(zhǎng)松 吳夢(mèng)楚 梁佳赟 雷雨田 楊 斌

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

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

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

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

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

1.1 1～9標(biāo)度

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

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

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

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

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

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

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

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

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

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

特征值法(EM)：

AW=λmaxW

(1)

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

方根法(RMS)：

(2)

和積法(SPM)[10]：

(3)

行和歸一化(NRA)：

(4)

最大特征值λmax的計(jì)算公式為：

(5)

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

3 實(shí)例計(jì)算與結(jié)果分析

3.1 建立比較矩陣

工業(yè)攪拌器實(shí)例來(lái)源于Scargiali F等的生物攪拌器研究[11]，表4為該文獻(xiàn)中列出的直葉圓盤(pán)渦輪攪拌器、開(kāi)啟渦輪攪拌器、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)則層對(duì)目標(biāo)層比較矩陣A，方案層對(duì)準(zhǔn)則層比較矩陣B1、B2、B3、B4[12]。

1～9標(biāo)度：

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

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

基于上述建立的比較矩陣，運(yùn)用EM法、RMS法、SPM法和NRA法解矩陣，得到總的權(quán)重系數(shù)[13，14]見(jiàn)表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)度計(jì)算得到的排序結(jié)果皆為A310型、斜葉開(kāi)啟渦輪、直葉圓盤(pán)渦輪和螺旋槳型攪拌器。結(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)度對(duì)應(yīng)效果較好，最終結(jié)果更符合實(shí)際。因此，應(yīng)用層次分析法進(jìn)行工業(yè)攪拌器選型可選擇90/9～99/9標(biāo)度建立比較矩陣。同樣由表5、6可知，4種計(jì)算方法求解權(quán)重向量的結(jié)果極為接近，精度基本一致。

3.3 分析計(jì)算耗時(shí)

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

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

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

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

4 結(jié)束語(yǔ)

針對(duì)實(shí)際工業(yè)攪拌器選型的層次分析法應(yīng)用，所選取的兩個(gè)標(biāo)度中，90/9～99/9標(biāo)度更符合實(shí)際，更適于層次分析法的攪拌器選型。在90/9～99/9標(biāo)度下比較綜合權(quán)重結(jié)果和耗時(shí)量可知，RMS法在同樣符合計(jì)算精度和排序一致性的前提下，相比其余3種矩陣計(jì)算方法耗時(shí)更少，更適于攪拌器選型的層次分析法應(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] 代麗,朱愛(ài)華,趙勻.應(yīng)用層次分析法計(jì)算分插機(jī)構(gòu)優(yōu)化目標(biāo)的權(quán)重[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(2):60～65.

[4] 王嘉駿,馮連芳,顧雪萍,等.層次分析法在攪拌器選型中的應(yīng)用[J].計(jì)算機(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)險(xiǎn)分析與對(duì)策研究[D].上海：同濟(jì)大學(xué),2006.

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

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

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

[10] 魏翠萍.層次分析法中和積法的最優(yōu)化理論基礎(chǔ)及性質(zhì)[J].系統(tǒng)工程理論與實(shí)踐,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é)報(bào),1999,19(4):502～505.

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

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

[15] 駱正清,楊善林.層次分析法中幾種標(biāo)度的比較[J].系統(tǒng)工程理論與實(shí)踐,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

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

吳長(zhǎng)松(1996-)，本科生，從事過(guò)程流體機(jī)械的研究。

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

TQ051.9

A

0254-6094(2017)04-0397-05

2017-03-18)