Sử dụng nguồn số liệu vệ tinh trích rút mực nước và độ cao sóng có nghĩa trên biển đông và kiểm định mô hình toán

ỦY BAN NHÂN DÂN THÀNH PHỐ HỒ CHÍ MINH SỞ KHOA HỌC VÀ CÔNG NGHỆ SỞ KHOA HỌC VÀ CÔNG NGHỆ TRUNG TÂM PHÁT TRIỂN KHOA HỌC VÀ CÔNG NGHỆ TRẺ CHƯƠNG TRÌNH KHOA HỌC VÀ CƠNG NGHỆ CẤP THÀNH PHỐ BÁO CÁO TỔNG HỢP KẾT QUẢ NHIỆM VỤ NGHIÊN CỨU KHOA HỌC VÀ CÔNG NGHỆ SỬ DỤNG NGUỒN SỐ LIỆU VỆ TINH TRÍCH RÚT MỰC NƯỚC VÀ ĐỘ CAO SĨNG CĨ NGHĨA TRÊN BIỂN ĐƠNG VÀ KIỂM ĐỊNH MƠ HÌNH TOÁN Cơ quan chủ trì nhiệm vụ: Trung tâm phát triển Khoa học Công nghệ trẻ Chủ nhiệm nhiệm vụ: ThS Ngô Nam Thịnh Thành phố Hồ Chí Minh - 20… ỦY BAN NHÂN DÂN THÀNH ĐỒN TP.HCM THÀNH PHỐ HỒ CHÍ MINH TRUNG TÂM PHÁT TRIỂN KHOA HỌC SỞ KHOA HỌC VÀ CÔNG NGHỆ VÀ CÔNG NGHỆ TRẺ CHƯƠNG TRÌNH KHOA HỌC VÀ CƠNG NGHỆ CẤP THÀNH PHỐ BÁO CÁO TỔNG HỢP KẾT QUẢ NHIỆM VỤ NGHIÊN CỨU KHOA HỌC VÀ CÔNG NGHỆ SỬ DỤNG NGUỒN SỐ LIỆU VỆ TINH TRÍCH RÚT MỰC NƯỚC VÀ ĐỘ CAO SĨNG CĨ NGHĨA TRÊN BIỂN ĐƠNG VÀ KIỂM ĐỊNH MƠ HÌNH TOÁN Chủ nhiệm nhiệm vụ: Ngơ Nam Thịnh Cơ quan chủ trì nhiệm vụ Đồn Thành Kim Thành phố Hồ Chí Minh- 2017 TRUNG TÂM PHÁT TRIỂN KHOA HỌC VÀ CÔNG NGHỆ TRẺ CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc ., ngày tháng năm 20 BÁO CÁO THỐNG KÊ KẾT QUẢ THỰC HIỆN NHIỆM VỤ NGHIÊN CỨU KH&CN I THÔNG TIN CHUNG Tên nhiệm vụ: Sử dụng nguồn số liệu vệ tinh trích rút mực nước độ cao sóng có nghĩa biển Đơng kiểm định mơ hình tốn Thuộc: Chương trình/lĩnh vực (tên chương trình/lĩnh vực): Khoa học tự nhiên Chủ nhiệm nhiệm vụ: Họ tên: Ngô Nam Thịnh Ngày, tháng, năm sinh: 14/04/1985 Nam/ Nữ: Nam Học hàm, học vị: Thạc sỹ Chức danh khoa học: Giảng viên Điện thoại: Tổ chức: 02838443006 Chức vụ Nhà riêng: Mobile: 0985101261 Fax: 02838449474 E-mail: namthinh2007@gmail.com Tên tổ chức công tác: Trường Đại học Tài nguyên Môi trường TP.HCM Địa tổ chức: 236B Lê Văn Sỹ, phường 1, Quận Tân Bình, TP.HCM Địa nhà riêng: I408, CC Đào Duy Từ, Hẻm 51 Thành Thái, Phường 14, Q.10, TP.HCM 3 Tổ chức chủ trì nhiệm vụ: Tên tổ chức chủ trì nhiệm vụ: Trung tâm Phát triển Khoa học Công nghệ trẻ Điện thoại: 38.233.363 Fax: 38244705 Website: www.khoahoctre.com.vn Địa chỉ: Số Phạm Ngọc thạch, P Bến Nghé, Quận Họ tên thủ trưởng tổ chức: ĐOÀN KIM THÀNH Số tài khoản: 3713.0.1083277.00000 Kho bạc: Nhà nước Quận Thành phố Hồ Chí Minh Tên quan chủ quản đề tài: Sở Khoa học Cơng nghệ TP.HCM II TÌNH HÌNH THỰC HIỆN Thời gian thực nhiệm vụ: - Theo Hợp đồng ký kết: từ tháng 10/2016 đến tháng 10/2017 - Thực tế thực hiện: từ tháng 10/2016 đến tháng 10/2017 - Được gia hạn (nếu có): Kinh phí sử dụng kinh phí: a) Tổng số kinh phí thực hiện: 80 tr.đ, đó: + Kính phí hỗ trợ từ ngân sách khoa học: 80 tr.đ + Kinh phí từ nguồn khác: ……………….tr.đ b) Tình hình cấp sử dụng kinh phí từ nguồn ngân sách khoa học: Số TT Theo kế hoạch Thực tế đạt Ghi Thời gian Kinh phí Thời gian Kinh phí (Số đề nghị (Tháng, (Tr.đ) (Tháng, năm) (Tr.đ) toán) năm) 10/2016 40 10/2016 40 10/2017 40 10/2017 40 … c) Kết quả sử dụng kinh phí theo khoản chi: Đối với đề tài: Đơn vị tính: Triệu đồng Số Nội dung TT khoản chi Theo kế hoạch Tổng NSKH Nguồn Thực tế đạt Tổng NSKH Nguồn khác Trả công lao động 69,722 69,722 khác 69,722 69,722 10,278 10,278 (khoa học, phổ thông) Nguyên, vật liệu, lượng Thiết bị, máy móc Xây dựng, sửa chữa nhỏ Chi khác Tổng cộng 10,278 80 10,278 80 - Lý thay đổi (nếu có): 80 80 Đối với dự án: Đơn vị tính: Triệu đồng Số Nội dung TT khoản chi Thiết bị, máy móc Theo kế hoạch Tổng NSKH Nguồn Thực tế đạt Tổng NSKH Nguồn khác khác mua Nhà xưởng xây dựng mới, cải tạo Kinh phí hỗ trợ cơng nghệ Chi phí lao động Nguyên vật liệu, lượng Thuê thiết bị, nhà xưởng Khác Tổng cộng - Lý thay đổi (nếu có): Các văn hành q trình thực đề tài/dự án: (Liệt kê định, văn quan quản lý từ công đoạn xét duyệt, phê duyệt kinh phí, hợp đồng, điều chỉnh (thời gian, nội dung, kinh phí thực có); văn tổ chức chủ trì nhiệm vụ (đơn, kiến nghị điều chỉnh có) Số Số, thời gian ban TT hành văn Tên văn Ghi … Tổ chức phối hợp thực nhiệm vụ: Số TT Tên tổ chức Tên tổ chức Nội dung Sản phẩm đăng ký theo tham gia thực tham gia chủ yếu Thuyết minh chủ yếu đạt Ghi chú* - Lý thay đổi (nếu có): Cá nhân tham gia thực nhiệm vụ: (Người tham gia thực đề tài thuộc tổ chức chủ trì quan phối hợp, khơng 10 người kể chủ nhiệm) Số Tên cá nhân Tên cá nhân Nội dung Sản phẩm T đăng ký theo tham gia tham gia chủ yếu đạt T Thuyết minh thực ThS Ngơ Nam Thịnh ThS Trần Thị Kim Nội dung 1,3,4,5,6,7,8, ThS Trần Thị Nội dung Kim ThS Võ Duy Long ThS Võ Duy Long ThS Ngô Nam Thịnh Nội dung Bộ số liệu báo cáo đáng tin cậy Báo cáo đánh giá thống kê Báo cáo đánh giá tin cậy Ghi chú* ThS Nguyễn Iêng Vũ ThS Nguyễn Trâm Anh Nội dung ThS Thái Tiểu Minh CN Phùng Thị Mỹ Diễm Nội dung 4,8 CN Nguyễn Thị Nhật Anh Nội dung CN Trương Thiên Phúc Nội dung CN Nguyễn Thị Thùy Trang CN Phan Thùy Linh Nội dung Nội dung Báo cáo đánh giá tin cậy Báo cáo đánh giá tin cậy Báo cáo đánh giá tin cậy Báo cáo đánh giá tin cậy Báo cáo đánh giá tin cậy Báo cáo đánh giá tin cậy - Lý thay đổi ( có): Một số cá nhân bận cơng việc Tình hình hợp tác quốc tế: Số TT Theo kế hoạch Thực tế đạt (Nội dung, thời gian, kinh (Nội dung, thời gian, kinh phí, địa điểm, tên tổ chức phí, địa điểm, tên tổ chức hợp tác, số đồn, số lượng hợp tác, số đoàn, số lượng người tham gia ) người tham gia ) - Lý thay đổi (nếu có): Tình hình tổ chức hội thảo, hội nghị: Ghi chú* Số TT Theo kế hoạch Thực tế đạt (Nội dung, thời gian, kinh (Nội dung, thời gian, phí, địa điểm ) Ghi chú* kinh phí, địa điểm ) - Lý thay đổi (nếu có): Tóm tắt nội dung, công việc chủ yếu: (Nêu mục 15 thuyết minh, không bao gồm: Hội thảo khoa học, điều tra khảo sát nước nước ngoài) Thời gian Số TT Các nội dung, công việc (Bắt đầu, kết thúc Người, chủ yếu - tháng … năm) quan (Các mốc đánh giá chủ yếu) Theo kế Thực tế hoạch đạt thực Thu thập liệu địa hình, khí 10/2016 – 10/2016 – Trần Thị tượng, hải văn, số liệu vệ tinh 12/2016 12/2016 Kim (sóng, mực nước) năm 2009, 2010 Đánh giá tổng quan loại 10/2016 – 10/2016 – Nguyễn Thị vệ tinh theo không thời 11/2016 11/2016 Nhật Anh Tải xử lý số liệu vệ tinh 10/2016 – 10/2016 – Trương phục vụ hiệu chỉnh kiểm 12/2016 12/2016 Thiên Phúc gian đo đạc tại biển Đông sóng, mực nước định với số liệu thực đo Xây dựng lưới tính địa 10/2016 – hình, loại số liệu biên khu 11/2016 10/2016 – Phùng Thị 11/2016 Mỹ Diễm vực nghiên cứu phục vụ mơ hình tốn Hiệu chỉnh kiểm định mơ 10/2016 – 10/2016 – Ngũn Thị hình tính tốn sóng Mike 21 12/2016 12/2016 Thùy Trang Ứng dụng mơ hình Mike 21 10/2016 – 10/2016 – Võ Duy SW (Spectral Waves) tính 12/2016 12/2016 Long Hiệu chỉnh kiểm định mơ 1/2017 – 1/2017 – Ngũn hình tính tốn dao động mực 5/2017 5/2017 Trâm Anh Ứng dụng mơ hình Mike 21 6/2017 – 6/2017 – Phùng Thị Flow Model FM tính tốn 8/2017 8/2017 Mỹ Diễm Đánh giá kết quả tính tốn 7/2017 – 7/2017 – Phan Thùy mơ hình tốn với số liệu 9/2017 9/2017 Linh SW (Spectral Waves) với số liệu thực đo tốn sóng tháng đặc trưng mùa gió (Đơng Bắc, Tây Nam) so sánh với số liệu vệ tinh nước Mike 21 Flow Model FM với số liệu thực đo mực nước tháng đặc trưng mùa gió (Đơng Bắc, Tây Nam) so sánh với số liệu vệ tinh 10 Hình 3.19 Kết tính tốn dao động mực nước khu vực Biển Đơng lúc 21 ngày 4/7/2009 Hình 3.20 Kết tính tốn dao động mực nước khu vực Biển Đông lúc 20 ngày 8/12/2009 70 3.2.4 Bộ thơng số mơ hình Kết quả tính tốn mơ hình sóng triều cho kết quả phù hợp với số liệu thực đo Vì thế, thơng số mơ hình lựa chọn sau: Bộ thông số Mike 21 SW Các tham số khác Số hướng Số tần số 16 25 Alpha = Gamma = 0.8 0.04 m Từ NOAA Tham số sóng vỡ Ma sát đáy Gió Điều kiện ban đầu Sigma a 0.07 Sigma b 0.09 Tham số đỉnh 3.3 Đà gió 100km Tần số đỉnh 0.4 Bộ thông số Mike 21 HD FM Từ mơ hình tồn cầu Mike 21 Biên mực nước Toolbox Manning 32 eddy 0.28 m2/s Gió Từ NOAA Mực nước ban đầu 0.5 m 3.3 So sánh giữa mơ hình tốn với số liệu vệ tinh 3.3.1 Kết so sánh giữa mơ hình sóng Mike 21 SW với số liệu vệ tinh Mơ hình sóng Mike 21 SW ứng dụng tính tốn độ cao sóng có nghĩa so sánh với số liệu vệ tinh tại trạm DK1 Để có nhiều liệu so sánh, mơ hình mơ sóng tháng gió mùa Tây Nam (tháng 6-8) tháng gió mùa Đơng Bắc (tháng 10 – 12) Kết quả so sánh thể sau: 71 Hình 21 Tương quan số liệu vệ tinh mơ hình độ cao sóng có nghĩa năm 2009 Hình 3.22 Kết so sánh số liệu vệ tinh mơ hình độ cao sóng có nghĩa tháng năm 2009 Số liệu vệ tinh sử dụng để so sánh số liệu lấy khoảng không gian với bán kính 0.50 thời gian khoảng Kết quả so sánh độ cao sóng có nghĩa số liệu vệ tinh mơ hình tại vị trí DK1 cho thấy hai chuỗi số liệu tương quan mức khá, với hệ số tương quan R2= 0.62 Sai số trung bình hai chuỗi số liệu khoảng 11% (Hình 3.22) Với hệ số tương quan sai số hai chuỗi số liệu cho thấy số liệu vệ tinh có khả sử dụng 72 việc hiệu chỉnh kiểm định mơ hình tốn, số liệu có thể thay số liệu thực đo biển sâu 3.3.2 Kết so sánh giữa mơ hình triều Mike 21 HD FM với số liệu vệ tinh Số liệu mực nước trung bình mà vệ tinh đo tại trạm hải văn Vũng Tàu cũng sử dụng để so sánh với kết quả tính tốn mơ hình triều Số liệu vệ tinh xử lý tại trạm Vũng Tàu cũng tương tự số liệu vệ tinh mục 3.2.2 Mơ hình tính tốn triều Biển Đơng cũng tính tháng gió mùa Tây Nam (tháng – 8) tháng gió mùa Đơng Bắc (tháng 10 – 12) Kết quả so sánh mực nước trung bình thể hình 3.23 3.24 Hình 3.23 Tương quan số liệu vệ tinh mơ hình độ cao mực nước trung bình năm 2009 73 Hình 3.24 Kết so sánh số liệu vệ tinh mô hình độ cao mực nước trung bình năm 2009 Kết quả so sánh mực nước trung bình số liệu vệ tinh mơ hình cho thấy chuỗi số liệu có tương quan tốt, hệ số tương quan khoảng 0.97 Tuy nhiên, xét độ cao có sai lệch đáng kể (Hình 3.24) Đây sai khác hệ cao độ hai chuỗi số liệu, điều cũng có thể khắc phục cách sử dụng hàm tương quan (hình 3.19) để đưa số liệu vệ tinh gần sát với số liệu thực đo trạm Vũng Tàu để so sánh thời điểm khác KẾT LUẬN VÀ KIẾN NGHỊ Các vệ tinh đo cao Jason Jason sử dụng để trích rút số liệu độ cao sóng có nghĩa độ cao mực nước khu vực Biển Đơng năm 2009 2010 Chương trình tải liệu vệ tinh đo cao độ cao sóng có nghĩa độ cao mực nước khu vực Biển Đông xử lý số liệu xây dựng cách tự động, qua có thể tiếp tục sử dụng làm tiền đề cho nghiên cứu sử dụng loại số liệu khác vệ tinh Đề tài ứng dụng mơ hình Mike 21 để tính tốn mơ sóng triều khu vực Biển Đông Đồng thời hiệu chỉnh kiểm định mơ hình với số liệu thực đo để tìm thơng số cho mơ hình sóng triều khu vực Biển Đơng Kết quả nghiên cứu cho thấy, số liệu vệ tinh độ cao sóng có nghĩa độ cao mực nước tương quan mức so với số liệu thực đo tại trạm giàn khoang DK1-7, Vũng Tàu trạm đo đạc Cà Mau Bán kính thu thập số liệu vệ 74 tinh tính từ tâm vị trí số liệu đo đạc lựa chọn từ 0.50 đến 1.50 độ cao sóng 0.50 mực nước Qua đó, trường hợp khơng có số liệu thực đo ngồi biển khơi độ cao sóng, độ cao mực nước có thể sử dụng hàm tương quan để hiệu chỉnh số liệu vệ tinh số liệu thực đo phục vụ hiệu chỉnh kiểm định mơ hình tốn cũng đánh giá khác nghiên cứu biển Công nghệ vệ tinh phát triển vượt bậc, vệ tinh dành cho nghiên cứu biển cũng quan tâm đầu tư cả nghiên cứu khoa học lẫn an ninh quốc phòng quốc gia Số liệu vệ tinh đáp ứng cập nhật số liệu, liên tục theo khơng gian thời gian Vì thế, việc sử dụng số liệu vệ tinh nghiên cứu nói chung nghiên cứu biển nói riêng mang lại hiệu quả định Hạn chế Hạn chế đề tài so sánh với số liệu thực đo độ cao sóng có nghĩa, độ cao mực nước trạm phía Nam Biển Đơng, cần bổ sung so sánh đánh giá nhiều trạm khác có số liệu thực đo để có kết quả đáng tin cậy Vệ tinh có quỹ đạo vận tốc định, khơng thể bao trùm hết Biển Đông cùng thời điểm mà vào thời gian khác nhau, số liệu vệ tinh cũng khơng liên tục tại vị trí cố định thường có sai lệch vị trí đo đạc với số liệu vệ tinh nên kết quả so sánh chưa thật mong đợi Độ cao bề mặt biển từ số liệu vệ tinh chưa chuyển cao độ chuẩn Việt Nam nên có chênh lệch đáng kể giá trị độ lớn, đúng xu hướng dâng rút Kiến nghị 75 Tiếp tục nghiên cứu bổ sung đánh giá nhiều trạm đo đạc với số liệu thực đo khu vực phía Bắc nhằm đảm bảo độ tin cậy số liệu vệ tinh Nghiên cứu sử dụng số liệu vệ tinh yếu tố khác như: gió, nhiệt độ,… TÀI LIỆU THAM KHẢO Tiếng Anh Ardhuin, F et al., 2007 Comparison of wind and wave measurements and models in the Western Mediterranean Sea Ocean Engineering, 34(3-4), pp.526–541 Carter D J T, 1993 Estimating extreme wave heights in the NE Atlantic from Geosat data Technical Report OTH 93 396, Health and Safety Executive Offshore Technology Report, HMSO, London Dimitrova, M., Kortcheva, A & Galabov, V., 2013 Validation of the operational wave model WAVEWATCH III against altimetry data from JASON-2 satellite Bulgarian Journal of Meteorology and Hydrology, 1(2), pp.4–17 George Galanis, Dan Hayes, George Zodiatis, Peter C Chu, Yu-Heng Kuo, George Kallos, 2011 Wave height characteristics in the Mediterranean Sea by means of numerical modeling, satellite data, statistical and geometrical techniques Marine Geophysical Research Queffeulou, P., 2003 Long-term Quality Status of Wave Height and Wind Speed Measurements from Satellite Altimeters In The Thirteenth (2003) International Offshore and Polar Engineering Conference pp 129 – 135 Tiếng Việt Nguyễn Duy Khang nnk, 2015, Khảo sát xu thay đổi mực nước biển vùng Nam Bộ sử dụng số liệu đo cao vệ tinh, hội thảo chương trình khoa học cơng nghệ vũ trụ giai đoạn 2012 – 2015 Nguyễn Kỳ Phùng, 2010, “Nghiên cứu trình tương tác biển - lục địa ảnh hưởng chúng đến hệ sinh thái ven bờ Đông bờ Tây Nam Bộ”, KC.09/06-10 76 Internet Trần Tuấn Dũng, Đo cao vệ tinh nghiên cứu cấu trúc địa chất biển Đông Việt Nam, link: http://www.vast.ac.vn/khoa-hoc-va-phat-trien/nghien-cuu-co-ban/979-docao-ve-tinh-trong-nghien-cuu-cau-truc-dia-chat-bien-dong-viet-nam (ngày truy cập 27/06/2017) Nguồn Wikipedia: https://vi.wikipedia.org/wiki/Sputnik_1 (ngày truy cập 24/03/2017) 10 Nguồn Wikipedia: https://vi.wikipedia.org/wiki/V%E1%BB%87_tinh (ngày truy cập 24/03/2017) 11 Nguồn http://www.aviso.altimetry.fr/: (ngày truy cập 14/07/2017) 12 http://www.aviso.altimetry.fr/en/techniques/altimetry/principle/pulses-andwaveforms.html(ngày truy cập 14/07/2017) 13 Nguồn http://www.aviso.altimetry.fr/: http://www.aviso.altimetry.fr/en/techniques/altimetry/principle/frequenciesused.html(ngày truy cập 15/07/2017) 14 Nguồn http://www.altimetry.info/: (ngày truy cập 15/07/2017) 15 http://www.altimetry.info/missions/past-missions/topexposeidon/tp-instruments/topexaltimeter-or-nra/ (ngày truy cập 16/07/2017) 16 Nguồn http://www.altimetry.info/: (ngày truy cập 16/07/2017) 17 http://www.altimetry.info/missions/past-missions/topexposeidon/tpinstruments/poseidon-1/ (ngày truy cập 16/07/2017) 18 Nguồn http://www.altimetry.info/: (ngày truy cập 18/07/2017) 19 http://www.altimetry.info/missions/past-missions/topexposeidon/tp-instruments/tmrtopex-microwave-radiometer/ (ngày truy cập 18/07/2017) 20 Nguồn http://www.altimetry.info/:http://www.altimetry.info/missions/pastmissions/topexposeidon/tp-objectives/ (ngày truy cập 19/07/2017) 21 Nguồn http://www.altimetry.info/missions/past-missions/jason-1/ (ngày truy cập 19/07/2017) 77 PHỤ LỤC Phụ lục Chương trình tải dữ liệu vệ tinh clear all; clc; ftp_client = ftp('avisoftp.cnes.fr'); cd(ftp_client,'AVISO/pub/jason-1/gdr_e/') for i=257:258 % Tu 257-331 formatSpec = 'cycle_%d'; str = sprintf(formatSpec,i)% Thiet lap string cycle_i cd(ftp_client,str) m = dir(ftp_client,'*.nc') filename = {m.name}; for k = 3:length(filename)% k chay tu vi tai 1, filename la NaN name = filename{k}(17:19) if isequal(name,'012') || isequal(name,'025')|| isequal(name,'038')|| isequal(name,'051') || isequal(name,'101')|| isequal(name,'114')|| isequal(name,'127')|| isequal(name,'190') || isequal(name,'203')|| isequal(name,'242'); % Cac file nam vung nghien cuu: 12, 25,38, 51,64,77,114,127,101,140,153,190,203,216,229,242 disp(' Downloading ') filename1 = filename{k} mget(ftp_client,filename1,'D:\congviec\2015\vuon uom\Download_vetinh') end end cd (ftp_client,' ')% return directory end close(ftp_client); 78 Phụ lục Chương trình xử lý tổng hợp số liệu vệ tinh Biển Đông clear all; clc; cd 'D:\cong viec\2015\vuon uom\Data daxuly\Jason2' fid = fopen('data1.txt','a'); fprintf(fid,'%s %s %s %s %s %s \n','Date','Time','Lon','Lat','Swh','Ssh') formatSpec = '%s %f %f %f %f\n'; m = dir; z={m.name}; for i=520:(length(z)-3) S = netcdf(z{i}) for j=1:length(S.VarArray(1,1).Data) Time(j)=((S.VarArray(1,1).Data(j))/86400)+730486; % Time tinh tu 01/01/2000 co datenum bang 730486 Vi vay Time = (S.VarArray(1,1).Data(j))/86400+730486 Date = datestr(Time(j),'dd/mm/yyyy HH:MM:SS');% Chuyen dinh dang duoi dang so sang dang if((S.VarArray(1,5).Data(j)*(1.0000e-006)101)&& S.VarArray(1,4).Data(j)*(1.0000e-006)4) fprintf(fid,formatSpec,Date,double(S.VarArray(1,5).Data(j))*(1.0000e006),double(S.VarArray(1,4).Data(j))*(1.0000e-006), double(S.VarArray(1,73).Data(j))*0.001,double(S.VarArray(1,147).Data(j))*0.0 01) end end end fclose(fid); 79 Phụ lục Chương trình xử lý xác định bán kính từ vị trí xác định trước clear all; clc; cd 'D:\cong viec\2015\vuon uom\Data daxuly\Jason2'; fid = fopen('data1.txt','r');% So lieu ve tinh -n = numel(textread('data1.txt','%1c%*[^\n]'))%Dem so hang file data1.txt cd % Xac dinh xem tinh SWH hay SSH -prompt1 = 'Co tinh SWH hay khong?(yes:1) '; askswh = input(prompt1); prompt2 = 'Co tinh SSH hay khong?(yes:1) '; askssh = input(prompt2); % -Diem can tinh la: -c=[110.62055 8.01773]; % Toa DK1 if (askswh==1) fid1 = fopen('Write_Swh_DK1.txt','r');% So lieu dac SWH -fid2 = fopen('WriteSWH.txt','w'); n1 = numel(textread('Write_Swh_DK1.txt','%1c%*[^\n]'))%Dem so hang file Write_Swh_DK1.txt end if (askssh==1) fid4 = fopen('Write_Ssh_DK1.txt','r');% So lieu dac SSH fid3 = fopen('WriteSSH.txt','w') end n2 = numel(textread('Write_Ssh_DK1.txt','%1c%*[^\n]'))%Dem so hang file Write_Ssh_DK1.txt name = textscan(fid,'%s %s %s %s %s %s \n',1); s = 1; v = 1; ns(1) = 733770.064756944; %00:00:00 01/01/2009 nv(1) = 733770.064756944; %00:00:00 01/01/2009 m = dir; for i = 1:(n-1) % -Merge cac diem tren track ve tinh vao tung nhom formatSpec = '%s %s %f %f %f %f\n'; [A1{i}] = textscan(fid,formatSpec,1); Date(i) = A1{i}(1); 80 Time(i) = A1{i}(2); Lon(i) = A1{i}(3); Lat(i) = A1{i}(4); dt1(i) = datenum((Date{i}), 'dd/mm/yyyy');% Chuyen dinh dang sang dinh dang so dt(i) = datenum((Time{i}), 'HH:MM:SS');% Chuyen dinh dang Gio sang dinh dang so ns(i) = datenum(strcat(Date{i},{' '},Time{i}),'dd/mm/yyyy HH:MM:SS');% ssh -nv(i) = datenum(strcat(Date{i},{' '},Time{i}),'dd/mm/yyyy HH:MM:SS');% swh -if (ns(i)-ns(1)>0.000162035) % -Khoang cach 15km-(16s) ssh -ns(1) = ns(i); s = s+1; end ss(i) = s; if (nv(i)-nv(1)>0.001)% -Khoang cach 100km -swh -nv(1) = nv(i); v = v+1; end vv(i) = v; for j=3:6 B(i,j) = A1{i}(j); A2(i,j) = cell2mat(B(i,j)); end % Hien thi phan tram tgian tinh toan per = i/(n-1)*100; disp(per);disp('%') end save('data.mat') ns(1) = 733770.064756944; nv(1) = 733770.064756944; [u ID] = unique(vv); % Loc cac diem 90s (100km) [u1 ID1]= unique(ss); % Loc cac diem vong 16 giay(15km) DK = [A2(:,3) A2(:,4)]% Lon Lat ve tinh % Nhap toa cac diem vao day-% - DK1(110.62055 8.01773); VungTau(107.16339 10.35477);PhuQuy(108.93333 %10.51666);BachLongVi(107.72808 20.13214) 81 % -Merge SWH for k=1:length(u) if (k == 1); swhM(k) = mean(A2((1:ID(k)),5)); % Tinh trung binh swh TMw(k) = nv(fix(mean(1:ID(k)))); LonM(k) = (A2((fix(mean(1:ID(k)))),3)) LatM(k) = (A2((fix(mean(1:ID(k)))),4)) else swhM(k) = mean(A2((ID(k-1)+1:ID(k)),5)); % % Tinh trung binh swh TMw(k) = nv(fix(mean(ID(k-1)+1:ID(k)))); LonM(k) = ( A2((fix(mean(ID(k-1)+1:ID(k)))),3)) LatM(k) = (A2((fix(mean(ID(k-1)+1:ID(k)))),4)) end end Lo = transpose(LonM);La = transpose(LatM) % -Merge SSH -for k=1:length(u1) if (k == 1); sshM(k) = mean(A2((1:ID1(k)),6)); % Tinh trung binh ssh TMs(k) = ns(fix(mean(1:ID1(k)))); LonN(k) = (A2((fix(mean(1:ID1(k)))),3)) LatN(k) = (A2((fix(mean(1:ID1(k)))),4)) else sshM(k) = mean(A2((ID1(k-1)+1:ID1(k)),6)); % % Tinh trung binh ssh TMs(k) = ns(fix(mean(ID1(k-1)+1:ID1(k)))); LonN(k) = ( A2((fix(mean(ID1(k-1)+1:ID1(k)))),3)) LatN(k) = (A2((fix(mean(ID1(k-1)+1:ID1(k)))),4)) end end Los = transpose(LonN);Las = transpose(LatN) % Nhap so lieu dac-SWH for jj = 1:n1 formatSpec1 = '%s %s %f\n'; [OBs{jj}] = textscan(fid1,formatSpec1,1); nOB(jj) = datenum(strcat(OBs{jj}{1},{' '},OBs{jj}{2}),'dd/mm/yyyy HH:MM:SS'); end % Nhap so lieu dac SSH for jj = 1:n2 formatSpec1 = '%s %s %f\n'; [OBss{jj}] = textscan(fid4,formatSpec1,1); 82 nOBs(jj) = datenum(strcat(OBss{jj}{1},{' '},OBss{jj}{2}),'dd/mm/yyyy HH:MM:SS'); end % -Tinh SWH (Tim diem lan can tinh SWH) if (askswh == 1) for zz = 1:n1 % -SWH c1= nOB(zz) % Date num cua Time OBs r1=0.08915; %0.08915 = tieng TTMw = transpose(TMw) % Chuyen vi TMw [iOB Dob]=rangesearch(c1,r1,TTMw)% Lay cac so lieu vong tieng % So sanh voi so lieu tai DK1 if isempty(iOB) disp(' -Khong tim thay iOB -') else r=1.5; % -1degree = 110,57km [idx D]=rangesearch(c,r,[Lo(iOB(:)) La(iOB(:))]) if isempty(idx) disp(' -Khong tim thay idx -') else %Tinh SWH swh(zz) = mean(swhM(iOB(idx(:)))); % Tinh trung binh swh TimeW(zz) = TTMw(fix(mean(iOB(idx(:))))); swh_obs(zz) = (OBs{zz}{3}); TimeOBsw(zz) = nOB(zz) fprintf(fid2,'%s %s %f %f \n',datestr(TimeW(zz),'dd/mm/yyyy HH:MM:SS'),datestr(TimeOBsw(zz),'dd/mm/yyyy HH:MM:SS'),swh(zz),swh_obs(zz))% end end end end % Tinh SSH if (askssh == 1) rs=0.15; %-(16.5km) 1degree = 110,57km [idxs Ds]=rangesearch(c,rs,[Los Las]) if isempty(idxs) 83 disp(' -Khong tim thay idxs -') else %Tinh SSH ssh = (sshM((idxs(:)))); % Lay nhung ssh ban kinh 16.5km TimeS = TMs(idxs(:)) fprintf(fid3,'%s %f \n',datestr(TimeS,'dd/mm/yyyy HH:MM:SS'),ssh)% end end fclose(fid); fclose(fid1); fclose(fid2); 84