• <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国产精品国产| 熟女人妻精品中文字幕| 80岁老熟妇乱子伦牲交| 亚洲一区二区三区欧美精品| 夫妻午夜视频| 51国产日韩欧美| 一个人免费看片子| 久久人妻熟女aⅴ| 久久这里只有精品19| 青春草亚洲视频在线观看| 9热在线视频观看99| 日韩制服丝袜自拍偷拍| 老司机影院毛片| 欧美日韩精品成人综合77777| 国产永久视频网站| 久久久久视频综合| 春色校园在线视频观看| 丁香六月天网| 99国产综合亚洲精品| 国产麻豆69| 日本爱情动作片www.在线观看| 国产精品嫩草影院av在线观看| 秋霞在线观看毛片| av视频免费观看在线观看| 亚洲av日韩在线播放| av免费在线看不卡| 热99国产精品久久久久久7| 亚洲少妇的诱惑av| 天堂中文最新版在线下载| 少妇人妻久久综合中文| 一区二区三区乱码不卡18| 三上悠亚av全集在线观看| 久久精品国产亚洲av涩爱| 欧美亚洲 丝袜 人妻 在线| 欧美另类一区| 激情视频va一区二区三区| 美女中出高潮动态图| 狠狠婷婷综合久久久久久88av| 国产精品成人在线| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩卡通动漫| 黄色 视频免费看| 久久久a久久爽久久v久久| 90打野战视频偷拍视频| 伦理电影大哥的女人| 中文欧美无线码| 久久影院123| 99热6这里只有精品| 国产精品国产av在线观看| 母亲3免费完整高清在线观看 | 人妻人人澡人人爽人人| 天堂8中文在线网| 亚洲精品国产色婷婷电影| 久久国产精品男人的天堂亚洲 | 女人被躁到高潮嗷嗷叫费观| 水蜜桃什么品种好| 黄色视频在线播放观看不卡| 久久久久精品性色| 18禁动态无遮挡网站| 91成人精品电影| 午夜福利,免费看| 久久鲁丝午夜福利片| 免费观看在线日韩| 午夜福利网站1000一区二区三区| 国产男女超爽视频在线观看| 最新的欧美精品一区二区| 亚洲一区二区三区欧美精品| 欧美日韩成人在线一区二区| 亚洲国产av新网站| 午夜久久久在线观看| 成人国语在线视频| 国产免费一区二区三区四区乱码| 丝袜美足系列| 麻豆乱淫一区二区| 丝袜喷水一区| 美女中出高潮动态图| 亚洲中文av在线| 亚洲一码二码三码区别大吗| 最黄视频免费看| 精品久久国产蜜桃| 99九九在线精品视频| 曰老女人黄片| 成人二区视频| 久久久a久久爽久久v久久| 80岁老熟妇乱子伦牲交| www.av在线官网国产| 日本免费在线观看一区| 在线看a的网站| 1024视频免费在线观看| 日本wwww免费看| 少妇的丰满在线观看| 午夜激情久久久久久久| 丰满乱子伦码专区| 亚洲精品av麻豆狂野| 五月天丁香电影| 亚洲欧美一区二区三区国产| 热re99久久精品国产66热6| 精品一区二区三区视频在线| 九色亚洲精品在线播放| 老熟女久久久| 久热这里只有精品99| 国产一区二区三区综合在线观看 | 五月玫瑰六月丁香| 久久97久久精品| 国语对白做爰xxxⅹ性视频网站| 久久 成人 亚洲| 欧美日韩视频高清一区二区三区二| 久久人人爽人人爽人人片va| 狂野欧美激情性xxxx在线观看| 免费观看在线日韩| 99热全是精品| 热re99久久国产66热| 热99国产精品久久久久久7| 久久久久久久大尺度免费视频| 亚洲图色成人| 免费在线观看完整版高清| 制服丝袜香蕉在线| 深夜精品福利| 天天躁夜夜躁狠狠躁躁| 十八禁网站网址无遮挡| 色吧在线观看| 亚洲精品国产av成人精品| 一个人免费看片子| 日韩人妻精品一区2区三区| 国产亚洲欧美精品永久| 久久精品国产亚洲av涩爱| 性高湖久久久久久久久免费观看| av天堂久久9| 国产精品一区二区在线不卡| 男男h啪啪无遮挡| 人人澡人人妻人| 亚洲精品aⅴ在线观看| 精品酒店卫生间| 亚洲成人一二三区av| 美女内射精品一级片tv| 婷婷色综合www| 韩国精品一区二区三区 | 日韩人妻精品一区2区三区| 亚洲一级一片aⅴ在线观看| 9色porny在线观看| 精品久久久精品久久久| 久久国产精品男人的天堂亚洲 | 黄色怎么调成土黄色| 99精国产麻豆久久婷婷| av视频免费观看在线观看| 亚洲第一区二区三区不卡| 国产精品一国产av| 夫妻午夜视频| www日本在线高清视频| 久久久a久久爽久久v久久| videosex国产| 久久国产精品大桥未久av| 欧美日韩国产mv在线观看视频| 成人毛片60女人毛片免费| 亚洲美女黄色视频免费看| 欧美精品av麻豆av| 精品一区二区三区视频在线| av黄色大香蕉| h视频一区二区三区| 97超碰精品成人国产| 成人毛片60女人毛片免费| 色吧在线观看| 黄色毛片三级朝国网站| 啦啦啦视频在线资源免费观看| 精品熟女少妇av免费看| 蜜桃在线观看..| 亚洲精品一二三| 涩涩av久久男人的天堂| 国产又色又爽无遮挡免| 人成视频在线观看免费观看| 亚洲综合色惰| av卡一久久| 香蕉精品网在线| a级片在线免费高清观看视频| 欧美97在线视频| 建设人人有责人人尽责人人享有的| 91久久精品国产一区二区三区| 国产精品成人在线| 久久这里只有精品19| 啦啦啦中文免费视频观看日本| 99久国产av精品国产电影| 亚洲av在线观看美女高潮| 国产在线免费精品| 欧美日韩精品成人综合77777| 欧美bdsm另类| 看免费成人av毛片| 新久久久久国产一级毛片| 国产一区有黄有色的免费视频| 欧美97在线视频| 亚洲欧美精品自产自拍| 五月伊人婷婷丁香| 91精品三级在线观看| 午夜福利,免费看| 满18在线观看网站| 国产精品偷伦视频观看了| 少妇人妻 视频| 啦啦啦中文免费视频观看日本| 久久午夜福利片| 成人二区视频| 成年人免费黄色播放视频| 欧美亚洲 丝袜 人妻 在线| 十分钟在线观看高清视频www| 亚洲国产av新网站| 欧美bdsm另类| 亚洲国产av新网站| 久久久久久久久久久久大奶| 99视频精品全部免费 在线| 免费看光身美女| 国产有黄有色有爽视频| 日韩制服骚丝袜av| 国产在视频线精品| 伊人久久国产一区二区| 久久精品人人爽人人爽视色| 国产一区二区在线观看日韩| 人人妻人人添人人爽欧美一区卜| 国产 精品1| 韩国高清视频一区二区三区| 18禁在线无遮挡免费观看视频| www.熟女人妻精品国产 | 五月天丁香电影| 一级黄片播放器| 国产成人欧美| 国产一区二区三区综合在线观看 | 黑人欧美特级aaaaaa片| 大码成人一级视频| 中文精品一卡2卡3卡4更新| 免费久久久久久久精品成人欧美视频 | 国产精品麻豆人妻色哟哟久久| 国产毛片在线视频| 母亲3免费完整高清在线观看 | 好男人视频免费观看在线| 一二三四中文在线观看免费高清| 亚洲综合色网址| 又粗又硬又长又爽又黄的视频| 搡老乐熟女国产| 99热这里只有是精品在线观看| 一区二区日韩欧美中文字幕 | 精品人妻一区二区三区麻豆| 成年av动漫网址| 久久久a久久爽久久v久久| 巨乳人妻的诱惑在线观看| 欧美精品高潮呻吟av久久| 你懂的网址亚洲精品在线观看| 日韩大片免费观看网站| 精品久久久精品久久久| 亚洲四区av| 高清欧美精品videossex| 免费观看性生交大片5| 亚洲精品日本国产第一区| 亚洲人与动物交配视频| 老熟女久久久| 亚洲第一区二区三区不卡| 欧美精品av麻豆av| 蜜桃国产av成人99| 少妇被粗大的猛进出69影院 | 青春草视频在线免费观看| av黄色大香蕉| 男女午夜视频在线观看 | 日韩制服骚丝袜av| 欧美日韩av久久| 国产精品不卡视频一区二区| 欧美精品人与动牲交sv欧美| 99香蕉大伊视频| 国产极品粉嫩免费观看在线| 一本—道久久a久久精品蜜桃钙片| 如日韩欧美国产精品一区二区三区| 国产欧美日韩综合在线一区二区| 毛片一级片免费看久久久久| a级片在线免费高清观看视频| 欧美人与性动交α欧美精品济南到 | 久久ye,这里只有精品| 一区二区av电影网| 久久热在线av| 亚洲国产最新在线播放| 午夜精品国产一区二区电影| 日本与韩国留学比较| 五月天丁香电影| 国产探花极品一区二区| 蜜臀久久99精品久久宅男| 国产精品一国产av| 亚洲天堂av无毛| 女人精品久久久久毛片| 美女xxoo啪啪120秒动态图| 久久国产亚洲av麻豆专区| 婷婷色av中文字幕| 亚洲精品视频女| 国产成人a∨麻豆精品| 国产又爽黄色视频| 十八禁网站网址无遮挡| 国产白丝娇喘喷水9色精品| 亚洲av中文av极速乱| 99香蕉大伊视频| 免费播放大片免费观看视频在线观看| 97在线视频观看| 久久精品国产鲁丝片午夜精品| 亚洲av日韩在线播放| 青春草亚洲视频在线观看| 丝袜在线中文字幕| 亚洲成人av在线免费| 菩萨蛮人人尽说江南好唐韦庄| 久久久国产一区二区| 91精品国产国语对白视频| 一级黄片播放器| 成年女人在线观看亚洲视频| 18在线观看网站| 精品一区在线观看国产| 高清欧美精品videossex| 丝袜在线中文字幕| av有码第一页| 啦啦啦在线观看免费高清www| 国产日韩欧美视频二区| 免费高清在线观看视频在线观看| 2021少妇久久久久久久久久久| av有码第一页| 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| 国产亚洲最大av| 精品一区二区三卡| 亚洲av福利一区| videos熟女内射| 久久精品久久久久久噜噜老黄| 亚洲人与动物交配视频| 22中文网久久字幕| 51国产日韩欧美| 男的添女的下面高潮视频| 老司机亚洲免费影院| 老司机影院毛片| 狂野欧美激情性xxxx在线观看| 久热久热在线精品观看| 国产高清不卡午夜福利| 亚洲情色 制服丝袜| 成人综合一区亚洲| 视频中文字幕在线观看| 亚洲av综合色区一区| 亚洲成人av在线免费| 国产精品女同一区二区软件| 91成人精品电影| av.在线天堂| 1024视频免费在线观看| 激情五月婷婷亚洲| 一级黄片播放器| 欧美变态另类bdsm刘玥| 波野结衣二区三区在线| 99热网站在线观看| 国产成人av激情在线播放| 如何舔出高潮| 人体艺术视频欧美日本| 国产成人a∨麻豆精品| 大香蕉久久网| 国产免费一区二区三区四区乱码| 尾随美女入室| 久久久久人妻精品一区果冻| 宅男免费午夜| 色视频在线一区二区三区| 伦精品一区二区三区| 久久综合国产亚洲精品| 欧美人与善性xxx| 成人18禁高潮啪啪吃奶动态图| 久热久热在线精品观看| 成人毛片60女人毛片免费| av免费在线看不卡| www.熟女人妻精品国产 | 国产精品蜜桃在线观看| 国产精品国产三级国产专区5o| 中文字幕制服av| 欧美日韩综合久久久久久| 日日爽夜夜爽网站| 成人免费观看视频高清| 性色avwww在线观看| 久久久久久久久久人人人人人人| 国产精品免费大片| 18禁动态无遮挡网站| 成年人免费黄色播放视频| 国产精品成人在线| 亚洲情色 制服丝袜| 久久久久网色| 国产不卡av网站在线观看| 久久久亚洲精品成人影院| 97在线视频观看| 国产爽快片一区二区三区| 哪个播放器可以免费观看大片| 国产老妇伦熟女老妇高清| 午夜福利网站1000一区二区三区| 男人爽女人下面视频在线观看| 大陆偷拍与自拍| 国产av国产精品国产| 亚洲人成网站在线观看播放| 免费人妻精品一区二区三区视频| av免费在线看不卡| 亚洲情色 制服丝袜| 国内精品宾馆在线| 成人午夜精彩视频在线观看| 一级,二级,三级黄色视频| 男男h啪啪无遮挡| 亚洲欧美成人精品一区二区| 日韩成人av中文字幕在线观看| 免费人成在线观看视频色| 欧美精品国产亚洲| 亚洲国产看品久久| 最黄视频免费看| 国产老妇伦熟女老妇高清| 少妇精品久久久久久久| 日本黄大片高清| 日韩三级伦理在线观看| 国产成人a∨麻豆精品| 人妻 亚洲 视频| 国产精品国产三级专区第一集| 亚洲内射少妇av| 自线自在国产av| 国产一区二区在线观看日韩| videos熟女内射| 久久久久久伊人网av| 国产精品蜜桃在线观看| 久久人妻熟女aⅴ| 亚洲伊人久久精品综合| videossex国产| 久久精品熟女亚洲av麻豆精品| 插逼视频在线观看| 18禁裸乳无遮挡动漫免费视频| 免费黄色在线免费观看| 制服丝袜香蕉在线| 一二三四在线观看免费中文在 | 成人18禁高潮啪啪吃奶动态图| 97超碰精品成人国产| 天堂8中文在线网| 亚洲精品美女久久av网站| 热re99久久国产66热| 免费观看av网站的网址| 少妇人妻精品综合一区二区| 欧美人与性动交α欧美精品济南到 | 在线精品无人区一区二区三| 欧美亚洲 丝袜 人妻 在线| 午夜免费鲁丝| 狠狠精品人妻久久久久久综合| 欧美xxⅹ黑人| 欧美性感艳星| 丁香六月天网| 69精品国产乱码久久久| 亚洲国产欧美日韩在线播放| 最新中文字幕久久久久| 天天影视国产精品| 男女国产视频网站| 久久精品国产亚洲av涩爱| 亚洲成人手机| 亚洲人成77777在线视频| xxxhd国产人妻xxx| 国产又爽黄色视频| 欧美精品一区二区免费开放| 美女内射精品一级片tv| 考比视频在线观看| 在线观看免费日韩欧美大片| 男男h啪啪无遮挡| 午夜视频国产福利| 极品人妻少妇av视频| 亚洲av综合色区一区| 一区二区三区精品91| 18禁观看日本| xxx大片免费视频| 天天影视国产精品| 日日撸夜夜添| 伊人久久国产一区二区| 国产精品成人在线| 美女国产高潮福利片在线看| 日本-黄色视频高清免费观看| 啦啦啦在线观看免费高清www| 美女国产高潮福利片在线看| 久久久欧美国产精品| 最新中文字幕久久久久| 亚洲人与动物交配视频| 亚洲,欧美,日韩| 亚洲成国产人片在线观看| 午夜激情av网站| 免费大片黄手机在线观看| 免费大片18禁| 亚洲第一av免费看| 夫妻午夜视频| 国产一区二区激情短视频 | 午夜av观看不卡| 久久久久久久精品精品| 少妇精品久久久久久久| 免费在线观看黄色视频的| 2021少妇久久久久久久久久久| 久久婷婷青草| 性高湖久久久久久久久免费观看| 搡女人真爽免费视频火全软件| 亚洲人成77777在线视频| 亚洲天堂av无毛| 丰满饥渴人妻一区二区三| 22中文网久久字幕| 人妻人人澡人人爽人人| 亚洲成色77777| 亚洲婷婷狠狠爱综合网| 天天躁夜夜躁狠狠久久av| 国产成人欧美| 少妇的丰满在线观看| a级片在线免费高清观看视频| 国产精品三级大全| 精品视频人人做人人爽| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品中文字幕在线视频| 国产免费一级a男人的天堂| 亚洲欧美日韩卡通动漫| 亚洲人成网站在线观看播放| 国产成人免费观看mmmm| 成年av动漫网址| 国产精品成人在线| 欧美成人午夜精品| 亚洲精品,欧美精品| 韩国高清视频一区二区三区| 精品一区二区三区四区五区乱码 | 亚洲av中文av极速乱| 蜜桃国产av成人99| 国产黄频视频在线观看| 麻豆精品久久久久久蜜桃| 精品酒店卫生间| 午夜久久久在线观看| 99热全是精品| 插逼视频在线观看| 91精品国产国语对白视频| 亚洲精品美女久久久久99蜜臀 | 久久97久久精品| 肉色欧美久久久久久久蜜桃| 韩国高清视频一区二区三区| 日韩大片免费观看网站| 伦精品一区二区三区| 中文字幕最新亚洲高清| 免费人成在线观看视频色| 成人黄色视频免费在线看| 精品少妇内射三级| 国产亚洲av片在线观看秒播厂| 国产免费一级a男人的天堂| 又黄又粗又硬又大视频| 中文字幕人妻丝袜制服| 亚洲成人手机| 久久久国产精品麻豆| videossex国产| 在线天堂最新版资源| 99久久精品国产国产毛片| 制服诱惑二区| 亚洲国产看品久久| av在线播放精品| 边亲边吃奶的免费视频| 免费黄网站久久成人精品| 国产亚洲欧美精品永久| 久久av网站| 高清不卡的av网站| 国产毛片在线视频| 有码 亚洲区| 久久这里有精品视频免费| 色婷婷av一区二区三区视频| 欧美xxxx性猛交bbbb| 狠狠婷婷综合久久久久久88av| 在线观看国产h片| 国产av一区二区精品久久| √禁漫天堂资源中文www| 久久久a久久爽久久v久久| 亚洲av男天堂| 久久久a久久爽久久v久久| 久久影院123| 精品国产露脸久久av麻豆| 建设人人有责人人尽责人人享有的| 国产极品粉嫩免费观看在线| 精品亚洲乱码少妇综合久久| 久久久久视频综合| 精品一区二区三卡| 高清av免费在线| 精品亚洲成a人片在线观看| 亚洲美女视频黄频| 欧美3d第一页| 蜜臀久久99精品久久宅男| 最近手机中文字幕大全| 大话2 男鬼变身卡| 国产 精品1| 91午夜精品亚洲一区二区三区| 欧美国产精品va在线观看不卡| 午夜福利在线观看免费完整高清在| 老司机亚洲免费影院| 观看av在线不卡| 午夜影院在线不卡| 国产免费视频播放在线视频| 久久久亚洲精品成人影院| 新久久久久国产一级毛片| 色视频在线一区二区三区| 日本免费在线观看一区| 狂野欧美激情性bbbbbb| 欧美bdsm另类| 少妇精品久久久久久久| 久久久精品94久久精品| 亚洲,欧美精品.| 人人澡人人妻人| 久久精品人人爽人人爽视色| 国产在视频线精品| 男的添女的下面高潮视频| 婷婷色av中文字幕| 国产成人免费观看mmmm| 两个人免费观看高清视频| 久久99热这里只频精品6学生| 欧美激情 高清一区二区三区| 色哟哟·www| 夜夜爽夜夜爽视频| 精品一区在线观看国产| 80岁老熟妇乱子伦牲交| 两个人看的免费小视频| 成年人免费黄色播放视频| 午夜91福利影院| 永久网站在线| a级毛色黄片| 免费播放大片免费观看视频在线观看| tube8黄色片| 男男h啪啪无遮挡| 黄片无遮挡物在线观看| 亚洲三级黄色毛片|