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

    工業(yè)用水量測算方法及應(yīng)用

    2017-02-23 02:14:58周申蓓劉亞靈鄭士鵬孫淑云
    水利經(jīng)濟 2017年1期
    關(guān)鍵詞:方法模型

    周申蓓,劉亞靈,鄭士鵬,2,3,孫淑云

    (1.河海大學(xué)商學(xué)院,江蘇 南京 211100; 2.浙江省水利河口研究院,浙江 杭州 310020;3.臺州市水利水電勘測設(shè)計院,浙江 臺州 318000;4.水利部綜合事業(yè)局,北京 100053)

    工業(yè)用水量測算方法及應(yīng)用

    周申蓓1,劉亞靈1,鄭士鵬1,2,3,孫淑云4

    (1.河海大學(xué)商學(xué)院,江蘇 南京 211100; 2.浙江省水利河口研究院,浙江 杭州 310020;3.臺州市水利水電勘測設(shè)計院,浙江 臺州 318000;4.水利部綜合事業(yè)局,北京 100053)

    工業(yè)用水總量測算是我國實行最嚴(yán)格水資源管理制度的一項重要工作。為了利用現(xiàn)有水利統(tǒng)計成果對工業(yè)用水總量進(jìn)行測算,文章選取了基于歷史數(shù)據(jù)的預(yù)測方法和計量結(jié)合抽樣的統(tǒng)計方法測算工業(yè)用水量,提出了“以計量結(jié)合抽樣法為主,趨勢外推法進(jìn)行對比校核”的工業(yè)用水量測算方法。結(jié)果表明:對10萬m3以上規(guī)模的工業(yè)采用全樣本計量,對規(guī)模以下的工業(yè)按一定比例進(jìn)行抽樣計量,測算值的相對誤差可控制在10%左右。對于技術(shù)和統(tǒng)計投入受限的省份,可利用2011年水利普查數(shù)據(jù)和其他歷史數(shù)據(jù),適當(dāng)采用趨勢外推法對原有數(shù)據(jù)進(jìn)行修正。

    工業(yè)用水量;統(tǒng)計;測算方法;水資源管理

    長期以來,水資源利用方式粗放、利用效率不高和水資源供需矛盾突出是制約我國可持續(xù)發(fā)展的主要瓶頸之一[1]。為解決我國日趨復(fù)雜的水資源問題,使水資源得以高效利用和有效保護(hù),我國在2013年提出了最嚴(yán)格水資源管理制度,其中用水總量是考核該制度落實情況的重要標(biāo)準(zhǔn)之一。隨著經(jīng)濟的迅速發(fā)展,工業(yè)規(guī)模的逐漸擴大,工業(yè)用水占據(jù)總用水量比重逐步提升。在水資源消耗中,工業(yè)用水僅次于農(nóng)業(yè)用水,用水量居第二位[2]。工業(yè)用水量核算在水資源管理考核工作中起著舉足輕重的作用,因此掌握工業(yè)行業(yè)的用水狀況,加強對工業(yè)行業(yè)用水量測算和統(tǒng)計的研究非常必要。

    為了應(yīng)對日益嚴(yán)重的水資源危機,人類在重視用水量統(tǒng)計的同時還開展了一系列用水量預(yù)測工作。美國在1956年首次進(jìn)行用水量預(yù)測,日本、英國、法國等國家也逐漸開展用水量預(yù)測工作。在預(yù)測方法上,可以分為兩大類:一類是基于工業(yè)用水歷史數(shù)據(jù)進(jìn)行預(yù)測,如徐洪福等[3]、黃勝[4]和郭法強[5]對灰色預(yù)測模型做了大量的研究,澳大利亞的Zhou等[6-7]利用時間序列預(yù)測模型分別預(yù)測了日用水量和時用水量,Levi等[8]建立了逐步回歸模型,向平等[9]建立了BP神經(jīng)網(wǎng)絡(luò)模型;另一類是基于影響因子進(jìn)行預(yù)測,其前提是具備一些詳細(xì)真實的經(jīng)濟和用水統(tǒng)計等相關(guān)資料,包括從業(yè)人員數(shù)量、產(chǎn)值、資產(chǎn)余額、取用水量、重復(fù)利用率等,預(yù)測的結(jié)果精度相對較高,如徐得潛[10]和翟春健[11]將生產(chǎn)函數(shù)引入城市用水預(yù)測,Leonid等[12]于1993建立了模式識別模型。

    在用水統(tǒng)計方面,相較于美國、加拿大以及歐盟等發(fā)達(dá)國家,我國水資源統(tǒng)計體系較為復(fù)雜,除了專門負(fù)責(zé)統(tǒng)計水利數(shù)據(jù)的統(tǒng)計部門外,還有在職責(zé)范圍內(nèi)進(jìn)行水利信息統(tǒng)計的其他部門,導(dǎo)致不同部門之間的信息不能進(jìn)行及時準(zhǔn)確地共享。當(dāng)前我國工業(yè)企業(yè)用水信息的獲取主要通過水利部水利綜合統(tǒng)計報表制度、工業(yè)統(tǒng)計報表制度以及水利普查3種途徑,統(tǒng)計口徑上往往存在不統(tǒng)一現(xiàn)象。目前,各省市及各行業(yè)一般依據(jù)典型企業(yè)用水平衡的數(shù)據(jù)進(jìn)行用水測算,但各省級行業(yè)的工作基礎(chǔ)和投入不同,導(dǎo)致行業(yè)用水統(tǒng)計數(shù)據(jù)存疑,需進(jìn)一步對相應(yīng)數(shù)據(jù)進(jìn)行梳理、校核和評估。

    筆者在不考慮利用其他經(jīng)濟統(tǒng)計資料的前提下,綜合考慮現(xiàn)有用水統(tǒng)計的工作基礎(chǔ)和數(shù)據(jù),討論了工業(yè)用水趨勢外推法與當(dāng)前計量及統(tǒng)計方法,提出了利用預(yù)測方法和計量結(jié)合抽樣的統(tǒng)計方法對工業(yè)用水進(jìn)行測算和校核的思路,以提高評估工業(yè)用水變化趨勢準(zhǔn)確率。以浙江省湖州市為例,利用一些常用的基于歷史數(shù)據(jù)法對用水量進(jìn)行預(yù)測并分析其適用性,為工業(yè)用水量測算提供依據(jù);然后利用計量結(jié)合抽樣的統(tǒng)計方法利用湖州市2011年水利普查數(shù)據(jù)對2011年工業(yè)用水總量進(jìn)行測算,并與利用2003—2010年數(shù)據(jù)預(yù)測的結(jié)果進(jìn)行比較分析。

    1 基于歷史數(shù)據(jù)的工業(yè)用水量測算方法回顧和適應(yīng)性分析

    1.1 方法回顧

    依據(jù)歷史觀察數(shù)據(jù)進(jìn)行用水量預(yù)測的方法主要有解釋性預(yù)測法和時間序列分析法。解釋性預(yù)測法主要包括回歸分析等;時間序列分析方法主要包括灰色預(yù)測、趨勢外推、指數(shù)平滑、自回歸(AR)、移動平均(MA)、自回歸移動平均(ARMA)及神經(jīng)網(wǎng)絡(luò)等方法[13]。時間序列預(yù)測指利用按照時間順序進(jìn)行觀測或記錄的過去至現(xiàn)在的用水量數(shù)據(jù)排列起來,不考慮外部的影響因素,比如天氣、價格等,找到記錄或觀測的數(shù)據(jù)和時間變化之間的關(guān)系,并延展所呈現(xiàn)的變化趨勢,以預(yù)測未來的用水情況[14]。

    趨勢外推法是基于歷史統(tǒng)計數(shù)據(jù),根據(jù)過去和現(xiàn)在的發(fā)展趨勢得到系列變化規(guī)律,探索研究對象發(fā)展變化相對于時間的函數(shù)關(guān)系進(jìn)行擬合并預(yù)測長期的變化趨勢的一種方法,主要應(yīng)用在科技、經(jīng)濟和社會發(fā)展的預(yù)測。Wan等[15]提出了趨勢外推模型對居民年需水量進(jìn)行預(yù)測;Pradhan[16]對該法進(jìn)行定義并預(yù)測了生活用水趨勢;Grouillet等[17]認(rèn)為趨勢外推法預(yù)測結(jié)果的準(zhǔn)確性很大程度上是由原始數(shù)據(jù)決定的。

    指數(shù)平滑法主要是通過對整個時間序列給予不同的權(quán)重而進(jìn)行加權(quán)平均的一種方法,在中期和短期預(yù)測中應(yīng)用較多。指數(shù)平滑法在商業(yè)、采礦業(yè)、水文學(xué)等方面都有廣泛應(yīng)用,適用于時間序列長期趨勢變動和水平變動實物的預(yù)測。Mun[18]認(rèn)為指數(shù)平滑法是一種比較簡單又比較重要的時間序列預(yù)測法,在預(yù)測方法中較為常用。在水文學(xué)方面,Aly等[19]和Caiado[20]利用季節(jié)性指數(shù)平滑法預(yù)測了城市用水量;姜彥立等[21]對季節(jié)性平滑指數(shù)法進(jìn)行了改進(jìn)處理并預(yù)測了北方地區(qū)的用水量;楊光燦等[22]在日用水量預(yù)測方面采用了改進(jìn)單指數(shù)平滑法;山紅翠等[23]采用該法預(yù)測了在規(guī)劃水平年的城鎮(zhèn)公共需水量;在采礦業(yè)方面,侯曉東等[24]為實現(xiàn)礦井防塵用水量的準(zhǔn)確預(yù)測,建立了自動調(diào)整平滑參數(shù)的單指數(shù)平滑法預(yù)測模型。

    灰色系統(tǒng)是指信息不完全的系統(tǒng)?;疑A(yù)測就是對含有不確定因素的系統(tǒng),基于過去和現(xiàn)在已知的及不確定的信息建立一個從過去延伸至未來的GM模型,進(jìn)而預(yù)測系統(tǒng)未來發(fā)展情況的一種方法[25]。通過分析各類用水指標(biāo)的影響因素可知,部分指標(biāo)的影響機制仍不明確,存在部分信息已知、部分未知的現(xiàn)象,具有“灰色”特征。利用灰色預(yù)測模型,在用水量預(yù)測方面已有很多研究成果。利用傳統(tǒng)的灰色預(yù)測模型,徐洪福等[3]和張鑫等[26]對城市中長期用水量進(jìn)行了預(yù)測,預(yù)測結(jié)果精度較高;向宇等[27]短期預(yù)測了全國用水量,呈現(xiàn)了模型的實用價值;梁學(xué)玉[28]和王春超等[29]在模型改進(jìn)方面也進(jìn)行了研究,提出了組合灰色預(yù)測模型以及多變量灰色模型?;疑P妥钌賰H需4個數(shù)據(jù)就可以進(jìn)行預(yù)測,可應(yīng)用在數(shù)據(jù)收集難度大成本高,數(shù)據(jù)十分少的情況。

    自回歸移動平均模型是在1976年由美國統(tǒng)計學(xué)家Box和英國統(tǒng)計學(xué)家Jenkins提出的,此后該模型在供、用水預(yù)測方面有著廣泛的應(yīng)用[30]。黃瀟瑩等[31]預(yù)測了成都市2009年各季度的城市供水量;孫平等[32]以深圳市水廠為例,提出了水廠逐日需水量過程預(yù)測方法;張忠國[33]提出了將ARIMA與ANN相結(jié)合的需水量預(yù)測模型,研究表明,組合模型預(yù)測精度更高。研究顯示,歷史數(shù)據(jù)越多,此法預(yù)測精度越高,但實際供用水量數(shù)據(jù)收集難度大,所得數(shù)據(jù)準(zhǔn)確率低。

    人工神經(jīng)網(wǎng)絡(luò)在用水量預(yù)測方面的應(yīng)用始于20世紀(jì)90年代,因其非線性逼近能力強、可以進(jìn)行大規(guī)模并行計算等特點受到學(xué)者的關(guān)注。周天佑等[34]利用反向傳播神經(jīng)網(wǎng)絡(luò)預(yù)測了日用水量;楊艷等[35]利用小波神經(jīng)網(wǎng)絡(luò)對城市用水量做了長期預(yù)測;孫月峰等[36]利用模糊神經(jīng)網(wǎng)絡(luò)預(yù)測了城市需水量;錢光興[37]利用RBP和泛化回歸神經(jīng)網(wǎng)絡(luò)對城市供水量進(jìn)行了預(yù)測,預(yù)測效果較好。

    1.2 各預(yù)測方法的適用性分析及方法選取

    工業(yè)用水測算時需考慮測算成本和測算精度兩個基本因素。測算成本主要由數(shù)據(jù)的獲得方式及獲取難易程度決定。傳統(tǒng)的預(yù)測方法對于數(shù)據(jù)獲取的要求不高,通常只需綜合性的數(shù)據(jù),不會更多地考慮工業(yè)產(chǎn)值、產(chǎn)業(yè)結(jié)構(gòu)和能源消耗等方面,因此測算成本較低。傳統(tǒng)預(yù)測方法的測算精度主要依賴于這些綜合性數(shù)據(jù)自身的趨勢變化,若在預(yù)測期間工業(yè)用水?dāng)?shù)據(jù)由于產(chǎn)業(yè)結(jié)構(gòu)調(diào)整或節(jié)水技術(shù)提升等方面的原因出現(xiàn)了較大波動,預(yù)測精度會受到較大的影響。因此,傳統(tǒng)預(yù)測方法只是按照歷史數(shù)據(jù)的變化規(guī)律進(jìn)行預(yù)判,并不能全面真實地測算工業(yè)用水情況。但許多學(xué)者的研究表明,一般情況下,測算精度不低于85%,或者說測算誤差小于等于15%,這樣的模型或方法都是可行的[38]。因此,在精度得到保障的前提下,筆者認(rèn)為傳統(tǒng)方法的測算是可行的。

    表1詳細(xì)對比了5種傳統(tǒng)預(yù)測方法,這5種預(yù)測方法都可以對中短期用水進(jìn)行預(yù)測,但是各方法的適用條件和特點存在較大差異。首先,ARMA和神經(jīng)網(wǎng)絡(luò)可以較好地處理非線性關(guān)系,預(yù)測精確度較高,但是對歷史數(shù)據(jù)要求較高,需要較長的時間序列數(shù)據(jù)。而其他3種方法雖然在非線性關(guān)系處理和預(yù)測精確度上不如上述兩者,但可以在數(shù)據(jù)相對匱乏的狀況下實現(xiàn)預(yù)測。從實際角度出發(fā),目前我國工業(yè)用水統(tǒng)計方面仍然存在著統(tǒng)計不夠完善和數(shù)據(jù)不足的問題,因此趨勢外推、指數(shù)平滑和灰色預(yù)測方法相對更加適用于工業(yè)企業(yè)用水測算。隨著統(tǒng)計的規(guī)范化和數(shù)據(jù)的積累,ARMA、神經(jīng)網(wǎng)絡(luò)方法甚至大數(shù)據(jù)和云計算技術(shù)將更多地應(yīng)用在用水測算領(lǐng)域。對趨勢外推、指數(shù)平滑和灰色預(yù)測方法進(jìn)一步比較,我們發(fā)現(xiàn)趨勢外推法可以實現(xiàn)線性和非線性關(guān)系的處理,并能保證一定的精確度,實用性和靈活性更強,所以本文選取趨勢外推法作為工業(yè)用水量測算的傳統(tǒng)預(yù)測方法。

    表1 幾種傳統(tǒng)預(yù)測方法適用性比較

    2 計量結(jié)合抽樣的工業(yè)用水量測算方法

    工業(yè)用水主要包括火(核)電工業(yè)用水及一般工業(yè)用水。在工業(yè)用水統(tǒng)計中,可以選擇將工業(yè)用水戶按規(guī)模以上和規(guī)模以下用水進(jìn)行分層計量和抽樣,即對規(guī)模以上用水戶采用完全計量的方法獲得用水?dāng)?shù)據(jù),對規(guī)模以下用水戶按統(tǒng)計精度要求進(jìn)行抽樣計量,使得工業(yè)用水統(tǒng)計結(jié)果的精度能夠滿足工業(yè)用水管理和評估所需。這種計量結(jié)合抽樣的方法的精度和成本取決于抽樣計量樣本數(shù)量和完全計量的用水戶規(guī)模。

    2.1 抽樣計量樣本量的確定

    由于企業(yè)工業(yè)用水量分布復(fù)雜,受行業(yè)結(jié)構(gòu)、規(guī)模特征和地域經(jīng)濟發(fā)展水平等因素的影響,無法預(yù)知工業(yè)企業(yè)用水總體的概率分布特征,因此可利用蒙特卡羅方法按不同樣本數(shù)量進(jìn)行隨機抽樣,從而確定較為經(jīng)濟可靠的計量樣本數(shù)量。

    筆者以浙江省湖州市用水規(guī)模10萬m3以下工業(yè)企業(yè)為例,利用Bootstrap方法對不同量級的樣本進(jìn)行重復(fù)抽樣,并根據(jù)重復(fù)抽樣的平均誤差確定抽樣比。Bootstrap重復(fù)抽樣方法如下:

    (1)以200為樣本量單位,抽取的樣本量分別為200,400,600,…,3 800,4 000,共20種不同類型的樣本量。

    (2)每一種樣本量,按照簡單不重復(fù)隨機抽樣的方式,每種方案均模擬100次,以減少樣本隨機性偏差。

    (3)每一種抽樣方式得到工業(yè)用水統(tǒng)計評估量、估計量方差和相對估計誤差,取平均值作為評估結(jié)果。

    由簡單隨機抽樣方案的估計結(jié)果可得樣本量與誤差的關(guān)系及樣本量與估計量方差的關(guān)系,如圖1和圖2所示。

    圖1 樣本量與誤差的關(guān)系

    圖2 樣本量與估計量方差的關(guān)系

    通過對圖1的樣本量和誤差的關(guān)系分析,不難發(fā)現(xiàn)隨著樣本量的增加,相對誤差逐步減少,但是減少的趨勢逐步放緩。具體來講,若用水量的測算精度要保證15%的精度要求,樣本量僅需要在400個左右,抽樣比不能低于7%;若測算精度要達(dá)到10%的精度要求,則樣本量規(guī)模要擴大4倍左右,即抽樣比不能低于27%。以上說明當(dāng)樣本量達(dá)到一定程度之后,樣本的規(guī)模對相對誤差的減少的貢獻(xiàn)并不顯著,然而樣本規(guī)模的增加意味著增加調(diào)查成本。從估計量方差也可以看出,當(dāng)樣本量達(dá)到400個左右的時候,其樣本波動性出現(xiàn)了大幅的降低。因此,考慮預(yù)測精度和調(diào)查成本之間的矛盾,7%的樣本抽樣比能夠滿足工業(yè)用水量的測算精度要求。

    2.2 一般工業(yè)用水分層標(biāo)準(zhǔn)的確定

    通常,工業(yè)用水大戶的數(shù)量不多,但其用水占比較大,對用水測量精度的影響十分明顯,所以為保障工業(yè)用水統(tǒng)計的精度,將規(guī)模以上用水工業(yè)企業(yè)全部納入用水監(jiān)控計量范圍是十分必要的。那么,工業(yè)企業(yè)用水分層標(biāo)準(zhǔn)的設(shè)定也決定著工業(yè)用水測算精度的高低。以浙江省湖州市2011年的工業(yè)企業(yè)用水為例,年取水量50萬m3及以上的企業(yè)共44家,所計量水量為5 901.14萬m3,占到總用水量的18.59%;10萬m3以上的企業(yè)共有252家,計量水量為1 0744.75萬m3,占到總用水量的33.75%。

    假設(shè)分別以用水量10萬m3和50萬m3作為工業(yè)用水分層標(biāo)準(zhǔn),根據(jù)Bootstrap重復(fù)抽樣的方法比較兩種不同分層標(biāo)準(zhǔn)條件下規(guī)模以下樣本量與誤差的關(guān)系,如下圖3所示。

    圖3 規(guī)模以下樣本量與誤差的關(guān)系

    通過對比發(fā)現(xiàn),年取水量10萬m3作為分層標(biāo)準(zhǔn)的計量精度要顯著高于年取水量50萬m3。所以,對于湖州市而言,將年取水量10萬m3及以上的工業(yè)企業(yè)全部納入用水監(jiān)控計量范圍,會有效地保障工業(yè)用水統(tǒng)計的精度。

    通過對工業(yè)用水量分層標(biāo)準(zhǔn)的確定,分別對規(guī)模以上的和規(guī)模以下的工業(yè)企業(yè)采用直接計量和以一定抽樣比進(jìn)行隨機抽樣的方法測算工業(yè)用水量。這種分層計量和隨機抽樣調(diào)查相結(jié)合的方法可以綜合考慮工業(yè)用水統(tǒng)計的精度要求及統(tǒng)計成本;且較傳統(tǒng)預(yù)測方法可以更加全面地掌握工業(yè)企業(yè)用水動態(tài),特別是對用水大戶;同時,相較于全面普查,該方法的應(yīng)用成本更低,且簡單易操作。

    2.3 對比分析兩類測算方法

    為了更加準(zhǔn)確地判斷分層計量和隨機抽樣相結(jié)合測算方法的精確度,筆者將其與傳統(tǒng)預(yù)測方法相比較。在傳統(tǒng)方法中,選取上文中工業(yè)用水預(yù)測效果最好的趨勢外推法;以浙江湖州市工業(yè)用水為研究對象,以湖州市2005—2010年工業(yè)用水量數(shù)據(jù)及湖州市2011年水利普查信息作為湖州市2011年工業(yè)用水量測算的基礎(chǔ)數(shù)據(jù)。根據(jù)湖州市2012年水資源公報的內(nèi)容,全市規(guī)模以上工業(yè)和規(guī)模以下工業(yè)用水量分別為2.21億m3和0.96億m3,即一般工業(yè)用水量為3.17億m3,火(核)電年取水量為0.217億m3,則工業(yè)用水總量為3.387億m3。本文規(guī)定測算的精確度以相對誤差為標(biāo)準(zhǔn),標(biāo)準(zhǔn)誤差計算如下:

    相對誤差

    利用趨勢外推法對湖州市2005—2010年工業(yè)用水量數(shù)據(jù)進(jìn)行計算,得到了一般工業(yè)用水量線性變化趨勢y=-0.374 2t+5.486 6,t∈[0,7],2011年湖州市一般工業(yè)用水量的估計值為2.87億m3,一般工業(yè)用水預(yù)測相對誤差為9.5%。

    同理,利用工業(yè)用水總量線性變化趨勢y=-0.366 1t+5.650 7,t∈[0,7],測算得到2011年湖州市工業(yè)用水總量的估計值為3.09億m3,那么相對誤差為8.8%。用水總量的估計值和其相對誤差作為比較分析的基準(zhǔn)值。將根據(jù)分層計量和隨機抽樣相結(jié)合測算方法對工業(yè)用水量再次核算,與基準(zhǔn)值相比較。

    根據(jù)分層計量和隨機抽樣相結(jié)合的測算方法對2011年湖州市5809家工業(yè)企業(yè)的用水量進(jìn)行測算,其中規(guī)模以上有252家企業(yè),其用水量直接計量的結(jié)果為1.07億m3;采用隨機抽樣的方法估計得到規(guī)模以下工業(yè)企業(yè)平均用水量為3.31萬m3,則規(guī)模以下企業(yè)的用水量1.84億m3。因此,湖州市一般工業(yè)用水量的估計值為2.91億m3。一般工業(yè)用水相對誤差為8.2%。

    直接計量2011年火(核)電年用水量為0.222億m3,則計量結(jié)合抽樣所得總用水量為3.132億m3??傆盟肯鄬φ`差為7.6%。

    表2描述和比較了上面兩種方法的測算結(jié)果。

    表2 各測算方法相對誤差比較

    從表2的結(jié)果比較可以看出,計量結(jié)合抽樣測算結(jié)果相對誤差更小,精度明顯優(yōu)于傳統(tǒng)的趨勢外推法。這主要因為簡單隨機抽樣使得一定規(guī)模以上水量直接計量,從而不產(chǎn)生統(tǒng)計誤差,主要誤差來源于數(shù)量眾多且分散、用水比重相對較小的規(guī)模以下用水隨機抽樣統(tǒng)計部分,此法對有一定計量基礎(chǔ)的一般工業(yè)用水量的評估較為適宜。所以,為了更加準(zhǔn)確地了解一個地區(qū)的工業(yè)用水量,計量結(jié)合抽樣測算的方法可以更加積極地得到推廣。同時,由于當(dāng)前階段我國的工業(yè)用水計量規(guī)章制度仍不完善,全面的工業(yè)用水計量還無法實現(xiàn),因此重點監(jiān)控,適量抽樣調(diào)查的工業(yè)用水量測算方法不僅成本相對便宜,而且也更加簡化、合理、實用。

    但是,應(yīng)該注意到,由于簡單隨機抽樣具有很大的隨機性,在實際操作過程中,通過反復(fù)地重復(fù)抽樣以消除偏差是很難操作的,因此還是需要采用傳統(tǒng)的預(yù)測方法校核計量結(jié)合抽樣測算的結(jié)果。如果這兩種方法計算的結(jié)果出現(xiàn)較大的差異,我們還需要進(jìn)一步考察和確定工業(yè)用水量的變化。

    3 結(jié) 語

    本文針對工業(yè)用水量測算問題進(jìn)行了初步的研究和探討,選取了基于歷史數(shù)據(jù)和計量結(jié)合抽樣的方法用于湖州市2011年的用水量測算。通過對比分析基于歷史數(shù)據(jù)的3種測算結(jié)果,指出3種預(yù)測方法的特點和適用性;并對比簡單隨機抽樣及分層抽樣的相對誤差和估計量方差,結(jié)果表明按取水量(10萬m3)分層抽樣方案能滿足測算要求。通過對工業(yè)用水量兩種測算方法的詳細(xì)研究,提出“以計量結(jié)合抽樣的工業(yè)用水量測算方法為主,趨勢外推法進(jìn)行對比校核”的工業(yè)用水量測算方法。測算結(jié)果表明兩種測算方法得到的一般工業(yè)用水量和工業(yè)用水總量非常相近,測算相對誤差均在10%以內(nèi),可達(dá)到精度要求。

    [ 1 ] 卞錦宇,劉恒,耿雷華, 等. 基于隨機前沿生產(chǎn)函數(shù)的我國工業(yè)用水效率影響因素研究[J]. 水利經(jīng)濟,2014,32(5):6-9,71.

    [ 2 ] 賈邵鳳. 水資源經(jīng)濟學(xué)[M]. 北京: 中國水利水電出版社, 2006: 146.

    [ 3 ] 徐洪福,袁一星,趙洪賓. 灰色預(yù)測模型在年用水量預(yù)測中的應(yīng)用[J]. 哈爾濱建筑大學(xué)學(xué)報, 2001, 34(4):51-54.

    [ 4 ] 黃勝. 灰色系統(tǒng)與多元逐步回歸耦合模型在年用水量預(yù)測中的應(yīng)用[J]. 節(jié)水灌溉, 2007(5):38-40.

    [ 5 ] 郭法強. 基于灰色GM(1,1)模型的城市工業(yè)用水量預(yù)測[J]. 地下水, 2010,32(3):108-110.

    [ 6 ] ZHOU S L, MCMAHON T A, WALTON A, et al. Forecasting daily urban water demand: a case study of Melbourne[J]. Journal of Hydrology, 2000, 236(3-4):153-164.

    [ 7 ] ZHOU S L, MCMAHON T A, WALTON A, et al. Forecasting operational demand for an urban water supply zone[J]. Journal of Hydrology, 2002, 259(1-4):189-202.

    [ 8 ] BREKKE L, TAKAICHI L. Suburban water demand modeling using stepwise regression[J]. Journal American Water Works Association, 2002, 94(10):65-75.

    [ 9 ] 向平,張蒙,張智,等. 基于BP神經(jīng)網(wǎng)絡(luò)的城市時用水量分時段預(yù)測模型[J]. 中南大學(xué)學(xué)報(自然科學(xué)版),2012, 43(8):3320-3324.

    [10] 徐得潛. 城市用水量預(yù)測[J]. 水資源與水工程學(xué)報, 1992, 3(3):73-76.

    [11] 翟春健,張宏偉,王亮. 工業(yè)生產(chǎn)函數(shù)法用于城市工業(yè)用水量預(yù)測的研究[J]. 天津工業(yè)大學(xué)學(xué)報, 2009, 28(5):79-81.

    [12] SHVARTSER L, SHAMIR U, FELDMAN M. Forecasting hourly water demand by pattern recognition approach[J]. Journal of Water Resources Planning & Management, 1993, 119(6):611-627.

    [13] DONKOR E A, MAZZUCHI T A, SOYER R, et al. Urban Water demand forecasting: review of methods and models[J]. Journal of Water Resources Planning & Management, 2014, 140(2):146-159.

    [14] BILLINGS R B, JONES, C V. Forecasting urban water demand[M]. Washington:American Water Works Association, 1996: 154.

    [15] WAN S N, KUCZERA G. Incorporating demand uncertainty in water supply headworks simulation[J]. Water Resources Research, 1993, 29(2):469-478.

    [16] BABEL M S, GUPTA A D, PRADHAN P. Multivariate econometric approach for domestic water demand modeling: an application to Kathmandu, Nepal[J]. Water Resource Management, 2007,21:573-589.

    [17] GROUILLE B, FABRE J, RUELLAND D, et al. Historical reconstruction and 2050 projections of water demand under anthropogenic and climate changes in two contrasted Mediterranean catchments[J]. Journal of Hydrology, 2015, 522:684-696.

    [18] MUN J. Modeling risk: applying Monte Carlo risk simulation, strategic real options,stochastic forecasting, and portfolio optimization[J]. Contact stress analysis, 1990,144:647-653.

    [19] ALY A H, WANAKULE N. Short-term forecasting for urban water consumption[J]. Journal of Water Resources Planning & Management, 2004, 130(5): 405-410.

    [20] CAIADO J. Performance of combined double seasonal univariate time series models for forecasting water demand[J]. Journal of Hydrologic Engineering, 2009, 15(3): 215-222.

    [21] 姜彥立, 周新華. 季節(jié)性指數(shù)平滑法在城市用水量預(yù)測中的應(yīng)用研究[J]. 科學(xué)技術(shù)與工程, 2007, 7(8):1666-1669.

    [22] 楊光燦,張鑫,史浩.改進(jìn)單指數(shù)平滑法在某市日用水量預(yù)測中的應(yīng)用[J]. 黑龍江交通科技,2013(11):156,158.

    [23] 山紅翠,盛東,郭文娟,等.湖南省生活需水量預(yù)測[J]. 水資源研究, 2016,5(1):52-58.

    [24] 侯曉東,司樹強,沈宏鵬. 基于自適應(yīng)單指數(shù)平滑法的礦井防塵用水量預(yù)測[J]. 煤炭與化工, 2015,38(01):61-64.

    [25] 鄧聚龍. 灰色系統(tǒng)基本方法[M]. 武漢:華中科技大學(xué)出版社, 2005.

    [26] 張鑫,趙晴,史成波. 灰色模型在中長期用水量預(yù)測中的應(yīng)用[J]. 交通科技與經(jīng)濟, 2011, 13(4):113-115.

    [27] 向宇,吳琴. 改進(jìn)的灰色預(yù)測模型在全國用水量預(yù)測中的應(yīng)用[J]. 湖北民族學(xué)院學(xué)報(自然科學(xué)版),2015, 33(1):25-30.

    [28] 梁學(xué)玉,張鑫,孫天青. 組合灰色預(yù)測模型在城市用水量預(yù)測中的應(yīng)用[J]. 人民黃河, 2010, 32(4):79-80.

    [29] 王春超,王麗萍,曹云慧,等. 改進(jìn)多變量灰色模型在城市用水量預(yù)測中的應(yīng)用[J]. 水電能源科學(xué), 2013, 31(2):27-29.

    [30] BOX G E P, JENKINS G M. Time series analysis: forecasting and control[M]. Holden:Holden-Day Incorporated, 1990:199-201.

    [31] 黃瀟瑩,張健. ARIMA模型在成都市供水量預(yù)測中的應(yīng)用[J]. 水資源與水工程學(xué)報, 2010, 21(4):111-113.

    [32] 孫平,王麗萍,陳凱,等. 基于時間序列模型ARMA的水廠逐日需水量過程預(yù)測方法[J].中國農(nóng)村水利水電,2013(11):139-142.

    [33] 張忠國. 基于ARIMA與ANN組合模型的需水量預(yù)測方法研究[J]. 地下水, 2015, 37(6):1-4.

    [34] 周天佐,羅大庸. MCPSO優(yōu)化BP-NN在日用水量預(yù)測中的應(yīng)用[J]. 微計算機信息, 2009, 25(9):151-153.

    [35] 楊艷,李靖,馬顯瑩,等. 基于小波神經(jīng)網(wǎng)絡(luò)的城市用水量長期預(yù)測研究[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報, 2010, 25(2):272-276.

    [36] 孫月峰,閆雅飛,張表志,等. 基于T-S模型的模糊神經(jīng)網(wǎng)絡(luò)城市需水量預(yù)測方法研究[J]. 安全與環(huán)境學(xué)報,2013, 13(2):136-139.

    [37] 錢光興,崔東文. RBF與GRNN神經(jīng)網(wǎng)絡(luò)模型在城市需水預(yù)測中的應(yīng)用[J]. 水資源與水工程學(xué)報, 2012, 23(5):148-152.

    [38] 高海英. 神經(jīng)網(wǎng)絡(luò)在陜西省工業(yè)用水量預(yù)測中的應(yīng)用[D].西安: 西安建筑科技大學(xué), 2007.

    10.3880/j.issn.1003-9511.2017.01.010

    周申蓓(1974—),男,副教授,博士研究生,主要從事水資源管理研究。E-mail:shenbei@hhu.edu.cn

    鄭士鵬(1991—),男,碩士研究生,主要從事水資源管理研究。E-mail:13777682662@163.com

    TU991.31

    A

    1003-9511(2017)01-0040-05

    2016-04-16 編輯:陳玉國)

    猜你喜歡
    方法模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    學(xué)習(xí)方法
    可能是方法不對
    3D打印中的模型分割與打包
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    欧美大码av| 校园春色视频在线观看| 亚洲国产日韩欧美精品在线观看 | 九九在线视频观看精品| 日本 av在线| av在线天堂中文字幕| 午夜精品在线福利| 免费观看精品视频网站| 国产极品精品免费视频能看的| 国产精品av视频在线免费观看| av在线天堂中文字幕| 无人区码免费观看不卡| 国产亚洲精品av在线| 三级国产精品欧美在线观看 | 欧美一级a爱片免费观看看| av天堂中文字幕网| 美女黄网站色视频| 欧美色视频一区免费| 精品久久久久久久久久久久久| 亚洲国产精品合色在线| 国产主播在线观看一区二区| 俺也久久电影网| 日韩欧美在线二视频| 日本熟妇午夜| 国产不卡一卡二| 国产黄色小视频在线观看| x7x7x7水蜜桃| 国产成+人综合+亚洲专区| 此物有八面人人有两片| 俄罗斯特黄特色一大片| 成熟少妇高潮喷水视频| 三级男女做爰猛烈吃奶摸视频| 国产精品香港三级国产av潘金莲| 午夜福利在线观看吧| 成熟少妇高潮喷水视频| 小说图片视频综合网站| 深夜精品福利| 精品熟女少妇八av免费久了| 99国产精品99久久久久| 村上凉子中文字幕在线| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产99精品国产亚洲性色| 久久精品国产综合久久久| 久久精品人妻少妇| 香蕉丝袜av| 亚洲人成伊人成综合网2020| 日韩国内少妇激情av| 国产野战对白在线观看| 亚洲国产精品成人综合色| 欧美一级毛片孕妇| 精品不卡国产一区二区三区| 免费高清视频大片| 非洲黑人性xxxx精品又粗又长| 一级毛片精品| 午夜视频精品福利| 俄罗斯特黄特色一大片| 蜜桃久久精品国产亚洲av| 国产黄片美女视频| 中文字幕人妻丝袜一区二区| 99热6这里只有精品| 99久久精品一区二区三区| 国产精品综合久久久久久久免费| 俄罗斯特黄特色一大片| 欧美中文综合在线视频| 精品国产亚洲在线| 成年女人永久免费观看视频| 美女高潮的动态| 国产亚洲精品av在线| 午夜影院日韩av| 偷拍熟女少妇极品色| 桃色一区二区三区在线观看| 婷婷亚洲欧美| 欧美黄色片欧美黄色片| 99在线视频只有这里精品首页| 精品国产乱子伦一区二区三区| 国产99白浆流出| 亚洲av成人av| 非洲黑人性xxxx精品又粗又长| 亚洲成av人片在线播放无| 亚洲av成人不卡在线观看播放网| 黄频高清免费视频| 国产精品 欧美亚洲| 国产亚洲av高清不卡| 美女黄网站色视频| 99久久国产精品久久久| 亚洲欧美日韩东京热| 97人妻精品一区二区三区麻豆| 亚洲欧美一区二区三区黑人| 久久午夜综合久久蜜桃| 18禁观看日本| 一夜夜www| 校园春色视频在线观看| 色哟哟哟哟哟哟| 久久香蕉国产精品| 国产成人欧美在线观看| 女同久久另类99精品国产91| 国产亚洲欧美在线一区二区| 精品日产1卡2卡| 免费无遮挡裸体视频| 久久久成人免费电影| 午夜福利视频1000在线观看| 窝窝影院91人妻| 好男人电影高清在线观看| 最近视频中文字幕2019在线8| 夜夜夜夜夜久久久久| 日本熟妇午夜| 亚洲精品色激情综合| 精品一区二区三区视频在线 | 此物有八面人人有两片| 日本a在线网址| 精品99又大又爽又粗少妇毛片 | 午夜激情欧美在线| 白带黄色成豆腐渣| 一二三四社区在线视频社区8| 狂野欧美白嫩少妇大欣赏| 久久精品国产综合久久久| 搡老岳熟女国产| 91久久精品国产一区二区成人 | 成人精品一区二区免费| 在线观看一区二区三区| 精品一区二区三区视频在线 | 色在线成人网| 男插女下体视频免费在线播放| 精品一区二区三区av网在线观看| 午夜福利在线观看免费完整高清在 | 亚洲专区字幕在线| av天堂中文字幕网| 免费看a级黄色片| 亚洲国产欧洲综合997久久,| а√天堂www在线а√下载| 别揉我奶头~嗯~啊~动态视频| 亚洲国产中文字幕在线视频| 欧美日本视频| 动漫黄色视频在线观看| 草草在线视频免费看| 欧美性猛交╳xxx乱大交人| 欧美乱妇无乱码| 法律面前人人平等表现在哪些方面| 好看av亚洲va欧美ⅴa在| 一区二区三区高清视频在线| 九色成人免费人妻av| 欧美乱妇无乱码| 黄色丝袜av网址大全| 少妇的逼水好多| АⅤ资源中文在线天堂| 真人一进一出gif抽搐免费| 亚洲成人久久爱视频| 麻豆一二三区av精品| 日本黄色视频三级网站网址| 免费av毛片视频| 一个人观看的视频www高清免费观看 | 99国产综合亚洲精品| 给我免费播放毛片高清在线观看| 成年女人毛片免费观看观看9| 国产精品野战在线观看| 精品国产超薄肉色丝袜足j| 男女下面进入的视频免费午夜| 香蕉丝袜av| 午夜激情欧美在线| 欧美最黄视频在线播放免费| 熟女电影av网| 一区福利在线观看| 欧美日韩乱码在线| 午夜激情福利司机影院| 国产精品亚洲一级av第二区| 日韩欧美免费精品| 亚洲欧洲精品一区二区精品久久久| 国产一区二区三区视频了| 黄色视频,在线免费观看| 校园春色视频在线观看| 亚洲av成人不卡在线观看播放网| 嫁个100分男人电影在线观看| 天堂av国产一区二区熟女人妻| 18禁黄网站禁片免费观看直播| 18禁裸乳无遮挡免费网站照片| 国产1区2区3区精品| 亚洲国产看品久久| 精品无人区乱码1区二区| 国产乱人伦免费视频| 亚洲av成人av| svipshipincom国产片| 国产高清视频在线播放一区| 老鸭窝网址在线观看| 中文字幕久久专区| 久久婷婷人人爽人人干人人爱| 在线观看日韩欧美| 在线观看免费视频日本深夜| 精品免费久久久久久久清纯| 99国产精品99久久久久| 99热只有精品国产| 蜜桃久久精品国产亚洲av| 午夜视频精品福利| 18禁裸乳无遮挡免费网站照片| 亚洲av成人一区二区三| 午夜免费观看网址| 成熟少妇高潮喷水视频| 国产男靠女视频免费网站| 丝袜人妻中文字幕| 人人妻人人看人人澡| avwww免费| 99精品久久久久人妻精品| 国产精品乱码一区二三区的特点| 欧美日韩中文字幕国产精品一区二区三区| 国产精品 欧美亚洲| 伊人久久大香线蕉亚洲五| 亚洲aⅴ乱码一区二区在线播放| 99热精品在线国产| 久久婷婷人人爽人人干人人爱| 黄色片一级片一级黄色片| 亚洲午夜精品一区,二区,三区| 久久久久久久久久黄片| 国产激情久久老熟女| 国产高清视频在线播放一区| 国产成人影院久久av| 国产亚洲精品av在线| 国产精品,欧美在线| 国产69精品久久久久777片 | 91麻豆av在线| 国产av在哪里看| 日韩高清综合在线| 亚洲国产精品合色在线| av在线天堂中文字幕| 国产 一区 欧美 日韩| 精品久久久久久,| 男人和女人高潮做爰伦理| 亚洲自偷自拍图片 自拍| 欧美黑人欧美精品刺激| 欧美中文综合在线视频| 国产精品久久久av美女十八| 亚洲精华国产精华精| 精品日产1卡2卡| 黄色 视频免费看| 哪里可以看免费的av片| 熟女人妻精品中文字幕| 久久精品91蜜桃| 一个人免费在线观看电影 | 一个人免费在线观看电影 | а√天堂www在线а√下载| 精品国内亚洲2022精品成人| 色吧在线观看| 欧美成人免费av一区二区三区| 午夜视频精品福利| 亚洲成av人片免费观看| 三级男女做爰猛烈吃奶摸视频| 黄色视频,在线免费观看| 99热这里只有是精品50| 国产伦精品一区二区三区视频9 | 亚洲人成网站高清观看| 国产精品久久久久久精品电影| e午夜精品久久久久久久| 亚洲中文字幕一区二区三区有码在线看 | 可以在线观看毛片的网站| 国产不卡一卡二| 很黄的视频免费| 欧美日韩黄片免| cao死你这个sao货| 国产亚洲精品av在线| 五月玫瑰六月丁香| 天堂影院成人在线观看| 男人舔女人下体高潮全视频| 免费看日本二区| 欧美一区二区精品小视频在线| 中亚洲国语对白在线视频| 国产午夜精品论理片| 丰满的人妻完整版| 久久天躁狠狠躁夜夜2o2o| 老司机午夜福利在线观看视频| 国产成+人综合+亚洲专区| 成人三级黄色视频| 高清毛片免费观看视频网站| 日韩大尺度精品在线看网址| 一区二区三区高清视频在线| 欧美一区二区国产精品久久精品| 国内精品久久久久久久电影| 动漫黄色视频在线观看| 99riav亚洲国产免费| 亚洲欧美日韩东京热| 国产一区二区激情短视频| 国产av一区在线观看免费| 精品久久久久久,| 国产美女午夜福利| 日韩中文字幕欧美一区二区| 中国美女看黄片| 成人永久免费在线观看视频| 极品教师在线免费播放| 啦啦啦韩国在线观看视频| 欧美午夜高清在线| tocl精华| 真人做人爱边吃奶动态| 日韩欧美精品v在线| 搡老妇女老女人老熟妇| 日韩欧美在线二视频| 色在线成人网| 国产亚洲欧美98| ponron亚洲| 国产亚洲精品久久久com| 国语自产精品视频在线第100页| 免费观看人在逋| 成人精品一区二区免费| 国产精品一区二区三区四区久久| 99热这里只有是精品50| 岛国视频午夜一区免费看| 99热这里只有精品一区 | 毛片女人毛片| 九九在线视频观看精品| 十八禁人妻一区二区| 国产毛片a区久久久久| 日韩欧美免费精品| 色综合婷婷激情| www.999成人在线观看| 国内久久婷婷六月综合欲色啪| 国产成人系列免费观看| 亚洲一区高清亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 欧美中文日本在线观看视频| 亚洲国产精品sss在线观看| 99久久国产精品久久久| 91字幕亚洲| 国产欧美日韩一区二区精品| 欧美日韩黄片免| 高清毛片免费观看视频网站| 一区二区三区激情视频| 毛片女人毛片| 三级国产精品欧美在线观看 | 国产av在哪里看| 熟女人妻精品中文字幕| 男人舔女人的私密视频| 非洲黑人性xxxx精品又粗又长| 日韩欧美 国产精品| 欧美乱妇无乱码| 久久香蕉精品热| 99久久无色码亚洲精品果冻| h日本视频在线播放| 亚洲欧美日韩卡通动漫| 国产毛片a区久久久久| 男女之事视频高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品自产拍在线观看55亚洲| 久久久久免费精品人妻一区二区| 久久中文字幕一级| 成人永久免费在线观看视频| 亚洲欧美一区二区三区黑人| 成人国产综合亚洲| 99久久无色码亚洲精品果冻| 一二三四在线观看免费中文在| 亚洲av熟女| 国产精品乱码一区二三区的特点| 99国产精品一区二区蜜桃av| 欧美精品啪啪一区二区三区| 成人午夜高清在线视频| h日本视频在线播放| 91久久精品国产一区二区成人 | 成年女人毛片免费观看观看9| 国产爱豆传媒在线观看| 亚洲成人久久爱视频| 午夜视频精品福利| 在线免费观看不下载黄p国产 | 国产精品久久电影中文字幕| 日本 av在线| 蜜桃久久精品国产亚洲av| 少妇丰满av| 熟女少妇亚洲综合色aaa.| 搡老岳熟女国产| 精品久久蜜臀av无| 午夜福利欧美成人| 日本三级黄在线观看| 在线观看一区二区三区| 久久久久精品国产欧美久久久| 国产真人三级小视频在线观看| 欧美中文日本在线观看视频| 少妇丰满av| 蜜桃久久精品国产亚洲av| 热99在线观看视频| 午夜久久久久精精品| 香蕉av资源在线| 天堂影院成人在线观看| 国产精品一区二区三区四区久久| 国产精品野战在线观看| 国产精品自产拍在线观看55亚洲| 免费在线观看视频国产中文字幕亚洲| 亚洲aⅴ乱码一区二区在线播放| 色老头精品视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品国产高清国产av| 亚洲欧美日韩卡通动漫| 又紧又爽又黄一区二区| 男人和女人高潮做爰伦理| 麻豆久久精品国产亚洲av| av女优亚洲男人天堂 | 成熟少妇高潮喷水视频| 久久久成人免费电影| 制服人妻中文乱码| 亚洲国产欧洲综合997久久,| 天堂av国产一区二区熟女人妻| 亚洲无线在线观看| 国产精品99久久久久久久久| 欧美zozozo另类| 老熟妇乱子伦视频在线观看| 老司机午夜十八禁免费视频| 日韩欧美精品v在线| 中出人妻视频一区二区| 亚洲片人在线观看| 亚洲天堂国产精品一区在线| 国产av不卡久久| 精品久久久久久久人妻蜜臀av| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品一区二区www| 一个人免费在线观看电影 | 毛片女人毛片| 伊人久久大香线蕉亚洲五| 日本撒尿小便嘘嘘汇集6| 免费看a级黄色片| 亚洲第一欧美日韩一区二区三区| 黄色女人牲交| 精品国产乱子伦一区二区三区| 国产精品九九99| 在线观看舔阴道视频| 亚洲人成网站在线播放欧美日韩| 男女视频在线观看网站免费| 九九在线视频观看精品| 99精品久久久久人妻精品| 波多野结衣高清无吗| 女人被狂操c到高潮| 男人舔女人下体高潮全视频| 九色国产91popny在线| 99在线人妻在线中文字幕| 小蜜桃在线观看免费完整版高清| 免费看美女性在线毛片视频| 国产乱人伦免费视频| 99riav亚洲国产免费| av在线蜜桃| 国产成人aa在线观看| 一进一出抽搐gif免费好疼| 免费搜索国产男女视频| netflix在线观看网站| 在线观看一区二区三区| 免费在线观看成人毛片| 午夜福利在线观看免费完整高清在 | 亚洲午夜精品一区,二区,三区| 色综合站精品国产| 国产伦一二天堂av在线观看| 国产成人av激情在线播放| 热99在线观看视频| 国产v大片淫在线免费观看| 草草在线视频免费看| 黑人欧美特级aaaaaa片| 欧美一区二区国产精品久久精品| 一级毛片高清免费大全| 波多野结衣高清作品| 欧美国产日韩亚洲一区| 亚洲国产精品sss在线观看| 给我免费播放毛片高清在线观看| 亚洲国产精品合色在线| 成人午夜高清在线视频| 小说图片视频综合网站| 熟女人妻精品中文字幕| 美女免费视频网站| 最近最新中文字幕大全电影3| 深夜精品福利| 在线国产一区二区在线| 免费人成视频x8x8入口观看| 国内精品美女久久久久久| 国产精品美女特级片免费视频播放器 | 天天躁日日操中文字幕| 国产一级毛片七仙女欲春2| 非洲黑人性xxxx精品又粗又长| 精华霜和精华液先用哪个| 一本精品99久久精品77| 成年女人永久免费观看视频| 757午夜福利合集在线观看| 高清毛片免费观看视频网站| 此物有八面人人有两片| 欧美成人一区二区免费高清观看 | 成年女人毛片免费观看观看9| 国产免费av片在线观看野外av| 午夜精品一区二区三区免费看| 嫩草影院精品99| 日韩精品中文字幕看吧| 亚洲五月婷婷丁香| 嫩草影视91久久| 国产v大片淫在线免费观看| 国产91精品成人一区二区三区| 国产精品野战在线观看| 国产成人精品无人区| 国产99白浆流出| 国产日本99.免费观看| 久久亚洲精品不卡| 免费在线观看日本一区| 午夜亚洲福利在线播放| 久久精品91蜜桃| 久久久久性生活片| 在线永久观看黄色视频| 国产欧美日韩精品一区二区| av视频在线观看入口| 亚洲精品在线美女| 美女高潮的动态| 免费av毛片视频| 一区二区三区国产精品乱码| 色播亚洲综合网| 悠悠久久av| 老司机午夜福利在线观看视频| 18禁国产床啪视频网站| 久久久久久久久中文| 99久久无色码亚洲精品果冻| ponron亚洲| 色综合欧美亚洲国产小说| 亚洲18禁久久av| 国产伦一二天堂av在线观看| 一本久久中文字幕| 国产野战对白在线观看| 观看免费一级毛片| 黄色视频,在线免费观看| 在线a可以看的网站| 无人区码免费观看不卡| 九九久久精品国产亚洲av麻豆 | 国产精品一及| 亚洲成av人片在线播放无| 久久精品91无色码中文字幕| 又粗又爽又猛毛片免费看| 欧美极品一区二区三区四区| 十八禁网站免费在线| 老司机深夜福利视频在线观看| 久久久久国内视频| 午夜免费成人在线视频| 制服人妻中文乱码| 久久中文字幕一级| 精品国产亚洲在线| 欧美丝袜亚洲另类 | 人妻丰满熟妇av一区二区三区| 91av网站免费观看| 亚洲中文字幕一区二区三区有码在线看 | 色哟哟哟哟哟哟| 又爽又黄无遮挡网站| 88av欧美| 免费看光身美女| 1024手机看黄色片| 亚洲国产看品久久| 国产伦精品一区二区三区四那| 欧美激情在线99| 国产爱豆传媒在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲无线在线观看| 欧美日韩综合久久久久久 | 亚洲专区中文字幕在线| 亚洲国产精品久久男人天堂| 露出奶头的视频| 国产精品亚洲一级av第二区| 97超级碰碰碰精品色视频在线观看| 国产免费av片在线观看野外av| 色综合婷婷激情| 日韩高清综合在线| 亚洲成人精品中文字幕电影| 欧美一级a爱片免费观看看| 91在线精品国自产拍蜜月 | 成年女人永久免费观看视频| 国产精品1区2区在线观看.| 精品无人区乱码1区二区| 国产欧美日韩一区二区三| 精品久久久久久久人妻蜜臀av| 亚洲精华国产精华精| 人人妻,人人澡人人爽秒播| 国产精品永久免费网站| 亚洲国产欧美人成| 一级黄色大片毛片| 欧美丝袜亚洲另类 | 特级一级黄色大片| 国产淫片久久久久久久久 | 亚洲中文字幕日韩| 免费看日本二区| 床上黄色一级片| 白带黄色成豆腐渣| 少妇熟女欧美另类| 黑人高潮一二区| 国内精品一区二区在线观看| 成人av在线播放网站| 久久亚洲精品不卡| 中国美白少妇内射xxxbb| 99久久精品国产国产毛片| av国产免费在线观看| 精品国内亚洲2022精品成人| 一个人免费在线观看电影| 寂寞人妻少妇视频99o| 欧美另类亚洲清纯唯美| 一级毛片电影观看 | 91在线精品国自产拍蜜月| 久久热精品热| 中文在线观看免费www的网站| 尤物成人国产欧美一区二区三区| 久久久久久久午夜电影| 欧美日韩综合久久久久久| 超碰97精品在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品合色在线| 在线免费观看不下载黄p国产| 国产精品,欧美在线| 亚洲最大成人中文| 少妇人妻精品综合一区二区| 最近最新中文字幕免费大全7| 伊人久久精品亚洲午夜| 狂野欧美白嫩少妇大欣赏| 亚洲一级一片aⅴ在线观看| 亚洲国产精品合色在线| 少妇熟女欧美另类| 国产男人的电影天堂91| 波野结衣二区三区在线| 韩国高清视频一区二区三区| 变态另类丝袜制服| 夫妻性生交免费视频一级片| 亚洲一区高清亚洲精品| 国产精品三级大全| 插逼视频在线观看| 成人亚洲精品av一区二区| 最近视频中文字幕2019在线8| 精品熟女少妇av免费看| 九九爱精品视频在线观看|