水利期刊網(wǎng)頁(yè)制作格式說(shuō)明

1、<p>　　文章編號(hào)：1672-3031(2012)03-0000-00</p><p><b>　　題目</b></p><p>　　姓名1,2，姓名1,2，姓名2</p><p>　　(1.單位，北京 100871； 2. 單位， 北京 100871)</p><p>　　摘要：泥沙運(yùn)動(dòng)特征的變化及其對(duì)港口

2、航道的影響是大型海岸工程實(shí)施時(shí)需要考慮的主要問(wèn)題。由于海岸生態(tài)系統(tǒng)往往對(duì)泥沙沖淤的變化非常敏感，生態(tài)評(píng)估也逐漸成為工程規(guī)劃及方案可行性論證研究的重要內(nèi)容。本文利用泥沙數(shù)學(xué)模型，以深圳灣為例研究了5種填海方案造成的海灣來(lái)沙量、含沙量和泥沙沖淤分布的變化，討論了泥沙淤積速率的變化對(duì)港口航道和紅樹(shù)林生態(tài)系統(tǒng)的影響，并通過(guò)岸線方案的對(duì)比指出了在分別考慮各種泥沙評(píng)價(jià)指標(biāo)時(shí)對(duì) 應(yīng)的五種填海方案的優(yōu)劣。定量分析表明，不同評(píng)價(jià)指標(biāo)對(duì)應(yīng)的方案優(yōu)劣次序有所

3、差異，協(xié)調(diào)多種指標(biāo)要求的綜合優(yōu)化岸線可以通過(guò)比較各種方案后給出。</p><p>　　關(guān)鍵詞：泥沙數(shù)學(xué)模型；填海工程；岸線；深圳灣</p><p>　　中圖分類(lèi)號(hào)： 文獻(xiàn)標(biāo)識(shí)碼：A </p><p><b>　　1 引言</b></p><p>　　填海

4、對(duì)當(dāng)?shù)鼗蛑苓叺貐^(qū)潮波特性的影響，已經(jīng)在許多填海工程的岸線確定中得到充分重視［1］。近年來(lái)，填海工程造成的區(qū)域泥沙運(yùn)動(dòng)特征的改變以及對(duì)港口、航道和海岸生態(tài)系統(tǒng)的影響也愈來(lái)愈受到廣泛的關(guān)注。泥沙淤積帶來(lái)的問(wèn)題尤其復(fù)雜。例如，香港維多利亞港附近由于填海造地使得一些區(qū)域潮流作用減弱，形成了有利于泥沙沉降的環(huán)境，沉積速率明顯增大［2］；廈門(mén)西海域的海岸工程和圍墾造地使納潮量減少了67%左右，潮汐對(duì)泥沙的輸運(yùn)能力因此大大減弱，使廈門(mén)港的淤積加重、通

5、航條件受到影響，每年必須花800多萬(wàn)元費(fèi)用清理航道；此外，沉積環(huán)境的改變還導(dǎo)致珍貴的文昌魚(yú)幾乎絕跡，海豚減少［3］。</p><p>　　過(guò)去關(guān)于填海后泥沙運(yùn)動(dòng)特征變化及其對(duì)港口航道影響的研究較多［4-6］，對(duì)海岸生態(tài)系統(tǒng)影響的研究卻較少。本文以深圳灣填海工程為例，嘗試?yán)媚嗌硵?shù)學(xué)模型對(duì)填海前后海灣的來(lái)沙量、沖淤分布規(guī)律、港口航道和生態(tài)保護(hù)區(qū)的影響進(jìn)行綜合考慮，通過(guò)對(duì)比分析典型填海岸線方案的優(yōu)劣為填海工程的優(yōu)化提

6、供依據(jù)。</p><p><b>　　2 研究區(qū)域概況</b></p><p>　　深圳灣是深圳和香港的交界水域，為珠江口東岸的一個(gè)半封閉型淺水灣(如圖1)，面積約為90km2，平均水深2.9m。注入深圳灣的河流中，深圳河是最大的河流。受南海不規(guī)則半日混合潮影響，海灣潮流方向?yàn)橥鶑?fù)流，年均潮差1.37m。根據(jù)近來(lái)的實(shí)測(cè)資料［7］和水動(dòng)力學(xué)模型計(jì)算結(jié)果，漲潮時(shí)深圳灣內(nèi)大

7、潮期流速為0.39～0.65m/s，小潮期為0.31～0.40m/s；落潮時(shí)，深圳灣內(nèi)大潮期流速為0.58～0.74m/s，小潮期為0.20～0.34m/s。</p><p>　　深圳灣具有重要的航運(yùn)功能，灣內(nèi)有蛇口港、東角頭港等主要港口。在灣頂附近的福田國(guó)家紅樹(shù)林自然保護(hù)區(qū)和米埔自然保護(hù)區(qū)分布著大面積的泥灘和紅樹(shù)林，是華南地區(qū)最重要的濕地生態(tài)系統(tǒng)。該區(qū)域?yàn)樵S多珍稀和瀕危生物物種提供了棲息地，尤其是以水鳥(niǎo)和候鳥(niǎo)為

8、主的鳥(niǎo)類(lèi)，其中包括一些在</p><p>　　收稿日期： 基金項(xiàng)目： </p><p>　　作者簡(jiǎn)介：第一作者姓名（出生年-），性別，籍貫，學(xué)歷，職稱(chēng)，主要從事（研究方向）研究。E_mail: </p><p>　　通信作者：姓名（出生年-），性別，籍貫，學(xué)歷，職稱(chēng)，主要從事（研究方向）研究。E_mail:</p><p>　　全球范圍內(nèi)受

9、到威脅的種類(lèi)［8］。</p><p>　　計(jì)劃在深圳灣北岸進(jìn)行的大規(guī)模填海造地工程，不可避免地會(huì)改變泥沙的運(yùn)動(dòng)規(guī)律，并對(duì)港口、航道和自然保護(hù)區(qū)造成影響。為減少填海對(duì)港口、航道的負(fù)面影響和避免生態(tài)環(huán)境惡化，有必要對(duì)深圳灣的填海工程的設(shè)計(jì)方案進(jìn)行全面規(guī)劃和綜合論證。擬議中的深圳灣北岸填海工程岸線方案共有五種(見(jiàn)圖)1，方案一至方案五的填海面積分別為8.8、13.2、10. 95、5.88、8.40km2。文獻(xiàn)［1］利

10、用水動(dòng)力學(xué)數(shù)學(xué)模型，預(yù)測(cè)了深圳灣填海后納潮量、流場(chǎng)、潮位等海灣水動(dòng)力條件的變化，并討論了流速、潮位變化對(duì)防洪、港口航道和紅樹(shù)林自然保護(hù)區(qū)的影響。本文將在水動(dòng)力學(xué)數(shù)學(xué)模型研究的基礎(chǔ)上，進(jìn)一步探討填海前后泥沙運(yùn)動(dòng)特征的變化及其影響。</p><p><b>　　3 泥沙模型 </b></p><p>　　3.1 基本方程 根據(jù)實(shí)測(cè)資料［7］分析，冬季深圳灣各測(cè)站實(shí)測(cè)平均

11、含沙量在0.0073～0.205kg/m3之間，而夏季在0.0137～0.2385kg/m3之間。由于深圳灣水體中含沙量不高，泥沙對(duì)水流的影響極小，泥沙數(shù)學(xué)模型與水動(dòng)力學(xué)模型可以采用非耦合的方式求解。深圳灣的床沙組成為砂層及亞黏土層，上覆海相淤泥層，底質(zhì)細(xì)、中值粒徑為0.004～0.006mm［9,10］，泥沙運(yùn)動(dòng)以懸移質(zhì)運(yùn)動(dòng)為主，推移質(zhì)運(yùn)動(dòng)可以忽略。由于缺乏足夠的懸移質(zhì) 泥沙級(jí)配資料，本研究將懸沙按均勻沙處理。懸移質(zhì)非平衡輸沙方程和河

12、床變形方程分別如下［11］</p><p><b>　　(1)</b></p><p><b>　　(2)</b></p><p>　　式中：S為垂線平均含沙量；u、v分別為沿x和y方向的垂線平均流速；α為恢復(fù)飽和系數(shù)；ω為懸沙沉降速度;H為實(shí)際水深；φ為水流挾沙力；KxKy分別為沿x和y方向的泥沙擴(kuò)散系數(shù)；η為河床變形厚

13、度(η<0為沖刷)；γm為泥沙干容重。水流挾沙力公式采用常見(jiàn)的半經(jīng)驗(yàn)公式［12］，即</p><p><b>　　(3)</b></p><p>　　式中：A和m為經(jīng)驗(yàn)系數(shù)；V為流速大??；g為重力加速度。</p><p>　　深圳灣的地形如圖1。泥沙模型的計(jì)算范圍從深圳河口到深圳灣口赤灣爛角嘴一線，共有168×132個(gè)網(wǎng)格，網(wǎng)格

14、大小為100m×100m，模型中的流速和實(shí)際水深，直接引用水動(dòng)力學(xué)模型計(jì)算結(jié)果。差分格式采用ADI（Alternative Direction Implicit）法［13］構(gòu)造，即計(jì)算時(shí)段分為兩個(gè)半時(shí)間步長(zhǎng)，在前半個(gè)時(shí)間步長(zhǎng)內(nèi)，只對(duì)空間x方向作變量運(yùn)算，而后半個(gè)步長(zhǎng)內(nèi)，只對(duì)空間y方向作變量運(yùn)算。差分格式中對(duì)流項(xiàng)采用迎風(fēng)格式。時(shí)間步長(zhǎng)經(jīng)調(diào)試取100s。計(jì)算域開(kāi)邊界條件根據(jù)泥沙濃度的監(jiān)測(cè)資料獲得；對(duì)于閉邊界，假設(shè)，為閉邊界外法線方

15、向的單位矢量；泥沙方程的初始值由實(shí)測(cè)含沙量給定。</p><p>　　模型驗(yàn)證 模型參數(shù)的率定采用1998年1月的實(shí)測(cè)水文資料［7］進(jìn)行。在深圳灣的水文流監(jiān)測(cè)點(diǎn)位中(如圖1)，VS3和VS4接近灣口開(kāi)邊界，利用這兩點(diǎn)的泥沙濃度時(shí)間序列，通過(guò)線性插值得到灣口開(kāi)邊界；河口的泥沙</p><p>　　假設(shè)枯水期(或豐水期)每個(gè)月的來(lái)沙量都等于模型預(yù)測(cè)的枯水期(或豐水期)月來(lái)沙量，以枯水期5個(gè)月

16、、豐水期7個(gè)月將12個(gè)月的來(lái)沙量累加，可以近似得到深圳灣的年來(lái)沙量。如表2所示，除方案二對(duì)應(yīng)的年來(lái)沙量減少外，其它方案的年來(lái)沙量都略有增加，其中方案四增加得最多，而方案三的年來(lái)沙量變化最小。因此，單就進(jìn)出深圳灣泥沙變化最小的要求來(lái)看，各方案的優(yōu)劣順序?yàn)椋悍桨溉⒁?、五、二、四?lt;/p><p>　　表2 填海前后深圳灣月來(lái)沙量和年來(lái)沙量變化 (單位：103t)</p>

17、<p>　　對(duì)各種填海岸線方案的影響進(jìn)行的分析表明：不同的考慮對(duì)象要求的泥沙評(píng)價(jià)指標(biāo)及大小不同，因而對(duì)應(yīng)岸線方案的優(yōu)劣次序排列也有所差異(如表4)，所以不能憑單個(gè)指標(biāo)來(lái)確定填海方案。綜合考慮各種泥沙評(píng)價(jià)指標(biāo)，可以認(rèn)為方案三、四較優(yōu)。</p><p><b>　　表4 方案優(yōu)劣比較</b></p><p>　　水動(dòng)力學(xué)模型和泥沙數(shù)學(xué)模型為評(píng)估填海前后水動(dòng)力學(xué)條

18、件和泥沙運(yùn)動(dòng)特征變化提供了基礎(chǔ)。填海工程對(duì)海灣水質(zhì)、環(huán)境容量和污染物排放方式方面的影響需要利用其它模型。</p><p><b>　　4 結(jié)論參考文獻(xiàn)：</b></p><p>　　林秉南，趙雪華，施麟寶. 河口建壩對(duì)毗鄰海灣潮波影響的計(jì)算(二維特征線理論法)[J]. 水利學(xué)報(bào)，1980(3):16-25.</p><p>　　潘少明，施曉冬，

19、王建業(yè)，等. 圍海造地工程對(duì)香港維多利亞現(xiàn)代沉積作用的影響[J]. 沉積學(xué)報(bào)，2000，18(1):22-28. </p><p>　　洪華生，陳宗團(tuán). 海岸帶綜合管理中面臨的科學(xué)問(wèn)題[J]. 海洋管理，1998(1)：28-31. </p><p>　　松林，丁平興. 湛江灣沿岸工程沖淤影響的預(yù)測(cè)分析Ⅱ. 沖淤的數(shù)值計(jì)算[J]. 海洋學(xué)報(bào)，1997，19(1):64-72. </p

20、><p>　　Han Z C, Cheng H P. Two-dimensional sediment mathematical model of Hangzhou Bay. Proc. 3rd Intern. Symp. on River Sedimentation[J]. The University of Mississippi, USA, 1986, 463-471. </p><p>

21、;　　Chien N, Wan Z H. Mechanics of Sediment Transport[D]. American Society of Civil Engineers, Virginia: ASCE Press, 1999. </p><p>　　Leendertse J J. Aspects of SYMSYS2D: a system for two dimensional flow comp

22、utation[M]. The Rand Corporation, 1987.</p><p>　　隋淑珍，張喬民. 華南沿海紅樹(shù)林海岸沉積物特征分析[J]. 熱帶海洋，1999，18(4):18-23.</p><p>　　倪晉仁，楊小毛，王光謙. 深港交界帶經(jīng)濟(jì)開(kāi)發(fā)過(guò)程中泥沙對(duì)生態(tài)環(huán)境的影響[J]. 地理學(xué)報(bào)，1998，53(4):350-355.</p><p&g

23、t;　　張喬民，于紅兵，陳欣樹(shù)，等. 紅樹(shù)林生長(zhǎng)帶與潮汐水位關(guān)系的研究[J]. 生態(tài)學(xué)報(bào)，1997, 17(3)：258-265.</p><p><b>　　題目（英文）</b></p><p>　　姓名（英文）1，2, 姓名（英文）1，2, 姓名（英文）2</p><p>　　(1.單位（英文）， Beijing 100871,China;

24、 2.單位（英文）， Beijing 100871,China)</p><p>　　Abstract: Attention should be paid on variations of sediment transport characteristics and their impacts on ports and navigation channels before reclaiming a bay. In

25、 general, a bay system would be sensitive to the change of sediment erosion or deposition, and thus the ecological assessment in terms of sediment models could be important for the feasibility study of the proposed optio