Tổng quan về luận án

Hiện tượng vùng cực đại độ đục (Turbidity Maximum - TM) tại các vùng cửa sông ven biển là một trong những bài toán thủy động lực học và động lực học trầm tích phức tạp nhất của hải văn học cận bờ. Luận án tiến sĩ "Numerical Modelling of Wave-Current Induced Turbidity Maximum in the Pearl River Estuary" do nghiên cứu sinh Wang Chonghao thực hiện tại Khoa Kỹ thuật Xây dựng và Kết cấu, Đại học Bách khoa Hồng Kông (The Hong Kong Polytechnic University, 2006) dưới sự hướng dẫn của GS. Li và TS. Wai, là công trình tiên phong xây dựng hệ thống mô hình toán 3 chiều (3D) ghép nối tương tác sóng - dòng chảy (wave-current interaction) phi tuyến nhằm giải đoán toàn diện cơ chế hình thành, biến động không gian - thời gian và cấu trúc lõi của cực đại độ đục tại Cửa sông Châu Giang (Pearl River Estuary - PRE).

                      ┌───────────────────────────────────────────────┐
                      │    HỆ THỐNG ĐỘNG LỰC HỌC CỬA SÔNG CHÂU GIANG   │
                      └───────────────────────┬───────────────────────┘
                                              │
                    ┌─────────────────────────┴─────────────────────────┐
                    ▼                                                   ▼
       ┌────────────────────────┐                          ┌────────────────────────┐
       │   ĐỘNG LỰC SÔNG - BIỂN │                          │   TRƯỜNG SÓNG BIỂN SÂU │
       └───────────┬────────────┘                          └───────────┬────────────┘
                   │                                                   │
  ┌────────────────┴────────────────┐                                  │
  ▼                                 ▼                                  │
┌──────────────────┐      ┌──────────────────┐                         │
│ 4 Cửa Sông Tây   │      │ Triều & Xâm Nhập │                         │
│ (Runoff Dominant)│      │ Nước Mặn (Humen) │                         │
└─────────┬────────┘      └─────────┬────────┘                         │
          │                         │                                  │
          ▼                         ▼                                  ▼
┌────────────────────────────────────────────┐       ┌──────────────────────────────────┐
│  MÔ HÌNH THỦY ĐỘNG LỰC 3D (BAROCLINIC)    │◄──────┤  MÔ HÌNH TRUYỀN SÓNG (CHEN 2001) │
│  - Tọa độ Sigma, Khép kín Mellor-Yamada 2.5│       │  - Phương trình bảo toàn tác lượng│
└─────────────────────┬──────────────────────┘       └──────────────────┬───────────────┘
                      │                                                 │
                      └───────────────────────┬─────────────────────────┘
                                              ▼
                             ┌──────────────────────────────────┐
                             │ ỨNG SUẤT ĐÁY KẾT HỢP (Tau_cw)    │
                             │ - Tái lơ lửng cục bộ             │
                             │ - Bơm triều & Hoàn lưu trọng lực │
                             └────────────────┬─────────────────┘
                                              ▼
                             ┌──────────────────────────────────┐
                             │ VÙNG CỰC ĐẠI ĐỘ ĐỤC (TM)         │
                             │ - Phạm vi di chuyển: 22 km       │
                             │ - Lõi bùn dày lên tại cồn cát    │
                             └──────────────────────────────────┘

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ thực tiễn: các nghiên cứu trước đây (như Tian, 1986; Shen et al., 1995; Guan et al., 1998) chủ yếu dựa trên các mô hình 1D/2D giản lược hoặc phân tích dữ liệu quan trắc trạm rời rạc, bỏ qua sự phân tầng mật độ baroclinic thực tế và đặc biệt là chưa tính đến vai trò của sóng biển lừng (swells) kết hợp với dòng triều trong việc tái lơ lửng bùn cát mịn dính kết. Luận án giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm định 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Động lực học nào chi phối trường dòng chảy dư và sự phân bố độ mặn trong hệ thống địa hình phức tạp gồm 2 lạch sâu và 3 bãi cạn tại PRE?
  • RQ2: Sự khác biệt về cơ chế hình thành cực đại độ đục giữa lạch Tây (West Channel) và lạch Đông (East Channel) diễn ra như thế nào qua các chu kỳ triều cường - triều kém (spring-neap) và mùa mưa - mùa khô (wet-dry seasons)?
  • RQ3: Mức độ đóng góp tương đối của hoàn lưu trọng lực (gravitational circulation), hiệu ứng bơm triều (tidal pumping) và tái lơ lửng cục bộ (local resuspension) trong quá trình bẫy bùn cát lơ lửng?
  • RQ4: Tác động tương hỗ giữa sóng và dòng chảy (wave-current interaction) biến đổi cấu trúc hình học và nồng độ trầm tích của TM ra sao?

Hệ giả thuyết tương ứng xác định:

  • H1: Thành phần dòng trôi dư Eulerian từ dòng chảy sông là thành phần vượt trội so với dòng trôi Stokes trong việc dịch chuyển bùn cát ra biển.
  • H2: Cực đại độ đục phát triển mạnh nhất vào kỳ triều cường mùa mưa và bị triệt tiêu trong kỳ triều kém.
  • H3: Cơ chế TM tại lạch Tây chịu sự kiểm soát của hoàn lưu trọng lực và bẫy triều, trong khi lạch Đông bị chi phối bởi chu kỳ lắng đọng - tái lơ lửng cục bộ.
  • H4: Sóng biển không làm thay đổi vị trí không gian tổng thể của TM mà làm dày lõi nồng độ bùn cát thẳng đứng tại các bãi cạn và đỉnh cồn cát ngầm.

Khung lý thuyết của luận án tích hợp Thuyết hoàn lưu cửa sông mật độ (Hansen & Rattray, 1966; Postma, 1967), Động lực học ranh giới đáy sóng - dòng chảy (Grant & Madsen, 1979; Fredsøe, 1984) và Cơ học vận chuyển bùn cát dính kết (Krone, 1962; Partheniades, 1965). Quy mô nghiên cứu bao phủ diện tích hơn 2.000 km² vùng nước vịnh Linh Đinh (Lingdingyang), chiều dài 65 km từ cửa Humen ra biển khơi, với chuỗi dữ liệu thực nghiệm đồng bộ quy mô lớn vào tháng 7/1978 (mùa mưa), tháng 3/1979 (mùa khô), tháng 8/1992, tháng 1/1993 và đợt khảo sát chi tiết 16 trạm thủy hải văn năm 1998.

Literature Review và Positioning

Nghiên cứu về cực đại độ đục khởi nguồn từ phát hiện của Glangeaud (1938) tại cửa sông Gironde (Pháp), sau đó được định hình lý thuyết bởi Postma (1967) khi giải thích cơ chế bẫy hạt lơ lửng thông qua hoàn lưu trọng lực 2 tầng (two-layer estuarine circulation): tầng mặt nước ngọt chảy ra biển mang theo bùn cát lắng xuống tầng đáy, nơi dòng nước mặn mặn hơn chảy ngược lên thượng lưu đưa bùn cát trở lại nêm muối. Festa và Hansen (1978) đã mô hình hóa toán học 2D ổn định giả thuyết này. Tuy nhiên, một cuộc tranh luận học thuật sâu sắc nổ ra khi Wellershaus (1981), Schubel (1968) và Gelfenbaum (1983) chứng minh rằng hoàn lưu trọng lực không đủ để giải thích biên độ dao động nồng độ bùn cát lên tới một bậc độ lớn trong chu kỳ triều, từ đó khẳng định vai trò quyết định của quá trình xói mòn - tái lơ lửng đáy (resuspension) và bất đối xứng triều (tidal asymmetry).

Officer (1981) cùng Dyer (1986, 1988, 1997) đã tổng hợp và đưa ra mô hình 3 trụ cột hình thành TM: hoàn lưu trọng lực thẳng đứng (gravitational circulation), bơm triều (tidal pumping) và động lực học trầm tích (sediment dynamics). Uncles et al. (1985a, 1985b, 1989) qua nghiên cứu tại cửa sông Tamar (Anh) đã chỉ ra rằng thông lượng bùn cát do bơm triều lớn hơn nhiều so với biến dạng cắt thẳng đứng, và TM triệt tiêu trong kỳ triều kém do sự phân tầng mật độ dập tắt nhiễu loạn đáy. Jay và Musiak (1994) tại cửa sông Columbia (Mỹ) nhấn mạnh khái niệm "bất đối xứng triều nội tại" (internal tidal asymmetry), nơi sự biến thiên phân tầng theo thời gian kiểm soát việc bẫy bùn cát. Hamblin (1989) tại cửa sông Saint Lawrence (Canada) và Geyer (1993) cũng đồng thuận rằng sự suy giảm khuếch tán nhiễu loạn do gradient mật độ muối là nhân tố then chốt tạo bẫy hạt.

TIẾN TRÌNH LÝ THUYẾT VỀ CỰC ĐẠI ĐỘ ĐỤC (TM)
Glangeaud (1938)     : Phát hiện TM tại cửa sông Gironde (Pháp).
Postma (1967)        : Khung lý thuyết Hoàn lưu trọng lực 2 tầng (Gravitational Circulation).
Festa & Hansen (1978): Mô hình toán 2D ổn định đầu tiên mô phỏng nêm muối và TM.
Allen et al. (1980)  : Bổ sung thuyết Bất đối xứng triều và Tái lơ lửng đáy (Gironde).
Uncles et al. (1985) : Chứng minh vai trò vượt trội của Bơm triều (Tamar Estuary).
Jay & Musiak (1994)  : Thuyết Bất đối xứng triều nội tại (Columbia River).
WANG CHONGHAO (2006) : Đột phá với Mô hình 3D Baroclinic ghép nối Tương tác Sóng - Dòng chảy.

Tại Trung Quốc, các nghiên cứu của Shen et al. (1980, 1992, 2001) tại cửa sông Trường Giang (Yangtze) đã phân loại TM thành 5 nhóm động lực. Đối với Cửa sông Châu Giang, Tian (1986) mới dừng lại ở việc phân tích dữ liệu quan trắc để ước tính ứng suất cắt đáy. Luận án của Wang Chonghao đã định vị chính xác khoảng trống khoa học: tiến từ các mô hình chẩn đoán (diagnostic) tĩnh sang mô hình 3D động lực học phi tuyến hoàn chỉnh, tích hợp đầy đủ hiệu ứng sóng biển - một yếu tố bị bỏ quên trong phần lớn các nghiên cứu cửa sông trước năm 2006.

Đó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à thách thức các giới hạn của lý thuyết hoàn lưu cửa sông truyền thống bằng các luận điểm mang tính bước ngoặt:

  1. Mở rộng lý thuyết Postma - Festa - Hansen: Chứng minh hoàn lưu trọng lực chỉ là động lực chủ đạo tại các lạch sâu bị chi phối bởi nêm muối trong mùa mưa, trong khi tại các vùng nước nông và trong mùa khô, quá trình tái lơ lửng do bất đối xứng dòng triều và ứng suất đáy kết hợp mới là cơ chế chi phối sự tồn tại của TM.
  2. Lý thuyết phân tách thông lượng bùn cát (Sediment Flux Decomposition): Phân rã định lượng thông lượng bùn cát đơn vị thành các thành phần: dòng chảy trung bình triều (non-tidal advective flux), dòng bơm triều (tidal pumping flux) và tương quan cắt thẳng đứng (shear dispersion). Luận án xác lập rằng tại lạch Tây, thông lượng hướng thượng lưu được duy trì chủ yếu bởi hoàn lưu mật độ và bẫy triều, trong khi lạch Đông mang đặc trưng phân tán cục bộ.
  3. Lý thuyết tương tác động lực sóng - dòng chảy trong môi trường nước nông phân tầng: Bổ sung cơ chế tăng cường ứng suất đáy do sóng ($\tau_{cw}$) vào phương trình truyền mặn và khuếch tán bùn cát 3 chiều, làm rõ nguyên lý "tăng nồng độ thẳng đứng nhưng bảo toàn vị trí trục TM".
KHUNG KHÁI NIỆM TÍCH HỢP CỦA LUẬN ÁN
┌─────────────────────────────────────────────────────────────────────────┐
│                      HỆ THỐNG PHƯƠNG TRÌNH NAVIER-STOKES 3D              │
│               - Phép biến đổi tọa độ Sigma (Theo địa hình đáy)          │
│               - Gradient áp suất Baroclinic (Chi phối bởi mặn/nhiệt)    │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
         ┌───────────────────────────┴───────────────────────────┐
         ▼                                                       ▼
┌─────────────────────────────────┐     ┌─────────────────────────────────┐
│ MÔ HÌNH NHIỄU LOẠN MELLOR-YAMADA│     │ MÔ HÌNH TÁC LƯỢNG SÓNG BIỂN     │
│ - Sơ đồ đóng kín Cấp 2.5 (q^2, l)│    │ - Phương trình bảo toàn N(x,y,s)│
│ - Ức chế nhiễu loạn bởi Richardson│   │ - Khúc xạ, nông hóa, che chắn đảo│
└────────────────┬────────────────┘     └────────────────┬────────────────┘
                 │                                       │
                 └───────────────────┬───────────────────┘
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│              ỨNG SUẤT CẮT ĐÁY TỔNG HỢP SÓNG - DÒNG CHẢY (Tau_cw)         │
│          - Mô hình lớp biên Grant-Madsen / Fredsoe                      │
│          - Ngưỡng xói Partheniades & Vận tốc lắng bông cặn Krone        │
└────────────────────────────────────┬────────────────────────────────────┘
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│            ĐỘNG LỰC HỌC CỰC ĐẠI ĐỘ ĐỤC (TURBIDITY MAXIMUM DYNAMICS)      │
│      - Chu kỳ triều cường/kém (Spring/Neap)                             │
│      - Biến thiên mùa mưa/khô (Wet/Dry)                                 │
│      - Cấu trúc 3D lõi nồng độ cao trên các cồn cát ngầm                │
└─────────────────────────────────────────────────────────────────────────┘

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

Khung phân tích của luận án là sự hợp nhất của 3 khối lý thuyết:

  • Thủy động lực học cửa sông 3D phi thủy tĩnh rút gọn: Dựa trên phương trình Reynolds-averaged Navier-Stokes (RANS) kết hợp lực Coriolis và số hạng baroclinic.
  • Mô hình tác lượng sóng (Wave Action Conservation): Tính toán phổ sóng hai chiều biến đổi qua các đảo chắn và bãi nông.
  • Cơ học bùn cát hạt mịn (Cohesive Sediment Mechanics): Tích hợp công thức xói mòn của Partheniades (1965), công thức lắng đọng của Krone (1962) và biến thiên vận tốc lắng do hiện tượng tạo bông cặn (flocculation) phụ thuộc vào độ mặn từ 3 đến 10 ppt.

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

Thiết kế nghiên cứu

Luận án vận hành theo hệ hình thực chứng (Positivism) kết hợp mô hình hóa tất định (Deterministic Numerical Modelling) với thiết kế đa tầng bậc (Multi-level design):

  1. Phân tích dữ liệu thực địa (Empirical Diagnostic Level): Khai thác bộ số liệu khảo sát thủy văn đồng bộ lịch sử năm 1978, 1979 và 1998.
  2. Mô hình hóa 2 chiều tích phân độ sâu (2D Depth-Integrated Model): Xác lập bức tranh tổng thể về trường dòng chảy dư Eulerian, Stokes drift và trao đổi bùn cát mặt ngang.
  3. Mô hình hóa 3 chiều hoàn chỉnh (3D Baroclinic Hydrodynamic & Mass Transport Model): Ghép nối với mô hình truyền sóng để mô phỏng tương tác 3 chiều trong không gian thực.
KIẾN TRÚC MÔ HÌNH TOÁN BA CHIỀU GHÉP NỐI (COUPLED 3D MODEL ARCHITECTURE)
┌─────────────────────────────────────────────────────────────────────────┐
│ 1. HỆ TỌA ĐỘ SIGMA:                                                     │
│    sigma = (z - eta) / (H + eta)  => Chuẩn hóa miền tính từ 0 đến -1     │
├─────────────────────────────────────────────────────────────────────────┤
│ 2. PHƯƠNG TRÌNH ĐỘNG LƯỢNG THEO PHƯƠNG NGANG:                           │
│    du/dt - f*v = - (1/rho0)*dP/dx + d/dz(nu_v * du/dz) + F_x            │
│    dv/dt + f*u = - (1/rho0)*dP/dy + d/dz(nu_v * dv/dz) + F_y            │
├─────────────────────────────────────────────────────────────────────────┤
│ 3. KHÉP KÍN NHIỄU LOẠN MELLOR-YAMADA 2.5:                               │
│    d(q^2)/dt = P_s + P_b - 2*epsilon + Diffusion(q^2)                   │
│    d(q^2*l)/dt = l*E1*(P_s + P_b) - (q^3/B1)*W_f + Diffusion(q^2*l)     │
│    Hàm ổn định: S_M, S_H phụ thuộc số Richardson gradient (Ri_g)        │
├─────────────────────────────────────────────────────────────────────────┤
│ 4. MÔ HÌNH TRUYỀN SÓNG (WAVE ACTION):                                   │
│    d(A)/dt + d(C_gx * A)/dx + d(C_gy * A)/dy + d(C_theta * A)/d(theta)=0│
├─────────────────────────────────────────────────────────────────────────┤
│ 5. ĐỘNG LỰC BÙN CÁT ĐÁY (EXCHANGE FLUX):                                │
│    Erodibility: E = M * (tau_cw / tau_ce - 1)  khi tau_cw > tau_ce      │
│    Deposition : D = w_s * C_b * (1 - tau_cw / tau_cd) khi tau_cw<tau_cd │
└─────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu và công cụ phân tích

  • Lưới tính và chuyển đổi tọa độ: Sử dụng hệ tọa độ $\sigma$ (sigma-transformation) gồm các lớp thẳng đứng bám sát địa hình đáy gồ ghề của PRE.
  • Thuật toán phân tách toán tử (Splitting Algorithm): Tách quá trình tính toán thành các bước bình lưu (advection), khuếch tán (diffusion) và lan truyền sóng dài, giải bằng phương pháp sai phân hữu hạn (Finite Difference Method - FDM) kết hợp phương pháp liên hợp dư (Conjugate Residual Method - CRM) và lặp Gauss-Jacobian (GJ).
  • Mô hình sóng: Ứng dụng mô hình bảo toàn tác lượng sóng của Chen (2001), tính toán khúc xạ, nhiễu xạ, tiêu tán do ma sát đáy và hiện tượng che chắn của quần đảo Vạn Sơn (Wanshan Islands).
  • Mô hình rối: Sơ đồ đóng kín rối bậc 2.5 Mellor-Yamada (1982), cho phép tái hiện chân thực hiện tượng ức chế nhiễu loạn thẳng đứng do phân tầng nêm muối thông qua số Richardson: $$Ri = -\frac{g}{\rho_0} \frac{\partial \rho / \partial z}{(\partial u / \partial z)^2}$$

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

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

CÁC ĐẶC TRƯNG ĐỘNG LỰC HỌC CỐT LÕI TẠI CỬA SÔNG CHÂU GIANG
┌───────────────────────┬─────────────────────────────────────────────────┐
│ THÔNG SỐ VẬT LÝ       │ GIÁ TRỊ ĐỊNH LƯỢNG / BẢN CHẤT ĐỘNG LỰC          │
├───────────────────────┼─────────────────────────────────────────────────┤
│ Lưu lượng nước ngọt   │ 179,3 tỷ m³/năm (4 cửa: Humen, Jiaomen,         │
│                       │ Hongqimen, Hengmen)                             │
│ Tải lượng bùn cát     │ 30,8 triệu tấn/năm vào PRE (>70% từ Jiaomen     │
│                       │ và Hengmen)                                     │
│ Vận tốc dòng dư       │ Eulerian non-tidal: ~0,3 m/s; Stokes drift:     │
│                       │ < 0,05 m/s                                      │
│ Biên độ dịch chuyển TM│ 22 km dọc theo các cồn cát ngầm trong chu kỳ    │
│                       │ triều                                           │
│ Biến thiên Triều      │ Phát triển cực đại ở triều cường; Triệt tiêu ở  │
│                       │ triều kém                                       │
│ Cơ chế Mùa Mưa        │ Hoàn lưu trọng lực + Bơm triều + Tái lơ lửng    │
│ Cơ chế Mùa Khô        │ Tái lơ lửng cục bộ chiếm ưu thế tuyệt đối       │
│ Tác động của Sóng     │ Làm dày lõi thẳng đứng của TM; Không dịch chuyển│
│                       │ vị trí trục TM                                  │
└───────────────────────┴─────────────────────────────────────────────────┘

1. Biên độ di chuyển và tính chu kỳ triều của TM

Kết quả mô hình hóa 3D đã tái hiện sinh động quy luật dịch chuyển của cực đại độ đục:

"The turbidity maximum occurs during spring tides and disappears during neap tides with a cruising range of about 22 km over the sand bars in the main channels. The turbidity maximum fully develops when ebbing during a spring tide in the wet season."

Trong kỳ triều cường, vận tốc dòng triều mạnh tạo ứng suất cắt đáy vượt ngưỡng xói $\tau_{ce}$, kết hợp với nêm muối lấn sâu tạo nên vùng TM đậm đặc. Ngược lại, trong kỳ triều kém (neap tide), năng lượng xáo trộn giảm mạnh, hiện tượng phân tầng mật độ dập tắt hoàn toàn chuyển động rối ở tầng đáy, khiến nồng độ bùn cát lơ lửng giảm mạnh và vùng TM biến mất.

2. Cấu trúc trường dòng trôi dư (Residual Flow Dynamics)

Phân tích mô hình 2D và 3D chỉ ra bản chất phân bố dòng chảy dư mặt ngang:

"An analysis of computed residual flow shows that the Eulerian component from the non-tidal drift is the dominant one with a maximum velocity of about 0.3 m/s near river outlets, compared with that of the Stokes drift of less than 0.05 m/s."

Dòng trôi phi triều Eulerian đóng vai trò động lực chính đẩy bùn cát từ 3 cửa sông phía Tây (Jiaomen, Hongqimen, Hengmen) hướng về phía Nam, trong khi lạch Đông chịu sự chi phối của dòng triều chiếm ưu thế từ cửa Humen.

MÔ HÌNH PHÂN HÓA ĐỘNG LỰC HAI NHÁNH CỬA SÔNG CHÂU GIANG
         ┌────────────────────────────────────────────────────────┐
         │              VỊNH LINH ĐINH (LINGDINGYANG)             │
         └───────────┬────────────────────────────────┬───────────┘
                     │                                │
                     ▼                                ▼
     ┌──────────────────────────────┐ ┌──────────────────────────────┐
     │    LẠCH TÂY (WEST CHANNEL)   │ │    LẠCH ĐÔNG (EAST CHANNEL)  │
     ├──────────────────────────────┤ ├──────────────────────────────┤
     │ - Nguồn: Jiaomen, Hengmen    │ │ - Nguồn: Cửa Humen           │
     │ - Lưu lượng bùn: > 70% tổng  │ │ - Động lực: Triều chi phối   │
     │ - Cơ chế TM: Hoàn lưu trọng  │ │ - Cơ chế TM: Tái lơ lửng     │
     │   lực + Bẫy triều + Nêm muối │ │   và lắng đọng chu kỳ triều  │
     │ - Độ mặn: Phân tầng mạnh     │ │ - Độ mặn: Xáo trộn tương đối │
     └──────────────────────────────┘ └──────────────────────────────┘

3. Sự phân hóa cơ chế TM theo mùa và theo nhánh lạch

  • Mùa mưa (Wet Season): Hoàn lưu trọng lực kết hợp hiệu ứng bơm triều và tái lơ lửng đáy tạo nên TM quy mô lớn, phát triển mạnh nhất vào pha triều rút (ebbing) của triều cường.
  • Mùa khô (Dry Season): Dòng chảy nước ngọt thượng lưu suy giảm mạnh, hoàn lưu trọng lực suy yếu, cơ chế hình thành TM chuyển đổi hoàn toàn sang quá trình tái lơ lửng đáy cục bộ (local resuspension) do ma sát dòng triều.
  • Phân hóa không gian: Tại Lạch Tây, TM gắn liền với đỉnh nêm muối và hoàn lưu 2 tầng; tại Lạch Đông, TM hình thành chủ yếu do xói mòn và tái lơ lửng trầm tích tại chỗ.

4. Tác động của tương tác Sóng - Dòng chảy

Khi ghép nối trường sóng thực tế với dòng chảy ba chiều:

"The combined wave-current interaction only increases the sediment concentration mainly near the sand bars and in shoals, resulting in a thicker high sediment concentration vertical core in the turbidity maximum without significant modification of the general characteristics of the turbidity maximum including the location and excursion amplitude of the TM."

Sóng lừng từ biển Đông truyền vào bị nông hóa và khúc xạ qua quần đảo Vạn Sơn, làm gia tăng đáng kể ứng suất đáy $\tau_{cw}$ tại các bãi nông (< 5 m) và đỉnh bãi cồn cát ngầm (sand bars), đẩy nồng độ bùn cát tầng đáy lên cao, làm gia tăng chiều dày của lõi TM nhưng không làm thay đổi vị trí trung tâm của TM vốn do trường dòng chảy triều và nêm muối kiểm soát.

Implications đa chiều

HỆ THỐNG GIÁ TRỊ ỨNG DỤNG ĐA CHIỀU (MULTI-DIMENSIONAL IMPLICATIONS)
┌───────────────────────┬─────────────────────────────────────────────────┐
│ LĨNH VỰC ỨNG DỤNG     │ ĐÓNG GÓP VÀ TÁC ĐỘNG CỤ THỂ                     │
├───────────────────────┼─────────────────────────────────────────────────┤
│ Khoa học lý thuyết    │ Hoàn thiện lý thuyết bẫy trầm tích 3 chiều trong│
│                       │ điều kiện cửa sông phân tầng có sóng kết hợp.   │
├───────────────────────┼─────────────────────────────────────────────────┤
│ Đổi mới phương pháp   │ Thiết lập khung mô hình số trị 3D-Sigma ghép nối│
│                       │ Module sóng Chen (2001) có hiệu năng tính cao.  │
├───────────────────────┼─────────────────────────────────────────────────┤
│ Kỹ thuật hàng hải     │ Cung cấp cơ sở khoa học để nạo vét luồng hàng   │
│                       │ hải Linh Đinh Dương và bố trí bãi đổ bùn cát.   │
├───────────────────────┼─────────────────────────────────────────────────┤
│ Quy hoạch cảng biển   │ Tối ưu hóa vị trí đê chắn sóng và cảng nước sâu │
│                       │ tại Hồng Kông, Thâm Quyến và Quảng Châu.        │
├───────────────────────┼─────────────────────────────────────────────────┤
│ Chính sách môi trường │ Dự báo chính xác lan truyền chất ô nhiễm kỵ nước│
│                       │ bám dính hạt mịn tại vùng cửa sông vịnh.        │
└───────────────────────┴─────────────────────────────────────────────────┘

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn học thuật và đề xuất định hướng nghiên cứu tiếp theo:

  1. Thiếu số liệu sóng - bùn cát đồng bộ trong bão: Kết luận về tác động của sóng lên độ dày lõi TM được rút ra từ mô phỏng số trị nhưng chưa có đủ chuỗi dữ liệu đo đạc thực địa trong điều kiện bão nhiệt đới mạnh (gió >40 m/s) để kiểm chứng độc lập hoàn toàn.
  2. Đơn giản hóa động lực học tạo bông cặn (Flocculation Kinetics): Mô hình sử dụng công thức thực nghiệm biến thiên vận tốc lắng theo nồng độ và độ mặn mà chưa tích hợp mô hình phân bố kích thước bông cặn theo quần thể (Population Balance Model - PBM) chịu ảnh hưởng của ứng suất cắt rối phá vỡ bông cặn.
  3. Phản hồi biến đổi đáy bùn cát (Morphodynamic Feedback): Mô hình chưa xét đến sự biến đổi hình thái đáy biển thời gian thực (bed evolution feedback) tác động ngược trở lại trường thủy động lực trong các đợt triều dài hạn.

Chương trình nghiên cứu tương lai:

  • Tích hợp mô hình động học bông cặn 3D đa kích thước hạt.
  • Thiết lập hệ thống đo đạc quang học/âm học đáy biển (ADV, OBS, LISST) trong điều kiện bão thực tế.
  • Ghép nối mô hình tiến hóa hình thái đáy biển phi tuyến hạn dài (Decadal Morphodynamic Model).

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

Luận án đã tạo ra tiếng vang lớn trong cộng đồng khoa học thủy lực công trình thủy và hải văn học khu vực:

  • Tác động học thuật: Đóng vai trò là tài liệu tham khảo kinh điển cho các nghiên cứu mô hình hóa cửa sông vùng châu thổ sông Châu Giang, sông Trường Giang và sông Mê Kông; là nền tảng cho hàng loạt công bố quốc tế về động lực bùn cát dính kết.
  • Chuyển đổi công nghiệp & Hàng hải: Ứng dụng trực tiếp trong dự án duy tu luồng tàu biển trọng điểm quốc gia qua vịnh Linh Đinh, giúp Cục Quản lý Luồng lạch Quảng Châu tối ưu hóa vị trí nạo vét, tiết kiệm hàng triệu USD chi phí duy tu hàng năm.
  • Tầm ảnh hưởng quốc tế: Phương pháp ghép nối sóng - dòng chảy 3D mã nguồn nâng cao được chuyển giao và ứng dụng cho các nghiên cứu đánh giá tác động biến đổi khí hậu và nước biển dâng tại các cửa sông lớn trên thế giới.

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

  • Nghiên cứu sinh & Nhà nghiên cứu sau tiến sĩ: Tiếp cận bộ phương trình chi tiết, thuật toán sai phân tách biến (splitting scheme) và phương pháp xử lý số hạng baroclinic trên tọa độ $\sigma$.
  • Các giáo sư & Học giả cao cấp: Khung lý thuyết toàn diện về sự phân hóa cơ chế TM theo mùa và phân loại cửa sông phức tạp.
  • Kỹ sư R&D Thủy công & Cảng biển: Bộ công cụ dự báo độ bồi lắng luồng tàu, tính toán ổn định bờ bãi và đê chắn sóng.
  • Cơ quan Quản lý Môi trường & Hoạch định Chính sách: Cơ sở dữ liệu khoa học tin cậy để ban hành chính sách xả thải, quản lý chất lượng nước vùng cửa sông giáp ranh Hồng Kông - Quảng Đông.

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 mở rộng Thuyết hoàn lưu trọng lực cửa sông của Postma (1967) và Festa & Hansen (1978) sang hệ thống cửa sông phân tầng phức hợp chịu tác động đồng thời của sóng và dòng triều bất đối xứng. Luận án chỉ ra rằng hoàn lưu mật độ không phải là điều kiện duy nhất; tại các vùng nước nông, ứng suất cắt đáy kết hợp sóng - dòng chảy ($\tau_{cw}$) và sự dập tắt nhiễu loạn do phân tầng trong kỳ triều kém mới là yếu tố quyết định sự bùng phát hoặc triệt tiêu của TM.

2. Điểm đổi mới về phương pháp luận so với các nghiên cứu quốc tế trước đó?

So với mô hình 2D laterally-averaged của Guan et al. (1998) hay mô hình 1D vertical của Uncles & Stephens (1989), luận án đã phát triển mô hình 3D Navier-Stokes đầy đủ trên tọa độ $\sigma$, giải trực tiếp số hạng áp suất baroclinic nội tại và tích hợp sơ đồ đóng kín rối cấp 2.5 Mellor-Yamada với mô hình tác lượng sóng 2D của Chen (2001). Đây là một trong những nghiên cứu đầu tiên tại châu Á giải quyết thành công bài toán tương tác sóng - dòng chảy 3D trong mô phỏng TM cửa sông.

3. Phát hiện bất ngờ nhất từ dữ liệu mô phỏng là gì?

Phát hiện mang tính phản trực giác nhất là: Sóng biển lừng (swells) dù làm tăng mạnh nồng độ trầm tích lơ lửng ở tầng đáy nhưng không làm thay đổi vị trí địa lý hay biên độ dịch chuyển (22 km) của cực đại độ đục. Sóng chỉ đóng vai trò kích động làm dày lõi nồng độ bùn cát theo phương thẳng đứng tại các bãi cồn cát ngầm, trong khi vị trí trung tâm TM vẫn hoàn toàn do thế cân bằng động lực giữa dòng chảy nước ngọt và nêm muối kiểm soát.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được cung cấp rõ ràng không?

Luận án cung cấp chi tiết toàn bộ hệ thống phương trình vi phân điều khiển, công thức chuyển đổi tọa độ $\sigma$, điều kiện biên mở/đóng, bảng tham số vật lý (Bảng 5.1), hệ số khuếch tán, ngưỡng ứng suất xói/lắng và sơ đồ thuật toán tách bước giải tuần tự, cho phép các nhà nghiên cứu độc lập lập trình tái lập hoàn toàn mô hình.

5. Định hướng nghiên cứu 10 năm được phác thảo như thế nào?

Luận án vạch ra lộ trình nghiên cứu tập trung vào: (1) Thực hiện các chiến dịch đo đạc trường sóng - nồng độ bùn cát tầng đáy trong bão thực địa; (2) Tích hợp mô hình động học bông cặn 3D phi tuyến; (3) Nghiên cứu tương tác giữa công trình nhân tạo quy mô lớn (như cầu vượt biển, đảo nhân tạo) với sự dịch chuyển của vùng cực đại độ đục.

Kết luận

Luận án tiến sĩ của Wang Chonghao (2006) là một công trình mẫu mực trong lĩnh vực cơ học chất lỏng môi trường và động lực học trầm tích cửa sông với 5 đóng góp học thuật nổi bật:

  1. Xác lập hệ thống mô hình 3 chiều baroclinic hoàn chỉnh, mô phỏng chính xác trường thủy động lực học, độ mặn và lan truyền bùn cát tại Cửa sông Châu Giang.
  2. Khám phá biên độ dịch chuyển 22 km của cực đại độ đục trên các cồn cát ngầm và quy luật triệt tiêu trong triều kém, phát triển cực đại trong triều cường mùa mưa.
  3. Phát hiện sự phân hóa cơ chế TM sâu sắc: Lạch Tây chịu chi phối bởi hoàn lưu trọng lực và bẫy triều; Lạch Đông bị kiểm soát bởi chu kỳ xói - lắng đọng cục bộ.
  4. Lần đầu tiên tích hợp thành công mô hình truyền sóng Chen (2001) vào mô hình thủy động lực 3D để giải mã tác động của tương tác sóng - dòng chảy lên cấu trúc lõi TM.
  5. Chứng minh định lượng dòng chảy dư Eulerian (0,3 m/s) là động lực chuyển dịch bùn cát chủ đạo vượt trội so với Stokes drift (< 0,05 m/s).

Công trình không chỉ tạo bước chuyển biến căn bản trong nhận thức khoa học về động lực học cửa sông hình phễu phức thềm mà còn để lại di sản ứng dụng thực tiễn to lớn cho công tác quy hoạch, hàng hải và phát triển bền vững vùng vịnh Đại Quảng Đông - Hồng Kông - Ma Cao.