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

    西北內(nèi)陸張掖盆地地下水溫度變化特征及其指示意義

    2011-04-19 12:36:52連英立張光輝聶振龍
    地球?qū)W報(bào) 2011年2期
    關(guān)鍵詞:承壓水潛水盆地

    連英立, 張光輝, 聶振龍, 艾 婷

    1)中國地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所, 河北石家莊 050061;

    2)河北工程大學(xué)經(jīng)濟(jì)管理學(xué)院, 河北邯鄲 056038

    西北內(nèi)陸張掖盆地地下水溫度變化特征及其指示意義

    連英立1), 張光輝1), 聶振龍1), 艾 婷2)

    1)中國地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所, 河北石家莊 050061;

    2)河北工程大學(xué)經(jīng)濟(jì)管理學(xué)院, 河北邯鄲 056038

    通過西北內(nèi)陸張掖盆地地下水溫度變化特征研究表明: 最近 20年(1990—2009)研究區(qū)內(nèi)潛水和承壓水溫度分別下降了0.78℃和1.17℃, 與氣溫升高導(dǎo)致山區(qū)冰雪融水和降水量增加有關(guān); 同時(shí), 人類開采取用地下水強(qiáng)度增大使得研究區(qū)內(nèi)潛水和承壓水溫度變化趨勢近同。盆地平原區(qū)地下水溫度變化與山區(qū)氣溫和降水量變化成反比關(guān)系, 即山區(qū)氣溫升高、降水量增加, 盆地地下水溫度降低; 山區(qū)氣溫降低、降水量減少, 盆地地下水溫度升高。由于盆地平原區(qū)地下水補(bǔ)給水源的70%~90%來自上游祁連山區(qū)降水和冰雪融水出山地表徑流, 水溫較低, 因此, 盆地平原區(qū)地下水溫度變化程度反映著其獲取有效補(bǔ)給水量多少, 具有一定的指示意義。

    張掖盆地; 地下水; 溫度; 影響因素; 指示意義

    Key words:Zhangye basin; groundwater; temperature; influencing factors; indicative significance

    地下水溫度是地下水動(dòng)態(tài)研究的一項(xiàng)重要內(nèi)容。地下水溫度的變化可以反映地下水補(bǔ)給水源組成、徑流和排泄條件的變化, 并從一定程度上反映氣候條件變化和人類活動(dòng)對地下水的影響。全球氣候變化國際間政府合作組織(IPCC)于2007年11月發(fā)布的第四次全球氣候變化報(bào)告指出, 最近 100年(1906—2005年)中, 全球氣溫普遍升高, 且在北半球高緯度地區(qū)氣溫升幅較大, 全球氣溫升高線性趨勢為0.74℃(IPCC, 2007)。根據(jù)中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)提供的數(shù)據(jù)計(jì)算得出, 張掖盆地最近20年(1990—2009)中的氣溫升高線性趨勢為 1.14℃, 其上游山區(qū)氣溫升高的線性趨勢為 1.27℃, 均大于全球平均值。在此背景下地下水溫度如何變化及其指示意義何在成為眾多學(xué)者關(guān)注的科學(xué)問題, 國內(nèi)外的有關(guān)學(xué)者就此開展了相關(guān)研究(Gosnold et al., 1997; Bodri et al., 2001; 方燕娜等, 2005; Miyakoshi et al., 2005; 方燕娜等, 2006; 林學(xué)鈺等, 2009; Taylor et al., 2009), 普遍認(rèn)為氣候變暖是導(dǎo)致地下水溫度升高的主要原因, 地下水溫度升高是地下水對全球氣候變暖的一種響應(yīng)。但這一研究結(jié)果是否適用于所有地區(qū)有待考究, 特別是在我國西北內(nèi)陸盆地,地下水主要來源于出山河水入滲補(bǔ)給, 且上游補(bǔ)給源區(qū)存在冰川, 其溫度變化規(guī)律及其指示意義如何尚不清楚, 仍需做大量深入的研究。

    本研究以河西走廊張掖盆地為例, 以多年地下水溫度動(dòng)態(tài)觀測數(shù)據(jù)和氣候變化觀測資料為基礎(chǔ),分析了近20年來的地下水溫度變化規(guī)律, 揭示了不同層位地下水溫度與各氣候要素之間的內(nèi)在聯(lián)系,并探討了影響其變化的主導(dǎo)因素, 進(jìn)而闡明地下水溫度變化的指示意義。

    1 研究區(qū)概況

    研究區(qū)張掖盆地位于河西走廊中部黑河流域中游, 東起民樂永固隆起, 西止高臺——南華榆木山構(gòu)造隆起地下水匯水線, 南、北界于祁連山與龍首山、合黎山之間, 是一個(gè)沿SE-NW向展布的半封閉的大型山間斷陷盆地, 面積5062 km2, 年降水量120 mm, 年水面蒸發(fā)量大于 2050 mm, 年平均氣溫 7.5℃, 屬于溫帶干旱氣候(張光輝等, 2005b; 賈艷琨等, 2008)。研究區(qū)內(nèi)的主要河流有黑河, 梨園河和山丹河(圖1)。

    圖1 研究區(qū)位置圖Fig.1 Location of the study area

    由于降水量小蒸發(fā)量大, 使得盆地自產(chǎn)水資源能力不足 10%, 盆地地下水主要依賴祁連山區(qū)出山口地表徑流補(bǔ)給, 包括冰川雪融水、降水和基巖裂隙水(基流)補(bǔ)給, 他們與祁連山區(qū)的氣候變化密切相關(guān)(張光輝等, 2006, 2005c)。研究區(qū)內(nèi)地下水按埋藏特征劃分, 南部為單層結(jié)構(gòu)潛水系統(tǒng)和北部為多層結(jié)構(gòu)潛水—承壓水系統(tǒng)(圖2)。由祁連山山緣至洪積扇扇緣的溢出帶, 為單層結(jié)構(gòu)的潛水系統(tǒng)。山前洪積扇頂部地下水水位埋深大于 200 m, 含水層由粗顆粒的砂礫卵石組成。至扇中地帶, 地下水水位埋深150~50 m, 含水層中含泥質(zhì)漸多。在扇緣地帶,含水層顆粒漸細(xì), 地下水水位埋深漸淺。至張掖—臨澤一帶, 地下水以泉的形式溢出, 含水層結(jié)構(gòu)由單一的潛水系統(tǒng)漸變?yōu)槎鄬拥臐撍袎核到y(tǒng)。在黑河—山丹河沿岸地帶, 地下水水位埋深小于 3 m, 溝壑和洼地內(nèi)有成片泉水出露。在山前地帶, 接受出山河水入滲補(bǔ)給, 地下水礦化度<1g/L, 屬于HCO3型水。

    在溢出帶以下的細(xì)土平原, 為多層結(jié)構(gòu)的潛水—承壓水系統(tǒng), 含水層以亞砂土、亞粘土和砂礫石互層為主, 含水層單層厚度 20~30 m, 第一承壓含水層頂板埋深10 m以上, 承壓水頭一般高于潛水水位 1~2 m, 在臨澤的農(nóng)場—小屯一帶為自流區(qū), 地下水水位高出地表0.5~3.0 m。潛水的礦化度>1 g/L,為 SO4-HCO3或 SO4型水, 承壓水的礦化度<1 g/L,為HCO3型水。

    圖2 張掖盆地水文地質(zhì)剖面圖Fig.2 Hydrogeological section of Zhangye basin

    2 地下水溫度變化特征

    此次研究選取了研究區(qū)內(nèi) 1990—2009年間近100眼觀測井的地下水溫觀測資料(個(gè)別井和年份溫度資料缺失), 和黑河上游山區(qū)和中游盆地的氣溫、降水資料, 取各觀測資料的多年(1990—2009)月平均值做出了研究區(qū)內(nèi)地下水溫度、氣溫、降水的年內(nèi)變化曲線(圖3), 同時(shí)取其年平均值做出了研究區(qū)內(nèi)氣溫、降水量、水位埋深和地下水溫度年際變化曲線(圖4)。

    2.1 地下水溫度的年內(nèi)變化特征

    研究區(qū)內(nèi)潛水年平均溫度為 10.9℃, 承壓水年平均溫度為10.7℃, 從(圖3)可以看出潛水和承壓水的年內(nèi)最低溫均出現(xiàn)在 3月份, 最低溫度分別為9.5℃和9.7℃, 盆地潛水的年內(nèi)最高溫度出現(xiàn)在9月份為12.5℃, 承壓水的年內(nèi)最高溫度出現(xiàn)在10月份為 11.8℃。從變化趨勢來看, 研究區(qū)內(nèi)潛水和承壓水年內(nèi)均呈現(xiàn)下降—上升—下降的變化趨勢, 但潛水溫度的年內(nèi)變化較承壓水劇烈, 其年內(nèi)最低和最高溫度均小于和大于承壓水, 其年內(nèi)變化幅度, 潛水約為3℃, 承壓水約為2℃。

    從(圖3)可以看出, 山區(qū)和盆地的氣溫、降水年內(nèi)變化存在較高的一致性, 均呈現(xiàn)上升—下降的變化趨勢, 氣溫和降水的年內(nèi)最大值在山區(qū)和盆地均出現(xiàn)在每年的 7月份, 山區(qū)和盆地的最低氣溫均出現(xiàn)在每年的 1月份, 降水量最小值山區(qū)出現(xiàn)在每年的12月份, 盆地出現(xiàn)在每年的2月份。研究區(qū)內(nèi)潛水溫度最大值分別滯后于山區(qū)和盆地的氣溫、降水2個(gè)月出現(xiàn); 承壓水溫度最大值分別滯后于山區(qū)和盆地的氣溫、降水 3個(gè)月出現(xiàn)。如果從多年連續(xù)變化來看, 每年的3—9月份是潛水溫度的上升期, 9月份至次年的3月份為潛水溫度的下降期, 其中, 5—8月份上升幅度最大, 10—12月份下降幅度最大; 每年的3—10月份是承壓水溫度的上升期, 10月份至次年的3月份為承壓水溫度的下降期, 其中4—7月份上升幅度最大, 10—12月份下降幅度最大。

    圖3 研究區(qū)氣溫與地下水溫度、降水量與地下水溫度年內(nèi)變化關(guān)系圖Fig.3 Variation curve of atmospheric temperature and groundwater temperature(a), precipitation and groundwater temperature(b) within a year in the study area

    圖4 研究區(qū)氣溫、降水、與地下水位、溫度年際變化曲線Fig.4 Variation curve of temperature, precipitation, water table and groundwater temperature between various years in the study area

    2.2 地下水溫度的年際變化特征

    相比年內(nèi)變化, 人們更關(guān)心地下水長時(shí)間的年際變化。從(表1)和(圖4)可以看出, 最近20年(1990—2009)中, 研究區(qū)內(nèi)潛水和承壓水的年平均溫度均呈現(xiàn)下降的趨勢, 分別下降了 0.78℃和 1.17℃, 下降速率分別為?0.04℃/a和?0.06℃/a, 承壓水溫度的下降速率是潛水的 1.5倍。由(圖 4)可以看出, 1990—1995年間研究區(qū)潛水和承壓水溫度呈現(xiàn)互逆變化,其分別經(jīng)歷了降低—升高、升高—降低的過程, 且在此過程中, 潛水和承壓水分別于 1992年和 1993年達(dá)到最大值和最小值, 期間承壓水溫度高于潛水溫度, 最大差值約為0.4℃。1995—1998年間潛水與承壓水溫度變化趨勢相同, 只有1996年承壓水溫度高于潛水溫度, 其余年份承壓水溫度均低于潛水溫度。自1998年之后研究區(qū)內(nèi)潛水和承壓水均呈現(xiàn)顯著下降趨勢, 其中2004年出現(xiàn)較大幅度回升, 但總體上呈現(xiàn)波動(dòng)下降的趨勢。在此過程中(1998—2009),潛水和承壓水的溫度波動(dòng)變化趨勢也趨于相同, 尤其是自2001年之后, 波動(dòng)變化趨勢完全一致, 且承壓水溫度低于潛水溫度, 二者溫度差值與波動(dòng)幅度大小成正比, 且自 2001至 2009年潛水、承壓水年平均溫度差值呈變小的趨勢, 到2009年二者差值幾乎為 0。在1998—2009年間分別于1998—2001年和2004—2006年兩次出現(xiàn)較大幅度下降, 且持續(xù)時(shí)間均為 2~3年。同時(shí), 在 2002年地下水溫度出現(xiàn)一次較大幅度回升。

    3 地下水溫度變化影響因素分析

    影響地下水溫度變化的因素主要有: 熱源(太陽輻射和地球內(nèi)部熱流)、地質(zhì)條件(包括地層巖性、構(gòu)造、地形等)、補(bǔ)給水源特征(包括河水補(bǔ)給和當(dāng)?shù)亟邓a(bǔ)給的溫度等特征)、當(dāng)?shù)貧夂蚝腿祟惢顒?dòng)影響等(方燕娜等, 2005; 方燕娜等, 2006; 林學(xué)鈺等, 2009)。對于同一地區(qū), 地質(zhì)條件短時(shí)期內(nèi)不會發(fā)生大的變化, 其對地下水溫度短期內(nèi)變化影響不大,盆地內(nèi)潛水平均水位埋深在 8~10 m, 太陽輻射和地球內(nèi)部熱流對其年均溫度變化影響微乎其微(巴建文等, 2009), 因此研究區(qū)地下水溫度變化主要受補(bǔ)給水源溫度以及當(dāng)?shù)貧夂蚝腿祟惢顒?dòng)的影響。

    表1 研究區(qū)氣溫、降水與地下水溫度變化表Table 1 Temperature, precipitation and groundwater temperature variation in the study area

    由(圖5)可以看出, 盆地地下水溫度的分布規(guī)律是由山前向盆地水溫逐漸升高, 由東南向西北水溫逐漸升高, 即與地下水流場方向一致; 同時(shí)與盆地內(nèi)距河道的距離成正比, 距離河道越近溫度越低,距離河道越遠(yuǎn), 溫度越高, 這也證明了盆地地下水主要補(bǔ)給來源為出山徑流且補(bǔ)給水源溫度低于盆地地下水溫度, 出山河水經(jīng)過長途運(yùn)移后吸收盆地內(nèi)的熱能, 溫度逐漸升高。地下水流場中水溫的變化屬于地下水熱量運(yùn)移范疇。地下水熱量運(yùn)移是一個(gè)非常復(fù)雜的過程, 其不僅與地下水及其所在的介質(zhì)有關(guān), 而且與地下水的運(yùn)動(dòng)狀態(tài)有密切的關(guān)系(王錦國等, 2001), 在水溫和流經(jīng)介質(zhì)相同的情況下, 流速越快吸熱越少, 溫度變化越小; 流速越慢吸熱越多, 水溫變化越大。(注: 其中個(gè)別不符合規(guī)律的高溫點(diǎn)大都位于泉水溢出帶, 地下水流速較慢, 受盆地輻射和氣候變化影響較大。)

    導(dǎo)致最近 20年來(1990—2009), 研究區(qū)內(nèi)潛水和承壓水溫度分別下降0.78℃和1.17℃的原因主要有以下3個(gè)方面:

    (1)研究區(qū)地下水補(bǔ)給量增大

    圖5 研究區(qū)內(nèi)采樣點(diǎn)潛水溫度分布圖(2007-9)Fig.5 Distribution of shallow groundwater temperatures of sampling points in the study area in September, 2007

    普遍認(rèn)為, 盆地 70%~90%的地下水補(bǔ)給來源于山區(qū)(張光輝等, 2005a), 山區(qū)降水量的增加直接導(dǎo)致出山徑流量增加, 進(jìn)而使得盆地地下水補(bǔ)給量的增加。由(表 1)可以看出最近 20年(1990—2009),山區(qū)和盆地的降水量分別增加了 1/4和 1/10, 分別為 95.25 mm 和 12.28 mm, 增加速率分別為 5.01 mm/a和0.65 mm/a。平均山區(qū)每增加1 mm降水, 盆地地下水(含潛水和承壓水)溫度下降0.01℃。研究區(qū)內(nèi)黑河與梨園河多年平均(1956—2009)流量為18.44億m/a, 2003—2009年平均流量為21.45億m/a, 為近60年來的最豐時(shí)期(圖6)。

    (2)補(bǔ)給水源組成發(fā)生變化, 補(bǔ)給水源溫度降低

    在研究區(qū)的上游山區(qū)存在大量冰川, 且普遍很小, 最大的也只有 2.48 km2, 絕大部分的冰川面積<1 km2, 所以它們對溫度的反應(yīng)很敏感(聶振龍等, 2001; 陽勇等, 2007)。由(表 1)可以看出最近 20年(1990—2009), 山區(qū)和盆地的氣溫分別升高了 1.27℃和1.14℃, 升高速率分別為0.07℃/a和0.06℃/a。根據(jù)中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)提供的數(shù)據(jù), 山區(qū)1990年4月和10月的平均溫度分別為-1.0℃和-1.8℃, 到2009年分別為3.9℃和0.8℃, 均達(dá)到0℃以上。由(圖7)可以看出山區(qū)氣溫升高速率1月、4月和 10月份最為明顯均超過了 0.1℃/a, 其中 10月份最高超過了0.12℃/a。4月和10月分別為冬春和秋冬季節(jié)的交替期, 該時(shí)段的氣溫升高直接導(dǎo)致山區(qū)融雪期延長, 融雪徑流過程提前, 融雪徑流量顯著增長, 徑流模數(shù)增加(王建等, 2005; 賀建橋等, 2008), 從而最終使得山區(qū)水源溫度趨于下降的趨勢。

    圖6 研究區(qū)內(nèi)河流出山口流量圖Fig.6 River discharge at the mountain-pass in the study area

    圖7 山區(qū)氣溫升高率(1990—2009)Fig.7 Temperature rise rate in the mountain area (1990—2009)

    (3)埋藏條件和補(bǔ)給途徑不同導(dǎo)致潛水和承壓水溫度下降速率不同

    研究區(qū)內(nèi)潛水埋藏深度較淺, 平均在 8~10 m左右, 其主要補(bǔ)給來源為引河灌溉和河流側(cè)滲, 在長距離的地表運(yùn)移過程中吸熱較多, 且存在強(qiáng)烈的蒸發(fā)影響(陳宗宇等, 2006); 而承壓水的補(bǔ)給是由山區(qū)來水在山前經(jīng)由粗顆粒砂礫石快速滲入地下(張光輝等, 2005a), 以潛流的方式補(bǔ)給, 水的運(yùn)移過程吸熱較少, 加之盆地承壓水埋深普遍大于 80 m, 因此在相同的補(bǔ)給水源條件下, 研究區(qū)內(nèi)潛水溫度下降速率比承壓水速率要小。

    1990—1995年間研究區(qū)潛水和承壓水溫度呈現(xiàn)互逆變化, 且承壓水溫度高于潛水溫度, 出現(xiàn)這種逆反趨勢的原因有以下2個(gè)方面:

    (1)潛水埋深較淺, 水溫受埋深和盆地氣候控制

    1990—1995年間統(tǒng)計(jì)的觀測井中, 75%的觀測井潛水埋深在10 m以淺, 使得研究區(qū)潛水溫度受盆地氣候的影響較大。雖然在此期間山區(qū)降水量和出山徑流量均經(jīng)歷了降低——升高的過程, 但是潛水并不是與山區(qū)降水量呈相反的響應(yīng)過程。由(圖4)可以看出盆地氣溫在 1990—1995年間經(jīng)歷了降低——升高的過程, 這一過程與研究區(qū)潛水溫度的變化趨勢相同。由(圖8)亦可以看出, 1990—1995年間, 研究區(qū)內(nèi)潛水溫度與潛水埋深之間的線性關(guān)系明顯, 即埋深增加溫度降低, 埋深減小溫度升高。因此, 在 1990—1995年間潛水溫度主要受控于埋深和盆地氣候。

    圖8 研究區(qū)內(nèi)潛水年均溫度、埋深散點(diǎn)圖Fig.8 Scatter diagram of shallow groundwater temperatures and depths

    (2)山區(qū)降水量減少是盆地承壓水溫度升高的主要因素

    由于盆地承壓水位埋深較大, 且其補(bǔ)給途徑主要由地下潛流補(bǔ)給, 因此其受盆地氣候變化的影響很小, 其溫度主要受控于山區(qū)降水量和補(bǔ)給水源溫度變化, 而補(bǔ)給水源溫度變化又與山區(qū)氣溫呈負(fù)相關(guān)。由圖4可以看出, 1990—1992年山區(qū)氣溫變化不大, 而同期山區(qū)降水量卻減少了將近100 mm, 山區(qū)降水量減少使得盆地承壓水接受補(bǔ)給量減少, 水量減少導(dǎo)致補(bǔ)給水在地下運(yùn)移過程中速度減慢吸熱較多水溫升高。因此, 山區(qū)降水量減少和運(yùn)移過程中吸熱升溫是導(dǎo)致同期盆地承壓水溫度上升的主要原因。同樣, 1992年之后山區(qū)氣溫和降水量均呈升高的趨勢, 這是導(dǎo)致同期盆地承壓水溫度下降的主要原因。

    1998年之后研究區(qū)內(nèi)潛水和承壓水溫度變化趨勢趨于相同, 尤其是自2001年之后, 變化趨勢完全一致, 這主要是因?yàn)?

    (1)2001年之后山區(qū)氣溫和降水量持續(xù)增加使得盆地地下水得到充分的補(bǔ)給, 潛水埋深增加, 受盆地氣候影響減小

    山區(qū)多年(1960—2009)平均氣溫為?0.94℃, 2001—2009年間平均氣溫為-0.09℃, 高于多年平均值0.85℃。山區(qū)多年(1960—2009)平均降雨量為409 mm/a, 其中2001—2009年平均降水量為440 mm/a,高于多年平均值 31mm, 2003—2009年更是達(dá)到了460 mm/a, 高于平均值51 mm。氣溫的升高導(dǎo)致山區(qū)融雪期延長(王建等, 2005), 冰川雪融水增加, 補(bǔ)給水源溫降低; 降水增加使得出山口河流量增大,盆地地下水得到充分補(bǔ)給。同時(shí), 潛水埋深增加(圖4)使得其受盆地氣候的影響減弱, 加之充足的補(bǔ)給最終使得潛水、承壓水溫度變化趨同。

    (2)黑河調(diào)水方案的實(shí)施使得研究區(qū)內(nèi)地下水開采量趨于增加, 原有的補(bǔ)給排泄方式被改變。

    隨著 2001年黑河向下游額濟(jì)納旗調(diào)水方案的實(shí)施, 為了保證下游額濟(jì)納旗生態(tài)用水, 使得盆地地表引水量受限, 導(dǎo)致盆地深層地下水開采量增加。由于開采量數(shù)據(jù)較難統(tǒng)計(jì), 因此沒有直接的證據(jù)證明盆地內(nèi)開采量增加, 但是, 由(圖9)可以看出,在第一次實(shí)施分水方案的2001年, 研究區(qū)內(nèi)開采井?dāng)?shù)量驟然增加了將近 900眼, 約為前 5年新增開采井?dāng)?shù)量的總和。同時(shí), 由于絕大部分開采井都是混合開采, 導(dǎo)致盆地內(nèi)潛水、承壓水水力聯(lián)系強(qiáng)度增加, 含水層之間連通性越來越好, 進(jìn)而使得各含水層的溫度變化趨于相同。

    圖9 研究區(qū)內(nèi)年新增開采井?dāng)?shù)量曲線圖Fig.9 Newly-increased mining wells every year in the study area

    眾所周知, 1998年由于厄爾尼諾—拉尼娜現(xiàn)象的出現(xiàn)導(dǎo)致全球氣候異常, 同年我國境內(nèi)降水量激增, 各地出現(xiàn)洪澇災(zāi)害, 而西部地區(qū)也不例外。由(圖4)可以看出, 在1998年和2003年山區(qū)和盆地降水量和氣溫均出現(xiàn)較大峰值, 其中1998年山區(qū)降水量約為多年平均值的 1.5倍, 短期內(nèi)大量降水注入盆地是導(dǎo)致盆地地下水溫度連續(xù) 2~3年下降的主要原因。同樣, 2002年的溫度升高與2000—2002年間持續(xù)干旱, 山區(qū)降水量減少有很大關(guān)系, 且這種影響也是在持續(xù)2~3年的干旱之后, 地下水溫度對氣候變化做出的響應(yīng)。因此, 可以說地下水對極端氣候的響應(yīng)具有2~3年的滯后。

    4 地下水溫度變化的指示意義

    上述分析表明, 盆地平原區(qū)地下水溫度與氣候變化和人類活動(dòng)的關(guān)系密切, 因此, 地下水溫度的變化可以被認(rèn)為是氣候變化和人類活動(dòng)強(qiáng)度的指示器。

    首先, 地下水溫度在空間上的分布反應(yīng)了研究區(qū)內(nèi)地下水補(bǔ)徑排條件, 即盆地內(nèi)地下水接受出山河水補(bǔ)給, 沿地下水流場和河道的水溫變化反應(yīng)了地下水補(bǔ)給過程中的徑流途徑。

    其次, 由于盆地平原區(qū)地下水補(bǔ)給水源的70%~90%來自上游祁連山區(qū)降水和冰雪融水出山地表徑流, 水溫較低, 因此, 盆地平原區(qū)地下水溫度變化程度反映著其獲取有效補(bǔ)給水量多少, 具有一定的指示意義。

    最后, 由于貫通井和開采量的增加使各含水層之間連通性越來越好, 導(dǎo)致區(qū)內(nèi)潛水和承壓水溫度變化趨于相同。因此, 研究區(qū)內(nèi)潛水和承壓水的溫度變化趨同程度可以反映人類活動(dòng)的強(qiáng)度, 具有一定指示意義。

    5 結(jié)論

    (1)研究區(qū)內(nèi)潛水和承壓水多年月平均溫度年內(nèi)變幅分別約為3℃和2℃, 最低溫度均出現(xiàn)在3月份,最高溫度分別出現(xiàn)在9月和10月; 潛水和承壓水年內(nèi)最高溫度分別滯后山區(qū)和盆地氣溫、降水最大值2個(gè)月和3個(gè)月出現(xiàn)。

    (2)最近20年(1990—2009)中, 研究區(qū)內(nèi)潛水和承壓水的年平均溫度均呈現(xiàn)下降的趨勢, 且其下降的線性趨勢分別為-0.78℃和-1.17℃, 下降速率分別為-0.04℃/a和-0.06℃/a, 承壓水溫度的下降速率是潛水的1.5倍。

    (3)山區(qū)氣溫升高降水量增加是導(dǎo)致盆地地下水溫度下降的主要因素, 且山區(qū)氣溫升高對研究區(qū)內(nèi)承壓水溫度變化的影響要大于對潛水溫度變化的影響, 降水量增加對研究區(qū)內(nèi)潛水溫度變化的影響大于對承壓水溫度變化的影響; 研究區(qū)內(nèi)河流出山口徑流增大和盆地平原區(qū)人類活動(dòng)是導(dǎo)致 2001年之后研究區(qū)內(nèi)潛水與承壓水溫度變化趨勢相同的主要原因。

    (4)研究區(qū)內(nèi)地下水溫度對極端氣候變化的響應(yīng)具有2~3年的滯后。

    (5)研究區(qū)內(nèi)地下水溫度變化程度可以反映其獲取有效補(bǔ)給水量的多少, 潛水和承壓水溫度變化的趨同程度可以反映人類活動(dòng)強(qiáng)度的強(qiáng)弱, 具有一定的指示意義。

    巴建文, 陳新民, 田遼西, 秦曉燕, 劉振華.2009.甘肅省張掖市北郊濕地環(huán)境地質(zhì)研究報(bào)告[R].甘肅張掖: 甘肅省地礦局水文地質(zhì)工程地質(zhì)勘察院.

    陳宗宇, 萬力, 聶振龍, 申建梅, 陳京生.2006.利用穩(wěn)定同位素識別黑河流域地下水的補(bǔ)給來源[J].水文地質(zhì)工程地質(zhì), (6): 9-14.

    方燕娜, 廖資生, 陳洪艷, 董維紅.2006.吉林中部平原區(qū)地下水溫動(dòng)態(tài)變化影響因素的關(guān)聯(lián)度分析[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版, 36(1): 66-72.

    方燕娜, 林學(xué)鈺, 廖資生, 陳洪艷.2005.吉林中部平原區(qū)淺層地下水溫度對全球氣候變暖的響應(yīng)[J].冰川凍土, 27(6): 833-837.

    賀建橋, 宋高舉, 蔣熹, 李全蓮, 武曉波.2008.2006年黑河水系典型流域冰川融水徑流與出山徑流的關(guān)系[J].中國沙漠, 28(6): 1186-1189.

    賈艷琨, 劉福亮, 張琳, 聶振龍, 申建梅, 陳宗宇.2008.利用環(huán)境同位素識別酒泉—張掖盆地地下水補(bǔ)給和水流系統(tǒng)[J].地球?qū)W報(bào), 29(6): 740-744.

    林學(xué)鈺, 方燕娜, 廖資生, 陳洪艷.2009.全球氣候變暖和人類活動(dòng)對地下水溫度的影響[J].北京師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 45(Z1): 452-457.

    聶振龍, 郭占榮, 焦鵬程, 董華.2001.西北內(nèi)陸盆地水循環(huán)特征分析[J].地球?qū)W報(bào), 22(4): 302-306.

    王建, 李碩.2005.氣候變化對中國內(nèi)陸干旱區(qū)山區(qū)融雪徑流的影響[J].中國科學(xué)(D輯), 35(7): 664-670.

    王錦國, 周志芳, 金忠青.2001.地下水熱量運(yùn)移模擬的BEM-FAM耦合法[J].水利學(xué)報(bào), (5): 71-76.

    陽勇, 陳仁升, 吉喜斌.2007.近幾十年來黑河野牛溝流域的冰川變化[J].冰川凍土, 29(1): 100-106.

    張光輝, 陳宗宇, 聶振龍, 劉少玉, 張翠云, 申建梅, 王金哲, 程旭學(xué), 張荷生.2006.黑河流域地下水同位素特征及其對古氣候變化的響應(yīng)[J].地球?qū)W報(bào), 27(4): 341-348.

    張光輝, 劉少玉, 謝悅波.2005a.西北內(nèi)陸黑河流域水循環(huán)與地下水形成演化模式[M].北京: 地質(zhì)出版社: 60-66.

    張光輝, 聶振龍, 劉少玉, 陳宗宇, 張翠云, 申建梅, 王金哲, 程旭學(xué), 張荷生.2005b.黑河流域走廊平原地下水補(bǔ)給源組成及其變化[J].水科學(xué)進(jìn)展, 16(5): 673-678.

    張光輝, 聶振龍, 張翠云, 申建梅.2005c.黑河流域走廊平原地下水補(bǔ)給變異特征與機(jī)制[J].水利學(xué)報(bào), 36(6): 715-720.

    References:

    BA Jian-wen, CHEN Xin-min, TIAN Liao-xi, QIN Xiao-yan, LIU Zhen-hua.2009.Northern Suburb Wetland Environmental and Geological Research Report in Zhangye City, Gansu Province[R].Gansu Zhangye: Gansu Bureau of Geology and Mineral Exploration Institute of Hydrogeology and Engineering Geology(in Chinese).

    BODRI L, CERMAK V, KUKKONEN I T.2001.Climate change of the last 2000 years inferred from borehole temperatures: data from Finland[J].Global and Planetary Change, 29(3-4): 189-200.

    CHEN Zong-yu, WAN Li, NIE Zhen-long, SHEN Jian-mei, CHEN Jing-sheng.2006.Identification of groundwater recharge in the Heihe Basin using environmental isotopes[J].Hydrogeology & Engineering Geology, (6): 9-14(in Chinese with English abstract).

    FANG Yan-na, LIAO Zi-sheng, CHEN Hong-yan, DONG Wei-hong.2006.Correlativity Analysis of the Affect Factors on Groundwater Thermal Regime in the Plain Area of Central Jilin[J].Journal of Jilin University(Earth Science Edition), 36(1): 66-72(in Chinese with English abstract).

    FANG Yan-na, LIN Xue-yu, LIAO Zi-sheng, CHEN Hong-yan.2005.Impact of Global Warming on the Groundwater Temperature of Shallow Aquifer in the Plain Area of Central Jilin Province[J].Journal of Glaciology and Geocryology, 27(6): 833-837(in Chinese with English abstract).

    GOSNOLD W D, TODHUNTER P E, SCHMIDT W.1997.The borehole temperature record of climate warming in the mid-continent of North America[J].Global and Planetary Change, 15(1-2): 33-45.

    HE Jian-qiao, SONG Gao-jü, JIANG xi, LI Quan-lian, WU Xiao-bo.2008.Relation between Glacial Meltwater Runoff and Mountainous Runoff in 2006 in Four Typical River Basins of Heihe River Water System[J].Journal of Desert Research, 28(6): 1186-1189(in Chinese with English abstract).

    IPCC.2007.Climate Change 2007: Impacts, Adaptation and Vul-nerability.Contribution of Working Group Ⅱto the Fourth Assessment Report of the Intergovernmental Panel on Climate Change[M].Cambridge, UK and New York, USA: Cambridge University Press.

    JIA Yan-kun, LIU Fu-liang, ZHANG Lin, NIE Zhen-long, SHEN Jian-mei, CHEN Zong-yu.2008.Recognition of Groundwater Recharge and Flow System in the Jiuquan-Zhangye Basin by Means of Environmental Isotopic Tracers[J].Acta Geoscientica Sinica, 29(06): 740-744(in Chinese with English abstract).

    LIN Xue-yu, FANG Yan-na, LIAO Zi-sheng, CHEN Hong-yan.2009.Impact of Global Warming and Human Activities on Groundwater Temperature[J].Journal of Beijing Normal University(Natural Science), 45(Z1): 452-457(in Chinese with English abstract).

    MIYAKOSHI A, TANIGUCHI M, OKUBO Y, UEMURA T.2005.Evaluations of subsurface flow for reconstructions of climate change using borehole temperature and isotope data in Kamchatka[J].Physics of The Earth and Planetary Interiors, 152(4): 305-313.

    NIE Zhen-long, GUO Zhan-rong, JIAO Peng-cheng, DONG Hua.2001.An Analysis of the Hydrological Cycle Characteristics in Interior Basins of Northwest China[J].Acta Geosicientica Sinica, 22(4): 302-306(in Chinese with English abstract).

    TAYLOR C A, STEFAN H G.2009.Shallow groundwater temperature response to climate change and urbanization[J].Journal of Hydrology, 39(3): 271-284.

    WANG Jian, LI Shuo.2005.Effect of climate change on snowmelt runoff in arid mountain areas of China[J].Science in China, Ser.D, 35(7): 664-670(in Chinese).

    WANG Jin-guo, ZHOU Zhi-fang, JIN Zhong-qing.2001.Coupled method of BEM-FAM for simulation of groundwater thermal transport[J].Journal of Hydraulic Engineering, (5): 71-76(in Chinese with English abstract).

    YANG Yong, CHEN Ren-sheng, JI Xi-bin.2007.Variations ofGlaciers in the Yeniugou Watershed of Heihe River Basin from 1956 to 2003[J].Journal of Glaciology and Geocryology, 29(1): 100-106(in Chinese with English abstract).

    ZHANG Guang-hui, CHEN Zong-yu, NIE Zhen-long, LIU Shao-yu, ZHANG Cui-yun, SHEN Jian-mei, WANG Jin-zhe, CHENG Xu-xue, ZHANG He-sheng.2006.Isotopic Characteristics of Groundwater and Its Response to Paleoclimatic Variation in Heihe River Basin[J].Acta Geoscientica Sinica, 27(4): 341-348(in Chinese with English abstract).

    ZHANG Guang-hui, LIU Shao-yu, XIE Yue-bo.2005a.Evolution Model of Water Cycle and Groundwater Formation in Northwest Inland Heihe River Basin[M].Beijing: Geological Publishing House: 60-66(in Chinese).

    ZHANG Guang-hui, NIE Zhen-long, LIU Shao-yu, CHEN Zong-yu, ZHANG Cui-yun, SHEN Jian-mei, WANG Jin-zhe, CHENG Xu-xue, ZHANG He-sheng.2005b.Characteristic and variation of groundwater recharge resources in the middle reaches of Heihe River basin[J].Advances in Water Science, 16(5): 673-678(in Chinese with English abstract).

    ZHANG Guang-hui, NIE Zhen-long, ZHANG Cui-yun, SHEN Jian-mei.2005c.Mechanism and characteristics of groundwater replenishment variation in middle reaches of Heihe River basin[J].Journal of Hydraulic Engineering, 36(6): 715-720(in Chinese with English abstract).

    中國地質(zhì)科學(xué)院2010年喜獲16項(xiàng)科技獎(jiǎng)

    2010年中國地質(zhì)科學(xué)院榮獲科技獎(jiǎng)16項(xiàng), 包括國家技術(shù)發(fā)明二等獎(jiǎng)1項(xiàng), 國土資源科學(xué)技術(shù)獎(jiǎng)一等獎(jiǎng)4項(xiàng)、二等獎(jiǎng)8項(xiàng), 其他省部級科技成果獎(jiǎng)3項(xiàng), 與2009年相比有較大幅度增長。

    礦產(chǎn)資源研究所張榮華研究員長期從事高溫高壓流體和反應(yīng)動(dòng)力學(xué)研究, 先后發(fā)明多種新的高溫高壓化學(xué)傳感器、有高溫高壓窗口的反應(yīng)裝置、高溫高壓檢測標(biāo)定實(shí)驗(yàn)平臺, 發(fā)展了高溫高壓環(huán)境原位直測技術(shù),在國際核心刊物發(fā)表多篇學(xué)術(shù)論文, 獲得6項(xiàng)國家發(fā)明專利和2項(xiàng)實(shí)用新型專利, 相關(guān)成果“高溫高壓流體和流動(dòng)反應(yīng)原位觀測裝置、方法、整合技術(shù)”榮獲2010年度國家發(fā)明二等獎(jiǎng)。

    地質(zhì)研究所楊經(jīng)綏研究員和許志琴院士團(tuán)隊(duì)“青藏高原地體拼合、碰撞造山及隆升機(jī)制”、國家地質(zhì)實(shí)驗(yàn)測試中心黃懷曾研究員團(tuán)隊(duì)“首都北京及周邊地區(qū)水、土環(huán)境污染機(jī)理與調(diào)控原理”、礦產(chǎn)資源研究所毛景文研究員團(tuán)隊(duì)“中國東部中生代多階段成礦的過程與背景”及地質(zhì)力學(xué)研究所參加完成的“長江三峽鏈子崖危巖體防治工程”分別榮獲2010年度國土資源科學(xué)技術(shù)獎(jiǎng)一等獎(jiǎng)。

    地球物理地球化學(xué)勘查研究所朱立新研究員團(tuán)隊(duì)“城市環(huán)境地球化學(xué)調(diào)查異常查證方法技術(shù)研究”、國家地質(zhì)實(shí)驗(yàn)測試中心饒竹研究員團(tuán)隊(duì)“多目標(biāo)地質(zhì)調(diào)查中主要有機(jī)物分析方法研究及應(yīng)用”、地質(zhì)研究所肖慶輝研究員團(tuán)隊(duì)“中國花崗巖重大地質(zhì)問題研究”、礦產(chǎn)資源研究所祝有海研究員團(tuán)隊(duì)“我國陸域永久凍土帶天然氣水合物資源遠(yuǎn)景調(diào)查”、礦產(chǎn)資源研究所楊建民研究員團(tuán)隊(duì)“天山銅礦帶找礦靶區(qū)優(yōu)選”、礦產(chǎn)資源研究所趙一鳴研究員團(tuán)隊(duì)“戰(zhàn)略性礦產(chǎn)單礦種編圖及成礦規(guī)律研究”、地質(zhì)研究所耿元生研究員團(tuán)隊(duì)“揚(yáng)子地臺西緣變質(zhì)基底演化”及礦產(chǎn)資源研究所參加完成的“冷水坑銀鉛鋅多金屬成礦系統(tǒng)與隱伏礦預(yù)測研究”分別榮獲2010年度國土資源科學(xué)技術(shù)獎(jiǎng)二等獎(jiǎng)。

    地質(zhì)研究所呂君昌完成的“河南省古生物化石地質(zhì)遺跡資源調(diào)查評價(jià)”獲河南省科學(xué)技術(shù)進(jìn)步獎(jiǎng) 2等獎(jiǎng); 巖溶地質(zhì)研究所完成的“中國水土流失與生態(tài)安全綜合科學(xué)考察”、“巖溶峰叢洼地水土保持研究與示范”分別獲得中國水土保持學(xué)會科學(xué)技術(shù)獎(jiǎng)一等獎(jiǎng)、二等獎(jiǎng)。

    本刊編輯部 采編

    Groundwater Temperature Variation in the Zhangye Inland Basin of Northwest China and Its Indications

    LIAN Ying-li1), ZHANG Guang-hui1), NIE Zhen-long1), AI Ting2)
    1) Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang, Hebei 050061;
    2) School of Economic Management, Hebei University of Engineering, Handan, Hebei 056038

    A study of groundwater temperature variation in Zhangye inland basin shows that the shallow groundwater and the confined groundwater decreased by 0.78℃ and 1.17℃ respectively in the study area from 1990 to 2009.This phenomenon was related to the increasing of melt water and precipitation caused by the air temperature rising in mountain areas.Meanwhile, the increase of groundwater exploitation led to the nearly same development trends of the temperature of the shallow groundwater and that of the confined groundwater in the study area.The groundwater temperature of the basin plain was in inverse proportion to the air temperature and precipitation in the mountain area.It means that, when the mountain temperature and precipitation increase the basin groundwater temperature decreases, and when the mountain temperature and precipitation decrease the basin groundwater temperature increases.70%-90% of the basin plain groundwater recharge which is at low temperature comes from the upstream of Qilian Mountains, and therefore the variation degree of basin plain groundwater temperature indicates the quantity of the effective recharge that the groundwater obtains, and has indicative significance to certain extent.

    P332; P641

    A

    10.3975/cagsb.2011.02.08

    本文由國家科技支撐計(jì)劃項(xiàng)目(編號: 2007BAD69B02; 2009BADA3B05)、國家自然基金項(xiàng)目(編號: 40602032)和河北省科技廳重點(diǎn)基礎(chǔ)研究項(xiàng)目(編號: 08966711D)聯(lián)合資助。

    2010-12-14; 改回日期: 2011-03-14。責(zé)任編輯: 閆立娟。

    連英立, 男, 1982年生。博士研究生。主要從事地下水資源可持續(xù)利用研究。E-mail: lianyingli123@163.com。

    張光輝, 男, 1959年生。研究員, 博導(dǎo)。長期從事區(qū)域水循環(huán)演化和地下水可持續(xù)利用研究。E-mail: Huanjing@heinfo.net。

    猜你喜歡
    承壓水潛水盆地
    潛水者
    中外文摘(2022年19期)2022-12-17 09:44:20
    地鐵深基坑承壓水控制研究
    深層承壓水污染途徑及防治研究
    盆地是怎樣形成的
    二疊盆地Wolfcamp統(tǒng)致密油成藏特征及主控因素
    快樂潛水
    潛水去
    鯨魚會得潛水病嗎?
    楚雄盆地扭動(dòng)構(gòu)造及其演化
    高承壓水上采煤可行性分析及安全開采評價(jià)
    金屬礦山(2013年7期)2013-03-11 16:54:07
    国产亚洲精品第一综合不卡| 人妻一区二区av| 国产成人a∨麻豆精品| 两个人看的免费小视频| 蜜桃国产av成人99| 久久久久久免费高清国产稀缺| 欧美另类一区| www.精华液| 爱豆传媒免费全集在线观看| 国产成人av教育| 人人妻人人澡人人看| 久久九九热精品免费| 成在线人永久免费视频| 久久国产精品大桥未久av| 久久久久精品人妻al黑| 日本精品一区二区三区蜜桃| 亚洲视频免费观看视频| 亚洲精品久久成人aⅴ小说| 好男人电影高清在线观看| 日韩一区二区三区影片| 亚洲一区中文字幕在线| 在线 av 中文字幕| 亚洲自偷自拍图片 自拍| 久久国产精品人妻蜜桃| www.精华液| 少妇的丰满在线观看| 亚洲伊人久久精品综合| 中国美女看黄片| 亚洲国产精品一区二区三区在线| 丝袜美腿诱惑在线| 午夜福利在线免费观看网站| 国产成人精品在线电影| 亚洲国产日韩一区二区| 国产91精品成人一区二区三区 | 岛国在线观看网站| 99九九在线精品视频| 夜夜骑夜夜射夜夜干| 亚洲欧美一区二区三区久久| 免费在线观看日本一区| 成人影院久久| 一级毛片电影观看| 91麻豆av在线| 2018国产大陆天天弄谢| 中文字幕色久视频| av电影中文网址| 久久精品亚洲熟妇少妇任你| 啦啦啦中文免费视频观看日本| av在线播放精品| 50天的宝宝边吃奶边哭怎么回事| 精品国产乱子伦一区二区三区 | 丝瓜视频免费看黄片| 人人妻,人人澡人人爽秒播| 男女午夜视频在线观看| 亚洲专区中文字幕在线| 丁香六月天网| 久久精品国产a三级三级三级| 欧美+亚洲+日韩+国产| 国产xxxxx性猛交| 一二三四社区在线视频社区8| a级片在线免费高清观看视频| 欧美激情极品国产一区二区三区| 国产亚洲欧美在线一区二区| 国产欧美亚洲国产| 日韩视频在线欧美| 久久九九热精品免费| 12—13女人毛片做爰片一| 女人久久www免费人成看片| 狠狠精品人妻久久久久久综合| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧洲精品一区二区精品久久久| 日本猛色少妇xxxxx猛交久久| 亚洲欧洲精品一区二区精品久久久| 国产高清国产精品国产三级| 高清在线国产一区| www.熟女人妻精品国产| 伦理电影免费视频| 久久久久久久久免费视频了| 色视频在线一区二区三区| 国产成人精品久久二区二区免费| 免费在线观看视频国产中文字幕亚洲 | 热99久久久久精品小说推荐| a级片在线免费高清观看视频| 日韩三级视频一区二区三区| 性色av乱码一区二区三区2| 亚洲精品美女久久久久99蜜臀| 亚洲一区二区三区欧美精品| 俄罗斯特黄特色一大片| 国产在线视频一区二区| 在线观看一区二区三区激情| 欧美+亚洲+日韩+国产| 亚洲 国产 在线| 国产精品麻豆人妻色哟哟久久| 在线亚洲精品国产二区图片欧美| 人人妻,人人澡人人爽秒播| 欧美av亚洲av综合av国产av| 午夜福利一区二区在线看| 十八禁网站网址无遮挡| 国产精品二区激情视频| 欧美日本中文国产一区发布| 大陆偷拍与自拍| 免费女性裸体啪啪无遮挡网站| 亚洲av国产av综合av卡| 国产精品成人在线| 99香蕉大伊视频| 老司机午夜福利在线观看视频 | 亚洲精品美女久久av网站| 高清视频免费观看一区二区| 五月开心婷婷网| 免费女性裸体啪啪无遮挡网站| 亚洲欧美一区二区三区黑人| 免费女性裸体啪啪无遮挡网站| 亚洲av国产av综合av卡| 国产极品粉嫩免费观看在线| 成年女人毛片免费观看观看9 | 国产精品国产三级国产专区5o| 成年av动漫网址| 精品亚洲乱码少妇综合久久| 巨乳人妻的诱惑在线观看| 久久久精品免费免费高清| 欧美精品人与动牲交sv欧美| 久久久久国内视频| 免费观看a级毛片全部| 国产精品秋霞免费鲁丝片| 首页视频小说图片口味搜索| av一本久久久久| 丝袜美足系列| 亚洲欧美精品综合一区二区三区| 婷婷成人精品国产| 青春草视频在线免费观看| 人妻一区二区av| 亚洲人成77777在线视频| 久久天躁狠狠躁夜夜2o2o| 欧美日韩视频精品一区| 国产一区二区在线观看av| 多毛熟女@视频| 久久影院123| 国产免费一区二区三区四区乱码| 成年人午夜在线观看视频| 日韩一区二区三区影片| 丝瓜视频免费看黄片| 精品国产一区二区三区四区第35| 久热这里只有精品99| 99re6热这里在线精品视频| 男女之事视频高清在线观看| 亚洲一区二区三区欧美精品| 久久香蕉激情| 国产亚洲午夜精品一区二区久久| xxxhd国产人妻xxx| 亚洲欧美清纯卡通| 精品少妇一区二区三区视频日本电影| 在线亚洲精品国产二区图片欧美| 波多野结衣av一区二区av| 欧美黄色片欧美黄色片| 俄罗斯特黄特色一大片| 久久免费观看电影| 免费不卡黄色视频| 国产精品欧美亚洲77777| 人人澡人人妻人| 91九色精品人成在线观看| 777久久人妻少妇嫩草av网站| 精品国产超薄肉色丝袜足j| 十八禁高潮呻吟视频| 亚洲精品久久成人aⅴ小说| 日韩三级视频一区二区三区| 母亲3免费完整高清在线观看| 最黄视频免费看| 午夜久久久在线观看| 九色亚洲精品在线播放| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区四区五区乱码| 成人国语在线视频| 一级毛片电影观看| 亚洲欧美清纯卡通| 色老头精品视频在线观看| 不卡一级毛片| 少妇精品久久久久久久| 国产精品二区激情视频| 国产成人精品久久二区二区91| 如日韩欧美国产精品一区二区三区| 日本wwww免费看| 日韩有码中文字幕| 波多野结衣一区麻豆| 亚洲国产av影院在线观看| 国产片内射在线| 欧美日韩国产mv在线观看视频| 亚洲五月色婷婷综合| 在线看a的网站| 中文字幕人妻丝袜制服| 亚洲av国产av综合av卡| 在线av久久热| 久久青草综合色| 超碰成人久久| 在线亚洲精品国产二区图片欧美| 国产淫语在线视频| 俄罗斯特黄特色一大片| 久久久国产成人免费| 亚洲avbb在线观看| 国产99久久九九免费精品| 一边摸一边抽搐一进一出视频| 韩国精品一区二区三区| 91麻豆精品激情在线观看国产 | 日本五十路高清| 自拍欧美九色日韩亚洲蝌蚪91| 老司机亚洲免费影院| 亚洲精品中文字幕在线视频| 欧美黑人精品巨大| 免费高清在线观看日韩| 成人黄色视频免费在线看| 久久久国产一区二区| 精品国产乱子伦一区二区三区 | 免费观看av网站的网址| 精品国产一区二区久久| 男女高潮啪啪啪动态图| 无限看片的www在线观看| 日韩制服丝袜自拍偷拍| 男女床上黄色一级片免费看| 亚洲黑人精品在线| 国产区一区二久久| 久久久精品94久久精品| 中文字幕精品免费在线观看视频| 亚洲av成人不卡在线观看播放网 | 男女高潮啪啪啪动态图| 色婷婷久久久亚洲欧美| 91成人精品电影| 欧美日韩中文字幕国产精品一区二区三区 | 国产国语露脸激情在线看| 欧美大码av| 丝袜脚勾引网站| 国产亚洲av片在线观看秒播厂| 午夜福利免费观看在线| 日本av免费视频播放| 青春草亚洲视频在线观看| av免费在线观看网站| 亚洲国产看品久久| 黄频高清免费视频| 操美女的视频在线观看| 国产免费视频播放在线视频| 亚洲国产欧美日韩在线播放| 久久av网站| 亚洲五月婷婷丁香| 制服人妻中文乱码| 九色亚洲精品在线播放| 精品视频人人做人人爽| 免费女性裸体啪啪无遮挡网站| 狂野欧美激情性xxxx| 国产成+人综合+亚洲专区| 午夜福利视频在线观看免费| 视频区欧美日本亚洲| 一区在线观看完整版| 久久综合国产亚洲精品| 一边摸一边抽搐一进一出视频| 亚洲欧美色中文字幕在线| 免费女性裸体啪啪无遮挡网站| e午夜精品久久久久久久| 亚洲情色 制服丝袜| 精品国产一区二区三区久久久樱花| a在线观看视频网站| 午夜成年电影在线免费观看| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品av久久久久免费| 操出白浆在线播放| 青青草视频在线视频观看| 老熟妇乱子伦视频在线观看 | 亚洲精品一二三| 满18在线观看网站| 日韩中文字幕视频在线看片| 热re99久久国产66热| 性色av一级| 宅男免费午夜| 最新的欧美精品一区二区| 正在播放国产对白刺激| 国产精品1区2区在线观看. | 最黄视频免费看| 午夜免费鲁丝| 欧美日韩亚洲高清精品| 日韩三级视频一区二区三区| 香蕉国产在线看| 中文字幕av电影在线播放| 在线亚洲精品国产二区图片欧美| 欧美日本中文国产一区发布| 一区在线观看完整版| 日本av手机在线免费观看| 久久久久久人人人人人| 亚洲精品粉嫩美女一区| 欧美精品一区二区免费开放| 色视频在线一区二区三区| 午夜福利影视在线免费观看| 韩国精品一区二区三区| 国产男女超爽视频在线观看| 日本一区二区免费在线视频| 免费一级毛片在线播放高清视频 | 99久久精品国产亚洲精品| 黄色毛片三级朝国网站| 啦啦啦啦在线视频资源| 美女高潮喷水抽搐中文字幕| 久久久久国内视频| 亚洲avbb在线观看| 一区二区三区乱码不卡18| 天天躁狠狠躁夜夜躁狠狠躁| 大陆偷拍与自拍| 黄色片一级片一级黄色片| 久久久精品免费免费高清| 国产男女超爽视频在线观看| 久久久久久亚洲精品国产蜜桃av| 一本久久精品| 搡老岳熟女国产| 亚洲一码二码三码区别大吗| 91麻豆av在线| 老司机深夜福利视频在线观看 | 久9热在线精品视频| 一区二区三区四区激情视频| 国产精品久久久久久精品古装| 亚洲,欧美精品.| 欧美精品一区二区免费开放| 久久久久久久精品精品| 欧美亚洲日本最大视频资源| 午夜福利在线免费观看网站| 久久久久国产一级毛片高清牌| 久久久久国产精品人妻一区二区| 国产av一区二区精品久久| 嫁个100分男人电影在线观看| 久久久久久久精品精品| 欧美+亚洲+日韩+国产| 俄罗斯特黄特色一大片| 午夜激情久久久久久久| 亚洲精品成人av观看孕妇| 日本精品一区二区三区蜜桃| 99久久国产精品久久久| 大片免费播放器 马上看| 欧美av亚洲av综合av国产av| 精品福利永久在线观看| 国产精品一区二区在线不卡| 午夜视频精品福利| 91av网站免费观看| 午夜成年电影在线免费观看| 日韩一区二区三区影片| 日韩,欧美,国产一区二区三区| 成人黄色视频免费在线看| 国产1区2区3区精品| 欧美日韩亚洲高清精品| 国产av精品麻豆| 少妇的丰满在线观看| 五月天丁香电影| 精品一品国产午夜福利视频| 黄片大片在线免费观看| 久久 成人 亚洲| 国产精品久久久av美女十八| 1024视频免费在线观看| 亚洲精品久久午夜乱码| 天天躁夜夜躁狠狠躁躁| 90打野战视频偷拍视频| 性少妇av在线| 欧美精品啪啪一区二区三区 | 欧美日本中文国产一区发布| 日本wwww免费看| 欧美性长视频在线观看| 男女无遮挡免费网站观看| 日韩大码丰满熟妇| 国产免费视频播放在线视频| cao死你这个sao货| 蜜桃国产av成人99| 女性生殖器流出的白浆| av又黄又爽大尺度在线免费看| 少妇 在线观看| 国产亚洲欧美在线一区二区| 亚洲精品国产一区二区精华液| 国产无遮挡羞羞视频在线观看| 美女脱内裤让男人舔精品视频| av又黄又爽大尺度在线免费看| 国产精品一二三区在线看| 亚洲av电影在线进入| av又黄又爽大尺度在线免费看| 欧美 日韩 精品 国产| 国产精品av久久久久免费| 欧美少妇被猛烈插入视频| 纯流量卡能插随身wifi吗| 伊人亚洲综合成人网| av在线播放精品| 日韩电影二区| 美女主播在线视频| 免费在线观看视频国产中文字幕亚洲 | 欧美黑人欧美精品刺激| 一本—道久久a久久精品蜜桃钙片| 国产亚洲一区二区精品| 国产男女超爽视频在线观看| 国产亚洲午夜精品一区二区久久| 少妇粗大呻吟视频| 天堂中文最新版在线下载| 我的亚洲天堂| 天天躁狠狠躁夜夜躁狠狠躁| 成人黄色视频免费在线看| 色婷婷久久久亚洲欧美| 捣出白浆h1v1| 欧美中文综合在线视频| 蜜桃国产av成人99| 最近最新中文字幕大全免费视频| 青青草视频在线视频观看| 午夜影院在线不卡| 免费观看a级毛片全部| 国产老妇伦熟女老妇高清| 亚洲国产日韩一区二区| 大码成人一级视频| 久久久久久久大尺度免费视频| 黄频高清免费视频| 亚洲国产精品一区三区| 无遮挡黄片免费观看| 成人手机av| 亚洲va日本ⅴa欧美va伊人久久 | 别揉我奶头~嗯~啊~动态视频 | cao死你这个sao货| 一区二区三区四区激情视频| 日韩电影二区| 欧美日韩福利视频一区二区| 久久人人爽人人片av| 免费高清在线观看日韩| 丰满人妻熟妇乱又伦精品不卡| 欧美午夜高清在线| 亚洲熟女毛片儿| 90打野战视频偷拍视频| 婷婷丁香在线五月| 热re99久久精品国产66热6| 90打野战视频偷拍视频| 久久久久精品国产欧美久久久 | 亚洲国产精品999| 桃红色精品国产亚洲av| 欧美乱码精品一区二区三区| 日韩欧美一区二区三区在线观看 | 黑人操中国人逼视频| h视频一区二区三区| 99久久国产精品久久久| 侵犯人妻中文字幕一二三四区| 19禁男女啪啪无遮挡网站| 国产亚洲午夜精品一区二区久久| 黑人巨大精品欧美一区二区蜜桃| 一级,二级,三级黄色视频| 91精品伊人久久大香线蕉| 国产欧美日韩一区二区三 | 高潮久久久久久久久久久不卡| 无遮挡黄片免费观看| 亚洲美女黄色视频免费看| 岛国在线观看网站| 捣出白浆h1v1| 国产精品 欧美亚洲| 午夜免费成人在线视频| 这个男人来自地球电影免费观看| 亚洲欧美日韩高清在线视频 | 欧美精品高潮呻吟av久久| 1024视频免费在线观看| 欧美乱码精品一区二区三区| 另类亚洲欧美激情| 桃花免费在线播放| 悠悠久久av| 男女免费视频国产| 亚洲中文日韩欧美视频| 欧美老熟妇乱子伦牲交| 夜夜夜夜夜久久久久| 丝袜在线中文字幕| av不卡在线播放| 国产在线视频一区二区| 老司机午夜福利在线观看视频 | 国产亚洲欧美精品永久| 12—13女人毛片做爰片一| 亚洲一码二码三码区别大吗| 亚洲 欧美一区二区三区| 国产日韩欧美在线精品| 一进一出抽搐动态| 亚洲中文字幕日韩| 老司机影院毛片| 黄色怎么调成土黄色| 久久狼人影院| 欧美精品人与动牲交sv欧美| 一级毛片电影观看| 青春草亚洲视频在线观看| 日本wwww免费看| 午夜福利在线观看吧| 日韩视频一区二区在线观看| 亚洲精品国产精品久久久不卡| 国产又爽黄色视频| 亚洲人成77777在线视频| 国产成人av教育| 国产亚洲一区二区精品| 亚洲人成77777在线视频| 91老司机精品| 日韩电影二区| 视频区图区小说| 水蜜桃什么品种好| 精品免费久久久久久久清纯 | 欧美日韩亚洲综合一区二区三区_| 国产片内射在线| 亚洲第一av免费看| 亚洲av男天堂| 成年av动漫网址| 亚洲五月婷婷丁香| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看| 中文字幕高清在线视频| www.999成人在线观看| 9色porny在线观看| 人妻 亚洲 视频| 亚洲国产av影院在线观看| 中文字幕制服av| 中文欧美无线码| 午夜两性在线视频| 免费高清在线观看视频在线观看| 精品福利观看| 高清黄色对白视频在线免费看| 一本久久精品| 俄罗斯特黄特色一大片| 精品少妇内射三级| 欧美日本中文国产一区发布| 老熟妇仑乱视频hdxx| 美女视频免费永久观看网站| 久久久久久久久久久久大奶| av电影中文网址| 精品久久久久久久毛片微露脸 | 丰满饥渴人妻一区二区三| 男女边摸边吃奶| 在线 av 中文字幕| 狠狠婷婷综合久久久久久88av| 美女福利国产在线| 岛国毛片在线播放| 亚洲情色 制服丝袜| 在线永久观看黄色视频| 法律面前人人平等表现在哪些方面 | 一本一本久久a久久精品综合妖精| 精品一区在线观看国产| 一区二区日韩欧美中文字幕| 一级片'在线观看视频| 天天影视国产精品| 亚洲综合色网址| 国产成人系列免费观看| av又黄又爽大尺度在线免费看| 久久久国产欧美日韩av| 女警被强在线播放| 老汉色∧v一级毛片| 欧美精品一区二区免费开放| 韩国高清视频一区二区三区| 国产黄色免费在线视频| 久久国产精品影院| 老鸭窝网址在线观看| 国产精品99久久99久久久不卡| 肉色欧美久久久久久久蜜桃| 老司机亚洲免费影院| 亚洲欧美一区二区三区黑人| 我要看黄色一级片免费的| 高清在线国产一区| 欧美97在线视频| 国产男人的电影天堂91| 巨乳人妻的诱惑在线观看| 国产成人欧美在线观看 | 国产高清国产精品国产三级| kizo精华| 欧美少妇被猛烈插入视频| 他把我摸到了高潮在线观看 | av在线播放精品| 国产精品一区二区在线观看99| 少妇猛男粗大的猛烈进出视频| 久久精品人人爽人人爽视色| 亚洲中文日韩欧美视频| 俄罗斯特黄特色一大片| 一本大道久久a久久精品| 欧美日韩亚洲综合一区二区三区_| 一本大道久久a久久精品| 在线看a的网站| 国产精品久久久久久精品电影小说| 亚洲av日韩在线播放| 欧美激情久久久久久爽电影 | 天天躁狠狠躁夜夜躁狠狠躁| 精品欧美一区二区三区在线| 男女无遮挡免费网站观看| 亚洲国产中文字幕在线视频| 国产精品香港三级国产av潘金莲| 法律面前人人平等表现在哪些方面 | 亚洲七黄色美女视频| 老司机影院成人| 1024视频免费在线观看| 另类精品久久| 波多野结衣一区麻豆| 91字幕亚洲| 久久久久久亚洲精品国产蜜桃av| 啦啦啦免费观看视频1| 国产精品 欧美亚洲| 亚洲男人天堂网一区| 国精品久久久久久国模美| 国产成人精品在线电影| 欧美 亚洲 国产 日韩一| 啪啪无遮挡十八禁网站| 亚洲熟女毛片儿| 成人黄色视频免费在线看| 韩国精品一区二区三区| 操美女的视频在线观看| 日韩精品免费视频一区二区三区| 久久午夜综合久久蜜桃| 一本一本久久a久久精品综合妖精| 亚洲精品av麻豆狂野| 亚洲性夜色夜夜综合| 免费av中文字幕在线| 蜜桃在线观看..| 亚洲色图综合在线观看| 精品久久久久久久毛片微露脸 | 国产日韩欧美亚洲二区| 亚洲avbb在线观看| 午夜福利,免费看| 99国产精品99久久久久| 久久人妻熟女aⅴ| 亚洲欧美精品综合一区二区三区| 后天国语完整版免费观看| 亚洲精品国产av蜜桃| 可以免费在线观看a视频的电影网站| 亚洲精品久久成人aⅴ小说| 黑人欧美特级aaaaaa片| 成人免费观看视频高清| 午夜福利影视在线免费观看| 欧美精品人与动牲交sv欧美| 不卡一级毛片|