Tổng quan về luận án
Bối cảnh biến đổi khí hậu toàn cầu, sự bùng nổ dân số và tốc độ tăng trưởng kinh tế nhanh chóng tại lưu vực sông Hồng – Thái Bình (LVSHTB) đang tạo ra áp lực chưa từng có lên an ninh nguồn nước. Bao phủ diện tích 169.000 km² trải dài qua ba quốc gia Việt Nam, Trung Quốc và Lào, LVSHTB nuôi dưỡng 29 triệu người (chiếm 95% dân số miền Bắc) và đóng góp xấp xỉ 20% GDP quốc gia. Dù tổng lượng dòng chảy trung bình đạt tới 135 tỷ m³ (với 82,54 tỷ m³, tương đương 61,1% hình thành trong lãnh thổ Việt Nam và 52,46 tỷ m³, chiếm 38,9% bắt nguồn từ Trung Quốc), sự phân bố dòng chảy không đều giữa các mùa và sự phụ thuộc nghiêm trọng vào 62 công trình thủy điện thượng nguồn phía Trung Quốc (tổng dung tích trên 4 tỷ m³) đặt ra thách thức sinh tử cho việc điều phối nguồn nước hạ du.
Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại nhiều thập kỷ qua là sự thống trị của cách tiếp cận kỹ thuật thuần túy (mô phỏng thủy văn - thủy lực) và sự thiếu vắng các mô hình phân bổ kinh tế liên thời gian có khả năng phản ánh giá trị kinh tế thực tế của nước dưới dạng hàm cầu động. Như McKinney và cộng sự (1999) đã tổng kết: "Bất kể tầm quan trọng to lớn của các biến kinh tế trong phân bổ và quản lý tài nguyên nước, các nghiên cứu tài nguyên nước nói chung vẫn bị thống trị bởi các nghiên cứu thủy văn, thủy lực, tức là quy hoạch tài nguyên nước theo quan điểm kỹ thuật". Các quy trình vận hành liên hồ chứa truyền thống tại Việt Nam (như các nghiên cứu của Hà Văn Khối và Lê Bảo Trung, 2003; Viện Quy hoạch Thủy lợi, 2005, 2007) chủ yếu giải quyết bài toán thủy lực tĩnh hoặc tách rời mùa kiệt và mùa lũ, chưa tích hợp hành vi kinh tế của người sử dụng nước trong điều kiện bất định.
Nghiên cứu của tác giả Đào Văn Khiêm (Chuyên ngành Quy hoạch và quản lý tài nguyên nước, Mã số: 62-62-30-01, Trường Đại học Thủy lợi, dưới sự hướng dẫn của GS.TS Nguyễn Quang Kim) giải quyết triệt để khoảng trống này thông qua việc thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:
- Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng và lượng hóa hàm mục tiêu kinh tế tích hợp đầy đủ hàm cầu động, hàm giá trị sử dụng nước (tưới nông nghiệp, thủy điện, sinh hoạt, du lịch) và cấu trúc hàm chi phí trong điều kiện thị trường có điều tiết tại Việt Nam?
- Câu hỏi nghiên cứu 2 (RQ2): Mô hình tối ưu hóa động ngẫu nhiên phi tuyến có thể giải quyết xung đột lợi ích liên thời gian giữa các ngành sử dụng nước tại các tiểu lưu vực phức tạp như thế nào dưới tác động của biến ngẫu nhiên dòng chảy và giá cả thị trường?
- Câu hỏi nghiên cứu 3 (RQ3): Hệ thống hồ chứa đa mục tiêu tại LVSHTB có tồn tại thuộc tính "độc quyền tự nhiên" (natural monopoly) hay không, và mức độ mạnh/yếu của cấu trúc này ảnh hưởng ra sao đến cơ chế can thiệp của Nhà nước?
Hệ thống giả thuyết nghiên cứu bao gồm:
- Giả thuyết 1 (H1): Việc chuyển đổi từ phân bổ cơ học sang mô hình tối đa hóa tổng lợi ích ròng xã hội (Net Economic Surplus) dựa trên thặng dư kinh tế biên sẽ gia tăng đáng kể hiệu quả kinh tế tổng thể mà vẫn đảm bảo các ràng buộc an toàn phòng lũ và mực nước giới hạn.
- Giả thuyết 2 (H2): Dòng chảy đến tự nhiên và cầu nước tưới nông nghiệp tuân theo các phân phối ngẫu nhiên (ARIMA/ngẫu nhiên thực nghiệm), và quỹ đạo tối ưu tìm được qua phương pháp quy hoạch phi tuyến ngẫu nhiên có độ ổn định cao trước các biến động cực đoan.
- Giả thuyết 3 (H3): Đặc tính chi phí biên giảm dần và chi phí vốn hạ tầng cố định khổng lồ tạo nên cấu trúc độc quyền tự nhiên yếu (weak natural monopoly) trong vận hành phân phối nước của hệ thống hồ chứa, đòi hỏi khung pháp lý điều tiết giá và phân bổ sản lượng tối ưu thay vì thả nổi thị trường tự do.
Khung lý thuyết của luận án được định hình từ Tuyên bố Dublin 1992 ("Nước có giá trị kinh tế và cần được coi như một hàng hóa kinh tế... ngành nước là một ngành nổi tiếng vì có nhiều 'thất bại thị trường'"), Lý thuyết Tối ưu hóa Kinh tế Phúc lợi Pareto, Lý thuyết Quy hoạch Phi tuyến và Ràng buộc Hỗn hợp Kuhn-Tucker (Sundaram, 1996), cùng Lý thuyết Độc quyền Tự nhiên của Baumol, Panzar và Willig (1982). Đóng góp đột phá của luận án được kiểm chứng thực nghiệm trên chuỗi dữ liệu dài hạn: 21 năm (1980–2000) đối với hệ thống hồ Núi Cốc (sông Công) và 25 năm (1980–2005) đối với liên hệ thống Sơn La – Hòa Bình (sông Đà) với tổng công suất phát điện lên tới 4.320 MW (Sơn La: 2.400 MW, Hòa Bình: 1.920 MW) và tổng dung tích phòng lũ kết hợp 7,0 tỷ m³.
Literature Review và Positioning
Tổng quan y văn quốc tế cho thấy sự tiến hóa rõ nét của ba dòng nghiên cứu chính:
- Dòng nghiên cứu mô hình hóa kỹ thuật - mô phỏng thủy văn: Chiếm ưu thế tuyệt đối trong thế kỷ 20 với các mô hình tiêu biểu như MIKE 11, HEC-ResSim, MIKE BASIN. Dòng nghiên cứu này mô phỏng chi tiết dòng chảy thủy lực, chuyển tiếp sóng lũ và cân bằng nước vật lý nhưng xem cầu nước là đại lượng cố định ngoại sinh (Loucks và Van Beek, 2005).
- Dòng nghiên cứu kinh tế học phân bổ tài nguyên nước: Khởi xướng bởi Douglas James và Lee (1971), Bowen và Young (1985), Dinar và Letey (1996), và Griffin (2006). Dòng nghiên cứu này nhấn mạnh việc thiết lập hàm sản xuất mùa vụ, đường cầu nước nông nghiệp phi tuyến dựa trên phương pháp bốc thoát hơi nước và mô hình quy hoạch toán học nông hộ để suy ra ý muốn thanh toán cận biên (Marginal Willingness to Pay - MWTP).
- Dòng nghiên cứu mô hình kinh tế - thủy văn tổng hợp (Integrated Hydro-Economic Modeling): Được định hình bởi Rosegrant và Meinzen-Dick (1996), McKinney, Ximing Cai và cộng sự (1999), Booker và Colby (1995) tại lưu vực sông Colorado, và hệ thống mô hình hóa Aquarius do nhóm nghiên cứu Rocky Mountain Station (Cơ quan Kiểm lâm Hoa Kỳ - USDA Forest Service) và Đại học Bang Colorado phát triển (Diaz, Brown, Sveinsson, 1997).
Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:
- Tranh luận 1 (Mô hình hóa cục bộ vs. Mô hình hóa hệ thống): Các tác giả như Howe và Orr (1974), Noel và Howitt (1982), Lee và Howitt (1996) ủng hộ kết nối lỏng (loose-coupling) giữa mô hình thủy văn độc lập và mô hình kinh tế đầu ra để bảo toàn độ chi tiết thủy lực. Ngược lại, Booker (1996), Henderson và Lord (1995) chứng minh rằng chỉ có tiếp cận mô hình hóa hệ thống đồng thời (unified system approach) mới giải quyết được tính liên hồi nội sinh giữa động thái trữ nước hồ chứa và biến động giá trị kinh tế biên.
- Tranh luận 2 (Cơ chế phân bổ hành chính vs. Cơ chế thị trường điều tiết): Trường phái điều phối hành chính truyền thống cho rằng nước là hàng hóa công cộng thuần túy cần cấp phát theo chỉ tiêu định mức. Ngược lại, trường phái kinh tế tài nguyên (Young, 2005; Griffin, 2006) khẳng định việc thiếu vắng tín hiệu giá và hàm lợi ích biên dẫn đến tổn thất vô ích (deadweight loss) nghiêm trọng trong phân bổ liên ngành.
[TIẾP CẬN KỸ THUẬT THUẦN TÚY] [TIẾP CẬN KINH TẾ ĐƠN LẺ]
- Mô phỏng MIKE 11 / MIKE BASIN - Đánh giá CVM/TCM tĩnh
- Ràng buộc thủy lực vật lý - Thiếu động học thủy văn
\ /
\ /
[KHOẢNG TRỐNG: THIẾU TỐI ƯU HÓA ĐỘNG NGẪU NHIÊN LIÊN HỒ VÀ
LƯỢNG HÓA CẤU TRÚC ĐỘC QUYỀN TỰ NHIÊN DƯỚI RÀNG BUỘC HỖN HỢP]
[VỊ THẾ ĐỘT PHÁ CỦA LUẬN ÁN ĐÀO VĂN KHIÊM (2017)]
- Tích hợp mô hình kinh tế - thủy văn Aquarius mở rộng
- Thuật toán phi tuyến ngẫu nhiên (Lingo 16/17 + ARIMA)
- Định vị cấu trúc Độc quyền tự nhiên yếu cấp lưu vực
Tại Việt Nam, các nghiên cứu trước đây như Đề tài cấp Nhà nước của GS.TS Lê Kim Truyền (2005–2006) hay nghiên cứu của Viện Quy hoạch Thủy lợi (2005, 2007) chỉ dừng lại ở quy hoạch cấp nước mùa kiệt hoặc ứng dụng phần mềm GAMS ở dạng tĩnh, chưa lượng hóa được hàm cầu dựa trên ý muốn chi trả của người sử dụng nước. Nghiên cứu của Claudia Ringler và Nguyễn Vũ Huy (2004) tại lưu vực sông Đồng Nai dù áp dụng tiếp cận phân bổ tối ưu nhưng dựa trên khung tối ưu tĩnh, bỏ qua các cấu phần động học và bất định ngẫu nhiên của dòng chảy.
Luận án của Đào Văn Khiêm định vị tiên phong tại giao điểm giữa kỹ thuật tài nguyên nước và kinh tế học công cộng hiện đại. So với hai nghiên cứu quốc tế điển hình là mô hình lưu vực sông Colorado của Booker và Colby (1995) và mô hình phân bổ nước ngầm - nước mặt Nam California của Dinar và Letey (1996), luận án không chỉ kế thừa khung phân tích hàm lợi ích tổng hợp từ hệ thống Aquarius mà còn phát triển đột phá thêm hai cấu phần chưa từng có: (1) Tích hợp hàm cấu trúc chi phí (hàm mũ và hàm bậc ba) vào hàm mục tiêu tối đa hóa lợi ích ròng, và (2) Thiết lập quy trình toán học kiểm định tính cộng dưới chi phí (subadditivity) để chứng minh đặc tính độc quyền tự nhiên của hệ thống liên hồ chứa.
Đó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 lý thuyết kinh tế tài nguyên và tổ chức ngành thông qua các đóng góp cụ thể:
- Mở rộng Lý thuyết Tối ưu hóa Kinh tế Tài nguyên Nước (Welfare Economics in Water Resources): Kế thừa tư tưởng của Douglas James, Lee (1971) và Griffin (2006), nghiên cứu đã phát triển khung xác định giá bóng (shadow price) liên thời gian của nước thông qua nhân tử Lagrange $\lambda(t)$ trong bài toán tối ưu hóa phi tuyến có ràng buộc hỗn hợp Kuhn-Tucker. Giá bóng phản ánh chính xác chi phí cơ hội biên của một đơn vị nước trữ trong hồ đối với các mục tiêu cạnh tranh (thủy điện vs. tưới tiêu hạ du).
- Tiên phong áp dụng Lý thuyết Độc quyền Tự nhiên (Theory of Natural Monopoly) vào Lưu vực sông: Dựa trên nền tảng lý thuyết cấu trúc thị trường của Baumol, Panzar và Willig (1982), luận án đã mở rộng khái niệm tính cộng dưới chi phí sang ngành nước. Nghiên cứu chứng minh rằng một hệ thống liên hồ chứa duy nhất vận hành tập trung sẽ tối thiểu hóa tổng chi phí xã hội so với việc phân chia quyền vận hành cho nhiều đơn vị độc lập.
- Phát triển Mô hình Cầu động Nước Nông nghiệp: Bổ cứu cho lý thuyết hàm sản xuất mùa vụ tĩnh của Hanks và Hill (1980), nghiên cứu thiết lập phương pháp ước lượng hàm cầu động đối với nước tưới, phản ánh sự biến thiên của giá trị sản phẩm nông nghiệp (giá lúa ngẫu nhiên) và các giai đoạn sinh trưởng nhạy cảm của cây trồng.
Hệ thống mệnh đề lý thuyết (Theoretical Propositions) được xác lập:
- Mệnh đề 1 (P1): Tổng lợi ích kinh tế ròng của lưu vực đạt cực đại khi và chỉ khi giá trị sản phẩm cận biên của nước giữa tất cả các điểm sử dụng ngoài dòng chảy (tưới, sinh hoạt, công nghiệp) và trên dòng chảy (thủy năng, sinh thái) cân bằng với giá bóng của tài nguyên nước tại từng thời đoạn $t$.
- Mệnh đề 2 (P2): Do tồn tại hiệu quả kinh tế theo quy mô và chi phí hạ tầng cố định cực lớn, hàm chi phí vận hành hồ chứa thỏa mãn điều kiện $C(\sum q_i) < \sum C(q_i)$, xác lập cấu trúc độc quyền tự nhiên yếu, đòi hỏi cơ chế định giá biên bậc hai (two-part tariff) có sự điều tiết của Nhà nước.
- Mệnh đề 3 (P3): Quỹ đạo tích - xả tối ưu của hồ chứa trong điều kiện ngẫu nhiên thủy văn không hội tụ về một điểm tĩnh mà tạo thành một dải phân phối xác suất ổn định, trong đó mức độ chiết khấu liên thời gian quyết định tỷ lệ đánh đổi giữa xả nước phát điện mùa khô và duy trì an toàn phòng lũ cuối mùa kiệt.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp liền mạch ba trụ cột lý thuyết: (1) Lý thuyết Mô hình hóa Hệ thống Thủy văn LVSHTB, (2) Lý thuyết Tối ưu hóa Phi tuyến Ngẫu nhiên Đa chiều, và (3) Lý thuyết Kinh tế học Điều tiết Độc quyền Tự nhiên.
Mô hình toán học cốt lõi mở rộng từ hệ thống Aquarius được thiết lập dưới dạng bài toán tối đa hóa Tổng Thặng dư Kinh tế Ròng (Net Economic Surplus - $TB_{Net}$):
$$\max TB_{Net} = \sum_{t=1}^{np} \sum_{u=1}^{nu} \int_{0}^{q_{ut}} f_{ut}(x),dx - \sum_{t=1}^{np} C_t\left(S_t, O_t, q_{ut}\right)$$
Phụ thuộc vào hệ thống ràng buộc hỗn hợp phi tuyến nghiêm ngặt:
- Ràng buộc cân bằng nước hồ chứa liên thời gian:
$$S(t+1) = S(t) + I(t) - O(t) - E(S(t)) - L(S(t))$$
Trong đó $S(t)$ là dung tích trữ, $I(t)$ là dòng chảy đến (ngẫu nhiên), $O(t)$ là tổng lượng xả quyết định, $E$ là bốc hơi mặt hồ, $L$ là tổn thất thấm/tràn.
- Ràng buộc đặc trưng kỹ thuật hồ chứa: Quan hệ phi tuyến giữa cao trình mực nước $H(t)$, dung tích $S(t)$, và diện tích mặt nước $A(t)$:
$$H(t) = \phi_1(S(t)), \quad A(t) = \phi_2(S(t))$$
- Ràng buộc công suất phát điện và dung tích phòng lũ:
$$N_{min} \le N_t = 9.81 \cdot \eta \cdot O_{power}(t) \cdot H_{net}(t) \le N_{max}$$
$$S(t) \le S_{max} - V_{flood}(t)$$
- Ràng buộc phi âm và biên giới hạn tài nguyên: $q_{ut} \ge 0, \quad O(t) \ge O_{min_env}$.
Điều kiện biên (Boundary conditions) của khung phân tích được xác định rõ: Khung mô hình áp dụng cho hệ thống hồ chứa có điều tiết năm và nhiều năm; dữ liệu kinh tế giả định thị trường sản phẩm nông nghiệp và điện năng vận hành theo cơ chế thị trường định hướng xã hội chủ nghĩa có sự điều tiết của Chính phủ; chuỗi thủy văn ngẫu nhiên tuân theo quy luật dừng có điều kiện sau khi khử xu thế mùa vụ.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án áp dụng thế giới quan Hiện thực Phê phán (Critical Realism) kết hợp Thực chứng Thực dụng (Pragmatic Positivism). Nghiên cứu thừa nhận các quy luật vật lý thủy văn tồn tại khách quan nhưng hành vi kinh tế của các tác nhân sử dụng nước bị chi phối bởi thể chế, giá cả và chính sách điều tiết. Thiết kế nghiên cứu phối hợp phương pháp hỗn hợp đa tầng (Multi-level Mixed Design):
[THIẾT KẾ PHƯƠNG PHÁP NGHIÊN CỨU ĐA TẦNG]
[CẤU PHẦN THỦY VĂN & VẬT LÝ] [CẤU PHẦN KINH TẾ & HÀNH VI]
- Thu thập chuỗi số liệu NCHMS (1980-2005) - Điều tra phỏng vấn CVM, TCM
- Dòng chảy đến tại Sơn Tây (118 tỷ m³/năm) - Ước lượng hàm cầu nước tưới, sinh hoạt
- Mô hình hóa ngẫu nhiên chuỗi thời gian ARIMA - Xây dựng cấu trúc chi phí (mũ & bậc ba)
[TỔNG HỢP VÀO MÔ HÌNH HỆ THỐNG AQUARIUS MỞ RỘNG]
- Xây dựng hàm mục tiêu Net Economic Surplus
- Tích hợp hệ thống ràng buộc hỗn hợp Kuhn-Tucker
[GIẢI BÀI TOÁN TỐI ƯU TRÊN MÔI TRƯỜNG LINGO 16/17]
- Thuật toán Non-Linear Programming (NLP/SLQP)
- Kiểm định Monte Carlo tính ổn định nghiệm (Robustness)
[PHÂN TÍCH ĐẶC TÍNH CẤU TRÚC ĐỘC QUYỀN TỰ NHIÊN]
- Đánh giá tính cộng dưới chi phí (Subadditivity)
- Đề xuất chính sách điều tiết cấp lưu vực
Quy mô không - thời gian và tập mẫu nghiên cứu được chuẩn hóa chi tiết:
- Tiểu lưu vực sông Công (Hệ thống Núi Cốc): Chuỗi thời gian phân tích 21 năm (1980–2000), phân giải theo bước thời gian tháng (252 tháng tính toán). Cung cấp nước tưới cho hàng chục ngàn hecta nông nghiệp Thái Nguyên, cấp nước sinh hoạt, du lịch hồ Núi Cốc và phát điện nhỏ.
- Tiểu lưu vực sông Đà (Liên hồ Sơn La – Hòa Bình): Chuỗi thời gian phân tích 25 năm (1980–2005), 300 bước thời gian tháng. Công suất lắp đặt liên hồ: 4.320 MW; điều tiết nguồn nước cho toàn bộ đồng bằng Bắc Bộ trong các đợt xả tăng cường vụ Đông Xuân.
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu được kiểm soát chất lượng qua các giao thức nghiêm ngặt:
- Thu thập và xử lý số liệu đầu vào: Dữ liệu thủy văn (lưu lượng đến, mực nước, lượng mưa) do Trung tâm Khí tượng Thủy văn Quốc gia (NCHMS) cung cấp. Số liệu nhu cầu nước mùa vụ do Cục Trồng trọt (DCP/MARD) và Ban Quản lý Khai thác Thủy lợi (DWR/MARD) cung cấp. Thông số kỹ thuật phát điện và giá điện do EVN và Cục Điều tiết Điện lực (ERAV) cung cấp.
- Kỹ thuật điều tra phi thị trường: Sử dụng phương pháp Đánh giá Ngẫu nhiên (Contingent Valuation Method - CVM) với kịch bản phỏng vấn cấu trúc để đo lường ý muốn chi trả (WTP) của các hộ dùng nước sinh hoạt và tưới tiêu; áp dụng Phương pháp Chi phí Du lịch (Travel Cost Method - TCM) do Harold Hotelling (1947) và Clawson - Knetsch (1966) phát triển để xây dựng hàm cầu giải trí/du lịch tại hồ Núi Cốc.
- Kỹ thuật Tam giác đạc (Triangulation): Kết hợp tam giác đạc dữ liệu (chuỗi quan trắc lịch sử + mô phỏng stochastic), tam giác đạc phương pháp (kinh tế lượng + quy hoạch toán học) và tam giác đạc lý thuyết (kỹ thuật thủy lợi + kinh tế học phúc lợi).
- Độ tin cậy và giá trị đo lường: Các phương trình hồi quy hàm cầu kinh tế lượng được kiểm định đa cộng tuyến, phương sai sai số thay đổi và tự tương quan, đạt hệ số tin cậy $R^2 > 0,85$ và hệ số Cronbach's $\alpha > 0,80$ đối với bảng hỏi điều tra CVM.
Data và phân tích
Đặc trưng dữ liệu thủy văn đầu vào thể hiện tính phức tạp cao độ của LVSHTB: Lưu lượng năm trung bình tại trạm thủy văn Sơn Tây đạt 118 tỷ m³ (tương đương lưu lượng bình quân $3.740\text{ m}^3/\text{s}$). Trong đó, sông Đà đóng góp tỷ trọng lớn nhất với 42% tổng lượng dòng chảy; sông Lô chiếm 25,4%; sông Thao có diện tích lưu vực tương đương sông Đà nhưng lượng dòng chảy chỉ chiếm 19%.
Công cụ và kỹ thuật số tiên tiến được huy động:
- Xử lý chuỗi thời gian ngẫu nhiên: Sử dụng mô hình Tự hồi quy Trung bình Trượt tích hợp (ARIMA) trên phần mềm Stata 13 và R Language để phân tích và sinh các kịch bản chuỗi dòng chảy đến ngẫu nhiên $I(t)$ và chuỗi biến động giá nông sản.
- Phần mềm tối ưu hóa chuyên dụng Lingo 16/17: Sử dụng solver phi tuyến (Nonlinear Solver) và kỹ thuật Quy hoạch Tuyến tính Toàn phương Ngẫu nhiên (Stochastic Linear Quadratic Programming - SLQP) có bản quyền của hãng LINDO Systems. Lingo vượt trội hơn phần mềm truyền thống GAMS nhờ khả năng xử lý mượt mà các ma trận ràng buộc phi tuyến quy mô lớn với hàng chục ngàn biến quyết định.
- Kiểm tra tính ổn định nghiệm (Robustness Checks): Thực hiện mô phỏng Monte Carlo với hàng trăm kịch bản nhiễu ngẫu nhiên về dòng chảy đến và nhu cầu dùng nước. Thử nghiệm so sánh độ nhạy của nghiệm tối ưu dưới hai cấu trúc chi phí khác nhau: mô hình chi phí hàm mũ $C(Q) = a \cdot e^{bQ}$ và mô hình chi phí bậc ba $C(Q) = c_0 + c_1 Q + c_2 Q^2 + c_3 Q^3$. Kết quả xác nhận các quỹ đạo trữ - xả tối ưu duy trì tính hội tụ cao với khoảng tin cậy 95% hẹp.
Phát hiện đột phá và implications
Những phát hiện then chốt
Phân tích định lượng từ các mô hình tối ưu hóa động tất định và ngẫu nhiên mang lại các phát hiện khoa học mang tính bước ngoặt:
- Gia tăng vượt trội của Tổng Lợi ích Ròng Kinh tế: Việc áp dụng quy tắc phân bổ tối ưu theo giá trị biên làm tăng tổng lợi ích ròng toàn hệ thống từ 18,5% đến 24,2% so với quy trình vận hành quy chuẩn tĩnh hiện hành. Mô hình chứng minh rằng việc cố định cơ học lượng xả tưới vụ Đông Xuân gây tổn thất phúc lợi kinh tế lớn cho ngành phát điện trong các tháng cao điểm mùa khô (tháng 1 đến tháng 3), trong khi cùng một lượng nước đó nếu được điều tiết linh hoạt theo hàm cầu động sẽ mang lại thặng dư kinh tế kép cho cả hai khu vực.
- Nhận diện thực nghiệm 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 chi phí cho hệ thống hồ Núi Cốc và liên hồ Sơn La – Hòa Bình xác nhận đường chi phí trung bình dài hạn (Long-run Average Cost - LAC) giảm dần trên toàn bộ miền sản lượng xả hiệu dụng. Hàm chi phí thỏa mãn tính cộng dưới (subadditivity):
$$C(Q_{tưới} + Q_{điện}) < C(Q_{tưới}) + C(Q_{điện})$$
Tuy nhiên, tại các tháng kiệt đỉnh điểm (tháng 3, tháng 4), tính cộng dưới bị suy giảm cục bộ do sự gia tăng đột biến của chi phí cơ hội biên, chứng minh hệ thống mang cấu trúc độc quyền tự nhiên yếu.
- Đặc tính Bền vững của Quỹ đạo Tối ưu Ngẫu nhiên: Khi đưa chuỗi dòng chảy đến ngẫu nhiên ARIMA và biến động giá lúa vào mô hình ngẫu nhiên, các quỹ đạo tối ưu (optimal trajectories) của mực nước hồ chứa không bị phá vỡ mà dao động xung quanh một dải quỹ đạo tiệm cận ổn định. Phân phối xác suất của tổng lợi ích ròng (Probability Density Function - PDF) có dạng đơn đỉnh lệch phải, khẳng định tính khả thi trong thực tiễn điều hành thời gian thực.
- Hiệu ứng Đánh đổi Liên thời gian (Inter-temporal Trade-off) và Tỷ lệ Chiết khấu: Khi áp dụng hệ số chiết khấu dòng tiền $r = 10%/\text{năm}$, hành vi tích - xả của hệ thống hồ Sơn La – Hòa Bình ưu tiên tạo thặng dư phát điện trong các năm đầu của chu kỳ quy hoạch, nhưng khi tính đến giá trị hiện tại ròng dài hạn (NPV) 25 năm, sự can thiệp điều tiết của Nhà nước nhằm duy trì mực nước tối thiểu cho hạ du giúp triệt tiêu ngoại ứng tiêu cực do sụt giảm dòng chảy môi trường.
Mực nước hồ (Z)
0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 T11 T12
[VỤ ĐÔNG XUÂN] [GIAI ĐOẠN LŨ] [TÍCH NƯỚC PHÁT ĐIỆN]
Implications đa chiều
- Về mặt lý thuyết: Luận án đặt nền móng cho việc ứng dụng lý thuyết kinh tế thị trường có điều tiết vào ngành nước Việt Nam; thống nhất khung phân tích kinh tế vi mô (hàm cầu, thặng dư tiêu dùng) với cơ học thủy lực lưu vực sông.
- Về mặt phương pháp luận: Thiết lập quy trình chuyển đổi bài toán phân bổ tài nguyên liên thời gian phức tạp thành bài toán phi tuyến giải được trên các phần mềm đại số ma trận hiện đại (Lingo, R), mở ra hướng tiếp cận có thể nhân rộng cho các lưu vực khác như sông Mê Kông, sông Đồng Nai hay sông Vu Gia - Thu Bồn.
- Về mặt ứng dụng thực tiễn: Cung cấp công cụ tính toán lịch xả tối ưu theo tháng và theo đợt cho Trung tâm Điều độ Hệ thống Điện Quốc gia (A0/EVN) và Tổng cục Thủy lợi (MARD), giúp chấm dứt tình trạng căng thẳng đàm phán cấp nước đổ ải vụ Đông Xuân hàng năm.
- Về mặt chính sách: Cung cấp cơ sở khoa học trực tiếp cho việc thực thi Luật Thủy lợi 2017 (chuyển từ cơ chế "thủy lợi phí" hành chính sang "giá dịch vụ thủy lợi" thị trường có sự kiểm soát trần giá của Nhà nước nhằm kiểm soát độc quyền tự nhiên).
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:
- Bất đối xứng thông tin xuyên biên giới: Chưa thể tích hợp đầy đủ số liệu vận hành thời gian thực của 62 nhà máy thủy điện tại phần thượng nguồn LVSHTB thuộc lãnh thổ Trung Quốc do rào cản chia sẻ dữ liệu quốc tế; các kịch bản dòng chảy đến biên giới phải dựa trên chuỗi quan trắc tại trạm nhập lưu Việt Nam.
- Phạm vi đối tượng sử dụng nước: Dù khung lý thuyết bao quát đa ngành, các mô hình chạy số số trị tập trung sâu vào hai ngành sử dụng nước chủ lực là tưới nông nghiệp và phát điện thủy điện; các ngành như giao thông thủy, thủy sản và hấp thụ ô nhiễm môi trường mới chỉ được lồng ghép dưới dạng các ràng buộc bất đẳng thức biên.
- Đơn giản hóa thủy lực mạng lưới hạ du: Để đảm bảo khả năng giải của bài toán tối ưu phi tuyến quy mô lớn, mô hình sử dụng quan hệ cân bằng nút và truyền triều dạng giản lược thay vì ghép nối trực tiếp mô hình thủy lực 2D phi tĩnh chi tiết.
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Hướng 1: Ứng dụng Trí tuệ Nhân tạo (AI), Mạng Nơ-ron Tích chập (CNN) và Học tăng cường sâu (Deep Reinforcement Learning) để giải bài toán tối ưu hóa điều phối thời gian thực (real-time adaptive dispatching) theo giờ/phút.
- Hướng 2: Mở rộng bài toán phân bổ tối ưu đa quốc gia dựa trên Lý thuyết Trò chơi (Game Theory - Nash Equilibrium) cho quản trị nguồn nước xuyên biên giới sông Hồng giữa Việt Nam và Trung Quốc.
- Hướng 3: Tích hợp mô hình kinh tế thủy văn với các kịch bản biến đổi khí hậu toàn cầu mới nhất (IPCC AR6 - SSP Scenarios) và biến động sử dụng đất lưu vực.
- Hướng 4: Lượng hóa chi tiết giá trị dịch vụ hệ sinh thái (Ecosystem Services Valuation) và dòng chảy môi trường (Environmental Flows) vào hàm mục tiêu kinh tế.
Tác động và ảnh hưởng
- Tác động học thuật: Luận án mở ra nhánh nghiên cứu Kinh tế Thủy văn Tổng hợp (Hydro-Economics) tại Việt Nam, đóng vai trò tài liệu tham khảo chuẩn mực cho các chương trình đào tạo tiến sĩ, thạc sĩ ngành Kỹ thuật Tài nguyên nước, Kinh tế Tài nguyên Môi trường tại Đại học Thủy lợi, Đại học Quốc gia Hà Nội và Viện Khoa học Thủy lợi Việt Nam. Ước tính công trình có tiềm năng tạo ra hàng trăm trích dẫn trong các nghiên cứu quản trị lưu vực sông tại Đông Nam Á.
- Chuyển đổi ngành Thủy lợi và Năng lượng: Cung cấp luận cứ định lượng giúp Tập đoàn Điện lực Việt Nam (EVN) và các Công ty TNHH MTV Khai thác Công trình Thủy lợi tối ưu hóa kế hoạch huy động nguồn lực, giảm thiểu lãng phí hàng trăm triệu mét khối nước mỗi năm trong các đợt điều tiết nước tăng cường.
- Ảnh hưởng chính sách quốc gia: Tác giả đã trực tiếp đóng góp kết quả nghiên cứu vào các diễn đàn chính sách cấp cao: Hội thảo "Định giá tài nguyên nước" do Ngân hàng Thế giới (World Bank) tài trợ, phản biện Đề tài "Kế toán Tài nguyên Nước" của Viện Chiến lược, Chính sách Tài nguyên và Môi trường (MONRE), và hoàn thiện các quy định điều tiết liên hồ chứa theo Quyết định số 1287/QĐ-TTg của Thủ tướng Chính phủ.
- Lợi ích kinh tế - xã hội: Giúp tăng sản lượng điện hàng năm hàng trăm tỷ đồng, đảm bảo an ninh lương thực cho hàng triệu hộ nông dân đồng bằng sông Hồng, đồng thời giảm thiểu rủi ro ngập úng mùa lũ và thiếu nước mùa khô.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận mẫu mực về mô hình hóa tối ưu hóa động phi tuyến ngẫu nhiên, giải quyết khoảng trống kết hợp giữa khoa học kỹ thuật và kinh tế lượng.
- Nhà hoạch định chính sách (Bộ Nông nghiệp & PTNT, Bộ Tài nguyên & Môi trường, Bộ Công Thương): Sở hữu công cụ khoa học vững chắc để ban hành biểu giá dịch vụ nước, quy trình vận hành liên hồ chứa đa mục tiêu và cơ chế giám sát độc quyền tự nhiên.
- Doanh nghiệp khai thác công trình thủy lợi và Thủy điện (EVN, A0, các Công ty Khai thác Thủy nông): Nâng cao hiệu suất sinh lời trên mỗi mét khối nước trữ, minh bạch hóa trách nhiệm và chi phí cơ hội trong việc phối hợp cấp nước hạ du.
- Cộng đồng dân cư lưu vực sông Hồng - Thái Bình: Thụ hưởng nguồn nước tưới ổn định, giảm thiểu rủi ro thiên tai lũ lụt và nâng cao chất lượng môi trường số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 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 (Natural Monopoly Theory) của Baumol, Panzar và Willig (1982) vào bài toán Quy hoạch và Quản lý Tài nguyên Nước cấp lưu vực sông tại Việt Nam. Luận án đã mở rộng lý thuyết này bằng cách xây dựng quy trình toán học kiểm định đặc tính cộng dưới chi phí (subadditivity) từ nghiệm tối ưu của mô hình tối ưu hóa động, xác lập rằng hệ thống hồ chứa đa mục tiêu vận hành theo cơ chế độc quyền tự nhiên yếu, từ đó luận giải tính tất yếu của việc Nhà nước phải can thiệp điều tiết giá và hạn mức phân bổ thay vì để thị trường tự do định đoạt.
2. Đổi mới phương pháp luận của luận án vượt trội hơn các nghiên cứu trước đây như thế nào?
So với các nghiên cứu trước đây tại Việt Nam (như Hà Văn Khối và Lê Bảo Trung, 2003 chỉ áp dụng quy hoạch động 2 chiều thuần túy kỹ thuật; Viện Quy hoạch Thủy lợi, 2005 dùng GAMS tĩnh; Claudia Ringler và Nguyễn Vũ Huy, 2004 tại sông Đồng Nai không có cấu phần động học), luận án của Đào Văn Khiêm đã tạo bước nhảy vọt:
- Mở rộng hệ thống mô hình hóa Aquarius bằng cách bổ sung hàm chi phí phi tuyến và thuật toán ngẫu nhiên.
- Sử dụng phần mềm Lingo 16/17 kết hợp ngôn ngữ lập trình thống kê R và Stata 13 để giải bài toán quy hoạch phi tuyến ngẫu nhiên (SLQP) với chuỗi thời gian kéo dài 21–25 năm, điều mà các công cụ trước đó không thực hiện được do bùng nổ số chiều tính toán.
3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là hiện tượng suy giảm tính độc quyền tự nhiên cục bộ theo mùa. Trong khi hầu hết các tháng trong năm hệ thống hồ chứa Núi Cốc và Sơn La – Hòa Bình thể hiện đặc tính độc quyền tự nhiên mạnh (chi phí biên giảm, tính cộng dưới chi phí thể hiện rõ), thì vào các tháng kiệt đỉnh điểm (tháng 3, tháng 4), tính cộng dưới chi phí bị suy yếu nghiêm trọng. Nguyên nhân kinh tế là do chi phí cơ hội của việc giữ nước phát điện tăng vọt khi nhu cầu tưới hạ du đạt đỉnh, tạo ra sự cạnh tranh gay gắt giữa hai đối tượng dùng nước, làm biến dạng cấu trúc chi phí truyền thống của hồ chứa.
4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) không?
Có. Luận án xây dựng quy trình nhân rộng 5 bước tường minh:
- Xác định sơ đồ cấu trúc hệ thống nguồn - cầu và các nút điều tiết lưu vực;
- Ước lượng các hàm cầu kinh tế động (hàm sản xuất mùa vụ, CVM, TCM);
- Thiết lập hàm mục tiêu thặng dư kinh tế ròng và ma trận ràng buộc kỹ thuật - thể chế;
- Mã hóa mô hình đại số trên phần mềm Lingo/GAMS và nạp chuỗi dữ liệu ngẫu nhiên (ARIMA);
- Chạy mô hình tối ưu hóa phi tuyến, kiểm tra tính ổn định nghiệm qua mô phỏng Monte Carlo và phân tích cấu trúc chi phí subadditivity.
5. Định hướng chương trình nghiên cứu 10 năm được phác thảo ra sao?
Chương trình nghiên cứu 10 năm được định hình theo 3 giai đoạn:
- Giai đoạn 1 (1-3 năm): Mở rộng khung mô hình cho toàn bộ 4 hồ chứa lớn nhất LVSHTB (Sơn La, Hòa Bình, Thác Bà, Tuyên Quang) với phân giải thời gian ngày/giờ phục vụ thị trường phát điện cạnh tranh.
- Giai đoạn 2 (4-6 năm): Tích hợp công cụ mô phỏng chất lượng nước và dòng chảy sinh thái vào bài toán tối ưu hóa đa mục tiêu; xây dựng nền tảng giao dịch hạn ngạch nước (water permit trading).
- Giai đoạn 3 (7-10 năm): Ứng dụng Trí tuệ Nhân tạo phân tán (Multi-Agent Systems) và Dữ liệu lớn viễn thám để thiết lập Hệ thống Hỗ trợ Ra Quyết định Thông minh (IDSS) cấp quốc gia trong điều kiện biến đổi khí hậu cực đoan.
Kết luận
Luận án Tiến sĩ Kỹ thuật của tác giả Đào Văn Khiêm đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp cốt lõi:
- Thiết lập thành công mô hình toán học tối ưu hóa động ngẫu nhiên cho bài toán phân bổ hiệu quả kinh tế tài nguyên nước tại hai tiểu lưu vực trọng điểm (hệ thống Núi Cốc và liên hệ thống Sơn La – Hòa Bình) thuộc LVSHTB.
- Xây dựng phương pháp ước lượng hàm cầu động đối với nước tưới nông nghiệp và phát điện thủy điện, chuyển đổi thành công các giá trị kinh tế vi mô thành hàm mục tiêu thặng dư phúc lợi xã hội ròng.
- Phát triển và mở rộng khung mô hình Aquarius, bổ sung các hàm cấu trúc chi phí phi tuyến (hàm mũ, hàm bậc ba) và xử lý bất định thủy văn bằng chuỗi thời gian ARIMA.
- Tiên phong chứng minh thực nghiệm cấu trúc Độc quyền tự nhiên yếu của hệ thống liên hồ chứa đa mục tiêu tại Việt Nam thông qua kiểm định tính cộng dưới chi phí (subadditivity).
- Đưa ra các giải pháp phân bổ nước tối ưu liên thời gian, chứng minh khả năng gia tăng tổng lợi ích kinh tế ròng từ 18,5% đến 24,2% so với quy trình vận hành tĩnh truyền thống.
- Cung cấp luận cứ khoa học và thực tiễn trực tiếp phục vụ công tác điều hành liên hồ chứa của Chính phủ, EVN, MARD và hỗ trợ thể chế hóa Luật Thủy lợi 2017 theo định hướng kinh tế thị trường có sự điều tiết của Nhà nước.
Công trình đánh dấu bước chuyển dịch mô thức (paradigm shift) từ quản lý cấp phát tài nguyên nước mang tính kỹ thuật hành chính thụ động sang quản trị kinh tế - thủy văn tổng hợp hiện đại, mở ra chân trời nghiên cứu liên ngành bền vững cho khoa học tài nguyên nước Việt Nam và khu vực.
