【正文】 8+770 324 667 62+500 364 575 80+000 883 750 63+000 414 582 80+780 358 850 63+500 458 643 81+150 802 1142 64+000 328 484 82+000 198 567 64+500 292 633 82+000 330 580 引黃濟青攔河閘 65+000 345 510 82+500 389 501 65+500 299 377 83+000 483 682 13 樁號 主河槽 寬度（ m） 堤距 （ m） 注釋 樁號 主河槽 寬度（ m） 堤距 （ m） 注釋 83+500 444 665 91+500 519 1050 84+000 440 616 92+000 360 630 84+420 435 734 93+000 326 550 南膠萊河匯入 84+575 432 582 濟青高速公路 93+500 248 920 河榮二村橋 85+000 506 758 93+680 183 膠濟鐵路橋 85+500 613 755 94+500 280 86+000 710 896 95+000 286 86+500 498 768 95+750 358 87+000 362 561 膠濟電氣化鐵路橋 96+000 316 204國道 87+500 417 538 96+500 238 88+000 388 601 97+000 288 88+500 403 786 97+500 428 89+030 272 551 98+000 366 89+500 366 570 98+820 492 89+985 383 605 99+500 404 90+000 345 589 100+000 327 90+000 408 580 賈疃橡膠壩 100+500 419 90+050 396 589 賈疃橋 101+000 358 850 90+500 212 676 101+500 802 1142 南莊橡膠壩 91+000 226 721 大沽河水文分析 水文氣候特征 大沽河流域地處膠東沿海，屬海洋性氣候。溫度一般在攝氏 22℃ ~38℃之間，流域內(nèi)溫差不大，終年無霜期約 200 天。 流域內(nèi)雨量較充沛，多年平均降水量為 ，徑流量為 108m3（ 1956年～ 2021 年）。降雨量年際變化較大，最大年降雨量為 （ 1964 年），最小年降雨量為 （ 1981 年）。降雨量年 內(nèi)分布不均，約有 80%集中在汛期 14 （ 6～ 9 月），尤以 8 月份最多，約占全年的 56%左右，年降雨量天數(shù)一般不足 90 天。 降雨量的年際變化導(dǎo)致了年徑流量的年際變化，且年徑流量的年際變化更為劇烈，最大年徑流量為 108m3（ 1964 年），最小年徑流量為 108m3（ 1968年），最大年徑流量是最小年徑流量的 328 倍。豐水年與枯水年一般交替出現(xiàn)，也有連豐連枯現(xiàn)象。徑流量年內(nèi)分配也不均勻，汛期洪水暴漲暴落，易形成水災(zāi)，枯水期徑流量很小，甚至斷流。 設(shè)計 洪水 大沽河流域上設(shè)有張家院、產(chǎn)芝、姜家許村、尹 府、葛家埠、嵐西頭、南村、閘子、山角底、大麻灣等多處水文站，南墅、孫受等 25 處雨量站，能滿足流域規(guī)劃治理、設(shè)計等各項工作對資料的需要。 根據(jù)《水利水電工程設(shè)計洪水計算規(guī)范》（ SL4493）第 條規(guī)定“??具有 30 年以上實測或插補延長洪水流量資料，并具有調(diào)查歷史洪水時，應(yīng)采用頻率分析法計算設(shè)計洪水”。南村水文站為大沽河干流控制站， 1951 年設(shè)站，有54 年實測流量資料，并且有調(diào)查歷史洪水資料，本次設(shè)計洪水按規(guī)范規(guī)定采用南村站實測流量資料頻率分析法計算大沽河流域現(xiàn)狀水利工程條件下設(shè)計洪水。 南村水文站洪 水頻率分析 現(xiàn)狀工程條件下南村水文站洪水過程的推求 南村水文站是大沽河干流的控制站，位于平度市南村鎮(zhèn)，控制流域面積3724km2。有 1951～ 2021 年的水文觀測資料，本次對南村站現(xiàn)狀工程條件下的洪水進行頻率分析。南村水文站以上集水面積內(nèi)有產(chǎn)芝、尹府兩座大型水庫及高 15 格莊、勾山、城子、北墅、黃同等五座中型水庫。 考慮到中型水庫無水文觀測資料，且汛期一般按控制運用方案運行，本次南村站頻率分析不對中型水庫的洪水過程進行處理。兩座大型水庫均設(shè)有水文站，自設(shè)立以來均有水文觀測資料，本次對兩座大型水庫的入庫洪水按汛期控 制運用方案重新進行調(diào)洪演算。 由產(chǎn)芝、尹府兩水庫的實測洪水資料，反推出歷年各次洪水的入庫洪水過程，按照本水庫的汛期控制運用方案，起調(diào)水位取當次洪水的實際起漲水位，對每次洪水的入庫洪水過程進行調(diào)洪演算，調(diào)算后的水庫泄流過程再進行河道匯流演算，演算到南村站。根據(jù)實測洪水資料，產(chǎn)芝水庫 1996 年、 1998 年有溢洪，尹府水庫 2021 年有溢洪。南村站實測洪水過程扣除經(jīng)河道匯流演算后的產(chǎn)芝、尹府實際溢洪過程，得出南村站以上產(chǎn)芝、尹府水庫以下區(qū)間洪水過程。南村站以上產(chǎn)芝、尹府水庫以下區(qū)間洪水過程與產(chǎn)芝、尹府水庫重新調(diào)洪演 算后到南村站的洪水過程 疊加 得出南村站現(xiàn)狀工程條件下的洪水過程。采用年最大值選樣方法選取南村站現(xiàn)狀工程條件下的年最大洪峰流量。水庫未建年份的洪峰流量按面積比的 2/3 次方折算求得。 頻率分析 將年最大洪峰流量系列采用數(shù)理統(tǒng)計法進行頻率分析，參數(shù)估計采用矩法計算出一組初值，用經(jīng)驗適線法擬定曲線、確定參數(shù)。 實測系列經(jīng)驗頻率計算公式 1111 ?? ??????? ????? Ln LmN aN aP m 式中： 16 Pm實測系列第 m項的經(jīng)驗頻率； N調(diào)查考證期； a在 N 年中連續(xù)順位的特大洪水項數(shù)； n實測 洪水系列項數(shù)； L實測洪水系列中抽出作特大值處理的項數(shù)； m實測洪水的序位。 均值 X 、變差系數(shù) Cv計算公式 ???????? ???? ? ?? ??aj nLi ij XLn aNXNX 1 11 ? ? ? ? ???????? ??????? ? ?? ??ajnLi ijXXLn aNXXNXCv1212111 式中： Xj歷史調(diào)查大洪水變量（ j=1， 2?? a）； Xi歷史調(diào)查大洪水變量（ i=L+1， L+2?? n）。 根據(jù)年最大洪峰流量系列和上述頻率計算公式及 PⅢ 型頻率曲線適線分析，采用 Cs=，得到南村站現(xiàn)狀工程條件下各頻率洪 峰流量如 表 22。 表 22 南村站現(xiàn)狀工程條件下年最大洪峰流量分析成果表 均值（ m3/s） Cv 各不同頻率設(shè)計值（ m3/s） 計算 適線 20% 10% 5% 3% 2% 1% 1153 1807 2508 3043 3476 4232 南村站設(shè)計洪水成果的合理性分析 將本次計算成果與兩次南村站實測流量資料頻率分析成果比較，比較 情況見表 23。 17 表 23 南村站實測流量頻率計算成果比較表 計算年份 重現(xiàn)期 P=20% P=10% P=5% P=3% P=2% P=1% 1998 2840 3330 3820 2021 1670 2220 2590 2960 3530 本次 1153 1807 2508 3043 3476 4232 成果比較（ %） 與 1998年 與 2021年 通過比較可以看出，本次計算成果，比 1998 年成果偏少 10%左右，比 2021年成果偏大 8%~17%左右。 1998 年進行頻率分析時，沒有將水庫未建年份的資料修正到現(xiàn)狀水利工程條件下，計算所依據(jù)的資料缺乏一致性；而 2021 年計算時將南村站以上的大中型水庫未建年份的資料都進行了修正，而本次計算只考慮了兩座大型水庫未建年份的資料修正，沒有考慮中型水庫的影響，這就造成了兩者 1971 年（北墅水庫修建）以前的資料不一致，本次計算這一 時期系列的均值是 2021 年分析計算時的 倍，由于本段時期屬豐水期，所以造成本次計算成果比 2021 年計算成果偏大，由于中型水庫沒有資料，強制修正容易造成較大誤差，且不考慮中型水庫的影響的設(shè)計成果偏安全，所以可以認為本次的設(shè)計洪水成果比較合理。 再將本次設(shè)計洪水的成果與南村站實測資料進行對照。南村站 1951 年建站，有 54 年實測資料，在 1958 年修建水利工程前的 7 年間，有 5 年最大流量超過2021m3/s，有 2 年超過 3000m3/s，而以后的 47 年間，最大流量沒有一年超過3000m3/s，超過 2021m3/s 有 3 年，超過 1000m3/s 的年份也僅有 11 年。再看看對水庫未建年份資料進行修正后的資料， 54 年中沒有一年超過 3000m3/s，最大流量為 1955 年的 2694m3/s，該流量的經(jīng)驗頻率為 55 年一遇，本次計算的 50 年一 18 遇洪峰流量值為 3476m3/s，是實測系列最大值的 倍，計算成果較符合本流域的實際情況。 綜合以上分析，本次南村站現(xiàn)狀水利工程條件下的設(shè)計洪水成果符合實際情況，合理可靠。 大沽河干流沿河各斷面洪水過程的推求 控制斷面及單元區(qū)劃分 本次分析計算，考慮氣候條件、下墊面條件的影響及流域干支流情 況，全河干流（產(chǎn)芝水庫以下）選取了 8 個控制斷面，分別為洙河口、小沽河口、五沽河口、落藥河口、南村站、流浩河口、南膠萊河口 和 入?？诎藗€控制斷面。 大沽河流域內(nèi)有大型水庫 2 座，中型水庫 8 座。本次考慮氣候條件，下墊面條件的影響及流域干支流情況，根據(jù)產(chǎn)匯流特性的差異，將全流域劃分了產(chǎn)芝、高格莊、北墅、黃同、尹府、洙河區(qū)間、小沽河區(qū)間、五沽河、落藥河、宋化泉、流浩河區(qū)間、王吳、南膠萊河區(qū)間、挪城、桃源河區(qū)間 15 個單元。 各斷面洪水過程的推求 各單元區(qū)平均雨量及全流域平均雨量的計算 根據(jù)大沽河雨量站的分布、資料系列情況 ，綜合考慮資料代表性、一致性和可靠性等，確定本次計算采用流域內(nèi)及流域邊界共 33 個雨量站進行流域平均雨量的計算。 由對歷史發(fā)生的較大洪水的分析，本著典型暴雨宜選用量大、強度大且對工程不利的原則，本次選用“ ”暴雨作為典型暴雨。根據(jù)選用的降雨量代表站，采用算術(shù)平均法計算全流域及各單元的流域平均雨量。 19 各控制斷面洪水過程的推求 分單元按暴雨法進行流域產(chǎn)流、匯流、河道匯流演算，各單元洪水 疊加 后求得各控制斷面洪水過程。本次采用典型暴雨法進行暴雨的時空分布分配，前期影響雨量采用本流域的設(shè)計值 45mm。 （ 1） 各單元產(chǎn)匯 流計算 流域內(nèi)產(chǎn)芝水庫有 1959~2021 年的水文觀測資料、尹府水庫有 1960~2021年的水文觀測資料、嵐西頭水文站有 1976~2021 年的水文觀測資料、南村水文站有 1951~2021 年的水文觀測資料，四站均有預(yù)報任務(wù)，有一套完整的產(chǎn)、匯流預(yù)報方案。另外王吳水庫有 1960~1985 年的水文觀測資料、吳家口水文站有1951~1970 年的水文觀測資料，王吳水庫當時預(yù)報任務(wù)，曾分析一套完整的產(chǎn)、匯流預(yù)報方案。本次根據(jù)吳家口水文站的水文觀測資料，分析了吳家口水文站以上王吳水庫以下區(qū)間流域的產(chǎn)流方案和匯流方案。 出 于對工程的安全考慮，采用以上各站降雨徑流點據(jù)的外包線作為本次各單元的產(chǎn)流方案。 由以上各站的暴雨洪水資料，分析了不同地區(qū)不同下墊面條件的多條時段單位線，最后綜合出適用于流域中不同單元的經(jīng)驗單位線。 根據(jù)各單元計算的流域平均雨量及確定的各單元產(chǎn)匯流方案推求出各單元的洪水過程。 （ 2） 水庫調(diào)洪演算 將產(chǎn)芝、尹府、高格莊、北墅、黃同、宋化泉、王吳、挪城八座大中型水庫推求的入庫洪水過程按照汛期控制運用方案進行調(diào)洪演算，以汛中限制水位作為起調(diào)水位。 20 （ 3） 河道匯流演算 采用華水匯流系數(shù)法，根據(jù)南村水文站洪水預(yù)報方案中的分析成果， 本次經(jīng)綜合分析了各河段的特性，確定出各河段的匯流系數(shù)。各單元推求的洪水過程經(jīng)河道匯流演算到各控制斷面。 （ 4） 各控制斷面洪水過程的合成 每一控制斷面的洪水過程，由該斷面以上的各單元洪水過程，經(jīng)河道匯流演算與控制斷面區(qū)間的洪水 疊加 ?；髁堪疵恳话倨椒焦?1m3/s 的流量計算。 設(shè)計洪水成果 根據(jù)以上計算方法，本次各支流匯入點流量成果如 表 24。 表 24 大沽河各支流匯入點 斷面洪峰流量 （ m3/s） 項目 產(chǎn)芝 水庫 洙河 河口 小沽 河口 五沽 河口 落藥 河口 南村站 流浩 河口 南膠萊 河口 入?？? 控制面積 （ km2） 879 1436 2605 3270 3661 3724 4047 5557 5年一遇 479 717 1178 1119 1165 1153 1198 1289 1279 20年一遇 479 1045 1906 2426 2530 2508 2493 2907 2715 50年一遇 479 1358 2725 3295 3474 3476 3571 4040 4000 起始水位的確定 在大沽河入海口段，先于河西屯設(shè)潮位站（ 1954～ 1966 年， 13 年），后遷站至山角底村處（ 1967～ 1986 年， 20 年），以測定大沽河入膠州灣口處的潮汐情況。由于地理、氣象、水文諸方面的綜合作用，該區(qū)域每年的高潮多發(fā)生在七、八月份，常有天文高潮位與河道的洪峰相遭遇，致使實測潮位成為潮汐、河口洪水共同作用的復(fù)合潮位。 21 據(jù) 1972 年山東省水文總站分析的兩站相關(guān)關(guān)系，將山角底站的潮位資料延長，形成從 1954～ 1985 年計 32 年的系列，并對之作 概率 分析，得 到在不同重現(xiàn)期，大沽河口段由潮汐、洪水共同形成的最高水位值。本規(guī)劃用之以作為推算大沽河各斷面水位的起始斷面水位。山角底至南莊的河道 2021 年進行了治理，推算時此段河道的斷面數(shù)據(jù)采用治理后的實際數(shù)據(jù)進行推算。 表 25 不同重現(xiàn)期之起始斷面水位 重現(xiàn)期（年） 平均 5年 10年 20年 30年 50年 水位（ m） 注：黃海高程系。 在推求水位時，大沽河的 洪水與起始斷面水位（入?？诙螐?fù)合潮位）采用同一重現(xiàn)期，使復(fù)合潮位與本規(guī)劃的設(shè)計洪水有著內(nèi)在的統(tǒng)一聯(lián)系。 大沽河行洪能力分析 大沽河通過近幾年的治理，河道行洪能力 得到 提高。 按天然河道水面線計算公式 推求大沽河 規(guī)劃 各 斷面 50 年一遇洪水位 見 表 26。 22 表 26 大沽河 規(guī)劃 斷面要素表 樁號 規(guī)劃河 底高程（ m） 主河槽 寬度 （ m） 堤距 （ m） 灘地高程 （ m） 堤頂高程 （ m） 洪峰流量（ m3/s） 洪水位 （ m） 50年超高 （ m） 注釋 左 右 左 右 50年一遇 50年一遇 左 右 0+000 479 同三高速公路橋 0+500 479 1+000 479 1+070 479 萊西市城西橋 1+500 479 1+770 479 濰萊高速公路橋 2+000 479 2+500 479 3+000 479 3+500 479 3+930 479 萊西城橋 4+000 479 4+500 479 5+000 840 479 5+300 920 479 上海路橋 5+800 730 479 6+500 980 479 7+000 1000 479 7+600 1100 479 23 樁號 規(guī)劃河 底高程（ m） 主河槽 寬度 （ m） 堤距 （ m） 灘地高程 （ m） 堤頂高程 （ m） 洪峰流量（ m3/s） 洪水位 （ m） 50年超高 （ m） 注釋 左 右 左 右 50年一遇 50年一遇 左 右 8+000 1050 479 8+575 667 479 8+780 716 479 橋 9+300 635 479 9+800 535 479 橋 10+000 535 479 10+500 680 479 11+000 700 479 11+500 700 479 12+000 515 479 12+500 900 479 13+000 1100 1358 潴河匯入 14+000 945 1358 14+500 1358 15+000 1358 15+500 1358 15+645 1060 1358 橋 16+000 1145 1358 16+500 1060 1358 17+000 1020 1358 17+500 1112 1358 24 樁號 規(guī)劃河 底高程（ m） 主河槽 寬度 （ m） 堤距 （ m） 灘地高程 （ m） 堤頂高程 （ m） 洪