matlab井函数,用MATLAB函数解承压含水层完整井非稳定流问题

用 MATLAB函数解承压含水层 完整井非稳定流 问题 李瑞 周德亮 (辽宁师范大学数学学院 辽宁大连 116029) 随着时代的进步，地下水数值模拟也开发出自己的研究软件， 用 MATLAB这个数学软件来研究地下水数值模拟也是可行的。 承压含水层完整井非稳定流可归纳为非线性偏微分方程的 求解问题 ，用 MATLAB编写的有限元或有限差分等算法的程序来 计算偏微分方程需要很大的工作量。所以对于一些偏微分方程的 模型，MATLAB提供了两种方法来解决这些问题。一种方法是用 PDE工具箱来解决特殊的偏微分方程求解问题，但它有较大的局 限性，只能求解部分的椭圆型、抛物型 、双曲型、特征值型的二阶 偏微分方程问题。另一种方法是采用 MATLAB中的函数来解决偏 微分方程求解问题，MATLAB已经编写好了这些函数的程序 ，在 脚本文件中直接调用就可以，这种方法可以解决不同类型的偏微 分方程 ，应用范围很广泛 。 1 pdepe函数 用法 根据承压含水层完整井非稳定流的方程模型，可采用 MAT— LAB中的函数 pdepe来求解。pdepe是用来求解一维偏微分方程 或方程组定解问题的命令函数，编写程序时，先布置结点，冉调用 函数。在解决定解问题时需先化成标准形式： c “， ax ) = u， ))+s “， )， ∈Q {u(x，t0)=u0(x) Ip ，f，“)+q ， 1舡， ，M， )=0,x∈an 1 m 其调用格式为 ： sol=pdepe(m，pdefun，icfun，bcfun，xmesh，tspan) m表 示标准方 程 中的指数 ，xmesh表 示标准 方程 中的 ，tspan 表示时间，其中 xmesh和 tmesh应取长度大于等于 3的单调递增 一 维数组。 pdefun是表示偏微分方程的函数，即生成， 成c “ )， 詈 ， )s ， ) 函数。编写格式为 ：【c，f,s]=pdefun(x⋯t u du)，其中，x,t是对应标准方 程中的与，u对应所求解 ，du对应的一阶导数，给定几个输入变量 即可表示 出 c,f,s这 三个 函数 。 icfun是表示偏微分方程初值条件的函数，定义初值 。 )。编 写格式 为 ：uO=icfun(x)。 bcfun是表示偏微分方程边界条件的函数，即生成p ，t，H)和 q ，t)函数。编写格式为 ：[pl，ql，pr，qr]=bcfun(xl，ul，Yd"，Ur，t)，其中 xl、ul和 XI"、ur分别对应左边界和右边界方程中的与所求解 ；pl、ql和 pr、qr 分别对应标准边界条件中左边界和右边界的函数。 对于偏微分方程定解问题来说，pdepe函数的输出参数 sol是 一 个二维数组，表示所求解 ，即 “( )对应的 ( )和 ￡ )的解为sol(i， i)；对于元偏微方程组定解问题来说 ，u=(uI’u：⋯， n)，输出参数 sol 是一个三维数组 ，solf ，n)表示方程组的解 ，即u 对应 ( )和 )时 的解为 sol(i,j，k)。 2承压 含水层 完整 井非 稳定流模 型 求解 2．1定降深模型 求解 考虑一均质、等向的承压含水层定流量完整井非稳定流问 题，井里的水位保持定降深。在这个模型中，用降深作为待求函 数，则模型的定解问题可以写成： 初始条件 边界条件 a‘S ． 1 aS — S aS ar2。r ar KM at S(r，t)l ：0 r S(r，f)I =0 lS(r，t)lr= =S 为厚度， 为渗透系数 ，Js 为贮水系数，5 为抽水前降深，R为 抽水井的