Tổng quan nghiên cứu

Ung thư gan hiện đang là một trong những thách thức y tế hàng đầu tại Việt Nam với gánh nặng bệnh tật và tỷ lệ tử vong rất cao. Theo báo cáo dịch tễ học, năm 2018 Việt Nam ghi nhận gần 165.000 ca mắc ung thư mới trên quy mô dân số 96,5 triệu người, và con số này ước tính vượt mức 189.000 ca vào năm 2020. Trong đó, ung thư gan dẫn đầu về tỷ lệ mắc với hơn 25.000 trường hợp ghi nhận trong năm 2018. Tỷ lệ sống sót sau 5 năm của bệnh nhân ung thư gan trung bình chỉ đạt khoảng 9%. Đáng chú ý, nếu tổn thương được phát hiện và can thiệp ở giai đoạn sớm, tỷ lệ sống trên 5 năm có thể đạt 19%, nhưng con số này sẽ giảm mạnh xuống 6,5% ở giai đoạn hai và chỉ còn khoảng 3,5% khi bệnh tiến triển sang giai đoạn muộn.

Phần lớn các ca ung thư gan phát triển trên nền bệnh lý xơ gan mãn tính. Do đó, việc tầm soát và chẩn đoán định lượng mức độ xơ hóa đóng vai trò then chốt giúp bác sĩ đưa ra phác đồ can thiệp kịp thời. Mặc dù sinh thiết mô gan được xem là tiêu chuẩn vàng truyền thống, kỹ thuật xâm lấn này tiềm ẩn nhiều biến chứng nghiêm trọng xảy ra ở 1% đến 5% bệnh nhân, với tỷ lệ tử vong được ghi nhận từ 1 trên 1.000 đến 1 trên 10.000 ca. Nhằm khắc phục triệt để các rủi ro trên, kỹ thuật tạo ảnh siêu âm đàn hồi sóng biến dạng đã ra đời như một giải pháp không xâm lấn tiên tiến.

Tuy nhiên, chất lượng tạo ảnh sóng biến dạng trong thực tế lâm sàng luôn bị suy giảm bởi nhiễu Gauss ngẫu nhiên và hiện tượng phản xạ sóng tại các mặt phân cách mô. Nghiên cứu này được thực hiện tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội vào năm 2020 nhằm giải quyết bài toán ước lượng mô đun shear phức trong môi trường mô mềm sinh học 2D có nhiễu. Mục tiêu trọng tâm là kết hợp phương pháp sai phân hữu hạn trong miền thời gian, bộ lọc thích nghi trung bình bình phương tối thiểu và thuật toán biến đổi ngược đại số Helmholtz để tái tạo chính xác độ đàn hồi và độ nhớt của mô xơ hóa, giảm thiểu sai số đo lường xuống dưới 8,02%.

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 dựa trên nền tảng vật lý truyền sóng âm và mô hình lưu biến học mô mềm. Đầu tiên, mô hình viscoelastic Kelvin-Voigt được áp dụng để biểu diễn mô đun shear phức CSM trong miền tần số. Mô hình này kết hợp đồng thời hai đặc tính cơ học cơ bản của mô mềm sinh học bao gồm độ đàn hồi tĩnh (ký hiệu là $\mu$, đơn vị Pascal) và độ nhớt động học (ký hiệu là $\eta$, đơn vị Pascal giây). Phương trình mô tả mối liên hệ giữa vận tốc truyền sóng biến dạng, mật độ khối của mô và các thành phần cơ học cho thấy sóng biến dạng lan truyền với vận tốc càng cao thì mô tổn thương có độ cứng càng lớn.

Khung lý thuyết thứ hai là hệ phương trình truyền sóng Navier và định luật phản xạ sóng âm Snell. Trong môi trường không đồng nhất, khi sóng biến dạng di chuyển qua mặt phân cách giữa hai tổ chức mô có trở kháng âm khác nhau, một phần năng lượng sóng sẽ bị phản xạ ngược lại với hệ số phản xạ R phụ thuộc vào mức độ chênh lệch trở kháng. Trở kháng âm của mô gan bình thường đạt khoảng 1,65 Mrayl, mô cơ đạt 1,7 Mrayl trong khi mô mỡ là 1,38 Mrayl.

Các khái niệm cốt lõi được vận dụng xuyên suốt bao gồm: Sóng biến dạng (Shear Wave) với đặc tính dao động vuông góc với phương truyền sóng và vận tốc chỉ bằng khoảng 60% sóng dọc; Kỹ thuật xung lực bức xạ âm (ARFI) và kích thích rung kim cơ học từ 50 Hz đến 450 Hz; Thuật toán lọc thích nghi trung bình bình phương tối thiểu (LMS); và Thuật toán biến đổi ngược đại số Helmholtz (AHI) dùng để giải trực tiếp bài toán ngược trường sóng.

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 kết hợp kiểm chứng định lượng trên môi trường tính toán khoa học chuyên sâu. Dữ liệu nghiên cứu được khởi tạo từ mô hình mô phỏng trường sóng 2D đại diện cho vùng mô mềm sinh học có kích thước không gian 120 nhân 120 mm. Bên trong môi trường nền có bố trí một khối xơ gan hình tròn tại tọa độ tâm (40 mm, 40 mm) với bán kính 15 mm.

Cỡ mẫu dữ liệu khảo sát bao gồm toàn bộ lưới không gian 2D với 120 nhân 120 điểm nút, tương đương 14.400 vị trí đo với bước lưới lấy mẫu đồng đều 1 mm theo cả hai trục tọa độ X và Y. Để đảm bảo độ tin cậy cao nhất, phương pháp chọn mẫu phân tích tập trung vào vùng khảo sát ROI có kích thước 75 nhân 75 mm tương ứng với các tọa độ từ 5 mm đến 79 mm, giúp loại bỏ vùng chuyển tiếp bất ổn định sát nguồn kích thích và vùng suy giảm biên độ quá mức ở rìa ngoài.

Phương pháp phân tích được thực hiện qua quy trình 3 giai đoạn chặt chẽ:

  1. Mô hình hóa quá trình lan truyền sóng biến dạng trong môi trường dị thể bằng phương pháp sai phân hữu hạn trong miền thời gian (FDTD) với nguồn kích thích rung kim tần số 200 Hz, biên độ 5 mm và mật độ mô 1.000 kg/m³.
  2. Ứng dụng bộ lọc thích nghi LMS với chiều dài bộ lọc xác định để triệt tiêu thành phần nhiễu Gauss cộng sinh ra trong quá trình thu nhận vận tốc hạt bằng siêu âm Doppler.
  3. Sử dụng thuật toán AHI chuyển đổi phương trình sóng từ miền thời gian sang phương trình Helmholtz trong miền tần số để ước lượng trực tiếp giá trị độ đàn hồi $\mu$ và độ nhớt $\eta$ tại từng điểm ảnh.

Lý do lựa chọn chuỗi giải thuật FDTD-LMS-AHI là nhờ khả năng mô phỏng chính xác sóng biến dạng trong môi trường không đồng nhất, loại bỏ nhiễu hiệu quả mà không làm méo tín hiệu gốc, đồng thời tiết kiệm đáng kể thời gian tính toán so với các thuật toán ước lượng hợp lẽ cực đại (MLEF), tạo điều kiện tối ưu để tích hợp vào các thiết bị siêu âm thời gian thực. Toàn bộ kịch bản mô phỏng và phân tích dữ liệu được hoàn thiện trong timeline nghiên cứu năm 2020.

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

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

Nghiên cứu đã ghi nhận ba phát hiện mang tính đột phá về khả năng tái tạo ảnh định lượng và kiểm soát sai số trong môi trường có nhiễu:

Thứ nhất, việc kết hợp bộ lọc LMS trước khi thực hiện thuật toán AHI mang lại hiệu quả vượt trội trong việc triệt tiêu nhiễu Gauss. Khi đánh giá trên toàn bộ vùng ảnh 75 nhân 75 mm, sai số chuẩn hóa của việc ước lượng độ đàn hồi đạt mức dưới 6,4%, trong khi sai số chuẩn hóa của độ nhớt đạt mức dưới 8,02%. Các giá trị này phản ánh độ tương đồng rất cao giữa bản đồ đàn nhớt ước lượng với mô hình lý tưởng ban đầu (độ đàn hồi nền 6.000 Pa, khối xơ 8.900 Pa; độ nhớt nền 1,2 Pa.s, khối xơ 2,7 Pa.s).

Thứ hai, nghiên cứu phát hiện sự phụ thuộc rõ rệt của độ chính xác ước lượng vào khoảng cách truyền sóng. Tại các vị trí cách nguồn kích thích vượt quá 80 mm, biên độ vận tốc hạt của sóng biến dạng bị suy giảm tự nhiên rất mạnh do tính chất tiêu tán năng lượng của mô mềm. Hiện tượng này làm tỷ số tín hiệu trên nhiễu giảm sút, dẫn đến sai số ước lượng tại các vùng xa tăng cao. Do đó, việc giới hạn khung ảnh tái tạo ở kích thước 75 nhân 75 mm là ngưỡng tối ưu giúp bảo toàn độ tin cậy chẩn đoán.

Thứ ba, nghiên cứu đã định lượng hóa cụ thể tác động tiêu cực của hiện tượng phản xạ sóng biến dạng tại bề mặt ranh giới khối xơ gan. Khi có hiện tượng phản xạ xảy ra, sóng tới và sóng phản xạ giao thoa làm biến dạng phân bố vận tốc hạt đo được. Kết quả đánh giá chỉ ra rằng: sai số chuẩn hóa của độ đàn hồi khi có phản xạ tăng từ 0,0430 (4,30%) lên 0,0754 (7,54%), tương ứng mức tăng tuyệt đối 3,24%; trong khi sai số chuẩn hóa của độ nhớt tăng vọt từ 0,1412 (14,12%) lên 0,2695 (26,95%), tức mức tăng tuyệt đối lên tới 12,83%.

Thảo luận kết quả

Các kết quả thực nghiệm mô phỏng có thể được minh họa trực quan thông qua các đồ thị phân bố vận tốc hạt cắt ngang không gian theo đường Line 10 (vùng mô lành không qua u) và Line 40 (đường cắt xuyên tâm khối xơ gan). Dữ liệu biểu diễn dạng biểu đồ đường cho thấy rõ sau khi đi qua bộ lọc thích nghi LMS, các dao động ngẫu nhiên do nhiễu Gauss trên tín hiệu vận tốc hạt tại tọa độ (40, 40) đã được làm mượt hoàn hảo, bám sát dạng sóng hình sin chuẩn. Trên các bản đồ nhiệt 2D, khối xơ gan hình tròn bán kính 15 mm được tái tạo sắc nét về cả vị trí không gian, đường biên kích thước lẫn giá trị độ cứng 8.900 Pa và độ nhớt 2,7 Pa.s.

Nguyên nhân dẫn đến sự gia tăng sai số mạnh ở tham số độ nhớt (tăng 12,83% khi có phản xạ) so với độ đàn hồi (chỉ tăng 3,24%) xuất phát từ bản chất toán học của mô hình Kelvin-Voigt trong thuật toán AHI. Độ nhớt $\eta$ được ước lượng thông qua phần ảo của số sóng phức và tỷ lệ nghịch trực tiếp với tần số góc dao động $\omega_0$. Do đó, hiện tượng giao thoa sóng phản xạ gây lệch pha cục bộ sẽ tác động trực tiếp và nặng nề hơn lên thành phần ảo so với thành phần thực đại diện cho độ đàn hồi $\mu$.

Khi so sánh với các nghiên cứu trước đây sử dụng bộ lọc tổ hợp hợp lẽ cực đại (MLEF) trong môi trường đồng nhất của nhóm tác giả Orescanin Marko hay các nghiên cứu mở rộng của tác giả Shigao Chen, phương pháp tiếp cận FDTD kết hợp LMS và AHI trong luận văn này thể hiện ưu thế vượt trội về mặt tài nguyên xử lý. Phương pháp AHI giải trực tiếp phương trình vi phân thay vì phải chạy các vòng lặp tối ưu hóa phi tuyến phức tạp của MLEF, giúp giảm thời gian tính toán xuống hàng chục lần, giải quyết bài toán nghẽn cổ chai trong việc xử lý tín hiệu siêu âm thời gian thực.

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

Dựa trên các kết quả đạt được, luận văn đưa ra bốn đề xuất hành động cụ thể nhằm thúc đẩy việc ứng dụng công nghệ vào thực tiễn y tế:

  1. Tích hợp khối giải thuật LMS-AHI vào phần cứng máy siêu âm: Các kỹ sư R&D thiết bị y tế tại các doanh nghiệp công nghệ cần lập trình nhúng trực tiếp thuật toán lọc LMS và phép biến đổi AHI lên các chip xử lý tín hiệu số DSP hoặc FPGA chuyên dụng. Mục tiêu cần đạt là rút ngắn thời gian xử lý ảnh 2D xuống dưới 0,5 giây mỗi khung hình với mức sai số đàn hồi duy trì dưới 6,5%, hoàn thành triển khai trong lộ trình 12 đến 18 tháng.

  2. Phát triển module thuật toán bù trừ sóng phản xạ đa hướng: Các nhóm nghiên cứu xử lý tín hiệu y sinh tại các trường đại học kỹ thuật cần nghiên cứu bổ sung thuật toán tách sóng phản xạ dựa trên biến đổi phân rã trường sóng (Wavefield Separation). Mục tiêu trọng tâm là giảm sai số chuẩn hóa độ nhớt trong môi trường có phản xạ từ mức 26,95% hiện tại xuống dưới ngưỡng 10%, thực hiện trong giai đoạn 2021 đến 2023.

  3. Thực nghiệm kiểm chuẩn trên Phantom sinh học đa lớp: Các phòng thí nghiệm kỹ thuật y tế cần phối hợp với các bệnh viện chuyên khoa để chế tạo các mẫu mô phỏng sinh học (tissue-mimicking phantoms) có độ cứng đa dạng từ 2 kPa đến 20 kPa tương ứng các giai đoạn xơ gan F0 đến F4 theo thang Metavir. Thử nghiệm trên cỡ mẫu tối thiểu 200 lượt đo thực nghiệm để đánh giá độ lặp lại và độ bền vững của thuật toán trước năm 2024.

  4. Chuẩn hóa quy trình chẩn đoán xơ gan không xâm lấn: Bộ Y tế và các hội đồng chuyên môn chẩn đoán hình ảnh cần xây dựng hướng dẫn kỹ thuật kết hợp siêu âm đàn hồi sóng biến dạng định lượng 2 tham số (độ cứng và độ nhớt) vào quy trình tầm soát bệnh gan mãn tính. Mục tiêu hướng tới là giảm từ 30% đến 40% các chỉ định sinh thiết gan không cần thiết trên toàn quốc trong giai đoạn 2023 đến 2025.

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

Nội dung và giải pháp kỹ thuật của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng trọng tâm:

  1. Bác sĩ chẩn đoán hình ảnh và chuyên gia nội gan mật: Luận văn cung cấp cơ sở vật lý chuyên sâu về cơ chế hình thành sóng biến dạng, giúp bác sĩ hiểu rõ ý nghĩa định lượng của độ đàn hồi (kPa) và độ nhớt (Pa.s). Từ đó, bác sĩ có thể ứng dụng dữ liệu siêu âm đàn hồi để đánh giá chính xác các giai đoạn xơ gan từ F0 đến F4, phân biệt u lành tính và ác tính mà không cần lạm dụng sinh thiết xâm lấn gây đau cho bệnh nhân.

  2. Kỹ sư kỹ thuật y sinh và lập trình viên xử lý ảnh y tế: Tài liệu là nguồn tham khảo chi tiết về thuật toán sai phân FDTD, giải thuật thích nghi LMS và kỹ thuật biến đổi ngược AHI. Đây là khung kiến trúc mẫu giúp các kỹ sư phát triển phần mềm xử lý tín hiệu trên các dòng máy siêu âm thế hệ mới.

  3. Học viên cao học và nghiên cứu sinh ngành Điện tử - Viễn thông: Cung cấp phương pháp luận chuẩn mực trong việc mô hình hóa toán học các bài toán truyền sóng trong môi trường sinh học không đồng nhất, xử lý nhiễu Gauss và đánh giá sai số định lượng trên phần mềm mô phỏng.

  4. Nhà quản lý bệnh viện và các đơn vị cung ứng thiết bị y tế: Giúp các nhà quản lý nắm bắt xu hướng công nghệ tạo ảnh y học tiên tiến, làm cơ sở khoa học để thẩm định, lựa chọn và đầu tư các hệ thống máy siêu âm đàn hồi sóng biến dạng hiện đại phục vụ nâng cao chất lượng khám chữa bệnh.

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

Siêu âm sóng biến dạng có điểm gì khác biệt so với siêu âm kiểu B truyền thống?
Siêu âm kiểu B truyền thống chỉ thể hiện cấu trúc hình thái học thông qua độ phản hồi âm xám của mô, khiến các khối u có cùng độ phản hồi rất khó phân biệt lành tính hay ác tính. Siêu âm sóng biến dạng bổ sung thông tin cơ học định lượng về độ đàn hồi (kPa) và độ nhớt (Pa.s), giúp đo lường trực tiếp độ cứng mô bệnh lý một cách khách quan.

Tại sao việc ước lượng đồng thời cả độ nhớt và độ đàn hồi lại quan trọng trong chẩn đoán xơ gan?
Mô mềm sinh học có tính chất đàn nhớt phức tạp, không chỉ thay đổi độ cứng mà còn thay đổi tính nhớt khi bị tổn thương cấu trúc ngoại bào. Ước lượng cả hai tham số theo mô hình Kelvin-Voigt giúp phản ánh toàn diện tình trạng biến đổi mô học qua từng giai đoạn F0 đến F4, nâng cao độ đặc hiệu trong chẩn đoán bệnh lý gan.

Bộ lọc thích nghi LMS giải quyết vấn đề gì trong bài toán tạo ảnh siêu âm?
Trong quá trình đo vận tốc hạt bằng chùm tia siêu âm Doppler, tín hiệu luôn bị pha tạp bởi nhiễu Gauss ngẫu nhiên làm biến dạng dạng sóng. Bộ lọc thích nghi LMS giúp ước lượng và triệt tiêu chính xác thành phần nhiễu này, hạ thấp sai số chuẩn hóa độ đàn hồi xuống dưới 6,4% và độ nhớt dưới 8,02% trước khi giải bài toán AHI.

Hiện tượng phản xạ sóng biến dạng tác động như thế nào đến độ chính xác của hình ảnh?
Khi sóng biến dạng gặp bề mặt phân cách giữa mô lành và khối xơ có trở kháng âm chênh lệch, sóng phản xạ sinh ra sẽ giao thoa với sóng tới làm thay đổi vận tốc hạt thực tế. Hệ quả là sai số chuẩn hóa độ đàn hồi tăng thêm 3,24% và sai số độ nhớt tăng mạnh thêm 12,83%, làm giảm độ phân giải biên khối u.

Thuật toán kết hợp FDTD-LMS-AHI có thể áp dụng cho việc tạo ảnh thời gian thực không?
Hoàn toàn khả thi. Khác với phương pháp ước lượng khả năng cực đại MLEF đòi hỏi lặp phi tuyến tốn kém bộ nhớ và thời gian, thuật toán AHI giải trực tiếp phương trình đại số Helmholtz trong miền tần số kết hợp bộ lọc LMS tối ưu, cho phép giảm thiểu khối lượng tính toán để tích hợp vào các chip DSP xử lý thời gian thực.

Kết luận

  • Luận văn đã xây dựng thành công quy trình tạo ảnh định lượng mô đun shear phức (CSM) 2D cho mô mềm sinh học thông qua việc tích hợp phương pháp FDTD, bộ lọc thích nghi LMS và biến đổi ngược đại số Helmholtz AHI.
  • Kiểm chứng mô phỏng chứng minh bộ lọc LMS kiểm soát xuất sắc nhiễu Gauss, duy trì sai số chuẩn hóa độ đàn hồi dưới 6,4% và độ nhớt dưới 8,02% trên vùng ảnh tối ưu kích thước 75 nhân 75 mm.
  • Khảo sát toàn diện và định lượng hóa ảnh hưởng của hiện tượng phản xạ sóng biến dạng, chỉ rõ mức tăng sai số 3,24% ở độ đàn hồi và 12,83% ở độ nhớt do hiện tượng giao thoa trường sóng tại ranh giới khối xơ.
  • Đóng góp giải pháp công nghệ có chi phí tính toán thấp, khắc phục hạn chế của các phương pháp tối ưu hóa lặp truyền thống, mở ra tiềm năng lớn cho các thiết bị siêu âm chẩn đoán y tế thời gian thực.
  • Kế hoạch phát triển tiếp theo tập trung vào việc thử nghiệm thuật toán trên mô hình 3D thực nghiệm và phantom đa tầng trong giai đoạn 2021-2023. Hãy liên hệ với các đơn vị nghiên cứu kỹ thuật y sinh chuyên sâu để chuyển giao và ứng dụng công nghệ tạo ảnh đàn nhớt tiên tiến này vào hệ thống thiết bị y tế lâm sàng ngay hôm nay!