Một số giải pháp trong khai thác nước ngầm bằng bãi giếng nhằm giảm thiểu hạ thấp mặt đất

Bài viết này trình bày giải pháp tối ưu hóa lưu lượng nước khai thác của các giếng trong bãi giếng khai thác nhằm giảm thiểu hạ thấp mặt đất ở khu vực TP. Hồ Chí Minh. Phương pháp nghiên cứu được sử dụng là giải bài toán quy hoạch tuyến tính để tìm lời giải tối ưu lưu lượng nước khai thác trong từng giếng trong bãi giếng khai thác.

Hội thảo Khoa học Quốc tế Phát triển Xây dựng bền vững điều kiện Biến đổi khí hậu khu vực đồng Sông Cửu Long SCD2021 MỘT SỐ GIẢI PHÁP TRONG KHAI THÁC NƯỚC NGẦM BẰNG BÃI GIẾNG NHẰM GIẢM THIỂU HẠ THẤP MẶT ĐẤT SOME SOLUTIONS IN GROUNDWATER EXPLOITATION BY GOOD YARDS FOR REDUCTION LOWERING THE GROUND Nguyễn Xuân Mãn, Nguyễn Duyên Phong ABSTRACT: Groundwater exploitation leads to a lowering of the groundwater level in a non-pressurized aquifer or a decrease in pressure in a pressurized aquifer If there is no effective and appropriate remedy, it will inevitably lead to the breakdown of the stratigraphic balance; As a result, the strata subsided The subsidence process spreads to the ground causing the ground to collapse Stratigraphic subsidence and ground subsidence cause damage to underground structures and aboveground structures (underground tunnels, roads, houses and civil works, and public buildings terrestrial industry) This report presents two solutions to minimize ground subsidence caused by groundwater extraction in the Mekong Delta and Ho Chi Minh City: The first solution: Build well yards to exploit groundwater near the area with replenishing water sources The second solution: Optimizing the exploitation water flow of the wells in the mining well yard KEYWORDS: Lowering the ground, groundwater extraction, lower pressure, wells yards, exploitation water flow TÓM TẮT: Khai thác nước đất dẫn đến hạ thấp mực nước ngầm tầng nước không áp làm giảm áp lực tầng nước có áp Nếu khơng có giải pháp khắc phục hiệu phù hợp tất yếu dẫn đến phá vỡ cân địa tầng; hậu địa tầng mặt đất bị sụp lún Sụp lún địa tầng sụt lún mặt đất gây tổn hại cho cơng trình xây dựng (cơng trình ngầm mặt đất) Hiện có hai giải pháp để giảm thiểu lún mặt đất khai thác nước ngầm: Xây dựng bãi giếng khai thác nước ngầm gần khu vực có nguồn nước bổ cập; Tối ưu hóa lưu lượng nước khai thác giếng bãi giếng khai thác Báo cáo trình bày giải pháp tối ưu hóa lưu lượng nước khai thác giếng bãi giếng khai thác nhằm giảm thiểu hạ thấp mặt đất khu vực TP Hồ Chí Minh Phương pháp nghiên cứu sử dụng giải tốn quy hoạch tuyến tính để tìm lời giải tối ưu lưu lượng nước khai thác giếng bãi giếng khai thác TỪ KHÓA: Khai thác nước ngầm, sụp lún, bãi giếng, hạ áp lực, hạ thấp mặt đất Nguyễn Xuân Mãn Học hàm, học vị: PGS.TS Khoa Xây dựng, Trường Đại học Mỏ - Địa chất Tel: 0903 010 864 Email: mannxdoky@gmail.com 75 SCD2021 International Conference on sustainable construction development in the context of climate change in the Mekong Delta Nguyễn Duyên Phong Học hàm, học vị: TS Khoa Xây dựng, Trường Đại học Mỏ - Địa chất Tel: 0967 318 556 Email: nguyenduyenphong@humg.edu.vn ĐẶT VẤN ĐỀ Để xây dựng công trình khai thác nước ngầm với việc bổ cập nước từ nguồn khác cần phải có điều kiện sau: có nguồn nước tự nhiên nhân tạo bên cạnh (sơng, hồ), có cấu trúc địa chất thủy văn quan hệ thủy lực nước mặt nước đất thuận lợi Tồn kiểu bổ cập quan hệ thủy lực nước mặt nước đất sau: Kiểu 1: Thường xảy đới dư ẩm dọc theo sông lớn Trong điều kiện tự nhiên phần lớn thời gian năm sơng, hồ nước đất cung cấp, dịng chảy ngầm hướng từ bờ phía sơng hồ Chỉ mùa lũ thời kỳ lũ nước đất tạm thời nước sông hồ cung cấp Sự cung cấp xảy đới ven bờ, làm cho dịng chảy ngầm có phương từ phía sơng hồ đới ven bờ Chiều rộng đới phụ thuộc vào: độ lớn sông; biên độ dao động mực nước sơng; tính thấm đất đá chứa nước Khi có cơng trình khai thác ven bờ với mực nước hạ thấp mực nước sơng hồ nước sơng cung cấp cho cơng trình khai thác Kiểu 2: Phổ biến vùng thiếu ẩm vùng có cấu trúc đặc biệt thuận lợi nước mặt quanh năm cung cấp cho nước đất Vào thời kỳ lũ giá trị cung cấp tăng lên Khi có cơng trình khai thác ven bờ cung cấp tăng Kiểu 3: Đặc trương cho trường hợp dao động mực nước sông không lớn Do độ nghiêng thủy lực nước đất đới ven sơng lớn nên khơng có cung cấp sơng Trong thời kỳ lũ, q trình thóat nước đất khơng khơng bị dừng lại mà cịn tăng lên gia tăng lượng cung cấp Các công trình khai thác ven bờ khơng hạ thấp mực nước đất xuống mực nước mặt, khơng nhận cung cấp từ phía sơng 76 Kiểu 4: Là kiểu quan hệ thủy lực nước mặt nước có áp nằm sâu ngăn cách lớp thấm nước yếu Trong đa số trường hợp mực áp lực tầng chứa nước cao mực nước mặt nên nước đất cung cấp cho nước mặt cách thấm xuyên từ lên qua lớp thấm yếu Nếu mực áp lực thấp mực nước mặt xảy tượng ngược lại Khi có cơng trình khai thác ven bờ, mực nước đất hạ thấp xuống mực nước mặt nước mặt cung cấp cho cơng trình khai thác cách thấm xuyên qua lớp thấm yếu nên lượng bổ sung không lớn Như việc xây dựng cơng trình khai thác nước thấm lọc ven bờ có hiệu từ nguồn nước mặt có quan hệ thủy lực với nước đất kiểu kiểu Đối với kiểu quan hệ phải đắp đập để nâng cao mực nước mặt Còn kiểu quan hệ khoan thêm lỗ khoan dẫn nước từ xuống Dưới trình bày sở lý thuyết để giải toán đặt GIẢI QUYẾT VẤN ĐỀ Đại lượng thấm từ nguồn nước mặt xác định công thức thủy động lực (Konopdianxep, E.H Iaxepva, 1983) sau: Q = K.M.B.I, (1) Trong (1): Q - lưu lượng thấm từ sông hồ, m /ngày; K - hệ số thấm đất đá tầng chứa nước lớp bùn sét lắng đọng đáy sông hồ, m/ngày; M - chiều dày tầng chứa nước, m; B - chiều dài tường bờ, m; I - độ nghiêng thủy lực, xác định công thức: (2) Trong (2): H1 - độ cao mực nước mặt, m; H2 - độ cao mực nước đất cơng trình khai thác, m; L - khoảng cách từ nguồn nước mặt đến Hội thảo Khoa học Quốc tế Phát triển Xây dựng bền vững điều kiện Biến đổi khí hậu khu vực đồng Sơng Cửu Long cơng trình khai thác, m Trong điều kiện tự nhiên cụ thể đại lượng K, M, B khơng thay đổi, để tăng lưu lượng thấm Q người ta phải làm tăng độ nghiêng thủy lực I cách: Cách 1: chuyển cơng trình khai thác gần sơng hồ tốt, chí đặt lịng sơng hồ (Hình 1) Hình Tăng độ nghiêng thủy lực dịng thấm từ sơng hồ cách chuyển cơng trình khai thác phía gần sơng Cách 2: hạ thấp sâu mực nước đất cơng trình khai thác, nhiên hạ thấp đến giới hạn cho phép (Hình 2) Hình Tăng độ nghiêng thủy lực dịng thấm từ sơng hồ hạ thấp mực nước đất cơng trình khai thác Cơng thức tổng qt để tính trữ lượng khai thác nước đất có dạng sau: (3) Trong (3): Qkt - trữ lượng khai thác tiềm năng, m /ng; Qtn - trữ lượng động tự nhiên, m3/ng; Vdh - trữ lượng tĩnh đàn hồi, m3; Vtl - trữ lượng tĩnh trọng lực, m3; α - hệ số xâm phạm vào trữ lượng tĩnh trọng lực cho phép (lấy 30% tầng chứa nước không áp); t - thời gian khai thác, thường hạn định 27 năm (104 ngày); Qct trữ lượng theo, m3/ng Các đại lượng Qtn, Vdh, Vtl hình thành điều kiện tự nhiên, trữ lượng theo (Qct) hình thành trình khai thác, mà bơm hút, mực nước đất bị hạ thấp lơi dịng chảy từ phía vào tầng chắn nước khai thác (thấm xuyên từ khối nước mặt xuống, thấm nghiêng từ tầng chứa nước kề cận vào ) Trong lượng thấm xuyên từ khối nước mặt số điều kiện chiếm tỉ trọng lớn Đó tiền đề để xây dựng giếng thấm kích thích (induced infiltration wells) ven bờ Theo số liệu Sở Tài nguyên Mơi trường, TP Hồ Chí Minh, 2005 (Tóm tắt đề án: Xây dựng mạng quan trắc lún đất khai thác nước ngầm vùng phía Nam TP Hồ Chí Minh) nước đất khu vực TP Hồ Chí Minh bổ cập từ ba nguồn chính: nước Kinh Đơng, nước mưa nước sơng Sài Gịn Sơng Sài Gịn, Kinh Đông, sông Đồng Nai, nước mưa, nguồn bổ cập cho tầng chứa nước đất Pleistocen, Pliocen Nước mưa nguồn bổ cập lớn cho tầng chứa nước Pleistocen gián tiếp cho Pliocen Chính vậy, mực nước ngầm dao động lớn theo mực nước dịng sơng, Kinh Đơng theo hai mùa: mùa khô mùa mưa Nước Kinh Đông cung cấp cho nước đất tầng Pleitocen Khi nước Kinh Đơng dâng cao mực nước đất vị trí cách kênh 20 ÷ 100 m dâng cao Khi xây dựng xong Kênh Đôi mực nước đất dâng lên ÷ 2,5 m so với chưa xây dựng Kênh Đôi Nước mưa nguồn bổ trợ nước đất lớn Khảo sát mực nước đất cho thấy, có trận mưa kéo dài ÷ mực nước đất Củ Chi tăng lên 0,1 ÷ 0,25 m Về mùa mưa mực nước đất cao mùa khô từ 1,5 ÷ 3,0 m Nước sơng Sài Gịn có quan hệ thủy lực với nước đất 77 SCD2021 SCD2021 International Conference on sustainable construction development in the context of climate change in the Mekong Delta Khi nước sông dâng cao bổ cập cho nước đất, ngược lại nước sơng cạn nước đất lại cung cấp cho sơng Theo tính tốn lượng nước bổ cập cho nước đất tầng Pleistocen từ nguồn bổ cập sau: + Từ Kinh Đơng là: QKD = q.L.b = 0,55 285000 1,0 = 156750 m /ngày Trong q lưu lượng đơn vị, m3/ngày (Theo Sở Tài ngun Mơi trường Tp HCM kinh Đông cung cấp cho nước đất với giá trị khoảng q = 0,55 m3/ngày); L - chiều dài phần kênh chưa bị xi măng hóa (L = 400 - 115 = 285 km = 285000 m; b - bề rộng tính tốn, lấy b = 1,0 m) + Từ sơng Sài Gịn là: QsSG = K.L.W.(Fsong - FAquifer)/M Với K - hệ số thấm lớp đất đá đáy sông, K = 1,0 m/ngày; L - chiều dài đoạn sông cung cấp nước đất, L = 22500 m; W - chiều rộng trung bình sơng, W = 80 m; M - chiều dày lớp trầm tích đáy sông, M = 4,0 m; Fsong - mực nước sơng, lấy trung bình theo trạm quan trắc Bình Dương 0,15 m; FAquifer - mực nước ngầm bên sơng lấy trung bình nhiều năm trạm quan trắc Q002 Bình Mỹ 0,0 m Thay số: QsSG = 67500 m3/ngày + Từ nước mưa là: Qm = F.W, m3/ngày F - diện tích hứng mưa khu vực TP Hồ Chí Minh; W - cường độ cung cấp nước mưa cho tầng Thay số Qm = 309532 m3/ngày Đề xuất xây dựng bãi giếng: Quy hoạch vị trí xây dựng nhà máy khai thác nên ưu tiên vị trí ven sơng, kênh, rạch hồ chứa nguồn cung cấp, bổ trợ lớn cho tầng chứa nước khai thác Giảm lưu lượng khai thác nước ngầm việc khai thác, xử lý nguồn nước mặt Trong vùng TP Hồ Chí Minh tồn hai tầng chứa nước lỗ hổng trầm tích Đệ Tứ bở rời có giá trị khai thác tốt: Tầng chứa nước Pleistocen Pliocen (trên, dưới) Hai tầng bổ cập nước mặt từ kênh Đơi, sơng Sài Gịn, sông 78 Đồng Nai, hệ thống sông, rạch khác nước mưa theo mùa mưa Kết quan trắc động thái nước đất xác định động thái thủy văn, tức vùng động thái nước đất bị chi phối chế độ thủy văn sông Sài Gịn, Kênh Đơi mùa mưa, mùa khơ Trên sở nghiên cứu luận khoa học đặc điểm địa chất thủy văn TP Hồ Chí Minh, đề xuất xây dựng bãi giếng khai thác nước thấm lọc dọc theo dải ven bờ sơng Sài Gịn, sơng Đồng Nai, sông khác kinh khu vực thành phố, Kinh Đông Bãi giếng đưa gần sơng độ nghiêng thủy lực lớn nhận lưu lượng thấm từ sông nhiều Các kết tính tốn cho thấy với diện tích đặt bãi giếng sát mép nước sông lưu lượng khai thác tăng 2,0 ÷ 3,0 lần so với bãi giếng cách xa mép nước sông (khoảng cách xa từ 350 ÷ 500 m) Tóm lại, cơng trình khai thác nước đặt giếng ven bờ dạng khai thác có nguồn bổ sung nhân tạo cần áp dụng rộng rãi thời gian tới TP Hồ Chí Minh để làm việc vùng ven bờ sơng nên quy hoạch thành cơng viên, khu du lịch, giải trí,… xây dựng cơng trình khai thác nước đất THIẾT LẬP BÀI TOÁN Như biết điều kiện địa chất thủy văn (ĐCTV) số loại mỏ nước đất không đơn giản bất đồng tính thấm, chứa nước Sơ đồ bố trí lỗ khoan khai thác trường hợp thường khơng có dạng hình học đặc trưng (đường thẳng, đường trịn, lưới, v.v ) mà có dạng phân bố dạng diện tích Vấn đề đặt nên khai thác từ lỗ khoan với lưu lượng để trị số hạ thấp mực nước nhỏ trị số hạ thấp mực nước cho phép mà tổng chúng đạt giá trị cực đại Đây nội dung tốn tối ưu lưu lượng khai thác lỗ khoan Thiết lập toán tối ưu lưu lượng nước lỗ khoan vị trí chúng biết Nghiên cứu tốn tối - ưu ĐCTV nói chung đặc biệt Hội thảo Khoa học Quốc tế Phát triển Xây dựng bền vững điều kiện Biến đổi khí hậu khu vực đồng Sông Cửu Long tốn tối - ưu hóa lưu lượng lỗ khoan cơng trình khai thác ấn định vị trí chúng J.K Gavitch F.M Botrever đề cập Nội dung toán xác định lưu lượng lỗ khoan khai thác nước Qi để cho tổng lưu lượng khai thác QT thời gian khai thác t đạt cực đại, mà trị số hạ thấp mực nước SJ lỗ khoan j không vượt trị số hạ thấp mực nước cho phép [S]J lưu lượng lỗ khoan thay đổi từ đến Qmax (giá trị Qmax phụ thuộc vào hệ số dẫn nước tầng chứa nước vị trí đặt lỗ khoan, cấu trúc đặc tính kỹ thuật máy bơm) Để đơn giản hóa, giả thiết chất lượng nước đảm bảo khơng phải đề cập đến giới hạn chất lượng nước (độ khóang hóa, thành phần hóa học, v.v ) Với điều kiện toán thiết lập sau: xác định lưu lượng lỗ khoan khai thác Qi(i = 1,2,3, , k), tổng chúng đạt cực đại: (4) Đồng thời tuân theo giới hạn nói trên: Sj(t) ≤ [S]j (5) ≤ Qi ≤ Qmax (6) Trong (4), (5) (6): F - hàm mục tiêu; i, j = 1, 2, , k số thứ tự số lỗ khoan khai thác nước Giá trị Sj phụ thuộc vào điều kiện ĐCTV vùng, lưu lượng lỗ khoan Qi thời gian khai thác t Đối với tầng chứa nước vô hạn giá trị Sj xác định theo công thức (4) hay (5) Đối với nước có áp mối quan hệ (5) thường mối quan hệ tuyến tính trị số hạ thấp Sj lưu lượng lỗ khoan đến tính thấm, chứa nước, loại biên, khoảng cách từ biên đến chúng khoảng cách lỗ khoan Về mặt thủy lực sức cản thủy lực nước đến lỗ khoan, xác định theo công thức: aji = fji/2T (9) Giá trị aj,i mặt vật lý trị số hạ thấp mực nước lỗ khoan thứ j ảnh hưởng khai thác lỗ khoan thứ i với lưu lượng đơn vị Biết aj,i cho phép xác định Sj(t) tất lỗ khoan tương tác Giới hạn (5), (6) biểu diễn dạng bất phương trình biểu thức Sj(t) viết hệ phương trình đại số tuyến tính Một số tác Botrever Ph.M.,(1978); I.K Gavitch, (1988) đề xuất cách xác định Si theo công thức sau đây: (10) Trong (10): k - hệ số thấm; m - chiều dày lớp đất đá, m; rij - khoảng cách từ lỗ khoan thứ i đến lỗ khoan thứ j, i = j rij bán kính ống lọc lỗ khoan Từ (5), (6) (8) nhận phương trình điều kiện giới hạn: F = Q1 + Q2 + Q3 + + Qk → Max (11) Các ràng buộc: (12) (13) (7) Trong đó: aji -Các hàm ảnh hưởng, xác định sau: (8) Trong (7) (8): fj,i - sức cản thủy lực không thứ nguyên; T - hệ số dẫn nước, m2/ngày; aj,i hàm ảnh hưởng lỗ khoan thứ i đến j có tính - Từ k lỗ khoan, thí dụ lỗ khoan thứ (i = 1) đặt lưu lượng nhất, lỗ khoan cịn lại xem khơng; - Xác định giá trị Sj cho tất lỗ khoan tương ứng với aji; 79 SCD2021 SCD2021 International Conference on sustainable construction development in the context of climate change in the Mekong Delta - Bằng phương pháp tương tự tính cho tất lỗ khoan cịn lại (i = 2, i = 3, …, i = k) tìm giá trị cịn lại hàm ảnh hưởng aj,i nhận Đối với tầng chứa nước có áp, vơ hạn, khai thác kéo dài, hàm ảnh hưởng aji tính theo cơng thức: (14) đây: T - hệ số dẫn nước, m2/ng; a - hệ số truyền áp, m2/ng; t - thời gian bơm khai thác, ngày-đêm; rji - khoảng cách từ lỗ khoan j đến lỗ khoan lại (m), j = i rji = rjj = đường kính ống lọc lỗ khoan hút nước (thường 0,2 ÷ 0,4 m) PHƯƠNG PHÁP GIẢI Bài toán thiết lập đặc trưng hàm mục tiêu F tất giới hạn tuyến tính với thơng số tối ưu Qi trình bầy tốn quy hoạch tuyến tính Nếu hàm mục tiêu hay giới hạn phi tuyến so với Q hệ phương trình (7), (8) toán quy hoạch phi tuyến Trong số trường hợp hàm phi tuyến thay hàm tuyến tính khoảng quy hoạch phi tuyến thay quy hoạch tuyến tính Giải tốn nhiều cách Một phương pháp giải dùng phương pháp simplex Để làm việc ta viết lại toán dạng: (15) Các ràng buộc: Lập bảng Simplex có dạng sau đây: Bảng Bảng sở t X1 X2 X3 Xm Xk G Q1 Q2 Q3 a11 a12 a13 a21 a22 a23 a31 a32 a33 t am1 am2 am3 ak1 ak2 ak3 C1=- C2=- C3=1 1 Qs a1S a2S a3S amS akS CS Qk B=[S] a1k B1 a2k B2 a3k B3 amm Bm akm Bk Ck=1 Biến đổi bảng sở theo nguyên tắc sau đây: - Chọn cột cho phép: cột chứa giá trị tuyệt đối Cs lớn Giả thiết cột s (cột tơ đậm) - Sau có cột sở s ta tính giá trị cột theo công thức: ajS = Bi/aiS ; j =1, 2, , k - Chọn hàng cho phép: hàng chứa giá trị nhỏ từ ajS tính Giả sử hàng m (hàng tô đậm) Biến đổi bảng sở theo bước sau: + tính yếu tố cho phép: amS thay a’mS = 1/amS; + yếu tố hàng cho phép: a’mi = ami/amS; + yếu tố cột cho phép: a’jS = - ajS /amS; + yếu tố lại bao gồm hàng tự B hàng G tính sau: + Thay đổi vị trí biến số QS cột s cho biến Xm hàng m Biến đổi bước bước tất yếu tố hàng G không âm (ứng với Xi) Khi giá trị Bki cột cuối ứng với Xi ≥ nghiệm tốn tối ưu (16) Bảng Simplex cuối ứng với giai đoạn biến đổi thứ k yếu tố hàng G Cki khơng âm có dạng: Bảng Bảng Simplex giai đoạn biến đổi thứ k Điều kiện: (17) 80 Q/X Q1 Q2 X1 ak11 ak21 X2 ak12 ak22 X3 ak13 ak23 Xk-1 Xk ak1k ak2k B=[S] Bk1 Bk2 Hội thảo Khoa học Quốc tế Phát triển Xây dựng bền vững điều kiện Biến đổi khí hậu khu vực đồng Sơng Cửu Long Q3 Qk G=-F ak31 akk1 Ck1 ak32 akk2 C k2 ak33 akk3 Ck3 Ckk-1 ak3k akkk Ckk Bk3 Bkk Nghiệm tối ưu tương ứng là: Q1 Q2 Qkt Dựa vào (14) số liệu cho dễ dàng xác định giá trị hàm ảnh hưởng (aji, i = 1, 2, 3; j = 1, 2, 3) ứng với lưu lượng đơn vị (lấy 1000 m3/ng) Bài tốn tối ưu hóa biểu diễn hệ phương trình sau: Bk1 Bk2 Bkk Trong trường hợp tốn nhiều biến việc giải lâu nhiều thời gian Chính mà người ta giải theo số thuật tốn lập trình sẵn Sau ứng dụng chương trình tính có sẵn để tìm nghiệm tối ưu Bài tốn thiết lập theo mơ hình chuẩn tốn quy hoạch tuyến tính sau: (18) X1, X2 X3 - biến số phụ không âm đưa vào theo phương pháp Simplex để bất phương trình dạng (12) trở thành phương trình dạng (16) Lập bảng sở (Bảng 3): Bảng Bảng sở tính tốn minh họa số X/Q X1 X2 X3 G = -F Q1 10,07 6,845 4,578 -1,0 Q2 6,845 10,07 4,414 -1,0 Q3 4,578 4,414 10,07 -1,0 B 50 50 50 Biến đổi theo bảng Simplex giai đoạn thứ cho ta kết Bảng 4: Bài tốn cịn giải Maple 9.5 Bảng Kết biến đổi bảng Simplex giai đoạn MINH HỌA LỜI GIẢI QUA VÍ DỤ CỤ THỂ 5.1 Giải phương pháp Simplex Sau nghiên cứu ví dụ: để cung cấp nước cho xí nghiệp người ta khoan lỗ khoan khai thác có đường kính ống lọc: r11 = r22 = r33 = 0,2 m Tầng chứa nước áp lực, vơ hạn có bề dày trung bình 30 m, áp lực mái tầng chứa nước 50 m Từ tài liệu hút nước thí nghiệm xác định hệ số dẫn nước T = 220 m2/ng, hệ số truyền áp a = 2,25.106 m2/ng Các lỗ khoan bố trí theo tam giác với khoảng cách sau: r12 = r21 = 300 m; r13 = r31 = 400 m r23 = r32 = 500 m Thời gian dự kiến khai thác t = 104 ngày đêm (tương đương 27 năm) Theo tài liệu thăm dò, trị số hạ thấp mực nước cho phép lỗ khoan lấy bằng: [S]1 = [S]2 = [S]3 = 50 m Q/X Q1 Q2 Q3 G=-F X1 0,036 -0,115 -0,440 0,039 X2 -0,118 0,191 -0,025 0,047 X3 -0,037 -0,031 0,129 0,027 B 2,046 2,204 3,129 -7,397 Lời giải nhận được: Q1 Q2 Q3 F=-G 2,046 2,204 3,129 +7,397 F = -G = +7379 m3/ng Như vậy, phương pháp Simplex cho ta giá trị lưu lượng khai thác tối ưu cho lỗ khoan: Q1 = 2046 m3/ng, Q2 = 2204 m3/ng Q3 = 3129m3/ng QT = Q1 + Q2 + Q3 = 7379 m3/ng = max 81 SCD2021 SCD2021 International Conference on sustainable construction development in the context of climate change in the Mekong Delta 5.2 Giải phần mềm Maple 9.5 KẾT LUẬN Dưới trích xuất kết giải theo Maple 9.5: Khai thác nước ngầm tất yếu làm cho mực nước ngầm tầng không áp giảm áp lực tầng chứa nước có áp giảm > restart: with(simplex); rb1:=10.07*Q1+6.845*Q2+4.578*Q3