Groundwater Modeling: Protocols and Numerical Methods - G. Oude Essink, Utrecht University

Mô hình hóa nước ngầm là phương pháp số hóa quá trình vận động, phân bố và chất lượng nước dưới lòng đất. Bài viết phân tích các phương pháp, phần mềm và ứng

• Tác giả: Gualbert Oude Essink
• Trường học: Utrecht University
• Chuyên ngành: Thủy văn
• Đề tài: Groundwater Modelling
• Loại tài liệu: Bài giảng
• Năm xuất bản: 2000
• Địa điểm: Utrecht

Tóm tắt

I. Tổng quan về mô hình hóa nước dưới đất

Mô hình hóa nước dưới đất (groundwater modeling) là quá trình sử dụng các mô hình toán học để mô phỏng và dự báo hành vi của nước trong các tầng chứa nước. Sự phát triển của lĩnh vực này gắn liền với bước đột phá trong công nghệ tính toán. Máy tính kỹ thuật số có khả năng lưu trữ, xử lý dữ liệu và thực hiện các phép tính phức tạp vượt xa năng lực vật lý của con người. Mô hình hóa giúp tích hợp lý thuyết hiện có với dữ liệu đo đạc thực tế. Quá trình xây dựng mô hình bao gồm nhiều bước: xác định mục đích, khái niệm hóa mô hình toán học, lựa chọn mã tính toán, chọn tham số và ứng suất thủy văn. Mỗi bước đòi hỏi sự hiểu biết sâu về hệ thống thủy văn nghiên cứu. Các mô hình được phân loại dựa trên thiết kế (một chiều, hai chiều, ba chiều), dựa trên quy trình trong chu trình thủy văn (dòng chảy, vận chuyển chất tan, vận chuyển nhiệt), và dựa trên ứng dụng thực tế. Sự kết hợp giữa lý thuyết vững chắc và dữ liệu chất lượng cao là yếu tố then chốt để xây dựng mô hình đáng tin cậy.

1.1. Lịch sử phát triển mô hình hóa thủy văn

Lịch sử mô hình hóa thủy văn bắt đầu từ nhu cầu tích hợp lý thuyết với dữ liệu đo đạc. Giai đoạn đầu, các phương pháp giải tích được sử dụng để giải các bài toán đơn giản. Sự ra đời của máy tính đánh dấu bước ngoặt quan trọng. Các nhà thủy văn có thể giải quyết bài toán phức tạp hơn với số lượng dữ liệu lớn. Từ thập niên 1960, các mã mô hình số đầu tiên xuất hiện. Phát triển phần cứng và phần mềm liên tục mở rộng khả năng ứng dụng. Ngày nay, mô hình hóa ba chiều và mô hình phân tán trở thành công cụ tiêu chuẩn trong quản lý tài nguyên nước.

1.2. Phân loại các mô hình toán học thủy văn

Mô hình toán học trong thủy văn được phân loại theo nhiều tiêu chí. Dựa trên thiết kế, mô hình chia thành một chiều, hai chiều và ba chiều. Mỗi loại phù hợp với bài toán cụ thể và mức độ phức tạp khác nhau. Dựa trên quy trình thủy văn, mô hình bao gồm mô hình dòng chảy, mô hình vận chuyển chất tan và mô hình vận chuyển nhiệt. Dựa trên ứng dụng, mô hình phục vụ quản lý tài nguyên nước, đánh giá ô nhiễm hoặc dự báo lũ lụt. Việc lựa chọn loại mô hình phù hợp là bước quan trọng đầu tiên trong quy trình mô hình hóa.

II. Thách thức trong mô hình hóa nước dưới đất

Mô hình hóa nước dưới đất đối mặt với nhiều thách thức phức tạp. Bài toán thu thập dữ liệu là rào cản lớn nhất. Thông tin về tham số nền đất như độ dẫn thủy lực, độ thấm, độ rỗng và hệ số trữ riêng thường không đầy đủ. Dữ liệu đầu vào và đầu ra bao gồm mưa, bốc hơi, dòng chảy bề mặt, thấm và dòng cơ sở cần được đo đạc chính xác. Điều kiện biên đặt ra vấn đề khác. Phần dưới mô hình thường dùng biên không thấm. Tuy nhiên, biên tự nhiên không phải lúc nào cũng có sẵn. Khi đó, biên nhân tạo phải được mô phỏng. Điều này có thể gây sai số đáng kể. Dữ liệu hóa địa chất cũng cần thiết: cation, anion, nhiệt độ, pH và đồng vị. Sai số trong dữ liệu đầu vào sẽ tích lũy qua quá trình tính toán. Hiệu chuẩn mô hình đòi hỏi dữ liệu quan trắc đủ dài và liên tục để đảm bảo độ tin cậy.

2.1. Bài toán thu thập và xử lý dữ liệu

Thu thập dữ liệu chất lượng cao là thách thức hàng đầu trong mô hình hóa nước dưới đất. Tham số nền đất phải được xác định qua quan trắc hiện trường và thí nghiệm. Độ dẫn thủy lực và độ thấm thay đổi không gian phức tạp. Dữ liệu thủy văn như lượng mưa, bốc hơi và thấm cần chuỗi thời gian dài. Hóa địa chất yêu cầu phân tích cation, anion, đồng vị và hợp chất hữu cơ. Dữ liệu không đầy đủ dẫn đến mô hình không đáng tin cậy. Phương pháp nội suy và ngoại suy được áp dụng để bổ sung dữ liệu thiếu.

2.2. Xác định điều kiện biên và ranh giới mô hình

Điều kiện biên quyết định độ chính xác của mô hình nước dưới đất. Ba loại biên chính được sử dụng: Dirichlet, Neumann và Cauchy. Biên tự nhiên như sông, núi hoặc tầng không thấm nên được ưu tiên sử dụng. Khi biên tự nhiên không có sẵn, biên nhân tạo phải được thiết lập cẩn thận. Ranh giới mô hình xác định kích thước và vùng quan tâm. Cần cân nhắc giữa diện tích mô phỏng và chi phí tính toán. Thiết lập biên sai lệch gây ra kết quả mô hình không phản ánh đúng thực tế hệ thống thủy địa chất.

III. Phương pháp mô hình hóa nước dưới đất hiệu quả

Phương pháp mô hình hóa nước dưới đất dựa trên các phương trình toán học mô tả quy trình thủy địa chất. Phương trình dòng chảy sử dụng định luật Darcy kết hợp với phương trình liên tục. Độ dẫn thủy lực và độ thấm là tham số vật lý cốt lõi. Phương trình vận chuyển chất tan mô tả quá trình đối lưu và phân tán. Phương trình vận chuyển nhiệt xem xét dẫn nhiệt và đối lưu nhiệt. Về mặt số, nhiều phương pháp được áp dụng. Phương pháp sai phân hữu hạn chia lưới đều đặn và giải phương trình tại các nút lưới. Phương pháp phần tử hữu hạn linh hoạt hơn với lưới không đều. Phương pháp giải tích phần tử kết hợp ưu điểm của cả hai cách tiếp cận. Phương pháp đặc trưng và phương pháp đi ngẫu nhiên phù hợp với bài toán vận chuyển chất tan. Mỗi phương pháp có ưu nhược điểm riêng. Lựa chọn phụ thuộc vào tính chất bài toán và tài nguyên tính toán hiện có.

3.1. Mô tả toán học quy trình thủy địa chất

Mô tả toán học quy trình thủy địa chất bắt đầu từ định luật Darcy. Định luật này liên hệ vận tốc dòng chảy với gradient thủy lực và độ dẫn thủy lực. Phương trình liên tục bảo toàn khối lượng chất lỏng trong hệ thống. Kết hợp hai phương trình tạo thành phương trình dòng chảy nước dưới đất. Phương trình vận chuyển chất tan (advection-dispersion) mô tả sự di chuyển của chất ô nhiễm trong tầng chứa nước. Phương trình vận chuyển nhiệt xem xét cả dẫn nhiệt và đối lưu. Các phương trình phi tuyến thường yêu cầu giải pháp số để xử lý.

3.2. Các phương pháp số phổ biến trong mô hình hóa

Phương pháp số là công cụ cốt lõi để giải các phương trình mô hình nước dưới đất. Phương pháp sai phân hữu hạn (FDM) xấp xỉ đạo hàm bằng hiệu phân tại các điểm lưới đều. Phương pháp phần tử hữu hạn (FEM) sử dụng lưới linh hoạt, phù hợp với hình học phức tạp. Phương pháp giải tích phần tử (AEM) kết hợp giải pháp giải tích và kỹ thuật số. Phương pháp đặc trưng (MOC) xử lý tốt bài toán vận chuyển chất tan. Phương pháp đi ngẫu nhiên mô tả phân tán bằng cách di chuyển các hạt. Mỗi phương pháp đòi hỏi cân nhắc về tính ổn định và sai số.

IV. Ứng dụng và triển vọng mô hình hóa nước dưới đất

Mô hình hóa nước dưới đất có ứng dụng rộng rãi trong quản lý tài nguyên nước và bảo vệ môi trường. Các mã mô hình phổ biến như MODFLOW, MOC và SEAWAT được sử dụng toàn cầu. MODFLOW mô phỏng dòng chảy nước dưới đất ba chiều với nhiều gói tính toán. MOC xử lý bài toán vận chuyển chất tan trong điều kiện mật độ khác nhau. SEAWAT kết hợp dòng chảy và vận chuyển trong môi trường nước mặn. Yêu cầu tính toán phụ thuộc vào tốc độ xử lý, kích thước mô hình và hiệu quả trình biên dịch. Hệ thống máy trạm UNIX với RAM lớn mở ra khả năng mô hình hóa ba chiều chi tiết. Các gói tính toán (packages) xử lý nguồn bên ngoài và điều kiện biên phức tạp. Tương lai của lĩnh vực hướng đến mô hình phân tán, tích hợp dữ liệu viễn thám và ứng dụng học máy. Mô hình hóa nước dưới đất tiếp tục là công cụ không thể thiếu trong quản lý bền vững tài nguyên nước ngầm toàn cầu.

4.1. Các mã mô hình nước dưới đất phổ biến

Nhiều mã mô hình nước dưới đất được phát triển và sử dụng rộng rãi trên toàn thế giới. MODFLOW là mã tiêu chuẩn mô phỏng dòng chảy ba chiều do USGS phát triển. MOC xử lý vận chuyển chất tan với khả năng thích ứng mật độ khác nhau. Các gói (packages) trong mã mô hình xử lý nguồn bên ngoài vào khối lưới hiệu quả. Gói SIP (Strongly Implicit Procedure) và SSOR (Slice-Successive Overrelaxation) giải hệ phương trình tuyến tính. Lựa chọn mã phù hợp phụ thuộc vào bài toán cụ thể và dữ liệu sẵn có.

4.2. Yêu cầu tính toán và triển vọng tương lai

Yêu cầu tính toán của mô hình nước dưới đất phụ thuộc vào nhiều yếu tố kỹ thuật. Tốc độ xử lý và dung lượng RAM ảnh hưởng trực tiếp đến thời gian thực thi. Kích thước mô hình quyết định dung lượng bộ nhớ cần thiết. Hệ thống máy tính hiện đại với RAM lớn cho phép mô hình hóa ba chiều chi tiết. Triển vọng tương lai bao gồm mô hình phân tán trên nền tảng đám mây. Tích hợp dữ liệu viễn thám và ứng dụng trí tuệ nhân tạo mở rộng khả năng dự báo. Công nghệ mới giúp giảm thời gian tính toán và tăng độ chính xác mô hình.

Nội dung chính

Tổng quan về giáo trình

Tài liệu bài giảng Groundwater Modelling I (Mã số: L4018/GWM I, phát hành tháng 9 năm 2000) do Tiến sĩ Gualbert Oude Essink biên soạn, trực thuộc Trung tâm Thủy văn Liên khoa Utrecht (Interfacultair Centrum voor Hydrologie Utrecht - ICHU), Viện Khoa học Trái đất và Khoa Địa vật lý thuộc Đại học Utrecht (Utrecht University), Hà Lan. Tập bài giảng này được phát triển dựa trên một phần tài liệu Hydrological models (f15D) của Đại học Công nghệ Delft (Delft University of Technology), phục vụ trực tiếp cho chương trình đào tạo chuyên ngành Thủy văn (study-path Hydrologie).

Mục tiêu đào tạo của giáo trình là trang bị cho người học hiểu biết định lượng và định tính về các quá trình động lực học nước dưới đất thông qua công cụ mô hình hóa số trị (numerical modelling). Giáo trình định hướng người học thiết lập và vận hành một quy trình mô hình hóa hoàn chỉnh: từ việc phân tích bài toán thực tế ngoài hiện trường, xây dựng mô hình khái niệm, lựa chọn mã nguồn tính toán, thiết kế lưới không gian và bước thời gian, thiết lập điều kiện biên và điều kiện ban đầu, đến quá trình hiệu chỉnh (calibration), kiểm định và phân tích độ tin cậy của kết quả mô phỏng.

Về mặt cấu trúc, tài liệu được phân chia thành hai phần chính:

• Phần I - Giao thức mô hình hóa (Modelling Protocol): Trình bày trình tự phương pháp luận từng bước để giải quyết một bài toán thủy văn bằng mô hình hóa số trị, đồng thời nhận diện các hạn chế và sai số kỹ thuật trong quá trình thực hiện.
• Phần II - Mô hình hóa nước dưới đất (Groundwater Modelling): Trình bày hệ thống các phương trình vi phân đạo hàm riêng mô tả quá trình dòng chảy chất lưu, vận chuyển chất tan và truyền nhiệt; các phương pháp giải số trị tương ứng và cơ chế hoạt động của các phần mềm mô phỏng chuyên dụng.

Điểm đặc thù của giáo trình là sự kết hợp chặt chẽ giữa cơ sở giải tích toán lý và quy trình kỹ thuật thực hành trên máy tính (sử dụng các mã nguồn như MODFLOW, MOC, MOC3D), gắn liền với dữ liệu thực nghiệm về các tầng chứa nước tại khu vực trũng ven biển của Hà Lan.

---

Nội dung kiến thức cốt lõi

Các chương và chủ đề chính

Nội dung giáo trình được cấu trúc qua 9 chương chuyên đề và hệ thống phụ lục kỹ thuật:

• Chương 1 - Giới thiệu (Introduction): Tổng thuật lịch sử phát triển ngành thủy văn từ thời cổ đại (Anaxagoras, Theophrastus), thời kỳ Phục hưng (Leonardo da Vinci) đến các phát minh cơ bản thế kỷ XIX-XX (Darcy, Manning, Sherman, Gumbel). Chương này nêu rõ vai trò của máy tính số, định nghĩa các khái niệm cốt lõi (hệ thống thủy văn, biến số, thông số) và phân định phạm vi giữa Địa lý tự nhiên, Địa hóa học, Địa vật lý, Địa chất học, Thủy văn học và Địa chất thủy văn.
• Chương 2 - Phân loại mô hình toán học (Classification of mathematical models): Phân chia mô hình theo thiết kế cấu trúc (giải tích đối chứng số trị; tất định đối chứng ngẫu nhiên; tập trung đối chứng phân bố; hộp đen, hộp xám đối chứng hộp trắng; thực nghiệm, khái niệm đối chứng dựa trên quy luật vật lý; dừng đối chứng không dừng; tuyến tính đối chứng phi tuyến; số chiều không gian từ 0D đến 3D); phân loại theo chu trình thủy văn và phân loại theo mục đích ứng dụng (nghiên cứu khoa học, thiết kế kỹ thuật, đánh giá/diễn giải, dự báo và quản lý tài nguyên nước).
• Chương 3 - Phương pháp luận mô hình hóa (Methodology of modelling): Khảo sát chi tiết chu trình mô hình hóa theo khung của Anderson & Woessner (1992), bao gồm: xác định mục tiêu, sơ đồ hóa và khái niệm hóa bài toán, lựa chọn mã nguồn máy tính, thiết kế lưới (grid design) và bước thời gian (temporal discretisation), cài đặt điều kiện biên (Dirichlet, Neumann, Cauchy), thiết lập điều kiện ban đầu, hiệu chỉnh thông số (phương pháp thử sai và phương pháp nghịch đảo tự động), đánh giá sai số (ME, MAE, RMS) và tiêu chí hội tụ lặp.
• Chương 4 - Thu thập dữ liệu (Data gathering): Hướng dẫn thu thập và xử lý các thông số hình học tầng chứa nước, hệ số thấm, độ truyền dẫn, hệ số tích nước, tải trọng thủy văn và dữ liệu địa hóa.
• Chương 6 - Mô tả toán học các quá trình thủy địa chất (Mathematical description of hydrogeologic processes): Thiết lập phương trình chuyển động của chất lưu (định luật Darcy), phương trình liên tục, phương trình dòng chảy tổng quát trong môi trường lỗ rỗng, phương trình trạng thái phụ thuộc mật độ; phương trình bình lưu - phân tán (advection-dispersion equation) mô tả vận chuyển chất tan và phương trình dẫn nhiệt - đối lưu nhiệt (conduction-convection equation).
• Chương 7 - Kỹ thuật giải số phương trình vi phân (Numerical solution techniques): Trình bày kỹ thuật khai triển chuỗi Taylor, phương pháp khử Gauss-Jordan, phương pháp sai phân hữu hạn (Finite Difference Method - FDM), phương pháp phần tử hữu hạn (Finite Element Method - FEM), phương pháp phần tử giải tích (Analytic Element Method - AEM), phương pháp đặc trưng (Method of Characteristics - MOC) và phương pháp bước đi ngẫu nhiên (Random Walk Method).
• Chương 8 & 9 - Khía cạnh số trị và các mã nguồn tính toán (Numerical aspects and selected groundwater codes): Phân tích tính ổn định số trị, sai số cắt cụt (truncation errors), sai số dao động (oscillation errors) và tiêu chí số Peclet; cấu trúc và nguyên lý vận hành của mã USGS MODFLOW (với các gói giải SIP, SSOR) cùng mã MOC/MOC3D thích ứng cho bài toán dòng chảy có mật độ thay đổi.

Logic phát triển nội dung đi từ việc thiết lập nền tảng nhận thức và phương pháp luận tổng quát (Phần I), tiếp nối bằng việc giải thức toán học các quá trình vật lý, rời rạc hóa số trị và ứng dụng các mã phần mềm máy tính tương ứng (Phần II).

Kiến thức nền tảng được xây dựng

• Các lý thuyết nền tảng: Định luật bảo toàn khối lượng, bảo toàn động lượng và bảo toàn năng lượng trong môi trường xốp; khái niệm Thể tích đại diện hữu hạn (Representative Elementary Volume - REV, theo Bear 1972) đóng vai trò chuyển tiếp từ cấp độ lỗ rỗng vi mô sang cấp độ môi trường liên tục; nguyên lý chồng chất (superposition) trong các hệ vi phân tuyến tính.
• Các nguyên lý cốt lõi: Định luật Darcy (1856) cho dòng chảy tầng trong môi trường rỗng; cơ chế dòng chảy phụ thuộc mật độ giữa nước ngọt và nước mặn; động học vận chuyển chất tan bao gồm quá trình bình lưu (advection), khuếch tán phân tử (molecular diffusion) và phân tán thủy lực (mechanical dispersion).
• Các khung lý thuyết chuẩn: Khung giao thức mô hình hóa chuẩn tắc phỏng theo Anderson & Woessner (1992); khung thống kê bình phương tối thiểu có trọng số (weighted least squares) và khung thống kê Fisherian dùng trong bài toán tối ưu hóa thông số nghịch đảo (inverse problem).

Kỹ năng phát triển

• Kỹ năng kỹ thuật: Thiết kế mạng lưới sai phân khối trung tâm (block-centered) và mắt lưới trung tâm (mesh-centered); thiết lập cấu trúc tệp dữ liệu đầu vào (input data files) cho phần mềm MODFLOW và MOC; cấu hình bộ nhớ mở rộng (Extended Memory RAM) thông qua các trình biên dịch FORTRAN (như Lahey Fortran F77L-EM/32) phục vụ mô phỏng 2D và 3D.
• Kỹ năng phân tích: Tính toán và đánh giá các chỉ số sai số định lượng trong hiệu chỉnh mô hình: Sai số trung bình (Mean Error - ME), Sai số tuyệt đối trung bình (Mean Absolute Error - MAE), Sai số toàn phương trung bình (Root Mean Squared Error - RMS); phân tích tiêu chí ổn định của bước thời gian và kích thước lưới $\Delta x$ qua số Peclet ($Pe \le 2 - 10$); kiểm tra cân bằng nước tổng thể (đạt ngưỡng sai số chấp nhận dưới 0.1% đến 1%).
• Kỹ năng thực hành: Đơn giản hóa cấu trúc địa tầng tự nhiên thành mô hình khái niệm nhiều tầng (như cấu trúc mặt cắt Holland profile gồm tầng chứa cát Pleistocene và lớp thấm yếu sét Holocene); gán đúng bản chất vật lý cho các điều kiện biên thủy văn tại biên biển, sông ngòi, đường phân chia dòng ngầm và hệ thống giếng khai thác.

---

Phương pháp giảng dạy và học tập

Giáo trình áp dụng phương pháp tiếp cận diễn dịch kết hợp thực nghiệm số trị: lý thuyết toán lý được trình bày gắn liền với phân tích bản chất cơ học chất lưu, tiếp nối bằng việc chuyển đổi các phương trình vi phân thành hệ đại số rời rạc và kiểm chứng thông qua các bài thực hành trên máy tính (computer practicals). Tài liệu nhấn mạnh quan điểm phản biện khoa học: mô hình chỉ là công cụ hỗ trợ phân tích, không đại diện hoàn toàn cho thực tế khách quan, và việc sử dụng các phần mềm đồ họa hiện đại không thể thay thế cho việc đánh giá tính đúng đắn của các giả định thủy địa chất.

Hệ thống bài tập và nghiên cứu tình huống (case studies) tập trung giải quyết các bài toán thủy văn có thực:

• Mô hình hóa hệ thống dòng chảy nước dưới đất tại vùng cồn cát ven biển của Công ty Cấp nước Đô thị Amsterdam (Gemeentewaterleidingen Amsterdam). Tình huống này phân tích sự biến đổi của dòng nước ngầm và hiện tượng nhiễm mặn trước và sau sự kiện cải tạo hồ Haarlemmermeer (năm 1853-1854), đối chiếu với các kịch bản khai thác giếng ngầm ở các năm 1956 và 1987.
• Mô phỏng tác động thủy triều lên hệ thống tầng chứa nước ven bờ nhằm chỉ rõ sai số khi áp dụng mô hình trạng thái dừng (steady-state) cho các quá trình biến đổi nhanh (transient).

Phương pháp đánh giá kết quả học tập dựa trên:

1. Độ chính xác trong việc lựa chọn điều kiện biên và rời rạc hóa lưới mô phỏng.
2. Khả năng lập trình và hiệu chỉnh mô hình đạt các mức chỉ tiêu (levels of calibration $\xi$).
3. Khả năng kiểm soát sai số cân bằng nước và giải trình nguyên nhân khi thuật toán lặp không hội tụ (do sai sót thiết kế lưới, tương phản độ dẫn thủy lực quá mức hoặc lỗi bộ nhớ máy tính).

Để tự học hiệu quả, người học được hướng dẫn theo dõi trình tự từng bước trong sơ đồ giao thức (Modelling Protocol flowchart), sử dụng bảng tổng hợp công thức (Formula sheet) ở trang 189 để kiểm tra phương trình, đồng thời đối chiếu các thuật ngữ chuyên ngành đã được định nghĩa chuẩn xác trong Chương 1 và Chương 2.

---

Điểm nổi bật và cập nhật

Tài liệu được cập nhật dựa trên bối cảnh khoa học và công nghệ tính toán tại thời điểm năm 2000:

• Kế thừa và mở rộng: Giáo trình phát triển trực tiếp từ các bài giảng Hydrological models (f15D) của Đại học Công nghệ Delft, bổ sung toàn diện phần lý thuyết toán học về dòng chảy mật độ không đồng nhất (non-uniform density distributions) và chuyển động của ranh giới mặn - ngọt.
• Tích hợp bước tiến công nghệ phần cứng: Cập nhật xu hướng chuyển dịch tính toán từ các hệ thống máy tính lớn (mainframe) sang máy tính cá nhân để bàn (Personal Computer) và máy trạm UNIX. Giáo trình dẫn chứng dữ liệu định lượng của Cutaia (1990) và đưa ra bảng đo kiểm hiệu năng thực tế (benchmark tests) của tác giả Oude Essink trên các bộ vi xử lý từ PC XT, Pentium Pro 200 MHz, Pentium III 733 MHz đến trạm làm việc UNIX Geophysics với bộ nhớ 1500 Mb EM RAM.

Hệ thống máy tính	Tốc độ tương đối (so với PC)	Cấu hình & Ghi chú
PC chuẩn	1.0	Nền tảng gốc DOS 640 Kb RAM
PC XT	1.2	Xử lý dữ liệu cơ bản
Pentium Pro 200 MHz	382.4	Đo kiểm thực nghiệm bởi Oude Essink
Pentium III 733 MHz	1461.3	Hỗ trợ mô phỏng không gian 3 chiều
Trạm làm việc UNIX (Geophysics)	Hiệu năng cao	Trang bị 495 Mb đến 1500 Mb EM RAM

• Xu hướng thuật toán mới: Giới thiệu các mã nguồn tự động hóa ước lượng thông số (Automated parameter estimation codes) dựa trên phương pháp giải bài toán nghịch đảo (inverse problem) và kỹ thuật mô phỏng ngẫu nhiên Monte Carlo để tạo lập chuỗi số liệu mưa nhân tạo.
• Gắn kết thực tiễn và nghiên cứu chính sách: Tài liệu kết nối trực tiếp với nghiên cứu Phân tích Chính sách Quản lý Nước của Hà Lan (PAWN-study - Policy Analysis for the Watermanagement of the Netherlands [Pulles, 1985]), đồng thời giải quyết các thách thức môi trường đương thời như ô nhiễm nước ngầm do mưa axit, nitrat hóa và các hợp chất hữu cơ khó tan trong nước (NAPLs - nonaqueous phase liquids).

---

Đối tượng sử dụng giáo trình

Giáo trình Groundwater Modelling I được thiết kế cho các nhóm đối tượng học thuật xác định:

• Sinh viên và Học viên cao học: Sinh viên các năm cuối bậc đại học và học viên cao học chuyên ngành Thủy văn học (Hydrology), Địa lý tự nhiên (Physical Geography), Địa chất học (Geology), Địa vật lý (Geophysics) và Địa hóa học (Geochemistry) thuộc Đại học Utrecht và các viện nghiên cứu liên kết trong khuôn khổ ICHU.
• Yêu cầu kiến thức tiên quyết (Prerequisites):
  • Đã nắm vững kiến thức cơ sở về địa chất thủy văn, bao gồm định luật Darcy và các bài toán dòng chảy nước dưới đất ở trạng thái dừng.
  • Có hiểu biết về toán giải tích, chuỗi Taylor, phép biến đổi Laplace và các phương pháp giải số cơ bản.
  • Sinh viên thuộc Khoa Địa lý Tự nhiên được khuyến nghị bắt buộc phải hoàn thành học phần Địa chất thủy văn (Groundwaterhydrology - HYDB); toàn thể sinh viên được khuyến nghị theo học trước học phần Các quá trình vận chuyển thủy văn (Hydrological Transport Processes - L3041/KHTP).
• Giảng viên và Nghiên cứu viên: Giảng viên sử dụng tài liệu làm đề cương chuẩn cho các khóa học lý thuyết kết hợp thực hành phòng máy tính; các nghiên cứu viên và kỹ sư địa kỹ thuật - môi trường có thể sử dụng làm tài liệu tham khảo để tra cứu phương trình toán lý, kỹ thuật gán điều kiện biên và tiêu chuẩn đánh giá sai số mô hình.

---

Câu hỏi thường gặp

1. Giáo trình này được biên soạn phục vụ cho đối tượng nào?

Giáo trình phục vụ sinh viên năm cuối đại học, học viên sau đại học và các nhà nghiên cứu thuộc các ngành Khoa học Trái đất, Địa vật lý, Địa lý tự nhiên và Thủy văn công trình cần sử dụng công cụ mô hình số để tính toán tài nguyên và chất lượng nước ngầm.

2. Cần chuẩn bị những kiến thức nền tảng nào trước khi học?

Người học cần hoàn thành các học phần toán giải tích, phương pháp số cơ bản, cơ học thủy khí môi trường xốp và địa chất thủy văn cơ sở (học phần HYDB và L3041/KHTP), nắm vững định luật Darcy và nguyên lý dòng chảy tầng chứa nước.

3. Giáo trình này có điểm khác biệt gì so với các tài liệu thủy văn đại cương?

Tài liệu tập trung chuyên biệt vào phương pháp luận mô hình hóa số trị (Modelling Protocol) 8 bước chặt chẽ, không dừng lại ở giải tích lý thuyết thuần túy mà hướng dẫn chi tiết các khía cạnh kỹ thuật máy tính: cấu trúc bộ nhớ FORTRAN, cơ chế giải lặp số trị, kiểm soát dao động số Peclet và cách xử lý dữ liệu thực tế trên các mã nguồn tiêu chuẩn như MODFLOW và MOC.

4. Làm thế nào để tự học và thực hành hiệu quả theo giáo trình?

Người học cần tuân thủ tuần tự các bước trong quy trình mô hình hóa ở Phần I; sau đó đối chiếu các phương trình toán lý ở Phần II khi xây dựng bài toán. Trong quá trình chạy thử mô hình, cần đặc biệt theo dõi tỷ lệ sai số cân bằng nước (yêu cầu dưới 1%) và đánh giá định lượng sai số qua ba đại lượng thống kê ME, MAE và RMS.

5. Tài liệu có cung cấp danh mục tài nguyên bổ trợ nào không?

Giáo trình cung cấp danh mục chi tiết các viện nghiên cứu và tổ chức phân phối mã nguồn quốc tế (như USGS, IGWMC, Scientific Software Group, WMO/HOMS, cơ sở dữ liệu STOWA/SAMWAT của Hà Lan), hệ thống thư mục tài liệu tham khảo học thuật (Anderson & Woessner 1992, Bear 1972, Dooge 1968, Konikow & Bredehoeft 1978/1992) và một bảng tổng hợp công thức toán học chuyên dụng ở trang 189.

---

Kết luận

Tập bài giảng Groundwater Modelling I (L4018/GWM I) của G. Oude Essink thiết lập hệ thống phương pháp luận quy chuẩn cho việc mô hình hóa các quá trình dòng chảy và lan truyền chất trong môi trường nước ngầm. Giá trị học thuật cốt lõi của tài liệu nằm ở sự liên kết nhất quán giữa phương trình toán lý cơ bản, cấu trúc giải thuật số trị và quy trình hiệu chỉnh dữ liệu thực địa.

Tiến trình học tập được định hình mạch lạc: bắt đầu từ việc hiểu rõ bản chất vật lý của hệ thống thủy văn, tiến hành khái niệm hóa bài toán, thiết kế mạng lưới tính toán, áp dụng các điều kiện biên và kết thúc bằng bước kiểm định sai số nghiêm ngặt trên phần mềm. Để mở rộng kiến thức, người học có thể khai thác các mã nguồn công cộng từ Cục Khảo sát Địa chất Hoa Kỳ (USGS), Trung tâm Mô hình hóa Nước dưới đất Quốc tế (IGWMC) và các công trình nghiên cứu chuyên sâu được liệt kê trong danh mục thư mục của giáo trình.

Trích đoạn nội dung tài liệu

L4018/GWM I Groundwater Modelling September 2000 G. Oude Essink Utrecht University Interfaculty Centre of Hydrology Utrecht Institute of Earth Sciences Department of Geophysics f. 15,- Contents General introduction v I Modelling Protocol 1 1 Introduction 3 1.1 Historical developments towards hydrologic modelling .3 Some drawbacks in modelling. 6 2 Classification of mathematical models 9 2.2 Based on the design of the mathematical model .3 Based on the processes in the hydrologic cycle .4 Based on the application of the model.

15 3 Methodology of modelling 19 3.1 Define purpose of the modelling effort .2 Conceptualisation of a mathematical model .3 Selection of the computer code .5 Preliminary selection of parameters and hydrologic stresses .1 Evaluating the calibration .3 First model execution .6 Model verification .8 Presentation of results .10 Why can things go wrong ?. 53 i ii 4 Data gathering 55 II Groundwater Modelling 59 5 Introduction 61 5.1 Classification based on the design of the model. 68 6 Mathematical description of hydrogeologic processes 69 6.1 Fluid flow: equation of motion and continuity .1 Equation of motion: Darcy’s law .3 Hydraulic conductivity and permeability .4 Density of groundwater .6 Equation of continuity .7 Groundwater flow equation .8 Equation of state .2 Solute transport: advection-dispersion equation .1 Equation of solute transport .3 Heat transport: conduction-convection equation .1 Taylor series development .3 Steady state methods .4 Non-steady state methods .4 Gauss-Jordan elimination .5 Finite difference method .6 Finite element method .7 Analytic element method .8 Method of characteristics .9 Random walk method. 116 iii 8 Numerical aspects of groundwater models 119 8.1 Stability analysis of the advection-dispersion equation .3 Analysis of truncation and oscillation errors.

128 9 Some selected groundwater codes 133 9.1 External sources into a block: packages .4 Strongly Implicit Procedure package (SIP) .5 Slice-Successive Overrelaxation package (SSOR) .5 MOC, (2D) adapted for density differences .1 Theoretical background of the groundwater flow equation .2 Theoretical background of the solute transport equation. 164 Continuı̈teitsvergelijking: niet stationair 177 References 181 Consulted literature. 187 Some distributors of computer codes. 187 Formula sheet 189 Index 191 iv General introduction The development of models has been the direct outcome of the need to integrate our existing theories with all physical and measured data.

The key factor in this development is the (still ongoing) breakthrough in computation technologies, because it is the digital computer which is capable to store, to manipulate data and to execute complex calculations beyond the physical ability of man, yet within his mental capacity. In order to avoid that you, as a hydrogeologist-in-spe, will get stranded in the fine art of modelling, these lecture notes1 are written to show you some ropes in the fantastic, tempting, and yet creepy world of groundwater modelling. The presence course, which comes under the ICHU2 , is called Groundwater Modelling I. The aim of this course is to gain more insight in the behaviour of groundwater processes, quantitatively as well as qualitatively, by means of numerical modelling.

These lecture notes are divided into two parts: I. Modelling protocol in this part, the procedure of modelling hydrologic processes is described. A hydro- geologist should advance this procedure, from coping with a hydrological problem towards solving the problem by means of (numerical) modelling, while skipping all kinds of traps on the track; and II. Groundwater modelling in this part, features of groundwater flow and solute transport are considered and numerical solution techniques of partial differential equations are discussed.

Moreover, some groundwater computer codes will be treated. During the computer practicals, attention will be paid to the modelling of standard problems by means of the codes MODFLOW and MOC(3D). Primarily knowledge should comprise basic knowledge on hydrogeologic processes such as the Darcy equation and stationary groundwater flow; as well as basic knowledge on discretisation techniques such as Taylor series and simple numerical solution techniques. For students of the Department of Physi- cal Geography, it is strongly recommended to have followed Groundwaterhydrology (HYDB).

For all students, the lectures Hydrological Transport Processes (L3041/KHTP) are also rec- ommended. Gualbert Oude Essink, September 2000 g.