Báo cáo nghiên cứu khoa học " Quy trình thử nghiệm dự báo trường dòng chảy, độ muối, nhiệt độ và mực nước tổng cộng cho khu vực Biển Đông bằng mô hình ROMS " potx

TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
362
_______
Quy trình thử nghiệm dự báo trường dòng chảy, độ muối,
nhiệt độ và mực nước tổng cộng cho khu vực Biển Đông
bằng mô hình ROMS
Nguyễn Minh Huấn
1,
*, Phạm Văn Sỹ
2
, Dương Hồng Sơn
2
1
Khoa Khí tượng Thủy văn và Hải dương học, Trường Đại học Khoa học Tự nhiên, ĐHQGHN,
334 Nguyễn Trãi, Hà Nội, Việt Nam
2
Viện Khoa học Khí tượng Thuỷ văn và Môi trường
Nhận ngày 11 tháng 8 năm 2010
Tóm tắt. Sự phát triển nhanh chóng của kinh tế biển và các hoạt động an ninh quốc phòng đảm
bảo chủ quyền trên biển đã đặt ra những vấn đề khoa học cấp thiết cần giải quyết đối với việc cung
cấp thông tin dự báo trường các yếu tố khí tượng thủy văn biển.
Trên thực tế, các yếu tố hải văn biển như dòng chả
y, nhiệt độ và độ muối và mực nước tổng
cộng là những yếu tố quan trọng cần thiết đáp ứng yêu cầu của các hoạt động kinh tế, an ninh quốc
phòng và nghiên cứu. Trong nghiên cứu này, nhóm tác giả đã xây dựng quy trình dự báo hạn ngắn
trường dòng chảy, nhiệt độ và độ muối bằng mô hình ROMS cho toàn bộ khu vực Biển Đông.

Mở đầu
∗
Sự phát triển nhanh chóng của kinh tế biển
và các hoạt động an ninh quốc phòng đảm bảo

chủ quyền trên biển đã đặt ra những vấn đề
khoa học cấp thiết cần giải quyết đối với việc
cung cấp thông tin dự báo trường các yếu tố khí
tượng thủy văn biển. Trong điều kiện nước ta
hiện nay, để có được các thông tin dự báo hạn
ngắn các trườ
ng yếu tố thủy văn biển có thể tiến
hành theo hai phương thức: Xây dựng hệ thống
thu nhận thông tin các trường khí tượng thủy
văn dự báo của các nước trên thế giới và khu
vực; Xây dựng, phát triển và ứng dụng hệ thống
các mô hình dự báo các trường khí tượng - thủy
văn biển. Theo phương thức thứ hai chúng ta sẽ
có được tính chủ động cao trong công tác dự
báo, các thông tin kết quả của hệ thống dự báo
nhận được sẽ đầy đủ, chi tiết và chính xác kịp
thời đáp ứng được yêu cầu của các hoạt động
kinh tế, an ninh quốc phòng và nghiên cứu.
∗
Tác giả liên hệ. ĐT: 84-4-38584943.
E-mail:
Trong các thông tin dự báo về các yếu tố
thủy văn biển trường dòng chảy, nhiệt độ và độ
mu
ối là các yếu tố quan trọng. Thông tin về
dòng chảy biển phục vụ trực tiếp cho công tác
hàng hải, tìm kiếm cứu nạn cứu hộ, dự báo lan
truyền ô nhiễm, thông tin về nhiệt độ và độ
muối góp phần vào việc xác định, dự báo ngư
trường đánh bắt cá. Trong nghiên cứu này,

nhóm tác giả đã xây dựng quy trình dự báo hạn
ngắn trường dòng chảy, nhiệt độ và độ muối
bằ
ng mô hình ROMS cho toàn bộ khu vực Biển
Đông.
N.M.Huấnvànnk./TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
363
1. Mô hình ROMS 3.0
ROMS là mô hình hoàn lưu đại dương, sử
dụng hệ phương trình nguyên thủy (primitive).
Là mô hình mã nguồn mở nên ROMS mang
tính cộng đồng rất cao, được rất nhiều các nhà
nghiên cứu sử dụng với nhiều qui mô không
gian và thời gian khác nhau: từ dải ven bờ tới
các đại dương thế giới; mô phỏng, dự báo cho
vài ngày, vài tháng và thậm chí tới hàng chục
năm. ROMS được xây dựng trên cơ sở các
nghiên cứu số trị bậ
c cao cùng với kỹ thuật tiên
tiến cho phép triển khai một cách có hiệu quả
các tính toán có độ phân giải cao. Mô hình giải
các phương trình thuỷ động lực thuỷ tĩnh và bề
mặt tự do cho các địa hình phức tạp trên hệ lưới
cong trực giao theo phương ngang và thích ứng
địa hình theo phương thẳng đứng [4].
1.1. Hệ toạ độ thích ứng địa hình theo phương
thẳng đứng
Hệ toạ độ thích ứng địa hình theo phương
thẳng đứng (lưới σ hoặc s) xấp xỉ địa hình đáy
biển và bề mặt tự do nhằm mục đích mô phỏng

các quá trình rối gần các bề mặt chất lỏng cũng
như các quá trình động lực vùng cửa sông ven
biển, mô phỏng trung thực hơn ảnh hưởng của
địa hình tới dòng chảy so với các mô hình sai
phân thông thường. Tuy nhiên, hệ toạ độ này
cũng có một số nhược
điểm nhất định, gây ra
sai số số học trong quá trình tính građien áp
suất tại các vị trí có độ dốc lớn, nhược điểm
này cũng đã được quan tâm trong nhiều năm
qua. Mặc dù các sai số này không thể loại bỏ
được hoàn toàn, nhưng có thể được giảm tới
mức có thể chấp nhận được qua phương pháp
tái tạo parabolic của Shchepetkin và
McWiliams, 2005 [2].
Hệ toạ độ thích ứng địa hình tổ
ng quát có
dạng:

(
)
()
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛

+
−
=
t,y,x)y,x(H
t,y,xz
ss
ς
ς
01 ≤≤− s
(1)
trong đó H là độ sâu, ζ là cao độ mực nước bề
mặt. Trong trường hợp s phụ thuộc tuyến tính
vào z, phép chuyển toạ độ trên sẽ trở thành hệ
toạ độ σ truyền thống. Các phép chuyển hệ toạ
độ này sẽ tạo ra hệ thống lưới không trực giao,
tuy nhiên tỷ số giữa kích thước lưới theo
phương ngang và phương thẳng đứng thường
được giả thiết là rấ
t lớn do vậy có thể bỏ qua
một số các số hạng liên quan tới hệ toạ độ cong
và phép chuyển toạ độ có thể được đơn giản
hoá thành:
zx
z
xx
z
∂
∂
∂
∂

−
∂
∂
=
∂
∂
.
σσ

Mối phụ thuộc không tuyến tính (1) cho
phép tăng cường độ phân giải của lưới tính vào
những miền có tính bất đồng nhất mạnh. Nhằm
mô phỏng tốt các quá trình động lực trong lớp
xáo trộn cũng như lớp thermocline, mô hình
ROMS sử dụng phép chuyển toạ độ dưới đây
(Haidvogel và nnk, 2000):
()(
sChhshz
ss
−
)
+
=
(2)
trong đó h
s
là độ sâu đặc trưng của lớp xáo trộn
và:
() ( )
(

)
()
()
[]
(
)
()
2/tanh2
2/tanh2/1tanh
sinh
sinh
1
θ
θ
θ
θ
θ
θ
θ
−+
+−=
ss
sC
bb
(3)
trong đó
θ, θ
b
là các các tham số không thứ
nguyên kiểm soát độ co giãn của lưới theo

phương thẳng đứng: độ phân giải tại bề mặt tỉ lệ
thuận với giá trị
θ
và tại đáy tăng lên khi
θ
b

tiến dần đến 1. Ưu điểm của phép giãn trục toạ
độ (2) và (3) là cho phép tăng độ phân giải tuỳ
ý tại những nơi độ sâu lớn (tăng độ phân giải
trong lớp nhảy vọt nhiệt độ, và lớp tà nhiệt)
trong khi duy trì độ phân giải tương đối đều tại
vùng nước nông.
N.M.Huấnvànnk./TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
364
1.2. Hệ toạ độ cong trực giao theo phương
ngang
Hệ toạ độ tính toán
ξ,η
là vuông góc và các
biên của miền tính trùng với các đường đẳng
ξ,η
. Khi hàm ánh xạ được xác định thì các hệ
số đo cũng được xác định. Các hệ số đo m và n
của hệ toạ độ cong trực giao liên kết các khoảng
cách sai phân theo hướng
ξ,η
với các cung thực
tế như sau:
()

m
ds
ξ
ξ
∂
=

()
n
ds
η
η
∂
=

(4)
do vậy, đoạn ds phải thoả mãn:
2
2
2
2
2
nm
ds
ηξ
∂
+
∂
=
(5)

hay trong toạ độ Đề các:
()(
()()
()
ηξηξ
ηξηξ
ηξηξηηξξ
ηξηξ
ddyyxxdyxdyx
dydydxdxdydxds
+++++=
+++=+=
2
222222
22
222
)
(6)
Với thành phần cuối cùng ở phương trình
trên bằng 0 trong hệ toạ độ trực giao, phép thế
hai phương trình trên dẫn đến (Wilkin and
Hedstrom, 1998):
()
()
2/1
22
2/1
22
−
−

+=
+=
ηη
ξξ
yxn
yxm
(7)
1.3. Hệ phương trình thuỷ động lực
Phương trình trạng thái
()
PTS ,,
ρ
ρ
=
và
giả thiết thuỷ tĩnh
g
z
P
ρ
=
∂
∂

Phương trình chuyển động
()
uu
zz
z
zzzz

DF
mn
H
g
zgP
n
H
vH
m
u
n
v
mn
f
mn
uH
sm
uvH
n
uH
mn
uH
t
++
⎟
⎟
⎠
⎞
⎜
⎜

⎝
⎛
∂
∂
+
∂
∂
+
∂
∂
−=
⎭
⎬
⎫
⎩
⎨
⎧
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
−
⎟
⎠
⎞
⎜

⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
−
⎟
⎠
⎞
⎜
⎝
⎛
Ω
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂

+
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
ξ
ς
ξρ
ρ
ξηξ
ηξ
0
2
11

(8)
()
vv
zz
z
zzzz
DF
mn
H
g
zgP
m
H
uH
m
u
n
v
mn
f
mn
vH
sm
vH
n
uvH
mn
vH
t
++

⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
∂
∂
+
∂
∂
+
∂
∂
−=
⎭
⎬
⎫
⎩
⎨
⎧
⎟
⎠
⎞
⎜
⎝
⎛
∂

∂
−
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
+
⎟
⎠
⎞
⎜
⎝
⎛
Ω
∂
∂
+
⎟
⎟

⎠
⎞
⎜
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
η
ς
ηρ
ρ

ηηξ
ηξ
0
2
11

(9)
Phương trình bảo toàn nhiệt độ và độ muối
()
TT
zzzzz
DF
mn
H
mn
TH
sm
vTH
n
uTH
mn
TH
t
+=
⎟
⎠
⎞
⎜
⎝
⎛

Ω
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂

ηξ
(10)
()
SS
zzzzz
DF
mn
H
mn
SH
sm
vSH
n
uSH
mn
SH
t
+=
⎟
⎠
⎞
⎜
⎝
⎛
Ω
∂
∂
+
⎟
⎠

⎞
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
+
⎟
⎠
⎞
⎜
⎝
⎛
∂
∂
ηξ
(11)
trong đó D
u
, D
v

, D
S
, D
T
là các thành phần nhớt
và khuếch tán rối đối với các biến động lượng
và nhiệt muối; F
u
, F
v
, F
S
, F
T
là các thành phần
ngoại lực và nguồn.
() ()
⎥
⎦
⎤
⎢
⎣
⎡
∂
∂
−
∂
∂
−
∂

∂
+−=
∂
∂
=Ω
ηξ
ηξ
z
nv
z
mu
t
s
sw
Ht
s
ts
z
1
1
,,,

(12)
() ()
ξξ
ς
ξ
∂
∂
+

∂
∂
+=
∂
∂
h
sCs
z
1
(13)
() ()
ηη
ς
η
∂
∂
+
∂
∂
+=
∂
∂
h
sCs
z
1
(14)
Với giả thiết áp suất thuỷ tĩnh cho phép tính
được thành phần vận tốc theo phương thẳng
đứng từ phương trình liên tục do đó các mô

hình loại này đôi khi còn được gọi là “tựa ba
chiều”, để phân biệt với các mô hình sử dụng
đầy đủ cả ba phương trình động lượng.
1.4. Sai phân hoá
Sai phân theo không gian

Sai phân trung tâm bậc hai trên lưới
Arakawa C áp dụng cho phương ngang (
ξ,η
)
với các điều kiện biên trượt tự do (toàn phần),
trượt một phần, hoặc điều kiện dính. Theo
phương thẳng đứng (trục s) sử dụng sai phân
xen kẽ bậc hai.

N.M.Huấnvànnk./TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
365
Cũng như các mô hình khác, tuy đã sử dụng
phương pháp giảm thiểu sai số gây ra bởi thành
phần gradient áp suất, ROMS sử dụng kỹ thuật
là trơn địa hình đáy để tăng độ ổn định của mô
hình như sau:

2/12/1
2/12/1
2
++
−+
+
−

=
∆
=
H
H
hh
h
h
r
(15)
Vì độ dày của các lớp nước thường nhỏ hơn
rất nhiều so với bước tính theo phương ngang
nên bước thời gian thường bị giới hạn bởi điều
kiện ổn định theo phương thẳng đứng
(Vreugdenhill 1994),
∆
t <
∆
z
2
/4Nv, hơn là
theo phương nằm ngang (chỉ tiêu Courant-
Friedrichs-Levy).
Sai phân theo thời gian
Do vậy, có hai cách sử dụng mô hình,
phương pháp thứ nhất là giải phương trình
chuyển động để tính vận tốc ngang, sau đó tích
phân phương trình liên tục từ mặt tới đáy để
tính thành phần vận tốc thẳng đứng và áp dụng
điều kiện biên động học tại mặt nước để tính sự

thay đổi của mực n
ước. Phương pháp thứ hai
hay kỹ thuật tách thời gian được sử dụng rộng
rãi trong các mô hình hoàn lưu đại dương nhằm
phân tách các sóng chính áp có tốc độ chuyển
động nhanh mà vẫn tiết kiệm thời gian tính:
Chế độ chính áp: giải hệ phương trình động
lượng hai chiều “trung bình theo độ sâu” với
bước thời gian ngắn, thoả mãn điều kiện
Courant-Friedrichs-Levy
Chế độ tà áp: giải hệ phương trình động
lượng ba chi
ều với bước thời gian tương đối
dài.
Tuy nhiên, sai số khác nhau giữa hai hệ
phương trình này đòi hỏi các phương pháp điều
chỉnh sao cho chúng cùng thỏa mãn phương
trình liên tục và bảo toàn các đại lượng vô
hướng.
2. Quy trình dự báo
Quy trình dự báo trường dòng chảy, nhiệt
độ, độ muối được viết trên mã nguồn mở
(UNIX) kết hợp với ngôn ngữ lập trình
MATLAB, cho phép thực hiện dự báo nghiệp
vụ
tự động, đạt hiệu quả cao (Hình 1).
Quy trình dự báo được chia làm 3 giai đoạn:
1. Giai đoạn tiền xử lý:
Download số liệu khí tượng, điều kiện ban
đầu và điều kiện biên… từ nguồn số liệu ECCO

dự báo đến 10 ngày.
Nội suy các yếu tố về miền lưới tính toán
Kiểm tra và xử lý lại kết quả đầu vào
2. Giai đoạn chạy mô hình
Ch
ạy mô hình ROMS với các thông số và
số liệu đã được ấn định
3. Giai đoạn xử lý kết quả
Trích xuất các yếu tố (dòng chảy, nhiệt độ,
độ muối) để tạo thành các file có dung lượng
nhỏ dạng netcdf
Thể hiện các trường dòng chảy, nhiệt độ, độ
muối dưới dạng ảnh và định dạng kml (định
dạng hiển thị trên google earth hay google map)
Xuất các yế
u tố (dạng file netcdf, dạng ảnh
và kml) lên bản tin dự báo biển
www.dubaobien.vn).

Hình 1. Sơ đồ quy trình dự báo dòng chảy, nhiệt
độ, độ muốivà mực nước tổng cộng bằng mô hình
ROMS cho khu vực Biển Đông.
N.M.Huấnvànnk./TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
366
Quy trình dự báo được thực hiện chạy tự
động thông một file script viết bằng trình thông
lệnh shell trên hệ điều hành Linux. Một số công
đoạn của quy trình đã được cải tiến nhằm phù
hợp với các điều kiện của Việt Nam bởi Trung
tâm Nghiên cứu Môi trường, Viện Khoa học

Khí tượng Thuỷ văn và Môi trường và Trường
Đại học Khoa học Tự nhiên trong khuôn khổ
của
đề tài cấp Nhà nước: “Nghiên cứu phát
triển và ứng dụng công nghệ dự báo hạn ngắn
trường các yếu tố thủy văn biển khu vực Biển
Đông” mã số KC.09.16/06-10.
Các kết quả đáng tin cậy của mô hình dự
báo đã được hoàn toàn khẳng định trong [1] nên
trong bài báo này không trình bày lại các phần
hiệu chỉnh và kiểm nghiệm mô hình.
3. Các điều kiện đầu vào của mô hình
3.1. Miền dự báo
Địa hình đáy biể
n được nội suy từ số liệu
địa hình đáy biển toàn cầu ETOPO-2 có độ
phân giải 2’ (khoảng 4km). Độ sâu lớn nhất của
miền tính khoảng 4500 m và độ sâu nhỏ nhất
được giới hạn bằng 10 m và được chia thành 12
lớp theo phưong thẳng đứng. Độ sâu nhỏ nhất
được giới hạn bằng 10 m nhằm tăng bước thời
gian giới hạn bởi hiệu ứng khuyếch tán rối theo
phương thẳng đứng
∆
t <
∆
z
2
/4Nv, quá trình
này chủ yếu xảy ra ở vùng nước nông gần bờ

(Vreugdenhill 1994).Với kính thước bước lưới
như trên, bước thời gian cho chế độ barocline là
1800 giây, và barotrop là 60 giây (Hình 2).
3.2. Điều kiện ban đầu, điều kiện biên và ngoại
lực
Số liệu đầu vào dự báo 10 ngày được cung
cấp từ mô hình hoàn lưu toàn cầu và số liệu
quan trắc của ECCO (dự án “Dự báo Hoàn lưu
và Khí tượng Đại dương”) do cơ
quan Nghiên
cứu Hải quân, Mỹ (ONR) cùng với rất nhiều
các tổ chức, cơ quan (NOAA, NASA…) đồng
thực hiện.

Hình 2. Bản đồ địa hình khu vực Biển Đông


Hình 3. Phân bố theo phương ngang và thẳng đứng
của trường nhiệt độ nước biển trên Biển Đông.

N.M.Huấnvànnk./TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
367
Lưu lượng của các sông lớn được lấy trung
bình cho hai mùa chính (mùa khô và mùa mưa):
Lưu lượng nước sông Mê Kông dao động
trong khoảng 37000 m
3
/s (tháng I) - 148000
m
3

/s (tháng VIII) và đổ ra biển tại hai cửa có toạ
độ (9.421
o
N;106.289
o
E) và (10.171
o
N;106.80
o
E)
Lưu lượng nước sông Hồng dao động trong
khoảng 3700 m
3
/s (tháng I) - 10000 m
3
/s (tháng
VIII) và đổ ra biển tại hai cửa có toạ độ
(20.385
o
N;106.643
o
E) và (20.654
o
N;106.914
o
E)
Trên cả hai hệ thống sông nhiệt độ dao
động trong khoảng 20
o
C-25

o
C và độ muối luôn
bằng 0 %
0.

Hình 4. Đẳng biên độ và đẳng pha sóng M2.
Các tham số của 6 sóng thuỷ triều, M2, S2,
N2, K2, K1, O1, được đưa vào các biên lỏng
của miền nghiên cứu. Các tham số này được lấy
từ mô hình thuỷ triều toàn cầu TPXO có độ
phân giải 1/4 độ kinh vĩ. Mô hình này sử dụng
phương pháp bình phương tối thiểu để xấp xỉ
phương trình thuỷ triều Laplace với số liệu cao
độ bề mặt biển thu được dọc theo 324 quỹ đạo
của vệ tinh TOPEX/Poseidon
4. Xu
ất bản tin dự báo
4.1. File dữ liệu
Để thuận tiện cho người dùng khi download
số liệu cũng như làm giảm dung lượng của file,
các file đầu ra được trích xuất cho riêng các yếu
tố: dòng chảy, nhiệt độ và độ muối trong khu
vực tính toán dưới định dạng netcdf (Hình 5).
Số liệu sau khi dowload, có thể sử dụng các
phần mềm hỗ trợ như ncbrowser để chuyển đổi
sang định dạng ASCII theo khu vự
c phù hợp
với yêu cầu của người sử dụng.

Hình 5. Giao diện download số liệu dự báo của các

yếu tố dòng chảy, nhiệt độ, độ muối và mực nước
tổng cộng khu vực Biển Đông.
4.2. Định dạng ảnh
Hình ảnh bản đồ phân bố các yếu tố: dòng
chảy, nhiệt độ và độ muối cho toàn khu vực
Biển Đông được chia thành 9 vùng nhỏ. Người
truy cập có thể phóng to các vùng bằng cách
kích con trỏ chuột lên từng khu vực (Hình 6).
N.M.Huấnvànnk./TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
368


Hình 6. Bản đồ dự báo phân bố trường vectơ dòng
chảy trên nền trường độ muối định dạng ảnh khu
vực Biển Đông và các vùng lựa chọn.
4.3. Định dạng kml – Google earth
Kết quả cũng có thể được thể hiện trên
Google earth và Google map bằng các file định
dạng klm để người dùng có thể sử dụng khả
năng có sẵn của Google earth và Google map
phóng to, thu nhỏ tuỳ ý các khu vực quan tâm .

Hình 7. Bản đồ dự báo phân bố trường dòng chảy
trên nền trường độ muối với định dạng kml khu vực
Biển Đông.
Kết luận
Quy trình dự báo trường dòng chảy, nhiệt
độ, độ muối và mực nước tổng cộng được xây
dựng trên cơ sở kế thừa các kết quả của cơ quan
Nghiên cứu Hải quân Mỹ (ONR) kết hợp với

các tổ chức, cơ quan khác như NOAA,
NASA…, là một quy trình khép kín được thực
thi tự động với nhiều các công đoạn phức tạp:
download số liệu d
ự báo toàn cầu, xử lý số liệu,
thực hiện chạy mô hình dự báo, xử lý kết quả
và xuất bản tin dự báo.
Kết quả của quy trình cho phép nhận được
thông tin dự báo toàn cảnh chi tiết trước nhiều
ngày về trường dòng chảy, nhiệt độ, độ muối và
mực nước tổng cộng cho khu vực Biển Đông,
các thông tin của hệ thống dự báo sẽ được cung
cấp đầ
y đủ, chi tiết và chính xác kịp thời đáp
ứng được yêu cầu của các hoạt động kinh tế,
hàng hải, tìm kiếm cứu nạn cứu hộ, an ninh
quốc phòng và nghiên cứu.

N.M.Huấnvànnk./TạpchíKhoahọcĐHQGHN,KhoahọcTựnhiênvàCôngnghệ26,Số3S(2010)362‐369
369
Tài liệu tham khảo
[1] Dương Hồng Sơn, Nguyễn Tài Hợi, “Dự báo
hạn ngắn trường dòng chảy và nhiệt độ Biển
Đông”, Đề tài Xây dựng mô hình dự báo các
trường khí tượng thuỷ văn vùng Biển Đông –
KC.09.04.
[2] A.F. Shchepetkin, J.C. McWiliams, The regional
oceanic modeling system (ROMS): a split-
explicit, free-surface, topography-following-
coordinate oceanic model. Ocean Model 9

(2005) 347-404.
[3]

[4]



Experimental forecast system for current, salinity, temperature
and water level fields in the South China Sea
Nguyen Minh Huan
1
, Pham Van Sy
2
, Duong Hong Son
2
1
Faculty of Hydro-Meteorology & Oceanography, Hanoi University of Science, VNU,
334 Nguyen Trai, Hanoi, Vietnam
2
Institute of Hydrology Meteorology and Environment

The rapid development of marine economy and national defense and security activities to ensure
maritime sovereignty posed scientific problem solving urgent need for providing forecast
meteorological and oceanographical information in the sea. In the country's current conditions, to
obtain the short-term forecasts of oceanographical factors in the sea can proceed in two ways: first
approach is build a system receiving information of meteorology and oceanography forecast from

countries and regions worldwide, second one is development and application of system models to
forecast the marine weather and hydrodynamics. By the second approach, we would have highly sense
in initiative in the marine forecast, the forecasting system will give full and accurate information and

timely to meet the requirements of economic , security, defense and research activities.
In the fact, oceanographical factors in sea such as currents, temperature, salinity and water level
are important ones, information about flows, water level and waves directly serve for high seas
navigation, the maritime search and rescue works, spreading pollution forecasting, information about
temperature and salinity contribute to the identification, forecasting fishing grounds. In this study, the
authors developed a short-term forecasting system for flows field, temperature, salinity and water level
by Regional Ocean Model System (ROMS) for the South China Sea.