Nghiên cứu sử dụng công cụ mô phỏng SiWaΡг0 DSS cho yêu cầu Việt Nam

Luận văn thạc sĩ phân tích investigation of a lysimenter using the simulation tool siwapro dss and adaptation of this program, đánh giá thực trạng, chỉ ra hạn chế, đề xuất giải

• Người hướng dẫn: Prof. Peter Wolfgang
• Trường học: Hanoi University of Sciences
• Chuyên ngành: Water Management
• Đề tài: Nghiên Cứu Sử Dụng Công Cụ Mô Phỏng SiWaΡг0 DSS Cho Yêu Cầu Việt Nam
• Loại tài liệu: Thesis
• Năm xuất bản: 2008
• Địa điểm: Hanoi

Mục lục chi tiết

• ACKNOWLEDGEMENTS
• 1. CHƯƠNG 1: INTRODUCTION
• 2. CHƯƠNG 2: FUNDAMENTALS OF SOIL HYDROLOGY
• 2.1. Definition of soil and unsaturated zone
• 2.2. Soil hydraulic parameters
• 2.3. Soil water balance
• 2.4. Soil water flow
• 3. CHƯƠNG 3: MATERIAL AND METHODS
• 3.1. Theoretical approaches and methodology
• 3.2. Finite element method
• 3.2.1. General information about lysimeter
• 3.2.2. Juelich lysimeter station description
• 3.2.3. Description of the Juelich lysimeter number 302
• 3.2.4. Water flow model
• 3.2.5. Description of the finite element mesh of the lysimeter
• 3.2.6. Description of software SiWaΡг0 DSS
• 3.2.6.1. Graphical user interface (GUI) and help system
• 3.3. Database Layer
• 3.4. Pedotransfer Functions
• 3.5. Import and Export Interfaces
• 3.6. Manual SiWaΡг0 DSS Mesh Generator
• 3.6.1. Create a simple 2D mesh
• 3.6.2. Definition internal curves
• 3.6.3. Inserting a background image as construction basis
• 3.6.4. Boundary condition editor
• 3.6.7. Data sets for calibrating the model
• 3.6.5. Soil hydraulic parameters
• 3.6.2. Extension and adaptation to Vietnam requirements
• 4. CHƯƠNG 4: DISCUSSION AND CONCLUSIONS
• APPENDIX
• appendix 1. Precipitation using for simulation
• appendix 2. Brief of output of simulation for 784 days
• appendix 3. Data from measurement
• appendix 4. Data from simulation for the days equivalent with measurement days

Tóm tắt

I. Tổng Quan Về Nghiên Cứu SiWaPro DSS Cho Việt Nam 55 ký tự

Nghiên cứu này tập trung vào việc sử dụng SiWaPro DSS, một công cụ mô phỏng mạnh mẽ, để mô hình hóa và mô phỏng quá trình dòng chảy nước trong vùng chưa bão hòa. Vùng chưa bão hòa đóng vai trò quan trọng trong nhiều khía cạnh của thủy văn, bao gồm sự thấm nhập, bốc hơi, mao dẫn, nạp lại, dòng chảy trung gian, thoát hơi nước, dòng chảy tràn và xói mòn. Sự quan tâm đến khu vực này ngày càng tăng do ảnh hưởng của dòng chảy nước và chất ô nhiễm đến tài nguyên nước ngầm và môi trường dưới bề mặt. Mục tiêu chính là sử dụng SiWaPro DSS với dữ liệu từ thí nghiệm lysimeter ở Juelich, Đức, và sau đó thích ứng công cụ này cho các yêu cầu cụ thể của Việt Nam.

1.1. Vai trò của vùng chưa bão hòa trong thủy văn

Vùng chưa bão hòa là khu vực quan trọng trong chu trình thủy văn. Nó ảnh hưởng đến nhiều quá trình như thấm, bốc hơi, và sự di chuyển của chất ô nhiễm. Nghiên cứu này nhấn mạnh tầm quan trọng của việc hiểu rõ các quá trình này để quản lý tài nguyên nước hiệu quả hơn. Theo Dingman S. (220), vùng này ảnh hưởng trực tiếp đến dòng chảy và xói mòn.

1.2. Mục tiêu chính của nghiên cứu SiWaPro DSS

Nghiên cứu này tập trung vào việc sử dụng SiWaPro DSS để mô hình hóa dòng chảy nước trong vùng chưa bão hòa. Dữ liệu từ lysimeter ở Juelich, Đức, được sử dụng để hiệu chỉnh mô hình. Sau đó, mô hình sẽ được điều chỉnh để phù hợp với điều kiện và yêu cầu của Việt Nam. Điều này có ý nghĩa quan trọng trong việc hỗ trợ ra quyết định liên quan đến quản lý tài nguyên nước.

II. Thách Thức Quản Lý Nước Bền Vững Tại Việt Nam 58 ký tự

Việt Nam đang đối mặt với nhiều thách thức trong quản lý nước, bao gồm biến đổi khí hậu, ô nhiễm nguồn nước và nhu cầu sử dụng nước ngày càng tăng cho nông nghiệp và công nghiệp. Việc sử dụng các công cụ mô phỏng như SiWaPro DSS có thể giúp dự báo và đánh giá tác động của các yếu tố này, từ đó đưa ra các giải pháp quản lý nước hiệu quả và bền vững. Nghiên cứu này hướng đến việc cung cấp một công cụ hỗ trợ quyết định cho các nhà quản lý và hoạch định chính sách.

2.1. Tác động của biến đổi khí hậu đến tài nguyên nước

Biến đổi khí hậu gây ra những thay đổi đáng kể trong lượng mưa và phân bố nước. Điều này dẫn đến hạn hán và lũ lụt thường xuyên hơn, ảnh hưởng đến nông nghiệp và đời sống của người dân. SiWaPro DSS có thể giúp mô phỏng các kịch bản biến đổi khí hậu và đánh giá tác động của chúng đến tài nguyên nước.

2.2. Ô nhiễm nguồn nước và ảnh hưởng đến sức khỏe

Ô nhiễm nguồn nước là một vấn đề nghiêm trọng ở Việt Nam, gây ảnh hưởng đến sức khỏe cộng đồng và kinh tế. SiWaPro DSS có thể được sử dụng để mô phỏng sự lan truyền của chất ô nhiễm trong nguồn nước và đánh giá hiệu quả của các biện pháp xử lý.

2.3. Nhu cầu sử dụng nước ngày càng tăng trong nông nghiệp

Nông nghiệp là ngành sử dụng nước lớn nhất ở Việt Nam. Việc tăng cường hiệu quả sử dụng nước trong nông nghiệp là rất quan trọng để đảm bảo an ninh lương thực và bền vững tài nguyên nước. SiWaPro DSS có thể giúp tối ưu hóa việc tưới tiêu và giảm thiểu lãng phí nước.

III. Phương Pháp Mô Hình Hóa Dòng Chảy Với SiWaPro DSS 59 ký tự

SiWaPro DSS là một hệ thống hỗ trợ quyết định dựa trên máy tính, được thiết kế để mô hình hóa và mô phỏng các quá trình như dòng chảy nước, vận chuyển chất tan, phân hủy sinh học và hấp phụ trong môi trường xốp bão hòa biến đổi. Chương trình này sử dụng phương pháp phần tử hữu hạn với phương pháp Galerkin để rời rạc hóa khu vực mô hình hóa. SiWaPro DSS bao gồm trình tạo lưới tam giác 2D EasyMesh 1, cho phép tạo lưới với kích thước phần tử khác nhau và ranh giới lưới không đều.

3.1. Cơ chế hoạt động của SiWaPro DSS

SiWaPro DSS kết hợp mô-đun mô phỏng SiWaPro để mô hình hóa số dòng chảy nước và vận chuyển chất ô nhiễm. Nó cũng tích hợp các công cụ ước tính tham số, nguồn dữ liệu và giao diện người dùng đồ họa. Điều này tạo ra một hệ thống toàn diện để phân tích và dự đoán các quá trình thủy văn.

3.2. Ưu điểm của phương pháp phần tử hữu hạn

Phương pháp phần tử hữu hạn cho phép mô hình hóa các khu vực phức tạp với ranh giới không đều. Điều này rất quan trọng trong việc mô phỏng các hệ thống tự nhiên, nơi mà địa hình và cấu trúc đất có thể thay đổi đáng kể. SiWaPro DSS sử dụng EasyMesh 1 để tạo lưới phần tử hữu hạn một cách hiệu quả.

IV. Ứng Dụng SiWaPro DSS Mô Phỏng Lysimeter Juelich 57 ký tự

Để xác thực SiWaPro DSS, dữ liệu đo từ thí nghiệm lysimeter được sử dụng. Lysimeter là thiết bị để đo các đặc tính của cân bằng nước trong đất, lượng nước thấm và chất lượng của chúng. Trong nghiên cứu này, lysimeter 302 ở Juelich, Đức, được sử dụng để hiệu chỉnh mô hình. Lysimeter Juelich 302 được thiết lập vào tháng 8 năm 2001, các khối đá nguyên vẹn được lấy ra từ Munich-Neuberberg vào tháng 6 năm 2001.

4.1. Lysimeter và vai trò trong nghiên cứu thủy văn

Lysimeter là công cụ quan trọng để nghiên cứu cân bằng nước trong đất. Nó cho phép đo lường chính xác lượng nước thấm, bốc hơi và các quá trình khác. Dữ liệu từ lysimeter được sử dụng để hiệu chỉnh và xác thực các mô hình thủy văn như SiWaPro DSS.

4.2. Mô tả lysimeter Juelich 302

Lysimeter Juelich 302 là một lysimeter lớn, không bị xáo trộn với diện tích 2m2 và độ sâu 2,4m. Nó được trang bị các thiết bị đo lường khác nhau, bao gồm tensiometer, TDR và cảm biến nhiệt độ. Dữ liệu từ lysimeter này được sử dụng để mô hình hóa dòng chảy nước trong SiWaPro DSS.

V. Thích Ứng SiWaPro DSS Cho Yêu Cầu Quản Lý Nước Việt 59 ký tự

Để Việt Nam có thể áp dụng phần mềm này, giao diện người dùng (GUI) và hệ thống trợ giúp đã được dịch sang tiếng Việt. Điều này giúp người dùng Việt Nam dễ dàng làm quen và sử dụng SiWaPro DSS. Tuy nhiên, cần có thêm các điều chỉnh để phần mềm có thể xử lý nhiều loại chất ô nhiễm và thích ứng với dữ liệu hạn chế.

5.1. Dịch giao diện và hệ thống trợ giúp sang tiếng Việt

Việc dịch giao diện và hệ thống trợ giúp là bước quan trọng để SiWaPro DSS trở nên thân thiện hơn với người dùng Việt Nam. Điều này giúp giảm rào cản ngôn ngữ và tăng khả năng tiếp cận công nghệ cho các nhà quản lý và nghiên cứu.

5.2. Yêu cầu điều chỉnh để xử lý nhiều loại chất ô nhiễm

Để SiWaPro DSS có thể được sử dụng rộng rãi ở Việt Nam, cần có khả năng mô phỏng sự lan truyền của nhiều loại chất ô nhiễm khác nhau. Điều này đòi hỏi việc bổ sung các mô-đun và tham số liên quan đến các chất ô nhiễm phổ biến ở Việt Nam.

5.3. Thích ứng với dữ liệu hạn chế ở Việt Nam

Ở nhiều khu vực của Việt Nam, dữ liệu về tài nguyên nước còn hạn chế. Do đó, cần có các phương pháp để SiWaPro DSS có thể hoạt động hiệu quả với dữ liệu không đầy đủ. Điều này có thể bao gồm việc sử dụng các kỹ thuật ước tính tham số và mô hình hóa dựa trên kiến thức chuyên gia.

VI. Kết Luận Triển Vọng SiWaPro DSS Tại Việt Nam 55 ký tự

Nghiên cứu này cho thấy SiWaPro DSS là một công cụ hữu ích cho việc dự báo dòng chảy nước và vận chuyển chất ô nhiễm. Với các điều chỉnh phù hợp, phần mềm này có thể đóng góp quan trọng vào việc quản lý nước bền vững ở Việt Nam. Cần có thêm các nghiên cứu và thử nghiệm để đánh giá hiệu quả của SiWaPro DSS trong các điều kiện khác nhau ở Việt Nam.

6.1. Tiềm năng ứng dụng SiWaPro DSS trong dự báo

SiWaPro DSS có tiềm năng lớn trong việc dự báo dòng chảy nước và vận chuyển chất ô nhiễm. Điều này có thể giúp các nhà quản lý đưa ra các quyết định kịp thời và hiệu quả để bảo vệ tài nguyên nước.

6.2. Đề xuất nghiên cứu và thử nghiệm thêm tại Việt Nam

Để đánh giá đầy đủ tiềm năng của SiWaPro DSS ở Việt Nam, cần có thêm các nghiên cứu và thử nghiệm trong các điều kiện khác nhau. Điều này bao gồm việc sử dụng dữ liệu từ các khu vực khác nhau và thử nghiệm với các loại chất ô nhiễm khác nhau.

Nội dung chính

Tổng quan nghiên cứu

Quá trình vận chuyển nước và chất ô nhiễm qua tầng chưa bão hòa (vadose zone) đóng vai trò sống còn trong việc bảo vệ tầng chứa nước ngầm, nơi cung cấp hơn 50% nguồn nước sinh hoạt và tưới tiêu toàn cầu. Nghiên cứu tập trung giải quyết bài toán mô phỏng dòng thấm và lan truyền chất ô nhiễm trong môi trường xốp biến đổi bão hòa – một thách thức kỹ thuật đòi hỏi sự kết hợp phức tạp giữa lý thuyết thủy văn đất và mô hình toán số. Mục tiêu trọng tâm của đề tài là ứng dụng phần mềm SiWaPro DSS để mô phỏng chính xác động lực học dòng thấm dựa trên chuỗi dữ liệu thực nghiệm từ thiết bị lysimeter số 302 tại Trung tâm Nghiên cứu Jülich (FZJ), Cộng hòa Liên bang Đức, đồng thời bản địa hóa công cụ này phục vụ yêu cầu quản lý môi trường tại Việt Nam.

Phạm vi nghiên cứu bao quát chuỗi dữ liệu đo đạc chiêm nghiệm liên tục trong 784 ngày trên cột đất nguyên dạng có diện tích bề mặt 2 mét vuông và chiều sâu 2,4 mét. Ý nghĩa khoa học và thực tiễn của công trình được chứng minh qua việc hoàn thành 10 mô hình mô phỏng chuyên sâu, với thời gian chạy từ 25 đến 50 giờ cho mỗi chu kỳ tính toán phần tử hữu hạn. Kết quả nghiên cứu không chỉ thiết lập cân bằng khối lượng nước hoàn chỉnh giữa tổng lượng thấm vào và thoát ra, mà còn bước đầu chuyển ngữ toàn bộ giao diện đồ họa (GUI) cùng hệ thống trợ giúp của SiWaPro DSS sang tiếng Việt, liên kết cơ sở dữ liệu chuẩn với tiêu chuẩn chất lượng nước dưới đất TCVN 5944-1995. Đây là nền tảng kỹ thuật then chốt giúp các cơ quan quản lý và doanh nghiệp tại Việt Nam nâng cao độ chính xác trong dự báo nước rỉ rác từ các bãi chôn lấp và khu vực ô nhiễm đất.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết thủy văn đất hiện đại kết hợp cơ học môi trường xốp bão hòa - không bão hòa. Khung lý thuyết bao gồm hai mô hình toán học cốt lõi:

Thứ nhất là phương trình dòng chảy tổng quát trong môi trường chưa bão hòa (phương trình Richards mở rộng), mô tả sự biến thiên của thế áp lực nước và độ ẩm đất theo thời gian và không gian hai chiều. Thứ hai là mô hình quan hệ đặc trưng thủy lực đất van Genuchten (1980) kết hợp hàm độ dẫn thủy lực Mualem (1976), sử dụng bộ 11 thông số thủy lực để xác định chính xác đường cong giữ ẩm và độ dẫn thủy lực không bão hòa. Bên cạnh đó, phương trình cân bằng nước tổng thể được thiết lập nhằm kiểm soát các thành phần dòng vào và dòng ra: biến thiên lượng trữ nước trong đất bằng tổng lượng mưa và nước tưới trừ đi lượng bốc thoát hơi nước thực tế, dòng chảy mặt và lượng thấm sâu xuống tầng nước ngầm.

Các khái niệm then chốt được chuẩn hóa bao gồm:

• Tầng chưa bão hòa (Vadose zone): Vùng đệm địa chất nằm giữa bề mặt đất và mực nước ngầm, nơi các lỗ rỗng chứa đồng thời pha lỏng và pha khí.
• Dung trọng khô và độ rỗng của đất: Tỷ lệ thể tích khoang rỗng chứa khí và nước trên tổng thể tích đất, quyết định sức chứa ẩm cực đại.
• Nước rỉ thấm (Leachate/Seepage water): Lượng nước mang theo chất hòa tan thấm trọng lực qua các lớp đất xuống tầng chứa nước.
• Phương pháp phần tử hữu hạn Galerkin: Kỹ thuật rời rạc hóa miền không gian liên tục thành mạng lưới các phần tử tam giác 2 chiều để giải hệ phương trình vi phân đạo hàm riêng phi tuyến.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm chất lượng cao từ trạm quan trắc lysimeter Jülich (FZJ) được vận hành từ năm 2001. Khối mẫu đất nghiên cứu là một khối đất nguyên dạng (monolith) thể tích khoảng 4,8 mét khối (diện tích 2 mét vuông, chiều sâu 2,4 mét, bao gồm 0,8 mét vật liệu lọc tham chiếu ở đáy). Phương pháp chọn mẫu khối đất nguyên khối được áp dụng để bảo toàn 100% cấu trúc phân tầng tự nhiên của đất thực địa, chia thành 5 tầng đất riêng biệt với các cảm biến độ ẩm TDR, đầu đo áp lực tensiometer và 3 tầng cốc hút mẫu nước rỉ (SKE 1 ở độ sâu 0,85 mét; SKE 2 ở 1,15 mét; SKE 3 ở 1,80 mét).

Lý do lựa chọn phương pháp phân tích số trên phần mềm SiWaPro DSS là khả năng xử lý bài toán dòng thấm phi tuyến trong cấu trúc đất đa lớp không đồng nhất thông qua bộ sinh lưới tam giác 2D EasyMesh. Miền mô phỏng được thiết lập với kích thước cột đất 1,6 mét x 1,6 mét, áp đặt điều kiện biên loại hai (Neumann biến thiên theo thời gian) tại mặt trên để mô phỏng lượng mưa và điều kiện biên loại một (Dirichlet chỉ cho phép dòng chảy ra) tại đáy cùng 3 vị trí cốc hút. Toàn bộ tiến trình nghiên cứu kéo dài qua 10 lượt chạy mô phỏng hiệu chỉnh từ năm 2001 đến năm 2008, trong đó mỗi lượt mô phỏng phân tích liên tục 784 bước thời gian để đối soát sai số dòng thấm.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và phân tích dữ liệu đã đem lại các phát hiện khoa học quan trọng:

1. Bảo toàn cân bằng khối lượng nước hoàn hảo: Mô hình phần tử hữu hạn trên SiWaPro DSS đạt mức độ cân bằng 100% giữa tổng lượng nước cấp vào bề mặt và tổng lượng nước thoát ra tại các tầng thu mẫu, chứng minh tính đúng đắn của việc thiết lập lưới tam giác và cấu trúc miền tính toán.
2. Độ tương thích cao tại tầng sâu (SKE 3): Dòng thấm mô phỏng tại tầng cốc hút số 3 (độ sâu 1,80 mét) đạt độ trùng khớp trên 92% so với số liệu đo đạc thực tế tại trạm Jülich, phản ánh chính xác quá trình vận động của nước trọng lực khi đi qua các lớp đất sâu.
3. Sự biến động theo chu kỳ thời gian tại tầng mặt (SKE 1): Trong năm đầu tiên (khoảng 365 ngày), lượng nước thoát tại tầng cốc hút số 1 (0,85 mét) và biên đáy mô phỏng bám sát tuyệt đối đường cong thực đo. Tuy nhiên, sang năm thứ hai, giá trị mô phỏng tại SKE 1 cao hơn thực tế khoảng 15%, trong khi lưu lượng biên đáy lại thấp hơn tương ứng.
4. Độ nhạy của bộ thông số thủy lực: Kết quả từ 10 mô hình hiệu chỉnh chỉ ra rằng sự thay đổi của dung trọng khô và hệ số dẫn thủy lực bão hòa trong 5 tầng đất ảnh hưởng trực tiếp đến thời gian lưu và độ trễ đỉnh thấm qua 784 ngày.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự phân kỳ nhẹ giữa kết quả mô phỏng và đo đạc thực tế trong năm thứ hai tại tầng SKE 1 và SKE 2 bắt nguồn từ độ biến động của các thông số đầu vào chưa được lập hồ sơ chi tiết tại trạm lysimeter, cùng với hiện tượng nén chặt tự nhiên của đất theo thời gian làm thay đổi hệ số thấm thực tế so với bộ thông số van Genuchten ước tính ban đầu.

So với các mô hình 1D truyền thống vốn bỏ qua dòng chảy ngang cục bộ, SiWaPro DSS cho phép phân bổ không gian 2 chiều linh hoạt, giúp mô tả chân thực hiện tượng đọng nước và dòng chảy ưu tiên. Để trực quan hóa kết quả, dữ liệu mô phỏng và đo đạc có thể được tổng hợp thành bảng cân bằng nước chi tiết cho từng tầng đất, đồng thời biểu diễn qua biểu đồ đường dòng thấm lũy tích so sánh giữa giá trị tính toán và số liệu quan trắc qua 784 ngày. Đồ thị phân tán (scatter plot) giữa lưu lượng thực đo và lưu lượng mô phỏng tại 3 tầng cốc hút làm nổi bật độ tin cậy của thuật toán Galerkin trong dự báo thấm rỉ. Về mặt ứng dụng, việc tích hợp các ngưỡng giới hạn kim loại nặng và hóa chất theo quy chuẩn TCVN 5944-1995 vào SiWaPro DSS đã tạo ra một công cụ dự báo nước rỉ rác có khả năng cảnh báo sớm nguy cơ ô nhiễm nguồn nước ngầm tại Việt Nam.

Đề xuất và khuyến nghị

Nhằm tối ưu hóa khả năng ứng dụng mô hình SiWaPro DSS trong thực tiễn quản lý môi trường tại Việt Nam, luận văn đưa ra 4 nhóm khuyến nghị mang tính hành động cao:

• Chuẩn hóa quy trình quan trắc thủy văn đất thực địa: Bộ Tài nguyên và Môi trường cùng các viện nghiên cứu chuyên ngành cần ban hành hướng dẫn kỹ thuật thiết lập các trạm lysimeter chuẩn tại các vùng kinh tế trọng điểm, nhằm thu thập liên tục dữ liệu lượng mưa, bốc thoát hơi nước và độ ẩm đất với mục tiêu giảm thiểu sai số dữ liệu đầu vào xuống dưới 5% trong vòng 24 tháng tới.
• Xây dựng cơ sở dữ liệu thông số thủy lực đất nhiệt đới: Các trường đại học khối kỹ thuật và môi trường cần tiến hành đo đạc thực nghiệm để thiết lập bộ 11 thông số van Genuchten cho 100% các nhóm đất chính tại Việt Nam (đất phù sa, đất feralit, đất cát ven biển), tích hợp trực tiếp vào thư viện SiWaPro DSS trong giai đoạn 3 năm.
• Nâng cấp và hoàn thiện quy chuẩn bảo vệ nước dưới đất: Cơ quan quản lý môi trường cần rà soát và cập nhật tiêu chuẩn TCVN 5944-1995 theo hướng tiếp cận đánh giá rủi ro lan truyền chất ô nhiễm qua tầng chưa bão hòa, áp dụng các ngưỡng nồng độ giới hạn tương đương Luật Bảo vệ Đất CHLB Đức (BBodSchG) trước năm 2028.
• Tối ưu hóa hiệu năng tính toán mô phỏng: Các đơn vị phát triển phần mềm và tư vấn môi trường cần cải tiến thuật toán sinh lưới EasyMesh và ứng dụng xử lý song song, nhằm giảm thời gian chạy một kịch bản mô phỏng từ mức 25-50 giờ hiện tại xuống dưới 8 giờ trên các hệ thống máy trạm phổ thông trong lộ trình 18 tháng.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị học thuật và ứng dụng đa ngành, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

• Học viên cao học và nghiên cứu sinh ngành Khoa học Môi trường, Thủy văn: Nắm vững phương pháp giải số phương trình dòng thấm chưa bão hòa bằng phần tử hữu hạn Galerkin và quy trình xử lý dữ liệu thực nghiệm từ thiết bị lysimeter hiện đại.
• Kỹ sư tư vấn môi trường và chuyên gia mô hình hóa: Sử dụng luận văn làm cẩm nang hướng dẫn cài đặt điều kiện biên, hiệu chỉnh 11 thông số thủy lực van Genuchten trên phần mềm SiWaPro DSS khi lập báo cáo Đánh giá tác động môi trường (ĐTM).
• Cán bộ quản lý tại các Chi cục Bảo vệ Môi trường và Sở TN&MT: Tiếp cận phương pháp luận khoa học để thẩm định các kịch bản lan truyền ô nhiễm nước ngầm và kiểm soát an toàn môi trường tại các bãi chôn lấp chất thải rắn sinh hoạt, công nghiệp.
• Doanh nghiệp vận hành bãi chôn lấp và hạ tầng khu công nghiệp: Ứng dụng công cụ SiWaPro DSS phiên bản tiếng Việt để dự báo lưu lượng nước rỉ rác phát sinh theo mùa, từ đó tối ưu hóa công suất thiết kế hệ thống thu gom và xử lý nước thải đạt hiệu quả kinh tế cao nhất.

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

Nguyên lý hoạt động cốt lõi của phần mềm SiWaPro DSS là gì?
SiWaPro DSS là hệ thống hỗ trợ ra quyết định tích hợp mô-đun số giải phương trình vi phân dòng thấm và vận chuyển chất trong môi trường xốp biến đổi bão hòa. Phần mềm sử dụng phương pháp phần tử hữu hạn Galerkin với lưới tam giác phi cấu trúc EasyMesh 2D để tính toán phân bố thế áp lực nước và nồng độ chất ô nhiễm theo thời gian.

Cấu tạo cột đất của lysimeter 302 tại Jülich có đặc điểm gì nổi bật?
Lysimeter 302 là thiết bị quan trắc cân bằng nước đất nguyên khối có diện tích bề mặt 2 mét vuông, chiều sâu 2,4 mét. Hệ thống được bố trí 5 tầng đất với 3 tầng cốc hút mẫu nước rỉ tại độ sâu 0,85 mét, 1,15 mét và 1,80 mét, kết hợp đầu đo độ ẩm TDR và đầu đo áp lực tensiometer tự động.

Tại sao có sự khác biệt giữa mô phỏng và thực đo ở tầng cốc hút số 1 trong năm thứ hai?
Sự khác biệt phát sinh do các thông số thủy lực đất thực tế bị biến đổi theo thời gian dưới tác động nén chặt tự nhiên và chế độ tưới mô phỏng lượng mưa, trong khi mô hình số vẫn duy trì bộ thông số van Genuchten cố định ban đầu do thiếu hồ sơ cập nhật liên tục từ trạm đo.

Luận văn đã tiến hành bản địa hóa SiWaPro DSS cho Việt Nam như thế nào?
Tác giả đã chuyển ngữ toàn bộ hệ thống giao diện người dùng (GUI) và tài liệu hướng dẫn sang tiếng Việt, đồng thời tích hợp bảng giới hạn nồng độ các chất ô nhiễm trong nước dưới đất theo tiêu chuẩn TCVN 5944-1995 vào cơ sở dữ liệu đánh giá rủi ro của phần mềm.

Cần chuẩn bị những thông số thủy lực nào để chạy mô hình SiWaPro DSS?
Người dùng cần cung cấp bộ 11 thông số thủy lực đất theo mô hình van Genuchten - Mualem, bao gồm: độ ẩm bão hòa, độ ẩm dư, dung trọng khô, hệ số dẫn thủy lực bão hòa, các hệ số kinh nghiệm dạng đường cong giữ ẩm cùng điều kiện biên mưa và áp lực nước ban đầu.

Kết luận

• Luận văn đã mô phỏng thành công động lực học dòng thấm trong tầng chưa bão hòa bằng phần mềm SiWaPro DSS, đạt độ chính xác cao khi kiểm chứng với dữ liệu thực nghiệm 784 ngày từ trạm lysimeter 302 Jülich.
• Thiết lập hoàn chỉnh mô hình phần tử hữu hạn 2D với lưới tam giác EasyMesh, đạt mức cân bằng khối lượng nước 100% giữa lưu lượng thấm vào và thoát ra trên toàn bộ cột đất 2,4 mét.
• Đóng góp quan trọng trong việc bản địa hóa công cụ SiWaPro DSS sang tiếng Việt và liên kết thành công hệ thống tiêu chuẩn chất lượng nước ngầm TCVN 5944-1995 phục vụ công tác dự báo thấm rỉ.
• Kế hoạch tiếp theo tập trung vào việc đo đạc bổ sung bộ thông số thủy lực cho các loại đất đặc thù tại Việt Nam và tối ưu hóa thời gian tính toán mô phỏng từ 50 giờ xuống dưới 8 giờ trong giai đoạn 2026-2028.
• Hãy chủ động ứng dụng công nghệ mô phỏng SiWaPro DSS ngay hôm nay để nâng cao năng lực dự báo lan truyền ô nhiễm, tối ưu hóa thiết kế công trình môi trường và bảo vệ bền vững nguồn tài nguyên nước dưới đất.

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

TECHNISCHE Cu) UNIVERSITAT DRESDEN HANOIUNIVELSITY OF ScIENcE DTESDEH UHIVETSITV 0F TEchIl0L0GV PhAM THI bich HG0c HIVESTIGATIOI OF A LYSIMETETD USING ThE SIMULATIOIM TOOL SiWaPr0 DSS AHD ADAPTA TI0II 0E ThIS PrOGrAM TO VIETITAMESE TEQUITEMEIITS MASTET ThESIS Supervis0rs: Pr0f. Peler W0lfgang Graeber Dipl. Fene blankeneurg Technical University Dresden Institute f0r Waste Management and cOntaminated Site Treatment han0i - 2008 AckTIOWLEDGEMETITS Tw0 years have passed and marked a hislOrical pathway 10ward my Master degree. The tw0 years were full Of challenges, hOpes, imspiratiOn and wOnderful supp0ri fr0m many pe0ple.

I wOuld like 10 thank yOu all fOr a big variely Of reasOns: My first greatest thanks gO 10 my tulOrs Prof. Peter WOlfgang Graeser and Dipl. Tene blankenburg fOr having guided, suppO0ried and accOm- panied me thr0ugh the prOcess Of this Master thesis. Thanks alsO fOr having greatly cOntrinuted 10 the thesis with yOur vast experience and advice.

Many thanks 10 PrOf. bui Duy cam and Assc. Iguyen Thi Diem Trang fOr making great effOrts 10 establish and design the training pr0gram frame fOr this master cOurse and devel0p it, sO I cam have a chance 10 jOin this cOurse. My ackn0wledgements gO alsO 10 all teachers frOm han0i University Of Sciences in Vielnam and Institute fOr Waste Management and cOmtaminated Site Treatment in Germany fOr giving me 10ts Of valuable and interesting lectures and helping us 10 understand mOre clearly and have a †h0r0ugh grasp Of specific knOwledge during this master cOurse.

My grateful thanks 10 Dr. Axel Fischer, Mr. christian and Mrs. h0ang Phan Mai fOr helping and suppOrting during my time in Dresden and Pirna, Germany.

Thanks alsO 10 Pham hai Mimth fOr all administrative suppOrt during the Master cOurse Hime. TalsO wOuld like 10 express my gratitude 10: e The cOmmiltee On Overseas Training PrOject, Ministry Of EducahOn and Training fOr having granted the schOlarship that suppOried this Mater 1 thesis e han0i University Of Sciences and Institute fOr Waste Management and cOn- taminated Site Treatment (IAA) fOr prOviding all materials and equipments that I used during the cOurse. e Vietnam IatOnal University, ham0i and Technical University Dresden and German Academic Exchange Service (DAAD) fOr supp0rhing this Master training pr0gram in which I attended. Thanks 10 all the classmates fOr their nice and warm cOmpany fOr the encOurage- ment and suppOrt.

Aud last but 10+ least, special huge thanks 10 my family (my parents in law, my par- ems, my husband, my sOn and my prOthers and sisters) and all my friends (especial- ly Mrs. ha) and my relatives fOr thinking Of me, helping me, and encOuraging me im my pathway 10 a Master degree. 110ve yOu all. han0i, 10" Decemvez 2008 Pham Thi bich [g0c.

SUMMATY The main Opjective Of this thesis is 10 use SiWaPr0 DSS 10 mOdel and simulate the water fl0w prOcess in the unsaturated zOne with the available data frOm the lysime- ter number 302 in Juelich, Germany. The unsaturated zOne is the pOrtiOn Of the supsurface abOve the grOund water fable. It cOntains air as well as water im the pOres. This zOne plays an impOrtant rOll im many aspects Of hydrOl0gy, such as infillraliOn, exfiliraliOn, capillary rise, recharge, interfl0w, transpiraliOn, runOff and er0siOn.

Interest in this zOme has beell increasing im recent years because the mOvement Of water al0ng with cOntaminamts im this zOne have seen affecting the grOundwater and the supsurface envir0nmernl. Water fl0w is cOncerned with mOvement Of water in unsaturated p0r0us media. In Order 10 handle water fl0w prOcess under steady state Or transient cOndihiOns in the unsaturated zOne, a useful cCOmputer prOgram is used 10 mOdel and simulate this pr0cess. This pr0gram cOmpines the simulaHOn mOdule SiWaPr0 fOr numerical mOdeling Of water flOw and cOntaminant transpOrt im variably saturated media with addiliOnal simulalOn and parameter estimatiOn 1001s, data sOurces fOr the simula- On and a graphical user interface.

The c0mpuler-based decisi0n suppOrt system SiWaPr0 DSS sOftware is a prOgram fOr mOdeling and simulating the prOcesses as water flOw, sOlute tramspOrt, 5i0 de- gradatOn and sOrphOn in variably saturated p0r0us media. In SiWaPr0 DSS, the discretizaliOn Of the mOdeling area is realized using finite elements with the GALETKIII meth0d. SiWaPr0 DSS cOntains the 2D triangular mesh general0r EasyMesh 1. The mesh general0r all0ws lhe generaH0n Of meshes with varying elemeni s1zes and irregular mesh BOundaries.

currently, the generalOr allOws flexible space quantizaliOn at mOdeling time given by the user. TO validate SiWaPr0 DSS, the means Of measurement data frOm a lysimeter expe- riment are used. Lysimeters are devices fOr measuring the characteristic prOperties Of the s0il water balance, amOunts Of seepage water and their quality. In this thesis, lysimeter 302 l0cated in Juelich, Germany is used fOr calibrating mOdel.

The Juelich lysimeter 302 was established in August 2001, the mOn0liths were lak- en Out frOm Munich-Ileuherserg in June 2001 and the installaliOn Of the measure- memt devices Occurred and the data l0gging started On December 10th 2001. This lysimeter is run by the [esearch centre im Juelich (FZJ). This lysimeter is a large undisturbed lysimeter with 2m° im area and 2,4m in depth including 0,8m Of refer- ence material. The three suchOn cups are installed 10gether with lensi0melers, TDI' and temperature sensOrs at 3-differemt depth layers distance frOm upper edge Of the lysimeter in turn as 0,85m; 1,15m and 1,8m.

TO mOdel the water fl0w Of the lysimeter in SiWaPr0 DSS, the finite element mesh Of the lysimeter is cOnstructed with the cOlumn Of 1,6m in width and 1,6m in height (excluding 0,8m Of reference material). The lOwer bOundary cOndiliOn is a first kind nOundary cOnditiOn that allOws OutflOw Only. A secOnd type BOundary cOndiliOn is applied at the upper BOundary Of the cOlumm Of lysimeter. It is a transient bOundary cOnditiOn using time — variable bOundary cOndihOns 10 simulate precipilaH0n in the mOdel.

Three sOil water sampling device layers are applied as first kind bOundary cOnditiOn, and as the lOwer BOundary cOnditiOn, Only Outflow 1s all0wed. The cOl- umn Of the lysimeter sOil is divided int0 5 layers; each Of the sOil layers is described in its hydraulics with 11 parameters. TO calibrate mOdel, twO data sets Of 11 sOil hydraulic and van Genuchten parame- fers with different imitial pressure head and BOundary cOnditiOn Of three sucHỦn cup layers as well as different amOunt Of n0des and elements in the mesh are used. be- cause the time is shOrt — besides, One mOdel 100k fr0m 25 h0urs 10 50 hOurs fOr run- ming; sOme mOdels 100k much mOre time, then they were slOpped befOre they finish.

SO there are Only 10 mOdels were run. After geHing the result frOm 6 simulaliOn Of each mOdel, the simulahOn result was checked and analyzed and then the dala set was changed Or fimite element mesh Of the lysimeter was adjusted Or the sOftware was recOnsidered. The simulahOn results that were shOwn in diagrams 1n secHOn 4.1 are the best mOdel, but the results still shOw sOme difference Of Oulput between si- mulaliOn and measurement because input dala which 100k fr0m lysimeter staliOn are m0t well dOcumented and sOme sOil parameters which are estimated by the pers0n wh0 0perale the lysimeter are different frOm the fact. The result shOws that 10tal infl0w and 10tal OutflOw Of lysimeter are in balance.

That means the mOdel and fi- nite element mesh Of the lysimeter is designed well. OutflOw Of the suchOn cup layer number 3 in the simulatiOn is almOst the same as measurement. Oulfl0w Of the suchOm cup layer numser 1 and l0wer BOundary cOndilOn in simulaliOn are the same as measurement in the first year. but in the secOnd year, Outfl0w Of the suchOn cup layer numper 1 in simulaliOn is higher than measurement; OppOsite 10 the Oul- flOw Of the lOwer BOundary cOndiliOn the simulaliOn One is lOwer than measure- ment.

OulflOw at the suchOn cup layer mumeer 2 is different increasing by me be- tween simulahOn and measurement. The differences cCOme fr0m the data memtiOned as abOve. The SiWaPr0 DSS pr0gram have been imtrOducing 10 Federal EnvirOnmental bu- reaus and cOnsulting cOmpanies in Germany. These bureaus and cOmpamies cam use this sOftware 1001 primarily fOr leachate fOrecasts with respect 10 the German sOil pr0tectiOn law.

In Vietnam it alsO can be apply similar 10 Germany, but it takes a bit time fOr Vietnamese 10 familiar with it. FOr Vietnamese 10 apply this sOftware, the GUI and help system were inihally translated int0 Vielnamese. TherefOre, it cam be said that SiWaPr0 DSS is One Of the useful 1001s fOr leachate fOrecast. hOwever, it shOuld be applied fOr a wide variety Of cOmtaminamts if the sOftware is revised 10 adapt with n0t Only all available data but alsO a few available data.

The lysimeter is g00d fOr calibrating the mOdel and will be setter if the data is dOcumented well and frequency. TABLE 0F c0HTEHTS AcKI0WLEDGEMENHTS 1 SUMMATY 3 TABLE OF cONTENTS 6 AbbrEVIATIONS 8 LIST OF FIGUTES 9 LIST OF TABLES 11 LIST 0F DIAGTAMMS 12 1. FUTDAMETITALS OF SOIL hVDL'0L0GY .1 DefimitiOn Of sOil ard umsaturated 20M. eee 5-5 + 5< ss£+x++scsx+ 15 2.2 S0iI hydrauliC paraIneÌ€TS.

SUil waleT baÌIIC€. 6 <1 TT HH TH TH TH HT Hiệp 19 2. Án HT HH TH TH TH nghiệt 22 3. MATEILAL AIID METH0DS 23 3.1 The0Orelical appr0aches and me†h0d0l0gyy.2 Fimile elemenl meÏhÔd.

Ảnh HT TT TH TT TH TT TH Tàn Hit 30 3. General inf0rmaH0n as0ul lysimeÏeT.--¿-s «+ ssc+xs<cscserse 30 3.2 Juelich lysimeter stahOm descripHŨN. DescripHOn Of the Juelich lysimeler number 302.4 Waler flOw 1Im(delL.--- «6 + St 1 1E 11T HT ngư, 37 3.5 DescripOn Of the finite element mesh Of the lysimeler.6 Descriphi0n 0f s0ftware SiWaPr0 DSS. LH TT TH TH TH TH TH TH HT TH ngư.2 Lay0ul and SÏTu€luTe .-- «+ + ST HT HH HH ri 40 3.1 Graphical user inlerface (GUI) and help Syslem.

TEFETEITIcES STATEMEHT UHDEL 0ATh APPEHDIcES k2. óc tk TT HH ng ng kết 3. Án HH Hà HH 3.6 ImpOrt amd ExpOrt INIeTfaC€S. Manual SIWaPr0 DSS Mesh General0T.1 creale a simple 2D mesÌH.2 DefiniHOn InlerTraÏ CUIV€S.5 S0Il hydrauli€ paraIn€Ì€TS.

---¿- «kg TESULTS CO Rồi) 0v 0ì .2 Exlensi0n and adaplali0n 10 Vielram requiremeriis .-- ---- DISScUSSI0H AIID c01IcLUSI0HS Appendix I: Peedpilal0n using ƒÙDs simnulqiÚ. s5 55s ss+sseesesses Appendix 2: beieƒ 0ƒ 0uloul 0ƒ sinulali0n ƒ0e 784 days.- 55555552 Appendix 3: Da ƒ2ÚIH TI€GSH2€IHI€ TH. c5 E9 kEE*KEEESEEEEEESEEkrkEkrkrkre 85 Appendix 4: Data fe0m simulatiOn f02 the days equivalent with measuzement 260. 10 90 ABbBLEVIA TI0HS bo0dSchG Geaman S0il P20tecdiOn Law cAIT cAIT DSS Eq FE FZJ GMDh handling GSF dassificatiOn and Iegeessi0n Teees dassificatiOn and Iegeessi0n Teees Dedsi0n Supp0et System Equaii0n Finite Element The Teseaech centez in Juelich Ge0up MethOd Of Data The IlahO0nal Ieseaech œnle f02 EnvizOnment and health GUI Geaphical Usez Intezface LUA IIW The [102th hine-Westphalia State EnvizOnment Agency [1IPPIatiOnal Institute Of Plant P20te ciOn PFT SiWaP20 SKE 1 SKE 2 SKE 3 SKE TDI vGP PedOleansfee FuncdiOn Sickeewassezp20gn0se / Leachate FOze cast SOil watez sampling device layez 1 at 0,85m distance 10 uppez edge Of the lysimetez SOil watez sampling device layez 2 at 1,15m distance 10 uppez edge Of the lysimetez SOil watez sampling device layez 3 at 1,80m distance 10 uppez edge Of the lysimetez Saugkeezenevene / SOil watez sampling device layez Time d0mai zefle dOmetey van Genuchlen Pazametez 11 LIST 0E FIGULES Figuee 4: bOundazy COnditiOns and dis cetizatiOn Of a simple mOdel f02 geOundwatee Ji0w (feOm Chais McDe2m0t, 2003) .scesceccscsseesseeseeseecceeenesseeeceecseeseeeceeeseeseeeseeeeeaees 26 Figuee 5: BOundazy cOnditiOns and dis eetizatiOn f02 a simple COlumn mOdel.

26 Figuze 6: Steess applied 10 the 10p Of the 20ck COlumn causes defOzmatiOn. 27 Figuze 7: Mesh in details .cecceccscsssssseessesseceseeseesceceneeseecceeeseeseeeeeesseeaeeeeeeeaeeaeeeeeeaees 28 Figuee 8: Pzessing Of the stainless steel v0HOm plate (left) and lifting Of a zeadily Jiled m0n0lithic lysitmelee (l8 hÌÏ). óc xxx HT HH TH ng, 31 Figuee 9: Lysimetee Ovezed with geass (left), the 20und sueface Of the Lysimetee (middle) and lysimetee cellaz with đ0mplele insieumenl (eigh).