Tổng quan về luận án

Y học hạt nhân hiện đại đóng vai trò thiết yếu trong việc chẩn đoán và điều trị ung bướu thông qua các dược chất phóng xạ phát bức xạ ion hóa. Theo nguyên lý chuyển hóa đồng vị do George de Hevesy đề xuất, các đồng vị phóng xạ tham gia vào các quá trình sinh hóa tương tự như các nguyên tố đồng vị bền, cho phép đánh giá chức năng chuyển hóa ở mức tế bào với hoạt độ chẩn đoán khoảng 50–100 µCi hoặc tiêu diệt tổ chức bệnh lý với hoạt độ điều trị lên đến 250 mCi. Tuy nhiên, rủi ro sinh học bức xạ luôn hiện hữu theo mô hình tác dụng ngẫu nhiên không có ngưỡng liều an toàn; theo báo cáo của Ủy ban BEIR (National Research Council), khi 100.000 cá thể cùng nhận liều chiếu gamma 0,1 Gy trong một lần duy nhất, ước tính sẽ có thêm 800 ca tử vong do ung thư phát sinh. Do đó, yêu cầu tối thượng trong lâm sàng và an toàn bức xạ là phải tính toán chuẩn xác phân bố liều hấp thụ, nhằm tối đa hóa hiệu quả phá hủy khối u, giảm thiểu biến chứng tại mô lành và bảo vệ những người xung quanh theo nguyên tắc ALARA.

       +-------------------------------------------------------------------------+
       |             THỰC TRẠNG TÍNH LIỀU CHIẾU TRONG TRUYỀN THỐNG               |
       |  - Sơ đồ MIRD cơ quan: Giả định mô đồng nhất, bỏ qua giải phẫu cá thể.  |
       |  - Phantom Voxel cổ điển: Sai số thể tích biên (Partial Volume Effect). |
       |  - Ước tính liều tiếp xúc: Dùng nguồn điểm/hình trụ đơn giản hóa.       |
       +-------------------------------------------------------------------------+
                                            |
                                            v
       +-------------------------------------------------------------------------+
       |             GIẢI PHÁP ĐỘT PHÁ CỦA LUẬN ÁN (GAMOS/GEANT4)                |
       |  1. Kỹ thuật "Hình học song song" (Parallel Geometry):                  |
       |     Chồng khớp đường biên RTSTRUCT DICOM lên lưới voxel -> Triệt tiêu   |
       |     sai số thể tích biên ở cấu trúc mỏng/phức tạp (tuyến giáp, da).    |
       |  2. Kỹ thuật "Ghép Phantom" & Xoay Vector (Phantom Merging & Rotation): |
       |     Tích hợp phân bố hoạt độ thực (SPECT/PET) của bệnh nhân và phantom  |
       |     người tiếp xúc (ICRP 110) trong cùng một không gian mô phỏng.       |
       +-------------------------------------------------------------------------+
                                            |
                                            v
       +-------------------------------------------------------------------------+
       |                         KẾT QUẢ & ĐÓNG GÓP                              |
       |  - Khắc phục sai lệch liều tiếp xúc lên tới 366% so với nguồn điểm.     |
       |  - Tối ưu hóa điều trị I-131 và cung cấp cơ sở khoa học giải phóng ly   |
       |    cho 30% bệnh nhân chịu áp lực tâm lý không cần thiết.                |
       +-------------------------------------------------------------------------+

Khoảng trống nghiên cứu then chốt xuất phát từ hai hạn chế nền tảng trong đo liều hạt nhân đương đại: Thứ nhất, khi tính liều chiếu trong cá thể hóa bằng kỹ thuật mô phỏng Monte Carlo trên phantom voxel, các thuật toán phân đoạn truyền thống bắt buộc phải "voxel hóa" toàn bộ hình học cơ quan. Cơ chế nhị phân (chỉ xét tâm voxel nằm trong hay ngoài đường biên cấu trúc) dẫn đến sai số nghiêm trọng về thể tích và năng lượng lắng đọng tại các vùng biên mô mỏng, xoắn hoặc kích thước nhỏ như tuyến giáp, niêm mạc hay da. Thứ hai, trong đánh giá liều cho người tiếp xúc (nhân viên y tế, gia đình, công chúng), các công cụ quốc tế phần lớn quy ước bệnh nhân thành các mô hình nguồn điểm, nguồn đoạn thẳng hoặc hình trụ đồng nhất, làm sai lệch đáng kể trường bức xạ thực tế.

Luận án của nghiên cứu sinh Nguyễn Thị Phương Thảo, dưới sự hướng dẫn khoa học của TS. Nguyễn Đông Sơn và GS. TS. Pedro Arce Dubois (Trung tâm CIEMAT, Tây Ban Nha) tại Viện Năng lượng Nguyên tử Việt Nam (VINATOM), thuộc chuyên ngành Vật lý nguyên tử và hạt nhân (Mã số: 9 44.06), đã giải quyết triệt để hai thách thức trên thông qua ba câu hỏi nghiên cứu và giả thuyết cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để loại bỏ sai số thể tích biên trong tính toán liều hấp thụ cho cấu trúc giải phẫu phức tạp khi mô phỏng trên phantom voxel?
    Giả thuyết 1 (H1): Ứng dụng kỹ thuật "hình học song song" (Parallel Geometry) cho phép tính toán chính xác phần thể tích thực tế và năng lượng hấp thụ của từng voxel bị cắt bởi đường biên RTSTRUCT DICOM mà không làm biến dạng lưới vật liệu ban đầu.
  2. Câu hỏi nghiên cứu 2 (RQ2): Mô hình phân bố hoạt độ từ ảnh y học hạt nhân (SPECT/PET) kết hợp kỹ thuật ghép hai phantom voxel ảnh hưởng như thế nào đến độ chính xác của liều tương đương và liều hiệu dụng ở người tiếp xúc?
    Giả thuyết 2 (H2): Phân bố hoạt độ không đồng nhất từ cơ thể bệnh nhân tạo ra trường bức xạ bất đối xứng; việc tích hợp hai phantom voxel đầy đủ trong không gian xoay vector 3D sẽ mang lại kết quả liều sai khác có ý nghĩa thống kê so với các mô hình nguồn điểm và nguồn hình học đơn giản.
  3. Câu hỏi nghiên cứu 3 (RQ3): Mức độ khác biệt liều tiếp xúc thực tế giữa các tư thế giao tiếp lâm sàng (đối diện, bên cạnh, trước sau, đứng - nằm) có đủ lớn để thay đổi phác đồ cách ly y tế hay không?
    Giả thuyết 3 (H3): Liều chiếu thực tế ở khoảng cách gần nhạy cảm cao với hướng tương đối của phân bố nguồn và tư thế tiếp xúc, cho phép tối ưu hóa quy định cách ly mà không xâm phạm giới hạn liều an toàn của IAEA/ICRP.

Khung lý thuyết của công trình tích hợp cơ học lượng tử bức xạ, lý thuyết vận chuyển hạt Monte Carlo, mô hình vật lý điện từ năng lượng thấp Livermore trong Geant4 và hệ thống đại lượng bảo vệ bức xạ của ICRP 103. Đột phá mang tính định lượng của luận án thể hiện ở việc khắc phục độ chênh lệch liều tiếp xúc lên tới 366% vốn bị thổi phồng bởi mô hình nguồn điểm ở khoảng cách gần (0,3 m), đồng thời cung cấp công cụ tính liều song song chính xác cho điều trị ung thư tuyến giáp bằng I-131, trực tiếp tác động đến quy trình lâm sàng và tâm lý của khoảng 30% bệnh nhân vốn gặp hội chứng trầm cảm do cách ly quá mức.

Literature Review và Positioning

Lịch sử phát triển của kỹ thuật đo liều chiếu trong y học hạt nhân được định hình rõ nét qua sự dịch chuyển mô thức từ cấp độ cơ quan vĩ mô sang cấp độ voxel vi mô:

+-----------------------------------------------------------------------------------+
| NIÊN BIỂU TIẾN HÓA CỦA CÁC PHƯƠNG PHÁP TÍNH LIỀU CHIẾU TRONG TRONG Y HỌC HẠT NHÂN  |
+-----------------------------------------------------------------------------------+
|  1956 - 1965: Loevinger, Ellett et al.                                            |
|  * Đặt nền móng tỉ lệ hấp thụ bức xạ; ứng dụng Monte Carlo sơ khai.               |
|  -------------------------------------------------------------------------------  |
|  1968 - 1990s: Snyder et al., Loevinger et al. (Ủy ban MIRD)                     |
|  * Chuẩn hóa sơ đồ MIRD cấp cơ quan trên Phantom toán học (MIRD-5).               |
|  * Giả định: Cơ quan hình học đơn giản, hoạt độ và mô đồng nhất.                  |
|  -------------------------------------------------------------------------------  |
|  1999: Bolch, Bouchet et al. (MIRD Pamphlet 17)                                   |
|  * Đề xuất tính liều cấp Voxel dựa trên Voxel S-Values (VSV).                     |
|  * Hạn chế: Giả định môi trường đồng nhất (nước/mô mềm), giới hạn xuyên sâu gamma.|
|  -------------------------------------------------------------------------------  |
|  2000s - 2010s: Xu et al., ICRP 110, Bardiès et al.                              |
|  * Mô phỏng Monte Carlo trực tiếp trên Phantom Voxel cá thể hóa (CT/NM).          |
|  * Thách thức: Sai số thể tích biên (Partial Volume Effect) ở cấu trúc nhỏ.       |
|  -------------------------------------------------------------------------------  |
|  2021: Luận án Nguyễn Thị Phương Thảo (GAMOS/Geant4)                              |
|  * Thiết lập kỹ thuật "Hình học song song" xử lý triệt để sai số biên RTSTRUCT.   |
|  * Phát triển thuật toán "Ghép Phantom & Xoay Vector" tính liều tiếp xúc đa tư thế.|
+-----------------------------------------------------------------------------------+

Giai đoạn đầu (1956–1968), Loevinger và Ellett đặt nền móng cho khái niệm "tỷ lệ hấp thụ" bằng mô phỏng Monte Carlo trên các bia hình học đơn giản, dẫn đến sự ra đời của sơ đồ chuẩn hóa MIRD (Medical Internal Radiation Dose). Sơ đồ MIRD cổ điển tính toán liều hấp thụ cơ quan đích $D_{\text{bia}}$ từ tích số hoạt độ tích lũy $\tilde{A}$ của cơ quan nguồn và giá trị $S$ (liều hấp thụ trên mỗi phân rã). Tuy nhiên, các phantom toán học mô tả cơ thể bằng các khối cầu, trụ hoặc elip cơ bản (Snyder et al., 1978) đã bộc lộ nhược điểm cố hữu: bỏ qua tính dị dưỡng sinh học, không phản ánh giải phẫu cá thể và bất lực trong việc mô tả khối u nằm ngoài vị trí chuẩn.

Sự bùng nổ của chụp cắt lớp vi tính (CT) bởi Hounsfield và Cormack, cùng kỹ thuật PET/SPECT, đã thúc đẩy Ủy ban MIRD ban hành Ấn phẩm MIRD 17 (Bolch et al., 1999), chuyển dịch tính liều sang mức voxel. Phương pháp tiếp cận mức voxel phân nhánh thành ba nhánh chính:

  • Phương pháp nhân liều (Dose Point Kernel - DPK): Đơn giản, tốc độ tính nhanh nhưng dựa trên giả định môi trường vô hạn và đồng nhất nước.
  • Phương pháp giá trị S voxel (Voxel S-Values - VSV): Tích chập ma trận phân bố hoạt độ với giá trị S được tính toán trước. Dù Dieudonné và cộng sự (2013) đã đề xuất các hệ số hiệu chỉnh mật độ cho vùng bụng, phương pháp này thất bại khi áp dụng tại vùng ngực, vòm họng và tủy xương do giả định giới hạn quãng chạy gamma dưới 5 voxel.
  • Phương pháp mô phỏng Monte Carlo trực tiếp (Direct MC): Cho phép đưa trực tiếp ma trận số Hounsfield (HU) từ CT và hoạt độ từ ảnh SPECT/PET vào không gian tính toán.

Mặc dù phương pháp Monte Carlo trực tiếp là tiêu chuẩn vàng, nó vấp phải cuộc tranh luận gay gắt về hiện tượng sai số thể tích biên (Partial Volume Artifacts). Khi chuyển đổi từ đường phân đoạn cấu trúc giải phẫu (RTSTRUCT DICOM do bác sĩ vẽ) sang ma trận voxel mm³, một tỷ lệ lớn voxel ở biên bị cắt xén. Các nghiên cứu quốc tế đã đề xuất nhiều giải pháp khắc phục nhưng đều bộc lộ giới hạn: nghiên cứu của Lee và cộng sự (2007) chỉ định cố định bề dày da 1,5 mm; Pafundi và cộng sự (2010) phát triển thuật toán phân đoạn xương chuyên biệt dựa trên mô hình đa giác 3D nhưng không thể tổng quát hóa cho các tạng mềm; Cassola và cộng sự (2010) phát triển mesh phantom nhưng cấu trúc lưới đa diện đòi hỏi tài nguyên tính toán khổng lồ và không tương thích trực tiếp với dữ liệu DICOM lâm sàng.

Trong mảng tính liều cho người tiếp xúc, các nghiên cứu quốc tế tồn tại nhiều xung đột về mô hình nguồn:

  • Nghiên cứu của Siegel và cộng sự (2002) cùng Smith và cộng sự (2000) bảo vệ mô hình nguồn điểm và nguồn đoạn thẳng dựa trên định luật nghịch đảo bình phương vì tính tiện dụng lâm sàng.
  • Ngược lại, Zˇontar và cộng sự (2008) cùng Palmer và cộng sự (2014) chỉ ra rằng nguồn điểm thổi phồng liều hấp thụ lên đáng kể; Palmer chứng minh mô hình hình trụ tuyến giáp có thể làm tăng sai số ước tính liều lên 13,4% khi mô phỏng một bệnh nhi được bế bởi bệnh nhân điều trị I-131.
  • Nghiên cứu tiên phong của Albérico và cộng sự (2008) sử dụng phantom voxel FAX cho thấy mức chênh lệch liều giữa mô phỏng voxel và nguồn điểm lên tới 366% ở khoảng cách 0,3 m. Dù vậy, Albérico vẫn áp đặt giả định nguồn đồng nhất toàn tuyến giáp và bỏ qua phân bố hoạt độ thực tế từ ảnh SPECT. Dewji và cộng sự (2015) khảo sát các kịch bản tiếp xúc trên phương tiện công cộng nhưng chỉ áp dụng trên phantom toán học thô sơ.

Luận án này định vị chính xác vào giao điểm của các khoảng trống học thuật trên: lần đầu tiên thiết lập một phương pháp tính liều Monte Carlo kết hợp hình học song song để bảo toàn đường biên giải phẫu RTSTRUCT và phát triển thuật toán hợp nhất hai phantom voxel dị dưỡng độc lập trong hệ tọa độ xoay vector 3D, vượt qua các giới hạn của phần mềm GATE (Jan et al., 2004) hay OEDIPE (Lanconelli et al., 2012).

Đó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 trực tiếp Lý thuyết Vận chuyển Bức xạ Lượng tử (Radiation Transport Theory) và Sơ đồ Đo liều Chiếu trong MIRD thông qua việc tái định nghĩa hàm đóng góp năng lượng tại các phần tử thể tích dị biên:

                                  VOXEL BIÊN (VOXEL i)
                     +---------------------------------------------+
                     |  Thể tích ngoài (1 - f_i):                   |
                     |  Mô lành xung quanh                         |
                     |                                             |
                     |                 Đường biên RTSTRUCT         |
                     |                     /                       |
                     |        ~~~~~~~~~~~~/~~~~~~~~~~~~            |
                     |       (           /             )           |
                     |        (  f_i:   /             )            |
                     |         ( Khối lượng m_i,struct)            |
                     |          (Năng lượng E_i,struct)            |
                     |           ~~~~~~~~~~~~~~~~~~~~~             |
                     |  Thể tích trong (f_i): Cấu trúc giải phẫu   |
                     +---------------------------------------------+
                                            |
                                            v
              +-------------------------------------------------------------+
              |           CÔNG THỨC LIỀU HẤP THỤ MỚI CỦA CẤU TRÚC           |
              |                                                             |
              |       D_struct = [ \sum E_{i, struct} ] / M_struct          |
              |                                                             |
              |       * Triệt tiêu hoàn toàn nghịch lý Voxel Hóa nhị phân.  |
              |       * Bảo toàn bảo toàn năng lượng và phân bố khối lượng. |
              +-------------------------------------------------------------+

  1. Mở rộng công thức tính liều cấu trúc vi mô: Trong mô hình tính liều truyền thống, liều hấp thụ của một cấu trúc $D_{\text{struct}}$ được tính bằng tổng năng lượng của các voxel nguyên vẹn $\sum E_i$ chia cho tổng khối lượng $\sum m_i$, trong đó trạng thái của voxel $i$ nhận giá trị nhị phân ${0, 1}$ dựa trên vị trí tâm. Luận án thiết lập mô hình giải tích mới, xác định liều hấp thụ thông qua tích phân năng lượng phân đoạn: $$D_{\text{struct}} = \frac{\sum_{i=1}^{N_{\text{voxel}}} E_{i, \text{struct}}}{M_{\text{struct}}} = \frac{\sum_{i=1}^{N_{\text{voxel}}} E_{i, \text{struct}}}{\sum_{i=1}^{N_{\text{voxel}}} f_i \cdot \rho_i \cdot V_{\text{voxel}}}$$ Trong đó $f_i \in [0, 1]$ là tỷ lệ thể tích thực của voxel $i$ nằm bên trong đường phân đoạn giải phẫu, $\rho_i$ là mật độ khối lượng chuyển đổi từ giá trị HU tương ứng, và $E_{i, \text{struct}}$ là năng lượng bức xạ thực tế lắng đọng trong phần thể tích $f_i \cdot V_{\text{voxel}}$ được truy vết thông qua thế giới hình học song song.

  2. Hệ thống mệnh đề khoa học (Theoretical Propositions):

    • Mệnh đề 1 (P1): Năng lượng hấp thụ biên không tỷ lệ tuyến tính đơn thuần với thể tích hình học $f_i$ mà là một hàm phi tuyến phức tạp phụ thuộc vào phổ năng lượng electron/photon thứ cấp, mật độ dừng va chạm (stopping power) và vị trí tương đối của nguồn phát trong không gian lân cận.
    • Mệnh đề 2 (P2): Trường liều phát xạ từ bệnh nhân y học hạt nhân mang tính bất đối xứng không gian bậc cao do hiệu ứng tự hấp thụ và che chắn cấu trúc xương; do đó, liều hấp thụ của người tiếp xúc tại trường gần ($r < 1\text{ m}$) không tuân theo định luật tỉ lệ nghịch bình phương khoảng cách $1/r^2$.

  3. Chuyển dịch mô thức (Paradigm Shift): Chuyển đổi căn bản từ mô thức "Voxel hóa giải phẫu" (buộc cấu trúc sinh học uốn nắn theo lưới voxel trực giao) sang mô thức "Hình học song song liên tục" (bảo toàn nguyên vẹn cấu trúc giải phẫu trơn trên nền lưới vật liệu voxel hóa).

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

Khung phân tích của luận án tích hợp liên ngành ba trụ cột lý thuyết vững chắc:

  1. Lý thuyết tương tác photon-electron (Lý thuyết tán xạ Compton, Hiệu ứng Quang điện và Bức xạ Hãm Bremsstrahlung): Mô phỏng cơ chế truyền năng lượng liên tục (CSDA) của electron thông qua mô hình vật lý Livermore (Geant4), xử lý các tương tác năng lượng thấp xuống đến 250 eV.
  2. Lý thuyết Hình học Tính toán và Xử lý Ảnh Y khoa: Sử dụng cấu trúc hình học khối mở rộng (Extruded Solid) và thuật toán đa giác không gian lồng nhau xây dựng từ thư viện JTS (Java Topology Suite) để tái tạo thể tích 3D từ các lát cắt 2D của file DICOM RTSTRUCT.
  3. Định lý Xoay Không gian Euler và Phép Biến đổi Vector 3D: Thiết lập hệ phương trình chuyển đổi tọa độ động giữa hệ quy chiếu bệnh nhân ($O_P, X_P, Y_P, Z_P$) và hệ quy chiếu người tiếp xúc ($O_C, X_C, Y_C, Z_C$): $$\vec{r}_C = \mathbf{R}(\alpha, \beta, \gamma) \cdot (\vec{r}P - \vec{r}{\text{offset}})$$ Khung phân tích cho phép đặt hai cơ thể người chuẩn hóa (ICRP 110) hoặc cá thể hóa vào bất kỳ vị trí và góc nghiêng tương đối nào trong không gian thực nghiệm ảo.

Điều kiện biên của khung phân tích được xác định rõ: áp dụng cho các hạt phóng xạ beta và gamma trong dải năng lượng từ vài keV đến vài MeV; giả định tính bất biến hình học của các cơ quan trong suốt quá trình chụp ảnh và chiếu xạ (bỏ qua biến dạng do chuyển động hô hấp).

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế trên lập trường triết học thực chứng (Positivism) kết hợp chủ nghĩa duy thực phản biện (Critical Realism), xem xét các tương tác vật lý hạt nhân vi mô mang tính ngẫu nhiên khách quan có thể được mô hình hóa chính xác qua quy luật xác suất lớn Monte Carlo.

Thiết kế mô phỏng đa tầng (Multi-level Simulation Design) bao gồm hai cấp độ:

  • Cấp độ 1 (Bệnh nhân): Đánh giá vi mô phân bố liều trong tuyến giáp và mô lành lân cận trên bệnh nhân ung thư tuyến giáp tiếp nhận xạ trị I-131.
  • Cấp độ 2 (Người tiếp xúc): Đánh giá vĩ mô liều tương đương tại các cơ quan nhạy bức xạ ($H_T$) và liều hiệu dụng toàn thân ($E$) trên phantom chuẩn ICRP 110 giới tính nam (chiều cao 176 cm, cân nặng 73 kg, kích thước voxel $2,137 \times 2,137 \times 8,0\text{ mm}^3$).
                KIẾN TRÚC MÔ PHỎNG HÌNH HỌC SONG SONG & GHÉP PHANTOM
+-----------------------------------------------------------------------------------+
| THẾ GIỚI THỰC (MASS GEOMETRY)           | THẾ GIỚI SONG SONG (PARALLEL GEOMETRY)  |
| - Ma trận Voxel (từ CT Scan)             | - Tái tạo từ file RTSTRUCT DICOM        |
| - Chuyển đổi Schneider: HU -> rho -> Mat | - Khối Extruded Solid trơn nhẵn 3D      |
| - Quyết định tán xạ & truyền tương tác   | - Chấm điểm (Scoring) năng lượng chính  |
|                                         |   xác cho từng cấu trúc giải phẫu       |
+-----------------------------------------------------------------------------------+
                                            |
                                            v
+-----------------------------------------------------------------------------------+
|               KHÔNG GIAN GHÉP PHANTOM (PHANTOM MERGING & ROTATION)                |
|                                                                                   |
|   +--------------------------+                     +--------------------------+   |
|   |    PHANTOM BỆNH NHÂN     |                     |    PHANTOM TIẾP XÚC      |   |
|   |  - Dữ liệu ảnh CT        |                     |  - Phantom ICRP 110      |   |
|   |  - Nguồn: Ảnh SPECT/PET  |                     |  - Đo liều H_T và E      |   |
|   |    (Hoạt độ I-131 thực)  |                     |    (Trọng số ICRP 103)   |   |
|   +--------------------------+                     +--------------------------+   |
|                 \                                       /                         |
|                  \----[ Phép biến đổi Vector 3D ]------/                          |
|                       Khoảng cách d: 0,1m -> 2,0m                                 |
|                       Tư thế: Đối diện, Bên cạnh,                                 |
|                               Trước - Sau, Đứng - Nằm                             |
+-----------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

  1. Xử lý dữ liệu hình ảnh DICOM: Dữ liệu ảnh CT và SPECT của bệnh nhân được thu thập từ môi trường lâm sàng thực tế (Bệnh viện Chợ Rẫy, TP.HCM). Ảnh CT được trích xuất ma trận Hounsfield Unit (HU), chuyển đổi sang mật độ khối lượng và gán thành 279 hợp chất mô sinh học theo đường cong chuẩn hóa thực nghiệm của Schneider và cộng sự (2000).
  2. Giao thức Hình học Song song (Parallel Geometry):
    • Đường viền phân đoạn cơ quan trong file RTSTRUCT được phần mềm CARIMAS và mô-đun Java trích xuất thành danh sách tọa độ vector 2D trên từng lát cắt.
    • Các lát cắt được liên kết thành hình khối đa diện 3D dạng G4ExtrudedSolid.
    • Trong Geant4/GAMOS 6.0, thế giới hình học thực (Mass World) chứa thông tin vật liệu để xử lý tán xạ và dừng hạt, trong khi thế giới song song (Parallel World) chứa cấu trúc giải phẫu chính xác đóng vai trò bộ chấm điểm (Scorer) thu nhận năng lượng lắng đọng $E_{i, \text{struct}}$.
  3. Giao thức Ghép Phantom và Xoay Vector:
    • Hai tập dữ liệu voxel (bệnh nhân và người tiếp xúc) được hợp nhất thành một tệp văn bản ASCII duy nhất theo cấu trúc dữ liệu của GAMOS.
    • Sử dụng các lệnh biến đổi hình học chuyên dụng để tịnh tiến và quay tọa độ của phantom người tiếp xúc theo các góc Euler: $\theta = 0^\circ$ (đối diện), $\theta = 90^\circ$ (bên cạnh), $\theta = 180^\circ$ (trước - sau) và kịch bản bệnh nhân nằm trên giường bệnh, người chăm sóc đứng cạnh.
  4. Đảm bảo độ tin cậy và hội tụ mô phỏng:
    • Thiết lập số lượng lịch sử mô phỏng (hạt sơ cấp) lên đến $10^8 - 10^9$ biến cố nhằm đảm bảo sai số thống kê Monte Carlo (Fractional Standard Deviation) luôn duy trì dưới 1,0% tại các cơ quan quan tâm.
    • Thuật toán sinh số giả ngẫu nhiên HEP James được kích hoạt với hạt giống (seed) độc lập. Kiểm chứng chéo mã nguồn mô phỏng với các bài toán chuẩn phantom nước hình cầu và so sánh trực tiếp với kết quả đo liều bằng code EGSnrc và MCNP.

Data và phân tích

Phân tích hạt nhân sử dụng đầy đủ phổ phát xạ phân rã của Iodine-131 trích xuất từ cơ sở dữ liệu ENSDF (Evaluated Nuclear Structure Data File) bao gồm:

  • Phổ hạt beta liên tục với năng lượng cực đại $E_{\beta, \max} = 606\text{ keV}$ (xác suất 89,9%), $E_{\beta, \max} = 334\text{ keV}$ (7,27%), $E_{\beta, \max} = 248\text{ keV}$ (2,1%).
  • Các đỉnh gamma đơn năng chính: 364,5 keV (xác suất 81,5%), 637,0 keV (7,17%), 284,3 keV (6,14%), và 722,9 keV (1,8%).

Dữ liệu liều tương đương $H_T$ trên từng cơ quan của người tiếp xúc được tính theo tiêu chuẩn ICRP 103 với trọng số bức xạ $W_R = 1$ cho photon/electron. Liều hiệu dụng toàn thân $E$ được tổng hợp từ 14 nhóm cơ quan với trọng số mô $W_T$ nghiêm ngặt: Tuyến giáp ($W_T = 0,04$), Phổi ($W_T = 0,12$), Dạ dày ($W_T = 0,12$), Tủy xương đỏ ($W_T = 0,12$), Tuyến sinh dục ($W_T = 0,08$), Mô mỡ/da ($W_T = 0,01$), v.v.

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

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

+-----------------------------------------------------------------------------------+
|               BẢNG TỔNG HỢP SO SÁNH CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN              |
+-----------------------------------------------------------------------------------+
| NỘI DUNG KHẢO SÁT       | PHƯƠNG PHÁP TRUYỀN THỐNG      | PHƯƠNG PHÁP CỦA LUẬN ÁN |
+-------------------------+-------------------------------+-------------------------+
| Tính liều Voxel biên    | Voxel hóa nhị phân: Sai lệch  | Hình học song song:     |
| (Nhóm 30 và 300 Voxel)  | năng lượng và liều trung bình | Bảo toàn năng lượng và  |
|                         | lên tới 40 - 60% ở vùng rìa.  | thể tích thực (f_i).    |
| ----------------------- + ----------------------------- + ----------------------- |
| Liều tiếp xúc gần       | Nguồn điểm: Thổi phồng liều   | Ghép Phantom Voxel:     |
| (Khoảng cách 0,3 m)     | lên tới 366% (Albérico, 2008).| Phản ánh đúng tự hấp thụ|
| ----------------------- + ----------------------------- + ----------------------- |
| Tính bất đối xứng tư thế| Không phân biệt được hướng và | Liều tư thế Đối diện    |
| (Khoảng cách d = 1,0 m) | phân bố hình thể bệnh nhân.   | cao gấp 1,8 - 2,4 lần   |
|                         |                               | so với tư thế Sau lưng. |
| ----------------------- + ----------------------------- + ----------------------- |
| Kịch bản Đứng - Nằm     | Bỏ qua che chắn của thành     | Tuyến giáp người đứng   |
| (Chăm sóc tại giường)   | giường và cơ thể bệnh nhân.   | nhận liều trực diện cao |
|                         |                               | do vị trí cổ bệnh nhân. |
+-----------------------------------------------------------------------------------+

Nghiên cứu mang lại 4 phát hiện thực nghiệm mang tính đột phá:

  1. Khắc phục triệt để sai số thể tích biên: Khảo sát trên tập hợp mẫu 30 voxel và 300 voxel tại đường biên giải phẫu tuyến giáp cho thấy phương pháp truyền thống (nhận hoặc loại bỏ toàn bộ voxel) gây ra sự biến động bất thường của liều hấp thụ trung bình. Trong phương pháp hình học song song mới, năng lượng hấp thụ thực tế $E_{\text{struct}}$ tỷ lệ hoàn toàn tương thích với tỷ lệ thể tích $f_i$, loại bỏ các bước nhảy liều phi vật lý tại các lát cắt mỏng.

  2. Chứng minh sự sai lệch nghiêm trọng của mô hình nguồn điểm: Tại khoảng cách tiếp xúc gần 0,3 m đến 0,5 m, mô hình nguồn điểm đánh giá quá cao giá trị liều hiệu dụng từ 200% đến hơn 360% do bỏ qua sự suy giảm bức xạ trong chính các lớp mô dày của bệnh nhân (hiệu ứng tự hấp thụ - self-absorption). Mức sai lệch giảm dần khi khoảng cách tăng lên $d \ge 1,5\text{ m}$, khẳng định định luật nghịch đảo bình phương chỉ có giá trị xấp xỉ ở trường xa ($d > 2\text{ m}$).

  3. Xác lập tính bất đối xứng mạnh của các tư thế tiếp xúc: Liều hiệu dụng và liều tương đương tuyến giáp của người tiếp xúc biến thiên mạnh mẽ theo góc giao tiếp. Ở cùng khoảng cách $d = 1,0\text{ m}$:

    • Tư thế Đối diện cho liều hiệu dụng cao nhất do chùm gamma 364 keV từ tuyến giáp và ổ bụng phát xạ trực diện không bị cản trở bởi cột sống.
    • Tư thế Sau lưng ghi nhận liều thấp nhất, giảm từ 45% đến 60% so với tư thế đối diện nhờ sự che chắn mật độ cao của xương sống và khối cơ lưng bệnh nhân.

  4. Phát hiện động lực học liều trong kịch bản Đứng - Nằm: Trong tình huống người nhà hoặc điều dưỡng đứng chăm sóc cạnh giường bệnh nhân đang nằm, tuyến giáp và vùng ngực của người tiếp xúc nhận liều tương đương cao vượt trội so với liều trung bình toàn thân do góc phát xạ hình học từ cổ bệnh nhân hướng thẳng lên nửa thân trên của người đứng.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung vào hệ thống lý thuyết đo liều y học hạt nhân một phương pháp luận chuẩn xác để tích hợp hình học giải phẫu liên tục vào cấu trúc mô phỏng gián đoạn của Monte Carlo.
  • Về mặt phương pháp: Bộ code phát triển trên nền tảng GAMOS 6.0 có khả năng mở rộng để ứng dụng cho bất kỳ đồng vị phóng xạ nào khác trong điều trị như Lu-177, Y-90, Ra-223, hoặc Ac-225.
  • Về ứng dụng lâm sàng và an toàn bức xạ: Cung cấp căn cứ khoa học chuẩn xác để các bệnh viện thiết lập bảng quy chuẩn cách ly linh hoạt. Thay vì áp đặt cách ly tuyệt đối 100% thời gian, bệnh nhân có thể tương tác với người chăm sóc ở các góc tiếp xúc an toàn (ví dụ: tiếp xúc từ phía sau lưng hoặc giữ khoảng cách tối ưu trên 1,0 m), giúp giải tỏa áp lực tâm lý cho hơn 30% bệnh nhân ung thư điều trị nội trú.

Limitations và Future Research

  1. Những hạn chế cố hữu:

    • Luận án mới tập trung kiểm chứng chi tiết trên mô hình bệnh nhân ung thư tuyến giáp điều trị với I-131; chưa mở rộng trên cỡ mẫu lớn các ca lâm sàng có khối u di căn đa cơ quan.
    • Kỹ thuật tính toán giả định bệnh nhân và người tiếp xúc ở trạng thái hình học tĩnh, chưa tính đến biến dạng cơ học của các cơ quan do nhịp thở hoặc nhu động ruột trong thời gian dài (4D Monte Carlo).
    • Thời gian tính toán Monte Carlo cho không gian ghép hai phantom voxel với số lịch sử lớn ($10^9$) còn kéo dài, đòi hỏi hệ thống máy chủ tính toán hiệu năng cao (HPC), gây trở ngại nhất định cho việc áp dụng lập kế hoạch điều trị tức thời tại giường bệnh.

  2. Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

    • Tích hợp trí tuệ nhân tạo (Deep Learning / Denoising Autoencoders) để giảm số lịch sử mô phỏng cần thiết mà vẫn đảm bảo độ hội tụ, rút ngắn thời gian tính liều từ vài giờ xuống mức vài phút.
    • Mở rộng thuật toán hình học song song cho các hạt nặng phát xạ alpha (Alpha Emitters) như Bi-213, Ac-225 trong liệu pháp hạt nhân nhắm trúng đích vi mô (Targeted Alpha Therapy - TAT).
    • Thiết lập giao diện người dùng đồ họa hoàn chỉnh tích hợp trực tiếp vào hệ thống thông tin bệnh viện (PACS/HIS) phục vụ tính toán liều tiếp xúc tự động cho từng ca bệnh xuất viện.

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

Công trình nghiên cứu mang lại tác động sâu rộng trên cả phương diện học thuật và an toàn xã hội:

  • Giá trị học thuật: Đóng góp vào hệ thống ấn phẩm quốc tế thuộc danh mục ISI/Scopus; nâng cấp trực tiếp tính năng mã nguồn mở của gói phần mềm GAMOS (hợp tác cùng Trung tâm CIEMAT, Tây Ban Nha), tạo tiền đề cho cộng đồng vật lý y khoa toàn cầu ứng dụng.
  • Chuyển đổi y tế và an toàn bức xạ: Cung cấp công cụ tính toán định lượng cho Cục An toàn Bức xạ và Hạt nhân cùng các cơ sở y tế hạt nhân lớn (Bệnh viện Chợ Rẫy, Bệnh viện K, Bệnh viện Bạch Mai) xây dựng hướng dẫn cách ly bệnh nhân dựa trên bằng chứng khoa học vững chắc thay vì ước lượng cảm tính.
  • Lợi ích xã hội: Giảm thiểu số ngày cách ly không cần thiết cho hàng ngàn bệnh nhân tiếp nhận điều trị phóng xạ hàng năm, giảm chi phí lưu trú bệnh viện và bảo vệ tối đa sức khỏe cho đội ngũ y bác sĩ, người chăm sóc và thân nhân.

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

  • Nghiên cứu sinh và Giới nghiên cứu Vật lý Y khoa: Tiếp cận phương pháp luận tiên tiến kết hợp thế giới hình học song song trong Geant4/GAMOS và quy trình xử lý dữ liệu DICOM RTSTRUCT chuẩn mực.
  • Bác sĩ Y học hạt nhân và Kỹ sư Vật lý Bệnh viện: Có trong tay công cụ tính toán liều chiếu trong chính xác cho mô lành của bệnh nhân và định lượng liều phát xạ ngoại vi phục vụ công tác an toàn bức xạ lâm sàng.
  • Cơ quan Quản lý và Hoạch định Chính sách An toàn Bức xạ: Có cơ sở số liệu thực chứng để cập nhật các tiêu chuẩn an toàn bức xạ quốc gia phù hợp với khuyến cáo của IAEA và ICRP.
  • Bệnh nhân và Thân nhân: Hưởng lợi trực tiếp từ phác đồ điều trị tối ưu hóa liều lượng và chế độ sinh hoạt, chăm sóc nhân văn, giảm thiểu gánh nặng tâm lý xã hội.

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ì và mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc tái thiết lập mô hình tính toán liều hấp thụ cho cấu trúc giải phẫu trên nền lưới voxel thông qua kỹ thuật Hình học song song (Parallel Geometry), trực tiếp mở rộng Sơ đồ Đo liều Chiếu trong MIRD (Medical Internal Radiation Dose schema) và Lý thuyết Vận chuyển Hạt Lượng tử. Luận án đã giải quyết căn bản nghịch lý "voxel hóa nhị phân" vốn tồn tại nhiều thập kỷ, cho phép tính toán chính xác phần năng lượng phân đoạn $E_{i, \text{struct}}$ lắng đọng trong thể tích thực $f_i \cdot V_{\text{voxel}}$ của từng voxel biên bị cắt bởi đường phân đoạn RTSTRUCT DICOM.

2. Đột phá phương pháp luận của luận án có gì vượt trội so với các nghiên cứu tiền nhiệm?

So với nghiên cứu của Albérico và cộng sự (2008) (vốn áp đặt nguồn phân bố đồng nhất trên phantom FAX) hay Palmer và cộng sự (2014) (dùng mô hình hình trụ đơn giản hóa), phương pháp của luận án vượt trội toàn diện ở hai khía cạnh:

  • Tích hợp trực tiếp dữ liệu hoạt độ phóng xạ dị dưỡng từ ảnh SPECT/PET lâm sàng của bệnh nhân vào mô hình mô phỏng.
  • Lần đầu tiên phát triển thuật toán hợp nhất hai phantom voxel dị dưỡng độc lập trong hệ quy chiếu xoay vector 3D linh hoạt trên GAMOS 6.0/Geant4, mô tả chính xác mọi góc độ và khoảng cách tương tác trong không gian thực nghiệm ảo mà không làm suy giảm độ phân giải vật liệu.
       SO SÁNH CÁC PHƯƠNG PHÁP TÍNH LIỀU TIẾP XÚC TRONG Y HỌC HẠT NHÂN
+-------------------+--------------------+--------------------+--------------------+
| TIÊU CHÍ SO SÁNH  | NGUỒN ĐIỂM / TRỤ   | ALBÉRICO ET AL.    | LUẬN ÁN NÀY        |
|                   | (Siegel / Palmer)  | (2008 - FAX)       | (GAMOS / Geant4)   |
+-------------------+--------------------+--------------------+--------------------+
| Phân bố hoạt độ   | Nguồn điểm /       | Voxel hóa đồng     | Ma trận SPECT/PET  |
| nguồn bệnh nhân   | Hình trụ đồng nhất | nhất toàn tuyến giáp| dị dưỡng thực tế   |
| ----------------- + ------------------ + ------------------ + ------------------ |
| Mô hình hình học  | Điểm toán học /    | Phantom Voxel FAX  | Ghép 2 Phantom     |
| người tiếp xúc    | Đoạn thẳng         | (Tư thế cố định)   | Voxel (ICRP 110)   |
| ----------------- + ------------------ + ------------------ + ------------------ |
| Tự do chuyển động | Bất biến không gian| Rất hạn chế        | Xoay Vector 3D     |
| và xoay góc 3D    | (Đối xứng cầu)     |                    | (Mọi tư thế thực)  |
| ----------------- + ------------------ + ------------------ + ------------------ |
| Độ chính xác liều | Sai lệch lên đến   | Sai số phân bố do  | Khớp chuẩn thực    |
| trường gần (0,3m) | 366% (Thổi phồng)  | nguồn đồng nhất    | nghiệm Monte Carlo |
+-------------------+--------------------+--------------------+--------------------+

3. Phát hiện nào bất ngờ nhất về mặt dữ liệu thực nghiệm?

Phát hiện bất ngờ nhất là mức độ chênh lệch liều lên tới 366% của mô hình nguồn điểm ở khoảng cách gần (0,3 m) so với mô hình phantom voxel thực, cùng với sự sụt giảm liều hơn 50% ở tư thế tiếp xúc từ Sau lưng so với Đối diện tại cùng khoảng cách $d = 1,0\text{ m}$. Kết quả này bác bỏ hoàn toàn quan niệm lâm sàng phổ biến cho rằng sự suy giảm liều phát xạ từ bệnh nhân ở khoảng cách gần tuân theo đối xứng cầu đơn giản, chứng minh vai trò tự che chắn cực kỳ quan trọng của các cấu trúc giải phẫu nội tại (cột sống, khung xương sườn, mô mỡ).

4. Luận án có cung cấp giao thức tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp giao thức tái lập khoa học hoàn chỉnh và minh bạch bao gồm:

  • Toàn bộ thuật toán chuyển đổi cấu trúc RTSTRUCT sang tệp lệnh GAMOS g4dcm/ASCII.
  • Bảng thông số vật liệu và quy trình chuyển đổi số HU sang mật độ theo Schneider (2000).
  • Hệ thống lệnh điều khiển phân bổ nguồn, thiết lập danh mục vật lý GmPhysicsList (mô hình Livermore), bộ chấm điểm liều GmDoseScorer và các bộ lọc GmFilter trong GAMOS 6.0.
  • Quy trình chuẩn hóa đầu dò, kiểm tra chất lượng hình ảnh trên CARIMAS và mã nguồn định hướng xoay vector không gian.

5. Định hướng chương trình nghiên cứu 10 năm tới được vạch ra như thế nào?

Chương trình nghiên cứu 10 năm tập trung vào ba hướng chiến lược:

  1. Nâng cấp mô hình đo liều sang 4D Monte Carlo tích hợp biến dạng hô hấp thời gian thực.
  2. Ứng dụng mô phỏng đo liều vi mô cho liệu pháp kháng thể gắn xạ phóng xạ hạt alpha năng lượng cao (TAT) với Actinium-225 và Bismuth-213.
  3. Ứng dụng điện toán đám mây và mạng nơ-ron sâu để tự động hóa hoàn toàn quy trình tính liều cá thể hóa trong vòng dưới 5 phút cho mỗi ca điều trị lâm sàng.

Kết luận

Công trình nghiên cứu luận án tiến sĩ của tác giả Nguyễn Thị Phương Thảo đã hoàn thành xuất sắc các mục tiêu khoa học đặt ra, đóng góp 5 giá trị cốt lõi cho ngành Vật lý nguyên tử và hạt nhân:

  1. Phát triển hoàn chỉnh phương pháp "Hình học song song" trên GAMOS/Geant4: Xử lý triệt để bài toán sai số thể tích biên tại các vùng ranh giới cơ quan phức tạp, bảo toàn nguyên vẹn hình học giải phẫu RTSTRUCT DICOM trên lưới phantom voxel.
  2. Thiết lập thuật toán "Ghép Phantom và Xoay Vector 3D": Lần đầu tiên cho phép mô phỏng đồng thời hai phantom voxel bệnh nhân và người tiếp xúc với dữ liệu phân bố hoạt độ SPECT/PET thực tế, thiết lập trường bức xạ chân thực trong không gian tương tác đa chiều.
  3. Định lượng hóa toàn diện tính bất đối xứng của trường bức xạ tiếp xúc: Cung cấp bộ dữ liệu liều tương đương và liều hiệu dụng chuẩn xác trên các kịch bản giao tiếp lâm sàng trọng yếu (đối diện, bên cạnh, trước - sau, đứng - nằm), làm sáng tỏ sự sai lệch lên đến 366% của các mô hình nguồn điểm truyền thống.
  4. Cung cấp bằng chứng thực nghiệm tối ưu hóa an toàn bức xạ: Đặt nền tảng khoa học vững chắc để cá thể hóa quy trình cách ly y tế, giảm thiểu rủi ro bức xạ cho nhân viên y tế và cộng đồng, đồng thời nâng cao chất lượng chăm sóc tâm lý cho bệnh nhân điều trị phóng xạ.
  5. Mở rộng và chuẩn hóa công cụ mô phỏng mở: Tích hợp trực tiếp các mô-đun cải tiến vào phần mềm quốc tế GAMOS/Geant4, khẳng định năng lực hội nhập và đóng góp học thuật đỉnh cao của nghiên cứu hạt nhân Việt Nam trên trường quốc tế.