考虑防火门易损性的震后火灾模拟

王海东 袁佳欣

摘   要：探討一种基于BIM技术的震后次生火灾模拟技术，定量评估防火门损伤对室内火灾温度场分布的影响. 通过选取3条场地类型及地震设计分组一致的地震波对9层Benchmark钢结构模型进行非线性结构分析，并以层间位移角作为非结构构件门损伤的评估指标，结合FEMA P-58中门易损性曲线建立防火门损伤模型，进而基于BIM且利用得到的防火门损伤模型快速准确地建立建筑火源信息模型，导入Pyrosim模拟分析其对火灾温度分布的影响，并得到基于燃烧物与着火点位置关联下的真实温度场分布. 研究结果表明：防火门损伤数量的增加直接影响框架柱的温度，且数量越多框架柱温度越高;除位于着火点及未燃烧的防火分区外，位于其他防火分区的框架柱考虑门损伤的温度明显高于未考虑门损伤的情况，最大达483.98 ℃;且着火点位于不同层数时，各层达到温度边界值的框架柱数量存在差异.

关键词：震后火灾;建筑信息模型;易损性;建筑火源信息模型;温度场分布

中图分类号：TU391;TU243.2                  文献标志码：A

Fire Following Earthquake Simulation Analysis

Considering of Fire Door Fragility

WANG Haidong1，2，YUAN Jiaxin2

（1. Key Laboratory of Building Safety and Energy Efficiency of Ministry of Education（Hunan University），Changsha 410082，China;

2. College of Civil Engineering，Hunan University，Changsha 410082，China）

Abstract：Post-earthquake fire simulation technology based on BIM was discussed to quantitatively evaluate the impact of fire door damage on the fire temperature distribution. Three seismic waves consistent with the site type and its seismic design group were selected to analyze the nonlinear structure of the 9-story Benchmark steel structure model，and the inter-story displacement angle  was used as the evaluation index of the non-structural component door damage，combined with the door vulnerability curve of FEMA P-58 to establish the door damage model. Then，using the fire door damage model obtained by BIM，the information model of the fire was quickly and accurately established. The Pyrosim was imported to simulate and analyze its impact on the fire temperature distribution，and the real temperature field distribution based on the correlation between the burning substance and the fire location was obtained. The results show that the increase in the number of door damages directly affects the temperature of the frame columns，and the larger the number，the higher the temperature of the frame columns; except for the fire and unburned fire compartments，the temperature of frame columns in other fire compartments with regard to door damage is significantly higher than that without considering door damage，with a maximum of 483.98 ℃. When the fire is located at different layers，the number of frame columns that reach the temperature boundary value of each layer is different.

Key words：post-earthquake fire;Building Information Modelling（BIM）;fragility;the information model of the fire;fire temperature distribution

地震过后往往易发生次生灾害，难以控制且损失巨大. 例如：1923年日本关东大地震、1999年台湾集集大地震、2008年汶川地震、2011年日本地震和2014年智利地震等，震后引发大火，火灾造成建筑物损坏甚至倒塌，因而震后火灾结构的性能已成为国内外学者研究的一个热点. Jelinek等[1]对钢框架在震后火灾下的反应与仅受火灾影响时结构的反应进行对比研究;Memari等[2]提出一种基于概率性能的地震火灾分析框架，并利用该框架分析钢结构在地震和火灾连锁危害下的反应，同时研究火灾荷载密度等因素对结构失效概率的影响;Suwondo等[3]研究了不同的火灾位置和地震破坏对荷载重分布路径和结构倒塌的影响;张文[4]通过模糊数学方法建立震后建筑火灾危險性评估模型，研究震后火灾对建筑结构的影响;毛小勇等[5]对轻钢-混凝土组合梁在标准升温下进行了抗火性能的研究;李强等[6]对先后经历地震和火灾作用的破坏形态、损伤机理进行分析，研究损伤指数对圆钢管混凝土柱震后耐火极限的影响. 然而，从已有的研究中发现，研究者对于火灾发生的位置、温度分布，尤其是与火灾温度密切相关的燃烧物数量、分布、防火分区等等都是根据自己的假定而并不与建筑本身的信息相关，缺乏一定的合理性. 与此同时，基于性能的抗震设计开始得到广泛的重视和发展[7]，笔者发现部分学者将建筑信息模型（Building Information Modelling，BIM）技术引入到火灾分析中，具有参考意义.

王婷等[8]将BIM模型导入火灾模拟软件PyroSim中，实现简化火灾建模;陆扬[9]将BIM 技术运用到性能化防火设计中，实现了消防性能相关模拟;杜二峰等[10]设计建造缩尺比为1 ∶ 4的门式刚架厂房模型，测量其模型内部各关键位置的热空气和钢构件温度. 但仅考虑消防性能并未将非结构构件考虑于结构有限元模型中，故由此评估建筑整体防火性能可能缺乏安全性. 且随着人们对建筑使用性的要求不断提高，利用基于FEMA P-58的抗震性能评估方法，开始专注于非结构构件的评估是抗震性能评估的趋势. 而防火门这一非结构构件是建筑在发生火灾时防止火灾蔓延的重要措施，地震发生将导致门有不同程度的破坏，从而影响结构温度场分布，最终影响结构的整体性能，因此通过有限元模型与该非结构构件进行联结，并对建筑整体的防火性能评估应值得关注. Xu等[11]利用BIM和美国抗震设计标准FEMA P-58开展了建筑消防喷淋系统的震害评价，表明其建筑抗火性能与未考虑地震作用时有明显区别. 但是，假定火灾位置进行火灾模拟，未将建筑内燃烧物与火灾位置进行关联，且未明确防火分区，难以得出震后火灾温度场的真实分布，可能缺乏合理性.

本文探讨一种利用BIM技术和FEMA P-58的考虑非结构构件破坏的次生火灾模拟方法，定量评估由于门损伤对火温分布的影响. 该方法基于BIM快速准确地建立建筑火源信息模型，并通过FEMA P-58的抗震性能评估方法评估震后门损伤的情况，进而通过Pyrosim模拟分析防火门破损程度对火灾温度分布的影响，得到基于燃烧物与火灾位置关联下的真实温度场分布，为后续震后火灾结构整体性能研究提供基础.

1 建模分析

1.1   建模思路

考虑震后防火门破损的火灾温度模拟主要包括2种不同的模型：①研究地震作用下非结构构件性能的有限元模型;②反映建筑物实际防火分区和燃烧物信息的火灾模拟模型.

如何将2种独立分割的模型信息统一，本文提出了基于BIM技术建立综合的建筑信息模型的方法，将以上2种模型有效联结，具体思路见图1. 该思路核心为建筑BIM信息模型，即为实现火源信息化的关键，并利用FEMA P-58得出非结构构件的易损性分析，同时结合有限元模型中结构时程分析得出门损伤概率模型，并利用火灾模拟软件FDS模拟震后火灾，分析得出该模型温度的影响，更有利于将地震及震后火灾结合考虑. 其中，有限元模型主要解决地震作用下非结构构件的破损分析问题;FDS火灾模型主要解决火灾模拟分析以及温度实时分析问题;BIM建筑信息模型主要提供结构信息、建筑防火分区、防火门破损位置、状态以及建筑火源分布情况，并贯穿于整个分析流程，为该建模思路的核心模型.

对于体量较大、非结构构件较多且需要建立建筑火源信息的模型，直接建立含有建筑火源信息的火灾数值模型的工作量庞大，故将建筑内火源信息建立于BIM模型中，从而为整体火灾模型的火源信息提供精细化的数据基础，进而通过格式的转换将其导入火灾模拟软件中得到建筑火源信息FDS模型. 将BIM模型转换为FDS模型，就是利用BIM软件将模型导出为fbx文件，再导入至图形化操作软件PyroSim中，软件将正确识别文件格式，并按构件完成分组，此时对模型设置网格数量并对各构件材料属性进行定义，最终运行.

该模型充分实现BIM技术“一模多用”的特点，减少了直接利用PyroSim建立FDS模型的复杂流程，而许镇等[12]，采取直接建模及转换生成FDS模型的方式分别进行火灾模拟，模拟得出同一位置的烟气浓度曲线图，两种方式得出的结果基本一致，证明该转换方式具备较高的准确性.

1.2   模型建立

1.2.1   BIM模型建立

本文将以9层Benchmark模型为例，利用Autodesk Revit2017进行建筑建模，图2为该建筑的三维BIM模型，模型体现建筑功能分区、建筑火源信息、结构全部信息，且各层采用相同的平面布局. 图3为该建筑的平面布置图. 根据美国防火规范NFPA 5000-2012将在该建筑每层设置6个防火分区，图4为各层平面防火分区分布图.

1.2.2   有限元模型创建

本文选用9层Benchmark模型，平面尺寸为45.73 m × 45.73 m，X方向、Y方向对称，各5跨，每跨9.15 m，首层高度为5.49 m，标准层高度为3.96 m，梁、柱单元材料分别采用的是屈服强度为248 MPa和345 MPa的工字钢，表1给出了9层Benchmark钢框架结构模型的构件尺寸，详细参数见文献[13]. 利用SAP2000有限元软件对9层Benchmark三维框架模型进行结构分析，阻尼比ζ = 0.05. 抗震设防烈度为8度（0.20g），场地类别为Ⅱ类（D类），场地分组为第1组.

梁和柱单元塑性铰定义为SAP2000中基于FEMA365中的默认铰属性的集中塑性铰模型. 其中柱单元采用的是考虑弯矩和轴力相互作用的P-M2-M3耦合塑性铰，梁单元采用的是考虑弯矩的M3铰，分别布置在梁柱0.1倍和0.9倍的长度位置处. 图5为该9层Benchmark模型的三维框架立面及平面布置图，图6为利用SAP2000建立的有限元模型.

1.2.3   火灾模拟模型创建

将该9层Benchmark建筑BIM模型导出为fbx文件，并将导出的文件导入至Pyrosim中，补充应有的信息并设置相应的参数可完成模型转换，图7为该9层Benchmark模型转换完成的火灾数值FDS模型.

2   地震后防火门破损评估分析

2.1   易损性数据

本文将利用基于FEMA P-58来考虑震后门损伤对火灾的影响. 在美国规范FEMA P-58中提供了门这一非结构构件的易损性情况的数据，其中包含门的易损性曲线以及两种不同破坏程度的说明. 表2为FEMA P-58中提供的门易损性数据表.

图8为FEMA P-58中提供的门的易损性曲线，将其与地震损伤状态相关联，该曲线中包括两种不同的破坏状态分别为DS1和DS2，其中横坐标为层间位移角数据，纵坐标为对应层间位移角下的损伤概率. DS1指门受到一定程度的变形，导致门的开合受到影响，而当结构变形消失时，门的变形也随之消失，为弹性变形，据研究显示，一般发生于层间位移角为0.001左右时;而DS2指门受到不可自行恢复的损伤且需要修理，据研究显示，一般发生于层间位移角为0.004～0.01时.

2.2   门损伤模型

2.2.1   门损伤概率模型

本文选取9层Benchmark三维框架模型，以最大层间位移角θmax作为结构的损伤判定指标，利用王海东等[14]结合我国抗震规范及美国HAZUS[15]耐震规范中对不同损伤的定义划分. 对应于FEMA P-58提供的门的易损性曲线，由于4种损伤状态下结构均已处于结构非线性阶段，故选用DS2为门损伤状态计算其门损伤概率. 表3为结合门的易损性曲线及各破坏阶段下对应的层间位移角限值，得出在不同损伤状态下建筑内门的整体损伤概率.

2.2.2   门损伤评估

采用非线性评估理论，依照Memari等[16]在研究过程中选取的地震波，并参考我国《建筑抗震设计规范》的标准，选取表4中3条场地类型为Ⅱ类，参照之前学者研究成果[17]，依照震中距及震级得出各地震波的相应设计分组，以便工程应用.

将以上选取的3条地震波按规范要求，做最大值归一化处理，按照抗震设防烈度为8度（0.20g）的罕遇地震情况，取PGA = 0.4g进行调幅. 作用于9层Benchmark三维框架模型中进行非线性结构分析，以各层的层间位移角这一结构评估值作为非结构构件门损伤的评估指数. 表5为最终得出的在各地震记录影响下模型各层层间位移角的平均值及层占比，从而得出建筑模型下，在不同损伤状态时，各层门的损伤情况.

2.2.3   9层Benchmark门损伤模型

该建筑火源信息模型中，门总数为289樘，其中每层防火门为11樘，其他为普通门，各层防火门位置一致，一层防火门位置如图9所示，其中防火门1-1～1-7防火门尺寸均为M1521，1-8～1-11均为M1821. 依照前文所述，将选取的地震记录作用于该建筑的有限元模型中，得到各层平均层间位移角并进行层占比计算. 通过门的易损性曲线对应各状态下门的损伤比例，得出不同状态下各层门损伤数量，从而得出该模型的门损伤模型，表6为地震记录下，结构不同损伤状态时各层门损伤数量.

其中，当结构处于严重破坏及倒塌状态时，门的损伤概率已到达100%，则门损伤数量不依照各层层间位移角平均值的层占比分布，故不考虑. 而结构在中等破坏状态下，考虑最不利布置原则，由于各层门损伤数量均大于11樘（各层防火门数量），则最不利情况为各层防火门均损伤，此时各层火灾均蔓延至各防火分区，均已达最不利状态，故本文不考虑.

本文后续将考虑在结构轻微破坏下建筑的火灾模拟分析，并考虑最不利布置原则，将建筑模型中各层损伤的门均设置为防火门，由此更大程度地扩大室内火灾区域，从而考虑各种火灾发生的最不利情况.

3   火灾模拟分析

本文以結构在轻微破坏状态下为例，利用PyroSim火灾模拟软件，分别将有无门损伤两种情况对火灾的影响进行模拟实验，模拟时间取为500 s. 由于同层门损伤状态下发生火灾，对其他层温度的影响较小，故对各层分别进行模拟实验. 本文选取模型的1层、5层及9层为模拟层数，得出最低、正中间层、最高层各柱处最高温度曲线.

3.1   火灾场景设置

图10中所示的起火点区域位于①防火分区（见图11），该区域在建筑设计中设置了茶水间，内有大功率用电设备，且与各防火分区均有接触面，当防火门损伤时，火灾具有最大的扩散面积，为最不利，且图10表示各层考虑门未损伤时的火灾场景. 图11为各层柱编号所对应框架柱所在位置及防火分区. 同时，由于该建筑功能为办公楼，考虑反应过程、燃烧产物等因素，设定与现实场景相符的起火场景及相应的化学反应. 本文则根据办公室的起火特性，选取办公室标准的Polyurethane GM37（聚氨酯）反应起火. 由于办公楼可能存在大量的沙发坐垫以及可能存放的纸箱，依据最不利原则，设定火灾为快速增长型，火灾增长系数α取0.046 89，火源最大热释放速率为1 000.0 kW，按照t2火模型计算得到火灾到达最大速率的时间为146 s，在本文中，由于设置的燃烧物有限，模拟火灾过程影响，发现各柱在500 s内均已达到火灾温度下降段，故选取该模型模拟时间为500 s.

本文主要研究该模型结构在轻微破坏下，1、5、9层在有无考虑门损伤状态两种情况下对结构温度的影响. 由前文可知，5层门损伤的数量为3，1、9层门损伤数量为4，由此进行模拟计算.

由于1、5、9层的门损伤数量以平面视角分类，则共有2种损伤类型，分别为门损伤数量为3和4. 考虑火灾现场设计的最不利原则，应将损伤的防火门设置在火源物质影响更大的位置处，故选择图12、图13所示的设置方式，门损伤数量分别为3和4.

3.2   温度数据

在以上火灾模拟中，共有6种工况，分别为工况1：9层4防火门损伤;工况2：9层门未损伤;工况3：5层3防火门损伤;工况4：5层门未损伤;工况5：1层4防火门损伤 ;工况6：1层门未损伤. 由于模型中各层框架柱共计36个，模拟500 s内各层柱中部温度变化情况，并记录其最高温度，按各防火分区进行统计，分别得出图14～图19各防火分区6种不同工况的最高温度情况，表7为各防火分区内包含的柱编号.

3.3   模拟结果分析

由各层柱在6种工况、不同门损伤情况下统计得到的最高温度可知，门损伤数量的增多将对框架柱温度有直接影响，且门损伤数量越多，框架柱温度越高. 统计可得，在火灾模拟500 s内，当框架柱最高温度高于300 ℃时，则100 s内基本达到构件的最高温度;当框架柱最高温度低于300 ℃时，则在模拟时间为300 s后基本达到构件的最高温度.

通过统计各工况下各框架柱的最高温度可知：框架柱位于着火点及在所有工况下均未燃烧的防火分区时，考虑地震作用和未考虑地震作用的框架柱燃烧温度基本吻合;而框架柱位于其他燃烧点的防火分区时，未考虑地震作用下门损伤情况的燃烧温度将明显低于考虑地震作用下门损伤情况的燃烧温度. 且以本文1、5、9层对应考虑地震作用与未考虑地震作用的所有工况下，统计可得考虑门损与未考虑门损的框架柱温度差值最大达483.98 ℃，该点为9层23号柱.

根据《建筑钢结构防火技术规范》[18]可知，当温度达到300 ℃以上时钢材的屈服强度将进行折减，故本文以计算屈服强度折减系数为1.0、0.95、0.9时钢材的温度作为统计边界值，分别为300 ℃、375 ℃、409 ℃. 图20为统计各工况下框架柱最高温度达到该三者统计边界值以上的數量，为后续分析温度对柱承载能力的影响提供基础. 且由各框架柱最高温度的统计结果可知，考虑地震作用门损伤时达到边界温度值的框架柱数量比不考虑地震作用门损伤时多达4.25～7.33倍;当考虑门损伤情况下，框架柱温度及达到统计边界值的数量排序均为：9层>1层>5层;当不考虑门损伤的情况下各层达统计边界值数量趋于稳定.

4   结   论

本文利用BIM技术和FEMA P-58提出了一种基于性能的考虑非结构构件破损的震后次生火灾模拟方法，定量评估由于门损伤对火灾温度场的影响，并反映于框架柱的温度，相关结论如下.

1）提出了基于BIM的建筑火源信息模型，并将其作为信息传递中心，将结构模型与火灾模型相结合.

2）将地震作用下的层间位移角与非结构构件门的易损性进行关联：

①当结构处于严重破坏或倒塌状态时，各层门损伤概率达100%，此时震后火灾温度场与防火门无关;

②当结构处于中等破坏时，门损伤概率为56%，此时各门为有概率的损伤，若考虑最不利情况为防火门全数损伤，火灾蔓延至全部防火分区，将达到火灾最不利影响;

③在结构处于轻微破坏时，门损伤的概率为10%，防火门是否损伤具有随机性，考虑地震作用下防火门损伤情况对震后火灾温度场分布具有显著的影响，故在中低地震设防区具有重要意义.

3）当结构处于轻微破坏时，通过对比6种工况下各框架柱的温度变化情况及最高温度统计可得：

①当框架柱位于非着火点及引燃的防火分区时，未考虑地震作用门损伤的燃烧温度将明显低于考虑地震作用门损伤的燃烧温度，且框架柱温度差值最大为483.98 ℃;

②通过统计钢材屈服强度折减系数为1.0、0.95、0.9所对应的边界温度值300 ℃、375 ℃、409 ℃的框架柱数量可知，考虑门损时达到边界温度值的框架柱数量将远远高于未考虑地震作用门损时的框架柱数量，其比值最高可达7.33倍;

③达统计边界值的数量具有差异但有一定规律性，在考虑地震作用门损伤时，各层达统计边界值数量为9层>1层>5层.未考虑地震作用门损伤时，各层达统计边界值数量趋于稳定.

参考文献

[1]    JELINEK T，ZANIA V，GIULIANI L. Post-earthquake fire resistance of steel buildings[J]. Journal of Constructional Steel Research，2017，138（8）：774—782.

[2]    MEMARI M，MAHMOUD H. Framework for a performance-based analysis of fires following earthquakes[J]. Engineering Structures，2018，171（5）：794—805.

[3]    SUWONDO R，CUNNINGHAM L，GILLIE M. Progressive collapse analysis of composite steel frames subject to fire following earthquake[J]. Fire Safety Journal，2019，103（1）：49—58.

[4]    张文. 基于模糊数学的震后建筑火灾危险性评估模型研究[J].地震工程学报，2018，40（6）：1372—1377.

ZHANG W. A model of fire risk assessment of post-earthquake buildings based on fuzzy mathematics[J]. China Earthquake Engineering Journal，2018，40（6）：1372—1377. （In Chinese）

[5]    毛小勇，肖岩. 标准升温下轻钢-混凝土组合梁的抗火性能研究[J]. 湖南大学学报（自然科学版），2005，32（2）：64—70.

MAO X Y，XIAO Y. Behavior of lightweight steel-concrete composite beams subjected to standard fire[J]. Journal of Hunan University（Natural Sciences），2005，32（2）：64—70. （In Chinese）

[6]    李强，王文达，毛文婧. 震损后圆钢管混凝土柱耐火性能有限元分析[J]. 建筑结构学报，2017，38（S1）：118—125.

LI Q，WANG W D，MAO W J. Finite element analysis of concrete-filled circular steel tubular columns subjected to post-earthquake fire[J]. Journal of Building Structures，2017，38（S1）：118—125. （In Chinese）

[7]    陆新征，曾翔，许镇. 建设地震韧性城市所面临的挑战[J]. 城市与减灾，2017（4）：29—34.

LU X Z，ZENG X，XU Z. The challenge of building a seismically resilient city[J]. City and Disaster Reduction，2017（4）：29—34. （In Chinese）

[8]    王婷，杜慕皓，唐永福，等. 基于BIM的火灾模拟与安全疏散分析[J]. 土木建筑工程信息技术，2014，6（6）：102—108.

WANG T，DU M H，TANG Y F，et al. Analysis on fire model and safety evacuation based on BIM[J]. Journal of Information Technology in Civil Engineering and Architecture，2014，6（6）：102—108. （In Chinese）

[9]    陆扬. 基于BIM的性能化分析手段在建筑防火设计中的研究与实践[J]. 土木建筑工程信息技术，2011，3（4）：63—71.

LU Y. BIM-based performance analysis tools in building fire safety design of the research and practice[J]. Journal of Information Technology in Civil Engineering and Architecture，2011，3（4）：63—71. （In Chinese）

[10]  杜二峰，舒贛平，吕晓. 1 ∶ 4缩尺厂房中功率火灾试验研究[J]. 湖南大学学报（自然科学版），2020，47（3）：82—88.

DU E F，SHU G P，L？譈 X.Medium power fire test study on a 1 ∶ 4 scaled portal frame building[J]. Journal of Hunan University （Natural Sciences），2020，47（3）：82—88. （In Chinese）

[11]  XU Z，ZHANG Z C，LU X Z，et al. Post-earthquake fire simulation considering overall seismic damage of sprinkler systems based on BIM and FEMA P-58[J]. Automation in Construction，2018，90（2）：9—22.

[12]  许镇，张宗才，陆新征，等. 基于BIM和FEMAP-58的建筑地震次生火灾模拟[J]. 消防科学与技术，2018，37（6）：745—749.

XU Z，ZHANG Z C，LU X Z，et al. Building post-earthquake fire simulation based on BIM and FEMA P-58[J]. Fire Science and Technology，2018，37（6）：745—749. （In Chinese）

[13]  OHTORI Y，CHRISTENSON R E，SPENCER B F J，et al. Benchmark control problems for seismically excited nonlinear buildings[J]. Journal of Engineering Mechanics，2004，130（4）：366—385.

[14]  王海東，武战战，刘巴黎. 地震动输入对中长周期结构抗倒塌性能评估的影响[J]. 地震工程与工程振动，2019，39（3）：176—184.

WANG H D，WU Z Z，LIU B L. Influence of ground motion input on evaluation of anti-collapse performance of medium-long period structure[J].Earthquake Engineering and Engineering Dynamics，2019，39（3）：176—184. （In Chinese）

[15]  Users manual：HAZUS99[S]. Washington D C：Federal Emergency Management Agency，l999：236—268.

[16]  MEMARI M，MAHMOUD H，ELLINGWOOD B. Post-earthquake fire performance of moment resisting frames with reduced beam section connections[J]. Journal of Constructional Steel Research，2014，103（9）：215—229.

[17]  冀昆，温瑞智，任叶飞. 中国抗震规范强震记录选取的初选条件研究[J]. 防灾减灾工程学报，2016，36（1）：44—49.

JI K，WEN R Z，REN Y F. Study on the initial of selection criteria strong motion records for Chinese seismic code[J]. Journal of Disaster Prevention and Mitigation Engineering，2016，36（1）：44—49. （In Chinese）

[18]  建筑钢结构防火技术规范：GB 51249—2017[S]. 北京：中国计划出版社，2017：19—20.

Code for fire safety of steel structures in buildings：GB 51249—2017[S]. Beijing：China Planning Press，2017：19—20. （In Chinese）