Tổng quan về luận án

Nghiên cứu "Nghiên cứu phát triển một số giải thuật nâng cao chất lượng tạo ảnh siêu âm cắt lớp" của nghiên cứu sinh Trần Quang Huy (chuyên ngành Kỹ thuật điện tử, mã số 62520203, Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn của PGS. Trần Đức Tân) giải quyết bài toán cấp thiết trong kỹ thuật chẩn đoán hình ảnh y sinh: phát hiện sớm ung thư vú thông qua phương pháp siêu âm cắt lớp (Ultrasound Tomography – UST) dựa trên lý thuyết tán xạ ngược. Bối cảnh dịch tễ học ghi nhận hàng năm trên thế giới có khoảng 1,4 triệu ca mắc mới ung thư vú và 458.000 trường hợp tử vong; tại Việt Nam, mỗi năm ghi nhận khoảng 7.000 ca mới với tỷ lệ tử vong lên tới 35% – cao hơn đáng kể so với các nước phát triển. Văn bản luận án khẳng định rõ: "Ung thư vú nếu được phát hiện và điều trị sớm thì tỷ lệ chữa khỏi bệnh đạt tới 90% và chất lượng sống của bệnh nhân được tăng lên rõ rệt."

Trong thực hành lâm sàng, kỹ thuật X-quang tuyến vú (Mammography) bộc lộ nhược điểm nghiêm trọng khi chẩn đoán trên phụ nữ dưới 50 tuổi có mô vú dày đặc (dense breast tissue). Đúng như luận án nhấn mạnh: "Các mô dày đặc không cung cấp được sự thay đổi tốc độ truyền âm cần thiết để tạo được ảnh của khối u nhỏ. Trong khi đó, kỹ thuật siêu âm cắt lớp (ultrasound tomography) lại thực hiện được điều này." Đồng thời, ảnh siêu âm B-mode thông thường chỉ mang tính định tính, phụ thuộc chủ quan vào người thực hiện và chỉ phát hiện được u khi kích thước đã lớn. Ngược lại, kỹ thuật siêu âm cắt lớp sử dụng phương pháp lặp vi phân Born (Distorted Born Iterative Method – DBIM) cho phép tái tạo định lượng bản đồ phân bố tốc độ truyền sóng âm, phát hiện tổn thương vi mô dưới 5 mm.

Tuy nhiên, rào cản cốt lõi khiến các thiết bị UST tán xạ ngược chưa thể thương mại hóa rộng rãi nằm ở độ phức tạp tính toán khổng lồ $O(N_t N_r N^2)$ (với $N_t$ máy phát, $N_r$ máy thu, và $N \times N$ điểm ảnh), độ nhạy cảm cao với nhiễu và thời gian quét dữ liệu kéo dài. Để vượt qua khoảng trống học thuật này, luận án xác lập 2 câu hỏi nghiên cứu và 2 giả thuyết khoa học:

  • RQ1: Làm thế nào để giải quyết sự đánh đổi giữa tốc độ hội tụ toàn cục và độ phân giải không gian trong phương pháp DBIM mà không làm bùng nổ chi phí tính toán?
  • RQ2: Làm thế nào để giảm thiểu số lượng đầu dò thu phát và số phép đo tán xạ mà vẫn bảo toàn độ chính xác tái tạo ảnh, đồng thời đảm bảo khả năng hiện thực hóa trên phần cứng hệ thống?
  • H1: Việc tích hợp kỹ thuật kết hợp hai tần số (Dual Frequency – DF) với điểm chuyển giao vòng lặp tối ưu sẽ mở rộng bán kính hội tụ ở tần số thấp và tối ưu hóa độ phân giải ở tần số cao, triệt tiêu hiện tượng tạo tác (artifacts).
  • H2: Việc ứng dụng kỹ thuật lấy mẫu nén giả ngẫu nhiên (Deterministic Compressed Sensing – DCS) dựa trên chuỗi hỗn loạn Logistic sẽ giảm mạnh số phép đo $N_t \times N_r$ so với chuẩn Nyquist mà không đòi hỏi phần cứng tạo số ngẫu nhiên phức tạp (HRNG/PRNG).

Phạm vi nghiên cứu tập trung vào mô hình hóa 2D sóng tán xạ với chùm Bessel bậc 0, trường tới đơn sắc, môi trường nền đồng nhất chứa đối tượng u có độ tương phản tốc độ âm $\Delta c$ từ 2% đến 5%, kiểm chứng thực nghiệm trên phantom hình cầu đường kính $0,3 \pm 0,05\text{ mm}$ ở tần số $0,64\text{ MHz}$.

Literature Review và Positioning

Lịch sử phát triển của siêu âm cắt lớp phân tách thành hai trường phái lý thuyết chính:

Trường phái thứ nhất dựa trên lý thuyết chiếu tia thẳng và hiệu chỉnh khúc xạ (Ray-based Tomography), tiêu biểu với các hệ thống thực nghiệm như CURE (Computerized Ultrasound Risk Evaluation) do Viện Ung thư Karmanos phát triển và hệ thống HUTT (High-resolution Ultrasonic Transmission Tomography) của Đại học Nam California. Mặc dù đạt tốc độ tái tạo nhanh, mô hình chiếu tia chỉ khôi phục chính xác các cấu trúc có kích thước lớn hơn nhiều bước sóng (từ 2 đến 5 bước sóng $\lambda$). Khi kích thước khối u nhỏ hơn bước sóng, hiệu ứng nhiễu xạ và tán xạ đa hướng chiếm ưu thế, khiến lý thuyết tia hoàn toàn mất hiệu lực.

Trường phái thứ hai dựa trên lý thuyết tán xạ sóng (Wave-based Inverse Scattering), giải trực tiếp phương trình sóng Helmholtz thông qua xấp xỉ Born bậc 1 hoặc xấp xỉ Rytov bậc 1. Luận án chỉ rõ sự tranh luận học thuật sâu sắc giữa hai phép xấp xỉ này. Xấp xỉ Rytov bậc 1 không giới hạn kích thước đối tượng nhưng mắc phải hiện tượng gập pha (phase wrapping) nghiêm trọng khi ước tính pha tổng từ áp suất tán xạ, gây ra sai số bùng nổ khi nghịch đảo ma trận. Ngược lại, xấp xỉ Born bậc 1 tuyến tính hóa bài toán ổn định hơn nhưng bị ràng buộc bởi điều kiện độ lệch pha cực đại $|\Delta \phi| < \pi$. Trích dẫn trực tiếp từ luận án chỉ ra: "Hạn chế này phù hợp với sự kỳ vọng của phương pháp xấp xỉ Born bậc 1 là phương pháp này sẽ bị mất hiệu quả khi $|\Delta \phi| > \pi$."

Trong nhánh nghiên cứu DBIM (Distorted Born Iterative Method), các công trình quốc tế như của Haddadin & Ebbini (1998) hay Lavarello & Oelze (2009) đã áp dụng cập nhật hàm Green phi tuyến theo từng bước lặp để tăng tốc hội tụ. Tuy nhiên, DBIM truyền thống rất nhạy cảm với nhiễu bậc cao và đòi hỏi số lượng đầu dò cực lớn. Khi so sánh với các hệ thống thương mại hóa và tiền lâm sàng hàng đầu thế giới:

  1. Hệ thống SoftVue (Delphinus Medical Technologies, Hoa Kỳ): Sử dụng mảng đầu dò vòng tròn gồm 2048 biến tử, 512 kênh phát và 512 kênh thu, hoạt động ở tần số $3\text{ MHz}$, độ phân giải ảnh B-mode $0,7\text{ mm}$, độ dày lát cắt $2,5\text{ mm}$, đường kính quét tối đa $22\text{ cm}$, thời gian quét mỗi lát cắt mất 15 giây.
  2. Hệ thống 3D USCT (Viện Công nghệ Karlsruhe – KIT, Đức): Sử dụng khẩu độ trụ với hàng nghìn đầu dò siêu âm 3D tĩnh, đối mặt với khối lượng dữ liệu khổng lồ hàng chục Gigabytes cho mỗi lần quét.

Nghiên cứu của NCS Trần Quang Huy đã định vị chính xác khoảng trống: Trong khi thế giới giải quyết bài toán bằng cách tăng số lượng phần cứng đầu dò (như SoftVue dùng 2048 biến tử), luận án đi sâu vào tối ưu hóa thuật toán xử lý tín hiệu bậc cao: phát triển giải thuật DF-DBIM xác định chính xác thời điểm chuyển giao tần số để phá vỡ giới hạn Born, và giải thuật DCS-DBIM ứng dụng chuỗi hỗn loạn để giảm thiểu số lượng đầu dò vật lý, tạo ra bước đột phá về tính khả thi phần cứng.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng nền tảng lý thuyết tán xạ ngược phi tuyến bằng cách tích hợp phương trình vi phân sóng Helmholtz không đồng nhất với lý thuyết điều khiển tối ưu và lý thuyết lấy mẫu nén:

Phương trình sóng mô tả sự lan truyền áp suất âm trong mô sinh học có dạng: $$(\nabla^2 + k_0^2) p(\vec{r}) = -O(\vec{r}) p(\vec{r})$$ Trong đó $k_0 = 2\pi f / c_0$ là số sóng trong môi trường nền đồng nhất (nước, $c_0 = 1484\text{ m/s}$), $k(\vec{r}) = 2\pi f / c(\vec{r})$ là số sóng trong môi trường mô, và hàm mục tiêu lý tưởng $O(\vec{r})$ được định nghĩa: $$O(\vec{r}) = k_0^2 \left( \frac{c_0^2}{c^2(\vec{r})} - 1 \right)$$

Thông qua biểu diễn tích phân sử dụng hàm Green hai chiều $G_0(|\vec{r} - \vec{r}'|) = -\frac{j}{4} H_0^{(1)}(k_0 |\vec{r} - \vec{r}'|)$ (với $H_0^{(1)}$ là hàm Hankel loại 1 bậc 0), trường tán xạ được mô hình hóa chính xác. Rời rạc hóa bằng phương pháp Moment (MoM) biến đổi bài toán liên tục thành hệ phương trình đại số ma trận: $$\bar{p} = (\bar{I} - \bar{C} \text{diag}(\bar{O}))^{-1} \bar{p}_{inc}$$ $$\bar{p}^{sc} = \bar{B} \text{diag}(\bar{O}) \bar{p}$$ Trong đó $\bar{B}$ và $\bar{C}$ là các ma trận hàm Green mô tả tương tác từ các phần tử lưới tới đầu dò thu và tương tác nội tại giữa các phần tử lưới.

Mô hình lý thuyết của luận án xác lập 2 mệnh đề then chốt:

  • Mệnh đề 1 (Chuyển giao đa tần DF-DBIM): Cho quá trình lặp vi phân Born với hai tần số $f_1 < f_2$, tồn tại một số vòng lặp tối ưu $N_{f1}^$ tại tần số $f_1$ sao cho ảnh khôi phục $\bar{O}^{(N_{f1}^)}$ nằm trong bán kính hội tụ của phép xấp xỉ Born bậc 1 ở tần số $f_2$, tức $|\Delta \phi(f_2)| < \pi$, giúp quá trình lặp tiếp theo tại $f_2$ đạt sai số thặng dư tương đối nhỏ nhất mà không bị phân kỳ.
  • Mệnh đề 2 (Đo nén xác định DCS-DBIM): Bằng việc phân bố vị trí máy thu theo chuỗi số giả ngẫu nhiên sinh từ ánh xạ Logistic $x_{n+1} = \mu x_n (1 - x_n)$ với tham số điều khiển $\mu = 4$ và $x_0 \in (0, 1)$, ma trận hệ thống $\bar{M}_{DCS}$ thỏa mãn tính chất đẳng cự cục bộ (Restricted Isometry Property – RIP) với xác suất tương đương phân bố ngẫu nhiên Gauss, cho phép giải bài toán tối thiểu hóa chuẩn $l_1$ ($l_1$-regularized least squares problem – $l_1$-LSP) với số lượng phép đo $M \ll N_t \times N_r$.

Khung phân tích độc đáo

Khung phân tích của luận án là sự hợp nhất của 4 trụ cột lý thuyết: (1) Lý thuyết sóng âm Helmholtz và hàm Green, (2) Xấp xỉ Born và phương pháp Moment MoM, (3) Chuẩn tắc hóa Tikhonov kết hợp tối ưu hóa $l_1$-LSP, và (4) Lý thuyết hệ động lực phi tuyến (bản đồ hỗn loạn).

Điều kiện biên của khung phân tích được xác định rõ ràng: đối tượng quan sát là các tổn thương dạng khối có đường kính $d \le \lambda$ đến $2\lambda$, độ tương phản tốc độ âm yếu đến trung bình ($\Delta c/c_0 \le 5%$), môi trường truyền sóng không xét đến tán xạ dị hướng bậc cao của xương hoặc khí (vốn có tốc độ âm $4090\text{ m/s}$ hoặc $330\text{ m/s}$).

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (Positivism), kết hợp chặt chẽ giữa suy diễn toán học hình thức, mô phỏng số học chuyên sâu trên máy tính và đối soát thực nghiệm vật lý. Thiết kế đa tầng bao gồm:

  1. Tầng 1 (Mô hình toán học): Xây dựng hệ phương trình tích phân tán xạ và các toán tử ma trận $\bar{B}, \bar{C}, \bar{M}$.
  2. Tầng 2 (Mô phỏng số học): Thiết lập các kịch bản kiểm thử trên MATLAB với các kích thước ma trận lưới $N = 16, 21, 22, 32, 64$ điểm ảnh; cấu hình thu phát từ 16 đến 22 máy phát/thu.
  3. Tầng 3 (Kiểm chứng thực nghiệm): Thu thập dữ liệu áp suất tán xạ thực tế từ hệ đo thí nghiệm tại bể nước để hiệu chuẩn mô hình mô phỏng.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được chuẩn hóa nghiêm ngặt:

  • Nguồn phát sóng tới: Sử dụng chùm Bessel bậc 0 đơn sắc: $$p_{inc}(\vec{r}) = J_0(k_0 |\vec{r} - \vec{r}t|) = \sum{m=0}^\infty \frac{(-1)^m}{m! \Gamma(m+1)} \left(\frac{k_0 |\vec{r} - \vec{r}_t|}{2}\right)^{2m}$$
  • Bộ thu tín hiệu tán xạ: Bố trí trên đường tròn bao quanh đối tượng. Quá trình phát – thu đơn phần tử tuần tự: tại mỗi thời điểm chỉ có 1 đầu dò phát, toàn bộ $N_r$ đầu dò thu ghi nhận đồng thời tín hiệu áp suất tổng $p(\vec{r}r)$, sau đó trừ đi áp suất trường tới $p{inc}(\vec{r}_r)$ để trích xuất tín hiệu tán xạ $p^{sc}(\vec{r}_r)$.
  • Triệt tiêu nhiễu và kiểm chuẩn: Dữ liệu thực nghiệm được đo đạc trong bể nước ở nhiệt độ ổn định $18,1^\circ\text{C}$, sử dụng bộ lọc trung bình (moving average filter) để khử nhiễu ngẫu nhiên trước khi đưa vào thuật toán nghịch đảo.

Data và phân tích

Dữ liệu mô phỏng và thực nghiệm được xử lý qua hai giải thuật cải tiến cốt lõi:

Các tham số thống kê mô phỏng:

  • Kịch bản DF-DBIM: Tần số $f_1 = 1\text{ MHz}$, $f_2 = 2\text{ MHz}$, số máy phát/thu $N_t = N_r = 22$, số điểm chia lưới $N = 22$ ($N^2 = 484$ pixel), tổng số vòng lặp $I_{total} = 8$. Khảo sát $N_{f1}$ biến thiên từ 1 đến 7.
  • Kịch bản DCS-DBIM: Lưới $N = 21$ ($N^2 = 441$ pixel), so sánh cấu hình đều truyền thống ($N_t = N_r = 20 \Rightarrow 400$ phép đo; $N_t = N_r = 22 \Rightarrow 484$ phép đo; $N_t = N_r = 30 \Rightarrow 900$ phép đo) với cấu hình DCS đề xuất ($N_t = N_r = 16 \Rightarrow 256$ phép đo; $N_t = N_r = 20 \Rightarrow 400$ phép đo).

Chỉ số đánh giá định lượng là Sai số thặng dư tương đối (Relative Residual Error – RRE): $$\text{RRE} = \frac{|\bar{O}{thực} - \bar{O}{khôi_phục}|2}{|\bar{O}{thực}|_2}$$

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Khám phá quy luật điểm chuyển giao tối ưu trong giải thuật DF-DBIM: Nghiên cứu chứng minh rằng việc kết hợp tần số không tự động mang lại kết quả tốt nếu không chọn đúng thời điểm chuyển giao $N_{f1}$. Khi khảo sát 4 kịch bản mô phỏng với các tỷ lệ số phép đo trên số biến khác nhau, kết quả chỉ ra:
    • Nếu chuyển giao quá sớm ($N_{f1} = 1$), thuật toán tại tần số cao $f_2 = 2\text{ MHz}$ bị phân kỳ hoặc tạo nhiều nhiễu đốm giả do ước lượng ban đầu chưa thỏa mãn điều kiện xấp xỉ Born.
    • Nếu chuyển giao quá muộn ($N_{f1} = 6$ hoặc $7$), sai số bị bão hòa ở mức cao của tần số thấp $f_1 = 1\text{ MHz}$, không tận dụng được số vòng lặp của tần số cao để làm mịn ảnh.
    • Điểm chuyển giao tối ưu đạt được tại $N_{f1} = 2$ hoặc $N_{f1} = 3$ (trên tổng số 8 vòng lặp), giúp sai số chuẩn hóa RRE giảm xuống mức thấp kỷ lục (dưới 0,15), vượt trội hoàn toàn so với việc chỉ chạy đơn tần $1\text{ MHz}$ hoặc $2\text{ MHz}$.
  Sai số RRE
     1       2       3       4       5       6       7       8
  1. Cắt giảm từ 47% đến 55% số lượng phép đo nhờ giải thuật DCS-DBIM: Trong cấu hình đo truyền thống, để tái tạo lưới $21 \times 21$, hệ thống cần $N_t = N_r = 22$ (484 phép đo). Phương pháp DCS-DBIM đề xuất chỉ cần $N_t = N_r = 16$ (256 phép đo) bố trí theo luật Logistic map đã đạt được chất lượng tương đương, thậm chí lỗi chuẩn hóa sau 8 vòng lặp còn nhỏ hơn cấu hình đo đều 484 phép đo. Khi so sánh ở quy mô 900 phép đo ($N_t = N_r = 30$), DCS-DBIM chỉ cần 400 phép đo ($N_t = N_r = 20$) để đạt cùng độ chính xác.
Phương pháp Cấu hình ($N_t \times N_r$) Tổng số phép đo Thời gian tạo ảnh (s) Sai số chuẩn hóa (RRE)
DBIM truyền thống $22 \times 22$ 484 Lớn ($O(N_t N_r N^2)$) 0,215
DCS-DBIM đề xuất $16 \times 16$ 256 (Giảm 47,1%) Giảm ~50% 0,184 (Tốt hơn)
DBIM truyền thống $30 \times 30$ 900 Rất lớn 0,162
DCS-DBIM đề xuất $20 \times 20$ 400 (Giảm 55,6%) Giảm ~58% 0,151 (Tốt hơn)

  1. Hiện thực hóa phần cứng lấy mẫu nén không cần mạch tạo số ngẫu nhiên: Phát hiện mang tính đột phá kỹ thuật là việc chứng minh ma trận đo giả ngẫu nhiên sinh bởi hàm phi tuyến Logistic hoàn toàn tương thích với phần cứng mảng đầu dò cố định. Kỹ sư chế tạo chỉ cần lập trình sẵn bảng vị trí đầu dò theo thuật toán đệ quy số học đơn giản mà không cần tích hợp bộ tạo số ngẫu nhiên thực HRNG cồng kềnh, đắt tiền và thiếu ổn định.

  2. Kiểm chứng thực nghiệm thành công với vi tổn thương $0,3\text{ mm}$: Dữ liệu thực nghiệm đo áp suất tán xạ của quả bóng đường kính $0,3\text{ mm}$ (ở bước sóng $\lambda \approx 2,3\text{ mm}$ tại $0,64\text{ MHz}$) cho thấy sự trùng khớp dạng sóng tán xạ giữa mô phỏng và thực tế sau khi qua bộ lọc trung bình, xác nhận mô hình toán học phản ánh chính xác bản chất vật lý của môi trường tán xạ.

Implications đa chiều

  • Về mặt học thuật: Mở ra hướng tiếp cận mới trong xử lý bài toán phi tuyến ngược bằng cách kết hợp tối ưu hóa đa tần và lý thuyết nén tín hiệu xác định, cung cấp công cụ toán học giải quyết các bài toán tán xạ điện từ và địa chấn.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn mực từ tính toán giải tích hàm Green, rời rạc hóa MoM, mô phỏng tham số đến hiệu chuẩn thực nghiệm bể nước.
  • Về mặt kỹ thuật y sinh: Cung cấp nền tảng thuật toán cốt lõi để phát triển các máy siêu âm cắt lớp thế hệ mới với số lượng kênh phần cứng giảm một nửa, hạ giá thành sản xuất thiết bị từ hàng trăm nghìn USD xuống mức bình dân.
  • Về mặt y tế công cộng: Tạo tiền đề cho các chương trình tầm soát ung thư vú đại trà an toàn (không nhiễm xạ tia X), không gây đau đớn (không ép ngực như mammography), phát hiện u ác tính ở giai đoạn khởi phát (< 5 mm), nâng tỷ lệ sống sót sau 5 năm của bệnh nhân lên trên 90%.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn khoa học:

  1. Giới hạn không gian 2D: Nghiên cứu giới hạn ở bài toán lát cắt phẳng 2D, giả định chùm sóng Bessel hình trụ. Trong thực tế lâm sàng, trường sóng âm trong cơ thể người là không gian 3D với hiệu ứng tán xạ ngoài mặt phẳng (out-of-plane scattering).
  2. Đơn tham số tốc độ truyền âm: Mô hình chỉ tập trung tái tạo độ tương phản tốc độ âm $\Delta c$, bỏ qua sự biến thiên đồng thời của độ suy hao âm $\alpha(\vec{r})$ và mật độ khối $\rho(\vec{r})$.
  3. Môi trường nền lý tưởng hóa: Mô phỏng giả định môi trường xung quanh u là đồng nhất, chưa xét đến tính bất đồng nhất phức tạp của các bè sợi và tuyến sữa trong vú thật.

Chương trình nghiên cứu tiếp nối 5 - 10 năm tới:

  • Mở rộng thuật toán DF-DBIM và DCS-DBIM sang không gian 3D đầy đủ (Full 3D Vector Wave Scattering).
  • Phát triển giải thuật tái tạo đồng thời đa tham số (Multi-parameter Reconstruction: Sound Speed, Attenuation, and Density).
  • Tích hợp mạng nơ-ron tích chập sâu (Deep CNNs) để học trước ma trận khởi tạo $\bar{O}^{(0)}$, rút ngắn số vòng lặp DBIM từ 8 vòng xuống còn 1 - 2 vòng lặp thời gian thực.
  • Thử nghiệm tiền lâm sàng trên phantom mô vú sinh học đa lớp (Tissue-mimicking Breast Phantoms) và thử nghiệm lâm sàng in-vivo.

Tác động và ảnh hưởng

  • Tác động học thuật: Đóng góp 02 bài báo ISI/Scopus và nhiều báo cáo hội nghị quốc tế chuyên ngành kỹ thuật y sinh (IEEE), mở ra các trích dẫn quan trọng trong lĩnh vực xử lý tín hiệu tán xạ ngược.
  • Tác động công nghiệp: Cung cấp giải pháp thiết kế mảng đầu dò tiết kiệm cho các công ty sản xuất thiết bị y tế trong nước và quốc tế, giảm chi phí chế tạo hệ thống mảng thu phát siêu âm cắt lớp.
  • Tác động kinh tế – xã hội: Giúp giảm gánh nặng chi phí y tế cho bệnh nhân ung thư vú tại Việt Nam; giảm tải cho các bệnh viện tuyến trên nhờ thiết bị sàng lọc không xâm lấn, giá thành thấp có thể triển khai tại các trung tâm y tế cơ sở.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học giả: Tiếp cận khung toán học chuẩn xác về phương pháp Moment, xấp xỉ Born và kỹ thuật lấy mẫu nén áp dụng cho bài toán sóng Helmholtz.
  • Kỹ sư R&D thiết bị y tế: Sử dụng trực tiếp thuật toán DCS-DBIM và sơ đồ chuỗi Logistic để thiết kế bố trí phần cứng đầu dò tối ưu trên máy chụp cắt lớp.
  • Bác sĩ chẩn đoán hình ảnh: Có thêm cơ sở tin cậy từ hình ảnh định lượng tốc độ truyền âm để chẩn đoán phân biệt u lành và u ác tính mà không cần sinh thiết xâm lấn.
  • Nhà hoạch định chính sách y tế: Có luận cứ khoa học để đầu tư công nghệ tầm soát sớm bệnh không lây nhiễm, nâng cao tuổi thọ trung bình và chất lượng dân số.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng Lý thuyết tán xạ ngược lặp Born (DBIM) thông qua việc tích hợp Cơ chế chuyển giao tần số tối ưu ($N_{f1}^*$) và Lý thuyết lấy mẫu nén xác định (DCS). Luận án đã chứng minh toán học và số học rằng việc chuyển giao tần số tại đúng vòng lặp thứ 2 hoặc 3 sẽ mở rộng miền hội tụ của xấp xỉ Born bậc 1 ($|\Delta \phi| < \pi$), khắc phục triệt để hiện tượng phân kỳ và nhiễu tạo tác mà lý thuyết DBIM đơn tần cổ điển của Chew & Wang (1990) mắc phải.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó? Trả lời:

  • So với nghiên cứu của Slaney et al. (1984) và Lavarello et al. (2009) (vốn chỉ khảo sát DBIM đơn tần hoặc đa tần liên tục đòi hỏi quét toàn bộ phổ tần số gây bùng nổ dữ liệu), luận án đề xuất cấu trúc DF-DBIM hai tần số với quy tắc ngắt vòng lặp tối ưu, giảm hơn 60% dung lượng bộ nhớ lưu trữ dữ liệu đo.
  • So với nghiên cứu của Provost & Lesage (2009) về lấy mẫu nén ngẫu nhiên trong siêu âm (đòi hỏi phân bố đầu dò ngẫu nhiên Gaussian/Bernoulli không thể chế tạo cơ khí chính xác), luận án đổi mới bằng phương pháp DCS sử dụng ánh xạ hỗn loạn Logistic, cố định hóa tọa độ đầu dò trên vành khuyên mà vẫn đảm bảo tính chất RIP để giải bài toán $l_1$-LSP.

3. Phát hiện bất ngờ nhất (counter-intuitive) có số liệu chứng minh là gì? Trả lời: Phát hiện bất ngờ nhất là Cấu hình đo ít đầu dò hơn lại cho sai số tái tạo nhỏ hơn cấu hình đo nhiều đầu dò khi ứng dụng DCS. Cụ thể, cấu hình DCS-DBIM với 256 phép đo ($N_t = N_r = 16$) cho sai số thặng dư tương đối $\text{RRE} = 0,184$, thấp hơn đáng kể so với cấu hình đo đều truyền thống 484 phép đo ($N_t = N_r = 22, \text{RRE} = 0,215$). Điều này giải thích bởi việc phân bố máy thu theo quy luật hỗn loạn Logistic đã triệt tiêu hiện tượng tương quan đồng pha (coherence) giữa các hàng của ma trận Green, giảm hiện tượng ma trận suy biến trong chuẩn tắc Tikhonov.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (replication protocol) không? Trả lời: Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình tái lập: (1) Bộ tham số mô phỏng gồm bước sóng, bán kính đối tượng, độ tương phản $\Delta c = 2%$, kích thước lưới $N$; (2) Công thức giải tích rời rạc hóa ma trận $\bar{B}, \bar{C}$ qua hàm Hankel và Bessel; (3) Mã giả chi tiết cho 4 giải thuật cốt lõi (Giải thuật DBIM chuẩn, Giải thuật DF-DBIM, Khảo sát $N_{f1}$, Giải thuật DCS-DBIM); (4) Thông số thực nghiệm bể nước ($T = 18,1^\circ\text{C}$, khoảng cách phát $104\text{ mm}$, khoảng cách thu $91\text{ mm}$, tần số $0,64\text{ MHz}$, bước quay $1^\circ$).

5. Lộ trình nghiên cứu 10 năm được phác thảo như thế nào? Trả lời: Lộ trình bao gồm: Giai đoạn 1 (1 - 3 năm): Nâng cấp thuật toán sang hệ tọa độ cầu 3D Vector DBIM; Giai đoạn 2 (4 - 6 năm): Tích hợp Deep Learning để giải bài toán ngược thời gian thực và chế tạo module đầu dò DCS 256 kênh; Giai đoạn 3 (7 - 10 năm): Hoàn thiện hệ thống máy chụp siêu âm cắt lớp vú 3D hoàn chỉnh và tiến hành thử nghiệm lâm sàng đa trung tâm tại các bệnh viện ung bướu.

Kết luận

  1. Xây dựng thành công mô hình toán học và giải thuật siêu âm cắt lớp tán xạ ngược 2D dựa trên phương trình sóng Helmholtz toàn phần và phương pháp Moment, được kiểm chứng nhất quán giữa mô phỏng số học và thực nghiệm vật lý trên đối tượng vi mô $0,3\text{ mm}$.
  2. Đề xuất giải thuật DF-DBIM đột phá, giải quyết triệt để sự đánh đổi giữa tốc độ hội tụ và độ phân giải không gian bằng cách xác định chính xác điểm chuyển giao vòng lặp tối ưu $N_{f1} = 2$ hoặc $3$, triệt tiêu sai số phân kỳ của xấp xỉ Born bậc 1.
  3. Phát triển thành công giải thuật DCS-DBIM ứng dụng chuỗi hỗn loạn Logistic, cắt giảm từ 47% đến 55% số lượng phép đo và thời gian tính toán mà vẫn đảm bảo độ chính xác tái tạo vượt trội so với cấu hình đo đều truyền thống.
  4. Giải quyết bài toán thắt nút cổ chai về phần cứng của kỹ thuật lấy mẫu nén trong siêu âm cắt lớp, mở đường cho việc chế tạo mảng đầu dò cố định giá thành thấp mà không cần mạch tạo số ngẫu nhiên phức tạp.
  5. Mở ra 3 hướng nghiên cứu học thuật mới: Tái tạo ảnh tán xạ 3D đa tham số (tốc độ âm, suy hao, mật độ), kết hợp học sâu (Deep Learning) tăng tốc nghịch đảo ma trận Green, và kỹ thuật chẩn đoán sớm ung thư vú không xâm lấn, đóng góp giá trị khoa học và thực tiễn to lớn cho nền y tế nước nhà và quốc tế.