周　斌， 黃國靜， 魯昌濤， 阮　航

(華南師范大學(xué)華南先進(jìn)光電子研究院，廣州 510006)

分布式光纖傳感中用于快速檢測的軟硬件設(shè)計(jì)

周斌*， 黃國靜， 魯昌濤， 阮航

(華南師范大學(xué)華南先進(jìn)光電子研究院，廣州 510006)

針對目前分布式光纖傳感系統(tǒng)檢測速度慢的缺點(diǎn)，利用可編程邏輯器件(FPGA)芯片為核心設(shè)計(jì)能快速檢測的數(shù)據(jù)采集系統(tǒng)，將1次全量程檢測時(shí)間從30 s甚至幾分鐘縮至10 s以內(nèi)，使得分布式光纖傳感接近實(shí)時(shí). 分布式光纖傳感技術(shù)是光纖傳感中最有應(yīng)用前景的技術(shù)之一，其最大的優(yōu)點(diǎn)是能實(shí)現(xiàn)大范圍超長距離檢測.介紹布里淵和拉曼這兩類分布式傳感的傳感原理以及其重復(fù)性、幀結(jié)構(gòu)的傳感數(shù)據(jù)特點(diǎn). 根據(jù)該特點(diǎn)設(shè)計(jì)基于FPGA內(nèi)部的多個(gè)先入先出(FIFO)的超長環(huán)形隊(duì)列進(jìn)行數(shù)據(jù)緩沖和求均算法操作，對信號進(jìn)行降噪處理，并設(shè)計(jì)基于USB2.0協(xié)議的CY7C68013A數(shù)據(jù)傳輸模塊將數(shù)據(jù)傳輸?shù)缴衔粰C(jī)顯示和存儲. 結(jié)果表明，該設(shè)計(jì)解決了分布式光纖傳感中因測量時(shí)間長，不能實(shí)時(shí)檢測的問題，實(shí)現(xiàn)準(zhǔn)實(shí)時(shí)檢測，大大提高分布式光纖傳感的性能，可應(yīng)用于需要準(zhǔn)實(shí)時(shí)監(jiān)測和大量傳感器的場合.

分布式光纖傳感； 數(shù)據(jù)采集； 光纖布里淵效應(yīng)； 光纖拉曼效應(yīng)

自低損耗光纖發(fā)明以來，光纖傳感得到迅猛發(fā)展，并應(yīng)用于各行各業(yè)中，如生物醫(yī)學(xué)[1-2]、建工[3-4]、能源[5]、交通[6-7]等領(lǐng)域. 使用分布式光纖傳感技術(shù)，可以讓光纖本身變成傳感器，每一段光纖都能檢測其沿線的物理量，非常適合需要海量傳感器的場合. 分布式光纖傳感利用光纖非彈性散射效應(yīng)，包括光纖布里淵效應(yīng)[8-9]和光纖拉曼效應(yīng)[10-11]，這些非彈性散射光的強(qiáng)度都非常弱，為了提高信噪比，需要將采集到的信號做幾千到上萬次的平均操作[12]. 目前多數(shù)分布式光纖傳感的數(shù)據(jù)采集系統(tǒng)采用商用設(shè)備和標(biāo)準(zhǔn)采集流程，數(shù)據(jù)的獲取和降噪平均算法. 因此1次全量程的傳感一般需要幾分鐘甚至更長時(shí)間，導(dǎo)致傳感過程的非實(shí)時(shí)，只能檢測到緩變的信號，很大程度上限制了分布式光纖傳感的應(yīng)用.

本文根據(jù)布里淵和拉曼這2類分布式傳感的特點(diǎn)，以可編程邏輯器件(FPGA)芯片為核心，專門針對分布式光纖傳感中1次全量程檢測耗時(shí)的問題，設(shè)計(jì)了能實(shí)現(xiàn)快速檢測的數(shù)據(jù)采集模塊系統(tǒng). 在數(shù)據(jù)采集架構(gòu)設(shè)計(jì)中，數(shù)據(jù)的獲取和降噪平均算法同時(shí)進(jìn)行；選用超大片內(nèi)RAM空間的FPGA芯片，采用超長隊(duì)列，將每組數(shù)據(jù)采集的時(shí)間間隔壓縮至最短；采用乒乓操作，將數(shù)據(jù)采集和數(shù)據(jù)上傳這2個(gè)流程獨(dú)立，將每2次全量程采集時(shí)間間隔壓縮至最短. 利用這個(gè)架構(gòu)，可以將1次完整的全量程傳感所需要的時(shí)間控制在10 s之內(nèi)，實(shí)現(xiàn)準(zhǔn)實(shí)時(shí)檢測，大大擴(kuò)大了分布式光纖傳感的應(yīng)用范圍.

FPGA采用并行的執(zhí)行模式、開發(fā)周期短、響應(yīng)時(shí)間快、穩(wěn)定性好，易于長期維護(hù)升級. 系統(tǒng)采用FPGA為主控芯片，利用內(nèi)部鎖相環(huán)提高工作頻率，控制AD信號的采集，利用超長隊(duì)列進(jìn)行數(shù)據(jù)緩存和處理，利用USB2.0協(xié)議將數(shù)據(jù)傳到上位機(jī)，達(dá)到準(zhǔn)實(shí)時(shí)監(jiān)測的目的.

1　研究方法

1.1分布式光纖傳感原理以及其傳感數(shù)據(jù)的特點(diǎn)

分布式光纖傳感包括拉曼傳感和布里淵傳感，目前普遍將時(shí)域OTDR(Optical Time Domain Reflectometer)[13]技術(shù)應(yīng)用于這2類分布式光纖傳感，分別稱為ROTDR(Raman OTDR)[14]和BOTDR(Brillouin OTDR)[15]. 其傳感數(shù)據(jù)形式非常類似，因此其數(shù)據(jù)采集模塊的架構(gòu)統(tǒng)一. 圖1分別是分布式光纖拉曼傳感和分布式光纖布里淵傳感的原理圖.

DAQ：數(shù)據(jù)采集模塊；AD：模數(shù)轉(zhuǎn)化；AOM:聲光調(diào)制器；WDM：波分復(fù)用器；APD：雪崩光電二極管；PD：光電二極管；EOM：電光調(diào)制器；PS：擾偏器.

圖1分布式光纖傳感技術(shù)的原理

Figure 1Principle of distributed optical fiber sensing technologies

聲光調(diào)制器(AOM)用于產(chǎn)生脈沖激光，波分復(fù)用器(WDM)用于分離拉曼散射的Stocks和anti-Stocks光，電光調(diào)制器(EOM)和濾波器用于產(chǎn)生泵浦光，擾偏器用于打亂光纖內(nèi)的偏振. 在上述系統(tǒng)中，探測光被調(diào)制成脈沖形式，為了實(shí)現(xiàn)高空間分辨率(如1 m)，脈沖光的脈寬通常在10 ns量級. 經(jīng)過傳感光纖后，信號光(背向Raman散射的Stocks和anti-Stocks光，以及被受激Brillouin放大的探針光)是連續(xù)的，經(jīng)過模數(shù)轉(zhuǎn)換之后，數(shù)據(jù)采集模塊采集這些時(shí)域上連續(xù)的信號光. 無論基于光纖Ramam效應(yīng)還是基于Brillouin效應(yīng)的分布式光纖傳感系統(tǒng)，采集的數(shù)據(jù)形式均為重復(fù)信號，定義為幀信號，每1次重復(fù)為1幀. 幀長度為1次全量程采集的數(shù)據(jù)長度，它取決于傳感光纖的長度和空間分辨率. 基于此，分布式光纖傳感的數(shù)據(jù)采集模塊架構(gòu)可以統(tǒng)一. 一般來說，光纖Raman分布式傳感的最大傳感距離為10～30 km，而光纖Brillouin分布式傳感的最大傳感距離能達(dá)到50～100 km，空間分辨率小于1 m. 因此在需要采集的數(shù)據(jù)格式中每1幀的長度非常大(104～105AD). 由于FPGA可支持的最大隊(duì)列長度≤131 072子幀，因此這類超長幀需要通過多塊存儲空間拼接才能實(shí)現(xiàn).

1.2采集模塊架構(gòu)

由于各類分布式光纖傳感系統(tǒng)數(shù)據(jù)格式的統(tǒng)一性，采集模塊架構(gòu)可統(tǒng)一設(shè)計(jì)，圖2是抽象后的采集模塊架構(gòu)，適用于ROTDR和BOTDR. 采集模塊的核心為FPGA芯片. 來自傳感光纖的周期性重復(fù)的多幀信號，經(jīng)AD芯片模數(shù)轉(zhuǎn)換后，進(jìn)入FPGA芯片，分布式光纖傳感通常多通道同時(shí)進(jìn)行，圖2是其中1個(gè)通道的工作流程. 進(jìn)入FPGA的周期性多幀傳感信號首先根據(jù)同步信號，被分割成單幀，然后進(jìn)入乒乓操作流程：開辟2塊相同大小的RAM，每塊RAM的大小為1幀傳感信號所需的空間大小，其中一塊RAM用于傳感數(shù)據(jù)的接收和幀與幀之間的累加，另一塊用于多幀傳感信號平均降噪算法處理以及和計(jì)算機(jī)之間的數(shù)據(jù)傳輸. 每當(dāng)數(shù)據(jù)接收完成，交換2塊RAM之間的功能. 平均降噪算法算法耗時(shí)較長，采用這種內(nèi)存的乒乓操作，將數(shù)據(jù)采集、平均降噪算法以及數(shù)據(jù)傳輸過程分開，能夠?qū)崿F(xiàn)2次全量程傳感之間零間隔.

每1幀數(shù)據(jù)按照環(huán)形隊(duì)列的數(shù)據(jù)格式存儲在FPGA的RAM中. 每采集1個(gè)新數(shù)據(jù)，隊(duì)列指針向后移1位，將當(dāng)前指針指向的數(shù)據(jù)取出并與采集的數(shù)據(jù)相加，累加值將指針位置的數(shù)據(jù)替換并壓入隊(duì)列. 上述的隊(duì)列更新流程如圖3A顯示，事實(shí)上，上述流程是傳感數(shù)據(jù)的累加過程，通過這個(gè)流程，可以將每2組數(shù)據(jù)采集的時(shí)間間隔壓縮為0. 在FPGA內(nèi)部，環(huán)形隊(duì)列由FIFO(先入先出)線性隊(duì)列組成，圖3B是首尾相連的、線性的環(huán)形隊(duì)列. Altera公司的FPGA芯片最大支持的FIFO隊(duì)列深度為131 072. 對于分布式光纖檢測來說，如果要實(shí)現(xiàn)很高的空間分辨率(如0.1 m)，即每隔0.1 m采集1個(gè)數(shù)據(jù)，對于100 km長度的檢測則需要106個(gè)數(shù)據(jù)，單個(gè)FIFO隊(duì)列已無法滿足分布式光纖采集所需要的數(shù)據(jù)量，需要將多個(gè)隊(duì)列拼接才能達(dá)到需要的深度.

圖3C顯示了如何利用多個(gè)FIFO隊(duì)列組成1個(gè)長環(huán)形隊(duì)列. 每個(gè)FIFO隊(duì)列至少都有1個(gè)狀態(tài)位和2個(gè)控制位. 狀態(tài)位FULL表示隊(duì)列是否已滿，2個(gè)控制位WRITE_EN和READ_EN用于控制隊(duì)列數(shù)據(jù)的壓入和彈出. 多個(gè)隊(duì)列組成環(huán)形隊(duì)列時(shí)，前一個(gè)隊(duì)列FULL狀態(tài)位連接到后一個(gè)隊(duì)列寫入控制位，同時(shí)控制本隊(duì)列的讀出控制位. 前一隊(duì)列的數(shù)

圖2　分布式光纖傳感技術(shù)中實(shí)現(xiàn)快速檢測的數(shù)據(jù)采集模塊架構(gòu)

圖3　快速檢測的數(shù)據(jù)采集模塊的數(shù)據(jù)結(jié)構(gòu)和算法

據(jù)輸出連接后一隊(duì)列的數(shù)據(jù)輸入. 這樣當(dāng)前一隊(duì)列滿時(shí)，自動(dòng)向下一隊(duì)列擴(kuò)展，實(shí)現(xiàn)多個(gè)隊(duì)列的拼接. 最后一個(gè)隊(duì)列的數(shù)據(jù)接入加法器，和來自AD轉(zhuǎn)換的傳感信號相加，再壓入第一個(gè)隊(duì)列，實(shí)現(xiàn)對多次重復(fù)測量的傳感信號的累加，用于平均降噪算法.

1.3FPGA的芯片選取和設(shè)計(jì)

采集模塊中FPGA芯片的選取主要考慮2個(gè)因素，即端口速率和內(nèi)部RAM大小. 由于光速非?？欤疚倪x取端口速率最高為250 MHz的FPGA芯片. 在數(shù)據(jù)采集模塊中，采用14位高精度模數(shù)轉(zhuǎn)換器，在降噪算法中采用10 000次以上的平均降噪，提高20 dB以上的信噪比[9]，為了提高運(yùn)算速率，平均算法采用移位來實(shí)現(xiàn)，因此平均次數(shù)必須是2的整數(shù)倍，設(shè)計(jì)為214.另外Altera公司FPGA支持最大的FIFO隊(duì)列長度為131 072，因此FIFO開辟的空間最大為131 072×28 bits. 典型的分布式光纖拉曼傳感和分布式光纖布里淵傳感的最大檢測距離相差5～10倍，因此其片內(nèi)RAM大小也不同，綜合成本和性能，對于分布式光纖拉曼傳感，選擇Arria II系列的EP2AGX65DF25C4N， 它擁有5 246 kbits的片內(nèi)存儲器. 而對于分布式光纖布里淵傳感，選擇了EP4SGX530HH35C4，它擁有27 376 kbits的片內(nèi)存儲器. 每個(gè)數(shù)據(jù)采集模塊，都采用4通道采集，實(shí)現(xiàn)分布式傳感的多通道復(fù)用. 上述2款芯片的存取空間能夠滿足4通道同時(shí)進(jìn)行乒乓操作和數(shù)據(jù)傳輸.

1.4USB數(shù)據(jù)傳輸模塊設(shè)計(jì)

USB數(shù)據(jù)傳輸過程是FPGA內(nèi)部乒乓操作的一部分. 內(nèi)部存儲空間RAM分為2部分(圖2)，接收數(shù)據(jù)時(shí)將其中一部分存儲空間RAM1交給接收模塊控制，組成環(huán)形隊(duì)列，進(jìn)行數(shù)據(jù)的接收與累加. 接收數(shù)據(jù)后，RAM1切換給USB模塊控制. 同時(shí)USB模塊放開對另一塊存儲空間RAM2的控制，將其釋放給數(shù)據(jù)接收模塊. 設(shè)計(jì)上位機(jī)(計(jì)算機(jī))的信號傳輸采用USB2.0協(xié)議. 對數(shù)據(jù)接收過程耗時(shí)Tr和USB信號傳輸過程耗時(shí)Tusb做如下估算：

(1)

(2)

式中，neff和L分別為傳感光纖的有效折射率和長度，k為降噪平均操作的次數(shù)，c光速；R為傳感中的空間分辨率，B為采集數(shù)據(jù)的精度，Susb為USB傳輸速率.

這里僅考慮USB凈荷的傳輸時(shí)間，忽略USB報(bào)文的頭部長度和協(xié)議建立時(shí)間，同時(shí)忽略平均降噪算法的耗時(shí)，相對于式(1)、式(2)，所忽略項(xiàng)影響非常小. 設(shè)計(jì)中，B=14 bits，采用14位的采樣精度；R=1，采用1 m的空間分辨率；Susb=48 Mbps，采用USB2.0的全速協(xié)議；neff為1.45. 據(jù)估算，只要重復(fù)采集的次數(shù)k>30，可滿足. 在實(shí)際應(yīng)用中，重復(fù)采集的次數(shù)從103至105量級變化會(huì)有幾千到一萬次，遠(yuǎn)大于上述條件. 因此采用乒乓操作是非常合理的一種設(shè)計(jì)，通過這個(gè)設(shè)計(jì)，可以做到零延遲的數(shù)據(jù)接收.

采用的USB芯片選擇賽普拉斯公司的CY7C68013A，這個(gè)芯片可以配置為普通端口模式(PORTS)、通用可編程接口模式(GPIF)和從器件FIFO(SLAVE)模式3種數(shù)據(jù)傳輸模式. FIFO(SLAVE)模式下，數(shù)據(jù)不需要經(jīng)過CY7C68013的微處理器控制，可以直接進(jìn)行數(shù)據(jù)傳輸，傳輸速度最快理論值能達(dá)480 Mbps. 設(shè)計(jì)中選取FIFO(SLAVE)模式，采用48 Mbps的傳輸速度. USB軟件設(shè)計(jì)分為2部分，一是CY7C68013的固件程序，固件程序通過執(zhí)行8051代碼完成相應(yīng)的功能，實(shí)現(xiàn)USB設(shè)備工作模式的選擇、中斷處理和數(shù)據(jù)的發(fā)送與接收的控制；二是FPGA程序設(shè)計(jì)，進(jìn)行邏輯和時(shí)序的控制. 這2部分共同協(xié)作完成整個(gè)系統(tǒng)數(shù)據(jù)的發(fā)送和接收.

2　系統(tǒng)仿真和調(diào)試

利用Quartus II軟件對FPGA芯片進(jìn)行設(shè)計(jì)、編譯、綜合. 同時(shí)利用Modelsim對各模塊的相關(guān)時(shí)序進(jìn)行仿真，縮短設(shè)計(jì)周期，同時(shí)確保達(dá)到設(shè)計(jì)要求. 同時(shí)也利用邏輯分析儀Sigbal Tap II進(jìn)行在線調(diào)試. 最終完成數(shù)據(jù)采集的相關(guān)功能設(shè)計(jì). 圖4是使用2 km的傳感光纖，利用數(shù)據(jù)采集模塊，在分布式光纖拉曼傳感系統(tǒng)中采集的信號. 數(shù)據(jù)通過USB2.0協(xié)議傳輸?shù)诫娔X，用Labview程序顯示最終采集到的信號波形.

采用2通道同時(shí)采集，圖4是采集的拉曼散射anti-Stokes信號和Stokes信號. 整個(gè)數(shù)據(jù)采集過程采用214次的降噪平均處理，所耗時(shí)間小于1 s，如果使用傳統(tǒng)的數(shù)據(jù)采集算法和架構(gòu)所耗時(shí)間會(huì)在30 s到幾分鐘. 本文的快速檢測軟硬件設(shè)計(jì)，能夠使這類分布式光纖傳感獲得準(zhǔn)實(shí)時(shí)的特性，實(shí)時(shí)性是光纖檢測中一個(gè)重要的特性，小于1 s的實(shí)時(shí)性，能夠使此模塊在火災(zāi)預(yù)警等實(shí)際應(yīng)用場合中的優(yōu)勢顯著.

3　結(jié)論

本文實(shí)現(xiàn)了一種能實(shí)現(xiàn)快速檢測的數(shù)據(jù)采集軟硬件系統(tǒng)，并應(yīng)用在分布式光纖傳感中，將全量程檢測檢測時(shí)間從30 s甚至幾分鐘縮短到10 s以內(nèi). 解決了目前分布式光纖傳感中測量時(shí)間長，不能實(shí)時(shí)檢測的問題，實(shí)現(xiàn)了準(zhǔn)實(shí)時(shí)檢測，大大提高了分布式光纖傳感的性能，有望拓寬其應(yīng)用領(lǐng)域.

圖4　2 km傳感光纖上采集的分布式拉曼傳感信號

[1]Candiani A, Sozzi M, Cucinotta A, et al. Optical fiber ring cavity sensor for label-free DNA detection[J]. IEEE Journal of Selected Topics in Quantum Electronics, 2012, 18(3): 1176-1183.

[2]Candiani A, Giannetti S, Sozzi M, et al. Microstructured optical fiber Bragg grating sensor for DNA detection[J]. Optical Fibers & Sensors for Medical Diagnostics & Treatment Applications XIII, 2013, 8576(3):343-348.

[3]Mokhtar M R, Owens K, Kwasny J, et al. Fiber-optic strain sensor system with temperature compensation for arch bridge condition monitoring[J]. Sensors Journal IEEE, 2012, 12(5): 1470-1476.

[4]Scott R H, Banerji P, Chikermane S, et al. Commissioning and evaluation of a fiber-optic sensor system for bridge monitoring[J]. Sensors Journal IEEE, 2013, 13(7): 2555-2562.

[5]Li L, Han Q, Liu T, et al. Reflective all-fiber current sensor based on magnetic fluids[J]. Review of Scientific Instruments, 2014, 85(8): 083107.

[6]Minardo A, Coscetta A, Porcaro G, et al. Structural health monitoring in the railway field by fiber-optic sensors[M]. Switzerland: Springer International Publishing, 2015: 359-363.

[7]Murayama H, Wada D, Igawa H. Structural health monitoring by using fiber-optic distributed strain sensors with high spatial resolution[J]. Photonic Sensors, 2013, 3:355-376.

[8]Kobyakov A, Sauer M, Chowdhury D. Stimulated Brillouin scattering in optical fibers[J]. Advances in Optics and Photonics, 2010, 2(1): 1-59.

[9]Nikles M, Thevenaz L, Robert P A. Brillouin gain spectrum characterization in single-mode optical fibers[J]. Journal of Lightwave Technology, 1997, 15(10): 1842-1851.

[10]Farahani M A, Gogolla T. Spontaneous Raman scattering in optical fibers with modulated probe light for distributed temperature Raman remote sensing[J]. Journal of Lightwave Technology, 1999, 17(8): 1379-1391.

[11]Dakin J P, Pratt D J, Bibby G W, et al. Temperature distribution measurement using Raman ratio thermometry[C]∥29th Annual Technical Symposium, International Society for Optics and Photonics. San Diego, USA, 1986: 249-256.

[12]Yuan C X, Wang Z N, Jia X H, et al. 94 km Brillouin distributed optical fiber sensors based on ultra-long fiber ring laser pumping[J]. Spectroscopy and Spectral Analysis, 2014, 34(5): 1191-1195.

[13]Barnoski M K, Jensen S M. Fiber waveguides: A novel technique for investigating attenuation characteristics[J]. Applied Optics, 1976, 15(9): 2112-2115.

[14]Dakin J P, Pratt D J, Bibby G W, et al. Distributed optical fibre Raman temperature sensor using a semiconductor light source and detector[J]. Electronics Letters, 1985, 21(13): 569-570.

[15]De Souza K, Lees G P, Wait P C, et al. Diode-pumped Landau-Placzek based distributed temperature sensor utilising an all-fibre Mach-Zehnder interferometer[J]. Electronics Letters, 1996, 32(23): 2174-2175.

[16]張悅, 張記龍, 李曉, 等. 光纖中拉曼散射回波超高靈敏度探測技術(shù)研究[J]. 光譜學(xué)與光譜分析, 2009, 29(5): 1300-1303.

Zhang Y, Zhang J L, Li X, et al. Research on the supersensitive detecting technology of backward Raman scattering signal in optic fiber[J]. Spectroscopy and Spectral Analysis, 2009, 29(5): 1300-1303.

【中文責(zé)編：譚春林英文責(zé)編：肖菁】

Hardware and Software Design for Rapid Detecting Property in the Distributed Fiber Optic Sensing System

Zhou Bin*， Huang Guojing, Lu Changtao, Ruan Hang

(South China Academy of Advanced Optoelectronics, South China Normal University, Guangzhou 510006, China)

A data acquisition system of the distributed sensing technologies is developed with rapid detection property using FPGA, for slow detection speed of a full-scale detecting of current distributed optical fiber sensing system. With our fast data acquisition system, the time consumption of a full-scale detection is suppressed to be less than ten seconds from 30s or even a few minutes, making it possible for the quasi real-time monitoring. Distributed fiber-optic sensing technology is one of the most promising fiber-optic sensor technologies. Its biggest advantage is the large range and long distance detection. Both of Brillouin and Raman distributed sensing principle and its repeatability, frame structure characteristics of the sensor data are introduced. According to these characteristics, long ring queue of a number of FIFO in FPGA are designed to buffer and average data, for signal noise reduction. Data transfer module based on USB2.0 of CY7C68013A is designed to transfer data to the PC for displaying and storage. Results showed that the design realizes quasi real-time detection, solving the lack of long time and not real-time detection. It greatly increases the performance of distributed optical fiber sensing. The proposed sensing system can be used in monitoring the power line systems, railway, subway tracks, tunnels, dams or landslide, where require real-time monitoring and a large number of sensors.

distributed fiber sensor; data acquisition; optical fiber Brillouin effect; optical fiber Raman effect

2015-04-12《華南師范大學(xué)學(xué)報(bào)(自然科學(xué)版)》網(wǎng)址：http://journal.scnu.edu.cn/n

國家自然科學(xué)基金項(xiàng)目(6130705)；廣東省引進(jìn)創(chuàng)新團(tuán)隊(duì)計(jì)劃資助項(xiàng)目(201001D0104799318)；中國博士后基金項(xiàng)目(2013M 531866)

周斌，講師，Email:zhoubin_mail@163.com.

O438.1

A

1000-5463(2015)05-0018-05