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

    基于大渦模擬與被動示蹤物模型的人工魚礁數(shù)值研究*

    2021-11-27 06:56:36王者也
    海洋與湖沼 2021年6期
    關(guān)鍵詞:上升流魚礁來流

    王者也 李 爽

    基于大渦模擬與被動示蹤物模型的人工魚礁數(shù)值研究*

    王者也 李 爽①

    (浙江大學(xué)海洋學(xué)院 舟山 316021)

    投放人工魚礁是解決海洋生態(tài)環(huán)境問題的重要措施之一。通過使用并行大渦模擬模式(the parallelizeda large-eddy simulation model, PALM)及被動示蹤物模型模塊, 研究了不同流速條件(0.1, 0.2, 0.4和0.6 m/s)下方型人工魚礁對流場形態(tài)、營養(yǎng)鹽的抬升作用、和湍流動能收支的影響。研究表明, 魚礁的存在使得其附近垂向速度增大, 產(chǎn)生上升流。受到上升流的抬升作用, 魚礁底部的示蹤物迅速進(jìn)入海洋上層, 之后遇到魚礁后方的背渦流, 示蹤物的抬升受阻, 高度逐漸降低。在上升流區(qū)域以及背渦流區(qū)域的共同影響下, 示蹤物抬升區(qū)域的最大高度與來流流速無關(guān)。不同上升流定義對應(yīng)的上升流區(qū)域的高度與來流流速均不相關(guān); 上升流區(qū)域的最大速度、平均速度與來流流速都成線性增加的關(guān)系; 然而隨著來流流速的增大, 不同上升流定義對應(yīng)的上升流區(qū)域面積的變化趨勢卻完全不同。投放魚礁后, 魚礁區(qū)域底部的混合增強, 區(qū)域底部的能量被輸運至上層。這說明, 魚礁的存在不僅對營養(yǎng)物質(zhì)具有抬升作用, 還能將能量向上輸運。

    大渦模擬; 人工魚礁; 被動示蹤物; 湍流動能收支

    隨著海洋生態(tài)環(huán)境逐步惡化, 許多養(yǎng)殖場不再適合養(yǎng)殖魚類、軟體動物和甲殼類動物。事實上, 因為大多數(shù)漁業(yè)資源已經(jīng)被過度開發(fā), 當(dāng)前漁業(yè)系統(tǒng)的可持續(xù)性正受到質(zhì)疑( Garcia, 2001)。聯(lián)合國糧食及農(nóng)業(yè)組織調(diào)查(FAO, 2020)顯示, 盡管水產(chǎn)養(yǎng)殖業(yè)依舊在發(fā)展, 海洋捕魚業(yè)似乎已經(jīng)到達(dá)上限。就資源開發(fā)的角度而言, 全球生態(tài)系統(tǒng)壓力過大, 生物多樣性減少, 并面臨崩潰的危險。因為資源的逐漸枯竭, 人造建筑物被用來增加魚類豐度以及多樣性, 改善棲息環(huán)境以及修復(fù)受損的魚礁(Lan, 2004)。目前, 存在著很多關(guān)于人工魚礁的定義。Seaman (2000)將人工魚礁定義為放置在海底的自然或人造物體, 其作用是調(diào)控海洋環(huán)境中的生物和物理參數(shù)。Grace (2001)研究了影響人工魚礁壽命的環(huán)境因素。大量具有不同材料成分和結(jié)構(gòu)的人工魚礁在世界各地被廣泛利用(Walker, 2002)。海洋是由眾多物種組成的生態(tài)系統(tǒng), 從簡單的浮游植物到更高營養(yǎng)層的魚類, 形成一個相互影響和相互依賴的海洋食物網(wǎng)。人工魚礁的投放可以改變流場、顆粒物分布, 有利于營養(yǎng)鹽的積累(Davis, 1982; Ambrose, 1990)。因此, 魚礁可以形成適宜魚類生存、繁殖, 躲避天敵的區(qū)域。魚群會被魚礁吸引, 然后逐漸形成人工生態(tài)系統(tǒng)。

    在開放的環(huán)境中投放魚礁將會改變原來穩(wěn)定的流場, 進(jìn)而促進(jìn)水體交換。魚礁的尺寸和結(jié)構(gòu)的多樣性導(dǎo)致不同地區(qū)的魚礁附近的流場不同。Vicente等(2008)對部署在不同地區(qū)的人工魚礁系統(tǒng)進(jìn)行營養(yǎng)循環(huán)和海底有機物富集方面的環(huán)境評估。事實上水體交換不僅可以為魚類提供不同的流體環(huán)境, 還可以促進(jìn)營養(yǎng)鹽的輸運。因此, 從修復(fù)海洋環(huán)境到豐富沿海生態(tài)系統(tǒng), 人工魚礁都是巖礁魚類的理想棲息地(Haro, 2004; Lan, 2004)。

    投放人工魚礁會影響魚礁區(qū)域的流速分布, 并產(chǎn)生不容忽視的湍流。研究人工魚礁對魚礁區(qū)域影響的方法主要分為兩種: 物理模型試驗和數(shù)值模擬。早期研究人員使用物理模型試驗對人工魚礁附近流場形態(tài)變化進(jìn)行初步分析。此后, 許多研究采用復(fù)雜的流體動力學(xué)模型, 通過數(shù)值模擬來探究流速的變化和渦流效應(yīng)。在投放人工魚礁前, 通過數(shù)值模擬對魚礁投放后產(chǎn)生的流場效應(yīng)進(jìn)行分析, 可以節(jié)省費用以及時間, 但模擬結(jié)果的可靠性需要驗證。Liu等(2006)利用有限體積法(finite volume method, FVM)對魚礁附近流場進(jìn)行了數(shù)值研究。Su等(2007)利用粒子成像測速儀(particle image velocimetry, PIV)對魚礁流場進(jìn)行測量與分析, 并對Liu等(2006)的數(shù)值結(jié)果進(jìn)行了進(jìn)一步的比較和驗證。Su等(2008)還研究了海底地形對人工魚礁內(nèi)外流場的影響。Jiang (2010)利用(renormalization group theory, RNG)-湍流模型和水槽實驗?zāi)P蛯误w魚礁產(chǎn)生的流場進(jìn)行分析, 得出模型結(jié)果可以代替試驗結(jié)果。

    經(jīng)魚礁提高的漁業(yè)資源的生態(tài)效率在一定程度上取決于人工魚礁的流場效應(yīng), 因此, 流場分析是研究人工魚礁生態(tài)效率的關(guān)鍵。上升流和背渦流的規(guī)模成為眾多研究中衡量魚礁流場效應(yīng)的主要標(biāo)準(zhǔn)。湍流的存在對于海洋上混合層的物理過程有著非常重要的影響。然而多數(shù)研究(Su, 2008; Liu, 2013; Li, 2017)并未分析湍流對于魚礁區(qū)域流場形態(tài)的影響。湍流在不同尺度上對浮游生物產(chǎn)生影響(Ki?rboe, 1993)。在近岸海域湍流能夠增加營養(yǎng)鹽通量(Dade, 1993)。一些研究數(shù)據(jù)表明湍流(即混合)與浮游植物群落組成之間存在聯(lián)系(Prairie, 2012; Acevedo-Trejos, 2015)。湍流可能彌補了它們?nèi)狈ψ晕彝七M(jìn)的器官的不足, 有利于它們遇到營養(yǎng)物質(zhì), 以及在真光層的持續(xù)存在(Margalef, 1978)?;旌夏軌蛟鰪姷撞繝I養(yǎng)鹽的擴散, 同時可能改變浮游植物的光照條件。本文將通過使用并行大渦模擬模式(the parallelized large-eddy simulation model, PALM)對魚礁區(qū)域的湍流動能收支、流場形態(tài)進(jìn)行分析。此外, 文章引入營養(yǎng)鹽模型, 以探究投放魚礁后的示蹤物分布變化。

    1 模型設(shè)計

    1.1 模式介紹

    本文使用的PALM基于Boussinesq近似, 經(jīng)過濾波的不可壓縮的非流體靜力學(xué)Navier-Stokes方程進(jìn)行求解。其基本方程如下:

    由于營養(yǎng)鹽與被動示蹤物在水體中的運動特性相似, 故將PALM中的被動示蹤物視為研究區(qū)域中的營養(yǎng)鹽, 其控制方程如下:

    其中, 上劃線表示濾波后去除次網(wǎng)格項的值; 雙撇號表示次網(wǎng)格尺度(subgrid-scale, SGS)變量; 尖括號表示SGS通量; 表1詳細(xì)列舉了方程的參數(shù)說明。上述方程通過離散笛卡爾網(wǎng)格上的控制方程對變量求平均值, 可以實現(xiàn)分辨尺度和SGS湍流模式的隱式分離, 其中采用Deardorff(1980)提出的方法建立SGS, 為方便起見, 本文使用了方程的連續(xù)形式。模式的詳細(xì)介紹見Huq等(2018)與Maronga等(2015)。

    本文通過計算魚礁上方各深度的上升流通量(Jiang, 2020)來分析不同來流流速下魚礁的抬升效果。

    其中,up表示上升流區(qū)域的面積;表示上升流區(qū)域的垂向速度。若流場中垂向速度大于或等于來流流速in的5%, 即≥0.05in, 則可稱該區(qū)域為上升流區(qū)域(Jiang, 2020)。

    1.2 礁體模型與計算區(qū)域

    人工魚礁結(jié)構(gòu)繁多, 如箱體型、三角型、組合型等(楊吝等, 2005)。在我國的東南部沿海海域, 魚礁多以米字型、立方體型為主(姜昭陽等, 2019)。在我國東部海域, 方型魚礁已被廣泛使用。李磊等(2018)通過研究箱體礁、三角形礁、框架礁對黑棘鯛的誘集作用, 發(fā)現(xiàn)箱體礁的誘集效果最佳。在近海區(qū)域, 付東偉等(2014)研究發(fā)現(xiàn), 開口小的單體魚礁產(chǎn)生的流場效應(yīng)更好。林軍等(2006)研究表明, 對于方型單礁來說, 其緩流區(qū)域體積與通透系數(shù)成反比。因此, 本次模擬選取邊長為3.0 m的實心方型人工魚礁作為研究目標(biāo)。趙海濤等(2006)認(rèn)為, 由于波浪對魚礁的沖擊作用, 魚礁應(yīng)投放在水深20—30 m的海域。結(jié)合人工魚礁實際案例(趙榮榮, 2019; 馮英明等, 2020), 本文的模擬區(qū)域水深為20 m。

    表1 參數(shù)列表

    Tab.1 List of parameters

    為避免人工魚礁流場遭受干擾(Ong, 2010), 本文將計算域的長度設(shè)置為100 m, 寬度設(shè)置為100 m。模型網(wǎng)格設(shè)置為200×200×80 (××), 網(wǎng)格大小設(shè)置為0.5×0.5×0.25 (m3), 即100 m×100 m×20 m (圖1)。參考Jiang等(2020)、Li等(2017)以及Liu等(2013)研究中, 模型的來流流速設(shè)置。本文設(shè)置四組理想實驗, 其來流流速為0.1、0.2、0.4和0.6 m/s; 并設(shè)置對照組, 即區(qū)域不設(shè)置魚礁, 其余的條件與實驗組一致。入流邊界條件為非循環(huán)邊界條件; 表層、底層邊界采用自由滑移邊界條件。

    曹欣中(1983)研究中的觀測數(shù)據(jù)詳細(xì)地展示了我國東部沿海海域表層營養(yǎng)鹽各深度的分布情況。為使此次研究的模擬結(jié)果更具代表性, 本文參考曹欣中(1983)的觀測結(jié)果, 設(shè)置了模擬區(qū)域各深度被動示蹤物的濃度分布。被動示蹤物濃度的分布見圖2。為與模型初始場使用的被動示蹤物濃度觀測數(shù)據(jù)的來源區(qū)域相對應(yīng), 此次研究設(shè)置的緯度為30°N。

    圖1 計算域示意圖

    圖2 初始被動示蹤物剖面

    注:表示被動示蹤物;表示深度;代表魚礁高度

    2 數(shù)值模擬結(jié)果與分析

    2.1 垂向流場形態(tài)

    2.1.1 上升流區(qū)域 魚礁對營養(yǎng)鹽的抬升作用依賴于魚礁區(qū)域的垂向流場。在模式穩(wěn)定時間內(nèi), 通過將魚礁區(qū)的垂向速度進(jìn)行水平()方向平均并與對照組相減得到圖3, 以研究魚礁區(qū)域垂向速度與上升流區(qū)域的垂向分布。從圖3可以看出, 魚礁的存在可以增大魚礁前方的垂向速度。在垂向上, 垂向速度增大區(qū)域的面積隨著來流流速的增大而增大。當(dāng)來流流速達(dá)到0.6 m/s時, 該區(qū)域的高度可達(dá)到水體表面, 實現(xiàn)了對全水深流場的影響。不同來流流速對應(yīng)的上升流部分在垂向上大體呈扇形且面積相近。當(dāng)來流流速為0.2 m/s時, 魚礁后方出現(xiàn)微弱的上升流, 其面積隨來流流速的增大而增大。在魚礁后方, 各來流流速對應(yīng)的垂向流速方向均被迫向下, 形成背渦流區(qū)域。

    圖3 來流流速為0.1 (a)、0.2 (b)、0.4 (c)、和0.6 m/s (d)時, 實驗組與對照組的垂向速度w在x-z方向上的差值

    注: 灰線包裹區(qū)域為上升流區(qū)域; 黑色方塊表示魚礁;代表模擬區(qū)域長度;表示深度;代表魚礁高度

    圖4展示了來流流速為0.2 m/s時, 不同深度下(-6.00、-5.67、-5.33、-5.00、-4.67、-4.33), 魚礁區(qū)的垂向速度與上升流區(qū)域在水平方向的分布。從圖4a可以看出, 上升流區(qū)域位于魚礁側(cè)后方且面積較小。對比圖4a—4f可知, 隨著深度的變淺, 上升流區(qū)域的水平分布面積在垂向上先增大后減小。魚礁后方垂向速度增大的部分與減小的部分相互脫離。在水平方向上, 魚礁上方的上升流部分形狀為圓形, 中心處的垂向速度最大, 越靠近外圍, 垂向速度越小。

    圖4 來流流速為0.2 m/s時, 實驗組與對照組在深度為-6.00L (a)、-5.67L (b)、-5.33L (c)、-5.00L (d)、-4.67L (e)、-4.33L (f)處, 垂向速度w在x-y方向上的差值

    注: 灰線包裹區(qū)域為上升流區(qū)域; 黑色方塊表示魚礁;代表模擬區(qū)域?qū)挾?代表模擬區(qū)域長度;代表魚礁高度

    為進(jìn)一步探究上升流區(qū)域的規(guī)模與強度, 文章通過使用參數(shù)a、max、up、up對上升流區(qū)域的規(guī)模和強度進(jìn)行描述。上升流最大速度(max)、上升流平均速度(a)被用來形容上升流的強度; 上升流高度(up)、上升流面積(up)被用來形容上升流的規(guī)模。

    為了使結(jié)論具有代表性, 增加兩組來流流速(0.3、0.5 m/s)。從圖5a中可以看出, 上升流的最大速度、平均速度均與來流流速(in)成線性關(guān)系, 其線性關(guān)系式分別為max=0.71in、a=0.08in, 相關(guān)系數(shù)均為1.00。從圖5b中可以看出, 隨著來流流速的增大, 上升流的高度基本不變, 均在2.42(7.25 m)附近。因此,up與in相關(guān)性不大。上升流面積與來流流速的關(guān)系可以用指數(shù)回歸多項式進(jìn)行表示。表達(dá)式為up=0.0013e14.1Um+32.92, 相關(guān)系數(shù)為0.97。

    2.1.2 上升流通量 除了關(guān)注上升流區(qū)域與垂向速度在垂向以及水平方向上的分布之外, 文章還關(guān)注了魚礁產(chǎn)生的上升流通量。本文選取魚礁上方不同高度來分析各來流流速對應(yīng)的上升流通量變化。

    圖5 上升流平均速度Wa及最大速度Wmax, 上升流高度Hup及面積Sup與來流流速Uin之間的關(guān)系

    從圖6中可以看出, 魚礁上方, 各來流流速對應(yīng)上升流通量在垂向上的變化趨勢相近, 即隨著深度的增大,先增大后減少, 與(Jiang, 2020)研究結(jié)果類似。隨著來流流速的增大, 魚礁附近的上升流通量逐漸增大。圖4中的上升流區(qū)域垂向面積的變化趨勢(先增大后減小)解釋了上升流通量的在垂向的變化趨勢。

    2.2 湍流動能收支

    根據(jù)公式(1)可得水平平均的SGS-TKE(turbulent kinetic energy)方程(Skyllingstad, 2000)為:

    其中,

    表示壓強,m表示次網(wǎng)格的渦黏系數(shù),表示重力加速度,表示湍動能, Δ表示耗散項。公式(7)的等號右邊, 第一項表示速度剪切項、第二項表示浮力項、第三項表示次網(wǎng)格耗散項、第四項和第五項分別表示湍流動能傳輸項和壓力項、第六項表示耗散項。下文中湍流動能的耗散項包含次網(wǎng)格耗散, 即兩項之和。

    人工魚礁的投放不僅會影響魚礁區(qū)域的流場分布, 還會產(chǎn)生不容忽視的湍流。圖7展示了模型的湍流動能收支情況。根據(jù)公式(8), 本文用紅線表示速度剪切項, 綠線表示湍流動能傳輸項, 藍(lán)線表示壓力項, 黑線表示耗散項。左圖為四種來流流速下的海洋湍流動能收支情況, 右圖均為左圖底部區(qū)域的局部放大圖。從圖7a、7c、7e、7g中可以發(fā)現(xiàn), 海洋湍流動能收支各項的數(shù)值在深度為-5.5處迅速增大, 然后迅速減小, 尤其是剪切項與耗散項的變化最為劇烈。此現(xiàn)象說明, 當(dāng)來流遇到魚礁后會產(chǎn)生強烈的湍流, 進(jìn)而增強混合。在圖7b中, 壓力項與湍流動能傳輸項的值在近海底處差別不大, 耗散項的最小值可以達(dá)到-1.14×10–6m2/s3; 在圖7d中, 剪切項與耗散項絕對值的最大值位置接近海底且耗散項的數(shù)值較大, 約為4.52×10–6m2/s3。在圖7f中, 耗散項在-6.33處達(dá)到最小值, 最小值為-1.89×10–5m2/s3, 在圖7h中, 對比剪切項和耗散項絕對值的最大值可知, 耗散項約為剪切項的1.5倍, 即4.46×10–5m2/s3。

    圖6 來流流速為0.1、0.2、0.4和0.6 m/s時, 上升流通量Q在不同深度處的分布

    注:表示深度;代表魚礁高度

    隨著來流流速的增加, 在近海底處的剪切、傳輸與壓力項必然增大, 導(dǎo)致向上傳遞的距離越遠(yuǎn), 所產(chǎn)生的耗散也相應(yīng)增大, 且隨著深度減小逐漸減弱。

    圖7 來流流速為0.1 (a, b)、0.2 (c, d)、0.4 (e, f)、和0.6 m/s (g, h)時水體底部的湍流動能收支

    注: 右圖均為左圖底部區(qū)域的局部放大圖;表示深度;代表魚礁高度

    2.3 被動示蹤物

    2.3.1 被動示蹤物分布 為驗證人工魚礁產(chǎn)生的上升流對被動示蹤物是否有抬升作用, 將實驗組減去對照組的被動示蹤物濃度得到圖8, 以探究魚礁被動示蹤物在垂向上的抬升效果。圖8a—8d中人工魚礁后方, 接近底部的被動示蹤物濃度為負(fù)值, 高于人工魚礁的部分被動示蹤物濃度為正值, 說明相較于未放置人工魚礁的對照組, 在人工魚礁的抬升作用下, 底部更多的被動示蹤物進(jìn)入上層水體。在人工魚礁上部, 被動示蹤物的抬升區(qū)域與上升流區(qū)域在垂向上均成扇形。

    圖8 來流流速為0.1 (a)、0.2 (b)、0.4 (c)、和0.6 m/s (d)時, 實驗組與對照組在x-z方向的被動示蹤物濃度S的差值

    注: 灰線包裹區(qū)域為上升流區(qū)域; 黑色方塊表示魚礁;代表模擬區(qū)域長度;表示深度;代表魚礁高度

    各來流流速對應(yīng)的被動示蹤物抬升區(qū)域, 在人工魚礁后方的不同位置處連續(xù)性明顯減弱且抬升高度降低。造成此現(xiàn)象的原因為, 在距人工魚礁一段距離后, 流場更加紊亂且垂向速度為負(fù)。

    圖9展示了來流流速為0.2 m/s時, 不同深度下(-6.00、-5.67、-5.33、-5.00、-4.67、-4.33), 被動示蹤物的水平分布。從圖9a可以看出, 受到在魚礁的抬升作用, 接近魚礁底部的營養(yǎng)鹽被水體裹挾進(jìn)入上層水體, 使得魚礁后方的營養(yǎng)鹽濃度明顯降低。從圖9b—9f中可以看出距離中心越遠(yuǎn)的位置, 營養(yǎng)鹽的抬升效果越弱。在水平方向上, 相較于圖4中上升流面積的變化趨勢, 魚礁上方營養(yǎng)鹽的抬升面積不隨深度的變化而變化。

    為進(jìn)一步探究魚礁對被動示蹤物的抬升效果, 在人工魚礁上方0以及后方5、10、15處, 將實驗組的被動示蹤物濃度與對照組相減得到圖10。隨著深度的加深, 被動示蹤物的抬升效果基本呈現(xiàn)先緩慢增大后迅速減小的趨勢。

    圖10a—10d中, 各來流流速對應(yīng)的被動示蹤物濃度分布在魚礁上方十分接近。在距離魚礁5處, 隨著來流流速的增大, 被動示蹤物濃度有所降低。在距離魚礁10處, 來流流速為0.1、0.2、0.4 m/s時對應(yīng)被動示蹤物濃度不斷降低; 當(dāng)來流流速達(dá)到0.6 m/s時, 被動示蹤物的抬升作用明顯增強, 相較于來流流速為0.4 m/s, 抬升的被動示蹤物濃度增大約2倍左右。在距離魚礁15處, 隨著來流流速的增大, 被動示蹤物的抬升作用先減弱后增強。對比圖10a—10d中的示蹤物分布, 可以發(fā)現(xiàn)當(dāng)來流流速大于0.2 m/s時, 魚礁的抬升效果有明顯的提升, 與圖8中被動示蹤物抬升區(qū)域的連續(xù)性相符合。

    2.3.2 被動示蹤物通量 文章還通過計算魚礁上方的上升流區(qū)域?qū)?yīng)的被動示蹤物通量up來衡量人工魚礁的抬升效果。

    其中,表示上升流通量;表示上升流區(qū)域up所含被動示蹤物濃度。研究選取魚礁上方不同高度來分析各來流流速對應(yīng)上升流區(qū)域的被動示蹤物通量的up變化。

    從圖11中可以看出, 魚礁上方, 各來流流速對應(yīng)被動示蹤物通量與上升流通量變化趨勢相近。隨著來流流速的增大, 魚礁附近的被動示蹤物通量up逐漸增大。

    3 討論

    3.1 上升流定義的對比

    黃遠(yuǎn)東等(2012)利用計算流體力學(xué)(computational fluid dynamics, CFD)技術(shù)研究不同流速下, 方型魚礁對應(yīng)的上升流和背渦流的分布情況, 發(fā)現(xiàn)方型魚礁的上升流高度和背渦流區(qū)寬度與來流流速無關(guān); 唐衍力等(2017)研究進(jìn)一步發(fā)現(xiàn), 不同人工魚礁結(jié)構(gòu)的上升流及背渦流面積的大小與來流流速不相關(guān), 最大上升流流速與來流流速呈線性關(guān)系。于定勇等(2019)發(fā)現(xiàn)隨著魚礁開口比的增大, 礁體產(chǎn)生上升流范圍以及垂向速度逐漸減少。

    圖9 來流流速為0.2 m/s時, 實驗組與對照組在深度為-6.00L (a)、-5.67L (b)、-5.33L (c)、-5.00L (d)、-4.67L (e)、-4.33L (f)處, 被動示蹤物S濃度在x-y方向上的差值

    注: 灰線包裹區(qū)域為上升流區(qū)域; 黑色方塊表示魚礁;代表模擬區(qū)域?qū)挾?代表模擬區(qū)域長度;代表魚礁高度

    圖10 來流流速為0.1 (a)、0.2 (b)、0.4 (c)、和0.6 m/s (d)時, 在距人工魚礁0L (黑線)、5L (黃線)、15L (藍(lán)線)、20L (紅線)處, 實驗組與對照組的被動示蹤物S濃度的差值

    圖11 來流流速為0.1、0.2、0.4和0.6 m/s時, 上升流區(qū)域?qū)?yīng)的營養(yǎng)鹽通量Fup在不同深度處的分布

    注:表示深度;代表魚礁高度

    與唐衍力等(2017)研究結(jié)果不同, 圖5中上升流面積與來流流速的關(guān)系為指數(shù)遞增。由于本文將垂向速度大于或等于來流流速in的5%的區(qū)域定義為上升流區(qū)域, 而唐衍力等(2017)則將垂向流速分量與來流流速之比等于或大于10%的區(qū)域定義為上升流區(qū)域。因此, 作者按照唐衍力等(2017)中的上升流的定義分別對上升流最大速度(max)、上升流平均速度(a)及上升流高度(up)、上升流面積(up)進(jìn)行統(tǒng)計得到圖12。圖12a中可以看出上升流的最大速度、平均速度均與來流流速成線性關(guān)系, 其線性關(guān)系式分別為max=0.71in、a=0.15in、相關(guān)系數(shù)均為1.00。從圖12b中可以看出隨著來流流速的增大, 上升流區(qū)域高度以及上升流區(qū)域面積基本不變, 均在2.42L(7.25m)以及1.44L2(13m2)附近, 與唐衍力等(2017)的研究結(jié)果一致。因此, 兩次研究結(jié)果出現(xiàn)差異的原因為對于上升流的定義不同。

    圖3與圖4表明, 在投放魚礁后, 魚礁附近的垂向流速明顯增大, 從而產(chǎn)生上升流。文章通過引用兩種定義方法(唐衍力等, 2017; Jiang, 2020)對上升流區(qū)域進(jìn)行描述, 分別得到圖5以及圖12。由于上升流區(qū)域的最大速度與上升流區(qū)域的定義無關(guān), 因此, 對比圖5與圖12可知, 在兩種不同的上升流定義下, 除了上升流區(qū)域的最大速度(max)相同之外, 上升流區(qū)域的平均速度(a)、高度(up)以及面積(up)的數(shù)值均不相同。值得注意的是, 兩種上升流定義對應(yīng)的上升流的最大流速(max)及平均流速(a)與來流流速均成線性增加的關(guān)系, 上升流的高度(up)均不隨來流流速的改變而改變, 隨著來流流速的增大, 上升流區(qū)域面積(up)的變化趨勢卻完全不同。

    圖12 改變上升流的定義后, 上升流平均速度Wa及最大速度Wmax, 上升流高度Hup及面積Sup與來流流速Uin之間的關(guān)系

    3.2 被動示蹤物的抬升效果

    在圖3中, 魚礁后方存在明顯的背渦流區(qū)域, 且隨著來流流速的增大, 魚礁產(chǎn)生背渦流隨之增強。這就使得魚礁對示蹤物的抬升效果受到限制。當(dāng)示蹤物遇到魚礁產(chǎn)生的上升流后, 會被迅速抬升至上層水體; 然后遇到魚礁后方的背渦流, 示蹤物的抬升受阻, 高度降低。因此, 從圖8中可以看出, 在垂向上各來流流速對示蹤物的最大抬升高度基本相同, 魚礁后方示蹤物的抬升高度有所降低, 示蹤物的抬升區(qū)域的面積基本相同。

    在圖10中, 魚礁底部對應(yīng)各流速抬升的示蹤物濃度均為負(fù)值, 即底部的示蹤物被抬升到魚礁上方。同時, 圖7中綠線代表的湍流傳輸項顯示, 區(qū)域底部的傳輸項在垂向上的數(shù)值變化均為先負(fù)后正, 即底部的能量向上運輸。此現(xiàn)象表明, 投置魚礁后, 區(qū)域底部的動能被帶往上層水體耗散。因此, 魚礁的存在使得底部能量向上輸運, 底層的示蹤物也隨之被帶往上層, 體現(xiàn)了魚礁對于投放海域底層的能量與營養(yǎng)物質(zhì)均有抬升作用。

    寧修仁等(2007)經(jīng)過現(xiàn)場觀測和取樣, 得到珠江口冬季真光層深度大約為2.9±1.5 m。陳雨(2014)發(fā)現(xiàn), 冬季長江口附近海域真光層深度在0.31—12.42 m之間變動, 且離岸越近, 真光層深度越小。從圖7中可以看出, 在近岸真光層范圍內(nèi)人工魚礁對于營養(yǎng)鹽均有抬升作用。

    4 結(jié)論

    為研究不同來流流速下, 方型人工魚礁對流場形態(tài)的影響以及對營養(yǎng)鹽的抬升作用, 本文設(shè)置四種來流流速(0.1、0.2、0.4和0.6 m/s), 分析不同流速下的流場形態(tài), 湍流動能收支各項的分布結(jié)果與被動示蹤物的分布狀態(tài), 并對比有無人工魚礁對被動示蹤物抬升效果的差別。主要結(jié)論如下:

    (1) 對于不同流速, 魚礁上方上升流通量隨著深度的變淺, 其變化趨勢均為先增大后減小; 隨著來流流速的增大, 魚礁上方的上升流通量隨之增大; 魚礁產(chǎn)生的上升流區(qū)域均呈扇形。魚礁的存在, 不僅使得其附近垂向速度增大, 還使其后方垂向速度減小。文章引用兩種定義方法(唐衍力等, 2017; Jiang, 2020) 分別對上升流區(qū)域進(jìn)行描述。結(jié)果發(fā)現(xiàn), 兩種上升流定義對應(yīng)的上升流區(qū)域的高度(up)均不隨來流流速的改變而改變, 上升流的最大流速(max)及平均流速(a)與來流流速均成線性增加的關(guān)系, 然而隨著來流流速的增大, 上升流區(qū)域面積(up)的變化趨勢卻完全不同。

    (2) 受到上升流的抬升作用, 魚礁底部的示蹤物迅速進(jìn)入海洋上層, 而后遇到魚礁后方的背渦流, 示蹤物的抬升受阻, 高度逐漸降低。上升流區(qū)域以及背渦流區(qū)域的共同影響下, 示蹤物抬升區(qū)域的最大高度與來流流速無關(guān); 且在魚礁后方, 示蹤物抬升區(qū)域的高度有所降低。各流速下, 魚礁上升流區(qū)域的營養(yǎng)鹽通量的變化趨勢與上升流通量相似。隨著與魚礁的距離增大, 魚礁的抬升作用先增強后減弱; 魚礁上方的抬升效果最佳。在近海區(qū)域的真光層, 魚礁可以實現(xiàn)對營養(yǎng)鹽的抬升作用。

    (3) 文章簡略探討了魚礁產(chǎn)生湍流對海洋上層的影響。剪切項、壓力項以及傳輸項為區(qū)域底部湍流動能的生成源, 平衡增大的耗散項主要靠剪切項。投放魚礁后, 魚礁區(qū)域底部的混合增強, 同時湍流剪切項與耗散項顯著增大。湍流傳輸項的分布顯示區(qū)域底部的能量被輸運至上層, 這說明, 魚礁的存在不僅對營養(yǎng)物質(zhì)具有抬升作用, 還能將能量向上層運輸; 此外, 它引起的上層混合也能使浮游生物的分布等進(jìn)行重新分配, 進(jìn)而形成特殊的生態(tài)區(qū)域。

    本文通過使用并行大渦模擬模式PALM及其被動示蹤物模塊, 初步探索了魚礁區(qū)域的流場形態(tài)與湍流過程對示蹤物空間分布的影響。潮流在我國近海十分顯著, 對近岸海洋環(huán)境的水動力過程影響較大, 未來工作重點便是在模式中加入周期性的潮流動力環(huán)境, 考慮潮流作用對魚礁區(qū)域的影響。

    于定勇, 楊遠(yuǎn)航, 李宇佳, 2019. 不同開口比人工魚礁體水動力特性及礁體穩(wěn)定性研究. 中國海洋大學(xué)學(xué)報, 49(4): 128—136

    付東偉, 陳 勇, 陳衍順等, 2014. 方形人工魚礁單體流場效應(yīng)的PIV試驗研究. 大連海洋大學(xué)學(xué)報, 29(1): 82—85

    馮英明, 許丙彩, 郝 義等, 2020. 日照市海洋牧場示范區(qū)人工魚礁選址適宜性分析. 山東國土資源, 36(1): 44—50

    李 磊, 陳 棟, 彭建新等, 2018. 3種人工魚礁模型對黑棘鯛的誘集效果研究. 海洋漁業(yè), 40(5): 625—631

    楊 吝, 劉同渝, 2005. 我國人工魚礁種類的劃分方法. 漁業(yè)現(xiàn)代化, (6): 22—23, 25

    陳 雨, 2014. 長江口鄰近海域冬季水體漫衰減特性及遙感反演初探. 上海: 華東師范大學(xué)碩士學(xué)位論文, 41—42

    林 軍, 章守宇, 2006. 人工魚礁物理穩(wěn)定性及其生態(tài)效應(yīng)的研究進(jìn)展. 海洋漁業(yè), 28(3): 257—262

    姜昭陽, 郭戰(zhàn)勝, 朱立新等, 2019. 人工魚礁結(jié)構(gòu)設(shè)計原理與研究進(jìn)展. 水產(chǎn)學(xué)報, 43(9): 1881—1889

    劉誠剛, 寧修仁, 蔡昱明等, 2007. 南海北部及珠江口細(xì)菌生產(chǎn)力研究. 海洋學(xué)報(中文版), 2007(2): 112—122

    趙榮榮, 2019. 長島擋浪島人工魚礁生態(tài)修復(fù)效果初步評價. 舟山: 浙江海洋大學(xué)碩士學(xué)位論文, 7—8

    趙海濤, 張亦飛, 郝春玲等, 2006. 人工魚礁的投放區(qū)選址和礁體設(shè)計. 海洋學(xué)研究, 24(4): 69—76

    唐衍力, 龍翔宇, 王欣欣等, 2017. 中國常用人工魚礁流場效應(yīng)的比較分析. 農(nóng)業(yè)工程學(xué)報, 33(8): 97—103

    黃遠(yuǎn)東, 姜劍偉, 趙樹夫, 2012. 方型人工魚礁周圍水流運動的數(shù)值模擬研究. 水資源與水工程學(xué)報, 23(3): 1—3

    曹欣中, 1983. 影響浙江近海營養(yǎng)鹽分布主要水動力因子的探討. 海洋科學(xué), 1983, 7(2): 36—38

    Acevedo-Trejos E, Brandt G, Bruggeman J, 2015. Mechanisms shaping size structure and functional diversity of phytoplankton communities in the ocean. Scientific Reports, 5(1): 8918

    Ambrose R F, Anderson T W, 1990. Influence of an artificial reef on the surrounding infaunal community. Marine Biology, 107(1): 41—52

    Dade W B, 1993. Near-bed turbulence and hydrodynamic control of diffusional mass transfer at the sea floor. Limnology and Oceanography, 38(1): 52—69

    Davis N, VanBlaricom G R, Dayton P K, 1982. Man-made structures on marine sediments: Effects on adjacent benthic communities. Marine Biology, 70(3): 295—303

    Deardorff J W, 1980. Stratocumulus-capped mixed layers derived from a three-dimensional model. Boundary-Layer Meteorology, 18(4): 495—527

    FAO, 2020. The state of world Fisheries and Aquaculture 2020. Rome: Food and Agriculture Organization of the United Nations.

    García S, Leiva Moreno I, 2003.Global overview of marine fisheries. Conference on Responsible Fisheries in the Marine Ecosystem, 103—123

    Grace R A, 2001. The factors and processes that influence artificial reef longevity. Marine Technology Society Journal, 35(2): 3—13

    Haro A, Castro-Santos T, Noreika J, 2004. Swimming performance of upstream migrant fishes in open-channel flow: a new approach to predicting passage through velocity barriers. Canadian Journal of Fisheries and Aquatic Sciences, 61(9): 1590—1601

    Huq S, De Roo F, Raasch S, 2019. Vertically nested les for high-resolution simulation of the surface layer in PALM (version 5.0). Geoscientific Model Development. 12, 2523—2538

    Jiang Z Y, Liang Z L, Tang Y L, 2010. Numerical simulation and experimental study of the hydrodynamics of a modeled reef located within a current. Chinese Journal of Oceanology and Limnology, 28(2): 267—273

    Jiang Z Y, Liang Z L, Zhu L X, 2020. Effect of hole diameter of rotary-shaped artificial reef on flow field. Ocean Engineering, 197: 106917

    Ki?rboe T, 1993. Turbulence, phytoplankton cell size, and the structure of pelagic food webs. Advances in Marine Biology, 29: 1—72

    Lan C H, Chen C C, Hsui C Y, 2004. An approach to design spatial configuration of artificial reef ecosystem. Ecological Engineering, 22(4—5): 217—226

    Li J, Zheng Y X, Gong P H, 2017. Numerical simulation and PIV experimental study of the effect of flow fields around tube artificial reefs. Ocean Engineering, 134, 96—104

    Liu T L, Liou S, Chen D W, 2006. Numerical flow simulation for the deployment analysis of artificial reefs. in Proceedings of the International Offshore and Polar Engineering Conference 533—537

    Liu Y, Zhao Y P, Dong G H, 2013. A study of the flow field characteristics around star-shaped artificial reefs. Journal of Fluids and Structures, 39: 27—40

    Margalef R, 1978. Life-forms of phytoplankton as survival alternatives in an unstable environment. Oceanologica Acta, 1(4): 493—509

    Maronga B, Gryschka M, Heinze R, 2015. The parallelized large-eddy simulation model (PALM) version 4.0 for atmospheric and oceanic flows: Model formulation, recent developments, and future perspectives. Geoscientific Model Development, 8(8): 2515—2551

    Ong M C, Utnes T, Holmedal L E, 2010. Numerical simulation of flow around a circular cylinder close to a flat seabed at high Reynolds numbers using a-model. Coastal Engineering, 57(10): 931—947

    Prairie J C, Sutherland K R, Nickols K J, 2012. Biophysical interactions in the plankton: a cross-scale review. Limnology and Oceanography: Fluids and Environments, 2(1): 121—145

    Seaman W Jr, 2000. Artificial reef evaluation: with application to natural marine habitats. London: CRC Press, 1—249

    Skyllingstad E D, Smyth W D, Crawford G B, 2000. Resonant wind-driven mixing in the ocean boundary layer. Journal of Physical Oceanography, 30(8): 1866—1890

    Su D T, Liu T L, Ou C H, 2007. A comparison of the PIV measurements and numerical predictions of the flow field patterns within an artificial reef. In: Proceedings of the 17th International Offshore and Polar Engineering Conference. Lisbon, Portugal: ISOPE, 2239—2245

    Su D T, Liu T L, Ou C H, 2008. Numerical investigation into effects of seabed topography on flows in and around artificial reefs. Fisheries Science, 74(2): 236—254

    Vicente M, Falc?o M, Santos M N, 2008. Environmental assessment of two artificial reef systems off southern Portugal (Faro and Olh?o): a question of location. Continental Shelf Research, 28(6): 839—847

    Walker B K, Henderson B, Spieler R E, 2002. Fish assemblages associated with artificial reefs of concrete aggregates or quarry stone offshore Miami Beach, Florida, USA. Aquatic Living Resources, 15(2): 95—105

    NUMERICAL SIMULATION OF CUBIC ARTIFICIAL REEF UNDER LARGE EDDY IN PASSIVE SCALAR MODEL

    WANG Zhe-Ye, LI Shuang

    ( Ocean College, Zhejiang University, Zhoushan 316021, China)

    Artificial reef is a widely used tool for marine ecosystem improvement. A numerical simulation was conducted in the parallel large eddy simulation model (PALM) in passive scalar mode to explore the effects of solid cubic artificial reef on flow field, turbulent kinetic energy (TKE) budget, and nutrient uplift under four inflow velocities (0.1, 0.2, 0.4, and 0.6 m/s). Results show that artificial reef significantly increases the vertical velocities near reef, and underpin inflow into upwelling. Scalars at the reef bottom are pushed up to the ocean surface quickly and then sunk by vortex at the reef rear. Under the joint effects of upwelling and vortex, the maximum height of uplift of passive scalar show no relation with inflow velocities. In addition, different definitions of upwelling present no influences on the heights and velocities of upwelling. Both average and maximum velocities of upwelling are increased linearly with inflow velocities. However, the area of upwelling change demonstrates a totally different trend. After reef is placed, mixing at the bottom is strengthened, dissipation and shear are remarkably enlarged, and energy at the bottom is transported to the upper layers. Therefore, artificial reef can not only uplift nutrients, but also transport energy upward.

    large eddy simulation; artificial reef; passive scalar; TKE (turbulent kinetic energy)

    * 國家自然科學(xué)基金項目, 41876003號; 國家重點研發(fā)計劃項目, 2017YFA0604101號。王者也, 碩士研究生, E-mail: 21934012@zju.edu.cn

    李 爽, 博士生導(dǎo)師, 副教授, E-mail: lshuang@zju.edu.cn

    2021-03-30,

    2021-07-28

    P731.26

    10.11693/hyhz20210300075

    猜你喜歡
    上升流魚礁來流
    夏季舟山上升流長期遙感觀測與分析
    人工魚礁區(qū)三種習(xí)見肉食性軟體動物營養(yǎng)生態(tài)位研究
    海洋通報(2022年5期)2022-11-30 12:04:06
    兩種典型來流條件下風(fēng)力機尾跡特性的數(shù)值研究
    能源工程(2022年2期)2022-05-23 13:51:48
    人工魚礁建設(shè)對海洋浮游植物群落變化的影響
    海洋通報(2020年5期)2021-01-14 09:27:04
    不同來流條件對溢洪道過流能力的影響
    遼東半島頂端海域上升流長期變化特征及影響因素*
    海洋與湖沼(2020年1期)2020-02-06 06:01:38
    人工魚礁在海洋牧場中的應(yīng)用
    珠江水運(2018年21期)2018-12-20 23:17:38
    瓊東上升流的年際變化及長期變化趨勢*
    海洋與湖沼(2016年1期)2016-01-15 03:50:49
    日照前三島人工魚礁對漁業(yè)資源影響的初步評價?
    彈發(fā)匹配驗證試驗系統(tǒng)來流快速啟動技術(shù)研究
    亚洲精品久久久久久婷婷小说| 高清欧美精品videossex| 日韩一区二区三区影片| 国产91精品成人一区二区三区 | 国产精品国产av在线观看| 天堂中文最新版在线下载| 国产国语露脸激情在线看| 91老司机精品| 婷婷色av中文字幕| 国产欧美亚洲国产| 99精品久久久久人妻精品| 久久影院123| 777久久人妻少妇嫩草av网站| 少妇精品久久久久久久| 狠狠婷婷综合久久久久久88av| 99香蕉大伊视频| 黑人欧美特级aaaaaa片| 欧美成人午夜精品| 91精品三级在线观看| 18在线观看网站| 亚洲欧洲日产国产| 亚洲精品成人av观看孕妇| 国产日韩欧美在线精品| 亚洲av欧美aⅴ国产| 亚洲精品乱久久久久久| 免费日韩欧美在线观看| 在线观看免费视频网站a站| 免费观看av网站的网址| av又黄又爽大尺度在线免费看| av电影中文网址| 亚洲专区中文字幕在线| 久久天堂一区二区三区四区| h视频一区二区三区| 日本猛色少妇xxxxx猛交久久| 欧美久久黑人一区二区| 啦啦啦中文免费视频观看日本| 99九九在线精品视频| 精品少妇一区二区三区视频日本电影| 精品亚洲乱码少妇综合久久| 欧美激情 高清一区二区三区| 永久免费av网站大全| 日本wwww免费看| 日本av手机在线免费观看| 国产精品久久久久久精品电影小说| 欧美+亚洲+日韩+国产| 亚洲国产精品999| 天堂中文最新版在线下载| 免费观看a级毛片全部| 黑人猛操日本美女一级片| 国产成人免费无遮挡视频| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 亚洲,欧美,日韩| 纵有疾风起免费观看全集完整版| 两性夫妻黄色片| xxx大片免费视频| 99久久综合免费| 在线亚洲精品国产二区图片欧美| 青春草视频在线免费观看| a级片在线免费高清观看视频| 国产av一区二区精品久久| 欧美精品人与动牲交sv欧美| 两个人看的免费小视频| 亚洲精品中文字幕在线视频| a级片在线免费高清观看视频| 国产精品成人在线| 在线观看免费高清a一片| 久久久久久久国产电影| 啦啦啦 在线观看视频| 欧美日韩一级在线毛片| 日日摸夜夜添夜夜爱| 午夜视频精品福利| 国产熟女欧美一区二区| 国产男女内射视频| 婷婷色av中文字幕| 亚洲成色77777| 国产一区二区在线观看av| 精品久久蜜臀av无| 99九九在线精品视频| 国产精品久久久久久精品古装| 午夜影院在线不卡| 麻豆乱淫一区二区| svipshipincom国产片| 精品久久蜜臀av无| 国产高清videossex| 久久久久久久精品精品| 国产在视频线精品| 亚洲精品在线美女| 美女午夜性视频免费| 蜜桃在线观看..| 女人久久www免费人成看片| 欧美在线一区亚洲| 成人手机av| 十分钟在线观看高清视频www| av天堂久久9| 国产精品熟女久久久久浪| 欧美日韩亚洲高清精品| 欧美激情高清一区二区三区| 欧美少妇被猛烈插入视频| 久久国产精品影院| 人妻一区二区av| 亚洲色图 男人天堂 中文字幕| 日日爽夜夜爽网站| 免费黄频网站在线观看国产| 亚洲激情五月婷婷啪啪| 免费一级毛片在线播放高清视频 | 国语对白做爰xxxⅹ性视频网站| 国产精品麻豆人妻色哟哟久久| 亚洲一码二码三码区别大吗| 精品一区二区三区av网在线观看 | 免费在线观看日本一区| 国产野战对白在线观看| 男女边摸边吃奶| 久久久精品94久久精品| 在线亚洲精品国产二区图片欧美| 男女免费视频国产| 青青草视频在线视频观看| 久久久久久人人人人人| 免费观看av网站的网址| 激情五月婷婷亚洲| 国产爽快片一区二区三区| 久久人妻福利社区极品人妻图片 | 亚洲成av片中文字幕在线观看| 69精品国产乱码久久久| 国产老妇伦熟女老妇高清| 考比视频在线观看| 日本一区二区免费在线视频| 亚洲国产精品成人久久小说| 一本一本久久a久久精品综合妖精| 别揉我奶头~嗯~啊~动态视频 | videos熟女内射| 亚洲欧美精品自产自拍| 国产精品一区二区免费欧美 | videosex国产| 欧美中文综合在线视频| 在线亚洲精品国产二区图片欧美| 交换朋友夫妻互换小说| 国精品久久久久久国模美| 新久久久久国产一级毛片| 大型av网站在线播放| 欧美大码av| 男人舔女人的私密视频| 色精品久久人妻99蜜桃| 男男h啪啪无遮挡| 国产成人系列免费观看| 男人操女人黄网站| 一本久久精品| kizo精华| 色婷婷av一区二区三区视频| 国产精品一二三区在线看| 看十八女毛片水多多多| 国产成人欧美在线观看 | 18禁黄网站禁片午夜丰满| 美国免费a级毛片| 80岁老熟妇乱子伦牲交| 国产成人一区二区三区免费视频网站 | 韩国高清视频一区二区三区| 美女脱内裤让男人舔精品视频| 欧美日韩亚洲国产一区二区在线观看 | a级片在线免费高清观看视频| 亚洲一区二区三区欧美精品| 18禁黄网站禁片午夜丰满| 亚洲av日韩在线播放| 亚洲国产欧美一区二区综合| 男的添女的下面高潮视频| 最新的欧美精品一区二区| 亚洲欧美一区二区三区黑人| 久久精品久久久久久噜噜老黄| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品国产av成人精品| 丝袜在线中文字幕| 午夜福利一区二区在线看| 亚洲人成77777在线视频| 亚洲五月婷婷丁香| 国产精品秋霞免费鲁丝片| 黑人欧美特级aaaaaa片| 亚洲精品第二区| 黄色一级大片看看| 成人三级做爰电影| 蜜桃国产av成人99| 人成视频在线观看免费观看| 精品国产乱码久久久久久小说| 国产av国产精品国产| tube8黄色片| 久久性视频一级片| 欧美精品啪啪一区二区三区 | av视频免费观看在线观看| 日本午夜av视频| 精品国产国语对白av| 国产主播在线观看一区二区 | 欧美av亚洲av综合av国产av| av福利片在线| a级毛片黄视频| 在线精品无人区一区二区三| 国产成人啪精品午夜网站| 老汉色∧v一级毛片| 你懂的网址亚洲精品在线观看| 久久午夜综合久久蜜桃| 搡老岳熟女国产| 如日韩欧美国产精品一区二区三区| 一本综合久久免费| 美女高潮到喷水免费观看| 一边亲一边摸免费视频| 十八禁高潮呻吟视频| 1024视频免费在线观看| 五月天丁香电影| 免费在线观看完整版高清| 色婷婷av一区二区三区视频| 女人被躁到高潮嗷嗷叫费观| 这个男人来自地球电影免费观看| 美国免费a级毛片| 欧美黄色片欧美黄色片| 欧美中文综合在线视频| 乱人伦中国视频| 色精品久久人妻99蜜桃| 如日韩欧美国产精品一区二区三区| 精品高清国产在线一区| 人人妻人人爽人人添夜夜欢视频| 多毛熟女@视频| 自线自在国产av| 一级,二级,三级黄色视频| 国产1区2区3区精品| 日韩欧美一区视频在线观看| 欧美黄色淫秽网站| 美女视频免费永久观看网站| 国产福利在线免费观看视频| 天天躁日日躁夜夜躁夜夜| 亚洲国产精品999| 狂野欧美激情性xxxx| 另类精品久久| 久久免费观看电影| 久久精品国产综合久久久| 国产91精品成人一区二区三区 | 亚洲国产欧美一区二区综合| 不卡av一区二区三区| 久久鲁丝午夜福利片| 亚洲精品国产一区二区精华液| 久久人人97超碰香蕉20202| 久久天躁狠狠躁夜夜2o2o | 男女边摸边吃奶| 欧美乱码精品一区二区三区| 黄色一级大片看看| 日本av手机在线免费观看| 精品一品国产午夜福利视频| 精品久久久久久电影网| 晚上一个人看的免费电影| 日韩人妻精品一区2区三区| 好男人视频免费观看在线| 久久av网站| 777米奇影视久久| 中文字幕制服av| a级毛片在线看网站| av电影中文网址| 亚洲精品美女久久av网站| 欧美国产精品一级二级三级| 亚洲国产看品久久| xxx大片免费视频| www.自偷自拍.com| 99热国产这里只有精品6| 狂野欧美激情性bbbbbb| 纵有疾风起免费观看全集完整版| 一本大道久久a久久精品| 捣出白浆h1v1| 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 少妇人妻久久综合中文| 欧美亚洲日本最大视频资源| 中国美女看黄片| 老司机深夜福利视频在线观看 | 极品人妻少妇av视频| 啦啦啦啦在线视频资源| 亚洲五月婷婷丁香| 亚洲精品第二区| 男女免费视频国产| 超碰成人久久| 肉色欧美久久久久久久蜜桃| av国产久精品久网站免费入址| 91精品三级在线观看| 一级毛片电影观看| 赤兔流量卡办理| 亚洲精品乱久久久久久| 免费黄频网站在线观看国产| 女人精品久久久久毛片| 乱人伦中国视频| 日本色播在线视频| 蜜桃国产av成人99| 后天国语完整版免费观看| 一二三四在线观看免费中文在| 精品免费久久久久久久清纯 | 亚洲一区二区三区欧美精品| 国产精品99久久99久久久不卡| 久久中文字幕一级| 欧美日韩国产mv在线观看视频| 精品少妇黑人巨大在线播放| 黄片播放在线免费| 人妻一区二区av| 久久久精品94久久精品| 赤兔流量卡办理| av网站在线播放免费| 久久久国产欧美日韩av| 在线av久久热| 极品少妇高潮喷水抽搐| 操美女的视频在线观看| 十八禁高潮呻吟视频| 免费在线观看黄色视频的| 精品国产超薄肉色丝袜足j| 久久天躁狠狠躁夜夜2o2o | 一区在线观看完整版| 久热爱精品视频在线9| 久久精品国产亚洲av高清一级| 大片免费播放器 马上看| 欧美激情极品国产一区二区三区| 一级毛片女人18水好多 | 好男人视频免费观看在线| 看免费成人av毛片| 黄色毛片三级朝国网站| 亚洲av日韩精品久久久久久密 | a级片在线免费高清观看视频| 亚洲国产中文字幕在线视频| 成年动漫av网址| 亚洲国产欧美在线一区| 国产又爽黄色视频| 日韩av免费高清视频| 国产成人系列免费观看| 777米奇影视久久| 免费黄频网站在线观看国产| av又黄又爽大尺度在线免费看| 精品亚洲成国产av| 亚洲精品乱久久久久久| 国产免费又黄又爽又色| 欧美日韩亚洲国产一区二区在线观看 | 欧美大码av| 精品一区二区三卡| 午夜激情久久久久久久| 男女无遮挡免费网站观看| 各种免费的搞黄视频| 国产亚洲av片在线观看秒播厂| 青草久久国产| 丰满饥渴人妻一区二区三| 又大又爽又粗| 啦啦啦啦在线视频资源| 日本午夜av视频| 免费在线观看影片大全网站 | 久热爱精品视频在线9| 国产精品久久久久久人妻精品电影 | 成人国产av品久久久| 欧美日韩av久久| 亚洲av日韩在线播放| 日本欧美国产在线视频| 国产精品久久久久久人妻精品电影 | 国产成人91sexporn| 成人国产一区最新在线观看 | 日韩精品免费视频一区二区三区| √禁漫天堂资源中文www| 精品一品国产午夜福利视频| 国产片内射在线| 80岁老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 人妻一区二区av| 黑人欧美特级aaaaaa片| 深夜精品福利| 一区二区三区四区激情视频| 五月开心婷婷网| 啦啦啦视频在线资源免费观看| 国产欧美日韩精品亚洲av| 黄片播放在线免费| 九草在线视频观看| 国产免费一区二区三区四区乱码| 最近中文字幕2019免费版| 国产免费一区二区三区四区乱码| 久久中文字幕一级| 国产片特级美女逼逼视频| 99国产精品一区二区蜜桃av | 国产色视频综合| 国产亚洲一区二区精品| 午夜av观看不卡| 飞空精品影院首页| 亚洲国产av新网站| 国产一区有黄有色的免费视频| 国产精品人妻久久久影院| 搡老岳熟女国产| 成人三级做爰电影| 国产精品久久久人人做人人爽| 国产视频首页在线观看| 精品少妇黑人巨大在线播放| 亚洲中文日韩欧美视频| 免费女性裸体啪啪无遮挡网站| 免费在线观看视频国产中文字幕亚洲 | 18在线观看网站| 日韩中文字幕欧美一区二区 | 男女午夜视频在线观看| 桃花免费在线播放| 亚洲国产看品久久| 久久久久精品国产欧美久久久 | 久久中文字幕一级| a级毛片黄视频| a级毛片在线看网站| 亚洲五月色婷婷综合| 免费人妻精品一区二区三区视频| 老司机靠b影院| 欧美日韩福利视频一区二区| videos熟女内射| 大码成人一级视频| 亚洲精品久久成人aⅴ小说| 欧美国产精品一级二级三级| 高清视频免费观看一区二区| 久久 成人 亚洲| 看十八女毛片水多多多| 国产精品久久久久成人av| 国产成人精品无人区| 久久久久国产精品人妻一区二区| 午夜av观看不卡| 成年av动漫网址| 国产成人91sexporn| 国产成人系列免费观看| av一本久久久久| 国产91精品成人一区二区三区 | 国产精品一区二区在线不卡| 精品人妻熟女毛片av久久网站| 在线观看免费高清a一片| 亚洲九九香蕉| 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说| 亚洲精品中文字幕在线视频| 日本vs欧美在线观看视频| 国产日韩欧美在线精品| av天堂久久9| 日本黄色日本黄色录像| 亚洲国产欧美网| 欧美激情高清一区二区三区| 一级a爱视频在线免费观看| 两性夫妻黄色片| 成年av动漫网址| 国产精品亚洲av一区麻豆| 精品一区二区三卡| 亚洲情色 制服丝袜| 19禁男女啪啪无遮挡网站| 99九九在线精品视频| 最近最新中文字幕大全免费视频 | 国产精品秋霞免费鲁丝片| 一个人免费看片子| 夫妻午夜视频| 捣出白浆h1v1| 一本大道久久a久久精品| 母亲3免费完整高清在线观看| 韩国高清视频一区二区三区| 久久这里只有精品19| 国产又爽黄色视频| 欧美久久黑人一区二区| 国产男女超爽视频在线观看| 中文乱码字字幕精品一区二区三区| 日韩一本色道免费dvd| 妹子高潮喷水视频| 国产麻豆69| 国产又爽黄色视频| 日本五十路高清| 成年女人毛片免费观看观看9 | 成人三级做爰电影| 人人妻人人爽人人添夜夜欢视频| 精品卡一卡二卡四卡免费| 精品免费久久久久久久清纯 | av片东京热男人的天堂| 午夜免费观看性视频| 亚洲av电影在线观看一区二区三区| 啦啦啦 在线观看视频| 亚洲欧美清纯卡通| 汤姆久久久久久久影院中文字幕| 久久久久久久大尺度免费视频| 操美女的视频在线观看| 国产午夜精品一二区理论片| 国产1区2区3区精品| 天堂中文最新版在线下载| 国产亚洲av高清不卡| 婷婷色综合大香蕉| 天天影视国产精品| 久久人人爽av亚洲精品天堂| 精品久久久久久久毛片微露脸 | 中文字幕亚洲精品专区| 亚洲国产欧美在线一区| 国产精品二区激情视频| 咕卡用的链子| 日韩av不卡免费在线播放| 久久久久久久大尺度免费视频| 日韩大片免费观看网站| 久久天堂一区二区三区四区| 丝袜美足系列| 久久久亚洲精品成人影院| 亚洲国产欧美日韩在线播放| 在线av久久热| 免费在线观看完整版高清| 王馨瑶露胸无遮挡在线观看| 日韩av不卡免费在线播放| 国产欧美日韩一区二区三区在线| 国产一区亚洲一区在线观看| 免费在线观看完整版高清| 好男人视频免费观看在线| 母亲3免费完整高清在线观看| 亚洲国产欧美日韩在线播放| 亚洲国产欧美在线一区| 亚洲精品美女久久久久99蜜臀 | 十八禁人妻一区二区| www.自偷自拍.com| 欧美少妇被猛烈插入视频| 大香蕉久久网| 丝袜脚勾引网站| 日本色播在线视频| 日韩视频在线欧美| 亚洲 欧美一区二区三区| 欧美xxⅹ黑人| 久久久久久免费高清国产稀缺| 日本vs欧美在线观看视频| 久久毛片免费看一区二区三区| 黄色毛片三级朝国网站| 男的添女的下面高潮视频| 午夜视频精品福利| 一边亲一边摸免费视频| 亚洲伊人色综图| svipshipincom国产片| 女性生殖器流出的白浆| 国产精品久久久久久精品电影小说| 国产熟女午夜一区二区三区| 国产在线观看jvid| 高潮久久久久久久久久久不卡| 国产亚洲欧美在线一区二区| 水蜜桃什么品种好| 女警被强在线播放| 赤兔流量卡办理| 精品卡一卡二卡四卡免费| 亚洲精品自拍成人| 亚洲精品一区蜜桃| 伦理电影免费视频| 亚洲九九香蕉| 免费看av在线观看网站| 国产男人的电影天堂91| 日韩,欧美,国产一区二区三区| 免费日韩欧美在线观看| 最黄视频免费看| 日韩 欧美 亚洲 中文字幕| 亚洲色图综合在线观看| 美女大奶头黄色视频| 欧美+亚洲+日韩+国产| 激情五月婷婷亚洲| 大陆偷拍与自拍| 国产精品国产三级国产专区5o| 国产97色在线日韩免费| 国产精品二区激情视频| 国产老妇伦熟女老妇高清| 在线观看免费午夜福利视频| 久久久精品免费免费高清| 好男人视频免费观看在线| 十八禁人妻一区二区| 亚洲成人免费av在线播放| 丰满人妻熟妇乱又伦精品不卡| 免费高清在线观看日韩| 国产精品一二三区在线看| 狂野欧美激情性xxxx| 日本黄色日本黄色录像| 午夜激情av网站| 亚洲av成人精品一二三区| 国产又色又爽无遮挡免| 女人爽到高潮嗷嗷叫在线视频| 欧美亚洲 丝袜 人妻 在线| 一个人免费看片子| av在线播放精品| 亚洲国产欧美一区二区综合| 日本五十路高清| 天堂8中文在线网| 亚洲国产精品一区二区三区在线| 999久久久国产精品视频| 激情五月婷婷亚洲| 欧美激情高清一区二区三区| 悠悠久久av| 欧美另类一区| 黑丝袜美女国产一区| 在线观看www视频免费| 大码成人一级视频| 成人免费观看视频高清| 日韩一区二区三区影片| 久久久久国产一级毛片高清牌| 大香蕉久久网| 欧美成人午夜精品| 精品欧美一区二区三区在线| 成人免费观看视频高清| 亚洲中文字幕日韩| 777米奇影视久久| 久久久国产精品麻豆| 久久人人爽人人片av| 欧美黄色淫秽网站| 岛国毛片在线播放| 国产日韩欧美视频二区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av日韩精品久久久久久密 | 国产精品 国内视频| 日本五十路高清| 校园人妻丝袜中文字幕| 亚洲av欧美aⅴ国产| 国产日韩欧美亚洲二区| 亚洲 欧美一区二区三区| 1024视频免费在线观看| 国产成人系列免费观看| 99国产精品99久久久久| 青草久久国产| 午夜福利乱码中文字幕| 天天躁日日躁夜夜躁夜夜| 国产熟女午夜一区二区三区| 国产精品久久久人人做人人爽| 欧美人与性动交α欧美软件| 日本一区二区免费在线视频| 国产伦理片在线播放av一区| 亚洲国产精品成人久久小说| e午夜精品久久久久久久| 在线天堂中文资源库| 亚洲欧美日韩另类电影网站| 可以免费在线观看a视频的电影网站| 女性生殖器流出的白浆|