專利名稱：一種在線優(yōu)化“德士古氣化爐”運(yùn)行工況的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種“急冷工藝”的德士古水煤漿加壓氣化工藝中“德士古氣化爐”運(yùn)行工況的控制方法，具體地說，涉及一種采用“急冷工藝”的德士古水煤漿加壓氣化工藝流程制備原料氣時(shí)，在正常工作狀態(tài)下，對其氣化爐的3個(gè)控制參數(shù)(包括從氣化爐頂部三流式燒嘴進(jìn)入氣化爐的氧氣總流量-后簡稱“氧氣流量”、通過氣化爐頂部的三流式燒嘴的中心管進(jìn)入氣化爐的氧氣流量-后簡稱“中心槍氧氣流量”和進(jìn)入氣化爐的急冷水流量-后簡稱“急冷水流量”)進(jìn)行在線優(yōu)化控制方法。
背景技術(shù)：
德士古水煤漿加壓氣化“急冷工藝”流程因其具有煤種適應(yīng)范圍較寬、設(shè)備相對簡單并易于實(shí)現(xiàn)大型化、氣化效率高、全過程污染輕微等特點(diǎn)，已成為大型甲醇生產(chǎn)廠家制備原料氣所采用的主流工藝之一。其工藝流程示意圖如圖1所示。圖1中，水煤漿經(jīng)煤漿給料泵加壓后送進(jìn)氣化爐頂部的三流式燒嘴中與空氣分離器送來的高壓氧氣，按照一定的氧煤比進(jìn)入氣化爐上部的燃燒室，在高溫高壓條件下進(jìn)行反應(yīng)，生成以氫氣(H2)、一氧化碳(CO)氣體、二氧化碳(CO2)氣體為主要成分的熱合成氣。 氧氣通過燒嘴的中心管和外環(huán)管，煤漿通過燒嘴的中環(huán)進(jìn)入氣化爐。煤中的灰份在高溫下熔融，熔渣與熱合成氣一起離開氣化爐上部的燃燒室進(jìn)入氣化爐下部的急冷室，被急冷水循環(huán)泵送來的急冷水淬冷后，將熱合成氣降溫至合適溫度。熔渣迅速固化并產(chǎn)生大量蒸汽， 被水蒸汽飽和并夾帶少量飛灰的合成氣從急冷室上部的合成氣出口排出。然后，該合成氣進(jìn)入文丘里洗滌器與來自灰水循環(huán)泵的洗滌水混合，使合成氣中所夾帶的固體顆粒尤其是細(xì)灰被完全潤濕，而后進(jìn)入碳洗塔，沿碳洗塔內(nèi)的下降管進(jìn)入塔底部的水浴中穿過水層，合成氣中摻雜的固體沉入水中，而氣體向上通過碳洗塔塔板，并在塔板上進(jìn)一步除塵后，經(jīng)碳洗塔頂部的除沫器除去其中夾帶的霧沫，即成為德士古水煤漿氣化系統(tǒng)的目標(biāo)產(chǎn)品，一般稱作“德士古合成氣”或“粗煤氣”，用于送至后續(xù)生產(chǎn)工段進(jìn)行下一步生產(chǎn)(參考文獻(xiàn)沈浚等.《合成氨》化學(xué)工業(yè)出版社，2001，p815)?！暗率抗藕铣蓺狻笔且环N混合性氣體，其有效成分是一氧化碳(CO)和氫氣(H2), 此外還含有一定數(shù)量的二氧化碳(CO2)、水蒸氣(H2O),微量的甲烷(CH4)、氮?dú)?N2)、硫化氫 (H2S)等。德士古氣化系統(tǒng)的“有效氣產(chǎn)率”是德士古氣化爐最重要的工況評判指標(biāo)和最重要的經(jīng)濟(jì)運(yùn)行指標(biāo)之一。“有效氣產(chǎn)率”是指單位質(zhì)量的標(biāo)準(zhǔn)煤經(jīng)氣化爐反應(yīng)后，生成的有效氣體(CCHH2)的產(chǎn)量。正常生產(chǎn)過程中，德士古氣化系統(tǒng)在一定的負(fù)荷(即一定的煤漿品質(zhì)和煤漿進(jìn)料流量，其中煤漿進(jìn)料流量由廠級調(diào)度人員根據(jù)全廠的裝置運(yùn)行狀況和生產(chǎn)目標(biāo)確定)下，操作人員的主要任務(wù)是根據(jù)裝置和原料情況，調(diào)整氣化爐的控制參數(shù)(氧氣流量、中心槍氧氣流量和急冷水流量)以優(yōu)化氣化爐的運(yùn)行工況，從而實(shí)現(xiàn)提高氣化系統(tǒng)的“有效氣產(chǎn)率”。正常生產(chǎn)過程中，在一定的煤漿品質(zhì)下，氣化系統(tǒng)的“有效氣產(chǎn)率”由氣化爐的進(jìn)料流量(包括煤漿流量、氧氣流量、中心槍氧氣流量、急冷水流量)、進(jìn)料溫度(包括煤漿溫
4度、氧氣溫度、急冷水溫度)、進(jìn)料壓力(包括煤漿壓力、氧氣壓力)和反應(yīng)條件(包括氣化爐爐膛溫度、爐膛壓力和煤漿燃燒狀況)共同決定。由于各種情況的制約，目前對氣化爐運(yùn)行工況的優(yōu)化主要通過調(diào)節(jié)進(jìn)入氣化爐的控制參數(shù)氧氣流量、中心槍氧氣流量、急冷水流量來實(shí)現(xiàn)。由于德士古氣化系統(tǒng)反應(yīng)機(jī)理復(fù)雜，具有嚴(yán)重的非線性、分布參數(shù)、大純滯后、時(shí)變、變量間強(qiáng)耦合、大規(guī)模以及不確定性等特點(diǎn)，至今國內(nèi)外均未見有準(zhǔn)確地用模型來計(jì)算不同生產(chǎn)條件下氣化爐的氧氣流量、中心槍氧氣流量、急冷水流量的最優(yōu)控制值的相關(guān)報(bào)導(dǎo)。操作人員往往根據(jù)自身的操作經(jīng)驗(yàn)，對當(dāng)前生產(chǎn)條件進(jìn)行判斷，并對上述3個(gè)控制參數(shù)進(jìn)行調(diào)節(jié)。不同的操作人員，對氣化爐的操作水平往往不同；而同一操作人員，也往往難以對各種各樣的生產(chǎn)條件均做出準(zhǔn)確判斷，并準(zhǔn)確調(diào)節(jié)控制參數(shù)。從而，氣化爐的運(yùn)行工況往往達(dá)不到最優(yōu)，氣化系統(tǒng)的“有效氣產(chǎn)率”也經(jīng)常上下波動(dòng)。鑒于此，如何能在正常生產(chǎn)條件時(shí)的一定的系統(tǒng)負(fù)荷下，在線優(yōu)化氣化爐的控制參數(shù)(氧氣流量、中心槍氧氣流量、急冷水流量)，從而優(yōu)化氣化爐的反應(yīng)工況，快速準(zhǔn)確地提高德士古氣化系統(tǒng)的“有效氣產(chǎn)率”就成為本發(fā)明需要解決的關(guān)鍵技術(shù)問題，也是各甲醇生產(chǎn)廠商所關(guān)注的重要指標(biāo)之一。

發(fā)明內(nèi)容
本發(fā)明的目的在于，提供一種在大規(guī)模甲醇生產(chǎn)中采用“急冷工藝”的德士古水煤漿加壓氣化工藝流程制備原料氣時(shí)，在一定的煤漿品質(zhì)和煤漿進(jìn)料流量下，通過在線優(yōu)化 “德士古氣化爐”的控制參數(shù)(如氧氣流量、中心槍氧氣流量和急冷水流量)，在線優(yōu)化“德士古氣化爐”運(yùn)行工況的方法，達(dá)到穩(wěn)定或/和提其“有效氣產(chǎn)率”之目的。本發(fā)明所述的在線優(yōu)化“德士古氣化爐”運(yùn)行工況的方法，包括如下步驟(1)數(shù)據(jù)的采集和預(yù)處理選取采用“急冷工藝”的德士古水煤漿加壓氣化工藝流程中，入爐(指“德士古氣化爐”，以下同)煤質(zhì)參數(shù)(包括水分，簡記為“M”;灰分，簡記為“A”;揮發(fā)分，簡記為“Va”; 固定碳含量，簡記為“C” ；灰熔點(diǎn)，簡記為“ST” ；共5個(gè))，入爐煤漿參數(shù)(包括煤漿濃度， 簡記為“Ce” ；煤漿流量，簡記為“Fe” ；煤漿壓力，簡記為“Pc” ；煤漿溫度，簡記為“Tc” ；共4 個(gè))，入爐氧氣參數(shù)(包括氧氣流量，簡記為“Fo”沖心槍氧氣流量，簡記為“Foe”;氧氣壓力，簡記為“Po”;氧氣溫度，簡記為“To”;共4個(gè))，入爐急冷水參數(shù)(包括急冷水流量，簡記為“Fw” ；急冷水溫度，簡記為“IV’ ；共2個(gè))，出爐合成氣參數(shù)(包括出爐合成氣流量， 簡記為"Fg"；合成氣中CO的含量，簡記為“CO” ；合成氣中H2含量，簡記為“H2” ；共3個(gè))， 共18個(gè)變量為“德士古氣化爐”工況在線優(yōu)化方法的輔助變量，其中6個(gè)變量(“M”、“A”、 “Va”、“C”、“ST” 和“Ce”)為人工分析變量，其余 12 個(gè)(“Fc”、“Pc”、“Tc”、“Fo”、“Foc”、 “P0”、“T0”、“Fw”、“l(fā)V，、“Fg”、“CO”、“H2”)為在線測量變量；采集以上18個(gè)輔助變量1個(gè)月的過程數(shù)據(jù)，應(yīng)用3 0準(zhǔn)則和一階數(shù)字濾波技術(shù)對在線測量變量數(shù)據(jù)進(jìn)行錯(cuò)誤數(shù)據(jù)的剔除和數(shù)據(jù)平滑后，從中選取穩(wěn)態(tài)工況下的對應(yīng)時(shí)間的數(shù)據(jù)組成歷史數(shù)據(jù)樣本，且歷史數(shù)據(jù)樣本不少于300組；(2)歷史工況評判和優(yōu)化區(qū)域的劃分本發(fā)明的發(fā)明人以氣化系統(tǒng)的“有效氣產(chǎn)率”為評判標(biāo)準(zhǔn)將歷史數(shù)據(jù)樣本中的歷史工況分為“優(yōu)”、“良”、“中”、“差”4個(gè)等級，利用K均值聚類方法將歷史數(shù)據(jù)樣本中的有效氣產(chǎn)率數(shù)據(jù)集分為4類(即4個(gè)等級)，獲得4個(gè)聚類中心(即4個(gè)類別的歷史數(shù)據(jù)樣本中的“有效氣產(chǎn)率”的平均值)，并根據(jù)聚類中心的值劃分“優(yōu)”、“良”、“中”、“差” 4種工況；(3)建立步驟(1)中所列“德士古氣化爐”的前15個(gè)輔助變量與出爐合成氣流量 Fg、合成氣CO含量和合成氣吐含量的3個(gè)BP (Back Propagation)神經(jīng)網(wǎng)絡(luò)關(guān)系模型采用3個(gè)結(jié)構(gòu)相同的BP神經(jīng)網(wǎng)絡(luò)分別建立“Fg”、“C0”和“H2”的BP神經(jīng)網(wǎng)絡(luò)模型；即，對由步驟(1)中獲得的歷史數(shù)據(jù)樣本中的18個(gè)輔助變量進(jìn)行歸一化，并對前15個(gè)輔助變量(除“Fg”、“C0”和"H2”外的15個(gè)變量)進(jìn)行主元分析(主元分析的目的是減小模型的規(guī)模)，然后，將8個(gè)主元變量作為3個(gè)BP神經(jīng)網(wǎng)絡(luò)的輸入變量，“Fg”、“C0”和“H2” 分別作為3個(gè)BP神經(jīng)網(wǎng)絡(luò)的輸出變量組成訓(xùn)練樣本，分別對3個(gè)BP神經(jīng)網(wǎng)絡(luò)模型進(jìn)行訓(xùn)練。(4)對工況進(jìn)行在線優(yōu)化對系統(tǒng)當(dāng)前工況進(jìn)行評判，計(jì)算當(dāng)前氣化系統(tǒng)實(shí)際的有效氣產(chǎn)率，按O)中所述方法評價(jià)系統(tǒng)當(dāng)前工況的等級(“優(yōu)”、“良”、“中”、“差”)，并根據(jù)工況等級判斷結(jié)果進(jìn)行相應(yīng)操作。如需要進(jìn)行優(yōu)化計(jì)算，則采用粒子群優(yōu)化算法求解優(yōu)化的控制參數(shù)(氧氣流量、中心槍氧氣流量、急冷水流量)，進(jìn)而實(shí)現(xiàn)“德士古氣化爐”運(yùn)行工況在線優(yōu)化。(5)優(yōu)化區(qū)域及模型的在線更新按表1的檢測頻率“滾動(dòng)”獲取、更新、保存最近1個(gè)月的步驟⑴中所述的18個(gè)輔助變量過程數(shù)據(jù)，當(dāng)步驟C3)中建立的某個(gè)BP神經(jīng)網(wǎng)絡(luò)模型的計(jì)算結(jié)果和對應(yīng)的變量測量結(jié)果的誤差超過限制時(shí)，按步驟( 中所述方法獲得新的優(yōu)化區(qū)域，同時(shí)對步驟C3)中超過限制的BP神經(jīng)網(wǎng)絡(luò)模型重新進(jìn)行訓(xùn)練，以保證優(yōu)化區(qū)域及模型不斷適應(yīng)系統(tǒng)特性的變化和生產(chǎn)工況的遷移。


圖1為本發(fā)明所述的采用“急冷工藝”的德士古水煤漿加壓氣化工藝流程示意圖。圖2為本發(fā)明所述的歷史數(shù)據(jù)樣本的獲取過程示意圖。圖3為本發(fā)明所述的Back Propagation(BP)神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練樣本的獲取過程
示意圖。圖4為本發(fā)明所述的“德士古氣化爐”的“Fg”的BP神經(jīng)網(wǎng)絡(luò)示意圖。圖5為本發(fā)明所述的“德士古氣化爐”的“C0”BP神經(jīng)網(wǎng)絡(luò)示意圖。圖6為本發(fā)明所述的“德士古氣化爐”的“H2”的BP神經(jīng)網(wǎng)絡(luò)示意圖。
具體實(shí)施例方式本發(fā)明所說的采用“急冷工藝”的德士古水煤漿加壓氣化工藝流程制備原料氣時(shí)， 系統(tǒng)在一定負(fù)荷下，在線優(yōu)化氣化爐運(yùn)行工況，提高“德士古合成氣”有效成分的產(chǎn)出量的方法，包括如下步驟(1)數(shù)據(jù)的采集和預(yù)處理(1-1)需采集數(shù)據(jù)的變量根據(jù)對采用“急冷工藝”的德士古水煤漿加壓氣化工藝流程的分析，本發(fā)明的發(fā)明人選擇了 6個(gè)人工分析變量(“M”、“A”、“Va”、“C”、“ST”和“Cc”)、12個(gè)在線測量變量(“ Fc ”、“ Pc ”、“ Tc ”、“ Fo ”、“ Foc，，、“ Po ”、“ To ”、“ Fw ”、“ Tw ”、“ Fg ”、“ CO ” 和 “ H2 ”)(請?jiān)斠姳?)。 表 權(quán)利要求
1. 一種在線優(yōu)化“德士古氣化爐”運(yùn)行工況的方法，包括如下步驟(1)選取采用“急冷工藝”的德士古水煤漿加壓氣化工藝流程中，入“德士古氣化爐” 的變量參數(shù)水分(“M”)；灰分(“A”)；揮發(fā)分(“Va”)；固定碳含量(“C”)；灰熔點(diǎn) (“ST”)，煤漿濃度(“Ce”)；煤漿流量(“Fe”)；煤漿壓力("Pc")；煤漿溫度(“Tc”)； 氧氣流量("Fo")；中心槍氧氣流量("Foe")；氧氣壓力(“Po”)；氧氣溫度(“To”)；急冷水流量(“Fw”)；急冷水溫度(“Tw”)；及出“德士古氣化爐”的變量參數(shù)出爐合成氣流量("Fg")；合成氣中CO的含量(“CO”)；和合成氣中H2含量("H2")；共18個(gè)變量參數(shù)為“德士古氣化爐”工況在線優(yōu)化方法的輔助變量；其中“M”、“A”、“Va”、“C”、“ST”和“Ce”為人工分析變量參數(shù)，其余為在線測量變量參數(shù)；采集以上18個(gè)輔助變量1個(gè)月的過程數(shù)據(jù)，應(yīng)用3 σ準(zhǔn)則和一階數(shù)字濾波技術(shù)對在線測量變量數(shù)據(jù)進(jìn)行錯(cuò)誤數(shù)據(jù)的剔除和數(shù)據(jù)平滑后，從中選取穩(wěn)態(tài)工況下的對應(yīng)時(shí)間的數(shù)據(jù)組成歷史數(shù)據(jù)，且歷史數(shù)據(jù)樣本不少于300組；(2)以氣化系統(tǒng)的“有效氣產(chǎn)率”為評判標(biāo)準(zhǔn)將歷史數(shù)據(jù)樣本中的歷史工況分為“優(yōu)”、 “良”、“中”、“差1個(gè)等級，利用K均值聚類方法將歷史數(shù)據(jù)樣本中的有效氣產(chǎn)率數(shù)據(jù)集分為 4個(gè)等級，獲得4個(gè)聚類中心，并根據(jù)聚類中心的值劃分“優(yōu)”、“良”、“中”、“差” 4種工況；(3)采用3個(gè)結(jié)構(gòu)相同的BP神經(jīng)網(wǎng)絡(luò)分別建立“Fg”、“C0”和"H2”的BP神經(jīng)網(wǎng)絡(luò)模型； 即，對由步驟(1)中獲得的歷史數(shù)據(jù)樣本中的18個(gè)輔助變量進(jìn)行歸一化，并對前15個(gè)輔助變量進(jìn)行主元分析，然后，將8個(gè)主元變量作為3個(gè)BP神經(jīng)網(wǎng)絡(luò)的輸入變量，“Fg”、“C0”和 "H2"分別作為3個(gè)BP神經(jīng)網(wǎng)絡(luò)的輸出變量組成訓(xùn)練樣本，分別對3個(gè)BP神經(jīng)網(wǎng)絡(luò)模型進(jìn)行訓(xùn)練；(4)計(jì)算當(dāng)前氣化系統(tǒng)實(shí)際的“有效氣產(chǎn)率”，按步驟O)中所述方法評價(jià)系統(tǒng)當(dāng)前工況的等級，并根據(jù)工況等級判斷結(jié)果進(jìn)行相應(yīng)操作；如需要進(jìn)行優(yōu)化計(jì)算，則采用粒子群優(yōu)化算法求解優(yōu)化的控制變量參數(shù)“你”、 “Foe”、和“Fw”，進(jìn)而實(shí)現(xiàn)“德士古氣化爐”運(yùn)行工況在線優(yōu)化；(5)按表1的檢測頻率“滾動(dòng)”獲取、更新、保存最近1個(gè)月的步驟(1)中所述的18個(gè)輔助變量過程數(shù)據(jù)表1類別序號符號中文描述單位變化范圍檢測頻率1"M"入爐煤水分wt%
每曰一次人工2"A"入爐煤灰分wt%[5,15]分3"Va"入爐煤揮發(fā)分wt%[30,50]每周一次及析變 量4"C"入爐煤固定碳含量wt%[30,60]更換煤種時(shí)5"ST"入爐煤灰熔點(diǎn)°C[1000,1400]6"Ce"入爐煤漿濃度%[55,70]每4小時(shí)一次7"Fe"入爐煤漿流量m3/h[40,80]8"Pc"入爐煤漿壓力MPa[5,7]9"Tc"入爐煤漿溫度0C[30,60]在線 測量 變量10"Fo"入爐氧氣流量m3/h[15000,35000]11 12 13 14"Foe" "Po" "To" "Fw"中心槍氧氣流量入爐氧氣壓力入爐氧氣溫度入爐急冷水流量m3/h MPa V m3/h[3000,6000] [7,8] [25,50] [250,400]在線實(shí)時(shí)測量每5秒保存一組當(dāng)前值15"Tw"入爐急冷水溫度。C[220,240]16"Fg"出爐合成氣流量m3/h[120000,250000]17"CO"合成氣CO含量V%[35,45]18"H2"合成氣H2含量V%[35,45]當(dāng)步驟(3)中建立的某個(gè)BP神經(jīng)網(wǎng)絡(luò)模型的計(jì)算結(jié)果和對應(yīng)的變量測量結(jié)果的誤差超過限制時(shí)，按步驟O)中所述方法獲得新的優(yōu)化區(qū)域，同時(shí)對步驟(3)中超過限制的BP 神經(jīng)網(wǎng)絡(luò)模型重新進(jìn)行訓(xùn)練，以保證優(yōu)化區(qū)域及模型不斷適應(yīng)系統(tǒng)特性的變化和生產(chǎn)工況的遷移。
2.如權(quán)利要求1所述的方法，其特征在于，其中步驟(3)中所建立的BP網(wǎng)絡(luò)均含有輸入層、隱含層和輸出層的三層網(wǎng)絡(luò)結(jié)構(gòu)，輸入層有8個(gè)節(jié)點(diǎn)，隱含層有8個(gè)節(jié)點(diǎn)，輸出層有1 個(gè)節(jié)點(diǎn)。
3.如權(quán)利要求1或2所述的方法，其特征在于，在步驟中計(jì)算當(dāng)前氣化系統(tǒng)實(shí)際的“有效氣產(chǎn)率”前，還包括判斷當(dāng)前裝置是否處于正?；蚍€(wěn)定狀態(tài)的步驟，具體包括如下步驟①各個(gè)測量變量均有有效的測量值，且近5分鐘時(shí)間段內(nèi)的變化范圍在該變量允許的變化范圍內(nèi)；②“Fe”和“Fo”在5分鐘時(shí)間段內(nèi)的波動(dòng)范圍均在士5%以內(nèi)；如果裝置當(dāng)前運(yùn)行狀態(tài)不正常或不穩(wěn)定，則優(yōu)化系統(tǒng)自動(dòng)切除，不對當(dāng)前工況進(jìn)行優(yōu)化。
4.如權(quán)利要求1或2所述的方法，其特征在于，在步驟(5)中，分別判斷步驟(3)中建立的“Fg”、“C0”和“H/’的BP神經(jīng)網(wǎng)絡(luò)模型計(jì)算結(jié)果與當(dāng)前測量值的相對誤差是否超限即連續(xù)10次相對誤差超過士 2%;如果有任何一個(gè)模型結(jié)果超此限，則系統(tǒng)按步驟(1)中所述方法，根據(jù)系統(tǒng)保存的過去1個(gè)月的歷史數(shù)據(jù)重新獲取歷史數(shù)據(jù)樣本，并按步驟(2)中所述方法，對歷史數(shù)據(jù)樣本中的歷史工況進(jìn)行K均值聚類，獲得新的優(yōu)化區(qū)域，同時(shí)，按步驟(3) 中所述方法，對歷史數(shù)據(jù)樣本進(jìn)行歸一化和主元分析后，對超限的BP網(wǎng)絡(luò)模型重新進(jìn)行訓(xùn)練，以獲得新的模型參數(shù)。
全文摘要
本發(fā)明涉及一種“急冷工藝”的德士古水煤漿加壓氣化工藝中“德士古氣化爐”運(yùn)行工況的控制方法。采用本發(fā)明所述的方法，能夠?qū)崿F(xiàn)在不增加生產(chǎn)成本的情況下，實(shí)時(shí)獲得“德士古氣化爐”的“有效氣產(chǎn)率”，并對當(dāng)前系統(tǒng)工況進(jìn)行評判，從而及時(shí)掌握裝置的運(yùn)行工況及生產(chǎn)狀態(tài)，并且通過優(yōu)化計(jì)算，在線調(diào)整控制參數(shù)(氧氣流量、中心槍氧氣流量、急冷水流量)，從而優(yōu)化反應(yīng)工況，快速準(zhǔn)確地提高德士古氣化系統(tǒng)“有效氣產(chǎn)率”，給生產(chǎn)廠商帶來更好的經(jīng)濟(jì)效益和更強(qiáng)大的競爭力。
文檔編號C10J3/72GK102399594SQ20111045921
公開日2012年4月4日 申請日期2011年12月30日 優(yōu)先權(quán)日2011年12月30日
發(fā)明者孫漾, 宋淑群, 崔晨, 張亞坤, 張凌波, 曹翠文, 顧幸生 申請人:兗礦國宏化工有限責(zé)任公司, 華東理工大學(xué)