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

    人工魚礁建設(shè)對福建霞浦海域營養(yǎng)鹽輸運的影響

    2016-07-06 10:30:07上海海洋大學(xué)海洋科學(xué)學(xué)院上海201306
    海洋科學(xué) 2016年2期
    關(guān)鍵詞:上升流礁體魚礁

    肖 榮, 楊 紅(上海海洋大學(xué) 海洋科學(xué)學(xué)院, 上海 201306)

    ?

    人工魚礁建設(shè)對福建霞浦海域營養(yǎng)鹽輸運的影響

    肖 榮, 楊 紅
    (上海海洋大學(xué) 海洋科學(xué)學(xué)院, 上海 201306)

    根據(jù)2012年—2013年對霞浦人工魚礁海域調(diào)查所獲得的實驗數(shù)據(jù), 應(yīng)用CFD軟件, 對礁區(qū)投放的三類不同結(jié)構(gòu)單位的魚礁組合在非定常流作用下的三維流場進行了數(shù)值模擬試驗, 分析了各魚礁組合之間的流場差異, 并估算了諸營養(yǎng)鹽垂直輸送通量。研究表明: 三種魚礁迎流面上升流最大高度為礁高的2.27~3.73倍, 上升流最大速度與來流速度之比為0.40~0.72, 上升流平均速度與來流速度之比為0.10~0.17; 計算域內(nèi)單位魚礁組合上升流平均速度為10–2cm/s量級, 三種單位魚礁組合區(qū)域營養(yǎng)鹽垂直通量大小關(guān)系為: 方型魚礁>船礁>梯形臺魚礁, 方型魚礁區(qū)PO-P、NO-N、SiO-Si的平均垂直通量分別為438.3 mg/(m2·d)、2212.5 mg/(m2·d)、7288.3 mg/(m2·d), 魚礁投放會改善該海域的水動力環(huán)境, 形成典型的上升流區(qū), 海域初級生產(chǎn)力得以明顯提升。

    人工魚礁; 流場效應(yīng); 數(shù)值模擬; 上升流; 營養(yǎng)鹽垂直輸運通量

    陸架海域的人工魚礁投放后, 受周期性潮流影響, 其周圍水體的流場隨之發(fā)生變化, 產(chǎn)生上升流和背渦流。陳勇等[1]研究表明, 上升流能夠?qū)⒑5椎某练e物和營養(yǎng)鹽帶到表層, 加快營養(yǎng)物質(zhì)循環(huán)速度,提高海域的基礎(chǔ)餌料水平, 從而提高海區(qū)初級生產(chǎn)力。人工魚礁產(chǎn)生的上升流和背渦流的強度可表征魚礁對水環(huán)境的物理作用的程度, 故可表征魚礁生態(tài)效果的優(yōu)劣。從海水化學(xué)要素分布來看, 上升流直接表現(xiàn)在營養(yǎng)鹽的區(qū)域性高值, 因而, 海水中氮、磷、硅營養(yǎng)鹽的分布成為上升流的化學(xué)指標之一[2]。目前關(guān)于人工魚礁的研究主要涉及生物和非生物因子[3]、生產(chǎn)力[4]、動植物區(qū)系時空變化[5]、三大效益[6]、選址及布局[7]、休閑漁業(yè)[8]、增殖資源[9]等角度進行研究, 對于有關(guān)人工魚礁建設(shè)產(chǎn)生的上升流對海域營養(yǎng)鹽輸運影響的研究鮮有報道[10], 而弄清和掌握上升流強度及分布對營養(yǎng)鹽含量及輸運的影響, 對科學(xué)的估算和預(yù)報投礁海域的生產(chǎn)力和營養(yǎng)現(xiàn)狀, 以及海洋環(huán)境保護, 綜合利用和管理意義重大[11]。本文通過對2012年、2013年霞浦人工魚礁投放海域的調(diào)查, 應(yīng)用計算流體力學(xué)和水環(huán)境化學(xué)相關(guān)原理, 研究霞浦投礁區(qū)三類單位魚礁組合上升流的分布及強度, 在數(shù)值模擬試驗的基礎(chǔ)上探討魚礁投放海域上升流的垂直速度大小, 估算諸營養(yǎng)鹽垂直輸送通量。

    1 材料和方法

    1.1 三維數(shù)值模型

    根據(jù)礁區(qū)建設(shè)的目標功能, 結(jié)合礁區(qū)的地形、水深和水文條件, 在完成適合該礁區(qū)的魚礁備選單體設(shè)計、并初步提出配置組合方案后, 在數(shù)值模擬計算區(qū)域里對不同魚礁組合配置進行優(yōu)化。本試驗選用在福建近岸海域人工魚礁區(qū)應(yīng)用廣泛、穩(wěn)定性良好的中空結(jié)構(gòu)梯形臺礁體、中空結(jié)構(gòu)方型礁體和船礁,對其礁體組合進行模擬研究。首先, 對魚礁單體進行CAD三維造型, 導(dǎo)入CAD文件到網(wǎng)格剖分軟件GAMBIT, 并復(fù)制多個單礁體、根據(jù)魚礁組合方案排列, 在GAMBIT中進行三維網(wǎng)格剖分, 采用三維非結(jié)構(gòu)化四面體網(wǎng)格; 其次, 將剖分后的網(wǎng)格數(shù)據(jù)導(dǎo)入FLUENT軟件, 在靠近礁體周圍進行網(wǎng)格加密處理, 湍流模型選用大渦模擬(large eddy simulation, LES)模型, 數(shù)值模擬不同單位魚礁組合的流場效應(yīng),以上升流區(qū)的規(guī)模和強度為主要指標, 對不同類型魚礁組合數(shù)值模擬結(jié)果進行分析, 據(jù)此選擇確定適合海域魚礁建設(shè)的配置組合方案。數(shù)值模擬技術(shù)采用三維雙精度數(shù)(3ddp)和時變(unsteady)模型、LES湍流模式, 求解器(solver)選用基于壓力隱式求解(pressure based), 流相采用常規(guī)液態(tài)水單項, 考慮重力作用的影響。數(shù)值計算區(qū)域來流方向設(shè)為流速入口(velocity-inlet), 計算區(qū)域出口設(shè)為自由出流條件(outflow), 計算區(qū)域的兩個側(cè)面和頂面均設(shè)為對稱邊界條件(symmetry); 計算區(qū)域的底面和魚礁壁面均設(shè)為無滑移邊界條件(wall)。

    湍流模型選擇是影響仿真精度的關(guān)鍵, 本文選用的湍流模式為大渦模擬法。大渦模擬是介于直接數(shù)值模擬(DNS)與Reynolds 平均法(RANS)之間的一種湍流數(shù)值模擬方法, 即用瞬時的N-S方法直接模擬湍流中的大尺度渦, 不直接模擬小尺度渦, 而小渦對大渦的影響通過近似的模型來考慮, 由此能模擬出所有大于網(wǎng)格尺度的渦的運動[12]。所用礁體模型的網(wǎng)格剖分精度可達數(shù)十厘米級, 足以達到對礁體的分辨率。為了計算效率和模擬精度并重, 一般在保證既能分辨礁體、又有足夠大的計算區(qū)域的前提下, 單體模型的網(wǎng)格數(shù)可控制在100萬個左右, 礁體組合的網(wǎng)格數(shù)在200萬~400萬個左右。以梯形臺魚礁計算域為例, 其縱向中軸斷面的網(wǎng)格劃分情況見圖1。

    圖1 梯形臺魚礁縱向中軸斷面網(wǎng)格剖面圖Fig. 1 The grid on the center longitudinal section of trapezoidal reef

    1.2 數(shù)值實驗

    基于流場效應(yīng)數(shù)值仿真研究人工魚礁不同礁型組合情況下的上升流強度, 實驗安排了三組數(shù)值,其一用于判定不同魚礁類型相同組合情況下上升流分布的差異, 其二用于比較典型鋼筋砼結(jié)構(gòu)礁型與船礁的流場效應(yīng)差異, 從而判定三個魚礁組之間上升流分布的差異, 為不同魚礁區(qū)營養(yǎng)鹽的垂直輸運提供流場基礎(chǔ)。以福建三沙根竹仔礁區(qū)建設(shè)中采用的A、B、C共三種單位魚礁組合為例來分析比較不同單位魚礁配置布設(shè)的差異和優(yōu)劣。A、B、C 3種單位魚礁組合的礁體個數(shù)分別為16、16、9, 組合魚礁礁距均為5m, 均勻分布。

    16礁組合: 分為中空結(jié)構(gòu)梯形臺魚礁(A型)和中空結(jié)構(gòu)方型魚礁(B型)兩組, 梯形臺礁體外形尺寸為3.0 m×3.0 m×3.0 m, 方型礁體外形尺寸為3.0 m× 3.0 m×3.0 m, 圖2為A、B、C3種魚礁模型的效果圖。魚礁的計算區(qū)域如圖3所示。來流方向從左至右, 坐標原點位于計算域進口底面中點, 計算域進口至魚礁組中心距離為28 m, 計算域的兩側(cè)面至魚礁組中心的距離為45 m, 整個計算域為90 m×90 m×15 m。礁體材料為鋼筋混凝土。

    9礁組合: 以廢舊漁船(C型)作為設(shè)置礁型, 魚礁效果圖見圖2, 計算域進口至船礁組中心距離為20 m,計算域的兩側(cè)面至船礁組中心的距離為45 m, 整個計算域為90 m×90 m×15 m。礁體材料為玻璃鋼質(zhì)。

    三種單位魚礁組合示意圖見圖4, 數(shù)值計算區(qū)域水深取為福建近岸魚礁投放區(qū)的典型水深15 m, 來流速度取島礁海域人工魚礁區(qū)的典型大潮流速0.5 m/s。

    1.3 營養(yǎng)鹽垂直通量的計算

    圖2 人工魚礁模型效果圖Fig. 2 Arendering of reef model

    圖3 魚礁計算區(qū)域Fig. 3 Reef model in the computational domain

    圖4 3種單位魚礁組合示意圖Fig. 4 Sketch of 3 types of artificial reef assemblages

    營養(yǎng)鹽的垂直通量指垂直方向營養(yǎng)鹽自底層向表層在單位時間內(nèi)通過單位面積的量, 即底層水體進入表層水體的營養(yǎng)鹽含量[14]。

    式中, α為單位換算常數(shù), C值為研究海區(qū)營養(yǎng)鹽平均濃度值(μmol/mL); ω為上升流上升速度。

    2 結(jié)果與分析

    2.1 單位魚礁組合流場效應(yīng)

    項目用海平面布置考慮施工投放的方便及礁群的互補、互助性, 人工魚礁平面布置為網(wǎng)格狀分組穿插布置, 除有利于增大礁體在礁區(qū)內(nèi)的迎流面積外,更有利于在礁區(qū)內(nèi)產(chǎn)生較多的渦流效應(yīng), 將礁區(qū)內(nèi)的不良物質(zhì)帶走, 補充礁區(qū)內(nèi)的營養(yǎng)物質(zhì), 調(diào)節(jié)礁區(qū)內(nèi)生態(tài)環(huán)境, 所形成的上升流將海底深層的營養(yǎng)鹽類帶到光照充足的表水面, 促進浮游生物的繁殖, 從而提高海洋的初級生產(chǎn)力。三沙根竹仔礁區(qū)礁體的設(shè)計應(yīng)充分考慮海區(qū)物種多樣性的需要, 魚礁材料和魚礁造型的選用和配置也應(yīng)考慮生物附著率的高低, 為其創(chuàng)造良好的棲息生長和繁衍場所。根據(jù)投礁區(qū)域地形、水深等條件, 為了在有限的投資下增加魚礁組的密集度, 3種單位魚礁組合均勻布置, 縱向每個魚礁組中心間距90 m, 橫向間距90 m, 形成網(wǎng)格狀錯開布置。

    無論是誘集漁業(yè)資源還是通過生態(tài)調(diào)控增殖漁業(yè)資源的人工魚礁建設(shè)項目, 首先是要讓魚類感知礁體的存在, 以及礁體存在所帶來的環(huán)境變異和餌料條件等差異。魚類對礁體的感知很大程度上是通過感知流場環(huán)境的變異(產(chǎn)生上升流和背渦流)而達成的, 流場環(huán)境的變異導(dǎo)致的營養(yǎng)鹽向上層水體的輸運將增加礁區(qū)餌料生物的供給量, 此外對于趨礁性魚類還通過內(nèi)部結(jié)構(gòu)的復(fù)雜化來達到增大棲息地空間和附著生物量, 從而提高誘集和增殖效果。

    基于上述數(shù)值方法, 模擬得到了3種單位魚礁組合的三維繞流流場, 圖5為不同魚礁組合y=0、y=4(過礁體中心)縱向斷面上升流速分布圖。取水流的垂向速度分量與來流速度V0之比大于或等于5%的水域作為上升流區(qū)域[15]。上升流規(guī)模用上升流的最大高度Ls與礁高H之比加以衡量, 而上升流強度用最大上升流速Vmax與來流速度V0之比和上升流域內(nèi)的上升流均值Vave與來流速度V0之比加以衡量。

    圖5 魚礁縱向斷面上升流速分布圖Fig. 5 Velocity contour on the center longitudinal section of the reef

    本試驗結(jié)果表明(表1): 在相同來流速度下, 3種單位魚礁組合上升流高度為礁高的2.3~3.7倍, 略大于張碩等[16]的研究結(jié)果(1.5~2.8倍), 其原因可能是由于試驗條件和尺度比不同所致。本試驗在進行數(shù)值計算時是以魚礁的實際尺寸來模擬, 采用的是投礁海域?qū)嶋H流速, 相比張碩[16]采用的礁型大得多,而且魚礁組合有交互效應(yīng)產(chǎn)生, 因而產(chǎn)生的上升流礁高比較大。3種單位魚礁組合上升流最大速度與來流速度之比為0.40~0.72, 上升流平均速度與來流速度之比為0.10~0.17, 與黃遠東等[17]的研究結(jié)果基本吻合。

    表1 不同組合情況下魚礁模型上升流規(guī)模和強度Tab. 1 The distribution and intensity of the upwelling in different combinations of reefs

    A型、B型單位魚礁組合為透水礁, 其絕對迎流面積(扣除礁迎流面中空部分)較小, 而C型單位魚礁組合為不透水礁, 其絕對迎流面面積較大, 因此, 同等程度下A型、B型單位魚礁組合上升流相對面積較C型要小。這也進一步佐證了魚礁透水性能對其流場變化, 尤其是對上升流影響區(qū)域的大小有直接的關(guān)系。A型、B型魚礁在不同來流方向上, 魚礁的迎流面積相對變化不大, 而C型魚礁長軸和短軸差異較大, 本文模擬的是長軸與來流方向垂直的情況,此時C型魚礁的絕對迎流面積達到最大, 隨著漲落潮往復(fù)流的變化, C型魚礁迎流面積在短軸與來流方向垂直時與本文模擬的情況之間變化, 相對來說C型魚礁的不同走向會對上升流產(chǎn)生明顯的影響, 單位迎流面積的C型魚礁產(chǎn)生的上升流相對面積隨著來流方向的改變會產(chǎn)生明顯的變化, 流場的穩(wěn)定性不如A型、B型魚礁。

    A型魚礁迎流面坡度的存在對自由來流的加速起到了一定的緩解作用(相對垂直放置的魚礁迎流面), 也表明A型魚礁的阻流作用要稍弱于B型魚礁; 而A型魚礁繞流的上升流平均速度與來流速度之比以及上升流最大速度與來流速度之比略大于B型魚礁。

    對于單位魚礁組合來說, 上升流主要形成在第一排礁體的迎流面, 魚礁組合之間的交互作用, 使得第二排礁體對第一排礁體產(chǎn)生影響, 由于A型、B型魚礁為透水礁, 水流通過第一排魚礁時形成一個緩流區(qū)域, A型、B型魚礁第二排礁體迎流面產(chǎn)生的上升流要比C型魚礁明顯, 隨著水流通過魚礁數(shù)量的增加, 魚礁迎流面產(chǎn)生的上升流逐漸減少, 整個計算域上升流速度取決于每個魚礁產(chǎn)生的上升流的綜合效應(yīng)。

    2.2 營養(yǎng)鹽輸運

    人工魚礁的流場效應(yīng)深刻影響著魚礁的物理環(huán)境功能及生態(tài)效應(yīng)的發(fā)揮。通常除碎浪帶外, 沿岸海域水體的垂向運動相對水平運動而言往往可以忽略。如果在潮流主流軸方向上投放人工魚礁, 可以生成很強的局部上升流, 其量值可以與水平流相當[6],從而促進表底層水體交換。通過這種水體的垂直交換功能, 上升流不斷將底層、近底層低溫、高鹽富營養(yǎng)的海水涌升至表層, 導(dǎo)致溫度、鹽度格局重新分布,使水文條件更適合于中、上層魚類棲息和集群活動的要求。另外, 餌料浮游生物高密度區(qū)主要出現(xiàn)在上升流區(qū), 這為中心漁場的形成創(chuàng)造了必要的條件[18]。

    人工魚礁產(chǎn)生的上升流將底層水帶入真光層的營養(yǎng)鹽含量可通過上升流海域營養(yǎng)鹽的垂直輸送通量計算。對于NO-N、PO-P、SiO-Si: (1)式中單位換算常數(shù)分別取αN=12096, αP=26758, αSi= 24270[19]; A型、B型、C型單位魚礁組合上升流平均流速ω取自CFD模擬值25.8×10–3、45.5×10–3、41.9×10–3cm/s, 東海上升流ω量級為10–3cm/s[20-23]、臺灣海峽上升流平均流速ω為4.4×10–3cm/s[24], 本文是以單位魚礁組合為研究對象, 流場調(diào)控距離取值為單位魚礁組合中心間距90 m, 根據(jù)人工魚礁區(qū)周邊水質(zhì)、生物資源等實際調(diào)查的結(jié)果表明, 人工魚礁區(qū)的影響范圍可達到其建設(shè)規(guī)模的近5~10倍之多,通常的面積范圍在100~101m量級的公頃數(shù)[25], 實際海域由于魚礁設(shè)置所產(chǎn)生的流場影響范圍在水平尺度上一般不超過魚礁規(guī)模的50倍[26], 所以數(shù)值計算得到的計算域內(nèi)上升流的平均速度要大于自然環(huán)境海域存在的上升流。營養(yǎng)鹽濃度C取2013年11月調(diào)查站位諸營養(yǎng)鹽含量平均值。

    表2 幾個上升流海域諸營養(yǎng)鹽的垂直輸送通量比較Tab. 2 The vertical transportation flux of nutrients in several upwelling areas

    3 結(jié)論

    基于ANSYS FLUENT12.1平臺, 數(shù)值模擬了來流速度為0.5m/s下的A型、B型、C型單位魚礁組合周圍的三維流場。通過分析過礁體中心的鉛垂平面上速度分布, 得出: (1)A型、B型單位魚礁組合上升流的最大高度與礁高之比為2.27~2.41、上升流最大速度與來流速度的比約為0.40~0.43, 上升流區(qū)域平均速度與來流速度的比約為0.10~0.13, C型魚礁最大高度與礁高之比為3.73、上升流最大速度與來流速度的比約為0.72, 上升流區(qū)域平均速度與來流速度的比約為0.17; (2)A型、B型魚礁外形尺寸相當、具有類似的空方比, 采用了相同的布設(shè)方案, 單位魚礁組合產(chǎn)生的上升流規(guī)模和強度差異不大; (3)C型單位魚礁組合為不透水礁, 其絕對迎流面面積較大, 因此, 同等程度下C型單位魚礁組合上升流相對面積較A型、B型要大; (4)單位魚礁組合流場效應(yīng)主要取決于礁體數(shù)量、礁體空方比、礁體迎流面積。

    計算域內(nèi)單位魚礁組合上升流平均速度為10–2cm/s量級, 相比自然海域上升流的速度要高, 主要由于本文數(shù)值模擬試驗以單位魚礁組合為基本研究對象, 水平尺度范圍為單位魚礁組合中心間距(90m), 實際海域單位魚礁組合流場調(diào)控范圍可以達到更大距離。通過對營養(yǎng)鹽的垂直輸運通量估算可知, 不管投放何種魚礁都能產(chǎn)生明顯的上升流, 營養(yǎng)鹽通量相比沿岸上升流處于較高的水平, 魚礁投放會改善該海域的水動力環(huán)境, 形成典型的上升流區(qū), 海域初級生產(chǎn)力得以明顯提升。今后, 可以在單位魚礁組合的基礎(chǔ)上來研究整個投礁海域生態(tài)調(diào)控范圍內(nèi)的流場效應(yīng), 評估礁區(qū)海域營養(yǎng)鹽的分布情況及水質(zhì)營養(yǎng)水平, 從定量的角度研究礁區(qū)營養(yǎng)鹽的垂直輸送通量。

    [1] 陳勇, 于長清, 張國勝, 等. 人工魚礁的環(huán)境功能與集魚效果[J]. 大連水產(chǎn)學(xué)院學(xué)報, 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 Fisheries University, 2002, 17 (1): 64-69.

    [2] Hung T C, Han B C, Hsu D S. Chemical and biomass studies: (2) Biological activities in upwelling off the southwestern coast of Taiwan including the Penghe area[J]. Acta Oceanogiaphics Taiwanica, 1987, (18): 62-74.

    [3] Haroun R J, Gomez M, Hernandez J J, et al. Environ-ment description of an artificial reef site in Gran Ganaria (Canary Islands, Spain) prior to reef placement [J]. Bulletin of Marine Science, 1994, 5: 932-938.

    [4] DeMartini E E, Barnett A M, Johnson T D, et al. Growth and production estimates for biomass-dominant fishes on a southern California artificial reef [J]. Bulletin of Marine Science, 1994, 55: 484-500.

    [5] Bortone S A, Van Tassell J, Brito A, et al. Enhancement of the nearshore fish assemblage in the Canary Islands with artificial habitats [J]. Bulletin of Marine Science, 1994, 55: 602-608.

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

    Lin Jun, Zhang Shouyu. Research advances on physical stability and ecological effectsof artificial reef[J]. Marine Fisheries, 2006, 28(3): 257-262.

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

    Zhao Haitao, Zhang Yifei, Hao Chunling, et al. Sitting and designing of artificial reefs[J]. Journal Of Marine Sciences, 2006, 24(4): 69-76.

    [8] 王偉定, 徐漢祥, 潘國良, 等. 浙江沿岸休閑生態(tài)型人工魚礁建設(shè)現(xiàn)狀與展望[J]. 浙江海洋學(xué)院學(xué)報(自然科學(xué)版), 2007, 26(1): 22-27.

    Wan Weiding, Xu Hanxiang, Pan Guoliang, et al. Current situation and prospect of recreational ecotypic artificial reef construction in Zhejiang[J]. Journal of Zhejiang Ocean University(Natural Science), 2007, 26(1): 22-27.

    [9] 段玉彤, 李志超. 利用人工魚礁增值河北漁業(yè)資源的探討[J]. 河北漁業(yè), 2007, 4: 43-44.

    Duan Yutong, Li Zhichao. Discussion on fishery resources multiplication of artificial in Hebei reefs[J]. Hebei Fisheries, 2007, 4: 43-44.

    [10] 章守宇, 張煥君, 焦俊鵬, 等. 海州灣人工魚礁海域生態(tài)環(huán)境的變化[J]. 水產(chǎn)學(xué)報, 2006, 30(4): 475-480.

    Zhang Shouyu, Zhang Huanjun, Jiao Junpeng, et al. Change of ecological environment of artificial reef waters in Haizhou Bay[J]. Journal of Fisheries of China, 2006, 30(4): 475-480.

    [11] 王珊珊, 曹志敏, 藍東兆, 等. 福建省興化灣灣口海域營養(yǎng)鹽的分布狀況及其影響因素[J]. 環(huán)境化學(xué), 2006, 25(2): 234-235.

    Wang Shanshan, Cao Zhimin, Lan Dongzhao, et al. The nutrients distribution and its influencing factors in Fujian Xinghua Bay[J]. Environmental Chemistry, 2006, 25(2): 234-235.

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

    Wang Fujun. Computational fluid dynamics analysis[M]. Beijing: Tsinghua University Press, 2004.

    [13] GB17387. 4-2007, 海洋監(jiān)測規(guī)范[S]. GB17387. 4-2007, The specification for marine monitoring[S].

    [14] 雷鵬飛. 浙江近海上升流速度及其營養(yǎng)鹽通量計算[J].海洋湖沼通報, 1984, (2): 22-26.

    Lei Pengfei. Evalution of upwelling rate and the nutrient flux in the zone of Zhejiang[J]. Transactions of Oceanology and Limnology, 1984, (2): 22-26.

    [15] 黃遠東, 姜劍偉, 趙樹夫. 方型人工魚礁周圍水流運動的數(shù)值模擬研究[J]. 水資源與水工程學(xué)報, 2012, 23(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, 23(3): 1-3.

    [16] 張碩, 孫滿昌, 陳勇. 不同高度混凝土模型礁上升流特性的定量研究[J]. 大連水產(chǎn)學(xué)院學(xué)報, 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 Fisheries University, 2008, 23(5): 353-358.

    [17] 黃遠東, 趙樹夫, 姜劍偉, 等. 多孔方型人工魚礁繞流的數(shù)值模擬研究[J]. 水資源與水工程學(xué)報, 2012, 23(5): 15-18.

    Huang Yuandong, Zhao Shufu, Jiang Jianwei, et al. Study on numerical simulation of water flow past a porous square artificial reef[J]. Journal of Water Resources and Water Engineering, 2012, 23(5): 15-18.

    [18] 蔡德華. 汕頭-臺灣淺灘漁場漁業(yè)自然條件及上升流區(qū)特征[J]. 水產(chǎn)科技, 2007(3): 22-23.

    Cai Dehua. Fishery natural conditions and characteristics of upwelling region in Shantou-Taiwan Bank[J]. Fisheries Technology, 2007(3): 22-23.

    [19] 陳水土, 阮五崎. 臺灣海峽上升流區(qū)氮、磷、硅的化學(xué)特性及輸送通量估算[J]. 海洋學(xué)報, 1996, 18(3) : 36- 44.

    Chen Shuitu, Ruan Wuqi. The chemical characteristics and transportation flux estimation of nitrogen, phosphorus and silicon in Taiwan strait upwelling area[J]. Acta Oceanologica Sinica, 1996, 18(3) : 36- 44.

    [20] 袁耀初. 東海三維海流的一個預(yù)報模式-黑潮調(diào)查研究論文選(五) [R]. 海洋出版社, 1993, 311-323.

    Yuan Yaochu. A three-dimensional current of the East China Sea forecast model—Kuroshio survey research papers election(五)[R].Ocean Press, 1993, 311-323.

    [21] 羅義勇. 東海沿岸上升流的數(shù)值計算[J]. 海洋湖沼通報, 1998, 3: 1-5.

    Luo Yiyong. Numerical modelling of upwelling in coastal areas of the East China Sea[J], Transactions of Oceanology and Limnology , 1998, 3: 1-5.

    [22] 羅義勇, 俞光耀. 風(fēng)和臺灣暖流引起東海沿岸上升流數(shù)值計算[J]. 青島海洋大學(xué)學(xué)報, 1998, 28(4): 536-542.

    Luo Yiyong, Yu Guangyao. Numerical studies of wind and twc- driven upwelling in coastal areas of the East China Sea[J]. Journal of Ocean University of Qingdao, 1998, 28(4): 536-542.

    [23] 朱建榮. 夏季長江口外水下河谷西側(cè)上升流產(chǎn)生的動力機制[J]. 科學(xué)通報, 2003, 48(23): 2488-2492.

    Zhu Jianrong. The produce dynamic mechanism of summer upwelling on the west side of the Yangtze river estuary underwater valley[J]. Chinese Science Bulletin, 2003, 48(23): 2488-2492.

    [24] 陳水土, 阮五琦. 臺灣海峽中北部夏季溶解氧飽和度分布與上升流的關(guān)系[J]. 臺灣海峽, 1991, 10(l): 16-24.

    Chen Shuitu, Ruan Wuqi. Summer upwelling and saturation dissolved oxygen distribution in the Central and Northern Taiwan Strait[J]. Journal of Oceanography in Taiwan Strait, 1991, 10(l): 16-24.

    [25] 中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所. 中國人工魚礁的理論與實踐[M]. 廣州: 廣東科技出版社, 2005.

    South China sea fisheries research institute, Chinese academy of fishery sciences. Theory and practice of Chinese artificial reefs[M]. Guangzhou: Guangdong Science and Technology Press, 2005.

    [26] 潘靈芝, 林軍, 章守宇. 鉛直二維定常流中人工魚礁流場效應(yīng)的數(shù)值實驗[J]. 上海水產(chǎn)大學(xué)學(xué)報, 2005, 14 (4) : 406 - 412.

    Pang Lingzhi, Lin Jun, Zhang Shouyu. A numerical experiment of the effects of artificial reef on vertical 2-dimensional steady flow field[J]. Journal of Shanghai Ocean University, 2005, 14(4): 406-412.

    [27] 吳麗云, 阮五琦. 閩南-臺灣淺灘漁場上升流區(qū)營養(yǎng)鹽的研究. 閩南-臺灣淺灘漁場上升流區(qū)生態(tài)系研究[M].北京: 科學(xué)出版社, 1991, 169-178.

    Wu Liyun, Ruan Wuqi. Study on the nutrient in the upwelling region in Minnan-Taiwan Bank fishing ground, Study of the ecosystems in the upwelling region in Minnan-Taiwan Bank fishing ground[M]. Beijing: Science Press, 1991, 169-178.

    [28] 王桂云. 東海海洋化學(xué)要素含量及其分布.黑潮調(diào)查研究論文集[C]. 北京: 海洋出版社, 1987, 267-284.

    Wang Guiyun. The content and distribution of chemical elements in the East China Sea, Kuroshio survey research papers election[C]. Beijing: Ocean Press, 1987, 267-284.

    [29] 韓舞鷹, 馬克美. 粵東沿岸上升流的研究[J].海洋學(xué)報, 1988, 10(1): 52-59.

    Han Wuying, Ma Kemei. Research on the coastal upwelling of Guangdong[J].Acta Oceanologica Sinica, 1988, 10(1): 52-59.

    [30] 楊紅, 李綱, 金成法. 長江口及鄰近海域磷酸鹽的再生和垂直通量[J]. 上海水產(chǎn)大學(xué)學(xué)報, 2005, 14(2): 162-167.

    Yang Hong, Li Gang, Jin Chengfa. The regeneration and vertical of phosphate in the Yangtze River Estuary and its adjacent area[J]. Journal of Shanghai Fisheries University, 2005, 14(2): 162-167.

    (本文編輯: 康亦兼)

    Influence of artificial reef construction on the transportation of nutrients in the off-shore area of Xiapu, Fujian

    XIAO Rong, YANG Hong
    (College of Marine Sciences, Shang Hai Ocean University, Shanghai 201306, China)

    Mar., 31, 2015

    artificial reef; effect of flow field; numerical simulation; upwelling; the vertical transportation flux of nutrients

    According to research performed in Xiapu, Fujian in 2012 and 2013, computational fluid dynamics was used to simulate the water flow pattern around three different artificial reefs in terms of unsteady flow effects, and the scale and intensity of the upwelling around the artificial reefs were simulated. In addition, the differences in flow characteristics were summarized, and the vertical transport fluxes of various nutrients were estimated. The results showed that the maximum height of the upwellings was approximately 2.27–3.73 times the reef height, and the maximum and average velocities of upwelling were approximately 0.40–0.72 and 0.10–0.17 times the incoming velocity. The upwelling velocity in a computational domain of three types of artificial reef assemblages was 10-2cm/s, and the magnitude of the relationship of nutrient vertical fluxes in different assemblages wassquare reefs>ship reefs>trapezoidal reefs. The average vertical fluxes of PO-P, NO-N, and SiO-Si were 438.3 mg/(m2·d), 2212.5 mg/(m2·d), and 7288.3 mg/(m2·d) in the square reef area.Setting artificial reefs could improve the hydrodynamic environment and form typical upwelling areas, significantly improving primary productivity.

    S953.1

    A

    1000-3096(2016)02-0094-08

    10.11759/hykx20150331002

    2015-03-31;

    2015-06-25

    農(nóng)業(yè)部2012年漁業(yè)資源保護和轉(zhuǎn)產(chǎn)轉(zhuǎn)業(yè)項目(D8006128012) [Foundation: The fishery resources protection and job transfer project in 2012 of Ministry of Agriculture (D8006128012)]

    肖榮(1988-), 男, 江西贛州人, 碩士研究生, 主要從事海洋環(huán)境學(xué)研究, 電話: 18817775220, E-mail: 1791198674@qq.com; 楊紅, 通信作者: E-mail: hyang@shou.edu.cn

    猜你喜歡
    上升流礁體魚礁
    夏季舟山上升流長期遙感觀測與分析
    人工魚礁區(qū)三種習(xí)見肉食性軟體動物營養(yǎng)生態(tài)位研究
    海洋通報(2022年5期)2022-11-30 12:04:06
    生態(tài)礁體布置形式對水流特征的影響
    水運工程(2022年3期)2022-03-25 03:15:06
    南麂列島海域兩種框型人工魚礁水動力性能試驗
    人工魚礁建設(shè)對海洋浮游植物群落變化的影響
    海洋通報(2020年5期)2021-01-14 09:27:04
    遼東半島頂端海域上升流長期變化特征及影響因素*
    海洋與湖沼(2020年1期)2020-02-06 06:01:38
    珠海萬山人工魚礁結(jié)構(gòu)設(shè)計
    人工魚礁在海洋牧場中的應(yīng)用
    珠江水運(2018年21期)2018-12-20 23:17:38
    海底新生代礁體下伏構(gòu)造的校正淺析
    海洋石油(2017年1期)2017-03-29 01:37:25
    瓊東上升流的年際變化及長期變化趨勢*
    海洋與湖沼(2016年1期)2016-01-15 03:50:49
    亚洲中文日韩欧美视频| 黄片小视频在线播放| 亚洲专区字幕在线| 一区二区三区国产精品乱码| 亚洲专区字幕在线| 97人妻精品一区二区三区麻豆 | 制服诱惑二区| 人人澡人人妻人| 国产精品免费视频内射| 亚洲欧美日韩无卡精品| 免费观看人在逋| 黄色片一级片一级黄色片| 久久人人爽av亚洲精品天堂| 国产一区二区三区视频了| 免费看a级黄色片| 国产午夜精品久久久久久| 午夜免费激情av| 欧美另类亚洲清纯唯美| 老熟妇乱子伦视频在线观看| 精品人妻在线不人妻| 日本欧美视频一区| 高清黄色对白视频在线免费看| 国内精品久久久久久久电影| 欧美乱码精品一区二区三区| 波多野结衣av一区二区av| 露出奶头的视频| 久9热在线精品视频| 91大片在线观看| 757午夜福利合集在线观看| 午夜影院日韩av| av天堂在线播放| 国产一区在线观看成人免费| 丝袜美腿诱惑在线| 国产成人欧美在线观看| 又黄又爽又免费观看的视频| 国产av一区二区精品久久| 最近最新中文字幕大全免费视频| 亚洲男人的天堂狠狠| 亚洲自偷自拍图片 自拍| 99国产极品粉嫩在线观看| 午夜福利免费观看在线| 女警被强在线播放| 在线观看舔阴道视频| www.自偷自拍.com| 在线观看免费视频网站a站| 性色av乱码一区二区三区2| 免费在线观看亚洲国产| 香蕉久久夜色| 国产精品香港三级国产av潘金莲| 69av精品久久久久久| 可以免费在线观看a视频的电影网站| 韩国av一区二区三区四区| 神马国产精品三级电影在线观看 | 精品人妻在线不人妻| 免费人成视频x8x8入口观看| 97超级碰碰碰精品色视频在线观看| 看片在线看免费视频| 中文字幕精品免费在线观看视频| 久久国产精品男人的天堂亚洲| 欧美日韩乱码在线| 老司机福利观看| 美国免费a级毛片| 一二三四在线观看免费中文在| 亚洲全国av大片| 两个人免费观看高清视频| 丰满的人妻完整版| 国产成人影院久久av| 中文字幕最新亚洲高清| 日韩中文字幕欧美一区二区| 在线天堂中文资源库| 好男人在线观看高清免费视频 | 12—13女人毛片做爰片一| 亚洲久久久国产精品| 成人手机av| 久久影院123| 黑人巨大精品欧美一区二区mp4| 午夜福利一区二区在线看| 中国美女看黄片| 制服人妻中文乱码| 丝袜在线中文字幕| 99国产综合亚洲精品| 欧美丝袜亚洲另类 | 中国美女看黄片| 天堂影院成人在线观看| 每晚都被弄得嗷嗷叫到高潮| av在线天堂中文字幕| 亚洲一区高清亚洲精品| 757午夜福利合集在线观看| 日韩精品免费视频一区二区三区| 满18在线观看网站| 国产精品一区二区在线不卡| 午夜福利高清视频| 俄罗斯特黄特色一大片| 欧洲精品卡2卡3卡4卡5卡区| 精品不卡国产一区二区三区| 国产欧美日韩精品亚洲av| 99香蕉大伊视频| 看片在线看免费视频| 亚洲伊人色综图| 日韩欧美免费精品| 级片在线观看| 亚洲精品中文字幕在线视频| 极品教师在线免费播放| 亚洲国产中文字幕在线视频| 亚洲中文字幕一区二区三区有码在线看 | av有码第一页| 久久草成人影院| 757午夜福利合集在线观看| 久热这里只有精品99| 啪啪无遮挡十八禁网站| 久久久久久久久免费视频了| 午夜免费成人在线视频| 午夜精品国产一区二区电影| 淫秽高清视频在线观看| 中文字幕色久视频| 国产精品99久久99久久久不卡| 宅男免费午夜| 亚洲精品国产精品久久久不卡| 99久久国产精品久久久| 国产精品国产高清国产av| 一区二区日韩欧美中文字幕| 后天国语完整版免费观看| 免费一级毛片在线播放高清视频 | 国产单亲对白刺激| 亚洲人成网站在线播放欧美日韩| 老汉色av国产亚洲站长工具| 成年女人毛片免费观看观看9| 日韩精品中文字幕看吧| 18禁国产床啪视频网站| 精品无人区乱码1区二区| 免费在线观看视频国产中文字幕亚洲| av视频免费观看在线观看| 成年人黄色毛片网站| 国产成人免费无遮挡视频| 午夜福利高清视频| 久久亚洲真实| 日韩大尺度精品在线看网址 | 午夜老司机福利片| 9热在线视频观看99| 欧美老熟妇乱子伦牲交| av片东京热男人的天堂| 一级黄色大片毛片| 啦啦啦 在线观看视频| 欧美色欧美亚洲另类二区 | 欧美乱妇无乱码| 91精品国产国语对白视频| 精品人妻1区二区| 亚洲aⅴ乱码一区二区在线播放 | 99久久综合精品五月天人人| 久久久久国产精品人妻aⅴ院| 一进一出好大好爽视频| 很黄的视频免费| 亚洲一区中文字幕在线| 天堂√8在线中文| 欧美色欧美亚洲另类二区 | 亚洲一卡2卡3卡4卡5卡精品中文| 搡老熟女国产l中国老女人| 久久亚洲精品不卡| tocl精华| 久久久久国产精品人妻aⅴ院| 国产精品免费一区二区三区在线| 久久久精品国产亚洲av高清涩受| 国产精品亚洲一级av第二区| 欧美日韩乱码在线| 国产精品一区二区在线不卡| 欧美av亚洲av综合av国产av| 1024视频免费在线观看| 欧美乱码精品一区二区三区| 日本精品一区二区三区蜜桃| 精品久久蜜臀av无| 丝袜在线中文字幕| a级毛片在线看网站| 精品国产乱子伦一区二区三区| 男女午夜视频在线观看| 精品欧美国产一区二区三| 国产亚洲精品综合一区在线观看 | 国语自产精品视频在线第100页| 国产精品一区二区三区四区久久 | www.熟女人妻精品国产| 18禁观看日本| 日日干狠狠操夜夜爽| 最新美女视频免费是黄的| 精品国产超薄肉色丝袜足j| 国产一区二区三区在线臀色熟女| 国产av一区在线观看免费| 男人舔女人的私密视频| 午夜福利在线观看吧| 欧美最黄视频在线播放免费| 天天添夜夜摸| 老司机午夜福利在线观看视频| 日韩精品中文字幕看吧| 亚洲伊人色综图| 国产精品久久视频播放| 精品高清国产在线一区| 宅男免费午夜| 亚洲av电影不卡..在线观看| 亚洲国产欧美日韩在线播放| 一级片免费观看大全| 亚洲色图综合在线观看| 国产精品久久久久久精品电影 | 日韩大尺度精品在线看网址 | xxx96com| 国产三级在线视频| 一级,二级,三级黄色视频| 亚洲 欧美一区二区三区| 美女免费视频网站| 三级毛片av免费| 老司机在亚洲福利影院| 黄色视频,在线免费观看| √禁漫天堂资源中文www| 老熟妇仑乱视频hdxx| 久久久久久久精品吃奶| 免费在线观看日本一区| 在线永久观看黄色视频| 无人区码免费观看不卡| 黑人巨大精品欧美一区二区mp4| 一个人免费在线观看的高清视频| 国产人伦9x9x在线观看| 久久狼人影院| 国产97色在线日韩免费| 午夜精品国产一区二区电影| 国产精华一区二区三区| 国产精品亚洲美女久久久| 国产成年人精品一区二区| 中文字幕色久视频| 人人妻,人人澡人人爽秒播| 又大又爽又粗| 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 久久人妻av系列| 天天躁狠狠躁夜夜躁狠狠躁| 国产极品粉嫩免费观看在线| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 热re99久久国产66热| 黄色女人牲交| 久久久久国产一级毛片高清牌| 国产一级毛片七仙女欲春2 | 久久国产精品人妻蜜桃| avwww免费| 国产国语露脸激情在线看| 欧美日韩精品网址| 搡老妇女老女人老熟妇| 亚洲午夜理论影院| 久久久水蜜桃国产精品网| 色av中文字幕| 无遮挡黄片免费观看| 99在线视频只有这里精品首页| 成人国产综合亚洲| 亚洲av电影在线进入| 日韩欧美国产在线观看| 亚洲国产精品sss在线观看| 欧美成狂野欧美在线观看| 一本综合久久免费| 99精品欧美一区二区三区四区| 大香蕉久久成人网| 国产亚洲精品av在线| 亚洲男人天堂网一区| 大码成人一级视频| 999久久久国产精品视频| 久久精品91无色码中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 日韩三级视频一区二区三区| 制服人妻中文乱码| 少妇被粗大的猛进出69影院| 色综合欧美亚洲国产小说| 97碰自拍视频| 久99久视频精品免费| 久久国产精品影院| 久久久久久久精品吃奶| 国产精品乱码一区二三区的特点 | 亚洲一区二区三区不卡视频| 啦啦啦 在线观看视频| 国产片内射在线| 丁香欧美五月| 国产97色在线日韩免费| 亚洲少妇的诱惑av| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一卡2卡三卡4卡5卡| 极品教师在线免费播放| 亚洲av成人一区二区三| 亚洲精品粉嫩美女一区| 国内精品久久久久精免费| 久久精品aⅴ一区二区三区四区| 国产亚洲欧美98| 国产成人欧美| 青草久久国产| 在线观看免费视频日本深夜| cao死你这个sao货| 两性午夜刺激爽爽歪歪视频在线观看 | 日本黄色视频三级网站网址| 日韩大码丰满熟妇| 久久国产精品影院| 午夜福利在线观看吧| 成人永久免费在线观看视频| 90打野战视频偷拍视频| 日本在线视频免费播放| 亚洲精品久久成人aⅴ小说| aaaaa片日本免费| 精品国产亚洲在线| 亚洲欧美一区二区三区黑人| 搞女人的毛片| 免费看美女性在线毛片视频| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区三区综合在线观看| 一级毛片女人18水好多| 亚洲激情在线av| 婷婷丁香在线五月| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人精品二区| 91字幕亚洲| 麻豆av在线久日| 不卡av一区二区三区| 国产精品爽爽va在线观看网站 | 黄色毛片三级朝国网站| 国产一区二区激情短视频| 亚洲情色 制服丝袜| 在线观看免费视频网站a站| 狂野欧美激情性xxxx| 午夜福利高清视频| 国产麻豆成人av免费视频| 免费在线观看影片大全网站| 久久九九热精品免费| 无限看片的www在线观看| 日韩三级视频一区二区三区| 91麻豆av在线| 久久狼人影院| 多毛熟女@视频| 久久国产精品男人的天堂亚洲| 国产97色在线日韩免费| 波多野结衣av一区二区av| 欧美色欧美亚洲另类二区 | 免费不卡黄色视频| 国产av在哪里看| 国产精品美女特级片免费视频播放器 | 啪啪无遮挡十八禁网站| 九色亚洲精品在线播放| a级毛片在线看网站| 激情在线观看视频在线高清| 亚洲精品中文字幕一二三四区| 国产熟女xx| 一个人免费在线观看的高清视频| 国产亚洲精品av在线| 淫秽高清视频在线观看| 欧美乱码精品一区二区三区| 国产亚洲精品久久久久5区| netflix在线观看网站| 美女高潮到喷水免费观看| 啦啦啦韩国在线观看视频| 国产精品二区激情视频| 亚洲 欧美一区二区三区| 天堂√8在线中文| 人妻丰满熟妇av一区二区三区| 国产精品,欧美在线| 国产片内射在线| 国产免费av片在线观看野外av| 法律面前人人平等表现在哪些方面| 美女免费视频网站| 亚洲国产中文字幕在线视频| 国产精品久久久人人做人人爽| 久久伊人香网站| 国产高清激情床上av| 老司机福利观看| 99re在线观看精品视频| svipshipincom国产片| 日本 av在线| 色av中文字幕| 少妇被粗大的猛进出69影院| 久久热在线av| 99在线人妻在线中文字幕| 欧美最黄视频在线播放免费| 久久久久久大精品| 好看av亚洲va欧美ⅴa在| 黄色视频,在线免费观看| 国产熟女xx| 免费看美女性在线毛片视频| 亚洲国产精品合色在线| 电影成人av| 国产精品免费视频内射| 青草久久国产| 一级a爱视频在线免费观看| 香蕉丝袜av| 午夜福利影视在线免费观看| 国产成人精品久久二区二区91| 又黄又粗又硬又大视频| 日日干狠狠操夜夜爽| 国产精品国产高清国产av| 一个人观看的视频www高清免费观看 | 国产精品久久久人人做人人爽| 国内久久婷婷六月综合欲色啪| 少妇被粗大的猛进出69影院| 黄片播放在线免费| 欧美激情极品国产一区二区三区| 一区二区三区国产精品乱码| 香蕉丝袜av| 亚洲va日本ⅴa欧美va伊人久久| 最近最新中文字幕大全免费视频| 黄色成人免费大全| 午夜精品在线福利| 好男人电影高清在线观看| 无人区码免费观看不卡| 午夜福利18| 日韩大码丰满熟妇| 一区二区三区高清视频在线| a在线观看视频网站| 国产一级毛片七仙女欲春2 | 色播在线永久视频| 国产亚洲精品久久久久5区| 女性被躁到高潮视频| 亚洲黑人精品在线| 免费无遮挡裸体视频| 久久天躁狠狠躁夜夜2o2o| 精品少妇一区二区三区视频日本电影| 99riav亚洲国产免费| 妹子高潮喷水视频| 人人妻人人澡人人看| 99国产极品粉嫩在线观看| 国产单亲对白刺激| 国产精品久久久久久精品电影 | 老司机靠b影院| 99精品久久久久人妻精品| 国产精品久久久久久亚洲av鲁大| 久久精品国产99精品国产亚洲性色 | 日日干狠狠操夜夜爽| 国产乱人伦免费视频| 十八禁人妻一区二区| 久久久久久免费高清国产稀缺| 久9热在线精品视频| 日韩精品青青久久久久久| 夜夜夜夜夜久久久久| 午夜福利欧美成人| 国产人伦9x9x在线观看| 亚洲熟妇熟女久久| 欧美日本视频| av网站免费在线观看视频| 一级黄色大片毛片| av福利片在线| e午夜精品久久久久久久| 精品久久久久久,| 校园春色视频在线观看| 深夜精品福利| 国产亚洲精品av在线| 狠狠狠狠99中文字幕| 曰老女人黄片| 亚洲五月婷婷丁香| 大型黄色视频在线免费观看| 亚洲一区高清亚洲精品| 免费一级毛片在线播放高清视频 | 在线永久观看黄色视频| 夜夜躁狠狠躁天天躁| 琪琪午夜伦伦电影理论片6080| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品无人区| 亚洲专区国产一区二区| av福利片在线| 一级a爱片免费观看的视频| www.精华液| 欧美激情 高清一区二区三区| 淫妇啪啪啪对白视频| 露出奶头的视频| 女性被躁到高潮视频| 十八禁人妻一区二区| 日本免费一区二区三区高清不卡 | 宅男免费午夜| 看黄色毛片网站| 两个人视频免费观看高清| 国产亚洲精品久久久久久毛片| 久99久视频精品免费| 波多野结衣av一区二区av| 亚洲第一电影网av| 国产精品日韩av在线免费观看 | 成人手机av| 乱人伦中国视频| 老司机福利观看| 99热只有精品国产| 精品国产国语对白av| 免费人成视频x8x8入口观看| 国产成人欧美| 少妇的丰满在线观看| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩高清在线视频| 久久中文字幕人妻熟女| 变态另类成人亚洲欧美熟女 | 国产真人三级小视频在线观看| 最新美女视频免费是黄的| 国产精品二区激情视频| 亚洲一区二区三区不卡视频| 日韩高清综合在线| 久久久久亚洲av毛片大全| 亚洲男人的天堂狠狠| 一二三四社区在线视频社区8| 啦啦啦韩国在线观看视频| 国产亚洲精品久久久久5区| 在线观看免费视频日本深夜| 久久这里只有精品19| 亚洲成人精品中文字幕电影| 最新美女视频免费是黄的| 90打野战视频偷拍视频| 亚洲欧美精品综合一区二区三区| 国产精品美女特级片免费视频播放器 | 操出白浆在线播放| 久久国产精品男人的天堂亚洲| 精品第一国产精品| 久久中文字幕人妻熟女| 国产亚洲av高清不卡| 久久这里只有精品19| 国产高清激情床上av| 国产激情欧美一区二区| 亚洲av第一区精品v没综合| 亚洲成人免费电影在线观看| 操美女的视频在线观看| 精品国产美女av久久久久小说| 中文字幕高清在线视频| 色播亚洲综合网| 老熟妇仑乱视频hdxx| 色av中文字幕| 国产精品九九99| 一边摸一边抽搐一进一出视频| 99在线视频只有这里精品首页| 亚洲国产精品999在线| 午夜久久久久精精品| 久久久国产精品麻豆| 777久久人妻少妇嫩草av网站| 亚洲一卡2卡3卡4卡5卡精品中文| 波多野结衣av一区二区av| 黑人欧美特级aaaaaa片| 免费女性裸体啪啪无遮挡网站| 日韩中文字幕欧美一区二区| 在线观看免费午夜福利视频| 嫩草影视91久久| 国产区一区二久久| 亚洲,欧美精品.| 精品卡一卡二卡四卡免费| 久久中文字幕一级| 日韩精品中文字幕看吧| 欧美成狂野欧美在线观看| 欧美av亚洲av综合av国产av| 我的亚洲天堂| 国产成人精品久久二区二区91| 每晚都被弄得嗷嗷叫到高潮| 校园春色视频在线观看| 琪琪午夜伦伦电影理论片6080| 777久久人妻少妇嫩草av网站| 18禁黄网站禁片午夜丰满| 在线观看一区二区三区| 成人精品一区二区免费| 在线观看午夜福利视频| 亚洲美女黄片视频| 欧美一级a爱片免费观看看 | 国产精品av久久久久免费| 麻豆成人av在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲视频免费观看视频| 欧美日本视频| 青草久久国产| 啦啦啦免费观看视频1| 可以在线观看毛片的网站| 电影成人av| 亚洲精品在线观看二区| 九色国产91popny在线| 日本免费a在线| 麻豆久久精品国产亚洲av| 精品欧美一区二区三区在线| 国产欧美日韩精品亚洲av| 亚洲午夜精品一区,二区,三区| 女警被强在线播放| 久久天堂一区二区三区四区| 国产精品亚洲美女久久久| 波多野结衣一区麻豆| 日韩 欧美 亚洲 中文字幕| 看片在线看免费视频| 久久久久国内视频| 久久国产乱子伦精品免费另类| 一级毛片女人18水好多| 久久亚洲真实| 国产男靠女视频免费网站| 成人18禁高潮啪啪吃奶动态图| 国产成人精品在线电影| 激情视频va一区二区三区| 少妇的丰满在线观看| 国产午夜精品久久久久久| 久久国产精品影院| 欧美中文日本在线观看视频| 国产一区二区三区综合在线观看| 午夜激情av网站| 老司机午夜十八禁免费视频| 国产精品免费一区二区三区在线| 日韩国内少妇激情av| 给我免费播放毛片高清在线观看| 午夜福利欧美成人| 99久久99久久久精品蜜桃| 成年人黄色毛片网站| 在线观看免费视频网站a站| 国产亚洲精品av在线| 真人一进一出gif抽搐免费| avwww免费| 韩国av一区二区三区四区| 男人操女人黄网站| 久久婷婷人人爽人人干人人爱 | 十八禁人妻一区二区| 国产又爽黄色视频| 色播在线永久视频| 制服诱惑二区| av视频免费观看在线观看| 中文字幕精品免费在线观看视频| 日韩有码中文字幕| 国产私拍福利视频在线观看| 欧美精品啪啪一区二区三区| 国产麻豆69| 亚洲av成人av| 精品少妇一区二区三区视频日本电影| 久久草成人影院| 91成年电影在线观看| 亚洲精品中文字幕在线视频| 色尼玛亚洲综合影院|