申请/专利权人：河海大学

申请日：2019-03-07

公开（公告）日：2019-07-09

公开（公告）号：CN109992839A

主分类号：G06F17/50(20060101)

分类号：G06F17/50(20060101);G01M10/00(20060101)

优先权：

专利状态码：有效-授权

法律状态：2022.09.16#授权;2019.08.02#实质审查的生效;2019.07.09#公开

摘要：本发明公开了一种避免抽水实验中水位数据噪音给反演计算带来误差的方法，通过基于一维、均质和各向同性的含水层中的地下水流的控制方程，表述关于降深的方程，结合实际的抽水状态，叠加引入水位恢复阶段的反向抽水；通过于公式推演，以降深的一阶导数表述的水力衰减，分析不同状态下水位恢复阶段和抽水阶段的降深Srec的衰减具有相关的一致性；基于上述的一致性，可用恢复阶段的降深数据来代替实际抽水阶段的降深数据进行反演；使用本发明的数据计算水力学走时，相对于使用水位下降数据计算的走时要精确很多，可有效避免抽水实验中的水位数据噪音给反演计算带来误差，具有很强的实用性和广泛的适用性。

主权项：1.避免抽水实验中水位数据噪音给反演计算带来误差的方法，其特征在于，包括以下步骤：S1、基于叠加原理，将抽水稳定后的降深Srec视为抽水和水位恢复的叠加，同时，将水位恢复视为反向抽水，抽水和反向抽水的降深分别为Spump1和Spump2，即Srec＝Spump1+Spump2；S2、基于一维、均质和各向同性的含水层中的地下水流的控制方程，表述关于降深Spump的方程，并以Spump的一阶导数表述水力衰减，二阶导数表述最大水力衰减，同时，获得最大值及最大值条件；S3、通过比较分别当抽水达到稳定或者近似稳定状态、达到准稳定状态时的降深Srec的取值，获得水位恢复阶段和抽水阶段的降深Srec的衰减具有相关的一致性；S4、基于上述的一致性，用水位恢复阶段的降深数据来代替实际抽水阶段的降深数据进行反演，以避免抽水实验中的水位数据噪音给反演计算带来误差。

全文数据：避免抽水实验中水位数据噪音给反演计算带来误差的方法技术领域本发明涉及一种避免水位数据噪音的方法，具体涉及一种避免抽水实验中水位数据噪音给反演计算带来误差的方法，属于水文地质与工程地质参数识别技术领域。背景技术近年来，水力层析法成为了水文地质反演中的一个热点，这种方法是建立在大量的地下水动力学实验和参数反演的基础上的。以抽水实验为例，在实验井中的某深度抽水的同时，在不同位置，不同深度的观测井中记录水头变化；实验结束后，改变实验点和观测点的位置，再次进行实验改变激发源井点和信号接收井点的位置；经多次实验后，就能够得到有一系列测量数据，这些信息犹如在空间组成的信息交叉网络，通过特定的反演算法，可以得到地下水含水层的水动力学参数的具体数值和分布。目前，在抽水试验过程中，由于设备设计和电路等对压力感应器的影响，使得抽水过程中所记录的水位变化的数据噪音较大，而这些噪音造成走时计算的不准确，最终导致反演不准确。因此，寻找一种减少或者避免噪音干扰的方法用于计算水力学走时，尤为重要。发明内容为解决现有技术的不足，本发明的目的在于提供一种避免抽水实验中水位数据噪音给反演计算带来误差的方法。为了实现上述目标，本发明采用如下的技术方案：避免抽水实验中水位数据噪音给反演计算带来误差的方法，包括以下步骤：S1、基于叠加原理，将抽水稳定后的降深Srec视为抽水和水位恢复的叠加，同时，将水位恢复视为反向抽水，抽水和反向抽水的降深分别为Spump1和Spump2，即Srec＝Spump1+Spump2；S2、基于一维、均质和各向同性的含水层中的地下水流的控制方程，表述关于降深Spump的方程，并以Spump的一阶导数表述水力衰减，二阶导数表述最大水力衰减，同时，获得最大值及最大值条件；S3、通过比较分别当抽水达到稳定或者近似稳定状态、达到准稳定状态时的降深Srec的取值，获得水位恢复阶段和抽水阶段的降深Srec的衰减具有相关的一致性；S4、基于上述的一致性，用恢复阶段的降深数据来代替实际抽水阶段的降深数据进行反演，以避免抽水实验中的水位数据噪音给反演计算带来误差。上述步骤S2中表述关于降深Spump的方程，步骤如下：A1、基于一维、均质和各向同性的含水层中的地下水流的控制方程，如下所示：式1中，K为渗透系数，SS为储水率，H为水头，r为任意点距离抽水井中心轴的水平距离，rp为P点距离抽水井中心轴的水平距离，Qrp，t为rp初始单位体积内的抽水率，t为抽水开始后的时间；A2、将方程式1，可以为描述关于降深的方程：式2中，s＝H-H0为降深，H0为初始水头；其中，表示了抽水前的流场散度，当初始水头散度稳定或者准稳定时，其值为零，则，方程式2可以写为：式3中，为水力扩散系数；当抽水率Q为常数，初始降深为零，并且区域无限大时，方程的解为：式中，Spump时抽水期间的降深，t为抽水开始后的时间。上述步骤S2中，以Spump的一阶导数表述水力衰减，二阶导数表述最大水力衰减，则：且，的最大值在为零时得出，为了满足这个条件，令此外，的最大值为：上述步骤S1中，基于叠加原理，将降深Srec视为抽水和反向抽水的叠加，即式8中，Spump1是当抽水率为常数Q时的降深，Spump2是虚拟的反向抽水中当抽水率为常数-Q时的降深；tp为停止抽水时的时间，t’为停止抽水后的时间，即，t＝tp+t’。上述步骤S3中当抽水达到稳定或者近似稳定状态时的降深Srec取值为：稳定或者近似稳定状态，指在所有的观测井中随时时间的进行降深的变化为零或者趋近于零；因此，当抽水达到稳定或者近似稳定时，Spump1r，t的一阶导数为0，则：当t'＝t时，因此，当时，一阶导数的最小值与相等，都等于上述步骤S3中当抽水达到准稳定状态时的降深Srec取值为：由于边界的影响，在实际抽水过程中很难达到稳定状态，要达到稳定状态需要大量充足的时间，相比，抽水状态达到某一规定程度时具有定量属性的准稳定状态更容易建立；降深Srecr，t’的一阶导数如下：于准稳定状态下，即，的相对变化小于1％的情况下，近似等于常数则，Srecr，t’的二阶导数如下：在时，达到最小值与抽水阶段的走时相同；的值计算如下：方程表明，水位恢复阶段的衰减与抽水阶段不完全相等，利用水位恢复阶段的衰减时，应增加近似常数本发明的有益之处在于：本发明的一种避免抽水实验中水位数据噪音给反演计算带来误差的方法，通过基于一维、均质和各向同性的含水层中的地下水流的控制方程，表述关于降深Spump的方程，结合实际的抽水状态，叠加引入水位恢复阶段的反向抽水；通过于公式推演，以Spump的一阶导数表述的水力衰减，分析不同状态下水位恢复阶段和抽水阶段的降深Srec的衰减具有相关的一致性；并基于上述的一致性，用恢复阶段的降深数据来代替实际抽水阶段的降深数据进行反演。本发明通过公式推导证明在抽水试验过程中，停止抽水后计算出的水位回升走时与抽水阶段水头下降的走时是一样的，因此使用水位回升的数据来做走时的计算，使用本发明的数据计算水力学走时，相对于使用水位下降数据计算的走时要精确很多，可有效避免抽水实验中的水位数据噪音给反演计算带来误差，具有很强的实用性和广泛的适用性。附图说明图1为本发明的降深-时间曲线图。图2为本发明的实施例的某次抽水试验的水位变化曲线图。具体实施方式以下结合附图和具体实施例对本发明作具体的介绍。在德国境内的实验基地上进行了上百次的抽水实验数据表明，所有的抽水阶段的水位降深数据都含有大量的噪音。如附图2所示，为其中一次抽水实验中的水位变化曲线图。从图中可见，抽水泵关闭前后的曲线有较大的差别；关闭前，因为降深本身绝对数值很小，所以导致信噪比非常差；关闭后，水位回升的曲线非常的平滑，不含噪音；可见抽水泵的运行，对数据有较大的干扰。本发明的一种避免水位数据噪音的方法，通过公式推演，得到上述的式10-式12，表明：在抽水试验过程中，停止抽水后计算出的水位回升走时与实际抽水阶段水头降深下降的走时是一样的；所以，可以使用水位恢复阶段的数据来做走时的计算。且，水位恢复阶段和抽水阶段的降深Srec的取值具有相关的一致性：水位恢复阶段的衰减与抽水阶段不完全相等，利用水位恢复阶段的衰减时，增加近似常数基于上述的一致性，用水位恢复阶段的降深数据来代替实际抽水阶段的降深数据进行反演，可以有效避免抽水实验中的水位数据噪音给反演计算带来误差。因此，使用此数据计算水力学走时，相对于使用水位下降数据计算的走时要精确很多。以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解，上述实施例不以任何形式限制本发明，凡采用等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

权利要求：1.避免抽水实验中水位数据噪音给反演计算带来误差的方法，其特征在于，包括以下步骤：S1、基于叠加原理，将抽水稳定后的降深Srec视为抽水和水位恢复的叠加，同时，将水位恢复视为反向抽水，抽水和反向抽水的降深分别为Spump1和Spump2，即Srec＝Spump1+Spump2；S2、基于一维、均质和各向同性的含水层中的地下水流的控制方程，表述关于降深Spump的方程，并以Spump的一阶导数表述水力衰减，二阶导数表述最大水力衰减，同时，获得最大值及最大值条件；S3、通过比较分别当抽水达到稳定或者近似稳定状态、达到准稳定状态时的降深Srec的取值，获得水位恢复阶段和抽水阶段的降深Srec的衰减具有相关的一致性；S4、基于上述的一致性，用水位恢复阶段的降深数据来代替实际抽水阶段的降深数据进行反演，以避免抽水实验中的水位数据噪音给反演计算带来误差。2.根据权利要求1所述的避免抽水实验中水位数据噪音给反演计算带来误差的方法，其特征在于，所述步骤S2中表述关于降深Spump的方程，步骤如下：A1、基于一维、均质和各向同性的含水层中的地下水流的控制方程，如下所示：式1中，K为渗透系数，SS为储水率，H为水头，r为任意点距离抽水井中心轴的水平距离，rp为P点距离抽水井中心轴的水平距离，Qrp，t为rp初始单位体积内的抽水率，t为抽水开始后的时间；A2、将方程式1，可以为描述关于降深的方程：式2中，s＝H-H0为降深，H0为初始水头；其中，表示了抽水前的流场散度，当初始水头散度稳定或者准稳定时，其值为零，则，方程式2可以写为：式3中，为水力扩散系数；当抽水率Q为常数，初始降深为零，并且区域无限大时，方程的解为：式中，Spump时抽水期间的降深，t为抽水开始后的时间。3.根据权利要求1所述的避免抽水实验中水位数据噪音给反演计算带来误差的方法，其特征在于，所述步骤S2中，以Spump的一阶导数表述水力衰减，二阶导数表述最大水力衰减，则：且，的最大值在为零时得出，为了满足这个条件，令此外，的最大值为：4.根据权利要求1所述的避免抽水实验中水位数据噪音给反演计算带来误差的方法，其特征在于，所述步骤S1中，基于叠加原理，将降深Srec视为抽水和反向抽水的叠加，即式8中，Spump1是当抽水率为常数Q时的降深，Spump2是虚拟的反向抽水中当抽水率为常数-Q时的降深；tp为停止抽水时的时间，t’为停止抽水后的时间，即，t＝tp+t’。5.根据权利要求1所述的避免抽水实验中水位数据噪音给反演计算带来误差的方法，其特征在于，所述步骤S3中当抽水达到稳定或者近似稳定状态时的降深Srec取值为：稳定或者近似稳定状态，指在所有的观测井中随时时间的进行降深的变化为零或者趋近于零；因此，当抽水达到稳定或者近似稳定时，Spump1r，t的一阶导数为0，则：当t'＝t时，因此，当时，一阶导数的最小值与相等，都等于6.根据权利要求1所述的避免抽水实验中水位数据噪音给反演计算带来误差的方法，其特征在于，所述步骤S3中当抽水达到准稳定状态时的降深Srec取值为：由于边界的影响，在实际抽水过程中很难达到稳定状态，要达到稳定状态需要大量充足的时间，相比，抽水状态达到某一规定程度时具有定量属性的准稳定状态更容易建立；降深Srecr，t’的一阶导数如下：于准稳定状态下，即，的相对变化小于1％的情况下，近似等于常数则，Srecr，t’的二阶导数如下：在时，达到最小值与抽水阶段的走时相同；的值计算如下：方程表明，水位恢复阶段的衰减与抽水阶段不完全相等，利用水位恢复阶段的衰减时，应增加近似常数

百度查询： 河海大学 避免抽水实验中水位数据噪音给反演计算带来误差的方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD