魏海滨 谷洪彪* 翟泽宇 朱柏正 徐邑荣

（1、防灾科技学院，河北 三河 065201 2、河北省地震动力重点实验室，河北 三河 065201）

1 概述

含水层参数表征了含水层介质的特征，其在水文地质领域是非常重要的研究课题，同时也是研究各类水文地质问题的基础[1]。而含水层参数的计算除了传统的抽水试验和微水试验以外，Hsieh等（1987）利用井水位与含水层孔压变化之间的相位差推导了反演含水层导水系数的方法[2]。基于此理论基础，日本学者M.Ishiguro，Y.Tamura，T.Sato和M.Ooe利用Fortran77语言改进开发了Baytap-G程序[3]。该程序利用贝叶斯建模程序来分析包含潮汐和其他变化的时间序列，以此可进行：（1）估算振幅因子和相位差；（2）确定趋势项和计算其频谱；（3）对缺失数据进行插值和对阶梯变化量进行估计；（4）粗略查找、识别异常数据；（5）计算ABIC值（Akaike's Bayesian Information Criterion）以此分析模型好坏[3]。将得到的相位差等参数手动输入至Matlab[4]中即可计算出导水系数。此方法理论简单易懂，但实际操作过程中需要手动输入数据占用了一定的工作量。本次研究针对传统计算方法步骤繁琐，不利于研究者使用的问题对原有程序进行了改进和嫁接，使改进程序实现以Excel表格配置简单参数，利用Matlab一键运行即可得到所求参数，简化操作步骤，提高研究者工作效率。

2 理论依据及原有计算方法

开放式承压含水井水位响应受地球固体潮扩张引起的压力水头扰动，导致地下水的流动，当含水层渗透率高时，那么井中的振荡相位与受迫潮汐应变几乎协调一致；当含水层渗透率低时，井中水位振荡就较大的滞后于潮汐应变。……