左程吉,王祝文,向　旻,周大鵬,劉宗利,楊　帆

(1.吉林大學(xué)地球探測(cè)科學(xué)與技術(shù)學(xué)院,吉林長(zhǎng)春130026;2.中國(guó)石油天然氣股份有限公司遼河油田分公司勘探開發(fā)研究院,遼寧盤錦124010)

基于電成像測(cè)井孔隙度分析技術(shù)的火山巖孔隙徑向非均質(zhì)性研究

左程吉1,王祝文1,向旻1,周大鵬1,劉宗利1,楊帆2

(1.吉林大學(xué)地球探測(cè)科學(xué)與技術(shù)學(xué)院,吉林長(zhǎng)春130026;2.中國(guó)石油天然氣股份有限公司遼河油田分公司勘探開發(fā)研究院,遼寧盤錦124010)

摘要:火山巖儲(chǔ)集層孔隙結(jié)構(gòu)復(fù)雜,孔隙發(fā)育類型多,產(chǎn)能預(yù)測(cè)困難。為此,根據(jù)電成像測(cè)井的方法原理,結(jié)合Archie公式和常規(guī)測(cè)井資料,將火山巖層的電成像數(shù)據(jù)轉(zhuǎn)換為核磁區(qū)間孔隙度類似的形式,通過與核磁區(qū)間孔隙度對(duì)比,得到火山巖儲(chǔ)層的孔隙區(qū)間特征。在此基礎(chǔ)上,針對(duì)火山巖復(fù)雜的孔隙類型,對(duì)孔隙區(qū)間展開進(jìn)一步分析,得到不同孔隙區(qū)間所占的比例,結(jié)合計(jì)算的孔隙度變異系數(shù),得到火山巖儲(chǔ)集層的孔隙區(qū)間特征及孔隙徑向非均質(zhì)性特征,進(jìn)而發(fā)現(xiàn)影響火山巖儲(chǔ)層產(chǎn)能的隔擋層的存在,為火山巖次生孔隙發(fā)育情況及火山巖儲(chǔ)集層產(chǎn)能預(yù)測(cè)提供了參考依據(jù)。

關(guān)鍵詞：電成像;孔隙度譜分析;火山巖;次生孔隙;產(chǎn)能預(yù)測(cè)

電成像測(cè)井覆蓋面積大,成果資料不僅直觀形象,而且具有較高的分辨率。近年來，電成像測(cè)井在地層精細(xì)結(jié)構(gòu)評(píng)價(jià)和不同巖性儲(chǔ)層產(chǎn)能預(yù)測(cè)等方面的應(yīng)用有了很大的進(jìn)展[1-2],將成像測(cè)井資料刻度為孔隙度就是其中一個(gè)方面[3-4]。目前國(guó)外的一些測(cè)井軟件已經(jīng)可以將電成像測(cè)井資料轉(zhuǎn)換為孔隙度頻譜進(jìn)行分析,但國(guó)內(nèi)的大部分軟件與其相比還存在一定的問題和差距。

火山巖層比其它巖層孔隙結(jié)構(gòu)更復(fù)雜,孔隙類型有不同成因的氣孔、收縮孔、破碎縫、溶蝕孔、晶間孔、構(gòu)造縫等,并且孔隙之間的連通性變化較大,火山巖孔隙的徑向非均質(zhì)性研究對(duì)于火山巖儲(chǔ)層產(chǎn)能預(yù)測(cè)具有重要意義。景建恩等[5]應(yīng)用電成像測(cè)井?dāng)?shù)據(jù)對(duì)裂縫型碳酸鹽巖儲(chǔ)層進(jìn)行了測(cè)井評(píng)價(jià),吳興能等[6]研究了將電成像測(cè)井資料變換為孔隙度分布圖像的方法,周彥球等[7]通過成像測(cè)井孔隙度頻譜與巖心孔隙分析資料的對(duì)比研究,證實(shí)了成像測(cè)井用于微觀解釋的可能性。我們利用Archie公式將電成像資料刻度為孔隙度,進(jìn)而得到每個(gè)采樣間隔的孔隙度分布譜,再對(duì)分布譜進(jìn)行譜分析,得到與核磁區(qū)間孔隙度相似的孔隙度分布區(qū)間。通過對(duì)比核磁區(qū)間孔隙度與電成像孔隙度分布區(qū)間,分析火山巖儲(chǔ)層的孔隙特征,再利用分布譜得到不同孔隙度所占的比例,結(jié)合孔隙度變異指數(shù),展示出孔隙度的徑向非均質(zhì)性,為火山巖的產(chǎn)能預(yù)測(cè)提供參考。

1孔隙度分布譜的計(jì)算

1.1電阻率刻度

電成像的紐扣電極得到的是地層電導(dǎo)率,需要將它表示為電阻率。已有研究成果表明,電成像測(cè)井電流流動(dòng)方式、聚集形式及探測(cè)深度范圍與常規(guī)測(cè)井中的淺側(cè)向測(cè)井相當(dāng)[8-9],可以選擇淺側(cè)向電阻率對(duì)電成像電導(dǎo)率進(jìn)行刻度[10-11]:

(1)

1.2孔隙度標(biāo)定

電成像測(cè)井及常規(guī)淺側(cè)向測(cè)井的探測(cè)深度大致為沖洗帶,故應(yīng)用經(jīng)典的Archie公式進(jìn)行孔隙度標(biāo)定:

(2)

式中:a,b為與巖性有關(guān)的系數(shù);m為膠結(jié)指數(shù);n為飽和度指數(shù);φ為孔隙度;Sxo為沖洗帶的含水飽和度;Rxo為沖洗帶電阻率;Rmf為泥漿溶液電阻率。

將利用刻度后的電成像電阻率Ri得到的孔隙度記作為φi,用常規(guī)淺側(cè)向電阻率RLLS得到的孔隙度記作φ0,代入公式(2)并整理得:

(3)

將公式(1)代入公式(3),整理得到每個(gè)紐扣電極經(jīng)過標(biāo)定后的孔隙度值:

(4)

根據(jù)公式(4)得到的孔隙度值消除了多個(gè)不確定性系數(shù)的影響,僅剩膠結(jié)指數(shù)m。m值的大小與多種因素有關(guān),研究層段的地層水電阻率、礦化度、溫度、壓力等基本一致,所以對(duì)m值的選取保持一致。由于火山巖的膠結(jié)物及膠結(jié)方式較復(fù)雜,所以選取的m值較大。在利用公式(4)得到每個(gè)紐扣電極孔隙度值的基礎(chǔ)上,用固定大小的窗長(zhǎng)及一定長(zhǎng)度的步長(zhǎng)對(duì)所有孔隙度值進(jìn)行統(tǒng)計(jì)分析[12-14],可以得到每個(gè)窗口內(nèi)的孔隙度頻譜(圖1)。圖1中,φmin為統(tǒng)計(jì)窗口內(nèi)的最小孔隙度;φph為主頻孔隙度,指示頻率最高點(diǎn)的孔隙度,一般與窗長(zhǎng)內(nèi)平均孔隙度接近;φw為主峰右側(cè)寬度,指示次生孔隙發(fā)育程度;φmax為統(tǒng)計(jì)窗口內(nèi)最大孔隙度[15]。圖2中第4道為遼河盆地東部凹陷Y井4450～4515m火山巖孔隙度分布譜。當(dāng)?shù)貙哟紊紫恫话l(fā)育、原生孔隙發(fā)育較均勻時(shí),孔隙度分布譜表現(xiàn)為前移的單峰;當(dāng)原生孔隙較小而次生孔隙發(fā)育較均勻時(shí),孔隙度分布譜表現(xiàn)為后移的單峰;當(dāng)原生孔隙發(fā)育且裂縫、溶蝕孔洞等次生孔隙發(fā)育不均勻時(shí),孔隙度分布譜表現(xiàn)為雙峰或多峰特征。因此,利用孔隙度分布譜可以在一定程度上對(duì)地層的孔隙分布情況進(jìn)行評(píng)價(jià)。為了更加精確地評(píng)價(jià)火山巖復(fù)雜的孔隙結(jié)構(gòu)及孔隙非均質(zhì)性,以下對(duì)孔隙度分布譜做進(jìn)一步分析。

圖1　一定窗長(zhǎng)內(nèi)的孔隙度分布譜及參數(shù)

圖2　遼河盆地東部凹陷Y井4450～4515m火山巖孔隙度分布譜

2孔隙度分析技術(shù)及應(yīng)用

2.1核磁區(qū)間孔隙度計(jì)算方法

核磁共振測(cè)井探測(cè)的是氫核的共振大小及其衰減,信號(hào)大小反映了地層的孔隙度,衰減速率反映了巖石的弛豫性質(zhì)。不同孔隙大小的流體中的氫核有不同的弛豫速率,通過核磁共振測(cè)井處理可以將不同弛豫速率與所對(duì)應(yīng)的信號(hào)大小反映出來,得到橫向弛豫時(shí)間T2的分布譜(圖3)。對(duì)T2分布譜從左到右、幅值從小到大按4,8,16,32,64,128,256,512,1024,2048ms等累加、排序,分別得到不同大小區(qū)間孔隙度；將各區(qū)間的孔隙度用不同顏色及符號(hào)逐項(xiàng)顯示,得到各個(gè)區(qū)間孔隙度(不同顏色區(qū)間)隨深度的變化趨勢(shì),即核磁區(qū)間孔隙度(圖2中第5道)。各區(qū)間孔隙度的累加和(區(qū)間孔隙度的包絡(luò))為有效孔隙度。

圖3　遼河盆地東部凹陷T2分布譜

2.2孔隙度分布譜分析技術(shù)

依據(jù)核磁區(qū)間孔隙度的計(jì)算方法,對(duì)電成像測(cè)井得到的孔隙度分布譜進(jìn)行類似處理。首先通過累加的方式統(tǒng)計(jì)出每個(gè)分布譜各個(gè)孔隙度區(qū)間的孔隙度曲線(分別代表各個(gè)采樣窗口內(nèi)孔隙度大小為0.03,0.05,0.07,0.10,0.13,0.15,0.17,0.20,0.50以下的數(shù)據(jù)構(gòu)成的區(qū)間組分),再用不同的顏色逐項(xiàng)顯示,得到電成像孔隙度分布區(qū)間圖(圖2中第6道)。對(duì)比分析核磁區(qū)間孔隙度與電成像孔隙度,可清晰直觀地看出不同孔隙度大小對(duì)總孔隙度的貢獻(xiàn)情況。

火山巖儲(chǔ)層孔隙結(jié)構(gòu)復(fù)雜,多由裂縫、孔洞發(fā)育構(gòu)成孔隙連通的通道,應(yīng)用上述方法可以看到不同大小的孔隙發(fā)育情況以及不同巖性的孔隙發(fā)育情況的差別。遼河盆地東部凹陷Y井4450～4500m為油斑粗面質(zhì)角礫巖及油跡粗面質(zhì)角礫巖,4500～4510m為熒光粗面巖,其中4450～4495m為火山巖儲(chǔ)層,產(chǎn)水72.1m3/d,產(chǎn)油17.9m3/d。從孔隙度分布譜上發(fā)現(xiàn),角礫巖儲(chǔ)集層段多呈雙峰或多峰分布,且孔隙度較大,說明大孔隙發(fā)育。從圖2核磁區(qū)間孔隙度圖與電成像區(qū)間孔隙度圖可以看出,該區(qū)角礫巖儲(chǔ)集層段的整體孔隙度較大,在15%左右,而且大孔隙所占比例較大,孔隙大小發(fā)育比較均勻,連通性也較好,所以儲(chǔ)集性能較好;粗面巖地層在孔隙度分布譜上多為單峰分布,且孔隙度較小,在區(qū)間孔隙度圖上顯示整體孔隙度較小,只有5%左右,而且都以小孔為主,所以粗面巖儲(chǔ)層物性和含油性較差。

以上分析發(fā)現(xiàn),可能成為火山巖儲(chǔ)集層的孔隙條件為:總孔隙度較大,孔隙度分布譜表現(xiàn)為后移的雙峰或多峰,在一定深度范圍內(nèi)大孔隙占較大比例且大孔隙發(fā)育程度較均勻。

值得注意的是,研究區(qū)4496～4505m層段為非儲(chǔ)集層,但孔隙度分布譜上多呈雙峰或多峰分布,整體孔隙度也較大且大孔隙有所發(fā)育。為解釋這一現(xiàn)象的原因,以便更加準(zhǔn)確地預(yù)測(cè)火山巖儲(chǔ)層的產(chǎn)能,下面對(duì)火山巖孔隙徑向非均質(zhì)性進(jìn)行研究。

3火山巖孔隙非均質(zhì)性研究

為了確定不同孔隙度所占的具體比例,我們對(duì)孔隙度分布譜進(jìn)行了進(jìn)一步分析。首先統(tǒng)計(jì)每個(gè)窗口孔隙度分布譜的最小孔隙度PorMin與最大孔隙度PorMax以及主頻孔隙度PorPeak,然后對(duì)每個(gè)窗口分布譜從左到右按孔隙度值從小到大的順序統(tǒng)計(jì)每個(gè)窗口內(nèi)累積頻率為0.2,0.4,0.6,0.8的臨界點(diǎn)的孔隙度值,分別得到孔隙度分布曲線Por0.2,Por0.4,Por0.6與Por0.8,在圖像上用不同顏色依次填充,反映不同孔隙度所占的比例,即孔隙徑向非均質(zhì)性(圖4中第4道)。

對(duì)每個(gè)窗口孔隙度分布譜進(jìn)行統(tǒng)計(jì)分析,可以得到頻率的最大值所對(duì)應(yīng)的孔隙度值及孔隙度的最大值,即圖1中的主頻孔隙度φph與最大孔隙度φmax,進(jìn)而分別得到孔隙度曲線PorPeak與PorMax。填充兩條曲線的中間區(qū)域,可以得到主頻右側(cè)寬度(圖4中第5道),該區(qū)域反映次生孔隙的發(fā)育尺度,不同大小孔隙越發(fā)育,該區(qū)域?qū)挾仍酱?進(jìn)而可以指示次生孔隙的非均質(zhì)性情況。在此基礎(chǔ)上,計(jì)算了孔隙度變異系數(shù)(圖4中第6道),它類似于油藏描述中的滲透率參數(shù),是每個(gè)采樣窗口內(nèi)所有數(shù)據(jù)點(diǎn)的孔隙度的均方差與孔隙度的均值的比值,指示了孔隙的非均質(zhì)性?？紫抖茸儺愊禂?shù)值越大,說明孔隙的徑向非均質(zhì)性越強(qiáng)。

分析圖4中的孔隙度分布曲線,可以發(fā)現(xiàn)遼河盆地東部凹陷Y井4450～4495m深度段的孔隙度值都較大,大孔隙占非常大的比例,并且主峰右側(cè)寬度較窄,孔隙度變異系數(shù)較小,說明次生孔隙發(fā)育程度均勻,反映孔隙大小均勻且連通性較好,所以此段儲(chǔ)集性能較好。進(jìn)一步觀察此層段發(fā)現(xiàn),在4475,4485m等深度點(diǎn)附近孔隙度變異系數(shù)較大,次生大孔隙發(fā)育尺度較多且小孔隙發(fā)育程度不均勻,可以推斷這些深度點(diǎn)附近存在隔擋層。

在4496～4505m深度段,對(duì)孔隙的非均質(zhì)性研究發(fā)現(xiàn),4496m處孔隙度較小,不同大小的孔隙所占的比例都很小,所以此處可能有隔擋層隔斷了儲(chǔ)集層的連通性。另外,此段的孔隙發(fā)育尺度較大,次生孔隙發(fā)育不均勻,且孔隙度變異系數(shù)普遍較大,說明此段孔隙的均勻性及連通性很差,即此段的儲(chǔ)集性能較差。

圖4　遼河盆地東部凹陷Y井火山巖孔隙非均質(zhì)性特征

以上分析可知,對(duì)火山巖孔隙的徑向非均質(zhì)性進(jìn)行研究,可以更加精確地評(píng)價(jià)孔隙發(fā)育特征,為評(píng)價(jià)火山巖儲(chǔ)層產(chǎn)能提供參考依據(jù)。

4結(jié)論

1) 利用Archie公式和常規(guī)測(cè)井資料,可將電成像測(cè)井資料標(biāo)定為孔隙度,進(jìn)而得到孔隙度分布譜,通過對(duì)譜峰的分布和參數(shù)進(jìn)行分析,可以對(duì)地層孔隙發(fā)育情況進(jìn)行有效解釋。

2) 對(duì)比分析電成像分布譜孔隙度與核磁區(qū)間孔隙度,可以得到火山巖儲(chǔ)層復(fù)雜的孔隙結(jié)構(gòu)特征:總孔隙度較大,孔隙度分布譜表現(xiàn)為后移的雙峰或多峰,在一定深度范圍內(nèi)大孔隙占較大比例,大孔隙發(fā)育程度較均勻。

3) 通過火山巖孔隙徑向非均質(zhì)性研究,可對(duì)一些孔隙度較大但儲(chǔ)集性能較差的現(xiàn)象做出解釋，判斷儲(chǔ)層中隔擋層的位置,進(jìn)而對(duì)火山巖儲(chǔ)層作出更精細(xì)的評(píng)價(jià)。

參考文獻(xiàn)

[1]王偉華,劉呈冰,于恒,等.微電阻率成像測(cè)井在孔洞型儲(chǔ)層定量評(píng)價(jià)中的應(yīng)用[J].國(guó)外測(cè)井技術(shù),1999,14(6):50-55

WANG W H,LIU C B,YU H,et al.The application of micro-resistivity imaging logging in quantitative evaluation of vug reservoir[J].Foreign Well Logging Technology,1999,14(6):50-55

[2]HORNBY B E,JOHNSTON D L.Fracture evaluation from the borehole stoneley wave[J].Geophysics,1989,54(10):1274-1288

[3]張守謙,顧純學(xué),曹廣華.成像測(cè)井技術(shù)及應(yīng)用[M].北京:石油工業(yè)出版社,1997:50-56

ZHANG S Q,GU C X,CAO G H.Imaging logging technology and its application[M].Beijing:Petroleum Industry Press,1997:50-56

[4]石油測(cè)井情報(bào)協(xié)作組.測(cè)井新技術(shù)應(yīng)用[M].北京:石油工業(yè)出版社,1998:5-10

PETROLEUM WIRELINE LOGGING INFORMATION COORDINATING GROUP.New technology in wireline logging[M].Beijing:Petroleum Industry Press,1998:5-10

[5]景建恩,魏文博,梅忠武,等.裂縫型碳酸鹽巖儲(chǔ)層測(cè)井評(píng)價(jià)方法:以塔河油田為例[J].地球物理學(xué)進(jìn)展,2005,20(1):78-82

JING J E,WEI W B,MEI Z W,et al.Method of well-logging interpretation for fracture reservoirs of carbonate rock:a case study in Tahe oil-field[J].Progress in Geophysics,2005,20(1):78-82

[6]吳興能,劉瑞林,雷軍,等.電成像測(cè)井資料變換為孔隙度分布圖像的研究[J].測(cè)井技術(shù),2008,32(1):53-56

WU X N,LIU R L,LEI J,et al.Study on converting electrical imaging log into porosity distribution image[J].Well Logging Technology,2008,32(1):53-56

[7]周彥球,李曉輝,范曉敏,等.成像測(cè)井孔隙度頻譜分析技術(shù)與巖心孔隙分析資料對(duì)比研究[J].測(cè)井技術(shù),2014,38(3):309-314

ZHOU Y Q,LI X H,FAN X M,et al.A comparative study of porosity spectrum techniques of micro-electric imaging logging and core analysis of porosity[J].Well Logging Technology,2014,38(3):309-314

[8]KUMAR T A,BHADURI A.井眼微電阻率成像測(cè)井在碳酸鹽巖儲(chǔ)層孔隙度分析中的應(yīng)用[J].羅景美,李愛華,趙文杰,譯.測(cè)井與射孔,2003,6(2):23-27

KUMAR T A,BHADURI A.Porosity analysis using borehole electrical images in carbonate reservoirs[J].LUO J M,LI A H,ZHAO W J,translater.Well Logging & Perforation,2003,6(2):23-27

[9]李曉輝,周彥球,緱艷紅,等.電成像測(cè)井孔隙度分析技術(shù)及其在碳酸鹽巖儲(chǔ)層產(chǎn)能預(yù)測(cè)中的應(yīng)用[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2012,42(4):929-930

LI X H,ZHOU Y Q,GOU Y H,et al.Porosity analysis of micro-electric imaging logging and its application in carbonate reservoir production capacity forecast[J].Journal of Jilin University(Earth Science Edition),2012,42(4):929-930

[10]張程恩,潘保芝,張曉峰,等.FMI測(cè)井資料在非均質(zhì)儲(chǔ)層評(píng)價(jià)中的應(yīng)用[J].石油物探,2011,50(6):630-633

ZHANG C E,PAN B Z,ZHANG X F,et al.Application if FMI data in evaluation of heterogeneity reservoir[J].Geophysical Prospecting for Petroleum,2011,50(6):630-633

[11]雍世和,張超謨.測(cè)井?dāng)?shù)據(jù)處理與綜合解釋[M].東營(yíng):石油大學(xué)出版社,2004:514-518

YONG S H,ZHANG C M.Data processing and comprehensive interpretation[M].Dongying:Petroleum University Press,2004:514-518

[12]張冰.利用測(cè)井資料研究沉積相與儲(chǔ)層非均質(zhì)性[D].長(zhǎng)春:吉林大學(xué),2010

ZHANG B.Research on sedimentary facies and reservoir heterogeneity using logging data[D].Changchun:Jilin University,2010

[13]嚴(yán)科,楊少春,任懷強(qiáng).儲(chǔ)層宏觀非均質(zhì)性定量表征研究[J].石油學(xué)報(bào),2008,29(6):870-879

YAN K,YANG S C,REN H Q.Research on quantitative characterization of macroscopic heterogeneity of reservoir[J].Acta Petrolei Sinica,2008,29(6):870-879

[14]LIU R L,LI N,FENG Q F,et al.Application of the tripe porosity model in well-log effectiveness estimation of carbonate in Tarim Oil-field[J].Journal of Petroleum Science and Engineering,2009,68(1/2):40-46

[15]楊億前.裂縫型低滲透儲(chǔ)層剩余油分布規(guī)律研究[D].荊州:長(zhǎng)江大學(xué),2012

YANG Y Q.Remanding oil distribution in fractured reservoir with low permeability[D].Jingzhou:Yangtze University,2012

(編輯：戴春秋)

The radial pore heterogeneity of volcanic reservoir based on the porosity analysis of micro-electric imaging logging

ZUO Chenji1,WANG Zhuwen1,XIANG Min1,ZHOU Dapeng1,LIU Zongli1,YANG Fan2

(1.CollegeofGeo-explorationScienceandTechnology,JilinUniversity,Changchun130026,China; 2.ResearchInstituteofExplorationandDevelopment,PetroChinaLiaoheOilfieldBranchCompany,Panjin124010,China)

Abstract:The volcanic reservoir is characterized by complex pore structures and various pore types; therefore,the productivity forecast is difficult.In terms of the principles of micro-electric imaging logging,combining with the Archie formula and conventional logging data,the micro-electric imaging data of volcanic formations is transformed into the form similar to the nuclear magnetic interval porosity.The porosity interval characteristics of volcanic reservoir are obtained by comparing with nuclear magnetic interval porosity.On this basis,aiming at the complex pore types existing in volcanic reservoir,we carry out further analysis on porosity intervals and achieve the proportion of different porosity intervals.Integrating with the calculated porosity variation coefficient,the characteristics of porosity intervals and radial pore heterogeneity are achieved.As a result,the interlayers in volcanic reservoir are identified.The research results provide reference for the development of secondary pores in volcanics and productivity forecast of volcanic reservoirs.

Keywords:micro-electric imaging,porosity analysis,volcanic,secondary pore,productivity forecast

收稿日期：2015-11-10;改回日期：2016-02-18。

作者簡(jiǎn)介：左程吉( 1991—),男,碩士在讀,主要從事地球物理測(cè)井新方法、新技術(shù)研究工作。

通訊作者：王祝文( 1961—),男,教授,博士生導(dǎo)師,主要從事地球物理測(cè)井新方法、新技術(shù)以及復(fù)雜巖性的地球物理測(cè)井解釋與評(píng)價(jià)、核地球物理、輻射與環(huán)境評(píng)價(jià)等方面的教學(xué)與研究工作。

基金項(xiàng)目：國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目(2009CB219303)及吉林大學(xué)研究生創(chuàng)新基金項(xiàng)目(2015131)共同資助。

中圖分類號(hào)：P631

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-1441(2016)03-0449-06

DOI:10.3969/j.issn.1000-1441.2016.03.016

This research is financially supported by the National Key Basic Research and Development Program of China (973 Program) (Grant No.2009CB219303) and the Graduate Innovation Fund of Jilin University (Grant No.2015131).