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

    基于貝葉斯時空模型的洋山港及其附近水域事故險情分布特征分析

    2024-06-05 13:25:25劉劍平杜佳欣尚麟
    航海 2024年3期
    關鍵詞:洋山港風險

    劉劍平 杜佳欣 尚麟

    摘 要:洋山深水港自2005年建港以來,已迅速崛起為中國最重要的國際集裝箱樞紐港。近期規(guī)劃調(diào)整將在北作業(yè)區(qū)設立新航道并新增錨地,這將對洋山深水港及附近水域的通航環(huán)境帶來實質(zhì)性改變。本文采用貝葉斯時空模型,考慮時空交互效應,對事故險情進行分析,提前研判水域風險。研究結(jié)果顯示,小洋山港區(qū)和東海大橋附近水域的風險最高,航行安全水平不穩(wěn)定。大戟洋和北鼎星附近水域風險次高,蘆潮港沿岸碼頭水域風險最低。事故險情主要發(fā)生在春夏兩季,特別是3月、8月和10月。這有助于制定更具針對性的監(jiān)管政策和安全措施,確保航行安全和港口高效運營。

    關鍵詞:貝葉斯時空模型;事故險情;風險;洋山港

    0 引 言

    上海港地處長江下游入海口、東海之濱,是世界上最繁忙的港口之一[1]。自然條件復雜,船舶種類多樣,船舶流量龐大。自2001年明確建設上海國際航運中心目標以來,上海港集裝箱吞吐量于2010年首次超越新加坡港成為世界第一,之后持續(xù)保持領先地位[2]。洋山深水港作為中國首個在海島建設的港口,位于杭州灣口[3],自2005年開港以來,其集裝箱吞吐量逐年增長,在上海港連續(xù)13年的集裝箱吞吐量排名第一中發(fā)揮著重要作用[4]。

    2016年12月,上海市人民政府與浙江省人民政府簽署了《關于共同推進小洋山區(qū)域開發(fā)等重大合作事項的框架協(xié)議》,旨在共同推進小洋山進一步綜合開發(fā)。新的小洋山北作業(yè)區(qū)項目的啟動將引起水域通航環(huán)境的實質(zhì)性變化,使得警戒區(qū)和附近水域的船舶交通管理變得更加復雜。因此,需要提前進行針對性研究,識別水域船舶航行風險,厘清水域事故險情分布特征,保障洋山深水港區(qū)及其附近水域的水上交通安全,確保航行安全和港口運營的可持續(xù)性[5]。

    貝葉斯時空模型是一種強大的統(tǒng)計學習方法,它綜合利用先驗信息和觀測數(shù)據(jù),以概率分布表示參數(shù)不確定性,從而實現(xiàn)更準確、可靠的建模,在探究時空變化或時空分布特征方面具有廣泛應用。鄒齊[6]利用貝葉斯時空分層模型,對新冠疫情傳播規(guī)律的時空相依性進行研究,為有針對性地進行疫情防控提供參考。董文錢[7]等人基于貝葉斯時空模型,探究街面秩序、市容環(huán)境等城管事件的時空演變特征及影響因素,為城市管理提供支持。徐冰[8]針對湖南省手足口病傳染問題,構建貝葉斯時空模型,研究手足口病的高發(fā)季節(jié)與時間,為建立傳染病早期預警系統(tǒng)奠定基礎。Yu[9]等人基于華北、朝鮮半島和日本的部分地區(qū)10年的碰撞報告,考慮空間和時間的相互影響,建立了貝葉斯時空模型來評估碰撞風險。

    貝葉斯時空模型通過分析時空分布特征,為風險防控提供了可靠的工具[10],使得有關部門能夠更準確地理解和預測潛在風險,從而采取針對性的措施。因此,本文利用貝葉斯時空模型,對洋山深水港及其附近水域事故險情分布特征進行建模分析,研究水域風險水平分布特點,為船舶交通安全組織管理提供參考。

    1 研究區(qū)域概況與數(shù)據(jù)來源

    1.1 研究區(qū)域

    本文研究區(qū)域為洋山海事局VTS管理服務區(qū),洋山海事局管轄洋山深水港及其附近超4 000 m2水域,VTS管理服務區(qū)示意圖如圖1所示[11]。岸線達123 km,轄區(qū)呈現(xiàn)“六多一線”的顯著特點,即:航道多,重要水工構筑物多;大型集裝箱船進出港艘次多,集裝箱吞吐量大;客運航線多,客船多;危險品船多;重點工程多;惡劣天氣多;東海大橋生命線。該水域船舶密度大、交通流復雜、種類多樣化、船舶大型化特征明顯,宏觀的東西向進出港與南北向過境船舶之間船長、吃水、航速差異較為明顯,總體交通流維度的交通環(huán)境較為復雜,安全管理難度較大。

    1.2事故險情報告數(shù)據(jù)及預處理過程

    為掌握水域事故險情發(fā)生特點,加強重點水域管理,對洋山海事局VTS管理服務區(qū)2017年1月至2023年4月事故險情報告進行分析,共收集134份有效記錄,在建模之前,需要對事故險情數(shù)據(jù)進行收集與預處理,事故險情數(shù)據(jù)的收集與預處理是進行基于貝葉斯時空模型的事故險情分布特征分析的關鍵步驟,準確、完整的數(shù)據(jù)收集以及有效的預處理過程對于后續(xù)模型建立和結(jié)果分析至關重要。

    根據(jù)洋山VTS管理服務區(qū)部分航路的位置,將該區(qū)域劃分為金山航道東西向西段至金山沿岸碼頭、蘆潮港沿岸碼頭、小洋山港區(qū)及東海大橋、大戟洋及北鼎星、主航道警戒區(qū)及大洋山以南5個區(qū)域,各區(qū)域的劃分及位置如圖2所示:

    本文提取事故險情報告中的時間與位置信息,根據(jù)經(jīng)緯度位置信息,將事故對應到各個區(qū)域中,對位置信息進行處理。事故險情發(fā)生位置屬于區(qū)域1時,位置點記為“1”;事故險情發(fā)生位置屬于區(qū)域2時,位置點記為“2”。以此類推,將所有的位置信息進行歸類,并記錄事故險情的發(fā)生年份,組成一組具有時空信息的數(shù)據(jù)樣本,并收集整理各個區(qū)域的交通流量,作為模型的輸入。

    2 研究方法

    根據(jù)貝葉斯理論,貝葉斯模型可以將樣本信息與先驗分布相結(jié)合來估計后驗分布,后驗分布表示影響因素導致事故險情發(fā)生的概率。由于典型的貝葉斯模型遵循泊松分布[12],因此似然函數(shù)可以表示為:

    L(y|μ)=i=1∏5? ?μiyiexp(-μi)/yi?。?)

    式中,yi表示服從泊松分布的船舶事故險情發(fā)生次數(shù),μi表示事故險情發(fā)生次數(shù)的期望值,i表示事故險情發(fā)生區(qū)域,i∈{1,2,3,4,5}。

    假定船舶事故險情發(fā)生是相互獨立的。因此,泊松分布μi的期望值可建模為 :

    E(yi)=μi=eiθi(2)

    式中ei為事故險情發(fā)生的預期次數(shù),θi為區(qū)域i的事故險情發(fā)生相對風險值。參數(shù)ei可表示為:

    ei=rk pi(3)

    式中,r為整個區(qū)域的事故險情發(fā)生概率,pi為區(qū)域i的船舶流量。r可以通過事故險情發(fā)生次數(shù)和每個區(qū)域i的交通流量表示為:

    r=5Σi=1 yi / 5Σi=1pi(4)

    在本文中,事故險情發(fā)生相對風險值θ是評估事故險情發(fā)生風險概率的關鍵參數(shù)。以空間和時間參數(shù)為橋梁,計算時空交互效應對潛在事故險情發(fā)生風險的影響。為了比較貝葉斯模型的性能,從考慮或不考慮時空交互效應出發(fā),本文提出了3種類型的貝葉斯時空模型,如模型(5)、(6)、(7)所示。在實踐中,通過DIC (Deviance Information Criterion)值來評價擬合優(yōu)度,DIC是對擬合優(yōu)度性能的評價標準,廣泛應用于多個候選模型之間的比較[13]。具體而言,DIC值越小代表模型的性能越好。

    θij=exp(a0+a1*timelj+ui+gi+di*timelj)(5)

    θij=exp(a0+a1*timelj+ui+vi+gj+ditimelj)(6)

    θij=exp(a0+a1*timelj+ui+vi+gj+psiij+ditimelj(7)

    其中參數(shù)j表示年份,α0為表示平均事故險情發(fā)生相對風險的截距參數(shù),由馬爾可夫鏈蒙特卡羅(MCMC)方法得到,α1為時間效應系數(shù),time1j為j年的時間效應,ui為由空間因素引起的隨機效應。在模型(6)和(7)中,vi分別表示非空間因素引起的隨機效應,gj表示自回歸時間效應,di表示位置效應,psiij表示時空交互效應。

    3 結(jié)果與分析

    經(jīng)求解,模型(5)的DIC值為130.6548,模型(6)的DIC值為129.9314,模型(7)的DIC值為115.9526。其中模型(7)的DIC值最小,因而表明時空交互效應psiij對船舶事故險情發(fā)生有較大影響,且模型(7)是評估事故險情風險的最佳模型。

    根據(jù)事故險情報告歸類得到的各區(qū)域發(fā)生的事故險情位置如圖 3所示,可以直觀地發(fā)現(xiàn),區(qū)域3和區(qū)域5發(fā)生事故險情的頻率較高,其次是區(qū)域4,區(qū)域1與區(qū)域2發(fā)生事故險情的頻率較低。從圖中也可以看出,在警戒區(qū)及警戒區(qū)以南事故險情發(fā)生的頻率較高,警戒區(qū)以北相對較低。

    利用所選模型(7),計算各研究區(qū)域的相對風險值θ,如表1所示。根據(jù)事故險情發(fā)生相對風險值,將5個區(qū)域劃分為低、中、高3個事故險情發(fā)生風險等級。

    根據(jù)表1中結(jié)果,將事故險情發(fā)生相對風險(RCR)值小于1.00的區(qū)域定義為低風險區(qū)域,RCR值在[1.00,1.50)之間的區(qū)域定義為中風險區(qū)域,RCR值大于1.50的區(qū)域定義為高風險區(qū)域。因此,區(qū)域3為高風險區(qū)域,區(qū)域4和5為中風險區(qū)域,區(qū)域1和區(qū)域2屬于低風險區(qū)域。

    圖4展示了各區(qū)域相對風險值逐年變化情況。觀察圖中可知,區(qū)域3、區(qū)域4和區(qū)域5的相對風險值均于2018年達到一個峰值,表明2018年可能是事故險情高發(fā)的一年,水域風險水平較高。區(qū)域3的相對風險水平波動較大,表明區(qū)域3事故險情發(fā)生的數(shù)量在歷年的分布并不均勻,可能是由于航行規(guī)則的改變,導致交通流量發(fā)生變化,進而影響事故險情的發(fā)生。除區(qū)域3外,其他區(qū)域的相對風險水平較為平穩(wěn),波動較小,表明這些水域較為穩(wěn)定,整體航行安全水平變化不大。

    2017年1月至2023年4月事故險情數(shù)變化情況如圖5所示,可以看出,事故險情數(shù)在2018年達到高峰,全年共發(fā)生54起。然而,隨后的幾年,即2019年至2021年,事故險情數(shù)量明顯減少。這可能是由于疫情導致了社會活動的限制以及各項安全措施的加強,從而降低了事故的發(fā)生率。而隨著疫情的穩(wěn)定與好轉(zhuǎn),2022年雖然事故數(shù)量未發(fā)生明顯變化,但險情數(shù)量明顯上升,導致事故險情數(shù)再次呈現(xiàn)上升趨勢。隨著經(jīng)濟的復蘇以及海上貿(mào)易的逐漸恢復,洋山海事局VTS管理服務區(qū)來往的交通流量也在逐步恢復疫情前的水平。因此,為防止事故險情數(shù)的進一步攀升,須繼續(xù)加強安全措施和管理,以確保航行安全和避免事故的發(fā)生。

    船舶事故險情在2017年至2023年間呈現(xiàn)出較為明顯的季節(jié)性分布趨勢,如圖6所示,從事故險情總量而言,春季和夏季是高發(fā)季節(jié),而秋季和冬季則相對較低,但在每一年里面季節(jié)性特點并不完全相同。每年春季洋山水域為大風及能見度不良天氣多發(fā)季節(jié),一定程度上導致了險情事故多發(fā)的客觀現(xiàn)象,而夏季主要受高溫、短時強對流天氣以及臺風影響,風險因素較多,海上風浪較大,增加了事故發(fā)生的可能性。

    為了更加詳細地總結(jié)險情事故在時間維度上的分布特點,本文進一步分析了2017至2023年來事故險情數(shù)在各月的分布情況,如圖7??芍鹿孰U情在3月發(fā)生的頻率較大,達到了月平均4.3起。其次為8月和10月,事故險情數(shù)為月平均2.3起。因此,應格外注意水域在春夏兩季的事故發(fā)生情況,特別是3月、8月與10月,采取一定措施預防事故的發(fā)生。

    4 結(jié) 論

    厘清事故險情發(fā)生時空分布特點在海事部門對水域進行針對性監(jiān)管方面具有重要意義。本文基于貝葉斯時空模型,考慮時空交互效應,對洋山深水港及其附近水域事故險情發(fā)生情況進行建模分析。研究結(jié)果顯示,小洋山港區(qū)以及東海大橋附近水域的相對風險水平最為突出。這一區(qū)域的風險水平呈現(xiàn)出逐年波動的趨勢,航行安全水平相對不穩(wěn)定。大戟洋及北鼎星附近水域的航行風險居次高,蘆潮港沿岸碼頭附近水域船舶航行風險最低。在洋山港及其附近水域,事故險情多發(fā)于春夏兩季,其中主要分布在3月、8月和10月。

    基于模型分析結(jié)果,針對高風險區(qū)域,可以采取更頻繁的巡航巡檢、加強導航警示設施建設等措施,以提高航行安全性。同時,對于低風險區(qū)域,可以適度調(diào)配監(jiān)管資源,優(yōu)化監(jiān)管策略,以確保資源的最優(yōu)利用。

    綜上所述,本文深入分析了洋山深水港及其周邊水域事故險情的時空分布特點,識別水域安全風險,為小洋山北作業(yè)區(qū)的規(guī)劃與建設后,海事部門水上交通安全組織提供有力保障。未來的研究可以考慮更多影響事故發(fā)生的因素,進一步提升預測模型的準確性和實用性。這將有助于在海域管理中更好地平衡航行效率與安全性,確保水域的可持續(xù)發(fā)展。

    參考文獻

    [1] 趙一飛.上海國際航運中心洋山深水港區(qū)國際中轉(zhuǎn)業(yè)務的新里程碑[J].集裝箱化,2021,32(12):1-2+19.DOI:10.13340/j.cont.2021.12.001.

    [2] 洋山深水港區(qū)四期工程上海國際航運中心[J].城鄉(xiāng)建設,2021(18):72-73.

    [3] 陳煒俊.上海港洋山深水港區(qū)發(fā)展對策[J].集裝箱化,2012,23(02):16-19.

    [4] 王桂英.發(fā)展國際集裝箱水水中轉(zhuǎn)助推上海國際航運中心建設——以洋山港為例[J].物流工程與管理,2012,34(05):1-2+7.

    [5] 黃少卿.洋山港——名副其實的東方大港[J].港口經(jīng)濟,2016(07):53-55.

    [6] 鄒齊. 基于貝葉斯分層時空模型的新冠疫情時空相依性研究[D].上海師范大學,2022.DOI:10.27312/d.cnki.gshsu.2022.000759.

    [7] 董文錢,董良,向琳等.基于貝葉斯時空模型的城管事件時空變化分析[J].地球信息科學學報,2020,22(05):1073-1082.

    [8] 徐冰. 基于貝葉斯網(wǎng)絡的傳染病時空預警模型研究[D].長安大學,2016.

    [9] Yu Y, Chen L, Shu Y, et al. Evaluation model and management strategy for reducing pollution caused by ship collision in coastal waters[J]. Ocean & Coastal Management, 2021, 203: 105446.

    [10] 梁永玉,田茂再.基于分層貝葉斯時空Poisson模型的流行病建模研究[J].系統(tǒng)科學與數(shù)學,2022,42(02):462-472.

    [11] 徐善雷,郭錦春.洋山港海事局:大數(shù)據(jù)應用讓水域更安全[J].中國海事,2018(11):70-71.DOI:10.16831/j.cnki.issn1673-2278.2018.11.023.

    [12] 趙昕. 基于貝葉斯時空模型的甘肅省猩紅熱發(fā)病時空分布及與氣象等因素的關聯(lián)研究[D].蘭州大學,2023.

    [13] 徐單. 基于貝葉斯新方法的時空獨立模型影響分析[D].東南大學,2016.董天翔. 江蘇省高速路網(wǎng)道路低溫時空分布及預警模型對比研究[D].南京信息工程大學,2018.

    猜你喜歡
    洋山港風險
    浙江省嵊泗縣洋山港海事局拆解涉漁“三無”船
    大型船舶進出洋山港四期碼頭交通組織
    中國航海(2018年4期)2019-01-07 03:34:32
    上海洋山港四期:智能碼頭的“中國芯”
    金橋(2018年4期)2018-09-26 02:24:50
    上海洋山港區(qū)四期全自動化碼頭開港試運行
    交通與港航(2018年1期)2018-02-27 07:57:07
    我國P2P網(wǎng)絡借貸的風險和監(jiān)管問題研究
    商(2016年27期)2016-10-17 06:18:10
    淺析應收賬款的產(chǎn)生原因和對策
    商(2016年27期)2016-10-17 05:41:05
    中國經(jīng)濟轉(zhuǎn)型的結(jié)構性特征、風險與效率提升路徑
    商(2016年27期)2016-10-17 05:33:32
    互聯(lián)網(wǎng)金融的風險分析與管理
    企業(yè)納稅籌劃風險及防范措施
    上海自貿(mào)區(qū)洋山港集裝箱碼頭
    18禁裸乳无遮挡动漫免费视频| 久久精品国产自在天天线| 最黄视频免费看| 亚洲成人中文字幕在线播放| 日本wwww免费看| 18禁在线播放成人免费| 久久精品国产亚洲网站| 人妻一区二区av| 永久免费av网站大全| 九九久久精品国产亚洲av麻豆| 在线观看国产h片| 丝袜脚勾引网站| 久久久久久久亚洲中文字幕| 免费av不卡在线播放| 国产精品偷伦视频观看了| 九色成人免费人妻av| 夜夜骑夜夜射夜夜干| 亚洲美女搞黄在线观看| 欧美最新免费一区二区三区| 精品国产乱码久久久久久小说| 欧美一区二区亚洲| 国产无遮挡羞羞视频在线观看| a级一级毛片免费在线观看| 免费观看在线日韩| 亚州av有码| 亚洲,一卡二卡三卡| 国产精品一二三区在线看| 亚洲精品国产av蜜桃| 亚洲人成网站高清观看| 国产白丝娇喘喷水9色精品| 国产av码专区亚洲av| 99久久精品热视频| 国语对白做爰xxxⅹ性视频网站| 午夜精品国产一区二区电影| 黄色日韩在线| 十分钟在线观看高清视频www | 中国三级夫妇交换| 最近中文字幕高清免费大全6| 伦理电影免费视频| av国产免费在线观看| 精品国产露脸久久av麻豆| 人人妻人人看人人澡| 国产毛片在线视频| 欧美最新免费一区二区三区| 日韩av不卡免费在线播放| 日本午夜av视频| 亚洲精品色激情综合| 国产精品秋霞免费鲁丝片| 亚洲欧美成人综合另类久久久| 国产免费福利视频在线观看| 久久久a久久爽久久v久久| 国产深夜福利视频在线观看| 涩涩av久久男人的天堂| 欧美变态另类bdsm刘玥| 全区人妻精品视频| 2018国产大陆天天弄谢| 国产中年淑女户外野战色| 美女中出高潮动态图| 大香蕉久久网| 日韩国内少妇激情av| 日韩,欧美,国产一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 美女xxoo啪啪120秒动态图| 少妇裸体淫交视频免费看高清| 欧美xxxx性猛交bbbb| 国产黄色视频一区二区在线观看| 三级国产精品片| 国产午夜精品一二区理论片| 亚洲精品自拍成人| 国产黄片视频在线免费观看| 日韩免费高清中文字幕av| 麻豆成人午夜福利视频| 91精品伊人久久大香线蕉| 欧美一区二区亚洲| 欧美日韩视频高清一区二区三区二| 高清av免费在线| 国产精品欧美亚洲77777| 国产精品精品国产色婷婷| 一区二区三区四区激情视频| 色婷婷久久久亚洲欧美| 国产v大片淫在线免费观看| 亚洲国产毛片av蜜桃av| 十八禁网站网址无遮挡 | 欧美少妇被猛烈插入视频| 狂野欧美白嫩少妇大欣赏| 欧美高清成人免费视频www| 黄色日韩在线| 十八禁网站网址无遮挡 | h视频一区二区三区| 大又大粗又爽又黄少妇毛片口| 久久毛片免费看一区二区三区| 内射极品少妇av片p| 国产免费福利视频在线观看| 麻豆国产97在线/欧美| 久久久久精品性色| 欧美人与善性xxx| 国产真实伦视频高清在线观看| 国产高清有码在线观看视频| av黄色大香蕉| 午夜福利网站1000一区二区三区| 伦理电影免费视频| 国产精品久久久久久精品古装| 在线精品无人区一区二区三 | 国产一级毛片在线| 日日啪夜夜爽| 女的被弄到高潮叫床怎么办| 天堂俺去俺来也www色官网| 少妇人妻精品综合一区二区| 狂野欧美激情性bbbbbb| 91精品国产九色| 最近最新中文字幕大全电影3| 欧美老熟妇乱子伦牲交| 精品国产露脸久久av麻豆| 国内少妇人妻偷人精品xxx网站| 街头女战士在线观看网站| 久久久久精品久久久久真实原创| 亚洲欧美精品专区久久| 国产精品一区二区三区四区免费观看| 欧美区成人在线视频| 纵有疾风起免费观看全集完整版| 一区二区三区精品91| 亚洲精品久久久久久婷婷小说| 99热国产这里只有精品6| 久久综合国产亚洲精品| av线在线观看网站| 欧美精品一区二区大全| 成人18禁高潮啪啪吃奶动态图 | 国产一区亚洲一区在线观看| 亚洲精品成人av观看孕妇| 日韩国内少妇激情av| 日韩制服骚丝袜av| 国产人妻一区二区三区在| 日韩一区二区视频免费看| 国产极品天堂在线| 欧美成人a在线观看| av黄色大香蕉| 欧美精品一区二区大全| 黄色怎么调成土黄色| 纯流量卡能插随身wifi吗| 大香蕉久久网| xxx大片免费视频| 日韩 亚洲 欧美在线| 国产精品人妻久久久久久| 国产精品偷伦视频观看了| av在线播放精品| 中文字幕亚洲精品专区| 久久久欧美国产精品| 伦理电影大哥的女人| 婷婷色麻豆天堂久久| 我的女老师完整版在线观看| 午夜激情久久久久久久| 久久综合国产亚洲精品| 男人爽女人下面视频在线观看| 国产v大片淫在线免费观看| 中文资源天堂在线| 日本黄大片高清| 免费不卡的大黄色大毛片视频在线观看| 久久人人爽av亚洲精品天堂 | 在线观看一区二区三区激情| 九草在线视频观看| 日韩大片免费观看网站| 国产精品一二三区在线看| 内射极品少妇av片p| 久久精品久久久久久噜噜老黄| 亚洲婷婷狠狠爱综合网| 国产精品国产av在线观看| 婷婷色综合大香蕉| av女优亚洲男人天堂| 亚洲无线观看免费| 少妇人妻久久综合中文| 建设人人有责人人尽责人人享有的 | 中国美白少妇内射xxxbb| av天堂中文字幕网| 国产午夜精品久久久久久一区二区三区| 成人黄色视频免费在线看| 国产一区二区在线观看日韩| 国产高清国产精品国产三级 | 久久精品国产亚洲av涩爱| 91午夜精品亚洲一区二区三区| 91精品国产九色| 另类亚洲欧美激情| 一区二区av电影网| 2021少妇久久久久久久久久久| 成年av动漫网址| 大又大粗又爽又黄少妇毛片口| 久久精品国产亚洲av涩爱| 国产男女超爽视频在线观看| 亚洲精品久久久久久婷婷小说| 色婷婷av一区二区三区视频| 亚洲精品国产成人久久av| 亚洲va在线va天堂va国产| 熟女人妻精品中文字幕| 九九爱精品视频在线观看| 亚洲国产精品999| 亚洲激情五月婷婷啪啪| av又黄又爽大尺度在线免费看| 日本av手机在线免费观看| 久热久热在线精品观看| 男女无遮挡免费网站观看| 国产探花极品一区二区| 午夜福利高清视频| 我的女老师完整版在线观看| 王馨瑶露胸无遮挡在线观看| 韩国高清视频一区二区三区| 联通29元200g的流量卡| 国产精品一区二区在线观看99| 日本-黄色视频高清免费观看| 国产女主播在线喷水免费视频网站| 在线观看一区二区三区| 在线观看免费高清a一片| 在线观看一区二区三区激情| 草草在线视频免费看| 狂野欧美白嫩少妇大欣赏| av视频免费观看在线观看| 黑丝袜美女国产一区| 男女边摸边吃奶| 欧美日韩视频精品一区| 搡女人真爽免费视频火全软件| 中文精品一卡2卡3卡4更新| 国产黄片视频在线免费观看| 麻豆精品久久久久久蜜桃| 少妇高潮的动态图| 国产精品精品国产色婷婷| 免费黄频网站在线观看国产| 美女主播在线视频| 国产精品久久久久久精品电影小说 | 久久精品久久久久久噜噜老黄| 精品人妻熟女av久视频| 99视频精品全部免费 在线| 成人18禁高潮啪啪吃奶动态图 | 亚洲av免费高清在线观看| 国产成人一区二区在线| 麻豆国产97在线/欧美| 亚洲不卡免费看| 亚洲精品自拍成人| 亚洲国产色片| 插逼视频在线观看| 偷拍熟女少妇极品色| 涩涩av久久男人的天堂| 男女免费视频国产| 日韩伦理黄色片| 国产黄片视频在线免费观看| 一区二区三区乱码不卡18| 国产v大片淫在线免费观看| 内射极品少妇av片p| 街头女战士在线观看网站| 国产一区亚洲一区在线观看| 日本爱情动作片www.在线观看| 99热6这里只有精品| 欧美极品一区二区三区四区| 午夜视频国产福利| 人体艺术视频欧美日本| 国内揄拍国产精品人妻在线| 精品久久久精品久久久| 久久这里有精品视频免费| 日韩欧美 国产精品| 亚洲一级一片aⅴ在线观看| 啦啦啦视频在线资源免费观看| 成年人午夜在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 成人免费观看视频高清| 一边亲一边摸免费视频| 少妇精品久久久久久久| 久久久久久伊人网av| 美女高潮的动态| 身体一侧抽搐| 99久久综合免费| 成人美女网站在线观看视频| 天堂中文最新版在线下载| 日韩伦理黄色片| 欧美另类一区| 国产精品熟女久久久久浪| 久久久久久人妻| 国产av码专区亚洲av| 日韩三级伦理在线观看| 最后的刺客免费高清国语| 夜夜爽夜夜爽视频| 亚洲真实伦在线观看| 久久久久人妻精品一区果冻| 性色av一级| 最近中文字幕2019免费版| 一级二级三级毛片免费看| 成人美女网站在线观看视频| 一本—道久久a久久精品蜜桃钙片| 嘟嘟电影网在线观看| 婷婷色麻豆天堂久久| 色婷婷久久久亚洲欧美| 亚洲欧美精品自产自拍| 性色avwww在线观看| 精品少妇黑人巨大在线播放| 欧美另类一区| 国产av一区二区精品久久 | 三级国产精品片| 边亲边吃奶的免费视频| 大香蕉97超碰在线| 国产有黄有色有爽视频| 亚洲高清免费不卡视频| 精品酒店卫生间| 久久久久人妻精品一区果冻| 夜夜爽夜夜爽视频| 日日撸夜夜添| 成年av动漫网址| 91aial.com中文字幕在线观看| 国产精品久久久久久精品电影小说 | 美女xxoo啪啪120秒动态图| 在线观看免费视频网站a站| 久久久久久久久大av| 亚洲av不卡在线观看| 国产精品偷伦视频观看了| 日韩 亚洲 欧美在线| 欧美xxⅹ黑人| 久久久色成人| av在线蜜桃| 欧美日韩国产mv在线观看视频 | 自拍欧美九色日韩亚洲蝌蚪91 | 老司机影院毛片| 日本av手机在线免费观看| 男人添女人高潮全过程视频| 五月玫瑰六月丁香| 少妇精品久久久久久久| 一本色道久久久久久精品综合| 久久久久久人妻| 一区在线观看完整版| 久久99热6这里只有精品| 日韩av不卡免费在线播放| 国产淫片久久久久久久久| 人体艺术视频欧美日本| 我要看日韩黄色一级片| 妹子高潮喷水视频| 国产欧美亚洲国产| 九九久久精品国产亚洲av麻豆| 免费观看的影片在线观看| 午夜免费观看性视频| 亚洲精品,欧美精品| 在线 av 中文字幕| 久久久久久九九精品二区国产| 毛片女人毛片| 精品熟女少妇av免费看| 久久国内精品自在自线图片| 肉色欧美久久久久久久蜜桃| 春色校园在线视频观看| 一级毛片黄色毛片免费观看视频| 大香蕉97超碰在线| 两个人的视频大全免费| 日韩 亚洲 欧美在线| 成人毛片a级毛片在线播放| 丰满少妇做爰视频| 插逼视频在线观看| 99久国产av精品国产电影| 欧美最新免费一区二区三区| 插阴视频在线观看视频| 欧美xxⅹ黑人| 99久久综合免费| 亚洲精品一二三| 欧美日韩视频精品一区| 夜夜骑夜夜射夜夜干| 精品少妇黑人巨大在线播放| 成人国产av品久久久| 亚洲欧美日韩无卡精品| 少妇人妻精品综合一区二区| 国产一区二区三区综合在线观看 | 男女啪啪激烈高潮av片| 亚洲自偷自拍三级| 中文乱码字字幕精品一区二区三区| 赤兔流量卡办理| 高清日韩中文字幕在线| 热99国产精品久久久久久7| 精品酒店卫生间| 日韩大片免费观看网站| 乱系列少妇在线播放| a级一级毛片免费在线观看| 亚洲精品日韩av片在线观看| 欧美日韩国产mv在线观看视频 | 在线观看免费日韩欧美大片 | 欧美精品一区二区大全| 六月丁香七月| 欧美+日韩+精品| 有码 亚洲区| 久久精品国产a三级三级三级| 国产精品偷伦视频观看了| 成年人午夜在线观看视频| 人妻制服诱惑在线中文字幕| av又黄又爽大尺度在线免费看| 亚洲精品456在线播放app| 国产精品一区二区性色av| 一二三四中文在线观看免费高清| 亚洲久久久国产精品| 国产精品精品国产色婷婷| 国产伦在线观看视频一区| 丰满迷人的少妇在线观看| 国产成人一区二区在线| 国产精品三级大全| 伊人久久国产一区二区| 亚洲av日韩在线播放| 最近中文字幕高清免费大全6| 国产av一区二区精品久久 | 久久 成人 亚洲| 国产大屁股一区二区在线视频| freevideosex欧美| 国模一区二区三区四区视频| 精品久久久噜噜| 美女主播在线视频| 国产精品嫩草影院av在线观看| 久久久久久久久久成人| 一级片'在线观看视频| 日韩电影二区| 成人18禁高潮啪啪吃奶动态图 | 亚洲av综合色区一区| 日本午夜av视频| 内地一区二区视频在线| 免费不卡的大黄色大毛片视频在线观看| 在线观看一区二区三区激情| 久久久精品94久久精品| 人人妻人人爽人人添夜夜欢视频 | 街头女战士在线观看网站| 91久久精品国产一区二区成人| 亚洲三级黄色毛片| 哪个播放器可以免费观看大片| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久久电影| 欧美xxxx性猛交bbbb| 身体一侧抽搐| 久久精品国产亚洲网站| 久久精品国产自在天天线| 九九爱精品视频在线观看| 成人国产av品久久久| 国产乱来视频区| 欧美精品一区二区免费开放| 女的被弄到高潮叫床怎么办| 五月伊人婷婷丁香| 美女高潮的动态| 亚洲精品aⅴ在线观看| 欧美zozozo另类| 精品久久久噜噜| 男女下面进入的视频免费午夜| 国产美女午夜福利| 日本-黄色视频高清免费观看| 精品一区二区三区视频在线| 自拍偷自拍亚洲精品老妇| 女的被弄到高潮叫床怎么办| 色吧在线观看| 在线观看免费视频网站a站| 我要看日韩黄色一级片| 色综合色国产| 99久久精品国产国产毛片| 22中文网久久字幕| 婷婷色av中文字幕| 国产精品秋霞免费鲁丝片| 国产精品麻豆人妻色哟哟久久| 国产亚洲午夜精品一区二区久久| 99re6热这里在线精品视频| 欧美zozozo另类| 国产有黄有色有爽视频| 狂野欧美白嫩少妇大欣赏| 成年免费大片在线观看| 九草在线视频观看| a级毛片免费高清观看在线播放| av线在线观看网站| 日本色播在线视频| 成年av动漫网址| 久久久久视频综合| 欧美变态另类bdsm刘玥| 欧美zozozo另类| 在线天堂最新版资源| 日韩 亚洲 欧美在线| 成人黄色视频免费在线看| 精品亚洲乱码少妇综合久久| 各种免费的搞黄视频| 丝袜喷水一区| 十分钟在线观看高清视频www | 日韩免费高清中文字幕av| 国产探花极品一区二区| 成人影院久久| 青春草视频在线免费观看| 日韩三级伦理在线观看| 男人狂女人下面高潮的视频| 亚洲av中文av极速乱| 水蜜桃什么品种好| 22中文网久久字幕| 99久国产av精品国产电影| 亚洲成人av在线免费| 国产免费视频播放在线视频| 亚洲欧美成人精品一区二区| 亚洲国产精品专区欧美| 男人狂女人下面高潮的视频| 天天躁夜夜躁狠狠久久av| 日本欧美国产在线视频| 久久热精品热| 欧美xxxx性猛交bbbb| 亚洲电影在线观看av| 国产精品久久久久久久电影| 观看美女的网站| kizo精华| 一级a做视频免费观看| 成人综合一区亚洲| 久久久久久久久大av| 亚洲图色成人| 日韩中字成人| 国产亚洲一区二区精品| 2022亚洲国产成人精品| 99热6这里只有精品| 亚洲av欧美aⅴ国产| 男人和女人高潮做爰伦理| 黑丝袜美女国产一区| 王馨瑶露胸无遮挡在线观看| 国产精品.久久久| 少妇裸体淫交视频免费看高清| 久久久久久久久久人人人人人人| 精品一区二区三区视频在线| 日本午夜av视频| 国产日韩欧美在线精品| 黄色怎么调成土黄色| 免费观看av网站的网址| 建设人人有责人人尽责人人享有的 | 国产精品秋霞免费鲁丝片| 成人综合一区亚洲| 亚洲色图综合在线观看| 午夜福利网站1000一区二区三区| 极品教师在线视频| 啦啦啦在线观看免费高清www| 亚洲在久久综合| 亚洲va在线va天堂va国产| 亚洲精品亚洲一区二区| 亚洲国产精品国产精品| 搡老乐熟女国产| 久久久久久久久久久免费av| 国内精品宾馆在线| 男人舔奶头视频| 亚洲美女搞黄在线观看| 国产亚洲午夜精品一区二区久久| 国产白丝娇喘喷水9色精品| 亚洲国产色片| 国内揄拍国产精品人妻在线| av在线app专区| 啦啦啦在线观看免费高清www| 欧美一级a爱片免费观看看| 国产精品人妻久久久影院| 精品亚洲成国产av| 在线观看三级黄色| av在线app专区| 成人黄色视频免费在线看| 嘟嘟电影网在线观看| 久久精品久久久久久噜噜老黄| 高清毛片免费看| 在线观看一区二区三区| 一级av片app| 日韩伦理黄色片| 纯流量卡能插随身wifi吗| 国产成人91sexporn| 最近最新中文字幕大全电影3| 成人特级av手机在线观看| 欧美最新免费一区二区三区| 亚洲精品成人av观看孕妇| 99热网站在线观看| 伦精品一区二区三区| 色婷婷久久久亚洲欧美| 26uuu在线亚洲综合色| 成人午夜精彩视频在线观看| 亚洲国产最新在线播放| 一区二区三区精品91| 99热这里只有是精品50| 亚洲,一卡二卡三卡| 欧美精品人与动牲交sv欧美| 99久久人妻综合| 亚洲国产精品一区三区| 天堂中文最新版在线下载| 国产精品蜜桃在线观看| 日韩av免费高清视频| 日本黄色日本黄色录像| 熟女人妻精品中文字幕| 美女国产视频在线观看| 久久av网站| 日韩亚洲欧美综合| 久久精品久久久久久噜噜老黄| 亚洲人成网站在线观看播放| 性色avwww在线观看| 欧美日韩视频精品一区| 秋霞伦理黄片| 一级毛片我不卡| 亚洲美女黄色视频免费看| 联通29元200g的流量卡| 好男人视频免费观看在线| 久久99蜜桃精品久久| 中文字幕久久专区| 99久久综合免费| 精品熟女少妇av免费看| 最近的中文字幕免费完整| 国产真实伦视频高清在线观看| 高清毛片免费看| 免费不卡的大黄色大毛片视频在线观看| 国产欧美日韩精品一区二区| 黄色怎么调成土黄色| 黑丝袜美女国产一区| 亚洲精品国产av成人精品| 夜夜看夜夜爽夜夜摸| 51国产日韩欧美| 日本午夜av视频| 国产亚洲av片在线观看秒播厂| 国产精品免费大片| 亚洲av二区三区四区| 免费在线观看成人毛片| 自拍欧美九色日韩亚洲蝌蚪91 | 国内揄拍国产精品人妻在线| 欧美老熟妇乱子伦牲交| 中文乱码字字幕精品一区二区三区| 亚洲av成人精品一二三区| 国产毛片在线视频| a 毛片基地| 久久久久久伊人网av| 久久鲁丝午夜福利片| 午夜激情福利司机影院| 五月玫瑰六月丁香| 亚洲中文av在线| 韩国av在线不卡| 大片电影免费在线观看免费|