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

    基于數(shù)值模擬的地表水體對地下水源地影響研究

    2020-04-06 04:10:36周瑞靜宋煒倪寶峰
    城市地質(zhì) 2020年4期
    關(guān)鍵詞:水質(zhì)

    周瑞靜 宋煒 倪寶峰

    摘 要:選取通州區(qū)納污河流運(yùn)潮減河和通州區(qū)區(qū)級水源地,運(yùn)用地下水水流模型Visual-Modflow和多組分溶質(zhì)運(yùn)移模型MT3DMSD耦合,以運(yùn)潮減河中濃度較高的耗氧量和氨氮為特征污染物,模擬了當(dāng)耗氧量和氨氮的初始濃度分別為20 mg/L、10 mg/L時,5年后耗氧量和氨氮在左岸潛水含水層的最大超標(biāo)范圍為153 m和176 m,對承壓水沒有影響。但在運(yùn)潮減河保持現(xiàn)狀水質(zhì)情況下,潛水污染面積會逐年增大并持續(xù)向下滲濾,對承壓水造成威脅,進(jìn)而影響水源地水質(zhì)質(zhì)量。

    關(guān)鍵詞:地表水;地下水;水源地;數(shù)值模擬

    Abstract: In order to explore the impact of inferior surface water on the water quality of groundwater drinking water, we take Tongzhou district sewage river and groundwater drinking water as an example, and use Visual MODFLOW and MT3DMSD to simulate the groundwater flow model. A groundwater flow model and a water quality model in the northern part of Tongzhou district are established. We take the chemical oxygen demand and ammonia nitrogen as the main research object, which contents are high in the Yunchaojian River. The maximum exceeding range of oxygen consumption and ammonia nitrogen in phreatic aquifer of the left bank is up to 94 m and 140 m after 5 years, under the initial concentration of 20 mg/L and 10 mg/L respectively, which has no impact on the confined water. However, under the condition of maintaining the current water quality of Yunchaojian River, the pollution area of phreatic aquifer will increase year by year. The continuous seepage of the polluted river water will threaten the confined aquifer, thus affect the quality of the underground source of drinking water in the confined aquifer.

    Keywords: surface water; groundwater; source of drinking water; numerical simulation

    0 引言

    首都北京是以地下水作為居民生產(chǎn)、生活主要供水水源的國際化大都市,2015年北京市將通州區(qū)定位為城市副中心,成為北京城市行政功能戰(zhàn)略轉(zhuǎn)移區(qū)域,肩負(fù)著承接中心城區(qū)的行政功能,并以此帶動其他功能疏解的使命。在為區(qū)域發(fā)展帶來機(jī)遇的同時,又對其可持續(xù)發(fā)展提出了更高的要求,也給該區(qū)地下水環(huán)境保護(hù)與管理帶來了更大的挑戰(zhàn)。目前,城市副中心區(qū)域第一含水層組地下水多個指標(biāo)出現(xiàn)超標(biāo),以無機(jī)指標(biāo)超標(biāo)為主,局部區(qū)域甚至呈現(xiàn)無機(jī)與有機(jī)組兩類污染雙重污染的狀況,地下水環(huán)境形勢嚴(yán)峻(郭高軒等,2014)。通州區(qū)地勢低凹,自古有“九河下梢”之稱,區(qū)內(nèi)地表水系水質(zhì)較差,多為劣Ⅴ類水質(zhì)(北京市水務(wù)局,2018)。運(yùn)潮減河作為通州區(qū)的納污河流,多年來河水來源主要是雨水直排、污水處理廠再生水排放等,水質(zhì)監(jiān)測結(jié)果多年為Ⅴ類,且河底沒有防滲襯層。

    李中陽等(2012)和馬闖等(2012)的研究結(jié)果顯示再生水灌溉會造成表層土壤重金屬濃度升高,深層土壤和含水層中重金屬濃度變化很小,表明符合灌溉標(biāo)準(zhǔn)的再生水不會對地下水造成重金屬污染,但是聚集在表層土壤中的重金屬對地下水的潛在污染風(fēng)險也不可忽視。陳衛(wèi)平等(2013)研究了不同回灌方式下再生水對地下水的污染風(fēng)險,重金屬的研究結(jié)果與李中陽和馬闖等的結(jié)論一致,均為對地下水的影響較小,而鹽分和硝態(tài)氮在地下水中都表現(xiàn)出明顯的增加趨勢。賈忠華等(2018)研究了雨水花園入滲雨水對地下水水質(zhì)的影響,結(jié)果表明地下水中總氮在短期內(nèi)顯著增加,氨氮增加不顯著,硝態(tài)氮在雨后有明顯升高。李賀強(qiáng)等(2015)用數(shù)值模擬法研究了深井雨水回灌對深層地下水的水質(zhì)影響,結(jié)果表明回灌水量和溶質(zhì)濃度均對溶質(zhì)的遷移擴(kuò)散有影響。杜新強(qiáng)等(2007)研究了地表水注入對地下水水質(zhì)的影響,并模擬了地表水進(jìn)入含水層后發(fā)生的水-巖相互作用,得出了優(yōu)質(zhì)地表水注入補(bǔ)給地下水可以改善劣質(zhì)地下含水層水質(zhì)。谷小溪(2012)利用人工試驗場研究了人工回灌對地下水水質(zhì)的影響,總結(jié)了人工回灌下地下水中8大離子和5種微量元素受混合作用、離子交換作用及含水層巖石溶解作用的影響程度,得出地下水具有一定的緩沖作用,自來水回灌地下水不會對地下水環(huán)境造成影響。

    本次研究使用地下水模擬領(lǐng)域廣泛應(yīng)用的Visual Modflow軟件系統(tǒng),該系統(tǒng)在整合了水流評價、平面和剖面流線示蹤分析、水文地質(zhì)參數(shù)估計與優(yōu)化和溶質(zhì)運(yùn)移評價等軟件基礎(chǔ)上,開發(fā)出了強(qiáng)大的圖形可視界面功能,是一款被國內(nèi)外各機(jī)構(gòu)認(rèn)可的地下水流場和溶質(zhì)運(yùn)移三維綜合評價軟件系統(tǒng)(王慶永等,2007;馮潔,2013;蔚東升等,2014;趙貝等,2015)。已有運(yùn)用Visual Modflow模擬再生水長期排入地表水體對當(dāng)?shù)氐叵滤|(zhì)量影響進(jìn)行評價的研究報道(崔瑜等,2018)。

    通州區(qū)唯一的區(qū)級飲用水水源地部分水源井傍河而建,取第四系潛水和微承壓水,其中兩眼水源井緊鄰運(yùn)潮減河,水源地的大規(guī)模取水也加劇了該區(qū)域地表水與地下水的相互影響。本研究為了探究劣質(zhì)地表水對地下水水源地水質(zhì)的影響,選取通州區(qū)納污河流運(yùn)潮減河和通州區(qū)區(qū)級水源地為例,運(yùn)用地下水水流模型Visual-Modflow和多組分溶質(zhì)運(yùn)移模型MT3DMSD耦合,以運(yùn)潮減河中濃度較高的耗氧量和氨氮為主要研究對象,模擬了區(qū)級水源地現(xiàn)狀取水量下運(yùn)潮減河中Ⅴ類地表水對區(qū)級水源地的影響程度和范圍。

    1 模擬區(qū)信息

    1.1 水源地概況

    通州區(qū)區(qū)級水源地位于通州區(qū)北部,有水源井28眼,其中第四系水源井13眼,位于北京與河北交界,鄰潮白河,基巖井15眼,位于城區(qū)北側(cè)。該水源地設(shè)計取水量1250萬t/a,供水對象包括通州城區(qū)4個街道辦事處全部區(qū)域以及永順鎮(zhèn)、梨園鎮(zhèn)和宋莊鎮(zhèn)部分區(qū)域,服務(wù)人口達(dá)25萬人。第四系水源井南側(cè)有運(yùn)潮減河,東側(cè)有潮白河,多個水源井傍河而建(北京市地質(zhì)工程勘察院,2016年)。

    1.2 第四系含水層特征

    模擬區(qū)第四系含水層特點(diǎn)是多層次,細(xì)顆粒,巖性以中、細(xì)、粉砂為主,局部地區(qū)以粗砂為主,夾有少量礫石。根據(jù)區(qū)域水文地質(zhì)條件,綜合分析隔水層的相對厚度及巖性、含水層組之間水力聯(lián)系和水質(zhì)等因素,可將模擬區(qū)含水層在垂向上概化為4個含水層組:

    潛水含水層組以中細(xì)砂為主,其次為細(xì)砂、粉砂,局部地區(qū)含粗砂、礫石。含水層一般有2~3層,累計厚度20~30 m,底板埋深40~50 m左右,地下水位埋深6~7 m。本含水層組富水性差,水位下降5 m時,單井出水量在500~1500 m3/d。

    第一承壓含水層組以細(xì)砂為主,其次為中砂、粗砂,局部地區(qū)含礫石。含水層一般有3~5層,累計厚度30~50 m,底板埋深80~100 m左右,地下水位埋深18~20 m。在通州區(qū)中北部地區(qū)該含水層富水性較好,水位下降5 m時,單井出水量在1500~3000 m3/d。

    第二承壓含水層組(底界深度150~180 m)和第三承壓含水層組(底界深度300 m)主要以細(xì)砂為主,中砂和粉砂次之,地下水水位埋深50 m左右。雖然富水性和補(bǔ)給條件較淺層地下水差,但由于埋藏較深,目前受到的污染較少,水質(zhì)較淺層地下水好。

    1.3 地下水質(zhì)量現(xiàn)狀

    區(qū)級水源地附近有區(qū)域地下水分層監(jiān)測井3眼,井深為41 m、103 m、225 m,分別監(jiān)測潛水、第一層承壓水、第二層承壓水。據(jù)監(jiān)測結(jié)果潛水含水層總硬度和溶解性總固體超標(biāo),耗氧量和氨氮不超標(biāo);而第一承壓含水層和第二承壓含水層地下水水質(zhì)中總硬度和溶解性總固體相對較好,耗氧量和氨氮不超標(biāo)(表1)。

    2 研究方法

    2.1含水層結(jié)構(gòu)概化

    通州區(qū)區(qū)級水源地位于潮白河沖洪積扇中下部,第四系沉積在模擬區(qū)范圍內(nèi)為300~500 m。由于本次研究地表水體對水源井水質(zhì)的影響,綜合考慮含水層概化需求及水井深度,模型在垂向上控制在300 m以內(nèi)。模擬區(qū)第四系含水層在300 m深度以內(nèi)有多個含水層,且各含水層之間存在厚度不等的粉質(zhì)黏土層,含水層巖性主要為多層薄層中細(xì)砂,局部含卵礫石層。垂向上將模型概化為3個含水層、2個弱透水層,共5層:第一層為潛水含水層,深度在40 m左右,巖性以細(xì)砂為主;第二層為弱透水層,厚度在10~20 m,巖性以粉質(zhì)黏土和黏土為主;第三層為淺層承壓含水層,深度在80~100 m,巖性主要為細(xì)砂,其次為中砂和粗砂;第四層為弱透水層,深度在100~120 m,主要為厚20~40 m的粉質(zhì)黏土與黏土互層,為該區(qū)域第一層厚度相對穩(wěn)定的隔水含水層;第五層為深層承壓含水層組,埋深在120 m以下,主要由細(xì)砂組成,中間夾粉質(zhì)黏土,由于農(nóng)村改水工程的開展,該層含水層開采量逐漸增加,現(xiàn)已成為主要的開采層。

    2.2 邊界條件概化

    側(cè)向邊界:模擬區(qū)東部為潮白河,通州區(qū)內(nèi)潮白河河段常年有水,河流底部未防滲,定義為河流邊界;西部為溫榆河和北運(yùn)河,通州區(qū)內(nèi)該河段常年有水,河流底部進(jìn)行了防滲處理,定義為河流邊界;根據(jù)地下水流場方向,沿地下水等水位線平行方向,將模擬區(qū)南部南劉各莊村一線,定義為側(cè)向流出邊界;北部以通州區(qū)區(qū)界與等水位線相切處為界,定義為零通量邊界;西北部定義為補(bǔ)給邊界(圖1)。

    垂向邊界:模擬區(qū)潛水埋深較淺,直接接受大氣降水和田間灌溉入滲補(bǔ)給,因此定義潛水自由水面為上邊界。模型底部為黏土層夾粉質(zhì)黏土層,為隔水含水層,定義為隔水邊界。

    2.3 水流模型建立

    模型中潛水含水層、隔水層、承壓水含水層的剖分均采用網(wǎng)格剖分法,并對水源地、農(nóng)村改水井所在位置進(jìn)行加密剖分,最終每層剖分為48765個單元格。

    模擬驗證期為2012年6月至2016年6月,以一個月作為一個時間段,識別模型后,地下水流場擬合圖和典型地下水動態(tài)曲線擬合圖結(jié)果見圖1和圖2,表明本次模擬結(jié)構(gòu)合理。在此基礎(chǔ)上對運(yùn)潮減河水質(zhì)維持現(xiàn)狀Ⅴ類的情況對地下水源地水質(zhì)的影響進(jìn)行模擬預(yù)測。

    2.4? MT3DMS溶質(zhì)運(yùn)移模型

    (1)模型構(gòu)建

    通過建立地下水溶質(zhì)運(yùn)移模型模擬污染物的運(yùn)移。此處考慮最不利情況,假定污染物已到達(dá)潛水含水層頂板,且達(dá)到最大濃度,以河水中各污染物的濃度作為初始源強(qiáng)計算,在水文地質(zhì)概念模型的基礎(chǔ)上預(yù)測污染物在地下水中的運(yùn)移規(guī)律。

    溶質(zhì)運(yùn)移模型的范圍與地下水流數(shù)值模型相一致,邊界分別對應(yīng)于水流模型的邊界。由于運(yùn)潮減河河道水質(zhì)為劣Ⅴ類,耗氧量和氨氮濃度均高于地下水,將其概化為定濃度線源補(bǔ)給邊界。

    (2)初始濃度設(shè)定

    運(yùn)潮減河位于通州區(qū)區(qū)級水源地二級保護(hù)區(qū)內(nèi),根據(jù)檢測結(jié)果(表2),運(yùn)潮減河河水水質(zhì)為劣Ⅴ類,對水源地存在污染的風(fēng)險。根據(jù)建立的數(shù)值模擬模型,對運(yùn)潮減河對水源地影響進(jìn)行預(yù)測評價。按照偏保守的評價原則,模擬地表水氨氮、耗氧量初始濃度分別設(shè)定為20mg/L、10mg/L。

    2.5 結(jié)果與分析

    根據(jù)預(yù)測結(jié)果可知,運(yùn)潮減河河水中超標(biāo)污染物可以進(jìn)入地下水環(huán)境中,隨時間增加污染物影響范圍逐漸擴(kuò)大,預(yù)測河水現(xiàn)狀水質(zhì)下,5年后耗氧量和氨氮在左岸潛水含水層的最大超標(biāo)范圍為153 m和176 m,對承壓水沒有影響,可見地水質(zhì)差的地表水體對地下水潛水含水層水質(zhì)有一定影響,短時間內(nèi)對承壓水無影響,見圖3至圖8。

    水源10#和水源19#臨運(yùn)潮減河左岸,水源10#取水層位為第二層承壓水,水源19#取水層位為第一層承壓水及第二層承壓水,由于水源地范圍內(nèi)潛水與第一層承壓水之間的弱透水層不夠連貫,不排除局部存在透鏡體的可能,因此潛水與第一承壓水之間水力聯(lián)系較密切,隨著水源地持續(xù)開采,地下水水位下降,井周邊水力梯度增加,將加劇地下水的徑流速度,從而加劇水體中污染物質(zhì)遷移速度。因此,運(yùn)潮減河河水對當(dāng)?shù)氐叵滤挠绊懖蝗莺鲆?,尤其需要重視其對水?9#的影響。需加強(qiáng)對水源地源水的監(jiān)測,密切關(guān)注監(jiān)測指標(biāo)的變化趨勢。同時應(yīng)對通州區(qū)內(nèi)水質(zhì)為劣Ⅴ類的地表水體進(jìn)行集中整治,必要時進(jìn)行河底淤泥研究治理,鋪設(shè)防滲襯層。

    3 結(jié)論

    (1)模擬區(qū)潛水取樣點(diǎn)水質(zhì)檢測結(jié)果顯示兩期潛水耗氧量均超標(biāo),氨氮含量未超標(biāo),按照地下水質(zhì)量分類屬于Ⅲ類,但是濃度從深層承壓水至潛水呈上升趨勢。

    (2)運(yùn)潮減河現(xiàn)狀水質(zhì)較差,按照耗氧量20 mg/L、氨氮10 mg/L的初始濃度進(jìn)行預(yù)測,5年后耗氧量和氨氮在左岸潛水含水層的最大超標(biāo)范圍為153 m和176 m,對承壓水含水層水質(zhì)沒有影響,可見水質(zhì)差的地表水體對潛水含水層水質(zhì)有一定影響,短時間內(nèi)對承壓水水質(zhì)無影響。

    (3)模擬區(qū)第一承壓水和潛水含水層之間的水力聯(lián)系緊密,運(yùn)潮減河保持現(xiàn)狀水質(zhì)情況下,潛水污染面積會逐年增大,甚至對承壓水造成污染,進(jìn)而影響水源地水質(zhì)。

    參考文獻(xiàn):

    北京市地質(zhì)工程勘察院,2016.北京市通州區(qū)城市地下飲用水水源地保護(hù)區(qū)調(diào)整與劃分[R].

    北京市水務(wù)局,2018.2017年北京市水務(wù)統(tǒng)計年鑒[R]..

    陳衛(wèi)平,呂斯丹,王美娥,等,2013.再生水回灌對地下水水質(zhì)影響研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報,24(5):1253-1262.

    崔瑜,吳泰然,周磊,等,2018.基于數(shù)值模擬的再生水對地下水質(zhì)影響實例[J].地質(zhì)通報,37(7):1365-1373.

    杜新強(qiáng),齊素文,廖資生,等,2007.人工補(bǔ)給對含水層水質(zhì)的影響[J].吉林大學(xué)學(xué)報(地球科學(xué)版)(2):293-297.

    谷小溪,2012.人工回灌條件下地下水中多組分溶質(zhì)的反應(yīng)遷移模擬:以上海市寶山顧村人工回灌試驗場地為例[D].吉林大學(xué):112.

    馮潔,2013.可視化地下水?dāng)?shù)值模擬軟件(VISUAL MODFLOW)在國內(nèi)的應(yīng)用[J].地下水,35(4):34-36.

    郭高軒,琚宜文,翟航,等,2014.北京市通州區(qū)地下水分層質(zhì)量評價及水化學(xué)特征[J].環(huán)境科學(xué),35(6):2114-2119.

    賈忠華,吳舒然,唐雙成,等,2018.雨水花園集中入滲對地下水水位與水質(zhì)的影響[J].水科學(xué)進(jìn)展,29(2):221-229.

    李中陽,樊向陽,齊學(xué)斌,等,2012. 再生水灌溉下重金屬在植物-土壤-地下水系統(tǒng)遷移的研究進(jìn)展[J].中國農(nóng)村水利水電(7):5-8.

    李賀強(qiáng),邢國平,王超,2015.雨水回灌對深層承壓含水層水質(zhì)影響的模擬研究[J].水土保持通報,35(1):139-142.

    馬闖,楊軍,雷梅,等,2012.北京市再生水灌溉對地下水的重金屬污染風(fēng)險[J].地理研究,31(12):2250-2258.

    王慶永,賈忠華,劉曉峰,等,2007.Visual MODFLOW及其在地下水模擬中的應(yīng)用[J].水資源與水工程學(xué)報(5):90-92.

    蔚東升,劉潔,馮兆東,2014.Visual MODFLOW在國內(nèi)的應(yīng)用[J].能源與節(jié)能(6):96-97,126.

    趙貝,劉建軍,2015.基于Visual Modflow地下水污染物溶質(zhì)運(yùn)移的模擬[J].石河子大學(xué)學(xué)報(自然科學(xué)版),33(6):787-792.

    猜你喜歡
    水質(zhì)
    水質(zhì)抽檢豈容造假
    環(huán)境(2023年5期)2023-06-30 01:20:01
    水質(zhì)檢測員——中華秋沙鴨
    水質(zhì)凈化廠提標(biāo)至一級A設(shè)計與運(yùn)行效果探討
    關(guān)于水質(zhì)監(jiān)測對環(huán)境保護(hù)的意義
    一月冬棚養(yǎng)蝦常見水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
    這條魚供不應(yīng)求!蝦蟹養(yǎng)殖戶、垂釣者的最愛,不用投喂,還能凈化水質(zhì)
    圖像識別在水質(zhì)檢測中的應(yīng)用
    電子制作(2018年14期)2018-08-21 01:38:16
    淺析黑臭水體成因、治理方法及水質(zhì)長效改善保持問題——水質(zhì)長效改善保持問題
    濟(jì)下水庫徑流水質(zhì)和垂向水質(zhì)分析及評價
    水質(zhì)的年輪——讀《時光的年輪》
    老师上课跳d突然被开到最大视频| 成人国产麻豆网| 午夜日本视频在线| 国产乱人视频| 欧美日本亚洲视频在线播放| 超碰av人人做人人爽久久| 熟女人妻精品中文字幕| 欧美极品一区二区三区四区| 国产久久久一区二区三区| 搞女人的毛片| 免费看a级黄色片| 久久久久久国产a免费观看| 国产精华一区二区三区| 成人综合一区亚洲| 亚洲欧美日韩东京热| 亚洲真实伦在线观看| 成年女人永久免费观看视频| 狂野欧美白嫩少妇大欣赏| 精品久久久噜噜| 国产老妇女一区| 我的老师免费观看完整版| 日本与韩国留学比较| 亚洲av成人精品一二三区| 免费播放大片免费观看视频在线观看 | 国产午夜精品久久久久久一区二区三区| 蜜臀久久99精品久久宅男| 人人妻人人澡欧美一区二区| 午夜a级毛片| 亚洲综合色惰| 亚洲伊人久久精品综合 | 免费观看的影片在线观看| 在线免费观看不下载黄p国产| 色视频www国产| 亚洲不卡免费看| 免费观看精品视频网站| 亚洲精华国产精华液的使用体验| 青春草视频在线免费观看| 亚洲无线观看免费| 最近中文字幕2019免费版| 欧美一区二区精品小视频在线| 亚洲欧美精品专区久久| 国产成人免费观看mmmm| 乱人视频在线观看| 97热精品久久久久久| 最近的中文字幕免费完整| 性色avwww在线观看| 午夜日本视频在线| 亚洲成人中文字幕在线播放| 男女下面进入的视频免费午夜| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费大片18禁| 男人和女人高潮做爰伦理| 麻豆乱淫一区二区| 级片在线观看| 日韩欧美在线乱码| 亚洲av二区三区四区| 天天躁夜夜躁狠狠久久av| 婷婷色综合大香蕉| 久久久久久久午夜电影| 欧美激情在线99| 少妇裸体淫交视频免费看高清| 欧美xxxx黑人xx丫x性爽| 亚洲成人久久爱视频| 亚洲精品乱码久久久v下载方式| 国产精品电影一区二区三区| 黄色一级大片看看| 亚洲精品aⅴ在线观看| 哪个播放器可以免费观看大片| 国产欧美另类精品又又久久亚洲欧美| 欧美一区二区亚洲| 国产精华一区二区三区| 国产精品久久久久久久久免| 久久这里有精品视频免费| 黄色一级大片看看| 亚洲av免费在线观看| 亚洲最大成人手机在线| 日韩一区二区视频免费看| 欧美一区二区精品小视频在线| 欧美一区二区国产精品久久精品| 91精品国产九色| 久久久久国产网址| eeuss影院久久| 久久韩国三级中文字幕| 高清毛片免费看| 亚洲成色77777| 美女国产视频在线观看| av女优亚洲男人天堂| 可以在线观看毛片的网站| 22中文网久久字幕| 日本av手机在线免费观看| 亚洲欧美精品综合久久99| 国产伦一二天堂av在线观看| www.色视频.com| 日韩高清综合在线| 看黄色毛片网站| 国产亚洲5aaaaa淫片| 午夜亚洲福利在线播放| 亚洲国产精品国产精品| 久久久久久九九精品二区国产| 久久久久久大精品| 成人亚洲精品av一区二区| 夫妻性生交免费视频一级片| 国产成人福利小说| 国产精品伦人一区二区| 波野结衣二区三区在线| 亚洲自偷自拍三级| 久久精品夜色国产| 麻豆成人av视频| 少妇人妻一区二区三区视频| 高清av免费在线| 亚洲va在线va天堂va国产| 欧美色视频一区免费| 国产精品美女特级片免费视频播放器| 中文字幕免费在线视频6| 直男gayav资源| 男女那种视频在线观看| 日韩,欧美,国产一区二区三区 | 丰满少妇做爰视频| 国产精品伦人一区二区| 天堂av国产一区二区熟女人妻| 色网站视频免费| 成人三级黄色视频| 久久久久久久午夜电影| 国产亚洲av嫩草精品影院| 国产91av在线免费观看| 18+在线观看网站| 国产极品精品免费视频能看的| 久久99热这里只有精品18| 日本免费a在线| 最近2019中文字幕mv第一页| 亚洲国产精品合色在线| 亚洲av成人av| 久久人人爽人人片av| 两性午夜刺激爽爽歪歪视频在线观看| 97超视频在线观看视频| 久久人人爽人人片av| 九草在线视频观看| 长腿黑丝高跟| 18禁动态无遮挡网站| 在线观看66精品国产| 爱豆传媒免费全集在线观看| 国产69精品久久久久777片| 国产毛片a区久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 黑人高潮一二区| 国产黄色视频一区二区在线观看 | 国模一区二区三区四区视频| 中文字幕av在线有码专区| 亚洲自拍偷在线| 国产精品久久视频播放| 亚洲精华国产精华液的使用体验| 99久久中文字幕三级久久日本| 麻豆成人午夜福利视频| 成年版毛片免费区| 综合色av麻豆| 亚洲图色成人| 国产精品国产三级专区第一集| 两个人视频免费观看高清| 国产大屁股一区二区在线视频| 国产女主播在线喷水免费视频网站 | 91在线精品国自产拍蜜月| av黄色大香蕉| 亚洲精品久久久久久婷婷小说 | 看黄色毛片网站| 看非洲黑人一级黄片| 国内精品宾馆在线| 大香蕉久久网| a级毛片免费高清观看在线播放| 老司机福利观看| 插逼视频在线观看| 我要搜黄色片| 国产色婷婷99| 日韩精品有码人妻一区| 精品欧美国产一区二区三| 91精品一卡2卡3卡4卡| 婷婷六月久久综合丁香| 国内精品宾馆在线| 观看免费一级毛片| 亚洲经典国产精华液单| 人人妻人人澡欧美一区二区| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久大av| 69av精品久久久久久| 一级av片app| 欧美日韩国产亚洲二区| 国产三级中文精品| 床上黄色一级片| 精品一区二区三区视频在线| 亚洲在线自拍视频| 欧美日韩一区二区视频在线观看视频在线 | 日韩av在线大香蕉| 精品久久久久久成人av| 国语自产精品视频在线第100页| 精品午夜福利在线看| 亚洲欧美一区二区三区国产| 欧美不卡视频在线免费观看| 国产亚洲最大av| 亚洲美女搞黄在线观看| 亚洲激情五月婷婷啪啪| 国内精品宾馆在线| 亚洲av二区三区四区| 97在线视频观看| 精品久久久久久成人av| 蜜臀久久99精品久久宅男| 美女脱内裤让男人舔精品视频| 黑人高潮一二区| 超碰97精品在线观看| 日本黄色视频三级网站网址| 观看美女的网站| 精品久久久久久久久久久久久| 欧美又色又爽又黄视频| 亚洲欧美日韩高清专用| 一个人看的www免费观看视频| 欧美一区二区国产精品久久精品| 蜜桃久久精品国产亚洲av| 欧美一区二区精品小视频在线| 日日摸夜夜添夜夜爱| 秋霞伦理黄片| 搡女人真爽免费视频火全软件| 最近的中文字幕免费完整| 波野结衣二区三区在线| 69人妻影院| 大又大粗又爽又黄少妇毛片口| 亚洲精品日韩av片在线观看| 韩国av在线不卡| 日本黄色视频三级网站网址| 性插视频无遮挡在线免费观看| 成人午夜精彩视频在线观看| 97超视频在线观看视频| 久久久久久久久久久丰满| 亚洲精品久久久久久婷婷小说 | 久久久久国产网址| 成人漫画全彩无遮挡| 国产色爽女视频免费观看| 国产真实伦视频高清在线观看| 毛片女人毛片| 国产精品一区二区在线观看99 | 男女视频在线观看网站免费| 欧美日韩国产亚洲二区| 成人特级av手机在线观看| 国产精品永久免费网站| av国产免费在线观看| 精品人妻熟女av久视频| 国产极品精品免费视频能看的| av视频在线观看入口| 可以在线观看毛片的网站| 草草在线视频免费看| a级毛色黄片| 国产极品精品免费视频能看的| 两个人视频免费观看高清| 天堂av国产一区二区熟女人妻| 久久久久久伊人网av| 又爽又黄无遮挡网站| 纵有疾风起免费观看全集完整版 | 久久久久网色| 国产老妇女一区| 精品免费久久久久久久清纯| 91精品一卡2卡3卡4卡| 久久久久久久久久黄片| 在线播放国产精品三级| 欧美日韩在线观看h| av女优亚洲男人天堂| 成年版毛片免费区| 青春草亚洲视频在线观看| 欧美一级a爱片免费观看看| 免费播放大片免费观看视频在线观看 | a级一级毛片免费在线观看| 亚洲国产精品专区欧美| 亚洲av电影不卡..在线观看| 深夜a级毛片| 狂野欧美白嫩少妇大欣赏| 夜夜看夜夜爽夜夜摸| 日日干狠狠操夜夜爽| 亚洲人成网站在线观看播放| 国产精品三级大全| 联通29元200g的流量卡| 精品午夜福利在线看| 日本三级黄在线观看| 亚洲自拍偷在线| 精品人妻熟女av久视频| 日韩欧美 国产精品| 日日摸夜夜添夜夜爱| 日日摸夜夜添夜夜添av毛片| 亚洲精品国产成人久久av| 一级爰片在线观看| 亚洲伊人久久精品综合 | 久久欧美精品欧美久久欧美| 亚洲成人久久爱视频| 全区人妻精品视频| 亚洲综合精品二区| 成人亚洲欧美一区二区av| 亚洲av一区综合| 成人一区二区视频在线观看| 日韩欧美精品v在线| 美女cb高潮喷水在线观看| 青春草视频在线免费观看| 亚洲av日韩在线播放| 国产黄a三级三级三级人| 一级黄片播放器| 成人高潮视频无遮挡免费网站| 七月丁香在线播放| 人妻少妇偷人精品九色| 亚洲,欧美,日韩| 国产成人免费观看mmmm| 国产真实乱freesex| 久久久久免费精品人妻一区二区| 男女下面进入的视频免费午夜| 亚洲四区av| 麻豆成人av视频| 美女国产视频在线观看| 亚洲人与动物交配视频| 国产精品国产高清国产av| 中文精品一卡2卡3卡4更新| 人妻系列 视频| 最近中文字幕2019免费版| 九九热线精品视视频播放| 日本三级黄在线观看| 男人舔奶头视频| 成人性生交大片免费视频hd| 91久久精品电影网| 两个人视频免费观看高清| 精品国内亚洲2022精品成人| 国产精品无大码| 一本久久精品| 人人妻人人看人人澡| 午夜a级毛片| 成人av在线播放网站| 日本熟妇午夜| 美女黄网站色视频| 国产精品伦人一区二区| 夫妻性生交免费视频一级片| 麻豆国产97在线/欧美| 女的被弄到高潮叫床怎么办| 你懂的网址亚洲精品在线观看 | 老司机福利观看| 国产熟女欧美一区二区| 国产国拍精品亚洲av在线观看| 三级国产精品片| 国产精品.久久久| 天堂√8在线中文| 一级av片app| 日本-黄色视频高清免费观看| 中文字幕熟女人妻在线| 成年av动漫网址| av.在线天堂| 一级二级三级毛片免费看| 直男gayav资源| 男人舔奶头视频| 18禁在线无遮挡免费观看视频| 秋霞伦理黄片| 久久这里只有精品中国| av女优亚洲男人天堂| 国产精品.久久久| 搡老妇女老女人老熟妇| 久久久久国产网址| 性插视频无遮挡在线免费观看| 99久久九九国产精品国产免费| 狠狠狠狠99中文字幕| 内射极品少妇av片p| 一级二级三级毛片免费看| 尾随美女入室| 久久精品久久精品一区二区三区| 白带黄色成豆腐渣| av国产久精品久网站免费入址| 看十八女毛片水多多多| 午夜老司机福利剧场| 又粗又爽又猛毛片免费看| 亚洲国产精品成人久久小说| 男的添女的下面高潮视频| 午夜爱爱视频在线播放| 午夜福利在线在线| 国产单亲对白刺激| 最新中文字幕久久久久| 内地一区二区视频在线| 99久久中文字幕三级久久日本| 国产又黄又爽又无遮挡在线| 久久人人爽人人爽人人片va| 成人午夜精彩视频在线观看| 99九九线精品视频在线观看视频| 亚洲最大成人中文| 好男人视频免费观看在线| 老司机影院毛片| 六月丁香七月| 在线a可以看的网站| 国产精品久久视频播放| 永久网站在线| 欧美xxxx性猛交bbbb| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 中文字幕av成人在线电影| 国产精品美女特级片免费视频播放器| av在线观看视频网站免费| 乱系列少妇在线播放| 亚洲激情五月婷婷啪啪| 欧美丝袜亚洲另类| 神马国产精品三级电影在线观看| 亚洲一区高清亚洲精品| 男人舔女人下体高潮全视频| 高清在线视频一区二区三区 | 欧美成人精品欧美一级黄| 亚洲成人中文字幕在线播放| 直男gayav资源| 最近最新中文字幕免费大全7| 丝袜喷水一区| 国产亚洲午夜精品一区二区久久 | 黄色欧美视频在线观看| 欧美性猛交黑人性爽| 日日摸夜夜添夜夜添av毛片| 日日摸夜夜添夜夜爱| .国产精品久久| 一区二区三区乱码不卡18| 久久精品91蜜桃| 国产精品嫩草影院av在线观看| 午夜福利高清视频| 超碰97精品在线观看| 国产不卡一卡二| 久久久久久久久久久丰满| 最近最新中文字幕免费大全7| 黄色欧美视频在线观看| 噜噜噜噜噜久久久久久91| 久久亚洲国产成人精品v| 日本av手机在线免费观看| 国产精品.久久久| 性色avwww在线观看| 精品一区二区三区人妻视频| 久久久久性生活片| 男人舔女人下体高潮全视频| 国产精品麻豆人妻色哟哟久久 | 秋霞在线观看毛片| 亚洲成人中文字幕在线播放| 亚洲天堂国产精品一区在线| 久久精品影院6| 中文字幕av在线有码专区| 我要搜黄色片| 国产高清不卡午夜福利| 国产精品嫩草影院av在线观看| 日韩一本色道免费dvd| 啦啦啦啦在线视频资源| 99久久精品热视频| 免费人成在线观看视频色| 最近视频中文字幕2019在线8| 亚洲,欧美,日韩| 99久久无色码亚洲精品果冻| 国产免费福利视频在线观看| 三级国产精品片| 欧美bdsm另类| 国产精品麻豆人妻色哟哟久久 | 中文字幕精品亚洲无线码一区| 性插视频无遮挡在线免费观看| 欧美潮喷喷水| 午夜亚洲福利在线播放| 三级男女做爰猛烈吃奶摸视频| 日本一二三区视频观看| 中文字幕av成人在线电影| 婷婷色av中文字幕| 亚洲综合色惰| 看片在线看免费视频| 一级二级三级毛片免费看| 色综合色国产| 美女cb高潮喷水在线观看| 久久久亚洲精品成人影院| 免费在线观看成人毛片| 久久精品国产自在天天线| 最近最新中文字幕大全电影3| 久久精品影院6| 尾随美女入室| 国模一区二区三区四区视频| 国产乱人视频| 亚洲欧美精品综合久久99| 国产精品一及| 看黄色毛片网站| av在线播放精品| 美女脱内裤让男人舔精品视频| 简卡轻食公司| 欧美日韩精品成人综合77777| 99久久成人亚洲精品观看| 99在线人妻在线中文字幕| 免费av观看视频| 国产免费视频播放在线视频 | 国产国拍精品亚洲av在线观看| 久久久久免费精品人妻一区二区| 国产精品久久视频播放| 日韩人妻高清精品专区| 中文字幕免费在线视频6| 国产成年人精品一区二区| 久久久久久久久中文| 最近中文字幕高清免费大全6| 欧美区成人在线视频| 亚洲成人久久爱视频| 精华霜和精华液先用哪个| 国产综合懂色| 亚洲最大成人av| 九草在线视频观看| 好男人在线观看高清免费视频| 波多野结衣高清无吗| 免费黄网站久久成人精品| 久久人人爽人人爽人人片va| 天天一区二区日本电影三级| 大香蕉97超碰在线| 亚洲成人精品中文字幕电影| 国产色爽女视频免费观看| 激情 狠狠 欧美| 熟妇人妻久久中文字幕3abv| 日本一二三区视频观看| 久久99蜜桃精品久久| av线在线观看网站| 欧美性猛交黑人性爽| 亚洲色图av天堂| 午夜爱爱视频在线播放| 久久精品国产自在天天线| 日韩视频在线欧美| 免费观看性生交大片5| 亚洲18禁久久av| 乱码一卡2卡4卡精品| 亚州av有码| 成人综合一区亚洲| 国产不卡一卡二| 婷婷六月久久综合丁香| 亚洲人成网站在线播| 美女国产视频在线观看| 国产色婷婷99| 国产黄a三级三级三级人| 亚洲成av人片在线播放无| 最近中文字幕高清免费大全6| h日本视频在线播放| 亚洲四区av| 国产乱人视频| 久久久国产成人免费| 国产乱人视频| 免费看av在线观看网站| 亚洲成人中文字幕在线播放| 午夜福利成人在线免费观看| 亚洲伊人久久精品综合 | 国产精品不卡视频一区二区| 国产高清国产精品国产三级 | www.色视频.com| 国产成人精品一,二区| 国产成人精品婷婷| 精品99又大又爽又粗少妇毛片| 看非洲黑人一级黄片| 丝袜喷水一区| 国产亚洲午夜精品一区二区久久 | 欧美丝袜亚洲另类| 日本色播在线视频| 女的被弄到高潮叫床怎么办| 七月丁香在线播放| 国产精品野战在线观看| 视频中文字幕在线观看| 在线a可以看的网站| 国产精品久久久久久精品电影小说 | 老女人水多毛片| 日本wwww免费看| 国产日韩欧美在线精品| 国国产精品蜜臀av免费| 美女被艹到高潮喷水动态| 久久久精品欧美日韩精品| 少妇熟女aⅴ在线视频| 一边摸一边抽搐一进一小说| 人妻夜夜爽99麻豆av| 在线免费观看不下载黄p国产| 九九在线视频观看精品| 老司机影院成人| 亚洲精品久久久久久婷婷小说 | 精品国产一区二区三区久久久樱花 | 国产黄片视频在线免费观看| 国产精品一区二区在线观看99 | 女人十人毛片免费观看3o分钟| 岛国毛片在线播放| 插阴视频在线观看视频| 精品国产露脸久久av麻豆 | 午夜精品在线福利| 2022亚洲国产成人精品| 亚洲国产成人一精品久久久| 日本wwww免费看| 久久久精品大字幕| 色哟哟·www| 欧美人与善性xxx| 六月丁香七月| 亚洲中文字幕一区二区三区有码在线看| 成人性生交大片免费视频hd| 国产午夜精品论理片| 搡女人真爽免费视频火全软件| 99久久成人亚洲精品观看| 免费观看a级毛片全部| 免费黄色在线免费观看| 男人舔奶头视频| 中文字幕亚洲精品专区| 91aial.com中文字幕在线观看| 国产久久久一区二区三区| 亚洲av成人av| 伊人久久精品亚洲午夜| 久久国内精品自在自线图片| 亚洲精品自拍成人| 麻豆精品久久久久久蜜桃| 亚洲精品乱码久久久久久按摩| 国产69精品久久久久777片| 成人国产麻豆网| 夜夜爽夜夜爽视频| 欧美性感艳星| .国产精品久久| 在线免费十八禁| 91av网一区二区| 最近最新中文字幕大全电影3| 久久久精品94久久精品| 成人高潮视频无遮挡免费网站| 国产在视频线精品| 午夜日本视频在线| 如何舔出高潮| 国产一区亚洲一区在线观看| 中文字幕av在线有码专区| 99国产精品一区二区蜜桃av| 午夜福利在线观看免费完整高清在| 乱系列少妇在线播放| 婷婷色综合大香蕉| 成人美女网站在线观看视频|