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

    遠海船舶噪聲分布初步研究

    2017-12-01 03:21:20薛睿超楊燕明黃二輝文洪濤
    聲學技術 2017年5期
    關鍵詞:船速航船環(huán)境噪聲

    薛睿超,楊燕明,黃二輝,文洪濤

    ?

    遠海船舶噪聲分布初步研究

    薛睿超,楊燕明,黃二輝,文洪濤

    (國家海洋局第三海洋研究所,福建廈門 361005)

    艦船數量、噸位以及航行距離的急劇上升導致了海洋環(huán)境噪聲級大幅度升高,對海洋環(huán)境的影響也日趨加重。但是當前的海洋環(huán)境噪聲研究所使用的歷史航船數據庫的信息卻非常匱乏且無法及時更新,不能反映出船舶噪聲真實的變化情況。此外,對于遠離岸邊的海域,岸基船舶自動識別系統(tǒng)(Automatic Identification System,AIS)與雷達無法監(jiān)測到,衛(wèi)星AIS技術的不成熟導致了船舶數據時效性和完整性的嚴重不足。因此,尚無可用的遠海大范圍的船舶數據。為了解決遠海船舶數據缺乏的問題,提出了一種基于多源AIS數據的數據融合方法,根據某海域在某個時間段內收到的所有AIS信息,處理后進行融合,融合結果可以得到任一時刻該海域的船舶分布情況,進而可計算出該海域的區(qū)域船舶噪聲級。

    海洋環(huán)境噪聲;船舶噪聲級;衛(wèi)星船舶自動識別系統(tǒng);遠海

    0 引言

    近年來,世界范圍內的大規(guī)模海上貿易日趨繁榮,大量的商船活動導致了海洋環(huán)境噪聲級的大幅升高。國內外對海洋環(huán)境噪聲的重視與日俱增[1],但是對海洋環(huán)境噪聲的認知卻停留在十幾年甚至幾十年前的測量數據上,海上聲學調查也只集中于軍事敏感區(qū)域,并且數據保密。人們對海洋環(huán)境噪聲實時變化的了解非常有限。有研究表明,1960~2000年,海洋環(huán)境噪聲呈上升趨勢,特別是低頻段的海洋環(huán)境噪聲以每10年3 dB的速度上升[2],這一趨勢直到近年來才逐漸減緩[3]。Ross[4-5]在1974年和1993年對上個世紀50~70年代的20年間,由于航船引起的低頻海洋環(huán)境噪聲的變化做了詳細的分析和總結。據Ross描述,到上個世紀70年代中期,海洋環(huán)境噪聲級的統(tǒng)計平均值已經明顯上升,不考慮氣候變化的影響,僅航船噪聲占重要地位的頻段(15~250 Hz)和地區(qū),平均噪聲級比Kunsen和Wenz時代增加了約3~5 dB,因為僅在1950年以后的25 年內,海上航船數量就增加了一倍多,并且船只的最大噸位和主機的推動力都大幅度增加。Piggott和Crouch[6-7]等人在Wenz測量的相同地點進行了海洋環(huán)境噪聲測量,結果發(fā)現,10~80 Hz范圍內海洋環(huán)境噪聲級增加了10 dB左右;高于80 Hz時頻段差異逐漸減小,在80~200 Hz之間僅增長了1~3 dB;高于200 Hz時海洋環(huán)境噪聲級變化不大。文獻[8]給出了1850~2000年間航船噪聲與自然噪聲的變化趨勢,如圖1所示??梢钥闯?,總的航船噪聲在上個世紀增長了接近40 dB。航船輻射噪聲是海洋環(huán)境噪聲低頻部分的主要來源,也是持續(xù)性噪聲的主要來源。在合適的傳播條件下,即使很遠處的航船輻射噪聲也能對聲接收點處的海洋環(huán)境噪聲強度產生影響。因此,掌握海面航船分布數據及其輻射噪聲源級對分析和掌握海洋環(huán)境噪聲級的變化趨勢十分重要。國外的航船數據已有多年積累,目前國內在此領域才剛剛起步,并缺乏有效的大規(guī)模測量手段[9]。

    圖1 航船噪聲和自然/生物噪聲發(fā)展趨勢

    船舶自動識別系統(tǒng)(Automatic Identification System,AIS),誕生于上世紀90年代初期,最初設計用于船舶避碰。AIS信息不僅包含了船舶時空上的信息,也給出了船舶當前的航行狀態(tài)。衛(wèi)星AIS的出現,對遠海船舶的航行具有革命性的意義[10-14],但是它的設計初衷并沒有考慮到從太空檢測信號的可能性,這導致了它在熱點海域的檢測存在盲點[15]。本文基于現有的AIS數據,提出了一種數據融合方法,可以在最大限度上還原遠海船舶的分布情況,根據AIS信息中的船速等信息可以計算出船舶的聲源級,進而計算出大范圍海域的航船噪聲分布。

    1 數據來源以及研究方法

    1.1 數據來源

    本文現有數據的主要來源為近年(2012~2015年)我國南海海域部分時段的航船數據,由國外商用衛(wèi)星AIS數據和部分岸基AIS數據組成。岸基AIS的覆蓋范圍僅為40 n mile(1 n mile=1852 m),因此無法滿足遠海用戶的使用需求。監(jiān)測大范圍海域的航船分布也首選衛(wèi)星AIS數據。中國南海是重要的海上運輸通道,航船密度高,目前也是海上爭端多發(fā)海域,事關我國海洋權益,戰(zhàn)略地位突出。如能有效掌握南海的航船分布,進而分析和掌握南海的航船輻射背景噪聲具有重要意義。因嚴重的時隙沖突,衛(wèi)星AIS在此類船舶高密度海區(qū)的船舶探測概率非常低[10-13],經過數據對比,發(fā)現很多區(qū)域的檢測概率低于20%,國外的商用衛(wèi)星AIS數據仍有大量遺漏和延時。

    1.2 研究方法

    為了解決商用衛(wèi)星AIS數據的遺漏和延時問題,本文考慮對其進行融合,具體方法和步驟如下:

    (1) 獲取一段連續(xù)時間的多源AIS數據,包括不同來源的星載AIS數據和岸基AIS數據(如有條件也可加入船載AIS數據);

    (2) 對所得到的數據進行統(tǒng)一的標準化處理和質量控制(因為數據來源不同,可能含有多種格式的數據且數據質量不一,需要統(tǒng)一數據格式以及剔除無效的數據);

    (3) 根據所得AIS數據中的時間、方位、船速、航向等信息,推算出每一時刻每一艘船舶的方位以及航行狀態(tài);

    (4) 通過自編的計算程序進行批量計算處理,繪制得到船舶分布圖和船舶軌跡圖。

    方法的實現流程如圖2所示。

    圖2 研究方法實現流程圖

    1.2.1 數據處理方法

    標準化后的數據為本文所需要的參數,包括水上移動通信業(yè)務標識碼(Maritime Mobile Service Identify,MMSI)、時間、船速、航向、經度、緯度等船舶相關參數;數據處理采用增加數據冗余的方法,即根據時間序列構造三維數據表,保證了數據處理的效率和準確度。本文中的時間分辨率設定為1 h,航行狀態(tài)包括船速、航向。另外,加入了時間間隔判斷條件,提高了結果的準確性;計算程序由Matlab語言編寫。

    1.2.2 計算模型

    本文提供了兩種船位計算模型,可以根據不同的需要進行選擇。

    模型一為線性模型。根據AIS數據中同一船舶不同時刻的經緯度信息,以固定的時間精度進行坐標位置的線性插值,結果如圖3所示。

    (a) 原始數據 (b) 計算結果

    圖3(a)中的兩個點表示同一船舶在不同時刻的船位,經過線性插值后得到圖3(b)的結果??梢杂檬?1)表示:

    模型一的計算方法簡單,速度快,但所得到的擬合結果精度較低,且當相鄰數據間隔時間太大會造成較大誤差。該模型適用于大范圍船舶分布的粗略統(tǒng)計以及其他對船位精度要求不高的應用。

    模型二為非線性(迭代)模型。根據AIS數據中同一船舶在不同時刻的經緯度、船速、航向等信息,進行綜合分析,采用迭代算法推算出每一時刻該船舶的船位以及相關航行狀態(tài)。示意圖見圖4。

    (a) 原始數據 (b) 計算結果

    2 船舶分布數據的計算分析

    本文的研究方法,利用AIS信息推算出不同時刻目標海域的船舶分布情況。首先,本文綜合了多源AIS數據,進行統(tǒng)一的標準化處理;其次,根據獲得的AIS信息,以時間為分組依據把所有參數分別進行分組;再根據方位、船速、航向信息進行擬合,構造連續(xù)時間的數組;最后,繪制船舶分布圖和軌跡圖。所獲取的原始數據的數量和質量越高,得到的結果越精確。下面對本文的研究方法進行說明。具體實施步驟如下:

    步驟1:根據獲得的岸基、衛(wèi)星數據進行標準化和質量控制,并提取所需參數,形成新的標準數據集。

    步驟2:處理數據。本文使用的多源星載AIS數據來源于星載AIS服務的供應商美國ORBCOMM公司和加拿大exactEarth公司的衛(wèi)星數據,岸基AIS數據來源于exactEarth公司,區(qū)域為中國南海部分海域的數據,數據時間包括2012年8月份和2015年5月份,共兩個月的數據(2012年數據不包含岸基數據)。

    步驟2-1:建立時間×船數×參數個數的三維數組,其中行為船舶的MMSI號,列為所獲數據的時間寬度(劃分到精度),頁為參數的個數,本文中選取的頁數為6,分別為時間、經度、緯度、船速、航向、衛(wèi)星數據標記;

    步驟2-2:把標準化后的數據按時間、MMSI號填入對應空格;

    步驟2-3:代入線性模型或者非線性迭代模型計算,所得結果填充表格數據,具體方法如下:

    (1) 線性模型:提取所有數據中的船舶的方位信息(經度、緯度),按照時間先后的順序進行排序,判斷相鄰數據點的時間間隔是否大于3天(可設置),若大于則在之間填充為“缺失”,符合小于條件的數據根據時間精度進行線性插值,得到中間所有時刻的船舶方位信息。

    (2) 非線性迭代模型:提取所有數據中的方位、航行狀態(tài)信息(包括經度、緯度、船速、航向),按照時間先后的順序進行排序,判斷相鄰數據點的時間間隔是否大于5天,若大于則在之間填充為“缺失”,符合小于條件的數據進行迭代計算:根據前后兩個數據的船速和航向數據,計算速度矢量在軸和軸方向的分量以及變化量,推算出中間時刻的方位、船速、航向信息。上述模型中加入時間間隔判斷條件是為了過濾部分中途停泊或者駛出區(qū)域范圍的船舶,以及惡劣天氣下(如臺風)錨泊的船舶。

    步驟2-4:去除異常點和超過最大允許時間間隔的數據;

    步驟2-5:存儲數據以便后續(xù)處理。

    步驟3:根據處理后的數據繪制相關示意圖。圖5(a)是由原始數據繪制得到的中國南海區(qū)域2012年8月份某時的船舶分布示意圖,不包括岸基數據,船舶數為264。圖5(b)為經本方法數據融合后得到的同一時刻同一區(qū)域的船舶分布示意圖,船舶數為1 629。圖5(c)是由原始數據繪制得到的南海區(qū)域2015年5月份某時的船舶分布示意圖,包括岸基數據,船舶數為2 092。圖5(d)為經本方法數據融合后得到的同一時刻同一區(qū)域的船舶分布示意圖,船舶數為7 887。通過對比可以看到,數據融合優(yōu)化前的船舶分布和數量嚴重失真,應用價值不大,優(yōu)化后的船舶分布和數量極大改善,接近真實情況。也可以看出,是否加入岸基數據,對結果中近岸船舶的分布影響很大。

    (a) 原始數據(2012年8月)

    (b) 融合后數據(2012年8月)

    (c) 原始數據(2015年5月)

    (d) 融合后數據(2015年5月)

    圖5 融合前后的船舶分布數據對比

    Fig.5 Ship data comparison before and after data fusion

    3 船舶噪聲的計算分析

    3.1 單一船舶聲源級計算[16-18]

    基于式(5)和融合后得到的南海航船分布數據,計算得到南海航船輻射噪聲的噪聲源級分布圖,如圖6所示。

    圖6 船舶聲源級分布

    3.2 區(qū)域船舶噪聲級計算

    考慮到某區(qū)域的海洋環(huán)境噪聲的影響因素有船舶、風、爆炸、打樁聲等等,各種因素疊加起來才得到該區(qū)域真實的海洋環(huán)境噪聲級。每一區(qū)域的聲場還受到該區(qū)域以及附近海域的海面海底反射以及介質散射的影響。如果考慮到這些,那么傳播模型需要精確的海洋環(huán)境參數數據,包括:海水聲速、海深、海底及海面邊界特性參數等等。計算相當復雜也缺少數據,因此本文只考慮船舶因素,并計算區(qū)域船舶噪聲級。

    區(qū)域船舶噪聲級是了解船舶噪聲空間分布以及其對海域內的海洋生物影響的重要參數。以往的區(qū)域船舶噪聲級計算方法使用的船舶數據主要來自歷史船舶數據庫(Historical Temporal Shipping,HITS)。它擁有全球最廣泛的船舶密度數據,能提供全球范圍內以1°*1°為單位的網格內船舶密度的數據,時間精度有月、季、年。其船舶類型包括商船、油輪、大型油輪、超大型油輪以及漁船。由于計算需要具體的船舶方位,所以根據數據庫的數據,得到相應海域的船舶密度,然后使用泊松分布函數得到離散的船舶數據。計算時還需要船舶的船速和長度數據,一般使用的是經驗數值,如表1所示。本文計算時若部分數據不完整,也采用表1的經驗數據。這種計算方法使用的船舶數據都是無法及時更新的以往數據和經驗數值,得到的船舶噪聲數據與真實測量值的誤差是不可避免的。

    表1 船長、船速的經驗數值

    為了計算區(qū)域的船舶噪聲源級,設計了如下計算模型:每一區(qū)域的聲源級取決于能對該區(qū)域產生影響的聲源,把這些聲源的貢獻疊加起來就能得到該區(qū)域的聲源級。該模型需要的數據有:不同的船舶序號,每艘船舶在不同時刻的坐標,每艘船舶的聲源級,每艘船舶在不同時刻的船速(船速不同,產生的聲源級也會不同)。具體處理步驟:

    (1) 把研究海域根據計算精度劃分成大小相同的網格(如1°*1°);

    (2) 計算海區(qū)內所有船舶的單獨聲源級;

    (3) 計算一定范圍內的船舶噪聲對每個網格的聲級貢獻,并去掉“奇點”。奇點指的是當前網格中心坐標1 km距離內的船舶,它會對計算結果產生非常大的影響。

    (4) 對結果圖像作局部平均以及平滑處理,平滑處理使用Matlab的shading interp函數。

    經過與美國HLS研究所的Michael Porter(聲學圖書館網站(http://oalib.hlsresearch.com/)的作者)的交流,計算方法的正確性得到了他的認同,并根據其意見在此基礎上作了一定改進(包括平滑處理等),結果如圖7所示。

    圖7 區(qū)域船舶噪聲級分布

    4 結論

    雖然從第二次世界大戰(zhàn)開始,軍事領域對于海洋低頻噪聲的研究和關注就沒有中斷過,然而所得到的有效數據還是非常匱乏。并且隨著海上貿易的發(fā)展,船舶數量、噸位和航行距離的增長,舊的數據庫已經不再適用。重新調查獲取數據所需的人力、物力成本對于任何國家、組織都是非常大的負擔,迫切需要一些較低成本的手段來獲取有效的數據。本文針對現有的衛(wèi)星AIS技術在熱點海域檢測概率不甚理想的問題,利用有限的AIS數據進行融合,融合后的數據能實時反映出海上任意范圍內船舶的位置分布和相關參數,提供了海上船舶分布的“快照”,且不需要安裝額外的設備。經過本文的方法優(yōu)化后,衛(wèi)星AIS檢測概率可以大幅度提高。本文把融合數據用于船舶噪聲預測,得到區(qū)域的船舶噪聲級分布圖,對于了解船舶噪聲時空分布有重要的參考價值。

    復雜的海洋環(huán)境中,需要綜合考慮的問題還有很多。研究結果目前存在不少缺陷:如只考慮了海區(qū)內船舶的噪聲,對高頻部分的噪聲還缺少相關的研究;沒有對淺、深海等不同的海洋條件進行區(qū)分;對于沒有安裝AIS設備的船舶缺少其數據;沒有考慮環(huán)境信息,特別是底質的聲學特性。所以本文給出的區(qū)域船舶噪聲級分布圖只是初步的結果。然而,我們相信,這些分布圖可以在說明船舶噪聲時提供一個重要的參考,以便將來做進一步的研究。

    [1] 郭新毅.海洋環(huán)境噪聲研究發(fā)展概述及應用前景[J]. 物理, 2014, 11(43): 723-731. GUO Xinyi. Overview of ocean ambient noise and application prospects[J]. Physics, 2014, 11(43): 723-731.

    [2] Mcdonald M A, Hildebrand J A, Wiggins S M. Increases in deep ocean ambient noise in the Northeast Pacific west of San Nicolas Island, California[J]. J. Acoust. Soc. Am., 2006, 120(2): 711-718.

    [3] Andrew R K, Howe B M, Mercer J A. Long-time trends in ship traffic noise for four sites off the north american west coast[J]. J. Acoust. Soc. Am., 2011, 129(2): 642-51.

    [4] 羅斯D [美]. 水下噪聲原理[M]. 北京: 海洋出版社, 1983. Ross D. Principle of underwater noise[M]. Beijing: China Ocean Press, 1983.

    [5] Ross D. On ocean underwater ambient noise[R]. Acoustics Bulletin January, 1993.

    [6] C L Piggott. Ambient sea noise at low frequencies in shallow water of the scotian shelf[J]. J. Acoust. Soc. Am., 1964, 36(11): 2152- 2163.

    [7] W W crouch, P J Burt. The logarithmic dependence of surface-generated ambient-sea-noise spectrum level on wind speed[J]. J. Acoust. Soc. Am., 1972, 51(3): 1066-1072.

    [8] 石杰, 張效民, 侯鐵雙, 等. 航船噪聲引起的低頻海洋環(huán)境噪聲級發(fā)展趨勢[J].魚雷技術, 2010, 18(2): 112-116. SHI Jie, ZHANG Xiaomin, HOU Tieshuang, et al. Noise level evolution of low-frequency ocean ambient noise caused by ship sources[J]. Torpedo Technology, 2010, 18(2): 112-116.

    [9] 衣雪娟, 林建恒, 孫軍平, 等. 海上航船分布及其對海洋環(huán)境噪聲的影響[J]. 海洋與湖沼, 2015, 46(6): 1270-1278. YI Xuejuan, LIN Jianheng, SUN Junping, et al. Ship distrbution in open sea and its impact on ocean ambient noise[J]. Oceanologia et Limnologia Sinica, 2015, 46(6): 1270-1278.

    [10] 程云, 陳利虎, 陳小前. 星載AIS檢測概率建模與仿真分析[J]. 國防科技大學學報, 2014, 36(3): 51-57. CHENG Yun, CHEN Lihu, CHEN Xiaoqian. Modeling and simulation analysis of detection probability for space-based AIS[J]. Journal of National University of Defense Technology, 2014, 36(3): 51-57.

    [11] 張昕, 李洪星, 陳萍, 等. 基于實際船舶分布的星載AIS仿真系統(tǒng)設計[J]. 無線電工程, 2014, 44(10): 15-17. ZHANG Xin, LI Hongxing, CHEN Ping, et al. Design of space- based AIS simulation system based on actual ship distribution[J]. Radio Engineering, 2014, 44(10): 15-17.

    [12] 宋果林, 鄒光南, 王艷峰, 等. 衛(wèi)星自動識別系統(tǒng)觀測模型及檢測概率分析[J]. 航天器工程, 2015, 24(4): 136-142. SONG Guolin, ZOU Guangnan, WANG Yanfeng, et al. Satellite-borne AIS observe model and detection probability analyze[J]. Spacecraft Engineering, 2015, 24(4): 136-142.

    [13] 劉暢. 船舶自動識別系統(tǒng)AIS關鍵技術研究[D]. 大連: 大連海事大學, 2013. LIU Chang. Study of key tennology of automatic identification system AIS[D]. Dalian: Dalian Maritime University, 2013.

    [14] 張哲, 馬桂山.衛(wèi)星探測AIS的分析及其發(fā)展的研究[J]. 航海技術, 2010(3): 38-40. ZHANG Zhe, MA Guishan. Analysis of satellite detection of AIS and study on its development[J]. Marine Technology, 2010(3): 38-40.

    [15] Yang M, Zou Y, Fang L. Collision and detection performance with three overlap signal collisions in space-based AIS reception[C]//Trust, Security and Privacy in Computing and Communications (TrustCom), 2012 IEEE 11th International Conference on. IEEE, 2012: 1641-1648.

    [16] Hamson R M. The modelling of ambient noise due to shipping and wind sources in complex environments[J]. Applied Acoustics, 1997, 51(3): 251-287.

    [17] 汪德昭, 尚爾昌. 水聲學M]. 2版. 北京: 科學出版社, 2013. WANG Dezhao, SHANG Erchang. Underwater acoustics[M]. 2ndedition. Beijing: Science Press, 2103.

    [18] 劉伯勝, 雷家煜. 水聲學原理[M]. 哈爾濱: 哈爾濱工程大學出版社, 2009: 235-248. LIU Bosheng, LEI Jiayi. Underwater acoustics principle[M]. Harbin: Harbin Engineering University Press, 2009: 235-248.

    Research on pelagic ship noise distribution

    XUE Rui-chao, YANG Yan-ming, HUANG Er-hui, WEN Hong-tao

    (Third Institute of Oceanography, State Oceanic Administration, Xiamen 361005, Fujian, China)

    As the number of ships continues to grow, the quantities of goods carried and the distances of ship traveled have made the sea ambient noise level dramatically rising, which impacts on the marine environment increasingly worse. Moreover, the ship history database can not reflect the real situation and is no longer applicable. Especially for the pelagic ships, there is a shortage of data available. With the aim of palliating this situation, this paper puts forward a data fusion method of multi-source AIS, which can get the every hour distribution of ships and then calculate the regional ship noise level for sound mapping.

    sea ambient noise; noise levels of ships; satellite-based AIS (Automatic Identification System); open sea

    OP733.22

    A

    1000-3630(2017)-05-0467-06

    10.16300/j.cnki.1000-3630.2017.05.012

    2017-01-13;

    2017-05-09

    國家海洋局第三海洋研究所基本科研業(yè)務專項(海三科2015016)、福建省科技計劃項目(2016H0018)

    薛睿超(1988-), 男, 福建惠安人, 碩士, 研究方向為船舶噪聲, AIS數據分析。

    楊燕明, E-mail: yangyanming@tio.org.cn

    猜你喜歡
    船速航船環(huán)境噪聲
    南湖航船指方向
    能效管理中的船速優(yōu)化
    環(huán)境噪聲智能檢測終端設計
    電子制作(2018年23期)2018-12-26 01:01:04
    基于云服務的環(huán)境噪聲智能監(jiān)測系統(tǒng)研究
    測控技術(2018年11期)2018-12-07 05:49:22
    淺談在強風條件下操縱大型LNG船靠泊天津臨港
    珠江水運(2018年5期)2018-04-12 16:48:04
    南海某海域環(huán)境噪聲與風速相關性研究*
    中國港航船企指數(CMEI)
    鐵路環(huán)境噪聲測量中的問題分析
    重載CAPESIZE船舶乘潮進連云港泊位實踐
    中國水運(2015年10期)2015-11-09 15:01:47
    2014中國港航船企收入榜單
    一本一本综合久久| 蜜桃在线观看..| 欧美日本中文国产一区发布| 五月伊人婷婷丁香| 综合色丁香网| 天堂俺去俺来也www色官网| 美女中出高潮动态图| 国产黄片美女视频| 亚洲国产精品专区欧美| 国语对白做爰xxxⅹ性视频网站| 亚洲国产欧美日韩在线播放 | 国产亚洲最大av| 国产在视频线精品| 一级爰片在线观看| 免费av不卡在线播放| 日本午夜av视频| 免费不卡的大黄色大毛片视频在线观看| 在线观看免费高清a一片| 免费在线观看成人毛片| 亚洲性久久影院| 美女大奶头黄色视频| 亚洲自偷自拍三级| 久久久久网色| 国产精品人妻久久久久久| 制服丝袜香蕉在线| 免费久久久久久久精品成人欧美视频 | 全区人妻精品视频| 大码成人一级视频| 精品一区二区三区视频在线| 一级毛片aaaaaa免费看小| 你懂的网址亚洲精品在线观看| 精品国产一区二区久久| 精品国产国语对白av| 在线亚洲精品国产二区图片欧美 | 久久国产乱子免费精品| 在线观看美女被高潮喷水网站| 免费在线观看成人毛片| 国产色婷婷99| 你懂的网址亚洲精品在线观看| 久久久久久久久大av| av卡一久久| 国产爽快片一区二区三区| 国产在线视频一区二区| 亚洲怡红院男人天堂| 毛片一级片免费看久久久久| 久久精品国产a三级三级三级| a级片在线免费高清观看视频| 日韩人妻高清精品专区| 丰满乱子伦码专区| 午夜激情福利司机影院| 99国产精品免费福利视频| 国产精品福利在线免费观看| 国精品久久久久久国模美| 国产一级毛片在线| 人妻制服诱惑在线中文字幕| 极品人妻少妇av视频| 国模一区二区三区四区视频| 免费少妇av软件| 国产免费福利视频在线观看| 亚洲精品国产成人久久av| 欧美精品亚洲一区二区| h视频一区二区三区| 久久久久久久久久久久大奶| a级毛片免费高清观看在线播放| 国产精品99久久99久久久不卡 | 婷婷色麻豆天堂久久| 国国产精品蜜臀av免费| 久久精品久久久久久久性| 在线播放无遮挡| 亚洲精品自拍成人| 校园人妻丝袜中文字幕| av黄色大香蕉| 中文天堂在线官网| 欧美另类一区| 91精品国产九色| 国产一区二区三区综合在线观看 | 免费人妻精品一区二区三区视频| 中国国产av一级| 欧美3d第一页| 免费av不卡在线播放| 久久久久久伊人网av| 一级毛片我不卡| 91成人精品电影| 午夜福利在线观看免费完整高清在| 亚洲怡红院男人天堂| 插阴视频在线观看视频| 男人添女人高潮全过程视频| 一级,二级,三级黄色视频| 97在线视频观看| 99热6这里只有精品| 国产成人一区二区在线| 九九久久精品国产亚洲av麻豆| 大片免费播放器 马上看| 黑人巨大精品欧美一区二区蜜桃 | 欧美人与善性xxx| 国产午夜精品久久久久久一区二区三区| 色视频在线一区二区三区| 国产精品伦人一区二区| 99视频精品全部免费 在线| 日韩中文字幕视频在线看片| 一级黄片播放器| 人人妻人人爽人人添夜夜欢视频 | 国产乱来视频区| 高清视频免费观看一区二区| 少妇人妻 视频| av国产久精品久网站免费入址| 日日啪夜夜爽| 你懂的网址亚洲精品在线观看| 精品午夜福利在线看| 下体分泌物呈黄色| 久久av网站| 一区在线观看完整版| 亚洲成人手机| 日韩欧美一区视频在线观看 | 国产精品蜜桃在线观看| 国产av精品麻豆| 五月玫瑰六月丁香| a级毛色黄片| 黄色配什么色好看| 色婷婷av一区二区三区视频| 国产精品一区www在线观看| 欧美高清成人免费视频www| 街头女战士在线观看网站| 啦啦啦视频在线资源免费观看| 老熟女久久久| 男女边摸边吃奶| 18禁在线播放成人免费| 国产亚洲av片在线观看秒播厂| 人妻夜夜爽99麻豆av| 丰满乱子伦码专区| 亚洲人成网站在线播| 欧美区成人在线视频| 欧美精品国产亚洲| 亚洲国产欧美日韩在线播放 | 51国产日韩欧美| 亚洲国产精品999| 欧美日韩精品成人综合77777| 亚洲成人手机| 国产午夜精品一二区理论片| 精品国产一区二区久久| 少妇熟女欧美另类| 99视频精品全部免费 在线| 日韩av免费高清视频| 国产亚洲一区二区精品| 插逼视频在线观看| 成人特级av手机在线观看| 成人二区视频| 成人毛片60女人毛片免费| 青春草亚洲视频在线观看| 纯流量卡能插随身wifi吗| 欧美日韩综合久久久久久| 2018国产大陆天天弄谢| 91精品国产九色| 久久国内精品自在自线图片| 在线播放无遮挡| 国内揄拍国产精品人妻在线| 国产 精品1| 亚洲久久久国产精品| 九九久久精品国产亚洲av麻豆| 麻豆成人午夜福利视频| 久久久久国产网址| 国产精品人妻久久久影院| 国产成人aa在线观看| 久久久久久久久久成人| 一区二区av电影网| 嫩草影院入口| 免费黄频网站在线观看国产| 成人毛片a级毛片在线播放| 精品久久久久久电影网| 日韩一区二区三区影片| 免费观看av网站的网址| 人人澡人人妻人| 亚洲av在线观看美女高潮| 国产精品一区二区性色av| 少妇的逼好多水| 超碰97精品在线观看| www.av在线官网国产| 国产精品无大码| 亚洲内射少妇av| 十分钟在线观看高清视频www | 少妇的逼好多水| 久久精品久久精品一区二区三区| 国产一级毛片在线| 成人亚洲精品一区在线观看| 3wmmmm亚洲av在线观看| 啦啦啦中文免费视频观看日本| 97精品久久久久久久久久精品| 精品视频人人做人人爽| 久久久久久久久久久免费av| 亚洲国产av新网站| 日韩成人av中文字幕在线观看| 水蜜桃什么品种好| 男人添女人高潮全过程视频| 欧美激情国产日韩精品一区| 在线观看av片永久免费下载| 国产欧美另类精品又又久久亚洲欧美| 一区二区三区免费毛片| 狂野欧美白嫩少妇大欣赏| 纯流量卡能插随身wifi吗| 日韩在线高清观看一区二区三区| 三级国产精品欧美在线观看| 国产精品国产三级专区第一集| 成人影院久久| 高清不卡的av网站| 女性被躁到高潮视频| 精品亚洲成国产av| 亚洲国产精品999| 插逼视频在线观看| 爱豆传媒免费全集在线观看| 久久国产亚洲av麻豆专区| 婷婷色av中文字幕| 国产精品99久久久久久久久| 男人舔奶头视频| 免费久久久久久久精品成人欧美视频 | 99久久精品国产国产毛片| 国产日韩欧美亚洲二区| 国产欧美亚洲国产| 国产午夜精品久久久久久一区二区三区| 有码 亚洲区| 国内精品宾馆在线| 97在线人人人人妻| 岛国毛片在线播放| 哪个播放器可以免费观看大片| 成人漫画全彩无遮挡| 丰满饥渴人妻一区二区三| 亚洲欧美成人精品一区二区| 在线观看av片永久免费下载| 少妇的逼好多水| 欧美日韩av久久| 欧美变态另类bdsm刘玥| 爱豆传媒免费全集在线观看| 自线自在国产av| 亚洲精品视频女| 亚洲高清免费不卡视频| 内射极品少妇av片p| 一级片'在线观看视频| 日韩一区二区三区影片| 在线观看av片永久免费下载| 一本一本综合久久| 熟女电影av网| 91精品国产国语对白视频| av视频免费观看在线观看| 99热网站在线观看| 中文字幕亚洲精品专区| av播播在线观看一区| 日韩亚洲欧美综合| 99视频精品全部免费 在线| 国产综合精华液| 日本91视频免费播放| 久久国产乱子免费精品| 下体分泌物呈黄色| av福利片在线观看| 一区二区三区四区激情视频| 少妇的逼水好多| 菩萨蛮人人尽说江南好唐韦庄| 麻豆精品久久久久久蜜桃| av.在线天堂| 中文字幕av电影在线播放| 午夜免费观看性视频| 亚洲av.av天堂| 久久精品久久久久久久性| 成人毛片60女人毛片免费| 大片免费播放器 马上看| 国产毛片在线视频| 久久午夜综合久久蜜桃| 赤兔流量卡办理| 免费不卡的大黄色大毛片视频在线观看| 如何舔出高潮| 最新的欧美精品一区二区| 亚洲av免费高清在线观看| 精品酒店卫生间| 成人特级av手机在线观看| 少妇精品久久久久久久| 99热网站在线观看| 乱人伦中国视频| 精品国产乱码久久久久久小说| 亚洲天堂av无毛| 精品一区二区三区视频在线| 日本-黄色视频高清免费观看| 亚洲,一卡二卡三卡| 中文乱码字字幕精品一区二区三区| 久久久欧美国产精品| 18禁裸乳无遮挡动漫免费视频| 精品99又大又爽又粗少妇毛片| 久久ye,这里只有精品| 婷婷色综合大香蕉| 日韩亚洲欧美综合| av.在线天堂| 国产精品久久久久久av不卡| www.av在线官网国产| 我要看日韩黄色一级片| 国产午夜精品一二区理论片| 精品午夜福利在线看| 久久久久久久久久久丰满| 日本色播在线视频| 另类精品久久| a级毛片在线看网站| 国产亚洲精品久久久com| 色94色欧美一区二区| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看 | 精品一品国产午夜福利视频| 久久久国产一区二区| 精品人妻一区二区三区麻豆| 国产欧美日韩综合在线一区二区 | 久久久亚洲精品成人影院| 亚洲四区av| 欧美最新免费一区二区三区| 18禁在线无遮挡免费观看视频| 99国产精品免费福利视频| 狂野欧美激情性xxxx在线观看| 亚洲第一av免费看| 爱豆传媒免费全集在线观看| 99热这里只有是精品50| 18禁裸乳无遮挡动漫免费视频| 亚洲不卡免费看| 观看av在线不卡| 99热全是精品| 亚洲真实伦在线观看| 亚洲国产精品成人久久小说| 亚洲av日韩在线播放| 午夜老司机福利剧场| 黄色一级大片看看| 少妇裸体淫交视频免费看高清| av视频免费观看在线观看| 国产亚洲欧美精品永久| 午夜影院在线不卡| 啦啦啦在线观看免费高清www| 91成人精品电影| 国产伦精品一区二区三区视频9| 极品教师在线视频| 免费观看a级毛片全部| 日韩一区二区三区影片| av播播在线观看一区| 久久精品国产亚洲网站| av在线观看视频网站免费| 在线免费观看不下载黄p国产| 91久久精品电影网| 久久久久国产精品人妻一区二区| 男女啪啪激烈高潮av片| 一级黄片播放器| 日韩av不卡免费在线播放| 老司机影院成人| 日本免费在线观看一区| 热99国产精品久久久久久7| 国产av码专区亚洲av| 亚洲婷婷狠狠爱综合网| 人体艺术视频欧美日本| 亚洲精品国产av成人精品| 欧美日韩综合久久久久久| 久久午夜福利片| 性高湖久久久久久久久免费观看| 久久国产亚洲av麻豆专区| 一个人看视频在线观看www免费| 18禁在线播放成人免费| 日日摸夜夜添夜夜添av毛片| 天堂8中文在线网| 91精品一卡2卡3卡4卡| 黑人猛操日本美女一级片| 国产伦精品一区二区三区四那| 亚洲国产欧美在线一区| 国产欧美日韩精品一区二区| 多毛熟女@视频| 久热这里只有精品99| 2018国产大陆天天弄谢| 少妇的逼水好多| 91久久精品国产一区二区三区| 日韩欧美一区视频在线观看 | 国产精品一区二区三区四区免费观看| 久久 成人 亚洲| 成年人免费黄色播放视频 | 日韩在线高清观看一区二区三区| 9色porny在线观看| 国产一区二区三区综合在线观看 | 国产高清国产精品国产三级| 久久亚洲国产成人精品v| 黑人猛操日本美女一级片| 少妇被粗大的猛进出69影院 | 国产精品偷伦视频观看了| 免费在线观看成人毛片| 国产 一区精品| 晚上一个人看的免费电影| 精品国产乱码久久久久久小说| 在线亚洲精品国产二区图片欧美 | 一区二区三区免费毛片| 久久人妻熟女aⅴ| 国产精品嫩草影院av在线观看| 国产免费福利视频在线观看| 最近2019中文字幕mv第一页| 2022亚洲国产成人精品| 亚洲精品中文字幕在线视频 | 在线免费观看不下载黄p国产| 亚洲av中文av极速乱| 国产黄色视频一区二区在线观看| 亚洲av日韩在线播放| 久久久久精品性色| 国产亚洲5aaaaa淫片| 另类精品久久| 3wmmmm亚洲av在线观看| 97在线人人人人妻| 一级片'在线观看视频| 日韩欧美一区视频在线观看 | 黄色日韩在线| 日韩免费高清中文字幕av| 欧美xxⅹ黑人| 人妻少妇偷人精品九色| 日产精品乱码卡一卡2卡三| 麻豆精品久久久久久蜜桃| 如日韩欧美国产精品一区二区三区 | 成人无遮挡网站| 免费观看在线日韩| 狂野欧美白嫩少妇大欣赏| av在线观看视频网站免费| 国产精品久久久久久精品古装| 精品久久久久久久久亚洲| 高清毛片免费看| 日本wwww免费看| 乱码一卡2卡4卡精品| 18禁裸乳无遮挡动漫免费视频| 夜夜爽夜夜爽视频| 久久午夜综合久久蜜桃| 大香蕉久久网| 一本久久精品| 交换朋友夫妻互换小说| 日本黄色片子视频| 少妇丰满av| 97在线人人人人妻| 精品久久久久久久久av| 中文资源天堂在线| 亚洲,一卡二卡三卡| .国产精品久久| 大香蕉97超碰在线| 国产黄色免费在线视频| 国产精品伦人一区二区| 美女国产视频在线观看| 在线观看免费视频网站a站| 久久午夜综合久久蜜桃| 女人久久www免费人成看片| 亚洲欧美日韩卡通动漫| 欧美区成人在线视频| 久久久久网色| 中文欧美无线码| 欧美成人午夜免费资源| 中文精品一卡2卡3卡4更新| 精品视频人人做人人爽| 国产综合精华液| 久久久久久久久久久免费av| 国产精品熟女久久久久浪| 亚洲精品亚洲一区二区| 亚洲国产成人一精品久久久| 成年av动漫网址| 精品久久久久久久久亚洲| 97在线人人人人妻| 午夜久久久在线观看| 国产精品国产av在线观看| 亚洲精品久久午夜乱码| 国产高清不卡午夜福利| 国产 精品1| 国产午夜精品久久久久久一区二区三区| 嘟嘟电影网在线观看| 国产精品秋霞免费鲁丝片| 一个人免费看片子| 国产亚洲最大av| 免费看光身美女| 久久99热这里只频精品6学生| 老司机影院毛片| 曰老女人黄片| 国产午夜精品久久久久久一区二区三区| 自线自在国产av| 亚洲av福利一区| 中文资源天堂在线| 王馨瑶露胸无遮挡在线观看| 我的老师免费观看完整版| 秋霞在线观看毛片| 亚洲美女视频黄频| 大又大粗又爽又黄少妇毛片口| 日韩强制内射视频| av有码第一页| 亚洲国产色片| 精华霜和精华液先用哪个| 亚洲图色成人| 777米奇影视久久| 亚洲无线观看免费| 久久 成人 亚洲| 寂寞人妻少妇视频99o| 三级国产精品片| 又大又黄又爽视频免费| 精品视频人人做人人爽| 久久精品久久精品一区二区三区| 亚洲第一区二区三区不卡| 麻豆乱淫一区二区| 一区二区av电影网| a级一级毛片免费在线观看| 欧美性感艳星| 色网站视频免费| 欧美xxxx性猛交bbbb| 精品国产乱码久久久久久小说| 欧美精品高潮呻吟av久久| 国产亚洲最大av| 不卡视频在线观看欧美| 亚洲av免费高清在线观看| 自拍偷自拍亚洲精品老妇| 高清毛片免费看| 有码 亚洲区| 久久久久国产网址| 成年人午夜在线观看视频| 国产精品无大码| 精品亚洲成国产av| 成人漫画全彩无遮挡| av天堂中文字幕网| 国产黄色视频一区二区在线观看| 国内揄拍国产精品人妻在线| 日韩欧美一区视频在线观看 | 亚洲国产精品一区二区三区在线| 色视频www国产| av免费在线看不卡| 久久av网站| 亚洲精品国产av成人精品| 免费高清在线观看视频在线观看| 精品人妻一区二区三区麻豆| 我要看日韩黄色一级片| 国产av精品麻豆| 久久人人爽人人片av| 九九久久精品国产亚洲av麻豆| 80岁老熟妇乱子伦牲交| 国产精品久久久久久av不卡| 男女边摸边吃奶| 亚洲av在线观看美女高潮| a级毛片免费高清观看在线播放| 国产黄片美女视频| 水蜜桃什么品种好| 国产精品秋霞免费鲁丝片| 国产在线男女| 大片免费播放器 马上看| 男女啪啪激烈高潮av片| 乱人伦中国视频| √禁漫天堂资源中文www| 国产一区二区三区综合在线观看 | 能在线免费看毛片的网站| 2022亚洲国产成人精品| 久久久国产欧美日韩av| 色视频在线一区二区三区| 亚洲欧美一区二区三区黑人 | 久久久久视频综合| 亚洲,欧美,日韩| 高清欧美精品videossex| 国产一区二区三区av在线| 免费少妇av软件| 久久97久久精品| 亚洲av欧美aⅴ国产| 亚洲美女搞黄在线观看| 黑人猛操日本美女一级片| 免费av不卡在线播放| 99热全是精品| 国产精品国产三级国产av玫瑰| 久久久久久久久久成人| 亚洲精品日本国产第一区| 欧美bdsm另类| 80岁老熟妇乱子伦牲交| 日韩欧美精品免费久久| 久久99精品国语久久久| 亚洲av在线观看美女高潮| 中国国产av一级| 精品熟女少妇av免费看| 国产无遮挡羞羞视频在线观看| 777米奇影视久久| 国产在线免费精品| 欧美日韩av久久| 日韩制服骚丝袜av| 国产毛片在线视频| 成年人午夜在线观看视频| 久久久久精品久久久久真实原创| 另类亚洲欧美激情| 国语对白做爰xxxⅹ性视频网站| 国产午夜精品久久久久久一区二区三区| 观看免费一级毛片| 久久久久久久精品精品| 两个人的视频大全免费| 久久精品国产亚洲av天美| 国产日韩一区二区三区精品不卡 | 最新的欧美精品一区二区| 一个人免费看片子| 久久6这里有精品| 人妻夜夜爽99麻豆av| 亚洲美女搞黄在线观看| 亚洲av电影在线观看一区二区三区| 成年美女黄网站色视频大全免费 | 99re6热这里在线精品视频| 春色校园在线视频观看| 街头女战士在线观看网站| 欧美97在线视频| av在线播放精品| 在线观看av片永久免费下载| av线在线观看网站| 亚洲av电影在线观看一区二区三区| 高清午夜精品一区二区三区| 久久99热6这里只有精品| 一级片'在线观看视频| 亚洲中文av在线| a级毛色黄片| 国产成人91sexporn| 欧美老熟妇乱子伦牲交| 日韩亚洲欧美综合| 精品久久久久久电影网| 内射极品少妇av片p| 亚洲人与动物交配视频| 一本色道久久久久久精品综合| 国产欧美亚洲国产| 少妇熟女欧美另类| 中文精品一卡2卡3卡4更新| 免费av不卡在线播放| 国产精品久久久久久精品古装| 国国产精品蜜臀av免费| 免费观看av网站的网址| 内地一区二区视频在线|