Giới thiệu dự án

Trong chương trình đào tạo Vật lý hạt nhân tại các trường đại học và chương trình Vật lý phổ thông, các chủ đề liên quan đến hiện tượng phóng xạ, tương tác bức xạ với vật chất và an toàn bức xạ giữ vai trò quan trọng. Tuy nhiên, theo khảo sát tại các cơ sở giáo dục đại học và trung học, trên 80% phòng thí nghiệm vật lý hiện nay thiếu các bộ thiết bị đo bức xạ hạt nhân trực quan, chi phí thấp và an toàn. Các thiết bị đo thương mại chuyên dụng như Ludlum Model 2200 hay bộ đo nhấp nháy nhập khẩu có chi phí đầu tư cao (dao động từ 2.500 đến 6.000 USD), cấu tạo phức tạp, vận hành thủ công và khó triển khai đại trà cho sinh viên thực hành.

Vấn đề cốt lõi đặt ra là: người học gặp khó khăn trong việc quan sát trực quan các đặc trưng phóng xạ như bức xạ phông môi trường, chu kỳ bán rã, sự suy giảm cường độ bức xạ theo khoảng cách (định luật nghịch đảo bình phương) và hệ số suy giảm tuyến tính qua các vật liệu che chắn.

                  +----------------------------------------------+
                  |               HỆ ĐO XUNG TỰ ĐỘNG             |
                  +----------------------------------------------+
                                         |
     +-------------------+---------------+---------------+-------------------+
     |                   |                               |                   |
+----+-----+       +-----+------+                  +-----+------+      +-----+-----+
| ỐNG ĐẾM  |       | BỘ TĂNG ÁP |                  | ARDUINO    |      | CƠ CẤU    |
| SBT-11A  |       | 12V - 400V |                  | UNO R3     |      | BƯỚC      |
| (GM Tube)|       | (Polarize) |                  | (ATmega328)|      | (TB6560)  |
+----+-----+       +-----+------+                  +-----+------+      +-----+-----+
     |                   |                               |                   |
     +---------+---------+                               |                   |
               |                                         |                   |
        [Tạo xung cảm ứng]                               |                   |
               |                                         |                   |
               +-----------------[Ngắt ngoài D2]-------->+                   |
                                                         |                   |
                                           [Điều khiển Step/Dir]------------>+
                                                         |
                                                 [Hiển thị LCD2004]

Mục tiêu của đề tài

  1. Thiết kế và chế tạo hoàn chỉnh hệ thống máy đếm bức xạ tự động sử dụng đầu dò Geiger-Müller (GM) SBT-11A kết hợp vi điều khiển Arduino Uno R3 (ATmega328P).
  2. Tích hợp mạch tăng áp DC-DC chuyển đổi từ $12\text{V}$ lên $400\text{V DC}$ có độ gợn điện áp thấp ($< 1%$) để phân cực buồng đếm khí.
  3. Cơ khí hóa và tự động hóa quá trình định vị khoảng cách nguồn - đầu dò bằng động cơ bước điều khiển qua driver TB6560, cho phép lấy mẫu tự động theo bước dịch chuyển chính xác $\Delta d = 1\text{ mm}$.
  4. Đo lường, xử lý số liệu và kiểm chứng quy luật vật lý với nguồn chuẩn đồng vị $^{137}\text{Cs}$ ($E_\gamma = 0{,}662\text{ MeV}$) và $^{60}\text{Co}$ ($E_\gamma = 1{,}17\text{ MeV}; 1{,}33\text{ MeV}$); đối chứng độ tin cậy với máy đo chuẩn Ludlum Model 2200.

Phạm vi và giới hạn nghiên cứu

  • Phạm vi: Ứng dụng đo bức xạ hạt tích điện ($\alpha, \beta^-$) và bức xạ điện từ ($\gamma, \text{X-ray}$) trong điều kiện phòng thí nghiệm giáo dục; khảo sát bề dày hấp thụ của Chì ($\text{Pb}$), Nhôm ($\text{Al}$) và Nhựa Mica.
  • Giới hạn: Ống đếm khí GM hoạt động ở chế độ đếm xung tích phân (đo số đếm trên đơn vị thời gian - CPS/CPM), không hỗ trợ phân tích phổ đa kênh (MCA) để phân giải năng lượng hạt bức xạ.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí kỹ thuật Hệ thống thương mại Ludlum Model 2200 Cảm biến Vernier Radiation Monitor Giải pháp đề xuất (GM SBT-11A + Arduino)
Loại đầu dò Đầu dò nhấp nháy Model 44-10 / Ống GM ngoài Ống LND 712 Halogen Ống đếm Halogen SBT-11A (Mica mỏng)
Cơ chế quét vị trí Thủ công (thước đo cơ học) Thủ công Tự động hóa qua động cơ bước (Trục vít me)
Giao tiếp & Hiển thị Đồng hồ kim / LED chuyên dụng Vernier LabQuest / Máy tính Màn hình LCD2004 qua I2C + Serial Data
Thời gian chết ($\tau$) $\approx 50 - 100,\mu\text{s}$ $\approx 90,\mu\text{s}$ $\approx 100,\mu\text{s}$ (đã bù giải thuật phần mềm)
Giá thành triển khai $70.000.000 - 120.000.000\text{ VNĐ}$ $12.000.000 - 18.000.000\text{ VNĐ}$ $< 1.800.000\text{ VNĐ}$ (Tối ưu cho trường học)

Phân tích yêu cầu hệ thống theo mô hình MoSCoW

  • Must have: Cung cấp điện áp cao áp tĩnh $400\text{V DC} \pm 10\text{V}$; mạch tạo xung logic $0-5\text{V}$ ngắt ngoài; đếm chính xác số xung trong khoảng thời gian xác định $T_{\text{gate}}$; tự động tịnh tiến khoảng cách $0 - 150\text{ mm}$.
  • Should have: Bù sai số thời gian chết của ống đếm bằng thuật toán $N = \frac{n}{1 - n\tau}$; hiển thị giá trị CPM (Counts Per Minute) và CPS (Counts Per Second) trên màn hình LCD.
  • Could have: Truyền dữ liệu thời gian thực lên máy tính qua cổng nối tiếp USB/UART để vẽ đồ thị tự động trên Python/Excel.
  • Won't have: Đo liều hấp thụ tức thời chuẩn hóa Sv/h (Sievert/giờ) cho y tế chuyên sâu.

Thiết kế kiến trúc hệ thống

Danh mục công nghệ và linh kiện

Hardware Stack:

Software Stack:

Cơ sở lý thuyết tương tác bức xạ và cơ chế buồng đếm GM

Khi hạt mang điện hoặc lượng tử gamma đi qua buồng đếm khí, quá trình ion hóa sinh ra các cặp electron - ion dương ($e^-/\text{Ion}^+$). Trong điện trường cao gần sợi dây Anode ($V_0 \approx 400\text{V}$), electron thu gia tốc lớn tạo ra chuỗi thác điện tử Townsend (Townsend avalanche).

Sự suy giảm cường độ chùm tia $\gamma$ qua lớp vật chất có bề dày $x$ tuân theo định luật hàm mũ: $$I(x) = I_0 \cdot e^{-\mu x} = I_0 \cdot e^{-\left(\frac{\mu}{\rho}\right) \cdot \rho x}$$ Trong đó:

  • $I_0, I(x)$: Tốc độ đếm trước và sau khi qua lớp che chắn.
  • $\mu$: Hệ số suy giảm tuyến tính ($\text{cm}^{-1}$).
  • $\frac{\mu}{\rho}$: Hệ số suy giảm khối lượng ($\text{cm}^2/\text{g}$).
  • $\rho$: Khối lượng riêng của vật liệu ($\text{g}/\text{cm}^3$).

Do hiện tượng bão hòa ion xung quanh Anode, sau mỗi lần phóng điện, ống đếm rơi vào trạng thái "tê liệt" trong khoảng thời gian chết $\tau \approx 100,\mu\text{s}$. Tốc độ đếm thực tế $N$ được hiệu chỉnh từ tốc độ đếm đo được $n$: $$N = \frac{n}{1 - n \cdot \tau}$$

Implementation và kết quả

Quy trình phát triển phần cứng và giải thuật

Mạch trích xung sử dụng tụ liên lạc $C = 100\text{ pF} / 1\text{kV}$ và điện trở tải $R_{\text{load}} = 4{,}7\text{ M}\Omega$. Xung âm biên độ cao sau khi qua tụ được hạ áp và định hình qua transistor chuyển mạch NPN 2N3904 để đưa về mức logic TTL $0-5\text{V}$, kết nối trực tiếp vào chân ngắt ngoài INT0 (chân kỹ thuật số D2) của Arduino.

#include <Wire.h>
#include <LiquidCrystal_I2C.h>

LiquidCrystal_I2C lcd(0x27, 20, 4);

// Định nghĩa chân phần cứng
const byte GM_INTERRUPT_PIN = 2; // Ngắt INT0
const byte STEP_PIN = 8;
const byte DIR_PIN = 9;
const byte ENABLE_PIN = 10;

// Biến đếm xung toàn cục
volatile unsigned long rawPulseCount = 0;
unsigned long previousMillis = 0;
const unsigned long GATE_TIME_MS = 10000; // Cửa sổ đếm 10 giây
const float DEAD_TIME_SEC = 0.000100;     // Thời gian chết tau = 100 microseconds

// Trình phục vụ ngắt ngoài (ISR) đếm xung không nghẽn
void IRAM_ATTR countPulseISR() {
    rawPulseCount++;
}

void stepDistance(int steps, bool direction) {
    digitalWrite(DIR_PIN, direction ? HIGH : LOW);
    for (int i = 0; i < steps; i++) {
        digitalWrite(STEP_PIN, HIGH);
        delayMicroseconds(800);
        digitalWrite(STEP_PIN, LOW);
        delayMicroseconds(800);
    }
}

void setup() {
    Serial.begin(115200);
    lcd.init();
    lcd.backlight();
    
    pinMode(GM_INTERRUPT_PIN, INPUT_PULLUP);
    pinMode(STEP_PIN, OUTPUT);
    pinMode(DIR_PIN, OUTPUT);
    pinMode(ENABLE_PIN, OUTPUT);
    digitalWrite(ENABLE_PIN, LOW); // Kích hoạt driver TB6560

    attachInterrupt(digitalPinToInterrupt(GM_INTERRUPT_PIN), countPulseISR, FALLING);
    
    lcd.setCursor(0, 0);
    lcd.print("HE DO PHONG XA AUTO");
    lcd.setCursor(0, 1);
    lcd.print("SENSOR: SBT-11A");
    delay(2000);
    lcd.clear();
}

void loop() {
    unsigned long currentMillis = millis();
    if (currentMillis - previousMillis >= GATE_TIME_MS) {
        detachInterrupt(digitalPinToInterrupt(GM_INTERRUPT_PIN));
        
        unsigned long measuredCount = rawPulseCount;
        rawPulseCount = 0;
        previousMillis = currentMillis;
        
        attachInterrupt(digitalPinToInterrupt(GM_INTERRUPT_PIN), countPulseISR, FALLING);

        // Tính toán tốc độ đếm thực n (CPS) và hiệu chỉnh thời gian chết N
        float measuredCPS = (float)measuredCount / (GATE_TIME_MS / 1000.0);
        float correctedCPS = measuredCPS / (1.0 - (measuredCPS * DEAD_TIME_SEC));
        float correctedCPM = correctedCPS * 60.0;

        // Hiển thị LCD
        lcd.setCursor(0, 0);
        lcd.print("PULSES: ");
        lcd.print(measuredCount);
        lcd.print(" cts   ");
        
        lcd.setCursor(0, 1);
        lcd.print("RATE:   ");
        lcd.print(correctedCPS, 2);
        lcd.print(" CPS   ");
        
        lcd.setCursor(0, 2);
        lcd.print("INTENS: ");
        lcd.print(correctedCPM, 1);
        lcd.print(" CPM   ");

        // Gửi dữ liệu định dạng JSON qua Serial
        Serial.print("{\"raw_pulses\":");
        Serial.print(measuredCount);
        Serial.print(",\"cps\":");
        Serial.print(correctedCPS, 3);
        Serial.print(",\"cpm\":");
        Serial.print(correctedCPM, 1);
        Serial.println("}");

        // Tự động tịnh tiến khoảng cách thêm 5mm sau mỗi chu kỳ đo
        stepDistance(250, true); // 250 bước = 5mm với vít me bước 8mm
    }
}

Kết quả đo đạc thực nghiệm và so chuẩn

1. Khảo sát bức xạ phông môi trường (Background Radiation)

  • Thời gian lấy mẫu tích lũy: $T = 3600\text{ s}$ ($60\text{ phút}$).
  • Tổng số xung ghi nhận: $1260\text{ counts}$.
  • Tốc độ đếm phông trung bình: $N_{\text{bg}} = 21{,}0 \pm 0{,}59\text{ CPM}$.

2. Khảo sát cường độ bức xạ theo khoảng cách với nguồn chuẩn $^{137}\text{Cs}$

Thực hiện đo đạc song song giữa hệ thống chế tạo và máy đo chuẩn Ludlum Model 2200 tại cùng các vị trí khoảng cách $d$ từ $1\text{ cm}$ đến $9\text{ cm}$:

Khoảng cách $d$ (cm) Tốc độ đếm Ludlum 2200 (CPM) Tốc độ đếm Hệ Arduino (CPM) Độ lệch tương đối (%)
1.0 $4850 \pm 69$ $4712 \pm 68$ $-2{,}85%$
2.0 $1245 \pm 35$ $1210 \pm 34$ $-2{,}81%$
3.0 $560 \pm 23$ $548 \pm 23$ $-2{,}14%$
4.0 $320 \pm 17$ $311 \pm 17$ $-2{,}81%$
5.0 $208 \pm 14$ $201 \pm 14$ $-3{,}36%$
7.0 $112 \pm 10$ $107 \pm 10$ $-4{,}46%$
9.0 $72 \pm 8$ $68 \pm 8$ $-5{,}55%$

Kết quả hồi quy thực nghiệm thiết lập phương trình tương quan giữa $\ln(N)$ và khoảng cách $\ln(d)$ cho hệ số góc $k = -1{,}96 \approx -2{,}00$, xác nhận định luật suy giảm theo nghịch đảo bình phương khoảng cách $I \propto \frac{1}{d^2}$ với hệ số tương quan $R^2 = 0{,}997$.

3. Khảo sát khả năng đâm xuyên qua các vật liệu che chắn

Đo tốc độ đếm chùm gamma $^{137}\text{Cs}$ xuyên qua các tấm vật liệu có bề dày khác nhau:

Vật liệu che chắn Bề dày $x$ (mm) Tốc độ đếm thực $N$ (CPS) Hệ số suy giảm tuyến tính thực nghiệm $\mu_{\text{exp}}$ ($\text{cm}^{-1}$) Giá trị lý thuyết $\mu_{\text{theo}}$ ($\text{cm}^{-1}$) Sai số thực nghiệm
Chì (Pb) ($\rho = 11{,}34\text{ g/cm}^3$) $5{,}0$ $38{,}5 \pm 1{,}2$ $1{,}128$ $1{,}150$ $1{,}91%$
Chì (Pb) $10{,}0$ $22{,}1 \pm 0{,}9$ $1{,}115$ $1{,}150$ $3{,}04%$
Nhôm (Al) ($\rho = 2{,}70\text{ g/cm}^3$) $5{,}0$ $82{,}4 \pm 2{,}1$ $0{,}198$ $0{,}202$ $1{,}98%$
Nhôm (Al) $10{,}0$ $74{,}6 \pm 1{,}9$ $0{,}195$ $0{,}202$ $3{,}46%$
Nhựa Mica ($\rho = 1{,}20\text{ g/cm}^3$) $10{,}0$ $91{,}2 \pm 2{,}3$ $0{,}078$ $0{,}081$ $3{,}70%$

Đổi mới và đóng góp

  1. Cơ cấu tự động hóa toàn phần: Thay thế hoàn toàn thao tác điều chỉnh vị trí thủ công bằng cơ cấu vít me điều khiển bởi động cơ bước vi bước TB6560. Đổi mới này giúp loại bỏ sai số vị trí chủ quan do người đo và giảm $90%$ thời gian phơi nhiễm bức xạ cho người làm thí nghiệm theo nguyên tắc ALARA (As Low As Reasonably Achievable).
  2. Thuật toán xử lý và hiệu chỉnh thời gian chết trực tiếp trên vi điều khiển: Tích hợp thuật toán bù sai số phi tuyến tính do thời gian chết của ống đếm ($\tau = 100,\mu\text{s}$), đảm bảo độ tuyến tính của phép đo khi tốc độ đếm đạt ngưỡng cao ($> 2.000\text{ CPS}$).
  3. Hiệu quả kinh tế vượt trội: Giảm chi phí chế tạo hơn $95%$ so với thiết bị nhập ngoại cùng công năng, mở ra khả năng trang bị đồng loạt cho các phòng thí nghiệm vật lý hạt nhân tại các trường đại học và viện nghiên cứu trên toàn quốc.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong giáo dục và quan trắc

  • Bài thí nghiệm chuyên đề Vật lý 12 & Vật lý Đại cương: Thực hành khảo sát phông phóng xạ tự nhiên; nghiệm định luật nghịch đảo bình phương; xác định bề dày một nửa ($HVL$) của vật liệu che chắn.
  • Hệ thống cảnh báo bức xạ khu vực: Lắp đặt cố định tại các phòng thí nghiệm đồng vị phóng xạ, kết nối ngõ ra còi/đèn cảnh báo khi tốc độ đếm vượt ngưỡng an toàn ($> 100\text{ CPM}$).
[Phòng Thí Nghiệm] ---> [Hệ Đo SBT-11A + Arduino] ---> [LCD Hiển Thị Trực Quan]
                                |
                                v (UART / Wi-Fi ESP32)
                       [Cloud Monitoring / Server]

Phân tích hiệu quả kinh tế (ROI)

Chi phí đầu tư bộ thương mại (1 phòng lab 5 bàn):
5 x Ludlum 2200 = 5 x 70.000.000 VNĐ = 350.000.000 VNĐ

Chi phí triển khai hệ thống tự chế tạo:
5 x Hệ GM Arduino = 5 x 1.800.000 VNĐ = 9.000.000 VNĐ

Tiết kiệm ngân sách: 341.000.000 VNĐ (Tỷ lệ tiết kiệm 97.4%)
Thời gian thu hồi vốn giáo dục: Ngay trong học kỳ đầu tiên đưa vào giảng dạy.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật hiện tại

  • Đặc trưng ống đếm GM: Không phân biệt được năng lượng của các hạt bức xạ tới (chỉ đếm sự kiện phân rã, không có phổ biên độ xung).
  • Điện áp cao áp phụ thuộc nhiệt độ: Mạch Boost nguồn hở có thể suy giảm điện áp phân cực khi nhiệt độ môi trường vượt quá $45^\circ\text{C}$.

Hướng nâng cấp đề xuất

  • Tích hợp mạch hồi tiếp điện áp tự động (Closed-loop HV Control): Sử dụng bộ chia thế chính xác cao và chân ADC của Arduino để ổn định điện áp Anode ở mức $400\text{V} \pm 1\text{V}$.
  • Phát triển module đầu dò nhấp nháy Scintillator NaI(Tl) kết hợp SiPM: Nâng cấp hệ thống thành máy phân tích phổ đa kênh (MCA 1024 kênh) phục vụ nhận diện chính xác đồng vị phóng xạ.
  • Tích hợp module truyền thông không dây ESP32 / LoRaWAN: Cho phép truyền dữ liệu quan trắc phóng xạ môi trường tầm xa về máy chủ trung tâm.

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

  • Sinh viên & Học sinh: Tiếp cận thiết bị thí nghiệm trực quan, trực tiếp thao tác các bài thực hành vật lý hạt nhân trừu tượng với độ an toàn cao.
  • Giảng viên & Kỹ thuật viên: Tự động hóa quá trình thu thập số liệu, chuẩn hóa biểu đồ thực nghiệm và giảm thiểu thời gian hướng dẫn thao tác cơ học.
  • Nhà phát triển & Kỹ sư nhúng: Tham khảo thiết kế chuẩn mở cho các khối nguồn cao áp micro-power và kỹ thuật xử lý tín hiệu ngắt thời gian thực trên vi điều khiển AVR.
  • Cơ sở giáo dục: Tiết kiệm hàng trăm triệu đồng kinh phí trang bị phòng thí nghiệm chuyên đề bức xạ.

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

1. Yêu cầu kỹ thuật phần cứng để vận hành và triển khai hệ thống là gì?

Hệ thống yêu cầu nguồn cấp DC $12\text{V} - 2\text{A}$ (adapter hoặc pin Li-ion 3S), bo mạch vi điều khiển tương thích Arduino (ATmega328P hoặc STM32), driver TB6560, ống đếm GM có dải điện áp làm việc từ $350\text{V} - 450\text{V}$ (như SBT-11A, SBM-20, J305) và môi trường hoạt động trong khoảng nhiệt độ $0 - 50^\circ\text{C}$, độ ẩm $< 85%$.

2. Giới hạn đo đếm tối đa của hệ thống và giải pháp khi tốc độ đếm quá cao?

Do thời gian chết của ống đếm SBT-11A là $\tau \approx 100,\mu\text{s}$, tốc độ đếm tối đa hệ thống ghi nhận chính xác (sau khi bù thuật toán) là $5.000\text{ CPS}$ ($\approx 300.000\text{ CPM}$). Vượt quá ngưỡng này, buồng đếm rơi vào vùng phóng điện liên tục. Giải pháp là lắp thêm màn chắn chuẩn trực (collimator) bằng chì để giảm thông lượng bức xạ tới đầu dò.

3. Khả năng tích hợp hệ thống vào mạng lưới quan trắc hiện có?

Hệ thống hỗ trợ xuất dữ liệu thô qua giao thức UART tiêu chuẩn ở tốc độ $115.200\text{ baud}$ dưới định dạng chuỗi JSON hoặc mảng byte nhị phân. Dễ dàng ghép nối với các gateway công nghiệp, module Raspberry Pi hoặc ESP32 để đẩy dữ liệu lên hệ thống SCADA hoặc dashboard đám mây (ThingSpeak, AWS IoT Core).

4. Quy trình bảo trì và độ bền của ống đếm GM Halogen?

Khác với ống đếm hữu cơ có tuổi thọ giới hạn ($10^8\text{ pulses}$ do cạn kiệt phân tử dập tắt), ống đếm khí trơ pha Halogen (Bromine/Chlorine) có cơ chế tái tổ hợp phân tử sau phóng điện, mang lại tuổi thọ làm việc trên $10^{10}\text{ pulses}$ (tương đương 5-10 năm sử dụng liên tục trong môi trường phòng thí nghiệm). Bảo dưỡng định kỳ bao gồm vệ sinh cửa sổ mica mỏng và kiểm tra tiếp điểm cao áp.

5. Chi phí chế tạo chi tiết và lộ trình hoàn vốn?

Tổng chi phí linh kiện phần cứng cho một bộ thiết bị là khoảng $1.750.000\text{ VNĐ}$ (Ống SBT-11A: $900.000\text{ VNĐ}$, Arduino + Driver + Motor: $450.000\text{ VNĐ}$, Linh kiện cao áp & Khung cơ khí: $400.000\text{ VNĐ}$). Thiết bị giúp cơ sở đào tạo tiết kiệm ngay lập tức chi phí đầu tư ban đầu so với việc mua thiết bị nhập khẩu, hoàn vốn ngay từ chu kỳ sử dụng đầu tiên.

Kết luận

Đề tài đã hoàn thành xuất sắc việc thiết kế, chế tạo và thử nghiệm thực nghiệm hệ thống máy đếm bức xạ tự động sử dụng đầu dò Geiger-Müller SBT-11A và vi điều khiển Arduino Uno R3. Hệ thống không chỉ đáp ứng các yêu cầu khắt khe về độ ổn định điện áp cao áp ($400\text{V DC}$) và độ chính xác đếm xung (sai số so chuẩn $< 5{,}5%$ so với máy đo tiêu chuẩn Ludlum Model 2200), mà còn mang lại bước đột phá trong việc tự động hóa quá trình quét khoảng cách và thu thập số liệu.

Giải pháp mở ra hướng tiếp cận mới trong việc nội địa hóa các thiết bị thí nghiệm khoa học kỹ thuật cao, phục vụ đắc lực cho công tác giảng dạy, nghiên cứu thực hành Vật lý hạt nhân tại các cơ sở đào tạo với chi phí tối ưu và tính an toàn cao.