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

    利用實(shí)測值估算斷面泥沙年負(fù)荷的方法比較

    2017-02-22 07:23:39蘇靜君李敘勇中國科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室北京00085交通運(yùn)輸部環(huán)境保護(hù)中心北京0003
    中國環(huán)境科學(xué) 2017年1期
    關(guān)鍵詞:泥沙通量徑流

    蘇靜君,李敘勇*,吳 震(.中國科學(xué)院生態(tài)環(huán)境研究中心,城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,北京 00085;.交通運(yùn)輸部環(huán)境保護(hù)中心,北京 0003)

    利用實(shí)測值估算斷面泥沙年負(fù)荷的方法比較

    蘇靜君1,李敘勇1*,吳 震2(1.中國科學(xué)院生態(tài)環(huán)境研究中心,城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,北京 100085;2.交通運(yùn)輸部環(huán)境保護(hù)中心,北京 100013)

    主要針對非固定頻率采樣、水質(zhì)監(jiān)測高度集中在汛期的情形,基于河北承德柳河子流域 1995~2011年間逐日徑流及泥沙含量數(shù)據(jù),比較了平均值法、比例法和回歸法估算泥沙年負(fù)荷的表現(xiàn),并提出此情況下泥沙年負(fù)荷估算流程.研究結(jié)果表明,比例法、流量加權(quán)平均濃度與時段平均流量之積、以及流量-泥沙濃度相關(guān)關(guān)系滿足一定條件的回歸法均可以獲得比較理想的結(jié)果.比例法在樣品量稀疏(采樣頻率≤30%)的情況下表現(xiàn)更優(yōu).基于瞬時流量平均與瞬時濃度平均之積、以及瞬時通量平均的算法對徑流量的時均離散項(xiàng)比較敏感,在低采樣頻率下往往高估時段污染負(fù)荷,因而適合于水質(zhì)數(shù)據(jù)比較充足情況下的負(fù)荷估算.研究結(jié)果可為干旱半干旱地區(qū)具有脈沖式徑流特征的流域泥沙含量監(jiān)測頻率的設(shè)置及泥沙年負(fù)荷估算方法的選擇提供一定依據(jù).

    斷面泥沙負(fù)荷估算;比例法;脈沖式徑流特征

    準(zhǔn)確估算河流斷面或流域出口污染物負(fù)荷或通量對河流污染物總量控制及流域水環(huán)境管理具有重要意義,也是率定和驗(yàn)證流域污染負(fù)荷模型的重要基礎(chǔ)[1].依據(jù)水量和水質(zhì)實(shí)測資料計(jì)算污染負(fù)荷是估算的最基本方法,其結(jié)果常被認(rèn)為最接近“真實(shí)值”[2].然而眾多研究表明,這個“真實(shí)值”估算的準(zhǔn)確性及精度受采樣頻率的設(shè)置以及估算方法的選擇影響很大[1,3-5].

    當(dāng)前大多數(shù)國家的水環(huán)境監(jiān)測體系因技術(shù)或成本等原因仍面臨水文數(shù)據(jù)分布連續(xù)(實(shí)時、逐日)且相對豐富、水質(zhì)數(shù)據(jù)分布離散(雙月、逐月、雙周、逐周)且相對稀缺的現(xiàn)實(shí)[1].隨著水質(zhì)樣品采集間隔增長,污染物年負(fù)荷估算的不確定性增加[5-6].在我國,還面臨水文與水質(zhì)數(shù)據(jù)采集時間不同步、地點(diǎn)不匹配的問題[7],這進(jìn)一步增加了準(zhǔn)確估算斷面污染負(fù)荷或通量的難度,同時也凸顯了估算方法選擇的重要性.

    國際上基于實(shí)測值進(jìn)行年污染通量或負(fù)荷估算的常用方法有3種:平均法、比例法及線性回歸法[1,3,6,8].平均法采用實(shí)測斷面平均濃度、樣品時間平均濃度、斷面瞬時流量、采樣期間平均流量、采樣代表平均流量、采樣期間平均負(fù)荷等概念,構(gòu)造了多種負(fù)荷估算方法[1,6-9].比例法由Beale提出[10],它運(yùn)用一個表征污染負(fù)荷與徑流量協(xié)方差的系數(shù)來調(diào)整平均法估算的年污染負(fù)荷[8],該方法在美國、加拿大、英國的多個流域(不同采樣頻率)得到成功的應(yīng)用[1,3,8].有研究者認(rèn)為平均法和比例法的估算效果因流域特性、采樣方法、頻率的不同而異[1,11];也有研究認(rèn)為比例法能夠提供更加準(zhǔn)確的估算結(jié)果,其尤為適用于徑流量數(shù)據(jù)相對豐富、污染物濃度數(shù)據(jù)相對短缺的情況[3-4,8].線性回歸法是利用徑流量與污染物含量或負(fù)荷之間的相關(guān)關(guān)系進(jìn)行的經(jīng)驗(yàn)估算,該方法因?qū)?shù)據(jù)量的需求較小、能夠提供方差較小的無偏估計(jì)而被廣泛應(yīng)用[12],其典型代表就是LOADEST[12-13].國內(nèi)在流域出口或河流控制斷面污染負(fù)荷估算方面也提出了適宜中國現(xiàn)狀的方法,例如平均濃度法、水文估算法等,考慮了季節(jié)和徑流路徑與非點(diǎn)源污染的內(nèi)在聯(lián)系并由此分割出點(diǎn)源和非點(diǎn)源污染負(fù)荷[2,14-18].在線性回歸估算污染負(fù)荷或通量方面,洪小康等[19]發(fā)現(xiàn)在次暴雨尺度單位面積總氮、總磷、COD負(fù)荷量與徑流量存在極強(qiáng)的正相關(guān)關(guān)系,并將其應(yīng)用于年徑流污染負(fù)荷的估算;李娜等[13]則利用LOADEST模型估算了寶象河氮、磷、懸浮物的入湖通量.

    然而眾多已有研究多針對固定頻率或日期采樣模式下污染負(fù)荷的估算.本研究擬針對非常規(guī)定期采樣如流量連續(xù)監(jiān)測、水質(zhì)監(jiān)測高度集中在汛期的情形,比較不同方法在泥沙年負(fù)荷估算的表現(xiàn),并提出此種情況下泥沙年負(fù)荷的估算流程.

    1 材料與方法

    1.1 流域出口或控制斷面實(shí)測污染負(fù)荷估算思想及方法概述

    某一時段內(nèi)流域出口或控制斷面的實(shí)測污染物(包括泥沙)負(fù)荷可以通過以下公式計(jì)算:

    式中:L代表某一時段的污染負(fù)荷;Q(t)和 C(t)分別是 t時刻的瞬時流量和污染物濃度.只要得到流域出口或監(jiān)測斷面的流量、污染物濃度及其相應(yīng)的監(jiān)測時間就可以計(jì)算某個時間段內(nèi)通過該出口或監(jiān)測斷面的污染負(fù)荷.由于污染物濃度和流量隨時間變化而變化,限于觀測條件、成本等多方原因,在實(shí)際情況下不可能完全做到長時間連續(xù)同步監(jiān)測.在我國,通過水文年鑒可以收集到研究流域內(nèi)水文站點(diǎn)的逐日流量,而大多數(shù)常規(guī)水質(zhì)監(jiān)測時間跨度大,多為按照固定頻率(雙月、逐月、雙周、逐周)采集瞬時樣品,年污染負(fù)荷則是以上采樣頻率對應(yīng)時段內(nèi)污染負(fù)荷的累加.

    本文針對的對象是非固定頻率采集、樣品分布高度集中(例如汛期集中采樣、非汛期樣品稀缺)情況下的泥沙年負(fù)荷估算.按月或者季節(jié)對時段泥沙負(fù)荷進(jìn)行估算難度較大,因?yàn)槟承┰路萆踔琳麄€非汛期都缺乏泥沙含量數(shù)據(jù)[20].因此需要一種新的方法將一個水文年劃分為若干可估算負(fù)荷的時間單元,且每個時間單元應(yīng)包括一定的污染物濃度數(shù)據(jù).本文采取的是對長時間成對逐日徑流-泥沙序列進(jìn)行聚類分割為若干徑流等級,并對缺乏泥沙數(shù)據(jù)的徑流序列依此分類標(biāo)準(zhǔn)進(jìn)行分級.屬于同一徑流等級的時間(如天數(shù))合并稱之為一個徑流層,其對應(yīng)的時段泥沙負(fù)荷根據(jù)表1列舉的方法進(jìn)行估算,年負(fù)荷則是各徑流層對應(yīng)時段泥沙負(fù)荷的累加.此估算思想基于 2個重要的假設(shè):①汛期高頻率的水質(zhì)監(jiān)測覆蓋了不同的水文條件(低、中、高流量);②隸屬于同一徑流層的時間單元,其徑流量、泥沙含量相似.

    表1列舉的方法1~4是平均法的不同表現(xiàn)形式.方法1和2中時段污染負(fù)荷是流量平均與污染物濃度平均之積在時段內(nèi)的累加,區(qū)別在于方法 1以離散的瞬時流量平均作為該時段的平均流量,而方法2則采用時段內(nèi)連續(xù)流量的平均.

    方法 3是一種基于瞬時污染通量平均的算法.方法4中時段污染負(fù)荷是流量加權(quán)平均濃度與時段流量平均乘積的累加,在多篇文獻(xiàn)中被用來代表污染負(fù)荷的“真實(shí)值”[1,22],本文也用其代表污染負(fù)荷的“真實(shí)值”.方法5是比例法,是在方法4的基礎(chǔ)上運(yùn)用Beale比例對估算結(jié)果做進(jìn)一步調(diào)整.Beale比例由Beale最先提出,是一個考慮了觀測污染負(fù)荷與流量協(xié)方差的修正參數(shù),曾在北美五大湖區(qū)域廣泛用于污染負(fù)荷估算.方法 6是一種經(jīng)驗(yàn)統(tǒng)計(jì)的方法,思路是基于流量與污染物濃度或負(fù)荷之間的關(guān)系建立回歸模型,根據(jù)連續(xù)的流量數(shù)據(jù)對缺失的污染物濃度或負(fù)荷進(jìn)行預(yù)測.

    表1 基于實(shí)測數(shù)據(jù)估算年污染負(fù)荷的方法一覽[1,6,8,21]Table 1 Pollutant load estimation methods

    1.2 研究區(qū)域

    柳河子流域是灤河流域的典型支流,位于河北承德市興隆縣和承德縣.該子流域多年平均氣溫和降雨量分別是 8.4℃和 650mm.受干旱半干旱季風(fēng)氣候影響,流域降水及徑流量年際及季節(jié)性差異顯著,汛期(6~9月)雨量約占全年降水量的60%~70%.流域海拔在745~1941m之間,主要土地利用類型為林地,占流域面積的55%.本文采用灤河水文年鑒中李營水文站 1995~2011年間逐日徑流資料及有記載逐日泥沙含量,李營水文站以上控制流域面積為620km2.

    1.3 徑流聚類分層

    對柳河李營水文站17a徑流-泥沙含量數(shù)據(jù)進(jìn)行初步分析顯示,共有356條同步逐日流量-泥沙含量記錄,占總時長的 5.7%,且大部分集中在汛期(圖 1).對同步的逐日徑流-泥沙含量數(shù)據(jù)集采用K-means方法進(jìn)行聚類,分成3類的組間距離平方和解釋了整體距離平方和的 85%,表明 3個徑流層可以清晰分割.表2顯示了聚類分析得出的各徑流層數(shù)據(jù)特征,以及將徑流全時間序列進(jìn)行分層后的數(shù)據(jù)情況.結(jié)果表明17a逐日徑流-泥沙匹配數(shù)據(jù)囊括了 17a間所有高流量情形和絕大多數(shù)的中流量情形,低流量時間段采樣頻率約為 4.8%.已有眾多研究表明,在一個水文年里少數(shù)幾場降雨-徑流事件可能貢獻(xiàn)絕大部分的年污染負(fù)荷[23].因此可以預(yù)期,基于本數(shù)據(jù)集的年負(fù)荷估算結(jié)果將不會因?yàn)檫z失主要降雨-徑流事件數(shù)據(jù)而偏離真實(shí)值.

    圖1 柳河子流域1995~2011年汛期、非汛期泥沙采集頻次、平均濃度及平均流量Fig.1 Sediment sampling times, mean sediment concentrations and mean flow in the flood season and non-flood season of year 1995~2011 in Liu river sub-watershed

    表2 聚類分析及應(yīng)用結(jié)果Table 2 Results of the cluster analysis

    1.4 回歸方程建立

    對同步流量-泥沙含量數(shù)據(jù)的分析發(fā)現(xiàn),流量和泥沙含量均呈近似對數(shù)正態(tài)分布,二者之間不存在明顯的相關(guān)關(guān)系,但流量與日泥沙負(fù)荷之間存在較強(qiáng)的對數(shù)回歸關(guān)系(R2=0.709,P<0.001),且此關(guān)系不受徑流分層的強(qiáng)烈影響(圖 2).因此,本文將不再對各徑流層分別建立回歸方程以進(jìn) 行預(yù)測.

    圖2 同步流量-泥沙含量數(shù)據(jù)特征分析Fig.2 The characteristics of paired flow-sediment data

    2 結(jié)果與討論

    2.1 平均法與比例法對時段泥沙負(fù)荷估算的比較

    表3列舉了1995~2011年不同徑流層對應(yīng)時段的總天數(shù)、監(jiān)測天數(shù)、瞬時泥沙濃度平均、瞬時流量平均、時段流量平均以及基于平均法及比例法估算的泥沙年負(fù)荷.采樣頻次在年、徑流層之間有所差異,例如在 1995、1996年相對密集,而在2000、2003~2004、2006~2007、2009~2010年較為稀疏;在高、中流量時段比較密集,而在低流量時段比較稀疏.以各徑流層采樣天數(shù)占總天數(shù)的比例表征采樣頻率,各方法估算時段泥沙負(fù)荷與方法4估算結(jié)果的比值(L/Lc)表征各方法估算結(jié)果與“真實(shí)值”的偏離程度,圖3a顯示了當(dāng)采樣頻率≤30%,方法 1235估算的時段污染負(fù)荷不同程度地偏離了其對應(yīng)的“真實(shí)值”.圖3b顯示方法1和3的估算結(jié)果偏差最大,分別是方法4估算結(jié)果的0.6~12.3和1~16.1倍;方法2和5估算結(jié)果偏離真實(shí)值的程度相對較小,分別為真實(shí)值的0.51~1.23和0.95~1.12倍.而當(dāng)采樣頻率繼續(xù)增加,方法1235計(jì)算的時段泥沙負(fù)荷逐漸趨近于“真實(shí)值”(L/Lc≈1)(圖3a).對于方法1和3在采樣頻率≤30%條件下高估時段泥沙負(fù)荷的現(xiàn)象,其原因在于這兩種算法沒有考慮瞬時流量平均與時段流量平均、瞬時通量平均與時段通量平均之間的差異,也就是徑流量的時均離散項(xiàng).圖 3c進(jìn)一步揭示了方法1和3估算結(jié)果偏離真實(shí)值的程度與瞬時流量平均與時段流量平均之比值呈現(xiàn)極顯著的正相關(guān),這說明方法1和3的結(jié)果極易受到徑流量的時均離散項(xiàng)變化的影響[24].Johnes等指出方法1在徑流量比較穩(wěn)定的情況下,如具有高基流指數(shù)的流域可以獲得合理的估算值;而在徑流量變化較大的情況,如具有脈沖式徑流特征或者低基流指數(shù)的流域可產(chǎn)生極大的誤差[1].富國、李懷恩和李層也提出方法1和3適合于點(diǎn)源占優(yōu)時的污染負(fù)荷估算[7,9].圖3c還表明方法2及 5的估算結(jié)果受徑流量的時均離散差異影響較小,較為適合采樣次數(shù)稀疏或非點(diǎn)源占優(yōu)時的污染負(fù)荷估算.但同時考慮兩種方法估算結(jié)果偏離真實(shí)值的程度,方法5的結(jié)果更接近無偏,且表現(xiàn)穩(wěn)健(圖3b).

    從表1方法5的算法來看,其不僅考慮了徑流量的時均離散差異,同時還利用 Beale比例考慮了時均泥沙通量和時均流量差異的復(fù)合影響,較之方法 4僅考慮時均流量差異更進(jìn)一步對估算結(jié)果進(jìn)行了修正.然而值得注意的是,有研究認(rèn)為方法 5的估算準(zhǔn)確性會隨著采樣頻率的增加而降低[4,8].本文估算結(jié)果沒有明顯支持這一結(jié)論,但圖3d顯示當(dāng)各時段采樣頻率>30%時Lrem/Lc比例的離散程度變大,這或許表明方法5表現(xiàn)的穩(wěn)定性會因采樣頻率增加而降低.此外,方法5要求每個估算時段單元內(nèi)應(yīng)該有至少兩個流量數(shù)據(jù)以及一個污染物濃度數(shù)據(jù)以避免無法計(jì)算的情況,例如表3列舉的1998年和2001年的高流量層、以及2008年的中流量層.因這3個流量層均只包括一天,取當(dāng)日泥沙通量為方法5在該時段的泥沙負(fù)荷以用于泥沙年負(fù)荷的估算.

    圖3 時段污染負(fù)荷估算值與真實(shí)值的比值隨采樣頻率的變化(a)、在不同方法間的比較(b)、與瞬時流量平均與時段流量平均比值之間的關(guān)系(c),以及比例法估算值與真實(shí)值比值隨采樣頻率的變化(d)Fig.3 Variations in the ratio of predicted sediment loads to true loads (L/Lc) with sampling frequency increasing (a); comparisons of L/Lc ratios among four different methods(b); correlations between L/Lc ratios and ratios of average instantaneous flow velocity to average flow velocity (c); and comparisons of Lrem/Lc ratios between different sampling frequency groups(d)

    表3 徑流分層后的采樣頻率及對應(yīng)流量、泥沙含量、泥沙年負(fù)荷Table 3 Characteristics of sampling,flow,sediment concentrations and predicted sediment loads under stratified flow conditions within each year

    注:a:總天數(shù);b:監(jiān)測天數(shù);:泥沙瞬時濃度平均(mg/L);:瞬時流量平均(m3/s);μqe:時段流量平均(m3/s).

    2.2 平均法、比例法與回歸法對泥沙年負(fù)荷估算的比較

    圖 4展示了在徑流分層的基礎(chǔ)上利用方法1~6估算的泥沙年負(fù)荷.對幾乎所有年份方法 5估算的泥沙年負(fù)荷與方法 4估算的“真實(shí)值”最為接近,方法6和方法2在多數(shù)年份也產(chǎn)生了較為接近“真實(shí)值”的結(jié)果(圖4).不同方法的估算值與真實(shí)值之間的線性回歸結(jié)果表明,方法 5和 6的估算結(jié)果與“真實(shí)值”擬合效果最佳,方法 2次佳,方法1和3最差(表4).然而需要注意的是,本文方法 6估算的泥沙年負(fù)荷是逐日泥沙負(fù)荷的累積,其時間序列仍保留了原有監(jiān)測數(shù)據(jù)計(jì)算的日泥沙負(fù)荷,缺失值由回歸方程計(jì)算插值.方法 6在本文滿意的估算效果很大程度上取決于所用泥沙監(jiān)測數(shù)據(jù)捕捉到了主要降雨-徑流事件對泥沙的輸移.另有眾多研究表明,依據(jù)流量-污染物濃度回歸關(guān)系估算年污染負(fù)荷往往低估真實(shí)的污染負(fù)荷[1,8].例如Walling和Web發(fā)現(xiàn),回歸法估算泥沙年負(fù)荷可能產(chǎn)生 23%~83%的低估[25-26]. Jones等發(fā)現(xiàn),即使基于逐日水質(zhì)數(shù)據(jù),對數(shù)回歸對流域出口總磷負(fù)荷的低估程度仍可達(dá)到真實(shí)值的 45%[1].本文對比了將回歸方程運(yùn)用至全部徑流時間序列以及缺失泥沙數(shù)據(jù)的徑流時間序列而得出的泥沙年負(fù)荷(Lreg1,Lreg),發(fā)現(xiàn)前者結(jié)果僅占真實(shí)值的 18.5%左右,而后者結(jié)果卻占99.1%(表4).Horowizt指出,回歸法容易低估高流量時段的污染濃度而高估低流量時段的污染物濃度[27],這可能是造成回歸法對年污染負(fù)荷整體低估的主要原因.因此在利用回歸法估算污染負(fù)荷時,在常規(guī)采樣(逐月或雙周)的基礎(chǔ)上增加高流量時段如汛期或者主要降雨-徑流事件的水質(zhì)數(shù)據(jù)可以獲得較為滿意的估算結(jié)果[28];而在水質(zhì)數(shù)據(jù)非正態(tài)分布或者樣品量少的情況下,回歸法表現(xiàn)欠佳[12].

    表4 不同方法估算年泥沙負(fù)荷與真實(shí)值之間的線型回歸結(jié)果Table 4 Summary of the correlationships between predicted annual sediment loads and measurements

    圖4 不同方法估算的泥沙年負(fù)荷比較Fig.4 Comparion of predicted annual sediment loads by different methods

    2.3 濃度數(shù)據(jù)高度集中分布情況下泥沙年負(fù)荷估算的建議

    圖5顯示了針對非常規(guī)固定采樣模式(濃度數(shù)據(jù)高度集中在汛期)下流域出口或河流斷面泥沙的估算流程.對于既定采樣頻率的歷史數(shù)據(jù),只能依賴于估算方法的選擇來降低泥沙負(fù)荷估算的不確定性和提高估算結(jié)果的準(zhǔn)確性.一般而言,估算方法的選擇受到采樣頻率、流域大小及河流基流指數(shù)、污染物類型以及人類活動(人口、點(diǎn)源排放等)的影響[1].本研究的區(qū)域柳河子流域是典型的干旱-半干旱山區(qū)流域,其年降雨-徑流量多集中分布在6~9月的幾次降雨事件中,同時流域面積較小(620km2)對徑流緩沖能力較弱,因此徑流呈現(xiàn)明顯季節(jié)性及脈沖式特征.弱化徑流影響的算法如方法1和3顯然不適合此種情形下泥沙年負(fù)荷的估算,可供選擇的有回歸法以及考慮了徑流量時均離散差異的方法5和4.

    回歸法能否應(yīng)用于污染負(fù)荷的估算以及估算的準(zhǔn)確性,首先取決于徑流量與污染物濃度或通量之間是否存在顯著的相關(guān)關(guān)系.然而并非所有的污染物其濃度或通量與流量存在良好的相關(guān)關(guān)系.目前這種相關(guān)關(guān)系多見于徑流與泥沙、以及與泥沙緊密相關(guān)的顆粒態(tài)磷、總磷、農(nóng)藥等污染物之間,而少見于流量與溶解態(tài)污染物如硝酸鹽、氨氮、氯離子等之間[8,19,28-29].Quilbé等[8]建議,在流量-泥沙含量相關(guān)關(guān)系的決定系數(shù)>0.5的情況下,可以優(yōu)先選擇利用回歸法進(jìn)行污染負(fù)荷的估算.考慮到本研究利用了流量與泥沙通量之間存在的較強(qiáng)對數(shù)回歸關(guān)系可以得到較為滿意的估算結(jié)果,同時參考其他文獻(xiàn)類似回歸關(guān)系的應(yīng)用[29],建議若流量-泥沙通量對數(shù)回歸關(guān)系的決定系數(shù)>0.65、且成對數(shù)據(jù)中徑流量覆蓋多種水文條件(如高、中、低等),可優(yōu)先使用回歸法;而若流量-泥沙含量相關(guān)關(guān)系的決定系數(shù)≤0.5、或者流量-泥沙通量對數(shù)回歸關(guān)系的決定系數(shù)≤0.65,可考慮選擇方法5和4.

    圖5 基于實(shí)測數(shù)據(jù)估算泥沙年負(fù)荷的流程Fig.5 The flow chart for annual sediment load estimation based on measured data

    雖然方法 4考慮了徑流量的時均離散差異,而且其估算結(jié)果在本文以及其它研究中被認(rèn)為可以代表河流斷面或流域出口泥沙或其它污染物負(fù)荷的“真實(shí)值”,但與方法 5相比缺少了對污染物通量(實(shí)質(zhì)上包括了流量和濃度)的時均離散差異的考量.方法5因同時修正了流量以及污染物通量的時均差異對估算結(jié)果的影響,使其尤其適用于徑流呈季節(jié)性及脈沖式特征的流域、或者事件尺度上污染物負(fù)荷的估算;而其在較低采樣頻率下更加穩(wěn)健的估算表現(xiàn),使得它成為濃度數(shù)據(jù)缺乏情況下河流斷面或流域出口泥沙及其它污染物負(fù)荷估算的優(yōu)先選擇.因此我們建議對于時段內(nèi)采樣頻率≤30%的數(shù)據(jù)集,優(yōu)先選擇方法5進(jìn)行估算;而對于時段內(nèi)采樣頻率>30%的數(shù)據(jù)集,可以基于操作簡單的原則優(yōu)先選擇方法 4.此外,針對泥沙含量高度集中在年中某一時段的情況,按徑流量大小將一個水文年劃分為若干時段單元(各單元包含一定數(shù)量水質(zhì)數(shù)據(jù)),并對各時段單元的泥沙負(fù)荷進(jìn)行估算,可以提高泥沙年負(fù)荷估算的準(zhǔn)確性和精度[8].

    3 結(jié)語

    本文研究了水質(zhì)數(shù)據(jù)高度集中情況下泥沙年污染負(fù)荷的估算方法,提出了基于徑流分層的流量加權(quán)平均濃度法或比例法可提供合理的估算結(jié)果;其中比例法在低采樣頻率(≤30%)情況下表現(xiàn)更加穩(wěn)健.應(yīng)用回歸法進(jìn)行污染物通量或負(fù)荷的估算需要充分考慮污染物類型、流域本身流量-污染物濃度之間的關(guān)系、以及汛期污染物濃度數(shù)據(jù)的豐度.其他基于瞬時流量平均與瞬時濃度平均之積、以及瞬時通量的算法在低采樣頻率下往往高估時段泥沙負(fù)荷,主要在于沒有考慮徑流的時均離散差異,適合于泥沙含量數(shù)據(jù)比較充足情況下的負(fù)荷估算.本文的研究結(jié)果可以為干旱半干旱地區(qū)具有脈沖式徑流特征的流域的泥沙含量監(jiān)測頻率的設(shè)置及污染負(fù)荷估算方法的選擇提供一定的依據(jù).

    [1] Johnes P J. Uncertainties in annual riverine phosphorus load estimation: Impact of load estimation methodology, sampling frequency, baseflow index and catchment population density [J]. Journal of Hydrology, 2007,332(1/2):241-258.

    [2] 陸海明,彭 輝,高俊杰,等.利用水量水質(zhì)實(shí)測資料估算流域非點(diǎn)源污染負(fù)荷方法 [C]//中國水文科技新發(fā)展-中國水文學(xué)術(shù)討論會論文集. 2012,948-951.

    [3] Dolan D M, Yui A K, Geist R D. Evaluation of river load estimation methods for total phosphorus [J]. Journal of Great Lakes Research, 1981,7(3):207-214.

    [4] Richards R P, Holloway J. Monte Carlo studies of sampling strategies for estimating tributary loads [J]. Water Resources Research, 1987,23(23):1939-1948.

    [5] Williams M R, King K W, Macrae M L, et al. Uncertainty in nutrient loads from tile-drained landscapes: Effect of Sampling frequency, calculation algorithm, and compositing strategy [J]. Journal of Hydrology, 2015,530:306-316.

    [6] Webb B W, Phillips J M, Walling D E, et al. Load estimation methodologies for British Rivers and their relevance to the Lois Racs(R) Programme [J]. Science of The Total Environment, 1997,194—195:379-389.

    [7] 李懷恩,李 層,河流控制斷面污染負(fù)荷計(jì)算方法比較 [J]. 水資源與水工程學(xué)報(bào), 2013,24(2):1-4.

    [8] Quilbé R, Rousseau A N, Duchemin M, et al. Selecting a calculation method to estimate sediment and nutrient loads in streams: Application to the Beaurivage River (Québec, Canada) [J]. Journal of Hydrology, 2006,326(1—4):295-310.

    [9] 富 國.河流污染物通量估算方法分析(I)—時段通量估算方法比較分析 [J]. 環(huán)境科學(xué)研究, 2003,16(1):1-4.

    [10] Beale E. Some uses of computers in operational research [J]. Industrielle Organisation, 1962,31(1):27-28.

    [11] Mukhopadhyay B, Smith E H. Comparison of statistical methods for estimation of nutrient load to surface reservoirs for sparse data set: Application with a modified model for phosphorus availability [J]. Water Research, 2000,34(12):3258-3268.

    [12] Park Y S, Engel B A. Analysis for regression model behavior by sampling strategy for annual pollutant load estimation [J]. Journal of Environmental Quality, 2015,44(6):1843-1851.

    [13] 李 娜,盛 虎,何成杰,等.基于統(tǒng)計(jì)模型loadest的寶象河污染物通量估算 [J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 2012,20(3):355-366.

    [14] 李懷恩.估算非點(diǎn)源污染負(fù)荷的平均濃度法及其應(yīng)用 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2000,20(4):397-400.

    [15] 陳友媛,惠二青,金春姬,等.非點(diǎn)源污染負(fù)荷的水文估算方法[J]. 環(huán)境科學(xué)研究, 2003,16(1):10-13.

    [16] 李強(qiáng)坤,李懷恩,胡亞偉,等.黃河干流潼關(guān)斷面非點(diǎn)源污染負(fù)荷估算 [J]. 水科學(xué)進(jìn)展, 2008,19(4):460-466.

    [17] 李強(qiáng)坤,李懷恩,孫 娟,等.基于有限資料的水土流失區(qū)非點(diǎn)源污染負(fù)荷估算 [J]. 水土保持學(xué)報(bào), 2008,22(5):181-185.

    [18] 李家科,李懷恩,董 雯,等.渭河關(guān)中段典型支流非點(diǎn)源污染監(jiān)測與負(fù)荷估算 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2011,31(7):1470-1478.

    [19] 洪小康,李懷恩.水質(zhì)水量相關(guān)法在非點(diǎn)源污染負(fù)荷估算中的應(yīng)用 [J]. 西安理工大學(xué)學(xué)報(bào), 2000,16(4):384-386.

    [20] Thomas R B, Lewis J. An evaluation of flow-stratified sampling for estimating suspended sediment loads [J]. Journal ofHydrology, 1995,170(1):27-45.

    [21] Littlewood I, Marsh T. Annual freshwater river mass loads from great britain, 1975—1994: Estimation algorithm, database and monitoring network issues [J]. Journal of Hydrology, 2005, 304(1):221-237.

    [22] Cassidy R, Jordan P. Limitations of instantaneous water quality sampling in surface-water catchments: comparison with near-continuous phosphorus time-series data [J]. Journal of Hydrology, 2011,405(1/2):182-193.

    [23] Withers P J A, Jarvie H P. Delivery and cycling of phosphorus in rivers: A review [J]. Science of The Total Environment, 2008, 400(1—3):379-395.

    [24] 郝晨林,鄧義祥,汪永輝,等.河流污染物通量估算方法篩選及誤差分析 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2012,32(7):1670-1676.

    [25] Walling D E, Webb B W. The Reliability of Suspended Sediment Load Data, Erosion and Sediment Transport Measurement [C]// Proceedings of the Florence Symposium. IAHS Publication, 1981,133:177-194.

    [26] Walling D E, Webb B W. The Reliability of rating curve estimates of suspended sediment yield: Some further comments [C]// Proceedings of the Porto Alegre Symposium, Sediment Budgets [C]. IAHS Publication, 1988,174:337—350.

    [27] Horowitz A J. An Evaluation of sediment rating curves for estimating suspended sediment concentrations for subsequent flux calculations [J]. Hydrological processes, 2003,17(17):3387-3409.

    [28] ParkYS, Engel B A. Use of pollutant load regression models with various sampling frequencies for annual load estimation [J]. Water, 2014,6(6):1685-1697.

    [29] Su J, van Bochove E, Thériault G, et al. Effects of snowmelt on phosphorus and sediment losses from agricultural watersheds in eastern canada [J]. Agricultural water management, 2011,98(5): 867-876.

    致謝:本文感謝承德市水文水資源勘測局的數(shù)據(jù)支持.

    Comparison of different methods estimating annual sediment loads in river cross sections based on irregularly measured data.

    SU Jing-jun1, LI Xu-yong1*, WU Zhen2
    (1.State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China;2.Environmental Protection Center, Ministry of Transport, Beijing 100013, China.). China Environmental Science, 2017,37(1):218~228

    By comparing the performance of averaging method, Beale ratio method and regression method, a framework focusing on estimating annual sedimentloads based on concentrated sampling during high flow period, was proposed. The main results were:1) the Beale ratio method and flow-weighted concentration method could provide more robust and accurate estimation results regardless of the sampling frequency. The Beale ratio method performed better when samples were sparse rather than when samples were sufficient. 2) The application of regression method was conditional, heavily depending on the significance of flow-sediment correlations. Increasing storm samples in the entire calculation dataset could improve the estimation accuracy. This study could provide a useful option in designing water sampling procedures and estimating pollutant loadings in watersheds characterized by pulsed runoffs.

    load estimation;Beale ratio method;pulsed runoff

    X522

    A

    1000-6923(2017)01-0218-11

    蘇靜君(1982-),女,湖北宜昌人,助理研究員,博士.主要從事非點(diǎn)源污染研究.發(fā)表論文10余篇.

    2016-05-16

    國家自然科學(xué)基金資助項(xiàng)目(41401590);水體污染控制與治理科技重大專項(xiàng)(2015ZX07203-005-01;2012ZX07203003)

    * 責(zé)任作者, 研究員, xyli@rcees.ac.cn

    猜你喜歡
    泥沙通量徑流
    泥沙做的父親
    冬小麥田N2O通量研究
    新疆多泥沙河流水庫泥沙處理措施
    土壤團(tuán)聚體對泥沙沉降速度的影響
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    泥沙滅火
    兒童繪本(2015年2期)2015-05-25 18:10:15
    探秘“大徑流”
    攻克“大徑流”
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    好男人在线观看高清免费视频| 国产色爽女视频免费观看| 成人欧美大片| 在线观看66精品国产| 免费看a级黄色片| 深夜精品福利| 亚洲五月婷婷丁香| 久久人妻av系列| 91狼人影院| 天天一区二区日本电影三级| 老熟妇乱子伦视频在线观看| 2021天堂中文幕一二区在线观| 久久这里只有精品中国| 日韩欧美精品v在线| 婷婷色综合大香蕉| 成人性生交大片免费视频hd| xxxwww97欧美| 国产又黄又爽又无遮挡在线| 国产成+人综合+亚洲专区| 国产精品一区二区三区四区免费观看 | 久久亚洲真实| 如何舔出高潮| 又黄又爽又免费观看的视频| 国产精品综合久久久久久久免费| 国产私拍福利视频在线观看| 免费搜索国产男女视频| 成年人黄色毛片网站| 亚洲第一欧美日韩一区二区三区| 好看av亚洲va欧美ⅴa在| 欧美最新免费一区二区三区 | 一本一本综合久久| 久久草成人影院| 天堂网av新在线| 亚洲经典国产精华液单 | 成年女人毛片免费观看观看9| 国产探花在线观看一区二区| av天堂中文字幕网| 搡老岳熟女国产| 久久久久久久午夜电影| 男人和女人高潮做爰伦理| 亚洲久久久久久中文字幕| 在线播放无遮挡| 欧美日韩黄片免| 丰满的人妻完整版| www.www免费av| 欧美极品一区二区三区四区| 亚洲无线观看免费| 好男人电影高清在线观看| 97超视频在线观看视频| 给我免费播放毛片高清在线观看| 国产高清三级在线| 精品不卡国产一区二区三区| www.色视频.com| 99国产综合亚洲精品| 麻豆成人午夜福利视频| 乱码一卡2卡4卡精品| 久久精品国产亚洲av涩爱 | 日韩亚洲欧美综合| 国产在线精品亚洲第一网站| 亚洲欧美日韩东京热| 亚洲精品在线观看二区| 国产精品三级大全| 熟女人妻精品中文字幕| 午夜视频国产福利| 2021天堂中文幕一二区在线观| 听说在线观看完整版免费高清| 性色avwww在线观看| 成人亚洲精品av一区二区| 人人妻人人看人人澡| 999久久久精品免费观看国产| 色综合亚洲欧美另类图片| 免费观看的影片在线观看| 亚洲综合色惰| 69人妻影院| 老司机福利观看| 黄色配什么色好看| 香蕉av资源在线| 国产精品久久视频播放| 如何舔出高潮| 在线a可以看的网站| av在线观看视频网站免费| 国产精品一区二区三区四区久久| 欧美日韩福利视频一区二区| 色5月婷婷丁香| 91字幕亚洲| 色哟哟·www| 国产精品嫩草影院av在线观看 | 首页视频小说图片口味搜索| 亚洲国产日韩欧美精品在线观看| 最好的美女福利视频网| 午夜亚洲福利在线播放| 午夜激情福利司机影院| 91九色精品人成在线观看| 国产午夜福利久久久久久| 精品一区二区免费观看| 嫩草影院入口| 亚洲av电影不卡..在线观看| 中文字幕av在线有码专区| av在线蜜桃| АⅤ资源中文在线天堂| 国产欧美日韩一区二区精品| 欧美潮喷喷水| 少妇高潮的动态图| 2021天堂中文幕一二区在线观| 欧美一区二区亚洲| 内地一区二区视频在线| 国产精品1区2区在线观看.| 久久欧美精品欧美久久欧美| 免费高清视频大片| 美女 人体艺术 gogo| 99久久99久久久精品蜜桃| 欧美最新免费一区二区三区 | 日韩 亚洲 欧美在线| 乱码一卡2卡4卡精品| 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 中亚洲国语对白在线视频| 欧美另类亚洲清纯唯美| 日本成人三级电影网站| 免费高清视频大片| 美女 人体艺术 gogo| 国产欧美日韩精品一区二区| 日本免费一区二区三区高清不卡| 露出奶头的视频| 最新在线观看一区二区三区| 国产高潮美女av| 97人妻精品一区二区三区麻豆| 天堂动漫精品| 在线免费观看不下载黄p国产 | 国产高清视频在线播放一区| or卡值多少钱| 国产精品久久电影中文字幕| 国产精品乱码一区二三区的特点| 香蕉av资源在线| 日韩欧美精品免费久久 | 深爱激情五月婷婷| 中文字幕免费在线视频6| 亚洲国产精品合色在线| eeuss影院久久| 两个人视频免费观看高清| 国产精品久久久久久久久免 | 国产精品精品国产色婷婷| 欧美3d第一页| 欧美日韩综合久久久久久 | 一个人免费在线观看电影| 精品久久久久久成人av| ponron亚洲| 超碰av人人做人人爽久久| 精品99又大又爽又粗少妇毛片 | 欧美精品国产亚洲| 精品午夜福利视频在线观看一区| 国产高清视频在线播放一区| 夜夜爽天天搞| 赤兔流量卡办理| 深夜a级毛片| 中国美女看黄片| 99精品在免费线老司机午夜| 久久精品国产自在天天线| 日本一二三区视频观看| 久久精品国产亚洲av涩爱 | 午夜日韩欧美国产| 一进一出好大好爽视频| 国产成人av教育| 在线十欧美十亚洲十日本专区| 国产精品,欧美在线| 99在线视频只有这里精品首页| 久久久久久久久中文| 国产一区二区在线观看日韩| 亚洲成人久久爱视频| 久久精品国产亚洲av涩爱 | 欧美日韩国产亚洲二区| 国产欧美日韩一区二区三| 成人av在线播放网站| 欧美+日韩+精品| 久久久久久久久大av| 最近在线观看免费完整版| 久久久久久久久久成人| 亚洲内射少妇av| av女优亚洲男人天堂| 成人欧美大片| 欧美一级a爱片免费观看看| 精品国产亚洲在线| 99国产精品一区二区蜜桃av| 国内久久婷婷六月综合欲色啪| 成人国产综合亚洲| 国产欧美日韩一区二区三| 51国产日韩欧美| 国产av麻豆久久久久久久| 麻豆久久精品国产亚洲av| 国产成年人精品一区二区| 日韩欧美精品免费久久 | 热99在线观看视频| 最后的刺客免费高清国语| 久久精品国产亚洲av天美| .国产精品久久| 国产三级黄色录像| 欧美色视频一区免费| 午夜精品一区二区三区免费看| 亚洲中文字幕一区二区三区有码在线看| 欧美黑人欧美精品刺激| 男人的好看免费观看在线视频| 91九色精品人成在线观看| 精品久久久久久,| 日韩欧美精品免费久久 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产精品成人综合色| 色尼玛亚洲综合影院| 日韩免费av在线播放| 国产精品爽爽va在线观看网站| 中亚洲国语对白在线视频| av黄色大香蕉| 亚洲精品亚洲一区二区| 亚洲最大成人手机在线| 在线国产一区二区在线| www日本黄色视频网| 啪啪无遮挡十八禁网站| 怎么达到女性高潮| 色综合亚洲欧美另类图片| 国产精华一区二区三区| 欧美激情久久久久久爽电影| 国产精品av视频在线免费观看| 观看美女的网站| 国内精品久久久久精免费| 老熟妇仑乱视频hdxx| 亚洲内射少妇av| 免费黄网站久久成人精品 | 亚洲av免费高清在线观看| 免费观看精品视频网站| 精品人妻视频免费看| 99国产精品一区二区蜜桃av| 99热这里只有是精品50| 日本 av在线| 午夜福利视频1000在线观看| 脱女人内裤的视频| 99久久精品热视频| 欧美xxxx性猛交bbbb| 波野结衣二区三区在线| 久久久久久久久久黄片| 欧美xxxx性猛交bbbb| 9191精品国产免费久久| 精华霜和精华液先用哪个| 精品一区二区三区视频在线| 天堂av国产一区二区熟女人妻| 亚洲精品一区av在线观看| 观看美女的网站| 校园春色视频在线观看| 日本黄色视频三级网站网址| 日本免费a在线| 性欧美人与动物交配| 国产探花极品一区二区| 久久精品久久久久久噜噜老黄 | 日韩欧美免费精品| 欧美中文日本在线观看视频| 成人国产综合亚洲| 人人妻人人看人人澡| 亚洲真实伦在线观看| 亚洲av日韩精品久久久久久密| 亚洲最大成人手机在线| 日韩欧美在线乱码| 少妇被粗大猛烈的视频| 午夜精品在线福利| 天堂av国产一区二区熟女人妻| 久久欧美精品欧美久久欧美| 嫩草影院精品99| 亚洲成av人片免费观看| 人妻制服诱惑在线中文字幕| 丰满的人妻完整版| 简卡轻食公司| 欧美日韩瑟瑟在线播放| 国产欧美日韩精品亚洲av| 久久久国产成人免费| 麻豆一二三区av精品| 成人国产一区最新在线观看| 在线观看美女被高潮喷水网站 | 特级一级黄色大片| 亚洲片人在线观看| 久久九九热精品免费| 亚洲三级黄色毛片| 色在线成人网| 亚洲七黄色美女视频| 内地一区二区视频在线| 亚洲av成人不卡在线观看播放网| 久久6这里有精品| 色综合婷婷激情| 五月伊人婷婷丁香| 免费观看的影片在线观看| 麻豆一二三区av精品| 国产男靠女视频免费网站| 麻豆成人午夜福利视频| 欧美高清成人免费视频www| 深夜a级毛片| 一个人观看的视频www高清免费观看| 国产精品一区二区免费欧美| 国产三级在线视频| 欧美色欧美亚洲另类二区| 看十八女毛片水多多多| 99在线人妻在线中文字幕| 亚洲人成伊人成综合网2020| 在线观看66精品国产| 少妇裸体淫交视频免费看高清| 国产精品永久免费网站| 大型黄色视频在线免费观看| 老鸭窝网址在线观看| 老女人水多毛片| 在线观看一区二区三区| 亚洲人与动物交配视频| 99久久九九国产精品国产免费| 高清日韩中文字幕在线| 亚洲第一电影网av| 一进一出抽搐gif免费好疼| 欧美性猛交╳xxx乱大交人| 成人亚洲精品av一区二区| 久久香蕉精品热| 国产精品日韩av在线免费观看| 黄色丝袜av网址大全| 日韩人妻高清精品专区| 高清毛片免费观看视频网站| 婷婷亚洲欧美| 亚洲国产欧美人成| 成人一区二区视频在线观看| eeuss影院久久| 久久精品国产99精品国产亚洲性色| 麻豆久久精品国产亚洲av| 两个人视频免费观看高清| 啦啦啦观看免费观看视频高清| 国产私拍福利视频在线观看| 国产精品伦人一区二区| 日本成人三级电影网站| 高清毛片免费观看视频网站| 亚洲国产精品合色在线| 亚洲七黄色美女视频| 国产探花在线观看一区二区| 亚洲成人中文字幕在线播放| 麻豆国产av国片精品| 人妻制服诱惑在线中文字幕| 女人十人毛片免费观看3o分钟| 国产精品电影一区二区三区| 国产av在哪里看| 精品人妻熟女av久视频| av天堂中文字幕网| 老鸭窝网址在线观看| 免费av观看视频| 免费av毛片视频| 高清在线国产一区| 无遮挡黄片免费观看| 99热这里只有是精品在线观看 | 哪里可以看免费的av片| 免费观看精品视频网站| 免费av毛片视频| 日本在线视频免费播放| 成人三级黄色视频| 老司机午夜福利在线观看视频| 我的老师免费观看完整版| 久久6这里有精品| 国内精品一区二区在线观看| 欧美性猛交黑人性爽| 夜夜躁狠狠躁天天躁| 最新中文字幕久久久久| 欧美成人a在线观看| 婷婷亚洲欧美| 亚洲最大成人av| 9191精品国产免费久久| 九色成人免费人妻av| 脱女人内裤的视频| 在线免费观看的www视频| 女人被狂操c到高潮| 最好的美女福利视频网| 一进一出抽搐动态| 好男人电影高清在线观看| 一边摸一边抽搐一进一小说| 免费av不卡在线播放| 超碰av人人做人人爽久久| 国产亚洲精品综合一区在线观看| www.999成人在线观看| 女人被狂操c到高潮| 成人三级黄色视频| 一进一出抽搐动态| 精品人妻1区二区| 精品欧美国产一区二区三| 亚洲人成网站在线播放欧美日韩| 欧美在线一区亚洲| 99热精品在线国产| 国产不卡一卡二| 男人和女人高潮做爰伦理| 少妇的逼水好多| 色在线成人网| 一级作爱视频免费观看| 中文字幕av成人在线电影| 在线观看午夜福利视频| 欧美三级亚洲精品| 18美女黄网站色大片免费观看| 老熟妇乱子伦视频在线观看| 国产精品女同一区二区软件 | 97超视频在线观看视频| 亚洲美女搞黄在线观看 | 亚洲国产欧洲综合997久久,| 欧美日韩亚洲国产一区二区在线观看| 国产免费av片在线观看野外av| 亚洲天堂国产精品一区在线| 国产av不卡久久| 成人精品一区二区免费| a级毛片免费高清观看在线播放| 九色国产91popny在线| 99热只有精品国产| 99热这里只有精品一区| 自拍偷自拍亚洲精品老妇| 一区二区三区四区激情视频 | 久久久成人免费电影| 精品乱码久久久久久99久播| 免费黄网站久久成人精品 | 亚洲国产欧美人成| 99精品在免费线老司机午夜| 久久精品国产亚洲av香蕉五月| 又爽又黄a免费视频| 午夜福利欧美成人| 欧美黄色淫秽网站| 欧美在线黄色| 麻豆av噜噜一区二区三区| ponron亚洲| 日本一二三区视频观看| 久久香蕉精品热| 好男人在线观看高清免费视频| 小说图片视频综合网站| 69av精品久久久久久| 精品不卡国产一区二区三区| 中亚洲国语对白在线视频| 人妻制服诱惑在线中文字幕| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品合色在线| 午夜亚洲福利在线播放| 国语自产精品视频在线第100页| 偷拍熟女少妇极品色| 久久久久久久久中文| 午夜亚洲福利在线播放| 亚洲,欧美,日韩| 男女床上黄色一级片免费看| 欧美高清成人免费视频www| 亚洲精品粉嫩美女一区| 啦啦啦韩国在线观看视频| 精品久久久久久久久av| 99久久精品一区二区三区| www.999成人在线观看| 麻豆一二三区av精品| 变态另类丝袜制服| 中文字幕久久专区| 最近视频中文字幕2019在线8| 丁香六月欧美| 少妇的逼好多水| 国产熟女xx| a级毛片免费高清观看在线播放| 精品国内亚洲2022精品成人| 最近视频中文字幕2019在线8| 国产欧美日韩精品一区二区| 国产av不卡久久| 91狼人影院| 国产一区二区在线av高清观看| 亚洲五月婷婷丁香| 757午夜福利合集在线观看| 日韩欧美精品免费久久 | 全区人妻精品视频| 亚洲欧美日韩高清专用| 久久欧美精品欧美久久欧美| 中文字幕免费在线视频6| 丰满人妻一区二区三区视频av| 精品国内亚洲2022精品成人| 免费观看人在逋| 亚洲av五月六月丁香网| 亚洲最大成人手机在线| 精品国产亚洲在线| 99热这里只有是精品在线观看 | 精品无人区乱码1区二区| 精华霜和精华液先用哪个| 18禁黄网站禁片午夜丰满| 麻豆久久精品国产亚洲av| 十八禁人妻一区二区| 在线观看av片永久免费下载| 精品人妻熟女av久视频| 国内少妇人妻偷人精品xxx网站| 欧美中文日本在线观看视频| 99国产极品粉嫩在线观看| 国产精品美女特级片免费视频播放器| 一级a爱片免费观看的视频| 日本撒尿小便嘘嘘汇集6| 成人av一区二区三区在线看| 欧美黄色淫秽网站| 无遮挡黄片免费观看| 精品日产1卡2卡| 97人妻精品一区二区三区麻豆| 国产中年淑女户外野战色| 亚洲 国产 在线| 成熟少妇高潮喷水视频| 国产高清激情床上av| 成人亚洲精品av一区二区| 国产精品一区二区三区四区久久| 精品一区二区三区av网在线观看| 国产精品一及| 9191精品国产免费久久| 亚洲av一区综合| 午夜两性在线视频| 欧美xxxx黑人xx丫x性爽| 亚洲 国产 在线| 亚洲av电影不卡..在线观看| 男人舔女人下体高潮全视频| 午夜精品一区二区三区免费看| 免费看美女性在线毛片视频| 天天一区二区日本电影三级| 一本一本综合久久| 日本黄色视频三级网站网址| 美女高潮的动态| 亚洲成人精品中文字幕电影| 乱码一卡2卡4卡精品| 一进一出好大好爽视频| 免费av观看视频| 亚洲最大成人av| 亚洲专区国产一区二区| 一边摸一边抽搐一进一小说| 人人妻,人人澡人人爽秒播| 精品熟女少妇八av免费久了| 欧美bdsm另类| 欧洲精品卡2卡3卡4卡5卡区| 欧美一区二区国产精品久久精品| 中文字幕高清在线视频| 日本 欧美在线| 全区人妻精品视频| 麻豆av噜噜一区二区三区| 国产一区二区三区视频了| 久9热在线精品视频| 村上凉子中文字幕在线| 成熟少妇高潮喷水视频| 久久精品影院6| 人人妻人人看人人澡| 一区二区三区免费毛片| 99国产精品一区二区三区| 久久久色成人| 18禁黄网站禁片免费观看直播| 亚洲第一欧美日韩一区二区三区| 十八禁网站免费在线| 国产av麻豆久久久久久久| 免费搜索国产男女视频| 九九在线视频观看精品| xxxwww97欧美| 精品久久久久久久久亚洲 | 人人妻人人看人人澡| 久久久精品欧美日韩精品| 国产精品亚洲一级av第二区| 美女cb高潮喷水在线观看| 亚洲欧美日韩无卡精品| 色综合欧美亚洲国产小说| 亚洲熟妇熟女久久| 久久久久久久久中文| 欧美成人一区二区免费高清观看| 波野结衣二区三区在线| 国产视频内射| 男人舔奶头视频| 久久午夜福利片| 国产一级毛片七仙女欲春2| 美女cb高潮喷水在线观看| 日韩欧美精品v在线| 色综合欧美亚洲国产小说| 国产成人福利小说| 中文字幕高清在线视频| 欧美性猛交╳xxx乱大交人| 久久精品91蜜桃| 精品人妻熟女av久视频| 日韩欧美一区二区三区在线观看| 久久伊人香网站| 欧美日本亚洲视频在线播放| 99热这里只有是精品50| 欧美性猛交黑人性爽| 欧美日韩综合久久久久久 | 超碰av人人做人人爽久久| 午夜福利在线观看免费完整高清在 | 女同久久另类99精品国产91| 午夜免费激情av| 日韩免费av在线播放| 尤物成人国产欧美一区二区三区| 香蕉av资源在线| 欧美丝袜亚洲另类 | 99国产综合亚洲精品| 熟女电影av网| 国产蜜桃级精品一区二区三区| 在线观看舔阴道视频| 国产一区二区三区在线臀色熟女| 精品国产亚洲在线| 亚洲成人久久爱视频| 亚洲综合色惰| 欧美一级a爱片免费观看看| 亚洲av熟女| 欧美激情久久久久久爽电影| 非洲黑人性xxxx精品又粗又长| 日本五十路高清| 成熟少妇高潮喷水视频| 午夜福利免费观看在线| 熟妇人妻久久中文字幕3abv| 熟女人妻精品中文字幕| 亚洲乱码一区二区免费版| 啦啦啦观看免费观看视频高清| 国产精品精品国产色婷婷| 免费在线观看影片大全网站| 亚洲成人久久爱视频| 国产精品一及| 欧美一级a爱片免费观看看| 性插视频无遮挡在线免费观看| 日本黄色视频三级网站网址| 国产v大片淫在线免费观看| 亚洲精品色激情综合| 亚洲午夜理论影院| 一区二区三区四区激情视频 | 免费在线观看影片大全网站| 18禁黄网站禁片午夜丰满| 在线免费观看的www视频| 色综合婷婷激情| 好男人在线观看高清免费视频| 国产野战对白在线观看| 亚洲精品456在线播放app |