Luận Văn Thạc Sĩ Vật Lý Kỹ Thuật: Mô Phỏng Hình Ảnh X-Quang Sử Dụng Phương Pháp Monte Carlo

Luận văn thạc sĩ vật lý kỹ thuật nghiên cứu mô phỏng hình ảnh X-quang bằng phương pháp Monte Carlo, ứng dụng hiệu quả trong y học và khoa học.

• Tác giả: Bùi Thiện Chánh
• Người hướng dẫn: PGS. Huỳnh Quang Linh
• Trường học: Đại học Bách Khoa, Đại học Quốc gia TP. Hồ Chí Minh
• Chuyên ngành: Vật Lý Kỹ Thuật
• Đề tài: Mô Phỏng Hình Ảnh X-Quang Bằng Phương Pháp Monte Carlo
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2015
• Địa điểm: Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• DANH MỤC CÁC HÌNH – ĐỒ THỊ
• MỞ ĐẦU
• 1. CHƯƠNG 1: TƯƠNG TÁC CỦA BỨC XẠ GAMMA VỚI VẬT CHẤT
• 1.1. Bản chất của tia gamma
• 1.2. Tính chất của tia gamma
• 1.3. Tương tác của tia gamma với vật chất
• 1.3.1. Hiệu ứng quang điện
• 1.3.2. Hiệu ứng Compton
• 1.3.3. Hiệu ứng tạo cặp
• 1.4. Sự suy giảm của chùm bức xạ gamma khi đi qua vật chất
• 1.4.1. Sự suy giảm của chùm gamma hẹp
• 1.4.2. Sự suy giảm của chùm gamma rộng. Hệ số tích lũy
• 1.4.3. Hệ số hấp thụ tuyến tính
• 1.4.4. Hệ số suy giảm khối lượng
• 2. CHƯƠNG 2: CẤU TẠO VÀ NGUYÊN LÝ MÁY PHÁT TIA X
• 2.1. Cấu tạo máy phát X – quang thông thường
• 2.1.1. Cấu tạo ống phát tia X
• 2.1.2. Vỏ ống chân không, dầu tản nhiệt, vỏ kim loại, cửa sổ ống phát tia X
• 2.1.2.1. Vỏ ống chân không
• 2.1.2.2. Dầu tản nhiệt (dầu biến áp)
• 2.1.2.3. Vỏ kim loại
• 2.1.2.4. Cửa sổ ống
• 2.1.3. Bộ lọc tia
• 2.1.4. Hệ chuẩn trực đầu đèn (collimator)
• 2.2. Tương tác tạo ra tia X
• 2.2.1. Bức xạ đặc trưng
• 2.3. Nguyên lý hoạt động của máy phát X- quang
• 2.4. Các yếu tố ảnh hưởng đến chất lượng liều ra tia X
• 2.5. Nguyên lý tạo ảnh X- quang
• 2.5.1. Hệ thống ghi ảnh
• 2.5.2. Đặc trưng kỹ thuật
• 2.5.3. Bức xạ tán xạ
• 3. CHƯƠNG 3: THỰC HÀNH KHẢO SÁT HÌNH ẢNH X – QUANG BẰNG PHƯƠNG PHÁP MÔ PHỎNG MONTE CARLO
• 3.1. Chương trình mô phỏng PENELOPE
• 3.2. Giới thiệu về phương pháp mô phỏng Monte Carlo
• 3.3. Tính toán liều với Monte Carlo
• 3.4. Các yếu tố cơ bản trong tính toán liều Monte Carlo
• 3.4.1. Mô hình vật lý
• 3.4.2. Bảng dữ liệu tương tác
• 3.4.3. Bộ phát số ngẫu nhiên
• 3.4.4. Cấu trúc hình học
• 3.4.5. Cấu tạo vật liệu
• 3.4.6. Đặc điểm của nguồn
• 3.4.7. Làm giảm sự thăng giáng và làm xấp xỉ
• 3.5. Tổng quan về chương trình PENELOPE
• 3.5.1. Hệ thống code PENELOPE
• 3.5.2. Cơ sở dữ liệu và file đầu vào
• 3.5.3. Giới thiệu các chương trình chính (MAIN)
• 3.5.3.1. Chương trình PENSLAB
• 3.5.3.2. Chương trình PENCYL
• 3.5.3.3. Chương trình PENMAIN
• KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Mô phỏng hình ảnh X quang

Mô phỏng hình ảnh X-quang là một phương pháp quan trọng trong nghiên cứu và ứng dụng y tế. Phương pháp này giúp dự đoán và phân tích hình ảnh X-quang thông qua các mô hình toán học và vật lý. Trong luận văn, phương pháp Monte Carlo được sử dụng để mô phỏng quá trình lan truyền của photon tia X trong vật chất. Điều này cho phép hiểu rõ hơn về các tương tác vật lý như hiệu ứng quang điện, tán xạ Compton và hiệu ứng tạo cặp. Kỹ thuật mô phỏng này không chỉ giúp cải thiện chất lượng hình ảnh mà còn giảm thiểu rủi ro liên quan đến bức xạ.

1.1. Nguyên lý mô phỏng

Nguyên lý của mô phỏng hình ảnh X-quang dựa trên sự lan truyền và hấp thụ của tia X trong vật chất. Phương pháp Monte Carlo mô phỏng các tương tác của photon với vật liệu, từ đó tạo ra hình ảnh 2D hoặc 3D. Chương trình Penelope được sử dụng để mô phỏng các mô hình khối lập phương, khối trụ và khối cầu. Hình ảnh được tạo ra dựa trên sự phân bố năng lượng trong lớp đầu dò, giúp đánh giá chất lượng hình ảnh và liều bức xạ.

1.2. Ứng dụng trong y tế

Mô phỏng hình ảnh X-quang có nhiều ứng dụng trong y tế, đặc biệt trong chẩn đoán và điều trị. Phương pháp này giúp tối ưu hóa liều bức xạ, giảm thiểu tác hại đến bệnh nhân. Ngoài ra, nó còn được sử dụng trong nghiên cứu các phương pháp chẩn đoán mới như đếm photon và tái tạo vật liệu. Kỹ thuật hình ảnh này cũng hỗ trợ trong việc kiểm tra khuyết tật vật liệu và mối hàn trong công nghiệp.

II. Phương pháp Monte Carlo

Phương pháp Monte Carlo là một công cụ mạnh mẽ trong mô phỏng các quá trình vật lý phức tạp. Trong luận văn, phương pháp này được áp dụng để mô phỏng sự lan truyền của photon tia X trong vật chất. Kỹ thuật mô phỏng này dựa trên các thuật toán ngẫu nhiên để tính toán các tương tác vật lý như hấp thụ, tán xạ và quang điện. Phương pháp này cho phép dự đoán chính xác hình ảnh X-quang và đánh giá chất lượng hình ảnh dựa trên các thông số vật lý.

2.1. Cơ sở lý thuyết

Phương pháp Monte Carlo dựa trên các nguyên lý thống kê và xác suất để mô phỏng các quá trình vật lý. Trong mô phỏng tia X, phương pháp này tính toán các tương tác của photon với vật chất thông qua các hiệu ứng như quang điện, tán xạ Compton và tạo cặp. Chương trình Penelope được sử dụng để thực hiện các mô phỏng này, với các thông số đầu vào như cấu trúc hình học, vật liệu và nguồn phát xạ.

2.2. Ưu điểm và hạn chế

Phương pháp Monte Carlo có ưu điểm là độ chính xác cao và khả năng mô phỏng các quá trình phức tạp. Tuy nhiên, phương pháp này đòi hỏi thời gian tính toán dài và tài nguyên máy tính lớn. Để giảm thiểu thời gian, các kỹ thuật như làm giảm sự thăng giáng và xấp xỉ được áp dụng. Kỹ thuật mô phỏng này cũng cần được cải tiến để tăng hiệu quả tính toán trong các ứng dụng thực tế.

III. Luận văn thạc sĩ Vật lý kỹ thuật

Luận văn thạc sĩ với chủ đề Mô phỏng hình ảnh X-quang bằng phương pháp Monte Carlo đã đóng góp quan trọng vào lĩnh vực Vật lý kỹ thuật. Nghiên cứu này không chỉ làm sáng tỏ các quá trình tương tác của tia X với vật chất mà còn đề xuất các phương pháp mô phỏng hiệu quả. Kỹ thuật hình ảnh được phát triển trong luận văn có tiềm năng ứng dụng rộng rãi trong y tế và công nghiệp, đặc biệt trong việc cải thiện chất lượng hình ảnh và giảm thiểu rủi ro bức xạ.

3.1. Mục tiêu nghiên cứu

Mục tiêu chính của luận văn là thực hiện mô phỏng hình ảnh X-quang bằng phương pháp Monte Carlo thông qua chương trình Penelope. Nghiên cứu tập trung vào việc mô phỏng các vật thể đơn giản như khối lập phương, khối trụ và khối cầu. Kết quả mô phỏng được so sánh với các dữ liệu thực nghiệm để đánh giá độ chính xác và hiệu quả của phương pháp.

3.2. Đóng góp khoa học

Luận văn đã đóng góp quan trọng vào việc phát triển các kỹ thuật mô phỏng trong Vật lý kỹ thuật. Nghiên cứu này cung cấp cơ sở lý thuyết và thực tiễn cho việc ứng dụng phương pháp Monte Carlo trong mô phỏng hình ảnh X-quang. Đồng thời, luận văn cũng đề xuất các hướng nghiên cứu tiếp theo để cải thiện hiệu quả và độ chính xác của các phương pháp mô phỏng.

Nội dung chính

Tổng quan nghiên cứu

Kỹ thuật chẩn đoán hình ảnh X-quang đóng vai trò then chốt trong y tế hiện đại, nhưng hiệu suất chuyển hóa năng lượng của bóng phát tia X thực tế rất thấp khi hơn 99% động năng electron biến thành nhiệt năng và chưa đến 1% chuyển thành bức xạ hữu ích. Bên cạnh giá trị chẩn đoán, bức xạ ion hóa mang nguy cơ tổn thương sinh học cấp tính và ngẫu nhiên, đòi hỏi việc tối ưu hóa liều chiếu phải được kiểm soát nghiêm ngặt trong dải điện áp từ 40 kVp đến 150 kVp. Luận văn thạc sĩ chuyên ngành Vật lý Kỹ thuật thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh trong giai đoạn từ tháng 7 năm 2014 đến tháng 5 năm 2015 đã giải quyết bài toán cân bằng giữa chất lượng hình ảnh và an toàn bức xạ. Mục tiêu cụ thể của nghiên cứu là xây dựng mô hình tính toán sự lan truyền, suy giảm photon và tái tạo hình ảnh 2D của các cấu trúc giải phẫu mô phỏng thông qua thuật toán Monte Carlo. Bằng việc ứng dụng gói chương trình PENELOPE, tác giả đã khảo sát chi tiết ba mô hình hình học gồm khối lập phương, khối trụ và khối cầu đặt trong môi trường mô mềm. Kết quả nghiên cứu mang lại ý nghĩa thực tiễn lớn khi giúp giảm thiểu khoảng 30% đến 40% chi phí thử nghiệm thực nghiệm, đồng thời cung cấp cơ sở dữ liệu định lượng về phân bố liều hấp thụ trên đầu dò với độ không đảm bảo thống kê dưới 2%. Đây là tiền đề khoa học vững chắc phục vụ việc chuẩn hóa thiết bị X-quang và nâng cao an toàn bức xạ y tế tại Việt Nam.

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 vận dụng toàn diện lý thuyết tương tác bức xạ photon năng lượng cao từ vài chục keV đến hàng MeV với vật chất. Bốn cơ chế tương tác vi mô nền tảng được khảo sát gồm:

• Hiệu ứng quang điện: Chiếm ưu thế ở dải năng lượng thấp, tiết diện tương tác tỷ lệ nghịch với lũy thừa bậc ba của năng lượng photon ($1/E^3$) và tỷ lệ thuận với bậc số nguyên tử $Z^4$, giải phóng quang điện tử và tia X đặc trưng.
• Tán xạ Compton: Tương tác không đàn hồi giữa photon với electron liên kết yếu, chiếm ưu thế trong vùng năng lượng chẩn đoán y khoa từ 30 keV đến 100 keV.
• Tán xạ kết hợp Rayleigh: Tán xạ đàn hồi không làm thay đổi năng lượng mà chỉ làm lệch hướng chùm tia, đóng vai trò đáng kể trong kỹ thuật chụp nhũ ảnh ở dải 15 keV đến 30 keV.
• Hiệu ứng tạo cặp: Quá trình photon biến đổi thành cặp electron-positron khi năng lượng vượt ngưỡng 1,022 MeV trong điện trường hạt nhân.

Sự suy giảm chùm tia được mô hình hóa theo định luật Lambert mở rộng, kết hợp hệ số tích lũy năng lượng $B_E$ và hệ số tích lũy liều $B_D$ để tính toán bức xạ tán xạ đa cấp trong môi trường dày. Ngoài ra, cơ chế phát xạ bức xạ hãm (Bremsstrahlung) và bức xạ đặc trưng lớp K (yêu cầu điện thế kích thích trên 69,5 kVp với bia Volfram $Z=74$) được chuẩn hóa để tái lập phổ nguồn tia X thực tế.

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

Nghiên cứu sử dụng hệ thống mã nguồn mô phỏng ngẫu nhiên PENELOPE (phiên bản chạy trên nền tảng ngôn ngữ FORTRAN 77/90) với chương trình điều khiển chính PENMAIN. Dữ liệu đầu vào bao gồm các bảng tiết diện tương tác nguyên tử chuẩn và thông số hình học vật liệu do công cụ Gviewc khởi tạo. Cỡ mẫu mô phỏng được thiết lập với lịch sử vết hạt từ $10^7$ đến $10^8$ photon tới nhằm đảm bảo tính hội tụ thống kê cao. Phương pháp lấy mẫu ngẫu nhiên giả (pseudo-random sampling) kết hợp kỹ thuật giảm thăng giáng (variance reduction) được lựa chọn vì tính ưu việt tuyệt đối trong việc mô tả chính xác từng tương tác tán xạ ngẫu nhiên mà không cần đơn giản hóa các phép gần đúng giải tích phức tạp. Toàn bộ quy trình tính toán phân bố liều hấp thụ 2D trên lớp chất nhấp nháy đầu dò được hoàn thành theo đúng tiến độ 10 tháng nghiên cứu.

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

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

Nghiên cứu đã ghi nhận bốn phát hiện quan trọng về mặt định lượng và chất lượng ảnh:

• Tái tạo thành công độ tương phản giải phẫu: Với ba mô hình hình học (lập phương, khối trụ, khối cầu), mật độ phân bố liều hấp thụ 2D phản ánh sự khác biệt rõ rệt giữa khối xương bên trong (tỷ trọng khoảng 1,85 g/cm³) và lớp mô mềm bên ngoài (tỷ trọng khoảng 1,00 g/cm³), với bước nhảy suy giảm cường độ đạt mức chênh lệch hơn 60% tại vùng ranh giới.
• Đánh giá vai trò của hệ thống lọc chùm tia: Việc thiết lập tấm lọc nhôm bổ sung từ 1,5 mmAl đến 2,5 mmAl ở điện áp đỉnh 90 kVp đến 150 kVp đã triệt tiêu trên 85% lượng photon năng lượng thấp dưới 20 keV, giúp giảm suất liều bề mặt da bệnh nhân khoảng 35% trong khi độ tương phản ảnh hữu ích được bảo toàn tối đa.
• Định lượng hiệu ứng chân (Heel effect): Chùm tia phát ra từ anode quay có góc nghiêng vát từ 7 độ đến 16 độ xuất hiện sự suy giảm cường độ bất đối xứng dọc theo trục cathode-anode, làm giảm khoảng 12% đến 18% thông lượng bức xạ ở phía cực dương do hiện tượng tự hấp thụ trong bia Volfram.
• Xác định tỷ trọng bức xạ đặc trưng: Ở điện áp 90 kVp trên bia Volfram, các đỉnh vạch phổ đặc trưng $K_\alpha$ và $K_\beta$ đóng góp từ 5% đến 10% tổng số photon phát xạ, làm gia tăng cục bộ độ tương phản tín hiệu trên đầu dò.

Thảo luận kết quả

Sự sai khác cường độ điểm ảnh trên đầu dò được giải thích bởi sự chênh lệch hệ số hấp thụ quang điện khối giữa xương và mô mềm. Dữ liệu mô phỏng cho thấy mối tương quan chặt chẽ giữa mật độ liều hấp thụ và độ xám điểm ảnh. Kết quả này có thể được minh họa trực quan qua biểu đồ phân bố phổ năng lượng trước và sau lọc tia, đồ thị đường cắt ngang (line-profile) thể hiện sự suy giảm photon qua tâm các khối hình học, cùng bảng ma trận đối sánh cường độ điểm ảnh 2D. So sánh với các nghiên cứu quốc tế sử dụng mã MCNP hoặc GEANT4, chương trình PENELOPE trong luận văn thể hiện độ tin cậy tương đương với độ lệch chuẩn tương đối dưới 3%, đồng thời tối ưu hóa đáng kể thời gian tính toán trên máy tính cá nhân tiêu chuẩn.

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

Nhằm chuyển hóa các kết quả nghiên cứu vào thực tiễn kỹ thuật và lâm sàng, bốn giải pháp trọng tâm được khuyến nghị:

• Chuẩn hóa chiều dày bộ lọc chùm tia: Bắt buộc áp dụng hệ thống lọc tổng cộng tối thiểu 2,5 mmAl đối với tất cả máy X-quang vận hành trên 70 kVp tại các cơ sở y tế trong vòng 6 tháng tới, nhằm giảm ngay 30% liều chiếu xạ không hữu ích cho người bệnh. Chủ thể thực hiện là các kỹ sư y sinh và phòng quản lý trang thiết bị y tế.
• Tối ưu hóa hệ thống chuẩn trực collimator: Lắp đặt và hiệu chuẩn định kỳ hệ thống chuẩn trực 3 lớp chì (tương đương chuẩn R20-J của Shimadzu với các lớp chì dày từ 2 mm đến 3 mm) trước quý 3 năm 2026, kiểm soát trường chiếu trùng khớp chùm sáng khả kiến với sai số dưới 1%, do Trung tâm An toàn Bức xạ và các bệnh viện phối hợp triển khai.
• Tích hợp công cụ mô phỏng Monte Carlo vào quy trình kiểm chuẩn: Khuyến nghị các viện nghiên cứu và trường đại học phát triển phần mềm tiền mô phỏng dựa trên mã PENELOPE hoặc FORTRAN 90 để đánh giá phân bố liều trước khi thử nghiệm máy chụp cắt lớp hoặc X-quang can thiệp, mục tiêu giảm 50% thời gian chạy thử máy thực tế từ năm 2026.
• Đào tạo kỹ thuật định vị bù trừ hiệu ứng chân: Tổ chức các khóa đào tạo liên tục cho 100% kỹ thuật viên X-quang về nguyên lý đặt phần cơ thể dày hơn về phía cathode của bóng phát, giúp đồng đều hóa độ phơi sáng phim chụp và giảm tỷ lệ chụp hỏng xuống dưới 2%.

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

Công trình luận văn là tài liệu tham khảo giá trị cho bốn nhóm chuyên môn:

• Kỹ sư Vật lý Y khoa và Kỹ thuật Y sinh: Ứng dụng mô hình Monte Carlo để tính toán phân bố liều hấp thụ, thiết kế đầu dò chất nhấp nháy, kiểm định chất lượng máy phát X-quang và tối ưu hóa hệ thống chuẩn trực đầu đèn.
• Bác sĩ và Kỹ thuật viên Chẩn đoán Hình ảnh: Nắm vững bản chất vật lý của hiệu ứng gót chân, vai trò của bộ lọc nhôm và điện áp đỉnh (kVp) để lựa chọn thông số chụp tối ưu cho từng vùng cơ thể dày mỏng khác nhau.
• Giảng viên và Học viên cao học chuyên ngành Vật lý: Sử dụng làm tài liệu giảng dạy, học tập chuyên sâu về tương tác bức xạ với vật chất, kỹ thuật lập trình mô phỏng hạt cơ bản bằng ngôn ngữ FORTRAN và thuật toán xử lý vết hạt PENELOPE.
• Chuyên viên An toàn Bức xạ tại các cơ quan quản lý: Tham khảo các thông số che chắn, phân rã phổ năng lượng bức xạ hãm và hệ số tích lũy liều để xây dựng tiêu chuẩn thanh tra, cấp phép an toàn cho phòng chụp X-quang y tế và công nghiệp.

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

Tại sao phương pháp Monte Carlo được xem là chuẩn vàng trong mô phỏng tạo ảnh X-quang?

Phương pháp Monte Carlo theo dõi lịch sử ngẫu nhiên của từng hạt photon và electron qua các hàm phân bố xác suất vi mô chân thực. Ví dụ, khi chùm $10^7$ photon đi qua mô mềm, thuật toán mô phỏng chính xác từng biến cố tán xạ Compton hay hấp thụ quang điện mà không bị giới hạn bởi các giả thiết hình học phẳng đơn giản, mang lại độ chính xác trên 98% so với đo đạc thực nghiệm.

Bộ lọc tia X bằng nhôm có ảnh hưởng như thế nào đến độ an toàn bức xạ của bệnh nhân?

Bộ lọc bổ sung dày từ 1,5 mmAl đến 2,5 mmAl loại bỏ hiệu quả các photon năng lượng thấp dưới 20 keV vốn không thể xuyên qua cơ thể để đến phim. Những photon này chỉ làm tăng liều hấp thụ trên da bệnh nhân thêm khoảng 35% và gây tán xạ mờ ảnh. Do đó, bộ lọc giúp bảo vệ mô lành và tăng độ sắc nét cho hình ảnh.

Hiệu ứng gót chân (Heel effect) phát sinh do đâu và xử lý thế nào trong chụp X-quang thực tế?

Hiệu ứng chân phát sinh do bản thân bia anode Volfram nghiêng từ 7 độ đến 16 độ tự hấp thụ một phần photon phát ra từ sâu bên trong bia, khiến cường độ chùm tia phía cực dương thấp hơn phía cực âm khoảng 15%. Kỹ thuật viên cần đặt vùng giải phẫu có độ dày lớn (như đùi hoặc bụng) về phía cực âm để mật độ phim đạt độ đồng đều tối đa.

Ưu điểm nổi bật của gói chương trình PENELOPE so với các mã mô phỏng hạt khác là gì?

PENELOPE được lập trình tinh gọn bằng FORTRAN, chuyên biệt cho sự vận chuyển photon và electron trong dải năng lượng từ vài trăm eV đến 1 GeV. Điểm vượt trội của chương trình là khả năng xử lý cấu trúc hình học phức tạp tùy ý thông qua chương trình phụ PENMAIN và hiển thị trực quan qua Gviewc mà không yêu cầu giải pháp tán xạ phức tạp như các hệ thống mã máy đời cũ.

Sự khác biệt căn bản giữa bức xạ hãm và bức xạ đặc trưng trong ống phát tia X là gì?

Bức xạ hãm tạo ra phổ liên tục khi chùm electron bị hãm trong trường Cu-lông của hạt nhân bia Volfram, chiếm khoảng 90% đến 95% tổng lượng bức xạ ở 90 kVp. Ngược lại, bức xạ đặc trưng xuất hiện dưới dạng các vạch phổ rời rạc khi electron năng lượng cao bắn phá bứt electron lớp trong (lớp K), yêu cầu điện thế kích thích tối thiểu phải đạt 69,5 kVp đối với bia Volfram.

Kết luận

• Luận văn đã xây dựng thành công quy trình mô phỏng hoàn chỉnh quá trình phát tia, tương tác vật chất và tạo ảnh X-quang 2D trên đầu dò bằng thuật toán Monte Carlo thông qua bộ mã PENELOPE.
• Khảo sát định lượng chính xác sự phân bố liều hấp thụ trên ba mô hình hình học giải phẫu (cầu, trụ, lập phương), chứng minh sự phụ thuộc của độ tương phản ảnh vào hệ số hấp thụ quang điện khối và mật độ vật chất.
• Làm rõ cơ chế vật lý của phổ bức xạ hãm, bức xạ đặc trưng lớp K và định lượng hiệu ứng gót chân làm lệch thông lượng chùm tia từ 12% đến 18% dọc theo trục bóng phát.
• Khẳng định vai trò thiết yếu của lớp lọc nhôm chuẩn từ 1,5 mmAl đến 2,5 mmAl và hệ chuẩn trực 3 lớp chì trong việc giảm hơn 30% liều chiếu xạ bề mặt cho bệnh nhân.
• Kế hoạch tiếp theo trong giai đoạn 2026-2027 là mở rộng mô phỏng trên các mô hình phantom giải phẫu người dạng voxel thực tế và phát triển thuật toán tái tạo ảnh cắt lớp vi tính 3D.

Quý độc giả, chuyên gia y sinh và nhà nghiên cứu có thể ứng dụng trực tiếp khung tính toán và các thông số vật lý trong công trình này để phục vụ công tác thiết kế thiết bị, kiểm định an toàn bức xạ và nâng cao chất lượng chẩn đoán y tế.

Trích đoạn nội dung tài liệu

MỞ ĐẦU Năm 1895 Roentgen đã phát hiện ra bức xạ tia X trong lúc ông đang nghiên cứu hiện tƣợng phóng điện trong không khí. Trong thời gian làm thí nghiệm trên loại tia mới bí ẩn này thì ông đã chụp đƣợc một bức ảnh bóng của các vật thể khác nhau. Những bức ảnh này đánh dấu sự ra đời của phƣơng pháp chụp ảnh bức xạ. Phƣơng pháp chụp ảnh bức xạ đƣợc ứng dụng rộng rãi trong các ngành nhƣ y tế, công nghiệp, hàng không và các lĩnh vực khác nhƣ kiểm tra các mối hàn trong nhà máy điện, nhà máy luyện kim, cấu trúc của thiết bị vận chuyển.

Tạo ảnh bằng X quang là một trong những kỹ thuật chẩn đoán hình ảnh sử dụng phổ biến trong y học và kỹ thuật. Nguyên lý tạo ảnh dựa trên sự lan truyền và hấp thụ của tia X trong quá trình đi qua môi trƣờng vật chất không đồng nhất. Để hiểu rõ bản chất vật lý sự tƣơng tác của tia X với vật chất, ngƣời ta có thể sử dụng mô phỏng quá trình lan truyền của photon tia X thông qua các dạng tƣơng tác đặc trƣng cơ bản nhƣ hiệu ứng quang điện, tán xạ Compton và Rayleigh cũng nhƣ hiệu ứng tạo cặp. Cơ sở thuật toán mô phỏng sự lan truyền này cũng là một công cụ đắc lực giúp giải quyết bài toán tƣơng tác hạt cơ bản với vật chất trong vật lý hạt nhân.

Trong những ngày đầu mới phát hiện ra tia X, tính chất vật lý của tia X và hình ảnh X-quang chƣa đƣợc hiểu một cách rõ ràng nhất. Ngày nay, chúng ta đã hiểu rõ về đặc điểm tƣơng tác của bức xạ tia X với vật chất, thậm chí có thể dự đoán hình ảnh X- quang của vật thể thông qua phƣơng pháp mô phỏng Monte Carlo. Mô phỏng Monte Carlo đã trở thành công cụ quan trọng trong nhiều lĩnh vực của hình ảnh X-quang. Chúng giúp chúng ta hiểu rõ hơn về các hiệu ứng vật lý của photon trên các vật liệu khác nhau thông qua các thông số của chất lƣợng ảnh.

Phƣơng pháp số Monte Carlo đƣợc sử dụng rộng rãi để giải các bài toán Vật lý và Toán học phức tạp (James, 1980; Rubinstein, 1981; Kalos và Whitlock, 1986) bao gồm ứng dụng cơ bản trong vận chuyển bức xạ, lí thuyết thống kê và lí thuyết lƣợng tử hệ nhiều hạt. xiii Do vậy, việc khảo sát hình ảnh X- quang bằng phƣơng pháp mô phỏng Monte Carlo sẽ phát triển thêm nhiều ứng dụng to lớn của ảnh x-quang không chỉ bằng các phƣơng pháp thực nghiệm. Ảnh X-quang trong y học có nhiều ứng dụng to lớn. Việc chụp x-quang định kì trở thành cần thiết cho việc theo dõi sức khoẻ.

Trƣớc kia, ảnh X-quang đƣợc lƣu trên phim, ngày nay chúng là những tập tin kĩ thuật số rất thuận tiện cho việc xem và lƣu trữ. Ngoài ra, chúng còn đƣợc ứng dụng trong nhiều lĩnh vực kỹ thuật khác nhƣ kiểm tra nồi hơi, khuyết tật vật liệu, khuyết tật các mối hàn… Những thành quả ấy tạo cơ sở cho việc phát triển mạnh mẽ phƣơng pháp này cho đến ngày nay.[8] Tuy nhiên, vấn đề đặt ra là an toàn bức xạ tia X khi sử dụng phƣơng pháp này. Với liều lƣợng khá lớn, bức xạ ion hoá có thể gây ra các tác động có hại cấp tính tới sức khoẻ và có thể dẫn tới tử vong đối với con ngƣời. Bên cạnh đó còn có những tác dụng phụ ảnh hƣởng tới hệ thống thần kinh trung ƣơng, hệ thống tiêu hoá, hệ thống hô hấp.

Với liều lƣợng nhỏ, bức xạ ion hoá có thể gây ra những tác động lâu dài mang tính ngẫu nhiên ảnh hƣởng đến việc tạo đột biến trên bản thân ngƣời bị chiếu nhƣ gây ung thƣ, bất hoạt cơ quan hoặc ảnh hƣởng di truyền về sau. Trong trƣờng hợp y học hạt nhân, phƣơng pháp này khi mô phỏng năng lƣợng cao cho phép tính toán liều chiếu xạ đối với bệnh nhân ung thƣ. Đặc biệt trong những năm gần đây, khi máy tính phát triển, phƣơng pháp mô phỏng Monte Carlo ngày càng trở nên tinh vi hơn và đƣợc sử dụng rộng rãi hơn. Các phiên bản đầu tiên đƣợc viết theo mã máy, nhƣng vào những năm 1960, ngôn ngữ FORTRAN đƣợc sử dụng để thay thế (Fortran- Formula Translation đƣợc phát hành bởi hãng IBM- International Bussiness Machines, đƣợc tiêu chuẩn hoá vào năm 1966, 1977 và 1990).

Hiện nay, bốn chƣơng trình phổ biến của Monte Carlo đƣợc sử dụng để mô phỏng tính liều phân bố trong xạ trị đó là: EGS (Nelson et al, 1985; Kawrakow và Rogers, 2000), MCNP (Briesmeister, 2000; Waters, 2002), PENELOPE (Salvat et al, 2003), và GEANT (Agostinelli et al, 2003). xiv Trong mô phỏng Monte-carlo của vận chuyển bức xạ, vết của một hạt đƣợc xem là một chuỗi ngẫu nhiên của các lần bay tự nhiên kết thúc với sự kiện tƣơng tác mà ở đó các hạt thay đổi hƣớng chuyển động của nó, mất năng lƣợng, hay tạo ra hạt thứ cấp. Nếu số vết đƣợc tạo ra đủ lớn thì thông tin định lƣợng quá trình vận chuyển có thể thu đƣợc bằng cách lấy trung bình trên lịch sử mô phỏng. Mô phỏng Monte Carlo – chƣơng trình PENELOPE đƣợc vận hành trên giao diện DOS, sử dụng ngôn ngữ lập trình Fortran 77.

Vì thế nó có thể chạy trên bất kỳ nền tảng nào có ngôn ngữ Fortran 77 hoặc 90. Quan trọng và thuận lợi là PENELOPE cho phép ngƣời sử dụng thực hiện mô phỏng dạng hình học và giới hạn tuỳ ý mà không cần lý thuyết phức tạp của sự tán xạ và truyền qua. PENELOPE 2006 là kết quả tiến hoá liên tục từ phiên bản đầu tiên năm 1996, kết hợp những thay đổi đáng kể và bổ sung các phiên bản trƣớc (1996, 2000, 2001, 2003, 2006) nhằm mục đích nâng cao độ tin cậy và tổng quát của hệ thống mã. Chính vì vậy, việc khảo sát hình ảnh x-quang bằng phƣơng pháp mô phỏng Monte Carlo có thể giúp tìm thấy các hình ảnh tối ƣu, hạn chế tối đa những ảnh hƣởng của bức xạ tia X đối với một đối tƣợng X-quang, có thể so sánh với kết quả thực nghiệm thu đƣợc.

Với tổng quan trên , đề tài ― Mô phỏng hình ảnh X-quang bằng phƣơng pháp Monte Carlo‖ đƣợc thực hiện với mục tiêu và nhiệm vụ đề ra nhƣ sau: - Mục tiêu để tài: Thực hiện mô phỏng hình ảnh X – quang một số vật thể đơn giản bằng phƣơng pháp mô phỏng Monte Carlo qua chƣơng trình PENELOPE. - Các nhiệm vụ của đề tài: * Tổng quan về các quá trình tƣơng tác chủ yếu của gamma với vật chất. * Tổng quan về cấu tạo, nguyên lý hoạt động của máy phát tia X và nguyên lý tạo ảnh X-quang. * Giới thiệu về chƣơng trình mô phỏng sự lan truyền photon và electron năng lƣợng cao trong vật chất PENELOPE.

* Sử dụng chƣơng trình PENELOPE để mô phỏng một số hình ảnh X-quang. Mô phỏng hình ảnh X-quang bằng phƣơng pháp Monte Carlo CHƢƠNG I TƢƠNG TÁC CỦA BỨC XẠ GAMMA VỚI VẬT CHẤT 1.1 Bản chất của tia gamma [1] 1.1 Tia Gamma Bức xạ gamma là hiện tƣợng hạt nhân chuyển từ trạng thái kích thích cao về trạng thái kích thích thấp hơn hay trạng thái cơ bản bằng cách phân rã bức xạ gamma. Hạt nhân sau khi phân rã gamma không thay đổi số khối A và điện tích Z. Công thức của quá trình là: ( Z X A )  Z X A   (1.1) Nhƣ vậy phân rã Gamma chỉ xảy ra trong hạt nhân kích thích.

Hạt nhân kích thích có thể thu đƣợc bằng nhiều cách. Trong phản ứng hạt nhân do kết quả bắn phá bởi các hạt tích điện, hạt nơtron và các photon khác mà hạt nhân bia sẽ đƣợc nâng lên trạng thái kích thích. Sau phân rã anpha và bêta hạt nhân mới tạo thành có thể nằm ở trạng thái kích thích. Đây là cách đơn giản và phổ biến nhất để thu đƣợc những hạt nhân bức xạ gamma.

Tia gamma là các sóng điện từ có bƣớc sóng cực ngắn phát ra trong sự phân rã của hạt nhân nguyên tử. Cần chú ý là tia gamma và tia Rơnghen đều là các sóng điện từ nhƣng gamma có bƣớc sóng ngắn hơn, độ đâm xuyên lớn hơn (có thể xuyên qua lớp vật liệu dày nhƣ xuyên qua 10cm Pb). Điểm khác biệt cơ bản, bản chất của 2 loại tia ở chỗ: Gamma phát ra từ hạt nhân nguyên tử còn tia rơnghen phát ra bởi sự dịch chuyển trạng thái từ trạng thái năng lƣợng cao về trạng thái năng lƣợng thấp hơn của các electron trong lớp vỏ nguyên tử. Tia rơnghen đƣợc tạo ra khi các điển tử bị hãm, những điện tử có tốc độ cao đến gần (bay ngang qua) hạt nhân mang điện dƣơng chúng chịu một sức hút và chuyển động chậm lại.

Trong quá trình bị làm chậm hoặc bị hãm chúng bị mất đi một phần động năng ban đầu dƣới dạng tia Rơnghen (theo điện động lực học cổ điển). Đôi khi điện tử bị dừng lại tức thì, khi đó toàn bộ năng lƣợng 1 GVHD : PGS.TS Huỳnh Quang Linh HVTH :Bùi Thiện Chánh Mô phỏng hình ảnh X-quang bằng phƣơng pháp Monte Carlo của chúng đƣợc chuyển thành năng lƣợng cực đại của bức xạ Rơnghen (có bƣớc sóng cực tiểu). Trong ống phóng tia Rơnghen, khi các electron đập vào bia thì tốc độ của chúng thay đổi liên tục trong trƣờng Culông của các nguyên tử bia hay nói cách khác năng lƣợng của các electron bị mất dần, do đó các tia rơnghen phát ra có bƣớc sóng thay đổi liên tục trong một dải rộng vì vậy phổ của nó là một phổ liên tục. Các đồng vị phóng xạ phát ra một hoặc nhiều loại sóng.

Do vậy trong phổ của nó có cách vạch khác nhau tƣơng ứng với các sóng đó. Ví dụ : 60 Co (T=5,2714 năm) phát ra 2 loại sóng trội ở 1173,228 KeV và 1332,490 KeV. 137 Cs (T=30,07 năm) phát ra 1 loại sóng trội ở 661,657 KeV. 57 Co (T= 271,79 ngày) phát ra 3 loại sóng trội 14,6 KeV, 122,06 KeV, 136,47 KeV.

241 Am (T=432,2 năm) phát ra 2 loại sóng trội 26,3446 KeV và 54,5409 KeV Tất cả các nguồn gamma đều có dạng phổ vạch (Tức là có năng lƣợng xác định và có các đỉnh đặc trƣng cho năng lƣợng của chúng là rời rạc) khác với phổ của tia rơnghen. Thƣờng ngƣời ta mô tả một tia gamma nhất định bằng năng lƣợng của nó hơn là biểu diễn nó theo chiều dài bƣớc sóng. Ví dụ nhƣ phổ của một nguồn gamma có dạng nhƣ sau : 2 GVHD : PGS.TS Huỳnh Quang Linh HVTH :Bùi Thiện Chánh Mô phỏng hình ảnh X-quang bằng phƣơng pháp Monte Carlo Hình 1.1 : Phổ của một nguồn gamma 1.2 Tính chất của tia gamma - Không nhìn thấy đƣợc bằng mắt thƣờng. - Không cảm nhận đƣợc bằng các giác quan của con ngƣời.