Tập bản đồ xác suất nguy hiểm động đất Việt Nam và biển đông

Bài báo này trình bày phiên bản mới nhất của tập bản đồ xác suất nguy hiểm động đất lãnh thổ Việt Nam và khu vực Biển Đông. Danh mục động đất cập nhật tới năm 2014 và các thông tin mới được công bố gần đây nhất về địa chấn kiến tạo và địa động lực khu vực Đông Nam Á được sử dụng để xác định ranh giới của 37 vùng nguồn chấn động trên lãnh thổ Việt Nam và khu vực Biển Đông giới hạn bởi kinh tuyến 1250 Đông.

Tạp chí Khoa học Cơng nghệ Biển; Tập 15, Số 1; 2015: 77-90 DOI: 10.15625/1859-3097/15/1/6083 http://www.vjs.ac.vn/index.php/jmst TẬP BẢN ĐỒ XÁC SUẤT NGUY HIỂM ĐỘNG ĐẤT VIỆT NAM VÀ BIỂN ĐÔNG Nguyễn Hồng Phương*, Phạm Thế Truyền Viện Vật lý Địa cầu-Viện Hàn lâm khoa học công nghệ Việt Nam * E-mail: phuong.dongdat@gmail.com Ngày nhận bài: 23-10-2014 TÓM TẮT: Bài báo trình bày phiên tập đồ xác suất nguy hiểm động đất lãnh thổ Việt Nam khu vực Biển Đông Danh mục động đất cập nhật tới năm 2014 thông tin công bố gần địa chấn kiến tạo địa động lực khu vực Đông Nam Á sử dụng để xác định ranh giới 37 vùng nguồn chấn động lãnh thổ Việt Nam khu vực Biển Đông giới hạn kinh tuyến 1250 Đơng Mơ hình tắt dần chấn động Toro cộng (1997) áp dụng cho vùng nguồn chấn động lãnh thổ thềm lục địa Việt Nam, mơ hình tắt dần chấn động xây dựng cho đới hút chìm Youngs, Chiou, Silva Humphrey (1997) áp dụng cho vùng nguồn máng biển sâu Manila Các đồ nguy hiểm động đất biểu thị phân bố không gian giá trị trung vị gia tốc cực đại (PGA) với xác suất bị vượt 10%, 5%, 2% 0,5% vòng 50 năm Các khu vực có độ nguy hiểm cao bao gồm 1) vùng Tây Bắc Việt Nam, nơi có vùng nguồn chấn động Điện Biên Lai Châu Sơn La, có giá trị PGA cực đại đạt tới 180 g 272 g tương ứng với chu kỳ thời gian từ 475 năm 9975 năm; 2) khơi Nam Trung Bộ, nơi có vùng nguồn chấn động kinh tuyến 1090 Cửu Long - Côn Sơn, có giá trị PGA cực đại đạt tới 118 g 285 g tương ứng với chu kỳ thời gian từ 475 năm 9975 năm Tập đồ xác suất nguy hiểm động đất cung cấp thông tin dự báo định lượng ngắn hạn, trung bình dài hạn độ nguy hiểm động đất lãnh thổ Việt Nam khu vực Biển Đơng sử dụng thiết kế kháng chấn nhiều ứng dụng địa chấn cơng trình Từ khóa: Bản đồ xác suất nguy hiểm động đất, gia tốc cực đại nền, vùng nguồn chấn động, mơ hình tắt dần chấn động MỞ ĐẦU Độ nguy hiểm động đất định nghĩa xác suất xuất giá trị rung động động đất gây điểm bề mặt Trái đất Các kết đánh giá độ nguy hiểm động đất sử dụng liệu đầu vào cho tính tốn đánh giá rủi ro động đất, kết hợp với kết đánh giá hiệu ứng địa phương (như khuếch đại rung động địa phương gây trầm tích bở rời bề mặt, cấu trúc địa chất hay địa hình địa phương) yếu tố dễ bị tổn thương (như thể loại, kết cấu hay tuổi thọ cơng trình chịu rủi ro) để ước lượng thiệt hại động đất gây khu vực xét Vì vậy, việc tính tốn vẽ đồ xác suất nguy hiểm động đất ln đóng vai trị quan trọng thiếu việc đánh giá độ nguy hiểm động đất cho khu vực nghiên cứu Ở phạm vi quốc gia hay khu vực, đồ nguy hiểm động đất thường xây dựng phương pháp xác suất biểu diễn phân bố không gian số tham số rung động gia tốc, vận tốc hay dịch chuyển Ở Việt Nam, tham số rung động thường lựa chọn để xây dựng đồ xác suất nguy hiểm động đất đại lượng gia tốc cực đại (viết tắt tiếng Anh PGA), 77 Nguyễn Hồng Phương, Phạm Thế Truyền đo đơn vị gal hay cm/s2 Giá trị thực tiễn tham số chỗ, thường sử dụng trực tiếp làm đầu vào cho tính tốn tải trọng cơng trình thiết kế kháng chấn Tập đồ xác suất nguy hiểm động đất Nguyễn Hồng Phương thành lập cho lãnh thổ Việt Nam năm 1993 [1] Dựa đồ phân vùng địa chấn kiến tạo, thuật tốn A C Cornell chương trình EQRISK R K McGuire sử dụng để tính tốn lập đồ gia tốc cực đại cho lãnh thổ Việt Nam Đây lần mô hình xác suất áp dụng tồn quy trình đánh giá độ nguy hiểm động đất lãnh thổ Việt Nam, từ khâu xử lý thống kê tài liệu động đất, ước lượng tham số cho vùng nguồn đến tính tốn thành lập đồ gia tốc cực đại Tuy nhiên, hạn chế số liệu động đất hiểu biết kiến tạo khu vực lúc giờ, toàn lãnh thổ Việt Nam có vùng nguồn chấn động xác định, tồn phần lãnh thổ miền Nam Việt Nam sử dụng vùng nguồn đơn lẻ Ngoài ra, ảnh hưởng chấn động từ phía Biển Đơng chưa xét đến cơng trình Nhược điểm nêu khắc phục đồ gia tốc cực đại thành lập cho lãnh thổ Việt Nam vùng lân cận công bố năm 1997 tác giả Trong cơng trình này, ảnh hưởng lan truyền chấn động từ quốc gia lân cận tới lãnh thổ Việt Nam tính đến việc bổ sung thêm vùng nguồn chấn động lãnh thổ Trung Quốc, Lào, Myanmar Thái Lan Ngoài ra, để tính đến ảnh hưởng lan truyền chấn động từ vùng biển tới lãnh thổ Việt Nam, ranh giới vùng nguồn chấn động Ấn Độ Dương vùng nguồn khác khu vực Thái Bình Dương xác định Tổng cộng có 17 vùng nguồn chấn động sử dụng cho tính toán vẽ đồ xác suất nguy hiểm động đất sử dụng cơng trình [2] Đến năm 2004, tập đồ tác giả cập nhật lại lần nữa, với việc hiệu chỉnh lại ranh giới vùng nguồn khu vực Biển Đông Việt Nam [3] Trên giới, vào năm 1998, nhân thập kỷ giảm nhẹ thiên tai toàn giới (IDNDR), Chương trình đánh giá độ nguy hiểm động đất tồn cầu Liên hợp quốc (GSHAP) cho đời đồ phân bố gia tốc thành lập cho 78 phạm vi toàn giới [4] Trên đồ này, khu vực Việt Nam chiếm vị trí không đáng kể, vùng nguồn chấn động có ảnh hưởng tới Việt Nam lại khơng phải nhà địa chấn Việt Nam thành lập Cũng cần nói thêm rằng, hiệu ứng chấn động biển lãnh thổ Việt Nam khơng tính đến cơng trình Từ năm 2000, nhiều đề tài, dự án với quy mô lớn triển khai thực lĩnh vực nghiên cứu địa chấn cho phép thu thập bổ sung thêm nhiều liệu động đất, địa chất vả địa vật lý, làm thay đổi đáng kể tranh địa chấn kiến tạo lãnh thổ Việt Nam khu vực Biển Đông Việc phát thêm chi tiết hóa hệ thống đới đứt gẫy sinh chấn lãnh thổ lãnh hải nước ta làm thay đổi đáng kể phân bố không gian vùng nguồn chấn động lãnh thổ vùng biển Việt Nam dẫn đến nhu cầu thiết yếu việc hiệu chỉnh đồ xác suất nguy hiểm động đất xây dựng Năm 2004, khuôn khổ đề tài nghiên cứu cấp Nhà nước Viện Vật lý Địa cầu, Nguyễn Đình Xuyên cộng sử dụng chương trình CRISIS99 để thành lập đồ PGA cho lãnh thổ Việt Nam với chu kỳ lặp lại chấn động 500 năm Ưu điểm đồ thành lập dựa đồ vùng nguồn chi tiết lãnh thổ thềm lục địa Việt Nam tính đến thời điểm Mặc dù vậy, cơng trình giới hạn khn khổ tính tốn lãnh thổ đới hẹp thềm lục địa Việt Nam mà khơng tính đến ảnh hưởng vùng nguồn chấn động khu vực Biển Đông Bản đồ sau sử dụng Tiêu chuẩn xây dựng Việt Nam TCXDVN 375:2006 ban hành vào năm 2006 [5] Năm 2010, khuôn khổ hai đề tài nghiên cứu cấp Nhà nước, Nguyễn Hồng Phương sử dụng chương trình CRISIS99 để cập nhật tập đồ xác suất nguy hiểm động đất cho lãnh thổ Việt Nam khu vực Biển Đông [6] Ưu điểm lớn đồ xây dựng sở sơ đồ vùng nguồn chấn động đầy đủ từ trước đến lãnh thổ Việt Nam lẫn khu vực Biển Đông Trên lãnh thổ Việt Nam, toàn vùng nguồn kế thừa từ đồ Nguyễn Đình Xuyên Nguyễn Hồng Phương năm 2004 [3, 5] Trên đồ này, vùng Tập đồ xác suất nguy hiểm động đất … nguồn khu vực Biển Đông vùng biển thuộc quốc gia lân cận Philippines, Malaysia, Thái Lan Đài Loan xây dựng chi tiết với độ tin cậy cao Tuy tập đồ đầy đủ từ trước đến nay, kết đồ năm 2010 công bố cho dải ven biển Việt Nam khu vực Biển Đông bao gồm hai quần đảo Hồng Sa Trường Sa Mục đích nghiên cứu hiệu chỉnh hoàn thiện tập đồ xác suất nguy hiểm động đất cho lãnh thổ Việt Nam khu vực Biển Đông xây dựng sở cập nhật số liệu động đất hiệu chỉnh ranh giới vùng nguồn chấn động có lưu ý tới tài liệu địa chấn kiến tạo khu vực Có thể coi cơng trình hiệu chỉnh tập đồ độ nguy hiểm động đất công bố năm 2010 [6] Do khuôn khổ viết, để tránh lặp lại thông tin công bố trước đây, tập trung mô tả thay đổi hiệu chỉnh áp dụng cho tập đồ xác suất nguy hiểm động đất thành lập CƠ SỞ LÝ THUYẾT Trong nghiên cứu này, phương pháp xác suất kinh điển đánh giá độ nguy hiểm động đất Cornell Esteva đề xuất lần năm 1968 áp dụng [7, 8] Phương pháp xác suất dựa giả thiết cho vùng nguồn chấn động khác xác định phạm vi khu vực nghiên cứu phát sinh trận động đất xảy phạm vi vùng nguồn có tính độc lập thống kê tần suất bị vượt λ(M) giá trị độ lớn động đất M bị vượt đơn vị thời gian tính cách sử dụng số liệu danh mục động đất Đại lượng λ(M) tính số trận động đất có độ lớn M bị vượt đơn vị thời gian đặc trưng cho tính địa chấn vùng nguồn chấn động Phương pháp xác suất giả thiết tất điểm nằm bên phạm vi vùng nguồn chấn động coi chấn tiêu động đất Khi đó, tần suất bị vượt giá trị gia tốc a gây nguồn điểm đơn lẻ (nguồn thứ i) tính theo công thức: vi  a   w j ij Mu M0   d dM M   Pr  A  a M , R  dM i ij (1) Ở M0 Mu giá trị cận cận độ lớn động đất, Pr(A>a|M, Rij) xác suất để gia tốc A điểm xét vượt giá trị a cho trước tác động trận động đất có độ lớn M có chấn tiêu nằm cách điểm tính khoảng Rij Trong chương trình tính tốn, Rij khoảng cách điểm tính đơn vị nguồn, tức thành phần nhỏ mà vùng nguồn chấn động chia Các trọng số wij gán cho đơn vị nguồn với giả thiết Σwij=1 Cuối cùng, ảnh hưởng tất N nguồn chấn động tới điểm xét tính công thức: v a   N i 1 vi  a  (2) Mơ hình tính địa chấn sử dụng nghiên cứu mơ hình Gutenberg-Richter cải tiến, theo tần suất bị vượt giá trị độ lớn động đất tính cơng thức:  M  M   M   0 e M  e  M u , M  M  M u e e (3) Trong λ0 tần suất bị vượt độ lớn động đất M0, β đại lượng tương đương với “tham số b” vùng nguồn chấn động Mu độ lớn động đất cực đại ước lượng cho vùng nguồn chấn động Với giả thiết phát sinh động đất vùng nguồn tn theo luật phân bố Pốt xơng, hàm mật độ xác suất đại lượng độ lớn động đất M cho bởi: p M     M d M   0  Me0  M u , M  M  M u dM e e (4) MƠ HÌNH NGUỒN CHẤN ĐỘNG Danh mục động đất Danh mục động đất tổng hợp thành lập cho khu vực nghiên cứu bao gồm 4.496 trận động đất, tập hợp từ nguồn số liệu sau đây: Danh mục động đất Việt Nam thời kỳ 114 - 2014, lưu trữ Trung tâm báo tin động đất cảnh báo sóng thần, Viện Vật lý Địa cầu, Viện Hàn lâm Khoa học Công nghệ Việt Nam Danh mục bao gồm số liệu 79 Nguyễn Hồng Phương, Phạm Thế Truyền động đất lịch sử số liệu đo máy, thu thập từ mạng lưới đài trạm quan trắc động đất Việt Nam trung tâm địa chấn giới ISC,USGS, NEIS, BEJ Danh mục động đất ghi nhận bên lãnh thổ Việt Nam vùng biển lân cận Việt Nam khoảng thời gian 1524 - 2014, giới hạn phạm vi φ=4,0 - 23,50N; λ=100 - 1250E) Các đới đứt gẫy hoạt động lãnh thổ Việt Nam Biển Đông Việc xác định đới đứt gẫy hoạt động có quy mơ, mức độ chế hoạt động khác có ý nghĩa lớn việc xác định vùng nguồn phát sinh động đất đánh giá mức độ nguy hiểm chúng vùng nghiên cứu Để phục vụ công tác đánh giá độ nguy hiểm động đất, đới đứt gẫy kiến tạo có biểu hoạt động đương đại tồn vùng nghiên cứu phân cấp theo vai trò phân chia yếu tố kiến trúc kiến tạo Việc phân cấp đứt gãy theo vai trò cấu trúc dựa luận thuyết kiến tạo mảng (mảng, vi mảng, khối cấu trúc ) tiến hóa địa động lực (ép, tách, trượt bằng) đứt gãy theo thời gian, độ sâu ảnh hưởng đứt gãy (vỏ, xuyên vỏ) Các đứt gãy cấp I ranh giới mảnh, vi mảng, đơn vị cấu trúc hệ địa động lực khác Kainozoi Các đứt gãy cấp II đứt gãy sinh kèm lông chim đứt gãy cấp I khu vực nghiên cứu, đóng vai trị ranh giới khối cấu trúc địa động lực, đới cấu trúc - địa động lực, phụ khối cấu trúc - địa động lực Các đứt gãy cấp III đứt gãy lông chim sinh kèm đứt gãy cấp II đóng vai trị ranh giới phụ khối - địa động lực, thường bị khống chế giới hạn khối - địa động lực miền cấu trúc - địa động lực Những đới đứt gẫy có quy mơ lớn (bậc I, II, III) với độ sâu ≥ 15 - 20 km đóng vai trị ranh giới khối kiến tạo có quy mơ khác đới có khả phát sinh động đất Các vùng kiến tạo lãnh thổ Việt Nam bị chia cắt thành đới nhỏ hơn, thường hẹp kéo dài theo phương cấu trúc địa chất, mạng lưới đới đứt gãy có kích thước khác Mức độ hoạt động đứt gãy đánh giá thực địa dựa dấu hiệu: địa chất, địa mạo, viễn thám, xuất lộ 80 nguồn nóng - nước khống, nứt - trượt đất, dị thường khí Radon, thủy ngân, mêtan, cacbonic, dị thường địa nhiệt Dấu hiệu tin cậy cao biểu hoạt động động đất thời kỳ đại hay biên độ dịch trượt tính từ đo lặp trắc địa xác Đơi biến vị thành tạo trầm tích Đệ tứ khi bị đứt gãy cắt qua cho phép xác định tuổi phương pháp huỳnh quang kích hoạt (OSL) Năm 2007, để xây dựng mơ hình nguồn tuyến phát sinh động đất phục vụ cho việc đánh giá độ nguy hiểm độ rủi ro động đất Việt Nam, Nguyễn Hồng Phương Phạm Thế Truyền thu thập, tổng hợp thành lập danh mục bao gồm 46 hệ thống đứt gẫy sinh chấn xác định toàn lãnh thổ thềm lục địa Việt Nam [9] Các đứt gẫy phân thành hai nhóm với cấp độ hoạt động khác nhau, cấp bao gồm đứt gẫy phân chia miền kiến tạo (mảng bậc 2) cấp bao gồm đứt gẫy đóng vai trị phân chia cấu trúc miền Từ năm 2010, nghiên cứu chi tiết đặc trưng địa chấn kiến tạo khu vực Đông Nam Á làm sáng tỏ chế phát sinh động đất sóng thần khu vực Biển Đông Việt Nam vùng biển lân cận [6, 10, 11] Danh mục đứt gẫy sinh chấn lãnh thổ Việt Nam nhờ chi tiết hóa bổ sung thêm tham số tính phân đoạn đứt gẫy hoạt động Ngoài ra, đới hút chìm đứt gẫy hoạt động phạm vi khu vực vùng Biển Đông Việt Nam nghiên cứu chi tiết bổ sung vào danh mục hệ thống đứt gẫy sinh chấn Bản đồ địa chấn kiến tạo lãnh thổ Việt Nam khu vực Biển Đơng minh họa hình cho thấy đới hút chìm đứt gẫy khu vực sau có khả phát sinh động đất mạnh ảnh hưởng tới Việt Nam: Đới hút chìm máng biển sâu Manila dài 1.150 km chạy dọc bờ tây quần đảo Philippines Cho đến nay, đới hút chìm Manila ghi nhận động đất có độ lớn vượt độ theo thang Mô men số trận sóng thần gây thiệt hại cho khu vực ven biển Philippines Đối với Việt Nam, đới hút chìm Manila coi nguồn phát sinh động đất mạnh nguy hiểm nhất, phát sinh sóng thần gây ảnh hưởng tới vùng bờ biển Việt Nam [10, 11] Tập đồ xác suất nguy hiểm động đất … Các đứt gẫy trượt có kích thước khu vực như: đứt gẫy kinh tuyến 1090 có chế trượt phải giai đoạn chạy dọc dải ven biển miền Trung Việt Nam, đứt gẫy vĩ tuyến Nam Hải Nam đứt gẫy Ba Chùa chạy theo phương tây bắc - đơng nam Hình Bản đồ địa chấn kiến tạo Việt Nam Biển Đông, số liệu động đất bao gồm động đất lịch sử, cập nhật đến năm 2014 Kèm theo đới đứt gẫy nêu cịn có số đới đứt gẫy có kích thước phạm vi nhỏ hơn, có khả sinh chấn Cửu Long - Côn Sơn, Thuận Hải - Minh Hải nằm thềm lục địa đông nam Việt Nam Borneo - Palawan nằm vùng biển Malaysia Philippines Các vùng nguồn chấn động Trên sở danh mục động đất thông tin đới đứt gẫy hoạt động biết, vùng nguồn chấn động xác định lãnh thổ Việt Nam Biển Đông Các vùng phát sinh động đất coi tổng cộng vùng cực động tất trận động đất cực đại có khả xảy đới phá hủy kiến tạo Đó hình chiếu mặt đứt gãy kiến tạo (kể từ ranh giới bên tầng hoạt động) lên mặt đất Tuy nhiên, nhiều trường hợp, điều kiện số liệu động đất không đầy đủ, việc đánh giá mức độ hoạt động đứt gẫy dựa số liệu địa vật lý, địa chất, chứng tự nhiên Ranh giới vùng nguồn xác định theo mật độ phân bố chấn tâm quan sát được, hay vào tổ hợp phân bố đứt gãy, 81 Nguyễn Hồng Phương, Phạm Thế Truyền cung núi lửa liên quan Ranh giới cuối nhận này, phản ánh trung thực đặc trưng địa chấn kiến tạo đới nằm, phương cấu trúc phân bố không gian chấn tâm, xác định vùng nguồn động đất khu vực nghiên cứu Hình Sơ đồ vùng nguồn phát sinh động đất lãnh thổ Việt Nam khu vực Biển Đông, (Các số đồ trùng với số thứ tự vùng nguồn liệt kê bảng 2) Áp dụng nguyên tắc trên, ranh giới 37 vùng nguồn chấn động vạch khu vực nghiên cứu gồm toàn lãnh thổ Việt Nam khu vực Biển Đơng Đây mơ hình vùng nguồn đầy đủ nhất, hiệu chỉnh sở kế thừa nghiên cứu trước Nguyễn Đình Xuyên (2004) Nguyễn Hồng Phương (2004, 2010) [3, 5, 6] Ở đây, giới hạn phía đơng khu vực nghiên cứu lấy tới sát bờ tây quần đảo Philippines Như vậy, 82 có tất 30 vùng nguồn chấn động nằm phần lục địa lãnh thổ Việt Nam vùng nguồn chấn động nằm Biển Đơng Trên hình minh họa đồ vùng nguồn chấn động sử dụng nghiên cứu Cần nhận xét hai vùng nguồn máng biển sâu Manila Palawan Biển Đông nằm cách bờ biển Việt Nam 1.000 km gần không gây ảnh hưởng chấn động tới lãnh thổ Việt Nam Hai vùng nguồn đưa vào tính Tập đồ xác suất nguy hiểm động đất … tốn với mục đích khơng bỏ sót chấn động tác động tới hai quần đảo Hoàng Sa Trường Sa Ước lượng tham số nguy hiểm động đất cho vùng nguồn chấn động Việc xử lý số liệu động đất đóng vai trị quan trọng quy trình đánh giá độ nguy hiểm động đất, đặc biệt trường hợp áp dụng phương pháp xác suất, nguyên lý quan trọng phương pháp đánh giá độ nguy hiểm động đất theo cách tiếp cận xác suất trận động đất sử dụng vào việc tính tốn phải kiện độc lập với mặt thống kê Vì thế, số liệu động đất sau đồng đơn vị đo độ lớn nhóm theo vùng nguồn phải tiếp tục xử lý để loại bỏ tất tiền chấn dư chấn, giữ lại rung động chuỗi số liệu động đất Chỉ có động đất có độ lớn từ 4,0 trở lên sử dụng cho tính tốn Ngun lý loại bỏ dư chấn (hoặc tiền chấn) biết đến rộng rãi Giả sử t thời điểm xảy động đất, h độ sâu chấn tiêu, M độ lớn động đất, i j số thứ tự hai trận động đất danh mục, j > i Trận động đất thứ hai coi dư chấn trận động đất thứ điều kiện sau thoả mãn: khoảng cách chấn tâm hai trận động đất nhỏ giá trị cho trước R(Mi), hj - hi ≤ H(M); Mj ≤ Mi, với T(M), R(M) H(M) hàm thực nghiệm [12] Thuật toán loại bỏ tiền chấn hoàn toàn tương tự Độ lớn động đất cực đại dự báo Mmax; Các tham số a b biểu thức tương quan độ lớn tần suất động đất (biểu thức Gutenberg - Richter) đại lượng suy diễn tương ứng λ β Chu kỳ lặp lại dự báo T(M) động đất mạnh vùng Các phương pháp hợp lý cực đại cực trị áp dụng đồng thời để ước lượng tham số nguy hiểm động đất Cơ sở lý thuyết hai phương pháp mơ tả chi tiết nhiều cơng trình trước [1, 3, 6, 10] Các nghiên cứu trước cho thấy phương pháp hợp lý cực đại cho kết ước lượng tham số nguy hiểm động đất xác thực so với phương pháp cực trị, đặc biệt giá trị Mmax Vì vậy, nghiên cứu này, phương pháp hợp lý cực đại ưu tiên áp dụng cho vùng nguồn chấn động Trong trường hợp số liệu động đất thu thập vùng nguồn q để áp dụng phương pháp xác suất - thống kê, phương pháp “tương đồng địa chấn kiến tạo” áp dụng Nguyên tắc phương pháp phát biểu sau: cho hai vùng nguồn động đất có điều kiện địa chấn kiến tạo địa động lực tương đồng, vùng nguồn khơng có đủ số liệu động đất để ước lượng tham số nguy hiểm động đất vùng nguồn có tham số nguy hiểm địa chấn ước lượng, cho phép gán tham số vùng nguồn cho vùng nguồn [5, 6] Trong cơng trình này, phần mềm Main.exe sử dụng để loại bỏ tự động tiền chấn dư chấn chuỗi số liệu động đất ứng với vùng nguồn chấn động Các chuỗi số liệu đưa vào tính tốn chứa tồn động đất để đảm bảo độ tin cậy cho kết tính tốn Phương pháp “tương đồng địa chất kiến tạo” áp dụng để ước lượng tham số nguy hiểm động đất cho vùng nguồn Tuy Hịa - Củ Chi Ngồi khơng ghi nhận số liệu động đất đủ lớn, vùng nguồn sơng Sài Gịn gán tham số vùng nguồn Sơng Vàm Cỏ Đơng, coi có điều kiện địa chất kiến tạo tương tự khu vực nghiên cứu Để phục vụ tính tốn vẽ đồ độ nguy hiểm động đất, tham số sau ước lượng cho vùng nguồn: Kết ước lượng tham số động đất cho vùng nguồn lãnh thổ Việt Nam Biển Đông liệt kê bảng 83 Nguyễn Hồng Phương, Phạm Thế Truyền Bảng Tham số nguy hiểm động đất vùng nguồn lãnh thổ Việt Nam Biển Đông TT Tên vùng nguồn 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36a 36b 36c 37 Sơn La Sông Mã - Pumaytun Đơng Triều - ng Bí Sơng Hồng - Sông Chảy Văn Sơn - Hà Giang Thái Nguyên - Bắc Cạn Cao Bằng - Tiên Yên Cẩm Phả Đông Bắc Trũng Hà Nội Sông Hiếu Sông Lô Quốc Lộ 13A Phong Thổ - Thanh Sơn Mường La - Bắc Yên Sông Đà Lai Châu - Điện Biên Mường Tè Mường Nhé Sông Cả - Khe Bố Rào Nậy Dakrong - Huế Trà Bồng Đà Nẵng Tam Kỳ - Phước Sơn Ba Tơ - Củng Sơn Tuy Hòa - Củ Chi Sơng Sài Gịn Nha Trang - Tánh Linh Sông Vàm Cỏ Đông Sông Hậu Kinh tuyến 109 Thuận Hải - Minh Hải Cửu Long - Côn Sơn Ba Tháp Bắc Biển Đông Máng biển sâu Manila Bắc Máng biển sâu Manila Trung Máng biển sâu Manila Nam Palawan 0,11 0,227 0,084 0,051 0,117 0,13 0,034 0,006 0,011 0,10 0,04 0,111 0,121 0,087 0,238 0,083 0,476 0,016 0,027 0,016 0,02 0,02 0,034 0,02 0,02 0,4 0,02 0,02 0,437 0,434 0,181 0,02 0,306 4,72 6,04 1,4 0,285 N Mmax_obs Mmax_HLCD M0 bHLCD H (km) 45 18 37 4 10 21 10 14 2 1 21 27 236 490 28 6,7 6,8 5,6 5,8 4,6 5,2 5,0 4,8 5,5 5,2 4,8 4,8 5,1 4,9 4,8 5,6 4,7 5,3 6,0 4,2 4,8 4,8 4,7 5,3 4,8 4,0 4,4 6,1 5,1 5,1 4,5 6,5 8,2 8,0 6,2 6,0 7,2 ± 0,54 7,3 ± 0,77 6,1 ± 0,54 6,3 ± 0,54 5,1 ± 0,82 5,7 ± 0,26 5,5 ± 0,92 7,3 ± 0,24 6,0 ± 0,54 5,7 ± 0,37 5,3 ± 0,21 5,3 ± 0,34 5,6 ± 0,23 5,4 ± 0,54 5,3 ± 0,54 6,5 ± 0,54 5,2 ± 0,45 5,8 ± 0,54 6,5 ± 0,54 6,0 ± 0,51 5,3 6,2 ± 0,54 5,3 5,2 5,8 ± 0,54 5,3 4,5 ± 1,56 4,5 ± 0,47 4,5 ± 1,56 4,9 ± 0,35 6,6 ± 0,28 5,6 ± 0,30 5,7 ± 0,38 5,1 7,0 ± 0,23 8,7 ± 0,93 8,5 ± 0,85 6,7 ± 0,28 6,4 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 4,0 5,0 5,0 5,0 4,0 0,49 0,59 0,48 1,00 1,21 0,33 1,18 0,16 0,63 0,61 0,25 0,58 0,29 0,44 1,30 0,32 0,79 0,66 0,74 0,58 1,0 0,63 1,0 1,0 0,14 1,0 1,16 0,59 1,16 0,36 0,28 0,32 0,51 22 22 22 17 10 10 12 12 12 12 12 10 10 12 12 12 12 12 17 12 12 12 12 12 12 10 10 10 10 12 12 12 10 10 33 35 35 35 30 0,30 0,65 0,88 0,56 Chú thích: - tần suất động đất có độ lớn M0 bị vượt hàng năm; N - số động đất ghi nhận phạm vi vùng nguồn; Mmax_obs - giá trị độ lớn động đất cực đại quan sát; Mmax_HLCĐ - giá trị độ lớn động đất cực đại tính phương pháp hợp lý cực đại; M0 - cận giá trị độ lớn động đất; bHLCD - giá trị b suy diễn từ kết phương pháp hợp lý cực đại; H - bề dầy tầng hoạt động vùng nguồn chấn động QUY LUẬT TẮT DẦN CHẤN ĐỘNG Việc xác lập biểu thức tắt dần chấn động cho khu vực nghiên cứu đóng vai trị quan trọng thường xem xét bước 84 độc lập quy trình đánh giá độ nguy hiểm động đất Ở Việt Nam, thiếu số liệu gia tốc từ trận động đất mạnh nên thời gian dài việc đánh giá nguy hiểm động đất phải sử dụng biểu thức tắt dần chấn động Tập đồ xác suất nguy hiểm động đất … xây dựng cho khu vực khác giới Mặc dù từ năm 2011, có hai biểu thức tắt dần chấn động xây dựng cho miền Bắc Việt Nam, chúng chưa áp dụng rộng rãi độ tin cậy biểu thức kiểm nghiệm [13, 14] Cho đến nay, có nhiều biểu thức tắt dần chấn động khác áp dụng để đánh giá độ nguy hiểm động đất cho lãnh thổ Việt Nam Bản đồ xác suất nguy hiểm động đất Nguyễn Hồng Phương thành lập năm 1993 sử dụng biểu thức tắt dần chấn động McGuire [1, 15] Tập đồ tác giả hiệu chỉnh vào năm 1997 2004, sử dụng biểu thức tắt dần chấn động xây dựng cho khu vực Vân Nam, Trung Quốc [2, 3, 25] Năm 2004, đồ xác suất nguy hiểm động đất lãnh thổ Việt Nam Nguyễn Đình Xuyên cộng thành lập áp dụng đồng thời hai biểu thức tắt dần chấn động Campbell (1997) Toro cộng (1997) [5] Năm 2010, Nguyễn Hồng Phương cộng sử dụng biểu thức tắt dần chấn động Campbell (1997) để thành lập đồ xác suất nguy hiểm động đất cho dải ven biển Việt Nam Biển Đông [6], đến năm 2014, tập thể tác giả sử dụng mơ hình NGA08 để thành lập đồ xác suất nguy hiểm động đất cho vùng Nam Trung Bộ Trong nghiên cứu này, biểu thức tắt dần chấn động Toro, Abrahamson Schneider (1997) sử dụng để thành lập tập đồ xác suất nguy hiểm động đất cho lãnh thổ Việt Nam Biểu thức xây dựng cho miền Đông Hoa Kỳ, khu vực bình ổn tính địa chấn so với miền Tây nước Mỹ theo phù hợp với điều kiện Việt Nam Biểu thức tắt dần chấn động Toro, Abrahamson Schneider (1997) có dạng tổng qt sau [23]: mơ hình lý thuyết, M độ lớn động đất theo thang mô men, RM khoảng cách chấn tiêu, Rjb khoảng cách gần theo phương nằm ngang tới đứt gẫy (còn gọi khoảng cách Joyner-Boor), đo km Cũng cần lưu ý rằng, việc áp dụng biểu thức tắt dần chấn động cho khu vực rộng lớn dẫn đến sai số kết ước lượng độ nguy hiểm động đất Vì vậy, vùng nguồn máng biển sâu Manila, nghiên cứu sử dụng quy luật tắt dần chấn động xây dựng cho đới hút chìm Youngs, Chiou, Silva Humphrey (1997) đề xuất [26]: LnY  aIF.GIF  aIS.GIS  1, 414M  b(10  M )3 (7)  c.ln  R  1, 782e 0,554M   0, 00607H Ở Y giá trị gia tốc cực đại (đo g), M độ lớn động đất theo thang mô men, R khoảng cách gần tới hình chiếu đứt gẫy mặt đất, aIF, aIS hệ số ứng với trường hợp động đất nội mảng rìa mảng, GIF, GIS số loại nguồn, GIF=1 động đất nội mảng GIF=0 trường hợp ngược lại, GIS=1 động đất rìa mảng GIS=0 trường hợp ngược lại, H độ sâu chấn tiêu, đo km TẬP BẢN ĐỒ ĐỘ NGUY HIỂM ĐỘNG ĐẤT LÃNH THỔ VIỆT NAM VÀ BIỂN ĐƠNG Chương trình CRISIS2007 Ordaz cộng áp dụng để tính tốn độ nguy hiểm động đất cho khu vực nghiên cứu, với số liệu đầu vào sau [16]: Bản đồ vùng nguồn chấn động khu vực nghiên cứu (hình 3) Các tham số nguy hiểm động đất tham số phục vụ tính tốn độ nguy hiểm động đất xác định cho vùng nguồn (bảng 2) ln Y  c1  c2  M    c3  M    c4 ln RM   c5  c6  max ln RM ,0   100   c6 RM   e   a   RM  R 2jb  C72 (5) (6) Ở Y giá trị gia tốc cực đại (đo g), hệ số từ c1 đến c7 xác định Bảng liệt kê toàn 37 vùng nguồn chấn động tham số tương ứng sử dụng việc tính tốn thành lập đồ độ nguy hiểm động đất cho lãnh thổ Việt Nam khu vực Biển Đông Đối với tất vùng nguồn, giá trị cận độ lớn động đất chọn M0=4,0 Gia tốc cực đại (đo đơn vị % gal) tính điểm mạng 85 Nguyễn Hồng Phương, Phạm Thế Truyền lưới 0,20 × 0,20 phủ lên tồn vùng nghiên cứu sử dụng để xây dựng đồ độ nguy hiểm động đất Hình Bản đồ gia tốc cực đại (PGA) lãnh thổ Việt Nam Biển Đơng tính cho đá với 10% xác suất bị vượt vòng 50 năm Trên hình 3, 4, minh họa tập đồ xác suất nguy hiểm động đất thành lập cho lãnh thổ Việt Nam khu vực Biển Đông với ranh giới phía đơng khu vực nghiên cứu thu nhỏ đến kinh tuyến 1250N Các đồ biểu thị phân bố không gian đại lượng gia tốc cực đại (PGA) đo đơn vị cm/s2 (hay gal) tính cho đá, tương ứng với 10%, 5%, 2% 0,5% xác suất bị vượt vịng 50 năm Với giả thiết phân bố Pốt xông phát sinh động đất, giá trị PGA thể đồ ứng với chu kỳ lặp lại 475, 975, 2475 9975 năm [3] Từ đồ kết đưa số nhận định sau đây: Phân bố không gian giá trị PGA đồ phản ánh rõ nét ranh giới vùng nguồn chấn động lãnh thổ Việt Nam khu vực Biển Đông Cường độ 86 chấn động mạnh tập trung vùng Tây Bắc lãnh thổ Việt Nam, nơi có đới đứt gẫy Điện Biên - Lai Châu Sông Mã thềm lục địa Nam Trung Bộ, nơi có đới đứt gẫy Kinh tuyến 1090 Hình Bản đồ gia tốc cực đại (PGA) lãnh thổ Việt Nam Biển Đông tính cho đá với 5% xác suất bị vượt vòng 50 năm Trên lãnh thổ Việt Nam, rung động mạnh quan sát thấy khu vực tây bắc đất nước, tập trung hai vùng nguồn Điện Biên - Lai Châu Sơn La Ứng với chu kỳ lặp lại 475, 975, 2475 9975 năm, gia tốc cực đại vùng nguồn Điện Biên - Lai Châu đạt giá trị 180 cm/s2, 202 cm/s2, 235 cm/s2 272 cm/s2, vùng nguồn Sơn La đạt giá trị 194 cm/s2, 210 cm/s2, 232 cm/s2 264 cm/s2 Cần lưu ý chu kỳ lặp lại khoảng 10.000 năm, giá trị rung động hai vùng nguồn nêu tương đương với cường độ chấn động bề mặt cấp IX theo thang MSK-64 Ngoài hai vùng chấn động mạnh nêu trên, lãnh thổ Việt Nam quan sát thấy Tập đồ xác suất nguy hiểm động đất … loạt vùng nguồn có rung động mạnh Ứng với chu kỳ lặp lại 9975 năm, vùng nguồn có rung động mạnh tới cấp VIII liệt kê bao gồm : Sông Hồng - Sông Chảy, Rào Nậy, Sông Cả - Khe Bố Trà Bồng đương với cường độ chấn động mặt tới cấp VIII theo thang MSK-64 Trên vùng thềm lục địa miền Trung Việt Nam, chấn động với gia tốc cực đại 122 gal (tương đương cường độ trung bình tới cấp VIII theo thang MSK-64) ứng với chu kỳ thời gian 9975 năm xuất vùng nguồn Bắc Biển Đông Khả rung động hai vùng quần đảo Hoàng Sa Trường Sa không cao Giá trị gia tốc cực đại hai khu vực dao động khoảng từ 28 đến 44 gal cho tất chu kỳ xét, tương đương với cường độ chấn động mặt tới cấp VI theo thang MSK-64 Hình Bản đồ gia tốc cực đại (PGA) lãnh thổ Việt Nam Biển Đơng tính cho đá với 2% xác suất bị vượt vòng 50 năm Giá trị gia tốc cực đại lớn quan sát thấy thềm lục địa Nam Trung Bộ, vùng nguồn Kinh tuyến 1090 Ứng với chu kỳ lặp lại 475, 975, 2475 9975 năm, gia tốc cực đại vùng nguồn Kinh tuyến 1090 đạt giá trị 118, 151 cm/s2, 200 cm/s2 285 cm/s2 Cần lưu ý chu kỳ lặp lại khoảng 10.000 năm, giá trị rung động vùng nguồn tương đương với cường độ chấn động bề mặt cấp IX theo thang MSK-64 Ngoài ra, thềm lục địa Nam Trung Bộ Việt Nam có hai vùng nguồn có rung động mạnh Cửu Long - Côn Sơn Thuận Hải - Minh Hải Ứng với chu kỳ lặp lại 9975 năm, rung động cực đại hai vùng nguồn đạt tới giá trị tương Hình Bản đồ gia tốc cực đại (PGA) lãnh thổ Việt Nam Biển Đông tính cho đá với 0,5% xác suất bị vượt vòng 50 năm Bảng liệt kê giá trị gia tốc cực đại (PGA) cường độ chấn động mặt I theo thang MSK-64 tương ứng tính số vùng nguồn chấn động tiêu biểu lãnh thổ thềm lục địa Việt Nam theo bốn chu kỳ thời gian khác 87 Nguyễn Hồng Phương, Phạm Thế Truyền Bảng Gia tốc cực đại cường độ chấn động mặt theo thang MSK-64 tính vùng nguồn chấn động lãnh thổ thềm lục địa Việt Nam N 31 16 33 20 32 19 22 35 Tên vùng nguồn Kinh tuyến 109 Điện Biên - Lai Châu Sơn La Cửu Long - Côn Sơn Rào Nậy Thuận Hải - Minh Hải Sông Hồng - Sông Chảy Sông Cả - Khe Bố Trà Bồng Bắc Biển Đông Chu kỳ thời gian (năm) 475 975 2475 9975 118/VIII 180/VIII 194/VIII 124/VIII 110/VII 98/VII 134/VIII 120/VIII 100/VII 85/VII 151/VIII 202/VIII 210/VIII 144/VIII 125/VIII 126/VIII 152/VIII 135/VIII 121/VIII 98/VII 200/VIII 235/VIII 232/VIII 178/VIII 150/VIII 150/VIII 179/VIII 160/VIII 142/VIII 118/VIII 285/IX 272/IX 264/IX 198/VIII 194/VIII 181/VIII 181/VIII 180/VIII 171/VIII 122/VIII Chú thích: Cột liệt kê mã vùng nguồn tương ứng với minh họa hình bảng Các số tử số biểu thị giá trị gia tốc cực đại tính gal, số mẫu số biểu thị cường độ chấn động mặt tương ứng KẾT LUẬN Nghiên cứu trình bày kết hiệu chỉnh tập đồ xác suất nguy hiểm động đất lãnh thổ Việt Nam khu vực Biển Đông xây dựng năm 2010 Trên sở cập nhật danh mục động đất đến năm 2014 sử dụng kết nghiên cứu kiến tạo - địa động lực khu vực Đông Nam Á, ranh giới 37 vùng nguồn chấn động toàn lãnh thổ Việt Nam khu vực Biển Đông tới kinh độ 1250E xác định Để tính tốn độ nguy hiểm động đất cho lãnh thổ Việt Nam, quy luật tắt dần chấn động Toro cộng (1997) xây dựng cho vùng có tính địa chấn ổn định lựa chọn Ngoài ra, ảnh hưởng rung động động đất Biển Đơng tính đến việc áp dụng quy luật tắt dần chấn động xây dựng cho đới hút chìm Youngs, Chiou, Silva Humphrey (1997) đề xuất cho vùng nguồn máng biển sâu Manila bờ tây Philippines Các kết trình bày dạng tập đồ biểu thị phân bố không gian đại lượng gia tốc cực đại (PGA) dự báo với xác suất bị vượt 10%, 5%, 2% 0,5% vòng 50 năm Tập đồ xác suất nguy hiểm động đất thành lập kết dự báo ngắn hạn, trung bình dài hạn độ nguy hiểm động đất lãnh thổ Việt Nam khu vực Biển Đông Ở phạm vi quốc gia khu vực, đồ cung cấp thông tin tham 88 chiếu quan trọng phục vụ công tác quy hoạch phát triển kinh tế, quốc phòng đất nước Ngoài ra, giá trị gia tốc cực đại nhận từ đồ sử dụng trực tiếp tính tốn thiết kế kháng chấn cho cơng trình xây dựng tải trọng động đất nghiên cứu liên quan đến phân vùng động đất TÀI LIỆU THAM KHẢO Nguyễn Hồng Phương, 1993 Đánh giá xác suất độ nguy hiểm động đất cho lãnh thổ Việt Nam Luận án tiến sĩ, Viện Vật lý Trái đất, Viện Hàn lâm Khoa học Liên bang Nga, Matxcơva (Tiếng Nga) Nguyen Hong Phuong, 1997 Probabilistic earthquake hazard assessment for Vietnam and adjacent regions Proceedings of the National Centre for Science and Technology of Vietnam, 9(1): 115-130 Nguyễn Hồng Phương, 2004 Bản đồ độ nguy hiểm động đất Việt Nam Biển Đơng Tạp chí Các khoa học Trái đất, 26(2): 97-111 GSHAP Final Report, 1997 Nguyễn Đình Xuyên, Lê Tử Sơn, 2009 Bản đồ phân vùng động đất lãnh thổ Việt Nam tiêu chuẩn xây dựng Việt Nam: TCXDVN 375-2006: Thiết kế cơng trình chiụ động đất Tạp chí Các khoa học Trái đất 30(4): 503-513 Tập đồ xác suất nguy hiểm động đất … Bùi Công Quế (Chủ biên), 2010 Nguy hiểm động đất sóng thần vùng ven biển Việt Nam Nxb Khoa học tự nhiên Công nghệ Chương V Tr 169-185 Cornell, C A., 1968 Engineering seismic risk analysis Bulletin of the Seismological Society of America, 58(5): 1583-1606 Esteva, L., 1968 Bases para la formulación de decisiones de diso sísmico Instituto de Ingeniería, Universidad Nacional Autónoma de México Nguyễn Hồng Phương, Phạm Thế Truyền, 2007 Xây dựng mơ hình nguồn tuyến đánh giá rủi ro động đất Việt Nam Tạp chí Khoa học Trái đất, 29(3): 228-238 10 Nguyen Hong Phuong, Bui Cong Que and Nguyen Dinh Xuyen, 2012 Investigation of earthquake tsunami sources, capable of affecting Vietnamese coast Natural hazards, 64(1): 311-327 11 Nguyen Hong Phuong, Bui Cong Que, Vu Ha Phuong and Pham The Truyen, 2014 Scenario-based tsunami hazard assessment for the coast of Vietnam from the Manila Trench source Physics of the Earth and Planetary Interiors, 236, 95-108 12 Keilis-Borok, V I., Knopoff, L., and Rotvain, I M., 1980 Bursts of aftershocks, long-term precursors of strong earthquakes Nature, 283, 259-263 13.Tran Viet Hung and Kiyomiya, O., 2012 Ground motion attenuation relationship for shallow strike-Slip earthquakes in Northern Vietnam Based on Strong motion records from Japan, Vietnam and adjacent regions Structural Eng./Earthquake Eng., 68(3): 509-525 14 Nguyen, L M., Lin, T L., Wu, Y M., Huang, B S., Chang, C H., Huang, W G., and Dinh, V T., 2012 The first peak ground motion attenuation relationships for north of Vietnam Journal of Asian Earth Sciences, 43(1): 241-253 15 McGuire, R K., 1976 FORTRAN computer program for seismic risk analysis (No 76-67) US Geological Survey 16 Ordaz, M., Aguilar, A., Arboleda, J., (?) Phần mềm CRISIS99 University of Mexico (UNAM) (TL tiếng Anh, tiếng Việt lẫn lộn) 17 Nguyen Hong Phuong, 1991 Probabilistic assessment of earthquake hazard in Vietnam based on seismotectonic regionalization Tectonophysics, 198(1): 81-93 19 Nguyen Hong Phuong, 2001 Probabilistic seismic hazard assessment along the Southeastern coast of Vietnam Natural hazards, 24(1): 53-74 22 Nguyen Hong Phuong, Pham The Truyen, 2014 Probabilistic seismic hazard assessment for the South Central Vietnam Vietnam Journal of Earth Scence (đang in) 23 Toro, G R., Abrahamson, N A., and Schneider, J F., 1997 Engineering model of strong ground motions from earthquakes in the central and eastern United States Seismological Research Letters, 68(1): 4157 24 Gusiakov, V K., 2005 Integrated tsunami database for the world ocean, version 5.15 Tsunami Laboratory, ICMMG SD RAS, Novosibirsk, Russia 25 Xiang, J., and Gao, D., 1994 The attenuation law of horizontal peak acceleration on the rock site in Yunnan area Earthquake Research in China, 8(4): 509-516 26 Youngs, R R., Chiou, S J., Silva, W J., and Humphrey, J R., 1997 Strong ground motion attenuation relationships for subduction zone earthquakes Seismological Research Letters, 68(1): 58-73 89 Nguyễn Hồng Phương, Phạm Thế Truyền PROBABILISTIC SEISMIC HAZARD MAPS OF VIETNAM AND THE EAST VIETNAM SEA Nguyen Hong Phuong, Pham The Truyen Institute of Geophysics-VAST ABSTRACT: In this paper, the probabilistic seismic hazard maps for the territory of Vietnam and the East Vietnam Sea revised from the existing 2010 ones are presented An earthquake catalog updated until 2014 and most recent seismotectonic and geodynamic information of South East Asia were used for delineation of 37 seismic source zones in the whole territory of Vietnam and the East Vietnam Sea area extended to 1250E While the Toro et al (1997) ground motion prediction equation was used for most seismic sources, the equation developed for the subduction zones by Youngs, Chiou, Silva and Humphrey (1997) was used for the Manila Trench source zone The hazard maps show distribution of the mean peak ground acceleration (PGA) with 10%, 5%, 2% and 0.5% probability of exceedance in 50 years The highest hazard areas were revealed in northwestern part of the country, in Dien Bien - Lai Chau and Son La source zones, with the maximum values of PGA for 475 - 9975 years of 180 - 272 gals; and off-shore Vietnam’s South Central coast, in the 1090Meridian Fault source zone, where the maximum values of PGA for 475 9975 years are 118 - 285 gals These PGA maps present both short - term and long - term forecasts of seismic hazard in Vietnam and the East Vietnam Sea and can be used as a reference for antiseismic design and many engineering applications Keywords: Probabilistic seismic hazard maps, peak ground acceleration, seismic source zones, ground motion prediction equation 90 ... thành lập đồ xác suất nguy hiểm động đất cho dải ven biển Việt Nam Biển Đông [6], đến năm 2014, tập thể tác giả sử dụng mơ hình NGA08 để thành lập đồ xác suất nguy hiểm động đất cho vùng Nam Trung... độ nguy hiểm động đất Ở Việt Nam, thiếu số liệu gia tốc từ trận động đất mạnh nên thời gian dài việc đánh giá nguy hiểm động đất phải sử dụng biểu thức tắt dần chấn động Tập đồ xác suất nguy hiểm. .. trình chiụ động đất Tạp chí Các khoa học Trái đất 30(4): 503-513 Tập đồ xác suất nguy hiểm động đất … Bùi Công Quế (Chủ biên), 2010 Nguy hiểm động đất sóng thần vùng ven biển Việt Nam Nxb Khoa