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

    云南南部地區(qū)深部電性結(jié)構(gòu)特征研究

    2014-09-25 02:15:20李冉湯吉董澤義肖騎彬詹艷
    地球物理學報 2014年4期
    關鍵詞:震區(qū)電性斷裂帶

    李冉,湯吉,2*,董澤義,肖騎彬,2,詹艷,2

    1中國地震局地質(zhì)研究所,北京 100029

    2地震動力學國家重點實驗室,北京 100029

    1 引言

    云南地區(qū)(21°N—29°N,97°E—106°E)位于青藏高原東南緣,是印度板塊與歐亞大陸北東向碰撞擠壓的前沿地帶和深部物質(zhì)逃逸地帶,由于青藏高原隆升產(chǎn)生的構(gòu)造運動及地殼深部物質(zhì)逃逸的影響,使這一地區(qū)的構(gòu)造活動非常復雜,地震活動也十分強烈,是中國大陸顯著的強震活動地區(qū)之一(周海濤等,2011).

    有研究表明,每年由印度板塊與青藏塊體碰撞作用向云南地區(qū)累計輸入約相當于一次6.8級地震的能量(皇甫崗等,2010),因此云南地區(qū)的孕震能力很強.據(jù)統(tǒng)計,僅從20世紀以來,云南有記錄的5級以上地震327次(部分數(shù)據(jù)來自中國地震臺網(wǎng)地震目錄),其中6.0~6.9級強震69次,7級以上大震13次(皇甫崗等,2000).云南地震資料顯示強震一般發(fā)生在地震的活躍期,兩次活躍期之間也有一段時間的平靜期,例如,1970年通海7.8級地震26a后,麗江發(fā)生7.0級地震;到目前為止,麗江地震發(fā)生已經(jīng)16a,云南地區(qū)至今還未發(fā)生7級以上地震,但6級以上地震逐漸突出.有趨勢表明云南即將或已進入強震活躍時段,未來面臨強震連發(fā)的嚴峻局面,因此對云南地區(qū)的孕震構(gòu)造特征研究以及未來發(fā)震形式的評估顯得尤為重要.

    云南地區(qū)的深大斷裂以南北向的怒江斷裂、瀾滄江斷裂、程海斷裂、元謀—綠汁江斷裂和小江斷裂及北西向的金沙斷裂和紅河斷裂為主.這些大斷裂把云南地區(qū)劃分為騰沖塊體、保山塊體、蘭坪—思茅弧后盆地、鹽源—麗江陸緣坳陷、滇中坳陷、康滇古隆起、滇東坳褶帶等構(gòu)造單元,組成了云南地區(qū)的構(gòu)造格架(蘇有錦等,1999).

    大地電磁測深方法是研究地下深部電性結(jié)構(gòu)的重要方法,其研究結(jié)果可以從電導率或電阻率隨深度的變化推測地球深部巖石圈結(jié)構(gòu)和物質(zhì)組成及狀態(tài),以尋找如斷層帶、高阻巖體或低阻局部熔融體等地球內(nèi)部特殊結(jié)構(gòu),為地球動力學、地質(zhì)災害防治、礦床成因研究等提供重要的支撐(金勝等,2010),為地球深部孕震環(huán)境提供科學依據(jù).

    在青藏高原及東邊緣地殼上地幔的電性結(jié)構(gòu)及地震成因的研究中,各國學者都取得了不少有價值的成果.在川西—藏東地區(qū)進行的大地電磁測深研究中,孫潔等(2003)推斷現(xiàn)今川滇菱形北部塊體處于熱狀態(tài),是近代活動塊體之一;湯吉等(2005)通過對青藏高原東北緣進行的大地電磁測深觀測,得到了該區(qū)地殼深部電性結(jié)構(gòu),趙國澤等(2004)對此提出了該區(qū)強烈構(gòu)造變形的動力學機制;詹艷等(2005)在海原弧形構(gòu)造區(qū)的深部電性結(jié)構(gòu)研究中發(fā)現(xiàn)了其殼內(nèi)低阻帶不連續(xù)發(fā)育,顯示了海原大地震的深部結(jié)構(gòu)和構(gòu)造背景;近年來有關青藏高原及其附近地區(qū)的地殼上地幔的低阻層的研究也得到了很多突破性的進展,魏文博等(2009)在青藏高原巖石圈導電性結(jié)構(gòu)與流變性的研究中討論了藏中、藏北和藏南巖石圈的流變性問題,提出了中下地殼的部分熔融理論;白登海等(2010)在對青藏高原東部大地電磁測深研究中,證明了地殼流模型;Unsworth等(2005)在對青藏高原南緣喜馬拉雅深部電性構(gòu)造研究中也證實大規(guī)模的“地殼流”的存在.這些研究對解釋已經(jīng)發(fā)生的構(gòu)造變形、地震活動、深部成礦及預測未來的地震活動等具有深遠意義.

    孟連—羅平大地電磁測深剖面位于青藏高原東南緣板塊活動帶,穿過云南耿馬—瀾滄、思茅—普洱和通海三個典型的強震活動區(qū)(見圖1).通過在該區(qū)的大地電磁測深研究,得到了云南南部地區(qū)近70km深度的地殼上地幔電性結(jié)構(gòu)特征及相互關系,為該區(qū)的強震孕震環(huán)境和地球動力學解釋提供電性依據(jù).

    2 數(shù)據(jù)采集、處理、分析與反演

    2.1 剖面位置

    大地電磁測深研究是喜馬拉雅項目“中國大陸深部孕震環(huán)境探查”第一期“南北地震帶南段”的研究內(nèi)容之一,布設了兩條大地電磁測深(MTMagnetotelluric)剖面.本研究區(qū)位于中國西南部地區(qū)的云南省境內(nèi),測線沿孟連—思茅—羅平布設,總長近569km.

    圖1剖面(方框中所示)即為本文所闡述的大地電磁測深剖面,該剖面西起孟連,向東穿過云南的主要構(gòu)造單元左貢—耿馬褶皺系、三江褶皺系和華南褶皺系,并跨過瀾滄江斷裂帶、哀牢山—紅河斷裂帶、楚雄斷裂和小江斷裂,經(jīng)過耿馬—瀾滄強震區(qū)、思茅—普洱強震區(qū)和通海強震區(qū).沿剖面布設了114個大地電磁測深點,測點間距平均在5km左右,但是在大型斷裂帶附近及主要強震區(qū)加密.

    2.2 野外觀測及資料處理

    圖1 云南地區(qū)地質(zhì)構(gòu)造與大地電磁測深剖面點位圖I保山地塊,II思茅地塊,III川滇菱形塊體,IV滇東南加里東褶皺帶.F1:怒江斷裂;F2:瀾滄江斷裂;F3:打洛—景洪斷裂;F4:木戛斷裂;F5:無量山斷裂;F6:阿墨江斷裂;F7:哀牢山斷裂;F8:紅河斷裂;F9:綠汁江斷裂;F10:普渡河斷裂;F11:小江斷裂;F12:彌勒—富源斷裂;F13:彌勒—師宗斷裂;F14:楚雄—建水斷裂.Fig.1 Yunnan regional tectonic map and location of MT sites I Baoshan block,II Simao block,III Sichuan-Yunnan rhombus-shape block,IV Southeastern Yunnan caledon fold belt.F1:Nujiang fault;F2:Lancangjiang fault;F3:Daluo-Jinghong falut;F4:Mujia falut;F5:Wuliangshan falut;F6:Amojiang fault;F7:Ailaoshan fault;F8:Red River fault;F9:Lüzhijiang fault;F10:Pudu river fault;F11:Xiangjiang fault;F12:Mile-Fuyuan fault;F13:Mile-Shizong fault;F14:Chuxiong-Jianshui fault.

    野外工作共完成114個有效大地電磁測深點的測量,野外測量采用了鳳凰公司生產(chǎn)的MTU-5A型大地電磁測深儀器采集的常規(guī)和長周期大地電磁測深資料,以及烏克蘭生產(chǎn)的超長周期大地電磁測深儀采集的超長周期大地電磁測深資料,其頻帶為0.00002~320Hz(測點詳情見圖1所示).在數(shù)據(jù)采集過程中,在離測線500km以外電磁干擾很小的地方布設了遠參考站.所有常規(guī)MT測點的野外記錄時間不小于20h,長周期MT測點的觀測時間不小于40h,超長周期的記錄時間不小于10天,但在后續(xù)的資料處理和反演解釋過程中,僅選取了常規(guī)周期的大地電磁測點資料進行處理.

    野外觀測資料為原始的時間序列數(shù)據(jù),以便資料處理時可以剔除受人文噪音干擾的信號.每個測點觀測完成后馬上進行現(xiàn)場的初步資料處理,若數(shù)據(jù)質(zhì)量不合格,立即進行復測或改變位置重新測量.資料處理采用鳳凰公司提供的SSMT2000數(shù)據(jù)處理軟件,觀測到的時間序列數(shù)據(jù)經(jīng)過傅氏變換,得到電磁場的自、互功率譜后,采用遠參考Robust資料處理技術(Egber and Booker,1986)估算大地電磁張量阻抗,以提高數(shù)據(jù)處理的質(zhì)量;再對各頻率的功率譜進行資料挑選,剔除因干擾噪音引起的飛點資料,經(jīng)編輯處理后得到各觀測點的視電阻率和阻抗相位、主軸方位角、二維偏離度和傾子等資料的數(shù)據(jù).為消除局部三維異常體對視電阻率和阻抗相位帶來的畸變影響,采用了Bahr張量阻抗分解技術(Karsten,1990)對114個常規(guī)的大地電磁測深點進行了張量阻抗分解.

    2.3 曲線分析

    每個測點的視電阻率和阻抗相位曲線反映了該測點地下電性隨深度的變化和該地區(qū)測點及附近地下介質(zhì)的電性不均勻性特征.通過對全部大地電磁測深點的視電阻率曲線和相位曲線分析發(fā)現(xiàn),沿剖面的曲線形態(tài)變化與地塊過渡具有很好的對應關系,不同地塊中的視電阻率和阻抗相位具有不同的曲線形態(tài),同一地塊的不同測點曲線具有相似性.圖2顯示了從南西至北東分布在保山地塊(001-012)、思茅地塊(013-042)、川滇菱形塊體(043-078)和滇東坳褶帶(079-114)的12個典型測點的視電阻率(ρs)和阻抗相位(φ)曲線,它們反映了不同區(qū)塊的電性結(jié)構(gòu)特征.

    圖2 大地電磁剖面典型測點視電阻率(Ωm)與相位(°)曲線圖Fig.2 Apparent resistivity(Ωm)and phase(°)curve of typical MT sites

    圖3 二維偏離度(Bahr張量分解)Fig.3 Two-dimensional skewness(Bahr-impedance tensor decomposition)

    在圖2中,003和010測點屬于保山地塊,視電阻率和相位曲線形態(tài)大致相同,視電阻率先隨深度的增加(頻率減?。┒鴾p小—增大—減小,在深部又有增大的趨勢,初步分析可能存在殼內(nèi)的相對低阻層;020、029和039測點位于思茅地塊,曲線結(jié)構(gòu)整體上有相似的特點,地下介質(zhì)的視電視阻率大致隨深度增加而減小,但在細節(jié)上還是存在一些差異,初步判斷殼內(nèi)深部存在低阻體,不同測點的視電阻率曲線形態(tài)的區(qū)別又可能是由于殼內(nèi)的低阻體規(guī)模存在差異;051、063、071號測點位于川滇菱形塊體上,因該塊體地質(zhì)結(jié)構(gòu)上就比較復雜,多條斷裂發(fā)育,在視電阻率曲線形態(tài)上表現(xiàn)為殼內(nèi)有相對高阻體的存在,在深部電阻率減??;082、092、105和109號測點屬于滇東坳褶帶,在圖2中可以看到該區(qū)測點的視電阻率曲線較其他地塊相對復雜,曲線形態(tài)波動較多,該區(qū)深部結(jié)構(gòu)可能也會相對復雜一些.

    2.4 區(qū)域維數(shù)分析和電性構(gòu)造走向

    大地電磁測深阻抗的二維偏離度是反應地下介質(zhì)電性結(jié)構(gòu)維數(shù)特征的參數(shù),它的取值范圍可判別地下介質(zhì)的維數(shù)特征,以了解地下介質(zhì)電性結(jié)構(gòu)的復雜程度.二維偏離度數(shù)值越小,說明地下介質(zhì)越趨于二維結(jié)構(gòu),在深部探測中,當二維偏離度小于0.3時,二維解釋是可以接受的.采用Bahr技術對114個測點進行了張量阻抗分解,結(jié)果(見圖3)顯示研究區(qū)內(nèi)大部分測點分解前后的視電阻率和阻抗相位變化不大,整體受畸變影響很小,大部分測點的整體二維偏離度很?。╯<0.3),可近似為二維結(jié)構(gòu).

    對114個測點的不同頻段的區(qū)域電性主軸方位角進行統(tǒng)計分析(圖4),發(fā)現(xiàn)在低頻段(0.0001~1Hz)、中頻段(1~10Hz)和高頻段(10~320Hz)的電性構(gòu)造主軸方位角在不同塊體及斷裂帶上有些波動,但除局部測點,總體上主軸方位角的變化在120°~160°之間.結(jié)合該區(qū)地質(zhì)和斷裂構(gòu)造走向多為北西30°,以及剖面方向北東60°的特點,判斷地下的電性構(gòu)造走向基本與剖面垂直,因此,將阻抗旋轉(zhuǎn)到測線方向(-30°)進行二維反演.研究區(qū)橫跨云南南部不同的區(qū)塊及深大斷裂,其橫向上結(jié)構(gòu)復雜,對比該剖面的TE模式和TM模式的反演結(jié)果,我們發(fā)現(xiàn)TM模式的反演結(jié)果在橫向上更直觀地反映地下實際結(jié)構(gòu),因此最終選擇TM模式對這條剖面全部測點進行二維非線性共軛梯度法反演.

    2.5 反演結(jié)果

    目前,在大地電磁測深反演運算中,比較常用的MT二維反演法主要有快速松弛(RRI)反演(Smith and Booker,1996)、Occam反演(Constable etal.,1987)和非線性共軛梯度反演(NLCG)(Rodi and Mackie,2001)等方法.在本研究中,二維反演選用TM模式,采用由陳小斌開發(fā)的MT-pioneer反演軟件(蔡軍濤等,2010),應用二維非線性共軛梯度法(NLCG),對沿剖面各測點觀測的視電阻率和相位進行了反演.

    將所有的觀測資料旋轉(zhuǎn)-30°,即到剖面方向(北東60°),這一方向基本與前面分析的區(qū)域電性構(gòu)造走向垂直.因此,將沿剖面方向的視電阻率和阻抗相位曲線規(guī)定為TM曲線(H極化),垂直剖面方向的曲線規(guī)定為TE曲線(E極化).二維NLCG反演中所選取的不均勻模型網(wǎng)格,在水平方向上遵循一個網(wǎng)格內(nèi)只能有一個測點,但相鄰兩個測點之間有1~3個網(wǎng)格的原則,因測點間距為5~10km,所以網(wǎng)格長度大約在1~3km,這樣的網(wǎng)格選取既不影響數(shù)據(jù)的擬合度,又能更好地反映地下的復雜結(jié)構(gòu);在垂直方向上,采用隨深度增加逐漸加大的方式選取網(wǎng)格間距.

    圖4 典型測點的主軸方位角(圖中數(shù)字單位為°,最下面一行數(shù)字為測點號)(a)0.00001~0.1Hz;(b)0.1~10Hz;(c)10~320Hz.Fig.4 Electrical regional structure strike angle of typical MT sites

    反演的理想效果是能夠使得反演模型較好地擬合數(shù)據(jù),同時又能得到光滑的模型.反演目標函數(shù)中的正則化因子影響了反演目標函數(shù)的擬合方向,即正則化因子越大,所得模型越平滑,正則化因子越小,數(shù)據(jù)擬合越好,但反演模型粗糙.L曲線分析(Hansen etal,1992)是選取最佳反演正則化因子的一種有效方法,本文選取1—500區(qū)間內(nèi)共13個正則化因子(500,300,200,100,80,70,60,50,40,30,20,10,1)進行多次反演,得到如圖5所示的L曲線,其縱軸為擬合均方差(RMS),表示數(shù)據(jù)擬合程度,橫軸為模型約束目標函數(shù)的平方根,評價模型光滑程度.從L曲線中看出,曲線的拐點對應正則化因子為50,再通過對比不同正則化因子的反演結(jié)果,最終選取τ=50的反演模型,此時的反演結(jié)果既保證了模型的光滑度,又能較好地擬合數(shù)據(jù)(RMS=3.62).

    圖5 不同反演參數(shù)下的正則化因子曲線(圖中數(shù)字為正則化因子)Fig.5 L-curve for regularization factor analysis in different inversion factors

    3 電性結(jié)構(gòu)分析

    根據(jù)二維反演處理的電性結(jié)構(gòu)以及視電阻率和阻抗相位曲線的特點,結(jié)合地表地質(zhì)和斷裂出露的情況以及活動地塊的劃分,繪制了電性構(gòu)造解釋圖(見圖6).其中,白色虛線為推斷的斷裂,由地震觀測結(jié)果推測莫霍面深度為40~50km(張中杰等,2005).總體來看,電性結(jié)構(gòu)表現(xiàn)出上地殼(0~20km)構(gòu)造比較復雜,低阻層和高阻層相互交錯,下地殼(20~40km)橫向上發(fā)育有不連續(xù)的低阻層,其中規(guī)模較大地區(qū)如圖6中A、B、C所示,而沿剖面的深部電性結(jié)構(gòu)特征顯示與區(qū)域地質(zhì)構(gòu)造有較好的對應關系.根據(jù)地質(zhì)構(gòu)造的地塊劃分來看,每個區(qū)塊的地下電性結(jié)構(gòu)都有不同程度的差異.

    反演模型見圖6b所示,其中藍色表示相對高電阻率,紅色表示相對低電阻率.由反演模型可以看出,沿測線地下介質(zhì)的電性結(jié)構(gòu)表現(xiàn)為橫向分塊,縱向分層的特點.

    圖7所示為MT剖面視電阻率與阻抗相位的觀測值與反演結(jié)果的擬合情況,由圖中可見視電阻率曲線及阻抗相位的擬合情況良好,表明二維反演結(jié)果具有較高的可信度.

    3.1 地塊電性特征

    二維電性結(jié)構(gòu)圖直觀地表現(xiàn)出電性結(jié)構(gòu)的分塊特征,自西南向東北分別為保山地塊(區(qū)塊I)思茅地塊(區(qū)塊II)、川滇菱形塊體(區(qū)塊III)、滇東坳褶帶(區(qū)塊IV),他們之間的邊界分別為瀾滄江斷裂帶、哀牢山—紅河斷裂帶以及小江斷裂.

    3.1.1 保山地塊

    保山地塊(001—012)位于瀾滄江斷裂南支以西,屬滇西南地區(qū).耿馬—瀾滄強震區(qū)位于研究區(qū)的瀾滄—勐遮斷裂帶與打洛—景洪斷裂帶交接處(云南省地質(zhì)礦產(chǎn)局,1990)(見圖1),該區(qū)地震頻發(fā),強震突出,地表發(fā)育有多條橫縱交切的活動斷裂帶,主要包括瀾滄江斷裂帶、木戛斷裂帶、打洛—景洪斷裂帶、瀾滄—勐遮斷裂帶等等.

    由資料處理得到的全部測點的視電阻率曲線表明,該區(qū)塊001—009號測點的視電阻率曲線形態(tài)基本一致,只有在靠近瀾滄江斷裂帶的測點曲線形態(tài)發(fā)生變化.圖2中所示的003號測點,兩個極化方式的視電阻率曲線基本重合.

    在二維反演結(jié)果圖(見圖6b)中顯示,地塊中電性結(jié)構(gòu)相對簡單,008號測點以西地區(qū)從地表到約5km深度,發(fā)育一相對高阻層,電阻率總體在100~300Ωm;5~15km的上地殼發(fā)育有一低阻層,電阻率總體在10Ωm以內(nèi),15km以下是高阻.而009—012號測點的地下0~25km深度卻發(fā)育有一規(guī)模較大的高阻區(qū),電阻率高達10000Ωm以上,高阻層以下表現(xiàn)為相對低阻.殼內(nèi)低阻層在構(gòu)造運動中起著重要的作用,它被普遍認為是塑性的軟弱層,活動性與可塑性強.耿馬—瀾滄強震區(qū)位于地殼低阻薄層與高阻體的交界處,因其西南部受到印度板塊擠壓,東部受到思茅塊體及川滇菱形塊體的阻擋,上部地殼便沿低阻層滑動或逆掩,形成了特殊的構(gòu)造帶(孫潔等,1989),導致該區(qū)的構(gòu)造活動十分強烈,地形受到強烈切割.保山地塊中的木戛斷裂、打洛—景洪斷裂和瀾滄—勐遮斷裂的交接地帶都曾發(fā)生過5.5級以上地震,因此推測該區(qū)地震多發(fā)與地殼中電性結(jié)構(gòu)的不均勻性以及多條斷裂帶交錯及相互作用有密不可分的聯(lián)系.

    3.1.2 思茅地塊

    思茅地塊(013—042)位于瀾滄江斷裂帶及哀牢山—紅河斷裂帶之間,隸屬唐古拉—昌都—蘭坪—思茅褶皺系,該地塊發(fā)育有一條較大的北北西走向的無量山斷裂帶,以及位于其東側(cè)的把邊江斷裂和阿墨江斷裂,把該地塊分成東西兩個部分.

    圖6 孟連—羅平剖面測點地質(zhì)構(gòu)造(a)及地殼上地幔電性結(jié)構(gòu)(b)解釋圖(上面為沿剖面各測點點位示意圖)白色實線、虛線:斷裂;黑色圓圈:強震區(qū)震源深度.LMF:瀾滄—勐遮斷裂;LCJF:瀾滄江斷裂;WLSF:無量山斷裂;ALSF:哀牢山斷裂;HHF:紅河斷裂;CJF:楚雄—建水斷裂;CTF:楚雄—通海斷裂;XJF:小江斷裂;MFF:彌勒—富源斷裂;MSF:彌勒—師宗斷裂.Fig.6 Electrical conductivity structure of crust and upper mantle and electrical tectonic analysis along the MT profile Menglian-Luoping(The diagram of sites along the profile for the above)Real and dashed line in White:Fault;Circle in black:Depth of strong earthquakes.LMF:Lancang-Mengzhe fault.LCJF:Langcangjiang fault;WLSF:Wuliangshan fault;ALSF:Ailaoshan fault;HHF:Red river fault;CJF:Chuxiong-Jianshui fault;CTF:Chuxiong-Tonghai fault;XJF:Xiaojiang fault;MFF:Mile-Fuyuan fault;MSF:Mile-Shizong fault.

    圖7 二維反演模型視電阻率、相位擬合斷面圖(a1)觀測的視電阻率;(a2)模型計算得到的視電阻率;(b1)觀測的阻抗相位;(b2)模型計算得到的阻抗相位.Fig.7 Apparent resistivity and phase curve pseudosection map of two-dimension inversion model(a1)Observation apparent resistivity;(a2)Computation apparent resistivity;(b1)Observation impedance;(b2)Computation impendace.

    從各測點的視電阻率曲線來看,整體上曲線的構(gòu)造形式有較大變化,曲線形態(tài)也直觀地顯示,在無量山斷裂附近地下介質(zhì)的高阻特點及東西兩側(cè)地下介質(zhì)中存在規(guī)模較大的低阻層.

    通過二維反演結(jié)果顯示,思茅地塊深部電性結(jié)構(gòu)表現(xiàn)為上下分層,橫向不連續(xù)的特點.無量山斷裂帶附近0~15km發(fā)育有一規(guī)模較大的高阻體,電阻率大約在10000Ωm以上,該高阻體將思茅地體分為東西兩部分,這兩部分在結(jié)構(gòu)上有相同點,也存在差異,如圖6b中A、B兩紅色區(qū)域,電阻率總體在1~100Ωm之間,它們規(guī)模都較大,但是在埋深和規(guī)模上有所區(qū)別.

    低阻體A的埋深在15~40km之間,東西跨度約60km,平均電阻率大約在10Ωm以內(nèi),其邊界有向兩側(cè)上地殼延伸的趨勢.在低阻體A以上的0~15km深度發(fā)育有一高阻蓋層,電阻率高達1000~10000Ωm.由該區(qū)地表區(qū)域地質(zhì)資料顯示的火山巖夾層出露情況來看(見圖6b),其地質(zhì)歷史時期火山活動較強烈,從普洱地區(qū)附近的溫泉分布情況來看,也說明該區(qū)地殼內(nèi)的地熱活動仍十分活躍.低阻體B的埋深大約在7~30km之間,東西跨度約60km,平均電阻率大約在10Ωm以內(nèi),形態(tài)表現(xiàn)為中間薄兩邊厚的蝶狀,其上0~10km發(fā)育兩個相對高阻層,分布于以阿墨江斷裂為界的兩側(cè).無量山斷裂帶是一條重要的地震活動斷裂帶,歷史上沿該斷裂帶曾發(fā)生多次MS6.0級以上強震,思茅—普洱強震區(qū)就位于無量山西麓,結(jié)合該區(qū)地質(zhì)資料和地熱資料等推斷其地震成因與深部構(gòu)造運動和熱驅(qū)動機制有關.

    3.1.3 川滇菱形塊體

    處于紅河斷裂以東,小江斷裂以西的云南地區(qū)(043—078)屬于川滇菱形塊體南端.該區(qū)域地況復雜,斷裂互相交錯,有北北西走向為主的楚雄—建水斷裂,南北走向的普渡河斷裂、小江斷裂、彌勒—富源斷裂等深大斷裂發(fā)育(見圖1).

    大地電磁測深結(jié)果顯示,該區(qū)地下5~50km的地殼上地幔低阻體規(guī)模較大(如圖6b低阻體C),電阻率總體在100Ωm以內(nèi),大部分電阻率值為1~10Ωm.該低阻體向東延伸部分離地表較近,且電阻率值略高于西側(cè).低阻體C的西側(cè)靠近地表附近有一組破碎低阻薄層,初步推斷其成因為紅河斷裂帶在歷史上經(jīng)歷從左旋水平位移轉(zhuǎn)變?yōu)橛倚呋?、伴隨不等量的垂直運動的過程中,地表巖石的大規(guī)模破碎及上下物質(zhì)的交換,實際上紅河斷裂帶東南段在現(xiàn)代構(gòu)造活動中的運動已大為減弱(柴天俊和楊繼武,1992).低阻體C以上的0~15km深度內(nèi)發(fā)育橫向上不連續(xù)的相對高阻塊體,電阻率約為100~1000Ωm.綜合該地塊地殼上地幔的電性結(jié)構(gòu)特征,結(jié)合其區(qū)域地質(zhì)構(gòu)造和板塊劃分,若向西北方向追蹤,我們發(fā)現(xiàn)該低阻體與羌塘地體中下地殼的管流體(Bai etal.,2010;葉高峰等,2007;Wei etal.,2001;Unsworth etal.,2004),以及川滇菱形塊體中下地殼的弱物質(zhì)層相對應(白登海等,2011),認為這是青藏高原地殼低阻管流通道的西南出口.本文中低阻體C的埋深與前人的研究結(jié)果存在一定的差異,是由于本文在川滇菱形塊體的大地電磁剖面測點的點距為2~5km,在二維反演時數(shù)據(jù)約束條件更多,反演得到的電性結(jié)構(gòu)模型反映的殼內(nèi)低阻層特征可能更加精細.

    3.1.4 滇東坳褶帶

    滇東坳褶帶是華南加里東褶皺系的西段,位于小江斷裂帶東南地區(qū)(079—114),區(qū)域地質(zhì)圖中可以看到地表出露多條錯綜復雜不同走向的斷裂(見圖1).

    大地電磁測深結(jié)果表明,該區(qū)中上地殼由一系列復雜的高阻體和低阻體組成(見圖6b),尤其在師宗地區(qū),地下5~10km處發(fā)育有一相對高阻體,電阻率約1000Ωm,周圍是放射狀的低阻體,電阻率約為10~100Ωm,電性結(jié)構(gòu)十分復雜.區(qū)域地質(zhì)資料顯示,與測線方向近平行的羅平、師宗、瀘西、朋普等地區(qū)屬于羅平—師宗斷褶束,地質(zhì)構(gòu)造以北東向斷裂為主,表現(xiàn)為斷裂密集的斷褶束特點.彌勒—師宗斷裂是該區(qū)的一條深大斷裂,由一系列的傾向NW、傾角40°~60°的逆沖斷層組成(董云鵬等,1999),盡管本研究所在的測線與彌勒—師宗斷裂的夾角很?。?5°左右),但從電性構(gòu)造主軸來看,其構(gòu)造走向與剖面西部基本一致,得到的地下介質(zhì)的電性結(jié)構(gòu)仍對褶皺帶地區(qū)有直觀的顯示.深部電性結(jié)構(gòu)表明上地殼區(qū)域發(fā)育有一相對高阻層,下覆的相對低阻體(100Ωm左右)向上地殼延伸,將高阻層以彌勒—師宗斷裂為界在其東西兩側(cè)分為兩個小的高阻體.在區(qū)域地質(zhì)構(gòu)造圖中能夠看出,該區(qū)與鄰區(qū)川滇菱形塊體出露的石炭系、二疊系的沉積相比,此處發(fā)育的三疊系沉積與川滇菱形塊體有較大區(qū)別(見圖6b).實際上,該區(qū)屬于穩(wěn)定的華南地體,在最近一個世紀內(nèi)并沒有較大的地震發(fā)生,地質(zhì)資料顯示滇東坳褶帶大約形成于古生代早期地殼運動的加里東運動,現(xiàn)代活動減弱.

    3.2 邊界斷裂帶分布

    3.2.1 瀾滄江斷裂帶

    瀾滄江斷裂總體近南北向延伸,呈北西向弧形彎曲,境內(nèi)長達800km,是青藏高原東南緣重要的地塊邊界斷裂帶.瀾滄江斷裂帶總體向東傾斜,而研究區(qū)所在區(qū)域為瀾滄江斷裂帶南段,斷裂帶在該段的傾角大致呈垂直向(鐘康惠等,2004).地質(zhì)物探資料顯示斷裂帶控制了巖漿活動和變質(zhì)作用,說明瀾滄江是一條深大斷裂帶.瀾滄江斷裂的走滑逆沖推覆構(gòu)造得到各界學者的一致共識,地震歷史資料和現(xiàn)代儀器記錄都表明,瀾滄江斷裂帶上地震活動比較弱,整條斷裂帶上僅僅在鳳慶地區(qū)發(fā)生過一次5.5級破壞性地震(李祥根,2010).地球深部電性結(jié)構(gòu)顯示,瀾滄江斷裂的上地殼(0~20km)中發(fā)育有規(guī)模較大的高阻體,電阻率高達10000~100000Ωm,高阻體在斷裂帶附近被阻斷,斷裂帶兩側(cè)地塊中都發(fā)育有連續(xù)性的規(guī)模較大的低阻層,并且都有向該斷裂延伸減薄的趨勢,從電性結(jié)構(gòu)特點上也可以看出瀾滄江斷裂是一條深大斷裂.

    3.2.2 哀牢山—紅河斷裂帶

    哀牢山—紅河斷裂帶是我國西南地區(qū)地質(zhì)歷史時期最為重要的一級構(gòu)造單元邊界斷裂帶之一,北起青藏高原東部,斜貫云南及越南北部,向東南延伸進入南海海域,總體呈北西—北北西走向,長達1000多公里,整體呈NE傾向,傾角較大,總體在50°~75°之間(柴天俊等,1992).該斷裂帶實為兩條近平行的深大斷裂的南段,分別為程海—哀牢山大斷裂及金沙江—紅河大斷裂(鄭慶鰲等,2006).

    根據(jù)地質(zhì)歷史資料顯示,哀牢山—紅河斷裂帶在現(xiàn)代地質(zhì)時期處于弱邊界作用(吳海威和楊繼武,1989).大地電磁測深結(jié)果顯示,哀牢山斷裂和紅河斷裂分別位于近地表高阻帶的東西兩端,均是低阻層與高阻體的交界,高阻層的電阻率大約為100~10000Ωm,兩側(cè)低阻層電阻率大約為1~100Ωm.區(qū)域地質(zhì)資料顯示哀牢山—紅河斷裂帶出露由一系列片巖、變粒巖、片麻巖、花崗質(zhì)混合巖及大理巖組成的深變質(zhì)巖系,而這種特殊的電性結(jié)構(gòu)也為該區(qū)地下地質(zhì)構(gòu)造的復雜性及其地質(zhì)歷史時期的陸塊碰撞作用提供了地電學的依據(jù).

    3.2.3 小江斷裂

    小江斷裂作為川滇菱形塊體與華南地體的邊界帶,備受地質(zhì)學家關注.根據(jù)地質(zhì)資料顯示,小江斷裂在加里東、華力西、印支、燕山和喜馬拉雅運動期一直都有活動,而且現(xiàn)今仍較活躍(云南省地質(zhì)礦產(chǎn)局,1990),已有的地震資料顯示,小江斷裂是一條具有強破壞性的地震頻發(fā)的活動斷裂帶(錢曉東和秦嘉政,2008).斷裂帶呈南北走向,全長約450km,在東川以南地區(qū)分成東西兩分支.深部電性結(jié)構(gòu)特征為大規(guī)模低阻體的發(fā)育,電阻率約10Ωm左右,深度為3~20km,上覆0~3km的相對高阻層,電阻率100~1000Ωm.小江斷裂在現(xiàn)代活動方式以左旋走滑為主,在印度板塊向北楔入及青藏高原隆起的過程中,斷裂帶西部的川滇菱形塊體向南東方向擠出,地球深部發(fā)育的上地殼低阻體及西側(cè)塊體中下地殼低阻體為這種理論提供了電性依據(jù).

    3.3 強震區(qū)電性特征

    大地電磁測深剖面穿過云南南部地區(qū)的三個強震區(qū):耿馬—瀾滄強震區(qū)、思茅—普洱強震區(qū)和通海強震區(qū),在地電模型中(見圖6b),我們可以觀察到它們的震源在埋深上有一個共同特點,就是都處于殼內(nèi)低阻體與高阻體交界變化的電性梯度帶上.這個現(xiàn)象表明,強震的發(fā)生,與殼內(nèi)低阻體的存在以及其與高阻體位置關系有密不可分的聯(lián)系.

    3.3.1 耿馬—瀾滄強震區(qū)

    該強震區(qū)位于保山地塊,由電性結(jié)構(gòu)模型(見圖6b)可以看出,其震源多位于地殼低阻薄層與高阻體的交界處,因保山地塊西南部受到印度板塊擠壓,東部受到思茅塊體及川滇菱形塊體的阻擋,上部地殼便沿低阻層滑動或逆掩,形成了特殊的構(gòu)造帶,導致該區(qū)的構(gòu)造活動十分強烈,地形受到強烈切割.保山地塊中的木戛斷裂、打洛—景洪斷裂和瀾滄—勐遮斷裂的交接地帶都曾發(fā)生過5.5級以上地震,因此推測該區(qū)地震多發(fā)與地殼中的電性結(jié)構(gòu)不均勻性以及多條斷裂帶交錯及相互作用有關.

    3.3.2 思茅—普洱強震區(qū)

    該強震區(qū)位于思茅地塊內(nèi),地處無量山西麓.無量山斷裂帶是一條重要的地震活動斷裂帶,歷史上沿該斷裂帶也曾發(fā)生多次MS6.0級以上強震,如圖6b中所示,思茅—普洱強震區(qū)處于低阻體A與無量山斷裂帶下高阻體的交界處,低阻體的電阻率為1~10Ωm,推測該低阻體可能與局部熔融有關.地表地質(zhì)資料顯示,該區(qū)火山巖出露廣泛,地質(zhì)歷史時期火山活動較強烈,普洱地區(qū)及附近溫泉廣泛分布,也說明該區(qū)地殼內(nèi)的地熱活動仍十分活躍.推測該地震區(qū)的強震活動可能與深部構(gòu)造運動和熱驅(qū)動機制有密切關系.結(jié)合該區(qū)地球化學結(jié)果和區(qū)域地質(zhì)資料及大地熱流值(100~150mV·m-2)(Hu etal.,2000),推測該地震多發(fā)區(qū)強震成因可能由深部構(gòu)造運動和熱驅(qū)動機制有密切關系.

    3.3.3 通海強震區(qū)

    通海強震區(qū)位于川滇菱形塊體南段,紅河—哀牢山斷裂帶以東,深部電性結(jié)構(gòu)模型(見圖6b)顯示,其震源區(qū)位于地殼低阻體C的東西兩部分交界處,上覆一系列規(guī)模較小的高阻體,電阻率約為100~1000Ωm.距離該強震區(qū)附近發(fā)育有小江斷裂這一活動的斷裂帶,近代地震記錄也表明,沿斷裂帶上地震頻繁的發(fā)生,我們有理由推測該區(qū)的強烈的地殼活動性與活動斷裂帶的發(fā)育及地球深部低阻熱物質(zhì)的流動等構(gòu)造運動有直接的聯(lián)系.

    4 結(jié)論和討論

    沿剖面分布的三個強震帶分別位于保山地塊、思茅地塊和川滇菱形塊體,從圖6b中可以看出,耿馬—瀾滄強震區(qū)、思茅—普洱強震區(qū)和通海強震區(qū)地下深部電性結(jié)構(gòu)都位于低阻層與高阻體的交界處,且附近斷裂密布,斷裂活動頻繁,這些都表明地殼的活動性是導致該地區(qū)發(fā)生大地震的重要因素.而相比之下,地震空白區(qū)在深部電性結(jié)構(gòu)上的特點是無大規(guī)模低阻體與高阻體的接觸,斷裂的活動性相對較弱.如剖面穿過的瀾滄江斷裂帶南段,該斷裂帶雖是一條深大斷裂,但多種資料表明,其近代活動性較弱;哀牢山—紅河斷裂帶也是一條深大斷裂,且在歷史構(gòu)造運動中,作為重要邊界帶,起到影響對該區(qū)地層發(fā)育、巖漿活動等的控制作用,但在現(xiàn)代構(gòu)造運動中已轉(zhuǎn)為弱邊界作用;滇東坳褶帶在現(xiàn)今的構(gòu)造中基本處于平靜的狀態(tài),雖然在地下深部電性結(jié)構(gòu)上顯示該區(qū)具有局部低阻體與高阻體相互交錯的特點,但因為其位于穩(wěn)定的華南地塊,目前并無大地震發(fā)生的可能.

    本研究開展的大地電磁測深剖面所穿過的三個強震區(qū),其地殼中普遍存在大規(guī)模的低阻異常體,那么是否就可以認為低阻異常是導致其強震發(fā)生的原因.一般認為,造成地殼導電率升高的主要因素包括地下導電體(如金屬礦體、石墨等)、地下流體、熔融體等.而在青藏高原及其周圍地區(qū)進行的大地電磁測深研究(趙國澤等,2004,2008;劉美等,2010;Bai etal,2010;葉高峰等,2007;魏文博等,2009;湯吉等,2005)表明,在青藏高原內(nèi)部以及其東邊緣殼內(nèi)存在連續(xù)的“管流”通道,即具有低黏滯度的部分熔融和含鹽流體等軟物質(zhì)流,其特點是容易流動.正是這些容易流動層起著能量傳遞的作用,使該區(qū)地殼內(nèi)保持活動性,也是強震孕震的條件.

    總體來說,通過本研究大致得到以下幾點結(jié)論:

    (1)沿孟連—羅平剖面的大地電磁測深反演結(jié)果得到了一個令人比較滿意的深部地電結(jié)構(gòu)模型,直觀地顯示了沿剖面上各測點地下的電性結(jié)構(gòu),清晰地反映出高阻體和低阻體的位置,大致與該區(qū)地表的區(qū)域地質(zhì)構(gòu)造相對應.

    (2)測線通過云南南部地區(qū)的三個強震區(qū),在深部電性結(jié)構(gòu)上存在著一個共同的特點,就是震源區(qū)附近的殼內(nèi)低阻體和局部高阻體的普遍存在,地震大多發(fā)生在高低阻變化的梯度帶上,這種對應關系在其他地震區(qū)也普遍存在.

    (3)根據(jù)青藏高原以及周圍地區(qū)的大地電磁測深的研究成果,發(fā)現(xiàn)在羌塘地體及青藏高原東緣中下地殼也有與川滇菱形塊體深部的低阻異常體相似的異常,認為是青藏高原殼內(nèi)熱物質(zhì)向其東南方向移動的管流通道的延伸.殼內(nèi)低阻體能量的傳遞和高阻體能量的積累作用有可能使云南南部地區(qū)地塊內(nèi)部斷裂維持活動狀態(tài),為該區(qū)三個強震區(qū)的孕震環(huán)境提供理論模型.

    對于研究區(qū)的部分測點地下電性結(jié)構(gòu)傾向局部三維構(gòu)造,二維反演在某種程度上不能滿足其電性結(jié)構(gòu)的全部要求,在該區(qū)開展三維反演的工作應是必不可少的,但由于某種原因本文不能一一闡述,在未來的研究中應加大力度對大地電磁測深三維反演工作進行更深一步的研究.

    致謝 本文在計算中應用了陳小斌研究員的大地電磁反演軟件,江漢石油管理局物探公司在野外資料采集中付出了辛勤的勞動,以及本文的評閱人為本文提出了非常有益的修改建議,在此作者對中國地震局地質(zhì)研究所電磁組的全體成員予以的幫助和指導一并表示衷心的感謝.

    Bai D H,Teng J W,Ma X B,etal.2011.Crustal flow beneath eastern Tibetan Plateau revealed by magnetotelluric measurements.China Basic Scienc(in Chinese),1:7-10,doi:10.3969/j.issn.1009-2412.

    Bai D H,Unsworth M J,Meju M A,etal.2010.Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging.Nature Geoscience,3(5):358-362,doi:10.1038/NGEO830.

    Bareau of Geology and Mineral Resources of Yunnan Province.1990.Regional Geology of Yunnan Province(in Chinese).Geological Publishing House.

    Cai J T,Chen X B,Zhao G Z.2010.Refined techniques for data processing and two dimensional inversion in magnetotelluric I:tensor decomposition and dimensionality analysis.Chinese J.Geophys.(in Chinese),53(10):2516-2526,doi:10.3969/j.issn.0001-5733.

    Chai T J,Yang J W.1992.Geological and geomorphological features of neotectonic movement of Red River fault zone.Yunnan Geology(in Chinese),11(3):260-267.

    Constable S C,Parker R L,Constable C G.1987.Occam′s inversion:apractical algorithm for generating smooth models from electromagnetic sounding data.Geophysics,52(3):289-300.

    Dong Y P,Zhu B Q.1999.The island arc-pillow type lava in Jianshui,southeastern Yunnan Province and its constraint on paleotethys of southern part of China.Chinese Science Bulletin(in Chinese),44(21):2323-2328.

    Egber G D,Booker J R.1986.Robust estimation of geomagnetic transfer functions.Geophys.J.Int.,87(1):173-194.

    Hansen P C.1992.Analysis of discrete ill-posed problems by means of the L-Curve.SIAM Review,34(4):561-580.

    Hu S B,He L J,Wang J Y.2000.Heat flow in the continental area of China:a new data set.Earth and Planetary Science Letters,179(2):407-419,doi:10.1016/S0012-821X(00)00126-6.

    Huang P G,Chen Y,Qin J Z,etal.2010.The Seismicity in Yunnan(in Chinese).Kunming:Yunnan Science &Technology Press.

    Huang P G,Shi S X,Su Y J.2000.Study on seismicity in Yunnan in the 20th century.Journal of Seismicity Research(in Chinese),23(1):2-9.

    Jin S,Zhang L T,Wei W B,etal.2010.Magnetotelluric method for deep detection of Chinese continent.Acta Geologica Sinica(in Chinese),84(6):808-817.

    Karsten B.1990.Geological noise in magnetotelluric data:a classification of distortion types.Phys.Earth Planet.,66(1-2):24-38.

    Li X G.2010.The Seismotectonics Movement of China(in Chinese).Seismological Press.

    Liu M,Bai D H,Xiao P F.2010.The electrical conductivity structure of the eastern Tibetan Plateau and its tectonic implications.Seismology and Geology(in Chinese),32(1):51-58,doi:10.3969/j.issn.0253-4967.

    Qian X D,Qin J Z.2008.Strong earthquake risk analysis of Xiaojiang fault zone and surrounding area.Journal of Seismological Researh(in Chinese),31(4):354-361.

    Rodi W,Mackie R L.2001.Nonlinear conjugate gradients algorithm for 2Dmagnetotelluric inversion.Geophysics,66(1):174-187,doi:10.1190/1.1444893.

    Smith J T,Booker J R.1996.Rapid inversion of two and three dimensional magnetotelluric data.J.Geophys.Res.,96(B3):3905-3922,doi:10.1029/90JB02416.

    Su Y J,Liu Z Y,Cai M J,etal.1999.Deep earth medium background of the strong earthquakes distribution in Yunnan region.Acta Seismologica Sinica(in Chinese),21(3):313-322.

    Sun J,Jin G W,Bai D H,etal.2003.Crustal and upper mantle eletrical structure of eastern Tibetan Plateau and its tectonic significance.Sci China Ser D-Earth Sci(in Chinese),33(Suppl.):173-180.

    Sun J,Xu C F,Jiang Z,etal.1989.The electrical structure of the crust and upper mantle in the west part of Yunnan Province and its relation to crustal tectonics.Seismology and Geology(in Chinese),11(1):1-11.

    Tang J,Zhan Y,Zhao G Z,etal.2005.Electrical conductivity structure of the crust and upper mantle in the northeastern margin of the Qinghai-Tibet plateau along the profile Maqin-Lanzhou-Jingbian.Chinese J.Geophys.(in Chinese),48(5):1205-1216.

    Unsworth M J,Jones A G,Wei W B,etal.2005.Crustal rheology of the Himalaya and Southern Tibet inferred from magnetotelluric data.Nature,438(7064):78-81,doi:10.1038/nature04154.

    Unsworth M J,Wei W B,Jones A G,etal.2004.Crustal and upper mantle structure of Northern Tibet imaged with magnetotelluric data.J.Geophys.Res.,109(B2),doi:10.1029/2002JB002305.

    Wei W B,Jin S,Ye G F,etal.2009.Conductivity structure and rheological property of lithosphere in Southern Tibet inferred from super-broadband magnetotulleric sounding.Sci China Ser D-Earth Sci(in Chinese),39(11):1591-1606,doi:10.1007/s11430-010-0001-7.

    Wei W B,Unsworth M J,Jones A G,etal.2001.Detection of widespread fluids in the Tibetan crust by magnetotelluric studies.Science,292(5517):716-719,doi:10.1126/science.1010580.

    Wu H W,Zhang L S,Ji S C.1989.The Red River-Ailaoshan fault zone——a Himalayan large sinistral strike-slip intracontinental shear zone.Scientia Geologica Sinica(in Chinese),1:1-10.

    Ye G F,Jin S,Wei W B,etal.2007.Research of conductive structure of crust and upper mantle beneath the south-central Tibetan Plateau.Earth Science-Journal of China University of Geosciences(in Chinese),32(4):491-498.

    Zhan Y,Zhao G Z,Wang J J,etal.2005.Crustal electric structure of Haiyuan arcuate tectonic region in the northeastern margin of Qinghai-Xizang Plateau,China.Acta Seismologica Sinica(in Chinese),27(4):431-440.

    Zhang Z J,Bai Z M,Wang C Y,ec al.2005.Crustal structure of Gondwana-and Yangtze-Typed blocks:An example by wideangle seismic profile from Menglian to Malong in Western Yunnna.Sci China Ser D-Earth Sci(in Chinese),34(5):387-392.

    Zhao G Z,Chen X B,Wang L F,etal.2008.MT evidence for the“pipe flow”in the crust of the northeastern margin of Qinghai-Xizang Plateau.Chinese Science Bulletin(in Chinese),53(3):345-350.

    Zhao G Z,Tang J,zhan Y,etal.2004.Electrical conductivity structure of the crust in the northeastern margin of the Qinghai-Tibet plateau and relationship of the block deformation.Sci China Ser D-Earth Sci(in Chinese),34(10):908-918.

    Zheng Q A,Yu G F,Wang W X.2006.The new knowledge of geological(fault)tectonic framework of Yunnan.Yunnan Geology(in Chinese),25(2):119-124.,doi:CN53-1041/P ISSN1004-1885.

    Zhong K H,Liu Z C,Shu L S,etal.2004.The Cenozoic strike-slip kinematics of the Lancangjiang fault zone.Geological Review(in Chinese),50(1):1-8.

    Zhou H T,Guo L Q.2011.Study on the current activity of main faults in Yunnan.Earthquake Research in Shanxi(in Chinese),146:29-32.

    附中文參考文獻

    白登海,騰吉文,馬曉冰等.2011.大地電磁觀測揭示青藏高原東部存在兩條地殼物質(zhì)流.中國基礎科學,13(1):7-10,doi:10.3969/j.issn.1009-2412.

    蔡軍濤,陳小斌,趙國澤.2010.大地電磁資料精細處理和二維反演解釋技術研究(一)——阻抗張量分解與構(gòu)造維性分析.地球物理學報,53(10):2516-2526,doi:10.3969/j.issn.0001-5733.

    柴天俊,楊繼武.1992.紅河斷裂帶新構(gòu)造活動的地質(zhì)地貌特征.云南地質(zhì),1992,11(3):260-267,224.

    董云鵬,朱炳泉.1999.滇東南建水島弧型枕狀熔巖及其對華南古特提斯的制約.科學通報,44(21):2323-2328.

    皇甫崗,陳颙,秦嘉政等.2010.云南地震活動性.昆明:云南科技出版社.

    皇甫崗,石紹先,蘇有錦.2000.20世紀云南地震活動研究.地震研究,23(1):1-9.

    金勝,張樂天,魏文博等.2010.中國大陸深探測的大地電磁測深研究.地質(zhì)學報,84(6):808-817.

    李祥根.2010.中國地震構(gòu)造運動.北京:地震出版社.

    劉美,白登海,肖鵬飛.2010.青藏高原東部巖石圈電性結(jié)構(gòu)特征及其構(gòu)造意義.地震地質(zhì),32(1):51-58,doi:10.3969/j.issn.0253-4967.

    錢曉東,秦嘉政.2008.小江斷裂帶及周邊地區(qū)強震危險性分析.地震研究,31(4):354-361.

    蘇有錦,劉祖蔭,蔡民軍等.1999.云南地區(qū)強震分布的深部地球介質(zhì)背景.地震學報,21(3):313-322.

    孫潔,晉光文,白登海等.2003.青藏高原東緣地殼、上地幔電性結(jié)構(gòu)探測及其構(gòu)造意義.中國科學D輯,33(增刊):173-180.

    孫潔,徐常芳,江釗等.1989.滇西地區(qū)地殼上地幔電性結(jié)構(gòu)與地殼構(gòu)造活動的關系.地震地質(zhì),11(1):1-11.

    湯吉,詹艷,趙國澤等.2005.青藏高原東北緣瑪沁—蘭州—靖邊剖面地殼上地幔電性結(jié)構(gòu)研究.地球物理學報,48(5):1205-1216.

    魏文博,金勝,葉高峰等.2009.藏南巖石圈導電性結(jié)構(gòu)與流變性——超寬頻帶大地電磁測深研究結(jié)果.中國科學D輯.39(11):1591-1606,doi:10.1007/s11430-010-0001-7.

    吳海威,張連生,嵇少丞.1989.紅河—哀牢山斷裂帶——喜山期陸內(nèi)大型左行走滑剪切帶.地質(zhì)科學,(1):1-8.

    葉高峰,金勝,魏文博等.2007.西藏高原中南部地殼與上地幔導電性結(jié)構(gòu).地球科學:中國地質(zhì)大學學報,32(4):491-498.

    云南省地質(zhì)礦產(chǎn)局.1990.云南省區(qū)域地質(zhì)志.北京:地質(zhì)出版社.

    詹艷,趙國澤,王繼軍等.2005.青藏高原東北緣海原弧形構(gòu)造區(qū)地殼電性結(jié)構(gòu)探測研究.地震學報,27(4):431-440.

    張中杰,白志明,王椿鏞等.2005.岡瓦納型和揚子型地塊地殼結(jié)構(gòu):以滇西孟連—馬龍寬角反射剖面為例.中國科學D輯,35(5):387-392.

    趙國澤,陳小斌,王立鳳等.2008.青藏高原東邊緣地殼“管流”層的電磁探測證據(jù).科學通報,53(3):345-350.

    趙國澤,湯吉,詹艷等.2004.青藏高原東北緣地殼電性結(jié)構(gòu)和地塊變形關系的研究.中國科學D輯:地球科學.34(10):908-918.

    鄭慶鰲,俞國芬,王維賢.2006.云南地質(zhì)(斷裂)構(gòu)造骨架新認識.云南地質(zhì),25(2):119-124,doi:CN53-1041/P ISSN1004-1885.

    鐘康惠,劉肇昌,舒良樹等.2004.瀾滄江斷裂帶的新生代走滑運動學特點.地質(zhì)評論.50(1):1-8.

    周海濤,郭良遷.2011.云南地區(qū)主要斷裂現(xiàn)今活動性研究.山西地震,(2):28-32.

    猜你喜歡
    震區(qū)電性斷裂帶
    流浪衛(wèi)星
    冷凍斷裂帶儲層預測研究
    依蘭—伊通斷裂帶黑龍江段構(gòu)造運動特征
    地震研究(2021年1期)2021-04-13 01:05:12
    民間引爆網(wǎng)絡事件的輿情特點——以“北電性侵事件”為例
    新聞傳播(2018年21期)2019-01-31 02:42:00
    用于燃燒正電性金屬的合金的方法
    色譜相關系數(shù)和隨鉆電性參數(shù)實時評價地層流體方法
    錄井工程(2017年3期)2018-01-22 08:40:08
    蘆山震區(qū)大田壩崩塌發(fā)育特征及其防治措施
    準噶爾盆地西北緣克-夏斷裂帶構(gòu)造特征新認識
    郯廬斷裂帶及兩側(cè)地區(qū)強震異常特征分析
    中國地震(2015年1期)2015-11-08 11:11:26
    帶電粒子在磁場中的多解問題
    久久精品亚洲精品国产色婷小说| 久久精品国产a三级三级三级| 欧美日韩亚洲国产一区二区在线观看 | 国产欧美日韩综合在线一区二区| 日韩三级视频一区二区三区| 国产xxxxx性猛交| 超碰成人久久| 久久久久精品国产欧美久久久| 在线观看免费日韩欧美大片| 黄色视频不卡| 亚洲第一av免费看| 欧美精品高潮呻吟av久久| 国产亚洲av高清不卡| 亚洲午夜理论影院| 午夜成年电影在线免费观看| 精品熟女少妇八av免费久了| 天天操日日干夜夜撸| 国产男女超爽视频在线观看| 久久久水蜜桃国产精品网| 身体一侧抽搐| 老司机深夜福利视频在线观看| 少妇猛男粗大的猛烈进出视频| 欧美乱色亚洲激情| 看黄色毛片网站| 脱女人内裤的视频| 欧美成人免费av一区二区三区 | 女性被躁到高潮视频| 变态另类成人亚洲欧美熟女 | 青草久久国产| 日韩有码中文字幕| 精品熟女少妇八av免费久了| 国产精华一区二区三区| 一边摸一边抽搐一进一小说 | 亚洲成国产人片在线观看| 母亲3免费完整高清在线观看| 国产蜜桃级精品一区二区三区 | 最近最新免费中文字幕在线| 国产亚洲欧美在线一区二区| 高清视频免费观看一区二区| 亚洲精品美女久久久久99蜜臀| 国产午夜精品久久久久久| 久久精品国产清高在天天线| 捣出白浆h1v1| 国产精品二区激情视频| 午夜福利在线免费观看网站| 亚洲欧美一区二区三区黑人| 超色免费av| 制服诱惑二区| 制服诱惑二区| 国产区一区二久久| 大码成人一级视频| 老司机午夜福利在线观看视频| 三级毛片av免费| 日韩有码中文字幕| 中国美女看黄片| 欧美激情高清一区二区三区| 久久久国产精品麻豆| 两个人免费观看高清视频| 亚洲五月色婷婷综合| 国产高清国产精品国产三级| 国产精品自产拍在线观看55亚洲 | 丝袜美腿诱惑在线| 国产精品影院久久| 亚洲欧美激情在线| 亚洲 欧美一区二区三区| 最新在线观看一区二区三区| 色94色欧美一区二区| 国产91精品成人一区二区三区| 国产亚洲欧美在线一区二区| 久久人人97超碰香蕉20202| 亚洲一区二区三区不卡视频| 老司机福利观看| 波多野结衣一区麻豆| 国产一区二区激情短视频| 亚洲欧美色中文字幕在线| 欧美 亚洲 国产 日韩一| 亚洲成国产人片在线观看| 国产一区二区三区综合在线观看| 日本黄色视频三级网站网址 | 动漫黄色视频在线观看| 又黄又粗又硬又大视频| 亚洲 欧美一区二区三区| 国产区一区二久久| 欧美日本中文国产一区发布| 日韩中文字幕欧美一区二区| 日韩欧美在线二视频 | 国产亚洲欧美精品永久| 日韩三级视频一区二区三区| 丰满饥渴人妻一区二区三| 91大片在线观看| 欧美精品一区二区免费开放| 99在线人妻在线中文字幕 | av网站免费在线观看视频| 在线天堂中文资源库| 男女午夜视频在线观看| 精品一区二区三区视频在线观看免费 | 我的亚洲天堂| 极品人妻少妇av视频| 90打野战视频偷拍视频| av天堂在线播放| 9191精品国产免费久久| 久久亚洲真实| 国产一卡二卡三卡精品| 国产成人欧美| av不卡在线播放| 精品第一国产精品| 亚洲免费av在线视频| 母亲3免费完整高清在线观看| 99re在线观看精品视频| 欧美不卡视频在线免费观看 | 无限看片的www在线观看| 少妇裸体淫交视频免费看高清 | 黑人巨大精品欧美一区二区mp4| 最新的欧美精品一区二区| 久久久久久久午夜电影 | 嫩草影视91久久| 亚洲专区字幕在线| 久久ye,这里只有精品| 亚洲成国产人片在线观看| 制服诱惑二区| 精品福利永久在线观看| 国产一区二区三区综合在线观看| 亚洲精品久久成人aⅴ小说| 亚洲少妇的诱惑av| 热re99久久精品国产66热6| 亚洲一区中文字幕在线| 看免费av毛片| 色尼玛亚洲综合影院| 少妇的丰满在线观看| 国产成人av教育| 久久精品成人免费网站| 国产av精品麻豆| 午夜免费观看网址| 在线观看一区二区三区激情| 12—13女人毛片做爰片一| 国产亚洲av高清不卡| 久久久久国产精品人妻aⅴ院 | 18禁国产床啪视频网站| 日韩制服丝袜自拍偷拍| 午夜激情av网站| 99国产精品一区二区三区| 黄色毛片三级朝国网站| 黄色丝袜av网址大全| 777米奇影视久久| 女人精品久久久久毛片| 18禁黄网站禁片午夜丰满| 一区二区三区精品91| 日韩人妻精品一区2区三区| 日本a在线网址| 俄罗斯特黄特色一大片| 日本黄色日本黄色录像| 国产91精品成人一区二区三区| 亚洲国产看品久久| 国产日韩欧美亚洲二区| 在线视频色国产色| 韩国精品一区二区三区| 好男人电影高清在线观看| 一区二区三区激情视频| 变态另类成人亚洲欧美熟女 | 女人久久www免费人成看片| 国产欧美日韩一区二区三| 精品视频人人做人人爽| 久久国产精品男人的天堂亚洲| 欧美精品高潮呻吟av久久| 91成人精品电影| 婷婷精品国产亚洲av在线 | 香蕉国产在线看| 在线观看免费日韩欧美大片| 久久精品成人免费网站| 这个男人来自地球电影免费观看| 天天操日日干夜夜撸| 母亲3免费完整高清在线观看| 黄片大片在线免费观看| 日韩人妻精品一区2区三区| 母亲3免费完整高清在线观看| 人成视频在线观看免费观看| 亚洲欧美一区二区三区黑人| 香蕉丝袜av| 国产亚洲精品久久久久5区| 国产一区二区三区在线臀色熟女 | www.999成人在线观看| 亚洲熟妇熟女久久| 久久久水蜜桃国产精品网| 亚洲久久久国产精品| 欧美黄色片欧美黄色片| 亚洲国产看品久久| 日本黄色视频三级网站网址 | 亚洲黑人精品在线| 久久久久久久精品吃奶| 交换朋友夫妻互换小说| 电影成人av| 久久性视频一级片| 99精品久久久久人妻精品| 久久香蕉激情| 成人黄色视频免费在线看| 一区在线观看完整版| 日韩三级视频一区二区三区| 日韩欧美三级三区| 一进一出抽搐gif免费好疼 | 欧美日韩黄片免| 丝袜在线中文字幕| 不卡一级毛片| 动漫黄色视频在线观看| 老汉色∧v一级毛片| 亚洲全国av大片| 亚洲精品在线观看二区| 婷婷精品国产亚洲av在线 | 人人妻人人澡人人看| 日本撒尿小便嘘嘘汇集6| e午夜精品久久久久久久| 久久久久视频综合| 久久久久视频综合| 高清av免费在线| 色尼玛亚洲综合影院| 黄色a级毛片大全视频| 成人精品一区二区免费| 亚洲国产欧美网| 天天躁日日躁夜夜躁夜夜| 欧美人与性动交α欧美软件| 免费日韩欧美在线观看| 欧美日韩精品网址| 国产亚洲精品久久久久久毛片 | 精品福利永久在线观看| 亚洲一区二区三区欧美精品| 成人手机av| 国产精品美女特级片免费视频播放器 | 国产精品综合久久久久久久免费 | 99re6热这里在线精品视频| 五月开心婷婷网| 男女之事视频高清在线观看| 国产成人精品久久二区二区91| 国产成人影院久久av| 一区在线观看完整版| 国产亚洲欧美98| 丰满人妻熟妇乱又伦精品不卡| 50天的宝宝边吃奶边哭怎么回事| 色播在线永久视频| 亚洲午夜精品一区,二区,三区| 久久久精品区二区三区| 久久国产精品影院| 国产99白浆流出| 新久久久久国产一级毛片| 人妻 亚洲 视频| 男女床上黄色一级片免费看| 日韩视频一区二区在线观看| 久久99一区二区三区| 18禁裸乳无遮挡免费网站照片 | 精品人妻熟女毛片av久久网站| 欧美日韩中文字幕国产精品一区二区三区 | 男女床上黄色一级片免费看| 国产一区二区三区在线臀色熟女 | 国产伦人伦偷精品视频| 男人的好看免费观看在线视频 | 村上凉子中文字幕在线| 久久亚洲精品不卡| 欧美激情 高清一区二区三区| 国产一区二区激情短视频| 天堂中文最新版在线下载| 这个男人来自地球电影免费观看| 人人妻人人澡人人看| 国产成人啪精品午夜网站| 精品国产乱子伦一区二区三区| 亚洲欧美一区二区三区久久| 国产精品免费视频内射| 国产又色又爽无遮挡免费看| 天堂动漫精品| 欧美激情高清一区二区三区| 一区二区三区精品91| 国产97色在线日韩免费| 国产aⅴ精品一区二区三区波| 男男h啪啪无遮挡| 曰老女人黄片| 黄色 视频免费看| 国产又爽黄色视频| 91精品三级在线观看| 女人高潮潮喷娇喘18禁视频| 久久精品熟女亚洲av麻豆精品| 欧美黑人欧美精品刺激| 国产有黄有色有爽视频| 亚洲五月婷婷丁香| 少妇的丰满在线观看| 亚洲欧美日韩高清在线视频| 欧美精品啪啪一区二区三区| 不卡一级毛片| 一a级毛片在线观看| 国产亚洲欧美在线一区二区| 夜夜爽天天搞| 精品欧美一区二区三区在线| 亚洲色图综合在线观看| 国产精品乱码一区二三区的特点 | 91国产中文字幕| 国产亚洲欧美98| 人妻丰满熟妇av一区二区三区 | 国产深夜福利视频在线观看| 欧美黑人精品巨大| 90打野战视频偷拍视频| 男人的好看免费观看在线视频 | 国产不卡一卡二| 国产三级黄色录像| 视频在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 久久久久精品国产欧美久久久| 一本综合久久免费| 每晚都被弄得嗷嗷叫到高潮| 一a级毛片在线观看| 91成年电影在线观看| 超色免费av| 夫妻午夜视频| 色婷婷久久久亚洲欧美| av中文乱码字幕在线| 亚洲中文av在线| 成人黄色视频免费在线看| 91大片在线观看| 久久国产精品影院| 人妻 亚洲 视频| 亚洲黑人精品在线| 一级片'在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 男女免费视频国产| 国产精品秋霞免费鲁丝片| 啦啦啦视频在线资源免费观看| 建设人人有责人人尽责人人享有的| 一区福利在线观看| 男女下面插进去视频免费观看| 中文字幕色久视频| 9191精品国产免费久久| 纯流量卡能插随身wifi吗| 亚洲男人天堂网一区| 91国产中文字幕| 9色porny在线观看| 黑人巨大精品欧美一区二区mp4| av欧美777| 中文字幕人妻丝袜制服| 亚洲成a人片在线一区二区| 看黄色毛片网站| 在线观看www视频免费| 欧美老熟妇乱子伦牲交| 1024视频免费在线观看| 欧美日韩亚洲国产一区二区在线观看 | 一区二区三区国产精品乱码| 纯流量卡能插随身wifi吗| 国产欧美亚洲国产| tocl精华| 50天的宝宝边吃奶边哭怎么回事| 久久久久久免费高清国产稀缺| 丰满的人妻完整版| 国产精品 国内视频| 成熟少妇高潮喷水视频| 老司机深夜福利视频在线观看| av中文乱码字幕在线| 久久人妻福利社区极品人妻图片| 亚洲成国产人片在线观看| 免费在线观看视频国产中文字幕亚洲| 日本五十路高清| 成人黄色视频免费在线看| 一边摸一边做爽爽视频免费| 国产精品.久久久| 国产片内射在线| 欧美日韩亚洲高清精品| 嫁个100分男人电影在线观看| 老熟妇乱子伦视频在线观看| 99国产综合亚洲精品| 黄色怎么调成土黄色| 亚洲国产看品久久| 久久这里只有精品19| 老汉色∧v一级毛片| 这个男人来自地球电影免费观看| 欧美性长视频在线观看| 国产一卡二卡三卡精品| av网站在线播放免费| 嫩草影视91久久| 久9热在线精品视频| 一级a爱视频在线免费观看| 久久人妻福利社区极品人妻图片| 多毛熟女@视频| 天天躁夜夜躁狠狠躁躁| 国产黄色免费在线视频| 咕卡用的链子| 天堂√8在线中文| 精品午夜福利视频在线观看一区| 黄色 视频免费看| 一夜夜www| 多毛熟女@视频| 久久国产乱子伦精品免费另类| 亚洲色图 男人天堂 中文字幕| 曰老女人黄片| 国产成人精品在线电影| 欧美另类亚洲清纯唯美| 久久久国产欧美日韩av| 国产精品秋霞免费鲁丝片| 一区二区三区激情视频| 日韩视频一区二区在线观看| 男女高潮啪啪啪动态图| 中文字幕制服av| av网站免费在线观看视频| 日日爽夜夜爽网站| 在线观看免费视频网站a站| 在线av久久热| 一级黄色大片毛片| 国产精品.久久久| 一级a爱视频在线免费观看| 日韩欧美一区视频在线观看| 99久久精品国产亚洲精品| 成人特级黄色片久久久久久久| 亚洲精品粉嫩美女一区| 国产精品免费大片| 一区二区三区激情视频| 免费高清在线观看日韩| 国产亚洲欧美精品永久| 女人被躁到高潮嗷嗷叫费观| 色94色欧美一区二区| 丰满迷人的少妇在线观看| 欧美老熟妇乱子伦牲交| 国产高清国产精品国产三级| 日韩欧美国产一区二区入口| 后天国语完整版免费观看| 又黄又粗又硬又大视频| 99在线人妻在线中文字幕 | 国产精品一区二区精品视频观看| 露出奶头的视频| 成人国语在线视频| 成人特级黄色片久久久久久久| 一夜夜www| 看黄色毛片网站| 亚洲,欧美精品.| 少妇的丰满在线观看| 老熟妇仑乱视频hdxx| 国产免费现黄频在线看| xxxhd国产人妻xxx| 精品亚洲成a人片在线观看| 欧美在线一区亚洲| 看片在线看免费视频| 午夜福利欧美成人| 在线观看日韩欧美| 国产精品98久久久久久宅男小说| 精品久久久久久久毛片微露脸| 亚洲精品久久午夜乱码| 久久国产乱子伦精品免费另类| 最新的欧美精品一区二区| 好看av亚洲va欧美ⅴa在| 人妻一区二区av| 国产精品亚洲av一区麻豆| 超色免费av| 亚洲中文日韩欧美视频| 欧美黑人精品巨大| 黄色片一级片一级黄色片| 欧美人与性动交α欧美软件| 国产午夜精品久久久久久| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| 日韩三级视频一区二区三区| 亚洲国产中文字幕在线视频| 涩涩av久久男人的天堂| 国产成人影院久久av| 老司机亚洲免费影院| 黄片小视频在线播放| 精品无人区乱码1区二区| 黄色成人免费大全| 欧美日韩亚洲综合一区二区三区_| 亚洲欧洲精品一区二区精品久久久| 91字幕亚洲| 一区二区三区精品91| 久久精品亚洲熟妇少妇任你| av有码第一页| 黄色a级毛片大全视频| 日本欧美视频一区| 久久香蕉激情| 一边摸一边做爽爽视频免费| 国产精品自产拍在线观看55亚洲 | 午夜免费观看网址| 欧美久久黑人一区二区| 欧美激情 高清一区二区三区| 午夜久久久在线观看| 纯流量卡能插随身wifi吗| 国产精品 国内视频| 王馨瑶露胸无遮挡在线观看| a级片在线免费高清观看视频| 国产区一区二久久| 男人舔女人的私密视频| 在线观看午夜福利视频| 国产精品国产高清国产av | 成人永久免费在线观看视频| 18禁美女被吸乳视频| 人人妻人人添人人爽欧美一区卜| 亚洲成人免费av在线播放| 一边摸一边抽搐一进一小说 | 国产精品一区二区在线观看99| 中文字幕最新亚洲高清| 亚洲国产精品sss在线观看 | 老汉色av国产亚洲站长工具| 午夜成年电影在线免费观看| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 欧美精品亚洲一区二区| 国产有黄有色有爽视频| 亚洲 欧美一区二区三区| 另类亚洲欧美激情| 成年女人毛片免费观看观看9 | 亚洲精品美女久久av网站| 欧洲精品卡2卡3卡4卡5卡区| 女人高潮潮喷娇喘18禁视频| 波多野结衣一区麻豆| 一本大道久久a久久精品| 国产一区二区激情短视频| 丰满迷人的少妇在线观看| 一a级毛片在线观看| 国产区一区二久久| 久久久久久久国产电影| 黄片大片在线免费观看| a级毛片黄视频| 久久人妻福利社区极品人妻图片| 国产激情欧美一区二区| 欧美老熟妇乱子伦牲交| 高清av免费在线| 亚洲少妇的诱惑av| 欧美日韩一级在线毛片| 中文亚洲av片在线观看爽 | 日韩免费高清中文字幕av| 在线天堂中文资源库| 国产aⅴ精品一区二区三区波| 精品一区二区三卡| 亚洲成a人片在线一区二区| 成人av一区二区三区在线看| 18禁裸乳无遮挡免费网站照片 | 我的亚洲天堂| 久久亚洲精品不卡| 婷婷丁香在线五月| 国产野战对白在线观看| 国产成人一区二区三区免费视频网站| 身体一侧抽搐| 国产精华一区二区三区| 亚洲国产精品合色在线| 香蕉丝袜av| 91在线观看av| 女人久久www免费人成看片| 精品国产乱子伦一区二区三区| 1024香蕉在线观看| 欧美国产精品一级二级三级| 啦啦啦免费观看视频1| 欧美黄色淫秽网站| 午夜免费成人在线视频| 欧美日韩一级在线毛片| 不卡一级毛片| 中国美女看黄片| 老司机深夜福利视频在线观看| 成在线人永久免费视频| 免费观看精品视频网站| 日韩欧美一区视频在线观看| 一级片'在线观看视频| 在线观看一区二区三区激情| 免费不卡黄色视频| 午夜免费成人在线视频| 久久这里只有精品19| 嫩草影视91久久| 99久久国产精品久久久| 亚洲全国av大片| 国产国语露脸激情在线看| 一进一出抽搐gif免费好疼 | 久久午夜亚洲精品久久| 精品一区二区三区视频在线观看免费 | 成年版毛片免费区| 69av精品久久久久久| 色老头精品视频在线观看| 亚洲av日韩在线播放| 久久久久国产精品人妻aⅴ院 | 欧美乱色亚洲激情| 少妇裸体淫交视频免费看高清 | 男女床上黄色一级片免费看| 久久人妻av系列| 黄片小视频在线播放| 久久人妻av系列| 亚洲av美国av| 国产成人av激情在线播放| 天堂动漫精品| 一级毛片女人18水好多| 日本黄色日本黄色录像| 脱女人内裤的视频| 日本黄色日本黄色录像| 国产激情欧美一区二区| 美女高潮到喷水免费观看| 国产成人免费无遮挡视频| 五月开心婷婷网| 热re99久久国产66热| 日韩视频一区二区在线观看| 精品久久久精品久久久| 中文字幕另类日韩欧美亚洲嫩草| 咕卡用的链子| 亚洲av第一区精品v没综合| 又黄又爽又免费观看的视频| 夜夜爽天天搞| 国产亚洲精品久久久久久毛片 | 变态另类成人亚洲欧美熟女 | 999久久久精品免费观看国产| 午夜精品在线福利| 中文字幕制服av| 精品久久久久久电影网| 老熟女久久久| 色婷婷av一区二区三区视频| 天天躁狠狠躁夜夜躁狠狠躁| 大片电影免费在线观看免费| 两个人免费观看高清视频| 人妻久久中文字幕网| 一区二区日韩欧美中文字幕| 久久久久精品国产欧美久久久| 精品乱码久久久久久99久播| 精品国产乱子伦一区二区三区| 国产三级黄色录像| 韩国av一区二区三区四区| 亚洲欧美一区二区三区黑人| 老司机在亚洲福利影院| 99riav亚洲国产免费| 999精品在线视频| 美女 人体艺术 gogo| 国产在线观看jvid| 国产精品免费大片|