劉宗利,王祝文,劉菁華,馮玉輝,曲 東

(1.吉林大學地球探測科學與技術(shù)學院,吉林長春130021;2.吉林大學地球科學學院,吉林長春130061;3.長城鉆探工程有限公司測井公司,遼寧盤錦124011)

遼河盆地東部凹陷中基性火成巖測井曲線特征研究

劉宗利1,王祝文1,劉菁華1,馮玉輝2,曲 東3

(1.吉林大學地球探測科學與技術(shù)學院,吉林長春130021;2.吉林大學地球科學學院,吉林長春130061;3.長城鉆探工程有限公司測井公司,遼寧盤錦124011)

受郯廬斷裂的影響,新生代以來,遼河盆地東部凹陷間斷發(fā)生了6次構(gòu)造運動,且均伴隨有強度不等的火山噴發(fā)活動,火成巖廣泛發(fā)育。以6口不同的井為例,揭示了東部凹陷內(nèi)中基性火成巖、致密玄武巖、氣孔玄武巖、玄武質(zhì)火山碎屑巖、粗面巖、粗面質(zhì)火山碎屑巖和輝綠巖的測井曲線組合特征。選取典型的巖性數(shù)據(jù)段總結(jié)出了該地區(qū)火成巖性測井響應特征,并制作了一套適合該地區(qū)的巖性識別圖版。其中，自然伽馬和深側(cè)向電阻率交會圖能很好地識別輝綠巖,中子和密度交會圖可以很好地區(qū)分玄武巖類;利用最小二乘支持向量機建立了火成巖巖性識別模型,識別結(jié)果與巖心、巖屑薄片定名對比,符合率較高。交會圖版和最小二乘支持向量機結(jié)合運用,很好地解決了研究區(qū)巖性識別的問題,為進一步開展油氣藏的解釋奠定了基礎(chǔ)。

遼河盆地東部凹陷;中基性火成巖;測井曲線;交會圖;巖性識別;最小二乘支持向量機

1 研究區(qū)概況

遼河盆地大地構(gòu)造位于華北板塊東北部[1],構(gòu)造上可分為沈北凹陷、大民屯凹陷、西部凸起、西部凹陷、中央凸起、東部凹陷和東部凸起7個單元。東部凹陷的東西兩側(cè)被東部凸起、中央凸起所夾持,北至大民屯凹陷[2],是一軸向約為北東43°展布的狹長含油氣凹陷(圖1)。

古近紀(65Ma)以來,太平洋板塊西緣俯沖方向由北北西轉(zhuǎn)變?yōu)楸蔽魑?。同時,印度板塊對中國大陸中部和東部產(chǎn)生北向推擠力。在這種情況下,郯廬斷裂從左旋活動轉(zhuǎn)向右旋活動[3-5]。在其影響下，新生代以來間斷發(fā)生了6次構(gòu)造運動且均伴隨有強度不等的火山噴發(fā)活動[6-7],火成巖廣泛發(fā)育,其中以沙河街組時期最具代表性。沙河街組縱向上發(fā)育大套玄武巖和粗面巖[8-9]。在東營組時期,巖漿侵入沙三段和沙一段地層,形成輝綠巖。溢流相發(fā)育致密玄武巖、氣孔玄武巖、粗面巖;火山通道相、爆發(fā)相多發(fā)育火山碎屑巖;侵入相發(fā)育輝綠巖[10-13]。

近年來,在對松遼、準噶爾、遼河盆地的勘探開發(fā)中,發(fā)現(xiàn)中基性火成巖具有良好的油氣藏特征[14-15],因此火成巖巖性的分類識別具有重要的現(xiàn)實意義。測井曲線特征的差異能很好地反應巖性特征。以常規(guī)曲線——自然伽馬(γ)、密度(ρ)、聲波時差(ΔtAC)、補償中子(φCNL)、深側(cè)向電阻率(RLLD)為基礎(chǔ),分析中基性火成巖巖性與測井曲線特征之間的關(guān)系,以期能更好地識別火成巖巖性,并為進一步的火成巖儲層評價提供幫助。

圖1 遼河盆地構(gòu)造分布

2 火成巖礦物組分及測井曲線特征

火成巖識別方法比較多,通常利用自然伽馬(γ)、密度(ρ)、聲波時差(ΔtAC)、補償中子(φCNL)、深側(cè)向電阻率(RLLD)等常規(guī)測井參數(shù)來識別火成巖巖性。圖2給出了火成巖的典型巖心及薄片顯微特征;圖3 到圖8給出了火成巖的測井曲線特征。

2.1 玄武巖

2.1.1 致密玄武巖

以X1井2876.5～2883.0m為例進行說明。致密玄武巖具有斑狀結(jié)構(gòu)、塊狀構(gòu)造。斑晶約占70%,主要由斜長石(75%)和輝石(25%)組成,斜長石較新鮮,輝石多已皂石化(圖2a),基質(zhì)為間粒-間隱結(jié)構(gòu),主要由斜長石、輝石以及玻璃質(zhì)組成,基質(zhì)蝕變中等。測井曲線γ在47～62API,ρ在2.45～2.64g/cm3,ΔtAC在61～72μs/ft(1ft=30.48cm),φCNL在19%～29%,RLLD在8～16Ω·m。表現(xiàn)為低自然伽馬、低電阻率、低聲波時差的特征(圖3)。

2.1.2 氣孔玄武巖

以X2井3682.5～3689.0m為例對氣孔玄武巖進行描述。其具有斑狀結(jié)構(gòu),氣孔-杏仁構(gòu)造,斑晶約占70%,主要由斜長石(80%)和輝石(20%)組成,部分斜長石可見聚片雙晶,輝石斑晶未見蝕變,基質(zhì)為間粒-間隱結(jié)構(gòu),主要由斜長石、輝石以及玻璃質(zhì)組成,基質(zhì)中見斜長石微晶、輝石小顆粒及玻璃質(zhì),蝕變較強。基質(zhì)中暗色礦物均已蝕變,多蝕變?yōu)樵硎?、綠泥石,氣孔內(nèi)充填方解石(圖2b),見氣孔中為方解石、硅質(zhì)兩期充填。此外溶蝕孔也較發(fā)育。

測井曲線γ在36～49API,ρ在2.31～2.54g/cm3,ΔtAC在59～87μs/ft(1ft=30.48cm),φCNL在20%～36%,RLLD在4～13Ω·m,深、淺側(cè)向電阻率范圍一致。與致密玄武巖相比,γ和ρ偏低,φCNL值偏大,深、淺側(cè)向電阻率偏小(圖4)。

2.2 玄武質(zhì)火山碎屑巖

X3井2467～2475m為玄武質(zhì)火山碎屑巖,凝灰結(jié)構(gòu),顆粒直徑均小于2cm,主要由漿屑、玻屑、晶屑、玄武巖巖屑(圖2c)組成,分選中等,磨圓較差、玻屑多呈弧面多角狀,晶屑、巖屑多為次棱角狀,火山灰膠結(jié)。巖石整體裂縫較為發(fā)育,多呈不規(guī)則狀,方解石充填。

測井曲線γ在24～39API,ρ在2.06～2.40g/cm3,ΔtAC在77～87μs/ft(1ft=30.48cm),φCNL在31%～35%,RLLD在3～4Ω·m,深、淺側(cè)向電阻率有分離現(xiàn)象。與玄武巖相比,自然伽馬值和電阻率明顯偏低,密度也偏小,聲波時差和補償中子值較大(圖5)。

2.3 粗面巖

以X4井2840～2926m對粗面巖巖性和測井曲線特征進行說明。巖心整體呈灰綠色,斑狀結(jié)構(gòu),斑晶主要為長石(圖2d),含量約為30%,斑晶有少數(shù)蝕變,有環(huán)帶結(jié)構(gòu),多呈長柱狀,表面裂紋發(fā)育,致密塊狀?；|(zhì)由堿性長石和少量斜長石組成,為間粒間隱結(jié)構(gòu)。

測井曲線γ在91～180API,ρ在2.11～2.59g/cm3,ΔtAC在52～93μs/ft(1ft=30.48cm),φCNL在4%～29%,RLLD在5～606Ω·m(圖6)。相比玄武巖自然伽馬和電阻率明顯呈現(xiàn)高值,粗面巖電阻率幅值變化大。可以用γ和雙側(cè)向電阻率曲線特征識別粗面巖和玄武巖巖性。

2.4 粗面質(zhì)火山碎屑巖

X5井3340.5～3398.5m為粗面質(zhì)火山碎屑巖,火山角礫結(jié)構(gòu),角礫直徑多大于2cm,主要由粗面質(zhì)巖屑(圖2e)、晶屑組成,分選、磨圓均較差,角礫多為棱角-次棱角狀,火山灰膠結(jié)。

測井曲線γ在69～115API,ρ在2.37～2.64g/cm3,ΔtAC在65～76μs/ft(1ft=30.48cm),φCNL在12%～23%,RLLD在16～61Ω·m(圖7)。相比粗面巖自然伽馬值偏低,密度值稍大,雙側(cè)向電阻率呈現(xiàn)低值。

2.5 輝綠巖

以X6井3638～3658m對輝綠巖進行巖性和測井曲線特征說明。輝綠結(jié)構(gòu),主要礦物為斜長石和輝石,結(jié)晶較粗,輝石多已蝕變綠泥石化,見少量黑色不透明磁鐵礦。發(fā)育較多淋濾孔,暗色礦物輝石蝕變?yōu)樵硎?圖2f)、綠泥石等粘土礦物后,經(jīng)過后期作用,粘土礦物被淋濾形成空隙,邊部仍可見綠泥石等,少量孔隙為斜長石溶蝕孔隙。

測井曲線γ在46～74API,ρ在2.55～2.71g/cm3,ΔtAC在58～75μs/ft(1ft=30.48cm),φCNL在16%～25%,RLLD在93～317Ω·m(圖8)。相對于玄武巖,電阻率值偏大,其余特征不太明顯。相對于粗面巖,γ明顯偏低,電阻率值偏小。

綜合以上描述可知,致密玄武巖與氣孔玄武巖和玄武質(zhì)火山碎屑巖相比,呈現(xiàn)高密度、低聲波時差、高電阻率的特征;粗面質(zhì)火山碎屑巖與粗面巖相比,呈現(xiàn)低電阻率特征。玄武巖與粗面巖相比,呈現(xiàn)低自然伽馬、低電阻率特征。輝綠巖呈低自然伽馬、高密度、高電阻率特征?；鸪蓭r巖性不同,測井曲線之間的組合特征也不同。

圖3 致密玄武巖的測井曲線特征

圖4 氣孔玄武巖的測井曲線特征

圖5 玄武質(zhì)火山碎屑巖的測井曲線特征

圖6 粗面巖的測井曲線特征

圖7 粗面質(zhì)火山碎屑巖的測井曲線特征

圖8 輝綠巖的測井曲線特征

3 火成巖巖性識別劃分

遼河盆地東部凹陷火成巖儲層巖性以玄武巖、粗面巖和輝綠巖為主,針對研究區(qū)實際需求和測井識別的可行性,將火成巖分為致密玄武巖、非致密玄武巖、粗面巖、非致密粗面巖、輝綠巖5類進行識別。其中非致密玄武巖為氣孔玄武巖以及包含玄武質(zhì)火山碎屑巖的玄武質(zhì)巖性的統(tǒng)稱,非致密粗面巖指的是包含粗面質(zhì)火山碎屑巖的粗面質(zhì)巖性的統(tǒng)稱。根據(jù)火成巖的巖性特征,總結(jié)出了該地區(qū)火成巖測井響應特征(表1),不同火成巖巖性測井響應特征不同。

3.1 交會圖版識別巖性

選取典型火成巖段,提取數(shù)據(jù)形成一套適合本區(qū)的火成巖識別圖版(圖9)。由γ-RLLD圖版可以準確識別輝綠巖、玄武巖類和粗面巖類,φCNL-ρ圖版可以區(qū)分致密玄武巖和非致密玄武巖,γ-RLLD和ΔtAC-γ圖版組合可以較為準確地識別和劃分粗面巖和粗面質(zhì)火山碎屑巖。

表1 遼河盆地東部凹陷火成巖測井響應特征

3.2 最小二乘支持向量機

3.2.1 最小二乘支持向量機原理及流程

支持向量機(SVM)[16-20]在解決小樣本、非線性及高維模式識別中表現(xiàn)出許多特有的優(yōu)勢,并能夠推廣應用到巖性分類識別等其它方面。

在特征空間Rn中,給定有N個樣本{(x1,y1),(x2,y2),…,(xN,yN)}的訓練集合,其中第i個輸入數(shù)據(jù)xi∈Rn,第i個輸出數(shù)據(jù)yi∈{-1,+1}是類別。支持向量機模型的目標是構(gòu)造一個如下形式的函數(shù):

(1)

式中:ω為超平面的法向量;b是超平面的截距;φ(x)為可將數(shù)據(jù)映射成線性可分的函數(shù)。

圖9 巖性識別交會圖版

最小二乘支持向量機(LS-SVM)是支持向量機的一種擴展形式,其損失函數(shù)采用誤差的二范數(shù),并用等式約束代替不等式約束,從而將二次規(guī)劃問題轉(zhuǎn)化為求解線性方程組的問題。

對應的目標函數(shù)為:

(2)

約束條件為:

(3)

其中,C≥0,為可調(diào)參數(shù),控制對誤差ε的懲罰程度,稱作懲罰因子。

用拉格朗日(Lagrange)方法解決這個優(yōu)化問題:

(4)

其中,α=[α1,α2,…,αN]為Lagrange算子,根據(jù)優(yōu)化條件得到:

(5)

其中,i=1,2,…,N。根據(jù)Mercer條件定義核函數(shù):

(6)

本文選擇高斯徑向基核函數(shù)

(7)

作為核函數(shù)。懲罰因子C和核函數(shù)寬度參數(shù)σ利用耦合模擬退火算法以及下山單純形法(Nelder-Mead)進行確定。最終得到?jīng)Q策函數(shù):

(8)

具體流程如圖10所示。

圖10 最小二乘支持向量機流程

3.2.2 最小二乘支持向量機巖性識別

將巖心/巖屑薄片深度進行歸位,與測井深度數(shù)據(jù)對應,并且薄片取自巖性較穩(wěn)定的井段;把薄片鑒定的結(jié)果作為巖性定名的標準。7口測試井最小二乘支持向量機的巖性識別與巖性樣本對比結(jié)果(表2)表明,不同井符合率有所差別,最高88.9%,最低76.1%,平均為79.5%,整體上符合較好。

表2 LS-SVM巖性測井識別符合率統(tǒng)計結(jié)果

4 結(jié)束語

1) 結(jié)合巖心和薄片特征,詳細描述和總結(jié)了研究區(qū)巖性的測井曲線特征。玄武巖低伽馬、低電阻率,致密玄武巖中高密度、中高中子,低聲波時差,氣孔玄武巖和玄武質(zhì)火山碎屑巖中低密度、中高中子和聲波時差;粗面巖高伽馬、高電阻率、中低密度和聲波時差、中子變化較大,粗面質(zhì)火山碎屑巖中伽馬、中低電阻率和密度、中高聲波時差和中子;輝綠巖中低伽馬和聲波時差、中高電阻率、中高密度和中子。

2) 利用最小二乘支持向量機法建立火成巖巖性識別模型,對遼河盆地東部凹陷火成巖巖性進行識別。識別出玄武巖、非致密玄武巖、粗面巖、非致密粗面巖、輝綠巖5種火成巖巖性,識別率較高;同時選取具有典型特征的穩(wěn)定數(shù)據(jù)段的層數(shù)據(jù),形成了一套適合本區(qū)的巖性圖版;最小二乘支持向量機程序識別與交會圖版結(jié)合應用,經(jīng)過實際生產(chǎn)驗證,能夠較好地解決本區(qū)巖性識別問題。

[1] 于小健.遼河盆地東部凹陷火山巖地層序列及其對儲層分布的影響研究[D].長春:吉林大學,2013 Yu X J.The research on volcanic sequence stratigraphy and its effect on reservoir in Eastern Sag of Liaohe Basin,NE China[D].Changchun:Jinlin University,2013

[2] 李軍.遼河坳陷東部凹陷沙三段沉積儲層特征與油氣關(guān)系[D].北京:中國地質(zhì)大學,2010 Li J.Reservoir characteristics and oil/gas relationship in the Es3sedimentation of the eastern depression of the Liaohe Rift[D].Beijing:China University of Geosciences,2010

[3] 陳全茂,李忠飛.遼河盆地東部凹陷構(gòu)造及其含油氣性分析[M].北京:地質(zhì)出版社,1998:12-14 Chen Q M,Li Z F.Structure and capability analysis of the Eastern Sag of Liaohe Basin[M].Beijing:Geological Publishing House,1998:12-14

[4] 廖興明,姚繼峰,于天欣.遼河盆地構(gòu)造演化與油氣[M].北京:石油工業(yè)出版社,1996:10-12 Liao X M,Yao J F,Yu T X.Tectonic evolution and oil and gas of Liaohe Basin[M].Beijing:Petroleum Industry Press,1996:10-12

[5] 趙澄林,孟衛(wèi)工,金春爽,等.遼河盆地火山巖與油氣[M].北京:石油工業(yè)出版社,1999:15-20 Zhao C L,Meng W G,Jin C S,et a1.Volcanic rocks and hydrocarbon in Liaohe Basin[M].Beijing:Petroleum Industry Press,1999:15-20

[6] 于小健,黃玉龍,李軍,等.遼河盆地東部凹陷古近系沙三段粗面巖成因及噴發(fā)模式[J].中國礦業(yè),2013,22(3):58-65 Yu X J,Huang Y L,Li J.Lithogenesis and eruption model of the eocene trachyte in the eastern sag of Liaohe basin[J].China Mining Magazine,2013,22(3):58-65

[7] 孫卉.遼河東部凹陷中新生代火山事件與構(gòu)造作用[J].中國石油勘探,2005(5):36-41 Sui H.Mesozoic-Cenozoic volcanic events and structure effects in the Eastern Sag in Liaohe depression[J].China Petroleum Exploration,2005(5):36-41

[8] 李思偉,王璞珺,丁秀春,等.遼河東部凹陷走滑構(gòu)造及其與火山巖分布的關(guān)系[J].地質(zhì)論評,2014,60(3):591-600 Li S W,Wang P J,Ding X C,et al.The relationship of strike-slip belts and volcanic rocks distribution in the eastern depression Liaohe basin[J].Geological Review,2014,60(3):591-600

[9] 蔡國鋼.遼河裂谷東部凹陷粗面巖成因機理探討[J].地球?qū)W報,2010,31(2):245-250 Cai G G.Genetic mechanism of trachyte in the eastern depression of Liaohe rift[J].Acta Geoscientica Sinica,2010,31(2):245-250

[10] 馮玉輝,黃玉龍,丁秀春,等.遼河盆地東部凹陷中基性火山巖相地震響應特征及其機理探討[J].石油物探,2014,53(2):206-215 Feng Y H,Huang Y L,Ding X C,et al.Relationship between intermediate-mafic volcanic facies and their corresponding seismic reflections:a case study from eastern depression of Liaohe basin[J].Geophysical Prospecting for Petroleum,2014,53(2):206-215

[11] 王巖泉,胡大千,蔡國剛,等.遼河盆地東部凹陷火山巖儲層特征與主控因素[J].石油學報,2013,34(5):896-904 Wang Y Q,Hu D Q,Cai G G,et al.Characteristics and controlling factors of cenozoic volcanic reservoirs in Liaohe basin,NE China[J].Acta Petrolei Sinica,2013,34(5):896-904

[12] 黃玉龍,單俊峰,邊偉華,等.遼河坳陷中基性火成巖相分類及儲集意義[J].石油勘探與開發(fā),2014,41(6):1-10 Huang Y L,Chan J F,Bian W H,et al.Facies classification and reservoir significance of the cenozoic intermediate and mafic igneous rocks in Liaohe depression,east China[J].Petroleum Exploration and Development,2014,41(6):1-10

[13] 楊光達.遼河坳陷東部凹陷中南部沙三段火山巖特征分析[J].石油天然氣學報,2011,33(8):54-57 Yang G D.Analysis of Es3volcanic rock property in the eastern sag of Liaohe depression[J].Journal of Oil and Gas Technology,2011,33(8):54-57

[14] 張占文,陳振巖,蔡國剛,等.遼河坳陷火成巖油氣藏勘探[J].中國石油地質(zhì),2005,10(4):16-22 Zhang Z W,Chen Z Y,Cai G G,et al.Exploration of volcanic oil and gas reservoir in Liaohe depression[J].China Petroleum Exploration,2005,10(4):16-22

[15] 劉嘉麒,孟凡超,崔巖,等.試論火山巖油氣藏成藏機理[J].巖石學報,2010,26(1):1-13

Liu J Q,Meng F C,Cui Y,et al.Discussion on the formation mechanism of volcanic oil and gas reservoir[J].Acta Petrologica Sinica,2010,26(1):1-13

[16] 于代國,孫建孟,王煥增,等.測井識別巖性新方法—支持向量機方法[J].大慶石油地質(zhì)與開發(fā),2005,24(5):93-95 Yu D G,Sun J M,Wang H Z,et al.A new method for logging lithology identification-SVM[J].Petroleum Geology & Oilfield Development in Daqing,2005,24(5):93-95

[17] 宋延杰,張劍風,閆偉林,等.基于支持向量機的復雜巖性測井識別方法[J].大慶石油學院學報,2007,31(5):18-20 Song Y J,Zhang J F,Yan W L,et al.A new identification method for complex lithology with support vector machine[J].Ournal of Daqing Petroleum Institute,2007,31(5):18-20

[18] 鐘儀華,李榕.基于主成分分析的最小二乘支持向量機巖性識別方法[J].測井技術(shù),2009,33(5):425-429 Zhong Y H,Li R.Application of principal component analysis and least square support vector machine to lithology identification[J].Well Logging Technology,2009,33(5):425-429

[19] 韓學輝,盧時林,支樂菲,等.應用最小二乘支持向量機識別J13井區(qū)杜家臺油層巖性[J].特種油氣藏,2011,18(6):18-21 Han X H,Lu S L,Zhi L F,et al.Identification of the lithology of Dujiatai formation in J13 well block with least squares support vector machine[J].Special Oil and Gas Reserviors,2011,18(6):18-21

[20] 韓學輝,支樂菲,劉榮,等.應用最小二乘支持向量機識別廣利油田沙四段儲層巖性[J].地球物理學進展,2013,28(4):1886-1892 Han X H,Zhi L F,Liu R,et al.A lithologic identification method in Es4reservoir of Guangli oil field with least square support vector machine[J].Progress in Geophysics,2013,28(4):1886-1892

(編輯：陳 杰)

The logging characteristics of the intermediate and basic igneous rock from the Eastern Depression of Liaohe Basin

Liu Zongli1,Wang Zhuwen1,Liu Jinghua1,Feng Yuhui2,Qu Dong3

(1.CollegeofGeoExplorationScienceandTechnology,JilinUniversity,Changchun130026,China；2.CollegeofEarthSciences,JilinUniversity,Changchun130061,China; 3.WellLoggingCompany,GWDC,Panjin124011,China)

Since the Cenozoic,because of the influence of the Tanlu Fault,six tectonic movements have occurred in the Eastern Depression of Liaohe Basin,which were accompanied by volcanic activity with different intensity; so igneous rock is widely developed.Based on six wells,the characteristics of log combination of the intermediate and basic igneous rock are analyzed,including tight basalt,stomatal basalt,basaltic volcanic clastic rock,trachyte,trachytic volcanic clastic rock and diabase.Selecting typical lithology data,the characteristics of igneous lithology logging response in this area are summarized; and a set of suitable lithology identification chart is established for the area.The crossplot of natural gamma versus deep lateral resistivity can well identify diabase,and the crossplot of neutron and density can well distinguish basalts.Using least squares support vector machine (LS-SVM),the igneous rock lithology identification model is constructed.Compared with core analysis and cutting description,the accuracy rate of the identification result reaches up to 79.5%.The combination of crossplot and LS-SVM well solves the problem of lithology identification and provides assistance for further interpretation of oil and gas reservoirs.

Eastern Depression of Liaohe Basin,intermediate and basic igneous rock,well logging curve,crossplot,lithology identification,least squares support vector machine (LS-SVM)

2015-04-01;改回日期：2015-06-26。

劉宗利(1990—),男,博士在讀,主要從事地球物理測井新方法、新技術(shù)研究工作。

王祝文(1961—),男,教授,博士生導師,主要從事地球物理測井、核地球物理、測井新方法新技術(shù)等方面的解釋理論和方法的教學與科研工作。

國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)項目(2012CB822002)和中國石油天然氣股份有限公司科學研究與技術(shù)開發(fā)項目(2012E3001)共同資助。

P631

A

1000-1441(2015)06-0787-09

10.3969/j.issn.1000-1441.2015.06.018