CN117077445B - 一种历史遗留矿区地表水环境治理界限及阈值确定方法 - Google Patents
一种历史遗留矿区地表水环境治理界限及阈值确定方法 Download PDFInfo
- Publication number
- CN117077445B CN117077445B CN202311333394.8A CN202311333394A CN117077445B CN 117077445 B CN117077445 B CN 117077445B CN 202311333394 A CN202311333394 A CN 202311333394A CN 117077445 B CN117077445 B CN 117077445B
- Authority
- CN
- China
- Prior art keywords
- river
- surface water
- water environment
- mining area
- water
- 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
- 239000002352 surface water Substances 0.000 title claims abstract description 68
- 238000005065 mining Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 102
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 22
- 231100000719 pollutant Toxicity 0.000 claims abstract description 22
- 238000013178 mathematical model Methods 0.000 claims abstract description 17
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 229910001385 heavy metal Inorganic materials 0.000 claims description 45
- 239000002245 particle Substances 0.000 claims description 20
- 238000009826 distribution Methods 0.000 claims description 19
- 229910052751 metal Inorganic materials 0.000 claims description 16
- 239000002184 metal Substances 0.000 claims description 16
- 238000003795 desorption Methods 0.000 claims description 12
- 238000004062 sedimentation Methods 0.000 claims description 12
- 238000001179 sorption measurement Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 239000000725 suspension Substances 0.000 claims description 9
- 238000009792 diffusion process Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 239000013049 sediment Substances 0.000 claims description 7
- 239000002351 wastewater Substances 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000008021 deposition Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 4
- 230000004907 flux Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 229910052717 sulfur Inorganic materials 0.000 claims description 4
- 238000004065 wastewater treatment Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 2
- 239000010865 sewage Substances 0.000 abstract description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 23
- 238000004422 calculation algorithm Methods 0.000 description 11
- 230000006399 behavior Effects 0.000 description 8
- 229910052742 iron Inorganic materials 0.000 description 6
- 238000002386 leaching Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 241000272168 Laridae Species 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- MBMLMWLHJBBADN-UHFFFAOYSA-N Ferrous sulfide Chemical compound [Fe]=S MBMLMWLHJBBADN-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 229910052683 pyrite Inorganic materials 0.000 description 2
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 2
- 239000011028 pyrite Substances 0.000 description 2
- 238000005295 random walk Methods 0.000 description 2
- 230000002378 acidificating effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 238000003911 water pollution 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
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Water Treatment By Sorption (AREA)
Abstract
本发明涉及地表水环境质量污染控制与管理技术领域,尤其为一种历史遗留矿区地表水环境治理界限及阈值确定方法,包括如下步骤:S1:获取历史遗留矿区地表水环境资料;S2:搭建地表水环境数学模型;S3:制定典型情景方案并进行水质达标距离模拟;S4:基于典型情景方案设置污染物河流达标断面;S5:基于地表水环境数学模型计算河流污染控制断面污染负荷;S6:设置河流污染控制断面阈值。本发明构造了历史遗留矿区地表水环境治理界限及阈值的确定方法,设置河流缓冲区,将污染物排放标准与地表水环境质量标准进行有机衔接。通过确定地表水治理界限及明确河流断面管控阈值,合理利用河流纳污能力,减少过度控污的投入。
Description
技术领域
本发明涉及地表水环境质量污染控制与管理技术领域,尤其是一种历史遗留矿区地表水环境治理界限及阈值确定方法。
背景技术
我国是矿产资源消费大国,自上世纪60年代以来,我国便大力开展了矿产资源开采和选矿等生产活动。由于早期矿山开采粗放、技术落后以及环境保护意识不强等问题,矿产资源开采殆尽后,废弃矿山遗留了大量环境问题未得到有效处理。以硫铁矿区为例,废弃矿渣在空气氧化和雨水淋溶等作用下形成的淋溶水和矿硐涌水普遍为酸性废水,多数含有毒有害重金属元素,迁移到附近水体中造成河道水质部分重金属严重超标,严重破坏区域水生态环境质量,对当地居民的健康造成威胁。
历史遗留矿区是指现状废弃,今后不再进行采矿活动,由政府承担治理恢复责任。我国历史遗留矿区多,治理难度较大。目前针对历史遗留区的生态环境问题缺乏明确的环境治理目标及有效的实施路径。以硫铁矿区为例,针对废水中主要污染指标铁来看,矿区污染物排放标准可参考《铁矿采选工业污染物排放标准》(GB28661-2012),铁的直接排放浓度限值为5mg/L;附近河流地表水水质控制标准要按照《地表水环境质量标准》(GB3838-2002),铁的浓度限值为0.3mg/L。可见,矿区污染物排放标准与地表水水环境质量标准未有效衔接,可能存在即使废弃矿区污染物治理排放达标情况下,地表水仍存在一段较长距离的超标河段。从另一方面考虑,若要求排放达到地表水标准,则会存在过度治理现象。
因此,考虑到针对目前历史遗留矿区地表水污染治理与管理的空白区,亟需提出一种既能满足地表水环境管理要求,又充分考虑矿区治理的经济可行性,在河流中设置一段缓冲区,明确该段水域地表水控制标准,以明确地表水治理界限及阈值。
发明内容
本发明的目的是通过提出一种历史遗留矿区地表水环境治理界限及阈值确定方法,以解决上述背景技术中提出的缺陷。
本发明采用的技术方案如下:
提供一种历史遗留矿区地表水环境治理界限及阈值确定方法,包括如下步骤:
S1:获取历史遗留矿区地表水环境资料;
S2:搭建地表水环境数学模型;
S3:制定典型情景方案并进行水质达标距离模拟;
S4:基于典型情景方案设置污染物河流达标断面;
S5:基于地表水环境数学模型计算河流污染控制断面污染负荷;
S6:设置河流污染控制断面阈值。
作为本发明的一种优选技术方案:所述历史遗留矿区地表水环境资料具体包括气象数据、水文数据、水质数据及污染源负荷数据。
作为本发明的一种优选技术方案:所述气象数据、水文数据、水质数据及污染源负荷数据通过数据拟合得到模型输入参数。
作为本发明的一种优选技术方案:所述S2中,分别搭建水动力模块、降雨径流模块、水质模块、重金属模块的地表水环境数学模型。
作为本发明的一种优选技术方案:所述水动力模块搭建如下:
其中,Q为流量,m3/s;x为沿水流方向空间坐标,m;t为时间坐标,s;q为旁侧入流流量,入流为正,出流为负,m2/s;A为主槽过水断面面积,m2;g为重力加速度,m/s2;C为谢才系数;R为水力半径,m;α为水流速度分布系数;
所述降雨径流模块搭建如下:
其中,Si为地表储水量,m3;Qi为降雨形成的径流量,m3;A为流域面积,m2;rave流域平均降雨量,m;f为产流系数,P为常数。
所述水质模块搭建如下:
其中,c为浓度,mg/L;D为扩散系数,m2/s;A为断面面积,m2;K为线性衰减系数;c2为源(汇)项浓度,mg/L;q为旁侧入流,m2/s;x为空间坐标,m;t为时间坐标,s。
所述重金属模块主要包括金属的吸附和解吸附,吸附的金属随悬浮颗粒沉降,底泥中的金属再悬浮,以及在水和沉积物交界面间金属的扩散运输。
其中:Aadss=kwKdSHMXSS
Ddess=kwXHM
式中:XHM为水体中吸附态重金属质量浓度,g/m3;Aadss为吸附过程中溶解态重金属浓度,g/(m3·d);Ddess为解吸过程中溶解态重金属浓度,g/(m3·d);Ssev为沉积过程吸附态重金属浓度,g/(m3·d);Rresv为起悬过程吸附态重金属浓度,g/(m3·d);t为计算时段,s;kw为水体中的解吸率,d-1;Kd为重金属在颗粒态和水之间的分配系数;SHM为水体中溶解态重金属质量浓度,g/m3;XSS为水体中悬浮颗粒SS质量浓度,g/m3;Vvsm为悬浮颗粒SS的沉降速度,m/d;z为计算层厚度,m;Rresrat为悬浮颗粒SS的再悬浮率,g/(m2·d);XHMS为底泥里吸附态重金属浓度,g/m2;XSED为沉降通量,g/m2。
作为本发明的一种优选技术方案:所述S3中,在流域长时间序列河流水质分析的基础上,结合历史遗留矿区地表水环境资料进行分析,确定典型设计水文条件,进行地表水环境数学模型的训练。
作为本发明的一种优选技术方案:所述S3中,基于已筛选的典型情景方案设计水文条件,开展河流水环境质量模拟,分析在典型情景方案下,河流污染物浓度的分布情况及确定污染物指标的达标距离,获得典型情景方案下污染物河流达标断面。
作为本发明的一种优选技术方案:所述S6中,将河流污染控制断面最大允许污染负荷与矿区废水按照排放标准达标排放的污染负荷进行对比,设置河流污染控制断面阈值。
作为本发明的一种优选技术方案:所述S6中,设置河流污染控制断面至河流达标控制断面为河段缓冲区,河段污染控制断面为矿区污染负荷输入的概化点,基于地表水环境质量标准获取河流污染控制断面阈值要求,设置矿区废水治理目标。
本发明提供的历史遗留矿区地表水环境治理界限及阈值确定方法,与现有技术相比,其有益效果有:
本发明构造了历史遗留矿区地表水环境治理界限及阈值的确定方法。其填补了历史遗留矿区所在流域地表水环境管理目标不明确的现实问题;设置河流缓冲区,有效地将污染物排放标准与地表水环境质量标准进行有机衔接。通过合理的确定地表水治理界限及明确河流断面管控阈值,在一定程度上合理利用河流纳污能力,减少过度控污的投入。总体而言,该方法为历史遗留矿区所在流域水环境管理提供了有力的技术支撑。
附图说明
图1为本发明优选实施例的方法流程图;
图2为本发明优选实施例中流程示意图;
图3为本发明优选实施例中W县地表水流向示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本实施例中的实施例及实施例中的特征可以相互组合,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1:
参照图1,本发明优选实施例提供了一种历史遗留矿区地表水环境治理界限及阈值确定方法,包括如下步骤:
S1:获取历史遗留矿区地表水环境资料;
S2:搭建地表水环境数学模型;
S3:制定典型情景方案并进行水质达标距离模拟;
S4:基于典型情景方案设置污染物河流达标断面;
S5:基于地表水环境数学模型计算河流污染控制断面污染负荷;
S6:设置河流污染控制断面阈值。
所述历史遗留矿区地表水环境资料具体包括气象数据、水文数据、水质数据及污染源负荷数据。
所述S2中,分别搭建水动力模块、降雨径流模块、水质模块、重金属模块的地表水环境数学模型。
所述气象数据、水文数据、水质数据及污染源负荷数据通过数据拟合得到模型输入参数。可以利用流量-水位关系曲线(流量-水位测验)估计水动力模型中的断面流量和水位的参数。降雨径流模型中的参数可以通过降雨和径流数据拟合得到。水质模型中的参数可以通过拟合实测水质数据得到。
所述水动力模块搭建如下:
其中,Q为流量,m3/s;x为沿水流方向空间坐标,m;t为时间坐标,s;q为旁侧入流流量,入流为正,出流为负,m2/s;A为主槽过水断面面积,m2;g为重力加速度,m/s2;C为谢才系数;R为水力半径,m;α为水流速度分布系数;
所述降雨径流模块搭建如下:
其中,Si为地表储水量,m3;Qi为降雨形成的径流量,m3;A为流域面积,m2;rave流域平均降雨量,m;f为产流系数,P为常数。
所述水质模块搭建如下:
其中,c为浓度,mg/L;D为扩散系数,m2/s;A为断面面积,m2;K为线性衰减系数;c2为源(汇)项浓度,mg/L;q为旁侧入流,m2/s;x为空间坐标,m;t为时间坐标,s。
所述重金属模块主要包括金属的吸附和解吸附,吸附的金属随悬浮颗粒沉降,底泥中的金属再悬浮,以及在水和沉积物交界面间金属的扩散运输。
其中:Aadss=kwKdSHMXSS
Ddess=kwXHM
式中:XHM为水体中吸附态重金属质量浓度,g/m3;Aadss为吸附过程中溶解态重金属浓度,g/(m3·d);Ddess为解吸过程中溶解态重金属浓度,g/(m3·d);Ssev为沉积过程吸附态重金属浓度,g/(m3·d);Rresv为起悬过程吸附态重金属浓度,g/(m3·d);t为计算时段,s;kw为水体中的解吸率,d-1;Kd为重金属在颗粒态和水之间的分配系数;SHM为水体中溶解态重金属质量浓度,g/m3;XSS为水体中悬浮颗粒SS质量浓度,g/m3;Vvsm为悬浮颗粒SS的沉降速度,m/d;z为计算层厚度,m;Rresrat为悬浮颗粒SS的再悬浮率,g/(m2·d);XHMS为底泥里吸附态重金属浓度,g/m2;XSED为沉降通量,g/m2。
所述S3中,在流域长时间序列河流水质分析的基础上,结合历史遗留矿区地表水环境资料进行分析,确定典型设计水文条件,进行地表水环境数学模型的训练。
所述S3中,基于已筛选的典型情景方案设计水文条件,开展河流水环境质量模拟,分析在典型情景方案下,河流污染物浓度的分布情况及确定污染物指标的达标距离,获得典型情景方案下污染物河流达标断面。
所述S6中,将河流污染控制断面最大允许污染负荷与矿区废水按照排放标准达标排放的污染负荷进行对比,设置河流污染控制断面阈值。
所述S6中,设置河流污染控制断面至河流达标控制断面为河段缓冲区,河段污染控制断面为矿区污染负荷输入的概化点,基于地表水环境质量标准获取河流污染控制断面阈值要求,设置矿区废水治理目标。
本实施例中,参照图2,以W县的为例进行模型的搭建和训练:
W县境内废弃矿山导致的水环境问题由来已久,上世纪60年代开始,硫铁矿长时间无序开采,现存大量的历史遗留矿山,对周边水体环境造成污染。但由于我国针对历史遗留矿区重金属的水环境污染治理目标一直不明确,治理界限不清晰,造成很长一段时间历史遗留矿区地表水环境治理受到制约,矿区周边河流由于长期受到矿堆淋溶水和矿硐涌水的影响,受到了严重污染。参照图3,本例以W县境内D河流域内历史遗留矿山为例,D河属于矿山废水直接排入的支流,D河汇入H河,H河汇入B河。从现状情况来看,B河已一定程度受到重金属污染。因此,根据本案技术方法,通过构建河流缓冲区,确定地表水环境治理界限及阈值,为历史遗留矿山治理明确目标。
收集W县所在区域1990~2019年近30年逐日降雨数据,2010~2019年近10年每年60min、180min、360min三个降雨历时逐小时过程数据,收集2015~2019年逐月流域考核断面水质数据,以及结合后续提出的不同典型情景方案下矿区污染负荷输出量。本实施例主要针对硫铁矿主要污染指标及流域地表水主要超标指标Fe举例说明。
分别搭建水动力模块、降雨径流模块、水质模块、重金属模块的地表水环境数学模型。
水动力模块搭建如下:
其中,Q为流量,m3/s;x为沿水流方向空间坐标,m;t为时间坐标,s;q为旁侧入流流量,入流为正,出流为负,m2/s;A为主槽过水断面面积,m2;g为重力加速度,m/s2;C为谢才系数;R为水力半径,m;α为水流速度分布系数;
降雨径流模块搭建如下:
其中,Si为地表储水量,m3;Qi为降雨形成的径流量,m3;A为流域面积,m2;rave流域平均降雨量,m;f为产流系数,P为常数。
水质模块搭建如下:
其中,c为浓度,mg/L;D为扩散系数,m2/s;A为断面面积,m2;K为线性衰减系数;c2为源(汇)项浓度,mg/L;q为旁侧入流,m2/s;x为空间坐标,m;t为时间坐标,s。
所述重金属模块主要包括金属的吸附和解吸附,吸附的金属随悬浮颗粒沉降,底泥中的金属再悬浮,以及在水和沉积物交界面间金属的扩散运输。
其中:Aadss=kwKdSHMXSS
Ddess=kwXHM
式中:XHM为水体中吸附态重金属质量浓度,g/m3;Aadss为吸附过程中溶解态重金属浓度,g/(m3·d);Ddess为解吸过程中溶解态重金属浓度,g/(m3·d);Ssev为沉积过程吸附态重金属浓度,g/(m3·d);Rresv为起悬过程吸附态重金属浓度,g/(m3·d);t为计算时段,s;kw为水体中的解吸率,d-1;Kd为重金属在颗粒态和水之间的分配系数;SHM为水体中溶解态重金属质量浓度,g/m3;XSS为水体中悬浮颗粒SS质量浓度,g/m3;Vvsm为悬浮颗粒SS的沉降速度,m/d;z为计算层厚度,m;Rresrat为悬浮颗粒SS的再悬浮率,g/(m2·d);XHMS为底泥里吸附态重金属浓度,g/m2;XSED为沉降通量,g/m2。
经率定验证,水动力模型流量结果率定较好,相对误差均在10%以内,水质模拟(Fe)相对误差基本在35%以内,模型能够较真实反映现实情况,可支撑后续情景方案模拟及控制断面阈值确定。
参照表1,利用1990~2019年逐日降雨数据,分析得到流域平均年降雨量约817.59mm,其中发生暴雨(降雨量≥50mm)的年均天数为2天,暴雨的年降雨均值为121.4mm,占年均总降雨量的14.85%;发生大雨(49.9mm≥降雨量≥25mm)的年均天数为7天,大雨的年降雨均值为240.6mm,占年均总降雨量的29.43%;发生中雨(24.9mm≥降雨量≥10mm)的年均天数为15天,中雨的年降雨均值为244.4mm,占年均总降雨量的29.89%;发生小雨(降雨量<10mm)的年均天数为126天,小雨的年降雨均值为211.2mm,占年均总降雨量的25.83%。
表1研究区域1990~2019年气象站降雨量统计情况
经分析,流域暴雨具有明显的季节性,暴雨出现的时序有一定的规律。暴雨最早发生在4月,最迟发生在11月,但量级及强度较大的暴雨一般在6~9月。流域大雨在3~11月份均有发生,主要集中在5~10月;中雨全年均有发生,主要集中在3~11月;小雨全年均有发生,各月份发生次数差异相对较小。
利用2010~2019年降雨过程,开展典型降雨过程曲线分析。将雨峰时段定位在出现几率最大的位置,雨峰时段在总雨量中所占的比例定义为各场降雨雨峰在总雨量中所占的比例平均值,其余各时段的具体位置和所占比例采用同样方法原理定义。
根据污染负荷调查情况及流域水环境质量评估,流域在平水期呈现出Fe的污染更为严重的现象,根据降雨在平水期发生的概率来看,各类雨都有可能在平水期发生,但考虑到降雨主要影响淋溶水对河流的污染,筛选暴雨、中雨开展典型情景方案研究,即确定的典型情景方案为平水期+暴雨、平水期+中雨。平水期选取流域平水期多年平均流量。
根据流域典型情景方案平水期+暴雨、平水期+中雨设计水文条件,基于已构建的水环境数学模型,计算得到D河汇入口、H河汇入口及B河汇入口处Fe的浓度值。参照表2,经计算,典型情景方案下以上断面Fe浓度均远远超过地表水环境质量目标要求,河段均受到一定程度污染。
表2典型情景方案下各断面Fe浓度
断面 | 平水期+中雨 | 平水期+暴雨 |
D河汇入口 | 8.50 | 37.33 |
H河汇入口 | 8.10 | 22.33 |
B河汇入口 | 1.32 | 3.81 |
达标控制断面即水质目标要求达到《地表水环境质量标准》(GB3838-2002)相关要求,即Fe的浓度限值为0.3mg/L。根据河流达标控制断面设置原则,D河上无水质考核断面,且矿区污染负荷汇入点距离D河汇入H河长度小于10km时,按照筛选原则D河汇入口作为达标控制断面。同时,典型情景方案模拟的达标断面在D河汇入口下游,因此,最终确定D河汇入口为达标控制断面。
基于已构建的水环境数学模型,分别计算在平水期+暴雨、平水期+中雨两种典型设计水文条件下,反算D河汇入口断面Fe浓度达到0.3mg/L时,河流污染控制断面最大允许污染负荷,即历史遗留矿山废水排入D河概化断面。经计算,D河河流污染控制断面Fe的最大允许负荷量为1625mg/s。
D河流污染控制断面Fe的最大允许负荷量为1625mg/s,按照《铁矿采选工业污染物排放标准》(GB28661-2012)计算得到D河片区历史遗留矿区达标排放的污染负荷为54mg/s,可见D河最大允许负荷量较历史遗留矿区达标排放的污染负荷量大,D河污染控制断面多年平均流量为0.05m3/s,计算得到污染控制断面浓度管控限值为32.05mg/L。由于计算的浓度管控限值大于《铁矿采选工业污染物排放标准》(GB28661-2012)中Fe的排放标准5mg/L,因此建议污染控制断面管控浓度按5mg/L进行设定。
实施例2
作为本发明的另外一种实施,与实施例1的不同之处在于:
本实施例中所述地表水环境的水动力模型、降雨径流模型、水质模型、重金属模型中,基于改进海鸥优化算法对各模型的参数进行寻优调整获得最优模型。
所述改进海鸥优化算法具体如下:
yI=δ×YI(τ)
其中,yI为第I只海鸥个体完成避免碰撞后的位置,δ为表示海鸥个体在搜索空间中的运动行为的附加变量,YI(τ)为第τ次迭代第I只海鸥个体的位置,R1为从2线性递减至0的随机数;
bI=B×(Yb(τ)-YI(τ))
B=2×δ×δ×R2
其中,Yb(τ)是第τ次迭代的最优海鸥个体,B为用于平衡算法的全局搜索和局部开发,bI为最优海鸥个体所在的方向,R2为[0,1]之间的随机数;
海鸥根据最优海鸥个体的位置信息更新位置:
ZI=|yI+bI|
其中,ZI为海鸥个体与最优海鸥个体之间的距离;
海鸥个体在攻击行为中以螺旋的飞行轨迹攻击捕食猎物,在x、y、z平面上,攻击行为为:
x′=r×cos q
y′=r×sin q
z′=r×q
r=ε×a×eqb
YI(τ)=(ZI×x′×y′×z′)+Yb(τ)+ε×θ
其中,r为海鸥个体螺旋飞行轨迹每一圈的半径,q为[0,2π]之间的一个随机数,a、b为用于定义螺旋飞行轨迹形状的常数,ε为动态收敛因子,θ为服从Levy分布的随机数向量;
其中,Levy(θ)表示参数θ服从Levy分布,r1、r2分别为服从标准差为σ1、σ2的正态分布随机数,其中,σ2=1,Γ(1+θ)为Gamma函数。
在具体实现过程中,在算法迭代初期,迁徙行为可以帮助海鸥种群更好地遍历整个搜索空间。在迁徙行为中,相邻海鸥个体之间为避免发生碰撞,通过下式更新位置:
yI=δ×YI(5)
其中,yI为第I只海鸥个体完成避免碰撞后的位置,δ为表示海鸥个体在搜索空间中的运动行为的附加变量,YI(5)为第5次迭代第I只海鸥个体的位置,R1为从2线性递减至0的随机数;
bI=B×(Yb(5)-YI(5))
B=2×δ×δ×R2
其中,Yb(5)是第5次迭代的最优海鸥个体,B为用于平衡算法的全局搜索和局部开发,bI为最优海鸥个体所在的方向,R2为[0,1]之间的随机数;
海鸥根据最优海鸥个体的位置信息更新位置:
ZI=|yI+bI|
其中,ZI为海鸥个体与最优海鸥个体之间的距离;
海鸥个体在飞行过程中会不断地改变速度和攻击角度。在攻击行为中,海鸥个体以螺旋的飞行轨迹攻击捕食猎物,在x、y、z平面上,攻击行为为:
x′=r×cos q
y′=r×sin q
z′=r×q
r=ε×a×eqb
YI(5)=(ZI×x′×y′×z′)+Yb(5)+ε×θ
其中,r为海鸥个体螺旋飞行轨迹每一圈的半径,q为[0,2π]之间的一个随机数,a、b为用于定义螺旋飞行轨迹形状的常数,ε为动态收敛因子,θ为服从Levy分布的随机数向量;
其中,Levy(θ)表示参数θ服从Levy分布,r1、r2分别为服从标准差为σ1、σ2的正态分布随机数,其中,σ2=1,Γ(1+θ)为Gamma函数,Γ(1+θ)=θ!;
在算法迭代初期,较大的攻击范围可以帮助海鸥个体更好地探索整个搜索空间,以提高种群的多样性。随着迭代次数的增加,在动态收敛因子ε的作用下攻击范围逐渐减小,海鸥个体围绕最优解进行精细地搜索,算法具备更强的局部开发能力。考虑到算法在迭代末期需要进行精细的搜索,为进一步提高算法的全局寻优能力,将Levy分布与动态收敛因子ε结合并应用到海鸥位置的更新过程。Levy分布中随机游走的特性可以提高算法跳出局部最优的能力,而动态收敛因子ε则帮助海鸥个体在迭代末期围绕最优解进行小范围的随机游走,以进一步强化算法的局部开发能力。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (4)
1.一种历史遗留矿区地表水环境治理界限及阈值确定方法,其特征在于,包括如下步骤:
S1:获取历史遗留矿区地表水环境资料;
S2:搭建地表水环境数学模型;
S3:制定典型情景方案并进行水质达标距离模拟;
S4:基于典型情景方案设置污染物河流达标断面;
S5:基于地表水环境数学模型计算河流污染控制断面污染负荷;
S6:设置河流污染控制断面阈值;
所述S2中,分别搭建水动力模块、降雨径流模块、水质模块、重金属模块的地表水环境数学模型;
所述降雨径流模块搭建如下:
其中,Si为地表储水量,m3;Qi为降雨形成的径流量,m3;Ab为流域面积,m2;rave流域平均降雨量,m;f为产流系数,K为线性衰减系数,P为常数;
所述S6中,将河流污染控制断面最大允许污染负荷与矿区废水按照排放标准达标排放的污染负荷进行对比,设置河流污染控制断面阈值;
所述S6中,设置河流污染控制断面至河流达标控制断面为河段缓冲区,河段污染控制断面为矿区污染负荷输入的概化点,基于地表水环境质量标准获取河流污染控制断面阈值要求,设置矿区废水治理目标;
所述S3中,在流域长时间序列河流水质分析的基础上,结合历史遗留矿区地表水环境资料进行分析,确定典型设计水文条件,进行地表水环境数学模型的训练;
所述S3中,基于制定的典型情景方案设计水文条件,开展河流水环境质量模拟,分析在典型情景方案下,河流污染物浓度的分布情况及确定污染物指标的达标距离,获得典型情景方案下污染物河流达标断面。
2.根据权利要求1所述的历史遗留矿区地表水环境治理界限及阈值确定方法,其特征在于:所述历史遗留矿区地表水环境资料具体包括气象数据、水文数据、水质数据及污染源负荷数据。
3.根据权利要求2所述的历史遗留矿区地表水环境治理界限及阈值确定方法,其特征在于:所述气象数据、水文数据、水质数据及污染源负荷数据通过数据拟合得到模型输入参数。
4.根据权利要求3所述的历史遗留矿区地表水环境治理界限及阈值确定方法,其特征在于,所述水动力模块搭建如下:
其中,Q为流量,m3/s;x为沿水流方向空间坐标,m;t为时间坐标,s;q为旁侧入流流量,入流为正,出流为负,m2/s;A为主槽过水断面面积,m2;g为重力加速度,m/s2;C为谢才系数;R为水力半径,m;α为水流速度分布系数;
所述水质模块搭建如下:
其中,c为浓度,mg/L;D为扩散系数,m2/s;A为主槽过水断面面积,m2;K为线性衰减系数;c2为源项浓度,mg/L;q为旁侧入流,m2/s;x为空间坐标,m;t为时间坐标,s;
所述重金属模块主要包括金属的吸附和解吸附,吸附的金属随悬浮颗粒沉降,底泥中的金属再悬浮,以及在水和沉积物交界面间金属的扩散运输;
其中:
Aadss=kwKdSHMXSS
Ddess=kwXHM
式中:XHM为水体中吸附态重金属质量浓度,g/m3;Aadss为吸附过程中溶解态重金属浓度,g/(m3·d);Ddess为解吸过程中溶解态重金属浓度,g/(m3·d);Ssev为沉积过程吸附态重金属浓度,g/(m3·d);Rresv为起悬过程吸附态重金属浓度,g/(m3·d);t为计算时段,s;kw为水体中的解吸率,d-1;Kd为重金属在颗粒态和水之间的分配系数;SHM为水体中溶解态重金属质量浓度,g/m3;XSS为水体中悬浮颗粒SS质量浓度,g/m3;Vvsm为悬浮颗粒SS的沉降速度,m/d;z为计算层厚度,m;Rresrat为悬浮颗粒SS的再悬浮率,g/(m2·d);XHMS为底泥里吸附态重金属浓度,g/m2;XSED为沉降通量,g/m2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311333394.8A CN117077445B (zh) | 2023-10-16 | 2023-10-16 | 一种历史遗留矿区地表水环境治理界限及阈值确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311333394.8A CN117077445B (zh) | 2023-10-16 | 2023-10-16 | 一种历史遗留矿区地表水环境治理界限及阈值确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117077445A CN117077445A (zh) | 2023-11-17 |
CN117077445B true CN117077445B (zh) | 2024-05-07 |
Family
ID=88706438
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311333394.8A Active CN117077445B (zh) | 2023-10-16 | 2023-10-16 | 一种历史遗留矿区地表水环境治理界限及阈值确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117077445B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102201034A (zh) * | 2011-04-26 | 2011-09-28 | 北京师范大学 | 一种临河村落非点源污染负荷的估算方法 |
CN107563139A (zh) * | 2017-09-06 | 2018-01-09 | 郑州大学 | 一种点源突发性水污染源事故溯源贡献程度计算方法 |
CN111581805A (zh) * | 2020-04-30 | 2020-08-25 | 沈阳理工大学 | 一种控制单元污染物输入与水质的响应分析方法 |
CN112287613A (zh) * | 2020-10-20 | 2021-01-29 | 河海大学 | 针对流域水环境控制断面的污染物削减方法 |
CN114997591A (zh) * | 2022-05-07 | 2022-09-02 | 河海大学 | 一种基于水环境数学模型的入河污染物削减及河道断面水质预测方法及装置 |
WO2023087995A1 (zh) * | 2021-11-16 | 2023-05-25 | 河海大学 | 一种分汊河道分流比与分污比及其演变规律的研究方法 |
CN116245042A (zh) * | 2022-12-22 | 2023-06-09 | 河海大学 | 基于水动力-水环境-水生生物的生态需水计算方法 |
CN116303815A (zh) * | 2023-01-04 | 2023-06-23 | 广东工业大学 | 一种河流断面的多级指纹图谱构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113689099B (zh) * | 2021-08-16 | 2022-05-06 | 河海大学 | 基于水量构成的环状河网水系水质考核断面污染负荷贡献率确定方法 |
-
2023
- 2023-10-16 CN CN202311333394.8A patent/CN117077445B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102201034A (zh) * | 2011-04-26 | 2011-09-28 | 北京师范大学 | 一种临河村落非点源污染负荷的估算方法 |
CN107563139A (zh) * | 2017-09-06 | 2018-01-09 | 郑州大学 | 一种点源突发性水污染源事故溯源贡献程度计算方法 |
CN111581805A (zh) * | 2020-04-30 | 2020-08-25 | 沈阳理工大学 | 一种控制单元污染物输入与水质的响应分析方法 |
CN112287613A (zh) * | 2020-10-20 | 2021-01-29 | 河海大学 | 针对流域水环境控制断面的污染物削减方法 |
WO2023087995A1 (zh) * | 2021-11-16 | 2023-05-25 | 河海大学 | 一种分汊河道分流比与分污比及其演变规律的研究方法 |
CN114997591A (zh) * | 2022-05-07 | 2022-09-02 | 河海大学 | 一种基于水环境数学模型的入河污染物削减及河道断面水质预测方法及装置 |
CN116245042A (zh) * | 2022-12-22 | 2023-06-09 | 河海大学 | 基于水动力-水环境-水生生物的生态需水计算方法 |
CN116303815A (zh) * | 2023-01-04 | 2023-06-23 | 广东工业大学 | 一种河流断面的多级指纹图谱构建方法 |
Non-Patent Citations (8)
Title |
---|
基于MIKE 11模型的引江济淮工程涡河段动态水环境容量研究;熊鸿斌;张斯思;匡武;吴蕾;朱慧娈;;自然资源学报;20170815(08);154-164 * |
基于一维水动力水质模型的纳污能力分段核定研究――以江门市江海区礼乐河为例;张家鸣;刘继艳;;人民珠江(07);85-88 * |
基于京杭运河五牧断面水质达标的水环境容量计算研究;韩梓流;逄勇;;水资源与水工程学报(04);81-87 * |
基于控制单元的襄河流域容量总量控制及水质达标对策;包琳琳;吴楠;;安徽农业大学学报(06);1084-1092 * |
基于龙王庙断面水质达标的城南河流域水环境容量;瞿一清;逄勇;;水资源保护(05);76-80 * |
太湖梅梁湾水环境监控预警体系研究;黄卫 等;环境监控与预警;6-9 * |
淀山湖流域平原河网水环境容量及控制断面水质达标方案研究;王子轩;逄勇;罗缙;王一舒;;水资源与水工程学报(06);61-65 * |
金属矿区重金属迁移对水体影响的数值模拟;轩晓博 等;水资源保护;30-35 * |
Also Published As
Publication number | Publication date |
---|---|
CN117077445A (zh) | 2023-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gholami et al. | A method of groundwater quality assessment based on fuzzy network-CANFIS and geographic information system (GIS) | |
Grande | Impact of AMD processes on the public water supply: hydrochemical variations and application of a classification model to a river in the Iberian Pyritic Belt SW Spain | |
Braghina et al. | Territorial management within the systems affected by mining. Case study the South-Western Development Region in Romania | |
CN112990739A (zh) | 基于多介质污染特征识别的磷矿区水环境风险评估方法 | |
Sarma et al. | Urban hydrology, watershed management and socio-economic aspects | |
CN117077445B (zh) | 一种历史遗留矿区地表水环境治理界限及阈值确定方法 | |
Laktionova et al. | Financing the greening of enterprises in industrial regions of Ukraine in the context of sustainable development | |
Akhter et al. | Climate modeling of Jhelum River basin-a comparative study | |
Uncles et al. | Turbidity in the Thames Estuary: How turbid do we expect it to be? | |
Wang et al. | Monitored and modeled correlations of sediment and nutrients with Chesapeake Bay water clarity | |
Omar et al. | Effectiveness of dredging and drains’ treatment on water quality of Rosetta branch | |
Fang et al. | Zero increase in peak discharge for sustainable development | |
Guan et al. | Probing the national development from heavy metals contamination in river sediments | |
CN114331787A (zh) | 以提升下游河流水质为目标的水库生态放流量的核算方法 | |
KR101859109B1 (ko) | 강우예측을 이용한 저영향개발시설의 제어방법 | |
CN111598757B (zh) | 一种平原河网区城市水系流动确定方法 | |
Turatsinze et al. | Quantification of sediments and delineation of their sources in Nyabugogo Stream, RWANDA: Relationship with economic development activities | |
Lim et al. | Prediction of pollution loads in agricultural reservoirs using LSTM algorithm: case study of reservoirs in Nonsan City | |
CN113158591A (zh) | 一种确定流域土地开发利用承载力的方法 | |
Braghină et al. | Territorial management within the systems affected by mining. Case study the South-Western Development Region in Romania. | |
CN113743760B (zh) | 河湖水系连通性与水体水质适配性分析方法 | |
Moon et al. | Determining the effective width of riparian buffers in Korean watersheds using the SWAT model. | |
Jha et al. | Analytical water quality model for biochemical oxygen demand simulation in River Gomti of Ganga Basin, India | |
Schastlivtsev et al. | Assessment of anthropogenic impact on river basins in southern Kuzbass mining areas | |
Babayan et al. | ASSESSING THE QUALITY OF WATER FROM SURFACE SOURCES OF DRINKING WATER SUPPLY TO TOWNS OF SOUTH ARMENIA. |
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 |