張小鋒，湛世輝，馮 霄

（中國石油大學(xué)（北京） 新能源研究院，北京 102249）

工業(yè)技術(shù)

采用減三線油回收裂解氣高溫余熱方案的研究

張小鋒，湛世輝，馮 霄

（中國石油大學(xué)（北京） 新能源研究院，北京 102249）

針對乙烯急冷過程中裂解氣420 ℃左右的熱量被直接降級至220 ℃使用，從而造成高品位熱量損失的現(xiàn)狀進(jìn)行了研究，提出了一種分等級回收裂解氣高溫余熱的新工藝流程，即引入煉廠減三線油直接噴淋來自廢熱鍋爐的高溫裂解氣，通過循環(huán)高溫減三線油去發(fā)生高壓蒸汽，與原工藝流程中急冷油循環(huán)去發(fā)生稀釋蒸汽組合，構(gòu)成分級多次回收高溫裂解氣余熱的新流程。使用Aspen Plus軟件模擬新流程和原流程，模擬結(jié)果表明，兩種流程所能提供的總熱量基本一致，新流程可以使裂解氣高溫余熱同時(shí)副產(chǎn)高壓蒸汽和稀釋蒸汽。以640 kt/a乙烯裝置為例進(jìn)行經(jīng)濟(jì)核算，可增加經(jīng)濟(jì)效益798萬元/a。

裂解氣；余熱回收；減三線油；乙烯裝置

裂解是乙烯裝置的第一道生產(chǎn)工序，能耗占整個(gè)裝置的50%～60%［1］，故降低裂解區(qū)的能耗意義重大。目前，研究者從兩個(gè)方面進(jìn)行了降低裂解區(qū)能耗的研究：一是提高裂解爐的熱效率；二是改善高溫裂解氣的熱量回收［2］。裂解氣的高溫余熱主要是指裂解氣由800～850 ℃降至165 ℃的熱量，其中，裂解氣從800～850 ℃降至接近露點(diǎn)（約420 ℃，裂解原料不同露點(diǎn)也不同）時(shí)，為了終止二次副反應(yīng)的發(fā)生，要在不足1 s內(nèi)完成冷卻，該過程需要特制的專用廢熱鍋爐才能實(shí)現(xiàn)。在該過程中，裂解氣的熱量用于發(fā)生超高壓蒸汽，超高壓蒸汽過熱后去裂解氣壓縮機(jī)和制冷壓縮機(jī)驅(qū)動(dòng)透平，能量得到合理利用。曹新波等［3］提出了將此部分余熱作為吸熱化學(xué)反應(yīng)的熱源，使余熱物理轉(zhuǎn)化為反應(yīng)產(chǎn)物所增加的化學(xué)的新方法，并將其用于裂解氣急冷與延遲焦化耦合工藝流程，提高了效率。由于焦化原料油的黏度性質(zhì)等原因，很難做到使焦化原料油與裂解氣在不足1 s的時(shí)間內(nèi)充分接觸、混合傳熱，故沒有工業(yè)裝置采用該流程。裂解氣從220 ℃到165 ℃左右的熱量現(xiàn)均用于發(fā)生稀釋蒸汽或替代低壓蒸汽直接作為加熱熱源。國內(nèi)的研究人員做過很多研究以改善這部分能量的利用效果。楊春生［4］改進(jìn)了急冷油塔系統(tǒng)的工藝流程，設(shè)置減黏塔以控制急冷油的黏度，達(dá)到發(fā)生更多稀釋蒸汽的目的，在實(shí)際應(yīng)用中取得了良好的效益。吳興松等［5］針對乙烯急冷裝置實(shí)際操作參數(shù)偏離設(shè)計(jì)值較大的現(xiàn)狀，提出了模擬組分分割、物性方法選取、單元簡化處理等一整套急冷系統(tǒng)的模擬方法，為系統(tǒng)改造提供了可靠依據(jù)，模擬結(jié)果與實(shí)際情況較吻合，提高了模擬計(jì)算熱量回收的準(zhǔn)確性。但裂解氣在420℃左右到220 ℃的范圍內(nèi)，由于涉及裂解氣的露點(diǎn)和裂解急冷油長時(shí)間循環(huán)易聚合結(jié)焦等問題，熱量無法通過換熱器傳熱被利用，關(guān)于該溫度區(qū)間熱量回收利用的研究極少。

本工作提出了采用減三線油回收裂解氣由420℃左右到220 ℃范圍內(nèi)熱量的工藝方案，使裂解氣的高溫余熱得到了更有效合理的利用。

1 乙烯裝置裂解氣急冷工藝流程

乙烯裝置裂解氣急冷工藝的流程見圖1。高溫（約850 ℃）裂解氣先經(jīng)廢熱鍋爐回收高品位熱能并副產(chǎn)超高壓蒸汽（壓力11.5 MPa），裂解氣冷卻至420 ℃；然后在急冷器中用循環(huán)急冷油直接噴淋降溫至220 ℃左右后進(jìn)入油洗塔塔釜，油洗塔塔釜急冷油少量作為裂解燃料油產(chǎn)品采出，其余先去稀釋蒸汽發(fā)生器作為熱源，以回收相當(dāng)于中壓蒸汽級位的熱量；冷卻后的急冷油分成兩股，一股送到急冷器以噴淋高溫裂解氣，另一股直接返回油洗塔中下部。

在油洗塔的中部設(shè)置有盤油循環(huán)，以回收相當(dāng)于低壓蒸汽級位的熱量；頂部用下游水洗塔冷凝下來的裂解汽油作為回流，以控制塔頂溫度在105 ℃左右；油洗塔塔頂?shù)牧呀鈿膺M(jìn)入水洗塔塔釜。

圖1 乙烯裝置裂解氣急冷工藝的流程圖Fig.1 Flowchart of quench process of ethylene plants.CG：cracking gas；LPS：low pressure steam；HPS：high pressure steam；DS：dilution steam，HTO：hydrogenation tail oil.

2 裂解氣急冷過程能量回收分析

以某640 kt/a乙烯裝置為例，在裂解氣從850 ℃至105 ℃的冷卻過程中，850 ℃至420 ℃之間的熱量用于發(fā)生超高壓蒸汽；420 ℃至220 ℃之間的熱量用于直接噴淋急冷油，循環(huán)急冷油去稀釋蒸汽發(fā)生器作為熱源，只能回收相當(dāng)于中壓蒸汽級位的熱量，造成熱量級位損失；220 ℃至185 ℃之間的熱量由循環(huán)急冷油取走，用于發(fā)生稀釋蒸汽；185 ℃至165 ℃之間的熱量由盤油取走，用于替代低壓蒸汽加熱工藝物流；165 ℃至65 ℃之間的熱量被循環(huán)急冷水取走，用于加熱工藝物流；65 ℃至40 ℃之間的熱量品質(zhì)過低，全部被循環(huán)水取走。

分析整個(gè)過程，裂解氣在420 ℃至220 ℃之間的熱量有40 MW左右降級使用，此部分熱量的品質(zhì)和總量都很可觀。文獻(xiàn)提出的方案基本是在420 ℃至220 ℃之間添加二級廢熱鍋爐和三級廢熱鍋爐，用來發(fā)生高壓蒸汽和中壓蒸汽，逐級利用熱量［6］。此類方案理論上是可行的，但由于裂解氣組成復(fù)雜，包含從H2到裂解燃料油等上百種烴類以及少量固體焦粒，420 ℃左右已接近裂解氣的露點(diǎn)溫度，如果添加二級廢熱鍋爐，裂解氣在二級廢熱鍋爐里降溫，會(huì)有少量重油冷凝出來從而結(jié)焦，使鍋爐很快失去作用，故只有以乙烷和丙烷為裂解原料的乙烯裝置采用了該技術(shù)。有些裝置采用了二級廢熱鍋爐，但裂解氣溫度仍在露點(diǎn)之上，沒有在裂解氣露點(diǎn)溫度之下使用廢熱鍋爐的例子?？梢钥紤]回收420 ℃至220 ℃之間的部分高溫?zé)崃坑脕砀碑a(chǎn)3.5 MPa的高壓蒸汽。

3 加入減三線油循環(huán)的新流程

本工作提出的方案是采用一種中間循環(huán)介質(zhì)將裂解氣420 ℃左右至220 ℃之間的熱量取走，即引入中間介質(zhì)直接噴淋來自廢熱鍋爐的高溫裂解氣，通過循環(huán)高溫中間介質(zhì)去發(fā)生高壓蒸汽，與原工藝流程中急冷油循環(huán)去發(fā)生稀釋蒸汽組合，構(gòu)成分級多次回收高溫裂解氣余熱的新流程。新工藝的流程如圖2所示。

3.1 中間介質(zhì)的選擇

對高溫裂解氣進(jìn)行噴淋冷卻的中間介質(zhì)，需滿足如下條件：1）中間介質(zhì)自身的構(gòu)成應(yīng)可以適應(yīng)操作條件下的循環(huán)運(yùn)行，外甩量適中；2）中間介質(zhì)組分不進(jìn)入或極少量進(jìn)入裂解氣分離系統(tǒng)。本工作使用Aspen Plus軟件對新工藝流程進(jìn)行模擬，中間介質(zhì)初步選擇了常三線油和減三線油，其餾程見表1。

圖2 新工藝的流程Fig.2 Flowchart of new process.

表1 常三線油與減三線油的餾程Table1 Distillation ranges of AGO and VGO

經(jīng)對比發(fā)現(xiàn)，使用常三線油，由于油品較輕，與高溫裂解氣直接噴淋混合后常三線油中大量的輕組分進(jìn)入裂解氣，影響后續(xù)分離系統(tǒng)；使用減三線油噴淋后，裂解氣中的部分高黏度重油冷凝，而減三線油只有極少一部分進(jìn)入裂解氣，對裂解氣的分離系統(tǒng)影響較小。

本模擬中使用的減三線油的餾程在399～517℃之間，這部分油品在常減壓裝置中被稱為減壓蠟油。減壓蒸餾塔可抽出減二線、減三線和減四線等餾程和組成不同的蠟油，且相鄰的油品之間有重合。不同的減壓裝置進(jìn)料不同，因此所產(chǎn)出的減壓三線蠟油也不同，但餾程基本為350～500 ℃，這類油品組成非常復(fù)雜，現(xiàn)在通常只能分析族組成和結(jié)構(gòu)族組成，分析所含鏈烷烴、環(huán)烷烴、芳烴、單環(huán)芳烴、雙環(huán)芳烴、多環(huán)芳烴和膠質(zhì)的組成，并需要借助大型儀器，如核磁共振儀、色譜和質(zhì)譜儀等，代價(jià)昂貴。減三線蠟油中普遍含烷烴和環(huán)烷烴較多，烷烴主要包括C20～36的正構(gòu)和異構(gòu)烷烴；環(huán)烷烴包括從單環(huán)到六環(huán)的帶有環(huán)戊烷環(huán)或環(huán)己烷環(huán)的環(huán)烷烴，主要以稠合類型為主；芳烴以單環(huán)、雙環(huán)和三環(huán)芳烴為主，同時(shí)含有一定量的四環(huán)及很少量高于四環(huán)的芳烴；多環(huán)芳烴也以稠合型為主。這樣的組成使減三線油中的輕組分進(jìn)入盤油，從而改變塔的熱負(fù)荷分布，但不改變塔頂裂解氣組成。急冷油中的高沸點(diǎn)餾分的組成與同樣餾程溫度的減壓蠟油有很大區(qū)別，主要為難以裂解的稠環(huán)多環(huán)芳烴（如萘類、茚類和苊類等），沸點(diǎn)溫度在500 ℃以上的組分占50%左右。這些不飽和烴類化合物在高溫下會(huì)發(fā)生縮聚反應(yīng)，生成大分子化合物，引起急冷油黏度增大。正是由于這兩種油品組成的差別，在第一次使用減三線油噴淋冷卻裂解氣的過程中，裂解氣中的重質(zhì)燃料油組分會(huì)冷凝進(jìn)入減三線循環(huán)油中，有益于改善目前乙烯裝置普遍存在的油洗塔塔釜急冷油循環(huán)過程中黏度大的問題。裂解重質(zhì)燃料油質(zhì)量約占裂解氣質(zhì)量的3%左右。使用的減三線油循環(huán)量約為裂解氣質(zhì)量的2倍，在循環(huán)過程中保持一定比例的減三線循環(huán)油的采出與補(bǔ)充更新，可以控制減三線循環(huán)油的黏度。此過程與乙烯裝置在急冷油黏度快速升高時(shí)引入調(diào)質(zhì)油降黏并加大燃料油采出的處理方法類似。

3.2 新工藝流程的模擬

根據(jù)國內(nèi)某640 kt/a乙烯裝置的設(shè)計(jì)數(shù)據(jù)，使用Aspen Plus（V7.2版本)軟件進(jìn)行工藝模擬計(jì)算，建立的新工藝流程模擬圖見圖3。

圖3 新流程急冷部分的Aspen Plus流程模擬示意圖Fig.3 Diagram of the quench process simulated by Aspen Plus. Oil-Rec：recycle oil.

該裝置有3種裂解原料，加氫尾油、石腦油和乙烷，分別使用3臺(tái)不同的裂解爐，3種原料產(chǎn)生的裂解氣組成也不同。設(shè)定3種原料產(chǎn)生的裂解氣從廢熱鍋爐出來的溫度均為420 ℃，此溫度的裂解氣在一級急冷器里與減三線油直接接觸冷卻至300℃，經(jīng)過閃蒸，液相進(jìn)入換熱器去副產(chǎn)高壓蒸汽，氣相進(jìn)入二級急冷器與急冷油直接接觸噴淋冷卻，冷卻到210 ℃后進(jìn)入油洗塔進(jìn)行分離。

新流程急冷部分的模擬結(jié)果見表2。由表2可看出，急冷需要循環(huán)油約1 050 t/h，循環(huán)油與裂解氣混合后將裂解氣冷卻到約300 ℃，而高溫循環(huán)油去換熱器發(fā)生高壓蒸汽，冷卻到260 ℃再次進(jìn)入循環(huán)。該過程給出熱量約30 MW，發(fā)生高壓蒸汽約61 t/h。

表2 新流程急冷部分的模擬結(jié)果Table 2 Calculated results of the simulated quench process

裂解氣與減三線循環(huán)油在一級急冷器中冷卻后，后續(xù)冷卻分離系統(tǒng)與原工藝一致。油洗塔塔釜循環(huán)急冷油與裂解氣在二級急冷器中噴淋冷卻至210 ℃左右進(jìn)入油洗塔，在油洗塔中繼續(xù)冷卻取熱。油洗塔底部循環(huán)急冷油取熱用來發(fā)生稀釋蒸汽，中部盤油循環(huán)取熱直接供給裝置的工藝用戶以替代低壓蒸汽，塔頂裂解氣冷卻至105 ℃左右進(jìn)入水洗塔。

使用Aspen Plus軟件對新工藝流程和原流程進(jìn)行模擬，模擬結(jié)果的對比見表3和表4。

從表3可看出，通過油洗塔塔頂去水洗塔的裂解氣流量和組成均改變很小，對后續(xù)分離系統(tǒng)幾乎沒有影響。

從表4可看出，原流程可發(fā)生稀釋蒸汽的熱量為52 MW，并提供替代低壓蒸汽熱量20 MW；而新流程可發(fā)生稀釋蒸汽的熱量僅有12 MW，并提供替代低壓蒸汽熱量28 MW，但新流程可同時(shí)提供發(fā)生高壓蒸汽的熱量30 MW。新工藝流程相當(dāng)于把原流程中用于發(fā)生稀釋蒸汽的熱量大部分用來發(fā)生了高壓蒸汽。

蒸汽品位和價(jià)格見表5。企業(yè)對稀釋蒸汽無價(jià)格核算，由于稀釋蒸汽替代的是中壓蒸汽，故取中壓蒸汽的價(jià)格174元/t，根據(jù)表4和表5的數(shù)據(jù)計(jì)算發(fā)生蒸汽的效益，原流程發(fā)生蒸汽的效益：25×102+90×174=18 210元/h；新流程發(fā)生蒸汽的效益：36×102+21×174+61×194=19 160元/h，新流程提高了高溫余熱的價(jià)值，減少了能量的降級利用損失。采用新流程可增加效益950元/h，按每年開工8 400 h計(jì)算，每年增加收益798萬元。由于本研究中所述蒸汽都是飽和蒸汽，其價(jià)格會(huì)略低些，但不影響結(jié)論的相對準(zhǔn)確性。

表3 新流程與原流程對后續(xù)系統(tǒng)的影響Table 3 Effects of the new and original processes on subsequent system

表4 新流程與原流程對油洗塔系統(tǒng)操作狀況的影響Table 4 Effects of the new and original processes on the operation parameters of oil washing tower

表5 蒸汽品位和價(jià)格Table 5 Prices of different steam grade

4 結(jié)論

1）針對乙烯裝置裂解氣余熱的回收，提出了增加減三線油對高溫裂解氣噴淋急冷，逐級降溫利用其熱量的新流程，可以合理利用高溫余熱間接發(fā)生高壓蒸汽。

2）新流程避開了裂解氣在露點(diǎn)處無法通過換熱器換熱回收高溫?zé)崃康膯栴}，對原流程改變較小，對裂解氣后續(xù)分離系統(tǒng)基本無影響，具有工業(yè)可操作性。

3）通過使用Aspen Plus軟件對急冷系統(tǒng)進(jìn)行模擬，對系統(tǒng)所產(chǎn)生的不同品位的蒸汽價(jià)值進(jìn)行核算，與原工藝流程相比，新流程可獲得顯著的經(jīng)濟(jì)效益。以640 kt/a乙烯裝置為例，每年可增加收益798萬元。

［1］ 肖雪軍，何細(xì)藕. 乙烯裝置裂解爐的節(jié)能技術(shù)［J］. 石油化工，2003，32（3）：254 - 257.

［2］ 孫曉. 乙烯裝置急冷系統(tǒng)能量回收工藝特點(diǎn)分析［J］. 乙烯工業(yè)，2003，15（4）：46 - 47.

［3］ 曹新波，湯志剛，段占庭，等. 乙烯急冷-延遲焦化集成節(jié)能流程［J］. 過程工程學(xué)報(bào)，2006，6（2）：267 - 271.

［4］ 楊春生. 乙烯裝置急冷油系統(tǒng)的設(shè)計(jì)和節(jié)能改造［J］. 中外能源，2011，16（1）：106 - 111.

［5］ 吳興松，楊春生. 乙烯裝置急冷系統(tǒng)流程模擬計(jì)算方法的探討［J］. 乙烯工業(yè)，2003，15（2）：9 - 15.

［6］ 朱懋斌，馮霄. 乙烯裝置汽油分餾塔系統(tǒng)中循環(huán)急冷油流量對產(chǎn)品的影響［J］. 計(jì)算機(jī)與應(yīng)用化學(xué)，2009，25（4）：399 -404.

（編輯 王 萍）

·最新專利文摘·

一種多產(chǎn)環(huán)己基苯的方法

該專利涉及一種多產(chǎn)環(huán)己基苯的方法，包括在含有MWW結(jié)構(gòu)分子篩催化劑存在下的苯與環(huán)己烯的接觸反應(yīng)得到環(huán)己基苯的過程Ⅰ，苯與多環(huán)己基苯接觸反應(yīng)進(jìn)一步得到環(huán)己基苯的過程Ⅱ，以及分離過程；過程Ⅰ和過程Ⅱ的產(chǎn)物經(jīng)分離后得到苯、甲基環(huán)戊基苯、環(huán)己基苯和多環(huán)己基苯，其中，分離過程得到的苯循環(huán)回用于過程Ⅰ和/或過程Ⅱ，過程Ⅱ中的多環(huán)己基苯來自分離過程。該方法具有較高的環(huán)己基苯收率。（中國石油化工股份有限公司；中國石油化工股份有限公司石油化工科學(xué)研究院）/CN 103848711 A，2014-06-11

一種脫氫催化劑及其制備方法

該專利涉及一種脫氫催化劑的制備方法，包括如下步驟：1）將氯鉑酸乙二醇溶液與含Sn氧化鋁載體混合，于60～120 ℃攪拌條件下回流1～6 h；2）將KOH乙二醇溶液加入步驟1）的體系中，升溫至150～200 ℃繼續(xù)回流5～120 min，結(jié)束后冷卻至室溫；3）分離步驟2）所得的固液混合物，經(jīng)洗滌至無氯離子后進(jìn)行真空干燥及惰性氣體氣氛下的焙燒；4）使用KNO3水溶液浸漬步驟3）所得固體，經(jīng)真空干燥、惰性氣體氣氛下焙燒后制得脫氫催化劑。該方法制備的催化劑能有效抑制SnPt合金的生成，提高貴金屬Pt的利用率及催化劑的穩(wěn)定性。（中國石油化工股份有限公司；中國石油化工股份有限公司撫順石油化工研究院）/CN 103785383 A，2014-05-14

首例可回收的熱固性樹脂新型聚合物

Chem Eng， 2014-07-01

美國IBM Almaden研究中心的研究人員發(fā)現(xiàn)了一種新的熱固性聚合物。這種聚六氫三嗪（PHT）新型材料，具有非常高的楊氏模量（高達(dá)14 GPa，與之相比典型的熱固性聚酰亞胺的楊氏模量為4 GPa），以及許多其他的特性，使其成為一種適宜的熱固性材料。但不同于傳統(tǒng)的熱固性塑料，PHT是可回收的， 當(dāng)將其放置在低pH的硫酸溶液中時(shí)，其可解聚回到可以被重新聚合的苯胺類單體。而形成PHT的反應(yīng)也很簡單，并且使用容易獲得且相對便宜的原料。IBM的團(tuán)隊(duì)認(rèn)為，該工藝可以擴(kuò)大到工業(yè)水平。

PHT的應(yīng)用包括典型的熱固性用途，如用于復(fù)合材料和電子設(shè)備，但其可回收性和強(qiáng)度大大拓展了它的潛力。盡管PHT仍處于開發(fā)階段的早期，且還沒有官方正式公布的相關(guān)信息，但在PHT熱固性樹脂方面已經(jīng)獲得了很多積極的商業(yè)利益。研究的下一步將是考察PHT背后的基本化學(xué)原理及其形成過程中的條件，并拓展新的可能的應(yīng)用領(lǐng)域。

一種丙烯酰胺聚合物及其制備方法和應(yīng)用

該專利涉及一種丙烯酰胺聚合物及其制備方法和應(yīng)用。具體步驟如下：在惰性氣體保護(hù)下，將含丙烯酰胺的單體水溶液與引發(fā)劑接觸，聚合初始溫度為-10～15 ℃，聚合時(shí)間為4～15 h，聚合過程中的最高溫度與最低溫度的溫差為10～50 ℃。采用該專利方法制備的聚丙烯酰胺聚合物的產(chǎn)品性能大幅優(yōu)于現(xiàn)有產(chǎn)品。（中國石油化工股份有限公司；中國石油化工股份有限公司北京化工研究院）/CN 103864986 A，2014-06-18

一種部分水解的丙烯酰胺聚合物及其制備方法和應(yīng)用

該專利涉及一種部分水解的丙烯酰胺聚合物及其制備方法和應(yīng)用。具體方法如下：在水解條件下，將丙烯酰胺聚合物膠體顆粒與無機(jī)堿性化合物接觸反應(yīng)，水解溫度為50～120 ℃，水解時(shí)間為0.5～6.0 h，氣相中氧氣的體積分?jǐn)?shù)為3%以下。該專利制得的丙烯酰胺聚合物具有較高的水解效率，且相對分子質(zhì)量較高。（中國石油化工股份有限公司；中國石油化工股份有限公司北京化工研究院）/CN 103864987 A，2014-06-18

A New Process of Waste Heat Recovery by Integrating Cracking Gas and Vacuum Gas Oil

Zhang Xiaofeng，Zhan Shihui，F(xiàn)eng Xiao
（Institute of New Energy，China University of Petroleum，Beijing，Beijing 102249，China）

Aimed at the unreasonable waste heat recovery of ethylene plants，namely that the heat of the cracking gas with high temperature（420 ℃) was degraded directly to lower temperature（220 ℃) during the recovery，a new process was proposed，in which vacuum gas oil was sprayed into the hot cracking gas to recover the heat of the cracking gas from waste heat boilers. In the new process，both high pressure steam and dilution steam can be generated，which results in multi-level heat recovery and reasonable energy utilization. The new and original processes were simulated by Aspen Plus software and the economic beneft was signifcant.

cracking gas； waste heat recovery； vacuum gas oil； ethylene plants

1000 - 8144（2014）10 - 1190 - 06

TQ 021.8

A

2014 - 04 - 18；［修改稿日期］ 2014 - 06 - 26。

張小鋒（1986—），男，新疆維吾爾自治區(qū)哈密市人，碩士生。聯(lián)系人：馮霄，電話 010 - 89739113，電郵 xfeng@ cup.edu.cn。

國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目（2012CB720500）。