张家垛油田阜三段湖相滩坝砂储层特征
2015-03-24丁晓琪张哨楠韩玫梅刘曦翔湛小红张咏梅
杨 鹏,丁晓琪,张哨楠,韩玫梅,刘曦翔,湛小红,张咏梅
(1.西南石油大学 油气藏地质及开发工程国家重点实验室,四川 成都 610500; 2.西南石油大学 地球科学与技术学院,四川 成都 610500; 3.成都理工大学 能源学院,四川 成都 610059)
张家垛油田阜三段湖相滩坝砂储层特征
杨 鹏1,2,丁晓琪1,3,张哨楠1,2,韩玫梅1,2,刘曦翔1,湛小红1,张咏梅1
(1.西南石油大学 油气藏地质及开发工程国家重点实验室,四川 成都 610500; 2.西南石油大学 地球科学与技术学院,四川 成都 610500; 3.成都理工大学 能源学院,四川 成都 610059)
湖相滩坝砂往往毗邻烃源岩,在高源-储压差条件下易聚集成藏。近年来在我国东部新生代地层中发现了大量的此类油气藏,成为勘探的热点。以张家垛油田阜三段为例,通过岩心观察、研究滩坝砂的分布规律;利用压汞资料、铸体薄片及扫描电镜,分析不同成因储层的孔喉分布和孔隙类型组合。研究结果表明:①滩坝砂中发育3种岩相,4种岩相组合类型,对应于4种测井相;②滩坝砂具有粒度细、单层厚度小、滩坝间互、叠合连片的特征,坝砂呈土豆状平行于湖岸线展布;③坝砂以细砂岩为主,杂基含量低;原生孔隙和次生孔隙发育,粘土中富高岭石而贫绿泥石和混层粘土;滩砂以粉砂-极细砂为主,杂基含量高,孔隙结构差,次生孔隙欠发育,粘土中富混层粘土而贫高岭石和绿泥石。
滩坝砂;储层;成岩作用;阜三段;张家垛油田;苏北盆地
湖相滩坝砂是陆相盆地中重要的一种储层类型[1],与其他类型的储层相比,滩坝砂具有以下特殊性:①粒度偏细,往往为细砂岩、粉砂岩,粒度粗的滩砂少见;②单层厚度小,多层可以纵向上相互叠置,分布规律较为复杂[2-3];③由于粘土含量高,储层物性较差,渗透率偏低;④滩坝砂往往与半深湖—深湖相的烃源岩伴生,具有优先捕集油气的能力[4]。由于滩坝砂的上述特征,导致滩坝砂普遍含油但非均质性强,形成圈闭的隐蔽性较强[5-6],油气勘探开发难度大[7]。本文以苏北盆地张家垛油田阜三段(阜宁组三段)下砂岩段(Ⅲ砂组)为例,通过岩心观察及实验分析,研究不同成因砂体剖面组合、分布规律及储层特征,为张家垛油田阜三段湖相滩坝砂的勘探和开发提供地质依据。
张家垛油田位于苏北盆地海安凹陷西部曲塘次凹的北部陡坡带[8](图1),处于断裂的下降盘,是一个由边界断层控制的大型鼻状构造[9-11],为油气富集的有利区带,具有北陡南缓的特点。油气主要富集于阜三段[12-14]。阜三段沉积期,湖盆底形平缓、生物及藻类繁盛,北部的泰州凸起为低幅的水下隆起,未提供物源,物源主要来自北部的三角洲。在湖浪和沿岸流的作用下,将来自北部三角洲的砂体再搬运、沉积,形成浅湖滩坝砂[15-17]。张家垛油田目前主要的勘探目的层为阜三段的滩坝砂岩性油藏,阜三段自下而上划分为4个砂层组,其中下砂岩段砂体发育,储层段埋深2 800~3 800 m[18]。
1 滩坝砂剖面组合特征
通过10余口井的岩心观察,根据岩相划分标准,在阜三段滩坝砂中归纳总结出3种岩相:平行层理细砂岩相(Sh),泥质纹层+生物钻孔粉-细砂岩相(Fr),生物扰动粉砂岩相(Fb)。
1.1 平行层理细砂岩相
该岩相主要发育在细砂岩中,主要的层理类型为平行层理,有时局部可过渡为脉状层理,断面可见少许炭屑。细砂岩具有杂基含量低、分选好、磨圆好的特点,说明沉积水动力较强,碎屑颗粒经历了长距离的搬运,细砂岩中可见藻屑。
1.2 泥质纹层+生物钻孔粉-细砂岩相
该岩相主要发育在粉砂岩中,少量在极细粒砂岩中,主要的层理类型为沙纹层理。由于强烈的生物扰动,泥质纹层呈不连续状,层面中可见大量的炭屑及生物钻孔。粉-细砂岩中杂基含量高和强烈的生物扰动反映出水体较浅但能量较弱。分选好、磨圆好说明经过了长距离的搬运。
1.3 生物扰动粉砂岩相
该岩相主要发育在粉砂岩和泥岩中,强烈的生物扰动使该岩相的原有层理特征不明显,粉砂岩和泥岩高度混合。当生物扰动特征强烈时,形成“似块状”的层理。该岩相粒度细、生物扰动强烈代表水体浅且能量弱的沉积环境。
1.4 剖面组合特征
以上3种岩相在剖面上常形成4种剖面组合类型,对应4种测井相(图2),反映出不同水动力条件。
1) 剖面组合Ⅰ(指形)
自下而上为Fr—Sh—Fr,3种岩相的厚度小,累计厚度不超过2 m。Sh岩相的厚度一般小于0.5 m,测井曲线形态为指形。沉积物供给有限,水动能量较弱,该剖面组合平面分布较稳定,往往代表滩砂。
图1 苏北盆地张家垛油田构造位置及阜三段Ⅲ砂组顶面构造
2) 剖面组合Ⅱ(似箱形)
自下而上为Fb—Fr—Sh—Fr,与剖面组合Ⅰ相比,该岩相的底部往往发育Fb,且Sh的厚度可超过1 m,最厚达到3 m,整个剖面组合的厚度可超过4 m,测井曲线形态为似箱形。沉积物供给稳定,水动能量较强,平面上呈透镜状分布,为典型的坝砂。
3) 剖面组合Ⅲ(漏斗形)
自下而上为Fb—Fr—Sh,与剖面组合Ⅱ相比,该岩相顶部不发育Fr,Sh与黑色泥岩直接接触。Sh的厚度变化较大(0.2~2 m),测井曲线形态为漏斗形。自下而上水动力变强,Sh厚度小时为滩砂,反之是坝砂。
4) 剖面组合Ⅳ(齿化钟形)
自下而上为Fb—Fr的多个叠加,但向上Fb的厚度增加而Fr厚度减小,形成低幅的齿化钟形。沉积物供给不充分,水动力能量弱,往往代表滩缘沉积。
2 滩坝砂分布规律
在单井沉积微相识别的基础上,进行了沉积微相的对比。图3为研究区一条东西向的微相对比图,可以看出下砂岩段的砂体主要为坝砂,其次为滩砂,并且坝砂向两侧逐渐相变为滩砂。滩坝砂在横向上稳定,可以在数千米范围内进行追踪对比,整个下砂岩段具有单层厚度小、滩坝间互和叠合连片的特征。
在利用相标志确定沉积相类型的基础上,通过对岩心的观察,利用单井相及剖面沉积相的识别结果,结合从地震反演中提取出的波阻抗值,绘制了阜三段下砂岩段(Ⅲ砂组)各单砂体的沉积微相平面分布(图4)。通过编制5个单砂体的沉积微相图,可以看出下砂岩段5个单砂体滩砂沉积微相分布面积大,而砂坝在不同的层段分布面积差异较大。总体来看,坝砂呈透镜状东西向平行于湖岸线展布。
图2 苏北盆地张家垛油田阜三段典型测井曲线岩相解释
图3 苏北盆地张家垛油田阜三段Ⅲ砂组沉积微相剖面
3 滩坝砂储集空间特征
作为油气储层评价中的重要指标,储层储集空间特征一直受到地质工作学者的广泛重视并取得了众多成果[19-21]。随着研究的深入,储层的储集空间特征的描述更为精细,这里选取6个典型的滩坝砂岩样品进行毛细管压力曲线测定(图5),同时进行薄片分析、粘土X-衍射和扫描电镜分析(表1)。根据6个样品的毛管压力曲线可以明显分为3类储层。
图4 苏北盆地张家垛油田阜三段Ⅲ砂组单砂体沉积微相平面分布
图5 不同储层类型的毛管压力曲线及孔喉分布
表1 粒度参数与粘土含量
Table 1 Parameters of particle size and content of clay
井号深度/m储层类型粒度分析X-衍射平均粒度/mm分选系数伊/蒙混层/%高岭石/%绿泥石/%伊利石/%张3X13634.4Ⅰ类0.181.240291318张32958.0Ⅱ类0.122.150182012张1012939.6Ⅱ类0.091.86511168张1012945.8Ⅲ类0.092.263131212苏2582750.8Ⅲ类0.082.37041313张1012945.3Ⅲ类0.072.67521211
3.1 Ⅰ类储层
Ⅰ类储层样品孔隙度大于15%,渗透率大于8×10-3μm2。喉道分选系数约2.4,排驱压力较低,一般小于0.5 MPa;中值压力在1 MPa以下,最小非饱和体积小于10%。孔喉的分布为非正态分布,粗孔喉占有较大的比例,中孔喉和细孔喉较少,次生孔隙发育。这类储层的微观特征表现为粒度较粗,主要为细砂岩,少量为粗粉砂岩,杂基含量较低。孔隙类型以粒间孔、粒内溶孔和铸模孔为主,微孔较少(图6a,b),其中次生孔隙占总孔隙的比例可达到1/3~1/2,甚至更高。由于溶蚀孔占有较大的比例,导致大孔隙多,喉道较粗。
由于杂基含量低,导致粘土矿物中伊/蒙混层很少,常见的粘土矿物为绿泥石和高岭石(图6b)。因此这种类型的储层存在一定的速敏,但水敏不强。根据沉积微相研究,认为这种砂体为坝砂成因,具有单层厚度大的特点,单层厚度超过3 m。结合该地区的试油成果,认为这种类型的储层具有自然产能,在不压裂的情况下即可获取工业油流。
3.2 Ⅱ类储层
Ⅱ类储层样品孔隙度在10%~15%,渗透率在(3~8)×10-3μm2。喉道分选系数为2.74,排驱压力在0.5~1 MPa,中值压力明显增加,在5 MPa左右,最小非饱和体积小于20%。孔喉仍非正态分布,有一部分孔喉偏粗,但比例不大,而细喉和微喉占有较大的比例,孔隙结构差。
图6 张家垛油田不同类型储层典型显微照片
这类储层的微观特征表现为粒度偏细,以泥质细砂岩和粉砂岩为主,杂基含量较高。孔隙类型有粒间孔、溶蚀孔和微孔,其中微孔占有较大的比例,可达20%~30%。由于高杂基含量,导致粘土矿物中的伊/蒙混层含量较高,且含有不同程度的高岭石和绿泥石(图6c,d)。因此,这种类型的储层既有水敏又有酸敏,但以水敏为主。结合沉积微相研究,认为这种砂岩为滩砂成因,单层厚度相对较小,主要分布在坝砂的周缘。从已有的试油试采资料来看,这种类型的储层不具有自然产能,需要对其进行压裂,改善储层的渗流能力,方可获取工业产能。
3.3 Ⅲ类储层
Ⅲ类储层样品孔隙度一般小于10%,渗透率小于3×10-3μm2。喉道分选系数为5左右,排驱压力迅速上升至2~4 MPa,中值压力达到50 MPa左右。从孔喉分布来看,基本不存在大的孔喉,全为微孔和细孔,孔喉明显较前两种储层偏细。
这种类型储层的微观特征表现为粒度细,以粉砂岩和泥质粉砂岩为主,杂基含量非常高,反应出水体能力弱,淘洗不干净。孔隙类型主要为微孔(图6e,f),其次为少量的粒间孔,溶蚀孔不发育。
由于杂基含量非常高,不利于孔隙水的流动,所以长石的溶蚀变弱,高岭石含量低[22]。这种类型储层中的粘土矿物几乎全为伊/蒙混层,水敏非常强。结合沉积微相的研究成果,认为这种类型的储层主要为滩砂及滩间沉积,弱的水动力导致粒度细,粘土含量高。目前尚未从这种类型储层获得工业产能,但从其它油田的勘探情况来看,不排除通过大规模的压裂,获得低产工业油流的可能。
4 结论
1) 坝砂主要发育平行层理细砂岩相,其次是泥质纹层+生物钻孔粉-细砂岩相,单层厚度往往超过2 m;而滩砂主要发育泥质纹层+生物钻孔粉-细砂岩相和平行层理细砂岩相,其次为生物扰动粉砂岩相,单层厚度一般小于2 m。
2) 坝砂分布面积较小,呈土豆状平行于湖岸线分布,研究区北部较南部发育。滩、坝砂体单层厚度小、滩坝间互、叠合连片。
3) 坝砂为水动力较强的产物,粒度一般为细砂岩,杂基含量低,原生孔隙和次生孔隙发育,孔隙结构好。坝砂向两侧均过渡为滩砂,形成滩砂的水动力相对较弱,粒度往往为粉-极细粒砂岩,原生孔隙发育而次生孔隙不发育,孔隙结构差、具有孔隙中等、渗透率差的特征。
[1] 操应长,王艳忠,徐涛玉,等.东营凹陷西部沙四上亚段滩坝砂体有效储层的物性下限及控制因素[J].沉积学报,2009,27(2):230-237. Cao Yingchang,Wang Yanzhong,Xu Taoyu,et al.The petrophysical parameter cutoff and controlling factors of the effective reservoir of beach and bar sandbodies of the upper part of the fourth member of the Shahejie Formation in west part of Dongying Depression[J].Acta Sedimentologica Sinica,2009,27(2):230-237.
[2] 邓宏文,马立祥,姜正龙,等.车镇凹陷大王北地区沙二段滩坝成因类型、分布规律与控制因素研究[J].沉积学报,2008,26(5):715-724. Deng Hongwen,Ma Lixiang,Jiang Zhenglong,et al.Sand bank generation types and its controls on their distribution,the second member of Shahejie Formation,Lower Tertiary,Dawangbei,Chezhen Depression[J].Acta Sedimentologica Sinica,2008,26(5):715-724.
[3] 刘迪,丁晓琪,韩玫梅.异常高压对储层成岩作用的控制[J].西南石油大学学报(自然科学版),2013,35(3):39-46. Liu Di,Ding Xiaoqi,Han Meimei.Govern of abnormal pressure on reservoir diagenesis[J].Journal of Southwest Petroleum University(Science & Technology Edition),2013,35(3):39-46.
[4] 邓宏文,高晓鹏,赵宁,等.济阳坳陷北部断陷湖盆陆源碎屑滩坝成因类型、分布规律与成藏特征[J].古地理学报,2010,12(6):737-747. Deng Hongwen,Gao Xiaopeng,Zhao Ning,et al.Genetic types,distribution patterns and hydrocarbon accumulation in terrigenous beach and bar in northern faulted-lacustrine-basin of Jiyang Depression[J].Journal of Palaeogeography,2010,12(6):737-747.
[5] 罗红梅,朱毅秀,穆星,等.渤海湾渤南洼陷深层湖相滩坝储集层沉积微相预测[J].石油勘探与开发,2011,38(2):182-190. Luo Hongmei,Zhu Yixiu,Mu Xing,et al.Seismic facies prediction of lacustrine beach and bar reservoirs in the deep zone of the Bonan Subsag,Bohai Bay Basin[J].Petroleum Exploration and Development,2011,38(2):182-190.
[6] 郭松,谭丽娟,林承焰,等.储层非均质性对油气成藏的影响——以博兴油田沙四上亚段滩坝相砂岩为例[J].石油实验地质,2013,35(5):534-538. Guo Song,Tan Lijuan,Lin Chengyan,et al.Influence of reservoir bed heterogeneity on hydrocarbon accumulation:Taking beach-bar sandstone of upper fourth member of Eocene Shahejie Formation in Boxing oilfield as an example[J].Petroleum Geology & Experiment,2013,35(5):534-538.
[7] 栾春华.湖相滩坝砂岩油藏储层预测研究[J].中国科技信息,2011,24:42-46. Luan Chunhua.The study of prediction of lacustrine beach bar sandstone reservoir[J].China Science and Technology in Formation,2011,24:42-46.
[8] 刘曦翔,杨正茂,丁晓琪,等.张家垛油田阜三段异常压力特征分析[J].断块油气田,2014,21(2):191-195. Liu Xixiang,Yang Zhengmao,Ding Xiaoqi,et al.Analysis on abnormal pressure characteristics of the third member of Funing Formation in Zhangjiaduo Oilfield[J].Fault-block Oil & Gas Field,2014,21(2):191-195.
[9] 魏祥峰,张廷山,梁兴,等.海安凹陷曲塘、海北次凹戴南组-三垛组一段沉积特征及演化[J].石油与天然气地质,2012,33(2):265-276. Wei Xiangfeng,Zhang Tingshan,Liang Xing,et al.Sedimentary characteristics and evolution of the 1st member of Dainan Formation and the 1stmember of Sanduo Formation in Qutang and Haibei Sub 2sags,Haian Sag[J].Oil & Gas Geology,2012,33(2):265-276.
[10] 焦里力.苏北盆地张家垛油田阜三段油气富集规律[J].石油与天然气地质,2012,33(2):166-173. Jiao Lili.Hydrocarbon enrichment patterns of the 3rd member of Funing Formation in Zhangjiaduo oilfield,Subei Basin[J].Oil & Gas Geology,2012,33(2):166-173.
[11] 鲁雪松,马淑芳,杨贵丽.苏北盆地金湖凹陷三垛组剥蚀厚度恢复及意义[J].西南石油大学学报(自然科学版),2009,31(4):13-17. Lu Xuesong,Ma Shufang,Yang Guili.Erosion thickness reconstruction of Sanduo Formation and its significance in Jinhu Sag,Subei Basin[J].Journal of Southwest Petroleum University(Science & Technology Edition),2009,31(4):13-17.
[12] 徐辉,田荣斌,李珊,等.张家垛阜三段储层敏感性和油水相渗特征[J].油气藏评价与开发,2011,01(6):48-51. Xu Hui,Tian Furong,Li Shan,et al.Reservoir sensitivity and oil-water relative permeability of the third member of Funing formation in Zhangjiaduo[J].Reservoir Evaluation and Development,2011,01(6):48-51.
[13] 黄方方,丁晓琪,万友利,等.张家垛油田阜三段储层成岩作用[J].断块油气田,2013,20(2):150-153. Huang Fangfang,Ding Xiaoqi,Wan Youli,et al.Diagenesis of Fu 3 Member reservoir in Zhangjiaduo Oilfield[J].Fault-block Oil & Gas Field,2013,20(2):150-153.
[14] 乐锦鹏,张哨楠,丁晓琪,等.苏北盆地张家垛油田阜三段储层特征及主控因素[J].东北石油大学学报,2013,37(1):71-77. Le Jinpeng,Zhang Shaonan,Ding Xiaoqi,et al.Reservoir characteristics of block 3 of Funing formation in Zhangjiaduo oilfield,Subei Basin and their main control factors[J].Journal of Northeast Petroleum University,2013,37(1):71-77.
[15] 冯胜斌,牛小兵,刘飞,等.鄂尔多斯盆地长7致密油储层储集空间特征及其意义[J].中南大学学报(自然科学版),2013,11:4574-4580. Feng Shengbin,Niu Xiaobing,Lu Fei,et al.Characteristics of Chang7 tight oil reservoir space in Ordos Basin and its significance[J].Journal of Central South University(Science and Technology),2013,11:4574-4580.
[16] 李国斌,姜在兴,陈诗望,等.利津洼陷沙四上亚段滩坝沉积特征及控制因素分析[J].中国地质,2008,35(5):911-921. Li Guobin,Jiang Zaixing,Chen Shiwang,et al.Sedimentary characteristics and controlling factors of beach bars in the Upper Submember of the Fourth Member of the Shahejie Formation in the Lijin subasin[J].Geology inChina,2008,35(5):911-921.
[17] 田继军,姜在兴.惠民凹陷与东营凹陷沙四上亚段滩坝沉积特征对比与分析[J].吉林大学学报,2012,42(3):612-623. Tian Jijun,Jiang Zaixing.Comparison and analysis of beach bars sedimentary characteristics of Upper Es4in Huimin and Dongying Depression[J].Journal of Jilin University(Earth Science Edition),2012,42(3):612-623.
[18] 陈倩倩,丁晓琪,万友利,等.张家垛油田阜三段Ⅲ砂组沉积微相特征与沉积模式[J].石油与天然气地质,2013,34(6):751-757. Chen Qianqian,Ding Xiaoqi,Wan Youli,et al.Microfacies characteristics and sedimentary models of the third member of Funing Formation in Zhangjiaduo oilfield[J].Oil & Gas Geology,2013,34(6):751-757.
[19] 俞亮亮.苏北盆地张家垛油田阜三段储层特征研究[J].中国石油和化工标准与质量,2013,5:164-165. Yu Liangliang.The study of reservoir characteristics in Fu III Interval,Zhangjiaduo Oilfield,Subei Basin[J].China Petroleum and Chemical Standard and Quality,2013,5:164-165.
[20] 张春明.山东惠民凹陷临南地区古近系沙三下亚段沉积相与储集性能研究[D].北京:中国地质大学,2012. Zhang Chunming.The sedimentary facies and reservoir characteristics of the lower part of the third member of Paleogene Shahejie Formation,in Linnan area,Huimin depression,Shandong province[D].Beijng:China University of Geosciences,2012.
[21] 赵军龙,闫博,赵靖舟,等.陇东环县油区M井区长8储集层特征及单井产能预测[J].石油与天然气地质,2013,34(5):694-699. Zhao Junlong,YanBo,ZhaoJingzhou,et al.Reservoir characteristics of Chang 8 and single-well productivity prediction of M wellblock in Huanxian oil district,Longdong area[J].Oil & Gas Geology,2013,34(5):694-699.
[22] 张哨楠,丁晓琪,万友利,等.致密碎屑岩中粘土矿物的形成机理与分布规律[J].西南石油大学学报(自然科学版),2012,34(3):174-182. Zhang Shaonan,Ding Xiaoqi,Wan Youli,et al.Formation mechanism and distribution of clay Minerals of deeply tight siliciclastic reservoirs[J].Journal of Southwest Petroleum University(Science & Techno-logy Edition),2012,34(3):174-182.
(编辑 张玉银)
Characteristics of lacustrine beach-bar sandstone reservoirs in the 3rdMember of the Funing Formation in Zhangjiaduo oilfield
Yang Peng1,2,Ding Xiaoqi1,3,Zhang Shaonan1,2,Han Meimei1,2,Liu Xixiang1,Zhan Xiaohong1,Zhang Yongmei1
(1.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,SouthwestPetroleumUniversity,Chengdu,Sichuan610500,China;2.SchoolofGeoscienceandTechnology,SouthwestPetroleumUniversity,Chengdu,Sichuan610500,China;3.CollegeofEnergyResources,ChengduUniversityofTechnology,Chengdu,Sichuan610059,China)
Lacustrine beach-bar sandstone is generally adjacent to source rocks,and it is favorable for hydrocarbon accumulation under high differential pressure between the source rock and reservoir.In recent years,a large number of oil and gas reservoirs have been discovered in the Cenozoic lacustrine beach bar sandstones in eastern China.Distribution of beach-bar sands in the 3rdMember of the Funing Formation was studied on the basis of the core observation.Mercury data,cast thin sections and scanning electron microscope were integrated to analyze the distribution of pore throat and the combinations of porosity types in reservoirs of different genetic types.The following results were obtained.①Three kinds of lithofacies were recognized in beach-bar sand and they can be organized into four kinds of lithofacies associations,which correspond to four kinds of logging facies;②The beach-bar sand has the characteristics of fine grain size,thin single layer,as well as interbedding,overlapping and lateral connection of the beach and barsands,and the potato-shaped bar sands are distributed in parallel with the lake shoreline;③The bar sands are dominated by fine sandstone with low matrix content;Primary pore and secondary pore are well developed,and clay minerals are rich in kaolinite and poor in chlorite and mixed layer;the beach sands are generally dominated by siltstone and very fine sandstone with high matrix content;Their pore structures are poor and secondary pores are underdeveloped,and the clay minerals are rich in mixed layer and poor in kaolinite and chlorite.
beach bar sand,reservoir,diagenesis,3rdMember of the Funing Formation,Zhangjiaduo oilfield,northern Jiangsu Basin
2014-12-23;
2015-03-25。
杨鹏(1986—),男,博士研究生,油气藏理论与技术(含油气经济评价)。E-mail:doctoryangpeng@qq.com。
国家自然科学基金项目(41302115);博士后科学基金项目(2012M511941);国家自然科学青年基金项目(41302115)。
0253-9985(2015)03-0456-07
10.11743/ogg20150314
TE122.2
A