资源描述:
第3 6 卷第3 期 爆破 V 0 1 .3 6N o .3 2 0 1 9 年9 月B L A S T I N G S e p .2 0 1 9 d o i 1 0 .3 9 6 3 /j .i s s n .1 0 0 1 4 8 7 X .2 0 1 9 .0 3 .0 1 9 爆破荷载下在役储油罐动态响应实验研究 杨德明““,熊祖钊a , b ,吴a , b 胡从骄“6 ,夏世林““ 武汉科技大学a .理学院工程力学系;b .中铁港航一武汉科技大学爆破技术研究中心,武汉4 3 0 0 6 5 摘要大型原油储罐的建设是国家经济与安全的重要保障,储罐地基的建设往往涉及到岩体爆破开挖作 业,这势必对在役储罐造成影响。为确保储罐建设工程安全实施,设计了一套缩尺模型试验方案,对储油罐 罐体与基础的动力响应进行了监测。分析了储油罐液高分别为1 /4 、1 /2 、3 /4 和满罐情况下,储油罐近区岩 体爆破开挖扰动下罐壁和底板的动态应变规律。结果表明罐体应变随罐高增加有减小趋势,最大应变出现 在罐底1 /5 位置,且罐壁与罐底连接部位是薄弱环节;随着储液高度增加,罐体和罐底应变呈增大趋势,且超 过1 /2 时增幅急剧变大。建议工程中对罐壁底部1 /5 以及与罐底结合部位进行重点监测,并且爆破时储液 量不宜超过罐体1 /2 。实验结果可为相关工程提供参考数据。 关键词储油罐;爆破振动;模型试验;动力响应 中图分类号T D 2 3 5文献标识码A文章编号1 0 0 1 4 8 7X 2 0 1 9 0 3 0 1 2 4 0 5 E x p e r i m e n t a lS t u d yo nt h eS t r u c t u r a lD y n a m i cR e s p o n s e o fO i lS t o r a g eT a n ku n d e rB l a s t i n gL o a d Y A N GD e m i n 9 8 一,X I O N GZ u z h a 0 8 一,删L i a n 9 8 一,H UC o n g - f i a 0 8 一,X I AS h i l i n 9 8 ,6 a .D e p a r t m e n to fM e c h a n i c s ,C o l l e g eo fS c i e n c e .b .C h i n aR a i l w a yP o r ta n dC h a n n e l E n g i n e e r i n gG r o u pC oL t d R e s e a r c hC e n t e ro fB l a s t i n gT e c h n o l o g y . W u h a n U n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y ,W u h a n4 3 0 0 6 5 ,C h i n a A b s t r a c t T h ec o n s t r u c t i o no fl a r g ec r u d eo i ls t o r a g et a n k si sa ni m p o r t a n tg u a r a n t e ef o rt h en a t i o n a le c o n o m y a n ds a f e t y .A n dt h ec o n s t r u c t i o no ft a n kf o u n d a t i o no f t e ni n v o l v e st h eb l a s t i n ge x c a v a t i o no fr o c km a s s ,w h i c hw i l li n e v i t a b l ya f f e c tt h ei n - s e r v i c es t o r a g et a n k s .I no r d e rt oe n s u r et h ei m p l e m e n t a t i o no ft h et a n kc o n s t r u c t i o np r o j e c ts a f e t y ,as c a l em o d e lt e s ts c h e m ew a sd e s i g n e dt om o n i t o rt h ed y n a m i cr e s p o n s eo ft h et a n kb o d ya n df o u n d a t i o n .T h ed y n a m i cs t r a i nl a wo ft a n kw a l la n df l o o ru n d e rt h ed i s t u r b a n c eo fb l a s t i n ge x c a v a t i o no fr o c km a s si nt h ev i c i n i t yo fo i l s t o r a g et a n kw a sa n a l y z e dw h e nt h et a n kl i q u i dh e i g h tw a s1 /4 ,1 /2 ,3 /4 a n df u l lt a n kr e s p e c t i v e l y .T h er e s u l t s s h o w e dt h a tt h es t r a i no ft a n kb o d yd e c r e a s e dw i t ht h ei n c r e a s eo ft a n kh e i g h tw i t ham a x i m u ms t r a i na tt h ep o s i t i o n o f1 /5o ft a n kh e i g h tf r o mb o t t o m ,a n dt h ec o n n e c t i o nb e t w e e nt a n kw a l la n dt a n kb o t t o mw a saw e a kl i n k .W i t ht h e i n c r e a s eo fl i q u i ds t o r a g eh e i g h t ,t h es t r a i no ft a n kb o d ya n dt a n kb o t t o mi n c r e a s e d ,a n dt h eg r o w i n gr a t ew e n tu p s h a r p l yw h e nt h el i q u i ds t o r a g eh e i g h te x c e e d e d1 /2 .I tw a ss u g g e s t e dt h a tm o n i t o r i n gs h o u l db et a k e na t1 /5l i q u i d s t o r a g eh e i g h tf r o mb o t t o mo ft h et a n kw a l la n dt h ej o i n tw i t ht h et a n kb o t t o m .T h el i q u i ds t o r a g es h o u l dn o te x c e e d 1 /2o ft h et a n kc a p a c i t yd u r i n gb l a s t i n g .T h ee x p e r i m e n t a lr e s u l t sc a np r o v i d er e f e r e n c ed a t af o rr e l a t e dp r o j e c t s . K e yw o r d s o i ls t o r a g et a n k ;b l a s t i n gv i b r a t i o n ;m o d e lt e s t ;d y n a m i cr e s p o n s e 万方数据 第3 6 卷第3 期杨德明,熊祖钊,吴亮,等爆破荷载下在役储油罐动态响应实验研究 近些年来,随着经济的快速发展,我国对石油需 求量不断增大。目前,储油罐逐渐向大型化发展,我 国也相继建成了不少直径6 0 ~1 5 0m 级别的大型储 油罐⋯。由于油罐规模庞大,且存放有易燃易爆有 害液体,一旦发生破坏,后果不堪设想。为确保罐区 基础爆破施工期间储油罐的安全运行,亟需对爆破 振动下储油罐结构的动应力应变、响应规律等问题 进行探讨[ 2 ] 。 在储油基地的建设过程中,由于地形的原因或 者出于安全的考虑,部分储油罐依山而建,在储油基 地的扩建过程中就涉及基础爆破开挖,而爆破振动 有可能影响既有储油罐结构,从而为后续运行带来 安全隐患。众所周知,爆破地震波具有幅值大、衰减 快、频率高、持续时问短的特点。根据已有的理论, 在基岩爆破开挖过程中邻近的储油罐在承受爆破荷 载时罐壁可能产生“象足”变形,这是因为罐壁轴向 压应力过大,罐壁局部失稳产生屈曲破坏,而且随着 储油罐液高的增加,罐壁质点振速都会增加[ 3 ] 。由 于罐内液体对液面以下罐壁质点振动有一定的缓冲 作用,导致罐壁质点振速的分布情况较为复杂,并且 罐壁质点最大质点振速与危险点的位置并不对应。 王开志等根据现场环境,设计了某部储油罐基础的 爆破开挖方案,优化了布孔方式和炮孔参数,并成功 实施了爆破,为同类的爆破工程提供了参考依 据㈠J 。TYw u 和GRL i u 等就罐底环形板与混凝 土环墙的设计方法进行了研究J 。许峰等运用A N S Y A /L S D Y N A ,研究了爆破荷载作用下大型储油罐 以及不同管道结构的动力响应规律J 。陈洋等以 实际工程为背景,运用A N S Y A /L S D Y N A 研究了爆 破振动下大型储油罐的响应规律[ 7 ] 。张美玲采用 A D I N A 有限元数值模拟方法,研究了储罐的静力作 用、动力特性以及储罐的地震动响应;并对比分析了 储液晃动、储液密度和深度变化对罐体带来的影响, 以及储液液面晃动波高变化等[ 8 ] 。路胜卓利用L S D Y N A 非线性有限元计算软件,通过缩比模型的破 坏模拟试验对数值模拟结果的可靠性进行了验 收稿日期2 0 1 9 0 4 2 2 作者简介杨德明 1 9 9 4 一 ,男,湖北孝感人,武汉科技大学在读硕 士,主要研究方向为爆炸理论及其应用, E m a i l 2 9 4 3 4 11 6 2 6 q q .c o m 。 通讯作者熊祖钊 1 9 6 3 一 ,男,博士、武汉科技大学教授,主要研究 方向为控制爆破、工程爆破, E .m a i l 9 3 7 7 0 2 9 7 8 q q . c o r n o 基金项目国家自然科学基金 编号5 1 7 7 9 1 9 3 爆破荷载作用下水泥 灌浆岩体的破坏机制;国家自然科学基金 编号 51 4 7 9 1 4 7 高应力区岩体的爆破破裂机制及其数值分析方 法;武汉科技大学2 0 1 7 年度研究生创新创业基金 编号 J C X 2 0 1 7 0 2 8 5 5 约束管道模型的简化及计算的研究 证[ 9 ] ,结果表明在爆炸作用下,拱顶储油罐呈现弹 塑性动力响应特性。 上述的研究都是针对储油罐进行的数值模拟, 较少涉及到模型试验,目前针对储油罐近区爆破开 挖对罐体的影响的实验研究相对较少。因此,设计 了一套模拟储油罐近区爆破实验方案,以此探究爆 破开挖对已建成的储油罐的影响,并为储油罐近区 爆破开挖提供了一套参考实验方案。 1实验部分 1 .1 实验目的 本实验通过在扇形水泥地基的炮孔里面装药起 爆从而模拟出爆破开挖扰动对附近既有储油罐的影 响,然后通过对储油罐上的仪器采集到的应变数据 进行分析,揭示爆破开挖对附近既有的大型储油罐 的影响规律,从而为既有储油罐附近的基础爆破开 挖施工提供理论依据,并为后续爆破开挖对附近既 有储油罐影响的研究提供实验方案。试验模型见 图1 。 形水泥地基 图1 实验模型 F i g .1E x p e r i m e n t a lm o d e l 1 .2 实验原理 通过对扇形水泥地基的炮孔装药填塞并起爆来 模拟储油罐近区的爆破开挖,在炸药起爆前点击动 态信号采集软件上的采集按钮就可以采集到爆破过 程中储油罐模型的应变变化规律,同时炸药起爆时, 动态应变仪会自动采集爆破时储油罐和扇形水泥的 振动数据。分析仪器记录下来的数据,得出结论。 背姆血 背爆而伽而迎爆而背爆而 .厂八一 ‘迫岁。j y 1- 底1 [ 6 】 【5J 伽I 仃i [ a l 【3 J [ 2 ] 【l J i I - 一 竺I _ 一 I _ 一 ● t _ “ - IL , 迎搬而 图2 储油罐应变片位置示意图 单位c m F i g .2D i a g r a mo ft a n ks t r a i ng a u g ep o s i t i o n u n i t c m 万方数据 1 2 6爆破2 0 1 9 年9 月 1 .3 实验方案 实验模型如图1 所示,圆形为已打好炮孔的水 泥基座,共分为十个扇形部分,每个部分之间都有隔 振沟,同时又各自与中心基座相连接,从而保证在不 破坏邻近扇形的前提下,将爆破产生的应力波传到、 至实验模型。实验时,在扇形水泥基座模型的炮孔 上装药,填塞和起爆,炸药起爆时会产生应力波,而 应力波会对储油罐模型造成一定的影响。储油罐模 型上贴有2 1 组应变片,炸药起爆时可以测到爆破时 储油罐的应变数据。储油罐应变片位置如图2 所 示,通过储油罐上仪器记录到的数据来分析爆破对 储油罐模型的影响。 1 .4 实验的优缺点 实验的优点 1 本实验的1 0 个扇形地基之间被隔断,众所 周知,当一个压缩波传播到自由面时,会反射一个往 回传播、大小相等的反射拉伸波。所以本实验保证 各个地基之间的独立性。 2 实验考虑了不同的液高和爆源不同的情 况,可以得到较为全面的结论。 实验的缺点 1 当一个扇形水泥地基做完实验准备进行下 一个水泥地基的实验时,都需要人工的转动储油罐 的方向,以保证每一次的迎爆面和背爆面都要一样, 这样才能保证实验数据的可比性。 2 由于实验模型的局限性,为保证罐底的安 全,每一次切换不同的水高时,都需要人工加水。 1 .5 试验步骤 试验主要分为3 部分 1 模型制作 在较为平整且土质较坚硬的地方浇筑圆形的水 泥地基,并用木板隔成十等分,待水泥初步凝固之 后,移除木板。待水泥地基达到实验所需要的强度 之后,在扇形水泥地基上距离中心分别为1 4 0c m 、 1 1 0c m 、8 0c m 、5 0c m 的地方打好炮孔,所需炮孔直 径为4c m ,孔深为1 2c m 。 2 测试系统布置与连接 在储油罐模型上预定位置位置贴好应变片,然 后将储油罐模型放置到水泥地基的中心位置,连接 好仪器,动态信号采集器需要接地。 3 荷载施加与数据采集 连接电脑,打开动态信号采集器以及相对应的 数据采集软件,将数据调试好。炮孑L 装药长度为 8c m ,填塞长度为4c m ,连接好起爆线路,准备起 爆,听到爆破倒计时口令马上点击数据采集。爆破 完成后停止数据采集,采集到爆破应变数据。调试 设备,准备下一轮的实验。 2 实验数据的分析 储油罐中心距爆源8 0c m ,液高为3 /4 时各测 点的应变值如表1 所示。 表18 0c m 爆距3 /4 液高各测点应变极值 T a b l e1M a x i m u ms t r a i nv a l u ea te a c hm e a s u r i n gp o i n tw i t h8 0c m e x p l o s i o nd i s t a n c ea n d3 /4l i q u i dh e i g h t 图3 分别为爆距8 0c m ,液高3 /4 工况下各测点 的Y 轴、名轴应变值与应变片布置高度的关系。结 果表明Y 轴和z 轴应变值随着高度的增加都呈现 出减小的趋势。迎爆面的应变值的变化十分明显, 万方数据 第3 6 卷第3 期杨德明,熊祖钊,吴亮,等爆破荷载下在役储油罐动态响应实验研究 1 2 7 并且在罐高1 /5 以下位置均出现较大的应变,这表 明罐高1 /5 以下为储油罐的危险区域。而环向应力 的最大值也基本出现在罐高1 /5 的高度以下,与轴 5 0 O , l 4 5 京4 ‘’ - N3 5 3 0 童2 5 重2 0 乳 ‘ 向应力的变化一致,并且与储油罐“象足”破坏理论 基本一致。 ∞ i 、 制 熠 厘 遗 .蔓 _ 羼 姑 簧 蛞 测点高度测点高度 a y 方向应变 I z 方向应变 a y - d i r e c t i o ns t r a i n h Jz - d i r e c t i o ns t r a i n 图38 0c m 爆距3 /4 液高时应变片布置高度与叭z 方向应变值的关系 F i g .3 T h er e l a t i o nb e t w e e ns t r a i ng a u g ea r r a n g e m e n th e i g h ta n ds t r a i nv a l u ei nYa n d 彳d i r e c t i o n a t3 /4l i q u i dh e i g h tw i t h8 0e me x p l o s i o nd i s t a n c e 图4 为爆距1 1 0c m 工况下,不同液高下测点应 变值的变化曲线图,由于测点较多,通过分析数据, 札 i 制 1 兰| 哥 H , 嵯 嗤 引 储液高度 ㈨也爆面 a 、E x p l o s i o n8 u “ l u e 储油罐的危险位置出现在1 /5 罐高以内,故选取迎 爆面和背爆面的1 号和2 号测点数据进行展示。 沁 i 制 崔 嚣 H , _ 嵯 嗤 乏 储液高度 f b ] 背爆面 f bJBa 1 .c ke x p l o s i o ns u r f a c e 图41 1 0c m 爆距工况下迎爆面、背爆面1 、2 测点的应变与储液高度的关系 F i g .4 T h er e l a t i o nb e t w e e ns t r a i na n dl i q u i ds t o r a g eh e i g h ta tm e a s u r i n gp o i n t s1a n d2o ft h e b l a s t i n gs u r f a c ea n db a c k b u r s ts u r f a c ea tt h ee x p l o s i o nd i s t a n c eo f11 0c m 分析图4 数据结果表明随着储液高度的增加, 每个测点的应变都逐渐增大,且从1 /2 到3 /4 液高 时应变急剧增加。其原因是因为随着储液增加,罐 壁液动压力也随之增加,罐体受到的剪力、弯矩增 大,此实验结果与文献[ 5 ] 、[ 6 ] 中的储罐罐壁应变 随液高增加而增大的结论基本吻合。 通过图5 数据分析表明底板径向应变从罐底边 缘沿径向逐渐减小,在罐底边缘应变为最大峰值。 环向应变从罐边缘沿径向有所增大但整体值较小。 其原因为罐底与罐壁交接处受到较大的力矩影响, 从而造成底板径向应变较大。这与曲晓建文献中的 现象基本一致0 | 。依据G B 5 0 3 4 1 油罐设计规 范1 1 j ,罐底边缘板厚度都比中腹板厚,这样的设计 要求增加了边缘板的强度,减小了其因边缘板径向 应力过大而破坏的几率。图6 为爆距1 4 0c m 工况 下储罐底部“底1 ”位置处应变与储液高度的关系, 从图可以看出随着储液高度增加其罐底边缘的应变 也随之增大,且储液从1 /2 到3 /4 时应变增幅较大, 满罐时罐底整体应变达到最大值。 万方数据 1 2 8爆破2 0 1 9 年9 月 0 土 坦 掣 崔 世 掣 8 0 6 0 4 0 0 底1底2底3 磷底测心编弓 图58 0c r n 爆距3 /4 液高底板各测点应变值 F i g .5 M a x i n m ms t r a i nv a l u ea te a c hm e a s u r i n gp o i n tw i t h 8 0c n le x p l o s i o nd i s t a n c ea n d3 /4l i q u i dh e i g h t c ‘ i 、 理 刽 崔 蘑 矗[ 世 4 1 3 l 2 0 l J 储液高度 图6 爆距为1 4 0a m 工况下底部测点应变与储液量的关系 F i g .6 T h er e l a t i o n s h i pb e t w e e nt h es t r a i no fb o t t o m m e a s u r i n gp o i n ta n dt h ev o l u m eo fl i q u i du n d e rt h ec o n d i t i o n o fe x p l o s i o nd i s t a n c eo f1 4 0e m 3结论 通过模型试验监测了不同储油罐液高条件下储 油罐近区岩体爆破开挖扰动下罐壁质点和底板质点 的应力变值变化过程,得到以下结论 1 在爆破施工中,减小储油罐储液量能够有 效减小爆破振动对储油罐罐壁与底板应变的影响。 爆破扰动下邻近储油罐应保持在半载状态以下。 2 在爆破地震波作用下,储罐底板与罐壁的 交界处为薄弱点,储液量由1 /2 增加至3 /4 时底板 应变急剧增加。本文所研究的储油罐危险点在储油 罐1 /5 高度处,爆破施工时应重点检测。 [ 1 ] [ 2 ] [ 2 ] 参考文献{ R e f e r e n c e s 贾庆山.储罐基础工程手册[ M ] .北京中国石化出版 社,2 0 0 2 2 1 - 2 3 5 . 路胜卓,王伟,张博一.大型浮顶式储油罐的爆炸破 坏机理实验[ J ] .爆炸与冲击,2 0 1l ,3 1 2 1 5 8 .1 6 4 . L US h e n g z h u o ,W A N GW e i ,Z H A N GB o .y i .E x p e r i m e n t a l r e s e a r c ho nd e s t r u c t i o nm e c h a n i s mo fl a r g e ..s c a l ef l o a t i n g .. r o o fo i lt a n ku n d e rb l a s tl o a d i n g [ J ] .E x p l o s i o na n dS h o c k W a v e s ,2 0 1 1 ,3 l 2 1 5 8 1 6 4 . i nC h i n e s e [ 3 ] 陈海阳.储罐地震动力响应分析及其地基处理[ D ] .东 营中国石油大学,2 0 1 2 . [ 3 ]C H E NH a i - y a n g .R e s e a r c ho fs e i s m i cr e s p o n s ea n ds o f t g r o u n dt r e a t m e n to fo i lt a n k s 【DJ .D o n g y i n g C h i n aU n i v e r s i t yo fP e t r o l e u m ,2 0 1 2 . i nC h i n e s e [ 4 ] 王开志,杨旭升,梁秋祥.某部储油罐基础开挖爆破设 谢J ] .工程爆破,2 0 1 6 ,2 2 4 7 7 8 1 . 【4JW A N GK a i z h i ,Y A N GX u - s h e n g ,L I A N GQ i u x i a n g . B l a s td e s i g no fo i ls t o r a g et a n kf o u n d a t i o ne x c a v a t i o n [ JJ . E n g i n e e r i n gB l a s t i n g ,2 0 1 6 ,2 2 4 7 7 - 8 1 . i nC h i n e s e 【5J TYW u ,GRL i u .C o m p a r i s o no fd e s i g nm e t h o d sf o ra t a n k b o t t o ma n n u l a rp l a t ea n dc o n e r e t er i n g w a l l [ JJ .I n t e r - n a t i o n a lJ o u r n a lo fP r e s s u r eV e s s e l sa n d P i p i n g , 2 0 0 0 9 5 1 1 .5 1 7 . [ 6 ] 许锋.爆破对罐区储油罐及地下管道结构的影响分 析[ D ] .武汉武汉科技大学,2 0 1 7 . [ 6 ]X UF e n g .A n a l y s i so fi n f l u e n c eo fs t o r a g et a n ka n du n d e r - g r o u n dp i p e l i n es t r u c t u r ei nt a n ka r e au n d e rb l a s t i n gl o a d i n dD ] .W u h a n W u h a nU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y ,2 0 1 7 . i nC h i n e s e [ 7 ] 陈洋,吴亮,许锋,等.爆破开挖振动下既有大 型储油罐的动力响应[ J ] .爆炸与冲击,2 0 1 8 ,3 8 6 1 3 9 4 .1 4 0 3 . 【7JC H E Nr a n g ,W UL i a n g ,X UF e n g ,e ta 1 .D y n a m i cr e s p o n s eo fe x i s t i n gl a r g eo i ls t o r a g et a n ku n d e rb l a s t i n ge x - c a v a t i o nv i b r a t i o n [ JJ .E x p l o s i o na n dS h o c kW a v e s ,2 0 1 8 , 3 8 6 1 3 9 4 1 4 0 3 . i nC h i n e s e [ 8 ] 张美玲.大型非锚固储油罐抗震性能研究[ D ] .成都 西南交通大学,2 0 1 8 . [ 8 ] Z H A N GM e i l i n g .E a r t h q u a k e r e s i s t a n c ea n a l y s i so fl a r g e u n a n c h o r e do i l s t o r a g et a n k 【DJ .C h e n g d u S o u t h w e s t J i a o t o n gU n i v e r s i t y ,2 0 1 8 . i nC h i n e s e [ 9 ] 路胜卓,陈卫东,王伟,等.拱顶储油罐爆炸作用下的动 力响应数值模拟[ J ] .油气储运,2 0 1 8 ,3 7 6 6 4 4 6 5 0 . 【9J L US h e n g z h u o ,C H E NW e i - d o n g ,W A N GW e i ,e ta 1 .N u - m e r i c a ls i m u l a t i o no nt h ed y n a m i cr e s p o n s eo ff i x e d r o o f o i ls t o r a g et a n ku n d e rt h ee f f e c to fb l a s tl o a d i n g [ J ] .O i l & G a sS t o r a g ea n dT r a n s p o r t a t i o n ,2 0 1 8 ,3 7 6 6 4 4 6 5 0 . i nC h i n e s e [ 1 0 ] 曲晓建,李玉坤,段冠,等.大型储罐罐底及边缘板应 力分析方法的对眦J ] .油气储运,2 0 1 1 1 2 9 1 9 - 9 2 2 . 【1 0JQ UX i a o - j i a n ,L IY u k u n ,D U A NG u a n ,e ta 1 .C o n t r a s to f s t r e s sa n a l y s i sm e t h o d sf o rt h eb o t t o ma n da n n u l a rp l a t e o fl a r g e s c a l et a n k [ JJ .O i l G a sS t o r a g ea n dT r a n s p o r t a t i o n ,2 0 1 1 1 2 9 1 9 .9 2 2 . i nC h i n e s e [ 11 ] G B 5 0 3 4 1 - - 2 0 1 4 立式圆筒形钢制焊接油罐设计规范 [ s ] .北京中国标准出版社,2 0 1 4 . [ 1 1 ] G B 5 0 3 4 1 - - 2 0 1 4C o d ef o rd e s i g no fv e r t i c a lc y l i n d r i c a l w e l d e ds t e e lo i lt a n k s [ s ] .B e i j i n g S t a n d a r dP r e s so f C h i n a ,2 0 1 4 . i nC h i n e s e 英文编辑黄刚 万方数据
展开阅读全文