Part 11 TÍNH TOÁN ỨNG DỤNG TRONG KTDK

©Copyright 2007

ỨNG DỤNG KHẢ NĂNG CỦA EXCEL 19

Bài toán mô hình Hillslope

3. Mô hình trong Excel

1

2

3

4

5

6

7

8

9

10

A B C D E F G H I J K L M NO P Q R S

NY EXPLANATION

i = 0 1 2 3 4 5 6 7 8 9 10 11 12 Dirichlech boundary condition

50.0 50.5 51.0 51.5 52.0 52.5 53.0 53.5 54.0 54.5 55.0 j = 1 No-flow boundary

2 Fictitious node

3 Inner cell

4 NX

5

6

7

2. Lưới sai phân

100 m

50 m

H = 50 55

Δx = Δy = 10 m

©Copyright 2007

ỨNG DỤNG KHẢ NĂNG CỦA EXCEL 20

Bài toán mô hình Hillslope

4. Kết quả phân bố cột áp thủy lực trong tầng chứa nước

H (m)

i = 1 2 3 4 5 6 7 8 9 10 11

j = 1 50.00 50.50 51.00 51.50 52.00 52.50 53.00 53.50 54.00 54.50 55.00

2 50.80 50.98 51.30 51.68 52.08 52.50 52.92 53.32 53.70 54.02 54.20

3 51.24 51.33 51.54 51.83 52.15 52.50 52.85 53.17 53.46 53.67 53.76

4 51.49 51.55 51.71 51.94 52.21 52.50 52.79 53.06 53.29 53.45 53.51

5 51.63 51.67 51.80 52.00 52.24 52.50 52.76 53.00 53.20 53.33 53.37

6 51.67 51.71 51.84 52.02 52.25 52.50 52.75 52.98 53.16 53.29 53.33

⇒ Kết quả tính toán thể hiện giá trị cột áp tại mỗi nút của

luới sai phân.

⇒ Sử dụng công cụ giải lặp của Excel

©Copyright 2007

ỨNG DỤNG KHẢ NĂNG CỦA EXCEL 21

Bài toán mô hình Hillslope

5. Biểu đồ thể hiện giá trị cột áp của trong tầng chứa nước

theo mặt cắt

48

49

50

51

52

53

54

(0, 50) (0, 40) (0, 30) (0, 20) (0, 10) (0, 0) (100, 0)

(x, y)

H (m) Analytic Numerical

©Copyright 2007

ỨNG DỤNG KHẢ NĂNG CỦA EXCEL 22

Giảm áp quanh một giếng bơm

‰ 1. Định nghĩa bài toán

ƒ Xét một tầng chứa nước hình

tròn có áp, chiều dày 30m, độ

dẫn suất thủy lực 13.3. Giả sử

biên của tầng thấm không ảnh

hưởng bởi lưu lượng giếng

bơm với lưu lượng 1000 m3/d

ở tâm của tầng chứa nước.

ƒ Tính toán sự phân bố cột áp

thủy lực trong tầng chứa.

1100 m

CONFINED AQUIFER

Well (Q = 1000 m3d-1)

©Copyright 2007

ỨNG DỤNG KHẢ NĂNG CỦA EXCEL 23

Giảm áp quanh một giếng bơm

3. Mô hình trong EXCEL

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

A B CDEFGH I J K LMN O P Q R S T

i = 0 1 2 3 4 5 6 7 8 9 10 11 12 EXPLANATION

30 30 30 30 30 30 30 30 30 30 30 30 j = 1

30 30 30 30 30 30 30 2 Dirichlech boundary condition

30 30 30 30 30 3 No-flow boundary

30 30 30 30 4 Fictitious node

30 30 30 5 Inner cell

30 30 6 NY Location of well

30 30 7

30 8

30 9

30 10

30 11

30 12

NX

2. Lưới sai phân

Giếng

©Copyright 2007

ỨNG DỤNG KHẢ NĂNG CỦA EXCEL 24

Giảm áp quanh một giếng bơm

H (m)

i = 1 2 3 4 5 6 7 8 9 10 11 12

30 30 30 30 30 30 30 30 30 30 30 30 j = 1

29.9149 29.9169 29.9232 29.9349 29.9558 30 30 30 30 30 30 30 2

29.8257 29.8295 29.8411 29.8605 29.8885 29.9246 29.9557 30 30 30 30 30 3

29.7287 29.7344 29.7511 29.7777 29.8129 29.8543 29.8982 29.9474 30 30 30 30 4

29.6204 29.6283 29.6512 29.6864 29.7309 29.7816 29.8354 29.8914 29.9474 30 30 30 5

29.4962 29.5074 29.5388 29.5858 29.6429 29.7056 29.7705 29.8354 29.8982 29.9557 30 30 6

29.3496 29.3662 29.4111 29.4750 29.5492 29.6274 29.7056 29.7816 29.8543 29.9246 30 30 7

29.1697 29.1966 29.2643 29.3538 29.4515 29.5492 29.6429 29.7309 29.8129 29.8885 29.9558 30 8

28.9360 28.9864 29.0955 29.2247 29.3538 29.4750 29.5858 29.6864 29.7777 29.8605 29.9349 30 9

28.6015 28.7176 28.9066 29.0955 29.2643 29.4111 29.5388 29.6512 29.7511 29.8411 29.9232 30 10

28.0346 28.3761 28.7176 28.9864 29.1966 29.3662 29.5074 29.6283 29.7344 29.8295 29.9169 30 11

26.7846 28.0346 28.6015 28.9360 29.1697 29.3496 29.4962 29.6204 29.7287 29.8257 29.9149 30 12

‰ Kết quả tính toán cột áp phân bố trong tầng chứa nước