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

    雄安新區(qū)旱澇趨勢(shì)分析

    2019-06-24 02:36周浩偉劉蘇峽李峰平
    南水北調(diào)與水利科技 2019年2期
    關(guān)鍵詞:趨勢(shì)分析雄安新區(qū)

    周浩偉 劉蘇峽 李峰平

    摘要:旱澇趨勢(shì)對(duì)雄安新區(qū)規(guī)劃建設(shè)至關(guān)重要。增補(bǔ)構(gòu)建了雄安新區(qū)1469年到2018年的連續(xù)550 a歷史旱澇等級(jí)數(shù)據(jù)。以每10 a中的偏澇(旱)年份頻次組成的旱澇等級(jí)數(shù)據(jù)序列(后簡(jiǎn)稱偏澇、偏旱序列)作為研究對(duì)象,采用Lee-Heghinian法、有序聚類法、滑動(dòng)T檢驗(yàn)法識(shí)別序列的突變點(diǎn),采用線性趨勢(shì)回歸檢驗(yàn)、Spearman秩次法識(shí)別序列的變化趨勢(shì),采用R/S分析法對(duì)序列進(jìn)行預(yù)測(cè)分析。結(jié)果表明:1469-2018年,雄安新區(qū)偏澇序列在1579-1588年間存在突變點(diǎn),突變點(diǎn)前、后的子序列分別呈現(xiàn)出不顯著的上升趨勢(shì)和顯著的上升趨勢(shì);偏旱序列在1919-1928年間存在突變點(diǎn),突變點(diǎn)前、后的子序列分別呈現(xiàn)出不顯著的下降趨勢(shì)和不顯著的上升趨勢(shì)。Hurst指數(shù)綜合預(yù)測(cè)結(jié)果顯示,在0.05顯著性水平下,未來偏澇序列呈從顯著下降到不顯著上升趨勢(shì),偏旱序列呈從不顯著上升到顯著下降趨勢(shì)。

    關(guān)鍵詞:雄安新區(qū);旱澇等級(jí);趨勢(shì)分析;突變分析;R/S分析

    中圖分類號(hào):P426文獻(xiàn)標(biāo)志碼:A

    Abstract:Trends of drought and flood are very important for the planning and construction of Xiong′an New Area.In this paper,550 years of historical drought and flood grade data from 1469 to 2018 in Xiong′an New Area were constructed.The data series of drought and flood grade with frequency of drought-prone and flood-prone years in every 10 years were our research object,which were abbreviated as drought-prone and flood-prone sequences in our study.The turning points of the sequences were identified by Lee-Heghinian method,ordered clustering and moving t-test.The trends of the sequences were identified with linear trend regression test and Spearman rank method.The sequences were predicted with a R/S analysis method.The results showed that there was a turning point in 1579-1588 for the flood-prone sequence.The subsequence showed an insignificant upward trend before this turning point,and a significant upward trend after it.There was a turning point in 1919 to 1928 for the drought-prone sequence.The subsequence showed an insignificant downward trend before this turning point,and an insignificant upward trend after it.At the significant level of 0.05,Hurst index forecasting showed that the flood-prone sequence would show a trend varying from significantly decreasing to insignificantly increasing in the future,while the drought-prone sequence would show a trend varying from insignificantly increasing to significantly decreasing.

    Key words:Xiong′an New Area;drought and flood grades;trend analysis;turning point;R/S analysis

    2017年4月1日,作為北京市非首都功能疏解集中承載地的河北雄安新區(qū)正式成立。2018年4月21日,《河北雄安新區(qū)規(guī)劃綱要》(后簡(jiǎn)稱《綱要》)正式發(fā)布,《綱要》明確指出,到2035年要將雄安新區(qū)建設(shè)成為“綠色低碳、信息智能、宜居宜業(yè)、具有較強(qiáng)競(jìng)爭(zhēng)力和影響力、人與自然和諧共生的高水平社會(huì)主義現(xiàn)代化城市”。此次《綱要》中特地提到了 “人與自然和諧共生”的建設(shè)目標(biāo),而要實(shí)現(xiàn)這一目標(biāo),最重要的一點(diǎn)就是要降低自然災(zāi)害對(duì)人們正常生活生產(chǎn)的影響,旱澇災(zāi)害是對(duì)雄安新區(qū)影響重大的主要自然災(zāi)害[1-2]。因此,研究雄安新區(qū)旱澇趨勢(shì)對(duì)實(shí)現(xiàn)2035年建設(shè)目標(biāo)具有重要意義。

    建國(guó)以來,由于大范圍水文站網(wǎng)的建設(shè),有關(guān)干旱和洪澇的研究數(shù)據(jù)得到質(zhì)的提升,利用這些數(shù)據(jù)開展的研究也數(shù)不勝數(shù)。然而,要對(duì)未來進(jìn)行準(zhǔn)確預(yù)測(cè),還需要挖掘更長(zhǎng)年限的歷史數(shù)據(jù),探明干旱和洪澇在長(zhǎng)時(shí)間范圍尺度下的演變特征,對(duì)未來旱澇情勢(shì)做出資料長(zhǎng)度更具代表性的判斷。

    近年來,已有不少學(xué)者利用歷史旱澇資料在全國(guó)開展了有關(guān)旱澇特征的研究。例如:畢碩本等[3]通過分析西北地區(qū)東部1470-1912年連續(xù)443年旱澇等級(jí)序列,發(fā)現(xiàn)該地區(qū)旱澇災(zāi)害具有明顯的階段性,且通過對(duì)比發(fā)現(xiàn)采用歷史文獻(xiàn)資料重建的旱澇序列與其他長(zhǎng)降水量序列的一致性較好。袁媛等[4]通過分析陜西省幾個(gè)具有代表性的氣象站1470-2008年旱澇等級(jí)序列,發(fā)現(xiàn)榆林和延安呈現(xiàn)從偏旱轉(zhuǎn)向正常的趨勢(shì),寶雞和西安呈現(xiàn)出從偏旱轉(zhuǎn)向偏澇的趨勢(shì),而漢中和安康主要以偏澇為主。劉東升等[5]通過重建并分析昆明市1322-2013年旱澇等級(jí)序列,發(fā)現(xiàn)歷史上昆明市旱澇災(zāi)害總體變化趨勢(shì)為先趨于澇,后趨于旱,且旱澇災(zāi)害的第一、二、三主周期分別為39 a,11 a和4 a。常奐宇等[6]通過分析北京市1470-2015年旱澇等級(jí)資料,發(fā)現(xiàn)北京市旱澇災(zāi)害整體上呈現(xiàn)出“澇-旱-澇-旱”的波動(dòng),且局部存在旱澇急轉(zhuǎn)的現(xiàn)象。YU等[7]通過分析內(nèi)蒙古地區(qū)近500年旱澇資料,發(fā)現(xiàn)近500年來內(nèi)蒙古地區(qū)旱澇災(zāi)害發(fā)生的頻率和強(qiáng)度均有所增加,且存在澇災(zāi)總是在旱災(zāi)發(fā)生后的4~5 a發(fā)生的一個(gè)滯后效應(yīng)。Zheng等[8]通過分析1736-2000年中國(guó)北方歷史旱澇資料,發(fā)現(xiàn)在這一時(shí)段內(nèi)中國(guó)北方共發(fā)生了29次特大旱災(zāi)和28次特大洪災(zāi),且在這些極端干旱(洪澇)事件中,大部分地區(qū)的降水量明顯減少(增加),尤其是夏秋兩季。Jiang等[9]通過分析長(zhǎng)江三角洲地區(qū)近1000年來的旱澇資料,發(fā)現(xiàn)1000-1600 年洪災(zāi)發(fā)生頻次呈負(fù)趨勢(shì),其后呈正趨勢(shì),且1300 年以后旱災(zāi)發(fā)生頻次增加。

    本文利用歷史旱澇等級(jí)資料研究雄安新區(qū)旱澇趨勢(shì),旨在為雄安新區(qū)防洪抗旱規(guī)劃建設(shè)提供參考依據(jù)。

    1 研究區(qū)概況

    雄安新區(qū)地處河北省保定市,同時(shí)也處于北京市、天津市和保定市腹地,其規(guī)劃范圍涵蓋了保定市雄縣、容城、安新三縣的行政轄區(qū),任丘市鄚州鎮(zhèn)、茍各莊鎮(zhèn)、七間房鄉(xiāng)以及高陽縣龍化鄉(xiāng),總規(guī)劃面積約1 770 km2。

    雄安新區(qū)位處白洋淀流域出口。白洋淀流域主要包含了白洋淀及其周圍的大清河水系。大清河水系在白洋淀流域內(nèi)呈扇形分布,而水系中匯入白洋淀的河流按來水方向分為南北中三支。其中,大清河南支和中支的河流直接注入白洋淀內(nèi),而北支的河流則是通過白溝引河最終引入白洋淀[10]。

    白洋淀流域?qū)儆谂瘻貛О敫珊禋夂騾^(qū),多年平均降水量約為570 mm。整個(gè)流域呈現(xiàn)西北高東南低的地勢(shì)特征,降水量在空間分布上極不均勻[11],西北部山區(qū)的降雨量多于東南部平原區(qū)。流域多年平均徑流量約為22.3億m,山區(qū)河流徑流量占總徑流量的80%以上。流域降水量時(shí)間分配不均勻,大部分降水集中在6月-9月,約占流域總降水量的80%。暴雨多集中于每年7月、8月內(nèi),歷史上發(fā)生的幾次對(duì)人類影響較大的暴雨均發(fā)生在這段時(shí)間內(nèi)[12]。

    由于白洋淀流域降雨時(shí)空分布不均勻的氣候特征和雄安新區(qū)的特殊位置,雄安新區(qū)位于易發(fā)生洪澇和干旱災(zāi)害的區(qū)域。

    2 研究方法

    2.1 雄安新區(qū)1469-2018年連續(xù)550 a的歷史偏澇(旱)序列生成方法

    《中國(guó)近五百年旱澇分布圖集》[13](后簡(jiǎn)稱《圖集》)及其續(xù)補(bǔ)[14]和再續(xù)補(bǔ)資料[15]刊印了1470-2000年全國(guó)120個(gè)站點(diǎn)旱澇等級(jí)分布圖?!秷D集》采用 5個(gè)等級(jí)表示各地降水情況,其中1級(jí)代表澇、2級(jí)代表偏澇、3級(jí)代表正常、4級(jí)代表偏旱、5級(jí)代表旱,各等級(jí)描述見表1?!秷D集》中每一個(gè)站點(diǎn)代表了現(xiàn)行政區(qū)劃的1~2個(gè)地區(qū)范圍,本文選取了保定站數(shù)據(jù)代表雄安新區(qū)歷史旱澇等級(jí)情況?;谏鲜觥秷D集》及其續(xù)補(bǔ)和再續(xù)補(bǔ)資料,首先讀出保定站1470-2000年旱澇等級(jí)。

    考慮到后面計(jì)算每10 a的偏澇(旱)年份頻次需要整10 a的旱澇等級(jí)數(shù)據(jù),還需把歷史旱澇等級(jí)數(shù)據(jù)往前推1 a。利用《中國(guó)氣象災(zāi)害大典》(河北卷)[16]中關(guān)于雄安新區(qū)所在區(qū)域的記載,對(duì)照表1中有關(guān)各旱澇等級(jí)的描述,延展出1469 年的旱澇等級(jí)數(shù)據(jù)。

    根據(jù)中國(guó)氣象數(shù)據(jù)網(wǎng)中的保定市2001-2018年月降水量資料(6月-9月),按照《圖集》采用的評(píng)定標(biāo)準(zhǔn),延展出2001-2018年的旱澇等級(jí)數(shù)據(jù)。旱澇等級(jí)計(jì)算標(biāo)準(zhǔn)[14]如下。

    3 結(jié)果分析

    3.1 突變點(diǎn)分析

    不同突變點(diǎn)檢驗(yàn)方法各有優(yōu)劣,為了使檢驗(yàn)結(jié)果更加準(zhǔn)確,本文分別采用有序聚類法、Lee-Heghinian法和滑動(dòng)T檢驗(yàn)法三種方法相互對(duì)比,對(duì)雄安新區(qū)偏澇(旱)總序列在55個(gè)時(shí)段中的突變點(diǎn)進(jìn)行檢驗(yàn)。檢驗(yàn)結(jié)果見圖2。

    根據(jù)突變點(diǎn)檢驗(yàn)結(jié)果,三種檢驗(yàn)方法對(duì)偏澇(旱)總序列突變點(diǎn)的檢驗(yàn)結(jié)果一致。其中偏澇總序列突變點(diǎn)均位于第11個(gè)時(shí)段(1569-1578年),偏旱總序列突變點(diǎn)均位于第45個(gè)時(shí)段(1909-1918年)。此外,三種檢驗(yàn)方法對(duì)偏澇(旱)總序列的檢驗(yàn)結(jié)果均通過了0.05顯著性水平檢驗(yàn),這表明突變點(diǎn)檢驗(yàn)結(jié)果可靠性較高。

    3.2 突變點(diǎn)前、后的偏澇(旱)子序列趨勢(shì)分析

    為了探究突變點(diǎn)前、后的偏澇(旱)子序列統(tǒng)計(jì)值特征,分別計(jì)算了各子序列的均值和方差。計(jì)算結(jié)果見表2。

    根據(jù)表2可知:突變點(diǎn)前、后的偏澇(旱)子序列平均值的差值為1.36次/(10 a)(1.96次/(10 a)),這表明突變點(diǎn)前、后偏旱子序列的反轉(zhuǎn)比偏澇子序列更加劇烈。此外,突變點(diǎn)前、后的偏澇(旱)子序列方差的差值為0.75(0.2),這表明突變點(diǎn)前、后的偏澇子序列穩(wěn)定性差異比突變點(diǎn)前、后的偏旱子序列穩(wěn)定性差異更大。

    采用線性趨勢(shì)回歸檢驗(yàn)法和Spearman秩次法分別對(duì)突變點(diǎn)前、后的偏澇(旱)子序列變化趨勢(shì)進(jìn)行識(shí)別。結(jié)果顯示,兩種趨勢(shì)檢驗(yàn)方法對(duì)4組序列的檢驗(yàn)結(jié)果一致:在0.05顯著性水平下,突變點(diǎn)前、后的偏澇(旱)子序列分別呈現(xiàn)出不顯著的上升趨勢(shì)(不顯著的下降趨勢(shì))和顯著的上升趨勢(shì)(不顯著的上升趨勢(shì))。

    3.3 突變點(diǎn)后的偏澇(旱)子序列未來趨勢(shì)判斷

    以突變點(diǎn)后的偏澇(旱)子序列作為研究對(duì)象,根據(jù)R/S分析方法,分別作出突變點(diǎn)后的偏澇(旱)子序列l(wèi)g(R/S)-lgτ雙對(duì)數(shù)關(guān)系圖,并采用最小二乘法進(jìn)行線性擬合。擬合結(jié)果見圖4。

    從圖4可看出,突變點(diǎn)后的偏澇(旱)子序列通過擬合得到的Hurst指數(shù)為0.3682(0.6081),[JP+1]小于(大于)0.5,表明突變點(diǎn)后偏澇(旱)子序列的未來變化趨勢(shì)與歷史變化趨勢(shì)情況呈現(xiàn)出負(fù)相關(guān)(正相關(guān))。根據(jù)3.2節(jié)中趨勢(shì)檢驗(yàn)結(jié)果,突變點(diǎn)后偏澇(旱)子序列呈現(xiàn)出顯著的上升趨勢(shì)(不顯著的上升趨勢(shì)),可預(yù)測(cè)未來偏澇(旱)子序列將呈現(xiàn)出顯著的下降趨勢(shì)(不顯著的上升趨勢(shì))。由于突變點(diǎn)后偏澇(旱)子序列擬合趨勢(shì)線的確定性系數(shù)分別達(dá)到了0.7922(0.9859),且擬合的趨勢(shì)線方程均通過了0.05顯著性水平檢驗(yàn),因此預(yù)測(cè)結(jié)果較為可靠。

    4 討論

    4.1 偏澇(旱)總序列與突變點(diǎn)前、后的偏澇(旱)子序列趨勢(shì)對(duì)比

    采用滑動(dòng)平均法作出偏澇(旱)總序列的滑動(dòng)平均曲線,見圖5。

    從圖5可看出。偏澇總序列整體波動(dòng)較大,且在前10個(gè)時(shí)段(1469-1568年)以及第40-45個(gè)時(shí)段(1859-1918年),偏澇年份的頻次處于整個(gè)歷史時(shí)期中的較高水平,而在第12-18(1579-1648年)、第22-30個(gè)時(shí)段(1679-1768年),偏澇年份的頻次則處于整個(gè)歷史時(shí)期中的較低水平;同樣,偏旱總序列整體波動(dòng)性也比較大,且在第12-18(1579-1648年)、第45-55個(gè)時(shí)段(1909-2018年),偏旱年份的頻次處于整個(gè)歷史時(shí)期中的較高水平,而在第22-30(1679-1768年)、第40-45個(gè)時(shí)段(1859-1918年),偏旱年份的頻次則處于整個(gè)歷史時(shí)期中的較低水平。

    此外,偏澇(旱)總序列的均值為2.91次/(10 a)(3.2次/(10 a)),方差為2.86(3.53),這表明在過去550年里,雄安新區(qū)偏旱年份出現(xiàn)的頻次要高于偏澇年份出現(xiàn)的頻次,且偏旱總序列的波動(dòng)性比偏澇總序列更大。

    采用線性趨勢(shì)回歸檢驗(yàn)法和Spearman秩次法分別對(duì)偏澇(旱)總序列的變化趨勢(shì)進(jìn)行識(shí)別。結(jié)果顯示,兩種趨勢(shì)檢驗(yàn)方法對(duì)兩組序列的檢驗(yàn)結(jié)果一致:在0.05顯著性水平下,偏澇(旱)總序列呈現(xiàn)出不顯著的下降趨勢(shì)(顯著的上升趨勢(shì))。

    由此可見,考慮突變點(diǎn)的偏澇(旱)子序列變化趨勢(shì)與不考慮突變點(diǎn)的偏澇(旱)總序列變化趨勢(shì)差異明顯。

    4.2 Hurst指數(shù)預(yù)測(cè)結(jié)果討論

    4.2.1 突變點(diǎn)前、后的偏澇(旱)子序列趨勢(shì)對(duì)Hurst指數(shù)預(yù)測(cè)效果的驗(yàn)證

    以突變點(diǎn)前的偏澇(旱)子序列作為研究對(duì)象,根據(jù)R/S分析方法,分別做出突變點(diǎn)前的偏澇(旱)子序列l(wèi)g(R/S)-lgτ雙對(duì)數(shù)關(guān)系圖,并采用最小二乘法進(jìn)行線性擬合。擬合結(jié)果見圖6。

    從圖6可看出,突變點(diǎn)前的偏澇(旱)子序列擬合得到的Hurst指數(shù)為0.2589(0.2927),小于0.5,表明未來變化趨勢(shì)與歷史變化趨勢(shì)情況呈負(fù)相關(guān)。由3.2節(jié)中趨勢(shì)檢驗(yàn)結(jié)果可知,突變點(diǎn)前的偏澇(旱)子序列呈現(xiàn)出不顯著的上升趨勢(shì)(不顯著的下降趨勢(shì)),可預(yù)測(cè)突變點(diǎn)后的偏澇(旱)子序列將呈現(xiàn)出不顯著的下降趨勢(shì)(不顯著的上升趨勢(shì))。這個(gè)結(jié)果與突變后的偏旱子序列趨勢(shì)一致,但與突變后的[CM(22]偏澇子序列實(shí)際趨勢(shì)不一致(呈現(xiàn)顯著上升趨勢(shì))。

    通過分析發(fā)現(xiàn),造成突變點(diǎn)前的偏澇子序列預(yù)測(cè)失敗的主要原因有兩個(gè):其一,突變點(diǎn)前的偏澇子序列樣本點(diǎn)分布散亂,使得最終擬合出的趨勢(shì)線方程的確定性系數(shù)僅為0.3454,且趨勢(shì)線方程未通過0.05顯著性水平檢驗(yàn),不能保證Hurst指數(shù)的準(zhǔn)確性;其二,突變點(diǎn)前的偏澇子序列用于預(yù)測(cè)的時(shí)滯長(zhǎng)度為9,其預(yù)測(cè)結(jié)果很難保證對(duì)序列長(zhǎng)度大于9的未來序列有效,而突變點(diǎn)后的偏澇子序列長(zhǎng)度為44,遠(yuǎn)大于9。

    由此可見,采用Hurst指數(shù)對(duì)未來趨勢(shì)進(jìn)行預(yù)測(cè)時(shí),首先應(yīng)保證Hurst指數(shù)擬合效果較好,因?yàn)镠urst指數(shù)的擬合效果將直接影響預(yù)測(cè)結(jié)果的準(zhǔn)確性。其次,還應(yīng)注意預(yù)測(cè)所用的時(shí)滯長(zhǎng)度,通常預(yù)測(cè)所用的時(shí)滯長(zhǎng)度不同,最終得到的預(yù)測(cè)結(jié)果也有所不同。

    4.2.2 偏澇(旱)總序列與突變點(diǎn)后的偏澇(旱)子序列未來趨勢(shì)預(yù)測(cè)對(duì)比

    以偏澇(旱)總序列為研究對(duì)象,根據(jù)R/S分析方法,分別作出偏澇(旱)總序列的lg(R/S)-lg τ雙對(duì)數(shù)關(guān)系圖,并采用最小二乘法進(jìn)行線性擬合。擬合結(jié)果見圖7。

    從圖7可看出,偏澇(旱)總序列通過擬合得到的Hurst指數(shù)為0.3362(0.2981),小于0.5,表明偏澇(旱)總序列未來變化趨勢(shì)與歷史變化趨勢(shì)呈負(fù)相關(guān)。據(jù)4.1節(jié)中趨勢(shì)檢驗(yàn)結(jié)果可知,偏澇(旱)總序列呈現(xiàn)出不顯著的下降趨勢(shì)(顯著的上升趨勢(shì)),可預(yù)測(cè)未來偏澇(偏旱)序列將呈現(xiàn)出不顯著的上升趨勢(shì)(顯著的下降趨勢(shì))。由于偏澇(旱)總序列擬合趨勢(shì)線的確定性系數(shù)達(dá)到了0.8122(0.8349),且擬合的回歸方程通過了0.05顯著性水平檢驗(yàn),因此預(yù)測(cè)結(jié)果也較為可靠。

    通過對(duì)比采用偏澇(旱)總序列與采用突變點(diǎn)后的偏澇(旱)子序列進(jìn)行預(yù)測(cè)的結(jié)果發(fā)現(xiàn),兩種預(yù)測(cè)方式的預(yù)測(cè)結(jié)果不同,但由于這兩種預(yù)測(cè)方式的擬合效果都較好,因此兩種預(yù)測(cè)結(jié)果也都比較可靠。通過分析推測(cè),造成兩種預(yù)測(cè)方式結(jié)果不同的主要原因是偏澇(旱)總序列用于預(yù)測(cè)的時(shí)滯長(zhǎng)度與突變點(diǎn)后的偏澇(旱)子序列用于預(yù)測(cè)的時(shí)滯長(zhǎng)度不同,導(dǎo)致兩種預(yù)測(cè)方式所能預(yù)測(cè)的序列長(zhǎng)度不同,所以最終產(chǎn)生了兩種不同的預(yù)測(cè)結(jié)果。由此可見,利用Hurst指數(shù)進(jìn)行預(yù)測(cè)時(shí),時(shí)滯長(zhǎng)度這一單獨(dú)因素對(duì)未來趨勢(shì)的預(yù)測(cè)結(jié)果也有著很大的影響。

    有關(guān)Hurst指數(shù)計(jì)算的匯總結(jié)果見表3。

    綜合前面的分析,最終可以預(yù)測(cè):在0.05顯著性水平下,未來偏澇序列呈從顯著下降到不顯著上升趨勢(shì),偏旱序列呈從不顯著上升到顯著下降趨勢(shì)。限于方法的限制,未來趨勢(shì)究竟在未來多少時(shí)限內(nèi)成立還需要做進(jìn)一步研究。

    5 結(jié)論

    [JP+1]本文通過對(duì)雄安新區(qū)1469-2018年偏澇(旱)序列進(jìn)行突變點(diǎn)、趨勢(shì)以及預(yù)測(cè)分析,得出以下結(jié)論。

    (1)雄安新區(qū)偏澇、偏旱序列的突變點(diǎn)分別出現(xiàn)在第11個(gè)時(shí)段(1569-1578年)和第45個(gè)時(shí)段(1909-1918年)。突變點(diǎn)前、后偏澇(旱)子序列分別呈現(xiàn)出不顯著的上升趨勢(shì)(不顯著的下降趨勢(shì))和顯著的上升趨勢(shì)(不顯著的上升趨勢(shì))。

    (2)從序列總體變化趨勢(shì)上看,雄安新區(qū)偏澇(旱)總序列呈現(xiàn)出不顯著的下降趨勢(shì)(顯著的上升趨勢(shì))。其中,偏澇年份頻次的較大值出現(xiàn)在第1-10(1469-1568年)、第40-45個(gè)時(shí)段(1859-1918年)[JP+1],較小值出現(xiàn)在第12-18(1579-1648年)、第22-30個(gè)時(shí)段(1679-1768年);偏旱年份頻次的較大值出現(xiàn)在第12-18(1579-1648年)、第45-55個(gè)時(shí)段(1909-2018年),較小值出現(xiàn)在第22-30(1679-1768年)、第40-45個(gè)時(shí)段(1859-1918年)。

    (3)Hurst指數(shù)預(yù)測(cè)結(jié)果顯示,在0.05顯著性水平下,未來雄安新區(qū)偏澇序列呈從顯著下降到不顯著上升趨勢(shì),而偏旱序列呈從不顯著上升到顯著下降趨勢(shì)。

    參考文獻(xiàn)(References):

    [1] 夏軍,張永勇.雄安新區(qū)建設(shè)水安全保障面臨的問題與挑戰(zhàn)[J].中國(guó)科學(xué)院院刊,2017,32(11):1199-1205.(XIA J,ZHANG Y Y.Problems and challenges in the construction of water security in Xiong′an New Area[J].Journal of the Chinese Academy of Sciences,2017,32(11):1199-1205.(in Chinese))

    [2] 王宏偉.雄安新區(qū)災(zāi)害風(fēng)險(xiǎn)與彈性城市建設(shè)[J].中國(guó)減災(zāi),2017,37(23):30-33.(WANG H W.Disaster risk and flexible urban construction in Xiong′an New Area[J].Disaster Reduction in China,2017,37(23):30-33.(in Chinese))

    [3] 畢碩本,錢育君,陳昌春,等.1470-1912年西北東部地區(qū)旱澇等級(jí)重建序列的特征及對(duì)比分析[J].干旱區(qū)地理,2016,39(1):12-21.(BI S B,QIAN Y J,CHEN C C,et al.Characteristics and comparative analysis of the reconstruction sequence of drought and flood grades in the eastern part of Northwest China from 1470 to 1912[J].Arid Region Geography,2016,39(1):12-21.(in Chinese))

    [4] 袁媛,黃少妮,郭春燕,等.近500年陜西旱澇序列重建與趨勢(shì)分析[J].干旱地區(qū)農(nóng)業(yè)研究,2015,33(5):232-236,244.(YUAN Y,HUANG S N,GUO C Y,et al.Reconstruction and analysis the trend of drought and waterlogging sequence in shaanxi province in recent 500 years[J].Agricultural Research in Arid Areas,2015,33(5):232-236,244.(in Chinese))[HJ]

    [5] 劉東升,張虎才.昆明1322年至2013年旱澇歷史的重建[J].云南地理環(huán)境研究,2015,27(6):8-14,32.(LIU D S,ZHANG H C.Reconstruction of drought and flood history in Kunming from 1322 to 2013[J].Yunnan Geographic Environment Research,2015,27(6):8-14,32.(in Chinese))

    [6] 常奐宇,翟家齊,趙勇,等.基于馬爾可夫鏈的北京市546年來的旱澇演變特征[J].南水北調(diào)與水利科技,2018,16(5):27-34.(CHANG H Y,ZHAI J Q,ZHAO Y,et al.Evolution characteristics of drought and flood in Beijing in 546 years based on Markov chain[J].South-to-North Water Transfer and Water Conservancy Science and Technology,2018,16(5):27-34.(in Chinese))

    [7] YU P,LONG S F,MA J W,et al.Temporal-spatial variability in correlations of drought and flood during recent 500years in Inner Mongolia,China[J].Science of The Total Environment,2018,633:484-491.DOI:10.1016/j.scitotenv.2018.03.200.

    [8] ZHENG J Y,YU Y Z,ZHANG X Z,et al.Variation of extreme drought and flood in North China revealed by document-based seasonal precipitation reconstruction for the past 300 years[J].Climate Change of Past,2018,14(8):1135-1145.DOI:10.5194/cp-14-1135-2018.

    [9] JIANG T,ZHANG Q,BLENDER R,et al.Yangtze Delta floods and droughts of the last millennium:Abrupt changes and long term memory[J].Theoretical and Applied Climatology,2005,82(3-4):131-141.DOI:10,1007/s00704-005-0125-4.

    [10] [ZK(#]劉丹丹.白洋淀水資源量變化及其原因分析[D].保定:河北農(nóng)業(yè)大學(xué),2014.(LIU D D.Variation of water resources in Baiyangdian and analysis of its causes[D].Baoding:Agricultural University of Hebei,2014.(in Chinese))

    [11] 高彥春,王金鳳,封志明.白洋淀流域氣溫、降水和徑流變化特征及其相互響應(yīng)關(guān)系[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2017,25(4):467-477.(GAO Y C,WANG J F,F(xiàn)ENG Z M.Characteristics of temperature,precipitation and runoff variations and their relationship in the Baiyangdian basin[J].Chinese Journal of Eco-Agriculture,2017,25(4):467-477.(in Chinese))

    [12] 王蕊,王中根,夏軍,等.氣候變化及南水北調(diào)中線工程對(duì)白洋淀流域水文的影響[J].水力發(fā)電學(xué)報(bào),2012,31(2):14-19.(WANG R,WANG Z G,XIA J,et al.Effects of climate change and middle route project of South-to-North Water Transfer on hydrology in Baiyangdian basin[J].Journal of Hydropower,2012,31(2):14-19.(in Chinese))

    [13] 中央氣象局氣象科學(xué)研究院.中國(guó)近五百年旱澇分布圖集[M].北京:地圖出版社,1981.(Chinese academy of meteorological science,yearly charts of dryness/wetness in China for last 500-year period[M].Beijing:Sino Maps Press,1981 (in Chinese))[HJ2.05mm]

    [14] 張德二,劉傳志.《中國(guó)近五百年旱澇分布圖集》續(xù)補(bǔ)(1980-1992年)[J].氣象,1993,19(11):41-45.(ZHANG D E,LIU C Z.Continuation of atlas of drought and flood distribution in China in the past 500 years (1980-1992) [J].Meteorology Monthly,1993,19(11):41-45.(in Chinese))

    [15] 張德二,李小泉,梁有葉.《中國(guó)近五百年旱澇分布圖集》的再續(xù)補(bǔ)(1993-2000年)[J].應(yīng)用氣象學(xué)報(bào),2003,14(3):379-388.

    (ZHANG D E,LI X Q,LIANG Y Y.A further supplement to the atlas of drought and flood distribution in China in the past 500 years (1993-2000) [J].Journal of Applied Meteorology,2003,14(3):379-388.(in Chinese))

    [16] 臧建升.中國(guó)氣象災(zāi)害大典 河北卷[M].北京:氣象出版社,2008.(ZANG J S.Hebei volume of Chinese meteorological disaster code [M].Beijing:Meteorological Press,2008.(in Chinese))

    [17] 馬德栗,劉敏,鞠英芹.長(zhǎng)江流域及三峽庫(kù)區(qū)近542年旱澇演變特征[J].氣象科技,2016,44(4):622-630.(MA D L,LIU M,JU Y Q.Evolution characteristics of drought and flood in the Yangtze River Basin and the Three Gorges Reservoir in recent 542 years [J].Meteorological Science and Technology,2016,44(4):622-630.(in Chinese))

    [18] 黃振平,陳元芳.水文統(tǒng)計(jì)學(xué)[M].北京:中國(guó)水利水電出版社,2011.(HUANG Z P,CHEN Y F.Statistical hydrology[M].Beijing:China Water Resources and Hydropower Press,2011.(in Chinese))

    [19] 呂琳莉,劉湘?zhèn)ィ芗t梅,等.雅魯藏布江中下游年徑流[JP+2]變化趨勢(shì)分析[J].人民黃河,2013,35(5):27-29.(LYU L L,LIU X W,ZHOU H M,et al.Trend analysis of annual runoff change in the middle and lower reaches of the Yarlung Zangbo River[J].People′s Yellow River,2013,35(5):27-29.(in Chinese))

    [20] 杜懿.水文時(shí)間序列的特性分析及預(yù)測(cè)研究[D].南寧:廣西大學(xué),2018.(DU Y.Characteristic analysis and prediction of hydrological time series[D].Nanning:Guangxi University,2018.(in Chinese))

    [21] 陳宏,尉英華,王穎,等.基于VIC水文模型的灤河流域徑流變化特征及其影響因素[J].干旱氣象,2017,35(5):776-783.(CHEN H,WEI Y H,WANG Y,et al.Runoff variation characteristics and influencing factors of Luanhe River basin based on VIC hydrological model[J].Drought Meteorology,2017,35(5):776-783.(in Chinese))

    [22] [ZK(#]周園園,師長(zhǎng)興,范小黎,等.國(guó)內(nèi)水文序列變異點(diǎn)分析方法及在各流域應(yīng)用研究進(jìn)展[J].地理科學(xué)進(jìn)展,2011,30(11):1361-1369.(ZHOU Y Y,SHI C X,F(xiàn)AN X L,et al.Advances in the analysis methods of hydrological sequence variation points and their application in various watersheds in China [J].Progress in Geographic Science,2011,30(11):1361-1369.(in Chinese))

    [23] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京:氣象出版社,2007.(WEI F Y.Modern climate statistical diagnosis and forecasting technology[M].Beijing:Meteorological Press,2007.(in Chinese))

    [24] 魯帆,嚴(yán)登華,王勇,等.中長(zhǎng)期徑流預(yù)報(bào)技術(shù)與方法[M].北京:中國(guó)水利水電出版社,2012.(LU F,YAN D H,WANG Y,et al.Medium and long term runoff forecasting techniques and methods [M].Beijing:China Water Resources and Hydropower Press,2012.(in Chinese))

    [25] 門寶輝,劉昌明,夏軍,等.R/S分析法在南水北調(diào)西線一期工程調(diào)水河流徑流趨勢(shì)預(yù)測(cè)中的應(yīng)用[J].冰川凍土,2005(4):568-573.(MEN B H,LIU C M,XIA J,et al.Application of R/S analysis method in prediction of runoff trend of water transfer river in the first stage of West Route of South-to-North Water Transfer Project [J].Glacier Frozen Soil,2005(4):568-573.(in Chinese))

    [26] 丁之勇,魯瑞潔,劉暢,等.環(huán)青海湖地區(qū)氣候變化特征及其季風(fēng)環(huán)流因素[J].地球科學(xué)進(jìn)展,2018,33(3):281-292.(DING Z Y,LU R J,LIU C,et al.Climate change characteristics and monsoon circulation factors around Qinghai Lake [J].Advances in Earth Science,2018,33(3):281-292.(in Chinese))

    [27] 陳建龍,劉永峰,錢鞠,等.R/S分析法與GM(1,1)灰色模型相結(jié)合的鴛鴦池水庫(kù)入庫(kù)徑流量預(yù)測(cè)[J].水資源與水工程學(xué)報(bào),2018,29(5):148-153+158.(CHEN J L,LIU Y F,QIAN J,et al.R/S analysis combined with GM(1,1) grey model [J].Journal of Water Resources and Water Engineering,2018,29(5):148-153+158.(in Chinese))

    [28] 陳彥光.基于Matlab的地理數(shù)據(jù)分析[M].北京:高等教育出版社,2012.(CHEN Y G.Analysis of geographic data based on MATLAB[M].Beijing:Higher Education Press,2012.(in Chinese))

    猜你喜歡
    趨勢(shì)分析雄安新區(qū)
    河北雄安新區(qū)發(fā)展SWOT分析及對(duì)策
    認(rèn)清建設(shè)雄安新區(qū)的根本戰(zhàn)略意圖
    唱衰雄安新區(qū)論調(diào)為什么站不住腳
    楊誼興:積極參與“雄安新區(qū)”機(jī)關(guān)企事業(yè)單位食堂社會(huì)化改造項(xiàng)目
    楊誼興:積極參與“雄安新區(qū)”機(jī)關(guān)企事業(yè)單位食堂社會(huì)化改造項(xiàng)目
    化工材料科學(xué)與工程的發(fā)展現(xiàn)狀及趨勢(shì)研究
    淺談核電項(xiàng)目設(shè)計(jì)質(zhì)量信息收集與分析
    高校微課教學(xué)現(xiàn)狀及未來發(fā)展趨勢(shì)研究
    從《遠(yuǎn)程教育》35年載文看遠(yuǎn)程教育研究趨勢(shì)
    基于知識(shí)圖譜的智慧教育研究熱點(diǎn)與趨勢(shì)分析
    久久人人精品亚洲av| 无人区码免费观看不卡| 亚洲第一电影网av| 国产精品1区2区在线观看.| 欧美日本亚洲视频在线播放| 久久中文字幕一级| 国产视频内射| 国产午夜精品论理片| 日本黄大片高清| 色噜噜av男人的天堂激情| 舔av片在线| 操出白浆在线播放| 国产精品野战在线观看| 久久久久亚洲av毛片大全| 亚洲精品美女久久久久99蜜臀| 99久久无色码亚洲精品果冻| 午夜福利成人在线免费观看| 欧美日韩瑟瑟在线播放| 国产一区在线观看成人免费| 国产精品精品国产色婷婷| 给我免费播放毛片高清在线观看| 在线播放国产精品三级| 久热爱精品视频在线9| 中文字幕人成人乱码亚洲影| 搡老妇女老女人老熟妇| 欧美日韩亚洲国产一区二区在线观看| 亚洲av五月六月丁香网| 色综合亚洲欧美另类图片| 两个人视频免费观看高清| 夜夜夜夜夜久久久久| 淫秽高清视频在线观看| 99久久国产精品久久久| 熟女少妇亚洲综合色aaa.| 又爽又黄无遮挡网站| 可以在线观看的亚洲视频| 黄片大片在线免费观看| 亚洲无线在线观看| АⅤ资源中文在线天堂| 国产高清视频在线观看网站| 视频区欧美日本亚洲| videosex国产| 亚洲性夜色夜夜综合| 男女之事视频高清在线观看| 夜夜爽天天搞| 亚洲精品久久国产高清桃花| 制服诱惑二区| 久久久精品国产亚洲av高清涩受| 狂野欧美激情性xxxx| 欧美大码av| 亚洲午夜精品一区,二区,三区| 不卡一级毛片| 久久中文看片网| 亚洲一码二码三码区别大吗| 天堂影院成人在线观看| 麻豆久久精品国产亚洲av| 国产高清视频在线观看网站| 一本一本综合久久| 看片在线看免费视频| 麻豆久久精品国产亚洲av| 亚洲精品色激情综合| 高清毛片免费观看视频网站| 看免费av毛片| 午夜免费成人在线视频| 亚洲av日韩精品久久久久久密| 在线观看免费午夜福利视频| 国产熟女xx| 午夜福利在线在线| 狠狠狠狠99中文字幕| 成年女人毛片免费观看观看9| 国产伦一二天堂av在线观看| 国产精品香港三级国产av潘金莲| 久久精品成人免费网站| 亚洲欧美日韩高清在线视频| 又粗又爽又猛毛片免费看| 欧美日韩中文字幕国产精品一区二区三区| bbb黄色大片| 天堂av国产一区二区熟女人妻 | 叶爱在线成人免费视频播放| 人妻久久中文字幕网| 天天一区二区日本电影三级| 国产精品,欧美在线| 日韩大码丰满熟妇| 国产激情偷乱视频一区二区| 亚洲精品在线观看二区| 亚洲aⅴ乱码一区二区在线播放 | www.999成人在线观看| 亚洲国产精品成人综合色| 国产精品 欧美亚洲| 变态另类成人亚洲欧美熟女| 亚洲精品美女久久久久99蜜臀| 国产av又大| 国产成人av激情在线播放| 麻豆久久精品国产亚洲av| 国产午夜精品久久久久久| av福利片在线| 国产免费男女视频| 国产91精品成人一区二区三区| 久久久久九九精品影院| 国产av麻豆久久久久久久| 又粗又爽又猛毛片免费看| 欧美性猛交╳xxx乱大交人| 两人在一起打扑克的视频| 人人妻人人看人人澡| 嫩草影视91久久| 亚洲性夜色夜夜综合| 久久久久久亚洲精品国产蜜桃av| 黄色成人免费大全| 亚洲精品美女久久av网站| 在线观看免费午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看 | 淫秽高清视频在线观看| 亚洲avbb在线观看| 999久久久精品免费观看国产| 国产高清视频在线观看网站| 日本撒尿小便嘘嘘汇集6| 亚洲一区二区三区不卡视频| av国产免费在线观看| 免费搜索国产男女视频| 一级毛片女人18水好多| 无遮挡黄片免费观看| 狠狠狠狠99中文字幕| 亚洲第一欧美日韩一区二区三区| 在线十欧美十亚洲十日本专区| 国产不卡一卡二| 亚洲黑人精品在线| 精品国产乱子伦一区二区三区| 国产午夜福利久久久久久| 亚洲精品粉嫩美女一区| 欧美黄色淫秽网站| 777久久人妻少妇嫩草av网站| 国产欧美日韩一区二区三| 18禁裸乳无遮挡免费网站照片| 久久国产精品人妻蜜桃| 日本免费一区二区三区高清不卡| 男女之事视频高清在线观看| 欧美久久黑人一区二区| 女人高潮潮喷娇喘18禁视频| 在线观看午夜福利视频| 99热这里只有是精品50| 亚洲精品中文字幕在线视频| 亚洲精品国产精品久久久不卡| x7x7x7水蜜桃| 欧美不卡视频在线免费观看 | 一区二区三区国产精品乱码| 村上凉子中文字幕在线| 国产成人影院久久av| 国产99久久九九免费精品| 日韩欧美国产一区二区入口| 欧美在线一区亚洲| 精品久久久久久久久久久久久| 国产一级毛片七仙女欲春2| 国产一区二区激情短视频| 美女 人体艺术 gogo| 国产午夜精品论理片| 久久精品夜夜夜夜夜久久蜜豆 | 欧美 亚洲 国产 日韩一| 丁香欧美五月| 亚洲18禁久久av| 真人一进一出gif抽搐免费| www.精华液| 欧美日本视频| 黄频高清免费视频| 久久草成人影院| 久久久精品大字幕| 免费一级毛片在线播放高清视频| 欧美 亚洲 国产 日韩一| 一边摸一边做爽爽视频免费| 国产成人系列免费观看| 宅男免费午夜| 欧美+亚洲+日韩+国产| 999久久久国产精品视频| 中文字幕人成人乱码亚洲影| 麻豆成人午夜福利视频| www日本在线高清视频| 久久99热这里只有精品18| 在线观看免费午夜福利视频| 在线a可以看的网站| 十八禁网站免费在线| 欧美乱妇无乱码| 淫妇啪啪啪对白视频| 亚洲av成人不卡在线观看播放网| 日本a在线网址| 高清在线国产一区| 999久久久国产精品视频| 香蕉丝袜av| 最新美女视频免费是黄的| 999久久久国产精品视频| 制服丝袜大香蕉在线| 中文字幕最新亚洲高清| 亚洲欧美日韩无卡精品| 欧美日韩福利视频一区二区| 在线a可以看的网站| 色综合站精品国产| 91大片在线观看| av中文乱码字幕在线| 18美女黄网站色大片免费观看| 99精品在免费线老司机午夜| 国产高清激情床上av| 精品欧美国产一区二区三| 久久久久免费精品人妻一区二区| 十八禁网站免费在线| 可以免费在线观看a视频的电影网站| 日本黄色视频三级网站网址| 亚洲国产精品成人综合色| 国产99久久九九免费精品| 巨乳人妻的诱惑在线观看| 日韩大码丰满熟妇| 国产在线精品亚洲第一网站| 亚洲真实伦在线观看| 亚洲黑人精品在线| 亚洲国产欧美网| 婷婷六月久久综合丁香| 女人被狂操c到高潮| 久热爱精品视频在线9| 亚洲国产欧洲综合997久久,| 国内久久婷婷六月综合欲色啪| 97碰自拍视频| 后天国语完整版免费观看| 熟妇人妻久久中文字幕3abv| 女生性感内裤真人,穿戴方法视频| 国产v大片淫在线免费观看| 熟女电影av网| 最近最新中文字幕大全免费视频| 日日爽夜夜爽网站| 久久精品aⅴ一区二区三区四区| 在线观看免费午夜福利视频| 亚洲 欧美 日韩 在线 免费| 欧美zozozo另类| 真人做人爱边吃奶动态| 欧美黑人欧美精品刺激| 色尼玛亚洲综合影院| 欧美又色又爽又黄视频| 最近在线观看免费完整版| www日本在线高清视频| 床上黄色一级片| 色av中文字幕| 亚洲乱码一区二区免费版| 一本久久中文字幕| av免费在线观看网站| 天天添夜夜摸| 99精品在免费线老司机午夜| 1024香蕉在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲美女黄片视频| 久久久精品国产亚洲av高清涩受| 12—13女人毛片做爰片一| 久久久久久久精品吃奶| 久久精品亚洲精品国产色婷小说| 九九热线精品视视频播放| 国产精品永久免费网站| 亚洲国产精品999在线| 成人av在线播放网站| 婷婷亚洲欧美| 日本免费a在线| 精品国内亚洲2022精品成人| 国产精品av视频在线免费观看| 99久久精品国产亚洲精品| 国产一级毛片七仙女欲春2| 久久久久久久精品吃奶| 免费在线观看视频国产中文字幕亚洲| 最近视频中文字幕2019在线8| 精品久久久久久久久久免费视频| 国内揄拍国产精品人妻在线| 国产v大片淫在线免费观看| 日本免费一区二区三区高清不卡| tocl精华| 成人高潮视频无遮挡免费网站| 久久久久久久午夜电影| 亚洲欧美日韩高清在线视频| 亚洲,欧美精品.| 精品久久蜜臀av无| 日韩欧美精品v在线| 国产成人aa在线观看| 小说图片视频综合网站| 欧美在线一区亚洲| 国产伦在线观看视频一区| АⅤ资源中文在线天堂| 夜夜躁狠狠躁天天躁| 久久精品国产亚洲av香蕉五月| 在线看三级毛片| 精品一区二区三区视频在线观看免费| a在线观看视频网站| 欧美激情久久久久久爽电影| 性色av乱码一区二区三区2| 可以在线观看毛片的网站| 1024香蕉在线观看| 色尼玛亚洲综合影院| 搡老岳熟女国产| 久久这里只有精品19| 日本黄大片高清| 亚洲精华国产精华精| 国产伦一二天堂av在线观看| 丝袜美腿诱惑在线| 午夜日韩欧美国产| 天堂√8在线中文| 日韩欧美国产在线观看| 国产免费av片在线观看野外av| 免费看日本二区| 亚洲av成人不卡在线观看播放网| 狂野欧美白嫩少妇大欣赏| 精品国产亚洲在线| 亚洲五月婷婷丁香| 全区人妻精品视频| 这个男人来自地球电影免费观看| 国产精品98久久久久久宅男小说| 欧美黑人巨大hd| 一本综合久久免费| 精品国产乱码久久久久久男人| 可以在线观看毛片的网站| 免费在线观看完整版高清| svipshipincom国产片| 亚洲成人中文字幕在线播放| 久久久国产成人精品二区| 国产一级毛片七仙女欲春2| 一个人免费在线观看电影 | 91字幕亚洲| 波多野结衣高清作品| 欧美极品一区二区三区四区| 国产精品一区二区三区四区免费观看 | 国产精品影院久久| 91av网站免费观看| 国产精品亚洲av一区麻豆| 一区二区三区激情视频| 国产伦人伦偷精品视频| 观看免费一级毛片| 人人妻人人澡欧美一区二区| 久久性视频一级片| 婷婷亚洲欧美| 亚洲无线在线观看| 国产成人av激情在线播放| 一区二区三区高清视频在线| 亚洲欧美一区二区三区黑人| 亚洲avbb在线观看| 悠悠久久av| 精品久久久久久久末码| 脱女人内裤的视频| 国产精品久久久av美女十八| 男女视频在线观看网站免费 | 免费在线观看完整版高清| 在线看三级毛片| 久久香蕉精品热| 国产精品自产拍在线观看55亚洲| 精品一区二区三区av网在线观看| 在线视频色国产色| 成年人黄色毛片网站| 搞女人的毛片| 亚洲av成人av| 人妻久久中文字幕网| 久久久久久人人人人人| 国产精品精品国产色婷婷| 观看免费一级毛片| 久久国产乱子伦精品免费另类| 男插女下体视频免费在线播放| 国产成+人综合+亚洲专区| 欧美黄色片欧美黄色片| 成人国语在线视频| 一区二区三区激情视频| 制服丝袜大香蕉在线| 午夜精品一区二区三区免费看| 少妇人妻一区二区三区视频| 老鸭窝网址在线观看| 国产精品一区二区免费欧美| 老司机靠b影院| 丝袜美腿诱惑在线| 欧美大码av| 国产男靠女视频免费网站| 久久精品成人免费网站| 欧美高清成人免费视频www| 国产精品久久久久久精品电影| 97超级碰碰碰精品色视频在线观看| 在线观看免费视频日本深夜| 国产精品自产拍在线观看55亚洲| 制服人妻中文乱码| 亚洲成人精品中文字幕电影| 成人手机av| www日本在线高清视频| 天堂√8在线中文| 亚洲男人的天堂狠狠| 国内少妇人妻偷人精品xxx网站 | 亚洲五月婷婷丁香| 欧美日本视频| 午夜a级毛片| 高清毛片免费观看视频网站| 人人妻人人看人人澡| 美女免费视频网站| 欧美成人免费av一区二区三区| 精华霜和精华液先用哪个| 亚洲第一欧美日韩一区二区三区| 一区福利在线观看| 午夜成年电影在线免费观看| 久久中文看片网| 淫妇啪啪啪对白视频| 一二三四社区在线视频社区8| 久久人妻av系列| 日本精品一区二区三区蜜桃| 天堂av国产一区二区熟女人妻 | 午夜精品久久久久久毛片777| 成人欧美大片| 欧美在线黄色| 麻豆av在线久日| 可以在线观看毛片的网站| 18禁黄网站禁片免费观看直播| 好男人在线观看高清免费视频| 白带黄色成豆腐渣| 午夜福利免费观看在线| 天堂影院成人在线观看| 成人三级黄色视频| 亚洲人成电影免费在线| 变态另类成人亚洲欧美熟女| 国产精品自产拍在线观看55亚洲| 怎么达到女性高潮| 国产主播在线观看一区二区| 成年免费大片在线观看| 12—13女人毛片做爰片一| 色精品久久人妻99蜜桃| 国产1区2区3区精品| 波多野结衣高清无吗| 精品免费久久久久久久清纯| xxx96com| av中文乱码字幕在线| 麻豆国产av国片精品| 国产亚洲av高清不卡| 亚洲aⅴ乱码一区二区在线播放 | 天堂√8在线中文| 最新美女视频免费是黄的| 精华霜和精华液先用哪个| 精品无人区乱码1区二区| 一区二区三区高清视频在线| 久久久久国内视频| 久久久国产精品麻豆| 欧美人与性动交α欧美精品济南到| 熟妇人妻久久中文字幕3abv| 变态另类丝袜制服| 日韩欧美国产在线观看| www日本黄色视频网| 精品一区二区三区四区五区乱码| 日本黄色视频三级网站网址| 一个人观看的视频www高清免费观看 | 我要搜黄色片| 亚洲av中文字字幕乱码综合| 国产精品,欧美在线| 精品久久久久久成人av| 在线视频色国产色| 久久久精品国产亚洲av高清涩受| 黄色 视频免费看| 好男人在线观看高清免费视频| 成人国产一区最新在线观看| 成人国语在线视频| 色综合亚洲欧美另类图片| 人妻丰满熟妇av一区二区三区| 视频区欧美日本亚洲| 好看av亚洲va欧美ⅴa在| 精品免费久久久久久久清纯| 国产亚洲精品久久久久久毛片| 麻豆成人午夜福利视频| 亚洲自拍偷在线| 午夜成年电影在线免费观看| 长腿黑丝高跟| 香蕉久久夜色| 一边摸一边抽搐一进一小说| 一进一出抽搐gif免费好疼| 国产亚洲欧美98| 一区福利在线观看| 精品国产乱子伦一区二区三区| or卡值多少钱| 午夜激情av网站| 欧美成狂野欧美在线观看| 亚洲av成人不卡在线观看播放网| 亚洲avbb在线观看| 亚洲av成人av| 特级一级黄色大片| 高清毛片免费观看视频网站| 99国产综合亚洲精品| 手机成人av网站| 两性夫妻黄色片| 亚洲专区中文字幕在线| 老司机福利观看| 亚洲av美国av| 91九色精品人成在线观看| 免费人成视频x8x8入口观看| 亚洲电影在线观看av| 久久久精品国产亚洲av高清涩受| av在线播放免费不卡| 国产成人av激情在线播放| 嫩草影院精品99| 国产精品一区二区精品视频观看| 在线观看日韩欧美| 国产一区在线观看成人免费| 高潮久久久久久久久久久不卡| av在线天堂中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 欧美人与性动交α欧美精品济南到| 久久人人精品亚洲av| www.熟女人妻精品国产| 久久香蕉国产精品| 久久久久国产精品人妻aⅴ院| 毛片女人毛片| 精品久久久久久,| 国产高清videossex| 琪琪午夜伦伦电影理论片6080| 757午夜福利合集在线观看| 久久性视频一级片| 久久人人精品亚洲av| 天堂影院成人在线观看| 中文字幕熟女人妻在线| 一级黄色大片毛片| 99热6这里只有精品| 中文字幕人成人乱码亚洲影| 国产探花在线观看一区二区| 国产精品自产拍在线观看55亚洲| 久热爱精品视频在线9| 18禁黄网站禁片午夜丰满| 国产不卡一卡二| 免费在线观看日本一区| 国产精品av久久久久免费| 亚洲精品一卡2卡三卡4卡5卡| 嫁个100分男人电影在线观看| 亚洲国产高清在线一区二区三| 欧美人与性动交α欧美精品济南到| 少妇人妻一区二区三区视频| 精品一区二区三区av网在线观看| 国产伦人伦偷精品视频| 亚洲午夜精品一区,二区,三区| 99在线视频只有这里精品首页| 国内少妇人妻偷人精品xxx网站 | 免费在线观看日本一区| 久久精品国产亚洲av高清一级| 不卡一级毛片| 欧美一级毛片孕妇| 国产av不卡久久| 757午夜福利合集在线观看| 真人一进一出gif抽搐免费| 天天躁夜夜躁狠狠躁躁| 亚洲精品国产一区二区精华液| 国产又色又爽无遮挡免费看| 国内精品久久久久精免费| 90打野战视频偷拍视频| 看黄色毛片网站| 亚洲一码二码三码区别大吗| 12—13女人毛片做爰片一| 老司机福利观看| 黄色视频不卡| 色老头精品视频在线观看| 亚洲人成电影免费在线| 欧美黑人精品巨大| 人人妻人人看人人澡| 18美女黄网站色大片免费观看| 伊人久久大香线蕉亚洲五| 亚洲精品国产一区二区精华液| 久久精品影院6| 欧美激情久久久久久爽电影| 亚洲成人久久爱视频| 久久国产乱子伦精品免费另类| 国产aⅴ精品一区二区三区波| 欧美日韩亚洲综合一区二区三区_| 日本黄色视频三级网站网址| 免费看日本二区| 免费高清视频大片| 亚洲国产日韩欧美精品在线观看 | 日韩三级视频一区二区三区| 岛国在线观看网站| 别揉我奶头~嗯~啊~动态视频| 国内揄拍国产精品人妻在线| 色综合婷婷激情| 十八禁人妻一区二区| 成人国产一区最新在线观看| 成人18禁在线播放| 久久久久久久午夜电影| 一级黄色大片毛片| 首页视频小说图片口味搜索| www.999成人在线观看| 又黄又爽又免费观看的视频| 国产麻豆成人av免费视频| 精华霜和精华液先用哪个| 制服人妻中文乱码| 亚洲av第一区精品v没综合| 精品国产超薄肉色丝袜足j| 91麻豆av在线| 色综合站精品国产| 国语自产精品视频在线第100页| 国内精品久久久久久久电影| 国产精品免费一区二区三区在线| netflix在线观看网站| 久久精品成人免费网站| 免费看a级黄色片| 50天的宝宝边吃奶边哭怎么回事| 91av网站免费观看| 日韩欧美 国产精品| 亚洲男人的天堂狠狠| 国产精品自产拍在线观看55亚洲| 亚洲人成网站在线播放欧美日韩| 狂野欧美激情性xxxx| 亚洲国产欧美人成| 亚洲专区国产一区二区| 欧美成狂野欧美在线观看| 午夜a级毛片| 韩国av一区二区三区四区| 香蕉国产在线看| 国产主播在线观看一区二区| 国产熟女xx| 别揉我奶头~嗯~啊~动态视频| 在线永久观看黄色视频| 欧美 亚洲 国产 日韩一| 国产真实乱freesex| 午夜免费观看网址| 特级一级黄色大片| 女警被强在线播放| 91麻豆精品激情在线观看国产| 亚洲最大成人中文| 午夜激情福利司机影院| 国产成年人精品一区二区| 精品福利观看| 国产精品亚洲av一区麻豆| 国产一区二区激情短视频| 成人特级黄色片久久久久久久|