田 旺, 秦 康, 李明豐, 胡元沖, 梁家林, 褚小立
(中國(guó)石化 石油化工科學(xué)研究院,北京 100083)
隨著原油逐漸重質(zhì)化和環(huán)保要求日趨嚴(yán)格,重油加氫改質(zhì)變得越來(lái)越重要[1]。蠟油加氫作為重要的加氫改質(zhì)工藝,可脫除硫、氮等雜質(zhì),降低芳烴和烯烴含量,為催化裂化或加氫裂化提供適宜的原料[2]。蠟油加氫工藝的核心是脫硫、脫氮催化劑,所用催化劑不僅決定產(chǎn)品性質(zhì),而且影響裝置的運(yùn)行周期。
蠟油加氫催化劑的中試評(píng)價(jià)主要關(guān)注加氫后蠟油的硫、氮質(zhì)量分?jǐn)?shù),目前的中試評(píng)價(jià)模型多為機(jī)理模型。機(jī)理模型的解釋性強(qiáng),過程較為清晰,但在同時(shí)預(yù)測(cè)加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)時(shí),機(jī)理模型需要考慮兩者之間的相互作用,數(shù)學(xué)表達(dá)較為復(fù)雜,動(dòng)力學(xué)模型參數(shù)求解困難,實(shí)際使用并不方便。
神經(jīng)網(wǎng)絡(luò)模型不需要考慮復(fù)雜的機(jī)理就可以實(shí)現(xiàn)相應(yīng)任務(wù)的預(yù)測(cè),成為研究熱點(diǎn)。周軼峰等[3]利用3層神經(jīng)網(wǎng)絡(luò),根據(jù)柴油餾分油性質(zhì)、工藝條件對(duì)加氫柴油中硫質(zhì)量分?jǐn)?shù)進(jìn)行預(yù)測(cè),得出不同工藝條件對(duì)加氫脫硫(HDS)反應(yīng)深度的影響大小。郭慶洲等[4]根據(jù)加氫處理潤(rùn)滑油基礎(chǔ)油的幾種性質(zhì),建立了預(yù)測(cè)其芳烴質(zhì)量分?jǐn)?shù)的神經(jīng)網(wǎng)絡(luò)模型。相比主成分回歸模型,神經(jīng)網(wǎng)絡(luò)模型的預(yù)測(cè)精度更高。胡元沖等[5]基于神經(jīng)網(wǎng)絡(luò)框架Keras,針對(duì)柴油產(chǎn)物中硫、氮、單環(huán)芳烴、多環(huán)芳烴質(zhì)量分?jǐn)?shù)分別建立了預(yù)測(cè)模型。李詔陽(yáng)等[6]采用長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(Long short-term memory,LSTM),以加氫裂化的操作參數(shù)為輸入變量,輕質(zhì)潤(rùn)滑油的4個(gè)性質(zhì)為輸出變量,建立了預(yù)測(cè)模型;通過隨機(jī)森林(Random forest,RF)算法對(duì)輸入特征進(jìn)行重要度排序,得到不同生產(chǎn)目標(biāo)對(duì)應(yīng)的優(yōu)化操作方案。田水苗等[7]利用BP神經(jīng)網(wǎng)絡(luò)建立數(shù)據(jù)驅(qū)動(dòng)模型,來(lái)預(yù)測(cè)石腦油、液化氣、燃料氣、精制蠟油流量以及精制蠟油中硫、氮質(zhì)量分?jǐn)?shù);仿真結(jié)果表明BP神經(jīng)網(wǎng)絡(luò)模型具有較高的預(yù)測(cè)精度。
盡管神經(jīng)網(wǎng)絡(luò)可同時(shí)預(yù)測(cè)多任務(wù),但是,多任務(wù)同時(shí)預(yù)測(cè)易出現(xiàn)各任務(wù)之間梯度不平衡、預(yù)測(cè)結(jié)果偏向某個(gè)任務(wù)的問題。為解決多任務(wù)同時(shí)優(yōu)化的平衡問題,Chen等[8]曾提出梯度歸一化算法,通過動(dòng)態(tài)調(diào)整梯度大小,自動(dòng)平衡深度學(xué)習(xí)模型的多任務(wù)優(yōu)化;但該方法只在特定的優(yōu)化器下才能表現(xiàn)出較好結(jié)果。Zhou等[9]針對(duì)多視圖分辨問題,考慮不同任務(wù)之間的方差不確定性,提出協(xié)同多目標(biāo)量子粒子群優(yōu)化算法(Cooperative multi-objective quantum particle swarm optimization, CMOQPSO),該算法即使在小樣本場(chǎng)景下也能得到較高的分類準(zhǔn)確率。Guo等[10]引入動(dòng)態(tài)任務(wù)優(yōu)先級(jí)概念,通過自適應(yīng)調(diào)整各任務(wù)損失函數(shù)的混合權(quán)重,自動(dòng)對(duì)難度較高的任務(wù)進(jìn)行優(yōu)先級(jí)排序;分類模型利用該算法辨別人類游泳姿態(tài)的準(zhǔn)確率甚至可以媲美單任務(wù)學(xué)習(xí)。Kendall等[11]考慮到視覺場(chǎng)景理解背景下的多任務(wù)學(xué)習(xí),涉及到不同單位和尺度的各種回歸和分類任務(wù),用各任務(wù)的同方差不確定性權(quán)衡多個(gè)損失函數(shù);該方法能夠自動(dòng)從數(shù)據(jù)中學(xué)習(xí)相對(duì)權(quán)重,且不受初始化權(quán)重的影響。Yu等[12]確定了3種任務(wù)同時(shí)優(yōu)化時(shí)各任務(wù)梯度相互干擾的客觀因素,開發(fā)了減輕任務(wù)梯度之間相互干擾的簡(jiǎn)單通用方法。該方法將某個(gè)任務(wù)的梯度投影到所有與其具有干擾的任務(wù)梯度平面,對(duì)一系列具有挑戰(zhàn)性的多任務(wù)監(jiān)督學(xué)習(xí)和多任務(wù)強(qiáng)化學(xué)習(xí)問題,顯著提高了模型的性能和效率,也明顯提高了模型的訓(xùn)練速度和預(yù)測(cè)效果。Ozan等[13]將多任務(wù)學(xué)習(xí)視為尋找Pareto最優(yōu)解的多目標(biāo)優(yōu)化,在界定出一個(gè)多目標(biāo)損失函數(shù)的上界并證明它可以有效地優(yōu)化之后,進(jìn)一步證明該上界優(yōu)化后能得到Pareto最優(yōu)解;該方法在數(shù)字分類、聯(lián)合語(yǔ)義分割、實(shí)例分割、深度估計(jì)、多標(biāo)簽分類等場(chǎng)景下能夠訓(xùn)練出高效的分類模型。
多任務(wù)學(xué)習(xí)中,對(duì)各任務(wù)預(yù)測(cè)的平衡問題研究,目前主要集中在計(jì)算機(jī)視覺、圖像分類和識(shí)別、自然語(yǔ)言處理等領(lǐng)域;對(duì)算法研究,主要關(guān)注算法的改進(jìn)和統(tǒng)計(jì)學(xué)的優(yōu)化;應(yīng)用對(duì)象主要針對(duì)分類問題。但是,從影響梯度平衡的本質(zhì)——不同任務(wù)數(shù)據(jù)分布的固有差異出發(fā),解決不同任務(wù)在優(yōu)化中出現(xiàn)梯度差異的回歸問題仍未得到充分的關(guān)注。
筆者以非線性多任務(wù)回歸問題為研究對(duì)象,根據(jù)加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)數(shù)據(jù)的分布特點(diǎn),提出了一種構(gòu)建損失函數(shù)的新方法,將其應(yīng)用于蠟油加氫催化劑中試評(píng)價(jià)數(shù)據(jù)模型。針對(duì)4種不同的蠟油加氫催化劑,基于深度學(xué)習(xí)框架,通過硬參數(shù)共享,建立了適用于4種催化劑的中試評(píng)價(jià)數(shù)據(jù)驅(qū)動(dòng)模型;同時(shí)通過全新的損失函數(shù),解決了模型訓(xùn)練過程中各任務(wù)梯度的平衡問題。
選用減壓蠟油和焦化蠟油的混合原料為研究對(duì)象,用以考察不同催化劑對(duì)原料的適用性。原料蠟油(記為樣品A、B、C、D)的主要性質(zhì)如表1所示。實(shí)驗(yàn)所用催化劑為中國(guó)石化石油化工科學(xué)研究院研制的4種催化劑,即2種NiMo劑(RS-2100和RN-410)、CoMo劑(RVS-420)和NiMoW劑(RN-32V)。

表1 蠟油加氫實(shí)驗(yàn)所用原料油主要性質(zhì)Table 1 Properties of raw oil for wax oil hydrogenation
蠟油加氫催化劑中試評(píng)價(jià)的工藝條件為:反應(yīng)溫度355~390 ℃,反應(yīng)壓力6.4~15 MPa,反應(yīng)空速0.6~1.6 h-1,氫/油體積比600~1200。其工藝流程如圖1所示,其中,加氫反應(yīng)器中裝填直徑2~3 mm 的待評(píng)價(jià)催化劑100~150 mL。原料蠟油經(jīng)加熱爐加熱并與H2混合后,進(jìn)入加氫反應(yīng)器,加氫反應(yīng)產(chǎn)物通過熱高分罐進(jìn)行分離,氣體產(chǎn)物經(jīng)洗滌后釋放,加氫蠟油進(jìn)入產(chǎn)品罐,離線分析加氫蠟油的硫、氮含量。

圖1 蠟油加氫催化劑中試評(píng)價(jià)工藝流程圖Fig.1 Pilot-scale evaluation flow chart of wax oil hydrogenation catalyst
蠟油加氫催化劑由Co、Mo、Ni、W 4種元素中的2種或3種混合制備而成,不同催化劑的脫硫、脫氮能力有所差異。為此,模型基于深度學(xué)習(xí)框架Keras,對(duì)不同的催化劑使用獨(dú)熱編碼,完成不同催化劑中試評(píng)價(jià)數(shù)據(jù)的并行訓(xùn)練,構(gòu)建適用于4種催化劑的數(shù)據(jù)模型。模型的輸入?yún)?shù)包括原料性質(zhì)、操作參數(shù)、催化劑牌號(hào);輸出參數(shù)為加氫蠟油硫、氮質(zhì)量分?jǐn)?shù)。建立的模型可以實(shí)現(xiàn)加氫蠟油硫、氮質(zhì)量分?jǐn)?shù)的同時(shí)預(yù)測(cè)。
利用Python編寫代碼,調(diào)用開源的深度學(xué)習(xí)框架Keras,搭建如圖2所示的模型架構(gòu)。模型設(shè)置4層網(wǎng)絡(luò)結(jié)構(gòu),分別為輸入層、第一隱含層、第二隱含層、輸出層。為避免出現(xiàn)過擬合現(xiàn)象,第一和第二隱含層的Dropout均設(shè)為0.1。隱含層傳遞函數(shù)分別為線性整流函數(shù)(Rectified linear unit,ReLU)和雙曲正切函數(shù)(Hyperbolic tangent,tanh),輸出層傳遞函數(shù)為線性函數(shù)linear。損失函數(shù)為均方誤差(Mean square error,MSE),采用Adam優(yōu)化器;模型評(píng)估指標(biāo)為平均絕對(duì)誤差(Mean absolute error,MAE)。為減少模型的隨機(jī)性,將批量處理規(guī)模batch_size設(shè)置為5,為保證模型可以復(fù)現(xiàn),固定初始化隨機(jī)種子。
調(diào)用機(jī)器學(xué)習(xí)庫(kù)Sklearn的數(shù)據(jù)集劃分包,按照70∶15∶15的比例將所有數(shù)據(jù)隨機(jī)劃分為訓(xùn)練集(Train set)、驗(yàn)證集(Validation set)和測(cè)試集(Test set)。為提高模型訓(xùn)練速度,使用Standard Scaler函數(shù)對(duì)輸入數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理,對(duì)輸出數(shù)據(jù)取對(duì)數(shù)。基于Keras進(jìn)行模型訓(xùn)練的迭代過程如圖2所示。

圖2 基于Keras的深度學(xué)習(xí)框架模型Fig.2 Deep learning framework model based on Keras
圖3為蠟油加氫中試數(shù)據(jù)模型訓(xùn)練集和驗(yàn)證集的MSE和MAE變化。由圖3(a)可見,隨著迭代次數(shù)(Epoch)的增加,訓(xùn)練集和驗(yàn)證集的均方誤差MSE迅速減小,并逐步逼近穩(wěn)定值,減小的幅度呈良好一致性,說(shuō)明模型訓(xùn)練過程合理。將平均絕對(duì)誤差MAE作為模型訓(xùn)練效果的評(píng)價(jià)指標(biāo),由圖3(b)可見,隨著Epoch的增加,訓(xùn)練集和驗(yàn)證集的MAE迅速減小,兩者之間的差值逐漸縮小并趨于相同。

MSE—Mean square error; MAE—Mean absolute error; Epoch—Iterative times圖3 蠟油加氫中試數(shù)據(jù)模型訓(xùn)練集和驗(yàn)證集的MSE和MAE變化Fig.3 MSE and MAE of training set and validation set in training process of pilot test data evaluation model for wax oil hydrogenation(a) MSE vs. Epoch; (b) MAE vs. Epoch
圖4為測(cè)試集樣本中加氫蠟油硫、氮質(zhì)量分?jǐn)?shù)的實(shí)測(cè)值和預(yù)測(cè)值對(duì)比。由圖4可見:在0~1000 μg/g區(qū)間內(nèi),模型對(duì)硫、氮質(zhì)量分?jǐn)?shù)的預(yù)測(cè)值與實(shí)測(cè)值吻合良好;在高于1000 μg/g區(qū)間,硫質(zhì)量分?jǐn)?shù)預(yù)測(cè)值與實(shí)測(cè)值的偏離明顯,而氮質(zhì)量分?jǐn)?shù)無(wú)實(shí)測(cè)數(shù)據(jù)。

圖4 測(cè)試集樣本中加氫蠟油硫、氮質(zhì)量分?jǐn)?shù)的實(shí)測(cè)值(wmea)和預(yù)測(cè)值(wpre)對(duì)比Fig.4 Comparison of the measured (wmea) andpredicted values (wpre) of sulfur and nitrogen massfractions of hydrogenated wax oil in the test set
為了進(jìn)一步定量測(cè)試集加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)預(yù)測(cè)效果的差異,筆者對(duì)模型的預(yù)測(cè)效果進(jìn)行了統(tǒng)計(jì)。表2為舊的損失函數(shù)模型對(duì)測(cè)試集中硫、氮質(zhì)量分?jǐn)?shù)預(yù)測(cè)效果的評(píng)價(jià)指標(biāo)。從表2可以看出,測(cè)試集加氫蠟油硫、氮質(zhì)量分?jǐn)?shù)預(yù)測(cè)效果的評(píng)價(jià)指標(biāo)差異較大。分析原因,主要是模型的損失函數(shù)對(duì)硫、氮質(zhì)量分?jǐn)?shù)進(jìn)行同時(shí)優(yōu)化時(shí),硫、氮質(zhì)量分?jǐn)?shù)差異大,導(dǎo)致兩者梯度更新時(shí)產(chǎn)生沖突。

表2 舊的損失函數(shù)模型對(duì)測(cè)試集輸出參數(shù)的預(yù)測(cè)效果Table 2 The prediction effect of the old loss function modelon the output parameters of the test set
模型訓(xùn)練過程中,對(duì)2個(gè)任務(wù)進(jìn)行優(yōu)化時(shí),如果一個(gè)任務(wù)的梯度比另一個(gè)任務(wù)的梯度大得多,那么它將主導(dǎo)平均梯度。此研究中,模型的損失函數(shù)表達(dá)式為均方誤差(MSE),即各任務(wù)預(yù)測(cè)值與實(shí)測(cè)值誤差的平方和再加和求平均,如式(1)所示。

(1)

上述模型數(shù)據(jù)集中,加氫蠟油的硫質(zhì)量分?jǐn)?shù)分布在0~6000 μg/g之間,而加氫蠟油的氮質(zhì)量分?jǐn)?shù)分布在0~1000 μg/g之間。模型在訓(xùn)練時(shí),會(huì)優(yōu)先滿足氮質(zhì)量分?jǐn)?shù)的預(yù)測(cè)誤差,忽略硫質(zhì)量分?jǐn)?shù)的預(yù)測(cè)誤差,從而出現(xiàn)硫、氮質(zhì)量分?jǐn)?shù)預(yù)測(cè)效果差異大的情況。
為解決非線性多任務(wù)學(xué)習(xí)模型預(yù)測(cè)過程中各任務(wù)預(yù)測(cè)誤差不平衡的問題,筆者對(duì)所建模型的損失函數(shù)進(jìn)行了重新設(shè)計(jì),分別給硫、氮質(zhì)量分?jǐn)?shù)的誤差平方和乘以不同的權(quán)重,該權(quán)重的作用是平衡模型訓(xùn)練過程中各任務(wù)的梯度,使各任務(wù)的梯度保持一致,同時(shí)達(dá)到最優(yōu)訓(xùn)練效果。假設(shè)第i個(gè)任務(wù)的損失函數(shù)(Li)由式(2)計(jì)算得到;多任務(wù)的損失函數(shù)(LMTL)由式(3)計(jì)算得到;式(4)中共享參數(shù)(Wsh,下標(biāo)sh是share的縮寫)的優(yōu)化受到所有任務(wù)損失函數(shù)的影響,并且不同任務(wù)損失函數(shù)對(duì)于共享參數(shù)的影響可以使用權(quán)重WFi進(jìn)行調(diào)節(jié),其中,Wsh+1為共享參數(shù)Wsh的下一步迭代,γ為參數(shù)更新過程中的步長(zhǎng)。
(2)
(3)
(4)
假設(shè)賦予加氫蠟油中氮質(zhì)量分?jǐn)?shù)的權(quán)重為WF1,賦予加氫蠟油中硫質(zhì)量分?jǐn)?shù)的權(quán)重為WF2, 訓(xùn)練模型數(shù)據(jù)集加氫蠟油的硫質(zhì)量分?jǐn)?shù)平均值為X,訓(xùn)模模型數(shù)據(jù)集加氫蠟油氮質(zhì)量分?jǐn)?shù)的平均值為Y,WF1和WF2與X和Y的關(guān)系如式(5)所示。
(5)
將式(5)代入式(3),得到模型的損失函數(shù)如式(6)所示。

(6)
對(duì)模型損失函數(shù)的代碼進(jìn)行重新編寫,模型訓(xùn)練時(shí),調(diào)用新的損失函數(shù),即式(6),模型訓(xùn)練完成后,對(duì)測(cè)試集加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)進(jìn)行同時(shí)預(yù)測(cè)。
舊、新?lián)p失函數(shù)模型對(duì)測(cè)試集加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)的預(yù)測(cè)效果對(duì)比,如圖5所示。修改損失函數(shù)前,舊模型的預(yù)測(cè)值為Predicted1;修改損失函數(shù)后,新模型的預(yù)測(cè)值為Predicted2。從圖5可以看出:相比Predicted1,在硫質(zhì)量分?jǐn)?shù)高于1000 μg/g區(qū)間,Predicted2與實(shí)測(cè)值的誤差明顯減小;在硫質(zhì)量分?jǐn)?shù)低于1000 μg/g區(qū)間,Predicted2與實(shí)測(cè)值的誤差略微減小。對(duì)于加氫蠟油氮質(zhì)量分?jǐn)?shù)的預(yù)測(cè)而言,Predicted2與實(shí)測(cè)值的誤差比Predicted1與實(shí)測(cè)值的誤差略微降低。這說(shuō)明新模型的損失函數(shù)可同時(shí)縮小加氫蠟油硫、氮質(zhì)量分?jǐn)?shù)的預(yù)測(cè)誤差,其中硫質(zhì)量分?jǐn)?shù)的預(yù)測(cè)誤差縮小尤其顯著。
修改損失函數(shù)前、后,模型對(duì)測(cè)試集的預(yù)測(cè)評(píng)估指標(biāo),如表3所示。從表3可以看出,構(gòu)建新的損失函數(shù)用于模型訓(xùn)練,加氫蠟油中硫質(zhì)量分?jǐn)?shù)的平均絕對(duì)誤差(MAE)、平均相對(duì)誤差(MRE)明顯降低;而加氫蠟油中氮質(zhì)量分?jǐn)?shù)的MAE、MRE略微降低。這說(shuō)明新的損失函數(shù)可改善模型對(duì)不同任務(wù)的預(yù)測(cè)效果。其中,對(duì)硫任務(wù)的調(diào)節(jié)作用明顯,對(duì)比舊、新?lián)p失函數(shù)模型,加氫蠟油硫質(zhì)量分?jǐn)?shù)的MAE從125.36 μg/g下降至49.89 μg/g,MRE從32.60%下降至9.98%,R2從0.98提升至0.99,說(shuō)明為硫任務(wù)損失函數(shù)設(shè)置的權(quán)重,對(duì)其梯度更新,起到了較好的修飾作用,使得模型在尋優(yōu)的過程中,逐漸向硫任務(wù)的全局最優(yōu)解靠攏。而模型對(duì)氮任務(wù)的調(diào)節(jié)作用較弱,加氫蠟油氮質(zhì)量分?jǐn)?shù)的MAE、MRE降低幅度較小,R2沒有變化,說(shuō)明為氮任務(wù)損失函數(shù)設(shè)置的權(quán)重,對(duì)其梯度更新也能起到修飾作用,但這種作用較弱,更多的表現(xiàn)為削弱硫任務(wù)梯度更新對(duì)氮任務(wù)梯度更新的影響。

圖5 新、舊損失函數(shù)模型對(duì)測(cè)試集輸出參數(shù)預(yù)測(cè)效果比較Fig.5 Comparison between old and new loss function models in predicting output parameters of the test set(a) w(Sulfur); (b) w(Nitrogen)

表3 新、舊損失函數(shù)模型對(duì)測(cè)試集預(yù)測(cè)評(píng)估指標(biāo)統(tǒng)計(jì)Table 3 New and old loss function models predict the evaluation index statistics of the test set
綜合來(lái)看,造成模型同時(shí)預(yù)測(cè)加氫蠟油硫、氮質(zhì)量分?jǐn)?shù)預(yù)測(cè)誤差差異大的關(guān)鍵原因是硫、氮質(zhì)量分?jǐn)?shù)分布的固有差異,這種差異導(dǎo)致模型在訓(xùn)練過程中,硫、氮任務(wù)梯度更新不一致,進(jìn)而出現(xiàn)氮任務(wù)逼近最優(yōu)解而硫任務(wù)未逼近最優(yōu)解的情況。因此,筆者基于加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)的分布規(guī)律,分別為硫、氮任務(wù)損失函數(shù)構(gòu)造不同的權(quán)重,來(lái)平衡兩者之間的差異。使用全新的數(shù)據(jù)進(jìn)行測(cè)試,結(jié)果表明新的損失函數(shù)基本解決了模型對(duì)2個(gè)任務(wù)同時(shí)預(yù)測(cè)時(shí)出現(xiàn)任務(wù)之間預(yù)測(cè)誤差不平衡的問題。
(1)基于深度學(xué)習(xí)框架Keras,利用硬參數(shù)共享,建立了適用于4種不同蠟油加氫催化劑中試評(píng)價(jià)的數(shù)據(jù)驅(qū)動(dòng)模型,實(shí)現(xiàn)了加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)的同時(shí)預(yù)測(cè)。
(2)根據(jù)加氫蠟油中硫、氮質(zhì)量分?jǐn)?shù)分布的差異,構(gòu)建新的損失函數(shù),解決了模型同時(shí)預(yù)測(cè)2個(gè)任務(wù)時(shí),出現(xiàn)的預(yù)測(cè)誤差不平衡的問題。相比舊的損失函數(shù),新的損失函數(shù)模型模擬結(jié)果表明,同時(shí)預(yù)測(cè)精制蠟油硫、氮質(zhì)量分?jǐn)?shù)的平均絕對(duì)誤差分別從125.36和45.95 μg/g降低至49.89和38.62 μg/g,平均相對(duì)誤差分別從32.6%和13.4%降低至9.98%和9.59%,基本滿足實(shí)際應(yīng)用的要求。
(3)蠟油加氫工藝中的脫硫和脫氮反應(yīng)是一種相互抑制的反應(yīng),傳統(tǒng)的機(jī)理模型構(gòu)建時(shí)需要考慮兩者之間的相互作用,加入一定的假設(shè),參數(shù)的求解過程較為復(fù)雜。筆者提出的非線性多任務(wù)同時(shí)預(yù)測(cè)的方法,為該問題的解決提供了另一種思路。