Nghiên cứu cơ bản về thủy âm và mô phỏng đánh giá các tham số chính

TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT Tập 8, Số 2, 2018 59–75 
59 
NGHIÊN CỨU CƠ BẢN VỀ THỦY ÂM 
VÀ MÔ PHỎNG ĐÁNH GIÁ CÁC THAM SỐ CHÍNH 
Phan Thanh Minha*, Nguyễn Ngọc Bìnha 
aKhoa Nghiệp vụ Viễn Thông, Trường Đại học Thông tin Liên lạc, Khánh Hòa, Việt Nam 
*Tác giả liên hệ: Email: phanthanhminh@tcu.edu.vn 
Lịch sử bài báo 
Nhận ngày 22 tháng 01 năm 2018 
Chỉnh sửa ngày 13 tháng 05 năm 2018 | Chấp nhận đăng ngày 14 tháng 05 năm 2018 
Tóm tắt 
Thông tin dưới nước cùng với các ứng dụng của nó là một lĩnh vực nghiên cứu đã và đang 
được phát triển nhanh chóng, mở rộng trong nhiều lĩnh vực như: Điều khiển xa trong ngành 
công nghiệp khai thác dầu mỏ ở ngoài biển; Tính toán ô nhiễm môi trường trong các hệ 
thống thuộc về môi trường; Truyền tiếng nói giữa các người nhái; Vẽ đáy đại dương để tìm 
ra các nguồn tài nguyên mới; Thông tin liên lạc giữa các thiết bị ngầm; Có hai cách thiết 
lập việc trao đổi thông tin giữa các thiết bị dưới nước: Cách thứ nhất là kết nối bằng cáp 
giữa máy phát và máy thu, cách này bảo đảm chất lượng tín hiệu tốt và giảm thiểu những 
tác động không mong muốn của môi trường. Tuy nhiên, chi phí cho việc triển khai bảo đảm 
liên lạc cao, công tác bảo quản và bảo dưỡng khó khăn, đặc biệt nếu việc trao đổi thông tin 
diễn ra ở độ sâu lớn, trong điều kiện cơ động thì đảm bảo thông tin theo kiểu này rất phức 
tạp. Cách thứ hai là thiết lập thông tin giữa các thiết bị bằng cách sử dụng nước như một 
môi trường truyền dẫn tín hiệu và kênh thông tin vô tuyến dưới nước như vậy được gọi là 
kênh thủy âm. Bài báo nghiên cứu những vấn đề cơ bản để thực hiện truyền tin dưới nước 
bằng thủy âm, trong đó trình bày cơ sở chung để thực hiện thông tin dưới nước, mô phỏng 
hệ thống thông tin thủy âm sử dụng kỹ thuật điều chế QPSK (Quadrature Phase Shift 
Keying), kết luận đánh giá và đề xuất hướng nghiên cứu tiếp theo. 
Từ khóa: QPSK; Thủy âm; Truyền dẫn. 
Mã số định danh bài báo:  
Loại bài báo: Bài báo nghiên cứu gốc có bình duyệt 
Bản quyền © 2018 (Các) Tác giả. 
Cấp phép: Bài báo này được cấp phép theo CC BY-NC-ND 4.0 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
60 
RESEARCH ON UNDERWATER COMMUNICATION AND 
SIMULATED ASSESSMENT OF MAIN PARAMETERS 
Phan Thanh Minha*, Nguyen Ngoc Binha 
aThe Faculty of Professional Telecommunications, Telecommunications University, Khanhhoa, Vietnam 
*Corresponding author: Email: phanthanhminh@tcu.edu.vn 
Article history 
Received: January 22nd, 2018 
Received in revised form: May 13th, 2018 | Accepted: May 14th, 2018 
Abstract 
Underwater communication and its applications is a research field that has been 
developing rapidly, extending in many fields such as remote control in the offshore oil 
industry, calculating environmental pollution, transmitting voice among frogmen, drawing 
the ocean floor to find new resources, communication between underground devices, etc. 
There are two ways of establishing information exchange between underwater devices. The 
simplest and most effective way is using a cable connection between the transmitter and the 
receiver, which ensures high quality of signal and minimizes unwanted effects of the 
environment. However, it has disadvantages such as high cost of deploying communication, 
difficult maintenance, and especially, if information exchanges take place at great depths 
and in mobile cases it would be very complex to ensure such an exchange. The second way 
is establishing information between devices using water as a signal transmission medium, 
which is called underwater communication channel. This article focuses on the basics of 
underwater communications including the general basis for conducting underwater 
communications, simulation of an underwater information system using QPSK modulation 
technique, and conclusions and recommendations for further researches. 
Keywords: QPSK modulation technique; Signal transmission; Underwater acoustic. 
Article identifier:  
Article type: (peer-reviewed) Full-length research article 
Copyright © 2018 The author(s). 
Licensing: This article is licensed under a CC BY-NC-ND 4.0 
Phan Thanh Minh và Nguyễn Ngọc Bình 
61 
1. CỞ SỞ CHUNG ĐỂ THỰC HIỆN THÔNG TIN DƯỚI NƯỚC 
1.1. Tại sao lại dùng âm thanh để truyền tin dưới nước? 
Bảng thống kê các tham số trong Bảng 1 cho thấy rằng sóng âm là tín hiệu thuận 
lợi nhất để truyền tin trong môi trường nước (suy hao truyền dẫn nhỏ, cự li liên lạc xa). 
Đại dương là một môi trường vô cùng phức tạp, tất cả đặc trưng riêng biệt của môi 
trường đại dương là tính tự nhiên không đồng nhất của nó. Có hai loại không đồng nhất: 
(i) Thông thường; và (ii) Ngẫu nhiên. Cùng với ảnh hưởng mạnh mẽ của trường sóng 
âm trong đại dương, sự thay đổi thông thường của tốc độ sóng âm so với các tham số 
của đại dương (nhiệt độ, độ sâu, độ mặn) đã dẫn tới hình thành “kênh sóng âm dưới 
nước” và điều đó là một hệ quả để truyền lan âm sóng đi xa. Sự không đồng nhất ngẫu 
nhiên của môi trường đại dương làm gia tăng sự tán xạ của sóng âm, do đó tạo ra những 
thay đổi thất thường trong trường sóng âm. 
Bảng 1. So sánh các loại tín hiệu sóng âm, điện từ, quang trong môi trường nước 
Tham số Acoustic (Sóng âm) Electromagnetic (Sóng điện từ) Optical (Sóng quang) 
Nominal (m/s) ~ 1,500  30,000,000  30,000,000 
Power Loss >0.1 dB/m/Hz ~28 dB/1Km/100MHz α Turbidity 
Bandwidth kHz MHz (10-150) MHz 
Frequency band kHz MHz THz 
Antenna size 0.1m 0.5m 0.1m 
Effective range Km 10m 10-100m 
Nguồn: Tác giả tổng hợp từ Internet. 
1.2. Ảnh hưởng của đại dương đến truyền sóng âm dưới nước 
1.2.1. Tốc độ âm thanh trong đại dương 
Sự thay đổi của tốc độ âm thanh trong đại dương c là tương đối nhỏ. Như là một 
qui luật, c nằm giữa 1450 và 1540m/s. Nhưng ngay cả các thay đổi nhỏ của c cũng ảnh 
hưởng đáng kể tới sự truyền lan của âm thanh trong đại dương. Một số thí nghiệm và đo 
lường hiện nay chỉ ra rằng tốc độ âm thanh tăng một cách phức tạp khi nhiệt độ, áp suất 
thủy tĩnh (hoặc độ sâu) và lượng muối hòa tan trong nước tăng. Tốc độ âm thanh c được 
đưa ra bởi Medwin và Clay (1998): 
c= 1449.2+4.6T-0.055T2+(1.34-0.01T)(S-35)+0.016z (1) 
Ở đây nhiệt độ T được tính bằng [0C], S độ mặn  000 , độ sâu z tính bằng [m], 
tốc độ âm thanh tính bằng [m/s]. 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
62 
1.2.2. Sự phụ thuộc của tốc độ âm thanh vào nhiệt độ, độ mặn, độ sâu 
Dạng đặc trưng của tốc độ âm thanh là một hàm của T, S và z là: 
( ) ( ( ), ( ), )c z c T z S z z= (2) 
Công thức (2) biểu thị mối quan hệ giữa T, S, z và c là quan trọng nhất đối với sự 
truyền lan của âm thanh trong đại dương. Đặc trưng của tốc độ âm thanh c theo độ sâu z 
là khác nhau đối với các vùng biển khác nhau, thay đổi theo thời gian (theo mùa). Tại 
các độ sâu dưới 1km, sự thay đổi của T và S thường nhỏ và tốc độ âm thanh tăng phần 
lớn do sự tăng của áp suất thủy tĩnh. Kết quả là tốc độ âm thanh tăng tuyến tính với độ 
sâu. 
Hình 1 chỉ ra nét đặc trưng thông thường của nhiệt độ tại bề mặt của biển cao 
hơn tại đáy biển. Ở đây có thể thấy, nhiệt độ giảm so với độ sâu đến độ sâu z=300m và 
sau đó đạt giá trị không đổi. Điều này phù hợp với đặc trưng mùa hè của một vùng biển 
thông thường. Tốc độ âm thanh thay đổi cùng với sự thay đổi của nhiệt độ, độ mặn và 
độ sâu. Sự ảnh hưởng của nhiệt độ và áp suất đối với tốc độ âm thanh được chỉ ra trong 
Hình 2. Ở đây nó được phân thành ba vùng: Vùng thứ nhất nhiệt độ là nhân tố chi phối 
đến tốc độ âm thanh; Vùng thứ hai là vùng chuyển đổi có độ sâu từ 200-400m, cả nhiệt 
độ và độ sâu đều chi phối đến tốc độ âm thanh; và Vùng thứ ba trên 400m, tốc độ âm 
thanh hoàn toàn phụ thuộc vào độ sâu. 
Hình 1. Nhiệt độ và độ sâu 
Phan Thanh Minh và Nguyễn Ngọc Bình 
63 
Hình 2. Tốc độ âm thanh và độ sâu 
Sự phụ thuộc của c vào độ mặn S được chỉ ra trong Hình 3. Ở đây, với sự gia 
tăng của S, tốc độ âm thanh c cũng tăng và giữ nguyên dạng đặc trưng của nó. 
Hình 3. Tốc độ âm thanh và độ mặn 
1.2.3. Suy hao truyền dẫn của âm thanh trong môi trường nước 
• Tán xạ hình cầu: Trong một môi trường đồng nhất và kéo dài vô tận, công 
suất phát ra bởi một nguồn được phát xạ theo tất cả các hướng trên bề mặt 
của một hình cầu. Điều này được gọi là phân bố hình cầu. 
2
0
0
( ) 10log 20logsphere
sphere
I r
G r
I r
= = 
 (3) 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
64 
Trong đó: r0 là khoảng cách tham khảo (=1m); I0 là cường độ âm thanh của 
nguồn tại khoảng cách r0; và I là cường độ âm thanh của nguồn tại khoảng cách r. 
• Tán xạ hình trụ: Tán xạ hình trụ tồn tại khi hai môi trường bị giới hạn bởi 
hai mặt phẳng phản xạ. Khoảng cách giữa các mặt phẳng được giả định là h 
> 10λ. Trong đó, λ biểu thị bước sóng của sóng âm. 
0
0
( ) 10log 10logcylinder
cylinder
I r
G r
I r
= = 
 (4) 
Năng lượng âm thanh của quá trình truyền sóng âm trong đại dương một phần bị 
hấp thụ do năng lượng bị chuyển thành nhiệt, một phần bị suy giảm do sóng âm bị tán 
xạ bởi sự không đồng nhất của môi trường. Nhưng trong các thí nghiệm ở đại dương 
thực, không thể phân biệt được giữa ảnh hưởng của sự hấp thụ và sự tán xạ. Cả hai hiện 
tượng trên góp phần làm suy giảm âm thanh trong nước biển. 
Trên cơ sở những thí nghiệm rộng rãi và thực nghiệm, công thức thực nghiệm 
sau đây đối với hệ số suy giảm trong nước biển được đưa ra. Công thức này thỏa mãn 
khi tần số nằm trong miền (a) như Hình 4. 
 
2 2
2 2
0.11 44
dB/km
1 4100
f f
f f
 = +
+ +
 (5) 
Trong đó : f là tần số [kHz]. 
Theo Brekhovskikh và Lysanov (2003) thì công thức này thỏa mãn khi tần số 
nằm trong miền (b) như Hình 4. 
( ) 
2 2
3 4
2
A
8.686*10 1 6.54*10 dB/kmT
T T
S f f Bf
P
f f
 −
= + − 
 (6) 
Trong đó:
6 62.34*10 , 3.38*10A B− −= = ; S là độ mặn [ppt]; P là áp suất thủy 
tĩnh [kg/cm2]; f là tần số [kHz]; và  6 1520/( 273)21.9*10 TTf kHz
− += là tần số nghỉ với T là 
nhiệt độ tính bằng [0C], với dải nhiệt độ từ 00 đến 300C, fT thay đổi xấp xỉ từ 59 đến 210 
kHz. Công thức này thỏa mãn khi tần số nằm trong miền (c) như Hình 4. 
( )
 
2
43
2 2
21 1 1 2 2 2
3 32 2 2 2
1 2 , Pure water
, Magnesium sulphate, Boric acid
dB/km
H O
MgSOB OH
A P f f A P f f
A P f
f f f f
 = + +
+ +
 (7) 
Trong đó: Boric axít B(OH)3 
Phan Thanh Minh và Nguyễn Ngọc Bình 
65 
1245
4
0.78 5 273
1 1 1
8.686
*10 ; 1; 2.8 *10
35
ph T
S
A P f
c
−
− += = = (8) 
Hình 4. Biểu đồ chỉ thị công thức thực nghiệm với các miền tần số khác nhau 
Magiê sunfat MgSO4 
( )
( )
( )
1990
8
273
4 9 2
2 2 ax ax 2
8.17*10
21.44* 1 0.025 ; 1 1.37*10 6.2*10 ;
1 0.0018 35
T
m m
S
A T P z z f
c S
−
+
− −= + = − + =
+ −
 (9) 
Tốc độ âm thanh xấp xỉ được cho bởi: 
ax1412 3.21 1.19 0.0167 mc T S z= + + + (10) 
Trong các miền chuyển tiếp nó tương ứng với f. Sự suy giảm tăng tương ứng với 
sự tăng của độ mặn và nhiệt độ, Hình 5 thể hiện sự suy giảm tăng tương ứng với tần số 
tăng. Khi nhiệt độ tăng thì độ suy giảm tăng (Hình 6). 
Hình 5. Hệ số suy giảm thay đổi theo độ mặn và nhiệt độ 200C 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
66 
Hình 6. Hệ số suy giảm thay đổi theo độ mặn và nhiệt độ 300C 
• Sự suy giảm của âm thanh trong lớp trầm tích: Sự suy giảm của âm thanh 
trong lớp trầm tích chủ yếu phụ thuộc vào sự thay đổi dạng của đáy. Cấu 
trúc của lớp đáy được tượng trưng bởi bt, biểu thị chất liệu lớp trầm tích của 
đại dương. Bảng 2 đưa ra các giá trị của bt cho mỗi dạng của lớp trầm tích. 
Công thức thực nghiệm sau được đưa ra để tìm sự suy giảm của sóng âm 
trong lớp trầm tích phụ thuộc vào bt. 
1 1
8.686 1kHz m
n
s
f
K 
= 
 (11) 
Trong đó: αs là sự suy giảm của sóng âm trong lớp trầm tích. 
Bảng 2. Các giá trị của bt cho mỗi dạng của lớp trầm tích 
Dạng của lớp trầm tích Giá trị bt 
Cát rất thô 0 
Cát thô 1 
Cát cỡ trung bình 2 
Cát mịn 3 
Cát rất mịn 4 
Bùn rất thô 5 
Bùn thô 6 
Bùn cỡ trung bình 7 
Bùn mịn 8 
Bùn rất mịn 9 
Đất sét 10 
Phan Thanh Minh và Nguyễn Ngọc Bình 
67 
Bảng 3 đưa ra các giá trị của K và n đối với bốn dạng của lớp trầm tích. 
Bảng 3. Các giá trị của K và n đối với bốn dạng của lớp trầm tích. 
Dạng của lớp trầm tích K n 
Bùn rất mịn 0.17 0.96 
Cát mịn 0.45 1.02 
Cát cỡ trung bình 0.48 0.98 
Cát thô 0.53 0.96 
1.2.4. Hiệu ứng Doppler 
Chuyển động tương đối giữa máy thu và máy phát gây nên thay đổi trong đáp 
ứng kênh truyền do hiệu ứng Doppler. Biên độ của hiệu ứng Doppler tỉ lệ với tỉ số a = 
v/c, trong đó v là chuyển động tương đối giữa máy thu và máy phát, c là vận tốc âm 
thanh trong nước. Do vận tốc của âm thanh trong nước là khá nhỏ khi so sánh với vận 
tốc của sóng điện từ trên không trung nên ảnh hưởng của hiệu ứng Doppler là rất lớn. 
Các thiết bị tự động dưới nước di chuyển với tốc độ khoảng vài m/s, tuy nhiên kể cả khi 
không có những chuyển động có mục đích thì những ảnh hưởng như sự trôi dạt gây nên 
bởi sóng, thủy triều cũng luôn tồn tại. Nói cách khác luôn có chuyển động tương đối 
giữa máy thu và phát, hệ thống thông tin thủy âm cần được thiết kế nhằm giải quyết vấn 
đề này. Vấn đề này có nhiều điểm tương đồng với kênh vô tuyến truyền dẫn qua vệ tinh. 
Méo tín hiệu gây ra do chuyển động của thiết bị tác động đến việc thiết kế thuật toán 
đồng bộ và ước lượng kênh truyền. 
Mức độ ảnh hưởng của méo lên tín hiệu phụ thuộc trực tiếp và giá trị của a. Làm 
phép so sánh như sau: Với hệ thống thông tin vô tuyến, với độ dịch chuyển tương đối 
giữa nguồn thu và phát là 160km/h thì 
71.5 10a −= , giá trị này đủ nhỏ để ảnh hưởng 
của hiệu ứng Doppler có thể được bỏ qua. Nói cách khác, việc xem xét ảnh hưởng độ 
méo của tín hiệu trong quá trình đồng bộ là không cần thiết. Xác xuất lỗi bit là rất nhỏ. 
Ngược lại, với kênh thông tin thủy âm, giả sử nguồn phát và thu chịu sự dịch chuyển 0.5 
m/s (1 knot), thì giá trị của a vào khoảng 
43 10a −= . Nếu như dịch chuyển tương đối 
lên tới vài m/s (tàu ngầm có thể dịch chuyển với vận tốc lớn hơn), giá trị của a vào 
khoảng 
310a −= , và đây là giá trị không thể bỏ qua được. 
Dịch chuyển Doppler và trải phổ Doppler sinh ra do chuyển động tương đối là 
một trong những yếu tố khác biệt giữa kênh thông tin thủy âm với kênh thông tin vô 
tuyến. Trải phổ Doppler gây lệch pha, trễ đồng bộ. Trong hệ thống thông tin thủy âm 
dùng đa sóng mang, hiệu ứng Doppler còn gây ra méo đặc biệt nghiêm trọng. Trong 
kênh thông tin vô tuyến trên không trung, sự nén, giãn theo thời gian là có thể bỏ qua, 
hiệu ứng Doppler là như nhau với tất cả các sóng mang con. Với kênh thông tin thủy 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
68 
âm, ảnh hưởng của dịch chuyển Doppler lên mỗi sóng mang con là khác nhau đáng kể, 
điều này gây nên méo Doppler không đồng bộ trên toàn bộ băng thông tín hiệu. 
2. MÔ PHỎNG HỆ THỐNG THÔNG TIN THỦY ÂM SỬ DỤNG KỸ 
THUẬT ĐIỀU CHẾ QPSK 
2.1. Sơ đồ khối hệ thống thông tin thủy âm 
Để tiến tới một hệ thống rời rạc theo thời gian mong muốn, xét hệ thống thông 
tin như Hình 7. Đây là hệ thống thông thường được thiết lập mà có thể đặc trưng cho 
bất kỳ hệ thống nào sử dụng kiểu điều chế biên độ cầu phương QAM (Quadrature 
Amplitude Modulation). Tóm tắt về hệ thống được mô tả như sau: Tại phía phát, các ký 
hiệu tuần tự d(n) được chuyển đổi sang tín hiệu băng gốc liên tục theo thời gian sbb(t) 
bởi một bộ điều chế biên độ xung (Pulse Amplitude Modulation-PAM). Trong đó các giá 
trị d(n) là các số phức rời rạc. Chuyển đổi tần số đường lên được thực hiện bằng cách 
nhân tín hiệu với
2 cj f te
, kết quả là tín hiệu băng thông s(t) được truyền qua kênh 
truyền như đã phân tích. 
Hình 7. Hệ thống mô phỏng kênh thủy âm 
Tại đầu thu, để tách sóng mang, tín hiệu thu r(t) được xử lý bằng một bộ đổi tần 
xuống mà có tín hiệu đầu ra phù hợp với tín hiệu tương đương băng gốc rbb(t). Sau bộ 
chuyển đổi đường xuống là bộ lọc thông thấp và bộ lọc phối hợp. Bộ tách sóng thực hiện 
ước lượng các ký hiệu đã được phát. Như đề cập ở trên, các đặc trưng của tín hiệu băng 
gốc sẽ có nhiều thuận lợi cho việc mô phỏng hệ thống. Hình 8 tượng trưng cho hệ thống 
tương đương băng gốc. 
Hình 8. Hệ thống tương đương băng gốc 
Phan Thanh Minh và Nguyễn Ngọc Bình 
69 
Mô phỏng thường được dựa trên hệ thống oversampled, có nghĩa là tốc độ lấy 
mẫu cao hơn tốc độ ký hiệu. Nhìn chung tốc độ lấy mẫu cao hơn sẽ phản ảnh chính xác 
hơn hệ thống nguyên bản. Tuy nhiên điều này dẫn tới thời gian mô phỏng dài hơn do 
nhiều mẫu hơn cần phải xử lý. Nhìn chung để sử dụng tốc độ oversample, tốc độ này là 
bội số của tốc độ ký hiệu. Trong đó Q là số bit trên một ký hiệu, Q là một số tự nhiên. 
2.2. Xây dựng mô hình mô phỏng 
Hệ thống mô phỏng được minh họa trong Hình 9. Trong đó bao gồm một nguồn 
bit, máy phát, kênh truyền, máy thu và một bộ tính tỉ số lỗi bit. Nguồn bit tạo ra các bit 
nhị phân tuần tự ngẫu nhiên và nó được phát đi bởi máy phát. Thông thường một nguồn 
bit ngẫu nhiên được dùng trong các mô phỏng và được sử dụng trong mô phỏng của bài 
báo. Máy phát chuyển các bit sang ký hiệu QPSK, đưa qua bộ tạo dạng xung và chuyển 
đổi đường lên được thực hiện để đưa đến tần số sóng mang yêu cầu. 
Tín hiệu đầu ra của máy phát được phát đi qua kênh thủy âm. Máy thu nhận tín 
hiệu này từ kênh truyền, ước lượng pha và định thời, giải điều chế các tín hiệu QPSK 
nhận được thành các bit tin và cung cấp cho khối tính toán tỉ số lỗi bit. Ở đây việc tính 
số lượng các bit lỗi xuất hiện để khảo sát hiệu năng của hệ thống. 
Hình 9. Hệ thống mô phỏng 
2.3. Kết quả mô phỏng 
Sử dụng lập trình Matlap, tiến hành đánh giá một số tham số như sau: 
2.3.1. Kết quả mô phỏng đối với các loại đáy 
Tiến hành đánh giá cho các loại đáy mà giá trị bt được trình bày trong Bảng 2. 
Các tham số khác được mô tả như sau: Vị trí nguồn (rs,zs)=(0,25) m; Vị trí máy thu (
1R
r ,
1R
z ) = (500,25)m; Tốc độ âm thanh c=1500m/s; Độ sâu của nước D=40m; Độ mặn 
S=35 ppt; Nhiệt độ của nước T=140C; và Độ pH là pH=8. 
Nhận thấy rằng tín hiệu phản xạ từ đáy tới bên thu có biên độ rất khác nhau cho 
những môi trường khác nhau. Với môi trường mà đáy có tính chất phản xạ tốt, tín hiệu 
thu do sự phản xạ từ đáy đóng vai trò đáng kể. Như vậy, khi tiến hành triển khai hệ 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
70 
thống trên thực tế, việc nghiên cứu tính chất phản xạ của mặt đáy là vô cùng quan trọng. 
Sự ảnh hưởng của tính chất đa đường đến tín hiệu trong môi trường đáy bùn và đáy đất 
sét hầu như không có sự khác biệt nào đáng kể. 
(a) 
(b) 
 Hình 10. Tín hiệu đa đường trong môi trường đáy cát và đáy bùn 
Ghi chú: a) Môi trường đáy là cát; b) Môi trường đáy là bùn. 
(a) 
(b) 
Hình 11. Tín hiệu đa đường trong môi trường đáy bùn và đáy là đất sét 
Ghi chú: a) Môi trường đáy là bùn; b) Môi trường đáy là đất sét. 
Phan Thanh Minh và Nguyễn Ngọc Bình 
71 
2.3.2. Kết quả mô phỏng đối với trường hợp khoảng cách giữa máy phát và máy 
thu thay đổi 
Tham số mô phỏng (Hình 12) bao gồm: Tần số sóng mang là 30 kHz; Dạng đáy 
bt= 1; Vị trí nguồn (rs,
1s
z )=(0, 20) m; Vị trí máy thu (
1R
r ,
1R
z )= (x,20) m, với x= 10, 
500, 1000, 1500. 
(a) 
(b) 
(c) 
(d) 
Hình 12. Tính chất đa đường tại các khoảng cách khác nhau với đáy cát 
Ghi chú: a) 10m; b) 500m; c) 1000m; và d) 1500m. 
2.3.3. Đặc tuyến BER 
Trong thông tin, mục tiêu là phải thu được tỉ số tín hiệu trên tạp âm cao nhất. 
Đối với thông tin dưới nước thì nhiễu có hai dạng, một là nhiễu do môi trường xung 
quanh tạo nên và thứ hai là hiện tượng truyền dẫn đa đường của chính tín hiệu đó. Do 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
72 
vậy, chúng tôi xem xét ảnh hưởng của hiện tượng truyền dẫn đa đường của kênh thủy 
âm và tất cả các ảnh hưởng đối với kênh truyền này. Sau đây là một số kết quả mô 
phỏng chỉ ra tỉ số lỗi bit (BER) đối với trường hợp chỉ có đường trực tiếp và đa đường 
cho hai trường hợp khi cấu trúc của dạng đáy khác nhau và khoảng cách giữa máy phát 
và máy thu khác nhau. 
• Đặc tuyến BER khi chưa xét đến ảnh hưởng của đa đường 
Tham số mô phỏng: Tốc độ âm thanh c=1500m/s; Độ sâu của nước D=40m; 
Độ mặn S=35 ppt; Nhiệt độ của nước T=140C; Độ pH =8; Tốc độ gió vw=6knots; Vị trí 
nguồn (rs,zs)=(0,x)m, x=10.35; Vị trí máy thu (
1R
r ,
1R
z )= (1000,x)m; và Dạng đáy bt= 
bùn thô. Hình 13 biểu thị giá trị BER của đường trực tiếp. Trong trường hợp này, chúng 
tôi chỉ xét đến đường trực tiếp do đó không có bất cứ một nhiễu nào, vì vậy chỉ có sự 
suy giảm về độ lớn của tín hiệu nên ta có BER của đường trực tiếp bằng không. 
Hình 13. BER của tia trực tiếp 
• Đặc tuyến BER khi xét đến ảnh hưởng của truyền dẫn đa đường 
Xét hai trường hợp, trường hợp một dạng đáy là bùn và trường hợp thứ hai dạng 
đáy là cát khi có tác động của hiện tượng truyền dẫn đa đường. Tham số mô phỏng: Vị 
trí nguồn (rs,zs)=(0,10) m; Vị trí máy thu (rR,zR)=(500,35)m; Tốc độ âm thanh 
c=1500m/s; Độ sâu của nước D=40m; Độ mặn S=35ppt; Nhiệt độ của nước T=140C; Độ 
pH =8; Tốc độ gió w 6 knots = ; và Dạng đáy: bt= 0,7. 
Hình 14 và Hình 15 cho thấy: Với môi trường có đáy phản xạ bt=7, thành phần 
phản xạ sẽ có biên độ nhỏ hơn so với trường hợp đáy có hệ số phản xạ là bt=0, do đó 
Phan Thanh Minh và Nguyễn Ngọc Bình 
73 
ảnh hưởng của thành phần phản xạ sẽ nhỏ hơn. Vì vậy, tỉ lệ lỗi bit của trường hợp bt=7 
sẽ nhỏ hơn tỉ lệ lỗi bit của bt=0. 
(a) 
(b) 
Hình 14. Biểu đồ BER đối với trường hợp bt=0 
Ghi chú: a) Thang tuyến tính; b) Thang log. 
(a) 
(b) 
Hình 15. Biểu đồ BER đối với trường hợp bt=7 
Ghi chú: a) Thang tuyến tính; b) Thang log. 
TẠP CHÍ KHOA HỌC ĐẠI HỌC ĐÀ LẠT [ĐẶC SAN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG] 
74 
3. KẾT LUẬN 
Thủy âm ngày càng được ứng dụng mạnh mẽ và đạt hiệu quả cao trong phát 
triển kinh tế cũng như bảo vệ chủ quyền lãnh thổ quốc gia. Những kết quả nghiên cứu 
trên của bài báo đã đưa ra kết quả về những tác động của các yếu tố đến quá trình truyền 
tin trong môi trường nước biển. Mô phỏng đưa ra kết quả cụ thể cho một mô hình hệ 
thống thông tin thủy âm sử dụng kỹ thuật điều chế QPSK. Kết quả có thể được ứng 
dụng trong việc chế tạo các thiết bị truyền tin dưới nước, trên cơ sở đó tổ chức xây dựng 
hệ thống thông tin trên biển ngày càng hoàn thiện phục vụ cho các mục đích của quốc 
gia. 
TÀI LIỆU THAM KHẢO 
Brekhovskikh, L. M., & Lysanov, Y. P. (2003). Fundamentals of ocean acoustics (3rd 
ed.). Berlin, Germany: Springer Press. 
Chengsheng, P., Liangchen, J., Ruiyan, C., & Yuanming, D. (2012). Modeling and 
simulation of channel for underwater communication network. International 
Journal of Innovative Computing, 8(3B), 2149-2156. 
Jensen, F. B., Kuperman, W. A., Porter, M. B., & Schmidt, H. (2000). Computational 
ocean acoustics. Berlin, Germany: Springer Press. 
John, G. P. (2001). Digital communications (4th ed.). New York, USA: McGraw-Hill. 
Johnny, R. J. (1996). Introduction to digital signal processing. New Delhi, India: 
Prentice-Hall of India. 
Kaya, A., & Yauchi, S. (2002). An acoustic communication system for subsea robot. 
Paper presented at The IEEE OCEANS Conference, USA. 
Kilfoyle, D., Preisig, J., & Baggeroer, A. (2005). Spatial modulation experiments in the 
underwater acoustic channel. IEEE Journal of Oceanic Engineering, 30(2), 406-
415. 
Kumar, R., Thakur, N., & Thakur, V. (2013). An overview of sonar and acoustic 
underwater communication. International Journal of Advanced Research in 
Electrical, Electronics and Instrumentation Engineering, 2(5), 1997-2003. 
Medwin, H., & Clay, C. S. (1998). Fundamentals of acoustical oceanography. 
Massachusetts, USA: Academic Press. 
Phùng, B. T. (2012). Nghiên cứu kênh thông tin dưới nước sử dụng kỹ thuật điều chế 
pha QPSK. (Luận văn Thạc sĩ), Học viện Bưu chính Viễn thông, Việt Nam. 
Stojanovic, M. (1996). Recent advances in high rate underwater acoustic 
communications. IEEE Journal of Oceanic Engineering, 21(2), 125-136. 
Phan Thanh Minh và Nguyễn Ngọc Bình 
75 
Stojanovic, M., Catipovic, J. A., & Proakis, J. G. (1993). Adaptive multichannel 
combining and equalization for underwater acoustic communications. The 
Journal of the Acoustical Society of America, 94(3), 1621-1631. 
Stojanovic, M. (2008). OFDM for underwater acoustic communications: Adaptive 
synchronization and sparse channel estimation. Paper presented at The IEEE 
International Conference on Acoustics, Speech and Signal Processing, USA. 
Suzuki, M., Sasaki, T., & Tsuchiya, T. (1992). Digital acoustic image transmission 
system for deep-sea research submersible. Paper presented at The Mastering the 
Oceans Through Technology Conference, USA. 
Urban, H. G. (2006). Handbook of underwater acoustic engineering. The Journal of the 
Acoustical Society of America, 120(5), 2393-2410.