Khóa Luận Tốt Nghiệp Về Mô Phỏng Tương Tác Tạo Ảnh X-Quang Bằng Phương Pháp Monte Carlo

Luận văn tốt nghiệp vật lý nghiên cứu tốt nghiệp vật lý mô phỏng quá trình tương tác tạo ảnh x quang bằng phương pháp monte carlo, điều tra thực trạng, phân tích số liệu, đề xuất

• Tác giả: Phùng Nhật Anh
• Người hướng dẫn: TS. Đỗ Xuân Hội
• Trường học: Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh
• Chuyên ngành: Vật Lý
• Đề tài: Mô Phỏng Quá Trình Tương Tác Tạo Ảnh X-Quang Bằng Phương Pháp Monte Carlo
• Loại tài liệu: luận văn tốt nghiệp
• Năm xuất bản: 2000-2004
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI NÓI ĐẦU
• 1. CHƯƠNG 1: TƯƠNG TÁC TIA X VỚI VẬT CHẤT
• 1.1. Tán xạ Rayleigh
• 1.2. Tán xạ Compton
• 2. CHƯƠNG 2: SỰ SUY GIẢM VÀ SỰ HẤP THỤ NĂNG LƯỢNG CỦA TIA X KHI ĐI QUA MÔI TRƯỜNG VẬT CHẤT
• 2.1. Hệ số suy giảm tuyến tính
• 2.2. Hệ số suy giảm khối lượng
• 2.3. Năng lượng hiệu dụng
• 2.4. Quang đường tự do trung bình
• 2.5. Phương pháp làm cứng chùm tia
• 2.6. Sự hấp thụ năng lượng tia X và tia gamma
• 2.7. Liều lượng sinh học của tia X
• 3. CHƯƠNG 3: CƠ SỞ CHỤP ẢNH X QUANG BẰNG PHIM
• 3.1. Các phương pháp chẩn đoán hình ảnh y khoa
• 3.2. Phương pháp X quang truyền qua (projection radiography)
• 3.3. Nguyên tắc hình học cơ bản của việc tạo ảnh
• 3.4. Băng đựng phim (screen - film cassette)
• 3.5. Bức xạ tán xạ trong ngành X quang truyền qua
• 4. CHƯƠNG 4: NGUỒN TẠO X - QUANG
• 4.1. Góc anốt và kích cỡ vệt chiếu
• 4.2. Bức xạ ngoài chùm (off-focus radiation)
• 4.3. Ống chuẩn trực (collimator)
• 4.4. Mạch điện trong thiết bị tạo X - quang
• 5. CHƯƠNG 5: PHỔ NĂNG LƯỢNG BỨC XẠ X - QUANG
• 5.1. Phổ bức xạ hãm (Bremsstrahlung)
• 5.2. Quang phổ tia X đặc trưng (characteristic X-Ray Spectrum)
• 6. CHƯƠNG 6: CƠ SỞ VỀ PHƯƠNG PHÁP MONTE-CARLO
• 6.1. Giới thiệu về phương pháp Monte-Carlo
• 6.2. Số ngẫu nhiên - Giả ngẫu nhiên
• 6.2.1. Phương pháp đồng dư tuyến tính
• 6.2.2. Các phương pháp mô phỏng một phân bố bất kì
• 6.2.2.1. Phương pháp chấp nhận - Loại bỏ (rejection technique)
• 6.2.2.2. Phương pháp hàm ngược
• 6.2.2.3. Phương pháp lấy mẫu (importance sampling)
• 7. PHẦN II: KẾT QUẢ THỰC HIỆN MÔ PHỎNG
• 7.1. Mục tiêu của chương trình
• 7.2. Xác định ngôn ngữ lập trình
• 7.3. Tạo dữ liệu mô phỏng
• 7.4. Lưu đồ của chương trình
• 7.5. Xây dựng chương trình mô phỏng
• 7.6. Mô tả các modul của chương trình mô phỏng
• 7.7. Kết quả và biện luận
• 8. PHẦN III: KẾT LUẬN
• PHỤ LỤC
• Phụ lục 1: Chương trình tạo dữ liệu mô phỏng
• Phụ lục 2: Thành phần và mật độ của một số vật chất sinh học
• Phụ lục 3: Chương trình mô phỏng
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về khóa luận tốt nghiệp vật lý mô phỏng

Khóa luận tốt nghiệp này tập trung vào việc mô phỏng quá trình tương tác tạo ảnh X quang bằng phương pháp Monte Carlo. Tia X đã trở thành một công cụ quan trọng trong chẩn đoán y khoa, và việc hiểu rõ quá trình tương tác của nó với vật chất là rất cần thiết. Mục tiêu của nghiên cứu là khảo sát các phương pháp mô phỏng để tối ưu hóa chất lượng ảnh X quang và đảm bảo an toàn cho bệnh nhân.

1.1. Lịch sử phát triển của tia X trong y học

Tia X được phát hiện vào năm 1895 bởi Wilhelm Conrad Roentgen. Phát hiện này đã mở ra một kỷ nguyên mới trong chẩn đoán y khoa, giúp bác sĩ có thể nhìn thấy bên trong cơ thể mà không cần phẫu thuật.

1.2. Vai trò của phương pháp Monte Carlo trong mô phỏng

Phương pháp Monte Carlo cho phép mô phỏng các quá trình ngẫu nhiên, giúp phân tích chính xác các tương tác của photon tia X với vật chất. Điều này rất quan trọng trong việc tối ưu hóa chất lượng ảnh X quang.

II. Vấn đề và thách thức trong quá trình tạo ảnh X quang

Mặc dù công nghệ chụp ảnh X quang đã phát triển mạnh mẽ, nhưng vẫn còn nhiều thách thức trong việc đảm bảo chất lượng ảnh và an toàn cho bệnh nhân. Các vấn đề như tán xạ không mong muốn và hấp thụ năng lượng cần được giải quyết.

2.1. Tán xạ và ảnh hưởng đến chất lượng ảnh

Tán xạ Rayleigh và Compton có thể làm giảm độ tương phản của ảnh X quang. Việc hiểu rõ các loại tán xạ này là cần thiết để cải thiện chất lượng ảnh.

2.2. An toàn cho bệnh nhân trong chẩn đoán

Liều lượng bức xạ cần được kiểm soát chặt chẽ để đảm bảo an toàn cho bệnh nhân. Việc sử dụng phương pháp mô phỏng giúp xác định liều lượng tối ưu cho từng loại chẩn đoán.

III. Phương pháp mô phỏng Monte Carlo trong tạo ảnh X quang

Phương pháp Monte Carlo là một công cụ mạnh mẽ trong việc mô phỏng quá trình tương tác của photon tia X với vật chất. Nó cho phép phân tích chi tiết các tương tác và tối ưu hóa các thông số chẩn đoán.

3.1. Nguyên lý hoạt động của phương pháp Monte Carlo

Phương pháp này sử dụng các số ngẫu nhiên để mô phỏng các quá trình tương tác phức tạp, từ đó giúp dự đoán kết quả chính xác hơn trong việc tạo ảnh X quang.

3.2. Ứng dụng của Monte Carlo trong nghiên cứu

Nghiên cứu đã chỉ ra rằng việc áp dụng phương pháp Monte Carlo có thể cải thiện đáng kể chất lượng ảnh X quang, đồng thời giảm thiểu liều lượng bức xạ cho bệnh nhân.

IV. Kết quả nghiên cứu và ứng dụng thực tiễn

Kết quả từ mô phỏng cho thấy rằng việc tối ưu hóa các thông số chẩn đoán có thể nâng cao chất lượng ảnh X quang. Nghiên cứu cũng chỉ ra rằng phương pháp Monte Carlo có thể được áp dụng rộng rãi trong các lĩnh vực khác nhau của y học.

4.1. Kết quả mô phỏng và phân tích dữ liệu

Dữ liệu thu được từ mô phỏng cho thấy sự cải thiện rõ rệt trong chất lượng ảnh X quang, với độ tương phản và độ phân giải cao hơn.

4.2. Ứng dụng trong chẩn đoán y khoa

Kết quả nghiên cứu có thể được áp dụng trong thực tiễn để nâng cao hiệu quả chẩn đoán, giúp bác sĩ đưa ra quyết định chính xác hơn.

V. Kết luận và triển vọng tương lai của nghiên cứu

Nghiên cứu đã chỉ ra rằng phương pháp Monte Carlo là một công cụ hữu ích trong việc mô phỏng quá trình tương tác tạo ảnh X quang. Tương lai của nghiên cứu này có thể mở ra nhiều hướng đi mới trong lĩnh vực chẩn đoán y khoa.

5.1. Hướng phát triển tiếp theo trong nghiên cứu

Cần tiếp tục nghiên cứu để cải thiện hơn nữa các thuật toán mô phỏng, nhằm nâng cao độ chính xác và hiệu quả trong việc tạo ảnh X quang.

5.2. Tác động của nghiên cứu đến ngành y tế

Nghiên cứu này không chỉ có ý nghĩa trong việc nâng cao chất lượng ảnh X quang mà còn góp phần bảo vệ sức khỏe cho bệnh nhân, giảm thiểu rủi ro từ bức xạ.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật chẩn đoán hình ảnh y khoa hiện đại, tia X đóng vai trò là công cụ cốt lõi phục vụ thăm khám giải phẫu và phát hiện bệnh lý. Tuy nhiên, việc sử dụng bức xạ ion hóa trong dải năng lượng từ 10 keV đến 300 keV luôn đặt ra bài toán tối ưu hóa đa mục tiêu: vừa nâng cao chất lượng tương phản và độ phân giải hình ảnh, vừa giảm thiểu liều hấp thụ bức xạ nhằm bảo vệ an toàn cho bệnh nhân theo nguyên lý ALARA. Đối với các kỹ thuật chuyên sâu như chụp nhũ ảnh ở dải năng lượng thấp 15 keV đến 30 keV hay chụp X-quang tổng quát từ 20 keV đến 150 keV, sự tương tác phức tạp của photon với các mô sinh học tạo ra bức xạ tán xạ làm suy giảm đáng kể chất lượng chẩn đoán.

Vấn đề cốt lõi đặt ra là làm thế nào để khảo sát định lượng chi tiết từng tiến trình tương tác vi mô của photon bên trong các môi trường vật chất không đồng nhất mà không cần thực hiện các phép chiếu chụp thử nghiệm lặp lại trên cơ thể người. Đề tài luận văn tập trung nghiên cứu, xây dựng thuật toán và chương trình mô phỏng quá trình tương tác tạo ảnh X-quang bằng phương pháp Monte Carlo. Nghiên cứu được triển khai tại Khoa Vật lý Trường Đại học Sư phạm Thành phố Hồ Chí Minh phối hợp cùng Trường Đại học Bách Khoa Thành phố Hồ Chí Minh.

Ý nghĩa thực tiễn của công trình thể hiện ở việc thiết lập mô hình tính toán số chuẩn xác, hỗ trợ xác định các thông số phát tia tối ưu cho từng vùng giải phẫu. Kết quả nghiên cứu cung cấp cơ sở khoa học để thiết kế các bộ lọc làm cứng chùm tia và lưới chống tán xạ, giúp giảm ước tính từ 20% đến 30% liều chiếu trên da bệnh nhân trong khi vẫn duy trì tỷ số tín hiệu trên nhiễu ở mức tối ưu cho việc chẩn đoán lâm sàng.

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ạ photon với vật chất trong cơ học lượng tử và quang học hạt nhân, kết hợp lý thuyết suy giảm hàm mũ chùm tia bức xạ và các quy chuẩn đo liều sinh học quốc tế.

Các cơ chế tương tác vi mô chủ đạo của chùm photon gồm:

• Tán xạ Rayleigh (tán xạ kết hợp): Xảy ra khi photon tương tác đàn hồi với toàn bộ nguyên tử mà không làm ion hóa, tiết diện vi phân tán xạ được xác định thông qua hàm phân bố Thompson nhân với thừa số hình thức nguyên tử.
• Tán xạ Compton (tán xạ không kết hợp): Tương tác giữa photon và electron lớp ngoài cùng, được mô tả bởi công thức tiết diện Klein-Nishina, làm photon mất một phần năng lượng và đổi hướng chuyển động.
• Hiệu ứng quang điện: Toàn bộ năng lượng photon truyền cho electron liên kết (chủ yếu lớp K với tỷ lệ khoảng 80%), xác suất hấp thụ tỷ lệ thuận với bậc ba số hiệu nguyên tử và tỷ lệ nghịch với bậc ba năng lượng photon, thể hiện các bước nhảy vọt tại ngưỡng hấp thụ K-edge như đối với Iot tại mức 33.2 keV.
• Phản ứng tạo cặp: Cần ngưỡng năng lượng tối thiểu 1.02 MeV nên hoàn toàn không xảy ra trong phạm vi X-quang chẩn đoán.

Các khái niệm đo lường quang học và liều lượng bao gồm: Hệ số suy giảm tuyến tính, hệ số suy giảm khối lượng, lớp bán trị, quãng đường tự do trung bình, thông lượng năng lượng, liều hấp thụ đo bằng đơn vị Gray và liều lượng hiệu dụng đo bằng Sievert theo tiêu chuẩn ICRP.

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

Luận văn ứng dụng phương pháp mô phỏng số ngẫu nhiên Monte Carlo để giải bài toán vận chuyển photon trong môi trường vật chất sinh học:

• Quy mô mẫu mô phỏng: Tập mẫu được thiết lập khảo sát từ 100.000 đến 1.000.000 lịch sử hạt photon đơn sắc và đa sắc nhằm khống chế sai số thống kê ngẫu nhiên xuống dưới mức 1.5%.
• Phương pháp chọn mẫu và phân tích: Chương trình sử dụng bộ tạo số giả ngẫu nhiên đồng dư tuyến tính chất lượng cao. Thuật toán tích hợp linh hoạt ba phương pháp lấy mẫu phân bố xác suất gồm kỹ thuật hàm phân bố nghịch đảo để xác định quãng đường tự do trung bình, phương pháp chấp nhận - loại bỏ để tạo phổ bức xạ liên tục Bremsstrahlung và góc tán xạ, kết hợp kỹ thuật lấy mẫu quan trọng nhằm tập trung khảo sát các photon truyền đến bề mặt nhận ảnh.
• Lý do lựa chọn: Monte Carlo là phương pháp số duy nhất cho phép theo dõi chi tiết từng biến cố va chạm vi mô của hạt, tính toán phân bố năng lượng lắng đọng trong cấu trúc giải phẫu nhiều lớp (mô mềm bao quanh xương) mà các phương trình giải tích cổ điển không thể giải quyết trọn vẹn. Quy trình mô phỏng được thực thi qua các module tính toán độc lập viết bằng ngôn ngữ lập trình chuyên dụng và kết xuất đồ họa qua môi trường MATLAB.

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 số và phân tích thực nghiệm đã mang lại những kết quả định lượng cụ thể sau:

• Đặc trưng phân bố các loại tương tác theo mức năng lượng: Trong môi trường mô mềm, ở mức năng lượng 30 keV, tán xạ Rayleigh chiếm tỷ trọng khoảng 12% tổng số tương tác, nhưng khi năng lượng chùm photon tăng lên trên 70 keV, tỷ lệ này giảm nhanh xuống dưới 5%. Tán xạ Compton bắt đầu chiếm ưu thế tuyệt đối khi năng lượng vượt ngưỡng 26 keV ở mô mềm và vượt ngưỡng 35 keV ở xương, khẳng định đây là nguồn gốc chính sinh ra bức xạ thứ cấp làm mờ ảnh.
• Mức độ suy giảm tương phản do bức xạ tán xạ: Khi tỷ số tán xạ trên sơ cấp đạt mức bằng 3, có tới 75% lượng photon đến đầu dò là photon tán xạ không mang thông tin cấu trúc giải phẫu. Kết quả tính toán cho thấy độ tương phản hình ảnh thực tế bị suy giảm từ mức ban đầu 50% xuống chỉ còn 12.5%, tức giảm 4 lần so với điều kiện không có tán xạ.
• Tương tác trong cấu trúc không đồng nhất mô mềm và xương: Khảo sát mô hình lớp xương lồng trong khối mô mềm chỉ ra hiệu ứng quang điện ở xương cao gấp 6 đến 8 lần so với mô mềm ở dải năng lượng dưới 50 keV do sự chênh lệch số hiệu nguyên tử hiệu dụng, tạo nên độ tương phản sắc nét trên biên giới hạn giữa xương và mô nhưng đồng thời tạo ra gradient tán xạ gây mờ vùng biên khoảng 15% đến 20%.
• Hiệu năng của các loại màn khuếch đại nhấp nháy: Bản khuếch đại đất hiếm Gadolinium Oxysulfide pha tạp Terbi đạt hiệu suất chuyển đổi năng lượng tia X thành ánh sáng khả kiến ở mức 15%, cao gấp 3 lần so với chất nhấp nháy cổ điển Canxi Vonfamat đạt 5%. Khi sử dụng hệ nhũ tương kép với độ dày màn 120 mg/cm², liều chiếu cần thiết để đạt mật độ quang học tiêu chuẩn từ 1.0 đến 1.2 giảm khoảng 60%.

Thảo luận kết quả

Cơ chế vật lý giải thích cho sự suy giảm chất lượng ảnh là do tiết diện tán xạ Compton giảm rất chậm theo năng lượng so với tiết diện hấp thụ quang điện. Khi tăng điện áp phát tia để tăng độ đâm xuyên, các photon tán xạ Compton mang năng lượng cao vẫn tiếp tục truyền thẳng hoặc lệch góc nhỏ đập vào phim, tạo thành phông bức xạ nền đồng nhất.

Dữ liệu mô phỏng được biểu diễn trực quan qua các ma trận phân bố mật độ quang học 2D và biểu đồ 3D trong MATLAB, thể hiện rõ vùng bóng của cấu trúc xương trên nền mô. Biểu đồ đường cong đặc trưng Hurter & Driffield cho thấy gradient tương phản trung bình đạt cực đại trong khoảng mật độ quang học từ 2.5 đến 3.5. Kết quả mô phỏng phân bố năng lượng và hệ số suy giảm khối lượng khớp hoàn toàn với các bảng dữ liệu thực nghiệm chuẩn của ICRP, chứng minh thuật toán Monte Carlo được xây dựng có độ tin cậy rất cao.

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

Dựa trên kết quả nghiên cứu mô phỏng, bốn giải pháp kỹ thuật cụ thể được khuyến nghị áp dụng trong thực tiễn chẩn đoán hình ảnh:

1. Chuẩn hóa thông số phát tia cho từng loại chẩn đoán lâm sàng:

   • Hành động: Thiết lập điện áp bóng phát và dòng chụp phù hợp; duy trì mức 25 kVp đến 30 kVp cho chụp nhũ ảnh và 80 kVp đến 120 kVp cho chụp ngực tổng quát.
   • Mục tiêu: Giảm 25% đến 30% liều hấp thụ bức xạ trên da bệnh nhân trong vòng 3 tháng áp dụng tại các khoa chẩn đoán hình ảnh.
   • Chủ thể thực hiện: Bác sĩ chẩn đoán hình ảnh và kỹ thuật viên X-quang.

2. Tối ưu hóa hình học chụp và sử dụng lưới chống tán xạ:

   • Hành động: Lắp đặt lưới chống tán xạ có tỷ số lá chì từ 8:1 đến 12:1 cho các bộ phận có độ dày giải phẫu trên 15 cm; đồng thời thu hẹp trường chiếu bằng ống chuẩn trực sát vùng cần khảo sát.
   • Mục tiêu: Kéo giảm tỷ số tán xạ trên sơ cấp xuống dưới 1.0, tăng độ tương phản hình ảnh thêm 40% ngay trong quý tiếp theo.
   • Chủ thể thực hiện: Kỹ sư vật lý y sinh và kỹ thuật viên vận hành máy.

3. Hiện đại hóa hệ thống nhận ảnh và tấm chuyển đổi quang học:

   • Hành động: Thay thế toàn bộ màn Canxi Vonfamat cũ bằng màn phosphor đất hiếm hoặc đầu dò kỹ thuật số dạng tấm phẳng.
   • Mục tiêu: Nâng cao hiệu suất hấp thụ lượng tử lên trên 65% theo lộ trình 12 tháng tại các cơ sở y tế.
   • Chủ thể thực hiện: Ban giám đốc bệnh viện và phòng trang thiết bị y tế.

4. Ứng dụng phần mềm mô phỏng vào nghiên cứu và đào tạo chuyên sâu:

   • Hành động: Triển khai mô hình tính toán Monte Carlo vào chương trình đào tạo Vật lý Y khoa và các dự án thiết kế thiết bị chẩn đoán.
   • Mục tiêu: Hoàn thiện bộ công cụ tính toán liều lượng tự động trong giai đoạn 2025-2026.
   • Chủ thể thực hiện: Các trường đại học, viện nghiên cứu và doanh nghiệp sản xuất thiết bị y tế.

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

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

• Bác sĩ chẩn đoán hình ảnh và kỹ thuật viên X-quang: Giúp hiểu rõ bản chất vật lý của sự tạo ảnh và tác động tiêu cực của bức xạ tán xạ, từ đó lựa chọn tối ưu điện áp phát tia và vị trí đặt lưới chống tán xạ nhằm thu được hình ảnh chẩn đoán sắc nét nhất với liều chụp thấp nhất.
• Kỹ sư Vật lý Y sinh và cán bộ an toàn bức xạ: Cung cấp cơ sở định lượng về liều hấp thụ, Kerma và liều lượng hiệu dụng theo tiêu chuẩn ICRP để thiết kế hệ thống che chắn phòng chụp và kiểm chuẩn chất lượng thiết bị định kỳ.
• Giảng viên, học viên cao học và nghiên cứu sinh ngành Vật lý: Cung cấp tài liệu tham khảo chi tiết về thuật toán mô phỏng Monte Carlo, kỹ thuật lấy mẫu biến ngẫu nhiên và mã nguồn xử lý tương tác bức xạ photon năng lượng thấp.
• Các nhà phát triển và sản xuất thiết bị y tế: Ứng dụng các thông số tính toán về phin lọc chùm tia, hiệu suất chuyển đổi của màn huỳnh quang đất hiếm và lưới chống tán xạ để tối ưu hóa thiết kế phần cứng của máy chụp X-quang thế hệ mới.

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

Tán xạ Rayleigh có gây ảnh hưởng đáng kể đến chất lượng ảnh X-quang không?
Dù chỉ chiếm từ 5% đến 12% tổng số tương tác ở dải năng lượng từ 30 keV đến 70 keV, tán xạ Rayleigh làm đổi hướng bay của photon mà không làm suy giảm năng lượng. Hiện tượng này khiến photon đến sai vị trí hình học trên đầu dò, làm mờ đường biên giải phẫu trong các kỹ thuật đòi hỏi độ phân giải cao như chụp nhũ ảnh.

Tại sao tán xạ Compton lại là nguyên nhân chính làm giảm độ tương phản hình ảnh?
Tán xạ Compton chiếm ưu thế tuyệt đối ở năng lượng trên 26 keV trong mô mềm. Tương tác này tạo ra các photon tán xạ ngẫu nhiên bay lệch quỹ đạo thẳng ban đầu. Khi tỷ số tán xạ đạt mức 3, có tới 75% photon đến đầu dò chỉ tạo thành phông nền mờ, triệt tiêu độ tương phản giữa các mô lân cận.

Phương pháp Monte Carlo có ưu thế gì so với đo đạc thực nghiệm trên mô hình vật lý?
Phương pháp Monte Carlo cho phép theo dõi chi tiết từng bước tương tác của hàng triệu photon trong cấu trúc giải phẫu 3D phức tạp mà không phát sinh liều nhiễm xạ thực tế. Phương pháp này giúp phân tách chính xác tỷ lệ đóng góp của từng cơ chế tương tác riêng biệt mà thực nghiệm không thể đo đạc độc lập.

Tại sao cần sử dụng tấm lọc nhôm để làm cứng chùm tia X-quang?
Chùm tia X phát ra từ bóng phát là chùm đa sắc chứa nhiều photon năng lượng thấp dưới 15 keV. Các photon mềm này sẽ bị hấp thụ hoàn toàn ở lớp da và mô nông của bệnh nhân mà không thể xuyên qua tới phim để tạo ảnh. Tấm lọc nhôm dày 1.5 mm đến 2.5 mm giúp loại bỏ các photon này, giảm liều chiếu trên da hiệu quả.

Màn khuếch đại phosphor đất hiếm mang lại lợi ích gì so với màn Canxi Vonfamat?
Màn phosphor đất hiếm như Gadolinium Oxysulfide có hiệu suất chuyển đổi năng lượng tia X thành ánh sáng khả kiến đạt 15%, cao gấp 3 lần so với mức 5% của Canxi Vonfamat. Nhờ đó, hệ thống chỉ cần liều chiếu tia X thấp hơn khoảng 60% mà vẫn đạt được độ đen phim tiêu chuẩn cho chẩn đoán.

Kết luận

• Luận văn đã hệ thống hóa toàn diện cơ sở vật lý tương tác của chùm photon năng lượng thấp từ 10 keV đến 300 keV với môi trường vật chất sinh học.
• Xây dựng thành công chương trình mô phỏng Monte Carlo theo dõi chính xác từng lịch sử hạt photon qua các môi trường mô mềm và xương.
• Định lượng rõ nét vai trò gây suy giảm độ tương phản của tán xạ Compton và hiện tượng mờ biên giải phẫu do tán xạ kết hợp Rayleigh.
• Khẳng định tính ưu việt của màn tăng sáng phosphor đất hiếm với hiệu suất chuyển đổi 15% trong việc cắt giảm liều bức xạ cho người bệnh.
• Đề xuất các giải pháp kỹ thuật khả thi về điện áp phát tia, bộ lọc làm cứng chùm tia và lưới chống tán xạ phục vụ tối ưu hóa chẩn đoán y khoa.

Đóng góp cốt lõi của công trình là đặt nền tảng tính toán số chính xác cho chuyên ngành Vật lý Y học tại Việt Nam, mở ra hướng ứng dụng tin học hóa trong mô phỏng tiền lâm sàng. Trong giai đoạn 6 đến 12 tháng tới, hướng phát triển tự nhiên của đề tài là mở rộng mô hình Monte Carlo trên cấu trúc cắt lớp vi tính 3D và tích hợp thuật toán xử lý ảnh số thông minh. Các chuyên gia y tế, kỹ sư vật lý y sinh và nhà nghiên cứu được khuyến khích khai thác ngay các dữ liệu và giải pháp từ công trình này để nâng cao chất lượng chẩn đoán và an toàn bức xạ trong thực tế.

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

Chương 1: TƯƠNG TÁC CỦA TIA X VỚI VẬT CHẤT Tia X nói chung là một phẩn của bức xạ gamma, bức xạ điện tử năng lượng cao với bước sóng A tương ứng với năng lượng trong khoảng (10 keV — 300 keV). Khi xem xét tương tác của tia gamma đối với vật chất, ta vận dụng thuyết photon của Einstein, xem tia gamma là tập hợp của các photon có năng hv lượng E=hv và động lượng =: Khi đi qua môi trường vật chất, các photon tia gamma có thể truyền qua, tán xạ hoặc bị hấp thụ. Có bốn loại tương tác chính của photon tia gamma với vật chất: (a) Tán xạ Rayleigh, (b) Tán xạ Compton, (c) Hiệu ứng quang điện, (d) Phản ứng tạo cặp. Trong đó có ba loại đầu tiên đóng vai trò quan trọng trong tương tác tia X trong ngành X - quang chẩn đoán và y học hạt nhân.

Quá trình suy giảm năng lượng của tia X khi qua môi trường vật chất được thể hiện cụ thể thông qua các tương tác sau: 1.1 Tán xạ Rayleigh: SVTH : Phùng Nhật Anh 4 Niên khoá: 2000 - 2004 LUẬN VAN TỐT NGHIỆP GVHD : TS. Đỗ Xuân Hội - TS. Huynh Quang Linh Trong tán xa Rayleigh, còn gọi là tán xa kết hợp (coherent), các photon tới tương tác và kích thích với toàn bộ nguyên tử, còn trong tan xạ Compton và hiệu ứng quang điện, photon tương tác và kích thích với từng clectron riêng lẻ. Sự tương tác này xảy ra chủ yếu ở những tia X dùng trong chẩn đoán có năng lượng thấp, như sử dung trong lĩnh vực Chup tia X vứ (15-30 keV) (mammography).

Trong quá trình tán xa Rayleigh, năng lượng cu@ photon tới truyền cho toàn bộ các electron trong nguyên tử bia dao động đồng bộ và bức xạ năng lượng, phat ra một photon có cùng năng lượng với photon tới nhưng theo một hướng khác. Ở sự tương tác này, electron không bị bứt ra ngoài, vì vậy không xảy ra sự ion hoá. Nói chung, góc tấn xạ càng lớn khi năng lượng của tia X giảm. Trong chụp ảnh y khoa, sự tan xạ tia X có ảnh hưởng xấu đến chất lượng ảnh.

Tuy nhiên, loại tương tác này có xác suất xảy ra thấp trong khoảng năng lượng chẩn đoán. Trong các mô mềm tương tác này chiếm ít hơn 5% tương tác của tia X có năng lượng trên 70 keV và chỉ chiếm khoảng 12% của tương tác khi photon ở mức năng lượng 30 keV. Tiết diện vi phân của tán xa Rayleigh: dơ 1 2 2 2 piled 5 Tol z + cos* 9)|FCX,Z)|ˆ I-l (1-1) r, là bán kính cổ điển của electron, r„= 2,8. F(x,Z) là thừa số hình thức (form factor), F(X,Z) hay F(hv,8,Z) phụ thuộc vào môi trường (Z) và năng lượng của photon tới thông qua thông số X.

9 Sin= X= — (1-2) Trong đó F(x,Z) là thừa số dang: * sin[(2/^)rs] F(X,Z) „2)= Jota ay + (1-3) 1-3 Với: - p(r) là mật độ electron. - _ r là khoảng cách từ electron đến hạt nhân. Ta có: dQ = sin9.dọ (1-4) ra: ta suy ) Từ(1-1 SVTH : Phùng Nhật Anh $ Niên khoá: 2000 - 2004 LUẬN VĂN TỐT NGHIỆP GVHD : TS. Đỗ Xuân Hội - TS.

Huynh Quang Linh =; - ú +cos 7 @)|F(X, zy 2nsin® (1-5) Công thức (1-5) có thé có được như sau: Trong đó “ene là hàm phân bố Thompson theo góc 0 và có giá tri: Tư1 r2(1 + cos2 6)2xsin6 (1-7) Tỷ số ma đóng vai trò là hàm phân phối f(6). Ta dùng công thức nay để mô phỏng tìm góc tán xạ 6 của photon ứng với từng giá trị của năng lượng photon tới hv (hay ứng với thông số X).Z)F sin8dôde (1-8) 00 g= mủ Jaxes 6)|F(x,ya i 2nd(cos@) (1-9) Tán xạ Rayleigh mặc dù có xác suất thấp nhưng góp phần đáng kể đến sự nhoè của ảnh X - quang. 12 Tán xạ Compton: SVTH : Phùng Nhật Anh 6 Nién khod: 2000 - 2004 LUẬN VĂN TỐT NGHIỆP GVHD : TS. Đỗ Xuân Hội ~ TS.

Huỳnh Quang Linh Tán xạ Compton, còn gọi là tán xạ không kết hợp (incohrent), là tương tác chủ yếu của photon tia X và tia gamma với mô mềm trong khoảng năng lượng chẩn đoán. Thực tế, tán xạ Compton không những chiếm ưu thế trong khoảng năng lượng chẩn đoán trên 26 keV trong các mô mềm mà còn chiếm ưu thế ở mức năng lượng chẩn đoán xấp xỉ 30 MeV. Sự tương tác này hầu như xảy ra giữa photon và các electron lớp ngoài cùng. Các electron thoát ra khỏi nguyên tử và photon bị tán xạ, nó bị mất bớt năng lượng.

Như tất cả các loại tương tác, cả năng lượng và xung lượng phải bị biến đổi. Vì vậy, năng lượng của photon tới (E) bằng tổng năng lượng của photon tán xạ (E) và động năng của electron thoát ra (E. Thể hiện trong công thức: E=E +E. Năng lượng liên kết của electron bứt ra là tương đối nhỏ và có thể bỏ qua.

E’ =hv’ là năng lượng của photon tan xạ. E=hvw là năng lượng của photon tới. Năng lượng của photon tán xạ được tính theo năng lượng của photon tới và góc tán xạ 6: 1 hv'=hv——— ar, 6) -10 (1-10) œ(1 — cos8) E. 1-11 ) a = hv/m,cŸ sE / 511keV là tỷ số năng lượng tia X và năng lượng nghỉ của electron.

Khi năng lượng của photon tới giảm, photon tấn xạ và electron tấn xạ hướng gần nhau hơn (hình 1. Trong chụp ảnh tia X, các photon được phát hiện bởi một máy thu. Vì vậy, giảm sự tương phản của ảnh. Nếu như góc tán xạ cố định, khi tăng năng lượng của photon tới thì phần năng lượng truyền photon tấn xạ giảm.

Do đó, đối với những photon tới có năng lượng cao, phan lớn năng lượng được truyền cho electron tấn xạ. Vi dụ, với một góc tán xạ 60°, nếu năng lượng của photon tới là 100 keV thì năng lượng của photon tán xạ chiếm 90%, nhưng nếu năng lượng của photon tới là 5 MeV thì phần năng lượng truyền cho photon tán xạ chỉ là 17%. Khi tán xa Compton xảy ra đối với tia X có năng lượng thấp, được ứng dụng trong chụp ảnh chẩn đoán (tit 18 đến 150 keV), phần lớn năng lượng của SVTH : Phùng Nhật Anh 7 Niên khoá: 2000 - 2004 LUẬN VĂN TỐT NGHIỆP GVHO : TS. Đỗ Xuân Hội - TS.

Huỳnh Quang Linh photon tới được truyền cho photon tán xạ. Ví dụ, tương tác Compton của một photon có năng lượng 80 keV, năng lượng nhỏ nhất của photon tán xạ là 61 keV. Vì vậy với sự mất năng lượng lớn nhất, photon tán xạ vẫn có năng lượng tương đối cao và có khả năng xuyên qua các mô. Ta cũng nhận thấy rằng năng lượng của tia X sau tấn xạ và động năng của electron đều phụ thuộc vào góc tán xạ.

Năng lượng của tia X thay đổi từ 0 đến hv và động năng của electron thay đổi từ E „. - 8 = 0": Năng lượng tia X không đổi, động năng của electron bằng 0. (1-14) l+a = sh = (1-15) em 1420 - 8= 180": Photon tia X tương tác trực tiếp với electron, electron bay theo hướng photon tới, photon tới bị tán xạ ngược lại một góc 180”. Sau tán xạ, on tia X có phot on nă có động năng lượng hv„„ va electr ng cực đại.

l BV nin =hv (1-16) I+2œ 1 KT” (1-17) Tiết diện vi phân của tán xạ Compton trên một electron được tính theo công thức O. Nishina [5]: (1-18) 2 ' số =4 tr) (+t sin*o) dQ 2 “(hv/ (hv'` hv 2 t, = ——— = 2,82.10°cm là bán kính cổ điển của electron.c Thay hv từ phương trình (2-10) vào phương trình(2-18) ta được: d seo. a (1—cos6) | (1-19) 2 2 đợ lại 1 SVTH : Phong Nhật Anh 8 Niên khoá: 2000 - 2004 LUẬN VAN TỐT NGHIỆP GVHO : TS. Đỗ Xuân Hội - TS.

Huỳnh Quang Linh Tỷ số se đóng vai trò như hàm phân phối g(6), hay hàm mật độ xác suất theo góc tan xạ 9. Ứng với mỗi giá trị năng lượng xác định, từ hàm phân phối xác suất g(6) ta mô phỏng để tìm giá trị góc tan xa 8. x10" Mat do zac sust Compton 7 ?* s~ \ \ \ > s- \ 1 \\ | - \ + 3- eet x -— — | = , ° os 1 us 2 25 3 3s Theta Hình 1-3: Mật độ xác suất Compton theo góc 9 ứng với năng lượng 200 keV. Đối với photon có năng lượng thấp, a<<1, phương trình trở thành: d.c Ị 2 — 2/0(1+cos 6) ti oes (1-20) - Tích phân phương trình (1-19) theo toàn bộ góc khối Q ta được tiết diện tán xạ Compton toàn phan: ° | a? l+2a a ad Ì*x — +%.| a l+œ|2(l+œ) =2 a 1 | ~—In(1+2a) | +—In(1 + 2a)- 1+3a 1-21 Phương trình này chỉ là tiết diện tấn xạ toàn phẩn trên một electron.

Dé có tiết diện tán xạ toàn phẩn trên một nguyên tử ta phải nhân tiết diện toàn phần của electron với số Z của nguyên tử. SVTH : Phùng Nhật Anh ọ Niện khoá: 2000 - 2004 LUẬN VĂN TỐT NGHIỆP GVHD : TS. Đỗ Xuân Hội - TS. Huỳnh Quang Linh 100 so = wo 7° 60 s0 40 i » 2 10}, °Q 30 S6 0 lãi 150 180 Scatter Angle (degrees) Hình 1-4: Dé thị minh hoạ xác Suất tương đối của tán xạ Compton theo góc tán xạ của photon có năng lượng 20, 80, 140 keV trong mô.

ata Năng lượng của photon tới phải cao hơn năng lượng liên kết ci#thi hiệu ứng Compton mới có thể xảy ra. Khi năng lượng của photon tới tăng, xác suất tán xạ Compton tăng so với tán xạ Rayleigh và hấp thụ quang điện. Xác suất tán xạ phụ thuộc vào mật độ electron (số electron/g*mat độ). Vì vậy xác suất của tán xạ Compton gần như chỉ phụ thuộc vào Z và xác suất tán xạ Compton trên đơn vị thể tích tỷ lệ với mật độ của vật liệu.

So sánh với những nguyên tố khác, do không có nơtron trong nguyên tử hiđro dẫn đến tăng gần như gấp đôi mật độ electron. Vì vậy, những vật liệu có hiđro có xác suất tán xạ Compton cao hơn so với vật liệu không có hiđro có cùng khối lượng. Tiết diện tấn xạ Compton, động năng của electron thay đổi theo năng lượng của photon tia X tới và số Z của nguyên tử vật chất. Tán xạ Compton là tác nhân chủ yếu gây nên sự nhoè của ảnh X - quang.3 Hiệu ứng quang điện: A: 0.5: Hiệu ứng quang điện SVTH : Phùng Nhật Anh 10 Nién khod: 2000 - 2004 LUAN VĂN TỐT NGHIỆP GVHD : TS.

Đỗ Xuân Hội - TS. Huỳnh Quang Linh Trong hiệu ứng quang điện, toàn bộ năng lượng của photon tới được truyền cho electron và electron này bứt ra khỏi nguyên tử.