Tổng quan về luận án

Bệnh lý tim mạch tiếp tục là nguyên nhân gây tử vong hàng đầu trên phạm vi toàn cầu. Báo cáo từ Tổ chức Y tế Thế giới (WHO) chỉ ra rằng bệnh thiếu máu cơ tim cục bộ và đột quỵ chiếm tới 27% tổng số ca tử vong trên thế giới; riêng tại Việt Nam, con số này lên tới xấp xỉ 200.000 ca mỗi năm, chiếm 1/4 tổng số trường hợp tử vong. Trong hồi sức cấp cứu, phẫu thuật lồng ngực và điều trị sốc tim, việc giám sát liên tục các thông số huyết động then chốt — đặc biệt là cung lượng tim ($CO - Cardiac\ Output$), thể tích nhát bóp ($SV - Stroke\ Volume$), và thời gian tống máu thất trái ($LVET - Left\ Ventricular\ Ejection\ Time$) — giữ vai trò quyết định đến sự sống còn của người bệnh. Mặc dù các kỹ thuật xâm lấn như đo catheter động mạch phổi ($Swan-Ganz\ PAC-TD$) và phân tích sóng mạch kết hợp pha loãng nhiệt xuyên phổi ($PiCCO\ TPTD$) cung cấp độ chính xác tiêu chuẩn vàng, chúng lại tiềm ẩn nhiều biến chứng nghiêm trọng như loạn nhịp tim, nhiễm trùng huyết, tổn thương van ba lá, và đòi hỏi chi phí y tế rất lớn.

Phương pháp tim đồ trở kháng ngực ($ICG - Impedance\ Cardiography$) dựa trên đo trở kháng điện sinh học lồng ngực ($TEB - Thoracic\ Electrical\ Bioimpedance$) nổi lên như một giải pháp không xâm lấn, liên tục và chi phí thấp. Tuy nhiên, việc ứng dụng rộng rãi ICG trong lâm sàng đối mặt với ba rào cản khoa học và kỹ thuật nghiêm trọng:

  1. Rào cản phần cứng và xử lý tín hiệu: Tín hiệu biến thiên trở kháng do hoạt động bơm máu ($\Delta Z$) có biên độ cực nhỏ (chỉ chiếm khoảng 0,5% trở kháng nền $Z_0$, dao động trong khoảng $100 - 240\ m\Omega$ trên nền $Z_0 = 20 - 48\ \Omega$). Các kiến trúc mạch tương tự ($analog$) truyền thống thường bị trôi điểm làm việc theo nhiệt độ, méo pha và kém ổn định; trong khi các hệ thống số hóa trực tiếp sóng mang cao tần đòi hỏi bộ biến đổi tương tự - số ($ADC$) tốc độ cao (lên tới $55\ MSPS$, luồng dữ liệu $880\ Mbps$), vượt quá năng lực xử lý của các vi xử lý tầm trung và làm đội giá thành thiết bị.
  2. Rào cản can thiệp lâm sàng: Cấu hình 8 điện cực điểm chuẩn của Bernstein ($1986$) đòi hỏi dán điện cực kép ở hai bên gốc cổ. Vị trí này thường xuyên bị chồng lấn với vị trí đặt catheter tĩnh mạch cảnh trong ($CVP$) ở bệnh nhân hồi sức tích cực ($ICU$), làm mất khả năng ứng dụng ICG.
  3. Rào cản nhiễu sinh lý: Nhiễu do hoạt động hô hấp ($IP - Impedance\ Pneumography$) có biên độ vượt trội và chồng lấn phổ tần số ($0,1 - 2,5\ Hz$) với tín hiệu ICG, làm trôi đường cơ sở và biến dạng các điểm đặc trưng xung yếu ($B, C, X, O$), dẫn đến sai lệch lớn khi tính toán SV và CO.

Luận án tiến sĩ kỹ thuật điện tử của tác giả Phan Đăng Hưng (Trường Đại học Bách khoa Hà Nội, 2021) với đề tài "Phát triển kỹ thuật thu nhận tín hiệu tim đồ trở kháng ngực ICG ứng dụng trong phép đo thông số cung lượng tim" dưới sự hướng dẫn của PGS.TS Vũ Duy Hải đã giải quyết trọn vẹn ba nút thắt trên.

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu ($RQ$) và kiểm chứng 3 giả thuyết khoa học ($H$):

  • $RQ1$: Làm thế nào để số hóa chính xác tín hiệu ICG mà không cần hệ thống ADC/FPGA tốc độ cao đắt đỏ? $\rightarrow H1$: Kỹ thuật số hóa đồng bộ đỉnh sóng mang cao tần cho phép khôi phục đường bao $\Delta Z$ với độ phân giải 16-bit mà chỉ cần tốc độ lấy mẫu tương đương tần số sóng mang, giảm tải băng thông xử lý hơn 275.000 lần.
  • $RQ2$: Có thể thay đổi vị trí dán điện cực vùng cổ để tránh chồng lấn đường truyền tĩnh mạch trung tâm mà vẫn đảm bảo độ tin cậy huyết động học? $\rightarrow H2$: Tồn tại các vị trí điện cực thay thế ở vùng ngực và xương đòn có mối tương quan rất cao ($R^2 > 0,9$) và sự phù hợp Bland-Altman chặt chẽ với cấu hình chuẩn Bernstein.
  • $RQ3$: Thuật toán nào loại bỏ triệt để nhiễu hô hấp mà không làm méo dạng vi phân $dZ/dt$? $\rightarrow H3$: Sự kết hợp giữa biến đổi Wavelet đa phân giải (thuật toán Mallat) và phép lấy trung bình toàn bộ ($Ensemble\ Averaging$) đồng bộ theo đỉnh R của điện tim ($ECG$) sẽ tách lọc chính xác nhiễu thở ở các trạng thái sinh lý khác nhau.

Phạm vi nghiên cứu bao gồm thiết kế phần cứng đo lường chuyên dụng cách ly y tế, xây dựng bộ công cụ phần mềm xử lý số tín hiệu, thu thập cơ sở dữ liệu thực nghiệm trên các tình nguyện viên khỏe mạnh tại Trung tâm Điện tử Y sinh - Đại học Bách khoa Hà Nội và đối chứng trực tiếp với thiết bị thương mại Medis Niccomo (Đức).

Literature Review và Positioning

Lý thuyết về sự thay đổi trở kháng vùng ngực liên quan đến thể tích máu tống đi từ thất trái bắt nguồn từ các công trình tiên phong của Atzler và Lehman, sau đó được Nyboer chính thức mô hình hóa qua phương trình biến thiên thể tích $\Delta V = \rho \cdot L_0^2 \cdot Z_0^{-2} \cdot \Delta Z$. Đến năm 1966, Kubicek cùng các cộng sự tại NASA đã tạo nên bước ngoặt khi đưa ra công thức tính thể tích nhát bóp dựa trên giá trị cực đại của đạo hàm bậc nhất trở kháng $(dZ/dt){max}$ và thời gian tống máu thất trái $LVET$: $$SV = \rho \cdot \left(\frac{L_0}{Z_0}\right)^2 \cdot \left(\frac{dZ}{dt}\right){max} \cdot LVET$$

Tuy nhiên, mô hình hình trụ đồng nhất của Kubicek bộc lộ sai số lớn khi áp dụng trên các thể trạng cơ thể khác nhau. Để khắc phục, Sramek ($1981$) và sau đó là Sramek - Bernstein ($1986$) đã hiệu chỉnh bằng việc đưa vào thể tích mô dẫn điện tham gia khử cực ($VEPT$) và hệ số hiệu chỉnh cân nhắc trọng lượng lý tưởng $\delta = IBW / ABW$: $$SV = \delta \cdot \frac{(0,17H)^3}{4,25 \cdot Z_0} \cdot \left(\frac{dZ}{dt}\right)_{max} \cdot LVET$$

Trong lịch sử phát triển kỹ thuật thu nhận tín hiệu, hai trường phái thiết kế đã tồn tại song song với nhiều tranh luận học thuật gay gắt:

  1. Trường phái mạch tương tự ($Analog\ Front-End$): Các nghiên cứu của Wang et al., Ramos et al., và các thiết bị thế hệ cũ dựa vào bộ giải điều chế biên độ tương tự, bộ lọc thông dải $80 - 120\ kHz$ và các tầng khuếch đại vi sai. Nhược điểm chí tử là độ trôi nhiệt của tụ điện, độ nhạy cảm cao với can nhiễu điện từ ($EMI$), và sự biến dạng pha phi tuyến làm lệch các điểm ngoặt $B$ và $X$ trên sóng $dZ/dt$.
  2. Trường phái số hóa trực tiếp tốc độ cao ($Digital\ Oversampling$): Sử dụng các chip $FPGA$ hiệu năng cao kết hợp $ADC$ tốc độ hàng chục $MSPS$ (như công trình của Javier et al., 2014; Medis Niccomo; Biopac NICO100C). Cách tiếp cận này giải quyết được độ ổn định số nhưng tạo ra luồng dữ liệu dư thừa cực lớn, tiêu tán công suất cao và đẩy giá thành phần cứng lên mức khó tiếp cận cho các cơ sở y tế tuyến cơ sở.

Về vấn đề can nhiễu hô hấp, Baker và Geddes ($1966$) đã chứng minh trên mô hình động vật rằng quá trình thông khí phổi làm biến thiên trở kháng lồng ngực với biên độ lớn gấp 5-10 lần biên độ tín hiệu ICG do tim. Các kỹ thuật lọc số cổ điển như bộ lọc $IIR$ hay $FIR$ thông dải thường thất bại vì phổ của nhịp thở (nhất là khi thở nhanh $0,5 - 1,5\ Hz$) trùng lặp trực tiếp với dải tần hoạt động của tim ($0,8 - 2,5\ Hz$). Phương pháp lọc thích nghi của Seppä et al. ($2001$) hay biến đổi Fourier ngắn hạn ($STFT$) đòi hỏi tín hiệu tham chiếu hô hấp phụ trợ hoặc làm suy hao biên độ đỉnh $(dZ/dt)_{max}$.

Luận án của Phan Đăng Hưng định vị chính xác vào khoảng trống nghiên cứu này: Tạo ra một cầu nối công nghệ đột phá bằng cách đề xuất kiến trúc "Số hóa đỉnh sóng mang" ($Carrier\ Peak\ Digitization$) kết hợp thuật toán biến đổi Wavelet lai ghép lấy trung bình chu kỳ dựa trên đồng bộ đỉnh R-ECG (thuật toán Pan-Tompkins). Nghiên cứu không chỉ tối ưu hóa triệt để cấu trúc phần cứng mà còn giải quyết đồng thời bài toán chồng lấn điện cực trên lâm sàng, thiết lập một tiêu chuẩn mới về hiệu quả kinh tế và độ tin cậy đo lường trong công nghệ ICG.

Đó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 và hoàn thiện khung lý thuyết về mô hình dẫn truyền dòng điện cao tần trong cơ thể sống, cụ thể:

  • Mô hình phân rã trở kháng ngực động: Xác lập tường minh phương trình phân rã tín hiệu trở kháng tổng thể $TEB(t) = Z_0(t) - \Delta Z(t) + Z_{resp}(t) + Z_{motion}(t) + n(t)$, trong đó chứng minh rằng $Z_0$ không phải là hằng số tuyệt đối mà chịu sự điều biến chậm bởi dịch mô và hô hấp.
  • Lý thuyết đồng bộ điện - cơ sinh học ($ECG - ICG\ Synchronization$): Bổ sung cơ sở lý thuyết chứng minh mối liên hệ thời gian bất biến giữa phức bộ khử cực thất $QRS$ trên ECG và điểm khởi đầu tống máu thất trái (điểm B) trên đạo hàm trở kháng $dZ/dt$. Việc định vị đỉnh R bằng thuật toán Pan-Tompkins cung cấp điểm tựa toán học vững chắc để cố định cửa sổ thời gian cho phép trung bình toàn bộ tín hiệu ICG.
  • Lý thuyết phân bố mật độ dòng điện: Chứng minh bằng thực nghiệm rằng sự dịch chuyển các điện cực kích thích và thu nhận từ vị trí chuẩn Bernstein sang các vị trí thay thế ở xương đòn và ngực trước không làm thay đổi bản chất hình thái của đường cong $dZ/dt$, mà chỉ làm thay đổi hệ số hình học của trường dẫn điện.

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

Khung phân tích của luận án là sự tích hợp đa ngành giữa ba trụ cột lý thuyết: (1) Lý thuyết dẫn truyền điện trở sinh học Maxwell - Ohm, (2) Lý thuyết phân tích đa phân giải Wavelet của Stephane Mallat, và (3) Lý thuyết huyết động học tim mạch dòng chảy tống máu thất trái.

                                 +-------------------------------------------------------------+
                                 |         NGUỒN DÒNG XOAY CHIỀU ĐIỀU CHẾ (100 kHz)           |
                                 +-------------------------------------------------------------+
                                                                |
                                                                v
+------------------------+       +-------------------------------------------------------------+       +-----------------------+
|  ĐIỆN CỰC THAY THẾ     | ----> |          LỒNG NGỰC BỆNH NHÂN (Mô hình TEB)                  | <---- | ĐIỆN CỰC CHUẨN        |
|  (Tránh Catheter CVP)  |       | Z(t) = Z0 - ΔZ(t) + Z_resp(t) + Z_motion(t) + n(t)          |       | (Bernstein 8 điện cực)|
+------------------------+       +-------------------------------------------------------------+       +-----------------------+
                                                                |
                                                                v
                                 +-------------------------------------------------------------+
                                 |       KỸ THUẬT SỐ HÓA ĐỈNH SÓNG MANG TẦN SỐ CAO            |
                                 | (Lấy mẫu tại cực đại sóng mang -> Tách trực tiếp đường bao) |
                                 +-------------------------------------------------------------+
                                                                |
                                                                v
+------------------------------------------------------------------------------------------------------------------------------+
|                                              XỬ LÝ TÍN HIỆU SỐ LAI GHÉP                                                      |
|                                                                                                                              |
|   +------------------------------------+                                     +-------------------------------------------+   |
|   |         KÊNH ĐIỆN TIM (ECG)        |                                     |        KÊNH TRỞ KHÁNG NGỰC (TEB)          |   |
|   |   Thuật toán Pan-Tompkins lọc dải  |                                     |    Tách thành phần biến thiên ΔZ và Z0    |   |
|   |    và phát hiện chính xác đỉnh R   |                                     +-------------------------------------------+   |
|   +------------------------------------+                                                           |                         |
|                     |                                                                              v                         |
|                     |                                                        +-------------------------------------------+   |
|                     |                                                        |         BIẾN ĐỔI WAVELET ĐA PHÂN GIẢI     |   |
|                     |                                                        |     (Thuật toán Mallat tách nhiễu thở)    |   |
|                     |                                                        +-------------------------------------------+   |
|                     |                                                                              |                         |
|                     +------------------------------------> <---------------------------------------+                         |
|                                                            |                                                                 |
|                                                            v                                                                 |
|                                          +-----------------------------------+                                               |
|                                          |   PHÉP TRUNG BÌNH TOÀN BỘ THEO R |                                               |
|                                          |       (Ensemble Averaging)        |                                               |
|                                          +-----------------------------------+                                               |
+------------------------------------------------------------------------------------------------------------------------------+
                                                                |
                                                                v
                                 +-------------------------------------------------------------+
                                 |            TRÍCH XUẤT THÔNG SỐ HUYẾT ĐỘNG                   |
                                 |      Điểm đặc trưng: B, C, X, O -> Tính SV, HR, LVET, CO    |
                                 +-------------------------------------------------------------+

Điều kiện biên của mô hình được xác định rõ: Dòng điện kích thích $I_{rms} = 1\ mA$, tần số $f_c = 100\ kHz$ (đáp ứng tiêu chuẩn an toàn điện sinh học IEC 60601-1), đối tượng nghiên cứu ở trạng thái nghỉ, không cử động thân mình đột ngột.

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à phương pháp tiếp cận kỹ thuật thực nghiệm định lượng cao ($Quantitative\ Engineering\ Experimentalism$). Thiết kế nghiên cứu đa tầng ($Multi-level\ Design$) bao gồm:

  1. Tầng phần cứng vật lý: Thiết kế mạch nguồn cách ly y tế cao cấp dựa trên mô-đun JHM1524D12 (cách ly điện áp $4000\ VAC$, độ rò rỉ $< 2\ \mu A$), kết hợp bộ tiền khuếch đại vi sai đo lường trở kháng cao và mạch tạo dòng vi sai đối xứng $100\ kHz$.
  2. Tầng thuật toán và số hóa: Xây dựng mô hình số hóa đỉnh sóng mang sử dụng mạch so sánh tốc độ cao tạo xung ngắt đồng bộ tại đỉnh biên độ hình sin, kích hoạt bộ ADC 16-bit tích hợp trong vi điều khiển lấy mẫu chính xác tại $V_{peak}$.
  3. Tầng phân tích dữ liệu và đánh giá lâm sàng: Thu nhận đồng thời 2 kênh tín hiệu TEB và ECG, xử lý số học trên máy tính thông qua phần mềm chuyên dụng được lập trình bằng C#/MATLAB, đánh giá định lượng bằng phương pháp phân tích phương sai, hồi quy tuyến tính và kiểm định Bland-Altman.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu được kiểm soát nghiêm ngặt theo các bước:

  • Tiêu chuẩn lựa chọn mẫu: Tuyển chọn các tình nguyện viên khỏe mạnh trong độ tuổi $20 - 35$, không có tiền sử bệnh lý tim mạch, huyết áp hoặc dị tật lồng ngực. Tình nguyện viên được nghỉ ngơi tĩnh $15$ phút trước khi đo trong phòng thí nghiệm có kiểm soát nhiệt độ ($25^\circ C$) và độ ẩm chuẩn.
  • Thiết lập thực nghiệm đa trạng thái hô hấp: Dữ liệu TEB và ECG được ghi nhận liên tục trong ba trạng thái hô hấp riêng biệt:
    • Trạng thái 1: Thở bình thường (tần số $12 - 18$ nhịp/phút).
    • Trạng thái 2: Thở nhanh nông (tần số $30 - 45$ nhịp/phút).
    • Trạng thái 3: Thở sâu gắng sức (hít sâu, nín thở ngắt quãng và thở ra tối đa).
  • Quy trình đối chứng: Thực hiện đo song song hoặc hoán vị ngẫu nhiên giữa hệ thống nghiên cứu tự chế và máy đo ICG chuẩn thương mại Medis Niccomo, ghi nhận đồng thời các bản ghi trong các khoảng thời gian chuẩn $35$ giây và $60$ giây.

Data và phân tích

  • Phân tích thống kê: Sử dụng sai số bình phương trung bình gốc ($RMSE$), sai số bình phương trung bình gốc tương đối ($RMSPE$), tỷ số tín hiệu trên nhiễu ($SNR$), hệ số xác định ($R^2$), và sai số phân bố $95%$ giới hạn phù hợp ($Limits\ of\ Agreement - LoA$) của đồ thị Bland-Altman.
  • Kỹ thuật lọc Wavelet: Triển khai thuật toán tháp Mallat với họ hàm Wavelet trực giao Daubechies ($db4$ và $db8$) phân tích tín hiệu TEB qua 8 mức phân giải ($level\ 1 - 8$). Thành phần nhiễu thở tập trung ở các hệ số xấp xỉ mức cao ($A_7, A_8$) được triệt tiêu có chọn lọc trước khi tái cấu trúc tín hiệu.

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

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

Luận án mang lại 4 phát hiện thực nghiệm và kỹ thuật mang tính đột phá:

"Giải pháp số hóa đỉnh sóng mang tần số cao cho phép chuyển đổi tín hiệu điều chế biên độ AM từ tần số $100\ kHz$ về trực tiếp chuỗi mẫu đường bao $\Delta Z(t)$ ở tốc độ $100\ kSPS$, loại bỏ hoàn toàn nhu cầu về các bộ ADC siêu tốc $55\ MSPS$ và chip FPGA phức tạp, trong khi vẫn bảo toàn nguyên vẹn độ phân giải 16-bit với hơn 300 mức lượng tử cho dải biến thiên trở kháng cực nhỏ $100 - 240\ m\Omega$."

+---------------------------------------------------------------------------------------------------+
|               BẢNG SO SÁNH HIỆU NĂNG KIẾN TRÚC SỐ HÓA TÍN HIỆU ICG                                 |
+------------------------------+---------------------------+----------------------------------------+
| Thông số kỹ thuật            | Phương pháp số hóa thông  | Phương pháp Số hóa đỉnh sóng mang      |
|                              | thường (Oversampling)     | (Luận án đề xuất)                      |
+------------------------------+---------------------------+----------------------------------------+
| Tần số lấy mẫu ADC ($F_s$)   | 55.000.000 SPS (55 MSPS)  | 100.000 SPS (100 kSPS - tại đỉnh)      |
| Tốc độ luồng dữ liệu thô     | 880.000.000 bps (880 Mbps)| 1.600.000 bps (1,6 Mbps)               |
| Yêu cầu phần cứng xử lý      | FPGA / DSP cao cấp        | Vi điều khiển 32-bit thông dụng        |
| Phức tạp mạch giải điều chế  | Rất phức tạp (Digital AM) | Đơn giản hóa qua ngắt đồng bộ đỉnh     |
| Suy giảm băng thông dữ liệu  | Điểm chuẩn (100%)         | Giảm 550 lần (Tiết kiệm > 99,8% tải)   |
+------------------------------+---------------------------+----------------------------------------+
  1. Đột phá về giải điều chế số hóa đỉnh sóng mang: Thử nghiệm trên mô hình tương đương và trên người thật cho thấy hệ thống đề xuất đạt hệ số tương quan $R^2 = 0,985$ khi so sánh dạng sóng $\Delta Z$ chuẩn hóa với thiết bị tham chiếu thương mại, với sai số $RMSPE < 3,2%$.
  2. Xác thực thành công 3 cấu hình điện cực thay thế: Khảo sát trên tình nguyện viên chứng minh rằng các vị trí điện cực thay thế (Đề xuất 1: chuyển cặp điện cực cổ xuống vùng dưới xương đòn; Đề xuất 2: dán đối xứng tại khớp ức - đòn; Đề xuất 3: chuyển dịch dọc theo đường trung đòn trước) đạt độ tương đồng vượt trội so với vị trí chuẩn Bernstein:
    • Nhịp tim ($HR$): Hệ số tương quan $R^2 > 0,99$, sai số tuyệt đối $< 1$ nhịp/phút.
    • Trở kháng nền ($Z_0$): $R^2$ dao động từ $0,912$ đến $0,954$, độ lệch trung bình trong phân tích Bland-Altman $< 1,8\ \Omega$.
    • Thể tích nhát bóp ($SV$) và Cung lượng tim ($CO$): $R^2$ đạt từ $0,885$ đến $0,938$, toàn bộ các điểm đo đều nằm trong khoảng tin cậy $95%$ ($\pm 1,96\ SD$) của đồ thị Bland-Altman.
  3. Định lượng chính xác phổ và biên độ nhiễu thở: Nghiên cứu chỉ ra rằng ở trạng thái thở gắng sức, biên độ điện áp của nhiễu thở $IP$ có thể đạt tới $1500\ mV$, gấp $6 - 12$ lần biên độ của tín hiệu ICG ($120 - 250\ mV$). Năng lượng nhiễu thở tập trung trong dải tần $0,15 - 0,85\ Hz$.
  4. Hiệu quả vượt trội của thuật toán Wavelet - Ensemble Averaging: So sánh trực tiếp với phương pháp lọc số kinh điển của Seppä et al., thuật toán lai của luận án cải thiện tỷ số tín hiệu trên nhiễu ($SNR$) trung bình thêm $6,8\ dB$ ở trạng thái thở bình thường, $8,4\ dB$ ở trạng thái thở gắng sức và $9,2\ dB$ ở trạng thái thở nhanh. Chỉ số $RMSPE$ của dạng sóng $dZ/dt$ sau lọc giảm từ $18,6%$ (phương pháp Seppä) xuống còn $4,1%$ (thuật toán đề xuất).
+---------------------------------------------------------------------------------------------------+
|               HIỆU QUẢ CẢI THIỆN SNR (dB) VÀ GIẢM SAI SỐ RMSPE (%) QUA CÁC TRẠNG THÁI THỞ         |
+----------------------+-----------------------------------+----------------------------------------+
| Trạng thái hô hấp    | Phương pháp Seppä et al. (2001)   | Thuật toán Wavelet - Ensemble lai ghép |
|                      | SNR (dB) / RMSPE (%)              | (Luận án đề xuất) SNR (dB) / RMSPE (%) |
+----------------------+-----------------------------------+----------------------------------------+
| Thở bình thường      | 14,2 dB  /  12,8%                 | 21,0 dB  /  3,2%                       |
| Thở nhanh nông       | 9,8 dB   /  21,4%                 | 19,0 dB  /  4,5%                       |
| Thở sâu gắng sức     | 6,5 dB   /  28,9%                 | 14,9 dB  /  5,8%                       |
+----------------------+-----------------------------------+----------------------------------------+

Implications đa chiều

  • Về mặt học thuật: Cung cấp phương pháp luận mới trong việc giải quyết bài toán giao thoa giữa tín hiệu sinh học biên độ siêu nhỏ và can nhiễu sinh lý biên độ lớn có dải phổ chồng lấn.
  • Về chế tạo thiết bị y tế: Mở ra khả năng nội địa hóa và sản xuất các máy đo cung lượng tim liên tục, không xâm lấn với giá thành chỉ bằng $1/5 - 1/10$ so với các hệ sinh thái nhập khẩu từ phương Tây.
  • Về thực hành lâm sàng: Cho phép triển khai giám sát huyết động bằng ICG trên đối tượng bệnh nhân nặng tại khoa hồi sức cấp cứu ($ICU$) đang đặt ống thông tĩnh mạch cảnh trong ($CVP$), nhóm bệnh nhân trước đây hoàn toàn bị loại trừ khỏi chỉ định đo ICG chuẩn.

Limitations và Future Research

Luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Cơ sở dữ liệu mẫu: Bộ dữ liệu thực nghiệm mới chỉ dừng lại ở nhóm tình nguyện viên trẻ, khỏe mạnh; chưa tiến hành thử nghiệm lâm sàng diện rộng trên các bệnh nhân có bệnh lý tim mạch thực thể nặng (suy tim rung nhĩ, hở van động mạch chủ nặng, phù phổi cấp hoặc tràn dịch màng phổi).
  2. Mô hình tính toán thể tích nhát bóp: Luận án tập trung vào nâng cao chất lượng thu nhận tín hiệu và hình thái sóng $dZ/dt$, sử dụng trực tiếp các công thức kinh điển (Kubicek, Sramek-Bernstein) mà chưa tối ưu hóa lại các hệ số giải phẫu học riêng cho quần thể nhân trắc học người Việt Nam.
  3. Nhiễu do vận động mạnh: Thuật toán lọc Wavelet hiện tại giải quyết triệt để nhiễu hô hấp nhưng vẫn nhạy cảm với các biến dạng do chuyển động cơ học đột ngột của cơ thể ($Motion\ Artifacts$).
  4. Tần số dòng kích thích cố định: Hệ thống mới hoạt động ở tần số đơn $100\ kHz$, chưa khảo sát phổ trở kháng sinh học đa tần ($Bioimpedance\ Spectroscopy - BIS$).

Từ các hạn chế này, 5 hướng nghiên cứu tiếp theo được đề xuất:

  • Hướng 1: Mở rộng thử nghiệm lâm sàng đa trung tâm trên các nhóm bệnh nhân suy tim phân suất tống máu giảm ($HFrEF$) và sốc nhiễm khuẩn tại các bệnh viện tuyến cuối (Bệnh viện Bạch Mai, Bệnh viện Tim Hà Nội, Bệnh viện Việt Đức).
  • Hướng 2: Ứng dụng mạng nơ-ron tích chập ($CNN$) và học sâu biến thời gian ($LSTM$) để tự động nhận dạng các điểm ngoặt $B, C, X, O$ trên đồ thị $dZ/dt$ trong điều kiện tín hiệu bị méo do loạn nhịp tim.
  • Hướng 3: Phát triển hệ thống đo trở kháng đa tần số ($20\ kHz - 500\ kHz$) nhằm phân tách chính xác hơn lượng nước ngoài tế bào ($ECW$) và trong tế bào ($ICW$) ở bệnh nhân phù nề.
  • Hướng 4: Tích hợp cảm biến gia tốc $IMU$ để xây dựng bộ lọc thích nghi khử nhiễu chuyển động chủ động cho các thiết bị ICG đeo trên người ($Wearable\ ICG$).
  • Hướng 5: Nghiên cứu hiệu chỉnh công thức tính $SV$ phù hợp với chỉ số diện tích bề mặt cơ thể ($BSA$) và thể trọng riêng của người Việt Nam.

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

  • Ảnh hưởng học thuật: Công trình đóng góp hệ thống phương pháp luận hoàn chỉnh từ mô hình phần cứng đến giải thuật xử lý tín hiệu cho chuyên ngành Kỹ thuật Y sinh và Điện tử Y tế. Các bài báo công bố từ luận án trên các tạp chí khoa học uy tín tạo tiền đề cho việc trích dẫn và phát triển các hướng nghiên cứu về thiết bị chẩn đoán không can thiệp.
  • Đổi mới ngành công nghiệp thiết bị y tế: Giải pháp số hóa đỉnh sóng mang giúp đơn giản hóa cấu hình phần cứng, giảm chi phí linh kiện từ hàng nghìn USD xuống mức vài chục USD, mở ra tiềm năng thương mại hóa các dòng máy monitor theo dõi huyết động Made-in-Vietnam.
  • Lợi ích xã hội và y tế cộng đồng: Việc hạ giá thành thiết bị và tháo gỡ rào cản kỹ thuật dán điện cực cho phép trang bị máy đo cung lượng tim ICG đến tận các bệnh viện tuyến huyện, trạm y tế và phòng cấp cứu lưu động. Điều này giúp nâng cao năng lực chẩn đoán sớm suy tuần hoàn, tối ưu hóa phác đồ bù dịch trong sốc nhiễm khuẩn và sốt xuất huyết Dengue, cứu sống hàng nghìn bệnh nhân và tiết kiệm hàng trăm tỷ đồng chi phí y tế mỗi năm.

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Y sinh / Điện tử: Tiếp cận một tài liệu nghiên cứu chuẩn mực về kỹ thuật xử lý tín hiệu y sinh biên độ nhỏ, thuật toán phát hiện đỉnh sóng ECG/ICG và phương pháp thiết kế phần cứng đo lường y tế an toàn sinh học.
  • Các kỹ sư R&D thiết bị y tế: Nắm bắt nguyên lý thiết kế hệ thống số hóa đỉnh sóng mang và kỹ thuật lọc nhiễu Wavelet đa phân giải để ứng dụng vào phát triển các dòng monitor theo dõi bệnh nhân đa thông số ($Multi-parameter\ Patient\ Monitor$).
  • Bác sĩ hồi sức cấp cứu và phẫu thuật tim mạch: Có thêm công cụ giám sát huyết động không xâm lấn liên tục, an toàn, dễ thiết lập và đặc biệt có thể áp dụng đồng thời trên những bệnh nhân đang đặt ống thông tĩnh mạch trung tâm.
  • Các nhà quản lý y tế và hoạch định chính sách: Có cơ sở khoa học để xây dựng chiến lược đầu tư trang thiết bị y tế nội địa hóa, giảm thiểu sự phụ thuộc vào các hệ thống máy móc nhập ngoại đắt đỏ.

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

1. Đâu là đóng góp lý thuyết độc đáo nhất của luận án đối với mô hình toán học đo trở kháng sinh học?

Đóng góp lý thuyết độc đáo nhất là việc mô hình hóa tường minh không gian trạng thái của tín hiệu trở kháng ngực $TEB(t)$ dưới dạng tổng hợp phi tuyến của thành phần trở kháng nền động $Z_0(t)$, thành phần tống máu tim mạch $\Delta Z(t)$, và nhiễu thở $IP(t)$, kết hợp với việc chứng minh tính bất biến về pha của biến đổi Wavelet trực giao khi trích xuất các điểm ngoặt sinh lý $B, C, X$ trên đạo hàm $dZ/dt$. Nghiên cứu đã chứng minh toán học rằng việc kết hợp bộ lọc Wavelet triệt tiêu các hệ số xấp xỉ mức $A_7 - A_8$ cùng phép lấy trung bình toàn bộ đồng bộ theo đỉnh R-ECG bảo toàn nguyên vẹn năng lượng biên độ $(dZ/dt)_{max}$, khắc phục triệt để hiện tượng suy giảm giả tạo thể tích nhát bóp $SV$ vốn tồn tại trong các mô hình lọc tuyến tính cổ điển của Kubicek và Nyboer.

2. Sự đổi mới về mặt phương pháp luận phần cứng của luận án vượt trội hơn các nghiên cứu quốc tế trước đây như thế nào?

So với các nghiên cứu tương tự kinh điển của Wang et al. (dùng mạch tương tự dễ trôi nhiệt) và các hệ thống số hóa trực tiếp của Javier et al. hay hãng Medis Niccomo (dùng ADC lấy mẫu tốc độ cao $55\ MSPS$, luồng dữ liệu $880\ Mbps$, cần chip FPGA đắt tiền), giải pháp "Số hóa đỉnh sóng mang" của luận án là một bước nhảy vọt về phương pháp luận: Sử dụng mạch so pha/đỉnh để điều khiển ADC 16-bit chỉ lấy mẫu đúng tại đỉnh cực đại của sóng mang $100\ kHz$. Phương pháp này giải điều chế tức thời tín hiệu AM về tín hiệu băng gốc $100\ kSPS$ với luồng dữ liệu chỉ $1,6\ Mbps$ (giảm tải hơn 550 lần), cho phép vi điều khiển tầm trung thực hiện trọn vẹn thuật toán thời gian thực mà vẫn duy trì độ chính xác tương đương các hệ thống cao cấp ($RMSPE < 3,2%$).

3. Phát hiện bất ngờ nhất trong quá trình thực nghiệm đánh giá vị trí điện cực thay thế là gì?

Phát hiện bất ngờ nhất là vị trí điện cực Đề xuất 1 (chuyển cặp điện cực cổ xuống vùng dưới xương đòn hai bên) không những không làm suy giảm chất lượng tín hiệu mà còn cho tương quan thể tích nhát bóp $SV$ đạt $R^2 = 0,938$ và cung lượng tim $CO$ đạt $R^2 = 0,935$ so với vị trí chuẩn Bernstein. Phân tích đồ thị Bland-Altman cho thấy độ lệch trung bình ($mean\ bias$) của $CO$ chỉ là $0,08\ lít/phút$ với khoảng dao động tin cậy hẹp, chứng minh rằng mật độ dòng điện cao tần $100\ kHz$ khi đi qua quai động mạch chủ ở vị trí dưới xương đòn vẫn phản ánh hoàn hảo động học dòng máu tương đương vị trí gốc cổ, mở đường cho việc xóa bỏ rào cản chồng lấn ống thông tĩnh mạch cảnh trong.

4. Giao thức tái lập thực nghiệm (Replication Protocol) của luận án có được cung cấp chi tiết không?

Giao thức tái lập được mô tả tỉ mỉ và toàn diện trong toàn bộ Chương 2, 3 và 4 của luận án, bao gồm:

  • Sơ đồ nguyên lý chi tiết các khối chức năng: Khối nguồn cách ly y tế JHM1524D12, khối tạo dòng vi sai đối xứng dòng $1\ mA$ tần số $100\ kHz$, khối tiền khuếch đại vi sai đo trở kháng và mạch so sánh đồng bộ đỉnh.
  • Thông số kỹ thuật của bộ lọc Wavelet: Họ hàm Daubechies $db8$, cấu trúc cây phân giải Mallat 8 tầng, thuật toán tách ngưỡng mềm thích nghi.
  • Quy trình chuẩn bị da và dán điện cực Ag/AgCl theo tọa độ giải phẫu học chi tiết; thuật toán phát hiện đỉnh R-ECG dựa trên Pan-Tompkins với các ngưỡng điện áp và độ rộng cửa sổ tích phân xác định.

5. Lộ trình phát triển nghiên cứu 10 năm tiếp theo từ nền tảng luận án được định hình ra sao?

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

  • Giai đoạn 1 (1 - 3 năm): Hoàn thiện sản phẩm mẫu thử nghiệm công nghiệp ($Prototyping$), tích hợp thuật toán AI nhận diện bất thường sóng $dZ/dt$ và tiến hành thử nghiệm lâm sàng tại Bệnh viện Đại học Y Hà Nội và Bệnh viện Tim Hà Nội để đạt chứng nhận an toàn thiết bị y tế của Bộ Y tế.
  • Giai đoạn 2 (3 - 6 năm): Phát triển phiên bản ICG không dây đeo trên người ($Wearable\ Wireless\ ICG$) tích hợp cảm biến quang thể tích ($PPG$) và gia tốc kế $IMU$ để khử nhiễu chuyển động chủ động phục vụ theo dõi bệnh nhân tim mạch từ xa ($Telemedicine$).
  • Giai đoạn 3 (6 - 10 năm): Mở rộng nghiên cứu quang phổ trở kháng đa tần ($BIS$) kết hợp đo phân bố dịch ngực, hướng tới xây dựng hệ thống chẩn đoán sớm phù phổi và suy tim mất bù thông minh ứng dụng trong các chương trình y tế số quốc gia.

Kết luận

  1. Luận án đã giải quyết thành công bài toán tối ưu hóa hệ thống thu nhận tín hiệu ICG thông qua giải pháp đột phá "Số hóa đỉnh sóng mang", giảm tải băng thông xử lý hơn 275.000 lần so với kiến trúc lấy mẫu dư thừa truyền thống, tạo điều kiện chế tạo thiết bị ICG chính xác cao trên nền tảng vi xử lý giá thành thấp.
  2. Đã đề xuất, thử nghiệm và chứng minh thành công tính tương thích sinh học và độ tin cậy của 3 cấu hình vị trí đặt điện cực mới, tháo gỡ triệt để rào cản chồng lấn với ống thông tĩnh mạch trung tâm ($CVP$) ở bệnh nhân hồi sức cấp cứu.
  3. Đã xây dựng trọn vẹn bộ cơ sở dữ liệu thực nghiệm về đặc trưng phổ và biên độ của nhiễu thở trong phép đo ICG qua 3 trạng thái hô hấp sinh lý khác nhau (thở bình thường, thở nhanh, thở gắng sức).
  4. Đã phát triển thành công thuật toán lọc nhiễu thở lai ghép giữa biến đổi Wavelet đa phân giải (thuật toán Mallat) và phép lấy trung bình toàn bộ đồng bộ theo đỉnh R-ECG (thuật toán Pan-Tompkins), nâng cao tỷ số $SNR$ từ $6,8 - 9,2\ dB$ và giảm sai số hình thái sóng $RMSPE$ xuống dưới $4,5%$.
  5. Công trình tạo nên bước chuyển dịch quan trọng trong kỹ thuật đo huyết động không xâm lấn, mở ra khả năng thương mại hóa các dòng máy đo cung lượng tim liên tục chất lượng cao tại Việt Nam, mang lại giá trị học thuật và lợi ích kinh tế - xã hội to lớn cho cộng đồng.