弟22卷第6動(dòng)力工程2088·2012年12月POWER ENGINEERINGDeC 2(002文章編號(hào):1000-6761(2002)062088-05氣化爐自動(dòng)建模系統(tǒng)研究白泉,李政,倪維斗清華大學(xué)熱能工程系,功力機(jī)械及工程研究所北京100084)摘要:為了減輕建模人員在建模公式推導(dǎo)、邏關(guān)系推導(dǎo)和端制仿真翟序的工作量,設(shè)計(jì)了氣化爐自動(dòng)建模系統(tǒng)。該系統(tǒng)在反應(yīng)物為多種氣體、多種國(guó)體(均考慮為單篩分)的條件下,按照小室檨型的擔(dān)架自動(dòng)列寫相應(yīng)各個(gè)物質(zhì)的質(zhì)量、能量平衡方程,并自動(dòng)生成相應(yīng)的計(jì)算機(jī)程序,方使仿真計(jì)算、由于該系魷避開了工程中常用的純數(shù)值建模的思路,采用了新的數(shù)學(xué)工具一以符號(hào)運(yùn)算為特長(zhǎng)的計(jì)算機(jī)代數(shù)系統(tǒng)為計(jì)算平臺(tái)才使得該系統(tǒng)具有強(qiáng)大的邏輯關(guān)系推導(dǎo)能力。系統(tǒng)在某氣化爐的建幞過(guò)程中得到了驗(yàn)證和示范,事實(shí)表明:該系統(tǒng)大大減少了氣化爐的建樸工作量,特別是減少由于化學(xué)反應(yīng)假定的不同帶來(lái)的建模復(fù)雜性,對(duì)于建模與仿真研究新方法的突破有所貢獻(xiàn)。關(guān)鍵詞:氣化爐;數(shù)學(xué)模型;邏輯關(guān)系;符號(hào)運(yùn)算中圖分類號(hào):TK229文獻(xiàn)標(biāo)識(shí)碼:A0引言和濃度。反應(yīng)動(dòng)力學(xué)模型包含了氣化爐中發(fā)生的多種化學(xué)反應(yīng)的化學(xué)動(dòng)力學(xué)信息,如:化學(xué)反應(yīng)路氣化爐是一個(gè)具有復(fù)雜物理化學(xué)過(guò)程的熱動(dòng)徑、化學(xué)反應(yīng)速率等,從這些信息可以推算出每種力系統(tǒng)裝置,氣化過(guò)程中發(fā)生的化學(xué)反應(yīng)的多樣反應(yīng)物的生成/消耗速率流動(dòng)模型和反應(yīng)動(dòng)力學(xué)性更加深了氣化爐的建模的復(fù)雜性。為了從最大模型為小室內(nèi)物質(zhì)能量平衡方程的列寫提供了程度上減少氣化爐的建模工作,特別是減少由于邊界條件。按照小室模型的建模思路將氣化爐沿化學(xué)反應(yīng)假定的不同帶來(lái)的建模復(fù)雜性,本文設(shè)軸向劃分為若干小室(圖1)后,在每個(gè)小室中建計(jì)了基于符號(hào)運(yùn)算的氣化爐自動(dòng)建模系統(tǒng)。立了煤氣組份及固體物質(zhì)的質(zhì)量和能量平衡方1氣化爐及其建模程。通過(guò)求解這些方程,可以得到氣化爐內(nèi)溫度反應(yīng)物濃度和含碳量的分布特性德士古氣化爐的工作過(guò)程是:水煤漿通過(guò)噴嘴在高速氧氣流作用下噴人氣化爐,氧氣和霧狀水煤漿在爐內(nèi)經(jīng)歷一系列復(fù)雜的物理、化學(xué)過(guò)程上拱嘎后生成以一氧化碳、氫氣、二氧化碳和水蒸氣為主射流段要成分的混煤氣及熔渣,二者一并離開反應(yīng)區(qū)后進(jìn)人底部的急冷室分離,上直筒度本課題組對(duì)氣化爐進(jìn)行過(guò)建模研究,在小室建模方法的基礎(chǔ)上建立了氣化爐整體模型。該模型大致可以分為3部分:流動(dòng)模型,反應(yīng)動(dòng)力學(xué)模型和小室內(nèi)物質(zhì)/能量平衡方程。氣化爐流動(dòng)模型由氣化爐的給料和兒何尺寸計(jì)算出氣化爐中各個(gè)小室內(nèi)氣固兩相物質(zhì)的流速收稿日期:2001-1120修訂日期:20020322作者簡(jiǎn)介:白泉,男,博十研究生,日前,主要從事新型熱動(dòng)圖1氣化爐小室的劃分力設(shè)備的計(jì)算機(jī)建模與j傷真妍究方面的工作中國(guó)煤化工CNMHG第6期動(dòng)力工程2089建立氣化爐自動(dòng)建模系統(tǒng)的動(dòng)機(jī)應(yīng)的反應(yīng)速率到各反應(yīng)物化學(xué)速率的映射關(guān)系,能夠結(jié)合流動(dòng)模型和其它模型提供的輸入?yún)?shù)在建模過(guò)程中,建模研究者遇到的一個(gè)實(shí)際推導(dǎo)出氣體、固體的物質(zhì)、能量平衡方程表達(dá)式問題是:氣化爐中同吋發(fā)生多種化學(xué)反應(yīng),由這些它還必須能夠處理多種氣體、多種固體反應(yīng)物和化學(xué)反應(yīng)推導(dǎo)各個(gè)物質(zhì)的質(zhì)量平衡方程和能量平有無(wú)能量平衡的情況。衡方程是一個(gè)繁瑣的過(guò)程,工作量大,重復(fù)性強(qiáng)程上最常用到的計(jì)算機(jī)語(yǔ)言(如C語(yǔ)言容易出錯(cuò)建立小室中氣體、固體的物質(zhì)能量平衡 FORTRAN語(yǔ)言》只注重純數(shù)值運(yùn)算,邏輯推導(dǎo)方程后,據(jù)此編制仿真程序代碼也是一個(gè)相刈簡(jiǎn)功能比較弱,要想完成上述工作,必須突破舊的思單、重復(fù)性強(qiáng)的過(guò)程,T作量大,而且容易出錯(cuò)。有路,采用更先進(jìn)的數(shù)學(xué)工具。以符號(hào)計(jì)算邏輯關(guān)時(shí)手1推導(dǎo)、編程和檢查工作加起來(lái)竟占了所有系推導(dǎo)和復(fù)雜計(jì)算見長(zhǎng)的計(jì)算機(jī)代數(shù)系統(tǒng)在此方建模工作量的50%以上。在添加/刪除某個(gè)化學(xué)面有著獨(dú)到的優(yōu)勢(shì)。因此本文采用計(jì)算機(jī)代數(shù)系反應(yīng)時(shí),推導(dǎo)和編程工作需要重新進(jìn)行,工作量較統(tǒng)作為工作平臺(tái)。3.2氣化爐自動(dòng)建模系統(tǒng)組成部分介紹建模研究者希望能夠有套工具幫助推導(dǎo)平衡(1)輸入?yún)?shù)方程表達(dá)式、幫助進(jìn)行仿真程序的編制,這樣研究氣化爐自動(dòng)建模系統(tǒng)的輸入可以分為3個(gè)部者可以將精力集中在模型的物理、化學(xué)假定是否分:反應(yīng)動(dòng)力學(xué)部分、流動(dòng)模型計(jì)算結(jié)果和其它建合理上面,而不是讓繁項(xiàng)的推導(dǎo)和編程工作占用小室模型必需的參數(shù)(圖2)。了過(guò)多的時(shí)間。①反應(yīng)動(dòng)力學(xué)部分的輸入本文設(shè)計(jì)的氣化爐自動(dòng)建模系統(tǒng)是為了滿足與反應(yīng)動(dòng)力學(xué)部分相關(guān)的輸入是指建模時(shí)研這種需求研制的究者假設(shè)的氣化爐中發(fā)生的所有化學(xué)反應(yīng)及其化3氣化爐自動(dòng)建模系統(tǒng)研究學(xué)反應(yīng)速率(表1)。同時(shí),還需要指明其中的固體3.1氣化爐自動(dòng)建模系統(tǒng)的設(shè)計(jì)思路反應(yīng)物是單篩分還是寬篩分,是否需要進(jìn)行能量氣化爐自動(dòng)建模系統(tǒng)的設(shè)計(jì)思路是:借助計(jì)平衡計(jì)算算機(jī)對(duì)輸入的化學(xué)反應(yīng)列表和化學(xué)反應(yīng)速率表達(dá)表1氣化爐自動(dòng)建模系統(tǒng)的輸入式進(jìn)行推導(dǎo),計(jì)算出各反應(yīng)物的生成/消耗速率對(duì)應(yīng)的化學(xué)反應(yīng)的化學(xué)反應(yīng)方程式再在此基礎(chǔ)上按照小室模型的框架列寫物質(zhì)、能速率的函數(shù)名量平衡方程,然后利用該系統(tǒng)將平衡方程表達(dá)式CO+1/202=CO2生成 FORTRAN程序通過(guò)將平衡方程代碼與流2H+1/20-HO動(dòng)模型和求解算法代碼一起編譯、連接,最終生成r=rutar能夠進(jìn)行仿真計(jì)算的氣化爐整體模型(圖2)。345H2C建立的自動(dòng)建模系統(tǒng)必須能夠處理各化學(xué)反C+I/rphiO2=(2-2/rphr=a3化學(xué)反應(yīng)列寰化學(xué)反應(yīng)的反速率L選動(dòng)攜型】[其它〔XO+(2/rphi1)COC-+CO,=2C化學(xué)反應(yīng)列表6 CO+H2O-CO2+H2化學(xué)反應(yīng)的反應(yīng)速事7 CO2+H2-CO+H2Or= rk4各小室內(nèi)色個(gè)物屬的是是量8C+H O=CO+H2C+2H,=CH,按平衡方程表達(dá)10 CO+H O=CO+H2rags11 CH,,HO-CO+3H2 CH+202CO2+2H:O分析計(jì)算納果圖2氣化爐整體梗型結(jié)構(gòu)和氣化爐自動(dòng)建模系統(tǒng)的作用范H中國(guó)煤化工率計(jì)算表達(dá)式都CNMHG提供相應(yīng)的函數(shù)時(shí)佚合應(yīng)速率用獨(dú)立的2090動(dòng)力工程第22卷上 ORTRAN函數(shù)計(jì)算。②流動(dòng)模型的計(jì)算結(jié)果值得指出的是:在某些情況下化學(xué)反應(yīng)方程自動(dòng)建模系統(tǒng)需要以流動(dòng)模型的計(jì)算結(jié)果作式中的系數(shù)并不是整數(shù),而是一個(gè)函數(shù)(如表1中為輸入,因此必須提供各個(gè)小室中的流動(dòng)參數(shù),在的4號(hào)反應(yīng)),這給系數(shù)的識(shí)別帶來(lái)了一定的難輸入的時(shí)候只要給出這些參數(shù)的函數(shù)名即可度由于采用符號(hào)計(jì)算工具本系統(tǒng)能夠成功識(shí)別③其它建立小室模型必需的參數(shù)該系數(shù),并將其處理成為一個(gè)特別的反應(yīng)方程式以小室模型為構(gòu)架列寫平衡方程時(shí),還需要系數(shù)其它一些變量,如揮發(fā)份分解產(chǎn)物的引入、從各小在建模的時(shí)候,并不一定所有的反應(yīng)物都參室排出爐朦的物質(zhì)的量等。如巢要列寫能量平衡與物質(zhì)、能量平衡計(jì)算,因此,需要給定進(jìn)行計(jì)算方程,還需要給出不同氣體/固體、不同溫度下對(duì)的氣體、固體物質(zhì)的名稱輸入時(shí)需要分別給定氣應(yīng)的比焓和受熱面吸熱量等。體物質(zhì)和固體物質(zhì)的名稱和種類。表2其它的參數(shù)列表本小室中質(zhì)量加入,排出爐物質(zhì)名稱給給人、帶出溫摩爾質(zhì)梃百分比度下的比焙CH√給料中C的百分含√所有固體總紿料量受熱面吸熱/化學(xué)反應(yīng)吸熱/散熱損失放熱綜上所述,所有其它輸入量在按照上述要求在自動(dòng)建模系統(tǒng)中我們沒有繼續(xù)使用手工推規(guī)范化”后可以制成表2。按照具體情況將表填導(dǎo)的思路,而是采用系數(shù)矩陣的方法從整體的角好、編制必要的子程序后,氣化爐自動(dòng)建模系統(tǒng)就度同時(shí)處理所有化學(xué)反應(yīng)方程式及反應(yīng)速率,思可以在這3部分輸入量的基礎(chǔ)上進(jìn)行自動(dòng)推導(dǎo)各路如下:個(gè)量之間的關(guān)系,生成質(zhì)量能量平衡表達(dá)式和仿首先,將各化學(xué)反應(yīng)方程式寫成反應(yīng)方程式真程序。向量,相應(yīng)的反應(yīng)速率寫成化學(xué)反應(yīng)的反應(yīng)速率(2)反應(yīng)動(dòng)力學(xué)模型中物質(zhì)關(guān)系的自動(dòng)推導(dǎo)向量,再?gòu)幕瘜W(xué)反應(yīng)方程式向量中識(shí)別并提取化自動(dòng)建模系統(tǒng)具有的第一個(gè)重要功能是由化學(xué)反應(yīng)涉及到的所有物質(zhì)的名稱,組成反應(yīng)物名學(xué)反應(yīng)方程式列表自動(dòng)推導(dǎo)反應(yīng)方程式和各物質(zhì)稱向量。然后,根據(jù)物質(zhì)名稱向量的順序從各個(gè)化的映射關(guān)系(即針對(duì)表1的輸入量進(jìn)行處理),生學(xué)反應(yīng)方程式向量中提取該物質(zhì)在各個(gè)反應(yīng)中的成各個(gè)反應(yīng)物由化學(xué)反應(yīng)引起的生成/消耗速率系數(shù)值(根據(jù)生成和消耗的不同將系數(shù)值處理成表達(dá)式正值或者負(fù)值),組成化學(xué)反應(yīng)系數(shù)矩陣。此系數(shù)研究人員手工推導(dǎo)時(shí),是按照同一種反應(yīng)物矩陣乘以反應(yīng)物名稱向量后能推出化學(xué)反應(yīng)方程的順序進(jìn)行的。過(guò)程是:先選定某種反應(yīng)物,再?gòu)氖较蛄?。此系?shù)矩陣轉(zhuǎn)置后乘以化學(xué)反應(yīng)的反應(yīng)多個(gè)化學(xué)反應(yīng)方程式中分布提取該物質(zhì)在各反應(yīng)速率向量后,就能得到綜合考慮所有化學(xué)反應(yīng)時(shí)中的系數(shù)然后乘以該反應(yīng)的反應(yīng)速率,這樣得到各個(gè)物質(zhì)的反應(yīng)速率向量。下面是一個(gè)用3個(gè)化該物質(zhì)的生成/消耗速率這種思路需要重復(fù)搜索學(xué)反應(yīng)說(shuō)明該過(guò)程的簡(jiǎn)單示例:多個(gè)化學(xué)反應(yīng)方程式,不但繁瑣、而且比較慢。中國(guó)煤化工CNMHG第6期動(dòng)力工程91C tO,cO-=C+c),zo系數(shù)矩陣轉(zhuǎn)置RO,+10RCO2CO+O2-2O2借助計(jì)算機(jī)代數(shù)系統(tǒng)的多項(xiàng)式符號(hào)運(yùn)算功則是:能,可以從化學(xué)反應(yīng)方程式向量中分離出物質(zhì)名外部供入該小室的質(zhì)量+上…個(gè)小室流入本稱向量和系數(shù)矩陣,借助符號(hào)運(yùn)算和矩陣運(yùn)算功小室的質(zhì)量一離開本小室的該物質(zhì)的量-本小室能,可以完成系數(shù)矩陣和表達(dá)式向量的邏輯與數(shù)流入到下一個(gè)小室的質(zhì)量+化學(xué)反應(yīng)生成/消耗值關(guān)系推導(dǎo),最終得到各個(gè)物質(zhì)由于化學(xué)反應(yīng)引的質(zhì)量=0起的生成/消耗速率表達(dá)式列寫小室的能量平衡方程式需要同時(shí)考慮所氣化爐自動(dòng)建模系統(tǒng)中考慮了8種氣體物質(zhì)有進(jìn)人、流出本小室的所有氣體、固體物質(zhì)攜帶的(氧氣、甲烷、一氧化碳、二氧化碳、硫化氫、氫氣,能量,同時(shí),再加上化學(xué)反應(yīng)產(chǎn)生/吸收的熱量、小氮?dú)?、水蒸?、』種固體物質(zhì)(碳)的物質(zhì)平衡。以室中所有受熱面的吸熱放熱量和散熱損失,其規(guī)表1中的輸入為例,9種物質(zhì)在發(fā)生12種化學(xué)反則是:應(yīng)的情況下,計(jì)算結(jié)果如圖3所示:外部供入該小室的所有物質(zhì)帶來(lái)的能量+上1. 165rvtat 2-2xk6個(gè)小室流入本小室的所有物質(zhì)帶來(lái)的能量離開本小室的所有物質(zhì)帶走的能量本小室流入到ch=a2-k5-xk6下一個(gè)小室的所有物質(zhì)帶走的能量+化學(xué)反應(yīng)產(chǎn)c=1+24-513(2-2}-xk1-xk3+xk4+太生/吸收的熱量+受熱面的吸熱/放熱量-散熱損co2"a1+A3+a3(-1+rphi! rvtar + xk1+ xk3在化學(xué)反應(yīng)生成/消耗速率表達(dá)式的基礎(chǔ)上結(jié)合流動(dòng)和其它輸入,自動(dòng)建模系統(tǒng)可以按照上述規(guī)則自動(dòng)選擇合適的參數(shù),“搭建”生成各小室h2. a1-2a2+a5-xk2+xk3-xk4+-3xk5內(nèi)各個(gè)物質(zhì)的質(zhì)量平衡表達(dá)式和小室能量平衡表b2o=--al-n5+0.34Srviar +xk2-xk3+ xk4-xk5+達(dá)式的形式。搭建過(guò)程中也需要針對(duì)物質(zhì)名稱問量和其它參數(shù)進(jìn)行邏輯關(guān)系推導(dǎo)。例如:在能量平衡方程式圖3化學(xué)動(dòng)力學(xué)筷型自模的結(jié)果中,自動(dòng)建模系統(tǒng)黹要根據(jù)某物質(zhì)流的物質(zhì)類別(3)建立小室內(nèi)物質(zhì)、能量平衡方程及表2中相應(yīng)屬性自動(dòng)將流量和比焓相乘生成流自動(dòng)建模系統(tǒng)完成的第二個(gè)重要功能是自動(dòng)動(dòng)物質(zhì)總焓值的表達(dá)式建立小室內(nèi)部各物質(zhì)的質(zhì)量能量平衡表達(dá)式。(4)計(jì)算機(jī)代碼的生成氣化爐某小室內(nèi)部的氣體、固體平衡關(guān)系示以前研究人員在推導(dǎo)得到物質(zhì)、能量平衡方意圖如圖4所示。在靜態(tài)條件下,針對(duì)某一個(gè)小室程的表達(dá)式后,需要手工將這些表達(dá)式編制成仿列寫的某一種氣體或者固體物質(zhì)的質(zhì)量平衡的規(guī)真程序代碼,才能進(jìn)行仿真計(jì)算本氣化爐自動(dòng)建模系統(tǒng)不但有自動(dòng)生成物上個(gè)小室的上個(gè)小室的質(zhì)、能量平衡表達(dá)式的功能,而且能夠自動(dòng)在表達(dá)氣體茂入回體入式的基礎(chǔ)上建立 FORTRAN語(yǔ)言程序代碼,在該代碼基礎(chǔ)上經(jīng)過(guò)部分加工后可以將這些代碼生成標(biāo)準(zhǔn)的 FORTRAN語(yǔ)言子程序模塊。與氣化爐其產(chǎn)生的揮發(fā)分它部分(如流動(dòng)模型和數(shù)學(xué)求解模型)放在一起編譯連接后,就能夠生成氣化爐整體模型的仿真程化學(xué)反應(yīng)序小室的本小室的生成耗中國(guó)煤化工成的第三個(gè)重要?dú)怏w霞比四體出功圖4小室模型示意圖HCNMH灬能夠根據(jù)輸入2092動(dòng)力工程第22卷的化學(xué)反應(yīng)列長(zhǎng)和各反廠的反應(yīng)速率自動(dòng)推導(dǎo)各學(xué)反應(yīng)假定提供了極人的方便,節(jié)省了大量推導(dǎo)、個(gè)反應(yīng)物的生成∥消耗速率,再在此基礎(chǔ)上按照小編程的工作量,節(jié)約了大量時(shí)間室模型的思路結(jié)合流動(dòng)等其它參數(shù)搭建多種氣事實(shí)表明:自動(dòng)建模系統(tǒng)的建立能夠使建模體、固體物質(zhì)的質(zhì)量能量平衡方程長(zhǎng)達(dá)式,最終研究人員減少在公式推演邏輯關(guān)系推導(dǎo)和編寫生成計(jì)算機(jī)髹序代碼。此代碼與其它部分結(jié)合組仿真程序等簡(jiǎn)單、低水平的細(xì)節(jié)工作上的投入,轉(zhuǎn)成氣化爐整體仿真程序而將更多的精力投入到模型的物理化學(xué)背景及模4自動(dòng)建模系統(tǒng)的驗(yàn)證型假定本身的正確與否上,有較大的實(shí)用價(jià)值參考文獻(xiàn)木文針對(duì)原有的德士古氣化爐整體模型進(jìn)行了驗(yàn)證和示范該模型將整個(gè)爐膛劃分為15個(gè)小{1李政,王人驕,等.1煤氣化爐數(shù)學(xué)模型研究一-建便室,每個(gè)小室中考慮了8種氣體、2種固體的質(zhì)量部分門].動(dòng)力工程,2001,21(2[2]李政,+天驕,等. Texaco煤氣化爐數(shù)學(xué)模型研究-計(jì)算平衡和小室能量平衡結(jié)果及分析[Jl動(dòng)力工程,2001,21(4)實(shí)踐表明:自動(dòng)建模系統(tǒng)生成的程序和手工131白泉李政,倪維斗,計(jì)算機(jī)代數(shù)在動(dòng)力系統(tǒng)建模與仿真編寫的程序計(jì)算結(jié)果一致,這驗(yàn)證了該系統(tǒng)推導(dǎo)中的應(yīng)用]動(dòng)力工程,2001,21(5):1469~147結(jié)果的正確性。在刪除某化學(xué)反應(yīng)時(shí),手工推導(dǎo)至4]SvAs. Caster M. Autonatic Diffcrentiation and imp少要花費(fèi)一天甚至幾天生成新的模型計(jì)算程序,thiation of Thermodynamic: Models Using a Computer Al自動(dòng)建模系統(tǒng)可以在幾秒鐘內(nèi)完成。gebra SystetmLC' L. European Symposium on Computer Aided5結(jié)論L時(shí)」陳宗海,過(guò)程系統(tǒng)建模與仿真「M],中國(guó)科學(xué)技術(shù)大學(xué)出版氣化爐自動(dòng)建模系統(tǒng)的使用不但能夠大大減6] Abraham1,bter.利用符號(hào)運(yùn)算進(jìn)行工程分析一一個(gè)突少該氣化爐模型第一次建模中遇到的繁瑣的表達(dá)破門.力學(xué)進(jìn)展,994,May25式推導(dǎo)工作,而且為模型建立以后修改模型的化【7黃新生,黃圳主,朱小謙多剛體系統(tǒng)計(jì)真機(jī)代數(shù)動(dòng)力學(xué)壁模研究|J1國(guó)防科技大學(xué)學(xué)報(bào),1998,02Study of Gasifier Auto-Modeling SystemBAI Quan, Li Zheng, Ni Wei-dou(Depl. of Thermal Engrg., Qinghua Univ, Beijing 100084, China)Abstract: In order to reduce the modeling work of model deduction logic relationship ircatment andwriting simulation programs in gasiler modeling procedure, a gasifier auto-modeling system is estab-lished. Such a system can deduct mass and cnergy balance equations based on cell model modelingmethod for multi gas and solid substances, and generate corresponding computer subroutine programIn contrast to pure numerical modeling method, computer algebra system, which is characterized assymbolic calculation, is used as basis. Demonstration and validation work showed the auto-modelingsystem could substitute hurnan's work to simplify gasifier modeling procedure, especially to reducedthe modeling complexity brought by the variation of chenical reaction assumption. Moreover,thegasifier auto-modeling system is a new sample of how to use computer to substitute human's modelingrk. Figs 4, tables 2 and rels 7.Key words: gasifier; mathematical model logic relationship; symbolic calculation中國(guó)煤化工CNMHG