Tổng quan nghiên cứu

Sự cố tràn dầu trên biển luôn là hiểm họa môi trường nghiêm trọng, đe dọa trực tiếp đến hệ sinh thái hải dương và gây thiệt hại kinh tế hàng triệu USD mỗi năm. Luận văn Thạc sĩ Toán học chuyên ngành Toán học tính toán của tác giả Nguyễn Thu Quyên, thực hiện năm 2013 tại Trường Đại học Khoa học Tự nhiên thuộc Đại học Quốc gia Hà Nội dưới sự hướng dẫn của Giáo sư Tiến sĩ Đặng Quang Á, tập trung giải quyết bài toán cấp bách này bằng phương pháp giải tích hiện đại. Vấn đề cốt lõi của nghiên cứu là xây dựng cơ sở lý thuyết và thuật toán số để xử lý đồng thời hai bài toán: bài toán thuận nhằm mô phỏng quá trình khuếch tán, lan truyền nồng độ dầu theo không gian - thời gian, và bài toán ngược nhằm định vị chính xác tọa độ nguồn phát cùng thời điểm xảy ra sự cố khi phát hiện chùm vết dầu trên biển.

Phạm vi nghiên cứu được triển khai thực nghiệm trên mô hình không gian một chiều và mở rộng sang không gian hai chiều tại vùng biển Đông thuộc khu vực Nam Trung Bộ Việt Nam, giới hạn từ kinh tuyến 108 độ Đông đến 112 độ Đông và vĩ tuyến 10 độ Bắc đến 14 độ Bắc. Với kịch bản giả định lượng dầu tràn lên tới 10.000 kg và theo dõi trong khung thời gian 60 giờ, luận văn đã chứng minh tính hữu hiệu vượt trội của phương pháp phương trình liên hợp. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc cắt giảm hơn 80% khối lượng tính toán so với các phương pháp quét lưới truyền thống, đồng thời cung cấp công cụ toán học tin cậy hỗ trợ các nhà quản lý môi trường ra quyết định kịp thời trong công tác cứu hộ bờ biển và bảo tồn các vùng sinh thái nhạy cảm.

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 phương trình đạo hàm riêng và giải tích hàm hiện đại. Trọng tâm lý thuyết là phương pháp phương trình liên hợp dựa trên đẳng thức Lagrange mở rộng, thiết lập mối quan hệ đối ngẫu giữa toán tử vi phân của bài toán ban đầu và toán tử liên hợp. Mô hình thủy động lực học lan truyền vật chất ô nhiễm được mô tả bởi phương trình truyền tải - khuếch tán - phân rã dạng Parabolic phi dừng, có tính đến trường vận tốc dòng chảy, hệ số khuếch tán dòng xoáy và tốc độ bốc hơi tự nhiên của dầu.

Khung lý thuyết của luận văn tích hợp 4 khái niệm cốt lõi: phiếm hàm độ nhạy xác định mối tương quan giữa nguồn phát và vùng quan trắc, không gian Hilbert các hàm suy rộng, định lý biểu diễn Riesz hỗ trợ chứng minh tính tồn tại của toán tử nghịch đảo bị chặn, và chuẩn năng lượng trong không gian Banach. Về mặt toán học thuần túy, tác giả đã sử dụng kỹ thuật ước lượng tiên nghiệm năng lượng để chứng minh chặt chẽ định lý về sự tồn tại, tính duy nhất và tính ổn định nghiệm suy rộng của bài toán tràn dầu hai chiều với các điều kiện biên hỗn hợp. Kết quả phân tích chứng minh rằng khi không còn nguồn xả mới, tổng lượng dầu tích phân trên toàn miền và chuẩn bậc hai của nồng độ dầu luôn suy giảm đơn điệu theo thời gian.

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

Nguồn dữ liệu tính toán của luận văn dựa trên các tham số vật lý và khí tượng thủy văn thực nghiệm, bao gồm trường vận tốc dòng chảy hai chiều thỏa mãn phương trình liên tục, hệ số khuếch tán dòng xoáy từ 0,9 mét vuông trên giây đến 3,0 mét vuông trên giây, cùng các giá trị hệ số phân rã bốc hơi từ 0 đến 0,1 trên giây.

Phương pháp phân tích số được lựa chọn là phương pháp sai phân hữu hạn kết hợp lược đồ luân hướng từng phần, hay còn gọi là phương pháp tách hướng không gian. Cỡ mẫu lưới sai phân được thiết lập cụ thể: không gian một chiều sử dụng lưới 301 điểm nút trên đoạn tọa độ từ -50 đến 50 với bước thời gian 0,5 giây; không gian hai chiều sử dụng mạng lưới không gian đồng nhất 201 nhân 201 điểm nút, tương đương 40.401 điểm tính toán, với bước thời gian tích phân là 300 giây. Lý do lựa chọn lược đồ sai phân luân hướng kết hợp toán tử làm trơn đơn điệu là nhằm đảm bảo lược đồ luôn ổn định tuyệt đối, triệt tiêu hiện tượng dao động số phi vật lý và đạt sai số chặt cụt cấp hai theo không gian và cấp một theo thời gian. Toàn bộ thuật toán giải số được lập trình và kiểm nghiệm qua các mốc thời gian mô phỏng từ 0 đến 60 giờ.

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

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

Nghiên cứu đã mang lại 4 phát hiện khoa học mang tính đột phá trong việc ứng dụng phương trình liên hợp để giải bài toán môi trường:

Thứ nhất, trong bài toán thuận một chiều, mô hình mô phỏng chính xác sự biến thiên của trường nồng độ dầu dưới ảnh hưởng của hệ số bốc hơi. Tại thời điểm 30 giây kể từ khi xuất hiện nguồn điểm, nồng độ cực đại đạt mức 0,052 đơn vị quy chuẩn tại tọa độ 26,3 mét khi hệ số bốc hơi bằng 0, và giảm mạnh hơn 35% khi đưa hệ số phân rã lên 0,1 trên giây.

Thứ hai, trong bài toán ngược một chiều, phương pháp liên hợp đã truy vết chính xác tuyệt đối vị trí nguồn phát ban đầu tại tọa độ gốc 0 mét. Khi thiết lập 5 điểm quan trắc nồng độ phân tán trên miền lan truyền tại các vị trí 20 mét, 24 mét, 26,3 mét, 30 mét và 33 mét, các đường đồng mức nồng độ liên hợp ngược thời gian đều giao nhau tại duy nhất một điểm gốc ban đầu với độ lệch sai số bằng 0.

Thứ ba, trong bài toán không gian hai chiều mô phỏng vùng biển Nam Trung Bộ, thuật toán đã tái hiện thành công quá trình phát tán của 10.000 kg dầu sau 60 giờ trôi dạt. Vùng nồng độ mở rộng theo hướng trường gió Tây Nam với diện tích bao phủ hàng trăm kilômét vuông mặt biển.

Thứ tư, khi giải bài toán ngược hai chiều từ dữ liệu chùm dầu thu nhận được tại mốc 60 giờ mà không cần biết chi tiết nồng độ bên trong, phương pháp liên hợp đã thu hẹp liên tiếp các vùng chứa nguồn phát tiềm năng qua các mốc thời gian ngược 12 giờ, 24 giờ, 42 giờ và định vị chính xác vị trí tàu gặp sự cố tại tọa độ 111 độ kinh Đông, 13 độ vĩ Bắc.

Thảo luận kết quả

Các kết quả đạt được khẳng định tính ưu việt vượt trội của việc kết hợp phương trình vi phân liên hợp với thuật toán sai phân luân hướng. Thay vì phải giải hàng nghìn bài toán thuận tương ứng với mọi điểm lưới khả dĩ trên toàn bộ diện tích biển 400 kilômét nhân 440 kilômét, kỹ sư tính toán chỉ cần giải từ 1 đến 5 bài toán liên hợp duy nhất theo hướng ngược thời gian. Điều này giúp giảm thiểu khoảng 85% thời gian tính toán của bộ vi xử lý và tiết kiệm đáng kể tài nguyên bộ nhớ.

Trong thực tế ứng dụng, các dữ liệu nồng độ và phân bố không gian được trình bày trực quan thông qua biểu đồ đường đẳng trị 2D và biểu đồ mặt cong nồng độ 3D. Đồ thị cho thấy các đường mức nồng độ co cụm dần thành một điểm cực trị đơn nhất khi lùi thời gian về thời điểm bắt đầu tràn dầu. So sánh với các nghiên cứu mô phỏng bằng phương pháp phần tử hữu hạn hay phương pháp hạt Monte Carlo truyền thống, phương pháp phương trình liên hợp của luận văn thể hiện độ ổn định cao hơn, không bị phụ thuộc vào số lượng hạt ngẫu nhiên và giải quyết triệt để bài toán thiếu hụt dữ liệu đầu vào ban đầu.

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

Dựa trên các kết quả lý thuyết và mô phỏng số đã đạt được, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể:

Thứ nhất, tích hợp thuật toán phương trình liên hợp vào hệ thống giám sát sự cố tràn dầu quốc gia của các cơ quan quản lý biển. Chủ thể thực hiện là Cục Biển và Hải đảo phối hợp với Trung tâm Ứng phó sự cố tràn dầu khu vực, hướng tới mục tiêu tự động hóa việc xác định tọa độ tàu xả thải lén lút với độ chính xác đạt trên 95% trong vòng 24 giờ kể từ khi nhận diện vệt dầu từ ảnh vệ tinh viễn thám.

Thứ hai, chuẩn hóa và số hóa cơ sở dữ liệu hải dương học thời gian thực. Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần triển khai hệ thống đo đạc dòng chảy tầng mặt và gió biển với mạng lưới phân giải cao 100 nhân 100 điểm nút trên toàn vùng biển kinh tế trọng điểm trong lộ trình 18 tháng tới, cung cấp trường vận tốc chuẩn xác làm đầu vào cho mô hình toán.

Thứ ba, phát triển gói phần mềm chuyên dụng hỗ trợ ra quyết định khẩn cấp. Các viện nghiên cứu toán học và công nghệ thông tin cần chuyển đổi mã nguồn tính toán sang các ngôn ngữ hiệu năng cao như C++ hoặc Python song song hóa, đặt chỉ tiêu xử lý bài toán ngược hoàn tất dưới 15 phút cho các sự cố có quy mô từ 5.000 đến 20.000 kg dầu, thực hiện trong vòng 12 tháng.

Thứ tư, tăng cường đào tạo liên ngành Toán ứng dụng và Khoa học môi trường biển. Các trường đại học khối khoa học tự nhiên cần tổ chức ít nhất 4 khóa tập huấn chuyên sâu hàng năm cho kỹ sư và nghiên cứu viên nhằm nâng cao năng lực ứng dụng toán học tính toán trong việc giải quyết các thách thức sinh thái biển.

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

Luận văn là tài liệu tham khảo giá trị cao cho 4 nhóm đối tượng trọng tâm:

Nhóm thứ nhất là các nhà nghiên cứu chuyên sâu về Toán ứng dụng và Toán học tính toán. Tài liệu cung cấp minh chứng chặt chẽ về việc thiết lập bài toán liên hợp, phân tích giải tích hàm trong không gian Hilbert và phương pháp giải số cho hệ phương trình Parabolic đa chiều.

Nhóm thứ hai là các kỹ sư môi trường, chuyên gia hải dương học và cán bộ kỹ thuật tại các trung tâm cứu hộ hàng hải. Nghiên cứu cung cấp công cụ mô phỏng trực quan giúp lập phương án quây chặn vết dầu, thả phao gom và sử dụng hóa chất phân tán đúng vị trí trọng yếu trên biển.

Nhóm thứ ba là các cơ quan quản lý nhà nước về tài nguyên, môi trường biển và cảnh sát biển. Luận văn cung cấp luận cứ khoa học có giá trị pháp lý vững chắc để truy vết, xác minh trách nhiệm pháp lý và yêu cầu bồi thường thiệt hại đối với các phương tiện vận tải gây ô nhiễm môi trường nước.

Nhóm thứ tư là học viên cao học, nghiên cứu sinh và sinh viên ngành Toán - Tin học, Vật lý địa cầu và Kỹ thuật tài nguyên nước. Đây là tài liệu tham khảo mẫu mực về cấu trúc nghiên cứu, phương pháp xây dựng mô hình toán học từ thực tiễn và kỹ thuật lập trình sai phân luân hướng.

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

Phương pháp phương trình liên hợp có ưu thế vượt trội gì so với phương pháp mô phỏng thuận thông thường? Thay vì phải thử hàng nghìn vị trí nguồn phát giả định rồi chạy mô phỏng xuôi thời gian, phương pháp liên hợp chỉ cần giải một bài toán ngược duy nhất từ dữ liệu quan trắc hiện tại. Điều này giúp giảm hơn 80% thời gian xử lý dữ liệu và tiết kiệm tài nguyên máy tính tối đa.

Mô hình có thể xác định được nguồn tràn dầu nếu chỉ biết nồng độ dầu tại một vài điểm rải rác không? Hoàn toàn có thể. Khi thu thập được nồng độ tại ít nhất 3 đến 5 vị trí quan trắc trong chùm dầu, phương pháp sẽ giải số các phương trình liên hợp tương ứng. Các đường đẳng trị nồng độ ngược thời gian sẽ giao nhau tại một giao điểm duy nhất, đó chính là tọa độ nguồn phát ban đầu.

Lược đồ sai phân luân hướng trong luận văn đạt bậc chính xác nào và có ổn định không? Lược đồ sai phân luân hướng được thiết lập với toán tử đơn điệu, đạt sai số chặt cụt cấp hai theo không gian và cấp một theo thời gian. Thuật toán đã được chứng minh ổn định tuyệt đối và không xuất hiện sai số phân tán số học khi tính toán trên lưới 201 nhân 201 điểm nút.

Hệ số bốc hơi của dầu ảnh hưởng như thế nào đến kết quả truy tìm vị trí nguồn phát? Hệ số bốc hơi làm giảm tổng lượng dầu trên mặt nước theo hàm mũ thời gian. Trong mô hình số một chiều, khi hệ số bốc hơi tăng từ 0 lên 0,1 trên giây, nồng độ cực đại giảm đi khoảng 35%, nhưng vị trí tọa độ giao điểm của các đường mức liên hợp vẫn hội tụ chính xác về nguồn gốc ban đầu.

Mô hình toán học trong luận văn có thể mở rộng cho các dạng ô nhiễm môi trường khác không? Mô hình có tính tổng quát hóa cao và có thể áp dụng ngay cho các bài toán lan truyền bụi mịn trong không khí, chất phóng xạ hoặc hóa chất độc hại trong sông ngòi và đầm phá với hệ số khuếch tán điều chỉnh từ 0,1 đến 5 mét vuông trên giây.

Kết luận

Luận văn Thạc sĩ của tác giả Nguyễn Thu Quyên đã giải quyết trọn vẹn và xuất sắc mục tiêu ứng dụng phương trình liên hợp vào bài toán môi trường thực tiễn với 5 đóng góp nổi bật:

  • Xây dựng thành công khung lý thuyết hoàn chỉnh cho bài toán thuận và bài toán ngược truyền tải khuếch tán chất ô nhiễm trong không gian một chiều và hai chiều.
  • Chứng minh chặt chẽ định lý về sự tồn tại, tính duy nhất và ổn định nghiệm suy rộng của bài toán tràn dầu hai chiều trong không gian Hilbert.
  • Thiết lập lược đồ sai phân luân hướng đơn điệu ổn định tuyệt đối với sai số bậc hai theo không gian và bậc một theo thời gian.
  • Định vị thành công nguồn phát 10.000 kg dầu tại tọa độ 111 độ kinh Đông, 13 độ vĩ Bắc trên vùng biển Đông sau 60 giờ phát tán mô phỏng.
  • Đề xuất quy trình tính toán tối ưu giúp giảm hơn 80% thời gian xử lý so với các phương pháp giải số truyền thống.

Trong giai đoạn tiếp theo, mô hình cần được phát triển tích hợp với dữ liệu gió biển theo thời gian thực và tự động hóa mã nguồn trên các nền tảng tính toán đám mây. Hãy tham khảo toàn văn luận văn để tiếp cận đầy đủ các chứng minh toán học chuẩn mực và thuật giải số chi tiết phục vụ công tác nghiên cứu và ứng dụng thực tế.