Tổng quan về luận án

Nghiên cứu phông phóng xạ tự nhiên và bức xạ môi trường là trụ cột nền tảng của vật lý hạt nhân ứng dụng và an toàn bức xạ y sinh. Theo báo cáo của Ủy ban Khoa học Liên Hiệp Quốc về Tác động của Bức xạ Nguyên tử (UNSCEAR) và Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA), các đồng vị phóng xạ nguyên thủy tồn tại từ khi hình thành Trái Đất—bao gồm Kali ($^{40}\text{K}$, chu kỳ bán rã $T_{1/2} = 1,26 \times 10^9$ năm), chuỗi Thori ($^{232}\text{Th}$, $T_{1/2} = 1,40 \times 10^{10}$ năm), và chuỗi Uran ($^{238}\text{U}$, $T_{1/2} = 4,47 \times 10^9$ năm)—là nguồn đóng góp tới hơn 80% tổng liều bức xạ hiệu dụng tự nhiên mà con người tiếp nhận hàng năm. Mặc dù các mô hình dịch tễ học về rủi ro bức xạ ion hóa liều thấp vẫn còn nhiều tranh luận học thuật, nguyên lý không ngưỡng tuyến tính (LNT) được chấp nhận rộng rãi xác nhận rằng khi liều hiệu dụng tích lũy tăng thêm $1\text{ Sv}$, xác suất phát sinh các bệnh lý ác tính, đặc biệt là ung thư biểu mô và ung thư bạch cầu, tăng thêm $5,5%$.

+--------------------------------------------------------------------------------------------------------+
|                                    BẢN ĐỒ PHÔNG PHÓNG XẠ TỰ NHIÊN                                      |
|                                                                                                        |
|  [Chuỗi 238U / 226Ra] --------+                                                                        |
|  [Chuỗi 232Th] ---------------> [Đất đá bề mặt] ---> [Bức xạ Gamma mặt đất] ---> [Chiếu ngoài: Raeq,  |
|  [Đồng vị đơn 40K] -----------+         |                                         Dair, AEDO, Hex]     |
|                                         v                                                              |
|                                  [Khí 222Rn / 220Rn] -> [Đường hô hấp] ---------> [Chiếu trong: Hin]   |
+--------------------------------------------------------------------------------------------------------+

Tại Cộng hòa Dân chủ Nhân dân Lào (CHDCND Lào), bức tranh dữ liệu về phông phóng xạ tự nhiên trong vỏ Trái Đất tồn tại khoảng trống nghiên cứu (research gap) nghiêm trọng. Trước công trình này, tài liệu quốc tế và bản địa chỉ ghi nhận duy nhất khảo sát sơ bộ của Sonexay và Lê Hồng Khiêm (2018) trên $10$ mẫu đất hạn hẹp tại huyện Thoulakhom, tỉnh Viêng Chăn dùng cho sản xuất vật liệu xây dựng. Chưa từng có một công trình nghiên cứu hệ thống nào đánh giá toàn diện hoạt độ phóng xạ tự nhiên, đặc tính địa chất và các chỉ số nguy hiểm bức xạ tại khu vực miền Trung nước Lào—nơi có mật độ dân số cao nhất cả nước, địa hình đan xen giữa đồng bằng và đồi núi phức tạp, đồng thời sở hữu các điểm dị thường khoáng sản và mỏ vàng lớn đã đưa vào khai thác quy mô công nghiệp (mỏ vàng Khamkeut ở tỉnh Bolikhamxay, mỏ vàng Sepone và Vilabuly ở tỉnh Savannakhet).

Luận án tiến sĩ của NCS. Somsavath Leuangtakoun với đề tài "Nghiên cứu ứng dụng phương pháp phổ gamma tự nhiên đánh giá tính phóng xạ đất đá trên bề mặt tại miền trung nước Lào", thực hiện tại Bộ môn Vật lý hạt nhân, Khoa Vật lý, Trường Đại học Khoa học Tự nhiên – ĐHQGHN dưới sự hướng dẫn của PGS. Bùi Văn Loát và TS. Phan Việt Cương, đã giải quyết triệt để khoảng trống khoa học này. Luận án giải quyết $4$ câu hỏi nghiên cứu và giả thuyết khoa học trọng tâm:

  1. RQ1 & H1: Làm thế nào để chuẩn hóa và nâng cao độ chính xác phép đo phổ kế gamma nhấp nháy $\text{NaI(Tl)}$ có độ phân giải năng lượng thấp khi phân tích đồng thời các đỉnh năng lượng chồng chập của $^{40}\text{K}$, $^{226}\text{Ra}$ và $^{232}\text{Th}$? Giả thuyết cho rằng việc xây dựng phương pháp ma trận phân tích phổ đa thông lượng kết hợp cửa sổ năng lượng đặc trưng sẽ khử được can nhiễu tán xạ Compton liên tục.
  2. RQ2 & H2: Sự sai khác về mật độ khối và thành phần chất nền (matrix effect) giữa mẫu chuẩn và mẫu phân tích môi trường thực tế gây ra sai số tự hấp thụ bức xạ gamma như thế nào? Giả thuyết khẳng định kỹ thuật xây dựng đường cong chuẩn nội hiệu suất ghi (internal standard efficiency curve) dựa trên các đỉnh gamma tự nhiên của $^{214}\text{Bi}$ kết hợp hiệu suất tuyệt đối tại vạch $1460,82\text{ keV}$ sẽ loại trừ hoàn toàn sai số hình học và hiệu ứng tự hấp thụ chất nền.
  3. RQ3 & H3: Trạng thái cân bằng phóng xạ thế kỷ (secular equilibrium) trong chuỗi $^{238}\text{U} \to {}^{226}\text{Ra} \to {}^{222}\text{Rn} \to {}^{214}\text{Pb} \to {}^{214}\text{Bi}$ có được bảo toàn trong các nhóm đất đá tự nhiên tại miền Trung Lào sau quy trình lưu mẫu kín $30$ ngày hay không?
  4. RQ4 & H4: Mức độ nguy hiểm bức xạ đối với sức khỏe cộng đồng dân cư sinh sống trên $30$ huyện thuộc $3$ tỉnh miền Trung Lào có vượt qua các ngưỡng khuyến cáo an toàn của UNSCEAR và IAEA hay không?

Khung lý thuyết của luận án tích hợp chặt chẽ giữa lý thuyết tương tác bức xạ gamma với vật chất (hấp thụ quang điện, tán xạ Compton, tạo cặp), động học phân rã phóng xạ chuỗi liên tiếp Bateman, lý thuyết cân bằng phóng xạ bền, và các mô hình đo liều bức xạ học (radiation dosimetry). Quy mô nghiên cứu bao quát $202$ mẫu đất đá bề mặt được thu thập theo mạng lưới chuẩn hóa trên toàn bộ $30$ đơn vị hành chính cấp huyện thuộc $3$ tỉnh miền Trung nước Lào: Bolikhamxay ($6$ huyện: $1$ đồng bằng, $5$ miền núi), Khammouane ($9$ huyện: $6$ đồng bằng, $3$ miền núi) và Savannakhet ($15$ huyện: $13$ đồng bằng, $2$ miền núi).

Literature Review và Positioning

Nghiên cứu bức xạ môi trường trên thế giới hình thành hai luồng quan điểm và trường phái lý thuyết lớn:

+-------------------------------------------------------------------------------------------------------+
|                                    TRANH LUẬN LÝ THUYẾT & ĐỊNH VỊ                                     |
|                                                                                                       |
|  [Trường phái Cân bằng thế kỷ Lý tưởng] <====== TRANH BIỆN ======> [Trường phái Bất cân bằng Môi trường]|
|  (Coi tỷ số hoạt độ 238U/226Ra = 1;                              (Tồn tại mất cân bằng địa hóa do     |
|   áp dụng hệ số chuyển đổi cố định)                                phong hóa, rửa trôi, thoát khí 222Rn)|
|                                           |                                                           |
|                                           v                                                           |
|                 [Định vị của Luận án: Phương pháp Chuẩn nội Hiệu suất ghi]                            |
|                 - Đo trực tiếp tỷ số 214Pb/214Bi và 238U/226Ra                                        |
|                 - Khử sai số hiệu ứng tự hấp thụ chất nền (Matrix effect)                             |
+-------------------------------------------------------------------------------------------------------+

Trường phái thứ nhất—đại diện bởi các nghiên cứu đo liều kinh điển của UNSCEAR (2000)—thường giả định trạng thái cân bằng phóng xạ thế kỷ lý tưởng được thiết lập hoàn hảo giữa hạt nhân mẹ đầu dãy $^{238}\text{U}$ và các đồng vị cháu con phát xạ gamma mạnh như $^{214}\text{Pb}$ và $^{214}\text{Bi}$. Do đó, nhiều tác giả quốc tế như Othman & Yassine (1996) tại Syria, Al-Amir et al. (2009) tại Jordan, hay Saleh et al. tại Malaysia đã ngoại suy trực tiếp hoạt độ $^{238}\text{U}$ thông qua hoạt độ của $^{226}\text{Ra}$ hoặc các vạch gamma của $^{214}\text{Bi}$.

Ngược lại, trường phái thứ hai—được nhấn mạnh bởi Malanca et al. (1996) tại Brazil, Taskin et al. (2009) tại Thổ Nhĩ Kỳ, và Ngô Quang Huy cùng cộng sự (2018) tại Việt Nam—chứng minh rằng trong môi trường đất đá bề mặt chịu tác động mạnh của quá trình phong hóa nhiệt đới, xói mòn, rửa trôi thủy văn và đặc tính thoát khí của Radon ($^{222}\text{Rn}$ và $^{220}\text{Rn}$), trạng thái cân bằng phóng xạ thường xuyên bị phá vỡ cục bộ. Nếu không có quy trình lưu giữ kín mẫu đủ thời gian (ít nhất $30$ ngày) và không kiểm soát tỷ số hoạt độ phóng xạ giữa $^{238}\text{U}$ và $^{226}\text{Ra}$, sai số định lượng hoạt độ riêng và liều hấp thụ có thể vượt quá $35%$.

Quốc gia / Vùng lãnh thổ Hoạt độ riêng $^{40}\text{K}$ ($\text{Bq/kg}$) Hoạt độ riêng $^{226}\text{Ra}$ ($\text{Bq/kg}$) Hoạt độ riêng $^{232}\text{Th}$ ($\text{Bq/kg}$) Suất liều không khí $D_{air}$ ($\text{nGy/h}$) Tài liệu tham khảo
Miền Trung Lào $272,0 \pm 15,2$ $37,8 \pm 2,4$ $54,2 \pm 3,1$ $62,5 \pm 4,2$ Luận án (2020)
Việt Nam (Trung bình 63 tỉnh) $292,6$ $28,6$ $50,7$ $54,5$ Ngô Quang Huy et al. (2018)
Thái Lan (Toàn quốc) $213,0$ $68,0$ $51,0$ $77,0$ UNSCEAR (2000), Prasit et al.
Malaysia (Toàn quốc) $310,0$ $66,0$ $82,0$ $92,0$ UNSCEAR (2000), Saleh et al.
Rio de Janeiro (Brazil) $114,0$ $32,0$ $74,0$ ($^{228}\text{Ra}$) $68,0$ Malanca et al. (1996)
Sakarya (Thổ Nhĩ Kỳ) $371,0 \pm 26,0$ $23,2 \pm 2,5$ $21,0 \pm 3,4$ $42,1$ Taskin et al. (2009)
Giá trị trung bình Thế giới $400,0$ $35,0$ $30,0$ $59,0$ UNSCEAR (2000)

Đặt trong bức tranh so sánh khu vực Đông Nam Á và quốc tế, luận án đã định vị chính xác vị trí học thuật của vùng địa chất miền Trung Lào: hoạt độ $^{232}\text{Th}$ ($54,2\text{ Bq/kg}$) và $^{226}\text{Ra}$ ($37,8\text{ Bq/kg}$) cao hơn mức trung bình toàn cầu ($30\text{ Bq/kg}$ và $35\text{ Bq/kg}$), tương đồng với các cấu trúc địa chất magma giàu khoáng sản tại Việt Nam, Thái Lan và Malaysia, trong khi hoạt độ $^{40}\text{K}$ ($272,0\text{ Bq/kg}$) thấp hơn trung bình thế giới ($400\text{ Bq/kg}$).

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và làm sâu sắc thêm các định luật vật lý hạt nhân nền tảng và mô hình định lượng bức xạ trong phân tích địa chất phóng xạ môi trường:

+--------------------------------------------------------------------------------------------------------+
|                                      KHUNG LÝ THUYẾT VÀ PHÂN TÍCH                                      |
|                                                                                                        |
|  [Quy luật Bateman & Phân rã liên tiếp] ===> [Trạng thái Cân bằng bền Thế kỷ: λ1N1 = λ2N2 = ... = λnNn]|
|                                                              |                                         |
|                                                              v                                         |
|  [Hàm phản ứng Detector: N(V) = ∫Φ(E)G(E,V)dE] ==> [Ma trận Deconvolution NaI(Tl)]                     |
|                                                              |                                         |
|                                                              v                                         |
|  [Tương tác Gamma: Quang điện + Compton + Tạo cặp] => [Đường cong chuẩn nội FEPE đỉnh 214Bi]           |
+--------------------------------------------------------------------------------------------------------+
  1. Lý thuyết phân rã phóng xạ chuỗi liên tiếp và phương trình Bateman: Nghiên cứu vận dụng hệ phương trình vi phân Bateman để mô tả biến thiên mật độ hạt nhân trong chuỗi $4n+2$ ($^{238}\text{U}$) và $4n$ ($^{232}\text{Th}$): $$\frac{dN_i(t)}{dt} = \lambda_{i-1}N_{i-1}(t) - \lambda_i N_i(t)$$ Khi hạt nhân mẹ có chu kỳ bán rã vượt trội ($T_1 \gg T_2, T_3, \dots, T_n$), hệ thống đạt trạng thái cân bằng bền thế kỷ (secular equilibrium): $$\lambda_1 N_1 = \lambda_2 N_2 = \dots = \lambda_n N_n \iff A_1 = A_2 = \dots = A_n$$ Chứng minh thực nghiệm của luận án khẳng định điều kiện lưu mẫu $30$ ngày trong hộp kín khí hoàn toàn thỏa mãn việc tái lập cân bằng giữa $^{226}\text{Ra}$ và các đồng vị cháu con phát xạ gamma $^{214}\text{Pb}$ ($295,22\text{ keV}, 351,93\text{ keV}$) và $^{214}\text{Bi}$ ($609,31\text{ keV}, 1120,3\text{ keV}, 1764,5\text{ keV}$).
  2. Lý thuyết hàm phản ứng và phân rã tín hiệu của đầu dò bức xạ: Thiết lập biểu thức giải tích liên hệ giữa phổ năng lượng photon tới $\Phi(E)$ và phổ biên độ xung đầu ra $N(V)$ thông qua hàm phản ứng $G(E,V)$: $$N(V) = \int \Phi(E) G(E,V) dE$$
  3. Mô hình định lượng liều lượng học bức xạ môi trường: Chuẩn hóa hệ thống phương trình tính toán các đại lượng đánh giá rủi ro bức xạ:
    • Hoạt độ Radi tương đương ($Ra_{eq}$): $$Ra_{eq} = A_{Ra} + 1,43 A_{Th} + 0,077 A_K \le 370\text{ Bq/kg}$$
    • Suất liều hấp thụ trong không khí ở độ cao $1\text{ m}$ ($D_{air}$): $$D_{air} = 0,0414 A_K + 0,461 A_{Ra} + 0,623 A_{Th}\text{ (nGy/h)}$$
    • Liều hiệu dụng chiếu ngoài hàng năm ($AEDO$): $$AEDO = D_{air} \times 8760\text{ h/năm} \times 0,2 \times 0,7\text{ Sv/Gy} \times 10^{-6}\text{ (mSv/năm)}$$
    • Chỉ số nguy hiểm chiếu ngoài ($H_{ex}$) và chiếu trong ($H_{in}$): $$H_{ex} = \frac{A_K}{4810} + \frac{A_{Ra}}{370} + \frac{A_{Th}}{259} \le 1$$ $$H_{in} = \frac{A_K}{4810} + \frac{A_{Ra}}{185} + \frac{A_{Th}}{259} \le 1$$

Khung phân tích độc đáo

Luận án đã tích hợp đồng thời $3$ trụ cột phương pháp luận:

+--------------------------------------------------------------------------------------------------------+
|                                     KIẾN TRÚC PHƯƠNG PHÁP NGHIÊN CỨU                                   |
|                                                                                                        |
|  [Hệ phổ kế NaI(Tl) 4"x4" + Buồng chì A340]        [Hệ phổ kế HPGe siêu phân giải (FWHM 1.8 keV)]      |
|                        |                                                 |                             |
|                        v                                                 v                             |
|       [Phương pháp Ma trận 3 Cửa sổ]                 [Kỹ thuật Đường cong Chuẩn nội Hiệu suất ghi]     |
|                        |                                                 |                             |
|                        +-----------------------+-------------------------+                             |
|                                                v                                                       |
|                                 [ĐỐI CHỨNG VÀ HIỆU CHỈNH CHÉO]                                         |
|                                 - Mẫu chuẩn quốc tế IAEA-375, RGU-1                                    |
|                                 - Hiệu chỉnh tự hấp thụ đỉnh 1460.82 keV                               |
|                                 - Đánh giá cân bằng 238U/226Ra                                         |
+--------------------------------------------------------------------------------------------------------+
  • Phương pháp ma trận phân giải phổ trên hệ $\text{NaI(Tl)}$: Thiết lập hệ phương trình ma trận tuyến tính $3 \times 3$ dựa trên tốc độ đếm tích hợp trong $3$ cửa sổ năng lượng đặc trưng: cửa sổ $^{40}\text{K}$ ($1360 - 1560\text{ keV}$ bao quanh đỉnh $1460,82\text{ keV}$), cửa sổ $^{226}\text{Ra}$ ($1660 - 1860\text{ keV}$ bao quanh đỉnh $^{214}\text{Bi}\text{ }1764,5\text{ keV}$), và cửa sổ $^{232}\text{Th}$ ($2500 - 2700\text{ keV}$ bao quanh đỉnh $^{208}\text{Tl}\text{ }2614,5\text{ keV}$). Các hệ số ma trận $a_{ik}$ được xác định từ các mẫu chuẩn đơn thành phần IAEA RGU-1, RGTh-1, RGK-1.
  • Kỹ thuật đường cong chuẩn nội hiệu suất ghi trên hệ bán dẫn HPGe: Sử dụng chính các vạch gamma đơn năng có cường độ phát xạ cao của $^{214}\text{Bi}$ phân bố đều từ dải năng lượng thấp đến cao ($609,31\text{ keV}, 768,36\text{ keV}, 934,06\text{ keV}, 1120,3\text{ keV}, 1238,1\text{ keV}, 1764,5\text{ keV}$) để xây dựng đường cong hiệu suất ghi tương đối thực tế của từng mẫu đo.
  • Kỹ thuật hiệu chuẩn kết hợp tuyệt đối - tương đối tại đỉnh $1460,82\text{ keV}$: Ghép nối đường cong hiệu suất tương đối với giá trị hiệu suất ghi tuyệt đối tại vạch $1460,82\text{ keV}$ của $^{40}\text{K}$ để khử hoàn toàn ảnh hưởng của sự sai khác mật độ mẫu từ $0,8\text{ g/cm}^3$ đến $1,6\text{ g/cm}^3$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng chuẩn tắc (positivism) với thiết kế thực nghiệm định lượng đa tầng, đối chứng chéo liên phòng thí nghiệm. Kế hoạch chọn mẫu được thiết kế đại diện tối ưu theo cấu trúc địa mạo và hành chính của $3$ tỉnh miền Trung Lào:

+--------------------------------------------------------------------------------------------------------+
|                                    QUY TRÌNH THỰC NGHIỆM ĐA TẦNG                                       |
|                                                                                                        |
|  [202 Điểm lấy mẫu tại 30 Huyện] --> [Tầng đất 0-30 cm / Bỏ tạp chất mặt]                             |
|                                                     |                                                  |
|                                                     v                                                  |
|  [Xử lý nhiệt 105-110°C / Nghiền mịn / Rây 0.2-1mm] -> [Đóng hộp kín 30 ngày (Đạt cân bằng thế kỷ)]    |
|                                                     |                                                  |
|                                                     v                                                  |
|  [Phổ kế HPGe Canberra (Genie 2000)] <=============+=============> [Phổ kế NaI(Tl) 4"x4" (Maestro-32)] |
|  - Độ phân giải FWHM 1.80 keV at 1332 keV                           - Buồng chì A340 nặng 1200 kg      |
|  - Chuẩn nội hiệu suất 214Bi                                        - Giải ma trận 3 cửa sổ năng lượng |
+--------------------------------------------------------------------------------------------------------+
  • Tỉnh Bolikhamxay ($6$ huyện): Gồm $1$ huyện đồng bằng (Paksan) và $5$ huyện miền núi (Khamkeut, Pakkading, Borikhane, Viengthong, Xaychamphone). Điểm lấy mẫu bao phủ cả khu vực mỏ vàng Khamkeut.
  • Tỉnh Khammouane ($9$ huyện): Gồm $6$ huyện đồng bằng (Thakhek, Nongbok, Hinboun, Mahaxay, Yommalath, Xebangfay) và $3$ huyện miền núi (Bualapha, Nakai, Gnommalath).
  • Tỉnh Savannakhet ($15$ huyện): Gồm $13$ huyện đồng bằng (Kaysone Phomvihane, Outhoumphone, Atsaphangthong, Phine, Songkhone, Champhone, Xonbuly, Xayboury, Virabury, Assaphone, Thapangthong, Phalanxay, Sayphouthong) và $2$ huyện miền núi (Sepone, Vilabuly). Điểm lấy mẫu bao trùm vùng khoáng hóa phức tạp xung quanh mỏ vàng Sepone và mỏ đồng-vàng Vilabuly.

Tổng số mẫu đất bề mặt nguyên dạng thu thập và xử lý phân tích đạt $202$ mẫu.

Quy trình nghiên cứu rigorous

Quy trình thu thập, tiền xử lý và đo đạc mẫu tuân thủ nghiêm ngặt hướng dẫn kỹ thuật của IAEA (1989):

  1. Giao thức lấy mẫu tại hiện trường: Tại mỗi điểm quan trắc đại diện, đất được lấy tại tầng bề mặt có độ sâu từ $0\text{ cm}$ đến $30\text{ cm}$ (hoặc từ $5\text{ cm}$ đến $25\text{ cm}$ ở các vị trí có cỏ rác phủ dày) theo phương pháp đường chéo $5$ điểm trên diện tích $100\text{ m}^2$, sau đó trộn đều theo kỹ thuật chia tư góc để lấy $1,5 - 2,0\text{ kg}$ mẫu đại diện.
  2. Tiền xử lý và tạo mẫu phân tích: Mẫu đất được sấy khô đồng nhất trong tủ sấy chuyên dụng ở nhiệt độ $105^\circ\text{C} - 110^\circ\text{C}$ trong $24$ giờ cho đến khi khối lượng không đổi để loại bỏ độ ẩm tự do. Sau đó mẫu được nghiền mịn bằng cối mã não và sàng qua rây tiêu chuẩn có kích thước mắt lưới từ $0,2\text{ mm}$ đến $1,0\text{ mm}$.
  3. Thiết lập cân bằng phóng xạ: Mẫu đất sau khi xử lý được cân chính xác bằng cân phân tích điện tử độ nhạy $0,1\text{ mg}$, nạp vào hộp hình trụ tiêu chuẩn (hoặc hộp Marinelli/hộp giếng thể tích $300 - 1000\text{ cm}^3$) và dán kín bằng keo chuyên dụng trong thời gian tối thiểu $30$ ngày nhằm ngăn chặn sự thất thoát của khí $^{222}\text{Rn}$ ($T_{1/2} = 3,82$ ngày) và $^{220}\text{Rn}$ ($T_{1/2} = 55,6$ giây), đảm bảo đạt trạng thái cân bằng phóng xạ thế kỷ tuyệt đối ($>99,5%$) giữa $^{226}\text{Ra}$ và các đồng vị con cháu $^{214}\text{Pb}$, $^{214}\text{Bi}$.
  4. Hệ thống thiết bị đo phổ gamma:
    • Hệ phổ kế gamma nhấp nháy $\text{NaI(Tl)}$: Đặt tại Bộ môn Vật lý hạt nhân, Trường ĐHKHTN – ĐHQGHN, sử dụng tinh thể $\text{NaI(Tl)}$ hình trụ kích thước $10,2 \times 10,2\text{ cm}$ ($4" \times 4"$), ống nhân quang điện ORTEC Model 276, điện áp cao thế $900\text{ V}$, hệ số khuếch đại $2,5$. Detector được bọc trong buồng chì che chắn Model A340 (North American Technology Services) có tổng khối lượng $1200\text{ kg}$, bề dày thân chì $10,2\text{ cm}$, bề dày nắp $13\text{ cm}$, đường kính trong $27,5\text{ cm}$, mặt trong lót lớp đồng nguyên chất dày $1\text{ mm}$ để triệt tiêu tia X huỳnh quang của chì. Hệ tích hợp phần mềm phân tích phổ MAESTRO-MCB25.
    • Hệ phổ kế gamma bán dẫn siêu tinh khiết HPGe: Sử dụng $4$ hệ phổ kế HPGe độ phân giải cao tại các viện nghiên cứu đầu ngành:
      • Hệ A (Viện Vật lý): Đầu dò ORTEC GEM20P4-70, hiệu suất tương đối $20%$, độ phân giải $FWHM = 1,80\text{ keV}$ tại đỉnh $1332\text{ keV}$ của $^{60}\text{Co}$, tỷ số Đỉnh/Compton $52:1$, cao thế $3800\text{ V}$, phần mềm Maestro-32.
      • Hệ B (Viện Khoa học và Kỹ thuật Hạt nhân): Đầu dò Canberra BE-6530, hiệu suất tương đối $60%$, độ phân giải $FWHM = 1,90\text{ keV}$, tỷ số Đỉnh/Compton $75:1$, cao thế $4500\text{ V}$, phần mềm Genie 2000.
      • Hệ C (Trung tâm Môi trường và Hóa học Quân sự): Đầu dò Canberra GC1518, hiệu suất $15%$, độ phân giải $FWHM = 1,80\text{ keV}$, tỷ số Đỉnh/Compton $44:1$, cao thế $3500\text{ V}$, phần mềm Genie 2000.
      • Hệ D (Viện Y học Phóng xạ và U bướu Quân đội): Đầu dò Canberra Model GC1518 (S/N 11037715), hiệu suất $15%$, $FWHM = 1,80\text{ keV}$, tỷ số Đỉnh/Compton $44:1$, cao thế $3500\text{ V}$, phần mềm Genie 2000.

Data và phân tích

Độ tin cậy của thuật toán phân tích được kiểm chứng thông qua mẫu chuẩn quốc tế IAEA-375, các mẫu chuẩn RGU-1, RGTh-1, RGK-1 và mẫu chuẩn thứ cấp TN1, TNK1. Thời gian đo phổ của mỗi mẫu kéo dài từ $24\text{ giờ}$ ($86.400\text{ s}$) đến $54.000\text{ s}$ nhằm thu thập phổ có diện tích đỉnh thuần (net peak area) đủ lớn, giảm độ bất định thống kê đếm xung xuống dưới $3% - 5%$.

+--------------------------------------------------------------------------------------------------------+
|                                    CHUỖI TRUYỀN SAI SỐ PHÉP ĐO                                         |
|                                                                                                        |
|  u(Net Area N) [Thống kê đếm xung Poisson] --------------------+                                      |
|  u(Hiệu suất ε) [Hiệu chuẩn đường cong & mật độ] -------------> [TỔNG HỢP SAI SỐ TOÀN PHẦN]           |
|  u(Tỷ số phân nhánh Iγ) [Dữ liệu hạt nhân chuẩn hóa] ---------> u_c(A) / A <= 5% - 8%                 |
|  u(Khối lượng m) [Cân phân tích điện tử 0.1 mg] --------------+                                      |
+--------------------------------------------------------------------------------------------------------+

Công thức truyền sai số tổng hợp $u_c(A)$ cho hoạt độ riêng $A$: $$\frac{u_c(A)}{A} = \sqrt{\left(\frac{u(N)}{N}\right)^2 + \left(\frac{u(\varepsilon)}{\varepsilon}\right)^2 + \left(\frac{u(I_\gamma)}{I_\gamma}\right)^2 + \left(\frac{u(m)}{m}\right)^2}$$

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả đo đạc thực nghiệm trên $202$ mẫu đất đá bề mặt tại $30$ huyện thuộc $3$ tỉnh miền Trung Lào mang lại các phát hiện khoa học quan trọng:

+--------------------------------------------------------------------------------------------------------+
|                             TƯƠNG QUAN HOẠT ĐỘ PHÓNG XẠ TỰ NHIÊN                                       |
|                                                                                                        |
|  Hoạt độ riêng (Bq/kg)                                                                                 |
|   400 |--------------------------------------------- [Trung bình Thế giới 40K = 400]                   |
|   300 |================== [40K Miền Trung Lào = 272.0]                                                 |
|   200 |                                                                                                |
|   100 |                                                                                                |
|    54 |================== [232Th Miền Trung Lào = 54.2] > [Thế giới 232Th = 30.0]                      |
|    38 |================== [226Ra Miền Trung Lào = 37.8] > [Thế giới 226Ra = 35.0]                      |
|     0 +------------------------------------------------------------------------>                       |
|                                    Đồng vị Phóng xạ Tự nhiên                                           |
+--------------------------------------------------------------------------------------------------------+
  1. Hoạt độ phóng xạ riêng của $^{40}\text{K}$, $^{226}\text{Ra}$ và $^{232}\text{Th}$:
    • Hoạt độ riêng của $^{40}\text{K}$ biến thiên trong dải từ $45,2\text{ Bq/kg}$ đến $680,5\text{ Bq/kg}$, với giá trị trung bình toàn vùng là $272,0 \pm 15,2\text{ Bq/kg}$, thấp hơn đáng kể so với trung bình thế giới ($400\text{ Bq/kg}$) và trung bình của Việt Nam ($292,6\text{ Bq/kg}$).
    • Hoạt độ riêng của $^{226}\text{Ra}$ dao động từ $12,1\text{ Bq/kg}$ đến $98,4\text{ Bq/kg}$, giá trị trung bình đạt $37,8 \pm 2,4\text{ Bq/kg}$, cao hơn mức trung bình toàn cầu ($35\text{ Bq/kg}$).
    • Hoạt độ riêng của $^{232}\text{Th}$ dao động trong khoảng từ $15,3\text{ Bq/kg}$ đến $124,6\text{ Bq/kg}$, với giá trị trung bình là $54,2 \pm 3,1\text{ Bq/kg}$, cao vượt trội ($1,81\text{ lần}$) so với giá trị trung bình toàn cầu ($30\text{ Bq/kg}$) và tương đương với mức nền của Việt Nam ($50,7\text{ Bq/kg}$) hay Thái Lan ($51,0\text{ Bq/kg}$).
  2. Dị thường phóng xạ tại các vùng khoáng hóa địa chất: Tại các huyện có mỏ vàng và thành tạo đá magma granit xâm nhập—cụ thể là huyện Khamkeut (tỉnh Bolikhamxay), huyện Sepone và huyện Vilabuly (tỉnh Savannakhet)—hoạt độ phóng xạ riêng của $^{226}\text{Ra}$ và $^{232}\text{Th}$ đạt cực đại, lần lượt vượt ngưỡng $95\text{ Bq/kg}$ và $120\text{ Bq/kg}$. Điều này xác nhận quy luật địa hóa phóng xạ: các nguyên tố nhóm Thori và Uran có xu hướng giàu hóa mạnh trong các pha magma axit muộn và quặng hóa đa kim cộng sinh.
  3. Tổng hợp hoạt độ Radi tương đương ($Ra_{eq}$): Giá trị $Ra_{eq}$ trên toàn miền Trung Lào dao động trong khoảng từ $48,6\text{ Bq/kg}$ đến $285,4\text{ Bq/kg}$, giá trị trung bình là $136,2 \pm 8,5\text{ Bq/kg}$. Tất cả các giá trị đo được đều nằm dưới giới hạn an toàn tối đa cho phép đối với vật liệu xây dựng và môi trường sống của UNSCEAR ($370\text{ Bq/kg}$).
  4. Suất liều hấp thụ và các chỉ số nguy hiểm bức xạ:
    • Suất liều hấp thụ trong không khí ở độ cao $1\text{ m}$ ($D_{air}$) dao động từ $22,4\text{ nGy/h}$ đến $128,7\text{ nGy/h}$, giá trị trung bình là $62,5 \pm 4,2\text{ nGy/h}$, cao hơn mức trung bình thế giới ($59\text{ nGy/h}$) và Việt Nam ($54,5\text{ nGy/h}$).
    • Suất liều hiệu dụng chiếu ngoài hàng năm ($AEDO$) biến thiên từ $0,028\text{ mSv/năm}$ đến $0,158\text{ mSv/năm}$, trung bình là $0,077 \pm 0,005\text{ mSv/năm}$, thấp hơn giới hạn khuyến cáo của ICRP cho công chúng ($1,0\text{ mSv/năm}$).
    • Chỉ số nguy hiểm chiếu ngoài ($H_{ex}$) dao động từ $0,13$ đến $0,77$ (trung bình $0,37 \pm 0,02$) và chỉ số nguy hiểm chiếu trong ($H_{in}$) dao động từ $0,16$ đến $0,98$ (trung bình $0,47 \pm 0,03$). Tất cả các chỉ số đều thỏa mãn điều kiện $H_{ex}, H_{in} < 1$.
+--------------------------------------------------------------------------------------------------------+
|                                  CÁC CHỈ SỐ AN TOÀN BỨC XẠ MIỀN TRUNG LÀO                              |
|                                                                                                        |
|  [Raeq: 136.2 Bq/kg]  =====================> [Giới hạn An toàn UNSCEAR: 370 Bq/kg]  (Đạt 36.8%)        |
|  [Hex: 0.37]          =====================> [Ngưỡng Nguy hiểm Tối đa: 1.00]        (Đạt 37.0%)        |
|  [Hin: 0.47]          =====================> [Ngưỡng Nguy hiểm Tối đa: 1.00]        (Đạt 47.0%)        |
|  [AEDO: 0.077 mSv/y]  =====================> [Giới hạn Công chúng ICRP: 1.00 mSv/y]  (Đạt  7.7%)        |
+--------------------------------------------------------------------------------------------------------+

Implications đa chiều

  • Đóng góp phương pháp luận: Khẳng định tính khả thi của việc kết hợp phương pháp ma trận trên hệ phổ kế $\text{NaI(Tl)}$ chi phí thấp với kỹ thuật chuẩn nội hiệu suất ghi trên hệ HPGe. Phương pháp này có thể chuyển giao cho các phòng thí nghiệm quốc gia đang phát triển để quan trắc phóng xạ môi trường diện rộng mà không đòi hỏi nguồn chuẩn trùng khớp tuyệt đối về mật độ chất nền.
  • Ứng dụng thực tiễn: Bộ số liệu của $30$ huyện cung cấp cơ sở định mức cho ngành công nghiệp khai khoáng và sản xuất vật liệu xây dựng (gạch, ngói, xi măng, đá ốp lát) tại Lào, giúp kiểm soát nguồn gốc đất sét và đá khai thác không vượt quá ngưỡng hoạt độ an toàn.
  • Khuyến nghị chính sách: Cung cấp cơ sở khoa học định lượng cho Bộ Khoa học và Công nghệ Lào, Bộ Tài nguyên và Môi trường Lào trong việc ban hành Quy chuẩn Kỹ thuật Quốc gia về an toàn bức xạ và xây dựng Bản đồ phông phóng xạ môi trường quốc gia.

Limitations và Future Research

  1. Giới hạn độ sâu lấy mẫu: Nghiên cứu chỉ tập trung vào lớp đất đá bề mặt tầng canh tác ($0 - 30\text{ cm}$), chưa tiến hành khoan thăm dò theo chiều sâu địa tầng để đánh giá gradient di chuyển của các hạt nhân phóng xạ theo nước ngầm.
  2. Giới hạn số lượng mẫu tại các điểm khoáng hóa: Mặc dù bao phủ $30$ huyện, mật độ lấy mẫu tại các vùng lõi mỏ vàng (Sepone, Khamkeut) mới dừng lại ở mức khảo sát ban đầu, chưa tạo thành lưới ô vuông đo đạc mật độ cao ($100 \times 100\text{ m}$).
  3. Đo trực tiếp Radon khí quyển: Chưa triển khai các thiết bị đo tích lũy Radon ($^{222}\text{Rn}$) trong nhà (indoor radon) bằng đầu dò vết hạt nhân CR-39 để đánh giá trực tiếp liều chiếu trong đối với các hộ dân sống gần vùng mỏ.
  4. Hướng nghiên cứu tương lai:
    • Mở rộng mạng lưới quan trắc phông phóng xạ cho các tỉnh miền Bắc và miền Nam Lào để hoàn thiện Bản đồ phóng xạ tự nhiên toàn quốc.
    • Ứng dụng mô hình tính toán Monte Carlo (phần mềm MCNP, Geant4) để mô phỏng chính xác hàm phản ứng của đầu dò và hiệu ứng tự hấp thụ cho các cấu trúc hình học mẫu phức tạp.
    • Nghiên cứu hệ số chuyển dịch sinh học (transfer factor) của $^{226}\text{Ra}$, $^{232}\text{Th}$, $^{40}\text{K}$ từ đất sang cây lương thực (lúa gạo, ngô) và chuỗi thức ăn tại các lưu vực sông miền Trung Lào.

Tác động và ảnh hưởng

Công trình của NCS. Somsavath Leuangtakoun đặt nền móng học thuật cho chuyên ngành Địa vật lý hạt nhân và Bức xạ môi trường tại CHDCND Lào:

  • Giá trị học thuật: Tạo ra bộ dữ liệu tham chiếu chuẩn hóa đầu tiên về hoạt độ phóng xạ tự nhiên của miền Trung Lào được cộng đồng khoa học quốc tế công nhận, cung cấp dữ liệu đầu vào cho các báo cáo định kỳ của UNSCEAR và IAEA.
  • Tác động kinh tế - xã hội: Giúp các cơ quan quản lý môi trường Lào giám sát tác động môi trường của các dự án thủy điện, dự án khai thác mỏ vàng, đồng, quặng bô xít và thạch cao dọc hành lang kinh tế Đông - Tây.
  • Hợp tác quốc tế: Thắt chặt liên kết nghiên cứu chuyên sâu giữa Trường Đại học Khoa học Tự nhiên – ĐHQGHN với các trường đại học, viện nghiên cứu tại Lào trong đào tạo nhân lực trình độ cao về kỹ thuật hạt nhân và an toàn bức xạ.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận quy trình giải bài toán ma trận phổ kế $\text{NaI(Tl)}$ và kỹ thuật chuẩn nội hiệu suất ghi không phụ thuộc chất nền.
  • Các nhà nghiên cứu địa chất và môi trường: Sử dụng số liệu nền để đối sánh sự tích tụ kim loại nặng và ô nhiễm phóng xạ do khai khoáng.
  • Doanh nghiệp khai khoáng và sản xuất vật liệu: Chủ động sàng lọc nguyên liệu đạt chuẩn $Ra_{eq} < 370\text{ Bq/kg}$.
  • Cơ quan hoạch định chính sách: Ban hành hạn ngạch liều bức xạ công cộng và xây dựng mạng lưới trạm quan trắc phóng xạ tự động quốc gia.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Đó là việc mở rộng lý thuyết hàm phản ứng đầu dò và chuẩn hóa phương pháp giải ma trận phân tách phổ chồng chập trên hệ phổ kế gamma nhấp nháy $\text{NaI(Tl)}$ kết hợp kỹ thuật chuẩn nội hiệu suất ghi trên hệ HPGe, khắc phục triệt để sai số do bất tương thích chất nền (matrix effect) giữa mẫu chuẩn và mẫu phân tích môi trường thực tế.
  2. Cải tiến phương pháp luận so với các nghiên cứu tiền nhiệm như thế nào? So với phương pháp chuẩn ngoại tuyệt đối truyền thống (vốn đòi hỏi mẫu chuẩn phải có cùng khối lượng riêng, hình học và thành phần hóa học chính xác với mẫu đất thực nghiệm—điều gần như bất khả thi đối với đất đá tự nhiên đa dạng), luận án sử dụng các đỉnh gamma tự thân của $^{214}\text{Bi}$ để xây dựng đường cong hiệu suất tương đối nội tại, sau đó chuẩn hóa tuyệt đối qua vạch $1460,82\text{ keV}$ của $^{40}\text{K}$, loại trừ hoàn toàn sai số tự hấp thụ chất nền.
  3. Phát hiện thực nghiệm gây bất ngờ nhất là gì? Hoạt độ riêng trung bình của $^{232}\text{Th}$ tại miền Trung Lào ($54,2\text{ Bq/kg}$) cao hơn $1,81\text{ lần}$ so với giá trị trung bình thế giới ($30\text{ Bq/kg}$), trong khi hoạt độ $^{40}\text{K}$ ($272,0\text{ Bq/kg}$) lại thấp hơn $32%$ so với trung bình thế giới ($400\text{ Bq/kg}$). Điều này phản ánh đặc tính vỏ phong hóa nhiệt đới sâu làm rửa trôi mạnh cation Kali dễ hòa tan, trong khi Thori có tính linh động địa hóa kém bị giữ lại và làm giàu trong các khoáng vật hạt mịn.
  4. Quy trình tái lập (replication protocol) có được cung cấp chi tiết không? Luận án cung cấp đầy đủ thông số hình học buồng đo, kích thước mẫu, quy trình tiền xử lý nhiệt $105^\circ\text{C}$, thời gian lưu mẫu $30$ ngày, điện áp làm việc, cấu hình phần mềm Genie 2000/Maestro, các ma trận hệ số $a_{ik}$ và phương pháp lan truyền sai số thống kê, cho phép mọi phòng thí nghiệm hạt nhân tái lập quy trình phân tích đạt độ chính xác tương đương.
  5. Chương trình nghị sự nghiên cứu 10 năm được định hình ra sao? Định hình lộ trình $3$ giai đoạn: (1) Số hóa bản đồ phông phóng xạ toàn bộ $18$ tỉnh/thành của Lào; (2) Thiết lập mạng lưới cảnh báo sớm bức xạ môi trường tại các vùng mỏ trọng điểm; (3) Tích hợp nghiên cứu mô hình hóa lan truyền phóng xạ sinh thái trong chuỗi nông nghiệp và nguồn nước sông Mê Kông.

Kết luận

  1. Luận án đã hoàn thiện cơ sở lý thuyết và phương pháp thực nghiệm nâng cao độ chính xác phép đo hoạt độ phóng xạ tự nhiên của $^{40}\text{K}$, $^{226}\text{Ra}$ và $^{232}\text{Th}$ trên cả hai hệ phổ kế gamma nhấp nháy $\text{NaI(Tl)}$ và bán dẫn HPGe.
  2. Thiết lập thành công phương pháp ma trận phân tích phổ đa cửa sổ trên hệ $\text{NaI(Tl)}$ và kỹ thuật đường cong chuẩn nội hiệu suất ghi kết hợp hiệu suất tuyệt đối tại đỉnh $1460,82\text{ keV}$ trên hệ HPGe, giải quyết dứt điểm sai số do hiệu ứng tự hấp thụ chất nền.
  3. Cung cấp bộ dữ liệu thực nghiệm chuẩn hóa đầu tiên về hoạt độ phóng xạ riêng của $202$ mẫu đất đá bề mặt tại $30$ huyện thuộc $3$ tỉnh miền Trung Lào: $^{40}\text{K} = 272,0 \pm 15,2\text{ Bq/kg}$, $^{226}\text{Ra} = 37,8 \pm 2,4\text{ Bq/kg}$, $^{232}\text{Th} = 54,2 \pm 3,1\text{ Bq/kg}$, và $Ra_{eq} = 136,2 \pm 8,5\text{ Bq/kg}$.
  4. Xác lập bộ bản đồ phóng xạ tự nhiên chi tiết cho $3$ tỉnh Bolikhamxay, Khammouane và Savannakhet, nhận diện các dị thường phóng xạ liên quan đến các thành tạo magma và mỏ vàng tại Khamkeut, Sepone, Vilabuly.
  5. Đánh giá toàn diện các chỉ số an toàn bức xạ: Suất liều hấp thụ $D_{air} = 62,5\text{ nGy/h}$, Liều hiệu dụng $AEDO = 0,077\text{ mSv/năm}$, Chỉ số nguy hiểm $H_{ex} = 0,37 < 1$ và $H_{in} = 0,47 < 1$, chứng minh môi trường đất đá bề mặt miền Trung Lào nằm trong ngưỡng an toàn bức xạ theo tiêu chuẩn quốc tế của UNSCEAR và IAEA.
  6. Mở ra các hướng nghiên cứu mới về trắc địa hạt nhân, quan trắc sinh thái phóng xạ và an toàn bức xạ công nghiệp, đóng góp trực tiếp vào sự nghiệp phát triển khoa học - công nghệ bền vững của CHDCND Lào.