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

    適用于中小河流受損岸坡修復(fù)的楔形丁壩群水力特性初探

    2018-09-18 07:08:36,,
    長江科學(xué)院院報(bào) 2018年9期
    關(guān)鍵詞:凹岸丁壩楔形

    ,,

    (1.交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究院,天津 300456;2.長江科學(xué)院 水力學(xué)研究所,武漢 430010)

    1 研究背景

    近年來,隨著河流治理工程理念的更新和變化,在滿足水資源開發(fā)利用和防洪安全的前提下,如何兼顧河流利用對(duì)生態(tài)環(huán)境的影響,促進(jìn)人與自然的和諧相處,實(shí)現(xiàn)河流利用的長期可持續(xù)性發(fā)展,已然成為現(xiàn)代水利科學(xué)關(guān)注的熱點(diǎn)。水利工程的生態(tài)影響問題已引起了社會(huì)各界的廣泛關(guān)注,同時(shí)也逐漸成為制約水利水電建設(shè)發(fā)展新的瓶頸[1]。鑒于此,河流再自然化日益成為人們新的必然的選擇[2-3]。

    新型的河道內(nèi)流動(dòng)控制結(jié)構(gòu)物,不同于傳統(tǒng)護(hù)岸單純以增強(qiáng)岸坡防護(hù)剛度為重點(diǎn),是通過主動(dòng)調(diào)整近岸處的流動(dòng)形態(tài),將高流速水流挑離至遠(yuǎn)離河岸位置,以達(dá)到保護(hù)堤岸免受沖蝕的目的。此外,經(jīng)河道內(nèi)流動(dòng)控制結(jié)構(gòu)物作用下流場(chǎng)調(diào)整所形成的沖坑還具有附帶的生態(tài)效益[4]。流動(dòng)控制結(jié)構(gòu)物作用下河道內(nèi)將形成新的流動(dòng)形態(tài)和沖淤格局,構(gòu)造出急流-緩流交替、深潭-淺灘交錯(cuò)的水流條件和多樣性地貌,使其更貼近天然河流的形態(tài)特征,為魚類或其他底棲生物提供適宜的棲息和產(chǎn)卵場(chǎng)地,從而有助于河流生物種群多樣性的提高。

    國內(nèi)外自20世紀(jì)90年代以來,已逐步開展河道內(nèi)新型流動(dòng)控制結(jié)構(gòu)物的模型試驗(yàn)和數(shù)值模擬研究工作,如導(dǎo)流屏[5-8]、堰壩[9-12]等。而楔形丁壩群以其阻水性小、穩(wěn)定性好和節(jié)省材料的優(yōu)點(diǎn)在中小河流岸坡修復(fù)領(lǐng)域正獲得廣泛的應(yīng)用[13-17],并初步形成了具有實(shí)踐指導(dǎo)意義的工程設(shè)計(jì)準(zhǔn)則,而國內(nèi)對(duì)于楔形丁壩目前還鮮見有相關(guān)的研究報(bào)道。

    2 河道內(nèi)流動(dòng)控制結(jié)構(gòu)物

    目前,國外在傳統(tǒng)丁壩形式基礎(chǔ)上衍生出了多種河道內(nèi)流動(dòng)控制結(jié)構(gòu)物類型,已廣泛用于河道生態(tài)修復(fù)領(lǐng)域,分類詳細(xì),多數(shù)已形成可用于指導(dǎo)工程實(shí)踐的設(shè)計(jì)指南,并積累了較為豐富的安裝、監(jiān)測(cè)和維護(hù)等工程經(jīng)驗(yàn)。而國內(nèi)對(duì)其劃分不甚明確,且多數(shù)還停留在河道控導(dǎo)工程等治河或通航領(lǐng)域,主要以工程應(yīng)用為導(dǎo)向,對(duì)其分類命名較為籠統(tǒng),且對(duì)生態(tài)效益的關(guān)注較少。

    為與傳統(tǒng)丁壩形式相區(qū)分,選取幾種常見的河道內(nèi)流動(dòng)控制結(jié)構(gòu)物(見圖1),就其各自的特點(diǎn)和適用性進(jìn)行介紹。

    圖1 常見的河道內(nèi)流動(dòng)控制結(jié)構(gòu)物類型Fig.1 Illustrations of typical instream flow control structures

    2.1 單一型堰(rock vane)

    單一型堰是一種獨(dú)臂導(dǎo)流結(jié)構(gòu)物,從河岸伸出突入河道當(dāng)中,突出部分約占河道寬度的1/3,并與來流構(gòu)成一定的角度,突入河道的一端直接插入河床,從而自河岸至另一端形成一定坡度,致使突入河道一端即使在小流量時(shí)依然被淹沒在水下,而其余部分則裸露在外。單一型堰通過將高流速水流挑離河岸,在河岸附近形成靜水區(qū)域,從而減弱水流對(duì)河岸的侵蝕,以達(dá)到保護(hù)河岸的目的。

    2.2 J型堰(J-hook vane)

    J型堰是對(duì)單一型堰的一種改進(jìn),在突入水中的一端額外布置留有間隔的漂石塊,從而在平面上呈“鉤”狀分布,“鉤”狀部分約占整個(gè)河道寬度的1/3?!般^”狀部分漂石間隔消散了來流部分能量,在下游形成漩渦,有利于泥沙輸移并為魚類提供掩蔽,相比單一型堰能形成更加寬大和更深的沖坑,從而更加利于底棲生物棲息。J型堰同單一型堰一樣,主要布置在河流彎道的凹岸一側(cè),通過降低河道坡降、流速和切應(yīng)力來防止河岸侵蝕。

    2.3 交叉堰(cross vane)

    交叉堰是將左右兩岸2個(gè)獨(dú)臂的單一型堰通過與來流相垂直的漂石體相連,三者各占河寬的1/3,垂直來流的中間堰體即使在小流量情況下依然被淹沒在水下。交叉堰主要用來控制坡降并將高流速水流挑離河岸,降低河床切應(yīng)力,從而減弱對(duì)河岸的侵蝕,形成穩(wěn)定的寬深比。交叉堰下游所形成的沖坑,在堰體結(jié)構(gòu)的掩蔽之下,能為底棲生物提供良好的棲息和產(chǎn)卵場(chǎng)地。

    2.4 W型堰(W-weir)

    W型堰結(jié)構(gòu)同交叉堰類似,同樣是一種跨河道的結(jié)構(gòu)物,通過兩列“V”型漂石體將左右兩岸2個(gè)單一型堰相連。從上游朝下游方向看去,該堰體外型呈“W”型,與交叉堰類似,該結(jié)構(gòu)物同樣起到控制坡降和挑離高流速水流作用以達(dá)到保護(hù)河岸的目的。與交叉堰不同的是在堰體下游構(gòu)成2條平行的河道深泓線,形成多樣化的流態(tài)特征,對(duì)改善魚類棲息地具有積極意義。

    2.5 導(dǎo)流屏(submerged vanes/Iowa vanes)

    導(dǎo)流屏是一種小型導(dǎo)流結(jié)構(gòu)物,垂直安裝在河床上,與水流構(gòu)成一定的角度,通過人工環(huán)流來重新分布水流與河道斷面內(nèi)的泥沙運(yùn)動(dòng),從而控制河床的沖淤形態(tài)。導(dǎo)流屏一方面激發(fā)次生環(huán)流,用以消除河流彎段離心力引起的橫向環(huán)流,改變近河床的切應(yīng)力分布,將含沙量較大的水流導(dǎo)入導(dǎo)流屏區(qū),促進(jìn)導(dǎo)流屏控制區(qū)的淤積;另一方面,它可增大局部水流阻力,減小導(dǎo)流屏控制區(qū)的流速和輸沙能力,有利于含沙水流與推移質(zhì)泥沙的淤積。

    2.6 丁壩(spur dike/groyne)

    丁壩是從河道岸邊伸出,在平面上與岸線構(gòu)成丁字形的河道整治結(jié)構(gòu)物。其主要作用有束窄河床、導(dǎo)水歸槽、調(diào)整流向、改變流速和導(dǎo)引泥沙等方面。按不同的設(shè)置目的、丁壩布置與水流不同的夾角,丁壩分為正丁壩、上挑丁壩和下挑丁壩。丁壩直接挑托水流遠(yuǎn)離堤岸,從而降低近岸流速和河床切應(yīng)力,促使泥沙在近岸淤積,以達(dá)到保護(hù)堤岸的目的。丁壩少數(shù)情況下也可能被完全淹沒,在壩頂形成過流,造成導(dǎo)流方向的改變。

    2.7 楔形丁壩(stream barb/bendway weir)

    楔形丁壩是一種低姿態(tài)的獨(dú)臂導(dǎo)流結(jié)構(gòu)物,主要用來保護(hù)堤岸免受沖刷。不同于丁壩單純挑托水流的功能,楔形丁壩壩體自河岸至突入河道一端形成一定坡度,分段淹沒在水下,從而在壩頂形成過流,即使在小流量情況下部分結(jié)構(gòu)物依然淹沒在水下,不但節(jié)省了材料,而且對(duì)水流阻礙作用小。壩體逆水流方向布置,與來流構(gòu)成一定角度(壩體軸線與沿壩根和岸堤交點(diǎn)切線的夾角),經(jīng)導(dǎo)流后的水流方向與壩體軸線相垂直,既不能因角度過小而使結(jié)構(gòu)物布置失去意義,也不能因角度過大而造成上游過分雍水并威脅結(jié)構(gòu)物自身安全。單個(gè)楔形丁壩的控制范圍有限,而彎道需要防護(hù)的岸坡則一直延伸至下游順直段,壩體通常沿彎道凹岸成群布置,分段保護(hù)岸坡免受沖蝕并發(fā)揮其聯(lián)合作用,最終將河道深泓線移至河道中心線的理想位置。

    楔形丁壩通過對(duì)彎道二次流結(jié)構(gòu)的調(diào)整,將高流速水流挑離河岸至河道中心線附近,繼而影響近岸的流速和河床切應(yīng)力分布。由于楔形丁壩被淹沒在水下,垂直壩頂軸線的過流與壩頭附近的收縮流交匯碰撞所形成的水流紊動(dòng)以及壩頂水躍均能就地削減河道內(nèi)的水流能量,避免將水流能量轉(zhuǎn)移至下游河段。合理的楔形丁壩布置會(huì)在上游岸線和結(jié)構(gòu)物附近形成一定范圍的回流區(qū)域,不僅能阻止高流速水流對(duì)河岸的淘刷,而且回水區(qū)穩(wěn)定的水流條件也有助于泥沙沉積。在楔形丁壩持續(xù)水力條件的作用下,丁壩壩頭和緊鄰結(jié)構(gòu)物下游形成沖坑,構(gòu)造出急流-緩流交替、深潭-淺灘交錯(cuò)的水流條件和地貌特征,能夠?yàn)榈讞锘螋~類提供良好的棲息庇護(hù)和覓食場(chǎng)所。此外,楔形丁壩群的聯(lián)合作用還能使河道深泓線從原先貼近凹岸岸線附近遷移至河道中心線的理想位置,從而有助于河勢(shì)的穩(wěn)定。

    3 楔形丁壩群導(dǎo)流特性

    3.1 193°彎曲水槽試驗(yàn)

    本節(jié)選取193°強(qiáng)彎曲水槽試驗(yàn)[18],參照國外楔形丁壩群的水力設(shè)計(jì)準(zhǔn)則,在該水槽凹岸布置1組(3個(gè))楔形丁壩群,通過三維數(shù)值模擬的方法對(duì)比結(jié)構(gòu)物布置前后水槽內(nèi)主要水力要素的變化,探討楔形丁壩群的水力特性,以此闡述楔形丁壩群在彎道凹岸岸坡治理的防護(hù)效果。該水槽試驗(yàn)的幾何尺寸和主要水力條件如下:寬度B=1.3 m的193°矩形彎曲水槽在其上游入口和下游出口分別連接9 m長和5 m長的矩形直段,彎段中心線的曲率半徑為Rc=1.7 m,水槽中心線總長為19.7 m,水槽底部為平底。入口流量Q=0.089 m3/s,平均水深H=0.159 m,平均流速V=0.43 m/s,弗勞德數(shù)Fr=0.34。水槽曲率半徑和寬度的比率為Rc/B=1.3,<3,屬強(qiáng)彎曲水槽,二次流效應(yīng)明顯。有關(guān)楔形丁壩的布置參照國外的相關(guān)設(shè)計(jì)準(zhǔn)則[13],結(jié)構(gòu)物突入寬度不超過河寬的1/4,突入角度控制在20°~30°之間,坡度≤1/5(入水一端高度為水深的一半),相鄰間距為長度的1~1.5倍,寬度通常為漂石直徑的2~3倍(本節(jié)考慮模擬需要,寬度取0.1 m)。193°強(qiáng)彎曲水槽內(nèi)楔形丁壩群的平面布置和剖面形狀如圖2所示。

    注:A,B,C表示楔形丁壩群編號(hào);B150表示150°彎段斷面;D0.5表示距彎段出口0.5 m下游順直段斷面,其他類推圖2 彎曲水槽內(nèi)楔形丁壩群布置Fig.2 Layout of stream barbs in channel bend

    3.2 三維數(shù)值模擬

    對(duì)水槽試驗(yàn)內(nèi)三維水流的模擬采用三維自由表面數(shù)值模型。模型求解以雷諾時(shí)均N-S方程和k-ε紊流模型作為控制方程的三維自由水面模型,數(shù)值方法見文獻(xiàn)[19-21],對(duì)雷諾時(shí)均N-S方程的求解采用壓力投影的方法。首先,對(duì)控制方程進(jìn)行分層積分,得到半離散化的方程,在前一時(shí)刻壓力值的基礎(chǔ)上,求解包含對(duì)流項(xiàng)和擴(kuò)散項(xiàng)的動(dòng)量方程,得到中間時(shí)刻的速度場(chǎng),其中,水平對(duì)流項(xiàng)和擴(kuò)散項(xiàng)采用Perot’s格式進(jìn)行離散求解;其次,求解由離散連續(xù)性方程和動(dòng)量方程聯(lián)立所得到的壓力Poisson方程,對(duì)中間時(shí)刻速度場(chǎng)進(jìn)行修正,得到新時(shí)刻的速度場(chǎng);然后,通過滿足自由水面邊界條件的水位演化方程計(jì)算新時(shí)刻的自由水面;最后,求解紊流封閉模型,計(jì)算渦黏系數(shù)。

    對(duì)水槽試驗(yàn)內(nèi)三維水流的模擬驗(yàn)證見文獻(xiàn)[19-21]。對(duì)布置有楔形丁壩群情況的模擬,除在彎道凹岸增設(shè)1組(3個(gè))楔形丁壩之外,其余如水槽尺寸、形態(tài)、網(wǎng)格劃分及相關(guān)水力參數(shù)等均保持不變。模擬采用標(biāo)準(zhǔn)k-ε紊流模型,時(shí)間步長取0.03 s,水平方向共劃分8 428個(gè)網(wǎng)格單元,垂向均勻分30層,糙率取值0.021,模擬時(shí)間150 s,直至水槽內(nèi)流動(dòng)達(dá)到恒定狀態(tài)為止。

    3.3 結(jié)果分析

    從水槽內(nèi)平面流場(chǎng)和流線分布、典型剖面上水面線和河床切應(yīng)力變化、典型斷面水位、垂向流速分布和二次流變化等幾個(gè)方面,對(duì)布置楔形丁壩群前后的模擬結(jié)果進(jìn)行分析和討論。

    圖3為楔形丁壩群作用下彎曲水槽內(nèi)的水位、局部流場(chǎng)和流線分布。從圖3中可以看出,受楔形丁壩群布置的影響,不論是水位還是流場(chǎng)或流線在結(jié)構(gòu)物附近及上方均有明顯的變化,水位在靠近凹岸一側(cè)的均勻梯度分布被打破,結(jié)構(gòu)物上游普遍形成雍水,水位抬升,過結(jié)構(gòu)物后水位下跌,從而水流能量得以部分削減。凹岸流線遇結(jié)構(gòu)物后方向發(fā)生改變,流向向凸岸一側(cè)偏轉(zhuǎn),在結(jié)構(gòu)物上方基本呈垂直結(jié)構(gòu)物走向流動(dòng),表明結(jié)構(gòu)物初步達(dá)到將主流挑離凹岸的目的。尤其是在楔形丁壩A和B之間還形成了明顯的局部回流區(qū)域,為泥沙在靠近岸坡一側(cè)沉積創(chuàng)造有利條件。

    圖3 布置楔形丁壩群后彎曲水槽內(nèi)的水位和流場(chǎng)(線)分布Fig.3 Contours of water level and streamlines or flow field in channel bend with stream barbs

    圖4 布置楔形丁壩群前后3個(gè)典型剖面上水面線與河床切應(yīng)力變化對(duì)比Fig.4 Water surfaces and bed shear stresses at three typical profiles along the channel bend with and without stream barbs

    圖4顯示布置楔形丁壩前后典型剖面上水面線和河床切應(yīng)力的變化對(duì)比,3條剖面從150°斷面起始,沿水槽軸線計(jì)算沿程距直至下游出口,分別距凹岸岸壁1/8,1/4,1/2河寬,在研究結(jié)構(gòu)物沿水流方向影響的同時(shí),也能注意其在橫向上的影響。與布置前相比,楔形丁壩群布置在結(jié)構(gòu)物上游形成雍水,過結(jié)構(gòu)物后水面迅速下跌,并且離岸壁距離越遠(yuǎn),對(duì)水面線的影響越弱,在水槽軸線剖面靠上游位置仍存在雍水現(xiàn)象,但水面線已相當(dāng)平緩。

    同樣,沿程河床切應(yīng)力變化表現(xiàn)出類似的現(xiàn)象,河床切應(yīng)力反映出近河床處的水流強(qiáng)度情況。從計(jì)算結(jié)果可以看出,在靠近岸壁1/8河寬位置,楔形丁壩布置高程相對(duì)較高,流速所受影響最大,雖然因?yàn)榻Y(jié)構(gòu)物的影響河床切應(yīng)力波動(dòng)較為劇烈,但切應(yīng)力大小相對(duì)結(jié)構(gòu)物布置前已有大幅減弱。1/4河寬剖面上的河床切應(yīng)力也存在微弱的波動(dòng),切應(yīng)力大小同樣也有所減弱,但在水槽軸線剖面上沿程河床切應(yīng)力明顯增強(qiáng),說明楔形丁壩群布置已將主流挑托遠(yuǎn)離岸壁的位置,主流向靠近水槽軸線位置發(fā)展,其長期作用能使河道深泓線穩(wěn)定在河道軸線位置。楔形丁壩群布置前岸壁附近的河床切應(yīng)力最大,靠近水槽軸線方向切應(yīng)力較小;而楔形丁壩群布置之后靠近岸壁處的河床切應(yīng)力明顯減弱,相應(yīng)地,水槽軸線位置的河床切應(yīng)力顯著增強(qiáng),說明楔形丁壩群布置使最大河床切應(yīng)力從凹岸岸壁處向水槽軸線位置偏移,從而起到保護(hù)堤岸免受沖蝕的目的。

    圖5為楔形丁壩群布置前后180°斷面上的水位和沿水深平均流速變化對(duì)比,圖5中左側(cè)是彎道凸岸,水面降低;而右側(cè)是凹岸,水面抬升。相較布置之前該斷面上水位均有所抬升,其中在靠近凹岸附近抬升較大。楔形丁壩群布置前斷面上流速大小保持均勻分布,布置后凹岸流速明顯減弱,流速梯度較強(qiáng);而凸岸流速隨著離凹岸距離越來越遠(yuǎn)則逐步增強(qiáng)。楔形丁壩群布置使該斷面水位雍高,凹岸流速降低,高流速水流被挑托至遠(yuǎn)離凹岸位置,但同時(shí)還應(yīng)防止凸岸一側(cè)被沖蝕的可能。

    圖5 布置楔形丁壩群前后180°斷面上的水位和水深方向平均流速對(duì)比Fig.5 Water level and depth-averaged velocity at section 180° with and without stream barbs

    圖6 布置楔形丁壩群前后180°斷面上典型測(cè)點(diǎn)垂線流速分布(縱向和橫向流速)對(duì)比Fig.6 Streamwise velocity and transverse velocity at typical points of section 180° with and without stream barbs

    圖6為布置楔形丁壩群前后180°斷面上典型測(cè)點(diǎn)位置垂向流速分布對(duì)比,2個(gè)典型測(cè)點(diǎn)均距離凸岸和凹岸0.477 m。從圖6中可以看出,布置楔形丁壩群之后,凹岸一側(cè)縱向流速大小有顯著減弱,垂向流速分布也有較大變化;凸岸一側(cè)縱向流速增大,垂向流速分布基本保持不變,可能帶來凸岸沖刷的不利影響。凹岸一側(cè)橫向流速分布基本與布置前呈對(duì)稱分布,表明表底層流速方向恰恰相反,即由原來的表層流向凹岸,底層向凸岸流動(dòng)的態(tài)勢(shì)轉(zhuǎn)變?yōu)楸韺酉蛲拱读鲃?dòng),而底層向凹岸流動(dòng)的情況??拷拱兑粋?cè)的橫向流速表層流速基本不變,中底層流速布置后有所增大。

    圖7為楔形丁壩群布置前后180°斷面上二次流流線分布變化,左側(cè)水位抬升為凹岸,右側(cè)水位降低為凸岸。從圖7中可以看出,布置楔形丁壩群之前整個(gè)斷面只存在一個(gè)逆時(shí)針方向的二次環(huán)流,表層向凹岸流動(dòng),底層流向凸岸。而布置楔形丁壩群之后,在斷面靠近凹岸一側(cè)新形成了一個(gè)順時(shí)針方向的環(huán)流渦,與布置前的二次渦流向正好相反,底層向凹岸流動(dòng),表層流向凸岸,從而削弱主渦環(huán)流將河床上的沙粒向凸岸一側(cè)輸移,達(dá)到保護(hù)凹岸岸坡免受沖蝕的目的。圖7中的二次流環(huán)流結(jié)果說明,楔形丁壩群布置能有效地改變彎曲水槽內(nèi)的二次流分布,進(jìn)而影響彎道內(nèi)的流速和切應(yīng)力分布。

    圖7 楔形丁壩群布置前后180°斷面二次流流線分布變化Fig.7 Secondary flow streamlines at section 180° with and without stream barbs

    4 結(jié) 論

    在當(dāng)前河流再自然化,強(qiáng)調(diào)適度人工干預(yù)措施進(jìn)行河流生態(tài)修復(fù)以恢復(fù)河流生態(tài)系統(tǒng)結(jié)構(gòu)和功能的背景下,應(yīng)拓展傳統(tǒng)丁壩水力學(xué)的生態(tài)效益,完善并推廣我國河流生態(tài)修復(fù)領(lǐng)域河道內(nèi)流動(dòng)控制結(jié)構(gòu)物的研究工作,進(jìn)行科學(xué)分類與命名,并研究其水力特性。

    (1)通過對(duì)彎道內(nèi)楔形丁壩群水力特性的初步研究,從水面形態(tài)、流速和切應(yīng)力及二次環(huán)流等幾方面變化對(duì)比分析可以發(fā)現(xiàn),在彎道凹岸岸坡布置楔形丁壩群能達(dá)到將高流速水流挑離凹岸,避免岸坡遭受沖蝕的目的,其作用機(jī)制主要是通過調(diào)整彎道內(nèi)的二次流結(jié)構(gòu),以發(fā)揮其挑托水流的目的。楔形丁壩群布置會(huì)造成彎道內(nèi)水位雍高,高流速水流向凸岸一側(cè)偏移,其雍水高度和偏移程度與楔形丁壩群的布置方式密切相關(guān)。

    (2)本文所采用的三維自由水面模型在成功模擬彎曲水槽三維水流結(jié)構(gòu)的基礎(chǔ)上,同樣適用于對(duì)河道內(nèi)流動(dòng)控制結(jié)構(gòu)物作用下水槽內(nèi)三維水力特性的模擬,該數(shù)值模型可以作為進(jìn)一步研究結(jié)構(gòu)物導(dǎo)流特性的有效工具。

    猜你喜歡
    凹岸丁壩楔形
    河道防洪治理工程丁壩設(shè)置應(yīng)注意的問題
    山東水利(2020年11期)2020-12-22 07:45:06
    History of the Alphabet
    鋼絲繩楔形接頭連接失效分析與預(yù)防
    Eight Surprising Foods You’er Never Tried to Grill Before
    河流凹凸岸芻議
    例談河流凹凸岸的判斷方法
    例談河流凹凸岸的判斷方法
    山區(qū)河流上下雙丁壩回流區(qū)水沙特性淺探
    江西建材(2018年1期)2018-04-04 05:26:28
    腹腔鏡下胃楔形切除術(shù)治療胃間質(zhì)瘤30例
    考慮水流紊動(dòng)的丁壩下游回流區(qū)水流挾沙力研究
    国内久久婷婷六月综合欲色啪| 久久精品91无色码中文字幕| 国产精品亚洲一级av第二区| 久久久久免费精品人妻一区二区 | 午夜成年电影在线免费观看| 欧美色视频一区免费| 久久久久久久午夜电影| 久久国产乱子伦精品免费另类| 午夜免费成人在线视频| 国产99白浆流出| 黄网站色视频无遮挡免费观看| 成年版毛片免费区| 久久精品亚洲精品国产色婷小说| 亚洲精品av麻豆狂野| 国产激情欧美一区二区| 精品久久久久久久人妻蜜臀av| 日本三级黄在线观看| 91成人精品电影| 老汉色∧v一级毛片| 色老头精品视频在线观看| 国产午夜福利久久久久久| 国内少妇人妻偷人精品xxx网站 | 国产人伦9x9x在线观看| 国产97色在线日韩免费| 香蕉国产在线看| or卡值多少钱| 日日夜夜操网爽| 两个人看的免费小视频| 最近最新中文字幕大全电影3 | 欧美日韩亚洲国产一区二区在线观看| 麻豆国产av国片精品| 久久天躁狠狠躁夜夜2o2o| 亚洲五月天丁香| 国产精品 国内视频| 夜夜爽天天搞| 少妇裸体淫交视频免费看高清 | 亚洲av成人不卡在线观看播放网| 97碰自拍视频| 老熟妇仑乱视频hdxx| 国产视频一区二区在线看| 久热爱精品视频在线9| 婷婷亚洲欧美| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三区av网在线观看| 亚洲avbb在线观看| 国产高清有码在线观看视频 | 久9热在线精品视频| 亚洲av成人av| 在线观看66精品国产| 久久伊人香网站| 又大又爽又粗| 亚洲激情在线av| 国产伦在线观看视频一区| 国产精品美女特级片免费视频播放器 | 午夜久久久在线观看| 少妇的丰满在线观看| 首页视频小说图片口味搜索| 国产视频内射| 老司机午夜十八禁免费视频| 国产亚洲欧美精品永久| 久久午夜综合久久蜜桃| 熟女少妇亚洲综合色aaa.| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产99精品国产亚洲性色| 99久久国产精品久久久| 国产成人av教育| 精品一区二区三区四区五区乱码| 中文字幕精品亚洲无线码一区 | 男女那种视频在线观看| 黄片播放在线免费| 国产人伦9x9x在线观看| 亚洲人成网站高清观看| 美国免费a级毛片| 婷婷精品国产亚洲av| 国产乱人伦免费视频| 亚洲色图av天堂| 一级黄色大片毛片| 亚洲人成网站在线播放欧美日韩| 成人特级黄色片久久久久久久| 亚洲五月色婷婷综合| 少妇 在线观看| 精品少妇一区二区三区视频日本电影| 妹子高潮喷水视频| 好男人电影高清在线观看| 国产三级在线视频| 搡老岳熟女国产| 欧美在线黄色| 岛国视频午夜一区免费看| 久久久久久大精品| 日本成人三级电影网站| 妹子高潮喷水视频| 此物有八面人人有两片| 国产精品香港三级国产av潘金莲| 精品免费久久久久久久清纯| 久久九九热精品免费| 女人被狂操c到高潮| 高潮久久久久久久久久久不卡| 免费看日本二区| 色播亚洲综合网| 亚洲专区中文字幕在线| 神马国产精品三级电影在线观看 | 老司机在亚洲福利影院| 日韩中文字幕欧美一区二区| 色综合亚洲欧美另类图片| 啪啪无遮挡十八禁网站| 熟女少妇亚洲综合色aaa.| 欧美不卡视频在线免费观看 | 久久久久久大精品| 非洲黑人性xxxx精品又粗又长| 国产精品国产高清国产av| 一边摸一边做爽爽视频免费| 亚洲专区中文字幕在线| 欧美日韩黄片免| 男女之事视频高清在线观看| 亚洲精品中文字幕一二三四区| 俺也久久电影网| 99久久精品国产亚洲精品| 亚洲成国产人片在线观看| 亚洲五月婷婷丁香| 久久久久免费精品人妻一区二区 | 色综合婷婷激情| 韩国av一区二区三区四区| 亚洲午夜精品一区,二区,三区| 老司机在亚洲福利影院| 亚洲中文字幕一区二区三区有码在线看 | 国产又爽黄色视频| 一a级毛片在线观看| 在线免费观看的www视频| 久久久国产精品麻豆| 黄频高清免费视频| www.www免费av| 黄色片一级片一级黄色片| 首页视频小说图片口味搜索| 国产aⅴ精品一区二区三区波| 真人做人爱边吃奶动态| 亚洲va日本ⅴa欧美va伊人久久| 欧美精品亚洲一区二区| 动漫黄色视频在线观看| 亚洲一区二区三区色噜噜| 韩国精品一区二区三区| 一区二区日韩欧美中文字幕| 一二三四社区在线视频社区8| 国产黄a三级三级三级人| 午夜福利免费观看在线| 大型av网站在线播放| 最近在线观看免费完整版| 精品久久蜜臀av无| 级片在线观看| 黄色毛片三级朝国网站| 成人永久免费在线观看视频| 亚洲人成网站在线播放欧美日韩| 欧美日本亚洲视频在线播放| 久久精品夜夜夜夜夜久久蜜豆 | 久久国产精品影院| 国产99白浆流出| 亚洲 国产 在线| 国产高清激情床上av| 国产av不卡久久| 日本 欧美在线| 国产成人一区二区三区免费视频网站| 精品久久蜜臀av无| 看免费av毛片| 欧美绝顶高潮抽搐喷水| 一区二区三区激情视频| 国产精品久久久久久亚洲av鲁大| 国产av一区二区精品久久| 久久久精品欧美日韩精品| 18禁裸乳无遮挡免费网站照片 | 中出人妻视频一区二区| 精品一区二区三区av网在线观看| 国产伦在线观看视频一区| 色哟哟哟哟哟哟| avwww免费| 成人免费观看视频高清| 男人的好看免费观看在线视频 | 后天国语完整版免费观看| 少妇的丰满在线观看| 中文字幕久久专区| 美女国产高潮福利片在线看| 日韩免费av在线播放| 男女下面进入的视频免费午夜 | 精品熟女少妇八av免费久了| 老司机午夜十八禁免费视频| 性欧美人与动物交配| 午夜久久久久精精品| 欧美日本视频| 身体一侧抽搐| 久久国产乱子伦精品免费另类| 国产精品,欧美在线| 亚洲精品粉嫩美女一区| 18禁美女被吸乳视频| 麻豆国产av国片精品| 精品久久久久久久久久久久久 | 伊人久久大香线蕉亚洲五| 一进一出好大好爽视频| 黄频高清免费视频| 亚洲欧美激情综合另类| 我的亚洲天堂| 又紧又爽又黄一区二区| 久久国产精品男人的天堂亚洲| 亚洲熟妇熟女久久| 亚洲精品久久国产高清桃花| 亚洲欧美精品综合久久99| 啦啦啦韩国在线观看视频| www日本在线高清视频| 此物有八面人人有两片| 亚洲第一青青草原| 国产黄片美女视频| 丝袜在线中文字幕| 久久精品夜夜夜夜夜久久蜜豆 | 日韩中文字幕欧美一区二区| av天堂在线播放| 免费av毛片视频| 欧美成人性av电影在线观看| 一本一本综合久久| 十分钟在线观看高清视频www| 校园春色视频在线观看| 精品卡一卡二卡四卡免费| 老汉色∧v一级毛片| 亚洲精品在线观看二区| 国产精品九九99| 日韩高清综合在线| 亚洲第一欧美日韩一区二区三区| 久久青草综合色| 久久香蕉精品热| 久久婷婷人人爽人人干人人爱| 久久人妻福利社区极品人妻图片| 久久午夜亚洲精品久久| 婷婷丁香在线五月| 真人做人爱边吃奶动态| 黄网站色视频无遮挡免费观看| 国产av一区二区精品久久| 成人午夜高清在线视频 | 精品国产亚洲在线| 久久精品国产亚洲av香蕉五月| 成在线人永久免费视频| 在线观看日韩欧美| 中文字幕最新亚洲高清| 欧美日韩瑟瑟在线播放| 国产精品影院久久| 亚洲精品美女久久av网站| a级毛片在线看网站| 首页视频小说图片口味搜索| 免费一级毛片在线播放高清视频| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 国产精品九九99| 久久久国产成人免费| 国产激情久久老熟女| 欧美久久黑人一区二区| 亚洲欧美精品综合一区二区三区| 男人舔女人下体高潮全视频| 免费女性裸体啪啪无遮挡网站| 在线av久久热| 国产亚洲精品综合一区在线观看 | 男人舔女人的私密视频| 亚洲成人国产一区在线观看| 亚洲国产欧美一区二区综合| 午夜福利成人在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美日韩无卡精品| 日韩欧美免费精品| 国产v大片淫在线免费观看| 亚洲最大成人中文| 欧美又色又爽又黄视频| 色综合站精品国产| 午夜老司机福利片| 精品国产国语对白av| 一进一出抽搐gif免费好疼| 哪里可以看免费的av片| 91国产中文字幕| 中文字幕最新亚洲高清| 精品不卡国产一区二区三区| 一区福利在线观看| 亚洲av熟女| 法律面前人人平等表现在哪些方面| 国内精品一区二区在线观看| 亚洲真实伦在线观看| av在线天堂中文字幕| 国产69精品久久久久777片| 国产精品亚洲美女久久久| 欧美日本亚洲视频在线播放| 亚洲性夜色夜夜综合| 桃色一区二区三区在线观看| 久久久久免费精品人妻一区二区| 午夜福利18| 天美传媒精品一区二区| 99久国产av精品| 欧美绝顶高潮抽搐喷水| 日韩中字成人| 免费人成视频x8x8入口观看| 韩国av在线不卡| 国产伦一二天堂av在线观看| 国产真实伦视频高清在线观看| 日本爱情动作片www.在线观看 | 五月伊人婷婷丁香| 一个人观看的视频www高清免费观看| 成熟少妇高潮喷水视频| 欧美不卡视频在线免费观看| 精品一区二区免费观看| 最后的刺客免费高清国语| 丝袜喷水一区| 悠悠久久av| 国产精品爽爽va在线观看网站| 麻豆精品久久久久久蜜桃| 亚洲性久久影院| 少妇丰满av| 性插视频无遮挡在线免费观看| 国产久久久一区二区三区| 国产精品久久久久久久电影| 久久精品影院6| 欧美xxxx性猛交bbbb| 国产激情偷乱视频一区二区| 欧美日韩国产亚洲二区| 国产三级中文精品| 亚洲无线在线观看| 久久国内精品自在自线图片| 国产精品嫩草影院av在线观看| 精品久久久噜噜| 亚洲成av人片在线播放无| 麻豆久久精品国产亚洲av| 亚洲自偷自拍三级| av女优亚洲男人天堂| 亚洲中文日韩欧美视频| 国产午夜精品久久久久久一区二区三区 | www日本黄色视频网| 一个人免费在线观看电影| 蜜桃亚洲精品一区二区三区| 网址你懂的国产日韩在线| 亚洲在线自拍视频| 亚洲精品一卡2卡三卡4卡5卡| 久久99热这里只有精品18| 色播亚洲综合网| a级毛色黄片| 日韩精品中文字幕看吧| 99久久精品热视频| 日本欧美国产在线视频| 精品久久久久久久人妻蜜臀av| 久久国产乱子免费精品| 超碰av人人做人人爽久久| 最新在线观看一区二区三区| 久久亚洲国产成人精品v| 无遮挡黄片免费观看| 国产一区二区三区在线臀色熟女| 老司机午夜福利在线观看视频| 成年版毛片免费区| 日韩欧美三级三区| 欧美色视频一区免费| 国产欧美日韩精品一区二区| 亚洲真实伦在线观看| av视频在线观看入口| 久久久久久久久久久丰满| 欧美日韩一区二区视频在线观看视频在线 | 男女下面进入的视频免费午夜| 国产精品久久久久久av不卡| 久久精品国产亚洲av天美| 亚洲精品一卡2卡三卡4卡5卡| 香蕉av资源在线| 国模一区二区三区四区视频| 精品一区二区三区av网在线观看| 12—13女人毛片做爰片一| 中文字幕熟女人妻在线| 久久草成人影院| 一级黄色大片毛片| 国产精品三级大全| 精品无人区乱码1区二区| 国产精品99久久久久久久久| 久久欧美精品欧美久久欧美| 亚洲综合色惰| 嫩草影院精品99| 国产精品一及| 中文字幕av成人在线电影| 天天一区二区日本电影三级| 人人妻人人看人人澡| 天堂网av新在线| 黄色日韩在线| 99久久精品热视频| 男插女下体视频免费在线播放| 真实男女啪啪啪动态图| ponron亚洲| 观看美女的网站| 久久久久精品国产欧美久久久| 中文字幕av成人在线电影| 免费大片18禁| 色尼玛亚洲综合影院| 91久久精品电影网| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 三级男女做爰猛烈吃奶摸视频| 久久久久久久久大av| 一个人观看的视频www高清免费观看| 日韩欧美国产在线观看| 亚洲中文字幕日韩| 丰满人妻一区二区三区视频av| 国产精华一区二区三区| 国产欧美日韩精品一区二区| 亚洲激情五月婷婷啪啪| 偷拍熟女少妇极品色| 成人亚洲精品av一区二区| 99视频精品全部免费 在线| 国产成人福利小说| 97超级碰碰碰精品色视频在线观看| 日韩精品中文字幕看吧| 我的女老师完整版在线观看| 美女cb高潮喷水在线观看| 可以在线观看毛片的网站| 国产乱人视频| 日韩成人av中文字幕在线观看 | av在线蜜桃| 国产一区二区激情短视频| 日韩 亚洲 欧美在线| 99热全是精品| 精品人妻一区二区三区麻豆 | av国产免费在线观看| 六月丁香七月| 成人欧美大片| 成人三级黄色视频| av卡一久久| a级一级毛片免费在线观看| 搡老岳熟女国产| 老司机影院成人| 午夜视频国产福利| 欧美bdsm另类| 最近在线观看免费完整版| 欧美激情国产日韩精品一区| 日韩制服骚丝袜av| 亚洲人与动物交配视频| 三级经典国产精品| 免费不卡的大黄色大毛片视频在线观看 | 变态另类丝袜制服| 一级黄片播放器| 免费av不卡在线播放| or卡值多少钱| 老师上课跳d突然被开到最大视频| av.在线天堂| 菩萨蛮人人尽说江南好唐韦庄 | 欧美不卡视频在线免费观看| 色哟哟·www| 国内久久婷婷六月综合欲色啪| 波多野结衣高清无吗| 亚洲五月天丁香| 中国美白少妇内射xxxbb| 亚洲四区av| 国产探花极品一区二区| 国产成人福利小说| 精品人妻一区二区三区麻豆 | 亚洲无线在线观看| 中文字幕人妻熟人妻熟丝袜美| 精品人妻视频免费看| 好男人在线观看高清免费视频| 国产精品一区二区三区四区免费观看 | 毛片一级片免费看久久久久| 欧美另类亚洲清纯唯美| 精品熟女少妇av免费看| 欧美性猛交黑人性爽| 亚洲人成网站在线播| 久久精品人妻少妇| 午夜免费激情av| 亚洲国产精品合色在线| 亚洲最大成人av| 久久久精品大字幕| 少妇丰满av| 久久人人爽人人片av| 悠悠久久av| 国产精品久久久久久av不卡| 中国美白少妇内射xxxbb| 深夜精品福利| 亚洲不卡免费看| 亚洲美女黄片视频| 好男人在线观看高清免费视频| 久久久久国产精品人妻aⅴ院| 我要看日韩黄色一级片| 欧美一级a爱片免费观看看| 亚洲美女视频黄频| 欧美丝袜亚洲另类| 欧美日韩一区二区视频在线观看视频在线 | 国产一区二区激情短视频| 国产在视频线在精品| 日本色播在线视频| 美女 人体艺术 gogo| 国产极品精品免费视频能看的| 男女视频在线观看网站免费| 91在线观看av| 国产欧美日韩精品亚洲av| 久久久久久久久久久丰满| av国产免费在线观看| 中文字幕精品亚洲无线码一区| 欧美日韩一区二区视频在线观看视频在线 | 国产精品一二三区在线看| 免费在线观看影片大全网站| 夜夜爽天天搞| 国产探花极品一区二区| 中国国产av一级| 日韩一区二区视频免费看| 国产男人的电影天堂91| 在线播放国产精品三级| 亚洲成人中文字幕在线播放| 国产日本99.免费观看| 两个人的视频大全免费| 亚洲欧美清纯卡通| 成人亚洲欧美一区二区av| 久久韩国三级中文字幕| 人妻制服诱惑在线中文字幕| 狂野欧美激情性xxxx在线观看| 可以在线观看毛片的网站| 十八禁网站免费在线| 亚洲天堂国产精品一区在线| 欧美bdsm另类| 给我免费播放毛片高清在线观看| 亚洲欧美成人精品一区二区| 国产亚洲91精品色在线| 男女那种视频在线观看| 国产精品av视频在线免费观看| 热99re8久久精品国产| 国产午夜精品论理片| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 在线天堂最新版资源| 国产精品精品国产色婷婷| 免费人成在线观看视频色| 长腿黑丝高跟| 日本一本二区三区精品| 精品久久久久久久久亚洲| 美女高潮的动态| 亚洲乱码一区二区免费版| 99热6这里只有精品| ponron亚洲| 久久午夜亚洲精品久久| 国产熟女欧美一区二区| 亚洲五月天丁香| 亚洲精品一卡2卡三卡4卡5卡| 天天躁日日操中文字幕| 久久久午夜欧美精品| 秋霞在线观看毛片| 自拍偷自拍亚洲精品老妇| 免费看av在线观看网站| 色在线成人网| 99热这里只有是精品在线观看| 自拍偷自拍亚洲精品老妇| 秋霞在线观看毛片| 蜜臀久久99精品久久宅男| 夜夜夜夜夜久久久久| 日产精品乱码卡一卡2卡三| 国产麻豆成人av免费视频| 日韩精品青青久久久久久| 51国产日韩欧美| 中文字幕av成人在线电影| 精品人妻熟女av久视频| 国产精品久久久久久久电影| 国产精品久久久久久av不卡| 亚洲精品一卡2卡三卡4卡5卡| 久久天躁狠狠躁夜夜2o2o| 国产视频一区二区在线看| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| 色播亚洲综合网| .国产精品久久| 一级毛片我不卡| 国产一级毛片七仙女欲春2| 一边摸一边抽搐一进一小说| 国产精品1区2区在线观看.| 美女内射精品一级片tv| 亚洲最大成人中文| 精品免费久久久久久久清纯| 亚洲最大成人中文| 欧美zozozo另类| 免费无遮挡裸体视频| 欧美成人一区二区免费高清观看| 国产熟女欧美一区二区| 久久午夜亚洲精品久久| 国产高潮美女av| 午夜爱爱视频在线播放| 国产爱豆传媒在线观看| 亚洲人与动物交配视频| 97碰自拍视频| 色视频www国产| 黄片wwwwww| 天堂影院成人在线观看| 黑人高潮一二区| 免费人成在线观看视频色| 久久久久久久午夜电影| 日日摸夜夜添夜夜添av毛片| 免费看美女性在线毛片视频| 日韩,欧美,国产一区二区三区 | 大香蕉久久网| 99国产精品一区二区蜜桃av| 亚洲av二区三区四区| 国模一区二区三区四区视频| 亚洲中文字幕一区二区三区有码在线看| 国产人妻一区二区三区在| 日韩精品青青久久久久久| 色视频www国产| 91av网一区二区| 成年av动漫网址| 国产精品,欧美在线| 亚洲三级黄色毛片| or卡值多少钱| 国产精品乱码一区二三区的特点| 亚洲欧美日韩卡通动漫| 亚洲专区国产一区二区| 国产亚洲精品久久久com| ponron亚洲| 成年免费大片在线观看| 精华霜和精华液先用哪个| 国产精品久久电影中文字幕| 色吧在线观看| 日日摸夜夜添夜夜添小说| 最近视频中文字幕2019在线8| 亚洲最大成人中文| 最后的刺客免费高清国语| 不卡一级毛片| 午夜激情欧美在线| av在线观看视频网站免费| 欧美极品一区二区三区四区| 亚洲av.av天堂|