资源描述:
第37卷第1期 2020年3月 Vo l . 37 No . 1 Ma r. 2020 bMg do i 10. 3963/j. issn . 1001 - 487 X. 2020.01 ・ 023 基于HAZOP与模糊理论的台阶爆破安全综合评价 叶海旷2,付威志李宁二王李管寫涛“2,王其洲“2 1.武汉理工大学资源与环境工程学院,武汉430070;2.矿物资源加工与环境湖北省重点实验室,武汉430070; 3.中南大学资源与安全工程学院,长沙410083 摘 要针对台阶爆破安全影响因素具有模糊性且难以量化的特点,提出了一种耦合危险性与可操作性分 析HAZOP和模糊理论的综合评价方法。首先通过建立HAZOP评价指标体系辨识台阶爆破中存在的危险 因素,分析危险因素产生的原因和造成的后果,运用层次分析法AHP确定各评价指标的权重;然后基于模 糊理论对台阶爆破作业安全进行定量描述,得到台阶爆破安全评分等级;最后以贵州省某土石方工程为例, 验证了该综合评价方法的科学性和实用性。 关键词台阶爆破;HAZOP;模糊算法;层次分析法;定量评价 中图分类号TD608 文献标识码A 文章编号1001 -487 X202001 -0152 - 07 Comprehensive uation of Bench Blasting Safety based on HAZOP 。 表2爆破管理单元HAZOP分析表 Table 2 HAZOP analysis for blasting operations manag ement unit n级单元引导词原因后果措施 技术水平No无证上岗 不了解爆破流程,不能正 确使用爆破设备及材料 加强作业人员安全教育 与培训 安全培训Less 三级教育、工种教育等不缺乏基本的安全常识,违实行岗前培训与考核制 符合要求章作业,造成事故度 安全责任制No 规章制度缺失、未按最新 要求制定 管理人员失职引发爆破 事故 落实安全责任制 检查制度Pa rt Of 未制定检查制度或检查 制度制定不合理 安全风险意识不足 制定详细严格的检查制 度,并落实 操作规程Pa rt Of安全操作规程不完善易引发事故加强作业人员操作培训 表3爆破设计单元HAZOP分析表 Table 3 Analysis of blasting operation desig n unit HAZOP II级单元 引导词 S s 蔽 安全设计依据 Pa rt Of 现场勘察与 数据采集 Pa rt 0 爆破参数设计 Reverse、 Ot her t ha n 预防保护措施 n T No、Less 安全核算 Pa rt Of 地质资料不详尽、不符合国家标 准 设计人员对爆破对象没有一定 的了解,现场勘察也是敷衍了事 延时、堵塞长度理、单响药量、起 爆方向、孔网设计参数、单耗选 取、抵抗线、起爆顺序不合理 预防保护措施设计缺少或不符 合要求 理论计算及公式的核实不到位 易发生爆破事故 爆破方案设计有误 造成爆破飞石事故 地震效应事故或产生有 毒气体和粉尘 爆破效果不好,容易发生 爆破事故 进行工程现场踏勘,严格按照 国标、行标与规范进行 详细掌握工程地质条件,进行 施工预处理 多方论证参数设置,制定严谨 施工方案 进行重复多次检验,确定预防 保护措施安全有效 做好理论演算和数据校准工 作 表4爆破施工单元HAZOP分析表 Table 4 HAZOP analysis for blasting production construction unit n级单元引导词原因后果措施 钻孔Less、Pa rt Of 钻孔长度或倾斜度不达标、打残眼、未对 环炮孔进谟测量验收 爆破飞石事故、残眼残留 炸药爆炸事故 严格按照设计方案逐孔 进行检查 装药 Less、Mo re、 Ot her t ha n 装药量过大或过小、装药工艺不合理或 未按规定装药,未严格检查起爆器材 地震效应、飞石 严格按照施工方案进 行,加强监管力度 堵塞Ot her t ha n堵塞长度、密实度不达标飞石、冲孔 严格按照施工方案,加 强炮孔堵塞检查 联网Ot her t ha n 未检查网路连接、起爆工艺不合理或未 按规定起爆,完成网路连接后,是否进 行再次检查并点算雷管、炸药用量 容易发生早爆、拒爆等爆 破事故 重复多次进行检验,确 保网路准确 防护警戒 Less、Pa rt Of、 Ot her t ha n 警戒范围不达标、未清场或清场不彻 底、爆信号不符合要求、未采取必要的 减震措施、未设置警示标志、警戒区内 未设置专职守卫人员、过早解除警戒 伤人事故、损坏机械设备 制定防护方案,指定负 责人,落实安全生产责 任制 第37卷第1期叶海旺,付威志,李宁,等基于HAZOP与模糊理论的台阶爆破安全综合评价155 12台阶爆破安全评价指标体系台阶爆破安全评价指标体系 基于HAZOP流程分析的基础上,运用AHP中 递阶层次结构理论,台阶爆破安全评价指标体系划 分为目标层、准则层、要素层三个层次,从上到下分 别用A、B、C表示三个层次,同一层次各要素用数字 1.2,3从左至右依次表示,准则层B]丄2、艮,要 素层C[ CC15 见表1 o 13评价指标权重计算评价指标权重计算 AHP是系统工程中常用的赋权方法,依据已建 立的台阶爆破安全评价指标体系,按照层次分析法 分析流程,构造两两比较判断矩阵,利用MATLAB 软件计算判断矩阵的最大特征值入““及其对应的特 征向量礦,作一致性检验,平均随机一致性指标RI 值依据表5选取,若检验系数C7 0. 1,则构造的矩 阵通过了 一致性检验,从四量化各指标的权 重mm , I级单元权重指标帀0. 294 0. 353 0. 353,爲 3,CR 0, H级指标权重,如表6所 ZJ* O 表5平均随机一致性指标RI Table 5 Averag e random consistency indicator RI n1234 5 6 7 891011 RI000.580.90 1.12 1.24 1.32 1.411.451.491.51 表6评价指标权重 Table 6 Weig ht of uation index 二级指标权重值陷特征值入碍检验系数CR 0. 307 C20. 231 C30. 1545.0000. 000 C40. 17 7 Cs0. 131 C60. 152 C70. 102 C80. 3965. 1440. 036 C90. 27 1 Go0. 07 9 Cn0. 138 c120. 260 130. 2605.0100. 002 C140. 260 C“0. 082 14模糊综合评判模糊综合评判 根据台阶爆破现场情况,评价因素U{U1 U2 ,其中E为评价指标体系中的各评价 因素,评语等级的论V { V2 v3 曲{优 良中差劣},因素S对评语片隶属度需要通过隶 属函数实现,要素层各评价因素的论域为实数集,选 择偏大型梯形模糊分布作为隶属度函数,如式1 所示,多个专家对要素层打分,将分值通过隶属度函 S 计算出其隶属度,其中%为专家打分,依据表7 中隶属度与评价指标的关系求出其评价等级,依次 得出要素层的评价向量。 .0 , 90 表7隶属度与评价等级的关系 Table 7 Relationship between membership deg ree and uation level 隶属度函数评价等级 S兀1优 2 亍WS() 1良 1 2 丁ws(x) y中 0 s()-差 S() 0劣 例如,设有k个专家参与打分,对某二级指标有 局人认为“优”某专家对该指标的打分为91分,由 S 可求得其隶属度为“优”,k2人认为“良”,, 有饥人认为“劣”,且弘 k,进而计算出该指标 ;1 的隶属度为{半 模糊综合评价数学模型的建立,一级指标 i 1,2,3,有/个二级指标,建立一级评价指标的 模糊关系矩阵R,,由模糊关系矩阵亿并结合二级指 标权重行列式Wc,依据式3进行模糊综合评价,得 到一级模糊评价矩阵N[N.,N2,N3r ,依据式4 对一级评价指标进行模糊综合评价得到台阶爆破评 价结果F。 156爆破2020年3月 R, 「讥加raj 2 “边 ra2 - raj \ . N 出。亿弘,他,弘 3 F Wbo N 4 式中R矩阵为专家评价打分构成;。为模糊算 子,进行模糊变换时要选择适宜的模糊合成算子,选 用在运算时兼顾各因素的M ,㊉算子[16471 o 为了能客观判定台阶爆破安全性等级,对评语 等级进行赋值,提高评价结果的科学性和可信度,令 S , 55 65 75 85 95,建立台阶爆破安全性 体系总得分数学模型,得到被评价对象的分值I则, 将分值与S对比,确定隶属等级。 5 叫丫化 s 5 [1 2应用实例 贵州省某土石方工程,整个工程飾约648万后, 填方约55万n ,石方爆破约500万m3,以该工程台 阶爆破为对象,建立台阶爆破安全评价体系,依据 15位专家学者对评价体系中各因素的打分计算各 因素的隶属度,对该工程石方爆破进行评价,汇总统 计结果,依据式1,计算二级指标隶属度,如表8 所示。 表8指标隶属度统计 Table 8 Indicator membership statistics 二级 指标 S i 劣 52 差 S s 中 s4 良优 000.30.60. 1 Q00.50.30. 10. 1 Q00.30.40.20. 1 c400.40.40. 10. 1 00.60.20. 10. 1 C60000.40.6 C700.40.40. 10. 1 Cs0000.40.6 C90.20.40.20. 10. 1 Cw00.60.20. 10. 1 C“00. 10. 10.60.2 C1200. 10.20.50.2 C“ 00.20.20.50. 1 C1400. 1 0.20.30.4 C1500.30.20.30.2 根据表8得岀二级评价矩阵坨 -000.30.60.1- 00.50.30.10.1 00.30.40.20.1 00.40.40.10.1 0.60.20.10.1- -0000.40. 00.4 0.40.10. R2 0000.40. 0.2 0.4 0.20.10. L 0 0.6 0.20.10. -00. 10. 10.60.2-| 00. 10.20.50.2 r3 00.20.20.50.1 00. 10.20.30.4 S 0.30.20.30.2J 根据表6,式3可得 Wo R1 00.311 0.3100.26 N2W2o R2 0.054 0.197 0.111 0.264 0.374 他二町。堆二0 0. 142 0. 186 0. 446 0. 226 田式4得 F Wo N 0.019 0.212 0.196 0.331 0.242 由式5得系统整体评分 5 U 若“F S; 采用加权平均法进行计算,得出系统整体评分 u为80. 65, I级指标单元得分7 5. 83、82. 07、 82.56,系统整体评分对应安全级别为“中”,该工程 台阶爆破存在风险,爆破管理单元得分较低,运用 HAZOP分析方法对爆破管理单元进行分析发现该 项目各类人员岗位职责和技术操作细则不完善,安 全标准不明确,检查制度未彻底落实,未建立安全监 督网络,对于管理人员和作业人员未进行经常性安 全教育。 3结论 1 针对台阶爆破的特点,引入了以系统工程 为基础的HAZOP危险性评价方法,将台阶爆破划 分为3个I级单元15个II级单元,利用引导词找出 爆破工艺中存在的危险因素,并对危险因素产生的 原因、导致的后果及采取的控制措施等进行分析。 2 基于HAZOP分析,建立了台阶爆破安全评 价指标体系,引用AHP确定各评价因素的权重,运 用模糊算法确定了台阶爆破安全等级,以事故预防 代替事故处理。 3 基于HAZOP的台阶爆破安全性模糊综合 第37卷第1期叶海旺,付威志,李 宁,等 基于HAZOP与模糊理论的台阶爆破安全综合评价157 评价,应用于某土石方爆破工程,能够识别台阶爆破 流程中潜在的风险,定量评价安全等级,具有很好的 准确性和实操性。 参考文献参考文献References [1] 屈金坡,孟军良,齐朝鹏预先危险分析方法在常峪铁 矿的应用[J].矿冶工程,2015,35424-26. [1 ] QU Jin -p o, MENG Jun -l ia n g, QI Zha o -p en g. Ap p l ic a t io n o p rel imin a ry ha za rd a n a l y sis in c ha n gy u iro n min e [ J ]. Min in g a n d Met a l l urgic a l En gin eerin g, 2015,35 4 24- 26. in Chin ese [2] 罗文贵,谢贤平,粟美维,等.爆破工程中电爆网络拒 爆事故分析及对策[J].安全与环境学报,2016, 162 20-24. [2] LUO Wen -gui,XIE Xia n -p in g,S U Mei-wei,et a l . Co un t er mea sures a ga in st t he p o wer-t ra n smissio n n et wo rk bl a st in g t hro ugh t he misfire a c c iden t a n a l y sis ba sed o n t he bl a st in g en gin eerin g t heo ry [ J ]. Jo urn a l o f S a fet y a n d En viro n men t ,2016,162 20-24. in Chin ese [3] 周令剑,王宏奇小型露天采石场爆破作业安全评价 方法的探讨[J].工业安全与环保,2006460-62. [3 ] ZHOU Lin g-Jia n , WANG Ho n g-q i. Disc ussio n s o n sa fet y a ssessmen t mea sures fo r bl a st in g a t sma l l st o n e q ua rries [J ]. In dust ria l S a fet y a n d En viro n men t a l Pro t ec t io n, 20064 60-62. in Chin ese [4] 崔晓荣.基于鱼骨图模型的露天煤矿高温爆破施工安 全研究[J].煤矿安全,2018,493 242-244,24 [4] CUI Xia o -ro n g. S a fet y a n a l y sis o high t emp era t ure bl a s t in g in o p en c a st c o a l min e ba sed o n fishbo n e mo del [J]. S a fet y in Co a l Min es ,2018,49 3 242-244,248. in Chin ese [5] 张景钢,安美秀.HAZOP在煤矿安全评价中的应用研 究[J].矿业安全与环保,2016,435111-114. [5] ZHANG Jin g-ga n g, AN Mei-xiu. Ap p l ic a t io n o f HAZOP in c o a l min e sa fet y eva l ua t io n [ J ]. Min in g S a fet y En viro n men t a l Pro t ec t io n,2016,435 111-114. in Chin ese [6] J S ECCATORE,C ORIGLIAS S O,G DE TOMI. Assessin g a risk a n a l y sis met ho do l o gy fo r ro c k bl a st in g o p era t io n s [C] //The 10t h In t ern a t io n a l S y mp o sium o n Ro c k Fra g men t a t io n by Bl a st in g. New Del hi ,2013 51-60. [7 ] 周荣义,李石林,刘何清 HAZOP分析中LOPA的应用 研究中国安全科学学报,2010,207 76-81. [7 ] ZHOU Ro n g-y i, LI S hi-l in , LIU He-q in g. S t udy o n a p p l i c a t io n o f LOPA in HAZOP [ J ]. Chin a S a fet y S c ien c e Jo urn a l ,2010,207 7 6-81. in Chin ese [8 ] DUNJO J, FTHENAKIS V, VILCHEZ JA, ARNALDOS J. Ha za rd a n d o p era bil it y HAZOP a n a l y sis. A l it era t ure re view [J ]. Jo urn a l o Ha za rdo us Ma t eria l s, 2010, 17 3 1- 3 19-32. [9]邱星栋,黄 坤,陈彰兵,等.甘醇脱水装置的HAZOP 分析[J].天然气工业,2016,309123-128. [9] QIU Xin g-do n g, HUANG Kun, CHEN Zha n g-bin g, et a l . HAZOP a n a l y sis o n a TEG dehy dra t io n un it A c a se st udy [J]. Na t ura l Ga s In dust ry , 2016,30 9 123-128. in Chin ese [10] 吴世斌,林宏,刘俊松,等浅谈白鹤滩大型水电工 程施工本质安全管理实践[J].中国安全生产科学技 术,2017 ,13S 2 139-143. [10] WU S hi-bin, LIN Ho n g, LIU Jun -so n g, et a l . Disc ussio n o n t he essen t ia l sa fet y ma n a gemen t p ra c t ic e fo r l a rge hy dro p o wer p ro jec t in Ba ihe Bea c h [ J ]. Jo urn a l o f S a fet y S c ien c e a n d Tec hn o l o gy, 2017 ,13 S 2 139-143. in Chin ese [11] 朱淑敏.基于AHP和S CL的露天矿山爆破安全评价 研究[D].武汉武汉理工大学,2014. [11 ] ZHU S hu-min . Resea rc h o n sa fet y eva l ua t io n o f o p en -p it min e bl a st in g ba sed o n AHP a n d S CL[ D]. Wuha n Wu ha n Un iversit y o f Tec hn o l o gy ,2014. in Chin ese [12] 梁光川,何慧娟,何 莎,等.HAZOP技术在CNG储 气井风险评价中的应用[J].石油与天然气化工, 20151 99-102. [12] LIANG Gua n g-c hua n, HE Hui-jua n, HE S ha ,et a l . Ap p l ic a t io n o f HAZOP t ec hn o l o gy in risk a ssessmen t o f CNG wel l s [ J]. Chemic a l En gin eerin g o f Oil a n d Ga s, 20151 99-102. in Chin ese [13] 赖连明.公路工程施工中的安全管理与风险控制研 究[D].南京东南大学,2017 . [13 ] LAI Lia n -min g. Resea rc h o n S a fet y Ma n a gemen t a n d Risk Co n t ro l in Highwa y En gin eerin g Co n st ruc t io n[ D]. Na n jin gS o ut hea st Un iversit y,2017 . in Chin ese [14] 王 涛,侯克鹏,郭振世,等.层次分析法AHP在尾 矿库安全运行分析中的应用[J].岩土力学,2008, 29S I 680-686. [14] WANG Ta o , HOU Ke-p en g, GUO Zhen -shi, et a l . Ap p l i c a t io n o f a n a l y t ic hiera rc hy p ro c ess t o t a il in gs p o n d sa fe t y o p era t io n a n a l y sis [ J ]. Ro c k a n d S o il Mec ha n ic s, 2008,29S I 680-686. in Chin ese [15] 刘淑金,朱国庆,张 磊.基于Ma t l a b层次分析法的 建筑火灾人员安全疏散系统研究[J].安全与环境工 程,2014,214 110-114. [15] LIU S hu-jin, NIU Guo -q in g, ZHANG Lei. Resea rc h o n t he p erso n n el eva c ua t io n sy st em fo r buil din g fire ba sed o n ma t l a b AHP[ J]. S a fet y a n d En viro n men t a l En gin eer in g,2014,21 4 110-114. in Chin ese [16] 郭丽杰,王 楠,康建新.石化装置HAZOP节点重要 度模糊综合评价研究中国安全科学学报,2015, 251 98-104. [16] GUO Li-jie,WANG Na n ,KANG Jia n -xin . Fuzzy c o mp re 158 爆破 2020年3月 hen sive eva l ua t io n o f HAZOP n o de imp o rt a n c e fo r p et ro c hemic a l p l a n t[ J]. Chin a S a fet y S c ien c e Jo urn a l ,2015 , 251 98-104. in Chin ese [17 ]王玉杰,黄平路,张惠聚.台阶爆破作业安全的模糊 评价[J].武汉理工大学学报,2004868-7 0,7 3. [17 ] WANG Yu-jie, HUANG Pin g-l u, ZHANG Hu-ju. Fuzzy Co mp rehen sive Eva l ua t io n o n t he S a fet y o f S o il a n d Ro c k Bl a st in g[ J]. Jo urn a l o Wuha n Un iversit y o Tec hn o l o gy ,2004 8 68-7 0,7 3. in Chin ese [18]王悦,宋文华,宋相杰,等.爆炸危险环境电气系统 安全评价方法研究[J].南开大学学报自然科学 版,2016,4957 9-83. [18] WANG Yue,S ONG Wen -hua,S ONG Xia n g-jie,et a l . S t udy o n sa fet y eva l ua t io n met ho d o f el ec t ric a l sy st em in exp l o sive a t mo sp heres [ J ]. Ac t a S c ien t ia rum Na t ura l ium Un iversit a t is Na n k a ien sis,2016,495 7 9-83. in Chin ese 英文编辑任高峰 (上接第151页) [7 ]黄润秋,李果,巨能攀.层状岩体斜坡强震动力响应 的振动台试验[J].岩石力学与工程学报,2013, 325 865-87 5 [7 ] HUANG Run -q iu, LI Guo, JU Nen g-p a n . S ha k in g t a bl e t est o n st ro n g ea rt hq ua k e resp o n se o f st ra t ified ro c k sl o p es [J]. Ro c k Mec ha n ic s a n d En gin eerin g, 2013,32 5 865-87 5. in Chin ese [8 ] LI J, DENG S , WANG M, et a l . Wea k dist urba n c e-t rig gered seismic even t s a n exp erimen t a l a n d n umeric a l in vest iga t io n [J ]. Bul l et in o f En gin eerin g Geo l o gy t he En viro n men t ,2019,7 8 4 2943-2955. [9] 肖卫国,王肖钧,劳 俊不同方式地下爆炸地震耦合 效应的数值模拟[J].计算物理,2011,286 797-802. [9 ] XIAO Wei-guo , WANG Xia o -jun, LAO Jun . Numeric a l simul a t io n o f seismic c o up l in g effec t s in un dergro un d ex p l o sio n s wit h differen t exp l o sio n mo des [J]. Co mp ut a t io n a l Phy sic s, 2011,286 7 97 -802. in Chin ese [10] 宗 琦,汪海波,周胜兵.爆破地震效应的监测和控 制技术研究[J].岩石力学与工程学报,2008,275 938-945. [10] ZONG Qi, WANG Ha i-bo, ZHOU S hen g-bin g. Resea rc h o n mo n it o rin g a n d c o n t ro l l in g t ec hn iq ues c o n siderin g effec t s o f ssmic sho c k[ J]. Ro c k Mec ha n ic s a n d En gi n eerin g ,2008,27 5 938-945. in Chin ese [11] 林大超,白春花,张 奇爆炸地震地面竖向振动的 幅值特性[J].北京北京理工大学学报,2001, 212 17 3-17 6. [11 ] LIN Da -c ha o , BAI Chun -hua , ZHANG Qi. Amp l it ude c ha ra c t erist ic s o vert ic a l gro un d surfa c e vibra t io n c a used by exp l o sio n seism[ J]. Tra n sa c t io n s o Beijin g In st it ut e o f Tec hn o l o gy ,2001,21 2 17 3-17 6. in Chin ese [12] BERGKVIS T N O,FERM R. Nuc l ea r exp l o sio n s 1945- 1998 [ R ]. Defen c e Resea rc h Est a bl ishmen t , 2000 14- 15. [13] 殷跃平,王 猛,李 滨,等汶川地震大光包滑坡动 力响应特征研究[J].岩石力学与工程学报,2013, 3110 1969-1982. [13] YIN Yue-p in g, WANG Men g, LI Bin , et a l . Dy n a mic re sp o n se c ha ra c t erist ic s o f Da gua n gba o l a n dsl ide t riggered by Wen c hua n ea rt hq ua k e[ J]. Ro c k Mec ha n ic s a n d En gin eerin g,2013 ,31 10 1969-1982. in Chin ese [14] 沈 蔚,徐全军,黄文华.用灰关联分析法确定爆破 震动参数的主要影响因素[J].工程爆破,2000, 64 8-13. [14] S HEN Wei, XU Qua n -jun, HUANG Wen -hua . Det ermi n a t io n o f ma in in fl uen c in g fa c t o rs o n bl a st in g vibra t io n p a ra met ers by grey c o rrel a t io n a n a l y sis [ J]. En gin eerin g Bl a st in g,2000,64 8-13. in Chin ese [15] 张 跃,皱寿平,宿芬.模糊数学方法及其应用[M]. 北京煤炭工业出版社,1992. [15] ZHANG Yue, ZHOU S ho u-p in g, S U Fen . Fuzzy ma t he ma t ic s met ho d a n d a p p l ic a t io n [ M ]. Beijin g Co a l In dus t ry Press, 1992. in Chin ese [16] 王俊平爆破地震波对周围物影响的分析[D].武汉 武汉理工大学,2005. [16] WANG Jun -p in g. An a l y sis o f t he in fl uen c e o f bl a st in g seismic wa ve o n n ea rby st ruc t ures [ D]. Wuha n Wuha n Un iversit y o Tec hn o l o gy,2005. in Chin ese [17 ] 郑 峰爆破地震效应影响因素的研究及工程应用 [D].武汉武汉科技大学,2007 . [17 ] ZHENG Fen g. S t udy a n d a p p l ic a t io n o n in fl uen c e fa c t o rs fo r t he bl a st in g seismic effec t[ D]. Wuha n Wuha n Un i versit y o f S c ien c e a n d Tec hn o l o gy,2007 . in Chin ese [18] 刘华丽,卢厚清,李宏伟,等.模糊综合评判法在边坡 稳定性分析中的应用[J].解放军理工大学学报学 报,2013,141 84-88. [18] LIU Hua -l i,LU Ho u-q in g,LI Ho n g-wei,et a l . Ap p l it io n o f n ew met ho d ba sed o n fuzzy c o mp rehen sive a ssessmen t fo r st a bil it y o f l a n dsl ide [ J]. PLA Un iversit y o S c ien c e a n d Tec hn o l o gy ,2013,14 1 84-88. in Chin ese 英文编辑陈东方
展开阅读全文