CN114626583A - 一种异质性污染场地地下水修复后反弹潜势预测评估方法 - Google Patents
一种异质性污染场地地下水修复后反弹潜势预测评估方法 Download PDFInfo
- Publication number
- CN114626583A CN114626583A CN202210187202.6A CN202210187202A CN114626583A CN 114626583 A CN114626583 A CN 114626583A CN 202210187202 A CN202210187202 A CN 202210187202A CN 114626583 A CN114626583 A CN 114626583A
- Authority
- CN
- China
- Prior art keywords
- rebound
- groundwater
- soil
- area
- concentration
- 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
- 239000003673 groundwater Substances 0.000 title claims abstract description 111
- 238000005067 remediation Methods 0.000 title claims abstract description 25
- 238000011156 evaluation Methods 0.000 title claims abstract description 13
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 148
- 231100000719 pollutant Toxicity 0.000 claims abstract description 148
- 239000002689 soil Substances 0.000 claims abstract description 122
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000035699 permeability Effects 0.000 claims abstract description 21
- 238000011835 investigation Methods 0.000 claims abstract description 11
- 238000012216 screening Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 96
- 238000005070 sampling Methods 0.000 claims description 46
- 238000005527 soil sampling Methods 0.000 claims description 25
- 230000008439 repair process Effects 0.000 claims description 23
- 238000012544 monitoring process Methods 0.000 claims description 12
- 239000011148 porous material Substances 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000011109 contamination Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 4
- 230000002829 reductive effect Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004445 quantitative analysis Methods 0.000 claims description 3
- 238000012800 visualization Methods 0.000 claims description 3
- 239000000356 contaminant Substances 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims 2
- 238000005259 measurement Methods 0.000 abstract description 3
- 230000003068 static effect Effects 0.000 abstract description 2
- MVPPADPHJFYWMZ-UHFFFAOYSA-N chlorobenzene Chemical compound ClC1=CC=CC=C1 MVPPADPHJFYWMZ-UHFFFAOYSA-N 0.000 description 52
- 239000004927 clay Substances 0.000 description 17
- 230000000694 effects Effects 0.000 description 9
- 238000011161 development Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000011897 real-time detection Methods 0.000 description 3
- 238000012502 risk assessment Methods 0.000 description 3
- 238000012954 risk control Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 239000003463 adsorbent Substances 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000003900 soil pollution Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Images
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/24—Earth materials
- G01N33/246—Earth materials for water content
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/80—Data visualisation
-
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Educational Administration (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Game Theory and Decision Science (AREA)
- Crystallography & Structural Chemistry (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Food Science & Technology (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Primary Health Care (AREA)
- Data Mining & Analysis (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明公开了一种异质性污染场地地下水修复后反弹潜势预测评估方法，包括构件模型、识别潜在污染反弹区域、潜在污染反弹区域污染调查、筛选反弹区域、反弹区低渗透土层污染物残留质量计算、计算反弹区域地下水含水层目标污染物残留质量、区域反弹潜势的预测和评估的步骤；本发明基于反弹区域地下水含水层及其临近区域低渗透地层中目标污染物质量守恒原理，创建了静态平衡情景下反弹区域含水层地下水中目标污染物浓度的预测模型；该模型为解析模型，求解较传统的数值模型简单，参数均可通过现场实测获取，易于推广应用，且模型预测结果的不确定性低。
Description
技术领域
本发明涉及属于场地污染调查与风险评估技术领域，具体涉及一种异质性污染场地地下水修复后反弹潜势预测评估方法。
背景技术
根据土壤法的要求，经污染状况调查与风险评估确定的污染风险超过可接受水平的场地，需开展风险管控或治理修复。同时，土壤法也明确要求，土壤污染责任人或土地使用权人应聘请第三方机构对场地风险管控或治理修复的效果开展评估，以确保达到了风险评估报告中确定的管控或修复目标。为此，生态环境部已制定并实施了《污染地块风险管控与土壤修复效果评估技术导则》(HJ25.5-2018)《污染地块地下水修复和风险管控技术导则》(HJ25.6-2019)，对场地修复效果评估工作的开展提出了相应的技术要求。其中，场地地下水修复效果评估需在修复设施停止运行且场地地下水流场基本稳定的前提下，连续开展2年的监测评估，只有在此期限内地下水中目标污染物浓度未出现明显的反弹升高且未超过相应的修复目标值，方可判断场地地下水修复达标。此技术要求的制定背景主要是考虑到场地水文地质条件的异质性，当场地地下水含水层区域存在非连续的粘土等低渗透透镜体时，如果这些透镜体中存在残留污染物，随着时间的推移，残留污染物可解吸扩散进入含水层造成地下水中目标污染物浓度反弹升高，甚至超过前期风险评估报告中确定的修复目标值。这一技术要求的实施，将场地修复效果评估的时限显著延长，在很大程度上阻碍了场地再开发建设的时序，对地方经济发展带来了不利影响。
目前一致认为，造成场地地下水修复设施停止运行后目标污染物浓度反弹升高的主要驱动机制是含水层低渗透土层中残留污染物缓慢解吸，在浓度梯度的驱动下扩散传输至含水层，导致地下水中目标污染物浓度逐渐升高反弹。这一过程受场地水文地质条件、残余污染物浓度与质量、目标污染物在低渗透地层以及含水层中的迁移转化规律等多重因素的影响，即使按HJ25.6中的技术要求在修复设施停止运行且场地地下水流场恢复稳定的情境下连续监测2年，依然可能捕捉不到反弹升高后地下水中目标污染物最高浓度，导致对反弹后是否超标做出假阴性的判断，造成场地再开发过程中可能暴露环境风险。
现在没有手段能够解决仅依靠两年地下水连续监测判断是否存在目标污染物升高反弹导致结论存在假阴性的技术问题，本发明提供一种异质性污染场地地下水修复后反弹潜势预测评估方法解决上述问题。
发明内容
本发明提供一种异质性污染场地地下水修复后反弹潜势预测评估方法，解决场地水文地质条件、残余污染物浓度与质量、目标污染物在低渗透地层以及含水层中的迁移转化规律等多重因素耦合影响下仅依靠两年地下水连续监测判断是否存在目标污染物升高反弹导致结论存在假阴性的技术问题。
本发明解决上述技术问题所采用的技术方案是：
一种异质性污染场地地下水修复后反弹潜势预测评估方法，包括以下步骤：
S1，构建场地污染概念模型：收集场地污染物浓度与空间分布数据、地层结构及其渗透性数据和地下水数据，利用三维可视化技术将上述数据进行叠加，构建场地污染概念模型，根据场地污染概念模型推断修复工程实施后残余污染物在场地中潜在的空间分布及其赋存介质；
S2，识别潜在污染反弹区域：根据场地污染概念模型，结合污染物的空间分布数据，识别出地下水含水层中低渗透土层，将其划定为会导致地下水污染物潜在反弹区域；
S3，潜在污染反弹区域污染调查：在潜在污染反弹区域进行土壤采样；
S4，筛选反弹区域：针对S3中在潜在污染反弹区域中采集的土壤样品进行污染物浓度检测及分析，
采用公式(1)计算目标污染物在土壤样品孔隙水中的浓度，并将其与潜在反弹区域地下水中对应目标污染物实测浓度进行比较：
式中，Clw为低渗透层土壤样品孔隙水中目标污染物浓度，mg/L；Cls为低渗透层土壤样品中目标污染物中浓度，mg/kg；ρls为低渗透土壤干密度，kg/L；θlw为土壤样品中水分体积含量，无量纲；θla为土壤样品中空气体积含量，无量纲；H为亨利常数，无量纲；kld为目标污染物平衡状态下土壤-水分配系数，L/kg；
根据浓度对比结果，确定污染反弹区域；
S5，反弹区低渗透土层污染物残留质量计算：针对S4中确定的污染反弹区域，采用泰森多边形插值划定每个采样点代表的反弹区域面积Alj；采用公式(2)计算每个反弹区域低渗透土层中目标污染物的残留质量：
式中，Mlk代表第k个反弹区域低渗透土层中目标污染物残留质量，mg；Clji为第k个反弹区域内水平方向第j个土壤采样点低渗透土层垂向上第i个采样位置土样中目标污染物浓度，mg/kg；Clj(i+1)为第k个反弹区域内水平方向第j个土壤采样点低渗透土层垂向上第i+1个采样位置土样中目标污染物浓度，mg/kg；Alj为第k个反弹区域内水平方向第j个土壤采样点代表的反弹区域低渗透土层面积，m2；Llji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间的距离，m；ρlji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间的土壤干密度，kg/L；n为第k个反弹区域内水平方向第j个采样点垂向上土壤采样位置数，无量纲；m为第k个反弹区域内低渗透土层水平方向上土壤采样点数量，无量纲；
S6，计算反弹区域地下水含水层目标污染物残留质量：针对S4中确定的污染反弹区域，采用公式(3)计算每个反弹区地下水含水层中目标污染物残留质量：
Mhk＝(Chw×Ahw×Lhw×θhw+khd×Chw×ρhs×Ahw×Lhw)×1000 (3)
式中，Mhk为修复设施停止运行后第k个反弹区地下水中目标污染物残留质量，mg；Chw为修复设施停止运行后第k个反弹区地下水中目标污染物平均浓度，mg/L；Ahw为每个反弹区域面积，m2；Lhw为每个反弹区域地下水含水层厚度，m；θhw为地下水含水层中土壤样品水分体积含量，无量纲；ρhs为地下水含水层土壤干密度，kg/L；khd为地下水含水层中目标污染物平衡状态下土壤-水分配系数，L/kg；其余符号标记同公式(1)(2)；
S7，预测区域反弹潜势：根据质量守恒原理，采用公式(4)预测平衡状况下各反弹区含水层地下中目标污染物的浓度：
式中，Ckwa为平衡状态下第k个反弹区域含水层地下水中目标污染物浓度，mg/L；θlji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间土壤样品水分体积含量，无量纲；Lhwm为反弹区域含水层地下水混合带高度，m，通常取值不大于2m；θhwm为反弹区域含水层地下水混合带土壤样品水分体积含量，无量纲；ρhsm为反弹区域含水层地下水混合带土壤干密度，kg/L；khdm为反弹区域含水层地下水混合带目标污染物平衡状态下土壤-水分配系数，L/kg；其余符号标记同公式(1)(2)(3)；
S8，评估区域反弹潜势：根据S7的计算结果将Ckwa与修复后场地地下水中目标污染物浓度Cw及其修复目标值进行比较，如Ckwa大于Cw且大于修复目标值，则表明未来该区域地下水中目标污染物浓度会反弹升高，且造成反弹后目标污染物浓度超标的风险较高；如Ckwa仅大于Cw，则表明未来该区域地下水中目标污染物浓度会反弹升高，但造成反弹后目标污染物浓度超标的风险较低；如Ckwa小于Cw，则表明未来该区域地下水中目标污染物浓度会反弹升高的风险较低。
进一步地，步骤S1中，场地污染物浓度与空间分布数据包括场地污染状况调查、修复设计、修复工程实施以及修复设施运行监测阶段中获得的污染浓度与空间分布数据，地下水数据包括地下水含水层分布及其流速与厚度数据。
进一步地，步骤S3中，土壤采样采用分步分阶段进行，具体方法如下：
先进行水平间距大于10m的低密度采样，对采集的土壤样品进行污染物实时检测，然后逐步缩小采样点之间的水平间距至不大于10m的高密度采样；
根据污染物实时检测结果对潜在反弹区域进行水平间距不大于10m的高密度采样，并确定出潜在反弹区域；
针对潜在反弹区域进行垂向间距不大于1m的采样。
进一步地，步骤S4中，当计算的孔隙水中目标污染物浓度低于潜在反弹区域地下水中对应目标污染物的实测浓度时，认为该样品所在局部低渗透土层中残留的污染物不会造成地下水中对应目标污染物浓度反弹；反之，认为样品所在局部低渗透土层中残留的污染物存在造成地下水中对应目标污染物浓度升高反弹的风险，根据后续的定量分析进行确定。
进一步地，步骤S4中，根据计算结果在S1中的场地污染概念模型中进一步进行反弹区域的细化筛选，确定反弹区域的具体分布、尺寸和土层的理化性质，以及每个反弹区域地下水含水层的厚度、含水层赋水介质的孔隙度和地下水流速。
进一步地，步骤S4的计算过程中，模型参数采用场地实测值而非经验参数。
本发明有益效果如下：
针对特定的污染场地，本发明在收集分析场地污染状况调查、修复设计、修复实施与修复设施运行的不同阶段的污染浓度与空间分布信息、地层结构与水文地质条件等信息的基础上，构建场地概念模型，识别潜在反弹区域，并通过污染精细刻画筛选反弹区域；在此基础上，基于质量守恒原理构建平衡状态下反弹区地下水中目标污染物浓度预测模型，据此评估污染物反弹潜势，结果足够保守，且不确定性较低，无需开展地下水中目标污染物浓度变化趋势连续监测，较传统方法有效的避免了地下水修复效果评估周期长、影响场地再开发利用进程、评估结论不确定性高导致再利用过程存在环境安全隐患的缺陷；
本发明中的反弹潜势预测方法基于反弹区域地下水含水层及其临近区域低渗透地层中目标污染物质量守恒原理，创建了静态平衡情景下反弹区域含水层地下水中目标污染物浓度的预测模型；该模型为解析模型，求解较传统的数值模型简单，参数均可通过现场实测获取，易于推广应用，且模型预测结果的不确定性低。
附图说明
图1为本发明的流程步骤示意图；
图2为本发明的场地概念模型示意图；
图3为本发明的反弹区域每个采样点代表的反弹面积划定示意图。
具体实施方式
下面将结合说明书附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本专利的描述中，需要理解的是，术语“中心”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本专利和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本专利的限制。
如图1所示，一种异质性污染场地地下水修复后反弹潜势预测评估方法，包括以下步骤：
S1，构建场地污染概念模型：收集场地污染物浓度与空间分布数据、地层结构及其渗透性数据和地下水数据，利用三维可视化技术将上述数据进行叠加，构建场地污染概念模型，根据场地污染概念模型推断修复工程实施后残余污染物在场地中潜在的空间分布及其赋存介质；
S2，识别潜在污染反弹区域：根据场地污染概念模型，结合污染物的空间分布数据，识别出地下水含水层中可能存在残留污染物的粉质粘土、粘土低渗透土层，将其划定为会导致地下水污染物潜在反弹区域；
S3，潜在污染反弹区域污染调查：在潜在污染反弹区域进行土壤采样；
S4，筛选反弹区域：针对S3中在潜在污染反弹区域中采集的土壤样品进行污染物浓度检测及分析，
采用公式(1)计算目标污染物在土壤样品孔隙水中的浓度，并将其与潜在反弹区域地下水中对应目标污染物实测浓度进行比较：
式中，Clw为低渗透层土壤样品孔隙水中目标污染物浓度，mg/L；Cls为低渗透层土壤样品中目标污染物中浓度，mg/kg；ρls为低渗透土壤干密度，kg/L；θlw为土壤样品中水分体积含量，无量纲；θla为土壤样品中空气体积含量，无量纲；H为亨利常数，无量纲；kld为目标污染物平衡状态下土壤-水分配系数，L/kg；
根据浓度对比结果，确定污染反弹区域；
S5，反弹区低渗透土层污染物残留质量计算：针对S4中确定的污染反弹区域，采用泰森多边形插值划定每个采样点代表的反弹区域面积Alj；采用公式(2)计算每个反弹区域低渗透土层中目标污染物的残留质量：
式中，Mlk代表第k个反弹区域低渗透土层中目标污染物残留质量，mg；Clji为第k个反弹区域内水平方向第j个土壤采样点低渗透土层垂向上第i个采样位置土样中目标污染物浓度，mg/kg；Clj(i+1)为第k个反弹区域内水平方向第j个土壤采样点低渗透土层垂向上第i+1个采样位置土样中目标污染物浓度，mg/kg；Alj为第k个反弹区域内水平方向第j个土壤采样点代表的反弹区域低渗透土层面积，m2；Llji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间的距离，m；ρlji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间的土壤干密度，kg/L；n为第k个反弹区域内水平方向第j个采样点垂向上土壤采样位置数，无量纲；m为第k个反弹区域内低渗透土层水平方向上土壤采样点数量，无量纲；
S6，计算反弹区域地下水含水层目标污染物残留质量：针对S4中确定的污染反弹区域，采用公式(3)计算每个反弹区地下水含水层中目标污染物残留质量：
Mhk＝(Chw×Ahw×Lhw×θhw+khd×Chw×ρhs×Ahw×Lhw)×1000 (3)
式中，Mhk为修复设施停止运行后第k个反弹区地下水中目标污染物残留质量，mg；Chw为修复设施停止运行后第k个反弹区地下水中目标污染物平均浓度，mg/L；Ahw为每个反弹区域面积，m2；Lhw为每个反弹区域地下水含水层厚度，m；θhw为地下水含水层中土壤样品水分体积含量，无量纲；ρhs为地下水含水层土壤干密度，kg/L；khd为地下水含水层中目标污染物平衡状态下土壤-水分配系数，L/kg；其余符号标记同公式(1)(2)；
S7，预测区域反弹潜势：根据质量守恒原理，采用公式(4)预测平衡状况下各反弹区含水层地下中目标污染物的浓度：
式中，Ckwa为平衡状态下第k个反弹区域含水层地下水中目标污染物浓度，mg/L；θlji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间土壤样品水分体积含量，无量纲；Lhwm为反弹区域含水层地下水混合带高度，m，通常取值不大于2m；θhwm为反弹区域含水层地下水混合带土壤样品水分体积含量，无量纲；ρhsm为反弹区域含水层地下水混合带土壤干密度，kg/L；khdm为反弹区域含水层地下水混合带目标污染物平衡状态下土壤-水分配系数，L/kg；其余符号标记同公式(1)(2)(3)；
S8，评估区域反弹潜势：根据S7的计算结果将Ckwa与修复后场地地下水中目标污染物浓度Cw及其修复目标值进行比较，如Ckwa大于Cw且大于修复目标值，则表明未来该区域地下水中目标污染物浓度会反弹升高，且造成反弹后目标污染物浓度超标的风险较高，即有较大概率导致目标污染物浓度超标，需要进一步实施修复或风险管控措施；如Ckwa仅大于Cw，则表明未来该区域地下水中目标污染物浓度会反弹升高，但造成反弹后目标污染物浓度超标的风险较低；如Ckwa小于Cw，则表明未来该区域地下水中目标污染物浓度会反弹升高的风险较低，低渗透土层将作为“吸附剂”进一步吸附含水层中流动相的目标污染物，降低目标污染物的迁移扩散风险。
进一步地，步骤S1中，场地污染物浓度与空间分布数据包括场地污染状况调查、修复设计、修复工程实施以及修复设施运行监测阶段中获得的污染浓度与空间分布数据，地下水数据包括地下水含水层分布及其流速与厚度数据。
进一步地，步骤S3中，土壤采样采用分步分阶段进行，具体方法如下：
先进行水平间距大于10m的低密度采样，对采集的土壤样品进行污染物实时检测，然后逐步缩小采样点之间的水平间距至不大于10m的高密度采样；
根据污染物实时检测结果对潜在反弹区域进行水平间距不大于10m的高密度采样，并确定出潜在反弹区域；
针对潜在反弹区域进行垂向间距不大于1m的采样。
优选的，采用膜界面半定量快速探测方式对采集的土壤样品进行污染物实时检测。
进一步地，步骤S4中，当计算的孔隙水中目标污染物浓度低于潜在反弹区域地下水中对应目标污染物的实测浓度时，认为该样品所在局部低渗透土层中残留的污染物不会造成地下水中对应目标污染物浓度反弹；反之，认为样品所在局部低渗透土层中残留的污染物存在造成地下水中对应目标污染物浓度升高反弹的风险，根据后续的定量分析进行确定。
进一步地，步骤S4中，根据计算结果在S1中的场地污染概念模型中进一步进行反弹区域的细化筛选，确定反弹区域的具体分布、尺寸和土层的理化性质，以及每个反弹区域地下水含水层的厚度、含水层赋水介质的孔隙度和地下水流速。
优选的，反弹区域的具体分布为水平和垂向位置数据，尺寸为面积和厚度数据，土层的理化性质为含水率、孔隙度和干密度数据。
进一步地，步骤S4的计算过程中，模型参数采用场地实测值而非经验参数以提高预测精度。
优选的，步骤S5中采用基于面积与厚度加权的计算方法，预测每个反弹区域低渗透土层中污染物残留质量，有效克服了污染物浓度与水文地质条件各向异性对残留质量预测精度的影响。
本发明的具体实施例，某氯苯污染场地地下水修复效果后反弹潜能预测与评估。场地土壤和地下水主要受氯苯污染，整个场地面积约18000m2，采用原位技术进行修复。修复设施运行约两年后实施单位运行监测结果显示地下水中氯苯浓度已低于修复目标值2.5mg/L。采用本发明的方法对上述场地地下水修复效果开展评估。
步骤一，收集数据构建场地概念模型；如图2所示，实施例场地地面以下18m深度范围内的土层由上而下依次是深度为0-3.5m的填土层，深度为3.5-10m的粉质粘土层，深度为10-12m的粘质粉土层，深度为12-18m的粉细砂层和深度为18-20m的粉粘层，其中，地下水主要赋存在粉细砂层。将污染物浓度及其空间分布数据与场地地层及水文地质条件数据进行三维叠加后发现，修复工程实施前，目标污染物主要集中在填土层之下的粉质粘土层以及地下水含水层上界面的粘质粉土层和下界面的粉质粘土层。修复设施运行过程中监测数据显示，地下水含水层上界面的粘质粉土层和下界面的粉质粘土层中氯苯的检出率与检出浓度较其它区域土层高。从平面分布看，主要集中在场地东南角区域。因此初步推测存在残留污染的区域主要是场地东南侧，污染物主要残留在这一区域地下水含水层上、下界面的低渗透土层。
步骤二，依据步骤一中构建的场地概念模型；初步识别案例场地地下水潜在的反弹区域主要位于场地东南角，其反弹的主要驱动力是这一区域地下水上覆粘质粉土层以及下伏粉质粘土层中残留的氯苯解吸扩散进入地下水含水层，造成这一局部区域地下水中氯苯浓度反弹升高。
步骤三，针对步骤二中识别出的场地地下水潜在反弹区域开展残留污染精细调查；通过高密度布点采样，并配合PID便携式设备快速检测技术，在潜在反弹区域布置土壤采样点31个采集土壤样品总计240个，然后进行氯苯浓度检测，同时设置17口地下水监测井进行氯苯浓度检测。
步骤四，采用公式(1)计算这240个土壤样品对应孔隙水中氯苯浓度，并将其与此区域17个地下水监测井地下水中氯苯按监测井代表的面积进行加权的平均浓度360μg/L进行对比，进一步筛选案例场地反弹区域。结果显示，仅11个点位地下水含水层下伏的粉质粘土层中氯苯孔隙水中浓度大于360μg/L，存在反弹的可能。具体点位信息如表1所示，这些采样点主要位于场地东南角区域，总面积约1653m2，可能造成地下水氯苯浓度反弹的下伏粉质粘土低渗透污染土层厚度约1～1.5m；反弹区域地下水含水层厚度6米，赋水介质为粉砂，因实施例的场地周围已建设止水帷幕，可认为场地内地下水流速为0。
表1土壤孔隙水中氯苯浓度预测结果
步骤五，如图3所示，针对步骤四筛选确定的反弹区域，采用多边形插值划定反弹区域每个采样点代表的反弹区域面积。
采用公式(2)计算反弹区域低渗透土层中氯苯的残留质量，如表2所示。
表2反弹区域低渗透土层氯苯残留质量
步骤六，采用公式(3)计算反弹区地下水含水层中氯苯残留质量，其中，反弹区域地下水中氯苯浓度(Chw)取值为此区域氯苯按面积加权的平均值0.36mg/L，反弹区域的面积为(Ahw)1653m2，每个反弹区域地下水含水层厚度(Lhw)为6m，地下水含水层中土壤样品水分体积含量(θhw)为0.4，地下水含水层土壤干密度(ρhs)为1.7kg/L，地下水含水层中目标污染物平衡状态下土壤-水分配系数(khd)2.05L/kg。将以上参数带入公式(3)计算的含水层中氯苯残留质量为13871314.8mg。
步骤七，采用公式(4)预测平衡状况下各反弹区域含水层中氯苯浓度，结果为2.03mg/L。对比可知，因局部区域地下水含水层下伏低渗透粉质粘土层中残留有氯苯，案例场地地下水修复设施停止后赋存在这些低渗透粉质粘土层中的氯苯将缓慢解吸，在浓度梯度的驱动下扩散传输进入地下水含水层，可能导致地下水中氯苯浓度反弹升高至2.03mg/L，但依然低于该场地地下水氯苯的修复目标值2.5mg/L。表明未来该区域地下水中氯苯浓度回反弹升高的可能性较低。
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内，不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
Claims (6)
1.一种异质性污染场地地下水修复后反弹潜势预测评估方法，其特征是，包括以下步骤：
S1，构建场地污染概念模型：收集场地污染物浓度与空间分布数据、地层结构及其渗透性数据和地下水数据，利用三维可视化技术将上述数据进行叠加，构建场地污染概念模型，根据场地污染概念模型推断修复工程实施后残余污染物在场地中潜在的空间分布及其赋存介质；
S2，识别潜在污染反弹区域：根据场地污染概念模型，结合污染物的空间分布数据，识别出地下水含水层中低渗透土层，将其划定为会导致地下水污染物潜在反弹区域；
S3，潜在污染反弹区域污染调查：在潜在污染反弹区域进行土壤采样；
S4，筛选反弹区域：针对S3中在潜在污染反弹区域中采集的土壤样品进行污染物浓度检测及分析，
采用公式(1)计算目标污染物在土壤样品孔隙水中的浓度，并将其与潜在反弹区域地下水中对应目标污染物实测浓度进行比较：
式中，Clw为低渗透层土壤样品孔隙水中目标污染物浓度，mg/L；Cls为低渗透层土壤样品中目标污染物中浓度，mg/kg；ρls为低渗透土壤干密度，kg/L；θlw为土壤样品中水分体积含量，无量纲；θla为土壤样品中空气体积含量，无量纲；H为亨利常数，无量纲；kld为目标污染物平衡状态下土壤-水分配系数，L/kg；根据浓度对比结果，确定污染反弹区域；
S5，反弹区低渗透土层污染物残留质量计算：针对S4中确定的污染反弹区域，采用泰森多边形插值划定每个采样点代表的反弹区域面积Alj；采用公式(2)计算每个反弹区域低渗透土层中目标污染物的残留质量：
式中，Mlk代表第k个反弹区域低渗透土层中目标污染物残留质量，mg；Clji为第k个反弹区域内水平方向第j个土壤采样点低渗透土层垂向上第i个采样位置土样中目标污染物浓度，mg/kg；Clj(i+1)为第k个反弹区域内水平方向第j个土壤采样点低渗透土层垂向上第i+1个采样位置土样中目标污染物浓度，mg/kg；Alj为第k个反弹区域内水平方向第j个土壤采样点代表的反弹区域低渗透土层面积，m2；Llji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间的距离，m；ρlji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间的土壤干密度，kg/L；n为第k个反弹区域内水平方向第j个采样点垂向上土壤采样位置数，无量纲；m为第k个反弹区域内低渗透土层水平方向上土壤采样点数量，无量纲；
S6，计算反弹区域地下水含水层目标污染物残留质量：针对S4中确定的污染反弹区域，采用公式(3)计算每个反弹区地下水含水层中目标污染物残留质量：
Mhk＝(Chw×Ahw×Lhw×θhw+khd×Chw×ρhs×Ahw×Lhw)×1000 (3)
式中，Mhk为修复设施停止运行后第k个反弹区地下水中目标污染物残留质量，mg；Chw为修复设施停止运行后第k个反弹区地下水中目标污染物平均浓度，mg/L；Ahw为每个反弹区域面积，m2；Lhw为每个反弹区域地下水含水层厚度，m；θhw为地下水含水层中土壤样品水分体积含量，无量纲；ρhs为地下水含水层土壤干密度，kg/L；khd为地下水含水层中目标污染物平衡状态下土壤-水分配系数，L/kg；其余符号标记同公式(1)(2)；
S7，预测区域反弹潜势：根据质量守恒原理，采用公式(4)预测平衡状况下各反弹区含水层地下中目标污染物的浓度：
式中，Ckwa为平衡状态下第k个反弹区域含水层地下水中目标污染物浓度，mg/L；θlji为第k个反弹区域内水平方向第j个采样点低渗透土层垂向上第i个和第i+1个土壤采样位置之间土壤样品水分体积含量，无量纲；Lhwm为反弹区域含水层地下水混合带高度，m，通常取值不大于2m；θhwm为反弹区域含水层地下水混合带土壤样品水分体积含量，无量纲；ρhsm为反弹区域含水层地下水混合带土壤干密度，kg/L；khdm为反弹区域含水层地下水混合带目标污染物平衡状态下土壤-水分配系数，L/kg；其余符号标记同公式(1)(2)(3)；
S8，评估区域反弹潜势：根据S7的计算结果将Ckwa与修复后场地地下水中目标污染物浓度Cw及其修复目标值进行比较，如Ckwa大于Cw且大于修复目标值，则表明未来该区域地下水中目标污染物浓度会反弹升高，且造成反弹后目标污染物浓度超标的风险较高；如Ckwa仅大于Cw，则表明未来该区域地下水中目标污染物浓度会反弹升高，但反弹后造成目标污染物浓度超标的风险较低；如Ckwa小于Cw，则表明未来该区域地下水中目标污染物浓度会反弹升高的风险较低。
2.根据权利要求1所述的一种异质性污染场地地下水修复后反弹潜势预测评估方法，其特征是：步骤S1中，场地污染物浓度与空间分布数据包括场地污染状况调查、修复设计、修复工程实施以及修复设施运行监测阶段获得的污染浓度与空间分布数据，地下水数据包括地下水含水层分布及其流速与厚度数据。
3.根据权利要求1所述的一种异质性污染场地地下水修复后反弹潜势预测评估方法，其特征是：步骤S3中，土壤采样采用分步分阶段进行，具体方法如下：先进行水平间距大于10m的低密度采样，对采集的土壤样品进行污染物实时测试，然后逐步缩小采样点之间的水平间距至不大于10m的高密度采样；
根据污染物实时测试结果对潜在反弹区域进行水平间距不大于10m的高密度采样，并确定出潜在反弹区域；
针对潜在反弹区域进行垂向间距不大于1m的采样。
4.根据权利要求1所述的一种异质性污染场地地下水修复后反弹潜势预测评估方法，其特征是：步骤S4中，当计算的孔隙水中目标污染物浓度低于潜在反弹区域地下水中对应目标污染物的实测浓度时，认为该样品所在局部低渗透土层中残留的污染物不会造成地下水中对应目标污染物浓度反弹；反之，认为样品所在局部低渗透土层中残留的污染物存在造成地下水中对应目标污染物浓度升高反弹的风险，根据后续的定量分析进行确定。
5.根据权利要求4所述的一种异质性污染场地地下水修复后反弹潜势预测评估方法，其特征是：步骤S4中，根据计算结果在S1中的场地污染概念模型中进一步进行反弹区域的细化筛选，确定反弹区域的具体分布、尺寸和土层的理化性质，以及每个反弹区域地下水含水层的厚度、含水层赋水介质的孔隙度和地下水流速。
6.根据权利要求1所述的一种异质性污染场地地下水修复后反弹潜势预测评估方法，其特征是：步骤S4的计算过程中，模型参数采用场地实测值而非经验参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210187202.6A CN114626583B (zh) | 2022-02-28 | 2022-02-28 | 一种异质性污染场地地下水修复后反弹潜势预测评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210187202.6A CN114626583B (zh) | 2022-02-28 | 2022-02-28 | 一种异质性污染场地地下水修复后反弹潜势预测评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114626583A true CN114626583A (zh) | 2022-06-14 |
| CN114626583B CN114626583B (zh) | 2022-09-06 |
Family
ID=81900052
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210187202.6A Active CN114626583B (zh) | 2022-02-28 | 2022-02-28 | 一种异质性污染场地地下水修复后反弹潜势预测评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114626583B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115639340A (zh) * | 2022-09-09 | 2023-01-24 | 生态环境部南京环境科学研究所 | 一种用于农药污染场地的异味类物质筛查方法 |
| CN116312824A (zh) * | 2023-03-20 | 2023-06-23 | 北京市生态环境保护科学研究院 | 针对地下水污染物非菲克弥散的含水层异质性识别方法 |
| CN117274014A (zh) * | 2023-09-26 | 2023-12-22 | 生态环境部环境规划院 | 一种多相抽提技术运行数据处理及管理系统及方法 |
| CN117973839A (zh) * | 2023-12-21 | 2024-05-03 | 北京市生态环境保护科学研究院 | 一种污染场地风险管控工艺参数确定及监测方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106485003A (zh) * | 2016-10-13 | 2017-03-08 | 北京市环境保护科学研究院 | 一种基于特定污染场地的地下水污染层次化评估方法 |
| CN110355193A (zh) * | 2019-07-19 | 2019-10-22 | 中国科学院南京土壤研究所 | 一种基于动态地下水循环的污染场地原位修复方法 |
| CN110449457A (zh) * | 2019-08-23 | 2019-11-15 | 中国科学院地理科学与资源研究所 | 污染土壤测土配方精准靶向修复方法及系统 |
| CN111400660A (zh) * | 2020-03-05 | 2020-07-10 | 北京市环境保护科学研究院 | 一种基于试算法的建设用地土壤修复目标值的计算方法 |
| US20200252283A1 (en) * | 2017-10-16 | 2020-08-06 | Chinese Research Academy Of Environmental Sciences | Multi-objective optimization method for groundwater pollution monitoring network |
| CN113379261A (zh) * | 2021-06-17 | 2021-09-10 | 上海市环境科学研究院 | 一种污染场地土壤地下水污染风险管控技术方法 |
| CN113772804A (zh) * | 2021-09-15 | 2021-12-10 | 宝航环境修复有限公司 | 地下水污染监测自然衰减修复预测方法、系统及装置 |
| CN113820466A (zh) * | 2021-10-09 | 2021-12-21 | 生态环境部南京环境科学研究所 | 一种基于工程安全性评估的修复后污染场地评估方法 |
-
2022
- 2022-02-28 CN CN202210187202.6A patent/CN114626583B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106485003A (zh) * | 2016-10-13 | 2017-03-08 | 北京市环境保护科学研究院 | 一种基于特定污染场地的地下水污染层次化评估方法 |
| US20200252283A1 (en) * | 2017-10-16 | 2020-08-06 | Chinese Research Academy Of Environmental Sciences | Multi-objective optimization method for groundwater pollution monitoring network |
| CN110355193A (zh) * | 2019-07-19 | 2019-10-22 | 中国科学院南京土壤研究所 | 一种基于动态地下水循环的污染场地原位修复方法 |
| CN110449457A (zh) * | 2019-08-23 | 2019-11-15 | 中国科学院地理科学与资源研究所 | 污染土壤测土配方精准靶向修复方法及系统 |
| CN111400660A (zh) * | 2020-03-05 | 2020-07-10 | 北京市环境保护科学研究院 | 一种基于试算法的建设用地土壤修复目标值的计算方法 |
| CN113379261A (zh) * | 2021-06-17 | 2021-09-10 | 上海市环境科学研究院 | 一种污染场地土壤地下水污染风险管控技术方法 |
| CN113772804A (zh) * | 2021-09-15 | 2021-12-10 | 宝航环境修复有限公司 | 地下水污染监测自然衰减修复预测方法、系统及装置 |
| CN113820466A (zh) * | 2021-10-09 | 2021-12-21 | 生态环境部南京环境科学研究所 | 一种基于工程安全性评估的修复后污染场地评估方法 |
Non-Patent Citations (1)
| Title |
|---|
| 钟茂生等: "修复达标土壤回填对地下水环境影响的层次化评估方法应用研究", 《环境科学》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115639340A (zh) * | 2022-09-09 | 2023-01-24 | 生态环境部南京环境科学研究所 | 一种用于农药污染场地的异味类物质筛查方法 |
| CN116312824A (zh) * | 2023-03-20 | 2023-06-23 | 北京市生态环境保护科学研究院 | 针对地下水污染物非菲克弥散的含水层异质性识别方法 |
| CN116312824B (zh) * | 2023-03-20 | 2023-08-22 | 北京市生态环境保护科学研究院 | 针对地下水污染物非菲克弥散的含水层异质性识别方法 |
| CN117274014A (zh) * | 2023-09-26 | 2023-12-22 | 生态环境部环境规划院 | 一种多相抽提技术运行数据处理及管理系统及方法 |
| CN117973839A (zh) * | 2023-12-21 | 2024-05-03 | 北京市生态环境保护科学研究院 | 一种污染场地风险管控工艺参数确定及监测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114626583B (zh) | 2022-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114626583B (zh) | 一种异质性污染场地地下水修复后反弹潜势预测评估方法 | |
| CN101556269B (zh) | 地下水污染模拟槽 | |
| KR101604772B1 (ko) | 다목적 개수로형 lid 실험장치 | |
| Charbeneau | Groundwater hydraulics and pollutant transport | |
| CN107831286B (zh) | 一种地下水污染快速诊断方法 | |
| CN104794298A (zh) | 一种铀尾矿库核素对库区地下水污染风险评价的方法 | |
| CN105678071A (zh) | 一种基于保护地下水的土壤修复目标的层次化制订方法 | |
| Beerten et al. | Geological, hydrogeological and hydrological data for the Dessel disposal site | |
| CN106485003A (zh) | 一种基于特定污染场地的地下水污染层次化评估方法 | |
| CN106517407A (zh) | 一种活性渗滤墙技术修复垃圾填埋场地下水污染的方法 | |
| Huang et al. | Field and numerical modeling study of turbidity current in Shimen Reservoir during typhoon events | |
| CN111008920B (zh) | 一种基于地下水位波动效应的污染场地调查方法 | |
| Follin et al. | Hydrogeological conceptual model development and numerical modelling using CONNECTFLOW, Forsmark modelling stage 2.3 | |
| CN113569330B (zh) | 一种建筑垃圾填埋堆置点评估方法、装置、设备及介质 | |
| Peter et al. | CO2 leakage test in a shallow aquifer for investigating the geochemical impact of CO2 on groundwater and for developing monitoring methods and concepts | |
| Abbaspour et al. | Chemical clogging and geotextile serviceability in subdrains adjacent to recycled concrete | |
| CN114878399A (zh) | 一种可视化土石混合体库岸边坡渗流侵蚀监测系统及方法 | |
| Divine et al. | Polyfluoroalkyl substances requiring a renewed focus on groundwater‐surface water interactions. | |
| Vincenzi et al. | Parametric and numerical modeling tools to forecast hydrogeological impacts of a tunnel | |
| AB | Äspö Hard Rock Laboratory Annual Report 2017 | |
| Gravesen et al. | Geological characterisation of potential disposal areas for radioactive waste from Risø, Denmark | |
| CN117324362A (zh) | 一种针对有色冶炼场地的地下水原位修复方法 | |
| Fordyce et al. | Assessing threats to shallow groundwater quality from soil pollutants in Glasgow, UK: development of a new screening tool | |
| Suuberg et al. | Modeling vapor intrusion processes and evaluating risks using subslab data | |
| CN208334365U (zh) | 一种用于模拟土体中溶质运移的试验装置 |
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 |