本發(fā)明涉及過程控制領(lǐng)域，尤其涉及一種動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法及裝置。

背景技術(shù)：

反應(yīng)器模型的設(shè)計對反應(yīng)器的設(shè)計、生產(chǎn)操作條件的優(yōu)化及工藝的改進具有重要影響，因此，建立能準確計算反應(yīng)速率的動態(tài)連續(xù)攪拌釜反應(yīng)器的建模模型意義重大，現(xiàn)有動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法存在計算速度慢、計算精度不高的問題。

技術(shù)實現(xiàn)要素：

為了克服現(xiàn)有技術(shù)的不足，本發(fā)明的目的在于提供一種動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法及裝置，其能解決現(xiàn)有技術(shù)計算速度慢、計算精度不高的問題。

本發(fā)明的目的采用以下技術(shù)方案實現(xiàn)：

一種動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法，包括：

獲取反應(yīng)器內(nèi)各進料流股的摩爾流量值、各氣相出料的摩爾流量值、各液相出料的摩爾流量值、各進料流股的摩爾焓值、各液相出料的摩爾焓值、各氣相出料的摩爾焓值、反應(yīng)器熱量值、液相出料流股中各組分的摩爾分數(shù)值、氣體出料流股中各組分的摩爾分數(shù)值、進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值；

獲取初始溫度值、初始壓力值、各反應(yīng)物質(zhì)的初始組分及各反應(yīng)物質(zhì)的化學計量關(guān)系；

根據(jù)閃蒸算法計算出反應(yīng)器中的初始總摩爾數(shù)及初始摩爾焓值；

根據(jù)所述初始溫度值及所述各反應(yīng)物質(zhì)的化學計量關(guān)系，計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值；

根據(jù)所述初始總摩爾數(shù)、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述反應(yīng)器的總體積值、所述各反應(yīng)物質(zhì)的當前反應(yīng)速率值計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)；

根據(jù)所述初始總摩爾數(shù)、所述初始摩爾焓值、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述各進料流股的摩爾焓值、所述各液相出料的摩爾焓值、所述各氣相出料的摩爾焓值、所述反應(yīng)器熱量值、所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值；

根據(jù)所述初始總摩爾數(shù)、所述各反應(yīng)物質(zhì)的初始組分、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述液相出料流股中各組分的摩爾分數(shù)值、所述氣體出料流股中各組分的摩爾分數(shù)值、所述進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分；

根據(jù)所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)、所述經(jīng)過預設(shè)時間間隔后的摩爾焓值、所述經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分計算出經(jīng)過預設(shè)時間間隔后的溫度值；

根據(jù)所述經(jīng)過預設(shè)時間間隔后的溫度值計算出經(jīng)過預設(shè)時間間隔后各反應(yīng)物質(zhì)的反應(yīng)速率值。

優(yōu)選的，所述根據(jù)所述初始溫度值及所述各反應(yīng)物質(zhì)的化學計量關(guān)系，計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值包括：

根據(jù)公式k＝a*exp(-e/r*t)*t^b計算出反應(yīng)速率r；

根據(jù)反應(yīng)速率r及所述各反應(yīng)物質(zhì)的化學計量關(guān)系計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值，其中，k、a、e、b為常量，t為初始溫度值，ri為第i個反應(yīng)物質(zhì)的反應(yīng)速率，i為i＝1到nc，nc為反應(yīng)器中反應(yīng)物質(zhì)的個數(shù)。

優(yōu)選的，所述根據(jù)所述初始總摩爾數(shù)、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述反應(yīng)器的總體積值、所述各反應(yīng)物質(zhì)的當前反應(yīng)速率值計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)包括：

根據(jù)公式mt+δt＝mt+δt*(∑fin-fv-fl+v*σri)計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)，其中，δt為預設(shè)時間間隔，mt為初始總摩爾數(shù)、fin為各進料流股的摩爾流量值、fv為各氣相出料的摩爾流量值、fl為各液相出料的摩爾流量值、v為反應(yīng)器的總體積值、ri為各反應(yīng)物質(zhì)的當前反應(yīng)速率值，mt+δt為經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)。

優(yōu)選的，所述根據(jù)所述初始總摩爾數(shù)、所述初始摩爾焓值、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述各進料流股的摩爾焓值、所述各液相出料的摩爾焓值、所述各氣相出料的摩爾焓值、所述反應(yīng)器熱量值、所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值包括：

根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值，其中，ht為初始摩爾焓值、hin為各進料流股的摩爾焓值、hl為各液相出料的摩爾焓值、hv為各氣相出料的摩爾焓值、q為反應(yīng)器熱量值，ht+δt為經(jīng)過預設(shè)時間間隔后的摩爾焓值。

優(yōu)選的，所述根據(jù)所述初始總摩爾數(shù)、所述各反應(yīng)物質(zhì)的初始組分、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述液相出料流股中各組分的摩爾分數(shù)值、所述氣體出料流股中各組分的摩爾分數(shù)值、所述進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分包括：

根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分，其中，ci(t)為第i個反應(yīng)物質(zhì)的初始組分，ci(l)為液相出料流股中各組分的摩爾分數(shù)值，ci(v)為氣體出料流股中各組分的摩爾分數(shù)值，為ci(in)進料流股中各組分的摩爾分數(shù)值、ci(t+δt)為經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分。

本發(fā)明還涉及一種動態(tài)連續(xù)攪拌釜反應(yīng)器的建模裝置，其特征在于，包括：第一獲取模塊、第二獲取模塊、第一計算模塊、第二計算模塊、第三計算模塊、第四計算模塊、第五計算模塊、第六計算模塊及第七計算模塊；

所述第一獲取模塊，用于獲取反應(yīng)器內(nèi)各進料流股的摩爾流量值、各氣相出料的摩爾流量值、各液相出料的摩爾流量值、各進料流股的摩爾焓值、各液相出料的摩爾焓值、各氣相出料的摩爾焓值、反應(yīng)器熱量值、液相出料流股中各組分的摩爾分數(shù)值、氣體出料流股中各組分的摩爾分數(shù)值、進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值；

所述第二獲取模塊，用于獲取初始溫度值、初始壓力值、各反應(yīng)物質(zhì)的初始組分及各反應(yīng)物質(zhì)的化學計量關(guān)系；

所述第一計算模塊，用于根據(jù)閃蒸算法計算出反應(yīng)器中的初始總摩爾數(shù)及初始摩爾焓值；

所述第二計算模塊，用于根據(jù)所述初始溫度值及所述各反應(yīng)物質(zhì)的化學計量關(guān)系，計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值；

所述第三計算模塊，用于根據(jù)所述初始總摩爾數(shù)、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述反應(yīng)器的總體積值、所述各反應(yīng)物質(zhì)的當前反應(yīng)速率值計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)；

所述第四計算模塊，用于根據(jù)所述初始總摩爾數(shù)、所述初始摩爾焓值、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述各進料流股的摩爾焓值、所述各液相出料的摩爾焓值、所述各氣相出料的摩爾焓值、所述反應(yīng)器熱量值、所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值；

所述第五計算模塊，用于根據(jù)所述初始總摩爾數(shù)、所述各反應(yīng)物質(zhì)的初始組分、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述液相出料流股中各組分的摩爾分數(shù)值、所述氣體出料流股中各組分的摩爾分數(shù)值、所述進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分；

所述第六計算模塊，用于根據(jù)所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)、所述經(jīng)過預設(shè)時間間隔后的摩爾焓值、所述經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分計算出經(jīng)過預設(shè)時間間隔后的溫度值；

所述第七計算模塊，用于根據(jù)所述經(jīng)過預設(shè)時間間隔后的溫度值計算出經(jīng)過預設(shè)時間間隔后各反應(yīng)物質(zhì)的反應(yīng)速率值。

優(yōu)選的，所述第二計算模塊，還用于根據(jù)公式

k＝a*exp(-e/r*t)*t^b計算出反應(yīng)速率r；

根據(jù)反應(yīng)速率r及所述各反應(yīng)物質(zhì)的化學計量關(guān)系計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值，其中，k、a、e、b為常量，t為初始溫度值，ri為第i個反應(yīng)物質(zhì)的反應(yīng)速率，i為i＝1到nc，nc為反應(yīng)器中反應(yīng)物質(zhì)的個數(shù)。

優(yōu)選的，所述第三計算模塊，還用于根據(jù)公式

mt+δt＝mt+δt*(∑fin-fv-fl+v*σri)計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)，其中，δt為預設(shè)時間間隔，mt為初始總摩爾數(shù)、fin為各進料流股的摩爾流量值、fv為各氣相出料的摩爾流量值、fl為各液相出料的摩爾流量值、v為反應(yīng)器的總體積值、ri為各反應(yīng)物質(zhì)的當前反應(yīng)速率值，mt+δt為經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)。

優(yōu)選的，所述第四計算模塊，還用于根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值，其中，ht為初始摩爾焓值、hin為各進料流股的摩爾焓值、hl為各液相出料的摩爾焓值、hv為各氣相出料的摩爾焓值、q為反應(yīng)器熱量值，ht+δt為經(jīng)過預設(shè)時間間隔后的摩爾焓值。

優(yōu)選的，所述第五計算模塊，還用于根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分，其中，ci(t)為第i個反應(yīng)物質(zhì)的初始組分，ci(l)為液相出料流股中各組分的摩爾分數(shù)值，ci(v)為氣體出料流股中各組分的摩爾分數(shù)值，為ci(in)進料流股中各組分的摩爾分數(shù)值、ci(t+δt)為經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分。

相比現(xiàn)有技術(shù)，本發(fā)明的有益效果在于：通過設(shè)定初始值，計算出預設(shè)時間間隔后反應(yīng)器中的總摩爾數(shù)、摩爾焓值及各反應(yīng)物質(zhì)的組分，計算速度快且計算精度高。

附圖說明

圖1為本發(fā)明實施例提供的動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法流程圖；

圖2為本發(fā)明實施例提供的動態(tài)連續(xù)攪拌釜反應(yīng)器的建模裝置示意圖。

具體實施方式

為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點更加清楚明白，以下結(jié)合附圖及實施例，對本發(fā)明進行進一步詳細說明。應(yīng)當理解，此處所描述的具體實施例僅僅用以解釋本發(fā)明，并不用于限定本發(fā)明。

如圖1所示，本發(fā)明實施例提供的動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法，包括：

步驟s101：獲取反應(yīng)器內(nèi)各進料流股的摩爾流量值、各氣相出料的摩爾流量值、各液相出料的摩爾流量值、各進料流股的摩爾焓值、各液相出料的摩爾焓值、各氣相出料的摩爾焓值、反應(yīng)器熱量值、液相出料流股中各組分的摩爾分數(shù)值、氣體出料流股中各組分的摩爾分數(shù)值、進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值。

具體的，反應(yīng)器模型與其他模型配合使用，上述變量從其他模型計算獲取或根據(jù)反應(yīng)物質(zhì)本身的特性獲取。

步驟s102：獲取初始溫度值、初始壓力值、各反應(yīng)物質(zhì)的初始組分及各反應(yīng)物質(zhì)的化學計量關(guān)系。

步驟s103：根據(jù)閃蒸算法計算出反應(yīng)器中的初始總摩爾數(shù)及初始摩爾焓值。

具體的，可以設(shè)定初始壓力和總摩爾數(shù)初值，通過閃蒸計算，對反應(yīng)器進行初始化，得到初始總摩爾數(shù)及初始摩爾焓值。所述閃蒸計算為現(xiàn)有技術(shù)，在此不再贅述。

步驟s104：根據(jù)所述初始溫度值及所述各反應(yīng)物質(zhì)的化學計量關(guān)系，計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值。

具體的，根據(jù)公式k＝a*exp(-e/r*t)*t^b計算出反應(yīng)速率r；

根據(jù)反應(yīng)速率r及所述各反應(yīng)物質(zhì)的化學計量關(guān)系計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值，其中，k、a、e、b為常量，可以選擇標準值，也可以重新設(shè)定，t為初始溫度值，ri為第i個反應(yīng)物質(zhì)的反應(yīng)速率，i為i＝1到nc，nc為反應(yīng)器中反應(yīng)物質(zhì)的個數(shù)。以aa+bb＝cc+dd反應(yīng)為例，a、b、c、d各反應(yīng)物質(zhì)的化學計量系數(shù)，

步驟s105：根據(jù)所述初始總摩爾數(shù)、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述反應(yīng)器的總體積值、所述各反應(yīng)物質(zhì)的當前反應(yīng)速率值計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)。

具體的，根據(jù)公式mt+δt＝mt+δt*(∑fin-fv-fl+v*σri)計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)，其中，δt為預設(shè)時間間隔，mt為初始總摩爾數(shù)、fin為各進料流股的摩爾流量值、fv為各氣相出料的摩爾流量值、fl為各液相出料的摩爾流量值、v為反應(yīng)器的總體積值、ri為各反應(yīng)物質(zhì)的當前反應(yīng)速率值，mt+δt為經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)。

步驟s106：根據(jù)所述初始總摩爾數(shù)、所述初始摩爾焓值、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述各進料流股的摩爾焓值、所述各液相出料的摩爾焓值、所述各氣相出料的摩爾焓值、所述反應(yīng)器熱量值、所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值。

具體的，根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值，其中，ht為初始摩爾焓值、hin為各進料流股的摩爾焓值、hl為各液相出料的摩爾焓值、hv為各氣相出料的摩爾焓值、q為反應(yīng)器熱量值，ht+δt為經(jīng)過預設(shè)時間間隔后的摩爾焓值。

步驟s107：根據(jù)所述初始總摩爾數(shù)、所述各反應(yīng)物質(zhì)的初始組分、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述液相出料流股中各組分的摩爾分數(shù)值、所述氣體出料流股中各組分的摩爾分數(shù)值、所述進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分。

具體的，根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分，其中，ci(t)為第i個反應(yīng)物質(zhì)的初始組分，ci(l)為液相出料流股中各組分的摩爾分數(shù)值，ci(v)為氣體出料流股中各組分的摩爾分數(shù)值，為ci(in)進料流股中各組分的摩爾分數(shù)值、ci(t+δt)為經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分。

步驟s108：根據(jù)所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)、所述經(jīng)過預設(shè)時間間隔后的摩爾焓值、所述經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分計算出經(jīng)過預設(shè)時間間隔后的溫度值。

具體的，采用閃蒸算法計算出經(jīng)過預設(shè)時間間隔后的溫度值。

步驟s109：根據(jù)所述經(jīng)過預設(shè)時間間隔后的溫度值計算出經(jīng)過預設(shè)時間間隔后各反應(yīng)物質(zhì)的反應(yīng)速率值。

具體的，反應(yīng)速率的計算方法與步驟s104相同，根據(jù)預設(shè)時間間隔后各反應(yīng)物質(zhì)的反應(yīng)速率值，重新計算出經(jīng)過下一個預設(shè)時間間隔后的總摩爾數(shù)、摩爾焓值、及各反應(yīng)物質(zhì)的組分，不斷重復計算，對反應(yīng)器動態(tài)特性進行模擬和分析，從而優(yōu)化生產(chǎn)。

本發(fā)明實施例提供的動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法通過設(shè)定初始值，預設(shè)時間間隔內(nèi)通過反應(yīng)速率動態(tài)計算反應(yīng)器中的總摩爾數(shù)、摩爾焓值及各反應(yīng)物質(zhì)的組分，計算速度快且計算精度高。

本發(fā)明實施例還提供一種動態(tài)連續(xù)攪拌釜反應(yīng)器的建模裝置，如圖2所示，裝置包括：

第一獲取模塊11、第二獲取模塊12、第一計算模塊13、第二計算模塊14、第三計算模塊15、第四計算模塊16、第五計算模塊17、第六計算模塊18及第七計算模塊19；

第一獲取模塊11，用于獲取反應(yīng)器內(nèi)各進料流股的摩爾流量值、各氣相出料的摩爾流量值、各液相出料的摩爾流量值、各進料流股的摩爾焓值、各液相出料的摩爾焓值、各氣相出料的摩爾焓值、反應(yīng)器熱量值、液相出料流股中各組分的摩爾分數(shù)值、氣體出料流股中各組分的摩爾分數(shù)值、進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值；

第二獲取模塊12，用于獲取初始溫度值、初始壓力值、各反應(yīng)物質(zhì)的初始組分及各反應(yīng)物質(zhì)的化學計量關(guān)系；

第一計算模塊13，用于根據(jù)閃蒸算法計算出反應(yīng)器中的初始總摩爾數(shù)及初始摩爾焓值；

第二計算模塊14，用于根據(jù)所述初始溫度值及所述各反應(yīng)物質(zhì)的化學計量關(guān)系，計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值；

第三計算模塊15，用于根據(jù)所述初始總摩爾數(shù)、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述反應(yīng)器的總體積值、所述各反應(yīng)物質(zhì)的當前反應(yīng)速率值計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)；

第四計算模塊16，用于根據(jù)所述初始總摩爾數(shù)、所述初始摩爾焓值、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述各進料流股的摩爾焓值、所述各液相出料的摩爾焓值、所述各氣相出料的摩爾焓值、所述反應(yīng)器熱量值、所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值；

第五計算模塊17，用于根據(jù)所述初始總摩爾數(shù)、所述各反應(yīng)物質(zhì)的初始組分、所述各進料流股的摩爾流量值、所述各氣相出料的摩爾流量值、所述各液相出料的摩爾流量值、所述液相出料流股中各組分的摩爾分數(shù)值、所述氣體出料流股中各組分的摩爾分數(shù)值、所述進料流股中各組分的摩爾分數(shù)值、反應(yīng)器的總體積值計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分；

第六計算模塊18，用于根據(jù)所述經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)、所述經(jīng)過預設(shè)時間間隔后的摩爾焓值、所述經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分計算出經(jīng)過預設(shè)時間間隔后的溫度值；

第七計算模塊19，用于根據(jù)所述經(jīng)過預設(shè)時間間隔后的溫度值計算出經(jīng)過預設(shè)時間間隔后各反應(yīng)物質(zhì)的反應(yīng)速率值。

優(yōu)選的，第二計算模塊14，還用于根據(jù)公式

k＝a*exp(-e/r*t)*t^b計算出反應(yīng)速率r；

根據(jù)反應(yīng)速率r及所述各反應(yīng)物質(zhì)的化學計量關(guān)系計算出各反應(yīng)物質(zhì)的當前反應(yīng)速率值，其中，k、a、e、b為常量，t為初始溫度值，ri為第i個反應(yīng)物質(zhì)的反應(yīng)速率，i為i＝1到nc，nc為反應(yīng)器中反應(yīng)物質(zhì)的個數(shù)。

優(yōu)選的，第三計算模塊15，還用于根據(jù)公式

mt+δt＝mt+δt*(∑fin-fv-fl+v*σri)計算出經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)，其中，δt為預設(shè)時間間隔，mt為初始總摩爾數(shù)、fin為各進料流股的摩爾流量值、fv為各氣相出料的摩爾流量值、fl為各液相出料的摩爾流量值、v為反應(yīng)器的總體積值、ri為各反應(yīng)物質(zhì)的當前反應(yīng)速率值，mt+δt為經(jīng)過預設(shè)時間間隔后的總摩爾數(shù)。

優(yōu)選的，第四計算模塊16，還用于根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的摩爾焓值，其中，ht為初始摩爾焓值、hin為各進料流股的摩爾焓值、hl為各液相出料的摩爾焓值、hv為各氣相出料的摩爾焓值、q為反應(yīng)器熱量值，ht+δt為經(jīng)過預設(shè)時間間隔后的摩爾焓值。

優(yōu)選的，所述第五計算模塊17，還用于根據(jù)公式

計算出經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分，其中，ci(t)為第i個反應(yīng)物質(zhì)的初始組分，ci(l)為液相出料流股中各組分的摩爾分數(shù)值，ci(v)為氣體出料流股中各組分的摩爾分數(shù)值，為ci(in)進料流股中各組分的摩爾分數(shù)值、ci(t+δt)為經(jīng)過預設(shè)時間間隔后的各反應(yīng)物質(zhì)的組分。

本實施例中的裝置與前述實施例中的方法是基于同一發(fā)明構(gòu)思下的兩個方面，在前面已經(jīng)對方法實施過程作了詳細的描述，所以本領(lǐng)域技術(shù)人員可根據(jù)前述描述清楚地了解本實施例中的裝置的結(jié)構(gòu)及實施過程，為了說明書的簡潔，在此就不再贅述。

為了描述的方便，描述以上裝置時以功能分為各種模塊分別描述。當然，在實施本發(fā)明時可以把各模塊的功能在同一個或多個軟件和/或硬件中實現(xiàn)。

通過以上的實施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本發(fā)明可借助軟件加必需的通用硬件平臺的方式來實現(xiàn)?；谶@樣的理解，本發(fā)明的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，該計算機軟件產(chǎn)品可以存儲在存儲介質(zhì)中，如rom/ram、磁碟、光盤等，包括若干指令用以使得一臺計算機設(shè)備(可以是個人計算機，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本發(fā)明各個實施例或者實施例的某些部分所述的方法。

本發(fā)明提供的動態(tài)連續(xù)攪拌釜反應(yīng)器的建模方法及裝置通過設(shè)定初始值，計算出預設(shè)時間間隔后反應(yīng)器中的總摩爾數(shù)、摩爾焓值及各反應(yīng)物質(zhì)的組分，計算速度快且計算精度高。

對本領(lǐng)域的技術(shù)人員來說，可根據(jù)以上描述的技術(shù)方案以及構(gòu)思，做出其它各種相應(yīng)的改變以及形變，而所有的這些改變以及形變都應(yīng)該屬于本發(fā)明權(quán)利要求的保護范圍之內(nèi)。