Tổng quan nghiên cứu

Trong bối cảnh trí tuệ nhân tạo và học máy phát triển vượt bậc, việc giải quyết các bài toán mô hình hóa hiện tượng tự nhiên đặt ra thách thức lớn về chi phí thu thập dữ liệu và năng lực tính toán. Bài toán xác định hệ số của phương trình đạo hàm riêng (PDEs), hay còn gọi là bài toán ngược, giữ vai trò then chốt trong nhiều lĩnh vực kỹ thuật, vật lý và định giá tài chính. Tuy nhiên, các phương pháp số truyền thống thường gặp khó khăn nghiêm trọng khi dữ liệu quan sát bị phân tán, có nhiễu hoặc khi không gian trạng thái mở rộng lên đến 100 chiều.

Luận văn thạc sĩ khoa học chuyên ngành Toán Tin thực hiện tại Viện Toán ứng dụng và Tin học thuộc Trường Đại học Bách khoa Hà Nội vào năm 2022 tập trung giải quyết bài toán ngược này thông qua việc nghiên cứu mạng học sâu và mạng nơ-ron thông tin vật lý (Physics-Informed Neural Networks - PINNs). Mục tiêu cốt lõi của nghiên cứu là xây dựng các thuật toán tối ưu cho phép ước lượng chính xác các tham số chưa biết của phương trình PDEs từ một lượng rất hạn chế các điểm quan sát thực nghiệm.

Phạm vi nghiên cứu bao quát từ các phương trình elliptic, parabolic kinh điển đến phương trình Navier-Stokes phi tuyến trong cơ học chất lưu và phương trình Black-Scholes 100 chiều có tính đến rủi ro vỡ nợ trong tài chính định lượng. Kết quả nghiên cứu mang lại ý nghĩa học thuật và ứng dụng to lớn khi giúp giảm thời gian tính toán định giá phái sinh đa chiều từ 140 phút xuống chỉ còn 36 phút (tiết kiệm khoảng 74,3% thời gian xử lý), đồng thời kiểm soát sai số tương đối ở mức dưới 0,2%, mở ra hướng tiếp cận đột phá cho lĩnh vực tính toán khoa học hiện đại.

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 nền tảng lý thuyết mạng nơ-ron nhân tạo và định lý xấp xỉ phổ quát do Kurt Hornik thiết lập năm 1991, sau đó được mở rộng cho mạng nơ-ron sâu bởi Zhou Lu vào năm 2017. Trọng tâm của khung lý thuyết là mô hình mạng nơ-ron thông tin vật lý do Maziar Raissi và cộng sự đề xuất năm 2017. Khác với học sâu truyền thống thuần túy hướng dữ liệu, mạng PINNs nhúng trực tiếp các định luật bảo toàn và toán tử đạo hàm vào cấu trúc huấn luyện.

Mô hình nghiên cứu triển khai 4 khái niệm cốt lõi:

  1. Hàm tổn thất tích hợp phần dư vật lý (Physics Residual Loss): Đóng vai trò ràng buộc đầu ra mạng phải thỏa mãn phương trình đạo hàm riêng tại các điểm collocation trong miền không-thời gian.
  2. Phương pháp vi phân tự động (Automatic Differentiation): Tính toán chính xác các đạo hàm riêng bậc 1 và bậc 2 của nghiệm xấp xỉ mà không gặp sai số làm tròn như phương pháp sai phân cổ điển.
  3. Phương trình vi phân ngẫu nhiên ngược (BSDEs) và biểu diễn Feynman-Kac phi tuyến: Thiết lập cầu nối chuyển đổi bài toán biên parabolic phi tuyến nhiều chiều sang hệ phương trình ngẫu nhiên có thể huấn luyện bằng mạng học sâu.
  4. Mạng nơ-ron dư (ResNet) và kỹ thuật điều chuẩn L2: Tăng cường độ ổn định cho quá trình huấn luyện sâu và ngăn chặn hiện tượng quá khớp (overfitting).

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng chuẩn hóa từ các công thức nghiệm giải tích và sơ đồ rời rạc hóa ngẫu nhiên Euler-Maruyama. Cụ thể, cỡ mẫu nghiên cứu bao gồm 1000 điểm dữ liệu lưới chia đều trên miền không gian hai chiều đối với phương trình Poisson và Navier-Stokes; đồng thời sinh ngẫu nhiên 64 minibatch mẫu tại mỗi bước lặp với 50 bước thời gian đối với mô hình Black-Scholes.

Phương pháp chọn mẫu áp dụng kỹ thuật lấy mẫu phân bố đều trên miền xác định và trên biên đối với các bài toán truyền nhiệt dừng, kết hợp với phương pháp lấy mẫu ngẫu nhiên chuyển động Brown đối với các bài toán tài chính đa chiều.

Lý do lựa chọn phương pháp phân tích tối ưu hóa kết hợp giữa Adam và L-BFGS là nhằm khai thác triệt để thế mạnh của từng thuật toán: thuật toán Adam (bậc 1) giúp thoát khỏi các điểm cực trị cục bộ nhanh chóng trong giai đoạn đầu, trong khi thuật toán L-BFGS (bậc 2 dựa trên ma trận Quasi-Newton) giúp quá trình tối ưu hội tụ cực nhanh về nghiệm chính xác ở giai đoạn tinh chỉnh. Toàn bộ quy trình thực nghiệm được triển khai lập trình trên nền tảng Python 3 và môi trường tính toán Google Colab trong suốt timeline thực hiện luận văn từ năm 2021 đến tháng 7 năm 2022.

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

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

Quá trình thực nghiệm số trên 5 lớp bài toán từ đơn giản đến phức tạp đã chứng minh năng lực vượt trội của mô hình:

  1. Ước lượng chính xác tham số phương trình Poisson và Navier-Stokes: Mạng học sâu xác định tham số khuếch tán của phương trình Poisson đạt giá trị xấp xỉ 1,05 (so với giá trị thực nghiệm 1,0) với hàm tổn thất giảm còn 0,0032. Đối với phương trình Navier-Stokes 2D phụ thuộc thời gian, hệ số nhớt được ước lượng chính xác tuyệt đối ở mức 0,01 với hàm tổn thất giảm sâu xuống mức 0,00006 và sai số chuẩn L2 của vận tốc đạt mức xấp xỉ 0,0000001.
  2. Tối ưu hóa xấp xỉ phương trình Black-Scholes 1 chiều: Mô hình đạt sai số tương đối chỉ 0,2% sau 3000 bước lặp trong thời gian 13 phút. Đặc biệt, việc ứng dụng kiến trúc mạng ResNet giúp khắc phục triệt để hiện tượng kẹt nghiệm tại giá trị 1,686 của mạng truyền thẳng cổ điển, đưa nghiệm hội tụ chuẩn xác về giá trị 1,450.
  3. Đột phá trong giải phương trình Black-Scholes phi tuyến 100 chiều: Khi định giá xác nhận quyền sở hữu kiểu châu Âu trên 100 tài sản cơ sở có rủi ro vỡ nợ, mô hình đạt sai số tương đối 0,05% sau 6000 bước lặp với thời gian chạy chỉ 36 phút, tìm ra giá trị định giá chuẩn xác 57,300.
  4. Xử lý thành công bài toán biên di động Stefan: Luận văn đã xác định đồng thời cả trường nhiệt độ và hàm vị trí biên chuyển động với sai số kiểm thử L2err xấp xỉ 0,0001 thông qua kiến trúc mạng nơ-ron kép song song.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương pháp đạt độ chính xác cao là nhờ việc nhúng trực tiếp toán tử vi phân vào hàm tổn thất, tạo ra cơ chế tự hiệu chỉnh vật lý mạnh mẽ. Mô hình không yêu cầu hàng triệu mẫu dữ liệu như học sâu truyền thống mà vẫn có thể khái quát hóa hoàn hảo từ vài chục đến vài trăm điểm quan sát.

Khi so sánh với các phương pháp giải số truyền thống, phương pháp mạng học sâu chứng minh tính ưu việt rõ rệt. Đối với bài toán định giá tài sản 100 chiều, phương pháp lặp Picard truyền thống phải mất đến 140 phút mới tìm ra kết quả, trong khi mạng học sâu chỉ cần 36 phút, giúp tiết kiệm tới 74,3% thời gian tính toán. Hơn nữa, nếu bỏ qua yếu tố rủi ro vỡ nợ, mô hình tuyến tính đơn giản sẽ định giá sai lệch nghiêm trọng lên mức 60,616, chứng tỏ tầm quan trọng của việc giải phương trình phi tuyến nhiều chiều.

Các kết quả nghiên cứu này có thể được trình bày một cách trực quan thông qua các biểu đồ đường thể hiện tốc độ suy giảm hàm tổn thất theo số vòng lặp từ bước 0 đến bước 6000, biểu đồ cột so sánh thời gian thực thi giữa phương pháp học sâu và phương pháp Picard, cùng bảng đối chiếu sai số tương đối giữa kiến trúc Feed-Forward thuần túy và kiến trúc ResNet.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 khuyến nghị hành động thiết thực:

  1. Tích hợp giải thuật PINNs vào hệ thống định giá tài chính tự động: Khối Quản trị rủi ro tại các ngân hàng thương mại và quỹ đầu tư cần áp dụng mô hình mạng học sâu để định giá các sản phẩm phái sinh phức tạp trên danh mục 100 tài sản, hướng tới mục tiêu cắt giảm 60% thời gian xử lý định giá rủi ro trong lộ trình triển khai 6 tháng.
  2. Xây dựng thư viện mã nguồn mở chuyên sâu về bài toán ngược PDEs: Các phòng thí nghiệm và nhóm nghiên cứu R&D công nghệ nên đóng gói thuật toán thành các module chuẩn hóa trên Python, hướng tới mục tiêu tối ưu độ trễ giải bài toán truyền nhiệt dừng xuống dưới 5 phút trong thời gian 3 tháng tới.
  3. Ứng dụng mô hình ResNet kết hợp sai phân hữu hạn trong mô phỏng vật liệu: Các viện nghiên cứu vật liệu và doanh nghiệp luyện kim cần áp dụng cấu trúc mạng nơ-ron kép để dự báo ranh giới chuyển pha trong bài toán Stefan, đảm bảo duy trì sai số biên di động dưới mức 0,1% trong khung thời gian thử nghiệm 9 tháng.
  4. Đổi mới chương trình đào tạo sau đại học ngành Toán - Tin ứng dụng: Các trường đại học kỹ thuật cần đưa chuyên đề Physics-Informed Machine Learning và kỹ thuật vi phân tự động vào giảng dạy, đặt mục tiêu 100% học viên cao học làm chủ công cụ giải bài toán ngược trong vòng 12 tháng.

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

Luận văn là nguồn tư liệu học thuật và thực tiễn giá trị cho 4 nhóm đối tượng chính:

  1. Chuyên gia phân tích định lượng (Quantitative Analysts) tại các định chế tài chính: Ứng dụng mô hình Black-Scholes 100 chiều để định giá chính xác các hợp đồng phái sinh đa tài sản có điều kiện rủi ro tín dụng và nguy cơ phá sản thực tế.
  2. Kỹ sư mô phỏng động lực học chất lưu và truyền nhiệt (CFD/Thermal Engineers): Tham khảo thuật toán ước lượng hệ số nhớt và hệ số dẫn nhiệt từ dữ liệu cảm biến rời rạc trong các hệ thống cơ khí và khí động học.
  3. Học viên cao học, nghiên cứu sinh ngành Toán ứng dụng, Khoa học dữ liệu và AI: Tiếp cận tài liệu chuẩn mực về lý thuyết BSDEs, liên hệ Feynman-Kac, cùng kỹ thuật kết hợp tối ưu hóa bậc 1 (Adam) và bậc 2 (L-BFGS).
  4. Kỹ sư phát triển phần mềm tính toán khoa học (Scientific Machine Learning Engineers): Khai thác cấu trúc mạng ResNet và phương pháp tích hợp vi phân tự động để phát triển các phần mềm mô phỏng kỹ thuật số thế hệ mới.

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

Mạng nơ-ron thông tin vật lý (PINNs) khác biệt gì so với mạng nơ-ron truyền thống?
Mạng nơ-ron truyền thống phụ thuộc hoàn toàn vào tập dữ liệu gắn nhãn lớn và dễ bị quá khớp khi thiếu dữ liệu. Ngược lại, PINNs tích hợp trực tiếp phương trình vi phân vào hàm tổn thất, đóng vai trò như một bộ điều chuẩn vật lý, giúp mô hình huấn luyện chính xác ngay cả khi chỉ có một lượng nhỏ điểm quan sát.

Mô hình trong luận văn giải quyết bài toán bùng nổ số chiều như thế nào?
Bằng cách chuyển đổi phương trình parabolic sang hệ phương trình vi phân ngẫu nhiên ngược (BSDEs) và xấp xỉ gradient thông qua mạng học sâu, mô hình đã giải quyết thành công bài toán Black-Scholes 100 chiều chỉ trong 36 phút, khắc phục hoàn toàn sự bùng nổ cấp số nhân của các phương pháp lưới cổ điển.

Tại sao cần kết hợp cả hai thuật toán tối ưu hóa Adam và L-BFGS?
Thuật toán bậc 1 Adam giúp mô hình thích ứng nhanh và vượt qua các điểm yên ngựa trong giai đoạn đầu. Sau đó, thuật toán bậc 2 L-BFGS sử dụng thông tin xấp xỉ ma trận Hessian để tăng tốc độ hội tụ, giúp sai số hàm tổn thất giảm sâu xuống mức 0,00006 mà không bị mắc kẹt tại cực trị địa phương.

Độ chính xác khi xác định hệ số phương trình Navier-Stokes đạt mức nào?
Trong thực nghiệm với 1000 điểm dữ liệu trên miền 2D, mô hình đã tìm ra chính xác hệ số nhớt thực nghiệm là 0,01 trùng khớp với giá trị lý thuyết, đưa hàm tổn thất toàn phần về mức 0,00006 và sai số chuẩn L2 của trường vận tốc gần như bằng 0.

Kiến trúc ResNet mang lại ưu thế gì so với mạng truyền thẳng thông thường?
Các đường nối tắt (skip connections) trong ResNet giúp dòng gradient lan truyền liên tục qua nhiều tầng ẩn, triệt tiêu hiện tượng triệt tiêu gradient (vanishing gradient). Thực nghiệm cho thấy ResNet giúp nghiệm hội tụ chuẩn xác về 1,450 thay vì bị dừng ở mức 1,686 như mạng truyền thẳng thông thường.

Kết luận

  • Luận văn đã xây dựng thành công khung phương pháp luận hoàn chỉnh sử dụng mạng học sâu và mạng PINNs để giải quyết bài toán xác định hệ số phương trình đạo hàm riêng.
  • Kiểm chứng thực nghiệm thành công trên 5 bài toán tiêu biểu: Poisson, Navier-Stokes, Black-Scholes 1 chiều và 100 chiều, cùng bài toán biên chuyển động Stefan.
  • Đạt bước tiến vượt bậc trong tính toán khoa học khi rút ngắn thời gian định giá rủi ro tài chính 100 chiều từ 140 phút xuống còn 36 phút với sai số tương đối chỉ 0,05%.
  • Đề xuất giải pháp kết hợp hiệu quả giữa mạng ResNet và hai thuật toán tối ưu Adam - L-BFGS giúp loại bỏ hiện tượng kẹt nghiệm cục bộ.
  • Cung cấp toàn bộ mã nguồn mở thực nghiệm, tạo tiền đề vững chắc cho việc ứng dụng AI vào mô phỏng vật lý và toán tài chính.

Trong giai đoạn 3 đến 6 tháng tiếp theo, hướng phát triển tự nhiên của đề tài là mở rộng giải thuật lên không gian trên 500 chiều và ứng dụng trên dữ liệu cảm biến công nghiệp thời gian thực. Hãy kết nối và tham khảo công trình nghiên cứu cùng mã nguồn thực nghiệm để ứng dụng ngay các giải pháp tính toán khoa học tiên tiến vào bài toán chuyên môn của bạn.