NGHIÊN CỨU DIỄN BIẾN BỒI XÓI VÙNG VEN BIỂN CẦN GIỜ BẰNG MÔ HÌNH TOÁN

NGHIÊN CỨU DIỄN BIẾN BỒI XĨI VÙNG VEN BIỂN CẦN GIỜ
BẰNG MƠ HÌNH TỐN
STUDY ON THE TENDENCY OF ACCRETION AND EROSION IN THE
CAN GIO COASTAL ZONE ON THE BASE MATHEMATICS MODEL
TS. Nguyễn Thị Bảy
(1)
, KS. Trần Thiện Tồn
(2)
(1) Ðại học Bách Khoa, Ðại Học Quốc Gia Tp Hồ Chí
Minh.
(2) Ðại học KHTN, Ðại Học Quốc Gia Tp Hồ Chí Minh.
TĨM TẮT
Bài báo này, các tác giả trình bày nghiên cứu về tính tốn dòng chảy, phù sa lơ lửng và khả
năng diễn biến bồi xói khu vực vùng ven biển Cần Giờ dưới ảnh hưởng của thủy triều và gió.
Cơ sở lý thuyết của bài tốn là hệ phương trình Reynolds, kết hợp với phương trình chuyển
tải bùn cát, lấy trung bình theo chiều sâu, có tính đến hàm số nguồn, mơ tả tốc độ bốc lên
hay lắng xuống của hạt, và phương trình liên tục bùn cát đáy. Các kết quả về nồng độ phù sa
lơ lửng và bồi xói đáy khu khu vực tương đối phù hợp so với thực đo.
ABSTRACT
In this paper, the authors present the studies on the numerical models to simulate the
currents, sediment transport by tide and wind. The results from the models are used to
estimate the tendency of accretion and erosion in the Can Gio coastal zone. The basic theory
in the numerical models is Reynolds equation system combined with the sediment transport
equation and the conservation equation for sediment mass. The computational results of the
mass concentration of sediment and the tendency of accretion and erosion in study area are
in agreement with the satellite photos.
I. VÙNG BIỂN CẦN GIỜ VÀ
VỊ TRÍ NGHIÊN CỨU
Cần Giờ được giới hạn trong tọa
độ đòa lý ở 10
0

15’- 10
0
40’N,
106
0
45’- 107
0
00’E, nằm trong
vùng cận xích đạo, thuộc vùng
cửa sông Đồng Nai. Đây là
vùng đất thấp, bằng phẳng, bò
ngập mặn quanh năm hoặc một
số thời gian trong năm. Cần Giờ
vừa mang đặc trưng khí hậu
biển vừa pha khí hậu lục đòa
Đông Nam Bộ (xem Hình 1).
Vùng biển Cần Giờ bao gồm:
vùng biển trước các cửa sông, vònh Gành Rái, vònh Đồng Tranh và vùng bãi triều Cần
Giờ.
Vùng biển trước cửa sông có bờ biển chạy theo hướng Đông Bắc – Tây Nam, chia làm
hai phần: từ Vũng Tàu lên Hàm Tân, phía Tây Nam từ Vũng Tàu đến Gò Công. Cửa
sông ở đây nông dần xuống phía Nam do ảnh hưởng bồi đắp cát từ đất liền.
Vònh Gành Rái ăn sâu vào đất liền, phía Đông giáp Vũng Tàu, phía Tây là Cần Giờ và
vùng bãi cạn, phía Nam là biển Đông, phía Bắc giáp đảo Long Sơn. Đổ nước vào vònh là
1
Hình1 : Vị trí vùng biển Cần Giờ và khu vực tính tốn
ba con sông lớn: sông Ngã Bảy, sông Thò Vải và sông Dinh. Đường bờ bao quanh vònh
khúc khuỷu và dốc.
Vònh Đồng Tranh, đổ vào vùng này là sông Soài Rạp và sông Đồng Tranh. Nhìn chung
đòa hình toàn vùng có hướng dốc từ Bắc xuống Nam, theo hướng các dòng sông và hướng

dốc từ Tây sang Đông, từ bờ ra biển. Đường bờ tương đối đơn giản, thoải, phần lớn là các
bãi bồi.
Dưới đây là một vài giá trò hàm lượng phù sa lơ lửng di đẩy và các cấp hạt đã được đo
đạc thực tế trong vùng. Đây là cơ sở để đưa các giá trò đầu vào, giá trò biên vào mô hình
tính toán. [8]
Phù sa lơ lửng
•
Sườn triều lên
- Tại vònh Gành Rái: Hàm lượng phù sa lơ lửng lớn nhất là 32.4mg/L và nhỏ
nhất là 16.9mg/L. Loại cấp hạt từ <0.038 ÷ 0.10mm chiếm tỉ lệ cao: 93 ÷ 98%, từ 0.1 ÷
1.0mm chiếm tỉ lệ thấp.
- Tại cửa Soài Rạp: Hàm lượng phù sa lơ lửng lớn nhất là 58.1mg/L và nhỏ nhất
là 17.3mg/L. Loại cấp hạt <0.038 ÷ 0.10mm chiếm tỉ lệ: 85 – 95%, từ 0.1 ÷ 1.0mm chiếm
tỉ lệ thấp.
•
Sườn triều xuống
- Tại vònh Gành Rái: Hàm lượng phù sa lơ lửng lớn nhất là 40.0mg/L và nhỏ
nhất là: 24.8mg/L. Loại cấp hạt <0.038 ÷ 0.1mm chiếm tỉ lệ 88 ÷ 93%, loại cấp hạt từ
>0.1 ÷ 1.0mm chiếm tỉ lệ thấp.
- Tại cửa Soài Rạp: Hàm lượng phù sa lơ lửng lớn nhất là 93.3mg/L và nhỏ nhất
là 10.8mg/L. Loại cấp hạt từ 0.038 ÷ 0.1mm chiếm tỉ lệ cao 93 – 97%, loại cấp hạt >0.1
÷ 1.0mm chiếm tỉ lệ thấp.
Phù sa di đẩy
•
Sườn triều lên:
- Tại vònh Gành Rái: Hàm lượng phù sa di đẩy lớn nhất là 1398.1mg/L và nhỏ
nhất là 225.9mg/L. Loại cấp hạt <0.038 ÷ 0.25mm chiếm tỉ lệ 98%, loại cấp hạt từ 0.5 ÷
1.0mm chiếm tỉ lệ 10%.
- Tại cửa Soài Rạp: Hàm lượng phù sa di đẩy lớn nhất là 1266.4mg/L và nhỏ nhất
là 386.7mg/L. Loại cấp hạt từ <0.038 ÷ 0.10mm chiếm tỉ lệ 6%.

•
Sườn triều xuống:
- Tại vònh Gành Rái: Hàm lượng phù sa di đẩy lớn nhất là 997.6mg/L và nhỏ
nhất là 91.5mg/L. Loại cấp hạt <0.038 ÷ 0.50mm chiếm tỉ lệ >90%, loại cấp hạt từ 0.5 ÷
1.0mm chiếm tỉ lệ <10%.
- Tại cửa Soài Rạp: Hàm lượng phù sa di đẩy lớn nhất là 2265mg/L và nhỏ nhất
là 106.4mm/L. Loại cấp hạt từ <0.038 ÷ 0.50mm chiếm tỉ lệ < 10%.
II. CƠ SỞ LÝ THUYẾT
Các phương trình cơ bản là các phương trình đã lấy trung bình theo chiều sâu như sau:
1. Hệ phương trình Reynolds
uA
h
vu
Ku
hx
gfv
y
u
v
x
u
u
t
u
sx
2
22
)(
)(
∇+

+
+
−
+
+
∂
−∂
−=−
∂
∂
+
∂
∂
+
∂
∂
ςςρ
τ
ςς
(1)
vA
h
vu
Kv
hy
gfu
y
v
v
x

v
u
t
v
sy
2
22
)(
)(
∇+
+
+
−
+
+
∂
−∂
−=+
∂
∂
+
∂
∂
+
∂
∂
ςςρ
τ
ςς
(2)

2. Phương trình liên tục:
2
[ ] [ ]
0
)()(
=
∂
+∂
+
∂
+∂
+
∂
∂
y
vh
x
uh
t
ςςς
(3)
3. Phương trình chuyển tải phù sa lơ lủng
Việc nghiên cứu sự chuyển tải phù sa, bùn cát dựa trên những giả thiết sau đây:
Hạt vật chất chia làm hai dạng: bùn cát đáy và bùn cát lơ lửng. Bùn cát đáy nằm
trong lớp đáy có bề dày a, ở độ cao z trong khoảng: - h < z < - h + a. Lớp bùn cát lơ lửng
nằm ở phía trên với cao độ: - h + a < z <ς. Theo thực nghiệm, giá trò “a” nhỏ nhất có thể
chấp nhận: a
min
= 0,01H (Vanrijn, 1986).
Sự phân bố bùn cát lơ lửng theo không gian và thời gian được mô tả bởi phương

trình hai chiều tổng quát sau:
H
S
y
y
HD
yHx
x
HD
xHyx
v
t
+
∂
∂
∂
∂
+
∂
∂
∂
∂
=









∂
∂
+
∂
∂
+
∂
∂
)(
1
)(
1 CCC
v
C
u
C
γ
(4)
Trong đó:
h: độ sâu tính từ mặt thoáng yên lặng đến đáy của từng điểm trong miền tính.
K: hệ số ma sát đáy:
ς : cao độ mực nước tại từng điểm trong miền tính
H = ς + h
u, v, _ Thành phần vận tốc trung bình dòng chảy theo phương x, y.
S
τ
: ứng suất tiếp gió, tính theo công thức thực nghiệm [18]:
saS
WWρCτ

10
=
, đơn vò N/m
2
3
10
10).067.075.0(
−
+=
WC
,
22
yx
WWW +=
,
W : vận tốc gió (m/s)
Ws = Wx, W
y
thành phần vận tốc gió theo các phương O
x
, O
y
a
ρ
- khối lượng riêng của không khí trên mặt thoáng [kg/m
3
].
C : Nồng độ phù sa lơ lửng trung bình theo chiều sâu.
γ
v

: hệ số phân bố vận tốc và nồng độ theo chiều sâu














−−















−−
=
−Z
v
H
a
Z
H
a
Z
γ
1
2.0
1).2.1(
1).1(
Với Z : Thông số lơ lửng, [3]:
)2(
* s
s
ωuχ
ω
Z
+
=
:
χ: Hằng số Karman.
u
*
: Vận tốc động lực của dòng chảy,
)(

22
2
2
*
vu
C
g
u +=
;
C : Hệ số Chezy,








=
s
K
H
C 12lg18
K
s
: Độ nhám tuyệt đối ở đáy, K
s
= 3D
90
, D

90
: đường kính hạt tương ứng với tần
suất 90%, [3]
D
x
, D
y
,

_ Hệ số truyền chất theo phương x, y, tính thực nghiệm của Elder

[11]:
3
1
93.5
−
=
CuhgD
x
;
1
93.5
−
=
CvhgD
y
S : Hàm số nguồn, mô tả sự bốc lên hay lắng xuống của hạt, tính thực nghiệm theo
Van Rijin [3].
S = E đối với τ
b

> τ
e
S = -D đối với τ
b
< τ
d
S = 0 đối với τ
d
≤ τ
b
≤ τ
e
Với:
E : Tốc độ xói.
D : Tốc độ bồi lắng.
τ
b
: Ứng suất tiếp đáy:
2
b
1
= kV
2
τ ρ
τ
d
: Ứng suất tiếp đáy tới hạn để các hạt lắng xuống.
τ
e
: Ứng suất tiếp đáy tới hạn để các hạt bốc lên.

ρ : khối lượng riêng nước biển, k là hệ số ma sát đáy ;
2 2
V u v= +
.
Hàm số E, D được tính thực nghiệm theo Van Rijn [3]:
e
eb
τ
ττ
ME
−
=

Với M=0,00001 kg/(sm
2
)
b
s,m
d
D 1 cw
 
τ
= −
 ÷
τ
 
Trong đó w
s,m
=(1- C)
γ

w
s
là độ thơ thuỷ lực của hạt trong hỗn hợp nước biển-bùn cát.
w
s
: độ thơ thuỷ lực trong nước trong; w
s
=(s-1)gd
2
/18ν
γ : hệ số dao động trong khoảng 2,3-4,6 tuỳ theo số Reynolds. Đối với hạt có đường
kính 0,05-0, mm, có thể chấp nhận γ = 4, [3].
S : Tỷ số giữa trọng lượng riêng của hạt và của nước.
3. Phương trình liên tục bùn cát đáy






∂
∂
+
∂
∂
+









∂
∂
∂
∂
+






∂
∂
∂
∂
+
−
=
∂
∂
y
q
x
q
y
C

HK
yx
C
HK
x
E
t
h
by
bx
yx
p
ε
1
1
Dòng bùn cát đáy q
b
= (q
bx
, q
by
) được xác đònh theo công thức [3]:
22
3.01.25.15.0
),(
))1((053.0
vu
vu
DTdgSq
mb

+
−=
−
∗
ε
p
: độ rỗng của lớp bùn cát đáy
D
*
: Thông số đường kính hạt.:
( )
3
1
2
1






−
=
∗
ν
Sg
dD
m
ν : Hệ số động học nhớt.
T : Thông số khả năng chuyển tải của hạt:







−
=
∗
∗∗
2
cr,
2
cr,
2
u
uu
T

cr,
u
∗
: Vận tốc tới hạn của dòng chảy, nếu lớn hơn vận tốc này, hạt sẽ bò dòng chảy bốc
lên hòa vào chuyển động của dòng, theo Shills :
15
1
15
8
15
9

15
8
cr,
)1(25.0u
−
∗
−=
νgdS
d : Đường kính trung bình của hạt; g: gia tốctrọng trường
III. PHƯƠNG PHÁP TÍNH TOÁN
4
Các phương trình cơ bản trên được giải bằng phương pháp sai phân hữu hạn sơ đồ
ADI, sai phân luôn hướng theo mỗi nửa bước thời gian.
Sau khi rời rạc hoá các phương trình vi phân trên từng nút của lưới sai phân, ta biến
đổi về hệ phương trình đại số tuyến tính dạng 3 đường chéo và sử dụng giải thuật truy
đuổi để giải ra nghiệm số u, z, v, C cho từng nút được quy đònh như trong lưới sai phân
(xem Hình 2)
1. Lưới sai phân:
2. Xử lý biên
a) Điều kiện ban đầu
Ở thời điểm ban đầu, cho vận tốc bằng 0, mực nước bằng mực nước tónh, nồng độ phù sa
bằng nồng độ nền thật nhỏ hoặc bằng 0
b) Điều kiện biên:
• Đối với bài toán tính dòng chảy, áp dụng điều kiện biên không thấm (u
n
=0) ở bờ
và dao động mực nước triều : ζ=f(t) ở biên lỏng.
• Đối với bài toán chuyển tải:
- Biên bờ:
0=

∂
∂
s
C
- Biên lỏng: Khi nước từ ngoài vào miền tính: C=C
0
với C
0
là nồng độ biên cho
trước. Khi nước chảy từ miền tính ra:
0
2
2
=
∂
∂
S
C
c) Tính toán đường bờ động
Việc tính toán các yếu tố thủy triều còn phụ thuộc vào lúc thủy triều lên hay
xuống.
Khi thủy triều xuống, làm cho diện tích mặt nước giảm xuống, những vùng cạn có
độ sâu thấp sẽ hiện ra và không bò ngập nước. Khi đó, ta sẽ không tính tại những vùng đo
(các ô lưới ở vùng đó được đóng lại. Giới hạn biên để tính toán sẽ được dời theo mực
nước rút của thủy triều, nghóa là mực nước thủy triều càng xuống thấp làm cho các vùng
cạn trơ ra càng nhiều thì giới hạn của vùng tính toán càng nhỏ lại.
Khi thủy triều lên, những vùng cạn này lại dần dần ngập trong nước, các ô lưới tính
toán của các vùng này lại được mở ra để tính. Giới hạn biên tính toán cũng sẽ lớn dần
theo mực nước lên của thủy triều. Khi những vùng cạn mới bắt đầu ngập nước, những
vùng này sẽ là những vùng biên mới thì ta cần có giá trò mực nước tại vùng này, ta sẽ

5
Chú thích:
Vận tốc u
Vận tốc v
, C
Độ sâu h
1/2 i=2 3 9/2
j=2
3
x
y
1
Hình 2 - Lưới sai phân
gán giá trò mực nước và nồng độ của nó bằng với giá trò mực nước và nồng độ của vùng
kế bên liền với nó.
IV. CÁC THÔNG SỐ SỬ DỤNG TRONG TÍNH TOÁN
Lưới tính tốn cho vùng biển
Cần Giờ bao gồm 280 ơ theo
chiều ngang và 116 ơ theo
chiều thẳng đứng, với
∆x=∆y=150m; tham số
Coriolis f=2ωsinϕ, với
ϕ=10
0
, suy ra f=2,53.10
-5
s
-1
hệ số rối ngang A=100m
2

/s,
hệ số ma sát đáy được tính
theo cơng thức:
2
2 1/ 3
g gn
k
C H
= =

Với n là hệ số nhám đáy.
Độ sâu h được lấy từ bản đồ
phân bố độ sâu (xem Hình
3), [7].
Các điều kiện biên về dòng
chảy:
Trên các biên lỏng, dao dộng mực nước được cho dưới dạng tổng hợp 11 sóng triều. Số
liệu về các sóng được lấy từ [6],[7],[8]. Sau đây là một vài bảng số liệu:
Các điều kiện về nồng độ
Nồng độ ban đầu trên toàn miền tính: 0.015 kg/m
3
Nồng độ biên lỏng: Khi dòng chảy hướng vào miền tính, ta cho nồng độ biên ở các biên
lỏng như sau [5],[7]:
- Cửa biển: 0.02 kg/m
3

- Đồng Tranh: 2.418 - 0.052 kg/m
3
- Ngã Bảy: 0.08 – 0.02 kg/m
3

- Cái Mép: 0.08 – 0.02 kg/m
3
- Soài Rạp: 0.3 – 0.075 kg/m
3
6
Bảng 1: Đặc điểm sóng tại trạm Sồi Rạp
STTT ên sóngBiên độ (m)Tần số (s
-1
)Pha ban
đầu
(rad)1M20.80.0001451.32N20.0890.0001450.96
323S20.310.0001451.67534K20.1530.0001451.
67535K10.6210.000070285.67516O10.4510.000
070285.67517P10.20.000070285.17518Q10.092
0.000070286.42019M40.0150.0002814.256710
M60.0020.0004211.277111MS60.0130.0002815
.2235
Bảng 2: Đặc điểm sóng tại trạm Đồng
Tranh
STTT ên sóngBiên độ (m)Tần số (s
-1
)Pha ban
đầu
(rad)1M20.780.0001451.392N20.0870.000145
1.05323S20.3120.0001451.76534K20.1530.00
01451.76535K10.6210.00007025.31356O10.4
510.00007025.76517P10.20.00007025.26518Q
10.0990.00007026.51019M40.0150.0002814.3
46610M60.0020.0004211.367111MS60.0130.
0002815.7651

Cần Giờ
Vũng Tàu
Hình 3 : Giới thiệu địa hình khu vực nghiên cứu
m
Các thông số khác:
Dựa vào thành phần các chất hữu cơ từ các sơng đổ ra biển Cần Giờ [4], theo tra cứu của
Van-Rijn [3] và tham khảo thêm trong [10], ta chọn:
τ
d
= 0.06 N/m
2
; τ
e
= 0.15 N/m
2
; M=0,00001 kg/sm
2
;
ν = 1.011E-6 m
2
/s
Hệ số γ trong công thức : w
s,m
=(1- C)
γ
w
s
được chọn = 4
Khối lượng riêng trung bình của hạt: ρ
s

= 2650 kg/m
3
Khối lượng riêng của nước biển: ρ
s
= 1025 kg/m
3
Thời gian tính
Mơ hình mơ phỏng và tính tốn sự thay đổi lòng dẫn trong vòng 3 tháng với bước thời gian
st 9=∆
.
V. KẾT QUẢ TÍNH TỐN
Các kết quả tính toán về mực nước, vận tốc nồng độ phù sa lơ lửng được trình bày
trong các Hình 4, 5, 6 và 7.
a) Phân tích mực nước triều
Kết quả tính toán (xem Hình 4, 5) cho thấy khi triều rút cực đại, dọc ven bờ biển Cần
Giờ mực nước rút ra xa trung bình khoảng 3 km, những khu vực mà mực nước rút xa nhất
là ở bãi biển Ngọc Điệp 5,5 km, bãi bồi ở khu vực ngoài cửa sông Soài Rạp (về phía Tây
Nam của vùng tính) 5 km. Ở bãi biển 30/4 khoảng 1,5 km.
Phân tích các kết quả tính toán cho ta thấy sự dao động đường bờ lúc nước lớn và
nước ròng trong một phạm vi khá lớn, điều này ảnh có hưởng nhiều đến tính toán dòng
chảy, phù sa lơ lửng và bồi xói lòng dẫn cho vùng này.
b) Phân tích phù sa lơ lửng
Nhìn chung nồng độ phù sa tập trung nhiều ở nhánh Đồng Tranh và Soài Rạp, có giá
trò trong khoảng 0.1 – 0.2 kg/m
3
. Trong khi đó nồng độ ở mũi Vũng Tàu tương đối thấp từ
0.025 – 0.05 kg/m
3
. Điều này được lý giải do lượng phù sa từ thượng lưu đổ vào vùng tính
Đồng Tranh và Soài Rạp cao, còn ở mũi Vũng Tàu chòu ảnh hưởng mạnh từ biên biển.

c) Tình hình diễn biến đáy biển.
Kết quả về tính toán khả năng bồi xói trong khu vực được trình bày trong Hình 8.
Sự bồi lắng và xói của vùng ven biển Cần Giờ là sự kết hợp của nhiều nguyên nhân
như: ảnh hưởng của gió, ảnh hưởng của sóng biển, dòng chảy, dòng chảy ven bờ, đòa
hình đáy, mật độ phù sa lơ lửng,…Trong đó vai trò chính là dòng chảy sát bờ và sóng.
Hình 8 cho thấy trong vùng tính, những khu vực bò xói thì xói sâu nhưng mà diện tích
bò xói lại hẹp, còn những khu vực được bồi thì bồi ít nhưng mà diện tích bồi rộng hơn.
Điều này chứng tỏ được sự cân bằng vật chất tương đối trong vùng tính. Nhưng nhìn một
cách tổng thể thì vùng biển Cần Giờ vẫn được xem là vùng bồi tụ trầm tích.
VI. KẾT LUẬN
Từ việc so sánh về kích thước bãi triều, phù sa trong vùng tính toán, số liệu thực đo
hàm lượng phù sa lơ lửng với kết quả tính được từ mô hình cho thấy rằng mô hình tính
toán có độ tin cậy và nếu bổ sung tính toán thêm ảnh hưởng của sóng thì ta có khả năng
áp dụng trong thực tế để tính toán. Đây là cách làm ít tốn kém và rất tiện lợi.
Lời cám ơn: Các tác giả xin chân thành cám ơn sự hỗ trợ về mặt tài chính của chương
trình nghiên cứu cơ bản nhà nước 2006.
7
C A ÀN G I Ơ Ø
V U ÕN G T A ØU
Hình 4 : Kết quả tính toán mực nước ven biển Cần Giờ lúc triều lên
C A ÀN G I Ơ Ø
V U ÕN G T A ØU
Hình 5 : Kết quả tính toán mực nước ven biển Cần Giờ lúc triều xuống
1
0
2
0
3
0
4

0
5
0
0
0 . 0 2 5
0 . 0 5
0 . 0 7 5
0 . 0 7 5
0 . 1
0 . 1 2 5
0 . 1 5
0 . 1 7 5
0 . 2
0 . 2 2 5
0 . 2 5
k g / m 3
Hình 6 : Sự tập trung phù sa ở vùng ven biển Cần Giờ (kết quả tính toán lúc triều lên)
8
0
0 . 0 2 5
0 . 0 5
0 . 0 7 5
0 . 0 7 5
0 . 1
0 . 1 2 5
0 . 1 5
0 . 1 7 5
0 . 2
0 . 2 2 5
0 . 2 5

k g / m 3
Hình 7 : Sự tập trung phù sa ở vùng ven biển Cần Giờ (kết quả tính toán lúc triều xuống)
m
- 0 . 0 0 5
- 0 . 0 0 1
0
0 . 0 0 1
0 . 0 0 5
0 . 0 1
0 . 0 1 5
0 . 0 2
0 . 0 2 5
0 . 0 3
0 . 0 3 5
C A ÀN G I Ơ Ø
C A ÀN G I Ơ Ø
Hình 8 : Kết quả tính toán sự thay đổi đáy biển Cần Giờ
TÀI LIỆU THAM KHẢO :
[1] Đặng Cơng Minh, Nguyễn Hữu Nhân; “Thủy triều biển Đơng”, Chương trình NC cấp
nhà nước KT.03, đề tài KT.03.03, (1993).
[2] Leo C. van Rijn, “Handbook- Sediment Transport by Currents and waves”.Delft
Hydaulic June 1989.
[3] Leo C. van Rijn, Principles of Sediment Transport in rivers, estuaries and coastal seas”.
Delft Hydaulic June 1993.
[4] Lê Trình “Nghiên cứu khả năng tác động của cơng nghiệp hố, đơ thị hố , giao thơng
thuỷ đến mơi trường Cần Giờ”. Đề tà i cấp sở KHCN & MT tp HCM(2001)
[5] Lê Quang Toại, Nguyễn Thế Hào, (1995), Tính tốn phân bố phù sa lơ lửng vùng cửa
sơng Sồi Rạp, Báo cáo đề tài khoa học cấp Tổng cục, Tổng cục khí tượng
[6]. Phan văn Hoặc, Nguyễn Hữu Nhân. “Nghiên cứu xâm nhập mặn trên sơng Đồng Nai
phục vụ nhà máy nước 100.000m

3
/ ngày”. Tổng cục Khí Tượng Thủy Văn, phân viện Khí
Tượng Thủy Văn tại TPHCM, 1993.
[7] Phan Văn Hoặc (2004), Báo cáo đề tài: “Nghiên cứu tương tác động lực học biển – sơng
ven biển Cần Giờ phục vụ xây dựng cơ sở hạ tầng cho du lịch TP. Hồ Chí Minh”, Sở Khoa
học và Cơng nghệ thành phố Hồ Chí Minh, Tp Hồ Chí Minh.
[8] Trung tâm Khí tượng thủy văn phía Nam (2000), “Vai trò của thuỷ triều trong vấn đề
ngập lụt tại tp HCM”, tp HCM 12/2000.
[9]. Ủy ban nhân dân tp HCM, Liên hiệp các hội khoa học và kỹ thuật. “Những luận cứ
khoa học làm cơ sở cho việc cải tạo bãi biển Cần Giờ thành khu bãi tắm, nghỉ ngơi, giải
trí”, tp. HCM, 2001.
9
Hình P16-Bờ biển và bãi triều Cần Giờ (Ảnh vệ tinh 12 / 7 /
1996)
[10]. U.Lumborg,A.Windelin. “Hydrography and cohesive sediment modelling:application
to the Romo Dyb tidal area”. J. of Marine systems, 38(2003) p.287-303.
[11] Bui Trong Vinh, Ichiro DEGUCHI (2004), The potential application of remote sensing
& gis and numerical models to investigate coastal process in Can Gio region (Saigon river
mouth – South Vietnam), Faculty of Geology & Petroleum, HoChiMinh City University of
Technology, Ho Chi Minh City 2004.
[12] Ủy ban nhân dân Tp.HCM, Liên hiệp các hội khoa học và kỹ thuật. “Những luận cứ
khoa học làm cơ sở cho việc cải tạo bãi biển Cần Giờ thành khu bãi tắm, nghỉ ngơi, giải
trí”, Tp.HCM, 2001.
10