Integrated resilience of urban rail transit network with active passenger flow restriction under major public health disasters
-
摘要: 分析了重大公共卫生灾害对城市轨道交通网络集成韧性的影响机理;基于韧性曲线模型对传统韧性测度方法进行了修正,构建了面向重大公共卫生灾害影响的城市轨道交通网络集成韧性测度方法;评估了城市轨道交通网络节点重要度水平,运用复杂网络构建了城市轨道交通网络拓扑模型,对节点客流进行了模拟分配;应用SEZIR传染病传播模型模拟了灾害传播过程,研究了城市轨道交通在重大公共卫生灾害背景下的集成韧性水平演化规律;以西安市疫情发展过程为研究对象,对主动客流限制下城市轨道交通网络的集成韧性水平进行了模拟和数值分析。研究结果表明:主动客流限制措施能够有效提高城市轨道交通网络对重大公共卫生灾害的阻断能力,当客流限制水平达到30%后,重大公共卫生灾害传播过程趋于平缓;主动客流限制措施会直接导致城市轨道交通网络运行效率降低,但能够提升城市轨道交通网络在重大公共卫生灾害影响下的集成韧性水平;当客流限制水平分别为70%、40%和20%时,城市轨道交通网络集成韧性水平的改善提升效果更加明显,累积改善效果分别可达到10.73%、46.87%和226.81%。Abstract: The influencing mechanism of major public health disasters on the integrated resilience of urban rail transit network was analyzed. The traditional resilience measurement method was modified by the resilience curve model, and an integrated resilience measurement method was constructed for the urban rail transit network affected by major public health disasters. The importance levels of urban rail transit network nodes were evaluated. A topological model of urban rail transit network was constructed by the complex network approach to simulate and assign the nodal passenger flow. The SEZIR infectious disease spread model was applied to simulate the spread process of disaster, and the evolution laws of the integrated resilience level of urban rail transit in the context of a major public health disaster were studied. The process of epidemic development in Xi'an was taken as the research object, the integrated resilience level of the urban rail transit network under active passenger flow constraints was simulated and numerically analyzed. Research results show that the ability of the urban rail transit network to interrupt the spread of major public health disasters can be effectively enhanced by active passenger flow restriction measures. The spread process of major public health disasters becomes gentle after the restriction level of passenger flow reaches 30%. Active passenger flow restriction measures are able to directly reduce the operational efficiency of the urban rail transit network, but the integrated resilience level of the urban rail transit network under the influence of major public health disasters can be improved. The improvement of integrated resilience level of the urban rail transit network is more significant when the passenger flow restriction level is 70%, 40%, and 20%, and the cumulative improvement is 10.73%, 46.87%, and 226.81%, respectively.
-
表 1 西安市城市交通主动限流措施
Table 1. Active traffic restriction measures in Xi'an
交通类型 主动限流措施 实现目标 城市轨道交通 风险车站“越站”运行 规避高风险区域 加开备用列车 控制列车满载率 调整线网运营服务时间 降低车站拥挤度,控制列车满载率 分批控制进站 降低车站拥挤度 关闭车站部分出入口 降低车站拥挤度 调整发车间隔 降低车站拥挤度,控制列车满载率 城市公交 车辆满载率不超过50% 控制车辆满载率 早晚高峰期间增加配车 控制车辆满载率 管控区域站点不停靠 规避高风险区域 客运站点关停 规避高风险区域 私家车 暂停机动车尾号限行措施 降低道路拥挤度,控制公共交通满载率 表 2 自然灾害与重大公共卫生灾害下城市轨道交通网络韧性研究对比
Table 2. Comparison of urban rail transit network resilience studies under major public health disasters and natural disasters
灾害类型 自然灾害 重大公共卫生灾害 结构层面影响 网络节点中断 网络结构未受到破坏 网络节点失效 网络线路中断 功能层面影响 城市轨道交通网络自身功能遭受灾害破坏,最大客流承载能力受损 城市轨道交通网络自身功能未遭到破坏,实际客流承载能力受到灾害防控措施限制 城市轨道交通网络运行影响 结构、功能受到灾害破坏,导致网络运行效率降低 灾害防控措施导致网络客流承载能力不能充分发挥,实际功能运行效率降低,对灾害的抵抗能力增加 城市轨道交通网络韧性内涵 网络遭受灾害干扰后,城市轨道交通网络系统性能受损-吸收干扰-系统性能恢复的过程 城市轨道交通网络在灾害干扰下引发实际客流承载功能、抵抗能力等综合系统运行性能恢复提升的过程 表 3 西安城市轨道交通网络线路
Table 3. Urban rail transit network lines of Xi'an
线路编号 起点 终点 站点数 1 沣河森林公园 纺织城 23 2 韦曲南 北客站 21 3 保税区 鱼化寨 26 4 航天新城 北客站(北广场) 29 5 西安东站 创新港 34 6 西安国际医学中心 西北工业大学 13 9 秦陵西 纺织城 15 14(1期) 贺韶 北客站(北广场) 9 表 4 西安城市轨道交通站点节点重要度分布
Table 4. Importance distribution of nodes at urban rail transit stations in Xi'an
节点编号 节点介数 节点度 重要度 节点编号 节点介数 节点度 重要度 0 0.000 0 0.25 0.500 0 78 0.182 6 0.50 0.889 7 1 0.030 9 0.50 0.738 0 79 0.184 3 0.50 0.891 4 2 0.061 4 0.50 0.768 5 80 0.151 1 0.50 0.858 2 3 0.091 5 0.50 0.798 6 81 0.126 4 0.50 0.833 5 4 0.121 2 0.50 0.828 3 82 0.103 6 0.50 0.810 7 5 0.150 5 0.50 0.857 6 83 0.082 2 0.50 0.789 3 6 0.179 4 0.50 0.886 5 84 0.061 7 0.50 0.768 8 7 0.207 9 0.50 0.915 0 85 0.041 7 0.50 0.748 8 8 0.236 0 0.50 0.943 1 86 0.028 0 0.50 0.735 1 9 0.263 7 0.50 0.970 8 87 0.033 3 0.50 0.740 4 10 0.291 0 0.50 0.998 1 88 0.191 5 0.50 0.898 6 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ 68 0.0309 0.50 0.738 0 146 0.317 9 0.50 1.025 0 69 0.061 4 0.50 0.768 5 147 0.344 4 0.50 1.051 5 70 0.091 5 0.50 0.798 6 148 0.370 5 0.50 1.077 6 71 0.121 2 0.50 0.828 3 149 0.000 0 0.25 0.500 0 72 0.150 5 0.50 0.857 6 150 0.030 9 0.50 0.738 0 73 0.179 4 0.50 0.886 5 151 0.061 4 0.50 0.768 5 74 0.207 9 0.50 0.915 0 152 0.073 4 0.50 0.780 5 75 0.236 0 0.50 0.943 1 153 0.068 4 0.50 0.775 5 76 0.087 4 0.50 0.794 5 154 0.074 9 0.50 0.782 0 77 0.211 4 1.00 1.211 4 155 0.085 4 0.50 0.792 5 表 5 Pearson相关系数结果
Table 5. Results of Pearson correlation coefficient
数据 实际新增确诊病例 模型模拟新增确诊病例 实际新增确诊病例 1.000 0.933 表 6 回归分析结果
Table 6. Regression analysis result
R2 F 自由度1 自由度2 显著性 0.870 294.935 1 44 0.000 表 7 ε取值与情景设置
Table 7. ε values and scenarios setting
情景 ε/% 描述 1 100 未采取客流限制措施 2 90 实际承载客流不超过最大可承载客流的90% 3 80 实际承载客流不超过最大可承载客流的80% 4 70 实际承载客流不超过最大可承载客流的70% 5 60 实际承载客流不超过最大可承载客流的60% 6 50 实际承载客流不超过最大可承载客流的50% 7 40 实际承载客流不超过最大可承载客流的40% 8 30 实际承载客流不超过最大可承载客流的30% 9 20 实际承载客流不超过最大可承载客流的20% 表 8 不同参数k下集成韧性模型效果分析
Table 8. Analysis of effect of integrated resilience model with different parameter k
k ε/% C C改善效率 R2 k=1 100 0.999 669 0.976 90 0.967 924 -0.031 756 80 0.928 448 -0.040 785 70 0.937 485 0.009 734 60 0.905 585 -0.034 028 50 0.870 427 -0.038 823 40 0.929 352 0.067 697 30 0.928 375 -0.001 051 20 1.213 994 0.307 655 k=2 100 0.999 669 0.987 90 1.006 206 0.006 539 80 0.999 405 -0.006 759 70 1.106 959 0.107 618 60 1.130 550 0.021 312 50 1.153 512 0.020 311 40 1.468 262 0.272 863 30 1.658 506 0.129 571 20 3.267 084 0.969 896 k=3 100 0.999 669 0.983 90 1.046 002 0.046 348 80 1.075 786 0.028 474 70 1.307 070 0.214 990 60 1.411 402 0.079 821 50 1.528 664 0.083 082 40 2.319 675 0.517 452 30 2.962 857 0.277 273 20 8.792 335 1.967 519 k=4 100 0.999 669 0.961 90 1.087 371 0.087 731 80 1.158 004 0.064 957 70 1.543 355 0.332 772 60 1.762 022 0.141 683 50 2.025 826 0.149 716 40 3.664 801 0.809 041 30 5.293 030 0.444 288 20 23.661 816 3.470 372 表 9 不同ε下城市轨道交通网络集成韧性水平
Table 9. Integrated resilience levels of urban rail transit network under different ε
ε/% V T C累积改善效果 100 0.000 000 0.999 669 90 0.038 788 0.931 099 0.006 539 80 0.073 646 0.862 528 -0.000 264 70 0.166 171 0.793 957 0.107 325 60 0.221 879 0.725 384 0.130 924 50 0.281 583 0.656 814 0.153 894 40 0.457 347 0.588 243 0.468 748 30 0.580 237 0.519 673 0.659 055 20 0.989 982 0.451 100 2.268 166 -
[1] 赵瑞东, 方创琳, 刘海猛. 城市韧性研究进展与展望[J]. 地理科学进展, 2020, 39(10): 1717-1731. doi: 10.18306/dlkxjz.2020.10.011ZHAO Rui-dong, FANG Chuang-lin, LIU Hai-meng. Progress and prospect of urban resilience research[J]. Progress in Geography, 2020, 39(10): 1717-1731. (in Chinese) doi: 10.18306/dlkxjz.2020.10.011 [2] 杨秀平, 王里克, 李亚兵, 等. 韧性城市研究综述与展望[J]. 地理与地理信息科学, 2021, 37(6): 78-84. https://www.cnki.com.cn/Article/CJFDTOTAL-DLGT202106012.htmYANG Xiu-ping, WANG Li-ke, LI Ya-bing, et al. Review and prospects of resilient city theory[J]. Geography and Geo-Information Science, 2021, 37(6): 78-84. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DLGT202106012.htm [3] TWUMASI-BOAKYE R, SOBANJO J O. Resilience of regional transportation networks subjected to hazard-induced bridge damages[J]. Journal of Transportation Engineering, Part A: Systems, 2018, 144(10): 04018062. doi: 10.1061/JTEPBS.0000186 [4] 刘振国, 姜彩良, 王显光, 等. 基于系统韧性提升交通运输疫情防控与应急保障能力对策[J]. 交通运输研究, 2020, 6(1): 19-23. https://www.cnki.com.cn/Article/CJFDTOTAL-JTBH202001004.htmLIU Zhen-guo, JIANG Cai-liang, WANG Xian-guang, et al. Countermeasures to improve epidemic prevention and control and emergency support capacity of transportation based on system resilience[J]. Transport Research, 2020, 6(1): 19-23. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JTBH202001004.htm [5] MURRAY-TUITE P M. A comparison of transportation network resilience under simulated system optimum and user equilibrium conditions[C]//IEEE. Proceedings of the 2006 Winter Simulation Conference. New York: IEEE, 2006: 1398-1405. [6] ZHANG X G, MAHADEVAN S, SANKARARAMAN S, et al. Resilience-based network design under uncertainty[J]. Reliability Engineering and System Safety, 2018, 169: 364-379. doi: 10.1016/j.ress.2017.09.009 [7] LIAO T Y, HU T Y, KO Y N. A resilience optimization model for transportation networks under disasters[J]. Natural Hazards, 2018, 93(1): 469-489. doi: 10.1007/s11069-018-3310-3 [8] 薛锋, 何传磊, 黄倩. 成都地铁网络的关键节点识别方法及性能分析[J]. 中国安全科学学报, 2019, 29(1): 93-99. https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201901016.htmXUE Feng, HE Chuan-lei, HUANG Qian. Identification of key nodes in Chengdu metro network and analysis of network performance[J]. China Safety Science Journal, 2019, 29(1): 93-99. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201901016.htm [9] 沈犁, 张殿业, 向阳, 等. 城市地铁-公交复合网络抗毁性与级联失效仿真[J]. 西南交通大学学报, 2018, 53(1): 156-163, 196. doi: 10.3969/j.issn.0258-2724.2018.01.019SHEN Li, ZHANG Dian-ye, XIANG Yang, et al. Simulation on survivability and cascading failure propagation of urban subway-bus compound network[J]. Journal of Southwest Jiaotong University, 2018, 53(1): 156-163, 196. (in Chinese) doi: 10.3969/j.issn.0258-2724.2018.01.019 [10] MA F, LIU F, YUEN K F, et al. Cascading failures and vulnerability evolution in bus-metro complex bilayer networks under rainstorm weather conditions[J]. International Journal of Environmental Research and Public Health, 2019, 16(3): 329. doi: 10.3390/ijerph16030329 [11] 张琳, 陆建, 雷达. 基于复杂网络和空间信息嵌入的常规公交-地铁复合网络脆弱性分析[J]. 东南大学学报(自然科学版), 2019, 49(4): 773-780. https://www.cnki.com.cn/Article/CJFDTOTAL-DNDX201904022.htmZHANG Lin, LU Jian, LEI Da. Vulnerability analysis of bus-metro composite network based on complex network and spatial information embedding[J]. Journal of Southeast University (Natural Science Edition), 2019, 49(4): 773-780. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DNDX201904022.htm [12] 马超群, 张爽, 陈权, 等. 客流特征视角下的轨道交通网络特征及其脆弱性[J]. 交通运输工程学报, 2020, 20(5): 208-216. doi: 10.19818/j.cnki.1671-1637.2020.05.017MA Chao-qun, ZHANG Shuang, CHEN Quan, et al. Characteristics and vulnerability of rail transit network based on perspective of passenger flow characteristics[J]. Journal of Traffic and Transportation Engineering, 2020, 20(5): 208-216. (in Chinese) doi: 10.19818/j.cnki.1671-1637.2020.05.017 [13] 王志如, 李启明, 梁作论. 城市地铁网络拓扑结构脆弱性评价[J]. 中国安全科学学报, 2013, 23(8): 114-119. doi: 10.3969/j.issn.1003-3033.2013.08.019WANG Zhi-ru, LI Qi-ming, LIANG Zuo-lun. Evaluation of urban metro network topological structure vulnerability[J]. China Safety Science Journal, 2013, 23(8): 114-119. (in Chinese) doi: 10.3969/j.issn.1003-3033.2013.08.019 [14] D'LIMA M, MEDDA F. A new measure of resilience: an application to the London underground[J]. Transportation Research Part A: Policy and Practice, 2015, 81: 35-46. doi: 10.1016/j.tra.2015.05.017 [15] ZHANG X, MILLER-HOOKS E, DENNY K. Assessing the role of network topology in transportation network resilience[J]. Journal of Transport Geography, 2015, 46: 35-45. doi: 10.1016/j.jtrangeo.2015.05.006 [16] ZHANG Dong-ming, DU Fei, HUANG Hong-wei, et al. Resiliency assessment of urban rail transit networks: Shanghai metro as an example[J]. Safety Science, 2018, 106: 230-243. doi: 10.1016/j.ssci.2018.03.023 [17] 王秋玲, 柯宇昊, 高诣民, 等. 灾害时空属性下高铁时序网络韧性分析[J]. 西安电子科技大学学报, 2022, 49(4): 156-166. https://www.cnki.com.cn/Article/CJFDTOTAL-XDKD202204018.htmWANG Qiu-ling, KE Yu-hao, GAO Yi-min, et al. Analysis of resilience of the high-speed rail temporal network under disaster temporal and spatial attributes[J]. Journal of Xidian University, 2022, 49(4): 156-166. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XDKD202204018.htm [18] 唐少虎, 朱伟, 程光, 等. 暴雨内涝下城市道路交通系统安全韧性评估[J]. 中国安全科学学报, 2022, 32(7): 143-150. https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK202207021.htmTANG Shao-hu, ZHU Wei, CHENG Guang, et al. Safety resilience assessment of urban road traffic system under rainstorm waterlogging[J]. China Safety Science Journal, 2022, 32(7): 143-150. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK202207021.htm [19] 李浩然, 王子恒, 杨起帆, 等. 复杂网络下地铁灾害链演化模型与风险分析[J]. 中国安全科学学报, 2021, 31(11): 141-147. https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK202111021.htmLI Hao-ran, WANG Zi-heng, YANG Qi-fan, et al. Evolutionary model and risk analysis of metro disaster chain under complex network[J]. China Safety Science Journal, 2021, 31(11): 141-147. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK202111021.htm [20] 光志瑞, 权经超, 党彤彤, 等. 常态化疫情防控下节假日客流特征及策略研究[J]. 都市快轨交通, 2022, 35(5): 56-61. https://www.cnki.com.cn/Article/CJFDTOTAL-DSKG202205009.htmGUANG Zhi-rui, QUAN Jing-chao, DANG Tong-tong, et al. Characteristics and optimization of holiday passenger flow in urban rail transit subject to coronavirus prevention and control measures[J]. Urban Rapid Rail Transit, 2022, 35(5): 56-61. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DSKG202205009.htm [21] 崔巍, 郝立宁, 木仁, 等. COVID-19疫情下一种新的地铁平峰运营策略[J]. 中国管理科学, 2022, 30(11): 321-332. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGK202211030.htmCUI Wei, HAO Li-ning, MU Ren, et al. A new balance strategy of urban subway operation under the COVID-19 epidemic[J]. Chinese Journal of Management Science, 2022, 30(11): 321-332. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGK202211030.htm [22] 刘斌, 丁波, 赵萌萌, 等. 新型冠状病毒肺炎疫情防控期间武汉地铁客运管理措施分析[J]. 城市轨道交通研究, 2020, 23(10): 5-9. https://www.cnki.com.cn/Article/CJFDTOTAL-GDJT202010003.htmLIU Bin, DING Bo, ZHAO Meng-meng, et al. Analysis of Wuhan metro passenger transport management measures during the period of COVID-19 prevention and control[J]. Urban Mass Transit, 2020, 23(10): 5-9. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GDJT202010003.htm [23] 何坚, 姜美利, 刘建荣. 疫情对地铁车站拥挤度评价的影响[J]. 南方建筑, 2021(2): 146-150. https://www.cnki.com.cn/Article/CJFDTOTAL-NFJZ202102022.htmHE Jian, JIANG Mei-li, LIU Jian-rong. The impacts of COVID-19 on congestion evaluation of metro station[J]. South Architecture, 2021(2): 146-150. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-NFJZ202102022.htm [24] NAN C, SANSAVINI G. A quantitative method for assessing resilience of interdependent infrastructures[J]. Reliability Engineering and System Safety, 2017, 157: 35-53. [25] 李瑞奇, 黄弘, 周睿. 基于韧性曲线的城市安全韧性建模[J]. 清华大学学报(自然科学版), 2020, 60(1): 1-8. https://www.cnki.com.cn/Article/CJFDTOTAL-QHXB202001001.htmLI Rui-qi, HUANG Hong, ZHOU Rui. Resilience curve modelling of urban safety resilience[J]. Journal of Tsinghua University (Science and Technology), 2020, 60(1): 1-8. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-QHXB202001001.htm [26] 黄莺, 刘梦茹, 魏晋果, 等. 基于韧性曲线的城市地铁网络恢复策略研究[J]. 灾害学, 2021, 36(1): 32-36. https://www.cnki.com.cn/Article/CJFDTOTAL-ZHXU202101007.htmHUANG Ying, LIU Meng-ru, WEI Jin-guo, et al. Research on urban metro network recovery strategy based on resilience curve[J]. Journal of Catastrophology, 2021, 36(1): 32-36. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZHXU202101007.htm [27] 吕彪, 高自强, 刘一骝. 道路交通系统韧性及路段重要度评估[J]. 交通运输系统工程与信息, 2020, 20(2): 114-121. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXT202002018.htmLYU Biao, GAO Zi-qiang, LIU Yi-liu. Evaluation of road transportation system resilience and link importance[J]. Journal of Transportation Systems Engineering and Information Technology, 2020, 20(2): 114-121. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSXT202002018.htm [28] 贾兴利, 周吴啸, 韩兴家, 等. 交通管控措施对市域新冠疫情传播的阻断效果分析[J]. 中国公路学报, 2022, 35(1): 252-262. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL202201022.htmJIA Xing-li, ZHOU Wu-xiao, HAN Xing-jia, et al. Blocking effects of traffic control measures on COVID-19 transmission in city territories[J]. China Journal of Highway and Transport, 2022, 35(1): 252-262. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL202201022.htm [29] 段佳勇, 郑宏达. 基于节点重要度的复杂网络脆弱性分析方法[J]. 控制工程, 2020, 27(4): 692-696. https://www.cnki.com.cn/Article/CJFDTOTAL-JZDF202004017.htmDUAN Jia-yong, ZHENG Hong-da. Vulnerability analysis method for complex networks based on node importance[J]. Control Engineering of China, 2020, 27(4): 692-696. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZDF202004017.htm [30] 耿辉, 徐安定, 王晓艳, 等. 基于SEIR模型分析相关干预措施在新型冠状病毒肺炎疫情中的作用[J]. 暨南大学学报(自然科学与医学版), 2020, 41(2): 175-180. https://www.cnki.com.cn/Article/CJFDTOTAL-JNDX202002012.htmGENG Hui, XU An-ding, WANG Xiao-yan, et al. Analysis of the role of current prevention and control measures in the epidemic of corona virus disease 2019 based on SEIR model[J]. Journal of Jinan University (Natural Science and Medicine Edition), 2020, 41(2): 175-180. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JNDX202002012.htm [31] TANG B, WANG X, LI Q, et al. Estimation of the transmission risk of 2019-nCoV and its implication for public health interventions[J]. Journal of Clinical Medicine, 2020, 9(2): 462. [32] 徐娟年, 马冉, 赵中睿, 等. 山东省新型冠状病毒肺炎疫情流行趋势初步预测[J]. 重庆师范大学学报(自然科学版), 2020, 37(2): 101-106. https://www.cnki.com.cn/Article/CJFDTOTAL-CQSF202002016.htmXU Juan-nian, MA Ran, ZHAO Zhong-rui, et al. Preliminary prediction of the epidemic tendency of 2019 novel coronavirus pneumonia in Shandong province[J]. Journal of Chongqing Normal University(Natural Science Edition), 2020, 37(2): 101-106. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CQSF202002016.htm