本發(fā)明實施例涉及石油化工
技術(shù)領(lǐng)域：
，具體涉及一種催化裂化反應(yīng)產(chǎn)物軟測量的方法及系統(tǒng)。
背景技術(shù)：
：催化裂化裝置(fccu,fluidcatalyticcrackingunit)是煉化企業(yè)的關(guān)鍵裝置，是重質(zhì)油深度加工、輕質(zhì)化的重要手段，隨著原油的重質(zhì)化、劣質(zhì)化加劇，催化裂化裝置的作用更加突出，在煉油工業(yè)中占有舉足輕重的地位。重質(zhì)原油通過催化裂化裝置被輕質(zhì)化，生成干氣(c1、c2、h2、h2s)、液化氣(c3、c4)、汽油(c5～c11)、柴油(c10～c20)、油漿(以稠環(huán)芳烴為主)、焦炭(縮合產(chǎn)物)等。受限于裝置尺寸、測量技術(shù)等諸多因素以及催化裂化反應(yīng)自身的復(fù)雜流動反應(yīng)體系高度耦合的特征，催化裂化反應(yīng)產(chǎn)物中的產(chǎn)品產(chǎn)率是無法使用儀表直接測量得到的，因此利用在線可測的輔助變量，如溫度、壓力、流量來實時估計難以檢測的產(chǎn)品產(chǎn)率的軟測量技術(shù)在催化裂化裝置中得到了廣泛應(yīng)用?，F(xiàn)有技術(shù)中，催化裂化軟測量技術(shù)大多使用集總動力學(xué)模型，通過集總動力學(xué)模型，獲得催化裂化反應(yīng)的產(chǎn)品產(chǎn)率。但是這種傳統(tǒng)的集總反應(yīng)動力學(xué)模型在建模過程中，只注重裂化反應(yīng)動力學(xué)過程，假定反應(yīng)氣體的狀態(tài)為等溫活塞流，而忽略了催化裂化系統(tǒng)中流動、傳熱、傳質(zhì)以及湍流脈動對反應(yīng)的影響，使得數(shù)值模擬存在一定的偏差，對催化裂化反應(yīng)軟測量結(jié)果的精度造成一定的影響。而且，傳統(tǒng)方法下得到的模型參數(shù)適應(yīng)性較差，往往是針對特定反應(yīng)條件和工況下得到的結(jié)果，因此當(dāng)反應(yīng)條件、催化劑性質(zhì)等因素變化時，需要通過現(xiàn)場采集數(shù)據(jù)來進(jìn)行模型校正與補(bǔ)償。因此，如何提出一種方案，能夠提高催化裂化反應(yīng)軟測量方法的適應(yīng)性，進(jìn)一步提高測量結(jié)果的精度，成為亟待解決的問題。技術(shù)實現(xiàn)要素：針對現(xiàn)有技術(shù)中的缺陷，本發(fā)明實施例提供了及一種催化裂化反應(yīng)產(chǎn)物軟測量的方法及系統(tǒng)。一方面，本發(fā)明實施例提供了一種催化裂化反應(yīng)軟測量輔助變量的選取方法，包括：采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，所述工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)；根據(jù)所述工況參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)；若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；根據(jù)所述優(yōu)化模型參數(shù)和所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，獲取催化裂化反應(yīng)不同工況對應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將所述優(yōu)化產(chǎn)物分布數(shù)據(jù)作為所述催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果。進(jìn)一步地，所述利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)，包括：若所述數(shù)值模擬為初次模擬，則設(shè)置模型參數(shù)初始值，并根據(jù)不同工況對應(yīng)的所述工況參數(shù)，設(shè)定流動模型參數(shù)；根據(jù)所述模型參數(shù)初始值和所述流動模型參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)所述不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)。進(jìn)一步地，所述若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，包括：若判斷獲知所述偏差總和大于預(yù)設(shè)閾值，則利用遺傳算法對所述模型參數(shù)進(jìn)行優(yōu)化。進(jìn)一步地，所述根據(jù)所述優(yōu)化模型參數(shù)和所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，獲取催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，包括：將所述優(yōu)化模型參數(shù)代入所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，對所述催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取所述優(yōu)化產(chǎn)物分布數(shù)據(jù)。進(jìn)一步地，所述催化裂化反應(yīng)動力學(xué)模型為十四集總動力學(xué)模型。進(jìn)一步地，所述工況參數(shù)包括：溫度數(shù)據(jù)、壓力數(shù)據(jù)和流量數(shù)據(jù)。進(jìn)一步地，所述實測產(chǎn)物分布數(shù)據(jù)包括：干氣收率、液化氣收率、汽油收率、柴油收率和焦炭收率。另一方面，本發(fā)明實施例提供一種催化裂化反應(yīng)產(chǎn)物軟測量的系統(tǒng)，包括：數(shù)據(jù)采集單元，用于采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，所述工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)；計算產(chǎn)物分布單元，用于根據(jù)所述工況參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)；模型參數(shù)優(yōu)化單元，用于若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；軟測量單元，用于根據(jù)所述優(yōu)化模型參數(shù)和所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，獲取催化裂化反應(yīng)不同工況對應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將所述優(yōu)化產(chǎn)物分布數(shù)據(jù)作為所述催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果。再一方面，本發(fā)明實施例提供一種用于催化裂化反應(yīng)軟測量輔助變量選取的電子設(shè)備，包括：至少一個處理器；以及與所述處理器通信連接的至少一個存儲器，其中：所述存儲器存儲有可被所述處理器執(zhí)行的程序指令，所述處理器調(diào)用所述程序指令能夠執(zhí)行上述的方法。又一方面，本發(fā)明實施例提供一種非暫態(tài)計算機(jī)可讀存儲介質(zhì)，所述非暫態(tài)計算機(jī)可讀存儲介質(zhì)存儲計算機(jī)指令，所述計算機(jī)指令使所述計算機(jī)執(zhí)行上述方法。本發(fā)明實時例提供的催化裂化反應(yīng)產(chǎn)物軟測量的方法和系統(tǒng)，利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，并對催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，獲得滿足條件的催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，重新利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，獲得最終的獲取催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，實現(xiàn)催化反應(yīng)產(chǎn)物軟測量。將催化裂化反應(yīng)動力學(xué)模型與計算流體力學(xué)法cfd相結(jié)合，通過優(yōu)化動力學(xué)模型參數(shù)，實現(xiàn)催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)產(chǎn)物軟測量結(jié)果的精度。此外，本發(fā)明實施例在催化裂化工況發(fā)生改變時，只需要相應(yīng)的調(diào)整催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，即可獲得不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)，實現(xiàn)不同工況的催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)軟測量方法的適應(yīng)性。附圖說明為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作一簡單地介紹，顯而易見地，下面描述中的附圖是本發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動的前提下，還可以根據(jù)這些附圖獲得其他的附圖。圖1為本發(fā)明實施例中催化裂化反應(yīng)產(chǎn)物軟測量方法的流程示意圖；圖2為本發(fā)明實施例中十四集總動力學(xué)模型的反應(yīng)網(wǎng)絡(luò)示意圖；圖3為本發(fā)明實施例中又一催化裂化反應(yīng)產(chǎn)物軟測量方法的流程示意圖；圖4為本發(fā)明實施例中利用計算流體力學(xué)法進(jìn)行催化裂化反應(yīng)數(shù)值模擬的流程示意圖；圖5為本發(fā)明實施例中又一催化裂化反應(yīng)產(chǎn)物軟測量的方法的流程示意圖；圖6為本發(fā)明實施例中進(jìn)行模型參數(shù)優(yōu)化時適應(yīng)度函數(shù)的曲線示意圖；圖7為本發(fā)明實施例中催化裂化反應(yīng)產(chǎn)物軟測量系統(tǒng)的結(jié)構(gòu)示意圖；圖8為本發(fā)明實施例中用于催化裂化反應(yīng)產(chǎn)物軟測量的電子設(shè)備的結(jié)構(gòu)示意圖。具體實施方式為使本發(fā)明實施例的目的、技術(shù)方案和優(yōu)點更加清楚，下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實施例是本發(fā)明一部分實施例，而不是全部的實施例。基于本發(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護(hù)的范圍。圖1為本發(fā)明實施例中催化裂化反應(yīng)產(chǎn)物軟測量方法的流程示意圖，如圖1所示，本發(fā)明實施例提供的催化裂化反應(yīng)產(chǎn)物軟測量的方法包括：s1、采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，所述工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)；具體地，根據(jù)需要采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)。其中工況參數(shù)包括：溫度數(shù)據(jù)、壓力數(shù)據(jù)和流量數(shù)據(jù)，實測產(chǎn)物分布數(shù)據(jù)包括：干氣收率、液化氣收率、汽油收率、柴油收率和焦炭收率。工況參數(shù)可以通過工業(yè)生產(chǎn)現(xiàn)場采集獲得實時工況參數(shù)或從工廠歷史數(shù)據(jù)庫中保存的數(shù)據(jù)中獲得歷史工況參數(shù)。本發(fā)明實施例將采集到的工況參數(shù)記為數(shù)據(jù)矩陣記為c，數(shù)據(jù)矩陣中包括但不限于：1)溫度數(shù)據(jù)：氣體出口溫度、新鮮原料溫度、回?zé)捰蜏囟?、回?zé)捰蜐{溫度、霧化蒸汽溫度、預(yù)提升蒸汽溫度以及催化劑溫度；2)壓力數(shù)據(jù)：反應(yīng)器出口壓力、反應(yīng)器壓降；3)流量數(shù)據(jù)：新鮮原料流量、回?zé)捰土髁?、回?zé)捰蜐{流量、催化劑的質(zhì)量流量、預(yù)提升蒸汽量、霧化蒸汽量。實測產(chǎn)物分布數(shù)據(jù)為現(xiàn)場檢測獲得的實時實測產(chǎn)物分布數(shù)據(jù)或歷史數(shù)據(jù)庫中存儲的歷史實測產(chǎn)物分布數(shù)據(jù)。即可以通過工業(yè)生產(chǎn)現(xiàn)場采集獲得實時實測產(chǎn)物分布數(shù)據(jù)或從工廠歷史數(shù)據(jù)庫中保存的數(shù)據(jù)中獲得歷史實測產(chǎn)物分布數(shù)據(jù)，本發(fā)明實施例可以將實測產(chǎn)物分布數(shù)據(jù)記為數(shù)據(jù)矩陣記為p。其中工況數(shù)據(jù)可以通過催化裂化裝置的監(jiān)控系統(tǒng)或催化裂化dcs系統(tǒng)獲得，或通過其他方式獲取，本發(fā)明實施例不作具體限定。s2、根據(jù)所述工況參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)；具體地，獲取到催化反應(yīng)的工況數(shù)據(jù)后，根據(jù)這些工況數(shù)據(jù)利用計算流體力學(xué)法cfd(computationalfluiddynamics)對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)。即根據(jù)不同工況對應(yīng)的工況數(shù)據(jù)獲取對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)。s3、若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；具體地，將獲取到的計算產(chǎn)物分布數(shù)據(jù)和預(yù)先采集到的實測產(chǎn)物分布數(shù)據(jù)作差，判斷計算產(chǎn)物分布數(shù)據(jù)和實測產(chǎn)物分布數(shù)據(jù)的偏差總和是否大于預(yù)設(shè)閾值，若大于則對催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；若計算產(chǎn)物分布數(shù)據(jù)和實測產(chǎn)物分布數(shù)據(jù)的偏差總和小于或等于預(yù)設(shè)閾值，則將此時的催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)作為優(yōu)化模型參數(shù)。s4、根據(jù)所述優(yōu)化模型參數(shù)和所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，獲取催化裂化反應(yīng)不同工況對應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將所述優(yōu)化產(chǎn)物分布數(shù)據(jù)作為所述催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果。具體地，根據(jù)優(yōu)化模型參數(shù)和催化裂化反應(yīng)動力學(xué)模型，利用計算流體力學(xué)法cfd，獲得催化裂化反應(yīng)不同工況對應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將優(yōu)化產(chǎn)物分布數(shù)據(jù)作為催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果，實現(xiàn)催化裂化反應(yīng)不同工況下的產(chǎn)物軟測量。需要說明的是，當(dāng)催化裂化反應(yīng)的工況發(fā)生改變時，相應(yīng)的調(diào)整工況數(shù)據(jù)，以及催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，重新執(zhí)行步驟s2-s4，獲得新的工況條件下對應(yīng)的催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，實現(xiàn)不同工況的催化裂化反應(yīng)產(chǎn)物軟測量。本發(fā)明實時例提供的催化裂化反應(yīng)產(chǎn)物軟測量的方法，利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，并對催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，獲得滿足條件的催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，重新利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，獲得最終的獲取催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，實現(xiàn)催化反應(yīng)產(chǎn)物軟測量。將催化裂化反應(yīng)動力學(xué)模型與計算流體力學(xué)法cfd相結(jié)合，通過優(yōu)化動力學(xué)模型參數(shù)，實現(xiàn)催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)產(chǎn)物軟測量結(jié)果的精度。此外，本發(fā)明實施例在催化裂化工況發(fā)生改變時，只需要相應(yīng)的調(diào)整催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，即可獲得不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)，實現(xiàn)不同工況的催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)軟測量方法的適應(yīng)性。在上述實施例的基礎(chǔ)上，所述利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)，包括：若所述數(shù)值模擬為初次模擬，則設(shè)置模型參數(shù)初始值，并根據(jù)不同工況對應(yīng)的所述工況參數(shù)，設(shè)定流動模型參數(shù)；根據(jù)所述模型參數(shù)初始值和所述流動模型參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)所述不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)。具體地，在初次進(jìn)行催化裂化反應(yīng)數(shù)值模擬時，需要設(shè)置催化裂化動力學(xué)模型的模型參數(shù)初始值，具體可以根據(jù)經(jīng)驗進(jìn)行設(shè)置，或采用其他方法設(shè)置，本發(fā)明實施例不作具體限定。設(shè)置好催化裂化動力學(xué)模型的模型參數(shù)初始值后，根據(jù)不同工況對應(yīng)的工況參數(shù)，設(shè)定流動模型參數(shù)。根據(jù)模型參數(shù)初始值和流動模型參數(shù)，利用計算流體力學(xué)法cfd對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)。其中，催化裂化反應(yīng)動力學(xué)模型屬于反應(yīng)模型，實際工業(yè)反應(yīng)器中還存在著復(fù)雜的氣固流動、傳熱、傳質(zhì)，描述氣固流動、傳熱的方程叫做流動模型(因為傳熱可以看作是熱量流動、傳質(zhì)可以看作組分流動，一般在化工中統(tǒng)稱為流動模型或者傳遞模型)。流動會影響反應(yīng)，反應(yīng)也會反過來影響流動，本發(fā)明實施例將流動模型和反應(yīng)模型耦合起來，即利用計算流體力學(xué)法cfd對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，以使獲得的催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)更加準(zhǔn)確，進(jìn)一步提高催化裂化反應(yīng)產(chǎn)物軟測量的測量結(jié)果的準(zhǔn)確性。其中，根據(jù)模型參數(shù)初始值和流動模型參數(shù)，利用計算流體力學(xué)法cfd對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬的過程如下：1)建立催化裂化提升管的三維幾何模型；2)對1)中三維幾何模型進(jìn)行網(wǎng)格劃分；3)使用預(yù)選采集的催化裂化反應(yīng)的工況參數(shù)設(shè)置各入口條件、邊界條件以及流場初始條件，設(shè)定流動模型參數(shù)；4)根據(jù)催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)初始值，設(shè)定反應(yīng)模型；5)使用求解器對催化裂化提升管進(jìn)行求解，得到計算產(chǎn)物分布數(shù)據(jù)。在上述實施例的基礎(chǔ)上，所述若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，包括：若判斷獲知所述偏差總和大于預(yù)設(shè)閾值，則利用遺傳算法對所述模型參數(shù)進(jìn)行優(yōu)化。具體地，獲取到催化裂化反應(yīng)的計算產(chǎn)物分布數(shù)據(jù)后，判斷計算產(chǎn)物分布數(shù)據(jù)和實測產(chǎn)物分布數(shù)據(jù)的偏差總和是否大于預(yù)設(shè)閾值。若大于，則利用遺傳算法對此時的催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至計算產(chǎn)物分布數(shù)據(jù)和實測產(chǎn)物分布數(shù)據(jù)的偏差總和小于或等于預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)。將所述優(yōu)化模型參數(shù)代入所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，對所述催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取所述優(yōu)化產(chǎn)物分布數(shù)據(jù)。即根據(jù)優(yōu)化獲得的優(yōu)化模型模型參數(shù)，重新設(shè)定催化裂化反應(yīng)的反應(yīng)模型，利用計算流體力學(xué)法cfd對催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲得催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)。在上述實施例的基礎(chǔ)上，所述催化裂化反應(yīng)動力學(xué)模型為十四集總動力學(xué)模型。具體地，本發(fā)明實施例采用的催化裂化反應(yīng)動力學(xué)模型為十四集總動力學(xué)模型，其是一種用集總方法簡化了的反應(yīng)動力學(xué)模型，反應(yīng)體系被分為14個集總，具體劃分步驟為：首先將原料油按照其流程范圍分為柴油層(l)、回?zé)捰?蠟油層(m)和減壓渣油-油漿層(h)三層，每一層按照結(jié)構(gòu)族組成又分成烷基碳(p)、環(huán)烷碳(n)和芳香碳(a)。由于減壓渣油和油漿中存在稠環(huán)芳香烴，進(jìn)一步通過在減壓渣油和油漿中增加新的稠環(huán)芳香烴(fah)集總。最后，考慮焦炭(ck)、汽油(go)、天然氣(lpg)以及干氣(dg)，得到原料油的十四集總劃分。圖2為本發(fā)明實施例中十四集總動力學(xué)模型的反應(yīng)網(wǎng)絡(luò)示意圖，如圖2所示，十四集總動力學(xué)模型共涉及49個反應(yīng)。十四集總動力學(xué)模型涉及49個待定動力學(xué)參數(shù)，即頻率因子k，k是維度為49的矢量，十四集總動力學(xué)模型各氣相組分(集總)的化學(xué)反應(yīng)速率表達(dá)式為如下公式(1)：式中：a0——表示催化劑活性校正因子；——表示催化劑的有效活性，kh為重芳環(huán)吸附常數(shù)；φ(cc)＝(1+0.51cc)-2.78——表示催化劑結(jié)焦的影響因素；cc——表示催化劑結(jié)焦百分產(chǎn)率；kj——表示集總j的反應(yīng)速率常數(shù)，k0表示指前因子；e——表示活化能；εg——表示空隙率；ραj——表示集總j的摩爾濃度。圖3為本發(fā)明實施例中又一催化裂化反應(yīng)產(chǎn)物軟測量方法的流程示意圖，如圖3所示，下面結(jié)合某煉廠催化裂化裝置提升管反應(yīng)器實際生產(chǎn)數(shù)據(jù)，具體介紹本發(fā)明實施例的技術(shù)方案：表101工況1的催化裂化原料性質(zhì)項目單位數(shù)值油品密度[kg/m3]914.8不同溫度下的粘度mm2/s33.57(50℃)/6.493(100℃)堿性氮含量[ppmwt]728總氮含量[ppmwt]2776硫含量[％]0.5412康氏殘?zhí)縖％]3.86釩含量[ppmwt]5鎳含量[ppmwt]15鈉含量[ppmwt]10鐵含量[ppmwt]10銅含量[ppmwt]0.3表102工況1的平衡催化劑性質(zhì)項目數(shù)據(jù)活性63比表面118孔體積0.3堆比0.85碳0.0560-20μm00-40μm8.90-80μm54.80-105μm75.30-149μm92.7aps75.4鐵5410鎳3760釩1800鎳釩比2.09鈉3800銻1960銻鎳比0.52鈣6150表103工況1的提升管操作條件項目單位數(shù)值原料油處理量[t/h]152霧化蒸汽質(zhì)量流率[t/h]9.6提升介質(zhì)體積/質(zhì)量流率[nm3/h]干氣/4574富氣質(zhì)量流率[nm3/h]41472煙氣冷卻水速率[kg/h]無煙氣冷卻水物料溫度[c]220霧化蒸汽溫度[c]290提升管出口溫度[c]520±2提升介質(zhì)溫度[c]37-38汽提氣溫度[c]400富氣溫度[c]35.5空氣鼓風(fēng)機(jī)排氣溫度[c]170環(huán)境溫度[c]20物料壓力[kpa]800霧化蒸汽壓力[kpa]1000提升介質(zhì)壓力[kpa]1000汽提氣壓力[kpa]1000富氣壓力[kpa]1090環(huán)境壓力[kpa]101反應(yīng)器壓力[kpa]170再生器壓力[kpa]200密相床堆積密度[kg/m3]480物料混合溫度[c]551.4物料預(yù)熱溫度[c]220提升管混合溫度[c]550提升管出口溫度[c]520再生器密相床溫度[c]680再生器旋分溫度[c]680、705煙氣溫度[c]647催化劑循環(huán)速率[t/h]1156再生催化劑中的焦炭含量[％]0.08劑油比7-8催化劑藏量[kg]180000新鮮催化劑的補(bǔ)充速率[t/h]0.3回?zé)挶?8～20％t1、采集數(shù)據(jù)。從催化裂化裝置生產(chǎn)現(xiàn)場采集到的生產(chǎn)數(shù)據(jù)、實驗室化驗數(shù)據(jù)，即采集工況數(shù)據(jù)，主要包括工況參數(shù)c以及實測產(chǎn)物分布數(shù)據(jù)p，某工況1下參數(shù)如表101～103所示。t2、建立催化裂化反應(yīng)的十四集總動力學(xué)模型。十四集總動力學(xué)模型的具體建立方法同上述實施例一致，此處不再贅述。t3、獲取計算產(chǎn)物分布數(shù)據(jù)。建立催化裂化反應(yīng)的三維數(shù)值模擬模型，將步驟t2中建立的十四集總動力學(xué)模型代入三維數(shù)值數(shù)值模擬模型，使用計算流體力學(xué)cfd方法進(jìn)行計算，得到計算產(chǎn)物分布數(shù)據(jù)y。根據(jù)三維數(shù)值模擬模型，利用計算流體力學(xué)cfd方法獲得計算產(chǎn)物分布數(shù)據(jù)y的方法如下：1)建立催化裂化提升管的三維幾何模型；2)對催化裂化提升管的三維幾何模型進(jìn)行網(wǎng)格劃分；3)使用步驟t1中采集的工況參數(shù)c設(shè)置各入口條件、邊界條件以及流場初始條件，設(shè)定流動模型參數(shù)；4)使用步驟t2中的十四集總動力學(xué)模型及模型參數(shù)k，設(shè)定反應(yīng)模型。對于該模型參數(shù)k，若為初次計算，使用模型參數(shù)初始值進(jìn)行計算，若非初次計算，則用步驟t4中優(yōu)化后的優(yōu)化模型參數(shù)進(jìn)行計算；5)使用求解器對該提升管進(jìn)行求解，得到計算產(chǎn)物分布數(shù)據(jù)y。t4、設(shè)置適應(yīng)度函數(shù)。將計算產(chǎn)物分布數(shù)據(jù)y與步驟t1中的實測產(chǎn)物分布數(shù)據(jù)p的偏差總和定義為適應(yīng)度函數(shù)，若大于預(yù)設(shè)閾值tol，則使用遺傳算法對模型參數(shù)k進(jìn)行優(yōu)化，獲得優(yōu)化模型參數(shù)。t5、模型參數(shù)的優(yōu)化。將優(yōu)化后的優(yōu)化模型參數(shù)代入步驟t3中的三維數(shù)值模擬模型，重復(fù)t3-t4，直到適應(yīng)度函數(shù)小于或等于預(yù)設(shè)閾值tol。對于計算流體力學(xué)cfd方法求解的具體過程不在本發(fā)明實施例的敘述范圍內(nèi)，可以用商業(yè)軟件、開源軟件或自己編寫代碼進(jìn)行計算。圖4為本發(fā)明實施例中利用計算流體力學(xué)法進(jìn)行催化裂化反應(yīng)數(shù)值模擬的流程示意圖，本發(fā)明實施例中可以使用fortran語言編寫代碼計算。表2原始工況下提升管反應(yīng)器數(shù)值模擬計算結(jié)果與工業(yè)裝置數(shù)據(jù)的對比(wt％)組成干氣液化氣汽油柴油焦炭工業(yè)值3.7516.3641.6529.037.53模擬值3.6016.1340.5228.227.55表3反應(yīng)網(wǎng)絡(luò)中頻率因子對應(yīng)的指前因子和活化能表4工況改變后計算結(jié)果和模擬結(jié)果產(chǎn)物干氣液化氣汽油柴油焦炭工業(yè)實際值4.0118.4237.7530.237.78ga-cfd模擬值3.7818.4337.7430.237.30相對誤差6.08％0.05％0.03％0.00％6.50％t6、軟測量結(jié)果獲取。當(dāng)工況不發(fā)生變化時，以步驟t3的三維數(shù)值模擬模型計算產(chǎn)物分布數(shù)據(jù)y作為軟測量結(jié)果。當(dāng)工況發(fā)生變化時，調(diào)用步驟t4到步驟t5，對模型參數(shù)進(jìn)行適應(yīng)調(diào)整。圖5為本發(fā)明實施例中又一催化裂化反應(yīng)產(chǎn)物軟測量的方法的流程示意圖，如圖5所示，當(dāng)工況變化時，圖1虛線部分將會發(fā)生作用，模型參數(shù)自動進(jìn)行調(diào)整。當(dāng)操作工況發(fā)生變化時，能在短時間內(nèi)獲得校正后的模型，適應(yīng)操作工況變化。本發(fā)明實施例通過耦合計算流體力學(xué)cfd方法和遺傳算法的三維的流動反應(yīng)模型對催化裂化反應(yīng)的十四集總動力學(xué)模型的模型參數(shù)進(jìn)行校正，使軟測量結(jié)果具有良好的精度和適應(yīng)性。工況1下，計算結(jié)果如表2所示，當(dāng)工況變化時，模型參數(shù)自動進(jìn)行調(diào)整。圖6為本發(fā)明實施例中進(jìn)行模型參數(shù)優(yōu)化時適應(yīng)度函數(shù)的曲線示意圖，如圖6所示，工況改變后，利用遺傳算法進(jìn)行模型參數(shù)的調(diào)整和優(yōu)化，在經(jīng)過400次迭代后，適應(yīng)度曲線趨于平穩(wěn)，此時49個模型參數(shù)達(dá)到最優(yōu)解，即獲得優(yōu)化模型參數(shù)。將優(yōu)化獲得的優(yōu)化模型參數(shù)代入催化裂化提升管反應(yīng)器的流動反應(yīng)模型中，經(jīng)過計算得到提升管出口處的產(chǎn)品收率，并將計算值與工業(yè)實際值對比，如表4。本發(fā)明實時例提供的催化裂化反應(yīng)產(chǎn)物軟測量的方法，利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，并對催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，獲得滿足條件的催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，重新利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，獲得最終的獲取催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，實現(xiàn)催化反應(yīng)產(chǎn)物軟測量。將催化裂化反應(yīng)動力學(xué)模型與計算流體力學(xué)法cfd相結(jié)合，通過優(yōu)化動力學(xué)模型參數(shù)，實現(xiàn)催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)產(chǎn)物軟測量結(jié)果的精度。此外，本發(fā)明實施例在催化裂化工況發(fā)生改變時，只需要相應(yīng)的調(diào)整催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，即可獲得不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)，實現(xiàn)不同工況的催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)軟測量方法的適應(yīng)性。圖7為本發(fā)明實施例中催化裂化反應(yīng)產(chǎn)物軟測量系統(tǒng)的結(jié)構(gòu)示意圖，如圖7所示，本發(fā)明實施例提供的催化裂化反應(yīng)產(chǎn)物軟測量系統(tǒng)包括數(shù)據(jù)采集單元71、計算產(chǎn)物分布單元72、模型參數(shù)優(yōu)化單元73和軟測量單元74，其中：數(shù)據(jù)采集單元71用于采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，所述工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)；計算產(chǎn)物分布單元72用于根據(jù)所述工況參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)；模型參數(shù)優(yōu)化單元73用于若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；軟測量單元74用于根據(jù)所述優(yōu)化模型參數(shù)和所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，獲取催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將所述優(yōu)化產(chǎn)物分布數(shù)據(jù)作為所述催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果。具體地，數(shù)據(jù)采集單元71根據(jù)需要采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)。其中工況參數(shù)包括：溫度數(shù)據(jù)、壓力數(shù)據(jù)和流量數(shù)據(jù)，實測產(chǎn)物分布數(shù)據(jù)包括：干氣收率、液化氣收率、汽油收率、柴油收率和焦炭收率。工況參數(shù)可以通過工業(yè)生產(chǎn)現(xiàn)場采集獲得實時工況參數(shù)或從工廠歷史數(shù)據(jù)庫中保存的數(shù)據(jù)中獲得歷史工況參數(shù)?？梢酝ㄟ^工業(yè)生產(chǎn)現(xiàn)場采集獲得實時實測產(chǎn)物分布數(shù)據(jù)或從工廠歷史數(shù)據(jù)庫中保存的數(shù)據(jù)中獲得歷史實測產(chǎn)物分布數(shù)據(jù)。獲取到催化反應(yīng)的工況數(shù)據(jù)后，計算產(chǎn)物分布單元72根據(jù)這些工況數(shù)據(jù)利用計算流體力學(xué)法cfd(computationalfluiddynamics)對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)。即根據(jù)不同工況對應(yīng)的工況數(shù)據(jù)獲取對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)。模型參數(shù)優(yōu)化單元73將獲取到的計算產(chǎn)物分布數(shù)據(jù)和預(yù)先采集到的實測產(chǎn)物分布數(shù)據(jù)作差，判斷計算產(chǎn)物分布數(shù)據(jù)和實測產(chǎn)物分布數(shù)據(jù)的偏差總和是否大于預(yù)設(shè)閾值，若大于則對催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；若計算產(chǎn)物分布數(shù)據(jù)和實測產(chǎn)物分布數(shù)據(jù)的偏差總和小于或等于預(yù)設(shè)閾值，則將此時的催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)作為優(yōu)化模型參數(shù)。軟測量單元74根據(jù)優(yōu)化模型參數(shù)和催化裂化反應(yīng)動力學(xué)模型，利用計算流體力學(xué)法cfd，獲得催化裂化反應(yīng)不同工況對應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將優(yōu)化產(chǎn)物分布數(shù)據(jù)作為催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果，實現(xiàn)催化裂化反應(yīng)不同工況下的產(chǎn)物軟測量。需要說明的是，當(dāng)催化裂化反應(yīng)的工況發(fā)生改變時，相應(yīng)的調(diào)整工況數(shù)據(jù)，以及催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，將新的工況環(huán)境對應(yīng)的工況數(shù)據(jù)和催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)輸入計算產(chǎn)物分布單元72、模型參數(shù)優(yōu)化單元73和軟測量單元74，重新進(jìn)行催化裂化反應(yīng)產(chǎn)物的軟測量，獲得新的工況條件下對應(yīng)的催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，實現(xiàn)不同工況的催化裂化反應(yīng)產(chǎn)物軟測量。本發(fā)明實施例通過耦合計算流體力學(xué)cfd方法和遺傳算法的三維的流動反應(yīng)模型對催化裂化反應(yīng)的十四集總動力學(xué)模型的模型參數(shù)進(jìn)行校正，使軟測量結(jié)果具有良好的精度和適應(yīng)性。其中，本發(fā)明實施例采用的催化裂化反應(yīng)動力學(xué)模型為十四集總動力學(xué)模型，其是一種用集總方法簡化了的反應(yīng)動力學(xué)模型，該模型的建立方法同上述實施例一致，此處不再贅述。本發(fā)明實施例提供的系統(tǒng)用于執(zhí)行上述方法，其具體實施方式同上述實施例一致，此處不再贅述。本發(fā)明實時例提供的催化裂化反應(yīng)產(chǎn)物軟測量的方法和系統(tǒng)，利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，并對催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，獲得滿足條件的催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，重新利用計算流體力學(xué)法cfd進(jìn)行催化裂化反應(yīng)的數(shù)值模擬，獲得最終的獲取催化裂化反應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，實現(xiàn)催化反應(yīng)產(chǎn)物軟測量。將催化裂化反應(yīng)動力學(xué)模型與計算流體力學(xué)法cfd相結(jié)合，通過優(yōu)化動力學(xué)模型參數(shù)，實現(xiàn)催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)產(chǎn)物軟測量結(jié)果的精度。此外，本發(fā)明實施例在催化裂化工況發(fā)生改變時，只需要相應(yīng)的調(diào)整催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)，即可獲得不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)，實現(xiàn)不同工況的催化裂化反應(yīng)產(chǎn)物軟測量，提高了催化裂化反應(yīng)軟測量方法的適應(yīng)性。圖8為本發(fā)明實施例中用于催化裂化反應(yīng)產(chǎn)物軟測量的電子設(shè)備的結(jié)構(gòu)示意圖，如圖8所示，所述裝置可以包括：處理器(processor)81、存儲器(memory)82和通信總線83，其中，處理器81，存儲器82通過通信總線83完成相互間的通信。處理器81可以調(diào)用存儲器82中的邏輯指令，以執(zhí)行如下方法：采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，所述工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)；根據(jù)所述工況參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)；若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；根據(jù)所述優(yōu)化模型參數(shù)和所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，獲取催化裂化反應(yīng)不同工況對應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將所述優(yōu)化產(chǎn)物分布數(shù)據(jù)作為所述催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果。此外，上述的存儲器82中的邏輯指令可以通過軟件功能單元的形式實現(xiàn)并作為獨立的產(chǎn)品銷售或使用時，可以存儲在一個計算機(jī)可讀取存儲介質(zhì)中?；谶@樣的理解，本發(fā)明的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻(xiàn)的部分或者該技術(shù)方案的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，該計算機(jī)軟件產(chǎn)品存儲在一個存儲介質(zhì)中，包括若干指令用以使得一臺計算機(jī)設(shè)備(可以是個人計算機(jī)，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本發(fā)明各個實施例所述方法的全部或部分步驟。而前述的存儲介質(zhì)包括：u盤、移動硬盤、只讀存儲器(rom，read-onlymemory)、隨機(jī)存取存儲器(ram，randomaccessmemory)、磁碟或者光盤等各種可以存儲程序代碼的介質(zhì)。本發(fā)明實施例提供一種非暫態(tài)計算機(jī)可讀存儲介質(zhì)，所述非暫態(tài)計算機(jī)可讀存儲介質(zhì)存儲計算機(jī)指令，所述計算機(jī)指令使所述計算機(jī)執(zhí)行上述各方法實施例所提供的方法，例如包括：采集催化裂化反應(yīng)不同工況對應(yīng)的工況數(shù)據(jù)，所述工況數(shù)據(jù)包括催化裂化反應(yīng)的實測產(chǎn)物分布數(shù)據(jù)和工況參數(shù)；根據(jù)所述工況參數(shù)，利用計算流體力學(xué)法對預(yù)選建立的催化裂化反應(yīng)動力學(xué)模型進(jìn)行數(shù)值模擬，獲取催化裂化反應(yīng)不同工況對應(yīng)的計算產(chǎn)物分布數(shù)據(jù)；若判斷獲知所述計算產(chǎn)物分布數(shù)據(jù)和所述實測產(chǎn)物分布數(shù)據(jù)的偏差總和大于預(yù)設(shè)閾值，則對所述催化裂化反應(yīng)動力學(xué)模型的模型參數(shù)進(jìn)行優(yōu)化，直至所述偏差總和小于或等于所述預(yù)設(shè)閾值，獲得優(yōu)化模型參數(shù)；根據(jù)所述優(yōu)化模型參數(shù)和所述催化裂化反應(yīng)動力學(xué)模型，利用所述計算流體力學(xué)法，獲取催化裂化反應(yīng)不同工況對應(yīng)的優(yōu)化產(chǎn)物分布數(shù)據(jù)，將所述優(yōu)化產(chǎn)物分布數(shù)據(jù)作為所述催化裂化反應(yīng)產(chǎn)物軟測量的結(jié)果。以上實施例僅用于說明本發(fā)明的技術(shù)方案，而非對其限制；盡管參照前述實施例對本發(fā)明進(jìn)行了詳細(xì)的說明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解：其依然可以對前述各實施例所記載的技術(shù)方案進(jìn)行修改，或者對其中部分技術(shù)特征進(jìn)行等同替換；而這些修改或替換，并不使相應(yīng)技術(shù)方案的本質(zhì)脫離本發(fā)明各實施例技術(shù)方案的精神和范圍。當(dāng)前第1頁12