- Luật Tài nguyên nước số 17/2012/QH13 ngày 21/6/2012.
- Nghị định 143/2003/NĐ-CP ngày 28 tháng 11 năm 2003 quy định chi tiết thi hành một số điều của Pháp lệnh khai thác và bảo vệ công trình thuỷ lợi.
- Nghị định 112/2008/NĐ-CP ngày 20 tháng 10 năm 2008 về quản lý, khai thác tổng hợp tài nguyên và môi trường các hồ chứa thuỷđiện, thuỷ lợi.
- Nghịđịnh 72/NĐ-CP ngày 07 tháng 02 năm 2007 về quản lý an toàn đập.
- Chỉ thị của Chính phủ tại văn bản số 21/CT-TTg ngày 14 tháng 10 năm 2013 về việc tăng cường quản lý đảm bảo an toàn hồ chứa nước.
- Thông tư 33/2008/TT-BNN ngày 04 tháng 02 năm 2008 về việc hướng dẫn một sốđiều thuộc NghịĐịnh 72/NĐ-CP.
- Chính sách hoạt động an toàn đập OP/BP 4.37 của Ngân hàng thế giới.
- Các quy chuẩn, tiêu chuẩn hiện hành khác có liên quan.
Các Quy chuẩn, tiêu chuẩn áp dụng để lập báo cáo an toàn đập
- Quy chuẩn kỹ thuật quốc gia QCVN 04-01:2010/BNNPTNT về Thành phần, nội dung lập Báo cáo đầu tư, Dự án đầu tư và Báo cáo kinh tế kỹ thuật công trình thủy lợi.
- Quy chuẩn kỹ thuật quốc gia QCVN 04-05:2012/BNNPTNT – Các quy định chủ yếu về thiết kế công trình thủy lợi.
- Tiêu chuẩn Việt nam TCVN 8216 : 2009- Thiết kế đập đất đầm nén.
- TCVN 8421:2010 Công trình thủy lợi – Tải trọng và lực tác dụng lên công trình do sóng và tàu.
- Các tiêu chuẩn Quy chuẩn hiện hành.
3.2.Tính toán kiểm tra an toàn công trình đầu mối.
3.2.1 Các thông số hồ chứa và công trình đầu mối theo điều kiện mới Tài liệu mực nước
- MNDBT= 33,2 m - MNLTK = 37,21 m - MNLKT = 38,15 m
- ∇chânđập = 1,4 m Tài liệu gió
- Hướng gió tính toán bất lợi nhất: Hướng gió vuông góc với trục đập.
- Vận tốc gió trung bình tại công trình: v = 1,5 m/s
+ Đà gió thổi ứng với MNDBT: L1= 368,56m = 0,369km + Đà gió thổi ứng với MNLTK: L2= 386,19m = 0,386km - Thời gian tác dụng của gió: t = 6 giờ
Cấp công trình
Tra theo QCVN 04-05-2012, cấp công trình là cấp I 3.2.2 Tính toán xác định cao trình đỉnh đập
Theo [2] cao trình đỉnh đập là cao trình lớn nhất xác định theo các trường hợp sau:
Trường hợp 1: Tính với MNDBT ở thượng lưu. Tốc độ gió tính toán được lấy theo tần suất p = 2%.
∇1 = MNDBT + Dh + hsl + a1
Trường hợp 2: Tính với MNLTK ở thượng lưu. Tốc độ gió tính toán lấy theo tần suất p = 25%.
∇2 = MNLTK+ Dh + hsl + a2 Trường hợp 3: Tính với MNLKT ở thượng lưu
∇3 = MNLKT + a3
Trong đó: a - chiều cao an toàn, xác định theo [2], công trình cấp I a1 = 1,2 m ứng với MNDBT
a2 = 1m ứng với MNLTK a3 = 0,3 m ứng với MNLKT
a. Xác định các tham số trong hai công thức trên
• Tính toán chiều cao nước dềnh do gió ∆h
∆h:Là chiều cao nước dềnh lên ứng với tốc độ gió được xác định theo công thức:
Trong đó :
- Kw: hệ số phụ thuộc vận tốc gió - V (m/s): Là vận tốc gió tính toán - L (m): Đà gió
- d (m): Cột nước trước đập
- a (độ): Góc kẹp giữa hướng gió tính toán với hướng vuông góc trục đập, α=0
• Xác định chiều cao sóng leo hsl
Theo [3], công thức chiều cao sóng dềnh được tính toán ứng với mức bảo đảm p = 1%.
Giả thiết sóng tạo dòng vùng nước sâu : H/λ > 0,5 Chiều cao sóng leo được xác định như sau:
hsl1% = Kr. Kp . Ksp. Krun .Ka. h1%
α cos ). 5 , 0 .(
. .
2
h d
g
L K V
h w
∆
= +
∆
Trong đó:
+ Kr, Kp: là các hệ số phụ thuộc vào đặc trưng lớp gia cố mái và độ nhám tương đối trên mái thượng lưu đập. Đập thiết kế có hệ số mái thượng lưu phía trên m = 3,25; Mái đập không được gia cố.
+ Ksp: là hệ số phụ thuộc vào vận tốc gió và hệ số mái dốc đập thượng lưu. + Krun: là hệ số phụ thuộc vào tỷ số λs/ hs1% và mái dốc đập.
- ls (m): Chiều dài trung bình của bước sóng, xác định theo công thức dưới đây:
Trong đó: với trị số .
+ Ka : là hệ số phụ thuộc vào góc kẹp giữa hướng gió tính toán với trục hồ.
• Xác định chiều cao hs1%
Chiều cao của sóng leo vùng nước sâu d > 0.5lsvới mức đảm bảo 1% được xác định theo công thức sau:
Trong đó:
- K1% : là hệ số có mức bảo đảm 1% phụ thuộc vào .2
v L
g tra theo đồ thị hình A2.
-hs Chiều cao trung bình của sóng
(m).
Trị số phụ thuộc vào
2 2
; . .
v d g v
L
f g được xác định theo đồ thị hình A1.
Kết quả tính toán được lập thành bảng sau:
π λ τ
.2
2 s =g
= v
g g v τ
τ .
= 2 .2 . , .
v d g v
L f g V gτ
s
s K h
h1% = 1%.
= 2 .2 v
h x g g hs v
2
. v
h g
Bảng 3.1 Kết quả tính sóng
Thông số Đơn vị MNDBT MNLTK
(P=0.5%)
∇ Ch©n ®Ëp TL m 1,40 1,40
d m 31,60 35,81
L m 1200,00 1500,00
V m/s 1,50 1,50
α độ 0,00 0,00
Kw*106 - 2,00 2,00
∆h m 0,00002 0,000019
g . L / V2 - 5232,00 6540,00
g .h / V2 - 0,0800 0,0880
g .τ / V - 3,90 4,20
hs m 0,02 0,02
τ s 0,60 0,64
λs m 0,56 0,64
d /λs - 56,89 55,58
gd/V2 - 137,776 156,132
k1% - 2,06 2,10
hs1% m 0,038 0,042
λs / hs1% - 14,70 15,20
Kr - 1,00 1,00
Kp - 0,90 0,90
Ksp - 1,10 1,10
Krun - 1,50 1,50
Kα - 1,00 1,00
hsl1% m 0,06 0,06
a m 1,20 1,00
Bảng 3.2 Kết quả tính toán cao trình đỉnh đập
Truờng hợp ∇MNTL (m) ∆h (m) hsl1% (m) a (m) ∇Đỉnh đập (m)
MNDBT 33,20 0,000017 0,056 1,20 34.26
MNLTK 37,21 0,000019 0,063 1,00 38.27
MNLKT 38,15 0,30 38,45
Kết Luận: Chọn cao trình đỉnh đập thiết kế
∇Đỉnh đập = +38.45 (m)<∇Đỉnh đập hiện trạng = +39.40 (m) Như vậy, giữ nguyên cao trình đỉnh đập hiện tại là +39,4m.
Nhận xét:Ở đây nhận thấy sau khi tính toán lại với các nhu cầu sử dụng nước đều tăng lên đáng kể so với nhu cầu hiện tại, tuy nhiên kết quả tính toán cao trình đỉnh đập lại cho thấy thấp hơn cao trình đỉnh đập hiện tại gần 1m, lý do là vì khi thiết kế hồ chứa nước Sông Mực thì có tính thêm cả dung tích tưới đủ cho cả hạ lưu của hồ Yên Mỹ sau này. Sau khi hồ Yên Mỹ xây dựng, đã giảm tải phần dung tích này cho hồ chứa Sông Mực, vì thế tuy các nhu cầu dùng nước đều lớn hơn so với hiện tại nhưng tính toán ra cao trình thấp hơn cao trình hiện trạng là hợp lý.
3.2.3 Lựa chọn giải pháp đảm bảo ổn định công trình
Hồ chứa nước Sông Mực được xây dựng năm 1977 và đưa vào khai thác năm 1981. Qua quá trình sử dụng cụm đầu mối công trình vẫn hoạt động bình thường.
Bề mặt đập qua quá trình sử dụng mặt đường bị hỏng, nay kiến nghị rải nhựa bề mặt đập dày 20cm, đảm bảo giao thông đi lại trên mặt đập. Mái thượng lưu đập được gia cố bằng bê tông cốt thép từ cơ ở cao trình +28,5m đến cao trình đỉnh đập. Mái hạ lưu giữ nguyên hiện trạng vẫn đang hoạt động bình thường.
Tràn xả lũ: Từ kết quả điều tiết lũ tác giả nhận thấy rằng bề rộng tràn nhỏ (∑B tràn = 8m) nên lưu lượng xả qua tràn khá nhỏ so với lưu lượng lũ đến, nhưng vẫn đảm bảo an toàn ứng với lũ thiết kế và lũ kiểm tra; hồ chứa vì thế không cần mở rộng bề rộng ngưỡng tràn.Tuy nhiên, MNDBT mới tăng lên 20cm so với MNDBT hiện trạng, vì thế đểđảm bảo giữ mực nước mới ở cao trình+33,2 thì ta tiến hành hàn thêm bản tap bằng thép ở đỉnh cửa van để nâng chiều cao của cửa van lên 20cm, cao trình đỉnh cửa van mới là +33,2m; đỡ phải thay mới cửa van tốn kém.
3.2.4 Kiểm tra ổn định công trình 3.2.3.1 Tài liệu tính toán
a. Mặt cắt tính toán
- Mặt cắt đập tại vị trí lòng sông, chiều cao đập lớn nhất.
b. Số liệu cơ bản
- Tài liệu vềđịa chất thủy văn và địa chất công trình.
- Tài liệu về công trình:
Bảng 3.3. Thông số hồ chứa
Hồ chứa Đập đất
Flv (km2)
MNDBT (m)
MNC (m)
MNLTK (m)
W (106m3 )
Ft (ha)
∇δ (m)
Hd (m)
Ld (m) 236 33,2 18 37,21 205,11 11335 39,4 38 470
- Các tiêu chuẩn thiết kế:
+ Các quy định chủ yếu về thiết kế: QCVN 04 – 05 : 2012/BNNPTNT.
+ Thiết kế đập đất đầm nén – TCVN 8216 : 2009.
+ Quy phạm TCVN 4253- 86: Nền các công trình thủy công.
+ Phần mềm địa kỹ thuật chuyên dụng có bản quyền Geostudio 2004.
Bảng 3.4. Hệ số thấm của các lớp đất
TT Loạt đất K (m/s)
1 Lớp 1 2,15 x10-5
2 Lớp 2 1,13 x10-6
3 Lớp 3 5,4 x10-5
4 Lớp 4 5,83x10-6
c. Trường hợp tính toán
- MNDBT ở thượng lưu +33,20 m, hạ lưu không có nước.
- MNLTK ở thượng lưu +37,21 m, mực nước tương ứng ở hạ lưu +32,0 m.
- MNLTK ở thượng lưu +38,15 m, mực nước tương ứng ở hạ lưu +32,0 m.
d. Phương pháp tính toán
- Phần mềm chuyên dụng Geostudio 2004 được sử dụng để kiểm tra ổn định thấm và ổn định trượt.
3.2.3.2 Tính toán ổn định thấm
a. Nguyên lý cơ bản của bài toán thấm
- Bài toán tính toán thấm tuân theo định luật của Darcy: q = k.i Trong đó:
+ q là lưu lượng thấm đơn vị + k là hệ số thấm
+ i là độ dốc thủy lực
- Tính toán thấm qua đập đất theo phương pháp phần tử hữu hạn xác định đường bão hòa ổn định, trường vận tốc và gradient thấm. Phương pháp này cho kết quả chính xác đặc biệt là trong trường hợp các vùng vật liệu phức tạp.
-Chương trình tính toán áp dụng là SEEP/W. Đây là một phần mềm PTHH được dùng để mô hình hoá các chuyển động và phân bố các áp lực
nước kẽ rỗng trong các môi trường rỗng như đất và đá - sản phẩm của GEO-Slope ternational Ltd – Canada.
Từ đó xác định tổng lưu lượng thấm nhằm kiểm tra lượng tổn thất do thấm có nằm trong phạm vi cho phép và kiểm tra độ bền thấm thân đập và nền.
- Tổng lưu lượng thấm: lưu lượng thấm qua thân đập được giới hạn như sau:
Wt<[W]cp= 1%Whồ Trong đó:
+ Wt là tổng lượng nước thấm trong thời đoạn 1 tháng, Wt = t.Qt
+ Qt là tổng lưu lượng thấm qua thân đập.
- Kiểm tra độ bền thấm cho thân đập và nền: Jkđ ≤[ ]J đ.
Trong đó: - [J]đ tra trong bảng 5 TCVN 8216:2009, [J]đ = 0,85.
b. Kết quả tính toán
Sau khi cập nhật các tài liệu khảo vào các mặt cắt thiết kế đập đất, tiến hành tính toán ổn định trượt mái và điều chỉnh lại mặt cắt thiết kế nếu chưa hợp lý. Qua kết quả tính toán kết luận về ổn định thấm:
- Các mặt cắt tính toán được chọn đều bất lợi về mặt ổn định thấm.
- Ở trường hợp mực nước dâng bình thường, ta thấy dòng bão hòa trong thân đập đổ vào đống đá tiêu nước, công trình làm việc bình thường; phù hợp với thực tế làm việc của công trình.
- Ở trường hợp mực nước lũ thiết kế và lũ kiểm tra có xuất hiện dòng thấm ra mái hạ lưu, gây nguy hiểm cho công trình. Kiến nghị áp mái tiêu nước hạ lưu trong phạm vi lòng sông cũ với chiều cao 2m tính từ đỉnh đống đá tiêu nước. (Kết quả tính toán xem phụ lục tính thấm 02)
- Về kiểm tra điều kiện Gradient trong thân đập: Jmax =0,4929 < [J] = 0,85; đảm bảo bé hơn Gradient giới hạn cho phép.
3.2.3.3 Tính toán ổn định mái đập
a. Nguyên lý cơ bản của bài toán tính ổn định
- Tính toán ổn định mái đập theo phương pháp Bishop, là một phương pháp được dùng phổ biến hiện nay.
- Chương trình tính toán áp dụng là SLOPE/W. Đây là một phần mền ứng dụng lý thuyết cân bằng giới hạn để xác định hệ số ổn định của các mái đất đá, được lập bởi GEO-Slope ternational Ltd–Canada.
-Hệ sốổn định cho phép ứng với các trường hợp tính toán.
Theo QCVN 04-05:2012/BNNPTNT, Trong đó:
K - hệ sốổn định mái của công trình.
nc - hệ số tổ hợp tải trọng, đối với tính toán theo TTGH thứ nhất.
nc =1,00 - đối với tổ hợp tải trọng cơ bản.
nc =0,90 - đối với tổ hợp tải trọng đặc biệt.
Ntt- tải trọng tính toán tổng quát.
R - sức chịu tải tổng quát.
m - hệ số điều kiện làm việc, tra phụ lục B thì đối với các mái dốc tự nhiên và nhân tạo m=1,00.
kn- hệ số bảo đảm được xét theo quy mô, nhiệm vụ công trình. Đối với các công trình cấp I lấy kn = 1,35.
- Hệ sốổn định mái cho phép của công trình [K] theo [2]
[K] = 1,35 - đối với tổ hợp tải trọng cơ bản.
[K] = 1,15 - đối với tổ hợp tải trọng đặc biệt.
b. Kết quả tính toán
Bảng 3.5. Kết quả tính toán ổn định mái đập
m k . ] n K N [
K R c n
tt
=
≥
=
Trường hợp tính toán Kminmin
[K]cp
(Tổ hợp cơ bản) MNDBT = +33,2m;
Hạ lưu không có nước 1,433
1,35 MNLTK = +37,21m;
Hạ lưu 32,00 1,363 MNLKT = +38,15m;
Hạ lưu 32,00 1,254 1,15
c. Kết luận, kiến nghị
- Các mặt cắt tính toán được chọn đều bất lợi về mặt ổn định trượt mái.
- Ởtrường hợp mực nước thượng lưu là mực nước dâng bình thường, hệ sốổn định mái k = 1,433 > [k] = 1,35 đảm bảo ổn định trượt của đập.
- Ở trường hợp mực nước thượng lưu là MNLTK thì hệ số k = 1,363 > [k] = 1,35; nên trường hợp này mái đập an toàn.
- Ở trường hợp mực nước thượng lưu là MNLKT thì hệ số k = 1,254 > [k] = 1,15; trường hợp này mái đập an toàn.
Như vậy, ta thấy các TH tính toán đều cho kết quả hệ số ổn định lớn hơn hệ số ổn định cho phép, công trình đảm bảo làm việc bình thường.
Kết luận Chương 3
Sau khi tính toán hồ chứa theo dung tích của nhiệm vụ mới, tác giả đã đưa ra được các thông số mực nước mới của hồ chứa. Từ đó tác giả tiến hành tính toán lại cao trình đỉnh đập, kết quả cho thấy cao trình đỉnh đập tính theo nhiệm vụ mới là +38,45m, thấp hơn cao trình đỉnh đập hiện tại là +39,4m, do nhiệm vụ mới của hồ chứa nước Sông Mực đã không còn tính đến dung tích mà hiện nay hồ chứa Yên Mỹ đã đảm nhận sau này; vì thế hồ Sông Mực đảm bảo đập an toàn trong mọi trường hợp về mực nước. Giải pháp đưa ra không phải gia cố, tôn cao đỉnh đập, chỉ rải nhựa lại bề mặt đỉnh đập dày 20cm.
Tràn xả lũ vẫn đảm bảo xả lũ với lưu lượng tính toán theo nhiệm vụ mới nên không cần mở rộng tràn mà chỉ sửa chữa cánh cửa van cung bằng cách hàn bản tab bằng thép trên đỉnh cửa van lên cao 20cm, với kích thước cửa van là nxbxh = 2x4x5,2m; đảm bảo hồ chứa tích nước với MNDBT theo nhiệm vụ mới là +33,2m.
Tác giả cũng đã kiểm tra lại điều kiện ổn định công trình về trượt, thấm ứng với các trường hợp mực nước hồ mới. Kết quả xác định là theo điều kiện mới và theo tiêu chuẩn hiện hành thì công trình là Cấp I; Kết quả cho thấy công trình đảm bảo ổn định về trượt khi hồ chứa làm việc theo điều kiện mới trong mọi trường hợp tính toán. Đối với ổn định thấm, kết quả tính toán cho thấy có xuất hiện dòng thấm ra mái với TH mực nước lũ thiết kế và mực nước lũ kiểm tra, tác giả kiến nghị làm áp mái tiêu nước hạ lưu trong phạm vi lòng sông cũ, với chiều cao 2m tính từ đỉnh đống đá tiêu nước. Như vậy với hồ Sông Mực để đáp ứng phục vụ theo nhiệm vụ mới thì cần kết hợp 2 giải pháp sửa chữa nâng cao cao trình cửa van lên 20cm và gia cố áp trúc mái hạ lưu thêm 2m tính từ đỉnh đống đá tiêu nước.
KẾT LUẬN VÀ KIẾN NGHỊ
1. Những kết quảđạt được của luận văn
Luận văn đã kế thừa các kết quả nghiên cứu trước đây trong xem xét đánh giá, tính toán dự báo nhu cầu sử dụng nước, đánh giá tiềm năng nguồn nước; xác định các vấn đề trong khai thác sử dụng tài nguyên nước; dự báo xu thế diễn biến tài nguyên nước trong kỳ quy hoạch trên địa vùng. Kết quả chính đạt được bao gồm:
- Đánh giá được đặc điểm hiện trạng các hồ chứa nước và hệ thống thủy lợi; tình hình quản lý khai thác hệ thống công trình thủy lợi tỉnh Thanh Hóa.
- Đánh giá hiện trạng hệ thống công trình hồ chứa nước Sông Mực, xác định được nhu cầu sử dụng nước hiện tại của công trình.
- Đã tính toán, xác định được nhu cầu sử dụng nước mới của các ngành nông nghiệp, sinh hoạt, thủy sản, công nghiệp… theo hướng phát triển kinh tế xã hội phụ thuộc trực tiếp từ nguồn nước hồ chứa sông Mực. Qua đó nhận thấy rằng đểđáp ứng được nhiệm vụ của hồ chứa theo nhiệm vụ mới thì cần thay đổi MNDBT từ +33,0 đến MNDBT đến 33,2m; MNLTK +37,21m;
MNLKT +38,15m.
- Với các mực nước trên, tác giả đã tiến hành tính toán lại xác định cao trình đỉnh đập của hồ Sông Mực theo nhiệm vụ mới được +38,45m nhỏ hơn cao trình đỉnh đập cũ +39,4m: Do đó không cần nâng cao cao trình đỉnh đập.
- Với các mực nước đã tính toán được tác giả đề ra giải pháp đó là kết hợp 2 phương pháp là sửa chữa nâng cao trình đỉnh cửa van cung lên 20cm và đắp áp trúc mái hạ lưu thêm 2m.
- Đánh giá kiểm tra an toàn hồ chứa nước sông Mực khi nâng cao dung tích hồ chứa theo nhiệm vụ mới.
2. Một sốđiểm tồn tại và hướng tiếp tục nghiên cứu 2.1. Một số điểm còn tồn tại
- Luận văn chưa tính toán được hết các nhu cầu dùng nước cụ thể cho một sốtrường hợp như phát điện, du lịch…
- Chưa áp dụng tính toán cụ thể các giải pháp khác nhau để từ đó so sánh, lựa chọn được phương án tối ưu nhất cho công trình.
- Tồn tại nữa của luận văn đó là tính toán điều tiết, cao trình, ổn định … của đập đều dựa theo các điều kiện địa chất, địa hình cũ khi thiết kế hồ chứa mà không có được các tài liệu hiện tại, khi mà sau một thời gian dài hoạt động hồ chứa đã bị bồi lắng, xói lở, thay đổi địa hình địa mạo nhiều. Tuy nhiên, trong phạm vi luận văn không có điều kiện tiến hành khảo sát lại được.
2.2. Hướng tiếp tục nghiên cứu của luận văn
- Tính toán chi tiết thêm các nhu cầu dùng nước của hồ chứa nước sông Mực, đảm bảo cấp nước cho chiến lược phát triển kinh tế xã hội của địa phương đến năm 2020, định hướng đến năm 2030.
- Áp dụng các phương án khác nhau để gia cố, đảm bảo an toàn cho công trình khi hồ chứa tích nước theo dung tích mới.
TÀI LIỆU THAM KHẢO
1. Bộ Nông nghiệp và PTNT (2002), Chương trình đảm bảo an toàn các hồ chứa nước, Nhà xuất bản xây dựng, Hà Nội.
2. Bộ Nông nghiệp và PTNT (2012), QCVN 04-05: 2012/BNNPTNT, Công trình thủy lợi - Các quy định chủ yếu về thiết kế.
3. Bộ Nông nghiệp và PTNT (2009), Tiêu chuẩn thiết kế đập đất đầm nén 8216-2009.
4. Bộ Nông nghiệp và PTNT (2012), TCVN4253 : 2012 Công trình thủy lợi – Nền các công trình thủy công – yêu cầu thiết kế, NXB Nông Nghiệp.
5. Bộ Thủy Lợi (1977), Quy phạm tính toán đặc trưng thủy văn thiết kế QP.TL. C-6-77.
6. Đinh Quang Dương (2014),Thủy lợi Thanh Hóa, Thanh Hóa.
7. Phạm Ngọc Quý (2008), Tràn sự cố trong đầu mối hồ chứa nước, Nhà xuất bản Nông nghiệp.
8. Nguyễn Xuân Trường (1972),Thiết kế đập đất, NXB Khoa học kỹ thuật.
9. Ngô Trí Viềng, Nguyễn Chiến, Nguyễn Văn Mạo, Nguyễn Văn Hạnh, Nguyễn Cảnh Thái (2004),Giáo trình thuỷ công, Trường Đại học Thuỷ lợi, Nhà xuất bản xây dựng, Hà Nội.
10. Viện quy hoạch thủy lợi (2012), Quy hoạch tổng thể thủy lợi Thanh Hóa đến 2020 và định hướng đến 2030
11. Viện kỹ thuật công trình (2014), Lập quy trình vận hành điều tiết hồ chứa nước Sông Mực.