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

黔西市各鄉(xiāng)鎮(zhèn)精細(xì)化氣溫預(yù)報(bào)訂正研究

2024-12-31 00:00:00姜軍志呂靜蔡軍楊鑫龍張玉逵
農(nóng)業(yè)災(zāi)害研究 2024年10期

摘 要:選取2020—2022年黔西市國(guó)家基本氣象站和27個(gè)鄉(xiāng)鎮(zhèn)區(qū)域自動(dòng)站日最高氣溫和日最低氣溫資料,利用最小二乘法得出各鄉(xiāng)鎮(zhèn)觀測(cè)站點(diǎn)與黔西國(guó)家基本站氣溫關(guān)系式,建立預(yù)報(bào)方法,采取2023年實(shí)況數(shù)據(jù)分別對(duì)方程組預(yù)報(bào)效果和經(jīng)驗(yàn)訂正預(yù)報(bào)效果進(jìn)行檢驗(yàn)。結(jié)果顯示,精細(xì)化氣溫預(yù)報(bào)訂正后的預(yù)報(bào)效果總體較好,對(duì)黔西市各鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)有較好的指導(dǎo)意義。

關(guān)鍵詞:氣溫預(yù)報(bào);預(yù)報(bào)效果;最小二乘法

中圖分類號(hào):P45 文獻(xiàn)標(biāo)志碼:B 文章編號(hào):2095–3305(2024)10–0-03

隨著我國(guó)經(jīng)濟(jì)社會(huì)的發(fā)展,社會(huì)各界對(duì)天氣預(yù)報(bào)的準(zhǔn)確率和精細(xì)化程度的要求越來(lái)越高。隨之而來(lái)對(duì)氣象預(yù)報(bào)精細(xì)化的需求日益增長(zhǎng),其中鄉(xiāng)鎮(zhèn)氣溫的預(yù)報(bào)是重點(diǎn)之一。以貴州省畢節(jié)市黔西市為例,以國(guó)家站為基準(zhǔn)的城區(qū)氣溫預(yù)報(bào)準(zhǔn)確率較高,但在開(kāi)展鄉(xiāng)鎮(zhèn)預(yù)報(bào)時(shí),通常以基本站為基準(zhǔn),采用經(jīng)驗(yàn)訂正的方式開(kāi)展氣象服務(wù),此方法無(wú)法滿足當(dāng)前社會(huì)需求。近年來(lái),鄉(xiāng)鎮(zhèn)區(qū)域站網(wǎng)布局越發(fā)完善,氣象數(shù)據(jù)更加全面準(zhǔn)確,為更好地開(kāi)展鄉(xiāng)鎮(zhèn)氣象服務(wù),基于陳云飛等[1]、黃治勇等[2]、邱學(xué)興等[3]對(duì)鄉(xiāng)鎮(zhèn)氣溫訂正預(yù)報(bào)方法的研究,采用SPSS軟件,對(duì)2020—2022年黔西國(guó)家站和各鄉(xiāng)鎮(zhèn)區(qū)域站日最高氣溫和日最低氣溫相關(guān)數(shù)據(jù)進(jìn)行分析和回歸計(jì)算,建立各鄉(xiāng)鎮(zhèn)區(qū)域站與國(guó)家基本氣象站相關(guān)聯(lián)的氣溫預(yù)報(bào)回歸方程,形成黔西市各鄉(xiāng)鎮(zhèn)氣溫精細(xì)化預(yù)報(bào)的訂正客觀方法,以提高黔西市各鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)準(zhǔn)確率,為農(nóng)業(yè)生產(chǎn)和人民生活提供幫助。

1 資料和方法

2020—2022年黔西市27個(gè)鄉(xiāng)鎮(zhèn)區(qū)域自動(dòng)站資料,

來(lái)源于貴州省氣象局提供的區(qū)域自動(dòng)站數(shù)據(jù)文件,黔西國(guó)家基本站資料來(lái)源于本站月報(bào)表數(shù)據(jù)文件。以下使用的氣象資料在計(jì)算前進(jìn)行了嚴(yán)格的質(zhì)量控制及均一性檢查,同時(shí)對(duì)缺測(cè)資料進(jìn)行了插補(bǔ)訂正。

選取2020—2022年黔西市27個(gè)鄉(xiāng)鎮(zhèn)區(qū)域自動(dòng)站和黔西國(guó)家基本氣象站的日最高氣溫及日最低氣溫資料,并通過(guò)統(tǒng)計(jì)分析各鄉(xiāng)鎮(zhèn)及黔西國(guó)家站的年平均溫度、年平均最高氣溫、年平均最低氣溫、月平均氣溫、月平均最高氣溫、月平均最低氣溫等氣象指標(biāo)的分布區(qū)域情況及變化趨勢(shì),發(fā)現(xiàn)鄉(xiāng)鎮(zhèn)氣溫與黔西國(guó)家站氣溫之間存在線性關(guān)系,并運(yùn)用最小二乘法擬合各鄉(xiāng)鎮(zhèn)區(qū)域自動(dòng)站與黔西國(guó)家站的最佳趨勢(shì)預(yù)報(bào)方程,形成不同鄉(xiāng)鎮(zhèn)的氣溫回歸模型,從數(shù)理統(tǒng)計(jì)的角度上客觀訂正各鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)。

根據(jù)中國(guó)氣象局于2005年7月1日下發(fā)的《中短期天氣預(yù)報(bào)質(zhì)量檢驗(yàn)辦法(試行)》,溫度預(yù)報(bào)檢驗(yàn)內(nèi)容包括最高氣溫和最低氣溫的預(yù)報(bào)誤差[4-5]。主要有以下3種公式:

平均絕對(duì)誤差:TMAE=|Fi-Oi|(1)

均方根誤差:TRMSE=(2)

預(yù)報(bào)準(zhǔn)確率:TTK=×100%(3)

在式(1)~式(3)中,F(xiàn)i為第i站(次)預(yù)報(bào)溫度,Oi為第i站(次)實(shí)況溫度。K值為預(yù)報(bào)準(zhǔn)確率檢驗(yàn)誤差值,K為1代表|Fi-Oi|≤1 ℃;K為2代表|Fi-Oi|≤2 ℃。NrK為預(yù)報(bào)正確的站(次)數(shù),NfK為預(yù)報(bào)的總站(次)數(shù)。

2 結(jié)果與分析

2.1 氣溫預(yù)報(bào)的訂正方法

利用2020—2022年共3年數(shù)據(jù)作為預(yù)報(bào)方程建立樣本。以黔西國(guó)家基本站為代表,采用最小二乘法擬合各鄉(xiāng)鎮(zhèn)區(qū)域站點(diǎn)與黔西站的每日最高氣溫和最低氣溫的關(guān)系式。設(shè)各鄉(xiāng)鎮(zhèn)的氣溫為因變量( y),黔西站的氣溫為自變量(x),建立了各鄉(xiāng)鎮(zhèn)區(qū)域站和黔西國(guó)家站的趨勢(shì)預(yù)報(bào)方程,所有趨勢(shì)預(yù)報(bào)方程的復(fù)相關(guān)系數(shù)均在0.970以上,回歸效果明顯。在實(shí)際運(yùn)用中發(fā)現(xiàn),使用趨勢(shì)方程對(duì)各鄉(xiāng)鎮(zhèn)氣溫訂正預(yù)報(bào)效果良好,這表明此方法建立的趨勢(shì)方程在預(yù)報(bào)鄉(xiāng)鎮(zhèn)氣溫時(shí)具有良好的實(shí)用性,對(duì)提高鄉(xiāng)鎮(zhèn)最低氣溫、最高氣溫預(yù)報(bào)質(zhì)量具有重大意義。各鄉(xiāng)鎮(zhèn)與黔西站的趨勢(shì)預(yù)報(bào)方程具體見(jiàn)表1。

2.2 鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)的效果檢驗(yàn)

通過(guò)2023年黔西市各鄉(xiāng)鎮(zhèn)逐日氣溫預(yù)報(bào)資料與同期實(shí)況氣溫進(jìn)行對(duì)比評(píng)估,選取的有效樣本數(shù)共計(jì)365 d,根據(jù)《中短期天氣預(yù)報(bào)質(zhì)量檢驗(yàn)辦法(試行)》的評(píng)分方法,統(tǒng)計(jì)得出黔西市各鄉(xiāng)鎮(zhèn)2023年每日氣溫預(yù)報(bào)的平均絕對(duì)方差、均方根誤差和預(yù)報(bào)準(zhǔn)確率。

從表2可以看出,黔西市各鄉(xiāng)鎮(zhèn)最高氣溫和最低氣溫的平均絕對(duì)誤差均≤3.1 ℃,均方根誤差≤4.1 ℃,

其中最高氣溫平均絕對(duì)誤差≤3 ℃的預(yù)報(bào)準(zhǔn)確率高于90%,最低氣溫平均絕對(duì)誤差≤2℃的預(yù)報(bào)準(zhǔn)確率高于60%,這一結(jié)果表明黔西市各鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)總體較為準(zhǔn)確。通過(guò)對(duì)比可發(fā)現(xiàn),最低氣溫預(yù)報(bào)準(zhǔn)確率明顯高于最高氣溫預(yù)報(bào)準(zhǔn)確率,最低氣溫的平均絕對(duì)誤差、均方根誤差等普遍低于最高氣溫,最高氣溫出現(xiàn)大值預(yù)報(bào)偏差的概率普遍高于最低氣溫。

一般情況下,當(dāng)實(shí)況氣溫與預(yù)報(bào)氣溫相近時(shí),公眾會(huì)認(rèn)為預(yù)報(bào)是準(zhǔn)確可信的。但當(dāng)實(shí)況氣溫與預(yù)報(bào)氣溫存在差異,且超過(guò)一定閾值時(shí),易損害氣象部門(mén)公信力[6-8]。根據(jù)天氣預(yù)報(bào)質(zhì)量評(píng)定相關(guān)辦法,溫度預(yù)報(bào)絕對(duì)誤差>4.5 ℃,即溫度預(yù)報(bào)評(píng)分60分作為大值的閾值[9-11]。從表2中可以發(fā)現(xiàn),黔西市氣溫預(yù)報(bào)≤4.5 ℃的預(yù)報(bào)準(zhǔn)確率高于75%,最高可達(dá)94.8%,表明黔西市各鄉(xiāng)鎮(zhèn)的氣溫趨勢(shì)預(yù)報(bào)方程可用性較高。

通過(guò)對(duì)2023年全年的各鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)效果檢驗(yàn)可以看出,此趨勢(shì)預(yù)報(bào)方程的建立,對(duì)黔西市各鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)具有一定的指導(dǎo)作用,對(duì)黔西市開(kāi)展精細(xì)化氣象服務(wù),提高預(yù)報(bào)服務(wù)質(zhì)量具有重要意義。

表1" 各鄉(xiāng)鎮(zhèn)最高、最低氣溫預(yù)報(bào)方程

鄉(xiāng)鎮(zhèn) 日最高氣溫預(yù)報(bào)方程 復(fù)相關(guān)系數(shù) 日最低氣溫預(yù)報(bào)方程 復(fù)相關(guān)系數(shù)

甘棠 y=1.040 7x+0.751 8 0.970 y=0.971 5x+0.527 2 0.969

大關(guān) y=0.997 1x+0.725 8 0.990 y=0.983 3x+0.726 6 0.993

定新 y=1.045 3x+0.562 9 0.972 y=0.960 6x+1.189 2 0.981

谷里 y=1.017 1x+0.530 5 0.991 y=1.002 7x+0.748 3 0.993

觀音洞 y=0.989 8x-0.420 1 0.986 y=0.993 9x-0.207 2 0.987

紅林 y=1.031 5x-0.787 6 0.984 y=0.939x-0.448 6 0.979

花溪 y=1.047x+0.048 4 0.973 y=0.978 2x+0.687 8 0.988

金碧 y=1.007 4x+0.599 1 0.991 y=0.985x+0.249 3 0.994

金蘭 y=0.990 7x+0.211 3 0.985 y=0.977x-0.001 1 0.988

金坡 y=1.001 4x-0.833 9 0.978 y=0.941 5x-0.704 7 0.964

錦星 y=1.028 7x+1.1 0.990 y=0.969 7x+0.678 0.974

林泉 y=1.006 2x-0.3364 0.992 y=0.965x-0.224 7 0.987

綠化 y=0.983 5x-0.05 0.992 y=0.966x+0.042 4 0.990

仁和 y=1.041 9x+0.063 9 0.985 y=0.963 4x+0.341 0.977

素樸 y=1.023 1x-0.040 2 0.986 y=1.004x-0.232 1 0.973

太來(lái) y=1.042 9x-0.345 8 0.983 y=1.002 2x+0.156 0.991

鐵石 y=1.014 4x+0.751 9 0.987 y=1.009 9x+0.663 7 0.987

五里 y=0.986 3x-0.414 6 0.985 y=0.989 9x-0.277 4 0.994

協(xié)和 y=1.014 8x+0.795 5 0.971 y=0.979 9x-0.150 5 0.971

新仁 y=0.988 7x+0.050 9 0.985 y=0.965 1x+0.591 4 0.979

永興 y=1.033 6x+0.597 6 0.988 y=0.974 3x+0.841 7 0.978

雨朵 y=0.991x+0.701 8 0.991 y=0.976 9x+0.566 1 0.987

長(zhǎng)堰 y=0.988 8x+0.460 8 0.993 y=0.940 7x+0.502 2 0.974

中建 y=1.039 9x+0.983 5 0.971 y=0.965 3x+0.879 8 0.968

中坪 y=1.020 4x+1.945 4 0.978 y=0.952 6x+2.104 5 0.971

鐘山 y=1.015 6x+0.441 1 0.989 y=0.981 8x-0.096 7 0.970

重新 y=1.006 3x-0.220 4 0.978 y=1.014 1x-0.562 9 0.967

表2" 黔西市各鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)評(píng)分

鄉(xiāng)鎮(zhèn) 平均絕對(duì)誤差/℃ 均方根誤差/℃ ≤2 ℃預(yù)報(bào)準(zhǔn)確率/% ≤4.5 ℃預(yù)報(bào)準(zhǔn)確率/%

最高氣溫 最低氣溫 最高氣溫 最低氣溫 最高氣溫 最低氣溫 最高氣溫 最低氣溫

甘棠 3.1 2.0 4.1 2.8 42.2 64.4 75.3 87.9

大關(guān) 2.8 1.8 3.7 3.9 46.3 66.3 78.6 92.3

定新 3.0 1.9 3.9 3.7 41.4 64.1 76.7 90.7

谷里 2.8 1.8 3.7 3.7 47.4 66.0 78.9 92.1

觀音洞 2.8 1.5 3.7 3.9 50.7 72.3 77.5 94.8

紅林 2.9 1.8 3.9 3.8 47.9 65.8 75.9 92.1

鄉(xiāng)鎮(zhèn) 平均絕對(duì)誤差/℃ 均方根誤差/℃ ≤2 ℃預(yù)報(bào)準(zhǔn)確率/% ≤4.5 ℃預(yù)報(bào)準(zhǔn)確率/%

最高氣溫 最低氣溫 最高氣溫 最低氣溫 最高氣溫 最低氣溫 最高氣溫 最低氣溫

花溪 2.9 1.8 3.8 3.7 44.4 66.3 79.2 92.3

金碧 2.8 1.8 3.7 3.7 50.1 66.0 79.5 92.1

金蘭 2.7 2.0 3.7 3.6 50.4 62.5 78.9 90.1

金坡 2.7 2.0 3.6 3.8 47.9 63.0 80.3 91.5

錦星 2.8 2.1 3.8 3.7 47.4 61.1 80.3 88.5

林泉 2.8 1.8 3.7 3.6 49.0 66.0 77.8 92.3

綠化 2.7 1.8 3.6 3.9 48.2 67.7 78.6 93.2

仁和 2.9 2.0 3.9 3.8 47.9 61.4 77.3 90.1

素樸 2.8 2.1 3.8 3.9 47.7 62.2 79.5 87.7

太來(lái) 2.8 1.8 3.9 3.8 48.8 66.6 80.0 92.3

鐵石 2.8 1.9 3.8 3.6 46.0 64.1 80.0 91.2

五里 2.7 1.7 3.6 3.8 50.4 67.7 79.2 94.0

協(xié)和 2.9 1.7 3.8 3.6 46.0 68.5 78.9 93.4

新仁 2.8 1.9 3.6 3.8 46.6 64.9 79.5 92.1

永興 2.8 2.0 3.8 3.6 48.5 63.6 78.1 90.1

雨朵 2.7 1.9 3.6 3.6 49.6 64.9 78.1 91.2

長(zhǎng)堰 2.7 2.0 3.6 3.9 49.9 63.6 79.2 90.1

中建 2.9 2.1 3.9 3.7 45.8 61.4 77.3 88.5

中坪 2.8 2.1 3.7 3.7 47.1 62.5 79.7 89.9

鐘山 2.7 2.1 3.7 3.8 49.9 62.5 79.7 86.8

重新 2.8 2.3 3.8 3.1 50.1 59.7 79.5 85.8

3 結(jié)論

(1)黔西市各鄉(xiāng)鎮(zhèn)區(qū)域自動(dòng)站與國(guó)家站氣溫的差值絕對(duì)值最大值均>2.0 ℃,由此可見(jiàn),各鄉(xiāng)鎮(zhèn)區(qū)域自動(dòng)站與國(guó)家站氣溫存在差異,單純用黔西國(guó)家站的氣溫預(yù)報(bào)代表鄉(xiāng)鎮(zhèn)溫度預(yù)報(bào),其準(zhǔn)確率較低。

(2)各鄉(xiāng)鎮(zhèn)最低氣溫預(yù)報(bào)準(zhǔn)確率明顯高于最高氣溫預(yù)報(bào)準(zhǔn)確率,意味著最高氣溫出現(xiàn)大值偏差的概率較大。通過(guò)對(duì)預(yù)報(bào)效果檢驗(yàn)可以發(fā)現(xiàn),各鄉(xiāng)鎮(zhèn)最高氣溫預(yù)報(bào)的絕對(duì)誤差≤2 ℃的預(yù)報(bào)準(zhǔn)確率為41.1%

~50.7%,最低氣溫預(yù)報(bào)的絕對(duì)誤差≤2 ℃的預(yù)報(bào)準(zhǔn)確率為59.7%~66.6%。

(3)因氣溫預(yù)報(bào)存在大值偏差,故采用溫度預(yù)報(bào)絕對(duì)誤差>4.5 ℃為大值預(yù)報(bào)誤差閾值。黔西市各鄉(xiāng)鎮(zhèn)最高氣溫絕對(duì)值誤差≤4.5 ℃的預(yù)報(bào)準(zhǔn)確率為75.3%~80.3%,最低氣溫絕對(duì)值誤差≤4.5 ℃的預(yù)報(bào)準(zhǔn)確率為85.8%~94.0%。

參考文獻(xiàn)

[1] 陳云飛,祁雁文,韓聰,等.基于城鎮(zhèn)預(yù)報(bào)的鄉(xiāng)鎮(zhèn)氣溫預(yù)報(bào)方法初探[J].內(nèi)蒙古氣象,2016(5):35-38.

[2] 黃治勇,張文,陳璇,等.湖北省鄉(xiāng)鎮(zhèn)溫度預(yù)報(bào)方法初探[J].氣象,2011,37(12):1578-1583.

[3] 邱學(xué)興,王東勇,朱紅芳.鄉(xiāng)鎮(zhèn)精細(xì)化最高最低氣溫預(yù)報(bào)方法研究[J].氣象與環(huán)境學(xué)報(bào),2013,29(3):92-96.

[4] 蔡軍,楊承梅,倪妮,等.威寧縣鄉(xiāng)鎮(zhèn)最高最低氣溫分析與預(yù)報(bào)方法探討[J].中低緯山地氣象,2019,43(5):29-35.

[5] 郝占宇.呼倫貝爾市氣候分區(qū)氣溫精細(xì)化訂正預(yù)報(bào)方法研究[J].農(nóng)業(yè)災(zāi)害研究,2024,14(1):212-214.

[6] 劉嘉慧敏,潘留杰,戴昌明,等.我國(guó)氣溫網(wǎng)格預(yù)報(bào)檢驗(yàn)及客觀訂正方法研究進(jìn)展[J].氣象科技進(jìn)展,2023,13(6):10-20.

[7] 李天江,楊曉玲,張占文,等.河西走廊東部智能網(wǎng)格溫度預(yù)報(bào)產(chǎn)品訂正方法探討[J].干旱氣象,2023,41(5):802-810.

[8] 羅晨藝,劉紅雙,張?jiān)魄铮?冬季黔北境內(nèi)云貴準(zhǔn)靜止鋒影響下最高氣溫預(yù)報(bào)訂正方法研究[J].中低緯山地氣象, 2023,47(5):24-28.

[9] 邱貴強(qiáng),時(shí)少英,王洪霞,等.2m氣溫集成訂正方法及在冬奧延慶賽區(qū)的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2023,34(4):400-412.

[10] 吳姍,姚浪,劉健平,等.不同天氣現(xiàn)象下畢節(jié)高溫特征及預(yù)報(bào)回歸檢驗(yàn)分析[J].中低緯山地氣象,2022,46(1):68-72.

[11] 李國(guó)翠,連志鸞,趙彥廠,等.石家莊溫度預(yù)報(bào)檢驗(yàn)及影響因子分析[J].氣象與環(huán)境學(xué)報(bào),2009,25(1):15-18.

收稿日期:2024-07-10

作者簡(jiǎn)介:姜軍志(1996—),男,山東梁山人,助理工程師,研究方向?yàn)榫C合氣象業(yè)務(wù)。

主站蜘蛛池模板: 欧美成人精品欧美一级乱黄| 天天综合网色中文字幕| 国产日韩欧美一区二区三区在线| 免费看的一级毛片| 亚洲天堂免费| 国产小视频a在线观看| 99视频国产精品| 亚洲日韩精品伊甸| 夜夜爽免费视频| 亚洲日韩AV无码精品| 亚洲欧美天堂网| 麻豆国产精品一二三在线观看| 日韩视频精品在线| 久爱午夜精品免费视频| 青青草国产精品久久久久| 麻豆精选在线| 婷婷99视频精品全部在线观看| 久久国产精品国产自线拍| 伊人久久大线影院首页| 国产高清在线精品一区二区三区| 91精品国产无线乱码在线| 91精品国产91久无码网站| 亚洲天堂网在线播放| 72种姿势欧美久久久久大黄蕉| 国产精品爽爽va在线无码观看| 狠狠色狠狠综合久久| 色悠久久综合| 久久久精品国产SM调教网站| 色悠久久久久久久综合网伊人| 456亚洲人成高清在线| 亚洲日韩欧美在线观看| 国模私拍一区二区| JIZZ亚洲国产| 国产精品免费入口视频| 国产不卡在线看| 啊嗯不日本网站| 国内精品九九久久久精品| 精品亚洲麻豆1区2区3区| 少妇精品网站| 天天躁日日躁狠狠躁中文字幕| 久久婷婷综合色一区二区| 国产一在线观看| 亚洲一区二区三区国产精华液| 久久综合五月婷婷| 国产成人免费观看在线视频| 精品国产一二三区| 欧美在线视频不卡| 伊人成人在线| 国产成人禁片在线观看| 免费a级毛片18以上观看精品| 久久天天躁狠狠躁夜夜躁| 99热这里只有精品免费| 久久精品嫩草研究院| 国产成人亚洲毛片| 国产乱子伦视频在线播放| 久久人搡人人玩人妻精品一| 日韩欧美综合在线制服| 老司机久久99久久精品播放| 黄色网页在线播放| 东京热高清无码精品| 亚洲人成在线精品| 亚洲中文无码av永久伊人| 在线视频一区二区三区不卡| 久久精品亚洲中文字幕乱码| 久久夜色撩人精品国产| 97av视频在线观看| 色哟哟色院91精品网站| 精品亚洲欧美中文字幕在线看| 午夜日韩久久影院| 伊人福利视频| 国产精品极品美女自在线| 亚洲全网成人资源在线观看| 日本午夜影院| 亚洲精品久综合蜜| 亚洲最大情网站在线观看| 久久国产亚洲偷自| 四虎国产在线观看| 九九热精品视频在线| 国产在线观看91精品亚瑟| 在线色国产| 丰满人妻中出白浆| 亚洲AV色香蕉一区二区|