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

    中國(guó)常用人工魚礁流場(chǎng)效應(yīng)的比較分析

    2017-10-13 06:27:27唐衍力龍翔宇王欣欣姜昭陽(yáng)張同征
    關(guān)鍵詞:上升流礁體場(chǎng)效應(yīng)

    唐衍力,龍翔宇,王欣欣,姜昭陽(yáng),程 暉,張同征

    ?

    中國(guó)常用人工魚礁流場(chǎng)效應(yīng)的比較分析

    唐衍力1,龍翔宇1,王欣欣1※,姜昭陽(yáng)2,程 暉1,張同征1

    (1. 中國(guó)海洋大學(xué)水產(chǎn)學(xué)院,青島 266003;2. 山東大學(xué)(威海)海洋學(xué)院,威海264209)

    不同結(jié)構(gòu)的人工魚礁在海中會(huì)產(chǎn)生不同的流場(chǎng)效應(yīng),為了對(duì)比分析不同結(jié)構(gòu)人工魚礁的流場(chǎng)效應(yīng)差異,該文利用CFX軟件,對(duì)6類礁型18種中國(guó)常用人工魚礁進(jìn)行了數(shù)值模擬。引入了2種相對(duì)評(píng)價(jià)方法(比礁高和比混凝土體積)和3個(gè)評(píng)價(jià)指標(biāo)(礁體中垂面上的上升流面積、背渦流面積和上升流高度)。首先分析了評(píng)價(jià)指標(biāo)適用的相對(duì)評(píng)價(jià)方法,然后對(duì)18種魚礁進(jìn)行比較研究,最后從每類礁型分別選出一個(gè)礁體為代表研究它們?cè)诓煌魉傧拢?.2, 0.4, 0.6, 0.8和 1.0 m/s)上升流和背渦流的差異。結(jié)果表明:研究上升流和背渦流面積時(shí)比混凝土體積法更科學(xué),研究上升流高度時(shí)比礁高法更有效;無(wú)論上升流還是背渦流,三角型礁的相對(duì)面積都為最大,復(fù)合型礁次之,框架性礁最??;上升流高度、上升流面積和背渦流面積都不隨來(lái)流速度變化;最大上升流流速與來(lái)流速度呈線性關(guān)系,其斜率在不同礁體間存在差異;該研究以期為不同海域不同要求下人工魚礁的選擇和設(shè)計(jì)提供了參考。

    魚礁;流場(chǎng);數(shù)值模擬;人工;比礁高,比混凝土體積

    0 引 言

    由于過(guò)度捕撈和環(huán)境污染,海洋生態(tài)環(huán)境遭到破壞,漁業(yè)資源瀕臨枯竭,利用人工魚礁改善和優(yōu)化海域生態(tài)環(huán)境已被廣泛采用[1]。美國(guó)、日本、韓國(guó)、澳大利亞等一些漁業(yè)發(fā)達(dá)國(guó)家,在利用人工魚礁修復(fù)和改善漁業(yè)生態(tài)環(huán)境、開發(fā)和保護(hù)漁業(yè)資源等方面,取得了諸多成效[2-5]。人工魚礁可為海洋生物提供索餌、繁殖、生長(zhǎng)等場(chǎng)所,達(dá)到修復(fù)漁業(yè)資源和改善海洋生態(tài)環(huán)境的目的[6]。人工魚礁投放后,會(huì)產(chǎn)生流場(chǎng)效應(yīng)、餌料效應(yīng)和避敵效應(yīng)。由于魚礁對(duì)水流的阻礙作用導(dǎo)致周圍的壓力場(chǎng)發(fā)生變化,形成了新的流場(chǎng)分布[7]。流體經(jīng)過(guò)礁體時(shí)產(chǎn)生的上升流和背渦流可促進(jìn)上下層海水的交換、加快營(yíng)養(yǎng)物質(zhì)循環(huán)、提高海域的初級(jí)生產(chǎn)力,從而誘集魚類。多數(shù)魚類喜棲息于流速緩慢的渦流區(qū),渦流還可造成浮游生物、甲殼類和魚類的物理性聚集。因此,上升流和背渦流規(guī)模是衡量人工魚礁流場(chǎng)效應(yīng)的重要指標(biāo),對(duì)人工魚礁物理環(huán)境功能造成技術(shù)的研究具有重要意義[8]。

    日本學(xué)者上北征男等對(duì)人工魚礁著底沖擊力和海底礁各種地形產(chǎn)生的涌升流進(jìn)行了研究[9-11]。鄭延璇等通過(guò)數(shù)值模擬和粒子圖像測(cè)速分析了3種疊放方式的圓管型人工魚礁周圍的流場(chǎng)分布[12]。劉洪生等通過(guò)風(fēng)洞試驗(yàn)和數(shù)值模擬研究了正方體、三棱柱及金字塔型人工魚礁周圍的流場(chǎng)[13]。崔勇等采用數(shù)值模擬方法研究了布設(shè)間距對(duì)人工魚礁流場(chǎng)效應(yīng)的影響,并與劉洪生等利用風(fēng)洞試驗(yàn)研究不同類型人工魚礁單體和不同組合正方體模型的流場(chǎng)效應(yīng)結(jié)果進(jìn)行了對(duì)比,二者確定2個(gè)礁體的最佳布設(shè)間距為礁體尺寸的1~1.5倍,結(jié)果一致[14-15]。先前對(duì)人工魚礁流場(chǎng)效應(yīng)的研究只局限于一種或幾種礁體,并且缺少不同形狀間的橫向?qū)Ρ取1疚倪x取了6類18人工魚礁進(jìn)行數(shù)值模擬,對(duì)比研究了其流場(chǎng)效應(yīng)的差異。

    本文基于計(jì)算流體力學(xué)理論,建立了單體人工魚礁周圍的湍流數(shù)值模型,采用有限體積法(finite volume method,F(xiàn)VM)對(duì)三維Navier-Stokes方程進(jìn)行離散求解。鑒于許多學(xué)者已經(jīng)證實(shí)了計(jì)算流體動(dòng)力學(xué)軟件CFX的可信度[16-19],運(yùn)用CFX分析了中國(guó)常用的6類礁型18種人工魚礁的流場(chǎng)分布特性。用上升流和背渦流作為人工魚礁流場(chǎng)效應(yīng)的衡量指標(biāo),采用比礁高(上升流高度、上升流面積或背渦流面積與礁高之比)和比混凝土體積(上升流高度、上升流面積或背渦流面積與礁體混凝土體積之比)2種方法,系統(tǒng)分析了不同人工魚礁的流場(chǎng)效應(yīng),以期為人工魚礁的設(shè)計(jì)、選型和生態(tài)效益評(píng)估提供參考依據(jù)。

    1 數(shù)值模擬方法

    1.1 控制方程

    假設(shè)流體為不可壓縮流動(dòng),三維笛卡爾坐標(biāo)系中,連續(xù)性方程和雷諾時(shí)均Navier–Stokes方程為[20]

    (2)

    式中為流體密度,kg/m3;,u(,=1,2,3,)為直角坐標(biāo)系,和3個(gè)方向的平均流速,m/s;為平均壓強(qiáng),Pa;為動(dòng)力粘度,Pa·s;為雷諾應(yīng)力,N。

    (4)

    1.2 數(shù)值模擬模型

    圖1為中國(guó)常用的18種人工魚礁(artificialreef,AR)模型[25],以形狀將其分為方型(AR01~AR05)、框架型(AR06~AR08)、棱臺(tái)型(AR09和AR10)、三角型(AR11和AR12)、圓柱型(AR13~AR15)和復(fù)合型(AR16~AR18)等6類礁型,其相關(guān)參數(shù)如表1所示。

    如圖2所示,將礁體靜置海底,豎直方向?yàn)榻阁w高度(m),沿水流方向?yàn)榻阁w長(zhǎng)度(m),與水流方向和豎直方向垂直的是礁體寬度(m)。對(duì)于側(cè)面相同的方型礁、框架型礁、棱臺(tái)型礁和復(fù)合型礁,選擇任一側(cè)面為迎流面,而AR18選擇船型側(cè)面為迎流面;對(duì)于三角型礁,選擇斜板面為迎流面;對(duì)于圓柱型礁,選擇管壁為迎 流面。

    1.3 計(jì)算域和邊界條件

    流場(chǎng)計(jì)算區(qū)域包括入口邊界(Inlet)、出口邊界(Outlet)、壁面(Wall)和開放面(Opening),為準(zhǔn)確反映礁體周圍流場(chǎng)情況,參考文獻(xiàn)[26],計(jì)算域設(shè)定長(zhǎng)度為魚礁前3倍礁長(zhǎng),后10倍礁長(zhǎng),寬度為5倍礁寬,高度為5倍礁高,即14×5×5,如圖3a所示。

    圖1 18種常用人工魚礁

    表1 18種常用人工魚礁相關(guān)參數(shù)

    文中所設(shè)定的邊界條件如下

    1)入口邊界為速度入口,設(shè)置來(lái)流速度,計(jì)算并給定邊界上湍動(dòng)能和湍動(dòng)耗散率;

    2)出口邊界為壓力出口,相對(duì)壓強(qiáng)設(shè)定為0 Pa;

    3)計(jì)算域的頂面和兩側(cè)面選擇開放面,設(shè)置為可移動(dòng)的壁面;

    4)計(jì)算域底面和礁體表面選擇壁面,設(shè)置為光滑和無(wú)滑移。

    1.4 網(wǎng)格劃分

    使用ANASYS Workbench Mesh模塊對(duì)人工魚礁和計(jì)算流域進(jìn)行四面體單元非結(jié)構(gòu)化網(wǎng)格劃分。由于礁體周圍流場(chǎng)變化梯度較大,為提高數(shù)值模擬的計(jì)算精度,在礁體周圍設(shè)置一個(gè)加密區(qū)域,范圍為5×3×3,即5倍礁長(zhǎng)、3倍礁寬和3倍礁高。對(duì)于加密區(qū)域,單元尺寸設(shè)置為0.1 m,其它區(qū)域使用相對(duì)稀疏的網(wǎng)格,如圖3b所示。網(wǎng)格無(wú)關(guān)性的預(yù)試驗(yàn)表明,當(dāng)網(wǎng)格數(shù)量在500萬(wàn)以上時(shí),網(wǎng)格數(shù)量對(duì)結(jié)果幾乎沒(méi)有影響,最終網(wǎng)格劃分的數(shù)量介于5~20×106。

    圖3 計(jì)算域,邊界條件和網(wǎng)格劃分

    1.5 流場(chǎng)效應(yīng)評(píng)價(jià)指標(biāo)和計(jì)算方法

    由于上升流和背渦流體積不易計(jì)算,選擇它們?cè)谥写姑嫔系拿娣e進(jìn)行研究。數(shù)值計(jì)算中,假定坐標(biāo)原點(diǎn)位于礁體底面中心,來(lái)流方向沿軸方向,軸為豎直方向。選取=0中垂面(平面)作為研究的基準(zhǔn)面(如圖4),將方向流速分量與來(lái)流速度之比等于或大于10%的區(qū)域定義為上升流區(qū)域[27]。由于背渦流沒(méi)有明確數(shù)量和方向上的定義,本文將在方向速度分量小于0的區(qū)域作為背渦流區(qū)域進(jìn)行簡(jiǎn)化計(jì)算[18]。上升流高度都是以礁體底部為零點(diǎn)進(jìn)行測(cè)量的。

    數(shù)值模擬中,用CFX對(duì)人工魚礁的流場(chǎng)效應(yīng)進(jìn)行分析。將18種礁體導(dǎo)入CFX進(jìn)行計(jì)算,計(jì)算結(jié)果導(dǎo)入CFD- Post后處理,上升流和背渦流面積采用Photoshop像素法計(jì)算[28],對(duì)所得的上升流和背渦流面積進(jìn)行比礁高和比混凝土體積2種相對(duì)化處理得到相對(duì)面積,然后進(jìn)行定量分析比較。從每種類型中選出一個(gè)代表礁體,進(jìn)行多流速模擬計(jì)算,流速分別為0.2, 0.4, 0.6, 0.8和1.0 m/s,研究流速對(duì)不同形狀人工魚礁的上升流和背渦流相關(guān)特性的影響。

    圖4 定義的中垂面

    2 結(jié)果與分析

    2.1 上升流的流場(chǎng)效應(yīng)

    圖5是18種礁體在中垂面上的上升流相對(duì)面積對(duì)比圖,可以看出18種礁體均存在上升流。采用比礁高評(píng)價(jià)法,上升流相對(duì)面積范圍0.74~46.24,AR06上升流相對(duì)面積最小,AR18最大(AR18的速度矢量圖如圖6a所示)。在相同礁高條件下,框架型礁體上升流面積均小于其他類型礁體,而復(fù)合型礁上升流面積最大,三角型礁次之。方型礁AR03、框架型礁AR07、棱臺(tái)型礁AR10、三角型礁AR11、圓柱型礁AR15和復(fù)合型礁AR18分別是同類型中上升流相對(duì)面積最大的礁體。采用比混凝土體積評(píng)價(jià)法,上升流相對(duì)面積范圍0.87~10.19,AR08上升流相對(duì)面積最小,AR13最大(AR13的速度矢量圖如圖6b所示)。相同混凝土體積條件下,框架型礁體上升流面積均小于其他類型礁體,方型礁AR02、框架型礁AR07、棱臺(tái)型礁AR10、三角型礁AR11、圓柱型礁AR13和復(fù)合型礁AR18分別是同類型中上升流相對(duì)面積最大的礁體。

    圖5 兩種上升流相對(duì)面積

    AR18由于迎流面寬度較大,用礁高比法分析其上升流相對(duì)面積明顯大于其它礁體;而用混凝土體積法進(jìn)行比較時(shí),差異明顯減小。所以在相同礁高條件下,混凝土體積的差異會(huì)對(duì)試驗(yàn)結(jié)果產(chǎn)生較大的影響,故認(rèn)為比混凝土體積法能更科學(xué)、準(zhǔn)確地評(píng)價(jià)礁體上升流效應(yīng)。

    綜上所述,比混凝土體積法得到的結(jié)果,框架型礁上升流相對(duì)面積最??;三角型礁上升流相對(duì)面積最大;復(fù)合型礁上升流相對(duì)面積僅次于三角型礁;方型礁、棱臺(tái)型礁和圓柱型礁的上升流相對(duì)面積差別不大。

    圖6是AR18、AR13和AR11在中垂面上的速度矢量圖??梢钥闯隽黧w經(jīng)過(guò)礁體,會(huì)在礁體上方產(chǎn)生上升流,后方產(chǎn)生背渦流,并且不同形狀的礁體上升流和背渦流的位置、規(guī)模等也會(huì)存在差異。

    圖6 礁體在中垂面上的速度矢量圖

    2.2 背渦流的流場(chǎng)效應(yīng)

    圖7是礁體在中垂面上的背渦流相對(duì)面積對(duì)比圖,可以看出18種礁體均存在背渦流。采用比礁高法,背渦流相對(duì)面積范圍0.19~31.80,AR07最小,AR18最大。在相同礁高條件下,框架型礁背渦流面積均小于其他類型礁體,而三角型礁、復(fù)合型礁均大于其他類型的礁體(除AR01)。方型礁AR01、框架型礁AR06、棱臺(tái)型礁AR10、三角型礁AR11、圓柱型礁AR14和復(fù)合型礁AR18分別是同類型中背渦流相對(duì)面積最大的礁體。采用比混凝土體積法,背渦流相對(duì)面積范圍0.27~10.52,AR08最小,AR11最大(AR11的速度矢量圖如圖6c所示)。相同混凝土體積條件下,框架型礁體背渦流面積均小于其他類型礁體(除AR09)。方型礁AR01、框架型礁AR06、棱臺(tái)型礁AR10、三角型礁AR11、圓柱型礁AR13和復(fù)合型礁AR16分別是同類型中背渦流相對(duì)面積最大的礁體。

    圖7 兩種背渦流相對(duì)面積

    通過(guò)比礁高和比混凝土體積分析不同礁體背渦流的流場(chǎng)效應(yīng),與上升流流場(chǎng)效應(yīng)同理,比較2種相對(duì)面積評(píng)價(jià)背渦流方法,比混凝土體積法仍是更科學(xué)準(zhǔn)確地評(píng)價(jià)背渦流面積方法。

    綜上所述,比混凝土體積法得到的結(jié)果,框架型礁背渦流相對(duì)面積最?。蝗切徒副硿u流相對(duì)面積最大;復(fù)合型礁背渦流相對(duì)面積僅次于三角型礁;方型礁、棱臺(tái)型礁和圓柱型礁的背渦流相對(duì)面積差別不大。

    2.3 流速對(duì)流場(chǎng)效應(yīng)的影響

    根據(jù)比混凝土體積法的結(jié)果,每種類型礁體以上升流與背渦流相對(duì)面積之和最大為標(biāo)準(zhǔn),選取的6個(gè)礁體為AR02、AR06、AR10、AR11、AR13和AR16。分析其在0.2、0.4、0.6、0.8和1.0 m/s 5種流速下上升流和背渦流的變化。

    如圖8所示,6個(gè)礁體在0.2、0.4、0.6、0.8和1.0m/s 5種流速下,上升流和背渦流面積雖然有波動(dòng)但變化很小,表明上升流和背渦流面積幾乎不隨來(lái)流速度的變化而變化。

    圖8 6個(gè)礁體在0.2、0.4、0.6、0.8和1.0m/s 5種來(lái)流速度下上升流和背渦流面積

    通過(guò)分析礁體中垂面上方向速度分量分布,得到礁體在不同來(lái)流速度下的最大上升流流速,最大上升流流速是在上升流區(qū)域中某一點(diǎn)處的最大值,相同礁體的最大上升流流速在相同點(diǎn)處取得,這與何文榮[29]的結(jié)論是一致的。圖9給出了6種礁體在不同來(lái)流速度下最大上升流速度的變化曲線。從圖中可以看出,所有礁體最大上升流流速與來(lái)流速度均為線性關(guān)系,擬合曲線2的范圍為0.9997~1。2越接近1,擬合度越高,即說(shuō)明礁體的最大上升流流速與來(lái)流速度的比幾乎不隨來(lái)流速度的變化而變化,此結(jié)論與何文榮[29]的金字塔型人工魚礁CFD研究中所得結(jié)論是一致的。從圖中進(jìn)一步分析得出不同礁體的最大上升流流速與來(lái)流速度比值存在差異,比值由大到小分別為:AR16>AR10>AR13>AR02>AR11> AR06。通過(guò)最大上升流流速與來(lái)流速度的比值排序可以看出,多數(shù)迎流面傾斜的礁體在相同的來(lái)流速度條件下,最大上升流流速更大。AR13由三個(gè)混凝土圓管組成,其迎流面也相當(dāng)于具有一定的傾角。AR11的比值比AR02小,可能是由于AR11迎流面上開口較多,影響了傾斜的迎流板該有的上升流效果。

    圖9 最大上升流流速與來(lái)流速度的關(guān)系

    圖10a為上升流高度比礁高與來(lái)流速度的關(guān)系,其上升流高度幾乎不隨來(lái)流速度變化,這與黃遠(yuǎn)東、劉彥的研究結(jié)論一致[30-31]。上升流高度比礁高由大到小依次為:AR11>AR16>AR02>AR10>AR13>AR06,礁高越大,其單位礁高的上升流高度越大(AR06除外);圖10b是上升流高度比混凝土體積與來(lái)流速度的關(guān)系,上升流高度比混凝土體積由大到小依次為:AR13>AR06>AR10> AR11>AR02>AR16,混凝土體積越小,其單位混凝土體積的上升流高度越大(AR10和AR11除外)。在研究上升流高度時(shí),比混凝土體積法傾向于混凝土體積較小的礁體,而比礁高法較為準(zhǔn)確的表征出礁體上升流高度。

    3 結(jié) 論

    本文運(yùn)用CFX對(duì)18種中國(guó)常用的人工魚礁進(jìn)行模擬,利用比礁高和比混凝土體積2種方法對(duì)中垂面上升流面積、背渦流面積和上升流高度進(jìn)行了分析,得出以下結(jié)論:

    1)將18種礁體分為方型礁、框架型礁、棱臺(tái)型礁、三角型礁、圓柱型礁和復(fù)合型礁等6大類,礁體的上升流和背渦流相對(duì)面積采用較好的比混凝土體積法評(píng)價(jià)。結(jié)果表明,無(wú)論上升流和背渦流相對(duì)面積,三角型礁為最好,復(fù)合型礁次之,而框架型礁為最差。方型礁、棱臺(tái)型礁和圓柱型礁相對(duì)面積相近且效果中等。

    2)選出的6個(gè)代表礁體的上升流面積和背渦流面積大小與來(lái)流速度不相關(guān);最大上升流流速與來(lái)流速度呈線性關(guān)系,其斜率在不同礁體間存在差異;礁體的上升流高度與來(lái)流速度不相關(guān)。

    3)18種礁體上升流和背渦流相對(duì)面積采用比混凝土體積法評(píng)價(jià),上升流高度采用比礁高法評(píng)價(jià)。在保證低成本的情況下,要得到更大的背渦流面積,應(yīng)選擇AR11、AR10和AR13;要得到更大的上升流面積,應(yīng)選擇AR13、AR11和AR18。AR11有著最大的背渦流相對(duì)面積和最大的相對(duì)影響面積(上升流相對(duì)面積與背渦流相對(duì)面積之和),而且在相同礁高的情況下,也有最大的上升流高度(其次是AR16和AR02),但其最大上升流流速較小。最大上升流流速最大的是AR16,其次是AR10和AR13。

    通過(guò)比礁高和比混凝土體積2種方法可以分析研究不同類型礁體的的上升流和背渦流相關(guān)特性,為不同海域、不同條件下人工魚礁礁型的選擇提供了一定的理論指導(dǎo),同時(shí)在人工魚礁設(shè)計(jì)上也有一定的參考價(jià)值。但相關(guān)參數(shù)隨礁高或混凝土體積的變化規(guī)律尚不明確,需今后進(jìn)一步研究。

    [1] 陳勇,于長(zhǎng)清,張國(guó)勝,等. 人工魚礁的環(huán)境功能與集魚效果[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào),2002,17(1):64-69. Chen Yong, Yu Changqing, Zhang Guosheng, et al. The environmental function and fish gather effect of artificial reefs[J]. Journal of Dalian Ocean University, 2002, 17(1): 64-69. (in Chinese with English abstract)

    [2] Branden K L, Pollard D A, Reimers H A. A review of recent artificial reef developments in Australia[J]. Bulletin of Marine Science -Miami-, 1994, 55(2/3): 982-994.

    [3] Clark S, Edwards A J. Use of artificial reef structures to rehabilitate reef flats degraded by coral mining in the maldives[J]. Bulletin of Marine Science -Miami-, 1994, 55 (2/3): 724-744.

    [4] Baine M. Artificial reefs: A review of their design, application,management and performance[J]. Ocean & Coastal Management, 2001, 44(3/4): 241-259.

    [5] 劉惠飛. 日本人工魚礁建設(shè)的現(xiàn)狀[J]. 現(xiàn)代漁業(yè)信息,2001(12):15-17. Liu Huifei. The status of construction of artificial fish reef in Japan[J]. Journal of Modern Fisheries Information, 2001(12): 15-17. (in Chinese with English abstract)

    [6] 王波,武建平,高峻,等. 關(guān)于青島建設(shè)人工魚礁改善近海生態(tài)和漁業(yè)環(huán)境的探討[J]. 海岸工程,2004(4):66-73. Wang Bo, Wu Jianping, Gao Jun, et al. An approach to construction of artificial reef to improve the offshore ecology and fishery environments in Qingdao[J]. Coastal Engineering, 2004(4): 66-73. (in Chinese with English abstract)

    [7] 張碩,孫滿昌,陳勇. 不同高度混凝土模型礁上升流特性的定量研究[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào),2008,23(5):353-358. Zhang Shuo, Sun Manchang, Chen Yong. Quantitative analysis of upwelling current features of a artificial concrete reef with different height[J]. Journal of Dalian Ocean University, 2008, 23(5): 353-358. (in Chinese with English abstract)

    [8] 賈曉平,陳丕茂,唐振朝,等. 人工魚礁關(guān)鍵技術(shù)研究與示范[M]. 北京: 海洋出版社,2011:30-33.

    [9] 上北征男. 水産増養(yǎng)殖施設(shè)の水力學(xué)的研究[R]. 水産工學(xué)研究所報(bào)告,1982(3):67-100. Yukio Uekita. Hydraulic studies on the aquacutural facilities[R]. Bulletin of National Research Institute of Fisheries Engineering, 1982(3): 67-100. (in Japanese with English abstract)

    [10] 秀島好昭,上北征男. 人工魚礁の著底衝擊力に関する研究[R]. 水産工學(xué)研究所報(bào)告,1983(4):59-71. Yoshiaki Hideshima, Yukio Uekita. Study on the impact force to the artificial fish scheme colliding with the sea bottom[R]. Bulletin of National Research Institute of Fisheries Engineering, 1983(4): 59-71. (in Japanese with English abstract)

    [11] 上北征男,中村充,秀島好昭. 海底礁による地形性湧昇流に関する研究[R]. 水産工學(xué)研究所報(bào)告,1984(5): 33-66. Yukio Uekita, Makoto Nakamura, Yoshiaki Hideshima. Observation of the upwelling causing behind the large reef in case of siomaki, Yamaguchi Pref[R]. Bulletin of National Research Institute of Fisheries Engineering, 1984(5): 33-66. (in Japanese with English abstract)

    [12] 鄭延璇,梁振林,關(guān)長(zhǎng)濤,等. 三種疊放形式的圓管型人工魚礁流場(chǎng)效應(yīng)數(shù)值模擬與PIV試驗(yàn)研究[J]. 海洋與湖沼,2014,45(1):11-19. Zheng Yanxuan, Liang Zhenlin, Guan Changtao, et al. Numerical simulation and experimental study on flow field of artificial reefs in three tube-stacking layouts[J]. Oceanologia et Limnologia Sinica, 2014, 45(1): 11-19. (in Chinese with English abstract)

    [13] 劉洪生,馬翔,章守宇,等. 人工魚礁流場(chǎng)風(fēng)洞實(shí)驗(yàn)與數(shù)值模擬對(duì)比驗(yàn)證[J]. 中國(guó)水產(chǎn)科學(xué),2009(3):365-371. Liu Hongsheng, Ma Xiang, Zhang Shouyu, et al. Validation and comparison between wind tunnel experiments and numerical simulation of flow field around artificial reefs[J]. Journal of Fishery Sciences of China, 2009(3): 365-371. (in Chinese with English abstract)

    [14] 崔勇,關(guān)長(zhǎng)濤,萬(wàn)榮,等. 布設(shè)間距對(duì)人工魚礁流場(chǎng)效應(yīng)影響的數(shù)值模擬[J]. 海洋湖沼通報(bào),2011(2):59-65. Cui Yong, Guan Changtao, Wan Rong, et al. Numerical simulation on influence of disposal space for effects of flow field around artificial reefs[J]. Transactions of Oceanology and Limnology, 2011(2): 59-65. (in Chinese with English abstract)

    [15] 劉洪生,馬翔,章守宇,等. 人工魚礁流場(chǎng)效應(yīng)的模型實(shí)驗(yàn)[J]. 水產(chǎn)學(xué)報(bào),2009(2):229-236. Liu Hongsheng, Ma Xiang, Zhang Shouyu, et al. Research on model experiments of effects of artificial reefs on flow field[J]. Journal of Fisheries of China, 2009(2): 229-236. (in Chinese with English abstract)

    [16] Kim K, Ahmed M R, Lee Y. Efficiency improvement of a tidal current turbine utilizing a larger area of channel[J]. Renewable Energy, 2012(48): 557-564.

    [17] Finnegan W, Goggins J. Numerical simulation of linear water waves and wave–structure interaction[J]. Ocean Engineering, 2012(43): 23-31.

    [18] Kim D, Woo J, Yoon H, et al. Wake lengths and structural responses of Korean general artificial reefs[J]. Ocean Engineering, 2014(92): 83-91.

    [19] Woo J, Kim D, Yoon H, et al. Characterizing Korean general artificial reefs by drag coefficients[J]. Ocean Engineering, 2014(82): 105-114.

    [20] 王福軍. 計(jì)算流體力學(xué)分析[M]. 北京:清華大學(xué)出版社,2004.

    [21] Yakhot V, Orzag S A. Renormalization group analysis of turbulence: basic theory[J]. Journal of Scientific Computation, 1986(1): 3-11.

    [22] Veersteg H K, Malalasekera W. An introduction to computational fluid dynamics: The finite volume method[J]. Pearson Schweiz Ag, 1995, 20(5): 400.

    [23] 郭鴻志. 傳輸過(guò)程數(shù)值模擬[M]. 北京:冶金工業(yè)出版社,1998.

    [24] 陶文銓. 數(shù)值傳熱學(xué)(第2版)[M]. 西安:西安交通大學(xué)出版社,2001.

    [25] 王磊,唐衍力,黃洪亮,等. 混凝土人工魚礁選型的初步分析[J]. 海洋漁業(yè),2009(3):308-315. Wang Lei, Tang Yanli, Huang Hongliang, et al. Primary analysis of model selection for artificial reefs[J]. Marine Fisheries, 2009(3): 308-315. (in Chinese with English abstract)

    [26] 鄭延璇,關(guān)長(zhǎng)濤,宋協(xié)法,等. 星體型人工魚礁流場(chǎng)效應(yīng)的數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(19):185-193. Zheng Yanxuan, Guan Changtao, Song Xiefa, et al. Numerical simulation on flow field around star artificial reefs[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(19): 185-193. (in Chinese with English abstract)

    [27] 關(guān)長(zhǎng)濤,劉彥,趙云鵬,等. 復(fù)合M型人工魚礁粒子圖像測(cè)速二維流場(chǎng)試驗(yàn)研究[J]. 漁業(yè)現(xiàn)代化,2010,37(1): 15-19. Guan Changtao, Liu Yan, Zhao Yunpeng, et al. Experimental study on two-dimensional flow field of the compound M artificial reef with particle image velocimetry[J]. Fishery Modernization, 2010, 37(1): 15-19. (in Chinese with English abstract)

    [28] 謝亮. Photoshop像素法在計(jì)算地圖面積中的應(yīng)用[J]. 電腦知識(shí)與技術(shù),2010,6(15):4021-4022.Xie Liang. Application of Photoshop pixel method in calculating map area[J]. Computer Knowledge and Technology, 2010, 6(15): 4021-4022. (in Chinese with English abstract)

    [29] 何文榮,黃遠(yuǎn)東,黃黎明,等. 金字塔型人工魚礁繞流的三維CFD模擬研究[J]. 水資源與水工程學(xué)報(bào),2013,24(5):71–76. He Wenrong, Huang Yuandong, Huang Liming, et al. Simulation of three-dimensional CFD of water flowat pyramid artificial reef[J].Journal of Water Resources & Water Engineering, 2013, 24(5): 71–76. (in Chinese with English abstract)

    [30] 黃遠(yuǎn)東,姜?jiǎng)?,趙樹夫. 方型人工魚礁周圍水流運(yùn)動(dòng)的數(shù)值模擬研究[J]. 水資源與水工程學(xué)報(bào),2012(3):1-3. Huang Yuandong, Jiang Jianwei, Zhao Shufu. Study on numerical model of water flows past a square artificial reef[J]. Journal of Water Resources and Water Engineering, 2012(3): 1-3. (in Chinese with English abstract)

    [31] Liu Y, Zhao Y P, Dong G H, et al. A study of the flow field characteristics around star-shaped artificial reefs[J]. Journal of Fluids & Structures, 2013, 39(5): 27-40.

    Comparative analysis on flow field effect of general artificial reefs in China

    Tang Yanli1, Long Xiangyu1, Wang Xinxin1※, Jiang Zhaoyang2, Cheng Hui1, Zhang Tongzheng1

    (1.,266003,; 2.,,264209,)

    Artificial reef (AR) is widely used to improve marine ecological environment. Some fisheries developed countries, such as America, Japan, Korea, Australia, New Zealand and some EU (European Union) states, have acquired a lot of achievements on the exploitation and protection of fishery resources by use of ARs. Flow field effect, bait effect and avoidance effect will be caused when the ARs are placed on the sea floor. Flow field effect is deemed to the main affecting mechanism of ARs, and the upwelling and back vortex are important indices to measure the flow field effect of ARs. So far, the study on the flow field effect of ARs has been limited to one or several reefs in China. However, no one has yet scientifically investigated the flow field effect of the Chinese general ARs. This study focused on the 18 types of general ARs, which were divided into 6 kinds, to investigate different flow field effects among them. And several representative reefs were selected for further study, which was the study on the effects of various flow velocity on the characteristics of upwelling and back vortex. Based on theory of computational fluid dynamics, numerical models of turbulent flow were built. Three-dimensional Navier-Stokes equations were solved by the finite volume method. The software CFX (computational fluid X) was used to study the performances of the flow field around reefs when flow passed through them. The 2 evaluation methods, the reef height ratio and the concrete volume ratio, were used to analyze the upwelling area, the back vortex area and the height of upwelling on the central plane. The reef height ratio is the ratio of the maximum upwelling height, upwelling area or back vortex area to reef height. And the concrete volume ratio is the ratio of the maximum upwelling height, upwelling area or back vortex area to concrete volume of reef. The paper was aimed to discuss the advantages and disadvantages of the 2 methods, and to analyze the different flow field effects among 18 types of ARs. Besides, the better reefs would be selected for diverse requirements. The results suggested the concrete volume ratio was a better criterion when studying the upwelling area and the back vortex area, and the height ratio was more appropriate when studying the height of upwelling. Regardless of the upwelling or the back vortex, triangle reefs had the largest relative area and complex reefs were the second, while frame reefs were the least. Then one reef was selected from each type to conduct the research to find out whether there were the effects of different flow velocity (0.2, 0.4, 0.6, 0.8 and 1.0 m/s) on upwelling area, back vortex area, the maximum upwelling velocity and maximum height of upwelling. The results demonstrated that the height of upwelling, the upwelling area, and the back vortex area fluctuated in an extremely narrow range under the 5 flow velocities. The maximum flow velocity of upwelling was linear with the flow velocity, and the slope showed the difference among different reefs. Moreover, reefs for different demands had been selected. AR11, AR10 and AR13 had a larger back vortex area with low cost. AR13, AR11 and AR18 should be chosen for the need of larger upwelling area. AR11 had not only the largest relative area of back vortex but also the largest relative influence area (sum of the upwelling area and the back vortex area). In addition, AR11 had the highest height of upwelling at the same reef height (followed by AR16 and AR02). The top 3 types of reefs with the maximum flow velocity were AR16, AR10 and AR13. The study provides the theoretical basis for the selection of AR, which can cater to various demands. Moreover, it also has certain reference value for the design of AR.

    reefs; flow fields; numerical simulation; artificial; the height ratio; the concrete volume ratio

    10.11975/j.issn.1002-6819.2017.08.013

    S96

    A

    1002-6819(2017)-08-0097-07

    2016-08-23

    2017-04-20

    海洋公益性專項(xiàng)(201305030);公益性行業(yè)(農(nóng)業(yè))科技專項(xiàng)(201203018)

    唐衍力,男,教授,博士,研究方向?yàn)槿斯~礁與海洋牧場(chǎng)、選擇性漁具漁法。青島 中國(guó)海洋大學(xué)水產(chǎn)學(xué)院,266003。 Email:tangyanli@ouc.edu.cn

    王欣欣,女,講師,博士,研究方向?yàn)闈O具設(shè)施水動(dòng)力學(xué)研究。青島 中國(guó)海洋大學(xué)水產(chǎn)學(xué)院,266003。Email:wxinxin@ouc.edu.cn

    唐衍力,龍翔宇,王欣欣,姜昭陽(yáng),程 暉,張同征.中國(guó)常用人工魚礁流場(chǎng)效應(yīng)的比較分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(8):97-103. doi:10.11975/j.issn.1002-6819.2017.08.013 http://www.tcsae.org

    Tang Yanli, Long Xiangyu, Wang Xinxin, Jiang Zhaoyang, Cheng Hui, Zhang Tongzheng. Comparative analysis on flow field effect of general artificial reefs in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(8): 97-103. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.08.013 http://www.tcsae.org

    猜你喜歡
    上升流礁體場(chǎng)效應(yīng)
    夏季舟山上升流長(zhǎng)期遙感觀測(cè)與分析
    生態(tài)礁體布置形式對(duì)水流特征的影響
    南麂列島海域兩種框型人工魚礁水動(dòng)力性能試驗(yàn)
    遼東半島頂端海域上升流長(zhǎng)期變化特征及影響因素*
    海洋與湖沼(2020年1期)2020-02-06 06:01:38
    場(chǎng)效應(yīng)晶體管短路失效的數(shù)值模型
    珠海萬(wàn)山人工魚礁結(jié)構(gòu)設(shè)計(jì)
    海底新生代礁體下伏構(gòu)造的校正淺析
    海洋石油(2017年1期)2017-03-29 01:37:25
    基于CH3NH3PbI3單晶的Ta2O5頂柵雙極性場(chǎng)效應(yīng)晶體管
    瓊東上升流的年際變化及長(zhǎng)期變化趨勢(shì)*
    海洋與湖沼(2016年1期)2016-01-15 03:50:49
    建筑學(xué)專業(yè)設(shè)計(jì)系列課程“場(chǎng)效應(yīng)”教學(xué)模式探索與實(shí)踐
    国产在线精品亚洲第一网站| 99久久人妻综合| 国产精品成人在线| 老汉色av国产亚洲站长工具| 一区二区三区激情视频| 黑人欧美特级aaaaaa片| 99久久国产精品久久久| 欧美日韩亚洲综合一区二区三区_| avwww免费| 日本黄色视频三级网站网址 | 91老司机精品| 大片免费播放器 马上看| 国产精品1区2区在线观看. | av国产精品久久久久影院| 精品第一国产精品| 熟女少妇亚洲综合色aaa.| 亚洲,欧美精品.| 蜜桃国产av成人99| 欧美国产精品一级二级三级| 免费观看人在逋| 美女午夜性视频免费| 欧美成人午夜精品| 午夜福利视频精品| 成人三级做爰电影| 丝袜人妻中文字幕| 视频区欧美日本亚洲| 99re在线观看精品视频| 日本黄色视频三级网站网址 | 精品国内亚洲2022精品成人 | 在线看a的网站| 在线观看免费日韩欧美大片| 国产精品免费一区二区三区在线 | 精品国产国语对白av| 一区二区av电影网| 日韩熟女老妇一区二区性免费视频| 狠狠精品人妻久久久久久综合| 亚洲成人手机| 亚洲avbb在线观看| 后天国语完整版免费观看| 欧美激情高清一区二区三区| 国产91精品成人一区二区三区 | 国产真人三级小视频在线观看| 丰满少妇做爰视频| 脱女人内裤的视频| 淫妇啪啪啪对白视频| a在线观看视频网站| 欧美 日韩 精品 国产| 日韩一区二区三区影片| 国产视频一区二区在线看| 精品国产亚洲在线| 一级黄色大片毛片| 免费日韩欧美在线观看| 中文字幕最新亚洲高清| 老鸭窝网址在线观看| 欧美乱妇无乱码| 国产97色在线日韩免费| 搡老乐熟女国产| 啦啦啦视频在线资源免费观看| 国产精品免费大片| 热99re8久久精品国产| 色精品久久人妻99蜜桃| 黑丝袜美女国产一区| 国产不卡av网站在线观看| 一个人免费看片子| 老熟妇仑乱视频hdxx| 欧美人与性动交α欧美软件| netflix在线观看网站| 2018国产大陆天天弄谢| 18禁裸乳无遮挡动漫免费视频| 午夜激情久久久久久久| 亚洲国产中文字幕在线视频| 国产高清国产精品国产三级| 2018国产大陆天天弄谢| 免费av中文字幕在线| 日本一区二区免费在线视频| 中亚洲国语对白在线视频| 大香蕉久久网| 日韩欧美一区视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 色老头精品视频在线观看| 亚洲色图综合在线观看| 国产男女内射视频| 别揉我奶头~嗯~啊~动态视频| 99精国产麻豆久久婷婷| 在线天堂中文资源库| 一进一出好大好爽视频| 亚洲国产看品久久| av线在线观看网站| 亚洲国产欧美日韩在线播放| 麻豆成人av在线观看| 中文字幕高清在线视频| 狂野欧美激情性xxxx| 亚洲国产av影院在线观看| 久9热在线精品视频| 又大又爽又粗| 交换朋友夫妻互换小说| 久久久国产精品麻豆| 宅男免费午夜| h视频一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 伊人久久大香线蕉亚洲五| 欧美亚洲日本最大视频资源| 国产亚洲av高清不卡| 97人妻天天添夜夜摸| 免费黄频网站在线观看国产| 欧美精品亚洲一区二区| 国产精品熟女久久久久浪| 在线观看舔阴道视频| 精品高清国产在线一区| 窝窝影院91人妻| 青青草视频在线视频观看| 高清av免费在线| 狂野欧美激情性xxxx| 成人国产av品久久久| 国产三级黄色录像| 大香蕉久久网| 天天操日日干夜夜撸| 满18在线观看网站| 巨乳人妻的诱惑在线观看| 国产精品二区激情视频| 免费观看人在逋| 欧美精品人与动牲交sv欧美| 免费人妻精品一区二区三区视频| bbb黄色大片| 欧美一级毛片孕妇| 亚洲精品av麻豆狂野| 亚洲专区国产一区二区| 在线观看一区二区三区激情| 国产欧美亚洲国产| 99riav亚洲国产免费| 久久久久久亚洲精品国产蜜桃av| 精品少妇久久久久久888优播| 97人妻天天添夜夜摸| 涩涩av久久男人的天堂| 啦啦啦 在线观看视频| 这个男人来自地球电影免费观看| e午夜精品久久久久久久| 91老司机精品| 亚洲欧美日韩另类电影网站| 国产亚洲欧美在线一区二区| 国产成人啪精品午夜网站| 日日夜夜操网爽| 12—13女人毛片做爰片一| 国产精品免费一区二区三区在线 | 久久人妻福利社区极品人妻图片| 久久人妻av系列| 天天躁夜夜躁狠狠躁躁| 亚洲精华国产精华精| 精品国产乱子伦一区二区三区| 乱人伦中国视频| 一夜夜www| 成人影院久久| 飞空精品影院首页| 日韩欧美免费精品| 2018国产大陆天天弄谢| 在线亚洲精品国产二区图片欧美| 69av精品久久久久久 | 热99久久久久精品小说推荐| 国产欧美日韩一区二区精品| 色播在线永久视频| 久久天堂一区二区三区四区| e午夜精品久久久久久久| 国产精品久久久久久精品电影小说| 人妻一区二区av| 精品一区二区三区四区五区乱码| 下体分泌物呈黄色| 免费高清在线观看日韩| 久久人人爽av亚洲精品天堂| 国产xxxxx性猛交| 中国美女看黄片| 两个人看的免费小视频| 亚洲国产毛片av蜜桃av| 久久香蕉激情| 在线观看66精品国产| 欧美激情极品国产一区二区三区| 中文字幕色久视频| 18禁观看日本| a在线观看视频网站| 色播在线永久视频| 国产免费av片在线观看野外av| 中文字幕人妻熟女乱码| 又紧又爽又黄一区二区| 俄罗斯特黄特色一大片| 日本av手机在线免费观看| 午夜福利视频在线观看免费| 69精品国产乱码久久久| 免费在线观看视频国产中文字幕亚洲| 久久中文看片网| 极品人妻少妇av视频| 精品福利观看| 亚洲性夜色夜夜综合| 国产精品免费大片| 自线自在国产av| 亚洲av欧美aⅴ国产| av国产精品久久久久影院| 十八禁高潮呻吟视频| 国产精品av久久久久免费| 97在线人人人人妻| 免费在线观看完整版高清| 超色免费av| 夫妻午夜视频| 国产精品久久久久久精品电影小说| 男人舔女人的私密视频| 五月开心婷婷网| 久久精品91无色码中文字幕| 亚洲精品乱久久久久久| 久久精品亚洲精品国产色婷小说| 精品亚洲成a人片在线观看| 免费高清在线观看日韩| 视频区图区小说| 热re99久久国产66热| 丝瓜视频免费看黄片| 国产精品一区二区在线观看99| 男人操女人黄网站| 脱女人内裤的视频| 三上悠亚av全集在线观看| 国产亚洲欧美在线一区二区| 亚洲国产av影院在线观看| 欧美乱码精品一区二区三区| 无遮挡黄片免费观看| 亚洲精品中文字幕一二三四区 | 国产高清视频在线播放一区| 天天操日日干夜夜撸| 黄色毛片三级朝国网站| 18禁国产床啪视频网站| 亚洲精品乱久久久久久| 飞空精品影院首页| 桃红色精品国产亚洲av| 激情在线观看视频在线高清 | bbb黄色大片| 在线观看66精品国产| 久久中文字幕人妻熟女| 日韩精品免费视频一区二区三区| 99精品久久久久人妻精品| 老司机在亚洲福利影院| 18禁美女被吸乳视频| 久久99热这里只频精品6学生| 国产又色又爽无遮挡免费看| 99热国产这里只有精品6| 又黄又粗又硬又大视频| 亚洲第一av免费看| 亚洲一区二区三区欧美精品| 婷婷丁香在线五月| 性少妇av在线| 欧美日韩国产mv在线观看视频| 国产亚洲av高清不卡| 在线观看66精品国产| 国产三级黄色录像| av网站在线播放免费| 99久久99久久久精品蜜桃| 精品一品国产午夜福利视频| 久久天躁狠狠躁夜夜2o2o| 久久国产亚洲av麻豆专区| 巨乳人妻的诱惑在线观看| 亚洲熟女精品中文字幕| 亚洲午夜理论影院| 色婷婷av一区二区三区视频| 99国产精品一区二区三区| 国产xxxxx性猛交| 中国美女看黄片| 国产成人av激情在线播放| 18禁观看日本| 日韩欧美免费精品| 亚洲人成电影免费在线| 一边摸一边抽搐一进一出视频| 捣出白浆h1v1| 免费不卡黄色视频| 欧美国产精品va在线观看不卡| 午夜福利欧美成人| 国产精品 国内视频| 50天的宝宝边吃奶边哭怎么回事| 中文字幕制服av| 中文字幕另类日韩欧美亚洲嫩草| 熟女少妇亚洲综合色aaa.| 亚洲欧美激情在线| 国产在线免费精品| 操出白浆在线播放| 国产在视频线精品| 一区福利在线观看| 男女之事视频高清在线观看| 极品人妻少妇av视频| av不卡在线播放| videos熟女内射| 一级,二级,三级黄色视频| 国产男女内射视频| 久久婷婷成人综合色麻豆| 国产午夜精品久久久久久| 性少妇av在线| 国产老妇伦熟女老妇高清| 日韩欧美三级三区| 99热国产这里只有精品6| 无限看片的www在线观看| 国产又色又爽无遮挡免费看| 免费人妻精品一区二区三区视频| 黄色毛片三级朝国网站| 久久狼人影院| 高清毛片免费观看视频网站 | 男女无遮挡免费网站观看| 男人操女人黄网站| 国产精品电影一区二区三区 | 精品乱码久久久久久99久播| 在线观看免费午夜福利视频| 在线十欧美十亚洲十日本专区| 国产精品免费一区二区三区在线 | 视频在线观看一区二区三区| 窝窝影院91人妻| 国产又色又爽无遮挡免费看| 欧美日韩视频精品一区| 一级毛片精品| 日韩 欧美 亚洲 中文字幕| 又大又爽又粗| 亚洲av片天天在线观看| 一区二区日韩欧美中文字幕| 这个男人来自地球电影免费观看| 国产亚洲精品一区二区www | 久久毛片免费看一区二区三区| 中文字幕人妻熟女乱码| 美国免费a级毛片| 国产免费av片在线观看野外av| 亚洲熟女精品中文字幕| 天天躁日日躁夜夜躁夜夜| 欧美变态另类bdsm刘玥| 99香蕉大伊视频| 五月天丁香电影| 老汉色av国产亚洲站长工具| 国产成人欧美在线观看 | 不卡一级毛片| 少妇裸体淫交视频免费看高清 | 大片免费播放器 马上看| 国产淫语在线视频| 婷婷丁香在线五月| 日韩 欧美 亚洲 中文字幕| 在线av久久热| 国产精品麻豆人妻色哟哟久久| av片东京热男人的天堂| 丝袜喷水一区| 日韩中文字幕视频在线看片| 搡老岳熟女国产| 欧美黄色淫秽网站| 丁香欧美五月| 熟女少妇亚洲综合色aaa.| 久久国产精品男人的天堂亚洲| 亚洲免费av在线视频| 国产深夜福利视频在线观看| 亚洲国产av影院在线观看| 嫩草影视91久久| 操出白浆在线播放| 韩国精品一区二区三区| 99热网站在线观看| 精品亚洲成a人片在线观看| aaaaa片日本免费| 人人妻,人人澡人人爽秒播| 搡老岳熟女国产| 午夜视频精品福利| 国产欧美日韩精品亚洲av| 桃花免费在线播放| 侵犯人妻中文字幕一二三四区| 一个人免费在线观看的高清视频| av一本久久久久| 精品亚洲成国产av| 欧美日韩成人在线一区二区| 精品久久久精品久久久| 国产精品国产av在线观看| 人人妻人人澡人人看| 777久久人妻少妇嫩草av网站| 亚洲av电影在线进入| 成人18禁高潮啪啪吃奶动态图| 国产精品98久久久久久宅男小说| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 国产在线一区二区三区精| 男女下面插进去视频免费观看| 成人影院久久| 99久久人妻综合| 电影成人av| 青青草视频在线视频观看| 十八禁高潮呻吟视频| 久久精品成人免费网站| 国产成人精品久久二区二区免费| 精品国产一区二区三区四区第35| 99riav亚洲国产免费| 亚洲专区字幕在线| 精品少妇黑人巨大在线播放| 色婷婷av一区二区三区视频| 久久久精品区二区三区| 亚洲av日韩在线播放| 啦啦啦在线免费观看视频4| 午夜福利乱码中文字幕| 亚洲国产av新网站| 亚洲熟女毛片儿| 纵有疾风起免费观看全集完整版| 美女福利国产在线| av线在线观看网站| 热99re8久久精品国产| 国产在线一区二区三区精| 夜夜夜夜夜久久久久| 亚洲成人免费av在线播放| 大码成人一级视频| 国产精品国产av在线观看| 久久影院123| 80岁老熟妇乱子伦牲交| 91老司机精品| 90打野战视频偷拍视频| 亚洲va日本ⅴa欧美va伊人久久| 在线看a的网站| 欧美精品亚洲一区二区| 欧美精品一区二区免费开放| 大型黄色视频在线免费观看| 一区二区三区精品91| 岛国在线观看网站| 他把我摸到了高潮在线观看 | 亚洲人成77777在线视频| 免费在线观看黄色视频的| 97人妻天天添夜夜摸| 极品教师在线免费播放| 一级毛片女人18水好多| 一区福利在线观看| 桃花免费在线播放| 久久人人97超碰香蕉20202| 嫁个100分男人电影在线观看| 在线 av 中文字幕| 国产男女超爽视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 少妇被粗大的猛进出69影院| 可以免费在线观看a视频的电影网站| 狠狠精品人妻久久久久久综合| 一区二区三区精品91| 水蜜桃什么品种好| 男女免费视频国产| 成人永久免费在线观看视频 | 精品高清国产在线一区| 久久久久精品人妻al黑| 最近最新中文字幕大全电影3 | 国产成人啪精品午夜网站| 蜜桃在线观看..| 久久99热这里只频精品6学生| 国产成人精品在线电影| 五月天丁香电影| 丝袜人妻中文字幕| 亚洲 国产 在线| 亚洲成国产人片在线观看| 亚洲精品久久成人aⅴ小说| 最近最新中文字幕大全电影3 | 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线观看免费 | 亚洲av成人一区二区三| 国产成人精品久久二区二区91| 亚洲专区国产一区二区| 亚洲人成电影免费在线| 搡老岳熟女国产| 亚洲精品国产色婷婷电影| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 19禁男女啪啪无遮挡网站| 一区二区三区乱码不卡18| cao死你这个sao货| 精品国产超薄肉色丝袜足j| 国产欧美日韩综合在线一区二区| 精品少妇内射三级| 亚洲精品一二三| 窝窝影院91人妻| 久久人妻福利社区极品人妻图片| 亚洲专区国产一区二区| 操出白浆在线播放| 欧美在线黄色| 我的亚洲天堂| 国产有黄有色有爽视频| 黑人巨大精品欧美一区二区mp4| 后天国语完整版免费观看| 十八禁高潮呻吟视频| 亚洲国产中文字幕在线视频| 捣出白浆h1v1| 一本久久精品| 十八禁人妻一区二区| 中文字幕制服av| 国产97色在线日韩免费| 国产欧美日韩一区二区三| 国产av国产精品国产| 成在线人永久免费视频| 黄色视频在线播放观看不卡| 亚洲色图 男人天堂 中文字幕| 精品人妻1区二区| 免费久久久久久久精品成人欧美视频| 精品人妻熟女毛片av久久网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲精华国产精华精| 女人被躁到高潮嗷嗷叫费观| 国产精品亚洲av一区麻豆| 久久天堂一区二区三区四区| 亚洲成人手机| 最近最新免费中文字幕在线| 国产高清视频在线播放一区| 人妻久久中文字幕网| 人人妻人人爽人人添夜夜欢视频| 亚洲午夜理论影院| 亚洲一区二区三区欧美精品| 黑人巨大精品欧美一区二区蜜桃| 最新在线观看一区二区三区| 欧美老熟妇乱子伦牲交| 午夜91福利影院| 亚洲精品自拍成人| av线在线观看网站| 国产片内射在线| 国产精品.久久久| 国产精品国产高清国产av | 人人澡人人妻人| 亚洲av日韩精品久久久久久密| 无限看片的www在线观看| 一级,二级,三级黄色视频| 悠悠久久av| 夜夜夜夜夜久久久久| 热99久久久久精品小说推荐| 精品卡一卡二卡四卡免费| 女人爽到高潮嗷嗷叫在线视频| 免费看十八禁软件| 精品人妻1区二区| 2018国产大陆天天弄谢| 亚洲第一青青草原| 国精品久久久久久国模美| 精品国产一区二区三区久久久樱花| 一级黄色大片毛片| 亚洲伊人久久精品综合| 91麻豆精品激情在线观看国产 | 高清在线国产一区| 91精品三级在线观看| 99国产精品一区二区三区| 日本av免费视频播放| 日韩成人在线观看一区二区三区| 国产精品一区二区免费欧美| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久精品古装| 亚洲精品乱久久久久久| 99国产精品免费福利视频| 日韩免费高清中文字幕av| 桃红色精品国产亚洲av| 亚洲美女黄片视频| 国产高清国产精品国产三级| 亚洲九九香蕉| 另类精品久久| 色综合欧美亚洲国产小说| 啦啦啦 在线观看视频| 曰老女人黄片| 一区二区av电影网| 国产不卡av网站在线观看| 日日摸夜夜添夜夜添小说| videosex国产| 又黄又粗又硬又大视频| 97人妻天天添夜夜摸| 黄频高清免费视频| 成年人午夜在线观看视频| 电影成人av| 日韩一区二区三区影片| 波多野结衣一区麻豆| av线在线观看网站| 日韩中文字幕视频在线看片| 一区福利在线观看| 精品第一国产精品| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 国产一区二区三区综合在线观看| 亚洲天堂av无毛| 12—13女人毛片做爰片一| 国产成人一区二区三区免费视频网站| 91成人精品电影| 桃花免费在线播放| 国产精品香港三级国产av潘金莲| av天堂久久9| 国产在线视频一区二区| 大片免费播放器 马上看| 国产伦人伦偷精品视频| 国产精品98久久久久久宅男小说| 高清毛片免费观看视频网站 | 极品教师在线免费播放| 免费日韩欧美在线观看| 一边摸一边抽搐一进一小说 | 午夜福利视频在线观看免费| 大码成人一级视频| 亚洲成人国产一区在线观看| 国产免费视频播放在线视频| 99在线人妻在线中文字幕 | 日日夜夜操网爽| 美女视频免费永久观看网站| 免费观看a级毛片全部| 日韩免费av在线播放| 极品人妻少妇av视频| 国产真人三级小视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 十八禁网站免费在线| 9191精品国产免费久久| cao死你这个sao货| 999久久久精品免费观看国产| 波多野结衣一区麻豆| 国产野战对白在线观看| 亚洲综合色网址| 国产精品免费一区二区三区在线 | 99香蕉大伊视频| 免费少妇av软件| 黑人欧美特级aaaaaa片| 日韩视频一区二区在线观看| 亚洲自偷自拍图片 自拍| 亚洲人成电影观看| 久久人人97超碰香蕉20202| 亚洲精品乱久久久久久| 1024香蕉在线观看| 亚洲自偷自拍图片 自拍| 国产熟女午夜一区二区三区| 久久精品国产亚洲av香蕉五月 | 成人av一区二区三区在线看| 少妇的丰满在线观看| 国产精品九九99| 日韩中文字幕欧美一区二区| 熟女少妇亚洲综合色aaa.| 久久中文字幕人妻熟女| 久久久久久人人人人人| 丰满迷人的少妇在线观看|