Tổng quan nghiên cứu

Theo thống kê từ GLOBOCAN năm 2012, toàn cầu ghi nhận hơn 7,149 triệu ca ung thư mới ở nam giới và hơn 4,030 triệu ca tử vong, trong đó Việt Nam có khoảng 125.000 ca mắc mới cùng 94.000 trường hợp tử vong mỗi năm. Nguy cơ mắc ung thư trước 75 tuổi tại Việt Nam ở mức 14,5%, đặt ra yêu cầu cấp thiết về việc nâng cao hiệu quả các kỹ thuật điều trị. Xạ trị truyền thống sử dụng chùm photon gặp nhược điểm lớn khi liều lượng bức xạ suy giảm theo hàm mũ, gây tổn thương diện rộng cho các tổ chức mô lành xung quanh khối u. Ngược lại, xạ trị proton mang lại bước đột phá nhờ khả năng tập trung năng lượng tối đa tại cuối quãng chạy (đỉnh Bragg), giúp tỷ lệ kiểm soát u sau 3 năm đạt trên 90% đối với ung thư phổi không tế bào nhỏ và xấp xỉ 100% đối với ung thư tiền liệt tuyến.

Nghiên cứu được thực hiện với mục tiêu làm rõ các đặc trưng tương tác vật lý - sinh học của hạt proton, đồng thời ứng dụng gói chương trình PENH trong bộ mã nguồn PENELOPE để mô phỏng phân bố liều lượng sâu của chùm proton trong các loại phantom mô phỏng sinh học ở dải năng lượng từ 100 MeV đến 250 MeV. Đề tài được triển khai trong phạm vi thời gian từ ngày 18/01/2018 đến ngày 14/07/2018 tại Khoa Khoa học Ứng dụng, Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh. Kết quả nghiên cứu đóng vai trò nền tảng vật lý tính toán, trực tiếp hỗ trợ việc tối ưu hóa kế hoạch xạ trị (Treatment Planning System) và đánh giá độ an toàn liều lượng trước khi Việt Nam đưa vào vận hành các trung tâm xạ trị proton và hạt nặng hiện đại.

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

Khung lý thuyết áp dụng

Khung lý thuyết của đề tài được xây dựng dựa trên các định luật vật lý hạt nhân và nguyên lý đo liều bức xạ y khoa:

  • Công thức dừng Bethe-Bloch: Xác định công suất dừng khối lượng dE/(ρdx) khi proton mang điện tích dương và khối lượng nghỉ 938 MeV/c2 (gấp 1836 lần electron) tương tác ion hóa với các electron nguyên tử trong vật chất. Công suất dừng tỷ lệ nghịch với bình phương vận tốc hạt, giải thích sự gia tăng đột biến của liều hấp thụ khi hạt chậm lại ở cuối quãng đường đi.
  • Lý thuyết đa tán xạ Coulomb (Multiple Coulomb Scattering - MCS): Mô tả độ lệch góc của chùm hạt theo phân bố Gauss khi xuyên qua lớp vật chất chắn, ảnh hưởng trực tiếp đến độ mở rộng biên trường chiếu với góc tán xạ của proton 200 MeV đạt khoảng 8,2x10^-4 radian.
  • Lý thuyết tương tác hạt nhân không đàn hồi: Phân tích cơ chế hấp thụ hạt nhân và tán xạ sinh hạt thứ cấp (neutron, proton thứ cấp), chiếm xác suất khoảng 1% ở mức năng lượng 25 MeV và tăng lên tới 25% ở mức năng lượng 200 MeV.
  • Hiệu ứng sinh học tương đối (RBE): Trong thực hành lâm sàng, tỷ số RBE giữa liều chuẩn photon và liều proton được thiết lập cố định ở mức 1,1, đồng thời kỹ thuật mở rộng đỉnh Bragg (Spread-Out Bragg Peak - SOBP) được ứng dụng để tạo trường liều phẳng đồng nhất trên toàn bộ thể tích bia u.

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

Nghiên cứu sử dụng phương pháp mô phỏng ngẫu nhiên Monte Carlo – tiêu chuẩn vàng trong tính toán đo liều xạ trị hiện đại. Dữ liệu hình học và vật liệu được thiết lập chi tiết trên 4 mô hình phantom đồng nhất gồm: phantom nước chuẩn (mật độ 1,00 g/cm3), phantom mô mềm (mật độ 1,05 g/cm3), phantom cơ vân (mật độ 1,05 g/cm3) và phantom xương đặc (mật độ 1,61 g/cm3). Các giá trị mật độ điện tử tương đối và thành phần khối lượng được chuẩn hóa theo dữ liệu chụp cắt lớp vi tính (CT).

Cỡ mẫu mô phỏng được cấu hình từ 10^6 đến 10^7 lịch sử hạt (histories) cho mỗi kịch bản năng lượng (100 MeV, 150 MeV, 250 MeV) và khoảng cách nguồn chiếu (từ 0 cm đến 40 cm). Phương pháp lấy mẫu ngẫu nhiên theo phân bố xác suất giúp tái hiện chân thực từng bước va chạm cứng và va chạm mềm của hạt. Thuật toán mô phỏng được thực thi trên gói phần mềm PENELOPE phiên bản mở rộng PENH chuyên dụng cho proton, sau đó toàn bộ ma trận dữ liệu phân bố liều được xử lý và phân tích thống kê thông qua phần mềm Origin. Lý do lựa chọn PENH là nhờ khả năng kết hợp mô hình tán xạ Eikonal chuẩn xác, vừa tối ưu hóa thời gian tính toán vừa đảm bảo sai số thống kê dưới 1,5%. Toàn bộ quy trình mô phỏng và xử lý số liệu được tiến hành liên tục trong 6 tháng (từ tháng 1/2018 đến tháng 7/2018).

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

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

Nghiên cứu đã xác lập các đặc tính phân bố liều sâu của chùm tia proton với những số liệu cụ thể:

  • Sự phụ thuộc của độ sâu đỉnh Bragg vào mức năng lượng chùm tia: Trong phantom nước, khi tăng năng lượng từ 100 MeV lên 150 MeV và 250 MeV, vị trí đỉnh Bragg dịch chuyển sâu dần từ 7,7 cm lên 15,8 cm và đạt cực đại ở 37,9 cm. Độ sâu xâm nhập ở mức 250 MeV tăng hơn 392% so với mức năng lượng 100 MeV, chứng minh tính linh hoạt trong việc kiểm soát độ sâu khối u.
  • Ảnh hưởng rõ rệt của mật độ vật liệu phantom: Cùng ở mức năng lượng 100 MeV, vị trí đỉnh Bragg trong phantom xương đặc chỉ đạt độ sâu 5,2 cm, ngắn hơn 32,5% so với phantom nước (7,7 cm) do mật độ khối của xương cao hơn đáng kể (1,61 g/cm3 so với 1,00 g/cm3). Trong khi đó, phantom mô mềm và phantom cơ vân có vị trí đỉnh Bragg tương đương nhau, lần lượt đạt 7,5 cm và 7,4 cm.
  • Đặc tính liều tại lối vào và tỷ lệ đỉnh liều: Tại mức năng lượng 200 MeV, chiều rộng nửa cực đại (FWHM) của đỉnh Bragg đạt khoảng 22 mm và liều giảm nhanh từ 90% xuống 10% chỉ trong khoảng cách 7 mm. Tỷ số giữa liều cực đại tại đỉnh Bragg và liều lối vào bề mặt dao động từ 3:1 đến 4:1.
  • Độ ổn định theo khoảng cách nguồn - phantom (SSD): Khi tăng khoảng cách từ 0 cm lên 10 cm, 20 cm và 40 cm ở năng lượng 250 MeV, độ sâu của đỉnh Bragg giữ nguyên ở mức 37,9 cm với sai lệch tọa độ dưới 0,5 mm, trong khi liều bề mặt giảm nhẹ từ 2% đến 5% do độ phân kỳ tự nhiên của chùm tia trong không khí.

Thảo luận kết quả

Các kết quả mô phỏng bằng PENH hoàn toàn tương thích với dữ liệu chuẩn của E. Sterpin và các cộng sự, với độ sai lệch vị trí đỉnh Bragg trong phantom nước dưới 1,2% ở cả hai mốc 100 MeV và 250 MeV. Cơ chế vật lý cốt lõi nằm ở công suất dừng tuyến tính Bethe-Bloch: khi động năng của proton giảm dần, tiết diện tương tác ion hóa tăng vọt, dẫn đến việc hơn 70% tổng năng lượng bị mất được giải phóng cục bộ tại vùng hẹp cuối hành trình.

Dữ liệu mô phỏng được biểu diễn trực quan thông qua hệ thống đồ thị đường cong Bragg (Dose-Depth Curves) và biểu đồ liều khối (Dose-Volume Histogram - DVH). Khác biệt căn bản giữa xạ trị photon và proton thể hiện rõ ở liều lối ra: chùm photon giải phóng liều cực đại ngay gần bề mặt rồi suy giảm chậm dần ra phía sau cơ thể, trong khi chùm proton triệt tiêu liều về 0 ngay sau đỉnh Bragg. Điều này mang ý nghĩa sinh học then chốt, cho phép tăng liều tiêu diệt khối u đích thêm 10% đến 20% mà vẫn bảo vệ tối đa các cơ quan lành nhạy cảm liền kề (như tủy sống, thân não, ốc tai hoặc tim).

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

Dựa trên kết quả nghiên cứu mô phỏng, 4 khuyến nghị chiến lược được đề xuất cho ngành vật lý y khoa và xạ trị ung bướu tại Việt Nam:

  1. Tích hợp giải thuật Monte Carlo PENH vào hệ thống lập kế hoạch xạ trị (TPS): Các bệnh viện chuyên khoa ung bướu và viện nghiên cứu cần ứng dụng module PENH để kiểm chuẩn độc lập các kế hoạch xạ trị proton phức tạp, mục tiêu đưa sai số liều tính toán thực tế xuống dưới 2% trong giai đoạn 2024 - 2026.
  2. Hoàn thiện đề án xây dựng Trung tâm Xạ trị Proton quốc gia: Ban Quản lý dự án của Bộ Y tế và Bệnh viện K cần đẩy nhanh việc đầu tư hệ thống máy gia tốc Cyclotron hoặc Synchrotron năng lượng 250 MeV, hướng tới mục tiêu điều trị thành công cho 500 đến 1.000 bệnh nhân ung thư nhi khoa và u nền sọ mỗi năm bắt đầu từ năm 2028.
  3. Mở rộng mô hình tính toán tương tác hạt nhân thứ cấp: Nhóm nghiên cứu vật lý hạt nhân cần phát triển thêm các thuật toán mô phỏng chi tiết dòng neutron tán xạ và electron thứ cấp năng lượng thấp (dưới 10 eV), giúp kiểm soát suất liều rò rỉ dưới 0,5 mGy trên mỗi Gray liều điều trị từ nay đến năm 2027.
  4. Đào tạo chuẩn hóa nhân lực Kỹ sư Vật lý Y khoa (Medical Physicist): Trường Đại học Bách Khoa và các trường đại học khối khoa học tự nhiên cần đẩy mạnh chương trình đào tạo chuyên sâu về xạ trị ion nặng và phần mềm mô phỏng hạt, phấn đấu cung cấp tối thiểu 50 kỹ sư đạt chứng chỉ đo liều lâm sàng quốc tế trong vòng 5 năm tới.

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

Nội dung luận văn mang lại giá trị học thuật và thực tiễn sâu sắc cho 4 nhóm đối tượng chính:

  • Kỹ sư Vật lý Y khoa (Medical Physicists): Nắm bắt phương pháp tính toán và thiết lập các thông số chùm tia, khai thác bộ công cụ Monte Carlo PENELOPE/PENH để kiểm chuẩn chất lượng máy gia tốc và tối ưu hóa phân bố liều 3D trên mô hình bệnh nhân.
  • Bác sĩ chuyên khoa Xạ trị Ung bướu: Hiểu rõ cơ sở vật lý của đỉnh Bragg, hiệu suất sinh học RBE = 1,1 và kỹ thuật SOBP nhằm đưa ra chỉ định lâm sàng chuẩn xác cho các khối u nằm sâu hoặc kháng tia xạ thông thường.
  • Học viên cao học và Nghiên cứu sinh ngành Vật lý Kỹ thuật, Kỹ thuật Hạt nhân: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về tương tác bức xạ với vật chất, kỹ thuật mô phỏng Monte Carlo và phương pháp xử lý dữ liệu vật lý bằng phần mềm Origin.
  • Nhà quản lý y tế và Chuyên gia Quy hoạch bệnh viện: Tham khảo các tiêu chuẩn kỹ thuật về máy gia tốc proton (Cyclotron, Synchrotron, Linac, ECRIPAC), hệ thống vận chuyển chùm tia và quy chuẩn che chắn bức xạ neutron phục vụ thẩm định các dự án đầu tư y tế công nghệ cao.

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

  • Ưu điểm vượt trội nhất của xạ trị proton so với xạ trị photon truyền thống là gì? Xạ trị proton sở hữu đặc tính đỉnh Bragg độc nhất, cho phép chùm hạt truyền năng lượng thấp ở lối vào, giải phóng liều cực đại tại vị trí khối u và triệt tiêu liều về 0 ngay sau đó. Ngược lại, photon suy giảm liều theo hàm mũ gây tổn thương kéo dài cho các mô lành phía sau khối u.

  • Tại sao chương trình PENELOPE với module PENH lại phù hợp cho mô phỏng liều proton? PENELOPE kết hợp thuật toán mô phỏng hỗn hợp giúp theo dõi chính xác tương tác của các hạt mang điện. Module PENH được bổ sung mô hình tán xạ Eikonal giúp mô tả hoàn hảo va chạm đàn hồi và ion hóa của proton trong môi trường vật chất với độ lệch vị trí đỉnh Bragg dưới 1,2%.

  • Đỉnh Bragg trải rộng (SOBP) có ý nghĩa như thế nào trong xạ trị lâm sàng? Đỉnh Bragg tự nhiên có bề rộng rất hẹp chỉ khoảng vài milimet đến 22 mm, không đủ bao phủ khối u lớn. Kỹ thuật SOBP chồng chập nhiều chùm proton ở các mức năng lượng khác nhau để tạo ra một vùng liều cao đồng nhất bao trọn toàn bộ thể tích hình học của khối u.

  • Mật độ của các loại mô sinh học ảnh hưởng ra sao đến quãng chạy của proton? Quãng chạy proton tỷ lệ nghịch với mật độ điện tử của môi trường. Ở cùng mức năng lượng 100 MeV, proton chỉ xuyên sâu 5,2 cm trong xương đặc có mật độ 1,61 g/cm3, trong khi có thể xuyên sâu tới 7,7 cm trong phantom nước chuẩn có mật độ 1,00 g/cm3.

  • Tại sao giá trị hiệu ứng sinh học tương đối (RBE) được cố định là 1,1 trong xạ trị proton? Giá trị RBE = 1,1 được chuẩn hóa quốc tế dựa trên các thử nghiệm lâm sàng, biểu thị rằng chùm proton có hiệu quả tiêu diệt tế bào cao hơn 10% so với liều photon tương đương. Việc cố định giá trị này giúp các bác sĩ dễ dàng chuyển đổi đơn vị liều từ phác đồ xạ trị photon truyền thống.

Kết luận

  • Luận văn đã xây dựng thành công quy trình ứng dụng chương trình PENELOPE và module PENH để mô phỏng chính xác phân bố liều sâu của chùm proton ở các mức năng lượng 100 MeV, 150 MeV và 250 MeV.
  • Xác thực độ tin cậy vượt trội của mô hình mô phỏng với độ sai lệch vị trí đỉnh liều dưới 1,2% so với dữ liệu công bố quốc tế của E. Sterpin trên phantom nước chuẩn.
  • Định lượng hóa ảnh hưởng của mật độ khối và mật độ điện tử trên 4 loại mô sinh học (nước, mô mềm, cơ vân, xương đặc), chứng minh quãng chạy giảm tới 32,5% khi đi qua mô xương đặc.
  • Thiết lập cơ sở dữ liệu vật lý y khoa quan trọng phục vụ công tác lập kế hoạch điều trị và chuẩn bị hạ tầng kỹ thuật cho các dự án trung tâm xạ trị proton tại Việt Nam giai đoạn 2025 - 2030.
  • Các đơn vị y tế và viện nghiên cứu cần tiếp tục mở rộng mô hình mô phỏng trên dữ liệu ảnh cắt lớp CT thực tế của bệnh nhân nhằm chuẩn hóa tối đa quy trình xạ trị công nghệ cao trong tương lai gần.