PHAST-SAFETI二次开发教程（07）：气象与地形工况——Weather、大气稳定度、地面条件与建筑阻塞区

PHAST-SAFETI二次开发教程07气象与地形工况——Weather、大气稳定度、地面条件与建筑阻塞区版本声明块工具/软件Phast™ / Safeti™版本主线 9.x具体以官方发布说明为准语言/环境Python ≥ 3.11官方库 PyPWS包名pypws桌面侧气象/地图设置以官方 User Manual 为准本文目标把「气象 × 地面 × 建筑/阻塞」这组环境输入整理成可枚举、可复现、可审计的工况矩阵一句话结论环境输入不是一个数而是一组维度——官方 PyPWS 把气象与地面抽象为Weather与Substrate两个实体、把大气稳定度抽象为AtmosphericStabilityClass、把地面类型抽象为SurfaceType与PoolSurfaceType批量场景的气象工况数等于各维度取值数的连乘例如风速 3 档 × 稳定度 6 档 × 地面 2 类即 36 个气象工况此处数值仅作组合逻辑示意因此「工况爆炸」必须在生成阶段就被算清并去重。本篇铁律适用三级标注取自系列登记的 10 条[硬约束] 铁律 4 单位与坐标系纪律风速、温度等一律 SI地图导入前先声明 GIS 坐标系基准。[硬约束] 铁律 2 闭源事实纪律实体属性名与枚举取值以官方 PyPWS 参考文档为准未确证字段在代码中显式标注待核对。[强约束] 铁律 6 批处理幂等 可续跑工况矩阵必须可重放命名唯一。[提示] 铁律 8 物理量不臆造稳定度分类、粗糙度等取值来自官方技术文档或项目法规依据不自造。〇、本篇要解决的认知问题一个Weather由哪些维度构成哪些是自由度、哪些是派生量大气稳定度在数据模型里是什么AtmosphericStabilityClass与风速是什么关系地面条件为什么要拆成Substrate、SurfaceType、PoolSurfaceType三层建筑尾流/室内释放在什么场景下必须打开对结果影响在哪一环地图上的阻塞区congested regions为什么是爆炸建模的输入而不是扩散的输入把气象做成批量维度时怎么避免组合爆炸与重复工况以下「认知问题回显FAQ」与上列问题一一对应、同序作答。一、机制解析1.1 环境输入的三个维度官方 PyPWS 的输入实体清单里环境输入占两个Weather与Substrate枚举清单里与它们相关的至少四个AtmosphericStabilityClass、SurfaceType、PoolSurfaceType以及与热辐射类结果相关的RadiationType。官方对泄漏计算的描述也点明「关联数据对象如 leaks、pipe ruptures、weathers 在调用计算时被消费」说明气象是被消费的独立数据对象而不是计算对象的一个散装字段。维度数据载体典型自由度作用环节风Weather风速、风向量扩散UDM、火焰倾角大气条件WeatherAtmosphericStabilityClass稳定度类别扩散的湍流参数地面/基底Substrate地表类型、粗糙度类特征扩散的近地面边界地表介质SurfaceType/PoolSurfaceType陆地、水面等液池所处表面液池扩展与蒸发建构筑物场景级设置尾流模型、室内释放开关扩散初始条件阻塞区地图/plot plan 几何阻塞区空间范围爆炸超压建模1.2 稳定度的位置风速 / 风向 │ ▼ ┌─────────────────┐ 稳定度类别AtmosphericStabilityClass │ Weather │◀──── 由风速与大气条件共同确定落在环境对象上 └────────┬────────┘ │ ▼ ┌─────────────────┐ 地面条件Substrate / SurfaceType / PoolSurfaceType │ 扩散UDM │◀──── 近地面边界与地表介质 └────────┬────────┘ │ ┌────────┴────────┐ ▼ ▼ 火灾/毒性结果 阻塞区几何 ──▶ 爆炸超压建模这张图要读出两条「不连线」的关系阻塞区不参与扩散的湍流参数它服务的是爆炸建模而地面条件参与扩散与液池不参与爆炸。把输入挂错环节是环境类错误中最常见的一类。1.3 组合爆炸的算术气象工况数 ∏(各维度取值数)。工程上必须先把这张乘积算出来再决定策略全因子所有组合都跑适合筛选与达标论证一次一因子固定基准、逐维变动适合敏感性归因。官方在 Safeti™ 侧使用sensitivity studies敏感性研究这一名称Phast™ 侧 9.0 引入的批量执行能力名为Batch runner——两者不要混用也不要生造别的功能名。1.4 环境输入的版本与证据纪律环境参数是最容易被「顺手抄」的一类输入上一份报告里的风速、别人模板里的稳定度、界面上残留的默认值都可能在你不知情的情况下进入工况。三条纪律把这类风险压到最低。第一条每个气象取值都要有出处。现场气象年、法规规定的设计工况或官方技术文档三者必有其一并在工况表里以来源列固定查不到出处的取值一律以官方文档为准不得凭经验填写铁律 8。第二条风速与稳定度必须成对登记。二者在物理上耦合只记风速的工况表在复现时无法确定湍流强度因此工况表的主键应包含稳定度去重也按「风速 稳定度」这一对来判断。第三条地面条件要区分「扩散边界」与「液池边界」。同一块地即使物理上只有一种形态在数据模型里也要分别声明扩散侧的地面/基底条件与液池侧的表面类型否则改了一处以为改了全部结果只动了半边。此外建筑尾流与阻塞区属几何与场景级设置它们的口径随版本演进较多落地前应回到官方 User Manual 与发布说明核对不要沿用旧项目的配置假设。二、完整操作与脚本逐段剖析2.1 第一段构造气象与地面实体说明以下实体类名、枚举名与导入路径取自官方 PyPWS 参考文档属性名以官方 PyPWS 参考文档为准铁律 2标注[核对]的行请在落地前按官方字段表确认名称与单位。# -*- coding: utf-8 -*-# 演示构造一个气象对象与一个地面基底对象frompypws.entitiesimportWeather,Substrate# 官方实体气象、基底frompypws.enumsimportAtmosphericStabilityClass# 稳定度枚举frompypws.enumsimportSurfaceType,PoolSurfaceType# 地面类型 / 液池表面类型frompypws.utilitiesimportgetAnalyticsApiTarget# 打印当前 API 目标print(Running against: %s%(getAnalyticsApiTarget()))# 先固定运行目标便于审计weatherWeather()# 声明气象对象weather.windSpeed5.0# [核对] 风速单位 m/sSIweather.stabilityClassAtmosphericStabilityClass.D# [核对] 大气稳定度类别示例取 D 类# 说明稳定度取值集合以官方枚举定义为准下方逐项打印用于确认枚举可用取值substrateSubstrate()# 声明地面/基底对象substrate.surfaceTypeSurfaceType.LAND# [核对] 地表类型示例取陆地substrate.poolSurfaceTypePoolSurfaceType.LAND# [核对] 液池所处表面类型substrate.surfaceRoughness0.03# [核对] 地表粗糙度单位 m按官方字段表核对量纲# 枚举自检确认当前安装版本提供了预期的枚举成员防止跨版本名称漂移members[mformindir(AtmosphericStabilityClass)ifm.isupper()]# 只取大写成员print(AtmosphericStabilityClass 成员,members)# 打印后填入你自己的取值白名单print(SurfaceType 成员,[mformindir(SurfaceType)ifm.isupper()])print(PoolSurfaceType 成员,[mformindir(PoolSurfaceType)ifm.isupper()])这段代码的两个设计意图一是用枚举自检把官方枚举成员打印出来避免凭记忆写错稳定度类别二是把所有[核对]行集中让代码评审时一眼看到哪些属性名还需要回到官方文档确认。2.2 第二段生成气象工况矩阵笛卡尔积 去重# -*- coding: utf-8 -*-# 演示把气象维度展开为可复现的工况矩阵并写出审计用的 CSVimportcsv# 标准库写 CSVimportitertools# 标准库笛卡尔积frompypws.entitiesimportWeather,Substrate# 实体frompypws.enumsimportAtmosphericStabilityClass,SurfaceType# —— 维度定义每一维的取值都应有出处法规/项目约定/官方文档——WIND_SPEEDS[1.5,5.0,9.0]# 风速档位 m/sSTABILITIES[AtmosphericStabilityClass.D,# 稳定度档位示例两档AtmosphericStabilityClass.F]SURFACES[SurfaceType.LAND,SurfaceType.WATER]# 地面类型档位示例两档defbuild_weather(wind_speed,stability):按风速与稳定度构造 Weather 实体。weatherWeather()# 新建气象对象weather.windSpeedfloat(wind_speed)# [核对] 风速 m/sweather.stabilityClassstability# [核对] 稳定度枚举returnweather# 返回可直接消费的气象对象defbuild_substrate(surface):按地表类型构造 Substrate 实体。substrateSubstrate()# 新建基底对象substrate.surfaceTypesurface# [核对] 地表类型returnsubstrate# 返回基底对象# —— 展开笛卡尔积并给出显式的工况标识 ——cases[]# 工况清单容器forindex,(wind_speed,stability,surface)inenumerate(itertools.product(WIND_SPEEDS,STABILITIES,SURFACES),start1):case_idW%03d%index# 唯一工况号幂等与续跑的基础铁律 6cases.append({case_id:case_id,# 工况号weather:build_weather(wind_speed,stability),# 气象对象substrate:build_substrate(surface),# 基底对象wind_speed:wind_speed,# 原始维度值便于审计stability:stability.name,# 枚举名便于落盘surface:surface.name,# 枚举名便于落盘})# —— 组合数自检乘积必须等于工况条数否则说明生成逻辑有漏铁律 6——expectedlen(WIND_SPEEDS)*len(STABILITIES)*len(SURFACES)# 各维连乘assertlen(cases)expected,组合数不符%d vs %d%(len(cases),expected)# —— 去重自检以风速, 稳定度, 地表三元组为主键 ——keys[(c[wind_speed],c[stability],c[surface])forcincases]# 主键列表assertlen(keys)len(set(keys)),工况矩阵存在重复主键# —— 落盘审计 CSV只写维度值不写不可序列化的对象 ——withopen(weather_cases.csv,w,newline,encodingutf-8)asfh:writercsv.DictWriter(fh,fieldnames[case_id,wind_speed,stability,surface])writer.writeheader()# 表头forcaseincases:# 逐工况写行writer.writerow({k:case[k]forkinwriter.fieldnames})# 只取标量字段print(生成工况数,len(cases))# 人工确认规模要点工况号case_id由生成顺序确定因此只要维度定义不变工况号就不变——这是断点续跑与结果对账的前提。2.3 第三段把气象工况接到泄漏计算并在批量中轮换# -*- coding: utf-8 -*-# 演示把气象/地面工况轮换注入泄漏计算收集排放记录首条速率frompypws.calculationsimportVesselLeakCalculation# 官方容器泄漏计算frompypws.enumsimportResultCode# 结果码defrun_case(case,vessel,leak,calculation_factoryVesselLeakCalculation):对单个气象工况执行一次泄漏计算返回 (case_id, 结果码, 首条记录)。 设计要点 1) 工厂函数注入便于单元测试时替换为假对象 2) 每次新建计算对象避免跨工况状态污染幂等要求。 calculationcalculation_factory()# 新建计算实例calculation.vesselvessel# 复用容器对象含上游相态calculation.leakleak# 复用泄漏对象calculation.weathercase[weather]# 注入本工况气象calculation.substratecase[substrate]# 注入本工况地面条件result_codecalculation.run()# 执行计算first_recordNone# 首条排放记录占位ifresult_codeResultCode.SUCCESS:# 仅在成功时取数recordscalculation.dischargeRecords# 排放记录列表first_recordrecords[0]ifrecordselseNone# 取首条稳态恒有 2 条else:calculation.print_messages()# 失败时打印诊断消息returncase[case_id],result_code,first_record# 返回结构化结果results[]# 结果容器forcaseincases:# 遍历上一段的工况矩阵case_id,code,recordrun_case(case,vessel,leak)# 执行单工况results.append({# 记录结构化结果case_id:case_id,result_code:code.name,# 结果码名称便于分组统计first_record:str(record),# 字段以官方参考章为准})failed[rforrinresultsifr[result_code]!SUCCESS]# 失败隔离print(总工况%d失败%d%(len(results),len(failed)))# 批量健康度一眼可见foriteminfailed:# 打印失败清单print(失败工况,item[case_id],item[result_code])# 便于针对性重跑工程要点失败工况不抛出中断整个批而是收集到failed清单配合 2.2 段的稳定case_id即可实现「只重跑失败工况」的断点续跑。2.4 建筑与阻塞区的处理位置建筑尾流模型与室内释放属于场景级环境设置地图上的阻塞区属于几何输入用于爆炸建模。二者在桌面侧配置以官方 User Manual 为准在 API 侧则以场景输入传递。落脚本时的纪律是几何入几何文件/配置不要塞进气象对象否则后续维护无法追溯哪一层改了。三、常见报错与排查现象 1枚举成员名写错脚本在构造处即抛AttributeError。根因跨版本枚举命名漂移或凭记忆写了字母类别而官方枚举另有命名体例。解法如 2.1 段那样用dir()打印实际成员并建立本地白名单把白名单与版本号一起入库。现象 2改风速后结果几乎不动。根因稳定度与风速在物理上耦合只改风速而稳定度保持在极端类别湍流参数被稳定度主导。解法把风速与稳定度成对设计为矩阵2.2 段而不是单维扫风速。现象 3液池蒸发/池火结果与现场经验差异大。根因Substrate与实际地表不符例如硬化地面按草地粗糙度、水面场景按陆地。解法先确认地表类型再确认液池所处表面类型最后核对粗糙度量纲。现象 4爆炸结果每次都一样改气象无反应。根因把气象当成爆炸的输入。阻塞区是爆炸建模的几何输入气象主要作用于扩散与火焰形态。解法修爆炸场景时改几何阻塞区与点火相关设置不要靠调气象。现象 5工况数远超预期批跑时间失控。根因维度连乘未预估3×6×2 一类组合在叠加孔径与物料维度后迅速上千。解法在生成阶段打印乘积并断言上限对筛选阶段改用一次一因子最终论证阶段再上全因子。现象 6同一气象在不同机器上结果不同。根因环境工况只存在于界面里没有随研究一起归档换机器后稳定度或粗糙度落回了默认值。解法把环境工况表含来源列与去重主键与研究文件一起纳入版本管理任何一次结果都标注所用的工况表版本。四、动手练习在 2.2 段的维度定义里加入第三个稳定度档位验证组合数与weather_cases.csv行数同步增长且已有case_id不发生错位。为run_case()写一个假计算对象返回固定ResultCode验证失败隔离逻辑在没有真实服务时也能被执行。挑一个你项目的现场气象按官方 User Manual 的口径把风速与稳定度落到枚举值并说明你选该稳定度的依据出处铁律 8。五、小结与下一篇预告本篇把环境输入整理为「气象 / 地面 / 建筑与阻塞区」三层明确了Weather、Substrate与三个地面相关枚举的位置给出了气象工况矩阵的生成、去重、审计与失败隔离写法并守住了一条容易做错的边界阻塞区服务爆炸建模不服务扩散。下一篇08 第一个批处理将进入桌面侧的官方批量路径Excel import快速建研究与Batch runner队列执行把本篇的工况矩阵落成可跑的研究集合。本篇认知问题回显FAQQ1一个Weather由哪些维度构成哪些是自由度、哪些是派生量答自由度是风速、风向与大气条件以稳定度类别表达等可独立设定的量它们在脚本侧落在Weather实体上地面/基底相关量则落在独立的Substrate实体上。派生量是由这些自由度经过扩散模型计算得到的湍流参数与浓度场——它们不写在环境对象上而是作为结果出现。完整的属性清单与期望单位以官方 PyPWS 参考文档为准且一律 SI 单位。Q2AtmosphericStabilityClass与风速是什么关系答它是把大气条件离散化为若干类别的枚举。稳定度与风速在物理上耦合同一风速下不同稳定度对应不同湍流强度因此工程上应把二者作为成对维度设计成矩阵而不是只扫风速。具体类别集合以官方枚举定义与官方技术文档为准。Q3地面条件为什么要拆成Substrate、SurfaceType、PoolSurfaceType三层答因为三者回答不同问题。Substrate是地面/基底对象承载地表类型与粗糙度等近地面边界特征作用于扩散SurfaceType表达地表介质如陆地、水面PoolSurfaceType表达液池所处的表面类型作用于液池扩展与蒸发。拆层是为了让「扩散的近地面边界」与「液池的热/质量交换边界」可分别设定与审计。Q4建筑尾流/室内释放在什么场景下必须打开影响哪一环答当泄漏发生在有显著建构筑物遮挡或处于室内的场景时尾流模型与室内释放设置会改变扩散的初始与近场条件必须在场景级环境设置中打开否则近场浓度会被系统性低估。这些设置属桌面侧场景配置入口与参数以官方 User Manual 为准。Q5阻塞区为什么是爆炸建模的输入而不是扩散的输入答阻塞区congested regions描述的是设备/结构密集造成的空间约束其物理作用是影响火焰加速与超压生成因此服务爆炸建模扩散侧的边界由地面条件与气象承担。二者不该互相替代这也是批处理中最常见的挂错环节错误。Q6把气象做成批量维度时怎么避免组合爆炸与重复工况答三步。第一生成前先算乘积并设上限断言第二用去重主键风速、稳定度、地表等标量组合校验唯一性第三用稳定且与生成顺序绑定的工况号如W001作为幂等键使得维度定义不变时工况号不变从而支持断点续跑与结果对账。策略上筛选阶段用一次一因子论证阶段再用全因子。