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ê (luan van thac si)

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ê (luan van thac si) 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ê (luan van thac si) 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ê (luan van thac si) 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ê (luan van thac si) 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ê (luan van thac si) 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ê (luan van thac si) 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ê (luan van thac si) 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ê (luan van thac si)

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 TÀI LIỆU THAM KHẢO 58 MỞ ĐẦU Những năm gần đây, thảm hoạ động đất, sóng thần toà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 mang tính thời qui mô toà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à 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 đấ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 3 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 toán sơ đồ khối chương trình dự báo động đất theo mô hình thống kê” Trong phầ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 toá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 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 xuyên xảy động đất tính toá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 toá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ố loà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 hoà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 đứt đoạn nằm sâu bên thạch Gần đây, chuyên gia khoa học Viện nghiên cứu Carnegie - Mỹ tìm cách để kiểm tra giám sát chiều dài đoạn đứt gãy, nhƣ dịch chuyển chúng vỏ Trái đất Phát phƣơng pháp đầy hữu ích, giúp cho việc dự báo trận động đất cách định vị xác đứt gãy có khả làm rung chuyển mặt đất gây trận động đất Trong đó, chuyên gia Viện Nghiên cứu vũ trụ thuộc Viện Hàn lâm khoa học Nga lại tiếp cận việc dự báo động đất sớm từ vũ trụ việc lắp đặt thiết bị dự báo động đất trạm không gian quốc tế Thiết bị ghi nhận biến đổi dòng điện tử proton có lƣợng trung bình khoảng không gian gần Trái đất Các nhà khoa học cho thay đổi dòng điện tử proton có liên quan tới trình địa vật lý Trái đất nhƣ dông tố, động đất, vv… Do vậy, ghi nhận đƣợc thay đổi này, dự báo động đất với độ xác cao Ngoài ra, ngƣời ta sử dụng hàng loạt phƣơng pháp nghiên cứu dự báo động đất khác dựa việc làm sáng tỏ dấu hiệu khác nhƣ: - Các dấu hiệu phát triển vùng nguồn (dị thƣờng mật độ đứt gãy sinh chấn), - Quan hệ vận tốc lan truyền sóng dọc sóng ngang, - Năng lƣợng giải phóng chuỗi động đất, hàng loạt dấu hiệu khác có đặc trƣng báo trƣớc động đất tƣơng lai (các động đất xảy ra) Trong loạt phƣơng pháp khác ngƣời ta lại sử dụng bất thƣờng đồng thời diễn biến số dấu hiệu có chất vật lý khác Một số công trình theo hƣớng chứa công thức thuật toán dự báo động đất Đó thuật toán dự báo động đất CN M8 Kết dự báo động đất theo kiểu nhƣ vùng động đất mạnh tƣơng lai, mà khoảng thời gian (thƣờng nhiều năm gọi dự báo trung hạn) xảy động đất Ƣu điểm hệ phƣơng pháp thuật toán dự báo trung hạn nhƣ làm sáng tỏ đƣợc dấu hiệu động đất tƣơng lai, có ý nghĩa vật lý cho phép đánh giá chúng phƣơng pháp thống kê Tuy nhiên, tính không xác định độ xác dự báo động đất trung hạn coi nhƣợc điểm phƣơng pháp Vì vậy, luận văn này, tác giả ứng dụng mô hình thống kê để thiết lập chƣơng trình thử nghiệm dự báo ngắn hạn động đất magnitude chúng, áp dụng lãnh thổ Việt Nam vùng lân cận Cách tiếp cận vừa đơn giản mà lại có khả loại bỏ đƣợc nhƣợc điểm phƣơng pháp nêu 1.2 Tình hình nghiên cứu dự báo động đất Việt Nam: Nghiên cứu dự báo động đất Việt Nam chủ yếu đƣợc tiến hành theo hai nhóm phƣơng pháp nhƣ sau:  Dự báo động đất cực đại sở tài liệu địa chất – địa vật lý  Dự báo động đất cực đại sở phân tích thống kê địa chấn * Dự báo động đất cực đại sở tài liệu địa chất – địa vật lý bao gồm số phương pháp sau: 1.2.1 Phương pháp tính Mmax theo quy mô vùng phát sinh động đất: Theo phƣơng pháp này, ngƣời ta dựa vào liên hệ kích thƣớc đoạn đứt gãy sinh chấn (L) nhƣ bề dày tầng sinh chấn (H) với M max động đất [37, 38, 43, 44]: Mmax ≤ 2lg L(km) + 1.77 Mmax ≤ 4lg H(km) + 0.48 Độ xác phƣơng pháp phụ thuộc vào việc xác định kích thƣớc đoạn đứt gãy sinh chấn bề dày tầng sinh chấn nhƣ Phƣơng pháp áp dụng tốt cho vùng có đứt gãy hoạt động nhƣng đồng thời lại dự báo đƣợc Mmax cho vùng khác đứt gãy hoạt động 1.2.2 Phương pháp kiến tạo vật lý: Theo phƣơng pháp cho rằng, vị trí, lƣợng động đất mạnh cực đại vùng phụ thuộc vào gradient vận tốc chuyển động kiến tạo thẳng đứng thời kì tân kiến tạo Trong gradient vận tốc lớn 10-8/ năm diện tích 1000 km2 xảy động đất cấp lần 1000 năm Sau gradient vận tốc tăng lên lần độ mạnh động đất tăng lên 1-2 cấp Song cần lƣu ý rằng, động đất trƣớc hết hệ chuyển động đại mà gradient vận tốc trung bình chu kỳ tân kiến tạo không phản ánh đƣợc, gradient vận tốc không phản ánh mức độ nguy hiểm động đất Thêm động đất mạnh không hệ chuyển động thẳng đứng mà thƣờng hệ chuyển động ngang Trong trƣờng hợp gradient vận tốc chuyển động thẳng đứng không nói lên điều quan trọng 1.2.3 Phương pháp đánh giá chuyên gia Cơ sở phƣơng pháp việc xác định mức độ phản ánh lớp đối tƣợng địa chất – địa vật lý đặc trƣng hoạt động động đất Mức độ phản ánh đƣợc biểu diễn qua trọng số tham số Sau phân loại đối tƣợng, gán trọng số cho đối tƣợng, tổng hợp ảnh hƣởng đối tƣợng để dự báo khu vực có khả phát sinh động đất mạnh nhƣ thành lập đồ dự báo chấn động cực đại [42] 1.2.4 Sử dụng kết hợp tài liệu địa chất – địa vật lý: Các phƣơng pháp bƣớc đầu đƣợc áp dụng Việt Nam số công trình sở sử dụng đặc trƣng cấu trúc vỏ [8, 39, 40] Ngoài ra, tác giả công trình [27] sử dụng số đặc trƣng cấu trúc vỏ Trái đất để áp dụng thử nghiệm dự báo động đất theo vecto dấu hiệu nhiều chiều Tuy nhiên, kết nhận đựơc công trình nêu sơ lƣợc Mặt khác chƣa có tác giả xây dựng đƣợc thuật toán thiết lập đƣợc chƣơng trình dự báo động đất phù hợp với điều kiện thực tế Việt Nam sở phƣơng pháp áp dụng * Các phương pháp dự báo động đất cực đại sở phân tích thống kê địa chấn sử dụng Việt Nam chủ yếu là: Bài toán phân bố cực trị Gumbel: Cơ sở lý thuyết hàm phân bố cực trị Gumbel đƣợc trình bày nhiều công trình nghiên cứu Việt Nam giới [52, 53, 10, 27, 28, 34] Nếu ta coi X biến ngẫu nhiên có hàm phân bố F(X) F(X) = P{X  x} xác suất x lớn n mẫu độc lập từ phân bố F(X) là: G(x) = P{ X1  x, X2  x, , Xn  x } = Fn (x) Đó hàm phân bố cực trị Nếu nhƣ ta biết đƣợc hàm phân bố ban đầu F(X) đơn giản để nhận đƣợc xác phân bố cực trị Nhƣng thông thƣờng ta đƣợc hàm phân bố ban đầu, nên cần phải xem đến dạng đƣờng tiệm cận phân bố cực trị Khi áp dụng lý thuyết phân bố cực trị Gumbel vào dự báo động đất phải tuân thủ giả thiết, [10]: - Các cực trị quan sát đƣợc khoảng thời gian cho trƣớc độc lập - Các điều kiện xảy khứ xảy tƣơng lai Gumbel xây dựng đƣợc loại hàm phân bố tiệm cận cực trị, hàm phân bố cực trị loại I có dạng: với  > G1 ( x)  exp[ e   ( xu ) ] đây: u - đặc trƣng cực trị, - hàm cƣờng độ cực trị ; u  tham số cần xác định Phân bố loại II có dạng: G ( x)  exp[ ( u  k ) ] với k > 0, x ≥ ε , u > ε ≥ x  Trong k – đại lƣợng hình dạng, ε – giá trị cận dƣới cực trị, u – đặc trƣng cực trị, u, k, ε tham số cần xác định Phân bố loại III có dạng: G ( x)  exp[ ( w x k ) ] với k > 0, x  w, u < w wu Trong đó, w – giá trị cận cực trị, u, k, w tham số cần xác định Từ hàm phân bố tiệm cận ta thấy phân bố tiệm cận loại II có tồn giá trị cận dƣới ε nên chúng không đƣợc sử dụng để đánh giá chấn cấp cực đại động đất Thông thƣờng dùng hàm lại để giải nhiệm vụ Nguyễn Kim Lạp Nguyễn Duy Nuôi [19] sử dụng hàm phân bố tiệm cận loại I Gumbel để tính độ nguy hiểm động đất cho vùng khu vực Đông Nam Á với chu kỳ khoảng số liệu cực trị tháng năm Nguyễn Hồng Phƣơng [28] sử dụng hàm phân bố loại III Gumbel kết hợp với hợp lý cực đại phân bố β để tính cho vùng nhỏ lãnh thổ Việt Nam Nguyễn Đình Xuyên Nguyễn Ngọc Thủy sử dụng hàm Gumbel loại I cải tiến để tính cho vùng thuộc khu vực thủy điện Sơn La [43] Trần Thị Mỹ Thành tính Mmax cho 23 vùng nguồn dựa theo hàm Gumbel III [34]… Phương pháp hợp lý cực đại tính đƣợc giới hạn chặn hai phía dãy trận động đất hoạt động địa chấn theo luật phân bố Poisson biểu thức Gutenberg-Richter Phƣơng trình biểu diễn quan hệ tần suất xuất động đất NM chấn cấp M phƣơng trình tiếng Gutenberg-Richter [53]: lg NM = a - bM Quy luật xuất động đất tuân theo quy luật phân bố Poisson Trong vùng nguồn, coi động đất kiện độc lập (loại bỏ tiền chấn dƣ chấn), xác suất P N để xảy N trận động đất có chấn cấp M ≥ m0, gây cƣờng độ chấn động I lớn mức i đó, toàn vùng nguồn khoảng thời gian t năm thỏa mãn phƣơng trình: e t (t ) n PN= P [ N = n ] = n! n = 0, 1,  vận tốc trung bình xuất động đất có chấn cấp M ≥ m0 Nếu xem pi xác suất xuất trận động đất có chấn cấp M ≥ m0, xác suất để xuất N trận động đất là: 10 Luận văn Thạc sĩ Bảng 4.1: Kết áp dụng chương trình thử nghiệm dự báo động đất cho khu vực Tây B c Việt Nam trường hợp TT (db) (db) Ttựa Tdb Tmin Tmax Mdb Mmin Mmax Tthực Thực Thực MThực ?T ?M Bảng 4.1.1: Trường hợp 1 20.98 102.44 19/8/2009 21.98 102.44 23/11/2011 12/6/2011 19/8/2009 6/5/2012 23/11/2011 5/10/2014 3.93 2.8 5.05 0 - - 2/3/2013 2.18 1.26 3.11 0 - - 22.98 102.44 13/2/2011 16/6/2012 13/2/2011 17/2/2014 2.86 1.73 3.99 0 - - 20.98 103.44 29/6/2011 18/12/2011 29/6/2011 8/12/2012 3.61 2.52 4.7 0 - - 21.98 103.44 4/11/2010 17/3/2011 4/11/2010 2/10/2011 2.48 1.44 3.52 0 - - 22.98 103.44 21/5/2011 13/6/2012 21/5/2011 8/8/2014 2.13 1.18 3.09 0 - - 20.98 104.44 23/2/2010 26/10/2010 23/2/2010 4/12/2011 2.5 1.62 3.38 0 - - 21.98 104.44 26/2/2006 30/4/2008 26/2/2006 15/6/2012 3.41 1.98 4.85 0 - - 20.98 105.44 20/6/2005 30/8/2006 20/6/2005 17/2/2009 2.44 1.49 3.38 0 - - Bảng 4.1.2: Trường hợp 21.98 102.44 28/11/2001 22.98 102.44 1/11/2000 18/5/2002 28/11/2001 3.9 3.07 4.74 15/3/2002 21.84 102.305 3.01 + - 1/11/2000 28/7/2004 3.36 2.36 4.35 24/2/2002 22.62 102.332 3.01 + + 20.98 103.44 13/11/2001 21.98 103.44 28/11/2001 20/6/2002 13/11/2001 4/7/2003 3.01 1.87 4.15 22/1/2002 20.79 103.69 3.9 + + 3/4/2002 28/11/2001 3/10/2002 3.79 2.71 4.88 14/4/2002 21.766 103.59 2.3 + - 22.98 103.44 15/8/2001 20.98 104.44 12/11/2001 26/2/2003 20.98 105.44 20/11/2001 6/9/2002 16/4/2003 15/8/2001 10/11/2005 2.93 1.94 3.91 26/7/2002 12/11/2001 5/4/2003 20/11/2001 3.3 + + 23/6/2002 21.13 104.525 3.01 + + 1.9 4.01 20/12/2001 20.8 105.326 3.01 + + 14/9/2003 3.53 2.66 4.39 8/1/2006 2.96 2/7/2002 22.796 103.639 Chú thích cho bảng 4.1: (db), (db) - tọa độ kiện dự báo; tmin - thời gian xảy sớm nhất; tmax - thời gian xảy muộn nhất; Mdb - magnitude kiện dự báo; Mmin - magnitude nhỏ kiện dự báo; Mmax - magnitude lớn kiện dự báo; Tthực, Thực, Thực, M thực- tương ứng thời gian, tọa độ, magnitude kiện xảy thực tế; ?T - thời gian dự báo có hay không; ?M - magnitude dự báo có hay không; “ +” : đúng; “ - “: không Trần Việt Phương 43 Luận văn Thạc sĩ Thử kiểm tra thời gian dự báo muộn (Tmax) so với thời gian dự báo (Tdb) (cả đại lƣợng phần mềm tính ra) thấy độ chênh lệch chúng tƣơng đối lớn, từ 183 ngày 1009 ngày Đồng thời, tác giả kiểm tra độ chênh lệch thời gian thực xảy động đất (Tthực) với thời gian dự báo (Tdb) thấy độ chênh lệch dao động từ -11 ngày đến 471 ngày Điều phản ánh thực tế mức độ tin cậy phƣơng pháp dự báo động đất theo mô hình thống kê Rõ ràng trƣờng hợp thử lùi danh mục khứ kiểm tra đƣợc mức độ tin cậy dự báo thông qua hiệu Tthực-Tdb Nhƣng áp dụng cho thời điểm tựa thời điểm tại, liệu để kiểm tra mức độ tin cậy, trận động đất dự báo nằm tƣơng lai Tuy nhiên, tác giả nhận thấy hiệu Tmax-Tdb vừa thử hiệu TthựcTdb có dấu hiệu tỷ lệ thuận với (xem hình 4.1) Bảng 4.2: So sánh độ chênh lệch Tmax-Tdb v i Tthực- Tdb Tdb stt 18/5/2002 6/9/2002 20/6/2002 3/4/2002 26/2/2003 26/7/2002 5/4/2003 Tmax 16/4/2003 28/7/2004 4/7/2003 3/10/2002 10/11/2005 14/9/2003 8/1/2006 Tmax-Tdb 333 691 379 183 988 415 1009 571.14 Tthực 15/3/2002 24/2/2002 22/1/2002 14/4/2002 2/7/2002 23/6/2002 20/12/2001 Tthực- Tdb 64.00 194.00 149.00 -11.00 239.00 33.00 471.00 162.71 Tthực- Tdb Tmax-Tdb Hình 4.1: So sánh độ chênh lệch thời gian Tmax-Tdb Tthực- Tdb Trần Việt Phương 44 Luận văn Thạc sĩ Thông qua quan hệ nhƣ vậy, dựa vào độ chênh lệch Tmax-Tdb để ƣớc lƣợng giá trị độ lệch Tthực- Tdb, từ cho phép dự đoán trƣớc mức độ xác phép dự báo chƣa có số liệu thực để đối chiếu Ngoài ra, thử nghiệm dự báo trƣờng hợp sử dụng danh mục động đất từ năm 1903-2011 thu đƣợc kết bảng 4.3 Bảng 4.3: Ví dụ kết dự báo đối v i số trận động đất TT (db) (db) Ttựa Tdb Tmin Tmax Mdb Mmin Mmax Tthực Thực Tthực - Tdb 21.98 22.98 20.98 21.98 22.98 20.98 21.98 20.98 102.44 102.44 103.44 103.44 103.44 104.44 104.44 105.44 11/28/2001 11/1/2000 11/13/2001 11/28/2001 8/15/2001 11/12/2001 4/20/2000 11/20/2001 7/31/2002 2/16/2003 5/21/2003 11/19/2002 4/21/2003 2/28/2004 5/30/2004 6/1/2004 11/28/2001 11/1/2000 11/13/2001 11/28/2001 8/15/2001 11/12/2001 4/20/2000 11/20/2001 6/24/2004 11/16/2005 4/28/2009 8/31/2005 2/9/2006 5/1/2010 10/28/2011 5/24/2010 3.87 3.36 2.96 3.78 2.95 3.47 3.86 2.96 2.97 2.27 1.29 2.67 2.61 2.56 2.01 4.77 4.45 4.64 4.9 3.89 4.33 5.16 3.92 3/15/2002 2/24/2002 1/22/2002 4/14/2002 7/2/2002 6/23/2002 8/1/2002 12/20/2001 21.84 22.62 20.79 21.766 22.796 21.13 21.94 20.8 -138 -357 -484 -219 -293 -615 -668 -894 Thực 102.305 102.332 103.69 103.59 103.639 104.525 104.023 105.326 MThực 3.01 3.01 3.9 2.3 3.3 3.01 2.6 3.01 ?T ?M + + + + + + + + + + + + + + + Từ bảng 4.3 thấy, độ chênh lệch thời gian bị âm (dự báo bị muộn) trị tuyệt đối lớn hẳn so với dự báo với danh mục tính từ năm 1970 Khi thu hẹp danh mục lại (1970 đến 2011), nhận thấy sai số có số âm số dƣơng (xem bảng 4.2), đồng thời sai số nhỏ trƣớc Độ chênh lệch thời gian đƣợc tính thời gian thực diễn trận động đất, trừ thời điểm dự báo động đất Nguyên nhân điều trƣớc năm 1970, số liệu động đất đƣợc ghi lại dựa vào ghi chép đơn truyền miệng, hệ thống đài trạm chƣa đầy đủ nên liệu bị thƣa nhiều Nhƣ vậy, tính đầy đủ tính liên tục số liệu động đất theo thời gian yếu tố quan trọng định độ xác phƣơng pháp dự báo theo mô hình thống kê Trần Việt Phương 45 Luận văn Thạc sĩ 4.3 Áp dụng thử nghiệm dự báo danh mục động đất Đông Nam Á Để kiểm nghiệm khả áp dụng chƣơng trình, tác giả áp dụng thử nghiệm chƣơng trình danh mục động đất Đông Nam Á Trƣớc xử lý số liệu, tác giả có khảo sát qua mức độ liên tục liệu, nhận thấy sau năm 2001 có bất thƣờng khoảng cách thời gian trận động đất liền (xem hình 4.2) 25 ∆T(ngày) 20 15 10 2006 2005 2001 2001 2001 2001 2000 2000 2000 2000 1999 1999 1999 1998 1998 1998 1998 1998 1997 1997 1997 1997 1997 Hình 4.2 : Biểu đồ thể tính liên tục dải số liệu thời gian xuất động đất Chú thích: ∆T - số ngày trận động đất liền Từ hình 4.2 thấy, từ năm 2002 đến năm 2007, khoảng cách thời gian trận động đất liền lớn hẳn so với trƣớc đó, nghĩa số lƣợng trận động đất danh mục rời rạc nhiều so với trƣớc năm 2002 Qua nhận xét sơ danh mục động đất nhƣ vậy, tác giả định sử dụng danh mục động giai đoạn 1997 - 2001 để áp dụng thử nghiệm chƣơng trình nhằm so sánh với kết áp dụng chƣơng trình sử dụng toàn danh mục (đến hết 2007) để kiểm tra ảnh hƣởng tính liên tục tính đầy đủ liệu đầu vào kết chƣơng trình Trần Việt Phương 46 năm Danh mục động đất Đông Nam Á bao gồm 37194 động đất giai đoạn từ 1/1/1997 đến 30/12/2007 Sau tách nhóm tiền chấn dƣ chấn, danh mục lại 15701 trận động đất độc lập Áp dụng chƣơng trình danh mục 15701 trận động đất độc lập với thông tin thiết đặt ban đầu nhƣ sau : - Độ xác dự báo (hay xác suất dự báo) Pg = 80%.; - Kích thƣớc tiểu vùng S: xλ = 0,5x0,5độ; Số kiện tối thiểu chứa tiểu vùng S 20 (trận) - Trƣờng hợp 1: Khi sử dụng danh mục động đất đến ngày 31/12/2005, chƣơng trình dự báo đƣợc 259 trận động đất, có 32 trận đƣợc dự báo (nghĩa thực tế xảy trận động đất nhƣ khoảng thời gian dự báo nói trên) - Trƣờng hợp 2: Khi lùi danh mục khứ trƣớc năm 2002, (giả sử lùi ngày 31/12/1999) chƣơng trình dự báo đƣợc 130 trận động đất, có 56 trận đƣợc dự báo khoảng thời gian, 99 trận magnitude, 45 trận thời gian lẫn magnitude động đất Để phù hợp với khuôn khổ luận văn, tác giả xin dẫn danh mục 130 trận động đất dự báo trƣờng hợp (Phụ lục 2) So sánh độ lệch thời gian xảy trận động đất dự báo trƣờng hợp với thời gian xảy kiện thực tế đƣợc trình bày hình 4.3 47 Hình 4.3: Phân bố độ chênh lệch thời gian dự báo thời gian đ xảy kiện thực tế (ngày) Chú thích: hình trên, trục tung số kiện dự báo, giá trị trục hoành độ chênh lệch tương ứng Ví dụ cột thứ hiểu : có kiện dự báo có độ chênh lệch khoảng từ -50 đến -30, cột thứ : có 16 kiện dự báo có độ chênh lệch từ -30 đến -10 Từ hình 4.3 ta thấy hầu hết độ chênh lệch số dƣơng, tức Tthực > Tdb Các kết hình 4.3 cho thấy độ chênh lệch trung bình thời gian dự báo thời gian xảy kiện thực tế 95.6 ngày Qua thấy kiện dự báo hầu hết có xu sớm so với thời gian xảy kiện thực tế Khoảng thời gian sớm trung bình 95 ngày Thời gian đủ để ngƣời chuẩn bị thực biện pháp sơ tán phòng tránh động đất xảy 48 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 Với mục tiêu so sánh kết chƣơng trình dự báo động đất với độ tin cậy khác nhau, chƣơng trình đƣợc áp dụng để dự báo động đất trƣờng hợp lựa chọn xác suất dự báo khác nhau: Pg = 70% Pg = 90% Các kết đƣợc trình bày bảng 4.4 Kết bảng 4.4 cho thấy việc thử nghiệm dự báo động đất với độ tin cậy (xác suất dự báo Pg) khác cho kết khác khoảng cách thời điểm sớm thời điểm muộn xảy kiện dự báo nhƣ độ chênh lệch magnitude cực đại magnitude cực tiểu kiện dự báo So sánh kết bảng 4.4.1 với kết bảng 4.4.2 cách trừ tất giá trị (phần tử) tƣơng ứng cột hàng bảng 4.4.1 cho bảng 4.4.2 Hiệu chúng đƣợc dẫn bảng 4.4.3 Các kết bảng 4.4.3 cho thấy độ chênh lệch thời điểm muộn (tmax) xảy kiện dự báo trƣờng hợp dao động từ 87 ngày đến 478 ngày Thời điểm sớm xảy kiện dự báo trƣờng hợp khác kiện dự báo thứ (116 ngày) Còn độ chênh lệch magnitude dự báo trƣờng hợp dao động từ 0.4 đến 0.54 đối magnitude cực tiểu từ (-0.4) đến (-0.54) magnitude dự báo cực đại Rõ ràng lựa chọn xác xuất tin cậy dự báo thấp (Pg = 70 %) chƣơng trình cho kết dự báo với magnitude cực đại, khoảng thời gian dự báo khoảng cách magnituge cực đại cực tiểu nhỏ so với trƣờng hợp lựa chọn xác xuất tin cậy dự báo cao (Pg = 90 %) 49 Bảng 4.4 So sánh kết dự báo động đất v i xác xuất tin cậy khác TT (db) (db) fiOp 21.98 22.98 20.98 21.98 22.98 20.98 20.98 21.98 22.98 20.98 21.98 22.98 20.98 20.98 0 0 0 Ttựa Tdb Tmin Tmax Mdb Mmin Mmax Tthực Thực Thực MThực ?T ?M Bảng 4.4.1 Trường hợp lựa chọn xác suất dự báo Pg = 70% 102.44 22.47 102.6 28/11/2001 18/5/2002 28/11/2001 11/2/2003 3.9 3.23 4.58 15/3/2002 21.84 102.305 3.01 + 102.44 22.81 102 1/11/2000 6/9/2002 25/2/2001 18/3/2004 3.36 2.55 4.16 24/2/2002 22.62 102.332 3.01 + + 103.44 21.39 103.046 13/11/2001 20/6/2002 13/11/2001 23/4/2003 3.01 2.09 3.93 22/1/2002 20.79 103.69 3.9 + + 103.44 22.45 103.03 28/11/2001 3/4/2002 28/11/2001 29/8/2002 3.79 2.92 4.67 14/4/2002 21.766 103.59 2.3 + 103.44 22.65 103.678 15/8/2001 26/2/2003 15/8/2001 5/5/2005 2.93 2.13 3.72 2/7/2002 22.796 103.639 3.3 + + 104.44 20.93 104.11 12/11/2001 26/7/2002 12/11/2001 27/6/2003 3.53 2.83 4.23 23/6/2002 21.13 104.525 3.01 + + 105.44 20.78 105.309 20/11/2001 5/4/2003 20/11/2001 29/6/2005 2.96 2.1 3.81 20/12/2001 20.8 105.326 3.01 + + Bảng 4.4.2 Trường hợp lựa chọn xác suất dự báo Pg = 90% 102.44 22.47 102.6 28/11/2001 18/5/2002 28/11/2001 19/7/2003 3.9 2.83 4.98 15/3/2002 21.84 102.305 3.01 + + 102.44 22.81 102 1/11/2000 6/9/2002 1/11/2000 8/2/2005 3.36 2.08 4.63 24/2/2002 22.62 102.332 3.01 + + 103.44 21.39 103.046 13/11/2001 20/6/2002 13/11/2001 20/10/2003 3.01 1.55 4.47 22/1/2002 20.79 103.69 3.9 + + 103.44 22.45 103.03 28/11/2001 3/4/2002 28/11/2001 24/11/2002 3.79 2.4 5.19 14/4/2002 21.766 103.59 2.3 + 103.44 22.65 103.678 15/8/2001 26/2/2003 15/8/2001 17/8/2006 2.93 1.66 4.19 2/7/2002 22.796 103.639 3.3 + + 104.44 20.93 104.11 12/11/2001 26/7/2002 12/11/2001 9/1/2004 3.53 2.42 4.64 23/6/2002 21.13 104.525 3.01 + + 105.44 20.78 105.309 20/11/2001 5/4/2003 20/11/2001 20/10/2006 2.96 1.6 4.31 20/12/2001 20.8 105.326 3.01 + + Bảng 4.4.3 So sánh kết dự báo động đất trường hợp lựa chon xác suất dự báo khác 0 0 -158 0.40 -0.40 0 0 0 0 116 -327 0.47 -0.47 0 0 0 0 -180 0.54 -0.54 0 0 0 0 -87 0.52 -0.52 0 0 0 0 -469 0.47 -0.47 0 0 0 0 -196 0.41 -0.41 0 0 0 0 -478 0.50 -0.50 0 0 λOp 50 4.5 Nhận xét Qua việc áp dụng thử nghiệm chƣơng trình dự báo động đất theo mô hình thống kê trƣờng hợp khác rút số nhận xét nhƣ sau: - Chƣơng trình cho kết khác lựa chọn tham số đầu vào khác Tùy theo mục đích yêu cầu nhiệm vụ nghiên cứu, chƣơng trình cho phép lựa chọn phƣơng án tính toán với độ xác (Pg) dự báo khác - Kết áp dụng thử nghiệm mô hình thống kê danh mục động đất khu vực Đông Nam Á giai đoạn 1997-2007 trƣờng hợp, đặc biệt lùi danh mục lại khứ năm, kiện dự báo (trùng với kiện xảy thực tế) chứng tỏ khả áp dụng chƣơng trình khu vực khác 51 KẾT LUẬN VÀ KIẾN NGHỊ Phƣơng pháp dự báo động đất theo mô hình thống kê có ƣu điểm lớn vừa xác định đƣợc vị trí chấn tâm, magnitude động đất, vừa xác định đƣợc thời gian xảy động đất mà lại độ tin cậy (xác suất Pg) dự báo Đã xây dựng đƣợc thuật toán sơ đồ khối chƣơng trình dự báo động đất theo mô hình thống kê Trên sở thuật toán sơ đồ khối mới, thiết lập đƣợc Chƣơng trình dự báo động đất theo mô hình thống kê Chƣơng trình dự báo động đất theo mô hình thống kê đƣợc thiết lập có ƣu điểm sau:  Giao diện thân thiện, dễ sử dụng,  Chƣơng trình cho phép lựa chọn tùy biến kích thƣớc tiểu vùng S độ tin cậy (hay xác suất Pg) dự báo tủy theo mục đích yêu cầu nhiệm vụ đặt nhƣ tùy thuộc vào đặc điểm đặc trƣng khu vực cần nghiên cứu  Không bị liệu : chƣơng trình có lỗi, ví dụ số liệu đầu vào lỗi, thông thƣờng chƣơng trình không lƣờng trƣớc tình phần mềm báo lỗi thoát, số liệu, tham số điền chƣơng trình có nguy bị chƣa kịp lƣu vào ổ cứng, nhƣng với chƣơng trình đƣợc thiết lập tác giả, hầu nhƣ loại lỗi không làm nội dung liệu đƣợc nhập vào  Dữ liệu đầu đƣợc hiển thị, bố trí thích hợp nên thuận tiện để sử dụng cho nhiều chƣơng trình khác nhau: Dữ liệu đầu đƣợc hiển thị dạng text, nhiên chúng đƣợc bó trí thích hợp nên cách bôi đen copy qua excel sử dụng đƣợc luôn, không cần phải cân chỉnh lại Áp dụng thử nghiệm chƣơng trình dự báo động đất theo mô hình thống kê trƣờng hợp khác cho phép rút số nhận xét nhƣ sau: - Chƣơng trình cho kết khác lựa chọn tham số đầu vào khác 52 Tùy theo mục đích yêu cầu nhiệm vụ nghiên cứu, chƣơng trình cho phép lựa chọn phƣơng án tính toán với độ xác (Pg) dự báo khác - Kết áp dụng thử nghiệm mô hình thống kê danh mục động đất khu vực Đông Nam Á giai đoạn 1997-2007 trƣờng hợp, đặc biệt lùi danh mục lại khứ năm, kiện dự báo (trùng với kiện xảy thực tế) chứng tỏ khả áp dụng chƣơng trình khu vực khác 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” Trong tƣơng lai, cần tiếp tục áp dụng thuƣnghiệm chƣơng trình đuợc thiết lập tác giả nhiều khu vực khác với nhiều phƣơng án khac để hoàn thiện chƣơng trình đƣa vào sử dụng phù hợp với điều kiện thực tiễn Kết dự báo động đất phƣơng pháp nêu túy mang chất xác xuất thống kê mà chƣa bao hàm chất vật lý kiện, nhiên lại cách tiếp cận với ƣu điểm bật nhƣ sau: Thuật toán sơ đồ khối chƣơng trình dự báo động đất sở ứng dụng mô hình thống kê vừa đơn giản mà lại có khả loại bỏ đƣợc nhƣợc điểm phƣơng pháp dự báo dài hạn trung hạn không xác định đƣợc độ xác dự báo động đất Đây khác biệt so với phƣơng pháp dự báo động đất khác ƣu điểm trội, vừa có ý nghĩa khoa học vừa có ý nghĩa thực tiễn cách tiếp cận Đồng thời nhƣợc điểm cách tiếp cận chƣa tính đến chất vật lý kiện khắc phục cách kết hợp phƣơng pháp dự báo động đất 53 mô hình thống kê với phƣơng pháp vật lý kiến tạo khác Ý tƣởng đƣợc tác giả luận văn tiếp tục thực tƣơng lai 54 TÀI LIỆU THAM KHẢO Tiếng Việt Allen K., Keilis-Borok V.I., Rotvain I.M nnk (1986), "Tập hợp tiền đề địa chấn dài hạn Kalifornia số khu vực khác", phƣơng pháp toán địa chấn địa động lực (Địa chấn tính toán V(19) Moscow, Nauka, 1987, tr 23-37 GS, TS Nguyễn Đình Xuyên (2004), Danh mục động đất Việt Nam (từ năm 114 đến năm 2003), báo cáo đề tài độc lập cấp nhà nước: “Nghiên cứu dự báo động đất dao động Việt Nam”, Hà Nội Dmitrieva O.E., Keilis-Borok V.I., Kosobocov V.G nnk, (1987), "Chẩn đoán chu kỳ xác suất cao động đất mạnh vùng nguy hiểm địa chấn Liên Xô (cũ) loạt nƣớc khác", Địa chấn tính toán, V(20) Moscow, Nauka, tr 99-111 Động đất miền Bắc Việt Nam, nha khí tƣợng xb, Hà Nội, 1968 Grishin A.P (2001) "Mô hình thống kê dự báo thời gian magnitude động đất", Địa chấn Núi lửa Viện hàn lâm khoa học Nga (4) Moscow, 2001, tr 60-65 Grishin A.P., Kondoskaya N.V., Levin L E., Solodinov L.N., Petrov A.L., Petrov O.M (2001)."Thử nghiệm dự báo thống kê vùng Kaspi (dự báo thời gian, tọa độ chấn tâm magnitude động đất)" Địa vật lý kỷ XXI Tuyển tập tài liệu đọc địa vật lý lần thứ mang tên Pheđưnski V.V Moscow, Thế giới khoa học, Tr 149-156 Grishin A.P., Kondorskaya N.V.(1997) "Thuật toán tách tiền chấn dƣ chấn từ danh mục động đất", Núi lửa địa chấn, (6) Viện hàn lâm khoa học Nga Matxcơva Tr 65-71 Gubin I E (1950) "Phƣơng pháp kiến tạo phân vùng động đất", công trình viện Vật lý địa cầu, Viện Hàn lâm khoa học Liên Xô, (13) Đặng Thanh Hải, Nguyễn Đức Vinh, Cao Đình Triều, (2002) “Dự báo dài hạn động đất khu vực Lai Châu – Điện Biên sở mô hình Thời gian – cấp độ mạnh , Tạp chí Khoa học Công nghệ, 40 (4), Hà Nội, tr.45 – 53 55 10 Đặng Thanh Hải (2003) Nghiên cứu đặc điểm cấu trúc sâu vỏ Trái đất phân vùng địa chấn kiến tạo miền Bắc Việt Nam, luận án Tiến sĩ Vật lý, 170 tr, Hà Nội 11 Đặng Thanh Hải (2006) Đặc trưng cấu trúc vỏ Trái đất xác định động đất cực đại vùng nguồn phát sinh động đất Việt Nam, báo cáo chuyên đề thuộc đề tài Nghiên cứu mã số 73.03.05, 23 tr, lƣu trữ Viện Vật lý địa cầu, Hà Nội 12 Đặng Thanh Hải, Cao Đình Triều (2006) “Đứt gãy hoạt động động đất miền Nam Việt Nam , Tạp chí Địa chất, loạt A(297) tr 11 - 23, Hà Nội 13 Keilis-Borok V.I., Kosobocov V.G (1986) "Các chu kỳ xác suất cao nảy sinh động đất mạnh giới" phƣơng pháp toán địa chấn địa động lực (Địa chấn tính toán), loạt V(19), Moscow, Nauka, 1987, tr 48-58 14 Kondorskaia N V., Bruk M G., Nosova O V (1986), "Xác định chấn tiêu dƣ chấn động đất Gazlisk 1984 phƣơng pháp nhóm chấn tiêu" Động đất Gazlisk Taskent FAN, C, tr 185-190 15 Kondorskaia N.V., Ngô Thị Lƣ, Niewiadomski J.(1994), "Về việc xác định tham số chấn tiêu động đất miền Bắc Việt Nam" Tin tức Viện hàn lân khoa học Nga, Tc Vật lý Địa cầu (4) Viện hàn lâm khoa học Nga, Moscow, 1994, tr 68-72 16 Kondorskaia N.V., Gorbunova I.V., Kiriev I.A Vandúieva N.V.(1993) " việc thành lập danh mục thống động đất mạnh khu vực Bắc á-âu theo tài liệu máy (1901-1990)" tính địa chấn phân vùng động đất khu vực Bắc á-âu Xbph Moscow 1993.:IFZ RAN Vật lý Địa cầu (4), tr 68-72 (Tiếng Nga) 17 Kosobocov V.G., Khili Dz., Dziuy Dz U I nnk.(1995) "Dự báo trung hạn sớm động đất Curill ngày 04 tháng 10 năm 1994 ngày tháng 12 năm 1995" Các vấn đề đại tính địa chấn địa động lực Trái đất (Địa chấn tính toán.), loạt V(28) Moscow, Nauka, 1996, tr 46-55 18 Kosobocov V.G., Mazhencov S.A.(1989), "Chẩn đoán chu kỳ độ nguy hiểm cao động đất mạnh Thiên Sơn theo thuật toán M8", Địa chấn tính toán, 56 Loạt V(22), Moscow, Nauka, 1989, tr 41-46 19 Nguyễn Kim Lạp, Nguyễn Duy Nuôi (1986) “ Độ nguy hiểm động đất khu vực Đông Nam Á , công trình khoa học Trung tâm nghiên cứu Vật lý địa cầu, Tập V(1985 - 1986), Viện Khoa học Việt Nam, Hà Nội, tr.7 – 13 57 [...]... chƣơng trình dự báo động đất bằng mô hình thống kê với việc sử dụng các thủ thuật toán học cần thiết 26 Chƣơng 3 THIẾT LẬP CHƢƠNG TRÌNH DỰ BÁO ĐỘNG ĐẤT THEO MÔ HÌNH THỐNG KÊ 3.1 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ê: Dựa trên cơ sở lý thuyết của phƣơng pháp đã trình bày ở chƣơng 2, chúng tôi đã xây dựng thuật toán của chƣơng trình dự báo động đất theo mô hình thống. .. động đất (magnitude) 25 Nhƣ vậy, trên cơ sở cách tiếp cận mới tác giả sẽ tiến hành xây dựng thuật toán và qui trình dự báo động đất theo mô hình thống kê Đây là một cách tiếp cận mới mang đặc trƣng thống kê Có thể áp dụng phƣơng pháp này phối hợp với các cách tiếp cận khác mang đặc trƣng vật lý kiến tạo sẽ cho kết quả khả quan trong nghiên cứu dự báo động đất Với thuật toán và sơ đồ khối đƣợc xây dựng, ... chúng (ti) Nhƣ vậy nhiệm vụ của mô hình thống kê là dự báo thời điểm xảy ra động đất Ti tại vùng S và độ mạnh (magnitude Mi) của nó Khi đó, thời điểm xảy ra động đất Ti và magnitude Mi của động đất dự báo đƣợc coi là sự kiện ngẫu nhiên Bài toán dự báo thời điểm xảy ra động đất Ti và magnitude Mi của nó chính là là bài toán xác định kỳ vọng toán học và các khoảng tin tƣởng của chúng với xác suất bảo... nghiệm dự báo ng n hạn động đất trên cơ s mô hình thống kê kết hợp sử dụng các phương pháp vật lý kiến tạo áp dụng đối v i l nh thổ Việt Nam và các vùng lân cận” nhằm giải quyết bài toán nan giải không chỉ đối với thực tế địa chấn Việt Nam nói riêng, mà cả trong thực tế địa chấn quốc tế nói chung Tác giả luận văn lựa chon đề tài: 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ô. .. thành thuật toán dự báo thời gian và magnitude của trận động đất xảy ra tiếp theo sự kiện tựa Các đại lƣợng và  Tpr  m pr đặc trƣng cho độ chính xác của việc dự báo trận động đất tiếp theo Dấu hiệu của độ tin cậy của dự báo động đất tiếp theo sau sự kiện tựa là các sự phụ thuộc sau: pr T  (Tmin ,Tmprax ) , pr m  (mmin , mmprax ) , Trong đó: T, m là thời điểm phát sinh và magnitude của động đất thực... các trận động đất [3], cũng nhƣ hàng loạt các dấu hiệu khác có đặc trƣng báo trƣớc các động đất tƣơng lai, và cũng sử dụng các bất thƣờng đồng thời trong diễn biến của một số các dấu hiệu có bản chất vật lý khác nhau [1, 18, 30, 31] Một số công trình theo hƣớng này chứa các công thức thuật toán dự báo động đất nhƣ: các thuật toán dự báo động đất KH và M8 [3, 13, 17] Kết quả dự báo đông đất theo các... đƣợc xác định tùy theo mục đích dự báo Trận động đất sau cùng này đƣợc gọi là trận động đất tựa (sự kiện tựa) và đƣợc ký hiệu là Zop Khi dự báo động đất mạnh, thƣờng thƣờng Mmin = 4-5 Trong các trƣờng hợp còn lại, giá trị Mmin đƣợc xác định bởi mức đại diện của danh mục động đất theo ngƣỡng magnitude Nhiệm vụ của mô hình là dự báo thời điểm xảy ra động đất và magnitude của trận động đất với chấn tâm... thúc và kiểm tra kết quả so với các trận t > topmax 7  xuất ra file danhmucTK.xml Nếu sốtrậntrongtiểuvùng(j) ... TRÌNH DỰ BÁO ĐỘNG ĐẤT THEO MÔ HÌNH THỐNG KÊ 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ê: Dựa sở lý thuyết phƣơng pháp trình bày chƣơng 2, xây dựng thuật toán chƣơng... chon đề tài: Xây dựng thuật toán sơ đồ khối chương trình dự báo động đất theo mô hình thống kê nhằm giải nội dung ý tƣởng nêu 14 Chƣơng PHƢƠNG PHÁP DỰ BÁO ĐỘNG ĐẤT THEO MÔ HÌNH THỐNG KÊ Trên giới,... cứu dự báo động đất Với thuật toán sơ đồ khối đƣợc xây dựng, thiết lập chƣơng trình dự báo động đất mô hình thống kê với việc sử dụng thủ thuật toán học cần thiết 26 Chƣơng THIẾT LẬP CHƢƠNG TRÌNH