吳 悠，林圣強(qiáng)，周偉星，羅忠敬，楊 斌

(1.清華大學(xué)能源與動(dòng)力工程系燃燒能源中心，北京 100084；2.鹽城工學(xué)院汽車工程學(xué)院，鹽城 224051；3.哈爾濱工業(yè)大學(xué)能源科學(xué)與工程學(xué)院，哈爾濱150001)

燃燒反應(yīng)動(dòng)力學(xué)機(jī)理診斷技術(shù)在燃燒反應(yīng)動(dòng)力學(xué)機(jī)理模型構(gòu)建、機(jī)理模型可靠性分析及機(jī)理模型簡(jiǎn)化上都扮演著重要的角色，特別是敏感性分析方法.在機(jī)理模型構(gòu)建過(guò)程中，敏感性分析方法用于調(diào)整機(jī)理模型中的速率系數(shù)，使機(jī)理模型具有更高的預(yù)測(cè)能力；在特定工況下，基元反應(yīng)的敏感性系數(shù)也是判斷機(jī)理合理性的一個(gè)重要依據(jù)；同時(shí)，敏感性方法在機(jī)理簡(jiǎn)化中也扮演重要角色.因此，發(fā)展可判別機(jī)理中各組分的重要性程度的機(jī)理模型診斷技術(shù)對(duì)燃燒動(dòng)力系統(tǒng)的認(rèn)識(shí)及應(yīng)用都有著深刻影響.傳統(tǒng)的基元反應(yīng)敏感性方法只能判別基元反應(yīng)在燃燒動(dòng)力學(xué)系統(tǒng)中的重要性，但無(wú)法判別機(jī)理中各組分的重要性大小，因此其所呈現(xiàn)的機(jī)理診斷信息及其在機(jī)理簡(jiǎn)化上的應(yīng)用都受到一定程度的限制[1].

目前燃燒數(shù)值計(jì)算還無(wú)法耦合龐大的燃燒反應(yīng)動(dòng)力學(xué)機(jī)理[2-3]，因此在燃燒數(shù)值計(jì)算中耦合簡(jiǎn)化機(jī)理是研究發(fā)動(dòng)機(jī)流場(chǎng)結(jié)構(gòu)的重要手段，也是發(fā)展機(jī)理簡(jiǎn)化工具的主要目的之一.機(jī)理診斷方法的可靠性也可以通過(guò)簡(jiǎn)化燃燒反應(yīng)動(dòng)力學(xué)機(jī)理驗(yàn)證；可靠性高的診斷方法能夠精準(zhǔn)判別機(jī)理中組分的重要性，從而得到比較緊湊的簡(jiǎn)化機(jī)理.目前已發(fā)展出很多機(jī)理簡(jiǎn)化方法，如敏感性分析法(SA)[4-5]、主成分分析法(PCA)[6]、計(jì)算奇異攝動(dòng)法(CSP)[7]、有向關(guān)系圖解法(DRG)[8-9]、基于誤差傳播的直接關(guān)系圖解法(DRGEP)[10-12]和路徑通量分析法(PFA)等[13-14].這些簡(jiǎn)化方法各有優(yōu)勢(shì)，但其中的圖解法通常作為簡(jiǎn)化大分子燃料機(jī)理首選的簡(jiǎn)化手段.因此在比較簡(jiǎn)化能力時(shí)采用圖解法作為比較對(duì)象，可較為顯著地突出組分敏感性方法在診斷機(jī)理中組分重要性程度的能力.而敏感性分析法可分為局部敏感性分析法和全局敏感性分析法，燃燒過(guò)程中組分燃料演化過(guò)程快速且劇烈，因此基元反應(yīng)速率系數(shù)很難精確測(cè)量出，燃燒動(dòng)力學(xué)機(jī)理也帶有一定的速率系數(shù)不確定性.目前燃燒動(dòng)力學(xué)系統(tǒng)的分析手段主要是依靠局部敏感性方法，但眾多研究表明，全局敏感性方法可得到更加精確的參數(shù)敏感性數(shù)值；而且應(yīng)用全局簡(jiǎn)化方法可以得到更加緊湊的機(jī)理.這些都表明，在分析燃燒動(dòng)力學(xué)時(shí)其自身所帶的參數(shù)不確定性特性是不可以忽略的，但基于組分的全局敏感性方法(STGSA)[15]計(jì)算量較大，因此本文提出基于組分的半全局敏感性方法，以降低求解敏感性時(shí)所需的樣本量，并用于簡(jiǎn)化USCMech Ⅱ[16]、DME[17]機(jī)理和JetSurF[18]機(jī)理，通過(guò)機(jī)理簡(jiǎn)化過(guò)程驗(yàn)證其可靠性，突出其在機(jī)理診斷中的實(shí)用性.

1 半全局敏感性方法

敏感性系數(shù)的定義方式有很多種，這里介紹常用基元反應(yīng)敏感性方法的計(jì)算方法，以了解用基元敏感性診斷反應(yīng)動(dòng)力學(xué)系統(tǒng)的方式.隨之，在基元反應(yīng)敏感性方法的基礎(chǔ)上，提出基于組分半全局敏感性方法的思路及其計(jì)算過(guò)程，并通過(guò)簡(jiǎn)化動(dòng)力學(xué)系統(tǒng)來(lái)驗(yàn)證該方法診斷動(dòng)力學(xué)系統(tǒng)的可靠性.

1.1 基元反應(yīng)敏感性方法

敏感性分析方法被廣泛應(yīng)用在反應(yīng)動(dòng)力學(xué)軟件，像CHEMKIN、Cantera 等，其目的在于監(jiān)測(cè)反應(yīng)動(dòng)力學(xué)系統(tǒng)的信息.敏感性系數(shù)定義如下：

式中：S 為敏感性系數(shù)；Y 表示未知數(shù)向量，這里指燃燒特性參數(shù)(點(diǎn)火延遲時(shí)間、火焰?zhèn)鞑ニ俣鹊?；x表示動(dòng)力系統(tǒng)中的監(jiān)測(cè)參數(shù)，這里指反應(yīng)動(dòng)力學(xué)速率系數(shù)中的指前因子；下標(biāo)i 表示第i 個(gè)基元反應(yīng).為使不同單位的敏感性系數(shù)有可比性，通常使用無(wú)量綱的歸一化敏感性系數(shù)：

式(1)～(3)被廣泛應(yīng)用于反應(yīng)動(dòng)力學(xué)系統(tǒng)，以求解基元反應(yīng)敏感性系數(shù).此方法在詳細(xì)模型中可較為精確地評(píng)價(jià)基元反應(yīng)的重要性程度.

1.2 基于組分半全局敏感性方法

1.1節(jié)主要推導(dǎo)如何應(yīng)用敏感性方法評(píng)價(jià)基元反應(yīng)的重要性程度，并可知通過(guò)擾動(dòng)基元反應(yīng)的指前因子可輕松獲得基元反應(yīng)的重要性程度，但應(yīng)用敏感性方法評(píng)判組分重要性大小并不容易實(shí)現(xiàn).不同于基元反應(yīng)的速率系數(shù)有3 個(gè)固定參數(shù)，組分濃度以中間產(chǎn)物存在于求解動(dòng)力學(xué)系統(tǒng)過(guò)程中，因此監(jiān)測(cè)此動(dòng)態(tài)參數(shù)在動(dòng)力學(xué)機(jī)理中的重要性程度是比較困難的.另一個(gè)難點(diǎn)在于，局部敏感性方法是單一時(shí)刻、單一擾動(dòng)的診斷方法，其從系統(tǒng)中獲得的信息也相對(duì)有限，而“全視野”的全局敏感性可更加全面地分析多參數(shù)變化對(duì)系統(tǒng)輸出的影響，為此筆者提出基于組分的半全局敏感性方法.

基于組分半全局敏感性方法監(jiān)測(cè)動(dòng)力學(xué)系統(tǒng)中組分的重要性分為兩個(gè)部分，如圖1 所示.第一步是通過(guò)模擬點(diǎn)火過(guò)程得到每個(gè)基元反應(yīng)的反應(yīng)進(jìn)程量，在評(píng)價(jià)組分重要性時(shí)擾動(dòng)所有與該組分(組分A)相關(guān)的基元反應(yīng)的反應(yīng)速率系數(shù)，而組分A 不參與的基元反應(yīng)的反應(yīng)速率系數(shù)保持不變.第二步用Sobol方法[17]進(jìn)行采樣，生成N×d 維矩陣.矩陣中每一個(gè)數(shù)代表在不確定性范圍內(nèi)對(duì)組分A 的相關(guān)基元反應(yīng)的反應(yīng)速率系數(shù)的擾動(dòng)，由此可以得到N 個(gè)在某一特定擾動(dòng)時(shí)組分A 的敏感性系數(shù).計(jì)算N 個(gè)樣本中組分A 的敏感性系數(shù)的平均值和方差：

圖1 基于組分半全局敏感性方法的示意Fig.1 Schematic of species-targeted semi-global sensitiveity analysis

式中：N 為樣本數(shù)量；d 是輸入?yún)?shù)個(gè)數(shù)，即組分A的相關(guān)基元反應(yīng)數(shù)量；表示擾動(dòng)k 后的點(diǎn)火延遲時(shí)間；Y0表示原機(jī)理計(jì)算的點(diǎn)火延遲時(shí)間；ki,Ad表示第i 個(gè)樣本中與組分A 相關(guān)的第d 個(gè)基元反應(yīng)的反應(yīng)速率系數(shù).然后可以得到組分A 的半全局敏感性 系數(shù)：

式中：XspeA和σspeA分別為N 個(gè)樣本中組分A 的敏感性系數(shù)的平均值和方差.本工作中樣本數(shù)量N 為64，反應(yīng)速率系數(shù)的不確定度為2.

2 結(jié)果和討論

為了驗(yàn)證基于組分半全局敏感性機(jī)理簡(jiǎn)化方法的有效性，本文應(yīng)用此方法分別簡(jiǎn)化乙烯詳細(xì)機(jī)理USCMech Ⅱ、二甲醚詳細(xì)機(jī)理和正庚烷詳細(xì)機(jī)理JetSurF，并與通過(guò)DRG 和DRGEP 方法得到的簡(jiǎn)化機(jī)理進(jìn)行比較.最后計(jì)算多工況下簡(jiǎn)化機(jī)理的點(diǎn)火延遲時(shí)間并與詳細(xì)機(jī)理的計(jì)算結(jié)果進(jìn)行比較.

2.1 敏感性方法診斷動(dòng)力學(xué)系統(tǒng)

為比較不同敏感性方法的診斷結(jié)果，本文分別應(yīng)用基元反應(yīng)敏感性方法、基于組分半全局敏感性方法和基于組分全局敏感性方法診斷相同工況下乙烯的點(diǎn)火燃燒特性，分析的初始條件設(shè)置為：T0＝1 000 K、p0＝0.1 MPa、φ＝1.0.

如圖2 所示，3 種方法診斷動(dòng)力系統(tǒng)的結(jié)果是不一樣的，組分敏感性方法是監(jiān)測(cè)哪些組分對(duì)燃燒過(guò)程起促進(jìn)或抑制作用，而基元反應(yīng)敏感性方法是監(jiān)測(cè)哪些基元反應(yīng)對(duì)燃燒過(guò)程起促進(jìn)或抑制的作用.

從圖2(a)中可知，對(duì)點(diǎn)火過(guò)程起重要作用的基元反應(yīng)有：(R1)：C2H3＋O2＝CH2O＋HCO；(R2)：C2H3＋O2＝CH2CHO ＋O；(R3)：C2H4＋HO2＝CH2OCH2＋OH；而且從基元敏感性系數(shù)中可知，C2H3對(duì)點(diǎn)火過(guò)程的促進(jìn)和抑制均起重要作用.圖2(b)顯示了各初始工況下半全局組分敏感性系數(shù)最大的20 個(gè)組分，因組分敏感性系數(shù)是通過(guò)擾動(dòng)與某一組分相關(guān)的多個(gè)基元反應(yīng)的速率系數(shù)得到組分的敏感性系數(shù)，某一組分的敏感性系數(shù)越大，則該組分對(duì)點(diǎn)火延遲時(shí)間的影響越大.但從圖2(b)中可知，組分HO2對(duì)燃燒過(guò)程的增強(qiáng)效果并不是非常有效，若僅調(diào)整(R3)(C2H4＋HO2＝CH2OCH2＋OH)并不能有效地提高機(jī)理預(yù)測(cè)能力，需綜合考慮該組分參與的數(shù)個(gè)基元反應(yīng).在調(diào)節(jié)模型中不能僅僅依靠基元反應(yīng)敏感性方法的監(jiān)測(cè)結(jié)果，通過(guò)組分敏感性手段則可以提供另一個(gè)視角有效地提高動(dòng)力學(xué)機(jī)理的預(yù)測(cè)能力，因此組分敏感性方法將為動(dòng)力學(xué)機(jī)理構(gòu)建提供有力的診斷手段.圖2(c)顯示了初始工況下全局組分敏感性系數(shù)最大的20 個(gè)組分，對(duì)比圖2(b)和圖2(c)可知，全局敏感性系數(shù)均為正值，不能判斷組分對(duì)點(diǎn)火延遲時(shí)間的影響是促進(jìn)還是抑制；而半全局敏感性系數(shù)可通過(guò)數(shù)值的正負(fù)判斷組分對(duì)點(diǎn)火延遲時(shí)間的正負(fù)貢獻(xiàn).

圖2 應(yīng)用不同敏感性方法診斷動(dòng)力學(xué)系統(tǒng)得到的敏感性系數(shù)Fig.2 Sensitivity coefficients obtained by different sensitivity methods for diagnosing kinetic systems

從基于組分半全局敏感性方法的定義上看，基于組分半全局敏感性系數(shù)是通過(guò)擾動(dòng)多個(gè)基元反應(yīng)的速率系數(shù)得到組分的敏感性系數(shù)，是比基元反應(yīng)方法更加宏觀的反應(yīng)動(dòng)力學(xué)系統(tǒng)的診斷手段，與全局組分敏感性系數(shù)對(duì)比，增加了組分對(duì)反應(yīng)動(dòng)力學(xué)系統(tǒng)的正負(fù)貢獻(xiàn)的判斷.從上述分析過(guò)程中可知，基于組分半全局敏感性方法對(duì)機(jī)理構(gòu)建、機(jī)理的合理性驗(yàn)證都是強(qiáng)有力的工具，但其敏感性系數(shù)的結(jié)果是否可代表組 分在系統(tǒng)中的重要程度是需要有效地驗(yàn)證.因此，為了驗(yàn)證基于組分半全局敏感性方法診斷效果的可靠性，本文應(yīng)用該方法對(duì)3 個(gè)燃燒動(dòng)力學(xué)機(jī)理進(jìn)行簡(jiǎn)化，并且與DRG 及DRGEP 方法進(jìn)行比較，突出基于組分半全局敏感性方法可準(zhǔn)確量化組分的重要性 程度.

2.2 乙烯燃燒機(jī)理簡(jiǎn)化

Wang 等[16]發(fā)展的USCMech Ⅱ機(jī)理包含111 種組分和784 個(gè)基元反應(yīng)，本文用該機(jī)理來(lái)預(yù)測(cè)乙烯的點(diǎn)火特性，用基于組分半全局敏感性方法對(duì)其進(jìn)行機(jī)理簡(jiǎn)化，并與其他機(jī)理簡(jiǎn)化方法DRG、DRGEP 對(duì)比.簡(jiǎn)化時(shí)考慮寬廣范圍下的多個(gè)工況，初始條件設(shè)置為：壓力范圍為0.1～3 MPa，當(dāng)量比范圍為0.5～2，溫度范圍為1 000～1 600 K.

從圖3 中可知，與DRG 和DRGEP 相比，基于組分半全局敏感性方法可以獲得組分?jǐn)?shù)最少、最緊湊的簡(jiǎn)化機(jī)理，其在簡(jiǎn)化機(jī)理上的效率明顯高于其他兩種方法.當(dāng)最大相對(duì)誤差為10%時(shí)，基于組分半全局敏感性方法得到的簡(jiǎn)化機(jī)理僅含32 個(gè)組分，而DRGEP和DRG 得到的簡(jiǎn)化機(jī)理分別含77 個(gè)組分和93 個(gè)組分.基于組分全局敏感性方法(STGSA)在相似的條件下簡(jiǎn)化乙烯的詳細(xì)機(jī)理USCMech Ⅱ，最終可得到一個(gè)32 組分的簡(jiǎn)化機(jī)理.與STGSA 相比，半全局敏感性方法的簡(jiǎn)化效果略低于STGSA，但可節(jié)省其96%以上的計(jì)算成本.

圖3 C2H4 簡(jiǎn)化機(jī)理中組分?jǐn)?shù)目與最大預(yù)測(cè)誤差之間的關(guān)系Fig.3 Maximum relative error in ignition delay times as a function of the number of species in the reduced mechanisms for ethylene/air

由于USCMech II 是高溫機(jī)理，簡(jiǎn)化機(jī)理的驗(yàn)證范圍設(shè)置為：壓力范圍為0.1～3 MPa，當(dāng)量比范圍為0.5～2，溫度范圍為1 000～1 600 K.圖4 比較了32組分的簡(jiǎn)化機(jī)理與詳細(xì)機(jī)理預(yù)測(cè)點(diǎn)火延遲時(shí)間，從圖中可知32 組分的簡(jiǎn)化機(jī)理可以精準(zhǔn)預(yù)測(cè)點(diǎn)火延遲時(shí)間，多工況條件下的點(diǎn)火延遲時(shí)間最大預(yù)測(cè)誤差為9%.從該算例中可知，基于組分半全局敏感性方法可在高溫工況下有效簡(jiǎn)化詳細(xì)機(jī)理，與傳統(tǒng)的圖解法DRG、DRGEP 相比可獲得更加精簡(jiǎn)的簡(jiǎn)化機(jī)理.

圖4 不同工況下C2H4 簡(jiǎn)化機(jī)理與詳細(xì)機(jī)理的點(diǎn)火延遲時(shí)間Fig.4 Ignition delay times for ethylene/air mixture over a range of initial temperatures at various pressures and equivalence ratios，calculated with the detailedand 32-species reduced mechanism，respectively

2.3 二甲醚燃燒機(jī)理簡(jiǎn)化

Wang 等[17]發(fā)展的DME 機(jī)理包含56 種組分和301 個(gè)基元反應(yīng)，本文用該機(jī)理來(lái)預(yù)測(cè)二甲醚的點(diǎn)火特性，用基于組分半全局敏感性方法對(duì)其進(jìn)行機(jī)理簡(jiǎn)化，并與其他機(jī)理簡(jiǎn)化方法DRG、DRGEP 對(duì)比.簡(jiǎn)化時(shí)考慮寬廣范圍下的多個(gè)工況，初始條件設(shè)置為：壓力范圍為0.1～3 MPa，當(dāng)量比范圍為0.5～2，溫度范圍為600～1 200 K.

從圖5 中可知，與DRG 和DRGEP 相比，基于組分半全局敏感性方法可以獲得組分?jǐn)?shù)最少、最緊湊的簡(jiǎn)化機(jī)理，其在簡(jiǎn)化機(jī)理上的效率明顯高于其他兩種方法.當(dāng)最大相對(duì)誤差為10%時(shí)，基于組分半全局敏感性方法得到的簡(jiǎn)化機(jī)理僅含32 個(gè)組分，而DRGEP和DRG 在最大相對(duì)誤差為30%得到的簡(jiǎn)化機(jī)理分別含39 個(gè)組分和42 個(gè)組分.

圖5 DME 簡(jiǎn)化機(jī)理中組分?jǐn)?shù)目與最大預(yù)測(cè)誤差之間的關(guān)系Fig.5 Maximum relative error in ignition delay times as a function of the number of species in the reduced mechanisms for dimethyl ether/air

由于該機(jī)理是中低溫氧化機(jī)理，簡(jiǎn)化機(jī)理的驗(yàn)證范圍設(shè)置為：壓力范圍為0.1～3 MPa，當(dāng)量比范圍為0.5～2，溫度范圍為600～1 200 K.圖6 比較了32 組分的簡(jiǎn)化機(jī)理與詳細(xì)機(jī)理預(yù)測(cè)點(diǎn)火延遲時(shí)間，從圖中可知32 組分的簡(jiǎn)化機(jī)理可以精準(zhǔn)預(yù)測(cè)點(diǎn)火延遲時(shí)間，多工況條件下的點(diǎn)火延遲時(shí)間最大預(yù)測(cè)誤差為7%.從該算例中可知，基于組分半全局敏感性方法可在中低溫工況下有效簡(jiǎn)化詳細(xì)機(jī)理，與傳統(tǒng)的圖解法DRG、DRGEP 相比可獲得更加精簡(jiǎn)的簡(jiǎn)化機(jī)理.

圖6 不同工況DME簡(jiǎn)化機(jī)理與詳細(xì)機(jī)理的點(diǎn)火延遲時(shí)間Fig.6 Ignition delay times for dimethyl ether/air mixture over a range of initial temperatures at various pressures and equivalence ratios，calculated with the detailed and 32-species reduced mechanism，respectively

2.4 正庚烷燃燒機(jī)理簡(jiǎn)化

Wang 等[18]發(fā)展的JetSurf 1.0 機(jī)理包含194 種組分和1 459 個(gè)基元反應(yīng)，此機(jī)理是正構(gòu)烷烴中至正十二烷的高溫反應(yīng)動(dòng)力學(xué)模型.為驗(yàn)證基于組分半全局敏感性方法對(duì)大分子燃料的機(jī)理簡(jiǎn)化的有效性，采用該機(jī)理預(yù)測(cè)正庚烷的點(diǎn)火特性，并與 DRG、DRGEP 方法對(duì)比.簡(jiǎn)化時(shí)考慮寬廣范圍下的多個(gè)工況，初始條件設(shè)置為：壓力范圍為0.1～3 MPa，當(dāng)量比范圍為0.5～2，溫度范圍為1 000～1 600 K.

從圖7 中可知，與DRG 和DRGEP 相比，基于組分半全局敏感性方法可以獲得組分?jǐn)?shù)最少、最緊湊的簡(jiǎn)化機(jī)理，其在簡(jiǎn)化機(jī)理上的效率明顯高于其他兩種方法.當(dāng)最大相對(duì)誤差為10%時(shí)，基于組分半全局敏感性方法得到的簡(jiǎn)化機(jī)理僅含56 個(gè)組分，而DRGEP和DRG 得到的簡(jiǎn)化機(jī)理分別含79 個(gè)組分和92 個(gè)組分.STGSA 在相似的條件下簡(jiǎn)化正癸烷的詳細(xì)機(jī)理Jetsurf 可得到一個(gè)48 組分的簡(jiǎn)化機(jī)理.與STGSA相比，基于組分半全局敏感性方法的簡(jiǎn)化效果略低于STGSA，但可節(jié)省其96%以上的計(jì)算成本.

圖7 n-C7H16 簡(jiǎn)化機(jī)理中組分?jǐn)?shù)目與最大預(yù)測(cè)誤差之間的關(guān)系Fig.7 Maximum relative error in ignition delay times as a function of the number of species in the reduced mechanisms for n-heptane/air

由于此機(jī)理是高溫機(jī)理，所以簡(jiǎn)化機(jī)理的驗(yàn)證范圍設(shè)置為：壓力范圍為0.1～3 MPa，當(dāng)量比范圍為0.5～2，溫度范圍為1 000～1 600 K.圖8 比較了56組分的簡(jiǎn)化機(jī)理與詳細(xì)機(jī)理預(yù)測(cè)的點(diǎn)火延遲時(shí)間，從圖中可知56 組分的簡(jiǎn)化機(jī)理可以精準(zhǔn)預(yù)測(cè)點(diǎn)火延遲時(shí)間，多工況條件下的點(diǎn)火延遲時(shí)間最大預(yù)測(cè)誤差為7%.從該算例中可知，基于組分半全局敏感性方法可有效簡(jiǎn)化大分子燃料的詳細(xì)機(jī)理，與傳統(tǒng)的圖解法DRG、DRGEP 相比可獲得更加精簡(jiǎn)的簡(jiǎn)化機(jī)理.

圖8 不同工況n-C7 H16 簡(jiǎn)化機(jī)理與詳細(xì)機(jī)理的點(diǎn)火延遲時(shí)間Fig.8 Ignition delay times for n-heptane/air mixture over a range of initial temperatures at various pressures and equivalence ratios，calculated with the detailed and 56-species reduced mechanism，respectively

3 結(jié)語(yǔ)

本文提出可評(píng)價(jià)組分重要性的半全局敏感性方法，并應(yīng)用在機(jī)理簡(jiǎn)化上.從簡(jiǎn)化USCMechⅡ、DME和JetSurF 機(jī)理的算例中可知，基于組分半全局敏感性方法可得到比DRG、DRGEP 方法更加緊湊的簡(jiǎn)化動(dòng)力學(xué)機(jī)理，同時(shí)也說(shuō)明該方法不僅可以診斷組分在反應(yīng)動(dòng)力學(xué)系統(tǒng)中的作用程度及影響效果，還可以分辨組分的正負(fù)貢獻(xiàn).