Tổng quan về luận án

Sự cố tràn dầu trên môi trường biển là một trong những thảm họa sinh thái và kinh tế nghiêm trọng nhất đối với các quốc gia ven biển. Biển Đông là vùng biển có mật độ giao thông hàng hải quốc tế nhộn nhịp thứ hai thế giới (sau Địa Trung Hải) với hơn 300–400 triệu tấn dầu thô được vận chuyển mỗi năm qua các tuyến hàng hải huyết mạch từ eo biển Malacca đến Đông Á. Bên cạnh đó, các bồn trũng dầu khí lớn như Nam Côn Sơn, Cửu Long, Sông Hồng, Malay – Thổ Chu với trữ lượng khai thác hàng tỷ mét khối và hàng trăm giàn khoan hoạt động liên tục tiềm ẩn rủi ro ô nhiễm rất cao. Thực tế lịch sử ghi nhận từ năm 1990 đến năm 2015 đã xảy ra hơn 50 sự cố tràn dầu lớn tại vùng biển Việt Nam, đặc biệt là hiện tượng dầu tràn vô chủ dạt vào bờ biển miền Trung và miền Nam trong hai năm 2007 và 2008 với khối lượng thu gom vượt trên 5.000 tấn, gây thiệt hại nghiêm trọng cho hệ sinh thái và kinh tế biển.

                     ┌───────────────────────────────────────────────┐
                     │   BÀI TOÁN TOÀN DIỆN VỀ LAN TRUYỀN DẦU TRÀN   │
                     └───────────────────────┬───────────────────────┘
                                             │
                    ┌────────────────────────┴────────────────────────┐
                    ▼                                                 ▼
      ┌───────────────────────────┐                     ┌───────────────────────────┐
      │   QUÁ TRÌNH XUÔI THỜI GIAN│                     │  QUÁ TRÌNH NGƯỢC THỜI GIAN│
      │    (Forward-in-Time)      │                     │     (Backward-in-Time)    │
      ├───────────────────────────┤                     ├───────────────────────────┤
      │ • Dự báo quỹ đạo loang    │                     │ • Truy tìm nguồn gốc thải │
      │ • Cảnh báo vùng nguy cơ   │                     │ • Tái thiết lập vị trí sự cố│
      │ • Mô phỏng biến đổi lý hóa│                     │ • Hỗ trợ pháp lý & ứng phó│
      └───────────────────────────┘                     └───────────────────────────┘

Trong bối cảnh đó, luận án tiến sĩ chuyên ngành Hải dương học (mã số: 62 44 02 28) 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" do nghiên cứu sinh Nguyễn Quốc Trinh thực hiện dưới sự hướng dẫn khoa học của PGS. Nguyễn Minh Huấn và PGS. Phùng Đăng Hiếu tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội (2018), đã đặt nền móng tiên phong trong việc giải quyết bài toán kép: dự báo lan truyền (xuôi thời gian) và truy tìm vị trí nguồn thải vô chủ (ngược thời gian) bằng phương pháp tiếp cận Euler.

Khoảng trống nghiên cứu (Research Gap)

Trước khi luận án được công bố, các mô hình mô phỏng dầu tràn trong nước và quốc tế chủ yếu tập trung vào bài toán dự báo xuôi thời gian dựa trên cách tiếp cận quỹ đạo hạt Lagrange (Lagrangian Particle-Tracking Models - LPTMs như GNOME của NOAA, MIKE21/31 SA/PA của DHI, hay OilSAS, OST-2D/3D). Đối với bài toán ngược thời gian nhằm xác định vị trí nguồn thải vô chủ:

  1. Các nghiên cứu quốc tế (như Batchelor, 1952; Peacock et al., 2010 với lý thuyết cấu trúc nhánh kết hợp Lagrangian - LCSs; hay BITT trong dự báo nguồn sinh vật phù du của Christensen et al., 2007) hầu như chỉ áp dụng Lagrange ngẫu hành không bảo toàn tính phân tán thực hoặc đạt hiệu suất hội tụ rất thấp (chỉ ~0,1% như kiểm nghiệm tại OMEX, Na Uy).
  2. Tại Việt Nam, các nghiên cứu của Vũ Thanh Ca và nnk. (2007), Trần Hồng Thái và nnk. (2008), hay Nguyễn Đình Dương và nnk. (2008) chỉ dừng lại ở phương pháp thử sai trực quan (trial-and-error) hoặc đảo dấu thời gian sơ khai trên mô hình hạt, thiếu một khung toán học chặt chẽ dựa trên phương trình vi phân đạo hàm riêng kết hợp đầy đủ các quá trình phong hóa phi tuyến tính.
  3. Thiếu sự đồng nhất về mặt phương pháp học giữa dữ liệu trường động lực môi trường (vốn xuất phát từ lưới Euler của mô hình khí tượng/thủy văn) với mô hình lan truyền chất ô nhiễm.

Câu hỏi nghiên cứu và Giả thuyết khoa học

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu cốt lõi:

  • RQ1: Làm thế nào để xây dựng một mô hình số trị Euler 2 chiều giải đồng thời phương trình truyền - khuếch tán dầu kết hợp hệ thống tham số hóa 8 quá trình phong hóa hóa lý trên bề mặt biển?
  • RQ2: Cơ sở toán - lý nào cho phép đảo ngược dòng thời gian ($t \rightarrow -t$) trong môi trường liên tục để tích tụ vết dầu tiêu tán ngược về tọa độ điểm xả thải ban đầu mà vẫn đảm bảo tính ổn định số và điều kiện nghiệm không âm?
  • RQ3: Mức độ chính xác và độ nhạy của mô hình Euler xuôi - ngược khi tích hợp với trường gió tái phân tích khí hậu CFSR, dòng chảy hải dương POM và sóng SWAN trên toàn bộ không gian Biển Đông là bao nhiêu?

Hệ thống giả thuyết kiểm định:

  • H1: Mô hình Euler tích hợp đầy đủ các trường lực (gió bề mặt $V_w$, dòng chảy tầng mặt $V_c$, sóng Stokes $V_s$) sẽ mô phỏng vận tốc trôi dạt tổng hợp $V_a$ chính xác hơn đáng kể so với việc chỉ xét dòng chảy hoặc gió đơn lẻ.
  • H2: Sơ đồ sai phân kết hợp đa tầng dưới lưới (sub-grid scheme bậc 2-4 và 2-8) sẽ triệt tiêu hiệu ứng khuếch tán số (numerical diffusion), cho phép bảo toàn khối lượng và độ tập trung nồng độ cực đại của vết dầu khi tích phân ngược thời gian.
  • H3: Vệt dầu vô chủ trôi dạt vào bờ biển Việt Nam năm 2007 và 2008 có thể xác định được vùng phát thải ban đầu ngoài khơi thông qua việc chạy mô hình ngược thời gian trong khoảng thời gian từ 5 đến 7 ngày trước khi phát hiện.

Khung lý thuyết và Phạm vi nghiên cứu

Nghiên cứu được xây dựng dựa trên sự tích hợp của Lý thuyết Thủy động lực học Hải dương học (Navier-Stokes/Ekman Dynamics), Định luật Khuếch tán Fickian, Định luật Bảo toàn Vật chất, và Nguyên lý Động học Phong hóa Dầu mỏ (Mackay, 1980; ASCE, 1996; Lehr et al., 2002).

  • Phạm vi không gian: Toàn bộ bồn trũng Biển Đông trong giới hạn tọa độ địa lý từ $1^\circ\text{N}$ đến $24^\circ\text{N}$ và từ $99^\circ\text{E}$ đến $121^\circ\text{E}$.
  • Phạm vi thời gian và kịch bản: Kiểm nghiệm trên điều kiện môi trường đồng nhất lý thuyết, trường môi trường trung bình tháng nhiều năm, và trường biến đổi thực tế theo các mốc quan trắc sự cố tràn dầu lịch sử (tiêu biểu là các mốc tháng 02/2007 và tháng 07/2008).

Literature Review và Positioning

Tổng quan các trường phái nghiên cứu dầu tràn

Lịch sử mô phỏng tràn dầu trên thế giới khởi đầu từ các nghiên cứu thực nghiệm bán giải tích của Fay (1969, 1971) và Hoult (1972) về ba giai đoạn loang dầu cơ học dưới tác dụng của trọng lực, quán tính, độ nhớt và sức căng bề mặt. Bước sang thập niên 1980, Mackay và nnk. (1980, 1982) đã đặt nền móng định lượng cho các quá trình phong hóa (weathering) bao gồm bay hơi, nhũ tương hóa và phân tán thẳng đứng.

Hệ thống mô hình hóa toán học sau đó phát triển thành hai trường phái đối nghịch:

  1. Trường phái tiếp cận Lagrange (Lagrangian Approach): Coi khối dầu tràn là tập hợp của hàng vạn chất điểm/hạt dầu độc lập di chuyển theo trường vận tốc và khuếch tán ngẫu hành Monte Carlo (Random Walk). Các mô hình tiêu biểu gồm GNOME (NOAA, Mỹ), MEDSLIK (Lardner et al., 1998, Síp), MIKE31 SA/PA (DHI, Đan Mạch), OSCAR (Reed et al., 1995, Na Uy). Ưu điểm của hướng này là tính toán nhanh, không phụ thuộc vào lưới không gian cố định, nhưng nhược điểm chí mạng là phải nội/ngoại suy trường môi trường từ lưới Euler sang vị trí hạt, gây mất tính liên tục và sai số lớn ở vùng ven bờ hoặc vùng gradient dòng chảy cao.
  2. Trường phái tiếp cận Euler (Eulerian Approach): Thiết lập phương trình truyền - khuếch tán nồng độ liên tục $C(x, y, t)$ và độ dày lớp dầu $h(x, y, t)$ trên hệ tọa độ không gian cố định. Tiêu biểu là các mô hình OILMAP (Spaulding et al., 1994), MOHID (Trancoso et al., 2003), GULFSPILL (Al-Rabeh et al., 1989). Ưu điểm vượt trội là sự đồng nhất tuyệt đối về cấu trúc dữ liệu với các mô hình khí tượng và thủy động lực học hải dương, bảo toàn trường thông lượng vật chất trên từng ô lưới.
                      SO SÁNH BÀI TOÁN NGƯỢC THỜI GIAN QUỐC TẾ
┌──────────────────────────────────┬──────────────────────────────────┐
│     TIẾP CẬN LAGRANGE (BITT)     │      TIẾP CẬN EULER (LUẬN ÁN)    │
├──────────────────────────────────┼──────────────────────────────────┤
│ • Mô hình: LPTMs, BITT           │ • Mô hình: Euler 2D ngược        │
│ • Kiểm nghiệm: OMEX (Na Uy)      │ • Ứng dụng: Biển Đông (Việt Nam) │
│ • Độ hội tụ: Rất thấp (~0.1%)    │ • Chỉ số: NASH > 0.85; R > 0.95  │
│ • Nhược điểm: Phân kỳ ngẫu hành  │ • Ưu điểm: Bảo toàn nồng độ &    │
│   khi đảo ngược thời gian        │   hóa lý phong hóa nghịch        │
└──────────────────────────────────┴──────────────────────────────────┘

Tranh luận học thuật về bài toán truy tìm nguồn gốc ngược thời gian

Vấn đề giải bài toán ngược thời gian (Backward-in-time inversion) luôn là thách thức lớn nhất trong động lực học chất lỏng do tính bất đối xứng nhiệt động lực học và entropy tăng (quá trình khuếch tán phân tán vật chất theo thời gian là quá trình không thuận nghịch tự nhiên).

  • Một số nhà nghiên cứu (như Dimou & Adams, 1993; Spivakovskaya et al., 2007) cho rằng mô hình Lagrange đảo ngược thời gian bằng cách đổi dấu vector vận tốc ($V = -V$) sẽ dẫn đến sự phân kỳ vô hạn của thuật toán ngẫu hành vì tham số ngẫu nhiên không thể đảo ngược chiều xác suất.
  • Ngược lại, trường phái giải bài toán biên ngược (Inverse Problem / Adjoint Model) của Marchuk (1995) và Bocquet (2005) trong khí tượng chứng minh rằng nếu sử dụng phương trình đạo hàm riêng Euler với toán tử đối ngẫu, nồng độ phân tán có thể được quy tụ ngược về nguồn nếu xử lý triệt để tiêu tán số.

Định vị đóng góp học thuật của luận án

So với các nghiên cứu quốc tế:

  • So sánh 1: Nghiên cứu của Christensen et al. (2007) sử dụng mô hình BITT (Backward-in-Time Trajectory) để tìm nguồn cá đẻ tại vùng biển Na Uy chỉ đạt xác suất hội tụ danh định 0,1% do hiệu ứng trôi dạt hỗn loạn không bảo toàn. Luận án của NCS. Nguyễn Quốc Trinh đã khắc phục nhược điểm này bằng cách sử dụng toán tử Euler kết hợp trường lực trôi dạt tổng hợp thực nghiệm, cho phép tái hội tụ trường nồng độ cực đại với hệ số Nash-Sutcliffe đạt trên 0,85.
  • So sánh 2: Nghiên cứu vật thể trôi của Breivik & Allen (2008) tại Biển Bắc chỉ áp dụng phương pháp Monte Carlo thuần túy không xét quá trình biến đổi khối lượng. Luận án tích hợp trực tiếp 8 quá trình phong hóa nghịch, khôi phục lại thể tích $V_0$ và mật độ $\rho_0$ ban đầu của vết dầu trước khi bị bốc hơi và nhũ tương hóa.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho hệ thống lý thuyết

Luận án đã mở rộng và hoàn thiện hệ thống lý thuyết lan truyền chất ô nhiễm trên biển qua các đóng góp cụ thể:

  1. Mở rộng lý thuyết trôi dạt bề mặt Ekman - Stokes: Thay vì sử dụng hệ số trôi dạt gió cố định (thường lấy 3% theo kinh nghiệm cổ điển), luận án thiết lập hàm phi tuyến tính động cho hệ số ma sát gió $\alpha_w$, ma trận quay góc lệch gió - dòng $\mathbf{D}_w$, hệ số dòng chảy $\alpha_c$ và hệ số sóng Stokes $\alpha_s$, phản ánh chính xác tương tác phi tuyến tại tầng biên đại dương - khí quyển.
  2. Thiết lập khung toán học bài toán ngược thời gian trên lưới Euler: Luận án đề xuất tiên đề bảo toàn ngược: Quá trình phân tán tiêu tán xuôi thời gian tương đương với quá trình tích tụ hội tụ ngược thời gian dưới tác dụng của trường lực liên hợp đảo chiều: $$\vec{V}_{\text{reverse}} = -\vec{V}_a, \quad \text{với } \vec{V}_a = \alpha_w \mathbf{D}_w \vec{V}_w + \alpha_c \vec{V}_c + \alpha_s \vec{V}_s$$
  3. Mô hình hóa động học phong hóa thuận nghịch: Xây dựng hệ phương trình cân bằng liên tục cho phép hiệu chỉnh ngược các đại lượng thể tích $V(t)$, khối lượng $M(t)$, khối lượng riêng $\rho(t)$, độ nhớt động học $\nu(t)$ và sức căng bề mặt $\sigma(t)$ từ trạng thái dầu phong hóa ngoài hiện trường trở về trạng thái nguyên khai tại nguồn xả.

Khung phân tích độc đáo và Hệ phương trình cấu vi mô

Mô hình toán học của luận án tích hợp 3 trụ cột lý thuyết chính: Thủy động lực học chất lỏng địa cầu, Động học hóa lý dầu mỏ và Giải thuật sai phân hữu hạn Euler.

                           HỆ THỐNG PHƯƠNG TRÌNH ĐỘNG LỰC HỌC
┌──────────────────────────────────────────────────────────────────────────────────┐
│  PHƯƠNG TRÌNH XUÔI THỜI GIAN (FORWARD EULER):                                    │
│  ∂C/∂t + ∂(Vx·C)/∂x + ∂(Vy·C)/∂y = ∂(Dx·∂C/∂x)/∂x + ∂(Dy·∂C/∂y)/∂y + Cm          │
│                                                                                  │
│  PHƯƠNG TRÌNH NGƯỢC THỜI GIAN (BACKWARD EULER):                                  │
│  ∂C/∂t + ∂(-Vx·C)/∂x + ∂(-Vy·C)/∂y = ∂(Dx·∂C/∂x)/∂x + ∂(Dy·∂C/∂y)/∂y + Cm        │
└──────────────────────────────────────────────────────────────────────────────────┘

1. Hệ phương trình lan truyền nồng độ và độ dày 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) + C_m$$

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

Trong đó:

  • $C(x, y, t)$: Nồng độ diện tích của lớp dầu trên mặt nước $[\text{kg}\cdot\text{m}^{-2}]$.
  • $h(x, y, t)$: Độ dày màng dầu $[\text{m}]$, liên hệ qua $C = \rho \cdot h$.
  • $V_x, V_y$: Các thành phần vận tốc trôi dạt tổng hợp theo phương kinh hướng và vĩ hướng $[\text{m}\cdot\text{s}^{-1}]$.
  • $D_x, D_y$: Hệ số khuếch tán rối ngang $[\text{m}^2\cdot\text{s}^{-1}]$, tham số hóa theo quy mô xoáy Okubo và độ dày vệt dầu.
  • $C_m$: Suất nguồn dầu xả thải $[\text{kg}\cdot\text{m}^{-2}\cdot\text{s}^{-1}]$.

2. Hệ phương trình tích tụ nồng độ ngược 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) + C_m$$

3. Mô hình tham số hóa 8 quá trình phong hóa:

Thể tích vết dầu $V(t)$ tại thời điểm $t$ bị suy giảm hoặc hiệu chỉnh theo công thức tổng quát: $$V = V_0 \left( 1 - F_e - F_d - F_{disc} - F_{sed} - F_{coast} - F_{oxy} - F_{bio} \right) (1 + F_w)$$

  • Tỷ phần bay hơi ($F_e$): Hàm phụ thuộc nhiệt độ sôi, áp suất hơi bão hòa, gió $V_w$ và diện tích màng dầu theo Mackay (1980).
  • Tỷ phần nhũ tương hóa ($F_w$): Tốc độ hấp thụ nước tạo nhũ tương nước-trong-dầu ("chocolate mousse"): $\frac{dF_w}{dt} = K_a (V_w + 1)^2 (1 - K_b F_w)$.
  • Biến đổi khối lượng riêng ($\rho$): $$\rho = F_w \rho_w + (1 - F_w) \rho_0 [1 + C_1 (T - T_e)][1 + C_2 F_e]$$ (với $C_1 = 8{,}0 \times 10^{-3}\text{ K}^{-1}$, $C_2 = 0{,}18$, $\rho_w$ là khối lượng riêng nước biển tầng mặt).
  • Biến đổi độ nhớt động lực học ($\mu$): $$\mu = \mu_0 \exp\left[ C_6 \left( \frac{1}{T} - \frac{1}{T_e} \right) + C_5 F_e + \frac{C_3 F_w}{1 - C_4 F_w} \right]$$
  • Sức căng bề mặt dầu - nước ($\sigma_w$) và dầu - không khí ($\sigma_A$): $$\sigma_w = \sigma_{w0}(1 + F_e), \quad \sigma_A = \sigma_{A0}(1 - F_e)$$
  • Phân tán thẳng đứng ($F_d$), hòa tan ($F_{disc}$), hấp thụ trầm tích đáy ($F_{sed}$), tương tác bờ bãi ($F_{coast}$), oxy hóa quang hóa ($F_{oxy}$) và phân hủy sinh học ($F_{bio}$) được số hóa thành các module thứ cấp phụ thuộc vào trường bức xạ mặt trời, độ đục trầm tích lơ lửng và đường bờ.

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

Thiết kế nghiên cứu và Triết lý khoa học

Nghiên cứu tuân thủ chặt chẽ triết lý Thực chứng (Positivism) kết hợp phương pháp mô phỏng số tất định (Deterministic Numerical Simulation). Toàn bộ hệ thống phương trình vi phân đạo hàm riêng phi tuyến được giải số dựa trên phương pháp sai phân hữu hạn (Finite Difference Method) trên lưới không gian 2D rời rạc hóa có lồng ghép lưới mịn dưới ô (nested sub-grid).

                      SƠ ĐỒ TÍCH HỢP HỆ THỐNG DỮ LIỆU ĐẦU VÀO
┌───────────────────────────┐  ┌───────────────────────────┐  ┌───────────────────────────┐
│        GIÓ KHÍ QUYỂN      │  │     DÒNG CHẢY BIỂN        │  │        SÓNG MẶT           │
│   Tái phân tích NCEP/CFSR │  │    Mô hình thủy động POM  │  │    Mô hình sóng SWAN      │
│     Độ phân giải 0.3°     │  │      Độ phân giải 1/12°   │  │    Độ phân giải 0.1°      │
└─────────────┬─────────────┘  └─────────────┬─────────────┘  └─────────────┬─────────────┘
              │                              │                              │
              └──────────────────────┬───────┴──────────────────────────────┘
                                     ▼
                      ┌─────────────────────────────┐
                      │    MÔ HÌNH DẦU TRÀN EULER   │
                      │  (Bước lưới: 0.01° - 0.06°) │
                      └─────────────────────────────┘

Quy trình nghiên cứu và Nguồn dữ liệu tích hợp

Hệ thống dữ liệu phục vụ mô phỏng quy mô toàn Biển Đông bao gồm:

  1. Dữ liệu địa hình đáy biển: Trích xuất từ cơ sở dữ liệu cao độ quốc gia kết hợp ETOPO1 với độ sâu từ $0\text{ m}$ đến trên $4.000\text{ m}$.
  2. Dữ liệu khí tượng tầng mặt: Trường gió 10m ($V_w$), nhiệt độ không khí ($T_a$) trích xuất từ dữ liệu tái phân tích hoàn lưu khí quyển toàn cầu CFSR (Climate Forecast System Reanalysis) của NCEP/NOAA với độ phân giải thời gian 1 giờ, không gian $0{,}3^\circ \times 0{,}3^\circ$.
  3. Dữ liệu hải dương học: Trường dòng chảy mặt ($V_c$), nhiệt độ nước biển ($T_w$), độ muối ($S$) thu nhận từ mô hình hoàn lưu đại dương Princeton Ocean Model (POM) thiết lập riêng cho Biển Đông với độ phân giải $1/12^\circ$.
  4. Dữ liệu sóng biển: Trường sóng mặt trích xuất từ mô hình mô phỏng sóng ven bờ SWAN (Simulating WAves Nearshore), cung cấp tham số chiều cao sóng có ý nghĩa $H_s$, chu kỳ đỉnh $T_p$ để tính dòng trôi Stokes: $$V_s \approx \frac{1}{16} \pi H_s^2 \frac{g}{T_p}$$
  5. Dữ liệu ảnh vệ tinh kiểm chứng: Bản đồ ảnh viễn thám quang học và radar khẩu độ tổng hợp (SAR) phát hiện vệt dầu nổi trên Biển Đông năm 2007 và 2008 của Cục Viễn thám Quốc gia và Viện Hàn lâm Khoa học và Công nghệ Việt Nam.

Kỹ thuật số trị và Đánh giá độ ổn định

  • Rời rạc hóa không gian: Thử nghiệm phổ bước lưới đa tỷ lệ $dx = dy \in {0{,}01^\circ; 0{,}02^\circ; 0{,}03^\circ; 0{,}04^\circ; 0{,}044^\circ; 0{,}05^\circ; 0{,}06^\circ}$, tương đương bước không gian từ $\sim 1{,}1\text{ km}$ đến $\sim 6{,}6\text{ km}$. Bước lưới tối ưu được lựa chọn cho tính toán toàn vùng là $dx = dy = 0{,}044^\circ$ ($\approx 5\text{ km}$).
  • Thuật toán xử lý dưới lưới (Sub-grid): Sử dụng sơ đồ sai phân lồng bậc 2-4 và bậc 2-8 nhằm giảm thiểu tiêu tán số, kiểm soát sai số cắt cụt (truncation error).
  • Bộ chỉ số thống kê định lượng đánh giá sai số:
    • Hệ số hiệu quả Nash-Sutcliffe (NASH): Đánh giá mức độ khớp giữa quỹ đạo mô phỏng và thực tế ($-\infty < \text{NASH} \le 1$): $$\text{NASH} = 1 - \frac{\sum (O_i - P_i)^2}{\sum (O_i - \bar{O})^2}$$
    • Phần trăm độ lệch (PBIAS): Đo lường xu thế đánh giá quá cao hoặc quá thấp của mô hình ($%$): $$\text{PBIAS} = \frac{\sum (O_i - P_i)}{\sum O_i} \times 100$$
    • Căn sai phương trung bình chuẩn hóa (RMSE): Đánh giá độ lệch tuyệt đối về tọa độ tâm và nồng độ vệt dầu: $$\text{RMSE} = \sqrt{\frac{1}{N}\sum (O_i - P_i)^2}$$
    • Hệ số tương quan tuyến tính Pearson ($R$).

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

Những phát hiện then chốt từ thực nghiệm mô phỏng

                       HIỆU SUẤT TRUY VẾT NGƯỢC THEO THỜI GIAN
┌───────────────────────────┬───────────────────────────┬───────────────────────────┐
│     SỰ CỐ THÁNG 02/2007   │     SỰ CỐ THÁNG 07/2008   │      CHỈ SỐ THỐNG KÊ      │
├───────────────────────────┼───────────────────────────┼───────────────────────────┤
│ • Gió: Đông Bắc (CFSR)    │ • Gió: Tây Nam (CFSR)     │ • NASH: 0.88 - 0.94       │
│ • Truy ngược: 5 - 7 ngày  │ • Truy ngược: 5 - 7 ngày  │ • PBIAS: < 6.5%           │
│ • Vị trí: Vùng mỏ ngoài   │ • Vị trí: Tuyến hàng hải  │ • R: 0.96 - 0.98          │
│   khơi Nam Trung Bộ       │   quốc tế Đông Nam Bộ     │ • RMSE: < 0.08 kg/m²      │
└───────────────────────────┴───────────────────────────┴───────────────────────────┘

1. Đóng góp vượt trội của cơ chế trôi dạt tổng hợp

Kết quả kiểm nghiệm cho thấy nếu chỉ sử dụng dòng chảy đơn thuần (mô hình thủy động thông thường), sai số vị trí tâm vệt dầu sau 48 giờ lên tới $42%$. Khi tích hợp đầy đủ công thức trôi dạt ba thành phần: $$V_a = \alpha_w \mathbf{D}_w V_w + \alpha_c V_c + \alpha_s V_s$$ với góc lệch $\theta = 45^\circ \exp(-V_w / 12{,}5)$ và $\alpha_w = 0{,}05[1 - \exp(-V_w / 17{,}5)]$, chỉ số NASH tăng từ 0,61 lên 0,93; hệ số tương quan $R$ đạt 0,97. Gió bề mặt và sóng Stokes đóng góp từ $55%$ đến $70%$ vận tốc dịch chuyển thực tế của màng dầu tầng mặt.

2. Kiểm soát thành công hiệu ứng khuếch tán số xuôi - ngược

Thử nghiệm trên các bước lưới từ $0{,}01^\circ$ đến $0{,}06^\circ$ cho thấy khi tăng kích thước ô lưới, diện tích vệt dầu bị khuếch tán nhân tạo mở rộng phi thực tế. Tuy nhiên, việc áp dụng sơ đồ sai phân kết hợp dưới lưới 2-8 đã triệt tiêu hiện tượng này:

  • Nồng độ cực đại của vệt dầu $C_{\max}$ được bảo toàn với sai số $< 4{,}8%$ sau 5 ngày mô phỏng liên tục.
  • Tỷ lệ quãng đường di chuyển giữa mô hình Euler và nghiệm giải tích phương trình bảo toàn đạt tỷ số lý tưởng xấp xỉ 1,00 (sai số $< 1{,}5%$).

3. Khả năng truy tìm chính xác nguồn phát thải vô chủ (Ngược thời gian)

Khi áp dụng giải bài toán ngược thời gian cho 2 sự cố tràn dầu thực tế trên Biển Đông:

  • Sự cố dầu tràn tháng 02/2007 (chế độ gió mùa Đông Bắc): Mô hình tính toán ngược từ ngày 25/02/2007 lùi 7 ngày về ngày 18/02/2007 đã tái hiện lại quá trình tích tụ từ vệt dầu loang rộng hàng ngàn $\text{km}^2$ ven bờ Nam Trung Bộ thu nhỏ dần về một điểm xả tập trung tại tọa độ ngoài khơi khu vực bể Cửu Long/Nam Côn Sơn, trùng khớp với tuyến hàng hải quốc tế tiếp cận vùng mỏ.
  • Sự cố dầu tràn tháng 07/2008 (chế độ gió mùa Tây Nam): Mô hình ngược thời gian từ ngày 25/07/2008 lùi 7 ngày đã chỉ ra nguồn phát thải xuất phát từ khu vực biển phía Tây Nam quần đảo Trường Sa, nơi có mật độ tàu dầu di chuyển với tải lượng lớn qua eo biển Bashi – Malacca.
  • Chỉ số kiểm định chất lượng quá trình ngược thời gian đạt: $\text{NASH} = 0{,}88 - 0{,}94$; $\text{PBIAS} < 6{,}5%$; $R > 0{,}96$.

4. Phân tích tính ổn định thuận nghịch (Euler Round-Trip Stability)

Luận án thực hiện phép thử kép: chạy mô hình xuôi thời gian 7 ngày từ điểm nguồn, sau đó lấy kết quả trường nồng độ cuối cùng làm điều kiện ban đầu để chạy ngược thời gian 7 ngày về quá khứ (và ngược lại). Kết quả chứng minh tâm khối lượng của vệt dầu quy tụ về điểm ban đầu với sai số khoảng cách dưới 1 bước lưới ($\Delta L < 5\text{ km}$ trên miền tính $2.500\text{ km}$), khẳng định tính đối xứng và độ tin cậy cơ học của thuật toán.

Ý nghĩa và Tác động đa chiều (Implications)

  • Ý nghĩa học thuật: Khẳng định tính khả thi của phương pháp Euler trong bài toán ngược động lực học chất ô nhiễm biển; thiết lập tiêu chuẩn tham số hóa phong hóa dầu tại vùng biển nhiệt đới Biển Đông (nhiệt độ cao đẩy mạnh tốc độ bay hơi trong 24 giờ đầu lên tới $35 - 45%$).
  • Ứng dụng thực tiễn: Cung cấp bộ công cụ phần mềm số trị độc lập, không phụ thuộc vào các gói thương mại đắt tiền nước ngoài (như MIKE của DHI), cho phép Trung tâm Quốc gia Ứng phó Sự cố Tràn dầu Việt Nam vận hành cảnh báo sớm và điều tra nguyên nhân tai nạn tràn dầu.
  • Chính sách và Pháp lý: Kết quả truy vết ngược thời gian là bằng chứng khoa học đanh thép phục vụ công tác đấu tranh pháp lý quốc tế, xác định trách nhiệm bồi thường thiệt hại môi trường đối với các tàu chở dầu quốc tế xả thải trộm hoặc gây tai nạn bỏ trốn trên vùng đặc quyền kinh tế (EEZ) của Việt Nam.

Limitations và Future Research

Hạn chế nghiên cứu

  1. Giả định mô hình 2 chiều tầng mặt: Mô hình hiện tại tập trung giải quyết phương trình 2D cho lớp màng dầu nổi trên mặt nước ($C$ và $h$), chưa mô hình hóa đầy đủ trường 3D liên tục của các giọt dầu phân tán chìm sâu vào cột nước dưới tác động của dòng chảy thẳng đứng $w_p$ và hiện tượng tái nổi (resurfacing).
  2. Độ phân giải dữ liệu môi trường đầu vào: Dữ liệu gió tái phân tích CFSR ($0{,}3^\circ$) và dòng chảy POM ($1/12^\circ$) tuy bao phủ toàn diện Biển Đông nhưng vẫn chưa phản ánh chi tiết hoàn lưu ven bờ, dòng triều cục bộ tại các vịnh nông và đầm phá phức tạp.
  3. Tham số hóa tương tác bãi bờ: Quá trình bám dính của dầu vào các dạng địa hình bờ khác nhau (rừng ngập mặn, bãi cát, rạn san hô, vách đá) mới dừng ở hệ số suy giảm thực nghiệm bậc một, chưa mô phỏng chi tiết cơ chế động lực sóng vùng sóng vỡ (surf zone).

Định hướng phát triển tương lai

  • Phát triển mô hình Euler-Lagrange 3D tích hợp: Kết nối mô hình màng dầu Euler tầng mặt với mô hình hạt Lagrange 3D trong cột nước để mô phỏng toàn diện số phận dầu từ đáy biển (rò rỉ giếng khoan sâu như thảm họa Deepwater Horizon) lên mặt biển.
  • Tích hợp đồng hóa dữ liệu vệ tinh thời gian thực (Data Assimilation): Ứng dụng kỹ thuật lọc Kalman (EnKF) để cập nhật trực tiếp diện tích và nồng độ vết dầu từ ảnh vệ tinh radar Sentinel-1, Cosmo-SkyMed vào mô hình số nhằm hiệu chỉnh tức thời điều kiện biên bài toán ngược.
  • Mở rộng thuật toán trí tuệ nhân tạo (AI/Physics-Informed Neural Networks - PINNs): Sử dụng mạng nơ-ron học sâu tích hợp vật lý để tăng tốc độ tính toán bài toán ngược, phục vụ tìm kiếm nguồn xả thải tức thời trong vòng vài giây hỗ trợ lực lượng Cảnh sát biển và Tìm kiếm cứu nạn (SAR).

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

                         MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG
┌───────────────────────────┬─────────────────────────────────────────────────────────┐
│ TRỤ CỘT TÁC ĐỘNG          │ KẾT QUẢ ĐỊNH LƯỢNG VÀ THỰC TIỄN                         │
├───────────────────────────┼─────────────────────────────────────────────────────────┤
│ 1. Khoa học & Học thuật   │ • Tiên phong giải thuật Euler ngược thời gian           │
│                           │ • Độ chính xác phục dựng nguồn thải đạt NASH > 0.88     │
├───────────────────────────┼─────────────────────────────────────────────────────────┤
│ 2. Công nghiệp & Khai thác│ • Kiểm soát an toàn sinh thái 100% bồn trũng dầu khí   │
│                           │ • Giảm thiểu thời gian khoanh vùng thu gom dầu tràn 30% │
├───────────────────────────┼─────────────────────────────────────────────────────────┤
│ 3. Quản lý Nhà nước       │ • Cung cấp công cụ điều tra xả thải trộm hàng hải       │
│                           │ • Cơ sở định lượng yêu cầu bồi thường theo UNCLOS 1982 │
├───────────────────────────┼─────────────────────────────────────────────────────────┤
│ 4. Kinh tế - Xã hội       │ • Bảo vệ nguồn lợi thủy sản, du lịch biển miền Trung    │
│                           │ • Tiết kiệm hàng chục tỷ đồng ngân sách dọn sạch bãi    │
└───────────────────────────┴─────────────────────────────────────────────────────────┘
  • Tác động học thuật: Luận án mở ra một hướng tiếp cận toán - tin mới trong ngành Hải dương học Việt Nam, thúc đẩy việc tự chủ công nghệ mô hình hóa môi trường biển và công bố các nghiên cứu chuyên sâu trên các tạp chí quốc tế chuyên ngành động lực học biển.
  • Chuyển đổi công nghiệp Dầu khí: Cung cấp cơ sở khoa học cho Tập đoàn Dầu khí Việt Nam (PVN) và các liên doanh khai thác dầu khí xây dựng kịch bản ứng phó sự cố rò rỉ đường ống đáy biển và sự cố giàn khoan tại các mỏ Bạch Hổ, Đại Hùng, Rồng, Lan Tây.
  • Bảo vệ cộng đồng và phát triển bền vững: Giảm thiểu thiệt hại kinh tế cho hàng triệu ngư dân nuôi trồng thủy hải sản ven biển và ngành du lịch biển tại các địa phương thường xuyên chịu ảnh hưởng của dầu tràn trôi dạt vô chủ (như Quảng Nam, Đà Nẵng, Bình Thuận, Bà Rịa - Vũng Tàu).

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận hoàn chỉnh về bài toán ngược Euler, mã nguồn và sơ đồ khối thuật toán sai phân hữu hạn đa tầng để mở rộng cho các bài toán lan truyền chất ô nhiễm khác (nhựa đại dương, phóng xạ, kim loại nặng, tảo nở hoa).
  • Cơ quan quản lý Nhà nước (Bộ Tài nguyên và Môi trường, Cục Biển và Hải đảo Việt Nam, Ủy ban Quốc gia Ứng phó sự cố, thiên tai và Tìm kiếm cứu nạn): Sở hữu công cụ khoa học chính xác để giám sát không gian biển, điều tra vi phạm hàng hải và ra quyết định điều động lực lượng ứng cứu kịp thời.
  • Các doanh nghiệp vận tải biển và Cảng biển quốc tế: Tối ưu hóa kế hoạch phòng ngừa sự cố tràn dầu, giảm thiểu rủi ro pháp lý và chi phí bảo hiểm hàng hải quốc tế.

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à thiết lập thành công khung toán học giải bài toán ngược thời gian bằng phương pháp Euler cho quá trình lan truyền và phong hóa dầu trên biển. Luận án đã mở rộng Phương trình Truyền - Khuếch tán liên tục (Advection-Diffusion Equation) và Lý thuyết Phong hóa Dầu mỏ của Mackay (1980) sang miền thời gian đảo nghịch ($t \rightarrow -t$). Bằng việc đảo chiều vector vận tốc tổng hợp liên hợp và thiết lập điều kiện nghiệm không âm kết hợp với sơ đồ sai phân khử khuếch tán số, luận án đã biến quá trình phân tán tiêu tán tự nhiên thành quá trình tích tụ hội tụ toán học, giải quyết triệt để hạn chế phân kỳ ngẫu hành vốn có của các mô hình hạt Lagrange BITT truyền thống.

2. Đổi mới phương pháp luận của luận án vượt trội hơn các nghiên cứu trước đây ở điểm nào?

So với mô hình BITT của Christensen et al. (2007) chỉ đạt tỷ lệ tìm nguồn danh định 0,1% và phương pháp thử sai trực quan của Vũ Thanh Ca et al. (2007), luận án vượt trội ở 3 khía cạnh:

  • Đồng nhất tuyệt đối không gian lưới tính với dữ liệu động lực đầu vào (CFSR, POM, SWAN), loại bỏ hoàn toàn sai số do phép nội suy Lagrange-Euler.
  • Tích hợp đồng thời 8 quá trình phong hóa hóa lý phụ thuộc nhiệt động lực học trong phương trình vi phân ngược, cho phép khôi phục chính xác thể tích và khối lượng riêng ban đầu của dầu tại thời điểm xả.
  • Áp dụng kỹ thuật sai phân lồng đa tầng dưới lưới bậc 2-8, duy trì tính ổn định số và bảo toàn khối lượng với sai số khoảng cách tâm vết dầu $< 5\text{ km}$ sau 7 ngày mô phỏng trên miền tính $2.500\text{ km}$.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và có bằng chứng định lượng đi kèm?

Phát hiện bất ngờ nhất là vai trò chi phối áp đảo của sóng Stokes và góc lệch gió tầng mặt so với dòng chảy đại dương thông thường. Trong điều kiện Biển Đông, dòng chảy địa phương tầng mặt chỉ đóng góp khoảng $30 - 45%$ vào tốc độ dịch chuyển của màng dầu; $55 - 70%$ động lực còn lại do ma sát gió và sóng Stokes dẫn dắt. Nếu không hiệu chỉnh ma trận góc lệch gió $\mathbf{D}_w$ và vận tốc Stokes $V_s$, mô hình xuôi sẽ lệch quỹ đạo thực tế tới $42%$ sau 48 giờ. Tuy nhiên, khi tích hợp đầy đủ công thức trôi dạt hiệu chỉnh, mô hình đạt độ tương quan $R = 0{,}97$ và chỉ số $\text{NASH} = 0{,}93$.

               ĐÓNG GÓP THÀNH PHẦN VẬN TỐC TRÔI DẠT DẦU TẦNG MẶT
   ┌────────────────────────────────────────┬───────────────────────────┐
   │ THÀNH PHẦN LỰC DẪN ĐỘNG                │ TỶ TRỌNG ĐÓNG GÓP (%)     │
   ├────────────────────────────────────────┼───────────────────────────┤
   │ Gió bề mặt 10m (Đã hiệu chỉnh góc lệch)│ ██████████████ 40 - 45%   │
   │ Sóng mặt (Dòng trôi sóng Stokes)       │ ████████ 20 - 25%         │
   │ Hoàn lưu dòng chảy đại dương tầng mặt  │ ███████████ 30 - 35%      │
   └────────────────────────────────────────┴───────────────────────────┘

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Có. Luận án cung cấp chi tiết toàn bộ quy trình tái lập bao gồm:

  • Toàn bộ hệ phương trình toán học giải tích, phương trình sai phân rời rạc hóa và các công thức tham số hóa phong hóa thực nghiệm (Phương trình 2.1 đến 2.30).
  • Hệ thống sơ đồ khối thuật toán chi tiết cho chương trình chính và các chương trình con tính toán động lực, phong hóa, sai phân bậc 1 và bậc 2.
  • Đặc tả thông số cấu hình miền tính lưới ($1 - 24^\circ\text{N}, 99 - 121^\circ\text{E}$), bước lưới $dx=dy=0{,}044^\circ$, bước thời gian thỏa mãn điều kiện ổn định Courant-Friedrichs-Lewy (CFL), và hướng dẫn trích xuất định dạng tệp đầu vào từ CFSR và POM.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year Research Agenda) được phác thảo như thế nào?

Chương trình nghiên cứu dài hạn bao gồm 4 lộ trình chiến lược:

  1. Năm 1–3: Nâng cấp mô hình lên cấu trúc 3 chiều hoàn chỉnh (3D Baroclinic Multiphase Model), tích hợp động lực học giọt dầu dưới sâu và tương tác trầm tích lơ lửng.
  2. Năm 4–5: Xây dựng hệ thống đồng hóa dữ liệu tự động với ảnh vệ tinh siêu phổ và Radar khẩu độ tổng hợp (SAR) đa nguồn.
  3. Năm 6–8: Ứng dụng kỹ thuật AI/Machine Learning để tối ưu hóa việc phân loại chữ ký hóa học của dầu kết hợp với mô hình số ngược thời gian nhằm xác định chính xác danh tính tàu xả thải từ dữ liệu AIS (Automatic Identification System).
  4. Năm 9–10: Xây dựng Hệ thống Chuyên gia Trợ giúp Quyết định Thời gian thực (Real-time Decision Support System) cấp quốc gia bảo vệ an ninh môi trường biển trên toàn bộ vùng biển Việt Nam và lân cận.

Kết luận

Công trình luận án tiến sĩ của NCS. Nguyễn Quốc Trinh là một bước tiến đột phá, đóng góp toàn diện cho ngành Hải dương học và công tác bảo vệ môi trường biển Việt Nam:

  1. Xây dựng thành công mô hình số trị Euler 2 chiều toàn diện: Giải quyết đồng thời bài toán dự báo lan truyền (xuôi thời gian) và tìm kiếm nguồn thải (ngược thời gian) cho vệt dầu tràn trên toàn bộ không gian Biển Đông.
  2. Đột phá lý thuyết bài toán ngược thời gian: Thiết lập khung toán học vững chắc cho bài toán ngược trong môi trường liên tục, chứng minh tính hội tụ và ổn định của phương pháp Euler đảo chiều vận tốc kết hợp bảo toàn tính chất hóa lý của khối dầu.
  3. Hoàn thiện hệ thống tham số hóa 8 quá trình phong hóa: Định lượng chính xác sự biến đổi phi tuyến tính của thể tích, mật độ, độ nhớt, sức căng bề mặt dưới điều kiện khí hậu nhiệt đới đặc thù của Biển Đông.
  4. Giải mã nguồn gốc các sự cố dầu tràn lịch sử: Phục dựng thành công vị trí phát thải của các đợt tràn dầu vô chủ nghiêm trọng năm 2007 và 2008 với độ chính xác cao ($\text{NASH} > 0{,}88$, $R > 0{,}96$), chứng minh tính ứng dụng thực tiễn vượt trội.
  5. Khẳng định tính tự chủ công nghệ mô hình hóa: Xây dựng hệ thống chương trình tính toán độc lập, tạo tiền đề vững chắc cho việc xây dựng hệ thống cảnh báo sớm và bảo vệ chủ quyền, an ninh môi trường biển đảo quốc gia trong kỷ nguyên kinh tế biển bền vững.