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

    雙殼油船破艙原油泄漏縮尺模型相似準(zhǔn)則

    2017-11-03 02:46:55盧金樹(shù)吳豪霄楊振波
    中國(guó)航海 2017年3期
    關(guān)鍵詞:油艙縮尺模型試驗(yàn)

    盧金樹(shù), 吳豪霄, 楊振波

    (浙江海洋大學(xué) 海運(yùn)與港航建筑工程學(xué)院, 浙江 舟山 316022)

    2017-04-15

    國(guó)家自然科學(xué)基金委員會(huì)面上項(xiàng)目(51079129);浙江省自然基金(LY18E090008)

    盧金樹(shù)(1974—),男,浙江象山人,教授,從事船舶安全與污染控制研究。E-mail:ljs_ljs@zjou.edu.cn

    1000-4653(2017)03-0103-06

    雙殼油船破艙原油泄漏縮尺模型相似準(zhǔn)則

    盧金樹(shù), 吳豪霄, 楊振波

    (浙江海洋大學(xué) 海運(yùn)與港航建筑工程學(xué)院, 浙江 舟山 316022)

    為使雙殼油船原油泄漏縮尺模型試驗(yàn)?zāi)苷鎸?shí)地反映事故原型,需確定各相似準(zhǔn)則的重要程度。針對(duì)油艙原油泄漏過(guò)程,基于相似理論,采用量綱分析的方法推導(dǎo)較完整的相似準(zhǔn)則數(shù)群,根據(jù)不同相似準(zhǔn)則的組合建立3種縮尺模型。進(jìn)行原型和不同相似準(zhǔn)則縮尺模型的數(shù)值試驗(yàn),并根據(jù)相似關(guān)系將縮尺模型試驗(yàn)的泄漏持續(xù)時(shí)間、泄漏總量等宏觀泄漏特征參數(shù)反推回原型,與原型數(shù)值試驗(yàn)結(jié)果相對(duì)比。采用黏性流體動(dòng)力學(xué)理論分析各種作用力對(duì)泄漏過(guò)程的影響。結(jié)果表明:弗汝德相似準(zhǔn)則和雷諾相似準(zhǔn)則是原油泄漏縮尺模型試驗(yàn)必須依據(jù)的準(zhǔn)則,韋伯相似準(zhǔn)則可不予考慮。在試驗(yàn)油品和原油密度相同的情況下,縮尺模型試驗(yàn)結(jié)果與原型最為接近。針對(duì)不同泄漏階段,應(yīng)采用不同的縮尺模型進(jìn)行試驗(yàn)。

    船舶工程;縮尺模型;相似準(zhǔn)則;原油泄漏;雙殼油船;黏性流體動(dòng)力學(xué)

    在原油泄漏研究中,原油泄漏規(guī)律的精準(zhǔn)預(yù)測(cè)是研究者關(guān)注的兩大工程問(wèn)題之一。[1-2]針對(duì)該問(wèn)題,通常采用理論分析、數(shù)值模擬和模型試驗(yàn)等方法進(jìn)行研究。由于油品泄漏涉及到的流動(dòng)行為復(fù)雜多樣,采用理論分析方法時(shí)需對(duì)泄漏預(yù)測(cè)模型進(jìn)行大量簡(jiǎn)化,導(dǎo)致模型精度不高。大尺度模型試驗(yàn)?zāi)艿玫捷^為真實(shí)的原油泄漏規(guī)律,但大尺度模型艙的制作成本較高,且原油泄漏會(huì)造成嚴(yán)重的環(huán)境污染。[3]因此,通常采用小尺度縮尺模型試驗(yàn)的方法研究原油泄漏規(guī)律。[4-7]縮尺模型試驗(yàn)的精度主要受以下3個(gè)因素影響:

    1) 確定原油泄漏過(guò)程的主導(dǎo)力,選取合理的相似準(zhǔn)則。

    2) 選取表征泄漏特征的特征長(zhǎng)度。

    3) 選取適當(dāng)大小的縮尺比。

    關(guān)鍵作用力的相似準(zhǔn)則缺失、表征泄漏特征的特征長(zhǎng)度選擇錯(cuò)誤及縮尺比過(guò)小都有可能引發(fā)縮尺效應(yīng)。根據(jù)建立縮尺模型所依據(jù)的相似準(zhǔn)則分類,采用的縮尺模型主要有弗汝德相似準(zhǔn)則模型(以下簡(jiǎn)稱Fr模型)和弗汝德-雷諾相似準(zhǔn)則模型(以下簡(jiǎn)稱Fr-Re模型)2種。

    Fr模型是早期的原油泄漏模型試驗(yàn)中采用的縮尺模型,主要用來(lái)探究雙殼油船和中層甲板油船的阻漏效果。[8-9]由于Fr模型只考慮重力的相似,忽略黏性力的相似,不可避免地會(huì)產(chǎn)生縮尺效應(yīng)。[8]此外,由于原油泄漏過(guò)程中存在流動(dòng)變態(tài)行為,在舷側(cè)破艙原油泄漏的密度差流階段[10-11],油品泄漏的動(dòng)力并非重力。在這種情形下,F(xiàn)r模型就不能真實(shí)地反映原型的泄漏過(guò)程。

    考慮到黏性力未相似會(huì)引發(fā)縮尺效應(yīng),文獻(xiàn)[12]針對(duì)單殼油船原油泄漏情形,提出同時(shí)考慮重力和黏性力相似的Fr-Re模型。該縮尺模型是在滿足Fr相似的前提下,選取特定黏度的油品進(jìn)行試驗(yàn),使模型的雷諾數(shù)與原型相等。當(dāng)縮尺比較大時(shí),難以找到黏度符合要求的油品[13],使得該模型具有一定的局限性,而且忽略了表面張力和浮力,采用該模型進(jìn)行試驗(yàn)時(shí)可能會(huì)引發(fā)縮尺效應(yīng)。

    基于上述情況,首先對(duì)原油泄漏過(guò)程進(jìn)行受力分析,確定影響原油泄漏的相關(guān)作用力?;谙嗨评碚摚瑢?duì)油艙原油泄漏過(guò)程進(jìn)行量綱分析,推導(dǎo)出較為完整的相似準(zhǔn)則數(shù)群,并據(jù)此建立3種不同的縮尺模型。在相同縮尺比下,通過(guò)對(duì)比分析3種縮尺模型的泄漏持續(xù)時(shí)間、泄漏速度和泄漏入海量等宏觀泄漏特征參數(shù),確定影響原油泄漏的主要作用力,最終得到建立縮尺模型必須依據(jù)的相似準(zhǔn)則。

    1 原油泄漏受力分析

    雙殼油船水下破艙原油泄漏是多種流動(dòng)行為復(fù)合的過(guò)程。在不同損傷情形下,原油泄漏的主導(dǎo)力不同,且作用強(qiáng)度會(huì)隨著泄漏過(guò)程的進(jìn)行而發(fā)生改變。這里以靜水環(huán)境中的雙殼油船底部破艙泄漏情形為例,分析原油泄漏過(guò)程中的相關(guān)作用力。

    圖1為雙殼油船原油泄漏示意,由于雙層底空間的存在,從油艙中泄漏出的原油不會(huì)直接泄漏到外界,而是先泄漏到雙層底中,隨后穿透雙層底過(guò)程中的水層泄漏到外界。[7]通常認(rèn)為原油從油艙泄漏到雙層底中的流動(dòng)屬于重力主導(dǎo)下的孔口出流,其流動(dòng)行為還受黏性力和油-氣表面張力的影響。[8]然而,在原油穿透雙層底中的水層泄漏到外界的過(guò)程中,其還會(huì)受到浮力的作用。在這些作用力中,重力是驅(qū)動(dòng)原油泄漏的力(動(dòng)力),而黏性力、表面張力和浮力均是阻礙原油泄漏的力(阻力)。當(dāng)某個(gè)作用力未按相同縮尺比縮小時(shí),就會(huì)導(dǎo)致動(dòng)力和阻力的相對(duì)大小發(fā)生變化,從而引發(fā)縮尺效應(yīng)。

    圖1 雙殼油船原油泄漏示意

    2 縮尺模型的建立

    目前已有研究采用的縮尺模型均未考慮表面張力和浮力。表面張力和浮力分別在油艙及雙層底中起作用。這里以油艙為研究對(duì)象,針對(duì)原油從油艙泄漏到雙層底的過(guò)程,采用量綱分析的方法推導(dǎo)出考慮表面張力的較為完整的相似準(zhǔn)則數(shù)群,并據(jù)此建立縮尺模型。

    2.1油艙油品泄漏量綱分析

    內(nèi)殼破孔處的原油泄漏平均速度v是衡量原油泄漏特征的重要指標(biāo)。相關(guān)參數(shù)有油艙油位H,破孔孔徑d,重力加速度g,油品密度ρ,油品動(dòng)力黏度系數(shù)μ和表面張力系數(shù)σ。描述油艙油品泄漏速度的一般函數(shù)式為

    f(v,H,ρ,g,d,μ,σ)=0

    (1)

    式(1)中:物理量的個(gè)數(shù)為n=7,采用[M],[L],[T]等3個(gè)基本量綱,m=3。根據(jù)π定理,這n個(gè)物理量的關(guān)系式可轉(zhuǎn)換成n-m=4個(gè)無(wú)量綱量的函數(shù)關(guān)系式,即

    f(π1,π2,π3,π4)=0

    (2)

    選擇H,ρ,g等3個(gè)物理量作為基本物理量,其余4個(gè)參數(shù)可表示為

    (3)

    根據(jù)量綱一致性原理得到4個(gè)無(wú)量綱數(shù),即

    (4)

    式(3)和式(4)中:π1為弗汝德數(shù);π2為幾何無(wú)量綱數(shù);π3和π4分別可轉(zhuǎn)換為雷諾數(shù)及韋伯?dāng)?shù)。弗汝德數(shù)、雷諾數(shù)和韋伯?dāng)?shù)分別是重力、黏性力及表面張力的相似準(zhǔn)則。根據(jù)滿足的相似準(zhǔn)則個(gè)數(shù)的不同,可建立Fr模型、Fr-Re模型和弗汝德-雷諾-韋伯模型(以下簡(jiǎn)稱Fr-Re-We模型)等3種縮尺模型。

    2.2原油泄漏模型參數(shù)相似關(guān)系的導(dǎo)出

    用m和p分別表示模型參數(shù)及原型參數(shù)。取幾何比尺為λ=lp/lm,根據(jù)推導(dǎo)的相似準(zhǔn)則導(dǎo)出模擬油艙原油泄漏過(guò)程的縮尺模型各參數(shù)與原型的相似關(guān)系。

    1) 速度相似關(guān)系為vp=vmλ1/2。

    2) 破孔孔徑相似關(guān)系為dp=dmλ。

    3) 油品動(dòng)力黏度相似關(guān)系為μp=μmλρλ3/2。

    4) 油品表面張力系數(shù)相似關(guān)系為σp=σmλρλ2。

    5) 當(dāng)模型的油品密度與原型相等時(shí),動(dòng)力黏度與表面張力系數(shù)相似關(guān)系分別簡(jiǎn)化為μp=μmλ3/2和σp=σmλ2。

    根據(jù)幾何相似關(guān)系和速度相似關(guān)系,可導(dǎo)出時(shí)間相似關(guān)系tp=tmλ1/2和油品體積相似關(guān)系Vp=Vmλ3。基于不同相似準(zhǔn)則建立的不同縮尺模型主要比尺見(jiàn)表1。

    表1 不同縮尺模型主要比尺

    3 不同相似準(zhǔn)則縮尺模型研究

    選取的相似準(zhǔn)則越合理,根據(jù)相似關(guān)系將縮尺模型試驗(yàn)結(jié)果反推回原型尺度得到的模型-原型結(jié)果與實(shí)際原型結(jié)果越接近。通過(guò)對(duì)比不同縮尺模型的油品泄漏特征參數(shù)與原型的偏差,即可判斷出各相似準(zhǔn)則的重要程度。參考相關(guān)研究[7],以某超大型油船的L型油艙邊艙為模型艙原型,以Lone Rock原油為油品原型,根據(jù)表1設(shè)計(jì)縮尺模型試驗(yàn)。采用數(shù)值試驗(yàn)的方法模擬原型和各縮尺模型的油品泄漏過(guò)程。通過(guò)對(duì)比模型自由泄漏階段持續(xù)的時(shí)間及油品泄漏入海量與原型的差異,探究油品泄漏過(guò)程中黏性力、表面張力和浮力的相對(duì)重要程度,進(jìn)而確定建立縮尺模型時(shí)應(yīng)采用的相似準(zhǔn)則。

    3.1縮尺模型試驗(yàn)設(shè)計(jì)依據(jù)

    在縮尺模型試驗(yàn)中,模型縮尺比確定之后,可通過(guò)改變特征流速v和選取具有特定黏度μ的液體作為流體介質(zhì)2種方法控制模型雷諾數(shù)Re=ρvd/μ。然而,在原油泄漏模型試驗(yàn)中,無(wú)法人為改變流速,只能通過(guò)選取具有特定黏度的油品控制模型的雷諾數(shù)。[13-14]同樣,模型的韋伯?dāng)?shù)也通過(guò)改變油品的表面張力系數(shù)來(lái)控制。

    3.2縮尺模型油品物性及模型幾何參數(shù)的確定

    小尺度縮尺模型試驗(yàn)的縮尺比通常為1:50~1:30,因此選取1:40作為縮尺模型的縮尺比。根據(jù)表1中的油品物性相似關(guān)系和模型幾何特征相似關(guān)系,確定縮尺模型試驗(yàn)對(duì)照組的油品物性和模型幾何參數(shù)(見(jiàn)表2)。此外,為探究浮力對(duì)油品泄漏的影響,設(shè)置試驗(yàn)油品密度和Fr-Re-We模型不同的對(duì)照組Fr-Re-We(D)。

    為保證縮尺模型與原型的幾何相似性,在原型和縮尺模型的幾何模型中均忽略船殼厚度。簡(jiǎn)化的模型艙示意見(jiàn)圖2。

    a)正視圖b)俯視圖

    圖2 簡(jiǎn)化的模型艙示意

    3.3研究方法及驗(yàn)證

    該研究采用文獻(xiàn)[15]給出的數(shù)值方法,經(jīng)物模試驗(yàn)驗(yàn)證,該方法能較為準(zhǔn)確地模擬油品泄漏過(guò)程。

    3.4不同相似準(zhǔn)則縮尺模型的試驗(yàn)結(jié)果

    相似準(zhǔn)則數(shù)反映各種作用力和慣性力的相對(duì)大小。根據(jù)Fr模型、Fr-Re模型、Fr-Re-We模型和Fr-Re-We(D)模型試驗(yàn)結(jié)果,分析黏性力、表面張力和浮力對(duì)原油泄漏的影響,從而確定縮尺模型試驗(yàn)應(yīng)采用的相似準(zhǔn)則。

    3.4.1不同相似準(zhǔn)則下的泄漏持續(xù)時(shí)間

    理論上,當(dāng)所有作用力都按相同比例縮小時(shí),根據(jù)相似關(guān)系將縮尺模型結(jié)果反推回原型得到的結(jié)果應(yīng)與原型相同。當(dāng)某個(gè)相似準(zhǔn)則未考慮(即某個(gè)作用力未按縮尺比縮小)時(shí),縮尺模型轉(zhuǎn)化到原型的結(jié)果會(huì)與原型產(chǎn)生偏差,該偏差就是該作用力未縮尺產(chǎn)生的縮尺效應(yīng)。不同縮尺模型油品自由泄漏階段持續(xù)時(shí)間見(jiàn)圖3。從圖3中可看出,Fr-Re模型和Fr-Re-We模型的泄漏持續(xù)時(shí)間與原型最為接近,F(xiàn)r-Re-We(D)模型次之,F(xiàn)r模型偏差最大。這說(shuō)明黏度偏大、密度偏小會(huì)導(dǎo)致原油自由泄漏階段持續(xù)時(shí)間偏短,黏性力和浮力是引發(fā)縮尺效應(yīng)的主要原因。

    表2 原油泄漏縮尺模型試驗(yàn)對(duì)照組

    圖3 不同縮尺模型油品自由泄漏階段持續(xù)時(shí)間

    Fr-Re-We模型是在Fr-Re模型的基礎(chǔ)上考慮韋伯相似準(zhǔn)則得到的縮尺模型。通過(guò)對(duì)比Fr-Re模型和Fr-Re-We模型可探究表面張力對(duì)原油泄漏的影響。圖3中,2個(gè)縮尺模型的油品泄漏速度和油品從外殼破孔泄漏入海的平均速度均基本相同,說(shuō)明在油品泄漏過(guò)程中表面張力的作用可忽略不計(jì)。

    從圖3~圖7中可看出,F(xiàn)r-Re模型和Fr-Re-We模型已能較為真實(shí)地反映原型,但仍存在偏差,這可能是未采用相應(yīng)的相似準(zhǔn)則對(duì)油水界面張力進(jìn)行縮尺導(dǎo)致的。

    Fr-Re模型是在Fr模型的基礎(chǔ)上考慮雷諾相似準(zhǔn)則得到的縮尺模型。通過(guò)對(duì)比Fr模型和Fr-Re模型可探究黏性力對(duì)原油泄漏的影響。圖4為不同縮尺模型油艙油品泄漏平均速度。從圖4中可看出,F(xiàn)r模型中油品從油艙破孔泄漏進(jìn)入雙層底的平均速度略小于Fr-Re模型。這是因?yàn)镕r模型僅保證重力的相似,而黏性力未按相同比例縮小。當(dāng)采用原油作為試驗(yàn)油品時(shí),由于其黏度遠(yuǎn)大于黏性力相似應(yīng)采用的油品黏度[10],導(dǎo)致油品在油艙中流動(dòng)的阻力偏大,油艙油品泄漏平均速度偏小。

    圖4 不同縮尺模型油艙油品泄漏平均速度

    在Fr模型中,油品進(jìn)入雙層底之后,由于油艙中的油品泄漏速度偏小,油射流所具備的動(dòng)能不足以克服自身的浮力和水層的阻力,導(dǎo)致油品泄漏入海率遠(yuǎn)低于原型。此外,由于油射流不能占據(jù)外殼破孔,使得外界海水大量進(jìn)入雙層底。在這2方面的綜合作用下,F(xiàn)r模型的雙層底液位上升速度遠(yuǎn)大于原型和其他縮尺模型,最終導(dǎo)致該模型自由泄漏階段持續(xù)時(shí)間遠(yuǎn)小于原型。因此,黏性力是原油泄漏過(guò)程中的重要作用力,當(dāng)黏性力未按相同比例縮小時(shí),會(huì)導(dǎo)致縮尺模型自由泄漏階段持續(xù)時(shí)間存在縮尺效應(yīng)。這也說(shuō)明在進(jìn)行縮尺模型試驗(yàn)時(shí),不能選取原油作為試驗(yàn)油品,而應(yīng)根據(jù)雷諾相似準(zhǔn)則選取特定黏度的油品。

    Fr-Re-We(D)模型是在Fr-Re-We模型的基礎(chǔ)上改變?cè)囼?yàn)油品密度得到的縮尺模型。通過(guò)對(duì)比Fr-Re-We(D)模型和Fr-Re-We模型可探究浮力對(duì)原油泄漏的影響。圖5為不同縮尺模型油品泄漏入海平均速度,2種模型的油艙油品泄漏速度基本相同,F(xiàn)r-Re-We(D)模型的油品泄漏入海平均速度小于Fr-Re-We模型,與原型的偏差更大。

    上述試驗(yàn)證明油品在穿透雙層底中的油水混合物時(shí)會(huì)受到浮力的影響。當(dāng)油品密度小于原油時(shí),油品與海水的密度差變大,油品在油水混合物中流

    圖5 不同縮尺模型油品入海泄漏平均速度

    動(dòng)時(shí)受到的凈浮力變大,油品泄漏入海的速度變小。泄漏入海速度變小意味著更多油品被保留在雙層底中,雙層底液位上升速度加快,自由泄漏階段持續(xù)時(shí)間變短(如圖3所示)。因此,浮力是油品泄漏過(guò)程中的另一個(gè)重要作用力,當(dāng)其未按相同比例縮小時(shí),會(huì)導(dǎo)致縮尺模型自由泄漏階段持續(xù)時(shí)間存在縮尺效應(yīng)。

    由于油品密度發(fā)生改變,F(xiàn)r-Re-We(D)模型的雷諾數(shù)并不與Fr-Re-We模型相同,但數(shù)值非常接近。然而,油品泄漏入海的速度與Fr-Re-We模型存在很大偏差,這說(shuō)明僅保證雷諾數(shù)相同不能準(zhǔn)確模擬原油在雙層底中的流動(dòng)行為。油品的密度和黏度是2個(gè)獨(dú)立的影響因子,通過(guò)控制油品的運(yùn)動(dòng)黏度來(lái)控制雷諾數(shù)相同是不可行的。

    3.4.2不同相似準(zhǔn)則下的油品泄漏入海量

    油品泄漏入海量是另一個(gè)重要的油品泄漏特征參數(shù)。不同縮尺模型自由泄漏階段油品泄漏入海量見(jiàn)圖6。從圖6中可看出,F(xiàn)r-Re模型和Fr-Re-We模型的泄漏入海量與原型最為接近,F(xiàn)r-Re-We(D)模型和Fr模型的偏差均較大。這說(shuō)明黏性力和浮力也會(huì)導(dǎo)致油品泄漏入海量存在縮尺效應(yīng)。

    圖6 不同縮尺模型自由泄漏階段油品泄漏入海量

    油艙油品泄漏量取決于自由泄漏階段持續(xù)時(shí)間和油品泄漏入海的平均速度。由于Fr模型未對(duì)黏性力進(jìn)行縮尺,采用原油作為試驗(yàn)油品會(huì)導(dǎo)致油品泄漏平均速度偏小,自由泄漏階段持續(xù)時(shí)間偏短。這2個(gè)參數(shù)偏小均會(huì)導(dǎo)致油艙油品泄漏量偏少(見(jiàn)圖7)。

    圖7 不同縮尺模型油艙油品泄漏量

    由于黏性力過(guò)大,使得油射流不能穿透雙層底中的水層泄漏到外界,直接改變雙層底中油水相互作用形式,導(dǎo)致模型油品泄漏入海量與原型存在巨大偏差(如圖6所示)。這進(jìn)一步證明要采用縮尺模型模擬原油泄漏過(guò)程,必須根據(jù)雷諾相似準(zhǔn)則選取低黏度的油品作為試驗(yàn)油品。

    在自由泄漏階段,油艙中的油品與雙層底中的油水混合物并未發(fā)生直接作用,浮力僅對(duì)雙層底中的油品產(chǎn)生作用,因此Fr-Re-We(D)模型的油艙油品泄漏速度與Fr-Re-We模型基本相同。由于Fr-Re-We(D)模型采用的油品與海水的密度差更大,油品在雙層底中流動(dòng)時(shí)所受凈浮力更大,油射流速度衰減更嚴(yán)重。因此,Fr-Re-We(D)模型的油品泄漏率偏小,雙層底中油品積累的速度偏大,自由泄漏階段持續(xù)時(shí)間偏短。雖然Fr-Re-We(D)模型油艙油品泄漏速度與Fr-Re-We模型相同,但Fr-Re-We(D)模型自由泄漏階段持續(xù)時(shí)間更短,導(dǎo)致油艙泄漏量偏少。在油艙泄漏量和油品泄漏入海率均減小的情況下,F(xiàn)r-Re-We(D)模型自由泄漏階段的油品泄漏入海量偏少(如圖6所示)。

    綜上所述,當(dāng)模型中的油品所受的黏性力和凈浮力偏大時(shí),會(huì)導(dǎo)致自由泄漏階段持續(xù)時(shí)間縮短,該階段的油品泄漏入海量減少。然而,在進(jìn)行縮尺物模試驗(yàn)時(shí),較難使浮力和黏性力的大小同時(shí)滿足相似準(zhǔn)則的要求,這勢(shì)必會(huì)導(dǎo)致自由泄漏階段結(jié)束時(shí)的油品泄漏特征參數(shù)存在縮尺效應(yīng)。若縮尺模型進(jìn)入原油受阻泄漏階段時(shí)的油艙油位和油艙油品泄漏速度(轉(zhuǎn)化到原型尺度)與原型存在偏差,則受阻泄漏階段的縮尺效應(yīng)可能會(huì)被進(jìn)一步放大。

    4 結(jié)束語(yǔ)

    針對(duì)原油自由泄漏階段的油艙原油泄漏過(guò)程,基于相似理論,采用量綱分析的方法推導(dǎo)出完整的相似準(zhǔn)則數(shù)群,并根據(jù)不同相似準(zhǔn)則建立3種縮尺模型。通過(guò)數(shù)值試驗(yàn)的方法模擬相同縮尺比下不同縮尺模型的油品泄漏過(guò)程,根據(jù)相似關(guān)系將縮尺模型數(shù)值試驗(yàn)結(jié)果轉(zhuǎn)換為原型尺度,將其與原型數(shù)值試驗(yàn)結(jié)果相對(duì)比,驗(yàn)證通過(guò)選取低黏度油品控制模型雷諾數(shù)的科學(xué)性。研究結(jié)果表明:

    1) 黏性力和浮力是影響原油泄漏過(guò)程的重要作用力,這2個(gè)作用力未按相同縮尺比縮小是引發(fā)縮尺效應(yīng)的主要原因。油-氣表面張力在原油泄漏過(guò)程中可忽略不計(jì)。

    2) 在進(jìn)行縮尺物模試驗(yàn)時(shí),應(yīng)根據(jù)傅汝德相似準(zhǔn)則建立縮尺模型,根據(jù)雷諾相似準(zhǔn)則選取適當(dāng)黏度的油品或盡量選取黏度小的油品,油品密度應(yīng)與原油相等或接近,以降低縮尺效應(yīng)的影響。

    [1] HYSING T, TORSET O. Reduction of Oil Outflows at Collisions and Groundings Through Improved Vessel Design Arrangement[J]. Marine Pollution Bulletin, 1994, 29(6): 368-374.

    [2] WIEL G V D, DORP J R V. An Oil Outflow Model for Tanker Collisions and Groundings[J]. Annals of Operations Research, 2011, 187(1): 279-304.

    [3] 蘭國(guó)新, 李穎, 劉丙新, 等. 基于多源遙感的海上溢油檢測(cè)研究[J]. 中國(guó)航海, 2013, 36(1): 115-120.

    [4] TAVAKOLI M T, AMDAHL J, LEIRA B J. Investigation of Interaction Between Oil Spills and Hydrostatic Changes[C]. ASME International Conference on Ocean, Offshore and Arctic Engineering, 2009: 803-811.

    [5] TAVAKOLI M T, AMDAHL J, LEIRA B J. Analytical and Numerical Modelling of Oil Spill from a Side Tank with Collision Damage[J]. Ships and Offshore Structures, 2012, 7(1): 73-86.

    [6] 趙衛(wèi)斌, 盧金樹(shù), 朱哲野, 等. 油船水下破艙原油泄漏模型試驗(yàn)設(shè)計(jì)[J]. 浙江海洋學(xué)院學(xué)報(bào)(自然科學(xué)版), 2013, 32(3): 255-259.

    [7] LU J S, LIU F C, ZHU Z Y. Effects of Initial Water Layer Thickness on Oil Leakage from Damaged DHTs[C]. The Twenty-Fourth International Ocean and Polar Engineering Conference, 2014: 618-625.

    [8] KARAFIATH G. Accidental Oil Spill Due to Grounding: Summary of Model Test Results. Summary Report, Jan-Jun 92[R]. Bethesda, MD Naval Surface Warfare Center, 1992: 1-61.

    [9] YAMAGUCHI K, YAMANOUCHI H. Oil Spills from the Double Hull Model Tanks[J]. Report of Ship Research Institute, 1992, 29: 1-38.

    [10] 盧金樹(shù). 油船水下破艙油品泄漏動(dòng)力學(xué)機(jī)理及特性研究[D]. 大連:大連海事大學(xué), 2015:49-51.

    [11] 盧金樹(shù), 劉楓琛, 吳文鋒, 等. 破孔高度對(duì)油船水下破艙油品泄漏的影響[J]. 大連海事大學(xué)學(xué)報(bào), 2015,41(2): 46-50.

    [12] SIMECEK-BEATTY D, LEHR W J, LANKFORD J. Leaking Tank Experiments for Heavy Oils[J]. Beatty, 2012(1): 127-131.

    [13] TAVAKOLI M T, AMDAHL J, LEIRA B J. Experimental Investigation of Oil Leakage from Damaged Ships Due to Collision and Grounding[J]. Ocean Engineering, 2011, 38(17): 1894-1907.

    [14] SIMECEK-BEATTY D, LEHR W J, LANKFORD J F. Leaking Tank Experiments With Orimulsion and Canola Oil[J]. National Oceanic and Atmospheric Administration, 2003(12): 30.

    [15] 盧金樹(shù), 劉楓琛, 朱哲野. 雙殼油船底部破艙水下原油泄漏過(guò)程三維數(shù)值模擬[J]. 江蘇科技大學(xué)學(xué)報(bào) (自然科學(xué)版), 2014, 28(5): 423-428.

    SimilarityCriteriaofPhysicalModelingforCrudeOilLeakagefromDoubleHullOilTanker

    LUJinshu,WUHaoxiao,YANGZhenbo

    (School of Maritime and Civil Engineering, Zhejiang Ocean University, Zhoushan 316022, China)

    Laboratory investigation of oil leakage from double hull oil tank requires the selection of an adequate similitude. In order to ensure the rationality of the scale model tests, the dimension analysis method is used to deduce the similarity criteria that are required in the scale model test. Based on these similarity criteria, three different kinds of scale models are presented. Numerical simulation for the prototype and its 1/40 scale models based on different similitude criteria are carried out. According to the similarity criteria, the result from the scale model is extrapolated to the prototype flow conditions and compared with those directly calculated from the prototype. Based on the viscous fluid dynamics, the effects of viscous force, surface tensional force and buoyant force on the crude oil leakage process are analyzed. The results indicate that the Froude and Reynolds similarity criteria are the dominant similarity criteria and Weber criteria is negligible. When the density of oil used in the model is equal to that of crude, the scale effect is minimized. In different oil leakage stage, different scale model is required to simulate the oil leakage process.

    ship engineering; scale model; similarity criteria; oil leakage; double hull oil tanker; viscous fluid dynamics

    X55;U674.133.1

    A

    猜你喜歡
    油艙縮尺模型試驗(yàn)
    爆炸荷載作用下鋼筋混凝土構(gòu)件縮尺效應(yīng)的數(shù)值模擬研究
    箱梁渦振的縮尺效應(yīng)及振幅修正研究
    尺度效應(yīng)對(duì)噴水推進(jìn)系統(tǒng)進(jìn)出口流場(chǎng)及推力影響分析
    殘油艙水分蒸發(fā)處理方式在IOPP證書(shū)附錄3.2.3項(xiàng)簽注問(wèn)題和記錄注意事項(xiàng)
    反推力裝置模型試驗(yàn)臺(tái)的研制及驗(yàn)證
    臺(tái)階式短加筋土擋墻行為特征的離心模型試驗(yàn)
    土石堆積體介質(zhì)力學(xué)特性縮尺試驗(yàn)研究
    巨厚堅(jiān)硬巖漿巖不同配比的模型試驗(yàn)研究
    電滲—堆載聯(lián)合氣壓劈烈的室內(nèi)模型試驗(yàn)
    撞擊不同貨油艙對(duì)被撞船舷側(cè)結(jié)構(gòu)碰撞性能影響
    国产淫语在线视频| 午夜激情av网站| 在线观看免费午夜福利视频| 国产精品 国内视频| 国产三级黄色录像| 国产91精品成人一区二区三区| 成年动漫av网址| 女人被躁到高潮嗷嗷叫费观| 一级黄色大片毛片| 日韩中文字幕欧美一区二区| 变态另类成人亚洲欧美熟女 | 国产1区2区3区精品| 欧美日韩亚洲国产一区二区在线观看 | 亚洲色图av天堂| 日韩免费高清中文字幕av| 一本大道久久a久久精品| 国产成人精品久久二区二区免费| 亚洲一区二区三区不卡视频| 一本综合久久免费| 好看av亚洲va欧美ⅴa在| 99精品在免费线老司机午夜| 精品国产国语对白av| 一区二区三区国产精品乱码| 精品卡一卡二卡四卡免费| 欧美乱码精品一区二区三区| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区mp4| ponron亚洲| 91字幕亚洲| 亚洲专区中文字幕在线| 国产深夜福利视频在线观看| 99久久人妻综合| 999久久久精品免费观看国产| 建设人人有责人人尽责人人享有的| 国产亚洲精品一区二区www | 国产乱人伦免费视频| 最近最新中文字幕大全电影3 | 欧美亚洲日本最大视频资源| 在线观看舔阴道视频| 91麻豆av在线| 两人在一起打扑克的视频| 最近最新免费中文字幕在线| 人妻 亚洲 视频| 精品国产美女av久久久久小说| 一区二区日韩欧美中文字幕| 午夜影院日韩av| 飞空精品影院首页| 搡老熟女国产l中国老女人| 久久久国产一区二区| 一级,二级,三级黄色视频| 黄频高清免费视频| 日韩视频一区二区在线观看| 亚洲视频免费观看视频| 久久久久精品人妻al黑| 午夜福利免费观看在线| av超薄肉色丝袜交足视频| 一区在线观看完整版| 一级作爱视频免费观看| 久久午夜亚洲精品久久| 99香蕉大伊视频| 极品人妻少妇av视频| 国产男女内射视频| 欧美大码av| 精品国内亚洲2022精品成人 | 美女福利国产在线| 日韩免费高清中文字幕av| 亚洲五月色婷婷综合| 老熟女久久久| 高清黄色对白视频在线免费看| av网站在线播放免费| 成年人午夜在线观看视频| 欧美在线黄色| 一边摸一边抽搐一进一小说 | 国产一区二区激情短视频| 麻豆国产av国片精品| 亚洲色图综合在线观看| 亚洲精品在线观看二区| 成在线人永久免费视频| 欧美一级毛片孕妇| 少妇 在线观看| 黄片小视频在线播放| 在线观看免费日韩欧美大片| 国产精品一区二区免费欧美| 国产精品综合久久久久久久免费 | 欧美+亚洲+日韩+国产| 一个人免费在线观看的高清视频| 老司机午夜福利在线观看视频| 99国产精品一区二区三区| 国产一区在线观看成人免费| 嫁个100分男人电影在线观看| 国产在线精品亚洲第一网站| 悠悠久久av| 十八禁网站免费在线| 国产精品久久久久久人妻精品电影| 高清毛片免费观看视频网站 | 亚洲第一青青草原| 老司机影院毛片| 中国美女看黄片| 日韩有码中文字幕| 国产精品一区二区免费欧美| 久久久久国产精品人妻aⅴ院 | 精品福利永久在线观看| 嫁个100分男人电影在线观看| 久久久精品国产亚洲av高清涩受| 女人精品久久久久毛片| 国产片内射在线| 91在线观看av| 欧美成人午夜精品| 国产亚洲精品第一综合不卡| 久久久久精品人妻al黑| 极品人妻少妇av视频| av不卡在线播放| 欧美黑人精品巨大| 欧美日韩一级在线毛片| 色播在线永久视频| 另类亚洲欧美激情| 国产精品成人在线| 亚洲av成人一区二区三| 午夜福利一区二区在线看| 亚洲五月婷婷丁香| 亚洲欧美一区二区三区久久| 国产不卡一卡二| 亚洲精品自拍成人| 欧美久久黑人一区二区| av中文乱码字幕在线| 日本a在线网址| 男女高潮啪啪啪动态图| 99久久人妻综合| 久久精品国产亚洲av高清一级| 欧美中文综合在线视频| 老司机靠b影院| 亚洲av电影在线进入| 女人被狂操c到高潮| 曰老女人黄片| 黄色视频不卡| 成熟少妇高潮喷水视频| 欧美日韩黄片免| 亚洲av日韩精品久久久久久密| 亚洲一区高清亚洲精品| 激情视频va一区二区三区| 久久久水蜜桃国产精品网| 亚洲人成伊人成综合网2020| 又紧又爽又黄一区二区| 免费黄频网站在线观看国产| 亚洲精品自拍成人| 黄色毛片三级朝国网站| 黄色女人牲交| 久99久视频精品免费| 久久精品亚洲精品国产色婷小说| 美女高潮喷水抽搐中文字幕| 中出人妻视频一区二区| 飞空精品影院首页| 久久亚洲真实| 免费在线观看日本一区| 国产亚洲精品第一综合不卡| 欧美亚洲 丝袜 人妻 在线| 亚洲精品av麻豆狂野| 亚洲成人免费av在线播放| 交换朋友夫妻互换小说| 高潮久久久久久久久久久不卡| 久久人人97超碰香蕉20202| 亚洲精品中文字幕一二三四区| 国产aⅴ精品一区二区三区波| 欧美日韩亚洲高清精品| 啪啪无遮挡十八禁网站| 色综合婷婷激情| 免费少妇av软件| 亚洲五月婷婷丁香| 欧美 亚洲 国产 日韩一| 国产av一区二区精品久久| 久9热在线精品视频| 91九色精品人成在线观看| 欧美黑人精品巨大| 亚洲成a人片在线一区二区| 九色亚洲精品在线播放| 亚洲,欧美精品.| 免费在线观看亚洲国产| 精品国产亚洲在线| 男女之事视频高清在线观看| 中文字幕av电影在线播放| 欧美精品高潮呻吟av久久| 亚洲精品自拍成人| 在线观看舔阴道视频| 少妇被粗大的猛进出69影院| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区mp4| 国产亚洲欧美在线一区二区| 亚洲综合色网址| 侵犯人妻中文字幕一二三四区| 国产欧美亚洲国产| 欧美一级毛片孕妇| 婷婷成人精品国产| av超薄肉色丝袜交足视频| 日日爽夜夜爽网站| 久久天堂一区二区三区四区| 国产麻豆69| 国产精品1区2区在线观看. | 人人妻人人添人人爽欧美一区卜| 午夜免费成人在线视频| 亚洲色图 男人天堂 中文字幕| 亚洲午夜理论影院| 99国产精品免费福利视频| 首页视频小说图片口味搜索| 午夜精品国产一区二区电影| 又大又爽又粗| www日本在线高清视频| 十八禁网站免费在线| 久久草成人影院| 不卡一级毛片| 欧美精品人与动牲交sv欧美| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| svipshipincom国产片| 欧美+亚洲+日韩+国产| 母亲3免费完整高清在线观看| 免费观看精品视频网站| 国产不卡一卡二| 欧美日韩福利视频一区二区| 日韩中文字幕欧美一区二区| 国产精品久久久久久人妻精品电影| 窝窝影院91人妻| 精品视频人人做人人爽| 99热国产这里只有精品6| 下体分泌物呈黄色| 精品国产乱码久久久久久男人| 国产高清国产精品国产三级| 国产真人三级小视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品免费视频内射| 51午夜福利影视在线观看| 看免费av毛片| 精品乱码久久久久久99久播| 多毛熟女@视频| 好男人电影高清在线观看| 久久亚洲真实| 久久精品亚洲熟妇少妇任你| 精品久久久久久久毛片微露脸| 久久久久精品国产欧美久久久| 午夜免费成人在线视频| 操出白浆在线播放| 少妇 在线观看| 国产精品av久久久久免费| 成人18禁高潮啪啪吃奶动态图| av中文乱码字幕在线| 大香蕉久久成人网| 久久国产精品影院| 人成视频在线观看免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品九九99| 国产精品 欧美亚洲| 极品少妇高潮喷水抽搐| 国产欧美日韩综合在线一区二区| 亚洲专区国产一区二区| 国产成人精品无人区| 一级毛片高清免费大全| 国产免费av片在线观看野外av| 亚洲成a人片在线一区二区| 国产无遮挡羞羞视频在线观看| 大型av网站在线播放| 在线观看免费高清a一片| 18禁国产床啪视频网站| 免费久久久久久久精品成人欧美视频| 国产在线一区二区三区精| 久久性视频一级片| 成年人午夜在线观看视频| 国产免费现黄频在线看| 久久热在线av| 18禁国产床啪视频网站| 人妻久久中文字幕网| 国产深夜福利视频在线观看| 久久午夜亚洲精品久久| 午夜激情av网站| 久久中文字幕人妻熟女| 亚洲成人手机| 日韩熟女老妇一区二区性免费视频| 久久ye,这里只有精品| 亚洲国产毛片av蜜桃av| 很黄的视频免费| 久久久精品区二区三区| 久久精品国产亚洲av高清一级| 看片在线看免费视频| 亚洲视频免费观看视频| 成人18禁在线播放| 欧美av亚洲av综合av国产av| 黄色片一级片一级黄色片| 成人特级黄色片久久久久久久| 欧美日韩亚洲综合一区二区三区_| 日日夜夜操网爽| 久久中文看片网| 成熟少妇高潮喷水视频| 女人高潮潮喷娇喘18禁视频| 精品一区二区三区四区五区乱码| 久久人人97超碰香蕉20202| 一区二区三区精品91| 久久久久国产精品人妻aⅴ院 | 亚洲精品一卡2卡三卡4卡5卡| 成在线人永久免费视频| av国产精品久久久久影院| 他把我摸到了高潮在线观看| 国内久久婷婷六月综合欲色啪| 在线观看一区二区三区激情| 十分钟在线观看高清视频www| 悠悠久久av| 亚洲人成电影观看| 国产乱人伦免费视频| 高清毛片免费观看视频网站 | 久久久精品区二区三区| 久久久久久久精品吃奶| 搡老熟女国产l中国老女人| 亚洲精品国产精品久久久不卡| 国产精品一区二区免费欧美| 99热只有精品国产| 国产精品 欧美亚洲| 搡老岳熟女国产| 99热网站在线观看| 国产精品综合久久久久久久免费 | 美女高潮喷水抽搐中文字幕| 国产aⅴ精品一区二区三区波| 亚洲片人在线观看| 老司机午夜福利在线观看视频| 99国产极品粉嫩在线观看| 麻豆国产av国片精品| 亚洲欧美激情在线| cao死你这个sao货| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品大桥未久av| 国产成人免费无遮挡视频| 久久狼人影院| tocl精华| 亚洲精品av麻豆狂野| 女人被躁到高潮嗷嗷叫费观| 丝袜美腿诱惑在线| 亚洲免费av在线视频| 动漫黄色视频在线观看| 免费女性裸体啪啪无遮挡网站| 国产精品免费大片| 国产99白浆流出| 大陆偷拍与自拍| 成人18禁高潮啪啪吃奶动态图| 十八禁高潮呻吟视频| 精品国产超薄肉色丝袜足j| 啦啦啦 在线观看视频| 国产熟女午夜一区二区三区| 精品国产国语对白av| 国产在线观看jvid| 亚洲人成电影观看| 亚洲av片天天在线观看| 99re在线观看精品视频| 日韩大码丰满熟妇| 好看av亚洲va欧美ⅴa在| 亚洲熟妇中文字幕五十中出 | 成人三级做爰电影| 身体一侧抽搐| 丝袜美腿诱惑在线| 午夜久久久在线观看| 欧美日韩亚洲高清精品| 大型av网站在线播放| 日日摸夜夜添夜夜添小说| 美女午夜性视频免费| 麻豆国产av国片精品| 又黄又爽又免费观看的视频| 欧美激情极品国产一区二区三区| 国产精品免费大片| 丝瓜视频免费看黄片| 电影成人av| 看黄色毛片网站| 1024香蕉在线观看| 久久久国产成人精品二区 | 在线观看www视频免费| 中文字幕人妻丝袜一区二区| 在线av久久热| av免费在线观看网站| 日韩熟女老妇一区二区性免费视频| 久久午夜综合久久蜜桃| 巨乳人妻的诱惑在线观看| 99久久国产精品久久久| 精品少妇一区二区三区视频日本电影| 一夜夜www| 国产激情欧美一区二区| a级片在线免费高清观看视频| 成人国语在线视频| 精品国产超薄肉色丝袜足j| 男女免费视频国产| 日日摸夜夜添夜夜添小说| 国产精品香港三级国产av潘金莲| 如日韩欧美国产精品一区二区三区| 9191精品国产免费久久| 国产高清激情床上av| 欧美成狂野欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲综合色网址| 亚洲精品中文字幕在线视频| 性色av乱码一区二区三区2| 高清欧美精品videossex| 叶爱在线成人免费视频播放| 久久久久国产精品人妻aⅴ院 | 久久婷婷成人综合色麻豆| 国产精品久久视频播放| 亚洲av成人一区二区三| 国产av一区二区精品久久| 欧美日韩福利视频一区二区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲男人天堂网一区| 久久久久国产一级毛片高清牌| 亚洲专区字幕在线| 人人妻人人爽人人添夜夜欢视频| 久久久久国产一级毛片高清牌| 一进一出抽搐gif免费好疼 | 人人妻人人澡人人爽人人夜夜| 搡老乐熟女国产| 成年动漫av网址| e午夜精品久久久久久久| 亚洲精品国产区一区二| 亚洲精华国产精华精| 精品国内亚洲2022精品成人 | 亚洲一区二区三区欧美精品| 国产免费男女视频| 午夜福利免费观看在线| 国产精品永久免费网站| 乱人伦中国视频| 日韩大码丰满熟妇| 激情在线观看视频在线高清 | 成人精品一区二区免费| 日韩精品免费视频一区二区三区| 精品熟女少妇八av免费久了| 建设人人有责人人尽责人人享有的| tocl精华| 日本五十路高清| 成年女人毛片免费观看观看9 | 国产区一区二久久| 欧美精品高潮呻吟av久久| 国产激情久久老熟女| 丰满的人妻完整版| 高潮久久久久久久久久久不卡| 美女高潮到喷水免费观看| 在线观看舔阴道视频| 男女免费视频国产| 久久国产精品影院| 一区在线观看完整版| 99久久人妻综合| 亚洲欧美一区二区三区黑人| 欧美日韩精品网址| 欧美日韩视频精品一区| 日韩欧美国产一区二区入口| 日韩三级视频一区二区三区| 亚洲精品久久午夜乱码| 国产精品av久久久久免费| 国产男女内射视频| 国产精品二区激情视频| 黑人操中国人逼视频| 久久99一区二区三区| 国产精品亚洲一级av第二区| av福利片在线| 欧美精品高潮呻吟av久久| 成人18禁在线播放| 国产欧美亚洲国产| 99国产精品一区二区三区| 久久精品91无色码中文字幕| 亚洲精品国产区一区二| 亚洲成a人片在线一区二区| 久久狼人影院| 人人妻人人添人人爽欧美一区卜| 黄色视频,在线免费观看| 色老头精品视频在线观看| 亚洲人成电影观看| 老熟女久久久| 日韩欧美在线二视频 | 亚洲成人手机| 天天操日日干夜夜撸| 午夜两性在线视频| 国产麻豆69| 99re在线观看精品视频| 在线播放国产精品三级| 国产av一区二区精品久久| 亚洲avbb在线观看| 巨乳人妻的诱惑在线观看| e午夜精品久久久久久久| 999久久久国产精品视频| 久久久国产成人免费| 超碰97精品在线观看| 中国美女看黄片| av免费在线观看网站| 亚洲av熟女| 亚洲自偷自拍图片 自拍| 色94色欧美一区二区| 欧美最黄视频在线播放免费 | 免费黄频网站在线观看国产| 老司机午夜十八禁免费视频| 好男人电影高清在线观看| 午夜精品久久久久久毛片777| 在线观看www视频免费| 超色免费av| 侵犯人妻中文字幕一二三四区| a级毛片在线看网站| 国产乱人伦免费视频| 国产成人精品无人区| 免费日韩欧美在线观看| 性色av乱码一区二区三区2| 美女国产高潮福利片在线看| 麻豆国产av国片精品| 日韩成人在线观看一区二区三区| 最近最新中文字幕大全电影3 | 亚洲精品av麻豆狂野| 亚洲精品中文字幕一二三四区| 一区二区三区激情视频| 免费不卡黄色视频| 丰满的人妻完整版| 亚洲中文日韩欧美视频| av超薄肉色丝袜交足视频| ponron亚洲| 成熟少妇高潮喷水视频| 久久久水蜜桃国产精品网| 久久精品亚洲熟妇少妇任你| 后天国语完整版免费观看| 久久国产乱子伦精品免费另类| 欧美精品av麻豆av| 亚洲性夜色夜夜综合| 亚洲五月色婷婷综合| x7x7x7水蜜桃| 丰满的人妻完整版| 午夜91福利影院| 又大又爽又粗| 涩涩av久久男人的天堂| 国产aⅴ精品一区二区三区波| 天天添夜夜摸| 侵犯人妻中文字幕一二三四区| 美女 人体艺术 gogo| 久久精品国产亚洲av香蕉五月 | 亚洲黑人精品在线| 亚洲专区中文字幕在线| 夫妻午夜视频| 露出奶头的视频| 又紧又爽又黄一区二区| 免费不卡黄色视频| 欧美黄色片欧美黄色片| 亚洲自偷自拍图片 自拍| 成年人黄色毛片网站| 欧美av亚洲av综合av国产av| 50天的宝宝边吃奶边哭怎么回事| 18禁美女被吸乳视频| 超碰97精品在线观看| 老司机午夜福利在线观看视频| 天堂动漫精品| 久久久久国产精品人妻aⅴ院 | 啦啦啦免费观看视频1| 脱女人内裤的视频| 啪啪无遮挡十八禁网站| 国产深夜福利视频在线观看| 搡老乐熟女国产| 欧美日韩亚洲综合一区二区三区_| 日韩中文字幕欧美一区二区| 亚洲精品乱久久久久久| 一二三四社区在线视频社区8| 欧美日韩瑟瑟在线播放| 久久久国产成人免费| 精品一品国产午夜福利视频| 精品高清国产在线一区| 最新的欧美精品一区二区| 午夜影院日韩av| 人人澡人人妻人| 天堂中文最新版在线下载| 黄色 视频免费看| 免费久久久久久久精品成人欧美视频| 国产免费现黄频在线看| 麻豆乱淫一区二区| 国产区一区二久久| 一二三四社区在线视频社区8| 又黄又爽又免费观看的视频| 黄色片一级片一级黄色片| 国产欧美日韩一区二区三| 久久精品亚洲精品国产色婷小说| 男人操女人黄网站| 一区福利在线观看| 色播在线永久视频| 91国产中文字幕| 搡老岳熟女国产| 91精品国产国语对白视频| 国产黄色免费在线视频| 亚洲国产精品sss在线观看 | 757午夜福利合集在线观看| 99热国产这里只有精品6| 精品国产美女av久久久久小说| 男男h啪啪无遮挡| 亚洲av第一区精品v没综合| 一边摸一边做爽爽视频免费| 久久香蕉激情| 91成人精品电影| 一级片'在线观看视频| 一进一出好大好爽视频| 淫妇啪啪啪对白视频| 19禁男女啪啪无遮挡网站| 精品卡一卡二卡四卡免费| 黄频高清免费视频| 亚洲专区国产一区二区| 午夜精品在线福利| av线在线观看网站| 大片电影免费在线观看免费| 国产精品久久久人人做人人爽| 一级毛片女人18水好多| 极品教师在线免费播放| 欧美不卡视频在线免费观看 | 女人久久www免费人成看片| 好看av亚洲va欧美ⅴa在| 欧美精品人与动牲交sv欧美| 亚洲精品国产区一区二| 亚洲成a人片在线一区二区| 国产高清激情床上av| 亚洲精品久久午夜乱码| 亚洲九九香蕉| 久热爱精品视频在线9| 久9热在线精品视频| 91成人精品电影| 多毛熟女@视频| 十八禁网站免费在线| 咕卡用的链子|