资源描述:
2 0 1 4 年第3 8 卷 第2 期 中国石油大学学报 自然科学版 J o u r n a lo fC h i n aU n i v e r s i t yo fP e t r o l e u m V 0 1 .3 8N o .2 A p r .2 0 1 4 文章编号1 6 7 3 - 5 0 0 5 2 0 1 4 0 2 - 0 0 9 9 4 3 4d o i 1 0 .3 9 6 9 /j .i s s n .1 6 7 3 - 5 0 0 5 .2 0 1 4 .0 2 .0 1 5 天然气水合物藏注热水开采敏感因素试验研究 李淑霞1 ’2 ,李杰1 ,徐新华1 ,李小森2 1 .中国石油大学石油工程学院,山东青岛2 6 6 5 8 0 ;2 .中国科学院可再生能源与天然气水合物重点实验室,广东广州5 1 0 6 4 0 摘要采用自制的一维天然气水合物 N G H 开采模拟试验装置,模拟海洋地质条件,在填砂模型 填砂管直径8c m , 长8 0c m 中生成天然气水合物,并通过注入热盐水进行热力开采的物理模拟试验,分析水合物藏地质因素和注热 参数对注热水开采能量效率 水合物分解所得甲烷气的热量与注入热量之比 的影响。结果表明影响能量效率的 因素由大到小依次为注热水温度、水合物饱和度、初始温度、注热水时间、注热水速度;在水合物藏地质因素中,水合 物饱和度越大,初始温度越高,注热水开采能量效率越高;在注热参数中,注热温度越高,能量效率越低;在注热开采 时需要合理优化注热水温度、注热水时间和注热水速度;试验条件下,注热水开采能量效率最高的试验组合为水合 物饱和度4 8 %,初始温度5 ℃,注热水温度4 0 ℃,注热水时间3 5 0m i n ,注热水速度1 2m L /m i n ,最高能量效率为 6 .7 4 。 关键词天然气水合物;注热盐水;水合物分解;能量效率 中图分类号T E3 7 5文献标志码A E x p e r i m e n t a ls t u d yo ni n f l u e n c i n gf a c t o r sf o rh y d r a t e d i s s o c i a t i o ni nah o tb r i n ei n j e c t i o np r o c e s s L IS h u .x i a l ”,L IJ i e l ,X UX i n .h u a l ,L IX i a o .s e n 2 1 .S c h o o lo fP e t r o l e u mE n g i n e e r i n gi nC h i n aU n i v e r s i t yo fP e t r o l e u m ,Q i n g d a o2 6 6 5 8 0 ,C h i n a ; 2 .K e yL a b o r a t o r yo fR e n e w a b l eE n e r g ya n dG a sH y d r a t e ,C h i n e s eA c a d e m yo fS c i e n c e s ,G u a n g z h o u51 0 6 4 0 ,C h i n a A b s t r a c t 7 r h ef o r m a t i o na n dd i s s o c i a t i o no fn a t u r a lg a sh y d r a t e N G H w a ss t u d i e du s i n gas p e c i a l l yd e s i g n e do n e d i m e n s i o n a le x p e r i m e n t a ls e t u pt os i m u l a t et h eh y d r a t eb e a r i n gs e d i m e n tb e l o wt h es e a b e d .N G HW a sf o r m e di nas a n d p a e km o d e l w i t h ad i a m e t e ro f8c ma n dl e n g ho f8 0c m .a n dt h e nh o tb r i n ew a si n j e c t e di n t ot h es a n d p a c kt od i s s o c i a t et h eg a sh y d r a t e .T h e i m p o r t a n tf a c t o r st h a ti n f l u e n c i n gt h ee n e r g ye f f i c i e n c yf o rh y d r a t ed i s s o c i a t i o nw e r ei n v e s t i g a t e d ,i n c l u d i n gt h et e m p e r a t u r eo f t h ei n j e c t e dh o tb r i n e ,t h es a t u r a t i o na n di n i t i a lt e m p e r a t u r eo fh y d r a t ef o r m e d ,a n dt h ei n j e c t i o nt i m ea n dr a t eo ft h eh o tb r i n e . T h ee n e r g ye f f i c i e n c yw a sd e f i n e da st h er a t i oo ft h eh e a to ft h em e t h a n eg a sr e l e a s e dt ot h ei n j e c t i o nh e a t .T h ee x p e r i m e n t a lr e s u i t ss h o wt h a t ,t h eo r d e rf r o mb i gt os m a l lo ff a c t o r st h a ti n f l u e n c i n gt h ee n e r g ye f f i c i e n c yf o rh y d r a t ed i s s o c i a t i o ni st e m p e r a t u r eo ft h ei n j e c t e dh o tb r i n e ,t h es a t u r a t i o n ,i n i t i a lt e m p e r a t u r e ,t h ei n j e c t i o nt i m ea n dr a t eo ft h eh o tb r i n e .T h eh i g h e rh y d r a t e s a t u r a t i o na n di n i t i a lt e m p e r a t u r e ,t h eh i g h e re n e r g ye f f i c i e n c y ,w h i l et h eh i g h e rt h eh o t b r i n et e m p e r a t u r e ,t h el o w e rt h e e n e r g ye f f i c i e n c y .T h et e m p e r a t u r e ,i n j e c t i o nt i m ea n di n j e c t i o nr a t eo fh o t b r i n ec anb eo p t i m i z e di no r d e rt oa c h i e v eah i I g h t h e r m a le f f i c i e n c yi nt h ep r o c e s so ft h eh y d r a t ed i s s o c i a t i o n .U n d e rt h ee x p e r i m e n t a lc o n d i t i o n s ,am a x i m u me n e r g ye f f i c i e n c yo f 6 .7 4i sa c h i e v e dd u r i n gt h ed i s s o c i a t i o no fh y d r a t ef o r m e da t5 ℃w i t has a t u r a t i o no f 4 8 %.w h e nh o t b r i n eo f 4 0 ℃i si n j e c t e d a tar a t eo f1 2m L /m i nf o r3 5 0m i n . K e yw o r d s n a t u r a lg a sh y d r a t e ;h o tb r i n ei n j e c t i o n ;h y d r a t ed i s s o c i a t i o n ;e n e r g ye f f i c i e n c y 天然气水合物的开采方法主要有注热法、降压 法、注化学剂法、C O 置换法以及固体开采法等‘1 4 1 。 在注热开采方面,K a m a t h [ 5 | 、U l l e r i c h ‘6 I 、W o n m o ‘7 引、 U c h i d a [ 9 1 等试验研究了注热开采动态唐良广‘1 0 ] 、 收稿日期2 0 1 3 1 0 1 2 基金项目国家自然科学基金项目 5 1 2 7 4 2 2 7 ;5 1 3 7 4 2 3 2 ;中国科学院可再生能源与天然气水合物重点实验室项目 y 2 0 7 k 2 作者简介李淑霞 1 9 7 0 一 ,女,教授,博士,主要从事油藏数值模拟、天然气水合物开采机制研究。E m a i l 二l i s h u x i a u p c .e d u .c n 。 万方数据 1 0 0 中国石油大学学报 自然科学版 万丽华[ ] 、孙建业等[ 1 2 ] 、李淑霞等[ 1 3 ] 、郝永卯等[ 1 4 ] 也进行了沉积物中天然气水和物 N G H 的注热试 验研究.但上述研究对注热开采的能量效率及敏感 因素分析不够。因此,笔者通过正交设计试验[ 15 | , 对水合物藏注热水开采的地质因素和注热参数进行 敏感性分析,为不同水合物藏的注热开采提供指导。 1 试验设备及方法 1 .1 试验设备 采用自制的天然气水合物开采模拟试验系统, 其具体试验流程见图1 。其中一维试验管尺寸为 驴8 0c m 8 0c m .工作压力2 0M P a 。沿试验管轴向 设计有1 1 套梯度取压、测温装置和电极系。填砂管 的注入端配有供液和供气设备,出口端设有回压装 置和气、液分离装置,利用数据采集处理软件自动采 集不同时刻的温度、压力、产气、产水等参数。 电子 天平 预 热 容 器 恒 温 箱 厂 压力传感器 1 放空 舶 氮 气 瓶 i 主主兰三妄妄四i 回压阀 温度传感器L l ;; 过滤器数据采集 图1天然气水合物开采模拟试验装置示意图 F i g .1 S k e t c hm a po fe x p e r i m e n t a ls y s t e m s i m u l a t i n gN G He x p l o i t a t i o n 1 .2 试验材料及步骤 试验中注入水为配置的质量分数为2 .0 %的盐 水,其中蒸馏水自制,N a C l 纯度≥9 9 .5 %;试验用 C H 。气纯度为9 9 .9 %。填砂管填人粒径为3 0 0 ~ 4 5 0 m 的石英砂,填砂后测得孑L 隙度为3 3 .4 %,渗 透率为1 .2t x m 2 。 试验操作步骤如下 1 N G H 等容生成。N G H 生成具体步骤参见 文献『1 6 ] 。当一次水合物的生成不能达到所需要 的饱和度时,需要进行二次注水或注气,然后继续降 温合成,直至达到设计的N G H 饱和度。 2 注热水开采。①加热预热罐中的盐水,使 其达到预定的注热温度,之后打开管线保温带,其设 定温度与预热罐的温度一致;②调节出口端回压阀 的压力,使其与试验管中的压力达到一致,之后打开 出口阀门③打开进水阀门,用平流泵按照预定的注 热速度向填砂管内注入热盐水,记录注热过程中的 产气产水量及压力、温度的变化;④当达到设定的注 热时间时,停止注热水,等待气体不再产出后,逐步 降低回压阀的压力到大气压,使填砂管中的剩余气 体产出。 2 注热水开采敏感因素试验 2 .1 正交设计试验结果 目前已发现的天然气水合物藏水合物饱和度差 别较大,从5 %到9 5 %[ 2 4 - 1 7 ] ,但渗透率随水合物饱 和度呈指数递减,在试验条件下,当水合物饱和度大 于5 0 %以后,渗透率已小于1 1 0 。I x m 2 ,使得注热 开采在技术上存在一定的难度,所以本组试验水合 物饱和度选取1 2 %~4 8 %实际水合物藏的初始温 度变化范围也比较大,从低于0 ℃到1 0 ℃以上,本 次试验选取一1 ~5o C 同时,根据实验室前期做过的 注热分解试验,确定注热量的原则是尽量使水合物 最大程度的分解且尽量减少注入的热量。根据以前 试验参数的合理范围,选取注热水的温度为4 0 ~ 1 0 0o C ,注热水的时间为1 5 0 ~4 5 0r a i n ,注热水速度 为9 ~1 8m L /m i n 。 采用五因素四水平正交设计试验,分别在不同 地质因素和注热参数下进行试验,研究水合物饱和 度A 、初始温度曰、注热水温度C 、注热水时间D 、注 热水速度E5 个参数对能量效率的影响,各因素取 值见表1 。 表1 正交试验因素水平表 T a b l e1P a r a m e t e r si no r t h o g o n a le x p e r i m e n t 水平A /% B P Cc /℃D /m i n E / m L - r a i n 。 1 214 0t 5 0 2 416 02 5 0 3 638 03 5 0 4 851 0 04 5 0 利用正交设计方案进行注热水开采试验。分析 各次试验的能量效率,并进行参数敏感性分析及最 优方案组合,结果见表2 。 表2 中均值表示每个因素每个水平下的平均能 量效率,极差表示每个因素各个水平平均能量效率 的最大值与最小值的差。可以看出,对注热水开采 能量效率影响由大到小的因素依次为注热水温度、 水合物饱和度、初始温度、注热水时间、注热水速度。 能量效率最大的组合为水合物饱和度4 8 %、初始温 度5 ℃、注热水温度4 0 ℃、注热水时间3 5 0r a i n 、注 热水速度1 2m L /m i n ,即为本组试验中的第1 6 次试 万方数据 第3 8 卷第2 期李淑霞。等天然气水合物藏注热水开采敏感因素试验研究 1 0 1 - 验,能量效率达到了6 .7 4 。 表2 正交设计试验的基本数据及敏感因素分析 T a b l e2B a s i cd a t aa n ds e n s i t i v i t yf a c t o r s a n a l y s i si no r t h o g o n a le x p e r i m e n t 2 .2 敏感因素 各个因素对能量效率影响的趋势见图2 。由图 2 可以看出 1 水合物饱和度越大,注热水开采的能量效率 越大。这是因为当其他条件相同时,水合物饱和度越 大,水合物分解吸热越多,热损失相对越小,注入同样 的热量时产出的气体越多,所以能量效率越大。 J I l 2 3 A 4 B l 占2B 3 8 4 C IC 2C 3 7 4D I D 2 D 3 口4 F l E 2 E 3E 4 因素水平 图2 能量效率随各因素变化趋势 F i g .2 V a r i a t i o no fe n e r g ye f f i c i e n c yw i t hv a r i o u sf a c t o r s 2 随着初始温度的增加,注热水开采的能量 效率增大。原因是当其他条件相同时,初始温度越 高,水合物分解所需的温差驱动力越小,从而用于水 合物分解的热量就相对越多,能量效率越大。 3 随着注热水温度的增加,能量效率逐渐减 小。因为当其他条件一致时,注入的热水温度越高, 在注热水过程中水合物已分解部位的温度升高得越 高,与周围环境的温差越大.向周围散失的热量越 多,能量效率越低。 4 对于注热水时间而言,在注热水前期 3 5 0 r a i n 以前 ,能量效率变化不大,在注热水后期 3 5 0 r a i n 以后 能量效率逐渐变小。这是由于注热前期 水合物没有完全分解,当其他条件不变时,水合物分 解产气量与注热时间成正比,而注入的热量也与注 热时间成正比,所以注热前期能量效率变化较小。 在注热水后期,水合物已经分解完毕,产气速率变得 很小,继续注热则能量效率逐渐变小。 5 随着注水速度的增加,能量效率先增加后 减小。这是因为随着注热水速度的提高,热前缘推 进速度加快,水合物能有效分解,从而能量效率增 加但注水速度太大使得温度升高加快,反应釜内外 温差加大,热损失加大,能量效率降低。 2 .3 注热水开采 综上分析可以看出考虑水合物藏的地质参数 时,水合物饱和度越大,初始温度越高,则能量效率 越大,但水合物饱和度的影响远大于初始温度的影 响。对实际水合物的开采,当水合物饱和度高到一 定程度后,水合物藏渗透率会急剧降低,从而会增大 注水压力,使注热水开采的可行性变差。建议对高 饱和度水合物藏首先进行相应的地层处理,然后再 进行注热水开采。 为了获得更大的能量效率。需要合理优化各种 注热参数,本组试验条件下,注热水温度对能量效率 影响最大,注热水温度越小,能量效率越大。但对水 合物藏的实际开采.当注热水温度比较低时会影响 水合物的分解速率,所以注入热水的温度不能过低 其次要选择合适的注热时间,使得注入的热量刚好 使水合物全部分解,避免注热量过多,引起能量效率 减少;当注热水速度过大或者过小时.都会影响水合 物分解产气量.所以需要合理优化注热水速度,使总 的能量效率达到最大。 3结论 1 对水合物藏注热水开采能量效率影响由大 到小的因素依次为注热水温度、水合物饱和度、初始 温度、注热水时间、注热水速度。 2 本试验条件下能量效率最大的参数组合为 水合物饱和度4 8 %、初始温度5 ℃、注热水温度4 0 ℃、注热水时间3 5 0r a i n 、注热水速度1 2r a L /m i n ,能 万方数据 1 0 2 中国石油大学学报 自然科学版 2 0 1 4 年4 月 量效率达到了6 .7 4 。 3 为了获得更高的能量效率,需要合理优化 注热参数在确保水合物分解速度前提下,适当降低 注热水温度;优化注热水时间,使得注入的热量刚好 川 使水合物完全分解;优化注热水速度,使产气量达到 最大。 参考文献 [ 1 ] 龚建明.M a l l i k 天然气水合物国际研讨会综述[ J ] .海 洋地质动态,2 0 0 4 ,2 0 2 2 8 3 1 . G O N GJ i a n m i n g .S u m m a r yo ft h ei n t e r n a t i o n a ls y m p o s i U Bo fM a l l i kn a t u r a lg a sh y d r a t e [ J ] .M a r i n eG e o l o g y L e t t e r s ,2 0 0 4 ,2 0 2 2 8 - 3 1 . [ 2 ] 李淑霞,陈月明,杜庆军.天然气水合物开采方法及数 值模拟研究评述『J ] .中国石油大学学报自然科学 版,2 0 0 6 ,3 0 3 1 4 6 1 5 0 . L IS h u x i a ,C H E NY u e m i n g ,D UQ i n g - j u n .C o m m e n t a r y o fp r o d u c t i o nm e t h o da n dn u m e r i c a ls i m u l a t i o no fn a t u r a l g a sh y d r a t e s [ J ] .J o u r n a lo fC h i n aU n i v e r s i t yo fP e t r o l e u m E d i t i o no fN a t u r a lS c i e n c e ,2 0 0 6 ,3 0 3 1 4 6 1 5 0 . [ 3 ] T A R U NG ,M O R I D I SGJ ,H O L D I T C HSA .A n a l y s i so f r e s e r v o i rp e r f o r m a n c eo ft h eM e s s o y a k h ag a sh y d r a t er e s e r v o i r [ R ] .S P E11 4 3 7 5 ,2 0 0 8 . [ 4 ]郭平,刘士鑫,杜建芬.天然气水合物气藏开发[ M ] . 北京石油工业出版社,2 0 0 6 . [ 5 ] K A M A T HVA ,G O D B O L ESP .E v a l u a t i o no fh o tb r i n e s t i m u l a t i o nt e c h n i q u ef o rg a sp r o d u c t i o nf r o mn a t u r a lg a s h y d r a t e s 『R ] .S P E1 3 5 9 4 .1 9 8 7 . [ 6 ] U L L E R I C HJW ,S E L I MMS ,S L O A NED .T h e o r ya n d m e a s u r e m e n to f h y d r a t ed i s s o c i a t i o n [ J ] .A I C h EJ , 1 9 8 7 ,3 3 5 7 4 7 - 7 5 2 . [ 7 ] W O N M OS ,H O S E O BL ,S U N J O O NK .E x p e r i m e n t a li n v e s t i g a t i o no fp r o d u c t i o nb e h a v i o r so fm e t h a n eh y d r a t es a t u r a t e di np o r o u sr o c k [ J ] .E n e r g yS o u r c e s ,2 0 0 3 ,2 5 8 4 5 .8 5 6 . [ 8 ] W O N M OS ,H O S E O BL ,H O J O O NY .A ne x p e r i m e n t a l s t u d yf o rh y d r a t ed i s s o c i a t i o np h e n o m e n aa n dg a sf l o w i n g a n a l y s i sb ye l e c t r i ch e a t i n gm e t h o di np o r o u sr o c k s [ J ] . K o r e a nC h e mE n gR e s ,2 0 0 4 ,4 2 1 1 1 5 1 2 0 . [ 9 ] U C H I D AT ,T A K E Y AS ,C H U V I L I NEM .D e c o m p o s i t i o no fm e t h a n eh y d r a t e si ns a n d ,s a n d s t o n e ,c l a y s ,a n d g l a s sb e a d s [ J ] .G e o p h y sR e s ,2 0 0 4 ,1 0 9 B 5 B 0 5 2 0 6 . [ 1 0 ] 唐良广,肖睿,李刚,等.热力法开采天然气水合物的 模拟实验研究[ J ] .过程工程学报,2 0 0 6 ,6 4 5 4 8 5 5 3 . T A N GL i a n g g u a n g ,X I A OR u i ,L IG a n g ,e ta 1 .E x p e r - [ 1 2 ] [ 1 3 ] [ 1 4 ] [ 1 5 ] [ 1 6 ] [ 1 7 ] i m e n t a li n v e s t i g a t i o no fp r o d u c t i o nb e h a v i o ro fg a sh y d r a t eu n d e rt h e r m a ls t i m u l a t i o n 『J ] .T h eC h i n e s eJ o u r - n a l o fP r o c e s sE n g i n e e r i n g ,2 0 0 6 ,6 4 5 4 8 - 5 5 3 . 万丽华,李小森,李刚,等.热盐水分解甲烷水合物实 验研究[ J ] .现代化工,2 0 0 8 ,2 8 7 4 7 5 0 . W A NL i h u a ,L IX i a o - s e n ,L IG a n g ,e ta 1 .A ne x p e f i m e n t a ls t u d yo nd i s s o c i a t i o nb e h a v i o ro fm e t h a n eh y d r a t e u n d e rh o tb r i n es t i m u l a t i o n 『J ] .M o d e r nC h e m i s t r yI n d u s t r y ,2 0 0 8 ,2 8 7 4 7 5 0 . 孙建业,业渝光,刘昌岭,等.沉积物中天然气水合物 合成及开采模拟实验研究[ J ] .中国海洋大学学报 自然科学版,2 0 0 9 ,3 9 6 1 2 8 9 .1 2 9 4 . S U NJ i a n y e ,Y EY u g u a n g ,L I UC h a n g l i n g ,e ta 1 . S i m u l a t i o ne x p e r i m e n to fg a sh y d r a t ef o r m a t i o na n de x p l o i t a t i o ni ns e d i m e n t s 『J ] .P e r i o d i c a lo fO c e a nU n i v e r s i t yo fC h i n a ,2 0 0 9 ,3 9 6 1 2 8 9 - 1 2 9 4 . 李淑霞.郝永卯,陈月明.多孔介质中天然气水合物 注热盐水分解实验研究[ J ] .太原理工大学学报, 2 0 1 0 ,4 1 5 6 8 0 - 6 8 3 . L IS h u x i a ,H A OY o n g m a o ,C H E NY u e m i n g .E x p e r i m e n t a ls t u d yo fh o t b r i n es t i m u l a t i o nf o rd i s s o c i a t i o no f N G Hi np o r o u sm e d i u m [ J ] .J o u r n a lo fT a i y u a nU n i v e r s i t yo fT e c h n o l o g y ,2 0 1 0 ,4 1 5 6 8 0 - 6 8 3 . 郝永卯,陈月明,李淑霞.天然气水合物注热开采实 验研究[ J ] .中国石油大学学报自然科学版,2 0 0 7 , 3 1 4 6 0 - 6 3 . H A OY o n g m a o ,C H E NY u e m i n g ,L IS h u x i a .E x p e f i m e r a ls t u d yo np r o d u c t i o no fn a t u r a lg a sh y d r a t eb yt h e r m a ls t i m u l a t i o n 『J ] .J o u r n a lo fC h i n aU n i v e r s i t yo fP e t r o l e u m ,2 0 0 7 ,3 1 4 6 0 6 3 . 方开泰,马长兴.正交与均匀试验设计[ M ] .北京 科学出版社,2 0 0 1 . 李淑霞,陈月明,王瑞和,等.初始压力对多孔介质 中气体水合物生成的影响lJ j .实验力学,2 0 0 9 ,2 4 4 3 1 3 .3 1 9 . HS h u x i a ,C H E NY u e m i n g ,W A N GR u i h e .O nt h e i n f l u e n c eo fi n i t i a lp r e s s u r eo nn a t u r a lg a sh y d r a t ef o r m a t i o ni nap o r o u sm e d i u m [ J ] .J o u r n a lo fE x p e r i m e n t a l M e c h a n i c s ,2 0 0 9 ,2 4 4 3 1 3 3 1 9 . 卢振权,祝有海,张永勤,等.青海祁连山冻土区天然 气水合物资源量的估算方法以钻探区为例[ J ] . 地质通报,2 0 1 0 ,2 9 9 1 3 1 0 1 3 1 8 . L UZ h e n q u a n ,Z H UY o u h a i ,Z H A N GY o n g q i n ,e t a 1 .E s t i m a t i o nm e t h o do fg a sh y d r a t er e s o u r c ei nt h eQ i l i a nM o u n t a i np e r m a f r o s ta r e a .Q i n g h a i .C h i n a ac a s e o ft h ed r i l l i n ga r e a [ J ] .G e o l o g i c a lB u l l e t i no fC h i n a , 2 0 1 0 ,2 9 9 1 3 1 0 1 3 1 8 . 编辑李志芬 万方数据 天然气水合物藏注热水开采敏感因素试验研究天然气水合物藏注热水开采敏感因素试验研究 作者李淑霞, 李杰, 徐新华, 李小森, LI Shu-xia, LI Jie, XU Xin-hua, LI Xiao-sen 作者单位李淑霞,LI Shu-xia中国石油大学石油工程学院,山东青岛266580;中国科学院可再生能源与天然气水合物 重点实验室,广东广州510640, 李杰,徐新华,LI Jie,XU Xin-hua中国石油大学石油工程学院,山东青岛 ,266580, 李小森,LI Xiao-sen中国科学院可再生能源与天然气水合物重点实验室,广东广州,510640 刊名 中国石油大学学报(自然科学版) 英文刊名Journal of China University of PetroleumEdition of Natural Science 年,卷期2014,382 本文链接
展开阅读全文