CN115329607B - 一种用于地下水污染的评估***及方法 - Google Patents

一种用于地下水污染的评估***及方法 Download PDF

Info

Publication number
CN115329607B
CN115329607B CN202211256594.3A CN202211256594A CN115329607B CN 115329607 B CN115329607 B CN 115329607B CN 202211256594 A CN202211256594 A CN 202211256594A CN 115329607 B CN115329607 B CN 115329607B
Authority
CN
China
Prior art keywords
groundwater
boundary
numerical
concentration
numerical simulation
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
Application number
CN202211256594.3A
Other languages
English (en)
Other versions
CN115329607A (zh
Inventor
张立川
王超
邹晔
秦志强
陈洪年
张宇飞
薛良方
朱世芳
张梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Lunan Geological Engineering Survey Institute of Second Geological Brigade of Shandong Geological Survey Bureau
Original Assignee
Shandong Lunan Geological Engineering Survey Institute of Second Geological Brigade of Shandong Geological Survey Bureau
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong Lunan Geological Engineering Survey Institute of Second Geological Brigade of Shandong Geological Survey Bureau filed Critical Shandong Lunan Geological Engineering Survey Institute of Second Geological Brigade of Shandong Geological Survey Bureau
Priority to CN202211256594.3A priority Critical patent/CN115329607B/zh
Publication of CN115329607A publication Critical patent/CN115329607A/zh
Application granted granted Critical
Publication of CN115329607B publication Critical patent/CN115329607B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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，并根据空间步长将实际研究场地的大小划分为
Figure DEST_PATH_IMAGE002
的三维区域，式中aL、bL和cL记为实际研究场地的长度、宽度和深度，并将计算区域内的坐标集合记为
Figure DEST_PATH_IMAGE004
；
待评估矿区地层性质和分类：用于设定数值计算潜水含水层的类型，并根据大部分的潜水含水层类型设定潜水含水层的性质和参数；
地表类型：用于设定数值模拟计算的上部边界，根据场地的开阔性和雨水渗入特性将上部边界设定为水量交换边界，记为
Figure DEST_PATH_IMAGE006
；
底层地下水的边界岩层性质：用于设定数值模拟计算区域的底部边界性质和参数，根据底部微渗透的特性，将研究区域的底部边界设定为零通量边界，记为
Figure DEST_PATH_IMAGE008
；
研究矿区四周的地下水流量以及水压信息：设定计算区域的上、下边界和左、右边界类型及参数，设定模拟区域的上边界和下边界均为不透水的零通量边界，分别将上边界和下边界记为
Figure DEST_PATH_IMAGE010
和
Figure DEST_PATH_IMAGE012
，将左边界和右边界设定为已知压力的边界，并分别将左边界和右边界记为
Figure DEST_PATH_IMAGE014
和
Figure DEST_PATH_IMAGE016
；
设定监测数据的取样时间间隔：作为数值模拟计算的监测时刻划分依据；
采样井信息：用于在数值模拟计算中对采样井编号并实施数值计算；
采样井水质信息中的SO4 2-浓度：通过采样的方式从待求的污染源位置SO4 2-浓度范围中抽取300组数据代入多相流数值模拟模型，并对其进行反演，取100组接近取样井SO4 2-浓度的数据进行核对，取准确度达到提前设定要求的数据进行统计计算。
作为本发明进一步的方案，所述数值追踪反演模块的数值建模信息包括：
初始条件：包括模拟区的初始压力分布、初始SO4 2-以及含水相以及初始SO4 2-的浓度，初始条件为：
Figure DEST_PATH_IMAGE018
；
Figure DEST_PATH_IMAGE020
；
Figure DEST_PATH_IMAGE022
；
式中：
Figure DEST_PATH_IMAGE024
记为模拟区域的初始压力分布函数，
Figure DEST_PATH_IMAGE026
记为已知的压力函数；
Figure DEST_PATH_IMAGE028
记为地下水的初始饱和度，
Figure DEST_PATH_IMAGE030
记为已知的地下水饱和度函数；
Figure DEST_PATH_IMAGE032
记为地下水中SO4 2-的初始浓度函数，
Figure DEST_PATH_IMAGE034
记为已知的SO4 2-初始浓度函数；
边界条件：根据水文地质条件设定已知的流量边界和压力边界，根据矿区的暴雨强度及周边河流信息设定源汇项以及SO4 2-浓度边界，边界条件为：
Figure DEST_PATH_IMAGE036
；
Figure DEST_PATH_IMAGE038
；
Figure DEST_PATH_IMAGE040
；
Figure DEST_PATH_IMAGE042
；
Figure DEST_PATH_IMAGE044
；
Figure DEST_PATH_IMAGE046
；
Figure DEST_PATH_IMAGE048
；
Figure DEST_PATH_IMAGE050
；
式中：
Figure DEST_PATH_IMAGE052
记为地下水中SO4 2-的浓度；
Figure DEST_PATH_IMAGE054
记为地下水的实际平均流速，
Figure DEST_PATH_IMAGE056
记为SO4 2-的饱和度，
Figure DEST_PATH_IMAGE058
记为地下水中SO4 2-的弥散系数张量，n记为边界的法线方向，
Figure DEST_PATH_IMAGE060
记为上部边界上的已知压力分布函数，
Figure DEST_PATH_IMAGE062
记为左边界的已知压力函数，
Figure DEST_PATH_IMAGE064
记为右边界的已知压力函数，
Figure DEST_PATH_IMAGE066
记为上部边界SO4 2-的已知浓度函数，
Figure DEST_PATH_IMAGE068
记为左边界的SO4 2-已知浓度函数，
Figure 640785DEST_PATH_IMAGE066
记为右边界的SO4 2-已知浓度函数；
时间步长：根据采样井的采样时间间隔设定数值模拟的计算时间步长T，并选取需要模拟计算的计算时刻t=10T，11T，12T，13T，14T；
多相流数值模拟计算偏微分方程：根据质量守恒定律建立多相流偏微分方程，多相流偏微分方程为：
Figure DEST_PATH_IMAGE070
，
Figure DEST_PATH_IMAGE072
；
式中：s记为组成物质种类，s=1记为组成物质为水，s=2记为组成物质为SO4 2-，
Figure DEST_PATH_IMAGE074
记为研究模拟区域地层的孔隙度，
Figure DEST_PATH_IMAGE076
记为地下水中SO4 2-的总浓度，
Figure DEST_PATH_IMAGE078
记为组成物质的密度，
Figure 571569DEST_PATH_IMAGE052
记为地下水中SO4 2-的浓度，
Figure DEST_PATH_IMAGE080
记为地下水中SO4 2-的饱和度，
Figure 776286DEST_PATH_IMAGE058
记为地下水中SO4 2-的弥散系数张量；
Figure DEST_PATH_IMAGE082
，
Figure 647290DEST_PATH_IMAGE072
；
式中：D记为SO4 2-在地下水中的分子扩散系数，
Figure DEST_PATH_IMAGE084
记为地层的扭曲系数，
Figure DEST_PATH_IMAGE086
记为克罗内克函数，
Figure DEST_PATH_IMAGE088
和
Figure DEST_PATH_IMAGE090
分别记为SO4 2-在地下水中的横向弥散度和纵向弥散度，
Figure DEST_PATH_IMAGE092
记为SO4 2-的渗流速度，
Figure DEST_PATH_IMAGE094
和
Figure DEST_PATH_IMAGE096
分别表示SO4 2-在i方向和在j方向上的渗流分速度；
Figure DEST_PATH_IMAGE098
，
Figure 733145DEST_PATH_IMAGE072
；
式中：
Figure DEST_PATH_IMAGE100
记为SO4 2-的相对渗透率，
Figure DEST_PATH_IMAGE102
记为SO4 2-的固有渗透率，
Figure DEST_PATH_IMAGE104
记为地下水的粘滞度，
Figure DEST_PATH_IMAGE106
记为地下水的压力，g记为重力加速度，
Figure DEST_PATH_IMAGE108
记为地下水的密度。
作为本发明进一步的方案，所述数值反演追踪模块连接有显示屏，所述显示屏用于显示数值模拟计算结果，并将计算结果显示在矿区三维地图上，标定可信度最高的串层污染井三维坐标，同时显示此处SO4 2-的浓度。
作为本发明进一步的方案，所述污染评估模块对所述数值反演追踪模块的计算结果可信度进行评估，所述可信度计算公式为：
Figure DEST_PATH_IMAGE110
；
式中：
Figure DEST_PATH_IMAGE112
记为数据可信度，
Figure DEST_PATH_IMAGE114
记为通过数值计算所得的SO4 2-浓度，
Figure DEST_PATH_IMAGE116
记为实际取样所得SO4 2-浓度。
本发明相比现有技术具有以下优点：该用于地下水污染的评估***及方法，通过数据采集模块的设置能够采集闭坑煤矿区域的地质条件信息和水文地质条件信息，使用SO4 2-的扩散代表矿坑积水反渗造成地下水污染的扩散，利用质量守恒和SO4 2-扩散规律建立数值模型，并对研究三维区域设置网格，结合已知的边界条件和初始条件求解污染源的确定位置坐标，计算结果得出后，在标定的目标处取样，分析数值模拟结果的可信度，并对可信度符合要求区域的串层污染井实施有效封堵或药剂投放，解决了现有的地下水污染评估***中没有结合数值模拟技术反演追踪煤矿串层污染井位置的不足，减轻了大量打井取样的人力投入，提高了煤矿串层污染井治理的精确度。
附图说明
图1是本发明的整体结构图。
具体实施方式
下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
如图1所示，本实施例提供一种技术方案：一种用于地下水污染的评估***，包括数据采集模块、数值反演追踪模块和污染评估模块，所述数据采集模块用于将研究矿区的地质和水文地质条件以及采样井地下水监测数据传输给所述数值反演追踪模块，所述数值反演追踪模块根据采集的数据信息设定水文地质条件、计算流场边界以及边界已知条件，建立多相流数值模拟模型，所述数值反演追踪模块连接有污染评估模块；
用于地下水污染的评估***通过数值模拟技术反演追踪闭坑煤矿区地下水污染源的方法如下：
步骤一：数据采集模块采集采样井地下水水质和水位信息、矿区的地质条件和水文地质条件以及串层污染井的范围，并将数据传输给数值反演追踪模块；
步骤二：数值反演追踪模块对研究的区域进行网格划分，对串层污染井范围内局部加密，数值反演追踪模块根据数据采集模块提供的数据，基于SO4 2-的扩散建立多相流数值模拟模型，根据数据采集模块提供的数据设定数值模拟初始条件、边界条件，圈定串层污染井的范围，求解污染源的三维坐标；
步骤三：根据数值模拟技术得到三维坐标，在坐标位置打孔取样，验证污染源位置的地下水样本中SO4 2-的浓度是否与数值模拟计算得到的计算浓度接近；
步骤四：在不同的串层污染井采用同样的方法进行数值模拟，收集数值模拟结果和采样结果；
步骤五：对数值模拟的计算结果可信度进行分析，利用数值模拟计算得出SO4 2-浓度减去实际采样地下水中SO4 2-浓度的绝对值，除以实际采样地下水中SO4 2-的浓度，得到的百分比作为可信度评估的依据，所得到的百分比数值大于80%的记为可信样本组，并对可信样本组对应的串层污染井进行有效封堵或者投放药剂进行治理。
通过数据采集模块的设置能够采集闭坑煤矿区域的地质条件信息和水文地质条件信息，使用SO4 2-的扩散代表矿坑积水反渗造成地下水污染的扩散，利用质量守恒和SO4 2-扩散规律建立数值模型，并对研究三维区域设置网格，结合已知的边界条件和初始条件求解污染源的确定位置坐标，计算结果得出后在标定的目标处取样分析数值模拟结果的可信度，并对可信度符合要求区域的串层污染井实施有效封堵或药剂投放，解决了现有的地下水污染评估***中没有结合数值模拟技术反演追踪煤矿串层污染井位置的不足，减轻了大量打井取样的人力投入，提高了煤矿串层污染井治理的精确度。
所述数据采集模块采集的矿区地质条件、水文地质条件以及地下水采样井实时监测数据包括：
实际研究场地的三维区域边界：用于划定数值模拟计算流场的三维计算区域，根据实际的场地大小设定计算的空间步长L，并根据空间步长将实际研究场地的大小划分为
Figure 523247DEST_PATH_IMAGE002
的三维区域，式中aL、bL和cL记为实际研究场地的长度、宽度和深度，并将计算区域内的坐标集合记为
Figure 543155DEST_PATH_IMAGE004
；
待评估矿区地层性质和分类：用于设定数值计算潜水含水层的类型，并根据大部分的潜水含水层类型设定潜水含水层的性质和参数；
地表类型：用于设定数值模拟计算的上部边界，根据场地的开阔性和雨水渗入特性将上部边界设定为水量交换边界，记为
Figure 749009DEST_PATH_IMAGE006
；
底层地下水的边界岩层性质：用于设定数值模拟计算区域的底部边界性质和参数，根据底部微渗透的特性，将研究区域的底部边界设定为零通量边界，记为
Figure 239771DEST_PATH_IMAGE008
；
研究矿区四周的地下水流量以及水压信息：设定计算区域的上、下边界和左、右边界类型及参数，设定模拟区域的上边界和下边界均为不透水的零通量边界，分别将上边界和下边界记为
Figure 872877DEST_PATH_IMAGE010
和
Figure 380082DEST_PATH_IMAGE012
，将左边界和右边界设定为已知压力的边界，并分别将左边界和右边界记为
Figure 170053DEST_PATH_IMAGE014
和
Figure 16786DEST_PATH_IMAGE016
；
设定监测数据的取样时间间隔：作为数值模拟计算的监测时刻划分依据；
采样井信息：用于在数值模拟计算中对采样井编号并实施数值计算；
采样井水质信息中的SO4 2-浓度：通过采样的方式从待求的污染源位置SO4 2-浓度范围中抽取300组数据代入多相流数值模拟模型，并对其进行反演，取100组接近取样井SO4 2-浓度的数据进行核对，取准确度达到提前设定要求的数据进行统计计算。
所述数值追踪反演模块的数值建模信息包括：
初始条件：包括模拟区的初始压力分布、初始SO4 2-以及含水相以及初始SO4 2-的浓度，初始条件为：
Figure 24056DEST_PATH_IMAGE018
；
Figure 487398DEST_PATH_IMAGE020
；
Figure 769475DEST_PATH_IMAGE022
；
式中：
Figure 533032DEST_PATH_IMAGE024
记为模拟区域的初始压力分布函数，
Figure 773520DEST_PATH_IMAGE026
记为已知的压力函数；
Figure 691536DEST_PATH_IMAGE028
记为地下水的初始饱和度，
Figure 777303DEST_PATH_IMAGE030
记为已知的地下水饱和度函数；
Figure 395366DEST_PATH_IMAGE032
记为地下水中SO4 2-的初始浓度函数，
Figure 806756DEST_PATH_IMAGE034
记为已知的SO4 2-初始浓度函数；
边界条件：根据水文地质条件设定已知的流量边界和压力边界，根据矿区的暴雨强度及周边河流信息设定源汇项以及SO4 2-浓度边界，边界条件为：
Figure 447953DEST_PATH_IMAGE036
；
Figure 337412DEST_PATH_IMAGE038
；
Figure 75561DEST_PATH_IMAGE040
；
Figure 861114DEST_PATH_IMAGE042
；
Figure 51924DEST_PATH_IMAGE044
；
Figure 745073DEST_PATH_IMAGE046
；
Figure 508368DEST_PATH_IMAGE048
；
Figure 58298DEST_PATH_IMAGE050
；
式中：
Figure 939666DEST_PATH_IMAGE052
记为地下水中SO4 2-的浓度；
Figure 170927DEST_PATH_IMAGE054
记为地下水的实际平均流速，
Figure 821351DEST_PATH_IMAGE056
记为SO4 2-的饱和度，
Figure 479866DEST_PATH_IMAGE058
记为地下水中SO4 2-的弥散系数张量，n记为边界的法线方向，
Figure 645268DEST_PATH_IMAGE060
记为上部边界上的已知压力分布函数，
Figure 414641DEST_PATH_IMAGE062
记为左边界的已知压力函数，
Figure 185151DEST_PATH_IMAGE064
记为右边界的已知压力函数，
Figure 811304DEST_PATH_IMAGE066
记为上部边界SO4 2-的已知浓度函数，
Figure 667265DEST_PATH_IMAGE068
记为左边界的SO4 2-已知浓度函数，
Figure 762301DEST_PATH_IMAGE066
记为右边界的SO4 2-已知浓度函数；
时间步长：根据采样井的采样时间间隔设定数值模拟的计算时间步长T，并选取需要模拟计算的计算时刻t=10T，11T，12T，13T，14T；
多相流数值模拟计算偏微分方程：根据质量守恒定律建立多相流偏微分方程，多相流偏微分方程为：
Figure 121738DEST_PATH_IMAGE070
，
Figure 387635DEST_PATH_IMAGE072
；
式中：s记为组成物质种类，s=1记为组成物质为水，s=2记为组成物质为SO4 2-，
Figure 262050DEST_PATH_IMAGE074
记为研究模拟区域地层的孔隙度，
Figure 904384DEST_PATH_IMAGE076
记为地下水中SO4 2-的总浓度，
Figure 118327DEST_PATH_IMAGE078
记为组成物质的密度，
Figure 351863DEST_PATH_IMAGE052
记为地下水中SO4 2-的浓度，
Figure 916836DEST_PATH_IMAGE080
记为地下水中SO4 2-的饱和度，
Figure 97282DEST_PATH_IMAGE058
记为地下水中SO4 2-的弥散系数张量；
Figure 165732DEST_PATH_IMAGE082
，
Figure 271966DEST_PATH_IMAGE072
；
式中：D记为SO4 2-在地下水中的分子扩散系数，
Figure 120973DEST_PATH_IMAGE084
记为地层的扭曲系数，
Figure 839530DEST_PATH_IMAGE086
记为克罗内克函数，
Figure 28066DEST_PATH_IMAGE088
和
Figure 806666DEST_PATH_IMAGE090
分别记为SO4 2-在地下水中的横向弥散度和纵向弥散度，
Figure 142970DEST_PATH_IMAGE092
记为SO4 2-的渗流速度，
Figure 665218DEST_PATH_IMAGE094
和
Figure 708260DEST_PATH_IMAGE096
分别表示SO4 2-在i方向和在j方向上的渗流分速度；
Figure 657762DEST_PATH_IMAGE098
，
Figure 215782DEST_PATH_IMAGE072
；
式中：
Figure 541721DEST_PATH_IMAGE100
记为SO4 2-的相对渗透率，
Figure 203385DEST_PATH_IMAGE102
记为SO4 2-的固有渗透率，
Figure 58208DEST_PATH_IMAGE104
记为地下水的粘滞度，
Figure 572366DEST_PATH_IMAGE106
记为地下水的压力，g记为重力加速度，
Figure 233155DEST_PATH_IMAGE108
记为地下水的密度。
所述数值反演追踪模块连接有显示屏，所述显示屏用于显示数值模拟计算结果，并将计算结果显示在矿区三维地图上，标定可信度最高的串层污染井三维坐标，同时显示此处SO4 2-的浓度。
通过显示屏的设置能够将数值计算的结果展示在显示屏上，并结合显示屏上的三维地图具体展示串层污染井的位置，便于研究人员对串层污染井的位置进行分析并讨论封堵方式和投药计量，有助于有效防治地下水污染。
所述污染评估模块对所述数值反演追踪模块的计算结果可信度进行评估，所述可信度计算公式为：
Figure 250789DEST_PATH_IMAGE110
；
式中：
Figure 542093DEST_PATH_IMAGE112
记为数据可信度，
Figure 543547DEST_PATH_IMAGE114
记为通过数值计算所得的SO4 2-浓度，
Figure 945710DEST_PATH_IMAGE116
记为实际取样所得SO4 2-浓度。
通过可信度公式的设置能够对数值计算的结果进行判别，进而能够对可信度达到要求的计算坐标实施串层污染井治理，避免数值模拟的计算结果不准确造成封堵材料和投放药剂的浪费。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (4)

1.一种用于地下水污染的评估***，包括数据采集模块、数值反演追踪模块和污染评估模块，其特征在于，所述数据采集模块用于将研究矿区的地质和水文地质条件以及采样井地下水监测数据传输给所述数值反演追踪模块，所述数值反演追踪模块根据采集的数据信息设定水文地质条件、计算流场边界以及边界已知条件，建立多相流数值模拟模型，所述数值反演追踪模块连接有污染评估模块；所述数值反演追踪模块的数值建模信息包括：
初始条件：包括模拟区的初始压力分布、初始SO4 2-以及含水相以及初始SO4 2-的浓度，初始条件为：
P(x，y，z，t)|t＝0＝f(x，y，z)，(x，y，z)∈Ψ；
Sl(x，y，z，t)|t＝0＝g(x，y，z)，(x，y，z)∈Ψ；
Csl(x，y，z，t)|t＝0＝h(x，y，z)，(x，y，z)∈Ψ；
式中：
P(x，y，z，t)|t＝0记为模拟区域的初始压力分布函数，f(x，y，z)记为已知的压力函数；
Sl(x，y，z，t)|t＝0记为地下水的初始饱和度，g(x，y，z)记为已知的地下水饱和度函数；
Csl(x，y，z，t)|t＝0记为地下水中SO4 2-的初始浓度函数，h(x，y，z)记为已知的SO4 2-初始浓度函数；
边界条件：根据水文地质条件设定已知的流量边界和压力边界，根据矿区的暴雨强度及周边河流信息设定源汇项以及SO4 2-浓度边界，边界条件为：
Figure FDA0004016170510000021
Figure FDA0004016170510000022
P(x，y，z，t)|Γzb＝f2(x，y，z，t)，(x，y，z)∈Γzb，t≥0；
Figure FDA0004016170510000028
Figure FDA0004016170510000023
Figure FDA0004016170510000024
Figure FDA0004016170510000025
Figure FDA0004016170510000026
式中：Csl记为地下水中SO4 2-的浓度；u1记为地下水的实际平均流速，S1记为SO4 2-的饱和度，
Figure FDA0004016170510000027
记为地下水中SO4 2-的弥散系数张量，n记为边界的法线方向，f1(x，y，z，t)记为上部边界上的已知压力分布函数，f2(x，y，z，t)记为左边界的已知压力函数，f3(x，y，z，t)记为右边界的已知压力函数，g1(x，y，z，t)记为上部边界SO4 2-的已知浓度函数，g2(x，y，z，t)记为左边界的SO4 2-已知浓度函数，g3(x，y，z，t)记为右边界的SO4 2-已知浓度函数；
时间步长：根据采样井的采样时间间隔设定数值模拟的计算时间步长T，并选取需要模拟计算的计算时刻t＝10T，11T，12T，13T，14T；
多相流数值模拟计算偏微分方程：根据质量守恒定律建立多相流偏微分方程，多相流偏微分方程为：
Figure FDA0004016170510000031
式中：l记为组成物质种类，l＝1记为组成物质为水，l＝2记为组成物质为SO4 2-，
Figure FDA0004016170510000032
记为研究模拟区域地层的孔隙度，
Figure FDA0004016170510000033
记为地下水中SO4 2-的总浓度，ρsl记为组成物质的密度，Csl记为地下水中SO4 2-的浓度，Sl记为地下水中SO4 2-的饱和度，
Figure FDA0004016170510000034
记为地下水中SO4 2-的弥散系数张量；
Figure FDA0004016170510000035
式中：D记为SO4 2-在地下水中的分子扩散系数，τ记为地层的扭曲系数，δij记为克罗内克函数，ahl和azl分别记为SO4 2-在地下水中的横向弥散度和纵向弥散度，vl记为SO4 2-的渗流速度，vi和vj分别表示SO4 2-在i方向和在j方向上的渗流分速度；
Figure FDA0004016170510000036
式中：Krl记为SO4 2-的相对渗透率，
Figure FDA0004016170510000041
记为SO4 2-的固有渗透率，μl记为地下水的粘滞度，Pl记为地下水的压力，g记为重力加速度，ρl记为地下水的密度；
所述数据采集模块采集的矿区地质条件、水文地质条件以及地下水采样井实时监测数据包括：
实际研究场地的三维区域边界：用于划定数值模拟计算流场的三维计算区域，根据实际的场地大小设定计算的空间步长L，并根据空间步长将实际研究场地的大小划分为aL×bL×cL的三维区域，式中aL、bL和cL记为实际研究场地的长度、宽度和深度，并将计算区域内的坐标集合记为Ψ；
待评估矿区地层性质和分类：用于设定数值计算潜水含水层的类型，并根据大部分的潜水含水层类型设定潜水含水层的性质和参数；
地表类型：用于设定数值模拟计算的上部边界，根据场地的开阔性和雨水渗入特性将上部边界设定为水量交换边界，记为Γs；
底层地下水的边界岩层性质：用于设定数值模拟计算区域的底部边界性质和参数，根据底部微渗透的特性，将研究区域的底部边界设定为零通量边界，记为Γd；
研究矿区四周的地下水流量以及水压信息：设定计算区域的上、下边界和左、右边界类型及参数，设定模拟区域的上边界和下边界均为不透水的零通量边界，分别将上边界和下边界记为Γsb和Γxb，将左边界和右边界设定为已知压力的边界，并分别将左边界和右边界记为Γzb和Γyb；
设定监测数据的取样时间间隔：作为数值模拟计算的监测时刻划分依据；
采样井信息：用于在数值模拟计算中对采样井编号并实施数值计算；
采样井水质信息中的SO4 2-浓度：通过采样的方式从待求的污染源位置SO4 2-浓度范围中抽取300组数据代入多相流数值模拟模型，并对其进行反演，取100组接近取样井SO4 2-浓度的数据进行核对，取准确度达到提前设定要求的数据进行统计计算。
2.根据权利要求1所述的一种用于地下水污染的评估***，其特征在于：所述数值反演追踪模块上设有显示屏，所述显示屏用于显示数值模拟计算结果，并将计算结果显示在矿区三维地图上，标定可信度最高的串层污染井三维坐标，同时显示此处SO4 2-的浓度。
3.根据权利要求1所述的一种用于地下水污染的评估***，其特征在于：所述污染评估模块对所述数值反演追踪模块的计算结果可信度进行评估，所述可信度计算公式为：
Figure FDA0004016170510000051
式中：QX记为数据可信度，CC记为通过数值计算所得的SO4 2-浓度，CR记为实际取样所得SO4 2-浓度。
4.根据权利要求1所述的一种用于地下水污染的评估***，其特征在于：用于地下水污染的评估***通过数值模拟技术反演追踪闭坑煤矿区地下水污染源的方法如下：
步骤一：数据采集模块采集采样井地下水水质和水位信息、矿区的地质条件和水文地质条件以及串层污染井的范围，并将数据传输给数值反演追踪模块；
步骤二：数值反演追踪模块对研究的区域进行网格划分，对串层污染井范围内局部加密，数值反演追踪模块根据数据采集模块提供的数据，基于SO4 2-的扩散建立多相流数值模拟模型，根据数据采集模块提供的数据设定数值模拟初始条件、边界条件，圈定串层污染井的范围，求解污染源的三维坐标；
步骤三：根据数值模拟技术得到三维坐标，在坐标位置打孔取样，验证污染源位置的地下水样本中SO4 2-的浓度是否与数值模拟计算得到的计算浓度接近；
步骤四：在不同的串层污染井采用同样的方法进行数值模拟，收集数值模拟结果和采样结果；
步骤五：对数值模拟的计算结果可信度进行分析，利用数值模拟计算得出SO4 2-浓度减去实际采样地下水中SO4 2-浓度的绝对值，除以实际采样地下水中SO4 2-的浓度，得到的百分比作为可信度评估的依据，所得到的百分比数值大于80％的记为可信样本组，并对可信样本组对应的串层污染井进行有效封堵或者投放药剂进行治理。
CN202211256594.3A 2022-10-14 2022-10-14 一种用于地下水污染的评估***及方法 Active CN115329607B (zh)

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 CN115329607A (zh) 2022-11-11
CN115329607B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116611274B (zh) * 2023-07-21 2023-09-29 中南大学 一种地下水污染运移可视化数值模拟方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109190280A (zh) * 2018-09-18 2019-01-11 东北农业大学 一种基于核极限学习机替代模型的地下水污染源反演识别方法
CN113239598A (zh) * 2021-06-08 2021-08-10 中国环境科学研究院 一种应用数值模拟的地下水污染源空间综合识别方法
WO2021208393A1 (zh) * 2020-04-15 2021-10-21 北京工业大学 一种大气污染物排放清单的反演估算方法
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模型的三维动态地下水污染模拟的安全饮用水区域确定方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107525907B (zh) * 2017-10-16 2019-12-31 中国环境科学研究院 地下水污染监测网多目标优化方法
CN112149871B (zh) * 2020-08-25 2022-09-13 河海大学 一种基于gis空间统计与随机模拟相结合的污染物点源解析方法
CN112697849B (zh) * 2020-11-20 2022-08-23 中科盛联(北京)新材料有限公司 动态污染源定位方法
CN114254572B (zh) * 2021-12-16 2024-01-02 西北工业大学太仓长三角研究院 考虑污染物沉积的航发压气机流场性能预测方法及***

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109190280A (zh) * 2018-09-18 2019-01-11 东北农业大学 一种基于核极限学习机替代模型的地下水污染源反演识别方法
WO2021208393A1 (zh) * 2020-04-15 2021-10-21 北京工业大学 一种大气污染物排放清单的反演估算方法
CN113239598A (zh) * 2021-06-08 2021-08-10 中国环境科学研究院 一种应用数值模拟的地下水污染源空间综合识别方法
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)

* Cited by examiner, † Cited by third party
Title
Source Inversion of a Sudden Pollution Accident Occurred in Plain River Network Based on Pattern Search Method;Gao Qi等;《IEEE Xplore》;20160121;全文 *
三维非稳定流数学模型在模拟废弃煤矿区地下水串层污染中的应用研究;杨之珍;《能源与环境》;20090430(第02期);全文 *
地下水DNAPLs污染源反演识别及不确定性分析研究;侯泽宇;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20181215;第2018年卷(第12期);全文 *
地下水水源保护区边界污染源防控的V-M数值模拟探析;杨佳立;《科技广场》;20171030(第10期);全文 *
基于贝叶斯理论的地下水DNAPLs污染源反演识别研究;闫雪嫚;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20220115;第2022年卷(第1期);全文 *

Also Published As

Publication number Publication date
CN115329607A (zh) 2022-11-11

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
Butters et al. Field scale transport of bromide in an unsaturated soil: 1. Experimental methodology and results
White Karst hydrology: recent developments and open questions
De Louw et al. Upward groundwater flow in boils as the dominant mechanism of salinization in deep polders, the Netherlands
Boggs et al. Field study of dispersion in a heterogeneous aquifer: 1. Overview and site description
Scheiber et al. Hydrochemical apportioning of irrigation groundwater sources in an alluvial aquifer
Perello et al. Excavation of the Aica-Mules pilot tunnel for the Brenner base tunnel: information gained on water inflows in tunnels in granitic massifs
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
Amiri et al. Hydrogeological assessment and estimation of groundwater inflow into the water transmission tunnel to Urmia Lake, Northwestern Iran
Hauwert Groundwater flow and recharge within the Barton Springs segment of the Edwards Aquifer, southern Travis and northern Hays counties, Texas
Gomo A groundwater-surface water interaction study of an alluvial channel aquifer
McKay et al. Lateral chloride migration from a landfill in a fractured clay‐rich glacial deposit
Lasher Application of fluid electrical conductivity logging for fractured rock aquifer characterisation at the University of the Western Cape's Franschhoek and Rawsonville research sites
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
Ellington Quantification of the impact of irrigation on the aquifer underlying the Vaalharts Irrigation Scheme
Lee et al. Investigating the origin and dynamics of salinity in a confined aquifer system in southeast Australia (Western Port Basin)
Schalla et al. Interim characterization report for the 300 Area process trenches
Lamontagne et al. Interconnection of surface and groundwater systems–River losses from losing/disconnected streams. Border Rivers site report. Water for a Healthy Country Flagship: Adelaide. 1CSIRO Land and Water, PMB 2, Glen Osmond SA 5064 2NSW Office of Water, 155–157 Marius Street, PO Box 550, Tamworth NSW 2340 This project is a National Water Commission (NWC) initiative, funded through the Raising National Water Standards Program. Additional funding was provided by NSW Office of Water
Middleton Aquifer–stream connectivity at various scales: Application of sediment–water interface temperature and vulnerability assessments of groundwater dependent streams
MacDonald et al. Groundwater systems
Alberti et al. Characterization of salt-water intrusion in the lower Esino Valley, Italy using a three-dimensional numerical model
Muldoon et al. Tracer study for characterization of groundwater movement and contaminant transport in fractured dolomite

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