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

    龍門山?jīng)_斷帶盆山體系地下水循環(huán)特征分析

    2024-06-12 18:53:13鄭童心漆繼紅許模譚謠謝杭
    人民長江 2024年5期
    關(guān)鍵詞:平原區(qū)水化學(xué)淺層

    鄭童心 漆繼紅 許模 譚謠 謝杭

    摘要:山區(qū)與平原區(qū)水資源聯(lián)系密切,探究山區(qū)-盆地地下水循環(huán)可為水資源開發(fā)和合理利用提供科學(xué)依據(jù)。龍門山?jīng)_斷帶及川西前陸盆地是中國典型的盆山系統(tǒng)代表,結(jié)合現(xiàn)場采樣及前人研究結(jié)果,采用水文地質(zhì)調(diào)查、水化學(xué)特征、氫氧同位素和聚類分析等技術(shù)方法,闡明了綿竹市境內(nèi)石亭江流域龍門山?jīng)_斷推覆體與前陸盆地形成的盆山體系的地下水循環(huán)特征。研究認為:大氣降水為山區(qū)及平原區(qū)淺層地下水的主要補給來源,其中J3l3地層還接受石亭江水的側(cè)向補給;地形對山區(qū)淺層地下水徑流具有主控作用,主要沿北西、南東方向發(fā)生徑流,對近平原區(qū)的淺層地下水給予少量補給;平原區(qū)淺層還接受地表水、河渠水補給,地下水自北西向南東發(fā)生運移,在玉泉鎮(zhèn)至板橋鎮(zhèn)附近以泉的形式進行排泄。NW向斷裂構(gòu)造溝通盆山體系,來源于研究區(qū)北西側(cè)更高山區(qū)的地下水沿斷裂帶形成深部徑流帶,形成具有較高礦化度的熱水。研究成果可為類似山區(qū)-盆地地下水循環(huán)系統(tǒng)研究提供參考。

    關(guān) 鍵 詞:地下水循環(huán);水化學(xué)特征;氫氧同位素;盆山系統(tǒng);石亭江流域;龍門山?jīng)_斷帶

    中圖法分類號:P641

    文獻標(biāo)志碼:A

    DOI:10.16232/j.cnki.1001-4179.2024.05.017

    0 引 言

    在自然界的水循環(huán)過程中,一般認為山區(qū)相比較于低海拔地區(qū)會接受更多的大氣降水補給和更少的蒸散發(fā)過程,山區(qū)水資源通常是平原或盆地水資源的重要補給來源[1-2。在水資源緊缺地區(qū),山區(qū)及平原區(qū)地下水流系統(tǒng)是維持該地區(qū)經(jīng)濟發(fā)展和日常生活的重要保障。因此,探究山區(qū)-盆地地下水循環(huán)系統(tǒng)是實現(xiàn)水資源合理開發(fā)的重要前提,對該地區(qū)堅持可持續(xù)發(fā)展理念具有重要意義。

    龍門山?jīng)_斷帶及川西前陸盆地是中國中西部陸內(nèi)盆山系統(tǒng)的典型代表[3-4,其形成和發(fā)展過程及構(gòu)造背景一直備受眾多地質(zhì)學(xué)家的重視和關(guān)注5-11。由于經(jīng)歷了多次構(gòu)造變革,龍門山盆山體系中山區(qū)斷裂縱橫、構(gòu)造豐富,從而產(chǎn)生了復(fù)雜的水文地質(zhì)情況;前陸盆地為沖洪積地層,富含大量可利用的淺層地下水資源。目前,針對盆地地下水系統(tǒng)和龍門山區(qū)地下水系統(tǒng),前人已經(jīng)進行了大量的調(diào)查研究。其中,對于盆地地下水系統(tǒng)的研究主要以淺層地下水為主,具有水資源利用價值。張敏12、高東東13等通過水文地質(zhì)條件、水化學(xué)特征等對德陽市淺層地下水進行環(huán)境質(zhì)量評價及可持續(xù)發(fā)展利用研究;陳盟等[14利用礦物風(fēng)化系統(tǒng)分析、相關(guān)性分析、主成分分析和PhreeQC 反向水文地球化學(xué)模擬等方法分析了廣漢市平原區(qū)淺層地下水水化學(xué)演化及控制因素,并對飲用水質(zhì)量進行了評價;唐金平等[15-16對湔江沖洪積扇淺層地下水化學(xué)及控制因素進行分析,并利用三維數(shù)值模擬建立淺層地下水循環(huán)模式,揭示地下水循環(huán)規(guī)律。對于山區(qū)地下水系統(tǒng)的研究主要以淺層地下水和深層熱水為主,龐練17通過現(xiàn)場水文地質(zhì)調(diào)查,對成蘭鐵路龍門山段淺層巖溶水文地質(zhì)條件進行研究并對隧道進行了涌水量評價;顏玉聰18、徐艷秋19等從地質(zhì)構(gòu)造及水文地球化學(xué)等方面入手,建立溫泉水文循環(huán)模型,探究深層熱水的成因模式。徐進勇20、尚英男21、孫光旭22等依據(jù)遙感地質(zhì)解譯、水化學(xué)和同位素等方法對麓棠、廣濟-洛水地區(qū)龍門山前緣熱水形成機制進行分析,得出該地區(qū)熱水主要來源于大氣降水沿龍門山區(qū)的裂隙通道進行深循環(huán)至龍門山山前排泄。

    龍門山?jīng)_斷帶及前陸盆地地下水聯(lián)系密切,但前人對該地區(qū)地下水的研究主要集中于山區(qū)或平原區(qū),沒有將二者結(jié)合進行系統(tǒng)的山區(qū)-盆地地下水循環(huán)研究。因此,本文在前人研究的基礎(chǔ)上,以綿竹市境內(nèi)石亭江流域龍門山?jīng)_斷推覆體與前陸盆地交接地帶為例,采用水化學(xué)和氫氧同位素技術(shù),闡明該區(qū)域山區(qū)與平原區(qū)地下水和地表水循環(huán)模式,探討山區(qū)-盆地地下水循環(huán)規(guī)律,以為類似的山區(qū)-盆地地下水循環(huán)系統(tǒng)研究提供借鑒。同時,由于山區(qū)及平原區(qū)不同的地形特征,該地區(qū)地下水循環(huán)情況復(fù)雜,對該地區(qū)地下工程建設(shè)具有一定的影響,因此對該地區(qū)地下水循環(huán)的研究,可為當(dāng)?shù)氐叵滤Y源的調(diào)查評價、合理開發(fā)及地下工程建設(shè)提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    研究區(qū)位于四川盆地西北部綿竹市內(nèi),地處成都平原與其西北側(cè)山區(qū)地段的交接地帶,以龍門山四大逆沖斷裂之一的關(guān)口隱伏斷裂為分界線,北西側(cè)為山地,南東側(cè)為平原盆地[22,區(qū)內(nèi)發(fā)育梅子溝隱伏斷層,與關(guān)口隱伏斷裂正交21。研究區(qū)內(nèi)地表水系發(fā)育,其中主要河流石亭江屬沱江上游主干支流,是長江的源頭之一。河流進入平原后切割砂礫卵石層,與地下水水力聯(lián)系密切。

    研究區(qū)內(nèi)山體大致呈北東-南西走向展布,主要出露侏羅系上統(tǒng)蓮花口組地層(J3l),為一套山前沖洪積相沉積,該地層整體傾向南東,為一單斜巖層,依據(jù)巖溶水文地質(zhì)結(jié)構(gòu)的劃分依據(jù),可將研究區(qū)內(nèi)水文地質(zhì)類型歸納為單斜山型水文地質(zhì)結(jié)構(gòu)23。依據(jù)巖性可將蓮花口組地層進一步分為4段,底部蓮花口組1段(J3l1)和3段(J3l3)巖性為鈣質(zhì)礫巖;蓮花口組2段(J3l2)地層為鈣質(zhì)砂巖,為相對隔水層;頂部蓮花口組4段(J3l4)覆蓋于第四系地層之下,出露較少,為長石砂巖、粉砂巖、泥巖互層。依據(jù)研究區(qū)地層巖性、含水巖組、溶蝕裂隙發(fā)育情況以及地下水的賦存形式,可將蓮花口組1段和3段劃分為含水層,2段和4段為相對隔水層,地下水類型主要為裂隙溶隙水。

    山前平原區(qū)為第四系地層,主要賦存松散巖類孔隙水,為一套沖積、沖洪積和冰水堆積及冰漬堆積層[24,除山前小部分地區(qū)分布中、下更新統(tǒng)(Q1+2)的砂質(zhì)黏土和含砂泥礫外,大部分地區(qū)均為全新統(tǒng)沖洪積層(Q4)地層,其上部為黃色砂土、黏質(zhì)砂土,下部為砂卵礫石層。

    2 樣品采集與分析

    基于對研究區(qū)地形地貌,地下水含水層類型的認識,本次研究于2022年10月對山區(qū)、山前、沖洪積扇扇頂、扇中和扇緣方向淺層地下水及地表水進行采樣,并參考了前人在該區(qū)域?qū)ι畈裤@井S1井的研究[21,采樣點及參考點的位置見圖1。山區(qū)地下水樣采自泉點或暗河,平原區(qū)水樣采自泉點、民井及灌溉井,地表水采自水庫水和河流,現(xiàn)場采用WTW Multi3410便攜式多參數(shù)水質(zhì)檢測儀對水樣的T(水溫)、pH、Eh和TDS進行及時測定。本次研究將部分水樣送至實驗室進行分析測試,測定的參數(shù)包括Na++K+、Ca2+、Mg2+、SO42-、Cl-、HCO-3、pH、TDS和δD、δ18O值等,檢測方法依據(jù)DZ/T0064-2021《地下水質(zhì)檢標(biāo)準(zhǔn)》,Na+、K+使用差減法得出,Ca2+、Mg2+、HCO-3采用滴定法測出,SO42-和Cl-使用離子色譜法測出,TDS由計算得出,氫氧同位素使用液態(tài)水同位素分析儀測出,測定結(jié)果以相對于V-SMOW標(biāo)準(zhǔn)的千分差表示。

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

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

    研究區(qū)采樣點水化學(xué)分析結(jié)果及參考研究結(jié)果如表1所列,山區(qū)、平原區(qū)淺層地下水、地表水和S1井地下熱水pH介于6.84~8.82之間,均值7.79,其中弱堿性水占92.6%,弱酸性水占7.4%且均在平原區(qū)。山區(qū)地下水TDS為204~481 mg/L,平原區(qū)地下水TDS介于403~600 mg/L,地表水TDS介于256~424 mg/L,均屬于淡水范疇;S1井地下熱水TDS高達980.6 mg/L,接近于微咸水。依據(jù)地下水特征、存在形式和出露地層,可將水樣分為山區(qū)地下水、平原區(qū)地下水、地表水和地下熱水4種類型,山區(qū)地下水依據(jù)含水層又可分為J3l1地下水、J3l2地下水和J3l3地下水。

    根據(jù)水樣檢測結(jié)果,繪制Piper三線圖,如圖2所示。從圖2中可以看出,研究區(qū)內(nèi)從山區(qū)到平原沖洪積扇扇緣部位,地下水基本符合一般的演化規(guī)律,淺層地下水與S1井地下熱水水化學(xué)類型呈現(xiàn)明顯差異。淺層地下水及地表水陽離子以Ca2+和Mg2+為主,陰離子以HCO3-和SO42-為主,根據(jù)舒卡列夫分類方法可分為HCO3-Ca、HCO3-Ca·Mg和HCO3·SO4-Ca·Mg等類型;而S1井地下熱水以Na+和HCO3-為主,且Na+含量較淺層地下水和地表水高出10倍左右,水化學(xué)類型為HCO3-Na型,表明該地區(qū)淺層地下水與S1井地下熱水的補給和徑流屬于完全不同的兩個系統(tǒng)。S1井地下熱水產(chǎn)自1 563.5 m深處的蓮花口組礫巖中[25,可排除人為污染的影響因素,Cl-濃度高達40.78 mg/L,水中的高氯化物是地下水通過深部循環(huán)發(fā)生水巖作用導(dǎo)致的[21,所以S1井地下熱水的水化學(xué)特征可反應(yīng)該地區(qū)深部地下水特征。

    平原區(qū)淺層地下水、地表水和山區(qū)的淺層地下水水化學(xué)具有差異性。山區(qū)地下水除個別泉點(S12、S13、S14)外均落在圖中①區(qū),屬于HCO3-Ca和HCO3-Ca·Mg型水。而平原區(qū)地下水、地表水和泉點(S12、S13、S14)均落在圖中②區(qū),水化學(xué)類型以HCO3·SO4-Ca·Mg為主,與山區(qū)地下水相比其差異主要體現(xiàn)在SO42-離子上。Piper三線圖的主菱形圖上顯示研究區(qū)水樣以堿土金屬和弱酸根離子為主[26,在陽離子三角圖中,山區(qū)地下水和平原區(qū)地下水離子的分布不存在明顯的差異,大都集中在Ca2+區(qū)域;在陰離子三角圖中,平原區(qū)地下水的SO42-占比落在18%~40%之間,山區(qū)地下水的SO42-占比均小于18%,進一步說明陰離子的濃度差異是影響山區(qū)和平原區(qū)地下水類型分異的關(guān)鍵因素。山區(qū)泉點S12、S13、S14為三箭水泉及其附近泉水,與山區(qū)其他泉點相比其補給來源除接受大氣降水外,還接受石亭江的側(cè)向補給,從而發(fā)生混合作用,其水化學(xué)成分不能代表山區(qū)淺層地下水化學(xué)特征。在Piper三線圖中,點S12、S13、S14落在②區(qū),更靠近平原區(qū)地表水和地下水,說明其接受石亭江地表水的補給量大于山區(qū)地下水的補給量。

    為了進一步探討山區(qū)淺層地下水、平原區(qū)淺層地下水和地表水的水化學(xué)特征和循環(huán)關(guān)系,對研究區(qū)內(nèi)所取得的26組水樣進行聚類分析,用來突出不同水樣點之間的空間關(guān)系,將樣本按照它們在性質(zhì)上的相似程度進行分類[27-29。本次研究選取26組水樣樣本,7個變量,分別為TDS、Na++K+、Ca2+、Mg2+、Cl-、SO42-和HCO3-。聚類分析前通過標(biāo)準(zhǔn)化使得變量的平均值為0,標(biāo)準(zhǔn)離差為1。本次聚類方法選取Ward(離差平方和)法,平方歐式距離,將樣本進行逐一歸類,具體計算由軟件SPSS 21.0進行,聚類結(jié)果見圖3和表2。由圖3可以看出,26組水樣分類效果較好,當(dāng)重心距離為5時,水樣可分為4類,同時根據(jù)各組內(nèi)水體樣品離子含量的平均值繪制Stiff圖。

    Ⅰ類水包括8個山區(qū)泉點,1個山區(qū)地表水,石亭江下游水和困牛山水庫水。其中,6個出露于J3l1和J3l3的鈣質(zhì)礫巖段泉點均屬于此類,地下水直接接受大氣降水補給,由于鈣質(zhì)礫巖段易發(fā)生溶蝕,落水洞、洼地較為發(fā)育,補給條件較好,地下水流動較快,且受地表污染較少,J3l1和J3l3的泉點具有礦化度較低、各種離子含量相差較小等特點。2個出露于J3l2的泉水也屬于此類,主要由于地下水切層流動,徑流途徑短,水中離子含量較低,所以屬于此類。山區(qū)地表水補給來源與泉點相似,因此具有相似的水化學(xué)性質(zhì),也都屬于Ⅰ類。困牛山水庫攔蓄地表水,其主要來源為石亭江江水,礦化度和各種離子含量也相對較低,都屬于Ⅰ類。

    Ⅱ類水包括5個山區(qū)泉點和1個山區(qū)地表水。這些泉點均出露于J3l2的鈣質(zhì)砂巖段內(nèi),相較于1段和3段地下水,由于這類地下水順層流動,徑流途徑長或在含水層中滯留時間長,使得水中各種離子含量偏高,且這些泉點都分布于云蓋村附近,補徑排途徑相似,因此具有相似的性質(zhì)。

    Ⅲ類水主要包括石亭江上游和下游水,困牛山水庫的引水渠水,S12三箭水泉和井水S26。此類水水化學(xué)成分相似,均為HCO3·SO4-Ca·Mg型。其中引水渠水來自于石亭江上游,與石亭江上游水具有相同的水化學(xué)性質(zhì)。三箭水泉流量較大,其主要補給來源大部分為石亭江引水,少部分來自于山區(qū)地下水,因此三箭水泉與石亭江上游水具有相似的水化學(xué)性質(zhì)。井水S26位于山區(qū)與平原的交界地帶,且接近困牛山水庫,而困牛山水庫和三箭水泉的主要補給來源均為石亭江水,因此S26可能接受困牛山水庫的滲漏補給和三箭水泉的補給,因此與三箭水泉同屬于一類。

    Ⅳ類水均為平原區(qū)井水和泉點,分布于沖洪積扇扇中和扇緣部位,其地下水水化學(xué)性質(zhì)與前三類水有很大區(qū)別。平原區(qū)地下水與山區(qū)地下水相比較,具有較高的TDS和離子含量,Na++K+含量達到22.7~52.4 mg/L,Cl-也高達14.1 mg/L,SO42-明顯高于山區(qū)地下水,與地表水相當(dāng)。除這些易受人為污染干擾的離子外,其Ca2+和Mg2+含量也與山區(qū)地下水不同。

    根據(jù)聚類結(jié)果及上述分析,研究區(qū)內(nèi)山區(qū)淺層地下水與部分平原區(qū)淺層地下水具有水力聯(lián)系,山區(qū)地下水又可分為J3l1和J3l3的鈣質(zhì)礫巖段含水層和J3l2的鈣質(zhì)砂巖段含水層。山區(qū)淺層地下水主要接受大氣降水的補給,J3l1和J3l3的鈣質(zhì)礫巖含水系統(tǒng)水化學(xué)特征反映的是補給徑流條件好,水循環(huán)較快;大部分J3l2的鈣質(zhì)砂巖段含水系統(tǒng)水化學(xué)特征反映的是地下水順層流動,徑流途徑較長,徑流環(huán)境較為開放,多為HCO-3型水。平原區(qū)地下水系統(tǒng)主要接受降雨入滲補給、河、渠水入滲補給和少量的山區(qū)地下水側(cè)向補給,徑流途徑長,并且存在人為作用對水化學(xué)的影響。

    3.2 氫氧同位素特征

    由于降水在蒸發(fā)和凝結(jié)過程中會發(fā)生同位素分餾,使得大氣降水的氫氧同位素之間存在線性關(guān)系。受當(dāng)?shù)亟邓亩握舭l(fā)、水蒸汽氣團的起源等影響,不同地區(qū)的大氣降水線會存在一定的差異。胡月等[30通過對成都地區(qū)113場降水樣品氫氧同位素的分析,得出成都的大氣降水線方程(CDLMWL):δD=7.83δ18O+8.79。

    山區(qū)地下水δ18O變化范圍為-9.98‰~-8.59‰,平均值為-9.39‰,δD變化范圍為-65.1‰~-58.25‰,平均值為-61.66‰;平原區(qū)地下水δ18O變化范圍為-9.81‰~-7.54‰,平均值為-8.63‰,δD變化范圍為-66.11‰~-49.59‰,平均值為 -57.06‰。根據(jù)研究區(qū)氫氧同位素分析測試及前人研究結(jié)果,繪制δD-δ18O關(guān)系圖(圖4)。從圖4中可以看出,研究區(qū)內(nèi)大部分水樣點分布在成都市大氣降水線的左上方,表明研究區(qū)地下水來源均為大氣降水且受蒸發(fā)作用影響微弱[31。S1井地下熱水位于關(guān)系圖的左下方,靠近大氣降水線,表明其來源主要為大氣降水,且與其他水體相比,其δD值要負得多。根據(jù)氫氧同位素高程效應(yīng)原理,地下水的δD值隨著補給高程的增加而呈現(xiàn)出更負的趨勢,因此可以確定S1井地下熱水的補給來源應(yīng)該是研究區(qū)北西側(cè)更高的山區(qū)[20。

    依據(jù)在研究區(qū)采集到的不同水體的氫氧同位素,分別對山區(qū)及平原區(qū)的δD和δ18O進行線性擬合,得到山區(qū)線性方程為δD=5.7δ18O-8.13(R2=0.87),平原區(qū)線性方程為δD=7.29δ18O+5.83(R2=0.99),與大氣降水線分布接近,再次表明研究區(qū)地下水均來源于大氣降水補給[32。大部分山區(qū)淺層地下水的樣品在δD-δ18O關(guān)系圖中同位素值接近,表明山區(qū)地下水有相似的補給來源[33。S12三箭水泉位于δD-δ18O 關(guān)系圖的左下方,呈現(xiàn)與其他山區(qū)地下水泉點相比較為偏負的現(xiàn)象,更接近于石亭江的δD、δ18O值,且S12三箭水泉和石亭江水水化學(xué)類型均為HCO3·SO4-Ca·Mg型,說明S12三箭水泉及其附近地下水大部分補給來源為石亭江的側(cè)向補給。平原區(qū)地下水在氫氧同位素散點圖上分布比較分散,山前中部平原由北向南地下水的同位素含量差別較大,表明平原區(qū)淺層地下水具有多種補給途徑,主要表現(xiàn)為北側(cè)井點S24、S25的δD、δ18O同位素值較大,且水化學(xué)類型均為HCO3-Ca·Na型;南側(cè)井點S23的δD、δ18O同位素值較小,其水化學(xué)類型為HCO3-Ca· Mg。結(jié)合水文地質(zhì)條件分析,從山區(qū)地下水、困牛山水庫水到平原區(qū)北側(cè)井點S24、S25,δ18O值越來越高,表明其補給高程相對較低,水庫水和近處高程較低的山區(qū)地下水存在共同補給S24、S25的可能,同時結(jié)合水化學(xué)成分及聚類結(jié)果分析,S24和S25的Na+含量較高,推測該點δD、δ18O值也可能受到人為因素影響。南側(cè)井點S23在δD-δ18O關(guān)系圖中更靠近石亭江水,推測其主要補給來源可能為石亭江的側(cè)向補給。

    山區(qū)淺層地下水氘盈余值介于10.46‰~14.75‰,平均值為13.46‰,平原區(qū)淺層地下水氘盈余值介于10.69‰~12.56‰,平均值為11.94‰(圖5),與全國氘盈余(d=8%~12%)[34-35接近,各泉點氘盈余變化較小,可能是由于地下水補給源構(gòu)成穩(wěn)定。大氣降水補給進入到地下后,隨著水巖作用的發(fā)生,同位素將發(fā)生交換。一般情況下,水巖作用不會影響水體的δD值,但可以導(dǎo)致水體相對富含δ18O,從而使d值降低,且水巖作用越強烈,地下水的d值就越小36。研究區(qū)山區(qū)淺層地下水氘盈余值略大于平原區(qū)淺層地下水,不排除山區(qū)地下水及石亭江水對平原區(qū)補給作用過程中發(fā)生水巖作用導(dǎo)致d值減小的可能。

    3.3 山區(qū)和平原區(qū)地下水循環(huán)模式

    綜合以上分析,研究區(qū)地下水循環(huán)模式為:大氣降水及地表水是研究區(qū)山區(qū)及平原區(qū)淺層地下水的主要補給來源,且蒸發(fā)作用微弱。山區(qū)淺層地下水受地形地貌、地層巖性、含水巖組空間組合特征的影響,在接受大氣降水補給后,分別向北東、南西沿垂向裂隙及層間溶隙發(fā)生徑流。其中J3l2段地下水在接受大氣降水補給后,少部分切層流動補給J3l3段地下水,大部分向北東、南西順層徑流形成較高TDS的地下水;J3l3段在接受大氣降水補給的同時,還接受石亭江的側(cè)向補給,沿垂向裂隙及層間溶隙發(fā)生徑流。山區(qū)深層地下水補給源為研究區(qū)北西側(cè)更高的山區(qū),在接受大氣降水后沿構(gòu)造帶循環(huán)至深部,經(jīng)熱儲層再經(jīng)機井開掘,形成溫泉井。平原區(qū)淺層地下水除接受大氣降水外,還接受地表水、渠水、和北西側(cè)山區(qū)裂隙溶隙水的少量補給,地下水自北西向南東運移形成地下徑流帶,在玉泉鎮(zhèn)至板橋鎮(zhèn)附近的溢出帶以泉的形式進行排泄(圖6)。

    4 結(jié) 論

    (1)基于水化學(xué)成分及聚類結(jié)果對研究區(qū)不同水體地下水進行分析,深層地下水和淺層地下水水化學(xué)組成呈現(xiàn)出顯著的差異,淺層地下水表現(xiàn)為北西側(cè)山區(qū)至南東側(cè)平原,地下水由HCO3-Ca和HCO3-Ca·Mg 型過渡為HCO3·SO4-Ca·Mg型;深層地下水主要為HCO3-Na型。三箭水泉(S12)及其附近泉點水化學(xué)類型呈現(xiàn)與石亭江地表水水化學(xué)類型相似,表明石亭江對三箭水泉具有一定的側(cè)向補給。

    (2)山區(qū)和平原區(qū)淺層地下水的氫氧同位素組成表明地下水接受大氣降水的補給,蒸發(fā)作用微弱;深層地下水主要為北西側(cè)龍門山脈更高的山區(qū)接受大氣降水并沿斷裂帶入滲循環(huán)至此。山區(qū)淺層地下水氫氧同位素值較為集中,其補給來源相似,其中三箭水泉的δD、δ18O值與其他山區(qū)淺層地下水相比,更靠近石亭江水,再次證明石亭江水對山區(qū)地下水存在側(cè)向補給;平原區(qū)淺層地下水δD、δ18O值較為分散,推測由于主要補給來源不同造成,證明平原區(qū)淺層地下水具有多種補給來源。同時山區(qū)與平原區(qū)的氘盈余值證明山區(qū)地下水對平原區(qū)地下水具有補給作用。

    (3)研究區(qū)山區(qū)地下水與平原區(qū)地下水補給過程受地質(zhì)條件、地形地貌等多種因素控制,僅通過水文地球化學(xué)及同位素技術(shù)難以定量化識別地下水循環(huán)模式,后期可結(jié)合數(shù)值模擬等方法建立概念模型,可進一步深入認識地下水循環(huán)規(guī)律。

    參考文獻:

    [1]董艷輝,符韻梅,王禮恒,等.甘肅北山-河西走廊-祁連山區(qū)域地下水循環(huán)模式[J].地質(zhì)科技通報,2022,41(1) :79-89.

    [2]梁杏,張婧瑋,藍坤,等.江漢平原地下水化學(xué)特征及水流系統(tǒng)分析[J].地質(zhì)科技通報,2020,39(1) :21-33.

    [3]李智武,宋天慧,王自劍,等.川西-龍門山盆山系統(tǒng)走向差異演化的變形、隆升和沉積記錄及關(guān)鍵構(gòu)造變革期討論[J].成都理工大學(xué)學(xué)報(自然科學(xué)版),2021,48(3):257-282.

    [4]劉樹根,羅志立,趙錫奎,等.中國西部盆山系統(tǒng)的耦合關(guān)系及其動力學(xué)模式:以龍門山造山帶-川西前陸盆地系統(tǒng)為例[J].地質(zhì)學(xué)報,2003,77(2):177-186.

    [5]YAN D P,ZHOU Y,QIU L,et al.The Longmenshan tectonic complex and adjacent tectonic units in the eastern margin of the Tibetan Plateau:a review[J].Journal of Asian Earth Sciences,2018,164:33-57.

    [6]JIN W Z,TANG L J,YANG K M,et al.Segmentation of the Longmen Mountains thrust belt,western Sichuan foreland basin,SW China[J].Tectonophysics,2010,485(1-4):107-121.

    [7]JIN W Z,TANG L J,YANG K M,et al.Tectonic evolution of the middle frontal area of the Longmen Mountain thrust belt,western Sichuan basin,China[J].Acta Geologica Sinica,2009,83(1):92-102.

    [8]LIAO T P,HU J J,ZHANG F R,et al.Relation between structural evolution of the Longmenshan orogenic zone and sedimentation of its foreland basin[J].Mining Science and Technology,2009,19(6):807-812.

    [9]李智武,陳洪德,劉樹根,等.龍門山?jīng)_斷隆升及其走向差異的裂變徑跡證據(jù)[J].地質(zhì)科學(xué),2010,45(4):944-968.

    [10]王二七,孟慶任,陳智樑,等.龍門山斷裂帶印支期左旋走滑運動及其大地構(gòu)造成因[J].地學(xué)前緣,2001,8(2):375-384.

    [11]賈東,陳竹新,賈承造,等.龍門山前陸褶皺沖斷帶構(gòu)造解析與川西前陸盆地的發(fā)育[J].高校地質(zhì)學(xué)報,2003,9(3):402-410.

    [12]張敏.德陽市平原區(qū)淺層地下水環(huán)境質(zhì)量研究[D].成都:成都理工大學(xué),2014.

    [13]高東東,吳勇,陳盟.德陽市地下水富集規(guī)律及可持續(xù)開發(fā)利用[J].科學(xué)技術(shù)與工程,2018,18(4):8-14.

    [14]陳盟,吳勇,高東東,等.廣漢市平原區(qū)淺層地下水化學(xué)演化及其控制因素[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2016,46(3):831-843.

    [15]唐金平,張強,胡漾,等.湔江沖洪積扇地下水化學(xué)特征及控制因素分析[J].環(huán)境科學(xué),2019,40(7):3089-3098.

    [16]唐金平.湔江沖洪積扇地下水化學(xué)特征與淺層地下水循環(huán)研究[D].成都:成都理工大學(xué),2020.

    [17]龐練.成蘭鐵路龍門山段水文地質(zhì)條件及其對隧道工程的影響研究:以鄧家坪隧道為例[D].成都:成都理工大學(xué),2010.

    [18]顏玉聰,劉峰立,郭麗爽,等.龍門山斷裂帶溫泉水文地球化學(xué)特征[J].地震研究,2021,44(2):170-184.

    [19]徐艷秋.四川盆地西部邊緣地下熱水水化學(xué)特征及成因研究[D].北京:中國地質(zhì)大學(xué)(北京),2021.

    [20]徐進勇,王多義,尹觀.綿竹三箭水S1井溫泉補給源的同位素示蹤研究[J].礦物巖石,2003,2(23):75-78.

    [21]尚英男.綿竹S1井及其鄰區(qū)溫泉地下水同位素及元素地球化學(xué)特征[D].成都:成都理工大學(xué),2003.

    [22]孫光旭.廣濟-洛水地區(qū)龍門山前緣溫泉形成機制及勘探前景研究[D].成都:成都理工大學(xué),2013.

    [23]李亞津,巫錫勇,蒲大攀,等.西南山區(qū)巖溶水文地質(zhì)結(jié)構(gòu)類型及地下水動力特征分析[J].路基工程,2020(6):13-18,24.

    [24]孫為奕.綿竹市平原區(qū)地下水脆弱性評價[D].成都:成都理工大學(xué),2017.

    [25]王多義.綿竹釀春池溫泉坡縷石的特征及成因探討[J].礦物巖石,2004,24(2):14-18.

    [26]黃靖宇,許模,許漢華,等.滇中高原巖溶關(guān)鍵帶典型含水系統(tǒng)水化學(xué)特征及成因分析[J].地質(zhì)科技通報,2022,41(5):347-356.

    [27]CLOUTIER V,LEFEBVRE R,THERRIEN R,et al.Multivariate statistical analysis of geochemical data as indicative of the hydrogeochemical evolution of groundwater in a sedimentary rock aquifer system[J].Journal of Hydrology,2008,353(3-4):294-313.

    [28]李舒,楊佳雪,李小倩,等.地下水化學(xué)組成的時空聚類分析與多級嵌套水流系統(tǒng)識別[J].地質(zhì)科技通報,2022,41(1):309-318.

    [29]王東輝,吳銘杰,李勝偉.岷江流域含石膏地層地下水水化學(xué)差異性研究[J].人民長江,2023,54(7):14-21.

    [30]胡月,劉國東,孟玉川,等.成都次降水穩(wěn)定氫氧同位素特征及水汽來源分析[J].環(huán)境科學(xué),2019,40(3):1179-1187.

    [31]PANG Z H,KONG Y L,LI J,et al.An isotopic geoindicator in the hydrological cycle[J].Procedia Earth and Planetary Science,2017,17:534-537.

    [32]蒲俊兵.重慶巖溶地下水氫氧穩(wěn)定同位素地球化學(xué)特征[J].地球?qū)W報,2013(6):713-722.

    [33]張世殊,劉希,施裕兵.四川省自一里水電站廠址區(qū)地下水來源識別[J].人民長江,2021,52(6):71-76.

    [34]宋洋.塔里木河流域東部降水氫氧穩(wěn)定同位素特征及其影響因素[D].蘭州:西北師范大學(xué),2021.

    [35]顧慰祖.同位素水文學(xué)[M].北京:科學(xué)出版社,2011.

    [36]尹觀,倪師軍,高志友,等.四川盆地鹵水同位素組成及氘過量參數(shù)演化規(guī)律[J].礦物巖石,2008,28(2):56-62.

    (編輯:劉 媛)

    Characteristics of groundwater circulation in basin-mountain system of Longmenshan

    alluvial belt:case of Mianzhu section of Shiting River Basin

    ZHENG Tongxin,QI Jihong,XU Mo,TAN Yao,XIE Hang

    (State Key Laboratory of Geohazard Prevention and Geoenvironment Protection,Chengdu University of Technology,Chengdu 610059,China)

    Abstract:The water resources of mountainous and plain areas are closely connected,and the investigation on the mountain-basin groundwater cycle can provide a scientific basis for the development and rational utilization of water resources.The Longmenshan alluvial belt and the western Sichuan foreland basin are typical representatives of the basin-mountain system in China.Combined with the field sampling and the results of previous studies,the groundwater circulation characteristics of the basin-mountain system formed by the Longmenshan alluvial overburden and the foreland basin in the Shiting River Basin in the territory of Mianzhu section were clarified by using hydrogeological investigation,hydrochemical characterization,hydrogen and oxygen isotopes,and cluster analysis methods.The study concludes that atmospheric precipitation is the main recharge source of shallow groundwater in the mountainous areas and plains,in which the J3l3stratum also receives lateral recharge from the Shiting River.Topography has a dominant role in controlling the shallow groundwater runoff in the mountainous areas,which mainly occurs along the north-west and south-east directions,and has a small amount of recharge of the shallow groundwater in the near-plain areas.The shallow layer in the plains also receives recharge from surface water and river and canal water,and the groundwater moves from northwest to southeast,finally it discharges in the form of springs along Yuquan Town to Banqiao Town.As the NW fracture structure connects with the basin-mountain system,the groundwater from the higher mountainous areas in the NW side of the study area forms a deep runoff zone along the fracture zone,forming hot water with high mineralization.The study results can be used as a reference for similar studies of mountain-basin groundwater circulation systems.

    Key words:groundwater circulation;hydrochemical characteristics;hydrogen and oxygen isotopes;basin-mountain system;Shiting River Basin;Longmenshan alluvial belt

    猜你喜歡
    平原區(qū)水化學(xué)淺層
    幸福湖泉域地下水水化學(xué)特征及其演化機制
    人民黃河(2025年1期)2025-01-16 00:00:00
    淺層換填技術(shù)在深厚軟土路基中的應(yīng)用
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評價
    基于淺層曝氣原理的好氧顆粒污泥的快速培養(yǎng)
    淺層地下水超采區(qū)劃分探究
    河北省平原區(qū)新近系熱儲回灌的可行性與前景分析
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    河北省平原區(qū)館陶組熱儲地下熱水動態(tài)特征
    保定市平原區(qū)淺層地下水水質(zhì)變化趨勢
    華北地質(zhì)(2015年3期)2015-12-04 06:13:29
    亚洲男人天堂网一区| 久久精品国产综合久久久| 一区二区三区激情视频| 欧美黄色淫秽网站| 久久精品成人免费网站| 啦啦啦免费观看视频1| 一级a爱视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 久久亚洲真实| 亚洲av美国av| 日韩欧美国产一区二区入口| 热re99久久国产66热| 免费在线观看黄色视频的| 成人影院久久| 男女下面插进去视频免费观看| 欧美av亚洲av综合av国产av| 欧美精品亚洲一区二区| 久久热在线av| 交换朋友夫妻互换小说| 午夜a级毛片| 亚洲熟女毛片儿| 在线观看舔阴道视频| 精品久久久久久电影网| 国产黄色免费在线视频| 88av欧美| 婷婷精品国产亚洲av在线| 一本综合久久免费| 法律面前人人平等表现在哪些方面| 国产精品美女特级片免费视频播放器 | 久久青草综合色| 亚洲男人的天堂狠狠| 久久中文字幕一级| 国产亚洲精品第一综合不卡| 亚洲一区中文字幕在线| av电影中文网址| 欧美日韩国产mv在线观看视频| 最近最新中文字幕大全电影3 | 三级毛片av免费| av福利片在线| 两性夫妻黄色片| 在线观看午夜福利视频| 一级作爱视频免费观看| 午夜福利在线观看吧| 免费高清在线观看日韩| 亚洲人成77777在线视频| 日韩三级视频一区二区三区| 精品乱码久久久久久99久播| 国产高清激情床上av| 国产精品香港三级国产av潘金莲| 在线观看午夜福利视频| 精品久久久久久久毛片微露脸| 日韩免费av在线播放| 精品久久蜜臀av无| 露出奶头的视频| 亚洲欧美精品综合一区二区三区| cao死你这个sao货| 成人免费观看视频高清| 国产精品国产av在线观看| 欧美日韩黄片免| 999精品在线视频| 亚洲五月色婷婷综合| 欧美成人午夜精品| 日日夜夜操网爽| 中文字幕高清在线视频| 一级毛片高清免费大全| 中文字幕av电影在线播放| 久久青草综合色| av国产精品久久久久影院| 熟女少妇亚洲综合色aaa.| 1024视频免费在线观看| 嫁个100分男人电影在线观看| 少妇裸体淫交视频免费看高清 | bbb黄色大片| 日本a在线网址| 啦啦啦 在线观看视频| 搡老岳熟女国产| 麻豆久久精品国产亚洲av | 黄色视频不卡| 国产精品国产av在线观看| 波多野结衣av一区二区av| 天堂俺去俺来也www色官网| 欧美大码av| 国产精品秋霞免费鲁丝片| 色精品久久人妻99蜜桃| 91精品三级在线观看| 三上悠亚av全集在线观看| 国产精品1区2区在线观看.| 无限看片的www在线观看| 色婷婷久久久亚洲欧美| 亚洲一区中文字幕在线| 在线观看免费视频日本深夜| 97人妻天天添夜夜摸| 一二三四在线观看免费中文在| 亚洲欧美日韩另类电影网站| 每晚都被弄得嗷嗷叫到高潮| 99国产极品粉嫩在线观看| 女同久久另类99精品国产91| 亚洲欧洲精品一区二区精品久久久| 精品欧美一区二区三区在线| 久久草成人影院| 成年女人毛片免费观看观看9| 最新在线观看一区二区三区| 久久精品91无色码中文字幕| 久久亚洲真实| 日韩人妻精品一区2区三区| 免费人成视频x8x8入口观看| 中文字幕人妻丝袜制服| 性色av乱码一区二区三区2| 亚洲精品中文字幕一二三四区| 久久久久久久精品吃奶| av超薄肉色丝袜交足视频| a级毛片在线看网站| 在线观看午夜福利视频| 香蕉国产在线看| 一区在线观看完整版| 在线免费观看的www视频| 欧美亚洲日本最大视频资源| 亚洲一区二区三区不卡视频| 少妇的丰满在线观看| 国产成人精品久久二区二区91| 另类亚洲欧美激情| 欧美在线黄色| 999精品在线视频| 九色亚洲精品在线播放| 99国产精品一区二区三区| 老司机在亚洲福利影院| 精品一区二区三区视频在线观看免费 | 男女高潮啪啪啪动态图| 黑人猛操日本美女一级片| 亚洲 欧美一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 女性被躁到高潮视频| 欧美日韩精品网址| 在线观看舔阴道视频| 中文字幕精品免费在线观看视频| 色综合婷婷激情| 丁香六月欧美| 大码成人一级视频| 波多野结衣一区麻豆| 欧美在线一区亚洲| 亚洲一码二码三码区别大吗| 亚洲av第一区精品v没综合| 免费观看人在逋| 免费人成视频x8x8入口观看| www.999成人在线观看| 日韩欧美一区二区三区在线观看| 久久国产乱子伦精品免费另类| 每晚都被弄得嗷嗷叫到高潮| 热re99久久精品国产66热6| 色播在线永久视频| 午夜成年电影在线免费观看| 国产精品一区二区在线不卡| 50天的宝宝边吃奶边哭怎么回事| 咕卡用的链子| 黄色成人免费大全| 欧美+亚洲+日韩+国产| 波多野结衣av一区二区av| 亚洲,欧美精品.| 国内毛片毛片毛片毛片毛片| 又大又爽又粗| av福利片在线| 一级片'在线观看视频| 别揉我奶头~嗯~啊~动态视频| 久久亚洲精品不卡| 亚洲成av片中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放 | 大香蕉久久成人网| 狂野欧美激情性xxxx| av国产精品久久久久影院| 午夜激情av网站| 神马国产精品三级电影在线观看 | 久久久久久免费高清国产稀缺| 岛国在线观看网站| 97人妻天天添夜夜摸| av免费在线观看网站| 看片在线看免费视频| av网站在线播放免费| 国产精品久久视频播放| 1024视频免费在线观看| 国产精品久久久久久人妻精品电影| 黄色视频不卡| 黄色视频不卡| 亚洲熟妇中文字幕五十中出 | 国产精品亚洲av一区麻豆| 精品午夜福利视频在线观看一区| 欧美丝袜亚洲另类 | 日韩有码中文字幕| 欧美中文日本在线观看视频| 中文字幕最新亚洲高清| 欧美乱色亚洲激情| 国产99白浆流出| 国产一区二区三区视频了| 国产男靠女视频免费网站| 国产主播在线观看一区二区| 亚洲熟妇熟女久久| 亚洲五月天丁香| 欧美成狂野欧美在线观看| 99香蕉大伊视频| 男人操女人黄网站| 国产精品亚洲av一区麻豆| 黑人欧美特级aaaaaa片| 国产亚洲精品综合一区在线观看 | 一级作爱视频免费观看| 欧美成狂野欧美在线观看| 欧美成狂野欧美在线观看| 亚洲色图 男人天堂 中文字幕| 99在线视频只有这里精品首页| tocl精华| 久久香蕉激情| 精品一区二区三区四区五区乱码| 久久精品成人免费网站| a在线观看视频网站| 欧美精品一区二区免费开放| 午夜精品在线福利| 国产极品粉嫩免费观看在线| 久久久精品国产亚洲av高清涩受| 在线观看日韩欧美| 如日韩欧美国产精品一区二区三区| 国产精品偷伦视频观看了| 后天国语完整版免费观看| 欧美日韩视频精品一区| 成人影院久久| 99国产精品99久久久久| 在线观看舔阴道视频| 9191精品国产免费久久| 757午夜福利合集在线观看| 久久精品91蜜桃| 无遮挡黄片免费观看| 欧美老熟妇乱子伦牲交| 男女做爰动态图高潮gif福利片 | 免费在线观看日本一区| 9色porny在线观看| netflix在线观看网站| 男人舔女人的私密视频| 欧美av亚洲av综合av国产av| 免费av毛片视频| 99在线人妻在线中文字幕| 天堂中文最新版在线下载| 99精品久久久久人妻精品| 国产精品美女特级片免费视频播放器 | 久久精品亚洲av国产电影网| 在线观看日韩欧美| 日本vs欧美在线观看视频| av网站在线播放免费| 搡老乐熟女国产| 可以在线观看毛片的网站| 好看av亚洲va欧美ⅴa在| 又紧又爽又黄一区二区| 久久亚洲精品不卡| 欧美日韩乱码在线| 午夜亚洲福利在线播放| 亚洲精华国产精华精| 制服诱惑二区| 精品一区二区三区四区五区乱码| 免费久久久久久久精品成人欧美视频| 91大片在线观看| 成人精品一区二区免费| 免费日韩欧美在线观看| 日本黄色视频三级网站网址| 韩国精品一区二区三区| 国产一区在线观看成人免费| 午夜福利,免费看| 99精国产麻豆久久婷婷| 久久性视频一级片| 亚洲成a人片在线一区二区| 50天的宝宝边吃奶边哭怎么回事| a级毛片黄视频| 女人高潮潮喷娇喘18禁视频| 岛国视频午夜一区免费看| 成人免费观看视频高清| 国产精品久久久久久人妻精品电影| 精品国产一区二区久久| av国产精品久久久久影院| 国产1区2区3区精品| 老司机深夜福利视频在线观看| 国产乱人伦免费视频| 欧美成人午夜精品| 神马国产精品三级电影在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 成人av一区二区三区在线看| 亚洲自拍偷在线| 亚洲第一青青草原| 久久久久久大精品| 最近最新免费中文字幕在线| 夜夜夜夜夜久久久久| 精品日产1卡2卡| 国产成人免费无遮挡视频| 成年版毛片免费区| 在线观看一区二区三区激情| 久久中文字幕人妻熟女| 国产国语露脸激情在线看| 少妇的丰满在线观看| 最新美女视频免费是黄的| 韩国av一区二区三区四区| 亚洲精品成人av观看孕妇| 久久久国产欧美日韩av| 99在线人妻在线中文字幕| 电影成人av| aaaaa片日本免费| 91老司机精品| 久久久久国产精品人妻aⅴ院| 麻豆成人av在线观看| 亚洲欧美日韩无卡精品| 在线看a的网站| 国产亚洲欧美98| 丝袜美足系列| 亚洲精品中文字幕在线视频| 美女大奶头视频| 午夜福利在线免费观看网站| 久久香蕉激情| 欧美大码av| 欧美日韩国产mv在线观看视频| 三上悠亚av全集在线观看| 女性生殖器流出的白浆| 操出白浆在线播放| 亚洲情色 制服丝袜| 一边摸一边抽搐一进一出视频| 亚洲一区二区三区色噜噜 | 无遮挡黄片免费观看| 亚洲九九香蕉| 精品欧美一区二区三区在线| 淫妇啪啪啪对白视频| 色老头精品视频在线观看| 成人国产一区最新在线观看| 1024香蕉在线观看| а√天堂www在线а√下载| 香蕉国产在线看| 操美女的视频在线观看| 黄片播放在线免费| 成年版毛片免费区| 国产男靠女视频免费网站| www.999成人在线观看| 啪啪无遮挡十八禁网站| 两人在一起打扑克的视频| 怎么达到女性高潮| 变态另类成人亚洲欧美熟女 | 一区二区三区激情视频| 亚洲第一青青草原| e午夜精品久久久久久久| 国产精品98久久久久久宅男小说| 亚洲精品国产精品久久久不卡| 国产在线精品亚洲第一网站| 午夜激情av网站| 琪琪午夜伦伦电影理论片6080| 91麻豆av在线| 老司机靠b影院| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩黄片免| 性色av乱码一区二区三区2| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区蜜桃| 天堂中文最新版在线下载| 国产激情欧美一区二区| 麻豆成人av在线观看| 国产免费现黄频在线看| 国产精品久久电影中文字幕| 在线观看www视频免费| 国产精品av久久久久免费| 免费av中文字幕在线| 国产高清视频在线播放一区| 亚洲中文av在线| 国产精品亚洲av一区麻豆| 国产精品国产高清国产av| 国产主播在线观看一区二区| 又大又爽又粗| 99精品在免费线老司机午夜| 亚洲 国产 在线| 欧美日韩乱码在线| 亚洲精品国产一区二区精华液| 九色亚洲精品在线播放| 亚洲片人在线观看| 亚洲一区二区三区欧美精品| 中文字幕最新亚洲高清| 每晚都被弄得嗷嗷叫到高潮| 淫秽高清视频在线观看| 久久久国产成人精品二区 | 一个人免费在线观看的高清视频| 亚洲欧美激情综合另类| 欧美乱色亚洲激情| 久久精品aⅴ一区二区三区四区| 亚洲一区二区三区欧美精品| 欧美久久黑人一区二区| 男人舔女人的私密视频| 欧美午夜高清在线| 免费搜索国产男女视频| 久久国产精品人妻蜜桃| 三上悠亚av全集在线观看| 精品国产乱码久久久久久男人| 国产精品久久久久久人妻精品电影| 国产精品永久免费网站| 亚洲欧美精品综合一区二区三区| 岛国视频午夜一区免费看| 亚洲成a人片在线一区二区| a级毛片黄视频| 1024香蕉在线观看| 欧美人与性动交α欧美精品济南到| 9色porny在线观看| 女人高潮潮喷娇喘18禁视频| 国产成人精品无人区| 淫妇啪啪啪对白视频| 日本五十路高清| 欧美在线一区亚洲| 精品国产乱子伦一区二区三区| 伊人久久大香线蕉亚洲五| 免费观看人在逋| av天堂久久9| 日韩免费av在线播放| 日日干狠狠操夜夜爽| 淫妇啪啪啪对白视频| 99久久人妻综合| 一级片免费观看大全| 久久人妻av系列| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 我的亚洲天堂| 91成人精品电影| 大香蕉久久成人网| 亚洲国产欧美网| 亚洲国产看品久久| 欧美色视频一区免费| 好看av亚洲va欧美ⅴa在| avwww免费| a在线观看视频网站| 看免费av毛片| 国产真人三级小视频在线观看| 91国产中文字幕| 国产1区2区3区精品| 精品国产一区二区三区四区第35| 欧美乱妇无乱码| 精品卡一卡二卡四卡免费| av视频免费观看在线观看| 国产熟女午夜一区二区三区| 精品久久久精品久久久| 久久中文看片网| 国产麻豆69| 少妇 在线观看| 亚洲性夜色夜夜综合| www.精华液| 欧美+亚洲+日韩+国产| 亚洲一码二码三码区别大吗| 伦理电影免费视频| 精品人妻1区二区| 日韩大码丰满熟妇| 丝袜美足系列| 黄片小视频在线播放| 黄色片一级片一级黄色片| 757午夜福利合集在线观看| 亚洲av成人av| 性欧美人与动物交配| 一边摸一边抽搐一进一出视频| 亚洲精品国产精品久久久不卡| 成人亚洲精品av一区二区 | 一夜夜www| 亚洲在线自拍视频| 亚洲精品在线观看二区| 亚洲五月婷婷丁香| av电影中文网址| 91精品三级在线观看| 精品欧美一区二区三区在线| 身体一侧抽搐| 国产伦人伦偷精品视频| 看免费av毛片| 巨乳人妻的诱惑在线观看| 高清毛片免费观看视频网站 | 天堂俺去俺来也www色官网| 99精国产麻豆久久婷婷| 日韩视频一区二区在线观看| 午夜免费成人在线视频| 亚洲精华国产精华精| 久久伊人香网站| 久久欧美精品欧美久久欧美| 黄片大片在线免费观看| 亚洲av片天天在线观看| 亚洲第一青青草原| 亚洲第一欧美日韩一区二区三区| 亚洲 欧美 日韩 在线 免费| 亚洲国产精品一区二区三区在线| 亚洲 欧美一区二区三区| 丁香欧美五月| 欧美日韩福利视频一区二区| 亚洲专区国产一区二区| 一个人免费在线观看的高清视频| 两个人看的免费小视频| 午夜福利欧美成人| 亚洲人成电影免费在线| 变态另类成人亚洲欧美熟女 | 精品一区二区三卡| 国产极品粉嫩免费观看在线| 9热在线视频观看99| 亚洲性夜色夜夜综合| 两个人看的免费小视频| 亚洲片人在线观看| 99久久精品国产亚洲精品| 久久人人97超碰香蕉20202| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美激情在线| 在线观看午夜福利视频| 99国产精品一区二区三区| 久久中文字幕一级| 91九色精品人成在线观看| 淫秽高清视频在线观看| 高潮久久久久久久久久久不卡| av欧美777| www国产在线视频色| 999久久久国产精品视频| 午夜福利在线观看吧| 在线十欧美十亚洲十日本专区| 色综合站精品国产| 国产日韩一区二区三区精品不卡| 美女午夜性视频免费| 精品久久久久久电影网| 日韩中文字幕欧美一区二区| 在线看a的网站| 一本综合久久免费| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品1区2区在线观看.| 99精国产麻豆久久婷婷| 日日爽夜夜爽网站| 欧美老熟妇乱子伦牲交| 一级片'在线观看视频| 亚洲伊人色综图| 男人操女人黄网站| 欧美成人午夜精品| 亚洲一区二区三区欧美精品| 久久精品亚洲熟妇少妇任你| 啦啦啦 在线观看视频| 50天的宝宝边吃奶边哭怎么回事| 日韩欧美一区二区三区在线观看| 欧美丝袜亚洲另类 | 亚洲熟女毛片儿| 视频在线观看一区二区三区| 欧美人与性动交α欧美软件| 90打野战视频偷拍视频| 日韩大尺度精品在线看网址 | 国产精品香港三级国产av潘金莲| 首页视频小说图片口味搜索| 999精品在线视频| 青草久久国产| 在线观看66精品国产| 日韩人妻精品一区2区三区| 淫妇啪啪啪对白视频| 欧美老熟妇乱子伦牲交| 精品日产1卡2卡| 女人高潮潮喷娇喘18禁视频| 美女 人体艺术 gogo| 18禁观看日本| 欧美色视频一区免费| 国产无遮挡羞羞视频在线观看| 成人黄色视频免费在线看| 超碰97精品在线观看| 高清毛片免费观看视频网站 | 亚洲国产欧美日韩在线播放| 日韩高清综合在线| 九色亚洲精品在线播放| 免费人成视频x8x8入口观看| 亚洲自偷自拍图片 自拍| 欧美黑人欧美精品刺激| 黄片播放在线免费| 日本黄色日本黄色录像| 脱女人内裤的视频| av有码第一页| 免费av中文字幕在线| 91老司机精品| 日韩精品免费视频一区二区三区| 中出人妻视频一区二区| 精品一品国产午夜福利视频| 一级片'在线观看视频| 国产精品影院久久| 亚洲人成伊人成综合网2020| 精品人妻1区二区| 亚洲avbb在线观看| 免费在线观看完整版高清| 一级毛片高清免费大全| 性色av乱码一区二区三区2| 成人永久免费在线观看视频| 一个人免费在线观看的高清视频| 久久青草综合色| 中文字幕人妻丝袜一区二区| 高潮久久久久久久久久久不卡| 亚洲 国产 在线| 他把我摸到了高潮在线观看| 人妻久久中文字幕网| 在线观看一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 神马国产精品三级电影在线观看 | 国产xxxxx性猛交| 国产成人欧美在线观看| 国产野战对白在线观看| e午夜精品久久久久久久| 欧美日本亚洲视频在线播放| 极品人妻少妇av视频| 成人国产一区最新在线观看| 日本一区二区免费在线视频| 欧美日韩亚洲高清精品| 后天国语完整版免费观看| 天堂俺去俺来也www色官网| 日韩大码丰满熟妇| 久久香蕉精品热| 大香蕉久久成人网| 亚洲精品在线美女| 十八禁人妻一区二区| 国产亚洲精品久久久久久毛片| 最近最新免费中文字幕在线| 九色亚洲精品在线播放| a级毛片在线看网站| 一级a爱视频在线免费观看| 亚洲国产欧美一区二区综合| 久久香蕉精品热| 真人一进一出gif抽搐免费| 色哟哟哟哟哟哟| 99香蕉大伊视频| 欧美日韩亚洲高清精品| 99热国产这里只有精品6| 亚洲专区中文字幕在线| 99香蕉大伊视频|