matlab 泰斯函数,Theis解析模型的数值解及其在MatLab中的实现

江苏地质，３１(３)，２４２—２４６，２００７

Ｔｈｅｉｓ解析模型的数值解及其在ＭａｔＬａｂ中的实现

李林子。朱国荣，江思珉，闵望

(南京大学地球科学系，江苏南京２１００９３)

摘要：讨论了如何应用数值积分方法求解含水层中地下水非稳定流抽水试验的解析模型以获得含水层的水文地质参数．并以Ｔｈｅｉｓ公式为例讨论求解原理。引入了ＭａｔＬａｂ的数值积分工具和优化工具，轻松实现了对水文地质参数的求解，此举既保证了解的可靠性和唯一性，也大大提高了方法的可操作性。最后以一个具体的承压水非稳定流抽水试验为例验证了这个方法，所得的结果与应用配线法得到的结果基本一致，从而肯定了方法的应用价值。

研究结果尽管是针对Ｔｈｅｉｓ模型进行的，但求解原理同样适用于其他各种模型。

关键词：ＭａｌＬａｂ；数值积分；Ｔｈｅｌｓ解析模型；水文地质参数；优化

中图分类号：Ｐ６２８文献标识码：Ａ文章编号：１００３－６４７４(２００７)０３－０２４２—０５

１研究思路

通过对含水层进行地下水非稳定抽水试验以求解水文地质参数是一种水文地质基础研究的通用方法。然而，在地下水非稳定流理论问世的４０余年里，对解析模型的求解一直没有比较简单易行且可靠的方法，教科书上要么采用配线法来获得水文地质参数，要么求解析式的近似解，如Ｔｈｅｉｓ解析模型的Ｊａｃｏｂ近似公式就是其中之一。

笔者在这里以Ｔｈｅｉｓ解析模型为例，提出了直接求解其解析解的方法，旨在克服以上的各种问题。

已知地下水非稳定流的Ｔｈｅｉｓ公式为

ｓ＝器ｒｅ(“)＝器ｆ争ｕ(１)其中

一丝

“一４ｎ

(２)式(１)和式(２)中，，为观测孔到抽水井的距离(量纲为Ｌ)，弘’为承压含水层贮水系数(无量纲)，ｒ为承压含水层导水系数(量纲为Ｌ２Ｔ“)，ｔ为抽水持续时间(量纲为Ｔ)，Ｑ为单位时间内从抽水井中抽出的水量(量纲为ｒＴ“)，ｓ为观测孔中对应于观测时间ｔ时的水位下降(量纲为Ｌ)。对式(１)的求解主要难点是处理积分项。求解过程存在以下两个问题。

①对定积分的求解，只要确定了积分上下限，就可采用数学上非常成熟的各种定积分近似计算公式来进行，例如矩形法、梯形法、抛物线法等，但积分区问无法事先给定。

②积分区间［ｕ，＋＊］的确定，根据(２)式，可知决定“的各要素中，ｒ和ｔ均为已知，丁和灿。为待求的水文地质参数。

第一个问题的答案将随着对第二个问题的求解而得到。解决第二个问题的思路是采用水文地质参数寻优方法确定ｒ和肛’。为此我们首先选定一组参数，将其代入(２)式获得ｕ，然后将“代人(１)式进行水位降深验证。再根据验证误差调整参数值，重复这个过程，最终求出导水系数和贮水系数，为模型的数值解提供必要的参数。

得益于计算机软件技术的进步，一个被誉为“工程师语言”的软件——ＭａｔＬａｂ出现了，该软件是一个集数值分析、符号运算、图形图像处理、系统控制和程序语言控制为一体的数学软件，特别是其工具箱中集成的各种数值计算工具包括解决问题所需要的数值积分工具和优化工具，为我们获得Ｔｈｅｉｓ模型的数值解提供了便利和可能。

收稿１３期：２００７—０４—１２；编辑：陆李萍

基金项目：南京大学“９８５工程”二期“教学创新和隶质培养”研究项目资助

作者简介：李林子(１９８７一)．女，云南昆明人，水文与水资源工程专业．