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

    一種計算航船對近海海洋環(huán)境噪聲貢獻的算法

    2020-11-10 11:54:16戎澤存胡長青
    聲學技術 2020年5期
    關鍵詞:接收點航船環(huán)境噪聲

    戎澤存,胡長青,趙 梅

    (1.中國科學院聲學研究所東海研究站,上海201815;2.中國科學院大學,北京100049)

    0 引 言

    海洋環(huán)境噪聲作為背景場,在很長一段時間里是被當作干擾項來看待,而事實上海洋環(huán)境噪聲中包含著豐富的信息,因此它具有重要的研究價值。海洋環(huán)境噪聲的組成比較復雜,學界普遍認為航船噪聲和風成噪聲是兩大主要組成部分,前者主要影響的是中低頻段噪聲(10~500 Hz),后者主要影響的是較高頻段噪聲(500~25 000 Hz)。除了以上兩大主要組成外,在低頻段(1~100 Hz),潮、涌、浪、大尺度湍流以及遠處的地震、風暴均對海洋環(huán)境噪聲產生貢獻[1]。

    作為一種隨機過程,海洋環(huán)境噪聲場不是一成不變的,尤其是隨著航運業(yè)的發(fā)展,交通航運噪聲對環(huán)境噪聲的影響更加顯著,使近海海域內的環(huán)境噪聲較之前發(fā)生較大變化,近海環(huán)境噪聲日益成為人們關注的重點。近年來的研究發(fā)現,低頻段海洋環(huán)境噪聲級呈上升趨勢,和經典的Wenz譜相比,在不同頻段內有著小到1~3 dB、大到10~12 dB的增加[2]。

    目前水聲設備向著低頻、遠程方向發(fā)展,作為低頻段海洋環(huán)境噪聲的主要噪聲源,海上航船對噪聲場的影響不容忽視。目前學界對風成海洋噪聲的研究較多,而對于低頻段的海洋環(huán)境噪聲涉足相對較少。因此本文針對淺海環(huán)境,提出一種計算航船噪聲對海洋環(huán)境噪聲貢獻的算法,主要關注頻段為 50~400 Hz,并利用實驗數據對該算法進行了驗證。

    1 理論推導

    1.1 聲傳播模型

    海洋環(huán)境下聲傳播模型主要包括波數積分模型、簡正波模型、射線模型、拋物方程模型、有限差分模型(有限元模型)[3]。本文選用了拋物方程模型對聲場進行計算。

    在研究低頻段海洋環(huán)境噪聲時,對于低頻、淺海條件,射線模型不再適用,因此本文采用了拋物方程(Parabolic Equations,PE)模型[3]。此方法可以較為精確地計算聲壓。

    拋物方程法是一種近似算法,將聲壓方程近似為拋物線型方程。采用柱坐標(r,φ,z)進行推導,等密度介質的三維亥姆霍茲方程為

    假設方位對稱,可得到:

    式中,p(r,z)是聲壓,k0=ω/c0是參考波數,n(r,z)=c0/ c(r,z)是折射率。

    根據Tappert的研究[4],可以假設式(2)的解滿足以下形式:

    將式(3)代入式(2),并作遠場近似,即 k0r?1,可得到:

    引入近軸近似:

    最終得到標準的窄角拋物線方程:

    為求解式(6)所示的PE方程,推導一系列基于算子形式的拋物波方程。首先定義如下算符:

    則式(6)可重寫為

    這里,求解方程的關鍵在于對Q算符的數學近似。距離相關聲學模型(Range-dependent Acoustic Model,RAM)程序作為一種聲場計算程序,它基于分步帕德(Padé)算法,對拋物方程組進行求解。在數學中,Padé近似值是一個給定階數的有理函數對一個函數的“最佳”近似值[9-10]。應用到實際場景時,將聲速梯度、海底聲速、海底底質密度和吸收系數等環(huán)境參數作為輸入,利用RAM程序對聲場進行計算,得到不同頻率聲波在特定海洋環(huán)境下的聲傳播特性。本文將航船視為點源,利用 RAM程序計算淺海中航船的傳播損失LT。在海洋環(huán)境中,聲波在傳播過程中的損失主要包括[5]:擴展損失、吸收損失、散射、海底反射損失。在這里假定海面為理想界面,不考慮海面的散射。

    通常深海環(huán)境下可不考慮海底反射損失。但在淺海環(huán)境下,除了擴展損失和海水吸收損失,海底反射損失是一項關鍵因素。本文關注的是淺海海洋環(huán)境,因此在利用 RAM進行計算時,必須將海底損失考慮在內,以得到聲波傳播過程中的總傳播損失。

    在利用RAM程序計算聲場時,還要注意以下兩點[3]:

    (1)首先要建立一個在深度上有限的求解區(qū)域(0≤z≤zmax),它包括實際物理域(海水和海底)及“海綿”層,“海綿”層是為了避免在下部計算邊界上產生附加反射,減少誤差;“海綿”層底部深度zmax=4H/3,其中H是傳播路徑上實際物理域的最大深度。

    (2)其次要選擇合適的網格尺寸(Δr,Δz),Δr為距離步長,Δz為深度步長。在淺海環(huán)境中,深度步長需滿足:Δz≤λ/4,其中λ是聲波波長;距離步長需滿足:2Δz≤Δr≤5Δz。為了獲得精確的數值解,唯一的辦法是按一定規(guī)則逐步減小Δr和Δz,對解(傳播損失)的收斂性仔細的檢驗。若Δr和Δz的值設置不合適,可能會導致不同頻率下的傳播損失計算值產生誤差,因此需要針對不同頻率分別進行收斂性檢驗。

    1.2 航船噪聲計算方法

    為了研究海洋環(huán)境噪聲場的水平方向分布,以接收點為中心,將海洋劃分為若干楔形區(qū)域,如圖1(a)所示。圖1中xOy平面代表海面;從海面開始往下直至接收點所在平面,接收點記為R,D為接收點的深度。在z=D處取垂直于z軸的截面,并繪制扇區(qū)分布如圖1(b)所示。航船位置按劃分的扇區(qū)分組,扇區(qū)均勻分布,每個扇區(qū)作為一個方向。例如,劃分間隔若為2°,整個區(qū)域將劃分為180個扇區(qū)。扇區(qū)劃分越密,誤差越小。扇區(qū)截面的半徑由選取的計算范圍決定,選取的范圍越大,半徑越大,黑色散點代表航船在扇區(qū)截面的投影。利用本文算法計算時,可提前計算目標海域的聲場,得到傳播損失隨傳播距離的變化,結合傳播損失和航船聲源級確定扇區(qū)的半徑。

    圖1 楔形區(qū)域及扇區(qū)劃分分布圖Fig.1 Division map of wedge-shaped area and sector area

    對某一個扇區(qū)的接收聲壓pR,可按式(9)和式(10)進行計算[6]:

    式(10)中,m表示航船序號,m=1,2,3,··,M,M為該扇區(qū)內船的數量,各個符號含義如下:

    (1)rm為該扇區(qū)內序號為m的航船到接收點的距離;

    (2)rc為該扇區(qū)內航船與接收點的最大距離,分別計算該扇區(qū)內的每艘航船到接收點的距離,取最大值即為rc;

    (3)Lsm是該扇區(qū)內航船 m在某一頻率下的聲源級,計算方法在1.3節(jié)討論;

    (4)ωm指該扇區(qū)內航船m對應的權重系數,由式(11)得到:

    按式(9)、(10)計算得到不同頻率的接收聲壓pR后,進而可得到接收點處水聽器的接收聲壓譜密度DRL(f)[7],計算公式為

    式中:參考聲壓pref=1μ Pa。對其他扇區(qū)可重復進行上述計算流程。

    1.3 航船噪聲源

    對航船噪聲源的模型,目前仍缺乏合適的數學與物理描述,一般仍采用經驗公式或半經驗公式描述。Ross[1]提出典型商船輻射噪聲的回歸公式,后來Hamson[7]更新了Ross的航船輻射噪聲的回歸方程,航船輻射噪聲功率譜密度如式(14)[1,7]所示:

    式(14)各參數代表的意義如下:

    (1)cv和cl分別代表航船速度和航船長度的羅斯冪律系數,cv=6,cl=2。

    (2)v,l分別是船舶的速度和長度。v0和l0分別是航船參考速度和參考長度,v0=12 kn(1kn=1.852 km·h-1)、l0=300 ft(1ft=0.3048 m)。

    (3)df是隨頻率變化的修正量,由式(15)給出,頻率f的單位是Hz[7]:

    (4)dl是隨航船長度變化的修正量,由式(16)給出:

    (5)Ls0(f)是僅與頻率相關的分量,由式(17)給出[7]:

    利用式(14),結合獲取的航船自動識別系統(tǒng)(Automatic Identification System,AIS)數據(AIS數據將在2.1節(jié)進行說明),可以得到某一頻率的航船輻射噪聲功率譜密度Ls。

    2 AIS數據介紹及預處理

    2.1 AIS數據介紹

    利用第1節(jié)構建的算法,結合AIS數據,進行了仿真計算。航船的位置、航速、船長等信息來源于船舶自動識別系統(tǒng)(AIS),以上數據作為算法的輸入,需要進行一定的預處理,在本節(jié)中將對此進行闡述。

    圖2所示的是仿真所考察的航船AIS數據分布范圍示意圖,其中中心點經緯度為[123.16°E,30.02°N],4個頂點的坐標分別為:[122.16°E,29.02°N]、[124.16°E,29.02°N]、[124.16°E,31.02°N]、[122.16°E,31.02°N]。該區(qū)域經度跨度、緯度跨度均為2°。

    航船AIS數據主要包括船名、海事移動服務識別參數(Maritime Mobile Service Identity,MMSI)、船舶類型、船舶狀態(tài)、船長、船寬、吃水、航向、速度、經度、緯度、時間戳等信息。文中船名和MMSI號已被隱去,必要時以手動編號代替。

    圖2 航船AIS數據分布范圍Fig.2 Distribution range of ship AIS data

    表1中代表的是某航船在兩個時刻對應的航船信息,包括船長、船速、吃水深度、經緯度和記錄時間等,在兩個時間點的航速和航向不同。表中UTC代表格林威治標準時間,與北京時間相差8個小時,“UTC+8”即為北京時間。

    表1 某航船AIS數據Table 1 AIS data of a ship

    2.2 數據預處理

    對于得到的AIS數據進行預處理,包括:

    (1)剔除存在明顯錯誤的數據,剔除航船中發(fā)動機已關閉的船舶;

    (2)篩選符合時間范圍內的數據,利用 Matlab軟件保存到數據庫;

    (3)利用航船的船速、航向對數據進行插值處理,作為數據的補充;

    (4)查看地圖可發(fā)現,島礁是一個必須考慮的因素,在某一個時刻下,如果航船與接收點的連線上出現了島礁,可能會出現衍射(繞射)?,F已證明,聲衍射的強弱和障礙物的大小與聲波波長的比值密切相關,其關系式為

    式中:k為聲波的波數;f、c、λ分別為聲波的頻率、速度和波長;a是代表障礙物尺度的量,如為球體a就代表球體的半徑。當ka?1時衍射很弱,聲波幾乎只是沿著直線傳播以至于被球體擋住傳播路徑而在其背面形成明顯的聲影。當島礁尺寸滿足ka?1時,則應將此艘航船予以剔除。

    經預處理得到航船分布和各艘船的信息后,依據算法,可獨立計算各扇區(qū)內船舶輻射噪聲級譜密度DRL(f)。

    3 仿真計算

    3.1 仿真水體環(huán)境

    水體環(huán)境參數參考某次海洋環(huán)境噪聲實驗。該次測量點位于圖2所示海域范圍內,水深約70 m,聲速梯度由電導率、溫度、深度(Conductivity,Temperature,Depth,CTD)系統(tǒng)測得,如圖3所示。用側掃聲吶測量海底,發(fā)現該海域海底較為平坦。同時通過對海底底質的采樣,可知其主要成分接近粉砂質砂。

    圖3 實驗海域聲速梯度剖面圖Fig.3 Acoustic velocity profile in the experimental area

    3.2 水平方向分布分析

    取某一時刻海域內的航船數據作為輸入,包括航船經緯度、長度、速度、吃水深度。預處理之后,運用1.3節(jié)的方法進行仿真計算。圖4即為該時刻預處理之后的航船位置分布圖,將此位置分布下的航船數據作為輸入,進行航船噪聲的仿真計算。

    按1/3倍頻程計,選定中心頻率,分別計算頻段內若干頻點的功率譜密度值,并取平均得到的結果如圖5所示。

    從圖5中可以看出:

    (1)在不同頻率下,航船噪聲功率譜密度水平分布圖形狀不同。這是由于不同頻率的聲源在實驗海域內的傳播損失曲線不同:給定兩個不同頻率、相同聲源級的聲源,兩者的傳播損失隨距離的變化不同,導致不同頻率下的水平分布圖不同。

    圖4 某一時刻的航船位置分布圖Fig.4 Distribution of ship positions at a certain time

    圖5 不同中心頻率航船噪聲功率譜密度的水平方向分布圖Fig.5 Horizontal distributions of the ship noise power spectral density at different center frequencies

    (2)圖5(a)~5(j)中:6°~12°、32°~36°、50°~52°、70°~72°、152°~154°、-152°~-122°、-118°~-30°、-18°~-16°和-12°~-8°方向的水平分布圖為“包絡”,其余為“線譜”。

    進一步分析數據,計算方法不變,從圖5(a)~5(h)中提取相應方向的數據形成的-116°方向、128°方向、-36°方向、34°方向上在這一時刻的航船信息,結果如表2所示。表2說明在這一時間點某一方向上有幾艘船,以及每艘船的船長、船速、吃水深度信息,在-116°方向上有兩艘船:A1,A2;在-36°方向上有四艘船:B1,B2,B3,B4;而在128°方向和34°方向上均只有一艘船,分別是D和E。

    4個方向上航船噪聲功率譜密度隨頻率的變化趨勢如圖6(a)~6(b)所示。

    依頻率提取多個方向上的數據,得出如下結論:

    (1)若扇區(qū)內存在多艘航船,各艘航船的航行狀態(tài)不同,受到航船噪聲疊加的影響,不同方向上航船噪聲功率譜密度隨頻率的變化趨勢差別較大,如圖6(a)、6(b)所示。

    (2)若扇區(qū)內僅存在一艘航船,則航船噪聲功率譜密度隨頻率的變化趨勢較為穩(wěn)定,圖形走勢基本一致,如圖6(c)、6(d)所示。

    表2 不同方向上的航船信息Table 2 Ship information in different directions

    3.3 功率譜密度分析

    采用本文算法,航船AIS信息輸入和水體環(huán)境參數與3.1節(jié)相同,得到了某一時刻的航船噪聲功率譜密度仿真結果。該時刻為2016年7月6日9時0分0秒,對不同方向的噪聲功率譜密度取平均,仿真結果如圖7中的星點所示,計算的頻點包括50、63、80、100、125、160、200、250、315、400 Hz。

    圖6 不同方向航船噪聲功率譜密度Fig.6 The ship noise power spectrum densities in different directions

    圖7 對某一時刻不同方向航船噪聲功率譜密度取平均的仿真結果Fig.7 Simulation result of the ship noise power spectral density averaged in different directions at a certain time

    對于海洋環(huán)境噪聲的研究,Wenz噪聲譜級圖是一個重要的參考。它是由Wenz等[8]總結出的經典海洋環(huán)境噪聲功率譜級圖。經典Wenz噪聲譜級圖中包含了航船、風、降雨、地震和爆炸等不同因素對環(huán)境噪聲的貢獻,如圖8所示:橫坐標代表頻率(以倍頻程記,頻率從1 Hz~100 kHz),縱坐標代表海洋環(huán)境噪聲譜級。Wenz噪聲譜級圖中,不同因素貢獻的噪聲頻段不同,圖中各曲線代表的影響因素如表 3 所示。航船對噪聲的貢獻主要是10 Hz~1 kHz頻段,圖8中描述了強航道、淺海通常航道、深海通常航道等條件下的噪聲譜級。

    圖8 Wenz噪聲功率譜級圖[8]Fig.8 Wenz noise power spectrum diagram[8]

    為了驗證本文算法的可靠性,將Wenz噪聲譜級圖中描述淺海通常航道的曲線取出,與圖7所示的仿真計算值同時繪制在圖9中。圖9中Wenz曲線 1、2分別代表淺海通常航道(Usual Traffic Noise-Shallow Water)的上、下限曲線。由圖9可以看出,仿真計算值與描述淺海通常航道的曲線2較為符合,在淺海通常航道的上、下限曲線之間,這也與水聽器和航船距離較遠這一事實符合。實驗所用水聽器為無指向性水聽器,布放位置為圖4中的接收點位置,深度為 41.03 m,水聽器與航船的距離從36.69~147.46 km不等。

    表3 Wenz噪聲功率譜級圖說明Table 3 Description of Wenz noise power spectrum diagram

    為了進一步驗證算法的可靠性,對實驗數據進行了分析。圖10中連續(xù)曲線代表的是在2016年7月6日8時59分30秒~9時30秒的時間段內,實測數據的1/3倍頻程功率譜密度。由圖10可見,在100~400 Hz頻段內,仿真值與實驗值較為符合。

    與圖 10不同,現選取 10個時間點(間隔為1 min),對10個時間點重復上述計算,并取平均,結果如圖11中星點所示。圖11中的實驗值是對共計10 min長度的實驗噪聲數據進行分析得到的。由圖10、11可以發(fā)現,其結果與圖10中星點代表的仿真計算值結果較為一致,在100~400 Hz頻段內兩者較為符合。

    圖9 某一時刻航船噪聲功率譜密度仿真結果(圖7)與典型淺海航道的Wenz譜曲線對比Fig 9 Comparison between the simulation result of ship noise power spectral density(Fig.7)and the Wenz spectrums in typical shallow sea channels

    圖10 某一時刻航船噪聲功率譜密度仿真結果(圖7)與實驗值對比Fig.10 Comparison between the simulation result of ship noise power spectral density(Fig.7)and the experimental values

    圖11 航船噪聲功率譜密度仿真值(某一時間段內取平均)與實驗值對比Fig 11 Comparison between the simulation result of ship noise power spectral density(averaging over a period of time)and the experimental value

    本文中 50~400 Hz頻段內的仿真結果與經典Wenz譜總體上較為符合,在Wenz曲線淺海通常航道的上下限曲線之內。但50~100 Hz頻段內的仿真結果略低于實測環(huán)境噪聲數據。分析實驗海域附近環(huán)境發(fā)現,由于實驗海域為近海,環(huán)境噪聲源種類繁多,噪聲源除了航船以外,還有工業(yè)活動、潮、涌、浪、遠處風暴等,并且這些聲源對低頻段的影響較大,這也是誤差產生的原因之一。

    4 結 論

    本文提出了一種計算航船對近海海洋環(huán)境噪聲貢獻的算法。利用本文算法,結合AIS數據和實際海域水文參數,可獲取接收點處航船噪聲的水平方向分布特點,并可初步定量分析航船噪聲對海洋環(huán)境噪聲的貢獻。本文主要關注的頻段為 50~400Hz,理論分析和數值計算均反映出了由于航船的非均勻分布,噪聲場在不同頻率下呈現明顯的非均勻分布特征這一特點。對噪聲功率譜密度進行了定量分析,并通過與經典Wenz譜和實測噪聲數據進行比對,在約50~400 Hz時仿真值和實驗值吻合較好,驗證了算法的可靠性。同時,在應用時需要盡可能多地獲取環(huán)境參數,以降低數值計算時的誤差。

    猜你喜歡
    接收點航船環(huán)境噪聲
    南湖航船指方向
    環(huán)境噪聲智能檢測終端設計
    電子制作(2018年23期)2018-12-26 01:01:04
    基于云服務的環(huán)境噪聲智能監(jiān)測系統(tǒng)研究
    測控技術(2018年11期)2018-12-07 05:49:22
    更正
    動態(tài)網絡最短路徑射線追蹤算法中向后追蹤方法的改進*1
    南海某海域環(huán)境噪聲與風速相關性研究*
    中國港航船企指數(CMEI)
    鐵路環(huán)境噪聲測量中的問題分析
    淺海波導界面對點源振速方向的影響?
    應用聲學(2015年3期)2015-10-27 02:52:49
    2014中國港航船企收入榜單
    欧美少妇被猛烈插入视频| av又黄又爽大尺度在线免费看| 色婷婷av一区二区三区视频| 久久久久久久精品精品| 另类精品久久| 不卡视频在线观看欧美| 午夜av观看不卡| 国产成人啪精品午夜网站| 毛片一级片免费看久久久久| 亚洲国产欧美日韩在线播放| 精品人妻在线不人妻| 丝瓜视频免费看黄片| 婷婷色av中文字幕| 男女无遮挡免费网站观看| 最近中文字幕2019免费版| 男人操女人黄网站| 亚洲精品久久午夜乱码| 纯流量卡能插随身wifi吗| 亚洲男人天堂网一区| 久久精品国产亚洲av高清一级| 性少妇av在线| 日韩免费高清中文字幕av| 嫩草影视91久久| 日本黄色日本黄色录像| 久久精品久久久久久噜噜老黄| 国产精品女同一区二区软件| 美女视频免费永久观看网站| 九草在线视频观看| 观看av在线不卡| 9色porny在线观看| 国产精品久久久久久久久免| 国产欧美亚洲国产| 国产麻豆69| 丁香六月欧美| 在线亚洲精品国产二区图片欧美| 日韩熟女老妇一区二区性免费视频| 免费观看a级毛片全部| 亚洲国产精品一区三区| 亚洲欧美激情在线| 亚洲精品美女久久久久99蜜臀 | 观看av在线不卡| 亚洲情色 制服丝袜| 侵犯人妻中文字幕一二三四区| 男人添女人高潮全过程视频| 日韩av在线免费看完整版不卡| 2018国产大陆天天弄谢| 精品少妇一区二区三区视频日本电影 | 国产成人午夜福利电影在线观看| 色精品久久人妻99蜜桃| 黄片播放在线免费| 性色av一级| 亚洲欧美成人精品一区二区| 一区在线观看完整版| 国产高清国产精品国产三级| 精品国产超薄肉色丝袜足j| 蜜桃在线观看..| 日韩人妻精品一区2区三区| 亚洲精品国产一区二区精华液| 一本大道久久a久久精品| 新久久久久国产一级毛片| 亚洲av电影在线进入| 久久精品久久久久久噜噜老黄| 国产1区2区3区精品| 午夜精品国产一区二区电影| 黑丝袜美女国产一区| 操美女的视频在线观看| 欧美国产精品va在线观看不卡| tube8黄色片| 久久99一区二区三区| 啦啦啦在线观看免费高清www| 亚洲四区av| 精品一区二区三卡| 永久免费av网站大全| 成人三级做爰电影| 免费黄色在线免费观看| 久久久久久久精品精品| 国产亚洲一区二区精品| 五月开心婷婷网| 国产欧美日韩综合在线一区二区| 91精品伊人久久大香线蕉| 国产爽快片一区二区三区| 国产一区二区激情短视频 | 亚洲精品自拍成人| 制服人妻中文乱码| 如何舔出高潮| 精品国产一区二区三区久久久樱花| 色婷婷av一区二区三区视频| 国产精品久久久久成人av| 一级爰片在线观看| 亚洲精品日韩在线中文字幕| av天堂久久9| 国产精品一二三区在线看| 亚洲,欧美精品.| 青草久久国产| 91精品伊人久久大香线蕉| 97在线人人人人妻| 中文天堂在线官网| 9色porny在线观看| 久久精品国产a三级三级三级| 黑丝袜美女国产一区| 狂野欧美激情性xxxx| 老司机亚洲免费影院| 精品国产一区二区三区四区第35| 欧美亚洲 丝袜 人妻 在线| 成人18禁高潮啪啪吃奶动态图| 青春草国产在线视频| 高清视频免费观看一区二区| 国产一区二区在线观看av| av电影中文网址| 国产成人系列免费观看| 久久久久久久久久久久大奶| 天美传媒精品一区二区| 男女午夜视频在线观看| a级毛片在线看网站| 精品第一国产精品| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久人妻精品电影 | 久久久久久人人人人人| 你懂的网址亚洲精品在线观看| 精品久久蜜臀av无| 国产乱来视频区| 欧美人与善性xxx| 在线观看免费午夜福利视频| 亚洲四区av| 久久毛片免费看一区二区三区| 91aial.com中文字幕在线观看| 国产99久久九九免费精品| 久久久久久久大尺度免费视频| 免费高清在线观看日韩| 欧美中文综合在线视频| 免费看av在线观看网站| 啦啦啦 在线观看视频| 亚洲av成人精品一二三区| 一区二区三区乱码不卡18| 国产又爽黄色视频| 国产一区二区三区av在线| 久久97久久精品| 亚洲精品一区蜜桃| 国产av码专区亚洲av| 日本猛色少妇xxxxx猛交久久| 中文乱码字字幕精品一区二区三区| 久久久精品免费免费高清| 一区二区三区精品91| 一区二区三区四区激情视频| 人人妻,人人澡人人爽秒播 | 少妇被粗大猛烈的视频| 下体分泌物呈黄色| 不卡视频在线观看欧美| 久久久国产一区二区| 丰满饥渴人妻一区二区三| 欧美激情极品国产一区二区三区| 国产亚洲午夜精品一区二区久久| 天天躁夜夜躁狠狠躁躁| 中文天堂在线官网| 日韩成人av中文字幕在线观看| 成年av动漫网址| 免费观看人在逋| 日本猛色少妇xxxxx猛交久久| 啦啦啦啦在线视频资源| 天堂中文最新版在线下载| 国产人伦9x9x在线观看| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲综合一区二区三区_| 中文精品一卡2卡3卡4更新| 侵犯人妻中文字幕一二三四区| 黄色视频不卡| 免费久久久久久久精品成人欧美视频| 伊人久久大香线蕉亚洲五| 免费女性裸体啪啪无遮挡网站| 国产精品欧美亚洲77777| 亚洲精品在线美女| 国产成人午夜福利电影在线观看| 国产女主播在线喷水免费视频网站| 2018国产大陆天天弄谢| 97人妻天天添夜夜摸| 两个人免费观看高清视频| 亚洲第一区二区三区不卡| 亚洲第一区二区三区不卡| 亚洲美女视频黄频| 制服人妻中文乱码| a 毛片基地| 搡老岳熟女国产| 国产在线免费精品| 国产精品偷伦视频观看了| 多毛熟女@视频| 亚洲欧美一区二区三区久久| 丝袜在线中文字幕| 两个人免费观看高清视频| 久久久久久久久久久久大奶| 国产97色在线日韩免费| 久久97久久精品| 超碰97精品在线观看| 女人久久www免费人成看片| av又黄又爽大尺度在线免费看| 亚洲精品在线美女| 满18在线观看网站| 国产免费一区二区三区四区乱码| 一本一本久久a久久精品综合妖精| 日本wwww免费看| a 毛片基地| 美女高潮到喷水免费观看| xxx大片免费视频| 国产精品一二三区在线看| 99精国产麻豆久久婷婷| 欧美日韩视频高清一区二区三区二| 日本91视频免费播放| 国产黄色免费在线视频| 国产在视频线精品| 日韩,欧美,国产一区二区三区| 中文字幕色久视频| 中文字幕制服av| 亚洲国产精品一区二区三区在线| 国产1区2区3区精品| 午夜福利乱码中文字幕| 99香蕉大伊视频| 国产av精品麻豆| 观看av在线不卡| 久久久久精品国产欧美久久久 | 亚洲中文av在线| 极品人妻少妇av视频| 中文字幕精品免费在线观看视频| 久久久欧美国产精品| 国产国语露脸激情在线看| 久久精品亚洲熟妇少妇任你| 久久国产精品大桥未久av| 男女高潮啪啪啪动态图| 国产精品久久久av美女十八| 亚洲国产欧美一区二区综合| 免费不卡黄色视频| 一级毛片电影观看| 久久国产亚洲av麻豆专区| 国产日韩欧美视频二区| 亚洲成av片中文字幕在线观看| 亚洲免费av在线视频| 青春草国产在线视频| 色网站视频免费| 国产片内射在线| 精品亚洲乱码少妇综合久久| 午夜av观看不卡| av女优亚洲男人天堂| videosex国产| 亚洲精品久久午夜乱码| 国产精品久久久人人做人人爽| 久久久久精品久久久久真实原创| 青春草国产在线视频| 欧美xxⅹ黑人| 王馨瑶露胸无遮挡在线观看| 久久久久精品人妻al黑| 成人影院久久| 日韩电影二区| e午夜精品久久久久久久| 亚洲熟女毛片儿| 在线天堂最新版资源| 精品一区二区三区av网在线观看 | 在线观看免费日韩欧美大片| 国产片内射在线| 韩国高清视频一区二区三区| 少妇被粗大的猛进出69影院| 黄色视频不卡| 美女主播在线视频| 亚洲精品久久午夜乱码| 中文字幕制服av| 亚洲精品美女久久av网站| 男女边摸边吃奶| 日韩中文字幕视频在线看片| 久久人人爽人人片av| 岛国毛片在线播放| 女人久久www免费人成看片| 极品少妇高潮喷水抽搐| 亚洲七黄色美女视频| 久久久久精品久久久久真实原创| 午夜福利视频精品| 国产野战对白在线观看| 在线看a的网站| 亚洲av福利一区| 色婷婷久久久亚洲欧美| 蜜桃在线观看..| 如日韩欧美国产精品一区二区三区| www.熟女人妻精品国产| 欧美日韩精品网址| 赤兔流量卡办理| 精品久久久久久电影网| 黄色 视频免费看| 久久久久久久久久久久大奶| 国产片特级美女逼逼视频| 国产男女超爽视频在线观看| 亚洲,一卡二卡三卡| 国产av一区二区精品久久| 各种免费的搞黄视频| 午夜精品国产一区二区电影| 久久性视频一级片| 国产亚洲av高清不卡| 51午夜福利影视在线观看| 大片免费播放器 马上看| 午夜福利,免费看| 欧美国产精品一级二级三级| 久久精品aⅴ一区二区三区四区| 午夜老司机福利片| 韩国精品一区二区三区| 哪个播放器可以免费观看大片| 秋霞在线观看毛片| 又大又爽又粗| 国产日韩欧美视频二区| 免费高清在线观看视频在线观看| 一边摸一边做爽爽视频免费| 久久久久久久久久久久大奶| 水蜜桃什么品种好| 交换朋友夫妻互换小说| 一区二区三区激情视频| 90打野战视频偷拍视频| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 黑人欧美特级aaaaaa片| 欧美黑人欧美精品刺激| 亚洲精品国产色婷婷电影| 国产又爽黄色视频| 又粗又硬又长又爽又黄的视频| 男女午夜视频在线观看| 日韩不卡一区二区三区视频在线| 黄色 视频免费看| 在线观看免费日韩欧美大片| 久久久精品区二区三区| 只有这里有精品99| 啦啦啦啦在线视频资源| 成人黄色视频免费在线看| 亚洲精品成人av观看孕妇| 久久久久久久久久久久大奶| 亚洲欧美精品综合一区二区三区| 亚洲精品视频女| 久久久久精品性色| 最近最新中文字幕大全免费视频 | 精品少妇一区二区三区视频日本电影 | 成人国产av品久久久| 成人手机av| 狂野欧美激情性bbbbbb| 亚洲国产欧美一区二区综合| 亚洲欧美一区二区三区国产| 老鸭窝网址在线观看| 99久国产av精品国产电影| 一级毛片黄色毛片免费观看视频| 大片免费播放器 马上看| 国产精品.久久久| 啦啦啦中文免费视频观看日本| 亚洲第一青青草原| 亚洲国产精品一区三区| 黄片小视频在线播放| 少妇人妻久久综合中文| 久久人人爽人人片av| 伦理电影大哥的女人| 色吧在线观看| 伊人久久国产一区二区| 日本wwww免费看| 自线自在国产av| 99香蕉大伊视频| 亚洲国产欧美日韩在线播放| 国产欧美日韩综合在线一区二区| 18禁观看日本| 日韩不卡一区二区三区视频在线| 韩国精品一区二区三区| 七月丁香在线播放| 巨乳人妻的诱惑在线观看| 久久久久国产一级毛片高清牌| 午夜影院在线不卡| 色综合欧美亚洲国产小说| 国产免费又黄又爽又色| 2018国产大陆天天弄谢| 亚洲人成网站在线观看播放| 国产1区2区3区精品| 啦啦啦 在线观看视频| 天堂8中文在线网| 亚洲视频免费观看视频| 精品第一国产精品| 精品视频人人做人人爽| 中文欧美无线码| 9191精品国产免费久久| av视频免费观看在线观看| 99久久精品国产亚洲精品| 久久人人爽人人片av| 丝袜美腿诱惑在线| 久久久久精品国产欧美久久久 | 在线观看国产h片| 久久毛片免费看一区二区三区| 精品一区在线观看国产| 久久亚洲国产成人精品v| 久久久久人妻精品一区果冻| 精品少妇黑人巨大在线播放| 毛片一级片免费看久久久久| 欧美激情高清一区二区三区 | 99国产综合亚洲精品| 欧美日韩国产mv在线观看视频| 人体艺术视频欧美日本| 亚洲精品日本国产第一区| 最近最新中文字幕大全免费视频 | 两个人免费观看高清视频| 激情五月婷婷亚洲| 亚洲,欧美精品.| 久久久久人妻精品一区果冻| 一本一本久久a久久精品综合妖精| 国产日韩一区二区三区精品不卡| 男女床上黄色一级片免费看| 欧美久久黑人一区二区| 日本黄色日本黄色录像| av国产久精品久网站免费入址| 香蕉丝袜av| 欧美人与性动交α欧美软件| 天天躁夜夜躁狠狠久久av| 曰老女人黄片| 日韩av在线免费看完整版不卡| 亚洲三区欧美一区| 日韩伦理黄色片| 久久精品国产a三级三级三级| 女人精品久久久久毛片| h视频一区二区三区| 国产精品欧美亚洲77777| 久久天堂一区二区三区四区| 我的亚洲天堂| 看十八女毛片水多多多| 超碰成人久久| 又大又爽又粗| 久久久久久久大尺度免费视频| 久久久精品国产亚洲av高清涩受| 午夜影院在线不卡| 日本猛色少妇xxxxx猛交久久| 男人舔女人的私密视频| 青草久久国产| 国产日韩欧美视频二区| 精品国产露脸久久av麻豆| 最近中文字幕高清免费大全6| 欧美黑人精品巨大| 99热国产这里只有精品6| 人人妻,人人澡人人爽秒播 | 一边亲一边摸免费视频| 欧美黄色片欧美黄色片| xxxhd国产人妻xxx| 国产精品av久久久久免费| 在线 av 中文字幕| 一区二区三区精品91| 久久久久久久大尺度免费视频| 精品一区二区免费观看| avwww免费| 亚洲国产精品一区三区| 国产精品.久久久| av电影中文网址| 欧美精品一区二区大全| 日韩中文字幕欧美一区二区 | 日日爽夜夜爽网站| 叶爱在线成人免费视频播放| 日韩 亚洲 欧美在线| 香蕉国产在线看| 午夜精品国产一区二区电影| 中文字幕色久视频| 亚洲成人国产一区在线观看 | 亚洲专区中文字幕在线 | 国产伦理片在线播放av一区| 丁香六月欧美| 久久精品国产亚洲av高清一级| 欧美精品亚洲一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 天天躁夜夜躁狠狠躁躁| 在线观看免费午夜福利视频| 晚上一个人看的免费电影| 女性被躁到高潮视频| 亚洲成人av在线免费| 黄色怎么调成土黄色| 男女无遮挡免费网站观看| 成人国产av品久久久| 精品一区二区三区av网在线观看 | 一区二区av电影网| 新久久久久国产一级毛片| 2018国产大陆天天弄谢| 亚洲精品美女久久久久99蜜臀 | 国产一区有黄有色的免费视频| 国产男女内射视频| 国产一区二区三区综合在线观看| 亚洲自偷自拍图片 自拍| 亚洲激情五月婷婷啪啪| 女人被躁到高潮嗷嗷叫费观| 久久久久国产一级毛片高清牌| 蜜桃在线观看..| 卡戴珊不雅视频在线播放| 亚洲欧美精品自产自拍| 亚洲欧美激情在线| 精品一区二区三区四区五区乱码 | 国产欧美亚洲国产| 热99国产精品久久久久久7| 麻豆av在线久日| 999久久久国产精品视频| 大码成人一级视频| 欧美人与善性xxx| 欧美乱码精品一区二区三区| 日韩av在线免费看完整版不卡| 18禁动态无遮挡网站| 最近最新中文字幕免费大全7| 亚洲人成网站在线观看播放| 国产成人a∨麻豆精品| 永久免费av网站大全| 亚洲av综合色区一区| 国产一区二区在线观看av| 日韩一区二区视频免费看| 人人妻人人添人人爽欧美一区卜| 肉色欧美久久久久久久蜜桃| 丰满饥渴人妻一区二区三| 性少妇av在线| 亚洲精华国产精华液的使用体验| 在线观看免费高清a一片| 纵有疾风起免费观看全集完整版| 如何舔出高潮| av又黄又爽大尺度在线免费看| 午夜福利,免费看| 大片免费播放器 马上看| 日韩人妻精品一区2区三区| 美女视频免费永久观看网站| av又黄又爽大尺度在线免费看| 一本一本久久a久久精品综合妖精| 亚洲国产最新在线播放| 男女边摸边吃奶| www.自偷自拍.com| 午夜久久久在线观看| 欧美激情极品国产一区二区三区| 美女国产高潮福利片在线看| 欧美日韩福利视频一区二区| 国产xxxxx性猛交| 日本av手机在线免费观看| 一级毛片黄色毛片免费观看视频| 一区二区三区激情视频| 亚洲一级一片aⅴ在线观看| 欧美激情 高清一区二区三区| 99九九在线精品视频| 电影成人av| 亚洲国产成人一精品久久久| 日本猛色少妇xxxxx猛交久久| 大香蕉久久网| 亚洲欧美一区二区三区国产| 精品一区二区三区四区五区乱码 | 久久久久久久久久久免费av| a级毛片黄视频| 日韩伦理黄色片| 午夜免费观看性视频| 中文字幕av电影在线播放| 在线看a的网站| 欧美日韩亚洲高清精品| 午夜久久久在线观看| 日韩熟女老妇一区二区性免费视频| 水蜜桃什么品种好| 80岁老熟妇乱子伦牲交| 日日爽夜夜爽网站| 日韩一本色道免费dvd| 狠狠精品人妻久久久久久综合| 国产精品免费大片| 90打野战视频偷拍视频| 久久人人爽人人片av| 国产97色在线日韩免费| 亚洲在久久综合| av不卡在线播放| 婷婷色综合大香蕉| 日韩电影二区| av卡一久久| 亚洲精品自拍成人| videos熟女内射| 成人18禁高潮啪啪吃奶动态图| 香蕉丝袜av| 美女视频免费永久观看网站| 日韩av免费高清视频| 欧美人与善性xxx| 成年人午夜在线观看视频| 美女主播在线视频| 精品人妻在线不人妻| 午夜91福利影院| 国产精品国产三级专区第一集| 午夜福利视频在线观看免费| 久久人人爽av亚洲精品天堂| 黄频高清免费视频| 亚洲国产欧美日韩在线播放| √禁漫天堂资源中文www| 老司机影院毛片| 在线观看www视频免费| 中文精品一卡2卡3卡4更新| 国产精品.久久久| 99热全是精品| 精品久久蜜臀av无| 免费在线观看完整版高清| 精品久久久精品久久久| 国产精品成人在线| 国产又爽黄色视频| 中文字幕高清在线视频| 国产精品偷伦视频观看了| 欧美97在线视频| 久久国产亚洲av麻豆专区| 日韩大码丰满熟妇| 国产免费又黄又爽又色| 毛片一级片免费看久久久久| 国产精品偷伦视频观看了| 美女中出高潮动态图| 午夜福利视频精品| 在线观看免费日韩欧美大片| 欧美97在线视频| 99久久99久久久精品蜜桃| 国产精品一国产av| 两个人免费观看高清视频| 最近手机中文字幕大全| 男人操女人黄网站| 赤兔流量卡办理| 亚洲四区av| 日本wwww免费看| 99国产精品免费福利视频| 超碰成人久久| 午夜老司机福利片| 男女床上黄色一级片免费看| 免费在线观看视频国产中文字幕亚洲 | 日韩中文字幕欧美一区二区 | 亚洲av电影在线观看一区二区三区| 中文字幕人妻丝袜制服| 女人被躁到高潮嗷嗷叫费观| av天堂久久9| 纵有疾风起免费观看全集完整版|