Tổng quan nghiên cứu

Trong bối cảnh tỷ lệ mắc các bệnh lý hiểm nghèo như ung thư ngày càng gia tăng, nhu cầu về các kỹ thuật chẩn đoán hình ảnh y sinh an toàn, chi phí hợp lý và độ chính xác cao trở nên cấp thiết hơn bao giờ hết. Các phương pháp chẩn đoán truyền thống hiện nay đều bộc lộ những hạn chế nhất định. Chụp cắt lớp vi tính (CT) mang lại hình ảnh sắc nét nhưng khiến cơ thể người bệnh phải hấp thụ liều tia X tương đương từ 6 tháng đến 3 năm phơi nhiễm bức xạ tự nhiên. Chụp cộng hưởng từ (MRI) có độ tương phản mô mềm vượt trội nhưng chi phí đắt đỏ và chống chỉ định với bệnh nhân mang dị vật kim loại. Trong khi đó, siêu âm B-mode thông thường an toàn tuyệt đối với vận tốc âm trung bình trong cơ thể khoảng 1540 m/s và dải tần lâm sàng từ 2 MHz đến 17 MHz, nhưng lại bị giới hạn về độ phân giải không gian và không thể phát hiện các dị vật hoặc vi tổn thương có kích thước nhỏ hơn bước sóng.

Để khắc phục rào cản trên, kỹ thuật chụp ảnh siêu âm cắt lớp (Ultrasound Computed Tomography - UCT) dựa trên nguyên lý tán xạ ngược đã được nghiên cứu nhằm tái tạo định lượng 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. Tuy nhiên, bài toán ngược phi tuyến trong siêu âm cắt lớp thường gặp hiện tượng phân kỳ hoặc mắc kẹt tại các cực tiểu cục bộ khi chỉ sử dụng một tần số sóng âm duy nhất. Nghiên cứu này tập trung giải quyết bài toán nâng cao chất lượng ảnh siêu âm cắt lớp thông qua việc đề xuất phương pháp kết hợp hai tần số (Dual-Frequency Distorted Born Iterative Method - DF-DBIM).

Mục tiêu cốt lõi của nghiên cứu là tối ưu hóa số bước lặp chuyển đổi giữa tần số thấp 1 MHz và tần số cao 2 MHz, từ đó tái tạo chính xác hàm mục tiêu của đối tượng tán xạ có đường kính 3.1 mm trong môi trường nhiễu Gaussian với 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 to lớn khi nâng chỉ số chất lượng ảnh tổng quát từ mức 0.8176 lên 0.9992, đồng thời giảm thiểu sai số tái tạo xuống mức tối ưu. Đây là bước đệm quan trọng nhằm thúc đẩy ứng dụng thương mại của hệ thống siêu âm cắt lớp an toàn, không nhiễm xạ trong tầm soát và chẩn đoán y khoa hiện đại.

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 trên nền tảng lý thuyết truyền sóng âm trong môi trường không đồng nhất và bài toán tán xạ ngược y sinh. Khi chùm sóng siêu âm tương tác với mô sinh học, trường áp suất tổng $p(r)$ được biểu diễn bằng tổng của sóng tới $p_{inc}(r)$ và sóng tán xạ $p_{sc}(r)$ thông qua phương trình tích phân Lippmann-Schwinger kết hợp hàm Green $G(r, r')$. Hai lý thuyết trụ cột được sử dụng để giải bài toán ngược tán xạ bao gồm phương pháp lặp Born (Born Iterative Method - BIM) và phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM). Trong đó, DBIM được lựa chọn làm giải pháp trọng tâm nhờ khả năng 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 vượt bậc so với BIM truyền thống.

Mô hình nghiên cứu định lượng vật thể thông qua hàm mục tiêu $O(r)$, phản ánh sự sai khác bình phương vận tốc truyền sóng giữa môi trường chuẩn (nước với vận tốc $c_0 = 1480$ m/s) và mô bệnh lý $c_1(r)$ (thường dao động từ 1460 m/s ở mô mỡ đến 1630 m/s ở mô cơ). Để rời rạc hóa các phương trình tích phân, phương pháp phần tử hữu hạn dạng Moment (Method of Moments - MoM) được triển khai kết hợp với xấp xỉ Born bậc một. Đánh giá độ chính xác của hình ảnh không chỉ dựa trên sai số trung bình bình phương (MSE) hay căn sai số toàn phương trung bình (RMSE), mà còn áp dụng Chỉ số chất lượng ảnh phổ quát $Q$ (Universal Image Quality Index). Chỉ số $Q$ biến thiên trong khoảng từ -1 đến 1 (với giá trị 1 là tái tạo hoàn hảo), tích hợp đồng thời ba yếu tố của hệ thống thị giác con người bao gồm: độ tổn hao tương quan, độ méo cường độ sáng và độ méo tương phản trên dải tần số khảo sát $f_1 = 1$ MHz và $f_2 = 2$ MHz.

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

Nghiên cứu sử dụng phương pháp mô phỏng số thực nghiệm toán tin trên không gian hai chiều với vùng quan tâm (ROI) được chia thành lưới kích thước $N \times N = 22 \times 22$, tương đương 484 điểm ảnh (pixel). Cỡ mẫu quan sát được khảo sát linh hoạt qua nhiều kịch bản thu phát khác nhau: cấu hình nhiều mẫu ($N_t = 44$ máy phát, $N_r = 22$ máy thu), cấu hình ít mẫu ($N_t = 15$ máy phát, $N_r = 7$ máy thu), và cấu hình mẫu trung bình ($N_t = 22$ máy phát, $N_r = 11$ máy thu). Phương pháp chọn mẫu bố trí các đầu thu phát đối xứng hình tròn xung quanh vật thể được lựa chọn nghiêm ngặt nhằm triệt tiêu hiện tượng dịch pha tín hiệu trong quá trình quét dữ liệu.

Lý do lựa chọn phương pháp phân tích bài toán ngược thông qua thuật toán Gradient liên hợp phi tuyến (Nonlinear Conjugate Gradient - NCG) kết hợp chuẩn hóa Tikhonov là nhằm khắc phục triệt để tính chất ill-posed (bài toán đặt không chỉnh) và ảnh hưởng của nhiễu nền. Tham số điều chuẩn được xác định tự động bằng phương pháp lũy thừa lặp với tỉ số xấp xỉ Rayleigh. Toàn bộ quy trình tính toán được thực hiện theo tiến trình 8 bước lặp tổng thể ($\text{sum_iter} = 8$), trong đó thuật toán đề xuất sẽ thực hiện $x$ bước lặp đầu tiên ở tần số thấp $f_1 = 1$ MHz để định vị vùng biên và hình học tổng thể, sau đó chuyển giao kết quả hàm mục tiêu làm giá trị khởi tạo cho $(8 - x)$ bước lặp tiếp theo ở tần số cao $f_2 = 2$ MHz dưới điều kiện nhiễu trắng Gaussian với SNR cố định ở mức 26 dB và độ chênh lệch vận tốc sóng từ 2% đến 5%.

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 và thực nghiệm số đã đem lại ba phát hiện mang tính bước ngoặt trong kỹ thuật tái tạo ảnh siêu âm cắt lớp:

Thứ nhất, nghiên cứu xác lập được mối quan hệ có tính quy luật giữa số lượng mẫu thu phát và số bước lặp tối ưu $x$ ở tần số thấp $f_1$. Khi số lượng mẫu đo lớn ($N_t = 44, N_r = 22$), thuật toán chỉ cần $x = 2$ bước lặp ở tần số 1 MHz là đủ dữ liệu để hội tụ; khi số mẫu đo giảm 3 lần ($N_t = 15, N_r = 7$), thuật toán đòi hỏi $x = 7$ bước lặp ở tần số thấp để tránh lỗi phân kỳ (giá trị sai số NaN ở các bước đầu). Với cấu hình mẫu trung bình ($N_t = 22, N_r = 11$), giá trị tối ưu được xác định chính xác tại $x = 4$, giúp cân bằng hoàn hảo giữa độ ổn định và chi phí tính toán.

Thứ hai, tần số thấp $f_1 = 1$ MHz thể hiện ưu thế vượt trội trong giai đoạn khởi tạo ban đầu. Tại bước lặp thứ nhất, chỉ số chất lượng ảnh $Q$ của tần số $f_1$ đạt 0.8338, cao hơn 50.18% so với giá trị $Q = 0.5552$ khi sử dụng đơn tần số cao $f_2 = 2$ MHz. Đến bước lặp thứ hai, chất lượng ảnh ở $f_1$ tiếp tục duy trì mức áp đảo với $Q = 0.9604$ so với $Q = 0.7101$ của $f_2$.

Thứ ba, sự kết hợp đa tần số DF-DBIM tạo ra bước nhảy vọt về chất lượng tái tạo ở các vòng lặp sau. Từ bước lặp thứ 5 (bắt đầu chuyển sang tần số $f_2 = 2$ MHz trên nền tảng hàm mục tiêu đã tối ưu), phương pháp DF-DBIM đạt chỉ số $Q = 0.9927$, vượt qua cả hai phương pháp đơn tần ($Q_{f1} = 0.9911$ và $Q_{f2} = 0.8046$). Tại bước lặp thứ 8 cuối cùng, DF-DBIM đạt chỉ số chất lượng cảm quan gần như tuyệt đối $Q = 0.9992$, trong khi phương pháp DBIM đơn tần chỉ dừng lại ở $Q = 0.9908$ (với $f_1$) và $Q = 0.8176$ (với $f_2$).

Thảo luận kết quả

Cơ chế vật lý đằng sau sự thành công của phương pháp DF-DBIM nằm ở sự tương hỗ hoàn hảo giữa hai dải bước sóng. Tần số thấp 1 MHz sở hữu bước sóng dài hơn, giúp giảm tính phi tuyến của hàm mục tiêu, mở rộng bán kính hội tụ của thuật toán và hạn chế tối đa nguy cơ rơi vào cực tiểu cục bộ. Sau khi cấu trúc đại thể của đối tượng được định hình vững chắc sau 4 bước lặp đầu tiên, việc đưa tần số cao 2 MHz vào 4 bước lặp tiếp theo sẽ tận dụng bước sóng ngắn để khôi phục các thành phần phổ tần số không gian cao, từ đó làm sắc nét đường biên và tái hiện chính xác độ tương phản âm học của vật thể.

Dữ liệu nghiên cứu khi được biểu diễn trực quan qua đồ thị hội tụ sai số và biểu đồ đường cắt 1D qua tâm đối tượng đã minh chứng rõ nét luận điểm này. Trên đồ thị sai số toàn phương, đường cong sai số của phương pháp đề xuất DF-DBIM luôn nằm thấp hơn hẳn so với đường sai số của hai phương pháp đơn tần số $f_1$ và $f_2$, chứng minh tốc độ hội tụ nhanh hơn và sai số tàn dư nhỏ hơn. Khi phân tích biểu đồ mặt cắt thẳng đứng đi qua trung tâm hàm mục tiêu, biên dạng tín hiệu khôi phục của DF-DBIM bám sát đường tín hiệu lý tưởng với độ lệch biên độ dưới 1%, trong khi các phương pháp đơn tần xuất hiện độ gợn sóng và méo biên dạng đáng kể. So với các hệ thống quét cắt lớp siêu âm thương mại như CURE, HUTT vốn thường bỏ qua hiện tượng nhiễu xạ bậc cao, phương pháp DF-DBIM mang lại độ chính xác định lượng vượt trội ngay cả trong điều kiện nhiễu nền 26 dB.

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

Dựa trên những kết quả thực nghiệm đạt được, luận văn đưa ra 4 nhóm giải pháp và khuyến nghị hành động mang tính ứng dụng cao:

  1. Chuẩn hóa quy trình tối ưu hóa số bước lặp đa tần số: Các nhóm nghiên cứu thuật toán xử lý ảnh y sinh cần áp dụng quy tắc xác định tham số $x$ linh hoạt theo mật độ đầu dò. Với hệ thống phần cứng 32 đến 64 kênh thu phát, cần thiết lập số bước lặp tần số thấp chiếm từ 40% đến 50% tổng chu kỳ lặp để đạt mục tiêu chỉ số $Q \ge 0.995$, hoàn thành tích hợp trong vòng 6 tháng tới.
  2. Nâng cấp thuật toán sang dải tần số liên tục (Multi-frequency DBIM): Các viện nghiên cứu điện tử viễn thông nên mở rộng thuật toán từ việc kết hợp 2 tần số rời rạc (1 MHz và 2 MHz) lên dải 3 đến 5 tần số quét liên tục (từ 1 MHz đến 5 MHz). Mục tiêu là nâng cao độ phân giải không gian của ảnh cắt lớp lên dưới 1 mm, triển khai thử nghiệm trong lộ trình 12 tháng.
  3. Phát triển phần cứng xử lý song song trên nền tảng GPU: Doanh nghiệp sản xuất thiết bị y tế cần phối hợp với các chuyên gia kỹ thuật máy tính để tăng tốc tính toán ma trận tán xạ Moment và thuật toán NCG bằng công nghệ CUDA hoặc OpenCL. Mục tiêu là giảm thời gian tái tạo một lát cắt $1024 \times 1024$ pixel từ vài phút xuống dưới 5 giây nhằm đáp ứng tiêu chuẩn hiển thị thời gian thực (tối thiểu 15 đến 24 khung hình/giây) trước năm 2026.
  4. Triển khai thử nghiệm lâm sàng trên mô hình mô nhân tạo (Phantom): Các bệnh viện và trung tâm chẩn đoán hình ảnh cần phối hợp với nhóm tác giả tiến hành đo đạc thực nghiệm trên các phantom mô phỏng u tuyến vú và tuyến giáp. Mục tiêu là kiểm chứng sai số đo vận tốc truyền âm dưới 1.5% trong môi trường lâm sàng thực tế trong giai đoạn 2024-2025.

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

Công trình nghiên cứu này mang lại giá trị học thuật và ứng dụng thực tiễn sâu rộng cho 4 nhóm đối tượng trọng tâm:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Điện tử - Viễn thông, Kỹ thuật Y sinh: Tài liệu cung cấp hệ thống cơ sở toán học chặt chẽ về bài toán ngược tán xạ sóng, phương pháp lặp vi phân Born và mã nguồn mô phỏng chi tiết bằng MATLAB trên lưới $22 \times 22$, phục vụ đắc lực cho các đề tài phát triển giải thuật xử lý tín hiệu y học.
  • Kỹ sư phát triển phần mềm và thuật toán tại các công ty thiết bị y tế: Nắm bắt kiến trúc tối ưu hóa kết hợp đa tần số và phương pháp điều chuẩn Tikhonov với tỉ số Rayleigh, giúp nâng cấp thuật toán tái tạo ảnh cho các dòng máy siêu âm thế hệ mới với chỉ số chất lượng thị giác $Q > 0.99$.
  • Bác sĩ và chuyên gia chẩn đoán hình ảnh: Có thêm cơ sở khoa học để đánh giá tiềm năng của kỹ thuật siêu âm cắt lớp định lượng, hiểu rõ sự khác biệt giữa hình ảnh B-mode quy ước và bản đồ phân bố vận tốc âm học (1460 m/s đến 4100 m/s) trong việc phát hiện sớm khối u ác tính.
  • Nhà quản lý y tế và các đơn vị đầu tư trang thiết bị bệnh viện: Tham khảo mô hình phân tích chi phí - hiệu năng của hệ thống siêu âm cắt lớp không nhiễm xạ, từ đó định hướng đầu tư các giải pháp chẩn đoán an toàn, tiết kiệm chi phí phần cứng từ 30% đến 40% so với việc lắp đặt hệ thống MRI hoặc CT đắt tiền.

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

Kỹ thuật kết hợp 2 tần số DF-DBIM vượt trội hơn DBIM đơn tần ở điểm nào?

DF-DBIM kết hợp ưu thế của cả hai dải sóng âm: tần số thấp 1 MHz giúp định hình cấu trúc đại thể và mở rộng vùng hội tụ, trong khi tần số cao 2 MHz làm sắc nét chi tiết biên. Nhờ đó, phương pháp đạt chỉ số chất lượng $Q = 0.9992$, vượt xa mức $0.9908$ và $0.8176$ của phương pháp đơn tần trong môi trường nhiễu 26 dB.

Tại sao việc lựa chọn số bước lặp x ở tần số thấp lại mang tính quyết định?

Nếu chọn sai giá trị $x$, quá trình tái tạo có thể bị phân kỳ hoặc cho chất lượng ảnh kém hơn cả phương pháp đơn tần. Với hệ thống 22 máy phát và 11 máy thu, việc chọn chính xác $x = 4$ trong tổng số 8 bước lặp giúp thuật toán đạt sai số tối thiểu và chất lượng cảm quan cao nhất.

Chỉ số chất lượng ảnh Q có ưu điểm gì so với sai số MSE truyền thống?

Sai số MSE chỉ đo lường sự sai khác toán học đơn thuần giữa các điểm ảnh mà bỏ qua cảm nhận thị giác. Chỉ số $Q$ biến thiên từ -1 đến 1, phản ánh chính xác 3 yếu tố thị giác con người gồm độ tương quan, độ sáng và độ tương phản, giúp đánh giá trung thực độ sắc nét của hình ảnh y sinh.

Cấu hình máy phát và máy thu ảnh hưởng như thế nào đến kết quả tái tạo?

Số lượng đầu dò quyết định mật độ dữ liệu thô. Khi có nhiều mẫu (44 máy phát, 22 máy thu), thuật toán hội tụ rất nhanh chỉ sau 2 bước lặp ở tần số thấp. Khi số mẫu ở mức trung bình (22 phát, 11 thu), cấu hình đối xứng giúp loại bỏ hiện tượng lệch pha và tái tạo hoàn hảo hàm mục tiêu.

Kỹ thuật siêu âm cắt lớp này có thể thay thế hoàn toàn chụp X-quang và CT không?

Siêu âm cắt lớp DF-DBIM là giải pháp bổ trợ và thay thế lý tưởng cho các mô mềm như tuyến vú, tuyến giáp nhờ ưu điểm không nhiễm xạ tia X và chi phí thấp. Tuy nhiên, đối với các cấu trúc xương có độ cản âm quá lớn (vận tốc trên 4000 m/s), phương pháp vẫn cần kết hợp thêm các kỹ thuật chẩn đoán chuyên sâu khác.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp kết hợp 2 tần số DF-DBIM nhằm nâng cao vượt bậc chất lượng ảnh siêu âm cắt lớp y sinh.
  • Tối ưu hóa thành công tham số bước lặp $x = 4$ ở tần số 1 MHz trong tổng số 8 chu kỳ lặp đối với cấu hình 22 máy phát và 11 máy thu.
  • Nâng chỉ số chất lượng ảnh phổ quát $Q$ lên mức kỷ lục 0.9992, vượt trội hoàn toàn so với việc sử dụng đơn tần số trong điều kiện nhiễu 26 dB.
  • Chứng minh tính ưu việt của việc tái tạo định lượng vận tốc âm học, mở ra hướng chẩn đoán khối u sớm không xâm lấn và an toàn tuyệt đối.
  • Xác lập lộ trình thử nghiệm trên mô hình thực nghiệm và tối ưu hóa phần cứng xử lý song song trong vòng 1 đến 2 năm tới.

Công trình nghiên cứu này đánh dấu bước tiến quan trọng trong lĩnh vực xử lý tín hiệu và chẩn đoán hình ảnh y tế. Để tìm hiểu chi tiết về thuật toán và ứng dụng giải pháp DF-DBIM vào hệ thống thiết bị thực tế, quý độc giả và các đơn vị nghiên cứu có thể liên hệ hợp tác chuyên môn và chuyển giao công nghệ ngay hôm nay.