長河段河床縱剖面分形特征研究

倪志輝，王 偉，吳立春

(1. 重慶交通大學(xué) 水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室 國家內(nèi)河航道整治工程技術(shù)研究中心，重慶 400074；2. 重慶交通大學(xué) 西南水運(yùn)工程科學(xué)研究所，重慶 400016；3.重慶第二師范學(xué)院，重慶 400067)

?

長河段河床縱剖面分形特征研究

倪志輝1,2，王 偉1，吳立春3

(1. 重慶交通大學(xué) 水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室 國家內(nèi)河航道整治工程技術(shù)研究中心，重慶 400074；2. 重慶交通大學(xué) 西南水運(yùn)工程科學(xué)研究所，重慶 400016；3.重慶第二師范學(xué)院，重慶 400067)

以長江干流水富至朱沱段為研究對(duì)象，建立一維非恒定分形數(shù)學(xué)模型，將研究河段分作3段，以岷江、沱江匯流口作為分界，計(jì)算各分段河道縱剖面分維數(shù)。研究表明：縱剖面分維數(shù)呈二階分維特性，且整體分維數(shù)小于分段分維數(shù)；整體分維數(shù)不等于各分段分維數(shù)算術(shù)平均值；縱剖面分維數(shù)隨流量增大而不斷減?。桓鞣侄魏拥揽v坡降隨縱剖面分形維數(shù)的變化形態(tài)類似波形。

航道工程；縱剖面；非恒定流；分形維數(shù)；變維分形；數(shù)值模擬

0 引 言

樞紐下泄的非恒定流使下游水流條件變得異常復(fù)雜，使下游河道的天然水沙過程發(fā)生了較大改變，將在一定程度上影響下游河道功能的正常發(fā)揮。目前針對(duì)樞紐下游長河段縱剖面分形特征的研究還很少，并且還處于理論探索分析階段，因此樞紐下游長河段縱剖面分維數(shù)的研究對(duì)下游城市的發(fā)展有著重要的戰(zhàn)略和現(xiàn)實(shí)意義。

20世紀(jì)70年代，本華·曼德博在《科學(xué)》上發(fā)表了《英國海岸線有多長》的文章，分形理論被提出，其作為一種新的概念和方法，正在許多領(lǐng)域開展應(yīng)用探索。許多學(xué)者應(yīng)用分形理論研究流域特征，對(duì)河流的分形研究獲得了大量成果[1-18]。

蜿蜒性是河流的重要特征，在研究河流的蜿蜒曲折程度等不規(guī)則性形態(tài)特征過程中誕生了河流的分形維數(shù)這個(gè)特征量，河流的分維由河長分維和河網(wǎng)分維兩種組成[19]。而河長的分維又分為整體分維和局部分維。對(duì)于整條河流分維的探討，即為整體分維，而將河流分段，探討各分段的分維，即為分段分維(局部分維)。河流的整體分維和局部分維是河流形態(tài)的重要特征量，二者有著緊密聯(lián)系。

樞紐下泄非恒定流中，水位、流量的頻繁急劇變化加劇了航道與岸灘的不穩(wěn)定性，且由于影響因素眾多，因此筆者利用分形維數(shù)代替多個(gè)影響因素，建立數(shù)學(xué)模型。研究對(duì)象為水富(橫江匯流口)至朱沱水文站共長約270 km的河段。如圖1，根據(jù)大量實(shí)測地形資料，建立一維非恒定分形數(shù)學(xué)模型，采用二階累計(jì)和變維分形方法，計(jì)算壩下長河段河道縱剖面分形維數(shù)，并初步探討河流長度分段分維數(shù)與整體分維數(shù)，坡降、流量與縱剖面分維數(shù)的相互關(guān)系。

圖1 水富至朱沱段河道示意Fig. 1 River course schematic of the river section from Shuifu to Zhutuo

1 研究方法

1.1 基本方程

一維水流運(yùn)動(dòng)和連續(xù)方程：

(1)

(2)

式中：A為過水面積；Q為流量；U為斷面平均流速；Z為水位；B為水面寬度；x為河段水平方向，沿下游為正；t為時(shí)間；ρ為水的密度；τb為河底切應(yīng)力。

由于所研究的問題為長河段、長時(shí)段內(nèi)發(fā)生的，在實(shí)際計(jì)算中對(duì)方程組進(jìn)行簡化。將整個(gè)計(jì)算時(shí)段劃分為小的計(jì)算時(shí)段，將長河段劃分為若干個(gè)短河段，可按恒定非均勻流考慮，水流運(yùn)動(dòng)變?yōu)?/p>

(3)

(4)

式中：Z2、Z1為計(jì)算段上、下游斷面水位；V2、V1為計(jì)算段上、下游斷面平均流速；α2、α1為計(jì)算段上、下游斷面的動(dòng)能修正系數(shù),一般取1.05；hf為沿程水頭損失；hj為局部水頭損失。

在流量、尾門位和水頭損失確定后，即可由式(4)算出河道斷面的各水力要素。

1.2 變維分形

按照B.B.MANDELBROT[20]的定義，分形分布滿足如下關(guān)系式:

s=Ay1-D

(5)

式中：S為歐氏長度；y為度量尺碼；D為分形維數(shù)；A為比例常數(shù)。

天然河道縱剖面分維，也存在統(tǒng)一的自相似性:

L=AQ1-D

(6)

式中：L為歐氏長度；Q為流量；D為分形維數(shù)；A為比例常數(shù)。

對(duì)式(5)兩端取自然對(duì)數(shù)得

lnL=lnA+(1-D)lnQ

(7)

把lnL繪于y軸，lnQ繪于x軸，用最小二乘法進(jìn)行直線擬合，它的斜率取為1-D，即可得出分維值D。

當(dāng)D為常數(shù)時(shí)，這種分形稱為常維分形。若D與特征線度呈函數(shù)關(guān)系，則稱之為變維分形。筆者采用由累計(jì)和序列獲得的分維模型，稱為“累計(jì)和分形”[21]。

1.3 計(jì)算河段的劃分

由于河流長度分維反應(yīng)了河流彎曲程度，而天然河流的不同分段的彎曲程度不盡相同，所以，筆者將河流分段，研究各分段與整體分維之間的關(guān)系。筆者以岷江、沱江匯流口為分界點(diǎn)，將水富(橫江匯流口)至重慶朱沱水文站的全長270km的長江干流河段分作3段進(jìn)行計(jì)算。第1段為水富(橫江匯流口)至宜賓(岷江匯流口)段，起點(diǎn)航道里程為1 073，終點(diǎn)航道里程為1 044。第2段為宜賓(岷江匯流口)至瀘州(沱江匯流口)段，起點(diǎn)航道里程為1 044，終點(diǎn)航道里程為912。第3段為瀘州(沱江匯流口)至朱沱水文站，起點(diǎn)航道里程為912，終點(diǎn)航道里程為805。筆者所述航道里程表示距宜昌距離，單位為km。為了保證計(jì)算結(jié)果的合理性與精度要求，研究河段布置了280個(gè)斷面，如圖2。金沙江江床地形圖為2008年3月測量；長江干流江床地形圖為2009年3月測量。

圖2 水富至朱沱段一維計(jì)算斷面布置Fig. 2 1-D arrangement of river section from Shuifu to Zhutuo

1.4 一維數(shù)學(xué)模型驗(yàn)證

1.4.1 水位驗(yàn)證

由于研究河段較長，限于篇幅，本次驗(yàn)證僅以有實(shí)測瞬時(shí)水面線資料的燒瓦沱至棧橋?yàn)┖佣?、桐梓林至黃蔥咀河段、宜賓至李莊河段為例，分別驗(yàn)證不同流量級(jí)下3個(gè)河段的水面線計(jì)算結(jié)果。燒瓦沱至棧橋?yàn)┖佣悟?yàn)證流量選取為洪水Q=10 322 m3/s、中水Q=8 660 m3/s和枯水Q=1 770 m3/s；桐梓林至黃蔥咀河段驗(yàn)證流量選取枯水Q=1 770 m3/s；宜賓至李莊河段選擇枯水與洪水兩級(jí)流量Q=2 960 、10 322 m3/s進(jìn)行驗(yàn)證。計(jì)算值與實(shí)測值繪于圖3。由圖3可知，模型與原型水位差值在±0.16 m以內(nèi)，誤差較小，計(jì)算結(jié)果與實(shí)測資料吻合度高。

(a)燒瓦沱至棧橋?yàn)┒嗡或?yàn)證

(b)桐梓林至黃蔥咀段水位驗(yàn)證

(c)宜賓至李莊段水位驗(yàn)證圖3 水位驗(yàn)證結(jié)果Fig.3 Water level verification resnlts

1.4.2 流速驗(yàn)證

限于篇幅，僅以有實(shí)測流速資料的馬皮包灘至和尚巖灘、棧橋?yàn)┲了槊诪┖佣巍⒁速e至李莊河段為例，分別驗(yàn)證不同流量級(jí)下3個(gè)河段的流速計(jì)算結(jié)果。燒瓦沱至棧橋?yàn)┖佣芜x取驗(yàn)證流量分別為洪水Q=10 322 m3/s、中水Q=8 680 m3/s和枯水Q=1 770 m3/s；棧橋?yàn)┲了槊诪┖佣芜x取洪、中、枯3級(jí)流量驗(yàn)證，驗(yàn)證流量分別為Q=10 322、8 680、1 770 m3/s；宜賓至李莊河段選取驗(yàn)證流量分別為枯水Q=2 960 m3/s和洪水Q=10 322 m3/s。計(jì)算值與實(shí)測值繪于圖4。

(a)燒瓦沱至棧橋?yàn)┒瘟魉衮?yàn)證

(b)棧橋?yàn)┲了槊诪┒瘟魉衮?yàn)證

(c)宜賓至李莊段流速驗(yàn)證圖4 流速驗(yàn)證結(jié)果Fig.4 Flow velocity verification results

由圖4可知，計(jì)算流速與實(shí)測流速的差值在±0.4 m/s以內(nèi)，誤差較小，計(jì)算結(jié)果與實(shí)測資料吻合度高。

1.5 糙率選取

河道糙率作為水力計(jì)算的主要參數(shù)，與河道糙率與河道水位、斷面形狀、坡降等參數(shù)相關(guān)，需要用實(shí)測資料進(jìn)行率定。長江上游河道糙率取值一般在0.03以上。根據(jù)數(shù)模計(jì)算結(jié)果，天然情況下研究河段糙率值變化很小，可采用線性內(nèi)插法進(jìn)行取值，將各分段天然狀態(tài)下部分流量級(jí)別下糙率值列于表1。

2 河道縱剖面分維數(shù)

依據(jù)第1節(jié)研究方法所述，以長江主干道實(shí)測地形圖為資料，統(tǒng)計(jì)特征流量下長江主干道的左右岸長度，結(jié)果見圖5。

從圖5可知，數(shù)據(jù)點(diǎn)分布趨勢不呈直線形式，說明研究河段長度呈現(xiàn)變維分形關(guān)系，需建立變維分形模型對(duì)其進(jìn)行處理。因此采用累計(jì)和變維分形的方法對(duì)長河段河流長度分形進(jìn)行計(jì)算，左岸計(jì)算過程見表2。右岸計(jì)算過程與左岸相同，計(jì)算結(jié)果見圖6及表3。

表2 水富至朱沱段左岸變維分維維數(shù)計(jì)算結(jié)果Table 2 Calculation results of variable fractal dimension of the left river bank of the river section from Shuifu to Zhutuo

圖6 水富至朱沱段2階累計(jì)和整體分維序列Fig.6 Second-order accumulation and integral fractal dimension sequence of the river section from Shuifu to Zhutuo

岸別直線相關(guān)方程分維值相關(guān)系數(shù)R2長江干流左岸2階y=1.1017x-6.96442.07380.9995長江干流右岸2階y=1.1057x-7.04272.07330.9995

注：y=lnQi，x=lnS2

由圖6可看出，經(jīng)2階累計(jì)和分形變換后，數(shù)據(jù)點(diǎn)與擬合直線吻合很好，說明水富至朱沱段河流分維呈現(xiàn)2階累計(jì)和分形關(guān)系。經(jīng)過計(jì)算，河段左岸分維數(shù)為D2=2.073 8，相關(guān)系數(shù)R2=0.999 5；右岸分維數(shù)D2=2.073 3，相關(guān)系數(shù)R2=0.999 5,均能通過F顯著性檢驗(yàn)。由此可見，水富至朱沱段河流長度具有二階分維特性。

3 分析與討論

筆者分別計(jì)算了研究河段的縱剖面分段分維數(shù)及整體分維數(shù)，將從分段分維數(shù)和整體分維數(shù)，縱剖面分維數(shù)與流量，縱剖面分維數(shù)與河道縱坡降3者的關(guān)系來探討河道縱剖面分形特性。

3.1 分段分維數(shù)

忽略較小支流的入?yún)R，以岷江、沱江匯流口作為分界點(diǎn)，將水富至朱沱270 km河段分作3段，即水富至宜賓段、宜賓至瀘州段、瀘州至朱沱段，按左、右岸分開計(jì)算各河段的分維數(shù)。計(jì)算結(jié)果見表4和圖7～圖9。

表4 分段計(jì)算數(shù)據(jù)直線相關(guān)方程及2階累計(jì)和分維值計(jì)算結(jié)果

注：y=lnQi，x=lnS2

圖7 水富至宜賓段2階累計(jì)和分段分維序列Fig.7 Second-order accumulation and subsection fractal dimension sequences of the river section from Shuifu to Yibin

圖8 宜賓至瀘州段2階累計(jì)和分段分維序列Fig.8 Second-order accumulation and subsection fractal dimension sequences of the river section from Yibin to Luzhou

圖9 瀘州至朱沱段2階累計(jì)和分段分維序列Fig.9 Second-order accumulation and subsection fractal dimension sequences of the river section from Luzhou to Zhutuo

由圖7～圖9可看出，水富至宜賓段左岸分維數(shù)為D2=2.078 2，相關(guān)系數(shù)R2=0.999 6，右岸分維數(shù)D2=2.077 9，相關(guān)系數(shù)R2=0.999 5；宜賓至瀘州段左岸分維數(shù)為D2=2.098 3，相關(guān)系數(shù)R2=0.999 2，右岸分維數(shù)D2=2.098 2，相關(guān)系數(shù)R2=0.999 5；瀘州至朱沱段左岸分維數(shù)為D2=2.101 8，相關(guān)系數(shù)R2=0.999 4，右岸分維數(shù)D2=2.101 6，相關(guān)系數(shù)R2=0.999 3，均能通過F顯著性檢驗(yàn)。說明數(shù)據(jù)點(diǎn)與擬合直線吻合地很好，即各分段具有較好的二階分維特征，得到的分形維數(shù)能夠很好地反映各分段河道的特征。

3.2 整體分維數(shù)與分段分維數(shù)的關(guān)系

表5為水富至朱沱段的整體和分段分維數(shù)結(jié)果。由表5可知，整體分維數(shù)小于各分段分維數(shù)。全河段左岸分維數(shù)為2.073 8，水富至宜賓段左岸分維數(shù)為2.078 2，比全河段左岸分維數(shù)增加了0.21%；宜賓至瀘州段左岸分維數(shù)為2.098 3，比全河段左岸分維數(shù)增加了1.18%；瀘州至朱沱段左岸分維數(shù)為2.101 8，比全河段左岸分維數(shù)增加了1.35%，由此可知分維數(shù)不斷增大，右岸亦如此，即整體分維數(shù)小于各分段分維數(shù)[22]。

表5 水富至朱沱段整體及分段分維數(shù)計(jì)算結(jié)果Table 5 Calculation results of integral and subsection fractal dimension of the river section from Shuifu to Zhutuo

由表5可知，整體分維數(shù)并不等于各河段分維數(shù)的算術(shù)平均值[23]。各河段左岸分維數(shù)算術(shù)平均值為2.092 8，并不等于全河段左岸分維數(shù)2.073 8；各河段右岸維數(shù)算術(shù)平均值為2.092 6，并不等于全河段右岸分維數(shù)2.073 3。所以全河段整體分維數(shù)不等于各河段分維數(shù)的算術(shù)平均值。

3.3 縱剖面分形維數(shù)與流量的關(guān)系

對(duì)研究河段3個(gè)分段的左、右岸縱剖面分維數(shù)進(jìn)行了計(jì)算，其相關(guān)系數(shù)都在0.99以上，說明數(shù)據(jù)點(diǎn)與擬合直線吻合得很好，流量與河道長度具有很好的相關(guān)性。將左、右岸Q=1 000～10 322 m3/s流量對(duì)應(yīng)的分維數(shù)繪于圖10。由圖10可知，縱剖面分形維數(shù)隨著流量不斷增大而不斷減小。

圖10 水富至朱沱段分形維數(shù)與流量關(guān)系Fig. 10 Relationship between fractal dimension and flow of the river section from Shuifu to Zhutuo

隨著河道來流量的加大，河岸逐漸被淹沒，河道彎曲程度逐漸變小[23]。因此，河道縱剖面分維數(shù)隨流量增大而減小。

3.4 縱剖面分形維數(shù)與河道縱坡降的關(guān)系

筆者統(tǒng)計(jì)了各河段在Q=1 000～10 322 m3/s流量下的河道縱坡降。圖11～圖13分別繪制了3個(gè)河段縱剖面分形維數(shù)與河道縱坡降關(guān)系曲線。由圖11～圖13可見，水富至宜賓段左岸，隨著流量增大，分維數(shù)減小，河道縱坡降不斷增大，分維數(shù)為2.099 7時(shí)坡降達(dá)到最大值0.249 7‰，對(duì)應(yīng)流量2 310 m3/s；分維數(shù)繼續(xù)減小，坡降開始減小，分維數(shù)為2.081 6時(shí)坡降達(dá)到最小值0.212 6‰，對(duì)應(yīng)流量3 969 m3/s；之后坡降隨分維數(shù)減小而增大。宜賓至瀘州段左岸分維數(shù)與河道坡降的變化規(guī)律與水富至宜賓段相類似。而瀘州至朱沱段左岸河道縱坡降與分形維數(shù)的變化曲線則與上述兩段不同，整體來說河道縱坡降隨著分維數(shù)減小而增大，分維數(shù)為2.123 5時(shí)坡降達(dá)到最大值0.260 7‰，對(duì)應(yīng)流量796 9 m3/s；之后坡降隨分維數(shù)減小而減小。各河段左、右岸縱坡降隨分形維數(shù)的變化規(guī)律基本相同，坡降隨分形維數(shù)的變化曲線類似波形。

圖11 水富至宜賓段分形維數(shù)與縱坡降關(guān)系Fig. 11 Relationship between fractal dimension and the gradient of the river section from Shuifu to Yibin

圖12 宜賓至瀘州段分形維數(shù)與縱坡降關(guān)系Fig. 12 Relationship between fractal dimension and the gradient of the river section from Yibin to Luzhou

圖13 瀘州至朱沱段分形維數(shù)與縱坡降關(guān)系Fig. 13 Relationship between fractal dimension and gradient of the river section from Luzhou to Zhutuo

水富至宜賓段位于向家壩下游，河床受樞紐下泄水流沖刷嚴(yán)重，河床床面形態(tài)變形加劇，床面高低不平，高程相差較大，且該河段地區(qū)地形高程相差較大，故坡降分維關(guān)系曲線呈現(xiàn)驟升驟降態(tài)勢,與簡諧波形狀相似。宜賓至瀘州段經(jīng)岷江、橫江匯流，流量較大，河床沖刷嚴(yán)重，坡降分維關(guān)系曲線亦驟升驟降。而瀘州至朱沱段地區(qū)地形較為平緩，水流較為平穩(wěn)，故坡降分維關(guān)系曲線較為平滑，與簡諧波形狀相似度最高。這進(jìn)一步揭示了流量、坡降、縱剖面分維數(shù)3者之間的關(guān)系，為城市防洪和航道整治等工程問題的解決提供一種新的依據(jù)。

4 結(jié) 論

1)筆者研究對(duì)象為向家壩下游約270 km的長河段，建立了一維非恒定-分形數(shù)學(xué)模型，該模型的計(jì)算結(jié)果與原型實(shí)測資料吻合良好，驗(yàn)證精度符合要求，表明該模型合理可靠。

2)研究河段左岸分維數(shù)為D2=2.073 8，相關(guān)系數(shù)R2=0.999 5；右岸分維數(shù)D2=2.073 3，相關(guān)系數(shù)R2=0.999 5，且均能通過F顯著性檢驗(yàn)。

3)水富至宜賓段左岸分維數(shù)為D2=2.078 2，相關(guān)系數(shù)R2=0.999 6，右岸分維數(shù)D2=2.077 9，相關(guān)系數(shù)R2=0.999 5；宜賓至瀘州段左岸分維數(shù)為D2=2.098 3，相關(guān)系數(shù)R2=0.999 2，右岸分維數(shù)D2=2.098 2，相關(guān)系數(shù)R2=0.999 5；瀘州至朱沱段左岸分維數(shù)為D2=2.101 8，相關(guān)系數(shù)R2=0.999 4，右岸分維數(shù)D2=2.101 6，相關(guān)系數(shù)R2=0.999 3，且均能通過F顯著性檢驗(yàn)。

4)研究河段縱剖面分形維數(shù)呈2階分維特性，且分段分維數(shù)大于整體分維數(shù)；各分段縱剖面分維數(shù)的算術(shù)平均值不等于整體分維數(shù)；縱剖面分維數(shù)隨流量增大而減小；各分段河道縱坡降隨縱剖面分形維數(shù)的變化曲線類似波形。

[1] 許模,王迪,漆繼紅,等.基于分形理論的卡斯特地貌形態(tài)分析[J].成都理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,38(3):328-333. XU Mu, WANG Di, Qi Jihong. Study on morphological characteristics of karst landform based on the fractal theory[J].JournalofChengduUniversityofTechnology(Science&TechnologyEdition, 2011, 38(3): 328-333.

[2] 范林峰,胡瑞林,張小艷,等.基于GIS和DEM的水系三維分形計(jì)盒維數(shù)的計(jì)算[J]. 地理與地理信息科學(xué),2012,28(6):28-30. FAN Linfeng, HU Ruilin, ZHANG Xiaoyan, et al. Calculation of 3D box dimension of river system based on GIS and DEM[J].GeographyandGeo-InformationScience, 2012, 28(6): 28-30.

[3] 倪志輝,吳立春,舒小紅.基于分形理論的挾沙水流摻混長度分析[J].人民黃河,2009, 31(9):32-33. NI Zhihui, WU Lichun, SHU Xiaohong. Analysis of water sediment mixing length based on fractal theory[J].YellowRiver, 2009, 31(9):32-33.

[4] 倪志輝,張緒進(jìn),胥潤生.長江黃河含沙量垂線分布的分形研究[J].人民長江,2011, 42(19):73-76. NI Zhihui, ZHANG Xujin, XU Runsheng. Fractal study on the vertical concentration distribution of sediment flow in Yangtze River and Yellow River[J].YangtzeRiver, 2011, 42(19):73-76.

[5] 江齊英,牛瑞卿.長江三峽庫岸幾何形態(tài)的分形分維研究[J].人民長江,2010,41(9):42-44. JIANG Qiying, NIU Ruiqing. Research on fractal and fractal dimension of Three Gorges Reservoir bank’s geometric shape[J].YangtzeRiver, 2010, 41(9): 42-44.

[6] 王衛(wèi)紅,徐鵬,田世民.分形理論在河型研究中的應(yīng)用探討[J].泥沙研究,2010(2):35-42. WANG Weihong, XU Peng, TIAN Shimin. Application of the fractal theory in the river pattern study[J].JournalofSedimentResearch, 2010(2):35-42.

[7] 高義,蘇奮振,周成虎,等.基于分形的中國大陸海岸線尺度效應(yīng)研究[J].地理學(xué)報(bào),2011,66(3):331-339. GAO Yi, SU Fenzhen, ZHOU Chenghu, et al. Scale effects of China mainland coastline based on fractal theory[J].ActaGeographicaSinica,2011, 66(3): 331-339.

[8] 周銀軍,陳立,歐陽娟,等.三峽蓄水后典型河段分形維數(shù)的變化分析[J].水科學(xué)進(jìn)展,2010,21(3):299-306. ZHOU Yinjun, CHEN Li, OUYANG Juan, et al. Changes in fractal dimensions of representative reaches after the impoundment of Three Gorges Project[J].AdvancesinWaterScience, 2010, 21(3): 299-306.

[9] 許光祥,鐘亮.河道剖面輪廓的一元分形插值模擬[J].人民長江,2012,43(3):20-23. XU Guangxiang, ZHONG Liang. Unary fractal interpolation simulation of channel profiles[J].YangtzeRiver, 2012, 43(3): 20-23.

[10] 郭碧云,傅旭東,張正峰.震后龍溪河流域河床縱剖面變化和溝谷分形[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2014,22(1):118-125. GUO Biyun, FU Xudong, ZHANG Zhengfeng. The Longitudinal change and gully bed fractality of Longxi river valley after the earthquake[J].JournalofBasicScienceandEngineering, 2014, 22(1):118-125.

[11] NIKORA V I, SAPOZHNIKOV V B, NOEVER D A. Fractal geometry of individual river channels and its computer simulation [J].WaterResourcesResearch, 1993, 29(29):3561-3568.

[12] SAPOZHNIKOV V, FOUFOULA-GEORGIOU E. Self-affinity in braided rivers[J].WaterResourcesResearch, 1996, 32(5): 1429-1439.

[13] FOUFOULA-GEORGIOU E, SAPOZHNIKOV V B. Anisotropic scaling in braided rivers: an integrated theoretical framework and results from application to an experimental river[J].WaterResourcesResearch, 1998, 34(4): 863-868.

[14] NYKANEN D K, FOUFOULA-GEORGIOU E, SAPOZHNIKOV V B. Study of spatial scaling in braided river patterns using synthetic aperture radar imagery[J].WaterResourcesResearch, 1998, 34(7):1795-1807.

[15] KINSNER W. A unified approach to fractal dimensions[J].InternationalJournalofCognitiveInformaticsandNaturalIntelligence, 2007, 1(4):26-46.

[16] 王欣凱,賈建軍. 浙江海島岸線分維數(shù)與其人工岸線比例的關(guān)系[J]. 海洋學(xué)研究, 2011, 29(4):25-31. WANG Xinkai, JIA Jianjun. Relationship between the fractal dimension and percentage of artificial coastlines of islands in Zhejiang Province[J].JournalofMarineSciences, 2011, 29(4): 25-31.

[17] ZHOU Z, YANG J, DENG Y, et al. Shortest-path fractal dimension for percolation in two and three dimensions[J].PhysicalReviewEStatisticalNonlinear&SoftMatterPhysics, 2012, 86(1): 1-5.

[18] NI Z H, WU L C, ZHANG X J, et al. A new model for the regulation width of waterway based on hydraulic geometry relation[J].MathematicalProblemsinEngineering, 2013(4):883-907.

[19] NI Z H, WU L C, WANG M H. The fractal dimension of river length based on the observed data[J].JournalofAppliedMathematics, 2013(1): 361-376.

[20] MANDELBROT B B. How long is the coast of britain? Statistical self-similarity and fractional dimension[J].Science, 1967, 156(3775):636.

[21] 吳立春,倪志輝.河流長度整體分維與局部分維關(guān)系研究——以長江重慶主城段為例[J].四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,36(1):147-151. WU Lichun, NI Zhihui. Study on the relationship between overall fractal dimension and partial fractal dimension of river length for example of the section of Chongqing City of Yangtze River[J].JournalofSichuanNormalUniversity(NaturalScience), 2013, 36(1): 147-151.

[22] 吳立春,倪志輝.長江重慶主城段河流長度分維數(shù)與洪水的關(guān)系[J].重慶交通大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,31(5):1042-1045. WU Lichun, NI Zhihui. Relationship between fractal dimension of the length in section of Chongqing city of Yangtze River and the flood[J].JournalofChongqingJiaotongUniversity(NaturalScience),2012, 31(5): 1042-1045.

[23] 倪志輝,吳立春,張緒進(jìn),等.潮流流速垂線分布的分類及其分形規(guī)律[J].人民黃河,2012,34(3):13-16. NI Zhihui, WU Lichun, ZHANG Xujin, et al. Patterns of vertical velocity profile of tidal current and their fractal law[J].YellowRiver, 2012, 34(3): 13-16.

(責(zé)任編輯：譚緒凱)

Fractal Characteristics of Longitudinal Section of River Bed

NI Zhihui1, 2, WANG Wei1, WU Lichun3

(1. Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education, National Engineering Research Center for Inland Waterway Regulation, Chongqing Jiaotong University, Chongqing 400074, P. R. China; 2. Southwestern Research Institute of Water Transportation Engineering, Chongqing Jiaotong University, Chongqing 400016, P. R. China; 3. Chongqing University of Education, Chongqing 400067, P. R. China)

Taking Shuifu to Zhutuo section in mainstream of the Yangtze River as the research object, an unsteady and fractal mathematics model with one-dimension was established. The studied section was divided into 3 segments, and the confluence of Minjiang River and Tuojiang River was taken as a demarcation. The fractal dimension of longitudinal section in river profile was calculated. The research indicates that: the fractal dimension of longitudinal section has two order fractal characteristics, and the overall fractal dimension is smaller than that of the piecewise. The value of overall fractal dimension is not equal to the arithmetic mean of fractal dimension of various sections. Fractal dimension of longitudinal section decreases as the flow increases. The longitudinal slope in various sections of river course is changing with the fractal dimension of longitudinal section, whose morphology is similar to the wave shape.

waterway engineering; longitudinal section; unsteady flow; fractal dimension; variable dimension fractal; numerical simulation

10.3969/j.issn.1674-0696.2017.07.10

2015-09-12；

2016-11-12

國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2016YFC0402104)；重慶市基礎(chǔ)與前沿研究計(jì)劃項(xiàng)目(cstc2016jcyjA0380)；省部共建水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室暨國家內(nèi)河航道整治工程技術(shù)研究中心開放基金項(xiàng)目(SLK2016B03)；內(nèi)河航道整治技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(NHHD-201514)

倪志輝(1980—)，男，湖南衡陽人，博士，副研究員，主要從事河流海岸水動(dòng)力學(xué)、環(huán)境及數(shù)值模擬方面的研究。E-mail：benny251@163.com。

TU612.3

A

1674-0696(2017)07-058-08