python元胞自动机模拟交通_文献分享 | 基于WCA2D与SWMM模型的城市暴雨洪涝快速模拟...

最新推荐文章于 2024-02-22 16:14:16 发布
最新推荐文章于 2024-02-22 16:14:16 发布
阅读量1.2k 收藏 2
点赞数 1
文章标签： python元胞自动机模拟交通

基于WCA2D与SWMM模型的城市暴雨洪涝

快速模拟

文献分享

第一作者：曾照洋 通讯作者：王兆礼 其他作者：赖成光、吴旭树、黄国如、胡庆芳 通讯单位：华南理工大学 土木与交通学院 亚热带建筑科学国家重点实验室 ff6d800cf49ecbeecbe2b8473b9c5b47.png

22ca84af40a128e120fd68c4220a11d1.png

[文章摘要] 为探讨计算高效的元胞自动机模型(WCA2D)与传统一维管网模型耦合的机制以及计算效果，尝试将WCA2D与SWMM模型耦合(SWMM/WCA2D)，以广州市长湴片区为例探究一种暴雨洪涝快速二维模拟技术，对比实测积水数据以及SWMM/LISFLOOD-FP模拟结果，得出如下结论：SWMM/WCA2D模拟结果与“20180607”实测积水数据相近，表明模型具有良好精度；通过与SWMM耦合，WCA2D和LISFLOOD-FP的模拟结果差异小(最大水深差值基本低于0.1m)、相关性强(相关系数基本超过0.7)，但前者计算效率是后者的3-5倍，表明WCA2D能够耦合SWMM且计算效率更高，为复杂城市化地区暴雨洪涝快速模拟提供了一种新方法。 [正文简述] 本研究区位于广州市天河区西北部的长湴片区(图1)，总面积约为1.57km²。该片区下垫面较为复杂，北部为山区，南部则为建筑比较密集的高度城市化地区。当暴雨发生时，大量山区的水流汇入长湴主干沟渠，容易造成该片区排水能力不足，从而导致地铁站附近出现积水。如在2018年6月8日，受长历时暴雨影响，天源路段最大水深将近1m，积水几乎没过地铁入口，严重威胁地下铁道系统的安全。本文以复杂城市化的广州市天河区长湴片区为例，构建SWMM/WCA2D模型，对2018年6月7日的场次暴雨(20180607)进行模拟，并采用实测积水数据进行验证，同时与SWMM/LISFLOOD-FP模型进行对比，旨在验证WCA2D的二维模拟效果以及模型耦合的可行性，以期为城市洪涝模拟、预警预报以及防涝减灾等工作提供一种新方法。

8b6b5091c5692609837f375385e1d3a4.png

图1 研究区位置与DEM8图

0843c291507b136ad738ac1f4293da39.png

图2 考虑建筑物的DEM5局部放大图 耦合WCA2D与SWMM模型实现城市暴雨洪涝的快速二维模拟，主要步骤如下：  (1)构建SWMM模型模拟降雨、地表产流、地面径流、管网汇流以及管网溢流过程；  (2)基于C# Winform平台，调用SWMM模型计算引擎的接口函数，提取节点溢流过程并保存为WCA2D模型所需的输入文件格式；  (3)生成模型配置文件，并链接地形数据与溢流文件；  (4)在C# Winform窗口驱动WCA2D模型计算引擎，模型通过读取配置文件并将模型参数、溢流过程、地形数据作为模拟条件，完成地表积水演进过程模拟；  (5)输出二维模拟结果文件，与SWMM/LISFLOOD-FP模型结果以及实测结果进行对比验证，并用ArcGIS与Python进行后处理。 根据是否考虑主干沟渠淤积深度以及二维模拟是否考虑建筑物对水流阻挡，再结合WCA2D与LISFLOOD-FP两种模型，共设置8种模拟情景(表1)。通过将主干沟渠的设计尺寸与现状断面对比，确定淤积深度约20cm；为模拟建筑物对地面积水的阻挡，将建筑物所在位置的高程进行加高处理(图2)。 表1 模拟情景汇总 0add0f3e757ccc4fe845b4ec003d8838.png 二维模拟的最大淹没水深和淹没范围如图3所示。结果表明，所有情景的积水范围主要分布于天源路以及长湴工业区，与实际调研结果相符，表明模拟结果能够很好的反映实际内涝情况。鉴于CA_B和LIS_B情景与下垫面的实际情况最接近，而积水深度普遍介于0.2-0.5m，与实际积水深度最为接近。因此，SWMM/WCA2D模拟结果与实测数据相近，且与SWMM/LISFLOOD-FP模拟结果相似，表明耦合模型具有较好精度。   640?wx_fmt=png 图3 所有情景最大积水深度空间分布图 分根据SWMM/WCA2D以及SWMM/LISFLOOD-FP模拟的积水深度相关系数(图4a)、积水面积(深度>0.05m，图4b)变化、计算效率(表2)等结果，可看出：随着降雨历时的增加，积水深度的相关系数先快速增大(至0.9以上)，在雨峰附近开始减小，最终趋于相对稳定(LIS_B与CA_B的相关系数接近0.9)；在2300min前，积水面积接近0，随着后雨峰影响，积水面积快速增大，而两种二维模型模拟的积水面积差异均比较小；从模拟时长角度分析， WCA2D模拟用时介于38-50s，而LISFLOOD-FP模拟用时达到170-200s，前者模拟效率是后者3-5倍。

23985083bfeee8fd5cdbe93f5b4ccd9f.png

(a)相关系数

426d3e77b8e4d6197aecf24734be0e1a.png

(b)积水面积 图4 二维模拟的积水深度相关系数变化以及积水面积变化图 表2 所有情景的模拟用时统计表(单位：秒) 36d4ad25eb2f4001a801737833c47e98.png 本研究获得国家重点研发计划资助项目(2018YFC1508205) ; 国家自然科学基金资助项目(51709117)的共同支持。 论文链接： http://journal16.magtechjournal.com/Jweb_skxjz/CN/10.14042/j.cnki.32.1309.2020.01.004 引用格式： 曾照洋, 赖成光, 王兆礼, 吴旭树, 黄国如, 胡庆芳. 基于WCA2D与SWMM模型的城市暴雨洪涝快速模拟[J]. 水科学进展, 2020, 31(1): 29-38.  ZENG Zhaoyang, LAI Chengguang, WANG Zhaoli, WU Xushu, HUANG Guoru, HU Qingfang. Rapid simulation of urban rainstorm flood based on WCA2D and SWMM model. Advances in Water Science, 2020, 31(1): 29-38. 编辑：邓梓锋 审核：王兆礼

点击 阅读原文 可进入文献网站！

weixin_39806779
关注
  • 1
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
【数学模型】基于Matlab实现洪水调度运算
matlab_dingdang的博客
10-05 2115
justw 水库调洪演算自编程序%水量平衡方程:(Q1+Q2)*dt/2-(q1+q2)*dt/2=V2-V1%水库蓄泄方程:q=m1*L*(z-z0)^1.5+m2*A*(z-z0)^0.5%zm--演算结果:[初始水位m 最高水位m 水位上涨m 初始库容1e4m3 最大库容1e4m3 调洪库容1e4m3 入库洪峰m3/s 最大泄量m3/s 洪峰削减m3/s]%Qt--洪水过程线:{[时间h 流量m3/s]}%Vz--水位库容曲线:[水位m 库容1e4m3]
python元胞自动机模拟交通_【数学建模】元胞自动机模拟交通流
weixin_39539002的博客
12-07 2038
元胞自动机(CA)是一种用来仿真局部规则和局部联系的方法。典型的元胞自动机是定义在网格上的，每一个点上的网格代表一个元胞与一种有限的状态。变化规则适用于每一个元胞并且同时进行。我们可以理解元胞自动机仿真需要理解三点。一是元胞，在matlab中可以理解为矩阵中的一点或多点组成的方形块，一般我们用矩阵中的一点代表一个元胞。二是变化规则，元胞的变化规则决定元胞下一刻的状态。三是元胞的状态，元胞...
元胞自动机模型代码matlab-SWMM-2DCA:此代号命名为雨水管理模型二维元胞自动机(SWMM-2DCA)，用于模拟城市排水系统与地表径
06-17
元细胞自动机模型代码matlab SWMM-2DCA 此代号命名为雨水管理模型二维元胞自动机 (SWMM-2DCA)，用于模拟城市排水系统与地表径流之间的洪水特性和流动相互作用。 代码是我的 MATLAB。 有关科学和技术信息，请参阅
CadFlood洪水分析软件使用手册
qq_36168428的博客
02-22 1017
计算河流洪水位是河湖划界、行洪论证等项目的关键工作之一。由于洪水分析专业性较强，目前用于洪水分析的软件较少，且一次仅能计算一条河流断面。在河湖划界等项目中，需计算数十、上百条河流断面洪水位，采用常规软件或手动计算，工作量十分繁琐。通过该软件进行洪水计算，可以提高计算准确性和减轻工作量，提高工作效率。
WCA-开源
04-27
1 / WCA解释一个文件，其中包含要分析的变量和表达式的定义，计算最坏情况分析并打印结果。 2 / WCADB扩展了MathCAD以根据从外部数据库获取的漂移来执行蒙特卡洛WCA。
lisflood-fp源码
08-20
lisflood-fp最新版源码，C++，推荐Visual Studio 2019打开，使用CMake组织编译
元胞自动机模拟交通路况.rar_交通元胞_交通模拟_元胞_元胞 交通_元胞自动机
07-15
元胞自动机，交通道路模拟，可修改参数，matlab源代码
元胞自动机模拟多晶长大的程序.rar_元胞自动机法_多晶_晶粒模拟_模拟晶粒长大_金属材料
07-14
元胞自动机（Cellular Automata，简称CA）是一种离散模型，广泛应用于复杂系统的研究，包括物理、生物、社会和计算机科学等多个领域。在金属材料科学中，元胞自动机法被用来模拟多晶材料的晶粒生长过程，这是一种...
元胞自动机模拟交通_matlab_元胞自动机_
10-01
交通元胞自动机模型通常包含以下几个关键要素： 1. **状态规则**：每个元胞根据其当前状态和相邻元胞的状态，按照一定的规则更新其下一时刻的状态。例如，在交通流中，车辆可能因前方有车而减速，也可能因道路空旷...
金属腐蚀的元胞自动机模拟代码 (1)_CA_元胞自动机_元胞_腐蚀_金属腐蚀元胞.rar
04-17
元胞自动机（Cellular Automata，简称CA）是一种离散模型，广泛应用于复杂系统的研究，包括物理、化学、生物等多个领域。在这个特定的场景中，我们关注的是它在金属腐蚀模拟中的应用。金属腐蚀是一个复杂的化学过程...
caesar-lisflood:CAESAR-Lisflood景观演化模型-开源
04-29
Caesar Lisflood是一个地貌/景观演化模型，将Lisflood-FP 2d水动力流模型（Bates等，2010）与CAESAR地貌模型相结合，模拟了河流集水区的侵蚀和沉积，并覆盖了从数小时到数千年的时间尺度。 这些sourceforge页面提供.exe的下载以及CAESAR-lisflood的源代码，并且可以在Wiki选项卡中找到说明和其他文档。 其中包括对CAESAR中使用的所有参数的完整说明，说明视频和其他信息。 可以在这里找到讨论和支持：https://groups.google.com/forum/#!forum/caesar-lisflood
Lisflood-FP 5软件包
最新发布
03-01
Lisflood-FP 5软件包 LISFLOOD-FP 欧洲联盟的欧洲环境局 2003年 流域尺度的洪水预测和损害评估模型，用于模拟洪水的形成、传播以及对环境和人类的影响。 地形数据、降雨数据、土地利用数据、土壤数据等。 洪水深度、洪水速度、洪水扩展范围、洪水对基础设施和人类活动的影响等。
flood-tiles:一个Web应用程序，用于在任何现代浏览器中动态模拟简单的淹没场景
05-12
洪水瓦 一个Web应用程序，用于在任何现代浏览器中动态模拟简单的淹没场景。 中HTML5画布层的实现在很大程度上受到启发，并从Mapbox的出色示例中被部分复制
lisflood-fp.zip
05-28
LISFLOOD-FP is a two-dimensional hydrodynamic model specifically designed to simulate floodplain inundation in a computationally efficient manner over complex topography. It is capable of simulating grids up to 106 cells for dynamic flood events and can take advantage of new sources of terrain information from remote sensing techniques such as airborne laser altimetry and satellite interferometric radar.
元胞自动机实现交通流模拟
04-24
包括换道系统，跟驰系统等功能。功能比较齐全，但是大部分程序都是加密的，只能拿来用，不能作为学习资源参考
元胞自动机.rar_交通流量_交通流量模型_元胞_元胞自动机_自动机
07-13
元胞自动机，可以用来建立交通流量模型，挺好用的。推荐下载
win10 运行lisflood_8踩过的坑
我的笔记
01-07 2795
1、安装cuda及VS2019 https://www.bilibili.com/read/cv7335974/ 2、netcdf.h头文件找不到问题 https://www.cnblogs.com/lyggqm/p/12926974.html https://blog.csdn.net/u011954296/article/details/54617048 3、编码问题 https://blog.csdn.net/qq_33485434/article/details/81196843 4、
swmm与lisflood-fp源码如何一起编译 CMake命令
我的笔记
08-09 1452
cmake_minimum_required(VERSION 3.13 FATAL_ERROR) include(CheckLanguage) project(lisflood-fp) set(_CONFIG config.default.cmake CACHE FILEPATH "CMake config file") include(${_CONFIG}) if (NOT CMAKE_BUILD_TYPE) set(CMAKE_BUILD_TYPE Release) endif() en.
Lisflood
警警的博客
01-26 815
文献：基于ＳＷＭＭ和ＬＩＳＦＬＯＯＤ－ＦＰ模型的居住区低影响开发模拟评价。问题：多个节点溢流同时，不行，待解决。目录下，cmd命令行。
元胞自动机与疏散模型的matlab算法,元胞自动机模拟多出口疏散模型的matlab实现...
05-29
关于元胞自动机与疏散模型的Matlab算法，可以参考以下步骤： 1. 定义元胞自动机的基本单位：每个元胞包含人员位置、速度、方向等属性。 2. 定义元胞自动机的演化规则：每个元胞根据周围元胞的状态，计算出自己的下一时刻状态，包括位置、速度、方向等。 3. 定义疏散模型的初始状态：包括人员位置、密度、出口位置等。 4. 基于元胞自动机的演化规则，模拟人员在疏散过程中的运动。 5. 根据疏散模型的结束条件，判断模拟是否结束。 6. 对模拟结果进行分析和可视化展示。 以下是一个简单的元胞自动机模拟多出口疏散模型的Matlab实现： ```matlab % 定义元胞自动机的参数 N = 200; % 元胞总数 L = 20; % 区域长度 W = 10; % 区域宽度 dt = 0.1; % 时间步长 T = 100; % 总模拟时间 v0 = 1.5; % 最大速度 tau = 0.5; % 安全时间 r = 0.3; % 感知半径 a = 2; % 舒适加速度 b = 1.5; % 减速加速度 % 定义疏散模型的参数 M = 4; % 出口数量 exit = [0, W/2; L/2, W; L, W/2; L/2, 0]; % 出口位置 density = 0.5; % 初始密度 N0 = round(N*density); % 初始人数 % 初始化元胞状态 pos = rand(N, 2)*[L, W]; % 位置 vel = zeros(N, 2); % 速度 dir = rand(N, 1)*2*pi; % 方向 exit_dist = zeros(N, M); % 距离出口的距离 exit_time = zeros(N, M); % 到达出口的时间 exit_chosen = zeros(N, 1); % 已选择出口的编号 for i = 1:N for j = 1:M exit_dist(i, j) = norm(pos(i, :)-exit(j, :)); exit_time(i, j) = exit_dist(i, j)/v0; end end % 模拟疏散过程 for t = 0:dt:T % 计算每个元胞的加速度 acc = zeros(N, 2); for i = 1:N % 感知周围元胞 nei = []; for j = 1:N if i ~= j && norm(pos(i, :)-pos(j, :)) < r nei = [nei, j]; end end % 计算期望速度和期望位置 v_des = exit_dist(i, exit_chosen(i))/exit_time(i, exit_chosen(i)); pos_des = pos(i, :) + v_des*exit(exit_chosen(i), :)/norm(exit(exit_chosen(i), :)); % 计算加速度 acc(i, :) = a*(v_des*exit(exit_chosen(i), :)/norm(exit(exit_chosen(i), :))-vel(i, :))/tau; for j = nei if norm(pos(i, :)-pos(j, :)) < r acc(i, :) = acc(i, :) + b*(vel(j, :)-vel(i, :))/(norm(pos(i, :)-pos(j, :))-r)^2; end end end % 更新每个元胞的状态 for i = 1:N vel(i, :) = vel(i, :) + acc(i, :)*dt; if norm(vel(i, :)) > v0 vel(i, :) = vel(i, :)*v0/norm(vel(i, :)); end pos(i, :) = pos(i, :) + vel(i, :)*dt; % 到达出口的处理 if exit_chosen(i) > 0 && norm(pos(i, :)-exit(exit_chosen(i), :)) < 0.1 vel(i, :) = [0, 0]; pos(i, :) = pos(i, :) + [L, W]; exit_chosen(i) = 0; end % 选择出口 if exit_chosen(i) == 0 min_time = Inf; for j = 1:M if exit_time(i, j) < min_time min_time = exit_time(i, j); exit_chosen(i) = j; end end end end % 可视化展示 clf; hold on; for i = 1:M plot(exit(i, 1), exit(i, 2), 'rx', 'MarkerSize', 10, 'LineWidth', 2); end for i = 1:N plot(pos(i, 1), pos(i, 2), 'bo', 'MarkerSize', 5, 'LineWidth', 1); end axis equal; axis([0, L, 0, W]); drawnow; end ``` 该代码实现了一个包含4个出口的疏散模型，模拟了人员在疏散过程中的运动，并可视化展示了模拟结果。其中，元胞自动机的演化规则基于人员的感知、期望速度和加速度，疏散模型的初始状态包括人员位置、密度和出口位置。模拟过程中，根据人员到达出口的时间选择最近的出口，当人员到达出口时，将其移动到对应区域的另一侧。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值