BÁO CÁO HỌC THUẬT ẢNH HƯỞNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN XÂM NHẬP MẶN TẦNG CHỨA NƯỚC HOLOCENE VÙNG BỜ BIỂN GIO LINH- QUẢNG TRỊ

ĐẠI HỌC MỎ - ĐỊA CHẤT KHOA MÔI TRƯỜNG BỘ MÔN ĐỊA SINH THÁI VÀ CÔNG NGHỆ MÔI TRƯỜNG BÁO CÁO HỌC THUẬT ẢNH HƯỞNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN XÂM NHẬP MẶN TẦNG CHỨA NƯỚC HOLOCENE VÙNG BỜ BIỂN GIO LINH- QUẢNG TRỊ Người thực hiện: Th.S Đỗ Cao Cường Hà Nội, tháng…/2019 MỞ ĐẦU Trong nghiên cứu bổ cập có nhiều cơng cụ khác để tiến hành nghiên cứu, yếu tố tự nhiên đóng vai trị quan trọng, ảnh hưởng trực tiếp đến nguồn hình thành nước đất Trên sở kết nghiên cứu điều kiện địa hình, khí hậu, thủy văn, cấu trúc địa chất, nguồn gốc hình thành trầm tích, đặc điểm kiến tạo, vận động nước đất, thông số địa chất thủy văn trạng nhiễm mặn khu vực kết hợp với yếu tố nhân tạo diễn biến khai thác tốc độ thị hóa, vv… đánh giá khả dự báo bổ cập diễn biến xâm nhập mặn tầng chứa nước Holocene Xâm nhập mặn vào nước đất nhạt đặc biệt vùng ven biển nhà khoa học giới nghiên cứu từ lâu Quá trình xâm nhập mặn xuất dọc theo khu vực bờ biển gây sơng nước biển chảy vào sâu đất liền tác động thủy triều bão Sự xâm nhập mặn vào tầng chứa nước khai thác Việt Nam phát từ kỷ trước Những năm 60 kỷ trước, tăng công suất khai thác nước đất từ 80.000m3/ng lên 160.000m3/ng TP Hồ Chí Minh ranh giới mặn nhạt xâm nhập vào đến gần 20km, sau phải chuyển sang khai thác nước mặt Sau gần 20 năm ngừng khai thác nước đất, ranh giới mặn nhạt di chuyển đến khoảng 2km Đoàn Văn Cánh (1996) tiến hành đề tài nghiên cứu “Tài nguyên môi trường nước đất vùng Nam Định - Hà Nam” Đề tài tập trung nghiên cứu chi tiết trữ lượng chất lượng môi trường nước đất khu vực Nam Định - Hà Nam Đây nghiên cứu tập trung vào tài nguyên nước đất khu vực Các kết tranh tổng thể trạng dự báo tài nguyên nước đất tương lai Đặng Hữu Ơn (1996) có dự báo khả nhiễm mặn nước đất Bà Rịa - Vũng Tàu thông qua việc xác định hệ số rỗng hữu hiệu dựa sơ đồ phễu hạ thấp mực nước cơng trình đưa vào hoạt động Tác giả xác định vận tốc dòng thấm trung bình theo hướng từ Biển vào cơng trình từ tính thời gian nước mặn xâm nhập vào cơng trình khai thác Ngun nhân biến đối khí hậu phức tạp đa dạng, bao gồm nguyên nhân tự nhiên nguyên nhân người Vì vậy, để giảm thiểu tác hại cần chủ động phối hợp, đề giải pháp đồng mang tính đa ngành, đa lĩnh vực để đối phó thích ứng với biến đổi khí hậu Trên sở định hướng nghiên cứu đó, đề tài “Ảnh hưởng biến đổi khí hậu đến xâm nhập mặn tầng chứa nước Holocene vùng bờ biển Gio Linh, Quảng Trị” nhằm mục đích giải số vấn đề tồn nghiên cứu đánh giá ảnh hưởng điều kiện tự nhiên như: cấu trúc địa chất, địa chất thủy văn, nguồn bổ cập trạng, nhu cầu khai thác diễn biến xâm nhập mặn, nhằm phục vụ việc khai thác bền vững nguồn tài nguyên quý giá vùng Chương Mối quan hệ thủy nước nước mưa, nước mặt với nước đất 1.1 Xác định mối quan hệ thủy lực nước mặt với nước đất 1.1.1 Thí nghiệm đổ nước hố đào a Cơ sở lý thuyết: Thí nghiệm đổ nước vào hố đào có vịng chắn, khống chế cột nước áp lực khơng đổi 0,1 m; tiến hành quan trắc thấm xác định lưu lượng thấm ổn định; tính toán hệ số thấm đất theo định luật Darcy, dựa ba giả thiết sau đây: - Dòng thấm từ đáy hồ đào vòng chắn hướng thẳng xuống đất, khơng chảy tản phía, nên tiết diện thấm tiết diện vòng chắn - Với chiều cao cột nước áp lực 10 cm nhỏ, nên lưu lượng thấm đạt đến ổn định, chiều sâu nước thấm xuống đất nhỏ vậy, có gradient thủy lực tương ứng 1; vận tốc thấm hệ số thấm - Áp lực mao dẫn không đáng kể b Thiết bị, dụng cụ thí nghiệm Vịng chắn thép, hình trụ có đường kính 50 cm, cao từ 20 cm đến 25 cm, thành dày từ mm đến mm đầu vát sắt mép ngồi Hình Dụng cụ đổ nước hố đào Thiết bị cấp nước, hai thùng đo định chuẩn có thơng số kỹ thuật Thùng đo định chuẩn chế tạo thép cứng, có dạng hình trụ trịn thẳng đứng khơng có lồi lõm, đường kính 400 mm, chiều cao 800 mm Ở phần thùng đo có hai van: van nối thơng với ống đo mực nước (2) làm thủy tinh, gắn bảng thang đo chia vạch mm đặt thẳng đứng cố định với thành thùng đo; van nối với ống dẫn nước Hai thùng đo định chuẩn đặt thẳng đứng giá đỡ vững gắn chặt chẽ với đai thép cứng (5), ống dẫn nước hai thùng kết nối với van ba nhánh (3) Nhánh thứ ba van nối với ống cứng (4) để dẫn nước kết nối với đầu ống phao nước Khi thí nghiệm, van ba nhánh cho phép đóng nước thùng vơi hết nước, đồng thời mở nước thùng thứ hai đổ đầy nước trước đó, đảm bao cho thí nghiệm liên tục Các dụng cụ thơng thường khác đồng hồ bấm giây, đồng hồ giờ, dụng cụ đào hố; gạt phẳng; thước đo có chiều dài khoảng 50 cm đến 100 cm với thang chia mm; thùng chứa nước dự trữ; vật liệu đất sét mềm dẻo, vật liệu sỏi hạt mm đến 10 mm; nhiệt kế chia độ đến 50 0C, số đọc xác đến 0,5 0C Hình Sơ đồ thí nghiệm đổ nước hố đào c Kết thí nghiệm: Chỉnh lý số đọc mực nước thang đo thùng cấp nước thời điểm quan trắc; tính lưu lượng thấm với khoảng thời gian q trình thí nghiệm, theo cơng thức 1: Q=V t (1) đó: Q lưu lượng thấm, lít/phút; V lưu lượng nước tiêu tốn thấm (l), thời gian đo từ t1 đến t2, tính theo cơng thức: V = V1 - V2, với V1 lượng nước thùng cấp nước thời điểm đo t V2 lưu lượng thùng thời điểm đo t2, xác định theo biểu đồ quan hệ số đo thang đo mực nước lượng nước thùng t khoảng thời gian hai lần đo: t = (t2 - t1), s (giây) Lập biểu đồ quan hệ thể tích V, lưu lượng Q thời gian t để xác định trị số lưu lượng ổn định (Qc) Tính hệ số thấm đất, theo cơng thức 2: Kth = Qc F (2) đó: Qc lưu lượng thấm ổn định, cm3/s; F tiết diện thấm, tiết diện vòng chắn, cm 2; tính theo cơng thức 2: F=   D2 (2') với: D đường kính vịng chắn trong, m;  số Pi, lấy 3,14 Bảng Các thơng tin vị trí thí nghiệm V/2 (lít) F (m2) b 1869200 136,73 70,49 10,60 5,30 0,20 1,04 725145 1880886 60,42 5,12 2,56 0,20 1,34 H06- Hiền Lương 718674 1880828 102,00 24,44 14,98 7,49 0,20 0,02 H01- Lâm Xuân 725966 1869987 102,00 16,17 13,10 6,55 0,20 0,02 H04-Gio An 714569 1872317 102,00 23,76 17,48 8,74 0,20 0,02 H03_Rừng cao su 720075,3 1866359 102,00 26,50 15,40 7,70 0,20 0,04 TT Tên hố đào X HĐ05- Cửa Tùng 728784 HĐ02- Lâm Xuân Y T (phút) 95,00 T1 (phút) V (lít) Bảng Kết tính tốn Tên hố đào TT X Y HĐ05- Cửa Tùng 728784 1869200 HĐ02- Lâm Xuân 725145 1880886 H06- Hiền Lương 718674 1880828 H01- Lâm Xuân 725966 1869987 H04-Gio An 714569 1872317 H03_Rừng cao su 720075,3 1866359 Mô tả thành phần thạch học Cát trắng nguồn gốc trầm tích gió biển (vmQ23) Cát trắng nguồn gốc trầm tích gió biển (vmQ23) Sét, sét pha (am Q23) Sét, sét pha (am Q23) Đất đá bazan phong hóa Đất đá bazan phong hóa Hệ số thấm (m/ngày) Hệ số thấm kiến nghị (m/ngày) 0,592 0,561 0,53 0,022 0,0205 0,019 0,025 0,0345 0,044 Nhận xét: Kết thí nghiệm đổ nước hố đào toàn khu vực cho phép phân vùng bổ cập cho tầng chứa nước Holocen, kết làm hiệu chỉnh mơ hình Wetspa khu vực Lượng bổ cập tính trung bình khu vực trầm tích có nguồn gốc gió biển trung bình 204,76 mm/năm, khu vực trầm tích sét, sét pha lượng bổ cập trung bình 7,48 mm/năm, khu vực đất đá bazan phong hóa lượng cập trung bình 12,59 mm/năm 1.1.2 Thí nghiệm thấm seepage a Cơ sở lý thuyết Bộ dụng cụ thấm rỉ seepage áp dụng từ năm 1940 nhà khoa học Israelson Reeve nghiên cứu nước kênh tưới (1944) phát triển vào năm 1970 để xác định lượng thấm rỉ cho hồ chứa (McBridge Plannkuch, 1975; Lee, 1977; John Lock; 1977; Lee Cherry, 1978) Dụng cụ thấm seepage sử dụng nhiều nghiên cứu lượng thấm sông, vùng ven biển hay hồ chứa lớn (Lee Hynes (1978), Woessner (1984); Libelo MacIntyre (1994); Cey n.n.k (1998) Landon nnk (2001) v.v Phương pháp xác định thấm rỉ seepage xác định cách dùng dụng cụ hình trụ, phần phía khơng có đáy, phía nắp kín nối với túi nhựa qua van phần rìa Khi tiến hành thí nghiệm, ấn từ từ dụng cụ hình trụ xuống đáy sông phần nắp dụng cụ cách đáy khoảng 2cm Bộ dụng cụ thí nghiệm lắp nghiêng xoay phần nắp van phía cao giúp cho khơng khí dụng cụ dễ dàng ngồi qua van Đánh số dụng cụ thí nghiệm túi nhựa, đưa vào túi nhựa lượng nước cho trước trước nối túi nhựa vào van dụng cụ thí nghiệm Sau thời gian thí nghiệm t, đóng van nối túi nhựa với dụng cụ seepage xác định biến đổi thể tích nước túi nhựa Lượng nước thấm xác định theo công thức: Q (V1  V0 ) t.F đó: Q: lượng thấm V0: thể tích nước cho vào túi nhựa; V1: thể tích nước túi nhựa sau thí nghiệm; t: thời gian tiến hành thí nghiệm; F: diện tích mặt cắt dụng cụ Thể tích nước (V1-V0) có dấu (-) lượng nước đi, có dấu (+) lượng nước bổ xung, tương tự với Q có dấu (-) lượng nước thấm xuống đáy sông mang dấu (+) lượng nước cấp lên túi Hầu hết sai số trình thí nghiệm thấm rỉ đáy sơng seepage chủ yếu bị ảnh hưởng trình lắp đặt dụng cụ thấm dó vào đáy sơng Trong q trình lắp đặt, dụng cụ không điều chỉnh cho phần đặt van vị trí cao khơng khí tích tụ giảm độ thấm rỉ Các kết đo thấm rỉ bị giảm xuống phần nắp dụng cụ thí nghiệm bị ấn sát với đáy sông, hồ Trong trường hợp này, kết đo lượng nước khơng cịn xác Đây dạng thí nghiệm đơn giản, không tốn sử dụng để xác định trực tiếp lượng nước cung cấp cho sông Dạng thí nghiệm có ưu điểm khơng cần làm thí nghiệm xác định tính thấm lớp trầm tích đáy sơng Lượng nước thu sau thí nghiệm sử dụng để phân tích thành phần hóa học để nghiên cứu thành phần nguồn nước cung cấp cho nước mặt Phương pháp thí nghiệm xác định thấm rỉ nhà khoa học kiểm chứng thí nghiệm phịng (Israelson Reeve (1944), McBridge Plannkuch (1975); Lee Cherry (1978), Cey n.n.k (1998) Landon nnk (2001)) Thí nghiệm kiểm chứng độ xác tiến hành bể có kích thước lớn Trong thời gian gần đây, Cey nnk (1998) Landon nnk tiến hành thí nghiệm thùng thí nghiệm có kích thước dài x rộng x cao 1,5mx1,5x1,5m Mẫu đất sử dụng thí nghiệm lấy đáy sơng vị trí tiến hành thí nghiệm seepage cần so sánh đối chứng Lượng nước cấp cho phần nước mặt cấp cho phần mẫu đất qua đáy thùng thí nghiệm đo đầu tự động (dataloger) đồng hồ đo Thí nghiệm kiểm chứng mức độ xác thí nghiệm seepage bố trí hình vẽ Hình Mơ tả thí nghiệm seepage phịng b Ngun lý cấu tạo dụng cụ thí nghiệm Trong chục năm qua, nhà khoa học đưa thiết kế khác dụng cụ thấm seepage nhằm hạn chế tối đa sai số gây dụng cụ giúp thuận lợi q trình làm thí nghiệm Tuy nhiên, hầu hết dụng cụ có cấu tạo giống thiết kế mà Lee (1977) đưa ra, bao gồm phần phần hình trụ hình hộp chữ nhật (Paulsen nnk, 2001) có nắp phía khoan lỗ nhỏ để lắp đặt đầu nối với ống nhựa, đầu ống nhựa nối với túi nhựa, ống lắp van để đóng mở lắp đặt thí nghiệm 10 Ngồi việc mơ hình hố địa tầng vùng nghiên cứu phải xây dựng cập nhật liệu địa hình, ĐCTV, giá trị bổ cập bốc hơi, điều kiện biên, điều kiện khai thác b Dữ liệu địa hình Bề mặt địa hình xây dựng dựa thông tin từ đồ địa hình số hố gán thơng tin sở đồ địa hình tỉ lệ 1/50.000 Bản đồ địa hình khu vực nghiên cứu thể hình 13 Hình 13 Bản đồ địa hình vùng c Dữ liệu Địa chất thuỷ văn 34 Ranh giới tầng bình đồ mặt cắt xây dựng từ liệu lỗ khoan khảo sát ĐC - ĐCTV khu vực nghiên cứu Các thông số ĐCTV cần nhập cho lớp gồm: hệ số thấm (thẳng đứng Kz nằm ngang Kx-y), hệ số nhả nước đàn hồi µ* hệ số nhả nước trọng lực µ Phân chia tính thấm, hệ số nhả nước tầng chứa nước lớp ngăn cách hay thấm nước yếu từ số liệu thực tế thí nghiệm lỗ khoan điểm khảo sát Các số liệu hệ số thấm, hệ số nhả nước tần chứa nước lấy đề án giai đoạn trước như: (Nguyễn Xuân Dương, 1978 [18]; Nguyễn Văn Thế, 1984 [17]; Nguyễn Văn Long, 1986 [16]; Lê Quang Mạnh, 1990 [6]; Nguyễn Trường Giang, 1995 [15]; Khổng Văn Bê, 2003 [3])… Các số liệu đặc trưng điển hình thể bảng: Bảng Hệ số thấm số lỗ khoan vùng nghiên cứu Tên lỗ khoan LK604 LK605 LK904 LK901 LK604B LK608 LK610 LK404 LK405 LK407 LK410 LK413B Chiều sâu khai thác, tuổi địa chất 15 – 45, 𝑎𝑚𝑄21−2 – 23, 𝑎𝑚𝑄21−2 - 11, 𝑎𝑚𝑄21−2 31 – 40, 𝛽𝑄21 22 – 58, 𝑎𝑚𝑄12−3 40 – 59, 𝑎𝑚𝑄12−3 – 22, 𝑎𝑚𝑄12−3 18 – 56, 𝑎𝑚𝑄12−3 44 – 63, 𝑎𝑚𝑄12−3 20 – 55, 𝑎𝑚𝑄12−3 20 – 58, 𝑎𝑚𝑄12−3 – 42, 𝑎𝑚𝑄12−3 q (l/s.m) 0.14 0.26 1.65 0.24 1.13 0.47 0.95 3.93 3.88 7.14 7.65 5.94 K (m/d) 0.47 1.93 16.31 0.57 2.89 2.04 4.65 10.59 23.74 18.57 37.35 25.81 Tên lỗ khoan G7 G9 G10 G11 LK602 LK603 LK432 LK408 LK906 LK905 LK907 LK908 LK424 G1 G2 G4 G6 23 – 86, 𝑎𝑚𝑄12−3 39 – 64, 𝑎𝑚𝑄12−3 41 – 59, 𝑎𝑚𝑄12−3 43 – 65, 𝑎𝑚𝑄12−3 55 – 72, 𝑎𝑚𝑄12−3 3.81 0.98 1.82 2.72 2.63 11.05 5.04 15.72 18.84 31.0 LK403 LK406 LK409 LK601 LK606 VB01 Chiều sâu khai thác, tuổi địa chất 39 – 63, 𝑎𝑚𝑄12−3 44 – 69, 𝑎𝑚𝑄12−3 43 – 68, 𝑎𝑚𝑄12−3 44 – 71, 𝑎𝑚𝑄12−3 26 – 72, 𝑁13 – 28, 𝑁13 61 – 115, 𝑁13 18 – 53, 𝐷2−3 19 – 27, 𝐷2−3 20 – 56, 𝐷2−3 19 – 43, 𝐷2−3 – 40, 𝑎𝑑𝑄11−2 & 𝐷2−3 , 𝑃2 11 – 58, 𝑂3 − 𝑆1 – 21, 𝑂3 − 𝑆1 10 – 58, 𝑂3 − 𝑆1 19 – 45, 𝑂3 − 𝑆1 16 – 38, 𝑂3 − 𝑆1 43-69, 𝑎𝑚𝑄12−3 q (l/s.m) 1.92 1.96 1.99 3.38 4.23 3.68 2.38 1.79 0.09 3.16 3.3 2.28 K (m/d) 15.6 24.76 34.72 15.6 6.23 37.69 8.06 5.07 0.26 28.52 15.74 1.25 0.88 0.78 0.36 0.09 0.58 1,98 2.4 5.95 0.72 0.98 2.55 17,7 Sau nhập thông số ĐCTV tiến hành phân vùng hệ số thấm hệ số nhả nước (xem hình 14 15) 35 Hình 14 Phân vùng hệ số thấm tầng chứa nước qh mô hình Hình 15 Phân vùng hệ số thấm tầng chứa nước qp mơ hình d Dữ liệu giá trị bổ cập bốc Bản đồ liệu giá trị bổ cập xác định sở tài liệu lượng mưa trạm khí tượng tồn vùng nghiên cứu (Cồn Cỏ, Cửa Việt, Đơng Hà, Khe Sanh, Gia Vịng), giá trị tính tốn mơ hình bổ cập Wetspa phần Bản đồ liệu bốc lấy tương tự trên, giá trị lượng bốc ngầm giới hạn chiều sâu 5m tính từ bề mặt đất Hình 16 Bản đồ lượng bổ cập nước đất tính trung bình theo ngày tính cho mùa khơ tháng I – VII năm 1981 (hình phải) diện tích phân bổ lượng bổ cập đến tầng chứa nước (hình phải) 36 e Các biên điều kiện biên mơ hình Dựa vào điều kiện thực tế vùng nghiên cứu, loại biên mơ hình bao gồm: Biên sông (River), biên mực nước không đổi (H = const), biên khơng dịng chảy (Q=0), biên tổng hợp (GHB) Trị số mực nước biên sông xác định theo tài liệu quan trắc điểm quan trắc nước mặt mặt cắt + Biên khơng dịng chảy đặt phía tiếp xúc đá gốc với tầng chứa nước; + Biên mực nước khơng đổi phía Tây Bắc – Đơng Nam mơ miền cấp từ bên vào vùng nghiên cứu + Biên mực nước không đổi lớp theo đường bờ mô tiếp xúc trực tiếp với nước biển Giá trị mực nước biên lấy theo giá trị trung bình quan trắc mực nước biển thực tế Trạm Khí tượng Hải văn Cồn Cỏ + Biên GHB (tổng hợp) bố trí theo đường bờ với lớp 2, mô tiếp xúc không trực tiếp với biển thông qua giá trị cản thấm Conductance Giá trị mực nước biên lấy theo giá trị trung bình quan trắc mực nước biển thực tế Trạm Khí tượng Hải văn Cồn Cỏ + Biên sơng bề mặt lớp mô tiếp xúc không trực tiếp sơng lớp mơ hình thơng qua giá trị cản thấm Conductance Trị số sức cản thấm C sông xác định theo tài liệu thí nghiệm seepage năm 2017 lớp bùn đáy sơng Hình 17 Các loại biên lớp – TCN qh Hình 18 Các loại biên lớp – TCN qp e Lưới sai phân bước thời gian chỉnh lý 37 Mơ hình dự báo xây dựng với lưới sai phân hữu hạn kích thước bước 250x250m phần trung tâm vùng nghiên cứu đan dày với lưới nhỏ có kích thước 50x50m Mơ hình có tổng số 133 hàng 134 cột (xem hình 19) Mơ hình dự báo sử dụng giá trị đầu vào tính tốn trước ảnh hưởng BĐKH&NBD Các giá trị mực nước biên dòng chảy, lượng bổ cập tính từ lượng mưa… Phân bố biên mặn thời điểm tầng chứa nước sử dụng làm giá trị ban đầu để mơ hình tính tốn q trình dịch chuyển Mơ hình dự báo chạy đến mốc năm 2100 với kịch phát thải (B1, B2, A2) Vị trí biên biển tác động nước biển dâng điều chỉnh sau khoảng thời gian 20 năm Hình 19 Lưới sai phân mơ hình 2.4 Kết mơ hình dự báo dịch chuyển biên mặn nước đất ứng với kịch biến đổi khí hậu Mơ hình trạng xây dựng nhằm giả lập lại biến đổi tài nguyên nước khứ Mô hình trải qua bước chỉnh lý Bước chỉnh lý ổn định nhằm chỉnh lý điều kiện biên thông số ĐCTV TCN mô hình Bước chỉnh lý khơng ổn định mơ hình thực nhằm chỉnh lý điều kiện khai thác bổ cập tầng chứa nước 38 Bước chỉnh lý không ổn định chạy: thời đoạn, chia 36 bước thời gian tính tốn Mỗi bước thời gian dài tương đương tháng Thực chất toán ngược thực nhiều lần nhằm giảm thiểu tối đa sai số đưa nghiệm toán tiến gần với giá trị thực tế quan sát vị trí quan trắc mực nước Sau hoàn tất việc nhập tất liệu đầu vào trình bày mục trên, việc hiệu chỉnh mơ hình thực theo hai bước: - Bước thứ giải toán ổn định (Steady State) với mục đích kiểm tra sơ lại thơng số ĐCTV điều kiện biên mơ hình (chủ yếu liệu thuộc tính) Bài tốn kết thúc mực nước ban đầu mơ hình thời điểm tháng 1/2015 xác lập với sai số so với thực tế nằm giới hạn cho phép - Bước thứ hai giải tốn khơng ổn định (Transient) với mục đích làm xác hóa thơng số ĐCTV, điều kiện biên biến đổi theo thời gian trị số hệ số nhả nước lớp Bài toán chỉnh lý kết thúc xác lập động thái mực nước theo thời gian với sai số so với thực tế bước tính toán đạt giá trị cho phép Thời gian chạy chỉnh lý không ổn định từ 1/1/2015 đến 31/12/2017 Kết hiệu chỉnh toán ổn định Hiệu chỉnh toán ổn định tiến hành qua giai đoạn: Giai đoạn 1: Trước thực toán hiệu chỉnh cần thiết tiến hành chỉnh sửa lỗi nhập liệu chủ quan lẫn khách quan Giai đoạn 2: Tiếp theo giải toán với điều kiện mực nước ban đầu tầng mực nước biên xác lập theo tài liệu liên quan từ đồ thủy đẳng áp thủy đẳng cao Lúc cần tiến hành kiểm tra thông số ĐCTV điều kiện biên, đặc biệt kiểm tra tính tương thích hệ thống liệu kèm theo mơ hình dịng chảy nước đất Việc giải toán thực nhiều lần với nghiệm lần dùng làm điều kiện mực nước ban đầu cho lần sau Việc thực toán lặp chấm dứt mực nước ban đầu cho mơ hình dịng chảy nước đất xác lập Kết hiệu chỉnh toán ổn định thời điểm 1/2015 xác lập trạng mực nước với loại sai số ghi nhận vị trí quan trắc 39 Tóm lại, tốn ổn định hoàn thành đạt yêu cầu cơng tác hiệu chỉnh mơ hình khu vực lớn Kết mực nước tầng chứa nước xác lập Quá trình thực hiệu chỉnh cho thấy điều kiện biên ranh giới mô hình khơng ảnh hưởng nghiệm tốn, đặc biệt khu vực trung tâm, nơi có cụm cơng trình khai thác Nhà máy nước Gio Linh Điều cho thấy việc sơ đồ hóa xác lập điều kiện biên phù hợp thực tế Sau xây dựng mơ hình dịng chảy nước đất, tiến hành dự báo dịch chuyển biên mặn tầng chứa nước thành tạo Đệ tứ vùng nghiên cứu ảnh hưởng BĐKH&NBD (ứng với kịch khác xây dựng phần trước); kết cho thấy xu mặn nhạt biến đổi phức tạp, cụ thể: a Đối với tầng chứa nước Holocene Kết mơ hình dịng chảy theo trạng năm 2015, 2018 Hình 20 Mực nước TCN qh năm 2015 Hình 21 Mực nước TCN qh năm 2018 40 Hình 22 Mực nước TCN qp năm 2015 Hình 23 Mực nước TCN qp năm 2018 Kết mơ hình trạng xâm nhập mặn TCN qh năm 2015, 2018 Hình 24 Xâm nhập mặn TCN qh năm 2015 Hình 25 Xâm nhập mặn TCN qh năm 2018 - Kết mơ hình dự báo cho thấy biến đổi nước đất tầng Holocene rõ rệt: Diện tích nước mặn biến đổi phức tạp qua giai đoạn, khu vực nhiễm mặn lục đại có 41 xu hướng giảm, xong vùng ven biển khu vực chịu ảnh hưởng sông sơng Bến Hải Thạch Hãn diện tích mặn lại tăng lên tác động ngập (xem hình 3.8) Những thập kỷ cuối kỷ 21, ranh giới mặn nhạt bị tác động lớn hơn, biến đổi nhanh tầng chịu ảnh hưởng mạnh mực nước biển dâng cao, phía trơng lục địa q trình rửa mặn diễn nhanh chóng; theo kịch phát thải cao (A2), đến năm 2100 diện tích nước mặn toàn vùng nghiên cứu 17,582 km2 Năm 2020 Năm 2040 Năm 2080 Năm 2100 42 *Ghi chú: màu đỏ chỉ phần nước đất có M>1g/l, màu xanh chỉ phần nước đất có M