CN115329607A - 一种用于地下水污染的评估***及方法 - Google Patents
一种用于地下水污染的评估***及方法 Download PDFInfo
- Publication number
- CN115329607A CN115329607A CN202211256594.3A CN202211256594A CN115329607A CN 115329607 A CN115329607 A CN 115329607A CN 202211256594 A CN202211256594 A CN 202211256594A CN 115329607 A CN115329607 A CN 115329607A
- Authority
- CN
- China
- Prior art keywords
- groundwater
- boundary
- numerical
- concentration
- calculation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000003911 water pollution Methods 0.000 title claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 80
- 238000005070 sampling Methods 0.000 claims abstract description 49
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 31
- 238000011156 evaluation Methods 0.000 claims abstract description 21
- 238000009792 diffusion process Methods 0.000 claims abstract description 18
- 239000003245 coal Substances 0.000 claims abstract description 14
- 238000005516 engineering process Methods 0.000 claims abstract description 10
- 238000003895 groundwater pollution Methods 0.000 claims abstract description 10
- 239000003673 groundwater Substances 0.000 claims description 70
- 238000004364 calculation method Methods 0.000 claims description 68
- 238000005065 mining Methods 0.000 claims description 23
- 238000012544 monitoring process Methods 0.000 claims description 16
- 238000011160 research Methods 0.000 claims description 13
- 230000009189 diving Effects 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 239000000470 constituent Substances 0.000 claims description 8
- 238000011109 contamination Methods 0.000 claims description 7
- 239000003814 drug Substances 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 239000004480 active ingredient Substances 0.000 claims description 3
- 239000008346 aqueous phase Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000008595 infiltration Effects 0.000 claims description 3
- 238000001764 infiltration Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000004080 punching Methods 0.000 claims description 3
- 239000002349 well water Substances 0.000 claims description 3
- 235000020681 well water Nutrition 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 2
- 239000012528 membrane Substances 0.000 claims description 2
- 238000005325 percolation Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000001223 reverse osmosis Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
- G01N33/182—Specific anions in water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Medicinal Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Sampling And Sample Adjustment (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及预测评估的数据处理领域,具体公开了一种用于地下水污染的评估***及方法,包括数据采集模块、数值反演追踪模块和污染评估模块。本发明通过数据采集模块的设置能够采集闭坑煤矿区域的地质条件信息和水文地质条件信息,使用SO4 2‑的扩散代表矿坑积水反渗造成地下水污染的扩散,利用质量守恒和SO4 2‑扩散规律建立数值模型,并对研究三维区域设置网格,结合已知的边界条件和初始条件求解污染源的确定位置坐标,在标定的目标处取样,分析数值模拟结果的可信度,解决了现有的地下水污染评估***中没有结合数值模拟技术反演追踪煤矿串层污染井位置的不足,提高了煤矿串层污染井治理的精确度。
Description
技术领域
本发明涉及预测评估的数据处理领域,具体涉及一种用于地下水污染的评估***及方法。
背景技术
由于地下水处于地表以下,且存在于岩石的空隙中,因此地下水被污染后难以被发现,地下水污染后无法将地表全部挖开寻找污染源,只能在地表开挖若干处取样井,通过取样井的地下水数据对地下水的污染及污染源进行评估,通过数值模拟的方法结合地下水采样监测反演识别地下水污染源就变得十分重要,尤其是在煤矿中通过数值模拟手段对地下水污染的评估更是一种亟待发展的技术。
在文献《煤矿区场地地下水污染防控技术研究进展及发展方向》中提及煤矿在经过长期大规模开采会闭坑停产,闭坑停产后矿坑水的停排引起地下水污染,但是煤矿在开挖的过程中会根据煤层所在的位置做复杂的开采,矿坑的走向和分布错综复杂,所以污染地下水的矿坑积水位置较难寻找,寻找引起地下水污染的串层污染井,进行有效封堵、投放物化材料或者药剂是防治矿区地下水污染加重的重要手段,但是现有的地下水污染评估***中还没有使用数值模拟的手段利用地下水采样监测数据反演识别污染地下水的串层污染井位置,针对以上不足,本发明提出一种用于地下水污染的评估***及方法,利用数值模拟手段结合地下水采样监测数据反演识别引起地下水污染的串层污染井位置,确定引起地下水污染的串层污染井位置后,有效封堵或者引起地下水污染的串层污染井,能够防止地下水污染加重。
发明内容
本发明所要解决的技术问题在于:如何解决现有的地下水污染评估***中没有利用地下水采样监测数据结合数值模拟方手段,反演识别污染地下水的串层污染井位置的不足,本发明提供了一种用于地下水污染的评估***及方法。
本发明是通过以下技术方案解决上述技术问题的,本发明包括数据采集模块、数值反演追踪模块和污染评估模块,所述数据采集模块用于将研究矿区的地质和水文地质条件以及采样井地下水监测数据传输给所述数值反演追踪模块,所述数值反演追踪模块根据采集的数据信息设定水文地质条件、计算流场边界以及边界已知条件,建立多相流数值模拟模型,所述数值反演追踪模块连接有污染评估模块;
用于地下水污染的评估***通过数值模拟技术反演追踪闭坑煤矿区地下水污染源的方法如下:
步骤一:数据采集模块采集采样井地下水水质和水位信息、矿区的地质条件和水文地质条件以及串层污染井的范围,并将数据传输给数值反演追踪模块;
步骤二:数值反演追踪模块对研究的区域进行网格划分,对串层污染井范围内局部加密,数值反演追踪模块根据数据采集模块提供的数据,基于SO4 2-的扩散建立多相流数值模拟模型,根据数据采集模块提供的数据设定数值模拟初始条件、边界条件,圈定串层污染井的范围,求解污染源的三维坐标;
步骤三:根据数值模拟技术得到三维坐标,在坐标位置打孔取样,验证污染源位置的地下水样本中SO4 2-的浓度是否与数值模拟计算得到的计算浓度接近;
步骤四:在不同的串层污染井采用同样的方法进行数值模拟,收集数值模拟结果和采样结果;
步骤五:对数值模拟的计算结果可信度进行分析,利用数值模拟计算得出SO4 2-浓度减去实际采样地下水中SO4 2-浓度的绝对值,除以实际采样地下水中SO4 2-的浓度,得到的百分比作为可信度评估的依据,所得到的百分比数值大于80%的记为可信样本组,并对可信样本组对应的串层污染井进行有效封堵或者投放药剂进行治理。
作为本发明进一步的方案,所述数据采集模块采集的矿区地质条件、水文地质条件以及地下水采样井实时监测数据包括:
实际研究场地的三维区域边界:用于划定数值模拟计算流场的三维计算区域,根据实际的场地大小设定计算的空间步长L,并根据空间步长将实际研究场地的大小划分为的三维区域,式中aL、bL和cL记为实际研究场地的长度、宽度和深度,并将计算区域内的坐标集合记为;
待评估矿区地层性质和分类:用于设定数值计算潜水含水层的类型,并根据大部分的潜水含水层类型设定潜水含水层的性质和参数;
研究矿区四周的地下水流量以及水压信息:设定计算区域的上、下边界和左、右边界类型及参数,设定模拟区域的上边界和下边界均为不透水的零通量边界,分别将上边界和下边界记为和,将左边界和右边界设定为已知压力的边界,并分别将左边界和右边界记为和;
设定监测数据的取样时间间隔:作为数值模拟计算的监测时刻划分依据;
采样井信息:用于在数值模拟计算中对采样井编号并实施数值计算;
采样井水质信息中的SO4 2-浓度:通过采样的方式从待求的污染源位置SO4 2-浓度范围中抽取300组数据代入多相流数值模拟模型,并对其进行反演,取100组接近取样井SO4 2-浓度的数据进行核对,取准确度达到提前设定要求的数据进行统计计算。
作为本发明进一步的方案,所述数值追踪反演模块的数值建模信息包括:
初始条件:包括模拟区的初始压力分布、初始SO4 2-以及含水相以及初始SO4 2-的浓度,初始条件为:
式中:
边界条件:根据水文地质条件设定已知的流量边界和压力边界,根据矿区的暴雨强度及周边河流信息设定源汇项以及SO4 2-浓度边界,边界条件为:
式中:记为地下水中SO4 2-的浓度;记为地下水的实际平均流速,记为SO4 2-的饱和度,记为地下水中SO4 2-的弥散系数张量,n记为边界的法线方向,记为上部边界上的已知压力分布函数,记为左边界的已知压力函数,记为右边界的已知压力函数,记为上部边界SO4 2-的已知浓度函数,记为左边界的SO4 2-已知浓度函数,记为右边界的SO4 2-已知浓度函数;
时间步长:根据采样井的采样时间间隔设定数值模拟的计算时间步长T,并选取需要模拟计算的计算时刻t=10T,11T,12T,13T,14T;
多相流数值模拟计算偏微分方程:根据质量守恒定律建立多相流偏微分方程,多相流偏微分方程为:
式中:s记为组成物质种类,s=1记为组成物质为水,s=2记为组成物质为SO4 2-,记为研究模拟区域地层的孔隙度,记为地下水中SO4 2-的总浓度,记为组成物质的密度,记为地下水中SO4 2-的浓度,记为地下水中SO4 2-的饱和度,记为地下水中SO4 2-的弥散系数张量;
式中:D记为SO4 2-在地下水中的分子扩散系数,记为地层的扭曲系数,记为克罗内克函数,和分别记为SO4 2-在地下水中的横向弥散度和纵向弥散度,记为SO4 2-的渗流速度,和分别表示SO4 2-在i方向和在j方向上的渗流分速度;
作为本发明进一步的方案,所述数值反演追踪模块连接有显示屏,所述显示屏用于显示数值模拟计算结果,并将计算结果显示在矿区三维地图上,标定可信度最高的串层污染井三维坐标,同时显示此处SO4 2-的浓度。
作为本发明进一步的方案,所述污染评估模块对所述数值反演追踪模块的计算结果可信度进行评估,所述可信度计算公式为:
本发明相比现有技术具有以下优点:该用于地下水污染的评估***及方法,通过数据采集模块的设置能够采集闭坑煤矿区域的地质条件信息和水文地质条件信息,使用SO4 2-的扩散代表矿坑积水反渗造成地下水污染的扩散,利用质量守恒和SO4 2-扩散规律建立数值模型,并对研究三维区域设置网格,结合已知的边界条件和初始条件求解污染源的确定位置坐标,计算结果得出后,在标定的目标处取样,分析数值模拟结果的可信度,并对可信度符合要求区域的串层污染井实施有效封堵或药剂投放,解决了现有的地下水污染评估***中没有结合数值模拟技术反演追踪煤矿串层污染井位置的不足,减轻了大量打井取样的人力投入,提高了煤矿串层污染井治理的精确度。
附图说明
图1是本发明的整体结构图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本实施例提供一种技术方案:一种用于地下水污染的评估***,包括数据采集模块、数值反演追踪模块和污染评估模块,所述数据采集模块用于将研究矿区的地质和水文地质条件以及采样井地下水监测数据传输给所述数值反演追踪模块,所述数值反演追踪模块根据采集的数据信息设定水文地质条件、计算流场边界以及边界已知条件,建立多相流数值模拟模型,所述数值反演追踪模块连接有污染评估模块;
用于地下水污染的评估***通过数值模拟技术反演追踪闭坑煤矿区地下水污染源的方法如下:
步骤一:数据采集模块采集采样井地下水水质和水位信息、矿区的地质条件和水文地质条件以及串层污染井的范围,并将数据传输给数值反演追踪模块;
步骤二:数值反演追踪模块对研究的区域进行网格划分,对串层污染井范围内局部加密,数值反演追踪模块根据数据采集模块提供的数据,基于SO4 2-的扩散建立多相流数值模拟模型,根据数据采集模块提供的数据设定数值模拟初始条件、边界条件,圈定串层污染井的范围,求解污染源的三维坐标;
步骤三:根据数值模拟技术得到三维坐标,在坐标位置打孔取样,验证污染源位置的地下水样本中SO4 2-的浓度是否与数值模拟计算得到的计算浓度接近;
步骤四:在不同的串层污染井采用同样的方法进行数值模拟,收集数值模拟结果和采样结果;
步骤五:对数值模拟的计算结果可信度进行分析,利用数值模拟计算得出SO4 2-浓度减去实际采样地下水中SO4 2-浓度的绝对值,除以实际采样地下水中SO4 2-的浓度,得到的百分比作为可信度评估的依据,所得到的百分比数值大于80%的记为可信样本组,并对可信样本组对应的串层污染井进行有效封堵或者投放药剂进行治理。
通过数据采集模块的设置能够采集闭坑煤矿区域的地质条件信息和水文地质条件信息,使用SO4 2-的扩散代表矿坑积水反渗造成地下水污染的扩散,利用质量守恒和SO4 2-扩散规律建立数值模型,并对研究三维区域设置网格,结合已知的边界条件和初始条件求解污染源的确定位置坐标,计算结果得出后在标定的目标处取样分析数值模拟结果的可信度,并对可信度符合要求区域的串层污染井实施有效封堵或药剂投放,解决了现有的地下水污染评估***中没有结合数值模拟技术反演追踪煤矿串层污染井位置的不足,减轻了大量打井取样的人力投入,提高了煤矿串层污染井治理的精确度。
所述数据采集模块采集的矿区地质条件、水文地质条件以及地下水采样井实时监测数据包括:
实际研究场地的三维区域边界:用于划定数值模拟计算流场的三维计算区域,根据实际的场地大小设定计算的空间步长L,并根据空间步长将实际研究场地的大小划分为的三维区域,式中aL、bL和cL记为实际研究场地的长度、宽度和深度,并将计算区域内的坐标集合记为;
待评估矿区地层性质和分类:用于设定数值计算潜水含水层的类型,并根据大部分的潜水含水层类型设定潜水含水层的性质和参数;
研究矿区四周的地下水流量以及水压信息:设定计算区域的上、下边界和左、右边界类型及参数,设定模拟区域的上边界和下边界均为不透水的零通量边界,分别将上边界和下边界记为和,将左边界和右边界设定为已知压力的边界,并分别将左边界和右边界记为和;
设定监测数据的取样时间间隔:作为数值模拟计算的监测时刻划分依据;
采样井信息:用于在数值模拟计算中对采样井编号并实施数值计算;
采样井水质信息中的SO4 2-浓度:通过采样的方式从待求的污染源位置SO4 2-浓度范围中抽取300组数据代入多相流数值模拟模型,并对其进行反演,取100组接近取样井SO4 2-浓度的数据进行核对,取准确度达到提前设定要求的数据进行统计计算。
所述数值追踪反演模块的数值建模信息包括:
初始条件:包括模拟区的初始压力分布、初始SO4 2-以及含水相以及初始SO4 2-的浓度,初始条件为:
式中:
边界条件:根据水文地质条件设定已知的流量边界和压力边界,根据矿区的暴雨强度及周边河流信息设定源汇项以及SO4 2-浓度边界,边界条件为:
式中:记为地下水中SO4 2-的浓度;记为地下水的实际平均流速,记为SO4 2-的饱和度,记为地下水中SO4 2-的弥散系数张量,n记为边界的法线方向,记为上部边界上的已知压力分布函数,记为左边界的已知压力函数,记为右边界的已知压力函数,记为上部边界SO4 2-的已知浓度函数,记为左边界的SO4 2-已知浓度函数,记为右边界的SO4 2-已知浓度函数;
时间步长:根据采样井的采样时间间隔设定数值模拟的计算时间步长T,并选取需要模拟计算的计算时刻t=10T,11T,12T,13T,14T;
多相流数值模拟计算偏微分方程:根据质量守恒定律建立多相流偏微分方程,多相流偏微分方程为:
式中:s记为组成物质种类,s=1记为组成物质为水,s=2记为组成物质为SO4 2-,记为研究模拟区域地层的孔隙度,记为地下水中SO4 2-的总浓度,记为组成物质的密度,记为地下水中SO4 2-的浓度,记为地下水中SO4 2-的饱和度,记为地下水中SO4 2-的弥散系数张量;
式中:D记为SO4 2-在地下水中的分子扩散系数,记为地层的扭曲系数,记为克罗内克函数,和分别记为SO4 2-在地下水中的横向弥散度和纵向弥散度,记为SO4 2-的渗流速度,和分别表示SO4 2-在i方向和在j方向上的渗流分速度;
所述数值反演追踪模块连接有显示屏,所述显示屏用于显示数值模拟计算结果,并将计算结果显示在矿区三维地图上,标定可信度最高的串层污染井三维坐标,同时显示此处SO4 2-的浓度。
通过显示屏的设置能够将数值计算的结果展示在显示屏上,并结合显示屏上的三维地图具体展示串层污染井的位置,便于研究人员对串层污染井的位置进行分析并讨论封堵方式和投药计量,有助于有效防治地下水污染。
所述污染评估模块对所述数值反演追踪模块的计算结果可信度进行评估,所述可信度计算公式为:
通过可信度公式的设置能够对数值计算的结果进行判别,进而能够对可信度达到要求的计算坐标实施串层污染井治理,避免数值模拟的计算结果不准确造成封堵材料和投放药剂的浪费。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (5)
1.一种用于地下水污染的评估***,包括数据采集模块、数值反演追踪模块和污染评估模块,其特征在于,所述数据采集模块用于将研究矿区的地质和水文地质条件以及采样井地下水监测数据传输给所述数值反演追踪模块,所述数值反演追踪模块根据采集的数据信息设定水文地质条件、计算流场边界以及边界已知条件,建立多相流数值模拟模型,所述数值反演追踪模块连接有污染评估模块;所述数值反演追踪模块的数值建模信息包括:
初始条件:包括模拟区的初始压力分布、初始SO4 2-以及含水相以及初始SO4 2-的浓度,初始条件为:
式中:
边界条件:根据水文地质条件设定已知的流量边界和压力边界,根据矿区的暴雨强度及周边河流信息设定源汇项以及SO4 2-浓度边界,边界条件为:
式中:记为地下水中SO4 2-的浓度;记为地下水的实际平均流速,记为SO4 2-的饱和度,记为地下水中SO4 2-的弥散系数张量,n记为边界的法线方向,记为上部边界上的已知压力分布函数,记为左边界的已知压力函数,记为右边界的已知压力函数,记为上部边界SO4 2-的已知浓度函数,记为左边界的SO4 2-已知浓度函数,记为右边界的SO4 2-已知浓度函数;
时间步长:根据采样井的采样时间间隔设定数值模拟的计算时间步长T,并选取需要模拟计算的计算时刻t=10T,11T,12T,13T,14T;
多相流数值模拟计算偏微分方程:根据质量守恒定律建立多相流偏微分方程,多相流偏微分方程为:
式中:s记为组成物质种类,s=1记为组成物质为水,s=2记为组成物质为SO4 2-,记为研究模拟区域地层的孔隙度,记为地下水中SO4 2-的总浓度,记为组成物质的密度,记为地下水中SO4 2-的浓度,记为地下水中SO4 2-的饱和度,记为地下水中SO4 2-的弥散系数张量;
式中:D记为SO4 2-在地下水中的分子扩散系数,记为地层的扭曲系数,记为克罗内克函数,和分别记为SO4 2-在地下水中的横向弥散度和纵向弥散度,记为SO4 2-的渗流速度,和分别表示SO4 2-在i方向和在j方向上的渗流分速度;
2.根据权利要求1所述的一种用于地下水污染的评估***,其特征在于:所述数据采集模块采集的矿区地质条件、水文地质条件以及地下水采样井实时监测数据包括:
实际研究场地的三维区域边界:用于划定数值模拟计算流场的三维计算区域,根据实际的场地大小设定计算的空间步长L,并根据空间步长将实际研究场地的大小划分为的三维区域,式中aL、bL和cL记为实际研究场地的长度、宽度和深度,并将计算区域内的坐标集合记为;
待评估矿区地层性质和分类:用于设定数值计算潜水含水层的类型,并根据大部分的潜水含水层类型设定潜水含水层的性质和参数;
研究矿区四周的地下水流量以及水压信息:设定计算区域的上、下边界和左、右边界类型及参数,设定模拟区域的上边界和下边界均为不透水的零通量边界,分别将上边界和下边界记为和,将左边界和右边界设定为已知压力的边界,并分别将左边界和右边界记为和;
设定监测数据的取样时间间隔:作为数值模拟计算的监测时刻划分依据;
采样井信息:用于在数值模拟计算中对采样井编号并实施数值计算;
采样井水质信息中的SO4 2-浓度:通过采样的方式从待求的污染源位置SO4 2-浓度范围中抽取300组数据代入多相流数值模拟模型,并对其进行反演,取100组接近取样井SO4 2-浓度的数据进行核对,取准确度达到提前设定要求的数据进行统计计算。
3.根据权利要求1所述的一种用于地下水污染的评估***,其特征在于:所述数值反演追踪模块上设有显示屏,所述显示屏用于显示数值模拟计算结果,并将计算结果显示在矿区三维地图上,标定可信度最高的串层污染井三维坐标,同时显示此处SO4 2-的浓度。
5.根据权利要求1所述的一种用于地下水污染的评估***,其特征在于:用于地下水污染的评估***通过数值模拟技术反演追踪闭坑煤矿区地下水污染源的方法如下:
步骤一:数据采集模块采集采样井地下水水质和水位信息、矿区的地质条件和水文地质条件以及串层污染井的范围,并将数据传输给数值反演追踪模块;
步骤二:数值反演追踪模块对研究的区域进行网格划分,对串层污染井范围内局部加密,数值反演追踪模块根据数据采集模块提供的数据,基于SO4 2-的扩散建立多相流数值模拟模型,根据数据采集模块提供的数据设定数值模拟初始条件、边界条件,圈定串层污染井的范围,求解污染源的三维坐标;
步骤三:根据数值模拟技术得到三维坐标,在坐标位置打孔取样,验证污染源位置的地下水样本中SO4 2-的浓度是否与数值模拟计算得到的计算浓度接近;
步骤四:在不同的串层污染井采用同样的方法进行数值模拟,收集数值模拟结果和采样结果;
步骤五:对数值模拟的计算结果可信度进行分析,利用数值模拟计算得出SO4 2-浓度减去实际采样地下水中SO4 2-浓度的绝对值,除以实际采样地下水中SO4 2-的浓度,得到的百分比作为可信度评估的依据,所得到的百分比数值大于80%的记为可信样本组,并对可信样本组对应的串层污染井进行有效封堵或者投放药剂进行治理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211256594.3A CN115329607B (zh) | 2022-10-14 | 2022-10-14 | 一种用于地下水污染的评估***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211256594.3A CN115329607B (zh) | 2022-10-14 | 2022-10-14 | 一种用于地下水污染的评估***及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115329607A true CN115329607A (zh) | 2022-11-11 |
CN115329607B CN115329607B (zh) | 2023-02-03 |
Family
ID=83914649
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211256594.3A Active CN115329607B (zh) | 2022-10-14 | 2022-10-14 | 一种用于地下水污染的评估***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115329607B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116611274A (zh) * | 2023-07-21 | 2023-08-18 | 中南大学 | 一种地下水污染运移可视化数值模拟方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109190280A (zh) * | 2018-09-18 | 2019-01-11 | 东北农业大学 | 一种基于核极限学习机替代模型的地下水污染源反演识别方法 |
US20200252283A1 (en) * | 2017-10-16 | 2020-08-06 | Chinese Research Academy Of Environmental Sciences | Multi-objective optimization method for groundwater pollution monitoring network |
CN112149871A (zh) * | 2020-08-25 | 2020-12-29 | 河海大学 | 一种基于gis空间统计与随机模拟相结合的污染物点源解析方法 |
CN112697849A (zh) * | 2020-11-20 | 2021-04-23 | 中科盛联(北京)新材料有限公司 | 动态污染源定位方法 |
CN113239598A (zh) * | 2021-06-08 | 2021-08-10 | 中国环境科学研究院 | 一种应用数值模拟的地下水污染源空间综合识别方法 |
WO2021208393A1 (zh) * | 2020-04-15 | 2021-10-21 | 北京工业大学 | 一种大气污染物排放清单的反演估算方法 |
CN114254572A (zh) * | 2021-12-16 | 2022-03-29 | 西北工业大学太仓长三角研究院 | 考虑污染物沉积的航发压气机流场性能预测方法及*** |
CN114386329A (zh) * | 2022-01-11 | 2022-04-22 | 南京大学 | 一种不确定性条件下dnapl污染场地修复的多目标优化方法 |
CN114492164A (zh) * | 2021-12-24 | 2022-05-13 | 吉林大学 | 基于多核极限学习机的有机污染物迁移数值模型替代方法 |
CN114707395A (zh) * | 2021-12-24 | 2022-07-05 | 吉林大学 | 基于超启发式-同伦算法的地下水有机污染源反演方法 |
CN114861502A (zh) * | 2022-05-27 | 2022-08-05 | 北京林业大学 | 基于Modflow模型的三维动态地下水污染模拟的安全饮用水区域确定方法 |
-
2022
- 2022-10-14 CN CN202211256594.3A patent/CN115329607B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200252283A1 (en) * | 2017-10-16 | 2020-08-06 | Chinese Research Academy Of Environmental Sciences | Multi-objective optimization method for groundwater pollution monitoring network |
CN109190280A (zh) * | 2018-09-18 | 2019-01-11 | 东北农业大学 | 一种基于核极限学习机替代模型的地下水污染源反演识别方法 |
WO2021208393A1 (zh) * | 2020-04-15 | 2021-10-21 | 北京工业大学 | 一种大气污染物排放清单的反演估算方法 |
CN112149871A (zh) * | 2020-08-25 | 2020-12-29 | 河海大学 | 一种基于gis空间统计与随机模拟相结合的污染物点源解析方法 |
CN112697849A (zh) * | 2020-11-20 | 2021-04-23 | 中科盛联(北京)新材料有限公司 | 动态污染源定位方法 |
CN113239598A (zh) * | 2021-06-08 | 2021-08-10 | 中国环境科学研究院 | 一种应用数值模拟的地下水污染源空间综合识别方法 |
CN114254572A (zh) * | 2021-12-16 | 2022-03-29 | 西北工业大学太仓长三角研究院 | 考虑污染物沉积的航发压气机流场性能预测方法及*** |
CN114492164A (zh) * | 2021-12-24 | 2022-05-13 | 吉林大学 | 基于多核极限学习机的有机污染物迁移数值模型替代方法 |
CN114707395A (zh) * | 2021-12-24 | 2022-07-05 | 吉林大学 | 基于超启发式-同伦算法的地下水有机污染源反演方法 |
CN114386329A (zh) * | 2022-01-11 | 2022-04-22 | 南京大学 | 一种不确定性条件下dnapl污染场地修复的多目标优化方法 |
CN114861502A (zh) * | 2022-05-27 | 2022-08-05 | 北京林业大学 | 基于Modflow模型的三维动态地下水污染模拟的安全饮用水区域确定方法 |
Non-Patent Citations (5)
Title |
---|
GAO QI等: "Source Inversion of a Sudden Pollution Accident Occurred in Plain River Network Based on Pattern Search Method", 《IEEE XPLORE》 * |
侯泽宇: "地下水DNAPLs污染源反演识别及不确定性分析研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
杨之珍: "三维非稳定流数学模型在模拟废弃煤矿区地下水串层污染中的应用研究", 《能源与环境》 * |
杨佳立: "地下水水源保护区边界污染源防控的V-M数值模拟探析", 《科技广场》 * |
闫雪嫚: "基于贝叶斯理论的地下水DNAPLs污染源反演识别研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116611274A (zh) * | 2023-07-21 | 2023-08-18 | 中南大学 | 一种地下水污染运移可视化数值模拟方法 |
CN116611274B (zh) * | 2023-07-21 | 2023-09-29 | 中南大学 | 一种地下水污染运移可视化数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115329607B (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Desaulniers et al. | Origin, age and movement of pore water in argillaceous Quaternary deposits at four sites in southwestern Ontario | |
Banks et al. | Assessing spatial and temporal connectivity between surface water and groundwater in a regional catchment: Implications for regional scale water quantity and quality | |
Schneider et al. | Towards improved instrumentation for assessing river-groundwater interactions in a restored river corridor | |
Berhane et al. | Water leakage investigation of micro-dam reservoirs in Mesozoic sedimentary sequences in Northern Ethiopia | |
Scheiber et al. | Hydrochemical apportioning of irrigation groundwater sources in an alluvial aquifer | |
Bonanno et al. | Flow directions of stream‐groundwater exchange in a headwater catchment during the hydrologic year | |
CN115329607B (zh) | 一种用于地下水污染的评估***及方法 | |
Ronchetti et al. | Tracer test to assess flow and transport parameters of an earth slide: The Montecagno landslide case study (Italy) | |
Li et al. | Groundwater circulation patterns in bedrock aquifers from a pre-selected area of high-level radioactive waste repository based on two-dimensional numerical simulation | |
Bourke et al. | A hydrological framework for persistent pools along non-perennial rivers | |
Rushton et al. | Significance of hydraulic head gradients within horizontal wells in unconfined aquifers of limited saturated thickness | |
Amiri et al. | Hydrogeological assessment and estimation of groundwater inflow into the water transmission tunnel to Urmia Lake, Northwestern Iran | |
CN109978273B (zh) | 一种基于地貌控水的矿井涌水量预测方法 | |
Vorlicek et al. | Quantitative hydrogeological studies of the Treviso alluvial plain, NE Italy | |
Michalski et al. | Characterization of transmissive fractures by simple tracing of in‐well flow | |
Hauwert | Groundwater flow and recharge within the Barton Springs segment of the Edwards Aquifer, southern Travis and northern Hays counties, Texas | |
Chiu et al. | The effect of hydrological conditions and bioactivities on the spatial and temporal variations of streambed hydraulic characteristics at the subtropical alpine catchment | |
Rickel et al. | Seasonal shifts in surface water-groundwater connections in a ferricrete-impacted stream estimated from electrical resistivity | |
Belan | Characterizating a fractured rock aquifer with hydraulic testing at a contaminated municipal well using flexible liner methods and depth discrete monitoring | |
YOKOYAMA et al. | A review of groundwater observation methods for slow-moving landslide | |
Wilkes et al. | Characterisation of the hydrogeology of the Augustus River catchment, Western Australia | |
Middleton | Aquifer–stream connectivity at various scales: Application of sediment–water interface temperature and vulnerability assessments of groundwater dependent streams | |
Schalla et al. | Interim characterization report for the 300 Area process trenches | |
Lee et al. | Investigating the origin and dynamics of salinity in a confined aquifer system in southeast Australia (Western Port Basin) | |
Mahindawansha et al. | Evaluation of Hydraulic Exchange in the Hyporheic Zone: A Depth-Oriented Analysis Focusing on Upstream and Downstream Conditions |
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 |