Tổng quan về luận án

Bối cảnh khoa học của nghiên cứu xuất phát từ thực tiễn gia tăng áp lực nhu cầu sử dụng nước và các biến động bất lợi của biến đổi khí hậu trên Lưu vực sông Hồng – Thái Bình (LVSHTB) – một lưu vực quốc tế có tổng diện tích 169.000 $km^2$ trải dài qua Việt Nam, Trung Quốc và Lào. Tại Việt Nam, vùng châu thổ sông Hồng tập trung 29 triệu dân (chiếm 95% dân số miền Bắc) và đóng góp gần 20% GDP quốc gia. Tổng lượng dòng chảy hàng năm tại trạm thủy văn Sơn Tây đạt khoảng 118 tỷ $m^3$, và khi tính toàn bộ châu thổ có thể lên tới 135 tỷ $m^3$ (trong đó 82,54 tỷ $m^3$, tương đương 61,1% hình thành nội địa và 52,46 tỷ $m^3$, tức 38,9% từ lãnh thổ Trung Quốc). Mặc dù hệ thống hạ tầng thủy lợi và thủy điện được đầu tư quy mô lớn với 4 hồ chứa đa mục tiêu chiến lược gồm Thác Bà (440 MW), Hòa Bình (1.920 MW), Tuyên Quang (750 MW) và Sơn La (2.400 MW) sở hữu tổng dung tích phòng lũ 8,45 tỷ $m^3$, dung tích hữu ích 19,86 tỷ $m^3$ và tổng công suất lắp đặt 6.246 MW, công tác quy hoạch truyền thống vẫn bộc lộ khoảng trống học thuật nghiêm trọng.

Khoảng trống nghiên cứu then chốt (research gap) nằm ở sự chia cắt giữa mô hình hóa thủy lực đơn thuần và phân tích kinh tế vi mô có tính đến bản chất động học và ngẫu nhiên. Các nghiên cứu trước đây (như Viện Quy hoạch Thủy lợi, 2005, 2007; Lê Kim Truyền, 2005) phần lớn sử dụng mô hình tối ưu tĩnh hoặc tách rời quy trình vận hành mùa lũ và mùa kiệt, thiếu vắng việc lượng hóa hàm cầu biên và hàm giá trị kinh tế của các đối tượng hưởng lợi.

Luận án "Áp dụng lý thuyết tối ưu hóa cho bài toán phân bổ hiệu quả tài nguyên nước ở lưu vực sông Hồng - Thái Bình" do NCS. Đào Văn Khiêm thực hiện dưới sự hướng dẫn khoa học của GS.TS. Nguyễn Quang Kim tại Trường Đại học Thủy lợi (2017/2018) đã giải quyết triệt để bài toán này. Nghiên cứu đặt ra hệ thống câu hỏi nghiên cứu (RQ) và giả thuyết khoa học (H):

  • RQ1: Làm thế nào để lượng hóa hàm giá trị kinh tế của các ngành dùng nước cạnh tranh (tưới nông nghiệp, thủy điện, sinh hoạt, du lịch) dưới dạng hàm mục tiêu tích hợp?
  • RQ2: Mô hình tối ưu hóa động ngẫu nhiên giải quyết sự bất định của dòng chảy đến và biến động giá cả nông sản như thế nào so với tiếp cận tất định truyền thống?
  • RQ3: Cấu trúc chi phí và lợi ích của hệ thống hồ chứa lớn có biểu hiện thuộc tính độc quyền tự nhiên (natural monopoly) hay không, và cơ chế điều tiết của Nhà nước cần can thiệp ở mức độ nào?
  • H1: Việc tích hợp chuỗi thời gian ngẫu nhiên (ARIMA/SARIMA) của dòng chảy và giá cả vào mô hình tối ưu hóa động sẽ làm thay đổi đáng kể quỹ đạo phân bổ nước tối ưu, nâng cao tổng thặng dư xã hội so với quy trình xả nước cứng nhắc.
  • H2: Hệ thống hồ chứa đa mục tiêu tại LVSHTB mang đặc tính độc quyền tự nhiên yếu (weak natural monopoly), đòi hỏi sự can thiệp và điều tiết công quyền thay vì để thị trường tự do vận hành thuần túy.

Khung lý thuyết của luận án được xây dựng dựa trên sự hợp nhất giữa Lý thuyết Phân bổ Tài nguyên Tối ưu Pareto (James & Lee, 1971), Kinh tế học Tài nguyên Nước Hiện đại (Loucks & Van Beek, 2017; McKinney et al., 1999) và Lý thuyết Tổ chức Công nghiệp về Độc quyền Tự nhiên (Baumol, Panzar, & Willig, 1982). Phạm vi thực nghiệm bao gồm chuỗi số liệu thủy văn và kinh tế 21 năm đối với hệ thống hồ Núi Cốc trên lưu vực sông Công (1980–2000) và 25 năm đối với hệ thống liên hồ chứa Sơn La – Hòa Bình trên sông Đà (1980–2004). Ý nghĩa đột phá của luận án là việc xây dựng thành công bộ công cụ mô hình hóa kinh tế thủy văn tổng hợp (Hydro-economic Modeling) giải bằng phần mềm Lingo, cung cấp luận cứ khoa học trực tiếp cho việc triển khai Luật Thủy lợi năm 2017 và điều chỉnh Quy trình vận hành liên hồ chứa theo Quyết định số 1287/QĐ-TTg của Thủ tướng Chính phủ.

Literature Review và Positioning

Tổng quan tài liệu quốc tế cho thấy sự tiến hóa rõ rệt của các trường phái quản lý tài nguyên nước cấp lưu vực. Theo báo cáo kinh điển của McKinney, Cai, Rosegrant, Ringler, và Scott (1999) cho Viện Quản lý Nước Quốc tế (IWMI), bản chất liên ngành của tài nguyên nước đòi hỏi sự tích hợp giữa các mô hình thủy văn - thủy lực và mô hình kinh tế phúc lợi xã hội:

"Bản chất liên ngành của những bài toán ngành nước yêu cầu những phương pháp mới để tổng hợp các khía cạnh kỹ thuật, kinh tế, môi trường, xã hội và pháp luật vào một khung mạch lạc."

Rosegrant và Meinzen-Dick (1996) nhấn mạnh nguyên lý nền tảng:

"Chỉ bằng cách nghiên cứu tất cả các cấu phần tương tác hưởng lợi từ và gây hại cho nguồn tài nguyên, sử dụng tối ưu từ quan điểm xã hội mới có thể được thiết lập."

   Trường phái Kỹ thuật Thủy văn        Trường phái Kinh tế Tài nguyên Nước
   (Loucks & Van Beek, 2017; MIKE)     (James & Lee, 1971; Booker & Colby, 1995)

Trong dòng chảy học thuật quốc tế, có hai cách tiếp cận chính để xây dựng mô hình kinh tế thủy văn tổng hợp:

  1. Tiếp cận mô hình hóa cục bộ (Compartmental Approach): Kết nối lỏng lẻo giữa mô hình thủy văn và mô hình kinh tế bằng cách truyền dữ liệu đầu ra qua lại (như Howe & Orr, 1974; Noel & Howitt, 1982; Lefkoff & Gorelick, 1990). Cách tiếp cận này duy trì độ chi tiết thủy lực nhưng gặp trở ngại lớn trong trao đổi thông tin lặp và hội tụ thuật toán.
  2. Tiếp cận mô hình hóa hệ thống (Holistic/System Approach): Hợp nhất toàn bộ phương trình cân bằng dòng chảy, quan hệ điều tiết hồ chứa và hàm giá trị kinh tế vào một hệ thống tối ưu hóa đơn nhất (Booker & Colby, 1995; Diaz, Brown, & Sveinsson với mô hình Aquarius).

Các tranh luận khoa học lớn trong tài liệu xoay quanh hai vấn đề:

  • Định giá nước nông nghiệp: Sự đối lập giữa mô hình sinh trưởng cây trồng - bốc thoát hơi nước (Crop Water Production Functions - Dinar & Letey, 1996; Hanks & Hill, 1980) và phương pháp quy hoạch toán học dẫn xuất hàm cầu nước tưới (Mathematical Programming - Griffin, 2006; Young, 2005; Bowen & Young, 1985).
  • Cơ chế điều tiết thị trường: Tranh luận giữa trường phái kinh tế thị trường tự do xem nước là hàng hóa tư nhân thuần túy và trường phái kinh tế thể chế nhấn mạnh các "thất bại thị trường" do ngoại ứng tiêu cực, hàng hóa công cộng và tính chất độc quyền tự nhiên của hạ tầng thủy công (Dublin Statement, 1992).

Tại Việt Nam, các nghiên cứu trước đây tạo ra những tiền đề quan trọng nhưng còn nhiều hạn chế. Dự án hợp tác giữa Đại học Tổng hợp Milan và Viện Quy hoạch Thủy lợi (2012–2016) đoạt giải thưởng EEARRI 2015 đã bước đầu tiếp cận quản lý lưu vực sông Hồng bền vững nhưng tập trung chủ yếu vào hỗ trợ ra quyết định đa mục tiêu. Nghiên cứu của Claudia Ringler và Nguyễn Vũ Huy (2004, 2007) trên lưu vực sông Đồng Nai đã ứng dụng hàm cầu kinh tế nhưng chỉ dừng lại ở bài toán tối ưu hóa tĩnh. Nghiên cứu của Viện Quy hoạch Thủy lợi (2005) sử dụng phần mềm GAMS nhưng chưa nội sinh hóa hàm cầu biên của người dùng nước.

Luận án của Đào Văn Khiêm định vị chính xác vào khoảng trống khoa học này: vượt qua giới hạn của mô hình Aquarius (vốn chỉ xét hàm cầu tĩnh, bỏ qua cấu trúc hàm chi phí hoàn chỉnh và chưa xét quá trình ngẫu nhiên của dòng chảy), luận án đã phát triển mô hình kinh tế thủy văn động ngẫu nhiên toàn diện, kết hợp kiểm định cấu trúc thị trường độc quyền tự nhiên. Nghiên cứu so sánh trực tiếp với hai điển cứu quốc tế lớn: mô hình phân bổ lưu vực sông Colorado của Booker và Colby (1995) và mô hình tối ưu hóa lưu vực sông ngòi Hoa Kỳ của Viện Nghiên cứu Rừng Rocky Mountain (USDA).

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

Đóng góp cho lý thuyết

Luận án đóng góp vào việc mở rộng và tái định hình các lý thuyết kinh tế tài nguyên nước trong điều kiện các nền kinh tế chuyển đổi có sự điều tiết của Nhà nước:

   Lý thuyết Phân bổ Tối ưu Pareto          Lý thuyết Tổ chức Công nghiệp
   (James & Lee, 1971; Douglas, 1971)       (Baumol, Panzar, & Willig, 1982)
  1. Mở rộng lý thuyết phân bổ Pareto của James Douglas (1971): Luận án thiết lập bài toán phúc lợi xã hội mở rộng $W$, tối đa hóa tổng thặng dư kinh tế ròng từ các đối tượng dùng nước phụ thuộc vào hệ ràng buộc cân bằng nước động: $$\max W = \sum_{t=1}^{T} \frac{1}{(1+\rho)^t} \left[ \sum_{i} \int_{0}^{x_{it}} P_i(u) du - C_t(S_t, R_t, Q_t) \right]$$ Trong đó:

    • $P_i(u)$ là hàm cầu nghịch đảo của ngành dùng nước $i$ (tưới, phát điện, cấp nước, du lịch).
    • $C_t(S_t, R_t, Q_t)$ là hàm chi phí vận hành, lưu trữ và xả nước của hệ thống tại chu kỳ $t$.
    • $\rho$ là tỷ suất chiết khấu kinh tế.
    • Điều kiện biên bậc nhất xác lập mối quan hệ: giá bóng tài nguyên (shadow price - nhân tử Lagrange $\lambda_t$) đo lường sự gia tăng phúc lợi xã hội biên khi dung tích hữu dụng tại thời điểm $t$ tăng thêm một đơn vị.
  2. Ứng dụng Lý thuyết Độc quyền Tự nhiên (Natural Monopoly Theory) của Baumol, Panzar, và Willig (1982): Luận án là công trình tiên phong tại Việt Nam và trên thế giới áp dụng khái niệm tính cộng dưới của hàm chi phí (Cost Subadditivity) vào hệ thống hồ chứa liên hoàn: $$C(q_1 + q_2) < C(q_1) + C(q_2)$$ Nghiên cứu chứng minh rằng do chi phí cố định xây dựng đập và hồ chứa cực kỳ lớn trong khi chi phí biên vận hành rất nhỏ, một hệ thống điều hành tập trung (như liên hồ Sơn La – Hòa Bình do Nhà nước quản lý thông qua A0/EVN và MARD) đạt hiệu quả kinh tế quy mô vượt trội so với cơ chế phân tán hoặc thị trường cạnh tranh tự do hoàn hảo.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu phần lý thuyết:

  • Lý thuyết Kinh tế Phúc lợi: Sử dụng khái niệm Thặng dư Người Tiêu dùng (Consumer Surplus - CS) và Thặng dư Nhà Sản xuất (Producer Surplus - PS) để xây dựng hàm mục tiêu tổng thể.
  • Lý thuyết Tối ưu hóa Động Ngẫu nhiên: Rời rạc hóa thời gian theo tháng ($t = 1, \dots, 12$ và mở rộng đến 25 năm) với các phương trình trạng thái bảo toàn khối lượng liên tục: $$S_{t+1} = S_t + I_t - R_t - E_t(S_t) - L_t$$ (với $S_t$: dung tích trữ, $I_t$: dòng chảy đến ngẫu nhiên, $R_t$: lượng nước xả sử dụng, $E_t$: tổn thất bốc hơi phụ thuộc diện tích mặt nước, $L_t$: thấm).
  • Lý thuyết Chuỗi thời gian Box-Jenkins: Mô hình hóa tính ngẫu nhiên của dòng chảy và giá cả qua các quá trình dừng ARIMA/SARIMA sau khi biến đổi Box-Cox và sai phân mùa.

Điều kiện biên (Boundary conditions) của mô hình được xác lập chặt chẽ: dung tích hồ không vượt quá dung tích chết $S_{\min}$ và dung tích thiết kế $S_{\max}$; lưu lượng xả hạ du tuân thủ lưu lượng tối thiểu duy trì dòng chảy môi trường và an toàn đê điều chống lũ: $R_{\min} \le R_t \le R_{\max}$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ nghiêm ngặt thế giới quan thực chứng (Positivism) và phương pháp luận định lượng hệ thống. Mô hình được thiết kế đa cấp (Multi-level Design):

  • Cấp vĩ mô: Toàn bộ lưu vực sông Hồng – Thái Bình với các tương tác thủy văn dòng chính và biên giới quốc tế.
  • Cấp trung mô: Hai tiểu lưu vực điển hình đại diện cho hai cấu trúc khai thác khác nhau:
    • Hệ thống hồ Núi Cốc (Sông Công): Đại diện cho hệ thống hồ chứa vừa và nhỏ đa mục tiêu (tưới nông nghiệp, nước sinh hoạt công nghiệp cho tỉnh Thái Nguyên, phát điện nhỏ và du lịch hồ).
    • Hệ thống liên hồ chứa Sơn La – Hòa Bình (Sông Đà): Đại diện cho hệ thống hồ chứa chiến lược siêu lớn (chiếm 42% tổng lượng dòng chảy tại Sơn Tây, kiểm soát 6.246 MW điện năng và điều tiết lũ, cấp nước vụ Đông Xuân cho toàn bộ hạ du đồng bằng Bắc Bộ).
   Số liệu Thủy văn 21-25 năm              Số liệu Kinh tế & Xã hội
   (Dòng chảy đến, Mưa, Bốc hơi)           (Giá nông sản, Điện lực, CVM/TCM)

Quy trình nghiên cứu rigorous

Quy trình xử lý số liệu chuỗi thời gian ngẫu nhiên tuân thủ phương pháp luận Box-Jenkins tiêu chuẩn:

  1. Kiểm tra tính dừng và biến đổi dữ liệu: Chuỗi dòng chảy đến thực tế $Q_t$ của hồ Núi Cốc (1980–2000) và hồ Sơn La (1980–2004) có tính phi dừng và biến động phương sai theo mùa. Tác giả áp dụng phép biến đổi logarit tự nhiên $Y_t = \ln(Q_t)$ và kỹ thuật Box-Cox để ổn định phương sai, sau đó lấy sai phân theo chu kỳ tháng $(1 - B^{12})$ để triệt tiêu tính chu kỳ mùa vụ, tạo ra chuỗi dừng $Z_t$.
  2. Nhận dạng và ước lượng mô hình ARMA/ARIMA: Sử dụng đồ thị tự tương quan (ACF) và tự tương quan riêng phần (PACF) để xác định bậc $(p, q)$ tiềm năng. Thuật toán tự động tìm kiếm mô hình tối ưu trong không gian $p \le 26, q \le 26$ bằng cách tối thiểu hóa tiêu chuẩn thông tin Akaike hiệu chỉnh (AICC): $$\text{AICC} = -2\ln(L) + 2k + \frac{2k(k+1)}{n - k - 1}$$
  3. Kiểm định tính độc lập và phân phối phần dư: Chuỗi phần dư $e_t$ được kiểm định tính ồn trắng (white noise) và tính độc lập đồng nhất (IID) thông qua kiểm định Ljung-Box $Q$-statistic. Kết quả xác nhận phần dư hoàn toàn ngẫu nhiên, không còn cấu trúc tự tương quan.
  4. Thu thập số liệu định giá phi thị trường: Sử dụng phương pháp Đánh giá Ngẫu nhiên (Contingent Valuation Method - CVM) phỏng vấn người dân về mức sẵn lòng chi trả (WTP) cho nước sinh hoạt, và phương pháp Chi phí Du lịch (Travel Cost Method - TCM) của Hotelling-Clawson để dựng đường cầu giải trí du lịch tại hồ Núi Cốc.

Data và phân tích

Hệ số / Thông số kỹ thuật Hệ thống hồ Núi Cốc (Sông Công) Hệ thống liên hồ Sơn La – Hòa Bình (Sông Đà)
Khoảng thời gian chuỗi số liệu 21 năm (1980 – 2000) 25 năm (1980 – 2004)
Dung tích toàn bộ / hữu ích ($10^6\ m^3$) $175,5\ /\ 168,0$ Sơn La: $9.260\ /\ 6.504$; Hòa Bình: $9.450\ /\ 5.600$
Công suất phát điện lắp đặt Trạm nhỏ thủy nông ($< 5\ \text{MW}$) Sơn La ($2.400\ \text{MW}$) + Hòa Bình ($1.920\ \text{MW}$)
Nhu cầu cấp nước chính Tưới lúa Thái Nguyên, cấp nước gang thép, sinh hoạt, du lịch Tưới Đông Xuân đồng bằng Bắc Bộ, phát điện lưới quốc gia, cắt lũ
Mô hình chuỗi thời gian tối ưu $\text{ARIMA}(p, d, q)$ với phần dư ồn trắng $\text{ARMA}(p, q)$ đa biến sau biến đổi Box-Cox
Công cụ phần mềm giải tối ưu Lingo 16/17, Stata 13, R-Language Lingo 16/17 (phi tuyến toàn cầu - Global Solver)

Các kiểm định độ vững (robustness checks) được thực hiện thông qua mô phỏng Monte Carlo chạy lặp đa tình huống (Multi-scenario simulations) với các kịch bản:

  • Kịch bản dòng chảy ngẫu nhiên kết hợp cầu tưới tất định.
  • Kịch bản dòng chảy ngẫu nhiên kết hợp cầu tưới ngẫu nhiên do giá lúa biến động.
  • Kịch bản thay đổi dạng hàm chi phí hệ thống từ hàm mũ ($C = a \cdot e^{b \cdot Q}$) sang hàm đa thức bậc ba ($C = a_0 + a_1 Q + a_2 Q^2 + a_3 Q^3$).

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

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

  1. LƯỢNG HÓA THẶNG DƯ PHÚC LỢI RÒNG
     Hồ Núi Cốc đạt 103 - 107 tỷ VNĐ/năm; Liên hồ Sơn La - Hòa Bình 
     đạt hàng nghìn tỷ VNĐ vượt trội so với quy trình vận hành tĩnh.
  2. GIẢI QUYẾT XUNG ĐỘT PHÁT ĐIỆN - TƯỚI NƯỚC THÁNG 1 ĐẾN THÁNG 3
     Mô hình tối ưu hóa điều phối linh hoạt dung tích trữ, chuyển nước
     sang các đợt xả tăng cường vụ Đông Xuân mà không suy giảm công suất đỉnh.
  3. XÁC LẬP THUỘC TÍNH ĐỘC QUYỀN TỰ NHIÊN YẾU (WEAK NATURAL MONOPOLY)
     Hàm chi phí trung bình nhiều sản lượng (Ray Average Cost - RAC) 
     giảm liên tục, khẳng định tính cộng dưới chi phí (Subadditivity).
  4. NGUY CƠ KHI PHI ĐIỀU TIẾT HÓA HOÀN TOÀN
     Không có sự can thiệp của Chính phủ dẫn đến phân bổ lệch về phát điện,
     gây thâm hụt nghiêm trọng nước tưới nông nghiệp vùng hạ du.
  1. Lượng hóa chính xác thặng dư kinh tế tổng thể: Mô hình tối ưu hóa động xác định tổng lợi ích ròng xã hội của hệ thống hồ Núi Cốc đạt mức trung bình 103 đến 107 tỷ VNĐ/năm tùy thuộc vào kịch bản dòng chảy và dạng hàm chi phí (Bảng 3.2, 3.6, 3.9 trong luận án). Đối với hệ thống Sơn La – Hòa Bình, tổng giá trị kinh tế kết hợp giữa phát điện và xả nước tưới hạ du đạt hàng chục nghìn tỷ đồng mỗi năm, vượt trội đáng kể so với việc điều hành theo các biểu đồ điều phối cố định truyền thống.
  2. Giải quyết xung đột thương mại giữa Thủy điện (EVN) và Nông nghiệp (MARD): Trong mùa kiệt (tháng 1 đến tháng 5), đặc biệt là 3 đợt xả nước phục vụ gieo cấy lúa vụ Đông Xuân, mô hình tối ưu động đã chỉ ra quỹ đạo mực nước hồ tối ưu giúp giảm thiểu tổn thất phát điện mà vẫn đáp ứng 100% nhu cầu nước cho các trạm bơm hạ du. Việc tính toán giá bóng của nước tại các thời điểm căng thẳng chỉ ra rằng giá trị kinh tế của một đơn vị nước dùng cho tưới ở các giai đoạn nhạy cảm mùa vụ cao hơn giá trị phát điện biên.
  3. Phát hiện cấu trúc Độc quyền Tự nhiên yếu (Weak Natural Monopoly): Kết quả chạy mô hình tối ưu hóa động với các dạng hàm chi phí bậc ba và hàm mũ chứng minh rằng hệ thống liên hồ chứa Sơn La – Hòa Bình thỏa mãn điều kiện tính cộng dưới của chi phí (Subadditivity) trên toàn bộ miền sản lượng thực tế (Chương 4). Đường chi phí trung bình nhiều sản lượng (Ray Average Cost) có xu hướng giảm dần khi quy mô khai thác kết hợp tăng lên.
  4. Minh chứng về sự thất bại thị trường nếu phi điều tiết: Khi chạy mô hình mô phỏng tình huống thị trường tự do không có sự can thiệp của Chính phủ (Mục 4.3.4), các nhà máy thủy điện có xu hướng tích nước để phát điện vào giờ cao điểm giá điện cao, dẫn đến sự thiếu hụt nghiêm trọng nước tưới cho hạ du trong các tháng 2 và 3, làm giảm mạnh tổng phúc lợi xã hội của toàn lưu vực.

Implications đa chiều

  • Về mặt Lý thuyết: Luận án mở ra hướng nghiên cứu mới kết hợp giữa Lý thuyết Tối ưu hóa Hiện đại và Kinh tế học Tổ chức Công nghiệp trong ngành nước, cung cấp khung giải tích để xác định giá bóng tài nguyên nước theo không gian và thời gian.
  • Về mặt Phương pháp luận: Quy trình tích hợp mô hình chuỗi thời gian ARIMA/SARIMA xử lý qua Box-Cox với ngôn ngữ mô hình hóa tối ưu Lingo có thể chuyển giao áp dụng trực tiếp cho các lưu vực sông lớn khác ở Việt Nam như sông Đồng Nai, sông Vu Gia – Thu Bồn hay sông Ba.
  • Về mặt Thực tiễn và Chính sách:
    • Cung cấp cơ sở định lượng để Bộ Nông nghiệp & PTNT và Bộ Công Thương đàm phán Kế hoạch Xả nước Hàng năm (ARP) và Kế hoạch Điều tiết Tháng (MRP).
    • Đóng góp trực tiếp luận cứ cho việc xây dựng khung giá dịch vụ thủy lợi theo cơ chế thị trường có điều tiết theo tinh thần Luật Thủy lợi 2017, thay thế cơ chế thủy lợi phí trước đây.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Giới hạn về dữ liệu xuyên biên giới: Khu vực thượng nguồn sông Hồng thuộc Trung Quốc có 62 nhà máy thủy điện vừa và nhỏ với tổng dung tích trữ trên 4 tỷ $m^3$ đang vận hành hoặc quy hoạch, nhưng thông tin thủy văn chi tiết bị hạn chế. Điều này tạo ra độ bất định lớn đối với dòng chảy đến biên giới trong mùa lũ và cuối mùa khô.
  2. Đơn giản hóa một số đối tượng dùng nước: Do giới hạn dung lượng tính toán, mô hình tập trung sâu vào hai ngành sử dụng nước lớn nhất là tưới lúa và phát điện; các nhu cầu giao thông thủy, thủy sản và chất lượng môi trường nước mới chỉ được đưa vào dưới dạng các ràng buộc bất đẳng thức biên thay vì xây dựng các hàm cầu nội sinh hoàn chỉnh.
  3. Giả định tính dừng cục bộ: Việc sử dụng các mô hình ARIMA/SARIMA dựa trên giả định chuỗi dừng sau biến đổi có thể gặp thách thức trong bối cảnh biến đổi khí hậu cực đoan dài hạn và hiện tượng hạ thấp lòng dẫn sông Hồng diễn ra dữ dội trong những năm gần đây.

Chương trình nghiên cứu tương lai (Future Research Agenda) cần tập trung vào 4 hướng:

  • Mở rộng mô hình tối ưu hóa động thích ứng (Adaptive Dynamic Programming) có kết nối dữ liệu viễn thám và dự báo khí tượng thủy văn thời gian thực.
  • Xây dựng mô hình cân bằng thủy lực - chất lượng nước tích hợp để tối ưu hóa đồng thời cả lượng nước và nồng độ ô nhiễm xả thải.
  • Nghiên cứu cơ chế hạn ngạch và thị trường trao đổi quyền sử dụng nước (Water Rights Trading) tại các tiểu lưu vực có mức độ cạnh tranh cao.
  • Phân tích tương tác chiến lược đa quốc gia sử dụng Lý thuyết Trò chơi (Game Theory) đối với việc điều tiết các hồ chứa thượng nguồn xuyên biên giới.

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

Nghiên cứu mang lại những ảnh hưởng sâu rộng trên nhiều bình diện:

  • Tác động Học thuật: Tạo nền tảng phương pháp luận chuẩn mực cho các nghiên cứu sinh và học viên cao học chuyên ngành Quản lý Tài nguyên Nước, Kinh tế Tài nguyên và Thủy văn Thủy lực tại Việt Nam; mở ra tiềm năng trích dẫn cao trong các tạp chí khoa học quốc tế uy tín thuộc danh mục ISI/Scopus về mô hình hóa kinh tế thủy văn (Hydro-economic Systems).
  • Chuyển đổi Ngành và Doanh nghiệp: Giúp Tập đoàn Điện lực Việt Nam (EVN), Trung tâm Điều độ Hệ thống Điện Quốc gia (A0) và các Công ty Khai thác Công trình Thủy lợi tối ưu hóa kế hoạch huy động tổ máy và lịch lấy nước, tiết kiệm hàng triệu mét khối nước xả từ các hồ chứa lớn.
  • Tác động Chính sách: Kết quả nghiên cứu được sử dụng làm tài liệu tham khảo chất lượng cao trong các hội thảo quốc tế lớn như Hội thảo "Định giá tài nguyên nước" của Ngân hàng Thế giới (World Bank) và phản biện Đề tài "Kế toán Nước" của Viện Chiến lược, Chính sách Tài nguyên và Môi trường (Bộ TN&MT).

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

  • Nghiên cứu sinh Tiến sĩ và Nhà khoa học: Tiếp cận một cấu trúc mô hình hóa hoàn chỉnh kết hợp kinh tế lượng chuỗi thời gian, lý thuyết tối ưu hóa toán học phi tuyến và kinh tế học vi mô ứng dụng.
  • Nhà quản lý và Hoạch định Chính sách (MARD, MOIT, MONRE, Văn phòng Chính phủ): Có công cụ khoa học chính xác để kiểm chứng tác động của các quyết định điều hành liên hồ chứa, hạn chế các quyết định mang tính hành chính cảm tính.
  • Tập đoàn Năng lượng và Doanh nghiệp Khai thác Thủy lợi: Tối đa hóa doanh thu phát điện và giảm chi phí vận hành bơm tưới thông qua việc dự báo chính xác các dòng chảy tối ưu.
  • Cộng đồng Xã hội: Hàng triệu hộ nông dân vùng đồng bằng sông Hồng được đảm bảo an ninh nguồn nước tưới trong vụ Đông Xuân, đồng thời toàn xã hội được hưởng lợi từ việc cấp điện ổn định và kiểm soát rủi ro ngập lụt mùa mưa bão.

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 lý thuyết độc đáo nhất là việc tích hợp thành công Lý thuyết Độc quyền Tự nhiên (Baumol et al., 1982) vào bài toán tối ưu hóa động phân bổ nguồn nước lưu vực sông. Luận án đã mở rộng mô hình phân bổ Pareto truyền thống của James Douglas (1971) bằng cách chứng minh tính chất cộng dưới của hàm chi phí (Cost Subadditivity) trong hệ thống hồ chứa đa mục tiêu, từ đó cung cấp cơ sở lý luận kinh tế học vững chắc khẳng định tính tất yếu của sự điều tiết công quyền của Nhà nước đối với cơ sở hạ tầng nước chiến lược.

2. Điểm đổi mới phương pháp luận so với ít nhất hai nghiên cứu tiền nhiệm?

So với mô hình Aquarius của Nhóm nghiên cứu Rừng Rocky Mountain (USDA) vốn chỉ dừng lại ở bài toán tĩnh và thiếu cấu trúc hàm chi phí, luận án đã phát triển mô hình tối ưu hóa động đa thời kỳ có tích hợp đầy đủ các hàm chi phí hệ thống (bậc ba và hàm mũ). So với nghiên cứu tối ưu hóa lưu vực sông Đồng Nai của Claudia Ringler và Nguyễn Vũ Huy (2004) sử dụng mô hình tĩnh trên GAMS, luận án đã thực hiện bước đột phá khi nội sinh hóa các quá trình ngẫu nhiên của dòng chảy đến và biến động giá nông sản thông qua chuỗi thời gian ARIMA/SARIMA được xử lý dừng hóa bằng biến đổi Box-Cox.

3. Phát hiện bất ngờ nhất (counter-intuitive) có số liệu chứng minh là gì?

Phát hiện bất ngờ nhất là trong điều kiện không có sự điều tiết của Chính phủ (mô hình cạnh tranh thị trường thuần túy), việc gia tăng dung tích trữ của các hồ chứa thượng nguồn không đồng nghĩa với việc nâng cao an ninh nguồn nước hạ du. Ngược lại, do tín hiệu giá điện giờ cao điểm chi phối hành vi tối đa hóa lợi nhuận của doanh nghiệp phát điện, lượng nước xả trong các tháng kiệt phục vụ gieo cấy lúa (tháng 1–3) bị sụt giảm nghiêm trọng, dẫn đến tổng thặng dư kinh tế toàn xã hội giảm sút so với phương án vận hành có ràng buộc điều tiết công cộng.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (replication protocol) không?

Có. Luận án mô tả chi tiết toàn bộ cấu trúc mã lệnh toán học trên môi trường phần mềm Lingo 16/17, trình bày đầy đủ bảng thông số các hàm cầu, hàm sản xuất nông nghiệp Cobb-Douglas, các phương trình hồi quy ARIMA với các hệ số tự hồi quy và trung bình trượt cụ thể, cùng hệ thống phụ lục chi tiết lưu trữ dữ liệu dòng chảy 21 năm (Núi Cốc) và 25 năm (Sơn La – Hòa Bình). Các nhà nghiên cứu khác hoàn toàn có thể sử dụng lại giao thức này để tái lập kết quả.

5. Khung chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Luận án định hình lộ trình nghiên cứu 10 năm tập trung vào 3 trụ cột: (i) Phát triển mô hình kinh tế thủy văn thời gian thực (Real-time Hydro-economic Forecasting) ứng dụng trí tuệ nhân tạo; (ii) Lượng hóa toàn diện giá trị kinh tế của dòng chảy môi trường và dịch vụ hệ sinh thái đất ngập nước; (iii) Thiết kế cơ chế phân bổ hạn ngạch nước và thị trường trao đổi quyền sử dụng nước có sự tham gia của khu vực tư nhân dưới sự giám sát của Nhà nước.

Kết luận

  1. Luận án đã thiết lập thành công khung lý thuyết và phương pháp luận hoàn chỉnh cho bài toán phân bổ hiệu quả kinh tế tài nguyên nước trên lưu vực sông Hồng – Thái Bình, kết hợp nhuần nhuyễn giữa kinh tế học phúc lợi và kỹ thuật tối ưu hóa động.
  2. Xây dựng và giải thành công mô hình tối ưu hóa động ngẫu nhiên trên phần mềm chuyên dụng Lingo 16/17 cho hai hệ thống tiêu biểu: hồ Núi Cốc (lưu vực sông Công) và liên hồ chứa Sơn La – Hòa Bình (lưu vực sông Đà).
  3. Đề xuất quy trình chuẩn hóa ứng dụng mô hình kinh tế lượng chuỗi thời gian ARIMA/SARIMA kết hợp biến đổi Box-Cox để mô phỏng chính xác tính ngẫu nhiên của dòng chảy thủy văn và biến động giá nông sản.
  4. Lượng hóa chi tiết hàm cầu và hàm giá trị kinh tế của các ngành dùng nước cạnh tranh chính (nông nghiệp, phát điện, cấp nước sinh hoạt, du lịch), xác lập giá bóng của nước tại các thời điểm căng thẳng mùa vụ.
  5. Chứng minh khoa học thuộc tính độc quyền tự nhiên yếu của hệ thống công trình thủy công liên hồ chứa, làm sáng tỏ cơ sở lý luận kinh tế cho việc Nhà nước duy trì điều tiết và điều hành thống nhất thay vì thả nổi thị trường.
  6. Cung cấp luận cứ khoa học thực tiễn vững chắc phục vụ đắc lực cho công tác thực thi Luật Thủy lợi 2017, hoàn thiện Quy trình vận hành liên hồ chứa theo Quyết định số 1287/QĐ-TTg, đặt nền móng vững chắc cho sự phát triển bền vững của ngành nước Việt Nam trong kỷ nguyên mới.