Tổng quan nghiên cứu

Hơn 70% các hệ thống đầm phá ven biển trên toàn cầu chịu sự chi phối mạnh mẽ của chế độ động lực học thủy triều và quá trình bồi lắng trầm tích phức tạp. Trong số đó, các hệ thống đầm phá có từ hai cửa biển trở lên thường xuyên đối mặt với nguy cơ mất ổn định hình thái, dẫn đến việc bồi lấp tự nhiên hoặc xói lở dịch chuyển luồng lạch với tỷ lệ biến động lòng dẫn lên tới 60% đến 80% sau các chu kỳ thiên tai cực đoan. Vấn đề cốt lõi trong kỹ thuật chỉnh trị bờ biển là xác định liệu một hệ thống nhiều cửa có thể cùng tồn tại ở trạng thái ổn định lâu dài hay không.

Trước đây, lý thuyết kinh điển của Escoffier (1940) và các phát triển của Van de Kreeke (1990) dựa trên mô hình tham số tập trung tuyến tính hóa từng kết luận rằng hệ thống hai cửa biển nối chung một vịnh nước dao động đồng nhất không thể đạt được trạng thái cân bằng mặt cắt ổn định vô điều kiện. Nhằm giải quyết điểm nghẽn học thuật này, luận văn thạc sĩ của tác giả C. Borsje thực hiện tại Đại học Công nghệ Delft (TU Delft) tập trung nghiên cứu sự phụ thuộc của các đường cong cân bằng dòng chảy vào đặc tính hình học cửa dẫn, đầm phá và chế độ thủy triều đại dương. Mục tiêu cụ thể là chứng minh sự tồn tại của trạng thái ổn định vô điều kiện (dạng Type E) thông qua các kịch bản vi sai biên độ triều ngoài biển và phân vùng nội vịnh (mô hình vách ngăn).

Phạm vi nghiên cứu bao gồm mô hình lý thuyết vịnh tham chiếu diện tích 100 km² với chu kỳ bán nhật triều 44.700 giây, kết hợp kiểm chứng thực nghiệm tại vịnh Matagorda (Mỹ) với diện tích 317 km² và ứng dụng thực tiễn tại đầm phá Tam Giang - Cầu Hai (Thừa Thiên Huế, Việt Nam). Nghiên cứu cung cấp cơ sở định lượng vững chắc với lưới tính toán 2.500 điểm ma trận thủy lực, hỗ trợ giảm thiểu chi phí nạo vét luồng hàng hải từ 25% đến 30% và nâng cao độ an toàn cho các công trình chỉnh trị cửa biển.

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 vận dụng lý thuyết ổn định mặt cắt cửa triều của Escoffier (1940), kết hợp với mối quan hệ kinh nghiệm giữa diện tích mặt cắt họng cửa $A$ và lăng đạo triều $P$ theo O'Brien ($A = C \cdot P^n$). Trạng thái cân bằng động lực học đạt được khi vận tốc dòng triều cực đại $\hat{u}$ bằng chính vận tốc cân bằng $u_{eq} = 1.0\text{ m/s}$, tương ứng với ứng suất tiếp đáy cân bằng $\tau_{eq}$ dao động trong khoảng hẹp từ $3.5\text{ N/m}^2$ đến $5.5\text{ N/m}^2$. Khi $\hat{u} > u_{eq}$, lòng dẫn ở trạng thái xói sâu; ngược lại khi $\hat{u} < u_{eq}$, lòng dẫn rơi vào trạng thái bồi cạn.

Luận văn thiết lập hệ thống phân loại 5 dạng cấu hình đường cong dòng chảy cân bằng trong mặt phẳng tọa độ $(A_1, A_2)$:

  1. Dạng Type A: Không tồn tại giao điểm ổn định, một hoặc cả hai cửa sẽ bị bồi lấp hoàn toàn.
  2. Dạng Type B: Có hai giao điểm nhưng đều nằm trên nhánh mất ổn định.
  3. Dạng Type C: Hai giao điểm, trong đó một giao điểm nằm trên nhánh ổn định của một cửa, dẫn đến việc chỉ duy nhất một cửa tồn tại.
  4. Dạng Type D: Xuất hiện giao điểm nằm trên nhánh ổn định của cả hai cửa nhưng mang tính ổn định có điều kiện; sau nhiễu động lớn, hệ thống có xu hướng chỉ giữ lại một cửa.
  5. Dạng Type E: Trạng thái ổn định vô điều kiện, trong đó các nhánh ổn định giao nhau tạo thành miền hội tụ khép kín, đảm bảo cả hai cửa cùng mở bền vững.

Hệ thống thủy động lực học được mô tả thông qua hệ phương trình vi phân phi tuyến gồm phương trình liên tục bảo toàn thể tích nước trong vịnh và phương trình động lượng 1D cho từng kênh dẫn, có xét đến hệ số ma sát đáy $f = 0.004$, hệ số tổn thất cục bộ cửa vào và cửa ra $m = 1.0$, cùng quán tính dòng chảy cục bộ.

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

Nguồn dữ liệu của nghiên cứu bao gồm tập thông số hình thái học thực địa từ vịnh Matagorda (Mỹ) với kênh Pass Cavallo dài 2.000 m, kênh Matagorda dài 2.880 m và hệ thống đầm phá Tam Giang - Cầu Hai với chuỗi số liệu đo đạc mặt cắt qua các cửa Thuận An, Tư Hiền và vết mở Hòa Duẩn.

Về cỡ mẫu và phương pháp chọn mẫu tính toán, tác giả thiết lập lưới quét tham số toàn diện (grid parameter sweep) ma trận $50 \times 50$, tạo ra 2.500 tổ hợp diện tích mặt cắt $(A_1, A_2)$ cho mỗi kịch bản mô phỏng trong miền giá trị từ $0\text{ m}^2$ đến $14.000\text{ m}^2$. Phương pháp chọn mẫu lưới đều này giúp bao quát toàn bộ không gian trạng thái thủy lực mà không bỏ sót các điểm kỳ dị chuyển tiếp.

Phương pháp phân tích được lựa chọn là giải tích kết hợp lặp số học (analytical-numerical iteration) để giải hệ từ 3 đến 5 phương trình vi phân phi tuyến. Phương pháp này được ưu tiên thay vì mô phỏng 3D phức tạp nhằm tối ưu hóa tốc độ xử lý hàng nghìn kịch bản tham số nhưng vẫn giữ được bản chất vật lý của tương tác lăng đạo triều. Thuật toán lặp tự điều chỉnh lưu lượng $Q_1, Q_2, Q_3$ và dao động mực nước $h_1, h_2$ cho đến khi độ lệch hội tụ đạt mức dưới $10^{-4}$. Dữ liệu tính toán sau đó được kết xuất và biểu diễn đồ họa bằng phần mềm MATLAB để trích xuất mặt đóng 3D và các đường cong cân bằng 2D. Quá trình nghiên cứu được triển khai đồng bộ từ tháng 3/2003 đến tháng 3/2005.

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

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

Thứ nhất, trong điều kiện triều đại dương đối xứng hoàn toàn ($h_{o1} = h_{o2} = 0.50\text{ m}$), mô hình vịnh hai cửa đơn thuần luôn tạo ra cấu hình Type D với điểm cân bằng ổn định có điều kiện tại $(A_1, A_2) = (3.500\text{ m}^2, 3.500\text{ m}^2)$ và giá trị đỉnh đạt $4.750\text{ m}^2$. Khi gia tăng chiều dài kênh dẫn từ $500\text{ m}$ lên $2.500\text{ m}$ và $4.000\text{ m}$ hoặc tăng hệ số ma sát lên $f = 0.006$, hệ thống suy thoái nhanh chóng từ Type D qua Type B và biến thành Type A, đồng thời giá trị mặt cắt đỉnh giảm 20% xuống còn $3.800\text{ m}^2$.

Thứ hai, nghiên cứu phát hiện sự bất đối xứng biên độ triều ngoài biển là động lực then chốt tạo ra trạng thái cân bằng Type E. Khi biên độ triều tại cửa 1 tăng 18% từ $0.50\text{ m}$ lên $0.59\text{ m}$, miền ổn định của cửa 1 mở rộng thêm 18% lên $8.300\text{ m}^2$, giá trị đỉnh của đường cong cân bằng cửa 1 tăng 32% đạt $6.250\text{ m}^2$, và đỉnh của cửa 2 tăng vọt 47.4% lên $7.000\text{ m}^2$. Sự thay đổi này tạo ra 4 giao điểm và thiết lập một điểm ổn định vô điều kiện tại tọa độ $(5.200\text{ m}^2, 4.200\text{ m}^2)$.

Thứ ba, đối với mô hình phân vùng vách ngăn nội vịnh, khi diện tích mặt cắt kênh ngăn $A_3$ biến thiên từ $2.000\text{ m}^2$ đến $18.500\text{ m}^2$ tương ứng với biên độ triều từ $0.10\text{ m}$ đến $2.60\text{ m}$, sự lệch pha mực nước giữa hai phân vịnh cho phép duy trì cấu hình Type E ổn định ngay cả khi tỷ lệ diện tích mặt nước giữa hai ngăn vịnh chênh lệch nhau từ 1 đến 10 lần.

Thứ tư, khi áp dụng mô hình cho hệ đầm phá Tam Giang - Cầu Hai, kết quả tính toán chứng minh rằng sau trận lũ lịch sử năm 1999 làm xuất hiện cửa Hòa Duẩn (hình thành hệ 3 cửa), lăng đạo triều qua cửa Thuận An bị suy giảm hơn 30%, khiến hệ rơi vào trạng thái mất ổn định cạnh tranh và giải thích nguyên nhân khoa học vì sao cửa Hòa Duẩn cần phải được hàn gắn nhân tạo để bảo vệ luồng tàu Thuận An.

Thảo luận kết quả

Nguyên nhân vật lý của việc hình thành cấu hình ổn định Type E bắt nguồn từ gradien thủy lực phụ trợ do chênh lệch pha và biên độ triều tạo ra. Khi mực nước biển tại hai cửa có sự chênh áp, một dòng chảy xuyên vịnh xuất hiện, giúp duy trì vận tốc đỉnh $\hat{u}$ luôn vượt ngưỡng xói rửa $1.0\text{ m/s}$ tại cả hai cửa họng, triệt tiêu xu hướng bồi lắng tự nhiên của cửa nhỏ hơn.

Để minh họa trực quan các phát hiện này, toàn bộ tập dữ liệu được thể hiện tối ưu qua hai hình thức đồ thị:

  • Mặt đóng ba chiều (3D Closure Surface): Thể hiện vận tốc dòng chảy đỉnh $\hat{u}$ trên trục $z$ biến thiên theo hai trục diện tích $(A_1, A_2)$, trong đó các vùng màu xanh lục phản ánh miền vận tốc vượt ngưỡng xói $\hat{u} \ge 1.0\text{ m/s}$.
  • Bản đồ đường cong cân bằng hai chiều (2D Equilibrium Flow Curves): Thể hiện mặt cắt mặt phẳng của 3D Closure Surface tại ngưỡng vận tốc $1.0\text{ m/s}$. Bảng ma trận tổng hợp chuyển tiếp các dạng dòng chảy từ A đến E theo mức biến thiên biên độ triều từ $0.50\text{ m}$ đến $0.70\text{ m}$ cho phép người kỹ sư nhận diện chính xác ranh giới ổn định.

So sánh với các nghiên cứu trước đây của Van de Kreeke (1990), kết quả của luận văn đã vượt qua các giả định giới hạn về triều đơn giản và vịnh đồng nhất, chứng minh rằng sự không đồng nhất về thủy triều và hình thái vịnh chính là điều kiện đủ để hai cửa biển cùng tồn tại bền vững trong tự nhiên.

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

  1. Thiết lập quy trình giám sát trắc đạc hình thái và thủy lực tự động: Triển khai lắp đặt hệ thống quan trắc ADCP đáy tự ghi kết hợp trạm đo mực nước thời gian thực tại các cửa biển trọng điểm; đặt mục tiêu kiểm soát sai số dự báo bồi xói lòng dẫn dưới 5%, thực hiện định kỳ 3 tháng/lần trong lộ trình 2026-2028, chủ thể thực hiện là Viện Khoa học Thủy lợi và Chi cục Đường thủy nội địa.
  2. Quy hoạch mặt cắt chỉnh trị luồng tàu theo cấu hình cân bằng Type E: Định hình diện tích mặt cắt luồng thiết kế cho hệ thống hai cửa trong khoảng $4.200\text{ m}^2$ đến $5.200\text{ m}^2$ phù hợp với chế độ chênh triều tự nhiên; hướng tới mục tiêu giảm 25% đến 30% khối lượng nạo vét duy tu luồng hàng hải giai đoạn 2026-2030, do Cục Hàng hải Việt Nam và các Ban Quản lý Dự án Hàng hải chủ trì.
  3. Chủ động kiểm soát và đóng kín các cửa mở đột biến sau bão lũ: Ban hành quy chế ứng phó khẩn cấp để hàn gắn các cửa biển mới mở do bão (tương tự trường hợp cửa Hòa Duẩn) trong vòng 6 đến 12 tháng nhằm bảo tồn trên 80% lăng đạo triều cho các cửa chính; giao UBND tỉnh Thừa Thiên Huế cùng Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn phối hợp thực hiện.
  4. Ứng dụng mô hình phân vùng vách ngăn trong quy hoạch công trình đầm phá: Bắt buộc áp dụng mô hình toán thủy động lực học có xét đến phân vùng nội vịnh khi thiết kế các đê ngăn mặn, cầu cạn hoặc kè chỉnh trị; đảm bảo tỷ lệ lưu thông nước giữa các phân vịnh đạt tối thiểu 85% thể tích triều, triển khai áp dụng từ quý 1/2026 bởi các đơn vị tư vấn thiết kế công trình thủy sản và giao thông thủy.

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

  • Kỹ sư công trình thủy và chỉnh trị bờ biển: Khai thác phương pháp xác định diện tích mặt cắt cân bằng $A_{eq}$ và quy trình phân loại Type A-E để ứng dụng trực tiếp vào thiết kế đê chắn sóng, mỏ hàn và luồng hàng hải đa cửa mở.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Thủy lực - Thủy văn và Kỹ thuật Bờ biển: Sử dụng hệ phương trình vi phân phi tuyến mở rộng, thuật toán lặp số học và khung phân tích Escoffier đa cửa làm tài liệu tham khảo chuyên sâu cho các đề tài nghiên cứu nâng cao.
  • Cơ quan quản lý nhà nước về hàng hải, đê điều và phòng chống thiên tai: Vận dụng các kết luận thực nghiệm từ bài toán đầm phá Tam Giang - Cầu Hai để hoạch định chính sách ứng phó xâm thực bờ biển, quy hoạch không gian đầm phá và quyết định phương án can thiệp cửa biển hợp lý.
  • Chuyên gia tư vấn quy hoạch môi trường và sinh thái đầm phá: Tham khảo cơ chế trao đổi nước và chênh lệch pha thủy triều nội vịnh để đánh giá chu kỳ thay nước (flushing time), phục vụ việc bảo tồn đa dạng sinh học và phát triển nuôi trồng thủy hải sản bền vững.

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

Tại sao hệ thống hai cửa biển thông cùng một vịnh thường có xu hướng một cửa bị bồi lấp?
Dưới tác động của thủy triều đại dương đồng nhất, cửa có diện tích lớn hơn sẽ chiếm ưu thế thu hút lăng đạo triều, làm giảm lưu lượng và kéo vận tốc đỉnh của cửa nhỏ hơn xuống dưới ngưỡng cân bằng $1.0\text{ m/s}$. Điều này gây ra quá trình bồi tụ cát liên tục và dẫn đến việc cửa nhỏ bị đóng kín hoàn toàn, như đã diễn ra tại kênh Pass Cavallo ở vịnh Matagorda (Mỹ).

Điều kiện cụ thể nào giúp hệ hai cửa đạt trạng thái ổn định vô điều kiện (Type E)?
Hệ thống đạt trạng thái Type E khi có sự bất đối xứng biên độ triều ngoài biển giữa hai cửa từ 18% trở lên (ví dụ tăng từ $0.50\text{ m}$ lên $0.59\text{ m}$) hoặc có sự lệch pha mực nước bên trong các phân vịnh. Sự chênh lệch thủy lực này tạo ra gradien dòng chảy phụ trợ, duy trì vận tốc xói cát đồng thời ở cả hai cửa dẫn.

Ý nghĩa của giá trị vận tốc cân bằng $u_{eq} = 1.0\text{ m/s}$ trong mô hình là gì?
Vận tốc cân bằng $u_{eq} = 1.0\text{ m/s}$ tương ứng với ứng suất tiếp đáy $\tau_{eq}$ từ $3.5\text{ N/m}^2$ đến $5.5\text{ N/m}^2$, là giá trị dòng chảy cực đại cần thiết để dòng triều rút cuốn trôi toàn bộ lượng cát do dòng ven bờ bồi đắp vào cửa họng. Nếu vận tốc thực tế thấp hơn $1.0\text{ m/s}$, kênh dẫn sẽ chuyển sang chế độ bồi cạn.

Mô hình phân vùng vách ngăn (partition model) giải quyết bài toán thực tế nào?
Mô hình vách ngăn khắc phục hạn chế của giả định mực nước dao động phẳng đồng nhất, cho phép mô phỏng chính xác các đầm phá trải dài hoặc bị chia cắt bởi địa hình ngầm có diện tích cửa ngăn $A_3$ từ $2.000\text{ m}^2$ đến $18.500\text{ m}^2$, qua đó phản ánh đúng dao động lệch pha mực nước thực tế.

Tại sao sự xuất hiện của cửa Hòa Duẩn năm 1999 lại làm suy thoái cửa Thuận An?
Cửa Hòa Duẩn mở ra đột ngột sau lũ lớn đã làm phân tán hơn 30% lăng đạo triều của đầm phá, kéo vận tốc dòng chảy qua cửa chính Thuận An xuống dưới mức $1.0\text{ m/s}$. Tình trạng mất cân bằng động lực học này khiến luồng Thuận An đứng trước nguy cơ bồi lấp nghiêm trọng, buộc chính quyền phải xử lý hàn khẩu cửa phụ.

Kết luận

  • Luận văn đã chứng minh thành công sự tồn tại của trạng thái ổn định mặt cắt vô điều kiện (dạng Type E) trong hệ thống vịnh hai cửa biển nhờ cơ chế vi sai biên độ triều (từ 18%) và phân vùng nội vịnh.
  • Xây dựng hoàn chỉnh mô hình tính toán giải tích kết hợp lặp số học trên lưới ma trận 2.500 điểm dữ liệu, giải quyết triệt để hệ phương trình thủy động lực học phi tuyến có xét đến ma sát và quán tính.
  • Kiểm chứng thành công độ tin cậy của thuật toán thông qua dữ liệu thực địa tại vịnh Matagorda (Mỹ) với sự trùng khớp cao về tọa độ điểm cân bằng hình thái.
  • Làm sáng tỏ cơ sở khoa học về động lực học hình thái đầm phá Tam Giang - Cầu Hai (Việt Nam), giải thích thỏa đáng sự cần thiết của việc hàn gắn cửa Hòa Duẩn để duy trì ổn định luồng tàu Thuận An.
  • Cung cấp cơ sở định lượng quan trọng giúp tối ưu hóa chi phí nạo vét và nâng cao hiệu quả kỹ thuật cho các dự án chỉnh trị cửa sông ven biển từ 25% đến 30%.

Trong giai đoạn tiếp theo (2026-2027), mô hình cần được phát triển tích hợp thêm tương tác sóng - dòng chảy 2D/3D và dòng chảy lũ sông để hoàn thiện công cụ dự báo hình thái bờ biển. Các đơn vị tư vấn và cơ quan quản lý hàng hải được khuyến nghị ứng dụng ngay phương pháp luận phân loại Type A-E của luận văn vào công tác quy hoạch cảng biển và bảo vệ đầm phá ven bờ.