我國河流水流泥沙數(shù)值模擬技術(shù)進展與應(yīng)用

張明進 ，張華慶 ，白玉川

（1.天津大學建筑工程學院，天津300072；2.交通運輸部天津水運工程科學研究所工程泥沙交通行業(yè)重點實驗室，天津300456）

河流水流泥沙數(shù)值模擬始于20世紀60年代，20世紀70年代以后逐步成熟。一維水沙數(shù)學模型發(fā)展最早，目前已建立了相對比較成熟的不同輸沙理論體系的模擬系統(tǒng)［1］；在20世紀80年代末也建立起平面二維水流泥沙數(shù)學模型和針對水流和泥沙側(cè)向變化較小的立面二維模型。20世紀90年代以來，隨著計算機的發(fā)展和水流泥沙研究領(lǐng)域的擴展和研究水平的不斷深入，三維水流泥沙數(shù)學模型在河流局部河段、湖泊泥沙問題研究方面得到了一定的發(fā)展與工程應(yīng)用［1-7］。目前，一、二維泥沙數(shù)學模型已相對比較成熟，三維模型也能應(yīng)用來解決一些具體問題。近幾十年、尤其是近20年來，我國水利水運工作者就內(nèi)河航道與港口泥沙輸移、河床演變與航道整治等進行了大量的研究，在數(shù)值模擬領(lǐng)域取得了許多重要的成果［8-14］。在模擬系統(tǒng)發(fā)展完善、數(shù)值計算方法、軟件系統(tǒng)集成等多方面取得了重要進展，在解決河床演變與整治、航道與港口工程泥沙等實際問題過程中發(fā)揮了重要的作用。

1 模擬系統(tǒng)技術(shù)進展與應(yīng)用

適用于不同目的，適應(yīng)不同水沙條件的模擬系統(tǒng)發(fā)展完善情況如下：

1.1 一維水沙數(shù)值模擬系統(tǒng)

適用于長河段和河網(wǎng)水沙模擬。一維水沙數(shù)值模擬是發(fā)展最早，理論基礎(chǔ)和數(shù)值離散格式都相對比較成熟的模擬系統(tǒng)。近10 a來，一維水沙數(shù)值模擬系統(tǒng)主要是在長河段和大型河網(wǎng)計算格式改進、非恒定水流模擬、泥沙模型改進等方面進行［1，8，13］。其中，計算格式的改進體現(xiàn)在水沙模型求解的汊點分組解法、半隱式有限元方法和桿系數(shù)據(jù)結(jié)構(gòu)等；在泥沙模型方面引入了非均勻沙不平衡輸沙理論，對不同的沖淤模式建立了不同的輸沙理論體系，并在非均勻沙挾沙能力、恢復(fù)飽和系數(shù)、混合層厚度、斷面沖淤面積分配等方面進行了改進［14-20］。

1.2 二維水沙數(shù)值模擬系統(tǒng)

包括平面二維和立面二維兩種模式。平面二維模擬系統(tǒng)適用于水平尺度遠大于垂向尺度的寬廣水域，如河流、湖泊、海岸等，水力參數(shù)（如流速、流向等）在垂向變化小于水平方向變化，其流態(tài)可用沿水深平均流動量來表示。立面二維模擬系統(tǒng)適用于窄深的水流通道，水力參數(shù)及有關(guān)變量（如流速、溫度、含沙量等）的垂向變化大于水平橫向變化，水流及其含有物的變化可用水平橫向平均值沿水深的變化來表示。

近20 a來，平面二維水沙模擬系統(tǒng)是發(fā)展最快、工程應(yīng)用最廣泛的模擬手段［21-35］。開發(fā)了適應(yīng)河道復(fù)雜邊界的正交曲線網(wǎng)格和非結(jié)構(gòu)多重網(wǎng)格模型，引入了高效的數(shù)值離散格式和并行計算技術(shù)，將彎道環(huán)流輸沙理論引入平面二維模擬系統(tǒng)中，建立了全沙輸移泥沙運移模型，引入了非均勻沙不平衡輸沙理論，灰色理論和拓撲技術(shù)得到了應(yīng)用，集成了GIS、計算可視化和過程模擬計算運行環(huán)境，開發(fā)了軟件化產(chǎn)品（如TK-2DC系統(tǒng)）等［36］。

立面二維模擬系統(tǒng)在異重流的入潛條件等方面進行了理論探索［1，34］，開發(fā)了異重流數(shù)學模型。水流泥沙運動方程中直接包含了密度差動力項，在模擬水庫異重流運動、干流向支流倒灌異重流、引航道泥沙淤積、水庫排沙等工程中得到了應(yīng)用和檢驗。

1.3 三維水沙數(shù)值模擬系統(tǒng)

三維水沙數(shù)值模擬系統(tǒng)一般由雷諾平均連續(xù)性方程和納維斯托克斯方程確定的流場及泥沙連續(xù)性方程和河床變形方程組成。主要發(fā)展體現(xiàn)在垂向空間離散坐標系的選擇、泥沙近底邊界條件的處理、自由表面的跟蹤計算、紊流模型等幾個方面［37-42］。

（1）垂向空間離散坐標系。在三維數(shù)值模型中，模型垂向空間的離散大多數(shù)均在笛卡兒坐標系或σ坐標系下進行。由于笛卡兒坐標系統(tǒng)存在適應(yīng)不規(guī)則邊界差和固定分層分辨率低的不足，三維水沙數(shù)值模擬多采用σ坐標變換。

（2）泥沙近底邊界條件。由于床面附近泥沙交換現(xiàn)象極其復(fù)雜，泥沙交換機理至今仍不十分清楚，三維水沙數(shù)值模擬中泥沙床面邊界條件大多為半理論半經(jīng)驗型公式，分為濃度型、梯度型和通量型。根據(jù)近底非均勻沙交換的隨機-力學模式，從理論上建立起了初步的非均勻沙床面平衡濃度公式，將床面平衡濃度表示為起動、不止動和懸浮概率以及單步運動平均時間比B的函數(shù)關(guān)系。殘差修正灰色系統(tǒng)模型和指數(shù)擬合關(guān)系用于懸移質(zhì)平衡濃度豎向分布實測資料獲得床面平衡濃度值，據(jù)此確定理論公式中的參數(shù)B，該計算方法結(jié)果相對較好。

（3）自由表面的計算。現(xiàn)有的計算自由表面的方法主要有剛蓋假定、水位函數(shù)法、MAC法、Lagrangian與ALE法、VOF法。目前已經(jīng)引入物理海洋學中的模式分離法和空氣動力學中的高效離散格式。自由表面直接由垂直積分后的方程組來計算，再將計算后的自由表面賦給三維計算模塊，這種處理方式快速穩(wěn)定。

（4）紊流模型。κ-ε雙方程模型是最簡便的精細紊流模型，也是目前應(yīng)用最廣泛的紊流模型，但其模型假定渦粘性系數(shù)是各向同性的標量，這導(dǎo)致水平和垂向擴散系數(shù)存在較大差異。內(nèi)河水沙模擬中已引入各向異性的紊流模型，在平面上紊動粘性系數(shù)采用大渦模型，垂向上紊動粘性系數(shù)采用κ-ε模型，其能夠較好地反映平面流動特點，又能夠模擬內(nèi)部紊動及與此密切相關(guān)的物理量沿垂線分布規(guī)律。

平面二維水沙數(shù)值模擬和三維水沙數(shù)值模擬在長江、珠江、海河、松花江等河流的航道整治、港口建設(shè)、通航論證等工程項目中得到了廣泛的應(yīng)用，為工程前期研究、初步設(shè)計、后期效果評估等提供了重要的技術(shù)支持。

立面二維水沙數(shù)學模型在水庫調(diào)度對溪洛渡電站下游水溫影響、三峽工程引航道淤積、劉家峽水庫及小浪底水庫異重流排沙等工程問題研究中得到了應(yīng)用。

2 數(shù)值模擬計算方法的一些創(chuàng)新及應(yīng)用

2.1 一維大型河網(wǎng)水沙計算快速求解方法

河網(wǎng)模型通常分為兩類，一類是一維圣維南方程求解法，另一類是為“組合單元法”。針對大型河網(wǎng)計算離散方程幾乎是滿秩的大型矩陣，求解效率低和在方案研究階段需要增加人工運河或某一河段增加斷面時，節(jié)點重新編碼繁瑣，耗時多，易出錯的特點。一些水運工程科研院所開發(fā)了基于半隱式有限元方法［20］、應(yīng)用汊點分組思想求解，模擬河網(wǎng)水流、泥沙動力學模型。

珠江三角洲水系為放射狀汊道河系（圖1），包括西、北江思賢滘以下和東江石龍以下河網(wǎng)水系及入注三角洲的中、小河流。該水系河網(wǎng)縱橫交錯，共有八大入海口門，分水分沙關(guān)系及河床演變十分復(fù)雜［20］。近十幾年來，由于珠江三角洲水系航道疏浚整治、圍墾造陸及大規(guī)模無序采沙等，改變了自然淤積的規(guī)律，河道普遍呈現(xiàn)加深下切的態(tài)勢。這些人類活動對河道行洪能力、納潮能力和輸沙能力等產(chǎn)生了較大影響。為了配合珠江河口規(guī)劃治理研究，研究西江、北江分水、分沙規(guī)律，以及網(wǎng)河區(qū)防洪規(guī)劃、洪水預(yù)報、河道演變規(guī)律及水資源綜合利用等，一些科研單位開發(fā)了采用半隱式有限元方法，建立了珠江三角洲河網(wǎng)水流、泥沙數(shù)學模型，模型具有計算穩(wěn)定性強、收斂快的特點［20］；采用有限元的桿系數(shù)據(jù)結(jié)構(gòu)，對于大型河網(wǎng)計算中增加河道和斷面數(shù)，數(shù)據(jù)準備簡單。模型經(jīng)1999年洪季和2001年枯季兩次大規(guī)模、長系列的水文泥沙資料的檢驗，各站水位、流量和含沙量計算結(jié)果與實測過程十分吻合。目前該河網(wǎng)模型在工程實際中得到了廣泛應(yīng)用。

2.2 并行計算技術(shù)的發(fā)展

并行計算（多個CPU同時工作完成一項計算）技術(shù)已經(jīng)廣泛應(yīng)用在航空、新材料研發(fā)、石油勘探、氣象預(yù)測等領(lǐng)域。水運行業(yè)已于2007年將最廣泛采用的并行實現(xiàn)方法——基于消息傳遞（MPI）的并行編程方法引入內(nèi)河港口與航道水流泥沙數(shù)值模擬中，成功開發(fā)了內(nèi)河水沙數(shù)值模擬軟件TK-2DC［43］。同時也開發(fā)出了基于共享內(nèi)存（OpenMP）并行編程方法的河道水沙數(shù)值模擬軟件Hydroinfo。

內(nèi)河并行計算技術(shù)在長江下游江心洲至烏江長河段航道系統(tǒng)整治前期工作中得到了應(yīng)用，并在長江下游安慶至南京長河段航道整治與維護計算平臺的開發(fā)中得到了應(yīng)用，為長河段水沙輸移規(guī)律及航道系統(tǒng)整治相關(guān)研究提供了技術(shù)支持，目前已經(jīng)服務(wù)于航道整治與維護工作中［43］。

2.3 模式分離技術(shù)的發(fā)展應(yīng)用

為了提高大尺度河道水沙三維計算數(shù)值穩(wěn)定性和縮短計算時間，在應(yīng)用并行計算方法的同時，物理海洋學中的模式分離技術(shù)已引入三維模擬中的自由表面計算。利用外模態(tài)和內(nèi)模態(tài)分別計算水位場和垂向流場，對于沒有劇烈波動的水流十分有效，避免自由表面跟蹤大量的計算，該技術(shù)應(yīng)用能揭示水流內(nèi)部結(jié)構(gòu)，又能較好地模擬自由表面［38-41］。

模式分離技術(shù)和非靜壓假設(shè)技術(shù)在河道水流泥沙三維數(shù)值模擬中已廣泛應(yīng)用。在長江上游匯合口河段航道整治、整治建筑物水毀機理的研究、橋墩三維水流泥沙數(shù)值模擬等方面取得了較好的應(yīng)用。

2.4 非靜壓假設(shè)下三維水流數(shù)學模型

采用半隱、分步的求解模式建立了模擬三維具有自由表面流動的非靜壓數(shù)值模型。模型中首先求解靜壓假設(shè)下的三維Navier-Stokes方程，然后計算非靜壓項，對靜壓的計算結(jié)果進行非靜壓修正以得到最終結(jié)果［18］。采用有限體積方法離散控制方程。采用滿足自由表面運動學邊界條件的水位控制方程來計算新時刻自由表面的位置。在表層采用對垂向動量方程積分的方式確保了自由表面處壓力邊界條件的精確給出。對流項和水平粘性項的計算采用半解析的semi-Lagrangian方法。通過封閉κ-ε紊流模型，形成了完整的三維流動數(shù)學模型，使模型具備模擬較復(fù)雜的三維水流泥沙的能力。

3 軟件化系統(tǒng)集成

河流水流泥沙數(shù)值模擬軟件已經(jīng)開始向軟件化發(fā)展，目前已經(jīng)開發(fā)了如TK-2DC［36］、Hydroinfo等軟件，這些軟件也在向品牌化、商品化方向發(fā)展。

3.1 前處理軟件系統(tǒng)開發(fā)

由于水運工程模擬內(nèi)河水域的邊界條件和地形變化相差很大，加之模擬方法的不同，在數(shù)值模擬中出現(xiàn)了多種多樣的網(wǎng)格剖分方法。目前已經(jīng)開發(fā)出二、三維正交曲線網(wǎng)格、矩形網(wǎng)格、非結(jié)構(gòu)四邊形和三角形網(wǎng)格生成軟件系統(tǒng)［36］。

3.2 后處理軟件系統(tǒng)開發(fā)

數(shù)值模擬結(jié)果需借助計算機仿真技術(shù)直觀展示。目前對于一維、二維和一些三維模擬結(jié)果已經(jīng)開發(fā)出二維和準三維后處理展示系統(tǒng)。對三維模擬數(shù)據(jù)的后處理，開發(fā)了能夠展示多種場景的軟件系統(tǒng)，包括表面場、截面流場、真三維立體電影方式。

3.3 河流水流泥沙數(shù)值模擬通用化軟件系統(tǒng)集成

河流水流泥沙數(shù)值模擬技術(shù)發(fā)展較快，已經(jīng)初步形成了一些通用化的軟件系統(tǒng)（如TK系列軟件），這些軟件包括完整的前處理、計算核心、后處理等模塊，已經(jīng)逐步從通用化、實用化、軟件化、可視化向品牌化、商品化方向發(fā)展［16-17］。

4 結(jié)語

隨著泥沙運動基本理論和河床演變等基礎(chǔ)理論學科的不斷發(fā)展，河流水流泥沙數(shù)值模擬技術(shù)已取得了很大的進展，但由于河流泥沙問題的復(fù)雜性，仍有很多問題待研究和解決。

（1）基礎(chǔ)理論的創(chuàng)新。數(shù)值模擬是依據(jù)對物理概念的合理科學認識基礎(chǔ)上的，只有理論上的創(chuàng)新才會帶動數(shù)值模擬基本上的技術(shù)創(chuàng)新。加強泥沙數(shù)值模擬中理論性關(guān)鍵問題的研究，如邊界阻力、紊動粘性系數(shù)、泥沙擴散系數(shù)、底部泥沙邊界的給定等。

（2）觀測技術(shù)的提高和觀測資料的系統(tǒng)化與公開化。數(shù)值模擬計算是以試驗和天然觀測資料為基礎(chǔ)的，是檢驗?zāi)M成果的重要依據(jù)，而當前不同觀測手段的差異及誤差仍制約著理論研究的深度與相關(guān)研究成果的可靠性。同時，觀測資料的系統(tǒng)化和公開化也是推進數(shù)值模擬技術(shù)發(fā)展的重要助力。

（3）數(shù)學模擬成套技術(shù)的標準化。河流水流泥沙數(shù)值模擬在強化基本物理圖景的同時，應(yīng)該使數(shù)值模擬成套技術(shù)向標準化、可視化和商業(yè)化發(fā)展，實現(xiàn)模型的開放性檢驗與應(yīng)用。

［1］胡春宏，王延貴，張燕菁.河流泥沙模擬技術(shù)進展與展望［J］.水文，2006，26（3）：37-41.HU C H，WANG Y G，ZHANG Y J.Development and Prospects of Simulating Technology for River Sedimentation［J］.Journal of China Hydrology，2006，26（3）：37-41.

［2］趙沖久.我國水運工程科學研究的現(xiàn)狀［J］.水道港口，2008，29（4）：229-232.ZHAO C J.Brief introduction of water transport engineering research of China［J］.Journal of Waterway and Harbor，2008，29（4）：229-232.

［3］王光謙.河流泥沙研究進展［J］.泥沙研究，2007（2）：64-81.WANG G Q.Advances in river sediment research［J］.Journal of Sediment Research，2007（2）：64-81.

［4］王光謙.中國泥沙研究述評［J］.水科學進展，1999，10（3）：337-344.WANG G Q.A Review of Sediment Research in China［J］.Advances in Water Science，1999，10（3）：337-344.

［5］曹叔尤.水力學及河流動力學基本問題研究的現(xiàn)狀與任務(wù)［J］.四川大學學報：工程科學版，2002（1）：1-5.CAO S Y.Status Quo and Tasks of Basic Research in Hydraulics and Fluvial Dynamics［J］.Journal of Sichuan University，2002（1）：1-5.

［6］萬遠揚，金中武，黃仁勇.泥沙模型研究述評與前景展望［J］.南水北調(diào)與水利科技，2006，4（1）：48-52.WAN Y Y，JIN Z W，HUANG R Y.Outlook and Prospect on Sedimentation Model Study［J］.South-to-North Water Transfers and Water Science&Technology，2006，4（1）：48-52.

［7］王兆印.泥沙研究的發(fā)展趨勢和新課題［J］.地理學報，1998，53（3）：245-255.WANG Z Y.Outlook for Sedimentation Research［J］.ACTA Geographica SINICA，1998，53（3）：245-255.

［8］陳雄波，俞國青，唐洪武.河道水流泥沙運動基本理論研究進展［J］.水利水電科技進展，2003，23（4）：59-64.CHEN X B，YU G Q，TANG H W.Advances on River Sedimentation Research［J］.Advances in Science and Technology of Water Resources，2003，23（4）：59-64.

［9］竇希萍，張幸農(nóng).南科院近 10 年泥沙研究工作綜述［J］.泥沙研究，1999（6）：17-20.DOU X P，ZHANG X N.Reviews of Sedimentation Research by NHRI in Recent Ten Years［J］.Journal of Sediment Research，1999（6）：17-20.

［10］胡春宏.我國江河治理與泥沙研究展望［J］.水利水電技術(shù)，2001（1）：50-52.HU C H.Outlook on River Regulation and Sediment Research［J］.Water Resources and Hydropower Engineering，2001（1）：50-52.

［11］劉興年，曹叔尤，方鐸.90 年代四川大學泥沙研究回顧與瞻望［J］.泥沙研究，1999（6）：33-36.LIU X N，CAO S Y，F(xiàn)ANG D.Reviews and Prospects of Sedimentation Research by Sichuan University in the 1990s［J］.Journal of Sediment Research，1999（6）：33-36.

［12］韓其為，胡春宏.50年來泥沙研究所主要研究進展［J］.中國水利水電科學研究院學報，2008，6（3）：170-182.HAN Q W，HU C H.Main research development by the Department of Sediment Research in the past 50 years［J］.Journal of China Institute of Water Resources and Hydropower Research，2008，6（3）：170-182.

［13］王光謙，張仁，惠遇甲，等.清華大學泥沙研究回顧［J］.泥沙研究，1999（6）：21-25.WANG G Q，ZHANG R，HUI Y J.Reviews of Sedimentation Research by Tsinghua University［J］.Journal of Sediment Research，1999（6）：21-25.

［14］談廣鳴，陳立.河床變形混交模型預(yù)測技術(shù)及其進展［J］.水利水電科技進展，2001（4）：14-19.TAN G M，CHEN L.Prediction of Riverbed Deformation Mix-Model and Its Development［J］.Advances in Science and Technology of Water Resources，2001（4）：14-19.

［15］謝鑒衡.江河演變與治理研究［M］.武漢：武漢大學出版社，2004.

［16］李孟國.海岸河口泥沙數(shù)學模型研究進展［J］.海洋工程，2006，24（1）：139-154.LI M G.A review on mathematical models of sediment in coastal and estuarine waters［J］.Ocean Engineering，2006，24（1）：139-154.

［17］秦崇仁，趙沖久.海岸泥沙運動研究綜述和展望［J］.水道港口，2003，24（2）：65-67.QIN C R，ZHAO C J.Review and Forecast of Coastal Sand Transport Research［J］.Journal of Waterway and Harbor，2003，24（2）：65-67.

［18］金生，倪漢根.水流和泥沙問題數(shù)值模擬與工程應(yīng)用［J］.大連理工大學學報，2000，40（S1）：78-82.JIN S，NI H G.Numerical analysis of flow and sediments and engineering application［J］.Journal of Dalian University of Technology，2000，40（S1）：78-82.

［19］胡海明，趙良奎，楊國錄.一維泥沙數(shù)學模型的發(fā)展方向［J］.水利電力科技，1995，22（1）：33-40.HU H M，ZHAO L K，YANG G L.The Forecast on 1-D Sediment Mathematical model［J］.Water Resources and Electric Power，1995，22（1）：33-40.

［20］張華慶，金生，沈漢堃，等.珠江三角洲河網(wǎng)非恒定水沙數(shù)學模型研究［J］.水道港口，2004，25（3）：121-128.ZHANG H Q，JIN S，SHEN H K，et al.Mathematical Model of Unsteady Flow and Sediment in Pearl River Delta［J］.Journal of Waterway and Harbor，2004，25（3）：121-128.

［21］趙明登，李義天.二維泥沙數(shù)學模型及工程應(yīng)用問題探討［J］.泥沙研究，2002（1）：66-70.ZHAO M D，LI Y T.Study on Two-Dimensional Sediment Mathematical Models and Problems in Application［J］.Journal of Sediment Research，2002（1）：66-70.

［22］陳界仁，沙撈·巴里，陳國祥.二維水庫水流泥沙數(shù)值模擬［J］.河海大學學報，2000，28（5）：11-15.CHEN J R，Salau O B E，CHEN G X.Two Dimensional Numulation of Flow and Sediment Transport in Reservoir［J］.Journal of Hohai University，2000，28（5）：11-15.

［23］夏軍強，王光謙.剖面二維懸移質(zhì)輸移方程的分步解法［J］.長江科學院院報，2000，17（2）：14-17.XIA J Q，WANG G Q.Fractional steps method for solving vertical 2-D suspended sediment transport equation［J］.Journal of Yangtze River Scientific Research Institute，2000，17（2）：14-17.

［24］余明輝，楊國錄.平面二維非均勻沙數(shù)值模擬方法［J］.水利學報，2000（5）：65-69.YU M H，YANG G L.Study on mathematical model for transportation of plane 2-D non-homogeneous sediment［J］.Journal of Hydraulic Engineering，2000（5）：65-69.

［25］張細兵，殷瑞蘭.平面二維水流泥沙數(shù)值模擬［J］.水科學進展，2002，13（6）：665-669.ZHANG X B，YIN R L.Planar 2-D flow and sediment mathematical modeling［J］.Adances in Water Science，2002，13（6）：665-669.

［26］張細兵，董耀華，殷瑞蘭.河道平面二維水沙數(shù)學模型的有限元法［J］.泥沙研究，2002（6）：60-65.ZHANG X B，DONG Y H，YIN R L.Finite Element Method for 2-D Flow and Sediment Mathematical Model in Channels［J］.Journal of Sediment Research，2002（6）：60-65.

［27］韋直林.河道水流泥沙問題的一種有限元解法［J］.武漢水利電力學院學報，1990，23（6）：77-86.WEI Z L.Finite Element Method to Solve River Flow and Sediment Problems［J］.Journal of Wuhan Hydraulic&Electric University，1990，23（6）：77-86.

［28］董文軍.永定新河二維輸沙有限元數(shù)值模擬［J］.泥沙研究，1996（4）：86-93.DONG W J.The 2-D Finite Element Method for Yongding New River Sediment Transport［J］.Journal of Sediment Research，1996（4）：86-93.

［29］郭慶超，韓其為，何明民.二維潮流及泥沙數(shù)學模型［J］.泥沙研究，1996（1）：48-55.GUO Q C，HAN Q W，HE M M.The 2-D Mathematical Model for Tide and Sediment Transport［J］.Journal of Sediment Research，1996（1）：48-55.

［30］張華慶，陸永軍.河道平面二維泥沙數(shù)學模型［J］.水道港口，1993（2）：36-43.ZHANG H Q，LU Y J.The 2-D Mathematical Model for River Sediment Transport［J］.Journal of Waterway and Harbor，1993（2）：36-43.

［31］張華慶，樂培九，楊細根.二維泥沙數(shù)學模型的改進——模型的建立［J］.水道港口，1996（4）：9-19.ZHANG H Q，YUE P J，YANG X G.Improvement on 2-D Mathematical Model of Flow and Sediment—Establishment of Model［J］.Journal of Waterway and Harbor，1996（4）：9-19.

［32］張華慶，樂培九，楊細根.擬合坐標系下平面二維水流泥沙數(shù)學模型及其應(yīng)用［J］.水利學報，1998（3）：39-47.ZHANG H Q，YUE P J，YANG X G.The Application of 2-D Mathematical Model for River Sediment Transport［J］.Journal of Hydraulic Engineering，1998（3）：39-47.

［33］白玉川，楊建民，黃本勝.二維水沙數(shù)學模型在復(fù)雜河道治理中的應(yīng)用［J］.水利學報，2003（9）：25-30.BAI Y C，YANG J M，HUANG B S.Application of 2-D mathematical model in training of complicated river channel［J］.Journal of Hydraulic Engineering，2003（9）：25-30.

［34］張紅武，趙連軍，王光謙，等.黃河下游河道準二維泥沙數(shù)學模型研究［J］.水利學報，2003（4）：1-7.ZHANG H W，ZHAO L J，WANG G Q，et al.Quasi-2D mathematical model for morphology of the Lower Yellow River［J］.Journal of Hydraulic Engineering，2003（4）：1-7.

［35］譚維炎，胡四一.二維淺水流動的一種普適的高性能格式——有限體積Osher格式［J］.水科學進展，1991，2（3）：154-161.TAN W Y，HU S Y.A Universal high-performance algorithm for two-dimensional unsteady shallow water flow computation-the finite-volume osher scheme［J］.Advances in Water Science，1991，2（3）：154-161.

［36］張明進，楊星.TK-2DC 軟件前處理模塊－貼體正交網(wǎng)格生成軟件開發(fā)［J］.水道港口，2008，29（2）：124-130.ZHANG M J，YANG X.Generation of orthogonal body-fitted meshes and software development［J］.Journal of Waterway and Harbor，2008，29（2）：124-130.

［37］賴錫軍，曲卓杰，周杰，等.非結(jié)構(gòu)網(wǎng)格上的三維淺水流動數(shù)值模型［J］.水科學進展，2006，17（5）：693-699.LAI X J，QU Z J，ZHOU J，et al.3-D hydrodynamic model for shallow water on unstructured grids［J］.Advances in Water Science，2006，17（5）：693-699.

［38］陸永軍，竇國仁，韓龍喜，等.三維紊流懸沙數(shù)學模型及應(yīng)用［J］.中國科學：E 輯，2004，34（3）：311-328.LU Y J，DOU G R，HAN L X，et al.Three-dimensional Mathematical Model of Turbulent Flow and Sediment Transport and Its Applications［J］.Science in China：Series E，2004，34（3）：311-328.

［39］陳國祥，陳界仁，沙撈·巴里.三維泥沙數(shù)學模型的研究進展［J］.水利水電科技進展，1998，18（1）：13-20.CHEN G X，CHEN J R，Salau O B E.Review of Development of 3D Sediment Mathematical Model［J］.Advances in Science and Technology of Water Resources，1998，18（1）：13-20.

［40］張新周，竇希萍，王惠民，等.三維水沙數(shù)值模擬中紊流隨機理論應(yīng)用初探［J］.泥沙研究，1998（3）：13-19.ZHANG X Z，DOU X P，WANG H M，et al.Turbulence stochastic model and its appliction in 3D numerical simulation of flow and sediment［J］.Journal of Sediment Research，1998（3）：13-19.

［41］李紹武，王建軍，郝品正，等.河道三維水流數(shù)值模型及應(yīng)用［J］.港工技術(shù)，2008（1）：1-4.LI S W，WANG J J，HAO P Z，et al.A 3-D Numerical Open Channel Flow Model and its Application［J］.Port Engineering Technology，2008（1）：1-4.

［42］李冰凍，李克鋒，李嘉，等.天然河段港口工程三維流場數(shù)值模擬［J］.水利水運工程學報，2006（2）：27-31.LI B D，LI K F，LI J，et al.Numerical simulation of 3-D flow for harbor engineering in natural rivers［J］.Hydro-Science and Engineering，2006（2）：27-31.

［43］王建軍，張明進.河道二維水沙數(shù)學模型并行計算技術(shù)研究［J］.水道港口，2009，30（3）：222-225.WANG J J，ZHANG M J.Parallel computer technology study on 20 numerical model of flow and sediment in rivers［J］.Journal of Waterway and Harbor，2009，30（3）：222-225.