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

    四川鮮水河-安寧河斷裂帶溫泉氫氧穩(wěn)定同位素特征*

    2021-04-17 01:31:16張磊郭麗爽劉樹文楊耀施得旸
    巖石學(xué)報 2021年2期
    關(guān)鍵詞:鮮水河氫氧康定

    張磊 郭麗爽 劉樹文 楊耀 施得旸

    1. 應(yīng)急管理部國家自然災(zāi)害防治研究院,北京 1000852. 中國地震局地殼動力學(xué)重點實驗室,北京 1000853. 北京大學(xué)地球與空間科學(xué)學(xué)院,北京 1008714. 四川省地震局,成都 6100415. 中國地震局地球物理研究所,北京 100081

    地殼表層水循環(huán)系統(tǒng)中的不同類型水體,如大氣降水、地表水和地下水具有不同的同位素組成特征(Clark and Fritz,1997)。2H和18O作為水分子的組成部分,它們的穩(wěn)定同位素是研究水的來源和水文循環(huán)過程的最有效的示蹤工具(Craig,1961;Gat,1996)。溫泉作為一種天然環(huán)保地?zé)豳Y源,是鏈接大氣降水、地表水、地下水和地下深部流體的水熱載體,是地球水循環(huán)過程中的重要組成部分, 因此, 對溫泉開展氫氧穩(wěn)定同位素研究,能夠有效識別溫泉的流體來源、水巖作用與水文循環(huán)過程等(周訓(xùn)等,2015;Lietal.,2018),為溫泉資源的開發(fā)與利用提供重要依據(jù)。

    四川西部高原地區(qū)分布著大量隆起山地型地?zé)豳Y源,主要受斷裂帶分布控制(王貴玲,2018)。尤其是位于川西青藏高原東緣大型活動斷裂帶-鮮水河和安寧河斷裂帶上分布著大量中高溫溫泉。李曉等(2018)通過水化學(xué)和氫氧穩(wěn)定同位素手段研究了鮮水河斷裂帶康定至道孚段的水熱活動區(qū)熱水的形成和聯(lián)系,提出了不同水源補(bǔ)給的熱水系統(tǒng)形成模式。Guoetal.(2017a)分析了康定高溫地?zé)崴乃瘜W(xué)組成、氣體組分和同位素組成,識別了影響熱水組成的水文地球化學(xué)過程(包括混合和脫氣作用)以及儲層溫度。Luoetal.(2017)利用地?zé)釡貥?biāo)、氫氧穩(wěn)定同位素和水文地球化學(xué)模擬等手段評價了康定南北部地?zé)衢_發(fā)潛力和優(yōu)先順序。也有學(xué)者通過對青藏高原東緣溫泉的氦同位素、重力、磁和地震數(shù)據(jù),獲得了巴塘到康定區(qū)域的地下熱結(jié)構(gòu)和熱液活動與地?zé)釀恿W(xué)的關(guān)系(Tangetal.,2017)。所有這些研究為探討青藏高原東緣溫泉的熱源和形成模式提供了依據(jù)。

    鮮水河-安寧河斷裂帶及其周邊大地震頻發(fā)(Wenetal.,2008),如發(fā)生在鮮水河斷裂帶上的1955年M7.5、1973年M7.6、1981年M6.9和2014年M6.3級地震。地震與構(gòu)造活動可引起泉水物理和水化學(xué)性質(zhì)發(fā)生變化,反之?dāng)嗔褞У叵滤绕涫菧厝刑N(yùn)含豐富的地下構(gòu)造活動信息(趙永紅等,2011;周曉成等,2020),對其深入研究有助于評價構(gòu)造活動對地下水資源的影響并可能獲得可靠的地震信息(Skeltonetal.,2014,2019)。Chenetal.(2014)和Lietal.(2019)等發(fā)現(xiàn)康定區(qū)域溫泉在大地震前后出現(xiàn)明顯的氫氧穩(wěn)定同位素變化,揭示了構(gòu)造活動引起了斷裂帶含水層滲透性的增加或減小。氫氧穩(wěn)定同位素技術(shù)為研究斷裂帶水循環(huán)過程及其與構(gòu)造活動之間的關(guān)系提供了證據(jù)。

    前人對青藏高原東緣鮮水河-安寧河斷裂帶溫泉的氫氧穩(wěn)定同位素研究多集中在鮮水河斷裂帶上康定區(qū)域,對鮮水河-安寧河斷裂帶上分布的眾多溫泉及地表水的氫氧穩(wěn)定同位素的系統(tǒng)對比研究還較缺乏。在解決鮮水河-安寧河斷裂帶不同位置溫泉的氫氧穩(wěn)定同位素分布差異原因、區(qū)域地表水對地?zé)崴难a(bǔ)給等方面仍需深入開展相關(guān)工作,這對于深入了解活動斷裂帶條帶分布地下水的循環(huán)特征有重要意義。為此本文以鮮水河-安寧河斷裂帶上的溫泉和地表水作為一個系統(tǒng)研究對象,開展氫氧穩(wěn)定同位素對比分析,同時結(jié)合溫泉水中化學(xué)組分特征,初步揭示活動斷裂帶上地?zé)崴耐凰氐厍蚧瘜W(xué)及水文循環(huán)特征。這對于研究青藏高原東緣斷裂帶上地?zé)崴乃难h(huán)過程及其與構(gòu)造活動的關(guān)系具有重要的意義。

    圖1 研究區(qū)域位置與采樣點分布示意圖Fig.1 The study area and the sampling locations

    1 研究區(qū)域概況

    研究區(qū)域(27.7°~31.2°N、101.1°~102.4°E,)主要位于四川西部道孚-康定-石棉-西昌一帶,地處青藏高原東緣和四川盆地西緣(圖1)。氣候類型主要為高原型季風(fēng)氣候和山地氣候,降雨集中在每年5月至9月。鮮水河斷裂帶位于四川甘孜至石棉一帶,為北西向的左旋走滑斷裂帶。安寧河斷裂帶位于四川石棉至西昌一帶,以左旋走滑為主,近南北走向,長約160km。

    強(qiáng)烈的構(gòu)造活動以及高山深切峽谷地形有利于溫泉的形成,研究區(qū)域大多數(shù)溫泉沿斷裂呈帶狀分布(羅來麟,1994)。鮮水河斷裂帶上的溫泉主要集中在道孚、中谷、康定和瀘定磨西等4個熱水區(qū)(羅來麟,1994),多為高中溫溫泉,屬于川西高原高中溫地?zé)釁^(qū),熱儲巖性主要為花崗巖、石灰?guī)r、大理巖、變質(zhì)砂板巖以及第四系松散層膠結(jié)等(王貴玲,2018)。鮮水河斷裂帶地下冷水主要是基巖裂隙水,其補(bǔ)給為大氣降水和高山冰雪融水(趙慶生,1984)。安寧河斷裂帶溫泉主要分布在喜德和西昌周邊,屬于川西南中低溫地?zé)釁^(qū),熱儲巖性一般為碳酸鹽巖(羅來麟,1994;王貴玲,2018)。

    2 樣品采集與測試

    本研究于2019年5月共采集了鮮水河-安寧河斷裂帶上38個樣品,包括溫泉水樣品24個,溫泉附近河水12個,冷泉水樣品1個,雪線附近積雪融水1個,采樣位置見圖1和表1。采樣點命名中,S代表溫泉水,S后數(shù)字代表采樣點編號,R代表溫泉附近的河流,R后數(shù)字為相應(yīng)位置的溫泉編號,CS1為冷泉水,M1為積雪融水(圖1、表1)。在溫泉主出水口處采集水樣,共采集了中高溫溫泉(60℃≤T<90℃)4處、中溫溫泉(40℃≤T<60℃)16處和低溫溫泉(25℃≤T<40℃)4處。采集溫泉附近的流動河水樣品,同時采集了道孚七美一處山頂雪線處積雪融水樣品,康定一處冷泉樣品。采樣現(xiàn)場水樣使用0.45μm過濾膜過濾,水樣裝滿后密封保存在高密度聚乙烯塑料瓶內(nèi)。在采樣現(xiàn)場,使用AZ8821型號數(shù)字溫度計測量水溫。

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

    3.1 氫氧穩(wěn)定同位素組成分布特征

    研究區(qū)域溫泉的海拔高度為1175~3946m,水溫變化范圍為37.29~74.81℃,水溫平均值為52.02℃(表1)。溫泉最高水溫為康定灌頂溫泉(74.81℃),最低水溫為西昌河西溫泉(37.29℃)。研究區(qū)域河水的海拔高度為1734~4197m,河水水溫變化范圍為6.48~23.82℃,水溫平均值為12.66℃(表1),呈現(xiàn)出海拔高水溫低、海拔低水溫高的特征,這與采樣點處的氣溫有關(guān)系。

    研究區(qū)溫泉、河水和積雪融水的δ18O和δ2H顯示出差異性的分布特征(圖2)。溫泉δ18O變化范圍為-19.04‰~-12.71‰,平均值為-16.42‰;δ2H變化范圍為-144.07‰~-88.63‰,平均值為-122.37‰(n=24)。在康定采集的1處冷泉(T=19.88℃),δ18O和δ2H分別為-13.66‰和-106.74‰。研究區(qū)域河水的δ18O變化范圍為-15.90‰~-10.85‰,平均值為-13.86‰;δ2H變化范圍為-118.21‰~-71.12‰,平均值為-98.99‰(n=12)。雪線處積雪融水的δ18O和δ2H分別為-10.27‰和-65.41‰,相比溫泉和河水樣品富集重同位素。

    通過水樣品的δ18O和δ2H與全球和區(qū)域大氣降水線對比,能夠了解水的來源與補(bǔ)給、水巖作用以及地下水混合作用等水文循環(huán)特征(Pangetal.,2017)。全球大氣降水線方程(Global Meteoric Water Line,GMWL)為δ2H=8δ18O+10(Craig,1961)。研究區(qū)域位于青藏高原東緣,采用青藏高原大氣降水線為區(qū)域大氣降水線(Local Meteoric Water Line,LMWL),降水線方程為δ2H=8.41δ18O+16.72(Kongetal.,2019)。δ18O和δ2H對比結(jié)果(圖2)顯示出研究區(qū)域的大部分溫泉和河水的δ18O和δ2H靠近GMWL和LMWL,說明補(bǔ)給源均為大氣降水成因。但是不同類型采樣點上同位素組成具有顯著的差異性,如24個溫泉樣品具有4個明顯的分布區(qū)域,河水R1~R17與R18~R23同樣具有顯著的同位素組成差異。上述氫氧穩(wěn)定同位素的離散分布特征,表明了大氣降水補(bǔ)給源的差別,這可能與不同位置大氣降水同位素分布差異因素有關(guān)。冷泉水同位素呈現(xiàn)出向右側(cè)偏離大氣降水線的特征(圖2)。冷泉水的氫氧穩(wěn)定同位素代表當(dāng)?shù)卮髿饨邓钠骄凰亟M成(王基華等,2000),其向右偏離大氣降水線的分布特征表明了補(bǔ)給過程中的蒸發(fā)作用引起了同位素的富集。

    表1 采樣點位置和氫氧穩(wěn)定同位素組成

    圖2 水樣δ18O-δ2H分布及其與全球和區(qū)域大氣降水線的關(guān)系Fig.2 Plot of δ18O vs. δ2H and their correlations with GMWL and LMWL

    溫泉與其附近的河流同位素值存在較大差異,溫泉比其附近河流δ18O和δ2H值低(圖2)。這種特征同樣出現(xiàn)在西藏地?zé)崽镏校缪虬司責(zé)崽?鄭淑蕙等,1982)。上述分析結(jié)果表明了近距離的河流和周邊大氣降水不是溫泉的主要補(bǔ)給源,溫泉的補(bǔ)給源為遠(yuǎn)距離的大氣降水下滲經(jīng)地下裂隙和管道運(yùn)移補(bǔ)給。但是S17溫泉與附近河流(R17)同位素值接近(圖2),表明了該溫泉的補(bǔ)給來源于附近的大氣降水。這與該溫泉為貢嘎山環(huán)繞的小區(qū)域氣候環(huán)境有直接關(guān)系。

    3.2 同位素高程效應(yīng)

    研究區(qū)域位于青藏高原東緣的不同位置,采樣點之間存在著明顯的海拔高度差。前人研究表明,高程效應(yīng)是影響大氣降水同位素的重要因素,表現(xiàn)出隨著海拔升高同位素值下降的特征(Clark and Fritz,1997)。高程效應(yīng)能夠區(qū)分不同海拔處降水對含水層的補(bǔ)給(Blasch and Bryson,2007),例如Songetal.(2006)通過懷沙河流域內(nèi)泉水的氧同位素高程效應(yīng)確定了地下水的補(bǔ)給范圍和水流路徑。溫泉的δ18O和δ2H與泉水出露高程之間存在著顯著的負(fù)相關(guān)關(guān)系(圖3a, b)。

    圖3 溫泉(a、b)和河水(c、d)采樣點海拔與δ18O和δ2H的關(guān)系Fig.3 Plot of altitude vs. δ18O and δ2H for hot springs (a, b) and for rivers (c, d)

    表2 溫泉和冷泉水化學(xué)組成結(jié)果(mg/L)

    雖然溫泉均是上升泉,但是其補(bǔ)給為區(qū)域內(nèi)大氣降水(圖2),而且沿著青藏高原東緣鮮水河-安寧河斷裂帶分布的研究區(qū)域溫泉高程變化明顯,因此溫泉的同位素高程效應(yīng)可能與大氣降水的同位素高程效應(yīng)有關(guān)。

    研究區(qū)域溫泉同位素分布特征受高程效應(yīng)控制,溫泉海拔與同位素值具有強(qiáng)負(fù)相關(guān)性(圖3),同位素高程效應(yīng)擬合公式為:

    δ18O=-0.0023h-10.011

    (1)

    δ2H=-0.0195h-67.813

    (2)

    式中,h為溫泉泉口出露高程(m)。因此得到氧同位素高程效應(yīng)為-0.23‰/100m,氫同位素高程效應(yīng)為-1.95‰/100m,即海拔每增加100m,δ18O值減少0.23‰,δ2H值減少1.95‰。Yuetal.(1984)通過對貴州-重慶-四川-西藏一線的水同位素研究,提出了大氣降水氧同位素高程效應(yīng)為-0.31‰/100m,氫同位素高程效應(yīng)為-2.6‰/100m,大于本文研究區(qū)域同位素效應(yīng),這與其研究區(qū)域海拔高度差(4980m)明顯大于本研究(2770m)有直接關(guān)系。本研究的同位素高程效應(yīng)更適用于鮮水河-安寧河斷裂帶上的四川道孚-康定-石棉-西昌區(qū)域。溫泉水同位素與溫泉高程的強(qiáng)相關(guān)性(圖3),表明了引起研究區(qū)域溫泉水同位素離散分布的主要因素為高程變化。河水的δ18O與海拔之間的R2值表明了斷裂帶河水存在著高程效應(yīng)(圖3c, d),但是河水的高程效應(yīng)R2值明顯小于溫泉,這可能與溫泉的補(bǔ)給和水汽來源地穩(wěn)定有關(guān),而河水的補(bǔ)給存在著以大氣降水為主的多種來源。

    3.3 溫泉氧同位素漂移

    地?zé)崴憩F(xiàn)出兩種典型的氧同位素分布特征,一類是氧同位素漂移并表現(xiàn)出水巖作用特征。由于補(bǔ)給熱儲的大氣降水與巖石之間同位素的不平衡,會發(fā)生水與巖石礦物之間氧同位素的交換,導(dǎo)致地?zé)崴赝凰馗患跬凰叵虼髿饨邓€右側(cè)漂移(顧慰祖等,2011)。氧同位素漂移可作為地?zé)嵯到y(tǒng)深部溫度的判定指標(biāo)(衛(wèi)克勤等,1983)。如西藏羊八井地?zé)崴跬凰仄x大氣降水線2‰,云南熱海地?zé)崽锏钠?‰~3‰,屬于與火山和巖漿有關(guān)的地?zé)崽镱愋?衛(wèi)克勤等,1983;Guoetal.,2017b)。另一種氧同位素靠近大氣降水線,未發(fā)生明顯的氧同位素漂移。本文報道的溫泉氫氧穩(wěn)定同位素靠近大氣降水線的分布特征,反映出溫泉的補(bǔ)給與循環(huán)速度快,并未發(fā)生明顯的水和巖石氧同位素交換作用(除S15康定龍頭溝溫泉)。S15溫泉泉口出露溫度72.17℃,地下水循環(huán)過程中與巖石發(fā)生了氧同位素交換作用,顯示出高溫作用下的地?zé)崴囊欢ǔ潭人畮r作用的特征。

    圖4 泉水樣品Na-K-Mg三角圖Fig.4 Na-K-Mg triangular diagram for springs

    3.4 水化學(xué)特征及水巖作用

    在Na-K-Mg三角圖上(圖4),大部分溫泉樣品和冷泉樣品落在Mg端元附近,表明了地下水的水巖作用尚未達(dá)到離子平衡狀態(tài),為循環(huán)相對較快的淺層冷水補(bǔ)給。S9、S21和S24溫泉處在部分平衡水區(qū)域,但是靠近Mg端元,顯示出了一定程度的水巖作用特征。地?zé)崴械奈⒘吭嘏c深部流體活動密切相關(guān),為深入研究影響流體組成的各種地質(zhì)過程提供了依據(jù)(Kaasalainenetal.,2015)。S12和S13溫泉與CS1冷泉同處康定同一區(qū)域,相隔距離近,溫泉Li和Sr含量略大于冷泉含量(表2),顯示出溫泉水巖作用弱的特征。來自氫氧穩(wěn)定同位素和水化學(xué)組成的結(jié)果,說明鮮水河-安寧河斷裂帶上的大部分溫泉以大氣降水補(bǔ)給為主且循環(huán)相對較快,小部分溫泉存在著一定程度的水巖反應(yīng)。

    3.5 溫泉循環(huán)深度

    同處康定區(qū)域的溫泉S9~S16與冷泉CS1有著明顯的同位素分布差異(圖2),顯示出溫泉同位素值較冷泉值低的特征。這種現(xiàn)象與我國東部郯廬斷裂帶和膠遼斷塊區(qū)溫泉與冷泉的同位素分布特征相同,即地?zé)崴摩?8O和δ2H值低于當(dāng)?shù)乩淙南鄳?yīng)同位素值,同位素差異與地下水深循環(huán)過程產(chǎn)生的同位素物理分餾有關(guān)(Shangguanetal.,1998)。Shangguanetal.(1998)發(fā)現(xiàn)地?zé)崴?H和δ18O的變化幅度與循環(huán)深度之間存在著強(qiáng)相關(guān)性,并能用來評價地?zé)崴难h(huán)深度,計算方法如下:

    d(δ18O)=-1.470Δδ18O+3.162

    (3)

    d(δ2H)=-0.397Δδ2H+1.064

    (4)

    式中,d為地?zé)崴畲笱h(huán)深度(km),Δδ18O為地?zé)崴c冷泉δ18O的差,Δδ2H為地?zé)崴c冷泉δ2H的差。公式(3)和(4)計算康定地區(qū)8口溫泉的d(δ18O)和d(δ2H)平均值估算得到溫泉循環(huán)深度為3.9~10.2km。使用上述方法估算的溫泉循環(huán)深度與δ18O和δ2H具有強(qiáng)相關(guān)性(圖5),表明隨著溫泉的循環(huán)深度增加,同位素值下降。由于研究區(qū)域多為高山峽谷區(qū),溫泉存在著強(qiáng)的高程效應(yīng)(圖3),溫泉同位素與循環(huán)深度存在強(qiáng)相關(guān)性(圖5),分析溫泉與冷泉同位素的差異性可能與補(bǔ)給區(qū)高程以及地下水深淺循環(huán)過程中的同位素分餾有關(guān)。

    圖5 溫泉循環(huán)深度與δ18O和δ2H的對比關(guān)系Fig.5 Plot of circulation depths vs. δ18O and δ2H for hot springs

    除上述氫氧穩(wěn)定同位素方法計算循環(huán)深度,還可以應(yīng)用地?zé)釡貥?biāo)計算溫泉的循環(huán)深度。常用的SiO2地?zé)釡貥?biāo)和陽離子地?zé)釡貥?biāo)等方法利用經(jīng)驗公式估算熱儲溫度(王瑩等,2007)。王逸凌等(2020)對鮮水河斷裂帶的溫泉采用了SiO2地?zé)釡貥?biāo)、陽離子地?zé)釡貥?biāo)以及多礦物平衡溫標(biāo)綜合得出了鮮水河斷裂帶不同區(qū)域的熱儲溫度,與本文研究康定區(qū)域類似范圍內(nèi)的熱儲溫度為180~ 250℃。如果溫泉的溫度受其循環(huán)所處的地?zé)峥刂?,可通過下式估算熱水循環(huán)深度(李娟等,2007):

    Z=G(TZ-T0)+Z0

    (5)

    式中Z為循環(huán)深度(m),G為地?zé)嵩鰷芈?,?0.4m/℃(張云輝,2018),TZ為熱儲溫度(℃),T0為補(bǔ)給區(qū)年平均氣溫,取7℃,Z0為年常溫帶深度,取30m。根據(jù)此式推算出康定區(qū)域溫泉的循環(huán)深度在3.6~5.0km。利用氫氧穩(wěn)定同位素計算的循環(huán)深度較利用熱儲溫度方法計算的值偏大,可能與同位素的高程效應(yīng)導(dǎo)致同位素值偏低有關(guān)。S9溫泉的Cl-和Li和Si含量較低(表2),也說明其循環(huán)深度可能比氫氧穩(wěn)定同位素計算的10.2km要小。該溫泉的海拔明顯較其他康定區(qū)域溫泉高,可能同位素高程效應(yīng)引起該溫泉與淺循環(huán)的冷泉氫氧同位素差值較大,進(jìn)而得到的循環(huán)深度偏高。對同一小范圍內(nèi)的溫泉S12、S13和冷泉CS1進(jìn)行計算得到的循環(huán)深度(圖5)也較利用公式(5)計算的偏高,可能與冷泉的蒸發(fā)效應(yīng)引起的同位素值偏大,引起Δδ18O和Δδ2H變化范圍偏大有關(guān)。因此,在利用δ2H和δ18O的變化幅度評價高山區(qū)域地?zé)崴难h(huán)深度時,應(yīng)考慮同位素高程效應(yīng)以及冷泉的蒸發(fā)效應(yīng)的影響??紤]到由不同經(jīng)驗公式估算循環(huán)深度可能存在不確定性,而且不同種方法之間存在著誤差,康定區(qū)域溫泉的循環(huán)深度可能由公式(5)估算的更加適用。

    3.6 溫泉成因與同位素監(jiān)測意義

    鮮水河斷裂帶為下切到地幔的深大活動斷裂帶,溫泉逸出氣體中氦同位素和碳同位素結(jié)果顯示出幔源流體作用(Zhouetal.,2015),但是氫氧穩(wěn)定同位素未顯示出明顯的向左或向右偏離大氣降水線的水巖作用和二氧化碳?xì)怏w作用特征(圖2),而且水化學(xué)揭示的大部分溫泉水巖作用特征不強(qiáng)烈。 這可能與補(bǔ)給份額以大氣降水補(bǔ)給為主有關(guān),地下水循環(huán)速度快、滯留時間短,而且深部補(bǔ)給份額少不足以引起水的同位素變化。溫泉和地表水氫氧穩(wěn)定同位素的分布特征,表明了溫泉的形成過程,即遠(yuǎn)距離補(bǔ)給的大氣降水和地表水通過地下裂隙運(yùn)移補(bǔ)給地下水,地下水在地下熱儲層加熱后通過斷裂通道上升到地表形成溫泉,這符合王貴玲和藺文靜(2020)提出的川西溫泉主要為經(jīng)典的對流型水熱系統(tǒng)成因類型特征。Tangetal.(2017)和Zhangetal.(2017)對川西巴塘-理塘-康定高溫水熱活動區(qū)進(jìn)行了地球物理(地震、航磁和重力)和氦同位素的觀測,結(jié)果顯示青藏高原東緣地?zé)嵯到y(tǒng)的成因為大氣降水和地表水沿斷裂帶滲入和深循環(huán),被地殼熱源加熱后經(jīng)循環(huán)上升到地表形成溫泉,熱源來自于逃逸的幔源熱流、花崗巖產(chǎn)生的放射性熱和斷層摩擦熱。Shietal.(2017)對川西溫泉的水化學(xué)聚類分析研究證明了康定區(qū)域溫泉水化學(xué)成因來源于淺層水混和。氫氧穩(wěn)定同位素結(jié)果與地球物理和水化學(xué)統(tǒng)計分析結(jié)果是一致的,而且氫氧穩(wěn)定同位素更易識別補(bǔ)給源的遠(yuǎn)近及其與地表水的補(bǔ)給關(guān)系。

    活動斷裂帶溫泉氫氧穩(wěn)定同位素監(jiān)測是研究構(gòu)造活動特征的重要手段之一。首先,氫氧穩(wěn)定同位素揭示出了溫泉的補(bǔ)給、循環(huán)與水巖作用特征,能夠評價是否適宜進(jìn)行地震觀測(Zhangetal.,2020)。其次,深循環(huán)地下水的同位素變化是穩(wěn)定的(Clark and Fritz,1997),近距離地表水的干擾較小,這在本研究溫泉和河流同位素區(qū)別分布上得到證實。該性質(zhì)對于排除小區(qū)域大氣降水和河流補(bǔ)給等環(huán)境干擾對溫泉氫氧穩(wěn)定同位素的影響十分重要,有利于識別與構(gòu)造活動有關(guān)的異常。再次,活動地震斷裂帶上溫泉的氫氧穩(wěn)定同位素值在強(qiáng)震活動前后會發(fā)生明顯的變化(Skeltonetal.,2014;Zhangetal.,2020)。已有的報道顯示出2008年汶川地震后,鮮水河-安寧河斷裂帶上部分溫泉出現(xiàn)明顯的同位素變化,比如位于康定區(qū)域的二道橋溫泉(S14)、龍頭溝溫泉(S15)和灌頂溫泉(S16)發(fā)生了明顯的氧同位素偏移(Chenetal.,2014;Lietal.,2019)。同位素變化與地殼應(yīng)力變化引起含水層性質(zhì)變化有關(guān),如含水層滲透性發(fā)生變化,致使不同類型的地下水混合或者補(bǔ)給比率發(fā)生變化,或者前震與同震促使含水層發(fā)生破裂,水與新鮮裂隙面接觸增加了水巖反應(yīng)程度(Thomas,1988;Claessonetal.,2007)。那么,定期開展活動斷裂帶上多口溫泉的氫氧穩(wěn)定同位素背景變化觀測,對于了解構(gòu)造活動強(qiáng)烈的鮮水河-安寧河斷裂帶的地下水循環(huán)特征、斷層/含水層破裂與愈合以及大地震的孕育過程有重要意義。

    4 結(jié)論

    (1)鮮水河-安寧河斷裂帶溫泉氫氧穩(wěn)定同位素空間變化明顯,其δ18O變化范圍為-19.04‰~-12.71‰,平均值為 -16.42‰,δ2H變化范圍為-144.07‰~-88.63‰,平均值為-122.37‰。氫氧穩(wěn)定同位素主要分布在全球和區(qū)域大氣降水線上,表明了溫泉的大氣降水成因。溫泉水缺少明顯的氧同位素漂移現(xiàn)象,表明了大部分溫泉弱的水巖作用特征,并在相關(guān)離子比值、Na-K-Mg三角圖、Li和Sr元素等指標(biāo)上得到證實。

    (2)研究區(qū)域溫泉同位素值具有明顯的同位素高程效應(yīng),溫泉同位素與海拔的擬合關(guān)系式為δ18O=-0.0023h-10.011,δ2H=-0.0195h-67.813。氧同位素高程效應(yīng)為 -0.23‰/100m,氫同位素高程效應(yīng)為-1.95‰/100m。

    (3)溫泉與其附近地表水的δ18O和δ2H具有明顯的差異性,表明了溫泉的補(bǔ)給源為遠(yuǎn)距離的大氣降水。氫氧穩(wěn)定同位素特征、水巖作用特征和循環(huán)深度表明溫泉的成因為遠(yuǎn)距離大氣降水運(yùn)移補(bǔ)給地下水,地下水在地下熱儲層加熱后通過斷裂通道上升到地表形成溫泉。

    致謝衷心感謝兩位審稿專家對本文提出的寶貴意見!

    謹(jǐn)以此文祝賀沈其韓院士百歲華誕,祝沈先生健康長壽!

    猜你喜歡
    鮮水河氫氧康定
    Three New Species of Diploderma Hallowell,1861 (Reptilia: Squamata:Agamidae) from the Shaluli Mountains in Western Sichuan,China
    自制液壓儲氣式氫氧燃料電池
    氫氧燃料電池演示實驗的改進(jìn)
    2014年四川康定MS6.3和MS5.8地震的應(yīng)力觸發(fā)研究
    地震研究(2021年1期)2021-04-13 01:04:48
    康定:作家們的原鄉(xiāng)記憶和故鄉(xiāng)神話——以“康定七箭”的鄉(xiāng)土小說為例
    阿來研究(2020年1期)2020-10-28 08:09:58
    汶川MS8.0大震前后鮮水河斷裂帶滑動及其西部地殼形變的動態(tài)變化
    鮮水河斷裂帶中南段同位素年代學(xué)及其地質(zhì)意義
    連鑄坯氫氧切割應(yīng)用與碳排放研究
    走進(jìn)康定
    鮮水河斷裂的幾何形態(tài)對地震發(fā)生的影響1
    伊人久久国产一区二区| 午夜精品一区二区三区免费看| 国产一区二区三区综合在线观看 | 亚洲精品乱码久久久v下载方式| 欧美日本视频| 日韩欧美 国产精品| 草草在线视频免费看| 观看美女的网站| 国产精品一区www在线观看| 极品少妇高潮喷水抽搐| 美女脱内裤让男人舔精品视频| 少妇被粗大猛烈的视频| 美女xxoo啪啪120秒动态图| 在现免费观看毛片| 黄色欧美视频在线观看| 边亲边吃奶的免费视频| 亚洲av成人av| 亚洲国产精品成人久久小说| 天堂av国产一区二区熟女人妻| 日本黄大片高清| 亚洲av.av天堂| 亚洲欧美日韩无卡精品| 热99在线观看视频| 日韩欧美三级三区| 麻豆国产97在线/欧美| 99九九线精品视频在线观看视频| 亚洲国产最新在线播放| 简卡轻食公司| 亚洲美女搞黄在线观看| 日韩欧美一区视频在线观看 | 色吧在线观看| 日韩av在线大香蕉| 韩国av在线不卡| 伊人久久国产一区二区| 少妇的逼水好多| 亚洲av日韩在线播放| 我的女老师完整版在线观看| 国内精品宾馆在线| 成人一区二区视频在线观看| 男女边摸边吃奶| 亚洲18禁久久av| 美女cb高潮喷水在线观看| 成年av动漫网址| 欧美+日韩+精品| 日韩av在线免费看完整版不卡| 亚洲av男天堂| 亚洲av电影在线观看一区二区三区 | 一级毛片 在线播放| 国产精品久久久久久av不卡| 色尼玛亚洲综合影院| 伦理电影大哥的女人| 3wmmmm亚洲av在线观看| 亚洲婷婷狠狠爱综合网| 乱码一卡2卡4卡精品| 成人亚洲精品一区在线观看 | 一个人看的www免费观看视频| 3wmmmm亚洲av在线观看| 日韩欧美精品免费久久| 亚洲av成人精品一区久久| 在线天堂最新版资源| 亚洲av男天堂| 亚洲国产色片| 成年女人在线观看亚洲视频 | 国产精品日韩av在线免费观看| 老司机影院成人| 天堂网av新在线| 永久网站在线| 丰满乱子伦码专区| 五月天丁香电影| 亚洲人成网站在线播| 亚洲最大成人中文| 中文欧美无线码| 观看美女的网站| 最近2019中文字幕mv第一页| 狠狠精品人妻久久久久久综合| 天堂网av新在线| 熟妇人妻久久中文字幕3abv| 中国美白少妇内射xxxbb| 国产激情偷乱视频一区二区| 中国美白少妇内射xxxbb| 国产男女超爽视频在线观看| 成人毛片60女人毛片免费| 街头女战士在线观看网站| 日韩大片免费观看网站| 亚洲最大成人av| 亚洲国产最新在线播放| 国产探花在线观看一区二区| 国产精品国产三级专区第一集| av一本久久久久| 成年免费大片在线观看| 熟妇人妻久久中文字幕3abv| 搡女人真爽免费视频火全软件| 男女视频在线观看网站免费| 亚洲在线自拍视频| 亚洲av国产av综合av卡| 一夜夜www| 亚洲国产成人一精品久久久| 插逼视频在线观看| 国产一区亚洲一区在线观看| 少妇熟女欧美另类| 大香蕉久久网| 亚洲av.av天堂| 男人和女人高潮做爰伦理| 夜夜爽夜夜爽视频| 2021天堂中文幕一二区在线观| 深爱激情五月婷婷| 亚洲精品乱码久久久v下载方式| 在线天堂最新版资源| 国产亚洲91精品色在线| 国产亚洲一区二区精品| 国产精品一区二区在线观看99 | 啦啦啦啦在线视频资源| 欧美日韩综合久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆乱淫一区二区| 亚洲成色77777| 一级毛片 在线播放| 久久久久网色| 日韩在线高清观看一区二区三区| eeuss影院久久| 亚洲精品影视一区二区三区av| 高清在线视频一区二区三区| av网站免费在线观看视频 | 国产一区有黄有色的免费视频 | 午夜老司机福利剧场| 一级爰片在线观看| 免费高清在线观看视频在线观看| 欧美成人a在线观看| 日本三级黄在线观看| 国产精品久久久久久久电影| 噜噜噜噜噜久久久久久91| 久久精品国产亚洲网站| 日日摸夜夜添夜夜添av毛片| 亚洲国产精品专区欧美| 男人和女人高潮做爰伦理| 国产又色又爽无遮挡免| 热99在线观看视频| 久久久久网色| 少妇熟女aⅴ在线视频| 久久精品夜夜夜夜夜久久蜜豆| 在线观看免费高清a一片| 乱码一卡2卡4卡精品| av网站免费在线观看视频 | 欧美变态另类bdsm刘玥| 插逼视频在线观看| 欧美bdsm另类| 欧美日韩视频高清一区二区三区二| 人妻少妇偷人精品九色| 十八禁国产超污无遮挡网站| 白带黄色成豆腐渣| 哪个播放器可以免费观看大片| 舔av片在线| 亚洲欧美中文字幕日韩二区| 亚洲人成网站高清观看| 免费av不卡在线播放| 91av网一区二区| 深爱激情五月婷婷| 99热这里只有精品一区| 肉色欧美久久久久久久蜜桃 | 亚洲精品成人av观看孕妇| 麻豆成人av视频| 黄片无遮挡物在线观看| 日韩亚洲欧美综合| 免费看光身美女| 亚洲综合色惰| 最近的中文字幕免费完整| 国产精品人妻久久久久久| 国产一区二区在线观看日韩| 简卡轻食公司| 成人亚洲精品一区在线观看 | 亚洲激情五月婷婷啪啪| 国产黄片视频在线免费观看| 精品国内亚洲2022精品成人| 91av网一区二区| 男插女下体视频免费在线播放| 成年女人在线观看亚洲视频 | 卡戴珊不雅视频在线播放| 国产大屁股一区二区在线视频| 中文天堂在线官网| 国产av国产精品国产| 国产色爽女视频免费观看| 久久久亚洲精品成人影院| 欧美激情在线99| 免费少妇av软件| 久久99精品国语久久久| 嫩草影院入口| 中国美白少妇内射xxxbb| 精品久久久久久电影网| 能在线免费看毛片的网站| 国产在视频线在精品| 欧美高清成人免费视频www| 亚洲精品国产av蜜桃| 一个人免费在线观看电影| 干丝袜人妻中文字幕| 国产又色又爽无遮挡免| 天天躁夜夜躁狠狠久久av| 久久久久久国产a免费观看| 国产片特级美女逼逼视频| 中文字幕av在线有码专区| 一区二区三区高清视频在线| 丰满人妻一区二区三区视频av| 美女cb高潮喷水在线观看| 啦啦啦韩国在线观看视频| 色网站视频免费| 国产精品av视频在线免费观看| 中国国产av一级| 久久久久网色| 日韩欧美一区视频在线观看 | 免费av观看视频| 国产精品无大码| 欧美激情在线99| 国产高清三级在线| 久久久久久久久久人人人人人人| 深夜a级毛片| 18禁在线无遮挡免费观看视频| 蜜桃久久精品国产亚洲av| 色视频www国产| 午夜视频国产福利| 欧美不卡视频在线免费观看| 国产成人精品久久久久久| 日韩中字成人| 日韩强制内射视频| 免费看av在线观看网站| 国产黄片美女视频| 亚洲熟妇中文字幕五十中出| 国产精品一二三区在线看| 国产精品一区二区在线观看99 | 91在线精品国自产拍蜜月| 亚洲精品第二区| 99久久精品一区二区三区| 我的老师免费观看完整版| 日韩欧美三级三区| 2021少妇久久久久久久久久久| 国产成人精品一,二区| 亚洲av日韩在线播放| 国产高清有码在线观看视频| 午夜免费男女啪啪视频观看| 日韩三级伦理在线观看| 高清毛片免费看| 日韩av在线免费看完整版不卡| 一区二区三区免费毛片| 伊人久久精品亚洲午夜| 有码 亚洲区| 亚洲精品乱码久久久v下载方式| freevideosex欧美| 亚洲欧美中文字幕日韩二区| 亚洲电影在线观看av| 在线观看美女被高潮喷水网站| 夜夜爽夜夜爽视频| 美女黄网站色视频| 超碰97精品在线观看| 国产一级毛片七仙女欲春2| 欧美激情在线99| 久久久欧美国产精品| 亚洲一级一片aⅴ在线观看| 99久国产av精品| 中文天堂在线官网| 亚洲av免费高清在线观看| 日韩,欧美,国产一区二区三区| 国产 亚洲一区二区三区 | 天堂俺去俺来也www色官网 | 2021少妇久久久久久久久久久| 日韩,欧美,国产一区二区三区| 有码 亚洲区| 国产在视频线在精品| 久久国产乱子免费精品| 一级黄片播放器| 你懂的网址亚洲精品在线观看| 久久精品久久精品一区二区三区| 人妻一区二区av| 中文精品一卡2卡3卡4更新| 国产91av在线免费观看| 91aial.com中文字幕在线观看| 亚洲最大成人手机在线| 亚洲色图av天堂| 免费大片18禁| 人人妻人人看人人澡| 丰满人妻一区二区三区视频av| 91久久精品国产一区二区三区| 亚洲国产最新在线播放| 亚洲av成人精品一区久久| 日韩av不卡免费在线播放| 美女xxoo啪啪120秒动态图| 中文字幕制服av| 九九在线视频观看精品| 日本猛色少妇xxxxx猛交久久| 超碰97精品在线观看| 亚洲成人av在线免费| 国产成人福利小说| 国产一级毛片在线| 三级男女做爰猛烈吃奶摸视频| 一边亲一边摸免费视频| 成年免费大片在线观看| 一个人看的www免费观看视频| 天堂√8在线中文| 韩国高清视频一区二区三区| 男的添女的下面高潮视频| 能在线免费看毛片的网站| 国内精品一区二区在线观看| 18禁在线播放成人免费| 国产伦理片在线播放av一区| 99久久中文字幕三级久久日本| 亚洲高清免费不卡视频| 最近最新中文字幕大全电影3| 精品人妻偷拍中文字幕| 青春草亚洲视频在线观看| 欧美丝袜亚洲另类| 波野结衣二区三区在线| 国产一区二区亚洲精品在线观看| 午夜免费观看性视频| 美女黄网站色视频| 免费黄频网站在线观看国产| 久久久久久久久大av| 亚洲精品国产av蜜桃| 国产精品无大码| 秋霞在线观看毛片| 国内少妇人妻偷人精品xxx网站| 91久久精品国产一区二区成人| 久久久久久久久久成人| 深爱激情五月婷婷| 欧美人与善性xxx| 激情 狠狠 欧美| 国产激情偷乱视频一区二区| 欧美另类一区| 80岁老熟妇乱子伦牲交| 男女国产视频网站| 成人综合一区亚洲| 一边亲一边摸免费视频| 久久这里有精品视频免费| 国产精品一区二区三区四区久久| 日韩欧美国产在线观看| 少妇丰满av| 成人二区视频| 婷婷色综合大香蕉| 女的被弄到高潮叫床怎么办| 91aial.com中文字幕在线观看| 免费在线观看成人毛片| 久久热精品热| 51国产日韩欧美| 亚洲精品aⅴ在线观看| 国产在视频线精品| 国产精品久久久久久久电影| 99久国产av精品国产电影| 全区人妻精品视频| 国产精品熟女久久久久浪| 欧美区成人在线视频| 欧美成人一区二区免费高清观看| 欧美不卡视频在线免费观看| 亚洲婷婷狠狠爱综合网| 美女脱内裤让男人舔精品视频| 亚洲高清免费不卡视频| 欧美日韩综合久久久久久| 国产精品久久视频播放| 国产精品一区二区在线观看99 | 搞女人的毛片| 国产黄色小视频在线观看| 中国美白少妇内射xxxbb| 激情 狠狠 欧美| 国产精品人妻久久久久久| 男女国产视频网站| 国产午夜精品论理片| 久久久精品欧美日韩精品| 午夜免费男女啪啪视频观看| 美女xxoo啪啪120秒动态图| 日韩成人伦理影院| 超碰av人人做人人爽久久| 亚洲,欧美,日韩| 免费大片18禁| 国产毛片a区久久久久| 日本爱情动作片www.在线观看| 不卡视频在线观看欧美| 女人十人毛片免费观看3o分钟| 内射极品少妇av片p| 视频中文字幕在线观看| 精品酒店卫生间| 久久韩国三级中文字幕| 九九在线视频观看精品| 亚洲精品日本国产第一区| 永久免费av网站大全| 成人性生交大片免费视频hd| 国产极品天堂在线| 欧美丝袜亚洲另类| 美女主播在线视频| 亚洲av在线观看美女高潮| 久久这里只有精品中国| 直男gayav资源| 亚洲综合色惰| 夜夜爽夜夜爽视频| 综合色av麻豆| 国产伦一二天堂av在线观看| 日韩av在线大香蕉| 国产av码专区亚洲av| 六月丁香七月| 欧美xxⅹ黑人| 久99久视频精品免费| 国产永久视频网站| 亚洲精品成人av观看孕妇| av在线观看视频网站免费| 亚洲成人久久爱视频| 在线观看av片永久免费下载| 成年免费大片在线观看| 美女主播在线视频| 精品人妻一区二区三区麻豆| 亚洲av成人av| 亚洲欧美成人精品一区二区| 一个人免费在线观看电影| 久久久久国产网址| 欧美xxxx黑人xx丫x性爽| 精品不卡国产一区二区三区| 精品一区二区三区视频在线| 午夜福利在线在线| 男人狂女人下面高潮的视频| 久久精品久久久久久噜噜老黄| 欧美高清性xxxxhd video| 3wmmmm亚洲av在线观看| 国产综合懂色| 色哟哟·www| 亚洲欧美成人精品一区二区| 久久97久久精品| 国产精品.久久久| 国模一区二区三区四区视频| 午夜福利网站1000一区二区三区| 性色avwww在线观看| 亚洲电影在线观看av| 亚洲欧美中文字幕日韩二区| 亚洲乱码一区二区免费版| 国产v大片淫在线免费观看| 热99在线观看视频| 精品久久久噜噜| 色哟哟·www| 久久久久免费精品人妻一区二区| 一区二区三区四区激情视频| 久久久久久久久久人人人人人人| 在线 av 中文字幕| 亚洲激情五月婷婷啪啪| 久热久热在线精品观看| 亚洲国产av新网站| 成人综合一区亚洲| 一级二级三级毛片免费看| 久久精品国产亚洲av涩爱| 欧美成人a在线观看| 日本-黄色视频高清免费观看| 美女黄网站色视频| 亚洲av成人精品一区久久| 国产一区二区亚洲精品在线观看| 国产精品99久久久久久久久| 成人午夜高清在线视频| 高清午夜精品一区二区三区| 成人性生交大片免费视频hd| 亚洲欧美日韩无卡精品| 亚洲人成网站在线播| 97超视频在线观看视频| 成人午夜精彩视频在线观看| 日本爱情动作片www.在线观看| .国产精品久久| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕| 久久久久久久国产电影| xxx大片免费视频| 少妇的逼水好多| 国产黄片视频在线免费观看| 一级毛片黄色毛片免费观看视频| 久久精品夜色国产| 国产午夜精品一二区理论片| 草草在线视频免费看| 女的被弄到高潮叫床怎么办| 国产午夜精品久久久久久一区二区三区| 国产黄色视频一区二区在线观看| 全区人妻精品视频| 能在线免费看毛片的网站| 国产老妇女一区| 春色校园在线视频观看| 成年女人看的毛片在线观看| 毛片女人毛片| 永久网站在线| 亚洲精品日本国产第一区| 国产精品一区二区性色av| 老司机影院成人| 国产女主播在线喷水免费视频网站 | 日日啪夜夜爽| 久久99热这里只有精品18| 久久精品国产亚洲av天美| 精品久久久久久久久亚洲| 日本黄色片子视频| 免费看a级黄色片| 欧美3d第一页| 亚洲欧美日韩无卡精品| 91精品国产九色| 深夜a级毛片| 如何舔出高潮| 男人舔女人下体高潮全视频| 中文字幕免费在线视频6| 国产一区二区在线观看日韩| 午夜福利成人在线免费观看| 美女主播在线视频| 内地一区二区视频在线| 久久久久精品久久久久真实原创| 国产91av在线免费观看| 亚洲av成人精品一二三区| 熟妇人妻不卡中文字幕| 最新中文字幕久久久久| 内地一区二区视频在线| 99久久精品热视频| av黄色大香蕉| 热99在线观看视频| 亚洲精品一区蜜桃| 国产伦精品一区二区三区视频9| 丝瓜视频免费看黄片| 日韩制服骚丝袜av| 国产久久久一区二区三区| 九九久久精品国产亚洲av麻豆| 精品久久久久久成人av| 在线观看一区二区三区| 亚洲图色成人| av免费在线看不卡| 亚洲av成人精品一区久久| 午夜激情久久久久久久| 中文字幕久久专区| 麻豆国产97在线/欧美| 久久久欧美国产精品| 乱人视频在线观看| 伦精品一区二区三区| 天美传媒精品一区二区| 成人高潮视频无遮挡免费网站| 国精品久久久久久国模美| 一个人看视频在线观看www免费| 91在线精品国自产拍蜜月| 免费大片黄手机在线观看| 亚洲国产精品成人综合色| 三级经典国产精品| 一区二区三区免费毛片| 夜夜看夜夜爽夜夜摸| 中文乱码字字幕精品一区二区三区 | 麻豆国产97在线/欧美| 欧美高清成人免费视频www| 综合色丁香网| 99久久九九国产精品国产免费| 91精品伊人久久大香线蕉| 三级男女做爰猛烈吃奶摸视频| 亚洲va在线va天堂va国产| 永久网站在线| 日韩 亚洲 欧美在线| 蜜臀久久99精品久久宅男| 色视频www国产| av国产免费在线观看| 丝瓜视频免费看黄片| av免费在线看不卡| 国产日韩欧美在线精品| 一边亲一边摸免费视频| 人人妻人人看人人澡| 欧美日韩视频高清一区二区三区二| 看免费成人av毛片| 18禁在线无遮挡免费观看视频| 熟女电影av网| 国产亚洲av嫩草精品影院| 精品一区二区三区人妻视频| 国产成人aa在线观看| 亚洲欧美精品专区久久| 中文字幕免费在线视频6| 婷婷六月久久综合丁香| 免费黄频网站在线观看国产| 国产国拍精品亚洲av在线观看| 亚洲在线观看片| 日韩强制内射视频| 不卡视频在线观看欧美| 嫩草影院入口| 寂寞人妻少妇视频99o| 色哟哟·www| 中文字幕免费在线视频6| 神马国产精品三级电影在线观看| 麻豆av噜噜一区二区三区| 日韩欧美精品v在线| 丝瓜视频免费看黄片| 国产精品人妻久久久影院| 欧美高清成人免费视频www| 青青草视频在线视频观看| 一个人看视频在线观看www免费| 国产精品日韩av在线免费观看| 国产成人午夜福利电影在线观看| 嫩草影院精品99| 国产午夜精品论理片| 青青草视频在线视频观看| 国内精品一区二区在线观看| 最后的刺客免费高清国语| 精品熟女少妇av免费看| 一本一本综合久久| 亚洲精品亚洲一区二区| 成年免费大片在线观看| 91av网一区二区| 中文字幕亚洲精品专区| 国产精品一及| 可以在线观看毛片的网站| 狂野欧美激情性xxxx在线观看| 日本熟妇午夜| 一级毛片我不卡| 亚洲国产高清在线一区二区三| 最近2019中文字幕mv第一页| 亚洲成人中文字幕在线播放| videos熟女内射| 丰满乱子伦码专区| 国内精品一区二区在线观看| 大香蕉97超碰在线| 日韩精品青青久久久久久| 亚洲av不卡在线观看| 亚洲精品乱久久久久久| 午夜视频国产福利| 久久97久久精品| 国产午夜福利久久久久久| 亚洲精品第二区| 日韩精品有码人妻一区| 免费少妇av软件| 偷拍熟女少妇极品色| 国产高清不卡午夜福利| 亚洲成人久久爱视频| 3wmmmm亚洲av在线观看|