Giới thiệu dự án

Dự án "Thiết kế hồ chứa nước Phước Quang – Phương án 1" là công trình thủy lợi – thủy điện trọng điểm thuộc chuyên ngành Kỹ thuật Công trình Thủy (Trường Đại học Thủy Lợi, Hà Nội). Đồ án tập trung giải quyết bài toán an ninh nguồn nước, cắt giảm lũ và cấp nước tưới tiêu bền vững cho vùng duyên hải Nam Trung Bộ.

                  LƯU VỰC SÔNG TRÀ BƯƠNG (F = 126 km²)
                                   │
                                   ▼
┌──────────────────────────────────────────────────────────────────────┐
│                  CỤM CÔNG TRÌNH ĐẦU MỐI PHƯỚC QUANG                  │
│                                                                      │
│   ┌─────────────────────┐  ┌──────────────────┐  ┌───────────────┐   │
│   │     ĐẬP ĐẤT CHÍNH   │  │   TRÀN XẢ LŨ     │  │ CỐNG LẤY NƯỚC │   │
│   │ L = 500 - 600 m     │  │ Tuyến yên ngựa   │  │ Đặt vai phải  │   │
│   │ Ngăn dòng chính     │  │ Ngưỡng tràn BTCT │  │ Không áp BTCT │   │
│   │ Thân đập đồng chất  │  │ Dốc nước tiêu    │  │ bc = 1.1 m    │   │
│   │ Chống thấm nền      │  │ năng sau tràn    │  │ Hc = 1.6 m    │   │
│   └──────────┬──────────┘  └────────┬─────────┘  └───────┬───────┘   │
└──────────────┼──────────────────────┼────────────────────┼───────────┘
               │                      │                    │
               │ (Chặn dòng & trữ)    │ (Tháo lũ thiết kế) │ (Qtk = 1.79 m³/s)
               ▼                      ▼                    ▼
        ┌──────────────┐       ┌──────────────┐     ┌──────────────┐
        │  LÒNG HỒ     │       │ HẠ LƯU SÔNG  │     │ KÊNH CHÍNH & │
        │  PHƯỚC QUANG │       │  TRÀ BƯƠNG   │     │ HỆ THỐNG N1, │
        │  W = 29.1    │       │  (Thoát lũ   │     │ N2 TƯỚI TỰ   │
        │  triệu m³    │       │   an toàn)   │     │ CHẢY 15000 ha│
        └──────────────┘       └──────────────┘     └──────────────┘

Bối cảnh và tính cấp thiết

Lưu vực sông Trà Bương (nhánh cấp 1 của sông Kỳ Lộ, tỉnh Phú Yên) có diện tích hứng nước $F_{lv} = 126\text{ km}^2$, chiều dài sông chính $L_s = 26{,}2\text{ km}$. Khu vực này chịu biến động khí hậu nhiệt đới gió mùa khốc liệt:

  • Mùa khô: Lượng mưa giảm sâu, lưu lượng dòng chảy kiệt chỉ đạt $Q_{85%} = 0{,}153\text{ m}^3/\text{s}$ (tháng 7-8), dẫn đến khô hạn nghiêm trọng trên 15.000 ha đất canh tác của hai xã Xuân Phước, Xuân Quang và thị trấn La Hai.
  • Mùa mưa: Lượng mưa năm đạt $1500 - 2000\text{ mm}$, thường xuyên hứng chịu các trận bão lớn với vận tốc gió cực trị $36\text{ m/s}$ (ghi nhận lịch sử $655\text{ mm}$/5 ngày), lưu lượng đỉnh lũ thiết kế $P = 1%$ lên tới $Q_{max} = 1032\text{ m}^3/\text{s}$ và lũ kiểm tra $P = 0{,}2%$ đạt $Q_{max} = 1313\text{ m}^3/\text{s}$, đe dọa trực tiếp đến tính mạng và cơ sở hạ tầng hạ du.

Tuyên bố bài toán (Problem Statement)

Khu vực hạ du sông Trà Bương đối mặt với 3 thách thức kỹ thuật lớn:

  1. Thiếu hụt nghiêm trọng công trình điều tiết dòng chảy năm ($W_0 = 77 \times 10^6\text{ m}^3$), gây lãng phí $85%$ lượng nước mặt đổ ra biển trong mùa lũ.
  2. Địa hình phân cách mạnh, nền địa chất thềm sông bồi tích á sét, á cát và cuội sỏi có tính thấm phức tạp, tiềm ẩn nguy cơ mất ổn định thấm và trượt cung trượt thân đập.
  3. Độ dốc thủy lực lớn đòi hỏi giải pháp tiêu năng đáy dốc nước và kết cấu tường chắn chịu lực ngang cực đại khi rút nước đột ngột.

Mục tiêu dự án

  1. Xác định thông số thủy văn - thủy năng: Tính toán điều tiết hồ chứa, xác định Mực nước chết ($MNC = +27{,}30\text{ m}$), Mực nước dâng bình thường ($MNDBT = +33{,}01\text{ m}$ đến $+35{,}18\text{ m}$), Mực nước lũ thiết kế ($MNLTK$) và lũ kiểm tra ($MNLKT$).
  2. Thiết kế kỹ thuật cụm công trình đầu mối Cấp II: Bao gồm đập đất đầm nén ngăn sông, tràn xả lũ thực dụng trên yên ngựa bờ tả và cống ngầm lấy nước không áp bờ hữu.
  3. Ứng dụng phần mềm chuyên ngành: Mô phỏng phần tử hữu hạn kết cấu tường bên dốc nước tràn bằng SAP2000 v14 và tính toán thấm qua đập bằng GeoSlope/GeoStudio.
  4. Lập dự toán và tối ưu hóa giải pháp kinh tế - kỹ thuật: Xác định chính xác khối lượng đào đắp, bê tông cốt thép M200, thép CII cho từng hạng mục công trình đầu mối.

Giải pháp và Phạm vi

  • Giải pháp lựa chọn (Phương án 1): Tuyến đập chính đắp bằng vật liệu địa phương (đất á sét đầm nén), tràn xả lũ đặt tại eo yên ngựa bờ tả tận dụng địa hình tự nhiên, cống lấy nước đặt ngầm qua thân đập ở vai phải dẫn vào kênh tưới tự chảy cao trình $Z_{kc} = +26{,}80\text{ m}$.
  • Phạm vi nghiên cứu: Thiết kế kỹ thuật cụm đầu mối công trình thủy lợi Cấp II theo quy chuẩn xây dựng Việt Nam QCVN 04-05:2012 và tiêu chuẩn TCVN 4118:2012, TCVN 9151:2012.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí so sánh Trạm bơm dã chiến / Đập dâng tạm Hồ chứa bê tông trọng lực Hồ chứa đập đất kết hợp tràn bên (Phương án 1)
Khả năng cấp nước Không ổn định ($Q < 0{,}5\text{ m}^3/\text{s}$), phụ thuộc mực nước kiệt Chủ động, dung tích lớn ($W > 25 \times 10^6\text{ m}^3$) Hoàn toàn chủ động theo mùa ($Q_{tk} = 1{,}79\text{ m}^3/\text{s}$, $W \approx 29{,}1 \times 10^6\text{ m}^3$)
Cắt lũ hạ du $0%$ khả năng cắt lũ Cắt lũ rất tốt qua tràn tự do/có van Cắt giảm đỉnh lũ từ $1313\text{ m}^3/\text{s}$ xuống ngưỡng an toàn
Vật liệu & Chi phí Chi phí ban đầu thấp, chi phí vận hành điện năng cao Chi phí vật liệu xi măng, cát đá rất cao Tận dụng $85%$ vật liệu đất đắp tại chỗ, kinh tế tối ưu
Địa chất thích ứng Không yêu cầu cao Yêu cầu nền đá cứng cường độ cao Phù hợp nền đất đệ Tứ, cuội sỏi, á sét có xử lý chống thấm
Độ bền công trình $1 - 3$ năm (dễ bị lũ cuốn trôi) $> 100$ năm $> 70$ năm (Công trình cấp II vĩnh cửu)

Phân tích yêu cầu theo mô hình MoSCoW

  • Must have (Bắt buộc): Đập ngăn sông an toàn thấm - trượt; Tràn tháo được lưu lượng kiểm tra $Q_{0{,}2%} = 1313\text{ m}^3/\text{s}$; Cống xả lưu lượng tưới $Q_{tk} = 1{,}79\text{ m}^3/\text{s}$ ở cao trình mực nước chết $+27{,}30\text{ m}$.
  • Should have (Nên có): Tháp van điều khiển thượng lưu bố trí van công tác và van sửa chữa; Thiết bị tiêu năng đáy dốc nước sau tràn hạn chế xói lở.
  • Could have (Có thể mở rộng): Tích hợp trạm thủy điện nhỏ chân đập tận dụng lưu lượng xả thừa mùa lũ.
  • Won't have (Loại trừ): Hệ thống cửa van tự động SCADA thời gian thực trong giai đoạn thiết kế kỹ thuật cơ sở ban đầu.

Thiết kế hệ thống và Cấu trúc kỹ thuật

Cụm công trình đầu mối được bố trí đồng bộ theo mặt bằng tổng thể:

  1. Đập chính: Chiều dài tuyến $L = 500 - 600\text{ m}$, đập đất không đồng chất hoặc đồng chất có chân khay chống thấm cắm sâu vào tầng đá gốc/lớp địa chất ít thấm, thiết bị thoát nước kiểu ống khói kết hợp lăng trụ hạ lưu.
  2. Tràn xả lũ: Tràn thực dụng đỉnh rộng không chân van đặt tại eo bờ tả, dốc nước bê tông cốt thép M200, bể tiêu năng sâu đáy gia cố tấm đan BTCT dày.
  3. Cống ngầm lấy nước: Cống tiết diện chữ nhật $b_c \times H_c = 1{,}1\text{ m} \times 1{,}6\text{ m}$, chiều dài toàn bộ $L = 64{,}54\text{ m}$ ($L_1 = 15\text{ m}$ từ cửa vào đến tháp van; $L_2 = 49{,}54\text{ m}$ từ tháp van đến cửa ra), độ dốc đáy cống $i = 0{,}00574$.

Technology Stack và Công cụ mô phỏng

  • Phân tích kết cấu (FEM): CSI SAP2000 v14.0 (Finite Element Method - Phần tử Shell kết hợp lò xo nền đàn hồi Winkler Area Springs).
  • Mô phỏng trường thấm: GeoStudio GeoSlope 2018 (Phân tích dòng thấm bão hòa và không bão hòa qua thân và nền đập).
  • Thiết kế hình học & Bản vẽ: AutoCAD Civil 3D 2022.
  • Tính toán thủy lực & Thủy năng: Numerical solver phát triển trên nền Python 3.10 / Excel VBA.

Cơ sở tiêu chuẩn kỹ thuật áp dụng

  • Quy chuẩn kỹ thuật quốc gia về công trình thủy lợi: QCVN 04-05:2012/BNNPTNT (Cấp công trình: Cấp II, hệ số tầm quan trọng $k_n = 1{,}15$).
  • Tiêu chuẩn thiết kế kênh và công trình trên kênh: TCVN 4118:2012.
  • Tiêu chuẩn tính toán thủy lực cống ngầm, cửa van: TCVN 9151:2012.
  • Tiêu chuẩn thiết kế kết cấu bê tông cốt thép thủy công: TCVN 5574:2018.

Phương pháp luận (Methodology)

Quy trình thiết kế tuân thủ mô hình tuyến tính nghiêm ngặt của kỹ thuật công trình (Waterfall Engineering Framework) kéo dài trong 14 tuần:

Tuần 1: Khảo sát địa hình, địa chất, khí tượng thủy văn (5%)
  └── Tuần 2-3: Tính toán điều tiết hồ chứa & xác định các mực nước đặc trưng (15%)
        └── Tuần 4-9: Thiết kế kỹ thuật Đập đất, Tràn xả lũ, Cống ngầm (45%)
              └── Tuần 10-14: Chuyên đề SAP2000, Tính nứt kết cấu & Lập dự toán (35%)

Ma trận đánh giá rủi ro và Biện pháp giảm thiểu

  • Rủi ro bồi lắng bùn cát: Nồng độ bùn cát $\rho_0 = 90\text{ g/m}^3$, dung trọng bồi lắng $\gamma_{bc} = 0{,}9\text{ T/m}^3$. Bố trí cao trình đáy cống $Z_{cv} = +26{,}25\text{ m}$ cao hơn cao trình bùn cát $Z_{bc} = +25{,}53\text{ m}$ một khoảng an toàn $> 0{,}7\text{ m}$.
  • Rủi ro xói lở kênh dẫn hạ lưu: Giới hạn lưu tốc thực tế $V_{max} = 0{,}49\text{ m/s}$ luôn nhỏ hơn lưu tốc không xói cho phép $[V_{kx}] = 0{,}66\text{ m/s}$ của nền đất á sét nhẹ.
  • Rủi ro mất ổn định lật và trượt tường ngực/tường bên tràn: Tính toán kiểm tra với 2 tổ hợp tải trọng bất lợi nhất: tải trọng thi công và tải trọng khi xả lũ xong rút nước đột ngột.

Implementation và kết quả

Quy trình phát triển và Thuật toán cốt lõi

1. Thuật toán xác định khẩu diện và tổn thất cột nước cống ngầm

Để đảm bảo tháo được lưu lượng thiết kế $Q_{tk} = 1{,}79\text{ m}^3/\text{s}$ ở mực nước bất lợi nhất ($MNC = +27{,}30\text{ m}$, mực nước đầu kênh $Z_{kc} = +26{,}80\text{ m}$, chênh lệch cho phép $[\Delta Z] \le 0{,}50\text{ m}$), tổng tổn thất qua cống được tính toán bằng phương pháp thử dần:

$$\sum Z_i = Z_1 + Z_p + Z_L + Z_v + Z_2 + J \cdot L_c \le [\Delta Z]$$

Trong đó:

  • Tổn thất cửa vào: $Z_1 = \frac{Q^2}{2g \varphi^2 \omega^2} - \frac{\alpha V_0^2}{2g}$
  • Tổn thất lưới chắn rác: $Z_L = \xi_L \frac{V_L^2}{2g}$ với $\xi_L = \beta \left(\frac{s}{b}\right)^{4/3} \sin\alpha = 1{,}67 \left(\frac{10}{50}\right)^{4/3} \sin(60^\circ) = 0{,}169$
  • Tổn thất khe phai và van: $Z_p = 2 \xi_p \frac{V_p^2}{2g}$; $Z_v = 2 \xi_v \frac{V_v^2}{2g}$
  • Tổn thất dọc đường: $J \cdot L_c = \frac{Q^2}{\omega^2 C^2 R} L_c$, với hệ số Chezy $C = \frac{1}{n} R^{1/6}$, nhám bê tông $n = 0{,}017$.

Sau khi tối ưu hóa giải tích, kích thước cống được chọn là $b_c = 1{,}10\text{ m}$, $H_c = 1{,}60\text{ m}$, đạt tổng tổn thất $\sum Z_i = 0{,}50\text{ m} = [\Delta Z]$.

import numpy as np

def calculate_culvert_backwater(Q=1.79, b=1.1, n=0.017, i=0.00301, g=9.81):
    """
    Tính toán đường mặt nước dâng C1 trong cống ngầm bằng phương pháp cộng trực tiếp.
    """
    h_k = (Q**2 / (g * b**2))**(1/3) # Độ sâu phân giới: 0.646 m
    # Giải tìm độ sâu dòng đều h0 theo Chezy-Manning
    h_0 = 1.01  # m (nghiệm phương trình lưu lượng)
    h_c = 0.14053  # Độ sâu co hẹp sau cửa van (m)
    
    # Khởi tạo bước tính
    h_steps = np.linspace(h_c, h_k, 50)
    L_cum = 0.0
    trajectory = [(0.0, h_c)]
    
    for idx in range(len(h_steps) - 1):
        h1, h2 = h_steps[idx], h_steps[idx+1]
        A1, A2 = b * h1, b * h2
        v1, v2 = Q / A1, Q / A2
        E1 = h1 + (v1**2) / (2 * g)
        E2 = h2 + (v2**2) / (2 * g)
        dE = E2 - E1
        
        R1 = A1 / (b + 2 * h1)
        R2 = A2 / (b + 2 * h2)
        C1 = (1 / n) * (R1**(1/6))
        C2 = (1 / n) * (R2**(1/6))
        
        J1 = (v1**2) / ((C1**2) * R1)
        J2 = (v2**2) / ((C2**2) * R2)
        J_avg = (J1 + J2) / 2
        
        # Khoảng cách giữa 2 mặt cắt
        dL = dE / (i - J_avg)
        L_cum += dL
        trajectory.append((L_cum, h2))
        
    return h_k, h_0, L_cum, trajectory

hk, h0, L_total, curve = calculate_culvert_backwater()
print(f"hk = {hk:.3f}m | h0 = {h0:.3f}m | Khoang cach den nuoc nhay = {L_total:.2f}m")

2. Phân tích nước nhảy trong cống

  • Kết quả chạy thuật toán cho thấy tại vị trí $L = 45{,}90\text{ m} < L_{tt} = 49{,}218\text{ m}$ (tính từ sau van), đường mặt nước dâng $C_I$ cắt đường phân giới $K-K$ ($h_r > h_k$).
  • Hiện tượng nước nhảy sinh ra ngay trong thân cống. Chiều cao nước nhảy liên hiệp $h_c'' = \frac{h_c}{2} \left(\sqrt{1 + 8 \frac{h_k^3}{h_c^3}} - 1\right)$ được tính toán chứng minh không chạm trần cống ($H_c = 1{,}6\text{ m} > h_c''$), bảo đảm an toàn khí động, không phát sinh chân không gây rung giật nguy hiểm cho khối bê tông.

Kiểm chứng và Đánh giá kết quả (Testing & Validation)

Kiểm tra ổn định tường chắn bên dốc nước tràn (Tổ hợp xả lũ xong, nước rút đột ngột)

Sử dụng các thông số địa chất: $\gamma_{tn} = 1{,}54\text{ T/m}^3$, $\varphi_{bh} = 18^\circ$, lực dính $C_{bh} = 1{,}6\text{ T/m}^2$, nền $C = 2{,}6\text{ T/m}^2$, hệ số ma sát $f = 0{,}65$.

                            SƠ ĐỒ TẢI TRỌNG TÁC DỤNG LÊN TƯỜNG BÊN
                                  │◄─── B = 6.5 m ───►│
                                  ┌──────────────────┐
                                  │      TƯỜNG       │▲
                       Áp lực đất │        P1        ││
                         chủ động │   ┌──────────────┤│
                          (Ecđ)   │   │              ││
                         ══════>  │   │  TL đất (P4) ││ H = 7.97 m
                         ══════>  │   │              ││
                         ══════>  │   │  TL nước(P5) ││
                                  ├───┴──────────────┤│
                       Áp lực     │  BẢN ĐÁY (P2,P3) │▼
                       thấm(Wth)  └───┬──────────┬───┘
                         ▲ ▲ ▲ ▲      │  Điểm A  │
                         └─┴─┴─┴──────┴──────────┴────────

  • Tổng hợp tải trọng:

    • Tổng lực thẳng đứng tính toán: $\sum G = 75{,}129\text{ T}$
    • Tổng lực đẩy ngang tính toán: $\sum P = 27{,}173\text{ T}$
    • Tổng mômen chống lật đối với tâm quay $A$: $M_{CL} = 210{,}048\text{ T}\cdot\text{m}$
    • Tổng mômen gây lật đối với tâm quay $A$: $M_{GL} = 75{,}188\text{ T}\cdot\text{m}$

  • Chỉ số an toàn trượt và lật:

$$K_L = \frac{M_{CL}}{M_{GL}} = \frac{210{,}048}{75{,}188} = 2{,}79 \ge [K_L] = 1{,}15 \quad (\text{Thỏa mãn điều kiện chống lật})$$

$$K_T = \frac{f \cdot \sum G + C \cdot B}{\sum P} = \frac{0{,}65 \times 75{,}129 + 2{,}6 \times 6{,}5}{27{,}173} = 2{,}42 \ge [K_T] = 1{,}15 \quad (\text{Thỏa mãn điều kiện chống trượt})$$

Phân tích nội lực phần tử hữu hạn bằng SAP2000 v14

Mô hình hóa kết cấu tường chắn dạng chữ L có sườn tăng cường bằng phần tử tấm vỏ 3D (Shell Elements), gán hệ số đàn hồi nền đất $K_0 = 1{,}4 \times 10^5\text{ kN/m}^3$:

  • Bản mặt tường:
    • Mômen uốn ngang $M_{11}$: $M_{max}^+ = +15{,}02\text{ kN}\cdot\text{m}$, $M_{max}^- = -1{,}41\text{ kN}\cdot\text{m}$.
    • Mômen uốn đứng $M_{22}$: $M_{max}^+ = +8{,}10\text{ kN}\cdot\text{m}$, $M_{max}^- = -71{,}37\text{ kN}\cdot\text{m}$ (tại vị trí ngàm với bản đáy).
  • Bản đáy tường chắn:
    • Mômen uốn ngang $M_{11}$: $M_{max}^+ = +237{,}57\text{ kN}\cdot\text{m}$ (tại tiết diện $y = 2{,}5\text{ m}$), $M_{max}^- = -29{,}84\text{ kN}\cdot\text{m}$.
    • Kiểm tra bề rộng khe nứt bê tông theo TCVN 5574:2018: $a_{crc} \le 0{,}2\text{ mm}$, kết cấu bảo đảm bền và chống thấm tuyệt đối.

Đổi mới và đóng góp

  1. Ứng dụng mô hình tương tác đất - kết cấu (Soil-Structure Interaction) 3D: Thay vì tính toán theo sơ đồ dầm dải 2D độc lập truyền thống, đồ án ứng dụng SAP2000 v14 mô hình hóa tường sườn và bản đáy trên nền đàn hồi Winkler ($K_0 = 1{,}4 \times 10^5\text{ kN/m}^3$). Giải pháp này giúp phân bố lại nội lực chính xác, giảm $18{,}5%$ khối lượng cốt thép chịu uốn tại nách liên kết so với phương pháp cổ điển.
  2. Định lượng không gian nước nhảy trong cống ngầm: Thiết lập thuật toán tích phân từng bước vẽ đường mặt nước $C_I$ và $b_I$, xác định chính xác tọa độ va đập thủy lực ($L = 45{,}90\text{ m}$), ngăn chặn triệt để nguy cơ cộng hưởng và tạo túi khí xâm thực thân cống ngầm dưới thân đập cấp II.
  3. Tối ưu hóa khẩu diện thủy lực kênh và cống: Tận dụng tối đa chênh lệch cột nước hạn chế ($\Delta Z = 0{,}50\text{ m}$), thiết kế tiết diện kênh hạ lưu $b_k = 2{,}30\text{ m}$, $h_h = 1{,}07\text{ m}$ đạt trạng thái thủy lực lợi nhất, khống chế lưu tốc xói lở $V_{max} = 0{,}49\text{ m/s} < 0{,}66\text{ m/s}$.
      PHƯƠNG PHÁP CỔ ĐIỂN (Dầm dải 2D)            MÔ HÌNH PTHH SAP2000 3D (Shell + Spring)
┌─────────────────────────────────────────┐     ┌─────────────────────────────────────────┐
│ Giả thiết ngàm cứng lý tưởng đáy móng    │     │ Xét biến dạng đàn hồi nền (K₀ = 1.4e5)  │
│ Thiên về quá an toàn, cốt thép dư thừa  │ VS  │ Tối ưu hóa biểu đồ mômen uốn            │
│ Không xét đến tương tác sườn - bản mặt  │     │ Tiết kiệm 18.5% lượng cốt thép chịu lực │
└─────────────────────────────────────────┘     └─────────────────────────────────────────┘

Ứng dụng thực tế và triển khai

Tình huống ứng dụng thực tế

  • Cấp nước tưới ổn định: Cung cấp nguồn nước tự chảy liên tục $Q_{tk} = 1{,}79\text{ m}^3/\text{s}$ cho 15.000 ha lúa, hoa màu và cây công nghiệp hai bên bờ sông Trà Bương, nâng hệ số quay vòng sử dụng đất từ 1 vụ lên 2-3 vụ/năm.
  • Phòng chống thiên tai: Cắt đỉnh lũ $P = 1%$ từ $1032\text{ m}^3/\text{s}$ và xả an toàn qua tràn dốc nước, bảo vệ toàn bộ dân sinh hạ lưu thị trấn La Hai khỏi ngập lụt cục bộ.

Lộ trình thi công và Hiệu quả đầu tư

  • Giai đoạn 1 (Tháng 1 - 3): Dẫn dòng thi công mùa khô qua cống ngầm kết hợp đê quai giai đoạn 1.
  • Giai đoạn 2 (Tháng 4 - 9): Đắp đập đất chính đến cao trình vượt lũ an toàn; thi công tràn xả lũ bê tông cốt thép M200 bờ tả.
  • Giai đoạn 3 (Tháng 10 - 12): Hoàn thiện tháp van, lắp đặt van công tác, nạo vét kênh chính và bàn giao vận hành.
  • Hiệu quả kinh tế: Dự án hoàn vốn đầu tư xã hội ước tính trong $6 - 8$ năm nhờ gia tăng giá trị sản xuất nông nghiệp và loại bỏ thiệt hại do bão lũ thường niên gây ra.

Hạn chế và hướng phát triển

  • Hạn chế: Mô hình SAP2000 sử dụng mô hình nền Winkler đơn giản hóa (lò xo đàn hồi tuyến tính), chưa mô phỏng phi tuyến đàn dẻo của đất nền khi chịu biến dạng trượt lớn dưới tác động của động đất.
  • Hướng phát triển:
    1. Ứng dụng phần mềm Plaxis 3D hoặc FLAC3D để kiểm tra tương tác động lực học (Dynamic Soil-Structure Interaction) khi có gia tốc động đất thiết kế cấp VII.
    2. Ứng dụng mô phỏng thủy lực 3D OpenFOAM / Flow-3D kiểm tra hiện tượng sóng đứng và xoáy cuộn trong bể tiêu năng tràn xả lũ.
    3. Lắp đặt hệ thống cảm biến đo áp lực nước lỗ rỗng (Piezometer) tự động truyền dữ liệu IoT về trung tâm giám sát an toàn đập.

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

                              CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI
                                            │
        ┌───────────────────┬───────────────┴───────────────┬───────────────────┐
        ▼                   ▼                               ▼                   ▼
┌───────────────┐   ┌───────────────┐               ┌───────────────┐   ┌───────────────┐
│  SINH VIÊN    │   │ KỸ SƯ THỦY LỢI│               │ CƠ QUAN QUẢN  │   │  NHÀ NGHIÊN   │
│   NGÀNH XÂY   │   │  & KẾT CẤU    │               │  LÝ & ĐỊA PHƯƠNG│ │     CỨU       │
│     DỰNG      │   │               │               │               │   │               │
│ Tài liệu mẫu  │   │ Quy trình mô  │               │ Cơ sở phê     │   │ Dữ liệu thủy  │
│ đồ án tốt     │   │ hình SAP2000  │               │ duyệt báo cáo │   │ văn, giải     │
│ nghiệp Thủy   │   │ & giải thuật  │               │ kỹ thuật kinh │   │ thuật thủy    │
│ công chuẩn    │   │ thủy lực cống │               │ tế hồ chứa    │   │ lực cống ngầm │
└───────────────┘   └───────────────┘               └───────────────┘   └───────────────┘
  • Sinh viên ngành Kỹ thuật Công trình Thủy / Xây dựng: Cung cấp tài liệu tham khảo hoàn chỉnh từ lý thuyết điều tiết thủy năng, thiết kế đập - tràn - cống đến mô hình hóa phần mềm.
  • Kỹ sư tư vấn thiết kế: Nắm bắt quy trình phân tích nội lực tường chắn bằng phần tử Shell trên nền đàn hồi và giải thuật kiểm soát nước nhảy trong cống ngầm.
  • Ban quản lý dự án & Chính quyền địa phương: Cơ sở khoa học vững chắc để thẩm định phương án kỹ thuật hồ Phước Quang, tối ưu ngân sách nhà nước.
  • Nhà nghiên cứu thủy văn - thủy lực: Dữ liệu chuẩn xác về đặc trưng thủy văn sông Kỳ Lộ và ứng dụng toán giải tích trong thủy lực công trình ngầm.

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

1. Yêu cầu kỹ thuật chính để triển khai cống ngầm dưới đập đất là gì?

Cống phải đặt trên nền đá gốc hoặc nền đất được xử lý gia cố; sử dụng bê tông cốt thép M200 trở lên. Mặt bằng cống đặt tại vai đập để giảm áp lực đất; bố trí các vòng cắt nước (khớp nối chống thấm) dọc thân cống để ngăn dòng thấm tập trung men theo thành cống.

2. Làm thế nào để giải quyết hiện tượng nước nhảy xuất hiện trong thân cống?

Kiểm tra chiều cao nước nhảy $h_c''$. Nếu $h_c'' < H_c$ (không chạm trần), chỉ cần gia cố bề mặt bê tông chống xâm thực và xử lý kín khít khớp nối. Nếu nước nhảy chạm trần tạo bọt khí gây rung giật, cần tăng độ dốc đáy cống $i$, dịch chuyển vị trí tháp van hoặc tăng chiều cao cống $H_c$.

3. Tại sao chọn phần tử Shell thay vì Solid khi mô hình tường chắn trong SAP2000?

Bản mặt và bản đáy tường chắn có chiều dày nhỏ hơn nhiều so với kích thước mặt bằng ($t \ll L, H$), phần tử Shell tích hợp cả trạng thái ứng suất phẳng trong mặt phẳng (Membrane) và uốn ngoài mặt phẳng (Plate Bending), cho kết quả mômen uốn và ứng suất chính xác mà không tốn tài nguyên tính toán như Solid.

4. Hệ số an toàn ổn định trượt và lật của công trình đạt bao nhiêu?

Theo QCVN 04-05:2012 đối với công trình Cấp II, hệ số cho phép $[K] = 1{,}15$. Tính toán thực tế đạt $K_L = 2{,}79$ (vượt $142{,}6%$) và $K_T = 2{,}42$ (vượt $110{,}4%$), bảo đảm an toàn tuyệt đối dưới mọi tổ hợp tải trọng.

5. Chi phí và thời gian thu hồi vốn của dự án như thế nào?

Dự án tận dụng đất đắp tại chỗ làm đập chính, giúp tiết kiệm $30 - 40%$ chi phí xây lắp so với đập bê tông. Nhờ cung cấp nước tưới cho 15.000 ha và bảo vệ an toàn sản xuất, thời gian thu hồi vốn đầu tư xã hội ước tính đạt $6 - 8$ năm.

Kết luận

Đồ án tốt nghiệp "Thiết kế hồ chứa nước Phước Quang – Phương án 1" của tác giả Nguyễn Quốc Khánh (Khoa Công Trình, Trường Đại học Thủy Lợi) là một công trình nghiên cứu ứng dụng mẫu mực, kết hợp nhuần nhuyễn giữa lý thuyết cơ sở công trình thủy công và công nghệ mô phỏng phần tử hữu hạn hiện đại (SAP2000, GeoSlope).

Công trình không chỉ đưa ra giải pháp toàn diện, tối ưu về mặt kết cấu và thủy lực cho cụm đầu mối Cấp II, mà còn mang ý nghĩa kinh tế - xã hội to lớn, giải quyết tận gốc bài toán thiếu nước tưới và phòng chống lũ lụt cho lưu vực sông Trà Bương, huyện Đồng Xuân, tỉnh Phú Yên. Đây là tài liệu kỹ thuật giá trị cao dành cho các kỹ sư, nhà nghiên cứu và sinh viên chuyên ngành Công trình Thủy.