明天你好 发表于 2013-5-12 20:23

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

本帖最后由 wangjie029 于 2013-5-13 21:04 编辑

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



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


附件


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




yue1988911 发表于 2013-5-13 09:04

拿而不留非礼也，感谢

左京生 发表于 2013-5-13 09:14

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

明天你好 发表于 2013-5-13 10:40

左京生 发表于 2013-5-13 09:14 static/image/common/back.gif
井点处降深不可能是负无穷大

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

左京生 发表于 2013-5-13 11:03

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:19

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


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

左京生 发表于 2013-5-13 11:58

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

r为开采井的半径

左京生 发表于 2013-5-13 12:32

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

liujun192616 发表于 2013-5-13 18:34

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

yyjoc 发表于 2013-5-13 19:11

感谢楼主，学习一下！

tlx1026 发表于 2013-5-13 20:14

明显感觉不一般:)

好好人 发表于 2013-5-14 07:07

谢谢分享。。。。。

lwl1961 发表于 2013-6-16 17:03

好东西，楼主厉害

tjjtjjtjja 发表于 2013-10-7 20:41

感谢楼主分享:)

liaoboquan 发表于 2014-3-19 16:34

下载学习一下
页: [1] 2
查看完整版本: 地下水导则F.6公式承压水非稳定计算程序1.10（测试版）修改

2022年环境影响评价师辅导