郭開華 王文靜 皇甫立霞 孫標(biāo)

中山大學(xué)工學(xué)院

LNG站場防火間距及安全性分析

郭開華 王文靜 皇甫立霞 孫標(biāo)

中山大學(xué)工學(xué)院

目前我國相關(guān)防火安全規(guī)范缺管道明確規(guī)定適用于工程事故的分析模型。針對(duì)LNG站場主要危險(xiǎn)源及防火間距的問題，依據(jù)NFPA59 A—2009標(biāo)準(zhǔn)推薦的DEGADIS、LNGFire3和PoFMISE模型原理，研發(fā)了適用于我國LNG工程設(shè)計(jì)的LNG蒸氣擴(kuò)散模型和池火熱輻射模型，并對(duì)模型的準(zhǔn)確性和可靠性進(jìn)行了對(duì)比驗(yàn)證：與Burro實(shí)驗(yàn)測定值對(duì)比，所研發(fā)的LNG蒸氣擴(kuò)散模型計(jì)算結(jié)果相對(duì)誤差為19.04%，優(yōu)于原DEGADIS模型32.88%的相對(duì)誤差。將池火熱輻射模型與PoFMISE和LNG Fire3模型的標(biāo)準(zhǔn)結(jié)果進(jìn)行了對(duì)比，結(jié)果表明：其兼具兩者的優(yōu)點(diǎn)，更趨安全可靠。據(jù)此開發(fā)的“液化天然氣（LNG）站場危險(xiǎn)性分析平臺(tái)”可用于主要危險(xiǎn)源及工藝設(shè)備的防火間距及LNG站場選址、規(guī)劃和設(shè)計(jì)過程中事故后果的分析評(píng)價(jià)。

LNG站場 防火間距 LNG安全性 重氣擴(kuò)散 池火災(zāi) 危險(xiǎn)性分析平臺(tái)

中國作為新興LNG進(jìn)口國，LNG產(chǎn)業(yè)發(fā)展雖然起步較晚，但近年來發(fā)展迅猛，LNG站場日益增多，站場的安全問題隨之凸顯。根據(jù)LNG低溫、易揮發(fā)、與空氣一定比例混合后可燃和易爆的特性，對(duì)LNG站場內(nèi)主要危險(xiǎn)源及無法規(guī)定具體防火間距的情形進(jìn)行有效約束，以滿足LNG工程設(shè)計(jì)和建設(shè)的應(yīng)用要求，是當(dāng)前亟待解決的問題，也是政府及企業(yè)關(guān)注的重點(diǎn)。由于LNG的大規(guī)模應(yīng)用，國際上對(duì)其安全性格外重視，在美國，LNG是唯一由美國聯(lián)邦法規(guī)（49 CFR，193部分）對(duì)其儲(chǔ)存設(shè)施選址和施工進(jìn)行詳細(xì)具體要求的可燃物質(zhì)［1］。美國NFPA 59A—2009《液化天然氣（LNG）生產(chǎn)、儲(chǔ)存和裝運(yùn)》標(biāo)準(zhǔn)規(guī)定［2］，對(duì)于站場內(nèi)主要危險(xiǎn)源及工藝設(shè)備的防火間距，除規(guī)定最小間距外，推薦采用可靠的事故模型進(jìn)行計(jì)算，確定事故后果的危害范圍，評(píng)價(jià)LNG站場的建筑紅線（Property Line）、與站場周邊公共區(qū)域的間距以及站場設(shè)施間距等的安全性。我國目前的相關(guān)標(biāo)準(zhǔn)，如GB 50183—2004《石油天然氣工程設(shè)計(jì)防火規(guī)范》［3］，亦提到采用事故模型計(jì)算主要危險(xiǎn)源的防火間距，但未明確規(guī)定適用于我國工程的事故數(shù)學(xué)模型。事故數(shù)學(xué)模型的建立是防火安全距離分析評(píng)估的重要技術(shù)手段。建立科學(xué)、可靠的事故模型，分析站場潛在的危險(xiǎn)，可保障已建站場的安全運(yùn)行，為新建、擴(kuò)建站場的安全設(shè)計(jì)提供依據(jù)。

筆者針對(duì)LNG站場主要危險(xiǎn)源及無法規(guī)定具體防火間距的情形，根據(jù)NFPA 59A—2009標(biāo)準(zhǔn)推薦采用的DEGADIS重氣擴(kuò)散模型原理和LNGFire3模型原理，以及Po FMISE大池火修正模型原理，研發(fā)了適用于我國國情的LNG蒸氣擴(kuò)散模型和池火熱輻射模型，并開發(fā)了“液化天然氣（LNG）站場危險(xiǎn)性分析平臺(tái)”，可用于LNG站場的選址、規(guī)劃和設(shè)計(jì)過程中，主要危險(xiǎn)源及工藝設(shè)備防火間距的分析評(píng)價(jià)。

1 國內(nèi)外標(biāo)準(zhǔn)對(duì)LNG站場防火間距規(guī)定

1.1 明確規(guī)定最小防火間距

根據(jù)LNG站場的安全特性與工藝特點(diǎn)，NF－PA59A—2009標(biāo)準(zhǔn)對(duì)儲(chǔ)罐間距（Container Spacing）、工藝設(shè)備間距（Process Equipment Spacing）、氣化器間距（Vaporizer Spacing）、裝卸設(shè)施間距（Loading and Unloading Facility Spacing）等，明確給出最小防火間距［2］。

國內(nèi)標(biāo)準(zhǔn)GB 50183—2004《石油天然氣工程設(shè)計(jì)防火規(guī)范》第10章“液化天然氣站場”、GB 50028—2006《城鎮(zhèn)燃?xì)庠O(shè)計(jì)規(guī)范》第9章《液化天然氣氣站》，對(duì)總儲(chǔ)量小于3 000 m3的小型LNG站場，明確給出相關(guān)的防火間距［3－4］。

1.2 模型計(jì)算確定防火間距

對(duì)于LNG站場內(nèi)主要危險(xiǎn)源及無法規(guī)定具體防火間距的情形，如LNG站場的建筑紅線（Property Line）、與站場周邊公共區(qū)域的防火間距、確定事故性泄漏的危害范圍、站場內(nèi)部無法明確規(guī)定防火安全間距的其他設(shè)施等，國內(nèi)外標(biāo)準(zhǔn)的相關(guān)規(guī)定分述如下。

1.2.1 NFPA 59A—2009標(biāo)準(zhǔn)規(guī)定

采用數(shù)學(xué)模型模擬事故發(fā)生來計(jì)算防火間距。首先設(shè)定LNG“溢出場景”和特定的大氣環(huán)境，要求在預(yù)設(shè)條件下，確保：①輻射熱流在站場的“建筑紅線”或最近的居住區(qū)不會(huì)超過特定值（臨界值）水平；②在大氣中LNG蒸氣的濃度不會(huì)超過燃燒下限的50%（此時(shí)LNG泄漏產(chǎn)生的蒸氣云在環(huán)境中擴(kuò)散，未被點(diǎn)燃）。對(duì)于重氣擴(kuò)散模型，NFPA 59A—2009標(biāo)準(zhǔn)推薦采用國際通用模型DEGADIS重氣擴(kuò)散基本模型，對(duì)火災(zāi)熱輻射模型的選取，NFPA 59A—2009標(biāo)準(zhǔn)在推薦原有的LNGFire3模型的基礎(chǔ)上，指出對(duì)于大尺寸池火宜選用PoFMISE模型。

1.2.2 GB 50183標(biāo)準(zhǔn)規(guī)定

GB 50183標(biāo)準(zhǔn)規(guī)定，LNG站場的區(qū)域布置除滿足標(biāo)準(zhǔn)規(guī)定的防火間距外，尚應(yīng)按“國際公認(rèn)的高濃度氣體擴(kuò)散模型和液化天然氣燃燒的熱輻射計(jì)算模型”進(jìn)行校核，但并未明確給出事故數(shù)學(xué)模型。

2 液化天然氣（LNG）站場危險(xiǎn)性分析平臺(tái)研究

筆者開發(fā)的液化天然氣（LNG）站場危險(xiǎn)性分析平臺(tái)（以下簡稱“平臺(tái)”）采用Visual Basic 6.0集成開發(fā)環(huán)境進(jìn)行面向?qū)ο笤O(shè)計(jì)，通過Fortran計(jì)算語言對(duì)LNG蒸氣擴(kuò)散模型和池火熱輻射模型［4－5］進(jìn)行編程計(jì)算，并封裝成動(dòng)態(tài)鏈接庫（DLL，Dynamic Link Library），為Visual Basic 6.0調(diào)用。平臺(tái)適用于LNG站場選址、規(guī)劃和設(shè)計(jì)過程中，站場主要危險(xiǎn)源及無法規(guī)定具體防火安全間距的情形。對(duì)于LNG站場潛在的主要危險(xiǎn)事故（如LNG泄漏擴(kuò)散和LNG火災(zāi)熱輻射），利用平臺(tái)的泄漏場景設(shè)計(jì)和事故模型計(jì)算功能，能直觀地在LNG站場布局圖上顯示事故的危害范圍，方便站場設(shè)計(jì)人員以及評(píng)估人員確定防火間距。

2.1 事故計(jì)算模型

平臺(tái)包括兩種事故模型：LNG蒸氣擴(kuò)散模型和池火熱輻射模型。這兩種事故模型的開發(fā)分別基于NFPA59A—2009標(biāo)準(zhǔn)推薦采用的、國際通用的DEGADIS重氣擴(kuò)散模型原理和LNGFire3池火熱輻射模型原理，詳細(xì)模型描述參見本文參考文獻(xiàn)［5－6］。在LNG重氣擴(kuò)散模型基礎(chǔ)上，依據(jù)LNG泄漏氣化率與空氣卷吸率之間的關(guān)系并考慮地面和大氣環(huán)境對(duì)擴(kuò)散源的影響因素，建立了LNG泄漏擴(kuò)散源模型，確定擴(kuò)散源尺寸、LNG氣化率和空氣實(shí)際卷吸率，使得模型更具科學(xué)性并切合應(yīng)用實(shí)際。對(duì)于池火熱輻射模型，確定了LNGFire3模型和PoFMISE模型與池火尺寸的對(duì)應(yīng)關(guān)系，修正兩種模型火焰表面輻射力，突破模型對(duì)火焰尺寸的限制，得到適用于大尺寸范圍的池火熱輻射計(jì)算；并在圓形池火基礎(chǔ)上，建立了對(duì)矩形池火熱輻射計(jì)算的方法，使模型能夠適合于計(jì)算各種不同形狀和尺寸的池火情況。

2.2 模型驗(yàn)證

2.2.1 LNG蒸氣擴(kuò)散模型驗(yàn)證

將研發(fā)的LNG蒸氣擴(kuò)散模型分別與DEGADIS重氣擴(kuò)散基本模型計(jì)算結(jié)果和Burro實(shí)驗(yàn)［7－8］測定值作比較，來驗(yàn)證模型的準(zhǔn)確性和可靠性。初始條件選取Burro系列實(shí)驗(yàn)條件［7］，如表1所示（B3～B9為實(shí)驗(yàn)序號(hào)）。驗(yàn)證結(jié)果示于表2。研發(fā)模型計(jì)算結(jié)果對(duì)實(shí)驗(yàn)測定值的平均相對(duì)偏差為19.04%，小于DEGADIS重氣擴(kuò)散基本模型計(jì)算結(jié)果的偏差32.88%。由此可見，研發(fā)的擴(kuò)散模型對(duì)于LNG泄漏蒸氣擴(kuò)散濃度的計(jì)算是比較準(zhǔn)確和可靠的。

2.2.2 池火模型驗(yàn)證

將研發(fā)的LNG池火熱輻射模型計(jì)算結(jié)果分別與LNGFire3模型和PoFMISE模型發(fā)表的標(biāo)準(zhǔn)計(jì)算結(jié)果［9］進(jìn)行比較，來驗(yàn)證模型的可靠性，驗(yàn)證結(jié)果見表3。研發(fā)模型綜合LNGFire3固體火焰模型原理的優(yōu)點(diǎn)和Po FMISE對(duì)大型池火表面熱輻射力的修正，獲得的結(jié)果介于兩種模型之間，對(duì)于小型池火較貼近LNGFire3，大型池火貼近PoFMISE，其相對(duì)誤差約為2.5%的正偏差，說明預(yù)測結(jié)果更趨安全。由此可知，研發(fā)的池火熱輻射模型綜合了LNGFire3和PoFMISE模型的優(yōu)點(diǎn)，是合理、可靠的。本模型還可對(duì)方形池火進(jìn)行模擬，表3中同時(shí)列出相同當(dāng)量直徑的方形池火的計(jì)算結(jié)果。

表1 Burro系列實(shí)驗(yàn)初始條件表［7］

表2 蒸氣擴(kuò)散濃度計(jì)算模型驗(yàn)證表（與Burro實(shí)驗(yàn)［8］和DEGADIS模型對(duì)比） 體積分?jǐn)?shù)

表3 池火熱輻射計(jì)算模型驗(yàn)證表（與LNGFire3和PoFMISE比較）

3 應(yīng)用舉例

3.1 站場描述及安全間距分析

設(shè)某LNG站場用于城市燃?xì)夤?yīng)，站場內(nèi)有2只LNG儲(chǔ)罐。場內(nèi)儲(chǔ)罐區(qū)、氣化區(qū)、裝卸車區(qū)等生產(chǎn)區(qū)均屬于危險(xiǎn)區(qū)，其中儲(chǔ)罐儲(chǔ)存大量LNG，屬重要危險(xiǎn)源，筆者僅對(duì)裝卸車區(qū)的泄漏擴(kuò)散事故和儲(chǔ)罐區(qū)的池火事故進(jìn)行防火間距計(jì)算，其他危險(xiǎn)源防火間距計(jì)算原理相同，以此類推。利用平臺(tái)對(duì)蒸氣擴(kuò)散事故和池火熱輻射事故進(jìn)行分析計(jì)算，確定站場與周邊環(huán)境（如居民區(qū)、公共聚集地等）之間的防火間距。

3.1.1 泄漏擴(kuò)散事故場景設(shè)定

設(shè)定蒸氣擴(kuò)散事故發(fā)生在裝卸車區(qū)，分為有集液池和無集液池兩種情況進(jìn)行模擬計(jì)算。LNG從卸料口泄漏，泄漏速率0.05 m3／s，持續(xù)10 min。泄漏的LNG流入集液池（有集液池）或地面（無集液池），發(fā)生閃蒸氣化，瞬時(shí)產(chǎn)生大量蒸氣，形成低溫的重氣云團(tuán)沿下風(fēng)向擴(kuò)散。依據(jù)NFPA 59A—2009標(biāo)準(zhǔn)，結(jié)合當(dāng)?shù)貧庀髼l件，設(shè)定重氣擴(kuò)散發(fā)生時(shí)的大氣及環(huán)境條件如表4。

表4 儲(chǔ)罐溢出場景大氣環(huán)境條件表［1］

3.1.2 池火熱輻射事故場景設(shè)定

依據(jù)NFPA 59A—2009標(biāo)準(zhǔn)，設(shè)定池火發(fā)生的大氣環(huán)境條件為：風(fēng)速10 m／s、大氣溫度30℃、相對(duì)濕度70%。假定儲(chǔ)罐內(nèi)的LNG完全泄漏到圍堰，并引發(fā)池火事故，一種情況為兩個(gè)儲(chǔ)罐共用一個(gè)方形圍堰，一種情況為單個(gè)儲(chǔ)罐建立圓形圍堰。

3.2 安全距離計(jì)算結(jié)果

3.2.1 泄漏擴(kuò)散事故安全距離計(jì)算結(jié)果

NFPA 59A—2009標(biāo)準(zhǔn)規(guī)定，重氣擴(kuò)散事故中站場建筑紅線處LNG蒸氣在空氣中平均濃度不超過甲烷燃燒下限的50%（即2.5%）。利用平臺(tái)對(duì)裝卸車區(qū)的蒸氣擴(kuò)散事故進(jìn)行計(jì)算分析，事故影響范圍如圖1。由圖1－a可知，由于設(shè)置了集液池，形成的LNG溢出范圍及重氣云團(tuán)尺寸較?。▓D中圓形區(qū)域所示），甲烷體積濃度2.5%的等值線的影響范圍有限；由圖1－b可知，由于未設(shè)置集液池，形成的LNG溢出范圍及重氣云團(tuán)尺寸（圖中圓形區(qū)域所示）遠(yuǎn)大于有集液池的情況，相應(yīng)甲烷體積濃度2.5%的等值線影響范圍非常大，遠(yuǎn)遠(yuǎn)超出廠區(qū)范圍。因此，在廠區(qū)的危險(xiǎn)區(qū)域設(shè)置集液池收集泄漏的LNG，可有效控制重氣云團(tuán)尺寸，減小LNG泄漏擴(kuò)散造成的危害范圍。

3.2.2 池火熱輻射事故安全距離計(jì)算結(jié)果

圖1 重氣沿下風(fēng)向擴(kuò)散體積濃度等值線分布圖

表5 池火熱輻射事故危害標(biāo)準(zhǔn)規(guī)定表［1］

NFPA59A—2009標(biāo)準(zhǔn)對(duì)熱輻射安全值的限定如表5所示，依據(jù)標(biāo)準(zhǔn)對(duì)火災(zāi)熱輻射強(qiáng)度安全值的規(guī)定，利用平臺(tái)分別對(duì)方形和圓形圍堰內(nèi)池火熱輻射進(jìn)行計(jì)算分析，熱輻射影響范圍如圖2。由圖2－a可知，兩個(gè)儲(chǔ)罐共用一個(gè)圍堰，導(dǎo)致池火面積較大，熱輻射波及廠區(qū)周圍較大范圍；由圖2－b可知，對(duì)單個(gè)儲(chǔ)罐建立與儲(chǔ)罐外罐相當(dāng)?shù)膱A形高圍堰（如雙容罐和全容罐），可大大減小池火面積，熱輻射波及范圍明顯縮小，且避免當(dāng)只有單個(gè)儲(chǔ)罐發(fā)生泄漏而引發(fā)與它共用一個(gè)圍堰的儲(chǔ)罐發(fā)生連帶事故。因此，為減小池火事故波及范圍，在設(shè)置儲(chǔ)罐圍堰時(shí)，應(yīng)考慮盡量減小圍堰面積，即減小池火面積，必要時(shí)可在儲(chǔ)罐周圍設(shè)置高圍堰或采用雙容罐和全容罐，一方面可將火災(zāi)控制在較小范圍，另一方面可盡量減小對(duì)附近其他危險(xiǎn)源的影響，以免發(fā)生連帶事故。這也正是為什么國際上普遍認(rèn)為采用雙容罐或全容罐安全等級(jí)高的原因。

4 結(jié)論

1）計(jì)算LNG站場內(nèi)主要危險(xiǎn)源和工藝設(shè)施的防火間距，應(yīng)采用可靠的、經(jīng)實(shí)際應(yīng)用驗(yàn)證準(zhǔn)確可靠的數(shù)學(xué)模型，預(yù)測事故后果影響的范圍，確定安全距離。

2）事故數(shù)學(xué)模型的建立是防火距離安全分析評(píng)估的重要技術(shù)手段。建立科學(xué)、準(zhǔn)確、可靠的事故模型，分析站場潛在的危險(xiǎn)，可保障已建站場的安全運(yùn)行，并為新建、擬擴(kuò)建站場的安全設(shè)計(jì)提供依據(jù)。

3）筆者研發(fā)的重氣擴(kuò)散模型和池火熱輻射模型及液化天然氣（LNG）站場危險(xiǎn)性分析平臺(tái)，可用于LNG站場選址、規(guī)劃和設(shè)計(jì)過程中，站場主要危險(xiǎn)源及工藝設(shè)施防火安全間距的分析和確定。

4）為減小LNG事故影響范圍，應(yīng)在LNG站場危險(xiǎn)區(qū)域采取相應(yīng)的防護(hù)措施，如在廠區(qū)危險(xiǎn)區(qū)域或危險(xiǎn)設(shè)備處設(shè)置集液池收集泄漏的LNG；在儲(chǔ)罐區(qū)設(shè)置圍堰控制火災(zāi)面積且應(yīng)考慮盡量減小圍堰面積，必要時(shí)可在儲(chǔ)罐周圍設(shè)置高圍堰或采用雙容罐和全容罐。

［1］RAJ P K，LEMOFF T.Risk analysis based LNG facility siting standard in NFPA 59A［J］.Journal of Loss Prevention in the Process Industries，2009，22：820－829.

［2］National Fire Protection Association.Standard for the production，storage，and handling of liquefied natural gas（LNG）［S］.NFPA 59A－2009，2009.

［3］中華人民共和國建設(shè)部，中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.GB 50183—2004石油天然氣工程設(shè)計(jì)防火規(guī)范［S］.北京：中國計(jì)劃出版社，2004.

［4］中華人民共和國建設(shè)部，中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.GB 50028—2006城鎮(zhèn)燃?xì)庠O(shè)計(jì)規(guī)范［S］.北京：中國建筑工業(yè)出版社，2006.

［5］孫標(biāo)，郭開華.LNG池火熱輻射模型及安全距離影響因素研究［J］.中國安全科學(xué)學(xué)報(bào)，2010，20（9）：51－55.

［6］孫標(biāo)，郭開華.LNG重氣擴(kuò)散安全距離及影響因素［J］.天然氣工業(yè)，2010，30（7）：110－113.

［7］KOOPMAN R P，BAKER J，CEDERWALL R T，et al.LLNL／NWC 1980 LNG spill tests（Volume I）Burro series data report［R］.CA（USA）：Lawrence Livermore National Lab.，1982.

［8］MOHAN M，PANWAR T S，SINGH M P.Development of dense gas dispersion model for emergency preparedness［J］.Atmosphere Environment，1995，29（16）：2075－2087.

［9］RAJ P K.Spectrum of fires in an LNG facility assessments，models and consideration in risk evaluations：final technical report［R］.Washington DC：U.S.Department of Transportation Pineline ＆Hazardous Materials Safety Administration，2006.

Fireproof spacing between equipments and safety analysis at LNG stations

Guo Kaihua，Wang Wenjing，Huangpu Lixia，Sun Biao
（School of Engineering，Sun Yat－Sen University，Guangzhou，Guangdong 510006，China）

NATUR.GAS IND.VOLUME 32，ISSUE 10，pp.90－94，10／25／2012.（ISSN 1000－0976；In Chinese）

This paper aims to get a better understanding of major hazards and determine the unsettled fireproof distance at LNG stations.There are neither specific fire protection codes nor any accident analysis models adaptable to engineering practices in present China.In view of this，the LNG vapor dispersion and pool－fire thermal radiation models are developed especially for LNG engineering projects based on the principles of DEGADIS，LNG Fire3，and PoFMISE models recommended in NFPA59A－2009.Then，the accuracy and reliability of the developed modes are validated and compared.（1）Compared with the field experiment of Burros data，the relative error in the calculation results with the LNG vapor dispersion model is 19.04%，better than the previous DEGADIS model with the relative error of 32.88%.（2）In contrast with the results between the PoFMISE and LNG Fire3 and the pool－fire thermal radiation model，the latter combines the merits of the former two models and is more secure and reliable.On this basis，the developed Hazard ＆Risk Analysis Platform for LNG Stations can be also applied to the analysis of major hazards and the fireproof spacing between equipments and to the evaluation of the consequences after accidents in the process of site selection，planning and design of an LNG station.

LNG station，fireproof spacing，LNG safety，diffusion of heavy gas，pool fire，Hazard ＆Risk Analysis Platform

郭開華等.LNG站場防火間距及安全性分析.天然氣工業(yè)，2012，32（10）：90－94.

10.3787／j.issn.1000－0976.2012.10.022

廣東省教育廳液化天然氣與低溫技術(shù)重點(diǎn)實(shí)驗(yàn)室資助項(xiàng)目（編號(hào)：39000－3211101），中山大學(xué)－BP液化天然氣中心資助項(xiàng)目（編號(hào)：99103－9390001）。

郭開華，1956年生，教授，美國密歇根大學(xué)博士；現(xiàn)任中山大學(xué)－BP液化天然氣中心主任，從事能源工程及液化天然氣技術(shù)研究工作。地址：（510006）廣東省廣州市新港西路132號(hào)中山大學(xué)工學(xué)院。電話：（020）39332893。E－mail：guokaih＠m(xù)ail.sysu.edu.cn

（修改回稿日期 2012－09－01 編輯 趙 勤）

DOI：10.3787／j.issn.1000－0976.2012.10.022

Guo Kaihua，professor，born in 1956，got his Ph.D degree from the University of Michigan.He has been engaged in research of en－ergy engineering and LNG technologies.

Add：No.132，West Xingang Rd.，Guangzhou，Guangdong 510006，P.R.China

E－mail：guokaih＠m(xù)ail.sysu.edu.cn