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

        ?

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

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

        摘 要:選取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)確率較高,但在開展鄉(xiāng)鎮(zhèn)預(yù)報(bào)時(shí),通常以基本站為基準(zhǔn),采用經(jīng)驗(yàn)訂正的方式開展氣象服務(wù),此方法無(wú)法滿足當(dāng)前社會(huì)需求。近年來(lái),鄉(xiāng)鎮(zhèn)區(qū)域站網(wǎng)布局越發(fā)完善,氣象數(shù)據(jù)更加全面準(zhǔn)確,為更好地開展鄉(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)方程具體見表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í),易損害氣象部門公信力[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ì)黔西市開展精細(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 ℃,由此可見,各鄉(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ù)。

        猜你喜歡
        最小二乘法
        基于壓電陶瓷的納米定位與掃描平臺(tái)模型辨識(shí)算法研究
        基于慣導(dǎo)角度量測(cè)的軌道平面最佳線形參數(shù)估計(jì)算法
        馬爾科夫鏈在市場(chǎng)預(yù)測(cè)中的應(yīng)用
        一種改進(jìn)的基于RSSI最小二乘法和擬牛頓法的WSN節(jié)點(diǎn)定位算法
        某型飛機(jī)副翼操縱系統(tǒng)建模研究
        最小二乘法基本思想及其應(yīng)用
        科技視界(2016年22期)2016-10-18 17:12:13
        基于改進(jìn)的最小二乘法的車道線檢測(cè)方法
        全國(guó)各省份經(jīng)濟(jì)發(fā)展影響因素的實(shí)證分析
        全國(guó)主要市轄區(qū)的房?jī)r(jià)收入比影響因素研究
        商(2016年20期)2016-07-04 01:23:26
        手動(dòng)求解線性回歸方程的方法和技巧
        考試周刊(2016年27期)2016-05-26 23:45:13
        国产乱妇无码大片在线观看| 亚洲一区二区女优av| 日本国产一区在线观看| 亚洲最大成人网站| 无码人妻久久一区二区三区不卡| 国产黄色免费网站| 国产精品人人爱一区二区白浆| 手机在线播放av网址| 欧美激情综合色综合啪啪五月| 99久久精品费精品国产一区二区| 无码日韩人妻AV一区免费| av一区二区不卡久久| 漂亮丰满人妻被中出中文字幕| 亚洲国产av无码专区亚洲av| 在线播放人成午夜免费视频| 中文字幕在线人妻视频| 国产成人高清在线观看视频 | 亚洲精品美女自拍偷拍| 国产影院一区二区在线| 欧美国产激情二区三区| 亚洲国产人在线播放首页 | 乱子伦视频在线看| 成年视频网站在线观看777| 免费在线视频亚洲色图| 玩弄丰满奶水的女邻居| 三上悠亚精品一区二区久久| 丰满人妻一区二区三区免费| 熟妇高潮一区二区三区在线观看| 性生交大全免费看| 亚洲日本在线va中文字幕| 中文字幕色资源在线视频| 久热re这里精品视频在线6| 久久香蕉免费国产天天看| 中文字幕二区三区在线| 国产av在线观看久久| 国产欧美日韩一区二区三区在线| 日韩毛片久久91| 久久综合国产精品一区二区| 精品人妻伦九区久久aaa片| 久久久久国产一级毛片高清版A| 中文片内射在线视频播放|