• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2017-08-01 00:00:55倪志輝吳立春
    關(guān)鍵詞:水富分維宜賓

    倪志輝,王 偉,吳立春

    (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

    猜你喜歡
    水富分維宜賓
    紹興市水富紡織器材有限公司
    木質(zhì)材料視覺感觀語義分維量化評(píng)價(jià)研究
    包裝工程(2023年16期)2023-08-25 11:36:32
    改進(jìn)的投影覆蓋方法對(duì)遼河河道粗糙床面分維量化研究
    The New Trends in Graphic Notation After 1945
    水富市總工會(huì):舉辦職工工間操培訓(xùn)班
    水富市總工會(huì):舉辦首屆職工籃球比賽
    宜賓面塑的保護(hù)、傳承與創(chuàng)新探討
    宜賓(外四首)
    岷峨詩稿(2017年4期)2017-04-20 06:26:33
    宜賓豆腐乳毛霉分離及應(yīng)用
    基于分形滲流模型的導(dǎo)電瀝青混凝土的分維計(jì)算
    久久中文看片网| 国产精品永久免费网站| 免费观看人在逋| 色5月婷婷丁香| 久久精品国产亚洲av天美| 中亚洲国语对白在线视频| 免费观看人在逋| 国产精品一区二区三区四区久久| 大型黄色视频在线免费观看| 老熟妇乱子伦视频在线观看| 女同久久另类99精品国产91| 天天一区二区日本电影三级| 欧美日韩国产亚洲二区| 亚洲精品久久国产高清桃花| 日韩欧美精品免费久久| 99热只有精品国产| 精品午夜福利在线看| 男人舔女人下体高潮全视频| 97碰自拍视频| 国产精品一区二区性色av| 在线国产一区二区在线| 欧美另类亚洲清纯唯美| 老司机午夜福利在线观看视频| 日本黄色片子视频| 亚洲国产精品sss在线观看| 永久网站在线| 欧美日本视频| 亚洲不卡免费看| 99在线视频只有这里精品首页| 欧美一区二区亚洲| 成人三级黄色视频| 人妻丰满熟妇av一区二区三区| 国产精品永久免费网站| 欧美成人a在线观看| 99九九线精品视频在线观看视频| 天堂av国产一区二区熟女人妻| 午夜爱爱视频在线播放| 精品人妻偷拍中文字幕| 国产乱人伦免费视频| 97热精品久久久久久| а√天堂www在线а√下载| 国产精品日韩av在线免费观看| 长腿黑丝高跟| 不卡一级毛片| 亚洲性久久影院| 最近最新中文字幕大全电影3| 国产精品久久久久久久久免| 久久人人精品亚洲av| 欧美三级亚洲精品| 18+在线观看网站| 欧美潮喷喷水| 国产亚洲精品久久久com| 久久精品综合一区二区三区| 国产 一区精品| 波多野结衣巨乳人妻| 免费黄网站久久成人精品| 久久人人爽人人爽人人片va| 午夜影院日韩av| 日韩欧美 国产精品| netflix在线观看网站| 日韩大尺度精品在线看网址| 1024手机看黄色片| 男人的好看免费观看在线视频| 亚州av有码| 九九久久精品国产亚洲av麻豆| 12—13女人毛片做爰片一| 网址你懂的国产日韩在线| 在线播放无遮挡| 99久国产av精品| 一本一本综合久久| 久久精品综合一区二区三区| 精品乱码久久久久久99久播| 网址你懂的国产日韩在线| 国产亚洲精品av在线| 最近最新中文字幕大全电影3| 九九爱精品视频在线观看| 免费看av在线观看网站| 老司机福利观看| 中文字幕免费在线视频6| 欧美zozozo另类| 岛国在线免费视频观看| 久久欧美精品欧美久久欧美| 韩国av一区二区三区四区| 婷婷丁香在线五月| 91麻豆av在线| 久久精品国产亚洲网站| 免费观看的影片在线观看| 99久久中文字幕三级久久日本| 最好的美女福利视频网| 久久国产乱子免费精品| 99久久久亚洲精品蜜臀av| 免费在线观看影片大全网站| 国内精品一区二区在线观看| 免费无遮挡裸体视频| 国内精品一区二区在线观看| 久久久久免费精品人妻一区二区| 国产伦精品一区二区三区四那| 久久人妻av系列| 啦啦啦韩国在线观看视频| 狂野欧美白嫩少妇大欣赏| 国产成人av教育| 国产精品精品国产色婷婷| 亚洲七黄色美女视频| a级毛片a级免费在线| 在线免费观看的www视频| 亚洲综合色惰| 久久热精品热| 婷婷精品国产亚洲av| 久久香蕉精品热| 欧美激情久久久久久爽电影| 不卡一级毛片| 两个人的视频大全免费| 日韩中字成人| 一区福利在线观看| 国产91精品成人一区二区三区| 在线免费观看的www视频| 极品教师在线免费播放| 91久久精品国产一区二区成人| 日本爱情动作片www.在线观看 | 一级黄色大片毛片| 国产老妇女一区| 欧美+亚洲+日韩+国产| 男人舔奶头视频| 丰满的人妻完整版| 亚洲成人精品中文字幕电影| АⅤ资源中文在线天堂| 嫩草影院新地址| 丰满乱子伦码专区| 欧美黑人巨大hd| 免费人成视频x8x8入口观看| 嫁个100分男人电影在线观看| 免费看美女性在线毛片视频| 男女之事视频高清在线观看| а√天堂www在线а√下载| 免费不卡的大黄色大毛片视频在线观看 | 熟妇人妻久久中文字幕3abv| 不卡视频在线观看欧美| 成人鲁丝片一二三区免费| 久久久久久久亚洲中文字幕| 国产精品久久视频播放| 99久久无色码亚洲精品果冻| 日本黄色片子视频| 欧美成人性av电影在线观看| 久久久久久久精品吃奶| 女人十人毛片免费观看3o分钟| 久久久久久久午夜电影| 极品教师在线免费播放| 久久亚洲真实| 日韩 亚洲 欧美在线| 狠狠狠狠99中文字幕| 亚洲国产色片| 如何舔出高潮| 淫妇啪啪啪对白视频| 国产成人av教育| 桃色一区二区三区在线观看| a在线观看视频网站| 国产精品美女特级片免费视频播放器| 给我免费播放毛片高清在线观看| www.www免费av| 亚洲国产色片| 国产真实乱freesex| 精品99又大又爽又粗少妇毛片 | 午夜福利欧美成人| 在线免费观看的www视频| 我的女老师完整版在线观看| av中文乱码字幕在线| 亚洲综合色惰| 精品午夜福利在线看| 一个人免费在线观看电影| 色播亚洲综合网| 真人一进一出gif抽搐免费| 最后的刺客免费高清国语| 国产69精品久久久久777片| 三级毛片av免费| 亚洲黑人精品在线| 成人性生交大片免费视频hd| 很黄的视频免费| 好男人在线观看高清免费视频| 免费搜索国产男女视频| 天美传媒精品一区二区| www.www免费av| 日本爱情动作片www.在线观看 | 亚洲国产欧洲综合997久久,| 91午夜精品亚洲一区二区三区 | 欧美色视频一区免费| av女优亚洲男人天堂| 成人性生交大片免费视频hd| bbb黄色大片| 最近最新中文字幕大全电影3| 精品一区二区三区av网在线观看| 天天一区二区日本电影三级| 欧美国产日韩亚洲一区| av在线亚洲专区| 91久久精品电影网| 亚洲成人中文字幕在线播放| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品久久久com| 啦啦啦啦在线视频资源| 久久久久久久亚洲中文字幕| 人妻制服诱惑在线中文字幕| 欧美激情久久久久久爽电影| 亚洲国产高清在线一区二区三| 国产私拍福利视频在线观看| 国产高清视频在线观看网站| 99国产精品一区二区蜜桃av| 别揉我奶头~嗯~啊~动态视频| 国产大屁股一区二区在线视频| 欧美中文日本在线观看视频| 亚洲第一电影网av| 欧美最新免费一区二区三区| 久久精品久久久久久噜噜老黄 | 日本a在线网址| 欧美成人性av电影在线观看| 成人鲁丝片一二三区免费| 直男gayav资源| 亚洲经典国产精华液单| 欧美+亚洲+日韩+国产| 日韩欧美 国产精品| 国产精品99久久久久久久久| 国产男靠女视频免费网站| 国产高清激情床上av| 国产一区二区三区在线臀色熟女| 免费人成视频x8x8入口观看| 淫秽高清视频在线观看| videossex国产| 亚洲精华国产精华精| 国产成人影院久久av| 亚洲国产精品久久男人天堂| 直男gayav资源| av专区在线播放| 久久久久国产精品人妻aⅴ院| 亚洲国产精品sss在线观看| 哪里可以看免费的av片| 国产男靠女视频免费网站| 听说在线观看完整版免费高清| 亚洲第一区二区三区不卡| 黄色女人牲交| 婷婷亚洲欧美| 国产精品国产高清国产av| 成人综合一区亚洲| 国产男靠女视频免费网站| 欧美精品国产亚洲| 欧美黑人巨大hd| 国产精品1区2区在线观看.| 一区二区三区四区激情视频 | 国产精品免费一区二区三区在线| 嫩草影视91久久| 免费av不卡在线播放| 啦啦啦韩国在线观看视频| 亚洲av五月六月丁香网| 97超级碰碰碰精品色视频在线观看| 午夜激情福利司机影院| 三级毛片av免费| 91麻豆av在线| 国产精品久久电影中文字幕| 十八禁国产超污无遮挡网站| 欧美高清成人免费视频www| 日韩在线高清观看一区二区三区 | 永久网站在线| 亚洲成人免费电影在线观看| 99久久久亚洲精品蜜臀av| 99久久中文字幕三级久久日本| 欧美一区二区亚洲| 亚洲精品影视一区二区三区av| 精品日产1卡2卡| 免费在线观看日本一区| 少妇的逼水好多| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产自在天天线| 亚洲av成人av| 欧美成人一区二区免费高清观看| 亚洲美女搞黄在线观看 | 久久精品国产鲁丝片午夜精品 | 亚洲成人免费电影在线观看| 亚洲男人的天堂狠狠| 日日摸夜夜添夜夜添av毛片 | 人人妻人人看人人澡| 91在线观看av| 国产伦在线观看视频一区| 看十八女毛片水多多多| 欧美成人一区二区免费高清观看| 一本一本综合久久| 91精品国产九色| 欧美成人一区二区免费高清观看| 国产精品av视频在线免费观看| 99久久久亚洲精品蜜臀av| 日韩国内少妇激情av| 久久久久精品国产欧美久久久| 人妻久久中文字幕网| 免费一级毛片在线播放高清视频| 亚洲avbb在线观看| 搞女人的毛片| 国产伦在线观看视频一区| 国产爱豆传媒在线观看| 亚洲综合色惰| 制服丝袜大香蕉在线| 日韩 亚洲 欧美在线| 亚洲一区高清亚洲精品| a级一级毛片免费在线观看| 久久热精品热| 亚洲 国产 在线| 国产91精品成人一区二区三区| 亚洲va在线va天堂va国产| 美女 人体艺术 gogo| 直男gayav资源| 三级国产精品欧美在线观看| 最好的美女福利视频网| 九九久久精品国产亚洲av麻豆| 最近视频中文字幕2019在线8| 美女大奶头视频| 国产国拍精品亚洲av在线观看| 亚洲熟妇熟女久久| 美女cb高潮喷水在线观看| 午夜精品在线福利| 能在线免费观看的黄片| 日本在线视频免费播放| 亚洲精华国产精华液的使用体验 | 深爱激情五月婷婷| 99在线人妻在线中文字幕| 久久久国产成人免费| 免费看日本二区| 老熟妇乱子伦视频在线观看| 一区二区三区高清视频在线| 免费电影在线观看免费观看| 听说在线观看完整版免费高清| 毛片女人毛片| 十八禁网站免费在线| 精品免费久久久久久久清纯| 男女边吃奶边做爰视频| 嫁个100分男人电影在线观看| 国产女主播在线喷水免费视频网站 | 岛国在线免费视频观看| 成人鲁丝片一二三区免费| 久久久久久久久久久丰满 | 极品教师在线视频| 国产亚洲欧美98| 国产麻豆成人av免费视频| 韩国av在线不卡| 成人二区视频| 人妻丰满熟妇av一区二区三区| 午夜免费男女啪啪视频观看 | 精品久久久久久久久av| 亚洲美女视频黄频| 精品一区二区免费观看| 有码 亚洲区| 亚洲欧美精品综合久久99| 老司机福利观看| 亚洲天堂国产精品一区在线| 在线播放国产精品三级| 十八禁国产超污无遮挡网站| 身体一侧抽搐| 国产免费av片在线观看野外av| 性色avwww在线观看| 两个人视频免费观看高清| 又爽又黄无遮挡网站| netflix在线观看网站| 亚洲在线观看片| 亚洲色图av天堂| 亚洲最大成人中文| 日日啪夜夜撸| 久久人人精品亚洲av| 国产精品爽爽va在线观看网站| 1024手机看黄色片| 日韩大尺度精品在线看网址| 色哟哟哟哟哟哟| 国产精品乱码一区二三区的特点| 最好的美女福利视频网| 男人舔奶头视频| 我要看日韩黄色一级片| 精品午夜福利视频在线观看一区| 欧美激情国产日韩精品一区| 1000部很黄的大片| 国产亚洲91精品色在线| av女优亚洲男人天堂| 可以在线观看毛片的网站| 91久久精品电影网| 日日撸夜夜添| 免费观看的影片在线观看| 真实男女啪啪啪动态图| 在线天堂最新版资源| 精品久久久久久久久av| 99精品久久久久人妻精品| 在线播放无遮挡| 久久99热这里只有精品18| 嫩草影院入口| 在现免费观看毛片| 日日夜夜操网爽| av天堂中文字幕网| 91午夜精品亚洲一区二区三区 | 桃色一区二区三区在线观看| 十八禁国产超污无遮挡网站| 国产老妇女一区| 亚州av有码| 国产中年淑女户外野战色| 成年版毛片免费区| 国产高清激情床上av| 国产伦人伦偷精品视频| 女的被弄到高潮叫床怎么办 | 国产熟女欧美一区二区| 久久久久久九九精品二区国产| 成人特级av手机在线观看| 性欧美人与动物交配| 大又大粗又爽又黄少妇毛片口| .国产精品久久| 免费观看在线日韩| 又紧又爽又黄一区二区| 天堂动漫精品| 成熟少妇高潮喷水视频| 网址你懂的国产日韩在线| 国产主播在线观看一区二区| 亚洲欧美日韩高清专用| 久久久久久久久中文| 老熟妇乱子伦视频在线观看| 国产黄a三级三级三级人| 极品教师在线免费播放| 亚洲av不卡在线观看| 伦理电影大哥的女人| 亚州av有码| 99热这里只有是精品在线观看| 国产精品亚洲一级av第二区| 亚州av有码| 久99久视频精品免费| 又紧又爽又黄一区二区| 欧美黑人欧美精品刺激| 99在线视频只有这里精品首页| 性欧美人与动物交配| 免费无遮挡裸体视频| av国产免费在线观看| 床上黄色一级片| 好男人在线观看高清免费视频| 国产精品国产高清国产av| 免费看a级黄色片| 日本a在线网址| 国产亚洲欧美98| 亚洲专区国产一区二区| 亚洲精华国产精华液的使用体验 | 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| 免费看a级黄色片| 亚洲综合色惰| 别揉我奶头~嗯~啊~动态视频| 亚洲国产欧美人成| 在线观看美女被高潮喷水网站| 午夜精品久久久久久毛片777| 一本精品99久久精品77| 午夜影院日韩av| 两个人视频免费观看高清| 一进一出抽搐动态| 一个人看的www免费观看视频| 国产精品三级大全| 九色国产91popny在线| 最好的美女福利视频网| 久久久色成人| 日韩国内少妇激情av| 国产精品人妻久久久影院| 免费看日本二区| 中文字幕av成人在线电影| 欧美激情在线99| 亚洲人成伊人成综合网2020| 日韩精品中文字幕看吧| 日日撸夜夜添| 欧洲精品卡2卡3卡4卡5卡区| 日韩一本色道免费dvd| 俺也久久电影网| 色综合婷婷激情| 免费av毛片视频| 好男人在线观看高清免费视频| 一个人看视频在线观看www免费| 久久精品综合一区二区三区| 中出人妻视频一区二区| 国产精品98久久久久久宅男小说| 搡老妇女老女人老熟妇| 精品久久久久久久久久免费视频| 久久精品夜夜夜夜夜久久蜜豆| 人妻少妇偷人精品九色| 免费观看在线日韩| 九九热线精品视视频播放| 精品午夜福利视频在线观看一区| 午夜免费成人在线视频| 成人无遮挡网站| 天天一区二区日本电影三级| 久久精品影院6| 最近最新中文字幕大全电影3| 日日啪夜夜撸| 国产美女午夜福利| 国产午夜精品论理片| 国产av不卡久久| 国产高清三级在线| 国产一区二区三区av在线 | 亚洲av中文字字幕乱码综合| 国产私拍福利视频在线观看| 国产美女午夜福利| 精品一区二区三区人妻视频| 午夜亚洲福利在线播放| 免费在线观看日本一区| 久久亚洲精品不卡| 又爽又黄无遮挡网站| 在线观看免费视频日本深夜| 美女高潮喷水抽搐中文字幕| 精品日产1卡2卡| 男人舔奶头视频| 99riav亚洲国产免费| av专区在线播放| 国产精品综合久久久久久久免费| 午夜福利视频1000在线观看| 亚洲无线观看免费| 高清日韩中文字幕在线| 欧美中文日本在线观看视频| 国产高潮美女av| 美女xxoo啪啪120秒动态图| 国产精品久久久久久久久免| 看免费成人av毛片| 国产av不卡久久| 亚洲av成人精品一区久久| 日韩精品有码人妻一区| 日本撒尿小便嘘嘘汇集6| 内射极品少妇av片p| 色综合亚洲欧美另类图片| 中文在线观看免费www的网站| 国产不卡一卡二| 国产欧美日韩一区二区精品| 国产精品一区二区性色av| 少妇猛男粗大的猛烈进出视频 | 久久精品国产亚洲网站| 三级男女做爰猛烈吃奶摸视频| 少妇裸体淫交视频免费看高清| 在线观看舔阴道视频| 国产男人的电影天堂91| 亚洲欧美日韩高清在线视频| 美女免费视频网站| 日本黄色片子视频| 欧美三级亚洲精品| 给我免费播放毛片高清在线观看| 精品福利观看| 久久久久久久精品吃奶| 又爽又黄a免费视频| 九九久久精品国产亚洲av麻豆| 日本黄色片子视频| 日本五十路高清| 亚洲无线在线观看| 99久国产av精品| 国产爱豆传媒在线观看| 国产精品美女特级片免费视频播放器| 波野结衣二区三区在线| 99在线人妻在线中文字幕| 国产激情偷乱视频一区二区| 亚洲综合色惰| 精品久久久噜噜| 少妇的逼好多水| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产三级国产av玫瑰| 国产精品久久电影中文字幕| 国产伦精品一区二区三区视频9| 蜜桃亚洲精品一区二区三区| 婷婷色综合大香蕉| 韩国av在线不卡| 日韩欧美精品v在线| 亚洲欧美日韩高清专用| 久久热精品热| 国产欧美日韩一区二区精品| 中文亚洲av片在线观看爽| av视频在线观看入口| 免费av不卡在线播放| 欧美性猛交╳xxx乱大交人| 伦精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 18+在线观看网站| 精品午夜福利在线看| 亚洲av中文字字幕乱码综合| 最近视频中文字幕2019在线8| 国产精品野战在线观看| 国产伦精品一区二区三区视频9| 乱码一卡2卡4卡精品| 婷婷六月久久综合丁香| 国产综合懂色| 国产主播在线观看一区二区| 国产在线精品亚洲第一网站| 哪里可以看免费的av片| 国产成人a区在线观看| 久久精品影院6| 亚洲三级黄色毛片| 日韩强制内射视频| 欧美极品一区二区三区四区| 久久久成人免费电影| 性插视频无遮挡在线免费观看| 国产黄片美女视频| 色精品久久人妻99蜜桃| 日韩精品中文字幕看吧| 啪啪无遮挡十八禁网站| 国产色爽女视频免费观看| 久久精品国产清高在天天线| 国产又黄又爽又无遮挡在线| 久久99热这里只有精品18| 国产精品久久电影中文字幕| 国产免费一级a男人的天堂| 18禁黄网站禁片午夜丰满| av福利片在线观看| 日韩亚洲欧美综合| 黄色日韩在线| 免费观看的影片在线观看| 婷婷精品国产亚洲av| 俺也久久电影网| 欧美最黄视频在线播放免费| 日韩强制内射视频| 亚洲精华国产精华液的使用体验 | videossex国产| 网址你懂的国产日韩在线| 最近在线观看免费完整版| 女同久久另类99精品国产91| 国产真实乱freesex| 亚洲人成网站高清观看| 99久久中文字幕三级久久日本| 午夜久久久久精精品| 高清日韩中文字幕在线| 久久精品国产亚洲av天美|