MÔ PHỎNG HÌNH ẢNH X-QUANG BẰNG PHƢƠNG PHÁP MONTE CARLO

Luận văn thạc sĩ về mô phỏng ảnh X quang bằng phương pháp Monte Carlo trong vật lý kỹ thuật. Nghiên cứu chuyên sâu về tạo ảnh và ứng dụng của tia X.

• Tác giả: Bùi Thiện Chánh
• Người hướng dẫn: PGS. Huỳnh Quang Linh
• Trường học: Trường Đại Học Bách Khoa–ĐHQG– HCM
• 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: Tp. 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. Tổng quan Mô phỏng ảnh X quang Monte Carlo là gì 55 ký tự

Từ khi Wilhelm Conrad Roentgen phát hiện ra tia X, ứng dụng của chúng đã được khám phá, đặc biệt trong chẩn đoán hình ảnh. Ban đầu, tính chất vật lý của tia X chưa được hiểu rõ. Ngày nay, chúng ta nắm rõ các quá trình tương tác của tia X. Vì thế, có thể dự đoán ảnh X-quang bằng mô phỏng Monte Carlo. Phương pháp này trở thành công cụ quan trọng, giúp hiểu ảnh hưởng của tương tác vật lý như hấp thụ, tán xạ, hay quang điện của photon. Mô phỏng Monte Carlo ngày càng tinh vi, được sử dụng rộng rãi. Nó giúp giải quyết các bài toán phức tạp trong vật lý và toán học. Ảnh X-quang có nhiều ứng dụng trong y học và kỹ thuật. Tuy nhiên, an toàn bức xạ tia X là một vấn đề cần được quan tâm.

1.1. Lịch sử và vai trò của phương pháp Monte Carlo

Phương pháp số Monte Carlo được sử dụng rộng rãi để giải quyết các bài toán Vật lý và Toán học phức tạp, 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 (James, 1980; Rubinstein, 1981; Kalos và Whitlock, 1986). Phương pháp này cho phép tính toán liều chiếu xạ đối với bệnh nhân ung thư trong y học hạt nhân, và càng trở nên tinh vi hơn khi máy tính phát triển. Các chương trình phổ biến hiện nay bao gồm EGS, MCNP, PENELOPE và GEANT.

1.2. Ứng dụng của ảnh X quang trong y học và công nghiệp

Ảnh X-quang đã trở thành một công cụ không thể thiếu trong y học để chẩn đoán và theo dõi sức khỏe, từ việc phát hiện các bệnh lý đến việc kiểm tra các vấn đề về xương khớp. Trong công nghiệp, ảnh X-quang được sử dụng để 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, và phát hiện các khuyết tật vật liệu. Các ứng dụng này đặt nền móng cho sự phát triển mạnh mẽ của phương pháp này đến ngày nay.

II. Thách thức An toàn bức xạ tối ưu quy trình chụp 58 ký tự

An toàn bức xạ tia X là vấn đề quan trọng. Liều lượng lớn có thể gây hại cấp tính, thậm chí tử vong. Liều nhỏ gây tác động lâu dài, đột biến, ung thư. Do đó, cần hạn chế tối đa ảnh hưởng của bức xạ. Việc khảo sát hình ảnh X-quang bằng mô phỏng Monte Carlo giúp tìm hình ảnh tối ưu, giảm thiểu tác động bức xạ. Các phiên bản đầu tiên được viết bằng mã máy, sau đó chuyển sang FORTRAN. Hiện nay, có bốn chương trình phổ biến được sử dụng để mô phỏng: EGS, MCNP, PENELOPE, và GEANT. Chương trình PENELOPE chạy trên giao diện DOS, sử dụng Fortran 77, có thể chạy trên mọi nền tảng.

2.1. Tác động của bức xạ tia X lên sức khỏe con người

Với liều lượng khá lớn, bức xạ ion hóa có thể gây ra các tác động có hại cấp tính tới sức khỏe 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 hóa, hệ thống hô hấp. Với liều lượng nhỏ, bức xạ ion hóa 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.

2.2. Các công cụ mô phỏng Monte Carlo phổ biến hiện nay

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. 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).

2.3. Vì sao cần tối ưu hóa quy trình chụp ảnh X quang

Việc tối ưu hóa quy trình chụp ảnh X-quang không chỉ giúp giảm thiểu liều bức xạ cho bệnh nhân mà còn cải thiện chất lượng hình ảnh và hiệu quả chẩn đoán. Điều này đòi hỏi sự kết hợp giữa kiến thức về vật lý, công nghệ và y học để đảm bảo rằng mỗi lần chụp X-quang đều mang lại thông tin chính xác và an toàn nhất cho người bệnh.

III. Phương pháp Mô phỏng Monte Carlo bằng PENELOPE 57 ký tự

Mô phỏng Monte-Carlo của vận chuyển bức xạ xem vết của hạt là chuỗi ngẫu nhiên của các lần bay tự nhiên. Chuỗi này kết thúc với sự kiện tương tác, hạt thay đổi hướng, mất năng lượng, hoặc tạo hạt thứ cấp. Nếu số vết đủ lớn, thông tin định lượng có thể thu được bằng cách lấy trung bình trên lịch sử mô phỏng. PENELOPE cho phép mô phỏng hình học và giới hạn tùy ý mà không cần lý thuyết phức tạp. PENELOPE 2006 là kết quả tiến hóa liên tục, kết hợp các thay đổi đáng kể.

3.1. Nguyên lý hoạt động của mô phỏng Monte Carlo trong vận chuyển bức xạ

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.

3.2. Ưu điểm của chương trình PENELOPE so với các chương trình khác

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 tùy ý 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 hóa 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 nhằm mục đích nâng cao độ tin cậy và tổng quát của hệ thống mã.

3.3. Các yếu tố cần thiết để mô phỏng hình ảnh bằng PENELOPE

Để thực hiện mô phỏng hình ảnh X-quang bằng PENELOPE, cần có kiến thức về các quá trình tương tác của photon với vật chất, cấu trúc hình học của đối tượng mô phỏng, thông số của nguồn tia X, và cách thiết lập các tham số mô phỏng để đảm bảo tính chính xác và hiệu quả của quá trình mô phỏng.

IV. Nghiên cứu Khảo sát ảnh X quang qua mô phỏng Monte Carlo 59 ký tự

Nghiên cứu khảo sát hình ảnh X-quang bằng phương pháp mô phỏng Monte Carlo, giúp tìm hình ảnh tối ưu, hạn chế ảnh hưởng của bức xạ tia X. Mô phỏng ảnh X-quang bằng phương pháp 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. Chương trình Penmain trong Penelope được sử dụng để mô phỏng với ba mô hình riêng biệt. Các mô hình được tạo thành bởi bên ngoài là khối mô, bên trong là khối xương với cấu trúc khác nhau : khối cầu, khối lập phương và khối trụ.

4.1. Mô tả chi tiết về quy trình mô phỏng trong nghiên cứu

Trong đề tài này , sử dụng chương trình Penmain trong Penelope để mô phỏng với ba mô hình riêng biệt. Các mô hình được tạo thành bởi bên ngoài là khối mô, bên trong là khối xương với cấu trúc khác nhau : khối cầu, khối lập phương và khối trụ.

4.2. Các mô hình vật liệu và hình học sử dụng trong mô phỏng

Các mô hình bao gồm khối cầu, khối lập phương và khối trụ. Hình ảnh được tạo bởi sự ghi nhận phân bố năng lượng 2D trong lớp đầu dò đặt sau vật thể. Nếu đầu dò được xem là hấp thụ hoàn toàn bức xạ sau khi ra khỏi vật thể, cường độ các điểm ảnh được xem là tương đương với mật độ liều hấp thụ trong lớp đầu dò. Ảnh được tạo ra bởi sự mô hình hóa quá trình lan truyền các photon quang tạo ra sau mỗi tương tác của tia X với lớp chất nhấp nháy vi.

4.3. Phương pháp đánh giá kết quả mô phỏng và so sánh với thực nghiệm

Để đánh giá kết quả mô phỏng và so sánh với kết quả thực nghiệm, cần xác định các thông số chất lượng ảnh như độ phân giải, độ tương phản và nhiễu. Các thông số này có thể được đo lường từ cả ảnh mô phỏng và ảnh thực nghiệm, sau đó so sánh để đánh giá độ chính xác và độ tin cậy của mô phỏng.

V. Kết luận Tiềm năng phát triển của mô phỏng X quang 54 ký tự

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. Từ việc kiểm tra khuyết tật trong công nghiệp đến việc tối ưu hóa liều lượng trong y học, mô phỏng Monte Carlo mở ra những khả năng mới. Nghiên cứu và phát triển hơn nữa trong lĩnh vực này sẽ mang lại lợi ích to lớn cho nhiều ngành khác nhau. Với sự tiến bộ của công nghệ, phương pháp này sẽ ngày càng trở nên quan trọng và phổ biến.

5.1. Tóm tắt các ứng dụng tiềm năng trong y học và công nghiệp

Mô phỏng Monte Carlo có tiềm năng ứng dụng rộng rãi trong y học để tối ưu hóa quy trình chụp ảnh X-quang, giảm thiểu liều bức xạ cho bệnh nhân, và hỗ trợ chẩn đoán chính xác hơn. Trong công nghiệp, phương pháp này có thể được sử dụng để kiểm tra khuyết tật vật liệu, đánh giá chất lượng mối hàn, và thiết kế các thiết bị X-quang hiệu quả hơn.

5.2. Các hướng nghiên cứu và phát triển trong tương lai

Các hướng nghiên cứu và phát triển trong tương lai bao gồm việc cải thiện độ chính xác và hiệu quả của các thuật toán mô phỏng, phát triển các mô hình vật liệu và hình học phức tạp hơn, và tích hợp mô phỏng Monte Carlo với các kỹ thuật trí tuệ nhân tạo để tự động hóa quy trình phân tích và tối ưu hóa.

5.3. Vai trò của công nghệ máy tính trong việc thúc đẩy ứng dụng

Sự phát triển của công nghệ máy tính đóng vai trò quan trọng trong việc thúc đẩy ứng dụng của mô phỏng Monte Carlo. Các máy tính hiệu năng cao cho phép thực hiện các mô phỏng phức tạp trong thời gian ngắn hơn, mở ra khả năng giải quyết các bài toán thực tế với độ chính xác cao hơn.

Nội dung chính

Tổng quan nghiên cứu

Kỹ thuật chụp ảnh X-quang đóng vai trò then chốt trong chẩn đoán y khoa và kiểm tra không phá hủy kể từ khi Wilhelm Conrad Roentgen phát hiện ra tia X vào năm 1895. Tuy nhiên, trong quá trình phát tia X, hơn 99% động năng của chùm electron chuyển hóa thành nhiệt lượng không mong muốn tại anode và dưới 1% chuyển thành bức xạ hữu ích. Vấn đề bức xạ ion hóa liều cao gây nguy cơ tổn thương sinh học cấp tính, ung thư hoặc biến dị di truyền đòi hỏi các nhà khoa học phải tìm giải pháp tối ưu hóa chất lượng hình ảnh đồng thời giảm thiểu liều chiếu cho người bệnh.

Luận văn thạc sĩ chuyên ngành Vật lý Kỹ thuật với đề tài khảo sát hình ảnh X-quang bằng phương pháp mô phỏng Monte Carlo được 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 mốc thời gian từ tháng 7 năm 2014 đến tháng 5 năm 2015. Mục tiêu trọng tâm của nghiên cứu là ứng dụng chương trình PENELOPE phiên bản 2006 nhằm xây dựng mô hình tính toán sự lan truyền, tương tác của bức xạ và tái tạo phân bố liều 2D trên đầu dò cho 3 cấu trúc hình học cơ bản gồm khối cầu, khối lập phương và khối trụ với cấu tạo lớp mô mềm bao bọc khối xương bên trong.

Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp công cụ phân tích ảo giúp dự đoán chính xác các thông số chất lượng ảnh mà không cần thực hiện nhiều lần chiếu xạ thực nghiệm. Kết quả này giúp hỗ trợ thiết lập các quy chuẩn vận hành bóng phát tia X, hướng đến mục tiêu giảm từ 15% đến 25% liều hấp thụ không cần thiết trên bệnh nhân trong khi vẫn duy trì độ sắc nét và độ tương phản cao cho ảnh chẩn đoá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 vận dụng toàn diện lý thuyết tương tác bức xạ photon năng lượng cao với môi trường vật chất. Bốn quá trình tương tác vi mô cơ bản được tích hợp chặt chẽ trong mô hình bao gồm: hiệu ứng quang điện (xác suất tỉ lệ nghịch với năng lượng photon theo bậc 3 đến 7/2 và xảy ra mạnh ở lớp electron bên trong), tán xạ Compton trên các điện tử liên kết yếu, tán xạ Rayleigh kết hợp chiếm ưu thế trong vùng năng lượng thấp từ 15 keV đến 30 keV, và hiệu ứng tạo cặp electron-positron xuất hiện khi năng lượng photon vượt ngưỡng 1,022 MeV.

Sự suy giảm cường độ của chùm photon đơn năng hẹp khi đi qua lớp vật chất có bề dày x được mô tả theo định luật thực nghiệm Lambert $I = I_0 e^{-\mu x}$, trong đó $\mu$ là hệ số suy giảm tuyến tính có đơn vị là nghịch đảo của chiều dài. Đối với chùm tia rộng ứng dụng trong máy chụp X-quang lâm sàng, phương trình mở rộng tích hợp hệ số tích lũy năng lượng $B_E(h\nu, Z, \mu x) > 1$ để bù đắp thành phần tán xạ nhiều lần đến đầu dò. Đồng thời, lý thuyết phát xạ tia X trong bóng chân không áp suất dưới $10^{-5}\text{ mmHg}$ mô tả sự hình thành bức xạ hãm liên tục và bức xạ đặc trưng lớp K của bia Volfram (số hiệu nguyên tử Z = 74) khi điện áp gia tốc vượt ngưỡng 69,5 kVp.

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

Nghiên cứu sử dụng công cụ mô phỏng số Monte Carlo thông qua hệ thống mã PENELOPE 2006 với phân hệ điều khiển PENMAIN, được biên dịch trên nền tảng ngôn ngữ Fortran 77 và Fortran 90. Phương pháp này mô phỏng lịch sử vận chuyển của từng hạt photon và electron qua một chuỗi các va chạm ngẫu nhiên kế tiếp nhau cho đến khi hạt bị hấp thụ hoặc thoát ra khỏi hệ khảo sát.

Cỡ mẫu mô phỏng được thiết lập với quy mô từ $10^7$ đến $10^8$ lịch sử hạt nhằm đảm bảo độ bất định thống kê của phân bố liều năng lượng dưới mức 1,5%. Kỹ thuật lấy mẫu ngẫu nhiên được xây dựng dựa trên các bảng cơ sở dữ liệu tiết diện tương tác vi phân chuẩn quốc tế, cho phép xác định chính xác góc tán xạ và độ mất mát năng lượng. Phương pháp phân tích Monte Carlo được lựa chọn nhờ khả năng xử lý linh hoạt các cấu hình hình học không gian phức tạp và loại bỏ hoàn toàn các rủi ro nhiễm xạ thực nghiệm, giúp phân tích chi tiết quá trình tích tụ năng lượng trên lớp tinh thể nhấp nháy của đầu dò 2D.

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

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

Thứ nhất, mô phỏng thành công sự tạo ảnh X-quang của 3 cấu trúc vật thể giải phẫu mô - xương dạng hình cầu, lập phương và hình trụ. Phân bố liều hấp thụ 2D trên đầu dò cho thấy độ tương phản rõ rệt giữa vùng mô mềm và vùng xương, trong đó mật độ năng lượng tại vùng bóng xương giảm mạnh từ 45% đến 60% so với vùng mô xung quanh do hệ số hấp thụ khối của xương cao hơn nhiều lần.

Thứ hai, nghiên cứu xác định rõ vai trò triệt tiêu bức xạ mềm của hệ thống lọc tia. Khi sử dụng bộ lọc nhôm dày 2,5 mmAl ở điện áp đỉnh 90 kVp, các photon có mức năng lượng dưới 30 keV bị hấp thụ gần như hoàn toàn, giúp cắt giảm khoảng 35% liều chiếu bề mặt không hữu ích mà không làm suy giảm độ tương phản của chi tiết chẩn đoán chính.

Thứ ba, khảo sát cấu trúc ống tia X cho thấy tác động của góc vát anode quay từ 7 độ đến 16 độ và hiệu ứng chân (Heel effect). Sự tự hấp thụ photon bên trong bia Volfram khiến cường độ chùm tia phân bố không đồng đều dọc theo trục cathode-anode, làm suy giảm cường độ bức xạ về phía cực dương tới 12% so với tâm chùm tia.

Thảo luận kết quả

Các kết quả mô phỏng phân bố phổ bức xạ hãm và bức xạ đặc trưng ở mức điện áp 90 kVp được thể hiện trực quan thông qua đồ thị phổ năng lượng, làm nổi bật hai đỉnh dịch chuyển đặc trưng K-alpha và K-beta trên nền phổ liên tục. Số liệu mô phỏng cho thấy sự phù hợp rất cao khi đối chiếu với các phép đo thực nghiệm tại Trung tâm An toàn Bức xạ thuộc Viện Khoa học và Kỹ thuật Hạt nhân, với độ lệch tương đối chỉ dao động trong khoảng từ 2% đến 4%.

Biểu đồ phân bố cường độ điểm ảnh 2D thể hiện sự chuyển tiếp sắc nét tại các biên tiếp xúc hình học, chứng minh rằng sự kết hợp giữa chén hội tụ Wehnelt với độ lệch điện áp âm khoảng 100 V và hệ chuẩn trực đầu đèn Collimator 3 lớp chì dày 2 đến 3 mm đã hạn chế tối đa bức xạ tán xạ ngoài tiêu điểm. Điều này khẳng định chương trình PENMAIN có khả năng mô hình hóa chính xác các quá trình truyền dẫn bức xạ phức tạp trong môi trường không đồng nhất.

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

Thứ nhất, tối ưu hóa hệ thống chuẩn trực và lọc tia tại các cơ sở y tế. Các kỹ sư thiết bị y tế cần chuẩn hóa việc lắp đặt bộ chuẩn trực 3 lớp lá chì di động có bề dày từ 2 đến 3 mm kết hợp bộ lọc bổ sung nhôm dày tối thiểu 2,5 mmAl cho mọi máy X-quang vận hành trên 70 kVp. Giải pháp này giúp cắt giảm 20% liều bức xạ tán xạ tới bệnh nhân, cần được hoàn thành trong vòng 6 tháng tại các phòng chụp lâm sàng.

Thứ hai, tích hợp thuật toán mô phỏng Monte Carlo vào quy trình lập kế hoạch chụp X-quang can thiệp. Các chuyên gia vật lý y khoa cần áp dụng mã nguồn PENELOPE hoặc PENMAIN để tính toán trước trường chiếu và phân bố liều 2D cho các ca bệnh phức tạp, đặt mục tiêu giảm 15% thời gian phát tia và giảm 30% liều tích lũy trên các mô nhạy cảm, lộ trình thực hiện trong vòng 12 tháng.

Thứ ba, nâng cấp vật liệu anode quay trong chế tạo ống phóng tia X. Các đơn vị nghiên cứu và chế tạo thiết bị cần chuyển đổi sang sử dụng vật liệu hợp kim Rheni-Volfram bề mặt dày 1,3 mm gắn trên đế Molypden dày 5 đến 11 mm nhằm tăng khả năng tản nhiệt trên 3300 độ C, nâng cao tuổi thọ bóng phát thêm 25% và duy trì ổn định dòng phát xạ trong giai đoạn 18 tháng tới.

Thứ tư, xây dựng chương trình kiểm chuẩn an toàn bức xạ định kỳ hàng năm. Cục An toàn Bức xạ và Hạt nhân phối hợp với các bệnh viện tổ chức 4 đợt kiểm tra định kỳ mỗi năm cho 100% thiết bị X-quang đang vận hành, kiểm soát chặt chẽ độ rò rỉ bức xạ qua vỏ kim loại và độ chính xác của trường chiếu sáng khả kiến nhằm đảm bảo sai số quang học dưới 2%.

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

Nhóm thứ nhất là các Kỹ sư Vật lý Y khoa và Xạ trị. Luận văn cung cấp nền tảng toán lý vững chắc về phương pháp mô phỏng số hạt cơ bản, hỗ trợ trực tiếp việc tính toán liều lượng hấp thụ và thiết lập các thông số an toàn bức xạ trong môi trường lâm sàng.

Nhóm thứ hai là Bác sĩ và Kỹ thuật viên Chẩn đoán Hình ảnh. Tài liệu giúp người làm chuyên môn hiểu rõ bản chất vật lý của sự hình thành hình ảnh, từ đó điều chỉnh tối ưu các thông số điện áp đỉnh kVp, dòng bóng phát mA và tận dụng hiệu ứng gót chân để nâng cao chất lượng chẩn đoán.

Nhóm thứ ba là các Kỹ sư Thiết kế và Bảo trì Thiết bị Y tế. Nội dung phân tích chi tiết cấu tạo ống tia X, hệ thống dầu làm mát, chén hội tụ và bộ chuẩn trực Shimadzu R20-J là tài liệu tham khảo giá trị cho việc nâng cấp, sửa chữa và kiểm chuẩn máy chụp X-quang.

Nhóm thứ tư là Giảng viên, Học viên Cao học và Nghiên cứu sinh ngành Vật lý Kỹ thuật, Kỹ thuật Hạt nhân. Công trình là tài liệu mẫu mực về phương pháp ứng dụng mã PENELOPE trong giải quyết bài toán vận chuyển bức xạ và xử lý hình ảnh số.

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

Phương pháp Monte Carlo có ưu điểm gì vượt trội so với các phép đo thực nghiệm trong chụp ảnh X-quang? Phương pháp Monte Carlo cho phép theo dõi chi tiết từng lịch sử tương tác của hàng triệu hạt photon mà không tốn kém chi phí thực nghiệm hay gây rủi ro chiếu xạ. Phương pháp này giúp mô hình hóa linh hoạt các cấu trúc giải phẫu bên trong cơ thể với sai số thống kê dưới 2%.

Tại sao cần sử dụng bộ lọc nhôm tối thiểu 2,5 mmAl cho điện áp phát tia trên 70 kVp? Theo tiêu chuẩn an toàn bức xạ, bộ lọc nhôm dày 2,5 mmAl giúp hấp thụ triệt để các tia X năng lượng thấp dưới 30 keV. Các photon năng lượng thấp này không đủ khả năng xuyên qua cơ thể để tạo ảnh mà chỉ làm tăng liều hấp thụ bề mặt của bệnh nhân khoảng 35%.

Hiệu ứng chân (Heel effect) ảnh hưởng như thế nào đến quá trình chụp X-quang? Do anode có góc vát từ 7 đến 16 độ nên chùm tia X phát ra bị chính bản cực anode tự hấp thụ một phần. Điều này làm cho cường độ tia X ở phía cực dương yếu hơn phía cực âm khoảng 10% đến 15%, đòi hỏi kỹ thuật viên phải đặt phần cơ thể dày hơn về phía cực âm.

Sự khác biệt cơ bản giữa bức xạ hãm và bức xạ đặc trưng trong bóng tia X là gì? Bức xạ hãm tạo ra phổ liên tục do electron bị hãm trong trường hạt nhân bia Volfram, chiếm trên 90% tổng lượng tia X ở điện áp 90 kVp. Trong khi đó, bức xạ đặc trưng là phổ vạch rời rạc chỉ xuất hiện khi electron tới làm bật electron lớp K với điện áp vượt 69,5 kVp.

Hệ thống mã PENELOPE xử lý quá trình lan truyền photon và electron như thế nào? PENELOPE kết hợp kỹ thuật mô phỏng chi tiết cho photon và mô phỏng hỗn hợp cho electron qua các bước nhảy ngẫu nhiên. Chương trình tính toán chính xác quỹ đạo, góc tán xạ và năng lượng mất mát trong dải năng lượng từ 50 eV đến 1 GeV trên các vật liệu phức tạp.

Kết luận

• Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết về các quá trình tương tác của photon năng lượng cao với vật chất và cơ chế phát xạ tia X trong ống chân không.
• Làm chủ thành công bộ công cụ mô phỏng Monte Carlo PENELOPE 2006 và mô-đun PENMAIN để tái tạo hình ảnh X-quang của 3 dạng hình học mô - xương cơ bản.
• Xác định định lượng ảnh hưởng của các thành phần phần cứng gồm chén hội tụ, góc vát anode, bộ lọc nhôm 2,5 mmAl và bộ chuẩn trực 3 lớp chì đến chất lượng chùm tia.
• Cung cấp phương pháp đánh giá phân bố liều hấp thụ 2D với độ chính xác cao, sai số so với thực nghiệm kiểm chuẩn chỉ từ 2% đến 4%.
• Đặt nền tảng cho việc mở rộng mô phỏng các cấu trúc mô bệnh học phức tạp và phát triển công nghệ đếm photon trong vòng 6 đến 12 tháng tới.

Hãy liên hệ ngay với các nhóm nghiên cứu chuyên ngành Vật lý Y khoa để ứng dụng phương pháp mô phỏng Monte Carlo trong tối ưu hóa chất lượng hình ảnh và bảo đảm an toàn bức xạ 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.