地下水导则F.6公式承压水非稳定计算程序1.10（测试版）修改

明天你好

本程序为《环境影响评价导则 地下水环境》（HJ 610-2011）中附录F中F.6的计算程序，非稳定流承压水开采井计算程序，本程序为测试版，仅供学习交流所用。程序经修改后界面如下：



使用方法见程序使用说明，附件带有实例。
计算结果可用surfer作图，见下：


附件
F.6承压水非稳定解析解（测试版）.rar (1.62 MB, 下载次数: 1504)

2013-5-13 21:04 上传

点击文件名下载附件

修改内容如下：
1 增加了 开采井的半径 输入；
2 解决了预测结果中出现无穷大的情况。

yue1988911

拿而不留非礼也，感谢

左京生

井点处降深不可能是负无穷大

明天你好

左京生 发表于 2013-5-13 09:14
井点处降深不可能是负无穷大

实际上不可能，但从导则推荐的这个模型，从数学求解的角度上来说是可以的。采用surfer制图的时候，这个无穷大是不参与插值的。

左京生

Sub 非稳定流承压水干扰井群法()
Dim tt, h, k, uu, t, rj, r, w As Double
Dim n, nj As Integer
Dim q(1 To 50), xc(1 To 50), yc(1 To 50), xj(1 To 50), yj(1 To 50), sj(1 To 50) As Double

  tt = Sheet2.Cells(2, 1).Value
  n = Sheet2.Cells(2, 2).Value
  h = Sheet2.Cells(2, 3).Value
  k = Sheet2.Cells(2, 4).Value
  t = Sheet2.Cells(2, 5).Value
  uu = Sheet2.Cells(2, 6).Value
  rj = Sheet2.Cells(2, 7).Value
  nj = Sheet2.Cells(2, 8).Value
     For i = 1 To n
        xc(i) = Sheet2.Cells(5 + i, 1).Value
        yc(i) = Sheet2.Cells(5 + i, 2).Value
        q(i) = Sheet2.Cells(5 + i, 3).Value
     Next i
     For i = 1 To nj
        xj(i) = Sheet2.Cells(5 + i, 5).Value
        yj(i) = Sheet2.Cells(5 + i, 6).Value
     Next i

       For i = 1 To 50
       Sheet2.Cells(i + 5, 8).Value = ""
       Next i
      
'--------------------Theis计算-------------------------
      For i = 1 To nj
        sq = 0
        For j = 1 To n
          xx = xj(i) - xc(j)
          yy = yj(i) - yc(j)
          r = Sqr(xx * xx + yy * yy) + rj
          Call wu(uu, t, tt, r, w)  '--调用Theis井函数--
          sq = sq + q(j) * w / (2 * k * 3.14)
        Next j
         
        sj(i) = sq / (2 * h)
     Next i
  '-----------Excel 输出-----------

       For i = 1 To nj
       Sheet2.Cells(i + 5, 8).Value = sj(i)
       Next i



End Sub
Sub 非稳定流潜水干扰井群法()
Dim tt, h, k, uu, t, rj, r, w As Double
Dim n, nj As Integer
Dim q(1 To 50), xc(1 To 50), yc(1 To 50), xj(1 To 50), yj(1 To 50), sj(1 To 50) As Double

  tt = Sheet2.Cells(2, 1).Value
  n = Sheet2.Cells(2, 2).Value
  h = Sheet2.Cells(2, 3).Value
  k = Sheet2.Cells(2, 4).Value
  t = Sheet2.Cells(2, 5).Value
  uu = Sheet2.Cells(2, 6).Value
  rj = Sheet2.Cells(2, 7).Value
  nj = Sheet2.Cells(2, 8).Value
     For i = 1 To n
        xc(i) = Sheet2.Cells(5 + i, 1).Value
        yc(i) = Sheet2.Cells(5 + i, 2).Value
        q(i) = Sheet2.Cells(5 + i, 3).Value
     Next i
     For i = 1 To nj
        xj(i) = Sheet2.Cells(5 + i, 5).Value
        yj(i) = Sheet2.Cells(5 + i, 6).Value
     Next i

       For i = 1 To 50
       Sheet2.Cells(i + 5, 8).Value = ""
       Next i
      
'--------------------Theis计算-------------------------
     For i = 1 To nj
        sq = 0
        For j = 1 To n
          xx = xj(i) - xc(j)
          yy = yj(i) - yc(j)
          r = Sqr(xx * xx + yy * yy) + rj
          Call wu(uu, t, tt, r, w)  '--调用Theis井函数--
          sq = sq + q(j) * w / (2 * k * 3.14)
        Next j
           sj(i) = h - Sqr(h * h - sq)
     Next i

  '-----------Excel 输出-----------

       For i = 1 To nj
       Sheet2.Cells(i + 5, 8).Value = sj(i)
       Next i

End Sub

明天你好

左京生 发表于 2013-5-13 11:03
Sub 非稳定流承压水干扰井群法()
Dim tt, h, k, uu, t, rj, r, w As Double
Dim n, nj As Integer

当r=0的时候，w(u)应该等于多少呢？
或者说r不可能等于0，r最小值为开采井的半径？

左京生

wangjie029 发表于 2013-5-13 11:19
当r=0的时候，w(u)应该等于多少呢？
或者说r不可能等于0，r最小值为开采井的半径？

r为开采井的半径

左京生

改进一下加个开采井的半径，把范围和井数放开，是非常好的程序供大家使用。

liujun192616

好东西，感谢楼主。学习一下

yyjoc

感谢楼主，学习一下！

tlx1026

明显感觉不一般

好好人

谢谢分享。。。。。

lwl1961

好东西，楼主厉害

tjjtjjtjja

感谢楼主分享

liaoboquan

下载学习一下