Tổng quan nghiên cứu
Theo thống kê y tế toàn cầu, ung thư hiện là một trong những nguyên nhân gây tử vong hàng đầu với hàng triệu ca mắc mới mỗi năm, trong đó việc phát hiện sớm qua chẩn đoán hình ảnh giúp tăng tỷ lệ điều trị thành công lên trên 80%. Trong các phương pháp hiện hành, siêu âm nổi bật nhờ ưu điểm không xâm lấn, chi phí thấp và không sử dụng bức xạ ion hóa nguy hại như chụp cắt lớp vi tính (CT) – vốn làm cơ thể hấp thụ liều bức xạ tương đương 6 tháng đến 3 năm trong tự nhiên. Tuy nhiên, kỹ thuật siêu âm phản xạ B-mode truyền thống vẫn tồn tại hạn chế lớn về độ phân giải không gian và không thể phát hiện các vi cấu trúc tổn thương có kích thước nhỏ hơn bước sóng âm.
Kỹ thuật chụp siêu âm cắt lớp (Ultrasound Computed Tomography - UCT) dựa trên tán xạ ngược ra đời nhằm giải quyết bài toán định lượng các thông số vật lý của mô như vận tốc truyền âm (1460 m/s ở mô mỡ đến 1580 m/s ở gan) và độ suy giảm âm. Mặc dù vậy, phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM) đơn tần số thường gặp hiện tượng phân kỳ hoặc mắc kẹt tại cực trị địa phương khi đối tượng có độ tương phản vận tốc âm lớn. Luận văn thạc sĩ chuyên ngành Kỹ thuật Điện tử tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội đã nghiên cứu phát triển kỹ thuật kết hợp hai tần số (Dual-Frequency DBIM - DF-DBIM) nhằm khắc phục triệt để nhược điểm này.
Nghiên cứu tập trung vào việc tối ưu hóa thuật toán tái tạo ảnh trên miền khảo sát kích thước $22 \times 22$ điểm ảnh, đường kính vùng tán xạ $4\lambda$ (tương đương 0,0031 m) với mức chênh lệch vận tốc truyền sóng từ 2% đến 5% dưới tác động của nhiễu Gaussian có tỷ số tín hiệu trên nhiễu (SNR) đạt 26 dB. Kết quả nghiên cứu có ý nghĩa thực tiễn quan trọng khi nâng chỉ số chất lượng ảnh phổ quát ($Q$) lên mức 0,9992 và giảm sai số khôi phục toán học qua 8 bước lặp, tạo tiền đề vững chắc cho việc ứng dụng siêu âm cắt lớp thương mại trong chẩn đoán ung thư sớm.
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 xây dựng trên nền tảng lý thuyết truyền sóng âm trong môi trường không đồng nhất, mô tả thông qua phương trình tích phân Helmholtz kết hợp hàm truyền Green. Khi chùm sóng tới $p^{inc}(r)$ tương tác với mô sinh học, sóng tán xạ $p^{sc}(r)$ được hình thành phụ thuộc vào hàm mục tiêu $O(r)$, trong đó $O(r)$ biểu diễn sự biến thiên số sóng giữa mô bệnh lý và môi trường đồng nhất chuẩn với vận tốc sóng trong nước $c_0 = 1480$ m/s.
Hai lý thuyết tái tạo ảnh phi tuyến giữ vai trò cốt lõi trong nghiên cứu gồm phương pháp lặp Born (BIM) và phương pháp lặp vi phân Born (DBIM). Trong đó, DBIM sử dụng thuật toán tối ưu phi tuyến bậc một để cập nhật liên tục hàm Green theo từng bước lặp, giúp tốc độ hội tụ nhanh hơn đáng kể so với BIM truyền thống. Phương pháp biến phân Tikhonov được tích hợp nhằm ổn định hóa bài toán ngược thiếu điều kiện, kết hợp giải thuật Gradient liên hợp phi tuyến (Nonlinear Conjugate Gradient - NCG) và phép xấp xỉ tỉ số Rayleigh để tối ưu hóa tham số điều hòa.
Chất lượng hình ảnh được lượng hóa bằng hai tiêu chuẩn khoa học: sai số trung bình bình phương (MSE/RMSE) và Chỉ số chất lượng ảnh phổ quát (Universal Image Quality Index - $Q$). Chỉ số $Q$ có miền giá trị từ -1 đến 1 (giá trị tối ưu bằng 1), phản ánh đồng thời ba thành phần: độ suy giảm tương quan, độ méo cường độ sáng và độ méo tương phản theo đặc tính thị giác con người.
Phương pháp nghiên cứu
Nghiên cứu áp dụng phương pháp mô phỏng số học trên phần mềm MATLAB với dữ liệu tán xạ thu nhận từ mô hình hệ đo tròn gồm các máy phát ($N_t$) và máy thu ($N_r$) phân bố đối xứng xung quanh vùng quan sát (ROI). Miền tái tạo được chia lưới không gian đa mức với kích thước $N \times N = 22 \times 22$ ô vuông (pixel). Quy trình nghiên cứu thực hiện qua 4 kịch bản đo lường để khảo sát tính phụ thuộc của bước lặp chuyển tiếp vào mật độ mẫu:
- Kịch bản số mẫu lớn: $N_t = 44$ máy phát và $N_r = 22$ máy thu.
- Kịch bản số mẫu nhỏ: $N_t = 15$ máy phát và $N_r = 7$ máy thu.
- Kịch bản số mẫu trung bình 1: $N_t = 22$ máy phát và $N_r = 11$ máy thu.
- Kịch bản số mẫu trung bình 2: $N_t = 27$ máy phát và $N_r = 14$ máy thu.
Dữ liệu mô phỏng được kiểm thử trong điều kiện thực tế với mức suy giảm công suất phát an toàn từ 1 mW đến 10 mW, bổ sung nhiễu trắng Gaussian có mức SNR cố định 26 dB. Quy trình thuật toán DF-DBIM đề xuất được chia thành 3 giai đoạn: xác định số bước lặp tối ưu $x$ ở tần số thấp $f_1 = 1$ MHz; khôi phục hàm mục tiêu sơ bộ $O_x$ sau $x$ bước lặp; sử dụng $O_x$ làm giá trị khởi tạo để tiếp tục lặp $(8 - x)$ bước ở tần số cao $f_2 = 2$ MHz. Phương pháp phân tích dữ liệu so sánh trực tiếp sai số toán học và tham số $Q$ giữa kỹ thuật đơn tần số và kỹ thuật đa tần số kết hợp.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Thứ nhất, nghiên cứu phát hiện số bước lặp tối ưu $x$ ở tần số thấp $f_1$ tỷ lệ nghịch với mật độ dữ liệu đo lường thu nhận được. Trong kịch bản 1 với mật độ mẫu cao ($N_t = 44, N_r = 22$), giá trị tối ưu đạt được tại $x = 2$. Ngược lại, ở kịch bản 2 với mật độ mẫu thưa ($N_t = 15, N_r = 7$), thuật toán cần tới $x = 7$ bước lặp ở tần số thấp để tránh phân kỳ. Với cấu hình mẫu trung bình ($N_t = 22, N_r = 11$), giá trị cân bằng lý tưởng được xác định chính xác tại $x = 4$ trên tổng số 8 bước lặp.
Thứ hai, thuật toán kết hợp hai tần số DF-DBIM vượt trội hoàn toàn về chất lượng khôi phục ảnh so với các phương pháp đơn tần số truyền thống. Tại bước lặp thứ 1, tần số $f_1 = 1$ MHz cho chỉ số $Q = 0,8338$, cao hơn 50,1% so với tần số $f_2 = 2$ MHz ($Q = 0,5552$). Khi hoàn thành toàn bộ chu trình 8 bước lặp, phương pháp DF-DBIM đạt chỉ số $Q = 0,9992$, tiệm cận giá trị lý tưởng tuyệt đối, vượt xa mức $Q = 0,9908$ của đơn tần số $f_1$ và $Q = 0,8176$ của đơn tần số $f_2$.
Thứ ba, đồ thị sai số toán học $err$ chứng minh tốc độ hội tụ của DF-DBIM vượt bậc ngay từ bước lặp thứ 5 (thời điểm kích hoạt tần số $f_2$). Trong khi phương pháp đơn tần số $f_2$ có mức sai số dừng ở mức cao 0,5198 do vướng cực tiểu cục bộ, phương pháp DF-DBIM liên tục kéo giảm sai số về mức tối thiểu. Dữ liệu mặt cắt thẳng đứng đi qua tâm của hàm mục tiêu cho thấy biên dạng vận tốc khôi phục của DF-DBIM khớp hoàn toàn với mô hình phân bố vận tốc thực nghiệm từ 2% đến 5%.
Thảo luận kết quả
Nguyên nhân cốt lõi giúp kỹ thuật kết hợp tần số thành công nằm ở cơ chế bù trừ sóng âm trong bài toán tán xạ phi tuyến. Tần số thấp ($f_1 = 1$ MHz) sở hữu bước sóng dài, ít bị ảnh hưởng bởi hiện tượng đảo pha và tán xạ phức tạp, giúp thuật toán nhanh chóng định vị cấu trúc tổng thể và bao hàm năng lượng chính của đối tượng. Khi chuyển tiếp sang tần số cao ($f_2 = 2$ MHz) ở nửa sau chu trình, các bước sóng ngắn sẽ bổ sung các thành phần tần số không gian cao, làm sắc nét đường biên và khôi phục các chi tiết có kích thước nhỏ.
Nếu dữ liệu được biểu diễn qua bảng so sánh đa thông số và đồ thị tiến trình hội tụ 8 bước lặp, người xem có thể quan sát rõ quỹ đạo suy giảm sai số dạng hàm mũ của DF-DBIM so với sự dao động bất ổn định của đơn tần số cao. So với các hệ thống quét lâm sàng quốc tế như CURE hay HUTT vốn bỏ qua hiện tượng nhiễu xạ bậc cao, việc áp dụng DBIM kết hợp điều hòa Tikhonov và đa tần số mang lại độ chính xác định lượng vật lý vượt trội, hạn chế tối đa ảnh giả và đảm bảo an toàn tuyệt đối cho người bệnh.
Đề xuất và khuyến nghị
Nhằm hoàn thiện và đưa kỹ thuật siêu âm cắt lớp kết hợp tần số vào ứng dụng y tế thực tế, các giải pháp cụ thể được đề xuất như sau:
- Phát triển thuật toán tự động thích ứng xác định điểm chuyển tiếp tần số $x$ dựa trên ước lượng ma trận thông tin Fisher, nhằm giảm thêm 25% thời gian tính toán mô phỏng trong vòng 6 đến 12 tháng tới do nhóm nghiên cứu thuật toán xử lý tín hiệu chủ trì thực hiện.
- Thiết kế và chế tạo mảng đầu dò siêu âm đa tần số băng rộng (tích hợp dải tần từ 1 MHz đến 5 MHz) sử dụng vật liệu gốm áp điện composite có độ nhạy cao, mục tiêu hoàn thành thử nghiệm phần cứng trong thời gian 12 đến 18 tháng bởi các trung tâm R&D thiết bị y sinh.
- Mở rộng mô hình giải bài toán ngược từ không gian 2 chiều (2D) lên không gian 3 chiều (3D) toàn phần, hướng tới mục tiêu kiểm soát sai số tái tạo thể tích dưới 1,5% trong lộ trình 18 đến 24 tháng do các phòng thí nghiệm trọng điểm công nghệ điện tử viễn thông đảm trách.
- Xây dựng quy trình thử nghiệm tiền lâm sàng trên các phantom mô vú nhân tạo có cấu trúc bệnh lý phức tạp (kích thước vi vôi hóa dưới 1 mm) với quy mô tối thiểu 50 mẫu thử, triển khai trong 24 tháng tới thông qua chương trình hợp tác giữa viện nghiên cứu và bệnh viện chuyên khoa ung bướu.
Đối tượng nên tham khảo luận văn
Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho nhiều nhóm độc giả chuyên sâu:
- Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Điện tử, Viễn thông và Kỹ thuật Y sinh: Cung cấp tài liệu tham khảo chuẩn mực về thuật toán DBIM, phương pháp tối ưu NCG và mã nguồn chương trình mô phỏng MATLAB phục vụ nghiên cứu luận văn, bài báo khoa học.
- Kỹ sư nghiên cứu và phát triển (R&D) tại các công ty thiết bị y tế: Cung cấp giải pháp công nghệ để tối ưu hóa kiến trúc thu phát phần cứng, cấu hình mảng đầu dò đa kênh và nâng cao độ phân giải cho các dòng máy siêu âm thế hệ mới.
- Bác sĩ chuyên khoa chẩn đoán hình ảnh và chuyên gia y tế: Giúp hiểu rõ cơ chế vật lý định lượng của siêu âm tán xạ ngược, từ đó có thêm góc nhìn khoa học trong việc phân tích hình ảnh và đánh giá sớm khối u tiền ung thư.
- Giảng viên và nhà khoa học trong lĩnh vực Toán ứng dụng và Xử lý ảnh: Bổ sung ca nghiên cứu thực tiễn điển hình về giải quyết bài toán ngược phi tuyến (Inverse Problems) và kỹ thuật điều hòa Tikhonov trong môi trường truyền sóng thực.
Câu hỏi thường gặp
Siêu âm cắt lớp (UCT) vượt trội hơn siêu âm B-mode truyền thống ở những điểm nào?
Siêu âm cắt lớp cho phép tái tạo định lượng chính xác các thông số vật lý của mô như vận tốc truyền âm và hệ số suy giảm âm thay vì chỉ hiển thị độ phản xạ định tính. Kỹ thuật này loại bỏ hoàn toàn nhiễu đốm (speckle), có khả năng phát hiện các tổn thương vi mô nhỏ hơn bước sóng âm và không phụ thuộc vào kỹ năng thao tác góc đặt đầu dò của bác sĩ.
Tại sao phương pháp lặp vi phân Born cần kết hợp hai tần số thay vì dùng đơn tần số?
Nếu chỉ sử dụng đơn tần số cao (2 MHz), thuật toán dễ rơi vào các cực tiểu địa phương khiến hình ảnh bị méo mó và sai số lớn (khoảng 0,5198). Khi kết hợp tần số thấp (1 MHz) ở các bước lặp đầu để tạo khung ảnh gốc ổn định, sau đó dùng tần số cao để tái tạo chi tiết, chỉ số chất lượng $Q$ được nâng lên tới 0,9992.
Việc lựa chọn số bước lặp tối ưu ở tần số thấp phụ thuộc vào yếu tố nào?
Số bước lặp $x$ phụ thuộc trực tiếp vào số lượng kênh thu phát (mật độ dữ liệu mẫu). Khi số mẫu lớn ($N_t = 44, N_r = 22$), hệ thống chỉ cần $x = 2$ bước lặp tần số thấp. Khi số mẫu trung bình ($N_t = 22, N_r = 11$), giá trị tối ưu là $x = 4$. Khi số mẫu quá ít ($N_t = 15, N_r = 7$), thuật toán cần tới $x = 7$ bước để đảm bảo hội tụ.
Chỉ số chất lượng ảnh phổ quát Q đánh giá hình ảnh y tế chính xác hơn chỉ số sai số MSE như thế nào?
Chỉ số sai số MSE chỉ tính toán độ lệch giá trị toán học thuần túy giữa các điểm ảnh mà bỏ qua tương quan không gian. Chỉ số $Q$ kết hợp đồng thời độ mất tương quan, độ biến dạng ánh sáng và độ tương phản quang học. Một ảnh có MSE cao vẫn có thể đạt cấu trúc giải phẫu rõ ràng, và chỉ số $Q$ phản ánh chính xác cảm nhận thị giác này.
Thuật toán DF-DBIM đề xuất xử lý ảnh hưởng của nhiễu môi trường như thế nào?
Nghiên cứu tích hợp phương pháp điều hòa Tikhonov cùng giải thuật Gradient liên hợp phi tuyến (NCG), sử dụng tỉ số Rayleigh để ước lượng tham số làm mịn tự động. Ngay cả khi làm việc trong môi trường chịu tác động của nhiễu trắng Gaussian với tỷ số tín hiệu trên nhiễu SNR = 26 dB, thuật toán vẫn khôi phục chính xác biên độ tán xạ.
Kết luận
- Luận văn đã xây dựng thành công thuật toán lặp vi phân Born kết hợp hai tần số (DF-DBIM), giải quyết triệt để bài toán phân kỳ và nâng cao chất lượng ảnh siêu âm cắt lớp y sinh.
- Tìm ra quy luật phụ thuộc và xác định chính xác số bước lặp chuyển tiếp tối ưu ($x = 4$ trên tổng số 8 bước) cho cấu hình đo lường 22 máy phát và 11 máy thu.
- Nâng cao vượt bậc chất lượng chẩn đoán hình ảnh với chỉ số chất lượng phổ quát $Q$ đạt 0,9992 và giảm thiểu tối đa sai số toán học so với các phương pháp đơn tần số.
- Chứng minh khả năng ứng dụng thực tế của kỹ thuật trong môi trường có nhiễu Gaussian 26 dB và độ tương phản vận tốc âm mô sinh học từ 2% đến 5%.
- Mở ra hướng nghiên cứu chuyển giao công nghệ mới trong chế tạo thiết bị chẩn đoán ung thư sớm không bức xạ với chi phí hợp lý cho ngành y tế Việt Nam.
Kế hoạch phát triển tiếp theo trong giai đoạn 12 đến 24 tháng tới là mở rộng thử nghiệm trên dữ liệu lâm sàng thực tế và tối ưu hóa phần cứng mảng đầu dò đa kênh. Các viện nghiên cứu, trường đại học và doanh nghiệp y tế quan tâm có thể khai thác trực tiếp kết quả luận văn để thúc đẩy thương mại hóa thiết bị siêu âm cắt lớp thế hệ mới.