CN105181895A - 利用海岸带多个观测孔潮汐效应地下水位信息确定含水层参数的方法 - Google Patents

Abstract

本发明公开了一种利用海岸带多个观测孔潮汐效应地下水位信息确定含水层参数的方法。描述滨海承压含水层地下水一维非稳定流动数学模型的解析解表明，在海潮影响下海岸带地下水位潮汐效率随距海岸的距离呈负指数函数衰减，地下水位与海潮的滞后时间随距海岸的距离呈线性函数增加。利用海岸带多个观测孔受海潮影响发生波动的地下水位观测资料，绘制潮汐效率与距海岸的距离和滞后时间与距海岸的距离的散点图，分别用负指数函数和线性函数拟合，确定衰减常数和斜率，从而求得承压含水层的储水系数与导水系数的比值。利用潮汐效应地下水位观测数据确定含水层水文地质参数，这在一定程度上可以减少利用抽水试验求含水层参数的费用。

Description

利用海岸带多个观测孔潮汐效应地下水位信息确定含水层参数的方法

1.技术领域

本发明涉及利用海岸带多个观测孔潮汐效应地下水位信息确定承压含水层水文地质参数的方法，可以用来研究海岸带地下水动力特征和进行地下水资源评价、开发和保护，属于海岸带地下水动力学技术领域。

2.背景技术

由于月球、太阳对地球的引力作用，致使海洋表面出现周期性涨落的现象，这就是海潮。因月球、太阳对地球各处引力的不同，海平面有时出现较大的涨落(大潮)，有时出现较小的涨落(小潮)，如此周而复始。在与海水有联系的滨海含水层的井孔中，可以观测到地下水位随着潮起潮落而出现的有规律的波动，这种现象称为潮汐效应，在潜水含水层和承压含水层中均可以出现。研究海岸带潮汐效应，包括以下问题：①滨海含水层地下水位随潮汐波动的幅度大小和地下水位波动幅度随距海岸的距离的增大而呈现的变化特点，以及潮汐影响的最远距离，②滨海含水层地下水位波动的周期和与海潮的时间滞后，③描述滨海含水层地下水位波动的数学模型及水位动态预测，④利用潮汐效应地下水位动态观测资料确定含水层参数等。

海潮引起滨海含水层的潮汐效应，是海岸带重要的自然现象，蕴含着有关滨海含水层的丰富信息。通过观测地下水位变化和定量研究这种潮汐效应，可以掌握海潮对滨海含水层地下水的影响程度，了解含水层多孔介质和地下水的诸多性质。利用观测的地下水动态资料可以求得含水层水文地质参数(如水力传导系数、储水系数等)；或者已知陆地某点处的承压水水头动态及其它必要的参数，可以反求承压含水层顶板延伸到海底的长度。在设计海岸附近开采井的开采降深和抽水设备时，地下水潮汐效应也是设计者必须考虑的因素。潮汐效应对滨海含水层海岸带咸、淡水界面或过渡带的变化有着显著的影响，从而影响到海岸带地质环境的变化，因而在研究海水入侵问题时，也必须对此予以足够的重视。在天然条件下海平面的周而复始的涨落对海岸带咸淡水界面和地下水位产生重要影响，人为开采海岸带地下水则直接影响到海岸带咸淡水界面的移动和地下水位的变化。

当含水层与海水有直接水力联系时，海水面升降变化引起的压力波会在含水层中传播，致使地下水位出现象海潮一样有规律的波动。显然，海岸带地下水位的波动幅度要小于海水面的波动，并且随着远离海岸，地下水位的波动逐渐减弱，到达一定距离后，就观测不到地下水位的这种波动了。可以建立描述海岸带地下水流动的数学模型，对潮汐效应现象进行模拟。以下考虑的是滨海承压含水层的情形。

如图1所示的均质各向同性滨海承压含水层，假定承压含水层初始地下水位水平并且与海平面一致，海岸边界为垂直边界，忽略海水密度与淡水密度的差异，海潮的波动可以用正弦函数描述。

选取图1所示的坐标，原点在平均海平面与海岸线交界处，向内陆方向为正。描述滨海承压含水层地下水一维非稳定流动的控制方程和初始条件及边界条件可以表示为：

$\frac{{\∂}^{2}H}{\∂x^2}=\frac{S}{T}\frac{\∂H}{\∂t}---\left(1\right)$

H(x，0)＝0(1a)

$H(0,t)=H_{0}sin\left(\frac{2\mathrm{\π}t}{t_0}\right)---\left(1b\right)$

H(∞，t)＝0(1c)

式中：H为以平均海平面为基准面的水位标高(m)，S为承压含水层的储水系数(无量纲)，T为承压含水层的导水系数(m²/h)，x为距海岸的距离(m)，t为时间(h)，H₀为海平面变幅的一半(即振幅)(m)，t₀为海平面波动周期(h)。上述定解问题的解为(Jacob1950；Ferris1951；Werner等，1951；Ingersoll等，1954)：

$H(x,t)=H_0{e}^{-x\sqrt{\frac{\mathrm{\π}S}{t_{0}T}}}sin\left(\frac{2\mathrm{\π}t}{t_0}-x\sqrt{\frac{\mathrm{\π}S}{t_{0}T}}\right)---\left(2\right)$

式(2)可以用来描述在海潮影响下滨海承压含水层距海岸距离为x的任意点在任意时刻t时的水位。

3.发明内容

本发明公开了一种利用海岸带多个观测孔潮汐效应地下水位信息确定承压含水层水文地质参数的方法。其基本内容如下。

将式(2)与式(1b)进行比较，可以看出承压含水层距海岸线x处的水位随海潮的波动而出现相应的波动，而且具有衰减效应和滞后效应。由式(2)可知距海岸线x处的承压含水层水位波动的振(H_x)为：

$H_x=H_0{e}^{-x\sqrt{\frac{\mathrm{\π}S}{t_{0}T}}}---\left(3\right)$

从式(1b)和式(3)可以求得潮汐效率(TE，即相同时间内地下水位变幅与海潮变幅之比)为：

$TE=e^{-x\sqrt{\frac{\mathrm{\π}S}{t_{0}T}}}---\left(4\right)$

由式(2)可以求得滞后时间(t_L，即地下水位高峰或低谷落后于相应的海潮高峰或低谷的时间)为：

$t_L=x\sqrt{\frac{t_{0}S}{4\mathrm{\π}T}}---\left(5\right)$

式(4)表明潮汐效率随距离x呈负指数衰减，式(5)表明滞后时间随距离x呈线性增加。

假设承压含水层的导水系数T＝31.25m²/h，储水系数S＝4.5×10^-4，海潮的振幅H₀＝2.5m，周期t₀＝24.7h，则可以依据式(2)求得不同距离x处地下水位随时间的变化，图2给出了x＝0，100，200，500，1000，1500，2500m处的水位变化曲线。可以看出，在海潮波动的诱导下，海岸带任意距离x处承压含水层地下水位随时间也出现相应的波动。随着x的增加，地下水位的变幅逐渐减小，出现波峰或波谷的时间逐渐增大。

4.附图说明

图1为滨海承压含水层示意图。

图1中，H₀为海潮振幅。

图2为滨海承压含水层距海岸不同距离处地下水位波动曲线。

图3是海岸带利用多个观测孔潮汐效应地下水位数据绘制的潮汐效率和滞后时间随距海岸距离的变化关系图。

图3中，A-潮汐效率与距海岸距离的负指数衰减关系；B-滞后时间与距海岸距离的线性增加关系。

5.具体实施方式

假设海岸带含水层符合图1中的条件。由式(4)描述的潮汐效率随距海岸距离x的负指数衰减变化关系可以求出衰减常数β为：

$\mathrm{\β}=\sqrt{\frac{\mathrm{\π}S}{t_{0}T}}---\left(6\right)$

由式(5)描述的滞后时间随距海岸距离x的线性变化关系可以求出斜率b为：

$b=\sqrt{\frac{t_{0}S}{4\mathrm{\π}T}}---\left(7\right)$

由式(6)和式(7)可知，如果通过潮汐效应野外水位观测资料能够确定潮汐效率的衰减系数β和滞后时间的斜率b，并已知潮汐波动周期t₀，则承压含水层储水系数与导水系数的比值S/T便可以确定，即：

$S/T=\frac{{\mathrm{\β}}^2{t}_0}{\mathrm{\π}}---\left(8\right)$

$S/T=\frac{4{\mathrm{\πb}}^2}{t_0}---\left(9\right)$

如果承压含水层的储水系数和导水系数中的一个为已知，则可以求得另外一个。例如，已知含水层的渗透系数和厚度，就可以求得储水系数。一般来说，利用潮汐效应观测资料，不能同时确定承压含水层储水系数和导水系数两个参数。为了简单和方便期间，有时把β称为潮汐传播系数。潮汐传播系数描述的是潮汐引起的压力波在含水层中传播的容易程度，与含水层的储水系数、导水系数和潮汐周期有关。

为了确定潮汐效率的衰减系数β和滞后时间的斜率b，需要在多个观测孔中观测承压含水层潮汐效应地下水位的变化。通过绘制潮汐效率与据海岸距离的散点图和滞后时间与据海岸距离的散点图，分别用负指数函数和线性函数来拟合，从而确定潮汐效率的衰减系数β和滞后时间的斜率。再根据式(6)和式(7)求得S/T。这两种方法分别称为潮汐效率法和滞后时间法。

例如，在1992年9月18-22日在广西北海市北部海岸曾在10个观测孔开展潮汐效应地下水位观测。在观测期间没有降雨，气压效应也可以忽略不计。观测孔地下水位变化曲线具有近似正弦曲线波动。

根据观测资料确定潮汐效率，并在潮汐效率与据海岸距离的关系图上作散点图(图3A)。可以看出，潮汐效率随距离大体上呈负指数衰减。实际观测数据与理论结果有一些差别，这可能是受一些因素的影响所致，例如，实际海岸边界并非垂直边界、含水层并非严格均质、可能存在少量越流、隔水顶板延伸到海里一定距离等。由图3A可以推测在北海市北部在距海岸3300m以远的地下水位已不受海潮波动的影响。

图3B表示的是滞后时间与距离的关系散点图，两者之间大体上呈线性关系，实际资料与理论结果也有一些偏差，其原因与潮汐效应与距离的关系类似。

从图3A和图3B可以确定潮汐传播系数(衰减系数)β＝0.00188921/m和斜率b＝0.0001227d/m。对于t₀＝24.7h，按式(8)的潮汐效率法求得的比率S/T＝1.169×10^-6d/m²，而按公式(9)的滞后时间法求得比率S/T＝1.83×10^-7d/m²。已知含水层的渗透系数为30m/d，厚度为25m，即T＝750m²/d＝31.25m²/h，用潮汐效率法求得含水层的储水系数S为8.7675×10^-4，用滞后时间法求得含水层的储水系数S为1.3725×10^-4，二者的数值较为接近。

利用潮汐效应地下水位观测数据可以确定含水层参数，这在一定程度上可以减少抽水试验求参数的费用。这种方法不能同时求得含水层的储水系数和导水系数，另外要求观测孔尽可能多、观测时间间隔尽可能短。

参考文献

Jacob，CE.1950.Flowofgroundwater.In：H.Rose(ed)，Engineeringhydraulics.JohnWiley&SonInc.，NewYork，321-386

Ferris，JG.1951.Cyclicfluctuationsofwaterlevelasabasisfordeterminingaquifertransmissibility.InternationalAssociationofScientificHydrology，Publication33：148-155

Ingersoll，LR，Zobel，OJ，Ingersoll，AC.1954.HeatConduction，withEngineeringandGeologicalandotherApplication.TheUniversityofWisconsinPress，Madison

Werner，PW，Noren，D.1951.Progressivewavesinnonartesianaquifers.TransactionofAmericanGeophysicalUnion，32：238-244

ZhouX.2008.DeterminationofaquiferparametersbasedonmeasurementsoftidaleffectsonacoastalaquifernearBeihai，China.HydrologicalProcess，22：3176-3180

周训，金晓媚，梁四海，沈晔，张红梅.2010.地下水科学专论.北京：地质出版社。

Claims (5)

1.一种利用海岸带多个观测孔潮汐效应地下水位信息确定含水层参数的方法，其特征是：描述滨海承压含水层地下水一维非稳定流动数学模型的解析解表明，在海潮影响下海岸带地下水位潮汐效率随距海岸的距离呈负指数函数衰减，地下水位与海潮的滞后时间随距海岸的距离呈线性函数增加。利用衰减常数和直线斜率，可以求得承压含水层的储水系数与导水系数的比值。

2.根据权利要求1所述的计算方法，其特征是：利用海岸带多个观测孔潮汐效应地下水位观测数据，绘制地下水位潮汐效率与距海岸的距离的关系散点图，用负指数函数拟合，确定衰减常数。

3.根据权利要求1所述的计算方法，其特征是：利用海岸带多个观测孔潮汐效应地下水位观测数据，绘制地下水位对海潮的滞后时间与距海岸的距离的关系散点图，用线性函数拟合，确定直线的斜率。

4.根据权利要求1所述的计算方法，其特征是：根据负指数函数的衰减常数，求得承压含水层的储水系数与导水系数的比值，称为潮汐效率法。

5.根据权利要求1所述的计算方法，其特征是：根据线性函数的斜率，求得承压含水层的储水系数与导水系数的比值，称为滞后时间法。