• <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ì)算
    videosex国产| 制服人妻中文乱码| 国产又爽黄色视频| 老司机在亚洲福利影院| 久久这里只有精品19| 色精品久久人妻99蜜桃| 国产精品爽爽va在线观看网站 | 欧美亚洲日本最大视频资源| 中国美女看黄片| 久久精品aⅴ一区二区三区四区| 少妇 在线观看| 亚洲国产日韩欧美精品在线观看 | 很黄的视频免费| 久久久久久久午夜电影| ponron亚洲| 亚洲欧美精品综合一区二区三区| 狠狠狠狠99中文字幕| 久久婷婷人人爽人人干人人爱| 日本五十路高清| 亚洲国产高清在线一区二区三 | 老司机在亚洲福利影院| 看片在线看免费视频| 欧美大码av| 欧美色欧美亚洲另类二区| 日韩大码丰满熟妇| 黄色视频不卡| 黄片小视频在线播放| 在线视频色国产色| 精品一区二区三区视频在线观看免费| 91老司机精品| 法律面前人人平等表现在哪些方面| 国产欧美日韩一区二区精品| 麻豆av在线久日| 亚洲七黄色美女视频| 国产免费av片在线观看野外av| 国产在线观看jvid| 日韩免费av在线播放| 日韩欧美国产在线观看| 久久婷婷人人爽人人干人人爱| 久久精品人妻少妇| 久久香蕉精品热| 国产精品乱码一区二三区的特点| netflix在线观看网站| 人人澡人人妻人| 高清在线国产一区| 亚洲天堂国产精品一区在线| 精品久久久久久久久久免费视频| 日日摸夜夜添夜夜添小说| 制服人妻中文乱码| 日韩大尺度精品在线看网址| 人人澡人人妻人| netflix在线观看网站| 日韩国内少妇激情av| 欧美绝顶高潮抽搐喷水| 欧美精品亚洲一区二区| 动漫黄色视频在线观看| 最好的美女福利视频网| 国产真人三级小视频在线观看| 国产私拍福利视频在线观看| 午夜福利欧美成人| 神马国产精品三级电影在线观看 | 午夜两性在线视频| 日韩欧美 国产精品| 亚洲精品一区av在线观看| 男人舔女人下体高潮全视频| 波多野结衣高清作品| 国产蜜桃级精品一区二区三区| 男人舔女人下体高潮全视频| АⅤ资源中文在线天堂| 极品教师在线免费播放| 最近最新中文字幕大全免费视频| 国产av又大| 久久久久国产一级毛片高清牌| 欧美日韩亚洲国产一区二区在线观看| 国内精品久久久久精免费| 熟妇人妻久久中文字幕3abv| 一区二区三区激情视频| 国产精品久久久久久亚洲av鲁大| 精品国产美女av久久久久小说| 国产精品,欧美在线| 久久久国产精品麻豆| 一个人观看的视频www高清免费观看 | 1024香蕉在线观看| 琪琪午夜伦伦电影理论片6080| 免费女性裸体啪啪无遮挡网站| 香蕉国产在线看| 亚洲免费av在线视频| 国产片内射在线| 在线观看舔阴道视频| 亚洲成人免费电影在线观看| 午夜福利在线在线| 草草在线视频免费看| 久久久水蜜桃国产精品网| 超碰成人久久| 国产精品乱码一区二三区的特点| 成人特级黄色片久久久久久久| 亚洲国产毛片av蜜桃av| 哪里可以看免费的av片| 69av精品久久久久久| 久久性视频一级片| 国产麻豆成人av免费视频| 精品午夜福利视频在线观看一区| 亚洲人成网站高清观看| 中文字幕高清在线视频| 久热爱精品视频在线9| 女人高潮潮喷娇喘18禁视频| 一区二区日韩欧美中文字幕| a级毛片a级免费在线| 亚洲欧美一区二区三区黑人| 一区二区三区精品91| 一a级毛片在线观看| 黄色视频,在线免费观看| 久久性视频一级片| 国产高清videossex| 少妇熟女aⅴ在线视频| 国产精品影院久久| 久久精品aⅴ一区二区三区四区| 哪里可以看免费的av片| 欧美精品啪啪一区二区三区| 91字幕亚洲| 一二三四社区在线视频社区8| 亚洲精华国产精华精| 午夜精品久久久久久毛片777| √禁漫天堂资源中文www| 亚洲第一电影网av| 色老头精品视频在线观看| 人人妻,人人澡人人爽秒播| 在线天堂中文资源库| 日韩高清综合在线| 熟女少妇亚洲综合色aaa.| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲aⅴ乱码一区二区在线播放 | 久久亚洲真实| 在线观看免费日韩欧美大片| 他把我摸到了高潮在线观看| 男女那种视频在线观看| 我的亚洲天堂| 国产精品98久久久久久宅男小说| 国产精品影院久久| 丁香欧美五月| 19禁男女啪啪无遮挡网站| 亚洲欧美一区二区三区黑人| 男人操女人黄网站| 成熟少妇高潮喷水视频| 天堂影院成人在线观看| 精品久久久久久久人妻蜜臀av| 亚洲五月色婷婷综合| 搡老岳熟女国产| 亚洲精品粉嫩美女一区| 免费搜索国产男女视频| 首页视频小说图片口味搜索| 欧美国产日韩亚洲一区| 超碰成人久久| 亚洲真实伦在线观看| 成人免费观看视频高清| tocl精华| 午夜久久久在线观看| 色综合婷婷激情| 日韩欧美 国产精品| 亚洲成人国产一区在线观看| 免费搜索国产男女视频| 国产亚洲欧美98| 50天的宝宝边吃奶边哭怎么回事| 在线十欧美十亚洲十日本专区| 久久欧美精品欧美久久欧美| 88av欧美| 在线观看午夜福利视频| 人妻久久中文字幕网| 一本综合久久免费| 人人妻,人人澡人人爽秒播| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩无卡精品| 神马国产精品三级电影在线观看 | 悠悠久久av| 国产一卡二卡三卡精品| 一级片免费观看大全| 久久香蕉激情| 国产av一区在线观看免费| 国产亚洲精品一区二区www| 国产高清视频在线播放一区| 久热这里只有精品99| 亚洲国产欧美一区二区综合| www.精华液| 欧美性长视频在线观看| 韩国精品一区二区三区| 男女床上黄色一级片免费看| 久久天躁狠狠躁夜夜2o2o| 熟女少妇亚洲综合色aaa.| 欧美性长视频在线观看| 久久久国产欧美日韩av| 欧美日韩中文字幕国产精品一区二区三区| 在线观看午夜福利视频| 色综合欧美亚洲国产小说| 老司机深夜福利视频在线观看| 999久久久国产精品视频| 国产精品二区激情视频| 18美女黄网站色大片免费观看| 成人亚洲精品av一区二区| 一本大道久久a久久精品| 午夜免费鲁丝| 亚洲成人免费电影在线观看| 亚洲,欧美精品.| 久久久久九九精品影院| 国产在线精品亚洲第一网站| 侵犯人妻中文字幕一二三四区| 亚洲片人在线观看| 18禁美女被吸乳视频| 精品午夜福利视频在线观看一区| 在线十欧美十亚洲十日本专区| 搡老岳熟女国产| 国产精品av久久久久免费| 非洲黑人性xxxx精品又粗又长| 此物有八面人人有两片| 黑丝袜美女国产一区| 精品熟女少妇八av免费久了| 看免费av毛片| 国内少妇人妻偷人精品xxx网站 | 国产成人精品久久二区二区免费| 国产99白浆流出| 欧美激情 高清一区二区三区| www.www免费av| 午夜免费鲁丝| 久久久久久九九精品二区国产 | 欧美激情极品国产一区二区三区| 在线观看日韩欧美| 国内精品久久久久精免费| 麻豆久久精品国产亚洲av| 亚洲av美国av| 俺也久久电影网| 久久热在线av| 国产精品免费视频内射| 国产不卡一卡二| 人成视频在线观看免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品一区二区www| 久久精品91蜜桃| 午夜福利在线在线| 免费高清在线观看日韩| 人成视频在线观看免费观看| 黄色片一级片一级黄色片| 韩国av一区二区三区四区| 在线视频色国产色| 久久人妻福利社区极品人妻图片| 桃色一区二区三区在线观看| 国产激情偷乱视频一区二区| 日本免费一区二区三区高清不卡| 亚洲精品美女久久久久99蜜臀| 1024手机看黄色片| 免费在线观看日本一区| 一区二区三区高清视频在线| 久久精品影院6| 男人操女人黄网站| 精品久久久久久久人妻蜜臀av| 亚洲真实伦在线观看| а√天堂www在线а√下载| 亚洲成国产人片在线观看| 中文字幕精品免费在线观看视频| www.自偷自拍.com| 国产黄a三级三级三级人| 琪琪午夜伦伦电影理论片6080| 国产麻豆成人av免费视频| 一级毛片高清免费大全| 国产精品综合久久久久久久免费| 亚洲中文av在线| 精品无人区乱码1区二区| 搞女人的毛片| 国产精品久久久av美女十八| 亚洲一码二码三码区别大吗| 精品午夜福利视频在线观看一区| 色婷婷久久久亚洲欧美| 日韩成人在线观看一区二区三区| 免费女性裸体啪啪无遮挡网站| 免费看日本二区| 久久久久国产精品人妻aⅴ院| 丰满的人妻完整版| 在线观看午夜福利视频| 夜夜躁狠狠躁天天躁| 色哟哟哟哟哟哟| 国产一区二区三区在线臀色熟女| 欧美色欧美亚洲另类二区| 超碰成人久久| 两个人免费观看高清视频| a在线观看视频网站| 黄色丝袜av网址大全| 精品少妇一区二区三区视频日本电影| 久久精品aⅴ一区二区三区四区| 黄片大片在线免费观看| 人妻丰满熟妇av一区二区三区| 日本免费一区二区三区高清不卡| 女性生殖器流出的白浆| 国产日本99.免费观看| 欧美性长视频在线观看| 他把我摸到了高潮在线观看| 国产熟女午夜一区二区三区| 看黄色毛片网站| 18禁黄网站禁片免费观看直播| 亚洲国产欧美一区二区综合| 国产精品一区二区免费欧美| 18禁观看日本| 欧美日韩乱码在线| 午夜成年电影在线免费观看| 不卡一级毛片| 中文字幕另类日韩欧美亚洲嫩草| 午夜久久久在线观看| 久久婷婷成人综合色麻豆| 女生性感内裤真人,穿戴方法视频| 亚洲成国产人片在线观看| 国产乱人伦免费视频| 哪里可以看免费的av片| 色综合欧美亚洲国产小说| 亚洲一区二区三区色噜噜| 欧美性猛交╳xxx乱大交人| 成熟少妇高潮喷水视频| 欧美三级亚洲精品| 俄罗斯特黄特色一大片| 午夜久久久在线观看| 国内少妇人妻偷人精品xxx网站 | 亚洲 欧美 日韩 在线 免费| 嫩草影院精品99| 午夜成年电影在线免费观看| 国内少妇人妻偷人精品xxx网站 | 成人国产一区最新在线观看| 99riav亚洲国产免费| 国产精品永久免费网站| 久久久久久人人人人人| 日本免费a在线| svipshipincom国产片| 一区二区三区高清视频在线| 性色av乱码一区二区三区2| 成人亚洲精品av一区二区| 亚洲成人久久性| 狠狠狠狠99中文字幕| 久久热在线av| ponron亚洲| 亚洲欧美一区二区三区黑人| 人人澡人人妻人| 亚洲国产毛片av蜜桃av| 一本久久中文字幕| 日韩国内少妇激情av| 正在播放国产对白刺激| 欧美日韩福利视频一区二区| 一区二区三区激情视频| 亚洲狠狠婷婷综合久久图片| 午夜久久久久精精品| 亚洲自拍偷在线| 中文字幕人妻熟女乱码| 欧美精品亚洲一区二区| 日韩视频一区二区在线观看| 怎么达到女性高潮| 两个人免费观看高清视频| 两人在一起打扑克的视频| 婷婷亚洲欧美| 国产av不卡久久| 亚洲免费av在线视频| 日本黄色视频三级网站网址| 国产成人精品久久二区二区91| 亚洲国产中文字幕在线视频| 丝袜在线中文字幕| 两个人视频免费观看高清| 91大片在线观看| 亚洲国产毛片av蜜桃av| 成年人黄色毛片网站| 手机成人av网站| 两个人免费观看高清视频| 精品日产1卡2卡| 午夜免费成人在线视频| 别揉我奶头~嗯~啊~动态视频| 麻豆成人av在线观看| 国产精品爽爽va在线观看网站 | 热re99久久国产66热| 哪里可以看免费的av片| 男女那种视频在线观看| 国产真人三级小视频在线观看| 久久香蕉激情| 黑人巨大精品欧美一区二区mp4| videosex国产| 欧美一级a爱片免费观看看 | 十八禁网站免费在线| 久久热在线av| 久久亚洲精品不卡| 美女 人体艺术 gogo| av免费在线观看网站| 黄频高清免费视频| 久久国产精品人妻蜜桃| 一级a爱视频在线免费观看| 久久精品国产亚洲av高清一级| 成人一区二区视频在线观看| av在线播放免费不卡| 久99久视频精品免费| 午夜激情av网站| 国产精品av久久久久免费| 91在线观看av| 深夜精品福利| 久久久久国产一级毛片高清牌| av天堂在线播放| 国产亚洲av嫩草精品影院| 亚洲自偷自拍图片 自拍| 欧美性猛交╳xxx乱大交人| 久久午夜综合久久蜜桃| 国产三级在线视频| 操出白浆在线播放| 在线观看66精品国产| 淫秽高清视频在线观看| 国产精品一区二区免费欧美| 老汉色∧v一级毛片| 一区二区三区精品91| 欧美乱色亚洲激情| 久久久久久久久久黄片| 精品欧美一区二区三区在线| 日本黄色视频三级网站网址| 18禁美女被吸乳视频| 中文在线观看免费www的网站 | 色综合站精品国产| 久久久久国产一级毛片高清牌| 天天躁夜夜躁狠狠躁躁| 午夜精品在线福利| 精品久久久久久久人妻蜜臀av| 国产午夜精品久久久久久| 欧美一级a爱片免费观看看 | av有码第一页| 黄色丝袜av网址大全| 成人特级黄色片久久久久久久| www.www免费av| 一级a爱视频在线免费观看| 日本a在线网址| 久久久久免费精品人妻一区二区 | 国内精品久久久久久久电影| 亚洲电影在线观看av| 香蕉丝袜av| 国产精品亚洲美女久久久| 黄色视频,在线免费观看| 欧美久久黑人一区二区| 日本 av在线| 国产视频内射| 欧美亚洲日本最大视频资源| 黄色视频不卡| 成人一区二区视频在线观看| 国产一区二区在线av高清观看| 看免费av毛片| 久久亚洲精品不卡| 国产成人av教育| av免费在线观看网站| 成在线人永久免费视频| 亚洲欧洲精品一区二区精品久久久| 国内少妇人妻偷人精品xxx网站 | 精品国内亚洲2022精品成人| 麻豆成人av在线观看| 久久久精品欧美日韩精品| 亚洲欧洲精品一区二区精品久久久| 嫩草影视91久久| 99riav亚洲国产免费| 香蕉久久夜色| 宅男免费午夜| 亚洲成a人片在线一区二区| 久久香蕉精品热| 亚洲人成77777在线视频| 黄色丝袜av网址大全| 午夜福利在线在线| 十八禁人妻一区二区| 国产精品亚洲美女久久久| 国产成年人精品一区二区| 精品久久蜜臀av无| 一进一出抽搐gif免费好疼| 欧美zozozo另类| 国产激情偷乱视频一区二区| 草草在线视频免费看| 免费在线观看成人毛片| 亚洲三区欧美一区| 亚洲欧美一区二区三区黑人| 欧美最黄视频在线播放免费| 国内揄拍国产精品人妻在线 | 欧美最黄视频在线播放免费| avwww免费| 三级毛片av免费| 国产黄片美女视频| 天天添夜夜摸| 91老司机精品| 国产不卡一卡二| 精品无人区乱码1区二区| 午夜两性在线视频| 亚洲精品中文字幕一二三四区| 脱女人内裤的视频| 国产蜜桃级精品一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产亚洲精品久久久久久毛片| 国产伦人伦偷精品视频| 啦啦啦观看免费观看视频高清| 99re在线观看精品视频| 一进一出抽搐动态| 亚洲国产高清在线一区二区三 | 99久久综合精品五月天人人| 精品电影一区二区在线| 日韩大码丰满熟妇| 欧美激情久久久久久爽电影| 亚洲人成77777在线视频| 50天的宝宝边吃奶边哭怎么回事| 午夜a级毛片| 女人被狂操c到高潮| 国内少妇人妻偷人精品xxx网站 | 国产精品综合久久久久久久免费| 中国美女看黄片| 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| 伦理电影免费视频| 精品国产亚洲在线| 黄频高清免费视频| 久久久久精品国产欧美久久久| 久久香蕉国产精品| 欧美久久黑人一区二区| 男女做爰动态图高潮gif福利片| 一卡2卡三卡四卡精品乱码亚洲| 国产三级在线视频| 亚洲成av人片免费观看| 国产午夜福利久久久久久| 香蕉丝袜av| 欧美日韩瑟瑟在线播放| 亚洲九九香蕉| 两个人看的免费小视频| 99精品欧美一区二区三区四区| 欧美最黄视频在线播放免费| 久久精品影院6| 久久国产精品人妻蜜桃| 午夜免费观看网址| 精品久久久久久久久久久久久 | 丰满的人妻完整版| 最新美女视频免费是黄的| 中文字幕人妻丝袜一区二区| 亚洲 欧美一区二区三区| 露出奶头的视频| 久久精品亚洲精品国产色婷小说| 69av精品久久久久久| 亚洲五月色婷婷综合| 欧美人与性动交α欧美精品济南到| 神马国产精品三级电影在线观看 | 高潮久久久久久久久久久不卡| 岛国视频午夜一区免费看| 成人一区二区视频在线观看| 欧美黑人欧美精品刺激| 国内久久婷婷六月综合欲色啪| 国产一区在线观看成人免费| 国产精品久久电影中文字幕| 巨乳人妻的诱惑在线观看| 国产成人欧美在线观看| 波多野结衣高清无吗| 日本撒尿小便嘘嘘汇集6| 久久精品人妻少妇| 欧美色视频一区免费| 亚洲 欧美 日韩 在线 免费| 欧美成人一区二区免费高清观看 | 一进一出抽搐动态| 久久人人精品亚洲av| 国产国语露脸激情在线看| 91成人精品电影| 免费在线观看亚洲国产| 精品国产超薄肉色丝袜足j| 国产成人精品久久二区二区91| netflix在线观看网站| 亚洲国产欧洲综合997久久, | 欧美中文综合在线视频| www.自偷自拍.com| 国产av又大| 精品不卡国产一区二区三区| 黄频高清免费视频| 中文字幕高清在线视频| 18禁美女被吸乳视频| 美女国产高潮福利片在线看| 国产成人影院久久av| 悠悠久久av| 欧美色欧美亚洲另类二区| 亚洲精品久久成人aⅴ小说| 国产高清videossex| 欧美激情极品国产一区二区三区| 波多野结衣巨乳人妻| 久久人妻福利社区极品人妻图片| 欧美激情久久久久久爽电影| 国产成人精品无人区| 女性被躁到高潮视频| 成人午夜高清在线视频 | 美女高潮到喷水免费观看| 欧美国产日韩亚洲一区| 精华霜和精华液先用哪个| 亚洲成a人片在线一区二区| 精品电影一区二区在线| 制服诱惑二区| 国产亚洲精品久久久久久毛片| 夜夜爽天天搞| 又黄又爽又免费观看的视频| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲精品综合一区在线观看 | 在线免费观看的www视频| 久久伊人香网站| 亚洲第一av免费看| 99久久精品国产亚洲精品| 亚洲 国产 在线| 在线观看日韩欧美| 999久久久精品免费观看国产| 久久精品成人免费网站| 曰老女人黄片| 两个人视频免费观看高清| 国产伦在线观看视频一区| 国产精品自产拍在线观看55亚洲| 看黄色毛片网站| a级毛片a级免费在线| netflix在线观看网站| 99国产精品一区二区蜜桃av| 国产乱人伦免费视频| 日韩高清综合在线| 国产99白浆流出| 日韩一卡2卡3卡4卡2021年| www日本黄色视频网| 国产精品久久电影中文字幕| 1024香蕉在线观看| 久久青草综合色| 欧美性猛交黑人性爽| 啦啦啦观看免费观看视频高清|