Tổng quan nghiên cứu

Hoạt động địa chấn tại Việt Nam thể hiện tính phân hóa sâu sắc theo không gian lãnh thổ, trong đó hơn 90% các trận động đất ghi nhận trong lịch sử tập trung chủ yếu tại khu vực Tây Bắc với trận động đất mạnh nhất đạt độ lớn M = 6,7 xảy ra vào năm 1983. Trái lại, khu vực Nam Trung Bộ và Nam Bộ thuộc khối nền Sunda tương đối ổn định với tần suất xuất hiện động đất thấp hơn nhiều. Tuy nhiên, việc xây dựng các công trình trọng điểm quốc gia như Nhà máy điện hạt nhân Ninh Thuận 2 tại xã Vĩnh Hải, huyện Ninh Hải, tỉnh Ninh Thuận đòi hỏi mức độ an toàn nghiêm ngặt ở cấp độ cao nhất. Các thảm họa địa chấn trên thế giới, tiêu biểu như sự cố tại nhà máy Kashiwazaki-Kariwa năm 2007 và Fukushima năm 2011, cho thấy dao động nền thực tế có thể vượt xa dự báo ban đầu nếu không đánh giá đầy đủ các đứt gãy kiến tạo ẩn sâu.

Mục tiêu cốt lõi của nghiên cứu là xác định chuẩn xác giá trị động đất thiết kế cực đại (MDE - Maximum Design Earthquake) và gia tốc rung động nền cực đại (PGA - Peak Ground Acceleration), làm cơ sở kỹ thuật trực tiếp cho việc tính toán kháng chấn và lựa chọn công nghệ lò phản ứng phù hợp. Phạm vi không gian của nghiên cứu bao quát toàn diện các cấu trúc địa chất trong bán kính 300 km tính từ tâm địa điểm xây dựng nhà máy, kết hợp phân tích chi tiết vùng nguồn chấn tiêu trong bán kính 40 km và 10 km lân cận. Nghiên cứu xử lý chuỗi số liệu địa chấn lịch sử kéo dài từ năm 1067 đến giai đoạn hiện đại, đặc biệt chú trọng các số liệu ghi đo bằng mạng lưới trạm quan trắc quốc gia sau năm 1975. Ý nghĩa của công trình thể hiện ở việc đảm bảo an toàn tuyệt đối cho chu kỳ vận hành 60 năm của nhà máy trước các kịch bản động đất có chu kỳ lặp lại lên tới 10.000 năm, đồng thời tối ưu hóa chi phí đầu tư xây dựng hạ tầng kỹ thuật hạt nhân cho Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết địa động lực học mảng, cơ học đứt gãy và lý thuyết xác suất thống kê các biến cố cực trị địa chấn. Ba khái niệm chuyên ngành then chốt được chuẩn hóa bao gồm:

  1. Động đất cực đại có thể xảy ra (MCE - Maximum Credible Earthquake) đại diện cho trận động đất lớn nhất mà một đới đứt gãy có khả năng phát sinh xét về mặt vật lý kiến tạo.
  2. Động đất thiết kế cực đại (MDE) là mức dao động nền chuẩn được quy định để kết cấu công trình chịu biến dạng mà không sụp đổ hay phát tán phóng xạ.
  3. Gia tốc rung động nền cực đại (PGA) phản ánh biên độ dao động lớn nhất của nền đá gốc tại địa điểm xây dựng.

Mô hình quan hệ Gutenberg-Richter dạng $\log N = a - bM$ được sử dụng để mô tả tần suất xuất hiện động đất theo độ lớn tại từng vùng nguồn. Kết hợp với đó là lý thuyết phân bố cực trị Gumbel loại I và loại III, cho phép ngoại suy xác suất xuất hiện và chu kỳ lặp lại $T(M)$ của các trận động đất mạnh từ chuỗi số liệu rời rạc. Về phương diện truyền sóng, nghiên cứu tích hợp các mô hình dự báo chuyển động đất (GMPEs - Ground Motion Prediction Equations) hiện đại, nổi bật là phương trình thực nghiệm của Campbell & Bozorgnia (2008) và Abrahamson & Silva (2008), giúp mô phỏng chính xác sự suy giảm biên độ gia tốc sóng địa chấn theo khoảng cách từ chấn tiêu tới vị trí đặt lò phản ứng.

Phương pháp nghiên cứu

Nguồn dữ liệu của luận văn được tổng hợp từ cơ sở dữ liệu động đất quốc gia với 10 trận chấn động tiêu biểu có độ lớn từ M = 5,0 đến M = 6,7, kết hợp tài liệu địa chấn nông phân giải cao và số liệu địa chấn dầu khí trên thềm lục địa. Cỡ mẫu nghiên cứu bao gồm toàn bộ 6 hệ thống đới đứt gãy lớn trong bán kính 300 km (như đới Nha Trang - Tánh Linh, Tuy Hòa - Củ Chi, Thuận Hải - Minh Hải) và 3 đới đứt gãy cục bộ trong bán kính 40 km (đứt gãy Đông Hòn Gió, Đông núi Bàu, sông Rừa). Phương pháp chọn mẫu cấu trúc địa chất được thực hiện theo tiêu chí đứt gãy hoạt động trong kỷ Đệ tứ (khoảng 10.000 năm trở lại đây) theo đúng khuyến cáo an toàn của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA).

Lý do lựa chọn phương pháp tất định (DSHA - Deterministic Seismic Hazard Analysis) làm phương pháp tiếp cận chủ đạo bắt nguồn từ thực tế chuỗi quan trắc địa chấn tại miền Nam Việt Nam còn ngắn, mạng lưới trạm đo dày đặc chỉ mới hình thành sau năm 1975. Trong điều kiện thiếu hụt chuỗi dữ liệu hàng ngàn năm, phương pháp xác suất (PSHA) dễ dẫn đến đánh giá thấp rủi ro thực tế. Phương pháp tất định cho phép thiết lập kịch bản sự cố xấu nhất (worst-case scenario) bằng cách đặt chấn tiêu động đất cực đại $M_{\max}$ tại vị trí đứt gãy gần nhà máy nhất. Quá trình tính toán mô phỏng được thực hiện liên tục trong mốc thời gian từ năm 2011 đến 2015, kết hợp kiểm chứng chéo với các kết quả khảo sát độc lập của các tổ chức tư vấn quốc tế.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phân tích cấu trúc địa chất kiến tạo khẳng định khu vực xây dựng Nhà máy điện hạt nhân Ninh Thuận 2 nằm trọn trong khối nền ổn định Sunda với độ lớn động đất cực đại giới hạn ở mức $M_{\max} \le 7,0$. Trong phạm vi bán kính 300 km, nghiên cứu đã phân vùng chi tiết 6 hệ thống đới đứt gãy chính có khả năng phát sinh chấn động ảnh hưởng đến địa điểm:

  1. Đới đứt gãy Nha Trang - Tánh Linh (cách nhà máy 65 km về phía Tây Bắc) gồm 5 phân đoạn, trong đó phân đoạn 1 và 2 hoạt động mạnh nhất.
  2. Đới cắt Tuy Hòa - Củ Chi (cách 75 km về phía Đông Bắc) gồm 4 phân đoạn.
  3. Đới đứt gãy Thuận Hải - Minh Hải trên thềm lục địa (cách 40 km về phía Đông Nam).
  4. Hệ thống đứt gãy Phan Thiết - Đèo Bảo Lộc (cách 100 km về phía Tây Nam).

Trong phạm vi bán kính 40 km, nghiên cứu phát hiện đứt gãy phía Đông núi Hòn Gió mang dấu hiệu chuyển động từ thế Pleistocen muộn, được xác định là nguồn sinh chấn kịch bản nguy hiểm nhất gây ảnh hưởng trực tiếp tới công trình. Áp dụng mô hình tắt dần Campbell & Bozorgnia (2008), gia tốc rung động nền do đứt gãy núi Bàu tạo ra tại địa điểm đạt 171 cm/s² (tương đương khoảng 0,174g). Khi mô phỏng kịch bản chấn động cực đại từ đứt gãy Đông Hòn Gió bằng phương pháp mô hình đứt gãy lai (Irikura - Miyake), đỉnh gia tốc rung động nền tính toán đạt 410 cm/s² (khoảng 0,418g). Để đảm bảo biên độ an toàn tuyệt đối, gia tốc rung động nền cực đại trong thiết kế cơ sở (SS) theo phương ngang được xác định ở mức 450 cm/s² (tương đương 0,459g) và theo phương thẳng đứng đạt 261 cm/s² (khoảng 0,266g).

Thảo luận kết quả

Khi đối chiếu giữa các phương pháp, tồn tại sự chênh lệch số liệu đáng kể giữa phương pháp xác suất và phương pháp tất định. Cụ thể, các nghiên cứu xác suất trước đây của các nhà khoa học trong nước chỉ ra rằng với chu kỳ lặp lại 1.000 năm, giá trị PGA chỉ đạt 0,029g (tương đương cấp V theo thang chấn cấp MSK-64), và với chu kỳ 10.000 năm, PGA đạt 0,063g (cấp VII MSK-64). Giá trị này thấp hơn nhiều so với kết quả 450 cm/s² (khoảng 0,459g) thu được từ phương pháp tất định do các chuyên gia quốc tế và luận văn thực hiện. Nguyên nhân chính là do phương pháp xác suất phân tán rủi ro theo chuỗi thời gian thống kê trên diện rộng, trong khi phương pháp tất định cô lập trận động đất giả định cực đại ngay trên đới đứt gãy Đông Hòn Gió tại khoảng cách gần nhất.

Toàn bộ các tham số động học này có thể được biểu diễn một cách trực quan thông qua đồ thị đường cong phổ phản ứng gia tốc (Spectral Acceleration Response Spectrum) ở các chu kỳ dao động từ 0,01 giây đến 10 giây. Đồ thị phổ cho thấy vùng năng lượng tập trung mạnh nhất rơi vào dải chu kỳ ngắn từ 0,05 giây đến 0,2 giây với đỉnh phổ gia tốc vượt ngưỡng 1.000 Gal. Bảng ma trận so sánh các tham số nguồn chấn cũng minh chứng rõ ràng rằng điều kiện địa chất nền đá gốc cứng tại xã Vĩnh Hải đóng vai trò triệt tiêu hiện tượng khuếch đại sóng địa chấn mặt, ngăn chặn nguy cơ cộng hưởng phá hủy kết cấu bê tông cốt thép của tòa nhà lò phản ứng.

Đề xuất và khuyến nghị

Dựa trên các kết quả tính toán động đất kịch bản, 4 nhóm giải pháp kỹ thuật và quản lý được khuyến nghị triển khai:

  1. Thiết lập chuẩn an toàn kháng chấn cơ sở: Cơ quan quản lý dự án và các đơn vị thiết kế phải áp dụng giá trị gia tốc rung động nền cực đại theo phương ngang tối thiểu 450 cm/s² (0,459g) và phương đứng 261 cm/s² (0,266g) làm tham số đầu vào bắt buộc cho thiết kế kháng chấn của các công trình nhóm cấp an toàn hạt nhân loại 1.
  2. Hiện đại hóa hệ thống quan trắc địa chấn chuyên dụng: Tập đoàn Điện lực Việt Nam phối hợp với Viện Vật lý Địa cầu lắp đặt mạng lưới quan trắc vi địa chấn cục bộ gồm tối thiểu 6 trạm kỹ thuật số 3 thành phần băng rộng trong bán kính 40 km quanh vị trí nhà máy, hoàn thành đo đạc liên tục tối thiểu 3 năm trước khi bắt đầu đổ bê tông móng lò.
  3. Thăm dò địa vật lý tầng sâu định kỳ: Thực hiện công tác khảo sát địa chấn nông phân giải cao kết hợp đo đạc vận tốc sóng cắt tại độ sâu 30 m ($V_{s30}$) với chu kỳ 2 năm một lần nhằm phát hiện sớm các chuyển dịch kiến tạo vi mô dọc theo đứt gãy Đông Hòn Gió và đứt gãy Thuận Hải - Minh Hải.
  4. Xây dựng quy chuẩn quốc gia về đánh giá địa điểm hạt nhân: Cục An toàn bức xạ và hạt nhân chủ trì ban hành văn bản hướng dẫn kỹ thuật kết hợp bắt buộc giữa phương pháp tất định (DSHA) và phương pháp xác suất (PSHA) trong thẩm định Báo cáo phân tích an toàn địa điểm (SAR), áp dụng đồng bộ từ giai đoạn 2016-2025.

Đối tượng nên tham khảo luận văn

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cho 4 nhóm đối tượng:

  1. Kỹ sư kết cấu và thiết kế công trình hạt nhân, thủy điện: Áp dụng trực tiếp bộ thông số gia tốc nền cực đại (PGA) và mô hình suy giảm sóng Campbell & Bozorgnia (2008) để tính toán tải trọng động lực học cho các bể chứa nhiên liệu, đập ngăn nước và tòa nhà lò phản ứng.
  2. Nhà nghiên cứu và giảng viên chuyên ngành Vật lý Địa cầu, Địa kỹ thuật: Sử dụng khung lý thuyết phân tích cực trị Gumbel, mô hình đứt gãy lai Irikura - Miyake và phương pháp giải bài toán thiếu hụt chuỗi dữ liệu địa chấn lịch sử làm tài liệu giảng dạy và nghiên cứu mở rộng.
  3. Cơ quan quản lý nhà nước về an toàn hạt nhân và môi trường: Sử dụng dữ liệu vùng nguồn chấn động trong bán kính 300 km làm căn cứ khoa học độc lập để thẩm tra, phê duyệt Báo cáo đánh giá tác động môi trường và Báo cáo an toàn địa điểm.
  4. Các đơn vị tư vấn khảo sát địa chất công trình ven biển Nam Trung Bộ: Kế thừa bản đồ phân bố các đới đứt gãy tại Ninh Thuận, Khánh Hòa, Bình Thuận nhằm phục vụ quy hoạch các dự án hạ tầng cảng biển nước sâu và khu công nghiệp năng lượng.

Câu hỏi thường gặp

1. Động đất thiết kế cực đại (MDE) khác gì so với động đất cực đại có thể xảy ra (MCE)? Động đất cực đại có thể xảy ra (MCE) là giá trị vật lý biểu thị trận động đất lớn nhất có thể phát sinh từ cấu trúc địa chất của vùng nguồn (khoảng M = 7,0 tại khối Sunda). Trong khi đó, Động đất thiết kế cực đại (MDE) là giá trị kỹ thuật quy chuẩn (như mức gia tốc ngang 450 cm/s²), được chọn làm mốc tải trọng để thiết kế các kết cấu hạt nhân duy trì trạng thái an toàn tuyệt đối khi xảy ra rung chấn MCE.

2. Vì sao luận văn ưu tiên phương pháp tất định (DSHA) hơn phương pháp xác suất (PSHA)? Khu vực miền Nam Việt Nam có lịch sử ghi nhận động đất ngắn, chuỗi số liệu trạm đo chuẩn xác chỉ bắt đầu từ sau năm 1975. Do thiếu dữ liệu chu kỳ lặp hàng ngàn năm, phương pháp xác suất cho kết quả gia tốc nền tương đối thấp (0,029g - 0,063g). Phương pháp tất định khắc phục điểm yếu này bằng cách giả định kịch bản xấu nhất xảy ra trên đứt gãy Đông Hòn Gió, mang lại hệ số an toàn cao hơn cho công trình hạt nhân.

3. Khu vực Ninh Thuận có nguy cơ xảy ra siêu động đất trên 8,0 độ Richter không? Theo các khảo sát địa động lực học của Cơ quan Khảo sát Địa chất Hoa Kỳ (USGS) và Viện Hàn lâm Khoa học và Công nghệ Việt Nam, Ninh Thuận nằm trọn trong mảng Sunda ổn định, không có các đới hút chìm sâu dưới lòng đất. Chiều dài và bề dày các đới đứt gãy như Nha Trang - Tánh Linh hay Tuy Hòa chỉ có khả năng tích lũy năng lượng để phát sinh động đất cực đại đạt mức $M_{\max} \le 7,0$.

4. Mô hình tắt dần của Campbell & Bozorgnia (2008) có ưu điểm gì trong tính toán rung động nền? Mô hình Campbell & Bozorgnia (2008) là phương trình thực nghiệm tiên tiến thuộc chương trình NGA-West, xét đến loại cơ cấu trượt của đứt gãy (trượt bằng, thuận, nghịch), độ sâu chấn tiêu, khoảng cách chấn đoạn ($R_{\text{seis}}$) và vận tốc truyền sóng cắt tầng nông ($V_{s30}$). Mô hình này cho phép mô phỏng chính xác sự suy giảm năng lượng sóng địa chấn trên nền đá gốc cứng đặc thù tại xã Vĩnh Hải.

5. Phạm vi bán kính khảo sát đứt gãy bắt buộc đối với nhà máy điện hạt nhân là bao nhiêu? Theo tiêu chuẩn an toàn của IAEA và quy chuẩn xây dựng hạt nhân quốc tế, phạm vi khảo sát địa chấn bắt buộc phải phân cấp thành ba vùng: bán kính 300 km để đánh giá các đới đứt gãy khu vực quy mô lớn, bán kính 40 km để xác định các đứt gãy gần có nguy cơ kích hoạt, và bán kính 10 km lân cận để kiểm tra hiện tượng đứt gãy mặt đất tại vị trí móng công trình.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở khoa học về bối cảnh địa chất kiến tạo và chế độ địa chấn tại khu vực Nam Trung Bộ, xác nhận tính ổn định tương đối của khối nền Sunda với ngưỡng động đất cực đại $M_{\max} \le 7,0$.
  • Xác định và phân loại chi tiết 6 hệ thống đới đứt gãy lớn trong bán kính 300 km và nhận diện đứt gãy Đông Hòn Gió trong phạm vi bán kính 40 km là nguồn sinh chấn kịch bản nguy hiểm nhất đối với công trình.
  • Chứng minh tính ưu việt của phương pháp tất định (DSHA) kết hợp mô hình tắt dần Campbell & Bozorgnia (2008) trong việc giải quyết bài toán thiếu hụt chuỗi dữ liệu quan trắc địa chấn lịch sử tại miền Nam.
  • Thiết lập thành công giá trị gia tốc rung động nền cực đại phục vụ thiết kế cơ sở (SS) với gia tốc ngang 450 cm/s² (0,459g) và gia tốc đứng 261 cm/s² (0,266g), đảm bảo an toàn tuyệt đối trước các chấn động cấp VII - VIII MSK-64.
  • Đóng góp luận cứ kỹ thuật chuẩn xác phục vụ trực tiếp công tác thẩm định an toàn hạt nhân và lựa chọn công nghệ lò phản ứng cho dự án Nhà máy điện hạt nhân Ninh Thuận 2.

Đóng góp lớn nhất của công trình là cung cấp một khung tính toán địa chấn kịch bản hoàn chỉnh, giải quyết hài hòa giữa yêu cầu an toàn hạt nhân tối đa và bài toán tối ưu chi phí xây dựng. Để ứng dụng hiệu quả kết quả nghiên cứu trong các giai đoạn tiếp theo, các cơ quan hữu quan cần tiếp tục mở rộng mạng trạm quan trắc địa chấn vi mô trước năm 2025 và cập nhật liên tục mô hình phản ứng nền đất. Hãy tham khảo chi tiết toàn văn tài liệu để nắm vững phương pháp mô hình hóa nguồn chấn và khai thác bộ tham số kỹ thuật phục vụ thiết kế các công trình công nghiệp trọng điểm quốc gia.