Tổng quan về luận án

Sự phát triển mạnh mẽ của kinh tế biển và công nghiệp dầu khí tại khu vực Biển Đông đã làm gia tăng đáng kể rủi ro ô nhiễm môi trường biển. Nằm trên tuyến hàng hải huyết mạch có lưu lượng vận tải dầu lớn thứ hai thế giới với khoảng 300 - 400 triệu tấn dầu thô lưu thông mỗi năm, Biển Đông đối mặt với nguy cơ thảm họa tràn dầu thường trực. Theo thống kê của Liên đoàn Ứng phó Ô nhiễm Biển Quốc tế (ITOPF), trong giai đoạn 1992 - 2011, Việt Nam đứng thứ 5 trên thế giới, thứ 4 tại Châu Á và thứ 3 trong khu vực Đông Nam Á về tần suất xảy ra các sự cố dầu tràn từ 7 tấn trở lên. Thực tế lịch sử ghi nhận hơn 50 sự cố dầu tràn nghiêm trọng từ năm 1990 đến năm 2015, điển hình là sự cố tràn dầu "bí ẩn" không rõ nguồn gốc trong các năm 2007 - 2008 làm dạt vào bờ biển Việt Nam hơn 5.000 tấn dầu trôi nổi, gây tổn thất nặng nề cho hệ sinh thái ven bờ, ngành du lịch và nuôi trồng thủy hải sản.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng các mô hình mô phỏng dầu tràn hiện nay trên thế giới và tại Việt Nam chủ yếu tiếp cận theo hệ quy chiếu Lagrange (Lagrangian Particle-Tracking Models - LPTMs) để dự báo xuôi thời gian. Việc giải bài toán ngược thời gian (Backward-in-Time Trajectory - BITT) để truy tìm nguồn gốc phát thải vẫn chủ yếu dựa trên kỹ thuật hạt đảo ngược đơn giản hoặc thử sai trực quan, thường bỏ qua tính tương thích về trường vật lý và các quá trình phong hóa phi tuyến tính. Luận án tiến sĩ chuyên ngành Hải dương học (mã số 62 44 02 28) của nghiên cứu sinh Nguyễn Quốc Trinh với đề tài "Nghiên cứu lan truyền dầu ở Biển Đông phục vụ cảnh báo và tìm kiếm nguồn thải" đã thiết lập hệ thống phương pháp tiếp cận Euler (Eulerian approach) toàn diện cho cả hai bài toán: cảnh báo lan truyền xuôi thời gian ($t > 0$) và truy nguyên nguồn phát thải ngược thời gian ($t < 0$).

Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập mạch lạc:

  • RQ1: Làm thế nào để mô hình hóa đồng bộ quá trình động lực học (bình lưu, khuếch tán) và 8 quá trình phong hóa hóa-lý của vệt dầu tràn trên mặt biển theo tiếp cận Euler?
  • RQ2: Cơ chế toán học và điều kiện biên nào cho phép nghịch đảo thời gian phương trình bình lưu - khuếch tán Euler nhằm xác định vị trí và thời điểm phát thải ban đầu với độ tin cậy cao?
  • H1: Mô hình Euler tích hợp sơ đồ sai phân dưới lưới (sub-grid) có khả năng triệt tiêu hiệu ứng khuếch tán số, nâng cao độ chính xác mô phỏng vệt dầu tương đương hoặc vượt trội so với các mô hình Lagrange truyền thống.
  • H2: Bài toán ngược Euler kết hợp nghịch đảo vận tốc trường động lực và bù trừ phong hóa cho phép tái cấu trúc không gian nguồn thải dầu trôi nổi trên Biển Đông trong phạm vi sai số chấp nhận được ($R > 0,85$; $NASH > 0,70$).

Nghiên cứu giới hạn phạm vi không gian toàn Biển Đông từ $1^\circ\text{N}$ đến $24^\circ\text{N}$ và $99^\circ\text{E}$ đến $121^\circ\text{E}$, sử dụng chuỗi dữ liệu tái phân tích khí quyển CFSR (NOAA), mô hình thủy động lực POM (Princeton Ocean Model) và sóng SWAN trong giai đoạn biến động thực địa 2007 - 2008.

Literature Review và Positioning

Các nghiên cứu mô phỏng sự cố tràn dầu trên thế giới khởi đầu từ các mô hình biến đổi phong hóa tĩnh của Blokker (1964) và Fay (1971), sau đó mở rộng thành các hệ thống tính toán động học đa chiều. Trong dòng nghiên cứu tiếp cận Lagrange, các mô hình thương mại và nghiên cứu nổi tiếng như GNOME (NOAA, Hoa Kỳ), MEDSLIK (Lardner, 1998), MIKE 31 SA/PA (DHI, Đan Mạch), và Delft3D-PART (Deltares, Hà Lan) sử dụng kỹ thuật ngẫu hành (random-walk) để theo dõi các hạt chất điểm. Nhóm mô hình này có ưu điểm tính toán nhanh nhưng tồn tại nhược điểm lớn là không sử dụng trực tiếp các trường môi trường liên tục mà phải thực hiện nội suy, ngoại suy không gian, dễ làm mất mát thông tin động lực quy mô nhỏ.

Ngược lại, dòng nghiên cứu tiếp cận Euler tập trung giải phương trình bảo toàn khối lượng và phương trình bình lưu - khuếch tán liên tục của chất ô nhiễm, tiêu biểu như mô hình OSCAR (Reed nnk., 1995), MOHID (Martins nnk., 2001), và OILMAP (Spaulding nnk., 1993). Ưu thế tuyệt đối của mô hình Euler là tính đồng nhất hoàn toàn về cấu trúc lưới và phương pháp giải số với các mô hình hoàn lưu đại dương (như POM, ROMS) và mô hình khí tượng (WRF).

Tại Việt Nam, các nghiên cứu tiên phong của Phạm Văn Ninh và Đinh Văn Mạnh (Viện Cơ học) với phần mềm OST-2D/3D, Vũ Thanh Ca nnk. (Viện Khoa học Khí tượng Thủy văn và Môi trường), Nguyễn Hữu Nhân nnk. với OilSAS, và Đinh Văn Ưu nnk. (Đại học Khoa học Tự nhiên) đã từng bước tiếp cận cả hai phương pháp. Tuy nhiên, về bài toán ngược thời gian truy tìm nguồn gốc phát thải, y văn quốc tế và trong nước còn tồn tại nhiều tranh luận học thuật sâu sắc:

Cuộc tranh luận thứ nhất diễn ra giữa quan điểm sử dụng cấu trúc nhánh kết hợp Lagrangian Coherent Structures (LCSs) của Peacock nnk. (MIT) đối lập với phương pháp quỹ đạo hạt ngược BITT (Backward-in-Time Trajectory) của Batchelder (2006) và Isobe (2009). Nhóm LCSs nhấn mạnh các ranh giới vận chuyển tàn dư (hyperbolic transport barriers), trong khi nhóm BITT lại vấp phải hạn chế nghiêm trọng về độ hội tụ (nghiên cứu thực nghiệm tại vùng OMEX Na Uy của Christensen nnk. chỉ đạt độ chính xác nguồn $0,1%$).

Cuộc tranh luận thứ hai nằm ở tính khả thi toán học của phương trình nghịch đảo thời gian Euler. Nhiều nhà nghiên cứu cho rằng bài toán ngược khuếch tán mang tính "ill-posed" (không chỉnh) do tính tiêu tán entropy không thuận nghịch. Nghiên cứu của Nguyễn Quốc Trinh đã định vị chính xác vị trí học thuật bằng cách đưa ra giải pháp giải số sai phân ổn định kết hợp điều kiện ràng buộc vật lý không âm ($C \ge 0$), đối chiếu trực tiếp với các thảm họa quốc tế như Ixtoc 1 (1979, Vịnh Mexico) và Deepwater Horizon (2010, BP).

Đó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à làm phong phú hệ thống lý thuyết thủy động lực học môi trường biển và lý thuyết lan truyền chất ô nhiễm thông qua các đóng góp đột phá:

  • Mở rộng định luật khuếch tán Fick và cơ học môi trường liên tục: Thiết lập công thức toán học mô tả sự nghịch đảo có điều kiện của quá trình tiêu tán vật chất. Trong bài toán xuôi, dầu tự loang và tiêu tán; trong bài toán ngược, quá trình thu gom và tích tụ vật chất về điểm kỳ dị (nguồn phát) được thiết lập thông qua toán tử bảo toàn khối lượng điều chỉnh.
  • Tích hợp tương tác phi tuyến dòng chảy - sóng - gió: Kế thừa lý thuyết dòng trôi Ekman (1905) và dòng sóng Stokes (1847), luận án phát triển công thức vận tốc trôi dạt tổng hợp $\mathbf{V}_a = \gamma_w D_w \mathbf{V}_w + \gamma_c \mathbf{V}_c + \gamma_s \mathbf{V}_s$, trong đó góc lệch $\theta$ và các hệ số trôi dạt $\gamma_w, \gamma_c, \gamma_s$ được tham số hóa phụ thuộc phi tuyến vào tốc độ gió mực 10m và cấu trúc trường dòng chảy thực tế thay vì sử dụng hệ số cố định như Stolzenbach nnk. (1977).
  • Hệ tiên đề bảo toàn và liên hệ pha: Xây dựng mối quan hệ hàm số xác định giữa ba đại lượng trạng thái bề mặt: nồng độ diện tích $C$ ($\text{kg}\cdot\text{m}^{-2}$), độ dày lớp dầu $h$ ($\text{m}$), và mật độ dầu $\rho$ ($\text{kg}\cdot\text{m}^{-3}$) qua biểu thức $C = \rho \cdot h$, giải quyết triệt để sự thiếu nhất quán trong việc mô phỏng đồng thời tính chất hóa học và động học của vệt dầu.

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

Khung phân tích của luận án tích hợp đồng bộ ba lý thuyết nền tảng: (1) Lý thuyết động lực học chất lỏng phi nén Navier-Stokes, (2) Lý thuyết phong hóa động lực học dầu khí của Mackay (1980), và (3) Lý thuyết khuếch tán rối đại dương Okubo (1971).

Hệ phương trình chủ đạo của bài toán xuôi thời gian:

$$\frac{\partial C}{\partial t} + \frac{\partial (V_x C)}{\partial x} + \frac{\partial (V_y C)}{\partial y} = \frac{\partial}{\partial x}\left(D_x \frac{\partial C}{\partial x}\right) + \frac{\partial}{\partial y}\left(D_y \frac{\partial C}{\partial y}\right) \pm C_m$$

Hệ phương trình bài toán ngược thời gian tái cấu trúc nguồn:

$$\frac{\partial C}{\partial t} - \frac{\partial (V_x C)}{\partial x} - \frac{\partial (V_y C)}{\partial y} = -\frac{\partial}{\partial x}\left(D_x \frac{\partial C}{\partial x}\right) - \frac{\partial}{\partial y}\left(D_y \frac{\partial C}{\partial y}\right) + C_m$$

Hệ thống phong hóa 8 thành phần được tham số hóa vi phân toàn phần, chuẩn hóa thể tích lớp dầu $V = V_0(1 - F_e - F_d - F_{disc} - F_{sed} - F_{coast} - F_{oxy} - F_{bio})(1 + F_w)$ và biến thiên độ nhớt $\mu$, sức căng bề mặt $\sigma_{w}$, $\sigma_{A}$ theo nhiệt độ và thời gian phơi nhiễm.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái nhận thức luận thực chứng (positivism) với phương pháp luận định lượng nghiêm ngặt kết hợp mô hình hóa số trị quy mô lớn. Cấu trúc mô hình được thiết kế phân tầng không gian đa cấp, kết hợp giữa lưới tính toán nền (coarse grid) và cấu trúc lồng ghép thích nghi dưới lưới (sub-grid refinement) theo các tỷ lệ $2^{-4}$ và $2^{-8}$.

Quy trình giải bài toán số trị sử dụng phương pháp sai phân hữu hạn (Finite Difference Method) trên lưới xen kẽ Arakawa-C với các bước rời rạc không gian được khảo sát chi tiết từ $0,01^\circ; 0,02^\circ; 0,03^\circ; 0,04^\circ; 0,044^\circ; 0,05^\circ$ đến $0,06^\circ$ (tương đương quy mô từ ~1,1 km đến ~6,6 km).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu kiểm định qua 5 bước nghiêm ngặt:

  1. Thu thập và tiền xử lý dữ liệu trường: Tích hợp dữ liệu khí tượng tái phân tích CFSR độ phân giải $0,3^\circ \times 0,3^\circ$ và $0,5^\circ \times 0,5^\circ$; trường hoàn lưu dòng chảy mặt từ mô hình 3D POM ($0,044^\circ$); trường sóng mặt từ SWAN.
  2. Kiểm tra tính hội tụ và triệt tiêu khuếch tán số: Sử dụng kỹ thuật sai phân ngược dòng (upwind scheme) bậc 1 kết hợp hiệu chỉnh bậc 2 Lax-Wendroff trên cấu trúc dưới lưới $2^{-8}$ nhằm cô lập và triệt tiêu sai số khuếch tán số học vốn là điểm yếu cố hữu của mô hình Euler.
  3. Triangulation dữ liệu: Đối soát chéo giữa kết quả mô phỏng số trị, chuỗi ảnh viễn thám quang học/SAR đa thời gian (ENVISAT, MODIS) năm 2007 - 2008 của Cục Viễn thám Quốc gia, và dữ liệu thực tế hơn 5.000 tấn dầu vón cục thu gom tại bờ biển miền Trung Việt Nam.
  4. Kiểm chuẩn độ ổn định thuận nghịch (Forward-Backward Loop): Thực hiện chu trình kép: mô phỏng xuôi từ nguồn giả định trong 7 ngày $\rightarrow$ lấy kết quả làm đầu vào mô phỏng ngược 7 ngày $\rightarrow$ đánh giá độ lệch tâm và biến dạng diện tích vệt dầu.

Data và phân tích

Đặc tính mẫu và thông số tính toán: Dữ liệu dầu thô chuẩn hóa (Bạch Hổ, DO, FO, Bunker C) với khối lượng riêng ban đầu $\rho_0 \in [800, 990]\text{ kg}\cdot\text{m}^{-3}$, độ nhớt $\mu_0 \in [5,0; 50,0]\text{ cP}$.

Các chỉ số thống kê định lượng đánh giá hiệu quả mô hình gồm:

  • Hệ số Nash-Sutcliffe Efficiency (NASH): Đánh giá độ khớp hình dạng và phân bố nồng độ.
  • Hệ số Percent Bias (PBIAS): Đo lường xu hướng ước tính thừa hoặc thiếu của thể tích dầu.
  • Sai số toàn phương trung bình (RMSE): Định lượng sai số vị trí tọa độ vệt dầu.
  • Hệ số tương quan Pearson ($R$): Đo lường tính đồng biến quỹ đạo di chuyển.

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

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

  • Triệt tiêu sai số khuếch tán số bằng kỹ thuật dưới lưới: Kết quả mô phỏng đối chứng chỉ ra rằng việc áp dụng lưới nền thông thường tạo ra diện tích loang dầu giả tạo lớn hơn $35% - 42%$ so với thực tế do khuếch tán số. Khi kết hợp cấu trúc dưới lưới $2^{-8}$, hình thái vết dầu co cụm chính xác, hệ số tương quan đạt $R = 0,94$, chỉ số NASH tăng từ $0,52$ lên $0,88$.
  • Xác định cơ chế trôi dạt đặc thù trên Biển Đông: Sự lan truyền vệt dầu trong điều kiện gió mùa Đông Bắc (tháng 2/2007) bị chi phối áp đảo bởi dòng chảy ven bờ chảy xiết về phía Nam kết hợp với ứng suất gió mặt ($V_w > 10\text{ m/s}$), làm vệt dầu dịch chuyển với tốc độ $0,45 - 0,65\text{ m/s}$. Ngược lại, trong gió mùa Tây Nam (tháng 7/2008), hoàn lưu xoáy nước ngoài khơi miền Trung làm phân tán vệt dầu ra vùng biển khơi với cơ chế khuếch tán chiếm ưu thế.
  • Nghịch đảo thời gian thành công xác định nguồn thải bí ẩn 2007 - 2008: Mô phỏng ngược thời gian trước 5 - 7 ngày từ các vệt dầu ghi nhận tại bờ biển Trung Bộ đã hội tụ thành công về các vùng tọa độ hẹp nằm trên tuyến hàng hải quốc tế trọng điểm ngoài khơi quần đảo Hoàng Sa và vùng biển ngoài khơi vịnh Bắc Bộ. Tọa độ nguồn phát thải ngược có độ lệch tâm cực đại $< 12,4\text{ km}$ sau 7 ngày mô phỏng ngược.
  • Tính toán định lượng quá trình phong hóa: Sau 120 giờ lan truyền, tỷ phần bay hơi $F_e$ đạt giá trị bão hòa $38% - 44%$ đối với dầu thô nhẹ, hàm lượng nước trong nhũ tương $F_w$ đạt ngưỡng cực đại $70% - 75%$, làm tăng thể tích khối chất lỏng lên gấp 3 lần và độ nhớt động lực học tăng đột biến hơn 2 bậc độ lớn ($> 1.200\text{ cP}$).
Kịch Bản Mô Phỏng Bước Lưới ($\Delta x, \Delta y$) Thời Gian Tính Hệ Số NASH PBIAS (%) RMSE (km) Hệ Số Tương Quan ($R$)
Đồng nhất (Xuôi) $0,044^\circ$ 120 giờ 0,92 -3,4% 3,2 0,96
Chế độ tháng (Xuôi) $0,044^\circ$ 72 giờ 0,84 +5,8% 8,7 0,89
Thực địa 05/02/2007 (Xuôi) $0,044^\circ$ (Dưới lưới $2^{-8}$) 168 giờ 0,86 +4,2% 7,4 0,91
Thực địa 05/07/2008 (Xuôi) $0,044^\circ$ (Dưới lưới $2^{-8}$) 168 giờ 0,83 -6,1% 9,1 0,88
Nghịch đảo 25/02/2007 (Ngược) $0,044^\circ$ (Dưới lưới $2^{-8}$) -168 giờ 0,79 +8,5% 12,4 0,85
Nghịch đảo 25/07/2008 (Ngược) $0,044^\circ$ (Dưới lưới $2^{-8}$) -168 giờ 0,76 -9,2% 14,1 0,82

Implications đa chiều

  • Đóng góp lý thuyết: Cung cấp bằng chứng thực nghiệm khẳng định tính ứng dụng thành công của mô hình Euler trong bài toán ngược thời gian dòng chảy địa vật lý, phá vỡ định kiến cho rằng chỉ có mô hình hạt Lagrange mới giải được bài toán truy tìm nguồn gốc.
  • Đổi mới phương pháp luận: Khung mô hình sai phân dưới lưới đa tỷ lệ mở ra hướng tiếp cận mới cho các bài toán phân tán chất ô nhiễm phóng xạ, vi nhựa và tìm kiếm cứu nạn hàng hải (SAR) trên vùng biển nhiệt đới.
  • Ứng dụng thực tiễn và chính sách: Cung cấp công cụ hỗ trợ ra quyết định (Decision Support System) trực tiếp cho Trung tâm Ứng phó Sự cố Tràn dầu Khu vực Miền Trung và Cục Hàng hải Việt Nam, chuyển đổi từ thế bị động thu gom rác ven bờ sang chủ động phát hiện, chặn bắt và xử lý pháp lý đối với các tàu xả thải trộm trên vùng đặc quyền kinh tế (EEZ).

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 hạn chế khoa học nội tại:

  1. Hiện tượng khuếch tán số tàn dư trong bài toán ngược dài hạn: Dù đã áp dụng kỹ thuật dưới lưới $2^{-8}$, khi thời gian mô phỏng ngược vượt quá 10 ngày ($-t > 240\text{ giờ}$), sai số khuếch tán tích lũy làm vùng nguồn ước tính bị mở rộng hình nan quạt, làm giảm độ sắc nét của tâm nguồn phát thải.
  2. Giả định 2D cho lớp dầu bề mặt: Mô hình hiện tại tập trung vào lớp dầu màng nổi 2 chiều, các tham số trao đổi phân tán thẳng đứng ($F_d$) và hòa tan ($F_{disc}$) chỉ được tham số hóa gián tiếp dạng thông lượng chìm mà chưa giải trực tiếp phương trình 3D hoàn chỉnh trong cột nước đại dương.
  3. Độ phân giải dữ liệu trường khí tượng ven bờ: Dữ liệu gió tái phân tích CFSR ($0,3^\circ - 0,5^\circ$) tuy thể hiện xuất sắc quy mô biển khơi nhưng chưa phản ánh đầy đủ hiệu ứng địa hình cục bộ phức tạp tại các vũng, vịnh kín ven bờ miền Trung.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng chiến lược:

  • Tích hợp mô hình 3D Euler - Lagrange lai (Hybrid 3D Modeling): Kết hợp tính toán trường nồng độ tổng thể bằng Euler và theo dõi các cụm váng dầu siêu nhỏ bằng hạt Lagrange trong cột nước sâu.
  • Đồng hóa dữ liệu thời gian thực (Real-time Data Assimilation): Ứng dụng bộ lọc Kalman mở rộng (EnKF) để đồng hóa liên tục dữ liệu ảnh viễn thám radar khẩu độ tổng hợp (SAR) và thiết bị trôi dạt đại dương (drifter) vào mô hình số.
  • Ứng dụng Trí tuệ Nhân tạo (AI/Machine Learning): Sử dụng mạng nơ-ron tích chập vật lý (PINNs) để giải bài toán nghịch đảo tham số nguồn phát thải trong thời gian thực dưới 1 phút.

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

Công trình nghiên cứu tạo ra những tác động học thuật và thực tiễn sâu rộng:

  • Tác động học thuật: Đặt nền móng cho trường phái mô hình hóa nghịch đảo Euler tại Việt Nam trong hải dương học ứng dụng, mở ra tiềm năng trích dẫn cao trong các tạp chí chuyên ngành quốc tế thuộc danh mục ISI/Scopus (như Marine Pollution Bulletin, Ocean Engineering, Journal of Marine Systems).
  • Chuyển giao ngành công nghiệp dầu khí: Cung cấp cơ sở dữ liệu tính toán nguy cơ tràn dầu cho Tập đoàn Dầu khí Quốc gia Việt Nam (PVN) và các liên doanh khai thác (Vietsovpetro) trong việc xây dựng báo cáo Đánh giá Tác động Môi trường (ĐTM) và kế hoạch ứng phó sự cố cho các mỏ Bạch Hổ, Đại Hùng, Rồng, Cá Ngừ Vàng.
  • Hoạch định chính sách quốc gia: Phục vụ trực tiếp cho Ủy ban Quốc gia Ứng phó Sự cố, Thiên tai và Tìm kiếm Cứu nạn trong việc phân bổ lực lượng tàu kéo, phao quây dầu chuyên dụng tại các vùng biển xung yếu.
  • Lợi ích xã hội và môi trường: Góp phần bảo vệ sinh kế cho hàng triệu ngư dân và bảo tồn đa dạng sinh học tại các khu bảo tồn biển trọng điểm như Cù Lao Chàm, Nha Trang, Phú Quốc.

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

  • Nghiên cứu sinh và Giảng viên ngành Hải dương học/Môi trường: Tiếp cận một cấu trúc mô hình số trị Euler hoàn chỉnh từ phương trình toán học đến sơ đồ thuật toán sai phân dưới lưới, làm tài liệu tham khảo chuyên sâu cho các bài toán phân tán chất.
  • Cơ quan Quản lý Hàng hải và Cảnh sát Biển: Sở hữu phương pháp luận khoa học để thực hiện điều tra pháp y môi trường biển (marine environmental forensics), truy vết và quy trách nhiệm pháp lý đối với các phương tiện vận tải quốc tế xả thải dầu lén lút.
  • Doanh nghiệp Vận tải Biển và Cảng Biển: Tối ưu hóa quy trình quản lý rủi ro và diễn tập ứng phó sự cố tràn dầu tại các trung tâm cảng biển lớn như Hải Phòng, Cái Mép - Thị Vải, Vân Phong.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc thiết lập thành công cơ sở toán - vật lý cho bài toán ngược thời gian lan truyền dầu theo tiếp cận Euler. Nghiên cứu đã mở rộng lý thuyết bình lưu - khuếch tán kinh điển và lý thuyết bảo toàn vật chất bằng cách chứng minh rằng: dưới điều kiện biên ràng buộc không âm và trường vận tốc động lực học đảo chiều $(-\mathbf{V}_a)$, mô hình Euler có khả năng tích tụ và tái lập điểm kỳ dị nguồn thải ban đầu mà không bị phá vỡ tính ổn định số trị.

2. Đổi mới phương pháp luận của luận án khi so sánh với các nghiên cứu quốc tế?

So với mô hình hạt ngược BITT (Batchelder, 2006; Christensen nnk., 2007) vốn chỉ theo dõi quỹ đạo hạt rời rạc với độ chính xác nguồn cực thấp ($0,1%$), và phương pháp LCSs của Peacock nnk. (MIT) chỉ xác định các ranh giới vận chuyển biểu kiến, mô hình Euler của luận án giải quyết trực tiếp bài toán trường nồng độ liên tục $C(x,y,t)$ kết hợp sơ đồ sai phân dưới lưới $2^{-8}$, cho phép tái tạo đồng thời cả vị trí nguồn, thời điểm phát thải và khối lượng dầu rò rỉ ban đầu.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là trong sự cố tràn dầu diện rộng năm 2007, hướng trôi dạt của dầu vào bờ biển Nam Trung Bộ không đơn thuần do gió mùa Đông Bắc đẩy thẳng vào bờ, mà do sự tương tác cộng hưởng giữa dòng chảy ven bờ tốc độ cao ép sát bờ biển miền Trung và cấu trúc xoáy nghịch ngoài khơi, tạo thành "hành lang vận chuyển siêu tốc" đưa dầu đi xa hơn 800 km chỉ trong vòng 10 ngày mà diện tích vệt dầu bề mặt không bị tan rã hoàn toàn nhờ quá trình nhũ tương hóa làm tăng độ nhớt lên hơn $1.200\text{ cP}$.

4. Luận án có cung cấp quy trình tái lặp (replication protocol) không?

Luận án cung cấp quy trình tái lặp chi tiết tuyệt đối: từ sơ đồ khối chương trình chính, hệ thống chương trình con tính toán sai phân bậc 1 và bậc 2 (Lax-Wendroff), bộ tham số chuẩn hóa 8 quá trình phong hóa, ma trận chuyển đổi góc lệch gió $D_w$, đến hệ tiêu chuẩn đánh giá sai số (NASH, PBIAS, RMSE) và các điều kiện biên cứng/lỏng được chuẩn hóa.

5. Khung chương trình nghiên cứu 10 năm được vạch ra như thế nào?

Chiến lược 10 năm tập trung vào 3 trụ cột: (1) Nâng cấp mô hình lên hệ thống tích hợp 3D Hydrodynamic-Spill-Ecosystem Model; (2) Tích hợp trí tuệ nhân tạo (AI-driven inverse solver) kết hợp viễn thám vệ tinh SAR thời gian thực; (3) Mở rộng phạm vi áp dụng cho toàn bộ vùng biển Đông Nam Á và thiết lập mạng lưới chia sẻ dữ liệu cảnh báo sớm quốc tế.

Kết luận

  1. Xây dựng thành công hệ thống mô hình số trị Euler 2D mô phỏng hoàn chỉnh quá trình lan truyền và 8 cơ chế phong hóa hóa-lý của dầu tràn trên biển, khắc phục triệt để sự thiếu đồng bộ giữa mô hình vết dầu và mô hình thủy động lực môi trường.
  2. Tiên phong giải quyết bài toán nghịch đảo thời gian theo tiếp cận Euler phục vụ truy tìm nguồn gốc ô nhiễm dầu, đạt độ chính xác cao về tọa độ nguồn thải ($R > 0,85$; sai số vị trí $< 14\text{ km}$ sau 7 ngày mô phỏng ngược).
  3. Phát triển kỹ thuật sai phân lồng ghép dưới lưới ($2^{-4}, 2^{-8}$), giải quyết thành công bài toán triệt tiêu sai số khuếch tán số học, nâng cao chỉ số hiệu quả mô hình NASH từ 0,52 lên 0,88.
  4. Giải mã khoa học nguồn gốc sự cố dầu tràn bí ẩn năm 2007 - 2008 trên vùng biển Việt Nam, chứng minh dầu phát tán từ tuyến vận tải hàng hải quốc tế ngoài khơi Biển Đông trong điều kiện gió mùa Đông Bắc và Tây Nam điển hình.
  5. Mở ra hướng nghiên cứu mới trong hải dương học ứng dụng và pháp y môi trường biển, cung cấp công cụ khoa học chuẩn xác phục vụ bảo vệ chủ quyền, an ninh môi trường và phát triển kinh tế biển bền vững của Việt Nam.