一种基于单井抽注的地下水污染物治理方法和装置

技术领域

本发明涉及环境治理领域，具体涉及一种基于单井抽注的地下水污染物治理方法和装置。

背景技术

目前，对地下水的利用过程中广泛存在单井抽注过程，例如，工业废水等污染物通过深井进入承压水层（含水层），并随着地下水的运动而扩散，地下水污染物治理基于单井抽注的地下水污染物治理过程又可以采用抽水的方式进行物理洗涤。然而，地下介质的非均匀结构使得菲克定律难以准确地描述地下水污染物即地下水溶质的运移规律，现有的手段难以了解地下水污染物扩散后的情况，无法准确地制定治理阶段抽水的各项参数。因此，如何有效地制定地下水污染物相应的治理方案是目前业界亟待解决的重要课题。

发明内容

有鉴于此，本发明实施例提供了一种基于单井抽注的地下水污染物治理方法和装置，以解决目前难以了解地下水污染物扩散后的情况导致无法准确地指定治理阶段抽水的各项参数的问题。

根据第一方面，本发明实施例提供了一种基于单井抽注的地下水污染物治理方法，所述方法包括：

确定待治理区域的第一参数、第二参数和第三参数；所述第一参数为表征待治理区域的地下水污染物在注水阶段的参数，所述第二参数为表征所述地下水污染物在静置阶段的参数；所述第三参数为表征所述待治理区域含水层的水力参数；

将所述第一参数、所述第二参数和所述第三参数输入至抽注三阶段溶质运移模型中，得到所述抽注三阶段溶质运移模型输出的第四参数；所述抽注三阶段溶质运移模型是基于样本区域含水层的样本水质参数以及所述样本区域的地下水污染物在所述注水阶段、所述静置阶段和抽水阶段采集的第一样本抽注参数训练得到的，所述注水阶段、所述静置阶段和所述抽水阶段依次顺序进行；所述第四参数为表征所述地下水污染物在所述抽水阶段的参数，所述第四参数包括所述抽水阶段的开始时刻、持续时长和抽水速度；

基于所述第四参数，制定待治理区域的治理方案。

结合第一方面，在第一方面第一实施方式中，所述确定待治理区域的第一参数、第二参数和第三参数，包括：

基于所述待治理区域的勘察数据，确定所述第三参数；所述第三参数包括反映非均匀性和各向异性的水质参数；

确定所述第一参数和所述第二参数；所述第一参数包括所述注水阶段的开始时刻、持续时长、注水速度和注水体积；所述第二参数为所述静置阶段的持续时长。

结合第一方面，在第一方面第二实施方式中，所述抽注三阶段溶质运移模型通过以下步骤建立得到：

基于所述第三参数，确定所述样本区域；

确定所述样本区域的所述样本水质参数；所述样本水质参数为反映非均匀性和各向异性的水质参数；

确定所述第一样本抽注参数；所述第一样本抽注参数包括地下水污染物在所述注水阶段的开始时刻、持续时长、注水速度及注水体积、所述静置阶段的持续时长、所述抽水阶段的开始时刻、持续时长及抽水速度；

将所述样本水质参数和所述第一样本抽注参数作为建立所使用的参数，建立抽注三阶段溶质运移模型；所述抽注三阶段溶质运移模型是基于所述注水阶段、所述静置阶段和所述抽水阶段三个阶段并分阶段建立的极坐标下的分数阶导数方程。

结合第一方面，在第一方面第三实施方式中，所述抽注三阶段溶质运移模型基于样本区域含水层的样本水质参数以及向所述样本区域注水的示踪剂在含水层的所述注水阶段、所述静置阶段和所述抽水阶段采集的第二样本抽注参数训练得到的；

相应的，所述抽注三阶段溶质运移模型通过以下步骤建立得到：

基于所述第三参数，确定所述样本区域；

确定所述样本区域的所述样本水质参数；所述样本水质参数为反映非均匀性和各向异性的水质参数；

向所述样本区域注水示踪剂，并确定所述第二样本抽注参数；所述第二样本抽注参数包括所述示踪剂在所述注水阶段的开始时刻、持续时长、注水速度及注水体积、所述静置阶段的持续时长、所述抽水阶段的开始时刻、持续时长及抽水速度；

将所述样本水质参数和所述第二样本抽注参数作为建立所使用的参数，建立抽注三阶段溶质运移模型；所述抽注三阶段溶质运移模型是基于所述注水阶段、所述静置阶段和所述抽水阶段三个阶段并分阶段建立的极坐标下的分数阶导数方程。

结合第一方面第二实施方式或者第一方面第三实施方式，在第一方面第四实施方式中，所述抽注三阶段溶质运移模型为：

式中，

结合第一方面第二实施方式，在第一方面第五实施方式中，所述抽注三阶段溶质运移模型还基于所述样本区域的工况参数训练得到的；所述工况参数包括地下水污染物的初始浓度和边界条件；

相应的，所述抽注三阶段溶质运移模型通过以下步骤建立得到：

基于所述第三参数，确定所述样本区域；

确定所述样本区域的所述样本水质参数和所述工况参数；

确定所述第一样本抽注参数；所述第一样本抽注参数包括所述注水阶段的开始时刻、持续时长、注水速度及注水体积、所述静置阶段的持续时长、所述抽水阶段的开始时刻、持续时长及抽水速度；

将所述样本水质参数、所述第一样本抽注参数和所述工况参数作为建立所使用的参数，建立抽注三阶段溶质运移模型；所述抽注三阶段溶质运移模型是基于所述注水阶段、所述静置阶段和所述抽水阶段三个阶段并分阶段建立的极坐标下的分数阶导数方程。

结合第一方面第五实施方式，在第一方面第六实施方式中，所述抽注三阶段溶质运移模型为：

式中，

第二方面，本发明实施例还提供一种基于单井抽注的地下水污染物治理装置，其特征在于，所述装置包括：

参数确定模块，用于确定待治理区域的第一参数、第二参数和第三参数；所述第一参数为表征待治理区域的地下水污染物在注水阶段的参数，所述第二参数为表征所述地下水污染物在静置阶段的参数；所述第三参数为表征所述待治理区域含水层的水力参数；

参数输入模块，用于将所述第一参数、所述第二参数和所述第三参数输入至抽注三阶段溶质运移模型中，得到所述抽注三阶段溶质运移模型输出的第四参数；所述抽注三阶段溶质运移模型是基于样本区域含水层的样本水质参数以及所述样本区域的及样本区域的地下水污染物在所述注水阶段、所述静置阶段和抽水阶段采集的第一样本抽注参数训练得到的，所述注水阶段、所述静置阶段和所述抽水阶段依次顺序进行；所述第四参数为表征所述地下水污染物在所述抽水阶段的参数，所述第四参数包括所述抽水阶段的开始时刻、持续时长和抽水速度；

方案确定模块，用于基于所述第四参数，制定待治理区域的治理方案。

第三方面， 本发明实施例还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述基于单井抽注的地下水污染物治理方法的步骤。

第四方面，本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述基于单井抽注的地下水污染物治理方法的步骤。

第五方面，本发明实施例还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述基于单井抽注的地下水污染物治理方法的步骤。

本发明提供的基于单井抽注的地下水污染物治理方法、装置、设备，基于单井抽注的前两个物理阶段以及待治理地区的地质特征，特别是待治理地区含水层的水力参数，水力参数能够一定程度上反映待治理地区地下介质的非均匀性以及各向异性，通过抽注三阶段溶质运移模型模拟地下水污染物即地下水溶质的运移过程，了解地下水污染物在待治理区域的分布情况，以此为依据评估地下水的污染程度，并借由抽注三阶段溶质运移模型预测得到的分布情况来预测地下水污染物在单井抽注的第三个阶段的各项参数，再基于预测得到的各项参数，制定待治理区域进行物理洗涤时的具体参数（指标），使得治理方案更能够适用待治理区域地下介质的物理结构，能够更好地评估和修复地下水污染，助力地下水动力学在生产生活中的应用，为地下水污染评估和修复提供相应的支撑。

附图说明

通过参考附图会更加清楚的理解本发明的特征和优点，附图是示意性的而不应理解为对本发明进行任何限制，在附图中：

图1示出了本发明提供的基于单井抽注的地下水污染物治理方法的流程示意图；

图2示出了本发明提供的基于单井抽注的地下水污染物治理方法中步骤S101具体的流程示意图；

图3示出了本发明提供的基于单井抽注的地下水污染物治理方法中建立抽注三阶段溶质运移模型的流程示意图之一；

图4示出了本发明提供的基于单井抽注的地下水污染物治理方法中建立抽注三阶段溶质运移模型的流程示意图之二；

图5示出了本发明提供的基于单井抽注的地下水污染物治理方法中建立抽注三阶段溶质运移模型的流程示意图之三；

图6示出了本发明提供的基于单井抽注的地下水污染物治理方法中CH3观测井的实验条件示意图；

图7示出了本发明提供的基于单井抽注的地下水污染物治理方法中CH3观测井中部分参数的示意图；

图8示出了本发明提供的基于单井抽注的地下水污染物治理方法中CH10观测井的实验条件示意图；

图9示出了本发明提供的基于单井抽注的地下水污染物治理方法中CH10观测井中部分参数的示意图；

图10示出了本发明提供的基于单井抽注的地下水污染物治理方法中CH3观测井的模型模拟结果示意图；

图11示出了本发明提供的基于单井抽注的地下水污染物治理方法中CH10观测井的模型模拟结果示意图；

图12示出了本发明提供的基于单井抽注的地下水污染物治理装置的结构示意图；

图13示出了本发明提供的电子设备的结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

下面结合图1描述本申请的基于单井抽注的地下水污染物治理方法，本申请是基于单井抽注过程来实现的，该方法包括：

S101、确定待治理区域的第一参数、第二参数和第三参数。

在本实施例中，第一参数为表征待治理区域的地下水污染物在注水阶段的参数，第二参数为表征地下水污染物在静置阶段的参数，第三参数为表征待治理区域含水层的水力参数。

对于单井抽注过程而言，其主要分为三个阶段，注水阶段、静置阶段和抽水阶段，并且，注水阶段、静置阶段和抽水阶段是依次顺序进行的，即，注水阶段之后方可进入静置阶段，静置阶段之后方可进入抽水阶段。

在复杂介质的运移过程中，非均质的含水层的结构决定着流体的流动方向和污染物的分布状况，复杂地质条件下同一地层的渗透系数范围相差可达数个数量级，地下水污染物在渗透系数较高的区域扩散快，而在渗透系数相对较低的区域迁移扩散慢。并且，当地下水污染物在地层分布状况不明时，现有钻探取样调查容易造成检测盲区，难以准确识别场地污染情况，严重阻碍场地的精准修复治理。在本申请中，第三参数包括反映非均匀性和各向异性的水质参数，通过了解待治理地区的非均匀性和各向异性，能够更好地模拟地下水污染物的运移过程和了解相应的运移规律。

优选的，第三参数可以包括渗透系数、贮水系数、导水系数、释水系数、给水度等水力参数。

S102、将第一参数、第二参数和第三参数输入至抽注三阶段溶质运移模型中，得到抽注三阶段溶质运移模型输出的第四参数。

在本申请中，抽注三阶段溶质运移模型是基于样本区域含水层的样本水质参数以及样本区域的地下水污染物在注水阶段、静置阶段和抽水阶段采集的第一样本抽注参数训练得到的。

在本实施例中，第一参数包括注水阶段的开始时刻、持续时长、注水速度和注水体积；第二参数为静置阶段的持续时长；第四参数为表征地下水污染物在抽水阶段的参数，第四参数包括抽水阶段的开始时刻、持续时长和抽水速度。

由于分数阶模型具有记忆性特点，可以刻画污染物运移过程中由于粘土或裂隙死穴带来的滞留现象，目前被广泛应用在地下水溶质运移模拟领域。但是现有的分数阶模型难以刻画地下水污染物在地下介质中注水阶段、静置阶段和抽水阶段三个物理过程的溶质运移规律。

本申请所建立并采用的抽注三阶段溶质运移模型是将注水阶段、静置阶段和抽水阶段这三个物理阶段间隔即分阶段单独描述、刻画和模拟，用以描述注水、静置、抽水三个不同的水力阶段发生的水力梯度大小和方向的瞬态变化，使得模拟得到的污染物溶质的迁移规律更加精确，能够模拟得到更为精确的地下水污染物在待治理区域的分布情况。

更具体的，抽注三阶段溶质运移模型是基于注水阶段、静置阶段和抽水阶段三个阶段并分阶段建立的极坐标下的分数阶导数方程。

S103、基于第四参数，制定待治理区域的治理方案。通过该治理方案进行地下水污染的相应治理。

本发明的基于单井抽注的地下水污染物治理方法，基于单井抽注的前两个物理阶段以及待治理地区的地质特征，特别是待治理地区含水层的水力参数，水力参数能够一定程度上反映待治理地区地下介质的非均匀性以及各向异性，通过抽注三阶段溶质运移模型模拟地下水污染物即地下水溶质的运移过程，了解地下水污染物在待治理区域的分布情况，以此为依据评估地下水的污染程度，并借由抽注三阶段溶质运移模型预测得到的分布情况来预测地下水污染物在单井抽注的第三个阶段的各项参数，再基于预测得到的各项参数，制定待治理区域进行物理洗涤时的具体参数（指标），使得治理方案更能够适用待治理区域地下介质的物理结构，能够更好地评估和修复地下水污染，助力地下水动力学在生产生活中的应用，为地下水污染评估和修复提供相应的支撑。

下面结合图2描述本发明的基于单井抽注的地下水污染物治理方法，步骤S101具体包括：

S1011、基于待治理区域的勘察数据，确定第三参数。在一些可能的实施例中，通过收集已有的研究区域的资料或者开展相关的地质调查例如开展单井抽注实验，测量待治理区域含水层沿不同方向的水力参数，分析该区域的地质结构，主要分析含水层的各向异性特点。

S1012、确定第一参数和第二参数。记录向待治理区域投入待治理区域时的各项指标参数，得到本申请所需的第一参数和第二参数。

下面结合图3描述本发明的基于单井抽注的地下水污染物治理方法，本申请中的抽注三阶段溶质运移模型通过以下步骤建立得到：

A101、基于第三参数，确定样本区域。

本申请在筛选样本区域时，基于已经得到的第三参数以及作为候选的区域已知的地下介质结构的参数，选择与待治理区域地下介质结构相似的候选区域作为样本区域。可以理解的是，可以选择至少一个候选区域作为本申请中的样本区域。

A102、确定样本区域的样本水质参数，同样的，样本水质参数为反映非均匀性和各向异性的水质参数。样本水质参数的数值会与第三参数的数值存在略微的差别，具体会根据所选择的样本区域，在参数的具体数值上存在不同，但都会反映相同的地下介质结构。例如，样本区域的样本水质参数的每项参数均为相应的第三参数±预设值。优选的，预设值可以为3%。

A103、确定第一样本抽注参数。

在本申请中，第一样本抽注参数包括地下水污染物在注水阶段的开始时刻、持续时长、注水速度及注水体积、静置阶段的持续时长、抽水阶段的开始时刻、持续时长及抽水速度。

A104、将样本水质参数和第一样本抽注参数作为建立所使用的参数，建立抽注三阶段溶质运移模型。

考虑到向样本区域投入污染物，会破坏区域的水质，造成对水资源的污染，在另一些可能的实施例中，抽注三阶段溶质运移模型基于样本区域含水层的样本水质参数以及向样本区域注水的示踪剂在含水层的注水阶段、静置阶段和抽水阶段采集的第二样本抽注参数训练得到的。选用对水质几乎无污染的示踪剂来代替污染物，并通过向样本区域投入示踪剂，之后获取示踪剂在含水层的注水阶段、静置阶段和抽水阶段采集的第二样本抽注参数，以此来建立抽注三阶段溶质运移模型。需要说明的是，此时，第二样本抽注参数包括示踪剂在注水阶段的开始时刻、持续时长、注水速度及注水体积、静置阶段的持续时长、抽水阶段的开始时刻、持续时长及抽水速度。

相应的，以示踪剂选用荧光素为例进行说明，下面结合图4描述本发明的基于单井抽注的地下水污染物治理方法，抽注三阶段溶质运移模型通过以下步骤建立得到：

A201、基于第三参数，确定样本区域。

本申请在筛选样本区域时，基于已经得到的第三参数以及作为候选的区域已知的地下介质结构的参数，选择与待治理区域地下介质结构相似的候选区域作为样本区域。可以理解的是，可以选择至少一个候选区域作为本申请中的样本区域。

A202、确定样本区域的样本水质参数，同样的，样本水质参数为反映非均匀性和各向异性的水质参数。

A203、向样本区域注水示踪剂，持续观察示踪剂在注水阶段、静置阶段和抽水阶段的浓度，整理、收集实验数据，并确定第二样本抽注参数。

A204、将样本水质参数和第二样本抽注参数作为建立所使用的参数，建立抽注三阶段溶质运移模型。

在这些实施例中，抽注三阶段溶质运移模型是基于注水阶段、静置阶段和抽水阶段三个阶段并分阶段建立的极坐标下的分数阶导数方程。

在上述实施例中，注水阶段和抽水阶段的水力梯度方向会发生变化，而静置阶段对流作用弱，静置阶段主要考虑扩散过程，这三个阶段都应考虑含水介质对溶质运移的滞留效果，因此，需要建立极坐标下的三阶段溶质运移的三阶段分数阶对流扩散模型也就是脉冲分数阶模型，得到抽注三阶段溶质运移模型。即，针对注水阶段、静置阶段和抽水阶段，采集样本区域的样本水质参数，以及第一样本抽注参数/第二样本抽注参数，并以此为建立模型所用的数据，建立适用于当前溶质运移过程的三阶段分数阶对流扩散模型。

需要说明得是，在某些条件下无静置阶段，即静置阶段的持续时长为0。

在上述实施例中，抽注三阶段溶质运移模型为：

式中，

上式中，分数阶阶数参数均为无量纲量，并适用于反映地下含水层的非均质性也就是地下介质的复杂程度，分数阶阶数低则含水层的滞留效果强，可能为裂隙发育不健全，死穴过多，反之，含水层的滞留效果弱，并且，注水阶段和抽水阶段的分数阶阶数的差异越大则反映了含水层的各向异性越强。

基于已确定的样本水质参数，以及第一样本抽注参数/第二样本抽注参数，率定上述公式中各个参数的具体值或者取值范围。

在另一些可能的实施例中，示踪剂还可以选用待治理区域的地下水中存在的物质，例如氯化钠（NaCl）。区别于上述的实施例，在这些实施例中，还需要考虑实际的工况，包括初始条件和边界条件，也就是抽注三阶段溶质运移模型还基于样本区域的工况参数训练得到的，工况参数包括地下水污染物的初始浓度和边界条件。

相应的，以示踪剂选用NaCl为例进行说明，下面结合图5描述本发明的基于单井抽注的地下水污染物治理方法，抽注三阶段溶质运移模型通过以下步骤建立得到：

A301、基于第三参数，确定样本区域。该步骤与步骤A101以及步骤A201类似，在此不做过多阐述。

A302、确定样本区域的样本水质参数和工况参数。同样可以基于单井抽注过程，确定样本水质参数和工况参数。

A303、确定第一样本抽注参数。在这些实施例中，第一样本抽注参数包括注水阶段的开始时刻、持续时长、注水速度及注水体积、静置阶段的持续时长、抽水阶段的开始时刻、持续时长及抽水速度。

A304、将样本水质参数、第一样本抽注参数和工况参数作为建立所使用的参数，建立抽注三阶段溶质运移模型。

在这些实施例中，抽注三阶段溶质运移模型是基于注水阶段、静置阶段和抽水阶段三个阶段并分阶段建立的极坐标下的分数阶导数方程。

在这些实施例中，抽注三阶段溶质运移模型为：

式中，

一般形式的截断分数阶导数定义为：

Riemann-Liouville定义下的分数阶导数定义为：

其中，

其中，

在本申请中，利用隐式有限差分方法求解抽注三阶段溶质运移模型，具体的：

当

当

当

当

当

当

当

当

当

以待治理区域的地质特征为风化断裂的花岗岩地质为例进行说明，通过本申请的方案，能够制定风化断裂的花岗岩地质区域的地下水污染物的治理方案。

花岗岩地质中结晶岩的裂隙发育很差，很少存在长的连通裂缝，因而，示踪剂溶液经过注水阶段、静置阶段、抽水阶段这三个不同的水力阶段发生的水力梯度大小和方向的瞬态变化，抽水阶段的溶质方向可能会产生超扩散过程。为了对风化断裂的花岗岩地质环境进行保护，可以在该地质环境处先进行地质探孔，并在探孔处开展本发明所需的单井示踪剂实验和实验期间对应的数值模拟与记录。可以理解的是，本申请可以选用多个观测井，具体的，本申请选了CH3和CH10两口观测井进行单井示踪剂实验，且，示踪剂均选用NaCl。

此时CH3和CH10两口观测井的实验条件分别如图6、7和图8、9所示，且收集实验数据后，得到的模型模拟结果图分别如图10和图11所示。

从图10和图11可以看出，本申请抽注三阶段溶质运移模型可以很好的刻画示踪剂粒子峰值的抵达和次扩散逐渐衰减的过程。当注水时间较短又无静置期时，模型的分数阶截断系数较高，表明溶质运移的滞留过程衰减更快；示踪剂粒子注入时间短且无静置期时，部分示踪剂会储存在离井较近的非连通裂隙中，当抽水时，粒子会沿裂隙通道快速回流到测试井，穿透曲线上表现为早峰的现象；示踪剂注入时间较长且存在静置期时，则不会观察到早峰现象，溶质回流过程早到达现象不明显，此时，模型的截断系数较高，这与抽水时，粒子容易找到回流通道有关。因此可以推测该区域的单井抽水阶段更容易洗出污染物溶质。

下面对本发明提供的基于单井抽注的地下水污染物治理装置进行描述，下文描述的基于单井抽注的地下水污染物治理装置与上文描述的基于单井抽注的地下水污染物治理方法可相互对应参照。

下面结合图12描述本申请的基于单井抽注的地下水污染物治理装置，本申请是基于单井抽注过程来实现的，该装置包括：

参数确定模块101，用于确定待治理区域的第一参数、第二参数和第三参数。

在本实施例中，第一参数为表征待治理区域的地下水污染物在注水阶段的参数，第二参数为表征地下水污染物在静置阶段的参数，第三参数为表征待治理区域含水层的水力参数。

对于单井抽注过程而言，其主要分为三个阶段，注水阶段、静置阶段和抽水阶段，并且，注水阶段、静置阶段和抽水阶段是依次顺序进行的，即，注水阶段之后方可进入静置阶段，静置阶段之后方可进入抽水阶段。

参数输入模块102，用于将第一参数、第二参数和第三参数输入至抽注三阶段溶质运移模型中，得到抽注三阶段溶质运移模型输出的第四参数。

方案确定模块103，用于基于第四参数，制定待治理区域的治理方案。通过该治理方案进行地下水污染的相应治理。

本发明的基于单井抽注的地下水污染物治理装置，基于单井抽注的前两个物理阶段以及待治理地区的地质特征，特别是待治理地区含水层的水力参数，水力参数能够一定程度上反映待治理地区地下介质的非均匀性以及各向异性，通过抽注三阶段溶质运移模型模拟地下水污染物即地下水溶质的运移过程，了解地下水污染物在待治理区域的分布情况，以此为依据评估地下水的污染程度，并借由抽注三阶段溶质运移模型预测得到的分布情况来预测地下水污染物在单井抽注的第三个阶段的各项参数，再基于预测得到的各项参数，制定待治理区域进行物理洗涤时的具体参数（指标），使得治理方案更能够适用待治理区域地下介质的物理结构，能够更好地评估和修复地下水污染，助力地下水动力学在生产生活中的应用，为地下水污染评估和修复提供相应的支撑。

图13示例了一种电子设备的实体结构示意图，如图13所示，该电子设备可以包括：处理器(processor)510、通信接口(Communications Interface)520、存储器(memory)530和通信总线540，其中，处理器510，通信接口520，存储器530通过通信总线540完成相互间的通信。处理器510可以调用存储器530中的逻辑命令，以执行基于单井抽注的地下水污染物治理方法，该方法包括：

确定待治理区域的第一参数、第二参数和第三参数；所述第一参数为表征待治理区域的地下水污染物在注水阶段的参数，所述第二参数为表征所述地下水污染物在静置阶段的参数；所述第三参数为表征所述待治理区域含水层的水力参数；

将所述第一参数、所述第二参数和所述第三参数输入至抽注三阶段溶质运移模型中，得到所述抽注三阶段溶质运移模型输出的第四参数；所述抽注三阶段溶质运移模型是基于样本区域含水层的样本水质参数以及所述样本区域的地下水污染物在所述注水阶段、所述静置阶段和抽水阶段采集的第一样本抽注参数训练得到的，所述注水阶段、所述静置阶段和所述抽水阶段依次顺序进行；所述第四参数为表征所述地下水污染物在所述抽水阶段的参数，所述第四参数包括所述抽水阶段的开始时刻、持续时长和抽水速度；

基于所述第四参数，制定待治理区域的治理方案。

此外，上述的存储器530中的逻辑命令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干命令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（ROM，Read-Only Memory）、随机存取存储器（RAM，Random Access Memory）、磁碟或者光盘等各种可以存储程序代码的介质。

另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各方法所提供的基于单井抽注的地下水污染物治理方法，该方法包括：

确定待治理区域的第一参数、第二参数和第三参数；所述第一参数为表征待治理区域的地下水污染物在注水阶段的参数，所述第二参数为表征所述地下水污染物在静置阶段的参数；所述第三参数为表征所述待治理区域含水层的水力参数；

将所述第一参数、所述第二参数和所述第三参数输入至抽注三阶段溶质运移模型中，得到所述抽注三阶段溶质运移模型输出的第四参数；所述抽注三阶段溶质运移模型是基于样本区域含水层的样本水质参数以及所述样本区域的地下水污染物在所述注水阶段、所述静置阶段和抽水阶段采集的第一样本抽注参数训练得到的，所述注水阶段、所述静置阶段和所述抽水阶段依次顺序进行；所述第四参数为表征所述地下水污染物在所述抽水阶段的参数，所述第四参数包括所述抽水阶段的开始时刻、持续时长和抽水速度；

基于所述第四参数，制定待治理区域的治理方案。

又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的基于单井抽注的地下水污染物治理方法，该方法包括：

确定待治理区域的第一参数、第二参数和第三参数；所述第一参数为表征待治理区域的地下水污染物在注水阶段的参数，所述第二参数为表征所述地下水污染物在静置阶段的参数；所述第三参数为表征所述待治理区域含水层的水力参数；

将所述第一参数、所述第二参数和所述第三参数输入至抽注三阶段溶质运移模型中，得到所述抽注三阶段溶质运移模型输出的第四参数；所述抽注三阶段溶质运移模型是基于样本区域含水层的样本水质参数以及所述样本区域的地下水污染物在所述注水阶段、所述静置阶段和抽水阶段采集的第一样本抽注参数训练得到的，所述注水阶段、所述静置阶段和所述抽水阶段依次顺序进行；所述第四参数为表征所述地下水污染物在所述抽水阶段的参数，所述第四参数包括所述抽水阶段的开始时刻、持续时长和抽水速度；

基于所述第四参数，制定待治理区域的治理方案。

以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干命令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行各个实施例或者实施例的某些部分所述的方法。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。