BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC NHA TRANG NGUYỄN THANH TÍN NGHIÊN CỨU XÁC ĐỊNH QUY MÔ ĐỘI TÀU KHAI THÁC HẢI SẢN PHÙ HỢP VỚI NGHỀ CÁ CỦA VIỆT NAM LUẬN VĂN THẠC SĨ KỸ THUẬT KHÁNH HÒA - NĂM 2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC NHA TRANG NGUYỄN THANH TÍN NGHIÊN CỨU XÁC ĐỊNH QUY MÔ ĐỘI TÀU KHAI THÁC HẢI SẢN PHÙ HỢP VỚI NGHỀ CÁ CỦA VIỆT NAM LUẬN VĂN THẠC SĨ KỸ THUẬT Ngành: Kỹ thuật khí động lực Mã số: 60520116 Quyết định giao đề tài: 502/QĐ-ĐHNT ngày 05/07/2016 Quyết định thành lập HĐ: 294/QĐ-ĐHNT ngày 15/03/2017 Ngày bảo vệ: 22/04/2017 Người hướng dẫn khoa học: PGS.TS Trần Gia Thái Chủ tịch Hội đồng: PGS.TS Phạm Hùng Thắng Khoa sau đại học: KHÁNH HÒA - NĂM 2017 LỜI CAM ĐOAN Tôi xin cam đoan luận văn hoàn thành trình nghiên cứu tài liệu, thực chuyến điều tra khảo sát hoạt động khai thác của tàu khai thác, tàu dịch vụ hậu cần hoạt động thu mua tại vựa hải sản tại cảng cá Phan Rí Cửa, cảng cá thành phố Phan Thiết tỉnh Bình Thuận Số liệu sử dụng luận văn toàn bộ kết quả, khảo sát thống kê của đề tài: “Nghiên cứu xác định quy mô đội tàu khai thác hải sản phù hợp với nghề cá Việt Nam” Số liệu luận văn hoàn toàn trung thực, xử lý theo phương pháp khoa học đảm bảo độ tin cậy Một số số liệu luận văn Lãnh đạo Chi Cục Thủy Sản Bình Thuận cung cấp cho phép sử dụng Kết nghiên cứu của luận văn mới, không trùng lặp với bất luận án bảo vệ học vị có trước Nha Trang, ngày 02 tháng 01 năm 2017 Tác giả luận văn Nguyễn Thanh Tín iii LỜI CẢM ƠN Tôi xin bày tỏ lòng cảm ơn sâu sắc tới PGS.TS Trần Gia Thái người trực tiếp hướng dẫn tạo điều kiện để hoàn thành luận văn Tôi xin trân trọng cảm ơn Lãnh đạo Chi Cục Thủy Sản Bình Thuận, KS Nguyễn Văn Cảnh phòng đăng kiểm cung cấp tạo điều kiện để sử dụng số liệu thực đề tài: “Nghiên cứu xác định quy mô đội tàu khai thác hải sản phù hợp với nghề cá Việt Nam” Xin cảm ơn KS Nguyễn Văn Quang, KS Hà Quốc Việt, KS Trần Anh Tuấn tận tình giúp đỡ chuyến điều tra khảo sát tại Bình Thuận Tôi xin gửi lời cảm ơn tới: Ban giám hiệu, Ban chủ nhiệm khoa sau đại học trường Đại học Nha Trang, Quý thầy khoa Kỹ Thuật Giao Thông truyền đạt kiến thức kinh nghiệm giúp hoàn thành luận văn Nha Trang, ngày 02 tháng 01 năm 2017 Tác giả luận văn Nguyễn Thanh Tín iv MỤC LỤC Lời cam đoan …………………………………… iii Lời cảm ơn iv Mục lục…………………………………………………………………………… …iii Danh mục bảng………………………………………………………………… … vii Danh mục hình …………………………………………………………………… x Trích yếu luận văn…………………………………………………………………… x MỞ ĐẦU…………………………………………………………………… ……… 1 Lý lựa chọn đề tài Ý nghĩa khoa học thực tiễn .3 2.1 Ý nghĩa khoa học 2.2 Ý nghĩa thực tiễn 3 Mục tiêu nghiên cứu 4 Phương pháp nghiên cứu .4 Đối tượng, phạm vi nội dung nghiên cứu .5 Chương TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU .5 1.1 Tình hình nghiên cứu nước .7 1.1.1 Tình hình nghiên cứu nước giới 1.1.2 Tình hình nghiên cứu nước .11 1.2 Tình hình nghiên cứu xây dựng mô hình tổ chức đội tàu .14 1.2.1 Các mô hình tổ chức đội tàu khai thác nước giới 14 1.2.2 Các mô hình tổ chức đội tàu khai thác nước ta 18 1.3 Đánh giá chung định hướng nghiên cứu 31 Chương XÂY DỰNG MÔ HÌNH QUY HOẠCH TỐI ƯU ĐỘI TÀU KHAI THÁC PHÙ HỢP VỚI NGHỀ CÁ VIỆT NAM……………………………….…………… 32 2.1 Mô hình toán tối ưu tổng quát .32 2.1.1 Khái niệm chung 32 2.1.2 Phân loại toán quy hoạch 34 2.2 Mô hình toán tối ưu đội tàu khai thác hải sản 35 2.2.1 Đầu vào – đầu của mô hình 36 v 2.2.2 Xây dựng hàm mục tiêu 49 2.2.3 Các điều kiện ràng buộc 50 2.2.4 Giải thuật xác định phân bố đội tàu khai thác theo khu vực mùa vụ .62 2.3 Mô hình toán tối ưu đội tàu mẹ - .65 2.3.1 Tối ưu hóa quy mô của đội tàu mẹ - .66 2.3.2 Tối ưu hóa lộ trình thu gom sản phẩm của tàu mẹ 76 2.4 Lập trình giải toán tối ưu 80 Chương MÔ HÌNH QUY HOẠCH TỐI ƯU ĐỘI TÀU KHAI THÁC TỈNH BÌNH THUẬN .82 3.1 Thực trạng nghề cá Bình Thuận 82 3.1.1 Đặc điểm ngư trường 82 3.1.2 Đặc điểm nguồn lợi mùa vụ khai thác vùng biển Bình Thuận 83 3.1.3 Năng lực đội tàu khai thác 85 3.1.4 Cơ cấu nghề của đội tàu khai thác xa bờ 90 3.2 Mô hình quy hoạch tối ưu đội tàu khai thác Bình Thuận .94 3.2.1 Xây dựng sở liệu đầu vào của mô hình 94 3.2.2 Kết đầu của mô hình 102 3.2.3 Điều chỉnh cấu của đội tàu khai thác Bình Thuận .110 3.2.4 Hỗ trợ việc định 116 3.3 Mô hình đội tàu mẹ - phù hợp nghề cá Bình Thuận 117 3.3.1 Xây dựng sở liệu đầu vào .117 3.3.2 Kết tối ưu hóa đội tàu mẹ - phù hợp với nghề cá Bình Thuận 122 3.3.3 Tối ưu hóa lộ trình của tàu mẹ .124 Chương KẾT LUẬN VÀ KHUYẾN NGHỊ 128 4.1 Kết luận……………………………………………………………………… 128 4.2 Khuyến nghị………………………………………………… 134 TÀI LIỆU THAM KHẢO 130 PHỤ LỤC 132 vi DANH MỤC CHỮ VIẾT TẮT CPUE: Năng suất khai thác (Catch Per Unit Effort) VND: Đơn vị tiền tệ CV: Đơn vị công suất máy DK1: Ngư trường khai thác nhà giàn DK1 DT: Đội tàu tập hợp tàu nghề khai thác KV: Khu vực khai thác MT: Loài mục tiêu khai thác vii DANH MỤC BẢNG Bảng 1.1: Thông số của một số mẫu tàu đánh cá quy mô công nghiệp 14 Bảng 1.2: Các mô hình tổ chức sản xuất biển Ninh Thuận 29 Bảng 2.1: Mùa gió số lượng loài, họ, giống thủy sản vùng biển Việt Nam 40 Bảng 2.2: Phân bố nguồn lợi vùng biển Việt Nam 42 Bảng 2.3: Tổng hợp chỉ số, biến số thông số của mô hình quy hoạch đội tàu 58 Bảng 2.4: Tổng hợp chỉ số, biến số thông số của mô hình mẹ-con 75 Bảng 3.1: Khả khai thác nhuyễn thể hai mảnh vỏ vùng biển Bình Thuận 84 Bảng 3.2: Tình hình phát triển đội tàu cá Bình Thuận giai đoạn (2013 – 2016) 85 Bảng 3.3: Tình hình phát triển đội tàu công suất từ 90 CV trở lên (2010 – 2016) .86 Bảng 3.4: Cơ cấu đội tàu khai thác có công suất máy 90 CV của tỉnh 92 Bảng 3.5: Phân nhóm đội tàu khai thác thủy sản của tỉnh Bình Thuận 94 Bảng 3.6: Tổng hợp thông số của đội tàu khai thác theo nhóm công suất máy 95 Bảng 3.7: Tổng hợp loài mục tiêu (k), khu vực khai thác (l) của đội tàu (i) 97 Bảng 3.8: Trữ lượng khả khai thác hải sản cho phép vùng biển Nam Bộ 99 Bảng 3.9: Ước tính sản lượng khai thác cho phép khai thác của loài hải sản .101 Bảng 3.10: Tổng số lượng tàu khai thác năm theo nghề theo mùa (chiếc) 103 Bảng 3.11: Bảng tổng hợp kết tối ưu với mục tiêu lợi nhuận lớn nhất 104 Bảng 3.12: Kết phân bố số lượng tàu, chuyến biển theo đội tàu (i), nhóm công suất (j) khai thác loài mục tiêu (k) tại khu vực (l), mùa vụ (m) 105 Bảng 3.13: Tổng sản lượng hải sản khai thác tối ưu theo loài (tấn) 109 Bảng 3.14: Phân bố cấu đội tàu theo nhóm công suất 113 Bảng 3.15: Phân bố cấu đội tàu theo nhóm nghề 114 Bảng 3.16: Kết phân bố lại lao động theo khu vực theo mùa vụ .116 Bảng 3.17: Sản lượng đánh bắt theo ngày CPUE (kg/ngày) của tàu đội tàu mẹ - (i), nhóm công suất (j), khu vực (l), mùa vụ (m) 120 Bảng 3.18: Giá thu mua sản phẩm khai thác tại cảng theo loài (s) (nghìn đồng) 120 Bảng 3.19: Số liệu điều tra thực tế thời gian chuyến biển (ngày) của tàu đội tàu mẹ - tàu 120 Bảng 3.20: Chi phí biến đổi cho ngày hoạt động (triệu đồng) 121 Bảng 3.21: Điều kiện đầu vào của chủ tàu 122 viii Bảng 3.22: Điều kiện đầu vào của thủy thủ .122 Bảng 3.23: Quy mô đội tàu mẹ tối ưu hoạt động theo khu vực mùa vụ 123 Bảng 3.24: Tọa độ của tàu mô hình 124 Bảng 3.25: Khoảng cách tàu mô hình 125 Bảng 3.26: Biến định hành trình di chuyển .125 Bảng 3.27: Trình tự di chuyển của tàu mẹ 126 ix DANH MỤC HÌNH Hình 1.1: Tàu lưới vây quy mô công nghiệp 17 Hình 1.2: Mô hình đội tàu luân phiên 19 Hình 1.3: Sơ đồ hoạt động của đội tàu mẹ - 20 Hình 1.4: Tàu Hải Vương mô hình tàu mẹ - 26 Hình 2.1: Cấu trúc của mô hình mối quan hệ toán học 35 Hình 2.2: Tuyến phân vùng khai thác thủy sản vùng biển Việt Nam .39 Hình 2.3: Cấu trúc của mô hình mối quan hệ toán học 61 Hình 2.4: Sơ đồ thuật toán giải thuật xác định phân bố của đội tàu theo khu vực mối quan hệ yếu tố mô hình quy hoạch 64 Hình 2.5: Sơ đồ thuật toán xác định quy mô của đội tàu khai thác dạng mẹ theo khu vực mùa vụ 67 Hình 2.6: Cơ cấu hoạt động của tàu mẹ tàu đội tàu mẹ - 76 Hình 2.7: Sơ đồ thuật toán tối ưu lộ trình tàu mẹ mô hình đội tàu mẹ - .79 Hình 3.1: Phạm vi ngư trường hoạt động của nghề cá Bình Thuận 82 Hình 3.2: Quá trình phát triển lực khai thác của đội tàu khai thác hải sản Bình Thuận giai đoạn 2010 - 2014 88 Hình 3.3: Biểu đồ phân bố số lượng tàu khai thác hải sản địa phương nghề cá của tỉnh Bình Thuận năm 2012 2014 89 Hình 3.4: Biểu đồ phân bố cấu nghề khai thác của đội tàu cá Bình Thuận .93 Hình 3.5: Biểu đồ phân bố số lượng tàu làm nghề khác từ năm 2013 đến 2016 93 Hình 3.6: Biểu đồ tổng sản lượng hải sản theo khu vực theo mùa vụ 109 Hình 3.7: Biểu đồ so sánh số lượng tàu của đội tàu trước sau chạy mô hình tối ưu đội tàu khai thác Bình Thuận 112 Hình 3.8: Phân bố đội tàu theo nhóm công suất 113 Hình 3.9: Điều chỉnh cấu đội tàu Bình Thuận theo nhóm nghề khai thác .114 Hình 3.10: Điều chỉnh cấu đội tàu khai thác Bình Thuận theo khu vực 115 Hình 3.11: Sự phân bố nguồn lực lao động theo không gian thời gian 116 Hình 3.12: Hình ảnh vị trí tàu mô hình đội tàu mẹ - khảo sát .124 Hình 3.13: Hành trình di chuyển ngắn nhất của tàu mẹ 126 Hình 3.14: Hành trình di chuyển của tàu mẹ theo cách thông thường 127 x LRONG(I,j,k,l,m) =E= DT(I,j,k,l,m)-(CPBD(i,j,k,l,'khathac')*TGKT(i,j,k,l,m) +CPBD(i,j,k,l,'dichuyen')* TGDC(i,j,k,l,m) +CPBD(i,j,k,l,'nambo')*TGNB(i,j,k,l,m)) ; * Phuong trinh xac dinh ty le danh bat, cung nhu tong san luong danh bat *Nang suat khai thac thuc te NANGSUATTT(I,j,k,l,m,S)$(CPUE(I,j,k,l,m,S) AND AL(m)) NSTT(I,j,k,l,m,S) =E= (A(l,m,S)*CPUE(I,j,k,l,m,S))$(TL(l,m,s)AND AL(m)); *Ty le danh bat thuc te TYLEDB(l,m,S)$(TL(l,m,s)AND AL(m)) TLDB(l,m,S)*TLKT*TL(l,m,s) =E= TSLKT(l,m,S); *He so tru luong khai thac HESO(l,m,S)$(TL(l,m,s)AND AL(m)) A(l,m,S) =E= - TLDB(l,m,S)**10 ; ****** Ham dieu kien thuc hien chuyen bien DKTHCB(I,j,k,l,m)$(SCBLN(I,j,k,l,m)and Al(m)) LRONG(I,j,k,l,m)*NLDB(I,j,k,l,m) =G= 0; *** Ham dieu kien dau vao cua chu tau: *Phan chia loi nhuan cua chu tau DOANHTHUCT(I,j) DTCT(I,j) =E= SUM((k,l,m)$(SCBLN(I,j,k,l,m)and Al(m)),(Chutau(i,j,'TLchiaLN')*LRONG(I,j,k,l,m)*NLDB(I,j,k,l,m))); *Chi phi co dinh, bao duong hang nam BAODUONG(I,j) CPHN(I,j) =E= Chutau(i,j,'CP-Codinh')*V3(I,j); * Chenh lech giua loi nhuan va chi phi co dinh CLECHCT(I,j) CLCT(I,j) =E= DTCT(I,j)-CPHN(I,j); *Dieu kien tri chuyen bien hang nam cua chu tau DKCT(I,j) CLTT(I,j)*V3(I,j) =G= 0; *** Ham dieu kien dau vao cua thuy thu: *Phan loi nhuan tra luong cho thuy thu LUONGTT(I,j) LTT(I,j) =E= (SUM((k,l,m)$(SCBLN(I,j,k,l,m) AND AL(m)),Thuythu(i,j,'TLchiaLN')*LRONG(I,j,k,l,m)*NLDB(I,j,k,l,m))); *Muc luong mong doi cua thuy thu LUONGMM(I,j) LMM(I,j) =E= SUM((k,l,m)$(SCBLN(I,j,k,l,m)and Al(m)),(Thuythu(i,j,'TN.ngay')*SLTT(i,j,k,l,m)*(TGDC(i,j,k,l,m)+TGKT(i,j,k,l, m)+TGNB(i,j,k,l,m))*NLDB(I,j,k,l,m)))/1000; *Chenh lech giua muc luong thuc te va luong mong muon CLECHTT(I,j) CLTT(I,j) =E= LTT(I,j)-LMM(I,j); *Dieu kien dau vao cua thuy thu DKTT(I,j) CLTT(I,j)*V3(I,j) =G= 0; *** Dieu kien rang buoc ve che chinh sach Hanche1 SUM((j),V3('Doi1',j))=L= 821; Hanche2 SUM((j),V3('Doi2',j)) =L= 753; *Hanche3 * SUM((i,j),V2(i,j,'VBAC')) =E= SUM((i,j),V2(i,j,'VNAM')); ********************************* * Gioi han tren duoi * ********************************* A.UP(l,m,S)=1; TLDB.UP(l,m,S)=1; MODEL HFMMA /ALL/; ********************************* * Option Definition * ********************************* $iftheni x%mode%==xbook $onecho >conopt.opt solvers conopt knitro minos snopt $offecho option nlp=conopt; OPTION ITERLIM = 1000000000000; OPTION RESLIM = 20000000000000; OPTION DOMLIM = 1000000; OPTION LIMROW = 1000; OPTION Rtmaxv=1.00e+50; HFMMA.OPTFILE = 1; *HFMMA.SCALEOPT = 1; HFMMA.WORKSPACE = 80.0; $endif ***************************************** * Initial Values * ***************************************** *V3.up(I,j) = 1000; *A.L(l,m,S) = 0.9; *NLDB.L(I,j,k,l,m) = 0; *V1.L(I,j,k,l,m) = 0; *SLKT.L(I,j,k,l,m,S) = 1; *LRONG.L(I,j,k,l,m) = 1; DKCT.L(I,J) = 0; DKTT.L(I,J) = 0; *************************************************  Giai lan thu nhat * ************************************************* SOLVE HFMMA USING NLP MAXIMIZING LNCT; ********************************* * Result Export * ********************************* * Define temporary parameter for eliminating values of NLDB * when corresponding VESSEL1 is zero; PARAMETER TTT(I,j,k,l,m) temporary parameter; TTT(I,j,k,l,m) = NLDB.L(I,j,k,l,m); NLDB.L(I,j,k,l,m) = TTT(I,j,k,l,m)$(V1.L(I,j,k,l,m)*(SCBLN(I,j,k,l,m))); *********************************************************  Giai lan thu hai * ********************************************************* V2.FX(I,j,m) = FLOOR(V2.L(I,j,m)); V2.UP(I,j,m) = FLOOR(V2.L(I,j,m)); V3.UP(I,j) = SMAX(m,FLOOR(V2.L(I,j,m))); ********************************************************* * Solve Statement (the 2nd Run) Squeeze=N * ********************************************************* SOLVE HFMMA USING NLP MAXIMIZING LNCT; ********************************* * Xuat ket qua * ********************************* TTT(I,j,k,l,m) = NLDB.L(I,j,k,l,m); NLDB.L(I,j,k,l,m) = TTT(I,j,k,l,m)$(V1.L(I,j,k,l,m)*(SCBLN(I,j,k,l,m))); FILE RESULTB /HFMM2B.TXT/; PUT RESULTB; RESULTB.PC = 5; parameter misc(*); misc("Loinhuan")=LNCT.l; misc("modelstat")=HFMMA.modelstat; Execute_Unload 'tmp.gdx',NLDB,V1,V2,V3,TSLKT,LTT,DTCT,A,NSTT,TLDB,TSL16,misc,DT; Execute 'GDXXRW.EXE tmp.gdx Squeeze=N O=Qhchung.xls par=misc rng=OUT!A1 Rdim=1 Cdim=0 var=V1 rng=OUT!A6 Rdim=4 Cdim=1 var=V2 rng=OUT!I2 Rdim=2 Cdim=1 var=V3 rng=OUT!O2 Rdim=2 Cdim=0 var=NLDB rng=OUT!S2 Rdim=4 Cdim=1 var=TSLKT rng=OUT!AA2 Rdim=2 Cdim=1 var=DT rng=OUT!AS2 Rdim=4 Cdim=1'; Phần 2: CHƯƠNG TRÌNH TỐI ƯU MÔ HÌNH SẢN SUẤT ĐỘI TÀU MẸ CON PHÙ HỢP VỚI ĐIỀU KIỆN KHAI THÁC CỦA TỈNH BÌNH THUẬN TRÊN GAMS ********************************* * Quy mo doi tau me con* ********************************* * Lua chon khu vuc khai thac $call =ask T=checklistbox M="Lua chon khu vuc khai thac" D="KV2|KV3" E="KV2|KV3" R="%s ' checked list box choice'" O=k.inc set k / $include k.inc/; display k; * Lua chon mua vu khai thac $call =ask T=checklistbox M="Lua chon mua vu khai thac" D="Vbac|Vnam" E="Vbac|Vnam" R="%s 'checked list box choice'" O=l.inc set l / $include l.inc/; display l; sets s Loai khai thac /L1*L14/ i Loai tau tham gia doi /Tme,Tcon/; parameter a He so suy giam nguon loi; a = 0.01; parameter TL Ty le phan chia loi nhuan giua chu tau va thuy thu; TL = 0.5; parameter TTRC Tai dang ki cua cac tau tham gia mo hinh (tan); TTRC = 25; parameter LCB Muc luong co ban cho lam viec cua thuy thu (Trieu dong); LCB = 0.3; parameter CP(i,*) chi phi bien doi cho mot hoat dong (trieu dong) *CP('Tme','CPNB') chi phi cho neo dau tai cang cua tau me (trieu dong) *CP('Tme','CPDC') chi phi cho di chuyen chuyen bien (trieu dong) *CP('Tme','CPGT') chi phi trao doi giua tau me va tau (trieu dong) *CP('Tme','CPCD') chi phi co dinh cua tau me nam (trieu dong) *CP('Tcon','CPNB') chi phi cho neo dau tai cang cua tau me (trieu dong) *CP('Tcon','CPDC') chi phi cho di chuyen chuyen bien (trieu dong) *CP('Tcon','CPGT') chi phi cho danh bat chuyen bien (trieu dong) *CP('Tcon','CPCD') chi phi co dinh cua tau nam (trieu dong) TG(i,k,l,*) Thoi gian hoat dong cua cac tau doi (Ngay) *TG('Tme',k,l,'TGNB') Thoi gian nam bo cua tau me giua cac chuyen bien (Ngay) *TG('Tme',k,l,'TGDC') Thoi gian di chuyen cua tau me chuyen bien (Ngay) *TG('Tme',k,l,'TGGT') Thoi gian giao tiep cua tau me voi tau doi (Ngay) *TG('Tcon',k,l,'TGNB') Thoi gian nam bo cua tau giua cac chuyen bien (Ngay) *TG('Tcon',k,l,'TGDC') Thoi gian di chuyen cua tau chuyen bien (Ngay) SLTT(i) So luong thuy thu tren moi tau (Nguoi) GIA(l,s) Gia thu mua ca (Nghin dong) CPUE(k,l,s) San luong danh bat cua loai S khu vuc K mua L tren khai thac (kg); $CALL GDXXRW Mcon.xls Cmerge=2 par=SLTT rng=Data!A25:B26 Rdim=1 par=CPUE rng=Data!A31:P35 Rdim=2 Cdim=1 par=CP rng=Data!A4:F6 Rdim=1 Cdim=1 par=TG rng=Data!A10:F18 Rdim=3 Cdim=1 par=GIA rng=Data!A42:O44 Rdim=1 Cdim=1; $GDXIN Mcon.gdx $LOAD CP,TG,CPUE,GIA,SLTT $GDXIN; display CP,TG,CPUE,GIA,SLTT; ********************************* * Khai bao cac bien * ********************************* VARIABLES LNCT Phan loi nhuan cua chu tau mot chuyen bien (Trieu dong) V1(k,l) # So luong tau me doi o khu vuc K mua L (Chiec) V2(k,l) # So luong tau doi o khu vuc K mua L (Chiec) V3(k,l) # Tong so luong tau toan doi (Chiec) CBM(k,l) Chieu dai chuyen bien cua tau me (Ngay) TTRM(k,l) Tai cua tau me(tan) TGTM(k,l) Thoi gian giua hai lan thu mua cua tau me (Ngay) SCBM(k,l) So luong chuyen bien cua tau me lan khai thac (Chuyen) CBC(k,l) Chieu dai chuyen bien cua tau (Ngay) TGKT(k,l) Thoi gian danh bat ca cua tau chuyen bien (Ngay) SLKT(k,l,s) San luong danh bat mot chuyen bien theo loai khu vuc K va mua L (Tan) DT(k,l) Tong doanh thu tren mot chuyen bien(Trieu dong) CPBD(k,l) Tong chi phi bien doi cho mot chuyen bien(Trieu dong) CPCD Tong chi phi bien doi cho mot chuyen bien(Trieu dong) LRONG Lai rong tren chuyen bien ((Trieu dong) DTCT(k,l) Phan chia loi nhuan cua chu tau chuyen bien(Trieu dong) CPHN(k,l) Chi phi co dinh co doi tau I chuyen bien:((Trieu dong)) CLCT(k,l) Loi nhuan cua chu tau sau tru chi phi co dinh:((Trieu dong)) LMM(k,l) Tong muc luong mong muon cua thuy thu tren doi tau i:(Trieu dong) LTT Muc luong cua thuy thu tren doi tau I:(Trieu dong) HS(k,l) Muc luong vuot qua du kien cua thuy thu tren doi tau I:(Trieu dong); POSITIVE VARIABLES V1, V2, V3, TQ, SLKT, REVENUE, TGTM, SCBM; POSITIVE VARIABLES LTT, DT, LRONG; ***************************************** * Khai bao cac ham * ***************************************** EQUATIONS * Ham muc tieu LOINHUAN Tong loi nhuan tu doi tau khai thac mot chuyen bien *Cac rang buoc ve dau vao chuyen bien CBIEN(k,l) So luong chuyen bien cua tau me lan khai thac cua tau (ngay) TOANDOI(k,l) Tong so luong tau toan doi TGCBM(k,l) Thoi gian chuyen bien cua tau me TGCBC(k,l) Thoi gian chuyen bien cua tau SANLUONGKT(k,l,S) San luong khai thac duoc chuyen bien DOANHTHU(k,l) CHIPHIBD(k,l) CHIPHICD LRCB Tong doanh thu chuyen bien Chi phi bien doi chuyen bien Chi phi co dinh chuyen bien Lai rong chuyen bien *Bieu thuc the hien dieu kien dau vao cua chuyen bien DKDAUVAO(k,l) Dieu kien dau vao cua chuyen bien *Bieu thuc the hien dieu kien dau vao cua chu tau LOINHUANCT(k,l) Tong loi nhuan cua chu tau BAODUONG(k,l) CLECHCT(k,l) Chi phi co dinh chuyen bien Chenh lech giua loi nhuan va chi phi co dinh dk1(k,l) jggjkh DK2(k,l) Dieu kien rang buoc ve tai tau TGTHUMUA(k,l) Khoang thoi gian giua hai lan thu mua cua tau me (Ngay) DKCT(k,l) Dieu kien dau vao cua chu tau *Bieu thuc the hien dieu kien dau vao cua thuy thu LUONGTT Muc luong thuc te cua thuy thu tren doi tau i LUONGMM(k,l) CLECH(k,l) Muc luong du kien cua thuy thu tren doi tai i Chenh lech giua muc luong duoc nhan va muc luong du kien cua thuy thu cua doi tau i DKTT(k,l) Dieu kien dau vao cua thuy thu; ***************************************** * Equation Definitions * ***************************************** * Ham muc tieu LOINHUAN LNCT =E= LRONG - CPCD - LTT; **Dieu kien rang buoc ve so luong chuyen bien va doi tau: *Chieu dai chuyen bien cua tau me TGCBM(k,l) CBM(k,l)=E=TG('Tme',k,l,'TGNB')+TG('Tme',k,l,'TGDC')+TG('Tme',k,l,'TGGT') *V2(k,l); *Chieu dai chuyen bien cua tau khai thac TGCBC(k,l) CBC(k,l) =E= TGKT(k,l)+TG('Tcon',k,l,'TGNB')+ TG('Tcon',k,l,'TGDC'); *Khoang thoi gian thu mua cua tau me (Ngay) TGTHUMUA(k,l) TGTM(k,l)=E= CBM(k,l)/V1(k,l); *So luong chuyen bien cua tau me dot khai thac cua tau CBIEN(k,l) SCBM(k,l)*TGTM(k,l)=E= TGKT(k,l); *So luong tau toan doi TOANDOI(k,l) V3(k,l) =E= V1(k,l)+V2(k,l); ** Dieu kien rang buoc ve tai *Dieu kien rang buoc ve tai tau me dk1(k,l) TTRM(k,l)=E= (Sum(s,CPUE(k,l,s)*TGTM(k,l))*V2(k,l))/1000; *Dieu kien rang buoc ve tai tau DK2(k,l) SUM(S,(CPUE(k,l,s)*TGTM(k,l)))/1000 =L= TTRC; *San luong khai thac duoc cua doi tau tren chuyen bien (tan) SANLUONGKT(K,L,S) SLKT(K,L,S)=E=(CPUE(K,L,S)-((V2(k,l)10)*a)*CPUE(K,L,S))*(TGKT(k,l) *V2(k,l))/1000; * Doanh thu cho chuyen bien (Trieu dong) DOANHTHU(K,L) DT(K,L) =E= SUM(S,(SLKT(K,L,S)*GIA(L,S))); * Chi phi bien doi mot chuyen bien (Trieu dong) CHIPHIBD(k,l) CPBD(k,l)=E=(CP('Tcon','CPDB')*TGKT(k,l)+CP('Tcon','CPDC')*TG('Tcon',k,l,'T GDC')+CP('Tcon','CPNB')*TG('Tcon',k,l,'TGNB'))*V2(k,l)+(CP('Tme','CPGT') *TG('Tme',k,l,'TGGT')*V2(k,l)+CP('Tme','CPDC')*TG('Tme',k,l,'TGDC')+ CP('Tme','CPNB')*TG('Tme',k,l,'TGNB'))*V1(k,l)*SCBM(k,l); *Chi phi co dinh mot chuyen bien CHIPHICD CPCD=E=sum((k,l),(CP('Tme','CPCD')/365*V1(k,l)*CBC(k,l)+CP('Tcon','CPCD')/ 365* V2(k,l)*CBC(k,l))); * Lai rong cua chuyen bien LRCB LRONG =E= SUM((k,l),(DT(K,L)-CPBD(k,l))); * Ham dieu kien dau vao cua chuyen bien DKDAUVAO(K,L) LRONG*V2(k,l) =G= 0; ***** Ham dieu kien dau vao cua chu tau: *Phan chia loi nhuan cua chu tau LOINHUANCT(k,l) DTCT(k,l) =E= TL*LRONG; *Chi phi bao duong tau hang nam chuyen bien BAODUONG(k,l) CPHN(k,l)=E=(CP('Tme','CPCD')/365)*CBC(K,L)*V1(k,l)+(CP('Tcon','CPCD')/36 5)*CBC(K,L)*V2(k,l); *Chenh lecnh giua loi nhuan va chi phi cua chu tau CLECHCT(k,l) CLCT(k,l) =E= DTCT(k,l)-CPHN(k,l); *Dieu kien dau vao cua chu tau DKCT(k,l) CLCT(k,l)*V2(k,l) =G= 0; ***** Ham dieu kien dau vao cua thuy thu: *Phan chia loi nhuan cho thuy thu LUONGTT LTT =E= (1-TL)*LRONG; *Muc luong mong muon cua thuy thu LUONGMM(k,l) LMM(k,l) =E= LCB*CBC(k,l)*(SLTT('Tme')+SLTT('Tcon')*V2(k,l)); *Chenh lech giua muc luong thuc te va luong mong muon CLECH(k,l) HS(k,l) =E= LTT-LMM(k,l); *Dieu kien tham gia cua thuy thu DKTT(k,l) HS(k,l)*V2(k,l) =G= 0; MODEL MECON /ALL/; V1.FX(k,l) = 1; V2.L(k,l) = 1; TGKT.l(k,l)=1; SLKT.L(K,L,S) = 1; LRONG.L = 1; CBM.Up(k,l)=12; CBC.UP(k,l)=60; DKCT.L(K,l) = 0; DKTT.L(k,l) = 0; SOLVE MECON USING NLP MAXIMIZING LNCT; V2.FX(k,l) = FLOOR(V2.L(k,l)); V2.UP(k,l) = FLOOR(V2.L(k,l)); SCBM.FX(k,l) = FLOOR(SCBM.L(k,l)); SCBM.UP(k,l) = FLOOR(SCBM.L(k,l)); SOLVE MECON USING NLP MAXIMIZING LNCT ; ********************************* * Result Export * ********************************* FILE RESULTB /HFMM2B.TXT/; PUT RESULTB; RESULTB.PC = 5; parameter misc(*); misc("LOINHUAN")=LNCT.l; misc("modelstat")=MECON.modelstat; Execute_Unload 'Mcon.gdx',V1,V2,V3,TTRM,SCBM,CBM,CBC,LNCT,DT,CPBD,CPCD,LTT,mis c; Execute 'GDXXRW Mcon.gdx o=Mcon.xls skipempty=0 zeroout=0 trace=2 index=OUT!B2'; Phần 3: CHƯƠNG TRÌNH TỐI ƯU HÀNH TRÌNH DI CHUYỂN CỦA TÀU MẸ TRONG SUỐT HÀNH TRÌNH THU MUA SẢN PHẨM KHAI THÁC TRÊN GAMS ********************************* * Hanh trinh di chuyen cua tau thu mua* ********************************* **Lua chon quang duong di chuyen cua tau me** set ii Cac vi tri den (cang ca va cac tau con) / i1*i8 / i(ii) subset of cities alias (ii,jj),(i,j,k); *Xac dinh so luong vi tri tau me muon di chuyen: i(ii) = ord(ii) 1 and ord(jj)>1; *Dieu kien rang buoc ve cac nhanh re: Nhanh(ij(i,j)) u(i) - u(j) + card(i)*x(i,j) =l= card(i)-1; model HTTM / muctieu, Hang, Cot, Nhanh /; option limcol=0,limrow=0; solve HTTM z using mip; display x.l; parameter misc(*); misc("Quangduong")=z.l; misc("modelstat")=HTTM.modelstat; Execute_Unload 'kcach.gdx',x,u,misc; Execute 'GDXXRW Kcach.gdx o=Kcach.xls trace=2 index=OUT!E1'; * skipempty=0 zeroout=0 Phụ Lục 3: MỘT SỐ HOẠT ĐỘNG KHẢO SÁT THỰC TẾ TẠI CẢNG CÁ TỈNH BÌNH THUẬN TRƯỜNG ĐẠI HỌC NHA TRANG CỘNG HOÀ - XÃ HỘI - CHỦ NGHĨA – VIỆT NAM Độc Lập – Tự Do – Hạnh Phúc Khoa Kỹ Thuật Giao Thông Nha Trang Ngày 10 Tháng Năm 2016 Báo Cáo Tổng Hợp (Kết khảo sát số liệu thực đề tài) NỘI DUNG - Khảo sát hoạt động của đội tàu khai thác thuỷ sản của tỉnh Bình Thuận VỊ TRÍ - Chi cục thuỷ sản Bình Thuận - Cảng cá Phan Rí Cửa – Huyện Tuy Phong – Bình Thuận - Cảng cá Phan Thiết – TP Phan Thiết – Bình Thuận THÀNH PHẦN THAM GIA STT Họ tên MSSV Nguyễn Thanh Tín 56CH212 Nguyễn Văn Quang 54131881 Hà Quốc Việt 54131786 Nguyễn Ngọc Quốc Việt 54131758 Nhiệm vụ Phỏng vấn chủ tàu, điền phiếu khảo sát; tổng hợp số liệu Thu thập nhật kí đánh bắt, tổng hợp số liệu Phỏng vấn chủ tàu, điền phiếu khảo sát; tổng hợp số liệu Thu thập nhật kí đánh bắt, tổng hợp số liệu KẾT QUẢ THỰC HIỆN - Thu thập 925 mẫu báo cáo nhật kí đánh bắt năm 2015 - Thu thập 214 mẫu báo cáo nhật kí đánh bắt năm 2016 - Phỏng vấn điền phiếu điều tra 150 phiếu (trong tại cảng cá Phan Rí Cửa 65 phiếu, cảng cá Phan Thiết 85 phiếu) ... THANH TÍN NGHIÊN CỨU XÁC ĐỊNH QUY MÔ ĐỘI TÀU KHAI THÁC HẢI SẢN PHÙ HỢP VỚI NGHỀ CÁ CỦA VIỆT NAM LUẬN VĂN THẠC SĨ KỸ THUẬT Ngành: Kỹ thuật khí động lực Mã số: 60520116 Quy t định giao... đề tài: Nghiên cứu xác định quy mô đội tàu khai thác hải sản phù hợp với nghề cá Việt Nam Số liệu luận văn hoàn toàn trung thực, xử lý theo phương pháp khoa học đảm bảo độ tin cậy Mô t số... Thủy Sản Bình Thuận, KS Nguyễn Văn Cảnh phòng đăng kiểm cung cấp tạo điều kiện để sử dụng số liệu thực đề tài: Nghiên cứu xác định quy mô đội tàu khai thác hải sản phù hợp với nghề cá Việt