CN109583035A - 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 - Google Patents
基于元胞自动机的城市面源污染物累积冲刷过程计算方法 Download PDFInfo
- Publication number
- CN109583035A CN109583035A CN201811307557.4A CN201811307557A CN109583035A CN 109583035 A CN109583035 A CN 109583035A CN 201811307557 A CN201811307557 A CN 201811307557A CN 109583035 A CN109583035 A CN 109583035A
- Authority
- CN
- China
- Prior art keywords
- cellular
- water
- pollutant
- amount
- rainfall
- 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
- 230000001413 cellular effect Effects 0.000 title claims abstract description 190
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- 230000008569 process Effects 0.000 title claims abstract description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 145
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 110
- 231100000719 pollutant Toxicity 0.000 claims abstract description 109
- 238000012546 transfer Methods 0.000 claims abstract description 40
- 238000009825 accumulation Methods 0.000 claims abstract description 23
- 208000035126 Facies Diseases 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 230000001186 cumulative effect Effects 0.000 claims abstract description 4
- 239000013618 particulate matter Substances 0.000 claims description 16
- 230000008595 infiltration Effects 0.000 claims description 15
- 238000001764 infiltration Methods 0.000 claims description 15
- 238000011160 research Methods 0.000 claims description 8
- 230000003628 erosive effect Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 238000009991 scouring Methods 0.000 claims description 5
- 239000002245 particle Substances 0.000 claims description 4
- 208000004998 Abdominal Pain Diseases 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 206010013990 dysuria Diseases 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 2
- 239000013049 sediment Substances 0.000 claims description 2
- 239000003403 water pollutant Substances 0.000 claims description 2
- 238000004064 recycling Methods 0.000 claims 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 abstract 1
- 230000035508 accumulation Effects 0.000 description 18
- 238000010586 diagram Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000004080 punching Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 229920001617 Vinyon Polymers 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003818 cinder Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Classifications
-
- 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
-
- 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
- G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
- G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于元胞自动机的城市面源污染物累积冲刷过程计算方法,包含:构建元胞地理场景,获得每个元胞的下垫面类型、高程和随时间变化的降雨量;根据下垫面下渗参数,进行产流计算和流向计算,进而通过汇流计算获得每个元胞的当前水深;根据污染物累积和冲刷相关系数和干期天数,计算每个元胞的污染物累计量,并进行污染物输移过程的模拟计算,获得每个元胞的陆相污染物转移量和水相污染物输移量,进而计算水相污染物浓度。本发明结合水利学和污染物累计冲刷原理,具有明确的物理含义和理论基础,能够有效模拟污染物的累计冲刷过程,为城市面源污染累积冲刷模拟提供了新方法,为城市面源污染负荷削减和控制提供指导和依据。
Description
技术领域
本发明涉及环境科学领域,尤其涉及一种基于元胞自动机的城市面源污染物累积冲刷过程计算方法。
背景技术
随着城市化的快速发展,城市不透水面不断扩大,降雨径流冲刷汇流所形成的面源污染已成为城市水环境质量下降的主要原因。面源污染模型能够对污染物的迁移转化过程进行定量模拟,已成为城市面源负荷来源解析和治理的重要研究内容。目前城市面源污染模型在模拟污染物冲刷时大多采用黑箱方法或灰箱方法,只能获得研究区的整体冲刷量而不能反映局部冲刷特征,模型参数缺乏物理含义,对污染物冲刷过程的机理描述尚不充分。
元胞自动机是一种时间、空间、状态都离散的网格动力学模型,具有模拟复杂***时空演化过程的能力,尤其适合变量存在局部非连续或者剧烈梯度变化的情景,在城市复杂地形面源污染累积与冲刷过程模拟中具有独特的优势。
因此,充分考虑污染物冲刷时的降雨溅蚀作用、径流冲刷作用和颗粒物再沉积等物理作用机制的基础上,结合水力学和污染物累积冲刷原理,构建城市面源污染物累积冲刷的元胞自动机模型,实现基于元胞自动机的城市降雨径流和污染物冲刷过程时空动态模拟,对于城市面源污染负荷削减和控制具有重要的指导意义,这也本发明的主要研究课题。
发明内容
发明目的:本发明目的是提供一种基于元胞自动机的城市面源污染物累积冲刷过程计算方法,以更机理的描述污染物的累积冲刷过程,为面源污染累积冲刷过程模拟提供新方法。
技术方案:本发明所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,主要包括如下步骤:
(1)构建元胞地理场景模型,包括定义地理场景模型元胞大小、领域关系、模型边界和出口,根据获取的研究区的下垫面类型、DEM高程和随时间变化的降雨量,获得每个元胞的下垫面类型、高程和随时间变化的降雨量;
(2)根据设置的下垫面下渗参数利用初损后损法进行每个元胞的产流和下渗量的计算得到元胞的净雨深,并根据中心元胞与邻域元胞高程的水位差进行流向计算,进而通过汇流计算获得每个元胞的当前水深;
(3)根据设置的污染物累积和冲刷相关系数和干期天数,采用指数累积公式计算每个元胞的干期污染物累积量,并考虑污染物的降雨溅蚀、径流冲刷和再沉积作用,对污染物输移过程进行模拟计算,得到径流冲刷强度、降雨冲刷强度和污染物再沉降强度,从而获得每个元胞的陆相污染物转移量和水相污染物输移量,进而计算水相污染物浓度。
进一步的,考虑到模型计算的结果有效性,所述方法还包括(4)结合实测的流量和污染物数据,利用纳什效率系数(EF)和皮尔逊相关系数(CC)进行计算结果的评价。
进一步的,所述步骤(2)中利用初损后损法将下渗简化为初损和后损两个部分,降雨开始到出现产流时的水量损失为初损量(I0);产流后的水量损失为后损,一般用后损平均入渗率(f)来表示;通过设置初损量和平均后损入渗率计算净雨深(h),也即每个元胞的水深,计算公式如下:
h=X-I0-ftc-Xn
式中,h为X形成的元胞的水深(净雨深),mm;X为总降雨量,mm;I0为初损量,mm;f为平均后损入渗率,mm/min;tc为超渗产流历时,min;Xn为后损阶段非超渗历时的降雨量,mm。
进一步的,所述步骤(2)中方法进行水流流向的判断,具体为:比较中心元胞与邻域元胞高程的水位差,排除高于平均高程的邻域元胞,循环该过程直到没有高于平均高程的邻域元胞可以被排除,确定可以进行水量转移的邻域元胞。
进一步的,所述步骤(2)中汇流计算具体为:在产流和流向确定的基础上,采用曼宁公式计算一个时间步长内中心元胞可以转移的水量,同时根据水位差得出可以转移的潜在最大水量,将二者相比,得出一个时间步长内中心元胞转移的真实水量。计算公式如下:
式中,T为中心元胞水深转移至邻居元胞所需时间,既水量/污染物转移时间步长;V为流速,m/s;h为中心元胞水深,m;s为中心元胞与第i个邻居元胞间的水面坡度;n为曼宁糙率系数,s/m1/3;t为模型时间步长;D为中心元胞与邻居元胞i的距离,m,对于四邻域而言,D为元胞边长,对于斜向邻域元胞,D为倍的元胞边长;fi为水深,表示中心元胞向邻居元胞i转移的水量。
比较时间步长t与计算出的T的大小,计算中心元胞向邻居元胞i转移的水量,如果T小于时间步长t,认为水量fi可以全部转移;否则只能转移相应比例的水量。根据元胞平均高程与邻居元胞的高差计算最多可以转移至的水量,每个邻居元胞转移水量fi求和可以得出转移总水深。并将转移总水深与中心元胞水深进行对比,确定每个元胞的实际转移水量。将上述规则应用于整个元胞空间,计算出每一时刻整个元胞空间内所有元胞的转移水量和剩余水量,并根据下一时刻步长的降雨和下渗量,获得下一时刻末的元胞水量,循环计算,直到整个模拟过程结束。
进一步的,所述步骤(3)中干期污染物累积量计算公式如下:
式中,M为干期沉积层污染物累积量,单位面积或者边沿长度的质量;td为干期时长,天;C1为最大增长可能,单位面积或者边沿长度的质量,C2为增长速率常数,1/天。
进一步的,所述步骤(3)中按照如下公式进行污染物冲刷过程计算,得到每个元胞陆相污染物转移量:
di=vsici
其中h0=0.33DR
Ωe=F(Ω-Ω0),其中Ω=ρwgSfq
Δpi=(rri+eri-di)*T*L2
式中,mT为单位面积的沉积层颗粒物总质量,kg/m2;mi为单位面积的沉积层中第i组颗粒物质量,kg/m2;P为降雨强度,m/s;ρs和ρw分别为沉积物密度和水密度,kg/m3;g为重力加速度常数,m/s2;vsi为第i组颗粒物在水中沉降速度,m/s;ci为第i组颗粒物在水体的浓度,kg/m3;ad为沉积层颗粒物可分离性,kg/m3;ad0为沉积物初始可分离性,kg/m3;q为单宽流量,m2/s;h为元胞水深,m;h0为临界水深,m;b为常数;DR为雨滴平均直径,mm;Ωe为有效径流功率,W/m2;Ω为计算的总径流功率,W/m2;Ω0为临界径流功率,W/m2;Sf为摩阻比降;F为比例系数;rri为径流冲刷强度;eri为降雨冲刷强度;di为污染物再沉降强度;△pi为每个元胞陆相污染物转移量;T为水量转移时间步长,也即模型污染物模拟时间步长,s;L为元胞边长,m。
进一步的,所述步骤(3)中根据每个元胞初始污染物累计量M和陆相污染物转移量△pi,利用初始污染物累积量减去陆相污染物转移量即可得出此时刻每个元胞陆相污染物总量,初始水相污染物量加上陆相污染物转移量即可得到水相污染物总量。根据每个元胞的水相污染物量除以元胞总的水量,即可以求出每个元胞的水体污染物浓度。在水量转移的基础上,按照每个元胞此时刻的污染物浓度,进行污染物量转移计算,得出整个元胞空间每个元胞当前时刻末的污染物量和污染物浓度,循环计算,直到整个模拟过程结束。
有益效果:与现有技术相比,本发明综合考虑了降雨溅蚀作用、径流冲刷作用和颗粒物再沉积作用等物理机制,建立了一种面源污染物累积冲刷过程元胞自动机模型。该模型考虑了颗粒物粒径对冲刷的影响,使用了真实的时间步长,相较于传统模型,本研究提出的模型具有明确物理含义和理论基础,不仅能够准确模拟面源污染物的累积冲刷输移过程,而且能够模拟不同级配颗粒物的冲刷规律和集水区出口的颗粒物粒径组成。
附图说明
图1是本发明的一个实施例的流程示意图;
图2是本发明的元胞边界条件示意图;
图3是流向判断算法流程示意图;
图4是污染物冲刷计算流程示意图;
图5是流量验证结果示意图;
图6是总悬浮物TSS验证结果示意图。
具体实施方式
如图1所示,本发明实施例公开的一种基于元胞自动机的城市面源污染物累积冲刷过程,包括以下步骤:
(1)获取研究区的下垫面类型、数字高程数据(DEM)和随时间变化的降雨量,设置元胞大小、时间步长、模型边界和出口,建立元胞基础地理场景。
该步骤具体包括:
(1-1)获取研究区的下垫面类型和数字高程数据(DEM)构建元胞基础地理场景,获得每个元胞的下垫面类型(例如:道路和屋顶等)和高程,设置模型的边界和模型的出口,定义每个元胞的大小、曼宁糙率系数、坐标、邻域关系,模型的时间步长和模拟时长,相关参数如表1所示。如图2所示,本发明中元胞的邻域关系为Moore型八邻域,边界元胞设置为极大值从而保证水量无法进入,本发明暂不考虑城市管网的影响和蒸散发损失。
表1地理元胞模型参数集示例表
参数 | 意义 | 说明 |
CA_length | 元胞边长 | 模型设置 |
n | 曼宁糙率 | 下垫面材质平均糙率 |
T_stop | 模拟时长 | 根据降雨过程进行设置 |
t | 时间步长 | 设置模拟时间的步长 |
(1-2)利用雨量计测定或通过气象站点获得研究区的降雨量数据,利用Arcgis软件进行站点降雨量数据的时空插值,获得每个元胞不同时刻的降雨量数据。
(2)根据设置的下垫面下渗参数(初损量、后损平均入渗率等),进行产流计算和流向计算,进而通过汇流计算获得每个元胞的当前水深。
该步骤具体包括:
(2-1)产流计算:利用初损后损法将下渗简化为初损和后损两个部分,降雨开始到出现产流时的水量损失为初损量(I0);产流后的水量损失为后损,一般用后损平均入渗率(f)来表示;通过设置初损量和平均后损入渗率计算净雨深(h),也即每个元胞的水深,计算公式如下:
h=X-I0-ftc-Xn
式中,h为X形成的元胞的水深,也即净雨深,mm;X为总降雨量,mm;I0为初损量,mm;f为平均后损入渗率,mm/min;tc为超渗产流历时,min;Xn为后损阶段非超渗历时的降雨量,mm。
比较降雨量和初损量,当累积降雨量达到初损量时,元胞开始具备产流条件。根据模型时间步长,分别计算出当前时间步长内各个元胞的降雨量和后损入渗量,并比较大小。若降雨量大于入渗量,该元胞产流,且产流量为降雨量与入渗量的差值;反之,该元胞不产流。将此规则用于整个元胞空间,结合元胞的大小,即可得到所有元胞当前时刻的水深。
(2-2)流向判断:根据DEM和元胞水深进行水流流向的判断,确定中心元胞与邻域元胞的平均高程,比较邻域元胞与元胞平均高程,排除高于平均高程的邻域元胞,循环该过程,直到没有高于平均高程的元胞可以被排除,确定可以进行水量转移的元胞,流向判断算法流程图如图3所示。
(2-3)汇流计算:在产流和流向确定的基础上,采用曼宁公式计算一个时间步长内中心元胞可以转移的水量,同时根据水位差,得出可以转移的潜在最大水量,将二者相比,可得出一个时间步长内中心元胞转移的真实水量。
具体的,假设中心元胞的水共流入k个邻居元胞,用第i个表示在k个邻居元胞中的某一个。首先,根据曼宁公式计算中心元胞水深h转移至第i个邻居元胞所需时间T,(也既水量/污染物转移时间步长),计算公式如下:
其中,V为流速,m/s;h为中心元胞水深,m;s为中心元胞与第i个邻居元胞间的水面坡度;n为曼宁糙率系数,s/m1/3;D为中心元胞与邻居元胞i的距离,m,对于四邻域而言,D为元胞边长,对于斜向邻域元胞,D为倍的元胞边长。
进而,比较时间步长t与计算出的T的大小,计算中心元胞向邻居元胞i转移的水量(用水深fi表示)。如果T小于时间步长t,认为水量可以全部转移,否则只能转移相应比例的水量。水深fi使用如下公式计算:
根据元胞平均高程与邻居元胞i间的高差计算最多可以转移至元胞i的水量,如果高差大于fi,则水量转移为fi;若高差小于fi,则转移水量的量使得两元胞高程相同,没有水位差。每个邻居元胞转移水量fi求和可以得出转移总水深(f0),将其与中心元胞水深进行对比,若中心元胞水深h小于f0,即中心元胞没有足够的水进行转移,每个邻居元胞水量转移量需按照水深进行缩小,重新计算f0;若中心元胞水深h大于f0,转移水量不变。将此规则应用于整个元胞空间,计算出每一时刻整个元胞空间内所有元胞的转移水量和剩余水量。
通过以上计算可以得到当前时刻末每个元胞的水量,当前时刻末每个元胞的水量作为下一时刻初始时元胞的水量,根据下一时刻步长的降雨和下渗量,获得下一时刻末的元胞水量,循环计算,直到整个模拟过程结束。
(3)根据设置的污染物累积和冲刷相关系数和干期天数,计算每个元胞的污染物累计量,并进行污染物输移过程的模拟计算,获得每个元胞的陆相污染物转移量、水相污染物输移量和水相污染物浓度。
该步骤具体包括:
(3-1)污染物累积过程计算:设置污染物累积干期天数,干期天数可以通过气象局网站公布的每日降水量数据进行确定。采用应用较多的SWMM软件中的指数累积公式计算每个元胞的污染物累积量;其中,干期污染物累积量计算公式如下:
式中,M为沉积层污染物累积量,单位面积或者边沿长度的质量;t为干期时长,天;C1为最大增长可能,单位面积或者边沿长度的质量,C2为增长速率常数,1/天。C1和C2为经验系数,可通过文献获得。
(3-2)污染物冲刷过程计算:结合元胞自动机的转换规则进行污染物输移过程的模拟计算,计算径流冲刷强度rri、降雨冲刷强度eri和污染物再沉降强度di,进而计算每个元胞陆相污染物转移量(△pi),具体计算公式如下:
di=vsici
Δpi=(rri+eri-di)*T*L2
其中,
Ωe=F(Ω-Ω0),Ω=ρwgSfq
式中,mT为单位面积的沉积层颗粒物总质量,kg/m2,可由(3-1)中计算的污染物累积量进行计算;mi为单位面积的沉积层中第i组颗粒物质量,kg/m2,可利用粒径分析仪测量粒径组分比例测定,并乘以总的污染物量M来获得;P为降雨强度,m/s;ρs和ρw分别为沉积物密度和水密度,kg/m3,为经验参数;g为重力加速度常数,m/s2;vsi为第i组颗粒物在水中沉降速度,m/s,可采用Stokes公式进行计算;ci为第i组颗粒物在水体的浓度,kg/m3;ad为沉积层颗粒物可分离性,kg/m3;ad0为沉积物初始可分离性,kg/m3,为模型经验参数;q为单宽流量,m2/s;h为元胞当前时刻水深,m;h0为临界水深,m,为模型经验参数;b为常数,单位为1,为模型经验参数;DR为雨滴平均直径,mm;Ωe为有效径流功率,W/m2;Ω为计算的总径流功率,W/m2;Ω0为临界径流功率,W/m2,为模型经验参数;Sf为摩阻比降;F为一个比例系数,反映了有效径流功率与超过临界径流功率部分的总径流功率的比值;T为水量转移时间步长,也即模型污染物模拟时间步长,s;L为元胞边长,m。
根据计算获得每个元胞的陆相污染转移量△pi,判断△pi的正负,△pi若为正,说明累积污染物从陆相被冲刷入水中量大于沉降量,污染物净进入水体;△pi若为负,说明污染物沉降量大于被冲刷进入水中的量,污染物从水体输出。
根据每个元胞陆相污染物转移量△pi,得出此时刻每个元胞水相污染物总量P水i(初始水相污染物量加上陆相污染物转移量)和陆相污染物总量P沉i(初始陆相污染物量减去陆相污染物转移量);由水相污染物量除以元胞总的水量,可以求出每个元胞水体中污染物浓度。在水量转移的基础上,按照每个元胞此时刻的污染物浓度,进行污染物量转移计算,得出整个元胞空间每个元胞当前时刻末的污染物量和污染物浓度,循环计算,直到整个模拟过程结束,污染物冲刷计算过程如图4所示。
(4)结合实测的流量和污染物数据,利用纳什效率系数(EF)和皮尔逊相关系数(CC)进行计算结果的评价。
使用流量计进行降雨过程流量的测量,获得每个时刻模型出口的实测流量数据。污染物数据获取时,利用经过校正的聚乙烯塑料瓶在出水口收集全过程的径流,考虑到降雨的初期冲刷效应以及雨量的大小,前期水样收集的时间间隔应设置相对较短(例如:每1-2min收集一次),后期水样收集时间间隔可适当加大,每次降雨过程一般收集15-20瓶水样,记录每瓶水样收集满所用时间。按照采样顺序进行编号并在进行样品的实验室分析。水质检测分析方法按照参照国家标准(《水和废水监测分析方法(第四版)》)进行,总悬浮物(TSS)检测采用重量法。
根据获得实测的流量和TSS数据,计算纳什效率系数(EF)和皮尔逊相关系数(CC)评价模拟结果的可靠性和适用性。其中,EF和CC为已有成熟的计算公式,不作具体介绍。流量验证结果如图5所示,污染物验证结果如图6所示。
以上所揭露的仅为本发明一种较佳实施例而已,不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (8)
1.基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于,该方法包括:
(1)构建元胞地理场景模型,包括定义地理场景模型元胞大小、领域关系、模型边界和出口,根据研究区的下垫面类型、DEM高程和随时间变化的降雨量,获得每个元胞的下垫面类型、高程和随时间变化的降雨量;
(2)根据设置的下垫面下渗参数利用初损后损法进行每个元胞的产流和下渗量的计算得到元胞的净雨深,并根据中心元胞与邻域元胞高程的水位差进行流向计算,进而通过汇流计算获得每个元胞的当前水深;
(3)根据设置的污染物累积和冲刷相关系数和干期天数,采用指数累积公式计算干期污染物累积量,并考虑污染物的降雨溅蚀、径流冲刷和再沉积作用,对污染物输移过程进行模拟计算,得到径流冲刷强度、降雨冲刷强度和污染物再沉降强度,从而获得每个元胞的陆相污染物转移量和水相污染物输移量,进而计算水相污染物浓度。
2.根据权利要求1所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于:所述步骤(2)中利用初损后损法将下渗简化为初损和后损两个部分,降雨开始到出现产流时的水量损失为初损量I0;产流后的水量损失为后损,用后损平均入渗率f表示;每个元胞的净雨深计算公式如下:
h=X-I0-ftc-Xn
式中,h为X形成的元胞的水深,即净雨深;X为总降雨量;tc为超渗产流历时;Xn为后损阶段非超渗历时的降雨量。
3.根据权利要求1所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于:所述步骤(2)中水流流向判断方法,具体为:比较中心元胞与邻域元胞高程的水位差,排除高于平均高程的邻域元胞,循环该过程直到没有高于平均高程的邻域元胞可以被排除,确定可以进行水量转移的邻域元胞。
4.根据权利要求1所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于:所述步骤(2)中汇流计算具体为:在产流和流向确定的基础上,采用曼宁公式计算一个时间步长内中心元胞可以转移的水量,同时根据水位差得出可以转移的潜在最大水量,将二者相比,得出一个时间步长内中心元胞转移的真实水量;计算公式如下:
式中,T为中心元胞水深转移至邻居元胞所需时间,既水量/污染物转移时间步长;h为中心元胞水深;s为中心元胞与第i个邻居元胞间的水面坡度;n为曼宁糙率系数;t为模型时间步长;D为中心元胞与邻居元胞i的距离,对于四邻域而言,D为元胞边长,对于斜向邻域元胞,D为倍的元胞边长;fi表示中心元胞向邻居元胞i转移的水量;
每个邻居元胞转移水量fi求和可以得出转移总水深,将转移总水深与中心元胞水深进行对比,确定每个元胞的实际转移水量;计算出每一时刻整个元胞空间内所有元胞的转移水量和剩余水量,并根据下一时刻步长的降雨和下渗量,获得下一时刻末的元胞水量,循环计算,直到整个模拟过程结束。
5.根据权利要求1所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于:所述步骤(3)中干期污染物累积量计算公式如下:
式中,M为干期沉积层污染物累积量,单位面积或者边沿长度的质量;td为干期时长;C1为最大增长可能,单位面积或者边沿长度的质量,C2为增长速率常数。
6.根据权利要求1所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于:所述步骤(3)中按照如下公式进行污染物冲刷过程计算,得到每个元胞陆相污染物转移量:
di=vsici
其中h0=0.33DR
Ωe=F(Ω-Ω0),其中Ω=ρwgSfq
Δpi=(rri+eri-di)*T*L2
式中,mT为单位面积的沉积层颗粒物总质量;mi为单位面积的沉积层中第i组颗粒物质量;P为降雨强度;ρs和ρw分别为沉积物密度和水密度;g为重力加速度常数;vsi为第i组颗粒物在水中沉降速度;ci为第i组颗粒物在水体的浓度;ad为沉积层颗粒物可分离性;ad0为沉积物初始可分离性;q为单宽流量;h为元胞水深;h0为临界水深;b为常数;DR为雨滴平均直径;Ωe为有效径流功率;Ω为计算的总径流功率;Ω0为临界径流功率;Sf为摩阻比降;F为比例系数;rri为径流冲刷强度;eri为降雨冲刷强度;di为污染物再沉降强度;△pi为每个元胞陆相污染物转移量;T为水量转移时间步长,也即模型污染物模拟时间步长;L为元胞边长。
7.根据权利要求1所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于:所述步骤(3)中根据每个元胞初始污染物累计量和陆相污染物转移量,利用初始污染物累积量减去陆相污染物转移量得出每个元胞陆相污染物总量,初始水相污染物量加上陆相污染物转移量得到水相污染物总量;根据每个元胞的水相污染物量除以元胞总的水量,求出每个元胞的水体污染物浓度;在水量转移的基础上,按照每个元胞的污染物浓度,进行污染物量转移计算,得出整个元胞空间每个元胞的污染物量和污染物浓度,循环计算,直到整个模拟过程结束。
8.根据权利要求1所述的基于元胞自动机的城市面源污染物累积冲刷过程计算方法,其特征在于:所述方法还包括:步骤(4)结合实测的流量和污染物数据,利用纳什效率系数和皮尔逊相关系数进行模拟计算结果的评价。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811307557.4A CN109583035B (zh) | 2018-11-05 | 2018-11-05 | 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811307557.4A CN109583035B (zh) | 2018-11-05 | 2018-11-05 | 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109583035A true CN109583035A (zh) | 2019-04-05 |
CN109583035B CN109583035B (zh) | 2022-10-25 |
Family
ID=65921493
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811307557.4A Active CN109583035B (zh) | 2018-11-05 | 2018-11-05 | 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109583035B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110442988A (zh) * | 2019-08-12 | 2019-11-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流流向计算方法及装置 |
CN111551467A (zh) * | 2020-06-19 | 2020-08-18 | 南京工业大学 | 一种沉积颗粒物冲刷污染量的测量及计算方法 |
CN112052545A (zh) * | 2020-08-25 | 2020-12-08 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 |
CN112270115A (zh) * | 2020-11-25 | 2021-01-26 | 同济大学 | 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 |
CN114970388A (zh) * | 2022-05-12 | 2022-08-30 | 河海大学 | 一种自适应时间步长的场地地表污染扩散方法 |
CN117010232A (zh) * | 2023-06-27 | 2023-11-07 | 西安理工大学 | 基于gpu加速技术的城市非点源污染全过程高分辨率模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243186A (zh) * | 2015-09-11 | 2016-01-13 | 广西大学 | 基于元胞自动机的扩散模拟方法 |
CN106202819A (zh) * | 2016-07-26 | 2016-12-07 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于元胞自动机的水域中藻类爆发模拟方法及*** |
CN108287950A (zh) * | 2017-12-27 | 2018-07-17 | 环境保护部环境规划院 | 基于控制单元水环境质量目标管理的水质模拟方法 |
-
2018
- 2018-11-05 CN CN201811307557.4A patent/CN109583035B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243186A (zh) * | 2015-09-11 | 2016-01-13 | 广西大学 | 基于元胞自动机的扩散模拟方法 |
CN106202819A (zh) * | 2016-07-26 | 2016-12-07 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于元胞自动机的水域中藻类爆发模拟方法及*** |
CN108287950A (zh) * | 2017-12-27 | 2018-07-17 | 环境保护部环境规划院 | 基于控制单元水环境质量目标管理的水质模拟方法 |
Non-Patent Citations (3)
Title |
---|
芮小平等: "基于元胞自动机的污染物潜水面扩散仿真算法", 《科学技术与工程》 * |
赖泽辉等: "基于元胞自动机的城市地表积水模拟研究", 《水土保持通报》 * |
陈桥等: "城市地表污染物累积和降雨径流冲刷过程研究进展", 《长江流域资源与环境》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110442988A (zh) * | 2019-08-12 | 2019-11-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流流向计算方法及装置 |
CN111551467A (zh) * | 2020-06-19 | 2020-08-18 | 南京工业大学 | 一种沉积颗粒物冲刷污染量的测量及计算方法 |
CN112052545A (zh) * | 2020-08-25 | 2020-12-08 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 |
CN112052545B (zh) * | 2020-08-25 | 2021-10-08 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 |
CN112270115A (zh) * | 2020-11-25 | 2021-01-26 | 同济大学 | 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 |
CN112270115B (zh) * | 2020-11-25 | 2022-06-14 | 同济大学 | 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 |
CN114970388A (zh) * | 2022-05-12 | 2022-08-30 | 河海大学 | 一种自适应时间步长的场地地表污染扩散方法 |
CN117010232A (zh) * | 2023-06-27 | 2023-11-07 | 西安理工大学 | 基于gpu加速技术的城市非点源污染全过程高分辨率模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109583035B (zh) | 2022-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109583035A (zh) | 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 | |
Noymanee et al. | Flood forecasting with machine learning technique on hydrological modeling | |
Fettweis et al. | The mud deposits and the high turbidity in the Belgian–Dutch coastal zone, southern bight of the North Sea | |
Licciardello et al. | Runoff and soil erosion evaluation by the AnnAGNPS model in a small Mediterranean watershed | |
CN111507375B (zh) | 一种城市内涝风险快速评估方法及*** | |
Vélez et al. | Parameter extrapolation to ungauged basins with a hydrological distributed model in a regional framework | |
Wu et al. | Simulation of soil loss processes based on rainfall runoff and the time factor of governance in the Jialing River Watershed, China | |
Barker et al. | Longitudinal distributions of river flood power: the combined automated flood, elevation and stream power (CAFES) methodology | |
Fonstad | Spatial variation in the power of mountain streams in the Sangre de Cristo Mountains, New Mexico | |
CN115391712A (zh) | 一种城市洪涝风险预测方法 | |
CN113409550B (zh) | 基于径流汇流模拟的泥石流灾害预警方法及*** | |
Kabir et al. | Process-based distributed modeling approach for analysis of sediment dynamics in a river basin | |
Chen et al. | Simulation of sediment yield from watershed by physiographic soil erosion–deposition model | |
Hrissanthou et al. | Estimate of sediment inflow into Vistonis Lake, Greece | |
Mustafa et al. | Application of swat model for sediment loads from valleys transmitted to Haditha reservoir | |
Herrero et al. | Modeling the sedimentary response of a large Pyrenean basin to global change | |
Vergnes et al. | Impact of river water levels on the simulation of stream–aquifer exchanges over the Upper Rhine alluvial aquifer (France/Germany) | |
Krajewski et al. | Operation of detention pond in urban area—example of Wyścigi Pond in Warsaw | |
Pan et al. | Study on the influence of temporal and spatial resolution of Rainfall Data on Watershed Flood Simulation performance | |
Zhu | GIS based urban flood inundation modeling | |
Haddad et al. | Extreme rainfall-runoff events modeling by HEC-HMS model for Koudiet Rosfa watershed, Algeria | |
Heidarnejad et al. | Estimation of sediment volume in Karaj Dam Reservoir (Iran) by hydrometry method and a comparison with hydrography method | |
Lee et al. | Geographic information system-coupling sediment delivery distributed modeling based on observed data | |
Hrissanthou et al. | Simulation modeling of runoff, groundwater flow and sediment transport into Kastoria Lake, Greece | |
Cieśla et al. | A new concept to forecast the process of suspended sediment accumulation in the bottom sediment of small reservoirs |
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 |