Tổng quan nghiên cứu

Trong hệ thống chẩn đoán hình ảnh y tế hiện đại, các kỹ thuật như chụp X-quang, chụp cắt lớp vi tính (CT), chụp cộng hưởng từ (MRI) và siêu âm đóng vai trò cốt lõi trong việc phát hiện sớm và điều trị các bệnh lý hiểm nghèo, đặc biệt là ung thư. Mặc dù chụp CT và MRI mang lại độ tương phản cấu trúc mô mềm chi tiết, nhưng chi phí kinh tế cao và nguy cơ nhiễm xạ ion hóa từ tia X vẫn là những rào cản đáng kể đối với việc tầm soát diện rộng. Siêu âm y khoa nổi lên như một giải pháp chẩn đoán an toàn tuyệt đối, không xâm lấn, mức năng lượng bức xạ thấp chỉ từ 1 mW đến 10 mW, và có chi phí chỉ bằng khoảng 10% đến 20% so với chụp MRI.

Tuy nhiên, phương pháp tạo ảnh siêu âm B-mode truyền thống dựa trên nguyên lý phản hồi âm thanh đơn thuần vẫn tồn tại nhược điểm lớn: độ phân giải không gian bị giới hạn bởi bước sóng, chất lượng hình ảnh phụ thuộc chủ quan vào kinh nghiệm của kỹ thuật viên và không thể nhận diện được các tổn thương vi thể nhỏ hơn bước sóng âm. Nhằm khắc phục hạn chế này, kỹ thuật chụp siêu âm cắt lớp (Ultrasound Computed Tomography - USCT) sử dụng nguyên lý tán xạ ngược đã được nghiên cứu sâu rộng. Phương pháp này cho phép tái tạo bản đồ phân bố định lượng vận tốc truyền âm trong mô sống (vốn dao động từ 1460 m/s ở mô mỡ đến 1580 m/s ở mô gan và đạt trung bình 1540 m/s trong mô người), từ đó bộc lộ rõ nét các khối u có kích thước dưới bước sóng.

Luận văn thạc sĩ kỹ thuật điện tử của tác giả Vũ Hoàng Tuấn Hiệp, dưới sự hướng dẫn khoa học của Phó Giáo sư Tiến sĩ Trần Đức Tân tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội (bảo vệ năm 2016), đã tập trung giải quyết bài toán tái tạo ảnh siêu âm cắt lớp thông qua thuật toán lặp vi phân Born (Distorted Born Iterative Method - DBIM) tích hợp kỹ thuật quyết định ngưỡng. Mục tiêu trọng tâm của công trình là giải quyết triệt để hai thách thức lớn nhất của kỹ thuật USCT: tốc độ tính toán chậm và tính bất ổn định trước tác động của nhiễu. Luận văn mở ra tiềm năng to lớn trong việc nâng cao tốc độ hội tụ thêm khoảng 30% đến 40% và duy trì chất lượng khôi phục ảnh định lượng vượt trội cho các ứng dụng sàng lọc ung thư vú trong thực tế 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

Khung lý thuyết của đề tài được xây dựng vững chắc trên nền tảng cơ học sóng âm phi tuyến và lý thuyết bài toán ngược tán xạ sóng. Khi sóng âm truyền qua môi trường cơ thể chứa 60% đến 70% nước, trường áp suất tổng cộng được mô tả qua phương trình sóng Helmholtz vi phân toàn phần. Dựa trên sự biến thiên trở kháng âm học và mật độ vật chất, hàm mục tiêu cần khôi phục biểu diễn sự chênh lệch bình phương số sóng giữa khối tổn thương và môi trường truyền dẫn cơ bản.

Nghiên cứu tập trung vào hai mô hình tái tạo bài toán ngược kinh điển:

  • Phương pháp lặp Born (Born Iterative Method - BIM): Cập nhật hàm truyền Green trong môi trường đồng nhất nền ở từng vòng lặp, có độ ổn định cao nhưng tốc độ hội tụ tương đối chậm.
  • Phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM): Cập nhật hàm Green trong môi trường không đồng nhất biến dạng theo từng bước lặp, giúp đẩy nhanh tốc độ hội tụ phi tuyến bậc một nhưng dễ phân kỳ khi dữ liệu bị nhiễm tạp âm.

Để giải quyết bản chất bài toán không chỉnh (ill-posed) khi nghịch đảo ma trận truyền đạt sóng âm, nghiên cứu ứng dụng mô hình chuẩn hóa Tikhonov kết hợp thuật toán tối ưu hóa Gradient liên hợp phi tuyến (Nonlinear Conjugate Gradient - NCG). Chuẩn hóa Tikhonov với tham số phạt thích nghi được xác định thông qua xấp xỉ tỉ số Rayleigh, giúp cân bằng hoàn hảo giữa độ mịn của nghiệm và năng lượng sai số chuẩn dư tương đối (Relative Residual Error - RRE). Khái niệm hàm mục tiêu âm học, vùng quan tâm (Region of Interest - ROI), và hàm Green hai chiều đóng vai trò là các trụ cột lý thuyết xuyên suốt quá trình thiết lập hệ phương trình tích phân Lippmann-Schwinger.

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

Nghiên cứu sử dụng phương pháp mô phỏng số học thực nghiệm trên môi trường phần mềm tính toán khoa học chuyên sâu MATLAB với dữ liệu mô phỏng trường sóng tán xạ đầy đủ.

  • Cỡ mẫu và cấu hình hệ đo: Không gian khảo sát ROI được chia lưới rời rạc hóa thành ma trận $N \times N$ phần tử ảnh (pixel), tương đương 4096 ô điểm ảnh trong lưới $64 \times 64$ với kích thước mỗi ô pixel đạt mức milimet. Mô hình hệ thống thu phát bao gồm 32 vị trí máy phát sóng siêu âm và 32 máy thu bố trí đối xứng hình tròn bao quanh đối tượng với góc quét toàn phần từ 0 độ đến 360 độ, tạo ra tổng cộng 1024 điểm đo áp suất tán xạ phức hợp.
  • Phương pháp chọn mẫu: Lựa chọn mô hình mẫu số (numerical phantom) đại diện cho mô bệnh học tiêu chuẩn với vật thể hình trụ tròn đặt trong môi trường nước (vận tốc truyền âm nền 1480 m/s). Vùng tổn thương mô phỏng có vận tốc truyền sóng 1540 m/s đến 1580 m/s, tương đương độ tương phản vận tốc âm thanh từ 4% đến 6.7%, đại diện cho các khối u xơ hoặc khối u ác tính giai đoạn sớm.
  • Lý do lựa chọn phương pháp phân tích: Thuật toán DBIM được chọn vì tính ưu việt về tốc độ tính toán so với BIM truyền thống khi xử lý trường sóng tán xạ phức tạp. Tuy nhiên, thay vì giữ nguyên toàn bộ các giá trị áp suất tán xạ đo được ở mọi cảm biến, tác giả tích hợp cơ chế quyết định ngưỡng thích nghi. Cơ chế này loại bỏ các thành phần tín hiệu tán xạ có biên độ thấp hơn ngưỡng sai số hoặc bị chi phối bởi nhiễu trắng Gaussian, giúp giải ma trận gradient liên hợp NCG với độ ổn định tuyệt đối và tránh hiện tượng bẫy cực tiểu địa phương.
  • Dòng thời gian nghiên cứu: Quá trình thiết kế mô hình toán học, lập trình thuật toán và kiểm thử đánh giá được triển khai chặt chẽ, nghiệm thu hoàn chỉnh vào tháng 8 năm 2016.

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

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

Quá trình kiểm nghiệm mô phỏng trên các tập dữ liệu có và không có nhiễu đã đem lại những phát hiện quan trọng mang tính đột phá về hiệu năng tái tạo ảnh:

  1. Tối ưu hóa thời gian tính toán và chu kỳ lặp: Thuật toán DBIM hỗ trợ bởi quyết định ngưỡng giúp hệ thống đạt trạng thái hội tụ nghiệm chỉ sau 5 đến 7 vòng lặp, giảm khoảng 35% số bước lặp so với DBIM nguyên bản (thường cần từ 10 đến 12 vòng lặp để đạt cùng mức sai số).
  2. Nâng cao độ chính xác định lượng hàm mục tiêu: Thuật toán tái tạo chính xác biên độ vận tốc âm của đối tượng với sai số biên độ cực đại dưới 3.2% so với giá trị lý tưởng, bộc lộ ranh giới hình học sắc nét của khối u mà không xuất hiện các vệt mờ biên thể (smearing artifacts).
  3. Khả năng kháng nhiễu vượt bậc: Khi đưa vào mô hình mức nhiễu tạp âm trắng Gaussian lên tới 10% đến 20% công suất tín hiệu, phương pháp đề xuất duy trì chỉ số sai số chuẩn dư tương đối RRE ở mức dưới 0.045, trong khi phương pháp DBIM truyền thống có RRE tăng vọt lên trên 0.12 và xuất hiện phân kỳ cục bộ.
  4. Năng lực khôi phục cấu trúc dưới bước sóng: Với tần số phát sóng siêu âm 2 MHz đến 5 MHz (tương ứng bước sóng khoảng 0.3 mm đến 0.77 mm trong nước), hệ thống có khả năng nhận diện và tái tạo rõ ràng các cấu trúc vi mô kích thước 0.5 mm, vượt qua rào cản nhiễu xạ truyền thống của các máy siêu âm B-mode thông thường.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương pháp DBIM kết hợp quyết định ngưỡng đạt được hiệu năng vượt trội nằm ở việc lọc bỏ các thành phần gradient nhiễu loạn trong không gian nghiệm. Khi tín hiệu áp suất tán xạ đi qua các kênh đo xa hoặc chịu suy giảm biên độ lớn, tỷ số tín hiệu trên nhiễu (SNR) giảm mạnh. Nếu đưa toàn bộ các giá trị này vào ma trận chuẩn hóa Tikhonov, thuật toán NCG sẽ mất rất nhiều bước lặp để triệt tiêu các thành phần sai lệch vô nghĩa. Bằng việc áp dụng ngưỡng biên độ, thuật toán chỉ giữ lại các vector trường sóng mang thông tin cấu trúc đặc trưng nhất, giúp đường hướng dịch chuyển của gradient tiến thẳng về nghiệm tối ưu toàn cục.

Trong quá trình phân tích học thuật, các tập dữ liệu định lượng này có thể được trực quan hóa vô cùng sinh động thông qua biểu đồ suy giảm hàm sai số RRE theo từng bước lặp và bảng so sánh các chỉ số chất lượng ảnh. Cụ thể, dữ liệu có thể được trình bày dưới dạng bảng thống kê chi tiết giá trị RRE và thời gian tính toán (giây) tại mỗi vòng lặp từ bước 1 đến bước 10 giữa ba phương pháp: BIM, DBIM cổ điển và DBIM có quyết định ngưỡng. Đồng thời, biểu đồ đường biểu diễn sự hội tụ sai số kết hợp cùng các bản đồ phân bố vận tốc âm 2D dạng ảnh màu trực quan qua từng bước lặp sẽ minh chứng rõ ràng sự tiệm cận nhanh chóng của ảnh tái tạo với hàm mục tiêu lý tưởng.

So với các công bố quốc tế cùng thời kỳ tập trung vào kỹ thuật lọc không gian tần số, cách tiếp cận quyết định ngưỡng trực tiếp trên miền áp suất tán xạ của tác giả vừa đảm bảo tính tinh giản về mặt thuật toán, vừa giảm thiểu tối đa chi phí tài nguyên phần cứng, tạo tiền đề thuận lợi để nhúng trực tiếp vào các vi xử lý tín hiệu số (DSP) thời gian thực.

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

Dựa trên những phát hiện khoa học từ luận văn, một số giải pháp công nghệ và lộ trình ứng dụng thực tiễn được đề xuất như sau:

  • Song song hóa thuật toán trên kiến trúc phần cứng GPU: Thực hiện lập trình tối ưu hóa các phép nhân ma trận và toán tử tích phân Green trên nền tảng CUDA/OpenCL nhằm nâng tốc độ xử lý khung hình lên mức trên 10 fps (khung hình/giây). Giải pháp này cần được triển khai trong thời gian 12 tháng bởi các nhóm kỹ sư xử lý tín hiệu số.
  • Thử nghiệm trên mô hình phantom sinh học đa cấu trúc: Tiến hành đo đạc thực nghiệm trên các phantom silicon và mô động vật có cấu trúc phức tạp nhằm kiểm chứng độ tin cậy với mục tiêu sai số định lượng vận tốc âm thanh dưới 5%. Kế hoạch dự kiến hoàn thành trong vòng 6 đến 12 tháng dưới sự chủ trì của các phòng thí nghiệm kỹ thuật y sinh.
  • Chế tạo và tích hợp đầu dò mảng vòng 3D đa phần tử: Phát triển hệ thống phần cứng đầu dò mảng pha hình khuyên với số lượng cảm biến từ 128 đến 256 kênh để quét toàn diện không gian 3 chiều, nâng độ phân giải không gian đa lát cắt trong lộ trình 24 tháng bởi các doanh nghiệp sản xuất thiết bị y tế.
  • Xây dựng bộ tiêu chuẩn lọc nhiễu thích nghi tiền xử lý: Chuẩn hóa quy trình tiền xử lý tín hiệu sóng âm thô trước khi nạp vào mô hình DBIM nhằm duy trì tỷ số SNR của hệ đo luôn đạt trên 15 dB. Mục tiêu này cần được nhóm nghiên cứu hoàn thiện trong thời gian 6 tháng tới.

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

Công trình luận văn này cung cấp hàm lượng tri thức chuyên sâu và mang lại giá trị thiết thực cho nhiều nhóm đối tượng:

  • Kỹ sư và nhà nghiên cứu thuật toán thị giác máy tính - xử lý ảnh y tế: Nắm vững cơ chế toán học giải bài toán ngược phi tuyến tính, ứng dụng thuật toán DBIM, NCG và chuẩn hóa Tikhonov vào các hệ thống tái tạo ảnh quang học, ảnh siêu âm và ảnh vi sóng.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Điện tử, Viễn thông, Kỹ thuật Y sinh: Sử dụng làm tài liệu tham khảo chuẩn mực về phương pháp lập mô hình tán xạ trường sóng, cách tổ chức mã nguồn MATLAB xử lý ma trận lớn và phương pháp phân tích hiệu năng hệ thống.
  • Bác sĩ và chuyên gia chẩn đoán hình ảnh: Tiếp cận nền tảng công nghệ siêu âm cắt lớp thế hệ mới, hiểu rõ giá trị định lượng của bản đồ vận tốc âm thanh trong việc sàng lọc ung thư biểu mô tuyến vú với kích thước nhỏ hơn 1 mm.
  • Đơn vị nghiên cứu phát triển và sản xuất thiết bị y tế: Khai thác cấu trúc hệ thống thu phát đối xứng và giải pháp phần mềm tối ưu để định hình thiết kế các dòng máy siêu âm cắt lớp thương mại thế hệ mới với giá thành cạnh tranh.

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

  1. Phương pháp chụp siêu âm cắt lớp USCT có ưu điểm gì vượt trội so với siêu âm B-mode thông thường? Siêu âm B-mode chỉ dựa trên cường độ phản hồi âm cục bộ và phụ thuộc vào góc đặt đầu dò, trong khi USCT thu nhận trường sóng tán xạ toàn phần $360^\circ$. Điều này cho phép tái tạo bản đồ vận tốc âm định lượng chính xác của từng voxel mô, đạt độ phân giải dưới bước sóng và loại bỏ hoàn toàn hiện tượng bóng cản âm thanh.

  2. Tại sao phương pháp lặp vi phân Born (DBIM) được ưu tiên sử dụng hơn phương pháp lặp Born (BIM)? Phương pháp DBIM cập nhật hàm Green phi tuyến theo nền môi trường thay đổi sau mỗi vòng lặp thay vì giữ nguyên nền đồng nhất như BIM. Nhờ đó, DBIM đạt tốc độ hội tụ bậc một nhanh hơn gấp 2 đến 3 lần, rút ngắn đáng kể thời gian tái tạo đối với các mô có độ tương phản âm học biến thiên liên tục.

  3. Kỹ thuật quyết định ngưỡng can thiệp vào giai đoạn nào để tăng tốc độ xử lý? Quyết định ngưỡng được áp dụng trực tiếp lên vector sai phân áp suất tán xạ đo được trước khi cập nhật gradient vi phân NCG. Kỹ thuật này loại bỏ các thành phần tín hiệu mang năng lượng nhiễu vượt trội, giảm độ phức tạp tính toán của ma trận và ngăn ngừa hiện tượng lặp phân kỳ do dữ liệu sai lệch.

  4. Sự khác biệt về vận tốc truyền âm trong các mô cơ thể người ảnh hưởng thế nào đến chẩn đoán? Vận tốc sóng âm thay đổi theo mật độ và độ đàn hồi mô: mô mỡ khoảng 1460 m/s, mô cơ 1545 m/s đến 1630 m/s, trong khi các khối u ác tính thường có mật độ tế bào cao khiến tốc độ âm thanh vượt ngưỡng 1560 m/s. Việc định lượng chính xác vận tốc âm giúp bác sĩ phân biệt bản chất mô lành và khối u ác tính một cách khách quan.

  5. Nghiên cứu này có thể ứng dụng trực tiếp cho các cơ quan nội tạng chứa khí hay xương không? Sóng siêu âm bị suy giảm và tán xạ cực mạnh khi gặp các bề mặt có độ chênh trở kháng âm quá lớn như xương (tốc độ lên tới 4100 m/s) hoặc khí trong phổi (tốc độ chỉ khoảng 650 m/s). Do đó, giải pháp này phù hợp tối ưu nhất cho việc tái tạo ảnh các tạng đặc và mô mềm đồng nhất như tuyến vú, gan hoặc tuyến giáp.

Kết luận

  • Luận văn đã giải quyết thành công bài toán tái tạo ảnh siêu âm cắt lớp độ phân giải cao bằng thuật toán DBIM cải tiến kết hợp cơ chế quyết định ngưỡng.
  • Giảm thiểu hiệu quả từ 30% đến 40% thời gian tính toán và chu kỳ lặp, mở đường cho khả năng xử lý hình ảnh y khoa theo thời gian thực.
  • Nâng cao rõ rệt độ bền vững của nghiệm trước tạp âm môi trường, giữ chỉ số sai số chuẩn dư tương đối RRE luôn ở mức an toàn dưới 0.045.
  • Đạt năng lực định lượng chính xác vận tốc truyền âm và bộc lộ hoàn hảo các cấu trúc tổn thương bệnh lý có kích thước nhỏ hơn bước sóng phát xạ.
  • Cung cấp mã nguồn thực nghiệm MATLAB hoàn chỉnh, đóng vai trò là tài liệu tham khảo học thuật và nền tảng kỹ thuật giá trị cao cho cộng đồng nghiên cứu xử lý tín hiệu y sinh.

Trong giai đoạn tiếp theo, việc tích hợp thuật toán trên phần cứng xử lý song song và thử nghiệm lâm sàng thực tế sẽ là bước đi quyết định để thương mại hóa công nghệ USCT. Quý độc giả, các nhà nghiên cứu và kỹ sư y sinh quan tâm có thể khai thác toàn diện các mô hình toán học và kết quả thực nghiệm trong luận văn để tiếp tục mở rộng và phát triển các hệ thống chẩn đoán hình ảnh tiên tiến phục vụ sức khỏe cộng đồng.