基于个人风险分布计算的泄漏事故疏散路线优化

陈晓坤 何春旭 史合 程方明 苏彬

摘 要：发生有毒有害氣体泄漏事故后，最优疏散线路的选择是降低人员伤亡的关键。本研究提出一种基于疏散路线个人风险分布特征分析的泄漏事故中疏散路线优化方法，考虑泄漏事故疏散路线安全性和疏散时间等影响因素，建立了疏散路径动态优化选择机制。将泄漏物料特性参数、泄漏源位置、存储状态参数、风向等事故发生时可快速获取的因素作为已知条件，并引入点火源、人员伤亡概率等不确定性因素的统计频率，给出疏散路线综合风险值的计算方法，进而实现了疏散路线安全性的量化表征。以某化工厂区液氨储罐泄漏事故为例，采用Safeti软件计算分析了各疏散路线个人风险，应用疏散路线动态优化选择机制，优选出针对该事故场景的最佳疏散路线。相关研究结果可为有毒有害气体泄漏事故应急疏散的科学决策提供理论指导。

关键词：泄漏事故;疏散;路线优化;个人风险中图分类号：X 936

文献标志码：A

文章编号：1672-9315（2021）01-0008-07

DOI：10.13800/j.cnki.xakjdxxb.2021.0102

Optimization of leakage evacuation route based on

personal risk distribution calculation

CHEN Xiaokun1，HE Chunxu1，3，4，SHI He2，CHENG Fangming1，3，4，SU Bin1，3

（1.School of Safety Science and Engineering，Xian University of Science and Technology，Xian 710054，China;

2.Xian Institute of Space Radio Technology，Xian 710100，China;

3.Shaanxi Engineering Research Center for Industrial Process Safety & Emergency Rescue，Xian 710054，China;

4.Urban Public Safety Research Institute，Xian University of Science and Technology，Xian 710054，China）

Abstract：After the occurrence of a toxic and harmful gas leakage accident，the choice of the optimal evacuation route is the key to reduce casualties.In this study，an evacuation route optimization method is proposed based on the analysis of individual risk distribution characteristics of evacuation routes.Considering the influencing factors such as leakage accident evacuation route safety and evacuation time，the dynamic optimization selection mechanism of evacuation path is established.The factors possible to be quickly obtained in accidents，such as leakage material characteristic parameters，leakage source location，storage state parameters，and wind direction，are taken as known conditions，and the statistical frequency of uncertain factors such as ignition source and casualty probability are introduced to design the evacuation route.The calculation method of the comprehensive risk value further realizes the quantitative characterization of the safety of the evacuation route.Taking the leakage accident of liquid ammonia storage tank in a chemical plant area as an example，the personal risk of each evacuation route was calculated and analyzed by Safeti software，and the dynamic optimization selection mechanism of evacuation route was applied to optimize the optimal evacuation route for the accident scene.The relevant research results can provide theoretical guidance for the scientific decision-making of emergency evacuation of toxic and harmful gas leakage accidents.

Key words：leak accident;evacuation;route optimization;individual risk

0 引 言

工业生产过程中泄漏事故发生后疏散路线选择是应急决策的重要问题之一，最优的疏散路线可以有效降低人员伤亡。然而，疏散路线受多种因素影响且关系复杂，优选过程不仅要考虑路线疏散时间和路线安全性，更要根据现场具体情况选择符合当前情况的疏散路线[1-5]。事故现场气象条件的变化、次生灾害事故的诱发及演变等都要求疏散路线需要随泄漏事故发展而不断变化的，疏散路线在空间和时间上均具有风险的不确定属性[6-10]。

目前，已有诸多学者针对危险化学品泄漏事故最优疏散路线选取开展研究。如赵婧璇等以某氯碱厂液氯泄漏为研究对象，采用计算流体力学Fluent 软件模拟计算得出区域疏散路线中各监测点的氯气扩散实时浓度，优选出在所有泄漏场景中人员累计中毒风险较小的最优疏散路线[11]。赵江平等为了了解在毒气泄漏情况下的中毒死亡概率及避难点的安全情况[12]，采用高斯模型计算毒气泄漏的浓度，并且使用Building EXODUS 软件模拟区域应急疏散过程。熊立春等引入时间因素改进优化了高斯扩散模型，并且建立风险矩阵、人口密度矩阵，为选取最佳疏散路线提供了依据[13]。郑杨炜等分析了氨的危害，通过MATLAB进行模拟，采用高斯烟团模型，从而划分出危险区域的范围[14]。LAW等通过使用计算流体力学CFD和ALOHA软件模拟马来西亚某工业园区周围的氯气泄漏的扩散情况，最终可知风向、风速、地表面貌会影响氯气的扩散模式，安全疏散路线应该根据季风季节而变化[15]。GAI等通过考虑在紧急决策中可能相互冲突的多个目标，针对有毒云雾释放的紧急情况优化疏散路线选择[16]。目的是将对个人和公众的潜在风险降到最低，并将个人疏散过程中沿着一条路径的总疏散时间降到最低。可见，目前相关研究多以泄漏后有毒有害气体浓度场的预测为基础，研究规划安全高效的疏散路线，多关注中毒等单一危害后果对疏散路线的影响，而对疏散过程火灾、爆炸等次生灾害事故的风险关注不多。因此，考虑泄漏发生后中毒、火災、爆炸等多种事故风险叠加对疏散路线安全性的影响，开展疏散路线优化研究对泄漏事故应急疏散过程的科学决策具有重要指导意义。

以泄漏事故人员疏散为研究对象，考虑路线疏散时间与路线安全性等主要因素，建立动态疏散路线优化选择机制。在疏散路径安全性方面，关注泄漏发生后多种事故风险叠加的影响，提出疏散路线综合风险值的概念，并给出计算方法，实现了疏散路线安全性的定量分析。最后，以某化工厂区液氨储罐泄漏事故为例，使用Safeti软件计算分析了不同疏散路线的综合风险值，采用动态疏散路径优化选择机制，确定了最优疏散路径，可为化工厂区泄漏事故疏散路径优选提供参考。

1 影响疏散路线选择的因素分析

1.1 疏散时间

疏散时间是指在泄漏事故发生后，处于多种事故风险叠加影响区域内的人群疏散到安全区域内所用的时间。若事故风险叠加影响区域内的人群疏散到安全区域内且人员无伤亡，则可以认为该疏散路线的选取是合理的。由于人员在多种事故风险叠加影响区域内的位置不同，不同人员对待泄漏事故的应急反应能力也不相同，泄漏事故发生时人员心里活动状态也有所差异，事故发生时路线的通畅程度也无法确定，这些都将影响整个疏散时间。因此，疏散时间对个体而言是随机变量，但在本研究中，将疏散时间看成是整体疏散时间，默认人群疏散位置、疏散能力等主观条件相同，是相对固定的变量[17-18]。

1.2 路线安全性

路线安全性是指人员在疏散路线上免受泄漏发生后多种事故风险叠加影响的能力[9]。确保疏散人员安全是疏路线选择的前提条件，泄漏事故疏散过程面临中毒、火灾、爆炸等多种事故风险，有毒有害气体浓度分布变化、潜在点火源位置等情况的不同导致事故灾区的区域风险在空间和时间上存在差异，不同疏散路线的安全性也是不同的。因此，选择风险较低的疏散路线，可提高疏散安全性，减少事故灾变造成的人员伤亡。

2 疏散路线安全风险

2.1 疏散路线安全性分析

路线安全性可采用路线风险大小进行衡量。风险是由事故后果与事故发生的频率确定的，泄漏事故疏散路线风险则由泄漏事故演化过程中的多种潜在事故后果与事故发生的频率共同确定。泄漏事故发生后，泄漏孔径、泄漏速率、风向、风速等诸多事故后果影响因素是能确定的，但点火源、大气稳定度等影响因素是很难直接确定的，且直接影响事故演变危害特性与风险大小。为了对路线风险进行定量分析，文中引入了疏散路线综合风险值这一概念来描述疏散过程各类事故风险的叠加，进而反映疏散路线的安全性。

2.2 疏散路线综合风险值

疏散路线综合风险值就是疏散个体在疏散路线上可能受到危险后果及某种程度伤害的频发程度。疏散路线综合风险值实质上就是个体在路线上的伤亡风险，体现了现场诸多事故后果危害不确定因素概率的叠加，其风险值由疏散人员在疏散路线上的个人风险与疏散路线长度所决定。

2.3 疏散路线综合风险值的计算

疏散路线综合风险值量化表示了路线的安全性，以某一泄漏事故场景为例，对疏散路径综合风险值的具体计算方法说明如下。

1）假设某一泄漏场景发生。

2）选择某一种风速V和大气稳定度D，从而给出这个天气等级M，以及该天气等级下的一种风向考虑风速、风向等条件能够迅速确定，而大气稳定度难以实时监测，大气稳定度频率为PD。

3）如果是可燃物泄漏，选择一个点火事件i并确定点火概率为Pi。如果考虑物质毒性影响，则不考虑点火事件。

4）泄漏事故发生后，在泄漏孔径、泄漏速率、风速、风向等能够快速确定的情况下，计算大气稳定度PD，点火事件（针对可燃物）条件下网格单元上的死亡概率Pd，计算中参考高度取1 m。

5）计算（D，i）条件下网格单元的个体风险（individual risk，IR）式（1）。

ΔIR=PD×Pi×Pd

（1）

6）对所有的天气事件重复2）～5）步的计算，对所有的点火事件重复3）～5）步的计算;则得到网格点处的叠加的个人风险，式（2）

IR=∑D∑iΔIR

（2）

7）计算出疏散线路上每个网格点处的个人风险值，即得到该疏散线路L上的个人风险分布曲线，并且对该个人风险分布曲线进行路线长度的积分，即可得到该路线上的综合风险值P，式（3）

P=∫LIR

（3）

在文中使用了风险定量分析软件Safeti，对目标疏散路线进行风险计算，从而确定多种事故影响的叠加风险。

3 疏散路径优化选择方法

不同事故场景下，不同的疏散路径有着不同的疏散时间，并且各条疏散路线的安全性也各不相同。因此，建立一种动态疏散路线优化选择机制是至关重要的。文中通过图1所示的疏散路线优化选择流程，表示综合考虑了事故发生时的疏散时间、路线安全性等因素下的动态选择机制，对每一条疏散路线进行综合筛选，最终选择出最优的疏散线路。

在泄漏事故发生后，应急决策人员采用动态疏散路线选择机制，选出在该泄漏事故发生时的最优疏散路线。从图1可知，应急决策人员应该迅速对每条疏散路线进行安全性分析，选择出其中安全性较高的疏散路线。然后在安全性较高的疏散路线中选择出疏散时间较小的疏散路线。

通过以上的疏散路线优化选择，应急决策人员能够根据具体泄漏事故的情况动态考虑路线安全性与疏散时间。突破了以往多关注疏散时间的优化方案，可有效避免选择二次事故及次生事故演化影响区域内的疏散路径，有效地减少人员伤亡，为应急疏散路线优选决策提供了一种方法[20]。

4 应用实例

4.1 厂区概况

某化工厂区平面图如图2所示，根据其场内人员分布可以分为3个区域：办公区（200人）、住宿区（450人）、操作区（50人）。厂区周边点火源可辨识为内部的若干條电线，厂区外的位于道路两侧的两条高压电线，厂区外道路上的车辆，同时认定人员所在区域为弱点火源。企业应急预案中确定了厂区内的疏散线路共有A，B，C，D等4条，分别去往四周的安全疏散点。已知液氨储罐为厂区主要危险源，共有6个含氨设备，其具体位置均以在厂区平面图中以设备标记形式表示，具体设备参数条件见表1。

4.2 疏散路线综合风险计算与分析

假设该企业某日风向为北风，1号液氨储罐发生了泄漏，据现场人员估计泄漏孔径约有100 mm。

根据储罐实时存量实时监测可知，事故发生时储罐内约有552.5 m3的液氨。并且场内风速监测仪显示实时风速为2.5 m/s。在泄漏事故发生后，泄漏孔径、风速、风向等因素能够迅速根据确定。然而还存在着点火概率、大气稳定度、死亡概率等诸多不能确定的参数[21]。在这种泄漏事故发生后，应急决策人员在选择疏散路线时应当考虑路线安全性，避免人员在疏散过程中误入危险区域。

通过使用Safeti来实现疏散路线 安全性计算，该软件通过输入相关工艺参数、设备参数、气象参数等，模拟特定事故情景，得出不同事故情景模式的影响范围，其不仅能构建单一设备设施的事故模型，并且对于不同设备之间的相互影响而造成的区域性影响也能进行模拟。从而能够使用计算机来对既定的事故模型进行影响范围及危害程度的计算，形成模拟分析图和详细的模拟结果报告。[22-25]

通过软件的风险计算情况能够得到如图4所示的事故风险等值线。图4明确表示了该事故下不同风险区域影响的位置。根据《GB36894—2018危险化学品生产装置和储存设施风险基准》中对于这种厂区的规定可容许个人风险值可知，当风险为3×10-6时，图中绿色区域是泄漏事故风险影响区域。考虑事故发生时，毒性物质影响区域相对较大，因此采用美国国家咨询委员会（national advisory committee，NAC）与国家研究委员会（national research council，NRC）针对国家、地方政府以及个人企业处理包括泄漏、灾难性暴露等紧急情况所制定的急性暴露标准AEGLS值。在文中中根据事故疏散长度选取AEGL-2在30 min的关注浓度220 ppm作为毒性事故影响区域的浓度值，如图中蓝色区域所示。可以看到风险区域及毒性后果影响区域均偏向南方，风险区域面积为498 680.15 m2，毒性后果影响区域面积为107 413 m2，这表明了在泄漏事故发生后仅仅关注事故影响区域是不够的，应急决策人员应当关注事故风险，从而选择出合理可行的路线。

在泄漏事故发生后，应急决策人员不仅要考虑事故后果影响区域，更要考虑事故风险区域，从而更加准确的根据事故情况决定人员疏散路线，因此文中引入了疏散路线综合风险值这一概念，通过模拟事故发生时情景，从而准确判断路线风险大小，从而为路线应急决策提供依据。根据前文的路线风险计算方法，结合软件计算本次事故各疏散线路的情况，从而得到图5所示的各疏散路线个人风险分布。

图5中横坐标为疏散路线行进的距离，纵坐标是风险值大小。可见，在泄漏事故发生后，人员所处的初始位置因为靠近事故储罐，因此初始风险都相对较大。从疏散路线A，B的风险趋势图上能看到，路线风险大小呈现波动变化趋势。这就说明事故发生后人员在该路线疏散时，会再次进入风险程度较高的区域内。

图4中绿色曲线代表风险值为3×10-6的风险区域，根据图4风险等高线可知，内部还分布诸多风险更高的区域。疏散路线A的路线不规则，在图4中位于风险区域与事故后果区内，并且能够看到随着路线进行疏散，路线跨越了不同的风险区域，因此造成其风险趋势波动变化。在路线B上进行疏散时，人员随着路线会再次进入高风险区域内，然后才能到达避难点，其路线的风险趋势图就是先增加在减小的趋势。

从风险趋势变化上证明疏散路线A，B在疏散过程中引发二次伤害可能性较大，不适合作为人员应急疏散的路线。疏散路线C，D在人员初始位置处风险处于最大，随着人员开始疏散，沿疏散路线个人风险值持续降低。对厂区已有的疏散线路综合风险值进行计算，疏散路线A，B，C，D的风险值分别为30.078，20.808 8，2.321 2，1.019 2。

4.3 疏散线路优化选择

根据已经建立的疏散路线动态选择机制，通过对路线安全性、疏散时间等参数指标的优化选择，从而可以分析选择出最优的疏散路线。考虑人员在应急疏散时，均采用最大移动速度，在文中中默认为1.5 m/s。本次泄漏事故中泄漏孔径为100 mm，风速为1.5 m/s，风向为北风的情况下的各条疏散线路综合评价见表2。

根据疏散路线动态优化选择机制，首先比较4条路线的综合风险值的大小，对路线综合风险值数值量级进行分类，发现疏散路线C，D在安全性上是满足要求的。其次考虑路线长度与疏散时间，疏散路线越长，人员所花费的时间也就越长，可能造成人员伤害的可能性也就越大。可见，疏散路线D所用疏散时间较短，因此确定为首选疏散路线。

5 結 论

1）针对泄漏事故疏散路线的科学优选问题，文中提出了一种基于疏散路线个人风险分布特征分析的疏散路线优化选择方法，该方法不仅关注单一泄漏事故的危害影响范围，而且综合考虑了现场环境条件下有害物料泄漏扩散过程可能发生的中毒、火灾、爆炸等多种事故的风险叠加影响，可有效提高疏散路线的安全性。

2）通过引入疏散路线综合风险值的概念来定量表征疏散过程安全性，并给出了具体计算方法，为疏散路线安全性定量分析提供了理论依据。

3）以某化工厂区液氨储罐泄漏安全疏散为具体应用实例，详细阐述了疏散路线综合风险计算及路线优化选择过程，可为有毒有害气体泄漏事故的应急疏散提供重要参考。

参考文献（References）：

[1] 潘旭海，蒋军成.重（特）大泄漏事故统计分析及事故模式研究[J].化学工业与工程，2002，19（3）：248-252.PAN Xuhai，JIANG Juncheng.Analysis on important release accidents and modes studying[J].Chemical Industry and Engineering，2002，19（3）：248-252.

[2]罗振敏，苏彬，程方明，等.液氨泄漏事故的后果模拟与风险评估[J].安全与环境工程，2017，24（6）：149-154.LUO Zhenmin，SU Bin，CHENG Fangming，

et al.Consequence simulation and risk assessment on leakage accident of liquefied ammonia[J].Safety and Environmental Engineering，2017，24（6）：149-154.

[3]潘东.液氨储罐事故后果模型分析及技术改造思路[J].过程工程学报，2018（1）：67-71.PAN Dong.Model analysis of liquid ammonia tank accident consequence and thinking of technical transformation[J].The Chinese Journal of Process Engineering，2018（1）：67-71.

[4]孙晖，张树海，吴明亮.液氨泄露蒸气云爆炸的风险分析[J].安全与环境工程，2010，17（4）：64-66.SUN Hui，ZHANG Shuhai，WU Mingliang.Risk analysis of liquid ammonie vapour cloud explosion[J].Safety and Environmental Engineering，2010，17（4）：64-66.

[5]李高艳，业成，支有冉，等.液氨泄漏事故伤害范围模拟及人员疏散研究[J].工业安全与环保，2015，41（5）：12-14.LI Gaoyan，YE Cheng，ZHI Youran，et al.Liquid ammonia leakage accident injury range simulation and evacuation[J].Industrial Safety and Environmental Protection，2015，41（5）：12-14.

[6]董晓强，李坤，刘玉婷，等.氨制冷机房氨气泄漏扩散的数值模拟[J].制冷学报，2017，38（6）：12-19.DONG Xiaoqiang，LI Kun，LIU Yuting，et al.Numerical simulation of the ammonialeakage diffusion in ammonia refrigeration room[J].Journal of Refrigeration，2017，38（6）：12-19.

[7]周宁，陈力，吕孝飞，等.大型LNG 储罐连续泄漏扩散过程影响因素的分析[J].化工进展，2019，38（10）：4423-4436.

ZHOU Ning，CHEN Li，LU Xiaofei，et al.Analysis on influencing factors of LNG continuous release & dispersion process from large scale tank[J].Chemical Industry and Engineerin Progress，2019，38（10）：4423-4436.

[8]王媛原.液氨泄漏事故的预防及应急救援对策[J].广东化工，2017，22（44）：114-115.WANG Yuanyuan.Study on the prevention and emergency rescue countermeasure of ammonia leakage accident[J].WWW.Gdchem.com，2017，22（44）：114-115.

[9]徐志胜，孔杰，游温娇，等.单洞双线铁路隧道火灾人员疏散安全性分析[J].中国安全生产科学术，2019，15（1）：181-186.XU Zhisheng，KONG Jie，YOU Wenjiao，et al.Analysis on safety of personnel evacuation in fire of single-hole double-line railway tunnel[J].Journal of Safety Science and Technology，2019，15（1）：181-186.

[10]王起全，王鸿鹏.基于情景构建的危化品事故应急疏散模拟研究[J].中国安全科学学报，2017，27（12）：147-152.WANG Qiquan，WANG Hongpeng.A scenario construction based study on emergency evacuation in hazardous chemicals accident[J].China Safety Science Journal，2017，27（12）：147-152.

[11]赵婧璇，赵云胜，郭颖.氯碱厂液氯泄漏事故人员疏散路径优化的数值模拟[J].安全与环境工程，2018，25（12）：150-153.ZHAO Jingxuan，ZHAO Yunsheng，GUO Ying.Numerical simulation on optimization of liquid chlorine leakage and evacuation pathway in a chlor-alkali plant[J].Safety and Environmental Engineering，2018，25（12）：150-153.

[12]趙江平，田海啸，宋倩文，等.氯气泄漏事故区域应急疏散研究——以西安某化工厂氯气储罐为例[J].中国安全科学学报，2014，24（1）：163-169.ZHAO Jiangping，TIAN Haixiao，SONG Qianwen，et al.Research on emergency evacuation in chlorine leak accident：taking a chlorine tank in a certain chemical factory in Xian as an example[J].China Safety Science Journal，2014，24（1）：163-169.

[13]熊立春，陈建宏，石东平.引入时间叠加的高斯液氨泄漏扩散模拟及人员疏散研究[J].中国安全生产科学技术，2013，11（11）：76-82.XIONG Lichun，CHEN Jianhong，SHI Dongping.Simulation on liquid ammonia leakage by Gaussian model introducing temporal superposition and personnel evacuation[J].Journal of Safety Science and Technology，2013，11（11）：76-82.

[14]郑杨炜，谢殿荣，何衍兴，等.基于MATLAB的液氨泄漏环境风险分析[J]安全与环境工程，2017，17（1）：41-44.ZHENG Yangwei，XIE Dianrong，HE Yanxing，et al.Environmental risk analysis of the leakage from liquid ammonia tank by MALABT[J]Safety and Environmental Engineering，2017，17（1）：41-44.

[15]LAW W P，ERAIN N，RAMLI N I.Assessment of chlorine leak dispersion around Gebeng industrial area and potential evacuation route[J].Atmospheric Research，2019，216：117-129.

[16]GAI Wenmei，DENG Yunfeng，

JIANG Zhongan.Multi-objective evacuation routing optimization for toxic cloud releases[J].Reliability Engineering & System Safety，2017，159：58-68.

[17]GAI W M，DENG Y F.Survey-based analysis on the diffusion of evacuation advisory warnings during regional evacuations for accidents that release toxic vapors：a case study[J].Journal of Loss Prevention in the Process Industries，2019，57：174-185.

[18]YUKIE L，MASAHIRO L.An evacuation route planning for safety route guidance system after natural disaster using multi-objective genetic algorithm[J].Procedia Computer Science，2016，96：1323-1331.

[19]CHEN P Z，CHEN G H，WANG L W.Optimizing emergency rescue and evacuation planning with intelligent obstacle avoidance in a chemical industrial park[J].Journal of Loss Prevention in the Process Industries，2018，56：119-127.

[20]DOU Z，MEBARKI A，CHENG Y.Review on the emergency evacuation in chemicals-concentrated areas[J].Journal of Loss Prevention in the Process Industries，2019，60：35-45.

[21]罗振敏，李青，程方明.液氨储罐泄漏事故后果的影响因素研究[J].消防科学与技术，2018，37（3）：410-413.LUO Zhenmin，LI Qing，CHENG Fangming.Study on the influencing factors of accident sequence in leakage of liquid ammonia tank[J].Fire Science and Technology，2018，37（3）：410-413.

[22]郭群.CASST-QRA定量风险评价法在重大危险源评估中的应用[J].煤炭技术，2015，34（6）：296-298.GUO Qun.Application of CASST-QRA quantitative risk assessment method in the assessment of major hazard sources[J].Coal Technology，2015，34（6）：296-298.

[23]吴宗之，多英全，魏利军，等.区域定量风险评价方法及其在城市重大危险源安全规划中的应用[J].中国工程科学，2006，8（4）：46-49.WU Zongzhi，DUO Yingquan，WEI Lijun，et al.A regional quantitative risk assessment method and its application in urban dangerous source safety planning[J].Engineering Science，2006，8（4）：46-49.

[24]罗艾民，陈平，多英全，等.定量风险评价中的危险辨识方法[J].中国安全生产科学技术，2010，6（1）：54-57.LUO Aimin，CHEN Ping，DUO Yingquan，et al.Hazard identification method in quantitative risk assessment[J].China Safety Science and Technology，2010，6（1）：54-57.

[25]任峰濤，胡卓华.化工园区安全的定量风险评价[J].化工管理，2013（22）：41-42.REN Fengtao，HU Zhuohua.Quantitative risk assessment of safety in chemical industry parks[J].Chemical Enterprise Management，2013（22）：41-42.

收稿日期：2020-03-13   责任编辑：刘 洁

基金项目：

国家重点研发计划项目子课题（2017-YFC-0804702）

通信作者：

陈晓坤，男，辽宁沈阳人，教授，博士生导师，E-mail：Chenxk@xust.edu.cn