Xâm nhập mặn ở đồng bằng sông cửu long do biến đổi khí hậu

Tạimột số tỉnh ven biển ĐBSCL, nước biển xâm nhập sâu vào các sông rạch khiến các dòng sông bịnhiễm mặn sớm, gây ảnh hưởng nghiêm trọng đến đời sống người dân và hoạt động nông nghiệp.Nh

ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA MÔI TRƯỜNG & TÀI NGUYÊN

- -BÁO CÁO PHÂN TÍCH HỆ THỐNG MÔI TRƯỜNG

ĐỀ TÀI:

XÂM NHẬP MẶN Ở ĐỒNG BẰNG SÔNG CỬU LONG DO BIẾN ĐỔI KHÍ HẬU

LỚP L01- NHÓM 05-HK 211-KHÓA 2019

GVHD: Nguyễn Hoàng Anh

MỤC LỤC

MỤC LỤC 1

LỜI CẢM ƠN 3

VAI TRÒ VÀ ĐÓNG GÓP CỦA THÀNH VIÊN TRONG NHÓM 4

DANH MỤC CÁC TỪ VÀ CÁC KÝ HIỆU VIẾT TẮT 6

DANH MỤC CÁC BẢNG 7

DANH MỤC CÁC BIỂU ĐỒ 9

DANH MỤC CÁC HÌNH VẼ 11

MỞ ĐẦU 12

1 ĐỀ TÀI NGHIÊN CỨU: 12

2 ĐỐI TƯỢNG NGHIÊN CỨU 12

3 PHẠM VI NGHIÊN CỨU 13

5 PHƯƠNG PHÁP NGHIÊN CỨU 13

6 KẾT CẤU CỦA ĐỀ TÀI 13

7 Ý NGHĨA CỦA ĐỀ TÀI 14

CHƯƠNG 1: TÍNH CẤP THIẾT CỦA ĐỀ TÀI XÂM NHẬP MẶN TẠI ĐỒNG BẰNG SÔNG CỬU LONG 15

1.1 HIỆN TRẠNG VẤN ĐỀ XÂM NHẬP MẶN TẠI ĐBSCL: 15

1.1.1 Khái niệm: 15

1.1.2 Diễn biến xâm nhập mặn tại Việt Nam do BĐKH: 15

1.1.3 Các yếu tố ảnh hưởng đến xâm nhập mặn ở ĐBSCL: 17

1.2 TRÌNH BÀY VÀ MÔ TẢ DỮ LIỆU: 28

1.2.1 Cơ sở lý thuyết kiểm định thống kê: 28

1.2.3 Thống kê mô tả, hệ số tương quan và kiểm định của các biến: 34

1.2.4 Ma trận hệ số tương quan: 69

1.2.5 So sánh kết quả hệ số tương quan và kiểm định giả thiết trong đánh giá tương quan các yếu tố với Độ mặn trung bình S: 71

1.3 MỐI QUAN HỆ GIỮA CÁC YẾU TỐ DẪN ĐẾN XÂM NHẬP MẶN: 72

1.3.1 Cơ sở lý thuyết sơ đồ nguyên nhân-hậu quả (sơ đồ xương cá): 72

1.3.2 Xác định và trình bày dưới dạng biểu đồ mối quan hệ của các thành phần: 74

CHƯƠNG 2: PHÂN TÍCH QUYẾT ĐỊNH & ĐỀ XUẤT GIẢI PHÁP 75

2.1 QUY TRÌNH VÀ CƠ SỞ LÝ THUYẾT CỦA CÁC PHƯƠNG PHÁP ÁP DỤNG PHÂN TÍCH RA QUYẾT ĐỊNH: 75

2.1.1 Xây dựng quy trình: 75

2.1.2 Cơ sở lý thuyết của các phương pháp áp dụng: 75

2.2 PHÂN TÍCH RA QUYẾT ĐỊNH VÀ ĐỀ XUẤT GIẢI PHÁP: 95

2.2.1 Áp dụng phân tích AHP vào vấn đề xâm nhập mặn ở ĐBSCL: 95

2.2.2 Phân tích SWOT: 102

2.2.3 Phân tích What-if: 106

2.2.4 Áp dụng phương pháp Trade-off Analysis cho đề tài: 112

2.2.5 Phân tích NPV - Ma trận ra quyết định: 114

2.2.6 Phân tích dịch vụ hệ sinh thái: 119

CHƯƠNG 3: TÓM TẮT NỘI DUNG VÀ KẾT QUẢ PHÂN TÍCH 126

CHƯƠNG 4: TỔNG KẾT 129

TÀI LIỆU THAM KHẢO 130

LỜI CẢM ƠN

Đầu tiên, nhóm chúng em xin gửi lời cảm ơn chân thành đến cô Nguyễn Hoàng Anh đã hướng dẫn và giúp đỡ trong suốt quá trình thực hiện và hoàn thành tốt đề tài này Cảm ơn đến tất

cả các bạn chung lớp vì đã góp ý, phản biện khi nhóm thuyết trình

Trong suốt học kì vừa qua, tuy rằng chưa từng được gặp cô trực tiếp, nhưng chúng em đều nhận thấy được cô là một giáo viên đầy tâm huyết với nghề Phương pháp dạy học qua tiểu luận các tuần của cô đã cho chúng em cơ hội cọ xát thực tế, trải nghiệm như thế nào là được tham gia một dự án từ đầu đến cuối Thực sự, chúng em đã gặp phải nhiều khó ăn khi bắt đầu, nhưng nhờ

có sự chỉ dẫn tận tình của cô, mà chúng em dần hoàn thiện được những thiếu sót của tiểu luận Côcho chúng em kiến thức, chúng em phải tự mày mò thực tế, cũng chính nhờ vậy mà những bài học trên lớp không còn khô khan, mà đầy tính trực quan bằng tiểu luận, hoạt động nhóm, phản biện,

Mặc dù đã có nhiều cố gắng trong suốt quá trình thực hiện đề tài, song, có thể còn có những mặt hạn chế, thiếu sót Nhóm chúng em rất mong nhận được ý kiến đóng góp và sự chỉ dẫncủa cô

Kính chúc cô luôn dồi dào sức khỏe và đạt nhiều thành công hơn nữa trong sự nghiệp dạy học và nghiên cứu của mình

Chúng em xin chân thành cảm ơn!

VAI TRÒ VÀ ĐÓNG GÓP CỦA THÀNH VIÊN TRONG NHÓM

Bài số 1: Làm Power Point + Thuyết trìnhTham gia tìm số liệu cả 10 bài và đảm nhận chínhphần tính toán vẽ biểu đồ phân tích dữ liệu

Bài số 10: Làm Power Point + Thuyết trìnhTổng hợp nội dung báo cáo

Trương Hồng Hạnh 1913280

Tham gia tìm số liệu cả 10 bài và tính toán vẽ biểu

đồ phân tích dữ liệuBài số 6: Làm Power Point + Thuyết trìnhBài số 10: Làm Power Point + Thuyết trìnhTổng hợp nội dung báo cáo

Trà Huỳnh Uyển Nhi 1911790

Tham gia tìm số liệu cả 10 bài và tính toán phântích dữ liệu

Bài số 6: Thuyết trìnhBài số 8: Thuyết trìnhBài số 10: Làm Power PointTổng hợp nội dung báo cáo

tích dữ liệu

Bài số 6: Thuyết trìnhBài số 8: Làm Power PointBài số 10: Thuyết trìnhTổng hợp nội dung báo cáo

DANH MỤC CÁC TỪ VÀ CÁC KÝ HIỆU VIẾT TẮT

BĐKH: Biến đổi khí hậu

DANH MỤC CÁC BẢNG

Bảng 1: Bảng số liệu các biến liên quan đến xâm nhập mặn 32

Bảng 2: So sánh kết quả hệ số tương quan và kiểm định giá thiết trong 69

Bảng 3: Bản chỉ số ngẫu nhiên RI với nhân tố n 79

Bảng 4: Strengths – Thế mạnh trong SWOT 81

Bảng 5: Weakness – Điểm yếu trong SWOT 83

Bảng 6: Opportunities – Cơ hội trong SWOT 84

Bảng 7: Threats – Rủi ro trong SWOT 85

Bảng 8: Bảng phân tích chi phí lợi ích 90

Bảng 9: Bảng so sánh thứ hạng quan trọng qua phân tích hệ số tương quan và phân tích AHP 100

Bảng 10: Bảng yếu tố - giải pháp 102

Bảng 11: Bảng SWOT cho giái pháp 103

Bảng 12: Bảng SWOT cho giải pháp 104

Bảng 13: Bảng SWOT cho giải pháp 105

Bảng 14: Bảng SWOT cho giải pháp 106

Bảng 15: Bảng SWOT cho giải pháp 107

Bảng 16: Bảng SWOT cho giải pháp Trồng rừng ngập mặn 108

Bảng 17: Bảng SWOT cho giải pháp 108

Bảng 18: Bảng số liệu về các phương án phòng chống xâm nhập mặn 110

Bảng 19: Chi phí của phương án 1 111

Bảng 20: What-if cho phương án 1 111

Bảng 21: Chi phí của phương án 2 112

Bảng 22: What-if cho phương án 2 112

Bảng 23: Chi phí của phương án 3 113

Bảng 24: What-if của phương án 3 113

Bảng 25: Bảng phân tích What-if 113 Bảng 26: Ma trận Trade-off Analysis cho các kịch bản giảm thiểu xâm nhập mặn tại ĐBSCL .114

Bảng 28: Bảng chi phí – lợi nhuận cho các kịch bản 117

Bảng 29: Bảng Chi phí – Lợi nuận cho Năng lượng điện 117

Bảng 30: Bảng Chi phí – Lợi ích cho nhà máy khử mặn 118

Bảng 31: Bảng Chi phí – Lợi ích cho xây dựng thêm tuyến đê biển 118

Bảng 32: Bảng Chi phí – Lợi ích trồng rừng ngập mặn dọc biển Đông và biển Tây 118

Bảng 33: Bảng Net benefits của các giải pháp ứng với kịch bản 118

Bảng 34: Ma trận trọng số chuẩn hóa dựa trên % Maximum 119

Bảng 35: Ma trận trọng số chuẩn hóa dựa trên % Range 119

Bảng 36: Ma trận trọng số chuẩn hóa dựa trên % Total 119

Bảng 37: Ma trận trọng số chuẩn hóa dựa trên Unit Vector 119

Bảng 38: Ma trận quyết định của từng giải pháp ứng với kịch bản và tổng thể kịch bản 119

Bảng 39: Đánh giá hiệu quả tức thời và lâu dài của các đề xuất 121

Bảng 40: Bảng các chỉ tiêu dịch vụ hệ sinh thái 124

Bảng 41: Bảng phân tích giá trị dịch vụ hệ sinh thái 125

Bảng 42: Bảng đánh giá các khu vực bị ảnh hưởng nghiêm trọng tại ĐBSCL 126

Bảng 43: Phân tích giá trị tổn thất (đơn vị tiền tệ/năm) của việc thất thoát các hệ sinh thái 127

Bảng 44: Bảng Tổng hợp DVHST và giá trị Lợi ích - Chi phí 127

Bảng 45: Bảng so sánh giá trị NPV và B/C trước và sau khi phân tích giá trị DVHST 128

DANH MỤC CÁC BIỂU ĐỒ

Biểu đồ 1: Biểu đồ lượng mưa và nhiệt độ trung bình qua các năm 18

Biểu đồ 2: Biểu đồ lượng bốc hơi 4 năm tại 4 địa điểm chính vùng ĐBSCL 19

Biểu đồ 3: Biểu đồ độ mặn qua các năm tại ĐBSCL 19

Biểu đồ 4: Dòng chảy trên dòng chính sông Mekong (xét tại trạm Kratie) 22

Biểu đồ 5: Dòng chảy trên dòng chính sông Mekong (xét tại trạm PrekDam) 23

Biểu đồ 6: Biểu đồ độ mặn lớn nhất đến cuối thnasg 2/2016 24

Biểu đồ 7: Biểu đồ nhu cầu nước ngọt giữa các năm tại đồng bằng sông cửu long 25

Biểu đồ 8: Biểu đồ xây dựng đập thủy điện trên sông Mekong từ 1966 - 2020 26

Biểu đồ 9: Bảng thống kê mô tả và biểu đồ boxplot của Độ mặn trung bình S 34

Biểu đồ 10: Diễn tiến của độ mặn trung bình S qua các năm 35

Biểu đồ 11: Biểu đồ histogram của Độ mặn trung bình S 36

Biểu đồ 12: Bảng thống kê mô tả và biểu đồ boxplot của Độ cao mực nước biển 37

Biểu đồ 13: Diễn tiến Độ cao mực nước biển qua các năm 38

Biểu đồ 14: Biểu đồ histogram của Độ cao mực nước biển 39

Biểu đồ 15: Biểu đồ scatter và biểu đồ tương quan tuyến tính giữa 40

Biểu đồ 16: Bảng thống kê mô tả và biểu đồ boxplot của Khai thác nước ngầm 42

Biểu đồ 17: Diễn tiến Khai thác nước ngầm qua các năm 43

Biểu đồ 18: Biểu đồ histogram của Khai thác nước ngầm 43

Biểu đồ 19: Biểu đồ scatter và biểu đồ tương quan tuyến tính giữa 44

Biểu đồ 20: Bảng thống kê mô tả và biểu đồ boxplot của Nhiệt độ trung bình năm 46

Biểu đồ 21: Diễn tiến Nhiệt độ trung bình năm qua các năm 47

Biểu đồ 22: Biểu đồ histogram của Nhiệt độ trung bình năm 48

Biểu đồ 23: Biểu đồ scatter và biểu đồ tương quan tuyến tính giữa 48

Biểu đồ 24: Bảng thống kê mô tả và biểu đồ boxplot của Lượng mưa trung bình năm 50

Biểu đồ 25: Diễn tiến Lượng mưa trung bình năm qua các năm 51

Biểu đồ 26: Biểu đồ histogram của Lượng mưa trung bình năm 52

Biểu đồ 28: Bảng thống kê mô tả và biểu đồ boxplot của Lượng bốc hơi 55

Biểu đồ 29: Diễn tiến Lượng bốc hơi qua các năm 56

Biểu đồ 30: Biểu đồ histogram của Lượng bốc hơi 57

Biểu đồ 31: Biểu đồ scatter và biểu đồ tương quan tuyến tính giữa 57

Biểu đồ 32: Bảng thống kê mô tả và biểu đồ boxplot của Lưu lượng thượng nguồn 59

Biểu đồ 33: Diễn tiến Lưu lượng thường nguồn qua các năm 60

Biểu đồ 34: Biểu đồ histogram của Lưu lượng thượng nguồn 61

Biểu đồ 35: Biểu đồ scatter và biểu đồ tương quan tuyến tính giữa 62

Biểu đồ 36: Bảng thống kê mô tả và biểu đồ boxplot của Số đập thủy điện 64

Biểu đồ 37: Diễn tiến Số đập thủy điện qua các năm 65

Biểu đồ 38: Biểu đồ histogram của Số đập thủy điện 66

Biểu đồ 39: Biểu đồ scatter và biểu đồ tương quan tuyến tính giữa 66

Biểu đồ 40: Biểu đồ xương cá (Biểu đồ nguyên nhân – kết quả) của 73

Biểu đồ 41: Biểu đồ Trade-off Analysis cho các kịch bản 116

DANH MỤC CÁC HÌNH VẼ

Hình 1: Miêu tả quá trình xâm nhập mặn 14

Hình 2: Bản đồ dự báo phân bố độ mặn khu vực ĐBSCL 15

Hình 3: Bản đồ phân tần địa hình ĐBSCL 16

Hình 4 Các yếu tố ảnh hưởng đến xâm nhập mặn ở ĐBSCL 17

Hình 5 Các kịch bản ảnh hưởng của nước biển dân tại ĐBSCL 18

Hình 6: Bản đồ hệ thống sông ngòi ĐBSCL 22

Hình 7: Hồ chứa thủy điện Cảnh Hồng (Trung Quốc) 28

Hình 8: Sơ đồ mối quan hệ giữa các yếu tố có liên quan đến sự ảnh hưởng của 33

Hình 9: Kết quả từ Excel tương quan giữa Độ cao mực nước biển và Độ mặn trung bình S 41

Hình 10: Kết quả từ Excel tương quan giữa Khai thác nước ngầm và Độ mặn trung bình S 46

Hình 11: Kết quả từ Excel tương quan giữa Nhiệt độ trung bình và Độ mặn trung bình S 50

Hình 12: Kết quả từ Excel tương quan giữa Lượng mưa trung bình và Độ mặn trung bình S 55

Hình 13: Kết quả từ Excel tương quan giữa Lượng bốc hơi và Độ mặn trung bình S 59

Hình 14: Kết quả từ Excel tương quan giữa Lưu lượng thượng nguồn và Độ mặn trung bình S 63

Hình 15: Kết quả từ Excel tương quan giữa Số đập thủy điện và Độ mặn trung bình S 68

Hình 16: Bảng ma trận hệ số tương quan 69

Hình 17: Cây phân cấp AHP 78

Hình 18: Các mức độ quan trọng tương đối của đánh giá AHP 79

Hình 19: Mẫu phiếu khảo sát đánh giá AHP 79

Hình 20: Cây phân cấp AHP 97

Hình 21: Ma trận so sánh theo cặp 98

Hình 22: Ma trận trọng số 99

Hình 23: Kết quả phân tích AHP 100

Hình 24: Ma trận hệ số tương quan 101

Hình 25: Các hệ sinh thái hay các quá trình tự nhiên bị phá hủy hoặc suy thoái 124

Hình 26:Bảng phân tích dịch vụ hệ sinh thái 125

MỞ ĐẦU

1 ĐỀ TÀI NGHIÊN CỨU:

Trong những năm gần đây Biến đổi khí hậu đã và đang gây ra những ảnh hưởng tiêu cựcđến cuộc sống của con người trên mọi lĩnh vực, cả về môi trường và kinh tế - xã hội, trở thànhthách thức lớn nhất đối với nhân loại trong thế kỷ 21 Biến đổi khí hậu kéo theo thiên tai và hiệntượng thời tiết cực đoan đã làm tăng tần suất và cường độ bão, nhiệt độ Trái Đất ấm dần lên vànước biển dâng, ngập lụt, hạn hán, mưa lớn, sạt lở đất, xâm nhập mặn…Trong bối cảnh biến đổikhí hậu, nước biển dâng sẽ mở rộng vùng nước biển xâm nhập vào đất liền và thu hẹp diện tíchvùng nước ngọt dẫn đến xâm nhập mặn tại nhiều địa phương trên cả nước

Đồng bằng Sông Cửu Long là khu vực chịu ảnh hưởng nặng nề nhất do xâm nhập mặn gây

ra tại nước ta Thông qua chuỗi số liệu quan trắc khí tượng thủy văn từ 1975 tới nay cho thấy:thời gian, cường độ và tần suất hoạt động của các cực đoan khí hậu-thời tiết ở Việt Nam nói riêng

và khu vực Đông Nam Á nói chung (trong đó có lưu vực sông Mekong) có xu thế gia tăng rõ rệt

Do đó mức độ hạn, kiệt, mặn trên ĐBSCL cũng ngày càng gia tăng cả về thời gian, không gian,cường độ, tần suất cùng với những thiệt hại kinh tế môi trường kèm theo do chúng gây ra Tạimột số tỉnh ven biển ĐBSCL, nước biển xâm nhập sâu vào các sông rạch khiến các dòng sông bịnhiễm mặn sớm, gây ảnh hưởng nghiêm trọng đến đời sống người dân và hoạt động nông nghiệp.Những năm gần đây, diễn biến xâm nhập mặn ở ĐBSCL ngày càng phức tạp và bất thường, cónăm đến sớm năm đến muộn so với cùng kỳ nhiều năm Đặc biệt, những tháng đầu năm 2016,diễn biến xâm nhập mặn tại ĐBSCL được đánh giá nặng nề nhất trong nhiều năm qua và dự báocòn diễn biến xấu hơn trong những năm tiếp theo

Trước diễn biến ngày càng khó lường và hậu quả nghiêm trọng của vấn đề xâm nhập mặn,

nhóm chúng em chọn đề tài “Nghiên cứu vấn đề xâm nhập mặn tại Đồng bằng Sông Cửu Long” để tìm hiểu rõ hơn về Biến đổi khí hậu và những nguyên nhân dẫn đến tình trạng xâm

nhập mặn để từ đó đề ra những phương án nhằm giảm thiểu xâm nhập mặn tại khu vực nghiên

2 ĐỐI TƯỢNG NGHIÊN CỨU

Đề tài nghiên cứu về Tình hình xâm nhập mặn tại Đồng bằng Sông Cửu Long ảnh hưởng

do Biến đổi khí hậu

3 PHẠM VI NGHIÊN CỨU

Khu vực Đồng bằng Sông Cửu Long với diện tích 40.577 km² giai đoạn từ năm 2000 đếnnăm 2020

4 MỤC TIÊU NGHIÊN CỨU

- Làm rõ mối quan hệ giữa các yếu tố tự nhiên và nhân tạo liên quan đến vấn đề xâm nhậpmặn

- Đề xuất và phân tích các biện pháp cụ thể nhằm giảm thiểu xâm nhập mặn

5 PHƯƠNG PHÁP NGHIÊN CỨU

- Phương pháp phân tích nguyên nhân - kết quả

- Thống kê mô tả

- Phân tích đa tiêu chí

- Phân tích tương quan

- Kiểm định

- Phân tích ma trận SWOT Analysis

- Phương pháp hỗ trợ ra quyết định What-if Analysis

- Phương pháp đánh đổi Trade-off Analysis

- Phân tích AHP

- Phân tích Lợi ích - Chi phí BCA

- Phân tích dịch vụ hệ sinh thái

- Phân tích giá trị tổn thất của việc thất thoát các HST

6 KẾT CẤU CỦA ĐỀ TÀI

Đề tài ngoài phần mở đầu, kết luận và tài liệu tham khảo thì nội dung đề tài có kết cấuthành 4 chương như sau:

Chương 2: Phân tích quyết định và đề xuất giải pháp.

Chương 3: Tóm tắt nội dung và kết quả

Chương 4: Tổng kết

7 Ý NGHĨA CỦA ĐỀ TÀI

Làm rõ cơ chế xâm nhập mặn tại Đồng bằng sông Cửu Long bị ảnh hưởng bởi Biến đổikhí hậu, đề xuất được các biện pháp khoa học và thích hợp từ đó kết quả nghiên cứu có thể triểnkhai áp dụng tại khu vực bị ảnh hưởng bởi xâm nhập mặn

CHƯƠNG 1: TÍNH CẤP THIẾT CỦA ĐỀ TÀI XÂM NHẬP MẶN

TẠI ĐỒNG BẰNG SÔNG CỬU LONG 1.1 HIỆN TRẠNG VẤN ĐỀ XÂM NHẬP MẶN TẠI ĐBSCL:

1.1.1 Khái niệm:

Xâm nhập mặn là quá trình thay thế nước ngọt trong các tầng chứa nước ở ven biển bằngnước mặn do sự dịch chuyển của khối nước mặn vào tầng nước ngọt

Hình 1: Miêu tả quá trình xâm nhập mặn

Biến đổi khí hậu (BĐKH) là sự thay đổi của khí hậu và của những thành phần liên quangồm đại dương, đất đai, bề mặt Trái đất, và băng quyển

1.1.2 Diễn biến xâm nhập mặn tại Việt Nam do BĐKH:

Từ 1961 – 2003:

- Nhiệt độ đại dương toàn cầu tăng 0.1oC

- Mực nước trung bình của các đại dương tăng lên ~ 1.8 mm/năm

- Xâm nhập mặn diễn ra tại hầu hết các địa phương ven biển, gây ảnh hưởng lớn đến hoạtđộng

- Hai đồng bằng rộng lớn của Việt Nam là Đồng bằng sông Hồng và Đồng bằng sông CửuLong là những nơi chịu ảnh hưởng lớn nhất

Hình 2: Bản đồ dự báo phân bố độ mặn khu vực ĐBSCL

- Cây lương thực: lúa, bắp, đậu, cam, quýt là nhóm cây trồng chịu mặn yếu (tối đa 2g/L)

- Tiêu chuẩn QC 01 được đưa ra bởi Bộ Y Tế, độ mặn cho phép trong nước sinh hoạt, ănuống ở khu vực bình thường là 0.25g/L trở xuống là an toàn; riêng khu vực ven biển và hảiđảo là 0.3g/L trở xuống.

- Hơn phân nửa diện tích ĐBSCL không đủ nước dùng cho nông nghiệp, và thậm chí là sinhhoạt

Ngoài ra xâm nhập mặn sẽ làm đất bị nhiễm phèn, nhiễm mặn, sạt lở không canh tác được,thiệt hại về người và của

1.1.3 Các yếu tố ảnh hưởng đến xâm nhập mặn ở ĐBSCL:

Hình 4 Các yếu tố ảnh hưởng đến xâm nhập mặn ở ĐBSCL

1.1.3.1 Các yếu tố tự nhiên do biến đổi khí hậu:

a) Ảnh hưởng nóng lên toàn cầu:

 Ảnh hưởng đến nước biển dâng

Hình 5 Các kịch bản ảnh hưởng của nước biển dân tại ĐBSCL

- ĐBSCL có 3 mặt đều là biển bao gồm rất nhiều quần đảo và đảo, với mặt bờ biển dài73.2km

- Địa hình ĐBSCL thấp và phẳng, cao độ bình quân là +1m so với mực nước biển, ít đồi núitrừ một số ở vùng Thất Sơn (An Giang), độ dốc bình quân 1 cm/km (1/100.000), có nhữngvùng trũng như vùng Đồng Tháp Mười, vùng từ giác Long Xuyên - Hà Tiên và một sốvùng trũng nhỏ ở U Minh ĐBSCL chịu ảnh hưởng của thuỷ triều và xâm nhập mặn hàngnăm với diện tích nhiễm mặn lên tới 1,7 triệu ha

- Gần 0,8 triệu ha nằm dọc theo bờ biển, các chủ yếu ở Bến Tre, Bạc Liêu, Cà Mau và KiênGiang và một số huyện ở Sóc Trăng, Trà Vinh là những vùng bị ảnh hưởng của xâm nhậpmặn nghiêm trọng nhất

- Nếu mực nước biển dâng cao 1 m, khoảng 40% diện tích đồng bằng sông Cửu Long, 11%diện tích đồng bằng sông Hồng và 3% diện tích của các tỉnh khác thuộc vùng ven biển bịngập, TP.HCM bị ngập trên 20% diện tích; khoảng 10-12% dân số nước ta bị ảnh hưởngtrực tiếp và tổn thất khoảng 10% GDP

Long lại thấp nên càng dễ dàng bị xâm nhập mặn, và nhất là khi mà lưu lượng dòng chảy từthượng nguồn sông Mekong giảm thấp.

 Ảnh hưởng đến chế nhiệt độ và lượng mưa

Biểu đồ 1: Biểu đồ lượng mưa và nhiệt độ trung bình qua các năm

Trong 1 năm, biên độ nhiệt và lượng mưa trung bình mỗi tháng của ĐBSCL có sự thay đổi

rõ rệt:

- Các tháng mùa mưa : tháng 5 đến tháng 10: lượng mưa vượt 1500mm

- Các tháng mùa khô: tháng 11 đến tháng 4: lượng mưa thấp có tháng thậm chí gần như hạnhán (tháng 2)

Ngược lại, nhiệt độ trung bình vào mùa nắng lại khá cao, khoảng 28.5 độ C, góp phần đẩymạnh sự bốc hơi nước, gây ra tình trạng thiếu nước ngọt vào mùa khô và lũ lụt vào mùa mưa

Nếu như không có BĐKH làm lưu lượng nước đổ về ĐBSCL bị suy giảm, thì việc nhiệt độcao vào mùa khô cũng không ảnh hưởng nhiều đến xnm, vì biên độ nhiệt không cao

Biểu đồ 2: Biểu đồ lượng bốc hơi 4 năm tại 4 địa điểm chính vùng ĐBSCL

Nhìn chung, các tỉnh đại diện ĐBSCL đều có lượng bốc hơi cao, trên 800mm/năm, và bốchơi không đồng đều qua các năm

Cà Mau có biên độ lượng bốc hơi cao nhất, Cần thơ biên độ ổn định nhất Điều này cũng làbằng chứng cho thấy Cà Mau hiện tại là tỉnh bị thiệt hại về xâm nhập mặn lớn nhất ở Đồng bằngSông Cửu Long Do biến đổi khí hậu khiến nhiệt độ trái đất nóng lên mà sự bốc hơi cũng tăngtheo, khiến tình trạng xnm ngày càng trầm trọng hơn

Đầu năm 2016, trong thời gian khu vực đang trải qua mùa khô, lượng nước con sông này

đã giảm sút nghiêm trọng gây nên tình trạng thiếu nước, dẫn đến hiện tượng nước mặn thâm nhậpngược dòng lên khu vực Nguyên nhân chính là do hiện tượng El Nino và các đập ngăn dòngMekong Dựa vào biểu đồ, ngoại trừ năm 2016 là năm có độ mặn lịch sử do đợt hạn hán nghiêmtrọng nhất 100 qua

Có thể thấy, năm 2019 so với cùng kỳ các năm còn lại, độ mặn đo tại các trạm đa phần códấu hiệu tăng dần về độ mặn, và hầu hết độ mặn ở trên mức ranh mặn 4g/L và trên rất nhiều lần

so với tiêu chuẩn nước sạch 0.25g/L của bộ y tế

Vì thế, trong tương lai, chúng ta cần phải cố gắng rất nhiều để có thể đề xuất ra các biệnpháp ứng phó với biến đổi khí hậu Chỉ có vậy thì đời sống của người dân ĐBSCL mới được đảmbảo

Kết luận: Trong những năm lượng mưa mùa khô dồi dào, nguồn sinh thủy lớn, lượng

nước ngọt trong hệ thống kênh, rạch và đồng ruộng lớn sẽ hạn chế mặn xâm nhập vào trong hệthống kênh, rạch nội động; nhưng nếu bốc hơi lớn thì sẽ làm cạn kiệt nguồn nước ngọt, tạo điềukiện thuận lợi để mặn xâm nhập vào trong nội đồng

b) Dòng chảy thượng nguồn:

Đồng bằng sông Cửu Long có hệ thống sông ngòi, kênh rạch khá phong phú, bao gồm hệthống sông thiên nhiên và kênh đào: Sông Tiền và sông Hậu, hệ thống sông Vàm Cỏ, hệ thốngsông Cái Lớn-Cái Bé, sông Giang Thành và các sông nhánh khác.

Tuy vậy, trong những năm gần đây, do BĐKH mà lưu lượng dòng chảy thượng nguồn ởĐBSCL đang dần chạm đến mức thấp kỷ lục

Biểu đồ 4: Dòng chảy trên dòng chính sông Mekong (xét tại trạm Kratie)

 Xét sự ảnh hưởng của dòng chảy sông Mekong đối với hạn mặn tại ĐBSCL vào mùa cạngiai đoạn từ trước đến năm 2016:

Trong thời gian từ 11/2015- 10/1/2016, dòng chảy tại Kratie rất thấp, đạt con số nhỏ nhấttrong hơn 30 năm Lưu lượng trung bình thời kỳ này nhỏ hơn giá trị trung bình của chuỗi số liệu1985-2015 khoảng gần 40%, và còn nhỏ hơn khoảng 28% so với năm hạn lịch sử 1998-1999 Từsau 15/1/2016, lưu lượng tại Kratie đạt mức trung bình nhiều năm

Biểu đồ 5: Dòng chảy trên dòng chính sông Mekong (xét tại trạm PrekDam)

 Xét sự ảnh hưởng của dòng chảy sông Mekong đối với hạn mặn tại ĐBSCL vào mùa cạngiai đoạn từ trước đến năm 2016:

o Do lưu lượng trên dòng chính sông Mekong quá thấp, thời điểm nước chảy ngược

từ sông chính vào hồ để trữ nước cho Biển Hồ muộn so với trung bình nhiều nămnhiều ngày, toàn bộ thời gian trữ nước cho hồ giảm gần ½ so với trung bình nhiềunăm Hệ quả, lượng nước dòng chính cấp cho Biển Hồ giảm 33%

o Mặt khác, dòng chảy bổ sung cho hồ từ lưu vực Biển Hồ giảm mạnh, do ít mưa.Điều này dẫn đến, lượng nước cấp cho ĐBSCL trong thời kỳ từ tháng 10/2015 đến22/2/2016 giảm đến 53% so với cùng kỳ trung bình nhiều năm Đặc biệt trong cáctháng 11/2015-2/2016, dòng chảy giảm 56%

Biểu đồ 6: Biểu đồ độ mặn lớn nhất đến cuối thnasg 2/2016

so với TBNN và CK Lịch sử khu vực ĐBSCL

 So sánh độ mặn của năm 2016 so với trung bình nhiều năm và độ mặn cao nhất trong lịch

sử, ta có thể nhìn ra độ mặn 2016 cao vượt trội so với trung bình nhiều năm và max lịch

sử, trở thành năm hạn mặn lịch sử tiếp theo trong vòng 100 năm

 Năm 2015, do ảnh hưởng của hiện tượng El-Nino nên mùa mưa đến trễ nhưng lại kết thúcsớm, tổng lượng mưa trên lưu vực thiếu hụt so với trung bình nhiều năm từ 20-50% Vìvậy, mùa lũ năm 2015 thuộc năm lũ nhỏ dần, dẫn đến dòng chảy trong mùa khô từ thượnglưu chảy về ĐBSCL xuống ở mức cực thấp

Kết luận: Do BĐKH, dòng chảy từ thượng nguồn đang tiến dần về mức cực thấp mà

lượng dòng chảy từ biển vào lại không đổi sẽ góp phần gây ra xâm nhập mặn

1.1.3.2 Các yếu tố nhân tạo:

Ngoài các yếu tố tự nhiên, chúng ta không thể không đề cập đến yếu tố nhân tạo trong cácyếu tố ảnh hưởng đến xâm nhập mặn do BĐKH

a) Nhu cầu sử dụng và khai thác nước ngọt:

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

Năm 2005 Năm 2015 Năm 2020

Chênh lệch nhu cầu giữa năm 2005 và 2020 khá lớn, từ năm 2015 đến 2020 có tăng nhẹ

Do biến đổi khí hậu nên xâm nhập mặn chủ yếu xảy ra vào mùa khô, gây nên tình trạng thiếunước ngọt dẫn đến nhu cầu sử dụng nước ngọt tăng

b) Đập thủy điện:

1966-1990 1990-2000 2000-2010 2010-20200

và chính sách bảo vệ ĐDSH; quan tâm bảo vệ môi trường,… cùng với sự xuất hiện của phongđiện, đã hạn chế các mặt có hại của việc xây dựng đập thủy điện và giảm số lượng đập đưa vào sửdụng nhưng vẫn đảm bảo được nguồn điện cung cấp cho người dân

Hình 7: Hồ chứa thủy điện Cảnh Hồng (Trung Quốc)

Cùng với vấn đề sử dụng nước, lưu lượng nước sông Mekong trong tương lai cũng khôngđược dự báo một cách chắc chắn Hiện tại trên thượng lưu sông Mekong có một số công trình đậpthủy điện đã hoàn thành, một số công trình đập thủy điện đang thi công, một số kế hoạch pháttriển trong đó có cả đập thủy điện tại thượng lưu (riêng ở lãnh thổ Trung Quốc đã có 4 đập quy

mô lớn và hơn 10 đập được quy hoạch trong năm 2011) Các công trình ở thượng lưu, đặc biệt làcác công trình đập thủy điện, sẽ có ảnh hưởng lớn đến chế độ dòng chảy ở hạ lưu sông Mekong.Các đập sẽ trữ nước mùa lũ và xả nước trong mùa khô để sản xuất điện và do đó sẽ làm tăng lưulượng hạ lưu sông Mekong Tuy nhiên, nếu khu vực đập thượng nguồn giảm xả nước, lại gây ratình trạng thiếu hụt lưu lượng, sẽ góp phần cổ vũ cho xâm nhập mặn

Ví dụ: Hồ chứa thủy điện Cảnh Hồng (Trung Quốc) giảm xả nước xuống hạ du từ ngày 5đến 24-1-2021 để bảo trì lưới điện Do đó, lưu lượng xả ra từ thủy điện Cảnh Hồng xuống hạ lưu

sẽ giảm gần 50%, xuống còn khoảng 1.000 m3/giây Việc giảm lưu lượng xả sẽ làm gia tăng tìnhtrạng xâm nhập mặn ở đồng bằng sông Cửu Long (ĐBSCL) từ cuối tháng 1 đến cuối tháng 2

Kết luận: Việc phát triển nông nghiệp và nhu cầu con người ngày càng tăng đã vô tình

góp phần gây ra BĐKH và ảnh hưởng trực tiếp đến tình trạng xâm nhập mặn ở vùng ĐBSCL.Trong tương lai, cần quy hoạch có kế hoạch đối với lưu lượng sông Mekong và cả ứng phó vớiBĐKH

1.2 TRÌNH BÀY VÀ MÔ TẢ DỮ LIỆU:

1.2.1 Cơ sở lý thuyết kiểm định thống kê:

a) Giá trị trung bình (Mean):

Giá trị trung bình là bình quân toán học đơn giản của một tập hợp gồm hai hoặc nhiều số.Giá trị trung bình của một tập hợp số đã cho có thể được tính theo nhiều cách, gồm có phươngpháp trung bình số học sử dụng tổng các số trong chuỗi và phương pháp trung bình hình học Tuynhiên, tất cả các phương pháp chính để tính trung bình đơn giản của một chuỗi số bình thườngđều tạo ra kết quả xấp xỉ bằng nhau

b) Giá trị trung vị (Median):

Trung vị là số nằm giữa trong một tập dữ liệu có các số được sắp xếp Để xác định giá trịtrung vị trong một chuỗi số, trước tiên các số phải được sắp xếp theo thứ tự giá trị từ thấp nhấtđến cao nhất hoặc cao nhất đến thấp nhất Trung vị có thể được sử dụng để xác định giá trị trungbình gần đúng hoặc giá trị trung bình, tuy nhiên không được nhầm lẫn trung vị với giá trị trungbình thực tế

c) Yếu vị (Mode):

Một tập hợp các số có thể có một hoặc nhiều hơn một yếu vị hoặc không có yếu vị nào cả

như vậy.

d) Độ lệch chuẩn (Standard deviation):

Độ lệch chuẩn được tính là căn bậc hai của phương sai, được tính bằng cách xác định sựchênh lệch giữa mỗi điểm dữ liệu so với giá trị trung bình Nếu một điểm dữ liệu nằm xa giá trịtrung bình, điểm đó có độ lệch cao trong tập dữ liệu, dữ liệu càng có độ dàn trải rộng thì độ lệchchuẩn càng cao

e) Phương sai mẫu (Sample Variance):

Phương sai được tính bằng cách xác định giá trị chênh lệch giữa mỗi số trong tập dữ liệuvới giá trị trung bình, sau đó bình phương các chênh lệch nhằm cho chúng mang giá trị dương vàkhông triệt tiêu lẫn nhau Cuối cùng chia tổng số lượng quan sát trong tập dữ liệu

f) Độ lệch (Skewness):

Nếu phân hối hình chuông bị dịch chuyển sang trái hoặc sang phải, nó được cho là bị lệch

Độ lệch được coi là một đại diện cho mức độ khác biệt của một phân phối nhất định so với phânphối chuẩn Một phân phối chuẩn có độ lệch bằng 0, trong khi phân phối phân phối xác suất logachuẩn sẽ có lệch về phía bên phải

g) Độ nhọn (Kurtosis):

Giống như độ lệch, độ nhọn là một biện pháp thống kê được sử dụng để mô tả các phânphối Trong khi độ lệch phân biệt các giá trị cực trị ở một đuôi so với đuôi kia, thì độ nhọn đolường các giá trị cực trị ở một trong hai đuôi

 Kiểm định phương trình hồi quy (kiểm định F):

Giả thiết H0: “Phương trình hồi quy tuyến tính không thích hợp”

Giả thiết H1: “Phương trình hồi quy tuyến tính thích hợp”

Miền bác bỏ H0: Trắc nghiệm F > F α,1,n-k

Ở đề tài này, ta có ⍺ = 0.05, n - k = 20 -2 = 18 nên F α,1,n-k = F 0.05,1,18= 4.4139

 Ước lượng đường hồi quy:

Phương trình hồi quy tuyến tính:

 Kiểm định hệ số a, b có nghĩa thống kê (a là hệ số tự do, b là hệ số góc)

Giả thiết H0: a, b không có ý nghĩa thống kê

Giả thiết H1: a, b có ý nghĩa thống kê

H0: ρ = 0 đối với H1: ρ ≠ 0 (hoặc ρ > 0 hoặc ρ < 0)

Người ta chứng minh được rằng với giả thiết H0, phân phối mẫu của R đối xứng; từ đó, thống kê

tqs = R√ n−2

1−R2~ Student (n – 2)Trắc nghiệm t được dùng trong trường hợp này

 Kiểm định hệ số hồi quy:

Giả thiết H0: Hệ số hồi quy không có ý nghĩa (=0)

Giả thiết H1: Hệ số hồi quy có ý nghĩa (≠0)

Miền bác bỏ H0: |t qs|>t ∝, n−2

Ở đây, do bảng số liệu đều có n = 20 và α=0.05 nên |t qs|>2.1009 → Bác bỏ H0

1.2.2 Thống kê số liệu quan trắc các năm ở Đồng bằng sông Cửu Long:

 Sơ đồ về mối quan hệ giữa các yếu tố có liên quan đến sự ảnh hưởng của Biến đổi khí hậuđến Xâm nhập mặn:

Hình 8: Sơ đồ mối quan hệ giữa các yếu tố có liên quan đến sự ảnh hưởng của

Biến đổi khí hậu đến Xâm nhập mặn

Từ sơ đồ trên, nhóm liệt kê ra 7 yếu tố sau đây để khảo sát tương quan ảnh hưởng trực tiếpđến xâm nhập mặn tại Đồng bằng sông Cửu Long:

- Lượng mưa trung bình năm (mm)

- Nhiệt độ trung bình năm (oC)

- Lượng bốc hơi trung bình năm (mm)

- Mực nước biển (mm)

- Lưu lượng nước thượng nguồn (m3/s)

- Khai thác mực nước ngầm (m3/ngày)

● Bảng thống kê số liệu quan trắc các năm ở Đồng bằng sông Cửu Long

Bảng 1: Bảng số liệu các biến liên quan đến xâm nhập mặn giai đoạn 2000-2020 tại Đồng bằng Sông Cửu Long

- Lưu lượng thượng nguồn (m3/s): Theo số liệu của Văn phòng thường trực, Ban chỉ đạotrung ương về Phòng chống thiên tai.

- Số đập thủy điện hoạt động (cái): Theo Bách khoa toàn thư

1.2.3 Thống kê mô tả, hệ số tương quan và kiểm định của các biến:

1.2.3.1 Độ mặn trung bình S:

a) Thống kê mô tả:

Biểu đồ 9: Bảng thống kê mô tả và biểu đồ boxplot của Độ mặn trung bình S

Nhận xét:

- Độ mặn có giá trị trung bình là 8.82 (g/L), giá trị trung vị là 5.65 (g/L)

- Độ mặn có xu hướng nghiêng về phía dưới giá trị trung bình

- Khoảng số phân tử = Q3 - Q1 = 10.2 - 2.32 = 7.88

- Khoảng cách giữa giá trị lớn nhất và nhỏ nhất là: 30.6 - 1.45 = 29.15

Nhìn chung độ mặn trung bình có sự dao động lớn, chênh lệch nhiều giữa các năm

Biểu đồ 10: Diễn tiến của độ mặn trung bình S qua các năm

Nhận xét: Từ năm 2000 đến năm 2020 độ mặn trung bình S(g/L) có sự thay đổi:

- Từ năm 2000 – 2003: độ mặn có xu hướng giảm dần từ 2.32 g/L đến 1.45 g/L

- Từ năm 2003 – 2004: độ mặn có xu hướng tăng dần từ 1.45 g/L đến 2.85 g/L

- Từ năm 2004 – 2005: độ mặn có xu hướng giảm dần từ 2.85 g/L đến 2.17 g/L

- Từ năm 2005 – 2006: độ mặn có xu hướng tăng dần từ 2.17 g/L đến 9.25 g/L

- Từ năm 2006 - 2009: độ mặn có xu hướng giảm dần từ 9.25 g/L đến 3.85 g/L

- Từ năm 2009 – 2010: độ mặn có xu hướng tăng dần từ 3.85 g/L đến 7.8 g/L

- Từ năm 2010 – 2012: độ mặn có xu hướng giảm dần từ 7.8 g/L đến 2 g/L

- Từ năm 2012 – 2016: độ mặn có xu hướng tăng dần từ 2 g/L đến 30.6 g/L

- Từ năm 2016 – 2017: độ mặn có xu hướng giảm dần từ 30.6 g/L đến 12.4 g/L

- Từ năm 2017 – 2019: độ mặn có xu hướng tăng dần từ 12.4 g/L đến 28.3 g/L

- Từ năm 2019 – 2020: độ mặn có xu hướng giảm dần từ 28.3 g/L đến 25.6 g/L

Nhìn chung, từ năm 2000 – 2020 độ mặn có xu hướng tăng từ thấp nhất 1.45 g/L (năm2003) đến cao nhất 30.6 g/L (năm 2016), sự tăng giảm độ mặn không đồng đều và có phần tănggiảm xen kẽ Đặc biệt, năm 2016 là năm có hạn mặn lịch sử trong hơn 100 năm qua, gây ra nhiềuthiệt hại nặng nề về con người và cả lao động sản xuất

Biểu đồ 11: Biểu đồ histogram của Độ mặn trung bình S

Nhận xét:

- Trong khoảng từ 0 đến 7 (g/L) : có 11 giá trị

- Trong khoảng từ 7 đến 14 (g/L): có 6 giá trị

- Trong khoảng từ 14 đến 21 (g/L): có 1 giá trị

- Trong khoảng từ 21 đến 28 (g/L): có 1 giá trị

- Trong khoảng từ 27 đến 35 (g/L): có 2 giá trị

Dữ liệu có xu hướng tập hợp phần lớn trong khoảng từ 0 - 7 (g/L) về phía bên trái, đây làdạng phân phối lệch trái

1.2.3.2 Độ cao mực nước biển

a) Thống kê mô tả:

Biểu đồ 12: Bảng thống kê mô tả và biểu đồ boxplot của Độ cao mực nước biển

Nhận xét:

- Độ cao mực nước biển có giá trị trung bình là 67.86 (mm), giá trị trung vị là 69 (mm)

- Độ cao mực nước biển chênh lệch ít với mức trung bình, có xu hướng nghiêng về phía trêngiá trị trung bình

- Khoảng số phân tử = Q3 - Q1 = 77.5 - 57 = 20.5

- Khoảng cách giữa giá trị lớn nhất và nhỏ nhất là: 86 - 45 = 41

Nhìn chung Độ cao mực nước biển có sự chênh lệch giữa các năm

Biểu đồ 13: Diễn tiến Độ cao mực nước biển qua các năm

Nhận xét: Từ năm 2000 đến năm 2020 độ cao mực nước biển(mm) ngày càng tăng dần:

- Từ năm 2000 - 2011, độ cao mực nước biển có xu hướng tăng từ 45mm lên 72mm

- Từ năm 2011- 2012, độ cao mực nước biển có xu hướng giảm từ 72mm xuống 71.5mm

- Từ năm 2012 - 2013, độ cao mực nước biển có xu hướng tăng từ 71.5mm lên 74.5mm

- Từ năm 2013 - 2014, độ cao mực nước biển có xu hướng giảm từ 74.5mm xuống 72mm

- Từ năm 2014 - 2017, độ cao mực nước biển có xu hướng tăng từ 72mm lên 84mm

- Từ năm 2017 - 2018, độ cao mực nước biển có xu hướng giảm từ 84mm xuống 81,5mm

- Từ năm 2018 - 2020, độ cao mực nước biển có xu hướng tăng từ 81.5mm lên 86mm.Nhìn chung, từ năm 2000- 2020, độ cao mực nước biển có xu hướng tăng từ thấp nhất là45mm (năm 2000) lên cao nhất là 86mm ( năm 2020) Mực nước biển ngày càng tăng cao làmdiện tích nước ngọt ngày càng bị thu hẹp, mở rộng hơn vùng xâm nhập mặn gây ảnh hưởng trựctiếp đến cuộc sống sinh hoạt cũng như sản xuất của người dân vùng Đồng bằng sông Cửu Long

Biểu đồ 14: Biểu đồ histogram của Độ cao mực nước biển

Nhận xét:

- Trong khoảng từ 40 đến 50 dữ liệu có 2 giá trị

- Trong khoảng từ 50 đến 60 dữ liệu có 4 giá trị

- Trong khoảng từ 60 đến 70 dữ liệu có 5 giá trị

- Trong khoảng từ 70 đến 80 dữ liệu có 5 giá trị

- Trong khoảng từ 80 đến 90 dữ liệu có 5 giá trị

Dữ liệu có xu hướng tập hợp phần lớn trong khoảng từ 60 - 90 (mm) về phía bên phải, đây

là dạng phân phối lệch phải

b) Tương quan giữa Độ cao mực nước biển đối với Độ mặn trung bình S: