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

    大汶河水質(zhì)狀況評價和污染源分析

    2021-10-08 08:57:42申洪鑫牛蓓蓓李富強李新舉
    安全與環(huán)境工程 2021年5期
    關(guān)鍵詞:大汶河水體流域

    申洪鑫,牛蓓蓓,李富強,王 英,李新舉*

    (1.山東農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,山東 泰安 271018;2.山東省煤田地質(zhì)局第三勘探隊,山東 泰安 250102)

    水資源是人類賴以生存和社會發(fā)展不可或缺的物質(zhì)資源。20世紀70年代以來,隨著人口的快速增加和經(jīng)濟的高速發(fā)展,水資源的需求量短時間內(nèi)急劇增長,同時城鄉(xiāng)污水大量排放污染了地表水及地下水,不可避免地形成了“水質(zhì)性缺水”等一系列水污染問題。我國70%江河湖泊的水質(zhì)受到各種污染物的影響,水污染不僅破壞了自然生態(tài)景觀,而且會影響經(jīng)濟發(fā)展和人類健康。因此,適時開展河流水質(zhì)評價和污染源解析,有助于對河流進行有效管理和污染防治,對于區(qū)域社會-經(jīng)濟-生態(tài)的可持續(xù)發(fā)展具有重要的意義。

    大汶河為黃河下游在山東境內(nèi)最大的支流,亦是山東省泰安市地表水資源的主要來源。它流經(jīng)東平湖注入黃河,作為南水北調(diào)東線工程重要的蓄水走廊——東平湖的唯一匯入河流,大汶河水質(zhì)狀況勢必會影響南水北調(diào)工程的調(diào)水水質(zhì)。隨著社會經(jīng)濟發(fā)展和城市化進程的加快,大汶河流域內(nèi)水資源短缺和水污染問題日漸凸顯。近年來,左欣等利用模糊數(shù)學(xué)法對2000年大汶河主要控制斷面水質(zhì)進行了評價,發(fā)現(xiàn)大汶河水質(zhì)污染嚴重,大部分河段水質(zhì)都不能夠滿足功能區(qū)的需求;何濤等通過分析大汶河沿線4個監(jiān)測點的水質(zhì)資料發(fā)現(xiàn),2007—2014年間受流域各項減排措施的影響,大汶河有機污染得到有效緩解,受攔蓄工程影響,河段內(nèi)TP含量明顯下降,TN含量總體上升,其水質(zhì)狀況整體有所好轉(zhuǎn);申祺等采用主成分分析和相關(guān)分析方法對2017年大汶河的采樣分析數(shù)據(jù)進行了水生態(tài)系統(tǒng)健康狀況評價,結(jié)果表明大汶河水生態(tài)系統(tǒng)健康狀況評價結(jié)果主要以一般和較差為主,分別占總采樣點的58.33%和20.83%。目前的已有研究主要偏重于常規(guī)水質(zhì)指標的分析評價,對于大汶河重金屬污染狀況和污染來源的分析研究較少。為此,本文基于2019年6月份和8月份大汶河21個斷面的兩次采樣分析數(shù)據(jù),采用單因子評價法、內(nèi)梅羅指數(shù)法和因子分析法等方法對大汶河水質(zhì)狀況進行了評價和污染源解析,以期為大汶河水體污染防治提供依據(jù)。

    1 研究區(qū)概況

    大汶河位于山東省中部泰山南麓,介于116°E~118°E、35.7°N~36.6°N之間,發(fā)源于山東沂源縣境內(nèi),自東向西流經(jīng)萊蕪、新泰、泰安、肥城、東平等縣、市,匯注東平湖,出陳山口后入黃河。大汶河全長208 km,流域面積為9 098 km,以大汶口、戴村壩為節(jié)點劃分為上、中、下游,東平縣馬口以下稱東平湖區(qū)。大汶河上游為主要集水區(qū),地形以低山丘陵為主,河岸穩(wěn)定,土壤肥沃;大汶河中、下游地形以平原為主,且土壤以沙土和沙壤土為主,透水通氣,保溫保濕性能良好,適合農(nóng)業(yè)發(fā)展。大汶河多年平均徑流量約為18.2億m,灌溉面積達290余萬畝。大汶河為沿岸居民的生產(chǎn)和生活提供用水保障,是山東省內(nèi)重要的河流之一。

    2 研究方法

    2.1 采樣點布設(shè)

    根據(jù)大汶河河流特征和污染源分布特征,在大汶河主要干流和支流布設(shè)了21個采樣斷面,編號依次為1~21號,采樣點分布如圖1所示。其中,1~9號采樣斷面位于大汶河北支的灜汶河、牟汶河;11、12號采樣斷面位于大汶河南支(柴汶河);大汶口以下大汶河干流采樣斷面5個(編號為10、16、17、19、20);海子河采樣斷面1個(編號為13);漕河采樣斷面2個(編號為14、15);匯河采樣斷面1個(編號為18);東平湖采樣斷面1個(編號為21)。

    圖1 大汶河監(jiān)測斷面分布圖Fig.1 Distribution map of monitoring sections in the Dawen River

    2.2 水樣采集與分析方法

    大汶河水質(zhì)檢測儀器見表1,具體檢測方法依據(jù)《地下水質(zhì)檢驗方法》(DZ/T 0064—1993)、《生活飲用水標準檢驗方法》(GB/T 5750—2006),并遵循《地表水環(huán)境質(zhì)量標準》(GB 3828—2002)的要求。

    表1 大汶河水質(zhì)檢測儀器Table 1 Water quality testing instruments in this study

    2.3 水質(zhì)評價方法

    2.3.1 單因子污染指數(shù)法

    單因子污染指數(shù)法是通過指數(shù)運算出在所有進行水質(zhì)評價的指標中最差的單項水質(zhì)指標,隨后通過此指標所屬的類別來確定水體的綜合水質(zhì)類別,實質(zhì)上是用水體各評價因子的監(jiān)測結(jié)果對照該項目的分類標準,確定其水質(zhì)類別。單因子污染指數(shù)法是進行水體水質(zhì)評價過程中應(yīng)用最為廣泛的評價方法之一,它有效地避免了確定指標權(quán)重的主觀隨意性。在此基礎(chǔ)上,徐祖信將其改進為單因子水質(zhì)標識指數(shù)法,并且已經(jīng)在水體水質(zhì)評價工作中得到了一定的應(yīng)用。

    單因子污染指數(shù)法的計算公式為

    (1)

    式中:

    P

    為采樣點污染物

    i

    的污染指數(shù);

    C

    為采樣點污染物

    i

    濃度的實測值(mg/L);

    S

    為采樣點污染物

    i

    的標準濃度值(mg/L)(取《地表水環(huán)境質(zhì)量標準》中的Ⅲ類水質(zhì)標準)。

    單因子污染指數(shù)法環(huán)境質(zhì)量評價分級標準,見表2。

    表2 單因子污染指數(shù)法環(huán)境質(zhì)量評價分級標準Table 2 Classification criteria for environmental quality assessment of single-factor pollution index method

    2.3.2 內(nèi)梅羅污染指數(shù)法

    內(nèi)梅羅污染指數(shù)最早是由美國敘拉古大學(xué)內(nèi)梅羅(N.L.Nemerow)教授在其所著的《河流污染科學(xué)分析》一書中提出的一種水污染指數(shù)。與單因子污染指數(shù)法相比,內(nèi)梅羅污染指數(shù)法在突出污染最嚴重的評價因子的同時,更兼顧了評價體系中其他因子的貢獻度,是一種比較綜合的評價方法。

    內(nèi)梅羅污染指數(shù)法的計算公式為

    (2)

    內(nèi)梅羅污染指數(shù)法環(huán)境質(zhì)量評價分級標準,見表3。

    表3 內(nèi)梅羅污染指數(shù)法環(huán)境質(zhì)量評價分級標準Table 3 Classification criteria for environmental quality assessment of Nemerow pollution index method

    2.4 水質(zhì)時空變化評價方法

    水質(zhì)隨時間和空間的變化程度分別根據(jù)其時間變化率

    T

    和空間變化率

    S

    來判斷,其計算公式如下:

    (3)

    (4)

    上式中:

    P

    P

    分別為比較時間內(nèi)起始時刻和終止時刻采樣斷面的內(nèi)梅羅污染指數(shù);

    P

    P

    分別為比較空間內(nèi)起始斷面和終止斷面采樣點的內(nèi)梅羅污染指數(shù)。

    2.5 因子分析法

    因子分析法是一種變量降維的方法,它假設(shè)原有變量的背后存在多個不可觀測的隱藏因子,它從研究原始變量的相關(guān)矩陣出發(fā),通過降維的方法將眾多原始變量歸納為少數(shù)幾個綜合因子,而這些因子能揭示出隱藏在原始變量之間的聯(lián)系,尤其是成因上的聯(lián)系,繼而解釋基于原始變量得到的旋轉(zhuǎn)后的成分矩陣、成分得分系數(shù)矩陣。

    在進行因子分析之前,需要對原始數(shù)據(jù)進行Z標準化處理,以消除不同污染物濃度量綱差異所帶來的影響。

    2.6 數(shù)據(jù)處理

    數(shù)據(jù)處理采用Excel、IBM SPSS Statistics 23軟件,作圖采用ArcGIS 10.2軟件。

    3 結(jié)果與分析

    3.1 大汶河水質(zhì)統(tǒng)計特征

    大汶河2期水質(zhì)指標統(tǒng)計結(jié)果,見表4。其中,各個水質(zhì)指標的平均值用于表征水體各水質(zhì)指標整體的含量水平;變異系數(shù)主要用來反映水體化學(xué)成分的變異程度,揭示水體水質(zhì)化學(xué)指標在空間上的差異。

    由表4可以看出:

    表4 大汶河豐、枯水期水質(zhì)指標統(tǒng)計結(jié)果Table 4 Statistical results of water quality indexes of the Dawen River in high flow and low flow period

    3.2 大汶河水質(zhì)評價

    3.2.1 單因子污染指數(shù)法水質(zhì)評價

    根據(jù)公式(1),對大汶河流域21個監(jiān)測斷面的16個監(jiān)測指標進行計算,得到大汶河流域各監(jiān)測斷面單因子污染指數(shù)見表5,其污染等級分類結(jié)果見圖2。

    表5 大汶河流域各監(jiān)測斷面水體水質(zhì)污染評價結(jié)果(單因子污染指數(shù)法)Table 5 Evaluation results of water quality pollution of monitoring sections in the Dawen River Basin(Single factor pollution index method)

    圖2 大汶河流域各監(jiān)測斷面水體水質(zhì)污染評價結(jié)果(單因子污染指數(shù)法)Fig.2 Evaluation results of water quality pollution of monitoring sections in the Dawen River Basin (single factor pollution index method)

    大汶河流域處水質(zhì)主要污染物為TN、Se和Hg。其中,只有在豐水期水體中發(fā)生Se含量超標情況,而含Se廢水主要是由顏料染料、硫酸制造和精煉銅等行業(yè)排放,經(jīng)調(diào)查發(fā)現(xiàn)大汶河沿岸有紡織業(yè)工廠的存在。Se中毒是一個緩慢的過程,對生態(tài)環(huán)境、人體健康都有極大的潛在危害,故有關(guān)部門需引起重視。大汶河上游的劉家廟村、前王莊,中游的石家門村和下游的匯河入大清河斷面處水體水質(zhì)污染嚴重,主要污染物為TN和Se,大部分斷面處水體中TN超標,存在較大的富營養(yǎng)化風(fēng)險。

    3.2.2 內(nèi)梅羅污染指數(shù)法水質(zhì)評價

    根據(jù)公式(2),對大汶河流域21個監(jiān)測斷面的16個監(jiān)測指標進行計算,得到基于內(nèi)梅羅污染指數(shù)法的大汶河流域各監(jiān)測斷面水體水質(zhì)污染評價結(jié)果,見表6。

    表6 大汶河流域各監(jiān)測斷面水體水質(zhì)污染評價結(jié)果(內(nèi)梅羅污染指數(shù)法)Table 6 Evaluation results of water quality pollution of monitoring sections in the Dawen River Basin(Nemerow pollution index method)

    結(jié)合表5和表6的數(shù)據(jù)可以看出,由于Se和Hg含量超標,大汶河流域大部分斷面水體中最大污染因子為Se和Hg,致使大汶河流域不同區(qū)域水體污染狀況以重度污染為主。

    大汶河流域水體污染程度由高到低可劃分為:大汶河上游>大汶河中游>大汶河下游。其中,大汶河上游區(qū)域,由于趙莊橋斷面、北店子橋斷面和北望泉林壩上游斷面位于中心城區(qū)河流下游,居民生產(chǎn)的生活用水量大,污水排放較多,且附近工廠較多,工業(yè)廢水集中排放,最終導(dǎo)致這3個斷面處水體水質(zhì)指標嚴重超標;大汶河中游區(qū)域,寧陽海子河斷面處水體水質(zhì)指標超標原因與大汶河上游區(qū)域相同,漕河肖家店橋斷面處水體水質(zhì)在此次監(jiān)測過程中超標最為嚴重,主要超標項為Hg。Hg在自然界中以三種形式存在,其中有機Hg對生態(tài)環(huán)境的危害最大,其具有穩(wěn)定性強、毒性高、生物積累性強等特點。水體中Hg的來源有很多,自然來源包括地表徑流、大氣降水等,其中地表徑流中的Hg大部分來源于可溶解土壤中的Hg;人為來源主要包括化工、冶金、輕工等多種行業(yè)排放的工業(yè)廢水。經(jīng)調(diào)查,漕河上游有化工、冶金等工廠的存在,大量工業(yè)廢水的排放導(dǎo)致水體中Hg的含量嚴重超標。

    采用單因子污染指數(shù)法和內(nèi)梅羅污染指數(shù)法對大汶河流域各監(jiān)測斷面水體水質(zhì)污染狀況進行評價的結(jié)果,見圖3。

    圖3 大汶河流域各監(jiān)測斷面水體水質(zhì)污染評價結(jié)果Fig.3 Water pollution assessment results of monitoring sections in the Dawen River Basin

    由圖3可見,兩種方法的評價結(jié)果基本一致,但內(nèi)梅羅污染指數(shù)法的評價更全面、更可靠。根據(jù)內(nèi)梅羅污染指數(shù)法水質(zhì)評價結(jié)果,流域內(nèi)枯水期Ⅴ類水質(zhì)斷面為14個,占比為66.7%;豐水期Ⅴ類水質(zhì)斷面為17個,占比為81%。整體來說,大汶河流域枯水期水體水質(zhì)明顯優(yōu)于豐水期,且從河流上游至下游污染程度逐漸加重。

    3.3 大汶河水質(zhì)時空變化分析

    根據(jù)公式(3)和(4),計算大汶河流域各監(jiān)測斷面水體水質(zhì)的時間變化率

    T

    和空間變化率

    S

    ,其結(jié)果見表7。

    表7 大汶河流域各監(jiān)測斷面水體水質(zhì)的時空變化Table 7 Spatiotemporal changes of water quality at monitoring sections in the Dawen River Basin

    由表7可知:大汶河豐水期與枯水期相比,監(jiān)測斷面2、8、10、16、18、20處的水質(zhì)明顯惡化(時間變化率平均值為-32.59%),監(jiān)測斷面1、5、6、17、19、21處的水質(zhì)輕微惡化(時間變化率平均值為-7.09%),監(jiān)測斷面3、4、7、9、11、12、13、14、15處的水質(zhì)改善(時間變化率平均值為18.44%)。其中,監(jiān)測斷面6、8、16、17均位于水壩上游,枯水期無補充水源,夏季用水高峰水庫水被調(diào)用,造成水量持續(xù)減少,水質(zhì)污染幾率降低。受2019年第9號臺風(fēng)“利奇馬”的影響,全省降水量較2018年同期偏多50%,雨水帶來的面源性污染是造成2019年大汶河豐水期水質(zhì)惡化的主要原因之一。

    本文選擇大汶河具有代表性的9處干流監(jiān)測斷面分析其水質(zhì)空間變化,并根據(jù)支流匯入的上、下游兩個斷面水質(zhì)的變化情況對支流水質(zhì)對干流水質(zhì)的影響進行了分析。

    3.4 因子分析

    為了消除不同變量單位差異而造成因子分析的誤差,在SPSS 23軟件中先對各采樣點的原始數(shù)據(jù)進行Z標準化處理,再進行因子分析。在進行因子分析前,須進行KMO檢驗和Bartlett球形度檢驗,以判斷數(shù)據(jù)是否適合進行因子分析。本研究中,KMO統(tǒng)計量顯示為0.57,一般認為KMO值大于0.5即說明適合做因子分析,因此本文數(shù)據(jù)適合做因子分析;Bartlett球形度檢驗的近似卡方值為259.724,自由度為105,概率值為0,即使球形假設(shè)被拒絕,即可以認為相關(guān)系數(shù)矩陣與單位矩陣有顯著差異,說明本文數(shù)據(jù)適合進行因子分析。

    運用主成分分析方法提取因子,采用最大方差法進行正交旋轉(zhuǎn),按照累計貢獻率80%以上和初始特征值大于1的原則抽取因子。由于本研究參與分析的樣本數(shù)量多達42個,因子數(shù)量為15個,樣本數(shù)與因子數(shù)都較大,故提取代表性的公因子的難度也較大,因此本次分析結(jié)果76.48%接近80%,可視為可以接受的范圍。本研究獲得6個公因子,旋轉(zhuǎn)前和旋轉(zhuǎn)后的累計方差貢獻率沒有發(fā)生變化,說明這6個公因子可以用來解釋大汶河流域15個原始變量的基本信息。旋轉(zhuǎn)后的因子載荷矩陣和各公因子的水質(zhì)主要控制指標,見表8和表9。

    表8 旋轉(zhuǎn)后的因子載荷矩陣Table 8 Factor loading matrix after rotation

    表9 2019年枯水期與豐水期提取的公因子Table 9 Common factors extracted from both low and high flow periods in 2019

    4 討 論

    4.1 大汶河水體污染來源分析

    大汶河流域水體污染主要分為兩個方面:一方面為水體富營養(yǎng)化(以TN污染為主);另一方面為水體重金屬污染(以Hg和Se污染為主)。

    大汶河沿岸土地利用類型以耕地為主,河流中TN與耕地呈正相關(guān)。這是由于耕地氮肥施加過多加重了水體中TN污染的程度,從而對流域水體形成面源性污染。大汶河水體中重金屬的主要來源是工業(yè)廢水排放,主要來自冶金、化工、紡織等行業(yè)。這是由于近年來該流域內(nèi)工業(yè)化進程的加快使得工業(yè)廢水排放量增加,從而對大汶河流域水體形成重金屬污染。迄今大汶河流域泰安市域內(nèi)仍有省控以上重點污水排放企業(yè)82家,涉及紡織、食品、造紙、化工、冶金、礦業(yè)等行業(yè)。此外,化工、冶煉及燃煤動力等也會產(chǎn)生含有Hg、Se等重金屬的廢氣,通過沉降在土壤和水體中富集。李瑞平等的研究證明泰安市6個縣市區(qū)農(nóng)田土壤Hg含量平均值是山東省土壤元素背景值的17.8倍,土壤中累積的重金屬會隨雨水進入水體,進一步對河流產(chǎn)生污染。

    4.2 大汶河流域管理對策

    水污染已嚴重威脅著我國的水資源安全,水污染的防治不僅僅是環(huán)境問題,更是關(guān)系到國家安全的戰(zhàn)略性問題,水污染的治理以及河道生態(tài)修復(fù)是我國可持續(xù)性發(fā)展的必要課題。

    大汶河流域上游以農(nóng)業(yè)面源污染為主,對于農(nóng)業(yè)面源污染的治理,可參考《重點流域農(nóng)業(yè)面源污染綜合治理示范工程建設(shè)規(guī)劃(2016—2020年)》,根據(jù)大汶河流域農(nóng)業(yè)面源污染組成的特征,因地制宜地建設(shè)農(nóng)田面源污染綜合防控工程,加強農(nóng)藥化肥的減量增效、農(nóng)業(yè)廢棄物的循環(huán)利用以及畜禽、水產(chǎn)養(yǎng)殖的污染治理;同時,要堅持生產(chǎn)生態(tài)的協(xié)調(diào)發(fā)展,堅持統(tǒng)籌推進,以地方政府為主體,以縣為基本單元,開展農(nóng)業(yè)面源污染的綜合治理。

    大汶河流域中、下游以水體重金屬污染為主,且在豐水期污染更為嚴重,表明存在某些企業(yè)在雨季增加不合格污水排放的情況。因此,管理部門應(yīng)加大監(jiān)管力度,提高監(jiān)管信息化水平,設(shè)計多個關(guān)鍵性指標,通過對這些指標的監(jiān)測來判斷工廠的污水處理是否合格,并責令不達標的企業(yè)停業(yè)整頓,同時將污染重的企業(yè)遷離重要水源地;同時,應(yīng)制定激勵政策,鼓勵企業(yè)開展污水治理技術(shù)創(chuàng)新,既保證了生態(tài)效益又保證了經(jīng)濟效益。此外,工業(yè)企業(yè)應(yīng)響應(yīng)國家號召,按照污染物總量控制的原則,擴建污水處理設(shè)施,增強污水處理的能力,提高水資源的利用率;與政府共同協(xié)商,大力推動清潔生產(chǎn),優(yōu)化工業(yè)園區(qū)廢水排放管理。

    除此之外,河長制是現(xiàn)階段水環(huán)境治理中的一種新型的行政管理體制,即由各級黨政主要負責人擔任“河長”,各自負責轄區(qū)內(nèi)河流的污染治理職責,并明確其責任、剛性化考核,在水環(huán)境治理中發(fā)揮其“領(lǐng)頭羊”的作用,同時還要加強監(jiān)督,構(gòu)建科學(xué)、有效的監(jiān)督管理體系。

    水體污染的治理需要多方力量的共同參與,政府有關(guān)部門要扮演好“服務(wù)者”的角色,企業(yè)是主體,個人是“監(jiān)督者”,要做到“共同治理,全民參與”。

    5 結(jié) 論

    大汶河流域水體整體水質(zhì)以Ⅴ類水為主,流域內(nèi)枯水期Ⅴ類水質(zhì)斷面占比為66.7%,豐水期Ⅴ類水質(zhì)斷面占比為81%。流域內(nèi)水體中主要污染物以TN、Hg和Se等為主。

    從時間變化來看,大汶河枯水期水質(zhì)明顯優(yōu)于豐水期;從空間變化來看,大汶河上游至下游水質(zhì)污染程度逐漸加重。

    大汶河流域水體整體以面源性污染為主,水體污染主要分為兩個方面:一方面為水體富營養(yǎng)化(以TN污染為主),主要與農(nóng)藥、化肥的使用有關(guān);另一方面為水體重金屬污染(以Hg和Se污染為主),主要與工業(yè)廢水排放有關(guān)。

    猜你喜歡
    大汶河水體流域
    由大汶河到泰山
    壓油溝小流域
    農(nóng)村黑臭水體治理和污水處理淺探
    生態(tài)修復(fù)理念在河道水體治理中的應(yīng)用
    堡子溝流域綜合治理
    羅堰小流域
    大汶河流域水資源質(zhì)量狀況及變化趨勢分析
    中國水運(2020年11期)2020-12-23 04:22:54
    大汶河梯級攔蓄水庫富營養(yǎng)化評估
    打造智慧流域的思路及構(gòu)想——以討賴河流域為例
    廣元:治理黑臭水體 再還水清岸美
    欧美日韩中文字幕国产精品一区二区三区 | 好男人电影高清在线观看| 在线观看免费午夜福利视频| 亚洲黑人精品在线| 高潮久久久久久久久久久不卡| 伊人久久大香线蕉亚洲五| 欧美久久黑人一区二区| 欧美乱妇无乱码| 国产av精品麻豆| www.999成人在线观看| 亚洲性夜色夜夜综合| 日韩高清综合在线| 午夜福利成人在线免费观看| 一本大道久久a久久精品| 午夜两性在线视频| 99国产极品粉嫩在线观看| 麻豆久久精品国产亚洲av| 高清在线国产一区| 欧美乱妇无乱码| 中文字幕人妻丝袜一区二区| 精品人妻1区二区| 色在线成人网| 女人高潮潮喷娇喘18禁视频| 亚洲av片天天在线观看| 99精品久久久久人妻精品| 国产xxxxx性猛交| 美女国产高潮福利片在线看| 欧美精品亚洲一区二区| 99久久综合精品五月天人人| 午夜亚洲福利在线播放| 少妇粗大呻吟视频| 99久久99久久久精品蜜桃| 国产精品电影一区二区三区| 黄色成人免费大全| 日本 av在线| 欧美成人性av电影在线观看| 欧美丝袜亚洲另类 | 两人在一起打扑克的视频| 久久久水蜜桃国产精品网| 午夜视频精品福利| 黄片小视频在线播放| 亚洲视频免费观看视频| 亚洲专区国产一区二区| 亚洲国产精品久久男人天堂| 午夜精品在线福利| 黑丝袜美女国产一区| 亚洲熟妇熟女久久| 满18在线观看网站| 亚洲avbb在线观看| 韩国av一区二区三区四区| 亚洲中文字幕一区二区三区有码在线看 | 亚洲美女黄片视频| 久久国产亚洲av麻豆专区| 亚洲欧美激情综合另类| 欧美人与性动交α欧美精品济南到| 非洲黑人性xxxx精品又粗又长| 国产熟女午夜一区二区三区| 亚洲激情在线av| 国产亚洲av高清不卡| 国产高清videossex| 精品无人区乱码1区二区| 午夜福利成人在线免费观看| 亚洲av成人一区二区三| 久久精品91蜜桃| 一区福利在线观看| 国产精品久久久人人做人人爽| 国产亚洲精品久久久久久毛片| 日韩欧美在线二视频| 亚洲av第一区精品v没综合| 他把我摸到了高潮在线观看| 天堂动漫精品| 午夜免费激情av| 中文字幕色久视频| 国产片内射在线| 大码成人一级视频| 亚洲精品中文字幕一二三四区| 美女午夜性视频免费| www.自偷自拍.com| 老熟妇仑乱视频hdxx| 纯流量卡能插随身wifi吗| 午夜福利高清视频| 欧美 亚洲 国产 日韩一| 国产亚洲精品综合一区在线观看 | 日本一区二区免费在线视频| 久久精品人人爽人人爽视色| 婷婷精品国产亚洲av在线| 久久久久久免费高清国产稀缺| 老司机靠b影院| 熟妇人妻久久中文字幕3abv| 色尼玛亚洲综合影院| 女警被强在线播放| 中亚洲国语对白在线视频| 一二三四社区在线视频社区8| 人人妻人人澡欧美一区二区 | 日韩欧美一区视频在线观看| 日日干狠狠操夜夜爽| 成人三级做爰电影| 亚洲av成人一区二区三| 久久国产亚洲av麻豆专区| 欧美人与性动交α欧美精品济南到| 久久亚洲真实| 免费高清在线观看日韩| 亚洲熟妇熟女久久| 视频区欧美日本亚洲| 女性生殖器流出的白浆| 欧美激情 高清一区二区三区| 亚洲专区中文字幕在线| 国产成人免费无遮挡视频| 少妇的丰满在线观看| 制服诱惑二区| 一边摸一边抽搐一进一小说| 国产精品一区二区三区四区久久 | 午夜精品国产一区二区电影| 一边摸一边做爽爽视频免费| 成人国产一区最新在线观看| 制服人妻中文乱码| 老鸭窝网址在线观看| 99国产综合亚洲精品| 久久久精品欧美日韩精品| av超薄肉色丝袜交足视频| 日本vs欧美在线观看视频| 97碰自拍视频| 波多野结衣巨乳人妻| 国内精品久久久久精免费| 精品久久久久久成人av| 在线观看免费午夜福利视频| 亚洲人成77777在线视频| 人人妻人人爽人人添夜夜欢视频| 国产黄a三级三级三级人| 久久狼人影院| 一级黄色大片毛片| 国产aⅴ精品一区二区三区波| 免费看美女性在线毛片视频| 欧美老熟妇乱子伦牲交| 欧美av亚洲av综合av国产av| 波多野结衣av一区二区av| 国产私拍福利视频在线观看| 亚洲精品久久成人aⅴ小说| 亚洲自拍偷在线| 亚洲全国av大片| 村上凉子中文字幕在线| 少妇被粗大的猛进出69影院| 亚洲五月天丁香| 欧美日韩福利视频一区二区| 成年人黄色毛片网站| 免费女性裸体啪啪无遮挡网站| 久久久精品欧美日韩精品| 人人妻人人爽人人添夜夜欢视频| 午夜a级毛片| 亚洲第一av免费看| 丝袜人妻中文字幕| 国产国语露脸激情在线看| 男人操女人黄网站| 1024视频免费在线观看| 久久久久久久午夜电影| 国产欧美日韩一区二区三| 亚洲熟妇熟女久久| 国产区一区二久久| 久久久久久免费高清国产稀缺| 久久久久久大精品| 欧美国产日韩亚洲一区| 99久久久亚洲精品蜜臀av| 99国产精品免费福利视频| 最近最新免费中文字幕在线| 国产成人精品无人区| 国产高清有码在线观看视频 | 亚洲一区二区三区色噜噜| 亚洲国产欧美日韩在线播放| 国产精品香港三级国产av潘金莲| 国产成人系列免费观看| 韩国精品一区二区三区| 日本五十路高清| 成人永久免费在线观看视频| 超碰成人久久| 国产亚洲精品一区二区www| 亚洲天堂国产精品一区在线| 中文字幕色久视频| 黄片大片在线免费观看| 亚洲精品美女久久av网站| 久久婷婷人人爽人人干人人爱 | 波多野结衣av一区二区av| 成人欧美大片| 亚洲av片天天在线观看| 日韩精品中文字幕看吧| 亚洲国产高清在线一区二区三 | 亚洲专区中文字幕在线| 国产精品98久久久久久宅男小说| 欧美日本视频| 亚洲av五月六月丁香网| 欧美黑人精品巨大| 青草久久国产| 久久青草综合色| 变态另类成人亚洲欧美熟女 | 午夜日韩欧美国产| 老司机福利观看| 18禁黄网站禁片午夜丰满| 麻豆国产av国片精品| 中亚洲国语对白在线视频| 欧美色视频一区免费| 欧美一级a爱片免费观看看 | 亚洲国产精品合色在线| 两性午夜刺激爽爽歪歪视频在线观看 | aaaaa片日本免费| 欧美日韩福利视频一区二区| 国产色视频综合| 色综合站精品国产| 老汉色av国产亚洲站长工具| 亚洲片人在线观看| 男女下面进入的视频免费午夜 | 亚洲国产精品sss在线观看| 一二三四在线观看免费中文在| 久热爱精品视频在线9| 亚洲五月色婷婷综合| 少妇粗大呻吟视频| 久久中文字幕一级| 黄色丝袜av网址大全| 精品少妇一区二区三区视频日本电影| 国产精品一区二区三区四区久久 | 久久香蕉国产精品| 在线视频色国产色| 亚洲第一av免费看| 久久狼人影院| 老汉色av国产亚洲站长工具| av福利片在线| 久久久久久久久免费视频了| 久久久久久久久免费视频了| 精品久久蜜臀av无| 人妻久久中文字幕网| 亚洲男人的天堂狠狠| 久久影院123| 亚洲欧美精品综合久久99| 国产亚洲欧美精品永久| av片东京热男人的天堂| 日本撒尿小便嘘嘘汇集6| 视频在线观看一区二区三区| 90打野战视频偷拍视频| 夜夜爽天天搞| 亚洲精品国产一区二区精华液| 精品国产乱子伦一区二区三区| 叶爱在线成人免费视频播放| 国产免费av片在线观看野外av| 精品人妻在线不人妻| 亚洲一区二区三区色噜噜| 午夜a级毛片| 国产一区二区三区在线臀色熟女| 日本精品一区二区三区蜜桃| 在线免费观看的www视频| 日本 av在线| 啦啦啦 在线观看视频| 中文字幕另类日韩欧美亚洲嫩草| 国产成人影院久久av| 国产精品野战在线观看| 国产精品精品国产色婷婷| 一区二区三区精品91| 亚洲色图综合在线观看| 91麻豆精品激情在线观看国产| www.www免费av| 看免费av毛片| 亚洲欧美精品综合一区二区三区| 天天躁夜夜躁狠狠躁躁| 日本 av在线| 亚洲熟妇中文字幕五十中出| av超薄肉色丝袜交足视频| 99久久精品国产亚洲精品| 欧美乱色亚洲激情| 精品国产一区二区三区四区第35| 国产欧美日韩一区二区精品| 咕卡用的链子| 12—13女人毛片做爰片一| 亚洲男人天堂网一区| 一边摸一边抽搐一进一小说| 天堂影院成人在线观看| 亚洲天堂国产精品一区在线| 亚洲国产日韩欧美精品在线观看 | 纯流量卡能插随身wifi吗| 国产高清激情床上av| 大型av网站在线播放| 国产成人啪精品午夜网站| 变态另类成人亚洲欧美熟女 | 一级片免费观看大全| 非洲黑人性xxxx精品又粗又长| 亚洲情色 制服丝袜| 成人国产综合亚洲| 亚洲国产精品999在线| 满18在线观看网站| 亚洲国产精品久久男人天堂| 91在线观看av| 精品少妇一区二区三区视频日本电影| 午夜免费成人在线视频| 久久中文看片网| 久久久久久久精品吃奶| 国产97色在线日韩免费| 亚洲自偷自拍图片 自拍| 午夜a级毛片| 亚洲黑人精品在线| 两人在一起打扑克的视频| 亚洲色图综合在线观看| 一区二区三区激情视频| 1024香蕉在线观看| 母亲3免费完整高清在线观看| 日韩一卡2卡3卡4卡2021年| av电影中文网址| 日韩欧美国产一区二区入口| 国产91精品成人一区二区三区| 国产成人精品无人区| 一进一出抽搐动态| 国产精品一区二区免费欧美| 亚洲熟女毛片儿| 久久精品人人爽人人爽视色| 成人欧美大片| 88av欧美| 亚洲欧洲精品一区二区精品久久久| 啦啦啦免费观看视频1| 日韩有码中文字幕| 欧美绝顶高潮抽搐喷水| 色在线成人网| 国产一区二区三区综合在线观看| 亚洲精品久久国产高清桃花| 性欧美人与动物交配| 国产亚洲精品一区二区www| 国产精品影院久久| 满18在线观看网站| 精品卡一卡二卡四卡免费| 神马国产精品三级电影在线观看 | 国产成人av教育| 多毛熟女@视频| 黑人巨大精品欧美一区二区mp4| 亚洲激情在线av| 午夜精品在线福利| 少妇 在线观看| 国产高清有码在线观看视频 | 99国产精品99久久久久| 黄网站色视频无遮挡免费观看| 91字幕亚洲| 99精品欧美一区二区三区四区| 欧美成人性av电影在线观看| 91麻豆av在线| 身体一侧抽搐| 日本免费a在线| 精品熟女少妇八av免费久了| 亚洲精品国产色婷婷电影| 咕卡用的链子| 91成人精品电影| 久久中文字幕人妻熟女| 精品一区二区三区四区五区乱码| 免费在线观看视频国产中文字幕亚洲| 一区二区日韩欧美中文字幕| 一个人观看的视频www高清免费观看 | 久久国产精品人妻蜜桃| 欧美成人午夜精品| 亚洲av日韩精品久久久久久密| 极品人妻少妇av视频| 成人欧美大片| 亚洲 欧美 日韩 在线 免费| 国产主播在线观看一区二区| a在线观看视频网站| 中文字幕久久专区| 丝袜在线中文字幕| av在线天堂中文字幕| 免费在线观看黄色视频的| 国产区一区二久久| www.熟女人妻精品国产| 亚洲一区二区三区不卡视频| 无遮挡黄片免费观看| 色播亚洲综合网| 欧美一级毛片孕妇| 欧美 亚洲 国产 日韩一| 精品一区二区三区视频在线观看免费| 久久亚洲精品不卡| 久久久久久久久久久久大奶| xxx96com| 97碰自拍视频| 国产精品精品国产色婷婷| 可以免费在线观看a视频的电影网站| 91在线观看av| 老熟妇仑乱视频hdxx| 亚洲熟妇熟女久久| 首页视频小说图片口味搜索| 久久国产精品影院| 黄片小视频在线播放| 国产区一区二久久| 亚洲国产高清在线一区二区三 | 国内精品久久久久精免费| 成人特级黄色片久久久久久久| 久久久水蜜桃国产精品网| 欧美中文日本在线观看视频| 男人的好看免费观看在线视频 | 女人高潮潮喷娇喘18禁视频| 久久久久久大精品| 欧美黄色淫秽网站| 曰老女人黄片| 中文字幕人妻丝袜一区二区| 两人在一起打扑克的视频| 国产精品亚洲一级av第二区| 亚洲七黄色美女视频| 欧美精品啪啪一区二区三区| 久久精品91无色码中文字幕| 国产一级毛片七仙女欲春2 | 别揉我奶头~嗯~啊~动态视频| 精品久久久精品久久久| 自线自在国产av| 欧美日韩中文字幕国产精品一区二区三区 | 极品教师在线免费播放| 可以免费在线观看a视频的电影网站| 国产高清视频在线播放一区| 国产精品久久久久久人妻精品电影| 9191精品国产免费久久| 国产主播在线观看一区二区| 老鸭窝网址在线观看| 欧美大码av| 国产一级毛片七仙女欲春2 | 免费av毛片视频| 真人一进一出gif抽搐免费| 老司机深夜福利视频在线观看| 一夜夜www| 久久久久国产精品人妻aⅴ院| 母亲3免费完整高清在线观看| 国产视频一区二区在线看| 亚洲成av人片免费观看| 嫩草影视91久久| 99香蕉大伊视频| 激情在线观看视频在线高清| 欧美黑人欧美精品刺激| 国产亚洲欧美98| 亚洲欧美精品综合一区二区三区| 91精品三级在线观看| 一区在线观看完整版| 亚洲自偷自拍图片 自拍| 1024香蕉在线观看| 老司机午夜十八禁免费视频| 纯流量卡能插随身wifi吗| x7x7x7水蜜桃| 亚洲精品国产区一区二| 亚洲精品中文字幕一二三四区| 欧美中文综合在线视频| 国产精品免费视频内射| 琪琪午夜伦伦电影理论片6080| 亚洲人成伊人成综合网2020| 777久久人妻少妇嫩草av网站| av有码第一页| 黄片播放在线免费| 日本精品一区二区三区蜜桃| 精品久久久久久久毛片微露脸| 久久天堂一区二区三区四区| 在线观看日韩欧美| 黑人操中国人逼视频| 欧美黑人精品巨大| 无人区码免费观看不卡| 97人妻天天添夜夜摸| 免费看美女性在线毛片视频| 高潮久久久久久久久久久不卡| 日韩欧美一区二区三区在线观看| 高清毛片免费观看视频网站| 最新美女视频免费是黄的| 日韩中文字幕欧美一区二区| 18禁国产床啪视频网站| 18禁裸乳无遮挡免费网站照片 | svipshipincom国产片| 亚洲激情在线av| 国产在线精品亚洲第一网站| 国产主播在线观看一区二区| 日韩一卡2卡3卡4卡2021年| 欧美精品啪啪一区二区三区| 精品乱码久久久久久99久播| 男人操女人黄网站| 最新美女视频免费是黄的| 亚洲第一av免费看| 一级作爱视频免费观看| 九色亚洲精品在线播放| 最近最新中文字幕大全免费视频| 中文亚洲av片在线观看爽| 在线观看日韩欧美| 巨乳人妻的诱惑在线观看| 久久精品国产99精品国产亚洲性色 | 精品国内亚洲2022精品成人| www.自偷自拍.com| 天堂√8在线中文| 久久香蕉激情| 国产精品亚洲一级av第二区| 国内毛片毛片毛片毛片毛片| 窝窝影院91人妻| 黄片小视频在线播放| 欧美最黄视频在线播放免费| 丝袜美腿诱惑在线| 国产精品一区二区精品视频观看| 免费看美女性在线毛片视频| 欧美久久黑人一区二区| 国产精品久久久av美女十八| 久久香蕉激情| 免费在线观看视频国产中文字幕亚洲| 国产精品香港三级国产av潘金莲| av超薄肉色丝袜交足视频| 中文字幕精品免费在线观看视频| 男女床上黄色一级片免费看| 国产精品 欧美亚洲| 国产真人三级小视频在线观看| 午夜成年电影在线免费观看| 一a级毛片在线观看| 成人三级做爰电影| 天天躁狠狠躁夜夜躁狠狠躁| 免费看十八禁软件| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 长腿黑丝高跟| 级片在线观看| 国产成人欧美| 亚洲精华国产精华精| 亚洲五月色婷婷综合| 亚洲人成电影免费在线| 在线观看一区二区三区| 亚洲人成伊人成综合网2020| 嫩草影视91久久| 日韩av在线大香蕉| 校园春色视频在线观看| 欧美在线黄色| 美女扒开内裤让男人捅视频| 国产又色又爽无遮挡免费看| 性欧美人与动物交配| 亚洲精品一区av在线观看| 久久久久久亚洲精品国产蜜桃av| 在线视频色国产色| 亚洲一区高清亚洲精品| 免费看a级黄色片| 欧美黄色淫秽网站| 视频区欧美日本亚洲| 国产极品粉嫩免费观看在线| 黄色 视频免费看| 老鸭窝网址在线观看| 精品熟女少妇八av免费久了| 亚洲五月色婷婷综合| 亚洲五月天丁香| 国产精品一区二区精品视频观看| 亚洲欧美一区二区三区黑人| 欧美绝顶高潮抽搐喷水| 精品免费久久久久久久清纯| 亚洲色图综合在线观看| 19禁男女啪啪无遮挡网站| 国产在线精品亚洲第一网站| 在线观看免费视频网站a站| 亚洲一区高清亚洲精品| 亚洲少妇的诱惑av| 可以在线观看的亚洲视频| 亚洲国产欧美日韩在线播放| 午夜精品久久久久久毛片777| 变态另类成人亚洲欧美熟女 | 久久久久久亚洲精品国产蜜桃av| 精品久久久久久,| 少妇的丰满在线观看| 国产精品香港三级国产av潘金莲| 法律面前人人平等表现在哪些方面| 国产精品久久久久久亚洲av鲁大| 嫁个100分男人电影在线观看| 国产精品99久久99久久久不卡| 亚洲七黄色美女视频| 精品久久久久久久人妻蜜臀av | 国产精品自产拍在线观看55亚洲| 搡老熟女国产l中国老女人| 精品免费久久久久久久清纯| 老司机在亚洲福利影院| 日韩高清综合在线| 天天躁夜夜躁狠狠躁躁| 日韩一卡2卡3卡4卡2021年| 精品久久久久久,| 久久久久国产一级毛片高清牌| 亚洲国产精品sss在线观看| 别揉我奶头~嗯~啊~动态视频| 精品国产超薄肉色丝袜足j| 1024视频免费在线观看| 久热爱精品视频在线9| 嫁个100分男人电影在线观看| 日本黄色视频三级网站网址| 一本久久中文字幕| 妹子高潮喷水视频| 亚洲人成电影免费在线| 亚洲少妇的诱惑av| 国产在线精品亚洲第一网站| av视频在线观看入口| 亚洲全国av大片| av天堂久久9| 少妇裸体淫交视频免费看高清 | 最近最新免费中文字幕在线| 日韩精品中文字幕看吧| 久久久久久久久久久久大奶| 久久国产精品影院| 在线观看免费日韩欧美大片| 一区二区三区激情视频| 午夜福利视频1000在线观看 | 69av精品久久久久久| 中文字幕精品免费在线观看视频| 亚洲一码二码三码区别大吗| 999久久久精品免费观看国产| 色av中文字幕| 亚洲av成人av| 好看av亚洲va欧美ⅴa在| 男女下面插进去视频免费观看| 国产精品98久久久久久宅男小说| 母亲3免费完整高清在线观看| 在线观看66精品国产| 国产成人影院久久av| 999精品在线视频| 中文字幕人妻熟女乱码| 黄色 视频免费看| 日韩一卡2卡3卡4卡2021年| 国产又爽黄色视频| 久久伊人香网站| 亚洲av美国av| 看片在线看免费视频| 国产亚洲精品久久久久5区| 欧美一区二区精品小视频在线| 伦理电影免费视频| 精品国产乱码久久久久久男人| 国产高清视频在线播放一区| 国产人伦9x9x在线观看| 日日干狠狠操夜夜爽| 乱人伦中国视频|