Xây dựng thuật toán và sơ đồ khối của chương trình dự báo động đất theo mô hình thống kê

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Trần Việt Phƣơng XÂY DỰNG THUẬT TỐN VÀ SƠ ĐỒ KHỐI CỦA CHƢƠNG TRÌNH DỰ BÁO ĐỘNG ĐẤT THEO MƠ HÌNH THỐNG KÊ LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2012 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Trần Việt Phƣơng XÂY DỰNG THUẬT TOÁN VÀ SƠ ĐỒ KHỐI CỦA CHƢƠNG TRÌNH DỰ BÁO ĐỘNG ĐẤT THEO MƠ HÌNH THỐNG KÊ Chun ngành : Vật Lý Địa Cầu Mã số : 604415 LUẬN VĂN THẠC SĨ KHOA HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: TSKH Ngô Thị Lƣ Hà Nội - Năm 2012 Luận văn Thạc sĩ môc lôc MỞ ĐẦU Chƣơng 1: TỔNG QUAN 1.1 Tổng quan tình hình nghiên cứu dự báo động đất giới 1.2 Tình hình nghiên cứu dự báo động đất Việt Nam 11 Chƣơng 2: PHƢƠNG PHÁP DỰ BÁO ĐỘNG ĐẤT THEO MÔ HÌNH THỐNG KÊ 19 2.1 Cơ sở lý thuyết phƣơng pháp 19 2.2 Qui trình dự báo động đất theo mơ hình thống kê 26 Chƣơng 3: THIẾT LẬP CHƢƠNG TRÌNH DỰ BÁO ĐỘNG ĐẤT THEO MƠ HÌNH THỐNG KÊ 31 3.1 Thuật toán sơ đồ khối chƣơng trình dự báo động đất theo mơ hình thống kê: 31 3.2 Ngôn ngữ lập trình 35 3.3 Chƣơng trình 36 3.4 Hƣớng dẫn sử dụng phần mềm 37 3.5 Các giao diện chƣơng trình 40 3.6 Code chƣơng trình 41 Chƣơng 4: ÁP DỤNG CHƢƠNG TRÌNH DỰ BÁO ĐỘNG ĐẤT THEO MƠ HÌNH THỐNG KÊ THỬ NGHIỆM DỰ BÁO CHO CÁC KHU VỰC CỤ THỂ 45 4.1 Các tài liệu sử dụng 45 4.2 Áp dụng thử nghiệm chƣơng trình khu vực tây bắc việt nam 45 4.3 Áp dụng thử nghiệm dự báo danh mục động đất đông nam 50 4.4 So sánh kết dự báo trƣờng hợp lựa chọn xác xuất tin cậy dự báo khác 53 4.5 Nhận xét 55 KẾT LUẬN 56 Trần Việt Phương Luận văn Thạc sĩ TÀI LIỆU THAM KHẢO 58 PHỤ LỤC 65 PHỤ LỤC 72 Bảng ký hiệu chữ viết tắt DMĐĐ Danh mục động đất DMĐĐĐL Trần Việt Phương Danh mục động đất độc lập Luận văn Thạc sĩ Các cơng trình cơng bố Trong trình học tập làm việc viện Vật lý địa cầu để hoàn thành luận văn, tác giả luận văn tham gia đồng tác giả cơng trình đƣợc cơng bố dƣới đây: Ngô Thị Lƣ, Trần Việt Phƣơng (2009), “Tách nhóm tiền chấn, dƣ chấn từ danh mục động đất khu vực Đông Nam Á (chu kỳ 1278-2008) phƣơng pháp cửa sổ không gian thời gian” Các khoa học Trái đất 31(1), Hà Nội, 2009.Tr 35- 43 Burmin V.Yu., Ngô Thị Lƣ, Trần Việt Phƣơng (2009), “Đánh giá tính hiệu hệ thống trạm địa chấn có Việt Nam” Các thiết bị địa chấn, Viện Hàn lâm khoa học Nga, 45(1), Moscow, 2009 Tr 44-61 (Tiếng Nga) Rodkin M.V.*, Ngo Thi Lu **, Pisarenko V.F.*, Tran Viet Phuong** and Vu Thi Hoan **(2010), “Change in the regime of growth of cumulative seismic energy with time: examination from the regional catalogue of Vietnam”, 8th General Assembly of Asian Seismological Commision (ASC 2010) (*International Institute of Earthquake Prediction Theory and Mathematical Geophysics Russian Academy of Sciences (IIEPT RAS); **- Institute of Geophysics VAST (Vietnam Academy of Science and Technology) V.Yu Burmin*, Ngo Thi Lu**, Tran Viet Phuong**(2010), “Design of an optimal network of seismic stations in North Vietnam” 8th General Assembly of Asian Seismological Commision (ASC 2010) *Schmidt Institute of Physics of the Earth, Russian Academy of Sciences, Moscow, Russia; **Institute of Geophysics of Vietnamese Academy of Science and technology, Hanoi Ngo Thi Lu, Nguyen Anh Quan, Tran Viet Phuong (2010), “Establishing a computer program for earthquake prediction on the vietnamese territory and adjacent regions by zoning of Earth’s crust types” Journal of Geology, series B(35-36/2010) tr 111-130) Trần Việt Phương Luận văn Thạc sĩ Ngô Thị Lƣ, Trần Việt Phƣơng (2012) “Thiết lập chƣơng trình dự báo động đất theo mơ hình thống kê”, Địa chất, Số 331-332; 5-8/2012, tr 40-49 ISSN 08667381 Ngô Thị Lƣ, Trần Việt Phƣơng (2012) “Cải biên thuật toán xây dựng sơ đồ khối chƣơng trình dự báo động đất cực đại phƣơng pháp vật lý kiến tạo”, Tc Địa chất 331-332; 5-8/2012, tr 50-58 ISSN 0866-7381 Ngô Thị Lƣ, Trần Việt Phƣơng (2012) “Về cách tiếp cận đề xây dựng thuật tốn quy trình dự báo động đất theo mơ hình thống kê”, Các khoa học Trái đất 34(3), Hà Nội, 2012, tr 3-7 Ngô Thị Lƣ (Chủ nhiệm đề tài), Trần Việt Phƣơng, Phùng Thị Thu Hằng, Nguyễn Hữu Tuyên nnk (2011) Đánh giá tiềm địa chấn lãnh thổ Việt Nam theo tổ hợp tài liệu địa chất-địa vật lý địa chấn, Nhiệm vụ hợp tác khoa học quốc tế hai viện Hàn lâm Khoa học Việt Nam viện HLKH Liên bang Nga theo Nghị định thƣ cấp Nhà nƣớc (giai đoạn 2008-2010) Hà Nội, 2011 Báo cáo tổng kết đề tài 163tr Trần Việt Phương Luận văn Thạc sĩ MỞ ĐẦU Những năm gần đây, thảm hoạ động đất, sóng thần tồn cầu ngày gia tăng, khu vực Đông Nam Á Do có vị trí địa lý nằm nơi tiếp giáp vành đai hoạt động địa chấn lớn liên quan với vành đai phá huỷ kiến tạo hoạt động mạnh mẽ hành tinh, nên Đông Nam Á chịu ảnh hƣởng hoạt động phá huỷ kiến tạo mạnh mà chịu độ nguy hiểm động đất sóng thần cao Đặc biệt, thảm hoạ động đất sóng thần Sumatra ngày 26.12.2004, thảm hoạ động đất Tứ Xuyên (12.05.2008) động đất (4.2010) (Trung Quốc), thảm họa động đất Tōhoku (Nhật Bản) vào ngày 11 tháng năm 2011 gây tổn thất vô nghiêm trọng ngƣời, về phá huỷ môi trƣờng Gần nhất, tƣợng động đất liên tục xảy khu vực đập thủy điện sông Tranh gây hoang mang dƣ luận nƣớc nói chung tỉnh Quảng Nam nói riêng Tình hình thực tế nêu làm cho vấn đề dự báo động đất, sóng thần, vốn vấn đề nan giải ln mang tính thời qui mơ tồn cầu, trở nên cấp thiết đƣợc nhà khoa học vô quan tâm Lãnh thổ VN tham gia vào thành phần ĐNA, nên nhiều chịu ảnh hƣởng vùng hoạt động kiến tạo vùng nguy hiểm địa chấn, đƣợc đặc trƣng tính địa chấn hoạt động kiến tạo tích cực nhƣ Do đó, nghiên cứu dự báo động đất lãnh thổ Việt Nam vùng lân cận vấn đề có ý nghĩa khoa học thực tiễn cao mà cịn nhiệm vụ có tính cấp thiết thực tế địa chấn Việt Nam Để giải vấn đề mang tính thời nhƣ lãnh thổ bất kỳ, nhiệm vụ quan trọng hàng đầu cần xây dựng thiết lập chƣơng trình, cho phép dự báo thời gian, địa điểm độ mạnh trận động đất xảy tƣơng lai gần nhằm xây dựng biện pháp phòng tránh, giảm thiểu tối đa thiệt hại Nhiệm vụ nhƣ có ý nghĩa bối cảnh Viện VLĐC thực Dự án tăng cƣờng trang thiết bị trạm, phục vụ báo tin động đất cảnh báo sóng thần Để giải vấn đề nêu trên, nhiệm vụ có tính cấp thiết cần tiến hành “Xây dựng thuật toán sơ đồ khối chương trình dự báo động Trần Việt Phương Luận văn Thạc sĩ đất theo mơ hình thống kê” Đó lý mà tác giả luận văn chọn tên đề tài nghiên cứu Mục tiêu luận văn Xây dựng thuật toán thiết lập chƣơng trình dự báo động đất theo mơ hình thống kê Áp dụng thử nghiệm chƣơng trình đƣợc thiết lập khu vực Tây Bắc Việt Nam khu vực Đông Nam Á để so sánh kết dự báo với với kiện xảy thực tế nhằm kiểm tra tính đắn chƣơng trình Nhiệm vụ luận văn Tìm hiểu phƣơng pháp dự báo động đất theo mơ hình thống kê sở tổ hợp tài liệu địa vật lý địa chấn Thu thập tài liệu địa vật lý địa chấn cho khu vực nghiên cứu, phân tích, lựa chọn chỉnh lý số liệu phục vụ hƣớng nghiên cứu Xây dựng thuật toán, sơ đồ khối, lựa chọn ngơn ngữ lập trình phù hợp thiết lập chƣơng trình phân loại vỏ Trái đất máy tính Áp dụng thử nghiệm chƣơng trình đƣợc thiết lập khu vực khác nhau, nhận xét đánh giá khả ứng dụng chƣơng trình Ý nghĩa khoa học thực tiễn luận văn Đã chuyển giao, áp dụng cách sáng tạo phƣơng pháp dự báo mơ hình thống kê, góp phần giải nhiệm vụ dự báo đƣợc thời gian, tọa độ magnitude động đất Việt Nam lân cận Thời gian gần có nhiều trận động đất mạnh xảy ra, gây nên tổn thất vô nặng nề ngƣời Vì vậy, xây dựng phần mềm có khả dự báo trƣớc đƣợc động đất yêu cầu vơ cấp thiết có ý nghĩa khoa học có ý nghĩa thực tiễn Trần Việt Phương Luận văn Thạc sĩ Những nội dung đƣợc thực luận văn góp phần thiết thực vào việc giải nhiệm vụ Đề tài độc lập cấp Viện khoa học Việt Nam, mã số: VAST.ĐL.09/11-12: “Xây dựng chương trình thử nghiệm dự báo ngắn hạn động đất sở mơ hình thống kê kết hợp sử dụng phương pháp Vật lý kiến tạo, áp dụng lãnh thổ Việt Nam lân cận” Cấu trúc luận văn Luận văn gồm phần mở đầu, chƣơng, kết luận, tài liệu tham khảo phụ lục Toàn nội dung nêu đƣợc trình bày 76 trang đánh máy khổ A4, với 10 hình vẽ bảng biểu minh họa (không kể phần phụ lục) Phần mở đầu gồm trang trình bày tính cấp thiết lý chọn đề tài: “Xây dựng thuật tốn sơ đồ khối chương trình dự báo động đất theo mơ hình thống kê” Trong phần cịn trình bày mục tiêu , nhiệm vụ, kết nhận đƣợc, điểm mới, ý nghĩa khoa học thực tiễn luận văn Chƣơng 1: gồm 10 trang, giới thiệu tổng quan tình hình nghiên cứu dự báo động đất giới Việt Nam, phƣơng pháp xác định magnutide cực đại động đất vùng phát sinh động đất mạnh Chƣơng 2: gồm 12 trang với hình vẽ trình bày sở lý thuyết phƣơng pháp quy trình dự báo động đất Chƣơng 3: gồm 14 trang với hình vẽ Nội dung chƣơng trình bày thuật tốn, sơ đồ khối, giới thiệu ngơn ngữ lập trình đƣợc sử dụng, giao diện chƣơng trình hƣớng dẫn sử dụng chƣơng trình Chƣơng 4: gồm 11 trang với hình vẽ bảng biểu trình bày kết áp dụng thử nghiệm chƣơng trình lập để dự báo động đất cho khu vực Tây Bắc Việt Nam khu vực Đông Nam Á Trên sở kết tác giả đƣa nhận xét đánh giá khả áp dụng ƣu điểm chƣơng trình lập Trần Việt Phương Luận văn Thạc sĩ Chƣơng TỔNG QUAN 1.1 Tổng quan tình hình nghiên cứu dự báo động đất giới: Từ lâu nhà khoa học thống kê tổng kết đƣợc nguyên nhân gây động đất, nơi thƣờng xun xảy động đất tính tốn đƣợc tham số trận động đất diễn Tuy nhiên, vấn đề dự báo sớm động đất xảy đâu, có độ lớn điều mà nhà khoa học hƣớng đến tốn chƣa có lời giải trọn vẹn Qua nhiều kỷ, ngƣời ta dựa sở khác nhau, từ hoạt động khác thƣờng số lồi vật tới hình thù kỳ lạ đám mây, biến đổi đột ngột mực nƣớc giếng, hay thay đổi hàm lƣợng radon hydro đất đá để tìm cách dự báo động đất Một lần ngƣời ta dự báo đƣợc xác trận động đất Haicheng, Trung Quốc, năm 1975 Lệnh di tản đƣợc phát ngày trƣớc trận động đất mạnh 7,3 độ Richter tàn phá thành phố Trong nhiều tháng trƣớc đó, ngƣời ta ghi đƣợc hàng loạt trận động đất nhỏ, với thay đổi mực nƣớc ngầm dâng lên bề mặt địa hình Tuy nhiên, trận động đất có dấu hiệu tiền báo nhƣ Sau thành công Haicheng, nhà địa chấn học Trung Quốc dự báo trận động đất với sức hủy diệt lớn Tangshan năm 1976 Với cƣờng độ 7,6 độ Richter, trận động đất cƣớp sinh mạng 250 nghìn ngƣời Ngày nay, với tiến khoa học kỹ thuật, cơng nghệ dự báo thiên tai ngày đƣợc hồn thiện Các trận động đất thƣờng kết chuyển động phận đứt gãy vỏ Trái đất, cấu tạo chủ yếu từ chất rắn Tuy chậm nhƣng mặt đất chuyển động động đất xảy ứng suất (nội lực phát sinh vật thể biến dạng tác nhân bên tác dụng) cao sức chịu đựng đất đá Các nhà khoa học nhận thấy việc đo thay đổi đoạn đứt gãy khó nhiều so với việc đo biến thiên ứng suất, đặc biệt Trần Việt Phương 10 Luận văn Thạc sĩ Nam lân cận, luận án Tiến sĩ Vật lý, Hà nội (161 tr) 35 Trần thị Mỹ Thành, Nguyễn Thị Cẩm nnk (2002) Quản lý danh mục động đất đồ chấn tâm động đất hệ thông tin địa lý, báo cáo đề tài nghiên cứu năm 2002 (Vì tiến phụ nữ) 36 Nguyễn Thị Kim Thoa tập thể tác giả (1999) Cơ sở liệu Vật lý Địa cầu – Số liệu địa chấn 1990-1998 (thông báo động đất danh mục động đất), đề án điều tra 1996-1998 37 Phạm Văn Thục (2007) Địa chấn học động đất Việt Nam, Bộ sách chuyên khảo Tài nguyên thiên nhiên môi trường Việt Nam, Nxb Khoa học tự nhiên công nghệ , Viện KH&CN Việt Nam, 378 tr 38 Cao Đình Triều (1999) “Về số quy luật hoạt động khả dự báo khu vực phát sinh động đất mạnh Việt Nam” Tạp chí Địa chất, loạt A(251), tr 14-21, Hà Nội 39 Cao Đình Triều, Nguyễn Hữu Tuyên, Thái Anh Tuấn (2006) “ Mối quan hệ đặc trƣng cấu trúc vỏ Trái đất hoạt động động đất Tây Bắc Việt Nam”, Tạp chí Các khoa học Trái đất, T.28, tr 155 – 164, Hà Nội 40 Cao Đình Triều, Ngơ Thị Lƣ, Mai Xn Bách nnk (2007) "Dự báo cực đại động đất phần đất liền lãnh thổ Việt Nam sở phân loại dạng vỏ Trái đất", Tuyển tập báo cáo khoa học Hội nghị Khoa học Kỹ thuật Địa Vật lý Việt nam lần thứ 5, NXB KH&KT, Hà Nội, năm 2007, tr 159-171 41 Cao Đình Triều, Rogozhin E.A., Yunga S.I., Ngô Thị Lƣ, Nguyễn Hữu Tuyên, Lê văn Dũng nnk (2009) "Một số kết bƣớc đầu khảo sát dấu vết nghi ngờ hoạt động động đất cổ để lại miền Tây Bắc Bộ, Việt Nam" Địa chất (311), Hà Nội, 2009, tr 1-10 42 Cao Đình Triều (2010) Tai biến động đất Việt Nam, Nhà xuất khoa học kỹ thuật, 303 tr, Hà Nội 43 Nguyễn Đình Xuyên, (chủ biên), Nguyễn Ngọc Thủy nnk (1996) Cơ sở liệu cho giải pháp giảm nhẹ hậu động đất Việt Nam, Báo cáo tổng kết đề tài độc lập cấp nhà nƣớc mã số KT-ĐL 92 – 07, tập, Viện Vật Trần Việt Phương 63 Luận văn Thạc sĩ lý địa cầu, Hà Nội 44 Nguyễn Đình Xuyên (2002) Động đất độ nguy hiểm động đất, tài liệu nội bộ, Lƣu trữ Viện Vật lý địa cầu, Hà Nội 45 Nguyễn Đình Xuyên, , Phạm Đình Nguyên, Phạm Quang Hùng, Nguyễn Ánh Dƣơng (2003) “Thử nghiệm dự báo động đất dài hạn theo quy luật hoạt động tiền chấn”, tạp chí Các khoa học Trái đất, 25(3), tr 193-200, Hà Nội 46 Nguyễn Đình Xuyên (1978) Đánh giá tính động đất lãnh thổ miền Bắc Việt Nam, luận án PTS toán-lý, Viện Vật lý Trái đất, Viện Hàn lâm Khoa học Liên Xô, Matscva 47 Nguyễn Đình Xuyên (1980) Sử dụng tài liệu thực địa về động đất phân vùng động đất lãnh thổ Việt Nam, kết nghiên cứu Viện Vật lý Địa cầu 1979, viện Khoa học Việt Nam, Hà Nội 48 Nguyễn Đình Xuyên, Nguyễn Ngọc Thủy tập thể tác giả (1996) Cơ sở liệu cho giải pháp giảm nhẹ hậu động đất Việt Nam,báo cáo tổng kết đề tài độc lập cấp nhà nƣớc KT-ĐL 92-07 (Phần danh mục động đất), Hà Nội 104 tr 49 Nguyễn Văn Yêm, Nguyễn Thị Cẩm nnk (1998) Động đất khu vực Tây Bắc Việt Nam năm 1998 việc xác định tính phù hợp thang magnitude động đất, viện Vật lý Địa cầu, Hà Nội 50 Zavialov A.D (1984) Độ nghiêng đồ thị lặp lại động đất  dấu hiệu báo trước động đất mạnh Dự báo động đất M ≥5, NXB Dolgin – Dushanbe- Moscow, Tr 173-184 Tiếng Anh 51 Earthquake Hazards program U S Geological Survey (USGS), National Earthquake Information Center World Data Center for Seismology, Denver.(1970-2011) 52 Gumbel E J, (1958) Statistics of Extremes, Columbia Univ Press 53 Gutenberg B Trần Việt Phương and Richter G F, (1954), Seismicity of the Earth and 64 Luận văn Thạc sĩ associated phenomena Princeton Univ Press 54 Ngo Thi Lu (1981), "Seismological Bulletin 1979-1980" National Center for Scientific Research of Vietnam, Hanoi, 1981 55 Ngo Thi Lu (1983), "Seismological Bulletin 1981-1982" National Center for Scientific Research of Vietnam, Hanoi, 1983 56 Ngo Thi Lu (1988), Bulletin of Vietnamese seismological stations in the period from 1981 to 1986 (near earthquakes), Institut of Geophysics, National Center for Scientific Research of Vietnam, Hanoi 57 Ngo Thi Lu (1988), Bulletin of Vietnamese seismological stations in the period from 1983 to 1984 (far earthquakes), Institut of Geophysics //National Center for Scientific Research of Vietnam, Hanoi 58 Ngo Thi Lu (1998) "Strong earthquakes on the territory of South-East Asia and peculiarities in their manifestation" Journal of earthquake prediction research, 7(2) 1998 PP 170 –185 59 Ngo Thi Lu, Gatinsky Y G., Kondorskaia N V (2000) "Seismicity and modem geokinematics of Southeast part of Asia" Reports of the Russian Acad Of Sci, 347(2) Moscow, 2000 C 247-251 (in Russian) 60 Seismological Bulletin of ISC (International Seismological Center) đến 2011 Trần Việt Phương 65 Luận văn Thạc sĩ PHỤ LỤC PHẦN CODE CHÍNH CỦA CHƢƠNG TRÌNH Lớp kiểu liệu trandongdatdubao public class trandongdatdubao { public double fi0, lamda0, dfi, fiop, lamdaop, top, tpr, tmin, tmax, mpr, mmin, mmax, treal, fireal, lamdareal, mreal; public char trt, trp, trm; } Lớp kiểu liệu trandongdat public class trandongdat { public int nam, thang, ngay, gio, phut, x, y; public double t, giay, fi, lamda, m; public void fullt2t() { DateTime d; try { d = new DateTime(nam, thang, ngay, gio, phut, (int)Math.Floor(giay), (int)((giay - Math.Floor(giay)) * 1000)); t = d.Ticks / 10000000.0 / 3600 / 24; } catch { t = -1; } } public void setXy(double lamdamin, double fimin, double dfi) { x = (int)((lamda - lamdamin - (lamda - lamdamin) % (2 * dfi)) / / dfi); y = (int)((fi - fimin - (fi - fimin) % (2 * dfi)) / / dfi); } } Lớp processor public class processor { public processor(Form1 f) { parent = f; } public Form1 parent; public string t2full(double t) { try { double ticks = t * 10000000.0 * 3600 * 24; DateTime d = new DateTime((long)ticks); return d.ToShortDateString(); Trần Việt Phương 66 Luận văn Thạc sĩ } catch { } return "err"; } public trandongdat[][][] tran; public bool isNo(char c) { foreach (int cc in new char[] { '0', '1', '2', '3', '4', '5', '6', '7', '8', '9', '.', '-' }) { if (cc == c) return true; } return false; } //Lệnh createArray, tạo mảng đối tượng trandongdat public trandongdat[] createArray(string text) { trandongdat[] kq; string s = ""; int i, k = -1, j = -1, rows = 0; int length = text.Length; for (i = 0; i < length; i++) { if (text[i].Equals('\n')) { rows++; } } kq = new trandongdat[rows]; for (i = 0; i < length; i++) { if (isNo(text[i])) s += text[i]; else if (s != "") { k++; k = k % 9; try { switch (k) { case 0: j++; kq[j] = new trandongdat(); kq[j].nam = Convert.ToInt32(s); break; case 1: kq[j].thang = Convert.ToInt32(s); break; case 2: kq[j].ngay = Convert.ToInt32(s); break; case 3: kq[j].gio = Convert.ToInt32(s); break; case 4: kq[j].phut = Trần Việt Phương 67 Luận văn Thạc sĩ Convert.ToInt32(s); break; case 5: kq[j].giay = Convert.ToDouble(s); break; case 6: kq[j].fi = Convert.ToDouble(s); break; case 7: kq[j].lamda = Convert.ToDouble(s); break; case 8: kq[j].m = Convert.ToDouble(s); break; } } catch { } s = ""; } } for (i = 0; i < kq.Length; i++) { kq[i].fullt2t(); } return kq; } //Lệnh process dự báo trận động đất (trandongdatdubao[]) public trandongdatdubao[] process(trandongdat[] toanbocactran, double topmax, double dfi, int nmin) { //tao mang dd chua cac tran dong dat co xay nho hon topmax int i = 0; foreach (trandongdat d in toanbocactran) { if ((d.t d.lamda) lamdamin = d.lamda; } // -tach cac tieuvung tran dong dat int xmax = 0, ymax = 0; foreach (trandongdat d in dd) { d.setXy(lamdamin, fimin, dfi); if (xmax < d.x) xmax = d.x; if (ymax < d.y) ymax = d.y; } tran = new trandongdat[xmax + 1][][]; int[][] sotrantrongoluoi = new int[xmax + 1][]; for (int x = 0; x độ lệch label9.Text = "% -> độ lệch "; try { dolech = dolechchuan(Convert.ToDouble(textBox3.Text) / 100); label9.Text += "" + Math.Round(dolech, 4); } catch { dolech = 0; }; } Trần Việt Phương 72 Luận văn Thạc sĩ PHỤ LỤC DANH MỤC CÁC SỰ KIỆN DỰ BÁO ĐÚNG (Trường hợp 2) TT (db) (db) 23.6 24.6 3.6 1.6 0.6 1.6 2.6 -5.4 -4.4 -6.4 10 -5.4 11 -6.4 12 -7.4 13 -6.4 14 -7.4 15 -8.4 16 -10.4 17 -9.4 18 -9.4 19 -8.4 20 -9.4 21 -8.4 22 -7.4 23 -10.4 Trần Việt Phương 94.67 94.67 95.67 97.67 98.67 98.67 98.67 102.67 102.67 103.67 103.67 104.67 105.67 105.67 106.67 112.67 113.67 113.67 114.67 114.67 115.67 115.67 115.67 116.67 T(tựa) 24/12/1999 26/12/1999 7/12/1999 2/12/1999 20/11/1999 22/9/1999 4/6/1999 22/9/1999 6/10/1999 13/10/1999 8/12/1999 9/8/1999 20/5/1999 21/12/1999 4/12/1999 17/11/1999 21/11/1999 18/12/1999 27/11/1999 19/12/1999 30/10/1999 23/12/1999 19/12/1999 25/12/1999 T(dựbáo) 17/2/2000 16/2/2000 25/1/2000 18/1/2000 30/12/1999 19/10/1999 1/7/1999 19/10/1999 20/11/1999 26/11/1999 12/1/2000 19/8/1999 29/6/1999 21/3/2000 22/5/1999 24/12/1999 1/7/2000 16/1/2000 16/1/2000 21/1/2000 19/11/1999 12/1/2000 18/1/2000 18/2/2000 tmin 24/12/1999 26/12/1999 7/12/1999 2/12/1999 20/11/1999 22/9/1999 4/6/1999 22/9/1999 6/10/1999 13/10/1999 8/12/1999 9/8/1999 20/5/1999 21/12/1999 4/12/1999 17/11/1999 21/11/1999 18/12/1999 27/11/1999 19/12/1999 30/10/1999 23/12/1999 19/12/1999 25/12/1999 tmax 14/2/2000 19/4/2000 4/8/2000 14/3/2000 18/2/2000 211/9/1999 8/8/1999 7/12/1999 4/1/2000 20/1/2000 16/2/2000 15/11/1999 23/8/1999 3/2/2000 17/9/1999 8/2/2000 26/2/2000 24/2/2000 4/4/2000 13/3/2000 29/12/1999 12/2/2000 26/2/2000 18/5/2000 73 Mdb Mmin Mmax 3.99 4.12 4.22 3.85 4.61 4.63 3.7 3.91 4.51 4.62 4.63 4.75 4.72 6.15 3.98 4.72 4.43 5.23 5.67 5.04 5.13 4.63 4.53 4.34 3.25 3.7 3.23 3.13 4.01 3.7 2.6 3.16 3.49 3.31 3.78 3.6 3.64 5.09 2.67 3.75 3.52 4.31 4.74 4.09 4.25 3.88 3.89 3.66 4.74 4.54 5.21 4.56 5.21 5.55 4.8 4.65 5.53 5.92 5.48 5.9 5.8 7.22 5.3 5.69 5.34 6.15 6.6 6 5.37 5.16 5.03 TThực 29/3/2000 31/3/2000 8/3/2000 6/1/2000 3/3/2000 14/3/2000 9/5/2000 31/1/2000 2/1/2000 19/2/2000 7/4/2000 6/4/2000 12/1/2000 25/2/2000 15/2/2000 1/7/2000 14/3/2000 16/3/2000 30/1/2000 12/1/2000 1/3/2000 6/1/2000 3/2/2000 9/1/2000 Thực Thực MThực ?T ?M 23.94 24.66 3.64 0.52 1.15 2.18 -5.76 -4.6 -6.72 -5.83 -6.76 -7.31 -6.24 -7.18 -8.32 -9.92 -9.89 -9.55 -8.2 -9.52 -8.32 -7.47 -10.57 94.57 94.46 96.07 98.09 98.55 98.63 99.14 102.33 102.24 104.15 104.01 105.06 105.18 105.17 106.68 113.03 113.42 113.92 114.65 114.75 115.18 116.08 116.05 116.45 4.2 4.1 4.6 4.6 3.3 3.3 4.5 4.4 4.5 3.4 4.2 4.7 4.4 4.6 4.4 5.2 4.4 4.6 5.2 4.9 4.6 4.8 + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + Luận văn Thạc sĩ TT (db) (db) 24 -9.4 25 -8.4 26 -11.4 27 -9.4 28 -8.4 29 -9.4 30 -8.4 31 -8.4 32 -4.4 33 -3.4 34 -2.4 35 -1.4 36 15.6 37 21.6 38 22.6 39 -4.4 40 -2.4 41 -1.4 42 0.6 43 13.6 44 17.6 45 21.6 46 22.6 47 23.6 48 24.6 49 0.6 50 9.6 51 19.6 52 21.6 Trần Việt Phương 116.67 116.67 117.67 117.67 117.67 118.67 118.67 119.67 119.67 119.67 119.67 119.67 119.67 119.67 119.67 120.67 120.67 120.67 120.67 120.67 120.67 120.67 120.67 120.67 120.67 121.67 121.67 121.67 121.67 T(tựa) 25/12/1999 26/12/1999 24/9/1999 24/12/1999 23/12/1999 25/7/1999 9/12/1999 30/12/1999 2/11/1999 5/11/1999 2/11/1999 5/10/1999 12/11/1999 10/11/1999 11/9/1999 28/9/1999 111/8/1999 27/12/1999 26/10/1999 4/12/1999 9/10/1999 1/12/1999 2/10/1999 20/9/1999 22/7/1999 111/9/1999 17/12/1999 3/11/1999 26/12/1999 T(dựbáo) 19/1/2000 22/1/2000 11/8/1999 22/1/2000 20/1/2000 31/8/1999 8/1/2000 4/2/2000 6/12/1999 22/11/1999 14/12/1999 13/6/1999 22/1/2000 9/12/1999 27/12/1999 4/11/1999 27/12/1999 12/1/2000 12/12/1999 26/12/1999 16/11/1999 18/12/1999 13/10/1999 5/10/1999 31/8/1999 9/1/2000 2/2/2000 8/12/1999 5/1/2000 tmin 25/12/1999 26/12/1999 24/9/1999 24/12/1999 23/12/1999 25/7/1999 9/12/1999 30/12/1999 2/11/1999 5/11/1999 2/11/1999 5/10/1999 12/11/1999 10/11/1999 16/11/1999 28/9/1999 111/8/1999 27/12/1999 26/10/1999 4/12/1999 9/10/1999 1/12/1999 2/10/1999 20/9/1999 2/8/1999 24/11/1999 17/12/1999 3/11/1999 26/12/1999 tmax 15/2/2000 29/3/2000 9/1/2000 13/3/2000 23/2/2000 13/10/1999 9/2/2000 21/3/2000 19/1/2000 16/12/1999 30/1/2000 21/7/1999 16/3/2000 11/1/2000 5/2/2000 31/12/1999 15/2/2000 14/3/2000 19/2/2000 22/1/2000 5/1/2000 10/1/2000 1/11/1999 25/10/1999 30/9/1999 24/2/2000 30/4/2000 8/2/2000 19/1/2000 74 Mdb Mmin Mmax 4.5 4.7 4.4 4.61 4.62 4.62 4.52 4.1 4.2 4.62 4.38 3.7 7.44 2.98 3.42 3.27 4.52 4.4 5.36 4.82 3.72 3.8 4.71 7.77 3.62 4.74 3.39 4.82 3.3 3.97 3.9 3.28 3.99 4.01 3.81 3.92 3.47 3.47 3.83 3.6 3.06 6.38 2.41 2.83 2.68 3.81 3.44 4.8 4.05 2.91 3.04 3.99 6.93 3.21 4.19 2.61 3.52 2.27 5.03 5.51 5.51 5.22 5.23 5.43 5.11 4.73 4.93 5.41 5.17 4.34 8.5 3.55 3.86 5.22 5.35 5.91 5.6 4.52 4.57 5.43 8.61 4.03 5.29 4.17 6.12 4.32 TThực 4/1/2000 31/12/1999 1/3/2000 6/1/2000 18/1/2000 28/1/2000 8/2/2000 2/2/2000 35/1/2000 7/2/2000 10/2/2000 16/5/2000 18/1/2001 1/7/2000 1/7/2000 5/4/2000 24/3/2000 22/2/2000 10/4/2000 25/1/2000 3/3/2000 11/7/2000 21/1/2000 22/4/2001 4/5/2001 30/1/2000 116/1/2000 1/3/2000 14/1/2000 Thực Thực MThực ?T ?M -9.2 -8.48 -11.6 -8.92 -8.63 -9.71 -8.74 -8.19 -4.47 -3.18 -2.55 -1.64 15.5 21.37 22.32 -4.71 -2.2 -1.83 0.72 13.62 17.37 22.04 22.2 23.96 24.46 0.31 9.7 19.83 21.33 116.47 116.24 117.75 117.7 117.53 118.54 118.33 119.96 119.61 119.8 120 119.94 119.28 120.12 120 120.2 120.84 120.92 120.94 120.82 120.45 121.07 121.01 121 121.01 121.57 121.84 121.49 121.81 4.8 4.5 4.9 4.5 4.6 5.4 4.6 4.6 4.2 4.5 4.6 3.9 3.9 3.8 3.1 3.1 4.6 4.6 3.9 4.1 3.9 3.1 3.2 3.7 3.2 4.6 4.4 3.8 3.3 + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + Luận văn Thạc sĩ TT (db) (db) 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 22.6 23.6 24.6 -0.4 22.6 23.6 24.6 -8.4 -0.4 6.6 23.6 24.6 -0.4 1.6 5.6 12.6 -0.4 0.6 1.6 2.6 3.6 4.6 5.6 9.6 10.6 11.6 12.6 24.6 -7.4 Trần Việt Phương 121.67 121.67 121.67 122.67 122.67 122.67 122.67 123.67 123.67 123.67 123.67 123.67 124.67 124.67 124.67 124.67 125.67 125.67 125.67 125.67 125.67 125.67 125.67 125.67 125.67 125.67 125.67 125.67 126.67 T(tựa) 13/12/1999 21/7/1999 13/12/1999 10/12/1999 24/10/1999 7/11/1999 23/12/1999 9/12/1999 30/12/1999 5/10/1999 19/12/1999 7/9/1999 13/12/1999 8/10/1999 25/12/1999 30/12/1999 25/9/1999 111/8/1999 25/11/1999 15/11/1999 23/12/1999 8/10/1999 4/12/1999 25/11/1999 31/10/1999 11/8/1999 15/11/1999 23/12/1999 6/8/1999 T(dựbáo) 21/12/1999 26/7/1999 19/12/1999 15/1/2000 22/11/1999 23/11/1999 31/12/1999 26/1/2000 3/2/2000 21/11/1999 20/1/2000 8/10/1999 22/1/2000 17/11/1999 20/2/2000 8/2/2000 5/11/1999 20/12/1999 23/12/1999 23/12/1999 9/2/2000 11/9/1999 27/12/1999 9/1/2000 21/12/1999 12/12/1999 22/12/1999 20/1/2000 19/9/1999 tmin 13/12/1999 21/7/1999 13/12/1999 10/12/1999 24/10/1999 7/11/1999 23/12/1999 9/12/1999 30/12/1999 5/10/1999 19/12/1999 7/9/1999 13/12/1999 8/10/1999 25/12/1999 30/12/1999 25/9/1999 20/11/1999 25/11/1999 15/11/1999 23/12/1999 8/10/1999 4/12/1999 25/11/1999 31/10/1999 12/11/1999 15/11/1999 23/12/1999 6/8/1999 tmax 1/1/2000 3/8/1999 30/12/1999 24/2/2000 8/1/2000 15/12/1999 15/1/2000 17/3/2000 15/3/2000 11/1/2000 22/2/2000 13/11/1999 11/3/2000 27/12/1999 25/3/2000 4/1/2000 4/1/2000 19/1/2000 19/2/2000 18/2/2000 14/2/2000 15/10/1999 22/1/2000 15/3/2000 28/2/2000 10/1/2000 10/2/2000 18/2/2000 4/11/1999 75 Mdb Mmin Mmax 3.7 3.3 3.5 4.11 3.27 3.4 3.59 4.53 3.64 4.43 3.9 4.43 5.47 4.43 4.08 3.99 4.08 5.75 4.11 4.75 3.87 4.31 4.21 3.65 4.62 3.97 4.41 3.21 3.8 2.58 2.77 3.2 2.47 2.44 2.47 3.46 2.27 3.39 2.97 3.55 4.19 3.78 2.88 3.38 4.58 3.25 3.89 3.37 3.75 3.56 2.49 3.68 3.47 3.77 2.44 2.87 4.4 4.02 4.23 5.02 4.07 4.37 4.71 5.59 5.01 5.48 4.84 5.31 6.74 5.09 5.17 5.1 4.78 6.92 4.96 5.61 4.37 4.88 4.85 4.81 5.56 4.48 5.06 3.97 4.74 TThực 5/3/2000 19/4/2001 3/2/2000 31/3/2000 2/1/2000 11/1/2000 31/12/1999 7/6/2000 1/7/2000 8/1/2000 21/1/2000 27/2/2000 5/3/2000 29/1/2000 2/1/2000 20/7/2000 17/3/2000 24/1/2000 19/6/2000 2/1/2000 4/4/2000 15/1/2000 11/7/2000 5/2/2000 30/1/2000 24/1/2000 26/3/2000 15/3/2000 15/5/2000 Thực Thực MThực ?T ?M 22.37 23.73 24.86 -0.11 22.28 23.4 24.8 -8.53 -0.12 6.9 23.38 24.81 -0.19 1.92 5.65 12.42 0.05 0.98 1.47 2.74 4.09 4.87 5.96 9.69 10.42 12.03 12.55 24.17 -6.93 121.82 121.63 122.14 123.06 122.66 122.99 122.25 123.96 123.29 123.76 123.76 123.18 124.38 124.3 125.16 125.14 125.95 125.35 125.66 126.01 125.59 125.54 125.83 125.29 126.14 125.95 125.87 125.2 127.16 4.3 4.2 3.6 3.2 3.5 3.5 3.3 3.4 3.1 4.4 3.8 4.8 3.9 3.7 4.3 3.8 5.7 4.1 4.3 3.5 4.4 4.2 3.8 3.6 4.6 + + + + + + + + + + + - + + + + + + + + + + + + + + + + + + + + + + Luận văn Thạc sĩ TT (db) (db) 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 -0.4 0.6 1.6 2.6 3.6 4.6 5.6 6.6 7.6 8.6 9.6 10.6 -7.4 -6.4 -4.4 -3.4 -2.4 -1.4 1.6 2.6 3.6 4.6 5.6 -8.4 -7.4 -6.4 -5.4 -3.4 -2.4 Trần Việt Phương 126.67 126.67 126.67 126.67 126.67 126.67 126.67 126.67 126.67 126.67 126.67 126.67 127.67 127.67 127.67 127.67 127.67 127.67 127.67 127.67 127.67 127.67 127.67 128.67 128.67 128.67 128.67 128.67 128.67 T(tựa) 7/7/1999 23/10/1999 14/12/1999 1/12/1999 26/12/1999 1/12/1999 18/12/1999 18/12/1999 31/10/1999 14/12/1999 10/12/1999 3/12/1999 4/12/1999 21/11/1999 9/12/1999 19/10/1999 31/5/1999 27/12/1999 15/11/1999 4/10/1999 4/12/1999 27/11/1999 28/12/1999 10/11/1999 24/12/1999 17/11/1999 5/9/1999 16/11/1999 12/12/1999 T(dựbáo) 17/8/1999 25/11/1999 10/1/2000 23/12/1999 14/1/2000 2/1/2000 12/1/2000 14/1/2000 27/11/1999 8/1/2000 30/12/1999 8/1/2000 2/1/2000 9/1/2000 19/1/2000 21/11/1999 5/7/1999 6/2/2000 6/12/1999 29/10/1999 1/7/2000 5/1/2000 12/1/2000 21/12/1999 18/1/2000 28/12/1999 23/10/1999 21/12/1999 16/1/2000 tmin 7/7/1999 26/10/1999 14/12/1999 1/12/1999 26/12/1999 1/12/1999 18/12/1999 18/12/1999 31/10/1999 14/12/1999 10/12/1999 3/12/1999 4/12/1999 21/11/1999 9/12/1999 19/10/1999 31/5/1999 27/12/1999 15/11/1999 4/10/1999 4/12/1999 27/11/1999 28/12/1999 10/11/1999 24/12/1999 17/11/1999 5/9/1999 16/11/1999 12/12/1999 tmax 4/10/1999 25/12/1999 3/3/2000 18/1/2000 7/2/2000 14/2/2000 1/3/2000 14/2/2000 31/12/1999 3/2/2000 31/1/2000 4/3/2000 5/2/2000 3/3/2000 5/3/2000 8/1/2000 8/10/1999 5/4/2000 2/1/2000 7/12/1999 20/2/2000 21/2/2000 17/4/2000 7/2/2000 31/3/2000 16/4/2000 22/12/1999 21/7/2000 1/3/2000 76 Mdb Mmin Mmax 3.78 4.5 4.51 4.82 4.3 3.8 4.49 3.07 4.51 4.12 4.1 4.3 4.73 3.78 4.14 4.3 4.02 4.74 3.6 5.44 5.86 4.11 4.1 3.15 3.98 4.01 4.54 3.79 5.44 3.08 3.65 3.46 3.84 3.35 3.06 3.47 2.05 3.82 3.22 3.39 3.57 3.94 3.13 3.24 3.54 3.23 3.82 2.62 4.75 3.3 3.34 2.33 2.79 2.95 3.3 3.26 4.78 4.48 5.35 5.56 5.8 5.26 4.53 5.51 4.09 5.19 5.02 4.81 5.03 5.51 4.42 5.03 5.06 4.81 5.66 4.57 6.14 6.72 4.91 4.86 3.97 5.16 5.07 5.77 4.31 6.09 TThực 4/8/2000 27/5/2000 18/1/2000 20/3/2000 31/1/2000 11/3/2000 25/1/2000 14/1/2000 4/1/2000 1/1/2000 1/3/2000 31/1/2000 28/3/2000 2/2/2000 18/1/2000 13/2/2000 25/1/2000 22/1/2000 15/2/2000 1/1/2000 23/2/2000 9/1/2000 17/8/2000 12/7/2000 25/1/2000 6/1/2000 30/5/2000 4/1/2000 10/7/2000 Thực Thực MThực ?T ?M 0.05 0.77 1.25 2.87 3.46 4.24 5.93 6.52 8.07 8.35 9.63 10.13 -7.63 -6.37 -4.16 -3.59 -2.67 -1.08 1.5 2.41 3.15 5.07 5.67 -7.94 -6.92 -5.99 -5.7 -3.16 -2.67 126.55 126.37 126.19 126.59 126.46 126.57 126.98 126.97 126.87 126.27 126.22 126.41 127.99 128.01 127.96 127.19 127.4 127.31 127.42 127.91 128.03 128.05 127.26 128.26 128.65 128.59 128.96 128.53 128.34 5.3 4.9 4.6 4.3 4.4 4.5 4.2 3.9 4.5 4.9 3.9 3.9 4.2 4.4 4.3 4.4 4.1 4.2 6.1 3.7 3.6 4.3 3.7 4.3 4.2 4.6 + + + + + + + + + + + + + + + - + + + + + + + + + + + + + + + + + + + + + - Luận văn Thạc sĩ TT (db) (db) 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 1.6 2.6 3.6 -8.4 -7.4 -6.4 -5.4 -4.4 -3.4 -7.4 -6.4 -5.4 -3.4 -3.4 -2.4 -1.4 -3.4 -2.4 -3.4 -2.4 Trần Việt Phương 128.67 128.67 128.67 129.67 129.67 129.67 129.67 129.67 129.67 130.67 130.67 130.67 130.67 134.67 134.67 136.67 138.67 138.67 139.67 139.67 T(tựa) 28/12/1999 26/12/1999 25/10/1999 21/7/1999 10/12/1999 26/12/1999 7/10/1999 1/12/1999 9/9/1999 18/5/1999 28/12/1999 27/12/1999 13/11/1999 19/10/1999 28/10/1999 3/12/1999 7/7/1999 11/8/1999 25/7/1999 18/12/1999 T(dựbáo) 10/2/2000 15/1/2000 211/9/1999 25/8/1999 5/1/2000 21/7/2000 22/11/1999 11/7/2000 8/10/1999 3/7/1999 11/7/2000 8/2/2000 1/7/2000 24/11/1999 3/12/1999 6/1/2000 19/8/1999 7/9/1999 28/8/1999 3/2/2000 tmin 28/12/1999 26/12/1999 25/10/1999 21/7/1999 10/12/1999 26/12/1999 7/10/1999 1/12/1999 9/9/1999 18/5/1999 28/12/1999 27/12/1999 13/11/1999 19/10/1999 1/11/1999 3/12/1999 7/7/1999 11/8/1999 25/7/1999 18/12/1999 tmax 14/1/2000 12/1/2000 9/1/2000 6/10/1999 17/2/2000 21/3/2000 24/1/2000 5/4/2000 25/11/1999 30/9/1999 12/2/2000 18/4/2000 4/3/2000 9/1/2000 4/1/2000 23/2/2000 8/10/1999 26/10/1999 24/10/1999 10/4/2000 77 Mdb Mmin Mmax 5.68 4.01 3.9 4.23 4.11 4.51 4.31 4.84 4.33 3.39 4.71 4.28 4.24 4.09 4.73 4.21 3.99 5.55 5.37 6.31 5.11 3.02 3.36 3.7 3.22 3.47 3.53 3.84 3.68 2.72 3.83 3.6 3.28 3.61 3.86 3.42 3.2 4.29 4.44 4.84 6.26 4.99 4.44 4.76 5.54 5.1 5.84 4.99 4.06 5.6 4.95 5.21 4.57 5.59 4.99 4.78 6.81 6.3 7.78 TThực 5/10/2000 1/3/2000 4/4/2000 13/6/2000 15/1/2000 22/1/2000 14/3/2000 9/9/2000 1/7/2000 21/2/2000 11/3/2000 14/3/2000 1/7/2000 1/3/2000 24/8/2000 5/5/2000 1/7/2000 10/1/2000 31/12/1999 22/2/2000 Thực Thực MThực ?T ?M 1.73 2.15 3.98 -8.78 -7.26 -6.59 -4.93 -4.79 -3.3 -6.95 -6.37 -5.53 -2.95 -3.73 -2.44 -1.51 -3.27 -2.18 -3.88 -2.01 128.44 128.8 128.35 129.43 129.43 129.78 130.1 129.45 129.54 130.54 131 131.11 130.25 134.48 134.32 136.74 138.5 139.16 139.58 139.29 4.2 4.2 4.7 4.4 3.9 4.3 5.1 4.3 3.9 4.7 4.9 3.4 4.5 3.8 4.4 4.2 3.9 4.6 5.1 4.7 + + + + + + + + + + + + + + + + + + + + + + -