Đang tải... (xem toàn văn)
đánh giá ảnh hưởng của rừng đến sự hình thành dòng chảy lũ trên một số lưu vực miền trung bằng mô hình sóng động học một chiều
ĐÁNH GIÁ ẢNH HƯỞNG CỦA RỪNG ĐẾN SỰ HÌNH THÀNH DÒNG CHẢY LŨ TRÊNMỘT SỐ LƯU VỰC MIỀN TRUNG BẰNG MÔ HÌNH SÓNG ĐỘNG HỌC MỘT CHIỀUNguyễn Thanh SơnTrường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà NộiASSESMENT OF FOREST IMPACTS ON FLOOD FLOW IN SOME RIVER BASINS INCENTRAL VIETNAM USING 1D KINEMATIC WAVE MODELNguyen Thanh SonCollege of Sciense, VNUTóm tắt. Việc mô phỏng lũ thường được giải quyết thông qua các bài toán về quá trình thấmvà quá trình tập trung nước trên lưu vực. Sử dụng mô hình sóng động học một chiều (KW-1D) phươngpháp phần tử hữu hạn kết hợp với phương pháp SCS (Soil Conservation Service) đã thu được nhữngthành công nhất định khi mô phỏng lũ trên một số lưu vực sông ngòi Miền Trung. Bài báo giới thiệukết quả áp dụng mô hình này để đánh giá tác động của các kịch bản thay đổi lớp phủ rừng đến sự hìnhthành đỉnh lũ, từ đó đề xuất việc khai thác rừng hợp lý và sử dụng mô hình như là một công cụ hỗ trợra quyết định cho các nhà quy hoạch lưu vực sông.1. Đặt vấn đềTheo [1], mô hình sóng động học một chiều dựa trên cơ sở xấp xỉ chi tiết không gian lưu vựcvà tích phân số trị các phương trình đạo hàm riêng mô tả các quá trình vật lý diễn ra trên lưu vực nhằmdiễn toán quá trình hình thành dòng chảy sông qua hai giai đoạn: dòng chảy trên sườn dốc và tronglòng dẫn. Mô hình cho phép đánh giá được tác động quy mô nhỏ của các yếu tố mặt đệm trên lưu vựcđến dòng chảy, mở ra một giai đoạn mới trong mô hình hoá các quá trình thuỷ văn.Dựa trên mô hình của Ross B.B và nnk, (Đại học Quốc gia Blacksburg, Mỹ) [6] dùng để đánhgiá ảnh hưởng của việc sử dụng đất đến quá trình lũ với mưa vượt thấm là đầu vào của mô hình,phương pháp phần tử hữu hạn kết hợp với phương pháp số dư của Galerkin được sử dụng để giải hệphương trình sóng động học của dòng chảy một chiều.Phương trình liên tục: (1)Phương trình động lượng (2)trong đó: Q: Lưu lượng trên bãi dòng chảy trên mặt hoặc trong kênh. q: Dòng chảy bổ sung ngang trênmột đơn vị chiều dài của bãi dòng chảy (mưa vượt thấm đối với bãi dòng chảy trên mặt và đầu ra củadòng chảy trên mặt đối với kênh dẫn). A: Diện tích dòng chảy trong bãi dòng chảy trên mặt hoặc trongkênh, x: khoảng cách theo hướng dòng chảy, t: thời gian, g: gia tốc trọng trường, S: độ dốc đáy của bãidòng chảy. Sf: độ dốc ma sát. y: độ sâu dòng chảyThuật giải hệ phương trình trên đã được trình bày trong [1], áp dụng cho lưu vực sông TràKhúc [5], Vệ [2],và Tả Trạch [4] kết hợp với phương pháp SCS để tính mưa hiệu quả làm đầu vào củamô hình đã cho kết quả mô phỏng dòng chảy khá tốt.Do tính chất của phương pháp SCS, khi thay đổi các điều kiện mặt đệm trên các lưu vực, sẽảnh hưởng tới việc hình thành dòng chảy mặt, nói chung và hình thành lũ, nói riêng. Việc thay đổi điềukiện sử dụng đất trên lưu vực dẫn tới sự thay đổi hệ số CN [7] và kéo theo là tính chất lũ. Bằng môhình KW -1D có thể khảo sát được mức độ ảnh hưởng của việc khai thác lớp phủ rừng đến sự hìnhthành lũ, từ đó có thể tư vấn cho các nhà quy hoạch về giới hạn nguy hiểm của mức độ sử dụng đất đốivới sự phát triển bền vững tài nguyên đất và nước.Qua khảo sát các lưu vực nghiên cứu, phân tích các điều kiện thuận lợi về địa hình, giao thông và các điều kiện dân cư, tác giả đã xây dựng các kịch bản sử dụng lớp phủ rừng, rút ra những kết luậnvề mức độ giới hạn sử dụng đất tác động đến sự hình thành lũ để các nhà quy hoạch có định hướngđúng, tránh gây các hiểm họa về thiên tai lũ lụt do sự quá tải khi khai thác lưu vực. Yếu tố sử dụng đấtđược xét đến là quá trình thay đổi lớp phủ rừng. Các bước phân tích và tính toán được tiến hành theocác kịch bản như sau:1. Phân tích các điều kiện địa lý tự nhiên và kinh tế xã hội trên các lưu vực sông nghiêncứu để xây dựng kịch bản sử dụng đất2. Thay đổi dần diện tích sử dụng đất trên các phần tử, lựa chọn lại hệ số CN và tính toánlại quá trình lũ bằng mô hình KW-1D3. Xây dựng quan hệ giữa diện tích sử dụng đất (%) theo các kịch bản với đỉnh và tổnglượng lũ4. Đưa ra kết luận về ảnh hưởng của quá trình sử dụng đất đến sự hình thành dòng chảylũ trên lưu vực.2. Số liệu Số liệu dùng để mô phỏng lũ và khảo sát ảnh hưởng của lớp phủ rừng là tài liệu mưa giờ vàtrích lũ từ năm 1998 - 2004 cùng với bản đồ số về lớp phủ rừng năm 2000, tỷ lệ 1: 50.000Số trận lũ được chọn để khảo sát như sau: Lưu vực Tả Trạch – Thượng Nhật: 9 trận; TràKhúc – Sơn Giang: 3 trận và Vệ – An Chỉ: 2 trận3. Đánh giá ảnh hưởng của lớp phủ rừng đến sự hình thành lũKhai thác rừng đang là một vấn đề xã hội cần quan tâm khi mà vai trò của rừng ngày càngđược khẳng định. Hiện nay về lý thuyết có thể hạn chế những nhu cầu về lâm sản bằng việc thay thếnhững vật liệu mới, nhưng thực tế không triệt tiêu được nhu cầu đó. Hơn nữa chính những nhu cầu ấycòn đem lại không ít việc làm cho người dân và vì vai trò của rừng rất quan trọng trong nền kinh tếquốc dân nên cần xác định mức độ khai thác rừng đảm bảo không vượt quá mức giới hạn, không làmtăng đột ngột quá trình lũ, gây ra những hậu quả nghiêm trọng đồng thời lại tận dụng được những tiềmnăng kinh tế của rừng.Qkhaithac(m3/s)%547 21.6544.9 25.3521 31.5493.8 35.8482.7 38.5477.9 41.8Qtrồng(m3/s)%477 44.2474.4 47.9472.5 50.6442.7 54.7434.6 58.6422.8 62.3410.5 64.9Wkhaithac(m3)%21333 21.621251 25.320319 31.519258 35.818825 38.518638 41.8Wtrông(m3)%18603 44.218501 47.918427 50.617265 54.716949 58.616489 62.316009.5 64.9 Hỡnh 1. nh hng ca lp ph rng n dũng chy l trờn sụng T Trch - Thng Nht trnl t 16h/21n 4h/23/X/2000Thc hin vic kho sỏt nh hng ca quỏ trỡnh khai thỏc rng trờn lu vc sụng T Trch -trm Thng Nht nhm xỏc nh ngng khai thỏc rng m khụng gõy t bin i vi quỏ trỡnh lbng kch bn ó nờu trờn. Tin hnh tng din tớch khai thỏc rng trờn lu vc dn ti vic xỏc lpli b thụng s ca mụ hỡnh KW 1D, c th l ch s CN v n cho cỏc phn t khi cú s thay i sdng t [7] (vi CN = 79, n = 0.3) . Sau ú tin hnh thay i cỏc thụng s CN v n trong cỏc filethụng s ca mụ hỡnh KW 1D tin hnh thc nghim s vi cỏc trn l ó la chn trờn lu vcsụng T Trch Thng Nht. S thay i din tớch rng trờn lu vc sụng T Trch vi kch bn t tng din tớch rnghin trng (khụng bn n chng loi rng) chim 44.3% din tớch lu vc, gim che ph theo cỏcmc tun t l 41.8% đ 38.5% đ 35.8% đ 25.3% đ 21.6% din tớch lu vc.Sau khi thay i CN v h s nhỏm trung bỡnh phn t trong cỏc file thụng s tin hnh tớnhtoỏn li cho 9 trn l ó chn trờn sụng T Trch. Th hin kt qu ca cỏc kch bn khai thỏc rng trờn lu vc sụng T Trch n dũng chy lc kho sỏt bng cỏc quan h Qmax~Frng, W~ Frng (Hỡnh 1)ỏnh giỏ nh hng ca vic khai thỏc rng n dũng chy l qua 9 trn l trờn sụng TTrch - trm Thng Nht, thy rng: khi din tớch khai thỏc rng tng thỡ nh v lng l tng lờn dos gim nhỏm v kh nng thm nc ca lu vc dn ti s tng dũng chy mt v nh l. Vicgim din tớch rng c th hin trong hỡnh 1 cho thy rng cú th khai thỏc rng nhng vn phi gidin tớch rng trờn lu vc vi che ph khong 35% din tớch lu vc thỡ khụng gõy t bin v l.Vic khai thỏc rng trờn lu vc sụng T Trch lm cho che ph di gii hn 35 % s lm cho lng dũng chy v nh l tng t ngt, gõy nguy him cho s phỏt trin kinh t xó hi.Vỡ vy, kin ngh vi cỏc nh quy hoch s dng t rng, khai thỏc rng hiu qu phi kthp vi vic tỏi sinh rng m bo hin trng rng khu vc n nh khụng cú hin tng nhl tng t bin gõy nguy him. Vic trng rng úng vai trũ ht sc quan trng trong quỏ trỡnh hỡnhthnh v iu tit dũng chy nht l i vi cỏc khu rng u ngun v rng phũng h. Vic tng dintớch rng s gúp phn iu hũa dũng chy trong nm, nú lm gim dũng chy mt v tng dũng chyngm, t ú dn n lm gim tng lng dũng mt, nh l. Ngoi ra rng cũn gúp phn vo vicchng xúi mũn v ra trụi trờn b mt lu vc, lm tng phỡ nhiờu ca t v lm cõn bng h sinhthỏi trờn lu vc. Ngoi lu vc sụng T Trch trm Thng Nht, ó tin hnh cỏc kho sỏt tngt i vi cỏc lu vc sụng Tr Khỳc trm Sn Giang, sụng V trm An Ch. Kt qu tng hpngng khai thỏc rng n nh l c trỡnh by bng 1.Bng 1. Kt qu kho sỏt ỏnh giỏ nh hng ca rng n dũng chy l trờn cỏc lu vc sụngnghiờn cuTT Lu vc nghiờn cuS trn l khosỏt thc nghimsNgng che ph rng gii hn (% din tớch lu vc)1 T Trch Thng Nht 9 352 Tr Khỳc Sn Giang 3 303 V An Ch 2 404. Kt lun1. Vi vic xp x chi tit khụng gian lu vc v tớch phõn cỏc phng trỡnh o hm riờng mụ tcỏc quỏ trỡnh vt lý din ra trờn lu vc, mụ hỡnh KW 1D cú kh nng ỏnh giỏ c nhngthay i trong phm vi nhng khụng gian nh trờn lu vc n quỏ trỡnh hỡnh thnh dũngchy. Tớnh bin ng theo khụng gian ca hỡnh dng lu vc, ca cỏc c tớnh thu vn vma cú th d dng c xột n trong mụ hỡnh trờn.2. Vi s liu u vo l ma vt thm v cỏc bn s v t, rng, a hỡnh vic ỏp dng mụhỡnh cú tớnh kh thi cao khi ỏnh giỏ tỏc ng s thay i ca cỏc yu t t nhiờn ti dũng chảy. Một sự biến động nào đó trên một phần tử sẽ có tác động đến toàn bộ hệ thống và ảnhhưởng đến dòng chảy trên sông. Có thể dùng phương pháp này để đánh giá các quy hoạch sửdụng đất đối với việc đảm bảo bền vững tài nguyên nước.3. Kết quả khảo sát các kịch bản sử dụng lớp phủ rừng cho thấy để đảm bảo tính chất tự nhiêncủa sự hình thành lũ, không nên khai thác rừng tới ngưỡng tới hạn độ che phủ của từng lưuvực. 4. Các kết luận trên được xây dựng qua việc khảo sát một số trận lũ, mang ý nghĩa về phươngpháp, chứng minh khả năng của mô hình KW - 1D và phương pháp SCS đối với việc đánh giácác điều kiện mặt đệm lên dòng chảy [3]. Tính thực tiễn và độ chính xác của nó cần đượckiểm chứng với số lượng các trận lũ lớn hơn. Ngoài ra, cần xây dựng thêm các kịch bản phốihợp giữa các yếu tố mặt đệm khác như lớp phủ thực vật, độ dốc, hướng sườn … để có các kếtluận chính xác và chi tiết hơn, sẽ được bàn tới trong các công bố tiếp theo.Các kết quả nghiên cứu công bố trong bài báo này được sự hỗ trợ kinh phí của Đề tài 705606thuộc Chương trình NCCB giai đoạn 2006 – 2008 và Đề tài QG. 07.15. Tài liệu tham khảo1. Nguyễn Thanh Sơn, Lương Tuấn Anh (2003), "Áp dụng mô hình thuỷ động học các phần tử hữuhạn mô tả quá trình dòng chảy lưu vực", Tạp chí khoa học Đại học Quốc gia Hà Nội. Khoa họcTự nhiên và Công nghệ, T.XIX, No1-2003, tr 90-992. Nguyễn Thanh Sơn, Ngô Chí Tuấn (2004) "Kết quả mô phỏng lũ bằng mô hình sóng động họcmột chiều lưu vực sông Vệ", Tạp chí khoa học Đại học Quốc gia Hà Nội. Khoa học Tự nhiên vàCông nghệ,T.XX, No3 PT-2004 tr. 44-50.3. Nguyễn Thanh Sơn (2006), "Áp dụng mô hình 1DKWM – FEM & SCS đánh giá tác động của quátrình đô thị hóa đến dòng chảy lũ trên một số sông ngòi Miền Trung" Tạp chí khoa học Đại họcQuốc gia Hà Nội. Khoa học Tự nhiên và Công nghệ, T.XXII, số 2B PT – 2006, tr. 149-157, HàNội4. Nguyễn Thanh Sơn (2006), "Kết quả mô phỏng lũ bằng mô hình 1DKWM – FEM & SCS lưu vựcsông Tả Trạch trạm Thượng Nhật" Tạp chí khoa học Đại học Quốc gia Hà Nội. Khoa học Tựnhiên và Công nghệ, T.XXII, No 3- 20065. Trương Quang Hải, Nguyễn Thanh Sơn (2006), “Kết quả mô phỏng lũ bằng mô hình 1DKWM–FEM&SCS lưu vực sông Trà Khúc – trạm Sơn Giang”. Tạp chí khoa học Đại học Quốc gia HàNội. Khoa học Tự nhiên và Công nghệ, T.XXII, số 4PT – 2006, Hà Nội. tr 59-686. Ross B. B, Contractor D.N and Shanholtz V. O., (1979), "Finite element model of overland andchannel flow for assessing the hydrologic impact of land - use change". Journal of Hydrology,(41), p.11-30.7. Chow YT, Maidment DR, Mays LW. (1988), Applied hydrology. NewYork: McGraw-Hill; . ĐÁNH GIÁ ẢNH HƯỞNG CỦA RỪNG ĐẾN SỰ HÌNH THÀNH DÒNG CHẢY LŨ TRÊNMỘT SỐ LƯU VỰC MIỀN TRUNG BẰNG MÔ HÌNH SÓNG ĐỘNG HỌC MỘT CHIỀUNguyễn Thanh. đất đến sự hình thành dòng chảylũ trên lưu vực. 2. Số liệu Số liệu dùng để mô phỏng lũ và khảo sát ảnh hưởng của lớp phủ rừng là tài liệu mưa giờ vàtrích lũ