mặt-nƣớc ngầm
Trong phần này, trình bày việc thiết lập giải một số bài toán bằng cách sử dụng mô hình MODFLOW và kiểm chứng với các kết quả thu được từ lời giải giải tích ở phần trên.
Sơ đồ hệ thống được thiết lập như hình 3.1.V ng tính được giới hạn như trong hình vẽ, Sông và tầng chứa nước được ng n cách bởi một lớp có độ dẫn thủy lực nhỏ, có thể coi như một tầng chứa không áp rất nhỏ. Tầng chứa nước bên dưới sông là có áp. Các đặc tính của dòng chảy được giả thiết là đồng nhất
tại mọi nơi, d ng chảy thẳng đứng trong tầng chứa có thể bỏ qua (giả thiết Dupuit). Có một giếng hút trong khu vực tính. Sông được giả thiết là có độ sâu bằng với độ dày tầng chứa. Khu vực tính toán được chia thành 3 vùng tính toán như trên hình 3.1.
Thiết lập mô hình
Mô hình được thiết lập với lưới không đều, chiều dài không giới hạn (y=∞) và chiều rộng giới hạn. Mô hình có 1 lớp. Các thông số của mô hình được cho như sau:
- K=20 m/ngày - S=4.3x10-4 m-1
- Độ rộng của sông: w=9m - Độ dày đáy sông b‟=0,5m
- Hệ số dẫn đáy sông K‟=0,5 m/ngày - Giếng hút, lưu lượng Q=210 m3/ngày
So sánh kết quả mô hình với kết quả từ lời giải giải tích
Mô hình tính được xây dựng cho từng bài toán theo Glover và Balmer, Hantush (1965), Hunt(1999). Sau khi thiết lập mô hình với các tham số như trên, chúng tôi chạy mô hình theo thời gian, lấy các kết quả về độ suy giảm dòng chảy (stream depletion). So sánh với các kết quả thu được từ giải tích với các cách tiếp cận của Glover và Balmer (1954), Hantush (1965) và Hunt (1999). Hình 3.5 biểu diễn kết quả tính toán từ mô hình MODFLOW với kết quả tính toán từ các bài toán của Glover và Balmer (1954), Hantush (1965) và Hunt (1999). Như trên hình, ta nhận thấy, những đường liền nét là kết quả tính toán từ nghiệm giải tích, những điểm chấm là kết quả tính toán của mô hình tương ứng với từng hệ thống được mô tả trong bài toán mẫu của Glover và Balmer (1954), Hantush (1965) và Hunt (1999); kết quả bài toán chạy bằng mô hình tương đối phù hợp với kết quả nghiệm giải tích của bài toán. Là cơ sở để tiến hành tính toán với các trường hợp sau.
Hình 3.5. Đồ thị so sánh kết quả mô hình với lời giải giải tích