SWMM模型是由美国环保署(EPA)资助开发的一个免费开源的城市雨洪管理模型，包括水量、水质、LID等多个模块[1]，已被世界各国的研究者广泛用于研究雨水问题，在中国也被广泛应用于海绵城市LID评估、排水管网复核、城市面源污染量估算等，并取得一系列的研究成果[2-3]。

然而，因为目前我国管网信息化程度较低，绝大部分管网数据以DWG格式存储，SWMM无法直接处理。若使用手工描绘方式进行建模，效率低且误差很大。所以，SWMM模型应用的核心问题是如何基于已有数据进行快速建模，通过计算机辅助建模可以大大提高建模效率，从而对SWMM模型的应用和推广有重要意义。

1　SWMM模型及文件结构

SWMM将城市排水管网系统中的水文和水力要素概化为管渠(Link)、节点(Node)和汇水区(Catchment)三种类型，用非线性水库模型模拟地表径流，用圣维南方程演算管网的输送过程，用累积冲刷模型模拟地表径流的污染[4]。SWMM模型主要用于城市区域降雨径流、合流制管道、污水管道和其他排水系统的规划设计、情景分析和方案评估等多个方面。资料翔实时也可应用于非城市区域的分析与模拟[5-6]。

2.1.1 正交试验 单体聚合反应中影响较大的因素有单体配比，温度，引发剂加量和pH值等。根据自由聚合原理，设定聚合反应的交联单体浓度均为30%，构建四因素四水平正交表分析缓凝剂性能影响因素的主次顺序[12,13]，正交表见表1。

SWMM的模型文件为INP格式文件，具有按块组织的特定文件格式。每一块都是从一个由中括号括着的关键词开始，关键词后为属性区，模型的各种参数反映了对应的管渠、雨量计、汇水区等的属性及相互关系，部分关键词如下。

[TITLE]工程标题

[OPTIONS]分析选项

因为提取过程需一次性完成，所以将步骤3融入到步骤2中从而得到完整的信息提取功能，之后再融入进步骤1的拓扑建立过程之中，即可在拓扑每增加一根管线同时，获得需要的管线、节点的各种信息。

[FILES]接口文件选项

图1为L2线自然电场法剖面曲线，可见3个频率的电位在140～250号点与400～600号点处均较低，推断该处裂隙较发育。

标高的提取流程如图4所示。

NOS1AP基因rs12742393位点多态性对瑞格列奈治疗我国汉族2型糖尿病患者疗效的影响 …………… 王 涛等（10）：1347

[HYDROGRAPHS]用于构造RDII进流量的单位水文过程线数据

[EVAPORATION]蒸发数据

本文选取河南省1989—2015年城镇化率、老龄人口抚养比和卫生机构数作为原始数据，运用GM(1,1)模型对河南省2016—2020年的城镇化率、老龄人口抚养比和卫生机构数进行预测，将输出结果作为多元线性回归预测方程的输入来计算河南省人口老龄化系数，实现灰色预测模型与多元线性回归方程的有机结合，最终取得较高精度的预测结果。

[TEMPERATURE]气温和雪融数据

招商部门在商场租赁谈判时，可以与租户进行沟通，在原商务谈判的租金总额、租赁使用面积不变的前提下，在合同中增加该铺位对应分摊的建筑面积，将租金按照建筑面积进行折算单位租金；以达到降低自用房建筑面积和原值、从而降低自用房房产税税额的目的。

煤矿安全生产过程中煤矿地质测量指的是对矿区的地形情况以及矿区地面的控制网进行分析和测绘，测量煤矿施工过程中相关地形地质状况，获取煤矿施工过程中矿区的地标移动和升降情况，保证煤矿安全生产过程中煤矿地质测量工作有效完成。以下是我国煤矿的主要分布图。

SWMM的INP文件中各节次序任意，并非所有节是必须的。分号(；)用于说明本行中，其后为注释而不是数据。数据项可以显示在行内任意列。各节关键词可以显示为混合小写和大写方式，仅仅前四个字符(加上起始中括号)用于区别关键词。

开展探究式实验的性别、年级和GPA并不影响学生对Q2和Q3的回答,但是女生比男生更相信探究式实验有助于强化她们对环境基本概念和应用的理解。对于Q4,三学年132位所调查的学生中,压倒性的92.4%学生是第一次接触探究式实验。仅有20个学生(占7.6%)以前接触过类似的探究式实验,这些不足8%的学生可能是由于外部课程(基础学科实验室或高中[12])而接触探究式实验。

2　建模实现过程

建模是进行SWMM模拟的第一步，基于VisualLISP快速提取管线信息并导入SWMM的过程可分为两部分，一是管线信息的提取，二是信息格式的转换与导出。

(2) 换乘模式。西安地铁2号线南稍门站位于十字路口南侧，车站结构形式为单柱双跨结构;5号线下穿友谊路,并下穿2号线区间。从车站改造难度及对地铁运营的影响考虑，5号线与2号线换乘首选通道模式，即通过在2号线车站站厅层侧墙上开洞，引出换乘通道与5号线车站进行换乘。

2.1　管线信息的提取

管线需提取的信息主要有以下内容：管径、管长、管线两端管内底标高、管线折点坐标、管线上下游与节点的对应关系、节点的坐标、节点内底高程、节点地面标高。为了提取到这些信息，有三个核心问题需要解决，一是建立管线的上下游拓扑关系，二是获取标高信息并与管线及节点对应，三是管径信息与管线的对应[7-8]。

2.1.1　建立管线的上下游拓扑关系

为了确定管线的上下游，首先自定义一个名为pck的块，需要用户通过点击下游管线末端插入到管线系统中，管线系统排出口位置即为块的基点。拓扑关系的建立流程如图1所示。

  

图1　管线拓扑关系建立流程图

 

Fig.1　Flow Chart of Pipelines Topology Establishment

如图2所示，管线拓扑建立的源代码，其中子程序“getnewaddnearobjss”承担根据搜索点表获取新增管线，并记录信息，同时更新搜索点表的任务；rcgx为管线端点的容差；gxename为管线的线形；gxlayer为管线图层。

2.1.5　可视化排错功能

  

图2　主程序源代码

 

Fig.2　Main Program Source Code

2.1.2　获取标高信息并与管线及节点对应

为了使提取后的标高管径信息能与管线及节点对应，需要合理的数据存储格式。这里使用表来存储数据，节点和管线的数据形式如下。

(节点坐标　地面标高　井底标高　节点编号)

(管线图元名　上游节点编号　下游节点编号上游管内底标高　下游管内底标高　管径　管线坐标　管线编号)

其中节点坐标和管线图元名可用来判断节点或管线是否已存在。这里讨论的标高标注形式为浙江样式，如图3所示。

  

图3　井管标注样式示意图

 

Fig.3　Sketch Map of Wells and Pipes Dimstyle

[RAINGAGES]雨量计信息

  

图4　管线特征标高识别及储存流程图

 

Fig.4　Flow Chart of Identification and Storage for Characteristic Elevation of Pipelines

需注意管线坐标保存时，应以下游节点坐标替换管线的下游坐标，以保证管线信息在提取后可以首尾相连，不受容差的影响；需使管线坐标按照上游至下游的顺序排列，以保证数据转换成SWMM格式后不至错乱。

三维管线高程信息可直接根据起终点Z值确定，可省略复杂的标注检索匹配流程[9]。

2.1.3　获取管径信息并与管线对应

服装企业设计、制造、销售所围绕的核心从根本上讲是服装产品，以产品数据管理作为主线可以将企业运营各方面数据贯穿为一个有机的整体。以往产品数据管理 （PDM）系统受制于行业需求以及供应商开发能力与服务水平，应用情况并不理想。产品数据管理 （PDM）系统在其他行业应用已成熟，有完善的实施过程且成效显著，然而服装行业整体所处的发展阶段，对PDM系统尚未产生足够需求，同时，考虑到PDM系统对于行业标准化的要求和业务流程适应性等因素，传统行业既有的成熟PDM系统很难嫁接应用于服装行业。

根据管线中点和设定的容差、图层等因素进行搜索，若搜索到多个标注则选择其中距离最近的那个。另外也可根据情况加入角度判断，保证标注和管线平行，降低出错几率。

水土流失造成土地退化，降低土地生产力，制约地区经济发展；加剧洪涝、干旱、风沙等自然灾害，制约生态文明建设；造成河道和水库淤积。威胁下游的防洪安全；带来面源污染，造成水质污染，威胁下游饮水安全。

2.1.4　功能融合

[REPORT]输出报告指令

采用SPSS 23.0 统计学软件进行数据分析，计数资料采用卡方检验，多个样本率的两两比较采用卡方分割法，检验水准α=0.05。

为了可视化检索提取信息的准确性，增加了标注功能；分图层分颜色将管线的中点与其所对应的管内底标高文字基点、管径文字基点连线，即可清楚地看出信息对应关系是否正确。同时标注管线及节点的编号，通过编号及信息连接线的标注进度也可知晓程序的运行情况。

2.2　数据的转换与导出

2.2.1　数据导出

将已有节点信息表和已有管线信息表，通过foreach函数导出并增加“ ”作为空格，即可得到规则的信息表格。图5为CAD原数据样式，表1和表2分别为提取出的节点和管段数据。

  

图5　CAD管线数据

 

Fig.5　CAD Pipelines Data

 

表1　节点信息提取

 

Tab.1　Node Information Extraction

  

节点信息(部分)坐标地面标高/m井底标高/m节点编号(5242900320176e+00600)211831951817(5242340320173e+00600)211811918616(5242610320175e+00600)20931934215(5242310320175e+00600)20911917414(5241890320181e+00600)207331911613(5241990320179e+00600)207311908512

 

表2　管段数据提取

 

Tab.2　Data Extraction from Pipe Sections

  

管段数据(部分)管线长度/m上游管内底标高/m下游管内底标高/m管径⁃长度/m⁃坡度管线节点坐标管线编号2934681951819342D400⁃293⁃60‰((5242900320176e+00600)(5242610320175e+00600))16121918619162D400⁃120⁃20‰((5242340320173e+00600)(5242320320175e+00600))15301934219162D400⁃300⁃60‰((5242610320175e+00600)(5242320320175e+00600))1461917419162D400⁃300⁃60‰((5242310320175e+00600)(5242320320175e+00600))13261911619077D400⁃260⁃15‰((5241890320181e+00600)(5241930320179e+00600))12

注：信息表中的坐标，为一或多对(x,y)形式的数据，对于有拐点的多段线同样适用

2.2.2　数据转换

将数据转换成为SWMM标准INP格式，是通过对已有节点信息表和已有管线信息表的操作来实现的。在数据输出时，分为[OUTFALLS]、[JUNCTIONS]、[CONDUITS]、[XSECTIONS]、[COORDINATES]、[VERTICES]6个部分。

  

图6　提取后自动生成INP文件

 

Fig.6　Automatically Generation of INP File after Extraction

通过nth、assoc等函数组合使用获取表中指定信息，经一定的对应、计算等处理后输出。SWMM中管内底标高是通过相对节点底部标高的偏移来实现的，所以需在获取管线内底标高后，根据其节点编号搜索对应的节点底部标高，经过计算后再输出，效果如图6所示。

3　应用案例

选取某公园来进行模型验证，管网基础数据为DWG格式管网规划图，汇水区基础数据为shp格式用地规划图，包括降雨、洪水位、蒸发等多种基础数据。根据本文上述建模过程所编写的CAD插件进行管道数据提取，成功提取为INP。汇水区由开源软件INPPINS转换为INP格式，将两个INP进行关键词对应即可合并，建立完整SWMM模型。输入模型参数后即可执行模拟，效果如图7所示。

分别输入芝加哥设计降雨数据和全年分钟降雨数据进行模拟，结果显示该公园规划海绵方案的年径流总量控制率在91.65%，达到90%的规划要求；规划排水管渠方案在P=5 a的情境下，无节点溢流，满足管网排水能力要求。本文快速提取方法大大降低了建模难度，提升了建模效率。

  

图7　管网及最终SWMM模型

 

Fig.7　Pipe Network and Final SWMM Model

4　结论

本文实现了将一定格式下的CAD管网数据快速提取的功能，可自动转换到SWMM模型中，大大减少建模人员的工作量，在规划、设计以及海绵城市建设效果评估中皆可得到应用。某公园的案例验证了方法的可行性。

本文所述方法仍有发展空间，如很多图纸并不是每个节点都有高程标注，缺失的部分仍需要手工来补充，后续可增加根据上下游标注自动补全标高的功能；另外由于管径标注形式多样，目前仅将标注整体提取，没有实现自动识别所需信息的功能；此外笔者开发了将汇水分区导入SWMM的程序，下一步有望将两个程序整合，实现更完整的功能，以期在海绵城市规划及设计评估等方面有更好的应用。

参考文献

[1]Zoppou C.Review of urban storm water models[J].Environmental Modelling & Software,2001,16(3):195-231.

[2]胡爱兵,任心欣,丁年,等.基于SWMM的深圳市某区域LID设施布局与优化[J].中国给水排水,2015,31(21):96-100.

[3]https://www.epa.gov/water-research/storm-water-management-model-SWMM[EB/OL].

[4]危唯.低影响开发技术在深圳某地区的应用研究[D].长沙:湖南大学,2014.

[5]吴亚男,熊家晴,任心欣,等.深圳鹅颈水流域SWMM模型参数敏感性分析及率定研究[J].给水排水,2015,41(11):126-131.

[6]任心欣,汤伟真,李建宁,等.水文模型法辅助低影响开发方案设计案例探讨[J].中国给水排水,2016,32(17):109-114.

[7]杨伟明,刘子龙,周玉文,等.基于CADTableConvert和雨水管网设计计算表的自动SWMM水力模型构建方法研究[J].给水排水,2016,42(4):124-127.

[8]张灵敏.排水管网水力计算及暴雨积水模拟方法研究[D].广州：华南理工大学,2015.

[9]陈淑珍.基于GIS和SWMM的城市排水管网建模方法研究[D].北京：北京工业大学,2014.