关于二维弥散解析法中连续污染源解析式的理解，求大神指点

lewisaa

最近在做一个地下水的二级评价项目，根据新导则用的是连续污染源解析式

根据公式的解释整个公式的前边部分都是定值，只有e的指数是和污染物的运移位置x有关，而解析解的后半部分是有关贝塞尔函数和第一类越流系统井函数有关，根据贝塔给出的公式，贝塞尔函数也是一个跟运移位置有关的函数，那就是在固定位置，他的函数值是不变的，最后是第一类越流系统井函数
我把公式中的u换成导则给出公式中的u2t/4DL，把r/b换成贝塔，进行计算，计算工具用的matlab。
情况说完了，现在问题来了！！！第一类越流井函数计算后1000天和2000天的结果都是小数点后十几为或者20几位，对整个解析解的影响很小，以硝酸盐为例，入渗浓度是38mg/l，分别计算100天、1000天和2000天的运移范围，最后发现100天和1000天还有2000天的超标范围基本都在一个位置，很困惑啊！！！我觉得是不是mt这个值不是公式说的单位时间注入的质量，或者说随时时间的推移每日注入示踪剂的质量是在变化的，本人是环境专业，地下水认识比较浅薄！！希望大神给指条明路！
我引用条件分别给水500m3/d,硝酸盐浓度38mg/l，M=45m，u=1.5m/d，n=0.3，DL=25m2/d ,DT=5m2/d

lewisaa

没有人给予解答吗？难道现在都是用软件做的？

dfan222

你好，第一类越流井函数用matlab怎么计算的

icespring

同样郁闷啊，刚刚用vba编程计算井函数，马上就出现了跟你完全一样的问题，W（u，β）根本就没有东西，计算半天就是一个定值，导则里面也没有给解释，好歹给个井函数的正规表达式啊，我现在看到了两个不同的井函数，另一个是单参数，用于抽水的，似乎也不能用吧。

wx_JKYY4XfR

...................

icespring

参考一下这篇文章的算法，看是不是好一点，https://wenku.baidu.com/view/08c ... l?qq-pf-to=pcqq.c2c

lewisaa

icespring 发表于 2018-11-13 14:42
参考一下这篇文章的算法，看是不是好一点，https://wenku.baidu.com/view/08cc5f3b376baf1ffc4fad01.html?q ...

非常感谢，我真觉得闭门造车太痛苦了