Tổng quan nghiên cứu

Bộ não con người là cơ quan điều khiển trung tâm tinh vi nhất với khối lượng trung bình khoảng 1,4 kg, chứa gần 1.000 tỷ tế bào thần kinh phối hợp chặt chẽ để điều hành toàn bộ hành vi thể chất và nhận thức. Trong y học hiện đại và kỹ thuật y sinh, việc tìm hiểu sự liên kết giữa các vùng chức năng trên vỏ não và hệ thống cơ bắp đóng vai trò then chốt trong chẩn đoán cũng như phục hồi chức năng. Vấn đề nghiên cứu trọng tâm của luận văn là xác định mối quan hệ tương hỗ giữa sự mất cân bằng cấu trúc cơ xương hàm (độ lệch quai hàm) và khả năng phát lực của hệ vận động chi trên (lực cánh tay). Trong thực tế, các tổn thương lệch khớp cắn thường dẫn đến sự suy giảm sức mạnh cơ bắp đối bên, nhưng cơ chế đáp ứng huyết động học trên vỏ não vận động vẫn chưa được định lượng rõ ràng bằng các phương pháp đo lường không xâm lấn tiện dụng.

Mục tiêu cụ thể của đề tài là ứng dụng kỹ thuật quang phổ cận hồng ngoại chức năng (fNIRS) 24 kênh để ghi nhận sự biến thiên nồng độ Oxy-Hemoglobin (Oxy-Hb), từ đó xây dựng giải thuật lọc nhiễu, phân tích đặc trưng tần số và hồi quy diện tích nhằm định lượng chính xác tác động của độ lệch quai hàm lên vùng vận động não bộ. Nghiên cứu được triển khai thực nghiệm tại Thành phố Hồ Chí Minh trên nhóm đối tượng thanh niên khỏe mạnh trong khoảng thời gian 50,05 giây cho mỗi chu kỳ vận động kết hợp giữa cắn miếng đệm (độ dày 0 mm, 2 mm, 3 mm) và nâng các mức tạ (từ 0 kg đến 7 kg). Ý nghĩa của nghiên cứu thể hiện qua việc giảm thiểu sai số nhận dạng tín hiệu vận động xuống dưới 5%, đồng thời cung cấp mô hình toán học thực nghiệm giúp các bác sĩ chuyên khoa và chuyên gia giao diện não - máy tính (BCI) có công cụ định lượng chuẩn xác phản ứng thần kinh trung ương.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên 2 nền tảng lý thuyết cốt lõi trong kỹ thuật y sinh và xử lý tín hiệu sinh học:

Lý thuyết quang học mô sinh học và Định luật Beer-Lambert mở rộng: Kỹ thuật fNIRS sử dụng ánh sáng trong dải bước sóng cận hồng ngoại từ 650 nm đến 950 nm có khả năng xuyên qua mô sọ đến độ sâu khoảng 3 cm tại vỏ não. Dựa trên sự hấp thụ quang phổ đặc trưng của Oxy-Hemoglobin (Oxy-Hb) và Deoxy-Hemoglobin (deOxy-Hb), mật độ quang (Optical Density) được tính toán để phản ánh mức độ tập trung huyết động học khi tế bào thần kinh hoạt động.

Lý thuyết liên kết thần kinh - mạch máu (Neurovascular Coupling): Khi một vùng não vận động bị kích thích bởi hoạt động thể chất, lưu lượng máu cục bộ tăng lên nhanh chóng để cung cấp oxy, dẫn đến sự gia tăng đột biến của Oxy-Hb và suy giảm tương ứng của deOxy-Hb.

Các khái niệm chính được sử dụng xuyên suốt bao gồm: Vỏ não vận động (Motor Cortex) tại bán cầu não trái; Giao diện Não - Máy tính (NIRS-BCI); Phân tích đa phân giải Wavelet rời rạc (MRA); Bộ lọc đa thức Savitzky-Golay và Tích phân diện tích đáp ứng huyết động.

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

Nguồn dữ liệu thực nghiệm được thu thập trực tiếp từ hệ thống máy quang phổ cận hồng ngoại chức năng FOIRE-3000 do hãng Shimadzu (Nhật Bản) sản xuất, vận hành ở 3 bước sóng chuyên biệt là 780 nm, 805 nm và 830 nm. Thiết bị được bố trí ma trận đo 4x4 bao gồm 8 nguồn phát quang và 8 đầu thu tạo thành 24 kênh tín hiệu bao phủ toàn bộ vùng điều khiển chuyển động ở bán cầu não trái.

Cỡ mẫu nghiên cứu gồm 5 đối tượng khỏe mạnh được tinh lọc kỹ lưỡng từ 10 ứng viên ban đầu thông qua phương pháp chọn mẫu có chủ đích, loại trừ các đối tượng không tuân thủ nghiêm ngặt quy trình đo hoặc bị nhiễu động do căng thẳng cơ học. Lý do lựa chọn cỡ mẫu này là nhằm bảo đảm tính đồng nhất sinh lý tối đa trong điều kiện thí nghiệm kiểm soát nghiêm ngặt.

Giao thức thí nghiệm khối (Block Design) được thiết lập với timeline 50,05 giây mỗi phiên: 20 giây đầu nghỉ ngơi thư giãn, 10 giây tiếp theo thực hiện đồng thời cắn miếng đệm và nâng tạ, 20 giây cuối cùng nghỉ ngơi để tín hiệu trở về trạng thái nền. Chu kỳ lấy mẫu đạt 0,130 giây tương ứng 386 mẫu dữ liệu cho mỗi kênh.

Phương pháp phân tích được lựa chọn gồm 3 bước logic:

  • Tiền xử lý: Sử dụng bộ lọc làm mịn đa thức Savitzky-Golay (kích thước khung cửa sổ n = 11, bậc đa thức M = 3) để loại bỏ nhiễu sinh lý tần số cao (nhịp tim, nhịp thở) mà vẫn bảo toàn đỉnh sóng huyết động.
  • Trích xuất đặc trưng: Áp dụng biến đổi Wavelet rời rạc (DWT) với họ Wavelet Biorthogonal (bior5.5) tại mức phân giải 6 kết hợp thuật toán ngưỡng thích nghi nhằm chọn lọc 6 đến 10 kênh tích cực nhất trong 24 kênh đo.
  • Định lượng tương quan: Thực hiện hồi quy đa thức bậc cao (từ bậc 3 đến bậc 20) trên tín hiệu Oxy-Hb làm phẳng và tính tích phân diện tích vùng vận động để so sánh định lượng.

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

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

Phân tích dữ liệu thực nghiệm từ 24 kênh quang phổ cận hồng ngoại đã đem lại 4 phát hiện khoa học quan trọng:

Hiệu quả lọc nhiễu vượt trội của cấu hình Savitzky-Golay: Với cửa sổ n = 11 và bậc đa thức M = 3, thuật toán đã triệt tiêu hơn 85% các thành phần gai nhiễu đột biến, duy trì dạng sóng mượt mà của Oxy-Hb trong khoảng thời gian kích hoạt 10 giây tốt hơn hẳn so với các cấu hình cửa sổ n = 17 hoặc n = 7.

Xác định chính xác kênh kích hoạt vận động qua Wavelet bior5.5: Biến đổi Wavelet tại mức phân tích thứ 6 với họ bior5.5 cho độ ổn định dạng sóng cao nhất so với họ Daubechies (db13) và Symlets (sym7), giúp phân tách rõ ràng dải tần số thấp mang thông tin huyết động và tự động định vị từ 6 đến 10 kênh tích cực trong khu vực vỏ não vận động.

Độ trễ sinh học của tín hiệu Oxy-Hb: Kết quả đo thực tế cho thấy nồng độ Oxy-Hb không tăng ngay lập tức tại giây thứ 20 mà có độ trễ sinh học rõ rệt. Vùng đáp ứng huyết động thực tế xuất hiện tập trung trong khoảng từ mẫu thứ 151 đến mẫu thứ 280 (tương đương mốc thời gian từ 19,58 giây đến 36,31 giây), phản ánh đúng thời gian trễ của dòng máu cục bộ.

Mối tương quan tuyến tính nghịch giữa độ lệch quai hàm và năng lượng kích hoạt não bộ: Khi nâng mức tải từ 4 kg lên 7 kg ở trạng thái bình thường (không cắn đệm), diện tích tích phân Oxy-Hb tăng dần từ 15% đến 35%. Tuy nhiên, khi chèn miếng đệm gây lệch quai hàm 2 mm và 3 mm, diện tích đáp ứng huyết động vùng vận động bị suy giảm tương ứng từ 18% đến 32%, kéo theo sự suy giảm rõ rệt lực nâng của cánh tay.

Thảo luận kết quả

Hiện tượng suy giảm diện tích Oxy-Hb khi quai hàm bị lệch có thể giải thích qua cơ chế tương tác ức chế thần kinh trung ương. Dây thần kinh sinh ba (Trigeminal nerve) từ cơ hàm truyền tín hiệu căng thẳng bất thường về thân não, gây cản trở đường dẫn truyền vỏ não - tủy sống (Corticospinal tract) điều khiển cánh tay đối bên. Khi đối tượng cắn miếng đệm không cân bằng, vỏ não vận động phải phân tán tài nguyên để điều chỉnh cân bằng cơ mặt, làm giảm lưu lượng máu và mức tiêu thụ oxy tập trung cho nhiệm vụ nâng tạ.

So sánh với các công trình nghiên cứu sử dụng tín hiệu điện cơ (EMG) trước đây của các tác giả quốc tế, kết quả của luận văn hoàn toàn tương đồng về mặt kết luận sinh cơ học nhưng vượt trội hơn nhờ cung cấp bằng chứng huyết động học trực tiếp từ bên trong vỏ não mà không bị ảnh hưởng bởi mồ hôi hay nhiễu điện cực da.

Dữ liệu nghiên cứu được biểu diễn trực quan và sinh động thông qua các đồ thị phân tách đa mức Wavelet và bảng ma trận diện tích đường cong hồi quy đa thức. Khi quan sát biểu đồ không gian 2D của 24 kênh đo, các kênh trung tâm như kênh 1, kênh 2, kênh 3, kênh 22, kênh 23 và kênh 24 thể hiện rõ nét sự co cụm của biên độ Oxy-Hb khi độ dày miếng đệm tăng từ 0 mm lên 3 mm.

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

Dựa trên các kết quả thực nghiệm và thuật toán đã xây dựng, luận văn đưa ra 4 khuyến nghị hành động cụ thể:

Thứ nhất, tích hợp module lọc Savitzky-Golay (n = 11, M = 3) và biến đổi Wavelet bior5.5 vào các phần mềm phân tích fNIRS trực tuyến. Đội ngũ kỹ sư y sinh cần tối ưu hóa mã nguồn thuật toán nhằm giảm độ trễ xử lý dữ liệu xuống dưới 50 mili-giây, hoàn thành triển khai thử nghiệm trong vòng 6 tháng tới.

Thứ hai, ứng dụng phương pháp đo diện tích Oxy-Hb vào quy trình chẩn đoán lâm sàng tại các bệnh viện Răng Hàm Mặt và trung tâm chỉnh hình. Bác sĩ chuyên khoa cần đưa chỉ số fNIRS vào đánh giá trước và sau phẫu thuật điều chỉnh khớp cắn, hướng tới mục tiêu phục hồi 90% cân bằng lực cơ cánh tay cho bệnh nhân sau 3 tháng can thiệp.

Thứ ba, phát triển các hệ thống giao diện Não - Máy tính (BCI) điều khiển chi giả dựa trên tín hiệu quang phổ hồng ngoại. Các nhóm nghiên cứu BCI cần áp dụng thuật toán hồi quy đa thức kết hợp ngưỡng thích nghi để tăng độ chính xác nhận dạng lệnh vận động đạt mức trên 95% trong lộ trình 12 tháng.

Thứ tư, mở rộng giao thức nghiên cứu đa trung tâm với cỡ mẫu từ 50 đến 100 người bao gồm cả đối tượng bệnh lý lệch khớp cắn thực tế và vận động viên chuyên nghiệp. Các trường đại học khối kỹ thuật y sinh cần chủ trì liên kết với các viện phục hồi chức năng để thiết lập bộ chỉ số chuẩn quốc gia trong vòng 18 tháng tới.

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

Nhà nghiên cứu và học viên cao học ngành Kỹ thuật Y sinh, Kỹ thuật Điện - Điện tử: Nắm bắt phương pháp xử lý tín hiệu không xâm lấn fNIRS, kỹ thuật lọc thích nghi đa thức và giải thuật phân tách đa phân giải Wavelet ứng dụng trên tín hiệu huyết động học thực tế.

Bác sĩ chuyên khoa Răng Hàm Mặt, Cơ Xương Khớp và Phục hồi Chức năng: Tham khảo cơ sở sinh lý thần kinh và bằng chứng định lượng về sự liên kết giữa sai lệch khớp cắn với suy giảm trương lực cơ ngoại vi để ứng dụng vào phác đồ điều trị chỉnh nha toàn diện.

Kỹ sư phát triển hệ thống Y tế thông minh và Giao diện Não - Máy tính (BCI): Khai thác cấu trúc thuật toán đặt ngưỡng và hồi quy đa thức 24 kênh để thiết kế các bộ giải mã tín hiệu thần kinh thời gian thực cho thiết bị hỗ trợ người khuyết tật.

Huấn luyện viên thể thao đỉnh cao và chuyên gia Khoa học Vận động: Áp dụng kiến thức về cân bằng khớp cắn để tối ưu hóa tư thế hàm cho vận động viên cử tạ, bắn súng và võ thuật nhằm tối đa hóa hiệu suất sinh lực chi trên.

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

Kỹ thuật quang phổ cận hồng ngoại (fNIRS) có ưu thế gì so với điện não đồ (EEG) và cộng hưởng từ chức năng (fMRI)? fNIRS cung cấp độ phân giải không gian vượt trội so với EEG và ít bị nhiễu do cử động cơ mặt khi đối tượng thực hiện động tác cắn hoặc nâng tải. So với fMRI, máy fNIRS FOIRE-3000 nhỏ gọn, chi phí vận hành thấp, cho phép đối tượng ngồi thực hiện vận động tự nhiên trong thời gian thực.

Tại sao bộ lọc Savitzky-Golay với cửa sổ 11 và bậc 3 lại được lựa chọn để làm phẳng tín hiệu Oxy-Hb? Bộ lọc Savitzky-Golay bảo toàn biên độ cực đại và độ rộng của đỉnh sóng huyết động tốt hơn các bộ lọc trung bình trượt thông thường. Thực nghiệm trên 386 mẫu cho thấy cấu hình n = 11, M = 3 loại bỏ triệt để nhiễu tần số cao mà không làm méo mó các đặc trưng vận động.

Họ Wavelet Biorthogonal (bior5.5) đóng vai trò gì trong việc xác định vùng vận động não bộ? Họ bior5.5 sở hữu tính đối xứng và đặc tính pha tuyến tính giúp tái tạo dạng sóng Oxy-Hb mượt mà và ổn định nhất ở mức phân giải 6. Thuật toán giúp trích xuất chính xác các thành phần xấp xỉ tần số thấp, hỗ trợ lọc bỏ các kênh không hoạt động trong 24 kênh đo.

Lệch quai hàm 3 mm làm suy giảm lực cánh tay như thế nào theo kết quả đo tín hiệu não? Độ lệch quai hàm 3 mm gây ức chế phản xạ thần kinh điều khiển, làm giảm từ 18% đến 32% diện tích tích phân Oxy-Hb tại vùng vận động bán cầu não trái. Sự suy giảm lưu lượng máu này tương ứng trực tiếp với hiện tượng mất sức mạnh cơ bắp cánh tay khi nâng mức tạ 4 kg đến 7 kg.

Ý nghĩa của việc tính diện tích đường cong hồi quy Oxy-Hb là gì? Diện tích tích phân phản ánh tổng năng lượng và mức tiêu thụ oxy của các nơ-ron tại vùng vỏ não vận động trong suốt chu kỳ kích hoạt 10 giây. Chỉ số diện tích giúp lượng hóa chính xác mức độ tham gia của não bộ thay vì chỉ dựa vào biên độ tức thời dễ bị sai lệch.

Kết luận

  • Luận văn đã xây dựng thành công quy trình xử lý tín hiệu fNIRS 24 kênh hoàn chỉnh, kết hợp bộ lọc Savitzky-Golay và biến đổi Wavelet bior5.5 ở mức phân tích thứ 6.
  • Định vị chính xác vùng vỏ não vận động tích cực tại bán cầu trái với khoảng 6 đến 10 kênh trọng điểm chịu trách nhiệm điều phối hoạt động cắn và nâng tải.
  • Chứng minh định lượng sự suy giảm diện tích nồng độ Oxy-Hb từ 18% đến 32% khi độ lệch quai hàm tăng từ 0 mm lên 3 mm trong các bài kiểm tra nâng tải từ 4 kg đến 7 kg.
  • Đóng góp giải thuật hồi quy đa thức thích nghi giúp nâng cao độ tin cậy trong phân tích tín hiệu huyết động học thần kinh không xâm lấn.
  • Mở ra hướng nghiên cứu liên ngành giá trị cao giữa chẩn đoán chỉnh hình nha khoa, kỹ thuật y sinh và hệ thống giao diện não - máy tính BCI.

Trong giai đoạn 6 đến 12 tháng tiếp theo, nhóm nghiên cứu dự kiến sẽ hoàn thiện thuật toán xử lý dữ liệu trực tuyến thời gian thực và tích hợp mạng nơ-ron nhân tạo (ANN) nhằm tự động hóa quy trình phân loại trạng thái vận động. Đề tài mở ra tiềm năng hợp tác sâu rộng với các cơ sở y tế chuyên khoa để thử nghiệm lâm sàng trên quy mô lớn, góp phần nâng cao chất lượng điều trị và chăm sóc sức khỏe cộng đồng.