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

    漾濞6.4級地震前后大理地區(qū)溫泉水化學(xué)組分特征分析

    2022-10-29 04:43:40李見張翔楊志堅李朝明段美芳
    華南地震 2022年3期
    關(guān)鍵詞:洱源溫泉水劍川

    李見,張翔,楊志堅,李朝明,段美芳

    (云南省地震局,昆明 650224)

    0 引言

    流體在地球動力學(xué)過程中有不可或缺的作用,水文地球化學(xué)在地震監(jiān)測預(yù)測方面得到廣泛應(yīng)用,地震常伴有流體各組分地球化學(xué)異常變化,地下流體異常對地震孕震區(qū)及其外圍區(qū)域介質(zhì)的變化有較為靈敏的響應(yīng),某些強(qiáng)震前會出現(xiàn)地下流體群體性異常變化[1-3]。2007年普洱地震震前后溫泉水中氟離子發(fā)生明顯變化[4],2008年汶川大地震伴有地下水位和流體地球化學(xué)異?,F(xiàn)象[5],2008年汶川8.0級地震前后的5月1日至7月4日,距震中約400 km的重慶北溫泉群所有溫泉的水溫下降了1℃左右,震后海拔較低泉點K+、Na+、Ca+濃度大幅降低[6],2015年阿左旗5.8級地震,周圍區(qū)域的水化學(xué)組成變化較大,對地震響應(yīng)較為敏感[7]。

    2021年5月21日21時48分37秒(北京時間),在云南省大理白族自治州漾濞縣(25.67°N,99.87°E)發(fā)生M6.4級地震,震源深度8 km。雷興林等[8]對漾濞地震特征和成因進(jìn)行研究,認(rèn)為漾濞地震受深部流體作用的影響明顯,漾濞地震后的流體特征需要進(jìn)一步研究。此次漾濞MS6.4地震發(fā)生于滇西北地區(qū),該區(qū)域歷史上M≥6.0地震頻發(fā),且大多數(shù)地震發(fā)生于主要斷裂及其附近區(qū)域,2013年3月以來,滇西北地區(qū)ML≥4.0地震顯著增強(qiáng),且主要集中分布在維西—喬后斷裂、程海斷裂及其附近地區(qū)。劉自鳳等[9]研究了滇西北地區(qū)b值以及其應(yīng)力張量方差的空間掃描,顯示出程海斷裂、龍蟠—喬后斷裂和維西—喬后斷裂南段及鄰區(qū)低值變化特征相對顯著,認(rèn)為強(qiáng)震發(fā)生的可能性大?;趶?qiáng)震危險性的較大背景,從2019年開始,對大理地區(qū)主體斷裂(維西—喬后斷裂、程海斷裂、紅河斷裂等)附近6個溫泉每季度定期進(jìn)行采樣觀測,主要測定斷裂附近溫泉水體的水氡、水汞、水體主要陰離子、固體CO2、電導(dǎo)率、pH值、水溫等地球化學(xué)量;此次漾濞6.4級地震后,也及時對區(qū)域內(nèi)溫泉開展了現(xiàn)場采樣觀測。共獲得6個溫泉震前震后8組連續(xù)流動觀測數(shù)據(jù)。因此,本文通過分析漾濞6.4級地震前后大理地區(qū)溫泉水化學(xué)組分,分析該地區(qū)溫泉水文地球化學(xué)變化規(guī)律,探索溫泉水地球化學(xué)特征變化與漾濞地震的關(guān)系。研究成果不僅對認(rèn)識該區(qū)溫泉水地球化學(xué)特征及其演化有重要意義,而且對探索和改進(jìn)利用該區(qū)溫泉水化學(xué)組分進(jìn)行地震預(yù)測的方法也具有重要意義。

    1 地質(zhì)背景

    從構(gòu)造上看,漾濞M6.4地震發(fā)生在川滇菱形塊體西側(cè)的印支塊體上,印支地塊自西向東分為三個次級構(gòu)造單元:臨滄花崗巖帶、瀾滄江斷裂帶和蘭坪—思茅盆地。蘭坪—思茅盆地內(nèi)部構(gòu)造復(fù)雜,發(fā)育了早期的逆沖斷層、推覆構(gòu)造與晚期的走滑構(gòu)造斷層[10-11]。此次地震發(fā)生在蘭坪—思茅盆地內(nèi),離東側(cè)維西—喬后—巍山斷裂約20 km(圖1)。維西—喬后—巍山斷裂帶,北起雪龍山東廘白濟(jì)汛一帶,經(jīng)維西、通甸、喬后,止于巍山盆地南端,走向北北西,傾向NE/NW,長約280 km。大致以玉獅場、平坡為界,可將斷裂分為北段、中段和南段3段,是一條走滑兼正斷的活動斷裂[12-13]。該區(qū)域是青藏高原物質(zhì)向東南側(cè)擠出的通道[14],也是特提斯—喜馬拉雅構(gòu)造域東南部弧形構(gòu)造轉(zhuǎn)折處,處于川滇菱形塊體西南邊界紅河斷裂帶與近南北向麗江—大理斷裂系統(tǒng)的交匯區(qū)及應(yīng)力集中帶[15]。區(qū)域內(nèi)強(qiáng)震頻發(fā),圖1b給出了1925年有儀器記錄以來的M≥6.0級地震M-T圖,包括漾濞M6.4級地震在內(nèi),共發(fā)生了22次6級以上地震,占同時期云南地區(qū)地震的34%,最近的一次強(qiáng)震頻發(fā)時段是1996—2003年,期間發(fā)生1次麗江7.0級地震,4次M≥6.0級地震。漾濞M6.4級地震與上一次地震相隔12年。大理地區(qū)位于紅河斷裂帶北端,構(gòu)造上表現(xiàn)為一系列近SN向及NNE向的活動斷裂,如離此次漾濞地震約20 km的維西—喬后—巍山斷裂,中新世以來的斷陷盆地從南北兩個方向逐漸向NW走向的紅河斷裂交會,構(gòu)成斷裂東北側(cè)典型的滇西北伸展裂陷區(qū)[16]。

    圖1 區(qū)域構(gòu)造與歷史地震震中分布Fig.1 Regional tectonics and epicenters distribution of the historical earthquakes

    此次定期觀測的6個溫泉為:洱源縣牛街鄉(xiāng)練渡村溫泉(以下簡稱洱源溫泉),永勝縣片角鎮(zhèn)熱河村溫泉(以下簡稱永勝溫泉),彌渡縣彌城鎮(zhèn)高芹溫泉(以下簡稱高芹溫泉),彌渡縣彌城鎮(zhèn)白總旗溫泉(以下簡稱白總旗溫泉),云龍縣漕澗下澡堂溫泉(以下簡稱云龍溫泉),劍川縣老君山新和村溫泉(以下簡稱劍川溫泉),位置如圖2所示,除劍川溫泉離漾濞6.4級地震震中距為104 km,其余溫泉震中距都在100 km內(nèi)(表1)。洱源溫泉位于紅河斷裂和麗江—劍川斷裂交互處,該區(qū)內(nèi)構(gòu)造體系呈SN向展布,為地下熱水提供循環(huán)通道;洱源盆地基底有大面積的多期巖漿侵入玄武巖,為地下熱水提供熱源;基底上覆中、新生界碎屑巖和沖湖積黏土層組成良好隔熱保溫層,三者構(gòu)成良好的儲熱構(gòu)造[17]。彌渡高芹溫泉、白總旗溫泉近紅河斷裂南段,彌渡盆地及周邊山區(qū)的碳酸鹽巖破碎帶地區(qū)是該區(qū)域地下水的補(bǔ)給區(qū),大氣降水沿巖溶裂隙及構(gòu)造裂隙向深部運(yùn)移獲得熱量形成深遠(yuǎn)高溫?zé)崴?,在盆地總體的由北至南的地下水流動方向上,運(yùn)移至溫泉出露區(qū)域,受致密巖體或者壓性斷裂帶阻攔,熱水沿構(gòu)造上升出露成泉[18]。云龍溫泉位于蘭坪—永平斷裂西側(cè),該區(qū)位于怒江及瀾滄江沿岸地下熱水密集區(qū),屬于中低溫水熱型地?zé)犷愋?,地殼地溫梯度?℃/100 m[19]。區(qū)域內(nèi)出露地層的巖性以紫紅色砂巖、粉砂巖、粉砂質(zhì)泥巖和泥巖為特征,統(tǒng)稱紅層[20]。劍川溫泉位于維西—喬后斷裂附近,出露于黃色砂巖的裂縫中,處于蘭坪—思茅盆地的邊緣[21],位于我國著名的滇藏地?zé)釒?nèi),地?zé)峄顒訌?qiáng)烈,構(gòu)造復(fù)雜,不同構(gòu)造體系的斷裂互相重合、交接、穿插,為地下熱水將深部熱量帶至地殼淺表、孕育高溫水熱系統(tǒng)、產(chǎn)生強(qiáng)烈水熱活動提供了重要的構(gòu)造條件[22]。永勝溫泉位于程?!e川斷裂紅巖坡段,程?!e川斷裂可分為八段,為最南段,區(qū)域內(nèi)活動斷裂發(fā)育,程海斷裂縱貫[23]。

    表1 溫泉基本信息表Table 1 Basic information sheet of hot spring

    圖2 溫泉取樣點分布圖Fig.2 Distribution map of hot spring sampling points

    2 水化學(xué)組分特征分析

    對6個溫泉點分別于2020年5、7、9、11月,2021年3、5、8、11月共8次觀測采樣,獲得溫泉水氡、水汞、固體CO2、陰離子、pH、電導(dǎo)率的連續(xù)8次觀測數(shù)據(jù),以及溫泉陽離子數(shù)據(jù)(無連續(xù)數(shù)據(jù))。水氡、水汞分別采用RAD7測氡儀和RA-915W測汞儀野外現(xiàn)場測定,其余取樣實驗室測定;溫泉固體CO2利用沉淀法獲取,取50 mL溫泉水樣,加入過量飽和Ba(OH)2溶液,將Ba CO3沉淀洗凈過濾、烘干,計算固體CO2值;用intergration離子色譜儀測定常量元素,儀器RSD%在5%以內(nèi),電導(dǎo)率為0.1 nS/cm,HCO3-用滴定法測定。數(shù)據(jù)見表2、表3。

    表2 溫泉采樣時間、水化學(xué)數(shù)據(jù)Table 2 Sampling time and water chemistry data of hot spring

    表3 溫泉水樣主要陰陽離子組分(2021年8月)Table 3 Main anion and cation components of hot spring water samples(August 2021)

    (接表2)

    Piper圖解和Na-K-Mg三角圖解常被用來評價水—巖化學(xué)平衡狀態(tài)、區(qū)分不同類型水樣以及判斷地下水循環(huán)深度[24]。從Piper圖可見(圖3),洱源溫泉屬于HCO3·SO4-Na型水,白總旗溫泉屬于HCO3·Cl-Na型水,高芹、劍川、永勝、云龍溫泉均屬于HCO3-Na型水。6個溫泉水常量離子組分測試結(jié)果(表2、表3)顯示,6個溫泉所含陰離子均以HCO3-占絕對優(yōu)勢,陽離子均以Na+占絕對優(yōu)勢,這可能和鈉的強(qiáng)溶解性和強(qiáng)遷移性有關(guān)。洱源,高芹,白總旗,劍川溫泉Na+占陽離子組分60%左右,且富含K+和Ca2+,反映出這4個溫泉可能受到巖漿巖和碳酸鹽巖的共同影響,大理地區(qū)廣泛分布石炭—二疊系海相沉積灰?guī)r,較高的HCO3-、Ca2+含量可能與灰?guī)r的溶解有關(guān),該區(qū)同時分布喜山期堿性巖漿巖,這或許導(dǎo)致了較高的Na+、K+;不同的是云龍和永勝溫泉主體為Na+,但K+,Ca2+,Mg2+等陽離子基本很少,這可能和離子交換作用有關(guān),K+,Ca2+,Mg2+吸附能力強(qiáng)于Na+,圍巖中的鈉離子置換鈣鎂離子。從樣品的Na-K-Mg平衡圖解(圖4)可見,永勝溫泉水屬于部分成熟水,表明其發(fā)生部分水巖反應(yīng),溫泉水與圍巖的反應(yīng)達(dá)到部分平衡狀態(tài);其余溫泉為未成熟水,且靠近Mg1/2端,表明其水巖反應(yīng)較強(qiáng),屬于淺層水區(qū)域。6個溫泉的pH為6.64~7.54之間,中性;電導(dǎo)率為0.22~1.34之間,6個溫泉pH和電導(dǎo)率多次測值變化不大,數(shù)據(jù)較穩(wěn)定。

    圖3 溫泉水的Piper圖解Fig.3 Piper diagram of hot spring water

    圖4 溫泉水的Na-K-Mg圖解Fig.4 Na-K-Mg ternary diagram of hot spring water

    3 各化學(xué)組分時空分布特征

    3.1 各化學(xué)化學(xué)組分空間分布特征

    根據(jù)不同的溫泉各組分空間分布圖(圖5)可以看出:K+、Ca2+、Mg2+、Na+、Cl-、SO42-、HCO3-、溫泉水溫(T)存在差異。區(qū)域內(nèi)溫泉都為中高溫溫泉(40℃~70℃),高溫溫泉(60℃~80℃)有洱源溫泉,白總旗溫泉和云龍溫泉,中溫溫泉有劍川溫泉,高芹溫泉,永勝溫泉,這和前人研究基本吻合[25]。洱源溫泉的陰陽離子含量、水溫(T)和固體CO2含量均最高,SO42-,HCO3-,Ca2+,Na+濃度分別達(dá)147.1 mg/L、627.5 mg/L、76.5 mg/L、219.3 mg/L、水溫74.0℃,固體CO2值最高達(dá)到262 mg/L。大理地區(qū)廣泛分布石炭—二疊系海相沉積灰?guī)r,同時洱源地區(qū)溫泉的CO2更多來源以地幔來源和變質(zhì)成因混合的深部來源CO2[24],極高的HCO3-和固體CO2值可能受碳酸鹽礦物溶解和深部來源共同影響。洱源盆地基底有大面積的多期巖漿侵入玄武巖,洱源溫泉可能同時受到巖漿巖和碳酸鹽巖的共同影響,且溫泉水溶蝕了沉積巖中的石膏或者其他硫酸鹽礦物;大面積的多期巖漿侵入玄武巖,以及上覆的中、新生界碎屑巖和沖湖積黏土層良好的儲熱構(gòu)造,可能是導(dǎo)致洱源溫泉水溫較高的主要原因。云龍和永勝溫泉各水化學(xué)組分均含量較低(圖5),較低的K+、Ca2+、Mg2+含量,可能和離子交換作用有關(guān),即圍巖中的鈉離子置換鈣鎂離子;彌渡白總旗溫泉和高芹溫泉HCO3-含量(413.3 mg/L)和固體CO2含量較高,彌渡盆地及周邊山區(qū)的碳酸鹽巖破碎帶地區(qū)是該區(qū)域地下水的補(bǔ)給區(qū)[18],這或許導(dǎo)致相對較高的HCO3-和固體CO2含量。各溫泉水氡,水汞差異性較大,位于蘭坪—永平斷裂西側(cè)的云龍溫泉的水氡含量最高,最高達(dá)到283 Bq/L,漾濞地震震中東南側(cè)紅河斷裂,程?!e川斷裂附近的洱源、白總旗、高芹、永勝溫泉水氡含量較低,基本在20 Bq/L以下。彌渡高芹溫泉水汞含量最高,最高達(dá)462.2 ng/L,但彌渡地區(qū)的白總旗溫泉和高芹溫泉同處紅河斷裂附近,其水汞含量相對較低,僅為15 ng/L左右,顯示出不同斷裂帶及同一斷裂帶上不同區(qū)位的溫泉水汞,水氡背景值各不相同,具有一定的區(qū)域性。

    圖5 溫泉各組分的空間分布圖Fig.5 Spatial distribution of components in hot springs

    3.2 水化學(xué)組分時間變化

    以漾濞6.4級地震前2021年3月采集的水樣為第一批水樣,計算6個溫泉后續(xù)各批水樣水化學(xué)參數(shù)相對于各自的第1批水樣各參數(shù)的變化率并作圖(圖6)。地震前后各化學(xué)參數(shù)的變化率分析發(fā)現(xiàn),漾濞地震后各斷裂帶溫泉水的化學(xué)參數(shù)出現(xiàn)了不同程度的變化。陰離子主要表現(xiàn)在Cl-、SO42-出現(xiàn)不同程度的變化,除云龍溫泉的Cl-、SO42-變化幅度為正,明顯上升,其余溫泉均表現(xiàn)出震后下降趨勢,且變化幅度較大,如劍川溫泉Cl-下降幅度最大,相對震前采集的水樣,下降幅度分別達(dá)到99%、63%、98%;永勝溫泉SO42-相對第一批水樣,其下降幅度達(dá)到98%。水氡,水汞相對震前水樣,部分溫泉為上升趨勢,除高芹水氡外,變化幅度不大,基本在20%以內(nèi);其余大部分溫泉變化幅度為負(fù),表現(xiàn)為顯著下降趨勢,如高芹溫泉水汞在震前達(dá)到高值,后持續(xù)下降,下降幅度約40%。

    圖6 溫泉各組分濃度的相對變化率Fig.6 The relative change rate of the concentration of each component in the hot spring

    根據(jù)漾濞地震前后大理地區(qū)各溫泉水化學(xué)參數(shù)隨時間變化圖(圖7)分析顯示。位于紅河斷裂附近的彌渡高芹,白總旗溫泉、維西—喬后斷裂附近的劍川溫泉、程?!e川斷裂附近的永勝溫泉的Cl-濃度存在同步變化趨勢,均在漾濞地震前出現(xiàn)明顯上升,后迅速下降發(fā)生地震,地震后降至上升前水平,劍川,白總旗溫泉變化最為顯著,如劍川溫泉Cl-從震前的4.6 mg/L降至0.05 mg/L(圖7a,7b)。永勝溫泉SO42-在漾濞地震前升至較高值后突然下降,由2021年3月23.8 mg/L下降至5月19日0.4 mg//L,下降幅度達(dá)到98%(圖7c),洱源溫泉具有季節(jié)性變化趨勢,漾濞地震前打破季節(jié)性趨勢,其余溫泉SO42-變化不明顯,彌渡高芹溫泉受季節(jié)性變化影響,長期趨勢穩(wěn)定。高芹溫泉水汞反應(yīng)較為靈敏,震前明顯升高,后持續(xù)下降發(fā)生地震,由震前462 ng/L下降295 ng/L;云龍溫泉水汞存在相同趨勢,但對比上一年度變化,可能和季節(jié)性有關(guān)。永勝水氡從2020年7月開始持續(xù)一年多下降,漾濞地震發(fā)生時處于較低值。6個溫泉的固體CO2表現(xiàn)出一定季節(jié)性變規(guī)律,變化不明顯。

    圖7 溫泉各化學(xué)組分濃度時間變化圖Fig.7 Time variation of chemical components in hot springs

    4 溫泉水的地球化學(xué)變化與漾濞6.4級地震的關(guān)系

    根據(jù)前文分析顯示,漾濞6.4級地震前后區(qū)域內(nèi)溫泉水的化學(xué)參數(shù)主要表現(xiàn)在Cl-和SO42-濃度以及溫泉水氡、水汞含量出現(xiàn)不同程度的變化,Cl-和SO42-變化最為顯著,如劍川溫泉Cl-,永勝溫泉SO42-,相對震前采集的水樣,下降幅度都達(dá)到98%以上,溫泉的水溫和固體CO2變化不明顯??傮w顯示出漾濞6.4級地震前大多數(shù)溫泉Cl-、SO42-、水氡、水汞濃度震前出現(xiàn)明顯上升,震后各組分含量顯著下降。許多學(xué)者對地下水前兆異常的機(jī)制進(jìn)行了研究,提出了許多解釋地下水前兆異常機(jī)制的模型[[27-30],其中含水層破裂/流體混合模型(AB/FM模型)認(rèn)為隔離含水層之間的水文屏障的前兆破裂,從而使不同含水不層發(fā)生混合,而這種前兆破裂是由地震前與應(yīng)力積累相關(guān)的地殼膨脹引起的。大理地區(qū)6個溫泉的水化學(xué)變化現(xiàn)象可以用AB/FM模型模型機(jī)制解釋。

    前人研究滇西北地區(qū)b值和應(yīng)力張量方差的空間掃描結(jié)果顯示程海斷裂、龍蟠—喬后斷裂和維西—喬后斷裂南段及鄰區(qū)低值變化特征相對顯著,這可能是區(qū)域應(yīng)力水平增強(qiáng)的表現(xiàn)[9],此次研究Cl-和SO42-震前升高,震后劇烈下降的劍川,永勝,白總旗,高芹等溫泉剛好在這些斷裂附近。劍川溫泉位于維西—喬后斷裂附近,永勝溫泉位于程海—賓川斷裂紅巖坡段,白總旗和高芹溫泉近紅河斷裂南段。既區(qū)域應(yīng)力積累增加,孕震區(qū)的巖石的內(nèi)部結(jié)構(gòu)經(jīng)歷巖石的彈性壓縮、原生裂隙的閉合、新微裂隙的產(chǎn)生、裂隙的彌散發(fā)展、裂隙的集結(jié)到巖石的宏觀破裂的演化過程[26],演化過程中深部流體擴(kuò)散上升、不同層位水體的混合以及孔隙流體的受擠壓流體巖石間的相互作用導(dǎo)致地下流體水化組分的異常變化[27-30]。洱源溫泉屬于HCO3·SO4-Na型水,白總旗溫泉屬于HCO3·Cl-Na型水,高芹、劍川、永勝、云龍溫泉均屬于HCO3-Na型水,雖然這些溫泉有一定含量的Cl-和SO42-,卻不是溫泉的主量化學(xué)元素,而SO42-及Cl-除少量的地表來源外,主要為中、深部來源[7],這些溫泉水較高的Cl-、SO42-濃度可能主要來源于深部高Cl-、SO42-濃度水的混入,地震孕育過程中,區(qū)域壓應(yīng)力增強(qiáng)導(dǎo)致高Cl-、SO42-濃度的深部水供給量的增加所致,地震后,隨著區(qū)域應(yīng)力的調(diào)整,深部供給的流體量減少,各溫泉的離子濃度下降。

    5 結(jié)論

    (1)大理地區(qū)主要活動斷裂附近6個溫泉主要為中高溫溫泉,洱源溫泉溫度最高,約74.5℃。溫泉水的常量離子分析表明,陽離子主要以Na+為主,Ca2+、Mg2+含量不高且差異較大,位于程?!e川斷裂的永勝溫泉和蘭坪—永平斷裂附近的云龍溫泉Ca2+、Mg2+含量極少,陰離子以HCO3-為主,較高的HCO3-和固體CO2含量可能受大理地區(qū)廣泛分布的石炭—二疊系海相沉積灰?guī)r影響。洱源溫泉屬于HCO3·SO4-Na型水,白總旗溫泉屬于HCO3·Cl-Na型水,高芹、劍川、永勝、漕澗溫泉均屬于HCO3-Na型水。永勝溫泉水屬于部分成熟水,其余溫泉為未成熟水。區(qū)域內(nèi)溫泉有較高水氡水汞背景值,不同斷裂帶及同一斷裂帶上不同區(qū)位的溫泉水汞,水氡背景值各不相同,具有一定的區(qū)域性。

    (2)漾濞6.4級地震前后大理地區(qū)主要斷裂附近溫泉水的Cl-、SO42-、水汞、水氡含量有明顯變化,白總旗溫泉Cl-,永勝、洱源溫泉SO42-以及高芹溫泉水汞含量在漾濞6.4級地震前明顯上升,震后含量顯著下降;可能是震前區(qū)域擠壓應(yīng)力增強(qiáng)過程高Cl-和SO42-,水體的混入,震后區(qū)域應(yīng)力的調(diào)整,中深部流體量補(bǔ)給量減少,各化學(xué)組分多數(shù)含量下降;溫泉水溫和固體CO2含量變化不明顯。

    猜你喜歡
    洱源溫泉水劍川
    千年繁華后的平淡劍川彌井村
    云南畫報(2021年5期)2021-07-22 08:45:18
    含西藏溫泉水的修護(hù)霜對敏感性皮膚輔助治療作用的臨床觀察
    劍川娃
    民族音樂(2019年3期)2019-08-14 01:05:18
    溫泉水輸送管道設(shè)計工程實例分析
    遇見最美洱源
    大理文化(2019年4期)2019-05-10 00:20:04
    喝溫泉水不如喝直飲水
    劍川花椒的主要栽培技術(shù)及病蟲害防治
    畫中有話
    洱源“雀之靈”的故鄉(xiāng)綻放上海
    云南畫報(2015年10期)2015-12-02 05:21:23
    洱源荷斯坦奶牛發(fā)情期各種陰道出血現(xiàn)象的診斷和治療
    99国产精品一区二区蜜桃av| 久久午夜福利片| 色播亚洲综合网| 国产一区二区在线av高清观看| 亚洲精品成人久久久久久| 亚洲精品在线观看二区| 久久香蕉精品热| 一区二区三区激情视频| 欧美+亚洲+日韩+国产| 亚洲专区中文字幕在线| 久久久久久九九精品二区国产| 一进一出抽搐gif免费好疼| 国产高潮美女av| 欧美日韩黄片免| 久久热精品热| 国产午夜福利久久久久久| 欧美最黄视频在线播放免费| 天堂av国产一区二区熟女人妻| 男女那种视频在线观看| 国产综合懂色| 亚洲av电影不卡..在线观看| 男人狂女人下面高潮的视频| 亚洲精品456在线播放app | 日韩欧美精品免费久久 | 欧美区成人在线视频| 精品一区二区三区人妻视频| 日韩免费av在线播放| 成人永久免费在线观看视频| 亚洲第一电影网av| 精品日产1卡2卡| 亚洲欧美精品综合久久99| 狂野欧美白嫩少妇大欣赏| 中文字幕精品亚洲无线码一区| 久久午夜亚洲精品久久| 欧美成人一区二区免费高清观看| 欧美xxxx黑人xx丫x性爽| 欧美乱妇无乱码| 国产亚洲av嫩草精品影院| 别揉我奶头~嗯~啊~动态视频| 在线免费观看不下载黄p国产 | 男人狂女人下面高潮的视频| 搞女人的毛片| 国产真实伦视频高清在线观看 | 成年免费大片在线观看| 亚洲成av人片在线播放无| 国产精品久久久久久亚洲av鲁大| 琪琪午夜伦伦电影理论片6080| 麻豆国产av国片精品| av在线老鸭窝| 成年女人看的毛片在线观看| 国产欧美日韩一区二区三| 天天躁日日操中文字幕| 内地一区二区视频在线| 精品福利观看| 高清日韩中文字幕在线| 午夜老司机福利剧场| 久久伊人香网站| 99久久九九国产精品国产免费| 18禁黄网站禁片免费观看直播| 一级毛片久久久久久久久女| 老司机午夜十八禁免费视频| 亚洲在线观看片| 成人欧美大片| 亚洲精品乱码久久久v下载方式| 老司机午夜福利在线观看视频| 国产日本99.免费观看| 欧美丝袜亚洲另类 | 老司机午夜十八禁免费视频| 国产精品99久久久久久久久| 国产伦在线观看视频一区| 少妇高潮的动态图| 日本成人三级电影网站| 国产免费一级a男人的天堂| 亚洲精品亚洲一区二区| 美女免费视频网站| 丁香欧美五月| 亚洲熟妇中文字幕五十中出| 免费大片18禁| 欧美成狂野欧美在线观看| 一级黄片播放器| 久久精品国产99精品国产亚洲性色| 久久久久国内视频| bbb黄色大片| 色尼玛亚洲综合影院| 午夜精品一区二区三区免费看| 亚洲午夜理论影院| 人妻夜夜爽99麻豆av| 久久国产乱子伦精品免费另类| 日本精品一区二区三区蜜桃| 性色av乱码一区二区三区2| 欧美乱妇无乱码| 久久草成人影院| 亚洲成人久久爱视频| 国内精品久久久久久久电影| 高潮久久久久久久久久久不卡| 国产伦在线观看视频一区| 日本黄色视频三级网站网址| 亚洲成人精品中文字幕电影| 欧美黑人欧美精品刺激| 亚洲黑人精品在线| 一进一出抽搐动态| 在线播放国产精品三级| 99久久99久久久精品蜜桃| 小说图片视频综合网站| 国产av在哪里看| 亚洲av成人av| 亚洲精品影视一区二区三区av| 午夜激情福利司机影院| 高清在线国产一区| 国产午夜精品论理片| 在线天堂最新版资源| 成人三级黄色视频| 九色成人免费人妻av| 18美女黄网站色大片免费观看| 亚洲国产色片| 夜夜夜夜夜久久久久| 亚洲欧美日韩高清在线视频| 亚洲av二区三区四区| 99久久久亚洲精品蜜臀av| av欧美777| 亚洲性夜色夜夜综合| 最后的刺客免费高清国语| 一进一出抽搐动态| 麻豆av噜噜一区二区三区| 精品国内亚洲2022精品成人| netflix在线观看网站| 91在线观看av| 亚洲成人免费电影在线观看| 日韩人妻高清精品专区| 亚洲国产精品sss在线观看| 91麻豆精品激情在线观看国产| 舔av片在线| 国产成人影院久久av| 久久99热6这里只有精品| 真实男女啪啪啪动态图| 在线播放无遮挡| 永久网站在线| 九九热线精品视视频播放| 97超级碰碰碰精品色视频在线观看| 免费看光身美女| 精品国产三级普通话版| 99在线视频只有这里精品首页| 真人做人爱边吃奶动态| 精品人妻偷拍中文字幕| 成人特级av手机在线观看| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| av国产免费在线观看| 最近最新免费中文字幕在线| 亚洲在线自拍视频| 国产av在哪里看| 日韩av在线大香蕉| 国内精品久久久久精免费| 久久精品人妻少妇| 亚洲av第一区精品v没综合| 亚洲美女搞黄在线观看 | 国产熟女xx| 丁香欧美五月| 嫁个100分男人电影在线观看| 级片在线观看| 久久精品人妻少妇| 99久久无色码亚洲精品果冻| 直男gayav资源| 国产精品98久久久久久宅男小说| 亚洲精品一区av在线观看| 男女床上黄色一级片免费看| 中亚洲国语对白在线视频| 国产美女午夜福利| 中文字幕人妻熟人妻熟丝袜美| 国产激情偷乱视频一区二区| 黄色一级大片看看| 97碰自拍视频| 天堂动漫精品| а√天堂www在线а√下载| 麻豆国产av国片精品| 国产又黄又爽又无遮挡在线| 国产亚洲精品综合一区在线观看| 一区二区三区高清视频在线| 亚洲av.av天堂| 国产淫片久久久久久久久 | 天堂动漫精品| 久久精品国产99精品国产亚洲性色| 国产视频一区二区在线看| 国产成人福利小说| 精品午夜福利视频在线观看一区| 精品国产亚洲在线| 在线观看美女被高潮喷水网站 | 国产中年淑女户外野战色| 一个人看视频在线观看www免费| 91麻豆精品激情在线观看国产| 少妇人妻精品综合一区二区 | 国产黄a三级三级三级人| 久久99热6这里只有精品| 久久国产精品影院| 极品教师在线免费播放| 黄色丝袜av网址大全| 老熟妇仑乱视频hdxx| 天美传媒精品一区二区| 免费在线观看日本一区| ponron亚洲| 亚洲内射少妇av| 久久久久久大精品| 日韩av在线大香蕉| 色综合欧美亚洲国产小说| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久热精品热| 欧美日韩中文字幕国产精品一区二区三区| 色播亚洲综合网| 午夜福利成人在线免费观看| 亚洲,欧美,日韩| 国产三级在线视频| 国产精品精品国产色婷婷| 日本撒尿小便嘘嘘汇集6| 少妇的逼水好多| 亚洲成人久久性| 丰满人妻一区二区三区视频av| 亚洲黑人精品在线| 欧美高清性xxxxhd video| 久久国产精品影院| 国产免费av片在线观看野外av| 日韩av在线大香蕉| 亚洲男人的天堂狠狠| 国产精品乱码一区二三区的特点| 男人舔奶头视频| 国产精品三级大全| 1024手机看黄色片| 国产色婷婷99| 亚洲中文字幕日韩| 97超视频在线观看视频| 在线观看午夜福利视频| 国产精品久久电影中文字幕| 免费无遮挡裸体视频| 日韩成人在线观看一区二区三区| 国产精品爽爽va在线观看网站| 一进一出抽搐动态| 淫妇啪啪啪对白视频| 国产一区二区三区视频了| 亚洲综合色惰| 国产成人影院久久av| 免费在线观看日本一区| 国产极品精品免费视频能看的| 香蕉av资源在线| 亚洲最大成人手机在线| 欧美国产日韩亚洲一区| 亚洲经典国产精华液单 | 免费搜索国产男女视频| 人妻久久中文字幕网| 免费无遮挡裸体视频| 久久亚洲精品不卡| 一边摸一边抽搐一进一小说| 中文字幕免费在线视频6| 国产毛片a区久久久久| 国产一区二区在线观看日韩| 啪啪无遮挡十八禁网站| 国产成人啪精品午夜网站| 免费无遮挡裸体视频| 婷婷六月久久综合丁香| a在线观看视频网站| 久久人妻av系列| 国产极品精品免费视频能看的| 亚洲av成人av| 非洲黑人性xxxx精品又粗又长| 好看av亚洲va欧美ⅴa在| 亚洲欧美日韩无卡精品| 女人十人毛片免费观看3o分钟| 亚洲色图av天堂| 91午夜精品亚洲一区二区三区 | 宅男免费午夜| 国产真实乱freesex| 成人精品一区二区免费| 国产野战对白在线观看| 91在线精品国自产拍蜜月| 亚洲男人的天堂狠狠| 国产精品久久久久久久电影| 国产精品一区二区三区四区久久| 亚洲人成电影免费在线| 国产在线男女| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 一卡2卡三卡四卡精品乱码亚洲| h日本视频在线播放| 在线观看一区二区三区| 桃红色精品国产亚洲av| 成人一区二区视频在线观看| 波多野结衣高清无吗| 国产精品女同一区二区软件 | 又爽又黄a免费视频| 国内精品美女久久久久久| 国产一区二区激情短视频| 两人在一起打扑克的视频| 亚洲国产精品合色在线| 久久精品国产清高在天天线| 午夜精品久久久久久毛片777| 18+在线观看网站| 日韩欧美精品v在线| 欧美成人a在线观看| 亚洲人成网站高清观看| 久久久久久九九精品二区国产| 欧美性感艳星| 欧洲精品卡2卡3卡4卡5卡区| 网址你懂的国产日韩在线| 久久精品国产清高在天天线| 午夜精品在线福利| 可以在线观看毛片的网站| 亚洲美女黄片视频| 中亚洲国语对白在线视频| 欧美最黄视频在线播放免费| 2021天堂中文幕一二区在线观| 色播亚洲综合网| 亚洲在线观看片| 精品人妻1区二区| 亚洲欧美日韩卡通动漫| 一卡2卡三卡四卡精品乱码亚洲| 熟女电影av网| 国产黄a三级三级三级人| 69av精品久久久久久| 成人亚洲精品av一区二区| 国产亚洲精品av在线| 村上凉子中文字幕在线| 亚洲国产精品成人综合色| 99久久成人亚洲精品观看| 亚洲成人久久性| 久久精品综合一区二区三区| 亚洲一区高清亚洲精品| 人妻制服诱惑在线中文字幕| 欧美黑人欧美精品刺激| 又爽又黄无遮挡网站| 精品一区二区三区视频在线观看免费| 变态另类丝袜制服| 99久国产av精品| 精品国产亚洲在线| 90打野战视频偷拍视频| 亚洲五月婷婷丁香| 欧美成人a在线观看| 夜夜躁狠狠躁天天躁| 老熟妇仑乱视频hdxx| 18禁在线播放成人免费| 亚洲不卡免费看| 老熟妇仑乱视频hdxx| 精品免费久久久久久久清纯| 国产高清三级在线| 久久人人精品亚洲av| 精品久久久久久久末码| 成人精品一区二区免费| 精品国产三级普通话版| 国产淫片久久久久久久久 | 亚洲精品在线美女| 免费高清视频大片| 美女 人体艺术 gogo| 国内精品一区二区在线观看| 有码 亚洲区| 久久人妻av系列| 2021天堂中文幕一二区在线观| 国模一区二区三区四区视频| 免费一级毛片在线播放高清视频| 91av网一区二区| 日本一二三区视频观看| 一本综合久久免费| 色综合亚洲欧美另类图片| 久99久视频精品免费| 欧美精品啪啪一区二区三区| 国产精品,欧美在线| 很黄的视频免费| 99视频精品全部免费 在线| 国产欧美日韩精品亚洲av| 色视频www国产| 精品不卡国产一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 欧美国产日韩亚洲一区| 国产乱人视频| 一边摸一边抽搐一进一小说| 国产精品人妻久久久久久| 色综合亚洲欧美另类图片| 很黄的视频免费| 日本成人三级电影网站| 18禁在线播放成人免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 十八禁网站免费在线| 国产成人欧美在线观看| 老熟妇乱子伦视频在线观看| 九色成人免费人妻av| 色av中文字幕| 两人在一起打扑克的视频| 天天一区二区日本电影三级| 五月玫瑰六月丁香| 国产中年淑女户外野战色| av天堂中文字幕网| 国产三级黄色录像| 日本三级黄在线观看| 国产免费一级a男人的天堂| 91在线精品国自产拍蜜月| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 精品熟女少妇八av免费久了| 网址你懂的国产日韩在线| 亚洲av电影在线进入| www.熟女人妻精品国产| 人人妻人人澡欧美一区二区| 日韩亚洲欧美综合| 看十八女毛片水多多多| 我要看日韩黄色一级片| 精品久久久久久久久久久久久| 少妇人妻一区二区三区视频| 最近在线观看免费完整版| 伦理电影大哥的女人| 久久久久久国产a免费观看| 国产高清三级在线| 亚洲国产欧美人成| 国产精品亚洲美女久久久| 91在线精品国自产拍蜜月| 国产白丝娇喘喷水9色精品| 毛片一级片免费看久久久久 | 欧美日韩瑟瑟在线播放| 一边摸一边抽搐一进一小说| 国产免费一级a男人的天堂| 国产精品美女特级片免费视频播放器| 自拍偷自拍亚洲精品老妇| 久久草成人影院| 一区福利在线观看| 亚洲人成网站高清观看| 亚洲无线在线观看| 亚洲最大成人中文| 精品久久久久久久末码| 1000部很黄的大片| 两人在一起打扑克的视频| 搡老岳熟女国产| 蜜桃久久精品国产亚洲av| 91麻豆精品激情在线观看国产| 久久久久精品国产欧美久久久| 五月伊人婷婷丁香| 久久久久九九精品影院| 麻豆国产av国片精品| 成人永久免费在线观看视频| 日本精品一区二区三区蜜桃| 久久精品国产亚洲av天美| 亚洲精品一区av在线观看| 简卡轻食公司| 亚洲av五月六月丁香网| 国产av不卡久久| 国产在视频线在精品| 亚洲成人中文字幕在线播放| 国产成人av教育| 精品国产亚洲在线| 亚洲精华国产精华精| 特级一级黄色大片| 噜噜噜噜噜久久久久久91| 成人欧美大片| 久久中文看片网| 欧美高清性xxxxhd video| 日韩成人在线观看一区二区三区| 国产人妻一区二区三区在| 亚洲国产高清在线一区二区三| 国产精品女同一区二区软件 | 国产在视频线在精品| 少妇被粗大猛烈的视频| 欧美区成人在线视频| 国产伦人伦偷精品视频| 国产精品野战在线观看| 男女床上黄色一级片免费看| 99久久成人亚洲精品观看| 91午夜精品亚洲一区二区三区 | 在线免费观看不下载黄p国产 | 神马国产精品三级电影在线观看| 久久午夜福利片| 黄色一级大片看看| 我要搜黄色片| 国产亚洲精品综合一区在线观看| 国产一区二区激情短视频| 亚洲在线观看片| 在线免费观看的www视频| 白带黄色成豆腐渣| 国产欧美日韩精品亚洲av| av在线老鸭窝| 色综合亚洲欧美另类图片| 极品教师在线视频| 午夜影院日韩av| 在线免费观看不下载黄p国产 | 久久人人精品亚洲av| 非洲黑人性xxxx精品又粗又长| 悠悠久久av| 精品欧美国产一区二区三| 成人特级黄色片久久久久久久| 亚洲精华国产精华精| 伦理电影大哥的女人| 岛国在线免费视频观看| 在线十欧美十亚洲十日本专区| 露出奶头的视频| 成人一区二区视频在线观看| 亚洲18禁久久av| av在线老鸭窝| 欧美午夜高清在线| 国产日本99.免费观看| 成年女人看的毛片在线观看| 亚洲精品影视一区二区三区av| 日韩精品中文字幕看吧| 狠狠狠狠99中文字幕| 91狼人影院| or卡值多少钱| 国产在线男女| 久久久久久久午夜电影| 成人欧美大片| 亚洲国产欧洲综合997久久,| 精品久久久久久久久久久久久| 午夜a级毛片| 国产成人aa在线观看| 在线a可以看的网站| 亚洲精华国产精华精| 欧美最新免费一区二区三区 | 色av中文字幕| 久久精品久久久久久噜噜老黄 | 亚洲人成伊人成综合网2020| 婷婷精品国产亚洲av在线| 国产精品亚洲av一区麻豆| 日本五十路高清| 一级黄片播放器| 午夜免费成人在线视频| 99久久精品国产亚洲精品| 最好的美女福利视频网| 亚洲一区高清亚洲精品| 无遮挡黄片免费观看| 成年女人毛片免费观看观看9| 亚洲乱码一区二区免费版| 久久久久性生活片| 亚洲成人久久性| 日韩 亚洲 欧美在线| 午夜福利在线在线| 成年免费大片在线观看| 亚洲av成人不卡在线观看播放网| 欧美xxxx性猛交bbbb| 国产乱人视频| 免费在线观看亚洲国产| 久久国产乱子伦精品免费另类| 成年女人毛片免费观看观看9| 一级a爱片免费观看的视频| 一卡2卡三卡四卡精品乱码亚洲| 久久精品夜夜夜夜夜久久蜜豆| 蜜桃久久精品国产亚洲av| 精品一区二区免费观看| 观看免费一级毛片| 男女床上黄色一级片免费看| 免费人成在线观看视频色| 丰满乱子伦码专区| 91狼人影院| 一区二区三区四区激情视频 | 久久精品夜夜夜夜夜久久蜜豆| 欧美国产日韩亚洲一区| 黄片小视频在线播放| 国产高清三级在线| 国产一区二区在线av高清观看| 亚洲午夜理论影院| 最近中文字幕高清免费大全6 | 99久久精品一区二区三区| 香蕉av资源在线| 欧美极品一区二区三区四区| 亚洲精品成人久久久久久| 欧美高清性xxxxhd video| 无遮挡黄片免费观看| 能在线免费观看的黄片| 午夜免费成人在线视频| 精品日产1卡2卡| 亚洲综合色惰| 国产探花在线观看一区二区| 国产精品日韩av在线免费观看| 国产成+人综合+亚洲专区| 毛片一级片免费看久久久久 | 亚洲国产精品合色在线| 日韩欧美精品免费久久 | 国产麻豆成人av免费视频| 精品久久久久久,| 99热这里只有精品一区| 欧美黄色片欧美黄色片| 麻豆一二三区av精品| 偷拍熟女少妇极品色| h日本视频在线播放| 男人和女人高潮做爰伦理| 大型黄色视频在线免费观看| 女人被狂操c到高潮| 国内揄拍国产精品人妻在线| 又黄又爽又刺激的免费视频.| 国产亚洲欧美在线一区二区| 国产精品亚洲美女久久久| 久久6这里有精品| 波野结衣二区三区在线| 青草久久国产| 在线观看av片永久免费下载| 亚洲一区高清亚洲精品| avwww免费| 国产色爽女视频免费观看| 成年女人永久免费观看视频| netflix在线观看网站| 精品免费久久久久久久清纯| 欧美最新免费一区二区三区 | 人人妻人人看人人澡| 一夜夜www| 日本成人三级电影网站| 久久久国产成人免费| 欧美日本亚洲视频在线播放| 91麻豆av在线| 国产精品免费一区二区三区在线| 真人做人爱边吃奶动态| 少妇裸体淫交视频免费看高清| 老熟妇乱子伦视频在线观看| 毛片女人毛片| 波多野结衣高清无吗| 成人亚洲精品av一区二区| 亚洲,欧美精品.| 人人妻人人看人人澡| 成人午夜高清在线视频| 小说图片视频综合网站| 日韩有码中文字幕| av在线观看视频网站免费| 99热这里只有精品一区| 久久久久久久久久黄片| 天堂av国产一区二区熟女人妻| 性色av乱码一区二区三区2|