黃河河口不同流路入海泥沙對(duì)下游影響二維數(shù)模分析

李東風(fēng) 張紅武 鐘德鈺 胡建永 馬梁超 江曹棟

摘?要：黃河河口入海泥沙在漲潮流和攔門沙的阻滯作用下對(duì)黃河下游泥沙運(yùn)動(dòng)產(chǎn)生重大影響。以已建立的河口海洋水沙二維數(shù)學(xué)模型為基礎(chǔ)，對(duì)規(guī)劃運(yùn)用初期的清水溝北汊流路、刁口河流路和十八戶流路典型水沙過程進(jìn)行了計(jì)算，分析了3條流路入海口的海洋動(dòng)力和泥沙運(yùn)動(dòng)，得到了入海泥沙對(duì)黃河下游的影響范圍為距濼口斷面90～215 km。研究了未來50 a進(jìn)入下游的水沙設(shè)計(jì)和流路規(guī)劃成果，采用清水溝、刁口河聯(lián)合運(yùn)用方案，計(jì)算了未來50 a入海水沙運(yùn)動(dòng)對(duì)下游的影響，結(jié)果表明：在此方案下，西河口水位維持在10.23～11.18 m。結(jié)合黃河三角洲高效生態(tài)經(jīng)濟(jì)區(qū)發(fā)展需求，清水溝和刁口河聯(lián)合運(yùn)用方案為最優(yōu)方案，十八戶流路海域封閉且海洋動(dòng)力比較弱，可以作為遠(yuǎn)景備用流路。

關(guān)鍵詞：黃河河口;二維數(shù)學(xué)模型;入海泥沙;黃河下游;清水溝流路;刁口河流路;十八戶流路

中圖分類號(hào)：TV142;TV882.1?文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2021.05.004

Abstract： The sediment flowing into the Yellow River estuary has a great influence to the sediment movement in the lower Yellow River under the blocking action of the tidal current and sandbars. Based on the established two-dimensional mathematical model of estuarine water and sediment， the typical water and sediment processes of North Branch estuarine of Qingshuigou， Diaokou River estuarine and Shibahu River estuarine in the early stage of planning and operation were calculated and the marine dynamics and sediment movement of the three channels into the estuary were analyzed. The results show that the influence range of sediment to the lower Yellow River is 90-215 km away from Luokou section. The results of water and sediment design and flow path planning in the next 50 years were studied. The combined operation scheme 4 of Qingshuigou and Diaokou rivers recommended by China Institute of Water Resources and Hydropower and Yellow River Engineering Consulting Co.， Ltd. was adopted to calculate the impact of water and sediment movement on the downstream in the next 50 years. The results show that under this scheme， the water level of Xihekou will be maintained between 10.23 m and 11.18 m. Combined with the development needs of the efficient ecological economic zone in the Yellow River delta， scheme 4 of the joint operation of Qingshuigou and Diaokou river is the optimal flow path and the sea area of Shibahu flow path is closed and the marine power is relatively weak， so it can be used as the future backup flow path.

Key words： Yellow River estuary; 2D numerical simulation; sediment into sea; Lower Yellow River; Qingshuigou flow path; Diaokouhe flow path; Shibahu flow path

1?引?言

黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展已上升為重大國家戰(zhàn)略[1]。有計(jì)劃地穩(wěn)定使用河口入海流路是黃河三角洲高效生態(tài)經(jīng)濟(jì)區(qū)發(fā)展的基礎(chǔ)。黃河河口治理與黃河下游治理密切相關(guān)，黃河口漲潮和攔門沙阻滯黃河入海水沙運(yùn)動(dòng)，對(duì)黃河下游防洪將產(chǎn)生不利影響。黃河河口入海流路穩(wěn)定與黃河下游治理特別是山東黃河密切相關(guān)[2]，2013年3月國務(wù)院批復(fù)的《黃河流域綜合規(guī)劃》在涉及“河口入海流路規(guī)劃”時(shí)指出：“入海流路規(guī)劃要遵循黃河河口自然演變規(guī)律，以保障黃河下游防洪安全為前提，以河口生態(tài)良性維持為基礎(chǔ)，充分發(fā)揮三角洲地區(qū)的資源優(yōu)勢，促進(jìn)地區(qū)經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展。綜合考慮各種因素，規(guī)劃期內(nèi)仍主要利用清水溝流路行河，保持流路相對(duì)穩(wěn)定，清水溝流路使用結(jié)束后，優(yōu)先啟用刁口河備用流路;馬新河和十八戶作為遠(yuǎn)景可能的備用流路。”同時(shí)，《黃河流域綜合規(guī)劃》在涉及“黃河長治久安的重大戰(zhàn)略措施”時(shí)還特別強(qiáng)調(diào)：“為了盡量減少河口淤積延伸對(duì)下游河道的反饋影響，要繼續(xù)進(jìn)行現(xiàn)行流路尾閭河道的治理，科學(xué)使用清水溝、刁口河、馬新河等入海流路，合理使用海域容沙空間;繼續(xù)進(jìn)行攔門沙治理，充分發(fā)揮海洋動(dòng)力輸沙能力。”

對(duì)于黃河河口流路穩(wěn)定和入海泥沙對(duì)黃河下游的影響研究，很多學(xué)者從分析實(shí)測資料歸納推理回歸得出經(jīng)驗(yàn)公式等方面進(jìn)行了深入研究，推動(dòng)了黃河河口治理研究。對(duì)于黃河河口治理研究手段的研究，2002年李國英提出建立“數(shù)字黃河”[3]和建立黃河數(shù)學(xué)模擬系統(tǒng)[4]，“黃河數(shù)學(xué)模擬系統(tǒng)建設(shè)要圍繞黃土高原土壤侵蝕模型、水庫調(diào)度模型、黃河下游河道水沙演進(jìn)模型、河口模型、水質(zhì)預(yù)警預(yù)報(bào)模型、寧蒙河段冰凌預(yù)報(bào)模型六大模型系統(tǒng)開展。要充分認(rèn)識(shí)黃河數(shù)學(xué)模擬系統(tǒng)建設(shè)的長期性和復(fù)雜性，持之以恒地進(jìn)行黃河數(shù)學(xué)模擬系統(tǒng)建設(shè)”。其中在論述河口模型時(shí)指出“黃河口不同于一般清水河流的河口，受上游來水來沙及海洋潮汐對(duì)河口泥沙的頂沖影響，其水沙運(yùn)動(dòng)規(guī)律較為復(fù)雜。應(yīng)加強(qiáng)機(jī)理研究，建立基于黃河口演變規(guī)律的數(shù)學(xué)模型，指導(dǎo)黃河口流路的規(guī)劃和治理”。

李澤剛[5]分析了黃河入海泥沙擴(kuò)散、黃河口攔門沙的形成與演變等規(guī)律，提出了建設(shè)西河口水沙控制工程、進(jìn)行河口段河道整治和攔門沙治理等穩(wěn)定黃河河口流路的對(duì)策。王愷忱[6]分析了黃河下游及河口河段沖淤特性，進(jìn)行了黃河河口延伸改道對(duì)下游河道的影響研究，用實(shí)測資料證實(shí)了利津至改道點(diǎn)河段短時(shí)段的沖淤主要受制于河口基準(zhǔn)面的狀況，同時(shí)進(jìn)一步證實(shí)了河口尾閭?cè)牒Ａ髀返慕^對(duì)長度是河口河段和下游河道沖淤幅度與發(fā)展趨勢的制約因素。胡春宏等[7-8]分析了黃河口水沙運(yùn)動(dòng)與演變基本規(guī)律，提出了黃河口治理的方向與措施。張紅武[9]2003年提出構(gòu)建二維河口數(shù)學(xué)模型及河口治理意見。馬睿等[10]基于1958年汛后—2000年汛后利津以下實(shí)測資料建立了河口段各部分累計(jì)淤積量與利津站累計(jì)來沙量的相關(guān)關(guān)系，利用準(zhǔn)二維水沙數(shù)學(xué)模型對(duì)黃河下游及河口段進(jìn)行耦合模擬，分析了黃河小浪底—河口段現(xiàn)狀治理模式和“兩道防線”治理模式對(duì)黃河河口各部分累計(jì)沖淤量、西河口以下河長延伸值等特征量的影響，以及兩種模式下黃河河口段未來50 a治理的效果。余欣等[11]界定了出汊、改道、流路穩(wěn)定及黃河河口穩(wěn)定的內(nèi)涵，系統(tǒng)分析了近期黃河河口水沙及流路、海岸的新變化，通過試驗(yàn)研究了尾閭河道出汊機(jī)制，提出出汊閾值、觸發(fā)條件以及河口流路穩(wěn)定綜合判別指標(biāo)，在此基礎(chǔ)上預(yù)判了清水溝清8汊河的穩(wěn)定性，提出近期流路安排應(yīng)優(yōu)先考慮北汊河。陳雄波等[12]對(duì)清水溝、刁口河流路聯(lián)合運(yùn)用方案進(jìn)行了計(jì)算比選分析。王崇浩等[13-15]利用一二維連接數(shù)學(xué)模型研究了黃河口潮流與泥沙輸移過程，對(duì)河口演變對(duì)下游河道反饋的影響進(jìn)行了分析，對(duì)黃河口清水溝流路沖淤發(fā)展與使用年限進(jìn)行了預(yù)測。王萬戰(zhàn)等[16]分析了黃河口海岸演變規(guī)律。王開榮等[17]對(duì)黃河現(xiàn)行清水溝流路汊河運(yùn)用方案進(jìn)行了探討。趙連軍[18]分析認(rèn)為，黃河下游河道演變與河口演變影響因素繁多，未來必須借助河工模型試驗(yàn)或數(shù)學(xué)模型計(jì)算等手段開展研究才可能徹底揭示黃河口演變對(duì)河道反饋影響的規(guī)律。李東風(fēng)等[19-23]利用差分方法計(jì)算了清水溝和刁口河流路入海泥沙運(yùn)動(dòng)和對(duì)港口的影響，采用有限元方法建立了河口二維數(shù)學(xué)模型并應(yīng)用于河口治理研究，研究了外海邊界條件的影響。上述成果是在老一輩黃河河口治理專家研究成果的基礎(chǔ)上取得的，為黃河河口治理研究打下了堅(jiān)實(shí)的基礎(chǔ)。

隨著高性能計(jì)算機(jī)的應(yīng)用，建立黃河下游長河段河道、海洋二維模型和應(yīng)用模型預(yù)測成為現(xiàn)實(shí)，但建立長河段濼口以下河道至海洋的二維水沙數(shù)學(xué)模型，分析入海泥沙對(duì)黃河下游影響的研究未見報(bào)道。為進(jìn)一步完善黃河河口二維水沙數(shù)學(xué)模型，發(fā)揮二維模型在黃河下游和河口治理研究中的優(yōu)勢，2016年國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2016YFC0402500）在第二課題中列入專題（2016YFC0402502），即未來50 a河口入海泥沙淤積延伸對(duì)黃河下游的影響，主要利用黃河下游和河口二維水沙數(shù)學(xué)模型研究未來50 a黃河入海泥沙對(duì)黃河下游的影響，本文為其部分成果。

2?基本理論

模型基本方程見文獻(xiàn)[19-22]。

（1）初始條件。計(jì)算域內(nèi)，在t=0時(shí)刻，給定水位、潮位、流速和含沙量為0。

（2）邊界條件。進(jìn)口邊界給定濼口斷面水流流量和含沙量。出口邊界用潮位變化控制計(jì)算域與海域相連接的水邊界，根據(jù)物理海洋學(xué)理論，潮高表達(dá)式為

開邊界兩個(gè)端點(diǎn)值由實(shí)測資料給定，其端點(diǎn)內(nèi)部的值由插值方法求得。在閉邊界上，認(rèn)為其法向流速為0，而沿切線方向的流速非0。泥沙運(yùn)動(dòng)邊界條件除了給出入口斷面各節(jié)點(diǎn)的含沙量外，在閉邊界還應(yīng)滿足法向輸沙通量為0。

（3）模型驗(yàn)證。模型驗(yàn)證見文獻(xiàn)[19-22]。

3?計(jì)算條件分析

根據(jù)水利部2010年審查通過的《黃河河口綜合治理規(guī)劃》，入海流路除了清水溝外，還包括刁口河備用流路、遠(yuǎn)景可能的十八戶與馬新河流路。規(guī)劃確定的改道標(biāo)準(zhǔn)為在設(shè)防流量10 000 m3/s時(shí)，西河口（二）站水位不超過12 m（大沽高程），未來規(guī)劃的流路是清水溝流路，清水溝流路行河完成后優(yōu)先啟用刁口河備用入海流路，馬新河和十八戶作為遠(yuǎn)景可能的備用流路;清水溝流路繼續(xù)使用清8汊河，滿足改道條件后改走清水溝北汊，再次滿足改道條件后改走1996年前清水溝流路原河道，使清水溝流路在50 a左右的時(shí)間內(nèi)保持相對(duì)穩(wěn)定。

3.1?計(jì)算范圍和地形

根據(jù)文獻(xiàn)資料分析得到的黃河河口入海泥沙河道淤積延伸對(duì)黃河下游的影響范圍，以及渤海灣電子海圖和黃河河道地形等資料，建立模型的范圍為：黃河下游河道濼口斷面為模型進(jìn)口，渤海灣足夠大的海域?yàn)槟Ｐ偷某隹?。鑒于馬新河流路入海泥沙進(jìn)入渤海灣后海域容沙體積和流路長度等方面的因素，僅布置了清水溝、刁口河和十八戶3條流路，建立黃河下游河道與海洋平面二維數(shù)學(xué)模型。利用2017年汛后大斷面數(shù)據(jù)對(duì)地形進(jìn)行校正，作為計(jì)算的初始地形。模型范圍和入海流路見圖1，各主要斷面距濼口斷面的距離見表1。

3.2?計(jì)算分析

3.2.1?水沙系列設(shè)計(jì)分析

關(guān)于水沙系列研究，陳雄波等[12]依據(jù)2012年通過水利部審查的《黃河水沙調(diào)控體系建設(shè)規(guī)劃》，對(duì)6個(gè)入海流路方案采用3個(gè)50 a水沙序列進(jìn)行計(jì)算，其中2020年后水沙序列均采用無古賢水庫情況，得到的利津站水沙量見表2[12]。

王崇浩等[15]參考小浪底水庫的設(shè)計(jì)水沙序列，并考慮有、無古賢水庫的調(diào)水調(diào)沙過程，經(jīng)下游河道沖淤計(jì)算得到艾山斷面未來80 a（2000年起算）水沙過程，見圖2。

綜合分析未來50 a水沙序列，本次計(jì)算使用文獻(xiàn)[15]考慮古賢水庫的成果，再根據(jù)《黃河水沙調(diào)控體系建設(shè)規(guī)劃》，按照不同時(shí)段和流量級(jí)（小于800 m3/s，大于800 m3/s、小于2 600 m3/s，大于2 600 m3/s、小于4 000 m3/s和大于4 000 m3/s）的天數(shù)得到濼口斷面的水沙過程。

3.2.2?入海流路計(jì)算方案的選擇

文獻(xiàn)[12]設(shè)計(jì)計(jì)算了6個(gè)入海流路方案，本次計(jì)算使用其推薦的方案4，即在清水溝按汊河、北汊1、北汊2和清水溝原河道行河的基礎(chǔ)上，將刁口河作為分洪通道方案。2030年后刁口河分洪通道平時(shí)關(guān)閉，利津流量在4 000 m3/s以上時(shí)刁口河開始過流;利津流量為4 000～7000 m3/s時(shí)維持清水溝流路4 000 m3/s的過流量，刁口河過流量為0～3 000 m3/s;利津流量大于7 000 m3/s以后，刁口河仍維持3 000 m3/s的過流量，其余通過清水溝下泄。

以上述條件為計(jì)算的初始和邊界條件，計(jì)算入海泥沙在海洋動(dòng)力作用下的輸運(yùn)擴(kuò)散和淤積沖刷情況。選取濼口以下主河槽內(nèi)和入海徑流方向上的655個(gè)點(diǎn)分析沿程水位和河床高程的變化，選擇重要斷面（見圖1）分析水位變化，水位分析點(diǎn)選擇在主槽位置。

3.3?流路入海泥沙淤積延伸影響分析

按照2018年水沙條件分析河流和河口入海泥沙運(yùn)動(dòng)情況。

3.3.1?北汊流路入?？诤Ｑ髣?dòng)力分析

圖3為濼口以下河段和黃河河口北汊入海后流速等值線與流速矢量疊加圖?？梢钥闯觯核鲝倪M(jìn)口流入河道后，在彎曲河道內(nèi)流動(dòng)，在寬河段也只在主槽內(nèi)流動(dòng)，包括孤東油田的北大堤和滯洪區(qū)附近的寬河段。水流流入河口海域后，泥沙在徑流作用下運(yùn)動(dòng)，一方面向深海方向運(yùn)動(dòng)，另一方面在海洋動(dòng)力作用下呈現(xiàn)出沿岸流的特性;黃河三角洲沿岸海域在三角洲頂點(diǎn)存在的高流速區(qū)域（垂線流速的平均值在0.45 m/s以上）和清水溝1996年沙嘴處的高流速區(qū)域都與實(shí)測資料一致。圖4為河道水位和海洋潮位等值線分布圖，可以看出，河道內(nèi)水位比降小，海域內(nèi)比降大，鑒于海域內(nèi)水位變幅比較大，海域潮位變幅只顯示了0 m等值線。河道水位變化受模型進(jìn)口濼口水沙過程和潮汐升降等海洋動(dòng)力條件的影響，同時(shí)攔門沙規(guī)模不斷擴(kuò)大和河口向外海延伸也對(duì)其造成影響。

3.3.2?泥沙運(yùn)動(dòng)含沙量分布分析

選擇濼口流量為500 m3/s、含沙量為8.0 kg/m3的典型水沙過程對(duì)泥沙在河道和河口中的運(yùn)動(dòng)進(jìn)行分析。圖5（a）為河道水流含沙量沿程分布，水流從進(jìn)口流入河道后在彎曲河道主槽內(nèi)流動(dòng)，流速比較大，含沙量也比較大，在比較窄的河道流速增大引起河床沖刷，水流含沙量增加到20 kg/m3。圖5（b）和（c）為剛開始落潮和漲潮時(shí)水流含沙量在口門的平面分布情況。水流流入河口海域后，泥沙一方面在徑流作用下向深海方向運(yùn)動(dòng)，由于流量比較小，漲潮動(dòng)力比較大，因此高渾濁度區(qū)域在口門附近;另一方面在海洋動(dòng)力作用下泥沙沿海岸運(yùn)動(dòng)，口門兩側(cè)水流含沙量較大，與沿岸流的特性一致。受北汊口門沙嘴的影響，落潮時(shí)期泥沙沙嘴偏向西北方向，漲潮時(shí)偏向東南方向，同時(shí)受漲落潮的影響，在口門附近含沙量最大，泥沙聚集在口門，而且含沙量等值線呈橢圓狀，中心含沙量最大。由圖5（c）可知，落潮結(jié)束、剛開始漲潮時(shí)口門中心最大含沙量為20 kg/m3;由圖5（b）可知，漲潮結(jié)束、剛開始落潮時(shí)挾沙水流動(dòng)能最小，含沙量變化有一個(gè)滯后過程，此時(shí)對(duì)應(yīng)含沙量最大，最大含沙量為32 kg/m3。泥沙在口門聚集是形成攔門沙的重要原因。另外，可見沿海岸和河道出現(xiàn)44 kg/m3的高含沙量區(qū)域，主要是在河道、海岸水動(dòng)力和邊界條件的影響下形成的。

3.3.3?淤積沖刷厚度分布分析

入海泥沙3.18億t/a條件下河道和河口泥沙淤積分布見圖6。可以看出，泥沙入海時(shí)在潮汐和攔門沙的作用下主要淤積在口門附近，淤積最大厚度為5.23 m，主要淤積在入海水流方向上6～8 km的范圍內(nèi)，而下游河道表現(xiàn)微沖和微淤的狀態(tài)，口門河道沖刷比較明顯，窄河段河床沖刷深度較大，沖刷淤積影響范圍為距離濼口斷面150～170 km處，即道旭斷面和利津斷面之間。主要原因是河道泥沙沖刷淤積的結(jié)果滯后于海洋動(dòng)力和攔門沙的影響。

3.3.4?水位和水位差分析

2018年水沙條件下濼口3 000 m3/s流量時(shí)汛前和汛后其以下縱剖面的水位和水位差見圖7?？梢钥闯觯涸诔毕蛿r門沙的作用下，下游河道和入?？陂T附近泥沙沖淤引起了水位變化，在3 000 m3/s流量時(shí)最大水位差為0.19 m;距河口的最遠(yuǎn)影響點(diǎn)為濼口以下約90 km處，即董家斷面和清河鎮(zhèn)斷面之間。主要原因是其他條件相同的情況下，3 000 m3/s大流量時(shí)，河床沖刷比較嚴(yán)重，沖刷后水位下降比較大，水位差也比較大，影響范圍距河口也比較遠(yuǎn)。

3.4?刁口河入海泥沙淤積對(duì)下游河道的影響

在進(jìn)口邊界條件即濼口流量為500 m3/s、水流含沙量為35 kg/m3下，得到刁口河流路1.5億t/a泥沙流入海域后淤積沖刷計(jì)算結(jié)果。選擇縱剖面水位分析河道水位的變化（選擇一日中同一時(shí)刻潮位條件下的水位）。從圖8可以看出河床淤積與沖刷的分界點(diǎn)為距濼口斷面185.33 km處，即利津斷面附近。該分界點(diǎn)以下河道淤積，海域淤積厚度較大。

3.5?十八戶流路入海泥沙運(yùn)動(dòng)及對(duì)下游水位的影響

十八戶流路入海位置特殊，作為遠(yuǎn)期備用流路需要進(jìn)行研究和探索，這里分析濼口流量為3 000 m3/s和含沙量為35 kg/m3條件下入海泥沙運(yùn)動(dòng)及對(duì)下游的影響。

3.5.1?入海泥沙運(yùn)動(dòng)分析

從十八戶入海泥沙運(yùn)動(dòng)、漲潮和落潮含沙量分布、一日漲潮落潮對(duì)河口和下游水位的影響圖可以看出，清水溝流路經(jīng)過30多a的使用，沙嘴突出萊州灣，使十八戶入海河口口門處于南北西三面陸地環(huán)繞，僅一面與萊州灣連通，使得入海泥沙運(yùn)動(dòng)擴(kuò)散的范圍減小。再由圖9可以看出，由于潮流等海洋動(dòng)力較弱，因此向外海輸送泥沙的能力減弱，造成入海泥沙淤積，而且淤積在口門附近。

3.5.2?漲潮流對(duì)下游水位的影響

由于入海口門附近海域海洋動(dòng)力較弱，因此漲潮流對(duì)河口入海徑流影響較小，明顯影響距離為入?？谝陨霞s15 km。

3.5.3?十八戶改道對(duì)改道點(diǎn)上下游水位的影響

十八戶流路改道位于距離濼口斷面約205.2 km位置。從圖10可以看出，改道點(diǎn)以下水位明顯下降，改道點(diǎn)以上水面比降遠(yuǎn)小于改道點(diǎn)以下水面比降。主要原因是入海水流流路在十八戶改道點(diǎn)附近的彎曲比較大，改道點(diǎn)以上水流在彎道附近的阻力大，水位下降緩慢，改道點(diǎn)距入?？诰嚯x短，水沙運(yùn)動(dòng)阻力小，水面比降大。

3.6?未來50 a入海泥沙對(duì)黃河下游的影響

根據(jù)3.2節(jié)設(shè)計(jì)的水沙過程和計(jì)算方案，考慮了各流路海域的狀況等，對(duì)各流路的使用年限進(jìn)行了分配，清水溝流路按照北汊河使用14 a、北汊1使用12 a、北汊2使用12 a、原河道使用12 a，刁口河作為分洪通道方案，計(jì)算未來50 a流路河流和河口入海泥沙淤積對(duì)河道的影響。使用Dell工作站并行計(jì)算方法，計(jì)算1 d相當(dāng)于模擬計(jì)算河流河口入海水沙運(yùn)動(dòng)1 a。

（1）河口入海泥沙淤積分析。圖11為河口泥沙淤積厚度平面分布圖，可以看出：一方面，入海泥沙主要淤積在口門附近，最大淤積厚度等值線呈橢圓狀，橢圓長軸與海岸方向一致，說明泥沙在沿岸流的作用下沿海岸方向往復(fù)運(yùn)動(dòng);另一方面，從泥沙淤積厚度0 m等值線的范圍看，入海泥沙向東南方向輸運(yùn)擴(kuò)散的距離和范圍大于向北部擴(kuò)散的范圍。圖12和圖13分別為河口口門方向上剖面淤積厚度和河床高程，可以看出：隨著入海泥沙的增加，最大淤積厚度不斷增大，流路使用初期其最大值增加比較快，向外海淤積的速度也比較快;隨著時(shí)間的增長，最大淤積厚度增加速度和向外海淤積的速度減小，泥沙主要向口門兩側(cè)輸運(yùn)和淤積。

（2）流路改道對(duì)下游的影響范圍。河口流路改道對(duì)下游的影響是一個(gè)復(fù)雜的過程，其主要表現(xiàn)為：在相同的水沙條件和河床比降條件下，流路改道初期，流路流程短、比降大，同流量下水流流速和挾沙能力增加，河道發(fā)生溯源沖刷;隨著流路使用時(shí)間的增長，溯源沖刷向上游發(fā)展的同時(shí)，河口入海流路向外海延伸，河道長度增長，溯源沖刷影響逐漸減弱直到結(jié)束。根據(jù)計(jì)算分析，入海泥沙河口淤積延伸對(duì)下游的影響范圍為距濼口斷面90～215 km，流路使用初期對(duì)下游的影響距離口門最遠(yuǎn)，隨著流路使用年限的增長對(duì)下游的影響越來越小，未來50 a入海泥沙對(duì)下游的影響范圍在董家斷面到利津斷面之間。其中：清水溝北汊河河口使用對(duì)下游的影響在董家斷面附近，北汊1河口使用對(duì)下游的影響在前左斷面附近，北汊2河口使用對(duì)下游的影響在十八戶斷面附近，原清水溝河道河口使用對(duì)下游的影響在利津斷面附近。

（3）對(duì)西河口水位的影響。根據(jù)入海流路改道標(biāo)準(zhǔn)，當(dāng)西河口流量為10 000 m3/s時(shí)水位達(dá)到12 m，因計(jì)算使用的設(shè)計(jì)水沙系列并沒有10 000 m3/s流量的，故確定西河口水位時(shí)，先計(jì)算出未來50 a不同流量的利津和西河口水位，再根據(jù)文獻(xiàn)[24-25]中的西河口和利津水位流量關(guān)系曲線，應(yīng)用水位平行抬高法插值推算分析，最終得到西河口10 000 m3/s流量對(duì)應(yīng)水位在流路使用初期和末期在10.23～11.18 m之間。其中：清水溝清8汊河河口使用14 a西河口水位在10.23～11.09 m之間，北汊1河口使用12 a西河口水位在10.39～11.18 m之間，北汊2河口使用12 a西河口水位在10.51～11.13 m之間，原清水溝河道河口使用12 a西河口水位在10.65～11.15 m之間。

4?結(jié)?論

以已建立和驗(yàn)證的黃河河口海洋水沙二維數(shù)學(xué)模型為基礎(chǔ)，計(jì)算分析了清水溝北汊流路、刁口河流路和十八戶流路入海泥沙運(yùn)動(dòng)對(duì)黃河下游的影響范圍，對(duì)清水溝流路按清8汊河、北汊1和北汊2、原河道行河，刁口河作為分洪通道方案計(jì)算分析了未來50 a河道和河口入海水沙運(yùn)動(dòng)對(duì)下游的影響，結(jié)果表明該方案清水溝流路可以使用50 a以上。3條入海流路中，十八戶流路入海口三面被陸地環(huán)繞再加上海洋動(dòng)力較弱和十八戶改道點(diǎn)附近河道彎曲水流阻力大，僅可以作為遠(yuǎn)景備用流路;刁口河流路入海口附近的海洋動(dòng)力強(qiáng)度和容沙能力處于清水溝流路和十八戶流路之間，宜作為清水溝流路行河結(jié)束后優(yōu)先選擇的流路;清水溝流路海洋動(dòng)力強(qiáng)，泥沙海域容沙能力大，結(jié)合工程措施可以將入海泥沙輸送到深海。

參考文獻(xiàn)：

[1]?習(xí)近平.在黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展座談會(huì)上的講話[J].中國水利，2019（20）：1-3.

[2]?張紅武，張羅號(hào)，景喚，等.山東對(duì)黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展的作用不可替代[J].水利水電技術(shù)（中英文），2021，52（1）：1-21.

[3]?李國英.建設(shè)“數(shù)字黃河”工程[J].中國水利，2002（2）：29-32，80.

[4]?李國英.建立黃河數(shù)學(xué)模擬系統(tǒng)[J].中國水利，2006（21）：36-38.

[5]?李澤剛.黃河近代河口演變基本規(guī)律與穩(wěn)定流路治理[M].鄭州：黃河水利出版社，2006：164-191.

[6]?王愷忱.黃河河口的演變與治理[M].鄭州：黃河水利出版社，2010：195-259.

[7]?胡春宏，曹文洪.黃河口水沙變異與調(diào)控I：黃河口水沙運(yùn)動(dòng)與演變基本規(guī)律[J].泥沙研究，2003，28（5）：1-8.

[8]?胡春宏，曹文洪.黃河口水沙變異與調(diào)控Ⅱ：黃河口治理方向與措施[J].泥沙研究，2003，28（5）：9-14.

[9]?張紅武.黃河河口數(shù)學(xué)模型及河口治理意見[C]//中國水利學(xué)會(huì)，黃河研究會(huì).黃河河口問題及治理對(duì)策研討會(huì)專家論壇文集.鄭州：黃河水利出版社，2003：124-126.

[10]?馬睿，韓鎧御，張紅武，等.不同治理模式下黃河河口段沖淤演變的數(shù)值模擬[J].水力發(fā)電學(xué)報(bào)，2018，37（4）：68-78.

[11]?余欣，吉祖穩(wěn)，王開榮，等.黃河河口演變與流路穩(wěn)定關(guān)鍵技術(shù)研究[J].人民黃河，2020，42（9）：66-70，129.

[12]?陳雄波，雷鳴，王鵬.清水溝、刁口河流路聯(lián)合運(yùn)用方案比選[J].海洋工程，2014，32（4）：117-123.

[13]?王崇浩，曹文洪，張世奇.黃河口潮流與泥沙輸移過程的數(shù)值研究[J].水利學(xué)報(bào)，2008，39（10）：1256-1263.

[14]?王崇浩，張世奇，曹文洪.黃河口演變對(duì)下游河道反饋影響的計(jì)算研究[J].水利水電技術(shù)，2007，38（9）：16-20，31.

[15]?王崇浩，安催花，曹文洪，等.黃河口清水溝流路沖淤發(fā)展與使用年限預(yù)測[J].水利學(xué)報(bào)，2008，39（6）：709-713.

[16]?王萬戰(zhàn)，張華興.黃河口海岸演變規(guī)律[J].人民黃河，2007，29（2）：27-28，32.

[17]?王開榮，凡姚申，韓沙沙，等.黃河現(xiàn)行清水溝流路汊河運(yùn)用方案探討[J].海洋科學(xué)，2020，44（10）：91-100.

[18]?趙連軍.黃河下游河道演變與河口演變相互作用規(guī)律研究[D].武漢：武漢大學(xué)，2005：111-143.

[19]?李東風(fēng)，李澤剛，張青玉.清水溝北汊流路入海泥沙對(duì)東營港影響的數(shù)值分析[J].黃渤海海洋，1998，16（1）：2-7.

[20]?李東風(fēng)，李澤剛.釣口河故道分洪入海泥沙對(duì)東營港影響的計(jì)算研究[J].人民黃河，1998，20（3）：3-5.

[21]?李東風(fēng)，張紅武，鐘德鈺，等.黃河河口水沙運(yùn)動(dòng)的二維數(shù)學(xué)模型[J].水利學(xué)報(bào)，2004，35（6）：1-6，13.

[22]?李東風(fēng)，張紅武，鐘德鈺，等.黃河河口潮流和泥沙淤積過程數(shù)值分析研究[J].水利學(xué)報(bào)，2004，35（11）：74-80.

[23]?鄒冰，李東風(fēng)，張紅武，等.河口導(dǎo)堤對(duì)水沙運(yùn)動(dòng)影響的二維數(shù)值模型[J].水利學(xué)報(bào)，2006，37（7）：880-885.

[24]?張翠萍.黃河下游利津至西河口河段河道水沙沖淤特性分析[C]//《水動(dòng)力學(xué)研究與進(jìn)展》編委會(huì).第十三屆全國水動(dòng)力學(xué)研討會(huì)文集.北京：海洋出版社，1999：410-416.

[25]?陳雄波，邱衛(wèi)國，劉娟.黃河河口“一水三流”行河方式研究[J].海岸工程，2019，38（3）：194-202.

【責(zé)任編輯?張華巖】