Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Bệ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 phạm vi toàn cầu cũng như tại Việt Nam. Hàng năm thế giới ghi nhận khoảng 1,4 triệu ca mắc mới và 458.000 trường hợp tử vong. Tại Việt Nam, số ca mắc mới mỗi năm khoảng 7.000 người với tỷ lệ tử vong ước tính khoảng 35%, đồng thời độ tuổi mắc bệnh có xu hướng trẻ hóa. Việc phát hiện khối u ở giai đoạn sớm—khi kích thước tổn thương có đường kính dưới 5 mm—cho phép nâng tỷ lệ điều trị thành công lên tới 90%. Kỹ thuật chụp X-quang tuyến vú hiện hành bộc lộ hạn chế rõ rệt đối với phụ nữ dưới 50 tuổi do mật độ mô tuyến vú dày đặc, không tạo ra sự khác biệt rõ rệt về độ cản quang để nhận diện tổn thương nhỏ. Phương pháp siêu âm B-mode truyền thống chỉ cung cấp hình ảnh định tính dựa trên sự thay đổi trở kháng âm, phụ thuộc lớn vào kỹ năng của người thao tác và khó phát hiện khối u ở giai đoạn khởi phát.

Kỹ thuật siêu âm cắt lớp dựa trên lý thuyết tán xạ cho phép xác định định lượng các thông số âm học của mô, đặc biệt là sự thay đổi tốc độ truyền âm. Trong đó, phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM) được đánh giá là công cụ toán học thích hợp để giải bài toán ngược phi tuyến nhằm tái tạo cấu trúc không đồng nhất. Tuy nhiên, rào cản lớn nhất khiến các thiết bị siêu âm cắt lớp chưa thể thương mại hóa rộng rãi là độ phức tạp tính toán rất cao ($O(N_t N_r N^2)$), dẫn đến thời gian xử lý kéo dài và nhạy cảm với nhiễu đo đạc. Bên cạnh đó, việc áp dụng lý thuyết lấy mẫu nén thông thường đòi hỏi cấu hình bố trí đầu dò ngẫu nhiên, gây khó khăn lớn khi thiết kế phần cứng thực tế. Do đó, khoảng trống nghiên cứu chính mà tác giả xác định là:

  • Thiếu giải pháp tối ưu hóa thời điểm chuyển giao tần số trong kỹ thuật kết hợp đa tần nhằm vừa đảm bảo tính hội tụ vừa cải thiện độ phân giải không gian.
  • Thiếu giải thuật lấy mẫu nén có cấu trúc tiền định (giả ngẫu nhiên) tương thích với phần cứng hệ thống đo để giảm số lượng đầu dò và khối lượng tính toán.

Mục tiêu và câu hỏi nghiên cứu

Mục đích nghiên cứu của luận án là xây dựng và phát triển các phương pháp, mô hình, giải thuật và công cụ nhằm tăng tốc độ lấy mẫu, giảm thiểu thời gian tính toán và nâng cao chất lượng tạo ảnh siêu âm cắt lớp thông qua việc tích hợp kỹ thuật kết hợp tần số và kỹ thuật lấy mẫu nén giả ngẫu nhiên vào phương pháp DBIM.

Các mục tiêu và nhiệm vụ cụ thể được đánh số theo văn bản luận án gồm:

  1. Nghiên cứu đề xuất kỹ thuật kết hợp tần số hiệu quả ứng dụng trong DBIM nhằm nâng cao chất lượng tạo ảnh siêu âm cắt lớp, cụ thể là tăng tốc độ hội tụ và cải thiện độ phân giải ảnh. Trọng tâm là xác định thời điểm chuyển giao tối ưu giữa tần số thấp và tần số cao.
  2. Nghiên cứu đề xuất kỹ thuật lấy mẫu nén giả ngẫu nhiên ứng dụng trong DBIM nhằm nâng cao chất lượng tạo ảnh siêu âm cắt lớp, cụ thể là giảm đáng kể số lượng đầu dò thu phát, giảm độ phức tạp tính toán mà vẫn đảm bảo độ chính xác khôi phục của ảnh.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu:
    • Sự thay đổi tốc độ truyền âm khi sóng siêu âm tương tác với môi trường không đồng nhất (mô bệnh lý, khối u).
    • Kỹ thuật tạo ảnh dựa trên trường tán xạ ngược của sóng siêu âm.
    • Phương pháp lặp vi phân Born (DBIM) giải bài toán tán xạ ngược.
    • Các kỹ thuật xử lý tín hiệu kết hợp: kỹ thuật kết hợp hai tần số (Dual Frequency - DF) và kỹ thuật lấy mẫu nén giả ngẫu nhiên (Deterministic Compressed Sensing - DCS).
  • Phạm vi nghiên cứu:
    • Không gian: Luận án giới hạn việc khảo sát trên mô hình lát cắt hai chiều (2D) trong hệ tọa độ trụ nhằm làm cơ sở dữ liệu để ghép nối thành ảnh ba chiều (3D). Vùng quan tâm (Region of Interest - RoI) chứa đối tượng hình trụ tròn đặt tại tâm hệ đo trong môi trường nền đồng nhất (nước).
    • Vật lý: Xét các đối tượng có kích thước tương đương hoặc nhỏ hơn bước sóng của sóng tới; mức độ biến thiên tốc độ truyền âm tương đối nhỏ ($\Delta c$ khoảng vài phần trăm), cho phép bỏ qua ảnh hưởng của phản xạ mạnh và khúc xạ tại ranh giới phân cách.
    • Dữ liệu và thực nghiệm: Sử dụng dữ liệu mô phỏng số và dữ liệu thực nghiệm đo đạc áp suất sóng tán xạ của vật thể hình cầu đường kính 0,3 mm trong bể nước ở tần số 0,64 MHz, nhiệt độ phòng thí nghiệm 18,1°C.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước trong y văn

Trong chương tổng quan, tác giả hệ thống hóa các hướng tiếp cận chính trong lịch sử phát triển kỹ thuật siêu âm cắt lớp y sinh:

  1. Siêu âm cắt lớp dựa trên lý thuyết chiếu tia (Ray-based Tomography): Được phát triển từ cuối những năm 1970 và đầu thập niên 1980 bởi Greenleaf và các cộng sự nhằm tái tạo lát cắt tốc độ âm và độ suy hao tuyến vú. Các hệ thống thực nghiệm tiêu biểu gồm hệ thống CURE của Viện Ung thư Karmanos (sử dụng phép chiếu tia thẳng và hiệu chỉnh khúc xạ), hệ thống HUTT của Đại học Nam California (khôi phục hệ số suy hao ở các tần số khác nhau để tái tạo ảnh 3D), và hệ thống siêu âm cắt lớp 3D của Viện Công nghệ Karlsruhe sử dụng mảng đầu dò hình trụ 2D. Zhang và cộng sự cũng mở rộng lý thuyết chiếu tia để tái tạo tham số phi tuyến từ sóng hài bậc hai. Tuy nhiên, các nghiên cứu mô phỏng và thực nghiệm chỉ ra rằng lý thuyết chiếu tia chỉ áp dụng chính xác cho các cấu trúc có kích thước lớn hơn bước sóng từ 2 đến 5 lần; khi cấu trúc có kích thước tương đương bước sóng, hiện tượng tán xạ chiếm ưu thế khiến mô hình chiếu tia bị sai lệch nghiêm trọng.
  2. Siêu âm cắt lớp dựa trên lý thuyết tán xạ (Scattering-based Tomography):
    • Xấp xỉ Rytov bậc 1: Không hạn chế kích thước đối tượng khảo sát nhưng gặp phải trở ngại lớn về hiện tượng gập pha (phase wrapping) khi ước tính pha tổng từ dữ liệu đo, dẫn đến sai số đột ngột trong phép nghịch đảo phương trình sóng.
    • Xấp xỉ Born bậc 1: Slaney và các cộng sự đã chứng minh phương pháp này hoạt động tin cậy khi độ dịch pha cực đại thỏa mãn điều kiện $|\Delta\phi| < \pi$. Khi vượt quá ngưỡng này, ảnh tái tạo xuất hiện nhiều thành phần giả chia cắt ở trung tâm đối tượng.
  3. Các cải tiến đối với phương pháp lặp vi phân Born (DBIM):
    • Kết hợp BIM và DBIM: Nhóm tác giả trong công trình tham khảo [94] đề xuất giải pháp kết hợp: sử dụng phương pháp lặp Born (BIM) trong những vòng lặp đầu tiên để chống nhiễu, sau đó chuyển sang DBIM để tăng tốc độ hội tụ.
    • Thông tin mật độ mô: Các nghiên cứu của Lavarello và cộng sự chỉ ra sự biến thiên mật độ có độ lớn tương đương biến thiên tốc độ truyền âm, việc khôi phục đồng thời mật độ cung cấp thêm thuộc tính cơ học của mô mềm.
    • Kỹ thuật nội suy: Ứng dụng phép nội suy láng giềng gần nhất (Nearest Neighbor) để tái tạo ảnh ở độ phân giải thấp trong các vòng lặp đầu, sau đó nâng độ phân giải ở các vòng lặp sau nhằm tiết kiệm chi phí tính toán.
    • Kỹ thuật đa tần (MF-DBIM và DF-DBIM): Sử dụng dải tần số từ thấp đến cao để khắc phục giới hạn xấp xỉ Born trong môi trường tán xạ mạnh. Tần số thấp giúp đảm bảo sự hội tụ toàn cục, tần số cao giúp cải thiện độ phân giải không gian.
    • Kỹ thuật lấy mẫu nén (Compressed Sensing - CS): Khởi xướng bởi Candes, Tao và Donoho (năm 2006). Các nghiên cứu áp dụng CS vào siêu âm cắt lớp bố trí máy thu phân bố ngẫu nhiên để khôi phục tín hiệu thưa với số lượng phép đo nhỏ hơn giới hạn Nyquist. Tại Việt Nam, nhóm nghiên cứu tại Trường Đại học Công nghệ - ĐHQGHN là đơn vị thực hiện các công bố quốc tế về hướng tiếp cận này.

Khoảng trống luận án lựa chọn giải quyết

Mặc dù các nghiên cứu trước đã ghi nhận tiềm năng của kỹ thuật kết hợp tần số và lấy mẫu nén, nhưng vẫn tồn tại hai hạn chế then chốt:

  • Trong kỹ thuật kết hợp tần số: Việc xác định số vòng lặp tối ưu ở tần số thấp trước khi chuyển sang tần số cao ($N_{f1}$) chưa có quy tắc định lượng cụ thể. Nếu chuyển giao sai thời điểm, chất lượng ảnh khôi phục có thể kém hơn so với việc chỉ sử dụng một tần số đơn lẻ. Luận án lựa chọn giải quyết bài toán tìm giá trị $N_{f1}$ tối ưu gắn liền với mối tương quan giữa số phép đo và số điểm ảnh.
  • Trong kỹ thuật lấy mẫu nén: Việc triển khai mảng đầu dò ngẫu nhiên đòi hỏi phần cứng tạo số ngẫu nhiên thực (HRNG) hoặc giả ngẫu nhiên (PRNG), gây phức tạp và khó khả thi khi chế tạo hệ thống đo vật lý. Luận án lựa chọn giải pháp lấy mẫu nén tiền định/giả ngẫu nhiên (DCS) dựa trên ánh xạ hỗn loạn logistic để cố định vị trí các đầu dò thu, giúp hệ thống vừa giảm được số phép đo vừa dễ dàng hiện thực hóa trên phần cứng.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung phân tích lý thuyết và toán học

Sự lan truyền của sóng siêu âm đơn sắc trong môi trường không đồng nhất hai chiều được mô tả thông qua phương trình sóng Helmholtz: $$\nabla^2 p(\vec{r}) + k^2(\vec{r}) p(\vec{r}) = 0$$

Trong đó áp suất âm tổng $p(\vec{r})$ là tổng của áp suất sóng tới $p_{inc}(\vec{r})$ và áp suất tán xạ $p_{sc}(\vec{r})$: $$p(\vec{r}) = p_{inc}(\vec{r}) + p_{sc}(\vec{r})$$

Số sóng trong môi trường nền đồng nhất là $k_0 = \omega / c_0$, và số sóng trong môi trường không đồng nhất là $k(\vec{r}) = \omega / c(\vec{r})$. Hàm mục tiêu âm học $O(\vec{r})$ đại diện cho sự sai khác về tốc độ truyền âm: $$O(\vec{r}) = k^2(\vec{r}) - k_0^2 = k_0^2 \left[ \left(\frac{c_0}{c(\vec{r})}\right)^2 - 1 \right]$$

Biểu diễn phương trình sóng dưới dạng tích phân thông qua hàm Green hai chiều $G_0(k_0 |\vec{r} - \vec{r}'|) = \frac{i}{4} H_0^{(1)}(k_0 |\vec{r} - \vec{r}'|)$, với $H_0^{(1)}$ là hàm Hankel loại một bậc 0: $$p_{sc}(\vec{r}) = \iint_{RoI} G_0(k_0 |\vec{r} - \vec{r}'|) O(\vec{r}') p(\vec{r}') d\vec{r}'$$

Để rời rạc hóa phương trình trên lưới $N \times N$ điểm ảnh, tác giả áp dụng phương pháp Moment (Method of Moments - MoM) với hệ hàm cơ bản dạng xung chữ nhật và phép lấy mẫu điểm Dirac. Hệ phương trình ma trận tương tác nội tại và tán xạ ra máy thu được xác lập: $$\bar{p} = (\bar{I} - \bar{B} \cdot \text{diag}(\bar{O}))^{-1} \bar{p}{inc}$$ $$\bar{p}{sc} = \bar{C} \cdot \text{diag}(\bar{O}) \bar{p}$$

Trong đó ma trận $\bar{B}$ biểu diễn tương tác giữa các điểm ảnh bên trong vùng chia lưới, ma trận $\bar{C}$ biểu thị tương tác từ các điểm ảnh tới các máy thu ngoài biên.

Áp dụng phép xấp xỉ Born bậc một tại mỗi vòng lặp vi phân ($n$), giả thiết độ biến thiên áp suất tán xạ phụ thuộc tuyến tính vào độ biến thiên hàm mục tiêu $\Delta \bar{O}$: $$\Delta \bar{p}_{sc} = \bar{M}_n \Delta \bar{O}_n$$

Do ma trận hệ thống $\bar{M}$ có điều kiện số xấu (ill-posed), bài toán nghịch đảo được giải thông qua phương pháp chuẩn tắc Tikhonov: $$\Delta \bar{O}n = \arg\min{\Delta \bar{O}} \left( |\bar{M}n \Delta \bar{O} - \Delta \bar{p}{sc}|^2 + \gamma |\Delta \bar{O}|^2 \right)$$

Hàm mục tiêu sau đó được cập nhật theo quy tắc lặp: $\bar{O}_{n+1} = \bar{O}_n + \Delta \bar{O}_n$.

Đối với phương pháp lấy mẫu nén giả ngẫu nhiên (DCS-DBIM), phương trình hiệu chỉnh được xây dựng dưới dạng bài toán bình phương tối thiểu chuẩn hóa $l_1$ ($l_1$-regularized least squares problem, $l_1$-LSP): $$\Delta \bar{O}n = \arg\min{\Delta \bar{O}} \left( |\bar{\Phi} \bar{M}n \Delta \bar{O} - \bar{\Phi} \Delta \bar{p}{sc}|_2^2 + \lambda |\Delta \bar{O}|_1 \right)$$

Ma trận lấy mẫu $\bar{\Phi}$ được thiết lập tiền định dựa trên chuỗi ánh xạ logistic hỗn loạn: $$x_{k+1} = \mu x_k (1 - x_k)$$ với $\mu$ là tham số điều khiển được chọn trong vùng hỗn loạn ($\mu = 4$).

Phương pháp nghiên cứu thực nghiệm và mô phỏng

Luận án kết hợp giữa phương pháp nghiên cứu lý thuyết giải thuật, mô phỏng số trên máy tính và kiểm chứng bằng dữ liệu đo đạc thực nghiệm.

  • Mô hình phát sóng: Nguồn phát được mô phỏng là chùm sóng Bessel bậc 0 lan truyền trong hệ tọa độ trụ: $$p_{inc}(\vec{r}) = J_0(k_0 |\vec{r} - \vec{r}_{Tx}|)$$
  • Hệ thống đo mô phỏng: Cấu hình tròn bao quanh đối tượng. Số lượng máy phát ($N_t$) và máy thu ($N_r$) thay đổi theo từng kịch bản thử nghiệm ($N_t, N_r \in {16, 20}$). Số điểm ảnh tái tạo trên lưới vuông $N \times N$ với $N \in {16, 21, 22, 32, 64}$.
  • Đánh giá định lượng: Độ chính xác được đánh giá qua hàm sai số chuẩn hóa (Normalized Error) và sai số thặng dư tương đối (Relative Residual Error - RRE) qua các vòng lặp: $$\text{Error} = \frac{|\bar{O}{true} - \bar{O}{recon}|2}{|\bar{O}{true}|_2}$$

Nguồn tư liệu và dữ liệu

Dữ liệu đầu vào của nghiên cứu bao gồm:

  1. Thông số vật lý tốc độ sóng âm trong các cơ quan, tổ chức sinh học của cơ thể người làm căn cứ thiết lập độ tương phản âm học:
Cơ quan / Mô sinh học Tốc độ truyền sóng âm (m/s)
Phổi 630
Mô mềm 1540
Não bộ 1550
Gan 1565
Máu 1570
Thận 1575
Cơ bắp 1600
Thần kinh thị giác 1615
Thấu kính mắt 1630
Xương 4090
  1. Thông số kỹ thuật của hệ thống thương mại SoftVue được dùng làm đối sánh yêu cầu thực tế:
    • Số biến tử mảng tròn: 2048 biến tử.
    • Số kênh thu / phát: 512 kênh thu, 512 kênh phát.
    • Tần số trung tâm: 3 MHz; độ phân giải tín hiệu: 14 bit.
    • Độ dày lát cắt: 2,5 mm; đường kính quét tối đa: 22 cm.
    • Thời gian quét một lát cắt: 15 giây.
  2. Dữ liệu thực nghiệm đo đạc áp suất sóng tán xạ: Thực hiện với đối tượng hình cầu đường kính 0,3 mm (sai số chế tạo 0,05 mm) đặt trong bể nước nhiệt độ 18,1°C. Nguồn phát sóng có tần số $f = 0,64$ MHz. Khoảng cách từ tâm hệ đo tới máy phát là 104 mm, tới máy thu là 91 mm. Máy thu quét với bước dịch chuyển $1^\circ$ trong góc mở không gian $120^\circ$.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 thiết lập bối cảnh khoa học và y học của đề tài. Tác giả phân tích gánh nặng bệnh tật của ung thư vú và chỉ rõ hạn chế của các phương pháp chẩn đoán hình ảnh truyền thống:

  • Chụp X-quang tuyến vú không hiệu quả ở phụ nữ trẻ do mô vú dày đặc gây cản trở phân biệt tổn thương.
  • Siêu âm B-mode phản xạ chỉ cho kết quả định tính và phụ thuộc chủ quan vào người khám.
  • Phương pháp siêu âm cắt lớp chiếu tia (truyền qua) chỉ tái tạo được cấu trúc kích thước lớn ($> 2-5\lambda$), hoàn toàn bỏ qua hiệu ứng nhiễu xạ và tán xạ của các tổn thương nhỏ.

Từ đó, tác giả khẳng định việc nghiên cứu siêu âm cắt lớp dựa trên lý thuyết tán xạ ngược sử dụng thuật toán DBIM là hướng đi tất yếu để phát hiện sớm các khối u kích thước dưới 5 mm. Tác giả khảo sát các công trình quốc tế liên quan đến xấp xỉ Born, xấp xỉ Rytov, các phương pháp kết hợp BIM-DBIM, kỹ thuật đa tần số và lấy mẫu nén. Trên cơ sở nhận diện các điểm nghẽn về tốc độ tính toán, sự bất ổn định khi chuyển giao tần số và tính bất khả thi khi triển khai phần cứng lấy mẫu ngẫu nhiên, tác giả xác định hai trọng tâm nghiên cứu: kỹ thuật kết hợp hai tần số (DF-DBIM) và kỹ thuật lấy mẫu nén giả ngẫu nhiên (DCS-DBIM).

Chương 2: Kỹ thuật tạo ảnh siêu âm cắt lớp

Chương 2 trình bày chi tiết cơ sở toán - lý và quy trình thực thi giải thuật DBIM chuẩn:

  • Xây dựng mô hình tán xạ 2D dựa trên phương trình tích phân Lippmann-Schwinger, hàm Green Hankel và phương pháp MoM.
  • Thiết lập giải thuật DBIM chuẩn (Giải thuật 2.1) với việc cập nhật ma trận hàm Green và giải bài toán nghịch đảo Tikhonov ở từng bước lặp. Độ phức tạp tính toán của mỗi vòng lặp được chỉ ra là $O(N_t N_r N^2)$.
  • Tiến hành thực nghiệm kiểm chứng tính đúng đắn của mô hình: So sánh tín hiệu tán xạ tính toán từ mô phỏng số với dữ liệu thực nghiệm đo từ vật thể hình cầu đường kính 0,3 mm đặt trong nước ở tần số 0,64 MHz ($R_{Tx} = 104$ mm, $R_{Rx} = 91$ mm, bước quét $1^\circ$ trên cung $120^\circ$). Kết quả sau khi lọc trung bình cho thấy biên độ và pha của dạng sóng mô phỏng tương thích chặt chẽ với số liệu thực nghiệm, xác lập độ tin cậy của mô hình toán để phát triển các thuật toán ở chương tiếp theo.

Chương 3: Thuật toán kết hợp tần số dùng trong tạo ảnh siêu âm cắt lớp

Chương 3 tập trung giải quyết bài toán nâng cao độ phân giải và tốc độ hội tụ bằng giải thuật DF-DBIM (Dual Frequency DBIM):

  • Tác giả lập luận việc sử dụng hệ thống đa tần số (MF-DBIM) đòi hỏi dung lượng bộ nhớ quá lớn và thời gian đo đạc kéo dài đối với bệnh nhân. Do đó, mô hình hai tần số DF-DBIM là sự dung hòa tối ưu giữa khối lượng tính toán và chất lượng tạo ảnh.
  • Quy trình giải thuật DF-DBIM: Sử dụng tần số thấp $f_1$ trong $N_{f1}$ vòng lặp đầu tiên để tái tạo cấu trúc tổng thể và đảm bảo tính hội tụ toàn cục; sau đó dùng tần số cao $f_2$ trong các vòng lặp còn lại ($N_{iter} - N_{f1}$, với tổng số vòng lặp $N_{iter} = 8$) để tái tạo chi tiết biên và nâng cao độ phân giải không gian.
  • Khảo sát sự phụ thuộc của chất lượng ảnh vào giá trị $N_{f1}$ qua 4 kịch bản mô phỏng trên đối tượng lý tưởng kích thước lưới $N = 22$:
    • Khảo sát các giá trị $N_{f1}$ từ 1 đến 7.
    • Phân tích mối quan hệ giữa số lượng phép đo ($N_t \times N_r$) và số biến cần tìm ($N^2$).
    • Kết quả chỉ ra rằng việc lựa chọn chính xác thời điểm chuyển giao $N_{f1}$ cho phép sai số chuẩn hóa giảm nhanh hơn và đạt giá trị thấp hơn rõ rệt so với thuật toán DBIM đơn tần số truyền thống. Nếu chuyển giao quá sớm khi tần số thấp chưa hội tụ đủ, ảnh sẽ bị phân kỳ hoặc xuất hiện nhiễu giả; nếu chuyển giao quá muộn, lợi ích cải thiện độ phân giải của tần số cao bị triệt tiêu.

Chương 4: Thuật toán lấy mẫu nén dùng trong tạo ảnh siêu âm cắt lớp

Chương 4 đề xuất giải thuật DCS-DBIM nhằm giảm thiểu số lượng phép đo và tăng tốc độ tái tạo ảnh:

  • Phân tích nhược điểm của kỹ thuật lấy mẫu nén ngẫu nhiên (CS-DBIM): Việc bố trí đầu dò ngẫu nhiên đòi hỏi mạch tạo số ngẫu nhiên phức tạp, gây trở ngại khi chế tạo đầu dò thực tế.
  • Đề xuất kỹ thuật lấy mẫu nén giả ngẫu nhiên (DCS-DBIM) sử dụng ánh xạ logistic để xác định tọa độ góc của các máy thu xung quanh đối tượng. Ma trận đo trở thành ma trận tiền định, hoàn toàn khả thi khi lập trình điều khiển phần cứng.
  • Chuyển đổi bài toán nghịch đảo từ chuẩn tắc Tikhonov sang bài toán bình phương tối thiểu có chuẩn tắc $l_1$ ($l_1$-LSP), khai thác tính thưa của ảnh vi phân mục tiêu $\Delta \bar{O}$.
  • Thực nghiệm mô phỏng trên lưới điểm ảnh $N = 21$:
    • Cấu hình đo truyền thống: $N_t = 20, N_r = 20$ (tổng cộng 400 hoặc 484 phép đo đối với các kịch bản đối chứng).
    • Cấu hình đề xuất DCS-DBIM: $N_t = 16, N_r = 16$ với máy thu phân bố giả ngẫu nhiên (tổng cộng 256 phép đo).
    • Kết quả so sánh cho thấy phương pháp DCS-DBIM với 256 phép đo đạt sai số chuẩn hóa và độ sắc nét tương đương phương pháp DBIM truyền thống cần tới 484 phép đo. Số lượng phép đo giảm gần 47%, dẫn đến thời gian tính toán và dung lượng bộ nhớ xử lý giảm đáng kể sau 8 vòng lặp.
    • So sánh giữa DCS-DBIM và CS-DBIM ngẫu nhiên chứng minh hiệu suất khôi phục của cấu hình giả ngẫu nhiên tương đương với cấu hình ngẫu nhiên lý thuyết, khẳng định tính đúng đắn của giải pháp thiết kế mảng thu tiền định.

Kết quả và những đóng góp mới

Đóng góp mới về mặt khoa học và lý luận

  1. Hoàn thiện phương pháp kết hợp hai tần số (DF-DBIM) trong siêu âm cắt lớp tán xạ: Luận án đã chỉ ra quy luật định lượng về sự phụ thuộc của chất lượng ảnh tái tạo vào thời điểm chuyển giao tần số $N_{f1}$. Tác giả chứng minh rằng việc kết hợp tần số thấp và tần số cao theo một lịch trình lặp tối ưu giúp khắc phục hiện tượng mất hội tụ của xấp xỉ Born trong môi trường tán xạ, đồng thời đạt độ phân giải không gian cao hơn so với việc chỉ sử dụng một tần số đơn lẻ.
  2. Xây dựng khung lý thuyết cho phương pháp lặp vi phân Born kết hợp lấy mẫu nén giả ngẫu nhiên (DCS-DBIM): Luận án là công trình tiên phong tại Việt Nam tích hợp chuỗi ánh xạ logistic hỗn loạn vào bài toán thiết kế cấu hình đo siêu âm cắt lớp, chứng minh khả năng áp dụng của bài toán $l_1$-LSP trong việc khôi phục ảnh tán xạ ngược với số lượng phép đo dưới ngưỡng lấy mẫu thông thường.

Đóng góp mới về mặt thực tiễn và ứng dụng

  1. Giảm thiểu chi phí phần cứng và số lượng phép đo: Giải pháp DCS-DBIM cho phép giảm số lượng đầu dò từ $N_t = N_r = 20$ xuống $N_t = N_r = 16$ (giảm từ 484 phép đo xuống 256 phép đo trong kịch bản khảo sát), giúp tinh gọn cấu trúc cơ khí của mảng đầu dò siêu âm.
  2. Tăng tốc độ tạo ảnh: Việc giảm số lượng phép đo trực tiếp làm giảm kích thước ma trận hệ thống $\bar{M}$, từ đó rút ngắn đáng kể thời gian tính toán của thuật toán lặp vi phân Born, tạo điều kiện kỹ thuật để nâng cao tốc độ quét trong chẩn đoán y tế.
  3. Tính khả thi trong chế tạo thiết bị: Khác với các mô hình lấy mẫu nén ngẫu nhiên thuần túy, cấu hình máy thu giả ngẫu nhiên tiền định của luận án có thể cài đặt trực tiếp trên các hệ thống vi điều khiển hoặc mảng cổng lập trình được (FPGA) của các máy siêu âm cắt lớp thế hệ mới mà không cần tích hợp mạch phát số ngẫu nhiên phức tạp.

Kiến nghị và giải pháp đề xuất

  • Ứng dụng mô hình chuyển giao tần số tối ưu và cấu hình đo DCS vào các quy trình xử lý tín hiệu của hệ thống siêu âm quét vú chuyên dụng nhằm nâng cao khả năng phát hiện tổn thương giai đoạn sớm.
  • Tiếp tục tối ưu hóa các tham số điều khiển của ánh xạ logistic ($\mu$) và tham số chuẩn hóa $\lambda$ trong bài toán $l_1$-LSP để thích ứng với các mức tỷ số tín hiệu trên nhiễu (SNR) khác nhau trong môi trường lâm sàng.

Hạn chế và hướng nghiên cứu tiếp

Các hạn chế được ghi nhận

Căn cứ vào phạm vi và nội dung được tác giả trình bày, luận án tồn tại một số giới hạn kỹ thuật:

  1. Mô hình hình học: Nghiên cứu mới dừng lại ở việc giải quyết bài toán tái tạo lát cắt 2 chiều (2D) trong hệ tọa độ trụ. Việc chuyển tiếp từ tập hợp các lát cắt 2D thành khối ảnh 3D hoàn chỉnh chưa được giải quyết trong khuôn khổ các giải thuật đề xuất.
  2. Đơn giản hóa đối tượng: Đối tượng khảo sát trong mô phỏng được lý tưởng hóa dưới dạng hình trụ tròn đồng chất đặt trong môi trường nền nước; các yếu tố phức tạp của mô tuyến vú thực tế (như các bè sợi, mô mỡ xen kẽ, cấu trúc giải phẫu không đều) chưa được đưa vào mô hình số.
  3. Tham số âm học: Luận án chủ yếu tập trung vào việc khôi phục duy nhất thông số tốc độ truyền âm ($\Delta c$) ở mức biến thiên nhỏ (khoảng vài phần trăm để thỏa mãn điều kiện xấp xỉ Born). Các thông số định lượng quan trọng khác như độ suy hao âm học và sự biến thiên mật độ mô ($\Delta \rho$) chưa được tích hợp đồng thời vào các giải thuật DF-DBIM và DCS-DBIM.

Hướng nghiên cứu tiếp

  • Mở rộng mô hình giải thuật toán học từ 2D sang không gian 3D đầy đủ để phản ánh chính xác trường sóng tán xạ trong thể tích khối ngực.
  • Nghiên cứu kết hợp khôi phục đa tham số: Tái tạo đồng thời tốc độ truyền âm, hệ số suy hao và mật độ khối u nhằm tăng cường độ nhạy và độ đặc hiệu trong chẩn đoán phân biệt u lành và u ác tính.
  • Nâng cấp mô hình mô phỏng với các phantom giải phẫu vú phức tạp và thử nghiệm giải thuật trên hệ thống đo đạc thực nghiệm lâm sàng đa kênh.

Giá trị tham khảo

Đối tượng thụ hưởng

  • Nghiên cứu sinh và học viên cao học: Luận án là tài liệu tham khảo có giá trị chuyên môn cao cho các ngành Kỹ thuật điện tử, Kỹ thuật viễn thông, Kỹ thuật y sinh và Xử lý tín hiệu số, đặc biệt trong các chủ đề về giải bài toán ngược phi tuyến, xấp xỉ Born và ứng dụng lý thuyết lấy mẫu nén.
  • Giảng viên và nhà nghiên cứu: Cung cấp tài liệu tham khảo bài bản về việc kết hợp phương pháp số (Moment) với các kỹ thuật tối ưu hóa hiện đại ($l_1$-LSP, chuẩn tắc Tikhonov) trong bài toán lan truyền sóng âm.
  • Kỹ sư thiết kế thiết bị y tế: Cung cấp cơ sở thuật toán và cấu hình phần cứng khả thi (bố trí máy thu DCS) cho việc nghiên cứu chế tạo các hệ thống chụp ảnh siêu âm cắt lớp thế hệ mới hỗ trợ tầm soát ung thư vú.

Các nội dung đáng tham khảo

  • Quy trình toán học chi tiết rời rạc hóa phương trình sóng Helmholtz bằng phương pháp Moment và thiết lập ma trận Green tương tác trong Chương 2.
  • Phương pháp luận khảo sát và lựa chọn thời điểm chuyển giao tần số $N_{f1}$ trong thuật toán DF-DBIM tại Chương 3.
  • Kỹ thuật xây dựng ma trận đo tiền định sử dụng ánh xạ logistic và quy trình giải bài toán $l_1$-LSP trong thuật toán DCS-DBIM tại Chương 4.

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

1. Tại sao kỹ thuật siêu âm cắt lớp tán xạ lại có khả năng phát hiện khối u nhỏ tốt hơn siêu âm B-mode thông thường?

Siêu âm B-mode hoạt động dựa trên nguyên lý ghi nhận sóng phản xạ từ các mặt phân cách âm học lớn, chỉ hiển thị định tính hàm trở kháng âm và gặp khó khăn khi phát hiện các tổn thương có kích thước nhỏ hơn bước sóng hoặc các tổn thương nằm sâu trong mô dày. Ngược lại, kỹ thuật siêu âm cắt lớp tán xạ ghi nhận toàn bộ trường sóng tán xạ theo mọi hướng ($0^\circ - 360^\circ$) khi sóng âm tương tác với vật thể có kích thước tương đương hoặc nhỏ hơn bước sóng. Dựa trên việc giải bài toán ngược phương trình sóng vi phân, phương pháp này khôi phục được định lượng giá trị tốc độ truyền âm tại từng điểm ảnh, giúp nhận diện chính xác vị trí và hình dạng của các khối u ngay từ giai đoạn khởi phát (dưới 5 mm).

2. Giới hạn vật lý nào quyết định tính khả dụng của phương pháp xấp xỉ Born bậc 1 trong siêu âm cắt lớp?

Theo nghiên cứu của Slaney và các cộng sự được dẫn trong luận án, phương pháp xấp xỉ Born bậc 1 chỉ hoạt động ổn định và chính xác khi áp suất sóng tán xạ rất nhỏ so với áp suất sóng tới. Điều kiện này tương đương với yêu cầu độ dịch pha cực đại của sóng khi truyền qua môi trường không đồng nhất phải thỏa mãn $|\Delta\phi| < \pi$. Khi độ dịch pha vượt ngưỡng này ($|\Delta\phi| \ge \pi$), phương pháp sẽ bị mất hiệu quả, khiến hình ảnh khôi phục xuất hiện nhiều thành phần nhiễu đốm giả tạo ở khu vực trung tâm đối tượng.

3. Ưu điểm nổi bật của giải thuật lấy mẫu nén giả ngẫu nhiên (DCS-DBIM) so với lấy mẫu nén ngẫu nhiên (CS-DBIM) là gì?

Kỹ thuật lấy mẫu nén ngẫu nhiên (CS-DBIM) yêu cầu các đầu dò thu tín hiệu phải được bố trí ở các vị trí ngẫu nhiên trên hệ đo, điều này đòi hỏi hệ thống phần cứng phải tích hợp bộ phát số ngẫu nhiên thực (HRNG) hoặc giả ngẫu nhiên (PRNG), gây phức tạp lớn trong thiết kế vi mạch và chế tạo cơ khí. Trong khi đó, giải thuật DCS-DBIM sử dụng chuỗi ánh xạ logistic có tính chất tiền định để phân bố vị trí máy thu. Điều này cho phép cố định cấu hình phần cứng đo đạc mà vẫn đảm bảo tính chất tán xạ thông tin cần thiết của lý thuyết lấy mẫu nén, giúp giảm số phép đo từ 484 xuống 256 mà chất lượng ảnh vẫn tương đương, tạo tính khả thi cao khi ứng dụng vào thiết bị thực tế.

4. Tại sao tác giả lại lựa chọn mô hình hai tần số (DF-DBIM) thay vì mô hình đa tần số (MF-DBIM)?

Về mặt lý thuyết, phương pháp đa tần số (MF-DBIM) với nhiều bước nhảy tần số kế tiếp sẽ mang lại chất lượng ảnh và tốc độ hội tụ tốt nhất. Tuy nhiên, việc phát và thu tín hiệu ở quá nhiều tần số đòi hỏi phải thu thập, truyền tải và lưu trữ một khối lượng dữ liệu khổng lồ, kéo theo yêu cầu dung lượng bộ nhớ hệ thống rất lớn và kéo dài thời gian quét trên cơ thể bệnh nhân. Phương pháp hai tần số (DF-DBIM) được tác giả lựa chọn như một giải pháp dung hòa kỹ thuật: tần số thấp $f_1$ đảm bảo giải thuật hội tụ đúng hướng mà không bị phân kỳ, trong khi tần số cao $f_2$ hoàn thiện độ phân giải không gian của ảnh, qua đó vừa đảm bảo chất lượng chẩn đoán vừa giảm thiểu độ phức tạp tính toán.

5. Dữ liệu thực nghiệm kiểm chứng trong luận án được thiết lập với các thông số vật lý nào?

Để kiểm chứng mô hình toán học của phương pháp DBIM trước khi phát triển các thuật toán mới, tác giả đã sử dụng tập dữ liệu thực nghiệm đo trường áp suất tán xạ của một vật thể hình cầu đặt trong môi trường nước. Các thông số cụ thể bao gồm:

  • Đường kính quả cầu: 0,3 mm (với sai số chế tạo 0,05 mm).
  • Tần số sóng âm phát: $f = 0,64$ MHz.
  • Môi trường truyền dẫn: Nước ở nhiệt độ ổn định 18,1°C.
  • Khoảng cách hình học: Từ tâm đối tượng đến máy phát là 104 mm; từ tâm đối tượng đến máy thu là 91 mm.
  • Quy cách lấy mẫu: Máy thu quét dịch chuyển từng bước $1^\circ$ trong không gian góc mở $120^\circ$.

Kết luận

Luận án của tác giả Trần Quang Huy đã giải quyết có hệ thống bài toán nâng cao chất lượng và tốc độ tạo ảnh siêu âm cắt lớp dựa trên phương pháp lặp vi phân Born (DBIM). Bằng việc đề xuất thành công giải thuật hai tần số DF-DBIM với thời điểm chuyển giao tối ưu và giải thuật lấy mẫu nén giả ngẫu nhiên DCS-DBIM dựa trên ánh xạ logistic, nghiên cứu đã chứng minh khả năng giảm thiểu đáng kể số lượng phép đo và chi phí tính toán mà vẫn duy trì độ chính xác khôi phục cấu trúc mô. Các kết quả mô phỏng số và dữ liệu thực nghiệm kiểm chứng cung cấp cơ sở khoa học và giải pháp kỹ thuật tin cậy, đóng góp thiết thực vào định hướng phát triển các hệ thống thiết bị siêu âm cắt lớp thế hệ mới phục vụ tầm soát sớm ung thư vú.