CN112381358B - 面向城市区域的绿化树木风灾破坏风险近实时评估方法 - Google Patents
面向城市区域的绿化树木风灾破坏风险近实时评估方法 Download PDFInfo
- Publication number
- CN112381358B CN112381358B CN202011165545.XA CN202011165545A CN112381358B CN 112381358 B CN112381358 B CN 112381358B CN 202011165545 A CN202011165545 A CN 202011165545A CN 112381358 B CN112381358 B CN 112381358B
- Authority
- CN
- China
- Prior art keywords
- wind speed
- situation
- target
- tree
- computational fluid
- 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.)
- Active
Links
- 230000006378 damage Effects 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 59
- 239000012530 fluid Substances 0.000 claims abstract description 125
- 230000002968 anti-fracture Effects 0.000 claims abstract description 26
- 238000004088 simulation Methods 0.000 claims description 46
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 6
- 238000012502 risk assessment Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Primary Health Care (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请提出了一种面向城市区域的绿化树木风灾破坏风险近实时评估方法,包括:获取目标区域的建筑数据和树木数据,根据建筑数据和树木数据建立目标区域的计算流体力学几何模型,根据流体力学几何模型建立目标区域的每个风向的计算流体力学模型;根据计算流体力学模型生成目标区域的情境库;采用基于情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速;计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,将树顶风速与抗折断极限风速、抗倾覆极限风速进行比较,生成目标区域中每棵树木的风灾破坏状态。根据本申请能够实现目标情境下目标区域内每棵绿化树木风灾破坏风险的近实时评估。
Description
技术领域
本申请涉及数据处理技术领域,尤其涉及一种面向城市区域的绿化树木风灾破坏风险近实时评估方法。
背景技术
风荷载引起的城市绿化树木破坏会造成直接或间接的经济损失,例如:城市在受到台风侵袭后,城市内树木受损,道路上倒伏和折断的树木对城市交通造成严重影响。因此,有必要在强风天气到来之前提前对城市区域的绿化树木进行针对性加固,从而减少风灾带来的损失。
合理、科学地评估风灾下城市区域绿化树木的破坏风险,可以帮助相关人员快速识别特定风灾情境下破坏风险较高的树木,从而提高城市树木防风工作开展的效率。然而,现有的针对树林系统的树木风灾破坏力学模型,适用于单个林分内树木风灾破坏的预测分析,而一个真实城市区域包含复杂的建筑群和多个林分对象,现有力学模型难以量化城市区域内不同林分之间以及林分与其周围建成环境之间的流体力学相互影响,导致难以准确评估城市内绿化树木的风灾破坏风险。此外,强风天气的多变性也对树木风灾破坏风险评估方法的效率提出了要求,气象预测信息随时间不断更新,风险评估结果需要随之快速更新,才能为城市防风工作提供有效的技术支持。
发明内容
本申请旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本申请提出一种面向城市区域的绿化树木风灾破坏风险近实时评估方法,以实现目标风灾情境下城市区域内绿化树木风灾破坏风险的评估。
本申请第一方面实施例提出了一种面向城市区域的绿化树木风灾破坏风险近实时评估方法,包括:
获取目标区域的建筑数据和树木数据,根据所述建筑数据和树木数据建立所述目标区域的计算流体力学几何模型,根据所述流体力学几何模型建立所述目标区域的每个风向的计算流体力学模型;
获取所述目标区域的情境库;其中,所述情景库通过如下方式得到:针对所述目标区域的每个风向的计算流体力学模型,改变入口边界的输入风速剖面、湍动能剖面、湍动能耗散率剖面,生成计算流体力学模型库;基于雷诺平均方程对所述模型库中的计算流体力学模型进行计算流体力学模拟,得到每个计算流体力学模型的平均风速模拟结果;根据所述每个计算流体力学模型的平均风速模拟结果,生成所述目标区域的情境库;
采用基于所述情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速;
计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,将所述树顶风速与所述抗折断极限风速、所述抗倾覆极限风速进行比较,生成所述目标区域中每棵树木的风灾破坏状态。
根据本申请实施例的面向城市区域的绿化树木风灾破坏风险近实时评估方法,可以实现目标情境下目标区域内每棵绿化树木风灾破坏风险的近实时评估,提高城市内绿化树木的风灾破坏风险评估准确性,为城市防风工作提供技术支持。
本申请第二方面实施例提出了一种面向城市区域的绿化树木风灾破坏风险近实时评估装置,包括:
建立模块,用于获取目标区域的建筑数据和树木数据,根据所述建筑数据和树木数据建立所述目标区域的计算流体力学几何模型,根据所述流体力学几何模型建立所述目标区域的每个风向的计算流体力学模型;
获取模块,用于获取所述目标区域的情境库;其中,所述情景库通过如下方式得到:针对所述目标区域的每个风向的计算流体力学模型,改变入口边界的输入风速剖面、湍动能剖面、湍动能耗散率剖面,生成计算流体力学模型库;基于雷诺平均方程对所述模型库中的计算流体力学模型进行计算流体力学模拟,得到每个计算流体力学模型的平均风速模拟结果;根据所述每个计算流体力学模型的平均风速模拟结果,生成所述目标区域的情境库;
生成模块,用于采用基于所述情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速;
处理模块,用于计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,将所述树顶风速与所述抗折断极限风速、所述抗倾覆极限风速进行比较,生成所述目标区域中每棵树木的风灾破坏状态。
本申请第三方面实施例提出了一种计算机设备,包括处理器和存储器;其中,所述处理器通过读取所述存储器中存储的可执行程序代码来运行与所述可执行程序代码对应的程序,以用于实现如第一方面实施例所述的面向城市区域的绿化树木风灾破坏风险近实时评估方法。
本申请第四方面实施例提出了一种非临时性计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如第一方面实施例所述的面向城市区域的绿化树木风灾破坏风险近实时评估方法。
本申请附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本申请的实践了解到。
附图说明
图1为本申请实施例所提供的一种面向城市区域的绿化树木风灾破坏风险近实时评估方法的流程示意图;
图2为本申请实施例所提供的一种面向城市区域的绿化树木风灾破坏风险近实时评估方法的框架图;
图3为本申请实施例所提供的一种真实复杂城市区域的计算流体力学模型示意图;
图4为本申请实施例所提供的一种真实复杂城市区域的计算流体力学模型库示意图;
图5为本申请实施例所提供的一种线性插值方法流程图;
图6为本申请实施例所提供的一种真实复杂城市区域在某具体风灾情境下的树木破坏状态示意图;
图7为本申请实施例所提供的一种面向城市区域的绿化树木风灾破坏风险近实时评估装置的结构示意图。
具体实施方式
下面详细描述本申请的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。
下面参考附图描述本申请实施例的面向城市区域的绿化树木风灾破坏风险近实时评估方法、装置及设备。
图1为本申请实施例所提供的一种面向城市区域的绿化树木风灾破坏风险近实时评估方法的流程示意图,如图1所示,该方法包括:
步骤101,获取目标区域的建筑数据和树木数据,根据建筑数据和树木数据建立目标区域的计算流体力学几何模型,根据流体力学几何模型建立目标区域的每个风向的计算流体力学模型。
具体地,建筑数据包括建筑底面多边形每个角点的坐标数据和建筑高度数据,树木数据包括胸径(Dbh,1.3m高度处的树干直径)、树木高度(H)、树冠深度(Dcrown)、树冠宽度(Wcrown)、树种信息以及所在土壤的类型。
本实施例中,将H数据相近、Dcrown数据相近、Wcrown数据相近、且品种相同、位置临近的树木进行组合形成林分对象,根据建筑数据和树木数据建立目标区域包含建筑和林分对象的计算流体力学几何模型。几何模型中可以忽略树冠以下的树干部分对风场的影响,通过距离地面预设高度的平行六面体表征每个林分对象。平行六面体的距地高度为一片林分所有树木树冠底部高度的平均值,顶部高度为树木高度的平均值,宽度为树冠宽度的平均值。
其中,计算流体力学几何模型包括表征建筑外形的建筑模型和表征树木影响区域的树木模型,对于计算流体力学几何模型中的树木模型,修正树木模型内部空间的Navier-Stokes方程和k-ε湍流模型。进而,根据修正后的计算流体力学几何模型,建立目标区域的每个风向的计算流体力学模型。
具体地,在本实施例中,考虑8种不同的风向,包括北风(N)、东北风(NE)、东风(E)、东南风(SE)、南风(S)、西南风(SW)、西风(W)和西北风(NW)。对于计算流体力学模型中非树木所在的流域空间,采用未进行修正的Navier-Stokes方程和k-ε湍流模型。
在本申请的一个实施例中,采用计算流体力学代码ANSYS FLUENT建立计算流体力学模型。计算流体力学模型包括:计算域、网格划分、计算流体力学参数设置和计算流体力学几何模型。
步骤102,获取目标区域的情境库。
其中,情景库通过如下方式得到:S1、针对目标区域的每个风向的计算流体力学模型,改变入口边界的输入风速剖面、湍动能剖面、湍动能耗散率剖面,生成计算流体力学模型库;S2、基于雷诺平均方程对所述模型库中的计算流体力学模型进行计算流体力学模拟,得到每个计算流体力学模型的平均风速模拟结果;S3、根据所述每个计算流体力学模型的平均风速模拟结果,生成目标区域的情境库。
对于上述步骤S1,具体地,风速剖面根据冯卡门常数、摩擦速度、地面粗糙度高度(根据目标区域所处地区的地面粗糙度而定)、距地面高度得到,湍动能剖面根据摩擦速度、参数cμ(取0.09)得到,湍动能耗散率剖面根据摩擦速度、冯卡门常数、地面粗糙度高度、距地面高度得到。
针对目标区域的每个风向的计算流体力学模型,入口边界采用7组输入风速剖面、湍动能剖面、湍动能耗散率剖面,三种剖面的形式由10m高处风速唯一确定,7组剖面对应的10m高处风速分别为0m/s、5m/s、10m/s、15m/s、20m/s、25m/s、30m/s。
具体地,在采用8种不同的风向和7种不同的入口边界条件后,共计得到8×7=56个计算流体力学模型,每个模型对应的情境采用Case-S-D的方式命名,其中,S表示风速大小,即0m/s、5m/s、10m/s、15m/s、20m/s、25m/s、30m/s;D表示风向,即N、NE、E、SE、S、SW、W、NW。
对于上述步骤S2,在形成计算流体力学模型库后,对计算流体力学模型库中的计算流体力学模型进行基于雷诺平均方程的计算流体力学模拟,得到每个计算流体力学模型的平均风速模拟结果。
对于上述步骤S3,保存每个计算流体力学模型的平均风速模拟结果,形成目标区域的计算流体力学模拟情境库。
步骤103,采用基于计算流体力学模拟情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速。
本实施例中,根据气象预测信息确定目标情境的入口边界风速和风向,在情境库中确定上部相邻情境和下部相邻情境。
在计算流体力学模拟情境库的所有情境中,上部相邻情境为与目标情境风向相同、入口边界风速不小于目标情境入口边界风速的至少一个情境中入口边界风速最小的情境。
在计算流体力学模拟情境库的所有情境中,下部相邻情境为与目标情境风向相同、入口边界风速不大于目标情境入口边界风速的至少一个情境中入口边界风速最大的情境。
作为一种示例,如以风向为南风、入口边界10m高处风速为18m/s的目标情境Case-18-S为例,则其上部相邻情境为Case-20-S,下部相邻情境为Case-15-S。
在获取目标情境的上部相邻情境和下部相邻情境后,对于目标情境下目标区域中的目标树木,根据上部相邻情境中目标树木的树顶风速模拟结果与下部相邻情境中目标树木的树顶风速模拟结果的插值,生成目标树木的树顶风速。
具体地,采用线性内插的方式,以目标情境Case-18-S下树木i的树顶风速计算为例,树木i在Case-18-S下的树顶风速Ui,Case-18-S采用其在Case-20-S下的树顶风速Ui,Case-20-S和其在Case-15-S下的树顶风速Ui,Case-15-S线性内插得到,为(3Ui,Case-20-S/5+2Ui,Case-15-S/5)。
根据基于所述计算流体力学模拟情境库的插值方法,生成目标情境下目标区域的所有树木的树顶风速,从而实现目标情境下目标区域所有树木的树顶风速计算。
步骤104,计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,将所述树顶风速与所述抗折断极限风速、所述抗倾覆极限风速进行比较,生成所述目标区域中每棵树木的风灾破坏状态。
具体地,采用GALES树木风灾破坏力学模型计算目标区域内每棵树木的抗折断极限风速和抗倾覆极限风速。
本实施例中,对于目标区域内的每棵树木,比较抗折断极限风速和抗倾覆极限风速的大小;将抗折断极限风速小于抗倾覆极限风速的树木标记为A类树木,将抗折断极限风速大于抗倾覆极限风速的树木标记为B类树木,将抗折断极限风速等于抗倾覆极限风速的树木标记为C类树木。
对于A类树木,将所述树顶风速与所述抗折断极限风速进行比较,若树顶风速小于抗折断极限风速,则预测树木在目标情境下不会破坏,否则,预测树木在目标情境下发生树干折断破坏;
对于B类树木,将所述树顶风速与抗倾覆极限风速进行比较,若树顶风速小于抗倾覆极限风速,则预测树木在目标情境下不会破坏,否则,预测树木在目标情境下可能发生整体倾覆破坏;
对于C类树木,将树顶风速与抗折断极限风速进行比较,若树顶风速小于抗折断极限风速,则预测树木在目标情境下不会破坏,否则,预测树木在目标情境下可能发生树干折断破坏或者整体倾覆破坏。
其中,抗折断极限风速为使树木发生树干折断破坏的最小树顶风速,抗倾覆极限风速为使树木发生整体倾覆破坏的最小树顶风速。
根据本申请实施例的方法,实现了目标情境下目标区域内每棵树木的风灾破坏风险的近实时评估,提高城市内绿化树木的风灾破坏风险评估准确性。
图2为本申请实施例所提供的一种面向城市区域的绿化树木风灾破坏风险近实时评估方法的框架图,如图2所示,本申请实施例的方法主要分为四个模块。
模块一:建立计算流体力学模型。
基于目标区域的建筑数据和树木数据,建立包含建筑和树木的计算流体力学几何模型,在ANSYS FLUENT中对树木占据的流域空间的流体控制方程进行修正,以模拟树木对城市风场的影响,由此建立不同风向情境的计算流体力学模型。
模块二:构建计算流体力学模拟情境库。
针对不同风向情境的计算流体力学模型,改变每个模型的入口边界条件,形成计算流体力学模型库。对每个计算流体力学模型进行基于雷诺平均方程的计算流体力学模拟,保存每个模型的平均风速模拟结果,构建目标区域的计算流体力学模拟情境库。
模块三:近实时预测目标情境下每棵树木的树顶风速。
依据气象预测信息得到目标情境的10m高处的风速数据和风向数据,从计算流体力学模拟情境库中获取目标情境的上部相邻情境和下部相邻情境。采用线性插值方法,近实时预测目标情境下目标区域内每棵树木的树顶风速。
模块四:评估目标情境下每棵树木的风灾破坏状态。
采用树木风灾破坏力学模型计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,并将抗折断极限风速、抗倾覆极限风速与模块三中预测得到的树顶风速进行比较,得到每棵树木具体的风灾破坏状态。
下面针对计算流体力学模型进行说明。
图3为本申请实施例所提供的一种真实复杂城市区域的计算流体力学模型示意图,如图3所示,展示了一个真实复杂城市区域的计算流体力学模型中的建筑模型和树木模型,以及对树木所占据的流域空间内流体的控制方程的修正方法。
具体地,在本申请的一个具体的实施例中,树木占据的流域空间的流体控制方程的修正方法分为两个步骤:
步骤1:对于树木占据的流域空间的雷诺平均方程中的动量守恒方程,添加源项Fd,Fd采用式(1)所示形式:
Fd=-ρ·cd·LAD·|U|·ui(1)
其中,ρ是空气密度;cd是机械阻力系数,根据树木冠层具体特征确定;LAD表示树木的叶面积密度,定义为单位群落体积的总植物单面叶面积,根据树木冠层具体特征确定;ui是风速;|U|是风速向量的模,可以通过式(2)进行计算:
|U|=(ui·ui)1/2(2)
步骤2:对于湍动能(k)和湍动能耗散率(ε)的输运方程,分别添加源项Sk和Sε,Sk和Sε分别采用式(3)和(4)所示形式。
Sk=ρ·cd·LAD·(βp|U|3-βd|U|k) (3)
其中,cε4、cε5、βp、βd为修正源项常数,根据树木冠层具体特征确定。
下面对计算流体力学模型库进行说明。
图4为本申请实施例所提供的一种真实复杂城市区域的计算流体力学模型库示意图。具体地,如图4所示,展示了本申请一个实施例的真实复杂城市区域的8种不同风向的计算流体力学模型的计算域尺寸、网格单元类型、网格单元数量。计算流体力学模型库中共计56个模型,相同风向、不同入口边界条件的模型在计算域尺寸、网格单元类型、网格单元数量上相同。对56个模型进行基于雷诺平均方程的计算流体力学模拟,即可得到每个模型的平均风速模拟结果,从而实现目标区域计算流体力学模拟情境库的构建。
下面对插值方法进行说明。
图5为本申请实施例所提供的一种线性插值方法流程图,如图5所示,线性插值方法的流程分为以下4个步骤。
步骤1:获取目标情境Case-S-D的入口边界风速信息S和风向信息D。
步骤2:判断目标情境的入口边界风速处于计算流体力学模拟情境库中哪两个情境的入口边界风速之间,从情境库中筛选出上部相邻情境Case-k-D和下部相邻情境Case-(k-5)-D。
步骤3:采用线性插值方法计算目标情境下第i棵树木的树顶风速Ui,目标情境,Ui,目标情境=(S-k+5)(Ui,Case-k-D-Ui,Case-(k-5)-D)/5+Ui,Case-(k-5)-D。
步骤4:对于目标区域内的所有树木,重复步骤3,得到目标区域内所有树木的树顶风速。
图6为本申请实施例所提供的一种真实复杂城市区域在某具体风灾情境下的树木破坏状态示意图。如图6所示,展示了根据本申请一个实施例的真实复杂城市区域在一西北风、10m高处风速为15.5m/s的风灾情境下的每棵绿化树木的破坏状态。
具体地,在本申请的一个具体的实施例中,树木的抗折断极限风速和抗倾覆极限风速采用GALES树木力学模型进行计算。图6中黑色块体代表建筑,灰色圆形代表不发生破坏的树木,标记为1的正三角形代表发生树干折断破坏的树木,标记为2的倒三角形代表发生整体倾覆破坏的树木。在真实复杂城市区域中,西部区域的一些树木由于长细比较大,在风灾情境下发生了树干折断破坏;中部区域的少数几棵树木将发生整体倾覆破坏。可以看到,通过本申请实施例提出的方法实现了面向城市区域尺度的绿化树木风灾破坏风险评估。由于在本申请实施例提出的方法中,计算耗时较长的计算流力学模拟发生在对目标情境分析之前,在对目标情境分析时只涉及初级的数学运算,因此本申请实施例的评估流程耗时不超过半小时。表明本申请实施例提出的方法能够实现面向城市区域尺度的绿化树木风灾破坏风险近实时评估,可以帮助城市管理者快速识别目标情境下目标区域内风灾破坏风险较高的树木,为城市树木防风工作提供技术支持。
综上,根据本申请实施例提出的面向城市区域的绿化树木风灾破坏风险近实时评估方法,通过建立目标区域包含建筑和树木的计算流体力学模型与修正树木所在空间内流体的Navier-Stokes方程和k-ε湍流模型,实现不同风向情境的计算流体力学模型的建立,该模型可以量化树木对城市风场的影响;通过改变不同风向情境的计算流体力学模型的入口边界条件,建立目标区域的计算流体力学模型库;通过对目标区域计算流体力学模型库的模型进行计算流体力学模拟,得到不同风向、不同入口边界条件的若干个情境的平均风速模拟结果,形成目标区域的计算流体力学模拟情境库;根据目标情境的入口边界风速信息和风向信息,提取目标情境在情境库中的上部相邻情境和下部相邻情境;通过对上部相邻情境和下部相邻情境树顶风速模拟结果进行插值,近实时预测目标情境下目标区域内每棵树木的树顶风速结果;通过树木风灾破坏力学模型计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,并通过比较树顶风速预测结果与抗折断极限风速和抗倾覆极限风速的大小来判断树木破坏状态。由此,可以实现目标情境下目标区域内每棵绿化树木风灾破坏风险的近实时评估,为城市防风工作提供技术支持。
为了实现上述实施例,本申请还提出一种面向城市区域的绿化树木风灾破坏风险近实时评估装置。
图7为本申请实施例所提供的一种面向城市区域的绿化树木风灾破坏风险近实时评估装置的结构示意图,包括:建立模块10,获取模块20,生成模块30,处理模块40。
其中,建立模块10,用于获取目标区域的建筑数据和树木数据,根据建筑数据和树木数据建立目标区域的计算流体力学几何模型,根据流体力学几何模型建立所述目标区域的每个风向的计算流体力学模型。
获取模块20,用于获取所述目标区域的情境库;其中,所述情景库通过如下方式得到:针对所述目标区域的每个风向的计算流体力学模型,改变入口边界的输入风速剖面、湍动能剖面、湍动能耗散率剖面,生成计算流体力学模型库;基于雷诺平均方程对所述模型库中的计算流体力学模型进行计算流体力学模拟,得到每个计算流体力学模型的平均风速模拟结果;根据所述每个计算流体力学模型的平均风速模拟结果,生成所述目标区域的情境库。
生成模块30,用于采用基于所述情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速。
处理模块40,用于计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,将所述树顶风速与所述抗折断极限风速、所述抗倾覆极限风速进行比较,生成目标区域中每棵树木的风灾破坏状态。
前述实施例对面向城市区域的绿化树木风灾破坏风险近实时评估方法的解释说明同样适用于本实施例的面向城市区域的绿化树木风灾破坏风险近实时评估装置,此处不再赘述。
根据本申请实施例的面向城市区域的绿化树木风灾破坏风险近实时评估装置,能够实现目标情境下目标区域内每棵绿化树木风灾破坏风险的近实时评估,为城市防风工作提供技术支持。
为了实现上述实施例,本申请还提出一种计算机设备,包括处理器和存储器;其中,处理器通过读取存储器中存储的可执行程序代码来运行与可执行程序代码对应的程序,以用于实现如前述任一实施例所述的面向城市区域的绿化树木风灾破坏风险近实时评估方法。
为了实现上述实施例,本申请还提出一种计算机程序产品,当计算机程序产品中的指令被处理器执行时实现如前述任一实施例所述的面向城市区域的绿化树木风灾破坏风险近实时评估方法。
为了实现上述实施例,本申请还提出一种非临时性计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如前述任一实施例所述的面向城市区域的绿化树木风灾破坏风险近实时评估方法。
在本申请的描述中,需要理解的是,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本申请的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本申请的限制,本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (9)
1.一种面向城市区域的绿化树木风灾破坏风险近实时评估方法,其特征在于,包括:
获取目标区域的建筑数据和树木数据,根据所述建筑数据和树木数据建立所述目标区域的计算流体力学几何模型,根据所述流体力学几何模型建立所述目标区域的每个风向的计算流体力学模型;
获取所述目标区域的情境库;其中,所述情境库通过如下方式得到:针对所述目标区域的每个风向的计算流体力学模型,改变入口边界的输入风速剖面、湍动能剖面、湍动能耗散率剖面,生成计算流体力学模型库;基于雷诺平均方程对所述模型库中的计算流体力学模型进行计算流体力学模拟,得到每个计算流体力学模型的平均风速模拟结果;根据所述每个计算流体力学模型的平均风速模拟结果,生成所述目标区域的情境库;
采用基于所述情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速;
计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,将所述树顶风速与所述抗折断极限风速、所述抗倾覆极限风速进行比较,生成所述目标区域中每棵树木的风灾破坏状态;
其中,所述采用基于所述情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速,包括:
根据气象预测信息或感兴趣的假想情境确定所述目标情境的入口边界风速和风向,在所述情境库中确定上部相邻情境和下部相邻情境;
对于所述目标情境下所述目标区域中的目标树木,根据上部相邻情境中所述目标树木的树顶风速模拟结果与下部相邻情境中所述目标树木的树顶风速模拟结果的插值,生成所述目标树木的树顶风速;
根据所述基于所述情境库的插值方法,生成所述目标情境下所述目标区域的所有树木的树顶风速。
2.如权利要求1所述的方法,其特征在于,所述计算流体力学几何模型包括表征建筑外形的建筑模型和表征树木影响区域的树木模型,所述根据所述流体力学几何模型建立所述目标区域的每个风向的计算流体力学模型包括:
对于所述计算流体力学几何模型中的树木模型,修正所述树木模型内部空间的Navier-Stokes方程和k-ε湍流模型;
根据修正后的计算流体力学几何模型,建立所述目标区域的每个风向的计算流体力学模型。
3.如权利要求1所述的方法,其特征在于,所述计算流体力学模型包括:计算域、网格划分、计算流体力学参数设置和计算流体力学几何模型。
4.如权利要求1所述的方法,其特征在于,所述上部相邻情境包括:在所述情境库的所有情境中,与所述目标情境风向相同、入口边界风速不小于目标情境入口边界风速的至少一个情境中入口边界风速最小的情境。
5.如权利要求1所述的方法,其特征在于,所述下部相邻情境包括:在所述情境库的所有情境中,与所述目标情境风向相同、入口边界风速不大于目标情境入口边界风速的至少一个情境中入口边界风速最大的情境。
6.如权利要求1所述的方法,其特征在于,所述将所述树顶风速与所述抗折断极限风速、所述抗倾覆极限风速进行比较,生成所述目标区域中每棵树木的风灾破坏状态,包括:
对于所述目标区域内的每棵树木,比较抗折断极限风速和抗倾覆极限风速的大小;
将抗折断极限风速小于抗倾覆极限风速的树木标记为A类树木,将抗折断极限风速大于抗倾覆极限风速的树木标记为B类树木,将抗折断极限风速等于抗倾覆极限风速的树木标记为C类树木;
对于所述A类树木,将所述树顶风速与所述抗折断极限风速进行比较,若所述树顶风速小于所述抗折断极限风速,则生成树木在所述目标情境下不会破坏的预测结果,否则,生成树木在所述目标情境下发生树干折断破坏的预测结果;
对于所述B类树木,将所述树顶风速与所述抗倾覆极限风速进行比较,若树顶风速小于抗倾覆极限风速,则生成树木在所述目标情境下不会破坏的预测结果,否则,生成树木在所述目标情境下发生整体倾覆破坏的预测结果;
对于所述C类树木,将所述树顶风速与所述抗折断极限风速进行比较,若树顶风速小于抗折断极限风速,则生成树木在所述目标情境下不会破坏的预测结果,否则,生成树木在所述目标情境下发生树干折断破坏或者整体倾覆破坏的预测结果。
7.如权利要求1所述的方法,其特征在于,所述抗折断极限风速为使树木发生树干折断破坏的最小树顶风速。
8.根据权利要求1所述的方法,其特征在于,所述抗倾覆极限风速为使树木发生整体倾覆破坏的最小树顶风速。
9.一种面向城市区域的绿化树木风灾破坏风险近实时评估装置,其特征在于,包括:
建立模块,用于获取目标区域的建筑数据和树木数据,根据所述建筑数据和树木数据建立所述目标区域的计算流体力学几何模型,根据所述流体力学几何模型建立所述目标区域的每个风向的计算流体力学模型;
获取模块,用于获取所述目标区域的情境库;其中,所述情境库通过如下方式得到:针对所述目标区域的每个风向的计算流体力学模型,改变入口边界的输入风速剖面、湍动能剖面、湍动能耗散率剖面,生成计算流体力学模型库;基于雷诺平均方程对所述模型库中的计算流体力学模型进行计算流体力学模拟,得到每个计算流体力学模型的平均风速模拟结果;根据所述每个计算流体力学模型的平均风速模拟结果,生成所述目标区域的情境库;
生成模块,用于采用基于所述情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速;
处理模块,用于计算目标区域每棵树木的抗折断极限风速和抗倾覆极限风速,将所述树顶风速与所述抗折断极限风速、所述抗倾覆极限风速进行比较,生成所述目标区域中每棵树木的风灾破坏状态;
其中,所述采用基于所述情境库的插值方法得到目标情境下目标区域每棵树木的树顶风速,包括:
根据气象预测信息或感兴趣的假想情境确定所述目标情境的入口边界风速和风向,在所述情境库中确定上部相邻情境和下部相邻情境;
对于所述目标情境下所述目标区域中的目标树木,根据上部相邻情境中所述目标树木的树顶风速模拟结果与下部相邻情境中所述目标树木的树顶风速模拟结果的插值,生成所述目标树木的树顶风速;
根据所述基于所述情境库的插值方法,生成所述目标情境下所述目标区域的所有树木的树顶风速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011165545.XA CN112381358B (zh) | 2020-10-27 | 2020-10-27 | 面向城市区域的绿化树木风灾破坏风险近实时评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011165545.XA CN112381358B (zh) | 2020-10-27 | 2020-10-27 | 面向城市区域的绿化树木风灾破坏风险近实时评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112381358A CN112381358A (zh) | 2021-02-19 |
CN112381358B true CN112381358B (zh) | 2022-09-02 |
Family
ID=74577326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011165545.XA Active CN112381358B (zh) | 2020-10-27 | 2020-10-27 | 面向城市区域的绿化树木风灾破坏风险近实时评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112381358B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114137636B (zh) * | 2021-11-11 | 2022-08-12 | 四川九通智路科技有限公司 | 一种环形压力传感器的区域气象监测管理方法及系统 |
CN114396975B (zh) * | 2022-01-05 | 2023-05-23 | 深圳市金阳盛城市服务集团有限公司 | 一种树木抗倒伏监控预警方法、装置、电子设备及介质 |
CN115935855B (zh) * | 2023-01-09 | 2023-05-09 | 北京科技大学 | 基于优化树木花粉浓度指标的城市绿化方法及装置 |
CN116523663B (zh) * | 2023-06-29 | 2023-09-01 | 北京科技大学 | 基于物理模型的城市建筑群风灾保险保费厘定方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101996283B (zh) * | 2010-11-26 | 2012-08-22 | 上海市浦东新区气象局 | 行道树城市街区风灾的动态预报方法 |
CN105956753A (zh) * | 2016-04-26 | 2016-09-21 | 同济大学 | 一种城市行道树风致破坏的快速评估方法 |
US11215597B2 (en) * | 2017-04-11 | 2022-01-04 | Agerpoint, Inc. | Forestry management tool for assessing risk of catastrophic tree failure due to weather events |
CN111553593B (zh) * | 2020-04-27 | 2023-07-18 | 上海应用技术大学 | 一种针对城市行道树的安全风险评估方法及系统 |
-
2020
- 2020-10-27 CN CN202011165545.XA patent/CN112381358B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112381358A (zh) | 2021-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112381358B (zh) | 面向城市区域的绿化树木风灾破坏风险近实时评估方法 | |
Vivoni et al. | Generation of triangulated irregular networks based on hydrological similarity | |
Brown et al. | The influence of lifting condensation level on low-level outflow and rotation in simulated supercell thunderstorms | |
CN111612806B (zh) | 一种建筑物立面窗户提取方法及装置 | |
CN115329691B (zh) | 一种基于cfd与gis的超大城市风环境模拟方法 | |
KR101943758B1 (ko) | 도시내 위험물질 확산에 따른 보행자환경에서의 위험영향 평가 방법 및 장치 | |
Etienne et al. | Spatial predictions of extreme wind speeds over Switzerland using generalized additive models | |
Fischer et al. | A spatial and seasonal climatology of extreme precipitation return-levels: a case study | |
CN105513051B (zh) | 一种点云数据处理方法和设备 | |
CN115470731A (zh) | 一种基于微气象的风场风电功率预测方法及系统 | |
CN110765694B (zh) | 一种基于简化型浅水方程组的城市地表水流数值模拟方法 | |
KR101463493B1 (ko) | 추계학적 태풍모의를 적용한 전지구 기후모델의 보완방법 | |
CN110852472A (zh) | 一种基于随机森林算法的陆地水储量预测方法及设备 | |
CN117973131A (zh) | 基于有限元模型的输电塔地质灾害失效预测方法及系统 | |
Walsworth et al. | Image modelling of forest changes associated with acid mine drainage | |
Lukač et al. | Wind resource assessment using airborne LiDAR data and smoothed particle hydrodynamics | |
WO2019128018A1 (zh) | 确定风资源的方法和装置 | |
CN108121888B (zh) | 沙滩潮间带剖面平衡坡度的分析方法 | |
Erlacher et al. | A GPU-based parallelization approach to conduct spatially-explicit uncertainty and sensitivity analysis in the application domain of landscape assessment | |
CN104732292A (zh) | 基于断面数据的主流速线自动规划方法 | |
Susiloningtyas et al. | Dynamic System Model of Land Use Affected by Sea Level Rise in the Coastal Area of Bengkulu City, Indonesia | |
Zhang et al. | Conformal adaptive hexahedral-dominant mesh generation for CFD simulation in architectural design applications | |
IŞINKARALAR | Spatio-temporal change of the morphology in west corridor development region of Ankara city and 2022-2039 growth estimation. | |
Ward et al. | An iterative method for the creation of structured hexahedral meshes over complex orography | |
Zhang et al. | Fast simulation method for ocean wave base on ocean wave spectrum and improved gerstner model with GPU |
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 |