CN115601520A - 一种基于二维预估和多尺度网格的火灾模拟方法及装置 - Google Patents
一种基于二维预估和多尺度网格的火灾模拟方法及装置 Download PDFInfo
- Publication number
- CN115601520A CN115601520A CN202211513592.8A CN202211513592A CN115601520A CN 115601520 A CN115601520 A CN 115601520A CN 202211513592 A CN202211513592 A CN 202211513592A CN 115601520 A CN115601520 A CN 115601520A
- Authority
- CN
- China
- Prior art keywords
- fire
- dimensional
- building group
- sub
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000007480 spreading Effects 0.000 claims abstract description 57
- 238000003892 spreading Methods 0.000 claims abstract description 57
- 230000007246 mechanism Effects 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims description 37
- 238000004891 communication Methods 0.000 claims description 22
- 239000012530 fluid Substances 0.000 claims description 14
- 230000005855 radiation Effects 0.000 claims description 14
- 238000012546 transfer Methods 0.000 claims description 13
- 238000009826 distribution Methods 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 230000001502 supplementing effect Effects 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 230000002265 prevention Effects 0.000 abstract description 5
- 230000009467 reduction Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 18
- 230000008569 process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000015654 memory Effects 0.000 description 6
- 239000000463 material Substances 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 239000000779 smoke Substances 0.000 description 3
- 239000004566 building material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 206010000369 Accident Diseases 0.000 description 1
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- XFXPMWWXUTWYJX-UHFFFAOYSA-N Cyanide Chemical compound N#[C-] XFXPMWWXUTWYJX-UHFFFAOYSA-N 0.000 description 1
- 208000000059 Dyspnea Diseases 0.000 description 1
- 206010013975 Dyspnoeas Diseases 0.000 description 1
- 206010019233 Headaches Diseases 0.000 description 1
- 206010028813 Nausea Diseases 0.000 description 1
- 208000012902 Nervous system disease Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 229910002091 carbon monoxide Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 231100000869 headache Toxicity 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008693 nausea Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/04—Architectural design, interior design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Mathematical Optimization (AREA)
- Remote Sensing (AREA)
- Computational Mathematics (AREA)
- Structural Engineering (AREA)
- Computer Graphics (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及土木工程防灾减灾技术领域，特别是指一种基于二维预估和多尺度网格的火灾模拟方法及装置，一种基于二维预估和多尺度网格的火灾模拟方法包括：导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型；根据所述建筑群二维模型进行三维建模，获得建筑群三维模型；基于所述建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围；根据所述二维火灾蔓延范围，基于所述建筑群三维模型，获得目标建筑群三维模型；基于所述目标建筑群三维模型，进行建筑群高效三维火灾模拟。本发明能够基于二维预估和可变多尺度网格，结合建筑群三维模型，能够在较短时间内完成建筑群精细化火灾模拟。
Description
技术领域
本发明涉及土木工程防灾减灾技术领域，特别是指一种基于二维预估和多尺度网格的火灾模拟方法及装置。
背景技术
近些年来，国内外发生了多起严重的建筑群火灾事故，造成严重的人员伤亡与巨大的经济损失。建筑中使用的许多可燃材料或高分子材料在燃烧时会释放出一氧化碳、氰化物等有毒烟气，人吸入这些烟气后，会出现呼吸困难、头痛、恶心、神经系统紊乱等症状，严重威胁人民的生命安全；在各地火灾造成的经济损失中，建筑火灾造成的经济损失居首位，建筑火灾产生的高温高热会对建筑结构造成破坏，甚至造成建筑的整体倒塌，建筑火灾发生后，建筑的修复重建、人员安置、生产经营停顿等，会造成巨大的间接经济损失；历史保护建筑和文化遗址一旦发生火灾，不仅会造成人员伤亡和财产损失，还会烧毁大量文物、经典书籍、古建筑等众多稀世珍宝，对人类文明成果造成不可弥补的损失。所以针对建筑群火灾蔓延模拟受到国内外学者高度重视，并已开展了较多研究。
目前，建筑群火灾蔓延模拟模型主要分为经验模型和物理模型。经验模型主要对火灾现场和计算机模拟得出的数据来进行汇总、测量和分析处理，总结经验公式，从而借助计算机的数据库、图形图像等功能，形象地分析建筑物中火灾的发展；物理模型即场模型，物理模型利用计算机，以流体力学为基础，利用连续性方程、动量方程、能量方程等数值求解的方法来得到火灾过程中典型状态参数的空间分布以及计算参数随时间的变化规律。
但在二维的经验模型或简化的物理模型中，难以考虑建筑三维空间特征及风场分布，也难以准确模拟火灾发展过程，可能造成火灾蔓延结果与实际存在较大偏差，因此，在现有火灾模拟技术中心，缺乏一种准确高效的建筑群火灾蔓延模拟方法。
发明内容
本发明实施例提供了一种基于二维预估和多尺度网格的火灾模拟方法及装置。所述技术方案如下：
一方面，提供了一种基于二维预估和多尺度网格的火灾模拟方法，该方法由电子设备实现，该方法包括：
导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型。
根据所述建筑群二维模型进行三维建模，获得建筑群三维模型。
基于所述建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围。
根据所述二维火灾蔓延范围，基于所述建筑群三维模型，获得目标建筑群三维模型。
基于所述目标建筑群三维模型，进行建筑群高效三维火灾模拟。
可选地，所述导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型，包括：
将建筑群卫星影像与几何模型加载到地理信息系统平台中，依据建筑物编号，为各建筑物匹配几何模型，补充层高、层数、建筑面积和结构防火等级的模型参数，进而完成二维模型的建立。
可选地，所述基于所述建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围，包括：
根据预设的建筑物间的火灾蔓延机制以及历史火灾统计数据，设置二维火灾模拟所需参数；其中，所述预设的火灾蔓延机制包括热辐射和热羽流，所述历史火灾统计数据包括起火点、起燃数量和气象条件；
根据预设的火灾蔓延经验公式、所述二维火灾模拟所需参数、以及所述建筑群二维模型，对建筑群的二维火灾蔓延范围进行预估。
可选地，所述基于所述目标建筑群三维模型，进行建筑群高效三维火灾模拟，包括：
S41、根据目标建筑群三维模型建立相应的流体动力学模型；
S42、在流体动力学模型中，通过递归坐标二分法将二维火灾蔓延范围初步划分出各子区域；再根据火源位置的影响范围，在各子区域中划分不同精细度的网格；
S43、计算各子区域的负载均衡指标，根据各子区域的负载均衡指标、基于各子区域的网格负载调整网格边界，计算网格调整后的各子区域的工作负载，直到各子区域的工作负载近似相等；
S44、开始模拟，设定时间步，读取各子区域的计算时间，若各子区域的计算时间分布不均衡，利用负载迭代函数动态调整各子区域的工作负载，直至模拟结束。
可选地，所述负载均衡指标为基于各子区域的计算时间测量值提出的量化指标，其中，所述负载均衡指标包括全域平均计算时间、子区域计算效率和累积不均衡负载。
可选地，所述网格负载为基于网格边界长度和网格数目通过负载描述函数计算得到，其中，所述网格负载包括网格拓扑结构负载和通信负载。
可选地，所述工作负载为计算结构平衡、通信平衡参数的量化指标，其中，所述工作负载包括结构负载平衡因子、通信负载平衡因子、火源辐射影响因子和可燃物含量因子。
可选地，所述利用负载迭代函数动态调整各子区域的工作负载，包括：
根据负载均衡指标中的累积不均衡负载，计算得到转移负载和网格移动方向；
根据所述转移负载和所述网格移动方向，通过贪婪算法进行负载迭代，调整各子区域的网格边界，直至各子区域的工作负载达到近似相等。
另一方面，提供了一种基于二维预估和多尺度网格的火灾模拟装置，该装置应用于一种基于二维预估和多尺度网格的火灾模拟方法，该装置包括：
二维模型构建模块，用于导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型；
三维模型构建模块，用于根据所述建筑群二维模型进行三维建模，获得建筑群三维模型；
范围预估模块，用于基于所述建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围；
目标获取模块，用于根据所述二维火灾蔓延范围，基于所述建筑群三维模型，获得目标建筑群三维模型；
火灾模拟模块，用于基于所述目标建筑群三维模型，进行建筑群高效三维火灾模拟。
可选地，所述二维模型构建模块，进一步用于：
将建筑群卫星影像与几何模型加载到地理信息系统平台中，依据建筑物编号，为各建筑物匹配几何模型，补充层高、层数、建筑面积和结构防火等级的模型参数，进而完成二维模型的建立。
可选地，所述范围预估模块，进一步用于：
根据预设的建筑物间的火灾蔓延机制以及历史火灾统计数据，设置二维火灾模拟所需参数；其中，所述预设的火灾蔓延机制包括热辐射和热羽流，所述历史火灾统计数据包括起火点、起燃数量和气象条件；
根据预设的火灾蔓延经验公式、所述二维火灾模拟所需参数、以及所述建筑群二维模型，对建筑群的二维火灾蔓延范围进行预估。
可选地，所述火灾模拟模块，进一步用于：
S41、根据目标建筑群三维模型建立相应的流体动力学模型；
S42、在流体动力学模型中，通过递归坐标二分法将二维火灾蔓延范围初步划分出各子区域；再根据火源位置的影响范围，在各子区域中划分不同精细度的网格；
S43、计算各子区域的负载均衡指标，根据各子区域的负载均衡指标、基于各子区域的网格负载调整网格边界，计算网格调整后的各子区域的工作负载，直到各子区域的工作负载近似相等；
S44、开始模拟，设定时间步，读取各子区域的计算时间，若各子区域的计算时间分布不均衡，利用负载迭代函数动态调整各子区域的工作负载，直至模拟结束。
可选地，所述火灾模拟模块，进一步用于：
所述负载均衡指标为基于各子区域的计算时间测量值提出的量化指标，其中，所述负载均衡指标包括全域平均计算时间、子区域计算效率和累积不均衡负载。
可选地，所述火灾模拟模块，进一步用于：
所述网格负载为基于网格边界长度和网格数目通过负载描述函数计算得到，其中，所述网格负载包括网格拓扑结构负载和通信负载。
可选地，所述火灾模拟模块，进一步用于：
所述工作负载为计算结构平衡、通信平衡参数的量化指标，其中，所述工作负载包括结构负载平衡因子、通信负载平衡因子、火源辐射影响因子和可燃物含量因子。
可选地，所述火灾模拟模块，进一步用于：
所述利用负载迭代函数动态调整各子区域的工作负载，包括：
根据负载均衡指标中的累积不均衡负载，计算得到转移负载和网格移动方向；
根据所述转移负载和所述网格移动方向，通过贪婪算法进行负载迭代，调整各子区域的网格边界，直至各子区域的工作负载达到近似相等。
另一方面，提供了一种电子设备，所述电子设备包括处理器和存储器，所述存储器中存储有至少一条指令，所述至少一条指令由所述处理器加载并执行以实现上述一种基于二维预估和多尺度网格的火灾模拟方法。
另一方面，提供了一种计算机可读存储介质，所述存储介质中存储有至少一条指令，所述至少一条指令由处理器加载并执行以实现上述一种基于二维预估和多尺度网格的火灾模拟方法。
本发明实施例提供的技术方案带来的有益效果至少包括：
本发明能够基于二维预估和可变多尺度网格，结合建筑群三维模型，能够在较短时间内完成建筑群精细化火灾模拟，可用于城市建筑群火灾蔓延过程的高效模拟，为城市火灾风险防控与消防规划提供准确高校的科学依据。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种基于二维预估和多尺度网格的火灾模拟方法流程图；
图2是本发明实施例提供的一种负载迭代函数的建筑群三维火灾模拟流程图；
图3是本发明实施例提供的一种负载计算示意图；
图4是本发明实施例提供的一种考虑负载的可变多尺度网格划分示意图；
图5是本发明实施例提供的一种基于二维预估和多尺度网格的火灾模拟装置框图；
图6是本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
本发明实施例提供了一种基于二维预估和多尺度网格的火灾模拟方法，该方法可以由电子设备实现，该电子设备可以是终端或服务器。如图1所示的一种基于二维预估和多尺度网格的火灾模拟方法流程图，该方法的处理流程可以包括如下的步骤：
S1、导入建筑群数据，根据建筑群数据，获得建筑群二维模型。
可选地，导入建筑群数据，根据建筑群数据，获得建筑群二维模型，包括：
将建筑群卫星影像与几何模型加载到地理信息系统平台中，依据建筑物编号，为各建筑物匹配几何模型，补充层高、层数、建筑面积和结构防火等级的模型参数，进而完成二维模型的建立。
一种可行的实施方式中，通过卫星采集建筑群的轮廓外形数据，根据建筑群轮廓外形分类简化为形容建筑群的多边形数据，再通过建筑群的多边形数据传输至地理信息系统平台，对卫星采集到的建筑群进行数据匹配，确定建筑群中建筑物的具体建筑高度、建筑采用的结构类型和建筑固有面积等与火灾相关的信息。
S2、根据建筑群二维模型进行三维建模，获得建筑群三维模型。
一种可行的实施方式中，根据建筑群的几何模型、建筑物楼层数和建筑物高度等火灾蔓延模拟相关数据的具体信息，在建筑群二维模型的基础上进行拉伸，建立建筑群三维模型。根据建筑群具体数据基于二维模型，建立建筑群三维模型是本领域的常用手段，此处不做赘述。
S3、基于建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围。
可选地，基于建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围，包括：
根据预设的建筑物间的火灾蔓延机制以及历史火灾统计数据，设置二维火灾模拟所需参数；其中，预设的火灾蔓延机制包括热辐射和热羽流，历史火灾统计数据包括起火点、起燃数量和气象条件；
根据预设的火灾蔓延经验公式、二维火灾模拟所需参数、以及建筑群二维模型，对建筑群的二维火灾蔓延范围进行预估。
一种可行的实施方式中，预设的火灾蔓延机制是指建筑物间的火灾蔓延，是通过热传播进行的，其形式与起火点、建筑材料、物质的燃烧性能和可燃物的数量等因素有关。在火场上燃烧物质所放出的热能，通常是以传导、辐射和对流三种方式传播，并影响火势蔓延扩大。预设的火灾蔓延经验公式是指根据起火点、建筑材料等关于火灾蔓延的相关数据计算火灾蔓延的范围。预设的火灾蔓延机制和预设的火灾蔓延经验公式是本领域常用技术手段，此处不做赘述。
S4、根据二维火灾蔓延范围，基于建筑群三维模型，获得目标建筑群三维模型。
一种可行的实施方式中，根据上述步骤得到在建筑群在二维平面上的建筑分布，根据上述步骤计算得到二维平面上的火灾蔓延范围，在建筑群三维模型中，选取二维火灾蔓延范围内的建筑群模型为目标建筑群三维模型，选取目标建筑群三维模型可以有效减少待模拟的建筑数量节省了模拟和计算时间，使计算计算响应更敏捷。
S5、基于目标建筑群三维模型，进行建筑群高效三维火灾模拟。
可选地，基于目标建筑群三维模型，进行建筑群高效三维火灾模拟，包括：
S51、根据目标建筑群三维模型建立相应的流体动力学模型；
一种可行的实施方式中，计算模拟建筑群火灾中的烟气和热传递过程，准确预测火灾发展的全过程，所以需要采用流体力学模型的环境下下进行模拟。
S52、在流体动力学模型中，通过递归坐标二分法将二维火灾蔓延范围初步划分出各子区域；再根据火源位置的影响范围，在各子区域中划分不同精细度的网格；
一种可行的实施方式中，在递归坐标二分法的基础上，按火源位置的影响范围，将计算域分为直接覆盖区，间接覆盖区以及辐射区，网格精度逐渐增加，其判别公式如下式（1）所示：
其中，d为网格边界到火源的距离，hr为火源的热释放速率。
S53、计算各子区域的负载均衡指标，根据各子区域的负载均衡指标、基于各子区域的网格负载调整网格边界，计算网格调整后的各子区域的工作负载，直到各子区域的工作负载近似相等；
可选地，负载均衡指标为基于各个子区域的计算时间测量值提出的量化指标，其中，负载均衡指标包括全域平均计算时间、子区域计算效率和累积不均衡负载。
一种可行的实施方式中，负载均衡指标中全域平均计算时间、子区域计算效率和累积不均衡负载的计算公式如下式（2）、（3）、（4）所示：
全域平均计算时间：
子区域计算效率：
累积不均衡负载：
其中，N为子区域数量；
为各子区域计算时间。
可选地，网格负载为基于网格边界长度和网格数目通过负载描述函数计算得到，其中，网格负载包括网格拓扑结构负载和通信负载。
一种可行的实施方式中，综合考虑包括但不限于网格数量、通信流量、可燃物等因素，提出负载描述函数，对初始状态的网格负载进行计算描述，负载描述函数中关于网格拓扑结构负载和通信负载计算公式如下式（5）、（6）所示：
网格拓扑结构负载：
通信负载：
其中，
为子区域的边界长度；
为网格精度；M为子区域内
的网格总数；T为子区域所有边界面上的网格总数，也是和邻近子区域通信负载的上限，只
有在所有边界面均为内边界时才取此值；f为映射函数，由选择的数值模拟求解器确定；μ为
常数，由选择的并行流场求解器决定。
可选地，工作负载为计算结构平衡、通信平衡参数的量化指标，其中，工作负载包括结构负载平衡因子、通信负载平衡因子、火源辐射影响因子和可燃物含量因子。
一种可行的实施方式中，调整网格边界将网格边界由网格负载大的区域向网格负载小的区域移动，该操作的本质是调整网格的疏密程度，使进行该操作后的在目标建筑群的范围内，工作负载在各个网格中达到均匀分布的状态。计算网格调整后的工作负载，通过调整网格边界使各子区域的工作负载达到近似相等。各子区域工作负载近似相等是指各子区域的工作负载误差不超过±1%。
工作负载中的结构负载平衡因子、通信负载平衡因子、火源辐射影响因子和可燃物含量因子的计算公式如下式（7）、（8）、（9）、（10）所示：
结构负载平衡因子：
其中，
和
分别为各子区域网格数的最大值与平均值。
通信负载平衡因子：
其中，
和
分别为各子区域通信表面边界网格数的最大值与平均值。
火源辐射影响因子：
其中，
和
分别为各子区域中心到火源距离的最大值与平均值。
可燃物含量因子：
其中，
和
分别为各子区域可燃物的数量与体积平均值。
工作负载计算公式如下式（11）所示：
工作负载：
S54、开始模拟，设定时间步，读取各子区域的计算时间，若各子区域的计算时间分布不均衡，利用负载迭代函数动态调整各子区域的工作负载，直至模拟结束。
可选地，如图2利用负载迭代函数动态调整各子区域的工作负载，包括：
根据负载均衡指标中的累积不均衡负载，计算得到转移负载和网格移动方向；
根据转移负载和网格移动方向，通过贪婪算法进行负载迭代，调整各子区域的网格边界，直至各子区域的工作负载达到近似相等。
一种可行的实施方式中，根据上述步骤计算结果，转移负载、移动方向和负载迭代的计算公式如下式（12）、（13）、（14）所示：
转移负载：
其中，
为子区域的分区负载；
为子区域的总负载。
移动方向：
迭代函数：
其中，
为惩罚因子，防止过度迭代，取1.3；误差
取0.01，当
时，迭代停
止，如图3、图4所示。
本发明能够基于二维预估和可变多尺度网格，结合建筑群三维模型，能够在较短时间内完成建筑群精细化火灾模拟，可用于城市建筑群火灾蔓延过程的高效模拟，为城市火灾风险防控与消防规划提供准确高校的科学依据。
图5是根据一示例性实施例示出的一种基于二维预估和多尺度网格的火灾模拟装置框图，该系统应用于实现一种基于二维预估和多尺度网格的火灾模拟方法。参照图5，该装置包括：
二维模型构建模块510，用于导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型；
三维模型构建模块520，用于根据所述建筑群二维模型进行三维建模，获得建筑群三维模型；
范围预估模块530，用于基于建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围；
目标获取模块540，用于根据二维火灾蔓延范围，基于建筑群三维模型，获得目标建筑群三维模型；
火灾模拟模块550，用于基于目标建筑群三维模型，进行建筑群高效三维火灾模拟。
可选地，二维模型构建模块510，进一步用于：
将建筑群卫星影像与几何模型加载到地理信息系统平台中，依据建筑物编号，为各建筑物匹配几何模型，补充层高、层数、建筑面积和结构防火等级的模型参数，进而完成二维模型的建立。
可选地，范围预估模块530，进一步用于：
根据预设的建筑物间的火灾蔓延机制以及历史火灾统计数据，设置二维火灾模拟所需参数；其中，预设的火灾蔓延机制包括热辐射和热羽流，历史火灾统计数据包括起火点、起燃数量和气象条件；
根据预设的火灾蔓延经验公式、二维火灾模拟所需参数、以及建筑群二维模型，对建筑群的二维火灾蔓延范围进行预估。
可选地，火灾模拟模块550，进一步用于：
S41、根据目标建筑群三维模型建立相应的流体动力学模型；
S42、在流体动力学模型中，通过递归坐标二分法将二维火灾蔓延范围初步划分出各子区域；再根据火源位置的影响范围，在各子区域中划分不同精细度的网格；
S43、计算各子区域的负载均衡指标，根据各子区域的负载均衡指标、基于各子区域的网格负载调整网格边界，计算网格调整后的各子区域的工作负载，直到各子区域的工作负载近似相等；
S44、开始模拟，设定时间步，读取各子区域的计算时间，若各子区域的计算时间分布不均衡，利用负载迭代函数动态调整各子区域的工作负载，直至模拟结束。
可选地，火灾模拟模块550，进一步用于：
负载均衡指标为基于各个子区域的计算时间测量值提出的量化指标，其中，负载均衡指标包括全域平均计算时间、子区域计算效率和累积不均衡负载。
可选地，火灾模拟模块550，进一步用于：
网格负载为基于网格边界长度和网格数目通过负载描述函数计算得到，其中，网格负载包括网格拓扑结构负载和通信负载。
可选地，火灾模拟模块550，进一步用于：
工作负载为计算结构平衡、通信平衡参数的量化指标，其中，工作负载包括结构负载平衡因子、通信负载平衡因子、火源辐射影响因子和可燃物含量因子。
可选地，火灾模拟模块550，进一步用于：
利用负载迭代函数动态调整各子区域的工作负载，包括：
根据负载均衡指标中的累积不均衡负载，计算得到转移负载和网格移动方向；
根据转移负载和网格移动方向，通过贪婪算法进行负载迭代，调整各子区域的网格边界，直至各子区域的工作负载达到近似相等。
本发明能够基于二维预估和可变多尺度网格，结合建筑群三维模型，能够在较短时间内完成建筑群精细化火灾模拟，可用于城市建筑群火灾蔓延过程的高效模拟，为城市火灾风险防控与消防规划提供准确高校的科学依据。
图6是本发明实施例提供的一种电子设备600的结构示意图，该电子设备600可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器（centralprocessing units，CPU）601和一个或一个以上的存储器602，其中，所述存储器602中存储有至少一条指令，所述至少一条指令由所述处理器601加载并执行以实现上述一种基于二维预估和多尺度网格的火灾模拟方法的步骤。
在示例性实施例中，还提供了一种计算机可读存储介质，例如包括指令的存储器，上述指令可由终端中的处理器执行以完成上述一种基于二维预估和多尺度网格的火灾模拟方法。例如，所述计算机可读存储介质可以是ROM、随机存取存储器（RAM）、CD-ROM、磁带、软盘和光数据存储设备等。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述方法包括：
导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型；
根据所述建筑群二维模型进行三维建模，获得建筑群三维模型；
基于所述建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围；
根据所述二维火灾蔓延范围，基于所述建筑群三维模型，获得目标建筑群三维模型；
基于所述目标建筑群三维模型，进行建筑群高效三维火灾模拟。
2.根据权利要求1所述的一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型，包括：
将建筑群卫星影像与几何模型加载到地理信息系统平台中，依据建筑物编号，为各建筑物匹配几何模型，补充层高、层数、建筑面积和结构防火等级的模型参数，进而完成二维模型的建立。
3.根据权利要求1所述的一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述基于所述建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围，包括：
根据预设的建筑物间的火灾蔓延机制以及历史火灾统计数据，设置二维火灾模拟所需参数；其中，所述预设的火灾蔓延机制包括热辐射和热羽流，所述历史火灾统计数据包括起火点、起燃数量和气象条件；
根据预设的火灾蔓延经验公式、所述二维火灾模拟所需参数、以及所述建筑群二维模型，对建筑群的二维火灾蔓延范围进行预估。
4.根据权利要求1所述的一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述基于所述目标建筑群三维模型，进行建筑群高效三维火灾模拟，包括：
S41、根据目标建筑群三维模型建立相应的流体动力学模型；
S42、在流体动力学模型中，通过递归坐标二分法将二维火灾蔓延范围初步划分出各子区域；再根据火源位置的影响范围，在各子区域中划分不同精细度的网格；
S43、计算各子区域的负载均衡指标，根据各子区域的负载均衡指标、基于各子区域的网格负载调整网格边界，计算网格调整后的各子区域的工作负载，直到各子区域的工作负载近似相等；
S44、开始模拟，设定时间步，读取各子区域的计算时间，若各子区域的计算时间分布不均衡，利用负载迭代函数动态调整各子区域的工作负载，直至模拟结束。
5.根据权利要求4所述的一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述负载均衡指标为基于各子区域的计算时间测量值提出的量化指标，其中，所述负载均衡指标包括全域平均计算时间、子区域计算效率和累积不均衡负载。
6.根据权利要求4所述的一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述网格负载为基于网格边界长度和网格数目通过负载描述函数计算得到，其中，所述网格负载包括网格拓扑结构负载和通信负载。
7.根据权利要求4所述的一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述工作负载为计算结构平衡、通信平衡参数的量化指标，其中，所述工作负载包括结构负载平衡因子、通信负载平衡因子、火源辐射影响因子和可燃物含量因子。
8.根据权利要求4所述的一种基于二维预估和多尺度网格的火灾模拟方法，其特征在于，所述利用负载迭代函数动态调整各子区域的工作负载，包括：
根据负载均衡指标中的累积不均衡负载，计算得到转移负载和网格移动方向；
根据所述转移负载和所述网格移动方向，通过贪婪算法进行负载迭代，调整各子区域的网格边界，直至各子区域的工作负载达到近似相等。
9.一种基于二维预估和多尺度网格的火灾模拟装置，其特征在于，所述装置包括：
二维模型构建模块，用于导入建筑群数据，根据所述建筑群数据，获得建筑群二维模型；
三维模型构建模块，用于根据所述建筑群二维模型进行三维建模，获得建筑群三维模型；
范围预估模块，用于基于所述建筑群二维模型，根据火灾蔓延机制与历史火灾统计数据，预估二维火灾蔓延范围；
目标获取模块，用于根据所述二维火灾蔓延范围，基于所述建筑群三维模型，获得目标建筑群三维模型；
火灾模拟模块，用于基于所述目标建筑群三维模型，进行建筑群高效三维火灾模拟。
10.根据权利要求9所述的一种基于二维预估和多尺度网格的火灾模拟装置，其特征在于，所述火灾模拟模块，进一步用于：
S41、根据目标建筑群三维模型建立相应的流体动力学模型；
S42、在流体动力学模型中，通过递归坐标二分法将二维火灾蔓延范围初步划分出各子区域；再根据火源位置的影响范围，在各子区域中划分不同精细度的网格；
S43、计算各子区域的负载均衡指标，根据各子区域的负载均衡指标、基于各子区域的网格负载调整网格边界，计算网格调整后的各子区域的工作负载，直到各子区域的工作负载近似相等；
S44、开始模拟，设定时间步，读取各子区域的计算时间，若各子区域的计算时间分布不均衡，利用负载迭代函数动态调整各子区域的工作负载，直至模拟结束。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211513592.8A CN115601520B (zh) | 2022-11-30 | 2022-11-30 | 一种基于二维预估和多尺度网格的火灾模拟方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211513592.8A CN115601520B (zh) | 2022-11-30 | 2022-11-30 | 一种基于二维预估和多尺度网格的火灾模拟方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115601520A true CN115601520A (zh) | 2023-01-13 |
| CN115601520B CN115601520B (zh) | 2023-06-09 |
Family
ID=84853292
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211513592.8A Active CN115601520B (zh) | 2022-11-30 | 2022-11-30 | 一种基于二维预估和多尺度网格的火灾模拟方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115601520B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116580327A (zh) * | 2023-07-12 | 2023-08-11 | 四川安信科创科技有限公司 | 基于无人机的火势发展预测方法及预测系统 |
| CN116957484A (zh) * | 2023-06-08 | 2023-10-27 | 中国消防救援学院 | 一种基于数字孪生的消防救援预案方法及系统 |
| CN118190939A (zh) * | 2024-05-17 | 2024-06-14 | 中国林业科学研究院木材工业研究所 | 一种既有建筑木构件耐火性预测方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009193551A (ja) * | 2008-02-15 | 2009-08-27 | Advance Soft Kk | 建造物内における火災及び有毒物質拡散のシミュレーション技術 |
| CN106815075A (zh) * | 2016-12-15 | 2017-06-09 | 上海交通大学 | 建筑火灾数值仿真的区域分解优化方法 |
| CN108564867A (zh) * | 2018-03-08 | 2018-09-21 | 清华大学 | 城市地震次生火灾模拟方法及系统 |
| CN108733876A (zh) * | 2018-04-04 | 2018-11-02 | 北京交通大学 | 火灾演变仿真方法及包括其的火灾疏散综合仿真方法 |
-
2022
- 2022-11-30 CN CN202211513592.8A patent/CN115601520B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009193551A (ja) * | 2008-02-15 | 2009-08-27 | Advance Soft Kk | 建造物内における火災及び有毒物質拡散のシミュレーション技術 |
| CN106815075A (zh) * | 2016-12-15 | 2017-06-09 | 上海交通大学 | 建筑火灾数值仿真的区域分解优化方法 |
| CN108564867A (zh) * | 2018-03-08 | 2018-09-21 | 清华大学 | 城市地震次生火灾模拟方法及系统 |
| CN108733876A (zh) * | 2018-04-04 | 2018-11-02 | 北京交通大学 | 火灾演变仿真方法及包括其的火灾疏散综合仿真方法 |
Non-Patent Citations (2)
| Title |
|---|
| 江列霖 等: "基于负载均衡分区法的建筑火灾并行数值模拟及应用", 《上海交通大学学报》 * |
| 许镇 等: "考虑地面高程的建筑群三维火灾蔓延模型", 《清华大学学报(自然科学版)》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116957484A (zh) * | 2023-06-08 | 2023-10-27 | 中国消防救援学院 | 一种基于数字孪生的消防救援预案方法及系统 |
| CN116580327A (zh) * | 2023-07-12 | 2023-08-11 | 四川安信科创科技有限公司 | 基于无人机的火势发展预测方法及预测系统 |
| CN116580327B (zh) * | 2023-07-12 | 2023-09-29 | 四川安信科创科技有限公司 | 基于无人机的火势发展预测方法及预测系统 |
| CN118190939A (zh) * | 2024-05-17 | 2024-06-14 | 中国林业科学研究院木材工业研究所 | 一种既有建筑木构件耐火性预测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115601520B (zh) | 2023-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115601520A (zh) | 一种基于二维预估和多尺度网格的火灾模拟方法及装置 | |
| Liu et al. | A planning-oriented resilience assessment framework for transmission systems under typhoon disasters | |
| US20220309201A1 (en) | Artificial-intelligence-assisted method for providing urban design form and layout with improved wind environment | |
| CN106649960B (zh) | 大气多污染物环境容量三维迭代计算方法 | |
| JP5922570B2 (ja) | 拡散物質の拡散状況予測装置、方法およびプログラム | |
| Ghisu et al. | An improved cellular automata for wildfire spread | |
| CN103544385A (zh) | 顾及空间异质性的多尺度空间负荷预测方法 | |
| HomChaudhuri et al. | Genetic algorithm based simulation–optimization for fighting wildfires | |
| CN104182343B (zh) | 操作系统的性能测试方法及装置 | |
| Bozorgmehr et al. | Utilizing dynamic parallelism in CUDA to accelerate a 3D red-black successive over relaxation wind-field solver | |
| Nale et al. | Fragility functions for local failure mechanisms in unreinforced masonry buildings: a typological study in Ferrara, Italy | |
| CN117371233A (zh) | 一种基于ue技术的水库防洪系统 | |
| CN116415754A (zh) | 一种森林火灾火势蔓延预测方法及系统 | |
| CN114818418A (zh) | 基于主动学习多元自适应回归样条的边坡可靠度分析方法 | |
| Baum et al. | Simulation of large industrial outdoor fires | |
| CN117113788A (zh) | 一种火灾和爆炸耦合下钢制常压储罐结构失效预测方法 | |
| Liang et al. | Fractal Design of Indoor and Outdoor Forms of Architectural Space Based on a Three‐Dimensional Box Dimension Algorithm | |
| Wong et al. | Integrating computational fluid dynamics (CFD) models with GIS: An evaluation on data conversion formats | |
| Chen et al. | Wargame Mapping and Implementation for Emergency Evacuation of Residents in Urban Waterlogging Disaster | |
| Parkhomenko et al. | Improved computing performance and load balancing of atmospheric general circulation model | |
| RU2611892C1 (ru) | Способ трехмерного моделирования заданного гидрогеологического объекта, реализуемый в вычислительной системе | |
| Wang et al. | A procedure to simulate spread of post-earthquake fire in urban area considering seismic damage to buildings | |
| Zhou et al. | Dispersion prediction of pollutants in settlements based on generative adversarial networks | |
| CN116957414B (zh) | 基于人工智能的村庄规划分析方法及装置 | |
| CN118094697A (zh) | 一种考虑建筑倒塌的震后火灾蔓延模拟方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |