Tổng quan nghiên cứu

Ung thư vú hiện là bệnh lý ác tính phổ biến hàng đầu ở nữ giới trên phạm vi toàn cầu. Theo các thống kê dịch tễ quốc tế, năm 2012 thế giới đã ghi nhận khoảng 1,67 triệu ca mắc mới (chiếm 25% tổng số các ca ung thư) và cướp đi sinh mạng của hơn 522.000 người. Tại khu vực Châu Á - Thái Bình Dương, tỷ lệ mắc bệnh có xu hướng tăng nhanh, đạt 65,7 ca trên 100.000 dân tại Singapore và lên tới 118 ca trên 100.000 dân tại Úc. Trước thực trạng đó, việc tầm soát sớm đóng vai trò then chốt nhằm nâng cao tỷ lệ sống sót sau 5 năm cho người bệnh.

Phương pháp chụp X-quang nhũ ảnh truyền thống (Mammography) từ lâu được xem là tiêu chuẩn vàng trong chẩn đoán sớm nhờ liều bức xạ thấp và độ nhạy cao với các vi vôi hóa. Tuy nhiên, kỹ thuật chụp chiếu 2D này bộc lộ hạn chế lớn đối với phụ nữ có mật độ mô vú dày đặc, khi hiện tượng chồng lấn giữa mô tuyến và mô mỡ có thể che khuất các khối u ác tính. Để khắc phục triệt để nhược điểm trên, kỹ thuật chụp cắt lớp tuyến vú kỹ thuật số (Breast Tomosynthesis) đã ra đời, cho phép thu nhận hình ảnh ở nhiều góc độ khác nhau để tái tạo không gian 3D.

Nghiên cứu này tập trung ứng dụng chương trình mô phỏng Monte Carlo PENELOPE kết hợp công cụ tạo ảnh PENEASY nhằm mô phỏng chi tiết quá trình truyền xạ trong kỹ thuật Mammography và Breast Tomosynthesis. Mục tiêu trọng tâm là khảo sát sự truyền photon, đánh giá liều tuyến vú trung bình và xác định cấu hình góc chụp tối ưu để phát hiện các khối vi vôi hóa có kích thước nhỏ dưới 4 mm. Đề tài thuộc chuyên ngành Vật lý Kỹ thuật, được hoàn thành tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh trong giai đoạn từ tháng 07/2016 đến tháng 01/2017. Kết quả nghiên cứu mang lại ý nghĩa thực tiễn to lớn trong việc cung cấp các thông số hình học và cơ sở vật lý chuẩn xác, phục vụ công tác thiết kế, tối ưu hóa thiết bị X-quang nhũ ảnh thế hệ mới với liều hấp thụ tối thiểu 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 nền tảng lý thuyết tương tác bức xạ với vật chất và các mô hình toán học trong vật lý y khoa:

  • Lý thuyết truyền xạ Monte Carlo qua mã PENELOPE: Chương trình PENELOPE mô phỏng chi tiết quỹ đạo và sự mất mát năng lượng của các hạt photon, electron và positron khi tương tác với môi trường vật chất bất kỳ trong dải năng lượng từ 100 eV đến 1 GeV. Hệ thống dữ liệu vật lý bao gồm các tiết diện tán xạ đàn hồi, tán xạ không đàn hồi, bức xạ hãm (Bremsstrahlung), hiệu ứng quang điện và tán xạ Compton.
  • Lý thuyết tạo ảnh X-quang và tái tạo cắt lớp: Nguyên lý suy giảm cường độ chùm tia X theo định luật truyền thẳng qua các mô sinh học có hệ số suy giảm khác nhau. Quá trình tái tạo ảnh 3D trong kỹ thuật Breast Tomosynthesis áp dụng thuật toán chiếu ngược lọc (Filtered Backprojection - FBP) kết hợp bộ lọc dốc (ramp filter), bộ lọc quang phổ Hanning và hàm lọc kiểm soát độ dày lát cắt (Hprofile) để loại bỏ ảnh giả out-of-plane.
  • Mô hình đánh giá liều lượng bức xạ: Đánh giá thông qua đại lượng liều tuyến vú trung bình (Mean Glandular Dose - Dg) và liều tương đối (Relative Glandular Dose - RGD), cho phép phân tích sự phân bố năng lượng bức xạ lắng đọng trong mô tuyến tại các góc chụp khác nhau.
  • Các khái niệm then chốt: Mô tuyến vú (glandular tissue), mô mỡ (adipose tissue), vi vôi hóa (microcalcifications - CaCO3), hàm chuyển đổi điều chế (Modulation Transfer Function - MTF) và hiệu suất lượng tử đầu dò (Detective Quantum Efficiency - DQE).

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

Nghiên cứu sử dụng phương pháp mô phỏng số thông qua chuỗi công cụ tích hợp chuyên dụng:

  • Nguồn dữ liệu và phổ tia X: Dữ liệu phổ phát xạ tia X được tạo lập từ phần mềm Spekcalc mô phỏng bóng phát tia anode Tungsten (W) với điện áp gia tốc 30 kVp, sử dụng cửa sổ lọc Berili dày 0,5 mm và bộ lọc Nhôm dày 0,5 mm. Cấu hình này mô phỏng chuẩn xác hệ thống thiết bị thực nghiệm CapeRay Soteria mammography với khoảng cách từ nguồn đến đầu dò là 65 cm.
  • Mô hình hóa mẫu vật (Phantom): Cỡ mẫu mô phỏng sử dụng phantom chuẩn CIRS Model 011A có kích thước hình học 12,5 x 10 x 4,5 cm. Phantom bao gồm 5 khối mô mềm kích thước 1 x 1 cm đại diện cho các tỷ lệ mật độ khác nhau (100% tuyến, 70% tuyến - 30% mỡ, 50% tuyến - 50% mỡ, 30% tuyến - 70% mỡ, 100% mỡ) và các khối cầu vi vôi hóa CaCO3 kích thước xấp xỉ 4 mm. Toàn bộ hình học được định nghĩa chính xác thông qua công cụ PENGEOM và PENVOX.
  • Hệ thống đầu dò và xử lý ảnh: Đầu dò phẳng CsI kích thước 20 x 20 cm, độ dày 1 mm, ma trận điểm ảnh 500 x 500 pixels với kích thước pixel đạt 0,04 cm (400 µm). Thuật toán ghép và tái tạo ảnh đa góc được lập trình thực hiện trên môi trường MATLAB.
  • Lý do lựa chọn phương pháp: Mô phỏng Monte Carlo là giải pháp tối ưu cho phép khảo sát sâu các tương tác lượng tử ở cấp độ vi mô, kiểm soát hoàn toàn các thông số vật lý mà không gây nguy cơ nhiễm xạ thực nghiệm trên cơ thể người. Tổng thời gian chạy mô phỏng cho mỗi cấu hình kéo dài liên tục 4 giờ để đảm bảo độ tin cậy thống kê cao. Toàn bộ tiến trình nghiên cứu được thực hiện trong thời gian 6 tháng (07/2016 – 12/2016).

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

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

Quá trình mô phỏng và phân tích định lượng đã mang lại 4 phát hiện khoa học quan trọng:

  1. Độ tương phản rõ rệt giữa các tỷ lệ mô vú: Trên hình ảnh mô phỏng PENEASY ở mức năng lượng 30 keV, hệ thống phân biệt sắc nét 5 cấp độ mật độ mô mềm. Khối mô tuyến 100% (mật độ 1 g/cm³) và khối mô mỡ 100% (mật độ 1 g/cm³) thể hiện mức độ suy giảm tia X chênh lệch rõ rệt, cho phép nhận diện chính xác ranh giới các vùng mô có tỷ lệ tuyến 70%, 50% và 30%.
  2. Khả năng phát hiện các tổn thương vi vôi hóa: Đầu dò CsI với độ phân giải ma trận 500 x 500 pixels ghi nhận rõ ràng 4 khối vi vôi hóa CaCO3 kích thước 4 mm. Tín hiệu độ xám của các hạt vi vôi hóa nổi bật hoàn toàn so với nền mô tuyến và mô mỡ xung quanh, chứng minh tính khả thi cao trong việc phát hiện các dấu hiệu ung thư giai đoạn sớm.
  3. Tối ưu hóa cấu hình góc quét trong Breast Tomosynthesis: Trong 4 cấu hình góc quét được khảo sát, phương án quét từ 70° đến 110° với 21 góc chiếu (bước quét 2°) mang lại chất lượng hình ảnh tái tạo lát cắt vượt trội nhất. Hình ảnh chồng lát có độ sắc nét biên cao hơn hẳn so với các cấu hình 11 góc (bước quét 4°) và cấu hình góc hẹp 80° đến 100° (6 góc hoặc 11 góc).
  4. Tương quan giữa số góc quét và chất lượng dữ liệu: Việc tăng số lượng góc chiếu giúp làm giảm hiện tượng ảnh giả sọc, nhưng nếu mở rộng góc quét quá lớn mà không tăng mật độ góc chụp sẽ làm suy giảm độ phân giải của lát cắt. Mức liều tuyến vú tương đối RGD duy trì độ biến thiên ổn định trong khoảng 0,5 đến 1,0 mGy đối với các góc chụp đối xứng qua trục 0°.

Thảo luận kết quả

Nguyên nhân giúp cấu hình 21 góc chiếu (dải góc 70° - 110°) đạt chất lượng tối ưu là do việc lấy mẫu không gian tần số Fourier trong vùng double-wedge domain được lấp đầy đồng đều hơn, giúp thuật toán Filtered Backprojection triệt tiêu hiệu quả các thành phần tần số gây nhiễu và hạn chế tối đa ảnh giả lan truyền giữa các lát cắt lân cận (out-of-plane artifacts). Khi bước góc quét giảm từ 4° xuống 2°, thông tin cấu trúc vi mô của khối vôi hóa được bảo toàn nguyên vẹn hơn khi thực hiện phép tích phân ngược.

So sánh với các nghiên cứu thực nghiệm quốc tế trên hệ thống CapeRay Soteria và nghiên cứu mô phỏng của Sechopoulos (2009), kết luận về sự đánh đổi giữa độ rộng góc và số góc chiếu hoàn toàn tương đồng. Trong thực tế, liều tuyến vú trung bình của kỹ thuật Breast Tomosynthesis (khoảng 1,30 mGy ở độ dày vú 5 cm) chỉ cao hơn khoảng 8% so với phương pháp Mammography truyền thống (khoảng 1,20 mGy), nhưng lại cung cấp khả năng bóc tách tổn thương vượt trội.

Các dữ liệu nghiên cứu này có thể được biểu diễn trực quan qua biểu đồ phổ năng lượng phát xạ Spekcalc từ 0 đến 30 keV và bảng ma trận đối chiếu chất lượng hình ảnh giữa các bước góc 2° và 4°. Những số liệu này khẳng định rằng việc lựa chọn bước quét nhỏ với dải góc trung bình là chìa khóa để cân bằng giữa độ phân giải hình ảnh và mức độ phơi nhiễm bức xạ của người bệnh.

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

Nhằm chuyển hóa các kết quả mô phỏng vào thực tiễn chế tạo và ứng dụng lâm sàng, 4 giải pháp trọng tâm được đề xuất như sau:

  • Nâng cấp độ phân giải vật lý của hệ thống đầu dò: Các đơn vị nghiên cứu và sản xuất thiết bị y tế cần cải tiến kích thước pixel từ mức 0,04 cm (400 µm) xuống dải tiêu chuẩn 70 - 85 µm sử dụng công nghệ chuyển đổi trực tiếp a-Se hoặc tấm phẳng CsI màng mỏng. Mục tiêu là nâng tỷ lệ nhận diện các cụm vi vôi hóa kích thước dưới 1 mm lên trên 95% trong lộ trình thực hiện 12 tháng.
  • Chuẩn hóa giao thức góc quét trên thiết bị nhũ ảnh cắt lớp: Các kỹ sư phát triển hệ thống cần thiết lập cố định chế độ quét từ 13 đến 21 góc chiếu trong phạm vi góc mở từ 30° đến 60° (tương đương biên độ ±15° đến ±30°) cho các dòng máy thương mại. Khuyến nghị các bệnh viện chuyên khoa triển khai thử nghiệm lâm sàng giao thức này trong vòng 6 tháng nhằm đánh giá hiệu quả chẩn đoán trên các phân nhóm mô vú dày.
  • Tối ưu hóa tham số điện áp phát tia và vật liệu bộ lọc: Bộ phận kỹ thuật hình ảnh tại các cơ sở y tế nên áp dụng phổ tia X từ bia Tungsten kết hợp bộ lọc Rhodium hoặc Bạc ở dải điện áp 28 - 34 kVp. Giải pháp này giúp giảm thiểu 20% đến 25% liều tuyến vú trung bình cho bệnh nhân mà vẫn duy trì tỷ số tín hiệu trên nhiễu (SNR) ở mức tối ưu.
  • Phát triển thuật toán tái tạo ảnh lặp nâng cao: Các nhóm nghiên cứu thuật toán y sinh cần đẩy mạnh ứng dụng kỹ thuật tái tạo lặp đồng thời (Simultaneous Iterative Reconstruction Technique) thay thế dần thuật toán chiếu ngược lọc cổ điển. Kế hoạch này cần hoàn thành trong vòng 18 tháng nhằm triệt tiêu hoàn toàn ảnh giả out-of-plane và tối ưu hóa thời gian xử lý dữ liệu lát cắt.

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

Công trình nghiên cứu mang lại giá trị học thuật và ứng dụng chuyên sâu cho 4 nhóm đối tượng chính:

  • Học viên, nghiên cứu sinh và giảng viên ngành Vật lý Y khoa / Vật lý Kỹ thuật: Luận văn là tài liệu tham khảo chuẩn mực về phương pháp lập trình mô phỏng Monte Carlo bằng PENELOPE/PENEASY, cung cấp hướng dẫn chi tiết cách mô hình hóa phantom 3D và tính toán liều lượng bức xạ hạt nhân trong dải năng lượng từ 100 eV đến 1 GeV.
  • Kỹ sư thiết kế và phát triển thiết bị y tế X-quang: Cung cấp bộ thông số kỹ thuật thực chứng về cấu hình bóng phát tia anode Tungsten, độ dày bộ lọc Berili/Nhôm và ma trận đầu dò CsI 500 x 500 pixels, hỗ trợ trực tiếp cho quá trình R&D các dòng máy nhũ ảnh kỹ thuật số thế hệ mới.
  • Bác sĩ chẩn đoán hình ảnh và kỹ thuật viên X-quang tuyến vú: Giúp đội ngũ y tế nắm vững bản chất vật lý hình thành ảnh, cơ chế xuất hiện ảnh giả và sự khác biệt về mặt quang học giữa ảnh chiếu 2D và lát cắt 3D, từ đó nâng cao độ chính xác khi đọc kết quả nhũ ảnh lâm sàng.
  • Chuyên gia an toàn bức xạ và các nhà quản lý y tế: Cung cấp cơ sở dữ liệu định lượng về liều tuyến vú trung bình (Mean Glandular Dose) để xây dựng các quy chuẩn kiểm định chất lượng máy chụp X-quang và thiết lập giới hạn liều an toàn trong các chương trình tầm soát ung thư vú cộng đồng.

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

Kỹ thuật Breast Tomosynthesis vượt trội hơn Mammography truyền thống ở những điểm nào?
Breast Tomosynthesis loại bỏ triệt để hiện tượng chồng lấn mô nhờ cơ chế quét đa góc và tái tạo không gian 3D dạng lát cắt. Kỹ thuật này giúp tăng tỷ lệ phát hiện khối u ác tính từ 25% đến 30% đối với phụ nữ có mô vú dày đặc so với hình ảnh chiếu 2D thông thường.

Chương trình PENELOPE đóng vai trò then chốt gì trong nghiên cứu mô phỏng nhũ ảnh?
PENELOPE sử dụng thuật toán Monte Carlo để mô phỏng chính xác đường đi và sự mất mát năng lượng của photon, electron trong dải từ 100 eV đến 1 GeV. Công cụ này cho phép đánh giá liều bức xạ và kiểm tra khả năng tạo ảnh trên máy tính mà không cần chiếu xạ thực nghiệm trên người.

Cấu hình góc quét nào mang lại chất lượng lát cắt nhũ ảnh tốt nhất?
Kết quả mô phỏng cho thấy cấu hình quét trong dải góc 70° đến 110° với 21 góc chiếu (bước quét 2°) mang lại độ sắc nét cao nhất. Cấu hình này giúp làm đầy dữ liệu không gian tần số Fourier và giảm thiểu tối đa hiện tượng mờ nhòe do ảnh giả out-of-plane.

Liều bức xạ của phương pháp Breast Tomosynthesis có nằm trong giới hạn an toàn không?
Có, liều tuyến vú trung bình của Breast Tomosynthesis trên mô vú dày 5 cm đạt khoảng 1,30 mGy, chỉ chênh lệch khoảng 8% so với mức 1,20 mGy của Mammography. Mức liều này hoàn toàn nằm trong giới hạn an toàn nghiêm ngặt của các tiêu chuẩn y tế quốc tế.

Làm thế nào để mô phỏng chính xác cấu trúc mô vú người trong nghiên cứu?
Nghiên cứu sử dụng phantom chuẩn CIRS Model 011A kích thước 12,5 x 10 x 4,5 cm, định nghĩa qua module PENGEOM với dữ liệu vật liệu ICRU. Mô hình phân chia 5 khối mô mềm theo tỷ lệ tuyến/mỡ khác nhau kết hợp các khối cầu vi vôi hóa CaCO3 kích thước 4 mm.

Kết luận

  • Làm chủ hoàn toàn kỹ thuật mô phỏng Monte Carlo thông qua bộ công cụ PENELOPE, PENEASY và Spekcalc trong nghiên cứu vật lý y khoa.
  • Mô hình hóa thành công cấu trúc phantom vú CIRS 011A với 5 cấp độ mật độ mô và nhận diện chính xác các cụm tổn thương vi vôi hóa 4 mm.
  • Xác lập cấu hình góc quét tối ưu gồm 21 góc chiếu trong biên độ 70° - 110° (bước quét 2°), tạo tiền đề cho thuật toán tái tạo lát cắt đạt độ tương phản cao.
  • Chứng minh tính ưu việt của hệ thống bóng phát tia anode Tungsten ở mức điện áp 30 kVp giúp cân bằng giữa độ nét hình ảnh và việc kiểm soát liều tuyến vú.
  • Đóng góp nguồn dữ liệu tham chiếu vật lý chuẩn xác, phục vụ trực tiếp cho định hướng tự chủ thiết kế và chế tạo máy nhũ ảnh cắt lớp 3D tại Việt Nam.

Đóng góp lớn nhất của luận văn là đã chuẩn hóa toàn bộ quy trình mô phỏng quang hình học từ nguồn phát, tương tác vật liệu đến đầu dò kỹ thuật số. Trong giai đoạn 2024 - 2026, các hướng nghiên cứu tiếp theo cần tập trung nâng cấp độ phân giải pixel đầu dò xuống dưới 100 µm và tích hợp các thuật toán tái tạo lặp thông minh. Các viện nghiên cứu, trường đại học và doanh nghiệp y tế nên khai thác triệt để kết quả của công trình này nhằm đẩy nhanh tiến độ nội địa hóa các dòng thiết bị chẩn đoán hình ảnh công nghệ cao.