Nghiên cứu xác định đáy chạy tàu hợp lý trong trường hợp có bùn loãng tại một số luồng hàng hải ở việt nam

BỘ GIÁO DỤC & ĐÀO TẠO BỘ NÔNG NGHIỆP & PT NÔNG THÔN VIỆN KHOA HỌC THUỶ LỢI VIỆT NAM NGUYỄN ANH TUẤN NGHIÊN CỨU XÁC ĐỊNH ĐÁY CHẠY TÀU HỢP LÝ TRONG TRƯỜNG HỢP CĨ BÙN LỖNG TẠI MỘT SỐ LUỒNG HÀNG HẢI Ở VIỆT NAM LUẬN ÁN TIẾN SĨ CHUN NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH THUỶ HÀ NỘI, NĂM 2020 BỘ GIÁO DỤC & ĐÀO TẠO BỘ NÔNG NGHIỆP & PT NÔNG THÔN VIỆN KHOA HỌC THUỶ LỢI VIỆT NAM ĐỀ TÀI LUẬN ÁN NGHIÊN CỨU XÁC ĐỊNH ĐÁY CHẠY TÀU HỢP LÝ TRONG TRƯỜNG HỢP CĨ BÙN LỖNG TẠI MỘT SỐ LUỒNG HÀNG HẢI Ở VIỆT NAM CHUN NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH THUỶ MÃ SỐ: 09-58-02-02 CÁN BỘ HƯỚNG DẪN NGHIÊN CỨU SINH Cán hướng dẫn GS TSKH Nguyễn Ngọc Huệ Cán hướng dẫn PGS TS Nguyễn Khắc Nghĩa Hà Nội, 2020 Nguyễn Anh Tuấn i LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng tôi, kết nghiên cứu đƣợc trình bày luận án trung thực, khách quan chƣa đƣợc bảo vệ học vị Tơi xin cam đoan thơng tin trích dẫn luận án đƣợc rõ nguồn gốc Hà Nội, ngày Tác giả luận án Nguyễn Anh Tuấn ii LỜI CÁM ƠN Tác giả xin chân thành cám ơn Viện Khoa học Thủy lợi Việt Nam đào tạo, hƣớng dẫn tạo điều kiện giúp đỡ tác giả suốt trình thực luận án Chân thành cảm ơn Cục Hàng hải Việt Nam gia đình tạo điều kiện thời gian hội tiếp cận, sử dụng nhiều nguồn số liệu phục vụ cho trình thực luận án Tác giả xin bày tỏ lòng biết ơn sâu sắc đến thầy hƣớng dẫn GS.TSKH Nguyễn Ngọc Huệ PGS.TS Nguyễn Khắc Nghĩa tận tình hƣớng dẫn, định hƣớng giúp tác giả hoàn thành luận án iii MỤC LỤC LỜI CAM ĐOAN LỜI CÁM ƠN DANH MỤC HÌNH VẼ DANH MỤC BẢNG BIỂU DANH MỤC CÁC TỪ VIẾT TẮT MỞ ĐẦU CHƯƠNG 1: TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU VỀ CHẠY TÀU TRONG LUỒNG CĨ BÙN LỖNG 1.1 Tổng quan nghiên cứu bùn loãng 1.1.1 Các vấn đề nghiên cứu 1.1.2 Phƣơng pháp nghiên cứu 1.2 Một số kết nghiên cứu giới 1.2.1 Nghiên cứu bùn loãng 1.2.2 Nghiên cứu đáy chạy tầu độ sâu chạy tầu 1.2.3 Nghiên cứu tính khả thi chạy tầu bùn loãng 1.2.4 Các ứng dụng xử lý lớp bùn loãng 1.3 Kết nghiên cứu Việt Nam 1.3.1 Kết nghiên cứu ứng dụng luồng Nam Triệu, Hải Phịng 1.3.2 Kết nghiên cứu luồng Sồi Rạp 1.3.3 Kết nghiên cứu luồng tầu nhà máy nhiệt điện Duyên Hải 1.3.4 Kết nghiên cứu luồng Sơng Hậu 1.4 Phân tích chung vấn đề nghiên cứu Việt Nam 1.5 Định hƣớng nghiên cứu luận án Kết luận chƣơng CHƯƠNG 2: CƠ SỞ LÝ THUYẾT, PHƯƠNG PHÁP NGHIÊN CỨU XÁC ĐỊNH LỚP BÙN LỖNG VÀ ĐÁY CHẠY TẦU TRONG LUỒNG CĨ BÙN LOÃNG 2.1 Cơ sở lý thuyết vấn đề nghiên cứu 2.1.1 Cơ sở lý thuyết chuyển động bùn dƣới tác động yếu tố động lực cửa sông ven biển 2.1.2 Cơ sở lý thuyết phép phân tích Fourier 2.1.3 Cơ sở lý thuyết mơ hình hồi quy đa biến 2.2 Cơ sở số liệu phục vụ cho nghiên cứu 2.2.1 Số liệu địa hình 2.2.2 Số liệu thủy hải văn, bùn cát iv 2.2.3 Ảnh viễn thám 2.3 Các phƣơng pháp nghiên cứu sử dụng luận án 2.3.1 Phƣơng pháp khảo sát trƣờng 2.3.2 Phƣơng pháp phân tích Fourier 2.3.3 Phƣơng pháp phân tích mơ hình hồi quy đa biến xây dựng tƣơng quan chiều dày lớp bùn loãng số yếu tố thủy hải văn, bùn cát Kết luận chƣơng CHƯƠNG 3: ĐIỀU KIỆN ÁP DỤNG CHẠY TẦU TRÊN BÙN LỖNG, PHÂN TÍCH FOURIER VÀ MƠ HÌNH HỒI QUY ĐA BIẾN XÁC ĐỊNH QUY LUẬT SA BỒI VÀ CHIỀU DÀY LỚP BÙN LOÃNG 3.1 Kết nghiên cứu điều kiện sử dụng công nghệ chạy tầu bùn loãng luồng tầu biển Việt Nam 3.1.1 Điều kiện khai thác, vận tải 3.1.2 Điều kiện tự nhiên 3.1.3 Thực tế bùn loãng cao trình đáy luồng số luồng hàng hải 3.2 Xây dựng phƣơng trình xác định chiều dày lớp bùn lỗng mơ hình hồi quy đa biến 3.2.1 Xác định tham số đầu vào 3.2.2 Tính tốn xây dựng công thức đặc trƣng xác định chiều dày lớp bùn lỗng mơ hình hồi quy đa biến 3.3 Xây dựng mơ hình tính tốn biến động sa bồi luồng tầu theo tháng 3.3.1 Dữ liệu đầu vào 3.3.2 Chƣơng trình tự động hóa tính tốn 3.3.3 Quy trình tính tốn 3.3.4 Kết xây dựng hàm chiều dày sa bồi theo tháng cho luồng Bạch Đằng Kết luận chƣơng CHƯƠNG 4: ỨNG DỤNG KẾT QUẢ NGHIÊN CỨU 4.1 Ứng dụng cơng thức đặc trƣng tính tốn độ dày lớp bùn loãng xác định đáy chạy tầu đáy nạo vét hợp lý trƣờng hợp có bùn lỗng luồng Sồi Rạp – Hiệp Phƣớc 4.1.1 Tính tốn thử nghiệm dự báo chiều dày lớp bùn loãng 4.1.2 Xác định đáy chạy tầu thiết kế đáy nạo vét 4.2 Xác định thời điểm nạo vét hợp lý theo tiêu chí hiệu khai v thác luồng Bạch Đằng Kết luận chƣơng KẾT LUẬN DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ TÀI LIỆU THAM KHẢO PHỤ LỤC vi Hình 1.1 Băng đo sâu tần số 210 kHz 33 kHz …………………… Hình 1.2 Kết mơ tƣơng tác sóng bùn lỗng ………… …… Hình 1.3 Tính chất lƣu biến bùn lỗng ……………… …………… Hình 1.4 Độ sâu luồng đƣợc đo thiết bị hồi âm đa tần (h.vẽ minh họa) Hình 1.5 Cảng Savannah vùng lân cận, bẫy bùn cát cửa ngăn triều … Hình 1.6 Giải pháp cơng trình cửa sơng Tapi: (a) phƣơng án kết hợp đê chìm hƣớng dịng dọc tuyến luồng; (b) phƣơng án đào bẫy bùn vũng quay tầu Hình 1.7 Vị trí lấy mẫu bùn cát lỗng tuyến luồng Sồi Rạp Hình 1.8 Luồng tầu nhà máy nhiệt điện Duyên Hải ………………….… Hình 2.1 Sơ đồ biến động lớp bùn lỗng …………………………… Hình 2.2 Biểu đồ mô tả khởi động bứt phá bùn lỗng ……………… Hình 2.3 Biểu đồ quan hệ S  Hình 2.4 Phân bố nồng độ kết bơng đƣợc lựa chọn ……………… Hình 2.5 Tính tốn tốc độ chìm lắng bị cản trở ………………………… Hình 2.6 Các trình vận chuyển bùn ….………………………… Hình 2.7 Kỹ thuật đo sâu hồi âm (hình vẽ minh họa) ……………………… Hình 2.8 Sơ đồ trình phân tích Fourier …… ……………………… Hình 2.9 Lƣợc đồ quy trình thiết lập phƣơng trình dự báo Hình 3.1 Cao độ đáy luồng tần số 33kHz 200kHz – Đoạn luồng Lạch Huyện … ……………………………………….……………… Hình 3.2 Cao độ đáy luồng tần số 33kHz 200kHz – Đoạn luồng Soài Rạp ……………………………………………………….……… Hình 3.3 Cao độ đáy luồng tần số 33kHz 200kHz – Đoạn luồng Quan Chánh Bố ……………………………………………………… Hình 3.4 Cao độ đáy luồng tần số 33kHz 200kHz – Đoạn luồng Vũng Tầu Thị Vải …………………………………………… ……… vii Hình 3.5 Bản đồ phân bố bùn lỏng tồn tuyến 2018………………………… Hình 3.6 Bản đồ phân bố bùn lỏng (mặt cắt 30-42) năm 2016-2017……… Hình 3.7 Bản đồ phân bố bùn lỏng (mặt cắt 30-42) năm 2018……….…… Hình 3.8 Lựa chọn khu vực lựa chọn tính tốn …………………………… Hình 3.9 Xây dựng mơ hình khu vực tính tốn (a) địa hình lịng dẫn; (b) lƣới tính; (c) điều kiện biên Hình 3.10 Kiểm định mực nƣớc, theo mực nƣớc theo số liệu thực đo …… Hình 3.11 Kiểm định vận tốc dịng chảy ………………………………… Hình 3.12 Kiểm định hƣớng dịng chảy … … ………………………… Hình 3.13 Đồng liệu hệ thống lƣới trực giao …… … Hình 3.14 Sơ đồ trình định dạng đồng liệu ……………………… Hình 3.15 Sơ đồ tính tốn xây dựng hàm hồi quy đa biến ………………… Hình 3.16 Sơ đồ rút gọn quy trình tính tốn hồi quy đa biến ……………… Hình 3.17 Scatter plot giá trị thu đƣợc từ phƣơng trình hồi quy giá trị quan trắc thực sử dụng nhân tố dự báo…………………………… Hình 3.18 So sánh kết tính tốn độ dày lớp bùn lỗng hàm hồi quy biến, biến biến, luồng Sồi Rạp Hình 3.19 Sơ đồ thuật tốn mơ hình tính tốn sa bồi tháng (SBThang) …… Hình 3.20 Sơ đồ giản lƣợc q trình tính tốn ……………………………… Hình 3.21 Biểu đồ độ dày sa bồi trung bình luồng Bạch Đằng…………… Hình 3.22 So sánh kết dự báo sa bồi luồng Bạch Đằng với chuỗi trung bình tháng sử dụng công thức BĐ1-BĐ3 hàm hồi quy bậc … …… Hình 3.23 Hình 3.20 So sánh lựa chọn công thức theo hệ số tƣơng quan R …… ……………………………………………………………… Hình 4.1 Đoạn luồng lựa chọn tính tốn thử nghiệm (Phao 49 đến phao 58 Luồng Sồi Rạp) …………………………………………………………… 109 Hình 4.2 Trƣờng vận tốc tính tốn trung bình năm (từ P49 đến P58) …… 110 Hình 4.3 Trƣờng hệ số kéo tính tốn …………………………… ……… 110 Hình 4.4 Trƣờng SPM (độ đục) trung bình đƣợc tính tốn từ ảnh Landsat8 111 Hình 4.5 Phân bố độ dày trung bình lớp bùn lỗng ………………………… 111 viii Hình 4.6 Độ dày lớp bùn lỗng đoạn từ P54 đến đoạn P56 ………………… 112 Hình 4.7 Độ dày lớp bùn loãng cho đoạn luồng từ phao 49 đến phao 54…… 112 Hình 4.8 Sơ đồ mật cắt dọc tuyến nghiên cứu thử nghiệm gồm 79 điểm … 113 Hình 4.9 Biểu đồ tổng giá trị hcl tổ chức nạo vét vào tháng khác luồng Bạch Đằng ………………… ………………………………… 122 endif enddo enddo i=n1,n+ng write(10,13) ns(i)-c(k1)*cos(2*pi*k1*(i-phi(k1))/n) enddo endif 13 format (3x,f12.4) 11 format (10x, a21) Endprogram II Giao diện chạy chương trình: Ảnh hình vẽ tả bác bước chạy Bước – Nhập file đầu vào Bước – Hiển thị thông số để xác định hàm lọc Bước – Lựa chọn hàm lọc ghi kết 28 PHỤ LỤC Chương trình hồi quy đa biến xác định chiều dày bùn lỗng (HQĐB) I Code chương trình PROGRAM HQĐB PARAMETER (NMAX=20000,MMAX=20,NMmax=NMAX*MMAX) REAL X(NMAX,MMAX),X1(NMmax) INTEGER IDX(MMAX) Integer Mbien CHARACTER*20 FN,FIN,FN1 Character*8 Ten(MMAX) Common Ten WRITE(*,'(A)')' CHO TEN FILE SO LIEU: ' READ(*,*) FN WRITE(*,'(A)')' CHO TEN FILE KET QUA: ' READ(*,*) FIN OPEN(1,FILE=FN,Status='Old') Open(3,File=FIn,Status='unknown') Read(1,*) Mbien Read(1,*) (Idx(i),i=1,Mbien) Read(1,777) (Ten(i),i=1,Mbien) 777 format(40A8) I=1 READ(1,*,END=100)(X(I,J),J=1,Mbien) I=I+1 IF(I-NMax) 5,5,100 100 N=I-1 if(n.le.Mbien+2)goto 2110 L=0 DO 160 J=1,Mbien DO 160 I=1,N L=L+1 X1(L)=X(I,J) 160 CONTINUE WRITE(*,'(A)')' CHO HANG SO EPSILON: ' READ(*,*) PRC WRITE(*,*)' EPSILON = ',PRC CALL HQ(N,Mbien,1,PRC,X1,IDX) 2110 Continue END SUBROUTINE HQ(N,M,NS,PCT,X,IDX) PARAMETER(NMAX=1000,MMAX=20,NMmax=NMAX*MMAX) C DIMENSION X(NMmax),XBAR(MMAX),STD(MMAX),D(MMAX) DIMENSION B(MMAX),IDX(MMAX),L(MMAX) DIMENSION RX(NMmax),R(NMmax),NSTEP(5),ANS(11),T(MMAX) Character *8 Ten(40) Common Ten Write(3,500) Ten(1) 500 Format(//10X,' HOI QUI GIUA ',A,' VA CAC YEU TO') Write(3,501) (Ten(i),i=2,M) 501 Format(5X,40A) 29 WRITE(3,2) N,M FORMAT(20X,'Hoi quy tung buoc') FORMAT(20X,'Dung luong mau',I5/10X,'So bien',7X,I5) IO=1 CALL CORRE(N,M,IO,X,XBAR,STD,RX,R,B,D,T) WRITE (3,5) FORMAT(/10X,'Ten bien',3x,'So TT bien',4X,'Trung binh', *4X,'Do lenh chuan') DO 105 I=1,M 105 WRITE(3,6) Ten(i),I,XBAR(I),STD(I) 6FORMAT (10X,A,5x,I4,4X,F10.2,6X,F10.2) IF(NS)135,135,140 135 WRITE(3,9) 9FORMAT(/20X,'So chon NS khong xac dinh CHAM DUT cong viec') GOTO 200 140 CALL MSTR(RX,R,M,0,1) NSEL=1 GOTO 150 145 CALL MSTR(R,RX,M,1,0) 150 WRITE(3,10)NSEL 10 FORMAT(//10X,'Lan chon thu',I10) N35=0 DO 155 K=1,M IF (IDX(K))152,153,153 152 WRITE(3,11)K,NSEL 11 FORMAT(20X,'Chi so bien thu',I5,2X,'Am o lan chon thu',I5) GOTO 185 153 IF(IDX(K)-3)155,154,152 154 N35=N35+1 155 CONTINUE IF (N35-1)156,157,156 156 WRITE(3,12)NSEL 12 FORMAT(20X,'Lan chon thu',I4,2X,'So bien phu thuoc khac 1') GOTO 185 157 CALL STEP(M,N,RX,XBAR,IDX,PCT,NSTEP,ANS,L,B,STD,TT,D,IER) 185 IF(NSEL-NS)190,200,200 190 NSEL= NSEL+1 GOTO 145 200 CO NTINUE return END SUBROUTINE CORRE(N,M,I0,X,XBAR,STD,RX,R,B,D,T) DIMENSION X(1),XBAR(1),STD(1),RX(1),R(1),B(1),D(1),T(1) DO 100 J=1,M B(J)=0 100 T(J)=0 K=(M*M+M)/2 DO 102 I=1,K 102 R(I)=0 FN=N L=0 IF(I0)105,127,105 30 105 DO 108 J=1,M DO 107 I=1,N L=L+1 107 T(J)=T(J)+X(L) XBAR(J)=T(J) 108 T(J)= T(J)/FN DO 115 I=1,N JK=0 L=I-N DO 110 J=1,M L=L+N D(J)=X(L)-T(J) 110 B(J)=B( J)+D(J) DO 115 J=1,M DO 115 K=1,J JK=JK+1 115 R(JK)=R(JK) +D(J)*D(K) GOTO 205 127 IF(N-M)130,130,135 130 K K=N GOTO 137 135 KK=M 137 DO 140 I=1,KK CALL DATA(M,D) DO 140 J=1,M T(J)=T(J)+D(J) L=L+1 140 RX(L)=D(J) FRR=KK DO 150 J=1,M XBAR(J)=T(J) 150 T(J)=T(J)/FRR L=0 DO 180 I=1,KK JK=0 DO 170 J=1,M L=L+1 170 D(J)=RX( L)-T(J) DO 180 J=1,M B(J)=B(J) +D(J) DO 180 K=1,J JK=JK+1 180 R(JK)=R(JK) +D(J)*D(K) IF(N-KK) 205,205,185 185 KK= N-KK DO 200 I=1,KK JK=0 CALL DATA(M,D) DO 190 J=1,M XBAR(J)=XBAR(J)+D(J) D(J)=D(J)-T(J) 190 B(J)=B(J)+D(J) 31 DO 200 J=1,M DO 200 K=1,M JK=JK+1 200 R(JK)=R(JK)+D(J)*D(K) 205 JK=0 DO 210 J=1,M XBAR(J)=XBAR(J)/FN DO 210 K=1,J JK=JK+1 210 R(JK)=R(JK)B(J)*B(K)/FN JK=0 DO 220 J=1,M JK=JK+J 220 STD(J)=SQRT(ABS( R(JK))) DO 230 J=1,M DO 230 K=J,M JK=J+(K*KK)/2 L=M*(J-1)+K RX(L)=R(JK) L=M*(K-1)+J RX(L)=R(JK) IF(STD(J)STD(K)) 225,222,225 222 R (JK)=1 GOTO 230 225 STD2=STD(J)*STD(K) R(JK)=R(JK)/STD2 C 225 R(JK)=R(JK)/(STD(J)*STD(K)) C WRITE(3,500) R(JK) C500 FORMAT(5X,F10.4) 230 CONTINUE FN=SQRT(FN-1.) DO 240 J=1,M 240 STD(J)=STD(J)/FN L=-M DO 250 I=1,M L=L+M+1 250 B(I)=RX(L) RETUR N END SUBROUTINE DATA(M,D) RETURN END SUBROUTINE MSTR(A,R,N,MSA,MSR) DIMENSION A(1),R(1) DO 20 I=1,N DO 20 J=1,N IF(MSR)5,10,5 IF(I-J)10,10,20 10 CALL EXTERNAL LOC(I,J,IR,N,N,MSR) IF (IR)20,20,15 15 R(IR)=0 CALL EXTERNAL LOC(I,J,IA,N,N,MSA) 32 IF(IA)20,20,18 18 R(IR)=A(IA) 20 CO NTINUE RETURN END SUBROUTINE EXTERNAL LOC(I,J,IR,N,M,MS) IX=I JX=J IF (MS-1)10,20,30 10 IRX=N*(J X-1)+IX GOTO 36 20 IF(IX-JX)22,24,24 22 IRX=IX+ (JX*JX-JX)/2 GOTO 36 24 IRX=JX+ (IX*IX-IX)/2 GOTO 36 30 IRX=0 IF(IXJX)36,32,36 32 IRX=IX 36 I R=IRX RETUR N END SUBROUTINE STEP(M,N,D,XBAR,IDX,PCT,NSTEP,ANS,L,B,S,T,LL,IER) DIMENSION D(1),XBAR(1),IDX(1),NSTEP(7),ANS(11),L(1),B(1),S(1) *,T(1),LL(1) IER=0 ONM=N-1 NFO=0 NSTEP(3)=0 ANS(3)=0 ANS(4)=0 NSTOP=0 DO 30 I=1,M LL(I)=1 IF(IDX(I))30,30,10 10 IF(IDX(I)-2)15,20,25 15 NFO=NFO+1 IDX(NFO)=I GOTO 30 20 NSTEP(3)=NST EP(3)+1 LL(I)=-1 GOTO 30 25 MY=I NSTEP(1)=MY LY=M*(MY-1) LYP=LY+MY ANS(5)=D(LYP) 30 CONTINU E NSTEP(2)=NFO MX=M-NSTEP(3)1 DO 140 NL=1,MX RD=0 33 IF(NL-NFO) 35,35,55 35 DO 50 I=1,NFO K=IDX(I) IF(LL(K))50,50,40 40 LYP=LY+K IP=M*(K-1)+K RE=D(LYP)**2/D(IP) IF(RD-RE) 45,50,50 45 RD=RE NEW=K 50 CO NTINUE GOTO 75 55 DO 70 I=1,M IF(I-MY) 60,70,60 60 IF(LL(I))70,70,62 62 LYP=LY+I IP=M*(I-1)+I RE=D(LYP)**2/D(IP) IF(RD-RE) 64,70,70 64 RD=RE NEW=I 70 CONTINUE 75 IF(RD) 77,77,76 76 IF(ANS(5)-(ANS(3)+RD))77,77,78 77 I ER=1 GOTO 150 78 RE=RD/ANS (5) IF(RE-PCT) 150,80,80 80 LL(NEW)=0 L(NL)=NEW ANS(1)=RD ANS(2)=RE ANS(3)=ANS(3)+RD ANS(4)=ANS(4)+RE NSTEP(4)=NL NSTEP(5)=NEW 85 ANS(6)=SQRT(ANS(4)) RD=NL RE=ONM-RD RE=(ANS(5)-ANS(3))/RE ANS(7)=(ANS(3)/RD)/RE C 90 ANS(8)=SQRT(RE) C C SUA LAI CAU LENH TREN BANG DOAN SAU C 90 IF(RE.GE.0.) THEN ANS(8)=SQRT(RE) ELSE ANS(8)=999.0 STOP 'DUNG VI KHONG TINH TIEP DUOC !!!' END IF IP=M*(NEW1)+NEW 34 RD=D(IP) LYP=NEW-M DO 100 J=1,M LYP=LYP+M IF(LL(J))100,94,97 94 IF(J-NEW) 96,98,96 96 IJ=M*(J-1)+J D(IJ)=D(IJ) +D(LYP)*D(LYP)/RD 97 D(LYP)=D(L YP)/RD GOTO 100 98 D(IP)=1/RD 100 CONTINUE LYP=LY+NEW B(NL)=D(LYP) IF(NL-1)112,112,105 105 ID=NL-1 DO 110 J=1,ID IJ=NL-J KK=L(IJ) LYP=LY+KK B(IJ)=D(LYP) DO 110 K=1,J IK=NLK+1 MK=L(IK) LYP=M*(MK-1)+KK 110 B(IJ)=B(IJ)-D(LYP)*B(IK) 112 ANS(9)=XBAR(MY) DO 115 I=1,NL KK=L(I) ANS(9)=ANS(9)B(I)*XBAR(KK) IJ=M*(KK1)+KK 114 S(I)=ANS(8)*SQRT(D(IJ)) 115 T(I)= B(I)/S(I) IP=M*(NEW1) DO 130 I=1,M IJ=I-M IK=NEW-M IP=IP+1 IF(LL(I)) 130,130,120 120 DO 126 J=1,M IJ=IJ+M IK=IK+M IF(LL(J)) 126,122,122 122 IF(J-NEW)124,126,124 124 D(IJ)=D(IJ)-D(IP)*D(IK) 126 CONTIN UE D(IP)=D(IP)/(RD) 130 CONTI NUE RD=NNSTEP(4) RD=ONM/RD 132 ANS(10)=SQRT(1-(1-ANS(6)**2)*RD) 134 IF(RD.GE.0) GOTO 905 35 ANS(11)=999 GOTO 906 905 ANS(11)=ANS(8)*SQRT(RD) 906 CONTINUE CALL INHQ(NSTEP,ANS,L,B,S,T,NSTOP) IF (NSTOP) 140,140,150 140 CONTINUE 150 R ETURN END SUBROUTINE INHQ(NSTEP,ANS,L,B,S,TT,NSTOP) DIMENSION NSTEP(5),ANS(11),L(1),B(1),S(1),TT(1) Character *8 Ten(40) Common Ten WRITE (3,1)NSTEP(4) FORMAT(5X,'Buoc thu ',I5) WRITE(3,4)ANS(6) FORMAT(/5X,'He so T.quan boi',F10.5) WRITE(3,6)ANS(8) FORMAT(5X,'Sai so chuan: ',F10.5) WRITE(3,7) FORMAT(5X,'Bien thu',5X,'He so hoi quy') K=NSTEP(4) DO 20 I=1,K 20 WRITE(3,8)L(I),B(I) FORMAT(5X,I7,6X,F10.5,10X,F15.5,10X,F10.5) WRITE(3,9)ANS(9) FORMAT(5X,'He so tu do',F10.5/) RETUR N END II Một số hình ảnh minh họa chạy chương trình - Bước – Gọi chương trình, nhập tên file đầu vào 36 - Bước – Nhập tên file ghi số liệu đầu Bước 3- Nhập tham số phân tích Bước – In kết phục vụ phân tích 37 ... trì luồng tàu làm sở xác định đáy chạy tàu thời điểm nạo vét hợp lý trường hợp có bùn loãng cho số luồng tàu Việt Nam 36 CHƢƠNG 2: CƠ SỞ LÝ THUYẾT, PHƯƠNG PHÁP NGHIÊN CỨU XÁC ĐỊNH LỚP BÙN LOÃNG... pháp xác định độ dầy lớp bùn loãng, ứng dụng cho số luồng tàu Việt Nam; (3) Xây dựng tiêu hiệu trì luồng tàu làm sở xác định đáy chạy tàu thời điểm nạo vét hợp lý trường hợp có bùn loãng cho số luồng. .. THUỶ LỢI VIỆT NAM ĐỀ TÀI LUẬN ÁN NGHIÊN CỨU XÁC ĐỊNH ĐÁY CHẠY TÀU HỢP LÝ TRONG TRƯỜNG HỢP CĨ BÙN LỖNG TẠI MỘT SỐ LUỒNG HÀNG HẢI Ở VIỆT NAM CHUYÊN NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH THUỶ MÃ SỐ: 09-58-02-02