亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于遙感技術(shù)的田灣核電溫排水分布研究

        2021-12-30 06:51:10石海崗章新益張恩董雙發(fā)張建永張春雷
        河北遙感 2021年4期
        關(guān)鍵詞:田灣排水口溫升

        石海崗,章新益,張恩,董雙發(fā),張建永,張春雷

        1.核工業(yè)航測(cè)遙感中心,河北 石家莊050002

        2.河北省航空探測(cè)與遙感技術(shù)重點(diǎn)實(shí)驗(yàn)室,河北 石家莊050002

        0 引言

        隨著公眾環(huán)保意識(shí)增強(qiáng)及國(guó)家產(chǎn)業(yè)結(jié)構(gòu)和能源結(jié)構(gòu)優(yōu)化,我國(guó)將持續(xù)穩(wěn)步推進(jìn)核電建設(shè)。核電作為清潔能源,運(yùn)行過程中不產(chǎn)生有害氣體(如SOx、NOx等)和溫室氣體CO2,但會(huì)形成大量的溫排水,造成周邊海域升溫[1],當(dāng)溫度過高且擴(kuò)散到核電取水口附近時(shí),會(huì)影響核電的正常運(yùn)行。此外,溫升較高時(shí)還會(huì)改變水體的質(zhì)量,影響到水生生物的繁殖和分布情況[2]。海岸線的變遷與海域使用、海岸保護(hù)等活動(dòng)密切相關(guān),對(duì)于海灣型的海域,岸線的變化影響到海域的納潮量和流場(chǎng)[3],間接影響到核電溫排水的分布情況。因此,核電站周邊海域環(huán)境監(jiān)測(cè),尤其是岸線、溫度場(chǎng)信息變化對(duì)核電的安全運(yùn)行和研究水域生態(tài)平衡有著重要意義[4]。遙感作為對(duì)地觀測(cè)綜合性技術(shù),在核電周邊海域環(huán)境監(jiān)測(cè)方面具有其它技術(shù)手段無法比擬的特點(diǎn)[5-6]。

        連云港作為歐亞大陸的東橋頭堡,沿海經(jīng)濟(jì)快速發(fā)展,核電周邊海域涉海工程逐年增加。田灣核電因自身機(jī)組擴(kuò)建也開展了系列工程建設(shè),導(dǎo)致了海域岸線形態(tài)和溫排水排水量變化顯著。海域岸線形態(tài)的改變和核電排水量的增加給水動(dòng)力和水環(huán)境帶來連鎖反應(yīng)[7],導(dǎo)致核電周邊海域環(huán)境發(fā)生顯著變化。以往海域環(huán)境變化研究多集中于分析岸線變遷驅(qū)動(dòng)力,為加強(qiáng)岸線管理提供依據(jù)。本次研究基于Landsat-5 存檔數(shù)據(jù)和Landsat-8 現(xiàn)勢(shì)數(shù)據(jù),利用遙感和地理信息技術(shù)提取田灣核電周邊海域海岸線數(shù)據(jù)和溫度場(chǎng)分布數(shù)據(jù),研究田灣核電周邊海域岸線變遷、核電擴(kuò)建對(duì)海域環(huán)境溫度場(chǎng)的影響。

        1 研究區(qū)與數(shù)據(jù)源

        1.1 研究區(qū)概況

        田灣核電廠址位于江蘇省連云港市高公島鄉(xiāng),廠區(qū)規(guī)劃建設(shè)8臺(tái)百萬千瓦級(jí)核電機(jī)組,一期1、2號(hào)機(jī)組和二期3、4號(hào)機(jī)組分別在2007年和2018年投入商運(yùn),三期5、6號(hào)機(jī)組分別在2020年9月和2021年5月投入商運(yùn),6臺(tái)機(jī)組單機(jī)容量均超過100萬千瓦。核電周邊海域向東開放,近岸區(qū)域?qū)儆倌喾植紖^(qū),南北兩側(cè)現(xiàn)均已被港口工程包圍(圖1,見中間彩頁(yè),下同)。核電站周邊為規(guī)則淺海半日潮型,一日經(jīng)歷2 次漲落潮,海域流場(chǎng)運(yùn)動(dòng)形式以往復(fù)流或接近往復(fù)流為主,歷時(shí)略有不同。落潮時(shí),海水自核電排水口沿導(dǎo)流堤向南流動(dòng),然后受潮汐影響向東北方向?yàn)a出,分別流經(jīng)排水口導(dǎo)流堤,取水明渠,連云港防波堤,至連島鎮(zhèn)北側(cè)。漲潮流大致與落潮流相反。

        1.2 遙感數(shù)據(jù)源

        本文數(shù)據(jù)源為2003 年12 月22 日、2005 年9 月6日和2008年2月19日的Landsat-5數(shù)據(jù),2013年11月15日、2017年2月27日和2020年3月23日的Land‐sat-8數(shù)據(jù)。為驗(yàn)證反演結(jié)果,獲取了2008年2月19日過境的MODIS 1KM-Level1B Calibrated Radi‐ances 數(shù)據(jù)和2013 年11 月15 日據(jù)海面測(cè)量數(shù)據(jù)。衛(wèi)星熱紅外波段參數(shù),衛(wèi)星過境時(shí)間、機(jī)組運(yùn)行工況,所處季節(jié)和潮態(tài)見表1。Landsat-5 衛(wèi)星攜帶的專題制圖儀(TM)、Landsat-8 的OLI 和TIRS 傳感器具體參數(shù)特征見文獻(xiàn)[8]。

        表1 田灣核電不同時(shí)相熱紅外數(shù)據(jù)概況表Table 1 Overview of thermal infrared data of Tianwan nuclear power plant in different time

        2 研究方法

        2.1 遙感數(shù)據(jù)預(yù)處理

        遙感影像預(yù)處理,主要是對(duì)Landsat-8 OLI 和TIRS、Landsat-5 TM 和MODIS 數(shù)據(jù)進(jìn)行輻射定標(biāo)、大氣校正、幾何精校正、海陸分離以及影像的裁剪和重采樣。幾何精校正通過5 萬地形圖數(shù)據(jù)為基準(zhǔn)進(jìn)行,使校正后的誤差控制在5m 之內(nèi);大氣校正通過ENVI 軟件的FLAASH 模塊進(jìn)行;海陸分離利用Landsat-8 OLI 和Landsat-5 TM 短波紅外波段來區(qū)分。對(duì)于MODIS 產(chǎn)品,幾何校正后,重采樣成120m。各影像按照研究區(qū)矢量進(jìn)行裁剪(圖1)。

        2.2 溫度場(chǎng)信息提取

        溫度場(chǎng)信息提取的基礎(chǔ)是對(duì)熱紅外波段開展溫度反演。中外學(xué)者均針對(duì)Landsat-8的2個(gè)熱紅外通道,開展了劈窗算法研究,并進(jìn)行了精度估算[9-10],但因Landsat-8 的B11 波段受條帶干擾嚴(yán)重,數(shù)據(jù)分發(fā)單位(https://glovis.usgs.gov/)經(jīng)過測(cè)試顯示,利用B11 波段的劈窗算法會(huì)帶來較大的誤差。針對(duì)Landat-5 熱紅外數(shù)據(jù),國(guó)內(nèi)外學(xué)者提出不同的反演算法,主要有單窗算法[11]、普適性單通道算法[12]和輻射傳輸方程法[13-14]。前兩種算法是通過對(duì)輻射傳輸方程算法中大氣相關(guān)參數(shù)近似簡(jiǎn)化,獲取地表溫度。為保證結(jié)果的一致性,對(duì)Landsat-8 和Landsat-5 均使用輻射傳輸方程算法開展溫度反演,并對(duì)反演結(jié)果開展驗(yàn)證。

        在無云情況下,不考慮大氣對(duì)電磁波的散射,水平大氣各種組分混合均勻,對(duì)于溫度Ts的條件下,地表的黑體輻射為∶

        式中,Ts是地物溫度,Lλ(Ts)是溫度為Ts時(shí)的黑體輻射,Lλ是傳感器接收到的大氣頂層輻射,可由傳感器輻射定標(biāo)獲取,ελ是地物的比輻射率,因海水接近黑體,比輻射率可以取0.995[11];Lλatm↓和Lλatm↑分別是大氣下行輻射和大氣上行輻射,τλ是波長(zhǎng)為λ時(shí)地表和傳感器之間的大氣透射率。

        (1)式經(jīng)變形,可以得到衛(wèi)星傳感器接收到的大氣頂層輻射Lλ:

        由(2)式可知,要求算地表溫度Ts,除海表比輻射率ελ外,還需計(jì)算5個(gè)參數(shù):大氣頂層輻射Lλ、溫度為Ts時(shí)的黑體輻射Lλ(Ts)、大氣透射率τλ、大氣上行輻射亮度Lλatm↑、大氣下行輻射亮度Lλatm↓。

        1.Lλ的計(jì)算

        大氣頂層輻射Lλ的計(jì)算是熱紅外傳感器輻射定標(biāo)的過程,即將傳感器觀測(cè)到的數(shù)據(jù)灰度(DN)值轉(zhuǎn)換成輻射亮度值的過程。公式如下:

        其中,a和b為定標(biāo)系數(shù),分別為圖像的增益和偏移,可以直接從各自元數(shù)據(jù)中獲取。

        2.大氣參數(shù)Lλatm↓,Lλatm↑和τλ的獲取

        根據(jù)衛(wèi)星過境時(shí)刻核電地面的平均溫度,氣壓,濕度及影像中心經(jīng)緯度,結(jié)合MODTRAN模塊,進(jìn)行在線(http∶//atmcorr.gsfc.nasa.gov/)大氣校正[11-12],獲取Lλatm↓,Lλatm↑和τλ參數(shù)。

        3.Lλ(Ts)和Ts的計(jì)算

        在獲取Lλatm↓,Lλatm↑和τλ后,獲取海域表層真實(shí)的Lλ(Ts),根據(jù)普朗克公式的反函數(shù),獲取海表真實(shí)溫度Ts,即

        式中,對(duì)于

        基于以上的算法和元數(shù)據(jù)中相關(guān)參數(shù),進(jìn)行波段運(yùn)算,獲得6 期數(shù)據(jù)海面的溫度場(chǎng)分布如圖2、表2所示。

        表2 田灣核電周邊海域不同時(shí)期溫度概況表Table 2 Temperature of sea region near Tianwan nuclear power plant at different periods

        2.3 岸線信息提取

        基于地理信息軟件Arcgis,對(duì)大氣校正后的Landsat-5、8多光譜數(shù)據(jù)進(jìn)行圖像融合(Landsat-8數(shù)據(jù))、波段組合和圖像增強(qiáng),采用目視解譯方法對(duì)6期遙感影像進(jìn)行解譯,提取核電周邊海域岸線情況(圖1、表3)。

        表3 田灣核電周邊海域不同時(shí)期岸線概況表Table3 Overview of coastline around Tianwan nucle‐ar power plant at different periods

        解譯結(jié)果顯示,在核電運(yùn)行前,2003~2008 年田灣核電周邊岸線未發(fā)生變化,2008~2017 年間,核電周邊發(fā)生了顯著變化:

        1.連島南側(cè)、進(jìn)水口北側(cè)的旗臺(tái)防波堤,從無到有,最終長(zhǎng)度達(dá)4.81km;

        2.取水明渠逐年擴(kuò)展,最終增加到4.54km;

        3.排水口導(dǎo)流堤向南擴(kuò)建約1.50km;

        4.南側(cè)徐圩港防波堤,從無到有,最終離岸距離達(dá)6.86km;

        5.核電排水口由原來相對(duì)開放,逐漸處于兩側(cè)防波堤環(huán)抱的人工海灣的灣底。

        2017-2020年岸線穩(wěn)定,未發(fā)生變化。

        3 結(jié)果與分析

        2013 年11 月15 日Landsat-8 衛(wèi)星過境期間開展了近同步海面溫度測(cè)量,海溫測(cè)量方式和溫度反演精度的分析研究,作者在文獻(xiàn)[15-16]中已有報(bào)道,本次研究不再贅述。

        3.1 Landsat-5 反演結(jié)果與MODIS 溫度數(shù)據(jù)對(duì)比分析

        因Landsat-5數(shù)據(jù)時(shí)相較早,未開展同步海面實(shí)測(cè)工作。為驗(yàn)證Landsat-5熱紅外波段溫度反演結(jié)果,以精度較高的MODIS數(shù)據(jù)的海表溫度二級(jí)產(chǎn)品[17](圖3)為基礎(chǔ),與Landsa-5反演溫度開展交叉驗(yàn)證[18]。

        因處于落潮末期,近岸灘涂出露面積較大,由于比熱容的差異,灘涂的溫度明顯高于海水溫度。MODIS 熱紅外數(shù)據(jù)空間分辨率為1km,海陸像元混合比Landsat-5要大,為減小誤差,30個(gè)隨機(jī)點(diǎn)位選取時(shí)多遠(yuǎn)離近岸灘涂。此外,2008年2月19日MODIS數(shù)據(jù)過境時(shí)南側(cè)閘口放水,選點(diǎn)時(shí)避開放水影響區(qū)域。將選取的數(shù)據(jù)擬合(表4、圖3),進(jìn)行交叉驗(yàn)證。Landsat-5偏差最大為3.61℃,最小小于0.1℃,平均偏差為0.64℃,兩組數(shù)據(jù)擬合關(guān)系式為y=0.7266x+1.2104,擬合后回歸系數(shù)的平方為0.8856,標(biāo)準(zhǔn)誤差為0.6663。

        表4 Landsat-5 反演結(jié)果與MODIS 反演結(jié)果對(duì)比表Table 4 Comparison of Landsat-5 retrieval results with MODIS temperature

        圖4 Landsat-5 反演結(jié)果與MODIS 溫度數(shù)據(jù)擬合圖Fig.4 Linear fitting of Landsat-5 retrieval results and MODIS temperature values

        數(shù)據(jù)擬合結(jié)果表明Landsat-5 反演結(jié)果與MODIS溫度數(shù)據(jù)之間線性特征明顯,具有很好的線性相關(guān)性和一致性,可以反映出Landsat-5 熱紅外波段溫度反演方法獲得的溫度場(chǎng)數(shù)據(jù)準(zhǔn)確可信。

        3.2 核電周邊溫度場(chǎng)分布特征

        2003年12月22日、2005年9月6日田灣核電站未運(yùn)行,附近海域溫度較為均勻(圖2a,2b),溫度分別在4.0~6.0℃和23.0℃~24.0℃之間,溫度較為均勻,排水口周邊海域無明顯溫度分異與溫升現(xiàn)象。

        2008年后獲取的4景數(shù)據(jù),核電已投入運(yùn)行,排水口周邊溫排水海域溫度層次分明,存在明顯的水溫分異現(xiàn)象,且距離排水口越近溫度越高,隨著距離增大,溫度逐步降低,到達(dá)一定距離后,溫度趨于穩(wěn)定。

        運(yùn)行后的4 期衛(wèi)星數(shù)據(jù),獲取時(shí)核電周邊海域均處在落潮時(shí)刻,溫排水形成的溫度場(chǎng)均具有沿海水落潮方向擴(kuò)散特征。2008 年2 月19 日周邊工程還未開始修建(圖2(c)),溫排水?dāng)U散較快,影響范圍相對(duì)較小。因無工程阻擋,進(jìn)水口處于明顯受到溫排水影響(4.9℃),處于高于周邊海水溫度(3.2℃)的范圍內(nèi)。2013年11月15日海域溫度場(chǎng)顯示(圖2(d)),由于旗臺(tái)防波堤和取水明渠向外海延伸,阻擋溫排水落潮擴(kuò)散通道,相較于2008 年2月19 日,溫度場(chǎng)整體向東偏移,進(jìn)水口平均溫度為15.4℃,略高于周邊海水溫度(15.0℃),溫排水影響的溫度場(chǎng)整體規(guī)模略有增大。

        2017 年2 月27 日溫度場(chǎng)結(jié)果(圖2(e))顯示,高溫?zé)崴芘潘趯?dǎo)流堤阻擋向南擴(kuò)散,過南側(cè)堤頭后,受潮汐拖曳作用,以約1.0km 寬度向東延伸2km,溫度降低2.0℃,繼續(xù)向東擴(kuò)散,因與外海海水混合,溫度迅速降低,2km 范圍內(nèi)降低4.0℃。因取水明渠阻擋,溫排水被限制在進(jìn)水口處南側(cè);向東南擴(kuò)散的溫排水羽跡,被徐圩港防波堤阻斷。溫排水展布進(jìn)一步向東、東南偏移,范圍增大明顯。

        2020 年3 月23 日核電周邊岸線未發(fā)生變化,核電二期3、4號(hào)機(jī)組投入商運(yùn),4 臺(tái)機(jī)組滿功率運(yùn)行,一期排水口西南側(cè)800m 處第二個(gè)排水口開始啟用,2 個(gè)排水口同時(shí)向外排放溫排水。因熱水混合較快,新排水口未形成獨(dú)立溫度場(chǎng)(圖2(f)),展布形態(tài)與2017年2月27日相似,但溫度場(chǎng)規(guī)模明顯擴(kuò)大。

        4 溫升統(tǒng)計(jì)與討論

        4.1 基準(zhǔn)溫度選取與溫升統(tǒng)計(jì)

        遙感監(jiān)測(cè)反映的是核電周邊海域?qū)嶋H水溫分布,是溫排水和自然因素綜合影響的結(jié)果。核電運(yùn)行后溫排水溫升信息提取,應(yīng)將自然水溫從遙感監(jiān)測(cè)的溫度場(chǎng)扣除,從而獲取核電站溫排水的影響范圍。但實(shí)際情況是核電周邊海域隨著季節(jié)、潮汐、核電運(yùn)行工況不同,絕對(duì)溫度不斷變化,基準(zhǔn)溫度不是一個(gè)固定值。如果沒有準(zhǔn)確的基準(zhǔn)溫度作為參考,就無法確定溫升分布范圍的大小及位置,無法進(jìn)行環(huán)境影響評(píng)價(jià)。目前技術(shù)條件下,無法把精確的本底溫度場(chǎng)從附近水體中完全剔除,但可根據(jù)一定的方法如臨近區(qū)域替代法、海灣平均溫度法或取水口法等選取基準(zhǔn)溫度[6]。

        經(jīng)運(yùn)行前溫度場(chǎng)對(duì)比分析,連島北側(cè)的溫度與溫排水區(qū)域溫度基本一致(圖2a、圖2b,表2),因此,本次研究的基準(zhǔn)溫度以連島北側(cè)的海水溫度作為參考,綜合考慮剔除溫排水影響區(qū)域后的核電附近海域平均溫度。將熱紅外反演的海域溫度場(chǎng)剔除基準(zhǔn)溫度,劃分溫升等級(jí)進(jìn)行編碼(圖5),統(tǒng)計(jì)各級(jí)溫升象元個(gè)數(shù),獲取核電溫排水形成的溫升分布結(jié)果(圖6)。

        溫升編碼圖(圖5a、b、c)顯示,隨著周邊岸線變化,取水口處溫升影響逐年降低,但排水口周邊溫排水海域各級(jí)溫升面積總體上均有增大(圖6)。因核電機(jī)組運(yùn)行功率相同,監(jiān)測(cè)季節(jié)及所處潮汐基本相同,可以認(rèn)為溫升面積增大,主要因?yàn)楹穗娭苓叞毒€環(huán)境改變導(dǎo)致。2020 年3 月23 日為核電二期3、4 號(hào)機(jī)組商運(yùn)后數(shù)據(jù)(圖6d),核電周邊工程建設(shè)穩(wěn)定,與2017 年相比,岸線無變化,溫升影響規(guī)模顯著增加,1℃以上溫升增加約12.11%,4℃以上溫升增加約6.98%(圖6)。排除其他影響因素,可推斷溫升面積增加為新增機(jī)組溫排水量增加引起的變化。

        4.2 討論

        隨著周邊工程的變化和核電的擴(kuò)建,核電站周邊海域環(huán)境發(fā)生了顯著的變化。核電周邊工程未建成時(shí),漲潮時(shí)刻,潮水從外海沿東北方向流進(jìn)排水區(qū)域,遇岸線阻擋,向東南方向流動(dòng),落潮時(shí)刻,潮水基本沿相反方向流出;工程建成后,取水明渠在漲潮時(shí)使潮水發(fā)生挑流,落潮時(shí)阻礙了潮水向東北方向的流動(dòng),而徐圩港防波堤阻擋了潮水向南擴(kuò)散。受此影響,溫排水海域漲落潮流速都有不同程度的減小,并且流向發(fā)生了偏轉(zhuǎn),不利于溫排水的擴(kuò)散,造成了溫排水影響的面積增大。2020 年3、4號(hào)機(jī)組投入運(yùn)行后,溫排水量增加明顯,進(jìn)一步改變了周邊海域環(huán)境,溫升也有了明顯的變化,但溫排水向外海擴(kuò)散到達(dá)一定距離后,海水流速,海水深度都明顯高于近岸,有利于熱量擴(kuò)散,因此,核電擴(kuò)建后溫排水影響面積并未按照機(jī)組增加幅度增加。

        如前所述,現(xiàn)有技術(shù)很難把自然溫度場(chǎng)完全分離出來,現(xiàn)有的基準(zhǔn)溫度扣除,會(huì)導(dǎo)致遙感監(jiān)測(cè)得到的溫升范圍可能與實(shí)際溫排水造成的溫升范圍不符,加之海陸混合像元的影響,也會(huì)導(dǎo)致遙感監(jiān)測(cè)得到的溫升范圍略大于實(shí)際溫排水造成的溫升范圍。為了扣除準(zhǔn)確的本底溫度場(chǎng),還需要開展季節(jié)、潮汐及近岸淺灘的影響因素研究。

        5 結(jié)論

        1.通過Landsat-8 反演溫度與近同步測(cè)量的測(cè)溫?cái)?shù)據(jù)對(duì)比、Landsat-5 與同日過境的MODIS 數(shù)據(jù)進(jìn)行交叉驗(yàn)證,顯示利用輻射傳輸方程算法對(duì)Land‐sat-8 和Landsat-5 的溫度反演結(jié)果可靠。田灣核電運(yùn)行前,排水口周邊海域溫度無明顯升溫現(xiàn)象。核電運(yùn)行后,海域溫升現(xiàn)象明顯。

        2.Landsat-5、8 數(shù)據(jù)岸線和溫度場(chǎng)信息提取結(jié)果顯示,田灣核電周邊海域環(huán)境變化明顯,周邊工程建設(shè),尤其是取水明渠的建設(shè)有效保護(hù)了核電冷卻水不受溫排水影響,阻礙了溫排水的擴(kuò)散速度、影響了溫排水?dāng)U散方向,造成溫排水溫升面積增大。核電二期擴(kuò)建后,因溫排水量增加,溫升面積進(jìn)一步擴(kuò)大,因外海流速、深度增大,有利于溫排水熱量擴(kuò)散,溫升幅度未按機(jī)組增加幅度增加。核電站需要時(shí)刻關(guān)注周邊海域環(huán)境變化,必要時(shí)進(jìn)行周邊海域環(huán)境信息提取、溫排水影響分析。

        3.Landsat-5、Landsat-8 數(shù)據(jù)能滿足田灣核電周邊海域環(huán)境監(jiān)測(cè)需求,為評(píng)估核電站周邊海域環(huán)境提供了存檔、現(xiàn)勢(shì)遙感數(shù)據(jù)。

        猜你喜歡
        田灣排水口溫升
        城市道路排水口安全自動(dòng)打開系統(tǒng)的設(shè)計(jì)
        電機(jī)溫升試驗(yàn)分析及無人值守電機(jī)溫升試驗(yàn)優(yōu)化
        電機(jī)溫升計(jì)算公式的推導(dǎo)和應(yīng)用
        田灣核電站運(yùn)行階段質(zhì)量保證分級(jí)管理
        田灣核電站運(yùn)行階段質(zhì)量保證分級(jí)管理
        田灣核電7、8號(hào)機(jī)組首臺(tái)預(yù)引入設(shè)備順利就位
        田灣核電6號(hào)機(jī)組首次并網(wǎng)成功
        排水口工程對(duì)京杭運(yùn)河通航條件的影響研究
        LED照明光源的溫升與散熱分析
        電子制作(2018年2期)2018-04-18 07:13:36
        可裁剪地漏過濾網(wǎng)
        全部孕妇毛片丰满孕妇孕交| 视频一区视频二区自拍偷拍| 亚洲国产一区二区三区| 欧美狠狠入鲁的视频777色| 黄又色又污又爽又高潮动态图 | 无码啪啪熟妇人妻区| 最近高清中文在线字幕观看| 国产精品反差婊在线观看| 久亚洲一线产区二线产区三线麻豆| 亚洲国产精品久久精品 | 日韩精品一区二区午夜成人版| 妇女性内射冈站hdwwwooo| 天堂av无码大芭蕉伊人av孕妇黑人| 毛片精品一区二区二区三区| 日本在线 | 中文| 中文字幕精品久久久久人妻红杏1| 久久精品国产亚洲AV香蕉吃奶| 国产青青草自拍视频在线播放| 手机久草视频福利在线观看| 国产女主播白浆在线观看| 亚洲国产成人精品无码区在线观看| 91自国产精品中文字幕| 国产实拍日韩精品av在线| 国精产品一区一区三区有限公司杨 | 91久久综合精品久久久综合| 国产精品18久久久| 亞洲綜合一區二區三區無碼| av免费在线手机观看| 精品综合久久88少妇激情| 国产国产人免费人成免费视频| 久久人人97超碰超国产| 免费高清视频在线观看视频| 精品国产三级a在线观看不卡| 国产爆乳美女娇喘呻吟| 中文字幕乱码免费视频| 欧洲亚洲色一区二区色99| 极品尤物在线精品一区二区三区 | 视频一区二区三区黄色| 国产免费av片在线播放| 国产艳妇av在线出轨| 日本国产一区二区在线|