999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

蕭山站遷站觀測(cè)資料對(duì)比評(píng)估

2018-04-23 05:36:03馬鳳華耿迪劉學(xué)軍鐘已承
浙江氣象 2018年1期
關(guān)鍵詞:風(fēng)速差異

馬鳳華 耿迪 劉學(xué)軍 鐘已承

(1.蕭山區(qū)氣象局,浙江 杭州 311201;2.陽(yáng)城縣氣象局,山西 晉城 048100)

0 引 言

蕭山國(guó)家一般氣象站始建于1954年5月1日,站址位于蕭山縣城東長(zhǎng)山公社山末址鄉(xiāng)村。1996年1月1日搬遷至北干街道柳橋社區(qū)北干山“山腰”。長(zhǎng)期以來,由于觀測(cè)場(chǎng)北側(cè)緊鄰山坡,西西北到東北方向有山脈遮擋,觀測(cè)場(chǎng)向外拓展空間受限,無(wú)法滿足現(xiàn)代氣象觀測(cè)業(yè)務(wù)發(fā)展需要。經(jīng)中國(guó)氣象局批準(zhǔn)同意,2016年1月1日遷移至新址,按照地面氣象觀測(cè)規(guī)范[1]和遷站規(guī)范[8]規(guī)定,臺(tái)站搬遷需要進(jìn)行為期一年的對(duì)比觀測(cè)。國(guó)內(nèi)已有許多文獻(xiàn)[2~7]對(duì)遷站觀測(cè)資料進(jìn)行了對(duì)比分析,提出遷站是產(chǎn)生氣候資料序列非均一性現(xiàn)象的最主要的原因之一。由于地理環(huán)境的不同和海拔的變化,使得蕭山國(guó)家一般氣象站新、舊站址相關(guān)氣象要素?cái)?shù)據(jù)產(chǎn)生差異。為了研究新址的氣象資料是否具有延續(xù)性和代表性,本文利用蕭山站2016年新、舊址的氣象要素觀測(cè)值,采用差值標(biāo)準(zhǔn)差、降水累計(jì)相對(duì)差值、風(fēng)向相符率、顯著性檢驗(yàn)等分析方法,進(jìn)行對(duì)比分析,尋找差異,并分析其產(chǎn)生的原因。

1 資料和方法

1.1 新、舊站址環(huán)境狀況對(duì)比

舊址位于蕭山城東北干山山腰,地理位置120°17′11″E、30°10′38″N,觀測(cè)場(chǎng)海拔為44.3 m,1996年1月1日正式啟用,西西北到東北方向有山脈阻擋,對(duì)日照和風(fēng)有一些影響,觀測(cè)場(chǎng)北面緊鄰山坡,西—北—東周圍10 m范圍內(nèi)有高于1 m的植物和樹木。新址位于舊址西北方向的北干山山頂,距離舊址198 m,地理位置120°17′03″E、30°10′40″N,觀測(cè)場(chǎng)海拔為96.5 m,周邊空曠,無(wú)高大障礙物遮擋,受城市等特定環(huán)境影響小。新、舊址的觀測(cè)儀器均為DZZ5型自動(dòng)氣象站,軟件均為臺(tái)站地面綜合觀測(cè)業(yè)務(wù)軟件,簡(jiǎn)稱ISOS。

1.2 資料來源

對(duì)比分析數(shù)據(jù)使用2016年1—12月蕭山站新址和舊址的逐日觀測(cè)資料,差異顯著性檢驗(yàn)使用1996—2015年舊址和2016年新址的觀測(cè)資料。對(duì)比分析評(píng)估時(shí),若新址或舊址觀測(cè)資料有缺測(cè),對(duì)應(yīng)的對(duì)比資料作缺測(cè)處理。

1.3 統(tǒng)計(jì)方法

1.3.1 差值及差值標(biāo)準(zhǔn)差

分別計(jì)算新、舊站址之間的日平均氣溫、日最高氣溫、日最低氣溫、日降水量、日平均相對(duì)濕度、日平均風(fēng)速的差值,并統(tǒng)計(jì)求取月、年差值平均值和差值標(biāo)準(zhǔn)差。差值平均值、差值的標(biāo)準(zhǔn)差計(jì)算公式分別為:

(1)

(2)

式中xi為新、舊站址第i次要素觀測(cè)值的差值,n為樣本容量。

1.3.2 降水量累計(jì)相對(duì)差值

降水量累計(jì)相對(duì)差值反應(yīng)了兩站降水總體差異情況。采用以下公式計(jì)算:

(3)

(3)式中,Ai為舊址第i次觀測(cè)降水量,xi為新舊址第i次降水量觀測(cè)值的差值;n為總降水次數(shù)。

1.3.3 風(fēng)向相符率統(tǒng)計(jì)

只有觀測(cè)風(fēng)速大于0.2 m/s時(shí),才統(tǒng)計(jì)風(fēng)向相符率。新址與舊址風(fēng)向角度差小于22.5°,即認(rèn)為兩者相符。采用以下公式計(jì)算:

(4)

1.3.4 顯著性檢驗(yàn)

利用舊址近20 a的觀測(cè)要素?cái)?shù)據(jù)月(年)平均值序列,對(duì)新址平行觀測(cè)期的月(年)平均值進(jìn)行顯著性檢驗(yàn)。采用以下公式計(jì)算:

假設(shè)x0為平行觀測(cè)期的某要素月(年)平均值,假定樣本容量為n月(年),這n年觀測(cè)的相應(yīng)要素月(年)平均值的樣本序列為x1,x2,…,xi,…,xn,統(tǒng)計(jì)量

(5)

在顯著水平a為0.05(0.01),自由度為n時(shí),當(dāng)|t|>ta,則表示平均值x0與其它年份平均值差異顯著。

2 新、舊站址各氣象要素差值分析

2.1 差值及差值標(biāo)準(zhǔn)差

分別計(jì)算新址與舊址之間的日平均氣溫、日最高氣溫、日最低氣溫、日降水量、日平均相對(duì)濕度、日平均風(fēng)速的差值,并統(tǒng)計(jì)求取1—12月、年差值平均值(見表1)和差值標(biāo)準(zhǔn)差(見表2)。在表1中,差值為正,代表新址的要素值高于舊址,反之,代表新址的要素值低于舊址。

表1 蕭山站新址與舊址各要素月(年)差值平均值表

表2 蕭山站新址與舊址各要素月(年)差值標(biāo)準(zhǔn)差σ

2.2 氣溫差異分析

經(jīng)過統(tǒng)計(jì)發(fā)現(xiàn)新、舊站址的氣溫存在著較大的差異。由表1可看出,新、舊站址平均氣溫的差值年平均為-0.4 ℃,即新址的年平均氣溫比舊址低0.4 ℃。1—12月的月平均氣溫偏差在-0.2~-0.5 ℃之間,其中1月、3月和8—12月差值為-0.4~-0.5 ℃,2月和4—7月差值為-0.2~-0.3 ℃,7月差值最小為-0.2 ℃,這說明新、舊址的月平均氣溫在春、夏季節(jié)偏差較小,秋、冬季節(jié)偏差較大。

新、舊址最高氣溫的年平均差值為-0.8 ℃,即新址最高氣溫的年平均值比舊址低0.8 ℃,相比于平均氣溫和最低氣溫,新舊址的最高氣溫相差較大。新、舊址最高氣溫的月平均差值在-0.4~-1.1 ℃之間,其中1—3月、8—12月的差值在-0.7~-1.1 ℃之間,其余月份差值在-0.4~-0.6 ℃之間,這同樣說明了新、舊址的最高氣溫在春、夏季節(jié)偏差小,秋、冬季節(jié)偏差大。

最低氣溫差值年平均為-0.2 ℃,相比于平均氣溫和最高氣溫,新舊址的最低氣溫相差較小。新、舊址最低氣溫的月平均差值在-0.1~-0.4 ℃之間,其中1月份新、舊址最低氣溫差值最大,為-0.4 ℃,其次為9、10、11月,差值均為-0.3 ℃。這也說明了新、舊址最低氣溫在春、夏季節(jié)偏差小,秋、冬季節(jié)偏差大的特征,但新、舊址最低氣溫偏差的季節(jié)特征不如平均氣溫和最高氣溫明顯。

由表2月(年)氣溫差值標(biāo)準(zhǔn)差表明,平均氣溫、最高氣溫、最低氣溫?cái)?shù)據(jù)離散程度較小,數(shù)據(jù)波動(dòng)小比較穩(wěn)定。經(jīng)過對(duì)比,2016年新、舊兩站的極端氣溫出現(xiàn)時(shí)間完全吻合,其中極端最高氣溫均出現(xiàn)在7月26日,極端最低氣溫均出現(xiàn)在1月25日;分析新、舊站平均氣溫隨時(shí)間的變化趨勢(shì)(圖1),兩者同樣具有極高的吻合度。

圖1 2016年新、舊址平均氣溫的逐月分布圖

2.3 相對(duì)濕度差異分析

由表1可見,新址全年各月相對(duì)濕度均大于舊址,新址年平均相對(duì)濕度比舊址高3%。2月、7月新、舊址的相對(duì)濕度月差值相對(duì)較小,9—12月是相對(duì)濕度差值較大的月份,其中10月份新舊址相對(duì)濕度差值為6%,為一年中相對(duì)濕度差值最大的月份。由表2可看出,月(年)相對(duì)濕度差值標(biāo)準(zhǔn)差較大,數(shù)據(jù)波動(dòng)較大。但兩站相對(duì)濕度在時(shí)間序列分布上呈現(xiàn)出高度一致性(圖2)。

圖2 2016年新、舊址平均相對(duì)濕度的逐月分布圖

2.4 降水量差異分析

2.4.1 降水量差值分析

從表1可看出,2016年新址的年降水量比舊址偏多110.3 mm,除2月份的降水量月差值為負(fù)值外,其余各月均為正值,這說明了只有2月份降水量新址少于舊址,其余各月份降水量新址均多于舊址。新舊址在4—6月和9—11月的降水量差異較大,差值在8.6~31.4 mm之間,1—3月、7—8月和12月差異相對(duì)較小,在-3.4~4.9 mm之間。3月份的降水量差值最小,為0.7 mm,代表3月新、舊址的降水量最接近。9月份差值為31.4 mm,是一年中降水量差值最大的月份。通過統(tǒng)計(jì),2016年新址的雨日比舊址偏少22 d。

由表2可看出,降水差值標(biāo)準(zhǔn)差較氣溫大,數(shù)據(jù)離散程度相對(duì)不穩(wěn)定,波動(dòng)較大。

2.4.2 降水量累計(jì)相對(duì)差值分析

計(jì)算2016年的逐月降水量和年降水量的累計(jì)相對(duì)差值,發(fā)現(xiàn)新、舊址的年降水量累計(jì)相對(duì)差值為7%,說明新舊站址間降水量的總量差異較大。由圖3可知:除2月份新址降水量少于舊址外,其余各月均大于舊址,3月、7月差異最小為1%,9月最大為16%,全年只有5個(gè)月(1月、3月、6—8月)累計(jì)相對(duì)差值小于等于4%。

圖3 2016年新、舊址的降水量累計(jì)相對(duì)差值的逐月分布圖

2.5 風(fēng)速、風(fēng)向差異分析

2.5.1 風(fēng)速差異分析

由表1可知,新、舊站址平均風(fēng)速差值年平均為2.0 m/s,月平均差值在1.3~2.7 m/s之間,其中4—7月的平均風(fēng)速差值在1.3~1.9 m/s之間,1—3月、8—12月在2.0~2.7 m/s之間,這說明2016年1—12月新址的平均風(fēng)速均大于舊址,新、舊址在春、夏季的風(fēng)速相差較小,秋冬季相差大。通過計(jì)算平均風(fēng)速日差值,發(fā)現(xiàn)平均風(fēng)速日差值最大為7.0 m/s、最小-0.2 m/s,2016年期間新址只有兩天的平均風(fēng)速較舊址偏小,其余均偏大,這說明新址的風(fēng)速基本上大于舊址的風(fēng)速,只有非常少數(shù)的情況下舊址的風(fēng)速大于新址,這與新舊址的地理位置有關(guān),舊址北側(cè)緊鄰山坡,西西北到東北方向有山脈遮擋,新址在山脈頂部的平地上。風(fēng)向?yàn)槲魑鞅钡綎|北方向時(shí),舊址受山脈阻擋,因此風(fēng)速會(huì)比新址小。

通過表2可以看出,平均風(fēng)速差值標(biāo)準(zhǔn)差較大,數(shù)據(jù)離散程度相對(duì)不穩(wěn)定,波動(dòng)較大。2016年最大風(fēng)速差值的年平均為3.6 m/s、月平均在3.0~4.6 m/s之間,日差值最大10.4 m/s、最小-0.7 m/s,新址全年只有一天較舊址偏小,其余均偏大;2016年極大風(fēng)速差值年平均為3.5 m/s、月平均在2.8~4.4 m/s之間、日差值最大23.4 m/s、最小-0.8 m/s,新址全年有5 d較舊址偏小,其余均偏大。由圖4可見,新、舊站址風(fēng)速差異較大,新址平均風(fēng)速、最大風(fēng)速、極大風(fēng)速差值平均都明顯大于舊址,平均風(fēng)速差值、最大風(fēng)速差值和極大風(fēng)速差值的最大值均出現(xiàn)在10月,并且發(fā)現(xiàn)風(fēng)速越大差值越大,當(dāng)出現(xiàn)大風(fēng)天氣現(xiàn)象時(shí),差異尤為顯著;新址全年出現(xiàn)大風(fēng)日數(shù)為15 d,而舊址全年未出現(xiàn)大風(fēng)天氣。新址全年靜風(fēng)出現(xiàn)回?cái)?shù)為8,而舊址為112。

圖4 新、舊址風(fēng)速差值的逐月分布圖

2.5.2 風(fēng)向差異分析

表3是新、舊址的最多風(fēng)向頻率風(fēng)向相符率統(tǒng)計(jì)表,可以看出新、舊址在1—12月和年最多風(fēng)向及頻率均相差較大。1月、8月和10月新、舊址的最多風(fēng)向比較接近,相差一個(gè)方位,4月最多風(fēng)向差異最大,相差8個(gè)方位。新址年最多風(fēng)向N,舊址年最多風(fēng)向WSW,相差5個(gè)方位。據(jù)統(tǒng)計(jì)年風(fēng)向頻率除WNW相同外,其余各風(fēng)向及靜風(fēng)都存在差異,其中N、WSW和NNW差別較大(圖略)。新址與舊址年風(fēng)向相符率只有42.5%,最大相符率在8月為53.2%,最小相符率在2月為36.5%,可見兩站風(fēng)向一致性較差。

2.6 顯著性檢驗(yàn)

為了進(jìn)一步了解遷站對(duì)觀測(cè)數(shù)據(jù)連續(xù)性的影響,對(duì)新址2016年和舊址近20 a(1996—2015年)的觀測(cè)資料進(jìn)行均值差異T檢驗(yàn),置信區(qū)間的置信度為95%,顯著性水平0.05(顯著性水平為0.05的標(biāo)準(zhǔn)值為2.086)(表4)。

表3 新、舊站址的最多風(fēng)向頻率風(fēng)向相符率統(tǒng)計(jì)表

表4 新、舊址各氣象要素的T檢驗(yàn)

結(jié)果顯示年平均氣溫、年降水量、年平均相對(duì)濕度,通過了T檢驗(yàn),年平均風(fēng)速未通過,表明平均氣溫、降水量和平均相對(duì)濕度連續(xù)性較好,平均風(fēng)速連續(xù)性相對(duì)較差。

由月平均值檢驗(yàn)結(jié)果可知:12月平均氣溫?cái)?shù)據(jù)有顯著性差異,其它月份數(shù)據(jù)顯著性差異較小,連續(xù)性較好;4月、5月的降水量差異較顯著,其次是9月和1月,其余月份的降水量差異不顯著;平均相對(duì)濕度除2月、10月、11月顯著外,其它月份的差異均不顯著;平均風(fēng)速除7月不顯著外,其他月份的差異均較顯著,年平均風(fēng)速的差異也較顯著。

3 結(jié)果分析

利用蕭山站站址遷移對(duì)比期(2016年1月1日至12月31日)新、舊址的氣溫、降水、相對(duì)濕度、風(fēng)向、風(fēng)速等氣象要素的逐日觀測(cè)數(shù)據(jù),采用差值標(biāo)準(zhǔn)差、降水累計(jì)相對(duì)差值、風(fēng)向相符率、顯著性檢驗(yàn)等統(tǒng)計(jì)分析方法,對(duì)以上要素進(jìn)行對(duì)比分析,結(jié)果表明:

1)受海拔高度和周圍地理環(huán)境差異的影響,2016年1—12月新址的平均氣溫、最高氣溫、最低氣溫均低于舊址,其年平均差值分別為-0.4 ℃、-0.8 ℃、-0.2 ℃,最高氣溫差異最大。但各類氣溫的差異具有一致的季節(jié)變化特征,均在春、夏季節(jié)比較接近,秋、冬季節(jié)相差較大。新、舊址的極端最高氣溫和最低氣溫的出現(xiàn)日期完全一致,平均氣溫隨時(shí)間變化趨勢(shì)吻合。

2)2016年1—12月新址的月平均相對(duì)濕度均大于舊址,新址年平均相對(duì)濕度比舊址高3%,其中9—12月相差較大。但兩者隨時(shí)間的變化趨勢(shì)基本一致。

3)2016年新址的年降水量比舊址偏多110.3 mm,但雨日比舊址偏少22 d。新、舊址年降水量的累計(jì)相對(duì)差值為7%。從降水量的差值、差值標(biāo)準(zhǔn)差及累計(jì)相對(duì)差值來看,新舊址的降水量觀測(cè)數(shù)據(jù)差異較大,主要體現(xiàn)在4—6月和9—11月,而1—3月、7—8月和12月差異相對(duì)較小。從日、月降水分析認(rèn)為,導(dǎo)致降水差異大的原因:一是地理位置不同,舊址海拔低,春夏時(shí)節(jié),副高北抬,南來暖濕空氣受山脈阻擋抬升,溫度下降幅度小,水汽凝結(jié)少,降水量相對(duì)偏少;而新址海拔高,暖濕空氣在被迫抬升過程中水汽凝結(jié)量多于舊址,造成新站降水多于舊站;秋冬多西北氣流,加之副高南撤,新站海拔高,曖濕空氣爬山冷卻凝結(jié)致山頂多雨,而舊站位于新站東南方向,海拔低,兩站水平相距只有198 m,受氣流越山后下沉,云滴蒸發(fā)加大影響,造成舊站降水偏少。二是地理環(huán)境不同,新舊站址所處地理環(huán)境不同,風(fēng)向和風(fēng)速出現(xiàn)了較大差異,因而對(duì)水汽輸送形成影響;三是季節(jié)影響,4—6月、9—11月局地對(duì)流性降水較多,反應(yīng)了降水的局地性特征。

4)2016年新址的平均風(fēng)速、最大風(fēng)速、極大風(fēng)速均比舊址偏大,年平均差值分別為2.0、3.6、3.5 m/s,其中春、夏季風(fēng)速相差較小,秋、冬季相差大,差值最大值均出現(xiàn)在10月。新站大風(fēng)日數(shù)多于舊站、靜風(fēng)出現(xiàn)次數(shù)少于舊站。兩站風(fēng)向出現(xiàn)的一致性較差,2016年的風(fēng)向相符率僅為42.5%,新址年最多風(fēng)向?yàn)楸憋L(fēng)(N),舊址年最多風(fēng)向是西西南(WSW),相差5個(gè)方位。主要原因是舊址的西西北到東北方向有連續(xù)障礙物阻擋,形成了“避風(fēng)港”;而新址地勢(shì)開闊,氣象觀測(cè)環(huán)境良好,全年各月風(fēng)向出現(xiàn)頻率符合季風(fēng)氣候特征,避免了舊站受山地屏障的影響,測(cè)得的風(fēng)速風(fēng)向更能代表本地氣候特征。

5)顯著性檢驗(yàn)結(jié)果表明,新、舊址的平均氣溫、平均相對(duì)濕度和降水量與舊址近20 a的觀測(cè)數(shù)據(jù)差異不顯著,各要素隨時(shí)間變化趨勢(shì)的一致性較好,根據(jù)其相對(duì)固定的數(shù)據(jù)差異,可在氣候分析和應(yīng)用中延續(xù)使用;受特殊地理環(huán)境影響,新、舊站風(fēng)要素差別較大,延續(xù)性使用價(jià)值不大,可根據(jù)不同的應(yīng)用需求,分別統(tǒng)計(jì)分析使用,新站風(fēng)要素更具有代表性。

參考文獻(xiàn):

[1] 中國(guó)氣象局.地面氣象觀測(cè)規(guī)范[M].北京:氣象出版社,2003:8.

[2] 宋莎莎.東平新舊站址氣象要素對(duì)比評(píng)估[B].山東氣象,2014,34(2):81-84.

[3] 黃愛榮,張莉.阿克蘇國(guó)家基準(zhǔn)站遷站資料對(duì)比分析[J].沙漠與綠洲氣象,2016,10(增刊):99-103.

[4] 周建平,孫照渤,倪東鴻,等.中國(guó)氣象臺(tái)站遷移對(duì)年平均氣溫均一性的影響[J].大氣科學(xué)學(xué)報(bào),2013,36(2):139-146.

[5] 吳利紅,駱月珍,孫莉莉,等.浙江省近34年年平均風(fēng)速序列均一性檢驗(yàn)研究[J].氣象科技,2008,36(5):661-665.

[6] 吳利紅,康麗莉,陳海燕,等.地面氣象站環(huán)境變化對(duì)氣溫序列均一性影響[J].氣象科技,2007,35(1):152-156.

[7] 胡蓉,楊長(zhǎng)登,鄧芩.玉屏新舊站址氣象要素對(duì)比分析[J].貴州氣象,2015,39(2):40-44.

[8] 中國(guó)氣象局.遷站規(guī)范[S].氣發(fā)〔2012〕93號(hào)文,2012:11.

猜你喜歡
風(fēng)速差異
相似與差異
音樂探索(2022年2期)2022-05-30 21:01:37
基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測(cè)
基于最優(yōu)TS評(píng)分和頻率匹配的江蘇近海風(fēng)速訂正
找句子差異
DL/T 868—2014與NB/T 47014—2011主要差異比較與分析
生物為什么會(huì)有差異?
基于GARCH的短時(shí)風(fēng)速預(yù)測(cè)方法
考慮風(fēng)切和塔影效應(yīng)的風(fēng)力機(jī)風(fēng)速模型
GE在中國(guó)發(fā)布2.3-116低風(fēng)速智能風(fēng)機(jī)
考慮風(fēng)速分布與日非平穩(wěn)性的風(fēng)速數(shù)據(jù)預(yù)處理方法研究
主站蜘蛛池模板: 精品一区二区久久久久网站| 日韩一区二区三免费高清 | 婷婷开心中文字幕| 亚洲天堂网在线视频| 亚洲天堂日韩av电影| 亚洲精品少妇熟女| 日韩成人在线视频| 性色一区| 婷婷亚洲视频| 在线无码九区| 精品国产aⅴ一区二区三区| 波多野结衣无码视频在线观看| 日韩一二三区视频精品| 国产成人福利在线视老湿机| 成人在线第一页| 亚洲一欧洲中文字幕在线| 91高清在线视频| 伊人久综合| 午夜无码一区二区三区在线app| 国产香蕉一区二区在线网站| 日本免费精品| 欧美精品高清| 久久这里只有精品国产99| 国产H片无码不卡在线视频| 国产成人高清在线精品| 亚洲三级片在线看| 欧美一区福利| 米奇精品一区二区三区| 亚洲三级视频在线观看| 欧美综合区自拍亚洲综合绿色 | 污视频日本| 国产精品美乳| 激情综合激情| 亚洲αv毛片| 国产视频一区二区在线观看| 一区二区在线视频免费观看| 婷婷色在线视频| 黄色国产在线| 99re在线视频观看| 国产性精品| 国产91特黄特色A级毛片| 人妻无码一区二区视频| 欧美一级专区免费大片| 中文字幕日韩视频欧美一区| 国产精品综合久久久| 久久香蕉欧美精品| 国产不卡在线看| 亚洲欧美国产高清va在线播放| 九九九精品视频| 国产精品一区在线观看你懂的| 日韩中文欧美| 麻豆国产原创视频在线播放| 亚洲人成网18禁| 欧美在线天堂| 欧美专区日韩专区| 精品视频在线观看你懂的一区| 午夜国产精品视频黄| 亚洲一区精品视频在线 | 福利小视频在线播放| 亚洲色婷婷一区二区| 极品国产在线| 亚洲精品中文字幕午夜| 在线无码九区| 国产无码性爱一区二区三区| 亚洲天堂自拍| 国产亚洲欧美日本一二三本道| 中文字幕人妻av一区二区| 国产午夜一级毛片| 午夜不卡视频| 亚洲人成亚洲精品| 91区国产福利在线观看午夜| 丁香婷婷久久| 亚洲av无码人妻| 欧美黄网在线| 精品一区二区久久久久网站| 色成人综合| 无码一区二区波多野结衣播放搜索| 国产凹凸视频在线观看| 亚洲91在线精品| 日韩av无码DVD| 精品欧美视频| 亚洲国产天堂久久综合|