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

    新疆博阿斷裂附近地表水和地下水的水化學(xué)和同位素特征及水質(zhì)評價

    2023-11-22 03:40:20楊明遠(yuǎn)趙佳怡馬超李鑫
    西北地質(zhì) 2023年6期
    關(guān)鍵詞:大氣研究

    楊明遠(yuǎn) ,趙佳怡 ,馬超 ,李鑫

    (1. 新疆維吾爾自治區(qū)地質(zhì)礦產(chǎn)勘察開發(fā)局第二區(qū)域地質(zhì)調(diào)查大隊,新疆 昌吉 831100;2. 中國地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所,河北 石家莊 050061;3. 自然資源部地?zé)崤c干熱巖勘察開發(fā)技術(shù)創(chuàng)新中心,河北 正定 050800)

    水資源緊缺一直是制約社會可持續(xù)發(fā)展的重要因素,尤其是中國西北干旱地區(qū)地下水資源匱乏,研究地下水的水文地球化學(xué)特征及水質(zhì)可為水資源的開發(fā)利用提供支持。

    水文地球化學(xué)的研究最早起始于地下水勘探開發(fā)利用的水文地質(zhì)研究,自20世紀(jì)60、70年代起,水文地球化學(xué)在地下水監(jiān)測及分析中應(yīng)用越來越廣泛(沈照理,1985)。20世紀(jì)60年代,開始研究地下水在地殼中的作用(任福弘等,1993),召開了相應(yīng)的水文地球與生物地球化學(xué)研討會,研究水文地球化學(xué)相關(guān)問題。1992年《水文地球化學(xué)》將此研究分為理論與應(yīng)用兩部分,從地下水形成、遷移等理論問題和飲用、地?zé)崴?、礦水和模擬預(yù)測的應(yīng)用兩方面進(jìn)行論述(葉思源等,2002)。近期水文地球化學(xué)圍繞油田水、海底礦藏、地?zé)崴?、環(huán)境等方面進(jìn)行特定的分析,研究手段也不斷增加(馬致遠(yuǎn)等,2017;劉瑞平等,2019)。

    新疆地質(zhì)構(gòu)造復(fù)雜多樣,地殼活動頻繁,地勢具有南高北低、西高東低少環(huán)山封閉盆地的特點(diǎn),地?zé)豳Y源豐富,各熱水區(qū)帶的水熱活動強(qiáng)度自北而南逐漸增強(qiáng),自西向東逐漸減弱(中國各省地?zé)豳Y源全解析,2020)。王書峰(1985)對新疆烏魯木齊河流域的水文地球化學(xué)進(jìn)行了初步研究,周金龍(2010)將新疆地下水劃分7個方面對其進(jìn)行整體分析研究,為新疆等西北干旱地區(qū)地下水資源與環(huán)境問題研究提供了可借鑒的方法和理念。顧新魯?shù)龋?015)對新疆地?zé)豳Y源成因類型及控?zé)崮J竭M(jìn)行了分析,將地?zé)岢梢蝾愋头譃閿嗔褬?gòu)造型、沉積盆地型和煤層自燃型,主要受區(qū)域斷裂構(gòu)造及次級構(gòu)造、大地?zé)崃鞅尘昂蜔嵩?、水源流體等因素的影響和控制。殷秀蘭等(2015)對烏魯木齊河流域北部的水文地球化學(xué)特征進(jìn)行了研究。蔣萬軍等(2016)、趙江濤等(2017)、雷米等(2020)分別對新疆吐—哈盆地、焉耆盆地平原區(qū)和昌吉州東部平原區(qū)的地下水行了水文地球化學(xué)分析,對地下水演化及分布特征進(jìn)行整體研究。

    新疆地質(zhì)演化歷史復(fù)雜,博羅科努-阿其克庫都克斷裂(簡稱“博阿斷裂”)是超巖石圈斷裂,有長期和多期活動的特點(diǎn),是良好的導(dǎo)熱導(dǎo)水通道及地下水容集通道,形成條件和巖性等都影響著其蓄水能力和地下水流場。因此,通過研究博阿斷裂附近的地表水及地下水的水化學(xué)等特征可判斷其對地下水演化、流場等影響??傮w來看,針對新疆單一地區(qū)的水文地球化學(xué)分析相對較多,但缺乏針對整體構(gòu)造的不同分布位置的水的水文地球化學(xué)研究。新疆屬于西北干旱地區(qū),水資源匱乏,通過研究博阿斷裂帶上泉水的水文地球化學(xué)特征可以更好分析斷裂構(gòu)造上對地下水化學(xué)演化、徑流的影響,為開發(fā)利用地下水及進(jìn)行水利工程建設(shè)等提供理論研究支撐。

    1 研究區(qū)概況

    博阿斷裂是一條巨大的NW向構(gòu)造,斷裂起始于哈薩克斯坦境內(nèi)沿阿拉胡西南側(cè)的阿拉套山山前,從阿拉山口由NW向SE進(jìn)入中國,在中國境內(nèi)長約為800 km(圖1;楊曉平等,2000)。博阿斷裂屬于超巖石圈斷裂(陳哲夫等,1997),屬準(zhǔn)噶爾-哈薩克斯坦古板塊與塔里木古板塊結(jié)合部位,次級斷裂及褶皺構(gòu)造十分發(fā)育,是褶皺山地斷裂型地?zé)豳Y源有利賦存地區(qū),以密集的NWW向大型右行壓扭性斷裂為主體,為天山山脈主要控?zé)針?gòu)造。從該斷裂控制的地層和侵入巖體看,它具有長期和多期活動的特點(diǎn),早古生代以前即已形成。中新生代以來,由于天山強(qiáng)烈隆起而導(dǎo)致斷裂重新活動,形成由南向北推擠的疊瓦狀斷裂帶,斷面向西南傾,傾角為50°~60°。沿斷裂發(fā)育破碎帶和動力變質(zhì)帶,其中片理、擠壓透鏡體長軸及小型褶皺曲軸面與主斷裂面成15°夾角,顯示斷裂的壓扭性質(zhì)(馮先岳,1997)。

    圖1 天山地區(qū)主要活動斷裂及采樣點(diǎn)位置(楊曉平等,2000)Fig. 1 Main active faults in Tianshan area and sampling point location

    研究區(qū)位于北天山優(yōu)地槽褶皺帶,西部位于依連哈比爾尕復(fù)背斜帶,受區(qū)域構(gòu)造控制,形成一系列NWW、SEE向壓性及壓扭性斷裂,NE—SW向和NW—SE向扭斷裂。置屬中山區(qū)地貌,山勢陡峻,溝谷縱橫交匯,河谷深切,呈“U”型谷。侵蝕作用強(qiáng)烈,溫泉出露部位地勢相對西高東低,南高北低,總體地勢是南高北低(陳鋒等,2016)。西部地區(qū)主要發(fā)育地層為中石炭統(tǒng)前峽亞組,由沉積作用、火山噴發(fā)作用形成的一套以淺海-濱海相正常沉積巖、火山碎屑沉積巖、火山碎屑巖建造,地層多呈單斜產(chǎn)出(陳首,2017)。中部位于昌吉市頭屯河水庫—三屯河努爾加水庫,地處天山北麓準(zhǔn)噶爾盆地南緣的頭屯河—三屯河沖洪積扇的山前傾斜平原上,地勢總趨勢西南高北東低,地表被兩河水系、沖溝侵蝕切割,切割深度自南向北逐漸變小。頭屯河發(fā)育有三至四級堆積階地,并有河漫灘發(fā)育,切割深度為25~35 m,河谷寬為400~600 m;三屯河發(fā)育有四級階地,切割深度為30~40 m,河岸寬為300~600 m,并有河漫灘發(fā)育。中部地區(qū)巖性單一,巖性為卵石、礫石、砂礫石、砂層等,結(jié)構(gòu)松散,透水性好。地形復(fù)雜,總體上呈西高東低,南高北低,海拔高程為1 095~1306 m,相對高差為275 m,山體地形走向以NW向的溝谷為主。受控于地質(zhì)構(gòu)造運(yùn)動、三屯河水系及頭屯河水系的徑流切割,以低山丘陵地貌為主,臺地地貌為次,局部地勢切割劇烈。東南地區(qū)位于烏魯木齊河流域,屬于天山北坡中段,準(zhǔn)噶爾盆地南緣,北部為傾斜平原,東望可見天山主峰博格達(dá)峰,南依天山中段天格爾峰,NW向準(zhǔn)噶爾盆地傾斜,區(qū)廣泛分布多層結(jié)構(gòu)的潛水、承壓水和深部自流水(周金龍,2010)。

    2 取樣與檢測方法

    在位于博阿斷裂帶的頭屯河水庫—三屯河努爾加水庫地區(qū)共采集地表水水樣3組,在同位于博阿斷裂帶的沙灣縣、呼圖壁縣及烏魯木齊市采集溫泉水樣6組。另外收集位于博阿斷裂帶的烏魯木齊河泉水水樣7組(表1;高朋等,2021)。

    采集的泉水樣品進(jìn)行了水質(zhì)全分析,部分溫泉樣品進(jìn)行了氫氧同位素分析,水質(zhì)全分析根據(jù)中華人民共和國國家標(biāo)準(zhǔn)飲用天然礦泉水檢驗方法(GB/T8538-2008)檢測,其中陽離子使用的ICP方法進(jìn)行離子濃度檢測,陰離子使用離子色譜進(jìn)行分析,陰陽離子平衡誤差控制在3 %以內(nèi)。氫氧同位素利用波長掃描—光腔衰蕩光譜法在溫度23 ℃,濕度50 %情況下進(jìn)行檢測。本次樣品均由國土資源部地下水礦泉水及環(huán)境檢測中心(中國地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所)檢測完成。

    3 結(jié)果分析與討論

    3.1 水化學(xué)特征分析

    (1)pH和TDS

    研究區(qū)溫泉水、地表水和冷泉水的水化學(xué)特征有明顯的區(qū)別。從溫泉水、地表水到地下水,pH值逐漸降低,分別為8.82~9.92、8.01~8.58和7.6~7.97,且均大于7.5,為堿性水。溫泉水、地表水和冷泉水的TDS含量為77~3 633.91 mg/L,且從地表水、溫泉水到冷泉水,TDS含量逐漸降低。其中,溫泉水和冷泉水的TDS分別為304.75~553.73 mg/L和77~170 mg/L,均小于1g/L,為淡水。然而,地表水的TDS為1 338.64~3 633.91 mg/L,大于1 g/L或 3g/L,為微咸水和咸水。

    (2)水化學(xué)特征及主要離子含量

    不同類型水的主要離子含量和水化學(xué)類型有明顯的區(qū)別(張帆等,2023)。地表水中的各種主要離子含量都明顯要高于溫泉水和冷泉水,而且溫泉水和冷泉水中各個主要離子的相對含量大小趨勢并不存在一致性(圖2)。由piper三線圖(圖3)可知地表水的水化學(xué)類型主要為SO4—Na、Cl·SO4—Na和HCO3·SO4—Ca·Na型,Na+為主要的陽離子,其含量為332.22~1 081.49 mg/L,SO42-為主要的陰離子,其含量為537.86~198.26 mg/L。溫泉水的水化學(xué)類型為SO4·HCO3—Na/HCO3·SO4—Na和HCO3·Cl—Na型,Na+為主要的陽離子,其含量為100.89~163.12 mg/L,HCO3—為最主要的陰離子,其含量為23.19~48.83 mg/L,此外,SO42—和Cl—為溫泉水的主要陰離子,含量分別為15.36~158.44 mg/L和14.32~96.63 mg/L。冷泉水的水化學(xué)類型為SO4·HCO3—Mg·Ca、HCO3—Ca、HCO3—Mg·Ca和SO4·HCO3-Ca型,SO42—和HCO3—為主要的陰離子,含量分別在13.78~28.79 mg/L和46.75~169.47 mg/L,Ca2+和Mg2+為主要的陽離子,含量分別為21.87~50.18 mg/L和3.39-8.7 mg/L。此外,溫泉水的Ca2+和Mg2+含量要顯著小于地表水和冷泉水,含量僅分別為1.6~4.01mg/L和0.24~0.73 mg/L。同時,冷泉水的Na+和Cl—含量明顯偏低,分別為0.55~3.68 mg/L和3.42~9.4 mg/L。

    圖2 研究區(qū)主要離子濃度圖Fig. 2 Major ion concentration in the study area

    圖3 研究區(qū)piper三線圖Fig. 3 Piper map of the study area

    (3)特征組分

    筆者主要對NO3-、NH4+、F-和偏硅酸進(jìn)行討論。由于并未收集到冷泉水的特征組分的數(shù)據(jù),因此只討論溫泉水和地表水中這些特征組分的含量特征。從圖4可以看到,溫泉水中的偏硅酸含量要顯著高于地表水,其在溫泉水和地表水中的含量分別為40.86~133.43 mg/L和16.75~37.66 mg/L。溫泉水的F-含量也要高于地表水,分別為2.43~10.97 mg/L和1.82~2.76 mg/L,并且DRT4和DRT5水樣的F-含量要顯著偏高分別為10.97 mg/L和9.67 mg/L。溫泉水中的NO3-含量顯著低于地表水,分別為0.14~5.74 mg/L和14.29~31.84 mg/L。然而,溫泉水的NH4+含量則明顯高于地表水,分別為0.39~5.7 mg/L和<0.04~2.2 mg/L。

    圖4 研究區(qū)特征組分濃度圖Fig. 4 Characteristic component concentration of study area

    3.2 水文地球化學(xué)過程

    地表水和地下水經(jīng)歷的不同的水-巖相互作用、混合作用、循環(huán)時間和路徑以及所處的氧化還原環(huán)境等因素都控制著水中的離子或元素的含量(張錫根,1988;張保健,2011;張未,2016; Wang et al.,2020)。因此,闡明地下水和地表水在循環(huán)過程中所經(jīng)歷的水文地球化學(xué)過程對了解地下水和地表水的化學(xué)組分的形成機(jī)理十分重要。

    (1)蒸發(fā)作用

    通過對河流、湖泊和雨水等樣本分析,Gibbs(1970)提出了地表水的3種離子來源:蒸發(fā)結(jié)晶、巖石風(fēng)化和降水(韓朝輝等,2023)。通過分析TDS vs. Na+/(Na++Ca2+)和TDS vs. Cl—/(Cl—+HCO3—)(圖5),可以發(fā)現(xiàn)地表水水樣點(diǎn)大部分位于蒸發(fā)結(jié)晶區(qū)域,冷泉位于巖石風(fēng)化區(qū)域,溫泉位于巖石風(fēng)化域和蒸發(fā)區(qū)域,說明研究區(qū)地表水受到了蒸發(fā)作用的強(qiáng)烈影響,冷泉和溫泉主要受到巖石風(fēng)化或水巖作用的影響。因此,地表水的主要離子含量與溫泉水和冷泉水相比顯著偏高。同時,地表水水樣點(diǎn)在TDS vs. Cl—/(Cl—+HCO3—)圖中位于Gibbs圖偏左,而在TDS vs. Na+/(Na++Ca2+)圖中位于Gibbs圖右側(cè),這可能是由于Ca2+發(fā)生了沉淀導(dǎo)致了其含量的降低,或者Na+與Ca2+發(fā)生了陽離子交換反應(yīng),導(dǎo)致Na+含量增加而Ca2+含量降低。

    圖5 研究區(qū)Gibbs圖Fig. 5 Gibbs map of the study area

    (2)礦物溶解與沉淀

    一般來說,地下水和地表水中的溶解組分的含量主要來自圍巖的溶濾、大氣降水以及人類活動等(Lasaga et al.,1994; Xing et al.,2013),礦物的溶解和沉淀會導(dǎo)致水化學(xué)類型發(fā)生變化(Gao et al.,2022)。

    運(yùn)用PHREEQC對溫泉水、地表水和冷泉水的部分礦物飽和指數(shù)(SI)進(jìn)行了模擬計算(圖6),發(fā)現(xiàn)硬石膏、螢石和石膏的SI值都<0,且大部分都<—1,說明研究區(qū)這些水中的硬石膏、螢石和石膏都處于溶解狀態(tài)。然而,方解石、白云石和文石的礦物飽和指數(shù)大部分為—1~1,表明研究區(qū)的這些水中的方解石、白云石和文石處于平衡和沉淀狀態(tài),有發(fā)生沉淀的趨勢。同時,溫泉水中的方解石、白云石和文石基本上都>0,處于沉淀狀態(tài),這也解釋了其水樣點(diǎn)在(Na+K)vs.HCO3圖中位于1∶1直線上方的原因。這些礦物對地下水水文地球化學(xué)過程起了相當(dāng)重要的作用。

    圖6 溫泉水(a)、地表水(b)和冷泉水(c)SI圖Fig. 6 SI map of (a) hot spring, (b) surface water and (c) cold spring

    通過Ca2+與Mg2+的比值分析,可更好了解方解石與白云石的溶解與沉淀(Gao et al., 2020)。作研究區(qū)Ca2+vs. Mg2+圖(圖7)顯示,水樣基本落在1∶1直線的下方,說明研究區(qū)中方解石比白云巖更占優(yōu)勢。

    圖7 研究區(qū)Ca2+ vs Mg2+離子關(guān)系圖Fig. 7 Diagram of Ca2+ vs Mg2+ relationship in the study area

    (3)陽離子交替吸附

    通過主要離子的毫克當(dāng)量數(shù)據(jù)點(diǎn)以及礦物溶解比直線可以揭示主要離子的主要來源。

    從HCO3—vs.(Ca2++Mg2+)、(Ca2++Mg2+)vs.(SO42—+HCO3—)、HCO3—vs. Ca2+、Na+vs. Cl—以及(Na++K+)vs.HCO3—的主要離子關(guān)系圖(圖8)可以看出,地表水樣點(diǎn)分布都十分分散,表明地表水中的這些元素來自于多種礦物的溶解。溫泉水水樣點(diǎn)在HCO3—vs. (Ca2++Mg2+)、(Ca2++Mg2+)vs.(SO42—+HCO3—)、HCO3—vs. Ca2+圖中都位于1∶1直線下方,說明有其他的HCO3—來源,而在Na+vs. Cl—以及(Na++K+)vs. HCO3—圖位于1∶1直線上方,說明Na+有其他的來源。通常在非碳酸鹽巖熱儲中的地下熱水中,由于深部CO2的影響,長石類礦物會發(fā)生溶解,導(dǎo)致地下熱水中有較高含量的Na+和HCO3—,同時由于少量鈣長石的溶解,因此地下熱水中具有少量的Ca2+。同時,含Mg2+的礦物很少,因此Mg2+的含量很低,低于地下冷水和地表水??梢杂萌缦禄瘜W(xué)方程式表示:

    圖8 主要離子關(guān)系圖Fig. 8 Major ion relationship map

    冷泉水樣點(diǎn)在(Ca2++Mg2+)vs.(SO42—+HCO3—)圖中剛好位于1∶1溶解比直線,這說明冷泉水中的Ca2+、Mg2+、SO42—和HCO3—主要來自于白云石、方解石以及石膏的溶解,可用如下化學(xué)方程式來表示:—vs. Ca2+以及SO4vs. Ca2+圖中都位于1∶1直線的上方,也說明了冷泉中這些元素的含量不止一種來源,表明離子交換作用對地下水化學(xué)成分影響較為強(qiáng)烈。

    (4)氧化還原環(huán)境

    NO3—和NH4+可以作為氧化還原環(huán)境的指示劑。在有氧條件下,亞硝酸細(xì)菌和硝酸細(xì)菌能夠?qū)H4+最終氧化成為NO3—,可用如下化學(xué)方程式來表達(dá):

    其數(shù)據(jù)點(diǎn)在HCO3—vs.(Ca2++Mg2+)、HCO3

    通過上述特征組分含量的分析可以知道,地表水中的NO3-含量遠(yuǎn)大于溫泉水,但是其NH4+含量卻小于溫泉水,這說明地表水處于相對氧化的環(huán)境,促進(jìn)了NH4+在亞硝酸細(xì)菌和硝酸細(xì)菌的作用下,被氧化成為NO3—。地表水與空氣直接接觸,空氣中的O2可以溶解進(jìn)入地表水中,因而使得地表水處于氧化狀態(tài)。同時,這也說明了溫泉水處于相對的還原狀態(tài)或者較弱的氧化狀態(tài),因此其NH4+含量要遠(yuǎn)高于地表水。溫泉水的循環(huán)路徑較長,在接受大氣降水補(bǔ)給之后,往往經(jīng)歷較長時間的深循環(huán),在循環(huán)過程中水中的大氣降水中的O2就會消耗導(dǎo)致含量降低,如果沒有其他O2的補(bǔ)給來源,就會導(dǎo)致溫泉水處于相對還原或者弱氧化的環(huán)境,因而較少的NH4+被亞硝酸細(xì)菌和硝酸細(xì)菌的氧化成為NO3—。

    3.3 氫氧同位素特征

    (1) 補(bǔ)給來源

    研究地下水氫氧同位素特征可以判斷地下水的起源(張帆等,2023;王新娟等,2023),確定地下水的補(bǔ)給條件、大氣降水與地表水和地下水的聯(lián)系程度,了解地下水的循環(huán)途徑(Song et al.,2017)。由于D和18O的蒸汽張力比較小,因此其在液相中富集,在氣相中貧化,從而導(dǎo)致不同的地下水循環(huán)氫氧同位素含量不同(張錫根,1988)。地下水中的δD的變化值除了少部分的混合作用的影響外,主要取決于補(bǔ)給溫度及補(bǔ)給高程,δ18O的變化則主要根據(jù)水-巖作用的交換程度和水、巖比值(Craig,1961)。

    將水樣點(diǎn)的氫氧穩(wěn)定同位素數(shù)據(jù)和全球大氣降水線以及當(dāng)?shù)卮髿饨邓€繪制在同一張圖中(圖9),判斷水樣點(diǎn)與大氣降水線的相對位置,可以揭示研究區(qū)水的補(bǔ)給來源。其中全球大氣降水線:δD=8δ18O+10(倪高倩等,2016),新疆烏魯木齊大氣降水線:δD=7.215δ18O+4.5(李暉等,2009)。由于缺乏研究區(qū)地表水的氫氧穩(wěn)定同位素數(shù)據(jù),因此只判斷研究區(qū)溫泉水和冷泉水的補(bǔ)給來源。

    圖9 研究區(qū)樣品δD-δ18O關(guān)系圖Fig. 9 Study area samples δD-δ18O diagram

    研究區(qū)溫泉水和冷泉水水樣點(diǎn)都落在大氣降水線附近,表明它們的主要補(bǔ)給來源為大氣降水。然而,溫泉水和冷泉水水樣點(diǎn)的分布位置有較大區(qū)別。溫泉水分布在全球大氣降水線和當(dāng)?shù)卮髿饨邓€的下方,且出現(xiàn)了“18O漂移”的現(xiàn)象,這主要是因為溫泉水與圍巖發(fā)生的水-巖相互作用導(dǎo)致溫泉水和圍巖發(fā)生了δ18O的交換,導(dǎo)致水中δ18O富集(Wang et al.,2019)。冷泉水有多數(shù)水樣點(diǎn)分布在全球大氣降水線和當(dāng)?shù)卮髿饨邓€的上方,這主要是因為發(fā)生了水汽再循環(huán),導(dǎo)致冷泉水中的δD富集。

    溫泉水相對于冷泉水表現(xiàn)出δ18O和δD貧化以及更高的溫度的特征,這說明溫泉水比冷泉水可能接受了更高的大氣降水補(bǔ)給,以及經(jīng)歷了更深的水循環(huán),其在含水層中滯留的時間更長,徑流相對較長,水-巖相互作用程度更加強(qiáng)烈。

    (2) 補(bǔ)給高程

    通常來說,大氣降水的δD和δ18O與氣溫呈線性關(guān)系,高程效應(yīng)指的就是海拔越高,溫度越低,δD和δ18O值同時也降低(孫占學(xué)等,1992)。主要原因是由于高程的增加導(dǎo)致氣溫降低,水汽冷凝同位素分餾,雨滴蒸發(fā)減少,因此伴有山體迎面風(fēng)的雨量增加,形成了山體迎面風(fēng)的超前降雨區(qū),從而隨著云團(tuán)上升,同位素含量減少,形成高程效應(yīng)(陳禮明,2019)。通過研究這種高程效應(yīng),可根據(jù)氫氧同位素推測地下水補(bǔ)給高程。根據(jù)資料可知地下水補(bǔ)給高程公式:

    式中:δG為樣品中δD的值;δP為大氣降水中δD的值;k為大氣降水中δD的高程梯度;h為取樣點(diǎn)高程(m)。根據(jù)李學(xué)禮等(2000)推理得出研究區(qū)大氣降水的同位素成分,大氣降水δD值為—44.61,高程梯度為—0.046‰/100 m。

    根據(jù)計算公式得到研究區(qū)溫泉水補(bǔ)給高程為3 064~4 161 m,冷泉補(bǔ)給高程為2 874~3 943 m(表2),南部天山山脈海拔3 800~4 200 m以上為常年積雪和現(xiàn)代冰川作用帶,研究區(qū)補(bǔ)給高程與附近山體融雪線高程較為一致,說明研究區(qū)泉水點(diǎn)的主要來源是周圍山脈的大氣降水與冰川融雪,通過斷裂,徑流、排泄等過程到達(dá)研究區(qū),經(jīng)過靜水壓力等作用下上升至地表形成泉水。

    表2 研究區(qū)補(bǔ)給高程表Tab. 2 Recharge elevation in the study area

    3.4 水質(zhì)評價

    F值評分法是中國水質(zhì)評價最常用的一種方法,適用于化學(xué)指標(biāo)分級明顯、無超標(biāo)或多項指標(biāo)均超標(biāo)等各類情況(王一凡等,2014)。利用F值評分法對研究區(qū)的溫泉水和地表水的水質(zhì)進(jìn)行評價。具體步驟如下(孟春芳等,2015):

    (1)首先對各個離子進(jìn)行單項組分評價,得到單項組分評價分值(表3)。

    表3 地下水質(zhì)量評分表Tab. 3 Groundwater quality scoring table

    (2)根據(jù)公式計算綜合評價分值F:

    (3)根據(jù)綜合評分值F確定水質(zhì)分級(表4)。

    表4 地下水質(zhì)量分級表Tab. 4 Groundwater quality classification table

    采用pH、TDS、Cl—、SO42—、NO3—、NO2—、NH4+、F等8種離子指標(biāo)進(jìn)行水質(zhì)分析,結(jié)果見表5。

    表5 研究區(qū)水質(zhì)F分值法評價結(jié)果表Tab. 5 Evaluation results of F-score method for water quality in the study area

    根據(jù)表5發(fā)現(xiàn),研究區(qū)溫泉水和地表水的水質(zhì)都大于7.2,屬于V類水,水質(zhì)均屬于極差等級,且地表水的水質(zhì)要比溫泉水更差。通過前面的分析可知,溫泉水循環(huán)深度較深,循環(huán)路徑較長,因此其水-巖相互作用時間較長,導(dǎo)致其水質(zhì)很差。地表水則由于強(qiáng)烈的蒸發(fā)濃縮作用,導(dǎo)致其離子和元素的含量很高,因此水質(zhì)很差。同時,人類活動對地表水產(chǎn)生的污染可能也是導(dǎo)致地表水F值很高、水質(zhì)較差的原因。

    4 結(jié)論

    (1)研究區(qū)不同類型的地表水和地下水的水化學(xué)類型有很大的區(qū)別。地表水的水化學(xué)類型主要為SO4-Na、Cl·SO4-Na和HCO3·SO4-Ca·Na型。溫泉水的水化學(xué)類型為SO4·HCO3-Na/HCO3·SO4-Na和HCO3·Cl-Na型。冷泉水的水化學(xué)類型為SO4·HCO3-Mg·Ca、HCO3-Ca、HCO3-Mg·Ca和SO4·HCO3-Ca型。

    (2)研究區(qū)地表水和地下水主要離子都主要來源于圍巖的溶濾或風(fēng)化,且地表水中主要離子的含量受到了強(qiáng)烈的蒸發(fā)作用的影響。研究區(qū)冷泉水中Ca2+、Mg2+、HCO3—和SO42—主要來源于白云石、方解石和石膏的溶解。溫泉水中的Na+、Ca2+、K+和HCO3—主要來源于長石類礦物的溶解。地表水中的主要離子則較復(fù)雜,具有多種來源。溫泉水和地表水中的NO—3和NH4+的相對含量主要受到了氧化還原環(huán)境的影響,而且地表水中很高的NO3—含量又受到了蒸發(fā)作用的影響。

    (3)研究區(qū)的溫泉水和冷泉水的補(bǔ)給來源都主要為大氣降水。由于水-巖相互作用,溫泉水表現(xiàn)出了明顯的“18O漂移”現(xiàn)象,而冷泉水的氫氧穩(wěn)定同位素受到了水汽再循環(huán)的影響。溫泉水和冷泉水的補(bǔ)給高程分別為3 064~4 161 m和2 874~3 943 m。

    (4)研究區(qū)的溫泉水和地表水的水質(zhì)極差,都屬于V類水,不適合直接飲用。

    猜你喜歡
    大氣研究
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    FMS與YBT相關(guān)性的實(shí)證研究
    太赫茲大氣臨邊探測儀遙感中高層大氣風(fēng)仿真
    2020年國內(nèi)翻譯研究述評
    遼代千人邑研究述論
    視錯覺在平面設(shè)計中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    新版C-NCAP側(cè)面碰撞假人損傷研究
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    一级毛片电影观看| 精品亚洲成国产av| 男人添女人高潮全过程视频| 夫妻性生交免费视频一级片| 日本与韩国留学比较| 国产成人91sexporn| 欧美日韩精品成人综合77777| 精品一区二区免费观看| 搡老乐熟女国产| 插阴视频在线观看视频| a 毛片基地| 色视频在线一区二区三区| 亚洲四区av| 久久99热6这里只有精品| 亚洲精品第二区| 国产91av在线免费观看| 色5月婷婷丁香| av.在线天堂| 中文欧美无线码| 爱豆传媒免费全集在线观看| 亚洲婷婷狠狠爱综合网| 麻豆乱淫一区二区| 精品久久久噜噜| 免费看光身美女| 亚洲欧美精品自产自拍| 色婷婷av一区二区三区视频| 免费观看无遮挡的男女| 国产无遮挡羞羞视频在线观看| 日日啪夜夜爽| 国产精品人妻久久久久久| 男的添女的下面高潮视频| 少妇丰满av| 一级毛片电影观看| 亚洲欧美中文字幕日韩二区| 午夜影院在线不卡| 99视频精品全部免费 在线| 国产亚洲5aaaaa淫片| 亚洲自偷自拍三级| 六月丁香七月| 亚洲精品国产av蜜桃| 日韩一区二区视频免费看| 色94色欧美一区二区| 男人舔奶头视频| 国产成人免费观看mmmm| 亚洲成人手机| 欧美成人精品欧美一级黄| 亚洲图色成人| 午夜福利,免费看| 黄色日韩在线| 亚洲精品,欧美精品| 日韩强制内射视频| 下体分泌物呈黄色| 精品久久久噜噜| 欧美最新免费一区二区三区| 最新中文字幕久久久久| 免费观看性生交大片5| 国产精品99久久99久久久不卡 | 91久久精品国产一区二区三区| 热re99久久国产66热| 狂野欧美白嫩少妇大欣赏| 精品人妻一区二区三区麻豆| 国产伦在线观看视频一区| 五月伊人婷婷丁香| 日韩中文字幕视频在线看片| 一本大道久久a久久精品| 新久久久久国产一级毛片| .国产精品久久| 精品国产乱码久久久久久小说| 日日爽夜夜爽网站| 乱码一卡2卡4卡精品| 99精国产麻豆久久婷婷| 免费人妻精品一区二区三区视频| 如日韩欧美国产精品一区二区三区 | 高清视频免费观看一区二区| 一区二区三区精品91| 97在线视频观看| 亚洲av免费高清在线观看| 国产成人精品一,二区| 寂寞人妻少妇视频99o| 中文精品一卡2卡3卡4更新| 性色avwww在线观看| 久久久久久久大尺度免费视频| 日韩,欧美,国产一区二区三区| 又爽又黄a免费视频| 日韩欧美一区视频在线观看 | 国产精品久久久久久av不卡| 亚洲成人一二三区av| 免费看光身美女| 男女边摸边吃奶| 777米奇影视久久| 免费大片18禁| 国产精品一二三区在线看| 国产淫语在线视频| 日本与韩国留学比较| 成人特级av手机在线观看| 欧美 亚洲 国产 日韩一| 亚洲国产毛片av蜜桃av| 久久久久久久久久久免费av| 日韩不卡一区二区三区视频在线| 人人妻人人添人人爽欧美一区卜| 91精品国产九色| 蜜桃在线观看..| 两个人的视频大全免费| 九色成人免费人妻av| 婷婷色综合www| 超碰97精品在线观看| 国产免费一区二区三区四区乱码| 中国国产av一级| 精品久久国产蜜桃| 久久精品国产自在天天线| 男女啪啪激烈高潮av片| 边亲边吃奶的免费视频| 纵有疾风起免费观看全集完整版| 交换朋友夫妻互换小说| 免费高清在线观看视频在线观看| 国产精品秋霞免费鲁丝片| 日韩三级伦理在线观看| 久久久精品94久久精品| 成人影院久久| av国产精品久久久久影院| 日本欧美国产在线视频| 乱系列少妇在线播放| 丰满少妇做爰视频| 精品一区二区三区视频在线| 国产一区有黄有色的免费视频| 午夜福利在线观看免费完整高清在| 精品久久久久久电影网| 亚洲色图综合在线观看| 日本欧美视频一区| 中文天堂在线官网| 午夜久久久在线观看| 日韩亚洲欧美综合| 在线观看人妻少妇| 高清黄色对白视频在线免费看 | 丝瓜视频免费看黄片| 免费看不卡的av| 亚洲精品乱码久久久久久按摩| 亚洲精品国产色婷婷电影| 丝袜脚勾引网站| 久久久久久久亚洲中文字幕| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| 狂野欧美激情性xxxx在线观看| 亚洲人与动物交配视频| 精品人妻熟女毛片av久久网站| 欧美亚洲 丝袜 人妻 在线| 欧美少妇被猛烈插入视频| 欧美丝袜亚洲另类| 黄色一级大片看看| 有码 亚洲区| 欧美精品亚洲一区二区| 久久av网站| 国内揄拍国产精品人妻在线| 亚洲成人手机| 免费黄频网站在线观看国产| 一级,二级,三级黄色视频| 免费久久久久久久精品成人欧美视频 | 精品酒店卫生间| 菩萨蛮人人尽说江南好唐韦庄| 2021少妇久久久久久久久久久| av在线app专区| 久久免费观看电影| 久久综合国产亚洲精品| 丝袜喷水一区| 亚洲精品乱码久久久久久按摩| 一级毛片电影观看| 中文字幕人妻熟人妻熟丝袜美| 色网站视频免费| 国内揄拍国产精品人妻在线| 亚洲国产精品成人久久小说| 国产成人午夜福利电影在线观看| 国产精品成人在线| 美女cb高潮喷水在线观看| 久久久精品免费免费高清| .国产精品久久| 日本vs欧美在线观看视频 | 亚洲av中文av极速乱| 国产高清有码在线观看视频| 免费看日本二区| 在线精品无人区一区二区三| 99久久中文字幕三级久久日本| 丰满少妇做爰视频| 美女中出高潮动态图| 久久久亚洲精品成人影院| 国产视频首页在线观看| 色视频www国产| 日本黄色日本黄色录像| 日韩伦理黄色片| 高清欧美精品videossex| 夜夜爽夜夜爽视频| 精品国产乱码久久久久久小说| 毛片一级片免费看久久久久| 老女人水多毛片| 男人舔奶头视频| 亚洲婷婷狠狠爱综合网| 午夜免费男女啪啪视频观看| 18禁裸乳无遮挡动漫免费视频| 三级国产精品片| 国产av精品麻豆| 久久ye,这里只有精品| 亚洲国产精品国产精品| 欧美精品一区二区免费开放| 黄色视频在线播放观看不卡| 欧美日韩一区二区视频在线观看视频在线| 国产av国产精品国产| 欧美一级a爱片免费观看看| 日本av免费视频播放| 午夜老司机福利剧场| 少妇的逼水好多| 国产精品福利在线免费观看| a 毛片基地| 99热这里只有是精品在线观看| 免费大片黄手机在线观看| 精品久久久久久久久亚洲| 99久久精品国产国产毛片| 久久这里有精品视频免费| 久久热精品热| 能在线免费看毛片的网站| 国产精品麻豆人妻色哟哟久久| 久久久午夜欧美精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产av成人精品| 黄色欧美视频在线观看| 天堂中文最新版在线下载| 2018国产大陆天天弄谢| 精品一区二区三区视频在线| 国产av码专区亚洲av| 日韩大片免费观看网站| 久久久久久久国产电影| 精品人妻偷拍中文字幕| 亚洲欧美日韩卡通动漫| 午夜免费观看性视频| 国产极品天堂在线| 欧美 亚洲 国产 日韩一| 国产成人freesex在线| av福利片在线| 久久99热6这里只有精品| 精品酒店卫生间| 国产国拍精品亚洲av在线观看| 噜噜噜噜噜久久久久久91| 青青草视频在线视频观看| 亚洲怡红院男人天堂| 美女主播在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产爽快片一区二区三区| 欧美变态另类bdsm刘玥| 国产av码专区亚洲av| 亚洲欧洲国产日韩| 国产伦精品一区二区三区视频9| 国产成人91sexporn| 91久久精品国产一区二区三区| 一区二区三区免费毛片| 人妻系列 视频| 成人美女网站在线观看视频| 免费大片黄手机在线观看| 欧美精品一区二区免费开放| 午夜福利影视在线免费观看| 人人妻人人爽人人添夜夜欢视频 | √禁漫天堂资源中文www| 麻豆乱淫一区二区| 三级经典国产精品| 一区二区av电影网| 91精品国产九色| 久久久久国产网址| 爱豆传媒免费全集在线观看| 婷婷色综合www| 欧美激情国产日韩精品一区| 亚洲av男天堂| 麻豆精品久久久久久蜜桃| 国产一区二区三区综合在线观看 | 国产高清不卡午夜福利| 精品国产乱码久久久久久小说| 欧美 亚洲 国产 日韩一| 日韩在线高清观看一区二区三区| 欧美日韩在线观看h| 六月丁香七月| 黄色视频在线播放观看不卡| 久久婷婷青草| 草草在线视频免费看| 老司机亚洲免费影院| 麻豆成人av视频| 七月丁香在线播放| 国产69精品久久久久777片| 免费看光身美女| 99热网站在线观看| 久热久热在线精品观看| 国产成人精品无人区| 人人妻人人添人人爽欧美一区卜| 黄色欧美视频在线观看| 26uuu在线亚洲综合色| 国产成人aa在线观看| 久久人妻熟女aⅴ| 啦啦啦视频在线资源免费观看| 亚洲精品一二三| 亚洲欧美成人综合另类久久久| 色视频www国产| 国产有黄有色有爽视频| 久久青草综合色| 国产成人精品一,二区| 少妇被粗大猛烈的视频| 亚洲欧美一区二区三区国产| 人体艺术视频欧美日本| 国产亚洲最大av| 国产成人aa在线观看| av有码第一页| 久久久精品94久久精品| 久久久a久久爽久久v久久| 久久99精品国语久久久| 欧美日韩国产mv在线观看视频| 亚洲欧美日韩东京热| 99九九线精品视频在线观看视频| 中文资源天堂在线| 我要看日韩黄色一级片| 国产成人午夜福利电影在线观看| 男女国产视频网站| 久久毛片免费看一区二区三区| 国产精品人妻久久久影院| 亚洲精品一区蜜桃| 在现免费观看毛片| 成人亚洲精品一区在线观看| 黄色怎么调成土黄色| 日日摸夜夜添夜夜添av毛片| 偷拍熟女少妇极品色| 成人午夜精彩视频在线观看| 男女国产视频网站| 18+在线观看网站| 国产精品三级大全| 两个人免费观看高清视频 | 少妇被粗大的猛进出69影院 | 男女免费视频国产| 亚洲欧美清纯卡通| 免费观看av网站的网址| 色婷婷久久久亚洲欧美| 大话2 男鬼变身卡| 久久久久网色| 国内少妇人妻偷人精品xxx网站| 老司机亚洲免费影院| 肉色欧美久久久久久久蜜桃| 春色校园在线视频观看| 日韩熟女老妇一区二区性免费视频| 又粗又硬又长又爽又黄的视频| 国产一级毛片在线| 国产成人精品久久久久久| 久久韩国三级中文字幕| 国产伦精品一区二区三区四那| 免费播放大片免费观看视频在线观看| 天堂俺去俺来也www色官网| 最近中文字幕2019免费版| 精品久久国产蜜桃| 一级毛片aaaaaa免费看小| 草草在线视频免费看| 免费在线观看成人毛片| 国产女主播在线喷水免费视频网站| 亚洲精品国产成人久久av| 午夜精品国产一区二区电影| 亚洲欧美精品自产自拍| 成人国产av品久久久| 国产高清有码在线观看视频| 最近2019中文字幕mv第一页| 亚洲三级黄色毛片| 日韩av不卡免费在线播放| 久久ye,这里只有精品| 国语对白做爰xxxⅹ性视频网站| 超碰97精品在线观看| 岛国毛片在线播放| 草草在线视频免费看| 午夜老司机福利剧场| 99re6热这里在线精品视频| 一本色道久久久久久精品综合| 另类精品久久| 夜夜看夜夜爽夜夜摸| 插阴视频在线观看视频| 亚洲欧美成人精品一区二区| 又粗又硬又长又爽又黄的视频| 九九爱精品视频在线观看| 亚洲精品乱码久久久久久按摩| 久久久久久久精品精品| 国产精品一区二区三区四区免费观看| 亚洲精品国产色婷婷电影| 免费少妇av软件| 人人妻人人澡人人看| 色视频www国产| 久久99一区二区三区| 国产高清有码在线观看视频| 中国美白少妇内射xxxbb| 中文精品一卡2卡3卡4更新| 亚洲精品国产av蜜桃| 高清黄色对白视频在线免费看 | 国产精品不卡视频一区二区| 国产成人免费观看mmmm| 国产91av在线免费观看| 美女大奶头黄色视频| 午夜91福利影院| 亚洲久久久国产精品| √禁漫天堂资源中文www| 久久久国产欧美日韩av| 国产一区二区三区av在线| 欧美成人精品欧美一级黄| 亚洲经典国产精华液单| 国产男女内射视频| 国产高清三级在线| 成人国产麻豆网| 超碰97精品在线观看| 最近中文字幕2019免费版| 最近最新中文字幕免费大全7| 亚洲精品456在线播放app| 蜜臀久久99精品久久宅男| 国产精品久久久久成人av| 国产成人精品福利久久| 久久青草综合色| 免费大片黄手机在线观看| 高清不卡的av网站| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 色94色欧美一区二区| 午夜久久久在线观看| 丁香六月天网| 自拍欧美九色日韩亚洲蝌蚪91 | 十八禁网站网址无遮挡 | 久久久久久久大尺度免费视频| 久久99精品国语久久久| 国产片特级美女逼逼视频| 插逼视频在线观看| 伦理电影免费视频| 成年av动漫网址| 精华霜和精华液先用哪个| a级毛色黄片| 美女大奶头黄色视频| 国内精品宾馆在线| 22中文网久久字幕| 男人狂女人下面高潮的视频| 亚洲综合色惰| 国产欧美日韩精品一区二区| 久久综合国产亚洲精品| 亚洲经典国产精华液单| 99久久精品国产国产毛片| 插逼视频在线观看| 王馨瑶露胸无遮挡在线观看| 成人毛片60女人毛片免费| 中文字幕精品免费在线观看视频 | 欧美成人精品欧美一级黄| 亚洲精品久久久久久婷婷小说| 老司机影院成人| 亚洲精品亚洲一区二区| 日韩 亚洲 欧美在线| 精品少妇黑人巨大在线播放| 日韩中文字幕视频在线看片| 自线自在国产av| 最黄视频免费看| 久热久热在线精品观看| 五月伊人婷婷丁香| 免费看av在线观看网站| 国产一区二区在线观看av| 久久精品国产鲁丝片午夜精品| 视频区图区小说| 欧美 日韩 精品 国产| 国产精品无大码| 天堂俺去俺来也www色官网| 69精品国产乱码久久久| 搡老乐熟女国产| 欧美 日韩 精品 国产| 亚洲国产精品一区二区三区在线| 偷拍熟女少妇极品色| 国产一区有黄有色的免费视频| 免费人妻精品一区二区三区视频| 日韩中字成人| 在现免费观看毛片| 成年人免费黄色播放视频 | 欧美亚洲 丝袜 人妻 在线| 国产亚洲最大av| av专区在线播放| 婷婷色综合大香蕉| 丝袜在线中文字幕| 午夜av观看不卡| 亚洲国产欧美日韩在线播放 | 少妇人妻精品综合一区二区| 一级毛片黄色毛片免费观看视频| 我要看黄色一级片免费的| av.在线天堂| 国产精品久久久久成人av| 国产av精品麻豆| 简卡轻食公司| 日韩av在线免费看完整版不卡| 最近2019中文字幕mv第一页| 日韩人妻高清精品专区| 女的被弄到高潮叫床怎么办| 国产极品粉嫩免费观看在线 | a 毛片基地| 我的老师免费观看完整版| 乱系列少妇在线播放| 日韩av不卡免费在线播放| 久久精品夜色国产| 边亲边吃奶的免费视频| 美女视频免费永久观看网站| 久久精品国产自在天天线| 午夜精品国产一区二区电影| 久久久久精品久久久久真实原创| 免费在线观看成人毛片| 国产精品秋霞免费鲁丝片| 又爽又黄a免费视频| 色吧在线观看| 香蕉精品网在线| 亚洲精品日韩av片在线观看| 黄色视频在线播放观看不卡| 久久精品国产亚洲av天美| 91午夜精品亚洲一区二区三区| 亚洲精品国产成人久久av| 97超碰精品成人国产| 丰满迷人的少妇在线观看| 亚洲欧美成人精品一区二区| 成人影院久久| 中国美白少妇内射xxxbb| 成人影院久久| 蜜桃在线观看..| av天堂久久9| 午夜免费观看性视频| 午夜久久久在线观看| 成人午夜精彩视频在线观看| 一级毛片电影观看| 午夜影院在线不卡| 又粗又硬又长又爽又黄的视频| 男人添女人高潮全过程视频| 国产精品免费大片| av视频免费观看在线观看| 久久ye,这里只有精品| 亚洲精品一二三| 午夜免费鲁丝| 一个人看视频在线观看www免费| av国产精品久久久久影院| 岛国毛片在线播放| 嫩草影院入口| 在线天堂最新版资源| 国产精品嫩草影院av在线观看| 久久99一区二区三区| 欧美xxxx性猛交bbbb| av在线观看视频网站免费| 超碰97精品在线观看| 最黄视频免费看| 又爽又黄a免费视频| 成人国产av品久久久| 九九爱精品视频在线观看| 国产亚洲午夜精品一区二区久久| 三级经典国产精品| 国产中年淑女户外野战色| 精品一区二区免费观看| 日韩成人伦理影院| 一本大道久久a久久精品| 欧美精品一区二区免费开放| 免费少妇av软件| 久久99热这里只频精品6学生| 高清黄色对白视频在线免费看 | 在线观看三级黄色| 九草在线视频观看| 国产 精品1| 亚洲av成人精品一二三区| 久久久欧美国产精品| 午夜91福利影院| 日本与韩国留学比较| 亚洲av成人精品一区久久| 啦啦啦中文免费视频观看日本| 交换朋友夫妻互换小说| 亚洲图色成人| 久久影院123| 久久国产精品大桥未久av | av.在线天堂| 乱系列少妇在线播放| 蜜臀久久99精品久久宅男| 三级国产精品片| 婷婷色麻豆天堂久久| 久久精品夜色国产| 在线观看免费日韩欧美大片 | 国产精品久久久久久精品古装| 国产极品粉嫩免费观看在线 | 久久人人爽人人片av| 国产精品成人在线| 亚洲人成网站在线播| 女人精品久久久久毛片| 国产淫语在线视频| 丰满迷人的少妇在线观看| 大香蕉久久网| 插阴视频在线观看视频| 欧美精品人与动牲交sv欧美| a级毛色黄片| 久久久久精品久久久久真实原创| 免费黄频网站在线观看国产| 国产亚洲午夜精品一区二区久久| 人人妻人人爽人人添夜夜欢视频 | av在线app专区| 男女边吃奶边做爰视频| 亚洲久久久国产精品| 午夜91福利影院| 成人毛片a级毛片在线播放| 免费看av在线观看网站| 九色成人免费人妻av| 国产高清不卡午夜福利| 亚洲精品中文字幕在线视频 | 少妇精品久久久久久久| 精品少妇久久久久久888优播| 伊人亚洲综合成人网| 99久久人妻综合| 18禁在线播放成人免费| 天天操日日干夜夜撸| 成年av动漫网址| 18禁在线播放成人免费| 成人毛片60女人毛片免费| 久久精品国产鲁丝片午夜精品| 综合色丁香网| 伊人亚洲综合成人网| 国产亚洲av片在线观看秒播厂| 午夜免费鲁丝| 亚洲精品自拍成人| 国产日韩欧美在线精品| 91精品伊人久久大香线蕉| 黄片无遮挡物在线观看| 久热这里只有精品99|