Luận văn Ứng dụng mô hình thủy lực một và hai chiều kết hợp hdm xây dựng bản đồ ngập lụt hạ lưu sông Cái Nha Trang

và chế biến thủy sản, sản xuất vật liệu xây dựng [5]. 2.2.2. Tình hình xã hội * Dân số: Theo số liệu niên giám thống kê tỉnh Khánh Hoà năm 2009 dân số trên lưu vực sông Cái Nha Trang là 552060 người (trong đó nam là 271602 người chiếm 49,2%, nữ là 280458 người chiếm 50,8%). Dân số ở thành thị là 317496 người chiếm 57,5%, nông thôn là 234564 người chiếm 42,5%; nơi có mật độ dân số cao nhất là Thành phố Nha Trang với 1540 người/km2. * Giáo dục : Trên lưu vực sông Cái Nha Trang có 82 trường tiểu học, chiếm 25,6% so với toàn tỉnh; 37 trường trung học cơ sở, chiếm 37%; 14 trường trung học phổ thông, chiếm 50%. Giáo viên tiểu học là 2010, chiếm 42,4% so với toàn tỉnh; 1756 giáo viên trung học cơ sở, chiếm 41,8%; 873 giáo viên trung học phổ thông, chiếm 47,9%. Trên lưu vực có 5 trường đại học, 3 trường cao đẳng, tập trung ở thành phố Nha Trang. * Y tế : Tổng số giường bệnh các huyện, thành phố trên lưu vực sông Cái Nha Trang là 1941, chiếm 64,5% so với toàn tỉnh, có 410 bác sĩ (chiếm 68,8%), 399 y sĩ (chiếm 60,1%), 600 y tá (chiếm 64,2%), 20 dược sĩ cao cấp (chiếm 76,9%), 124 dược sĩ trung cấp (chiếm 76%), 39 dược tá (chiếm 52%) [1].  27 CHƯƠNG 3 CƠ SỞ LÝ THUYẾT MÔ HÌNH HDM 3.1 CƠ SỞ LÝ THUYẾT Đây là mô hình có nguồn gốc từ Hoa Kỳ, được Viện Khoa học Khí tượng Thủy văn & Môi trường Hà Nội cải tiến để dễ sử dụng, phù hợp với điều kiện của Việt Nam và tích hợp thêm các công cụ tính toán. Điển hình là đã ứng dụng thành công cho nhiều đề tài như đề tài cấp nhà nước về diễn toán ngập lụt hạ lưu hồ Hòa Bình khi hệ thống hồ Sơn La - Hòa Bình bị vỡ đập, đề tài cấp tỉnh về xây dựng bản đồ nguy cơ ngập lụt sông Cà Ty và sông La Ngà tỉnh Bình Thuận, xây dựng bản đồ nguy cơ ngập lụt tỉnh Bình Định, bổ sung đặc điểm khí tượng thủy văn Phú Yên và xây dựng bản đồ ngập lụt sông Kỳ Lộ, đề tài cấp bộ về nghiên cứu xây dựng hệ thống phân tích, giám sát, cảnh báo và dự báo lũ, ngập lụt và hạn hán cho hệ thống sông Ba, dự án cấp tỉnh: Lập bản đồ ngập lụt lưu vực sông Dinh Ninh Hòa và sông Cái Nha Trang. Mô hình HDM được PSG.TS Lương Tuấn Anh phát triển từ mô hình DHM, mô hình DHM nguyên bản chỉ là mô hình thủy lực một chiều trong kênh hình chữ nhật đã được phát triển thành mô hình HDM là mô hình một chiều trong sông tự nhiên kết hợp với hai chiều trong bãi ngập. Ngoài ra mô hình HDM còn được thêm các công cụ tính toán và xử lý các công trình tác động đến diễn biến ngập lụt. Khi diễn toán trong sông mô hình là một chiều đối với khu vực tràn từ sông vào bãi ngập mô hình tự động chuyển thành hai chiều. Lượng nước tràn qua công trình được tính theo công thực đập tràn đỉnh rộng. Mô hình HDM (Hydro Dynamic Model) được phát triển từ mô hình DHM (Diffusion Hydrodynamic Model), là sự kết hợp giữa mô hình dòng một chiều trong sông và dòng hai chiều trên bãi sông cơ sở: - Thêm 01 thành phần quán tính vào phương trình chuyển động của dòng chảy trong hệ thống sông tạo ra khả năng mô hình có thể mô phỏng dòng chảy trong vùng có ảnh hưởng của thuỷ triều, nước dâng ven biển. - Bổ sung các điều kiện biên trong bao gồm đập tràn, cống có cửa điều khiển vào hệ phương trình 01 chiều và đập tràn, cống đồng thời đê, đường phía trên, cống  28 phía dưới đối với hệ phương trình 02 chiều mô tả chuyển động của dòng ngập lụt trên bãi sông, cho phép mô hình có thể mô phỏng đê, đường, cống qua đê ... - Mô hình có thể ứng dụng cho lòng sông có hình dạng bất kỳ thay thế mặt cắt hình chữ nhật như trong mô hình DHM. Hệ phương trình sóng động lực Saint - Venant được áp dụng trong tính toán truyền lũ một chiều trong hệ thống sông và xấp xỉ sóng khuyếch tán bằng cách bỏ qua các thành phần gia tốc trong phương trình bảo toàn động lượng trong hệ phương trình 02 chiều mô tả quá trình ngập lụt. Mô hình 1 chiều và 2 chiều được kết nối với nhau dựa trên nguyên tắc liên tục và chảy tràn bờ, tràn qua đê, đường và các công trình khác. Hệ phương trình cơ bản của mô hình có dạng như sau. 3.1.1 Mô hình một chiều trong sông Dòng chảy trong hệ thống sông được mô tả bằng phương trình sóng động lực bao gồm: Phương trình liên tục: q t A x Q =∂ ∂+∂ ∂ (3.1) Phương trình động lực: 0)( 2 =+∂ ∂+∂ ∂+∂ ∂ fSx Zg A Q xt Q (3.2) Độ dốc ma sát được xác định theo phương trình Manning: 3/42 2 f RA QQn S = Giải hệ phương trình trên có thể được dựa trên cơ sở quan niệm đoạn sông như một thể tích hữu hạn và xấp xỉ phương trình liên tục, phương trình động lượng bằng các quan hệ: * Quan hệ dạng sông Áp dụng cách xử lý tương tự cách xử lý của GS. Nguyễn Ân Niên trong mô hình KOD bằng việc xấp xỉ phương trình liên tục và động lượng dưới dạng: dt dZFQQ ii .1 =− −  29 ),( 1 iii ZZfQ −= Hình 3. Mô hình 1 chiều trong hệ thống sông Trong phương trình (3.2) gia tốc đối lưu có trị số nhỏ nhất và không đáng kể trong trường hợp dòng chảy trong kênh dẫn hoặc sông tự nhiên có hình dạng mặt cắt ít thay đổi. Đối với các kênh dẫn thành phần này có thể bỏ qua và phương trình động lượng có dạng: 03/42 2 =+∂ ∂+∂ ∂ RA QQn gA x ZgA t Q Áp dụng sơ đồ sai phân hiện tăng cường (enhanced explicit): 0 )1(3/41 121 1 1 1 1 =+∆ −+∆ − ++ ++ − + + + jj jj i J i j ij j i j i RA QQn g x ZZAg t QQ Đặt x ZZGRAD j i j i ∆ −= + − + 1 1 1 ; Ta có: )1/()1(),( )1(3/41 2 11 1 11 ++ ++ − ++ +∆−∆== jj j jj i j i j i j i RA Qn g t GRADAgQ t ZZfQ * Quan hệ dạng đê, đập tràn - Trường hợp chảy tự do: 1 1 1 *3/2 +− + j i j i hh p 1 1 1 1 11 2**35.0),( +− + − ++ == jijijiji ghBZZfQ  30 - Trường hợp chảy ngập: 1 1 1 *3/2 +− + j i j i hh f )(2***91.0),( 111 11 1 11 ++ − ++ − ++ −== jijijijijiji hhgBhZZfQ Trong đó: o j i j i ZZh −= +−+− 1111 và oZ là cao độ ngưỡng tràn * Quan hệ cống có cửa điều khiển - Trường hợp 1 1 1 + − + j i j s hh f và 11 ++ j i j s hh f (chảy kiểu đập tràn): phương pháp tính như trường hợp (b) - Trường hợp 1 1 1 + − + j i j s hh p và 11 ++ j i j s hh p (chảy kiểu cống hoàn toàn): )(**2***7.0 111 11 ++ − ++ −= jijijsji hhghBQ - Trường hợp 1 1 11 + − ++ j i j s j i hhh pp (trường hợp chảy nửa cống): 21 1 QQQ ji +=+ [ ]2/31112/31111 )()(****2*3/2 ++−++− −−−= jsjijiji hhhhBgQ µ )(**2*** 111 1 2 ++ − + −= jijiji hhghBQ µ ; .61.0=µ Ký hiệu h, độ mở của cống. 3.1.2 Mô hình hai chiều trong vùng ngập Dòng chảy hai chiều trong vùng bãi ngập được mô tả bằng phương trình sóng khuyếch tán hai chiều: - Phương trình liên tục: q t h y vh x uh =∂ ∂+∂ ∂+∂ ∂ (3.3) - Phương trình động lượng theo phương x và theo phương y: 0)S x Z(g y uv x uu t u fx =+∂ ∂+∂ ∂+∂ ∂+∂ ∂ (3.4) 0)S x Z(g y vv x vu t v fy =+∂ ∂+∂ ∂+∂ ∂+∂ ∂ (3.5) Khi bỏ qua thành phần gia tốc theo xấp xỉ của sóng khuyếch tán thì phương trình động lượng theo phương Ox và theo phương Oy trở thành:  31 fxx0 SSx h −=∂ ∂ (3.6) fyy0 SSy h −=∂ ∂ (3.7) - Độ dốc ma sát được xác định theo phương trình Manning theo phương Ox và theo phương Oy: 3/42 2 fx RA QQn S = và 3/42 2 fy RA QQn S = (3.8) Trong đó: h = độ sâu. z = cao độ mực nước. vx và vy = các thành phần vận tốc theo phương Ox và Oy. t = thời gian. g = gia tốc trọng trường. Sox và Soy = độ dốc đáy theo phương Ox và Oy. Sfx và Sfy = thành phần độ dốc ma sát theo phương Ox và Oy. nfx và nfy = hệ số Manning theo phương Ox và Oy. 3.1.3 Ghép nối mô hình một chiều và hai chiều Kết nối giữa mô hình dòng chảy một chiều trong hệ thống sông và mô hình dòng chảy hai chiều trong vùng ngập. Mô hình một chiều dòng chảy trong sông và mô hình hai chiều dòng chảy trong vùng ngập được liên kết qua sự trao đổi dòng chảy vùng ngập và dòng chảy Hình 4. Liên kết mô hình 1 chiều và 2 chiều  32 trong sông. Dựa trên cơ sở nguyên tắc liên tục và trao đổi dòng chảy trong sông và vùng ngập lụt thông qua điều kiện biên trong dưới dạng công trình như bờ sông, đê bao, đập tràn trên cơ sở lý thuyết được thể hiện trong công trình nghiên cứu. Mô hình 1 chiều dòng chảy trong sông và mô hình chảy tràn được kết nối dựa trên phương trình trao đổi nước giữa sông và ô lưới. tOIVV tttl t l ∆−+= +++ )( 111 Trong đó: Vl: thể tích nước trong ô tại thời điểm t+1 hoặc t, phụ thuộc vào mực nước. I: dòng chảy từ đoạn sông vào ô lưới gần kề. O: dòng chảy từ ô lưới ra đoạn sông gần kề. Dựa theo công thức đập tràn đỉnh rộng, dòng chảy vào hoặc ra ô lưới được tính theo công thức sau: LhhscI wrbf ∆−= 2/3)( nếu wr hh > và fpr hh > LhhscO wfpbf ∆−= 2/3)( nếu wfp hh > và rfp hh > Trong đó: cf: hệ số lưu lượng. hfp: mực nước trong khu ngập. ∆L: chiều dài đoạn sông. hr: mực nước trong sông. hw: cao trình bờ. sb: hệ số dòng chảy ngập. 3.2 YÊU CẦU CỦA MÔ HÌNH 3.2.1 Yêu cầu về mô phỏng Miền tính được mô phỏng là các ô lưới hình vuông có kích thước như nhau, mỗi ô lưới được đánh số thứ tự từ trái sang phải, từ trên xuống dưới. Ô lưới chia càng nhỏ thì độ phân giải càng cao, mô phỏng càng sát với thực tế. Tuy nhiên ô lưới càng nhỏ thì sô lượng ô lưới càng lớn khi đó số lượng phép tính quá lớn dẫn đến tình trạng tràn bộ nhớ trong quá trình tính toán.  33 Các công trình được mô phỏng theo biên của các ô, sông suối được mô phỏng theo các đường gấp khúc đi qua trung tâm của các ô. Vì vậy ô lưới càng nhỏ thì mô phỏng công trình và sông suối càng chính xác. Chiều rộng của sông lớn nhất không được lớn hơn kích thước của ô nên cần xác định kích thước ô lưới phù hợp để đảm bảo số lượng ô lưới không quá nhiều, sông suối nằm trọn trong ô lưới và có độ phân giải phù hợp để đảm bảo độ chính xác cho mô phỏng công trình, mô phỏng sông. Hệ thống ô lưới được bao hết cả vùng có khả năng ngập dùng để diễn toán ngập lụt. Các ô lưới được tạo liên kết hệ thống với các ô lưới xung quanh thông qua số thứ tự của ô lưới đã được đánh dấu. Các công trình như đường giao thông, đê được tạo liên kết công trình để xác định lượng nước tràn bờ giữa các ô lưới, đánh giá tác động của công trình đến diễn toán ngập lụt. Cao độ ô lưới được tính toán từ bản đồ địa hình, mỗi ô lưới được xác định một trị số cao độ nền của ô được tính từ cao độ nền địa hình. Hệ số nhám ô lưới được xác định vào điều kiện lớp phủ và thổ nhưỡng bề mặt của vùng diễn toán ngập lụt. Ngoài tạo liên kết hệ thống ô lưới còn phải tạo liên kết sông với ô lưới mà sông chảy qua, đối với các ô lưới này hệ số nhám được xác định là hệ số nhám trung bình của lòng sông và địa hình được xác định trên cơ sở số liệu mặt cắt ngang. Tổng cao trình đáy sông và độ sâu bằng cao trình nền của ô lưới có sông chảy qua nên sau sau khi tính toán các đặc trưng của mặt cắt thì cần phải kiểm tra lại. 3.2.2 Yêu cầu về số liệu Địa hình nền ảnh hưởng rất lớn đến diễn toán ngập lụt nên số liệu địa hình sử dụng trong mô hình càng chi tiết càng tốt. Bản đồ địa hình có thể sử dụng là bản đồ giấy, tuy nhiên loại bản đồ này rất bất tiện khi sử dụng. Để thuận tiện tính toán số liệu địa hình thì bản đồ địa hình cần được số hóa và có thuộc tính cao độ của đường bình đồ và điểm cao độ. Tạo thêm các lớp về lưới tính và số ô lưới trên bản đồ số để dễ dàng biên tập, xử lý và tính toán. Mỗi ô lưới nên có ít nhất một điểm cao độ và giá trị cao độ nền ô lưới được tính trung bình của các điểm trong ô, các ô lưới không có điểm cao độ thì được nội suy từ các ô lưới lân cận.  34 Số liệu mặt cắt đo đạc được dẫn cao độ với hệ cao độ và tọa độ thống nhất với bản đồ địa hình nền. Số liệu mặt cắt sông càng nhiều thì càng tốt, tuy nhiên vì điều kiện kinh tế không thể đo mặt cắt ngang nhiều vì vậy để đảm bảo độ chính xác của trong quá trình diễn toán ngập lụt thì cần đo mặt cắt tại những vị sông cong, sông co thắt hoặc mở rộng đột ngột, những đoạn sông có thay đổi nhiều về các yếu tố thủy lực. Mỗi ô lưới có sông chảy qua đều phải tính toán mặt cắt, tuy nhiên không phải ô lưới nào cũng có mặt cắt được đo vì vậy các ô lưới không được đo mặt cắt thì được nội suy và ngoại suy từ số liệu mặt cắt thực đo. Các số liệu mặt cắt nội suy và ngoại suy sau khi tính toán thì được kiểm tra tính hợp lý với cao độ của bản đồ địa hình nền. Số liệu Khí tượng Thủy văn dùng để kiểm định mô hình là số liệu của một trận mưa lũ quan trắc đồng bộ. Để mô phỏng đúng với thực tế, ngoài các trạm Thủy văn đang sử dụng cần phải bổ sung thêm các trạm mới, khống chế tất cả các nhánh sông gây ngập lụt cho hạ lưu. Các số liệu quan trắc đều được chỉnh biên tại Đài Khí tượng Thủy văn khu vực Nam Trung Bộ, đảm bảo số liệu sử dụng là chính xác, theo đúng quy trình quy phạm của ngành Khí tượng Thủy văn. Các số liệu đo đồng bộ trên lưu vực sông Cái Nha Trang gồm lưu lượng của sông Thác Ngựa, sông Chò, sông Suối Dầu, số liệu thủy triều tại cầu Bình Tân, cầu Trần Phú, số liệu mực nước kiểm tra tại 5 trạm mực nước từ biên trên xuống biên dưới. Các trạm mực nước kiểm tra được bố trí trên, dưới Quốc lộ 1 và đường sắt, bên trái và phải đường 23/10 vì đây là các công trình ảnh hưởng nhiều đến diễn toán ngập lụt ở hạ lưu. Số liệu khí tượng, thủy văn, hải văn dùng để chạy cho các kịch bản ngập được thống kê và tính theo các tần suất thiết kế là 1%, 3%, 5%, 10%. Các số liệu thống kê để tính tần suất càng dài thì càng tốt, đối với lưu vực sông Cái Nha Trang sử dụng số liệu thống kê của trạm Khí tượng Nha Trang, trạm Thủy văn Đồng Trăng và điểm đo mưa Khánh Vĩnh. Vì sông đổ ra biển nên biên dưới của mô hình cần phải tính toán số liệu mực nước triều ứng với các kịch bản ngập. Số liệu triều của các mô hình dự báo đều tính theo “0” Hải đồ nên phải chuyển cao độ số liệu mực nước triều cùng hệ cao độ với bản đồ địa hình nền.  35 CHƯƠNG 4 ỨNG DỤNG MÔ HÌNH HDM XÂY DỰNG BẢN ĐỒ NGẬP LỤT HẠ LƯU SÔNG CÁI NHA TRANG 4.1 TÍNH TOÁN SỐ LIỆU ĐẦU VÀO 4.1.1 Tính toán số liệu địa hình * Địa hình Dữ liệu địa hình là dữ liệu rất quan trọng của mô hình thủy lực, số lượng tính toán rất lớn, đòi hỏi độ chính xác cao. Dữ liệu địa hình của lưu vực sông Cái Nha Trang được tính toán từ bản đồ tỷ lệ 1/10.000 hệ tọa độ VN 2000 do Bộ Tài nguyên Môi trường xây dựng năm 2006 bằng công nghệ hàng không và được lưu ở dạng Microstation. Trên lưu vực sông Cái Nha Trang thu thập bản đồ dùng để tính số liệu địa hình gồm có 46 xã, phường, thị trấn, tập trung ở huyện Diên Khánh, Cam Lâm và thành phố Nha Trang. Bản đồ được lưu thành các mảnh theo các xã riêng biệt, xã ít nhất có 1 mảnh bản đồ, xã nhiều nhất có 12 mảnh bản đồ. Để có dữ liệu địa hình cho lưu vực sông Cái Nha Trang thì cần phải chuyển đổi và ghép các mảnh bản đồ thành bản đồ lưu vực hoàn chỉnh bao gồm các lớp thông tin là đường bình đồ, giao thông, sông suối, điểm cao độ và tên địa danh. Đã dùng phần mềm MapInfo 11.0 để chuyển đổi bản đồ từ Microstation sang MapInfo, hệ quy chiếu và hệ tọa độ VN 2000 - Khánh Hòa múi 3. Cao độ địa hình đầu vào của mô hình là cao độ của các ô lưới, các cao độ này được tính từ bản đồ địa hình trên. Từ các đường bình đồ và các điểm cao độ trên bản đồ tính toán được cao độ của các ô lưới. Các điểm cao độ và các đường bình đồ đều chưa được cấy các giá trị, sau khi chuyển đổi từ bản đồ gốc ở dạng Microstation bản đồ mới chỉ ở dạng số, chưa có thuộc tính về cao độ. Để tính toán cao độ ô lưới được dễ dàng và xử lý kết quả chạy mô hình sau này thì cần phải gán cao độ cho các đường bình đồ và điểm cao độ. Các đường bình đồ sau khi gán cao độ dược chuyển thành các điểm bằng phần mềm Vertical Maper 3.0 dùng để tính toán số liệu đầu vào và xử lý kết quả đầu ra của mô hình. Các điểm cao độ được cập nhật tọa độ tự động bằng phần mềm MapInfo 11.0.  36 Ngoài ra do bản đồ địa hình được xây dựng từ năm 2006, sau 7 năm đã có sự thay đổi ở một số vùng phía tây thành phố Nha Trang do xây dựng các công trình. Việc đô thị hóa của thành phố đã xây dựng nhiều khu đô thị, kéo theo đó là việc san lấp mặt bằng để xây dựng công trình đã ảnh hưởng đến ngập cục bộ của các vùng này. Các khu đô thị được cập nhật thay đổi cao độ nền so với trước đây là khu đô thị Phước Long, Mỹ Gia, Vĩnh Điềm Trung, Nam Sông Cái, tăng cao độ nền thêm từ 1,5 - 2,0m. Những thay đổi cao độ địa hình tại các khu đô thị được cập nhật bổ sung trên nền bản đồ địa hình và sử dụng để tính toán số liệu đầu vào cho mô hình. Hình 5. Miền tính mô hình HDM * Mặt cắt Các mặt cắt được đo tại vị trí sông cong, lòng sông bị biến đổi nhiều, có sự thay đổi về thủy lực. Trên lưu vực sông Cái Nha Trang được bố trí 32 mặt cắt ngang, trong đó 3 mặt cắt trùng với trạm đo lưu lượng, 4 mặt cắt trùng với trạm đo  37 mực nước kiểm tra và 2 mặt cắt tại cửa ra trùng với trạm đo thủy triều. Mặt cắt được dẫn cao độ Quốc gia, đã được kiểm tra thẩm định và kế thừa từ dự án “Lập bản đồ lưu vực sông Dinh Ninh Hòa và sông Cái Nha Trang” do Đài Khí tượng Thủy văn khu vực Nam Trung Bộ thực hiện năm 2010. Các mặt cắt ngang sông được bố trí trên sông Thác Ngựa, sông Chò, sông Cái, sông Suối Dầu và sông Quán Trường trên lưu vực sông Cái Nha Trang. Ngoài ra còn xác định thêm các mặt cắt ở mương Bầu Sấu, sông Tháo và sông Tắc để mô phỏng phù hợp với thực tế. Ngập lụt ở vùng hạ lưu do sông Cái Nha Trang và sông Suối Dầu gây ra vì vậy mặt cắt của trên các sông này có vai trò rất quan trọng trong mô hình, đối với mương Bầu Sấu, sông Tháo và sông Tắc không gây ngập lụt nhưng lại ảnh hưởng rất nhiều đến quá trình thoát lũ. Hình 6. Sơ đồ thủy lực mô hình HDM Các số liệu mặt cắt đo từ năm 2010, đến nay tại một số đoạn sông đã có sự thay đổi của việc xây dựng các công trình thủy lợi. Xác định vị trí thay đổi do xây dựng kè ở khu vực thị trấn Diên Khánh và thành phố Nha Trang để đo lại mặt cắt ngang. Vùng hạ lưu sông Cái Nha Trang có 03 mặt cắt được đo lại và đo thêm 02 mặt cắt vào giữa năm 2013 để kịp thời cập nhật sự thay đổi. Ngoài ra do dự án chỉnh trị sông Tắc và sông Quán Trường đã xây dựng kè 2 bờ bằng đá hộc với mặt  38 cắt được mở rộng gấp khoảng 3 lần trước đây và cao trình bờ sông được nâng cao gấp khoảng 2 lần. Sông được chỉnh trị ngắn và thẳng hơn trước đã tăng khả năng thoát lũ cho phía tây thành phố Nha Trang. Những thay đổi do chỉnh trị sông Tắc và sông Quán Trường được cập nhật dựa trên hồ sơ thiết kế thi công và bản đồ quy hoạch của dự án này (hình 42 - phụ lục). Hình 7. Tính toán các yếu tố mặt cắt ngang Hình 8. Phân bố các mặt cắt ngang trên hệ thống sông Cái Nha Trang Mỗi mặt cắt ngang được tính toán các yếu tố về độ sâu, cao trình đáy sông, chiều ngang sông ứng với các mức độ sâu khác nhau, hệ số nhám lòng sông. Mặt  39 cắt ngang tính cho từng ô lưới có sông chảy qua, đối với ô lưới không được đo mặt cắt thì nội suy và ngoại suy từ các ô lưới có số liệu đo. * Số liệu công trình Các công trình giao thông thủy lợi ảnh hưởng đến diễn biến ngập lụt ở hạ lưu vì vậy cần phải nhập dữ liệu công trình vào mô phỏng ngập lụt. Dữ liệu công trình bao gồm cao độ của đường giao thông, đập tràn, thông số của cống ngầm. Các dữ liệu này được đo đạc ở ngoài thực địa, hồ sơ kỹ thuật các công trình và trên bản đồ số hóa. Các công trình được mô phỏng theo đường biên của các ô lưới, cao độ công trình sẽ quyết định đến việc lượng nước tràn từ ô này sang ô khác. Cao độ của một đường giao thông trên thực tế ở các đoạn đường khác nhau không cùng một trị số, số điểm cao độ của đường giao thông trên bản đồ là không nhiều, chính vì vậy cần phải bổ sung thêm cao độ từ các hồ sơ kỹ thuật xây dựng và điều tra khảo sát thực địa. Các công trình ảnh hưởng đến quá trình ngập lụt hạ lưu sông Cái Nha Trang bao gồm đường Quốc lộ 1A, đường sắt, đường 23 tháng 10. Các đường liên huyện, liên xã và đường làng ít ảnh hưởng, có những đoạn đường không ảnh hưởng đến quá trình ngập lụt. Những đường và đoạn đường không ảnh hưởng đến ngập lụt là những đường và đoạn đường có cao trình thấp hơn cao độ vùng đất hai bên đường. Các con đường ảnh hưởng nhiều đến ngập lụt hạ lưu sông Cái Nha Trang có cao trình từ 4,5 đến 5,0 m và cao hơn vùng đất hai bên đường từ 0,5 đến 2,5 m. Đường Quốc lộ 1A và đường sắt chạy theo hướng bắc - nam nên ảnh hưởng rất nhiều đến thoát lũ hạ lưu, gây ngập lụt sâu về phía thượng lưu các con đường này. Đối với đường 23 tháng 10 ảnh hưởng đến việc trao đổi nước giữa sông Cái Nha Trang và sông Quán Trường. Trên các con đường này có nhiều cầu cống bắc qua các kênh, mương và sông nhỏ, đáng chú nhất là cầu Lùng trên Quốc lộ 1A không bắc qua sông mà qua một cái hồ nhỏ và có ảnh hưởng rất nhiều đến quá trình thoát lũ hạ lưu sông Cái Nha Trang. Khi lũ trên sông Cái lớn hơn sông Suối Dầu thì lượng nước trên sông Cái tràn chảy ngược vào sông Suối Dầu khiến mực nước phía tây nam thị trấn Diên Khánh dâng cao lúc đó sẽ có lượng nước đáng kể qua cầu Lùng xuống phía hạ lưu đường Quốc lộ 1A và nhập vào sông Quán Trường. Các cầu lớn khác ít  40 ảnh hưởng đến quá trình ngập lụt hạ lưu sông Cái Nha Trang như cầu Phú Cốc, cầu Sông Cái, cầu Suối Cát, cầu Hà Dừa, cầu Vĩnh Phương, cầu Bình Tân và cầu Trần Phú vì các cầu này có các mố cầu nằm trên hai bờ, các yếu tố thủy lực ảnh hưởng không đáng kể đến lưu lượng chảy qua chân cầu. Ngoài cầu Lùng còn có cầu khác ảnh hưởng nhiều đến ngập lụt hạ lưu sông Cái Nha Trang đó là cầu tràn sông Cái, làm cản trở dòng chảy. Hạ lưu sông Cái Nha Trang có hai công trình thủy lợi ảnh hưởng đến ngập lụt đó là đập ngăn mặn Vĩnh Phương và đập Suối Cát, đây là hai đập tràn. Đập Suối Cát nằm trên sông Suối Dầu ảnh hưởng đến phân phối lưu lượng giữa sông Suối Dầu và sông Quán Trường, đập ngăn mặn Vĩnh Phương làm chậm quá trình thoát lũ trên sông Cái. Tuy nhiên trong quá trình khảo sát đã đặt mặt cắt tại vị trí các tuyến đập nên mức độ ảnh hưởng của các đập cũng đã có trong tài liệu mặt cắt ngang, mặt khác trên đập Suối Cát có nhiều cửa mở khi lũ về và đập Vĩnh Phương được làm bằng đá xếp, nước chảy được qua thân đập. + Mô phỏng công trình hiện trạng Hình 9. Mô phỏng sông suối và công trình trên miền tính Để xây dựng dữ liệu công trình cần phải xác định cao độ của từng đoạn đường giao thông và mô phỏng trên hệ thống lưới tính. Các công trình được mô Mô phỏng sông suối Mô phỏng công trình giao thông  41 phỏng theo biên của các ô lưới, liên kết công trình được xác định giữa các ô lưới liền kề nhau theo quy luật ô lưới đó với ô lưới phía trước, phía phải, phía sau và phía trái. Quá trình liên kết công trình được thực hiên ở tất cả các ô lưới, nếu giữa các ô không có liên kết công trình thì cao độ công trình được xác định là “0”, nếu có liên kết công trình thì cao độ được xác định bằng cao độ của đoạn công trình đó. Tất cả các cao độ này được quy đổi thống nhất sang cao độ Quốc gia. Giữa hai ô lưới vừa có liên kết công trình, vừa có sông chảy qua thì cần phải xác định đặc điểm của công trình bắc qua sông. Nếu đó là cống hoặc cầu làm ảnh hưởng đến mặt cắt ngang sông thì sẽ gây ảnh hưởng nhiều đến diễn biến dòng chảy trong sông, nếu là cầu không ảnh hưởng đến mặt cắt ngang sông thì không mô phỏng, lúc đó dòng chảy qua cầu như dòng chảy trong sông tự nhiên. Khác với mô phỏng công trình, sông được mô phỏng theo trung tâm các ô lưới, các nhánh sông cách nhau ít nhất một ô lưới, trừ các ô lưới chéo nhau. + Chia lưới ô vuông Sau khi chuyển đổi định dạng bản đồ và hệ tọa độ, tiến hành chia lưới ô vuông. Đối với sông Cái Nha Trang đã chia lưới với kích thước 300 x 300m khống chế toàn bộ vùng có khả năng ngập bao gồm 2250 ô lưới. Việc chia ô lưới được căn cứ trên cơ sở xác định sơ bộ vùng có khả năng ngập để xác định phạm vi vùng chia ô lưới. Ngoài ra việc chia ô lưới còn căn cứ vào yêu cầu mô phỏng của mô hình để xác định kích thước ô lưới cho phù hợp. Sử dụng bản đồ địa hình để tính toán cao độ các ô lưới, đánh số ô lưới, cập nhật số ô và kinh vĩ độ của trọng tâm các ô. Các ô liên kết với nhau theo một trật tự nhất định, để tạo dữ liệu đầu vào cho mô hình cần phải tạo liên kết hệ thống. Việc tạo liên kết hệ thống theo cách thủ công mất rất nhiều thời gian và dễ bị nhầm lẫn, để khắc phục dự án đã mô phỏng lưới ô vuông theo ma trận hai chiều và lập chương trình để tạo liên kết bằng ngôn ngữ lập trình Fortran. Mỗi ô lưới được tạo liên kết với ô lưới phía trước, phía phải, phía sau và phía trái, đối với các ô ở biên thì phía nào giáp biên thì được liên kết với ô số “0”. Việc đánh số các ô theo một quy luật nhất định, các số sử dụng là số tự nhiên được bắt đầu bằng 1. Chia ô lưới là bước  42 đầu tiên của việc tạo số liệu địa hình cho mô hình, ảnh hưởng tới tất cả các quá trình tính toán sau này. Các thông tin cần xác định với một ô lưới bao gồm: liên kết hệ thống, liên kết công trình, cao độ ô lưới, cao độ mốc so với mặt chuẩn, kinh độ, vĩ độ và hệ số nhám. Để xác định các thông tin này cần đến sự hỗ trợ của phần mềm MapInfo 11.0 và ng