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

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

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

    周瑞靜 宋煒 倪寶峰

    摘 要:選取通州區(qū)納污河流運(yùn)潮減河和通州區(qū)區(qū)級(jí)水源地,運(yùn)用地下水水流模型Visual-Modflow和多組分溶質(zhì)運(yùn)移模型MT3DMSD耦合,以運(yùn)潮減河中濃度較高的耗氧量和氨氮為特征污染物,模擬了當(dāng)耗氧量和氨氮的初始濃度分別為20 mg/L、10 mg/L時(shí),5年后耗氧量和氨氮在左岸潛水含水層的最大超標(biāo)范圍為153 m和176 m,對(duì)承壓水沒(méi)有影響。但在運(yùn)潮減河保持現(xiàn)狀水質(zhì)情況下,潛水污染面積會(huì)逐年增大并持續(xù)向下滲濾,對(duì)承壓水造成威脅,進(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)、生活主要供水水源的國(guó)際化大都市,2015年北京市將通州區(qū)定位為城市副中心,成為北京城市行政功能戰(zhàn)略轉(zhuǎn)移區(qū)域,肩負(fù)著承接中心城區(qū)的行政功能,并以此帶動(dòng)其他功能疏解的使命。在為區(qū)域發(fā)展帶來(lái)機(jī)遇的同時(shí),又對(duì)其可持續(xù)發(fā)展提出了更高的要求,也給該區(qū)地下水環(huán)境保護(hù)與管理帶來(lái)了更大的挑戰(zhàn)。目前,城市副中心區(qū)域第一含水層組地下水多個(gè)指標(biāo)出現(xiàn)超標(biāo),以無(wú)機(jī)指標(biāo)超標(biāo)為主,局部區(qū)域甚至呈現(xiàn)無(wú)機(jī)與有機(jī)組兩類污染雙重污染的狀況,地下水環(huán)境形勢(shì)嚴(yán)峻(郭高軒等,2014)。通州區(qū)地勢(shì)低凹,自古有“九河下梢”之稱,區(qū)內(nèi)地表水系水質(zhì)較差,多為劣Ⅴ類水質(zhì)(北京市水務(wù)局,2018)。運(yùn)潮減河作為通州區(qū)的納污河流,多年來(lái)河水來(lái)源主要是雨水直排、污水處理廠再生水排放等,水質(zhì)監(jiān)測(cè)結(jié)果多年為Ⅴ類,且河底沒(méi)有防滲襯層。

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

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

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

    1 模擬區(qū)信息

    1.1 水源地概況

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

    1.2 第四系含水層特征

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

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

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

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

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

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

    2 研究方法

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

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

    2.2 邊界條件概化

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

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

    2.3 水流模型建立

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

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

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

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

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

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

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

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

    2.5 結(jié)果與分析

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

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

    3 結(jié)論

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

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

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

    參考文獻(xiàn):

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    猜你喜歡
    水質(zhì)
    水質(zhì)抽檢豈容造假
    環(huán)境(2023年5期)2023-06-30 01:20:01
    水質(zhì)檢測(cè)員——中華秋沙鴨
    水質(zhì)凈化廠提標(biāo)至一級(jí)A設(shè)計(jì)與運(yùn)行效果探討
    關(guān)于水質(zhì)監(jiān)測(cè)對(duì)環(huán)境保護(hù)的意義
    一月冬棚養(yǎng)蝦常見(jiàn)水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
    這條魚供不應(yīng)求!蝦蟹養(yǎng)殖戶、垂釣者的最愛(ài),不用投喂,還能凈化水質(zhì)
    圖像識(shí)別在水質(zhì)檢測(cè)中的應(yīng)用
    電子制作(2018年14期)2018-08-21 01:38:16
    淺析黑臭水體成因、治理方法及水質(zhì)長(zhǎng)效改善保持問(wèn)題——水質(zhì)長(zhǎng)效改善保持問(wèn)題
    濟(jì)下水庫(kù)徑流水質(zhì)和垂向水質(zhì)分析及評(píng)價(jià)
    水質(zhì)的年輪——讀《時(shí)光的年輪》
    午夜激情av网站| 亚洲成人手机| 久久99精品国语久久久| 欧美在线一区亚洲| a级毛片在线看网站| 中文欧美无线码| 99香蕉大伊视频| av网站在线播放免费| 少妇的丰满在线观看| 欧美日韩视频精品一区| 中文字幕人妻熟女乱码| 欧美日韩视频精品一区| 一区二区av电影网| 亚洲免费av在线视频| 大码成人一级视频| svipshipincom国产片| 青春草亚洲视频在线观看| xxxhd国产人妻xxx| 亚洲精品久久午夜乱码| 精品第一国产精品| 亚洲欧美中文字幕日韩二区| 国产一级毛片在线| 国产成人精品久久久久久| 国产精品熟女久久久久浪| 亚洲精品av麻豆狂野| 777米奇影视久久| 一级毛片 在线播放| 十分钟在线观看高清视频www| 视频在线观看一区二区三区| 国产精品无大码| 国产欧美日韩综合在线一区二区| 操美女的视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av电影在线进入| 日韩伦理黄色片| 精品少妇久久久久久888优播| 麻豆精品久久久久久蜜桃| av视频免费观看在线观看| 欧美精品一区二区大全| 国产精品一区二区在线不卡| 老汉色∧v一级毛片| 尾随美女入室| av福利片在线| 精品国产国语对白av| 韩国精品一区二区三区| 国产高清不卡午夜福利| 亚洲av福利一区| 在现免费观看毛片| 亚洲第一av免费看| 纵有疾风起免费观看全集完整版| 又粗又硬又长又爽又黄的视频| 免费在线观看视频国产中文字幕亚洲 | 一二三四中文在线观看免费高清| 丝瓜视频免费看黄片| 亚洲中文av在线| 精品国产乱码久久久久久小说| 日本爱情动作片www.在线观看| 亚洲精品在线美女| 女人久久www免费人成看片| 一区在线观看完整版| 综合色丁香网| 三上悠亚av全集在线观看| 国产成人精品福利久久| 一级毛片 在线播放| 婷婷色综合www| 亚洲国产精品999| 中文天堂在线官网| 伦理电影免费视频| 国产精品久久久久久久久免| 亚洲精品美女久久av网站| 国产熟女午夜一区二区三区| 好男人视频免费观看在线| 亚洲美女搞黄在线观看| 天堂中文最新版在线下载| 精品午夜福利在线看| 国产精品嫩草影院av在线观看| 精品国产一区二区三区四区第35| 亚洲精品日韩在线中文字幕| 久久热在线av| 乱人伦中国视频| 夫妻性生交免费视频一级片| 2018国产大陆天天弄谢| 久久精品久久久久久久性| 欧美黑人欧美精品刺激| 一区二区av电影网| 欧美乱码精品一区二区三区| 精品一品国产午夜福利视频| 国产极品天堂在线| 国产av国产精品国产| 久久人人97超碰香蕉20202| 成人亚洲精品一区在线观看| 精品午夜福利在线看| kizo精华| 国产黄频视频在线观看| 美女福利国产在线| 国产xxxxx性猛交| 欧美日韩一区二区视频在线观看视频在线| 看非洲黑人一级黄片| 欧美日韩视频高清一区二区三区二| 在线 av 中文字幕| 国产毛片在线视频| 99热网站在线观看| 中文字幕高清在线视频| 一区二区三区激情视频| 欧美黑人欧美精品刺激| 国产精品久久久久久人妻精品电影 | 妹子高潮喷水视频| 国产午夜精品一二区理论片| h视频一区二区三区| 亚洲精品日本国产第一区| 999精品在线视频| 波多野结衣av一区二区av| 国产成人一区二区在线| 天天添夜夜摸| 一级片免费观看大全| 亚洲美女视频黄频| 欧美日韩亚洲国产一区二区在线观看 | 久热爱精品视频在线9| 999久久久国产精品视频| 成人午夜精彩视频在线观看| 国产高清不卡午夜福利| 哪个播放器可以免费观看大片| 夫妻午夜视频| 一本久久精品| 在线观看人妻少妇| 精品国产露脸久久av麻豆| 街头女战士在线观看网站| 亚洲精品av麻豆狂野| 精品人妻一区二区三区麻豆| 亚洲精品久久久久久婷婷小说| 黄片播放在线免费| 丝袜喷水一区| 两性夫妻黄色片| 美女福利国产在线| 男女无遮挡免费网站观看| 亚洲一区二区三区欧美精品| av卡一久久| av在线app专区| 精品福利永久在线观看| 丁香六月欧美| 午夜福利乱码中文字幕| 亚洲欧美成人综合另类久久久| 又黄又粗又硬又大视频| 一本一本久久a久久精品综合妖精| 99香蕉大伊视频| 最近手机中文字幕大全| 国产精品二区激情视频| 成人国产av品久久久| 男女下面插进去视频免费观看| 国产视频首页在线观看| 十八禁高潮呻吟视频| 人成视频在线观看免费观看| 丁香六月欧美| 考比视频在线观看| 国产欧美日韩综合在线一区二区| 最近2019中文字幕mv第一页| 午夜福利在线免费观看网站| av一本久久久久| 久久久久久人人人人人| 电影成人av| 国产99久久九九免费精品| 极品少妇高潮喷水抽搐| 在线观看免费日韩欧美大片| 高清欧美精品videossex| 国产精品久久久人人做人人爽| 哪个播放器可以免费观看大片| 国语对白做爰xxxⅹ性视频网站| 精品视频人人做人人爽| 免费观看a级毛片全部| 欧美乱码精品一区二区三区| 亚洲,欧美,日韩| 性少妇av在线| 国产精品国产三级国产专区5o| 国产精品.久久久| 色综合欧美亚洲国产小说| 天堂俺去俺来也www色官网| 性高湖久久久久久久久免费观看| 国产精品欧美亚洲77777| 亚洲国产欧美日韩在线播放| 久久久亚洲精品成人影院| 久久 成人 亚洲| 亚洲国产看品久久| 国产极品粉嫩免费观看在线| 99久国产av精品国产电影| 亚洲精品一二三| 亚洲国产av影院在线观看| 欧美日韩av久久| 久久精品国产亚洲av涩爱| 一本色道久久久久久精品综合| 日韩一区二区三区影片| 90打野战视频偷拍视频| 成人18禁高潮啪啪吃奶动态图| 九九爱精品视频在线观看| 在线看a的网站| 久久人人爽av亚洲精品天堂| 青春草视频在线免费观看| 亚洲美女黄色视频免费看| 精品久久久精品久久久| 国产成人免费观看mmmm| 欧美 亚洲 国产 日韩一| 高清视频免费观看一区二区| 国产男女内射视频| 悠悠久久av| 精品免费久久久久久久清纯 | 日韩一区二区三区影片| 亚洲第一区二区三区不卡| 国产男人的电影天堂91| 婷婷色av中文字幕| 久久人妻熟女aⅴ| 伦理电影免费视频| 国产在线一区二区三区精| 夫妻性生交免费视频一级片| 中国三级夫妇交换| 观看美女的网站| 汤姆久久久久久久影院中文字幕| 爱豆传媒免费全集在线观看| 丰满乱子伦码专区| 亚洲国产精品一区三区| 侵犯人妻中文字幕一二三四区| 精品一区二区三卡| 国产熟女欧美一区二区| 亚洲一区二区三区欧美精品| 久久久久视频综合| a 毛片基地| 一区二区日韩欧美中文字幕| 高清黄色对白视频在线免费看| 国产精品亚洲av一区麻豆 | 日韩一区二区三区影片| 搡老岳熟女国产| 精品酒店卫生间| 日韩制服骚丝袜av| 欧美人与性动交α欧美软件| 亚洲第一av免费看| 永久免费av网站大全| 2018国产大陆天天弄谢| 国产精品99久久99久久久不卡 | 国产成人免费观看mmmm| 超碰成人久久| 久久精品熟女亚洲av麻豆精品| 欧美精品一区二区大全| 夫妻午夜视频| 久久人人爽人人片av| 亚洲男人天堂网一区| av福利片在线| 国产免费一区二区三区四区乱码| 精品亚洲乱码少妇综合久久| 老司机影院毛片| 亚洲av欧美aⅴ国产| 女人高潮潮喷娇喘18禁视频| 日韩熟女老妇一区二区性免费视频| 亚洲av成人精品一二三区| 日本一区二区免费在线视频| 亚洲情色 制服丝袜| 人人妻,人人澡人人爽秒播 | 午夜福利影视在线免费观看| 精品亚洲乱码少妇综合久久| 香蕉丝袜av| xxx大片免费视频| 下体分泌物呈黄色| 秋霞伦理黄片| 麻豆av在线久日| 精品一区在线观看国产| 亚洲三区欧美一区| 亚洲成av片中文字幕在线观看| 欧美日韩视频精品一区| 亚洲精品av麻豆狂野| 丁香六月天网| a级毛片在线看网站| 久久影院123| 亚洲国产成人一精品久久久| 日韩精品免费视频一区二区三区| 99久久精品国产亚洲精品| 蜜桃国产av成人99| 亚洲av国产av综合av卡| 欧美日韩av久久| 国产熟女欧美一区二区| 亚洲成人手机| 久久久国产精品麻豆| 久久久久人妻精品一区果冻| 久久久欧美国产精品| 色综合欧美亚洲国产小说| 亚洲图色成人| 久久久久精品久久久久真实原创| 丝袜在线中文字幕| 国产一区二区三区av在线| 午夜激情av网站| 男女之事视频高清在线观看 | 2021少妇久久久久久久久久久| 大香蕉久久网| xxxhd国产人妻xxx| 亚洲成av片中文字幕在线观看| 欧美另类一区| 男人添女人高潮全过程视频| 国产一区亚洲一区在线观看| 欧美精品高潮呻吟av久久| 久久久久精品久久久久真实原创| 中文字幕制服av| 超色免费av| av女优亚洲男人天堂| 欧美日韩亚洲综合一区二区三区_| 一区二区三区四区激情视频| 亚洲久久久国产精品| 亚洲精品国产色婷婷电影| 亚洲欧美色中文字幕在线| 国产一区二区在线观看av| 亚洲国产欧美日韩在线播放| 日本午夜av视频| 日韩一区二区三区影片| 又大又爽又粗| 青春草国产在线视频| 一级片免费观看大全| 欧美国产精品va在线观看不卡| 免费高清在线观看日韩| 亚洲国产最新在线播放| 亚洲av日韩精品久久久久久密 | 夜夜骑夜夜射夜夜干| 看十八女毛片水多多多| 一二三四中文在线观看免费高清| 在线天堂中文资源库| 9色porny在线观看| 99re6热这里在线精品视频| 国产视频首页在线观看| av视频免费观看在线观看| 高清av免费在线| 欧美精品高潮呻吟av久久| 亚洲av男天堂| 新久久久久国产一级毛片| 人人妻人人澡人人看| 如何舔出高潮| 最黄视频免费看| 国产一级毛片在线| 少妇被粗大猛烈的视频| 熟妇人妻不卡中文字幕| 国产熟女午夜一区二区三区| 国产欧美日韩一区二区三区在线| 操出白浆在线播放| 亚洲专区中文字幕在线 | 国产精品久久久av美女十八| 成年女人毛片免费观看观看9 | 人妻一区二区av| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷成人精品国产| 在线观看免费午夜福利视频| 伊人久久大香线蕉亚洲五| 观看av在线不卡| 亚洲精品中文字幕在线视频| 99香蕉大伊视频| 国产精品久久久久久久久免| 99热国产这里只有精品6| 另类精品久久| 日本色播在线视频| 中文字幕制服av| 欧美人与性动交α欧美精品济南到| 韩国高清视频一区二区三区| 两性夫妻黄色片| 精品国产超薄肉色丝袜足j| 国产精品偷伦视频观看了| 国产精品国产三级国产专区5o| 亚洲精品久久久久久婷婷小说| 亚洲五月色婷婷综合| 1024视频免费在线观看| 国产亚洲一区二区精品| 亚洲美女黄色视频免费看| 精品国产超薄肉色丝袜足j| 大香蕉久久成人网| 欧美 日韩 精品 国产| 观看av在线不卡| 女人被躁到高潮嗷嗷叫费观| 久久久精品国产亚洲av高清涩受| 十八禁高潮呻吟视频| 日日撸夜夜添| 国产精品麻豆人妻色哟哟久久| 国产成人欧美在线观看 | 亚洲av日韩在线播放| 欧美国产精品va在线观看不卡| 国产免费福利视频在线观看| 男人爽女人下面视频在线观看| 国产女主播在线喷水免费视频网站| 中文乱码字字幕精品一区二区三区| 国产熟女欧美一区二区| 不卡av一区二区三区| 99re6热这里在线精品视频| 9色porny在线观看| 不卡视频在线观看欧美| 啦啦啦在线免费观看视频4| 在线精品无人区一区二区三| 亚洲 欧美一区二区三区| 不卡视频在线观看欧美| 国产99久久九九免费精品| 岛国毛片在线播放| 色播在线永久视频| 欧美xxⅹ黑人| 最近的中文字幕免费完整| 日韩大码丰满熟妇| 亚洲av综合色区一区| 制服诱惑二区| 自线自在国产av| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 久久这里只有精品19| 男女午夜视频在线观看| 国产精品无大码| 2021少妇久久久久久久久久久| 久久精品熟女亚洲av麻豆精品| 国产精品av久久久久免费| 夫妻午夜视频| 精品一品国产午夜福利视频| 久久午夜综合久久蜜桃| 国产在线视频一区二区| 卡戴珊不雅视频在线播放| 亚洲人成77777在线视频| 人妻人人澡人人爽人人| 男的添女的下面高潮视频| 97精品久久久久久久久久精品| 久久99精品国语久久久| 狠狠精品人妻久久久久久综合| 精品一区二区免费观看| 国产xxxxx性猛交| 国产精品 欧美亚洲| 国产不卡av网站在线观看| 亚洲精品国产色婷婷电影| 久久亚洲国产成人精品v| 一级毛片电影观看| 男女床上黄色一级片免费看| 色播在线永久视频| 天堂中文最新版在线下载| 午夜激情av网站| 在线观看免费日韩欧美大片| 如何舔出高潮| 中国三级夫妇交换| 国产av一区二区精品久久| 亚洲精品国产一区二区精华液| 国产精品久久久av美女十八| 色婷婷av一区二区三区视频| 亚洲熟女毛片儿| 叶爱在线成人免费视频播放| 日韩中文字幕视频在线看片| 国产亚洲精品第一综合不卡| 久久久精品国产亚洲av高清涩受| 国产成人欧美| 日韩av不卡免费在线播放| 欧美国产精品va在线观看不卡| 国产精品.久久久| 丰满迷人的少妇在线观看| 两个人免费观看高清视频| 久久影院123| 国产在线免费精品| 美国免费a级毛片| 伊人亚洲综合成人网| 青春草国产在线视频| av福利片在线| 搡老岳熟女国产| 国产亚洲最大av| 久久精品国产亚洲av涩爱| 色播在线永久视频| 亚洲av欧美aⅴ国产| 免费黄色在线免费观看| 久久精品国产亚洲av高清一级| 中文乱码字字幕精品一区二区三区| 日韩人妻精品一区2区三区| 日韩 亚洲 欧美在线| 日韩av在线免费看完整版不卡| 一边摸一边做爽爽视频免费| 国产精品一国产av| 精品亚洲成国产av| 国产精品 欧美亚洲| 精品人妻熟女毛片av久久网站| av国产久精品久网站免费入址| 亚洲国产精品一区二区三区在线| 性色av一级| 午夜老司机福利片| 男女无遮挡免费网站观看| 亚洲情色 制服丝袜| 人人妻人人澡人人爽人人夜夜| 色婷婷久久久亚洲欧美| 91aial.com中文字幕在线观看| 亚洲精品一区蜜桃| 国产成人啪精品午夜网站| 亚洲久久久国产精品| 99九九在线精品视频| 亚洲av欧美aⅴ国产| 亚洲精品中文字幕在线视频| 亚洲一级一片aⅴ在线观看| 另类亚洲欧美激情| 看免费av毛片| 欧美最新免费一区二区三区| 五月天丁香电影| 街头女战士在线观看网站| www.自偷自拍.com| 夫妻午夜视频| 伦理电影大哥的女人| www日本在线高清视频| 91国产中文字幕| 精品人妻一区二区三区麻豆| 欧美精品av麻豆av| 成人亚洲精品一区在线观看| 黄色视频在线播放观看不卡| 欧美人与性动交α欧美精品济南到| 波多野结衣av一区二区av| 久久久久人妻精品一区果冻| 国产高清不卡午夜福利| 一级黄片播放器| 日本wwww免费看| 秋霞伦理黄片| 毛片一级片免费看久久久久| 久久久久久久久久久久大奶| 狠狠精品人妻久久久久久综合| 成人黄色视频免费在线看| 在线观看人妻少妇| 国产精品 欧美亚洲| 国产在线免费精品| 亚洲欧美中文字幕日韩二区| av免费观看日本| 一个人免费看片子| 亚洲情色 制服丝袜| 国产有黄有色有爽视频| 一区二区三区激情视频| 日日撸夜夜添| 久久av网站| 日韩免费高清中文字幕av| 超色免费av| 精品亚洲成a人片在线观看| 搡老乐熟女国产| 色网站视频免费| 国产亚洲av高清不卡| 精品国产乱码久久久久久小说| 国产色婷婷99| 亚洲色图综合在线观看| 免费观看av网站的网址| 亚洲三区欧美一区| 日日撸夜夜添| 深夜精品福利| a级毛片在线看网站| 亚洲专区中文字幕在线 | 午夜免费观看性视频| 在线观看免费高清a一片| 电影成人av| 亚洲图色成人| 搡老乐熟女国产| 美女脱内裤让男人舔精品视频| 久久人人爽人人片av| 99久久精品国产亚洲精品| 国产不卡av网站在线观看| 精品久久久精品久久久| 国产精品国产三级专区第一集| 涩涩av久久男人的天堂| 天堂8中文在线网| 国产有黄有色有爽视频| 亚洲精品中文字幕在线视频| 又粗又硬又长又爽又黄的视频| 深夜精品福利| 久久久久网色| 欧美日韩一级在线毛片| 巨乳人妻的诱惑在线观看| 久久97久久精品| 制服丝袜香蕉在线| 久久久久久久大尺度免费视频| 七月丁香在线播放| 人妻一区二区av| 亚洲四区av| 在线观看免费午夜福利视频| 国产色婷婷99| 日本爱情动作片www.在线观看| 99精国产麻豆久久婷婷| 老司机影院成人| 电影成人av| 97在线人人人人妻| 国产成人av激情在线播放| 久久久久久久久久久免费av| 久久人人爽人人片av| 国产精品亚洲av一区麻豆 | 国语对白做爰xxxⅹ性视频网站| 亚洲人成电影观看| 亚洲第一av免费看| 91精品国产国语对白视频| 亚洲av国产av综合av卡| 新久久久久国产一级毛片| 欧美日韩视频高清一区二区三区二| 香蕉丝袜av| 亚洲av综合色区一区| 免费久久久久久久精品成人欧美视频| videosex国产| 日韩中文字幕欧美一区二区 | 欧美日韩视频高清一区二区三区二| 老汉色av国产亚洲站长工具| 水蜜桃什么品种好| 日本vs欧美在线观看视频| 99国产精品免费福利视频| 黄色视频在线播放观看不卡| 亚洲国产av影院在线观看| 尾随美女入室| 国产精品熟女久久久久浪| 人人妻人人澡人人看| 亚洲视频免费观看视频| 日本欧美国产在线视频| 成人国产麻豆网| 亚洲视频免费观看视频| 久久精品亚洲熟妇少妇任你| 欧美日韩福利视频一区二区| 蜜桃在线观看..| 一区二区三区乱码不卡18| 别揉我奶头~嗯~啊~动态视频 | 日韩 欧美 亚洲 中文字幕| 91精品国产国语对白视频| 国产福利在线免费观看视频| 亚洲精品一区蜜桃| 校园人妻丝袜中文字幕| 欧美精品人与动牲交sv欧美| 国产亚洲av片在线观看秒播厂| 久久久久久久精品精品| 久久天躁狠狠躁夜夜2o2o | 国产欧美亚洲国产| 国产av国产精品国产| 国产精品av久久久久免费| 国产片特级美女逼逼视频|