Tổng quan về luận án

Nghiên cứu lựa chọn quy mô công trình dẫn dòng thi công trong xây dựng công trình thủy lợi, thủy điện ở Việt Nam là một luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng Công trình Thủy (Mã số: 9580202) do nghiên cứu sinh Mai Lâm Tuấn thực hiện tại Trường Đại học Thủy lợi dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Quang Cường và PGS.TS. Lê Văn Hùng (2019). Luận án giải quyết trực tiếp một trong những bài toán phức tạp, tốn kém và tiềm ẩn rủi ro sự cố cao nhất trong kỹ thuật xây dựng công trình đầu mối: tối ưu hóa kinh tế - kỹ thuật quy mô các hạng mục dẫn dòng xả lũ tạm thời trong suốt chu kỳ xây dựng đập.

                                  QUY HOẠCH DẪN DÒNG THI CÔNG
   ĐẶC TRƯNG KHÍ TƯỢNG - THỦY VĂN                                      PHƯƠNG ÁN DẪN DÒNG LIÊN HỢP
  (5 vùng khí hậu: Bắc Bộ, Bắc Trung Bộ,                            (Cống / Hầm kết hợp Đập dở dang)
   TẦN SUẤT THIẾT KẾ DẪN DÒNG P(%)                                  (Xử lý chuyển tiếp Chảy có áp/Không áp,
   (Phân cấp theo H cột nước & W hồ)                                  Chảy ngập/Tự do, Đường mặt nước & Vmax)
                                   TỐI ƯU HÓA KINH TẾ - KỸ THUẬT
                         (Cân bằng chi phí công trình tạm vs Rủi ro ngập lụt)

Khoảng trống nghiên cứu (Research gap) then chốt mà luận án nhận diện xuất phát từ các bất cập mang tính hệ thống trong các quy chuẩn, tiêu chuẩn hiện hành tại Việt Nam, cụ thể là Quy chuẩn Kỹ thuật Quốc gia QCVN 04-05:2012/BNN&PTNT và Tiêu chuẩn Quốc gia TCVN 9160:2012 ("Công trình thủy lợi - Yêu cầu thiết kế dẫn dòng trong xây dựng"). Các văn bản này quy định tần suất lưu lượng thiết kế dẫn dòng ($P%$) thuần túy dựa vào cấp của công trình tạm thời, chưa phản ánh đầy đủ và định lượng được mức độ rủi ro khi công trình chính tham gia chắn nước hoặc xả lũ dở dang qua thân đập. Tiêu chuẩn hiện hành thiếu cơ sở khoa học định lượng để người kỹ sư nâng hoặc hạ cấp tần suất thiết kế căn cứ vào chiều cao cột nước dâng ($H$), dung tích phòng lũ tức thời của lòng hồ ($W$), và mức độ thiệt hại kinh tế - xã hội ở vùng hạ du khi xảy ra sự cố vỡ đê quai hay vỡ đập dở dang. Đồng thời, về mặt thủy lực, các phương pháp tính toán truyền thống trong TCVN 9160:2012 chỉ áp dụng cho từng công trình dẫn dòng độc lập, thiếu thuật toán liên hợp giải quyết đồng thời bài toán thủy lực qua cống/hầm và tràn qua đỉnh đập đang xây dựng có xét đến quá trình điều tiết lũ và hiện tượng chuyển tiếp chế độ dòng chảy.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ sở khoa học nào cho phép định lượng hóa việc lựa chọn tần suất thiết kế ($P%$) và phân chia thời đoạn thi công tối ưu thích ứng với đặc thù phân phối dòng chảy lũ - kiệt của 5 vùng khí hậu thủy văn Việt Nam?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật diễn biến trường lưu tốc lớn nhất ($V_{\max}$) trên bề mặt dốc nước của đập đá đổ đang xây dựng chịu ảnh hưởng như thế nào bởi các thông số hình học (hệ số mái dốc $m$, chiều dài dốc nước $L$, chênh lệch mực nước hạ lưu $H_h$)?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xây dựng một thuật toán thủy lực liên hợp giải quyết triệt để hiện tượng bước nhảy số học tại điểm chuyển tiếp giữa các chế độ chảy ngập - không ngập, chảy không áp - có áp khi đồng thời xả lũ qua cống đáy và tràn qua thân đập dở dang?
  • Giả thuyết khoa học 1 (H1): Việc nâng/hạ cấp tần suất dẫn dòng phải là một hàm phụ thuộc trực tiếp vào dung tích hồ và chiều cao công trình chắn dòng tại thời điểm dẫn dòng, thay vì cố định theo cấp công trình tạm.
  • Giả thuyết khoa học 2 (H2): Phân bố lưu tốc cực đại $V_{\max}$ trên bề mặt đập đá đổ cho nước tràn qua có thể kiểm soát và định lượng chính xác thông qua mô hình thủy lực 1D cải tiến tích hợp điều kiện biên 2D/3D, làm cơ sở tối ưu hóa kết cấu gia cố tấm bê tông cốt thép kết hợp neo.
  • Giả thuyết khoa học 3 (H3): Mô hình tính toán điều tiết lũ liên hợp cống - đập có thể giải quyết hội tụ trơn tru các phương trình thủy lực chuyển tiếp, cho phép giảm quy mô kích thước công trình dẫn dòng ngầm đắt tiền mà vẫn đảm bảo tuyệt đối an toàn thi công.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng Thủy lực học công trình kinh điển (phương trình năng lượng Bernoulli, phương trình động lượng, phương trình liên tục), Lý thuyết dòng chảy không ổn định một chiều và hai chiều (hệ phương trình Saint-Venant), kết hợp với Lý thuyết phân tích kinh tế - rủi ro công trình tạm. Luận án mang lại tác động định lượng đột phá: đề xuất khung tần suất thiết kế mới cho đập tràn dở dang với dải tần suất $P = 1% - 5%$ tùy thuộc vào dung tích trữ và chiều cao đập; phát triển phần mềm chuyên dụng tính toán thủy lực dẫn dòng liên hợp giúp rút ngắn thời gian tính toán từ nhiều tuần xuống vài phút, cho phép phân tích so sánh hàng chục kịch bản kinh tế - kỹ thuật; tối ưu hóa hàng trăm tỷ đồng chi phí xây dựng công trình dẫn dòng và gia cố bảo vệ cho các cụm công trình thủy điện lớn tại Việt Nam.

Phạm vi nghiên cứu bao quát chuỗi số liệu thực tế tại các công trình trọng điểm quốc gia như Thủy điện Sơn La ($2.400\text{ MW}$, $Q_{\text{dd}} = 16.044\text{ m}^3/\text{s}$), Lai Châu ($1.200\text{ MW}$, $Q_{\text{dd}} = 10.388\text{ m}^3/\text{s}$), Tuyên Quang ($342\text{ MW}$, $Q = 5.036\text{ m}^3/\text{s}$), Cửa Đạt ($97\text{ MW}$), Hòa Bình ($1.920\text{ MW}$), Huội Quảng ($520\text{ MW}$), Bắc Hà ($90\text{ MW}$), kết hợp khảo sát 5 vùng thủy văn trên phạm vi toàn lãnh thổ Việt Nam.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về công tác dẫn dòng thi công cho thấy ba dòng nghiên cứu (research streams) chính:

                               CÁC DÒNG NGHIÊN CỨU TRONG LITERATURE
  STREAM 1: XẢ LŨ QUA ĐẬP ĐÁ ĐỔ         STREAM 2: THỦY LỰC XẢ KẾT HỢP         STREAM 3: CỬA VÀO CÔNG TRÌNH NGẦM
  - Henry Olivier (1967): q, m, d       - Hoàng Văn Tân (1999): Cơ chế dòng   - Nguyễn Danh Oanh (2003): Áp suất
  - Vũ Trọng Hồng (1972): Neo mái         tràn mặt + xả đáy qua cống            chân không, trụ nhô giảm áp
  - Trần Quốc Thưởng (2016): Kết cấu    - Lê Bá Sơn (1994): Tổn thất va chạm  - Sổ tay Thủy lực Nga (1985):
    tấm BTCT kết hợp neo                  dòng xả (chiếm 20% - 50%)              Tính toán công trình độc lập
  1. Dòng nghiên cứu dòng chảy tràn qua đập đá đổ và gia cố mái hạ lưu: Henry Olivier (1967) đặt nền móng lý thuyết thiết lập mối quan hệ giữa lưu lượng đơn vị $q$, độ dốc mái $m$ và đường kính hạt đá $d$ khi cho dòng chảy tràn qua đập đá đổ tự nhiên. Tại Việt Nam, luận án của Vũ Trọng Hồng (1972) chứng minh rằng khi mái hạ lưu đập đá đổ dốc hơn $1:3$, lực thủy động của dòng tràn sẽ phá vỡ liên kết trọng lực của khối đá, bắt buộc phải sử dụng kết cấu neo thép đặc biệt để giữ ổn định. Kế thừa hướng đi này, đề tài cấp nhà nước của Trần Quốc Thưởng (2016) đã xác định các thông số kết cấu gia cố bề mặt đập đá đổ bằng hệ thống tấm bê tông cốt thép kết hợp neo chốt sâu vào thân đập và rọ đá bảo vệ mái hạ lưu.
  2. Dòng nghiên cứu thủy lực xả kết hợp giữa dòng tràn mặt và tháo sâu: Luận án của Hoàng Văn Tân (1999) nghiên cứu thực nghiệm khả năng tháo kết hợp giữa tràn đỉnh rộng và cống đáy, chỉ ra rằng tương tác giữa dòng chảy xiết dưới đáy và dòng tràn tự do trên mặt làm thay đổi hệ số lưu lượng của từng công trình so với trạng thái làm việc độc lập. Đồng thời, nghiên cứu của Lê Bá Sơn (1994) lượng hóa tổn thất năng lượng do sự giao thoa, va đập giữa dòng phun tự do và dòng tháo sâu đáy, chứng minh tổn thất cục bộ do tương tác dòng chảy chiếm từ $20%$ đến $50%$ tổng tổn thất năng lượng từ thượng lưu về hạ lưu.
  3. Dòng nghiên cứu thủy động lực học công trình ngầm và cống dẫn dòng: Nguyễn Danh Oanh (2003) giải quyết bài toán thủy động lực học tại cửa vào đường hầm dẫn dòng thi công, đề xuất kết cấu cửa vào có trụ nhô giúp triệt tiêu áp suất chân không và hiện tượng khí thực (cavitation) nguy hiểm trên trần cửa hầm.

Trong y văn tồn tại cuộc tranh luận học thuật sâu sắc giữa hai quan điểm đối nghịch về chiến lược dẫn dòng mùa lũ:

  • Trường phái An toàn thụ động (Conservative Structural Isolation): Đại diện bởi các quy chuẩn xây dựng truyền thống của Liên Xô cũ và các sổ tay thủy lực công trình cổ điển (như Sổ tay Thủy lực Nga, 1985; Fisher & Ramsbottom, 2001). Quan điểm này ưu tiên tuyệt đối việc bảo vệ thân đập chính bằng cách cô lập hoàn toàn hố móng trong mùa lũ, chấp nhận đào các đường hầm xả lũ tạm thời kích thước khổng lồ xuyên qua núi đá hoặc đắp đê quai vượt lũ cực cao. Hạn chế cốt tử là đẩy chi phí công trình tạm lên mức $15% - 25%$ tổng mức đầu tư dự án và để lại gánh nặng hoành triệt vô cùng phức tạp.
  • Trường phái Tối ưu hóa chủ động (Dynamic Overtopping & Hybrid Diversion): Được thúc đẩy mạnh mẽ qua các thực tiễn xây dựng đập lớn tại Trung Quốc (như đập CFRD Thiên Sinh Kiều cao $178\text{ m}$, Thủy Bố Á cao $233\text{ m}$, Giao Khê cao $132,5\text{ m}$) và một số công trình tại Mỹ. Trường phái này chủ trương tận dụng chính thân đập đá đổ hoặc bê tông đang thi công dở dang làm công trình tràn xả lũ mùa lũ kết hợp với cống/hầm dẫn dòng đáy.

Nghiên cứu của Mai Lâm Tuấn định vị chính xác tại điểm giao cắt giữa hai trường phái: khoa học hóa giải pháp dẫn dòng kết hợp thông qua việc xây dựng mô hình toán thủy lực liên hợp có khả năng mô phỏng quá trình điều tiết lũ, phân tích trường lưu tốc bề mặt để đưa ra giải pháp gia cố neo tối ưu, xóa bỏ sự mơ hồ định tính trong các quy định nâng hạ cấp tần suất thiết kế dẫn dòng của QCVN 04-05:2012.

So sánh với các nghiên cứu và thực tiễn công trình quốc tế tiêu biểu:

  • Công trình đập Braunla (Sông Snake, Mỹ): Đập cao $88,5\text{ m}$, dòng lũ thiết kế $2.000\text{ m}^3/\text{s}$. Người Mỹ đã lựa chọn giải pháp không đào thêm hầm thứ hai mà cho lũ tràn qua đập đá đổ đang xây dựng ở cao trình $+552,0\text{ m}$ với lớp nước tràn dày tới $6,1\text{ m}$, bề mặt đập được bảo vệ bằng lớp đá đầm nén chọn lọc dày $2,4\text{ m}$.
  • Công trình đập Ord (Úc): Đập đá đổ lõi giữa cao $98\text{ m}$, dẫn dòng mùa lũ cho tràn qua thân đập ở cao trình $+41,0\text{ m}$ và $+48,8\text{ m}$. Bề mặt đỉnh đập được gia cố bằng bê tông dày $1,0\text{ m}$, mái hạ lưu gia cố bằng hệ neo thép $\varnothing 25\text{ mm}$ (chiều dài neo lên tới $10\text{ m}$, mật độ $1,4 \times 0,45\text{ m}$/neo) kết hợp lưới rọ đá $152 \times 152\text{ mm}$, chịu được lưu tốc dòng tràn thực tế đạt gần $4,5\text{ m}/\text{s}$ an toàn tuyệt đối.

Đó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 trực tiếp lý thuyết thủy lực công trình và lý thuyết quản trị rủi ro công trình thủy:

  1. Mở rộng lý thuyết thủy lực dốc nước trên bề mặt nhám phi quy ước: Công trình hoàn thiện cơ sở lý thuyết tính toán phân bố lưu tốc lớn nhất ($V_{\max}$) dọc theo chiều dài dốc tràn bề mặt đập đá đổ. Thiết lập mối quan hệ giải tích phụ thuộc giữa $V_{\max}$ với 4 biến độc lập: lưu lượng đơn vị $q$ ($\text{m}^3/\text{s/m}$), hệ số mái dốc $m$, chiều dài dốc nước $L$ ($\text{m}$) và độ chênh lệch mực nước cuối dốc so với mực nước hạ lưu $H_d - H_h$.
  2. Khái niệm hóa lại rủi ro trong lý thuyết tần suất thiết kế công trình tạm: Luận án thách thức giả định truyền thống xem công trình dẫn dòng là "công trình tạm thuần túy". Đưa ra luận điểm khoa học chứng minh khi thân đập chính hoặc đê quai cao tham gia chắn nước, tính chất rủi ro chuyển dịch từ "rủi ro kinh tế cục bộ hố móng" sang "rủi ro an toàn đập vùng hạ du", đòi hỏi tần suất thiết kế phải tiệm cận tần suất của công trình chính.
                           KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN

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

Khung phân tích của luận án là sự tích hợp hữu cơ giữa ba trụ cột:

  • Trụ cột Thủy lực học: Tích hợp mô hình phương trình bảo toàn năng lượng và phương trình liên tục cho hệ thống hai hay nhiều công trình cùng làm việc. Thuật toán phân tách dòng chảy thành các phân đoạn, tự động phát hiện và chuyển đổi trạng thái dòng chảy: từ chảy tự do đỉnh rộng sang chảy ngập, từ cống hở chảy không áp sang cống ngập chảy có áp với các hệ số lưu lượng $\mu$, $\varphi$, hệ số ngập $\sigma_n$ được hiệu chỉnh động.
  • Trụ cột Cơ học công trình: Xác định ứng suất kéo và lực đẩy nổi mạch động thủy lực tác động lên các tấm bê tông lát mái hạ lưu, làm cơ sở tính toán đường kính, chiều dài và mật độ bố trí neo thép chịu lực.
  • Trụ cột Kinh tế - Rủi ro: Hàm mục tiêu tối thiểu hóa tổng chi phí $C_{\text{total}} = C_{\text{xd}} + C_{\text{gc}} + C_{\text{ht}} + E(\text{Damage}) \to \min$, trong đó cân bằng chi phí xây dựng công trình dẫn dòng ($C_{\text{xd}}$), chi phí gia cố đập ($C_{\text{gc}}$), chi phí hoành triệt ($C_{\text{ht}}$) và tổn thất rủi ro kỳ vọng $E(\text{Damage})$.
  • Điều kiện biên áp dụng (Boundary conditions): Áp dụng tối ưu cho đập đá đổ bản mặt bê tông (CFRD), đập bê tông đầm lăn (RCC) và đập đất đá hỗn hợp có địa hình tuyến hẹp chữ V hoặc chữ U, nơi lưu lượng lũ mùa hè vượt quá $1.000\text{ m}^3/\text{s}$ và chênh lệch lưu lượng mùa kiệt - lũ lớn hơn 5 lần.

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

Thiết kế nghiên cứu

Luận án vận dụng triết lý nghiên cứu thực chứng (Positivism) kết hợp mô hình hóa giải tích - số (Analytical-Numerical Modeling) và kiểm chứng thực nghiệm (Empirical Validation). Thiết kế nghiên cứu đa cấp độ (Multi-level research design) được triển khai qua 3 tầng phân tích:

  1. Cấp độ vĩ mô: Phân tích chuỗi số liệu khí hậu thủy văn 5 vùng địa lý toàn quốc để xác định phân phối mùa kiệt - lũ và quy luật xuất hiện lũ tiểu mãn.
  2. Cấp độ trung mô: Mô phỏng thủy lực tương tác hệ thống công trình dẫn dòng (cống ngầm + đập dở dang + hồ chứa điều tiết).
  3. Cấp độ vi mô: Phân tích trường lưu tốc thủy động lực 2D/3D trên bề mặt mái dốc đập đá đổ và vùng tiêu năng hạ lưu.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ nghiêm ngặt các bước logic khoa học:

  • Chiến lược chọn mẫu dữ liệu: Toàn bộ dữ liệu khí tượng thủy văn, mặt cắt địa hình và thông số kỹ thuật được trích xuất từ các hồ sơ thiết kế kỹ thuật chính thức và báo cáo hoàn công của các công trình thủy điện lớn cấp đặc biệt và cấp I tại Việt Nam (Sơn La, Lai Châu, Tuyên Quang, Cửa Đạt, Hòa Bình, Bản Chát, Huội Quảng).
  • Tam giác đạc phương pháp (Methodological Triangulation): Kết quả tính toán của chương trình máy tính tự lập được đối chứng chéo đồng thời với 3 nguồn kiểm chứng độc lập:
    1. Số liệu đo đạc thực tế tại hiện trường các đợt xả lũ qua công trình thật (Thủy điện Sơn La năm 2009, Thủy điện Tuyên Quang năm 2006, Thủy điện Hòa Bình năm 1986).
    2. Kết quả thí nghiệm mô hình vật lý thủy lực (Physical scale model tests) tỷ lệ $1/50$ và $1/100$ thực hiện tại phòng thí nghiệm trọng điểm quốc gia.
    3. Mô phỏng số 3 chiều chi tiết sử dụng phần mềm thủy động lực học thương vi tính cao cấp Flow-3D.
                           TAM GIÁC ĐẠC PHƯƠNG PHÁP CỦA LUẬN ÁN

Data và phân tích

Luận án tiến hành xây dựng các thuật toán chuyên dụng viết trên nền tảng ngôn ngữ lập trình kỹ thuật:

  • Thuật toán tính thủy lực cống tự động: Xử lý điều kiện chuyển tiếp $Q - H_0$ giữa 4 trạng thái: cống không áp chảy tự do, cống không áp chảy ngập, cống có áp chảy tự do và cống có áp chảy ngập. Khắc phục hoàn toàn hiện tượng phân kỳ hay gián đoạn số học (numerical oscillation) tại các điểm uốn của đường cong quan hệ lưu lượng - mực nước.
  • Thuật toán điều tiết lũ thi công: Giải phương trình cân bằng nước vi phân $\frac{dW}{dt} = Q_{\text{vào}}(t) - Q_{\text{ra}}(Z)$ theo phương pháp sai phân hữu hạn từng bước thời gian $\Delta t$, tự động điều chỉnh lưu lượng tháo liên hợp $Q_{\text{ra}} = Q_{\text{cống}}(Z) + Q_{\text{đập}}(Z)$.
  • Mô phỏng 3D bằng Flow-3D: Giải hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp mô hình rối $k-\varepsilon$ hai phương trình và thuật toán Volume of Fluid (VOF) để theo dõi bề mặt tự do của dòng tràn trên dốc đập đá đổ. Kết quả so sánh giữa chương trình tự lập và thí nghiệm mô hình vật lý tại công trình Thủy điện Sơn La và Tuyên Quang cho thấy sai số lưu lượng $\Delta Q < 3,5%$, sai số mực nước thượng lưu $\Delta Z < 2,0%$, chứng minh độ tin cậy và độ vững chắn (robustness) tuyệt đối của thuật toán.

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

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

                               4 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN
  1. Định lượng hóa cơ sở nâng/hạ tần suất thiết kế dẫn dòng khi đập chính tham gia: Bác bỏ quan điểm áp dụng cơ học tần suất công trình tạm ($P = 5% - 10%$) khi thân đập chính đã thi công lên cao. Luận án chỉ ra rằng khi chiều cao đập chắn nước vượt quá $50\text{ m}$ và dung tích hồ trữ tạm thời vượt trên $100 \text{ triệu m}^3$, rủi ro vỡ đập sẽ gây thảm họa khôn lường cho hạ du. Do đó, tần suất thiết kế bắt buộc phải nâng cấp lên $P = 1% - 2%$ (tương đương cấp công trình chính). Ngược lại, ở giai đoạn hố móng ban đầu, chấp nhận hạ cấp xuống $P = 10% - 20%$ và cho nước tràn ngập hố móng có kiểm soát sẽ đem lại hiệu quả kinh tế tối ưu.
  2. Quy luật phi tuyến tính của trường lưu tốc cực đại $V_{\max}$ trên mái đập: Dữ liệu tính toán và mô phỏng chỉ ra rằng khi hệ số mái dốc $m \ge 8$, dòng chảy có xu hướng tiêu tán năng lượng nội tại rất lớn dọc theo dốc tràn, làm giảm đáng kể tốc độ gia tăng của $V_{\max}$. Khi chênh lệch độ cao cuối dốc nước so với mực nước hạ lưu $H_d - H_h \le 0$ (trạng thái ngập hoàn toàn), $V_{\max}$ giảm từ $35%$ đến $50%$ so với trạng thái chảy tự do vào không khí.
  3. Giải mã nguyên nhân sự cố vỡ đập dở dang tại công trình Cửa Đạt (2007): Bằng mô hình thủy lực, luận án chỉ rõ nguyên nhân thất bại tại Cửa Đạt khi cho lũ tràn qua đập đá đổ dở dang: lựa chọn cao trình ngưỡng tràn quá cao ($+50,0\text{ m}$), tạo độ chênh cột nước sau đập lên tới $17,5\text{ m}$, vượt quá khả năng chịu lực của thảm rọ đá dày $0,5\text{ m}$. Lực đẩy nổi mạch động và hiện tượng thấm rối qua khối đá đổ không được cố kết bằng neo sâu đã phá hủy thảm rọ đá từ chân khay, dẫn đến xói trôi thân đập. Ngược lại, bài học thành công tại Thủy điện Tuyên Quang là do đặt cao trình dẫn dòng thấp ($+48,0\text{ m}$), khống chế cột nước rơi chỉ $\approx 4\text{ m}$, kết hợp phủ bản mặt bê tông cốt thép.
  4. Cơ chế tương tác thủy lực cống đáy - tràn mặt: Khi cống dẫn dòng và tràn đập cùng hoạt động, dòng xả đáy tốc độ cao tạo ra hiệu ứng hút (sink effect), làm hạ thấp đường mặt nước thượng lưu trước ngưỡng tràn đập dở dang, giúp duy trì chế độ chảy ngập trên đỉnh đập lâu hơn, giảm vận tốc va đập và giảm áp lực xói mòn lên kết cấu gia cố đỉnh đập.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình toán học tổng quát giải quyết bài toán thủy lực liên hợp đa công trình có xét đến biến đổi trạng thái chảy và quá trình điều tiết lũ hồ chứa.
  • Về mặt phương pháp luận: Thiết lập quy trình tính toán tự động hóa hoàn chỉnh thay thế các phương pháp tra bảng và tính toán thủ công rời rạc, tạo chuẩn mực phương pháp luận mới cho tính toán dẫn dòng thi công.
  • Về mặt ứng dụng thực tiễn: Cung cấp giải pháp kỹ thuật thiết kế kết cấu gia cố bề mặt đập đá đổ bằng tấm bê tông cốt thép kết hợp neo chốt sâu (chiều dài neo $L_{\text{neo}} = 4 - 10\text{ m}$, đường kính neo $d \ge 25\text{ mm}$), cho phép xả an toàn lưu tốc dòng tràn lên tới $V = 12 - 15\text{ m}/\text{s}$.
  • Về mặt chính sách và quy chuẩn: Đóng góp luận cứ khoa học trực tiếp để Bộ Nông nghiệp và PTNT, Bộ Xây dựng rà soát, sửa đổi bổ sung QCVN 04-05:2012 và TCVN 9160:2012, đưa bảng phân cấp tần suất thiết kế dẫn dòng theo chiều cao cột nước và dung tích hồ vào quy phạm bắt buộc.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu khoa học:

  1. Giới hạn về mô hình thấm rối: Chưa tích hợp bài toán thấm rối phi tuyến ba chiều (Non-Darcy turbulent seepage) qua khối đá đổ thân đập trong điều kiện mực nước dâng nhanh đột ngột ở thượng lưu.
  2. Giới hạn về phân tích xung áp mạch động: Mặc dù đã tính đến lưu tốc trung bình $V_{\max}$, nghiên cứu chưa mô hình hóa chi tiết trường áp suất mạch động ngẫu nhiên (turbulent pressure fluctuation spectrum) tác động lên từng liên kết mối nối giữa các tấm bê tông lát mái dốc.
  3. Giới hạn hình thái thung lũng sông: Mô hình tập trung chủ yếu vào dạng thung lũng sông hẹp và vừa (đặc trưng miền núi Việt Nam), chưa mở rộng cho các hình thái lòng sông phân nhánh phức tạp vùng đồng bằng.

Chương trình nghiên cứu tương lai (Future research agenda):

  • Phát triển mô hình ghép đôi chất lưu - kết cấu đa pha (Fluid-Structure Interaction - FSI) sử dụng phương pháp hạt mịn mượt SPH (Smoothed Particle Hydrodynamics) để mô phỏng chính xác sự phá hủy của từng thanh neo và tấm bê tông dưới tác động của dòng tràn vận tốc siêu cao ($V > 18\text{ m}/\text{s}$).
  • Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Learning) để dự báo thời gian thực chế độ đóng mở cống kết hợp xả tràn dựa trên chuỗi dự báo mưa vệ tinh phân giải cao, tối ưu hóa quy trình điều hành dẫn dòng thi công thích ứng với biến đổi khí hậu cực đoan.
  • Mở rộng nghiên cứu giải pháp gia cố mái đập dở dang bằng vật liệu mới như bê tông tự lèn sợi polyme (FRC) và thảm bê tông địa kỹ thuật tự đông cứng dưới nước.

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

  • Ảnh hưởng học thuật (Academic Impact): Luận án là tài liệu tham khảo cốt lõi, mở ra hướng nghiên cứu chuyên sâu về thủy lực công trình tạm thời trong các chương trình đào tạo sau đại học ngành Kỹ thuật Thủy lợi, Thủy điện. Ước tính công trình tạo tiền đề cho hàng chục bài báo khoa học trên các tạp chí chuyên ngành uy tín trong nước và quốc tế.
  • Tái thiết ngành công nghiệp xây dựng thủy năng (Industry Transformation): Tạo bước ngoặt công nghệ cho các tổng thầu xây lắp (Song Da Corporation, Vinaconex, EVN, Trường Sơn) trong việc chủ động chuyển đổi từ phương án đào hầm dẫn dòng tốn kém sang giải pháp dẫn dòng qua đập dở dang kết hợp gia cố neo. Ứng dụng trực tiếp cho các dự án mở rộng thủy điện quy mô lớn như Thủy điện Hòa Bình mở rộng ($480\text{ MW}$), Thủy điện Ialy mở rộng ($360\text{ MW}$), Thủy điện Trị An mở rộng và các dự án thủy điện của Việt Nam đầu tư tại Lào (Thủy điện Luang Prabang, Xekaman).
  • Lợi ích kinh tế - xã hội định lượng: Tiết kiệm hàng trăm tỷ đồng chi phí đào hầm và đổ bê tông cống tạm; rút ngắn tiến độ thi công toàn dự án từ 6 đến 12 tháng; triệt tiêu nguy cơ vỡ đê quai gây ngập lụt, bảo vệ an toàn tính mạng và tài sản cho hàng triệu người dân sinh sống tại lưu vực hạ du các công trình thủy lợi, thủy điện.

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

                               CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI

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? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc xây dựng mô hình giải tích - số mô tả tương tác thủy lực chuyển tiếp phi tuyến giữa cống xả đáy và dòng tràn qua đập đá đổ dở dang có xét đến quá trình điều tiết lũ hồ chứa. Nghiên cứu mở rộng trực tiếp lý thuyết dòng chảy không ổn định và lý thuyết dốc nước của Henry Olivier (1967) và Vũ Trọng Hồng (1972) từ trạng thái tĩnh đơn lẻ sang trạng thái động tương tác liên hợp đa công trình.
  2. Điểm đổi mới phương pháp luận so với ít nhất 2 nghiên cứu đi trước là gì? Trả lời: So với Sổ tay Thủy lực Nga (1985) và TCVN 9160:2012 (chỉ tính riêng rẽ từng công trình độc lập và sử dụng bảng tra tĩnh), phương pháp luận của luận án giải quyết đồng thời phương trình liên tục và năng lượng của hệ ghép nối, tự động khử bước nhảy số học khi chuyển tiếp chế độ chảy có áp/không áp. So với các mô hình 1D như Saint-Venant thuần túy (MIKE 11, HEC-RAS không chuyên cho thủy lực công trình dẫn dòng), chương trình tự lập của luận án tích hợp trực tiếp thuật toán điều tiết lũ với đường cong đặc tính công trình tháo, giảm thời gian thiết lập mô hình từ nhiều ngày xuống vài phút.
  3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là vai trò quyết định của mực nước ngập hạ lưu $H_h$ và hệ số mái $m \ge 8$ trong việc triệt tiêu lưu tốc cực đại $V_{\max}$ trên mái đập đá đổ. Dữ liệu chỉ ra rằng việc hạ thấp cao độ ngưỡng tràn để tạo chế độ chảy ngập ($H_d - H_h \le 0$) làm giảm lưu tốc $V_{\max}$ từ $15\text{ m}/\text{s}$ xuống dưới $8\text{ m}/\text{s}$, làm giảm tới $70%$ khối lượng và chi phí neo thép gia cố mái dốc mà vẫn đảm bảo an toàn tuyệt đối trước xói trôi.
  4. Luận án có cung cấp quy trình tái lập (Replication protocol) hoàn chỉnh không? Trả lời: Có. Luận án cung cấp đầy đủ sơ đồ khối thuật toán (Flowcharts), hệ thống phương trình toán học thiết lập điều kiện biên, bảng thông số hình học, hệ số nhám $n$, hệ số lưu lượng $\mu$, cũng như quy trình chi tiết từng bước tính toán thủy lực cống, thủy lực đập dở dang và điều tiết lũ hồ chứa, cho phép bất kỳ kỹ sư hoặc nhà nghiên cứu nào lập trình tái lập chính xác kết quả.
  5. Chương trình nghiên cứu 10 năm tới được phác thảo ra sao? Trả lời: Lộ trình 10 năm tập trung vào 3 trọng tâm: (1) Số hóa và tự động hóa thiết kế dẫn dòng bằng công nghệ BIM thủy công tích hợp AI; (2) Phát triển các tiêu chuẩn tính toán kết cấu gia cố đập bằng vật liệu polyme nanocomposite chịu mài mòn xâm thực; (3) Mở rộng khung phân tích dẫn dòng thích ứng với các kịch bản biến đổi khí hậu cực đoan và lũ quét lưu vực sông quốc tế xuyên biên giới.

Kết luận

  1. Luận án đã hệ thống hóa và làm sáng tỏ cơ sở khoa học toàn diện về công tác dẫn dòng thi công trong xây dựng công trình thủy lợi, thủy điện ở Việt Nam, gắn kết chặt chẽ điều kiện khí hậu thủy văn 5 vùng địa lý với quy mô công trình xả tạm.
  2. Đã lượng hóa tiêu chí nâng/hạ cấp tần suất thiết kế dẫn dòng $P(%)$, đề xuất khung tần suất khoa học dựa trên chiều cao cột nước dâng và dung tích hồ chứa khi đập chính tham gia dẫn dòng, khắc phục khoảng trống tồn tại trong QCVN 04-05:2012 và TCVN 9160:2012.
  3. Khám phá quy luật phân bố trường lưu tốc lớn nhất $V_{\max}$ trên bề mặt đập đá đổ đang xây dựng phụ thuộc vào các thông số hình học $m$, $L$, $H_d - H_h$, cung cấp công thức giải tích chuẩn xác phục vụ thiết kế kết cấu gia cố mái hạ lưu bằng tấm bê tông cốt thép kết hợp neo chốt.
  4. Xây dựng thành công thuật toán và chương trình phần mềm tính toán thủy lực dẫn dòng liên hợp cống - đập dở dang tích hợp điều tiết lũ hồ chứa, xử lý triệt để bài toán chuyển tiếp chế độ dòng chảy không gây bước nhảy số học.
  5. Kiểm chứng thành công độ tin cậy của mô hình qua số liệu thực tế tại các đại công trình thủy điện Sơn La, Tuyên Quang, Hòa Bình, Cửa Đạt, chứng minh tiềm năng tiết kiệm hàng trăm tỷ đồng và rút ngắn tiến độ thi công từ 6 đến 12 tháng.
  6. Mở ra ba hướng nghiên cứu mới mang tầm quốc tế: mô phỏng tương tác động lực chất lưu - kết cấu (FSI) tấm neo gia cố đập dở dang, điều hành dẫn dòng thông minh thời gian thực ứng dụng AI, và chuẩn hóa khung quy phạm thiết kế dẫn dòng thi công công trình thủy thích ứng với biến đổi khí hậu toàn cầu.