Tổng quan nghiên cứu

Ung thư vú hiện là một trong những bệnh lý ác tính phổ biến hàng đầu và là nguyên nhân gây tử vong cao nhất do ung thư ở phụ nữ trên toàn cầu cũng như tại Việt Nam. Việc phát hiện sớm các khối u có kích thước nhỏ đóng vai trò quyết định trong việc nâng cao tỷ lệ điều trị thành công lên trên 90%. Tuy nhiên, các kỹ thuật chẩn đoán hình ảnh truyền thống vẫn bộc lộ nhiều hạn chế: chụp X-quang tuyến vú tiềm ẩn nguy cơ phơi nhiễm bức xạ ion hóa, chụp cộng hưởng từ (MRI) có chi phí đắt đỏ và thời gian kéo dài, trong khi siêu âm phản xạ 2D (B-mode) thông thường chỉ cung cấp thông tin định tính với độ phân giải giới hạn bởi bước sóng, khó phân biệt được ranh giới giữa mô lành, mô xơ và mô mỡ khi có sự hòa lẫn cấu trúc vi mô.

Trước thực trạng đó, công trình nghiên cứu tập trung giải quyết bài toán tái tạo hình ảnh định lượng mô học thông qua kỹ thuật siêu âm cắt lớp tán xạ ngược (Ultrasound Tomography). Mục tiêu trọng tâm của đề tài là xây dựng mô hình toán học tích hợp biến đổi mật độ mô và đề xuất giải thuật lặp Born biến dạng kết hợp tần số kép (Dual-Frequency Distorted Born Iterative Method - DF-DBIM). Nghiên cứu được triển khai mô phỏng trên vùng tán xạ kích thước đường kính 10 mm với lưới chia ma trận 40x40 pixel, khảo sát dải tần số từ 0,64 MHz đến các tần số bậc cao, sử dụng hệ thống đo từ 12 máy thu đến 21 máy phát bố trí đối xứng.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc trích xuất đồng thời hai thông số vật lý quan trọng: tốc độ truyền âm và mật độ vật chất của mô. Kết quả thực nghiệm mô phỏng ghi nhận giải thuật DF-DBIM giúp sai số chuẩn hóa (Relative Residual Error - RRE) giảm nhanh chóng từ mức 0,72 xuống còn 0,12 chỉ sau 4 vòng lặp đầu tiên, đồng thời rút ngắn 35% đến 45% thời gian xử lý so với các phương pháp đơn tần số truyền thống. Đây là tiền đề công nghệ quan trọng để phát triển các thiết bị siêu âm cắt lớp an toàn, công suất phát thấp từ 1 mW đến 10 mW, hỗ trợ sàng lọc ung thư vú không xâm lấn hiệu quả tại các cơ sở y tế.

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 sóng âm trong môi trường không đồng nhất, lý thuyết tán xạ ngược phi tuyến và giải thuật tối ưu hóa số học. Khung lý thuyết bao gồm các mô hình và khái niệm nòng cốt sau:

Thứ nhất, phương trình sóng Helmholtz tổng quát cho môi trường biến đổi mật độ: Khi sóng siêu âm lan truyền trong mô sinh học, sự không đồng nhất về mật độ $\rho(\vec{r})$ và tốc độ âm $c(\vec{r})$ làm xuất hiện cả thành phần tán xạ đơn cực (phụ thuộc vào hệ số nén $\kappa$) và tán xạ lưỡng cực (phụ thuộc vào mật độ $\rho$). Phương trình vi phân mô tả áp suất âm $p(\vec{r})$ kết hợp số sóng $k_0$ trong môi trường nền đồng nhất (như nước có trở kháng $Z_0 = 1,48 \times 10^6 \text{ kg/m}^2\text{s}$) và hàm nguồn được biểu diễn dưới dạng tích phân liên kết hàm Green.

Thứ hai, hàm mục tiêu định lượng (Object Function): Để tái tạo cấu trúc hình học và đặc tính vật lý của khối u, hàm mục tiêu $O(\vec{r})$ được thiết lập bao gồm sự sai khác bình phương số sóng và toán tử gradient mật độ theo công thức: $$O(\vec{r}) = \left[ k^2(\vec{r}) - k_0^2 \right] - \rho^{1/2}(\vec{r}) \nabla \cdot \left[ \rho^{-1/2}(\vec{r}) \nabla \right]$$ Trong đó, việc xét đến yếu tố mật độ mang lại độ tương phản âm học vượt trội so với mô hình chỉ xét tốc độ âm thanh.

Thứ ba, nguyên lý kết hợp tần số kép (DF-DBIM): Phương pháp giải bài toán tán xạ ngược thông qua hai tần số $f_1$ (tần số thấp, ví dụ 0,64 MHz) và $f_2$ (tần số cao). Tần số thấp đảm bảo tính lồi và độ hội tụ toàn cục của giải thuật về vùng lân cận nghiệm thực tế, trong khi tần số cao đảm nhận việc tái thiết các biên phân giải không gian chi tiết mà không gây phân kỳ.

Thứ tư, kỹ thuật tối ưu hóa và khử nhiễu: Áp dụng phương pháp điều chuẩn Tikhonov nhằm ổn định nghiệm của bài toán không chỉnh (ill-posed problem), kết hợp với thuật toán Gradient liên hợp phi tuyến (Nonlinear Conjugate Gradient - NCG) và xấp xỉ tỉ số Rayleigh để tối ưu hóa bước nhảy ma trận trong không gian đa chiều.

+-------------------------------------------------------------+
|               TRƯỜNG SÓNG SIÊU ÂM TỚI                      |
|       (Tần số thấp f1 = 0,64 MHz + Tần số cao f2)           |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|                BÀI TOÁN TÁN XẠ NGƯỢC                        |
|       Trích xuất tín hiệu áp suất phân tán p_sc(r)          |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|            MÔ HÌNH HÓA VÀ RỜI RẠC HÓA (MoM)                 |
|       Thiết lập ma trận hàm Green G_r, G_s (40x40 pixel)    |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|           GIẢI THUẬT TẦN SỐ KÉP (DF-DBIM) + NCG             |
|   - Tần số thấp f1: Định vị vùng hội tụ nghiệm gốc          |
|   - Tần số cao f2: Tinh chỉnh độ phân giải vi mô            |
|   - Chuẩn hóa Tikhonov khử nhiễu & tối ưu hóa Rayleigh      |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|          HÌNH ẢNH MẬT ĐỘ VÀ TỐC ĐỘ ÂM TÁI TẠO               |
|      (Phát hiện u kích thước < 2 mm, sai số RRE = 0,12)     |
+-------------------------------------------------------------+

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 kết hợp giải thuật nghịch đảo trên dữ liệu tán xạ trường sóng với các thông số kỹ thuật được kiểm soát chặt chẽ:

  • Nguồn dữ liệu và thiết lập cấu hình đo: Dữ liệu tán xạ sóng âm được tạo lập từ các mô hình hình học và mô phỏng thực nghiệm đối tượng hình trụ đồng tâm đại diện cho cấu trúc khối u tuyến vú. Vùng quan tâm (Region of Interest - ROI) có kích thước đường kính 10 mm, được chia lưới rời rạc thành ma trận $40 \times 40$ pixel (tương đương 1600 phần tử tính toán). Kích thước mỗi phần tử pixel $h$ được thiết lập nhỏ hơn một phần mười bước sóng của chùm tia tới nhằm thỏa mãn điều kiện xấp xỉ liên tục.
  • Phương pháp chọn mẫu và phân bố cảm biến: Hệ thống mô phỏng bố trí theo cấu hình tròn đối xứng bao gồm 21 vị trí máy phát ($N_t = 21$) và 12 đến 16 vị trí máy thu ($N_r = 12 - 16$). Việc lựa chọn cấu hình đối xứng xoay quanh tâm vật thể giúp triệt tiêu hiện tượng lệch pha tín hiệu, đảm bảo thu thập đầy đủ dữ liệu trường tán xạ đa hướng $360^\circ$.
  • Phương pháp phân tích số học: Sử dụng phương pháp Moment (Method of Moments - MoM) để chuyển đổi phương trình vi tích phân phi tuyến sang hệ phương trình ma trận tuyến tính cục bộ. Lý do lựa chọn giải thuật DF-DBIM kết hợp phương pháp NCG là khả năng giải quyết triệt để tính chất phi tuyến cao của tán xạ ngược, giúp thuật toán vượt qua các cực tiểu cục bộ và đạt tốc độ hội tụ nhanh hơn gấp 2 đến 3 lần so với phương pháp lặp Born cổ điển (BIM). Toàn bộ quá trình tính toán lặp được thực hiện qua chuỗi 8 chu kỳ vòng lặp liên tiếp để đánh giá tính ổn định của nghiệm.

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à phân tích định lượng trên các kịch bản thực nghiệm đã mang lại những phát hiện quan trọng về hiệu năng của giải thuật DF-DBIM trong việc tạo ảnh mật độ:

Thứ nhất, việc bổ sung tham số biến đổi mật độ $\rho(\vec{r})$ vào hàm mục tiêu đã giải quyết được điểm nghẽn của các mô hình siêu âm truyền thống. Khi kiểm tra trên mô hình mô phỏng với chênh lệch tốc độ truyền sóng 20% và mật độ thay đổi liên tục, ảnh tái tạo mật độ phản ánh chính xác ranh giới mô xơ và mô mỡ, trong khi các phương pháp bỏ qua yếu tố mật độ làm biến dạng hình học khối u lên tới 28%.

Thứ hai, tốc độ hội tụ và độ chính xác của giải thuật DF-DBIM vượt trội so với DBIM đơn tần số. Trong kịch bản kiểm thử với lưới $40 \times 40$ pixel, sai số chuẩn hóa (RRE) giảm mạnh từ giá trị ban đầu 0,72 xuống 0,38 ở vòng lặp thứ nhất, tiếp tục giảm xuống 0,21 ở vòng lặp thứ hai, đạt 0,15 ở vòng lặp thứ ba và ổn định ở mức 0,12 sau vòng lặp thứ tư (tương ứng mức giảm sai số tổng thể 83,3%).

Thứ ba, hiệu quả phân giải không gian và khả năng chống nhiễu được cải thiện rõ rệt khi kết hợp hai tần số. Với cấu hình 20 máy phát và 16 máy thu ($N_t = 20, N_r = 16$), sử dụng riêng tần số $f_1 = 0,64\text{ MHz}$ cho ảnh có đường bao mờ và sai số dư sau 8 vòng lặp là 0,28; sử dụng riêng tần số cao $f_2$ khiến giải thuật mất 14 vòng lặp và dễ rơi vào nghiệm cục bộ khi độ lệch pha $\Delta\phi > 0,5\pi$. Khi kết hợp đồng thời $f_1 + f_2$, sai số giảm về mức 0,08 sau 8 vòng lặp, khôi phục hoàn chỉnh các cấu trúc có kích thước nhỏ hơn 1,5 mm.

Thứ tư, tính ổn định của giải thuật khi thay đổi tỷ số tần số $f_{\min}/f_0$ trong dải từ 0,1 đến 0,9. Nghiên cứu xác nhận rằng khi thiết lập ngưỡng dung sai DBIM ở mức 0,01 kết hợp thuật toán NCG, sai số tái tạo bình phương trung bình (RMSE) luôn duy trì ở mức thấp và ổn định đối với các bán kính hình trụ từ $\lambda_0$ đến $4\lambda_0$.

Thảo luận kết quả

Cơ chế nền tảng giúp phương pháp DF-DBIM đạt hiệu quả vượt trội nằm ở sự phối hợp tương hỗ giữa hai miền tần số. Sóng âm tần số thấp $f_1$ thâm nhập sâu vào môi trường với độ suy giảm thấp, cung cấp thông tin nền tảng về dạng phân bố không gian tổng thể của hàm mục tiêu, đưa giá trị khởi tạo của vòng lặp vào vùng hội tụ của nghiệm thực. Tiếp đó, sóng âm tần số cao $f_2$ với bước sóng ngắn sẽ thực hiện nhiệm vụ khôi phục các thành phần tần số không gian cao, làm sắc nét các đường biên của khối u.

Dữ liệu thực nghiệm của nghiên cứu có thể được trực quan hóa tối ưu thông qua bảng ma trận so sánh sai số và biểu đồ phân rã RMSE qua các vòng lặp:

Thông số kịch bản Đơn tần số thấp ($f_1$) Đơn tần số cao ($f_2$) Kết hợp tần số ($f_1 + f_2$)
Cấu hình phát/thu ($N_t \times N_r$) $21 \times 12$ $21 \times 12$ $20 \times 16$
Sai số chuẩn hóa ban đầu (Vòng 1) 0,72 0,85 0,38
Sai số chuẩn hóa cuối (Vòng 8) 0,28 0,22 0,08
Thời gian tính toán trung bình 120 giây 310 giây 75 giây
Khả năng phân tách biên dạng Thấp (đường biên nhòe) Trung bình (dễ nhiễu hạt) Rất cao (trùng khớp biên thực)

Trên biểu đồ suy giảm sai số, đường biểu diễn của DF-DBIM có độ dốc đứng từ vòng lặp 1 đến vòng lặp 4 và tiệm cận đường lý tưởng nhanh hơn 2,5 lần so với DBIM đơn tần. So với các công trình nghiên cứu trước đây vốn thường bỏ qua thành phần mật độ hoặc chỉ áp dụng lý thuyết tán xạ Born bậc 1 tuyến tính, phương pháp DF-DBIM đề xuất đã giải quyết thành công phương trình vi phân phi tuyến đầy đủ. Điều này giúp loại bỏ hoàn toàn các giả định tuyến tính hóa phi thực tế, mở rộng phạm vi áp dụng cho các cấu trúc mô có độ tương phản âm học lớn.

Ý nghĩa thực tiễn của phát hiện này đặc biệt quan trọng trong việc ứng dụng vào các hệ thống máy tính chẩn đoán lâm sàng tại bệnh viện. Việc giảm số lượng điểm ảnh cần xử lý trung gian và rút ngắn thời gian lặp giúp các máy tính cấu hình tiêu chuẩn có thể dựng lại hình ảnh cắt lớp 2D/3D trong thời gian thực, đáp ứng yêu cầu khắt khe của quá trình thăm khám lâm sàng.

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

Dựa trên các kết quả đạt được, nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật và lộ trình triển khai nhằm đưa giải thuật siêu âm cắt lớp mật độ vào ứng dụng thực tế:

  • Thiết kế và chuẩn hóa phần cứng đầu dò mảng pha đa tần số: Cần phát triển các đầu dò siêu âm mảng ma trận tròn (Circular Ring Arrays) tích hợp dải tần rộng từ 0,5 MHz đến 5,0 MHz với số lượng kênh tối thiểu 32 máy phát và 32 máy thu. Target metric: Nâng tỷ số tín hiệu trên nhiễu (SNR) của hệ đo lên trên 25 dB tại chiều sâu khảo sát 50 mm; Timeline: 12 đến 18 tháng; Chủ thể thực hiện: Các kỹ sư phần cứng thiết bị y tế phối hợp cùng các phòng thí nghiệm kỹ thuật điện tử.
  • Tối ưu hóa thuật toán DF-DBIM trên nền tảng tính toán song song GPU: Chuyển đổi toàn bộ các phép nhân ma trận hàm Green và bước lặp Gradient liên hợp phi tuyến sang môi trường lập trình song song CUDA hoặc OpenCL. Target metric: Giảm thời gian tái tạo hình ảnh cho khối thể tích xuống dưới 3 giây/mặt cắt; Timeline: 6 đến 12 tháng; Chủ thể thực hiện: Nhóm chuyên gia xử lý tín hiệu và lập trình hệ thống nhúng y sinh.
  • Xây dựng cơ sở dữ liệu đối chuẩn mô học tuyến vú đa thông số: Thiết lập bộ dữ liệu chuẩn về trở kháng âm, hệ số suy giảm âm và mật độ mô thực tế của phụ nữ Việt Nam thông qua các mẫu sinh thiết giải phẫu bệnh. Target metric: Thu thập và chuẩn hóa dữ liệu từ ít nhất 500 ca bệnh lâm sàng; Timeline: 24 tháng; Chủ thể thực hiện: Các bệnh viện chuyên khoa ung bướu phối hợp với các viện nghiên cứu y sinh học.
  • Ban hành quy chuẩn an toàn năng lượng và quy trình chẩn đoán lâm sàng: Hoàn thiện khung hướng dẫn vận hành thiết bị siêu âm cắt lớp, bảo đảm công suất chùm tia $I_{\text{SPTA}} < 100\text{ mW/cm}^2$ và mật độ năng lượng dưới $50\text{ J/cm}^2$ theo tiêu chuẩn an toàn của Viện Siêu âm Y khoa Hoa Kỳ (AIUM). Target metric: Đạt 100% tiêu chuẩn kiểm định an toàn bức xạ phi ion hóa; Timeline: 18 tháng; Chủ thể thực hiện: Cơ quan quản lý trang thiết bị y tế và Hội đồng đạo đức y sinh học.

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

Nội dung và giải pháp kỹ thuật trong luận văn mang lại giá trị thực tiễn cao cho nhiều nhóm đối tượng chuyên môn:

  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Y sinh, Điện tử - Viễn thông: Cung cấp tài liệu chuyên sâu về phương pháp giải bài toán ngược phi tuyến, kỹ thuật xử lý tín hiệu âm học và thuật toán lặp vi phân Born. Use case: Áp dụng khung lý thuyết toán học của đề tài để phát triển các giải thuật tái tạo ảnh cắt lớp điện trở (EIT) hoặc cắt lớp vi sóng (Microwave Tomography).
  • Kỹ sư nghiên cứu và phát triển (R&D) thiết bị chẩn đoán hình ảnh: Nắm bắt kiến trúc hệ thống đo tán xạ siêu âm và phương pháp tối ưu hóa ma trận hàm Green. Use case: Ứng dụng mô hình bố trí đầu dò 20 phát - 16 thu để thiết kế hệ thống phần cứng máy siêu âm cắt lớp vú thế hệ mới.
  • Bác sĩ chẩn đoán hình ảnh và chuyên gia ung bướu: Hiểu rõ bản chất vật lý của sự tương phản mô dựa trên mật độ và tốc độ âm thanh, bổ sung góc nhìn định lượng bên cạnh hình ảnh B-mode truyền thống. Use case: Phối hợp đánh giá các tổn thương xơ hóa vi mô có kích thước dưới 2 mm khó quan sát trên ảnh siêu âm thông thường.
  • Các nhà quản lý y tế và đơn vị đầu tư công nghệ bệnh viện: Cung cấp luận cứ khoa học về hiệu quả kinh tế và tính an toàn của công nghệ siêu âm cắt lớp không sử dụng tia X. Use case: Xây dựng dự án đầu tư trang thiết bị tầm soát sớm ung thư vú quy mô lớn tại các trung tâm y tế dự phòng.

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

Tại sao việc tái tạo thông số mật độ lại quan trọng hơn việc chỉ tái tạo tốc độ truyền âm?
Tốc độ truyền âm giữa mô lành và mô mỡ đôi khi có sự tương đồng lớn, đặc biệt khi các mô xơ hòa lẫn với mô mỡ trong tuyến vú. Việc tái tạo đồng thời mật độ $\rho(\vec{r})$ giúp xác định chính xác trở kháng âm $Z = \rho \cdot c$, tạo ra độ tương phản âm học vượt trội hơn 30%, giúp phân biệt rõ ràng giữa khối u ác tính và các mô mềm xung quanh.

Phương pháp DF-DBIM khắc phục nhược điểm gì của giải thuật DBIM truyền thống?
DBIM truyền thống sử dụng đơn tần số thường gặp mâu thuẫn: tần số thấp cho độ phân giải kém, trong khi tần số cao dễ mắc kẹt ở cực tiểu cục bộ và đòi hỏi thời gian tính toán kéo dài. DF-DBIM kết hợp tần số thấp (định vị vùng hội tụ) và tần số cao (tăng độ sắc nét), giúp giảm 83,3% sai số sau 4 vòng lặp và rút ngắn 40% thời gian xử lý.

Thiết bị siêu âm cắt lớp tán xạ có gây tác dụng phụ sinh học nguy hiểm cho bệnh nhân không?
Công nghệ này hoàn toàn an toàn do sử dụng sóng siêu âm chẩn đoán với công suất phát chỉ từ 1 mW đến 10 mW, cường độ $I_{\text{SPTA}} < 100\text{ mW/cm}^2$, thấp hơn rất nhiều so với ngưỡng gây tổn thương nhiệt hoặc hiện tượng tạo lỗ vi mô. Thiết bị không sử dụng bức xạ ion hóa, cho phép thực hiện thăm khám định kỳ an toàn cho cả phụ nữ mang thai.

Cấu hình bao nhiêu đầu phát và thu là tối ưu cho việc tái tạo ma trận ảnh $40 \times 40$ pixel?
Qua các kịch bản thực nghiệm, cấu hình gồm 20 đến 21 máy phát ($N_t$) kết hợp 12 đến 16 máy thu ($N_r$) bố trí đối xứng xung quanh vùng khảo sát 10 mm là tối ưu nhất. Cấu hình này cung cấp đủ số lượng phương trình độc lập để giải ma trận 1600 phần tử, vừa đảm bảo triệt tiêu hiện tượng dịch pha vừa tiết kiệm tài nguyên tính toán.

Làm thế nào để thuật toán xử lý được nhiễu đo lường trong điều kiện lâm sàng thực tế?
Nghiên cứu tích hợp phương pháp điều chuẩn Tikhonov và giải thuật Gradient liên hợp phi tuyến (NCG). Bằng cách thiết lập tham số giảm chấn Rayleigh và khống chế ngưỡng dung sai ở mức 0,01, thuật toán loại bỏ hiệu quả các thành phần nhiễu đo ngẫu nhiên, giúp đường sai số hội tụ ổn định và mịn màng qua chuỗi 8 vòng lặp tính toán.

Kết luận

  • Xây dựng thành công mô hình toán học giải bài toán tán xạ ngược siêu âm có xét đến sự biến thiên mật độ mô sinh học, vượt qua giới hạn của các phương pháp tạo ảnh phản xạ truyền thống.
  • Đề xuất và kiểm chứng hoàn thiện giải thuật tần số kép DF-DBIM kết hợp kỹ thuật tối ưu hóa NCG, giải quyết triệt để sự đánh đổi giữa tốc độ hội tụ và độ phân giải không gian.
  • Kết quả mô phỏng trên lưới ma trận $40 \times 40$ pixel khẳng định sai số chuẩn hóa giảm ấn tượng từ 0,72 xuống 0,12 chỉ sau 4 vòng lặp, nâng cao độ chính xác nhận dạng khối u kích thước dưới 2 mm.
  • Khẳng định tính an toàn sinh học vượt trội của hệ thống siêu âm công suất thấp (1 - 10 mW), mở ra giải pháp sàng lọc ung thư vú không bức xạ chi phí thấp.
  • Thiết lập giải pháp kỹ thuật tối ưu hóa phần cứng và giảm thời gian tính toán, tạo điều kiện thuận lợi cho việc thương mại hóa trên các máy tính y tế thông thường.

Luận văn đã đóng góp một bước tiến khoa học quan trọng trong lĩnh vực xử lý tín hiệu và chẩn đoán hình ảnh y sinh học. Trong giai đoạn tiếp theo từ năm 2024 đến năm 2026, các nhóm nghiên cứu cần tập trung mở rộng giải thuật sang mô hình không gian 3 chiều với kích thước ma trận $512 \times 512$ pixel và triển khai thử nghiệm trên mẫu sinh học thực nghiệm. Các nhà khoa học và kỹ sư y sinh hãy đẩy mạnh ứng dụng giải thuật DF-DBIM vào các hệ thống siêu âm cắt lớp thế hệ mới nhằm nâng cao năng lực phát hiện sớm ung thư vú vì sức khỏe cộng đồng.