ÁP DỤNG PHẦN MỀM MẠNG THẦN KINH NHÂN TẠO DỰ BÁO LŨ HẠ LƯU SÔNG THU BỒN – VU GIA ppsx

Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008

137
ÁP DỤNG PHẦN MỀM MẠNG THẦN KINH NHÂN TẠO
DỰ BÁO LŨ HẠ LƯU SÔNG THU BỒN – VU GIA
APPLYING ARTIFICIAL NEURAL NETWORKS SOFTWARE FOR FLOODING
PREDICTION IN THU BON – VU GIA RIVERS

SVTH: PHAN TẤN PHÁT
Lớp: 03X2C, Khoa XDTLTĐ, Trường ĐHBK,ĐHĐN
DHKH: GS.TS NGUYỄN THẾ HÙNG.
Khoa XDTLTĐ, Trường Đại học Bách khoa, ĐHĐN

TÓM TẮT:
Đề tài này nghiên cứu, áp dụng phần mềm mạng thần kinh nhân tạo (ANNs) để dự báo lũ ở hạ
lưu sông Thu Bồn – Vu Gia. Phầ n mề m ANNs có thể đượ c sử dụ ng tố t như mộ t công cụ trong
việ c dự bá o lũ mà căn nguyên củ a kỹ thuậ t nà y đượ c m ô phng theo hoạt đng ca b no
con ngườ i chứ a hà ng tỉ neurons v ới các mi quan h ni liền nhau ; khả năng dự báo, mô
phng ca mô hình là dựa vào khả năng học những kinh nghim trong quá khứ.
Đề tài so sánh s liu tính toán và dự báo lũ trên hạ lưu sông Thu Bồn – Vu Gia cho thấy có
sự phù hợp tt; và đề nghị áp dụng phần mềm ANNs trong vic dự báo lũ trên các sông ca
khu vực miền Trung.
ABSTRACT:
In this paper, author applly Artificial Neural Networks (ANNs) software to simulate and predict
flood in downstream Thu Bon – Vu Gia Rivers. ANNs softwares can be used as a good tool to
predict flood in rivers that the origin of this technique can be traced to the functioning of the
human brain which contains billions of neurons and their interconnection based on learning
abilities from past experiences. This paper compare the calculation results by data collections
in downstream Thu Bon – Vu Gia rivers with a best agreement; and author suggest to use
ANNs software for predicting flood in many other rivers of central Viet Nam.

1. Đặt vấn đề
1.1. Tổng quan các phương pháp dự báo lũ trên sông hiện nay
Hiện nay có nhiều mô hình dự báo lũ trên mạng lưới sông; trong số đó có thể chia làm
hai nhóm chính là:
- mô hình thuỷ lực thường chính xác hơn, nhưng cần phải có đầy đủ các số liệu đầu vào như
địa hình, độ nhám lòng dẫn
- mô hình thuỷ văn không đòi hỏi số liệu đầu vào đầy đủ như mô hình thuỷ lực, nhưng bị hạn
chế khi có sự thay đổi địa hình, độ nhám lòng dẫn mà mô hình chưa thể hiệu chỉnh về ảnh
hưởng của nó; cũng như tác động của biên hạ lưu.
1.2. Lựa chọn mô hình toán dự báo
Nhiều lúc trong tính toán dự báo, không có được đầy đủ các số liệu địa hình, nên
không thể áp dụng các mô hình thuỷ lực. Do đó phải sử dụng các mô hình thuỷ văn để dự báo.
Trong bài báo nầy, tác giả thử áp dụng mạng ANNs để dự báo mực nước lũ hạ lưu
sông.
Với mô hình mạng thần kinh nhân tạo có một số ưu điểm và nhược điểm như sau:
a. Ưu Điểm:
- ANNs có khả năng “họ c” cá c mố i quan hệ giữ a cá c biế n đầ u và o và đầ u ra khi cá c
quy luậ t tự nhiên đượ c ưu tiên là không đượ c biế t đế n hoặ c không đượ c biế t mộ t cá ch chính
xác.
- ANNs là mộ t công cụ hữ u í ch trong việ c mô phỏ ng cá c quá trình phc tạp.
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008


138
- Các thut toán thì đơn giản và không cần giải các phương trình vi phân tng phần
phứ c tạ p vớ i cá c bà i toá n kè m theo giố ng như tính không ổ n định củ a cá c thuậ t toá n.
- ANNs là m việ c tố t mộ t khi cá c tậ p hợ p luyệ n là không đầ y đủ hoặ c chứ a nhiễ u.
- ANNs có khả năng phỏ ng theo cá c giả i phá p vượ t thờ i gian.
- Mộ t khi mạ ng đã đượ c luyệ n thì nó đượ c sử dụ ng mộ t cá ch dể dà ng.
b. Nhược Điểm:

- ANNs yêu cầ u dữ liệ u phả i đầ y đủ cả về số lượ ng lẫ n chấ t lượ ng . Đây là yêu cầ u
quan trọ ng vớ i tấ t cả kỹ thuậ t mô phỏ ng và ANNs không thể bị loạ i bỏ .
- Các cách hướng dẫn để chn cấu trúc mạng cho ph hợp với các bài toá n cũ ng không
đượ c tìm thấ y.
2. Cơ sở lý thuyết của mô hình
2.1. Lan truyề n tiế n củ a ANNs
Quá trình lan truyền tiến tính giá trị các nút xuất t mẫu nhp vào mạng. Tiế n trình nà y
đượ c sử dụ ng trong hai tình huố ng: khi luyệ n mạ ng và khi sử dụ ng mạ ng. Khi luyệ n mạ ng, lan
truyề n tiế n đượ c sử dụ ng lặ p đi lặ p lạ i từ mẫ u nà y đế n mẫ u khá c từ thế hệ nà y sang thế hệ
khác cho đến khi trng số đạt đến giá trị thích hợp . Trái lại, khi sử dụ ng mạ ng lan tru yề n tiế n
chỉ được thực thi một lần cho tng mẫu nhp Tuy nhiên, d trong trạng thái luyện mạng hay
sử dụ ng mạ ng, các thao tác trong thủ tục lan truyến tiến là như nhau.
Một mạng lan truyền tiến của ANNs có một lớp nhp (lớp input), một lớp xuất (lớp output) và
một hoặc nhiều lớp ẩn ở giữa lớp nhp và lớp xuất. Mỗi neuron trong một lớp thì kết nối đến
tất cả các neuron của lớp kế tiếp ngay sau nó và các neuron trong cng một lớp không liên kết
với nhau.
Một neuron đặc trưng được thể hiện trong hình 2. Mỗi neuron nhn tín hiệu t mỗi
neuron trong lớp liền trước. Mỗi neuron trong lớp ẩn và lớp xuất, mỗi nút nhp thì được nhân
với các trng số của cung liên kết giữa nút nhp và nút ẩn, tích số của chúng được cộng dồn
lại. Kế tiếp, một hàm truyền được áp dụng trên tổng trng hóa này cng với một ngưỡng của
nút ẩn đó để cho ra giá trị thực của nút ẩn. Các trng số được đưa đến neuron th J
th
trong một
lớp có dạng vectơ trng W
j
= (w
1j
,…w
ij
,…,w

nj
), với w
ij
là trng số của cung liên kết t neuron
th i
th
trong lớp liền trước đến neuron trong lớp hiện tại.
Sau khi nén tổng trng hóa của nó, đến lượt mình, mỗi nút ẩn sẽ gửi kết quả đến tất cả
các nút xuất.

Output Variables
Output Layer
Hidden Layer
Input Layer
Input Variables

j
y
j
b
j
Xn
Xi
X
2
X
1

Hình 1 Hình 2
Hình 1: Biểu đồ thể hiện mạng lan truyền tiến 3 lớp của ANNs

Hình 2: Quá trình hoạt động của một nút đặc trưng của ANNs
Mỗi nút xuất thực hiện các thao tác tương tự như đã thực hiện trong nút ẩn để cho ra
giá trị kết xuất của nút xuất. Gía trị của các nút xuất chính là giá trị thực, nghĩa là giá trị của
các biến phụ thuộc cần xác định.
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008

139
Hầu hết người ta thường sử dụng hàm truyền mà về mặt hình hc, đồ thị của nó có
dạng đường cong dạng chữ S tăng đều đặn, được gi là hàm Sigmoid hay hàm logistic. Hàm
logistic có các đặc điểm sau:
+ Hàm truyền là hàm bị chặn.
+ Hàm truyền là hàm đơn điệu tăng.
+ Hàm truyền là hàm liên tục và trơn.
Về mặt toán hc, hàm logistic được định nghĩa như sau:

1
y
1
1
j
z
e


(1)
Với z = w
ij

Chàûn dæåïi
Chàûn trãn

0.5
f(x)
1
0
-15
-10 -5
15105

Hình 3 : Hàm logistic

2.2. Sự chuẩn hóa dữ liệu
Trước khi áp dụng ANNs, dữ liệu đầu vào phải được chuẩn hóa rơi vào trong khoảng
[0,1]. Một biến tượng trưng là lưu lượng Q, biên này có thể thay đổi t 0 đến một vài giá trị
lớn nhất Q
max
có thể được chuẩn hóa bởi công thưc sau :
Q
s
= Q/Q
max
Với Q
s
là lưu lượng chuẩn hóa.
2.3. Thuật toán lan truyền ngược sai số
Thut toán lan truyền ngược sai số (BP) dựa vào việc khái quát các quy lut delta đã
được đưa ra bởi Rumelhart (1986) để hiệu chỉnh các trng số của các cung liên kết trong quá
trình luyện mạng.
Lan truyền ngược là một phương pháp cho phép xác định tp trng tốt nhất của mạng
giải một bài toán được cho. Việc áp dụng phương pháp lan truyền ngược là một quá trình lặp
đi lặp lại nhiều lần hai tiến trình chính: lan truyền tiến để thực hiện ánh xạ và lan truyền ngược

sai số để cp nht các trng.
Các trng của mạng là các hệ số của mô hình. Phương pháp giảm gradient được dng
để cp nht những hệ số này sao cho giảm thiểu sai số của mô hình. Sai số được đo bằng
phương pháp sai số trung bình bình phương là phương pháp thường được sử dụng để xây dựng
các mô hình.
Sai số trung bình bình phương thường được sử dụng để đo lường sự trng khớp giữa
ánh xạ cần xây dựng với hàm đích được cho trước (qua tp mẫu).

Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008


140

2
11
()
Nm
pn pn
pn
E T O



(3)

Với T
pn
là giá trị đích của neuron th n
th
của mẫu th p

th
;
O
pn
là giá trị đầu ra của neuron th n
th
của mẫu th p
th
;
N là tổng số mẫu ;
m là tổng số neuron đầu ra.
Tăng trng số của các nút kết nối t i đến j ở bước lan truyền th n, Δw
ij
(n) được tính
bởi công thc:

**
ij ij
ij
w ( ) w ( 1)
w
E
nn


     

(4)
Với α và ε là hệ số momen và tốc độ hc.
* Bắt đầu tính toán với giá trị ban đầu của các trng số.Các bước của thut toán như sau:

+ Bước 1: Thực hiện các bước t 2-9 cho đến khi ngng các điều kiện là ph hợp.
+ Bước 2: Đối với mỗi cặp luyện của tp dữ liệu, tiến hành các bước t 3-8.
(1) Lan truyền tiến:
+ Bước 3: Mỗi neuron đầu vào (X
i
, I = 1, 2, …, n) nhn tín hiệu đầu vào x
i
và gửi nó đến
tất cả các nút trong lớp kế tiếp (lớp ẩn).
+ Bước 4: Mỗi neuron trong lớp ẩn (Z
j,
j = 1, 2, …, p) cộng tổng các trng số của tín hiệu
đầu vào.
Zin
j
= v
oj
+ ∑x
i
v
ij
với i = 1, 2, ….,n (5)
Với v
ij
các trng số của cung liên kết
V
oj
là giá trị nghiêng
Hàm hoạt động thì được áp dụng để tính giá trị đầu ra.
Z

j
= f(Zin
j
) (6)
+ Bước 5: Mỗi neuron trong lớp xuất (Y
k
, k = 1, 2, ….,m) cộng tổng các trng số của tín
hiệu đầu vào.
Yin
k
= w
ok
+ ∑z
i
w
ij
với j = 1, 2, , p (7)
Tiếp tục hàm hoạt động được áp dụng để tính giá trị đầu ra.
Y
k
= f(Yin
k
) (8)
(2) Lan truyền ngược sai số
+ Bước 6: Mỗi neuron trong lớp xuất (Y
k
, k = 1, 2, …, m) tính sai số sử dụng mẫu đích
ph hợp với các mẫu luyện đầu vào.
δ
k

= (t
k
– y
k
).f(Yin
k
) (9)
Tính số hạng hiệu chỉnh trng số (để cp nht w
ij
)

w
jk k j
Z

  
(10)
Tính số hạng hiệu chỉnh độ dốc (để cp nht w
ok
cuối cng)

w.
ok k

  
(11)
Và gửi δ
k
đến các nút trong lớp liền trước.
+ Bước 7: Mỗi nút trong lớp ẩn (Z

j
, j = 1, 2, …, p) cộng tổng delta đầu vào (t các nút
trong lớp kế tiếp)
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008

141

w
j k jk
in



với k = 1, 2, …, m (12)
Được nhân bởi đạo hàm của hàm hoạt động để tính số hạng sai số.

. ( )
j j j
in f Zin


(13)
Tính số hạng hiệu chỉnh trng số (được sử dụng để cp nht v
ij
sau):

ij
.
jj
vx


  
(14)
Và tính số hạng hiệu chỉnh độ giảm dốc (được sử dụng để cp nht v
j
sau):

oj
.
j
v


(15)
(3) Cp nht các trng số và độ dốc:
+ Bước 8: Mỗi neuron trong lớp xuất, Y
k
, k = 1, 2, …, m, cp nht độ giảm dốc và trng
số (j = 0, 1, …, p):

w ( ) w ( ) w
jk jk ik
new old  
(16)

Mỗi nút ẩn (Z
j
, j = 1, 2, …, p) cp nht độ giảm dốc và trng số (i = 1, 2, …, n):

ij ij

( ) ( )
ij
new v old vv   
(17)
+ Bước 9. Kiểm tra lại nếu các điều kiện ngng là an toàn
2.4. Các yêu cầu chung cho các file đầu vào
Nhóm tài liệu thủy văn:
1. Thời đoạn tính toán: liệt tài liệu thu thp được và thời gian đo đạt mực nước hạ lưu
sông cần dự báo.
2. Các điều kiện biên gồm:
a. Các điều kiện biên trên:
- Lưu lượng thượng nguồn sông
- Lượng mưa thượng nguồn sông
b. Các điều kiện biên dưới: mực nước tại vị trí cần dự báo
2.5. Kiểm soát kết quả đầu ra của mô hình
Đầu ra của mô hình là mực nước dự báo lũ do mô hình tính tại hạ lưu sông Thu Bồn và
sông Vu Gia
3. Áp dụng mô hình mạng thần kinh nhân tạo (ANNs) để dự báo lũ trên sông Thu Bồn –
Vu Gia
3.1. Điều kiện địa lý tự nhiên lưu vực sông Thu Bồn và sông Vu Gia
Lưu vực hệ thống sông Thu Bồn – Vu Gia có vị trí địa lý nằm trong khoảng 14
o
56 đến
16
o
05, 107
o
13 đến 108
o
28 kinh độ Đông, diện tích lưu vực khoảng 10350km

2
, hệ thống sông
có hai cửa là Cửa Đại ở Hội An thuộc Quảng Nam và Của Hàn ở Đà Nẵng.
Địa hình chủ yếu là đồi núi chiếm 2/3 diện tích, phần còn lại là đồng bằng nhỏ hẹp
chạy dc ven biển. Điah hình thấp dần t Tây sang Đông.
3.2. Dự báo lũ trên sông Thu Bồn và sông Câu Lâu bằng mô hình mạng thần kinh nhân
tạo
3.2.1. Điều kiện biên
a. Biên lưu lượng
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008


142
Lưu lượng nước thượng nguồn sông Thu Bồn và sông Câu Lâu tại trạm đo Nông Son
và trạm Thành Mỹ.
b. Biên lượng mưa
Lượng mưa thượng nguồn sông Thu Bồn và sông Câu Lâu tại trạm đo Thành Mỹ và
trạm Nông Sơn trong thời đoạn tính toán.
c. Biên mực nước
Mực nước hạ lưu sông tại vị trí cần dự báo tương ng với lưu lượng và lượng mưu ở
thượng nguồn sông.


Hình 7:Cột nước đo tại hạ lưu sông Vu Gia Hình 8:Cột nước đo tại hạ lưu sông Thu Bồn
3.2.2. Kết quả chạy mô hình

Hình 9 Hình 10
Hình 9: Cột nước giữa thực đo và mô hình tính tại hạ lưu sông Vu Gia;
Hình 10: Cột nước giữa thực đo và mô hình tính tại hạ lưu sông Thu Bồn.

4. Kết luận và kiến nghị
4.1. Kết luận
Việc ng dụng phần mềm ANNs vào trong dự báo lũ gặp rất nhiều thun lợi. Nó có
khả năng dự báo lũ ở hạ lưu sông tương đối chính xác chỉ dựa vào tài liệu thủy văn mà không
cần tài liệu về địa hình địa chất.
Qua kết quả chạy mô hình ta thấy, đường quá trình mực nước lên xuống tại hạ lưu sông
Thu Bồn và sông Vu Gia giữa thực đo và tính toán là tương đối ph hợp.
4.2. Kiến nghị
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008

143
Nước ta là một nước nông nghiệp, luôn chịu ảnh hưởng của thiên tai, hạn hán, lũ lụt
liên tục; vì vy “thủy lợi” là vấn đề quan trng trong sự nghiệp phát triển của đất nước.
Trong những năm qua ngành Thủy Lợi đã có những thành tích đáng kể trong nghiên
cu và áp dụng các tiến bộ khoa hc kỹ thut trong xây dựng công trình. Tuy nhiên, một số
lĩnh vực tính toán vẫn còn yếu. Công nghệ khảo sát, dự báo lũ của ta còn nhiều hạn chế. Vì
vy cần đưa phần mềm ANNs vào dự báo lũ trong điều kiện thiếu các tài liệu về địa hình, độ
nhám lòng sông… để phục vụ cho công tác xây dựng công trình và phòng chống thiên tai.

TÀI LIỆU THAM KHẢO

[1] TS Nguyễn Đình Thúc, Hoàng Đc Hải, Giáo trình mạng trí tuệ nhân tạo mạng Nơron
phương pháp và ứng dụng, Nhà xuất bản giáo dục năm 2000
[2] Trung tâm quân đội Mỹ, Water resources systems planning and management