BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN VĂN GIÁP NGHIÊN CỨU CƠ CHẾ BỒI, XÓI VÀ ĐỀ XUẤT GIẢI PHÁP CHỈNH TRỊ KHU VỰC CỬA LẤP TỈNH BÀ RỊA - VŨNG TÀU Chuyên ngành: Xây dựng công trình thủy Mã số: 62-58-40-01 NGƯỜI HƯỚNG DẪN KHOA HỌC: 1. PGS.TS. Nguyễn Trung Việt 2. PGS.TS. Trương Văn Bốn HÀ NỘI, NĂM 2015 LỜI CAM ĐOAN Tác giả xin cam đoan công trình nghiên cứu tác giả. Các kết nghiên cứu kết luận luận án trung thực, không chép từ nguồn hình thức nào. Việc tham khảo nguồn tài liệu thực trích dẫn ghi nguồn tài liệu tham khảo theo quy định. Tác giả luận án Nguyễn Văn Giáp i LỜI CẢM ƠN Để hoàn thành luận án này, tác giả nhận giúp đỡ, hướng dẫn chu đáo, tận tình giáo viên hướng dẫn, quan tâm, sát cánh gia đình, quan đồng nghiệp. Đặc biệt, tác giả nhận nhiều điều kiện thuận lợi từ Ban Giám hiệu nhà trường, Phòng Đào tạo ĐH&SĐH, Khoa Công trình, Bộ môn Thủy công - Khoa Công trình trình học tập nghiên cứu. Tác giả xin trân trọng gửi lời cảm ơn sâu sắc tới PGS.TS Nguyễn Trung Việt, PGS.TS Trương Văn Bốn quan tâm hướng dẫn, giúp đỡ để tác giả hoàn thành luận án. Tác giả xin bày tỏ lòng biết ơn tới quan, Ban Giám hiệu nhà trường, Phòng Đào tạo ĐH&SĐH, Khoa Công trình, Bộ môn Thủy công - Khoa Công trình quan tâm giúp đỡ, tạo điều kiện để tác giả hoàn thành luận án này. Tác giả xin bày tỏ lòng biết ơn tới nhà khoa học quan tâm chia sẻ, góp ý bổ sung cho tác giả nhiều thông tin bổ ích. Cuối tác giả xin trân trọng cảm ơn gia đình, đồng nghiệp sát cánh động viên tác giả vượt qua khó khăn thực luận án. Tác giả luận án ii MỤC LỤC MỞ ĐẦU .1 1. Tính cấp thiết luận án 2. Mục tiêu nghiên cứu 3. Đối tượng phạm vi nghiên cứu .2 4. Nội dung nghiên cứu .2 5. Cách tiếp cận phương pháp nghiên cứu 6. Ý nghĩa khoa học thực tiễn .4 7. Cấu trúc luận án CHƯƠNG TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU VỀ CỬA SÔNG .7 1.1 Tổng quan chung nghiên cứu cửa sông .7 1.2 Tổng quan nghiên cứu phân loại cửa sông .9 1.2.1 Các nghiên cứu phân loại cửa sông giới .9 1.2.2 Các nghiên cứu phân loại cửa sông Việt Nam 1.2.3 Các nghiên cứu phân loại cửa sông Cửa Lấp .11 1.3 Tổng quan nghiên cứu động lực học cửa sông .11 1.3.1 Các nghiên cứu động lực học cửa sông giới .11 1.3.2 Các nghiên cứu động lực học cửa sông Việt Nam 13 1.3.3 Các nghiên cứu động lực học cửa sông Cửa Lấp .15 1.4 Tổng quan nghiên cứu công trình chỉnh trị cửa sông 17 1.4.1 Các nghiên cứu công trình chỉnh trị cửa sông giới .17 1.4.2 Các nghiên cứu công trình chỉnh trị cửa sông Việt Nam 18 1.4.3 Các nghiên cứu công trình chỉnh trị cửa sông Cửa Lấp 19 1.5 Tổng quan phương pháp nghiên cứu cửa sông .20 1.5.1 Đặc điểm phương pháp nghiên cứu cửa sông .20 1.5.2 Các phương pháp nghiên cứu truyền thống 21 1.5.3 Các phương pháp nghiên cứu đại 22 1.6 Nhận xét chung thành tựu nghiên cứu cửa sông .23 1.7 Các tồn nghiên cứu cửa sông 23 1.7.1 Các tồn chung .24 1.7.2 Những tồn nghiên cứu Cửa Lấp .25 iii 1.8 Đặt vấn đề nghiên cứu luận án 26 1.9 Kết luận chương 26 CHƯƠNG CƠ SỞ LỰA CHỌN MÔ HÌNH THỦY THẠCH ĐỘNG LỰC CHO KHU VỰC CỬA LẤP .28 2.1 Đặt vấn đề 28 2.2 Cơ sở lý luận mô hình tính toán thủy động lực học hình thái cửa sông ven biển .28 2.3 Lựa chọn mô hình tính toán cho khu vực nghiên cứu .30 2.3.1 Giới thiệu mô hình số đại .30 2.3.2 Lựa chọn mô hình tính toán cho khu vực nghiên cứu .34 2.4 Cơ sở lý thuyết mô hình EFDC 34 2.5 Cơ sở nghiên cứu ảnh hưởng chân vịt tàu thủy tới tượng bồi, xói luồng tàu 42 2.5.1 Cơ sở lý thuyết nghiên cứu ảnh hưởng chân vịt tàu thủy .42 2.5.2 Cơ sở tính ảnh hưởng chân vịt tàu thủy mô hình EFDC .47 2.6 Kết luận chương 52 CHƯƠNG NGHIÊN CỨU CHẾ ĐỘ THỦY ĐỘNG LỰC HỌC VÀ BỒI, XÓI KHU VỰC CỬA LẤP .53 3.1 Xây dựng sở liệu mô hình tính toán 53 3.1.1 Xây dựng sở liệu khu vực tính toán 53 3.1.2 Lựa chọn xác định miền tính toán 57 3.1.3 Kết hiệu chỉnh mô hình cho số liệu khảo sát tháng 10 năm 2010 .60 3.1.4 Xây dựng tương quan để tính toán mở rộng cho số liệu năm 2013 .64 3.2 Nghiên cứu, đánh giá chế độ thủy động lực học trình bồi, xói 67 3.2.1 Chế độ thủy động lực học khu vực Cửa Lấp 67 3.2.2 Đánh giá ảnh hưởng sóng đến dòng chảy ven bờ vận chuyển bùn cát khu vực Cửa Lấp 74 3.2.3 Cơ chế vận chuyển bùn cát khu vực Cửa Lấp 80 3.2.4 tàu Đánh giá ảnh hưởng chân vịt tàu tới tượng bồi lấp trở lại luồng 86 3.2.5 Phân tích, đánh giá kết tính toán 89 3.3 Kết luận chương 91 iv CHƯƠNG ỨNG DỤNG KẾT QUẢ NGHIÊN CỨU ĐÁNH GIÁ HIỆU QUẢ KỸ THUẬT GIẢI PHÁP CHỈNH TRỊ CỬA LẤP .93 4.1 Cơ sở để bố trí không gian công trình chỉnh trị Cửa Lấp 93 4.1.1 Mục tiêu chỉnh trị 93 4.1.2 Tham số thiết kế luồng tàu 93 4.1.3 Tham số thiết kế công trình bảo vệ tôn tạo bờ biển .96 4.2 Các phương án bố trí công trình chỉnh trị 97 4.2.1 Nhiệm vụ công trình chỉnh trị khu vực Cửa Lấp 97 4.2.2 Các phương án chỉnh trị đề xuất từ nghiên cứu có 97 4.2.3 Phân tích lựa chọn phương án công trình chỉnh trị khu vực Cửa Lấp 100 4.3 Đánh giá hiệu kỹ thuật công trình chỉnh trị 102 4.3.1 Chế độ thủy động lực diến biến hình thái theo phương án chỉnh trị 102 4.3.2 Đánh giá trình bồi lấp trở lại luồng tàu ảnh hưởng chân vịt tàu thủy lúc di chuyển .104 4.4 Kết luận chương 107 KẾT LUẬN VÀ KIẾN NGHỊ .109 I. Những kết đạt luận án .109 II. Những đóng góp luận án 111 III. Những tồn kiến nghị hướng phát triển .111 DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ 113 TÀI LIỆU THAM KHẢO .114 CÁC PHỤ LỤC .120 v DANH MỤC CÁC HÌNH Hình 2-1 Phương pháp tiếp cận vấn đề nghiên cứu 29 Hình 2-2 Sơ đồ vận chuyển bùn cát dòng chảy 37 Hình 2-3 Bùn cát bị xáo trộn ảnh hưởng chân vịt tàu thủy 43 Hình 2-4 Các yếu tố ảnh hưởng tới vận chuyển bùn cát chân vịt tàu thủy 43 Hình 2-5 Mô tả trình tính toán sa bồi dòng chảy chân vịt tàu gây 44 Hình 2-6 Định nghĩa vùng vùng cho phân bố vận tốc sau tàu .45 Hình 2-7 Chia lưới cho tàu mô hình thủy lực .48 Hình 2-8 Sơ đồ khối tính toán vận chuyển bùn cát tàu mô hình EFDC .49 Hình 2-9 Sơ đồ khối chi tiết mô đun tính toán vận chuyển bùn cát tàu .50 Hình 2-10 Giao diện mô đun tính toán ảnh hưởng chân vịt tàu thủy .51 Hình 2-11 Thông tin quyền .51 Hình 3-1 Tổng hợp địa hình khu vực Cửa Lấp (a) nội suy số liệu địa hình khoảng cách 5x5 m (b) .54 Hình 3-2 Hoa gió trạm Vũng Tàu từ 01/2012 đến 12/2013 54 Hình 3-3 Vị trí trạm đo sóng dòng chảy .55 Hình 3-4 Mực nước thực đo vị trí khu vực tính toán .55 Hình 3-5 Lưu lượng nước thực đo vị trí Cửa Lấp 56 Hình 3-6 Vận tốc thực đo trạm kiểm tra Cửa Lấp 56 Hình 3-7 Vị trí lấy mẫu bùn cát đợt khảo sát năm 2010 (a), kết phân tích thành phần hạt đưa vào mô hình (b) .57 Hình 3-8 Sơ hoạ phạm vi mô 57 Hình 3-9 Kết so sánh mực nước thực đo mực nước tính toán trạm Vũng Tàu .59 Hình 3-10 Lưới mịn lấy điều kiện biên thủy lực cách trích xuất kết từ mô hình lưới thô 59 Hình 3-11 Mực nước thực đo tính toán trạm Cửa Lấp năm 2010 61 Hình 3-12 Vận tốc thực đo tính toán vị trí 0,8 H 61 Hình 3-13 Vận tốc thực đo tính toán vị trí 0,6 H 62 Hình 3-14 Vận tốc thực đo tính toán vị trí 0,2 H 62 Hình 3-15 Kết so sánh hướng sóng vị trí Cửa Lấp 62 Hình 3-16 Kết so sánh độ cao sóng vị trí Cửa Lấp 63 Hình 3-17 Các phân đoạn tính tương quan 65 Hình 3-18 Tương quan mực nước Cửa Lấp trạm triều Vũng Tàu theo đoạn kỳ triều với số liệu đo năm 2013 .65 Hình 3-19 Tương quan lưu lượng mực nước Cửa Lấp với số liệu đo năm 2013 .66 Hình 3-20 Mực nước thực đo tính toán trạm Cửa Lấp năm 2013 66 Hình 3-21 Trường vận tốc dòng chảy: (a) triều lên, (b) triều xuống 67 Hình 3-22 Vị trí khảo sát vận tốc dòng chảy .68 Hình 3-23 Biến thiên độ lớn vận tốc dòng chảy vi trí khác (2013) 68 vi Hình 3-24 Trường dòng chảy dư tổng cộng: a) gió mùa Đông Bắc; b) gió mùa Tây Nam .69 Hình 3-25 Chế độ dòng chảy khu vực Cửa Lấp .71 Hình 3-26 Vị trí mặt cắt khảo sát lưu lượng dòng chảy 72 Hình 3-27 Lưu lượng qua: (a) Mặt cắt 1, (b) Mặt cắt 2, (c) Mặt cắt 72 Hình 3-28 Vị trí hoa tần suất vận tốc dòng chảy 73 Hình 3-29 Trường dòng chảy khu vực Cửa Lấp sóng gió mùa Đông Bắc 75 Hình 3-30 Trường dòng chảy khu vực Cửa Lấp sóng gió mùa Tây Nam 77 Hình 3-31 Điểm khảo sát ứng suất đáy sóng gió hướng Đông Bắc Tây Nam 79 Hình 3-32 Kết so sánh ứng suất đáy sóng gió hướng Đông Bắc Tây Nam 79 Hình 3-33 Phân bố hàm lượng bùn cát kết dính (a), bùn cát không dính (b) .81 Hình 3-34 Ứng suất đáy: (a) dòng chảy, (b) sóng gió mùa Đông Bắc 81 Hình 3-35 Ứng suất đáy: (a) dòng chảy, (b) sóng gió mùa Tây Nam .81 Hình 3-36 Vị trí khảo sát ứng suất đáy 83 Hình 3-37 Kết so sánh ứng suất đáy sóng dòng chảy vị trí khảo sát.83 Hình 3-38 Biến đổi tổng lượng bùn cát lơ lửng điểm khảo sát mô hình toán .84 Hình 3-39 Ảnh chụp vệ tinh khu vực Cửa Lấp tháng 8/2010 tháng 10/2013 84 Hình 3-40 Kết mô diễn biến hình thái đáy khu vực Cửa Lấp năm 2013 .85 Hình 3-41 Vị trí neo đậu tàu thuyền (a) tuyến luồng tàu vào cảng Phước Tỉnh (b) .86 Hình 3-42 Cắt dọc địa hình đáy luồng tàu .86 Hình 3-43 Vận tốc, ứng suất, nồng độ bùn cát tàu: (a) Trường vận tốc; (b) Mặt cắt vận tốc dòng chảy, ứng suất đáy nồng độ bùn cát lơ lửng vị trí 1m từ tim tàu .88 Hình 3-44 Trường ứng suất đáy bùn cát lơ lửng gây chân vịt tàu .88 Hình 3-45 Mô nồng độ bùn cát lơ lửng có tàu tàu chạy qua .89 Hình 4-1 Tuyến lạch sâu vùng Cửa Lấp 94 Hình 4-2 Mặt cắt ngang luồng đào kênh chạy tàu .96 Hình 4-3 Sơ họa bố trí cụm công trình chỉnh trị Cửa Lấp PA1, PA2 .98 Hình 4-4 Biến đổi mức đáy khu vực Cửa Lấp theo phương án, gió Tây Nam 100 Hình 4-5 Biến đổi mức đáy khu vực Cửa Lấp theo phương án, gió Đông Bắc .100 Hình 4-6 Mặt sơ họa bố trí cụm công trình chỉnh trị Cửa Lấp (PA3) 102 Hình 4-7 Chế độ thủy động lực khu vực Cửa Lấp: a) không nạo vét; b) có nạo vét 103 Hình 4-8 Diến biến hình thái khu vực Cửa Lấp: a, b) sau tháng nạo vét không nạo vét; c, d) sau 11 tháng nạo vét không nạo vét 103 Hình 4-9 Diễn biến bồi lấp trở lại luồng tàu khu vực Cửa Lấp .105 vii DANH MỤC CÁC BẢNG Bảng 1.1 Các cách phân loại cửa sông giới 10 Bảng 1.2 Các cách phân loại cửa sông Việt Nam 11 Bảng 1.3 Thống kê số công trình chỉnh trị cửa sông Miền Trung .20 Bảng 3.1 Vị trí trạm hiệu chỉnh .55 Bảng 3.2 Các thông số lưới thô lưới mịn mô hình 58 Bảng 3.3 Chỉ tiêu đánh giá sai số thực đo tính toán trạm Cửa Lấp .61 Bảng 3.4 Các kịch tính toán 75 Bảng 3.5 Ứng suất tới hạn tương ứng với kích cỡ bùn cát (theo Julien - 1998) .82 Bảng 3.6 Thống kê thông số tàu điển hình khu vực Cửa Lấp .87 Bảng 3.7 Ảnh hưởng chân vịt tàu thủy theo độ sâu chạy tàu .89 Bảng 4.1 Kích thước loại tàu điển hình 93 Bảng 4.2 Tọa độ điểm tuyến luồng tàu số 94 Bảng 4.3 Bảng tính độ sâu luồng tàu thiết kế cho loại tàu điển hình 95 Bảng 4.4 Tần suất mực nước trạm Vũng Tàu .95 Bảng 4.5 Mực nước chạy tàu thiết kế cho loại tàu 96 Bảng 4.6 Tọa độ tuyến chỉnh trị bờ biển phía Tây Cửa Lấp .96 Bảng 4.7 Thống kê công trình cụm Cửa Lấp PA1 98 Bảng 4.8 Thống kê công trình cụm Cửa Lấp PA2 99 Bảng 4.9 Thống kê cụm công trình chỉnh trị khu vực Cửa Lấp PA3 101 Bảng 4.10 Mức độ bồi lấp trở lại luồng tàu .106 viii CÁC KÝ HIỆU VIẾT TẮT VÀ GIẢI THÍCH CÁC THUẬT NGỮ CV Mã lực (Cheval Vapeur) DWT Đơn vị đo lực vận tải an toàn tàu tính ĐH Đại học EFDC Mô hình thủy động lực học môi trường chiều EFDC (Environmental Fluid Dynamics Code) GDP Tổng sản phẩm quốc nội (Gross Domestic Products) KH-CN Khoa học công nghệ KTTV Khí tượng thủy văn KTXH Kinh tế xã hội MC Mặt cắt PA Phương án ix If (UStar = 0#) Then GoTo 25 dgr = D50 * ((GRAV * 1.65 / VNU ^ 2) ^ 0.3333) If (dgr = 21 Then RMILE = CSng(OVal(1)) CDUM1 = OVal(2) SPEED = OVal(3) SIZE = OVal(4) CDUM1 = OVal(5) CDUM1 = OVal(6) STAGE = OVal(7) CDUM1 = OVal(8) TRAFFIC = OVal(9) STGSL = OVal(10) CELLID = OVal(11) DEPTH = CDbl(OVal(12)) CDUM1 = (OVal(13)) CDUM1 = (OVal(14)) DTSL = CDbl(OVal(15)) HMAX = CDbl(OVal(16)) CDUM1 = (OVal(17)) TPSHR = CDbl(OVal(18)) CDUM1 = (OVal(19)) CDUM1 = (OVal(20)) TAMB = CDbl(OVal(21)) End If nReads = nReads + If (depthx = 12 Then RTEST = CDbl(OVal(1)) CELLID1 = OVal(2) CDUM1 = OVal(3) CDUM1 = OVal(4) CDUM1 = OVal(5) CDUM1 = OVal(6) CDUM1 = OVal(7) CDUM1 = OVal(8) GS = CDbl(OVal(9)) VS = CDbl(OVal(10)) COHSED = CInt(OVal(11)) SEDTYP = CInt(OVal(12)) End If RMILED = Abs(RMILE - RTEST) If (RMILED < 0.01) Then GoTo 120 Next 120 'Backspace (60) !NMD NEED MORE WORK ON THIS RMILE1 = RTEST IVLEN = For jj = To 100000 If jj > Then If EOF(60) Then Exit For Line Input #60, Lin$ Call Parse_Line(Lin$, OVal(), nTokens) If nTokens >= 12 Then RTEST = CDbl(OVal(1)) CELLIDV(jj) = OVal(2) CDUM1 = OVal(3) CDUM1 = OVal(4) CDUM1 = OVal(5) CDUM1 = OVal(6) CDUM1 = OVal(7) CDUM1 = OVal(8) GSV(jj) = CDbl(OVal(9)) VSV(jj) = CDbl(OVal(10)) COHSEDV(jj) = CInt(OVal(11)) SEDTYPV(jj) = CInt(OVal(12)) End If 456 If (HMAX < 0#) Then If (HMAX -9#) Then HMAX = Abs(HMAX) End If End If DTSL = Abs(DTSL) RMILED = Abs(RMILE - RMILEB) If (RMILED > 0.01 Or TRAFFIC TRAFFICB _ Or STGSL STGSLB) Then Exit For If (DTSL < DTSLMIN) Then DTSLMIN = DTSL TPSHB = TPSHR DASL = DEPTH CELLIDC = CELLID End If Next 110 IHAVE = Close #70 Open NAVFOUT For Input As #70 Line Input #70, Lin$ nReads = nReads - nLines For jj = To nReads Line Input #70, Lin$ 'Backspace (70) !NMD NEED TO WORK ON THIS Next If EOF(70) Then Exit Do Line Input #70, Lin$ Call Parse_Line(Lin$, OVal(), nTokens) If nTokens >= 21 Then RMILE = CSng(OVal(1)) CDUM1 = OVal(2) SPEED = OVal(3) SIZE = OVal(4) CDUM1 = OVal(5) CDUM1 = OVal(6) STAGE = OVal(7) CDUM1 = OVal(8) TRAFFIC = OVal(9) STGSL = OVal(10) CELLID = OVal(11) DEPTH = CDbl(OVal(12)) CDUM1 = CDbl(OVal(13)) CDUM1 = CDbl(OVal(14)) DTSL = CDbl(OVal(15)) HMAX = CDbl(OVal(16)) CDUM1 = (OVal(17)) 144 Else CELLIDV(jj) = CELLID1 GSV(jj) = GS VSV(jj) = VS COHSEDV(jj) = COHSED SEDTYPV(jj) = SEDTYP End If RMILED = Abs(RMILE - RTEST) If (RMILED > 0.01) Then GoTo 130 IVLEN = IVLEN + Next 130 'continue IRELEV = End If For jj = To IVLEN If (CELLIDV(jj) = CELLID) Then GS = GSV(jj) VS = VSV(jj) COHSED = COHSEDV(jj) SEDTYP = SEDTYPV(jj) For i_bw = To i_bwcnt If (Trim(CELLID) = Trim(bw_cellid(i_bw))) Then COHSED = bw_cohsed(i_bw) SEDTYP = bw_sedtyp(i_bw) GoTo 1234 End If Next 1234 End If If (CELLIDV(jj) = CELLIDC) Then GSC = GSV(jj) VSC = VSV(jj) End If Next GS = GS / 1000# VS = VS / 100# GSC = GSC / 1000# VSC = VSC / 100# If (SPEED = "S") Then VBRG = VBRGS ElseIf (SPEED = "M") Then VBRG = VBRGM ElseIf (SPEED = "F") Then VBRG = VBRGF End If If (SIZE = "L") Then BLENGTH = BLENGTHL BWIDTH = BWIDTHL ElseIf (SIZE = "S") Then BLENGTH = BLENGTHS BWIDTH = BWIDTHS ElseIf (SIZE = "M") Then BLENGTH = BLENGTHM BWIDTH = BWIDTHM ElseIf (SIZE = "B") Then BLENGTH = BLENGTHB BWIDTH = BWIDTHB End If If (STAGE = "H") Then IDSF = ElseIf (STAGE = "M") Then IDSF = ElseIf (STAGE = "L") Then IDSF = End If iDS = iDS + IDSF If (iDS = 11) Then DSNSBS = "p26hq" GSM = 0.0005 VSM = 0.06 TSHIFT = 140# NPASS = ElseIf (iDS = 12) Then DSNSBS = "p26mq" GSM = 0.0005 VSM = 0.06 TSHIFT = 140# NPASS = ElseIf (iDS = 13) Then DSNSBS = "p26lq" GSM = 0.0005 VSM = 0.06 TSHIFT = 140# NPASS = ElseIf (iDS = 21) Then DSNSBS = "p8hqbs" GSM = 0.0005 VSM = 0.06 TSHIFT = 167# NPASS = ElseIf (iDS = 22) Then DSNSBS = "2p8mqbs" GSM = 0.0005 VSM = 0.06 TSHIFT = 140# NPASS = ElseIf (iDS = 23) Then DSNSBS = "p8lqbs" DSNSRC = "p8lqrc" GSM = 0.0005 VSM = 0.06 TSHIFT = 103# NPASS = ElseIf (iDS = 31) Then DSNSBS = "lghqsbs" DSNSRC = "lghqsrc" DSNFBS = "lqhqfbs" DSNFRC = "lghqfrc" GSM = 0.0001 VSM = 0.00658 TSHIFT = 213# NPASS = ElseIf (iDS = 32) Then DSNSBS = "lgmqsbs" DSNSRC = "lglqsrc" DSNFBS = "lgmqfbs" DSNFRC = "lgmqfrc" GSM = 0.0001 VSM = 0.00658 TSHIFT = 219# NPASS = ElseIf (iDS = 33) Then DSNSBS = "lglqsbs" DSNSRC = "lglqsrc" DSNFBS = "lglqfbs" DSNFRC = "lglqfrc" GSM = 0.0001 VSM = 0.00658 TSHIFT = 175# NPASS = End If INOREAD = If (iDS = IDSOLD) Then INOREAD = 145 IDSOLD = iDS iDS = iDS - IDSF If (GS > 9#) Then GS = GSM If (GSC > 9#) Then GSC = GSM If (VS > 90#) Then VS = VSM If (VSC > 90#) Then VSC = VSM If (iDS >= 31 And iDS = NSL(L) Then For i = To NSL(L) SL(i, L) = CDbl(OVal(i)) Next End If Line Input #10, Lin$ Call Parse_Line(Lin$, OVal(), nTokens) If nTokens >= NSL(L) Then For i = To NSL(L) TPSHM(i, L) = CDbl(OVal(i)) Next End If Line Input #10, Lin$ Call Parse_Line(Lin$, OVal(), nTokens) If nTokens >= NSL(L) Then For i = To NSL(L) CFAMB(i, L) = CDbl(OVal(i)) Next End If For i = To NST(L) Line Input #10, Lin$ Call Parse_Line(Lin$, OVal(), nTokens) ST(i, L) = CDbl(OVal(1)) If nTokens >= NSL(L) + Then For j = To NSL(L) CONM(i, j, L) = CDbl(OVal(i + 1)) Next End If ST(i, L) = ST(i, L) + TSHIFT If (L = 3) Then For LL = To NSL(L) CONM(i, LL, L) = CONM(i, LL, L) * 0.67 Next End If If (L = 4) Then For LL = To NSL(L) CONM(i, LL, L) = CONM(i, LL, L) * 0.67 Next End If Next If (DFNAME = DSNAME) Then GoTo 10 Next 10 'continue 'REWIND (10) !NMD NEED MORE WORK ON THIS Close #10 Open Tmp$ For Input As #10 Next End If For i = To 1000 COUT(i) = 0# Next For L = To NPASS If (L = Or L = 4) Then For i = To NSL(L) XDAT(i) = SL(i, L) YDAT(i) = TPSHM(i, L) Next Call Value(NSL(L), DTSL, XDAT, YDAT, 0#, 1#, 1#, TPSHM1) Else TPSHM1 = TPSHM(1, L) End If For i = To NST(L) For j = To NSL(L) XDAT(j) = SL(j, L) YDAT(j) = CONM(i, j, L) Next Call Value(NSL(L), DTSL, XDAT, YDAT, 0#, 1#, 1#, CASD(i)) Next If (L = 1) Then ratio = CSAMB(TPSHB, GSC, VSC) / CSAMB(TPSHM1, GSM, VSM) RPD = 1# ElseIf (L = 2) Then ratio = CSAMB(TPSHR, GSC, VSC) / CSAMB(TPSHM1, GSM, VSM) DTSLC = (BWIDTH + 20#) / 2# If (DTSL < DTSLC) Then ratio = 0# RPD = 1# ElseIf (L = 3) Then TC = 0.3 If (IDSF = 1) Then TC = 1.2 ratio = (TPSHB - TC) / (TPSHM1 - TC) If (ratio < 0#) Then ratio = 0# RPD = dmax((1# - TAMB / TC), 0#) Else TC = 0.3 If (IDSF = 1) Then TC = 1.2 ratio = (TPSHR - TC) / (TPSHM1 - TC) If (ratio < 0#) Then ratio = 0# If (DTSL < 25#) Then ratio = 0# RPD = aMax((1# - TAMB / TC), 0#) End If If (L = 1) Then Camb = CSAMB(TAMB, GSC, VSC) ElseIf (L = 3) Then For j = To NSL(L) XDAT(j) = SL(j, L) YDAT(j) = CFAMB(j, L) Next Call Value(NSL(L), DTSL, XDAT, YDAT, 0#, 1#, 1#, CFINES) Camb = CFINES Else Camb = 0# End If For i = To NST(L) CASD(i) = CASD(i) * ratio Next For i = To NST(L) XDAT(i) = ST(i, L) YDAT(i) = CASD(i) 146 Next For i = To nPts VS1 = VSC If (L > 2) Then VS1 = 0.0025 Call Value(NST(L), CSng(tOut(i - 1)), XDAT, YDAT, VS1, DEPTH, RPD, COUT1) COUT(i) = COUT(i) + COUT1 + Camb Next Next If (depthx > 0#) Then depth_orig = DEPTH For i = To i_bwcnt If (Trim(CELLID) = Trim(bw_cellid(i))) Then i_writeout = -1 If (DEPTH > depthx) Then DEPTH = depthx GoTo 444 End If End If Next Else i_writeout = -1 End If 444 tvrmax = 1# If (DEPTH 3000#) Then YDAT(i) = 3000# ElseIf (SEDTYP = 2) Then If (YDAT(i) > 1500#) Then YDAT(i) = 1500# End If Next COUT1 = 0# For i = To nPts Call Value(numpnts, CSng(tOut(i - 1)), XDAT, YDAT, VS, DEPTH, 1#, COUT1) If (depthx > 0#) Then If (depth_orig > depthx) Then COUT1 = COUT1 * (depthx ^ / depth_orig ^ 2) End If End If COUT(i) = COUT(i) + COUT1 Next End If NLNS = nPts / 10 nRM = nPts Mod 10 If (IONCE = 0) Then IONCE = Print #5, "START TIME END TIME NO. TIME STEPS" Print #5, PushRight(Format(TIMES, "0.00"), 14) + PushRight(Format(TIMEE, "0.00"), 14) + PushRight(nPts, 10) Print #5, "TIME VALUES (SECONDS)" 'WRITE(5,13) TIMES,TIMEE,NPTS k=1 For i = To NLNS Lin$ = "" For j = k To k + Lin$ = Lin$ + PushRight(Format(tOut(j - 1), "0.00"), 14) Next Print #5, Lin$ 'WRITE(5,15) (TOUT(J),J=K,K+9) k = k + 10 Next If (nRM > 0) Then NRMM1 = nRM - Lin$ = "" For j = k To k + NRMM1 Lin$ = Lin$ + PushRight(Format(tOut(j - 1), "0.00"), 14) Next Print #5, Lin$ 'WRITE(5,15) (TOUT(J),J=K,K+NRMM1) End If End If If (i_writeout 0) Then i_writeout = Print #5, CELLID + " " + TRAFFIC + " " + STGSL 'WRITE(5,14) CELLID,TRAFFIC,STGSL iwrite = For i = To nPts If ((COUT(i) - COUT(1)) >= 5#) Then iwrite = Next Print #5, CStr(iwrite) + PushRight(Format(COUT(1), "0.00"), 14) _ + PushRight(Format(DEPTH, "0.00"), 14) 'WRITE(5,16) IWRITE, COUT(1),DEPTH If (iwrite = 1) Then k=1 For i = To NLNS Lin$ = "" For j = k To k + Lin$ = Lin$ + PushRight(Format(COUT(j), "0.00"), 14) Next Print #5, Lin$ ' WRITE(5,15) (COUT(J),J=K,K+9) k = k + 10 Next If (nRM > 0) Then NRMM1 = nRM - Lin$ = "" For j = k To k + NRMM1 Lin$ = Lin$ + PushRight(Format(COUT(j), "0.00"), 14) Next Print #5, Lin$ 'WRITE(5,15) (COUT(J),J=K,K+NRMM1) End If End If End If End If If (INODAT = 0) Then IRECORD = IRECORD + Loop 'MAIN LOOP END Close #5 Close #10 Close #50 Close #60 Close #70 frmPropwash.MousePointer = vbDefault frmPropwash.Refresh ExitSub: Exit Sub NavsedError: MsgBox Error(Err) GoTo 5000 End Sub 147 Private Sub Value(NWP As Integer, DLOC As Double, XDAT() As Double, YDAT() As Double, _ VS As Double, DEPTH As Double, RPD As Double, VOUT As Double) Dim TC As Double Dim VALUET As Double Dim i As Integer TC = -VS * RPD / DEPTH VALUET = YDAT(1) If(XDAT(1)>DLOC) Then VALUET=YDAT(1)*Exp(TC*(XDAT(1)-DLOC)) For i = To NWP - If (XDAT(i) DLOC) Then VALUET = (DLOC - XDAT(i)) * (YDAT(i + 1) - YDAT(i)) / (XDAT(i + 1) - XDAT(i)) + YDAT(i) End If Next If (XDAT(NWP) 1.5) Then h0 = 1.5 End If tbp = tl / tvelg For i = To RATOT '!total number of ratios defined in TRATIO() and VRATIO() Call calratio(TRATIO(i), tbp, tt) '!tt is the value "returned" Call calratio(VRATIO(i), tvrmax, vr) 'vr is the value "returned" Next Close (1) If (SEDTYP = 1) Then RHOM = RHOSFT tcp1 = TCP1SF tcp2 = TCP2SF ERP1 = ERP1SFT ERP2 = ERP2SFT Call wavavg(HMAX, dtrun, waveperiod, r4slopeF, filewave, numpnts) Call vestuns(h0, filewave, vestout, time, conc) ElseIf (SEDTYP = 2) Then RHOM = RHOMED tcp1 = TCP1MD tcp2 = TCP2MD ERP1 = ERP1MED ERP2 = ERP2MED Call wavavg(HMAX, dtrun, waveperiod, r4slopeF, filewave, numpnts) Call vestuns(h0, filewave, vestout, time, conc) ElseIf (SEDTYP = 3) Then RHOM = RHOHRD tcp1 = TCP1HD tcp2 = TCP2HD ERP1 = ERP1HRD ERP2 = ERP2HRD '* --- vestuns+acker Call wavavg(HMAX, dtrun, waveperiod, r4slopeF, filewave, numpnts) Call vestuns(h0, filewave, vestout, time, conc) End If '* stop End Sub Private Sub calratio(ratio As Double, tbp As Double, ret As Double) ret = tbp * ratio Private Sub wavavg(HMAX As Double, dtrun As Double, waveperiod As Integer, r4slopeF As Integer, filewave As String, noi As Integer) Dim i, j, k, nPts, int0, xint, iX, iw1, iw2, iw3, iw4 As Integer Dim SUM, irang1, irang2, irang3, irang4, iPos As Integer Dim h(3000), H1, h3, h4, h5, xm1, xm2, xm3, xm4, Avg As Double 'open(unit=6,file=FILEWAVE,status='unknown') Open filewave For Output As #6 int0 = iX = '* npts = 1915 nPts = dtrun * 30 ' start of specified values iw1 = 25 iw2 = 75 iw3 = 200 iw4 = nPts '* units: cms H1 = 0# h3 = HMAX * 0.2 '* if (h3 = tauero) Then fs = ERSMAX * Exp(-ERP1 * tauero ^ ERP2) * (taubed - tauero) Else fs = 0# End If dzbt = dzb + dt * fs / rhob If (dzbt 0.000001) Then qkk = 0.5 * (qk + qkk) GoTo End If wnu = 0.5 * (qk + qkk) wlength = 2# * 3.14159265 / wnu wfr = 2# * PI / wperiod For i = To NGRIDS + Z = (Float(i - 1)) * dz dvel(i) = wheight / 2# * wfr * wnu * SinH(wnu * Z) / SinH(wnu * h0) kn(i) = WDIFFK * wfr * wheight ^ 2# * SinH(wnu * h0 / (h0 + RD) * (Z + RD)) * cosh(wnu * Z) / (2# * (SinH(wnu * h0)) ^ 2#) Next rig = g * (Visc * wperiod) ^ 0.5 * (RHOM - RHOW) / RHOW / (dudb * dudb) ub = sig / 2# * wheight / (SinH(wnu * h0)) End Sub Private Sub hydrowave2(h0, dvel, kn, taubed) Dim dvel(1), kn(1) As Double sig = 2# * 3.14159265 / wperiod qkk = (sig * sig) / 9.8 qk = sig ^ 2# / (9.8 * TanH(qkk * h0)) If (Abs(qk - qkk) > 0.000001) Then qkk = 0.5 * (qk + qkk) GoTo End If wnu = 0.5 * (qk + qkk) wlength = 2# * 3.14159265 / wnu wfr = 2# * PI / wperiod For i = To NGRIDS + Z = (Float(i - 1)) * dz dvel(i) = wheight / 2# * wfr * wnu * SinH(wnu * Z) / SinH(wnu * h0) kn(i) = WDIFFK * wfr * wheight ^ 2# * (SinH(wnu * Z)) * cosh(wnu * Z) / (2# * (SinH(wnu * h0)) ^ 2#) Next taubed= RHOW * WFRIC / 2# *(wheight / 2# * wfr / SinH(wnu* h0)) ^ 2# End Sub Private Sub wavecurrent1(h0, dvel, kn, RD) 152 Dim dvel(1), kn(1) As Double Dim knc(500), knw(500), RD As Double Dim dvelc(500), dvelw(500) As Double kkss = 3.15 * 100000# * 9.8 ^ 3# * QN ^ 6# Y0 = kkss / 30.1 For i = To NGRIDS + Z = (Float(i - 1)) * dz aatt = alog(h0 / Y0) dvelc(i) = uu / Z / (aatt - 1#) knc(i) = (KK * KK * uu) / (aatt - 1#) * Z * (h0 - Z) / h0 Next sig = 2# * 3.14159265 / wperiod qkk = (sig * sig) / 9.8 qk = sig ^ 2# / (9.8 * TanH(qkk * h0)) If (Abs(qk - qkk) > 0.000001) Then qkk = 0.5 * (qk + qkk) GoTo End If wnu = 0.5 * (qk + qkk) wlength = 2# * 3.14159265 / wnu wfr = 2# * PI / wperiod For i = To NGRIDS + Z = (Float(i - 1)) * dz dvelw(i) = wheight / 2# * wfr * wnu * SinH(wnu * h0 / (h0 + RD) * (Z + RD)) / SinH(wnu * h0) knw(i) = WDIFFK * wfr * wheight ^ 2# * SinH(wnu * h0 / (h0 + RD) * (Z + RD)) * cosh(wnu * Z) / (2# * (SinH(wnu * h0)) ^ 2#) dvel(i) = dvelw(i) + dvelc(i) kn(i) = knw(i) + knc(i) Next rig =g *(Visc * wperiod) ^ 0.5 *(RHOM - RHOW) / RHOW / (dudb * dudb) ub = sig / 2# * wheight / (SinH(wnu * h0)) End Sub Private Sub wavecurrent2(h0, dvel() As Double, kn() As Double, taubed) Dim knc(500), knw(500), RD As Double Dim dvelc(500), dvelw(500) As Double ReDim Preserve dvel(1), kn(1) kkss = 3.15 * 100000# * 9.8 ^ 3# * QN ^ 6# Y0 = kkss / 30.1 For i = To NGRIDS + Z = (Float(i - 1)) * dz aatt = alog(h0 / Y0) dvelc(i) = uu / Z / (aatt - 1#) knc(i) = (KK * KK * uu) / (aatt - 1#) * Z * (h0 - Z) / h0 Next ff = 8# * 9.8 * QN ^ 2# / h0 ^ (1# / 3#) taubedc = RHOW * ff * uu ^ 2# / 8# sig = 2# * 3.14159265 / wperiod qkk = (sig * sig) / 9.8 qk = sig ^ 2# / (9.8 * TanH(qkk * h0)) If (Abs(qk - qkk) > 0.000001) Then qkk = 0.5 * (qk + qkk) GoTo End If wnu = 0.5 * (qk + qkk) wlength = 2# * 3.14159265 / wnu wfr = 2# * PI / wperiod For i = To NGRIDS + Z = (Float(i - 1)) * dz dvelw(i) = wheight / 2# * wfr * wnu * SinH(wnu * Z) / SinH(wnu * h0) knw(i) = WDIFFK * wfr * wheight ^ 2# * SinH(wnu * Z) * cosh(wnu * Z) / (2# * (SinH(wnu * h0)) ^ 2#) dvel(i) = dvelw(i) + dvelc(i) kn(i) = knw(i) + knc(i) Next taubedw=RHOW*WFRIC/2# *(wheight/2#*wfr/SinH(wnu*h0)) ^ 2# taubed = taubedw + taubedc End Sub Private Sub result(h0, tt, c, WS, cbar) '* cbar is returned Dim tt As Double Dim c(1), WS(1) As Double suspmas=0# wssum = 0# For i = NGRIDS To Step -1 zz = (Float(i) - 0.5) * dz suspmas = suspmas + c(i) * dz wssum = wssum + WS(i) Next cbar = suspmas / h0 * 1000 wbar = wssum / NGRIDS Print #11, PushRight(Format(tt, "0.0000"), 16) + PushRight(Format(wheight, "0.0000"), 12) _ + PushRight(Format(wperiod, "0.0000"), 12) + PushRight(Format(cbar, "0.000000"), 14) Private Sub prowmita(PWD1, PWD2, PDEN1, PDEN2, PWP, PAMP1, PMU1, Puu, PGG1, PGG2, iP, RRR, rr2) Dim iP As Integer Dim PWD1, PWD2, PDEN1, PDEN2, PWP, PAMP1, PMU1, Puu, PGG1, PGG2, RRR As Double Dim DEN1, DEN2, WP, AMP1, SIGN, SIGN2, PI As Double Dim WD1, WD2, CRIT, WDY1, WDY2 As Double Dim MU1, NU1 As Double Dim AQ1, BQ0, BQ1 As Double Dim LSINH, LCOSH As Complex Dim A2, B2, C2, D2, E2, F2, G2, H2, AMP2, BAMP2 As Complex Dim MU2, NU2 As Complex Dim A1, B1, C1, D1, E1, F1, G1, H1, BAMP1, WN As Complex Dim UY1, UY2 As Complex Dim ii, LAM1, LAM2, BET1, BET2 As Complex PI = DACOS(-1#) ii = Complex(0#, 1#) SIGN = 1# SIGN2 = 1# CRIT = 0.000000000000001 WD1 = dble(PWD1) WD2 = dble(PWD2) DEN1 = dble(PDEN1) DEN2 = dble(PDEN2) WP = dble(PWP) AMP1 = dble(PAMP1) MU1 = dble(PMU1) gg1 = dble(PGG1) gg2 = dble(PGG2) uu2 = dble(Puu) '***** MUD RHEOLOGICAL MODEL ********* If (PGG1 >= 0#) Then AQ1 = uu2 / (gg1 + gg2) BQ0 = (gg1 * gg2) / (gg1 + gg2) BQ1 = (uu2 * gg1) / (gg1 + gg2) MU2 = ((6.2832 / WP * (BQ1 - BQ0 * AQ1)) _ + ii * (BQ0 + AQ1 * BQ1 * (6.2832 / WP) ^ 2#)) _ / (6.2832 / WP * (1# + (AQ1 * 6.2832 / WP) ^ 2#)) Else MU2 = uu2 + ii * (gg2 * dble(WP / (2# * 3.14159265))) End If NU1 = MU1 / DEN1 NU2 = MU2 / DEN2 '***** CALL JIANG'S (1993) SUBROUTINE ****************** Call CWAVE2ND(WD1, WD2, DEN1, DEN2, MU1, MU2, WP, AMP1, SIGN, _ SIGN2, WN, A1, B1, C1, D1, E1, F1, G1, H1, BAMP1, _ A2, B2, C2, D2, E2, F2, G2, H2, BAMP2, AMP2, CRIT) 153 LAM1 = SIGN * (WN * WN - ii * (6.2832 / WP) / NU1) ^ 0.5 LAM2 = SIGN * (WN * WN - ii * (6.2832 / WP) / NU2) ^ 0.5 BET1 = SIGN2 * (4# * WN * WN - 2# * ii * (6.2832 / WP) / NU1) ^ 0.5 BET2 = SIGN2 * (4# * WN * WN - 2# * ii * (6.2832 / WP) / NU2) ^ 0.5 iP1 = iP WDY1 = DSQRT(NU1 * WP / 2# / PI / 2#) WDY2 = WD2 UY1 = ii * (A1 * LCOSH(WN * WDY1) + B1 * LSINH(WN * WDY1) _ + C1 * LAM1 / WN * CDEXP(LAM1 * (WDY1 - WD1)) _ - D1 * LAM1 / WN * CDEXP(-LAM1 * WDY1)) UY2 = ii * (E1 * LCOSH(WN * WDY2) + F1 * LSINH(WN * WDY2) _ + G1 * LAM2 / WN * CDEXP(LAM2 * (WDY2 - WD2)) _ - H1 * LAM2 / WN * CDEXP(-LAM2 * WDY2)) rr2 = SNGL(CDABS(BAMP1)) WDY1 = DSQRT(NU1 * 2# * WP / 2# / PI) UY1 = ii * (A1 * LCOSH(WN * WDY1) + B1 * LSINH(WN * WDY1) _ + C1 * LAM1 / WN * CDEXP(LAM1 * (WDY1 - WD1)) _ - D1 * LAM1 / WN * CDEXP(-LAM1 * WDY1)) RRR = SNGL(CDABS(UY1 - UY2)) End Sub Private Function LSINH(X) As Complex Dim LSINH As Complex Dim X As Complex LSINH = (CDEXP(X) - CDEXP(-X)) / 2# End Function Private Function LCOSH(X) As Complex Dim LCOSH As Complex Dim X As Complex LCOSH = (CDEXP(X) + CDEXP(-X)) / 2# End Function Private Sub CWAVE2ND(WD1, WD2, DEN1, DEN2, MU1, MU2, WP, AMP1, SIGN, _ SIGN2, WN, A1, B1, C1, D1, E1, F1, G1, H1, BAMP1, _ A2, B2, C2, D2, E2, F2, G2, H2, BAMP2, AMP2, CRIT) Dim Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8 As Complex Dim W1, W2, W3, W4, W5, W6, W7, W8, W9, W10 As Complex Dim W11, W12, W13, W14, W15, W16, W17, W18, W19, W20 As Complex Dim W21, W22, W23, W24, W25, W26, W27 As Complex Dim A2, B2, C2, D2, E2, F2, G2, H2, AMP2, BAMP2 As Complex Dim BETA1, BETA2, SH1, SH2, CH1, CH2, SINH1, SINH2, COSH1, COSH2 As Complex Dim AN1, AN2, DA1, DA2 As Complex Dim A1, B1, C1, D1, E1, F1, G1, H1, BAMP1, WN, DWN As Complex Dim DEN1, DEN2, WP, SMA, AMP1, SIGN, SIGN2 As Double Dim WD1, WD2, WD, DWD1, DWD2, CRIT As Double Dim MU1, NU1, RBAMP1, RBAMP2, RAMP2 As Double Dim MU2, NU2 As Complex WD = WD1 + WD2 NU1 = MU1 / DEN1 NU2 = MU2 / DEN2 SMA = 6.28 / WP Call WN1ST1(WD1, WD2, DEN1, DEN2, MU1, MU2, WP, AMP1, SIGN, _WN, A1, B1, C1, D1, E1, F1, G1, H1, BAMP1, CRIT) DWN = Complex(2#, 0#) * WN BETA1 = SIGN2 * CDSQRT(DWN ^ - Complex(0#, 2#) * SMA / NU1) BETA2 = SIGN2 * CDSQRT(DWN ^ - Complex(0#, 2#) * SMA / NU2) AN1=Complex(0#,1#)*DEN1*SMA/WN-Complex(4#,0#)*WN*MU1 AN2=Complex(0#, #)*DEN2*SMA/WN-Complex(4#, 0#)*WN*MU2 DWD1 = 2# * WD1 DWD2 = 2# * WD2 Call SINHCOSH(WN, DWD1, SH1, CH1) Call SINHCOSH(WN, DWD2, SH2, CH2) Call SINHCOSH(WN, WD1, SINH1, COSH1) Call SINHCOSH(WN, WD2, SINH2, COSH2) DA1 = SIGN * CDSQRT(WN ^ - Complex(0#, 1#) * SMA / NU1) DA2 = SIGN * CDSQRT(WN ^ - Complex(0#, 1#) * SMA / NU2) Y1 = AMP1 * (WN * A1 * COSH1 + WN * B1 * SINH1 + DA1 * C1) Y2 = AMP1 * MU1 * ((Complex(2#, 0#) + DA1 ^ / WN ^ 2) * (WN ^ * ( _ A1 * SINH1 + B1 * COSH1) + DA1 ^ * C1) - (WN ^ * ( _ A1 * SINH1 + B1 * COSH1) + DA1 ^ * C1 / WN ^ 2)) Y3 = -Complex(2#, 0#) * AMP1 * (Complex(2#, 0#) * WN ^ * (A1 * COSH1 + _ B1 * SINH1) + DA1 * (WN ^ + DA1 ^ 2) * C1) Y4 = BAMP1 * (WN * A1 - DA1 * D1) Y5 = Complex(2#, 0#) * BAMP1 * (WN ^ * (E1 * SINH2 + F1 * COSH2) + DA2 ^ * G1 - WN ^ * B1 - DA2 ^ * D1) Y6=AMP1*(WN*(E1*COSH2 -F1*SINH2)+DA2*G1-WN*A1+DA1*D1) Y7 = AMP1 * (-MU2 * (WN ^ + DA2 ^ 2) * (E1 * SINH2 + F1 * COSH2) MU2 * Complex(2#, 0#) * DA2 ^ * G1 + MU1 * (WN ^ + DA1 ^ 2) * B1 + Complex(2#, 0#) * MU1 * DA1 ^ * D1) Y8 = AMP1 * (MU2 * (Complex(4#, 0#) * WN ^ * (E1 * COSH2 + F1 * SINH2)+_Complex(2#,0#)*DA2*(DA2^2+WN^2)*G1)-MU1*(Complex( 4#, 0#)*_WN^3*A1 –Complex (2#,0#)*DA1*(DA1^2+WN^2)*D1)) W1 = BETA2 * SH1 / DWN - CH1 W2 = BETA2 * CH2 - DWN * SH2 W3=AN1*CH1-Complex(0#,1#)*DEN1*9.8*SH1/SMA/Complex(2#,0#) W4=AN1*SH1-Complex(0#,1#)*DEN1*9.8*CH1/SMA/Complex(2#,0#) W5 = DWN * (Complex(2#, 0#) * MU1 * BETA1 + Complex(0#, 1#) * _ DEN1 * 9.8 / SMA / Complex(2#, 0#)) W6 = (W2 * W3 + DWN * W1 * W4) / W5 W7 = BETA2 * W3 / W5 W8 = (BETA1 * W3 - DWN * W4) / W5 W9 = (W3 * Y5 + DWN * W4 * Y6 - DWN * (Y2 + _ Complex(0#, 1#) * DEN1 * 9.8 * Y1 / SMA / Complex(2#, 0#))) / W5 W10 = Complex(8#, 0#) * WN ^ * (W2 * SH1 / DWN + W1 * CH1) + _ W6 * (BETA1 ^ + DWN ^ 2) W11 = Complex(8#, 0#) * WN ^ * BETA2 * SH1 / DWN + _ W7 * (BETA1 ^ + DWN ^ 2) W12 = Complex(8#, 0#) * WN ^ * (BETA1 * SH1 / DWN - CH1) + _ W8 * (BETA1 ^ + DWN ^ 2) W13 = Y3 - Complex(8#, 0#) * WN ^ * (SH1 * Y5 / DWN + CH1 * Y6) - W9 * (BETA1 ^ + DWN ^ 2) W14 = W2 * AN1 / DWN - BETA2 * AN2 * CH2 / DWN + AN2 * SH2 + Complex(0#, 1#) * (DEN2 - DEN1) * 9.8 * W1 / SMA / Complex(2#, 0#) W15 = BETA2 * AN1 / DWN + Complex(2#, 0#) * MU2 * BETA2 * AN2 W16 = BETA1 * AN1 / DWN + Complex(2#, 0#) * MU1 * BETA1 W17 = Y7 - Complex(0#, 1#) * (DEN2 - DEN1) * 9.8 * (Y6 + Y4) / SMA / Complex(2#, 0#) - AN1 * Y5 / DWN W18 = Complex(8#, 0#) * MU1 * WN ^ * W1 - Complex(8#, 0#) * MU2 * WN ^ * (BETA2 * SH2 / DWN - CH2) W19 = -Complex(8#, 0#) * MU2 * WN ^ * (BETA2 ^ + DWN ^ 2) W20 = MU1 * (BETA1 ^ - DWN ^ 2) W21 = Y8 - Complex(8#, 0#) * MU1 * Y6 * WN ^ W22 = W11 * W14 - W15 * W10 W23 = W12 * W14 - W16 * W10 W24 = W13 * W14 - W17 * W10 W25 = W11 * W18 - W19 * W10 W26 = W12 * W18 - W20 * W10 W27 = W13 * W18 - W21 * W10 D2 = (W24 * W25 - W27 * W22) / (W23 * W25 - W26 * W22) G2 = (W24 - W23 * D2) / W22 H2 = (W13 - W11 * G2 - W12 * D2) / W10 F2 = -H2 E2 = BETA2 * H2 / DWN BAMP2=Complex(0#, 1#)*(W1*H2+Y6+Y4)/SMA/Complex(2#, 0#) B2 = W1 * H2 - D2 + Y6 A2 = (H2 * W2 + BETA2 * G2 + BETA1 * D2 + Y5) / DWN C2 = W6 * H2 + W7 * G2 + W8 * D2 + W9 AMP2 = Complex(0#, 1#) * (A2 * SH1 + B2 * CH1 + C2 + Y1) / SMA / Complex(2#, 0#) RBAMP1 = dReal(BAMP1) RBAMP2 = dReal(BAMP2) 154 RAMP2 = dReal(AMP2) End Sub Private Sub WN1ST1(WD1, WD2, DEN1, DEN2, MU1, MU2, WP, AMP1, SIGN, _ WN, A1, B1, C1, D1, E1, F1, G1, H1, SAMP1, CRIT) Dim WFUN1, WFUNX, WN1, WN2, WN, dX As Complex Dim X, FOLD, FNEW, TED, TEG, TEH As Complex Dim A1, B1, C1, D1, E1, F1, G1, H1, SAMP1 As Complex Dim DA1, DA2, SINH1, COSH1, SINH2, COSH2 As Complex Dim WP, SMA, WD1, WD2, DEN1, DEN2, WDTEP As Double Dim WD, MU1, NU1, AMP1, SIGN, CRIT As Double Dim MU2, NU2 As Complex WD = WD1 + WD2 '! TOTAL DEPTH NU1 = MU1 / DEN1 NU2 = MU2 / DEN2 SMA = 6.28 / WP WDTEP = WD1 + WD1 '! FOR ITERATION PURPOSE Call AIRY(WP, WD1, WN1) '!WN1>WN2 Call AIRY(WP, WDTEP, WN2) dX = WN2 - WN1 X = WN2 Call FUNC1ST1(WD1, WD2, DEN1, DEN2, MU1, MU2, WP, SIGN, _ AMP1, WN1, WFUN1, TED, TEG, TEH, H1, _ DA1, DA2, SINH1, COSH1, SINH2, COSH2) FOLD = WFUN1 n=0 100 Call FUNC1ST1(WD1, WD2, DEN1, DEN2, MU1, MU2, WP, SIGN, _ AMP1, X, WFUNX, TED, TEG, TEH, H1, _ DA1, DA2, SINH1, COSH1, SINH2, COSH2) n=n+1 If (n > 100) Then GoTo 200 FNEW = WFUNX dX = -FNEW * dX / (FNEW - FOLD) X = X + dX If (CDABS(dX) > CRIT) Then FOLD = FNEW GoTo 100 End If 200 WN = X D1 = TED * H1 G1 = -TEG * H1 F1 = -H1 E1 = DA2 * H1 / WN SAMP1 = Complex(0#, 1#) * (E1 * SINH2 + F1 * COSH2 + G1) / SMA B1 = E1 * SINH2 + F1 * COSH2 + G1 - D1 A1=(E1*WN*COSH2+F1*WN*SINH2+G1*DA2+D1*DA1) / WN C1 = -Complex(2#, 0#) * NU1 * AMP1 * WN ^ End Sub Private Sub AIRY(PERIOD, DEPTH, WN) Dim WN As Complex Dim GRAV,SUM,PERIOD,DEPTH,Y,CGM,DUM,D(10), KH As Double GRAV = 9.8 CGM = 6.2832 / PERIOD Y = CGM ^ * DEPTH / GRAV D(1) = 0.66666666 D(2) = 0.35555555 D(3) = 0.1608465608 D(4) = 0.0632098765 D(5) = 0.0217540484 D(6) = 0.0065407983 SUM = 1# For j = To DUM = Y ^ j SUM = SUM + D(j) * DUM Next KH = Sqr(Y ^ + Y / SUM) WN = KH / DEPTH End Sub Private Sub FUNC1ST1(WD1, WD2, DEN1, DEN2, MU1, MU2, WP, SIGN, AMP1, WN, WFUN, TED, TEG, TEH, H1, _ DA1, DA2, SINH1, COSH1, SINH2, COSH2) Dim DA1, DA2, s2, Q2, TED, TEG, TEH, WN As Complex Dim AM1, AM2, SINH1, SINH2, COSH1, COSH2 As Complex Dim WFUN, H1 As Complex Dim WP, SMA, WD1, WD2, DEN1, DEN2, NU1 As Double Dim MU1, AMP1, SIGN As Double Dim MU2, NU2 As Complex NU1 = MU1 / DEN1 NU2 = MU2 / DEN2 SMA = 6.28 / WP Call SINHCOSH(WN, WD1, SINH1, COSH1) Call SINHCOSH(WN, WD2, SINH2, COSH2) DA1=SIGN*CDSQRT(WN^2-Complex(0#,1#)*SMA / NU1) '! `+' OR `-' DA2 = SIGN * CDSQRT(WN ^ - Complex(0#, 1#) * SMA / NU2) s2 = SINH2 - COSH2 * DA2 / WN Q2 = COSH2 - SINH2 * DA2 / WN AM1 = Complex(0#, 1#) * SMA * DEN1 / WN - 2# * WN * MU1 AM2 = Complex(0#, 1#) * SMA * DEN2 / WN - 2# * WN * MU2 TEG = (2# * DA1 * (MU2 * Q2 - MU1 * Q2) + _ (AM2 - AM1) * s2 - Complex(0#, 1#) * (DEN2 - DEN1) * 9.8 * Q2 _ / SMA) / (2# * MU1 * DA1 + DA2 * AM1 / WN _ + 2# * MU2 * DA2 + Complex(0#, 1#) * (DEN2 - DEN1) * 9.8 / SMA _ - DA1 * MU2 * (DA2 ^ + WN ^ 2) / WN / WN) TED = (2# * WN ^ * (MU2 * Q2 - MU1 * Q2) - TEG * _ (2# * MU1 * WN ^ - MU2 * (DA2 ^ + _ WN ^ 2))) / SMA / DEN1 / Complex(0#, 1#) TEH = (DA1 * COSH1 / WN - SINH1) * TED - _ TEG * (DA2 * COSH1 / WN + SINH1) - _ s2 * COSH1 - Q2 * SINH1 H1 = (DEN1 * 9.8 * AMP1 - 4# * DEN1 * (NU1 ^ 2) * AMP1 * _ (WN ^ 2) * DA1) / AM1 / TEH WFUN = (DA1 * SINH1 / WN - COSH1) * TED * H1 - _ (DA2 * SINH1 / WN + COSH1) * TEG * H1 - _ (s2 * SINH1 + Q2 * COSH1) * H1 + _ AM1 * WN * AMP1 / DEN1 End Sub '************************************' Private Sub SINHCOSH(WN, WD, SINHX, COSHX) Dim XKH, WN, TEMP1, TEMP2, SINHX, COSHX As Complex Dim WD As Double XKH = WN * WD TEMP1 = CDEXP(XKH) TEMP2 = CDEXP(-XKH) SINHX = (TEMP1 - TEMP2) / 2# COSHX = (TEMP1 + TEMP2) / 2# End Sub 155 Phụ lục 04: Bảng tần suất vận tốc dòng chảy ứng với vị trí khảo sát V(m/s) - .2 .2 - .4 .4 - .6 .6 - .8 .8 - >= Tổng N 16.88 16.36 9.09 3.73 1.8 1.26 49.12 NE 1.91 0 0 1.91 E 1.45 0 0 1.45 V(m/s) - .2 .2 - .4 .4 - .6 .6 - .8 .8 - >= Tổng N 0.06 0 0 0.06 NE 29.78 16.8 6.68 3.57 2.2 0.63 59.65 E 0.53 0 0 0.53 V(m/s) - .1 .1 - .2 .2 - .3 .3 - .4 .4 - .5 >= .5 Tổng N 2.07 0 0 2.07 NE 26.47 33.25 0.37 0 60.09 E 8.42 0.33 0 0 8.75 V(m/s) - .1 .1 - .2 .2 - .3 .3 - .4 .4 - .5 >= .5 Tổng N 10.03 2.48 0.01 0 12.53 NE 10.18 0.2 0 0 10.38 E 19.13 9.45 0.01 0 28.59 V(m/s) - .1 .1 - .2 .2 - .3 .3 - .4 .4 - .5 >= .5 Tổng N 0.61 0.01 0 0 0.62 NE 0.87 1.8 0.97 0.33 0.1 4.07 E 1.09 5.42 15.87 23.59 17.97 9.46 73.41 V(m/s) - .2 .2 - .4 .4 - .6 .6 - .8 .8 - >= Tổng N 0.06 0 0 0.06 NE 2.15 8.14 23.53 21.83 9.58 1.58 66.81 E 3.08 5.92 2.51 0.19 0 11.7 V(m/s) - .2 .2 - .4 .4 - .6 .6 - .8 .8 - >= Tổng N 0.68 0 0 0.68 NE 3.26 1.63 0.01 0 4.9 E 8.5 35.77 26.13 0.88 0 71.28 VỊ TRÍ SE S 5.99 15.15 18.6 7.12 0.25 0 0 5.99 41.12 VỊ TRÍ SE S 0.26 0.43 0 0 0 0 0 0.26 0.43 VỊ TRÍ SE S 1.31 2.2 0 0 0 0 0 1.31 2.2 VỊ TRÍ SE S 18.16 10.46 7.31 1.22 0 0 0 0 25.48 11.69 VỊ TRÍ SE S 0.76 0.75 1.47 0.68 0.01 0 0 0 2.25 1.43 VỊ TRÍ SE S 1.28 1.6 0.25 0 0 0 0 1.28 1.85 VỊ TRÍ SE S 0.73 0.49 0 0 0 0 0 0.73 0.49 156 SW 0.04 0 0 0.04 W 0.04 0 0 0.04 NW 0.33 0 0 0.33 Tổng 41.79 34.96 16.21 3.98 1.8 1.26 100 SW 19.06 13.54 5.52 0.7 0.08 0.02 38.92 W 0.07 0 0 0.07 NW 0.07 0 0 0.07 Tổng 50.26 30.34 12.2 4.27 2.28 0.65 100 SW 9.65 9.42 0.32 0 19.39 W 4.67 0.13 0 0 4.8 NW 1.39 0 0 1.39 Tổng 56.18 43.13 0.69 0 100 SW 3.28 0.06 0 0 3.34 W 2.82 0.61 0 0 3.42 NW 3.77 0.81 0 0 4.58 Tổng 77.83 22.15 0.02 0 100 SW 0.38 1.24 1.03 0.65 0.14 3.45 W 0.55 2.25 2.06 2.52 3.09 3.22 13.68 NW 0.89 0.21 0 0 1.11 Tổng 5.9 13.08 19.95 27.09 21.3 12.68 100 SW 3.21 6.4 7.8 0.81 0 18.21 W 0.07 0 0 0.07 NW 0.01 0 0 0.01 Tổng 11.47 20.71 33.83 22.82 9.58 1.58 100 SW 2.54 2.48 2.6 0.29 0 7.92 W 2.83 3.83 5.5 1.2 0.06 13.42 NW 0.59 0 0 0.59 Tổng 19.61 43.71 34.25 2.37 0.06 100 [...]... biển Cửa Lấp; - Đề xuất được giải pháp KH-CN khả thi cho việc ổn định luồng tàu, bảo vệ và tôn tạo bờ biển đang bị sạt lở ở khu vực Cửa Lấp 3 Đối tượng và phạm vi nghiên cứu - Đối tượng nghiên cứu: Vùng biển và cửa sông khu vực Cửa Lấp; - Phạm vi nghiên cứu: Phạm vi nội dung: Tập trung vào hai vấn đề liên quan đến sạt lở bờ biển và bồi lấp luồng tàu; phạm vi không gian: Luồng tàu từ cảng cá Phước Tỉnh. .. tàu; Phạm vi nghiên cứu hạn chế trong khu vực Cửa Lấp 1.9 Kết luận chương 1 Chương I trình bày tổng quan một số tình hình nghiên cứu về cửa sông trên thế giới và trong nước, qua đó tìm vấn đề nghiên cứu cho luận án Các thành tựu nghiên cứu được trình bày qua các vấn đề: Phân loại cửa sông, động lực học và diễn biến cửa sông, công trình chỉnh trị cửa sông và phương pháp nghiên cứu cửa sông Nội dung thảo... Các nghiên cứu về công trình chỉnh trị cửa sông ở Cửa Lấp Trong đề tài KH-CN cấp nhà nước của Trương Văn Bốn (2013), các giải pháp KHCN chỉnh trị tổng thể chống bồi, xói để ổn định khu vực Cửa Lấp đã được đề xuất: Nghiên cứu đã kiến nghị và đưa ra được 2 phương án chỉnh trị nhằm thích nghi với hệ thống tự nhiên của khu vực và khả năng đáp ứng của con người Nâng cao sự ổn định của đoạn 19 bờ có nguy cơ. .. chỉnh trị Tháng 5/2013, Trương Văn Bốn đã bảo vệ thành công đề tài KH-CN cấp nhà nước: Nghiên cứu nguyên nhân, cơ chế diễn biến hình thái và đề xuất các giải pháp KHCN nhằm ổn định vùng cửa biển Lộc An, Cửa Lấp tỉnh Bà Rịa - Vũng Tàu [48] Đề tài đã giải quyết được một số nội dung liên quan đến khu vực Cửa Lấp như sau: - Nghiên cứu tính toán dòng chảy ven bờ: Kết quả tính toán tại sườn triều lên và triều... mỗi cửa sông đều là một trường hợp cá biệt cần nghiên cứu, không thể áp dụng những giải pháp chung giống nhau Riêng Cửa Lấp, một cửa ngõ quan trọng của thành phố Vũng Tàu, trải qua quá trình nghiên cứu khoảng 15 năm (tính từ năm 2000, Trương Đình Hiển tiến hành nghiên cứu đề tài "Nghiên cứu các điều kiện động lực học khu vực cửa Lộc An và phụ cận nhằm phục vụ quy hoạch phát triển kinh tế xã hội tỉnh Bà. .. vấn đề ô nhiễm môi trường vùng cửa sông được quan tâm hơn cả, bởi vì sự ô nhiễm ở vùng này thường gây tác hại nặng nề, dễ lan truyền đến các vùng miền khác và rất khó giải quyết Trong phần tổng quan này, luận án đề cập đến bốn vấn đề cơ bản về nghiên cứu cửa sông: Phân loại cửa sông, động lực học cửa sông, công trình chỉnh trị cửa sông và phương pháp nghiên cứu về cửa sông 1.2 Tổng quan các nghiên cứu. .. nhiều vấn đề phải xem xét, nghiên cứu; - Các nghiên cứu đã chỉ ra sự phức tạp của thủy động lực khu vực Cửa Lấp, tuy nhiên kết quả nghiên cứu chưa phân tích được vai trò của từng yếu tố và mối tương tác giữa sóng và dòng chảy trong cơ chế chung về xói bồi của từng khu vực Kết quả nghiên cứu mới chỉ đưa ra được bản đồ minh họa cho sự bồi xói trong một tháng trong thời kỳ gió mùa Đông Bắc và gió mùa... hội tỉnh Bà Rịa - Vũng Tàu" ) [1], đã có một số đề tài KH-CN, dự án đầu tư tuy đã được nghiệm thu, nhưng các kết quả nghiên cứu vẫn 1 còn nhiều ý kiến chưa đồng thuận Giải pháp được Sở Khoa học và Công nghệ tỉnh Bà Rịa - Vũng Tàu đánh giá là thành công nhất khi chỉnh trị bờ biển khu vực cửa Lộc An (năm 2005) là sử dụng công nghệ stabiplage của Pháp cũng chưa thể áp dụng rộng ra cho khu vực khác Việc... cứu về phân loại cửa sông Những đề tài nghiên cứu cơ bản về cửa sông trước hết là nghiên cứu sự hình thành và phát triển cửa sông, châu thổ cửa sông, phân loại và phân đoạn cửa sông Trong luận án này quan tâm đến vấn đề phân loại cửa sông 1.2.1 Các nghiên cứu về phân loại cửa sông trên thế giới Trong giáo trình “Động lực học và công trình cửa sông” của Lương Phương Hậu (2005) [5] và luận văn thạc sĩ... án này lấy Cửa Lấp làm đối tượng nghiên cứu góp phần giải quyết một số vấn đề đang có tính thời sự cấp thiết nêu trên 7 Cấu trúc của luận án Ngoài phần mở đầu, kết luận, nội dung luận án gồm 4 chương: Chương 1: Tổng quan tình hình nghiên cứu về cửa sông Chương 2: Cơ sở lựa chọn mô hình thủy thạch động lực cho khu vực Cửa Lấp Chương 3: Nghiên cứu chế độ thủy động lực học và bồi, xói khu vực Cửa Lấp Chương . VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN VĂN GIÁP NGHIÊN CỨU CƠ CHẾ BỒI, XÓI VÀ ĐỀ XUẤT GIẢI PHÁP CHỈNH TRỊ KHU VỰC CỬA LẤP TỈNH BÀ RỊA - VŨNG. luồng tàu, bảo vệ và tôn tạo bờ biển đang bị sạt lở ở khu vực Cửa Lấp. 3. Đối tượng và phạm vi nghiên cứu - Đối tượng nghiên cứu: Vùng biển và cửa sông khu vực Cửa Lấp; - Phạm vi nghiên cứu: . bốn vấn đề cơ bản về nghiên cứu cửa sông: Phân loại cửa sông, động lực học cửa sông, công trình chỉnh trị cửa sông và phương pháp nghiên cứu về cửa sông. 1.2 Tổng quan các nghiên cứu về phân