Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình

Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình Nghiên cứu ứng dụng mô hình HEC-6 để mô phỏng diễn biến lòng sông Đà phía hạ lưu đập Hòa Bình

M Ụ C L Ụ C

LỜI MỚ Đ Ẩ U 2

Chương 1 Đ I Ể U KI ÊN ĐỊA LÝ T ự N H IÊ N LƯU v ự c S Ô NG HỒN G VÀ ĐẶC Đ IỂ M Đ O Ạ N M Ạ N G S Ô N G N G H I Ê N c ứ u 3

1.1 Vị trí địa l ý 3

1.2 Đ ị a h ì n h 4

1.3 Cấu tạo địa chất và đ ấ t 4

1.4 Lớp phủ thực v ậ t 5

1.5 Khí h ậ u

1.6 M ạng lưới sóng suối, đặc điế m đoạn m ạ n g s ồng nghiên cứu va tình hình tài liệu thúy văn đã thu t h ậ p 7

1.7 Q u y trình vận hành và một số thông số cơ bán của hô H oà b ì n h 1 2 Chương 2 c ơ SỞ LÝ T H U Y Ế T CỦA M Ô HÌ N H H E C - 6 14

2.1 Giới thiệu m ỏ hình H E C - 6 14

2.2 Cơ sớ lý thuyết tính toán thúy l ự c 1?

2.3 Cơ sớ lý thuyết tính toán bùn c á t 20

2.4 Các yêu cầu về số liệu đầu v à o 2C> 2.5 Kiếm soát đẩu r a 3 I Chương 3 ỨNG DUNG MỔ HÌNH HEC-6 TÍNH TOÁ N VẢ D ự BẢO DIKN BIẾN Đ OAN SÔNÍỈ h ổ n g T Ừ H O À b ì n h đ ế n h a n ổ i 33

3.1 Cơ sớ s ố l i ệ u 33

3.2 Sư đồ t í n h 34

3.3 Phân tích, xử lý số liệu đầu v à o 34

3.4 Đ i ểu kiện ban đầu và điều kiện b i ê n 39

3.5 Hi ệu chi nh và k iể m n gh iệ m m ô hình thúy l ực 39

3.6 Tí nh toán, lựa chọn côn g thức sức tái c á t 43

3.7 M ô p h ỏ n g diễn biến l ò ns sông H ồ n e đoan n ghi ên cứu đến nam 1997 ' ĩ 44

3.8 D ư báo th ử diễn biến lòng s ô n a H ổ n e đ o an n s h i é n cứu đến năm 1997 T ’ 46

3.9 Dụ báo d iễn biến lòng sông H ó n g đ o an n g h i ê n cứu đến năm 2052 48 KẾT L l l Ậ N 51

TÀI LIỆU T H A M K H Ả O 54

L Ò I M Ở Đ Ẩ U

Với hai n h iệ m vụ q uan trọng h àn g đầu là phát điện và p h ò n s lũ cho ha du từ khi vận hành, hổ Hoà Bình đã và đang đ em lại cho con người nhii ne lọi ích hết sức lo lớn nhưng c ũ ng đã và đ an g gây ra cho con người nh ững thay đối bát lợi như quá trình hổi lắng lòng hổ, qu á trình xói lớ bờ hố, quá trình xói lớ cục bô 0' ha lưu sau chân đập và q u á trình xói lớ lan truyền vé hạ du Bơi vậy đánh giá m ỏ phónsi và tlư báo các thay đổi sau khi xây d ựne cô ng trình đ ập H o à Bình nói l i e n s và các cõim trình xây tlưng Irẽn sõng nói c hun g nhằin để xuất các biện pháp hạn c h ế các ihicí hại luón lá dê tai có tính cap Ihiết và có ý nghĩa thực tiền lớn C õ n e cụ hữu hiện và kinh

tế nhất đế giai qu yết các đề tài này là mô hình toán

H E C - 6 là m ô hình số thủy đ ộng lực một chiều trong lòng dán hớ có biên di

các thay đổi trong trắc diện dọc sông do xói hoặc bồi t r o n s các thời đoạn vừa

Đê tài này n ghiê n cứu ứng d u n s mò hình H E C - 6 đê’ m õ p h ó n e cliỏn biên lòn Li sõng Đà phía hạ lưu đập Hòa Bình (đoạn từ Hòa Bình đôn niũi ha Th no -Đ à) và lòim sông H ỏ ng (đoan từ ngã ha T h a o - Đ à đến Hà Nội) Báo cáo đê tài đưực bó cục ihành

ba ch ươn g ( k h ôn g kế phấn m ở đáu và kết luận):

Ch ươn g 1 Điêu kiện địa lý tự nhiên lưu vực s ô n e Hỏn li và tlặc đ ié m (.loạn

mạ ng s óng nghiê n cứu

Chư ơn g 2 Cư sớ lý thuyết của mò hình I IEC-6

Chư ơn g 3 line d u n s m õ hình H EC- 6 tính toán và dụ' bán dicn hicn đoan SPI1L!

H ổ n s từ Hò a Bình đến Hà Nội

Lưu vực s ô ng T hao tính đến Việt Trì có diện tích 51 800 krrr trong đó phán ihuộc địa phận Việt Nam chi 12 00 0 krrr N ằ m ớ vị trí c hu y ê n tiếp từ đó n u hắc saiiii táy bác cua Bác Bộ lưu vực sòng T hao kéo dài k h o án a 250 km từ tây bãc xLìônu dóng nam với chiều r ô n e trung hình khoang 50 km Giới hạn phía d o n e là dãy núi Con Voi (đường phán nước giữa hai lưu vực s õng T hao va s õng C ha y ); gioi han phía tây la dãy núi cao Ho àn g Liên Sơn-Pu L uô ng (đường phán nước mữa hai Ill'll vực sỏne T hao và sông Đà) với đinh F an- xi- pan cao nhất nước ta (3143 m ) [3][6].

Lun vực sông Đ à có tổng diên tích là 52 9 00 k n r trontì đi') phan thuộc dĩa phán Việt N am k h o á n g 26 800 k m 2, N ã m sáu trong lục địa thuộc khu táv hãc Bác

Bộ, lưu vực sông Đ à kéo dài tới 380 km theo hướnti tãv b ác- đôn II nam \'ứi c h i c II rónti trung hình k h o á n g 80 km Giới hạn phía đ ô n s là dãy núi HoànII Liẽn Sơn-Pu Lu ó na (đường phân nước giữa hai lưu vực s ông Đ à và s ô n ẹ T h a o ) với các đỉnh cao từ 250 0

m đến trên 3 00 0 m Giới hạn phía tày là các dãy núi cao như: Pu-đen-đinh( 1886 m) Phu-hiiỏi -Lonẹ ( 21 78 m) Phu-la-ina (1801 m) Phu-Tur m ( 1 4 86 nu \'à Phu-Sanu ( i 5 ] s 111) C hu n g cunii chính la đưưng phan nirớc mil'll Ill'll vực sòn e Đ a VỨ1 các lưu vực sòng Mẽ K ô n g và sông Mã Phía cực bác có dãy núi cao Pu- si -lune ( 3076 m ) \'à Ngũ-đài- sơn ( 30 48 m) Phía đông nam là vùn tỉ núi t hấp Ba Vì (1287 m) \ a Vicn

N am ( 1029 m) - đư ờn g phân nước giữa hai lưu vực s ô n e Đ à và s óng Đ á y [5][6]

Lưu vực Sỏns Lô có t o n s diện tích là 39 0 0 0 k n r tronsỉ dó phan nam tròn địa phận Việt Nam k h o an g 22 600 krrr Giói hạn phía đ ổ n s là cánh c u n<2 N e ãn Son \ a canh c un g s ỏ n e G â m , phía d ò n g nam là dãy Ill'll T a m Đ a o \'à phía lãv là dãy núi Con Voi [6],

V ù n u d ồ n g bằng châu thố do phù sa sóng HỔI12 bói đáp nén từ Viội Trì mo' rộng ra tới biển V ù n g đ ồ n g b ằng này rộng tới 1 1 3 00 k n r có d a n s tam giác với đỉnh là Việt Trì [4],

3

Đặc đi ể m nổi bật cu a địa hình trong lưu vực sông Đ à là các d ạ n e địa hình mil

và cao ngu y ên đều cao bị chia cắt m ạ n h theo chi ều thảng đ ứ n s và đ a n s được 11 ũ ne lén mạ nh do chịu ánh hưởng của hoạt đ ộn g tân kiến tạo Các địa hình núi cao nguyên và t hung lũng sắp xếp song song theo hướng tây b á c - đ ỏ n s nam với đo cao

t r u n g h ì n h l ư u v ự c đ ạ t 9 6 5 m [ 5 ] [ 6 ] ,

Lưu vực sổng Lô có dạng địa hình phổ biên là đổi núi thấp Các (lũv núi I'.'n trong Ill'll vực đểu quy tụ về phía đ ông nam và m ớ rộ ng về phía bắc tao cho Ill'Ll vực

có dạng hình nan quạt Đ ộ cao bình qu ân lưu vực đạt 884 m, tr one đó độ cao từ 200

m trờ lén ch iế m tới 7 0- 8 0 % diện tích lưu vực [5][6]

1.3 Cáu tao địa chất và đát

Các vùng núi cao trong lưu vực của hệ thốn 2 s óng Hỏ n tí được cấu lạo hời các loai đá như: granit đá phiên, sa diệp thạch, phiến thạch, sa thạch, cát kct, cuội két va

đá vôi Dãy H o à n g Liên Sơn và các dãy núi O' v ùng biên giói Vi ệt -Lào được call lạo hói các đá có n g u ồ n gốc m a cma Khói núi vòm s óng Ch ay ỏ' thượng i m u m Sony Cháv với đính T â y Cô n Lĩnh cao 2418 m là khối núi granit lớn nhất, diện lích tới

2500 k n r Phía nam khối núi này là vùng đồi núi t hấp ( bao ỉiổm thung lũne snnu Cháy, sòn í Lõ hạ lưu s ông G â m và s ỏ na Phó Đ á y ) được cấu tạo chu yếu hăim ilá phiên, dá cát kèl đá phiên kết tinh và các loại đá biến chát khác Cánh cu nu sonii

G á m được cấu tạo băng đá phiến thạch anh và đá cát kết đá vòi Đá vôi phán bu rộng rãi trên các cao n g u yê n đá vôi Tá phình-Sín Chai Sơn L a - M ộ c C hâ u ironii Ill'll vực sông Đà các cao n g u y ên đ ồn g Vãn Q u á n Bạ trong lưu vực sôniỉ Lo va cánh

eu n c N e â n Sơn, Y ê n L ạ c [7],

P h á n b ă c v à đ ô i ì i ỉ b á c l ư u v ự c t h u ộ c v ù n t i đ á t n ú i p h á t t r i c n t r ô n đ á VOI h i c i n trtV, lì đ ấ t h ã n g , c ó r ừ n a c h e p h ủ , đ ấ t p h á t t r i ế n t r ẽ n d i ệ p t h ạ c h , s a i h a c h \ ' à đ á v ù i

D o v ậ y , l i r ự ni i b ù n c á t hi x ó i m ò n đ ê c u n a c ấ p c h o SỎI12 k h ó n s n h i ê u n ó n sõiiii [,i'i

m a n s rât ít b ù n c á t V ù n e đ ấ t t h u ộ c d ã y n ú i F a n xi p a n - P u l u ỏ n t i c ó d i ẹ n l í c h r ọ n e

độ cao và địa hình thay đổi lớn được phủ bới các loại đất rất k hác nhau Đủi O' vùnt: này được phát triển trên các loại đá gốc như diệp thạch tinh thè hoa ciro’Hi’ Cai i’j

4

xuố ng phía tây n am , diệp thạch và đá vôi càng nhiều còn ớ phía đ ô n g nam la diệp thạch và hoa cương Đ â y là khu vực cun g cấp bùn cát q u a n trọng cho s õ n s Đ à và đón g gó p phấn c hủ yếu vào bùn cát sông Hồng Khu vực bên phải cua s o n s Đa co cao n g u y ê n đá vói k é o dài từ tây bắc x u ố n g đ ỏ n g n am , đát phát triên trên đá vói có

độ mịn lớn Ngoài ra, còn có các loại đất phát triển trên d iệp thạch, sa thạch, hoa cương T h á m thực vật ớ đáy bị phá hủy n gh iê m trọng nén rất thuận lọi cho xói mon[8], Các vùng đ ồ n g báng xen đồi được bồi lấp bới phù sa hiên đại hoặc bới cac vật liệu từ các đổi đưa xuống Vù ng đồng bằng ha lưu s ông H ỏ n s được bói tích bơi phu

sa sóng với lớp trầm tích dày [6],

1.4 Lớp phủ thực vật

Trước nh ữn g năm 80, lớp phú rừng cua lưu vực hệ t h òn e sòng H o n s còn kha

tý lệ phú rừng của một số lưu vực sóng vừa và nho c ủ a sóntỉ Đ à từ 30 đến 5 5 ' ; cua sông T hao từ 50 đến 6 5 % và cua sông Lô từ 45 đến 9 2 % ( b án g 1.1)

B ảng 1.1

Tỷ lệ che phủ rừng của m ộ t sô lưu vực s ô n g vừa và nhỏ tro n g lưu vực hệ th ố n g sô n g H ống năm 1980 [8]

TTThuộc lưu vực son2

N hưi m từ n h ũ n g n ă m siữa thập ky 80 lớp phú thực vật trẽn lưu vực chu yêu la

r ùng dã bị tàn phá nãng T h e o tài liệu cua Viện Đ i êu tra Q u y h oạ ch rừne năm 1990[9], tv lệ rừng p hò n g hộ trên thượng n gu ồn s òn g Thao Đ à và Ló chi còn khoanti lừ

8 r’f đến 2 7 r/c ; ít hơn nhi ều so với tiêu ch uấ n đô che phu rìnm đáu nìiuón p h o n e hộ an loàn q u y định c ua F A O (-10%) T hư ợng n gu ồ n s ó n e T h a o (phán Việt N a m ) là nơi 1\

5

lệ rừng còn ít nhất, chỉ từ 8, 5% đến 15% ( bảng 1.2) Đ â v là một trong nhũ ng nguyên nhân chí nh làm tăng lượng dò ng chảy bùn cát trên thượng n g u ồ n các sõng Đà Thao

Lô trong n hữ ng n ã m từ 1986 đến nay Sự tãng diện tích c ác vùng đất trốn a đổi trọc

đã khiến đất bị xói m ò n n g hi ê m trọng, làm tăng lương bùn cát trong sõng vào m ùa

mư a lũ gây hiện tượng bổi lấp lòng sông ờ hạ lưu, ảnh hư ớn g n g hi ê m irọntỉ đến VÍCL tiéu thoát lũ và gây trớ ngại cho giao thông thủy

Diện tích đất có rừnsi (ha)

Ty lệ rừns so với diện tích đat tư nhicn(rv )

1.5 Khí hậu [7]

Khí hậu trong lưu vực hệ thống sõng Hổn g ớ phán lãnh thổ Viél Nam la khí hậu nhiệt đới gió mùa Các vếu lố khí hậu chịu tác đ ộ n g c úa địa hình, bicn doi mạnh theo k hõ ng gian và thời gian

Bức xa tổng c ộn g trung bình năm đạt k h o á n g 80 k c a l / c n r ớ vùnu núi cao va tăng lên tới 120 k c a l / c n r ờ vùng đ ồn g bằng Cân bằng bức xạ Irung bình nam cìmti hiến đổi từ dưới 4 0 k c a l / c m 2 ớ vùng núi cao đến trên 70 k c a l / c i r r O' vùng đ òn g Mím Can hãng hức xạ trung bình tháng t ư ơ ns đối cao trong các thán LI mù a hè (7-9

k cu l /c m2) và tương đối tháp tron li các thán a m ù a đ ôn g (3-6 k c a l / c n r )

Sỏ giò' n ắ n g t rung binh n ăm biến đổi trong p h ạ m vi từ dưới 1400 ụiờ ơ vìnm núi cao đến hơn 20 00 g iờ ớ các thung lũng

Nhiệt độ k h ô n g k hí trung bình n ăm cũ ng biến đổi m ạ n h theo k h ỏ n a eian va theo mùa T r o ng thời kỳ gió m ù a m ù a hạ nhiệt độ k h ô n g khí t r u n s bình t h á n ” khoánti 15-20"C ớ v ù n a núi, 20-30"C ớ vù na trung dll và đổriiỉ b a n s Troim [hoi k\ LLÌÓ mùa m ùa u ổ n g , nhiệt độ kh ôn g khí t ru n s bình thá ng k h oa i m 10-15’ C ớ \ùnti núi ] 5 - 2 0 ’C ứ \ Cintí, trung du và đ ồ n g bang

Nhiệt độ k h ô n g khí cao nhất tuyệt đối có thế trẽn 40"C’ f(42.fi"C lai Lai Cháu (5/1928) 4 1 r c tại Y e n Châu (13/5/ 1966) 40.7 ’C tại Phú Ho (5/1931 ) va 42.X c tại Láníi ( 5/1926)]

6

Nhiệt độ k h ô n g khí thấp nhất tuyệt đối có thế hạ t hấp dưới 0"C o vùng núi cao [-4,5"C tại Sìn H ổ (3/1/ 1974) , - 4 , 7 ° c tại Cò Nòi (2/1/1974) -5.7"C tại Hoang Llen Sơn ( 14 /12/ 1975) ] và dưới 5"C ớ trung du và đ ồ n g bằng.

Đ ộ ám tương đối củ a khô ng khí trung bình n ă m vào k ho áng 8 0 - 8 5 r f , biên đoi không lớn theo k h ô n g gian trong lưu vực nh ưn g thay đổi m ạ n h theo mùa tương đoi cao trong m ù a m ư a và tương đối thấp trong m ù a khỏ

Lượng m â y tống q u an trung bình n ăm biến đổi trong ph ạm vi 6-cS.5 phan mười bầu trời, tương đối thấp ở lưu vực sồng Đ à và tương đối cao ờ lưu vực sònn Ló

Tố c độ gió trung hình năm biến đổi trong p hạm vi r ộng, từ dưới 1 m/s ớ thnnti lũng và những sườn núi khuất gió đến 3-4 m/s ớ đ ổ n g băng ven biến, thậm chí có the lẽn tứi trẽn 4 m/s o vùng níu cao như ớ H o à n g Liên Sơn

Lương bốc hơi trung bình n ăm (đo b ằn g ố ng Piche) c ũ ng hiến đổi mạ n h tro nu

không gian, có xu t h ế lớn ở vùng trung du và đ ồn g b ằng ( 9 0 0 - 1 0 0 0 n u n ) và nho ớ

vùng núi ( 5 0 0 - 7 00 mill)

Do ánh hướng c ua địa hình nén lư ợn” m u a nãm trung bình nhicu nám ( X ) phân bó rất k h ổ n g đêu trong lưu vực từ 1100-1200 m m ớ những nơi khnấi nió mua

ám đến trẽn 4 0 0 0 m m ớ sườn núi đón gió m ù a ám Nơi m ư a nhiêu nhát co lượim

mu a hàng năm lớn s ấ p 5 lán nơi m ưa ít nhất

Tr on g lưu vực hệ thống sông Hỏng, khu vực Bắc Q u a n a thuộc sườn núi Tày Côn Lĩnh có X , lớn nhất, đạt tới gần 50 00 mm Ngoài ra, còn một sỏ trung tâm inu'a lớn (X ,>3000 m m ) xuất hiên ứ vùng núi cao Hoàn Si Liên Sơn khu vực núi cao ó' hicn giói Việt- Tr ung, ta ngạn s õng Đà T h u n g l ũne Yen Cháu- Mộ c Châu O' hữu nu.; Ill sông Đà và t h ung lũng thượng ng uồ n s ông Ch ảy (Sin M ầ n - N ậ m Vai )la nluìnu míi L'i)

x„ nhỏ nhất (1 100- 1200 mm ) Ngoài ra, t h u n a lũn ti s ỏ n e Mi ên ứ c ao niHiyen Đonii

Vãn trong lưu vực s ôn g G â m cũng có X u tương đói n hỏ ( 1 2 0 0 - 1 3 0 0 m m | , Vimiỉ ha

Ill'Ll thuộc đ ồ n g bằng châu thổ Bắc Bộ có X 0 biến đổi tr ona p hạm VI 1300-2000 mm

tháp nhấl ỏ' khu vực Phúc Yên - Kim A nh thuộc lưu vực s o n s Cà Lò ( 1300-1400

m i n )

C ũn g n hư các yếu tô khí tượng khác, m ưa cũn g biên đổi theo mua Múa mua Irons’ lưu vực hệ t h ốn g s õ n e H ồ n g thường kéo dài từ t h á n ” V đến tháim IX X (chi >

6 t ha ne) n h ưn g có lượng m ư a chi ếm tới 7 5- 8 5 % lượng m ư a ca năm Trotm khi do

m ùa khô kéo dài tới 6-7 tháng nh img có t o n s lương m ư a chí c h i ế m l > 2 V í lượn ti nura cá năm

1.6 Mạ ng lưới s òng suối, đặc điểm đoan ma ng sông nghiên cứu và tình hình tai liệu thúy văn đã thu thập [7]

1.6.1 M a n g ìưóỉ sóng suôi

M a n e lưới sõniỉ suối trong hệ t h ốn e s ô n s H ỏ n e phát t n ế n IvhónL2 đẽ 11 \ ớ i mịii

dò lưứi s òng từ 0, 25 đ ến 0.5 k m / k n r O' nhữnii cao noLivên đá vôi kh ỏ han (Sơn

La-7

Mộc Châu, Đ ồ n g V ã n ) cho đến trên 1,55 k m / k m 2 ớ nh ưn g nơi m ư a nhièu và có địa hình bị chi a cắt mạnh.

Hệ thố ng sông H ố ng do ba sông: Thao Đ à và Ló hợp thành t ron" đó sỏni:

T hao được coi là d òn g chính (hình 2.4).

D òn g ch ính s ô ng H ổ n g bất ng uồn từ hồ Đại Lý ớ độ cao eân 2000 m tron đính Ngụ y Sơn thuộc tính Vân N a m - T r u n g Q uốc , chảy theo h ư ớ ne tây bă c- đon e nam q u a tính Vân N a m rồi đố vào nước ta tại biên giới V i ẽ t - T r u n s thuộc Phô Lu huyện Bát Xát tính Lào Cai; tiếp tục cháv q u a Lào Cai Yên Bái, Hà Nội Hưnti Yen rối đổ ra vịnh Bắc Bộ tại cửa Ba Lạt Phần n ằm trẽn lảnh thổ T r u ng Q u ố c dược ÍỈỌÌ la sòng Ng uyên, đoạn trung lưu từ biên giới đến Việt Trì được íiọi là s ô na T hao và ŨI' Việt Trì đến Ba Lạt được gọi là sông Hồng Tính từ n g u ô n đến cửa Ba Lại dò nu chính s ông H ồn g dài 1126 km, trong đó 556 k m cháy trong lãnh thổ nưóc ta ơ ha lưu, s ông H ổ n g có m ột số phân lưu n hư các sông Đ á y, Đ u ô n g Luộc, Trà Lý, Đ ào và Ninh Cơ Sông Đ á y vốn là con sông tự nhiên có cửa Hát M ô n đưa nước từ sõn e Hông cháy vào sông Đáy T ừ sau khi đ ập Đ á y được xây dựn g (vào năm 1937) CƯU song Hat Món bị bôi lap nên sống Đáy dã trớ ihanh con s ô nạ tiẽu nước chI khi lũ sõng Hỏ ng đặc biệt lớn đế đ ảm báo an toàn cho thủ đô Hà Nội và đ ỏ n s hãng SỎHL! Hồng Ngoài ra, các s ông Cà Lổ Ngũ H u y ệ n Khê trước đáy c ũ n e là phán lưu sons: Hổng ớ phía bờ trái n h u n g nay, hai cứa sòng nàv dã bị bói láp nên chi con la sòIiii tiêu nước cho vùn g Vĩn h Phúc và Bãc Ninh

Sòng Đ à và sông Lô cũn g bát n gu ồn từ tính Vân N a m TrunSI Ọuóc Phun Ihương Ill'll năm trong lãnh thọ Trung Q u ố c cứa s ó n s Đ à đitợc tiọi la s ó ne Lý Ticn còn của s ông Ló được ơọi là sông Bàn Long

Dòn g chí nh s ô ng Đ à cháy vào nước ta tại xã Ka Lon tỉ huyện Muờníi Tc tinh Lai Châu, tiếp tục c h á y q u a ihị trấn Mường Tè và các thị xã Lai Chau, Sơn La vu I lòa Bình rồi đổ vào s ô ng T h a o tại T run g Hà Từ n e u ồ n đến T r u n g Hà, sỏ nc Đ a dai I

010 kin diện tích lưu vực 52 900 k n r trong đó phán n ằ m tron í: lãnh thổ nước la dài

17 200 k m : , trong đỏ ph ấn trong lãnh thố nước ta dài 217 k m diện lích Ill'Ll vực

9 78 0 k n r Sông C h ay bát n gu ồn từ sườn phía táy bắc dãy T á y Co n Linh thuộc huyện

Su Phì tinh Hà G i an g, ch áy theo hướng tâv b ắ c - đ ỏ n s n a m đ ến phía dưới Si- Ma- Cao

h uyện Bắc Hà tinh Lào Cai thì c hu yển hướng đ ô n g bắc- tây n am , chay q u a Bao Nhai Rao Yên ròi đổ vào hổ chứa Thác Bà Sau dó phẩn ha lưu s ô n e C ha y đo vào s on e Lo tại Đ o a n Hùn g Sông Ch áy dài 319 km diện tích lưu vực 6 500 k n r trong đỏ có

1920 k n r n am trên lãnh thổ T r un g Quốc

8

T ổ n g diện tích lưu vực của hệ t hống s ống H ổ n g k h o án g 155 0 0 0 k n r trong

đó chí có 72 7 0 0 k m 2 ( 4 7% diện tích) n ầ m trong lãnh thỗ nước ta Diện tích lưu \ ực sóng H ổ n g tính đến Sơn T â y k ho ản g 143 700 k m 2 trong đó lưu vực s ô n s Đà chiêm 36% s ông Thao: 3 7 , 5% và sóng Lô: 26,5% Tí nh đén Sơn Tây sônu Hỏnii co 615 sõng suối các loại có chiều dài íừ 10 km trớ lẽn trong đó có 57 s õ na n hánh cáp 1

195 sông n h á nh cấp 2 245 sông n há nh cấp 3, 97 sông n h án h cấp 4 18 sons: nhánh cấp 5 và 2 s ó n g n h á nh cấp 6

1.6.2 Đãc điểm đoạn m ạn g sông nghiên cứu

Đoạ n m ạ n g song nghiên cứu được tính lừ các trạm thúy vãn Hoa Bình-Yen Bái-Vụ Q u a n g đến các trạm thúy văn Hà Nội -Thượn g Cát; bao m')in: đoạn son e Đa

từ Hoà Bình đến n gã ba T hao- Đà, đoạn sòng T h a o từ Y ê n Bái đến Việt Tri đoan sông Lô từ Vụ Q u a n g đến Việt Trì đoạn sòng H ổ n g từ Việt Trì đến Hà Nội va đoan sông Đ u ố n g từ ngã ba H ồ n g - Đ u ố n g đến Th ượng Cát (hình 1.2)

Phía t h ư ợ n e lưu trạm Hoà Bình, sông Đ à ch áy theo hướnn t â y -đ ô ne n hư n s từ

I loa Bình dẽn T r u n g Hà do ánh hướne của các ngọn núi rán chăc nhu Ra Vi Yicii Nam Đ ò i Thỏi làm s ô ng Đà đổi hướng, chay theo hướrm nam-bá c Lòn tí sỏne (í đoạn này m ớ rộng rõ rệt ( khoáng 200 m t r o n 2 m ù a can) nên thuận lợi cho iinio thrum thủy Độ dốc lòng sông eiáin nhưng khô ng nhiéu, còn k h o á n a 0.42 %"

Đoạ n s ông T h a o từ Yên Bái đến n e ã ba T h a o - Đ à van tiếp tục chay theo huứni! lây h ắc -đ ôn a n am n h ư n e riêmi đoạn lừ ngã ba T h a o - Đ à đến Việt Trì thì đổi hướnu

1'hay theo hướng lâv n am -đ ỏi m bác Lòng sóng đoan nay mó' r ộn" tới 300- 400 m nen vặn tái trẽn sónt! thu ận tiện hơn Địa hình hai bên hờ ha tháp x u õ n a đưói 25 111 Bãi nổi \ LI at hiện nhiên và thườn £ bị nước lũ tràn neập Tại dây đã xuất hiên nhữnu CỈIKII1

đô đáu tién cua sóng Hổng

Từ Vụ Q u a n g đến Việt Trì sau khi nhận th êm nước cua sõníí Ciãm sỏny Lo vần cháy theo hướng tây b ắc- đỏn g nam Trước khi đổ vào Việt Trì s ó n s Lu u m nhàn them một n h á nh bên trái nữa là sồng Phó Đáv L ò n g s ô n e đoan từ Vụ Q u a n u đến Việt Trì mỏ' rộn g tới trên 200 m nên t huyén bè đi lại khá thuận tiện, thác yhẽnh không còn chi xuất hiện các bãi nổi

Kế từ Việt Trì s ôn g H ồ n g bắt đầu cháy vào v ù n2 đ õ n s b a n s n h u n s \'ăn tiep tục theo hư ớng tây b á c - đ ô n g n am ra tới tận cứa Ba Lạt Dưới thị xã Sơn Táv chừnii

I 2 km SÔI12 HỒ112 tách ra mộ t phân lưu là s õ na Đáy Phân lưu sónii Đ á y chi phan lũ khi có kì lớn đế d a m bao an toàn cho thu đô Hà Nội và đ ó n s b a n e sónii Hóng Tại hác bãi C a m Lò có s ỏ n s Cà Lô nối s ông H ổ n a \'ới s ông c á u n h ư n s h iện đã bị cát hỏi Ví) đò vít kín K h o á n g 8 k m trước khi tói Hà Nội m ộ t ph án lưu nữa lá sõnu Đuonìi nòi s ỏ ne H ổ n u với SÔI12 Thái Bình Trước đáv, ph án lưu s o n s Đuốnii rất nho n hu nu

nì l e)56, cứa s ò n g Đ u ô n g được m ớ r ộng \ à s óng Đ u ố n g irớ thánh phán Ill'll LỊLKtn

I r o n u n h á t c u a SÕ1112 H ổ i m

9

H ình 1.2 Sơ đổ đoạn m ạng sò n g n g h iên cứu

LỔ.3 Tình hình nghiên cứu thuỷ ván và các s ố liệu đ ả thu thập

Việc q u a n trắc thuý văn đã được tiến hàn h tại mộ t số đ i ể m trẽn sỏrm Đa nhu Hoà Bình (từ 1902) Lai C h âu (từ 1927 ) ; trên s ông T h a o n hư Yen Bái Tu y cn Ọuaim, Phú T h ọ ; trên sông Ló n hư G h ề n h Gà Vụ Q u a n g T u y nhiên, lại da so các trạm truức dây chí tập trung q u an trắc mực nước Lưu l ư ợ n ” nưức khónti đuuv quan irăc liôn lục và thường chi được tiến h àn h trong giai đ o ạn m ù a lũ

Sau n ă m 1956 việc q uan trác liên tục c h o hê t hố n g các yêu tở ihuv vãn chính như mực nước Ill'll lượng nước, bùn cát, nhiệt độ n ư ớ c được tiến h àn h ớ một su vị trí q u a n trọng; trên s ôn g Đ à n hư H o à Bình T ạ Bú Lai C h á u trẽn sỏ ne T hao như

10

Yên Bái Phú T h ọ , trẽn sòng Lô n hư G h ẽ n h Gà V ụ Q u a n g và trẽn s ôn c Hỏni: như Sơn Tây H à N ộ i

Cô n g trình thủy điện H o à Binh được khởi c ỏ n e xây d ự n s n gày 6 / X / \ L) l L>: lap

sóng đợt I vào t há n g X I I / 1983; lấp dẫn dòn g thi c ô n g vào tháng 1/1986 và hoàn thành vào thá ng XI I/1994 N ga y từ khi khởi c ô ng xây dựng, c h ế độ thuy vãn cua sóng Đ à dã có nh ững thay đổi nhất định Tu y nhiên, n h ữn g thay đỏi rõ rệt cua đòn tí cháy nước và d ò n g cháy bùn cát được ghi nhận ké từ thá ng 1/1986 khi lấp dần clnnu

đế thi cô ng và đặc biêt là từ 1988, khi hồ tích nước đê c hạ y tổ m á y số 1 Cĩìnu lừ năm 1988, do vị trí trạm thuý vãn Ho à Bình ỏ q u á gần đ ậ p chính, chê độ đ ò ne chay chịu ánh hư ớng trực tiếp từ cóng trình, kh ô ng đ ảm bảo cho c ỏ n s lác đo đạc vì the

tram c hu yể n về phía hạ lưu k ho án g 4,5 k m và m a n g tên mới là Bến Naọc

Bảng 1.3

T h ô n g kè tin h hinh s ố liệu th ủ y văn th ự c đo đã thu th ậ p

Vi trí tram

Thơi kì quan irác

1

Vu Ọ u a n s

■ 105"17 20 21"30'50"

Ọ p, t" nu ức bình quàn ns,ày i m - 2003

+ Thời kỳ trước khi có hố H o à Bình: từ 1985 trớ vé trước (coi anh hirớnti cua

hò tronu thời kỳ 19 79 - 1985 là kh ô ng đ á n s kê)

+ Thời kỳ sau khi có hồ Hoà Bình: từ 1986 đến nay

Đ ế p h ụ c v ụ đ é t ài n e h i ê n c ứ u , b á o c á o d ã t h u t h ậ p c á c s ố l i ệ u YC doDLi c h u \

nước nhiệt độ nước và lượng n g âm cát c ủa 6 trạm là: H o à Bình Yen Bái Vụ Ọ ua nu Son Tây, H à Nội T h ư ợ n g Cát thò'i kỳ từ sau khi có anh hưứnìi cua ho Hoa Binh (tư nam 1986) Tì nh hình sô liệu đã thu thập đươc t h ô n s kê tronti b a n s 1.3

1.7 Quy trình vận hành và một sỏ thõng sô cơ bàn của hó Hoà binh [5]

1.7.1 Q uy trinh vận hành hổ Hòa Bình

Hồ H ò a Bình được thiết k ế đế có thế cắt được những trận lũ lớn nhu lũ lịch SƯ' tháng VIII n ãm 1971, đ ám bảo cho mực nước H à Nội k h ô n g vượt quá 13,3 m Đê nâng cao hiệu q u á c h ố n g lũ cho hạ du và phát điện, g iả m thiếu mâ u thuần siiữa các ngành d ùn g nước, cần phải có quy trình vận h à nh hợp lý N g à y 12/VI/1997 Bộ trưởng Bộ N ó n g n g h i ệ p và Phát triển nô ng thôn Tr ường ban Phòng chốniỉ lụt bào

Tr ung ương N g u y ễ n C ốn g Tạn đã ký quyết định số 57 P C L B T Ư / Q Đ vé qu y trình vận hành hổ H ò a Bình [7] với một sỏ nội d un g có thế tó m tát n h ư sau:

Điều 6 q uy định 3 thời kỳ vận hành: (1) Thời kỳ lũ sớm lừ 1 5/VI - 15/VII (2 I Thời kỳ lũ chí nh vụ từ 16/VII - 25/VIII, (3) Thời kỳ lũ m u ộ n từ 26/VIII - 15/1X

Đ iéu 7 q u y định khi hồ chứa làm nh iệ m vụ cát lũ cho ha du mực nước tron Lì

hổ từng giai đo ạn cho phép duy trì ở mức tối đa như sau:

/ ) T l i ờ i k ỳ l ũ s ớ m l ừ J 5 / V I - 1 5 Ỉ V I I :

+ ) Từ 1 5/VI - 15/VII: mực nước hồ được điều tiết ớ cao trinh 91 = 2m

+ ) Trước n g ày 25/VI do yêu cầu cắt lũ tiếu mã n đ á m báo khòm: xa xuốnu ha Ill'Ll Ill'Ll lưựng lớn hơn 4 0 0 0 in' /s thi có thê nãn g mứ c nước hổ cao hon nhmiii khoiiLi quá cao trình + 1 0 0 m và phái đua về mức 91 X 2 m trước n gày 1/VII

2) Thời kỳ lũ chí nh vụ từ 16/VII - 2 5 /VIII:

Mực nước hổ được điều tiết ớ cao trình 9 1 + 2 m

j I I hời kỳ lũ niiiộn ỉừ 2(ìi\ I I I - 15/1X:

Tuỳ theo tình hình thời tiết cuối tháng VIII dự báo đáu t h á n e IX v;i kẽì qiui tích nuúc làm nh iệ m vụ cắt lũ của hổ t ro ne thời kv lũ chí nh vụ, Ban chi dạn PCLBTU sẽ trao đối với T ổ n g Cóng ty Đi ện lực Việt N a m và Tổn tỉ cục Khí tưony

Th uý văn đê q u yế t định

1.7.2 M ọt s ố thôn g so co bán cùa cóng trình thủy điện lìũà bình

Mộ t số t h ông số cơ bản cúa c ô n s trình thủv điện H ò a Binh đựơc trình bày trong bang 1.4

B áng 1.4 M ột số th ô n g sô cơ bàn của hố Hòa B inh [8]

1 Nhà máy:

- Kích thước cửa xả đáy m X m 6 X 10

C h ư ơ n g 2

C ơ S Ớ L Ý T H U Y Ế T C Ú A M Ô H Ì N H H E C - 6 [11]

2.1 Giới thiệu mô hình HEC-6

Wi l li am A T h o m a s là người đầu tiên có côn g thiết kế xây dựng và phái I] lên

mỏ hình HEC- 6 Phiên bản đầu tiên là H E C- 6 2.7 được xây d u n e tù nãm 1476 tại Trung tám Kỹ thuật T h ú y văn (Hydr ologic En gi n ee ri ng Cent er - HEC) cua Hoa Kỳ

Từ đó đến nay, T h o m a s và các cộng sự đã liên tục phát t n ế n m ó hình và cho ra đời nhiều phiên bán k h á c c ua m ó hình với sự bổ sun e m ột số tính năng Báo cáo này sư dụng H E C - 6 phiên bản 4.1 cồ ng bô năm 1993 Phiên bán này dã được bổ sunti tính năng m ô ph ỏng tính toán vận c huyển bùn cát trên m ạ n g sóng và c ho các hạt bùn cát

cổ đường kính lớn

2.1.1 M a c đích và nguyén lý của mớ hình

H E C - 6 là m ỏ hình số thúy động lực một chiêu trong l ò n s dẫn hơ có bién động, được thiết k ế để m ỏ p hó ng và d ự báo những biến d ạ ng c ủa trãc diện dọc sone xay ra do xói ho ặc bồi Irong thời đoạn tính toán vừa Số liệu quan trác d ònn cha> lid’ll tục được phân th ành chuồi các dòn g chay ổn định có lưu lượnẹ va thòi tiian đuv trì thay đoi Đ ườ n g mặt nước được tính toán đối với mỏi Ill'Ll lương va tứ đo cho đn dốc năng lương, lốc đỏ đỏ sâu lại mỗi mặt căl nuanii Sail đo sức tai cái đươc tính toán tại mỗi mặ t cát Sức tái cát này kết hợp với thời gian d uy trì dnniỉ chay cho phép tính thê tích bù n cát trong mỗi đoạn sổng Ti êp theo, lượng hỏi hoặc xoi tại mỏi mat cát dược tính toán và mặ t cát n g an g được hiệu chì nh lại Sau đỏ các tính toán này được thực hiện đối với lưu lượng k ế tiếp trong dãy và chu kỳ được lãp lại lừ đau với những t h òn g tin hình học mới nhất Các tính loan bùn cát đưực tiên hanh ihen ùrnt: phấn kích thước hạt và băng cách đó cho phép m ò phỏ ng thúy lực quá trinh \ Ỏ1 hoặc thô hoá lòng sõng

2.1.2 K h ả năng ứng dụ n g và các hạn chê của mó hình H E C -6

Mó hình H E C - 6 có thê được ứng d ụn g đế:

thiết k ố m ò p h o n e xu thế bổi xói l òna dẫn t ro na thòi hạn dài do kết qua cua SU' t h a y đ ổ i c ó t í n h t h ư ờ n g x u y ê n v à c h u k y c ủ a Ill'Ll l ư ợ n 2 h o ặ c d o s ư t h a y đo i hình d ạn g l òne dẫn

- đ á nh íiiá hổi láne trong hổ chứa

ihiết k ế thu hẹp l òn g dẫn, đáp ứng yêu cáu d u y trì chi êu sáu vận tài thủy

d ự báo ánh h ưứng cua nao vét đen tốc độ hổi l áns

L1'ỚC tính kh a n ã n s xói lớn nhất t ro n s tr ư ờn 2 h ợp 10 lớn và đánh íiiá quá trinhbổi xói trong lòng d ẫn co định

14

kh ổn g tính được vận c hu yên bùn cát theo phương n gang nén khóntì có kha năng tính xói lớ bờ và m ô phỏ ng sự phát triến c ú a n hững đ o an s ó na con tí.

- chỉ áp d ụ n g được cho s ông k h ô n g chịu ánh hư ớn g triều,

- 2 »

trong đó: g là gia tóc trọng lực; lic ỉà tổn thài n ăng lượng; Vị v : là tốc đô t run li hình

tại đấu và cuối đ o ạn sông; WS |, W S 2 là các cao trình mặt nước tai đáu vù cuối đoạn sông; 0C|, Gt- là hệ s ố phân bô tốc độ đối với dò ng c há y tại đầu và cuối đoan snniz

" ^ O ư ơ n g đ ỏ d o c r E n g l t Ị é n g -I— - Đ ư ơ n g m â t n ư ơ c I

D on g c h à y t * ■

Đ a y lo n g d ẫ n

H inh 2.1 Các sò hạng của

ph ư ơ n g trin h năng lượng.

2.2.2 Tôn thất thúy lưc

2 2 2 / T ố n tlníỉ (lo m a sát

Tổn thất n ăn g lượng l \ bao e ổm tổn thất do m a sát hj và tổn thất do hình đạn ổh„:

Tron li số hạn a tổn thất do hình dạng, chi xét các tốn thất do co h ẹ p và m ớ rộn II

Đ ể x ấp xi phân phối d ò n s cháy theo p h ư ơ n s n e a n e s ỏ ne được phan chia thành các dái t ư ơ n s tụ' nhau vê tính chất thiiy lực t heo p h ư ơ n s clone chay Mói mặt

15

cắt n g a ng được ch ia n h ỏ thành các tiểu m ậ t cắt Tổ n thất do m a sát được tính theo công thức sau:

R, + R,

’ I n i J L

-trong đó: A | , A-, tương ứng là diện tích của dò n g ch áy thõng t h ư ờ n2 tại các mật cài thượng và hạ lưu; NSS là tống số các tiếu mặt cắt của mặt cát n s a n g : K là mô d u n lưu lượng của tiếu mặ t cắt: Lị là chiều dài c ủa dái thứ j giữa các tiếu mặt cát: n là hẹ

sổ nh ám Mann in g: Q là lưu lượng nước; R| và R ; tương line là bán kính thúy lực phía thượng và hạ lưu

r õ 11«.

2.2.3 Tính toán các yếu tó thủy lực

Mỗi mặt căt n g a n g được định nghĩa bới các toạ độ ( X Y) n hư thé hiện troniĩ hình 2.2 Đ ế t hu ậ n lợi cho việc ấn định các giá trị n, chi ều dài sòng.,.: mỗi mặt cat ngang được phán c h i a thành các tiểu mặt cắt, thường g ồ m m ột lòng chí nh với các hãi trái và phái

16

Diện tích, chu vi ướt bán kính thúy lực m ô đu n lưu lượng, hệ sổ phán bò tốc

độ độ rộng và độ sáu hiệu q u á c ủ a mỏi tiêu mặt cắt được xác định từ các phan 10 hình thang n hỏ hơn, được giới ‘hạn bới mặt nước và các đ iế m toa đổ ké 1 ién cua mat cất ngang

Sức tái cát c ủ a d ò n g cháy trong lòng dẫn có mặ t cắt n a a n s kh ôn g phai hình chữ nhật được tính n h ờ sử d u n g độ sâu hiệu quá E F D và độ r ông hiệu qua EFYV lính

từ độ sâu hiệu q uả đế báo toàn A ( D :/:<) cho mặt cắt ngang Đ ộ sâu hiệu q ua dược tinh theo cóng thức :

Việc lính sức tái cát chi dựa trên các tha m sử thúy lực c ua chi lone dẫn chính nên các yêu tố thúy lực c ũ n g chí được xác định theo tài liệu địa hình eiới han tron" phạm vi của lòng đẫn chính

Tai mồi mặt cát n s a n s H E C -6 kiếm tra các tha m sô thúy lực thích hop de xác (.tịnh \ c m đ o an SÕ11£ là đ o ạn dòim chúv êm hay xiet Néu c h á \ Om các tính loan tiến hành theo c hi ên ngược d o n e còn nêu cháy xiêi H E C - 6 xấp XI hình d a n s hình học lòna dần n h ờ độ sâu và độ rộng hiệu q uá và xác dinh cao trình mặi nước dưa trẽn độ sâu c há y xiết t h ô n c thường

Ba bước chí nh đươc sứ du ng đế hội tụ các tính toán thử sai trong khi tính cao trình mặt nước c ủ a m ặ t cắt n gan g thượng lưu Hình 2 4 thê hiện dãy các thư sai liên tiếp đế hội tụ p hương ph áp bậc thang t hông dụng Q u y trình tính toán nhu sau:

* Th ử sai lần 1: Dựa trên cao trình mặ t nước ban đầu

* Th ứ sai lần 2: Giả thiết lượng thav đối bàng 9 0% c ủa AY,

* Th ứ sai lấn 3: Các cao trình cúa bước thử sai lần 1 và 2 đã gia thiết được s h é p nốivới nhau b ăn g mộ t đường thẳng và các đáp án c ua các tính toán thứ sai lán I và 2cũng được g h é p nối với nhau bằng một đường thắng G i ao điếm c ua hai đườrm thane này trở thành giá trị giả thiết c ủa thử sai lần 3

Q u á trình này tiếp tục cho tới khi các giá trị gia thiết và tính toán cua cao trình mặt nước n ăm trong phạm vi dung sai sai số có thế cho phép Khi đó cao n i n h mật nước tính toán trớ thành nghiêm hôi tụ

Nêu sau 20 lần thứ sai mà sai số vẫn lớn hơn sai số cho phép thì cao trình mật nước tính toán cuối c ù ng được sử dụng

0 WSE tính toán 1

* W5E giả thiẻt I

-Dung sai sai sổ cò thế

Só thử sai

H inh 2.4 Sự hội tụ của các cao trin h m ặt nước

aiả th iế t và tính toán

2.2.4 Các tham sô thủy lực điển hình sử dung trong các tính toán bùn cát

Các t ha m NÒ thuy lực được biến đối thanh các eiá trị điên hình dối với m ỗ i

đoan s ồng trước khi tính sức tái cát Các phương trinh tòn ti q uát dược trinh bày tlưói đây Có thê sửa đổi các trọng số cù ng với các sò liệu đầu vào

B àng 2.1

Một vài trọng số khác được nghiên cứu trong khi xây d ự n s c ô n e thức cua NO'

đồ tính Báng 2.1 trên trình bày bô hệ số cho tính toán gần n hư ổn định nhát và nho

đó cho phép các bước thời gian dài nhất (sơ đổ l ) và bộ hệ số n hay nhất đói với việc thay đổi c ao trình đáy nhưng lại véu cầu các bước thời gian n a án hơn đe được on

định (sơ đổ 2) Sơ đồ l t h ư ờ n s được chọ n nhiêu nhất bới vì độ dốc nãny lƯơnỊ 1 tính toán cổ thê thay đổi m a n h lừ mặt cắt này đến mãt cắt khác Irons khi cliẻn biến sônìi ihựt lè có thể phụ thuộc \ à o các đặc t runs cua đoạn

2.2.5 Độ nhám thủy lực

Đ ộ n h á m biên c ủa s ôn g bồi tích có q u an hê chặt chẽ với vận c h u y e n bùn cát

và c huy ên đ ộ n e cíia vật liệu đáy Cấc tổn thất n ã n s l ư ợ n2 cua đirờntí mặt nu(íc phai bao gổin các ánh hiro'ng c u a tất cá các tổn thát như: tốn thát do độ n h á m hat cua đ á \

đ ộ n g ; tốn that do SƯ cán trớ gây bới c á c hình d ạ n 2 đ áy n hư s ó n a HỢn đụn cái l í n h

k h ổ n g đẽu c u a bờ, thực vật; các tòn thát do thu hẹp hoặc m ơ rộng: các lỏn thát do

k húc s ô n g c o n g và các tổn thất do các chỗ nổi T rừ tổn thất do thu hep va m ớ r ó n V

tất cả các tổn thất kh ác đổu dược bao g ồ m trong th a m số độ n h á m đơn Manniniỉ n

19

2.3 Cơ sở lý thuyết tính tơán bùn cát

2.3.1 Phương trình liên tục của vật liệu bùn cát

2.3.1.1 K h á i n iệ m t h ể tích k h ốn g c h ế

Mỗi mật cắt n g an g m ó tả một thế tích kh ốn g chế Độ rộng c ủa thê tích khóntỉ chê thường b ãn g độ rộn g lòng động và độ sâu cua nó kéo dài từ mặt nước đèn dinh cua đá góc hoặc đ ến điếm k hố n g c h ế địa chất khác n ă m dưới hê mặt đ áy sons Tr on e các diên tích k hô n g có đá gốc, mộ t siới han tuỳ ý (gọi là “ đáy mô hình) được ấn định (hình 2.5) T h ế tích kh ốn g c h ế đối với mật cắt n g an g 2 được mô ta bới các đường đứt nét to Các thể tích khố ng c h ế đối với các mật cắt 1 và 3 được noi với nhau bới mặt cát 2 Phương trình liên tục bùn cát được viết cho thê tích kiêm soát này còn p hương trình n ăn g lượng lại được viết giữa các mặt cãt n ean c

T ố n g lượng bùn cát trên lòng sông được xác dinh h ăne cách xáp XI diện tích giới han trung bình n hư sau:

trong đó: lị, là độ r ộng lòng động: L ư Lứ tương ứng là chiêu dài c úa (.loạn thương và

hạ lưu sử d ụ n g Irong tính toán thể tích k hố ng chế, \ ’SO| thế tích hun cát trong thê tích không c h ế và Kv là đ ộ s â u c u a b ù n c á t t r o n g t hế t í ch k h ô n g chõ

Đôi với mỗi đ ộ sâu D thể tích chất l ỏna trong cột nước là:

Các th am sỏ thiiy lực độ rộng B„ vù độ sáu D sứ d u n g để giai phương trình nàng

lượng được tính b a n s c ác h láv trung bình trẽn k h o a ng k h ô n g gian n hư nhau

Các th a m số thúy lực cấp phối vặt liệu đáy và sức tai cát lính toán đirơc gia thicí là đ ổ n e nhát trong ihẽ tích kh ủn g chẽ L.irnìiy hun cát cha_\ vào đưiíL gia ihici la

2 0

xáo trộn đều với b ù n cát hiện có trong thể tích k h ố n g c h ế và sư k h uế ch tán cua tất ca các lớp kích thước hạt đươc xem n hư xáy ra ngay lập tức.

Sõng trao đổi hùn cát với các biên cùa nó theo cá hai chiêu t hã n s đứ n° cũ na như năm n g an g và thay đối hình dạng c ua nó băng c ác h hình thành các lòne dẫn các

đê tư nhiên, các khúc cong, các đáo và các hình d ạ ng khác trên mặ t bãntỉ T u \ nhien

H EC- 6 chi m ó hình hoá được sự trao đổi bùn cát c ù ng sự thay đói 1ÒI12 s on e iheo chiếu thảng đứn g c ò n độ rộng và độ sâu của nó do ngưò'i sử d ụ n s x á c đ ị n h

H E C -6 tính đên hai nguồn bùn cát: bùn cát tr one nước chay vào \ à hùn cat long sông Lượng bùn cát cháy vào là điêu kiện biên và được m o ta L'Cmii với sô liệu đầu vào Thê tích k h ố n g c h ế bùn cát lòng sông c u ng cấp thành phán bẽ phót nụuòn

và cũng được m ỏ tả tr ong số liệu đầu vào

Lý thuyết vận c h u y ể n cát và các kích thước lớn hơn liên kết quan đến suất vặn chuyển c ù ng với cấp phối hạt bùn cát trên bề mặ t lòng sông và thúy lực dòng cháy Các tính toán hoá thó yêu cầu cấp phối hạt vật liêu bén dưới bé mặt lònu Nona Độ sáu cưa đá góc hoặc cua mộ t sổ vật liệu khác có the nuãn ngừa suy thoái cũ n e phai dược nhận d ạ n g dế giới hạn quá trình xói Các Yêu cấu này được vào HLỈC-6 hãn ì: cách tính toán tách riêng cấp phối hạt hê mặt lòng sông và dưới mặt lòn tỉ Sony

Các toa độ được nối bới các đường liên nét t ron a hình 2.6 xác định hình đ a n e mặt cắt n e a n a ban đầu tại lúc bát đầu mô phỏnti Đối với các điêu kiện X Ó I, hiệu iiiừa lượn ti bùn cát c h á y vào và khá n ã n a vận c huy ến c ua đ oạn sònn được biên thanh the tích xói Nêu c ao trình đáy m ó hình k h ỏ n s được xác đinh trong các diêu kiện ban dầu thì giá trị m ạ c đinh b a n s 10 ft sẽ được sư d ụ n s sau đó trư thành đó sâu cực dại

đê xói hiện có của vật liệu lòng sông

2.3.1.2 Plìiùniii trình liên lục của bùn cát

Phươnn trình liên tục cua vật liệu bùn cát ( phương trình I : \ n e n la co SO' đc mó phong c h u y ê n d ỏ n g t heo ph ương i hã ne đứng c ua lòng sõng:

L ò n g d ẳ n

-CCorao

K hoả ng cách

H in h 2.6 Vật liệu bùn cát tro n g đáy sông

21

trong đó: ổ c, là độ rộn g lòng động, t là thời gian G là lun lượng bùn cát (iV/s) trung bình trong hước thời gian At X l à k h o ản g cách dọc theo lòng dẫn và ì \ là đỏ sâu bùn

Hmh 2 7 Lưỡi tính toán

Đ ô sâu vật liệu đ áy ban đấu tại điế m p định n gh ĩa 2Ĩá trị ban đáu cua YLưu lượng bùn cát G u là t on e lượng bùn cát theo kích thước hạt vào the tích kiếm soát từ thể tích k iể m soát phía thượng lưu Đối với đ oạn s õng xa nhất phía thương lưu đây là điều kiện biên lưu lượn 2 bùn cát cháy vào do n lĩ ười sư d u n e cung cáp Bùn cát rời khỏi thó tích kiêm soát G trớ thành G u đối với thó tích kiêm soát tiép theo phía hạ lưu

Lưu lượng bùn cát Gj được tính b a n s cách x em xét k ha n ă n e vận chuyên tại

đ iểm p, d ò n e bùn cát vào vật liệu hiện có trẽn l ò n e s ô n s và sư hoá thỏ Sư khác nhau ai lìa Cij và Ciu là tổng \'ật liệu hồi hoặc xói trong mõi đo ạn s õ na yhi nhãn " \ ' ì m2 tính t dá n " trẽn hình (2.7 I và được biên đổi thành sư thay đổi cùa cao trình đ á \ nh ơ SƯ dụng phương trình (2.20)

Khá n ãn u vận c h u y ế n c ủa mỏi kích thước hạt được tính toán đối Ml] các dicu kiện thủv lực tại đáu thời đ o ạn và k hô n g được tính lại Irons: thòi đoan nay Bin váy ihời đ oạ n tính t oán phái đu n s ă n đế các thay đói can trình đ á \ do xói hoặc bói iron 11

thời đoạn k h ỏ n g ảnh hướng đ á ng kể đến khả n ă ng vận c hu y ê n tại cuối thời đoan Bước thời gian điên hình dối với các lưu lương lớn là phần cùa ng ày còn đối với dònn cháy kiệt là vài ng ày hoặc thậ m chí vài tháng vẫn có thò thoa mãn Ván đẽ quan trọng cần đánh giá là tống thay đối cao trình đáy có thẽ được chấp nhận tron£ moi bước thời gian N ó phải k hông vươt quá lượng thav đối đáy cho phép i r o n s một thòi đoạn tính toán là 1 ft hoặc 10% độ sâu dò ng nước Tuy nhiên, cấp phối vật liệu day được tính lại trong thời đoạn vì tổng lượng; vật liệu vận c h u v ế n rất nh ay đối với cấp phổi hạt vật liệu đáy.

2.3.1.4 Tính toán lại cấp p h ố i hạt vật liệu dá x

H E C - 6 giái phương trình Ex n er đối với SU' liên tục cùa bùn cát Nêu kha rũuiii vận c hu y ể n lớn hơn lượng bùn cát vào thê tích k iém soát, bùn cát hiện có được di chuyến từ đáy đế t hoá m ã n lính liên tục Vì kha n ã n e vận c hu yế n đối với một kích thước đã cho phu thuộc vào phán kích thước hat trẽn đáy nên cán thưòìm xiivên tính lại các phán kích thước hiện có do bùn cát trao đổi với đáy Vẽ mặt lý luận, so các lượng gia trao đổi SPI trong một bước thời gian có q uan hê với chiêu dài bước thòi gian At với tốc độ trong một đoạn sông V và chiều dài đ oan sônti L hơi:

Thường thường, lượng gia trao đối có thế nhó hơn aiá trị này mù k h ô n c phát sinh các vấn đế s ố đán í! kê SPI han đáu phái được đạt bằni! 0 và một sự kiện ihuy văn cực đoan đưực m ỏ phỏng Đây phái là trường hợp ổn dinh nhất Sail dó bãt dâu lừ

SPI = 5 0 h o ặ c l ớ n h o n g i á trị S P I p h á i ui á m t h e o l i r ợi m ííia c u a lt) c h o d ẽ n k h i c ú c

kết quá trớ nên khác đ á n g kể các kết quá với SPI = 0 Sử sụn a SPI nho nhất lie cho một lời giai sát với SFI = 0

2.3.2 Xác định các lớp hoạt dóng và không hoạt dông.

H i iC - 6 su d u n u khái niệm lớp hoạt d o n e và khỏim hoat đ ộ n ” Lớp huại dnnti được gia thiết được x áo trộn liên tục với dòng chav n h ư n e có thế có mọt bé mặt các hat c hu yến đ ộ n g c h ậ m báo vệ các hạt mịn hơn khỏi bị cu ốn trôi vào cỉònn chay Hai quá trình k há c nh au được m ỏ phỏng: (1) xáo trộn xảy ra giữa các hat hùn cát lònti sòng và hồn h ợ p chất l ỏ n s - b ù n cát do năng lượng trong chất lóng c hu y ên độníi va (2) xáo trôn xáy ra giữa lớp hoại d ộ n s và lớp k h ô n e hoạt đ ộ n 2 do c hu y ế n đợnti cua he­

rn ặt lòng sòng Các cơ c h ế xáo trộn được quy vào sự rỏi dó n g quy mỏ lớn và ứnii suất tiếp đ áy do nước c h u y ế n động Chiêu sâu xáo trộn (aọi la "ch iêu sáu cán bũmi” ) đirực biếu diễn bàng m ột h àm của cường độ d ò n ° chay (kru l ượ n” đưn vị), dỏ dốc năng lượng và kích thước hạt

2 3 2 Ì C h i ê u sáu câ n băníỊ

Đ i é u k i ệ n t h ủ y l ự c n ã n g l ư ợ n a c u e t i ê n m à lại đ ó v ừ a till c h o m ộ t k í c h thưtVc

h ạ t r i ê n g h i ệ t k h ô n g c h u y ế n đ ộ n g t r ò n b ẽ m ặ t l ò n n s ó n t : d ư ợ c l í n h h ã n g c á c h kCi h o p

tirơnẹ ứns các phươnc trình M anning Stricklcr \ ’à Einstein:

trong đó: d là đườn g kính hat; D là độ sâu đ ò n e nước; V là tốc độ d ò n e nước: p là

khỏi lượng riêng c ủ a bùn cát; p, là khối lượng riêng cứa nước: y\i là arờiiii J ộ vận

chuyển trong h àm bùn cát đáy của Einstein, có qu an hc với số nạh ị ch đao cua tham

Vì vây chiêu sâu cán băng đối với kích thước hạt và lưu lượm: đơn vị d ã cho hãn LI:

trong đó Dc là c hi ểu sâu d òn g nước cực tiếu đối với vận c h u y ế n bùn cát khỏim đánL!

kế (tức chiểu sâu cân b a n s ) đối với kích thước hat <j

2.3.3 S ư thay đổi cao trình đáy

2 3 3 1 Lons’ d i m tĩúx d o n ”

Khi xói h oặc bồi xảy ra trong mộ t bước thời sian H E C - 6 điếu chinh các cao trình mặt cắt n c a n SI trong p hạ m vi phấn 1ÒI12 đ ộ n s cua mặt cắt n s a n c Đối với xói, phan l òne sõnti chí được c h u y ê n đ ộ n ẹ theo ph ư ơn s t h á n s đ ứ n s nốLt nó ư tnuiii pham

vi đáy đ ộ n a và ớ dưới mặt nước {nghĩa là ướt) Bổi l á n2 cho phép ớ rmoài ÌỈIỚI han vận c h u y ế n và ơ dưới mặt nước Một khi các giới hạn xói hoặc hôi được xác định,thể tích xói hoặc bổi được phân chia bới độ rộng hiệu q uá và chiéu dài cua thế líchkiểm soát đế nhận được sự thay đổi cao trình đáy Sau đó các thanh phán i hă ne đứnu của toan độ mặt cãt n e a n e trong phạm vi các eiói hạn xói hoặc hỏi n à \ dil'o'L diéu chinh n hư the hiện trong hình (2.8) và (2.9)

hiện b ăn g c ác h định rõ đ áy m ô hình E M B nhó hưn m ột lượng nh o so với cao trình

đ iếm sâu nhất ban đầu K h ô n g có bùn cát đ ón g g óp vào d o n e cháy bùn cát tại mặt cắt ngang

2.3.4 T rong lượng bổi lắng đơn vị

2 3 4 1 TrọiiíỊ -lượniỊ đơn vị ba n đuìt

T r ọ n g lượng đơn vị là trọng lượng của mỗi đơn vị thê tích chai bôi lane biéu diễn bâng trọng lượng khô

trong đó: Pc| là độ rỗn g của chất bổi lắng SG là trọng lực đơn vị cua c á c hạt hùn cát

7 là trọng lượng đơn vị c ua nước v à 7 là trọng lượng đơn VỊ cu a bùn cát

25

trong đó: B !à hệ số có két đối với hùn hoặc sét, T là thời eian tích I LI V troniỉ các nam

và 7, là trọ ng lượng đơn vị ban đau cúa bùn cát bôi l ãns, thường sal! 1 năm c ố két

T r o ng lượng đơn vị c ố kêi irung binh qua k ho án g thời aian T yêu câu lích hop qua thòi iíian T r ọ ng lượng đơn vị nàv được tính n hờ sứ du nu quan hẹ do Muller (1953) phát t n c n sail đây:

" • 11 T - ]Trọn*: lượm: đon vị này được sứ dung đê hiến đổi t rọna lượne bùn cái thanh the tích

đ ế t í n h t o á n s ự t h a y đ ổ i c a o t r ì n h đ á y s ó n g

2.3.5 Các đặc trưng của hạt bùn cát

Bốn đạc trưng c ơ bản của bùn cát: kích thước, hê số hình dạ ne t r o n ” lực n e n ”

và độ thô thúy lực rất q u a n trọng trong dự báo vận c hu y ến bùn cát

Các lóp kích thước hạt được ấn định i r o n s H E C - 6 theo hệ t h ốn e phán chia cua Hội Địa lý Mỹ ( A G U ) nhu' t r o n s b á n s 2.2 sau:

Bang 2.2 Phân loại bùn cát theo kích thước hạt

trong đó: a b c tương ứne là chi ều dài các true Irực giao lẫn n han dai nhái, trunii binh và n gãn nhất c ủ a các hạt bùn cát.

Hệ số hình d ạ n g hạt b ãn g 1.0 đối với hình cáu lý t ư ớ ns và có thố nhỏ hơn 0,1 đối với các hạt có hình d ạng rất bất quy tác H E C - 6 sứ d ụ n c một hệ số hình dạnn mặc đinh b ãn g 0 6 6 7 n h ưn g nó có thế được xác định bới người sư dung Nếu “đường kính bói l ã n g ” dược sứ dụng, nó sẽ được xác định b ăn g đặc trưng tốc đỏ rơi cua hạt

và hệ số hình d ạ n g hạt b ăng 1,0 được sứ dụng N ếu đ ư ò n a kính s à n2 thực té được su dụng thì hệ sỏ hình d ạ n g thực phai được sư dụng

T r ọ n g lực r iê ne c ủ a mộ t hạt bị chi phối hơi th à nh phấn k ho án g c ua nó Trong các hệ t h òng s ò ng tụ' nhiên, vật liêu đáy c hi ếm ưu t h ế là thạ ch anh có trọng lực rióng bãnti 2,65 H E C - 6 sử d ụ n g 2.65 như một mặ c định: tuy nhiên, các tỉiá trị tronư lực rièng c ủa cát bìm và sét có thê được vào

Hai kỹ thuật để tính độ thô thúy lực c ủa hạt có thế SU' d u n g trong HEC-6 Kỹ thuật thứ nhất dựa trên độ thô thúy lực xác định bơi Toffaleti (19 66 ) \'à tương tự với phương pháp R u b e y ( V an on i (1975) Phương p háp này gia thiết hệ số hình dạng bàng 0.9 Kỹ thuật thứ hai đưa vào xem xút hệ số hình d an g cu a hat d ùn g quy trinh

m ô la trong ICVVR ( 19 57 ) vù được m ỏ ta chi tiết bơi Wi l li am í 1980) Phương pháp thứ hai được m ặ c định

27

Khi ứng suất tiếp đáy tãng, xói hạt c h u y ể n th àn h xói khói và tốc dỏ xói tăn tí lẽn Về lý thuyết, vì tóc độ xói khối có thê là vỏ định nên Ar iathur ai và Krone í 1976)

dã dò nehị sứ d ụ n g mộ t "thời man xói đặc t r i m s ” T Với mộ t thời đ oan tính toán At.phương trình xói khối tro thành:

dưới các đ iều kiện xói HEC-óTciểm soát q uá trình này b ằng các h ki ế m tra ty lẽ phân trãm cúa sét trong đáy N ế u tý lệ này trên 10% thì tốc độ c u ốn trói sét hạn chê tỏc đỏ cuốn trôi hùn cát và sỏi c ũ ng n hư ngãn ngừa xói bùn, cát và soi trước khi xói sét cho

dù ứng suất tiếp đáy đủ đế xói các hạt nàv nhưng k h ô n g đu đê xói sét dính

Vì khái n iệ m “ n ồ n g đó cán b à n g ” c ủa Einstein được sứ d ụ n s đối với chái tai khô ng dính nên k h ô n g có sự c h ế ngư th em nào được yêu cáu đế niái hạn nontz đo cát

và sỏi Tu y nhiên, khi bùn cát dính được bao g ồ m thì sẽ k h ô n s có n o n e đỏ cân hãnti HEC- 6 giá thiết r ằn g xói và cuốn trôi các hạt m ị n bị han c h ế bơi một " n ò n e độ dims: cháy bùn cực đ ạ i ” N ồ n g độ dòng cháy bùn cực đai sư dụnsỉ bới H E C - 6 là 800 000ppm Nếu nô ng độ chất mị n (nghĩa là bùn và sét) tại mậ t cát ngunsi báì kỳ vượt quá

50 0 0 0 p pm, sự phán tác dụ ng tăng lên và một lời nhăn sẽ được in ra nói rõ t one sô lán các n ón g đ ó cao đã được phát hiện ra Khi nò n g độ Yươt quá 800 0 00 ppm mỏi nóng độ kích thước hạt được giam di mót cách tương line đé nònsd do tone cộntí la

h àm kết hợp c u a Toffalcli và Meyer -P et er và Muller,

q u an hộ cua Laur sen (1958) do M a d d e n ( 1985) sưa đòi

(1 96 2) dối vó'i bồi bùn cát dính do Ariathurai và Kr one ( 197 6) sứa đối

- q u an hệ L a u r s e n (Copel and và T h om a s 1989) do Copelancỉ (1990) sưa đối.người sử d ụ n g xác đ ịn h các hê số vận c hu yể n dựa trẽn các số liệu q u an trắc

2.4 Các yêu cáu vé sỏ liệu đáu vào

Các số liệu dầu vào được phàn thành các n h ó m số liệu hình hoc bùn cat th u\ văn và các lệnh đặc biệl

2.4.1 N h ó m sô liệu hình học

I ) C á c m ă l c á t I i i i a n g ( c á c b a n iilii X 1 X 3 CìR ).

29

2) Các giá trị hệ s ố n h á m M a n n i n g n, các hệ số tổn thất do co hẹp hay m ớ ròng (các bán ghi NC, N V, $ K L , $KI).

3) Các giới han đ áv độ n g (các bán ghi H HD)

4) N ao vét (các bán ghi H, H D , $ D R E D G E $ N O D R E D G E )

5) Diện tích d ò n g c h ả y kh ô ng hiệu q u ả (bản ghi X3)

6) Các giới hạn vận c h u y ế n (bán ghi XL)

2.4.2 N hóm sô liệu bùn cát

]) Các đặc trưng cửa bùn cát (các bản ghi II 12 13 14)

2} Lưu lượng và thành phần hat của bùn cát chãv vào (các bản chi LQ LT LF)

3) T hà n h phan hạt cửa vật liệu bùn cát trẽn đáy sõnti ( ban d l l PF)

a) Các điều kiện biên t hượna lưu:

- Lưu lượng và nhiệt độ nước (các ban ghi Q T )

- Lưu lượng b ùn cát (các bản ghi LQ, LT LF)

hì Các điếu kiên biên hạ Ill'Ll (các ban ghi SRATINCi RC R S).

c ) Các điều kiện nội biên (các bản ghi QT X5 R )

2.5 Kiểm soát đầu ra

H E C - 6 c h o phép neirời sử dụ ng có thè lựa ch on các m ứ c thõ ng Ún dâu ra vó các n h ó m sỏ liệu hình học bùn cát thúy văn với mức độ chi tiết khác nhau Các mức dấu ra này được t ó m tát tron ạ b ản g 2.3 dưới đày

30

B à n g 2.3 Tóm tắt các mức đầu ra của m ô hình HEC-6

thời đoạn dòng chảy, các tham sỗ thủv lực thòns thường cho mỏi mặt cắt ngang

B Đưa ra thông tin về hình học ban đẩu phán hò các tham sỏ iluiv luc qua các lie’ll mặt cát

E Mó tá thông tin một cách chi tiết: các so lieu thúy lực doi với mỏi sỏ gia diện tích, mỏi cao trình thứ n°hiệm tro 112 các tính toán nuớc vạl tại lừng mặt cắt ngang

A Đưa ra thon” tin vẽ the tích hun cát vào và ra khoi mỏ hình, h ill'll

B Đưa ra các thỏna tin vé các ihay dối cao trinh đáy các cao timh mat nước, cao trình đườns sâu nhất Ill'll lượn tỉ bùn cát ra klioi mo hình

c Đau ra chi net bad íióm sứa lai cál, lưu l ư Ợ M i i hun cát va t-ãp plu'1 hạt cua bùn CÚI đáy

D Mò ta thông tin chi tiẽt

SDREDGE A - E Các mức A - E cung cap đâu ra từ các d i ẻ n toan nao vét từ mức lap

lai các số liêu đơn gián (mức A) đến mức mò ta chi nét các thong

A Cung cấp đầu ra đỏi với tất cả các mãt cát ngang tai mức đáu ra cuaCột 8

bủn ghi Q

lựa chọn SPRT

í

SVOl - Đưa ra thòng till VC IIIV tích các tha\ đổi đá\ vu the tích

$VOL

Cót 8

A Đưa ra thóng tin về trọng lượng bùn cát luỹ tích chuyên qua mỏi mặt cắt ngang đối với mỗi lớp kích thước bùn cát

dụng cùng với bản ghi SVOL

32

C h ư ơ n g 3

Ứ N G D Ụ N G M Ô H Ì N H H E C - 6 T Í N H T O Á N V À D ự B Ả O

D I Ễ N B I Ế N Đ O Ạ N S Ô N G H Ổ N G T Ừ H O À B Ì N H Đ Ế N H Ả N Ộ I

3.1 Cơ sở sô liệu

Dựa t heo các yêu cầu về số liệu đẩu vào c ủa m ỏ hình H E C - 6 và sỏ liệu thực tè

do đạc hiện có để tài đã thu thập được các sỏ' liệu sau đáy:

* Sô liệu địa hình: do Đ oà n kháo sát sông H ỏ n g đo đạc Số liệu này hao nôm:

- 38 mặt cất n g a n g (từ sô 10 đên sổ 47) trẽn s ỏ ne Đ à đoan từ Hoà Bình đến ngã ba T h a o - Đ à đo n ă m 1992 và 1997

- 17 mặt cắt n g a n g (từ số 48 đên sô 64) trên s õng H ô n g đoan từ ngã ha T hao

Đà đến Hà Nội đo n ă m 1992 và 1997

- 10 mậ t căt n g a ng (từ số 1 đến số 10) trên sòng T hao đoạn từ Yên Bái đen ngã ba Th a o - Đ à đ o năm 1992

- 6 mặt cắt n g a n g (từ sổ 1 đến số 6) trẽn sóng Lô đoạn từ Vụ Qimnú đốn nszã

ba Lô - H ồ n g đo n ăm 1992

- Trác diện dọc đáy sông từ Ho à Bình đến Hà Noi đo n ăm 1992 và 1997

* S ỏ l i ệ u t h ú y v ã n : c u a T r u n g t á m K h í t ư ợ n g T h ủ v v a n ọ LI ố c ilia, s ỏ l i ệ u n a )

bao gôm:

- Lưu lượng nước và nhiệt độ nước bình quân n g àv thực đo tại các tram Hoa Bình Yên Bái Vụ Q u a n g và Thượng Cát thời kỳ từ sau khi có hổ Hoa Binh (lừ 1986 đến 2003),

- Mực nước bình q u á n ngày thực đo lại hai trạm Sun Táy và Ila Nội thời ky tư sau khi có hò H oà Rình (từ 1986 đen 2003)

* Số liệu bùn cát:

- Đ ộ đục nước s ông bình q uà n ngày thực đo tại các trạm H o à Binh Yen Bái

Vụ Q u a n g và T h ư ợ n g Cát thời kỳ từ sau khi có hổ H ò a Bình (từ 1986 đến 2 00 3) cua

Tr un g tâm K h í tượng T h ủ y văn Q u ố c gia

- Số liệu phán tích thành phần hạt bùn cát lơ l ử n s n g ày cu a 137 lãn đo tại tram H oà Bình và 189 lán đo tai tram Yên Bái các n ă m 70-73 cu a phònii MỎI trường thuộc Viện K h o a học T h ú y lợi

- Sô liệu phân tích thành ph ần hạt bùn cát lơ lúng t r u ne hình t h á n e các nám

70 + 74 tại t rạm T h ư ợ n g Cát và trung bình t h á n e tru ng hình nhiều n ăm tại trạm Vu

Quang

- Sô liệu phân tích thành phán hạt vật liệu đáy dọc s ó n c Đa và s ó ne Hóm: đo dạc n ăm 1996 cửa Viên K ho a hoc Th ủ y loi

33

3.2 Sơ đồ tính

Đ o ạ n m ạ n g s ồn g nghiên cứu lựa chọ n là m ạ n e sông H ổ n e từ Hoà Bình-Yẽn Bái-Vu Q u a n g đến Hà NóI -Thương Cát Vì để tài chu yếu n eh i ẽn cứu diẻn hiên l ò nl: sòng hạ lưu đ ập Hoà Bình nén đoạn s ông Đ à từ Ho à Bình đến n e ã ba T h a o - Đ à vùđoạn s ông H ổ n g từ n gã ba T h a o - Đ à đến Hà Nội được coi là d ò n e chính Trẽn đoạndòng ch ín h này có hai n h án h nh ập lưu và mộ t n h á nh phân lưu N h á n h n háp lưu thú' nhất là s ông T h a o tai n g ã ba T h a o - Đ à và n há nh n hậ p lưu thứ hai là sông Lõ tại neã

ba L ô- H ồn g Phân lưu du y nhất là sông Đ u ố n g từ n gã ba H ồ n a - Đ u ố n s Dựa trẽn cơ

sớ số liệu địa hình đã thu t hập được,

đoạn sông n g h ié n cứu trên đã được

sơ đổ hóa n h ư trong hình 3.1 Đế

gián hoá, riêng phân lưu sông Đ u ố n g

được coi là đ i ê m phân lưu cục bộ từ

tham số thủy lực n h ư các hệ sô n há m

Manning của lòng chí nh, c ủa bãi trái

và hài phái; các hệ số tốn thât đo co

• đèn sổ 6) trên s ô ng Lô và k h o á n g cách giũa chúng

Cán cứ vào hình d a n g trên hình vẽ của từng mặ t cát n ga ng, tiến h àn h chia tnỏi mặt cát n g a n g th ành ba phân khu có hệ số n h á m kh ác nhau: bãi trái, lòng chính va hãi phai Các hệ số n h á m M a n n i n g n c ủa lòn SI chí nh và hai hãi s óng được lựa chon

Hệ sô tổn thất do co h ẹp lấy cô định bằng 0.1 và hệ sỏ tốn thất do m ớ rộng láy

Các th am sỏ vận c h u yê n bùn 'cát: Chọn sỏ lần tính toán lai thành phán vật liêu trên

đáy bàn g 15: chọ n p hư ơn g pháp tính ứng suất tiếp đáy hàng yDS đối với xói và hỏi sét hoặc bùn; c h ọ n ph ương pháp tính vận c hu y ển bùn cát dính là phirơna pháp Krone nếu là bổi và p hư ơn g p háp Ariathurai nếu là xói; các th a m số khác láy theo mặ c định trong HEC-6

C ác qu a n hệ íỊiữu Ill'll lưựiiỊỉ nước và Ill'll híơnạ bùn cút lõiií’ CÕIIÍ’: Vì khônII có số

liệu đo dạc bùn cát đáy nén lưu lượng bùn cát tổng cộ ng được lây theo kinh niihiein bằng 1.25 lần lưu lượng bùn cát lơ lứng Q ua n hệ giữa lưu lượna nước và lưu lượn id hùn cát lổng c ộ n g c há y vào đoạn sông nghiên cứu q u a mật cát xa nhất phía thượnn lưu (mặt cất tram Ho à Bình), quan hệ giữa lưu lượng nước và lưu lượng hùn cát chay vào đo ạn s ô n g n gh iê n cứu của nhánh nháp lưu sông T h a o (tram Yen Bái) L'ĩint> như

hệ giữa lưu lượng nước và lưu lượn2 bùn cát chay vào đoạn sóng nuhiOn cứu cua nhánh n h ập lưu s ô ng Lõ (trạm Vụ Q u a n g ) được xây dưn g dưới cỉạnti đướim conn quan hệ giữa lưu lượng bùn cát tổng cộng Q, ( t á n / n s à y ) và lưu lượmỉ nước Q (ft / s ) Các đường con g q u a n hệ này được xây dựng theo số liệu trung bình ngày trung bình nhiều năm thời kỳ từ sau khi có hổ Hò a Bình (từ n ă m 1986 » C h ú n e đéu có dạn tỉ hàm m ũ khá chặt chẽ với hệ số tương qu an R khá lớn (trên 0.95) Kết quá cu thó dược t hõng kẽ trong báng 3.1 và được thế hiện trẽn các hình 3.2; 3.3 và 3.4

B àng 3.1

P h ư ơ n g trìn h hối q u y và hệ sô tư ơ n g của các quan hệ Qs - Q tại các trạm

T h à n h p lú in kích thước hạt của lưu lượng b ùn cát t ónx cộng: Sò liệu phán tích

thành phần hạt h ùn cát lơ lửng của 137 lán đo tại t rạm H o à Bình và 189 lán đo tại trạm Yen Bái các n ă m 70-73 của phòng Mỏi trường t h uộ c Vi ện K hoa hoc T h ú y lợi

vù số liệu phân tích thà nh phần hạt hùn cát lơ lừng t r u n a bình t h á n e Irunc binh nhiẽu năm tại trạm Vụ Q u a n g được sử d u n e đê lựa c h ọn 9 đ ườ ng c o n e cap phối thành phán

35

kích thước hạt c ù a lưu lượng bùn cát tổng cộ ng tương ứng với 9 c ấp lưu lượng nước

và lưu lượng bùn cát t ổng cộ ng khác nhau tại từng trạm Việc lựa ch ọn dưa trẽn

ngu yên tắc c h ọn c ác đường cong c ấp phối hat có d ạ ng phổ biến và lưu ì ươn íi nước

tương ứng với các lần lấy mẫ u phán tích hạt ấy phai bao trùm các cấp lưu lượng từ nhỏ nhất đến lớn nhất

H ình 3.2 Q uan hệ Qs - Q tại trạm Hoà Binh

(tru n g bình ngáy tru n g bình nhiếu năm th ờ i ký 1986-2002)

H ình 3.3 Q uan hệ Q s • Q tại trạm Yên Bái

(tru n g b in h n gày tru n g binh nhiếu năm th ờ i kỳ 1986-2002)

36

Qs (tấn/ngày)

H ình 3.4 Quan hệ Qs - Q tại trạm Vụ Q uang

(tru n g b in h ngày tru n g binh nhiếu năm th ờ i kỷ 1986-2002)

Từ các đư ờn g c o n s cap phối hai đã lựa chọn, tiến hành phân cáp các lóp kíchihước hat theo hệ t h ống phân chia c ủa Hôi Đ ịa lý M ỹ ( A G U ) như H EC- 6 dã ấn định.Kết q uá phân c ấp các lớp kích thước hạt ứng với 9 cấp lưu lượng nước và Ill'Ll ILi'trnii bùn cát tống c ộ n g tai các train Hòa Bình, Yên Bái Vu Ọ u a n g dược t h o n " ké tron lí các h án e 3.2 3.3 và 3.4

Bang 3.2

Các lớ p kích th ư ớ c hạt theo hệ th ố n g phân chia của Hội Địa lý Mỹ (AGU) ứng vớ i 9 cãp lưu lư ợng nước và lưu lượng bùn cát tô n g cộ n g tạ i trạm Hoa Binh

I Q d ì V s ) 12 l 13 4 2 3 7 8 7 1 3 3 6 1 4 1 9 6 5 1808 1 1 2 3 3 0 7 7 3 2 5 2 4 8 3 5 6 6 7 X 4 4 ] 4 vì Q,

So liẹu p h í m tích tliủnli p h â n hạt vái liệu cỉớỵ được lâv theo số liệu đo đạc

nam 1996 c ủ a Vi ện K h o a hoc T h ú y lợi

3.3.3 Sô liệu thúy văn

Ch ọ n thời đ o ạn tính toán At bằng 1 ngày; q u á trình lưu l ư ợ n2 nước va nhiệt do nước s ô n g đến đ o a n sõng nghiê n cứu lấy quá trình Ill'Ll lươn tỉ nước và nhiệt (In nun'L

3 8