Tổng quan nghiên cứu

Theo thống kê từ Tổ chức Y tế Thế giới, ung thư hiện là một trong những nguyên nhân gây tử vong hàng đầu toàn cầu với khoảng gần 10 triệu ca mắc mới và hơn 8 triệu ca tử vong mỗi năm. Tại Việt Nam, mỗi năm ghi nhận khoảng 150.000 ca mắc mới cùng hơn 50.000 trường hợp tử vong do các bệnh lý ác tính. Trong số các kỹ thuật can thiệp hiện đại, xạ phẫu định vị lập thể bằng dao gamma là giải pháp tối ưu cho các tổn thương nội sọ phức tạp, thu hút hơn 500.000 bệnh nhân trên toàn thế giới lựa chọn điều trị kể từ năm 1968.

Vấn đề nghiên cứu cốt lõi xuất phát từ thực tiễn vận hành hệ thống xạ phẫu dao gamma quay GammaART-6000TM tại Trung tâm Y học hạt nhân và Ung bướu, Bệnh viện Bạch Mai. Khác với hệ thống dao gamma cổ điển Leksell cố định 201 nguồn, hệ thống xoay RGS sử dụng 30 nguồn Cobalt-60 với tổng hoạt độ ban đầu xấp xỉ 6.000 Ci. Việc tích hợp cơ chế quay đồng bộ giữa khối nguồn và ống chuẩn trực thứ cấp đòi hỏi một quy trình đảm bảo chất lượng (QA) nghiêm ngặt nhằm kiểm soát độ chính xác cơ học và liều lượng bức xạ.

Mục tiêu cụ thể của luận văn là thiết lập quy trình kiểm chuẩn QA toàn diện theo tiêu chuẩn quốc tế và ứng dụng phương pháp mô phỏng Monte Carlo qua chương trình MCNP5 để tính toán phân bố liều 3D. Nghiên cứu được thực hiện tại Bệnh viện Bạch Mai và Đại học Khoa học Tự nhiên Hà Nội trong giai đoạn 2015–2016. Kết quả nghiên cứu mang ý nghĩa then chốt khi kiểm soát sai số đồng tâm dưới 0,3 mm, sai lệch liều thực nghiệm dưới 5%, đồng thời cung cấp công cụ tính toán độc lập giúp tối ưu hóa kế hoạch xạ phẫu và nâng cao an toàn bức xạ cho bệnh nhân.

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 dựa trên lý thuyết tương tác bức xạ gamma với vật chất sống và lý thuyết vận chuyển hạt ngẫu nhiên Monte Carlo. Năng lượng photon phát ra từ nguồn Cobalt-60 ở hai mức đỉnh 1,17 MeV và 1,33 MeV tạo ra chuỗi tương tác ion hóa trực tiếp và gián tiếp. Trong cơ thể người với hơn 85% là nước, khoảng 80% tác động sinh học phá hủy tế bào u bắt nguồn từ cơ chế gián tiếp thông qua các gốc tự do peroxy và ion hoạt hóa.

Mô hình nghiên cứu kết hợp giữa đo đạc thực nghiệm theo nghị định thư quốc tế IAEA TRS-398, IAEA TRS-227 và báo cáo AAPM TG-42, TG-51 của Hiệp hội Vật lý Y học Hoa Kỳ. Ba khái niệm nền tảng trong hệ thống RGS bao gồm:

  • Suất liều hấp thụ tuyệt đối: Đại lượng đo lường năng lượng bức xạ hấp thụ trên một đơn vị khối lượng mô mô phỏng tại tâm hình học, đạt giá trị lớn hơn hoặc bằng 3 Gy/phút khi nguồn mới nạp.
  • Hệ số trường chiếu (Output Factor - OF): Tỷ số giữa suất liều tại kích thước trường chuẩn trực cụ thể so với trường chuẩn 18 mm.
  • Độ rộng nửa cực đại (FWHM): Chỉ số đặc trưng cho độ phân bố không gian của chùm tia bức xạ và vùng chuyển tiếp liều (penumbra).

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

Nguồn dữ liệu thực nghiệm được thu thập trực tiếp trên hệ thống GammaART-6000TM tại Bệnh viện Bạch Mai. Nghiên cứu sử dụng buồng ion hóa thể tích siêu nhỏ A1SL (thể tích 0,053 cm3, hệ số chuyển đổi 55,8 R/nC) kết hợp liều kế Electrometer MAX-4000, đặt trong phantom cầu Polystyrene đường kính 16 cm (mật độ khối lượng 1,16 g/cm3) tương đương mô não người trưởng thành. Cỡ mẫu thực nghiệm bao gồm 5 khoảng thời gian phát xạ (60s, 120s, 180s, 240s, 300s) được lặp lại 3 lần cho mỗi kích thước collimator (4 mm, 8 mm, 14 mm, 18 mm). Phương pháp chọn mẫu đo đạc được thiết lập theo quy trình chuẩn hóa của IAEA nhằm loại bỏ sai số ngẫu nhiên do dòng rò và biến thiên nhiệt độ, áp suất môi trường.

Về phương pháp phân tích mô phỏng, tác giả lựa chọn phần mềm MCNP5 (Monte Carlo N-Particle) dựa trên cơ sở dữ liệu hạt nhân ENDF/B và ACTL. Lý do lựa chọn MCNP5 là khả năng giải quyết các bài toán phân bố bức xạ không gian 3D phức tạp mà các phương pháp giải tích thông thường không thể đáp ứng chính xác. Mô hình tính toán được thiết lập qua hai giai đoạn: mô phỏng nguồn đơn kênh và mô phỏng tổ hợp 30 chùm tia xoay đồng quy tại góc quét 360 độ. Quá trình thu thập và xử lý số liệu kéo dài liên tục trong 12 tháng từ tháng 8/2015 đến tháng 8/2016.

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

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

Quá trình kiểm chuẩn thực nghiệm và mô phỏng tính toán đã mang lại 4 phát hiện quan trọng về đặc tính vật lý của máy dao gamma quay:

  • Độ chính xác hình học và hệ thống điều trị tự động (ART): Hệ thống định vị lập thể đạt sai số cơ học chỉ từ 0,0 mm đến 0,1 mm trên các trục tọa độ X, Y, Z. Độ sai lệch vị trí chuyển động của hệ thống ART luôn duy trì dưới 0,2 mm, thấp hơn nhiều so với ngưỡng cảnh báo 0,5 mm của nhà sản xuất.
  • Hệ số trường chiếu (Output Factor): Kết quả đo đạc thực nghiệm cho thấy hệ số OF của các ống chuẩn trực đạt độ tương thích cao với thông số nhà sản xuất ban hành: collimator 18 mm đạt 1,000; collimator 14 mm đạt 0,982 (sai lệch 0,2% so với giá trị chuẩn 0,980); collimator 8 mm đạt 0,910 (trùng khớp tuyệt đối 0,910); và collimator 4 mm đạt 0,569 (sai lệch 0,17% so với giá trị chuẩn 0,570).
  • Độ chính xác phân liều tại tâm điều trị: Tại tâm phantom cầu với liều chỉ định 10 Gy (tương ứng đường đồng liều 50% bao quanh bia 5 Gy), liều đo được ở các collimator 18 mm, 14 mm và 8 mm lần lượt là 9,90 Gy (sai lệch 0,96%), 9,81 Gy (sai lệch 1,86%) và 9,79 Gy (sai lệch 2,12%). Cả 3 kích thước này đều đáp ứng xuất sắc yêu cầu lâm sàng với sai số dưới 5%. Riêng trường chiếu 4 mm ghi nhận liều 8,89 Gy (sai lệch 11,13%) do hiệu ứng trường bức xạ siêu nhỏ làm giảm khả năng cân bằng điện tử trong buồng ion hóa.
  • Phân bố liều và đường đồng liều trên phim Gafchromic RTQA: Kết quả chiếu xạ phim thực nghiệm ghi nhận độ rộng trung bình của đường đồng liều 50% tương ứng là 5,7 mm (collimator 4 mm), 9,8 mm (collimator 8 mm), 13,5 mm (collimator 14 mm) và 22,8 mm (collimator 18 mm). Tỷ lệ vùng tối chuyển tiếp liều nhỏ giúp giảm thiểu liều bức xạ chiếu vào mô não lành xung quanh khối u.

Thảo luận kết quả

Khi phân tích dữ liệu phân bố liều, các kết quả có thể được trực quan hóa thông qua biểu đồ biến thiên liều dọc trục không gian 3D (Ox, Oy, Oz) và bảng so sánh độ rộng FWHM giữa phần mềm mô phỏng MCNP5 với hệ thống lập kế hoạch thương mại Osirix. Dữ liệu biểu đồ profile liều cho thấy độ dốc phân bố liều giảm cực kỳ nhanh chóng bên ngoài thể tích đích điều trị, khẳng định ưu thế vượt trội của cơ chế 30 chùm tia xoay.

Nguyên nhân dẫn đến sự sai lệch liều trên 10% tại trường chiếu 4 mm là do kích thước buồng ion hóa A1SL (đường kính trong 4,1 mm) lớn hơn kích thước trường bức xạ thực tế, dẫn đến hiện tượng lấy tích phân trung bình thể tích làm suy giảm giá trị đo. Điều này hoàn toàn phù hợp với các công bố quốc tế của Joel Y.C Cheung và Al-Dweiri khi nghiên cứu trường chiếu nhỏ trên máy Leksell Gamma Knife.

So sánh với các nghiên cứu mô phỏng bằng code EGS4 hoặc PENELOPE, chương trình MCNP5 trong luận văn này đã tái hiện chính xác cấu trúc hình học phức tạp của 30 kênh chuẩn trực vonfram xoay. Sai số thời gian phát xạ chỉ 1 giây trên tổng chu kỳ 600 giây (sai số 0,16%), và suất liều rò rỉ bề mặt đạt dưới 10 microSv/h, chứng minh hệ thống RGS vận hành hoàn toàn ổn định và an toàn cho quy trình xạ trị lâm sàng.

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

Dựa trên các kết quả đo đạc thực nghiệm và dữ liệu mô phỏng Monte Carlo, nghiên cứu đưa ra 4 giải pháp thực tiễn nhằm hoàn thiện quy trình xạ phẫu:

  • Chuẩn hóa quy trình QA định kỳ cho hệ thống cơ khí và định vị ART: Tiến hành kiểm chuẩn tọa độ không gian 3D với tần suất 1 lần/tuần, duy trì mục tiêu sai số dịch chuyển cơ học dưới 0,1 mm. Đơn vị thực hiện: Kỹ sư Vật lý Y khoa phối hợp với bộ phận kỹ thuật thiết bị xạ trị.
  • Ứng dụng hệ số hiệu chỉnh trường chiếu nhỏ cho collimator 4 mm: Áp dụng bổ sung hệ số hiệu chỉnh thể tích buồng đo hoặc chuyển sang sử dụng detector diode silic/phim xạ trị chuyên dụng để kiểm soát sai số liều dưới 3% trước khi lập kế hoạch cho các tổn thương siêu nhỏ. Timeline áp dụng: Ngay trong quý đầu tiên của chu kỳ bảo dưỡng hệ thống. Chủ thể thực hiện: Nhóm Vật lý Y học hạt nhân.
  • Tích hợp code mô phỏng MCNP5 vào hệ thống kiểm tra chéo kế hoạch điều trị: Sử dụng MCNP5 như một công cụ thẩm định độc lập song song với phần mềm Osirix đối với các ca bệnh u não có hình thái bất thường hoặc gần các cơ quan trọng yếu như giao thoa thị giác và thân não. Thời gian hoàn thành thử nghiệm: Trong vòng 6 tháng. Chủ thể: Nhóm nghiên cứu vật lý nguyên tử và bác sĩ xạ phẫu.
  • Tái nạp hoặc hiệu chỉnh chu kỳ nguồn phóng xạ Cobalt-60: Lập lịch trình bù trừ thời gian phát xạ tự động theo chu kỳ bán rã 5,26 năm của Cobalt-60, đảm bảo suất liều điều trị thực tế không suy giảm quá 50% sau mỗi 5 năm vận hành. Chủ thể thực hiện: Ban lãnh đạo Trung tâm Ung bướu và Kỹ sư trưởng an toàn bức xạ.

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

  • Kỹ sư Vật lý Y khoa và Cán bộ An toàn bức xạ: Tiếp cận quy trình chi tiết về kiểm chuẩn chất lượng (QA/QC), phương pháp chuyển đổi giá trị đọc trên Electrometer, tính toán hệ số kTP và kỹ thuật chuẩn liều tuyệt đối bằng buồng ion hóa A1SL.
  • Bác sĩ chuyên khoa Ung bướu và Xạ phẫu thần kinh: Nắm vững nguyên lý phân bố đường đồng liều 50%–90%, đặc tính suy giảm liều ngoài rìa khối u để tối ưu hóa chỉ định liều và bảo tồn tối đa các cấu trúc thần kinh lành lanh quanh tổn thương.
  • Học viên cao học và Nghiên cứu sinh ngành Vật lý hạt nhân / Vật lý nguyên tử: Sử dụng tài liệu như một cẩm nang mẫu về phương pháp xây dựng input file MCNP5, định nghĩa mặt cong (surface card), ô mạng (cell card), tally năng lượng F6/F8 và kỹ thuật xử lý phổ photon Cobalt-60.
  • Đơn vị quản lý trang thiết bị y tế và Bệnh viện chuyên khoa: Tham khảo các tiêu chuẩn kỹ thuật kiểm định máy xạ phẫu dao gamma quay, lập kế hoạch bảo trì phòng hộ bức xạ và đánh giá hiệu năng kinh tế - kỹ thuật của thiết bị RGS.

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

Dao gamma quay RGS có điểm gì vượt trội so với dao gamma Leksell cổ điển?

Hệ thống RGS sử dụng 30 nguồn Cobalt-60 xoay đồng bộ thay vì 201 nguồn tĩnh như hệ thống Leksell cũ. Cải tiến này giúp giảm đáng kể chi phí nạp nguồn, đồng thời tạo ra 30 cung chiếu xoay 360 độ giúp phân tán liều bức xạ qua thể tích mô lành lớn hơn, giảm nguy cơ tổn thương não lành và sai số cơ học luôn dưới 0,3 mm.

Tại sao cần phải thực hiện quy trình QA định kỳ trước khi xạ phẫu?

Quy trình QA đảm bảo tính chính xác tuyệt đối về mặt cơ học và liều lượng bức xạ. Với nguồn Cobalt-60 có chu kỳ bán rã 5,26 năm, suất liều thay đổi liên tục theo thời gian. QA giúp kiểm soát sai số đồng tâm dưới 0,5 mm và sai lệch liều thực tế dưới 5%, đảm bảo tiêu diệt khối u chính xác mà không gây biến chứng cho bệnh nhân.

Tại sao có sự sai lệch liều lên tới 11,13% khi đo ở ống chuẩn trực 4 mm?

Sự sai lệch này xuất phát từ hiệu ứng trường chiếu nhỏ (small field dosimetry). Buồng ion hóa A1SL có đường kính 4,1 mm, lớn hơn đường kính chùm tia 4 mm, dẫn đến hiện tượng thiếu hụt cân bằng điện tử tích cực và tích phân liều không đầy đủ. Điều này được khắc phục bằng cách sử dụng phim Gafchromic hoặc hệ số hiệu chỉnh đặc thù.

Chương trình mô phỏng Monte Carlo MCNP5 đóng vai trò gì trong xạ phẫu?

MCNP5 là chương trình mô phỏng số giải bài toán vận chuyển hạt photon và electron ngẫu nhiên trong không gian 3D. Phần mềm cho phép tính toán chính xác phổ năng lượng photon (1,17 MeV và 1,33 MeV) và phân bố liều trong môi trường không đồng nhất, đóng vai trò công cụ kiểm chứng độc lập cho phần mềm lập kế hoạch điều trị Osirix.

Mức độ an toàn bức xạ xung quanh phòng máy RGS được kiểm soát như thế nào?

Hệ thống RGS được trang bị cơ cấu che chắn kép gồm thân máy bằng gang dày và cửa đóng mở dạng hàm cá sấu. Kết quả kiểm tra rò rỉ bức xạ thực tế cho thấy suất liều bề mặt vỏ máy luôn duy trì dưới 10 microSv/h, hoàn toàn tuân thủ các quy chuẩn an toàn bức xạ nghiêm ngặt của IAEA và Bộ Y tế Việt Nam.

Kết luận

  • Luận văn đã xây dựng thành công quy trình đảm bảo chất lượng (QA) toàn diện cho hệ thống xạ phẫu dao gamma quay GammaART-6000TM tại Bệnh viện Bạch Mai, kiểm soát sai số hình học dưới 0,1 mm.
  • Thực nghiệm đo liều hấp thụ khẳng định độ chính xác cao với sai lệch liều tại tâm chỉ từ 0,96% đến 2,12% đối với các ống chuẩn trực tiêu chuẩn từ 8 mm đến 18 mm.
  • Ứng dụng thành công chương trình MCNP5 mô phỏng chính xác tương tác chùm tia Cobalt-60 và phân bố liều không gian 3D trong phantom mô phỏng.
  • Đã phân tích chi tiết hiệu ứng trường chiếu nhỏ trên collimator 4 mm, cung cấp cơ sở dữ liệu quan trọng để hiệu chỉnh kế hoạch xạ phẫu lâm sàng.
  • Thiết lập quy trình kiểm tra chéo độc lập giữa mô phỏng Monte Carlo và phần mềm Osirix, nâng cao chất lượng điều trị ung bướu nội sọ.

Kế hoạch tiếp theo định hướng mở rộng mô phỏng trên mô hình phantom đầu người không đồng nhất (kết hợp dữ liệu voxel CT/MRI thực tế) trong vòng 12 tháng tới. Các trung tâm xạ trị và nhà nghiên cứu quan tâm có thể khai thác toàn bộ tham số kỹ thuật và file cấu hình mô phỏng chi tiết trong tài liệu luận văn để chuẩn hóa quy trình điều trị tại cơ sở.