Nghiên cứu lan truyền dầu ở biển đông phục vụ cảnh báo và tìm kiếm nguồn thải

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Quốc Trinh NGHIÊN CỨU LAN TRUYỀN DẦU Ở BIỂN ĐÔNG PHỤC VỤ CẢNH BÁO VÀ TÌM KIẾM NGUỒN THẢI LUẬN ÁN TIẾN SĨ HẢI DƯƠNG HỌC Hà Nợi - 2018 ĐẠI HỌC Q́C GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Quốc Trinh NGHIÊN CỨU LAN TRUYỀN DẦU Ở BIỂN ĐƠNG PHỤC VỤ CẢNH BÁO VÀ TÌM KIẾM NGUỒN THẢI Chuyên ngành: Hải Dương học Mã số: 62 44 02 28 LUẬN ÁN TIẾN SĨ HẢI DƯƠNG HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC HDC: PGS.TS Nguyễn Minh Huấn HDP: PGS.TS Phùng Đăng Hiếu XÁC NHẬN NCS ĐÃ CHỈNH SỬA THEO QUYẾT NGHỊ CỦA HỘI ĐỒNG ĐÁNH GIÁ LUẬN ÁN TM Tập thể hướng dẫn khoa học Chủ tịch Hội đồng đánh giá Luận án Tiến sĩ PGS TS Nguyễn Minh Huấn GS TS Đinh Văn Ưu Hà Nội - 2018 LỜI CAM ĐOAN Tôi xin cam đoan, cơng trình nghiên cứu riêng tơi Các số liệu tài liệu sử dụng luận án trung thực Các kết nghiên cứu nêu luận án trung thực chưa công bố cơng trình nghiên cứu khác Hà Nội, ngày tháng 03 năm 2018 Tác giả Nguyễn Quốc Trinh i LỜI CẢM ƠN Luận án hồn thành khoa Khí tượng, Thuỷ văn Hải dương học thuộc Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội giai đoạn từ năm 2011 đến 2016 Để luận án hoàn thành đạt kết này, nghiên cứu sinh gửi lời cảm ơn chân thành đến thầy hướng dẫn trực tiếp PGS TS Nguyễn Minh Huấn PGS TS Phùng Đăng Hiếu, thầy khoa Khí tượng, Thuỷ văn Hải dương học nói riêng Trường Đại học Khoa học Tự nhiên nói chung, bạn bè đồng nghiệp tạo điều kiện tốt vật chất tinh thần cho nghiên cứu sinh học tập, trao đổi lĩnh hội kiến thức Nghiên cứu sinh trân trọng cảm ơn đến Hội đồng Khoa học góp ý nhiều để nội dung luận án hồn thiện tốt Để có ngày hôm nay, nghiên cứu sinh cảm ơn bố mẹ người thân gia đình ln thương yêu, động viên khích lệ vật chất tinh thần Nghiên cứu sinh trân trọng cảm ơn tất người ủng hộ, tạo điều kiện động viên khích lệ tinh thần để luận án hoàn thành Nghiên cứu sinh Nguyễn Quốc Trinh ii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CHỮ VIẾT TẮT vi DANH MỤC BẢNG BIỂU vii DANH MỤC HÌNH VẼ viii MỞ ĐẦU 1 Tính cấp thiết đề tài Mục tiêu nghiên cứu Nội dung nhiệm vụ nghiên cứu Phạm vi nghiên cứu Phương pháp nghiên cứu Ý nghĩa khoa học thực tiễn Các đóng góp luận án Bố cục luận án CHƯƠNG TỔNG QUAN VỀ THỰC TRẠNG SỰ CỐ DẦU TRÀN VÀ CÁC NGHIÊN CỨU DẦU TRÀN TRÊN BIỂN 1.1 Thực trạng cố tràn dầu biển giới, khu vực Biển Đông biển Việt Nam 1.1.1 Hiện trạng cố dầu tràn giới 1.1.2 Hiện trạng cố dầu tràn khu vực Biển Đông biển Việt Nam 1.1.3 Nguy xảy cố dầu tràn khu vực biển Việt Nam Biển Đông 1.2 Tình hình nghiên cứu mơ dầu tràn biển 11 iii 1.2.1 Tình hình nghiên cứu mơ dầu tràn biển giới 12 1.2.2 Tình hình nghiên cứu mô dầu tràn biển nước 14 1.3 Tình hình nghiên cứu khả tìm kiếm nguồn thải 17 1.3.1 Tình hình nghiên cứu khả tìm kiếm nguồn thải giới 17 1.3.2 Tình hình nghiên cứu khả tìm kiếm nguồn thải dầu nước 18 1.4 Kết luận chương 19 CHƯƠNG NGHIÊN CỨU QUÁ TRÌNH LAN TRUYỀN DẦU VÀ KHẢ NĂNG TÌM KIẾM NGUỒN THẢI DẦU 21 2.1 Phương pháp nghiên cứu lan truyền biến đổi dầu 21 2.1.1 Cơ sở lý thuyết lan truyền biến đổi dầu xuôi thời gian 23 2.1.2 Cơ sở lý thuyết lan truyền biến đổi dầu ngược thời gian 29 2.2 Tham số hóa q trình động lực phong hố dầu tràn 34 2.2.1 Các trình động lực của dầu tràn 35 2.2.2 Các q trình phong hóa dầu của dầu tràn 39 2.3 Phương pháp rời rạc hoá, tiêu đánh giá sơ đồ thuật toán 51 2.3.1 Phương pháp giải số 51 2.3.2 Phương pháp rời rạc hoá với sơ đồ sai phân hữu hạn 53 2.3.3 Các sơ đồ cấu trúc thuật toán 64 2.3.4 Chỉ tiêu đánh giá sai số 69 2.4 Kết luận chương 72 CHƯƠNG ĐÁNH GIÁ, PHÂN TÍCH VÀ THẢO LUẬN KẾT QUẢ NGHIÊN CỨU ÁP DỤNG CHO KHU VỰC BIỂN ĐÔNG 73 3.1 Các bước thực thơng tin liệu phục vụ tính toán 73 3.1.1 Các bước thực 73 3.1.2 Thông tin liệu phục vụ làm đầu vào cho tính tốn 74 iv 3.2 Kết nghiên cứu đánh giá kiểm nghiệm xuôi thời gian 82 3.2.1 Ðánh giá hiệu ứng của sơ đồ sai phân luới kết hợp luới 82 3.2.2 Đánh giá kiểm ngiệm hiệu ứng của thành phần tham gia 84 3.2.3 Đánh giá kiểm nghiệm hiệu ứng khuếch tán số xuôi thời gian 88 3.3 Kết nghiên cứu đánh giá thử nghiệm xuôi thời gian 92 3.3.1 Đánh giá thử nghiệm điều kiện môi trường đồng xuôi thời gian 92 3.3.2 Đánh giá thử nghiệm điều kiện mơi trường trung bình tháng xi thời gian 95 3.3.3 Đánh giá thử nghiệm điều kiện môi trường thực xuôi thời gian 99 3.4 Kết nghiên cứu khả tìm kiếm nguồn thải ngược thời gian 105 3.4.1 Một số kết mô ngược thời gian giới 105 3.4.2 Đánh giá thử nghiệm ngược thời gian điều kiện môi trường đồng 107 3.4.4 Đánh giá thử nghiệm điều kiện môi trường biến đổi thực ngược thời gian 114 3.4.5 Đánh giá kiểm nghiệm hiệu ứng khuếch tán số ngược thời gian 122 3.5 Đánh giá độ ổn định mơ hình xi – ngược theo thời gian 124 3.6 Đánh giá độ ổn định mơ hình ngược - xi thời gian 128 3.7 Kết luận chương 131 KẾT LUẬN VÀ KIẾN NGHỊ 132 Kết luận 132 Kiến nghị 132 DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC CỦA TÁC GIẢ CÓ LIÊN QUAN ĐẾN LUẬN ÁN 134 TÀI LIỆU THAM KHẢO 135 PHỤ LỤC 150 v DANH MỤC CHỮ VIẾT TẮT API: Viện Dầu khí Hoa Kỳ (American Petroleum Institute) NOAA: Cục Khí tượng Đại dương Quốc gia Hoa Kỳ (National Oceanic and Atmospheric Administration) ITOPF: Liên đồn Ứng phó nhiễm biển Quốc tế (International Tanker Owners Pollution Federation Ltd) ASCE: Hiệp hội kỹ sư dân dụng Mỹ (American Society of Civil Engineers) CFSR: Hệ thống liệu khí hậu tái phân tích (Climate Forecast System Reanalysis) POM: Mơ hình Đại dương Trường Đại học Princeton (Princeton Ocean Model) SWAN: Mô hình hình mơ sóng ven bờ (Simulating WAves Nearshore) LPTMs: Mơ hình Lagrangian theo dõi quỹ đạo hạt (Lagrangian Particle-Tracking Models) BITT: Quỹ đạo ngược thời gian (Backward-in-Time-Trajectory) LCSs: Cấu trúc nhánh kết hợp Lagrangian (Lagrangian Coherent Structures) BP : Dầu khí Anh quốc (British Petroleum) nnk: Nhiều người khác vi DANH MỤC BẢNG BIỂU Bảng 1.1 Thống kê số vụ tràn dầu lớn hàng hải từ 1970 đến 2015 Bảng 2.1 Đặc trưng mật độ 15oC loại dầu 25 Bảng 2.2 Tổng hợp thông số đặc trưng tham số hố q trình động lực, phong hố tính chất dầu 50 Bảng 2.3 Hệ thống trình hố mơ hình dầu tràn vật thể trơi biển 71 Bảng 3.1 Thông tin liệu đầu vào phục vụ tính tốn lan truyền biến đổi dầu tràn 80 Bảng 3.2 Kết tính tốn theo số NASH, PBIAS, RMSE hàm tương quan trường hợp thành phần tham gia 85 Bảng 3.3 Kết mơ diện tích lớp dầu tràn mặt biển với kịch bước lưới khác 89 Bảng 3.4 Kết tính tốn theo số NASH, PBIAS, RMSE hàm tương quan mô dầu tràn theo thời gian với yếu tố môi trường chế độ tháng 96 Bảng 3.5 Kết tính toán theo số NASH, PBIAS, RMSE hàm tương quan lan truyền dầu với yếu tố môi trường (CFSR POM) theo thời gian thực 100 Bảng 3.6 So sánh đánh giá kết lan truyền dầu thải ngược thời gian theo số NASH, PBIAS, RMSE hàm tương quan 108 Bảng 3.7 Kết tính tốn theo số NASH, PBIAS, RMSE hàm tương quan lan truyền dầu theo thời gian với yếu tố môi trường chế độ tháng 111 Bảng 3.8 Kết tính tốn theo số NASH, PBIAS, RMSE hàm tương quan lan truyền dầu với yếu tố môi trường (CFSR POM) theo thời gian thực 115 vii DANH MỤC HÌNH VẼ Hình 1.1 Một số hình ảnh tàu bị tai nạn gây cố tràn dầu giới Hình 1.2 Đồ thị phân bố số vụ tràn dầu quốc gia có cố dầu tràn nhiều giới giai đoạn 1992-2011 Hình 1.3 Sơ đồ phân bố tàu hoạt động có phát báo Biển Đông lân cận giai đoạn 1995-2005 Hình 1.4 Sơ đồ phân bố cố tràn dầu khu vực Biển Đông lân cận giai đoạn 1976-2005 Hình 1.5 Sơ đồ phân bố mật độ dầu (triệu tấn/năm) vận chuyển qua tuyến hàng hải khu vực Biển Đông lân cận năm 2005 Hình 1.6 Bản đồ phân bố vệt dầu tràn từ ảnh vệ tinh Biển Đông năm 2007 2008 Hình 1.7 Phân bố tiềm dầu khí khu vực Biển Đông lân cận giai đoạn 2006-2007 10 Hình 1.8 Giao diện website quản lý tàu vận tải dầu Cục Hàng hải Việt Nam 10 Hình 1.9 Hệ thống giao thơng hàng hải quốc tế vận chuyển dầu biển 11 Hình 2.1 Các trình vật lý, phong hóa biến đổi dầu tràn biển 21 Hình 2.2 Sơ đồ nguyên lý trình lan truyền biến đổi tính chất dầu tràn 23 Hình 2.3 Sơ đồ mơ tả q trình tác dụng lực làm di chuyển xi- ngược theo thời gian 29 Hình 2.4 Sơ đồ mơ tả q trình khuếch tán vật chất xi- ngược theo thời gian 30 Hình 2.5 Sơ đồ mơ tả q trình bảo tồn vật chất xi- ngược theo thời gian 30 Hình 2.6 Sơ đồ trình động lực phong hoá biến đổi theo thời gian 34 Hình 2.7 Biến đổi hệ số khuếch tán ngang theo thời gian 38 Hình 2.8a Hình ảnh nhũ tương biển 39 Hình 2.8b Biến đổi tỷ lệ phần dầu bị nhũ tương (Fw) theo thời gian 39 Hình 2.9a Tỷ lệ bay nhiệt độ 15oC cho số loại dầu khác 40 Hình 2.9b Biến đổi tỷ phần bay (Fe) theo thời gian 40 viii Phụ lục 2.1 Cấu trúc chương trình tính tốn lan truyền, phong hố biến đổi dầu xi theo thời gian: + Chương trình viết sau: !************************************************************* Program Oil_Euler Implicit none Include 'Oil.h' !=========================================================== 1000 write(*,*) 'Oil Spill Running .' ! Initialise: 1001 Call Initialise ! .Input Information Reading: 1002 Call Input_Read ! Input of constants: 1003 Call Constants_Input ! Inputs of oil' 1004 Call setup_oil ! Setup time 1005 CALL SET_TIME CALL Julian ! Update time 1006 CALL SET_TIME CALL Julian ! Initialise Condition: 1007 ttime = dem = in_month = imonth ! Begining for Computed 13579 CALL SET_TIME CALL Julian ! - Input of filenames and constants: 154 ! 'Environments_Compute 1008 Call Input_Data ! Oil_Spill_Compute 1018 Call Oil_Compute ! Output Oil Spill 1019 Call Outputs c 1020 isecond = isecond + dto ttime = ttime + dto if(INday.lt.Ndaye) goto 13579 C End !************************************************************* Phụ lục 2.2 Cấu trúc chương trình tính tốn lan truyền, phong hố biến đổi dầu xi theo thời gian !**************************************************************** Subroutine Oil_Compute 1000 Include 'Oil.h' ! - Xử lý trình sai phân bậc cho nồng độ độ dày lớp dầu 1001 Call Deflnite_Scheme(im,jm,Coil ,dx,dy,doilx ,doily ) 1002 Call Deflnite_Scheme(im,jm,Ho ,dx,dy,dhox ,dhoy ) ! - Tính tốn q trình bình lưu -1003 Call Advecsion ! - Tính tốn q trình khuếch tán 1004 Call Disffusion ! - Tính tốn q trình xáo trộn thẳng đứng 1005 Call Vertical_Velocity (nếu có) ! - Tính tốn q trình lan truyền học 1006 Call Spreading ! - Xử lý trình sai phân bậc cho nồng độ độ dày lớp dầu - 155 1007 Call Deflnite_Scheme(im,jm,doilx,dx,dy,doilxx,doilxy) 1008 Call Deflnite_Scheme(im,jm,doily,dx,dy,doilyx,doilyy) 1009 Call Deflnite_Scheme(im,jm,dhox,dx,dy,dhoxx,dhoxy) 1010 Call Deflnite_Scheme(im,jm,dhoy,dx,dy,dhoyx,dhoyy) ! - Tính tốn q trình nguồn dầu (nếu có) -1011 Call Oil_Source ! - Xác định lại mật độ dầu (nếu cần) -1012 Call Update_Density ! - Tính tốn lan truyền dầu theo thời gian sai phân trung tâm 1013 Call Deflnite_Scheme_Time ! - Xác định mật độ dầu -1014 Call Update_Density ! - Xác định nhiệt độ dầu 1015 Call Oil_Temperature ! - Xác định độ nhớt dầu -1016 Call Oil_Viscosity ! - Xác định sức cắng bề mặt dầu -1017 Call Oil_Surface_Tension ! - Xác định tỷ phần bay dầu 1018 Call Evaporation ! - Xác định tỷ phần nhũ tương dầu 1019 Call Emulsification ! - Xác định tỷ phần dầu sóng -1020 Call Dispersion_wave ! - Xác định tỷ phần phân tán thẳng đứng dầu -1021 Call Dispersion ! - Xác định tỷ lệ hoà tan dầu -1022 Call Dissolution ! - Xác định tỷ lệ bị phân tán trầm tích dầu (nếu có) -1023 Call Sedimention ! - Xác định tỷ phần tương tác bờ bãi biển (nếu có) - 156 1024 Call Shoreline ! - Xác định tỷ phần ơ-xy hố lớp dầu 1025 Call Oil_Oxy ! - Xác định tỷ phần phân huỷ sinh học lớp dầu -1026 Call Biological ! - Xác định biết đổi tính chất dầu mật độ, thể tích, độ nhớt, ứng suất bề mặt dầu -1017 Call Oil_Properties ! - Xác định lại trạng thái cho qua stình tính tốn dầu -1018 Call Convert_Contants ! Kết thúc vịng tính Return 1019 End ****************************************************************** Phụ lục 2.3 Cấu trúc file điều khiển đầu vào tính tốn lan truyền, phong hố biến đổi dầu xuôi theo thời gian Caculation Zone: 99 ! West 121 ! East ! South 23 ! North Start time: 2007 ! years 02 ! months 01 ! days 00 ! hours 00 ! mints 0.0 ! sects End time: 2007 ! yeare 02 ! monthe 20 ! daye 00 ! houre 00 ! minte 0.0 ! secte Step time: 600 ! dto [second]: External (2-D) time step Outputs Time: ! itime_hotstart [Integer]: Output write to hot_start time step (hour) 157 ! itime_grid [Integer]: Output write to fine_grid time step (hour) ! itime_tecpl [Integer]: Output write to fine_grid time step for tecplot form (hour) ! igrid Outputs File: /Outputs/Test1_realfb Caculation kinds: ! Index to indicate whether run with wind (0/1: no/yes oil input file) ./Forward/Vetdau2007p.txt ! Index to indicate whether run with wind (0/1: no/yes wind input file) ./Forward/Wind_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Pressure input file) ./Forward/Pressure_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Temperature Air input file) ./Forward/TempAir_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Long Wave Radiation input file) ! Index to indicate whether run with wind (0/1: no/yes Sorf Wave Radiation input file) ! Index to indicate whether run with wind (0/1: no/yes Rain input file) ./Forward/Rain.inc ! Index to indicate whether run with wind (0/1: no/yes Cloud input file) ./Forward/Cloud.inc ! Index to indicate whether run with wind (0/1: no/yes Relative Humidity input file) ./Forward/Rh2m.inc ! Index to indicate whether run with tidal (0/1: no/yes tidal input file) ! Index to indicate whether run with level waterl (0/1: no/yes level input file) ./Forward/Level_POM.inc ! Index to indicate whether run with wind (0/1: no/yes Current input file) ./Forward/Currents_POM.inc ! Index to indicate whether run with wind (0/1: no/yes Temperature Water input file) ./Forward/Temps_POM.inc ! Index to indicate whether run with wind (0/1: no/yes Salt input file) ./Forward/Salt_POM.inc ! Index to indicate whether run with wind (0/1: no/yes Wave input file) ! Index to indicate whether run with wind (0/1: no/yes River input file) ! Index to indicate whether run to start from restart file (0/1: no/yes restart input file) Phụ lục 2.4 Cấu trúc file dầu tràn đầu vào tính tốn lan truyền, phong hố biến đổi dầu xi theo thời gian (“./Forward/Vetdau2007p.txt”) id Year mon day hh sec 229 2007 01 16 15 07 34.000 229 2007 01 16 15 07 34.000 229 2007 01 16 15 07 34.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 longitude latitude Volume(m3) Time(hour) 106.918785 8.132027 0.321688294E+00 106.915619 8.134459 0.321688294E+00 106.915619 8.135269 0.321688294E+00 109.896767 13.923503 0.172506022E+01 109.892715 13.926204 0.172506022E+01 109.891022 13.926879 0.172506022E+01 109.895416 13.858673 0.172506022E+01 109.853668 13.690890 0.474384129E+00 109.859909 13.704406 0.474384129E+00 109.862686 13.708912 0.474384129E+00 158 Phụ lục 2.5 Cấu trúc chương trình tính tốn lan truyền, phong hố biến đổi dầu ngược theo thời gian: + Chương trình viết sau: !************************************************************* Program OilBack_Euler Implicit none Include 'Oil.h' !=========================================================== ! Back Computed 2000 write(*,*) 'Oil Spill Back Running .' ! Initialise: 2001 Call Initialise ! .Input Information Reading: 2002 Call InputBack_Read ! Input of constants: 2003 Call Constants_Input ! Setup oil back : 2004 Call setup_oil_Back ! Setup time back : 2005 CALL SET_TIME CALL Julian ! Update time 2006 CALL SET_TIME CALL Julian ! Initialise Condition: 2007 ttime = dem = in_month = imonth ! Begining for Computed 97531 CALL SET_TIME CALL Julian 159 ! 'Environments_Compute 2008 Call Input_Data ! Oil_Spill_Back_Compute 2018 Call Oil_BackCompute Output Oil Spill Back 2019 Call Outputs_Back c 2020 isecond = isecond + dto ttime = ttime + dto if(INday.gt.Ndays) goto 97531 C End !************************************************************* Phụ lục 2.6 Cấu trúc chương trình tính tốn lan truyền, phong hố biến đổi dầu ngược theo thời gian !**************************************************************** Subroutine Oil_BackCompute 2000 Include 'Oil.h' ! - Xử lý trình sai phân bậc cho nồng độ độ dày lớp dầu 2001 Call Deflnite_Scheme(im,jm,Coil ,dx,dy,doilx ,doily ) 2002 Call Deflnite_Scheme(im,jm,Ho ,dx,dy,dhox ,dhoy ) ! - Tính tốn q trình bình lưu -2003 Call Advecsion ! - Tính tốn q trình khuếch tán 2004 Call Disffusion ! - Tính tốn q trình xáo trộn thẳng đứng 2005 Call Vertical_Velocity (nếu có) ! - Tính tốn q trình lan truyền học 2006 Call Spreading ! - Xử lý trình sai phân bậc cho nồng độ độ dày lớp dầu - 160 2007 Call Deflnite_Scheme(im,jm,doilx,dx,dy,doilxx,doilxy) 2008 Call Deflnite_Scheme(im,jm,doily,dx,dy,doilyx,doilyy) 2009 Call Deflnite_Scheme(im,jm,dhox,dx,dy,dhoxx,dhoxy) 2010 Call Deflnite_Scheme(im,jm,dhoy,dx,dy,dhoyx,dhoyy) ! - Tính tốn q trình nguồn dầu (nếu có) -2011 Call Oil_Source ! - Xác định lại mật độ dầu (nếu cần) -2012 Call Update_Density ! - Tính toán lan truyền dầu theo thời gian sai phân trung tâm 2013 Call Deflnite_Scheme_Time ! - Xác định mật độ dầu -2014 Call Update_Density ! - Xác định nhiệt độ dầu 2015 Call Oil_Temperature ! - Xác định độ nhớt dầu -2016 Call Oil_Viscosity ! - Xác định sức cắng bề mặt dầu -2017 Call Oil_Surface_Tension ! - Xác định tỷ lệ bay dầu -2018 Call Evaporation ! - Xác định tỷ lệ nhũ tương dầu 2019 Call Emulsification ! - Xác định tỷ lệ phân tán dầu sóng 2020 Call Dispersion_wave ! - Xác định tỷ lệ phân tán thẳng đứng dầu 2021 Call Dispersion ! - Xác định tỷ lệ hoà tan dầu -1022 Call Dissolution ! - Xác định tỷ lệ bị phân tán trầm tích dầu (nếu có) -2023 Call Sedimention ! - Xác định tỷ lệ bị phân tán bời bãi dầu (nếu có) 2024 Call Shoreline 161 ! - Xác định tỷ lệ bị ơ-xy hố dầu -2025 Call Oil_Oxy ! - Xác định tỷ lệ bị phân huỷ sinh học dầu 2026 Call Biological ! - Xác định biết đổi tinh chất dầu mật độ, thể tích, độ nhớt, ứng suất bề mặt dầu -2017 Call OilBack_Properties ! - Xác định lại trạng thái cho qua stình tính tốn dầu -2018 Call Convert_Contants ! Kết thúc vịng tính Return 2019 End ****************************************************************** Phụ lục 2.7 Cấu trúc file điều khiển đầu vào tính tốn lan truyền, phong hoá biến đổi dầu ngược theo thời gian Caculation Zone: 99.0 ! West 121.0 ! East 1.5 ! South 23.5 ! North Start time: 2007 ! years 02 ! months 10 ! days 00 ! hours 00 ! mints 0.0 ! sects End time: 2007 ! yeare 03 ! monthe 15 ! daye 00 ! houre 00 ! minte 0.0 ! secte Oil Time 2007 ! year0 03 ! month0 15 ! day0 00 ! hour0 162 00 ! minute0 ! second0 Step time: 600 ! dto [second]: External (2-D) time step Outputs Time: ! itime_hotstart [Integer]: Output write to hot_start time step (hour) ! itime_grid [Integer]: Output write to fine_grid time step (hour) ! itime_tecpl [Integer]: Output write to fine_grid time step for tecplot form (hour) ! igrid Outputs File: /Outputs/Real02b Caculation kinds: ! Index to indicate whether run with wind (0/1: no/yes oil input file) ./Backward/Vetdau2007p.txt ! Index to indicate whether run with wind (0/1: no/yes wind input file) ./Backward/Wind_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Pressure input file) ./Backward/Pressure_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Temperature Air input file) ./Backward/TempAir_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Long Wave Radiation input file) ! Index to indicate whether run with wind (0/1: no/yes Sorf Wave Radiation input file) ! Index to indicate whether run with wind (0/1: no/yes Rain input file) ./Backward/Rain.inc ! Index to indicate whether run with wind (0/1: no/yes Cloud input file) ./Backward/Cloud.inc ! Index to indicate whether run with wind (0/1: no/yes Relative Humidity input file) ./Backward/Rh2m.inc ! Index to indicate whether run with tidal (0/1: no/yes tidal input file) ! Index to indicate whether run with level waterl (0/1: no/yes level input file) ./Backward/Level_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Current input file) ./Backward/Currents_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Temperature Water input file) ./Backward/Temps_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Salt input file) ./Forward/Salt_CFSR.inc ! Index to indicate whether run with wind (0/1: no/yes Wave input file) ! Index to indicate whether run with wind (0/1: no/yes River input file) ! Index to indicate whether run to start from restart file (0/1: no/yes restart input file) Phụ lục 2.8 Cấu trúc file dầu tràn đầu vào tính tốn lan truyền, phong hố biến đổi dầu xuôi theo thời gian (“./Backward/Vetdau2007p.txt”) id Year mon day hh sec 229 2007 01 16 15 07 34.000 229 2007 01 16 15 07 34.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 2007 01 19 02 32 47.000 longitude 106.918785 106.915619 109.896767 109.892715 109.853668 109.862686 163 latitude 8.132027 8.134459 13.923503 13.926204 13.690890 13.708912 Volume(m3) 0.321688294E+00 0.321688294E+00 0.172506022E+01 0.172506022E+01 0.474384129E+00 0.474384129E+00 Phụ lục 2.9 Cấu trúc file địa hình phục vụ tính tốn (“Depth.inc”) Depth Information: / / /TOPOs/ EVS_64 grd longmin longmax latmin 90.0 140.0 -10.0 40.0 3201 ! id_style ! id_title latmax nlong nlat 3201 Phụ lục 2.9a Cấu trúc định dạng file địa hình phục vụ tính tốn (“EVS_64.grd”) Phụ lục 2.10 Cấu trúc file gió đầu vào phục vụ tính tốn (“Wind_CFSR.inc”) / / /DATA/meteorological/CFSR/Wind/CFSRWind 00 2007 01 01 00 00 00 ! to start time of wind data: year month day hour minute second 2009 01 01 00 00 00 ! to end time of wind data: year month day hour minute second hour ! to step time of wind field data longmin longmax latmin latmax nlong nlat 95.0 145.0 -15.0 35.0 101 101 164 Phụ lục 2.10a Cấu trúc định dạng file gió đầu vào phục vụ tính tốn thời điềm (“CFSRWind20070215000000”) Phụ lục 2.11 Cấu trúc file nhiệt độ khơng khí đầu vào phục vụ tính tốn (“Temp_CFSR.inc”) / / /DATA/meteorological/CFSR/Tema/CFSRTema 00 2007 01 01 00 00 00 ! to start time of Temp data: year month day hour minute second 2009 01 01 00 00 00 ! to end time of Temp data: year month day hour minute second hour ! to step time of temp field data longmin longmax latmin latmax nlong nlat 95.0 145.0 -15.0 35.0 101 101 Phụ lục 2.11a Cấu trúc định dạng file nhiệt độ khơng khí đầu vào phục vụ tính tốn thời điềm (“CFSRTema20070215000000”) 165 Phụ lục 2.12 Cấu trúc file dịng chảy đầu vào phục vụ tính tốn (“Curent_POM.inc”) / / /Models/POM/POM_EVS/Outputs/FLOW3Ds 00.txt 2007 01 01 00 00 00 ! to start time of current data: year month day hour minute second 2009 01 01 00 00 00 ! to end time of current data: year month day hour minute second hour ! to step time of current field data longmin longmax latmin latmax nlong nlat id(1,2,3,4: colum; 5,6,7,8:table) 98.750 125.000 -3.000 27.000 316 361 Phụ lục 2.12a Cấu trúc định dạng file dòng chảy đầu vào phục vụ tính tốn thời điềm (“FLOW3Ds20070215000000.txt”) Phụ lục 2.13 Cấu trúc file nhiệt đô nước đầu vào phục vụ tính tốn (“Temp_POM.inc”) / / /Models/POM/POM_EVS/Outputs/TEMP3Ds 00.txt 2007 01 01 00 00 00 ! to start time of Temp data: year month day hour minute second 2009 01 01 00 00 00 ! to end time of Temp data: year month day hour minute second hour ! to step time of Tem field data longmin longmax latmin latmax nlong nlat id(1,2,3,4: colum; 5,6,7,8:table) 98.750 125.000 -3.000 27.000 316 361 166 Phụ lục 2.13a Cấu trúc định dạng file nhiệt độ bề mặt nước đầu vào phục vụ tính tốn thời điềm (“TEMAs20070215000000.txt”) Phụ lục 2.14 Cấu trúc dạng file đầu dạng điểm (“Oil_real_Pmax_f.txt”) Phụ lục 2.15 Cấu trúc dạng file đầu dạng trường theo thời điểm (“Test_Real07f_aa1f_20080708120000.txt”) 167 Phụ lục 2.16 Cấu trúc dạng file đầu dạng trường cực đại theo thời điểm (“Test_Real07f_aa1max_f_20080708120000.txt”) Phụ lục 2.17 Cấu trúc dạng file đầu dạng trường theo thời điểm khác (“Test_Real07f_aa1f_20080708120000.plt”) 168 ... Nghiên cứu trình lan truyền, phong hố biến đổi tính chất dầu; 4, Nghiên cứu khả tìm kiếm nguồn phát thải dầu gây ô nhiễm; 5, Mô lan truyền dầu khả tìm kiếm nguồn thải dầu Phạm vi nghiên cứu Khu... khả tìm kiếm nguồn thải dầu (bài tốn ngược thời gian) 20 CHƯƠNG NGHIÊN CỨU QUÁ TRÌNH LAN TRUYỀN DẦU VÀ KHẢ NĂNG TÌM KIẾM NGUỒN THẢI DẦU 2.1 Phương pháp nghiên cứu lan truyền và biến đổi dầu. .. hình nghiên cứu khả tìm kiếm nguồn thải dầu nước 18 1.4 Kết luận chương 19 CHƯƠNG NGHIÊN CỨU QUÁ TRÌNH LAN TRUYỀN DẦU VÀ KHẢ NĂNG TÌM KIẾM NGUỒN THẢI DẦU 21 2.1 Phương pháp nghiên