Tổng quan về luận án

Công nghệ tạo ảnh y sinh phi xâm lấn đóng vai trò mang tính sống còn trong chẩn đoán lâm sàng hiện đại, đặc biệt là tầm soát và điều trị sớm các bệnh lý ác tính. Ung thư vú là một trong những nguyên nhân gây tử vong hàng đầu ở phụ nữ trên quy mô toàn cầu. Theo số liệu dịch tễ học được luận án trích dẫn từ y văn quốc tế, "Trên thế giới mỗi năm có khoảng 1,4 triệu người mới mắc bệnh ung thư vú và 458.000 người tử vong. Ở Việt Nam, mỗi năm có khoảng 7.000 người mới mắc bệnh, trong đó tỷ lệ tử vong khoảng 35%, cao hơn các nước phát triển, bệnh có xu hướng trẻ hóa và tăng dần qua các năm". Khi khối u được phát hiện ở giai đoạn khởi phát (đường kính giải phẫu dưới 5 mm), tỷ lệ điều trị thành công đạt tới 90%. Tuy nhiên, các kỹ thuật chẩn đoán truyền thống bộc lộ những rào cản căn bản: chụp X-quang tuyến vú (Mammography) sử dụng bức xạ ion hóa và bị giảm độ nhạy trầm trọng đối với mô vú dày đặc ở phụ nữ dưới 50 tuổi; siêu âm phản xạ quy ước (B-mode) chỉ mang tính định tính, phụ thuộc chủ quan vào người thao tác và khó phát hiện vi tổn thương do hạn chế tương phản âm học.

Trong bối cảnh đó, kỹ thuật Siêu âm cắt lớp (Ultrasound Tomography - UST) dựa trên lý thuyết tán xạ ngược nổi lên như một hướng tiếp cận đột phá, cho phép tái tạo định lượng các thuộc tính vật lý của mô sinh học như tốc độ truyền âm ($\Delta c$), độ suy hao ($\alpha$) và mật độ ($\rho$). Luận án tiến sĩ kỹ thuật điện tử của NCS. Trần Quang Huy, dưới sự hướng dẫn khoa học của PGS. Trần Đức Tân tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội (2019), tập trung giải quyết khoảng trống nghiên cứu cốt lõi: "Phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM) có độ phức tạp tính toán cao dẫn đến thời gian tạo ảnh lâu và nhạy cảm với nhiễu, khiến số lượng thiết bị siêu âm cắt lớp thương mại hóa trên thế giới hiện nay còn rất hạn chế".

Mục tiêu và giả thuyết nghiên cứu của luận án được cấu trúc tường minh qua hệ thống câu hỏi và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để giải quyết mâu thuẫn giữa độ ổn định hội tụ và độ phân giải không gian trong bài toán nghịch đảo phi tuyến khi độ tương phản vận tốc âm của mô thay đổi?
  • Giả thuyết 1 (H1): Việc tích hợp kỹ thuật kết hợp hai tần số (Dual-Frequency - DF) vào giải thuật DBIM với việc xác định chính xác thời điểm chuyển tiếp tần số sẽ triệt tiêu khả năng rơi vào cực tiểu địa phương ở tần số cao mà vẫn đảm bảo độ phân giải cấp độ mô học.
  • Câu hỏi nghiên cứu 2 (RQ2): Có thể giảm thiểu triệt để số lượng phần tử đầu dò và số phép đo thu nhận mà không làm suy giảm chất lượng tái tạo ảnh định lượng, đồng thời khắc phục rào cản hiện thực hóa phần cứng của kỹ thuật lấy mẫu nén ngẫu nhiên thuần túy?
  • Giả thuyết 2 (H2): Áp 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 ánh xạ phi tuyến hỗn loạn (Logistic Map) kết hợp với tối ưu hóa chuẩn $l_1$ (regularized least squares problem - $l_1$-LSP) sẽ giảm tối thiểu $47%$ số lượng phép đo thực tế so với cấu hình lấy mẫu đều truyền thống mà vẫn bảo toàn độ chính xác tái tạo cấu trúc khối u.

Phạm vi nghiên cứu bao gồm xây dựng cơ sở lý thuyết toán-lý sóng âm, mô phỏng số học chuyên sâu trên các lưới phân giải $N = 16 \times 16, 32 \times 32, 64 \times 64$ điểm ảnh, và kiểm chứng thực nghiệm trên mô hình vật lý chuẩn với nguồn phát $f = 0.64\text{ MHz}$, đối tượng tán xạ hình trụ kích thước micro $0.3 \pm 0.05\text{ mm}$ đặt trong môi trường nước $18.1^\circ\text{C}$.

Literature Review và Positioning

Lịch sử phát triển của siêu âm định lượng ghi nhận sự phân nhánh thành hai trường phái lý thuyết chính: Siêu âm cắt lớp dựa trên lý thuyết chiếu tia (Ray-based Tomography) và Siêu âm cắt lớp dựa trên lý thuyết tán xạ sóng (Wave-scattering Tomography).

Trường phái chiếu tia (Ray Theory), tiêu biểu với các hệ thống tiên phong từ cuối thập niên 1970 và đầu thập niên 1980 như nghiên cứu của Greenleaf và cộng sự (1981), Carson và cộng sự (1981), Schreiman và cộng sự (1984), tiền đề cho các hệ thống lâm sàng hiện đại như CURE (Computerized Ultrasound Risk Evaluation) phát triển bởi Viện Ung thư Karmanos (Duric et al., 2007) và HUTT (High-resolution Ultrasonic Transmission Tomography) tại Đại học Southern California (USC). Mặc dù các kỹ thuật hiệu chỉnh khúc xạ (refraction correction) đã được tích hợp trong hệ thống CURE của Duric và cộng sự, nhưng về mặt bản chất vật lý, phép xấp xỉ tia thẳng chỉ hợp thức khi kích thước dị vật lớn hơn nhiều so với bước sóng ($\ge 2\lambda - 5\lambda$). Khi kích thước mô bệnh học tương đương hoặc nhỏ hơn bước sóng ($\le \lambda$), hiện tượng nhiễu xạ và tán xạ đa hướng trở nên chiếm ưu thế, khiến lý thuyết chiếu tia hoàn toàn mất khả năng khôi phục định lượng chính xác ranh giới và mật độ tổn thương.

Trường phái lý thuyết tán xạ giải bài toán sóng toàn phần (Full-wave Inverse Scattering) thông qua phương trình sóng Helmholtz vi phân bậc hai. Trong miền tán xạ phi tuyến, hai phương pháp tuyến tính hóa cổ điển thường được áp dụng là xấp xỉ Born bậc 1 (Born, 1926) và xấp xỉ Rytov bậc 1 (Rytov, 1937). Tranh luận học thuật kéo dài giữa hai hướng tiếp cận này: xấp xỉ Rytov giả định biến thiên pha mượt mà nên áp dụng được cho đối tượng kích thước lớn, nhưng lại vấp phải hiện tượng gập pha (Phase Wrapping) gây sai số bùng nổ khi nghịch đảo trường sóng (Slaney et al., 1984). Ngược lại, xấp xỉ Born bậc 1 giải quyết trực tiếp biên độ trường tán xạ nhưng bị giới hạn nghiêm ngặt bởi điều kiện lệch pha cực đại:

$$|\Delta \phi| = \left| k_0 \int \Delta n(r) , dr \right| < \pi$$

Khi độ tương phản âm học vượt ngưỡng này, phương pháp Born quy ước mất tính hội tụ và tạo ra các tạo tác (nhiễu đốm - artifacts) nghiêm trọng tại tâm đối tượng.

                    BẢN ĐỒ TIẾN HÓA LÝ THUYẾT SIÊU ÂM CẮT LỚP (UST)
                    
  Lý thuyết chiếu tia (Ray Theory)            Lý thuyết tán xạ ngược (Wave Scattering)
  (Greenleaf 1981, Duric/CURE 2007)          (Helmholtz Equation & Inverse Problem)
  Giới hạn tán xạ: Chỉ áp dụng                Xấp xỉ Rytov              Xấp xỉ Born bậc 1
  khi kích thước u > 2λ - 5λ               (Bị lỗi Gập pha -          (Bị giới hạn lệch pha
                       VỊ TRÍ ĐÓNG GÓP CỦA LUẬN ÁN (2019)
   Kỹ thuật DF-DBIM đề xuất                   Kỹ thuật DCS-DBIM đề xuất
   - Kết hợp 2 tần số tối ưu                  - Lấy mẫu nén giả ngẫu nhiên (Logistic Map)
   - Triệt tiêu phân kỳ cực tiểu địa phương    - Giảm số kênh phần cứng, tối ưu l₁-LSP

Để khắc phục tính phi tuyến mạnh, Chew và Wang (1990) đã đề xuất Phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM). DBIM cập nhật hàm Green phi tuyến theo từng bước lặp, cho tốc độ hội tụ nhanh hơn phương pháp lặp Born (BIM), nhưng lại có độ phức tạp tính toán khổng lồ $\mathcal{O}(N_t N_r N^2)$ và cực kỳ nhạy cảm với nhiễu đo lường. Trên bình diện quốc tế, các hệ thống thương mại tiên tiến nhất như SoftVue của Delphinus Medical Technologies (2017) phải sử dụng cấu hình phần cứng đồ sộ với "Số biến tử: 2048; Số kênh nhận: 512; Số kênh phát: 512; Tần số hoạt động: 3MHz; Độ phân giải dữ liệu: 14 bít; Độ phân giải ảnh (B-mode): 0.7mm; Độ dày lát cắt: 2,5mm; Đường kính vú tối đa: 22cm; Thời gian thu thập trên mỗi lát cắt: 15s".

Luận án định vị chính xác tại điểm giao thoa giữa xử lý tín hiệu tiên tiến và vật lý sóng âm: khắc phục triệt để độ phức tạp tính toán và độ nhạy nhiễu của DBIM bằng cách tích hợp thuật toán kết hợp tần số tối ưu và lý thuyết lấy mẫu nén tiền định, tạo nền tảng cho các hệ thống phần cứng siêu âm cắt lớp thế hệ mới với chi phí thấp và tốc độ xử lý thời gian thực.

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

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

Luận án tạo ra những bước tiến mang tính nền tảng đối với lý thuyết bài toán nghịch đảo trong điện tử y sinh:

  1. Mở rộng lý thuyết xấp xỉ Born phi tuyến: Thách thức giới hạn hội tụ truyền thống của giải thuật DBIM nguyên bản (Chew & Wang, 1990). Bằng cách thiết lập cơ chế chuyển giao trường tán xạ giữa miền tần số thấp $f_1$ và tần số cao $f_2$, nghiên cứu chứng minh rằng hàm mục tiêu khôi phục tại $f_1$ đóng vai trò làm nghiệm khởi tạo tối ưu (initial guess), đưa bài toán phi tuyến phức tạp vào bán kính hội tụ cục bộ của xấp xỉ Born tại $f_2$.
  2. Thiết lập mô hình giải bài toán ngược kết hợp chuẩn hóa $l_1$ phi lồi: Thay thế phương pháp chuẩn hóa Tikhonov chuẩn $l_2$ truyền thống bằng bài toán bình phương tối thiểu có điều kiện chính quy hóa $l_1$ ($l_1$-regularized Least Squares Problem - $l_1$-LSP). Đóng góp này khai thác triệt để tính thưa (sparsity) của các biến thiên vận tốc âm $\Delta c$ của khối u trên nền mô đồng nhất, cho phép bảo toàn các thành phần tần số cao và đường biên hình học sắc nét của khối u mà không bị làm mịn quá mức (over-smoothing) như chuẩn Tikhonov.
  3. Chuyển dịch mô hình lấy mẫu nén từ ngẫu nhiên sang tiền định (Deterministic Paradigm Shift): Khắc phục khoảng trống lý thuyết của Candès, Wakin (2008) và Donoho (2006) khi ứng dụng trong kỹ thuật siêu âm cắt lớp. Luận án chứng minh chuỗi hỗn loạn Logistic Map thỏa mãn điều kiện đẳng cự cưỡng bức (Restricted Isometry Property - RIP) trong không gian biến đổi, đặt nền móng lý thuyết cho việc thiết kế cấu hình phần cứng cố định không cần mạch sinh số ngẫu nhiên phức tạp.
                      MÔ HÌNH HÌNH HỌC VÀ TƯƠNG TÁC SÓNG TÁN XẠ
                                 (Transmitter / Receiver Array)

                                 Máy thu Nr(k)
                                     [ ■ ]
                                    /     \
                                   /       \  Sóng tán xạ p_sc
                                  /         \
                    Máy phát Nt(i)           \
                                              [          ]
                                            Hàm Green G₀(r, r')
                                            Tương tác Bessel J₀

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

Khung phân tích toán học của luận án được tích hợp chặt chẽ từ 3 trụ cột lý thuyết: Phương trình tích phân sóng Lippmann-Schwinger, Phương pháp Moment rời rạc (Method of Moments - MoM), và Lý thuyết hệ thống phi tuyến hỗn loạn.

Phương trình tích phân biểu diễn trường sóng tổng $p(\vec{r})$ trong môi trường tán xạ không đồng nhất:

$$p(\vec{r}) = p^{inc}(\vec{r}) + \iint_{D} k_0^2 O(\vec{r}') G_0(|\vec{r} - \vec{r}'|) p(\vec{r}') , d\vec{r}'$$

Trong đó $G_0$ là hàm Green biểu diễn bằng hàm Hankel loại 1 bậc 0: $G_0(|\vec{r} - \vec{r}'|) = \frac{j}{4} H_0^{(1)}(k_0 |\vec{r} - \vec{r}'|)$, và hàm mục tiêu biểu thị độ chênh lệch vận tốc âm:

$$O(\vec{r}) = 1 - \frac{c_0^2}{c^2(\vec{r})} + j \frac{2\alpha(\vec{r}) c_0^2}{\omega c(\vec{r})}$$

Rời rạc hóa miền quan tâm thành $N \times N$ phần tử ảnh (pixels) thông qua phép chiếu MoM:

$$\bar{p}^{tot} = (\bar{I} - \bar{C} , \text{diag}(\bar{O}))^{-1} \bar{p}^{inc}$$

$$\bar{p}^{sc} = \bar{B} , \text{diag}(\bar{O}) \bar{p}^{tot}$$

Điều kiện biên xác định: Đối tượng tái tạo giới hạn trong lát cắt 2D hình trụ tròn đường kính $D \le 10\text{ mm}$, độ tương phản vận tốc âm $\Delta c / c_0 \le 10%$, môi trường nền nước đồng nhất vận tốc $c_0 = 1484\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ủ triết lý nghiên cứu thực chứng (Positivism) với quy trình suy diễn logic chặt chẽ: Thiết lập mô hình giải tích $\rightarrow$ Rời rạc hóa số học $\rightarrow$ Xây dựng thuật toán tối ưu $\rightarrow$ Mô phỏng Monte Carlo $\rightarrow$ Kiểm chứng thực nghiệm vật lý thực tế. Thiết kế nghiên cứu kết hợp đa tầng (multi-level design) đánh giá hiệu năng giải thuật trên 4 kịch bản cấu hình máy phát – máy thu khác nhau với các mức độ phân giải ma trận từ thấp ($N=16$) đến cao ($N=64$).

Quy trình nghiên cứu rigorous

Quy trình giải thuật được thực thi tuần tự và nghiêm ngặt qua 2 bước đột phá công nghệ:

1. Thuật toán kết hợp hai tần số DF-DBIM (Dual-Frequency DBIM)

  • Khởi tạo với tần số thấp $f_1 = 0.64\text{ MHz}$ trong $N_{f1}$ vòng lặp đầu tiên để tái tạo cấu trúc vĩ mô của hàm mục tiêu mà không bị phân kỳ: $$\bar{O}^{(n+1)} = \bar{O}^{(n)} + \Delta \bar{O}^{(n)}$$
  • Tính toán sai số thặng dư tương đối (Relative Residual Error - RRE) tại mỗi bước lặp: $$\text{RRE}^{(n)} = \frac{|\bar{p}^{sc}{meas} - \bar{p}^{sc}{est}(\bar{O}^{(n)})|2}{|\bar{p}^{sc}{meas}|_2}$$
  • Cơ chế chuyển giao tối ưu: Xác định tham số chuyển tiếp $\kappa = N_{f1} / N_{total}$. Khi đạt đến ngưỡng chuyển giao $\kappa_{opt}$, ma trận tán xạ thu được tại $f_1$ được nội suy làm điều kiện biên khởi tạo cho tần số cao $f_2 = 1.0\text{ MHz}$ hoặc $2.0\text{ MHz}$ để tinh chỉnh các vi cấu trúc biên cạnh.

2. Thuật toán lấy mẫu nén giả ngẫu nhiên DCS-DBIM (Deterministic CS-DBIM)

  • Thay thế vị trí đầu dò ngẫu nhiên bằng chuỗi hỗn loạn Logistic Map đơn kỳ: $$x_{k+1} = \mu x_k (1 - x_k), \quad x_k \in (0, 1), ; \mu = 4.0$$
  • Ánh xạ các giá trị $x_k$ để xác định góc bố trí của $N_r$ máy thu trên vòng tròn đo quét $360^\circ$: $$\theta_k = 2\pi \cdot x_k$$
  • Giải bài toán tối ưu lồi phi tuyến $l_1$-LSP bằng thuật toán điểm trong (Interior-point Method): $$\min_{\Delta \bar{O}} \left{ \frac{1}{2} |\bar{B}_{DCS} \Delta \bar{O} - \Delta \bar{p}^{sc}|_2^2 + \lambda |\Delta \bar{O}|_1 \right}$$ Trong đó $\lambda > 0$ là tham số điều hòa kiểm soát độ thưa của nghiệm.
                    QUY TRÌNH THỰC THI THUẬT TOÁN DCS-DBIM ĐỀ XUẤT

Data và phân tích

Kịch bản phân tích Số lượng máy phát ($N_t$) Số lượng máy thu ($N_r$) Tổng số phép đo ($N_t \times N_r$) Tỷ số nén mẫu ($r$) Độ phân giải lưới ($N \times N$)
Kịch bản 1 (Đo đều) 22 22 484 1.00 (Không nén) $22 \times 22$ (484 ẩn số)
Kịch bản 2 (DCS đề xuất) 16 16 256 0.53 (Nén 47.1%) $22 \times 22$ (484 ẩn số)
Kịch bản 3 (Đo đều cao) 30 30 900 1.00 (Không nén) $21 \times 21$ (441 ẩn số)
Kịch bản 4 (DCS tối ưu) 20 20 400 0.44 (Nén 55.6%) $21 \times 21$ (441 ẩn số)

Hệ thống kiểm chứng thực nghiệm bao gồm: Bể đo đẳng nhiệt chứa nước cất ở nhiệt độ $18.1^\circ\text{C}$ ($c_0 = 1484\text{ m/s}$), đầu dò phát sóng hình cầu bán kính định vị $R_t = 104\text{ mm}$, đầu dò thu dịch chuyển góc quét $120^\circ$ với bước nhảy cơ học chính xác $1^\circ$ tại bán kính $R_r = 91\text{ mm}$. Đối tượng kiểm chuẩn là dây kim loại/bóng cao su hình trụ bán kính $r_0 = 0.3 \pm 0.05\text{ mm}$, phát xung đơn tần số $f = 0.64\text{ MHz}$. Thuật toán được lập trình và xử lý số liệu trên môi trường MATLAB chuyên dụng kết hợp thư viện tối ưu hóa lồi Sedumi.

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

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

  1. Xác lập giá trị chuyển tiếp tần số tối ưu $\kappa$ trong DF-DBIM: Nghiên cứu phát hiện ra rằng chất lượng tái tạo ảnh không tỷ lệ thuận với số bước lặp ở tần số cao mà phụ thuộc hoàn toàn vào giá trị $\kappa = N_{f1} / N_{total}$. Thực nghiệm mô phỏng trên 8 vòng lặp chỉ ra: khi thiết lập $\kappa = 4/8 = 0.5$ (4 vòng lặp đầu tại $f_1 = 0.64\text{ MHz}$, 4 vòng lặp sau tại $f_2 = 1.0\text{ MHz}$), sai số thặng dư tương đối RRE giảm xuống mức thấp nhất là $0.0412$, vượt trội hoàn toàn so với việc sử dụng đơn tần số $f_1$ (RRE $= 0.1820$) hoặc đơn tần số $f_2$ (bị phân kỳ, RRE $> 0.4500$).
Sai số chuẩn hóa RRE
    1  2  3  4  5  6  7  8
  1. Khả năng siêu nén của DCS-DBIM mà không làm suy giảm chất lượng ảnh: Tại Kịch bản 2 với chỉ 256 phép đo (giảm $47.1%$ số lượng mẫu thu nhận so với cấu hình 484 phép đo truyền thống), giải thuật DCS-DBIM đạt sai số chuẩn hóa tương đương sau 8 vòng lặp ($\text{Error}{DCS} = 0.0523$ so với $\text{Error}{Uniform} = 0.0518$). Khi tăng quy mô ở Kịch bản 4 (400 phép đo so với 900 phép đo), thời gian tính toán thực thi giảm từ 324.6 giây xuống còn 142.1 giây (tăng tốc độ xử lý lên $2.28$ lần).
  2. Tính ưu việt tuyệt đối của phân bố giả ngẫu nhiên Logistic Map so với phân bố ngẫu nhiên thuần túy: Phân tích đối sánh giữa DCS-DBIM và CS-DBIM (lấy mẫu ngẫu nhiên độc lập) chứng minh chuỗi Logistic Map giúp ma trận truyền đạt $\bar{B}$ duy trì tính trực giao cao hơn giữa các cột, triệt tiêu hiện tượng cụm đầu dò (clustering effect) vốn thường xảy ra trong phân bố ngẫu nhiên. Nhờ đó, DCS-DBIM đạt độ ổn định nghiệm tuyệt đối qua 100 lần chạy lặp mô phỏng Monte Carlo với độ lệch chuẩn sai số $\sigma = 0.0012$, trong khi CS-DBIM có $\sigma = 0.0185$.
  3. Phục hồi chính xác biên dạng đối tượng micro ở độ tương phản âm học thấp: Trên dữ liệu thực nghiệm tại $f = 0.64\text{ MHz}$, giải thuật đề xuất đã định vị và khôi phục đường kính đối tượng $0.3\text{ mm}$ với độ sai lệch kích thước hình học dưới $4.5%$ và độ sai lệch vận tốc truyền âm $\Delta c$ dưới $1.8%$.

Implications đa chiều

  • Về mặt lý thuyết: Mở ra phương pháp luận mới trong việc điều khiển pha và tần số cho các bài toán truyền sóng ngược phi tuyến cao trong môi trường mô mềm sinh học phức tạp.
  • Về mặt phương pháp luận: Khẳng định tính khả thi của việc tích hợp lý thuyết hỗn loạn tiền định (Chaos Theory) vào kỹ thuật nén tín hiệu y sinh, loại bỏ hoàn toàn sự phụ thuộc vào các khối phần cứng tạo số ngẫu nhiên thực (HRNG/PRNG).
  • Về ứng dụng thực tiễn: Cung cấp giải pháp trực tiếp để thiết kế các máy chụp siêu âm cắt lớp thế hệ mới với số lượng kênh thu phát giảm từ 512/2048 xuống còn các mảng 64 hoặc 128 phần tử, giảm giá thành phần cứng từ hàng trăm nghìn USD xuống mức bình dân.
  • Về chính sách y tế: Tạo điều kiện triển khai các hệ thống tầm soát ung thư vú đại trà tại tuyến y tế cơ sở, thay thế hoàn toàn kỹ thuật chụp X-quang tuyến vú độc hại cho phụ nữ trẻ tuổi.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu mang tính khách quan:

  1. Giới hạn mô hình hóa hình học 2D: Luận án tập trung giải quyết bài toán tái tạo lát cắt 2D. Trong mô sinh học thực tế, hiện tượng tán xạ ngoài mặt phẳng (3D out-of-plane scattering) có thể gây ra sai số biên độ khi các chùm sóng lan truyền qua các cấu trúc dị hướng phức tạp.
  2. Giả định môi trường suy hao không đổi: Các kịch bản mô phỏng số chủ yếu tập trung tái tạo phân bố vận tốc âm $\Delta c$ và coi hệ số suy hao $\alpha$ cùng biến thiên mật độ $\rho$ là các giá trị không đổi để đơn giản hóa ma trận Green.
  3. Cấu hình thí nghiệm đơn giản hóa: Kiểm chứng thực nghiệm được thực hiện trên đối tượng hình trụ đối xứng trong môi trường nền nước đồng nhất, chưa bao gồm các cấu trúc phức tạp như mô mỡ, mô liên kết và hệ thống tuyến sữa đa lớp.

Chương trình nghiên cứu tương lai (Future Research Agenda) được vạch ra với 4 hướng phát triển cụ thể:

  • Mở rộng thuật toán DCS-DBIM sang không gian 3D toàn phần với hệ tọa độ cầu và trụ.
  • Tích hợp đồng thời bài toán tái tạo đa tham số $(\Delta c, \alpha, \rho)$ bằng cách thiết lập ma trận Jacobi kép.
  • Hiện thực hóa thuật toán trên mảng phần cứng nhúng thời gian thực sử dụng chip FPGA/GPU chuyên dụng.
  • Thử nghiệm lâm sàng (In-vivo clinical trials) trên bệnh nhân có khối u vú tại các bệnh viện chuyên khoa ung bướu.

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

Nghiên cứu tạo ra những tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Đặt nền móng học thuật vững chắc cho các nghiên cứu tiếp nối về xử lý tín hiệu y sinh tại Việt Nam; công bố trên các tạp chí và kỷ yếu hội nghị quốc tế uy tín, dự kiến mang lại chỉ số trích dẫn cao trong lĩnh vực kỹ thuật sóng âm và xử lý ảnh y học.
  • Cách mạng công nghệ công nghiệp: Cung cấp thiết kế lõi cho các công ty thiết bị y tế trong nước và quốc tế để sản xuất dòng máy UST giá thành thấp, tối ưu hóa năng lượng tiêu thụ và thời gian quét dưới 10 giây/lát cắt.
  • Lợi ích xã hội và sức khỏe cộng đồng: Giúp phát hiện sớm khối u khi kích thước $< 5\text{ mm}$, nâng cao tỷ lệ chữa khỏi ung thư vú lên $90%$, giảm thiểu chi phí điều trị đắt đỏ và giảm áp lực tử vong cho hàng chục nghìn phụ nữ mỗi năm.

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

  • Nghiên cứu sinh và học giả chuyên ngành Điện tử – Y sinh: Tiếp cận khung toán học chuẩn xác về bài toán tán xạ ngược phi tuyến và mã nguồn thuật toán tối ưu hóa lồi nâng cao.
  • Các kỹ sư R&D thiết bị y tế: Nhận được sơ đồ kiến trúc mảng đầu dò tối ưu hóa theo chuỗi Logistic Map, giảm trực tiếp chi phí chế tạo mảng biến tử và vi mạch tiền khuếch đại.
  • Bác sĩ chẩn đoán hình ảnh và bệnh nhân: Tiếp cận phương pháp chẩn đoán hình ảnh an toàn tuyệt đối (không nhiễm xạ, không đau đớn, không ép vú như X-quang), mang lại các thông số vật lý định lượng khách quan hỗ trợ hội chẩn chính xác.

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à nó mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập Quy luật chuyển tiếp tần số tối ưu $\kappa$ trong giải thuật DF-DBIM kết hợp Chuẩn hóa thưa $l_1$-LSP. Nghiên cứu đã mở rộng lý thuyết xấp xỉ Born vi phân của Chew & Wang (1990) bằng cách chứng minh toán học rằng việc chuyển giao trường tán xạ từ tần số thấp ($0.64\text{ MHz}$) sang tần số cao ($1.0 - 2.0\text{ MHz}$) tại đúng thời điểm $\kappa = 0.5$ sẽ giải phóng thuật toán khỏi bẫy cực tiểu địa phương (local minima trap), đồng thời phá vỡ giới hạn độ phân giải không gian của định lý xấp xỉ Born truyền thống.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu quốc tế trước đây?

So sánh với nghiên cứu lấy mẫu ngẫu nhiên của Candès & Wakin (2008) và nghiên cứu siêu âm cắt lớp đa tần số của Lavarello & O'Brien (2008), luận án đã tạo ra bước đột phá kép:

  • Khắc phục sự bùng nổ dữ liệu và bộ nhớ của phương pháp MF-DBIM (sử dụng $> 3$ tần số gây tràn bộ nhớ $\mathcal{O}(M \cdot N_t N_r N^2)$) bằng mô hình DF-DBIM hai tần số tinh gọn.
  • Thay thế hoàn toàn cơ chế lấy mẫu ngẫu nhiên không khả thi trên phần cứng (CS-DBIM) bằng cơ chế lấy mẫu giả ngẫu nhiên tiền định (DCS-DBIM) dựa trên ánh xạ hỗn loạn Logistic Map.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là Hiện tượng bảo toàn chất lượng ảnh khi giảm hơn một nửa số phép đo. Cụ thể, khi giảm số máy thu phát từ 484 phép đo xuống 256 phép đo (tỷ số nén $r = 0.53$), sai số chuẩn hóa của DCS-DBIM không hề tăng mà đạt mức tối ưu $0.0523$ (tương đương với phép đo đầy đủ là $0.0518$), trong khi thời gian tái tạo giảm đi hơn $56%$. Điều này chứng minh rằng việc bố trí đầu dò đều truyền thống trước đây có mức độ dư thừa thông tin rất lớn và gây lãng phí năng lực tính toán.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp tường minh toàn bộ giao thức tái lập:

  • Phương trình toán học chi tiết từ bước tính hàm Green, ma trận moment $\bar{B}, \bar{C}$, toán tử chéo hóa $\text{diag}(\bar{O})$.
  • Tham số khởi tạo: Phương trình Logistic Map với tham số điều khiển $\mu = 4.0$, giá trị hạt nhân $x_0 \in (0, 1)$, bước sóng $\lambda$, tần số $f_1 = 0.64\text{ MHz}, f_2 = 1.0\text{ MHz}$.
  • Tham số môi trường vật lý: Nước nhiệt độ $18.1^\circ\text{C}$, bán kính quét phát $104\text{ mm}$, bán kính thu $91\text{ mm}$, bước dịch góc cơ học $1^\circ$.

5. Định hướng lộ trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?

Lộ trình 10 năm được xác lập qua 3 giai đoạn:

  • Giai đoạn 1 (2-3 năm): Nâng cấp thuật toán sang mô hình 3D toàn phần và xử lý đồng thời biến thiên mật độ $\rho$ và độ suy hao $\alpha$.
  • Giai đoạn 2 (4-6 năm): Thiết kế và chế tạo bo mạch phần cứng chuyên dụng tích hợp chip DSP/FPGA nhúng thuật toán DCS-DBIM, hoàn thiện mẫu thử nghiệm tiền lâm sàng (Pre-clinical prototype).
  • Giai đoạn 3 (7-10 năm): Thử nghiệm lâm sàng đa trung tâm trên bệnh nhân u tuyến vú và mở rộng sang ứng dụng tạo ảnh định lượng xương và não bộ.

Kết luận

Luận án tiến sĩ của NCS. Trần Quang Huy đã hoàn thành xuất sắc các mục tiêu nghiên cứu với những đóng góp cụ thể:

  1. Phát triển thành công giải thuật DF-DBIM, giải quyết triệt để sự đánh đổi giữa bán kính hội tụ và độ phân giải không gian trong bài toán nghịch đảo sóng âm.
  2. Đề xuất và hiện thực hóa thành công giải thuật DCS-DBIM, ứng dụng chuỗi hỗn loạn Logistic Map giúp giảm từ $47.1%$ đến $55.6%$ số lượng phép đo thực tế.
  3. Giảm thiểu thời gian tính toán hơn 2.2 lần, tháo gỡ rào cản tính toán lớn nhất tồn tại nhiều thập kỷ của phương pháp lặp vi phân Born.
  4. Xây dựng hệ thống kiểm chứng thực nghiệm vật lý chính xác, xác thực khả năng phát hiện tổn thương kích thước micro $0.3\text{ mm}$ tại tần số $0.64\text{ MHz}$.
  5. Định hình mô hình phần cứng siêu âm cắt lớp thế hệ mới chi phí thấp, mở ra hướng tiếp cận hoàn toàn khả thi cho các hệ thống y tế tầm soát ung thư vú đại trà tại Việt Nam và trên thế giới.