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

    孕震構(gòu)造勢位(LDLEPS)及其相關(guān)的地震地質(zhì)信息分析

    2012-12-07 07:13:26黃秀銘
    地震地質(zhì) 2012年3期

    黃秀銘

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

    孕震構(gòu)造勢位(LDLEPS)及其相關(guān)的地震地質(zhì)信息分析

    黃秀銘

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

    在全變率(NDR)指標的基礎(chǔ)上,結(jié)合GIS技術(shù)進一步開發(fā)了孕震構(gòu)造勢位(LDLEPS)分析法,以便定量評估地震危險性。在充分利用已有資料基礎(chǔ)上繪制了華北地區(qū)LDLEPS圖和地震分布圖,從圖中發(fā)現(xiàn),孕震構(gòu)造勢位與地震空間分布的輪廓、位置、震級大小等是完全匹配的,即LDLEPS圖可以用來預測未來震情。為此,在2006年繪制了四川、青海等地LDLEPS圖,并據(jù)圖對該地未來震情進行了試驗性預測,后來的結(jié)果表明預測基本上是正確的。文中進一步討論了孕震構(gòu)造的形成機制:華北幔柱亞熱柱活動后期對地殼的熱作用與構(gòu)造再造作用而衍生出孕震的二元結(jié)構(gòu)模型,即下層半塑性流變層與上層脆性層,下層動力因素控制上層孕震構(gòu)造的形成與發(fā)展。根據(jù)LDLEPS圖把華北地區(qū)劃分出甲~己6個孕震構(gòu)造異常區(qū),其中乙異常區(qū)(京西北)地熱資源豐富,是深部熱物質(zhì)溢出地表的表現(xiàn)。其地震成因是:下層流變層的流體沿斷層面涌入上層,斷層復活致震,這是下層動力因素對上層間接作用的結(jié)果,其發(fā)震機制類似于水庫誘發(fā)地震。乙異常區(qū)僅產(chǎn)生中等震級以下的地震。甲孕震構(gòu)造異常區(qū)(渤海及其邊緣帶)是下層流變層強勁動力因素直接控制上層孕震構(gòu)造的形成與發(fā)展,從而產(chǎn)生中頻強震以下的地震。另外對探槽剖面及海相沉積物震積巖的古地震遺跡研究發(fā)現(xiàn),甲區(qū)史前曾發(fā)生大量古地震事件與古海嘯事件。現(xiàn)今該區(qū)高勢位孕震構(gòu)造和強震事件是過去事件的延續(xù)和重復。同時指出,華北南緣之SN地塊在構(gòu)造地形與震情上存在明顯差異。

    GIS 地殼運動 地震危險性 孕震構(gòu)造勢位(LDLEPS) 古地震

    0 引言

    近年來國內(nèi)日新月異的計算機技術(shù)促進了地學研究迅速進展。隨著圖形技術(shù)與人機接口技術(shù)的發(fā)展,人們可以更多地同圖形信息進行交互,可以利用高級的人機接口設(shè)備沉浸到由計算機生成的多種信息源構(gòu)成多維信息空間中探索和分析問題。由虛擬現(xiàn)實系統(tǒng)生成的多維地質(zhì)信息空間模擬地質(zhì)環(huán)境。如3D地表環(huán)境污染場、構(gòu)造運動應(yīng)力場和本文的LDLEPS分布場等。通過GIS手段,可以使抽象的數(shù)據(jù)模型化,用3D圖像語言把各種復雜地質(zhì)現(xiàn)象變化的規(guī)律表現(xiàn)得淋漓盡致,一目了然。

    地震是地殼運動突變事件之一。地殼運動含有豐富的孕震信息,充分了解地殼運動規(guī)律,是研究地震事件必經(jīng)之路。筆者曾指出:速度(V)、地變率(Az)與全變率(NDR)分別表示地殼運動Ⅰ、Ⅱ和Ⅲ型的指標;地震事件與前者沒有直接聯(lián)系(黃秀銘,1988,1994);并通過具體震例證明:后二者具有定量評估地震危險性的功能。近年來,筆者在上述研究基礎(chǔ)上,運用GIS技術(shù),創(chuàng)制了孕震構(gòu)造勢位(LDLEPS)工作法,對地震地質(zhì)研究與震情預測工作,在可視性、可操作性、有效性和定量化方面,作了新的探索與嘗試。

    所謂LDLEPS是指能夠反映處于殼中震源層高應(yīng)變狀態(tài)的孕震體孕育勢態(tài)的地殼運動數(shù)據(jù)模型,用3D圖像語言使之轉(zhuǎn)換成可定量評估地震危險性的虛擬地質(zhì)環(huán)境的動態(tài)指標。

    1 LDLEPS制作

    這里Mt為測區(qū)單元的地殼運動異??倧姸?S為測區(qū)單元地殼變形彎曲度;τ為區(qū)內(nèi)單元地殼變形反差度;K為因構(gòu)造地質(zhì)環(huán)境差異而引入的歸一化的修正系數(shù)。

    LDLEPS工作原理類同NDR(黃秀銘,1996)。依照上述思路,LDLEPS設(shè)制如下:

    (1)NDR制作。在標有地殼運動參量值的地區(qū),做以下計算

    (2)LDLEPS制作。在NDR工作基礎(chǔ)上,利用GIS技術(shù),模擬地下孕震構(gòu)造勢態(tài),把能夠反映震前不同的震情勢態(tài)的NDR數(shù)值場,用其不同的幅度及其分布范圍的峰狀體圖像分別表示不同級別的孕震構(gòu)造體的存在、分布范圍、危險程度及其發(fā)展趨勢,從而作出LDLEPS圖(圖1)。然后把處于同一投影坐標系的地震震中疊加到圖上。其中1966—1998年華北(或1950—1998年青海、四川等地)地震震中資料主要取自鄧志輝和楊竹轉(zhuǎn)提供的“mapsis”應(yīng)用軟件。2000—2010年華北(或青海、四川等地)地震目錄選自網(wǎng)絡(luò)資料。因篇幅有限,LDLEPS作圖過程與操作步驟請參閱本文參考文獻(羅云啟,2003;Duchham et al.,2003;吳信才,2004;Christopher,2008;Alias,2008,Chang,et al.,2009)。

    圖1 華北地區(qū)1966—2010年LDLEPS 3D圖與1966—1998年及2000—2010年地震震中3D分布圖Fig.1 A 3D map of both 1966-2010 LDLEPS and earthquakes during 1966-1998 and 2000-2010 in the North China region.

    2 LDLEPS的應(yīng)用條件

    此方法適合于板內(nèi)震前地殼運動參量變化場,尤其適用于淺源型多震區(qū);地下孕震構(gòu)造變動的高剪切應(yīng)變場可徑向傳遞到地表,可充分表現(xiàn)在地表垂直差異運動上;本法適用于中長期地震預測工作。

    該方法不適于震后短期地殼運動參量變化場。這是因為:大震后由于震源體應(yīng)力充分釋放、彈性回跳等,造成地殼應(yīng)力重新調(diào)整,地表塊體處于松弛狀態(tài),即震中附近很大范圍的地表塊體呈“四分五裂”狀態(tài)(Riccardo,2004)。此時獲得的LDLEPS信息混雜著此類地表無序運動的信息,因而不能反映地下孕震構(gòu)造變動的真實情況。

    3 華北地區(qū)LDLEPS特征及其相關(guān)的地震地質(zhì)信息分析

    圖1包含兩部分內(nèi)容,華北1966—2010年 LDLEPS分布和1966—1998年、2000—2010年震中分布。下面對二者分布的空間關(guān)系,LDLEPS內(nèi)容與該區(qū)地震地質(zhì)信息的關(guān)系,以及LDLEPS功能等方面,進行論述。

    3.1 圖1孕震構(gòu)造及其勢位的特征

    圖1中LDLEPS分布圖顯示:既有高矗峰狀的高勢位孕震構(gòu)造區(qū),也有呈平坦或低矮細小條形峰狀的低位勢位值地區(qū),還有介于二者之間的不同幅度的勢位值居中的地區(qū)。有的孕震構(gòu)造圖像彼此關(guān)系緊密,成叢或成串地組合在一起,有的稀疏、松散,彼此互不相聯(lián)。據(jù)此,將圖1中的LDLEPS劃分成甲~己6區(qū),某些區(qū)再進一步劃分為幾個次區(qū)。甲異常區(qū):峰狀體成叢地組合,面積較大,如邢臺(甲1)、渤海西側(cè)(甲2)、海城(甲3)和唐山地區(qū)(甲4);乙異常區(qū):在一定范圍內(nèi)呈多條成串的峰狀鏈,如從京西北至包頭五原一帶;丙異常區(qū):偶有中等幅度峰狀體及多條低矮的峰林鏈的圖像,如銀川附近;丁異常區(qū):峰狀體圖像成叢有度、峰幅不高,如鄭州以東、東北和魯西南一帶;戍異常區(qū):呈規(guī)模小、幅度不大的孤立峰狀體分布,如石家莊西南、太原東南與鄭州以北交會帶,即昔陽、和順及平定等地(戍a區(qū)),太原附近(戍b區(qū)),蒲縣、霍州附近(戍c區(qū))和運城一帶(戍d區(qū))等;己正常區(qū):孤立的峰狀體、低矮呈小鏈條的小峰狀體與平坦的地區(qū)。

    3.2 圖1中與LDLEPS圖相對應(yīng)的1966—1998年及2000—2010年華北地震活動

    此期華北大震頻頻,中小震不斷,地震較密集地分布在甲至戍的LDLEPS異常區(qū)。甲1-甲4區(qū)出現(xiàn)強震及其余震(編號見圖1):1)1966年M7.2邢臺地震(包括1966年M6.8河間地震、M6.7寧晉東南地震、M6.2束鹿南地震,1967年M6.3河間大城地震、M5.7束鹿南地震);2)1969年M7.4渤海地震;3)1975年M7.3海城地震;4)1976年M7.8唐山地震;5)1976年M7.1灤縣地震;乙區(qū):京西北一帶,如1990年M4.9大海坨地震(京北);6)1976年M6.2和林格爾地震;7)1979年M6五原地震;8)1981年M 5.8豐鎮(zhèn)地震;9)1989年M6.1大同堡村地震;10)1996年M6.4包頭地震;11)1998年M6.2張北地震;丙區(qū):發(fā)生大量中小地震,12)1976年M6.2內(nèi)蒙古阿拉善左面旗北地震、1982年M5.7寧夏海源地震等;丁區(qū):1998年M5.9菏澤地震,近20年來在河南鄭州東側(cè)與東北側(cè)小震頻仍,如2008年M4.3封丘-蘭考地震和M3濮陽地震、2010年M4.7周口市(太康縣、扶林縣和西華縣交界)地震等;戍區(qū):近年來小震頻繁,如戍a區(qū)2010年M3.1昔陽地震、M3.2和順地震及M3.2平定地震;戍b區(qū):1989年ML4.5太原附近地震、1992年ML4.3榆次地震、2010年M4.6陽曲地震;戍c區(qū):2010年M3.1霍州地震、M3.7蒲縣地震和M4.7洪洞地震;戍d區(qū):2010年M3.5和M3.1萬榮地震、M4.8河津地震;己區(qū):幾乎無地震。

    上述分析表明,圖1各異常區(qū)的孕震構(gòu)造圖像類型及其勢位特征與其相應(yīng)的地震分布范圍及規(guī)模是完全匹配對應(yīng)的。即孕震構(gòu)造異常區(qū)的勢位強度甲d>乙d>丙d>丁d>戊d>己d,異常區(qū)地震強度甲e>乙e>丙e>丁e>戊e>己e,相呼對應(yīng)。換言之,圖1的LDLEPS圖本身就具有預測未來震情的功能。例如,2009年發(fā)生在西安附近的M4.4小震和M3陜西榆林府谷神木小震,一度引起人們恐慌,當時許多人為正在籌備世界園藝博覽會的西安安全問題心焦如焚。其實從圖1的LDLEPS判斷,西安地區(qū)充其量只會發(fā)生M≤5的地震。神木地區(qū)(是礦震)亦然,是安全的。這2個地區(qū)以后的震情完全印證了這個判斷的正確性。

    3.3 與圖1中LDLEPS異常區(qū)相對應(yīng)的地震地質(zhì)信息分析

    3.3.1 圖1甲孕震構(gòu)造異常區(qū)的地震地質(zhì)信息

    3.3.1.1 包括甲區(qū)在內(nèi)的渤海及周緣古地震遺跡

    圖1中的甲LDLEPS異常區(qū)涵蓋邢臺、西渤海、海城及唐山地區(qū)(1976年震源體下方為低速體),它們位于渤海斷陷區(qū)邊緣。此乃中國東部新構(gòu)造運動最強烈的地區(qū)。這里地震頻發(fā),古地震遺跡很多,如探槽剖面揭示的唐山古地震遺跡和海相沉積物的震積巖(Seismite)(Seilacher,1969)。震積巖是M>6地震時海底陡坡沉積物滑坡與崩積的產(chǎn)物,被認為是古地震的遺跡。一般標志有:斷裂遞變層、震裂縫、微同沉積斷裂、振動液化卷曲變形構(gòu)造和火焰狀構(gòu)造、液化砂、(泥)巖脈、假結(jié)核、層內(nèi)褶皺等。近年來在渤海地區(qū)發(fā)現(xiàn)大量震積巖。如濟陽拗陷古近紀地層(陳世悅等,2003)和山東惠明凹陷東營古近紀三角洲(袁靜,2004)中的震積巖,山東臨朐紅絲石層中古地震遺跡(田洪水等,2006a);山東安丘地區(qū)郯廬斷裂帶古近紀沖積物震積巖(田洪水等,2006b);惠明凹陷地層的震積現(xiàn)象(張小莉等,2006);東營凹陷現(xiàn)河地區(qū)沙三段震積巖(魏垂高等,2006);濟陽坳陷陽信洼陷古近紀震積巖(王克等,2008)等等;魯中臨朐山旺盆地古地震遺跡:中新世晚期山旺硅藻土頁巖中蘊藏大量動植物化石,其中顯現(xiàn)栩栩如生的動物活動景象與動物臨死掙扎景象并存的大量動物被掩埋的場面,這無疑是突如其來的古地震事件造成的。渤海沿岸還發(fā)現(xiàn)多處古地震引發(fā)的海嘯海侵的古地震遺跡,如兩千多年前萊州灣南岸濱海地帶,漢武帝曾視察過現(xiàn)今昌邑城北瓦城村莊,如今僅剩一片瓦礫,城址與古村址均被海砂埋積了。西漢末期萊州灣海侵,井水變咸,迫使壽光、益都和廣饒縣城南遷,并有地震記錄佐證。綜上可見,渤海地區(qū)過去一直是強震多發(fā)區(qū)。

    3.3.1.2 圖1中甲1 LDLEPS異常區(qū)——邢臺地區(qū)的古地震遺跡

    甲1邢臺LDLEPS異常區(qū)是圖1中甲孕震構(gòu)造異常區(qū)的重要組成部分。這里以甲1邢臺古地震遺跡為例進行說明。邢臺地區(qū)位于河北滏陽河與小漳河之間,即束鹿斷陷盆地。1966年地震高發(fā)區(qū)是沿著滏陽河走向延伸到渤海邊。其震源體下方為低速體,極震區(qū)被若干組深淺不同的斷層切割成破碎帶。邢臺地區(qū)古地震遺跡很多。現(xiàn)今寧晉泊是古代大面積廣阿沼澤的一部分。在這樣地表坑坑洼洼易被洪水泛濫的泥沙淤平的環(huán)境中,出現(xiàn)大型湖泊斷陷地,當屬地震構(gòu)造產(chǎn)物。據(jù)鉆孔資料揭示,在邢臺寧晉地面以下1m左右,廣泛分布著湖相淤泥黑色土層,在地面以下50m深度之內(nèi),它們共有3層:1~13m;18~26m和30~47m深度。這說明過去這里經(jīng)歷過地震構(gòu)造斷陷與河流泥沙淤平的重復地質(zhì)過程。由此推測歷史上這里可能發(fā)生3次古強震事件,表明邢臺地區(qū)過去也是強震多發(fā)的地區(qū)。

    圖1中甲LDLEPS異常區(qū)和強震出現(xiàn)的部位正處于渤海裂谷地震構(gòu)造分布范圍,它們與那里的古地震信息是一脈相承的,不是偶發(fā)事件,是有其構(gòu)造活動歷史淵源的。

    3.3.2 圖1中乙孕震構(gòu)造異常區(qū)地震地質(zhì)信息分析

    乙孕震構(gòu)造異常區(qū)位于燕山山前與陰山山前的地帶,起自京北,經(jīng)涿鹿、大同、包頭直至五原附近,由許多短小斷裂構(gòu)成的NW向斷裂帶是深部構(gòu)造變動較為劇烈的地帶。

    乙孕震構(gòu)造異常區(qū)的地熱資源豐富,如北京地熱異常區(qū),下花園-懷來-張山營一帶有4個地熱異常區(qū):后郝窯-萬家窯、暖泉-沙營、奚家堡-西辛堡和塘子廟-胡家營,赤城塘子溝,涿鹿湯寧寺,宣化煙筒鐵礦區(qū)溫泉,蔚縣暖泉,廣靈作疃溫泉,三馬坊溫泉,渾源湯頭溫泉、靈丘泉與門頭泉,天鎮(zhèn)溫泉,西冊田南水域,大同縣瓜園李汪澗村地熱溫泉,包頭濱河新區(qū)(萬水泉鎮(zhèn))地熱溫泉等。此區(qū)如此豐富的地熱資源是該區(qū)深部熱流物質(zhì)經(jīng)殼中斷裂裂隙向地表外溢的表現(xiàn)。這與本區(qū)廣泛分布玄武巖、火山口及有色金屬礦體的地質(zhì)事實一樣,意味著該區(qū)深部蘊藏著豐富的熱源流變物質(zhì)。

    據(jù)深部地球物理探測資料(王帥軍等,2011),乙區(qū)上、中、下地殼均有低速層,如張北,涿鹿;懷來上、中地殼存在低速層,天鎮(zhèn)下地殼有低速層,并且地殼出現(xiàn)高、低速層相間的結(jié)構(gòu)(俗稱“三明治”)。低速層的存在意味著那里有半塑性流變層。乙區(qū)豐富的地熱資源表明深部熱物質(zhì)的存在與通過地殼某些構(gòu)造斷裂孔隙向地表外溢的事實,這說明區(qū)內(nèi)低速層-流變層可能存在較強的動力因素。該帶地震活動特點是:小震頻頻,M4~5地震常見,M5~6地震較多,震級上限約為M6.4?;诖?,推測該區(qū)地震成因可能是:軟流層流體上涌到上層-脆性層的斷層面,使其孔隙水壓力增加,或斷層面礦物被溶解掉,摩擦力因之減弱,觸發(fā)斷層復活,孕震構(gòu)造得以形成與發(fā)展,當斷裂面上凹凸障礙體破裂之時,即發(fā)地震。鑒于以上事實與本區(qū)地震活動特點,一般認為乙區(qū)地震成因類似于水庫誘發(fā)地震機制,是流變層動力因素間接誘發(fā)上層地震,其特點是:這里頻發(fā)中等震級以下的地震,不易產(chǎn)生大震;僅在局部出現(xiàn)分布范圍不大、幅度中等的孕震構(gòu)造的峰狀體。

    乙異常區(qū)的高勢位孕震構(gòu)造及地震,均分布在包頭、張北、和林格爾、豐鎮(zhèn)、大同等地,它們正處于地殼厚度突變的和M面陡坡帶之上的地質(zhì)環(huán)境部位。力學分析認為:上述地質(zhì)環(huán)境是不穩(wěn)定的。原因是:假定依附于M陡坡上的是巖體A,其重量為G,垂直向下的重力矢量,可分解成兩個分量g1及g2。g1為垂直于M陡坡面且方向向里的分量,它被陡坡反向支撐力抵消;g2為方向向下且平行陡坡的分量,其中小部分被陡坡面的摩擦力消耗掉,但大部分下滑重力分量仍作用于陡坡上。換言之,任何一塊出現(xiàn)在陡坡上,定然受到與陡坡面平行且下滑的重力分量作用。這樣巖體A無論是靜態(tài)還是動態(tài),均處于不穩(wěn)定的失衡狀態(tài),從而使附近空間存在趨于“變平”勢態(tài),地體易產(chǎn)生雪崩式的滑塌現(xiàn)象,并波及其上巖層,引致孕震構(gòu)造形成。顯然它也成為乙區(qū)形成孕震構(gòu)造的深部構(gòu)造背景條件之一(圖6)。

    簡言之,乙孕震構(gòu)造異常區(qū)出現(xiàn)是與其深部構(gòu)造變動的背景密切相關(guān)的。

    3.4 圖1各處LDLEPS峰狀體組合圖案差異的地震地質(zhì)背景

    由圖1可見,華北各地的LDLEPS峰狀體組合圖案顯著不同。北側(cè)峰狀鏈平直,呈NWW向,這是因為它們受北部內(nèi)蒙古地軸與南部華北平原沉降區(qū)的大地構(gòu)造分界影響。中間地區(qū)的LDLEPS峰狀體組合圖案呈海星放射狀、環(huán)形或流線型弧形鏈狀的形式,鮮見平直峰鏈。這些形態(tài)可能顯示流變層頂面的流體運移的痕跡(下述),或者它們對上層沖擊造就LDLEPS樣貌于平面投影的軌跡。華北南側(cè)的峰狀體幅度相近,彼此組合圖案相似,分布較為零亂,但總體仍隱隱約約顯示EW方向的趨勢,與前述樣貌不同。這是由于它的大地構(gòu)造背景使然。南側(cè)以北為華北地塊,以南為華南地塊。前者由太古代基底巖構(gòu)成,后者以元古代基底為主,兼有少量太古代巖石。它們在元古代末同為超大陸聯(lián)合體,約700Ma前發(fā)生分裂。古生代期間,二者各自獨立運動。經(jīng)歷古生代,原位于南半球的華南地塊開始向北飄移,在三疊紀中期(約230Ma前)與華北地塊碰撞,潛入后者南緣之下。西自西安秦嶺以南,南至大別山,東至蘇魯接壤地帶均有其超高壓變質(zhì)巖帶蹤跡。穩(wěn)定的華南古陸地塊,因其形成久遠的碰撞縫合帶緣故,主要特點是:二者碰撞帶在構(gòu)造與地形上形成諸如秦嶺、中條山與魯南的EW向山脈,一改晉西南風陵渡黃河中上游水系流向多變的習性,至此以下(入黃海的)黃河呈西東流向,那里呈蜈蚣狀的地表主干水系走向EW,NS成行、如渭河、丹江口以上的漢水、入黃海時黃河(風陵渡以東)下游河段、淮河、長江等,其地震活動頗具特色:地震頻度中等,一般僅發(fā)生M6.4以下地震。震中常呈孤狀分布,1900年前主要分布在長江沿岸、南黃海至鄭州環(huán)圍地帶等,1900年之后分布在武漢北側(cè)、九江、溧陽至南黃海環(huán)圍地帶等。這與華北地區(qū)下述特點大相徑庭:即盆嶺構(gòu)造、水系順勢而行,方向多樣與頻發(fā)中、強地震的特點。

    應(yīng)當指出,圖1時段在華南地塊潛入華北地塊南緣的復合帶地震活動并不強烈,但潛伏帶附近獨具特色:具有該地域近EW向強震帶的構(gòu)造背景。歷史上那里強震頻發(fā):如甘南禮縣東北1654年M8地震,武都西南1879年M8地震,天水東734年M7地震及甘谷附近143年M7地震,陜南華縣1556年M8地震及1501年M7朝邑地震,晉南1303年M8洪洞地震及1695年M8臨汾地震,冀南1830年M7.5磁縣西地震,魯西南1937年M7菏澤地震和1668年M8.5魯南湯頭東(郯城—莒縣之間)地震。人們發(fā)現(xiàn):晉南到魯南一帶呈現(xiàn)b值低和LURR峰值的異常帶,即該帶邊界附近可能潛伏地震異常信號(Zhang Lang-ping et al.,2011)。

    上面3個方面的論述表明,LDLEPS的內(nèi)涵彰顯它在洞察一個地區(qū)地震地質(zhì)內(nèi)在矛盾與運行規(guī)律上是富有成效的,同時在定量與可視性方面也能體現(xiàn)其新的思路與風格。更為重要的是,它有助于了解LDLEPS形成的主要原因即下述的二元結(jié)構(gòu)致震模式。

    4 利用LDLEPS方法進行未來震情預測的試驗

    為了更好地了解LDLEPS的地震預測功能,2006年我們又編制了圖2及圖3中的2000—2020年LDLPES 3D圖,2011年又把該區(qū)2000—2010年M8.2震中3D圖疊加到前圖上,則成現(xiàn)在的圖3,對青海、四川等地專門進行未來震情預測的試驗。

    4.1 1950—1998青海、四川等地LDLEPS圖與M4.7~7地震震中圖分析

    4.1.1 圖2孕震構(gòu)造分析

    此圖2可見,孕震構(gòu)造異常區(qū)分布在下列地區(qū):A1區(qū)——從105°E與35°N交會處附近,向NWW方向延伸,經(jīng)蘭州附近,并延及武威與西寧附近的地區(qū);A2區(qū)——從西寧向SW延伸的地區(qū);A3區(qū)——蘭州向NE延及40°N附近的地區(qū);A4區(qū)——100°E東西兩側(cè)與30°N以南的地區(qū);A5區(qū)——昆明東北側(cè)與西南側(cè)的地區(qū)。

    圖2 1950—1998年青海、四川等地LDLEPS 3D圖與M4.7~7地震3D分布圖Fig.2 A 3D map for both LDLEPS and earthquake epicenters with M4.7~7 in Qinghai and Sichuan Provinces in the western part of China during 1950—1998.

    4.1.2 圖2地震分析

    1950—1998年本區(qū)發(fā)生的地震:A1區(qū)8次M4.7~7地震;A2區(qū)3次M6~7地震;A3區(qū)10次M4.7~6.2地震;A4區(qū)7次M6~7地震;A5區(qū)2次M6.3~6.6地震。值得指出,在成都西北側(cè)與西側(cè)的廣大地區(qū)未見孕震構(gòu)造,也未發(fā)生地震。

    4.2 圖3中青海、四川等地2000—2020年LDLEPS圖與2000—2010年M6~8.2地震震中圖分析

    由圖3可見,這里顯示兩處極其明顯的孕震構(gòu)造異常區(qū):1)西北部B區(qū)范圍為90°~100°E與33°~40°N,在圓環(huán)狀基座背景上矗立著許多高勢位孕震構(gòu)造峰狀(叢)體區(qū)域,并在其西北角與南邊分別出現(xiàn)井噴式峰狀體;2)東南部C區(qū)范圍為100°~103°E與30°N交會處,高勢位孕震構(gòu)造的峰狀叢體,向NE方向逐然變窄變矮,朝105°E與35°N相交處附近延伸著,且在成都西側(cè)與西北側(cè)也出現(xiàn)另一井噴式峰狀體。另外,還存在較小的勢位值的D區(qū)范圍:在武威外側(cè)附近、北面40°N附近、40°N與105°E相交附近,以及100°E東西側(cè)與30°N以南的小型孕震構(gòu)造區(qū)。此外還發(fā)現(xiàn):圖3中C區(qū)大型孕震構(gòu)造異常區(qū),特別是成都西側(cè)和西北側(cè)呈井噴式峰狀體恰好位于該區(qū)地殼厚度突變和M面陡傾斜坡(地震危險標志)之上的地質(zhì)環(huán)境部位(圖6)?;趫D3中B和C孕震構(gòu)造異常區(qū)(包括井噴式峰狀體)突然出現(xiàn)的反?,F(xiàn)象與圖6中的深部異常構(gòu)造背景,2006年作者曾預測本區(qū)未來15年震情:B和C區(qū)可能分別發(fā)生M7.3以上與M7.5以上地震。結(jié)果該區(qū)相繼發(fā)生(圖3):B區(qū)2010年M 7.1青海玉樹地震(包括2001年M 8.2青海地震或稱昆侖山山口西MS8.1地震)和C區(qū)2008年M8四川汶川地震(該震源層下方為低速體)及余震。這3處大震震中均分別位于上述大型孕震構(gòu)造異常區(qū)井噴式峰狀體范圍內(nèi)。

    這次做的預測該區(qū)未來震情試驗,在空間位置與時間方面尚可,但震級預測略有偏差。但已經(jīng)證明:LDLEPS確有預測未來震情的功能。

    圖3 青海、四川等地2000—2020年LDLEPS 3D圖與2000—2010年M 6~8.2地震震中3D分布圖Fig.3 A 3D map for both 2000-2020 LDLEPS and earthquake epicenters with M6~8.2 in Qinghai and Sichuan Provinces during 2000—2010.

    5 討論與結(jié)語

    5.1 孕震構(gòu)造形成與發(fā)展的主要原因

    5.1.1 二元結(jié)構(gòu)模型與孕震構(gòu)造成因

    圖4 由下層半塑性流變動力因素引致上層巖體脆性變形的地質(zhì)遺跡Fig.4 Photos showing the geologic relics resulting from the rheologic deformation of the lower rock mass and the brittle deformation of the upper rock mass.

    圖5 剖面上顯示的巖石變形現(xiàn)象Fig.5 Rock deformation revealed on the cross sections.

    幔柱構(gòu)造理論(Manuyama,1994)認為,地幔熱物質(zhì)可直通巖石圈,形成幔柱構(gòu)造。據(jù)李紅陽等(2002)的研究成果、地球物理深部探測資料、地熱異常以及地震地質(zhì)資料,華北地下發(fā)育很多不同級別的幔柱亞熱柱構(gòu)造,如新生代時期,渤海西部深度50~70km出現(xiàn)頂部呈蘑菇狀亞熱柱隆升構(gòu)造,其頂部之上厚度約為20km的低速層分別向周圍山體之下延伸。深部熱物質(zhì)涌入地殼,殼內(nèi)增溫,中下地殼塑性增強。同時隆頂?shù)闹小⑾碌貧は蛲鈧?cè)向流變,頂部中下地殼減薄,向外圍增厚,產(chǎn)生側(cè)向拆離滑脫。地震地質(zhì)觀點強調(diào),幔柱亞熱柱活動后期對地殼的熱作用與構(gòu)造再造作用,衍生出二元結(jié)構(gòu):下層半塑性的流變與上層脆性層,亦即下地殼呈半塑性,中地殼具有高速、低速相間的特征,而上地殼呈現(xiàn)脆性特性。通過火山噴發(fā)、間歇性泥火山噴泉、升至地表的深成巖漿巖與變質(zhì)核雜巖的韌性剪切帶、金屬礦床以及大氣、海洋流體學和地熱異常等資料,人們逐漸了解:殼中流變層物質(zhì)呈4態(tài),包括固體以及流體的3態(tài)——液態(tài)、氣態(tài)和等離子體;存在流變介質(zhì)的對流、擴散以及相伴的物理、化學過程,導致質(zhì)量、動量、能量輸送的動力學現(xiàn)象出現(xiàn)。它包括:1)若輕的流體在下面,重的在上面,就會發(fā)生上下翻滾的對流運動;2)擴散,指2個相鄰的、溫度或濃度不同的流體團間的質(zhì)量、能量、動量交換,湍流擴散可以是分子擴散的幾百、上千倍;3)黏彈性流體,當它發(fā)生旋轉(zhuǎn)時,會以“龍卷風”似的方式沿著中間旋轉(zhuǎn)軸向上爬升;4)流體具有在下層圍壓下向上層巖隙強烈增壓噴瀉作用與動力回抽減壓下巖體坍塌作用;5)對上層產(chǎn)生不均勻浮力的振蕩與搖晃作用;6)流變層頂部大型強勢的水平流動對上層底部產(chǎn)生層間水平切變作用。下層強大動力因素無疑會對上層產(chǎn)生不斷沖擊作用:撕裂上層巖體、流體向上層脆性層的斷裂裂隙滲透與擠壓或令其瓦斯爆炸(類似礦震)等,引發(fā)構(gòu)造變形。圖4為原位于中地殼層位而后隆升至地表的北京云蒙山北部天仙溪黃土良村的地質(zhì)剖面。上下層構(gòu)造形態(tài)迥然不同。圖4c,d,e,f處于中地殼范圍的上層,其巖脈前身乃為不同時期形成的不同方向的裂縫、斷裂及其位移,且后者穿插前者的斷裂裂隙網(wǎng)絡(luò)系。它是歷史積淀的現(xiàn)象,是中地殼范圍近鄰的下層(圖4 a)各期不同方向流變動能對上層控制作用的結(jié)果,之后流體物質(zhì)才上涌并分別充填各自前期的裂隙中,形成巖脈。圖5a,b表示曾處于殼中上層脆性層的兩處相距甚遠的長城系大理巖,共同受到下層(花崗巖質(zhì)流變層)流變動力作用,導致它們同時發(fā)生類似的構(gòu)造變形、碎裂,并部分墜落到下層流變層的花崗巖層中,成為后者的捕虜體。其具體過程可能是:上層脆性層大理巖先發(fā)生構(gòu)造事件-沿葉理面斷裂位移,繼而碎裂,終至個別碎塊沿著撕裂的裂隙下墜到下層流變層中。圖5d表示下層流變層花崗巖質(zhì)流體沿著裂隙侵入上層脆性層閃長巖的地質(zhì)剖面?;鹕降卣鹉P?上部脆性發(fā)震層與下部塑性巖漿層,且二者層位與時同步變淺)凸顯二元結(jié)構(gòu)致震的真諦。地震成因是由下層流變活動引致上層脆性層孕震構(gòu)造形成,致震。其類型包括圖1中甲和乙異常區(qū)的事件以及上述事件;也包括上地殼上層(或中地殼內(nèi)上層)的底部層間滑脫事件,如均為低角度正斷層的北京密云云蒙山東南側(cè)河防口-懷北莊拆離滑脫(片理走向SE)斷層(朱大崗等,2006)和云蒙山石城斷層帶(含有假玄武巖玻璃),二者生成的環(huán)境均是殼中韌性剪切與彈塑性的過渡層即地殼震源層部位(史蘭斌等,1984),以及意大利中部Elba的isle中Zuccale拆離型震源斷層(其性狀:低角度正斷層、斷層角礫巖和碎裂巖、滲透性強、富含CO2含水流體灌入及葉片狀斷層核存有水破裂作用的碳酸鹽巖脈)。

    顯然,下層動力因素使上層應(yīng)變場增強,這是板內(nèi)孕震構(gòu)造形成與發(fā)展的主要原因。

    5.1.2 二元結(jié)構(gòu)模型與孕震構(gòu)造的關(guān)系

    利用地熱異常、玄武巖、有色金屬礦床、剖面低速層以及地震地質(zhì)等資料,編制出由華北幔柱亞熱柱活動后期對地殼熱作用與構(gòu)造再造作用衍生的(部分)二元結(jié)構(gòu)分布圖(圖7)。從古地震遺跡與歷史地震記錄可見,這些二元結(jié)構(gòu)區(qū)活性是變化的。如圖7中4區(qū)自渤海幔柱亞熱柱出現(xiàn)以來一直存在二元結(jié)構(gòu),大震不斷。近300a來其活性集中在其西北、西南、東和北的邊緣。1區(qū)自中衛(wèi)南經(jīng)銀川向NNW至中蒙邊界,并與蒙古和貝加爾裂谷的亞熱柱遙相呼應(yīng)。其活性歷來強勢。16、18、20世紀最活躍,至今余威不減。5區(qū)活性在4ka前北邊強勢,17世紀東南活躍。之后適度轉(zhuǎn)向西南。3區(qū)1ka前LDLEPS值在定襄、忻縣一帶增大。17區(qū)在14、17世紀末活性最強。6、7區(qū)是圖7中最弱的。顯然,這些二元結(jié)構(gòu)的上層的孕震構(gòu)造活化與平靜期的習性,與其垂向相對應(yīng)的下層的相應(yīng)空間部位附近的流變層動力因素活躍與平靜期的屬性相匹配的。分析圖7發(fā)現(xiàn),二元結(jié)構(gòu)活躍的地區(qū),大多出現(xiàn)在以Qh斷裂為界的地塊邊緣帶,如鄂爾多斯塊體邊緣帶,而很少出現(xiàn)在腹部地帶。因為腹部斷裂少,地層水平(傾角不到1°),LDLEPS值低。一般認為,二元結(jié)構(gòu)的下、上層和孕震構(gòu)造呈如下關(guān)系:若下層動力因素弱(如鄂爾多斯地塊下地殼存有含CO2鹽水的流體層,電阻率低,且康氏面附近可能有不透水巖層阻隔(Zhao et al.,2011),流體無力向上層滲透,上層整體性好,則不易產(chǎn)生較大孕震構(gòu)造;若下層流體動力較強,中層整體性差,流體容易上涌,令上層斷裂復活,容易產(chǎn)生孕震構(gòu)造,即圖1中乙異常區(qū)的情形。若下層極強勁,中層整體性一般或較差時,會撕裂上層巖體,易產(chǎn)生類似圖1中甲異常區(qū)的孕震構(gòu)造。

    圖6 圖3中D孕震構(gòu)造異常區(qū)和2008年M 8四川汶川地震及其余震區(qū)之下的M面3D示意圖Fig.6 A 3D map of M surface under both the anomalous LDLEPS region D in Fig.3 and the 2008 M8 Wenchuan epicentral area and its aftershock region.

    圖7 華北地區(qū)幔柱亞熱柱活動后期由地殼熱作用和構(gòu)造再造作用而衍生的(部分)二元結(jié)構(gòu)區(qū)分布Fig.7 The distribution map of bilayer structures generated by both geothermal process and tectonic reworking of the crust during the late-stage evolution of sub-mantle plume and mantle branch in North China.

    5.1.3 二元結(jié)構(gòu)模型有助于詮釋殼中地震地質(zhì)現(xiàn)象

    地殼結(jié)構(gòu)是分層的。殼中主震源層分布在中層位,即花崗巖質(zhì)及其上下界面附近,一般深度在8~22km,并以層狀“三明治”形式交替出現(xiàn)。由于二元結(jié)構(gòu)模型的存在,地殼深層處于半流變狀態(tài),能量便于自由傳遞與轉(zhuǎn)化,應(yīng)力在流動過程中得以快速釋放,卻不利于彈性應(yīng)變能存儲積累和突發(fā)式釋放,故地震較少。地殼表層因距離太遠,下層動力因素對地殼表層的直接作用或影響不甚顯著,一般能夠波及或分配到地殼表層的能量相對是間接緩慢的。此外地殼表層為開放性剛體,隸屬另一構(gòu)造動力學系統(tǒng)。巖石主要表現(xiàn)為脆性,淺部圍壓低,巖石強度低,在應(yīng)力作用下巖體易產(chǎn)生裂隙,應(yīng)力迅速釋放,有點能量可以通過調(diào)節(jié)表層塊體運動的途徑得以緩慢釋放,它的能量存儲積累和突發(fā)性釋放的條件較差,故該層位地震活動亦少。唯有殼中主震源層即二元結(jié)構(gòu)的上層范圍,既易從下層獲得能量,本身又是封閉性剛體,巖體不易生成裂隙,那里又沒有較強的塑性流動,故應(yīng)力便于積累,既具有有利于能量以彈性應(yīng)變能存儲積累又難以快速釋放的特點。那里儲能一旦超過本身負荷就會突發(fā)性地釋放,故該層發(fā)震頻繁。地震時此層地震波向上傳播,沿途幾無耗能,僅于空-地分界面,傳遞受阻,大量能量駐積于地表附近,且于數(shù)百秒內(nèi)全部釋放,并以剩余變形與地面振動形式對建筑物施加破壞。地震問題實際上就是堆積于地表至深度70m左右范圍的能量如何處置的問題。是否可以對這些能量通過反射屏反射、波場干涉、涂料吸收、與建筑物隔絕聯(lián)系等,予以遏抑、轉(zhuǎn)化、削弱,實現(xiàn)減輕震災目標。

    5.2 初步結(jié)論

    綜上所述,可以得出以下初步結(jié)論:

    板內(nèi)幔柱亞熱柱活動后期對地殼的熱作用與構(gòu)造再造作用,衍生出二元結(jié)構(gòu)。LDLEPS是定量評估地震危險性的地質(zhì)指標,適用于板內(nèi)中長期地震預測工作。它兼具新穎性、可視性、可操作性與有效性,工作可以計算機化。它在地學研究、地震預測與工程建設(shè)方面,有著廣闊的應(yīng)用前景。

    工作區(qū)LDLEPS異常區(qū)的出現(xiàn),可能與地殼中二元結(jié)構(gòu)模型密切相關(guān)。下層流變層的動力因素可以對上層脆性層產(chǎn)生控制作用,導致包括二者界面之上底層層間低角度滑脫斷層在內(nèi)的孕震構(gòu)造形成發(fā)展,致震。

    由包括第四紀地殼運動項目在內(nèi)的綜合地震地質(zhì)研究和地球物理深部探測獲取的信息源,是分析二元結(jié)構(gòu)模型的有效工具。

    致謝 在工作過程中,中國科學院地質(zhì)與地球物理研究所趙希濤及高浩中提供云蒙山野外工作條件,原地質(zhì)部系統(tǒng)諸單位同仁野外工作中給予賜教,密云縣旅游局關(guān)昊和優(yōu)山美地公司(北京)馬少輝在野外工作中給予關(guān)照,優(yōu)山美地公司工程師梅友堯無償提供計算應(yīng)用軟件,在此一并表示真誠謝意。

    陳世悅,袁文芳,鄢繼華.2003.濟陽坳陷早第三紀震積巖的發(fā)現(xiàn)及其意義[J].地質(zhì)科學,38(3):377—384.

    CHEN Shi-yue,YUAN Wen-fang,YAN Ji-hua.2003.Discovery and significance of earthquake event deposits of early Tertiary in the Jiyang depression[J].Chinese Journal of Geology,38(3):377—384(in Chinese).

    黃秀銘.1979.關(guān)于古地震地質(zhì)研究[J].地震地質(zhì)譯叢,1(4):34—44.

    HUANG Xiu-ming.1979.Research of palaoeseismic geology[J].Translated Works on Seismogeology,1(4):34—44(in Chinese).

    黃秀銘.1988.地變率指標與華北地震區(qū)劃[J].地震地質(zhì),10(3):33—44.

    HUANG Xiu-ming.1988.Indicators of topographical deformation ratio and seismic zonation of North China [J].Seismology and Geology,10(3)::33—44(in Chinese).

    黃秀銘.1994.地變率指標初探[A].見:劉若新主編.現(xiàn)代地球動力學研究及其應(yīng)用.北京:地震出版社.294—300.

    HUANG Xiu-ming.1994.Preliminary discussion on the indicator system of net deformation ratio(NDR)[A].In:LIU Ruo-xin(ed).Researches and applications of modern earth's dynamics.Seismological Press,Beijing.294—300(in Chinese).

    黃秀銘.1996.地震運動全變率指標(NDR)系與地震事件[J].地震地質(zhì),18(4):427—435.

    HUANG Xiu-ming.1996.The indicator system of net deformation ratio(NDR)of crustal movement and seismic event[J].Seismology and Geology,18(4):427—435(in Chinese).

    黃云波.1998.北京密云典型韌性剪切帶特征[J].礦物巖石地球化學通報,17(3):179—182.

    HUANG Yun-bo.1998.Features of a typical ductile shear zone in Miyun,Beijing[J].Bulletin of Mineralogy Petrology and Geochemistry,17(3):179—182(in Chinese).

    李紅陽,牛樹銀,王立峰,等編著.2002.幔柱構(gòu)造[M].北京:地震出版社.

    LI Hong-yang,NIU Shu-yin,WANG Li-feng,et al.2002.Mantle Plume Structure[M].Seismological Press,Beijing(in Chinese).

    羅云啟.2003.數(shù)字化地理信息系統(tǒng)建設(shè)與Mapinfo高級應(yīng)用[M].北京:清華大學出版社.

    LUO Yun-qi.2003.Digital Geography Information System Construction and Advanced Application of MapInfo[M].Tsinghua University Press,Beijing(in Chinese).

    牛樹銀,孫愛群,邵振國,等.2001.地幔熱柱多級演化及其成礦作用[M].北京:地震出版社.

    NIU Shu-yin,SUN Ai-qun,SHAO Zhen-guo,et al.2001.The Multiple Evolution of Mantle Plume and Its Mineralization[M].Seismological Press,Beijing(in Chinese).

    牛樹銀,史中華,張國軍.2004.華北東部盆嶺區(qū)的形成與演化[J].河北地質(zhì)礦產(chǎn)信息,(2):2—5.

    NIU Shu-yin,SHI Zhong-hua,ZHANG Guo-jun.2004.Formation and evolution of the basin-ridge province in the eastern North China[J].Information of Geology and Mineral Resources in Hebei Province,(2):2—5(in Chinese).

    史蘭斌,陳孝德.1984.北京北石城斷層帶與與河防口斷層帶巖石顯微構(gòu)造研究[J].地震地質(zhì),6(3):1—13.

    SHI Lan-bin,CHEN Xiao-de.1984.Microstructures of fault rocks in Beishicheng and Hefangkou Fault zones,Beijing[J].Seismology and Geology,6(3):1—13(in Chinese).

    田洪水,張增奇,張邦花,等.2006a.山東臨朐紅絲石層中古地震事件記錄[J].中國地質(zhì),33(5):1137—1142.

    TIAN Hong-shui,ZHANG Zeng-qi,ZHANG Bang-hua,et al.2006a.Records of paleoseismic events in the red-silk stone bed in Linqu,Shandong[J].Geology in China,33(5):1137—1142(in Chinese).

    田洪水,張增奇.2006b.山東安丘地區(qū)郯廬斷裂帶古近紀沖積物震積巖序列[J].地質(zhì)科學,41(2):217—225.

    TIAN Hong-shui,ZHANG Zeng-qi.2006b.Seismite succession of the Paleogene alluvium along the Tancheng-Lujiang Fault zone in Anqiu area,Shandong[J].Chinese Journal of Geology,41(2):217—225(in Chinese).

    王克,劉顯陽,趙衛(wèi)衛(wèi),等.2008.濟陽坳陷陽信洼陷古近紀震積巖特征及其地質(zhì)意義[J].巖性油氣藏,20(2):54—59.

    WANG Ke,LIU Xian-yang,ZHAO Wei-wei,et al.2008.Characteristics and geological significance of seismites of Paleogene in Yangxin subsag of Jiyang depression[J].Lithologic Reservoirs,20(2):54—59(in Chinese).

    王帥軍,王夫運,張成科,等.2011.北京及附近地區(qū)地殼速度結(jié)構(gòu)與構(gòu)造[[J].華北地震科學,29(2):6—12.

    WANG Shuai-jun,WANG Fu-yun,ZHANG Cheng-ke,et al.2011.Study on crust velocity structure and tectonic of the capital area[J].North China Earthquake Sciences,29(2):6—12(in Chinese).

    魏垂高,張世奇,姜在興,2006.東營凹陷現(xiàn)河地區(qū)沙三段震積巖特征及其意義[J].沉積學報,24(6):798—805.

    WEI Chui-gao,ZHANG Shi-qi,JIANG Zai-xing.2006.Discovery and significance of seismite of Silurian in member Ⅲat Shahejie formation in Xianhe area of Dongying sag [J].Acta Sedimentologica Sinica,24(6):798—805(in Chinese).

    吳信才.2004.MAPGIS 6.5地理信息系統(tǒng)[M].北京:電子工業(yè)出版社.

    WU Xin-cai.2004.MAPGIS 6.5 Geographic Information System[M].Electronic Industry Press,Beijing(in Chinese).

    袁靜.2004.山東惠明凹陷古近紀震積巖特征及其地質(zhì)意義[J].沉積學報,22(1):41—45.

    YUAN Jing.2004.The property and geological significance of seismites of Paleogene in Huimin sag,Shandong Province[J].Acta Sedimentologica Sinica,22(1):41—45(in Chinese).

    張小莉,查明,楊劍平.2006.惠明凹陷的震積作用及低電阻率油層[J].沉積學報,24(6):829—833.

    ZHANG Xiao-li,ZHA Ming,YANG Jian-ping,et al.2006.Seismic deposition and low contrast resistivity oil layers in Huimin sag[J].Acta Sedimentologica Sinica,24(6):829—833(in Chinese).

    朱大崗,崔盛芹,吳珍漢,等.2006.北京云蒙山地區(qū)擠壓-伸展體系構(gòu)造特征及其巖石組構(gòu)的動力學分析[A].見:燕山陸內(nèi)造山帶造山過程及動力機制.北京:地震出版社.150—158.

    ZHU Da-gang,CUI Sheng-qin,WU Han-zhen,et al.2006.Tectonic features of the compressional-extensional system and dynamic analysis of rock fabric in Yunmengshan area,Beijing[A].In:Orogeny and Dynamic Mechanism of the Yanshan Intracontinental Orogenic Belt.Seismological Press,Beijing.150—158(in Chinese).

    Alias Abdul-Rahman,Morakot Pilouk.2008.Spatial Data Modeling for 3D GIS[M].Berlin:Springer.

    Chang Kang-tsung.2009.Introduction to Geographic Information Systems[M].The McGraw-hill Education(Asia)Companies Ins.

    Christopher P.2008.Landscape Analysis and Visualization:Spatial Models for Natural Resource Management and Planning[M].Berlin:Springer.

    Duckham M,Geodchild M F,Worboys M F,et al.2003.Foundation of Geographic Information Science[M].Taylor&Francis,London.1—18.

    Framk A.1998.Introduction to Geomorphology[M].London:Arnold.

    Hans-Ulrich Schmincke.2010.火山學[M].古今書院.

    Maruyama S P.1994.Tectonics[J].Jour Geol Soc Japan,100(1):24—49.

    Riccardo E M.2004.Crustal Rheology and Postseismic Deformation:Modeling and Application to the Apennines[M].Delft University Press.

    Seilacher A.1969.Fault-graded beds interpreted as seismites[J].Sedimentology,13:155—159.

    Takeshi Kuritani.2007.水とマグマ [J].地學雜誌,116(1):133—153.

    Zhang Lang-ping,Li Zhi-xiong,Shao Zhi-gang,et al.2011.Research on strong earthquake tendency on active tectonic block boundaries in the Chinese mainland[J].Earthquake Research in China,25(2):189—200.

    Zhao Guo-ze.2000.Evidence for the location and cause of large crustal earthquake in Japan [J].J Geophys Res,105,13,579—13,594,20.

    THE METHOD OF LATENT DEVELOPMENT LEVEL OF EARTHQUAKE-PREGNANT STRUCTURES AND THE ANALYSIS OF RELEVANT SEISMOGEOLOGICAL INFORMATION

    HUANG Xiu-ming

    (Institute of Geology,China Earthquake Administration,Beijing 100029,China)

    The approach of latent development level of earthquake-pregnant structures in the crustal movement(LDLEPS)is a seismogeologic dynamic indicator for evaluating seismic risk quantitatively,developed based on the net deformation ratio using GIS techniques.The following concerns are addressed in this paper,including the working principle and applying conditions for LDLEPS,the drawing of LDLEPS maps,the relationship between the distribution characteristic of LDLEPS and earthquake in the North China region,and the seismogeologic analysis of the LDLEPS anomaly areas in North China.An example is presented of the use in forecasting earthquakes with LDLEPS method,and more discussions are made on the origin of formation of LDLEPS.Contrasting the spatial distribution of LDLEPS with earthquakes of 1966 -1998 in the North China region in Fig 1,we found that they match and correspond completely to each other in terms of lineament,location,and size.This demonstrates that the LDLEPS method can be used in quantitative evaluation of seismic risk.To verify this conclusion,we drawn the LDLEPS maps of Sichuan and Qinghai regions and used them to forecast experimentally the earthquake trend in these regions in 2006,and the afterward outcomes have proved that the forecast is basically correct,as shown in Fig.2 and 3.As a consequence,the 2008 M8 Wenchuan earthquake in Sichuan Province and the 2010 M7.1 Yushu earthquake in Qinghai Province just happened in the expected abnormal LDLEPS regions and times,respectively.Formation mechanism of earthquake-pregnant structures is discussed in detail in the paper.According to the geophysical exploration and seismogeologic survey data,a bilayer-structured earthquake-generation model is proposed,derived from the crustal thermal processes and tectonic reworking resulting from the late-stage evolution of sub-mantle plumes and mantle-branch structures in North China,i.e.the semi-plastic rheologic layer or lower layer in the middle crust and the rigid layer or upper layer in the upper crust,and that the dynamic of the lower layer induces the formation and development of the seismogenic structures in the upper layer,causing tectonic deformation and generating earthquakes.Based on LDLEPS maps,we divide the North China region into 6 seismogenic structure anomaly zones,as shown in Fig.1.Of them,the anomalous seismogenic structure zone B,located in the northwestern Beijing,is rich in geothermal resources,where hot fluids upwell into the fault planes of the upper layer from the lower rheologic layer,reactivating the faults and generating the moderate and lower-magnitude earthquakes,with a causative mechanism similar to that of reservoir-induced earthquake.In the anomaly zone A of seismogenic structures(Bohai Sea area and its margins),the powerful driving force from the lower rheologic layer controls directly the generation and development of the upper seismogenic structures,producing moderate and lower-size earthquakes.Researches on palaeoseismic traces from trench profiles and seismites in marine sediments reveal that there are many strong seismic events with magnitude≥6 and palaeotsunami events in this region in the prehistoric times.The nowadays high level of LDLEPS and strong earthquake events are the continuation and repetition of their past activities.It is also noted that there are obvious differences in structural landforms and earthquake trends between the South and North China blocks owing to difference of their tectonic background.Having the advantage of novelty,uniqueness,visibility,operability,effectiveness and computerization,the LDLEPS method can be widely applied to areas of earth science research,earthquake prediction and engineering.

    GIS,crustal movement,seismic risk,latent development level of earthquake-pregnant structure(LDLEPS),paleoseismology

    P315.2

    A

    0253-4967(2012)03-0500-16

    10.3969/j.issn.0253 - 4967.2012.03.011

    2011-07-05收稿,2011-12-29改回。

    黃秀銘,男,1940年生,1966年畢業(yè)于華東師范大學地貌專業(yè),副研究員,主要從事地震地質(zhì)與新構(gòu)造運動研究,電話:010-62365465,E -mail:HXM1940@Yahoo.com.cn。

    久久久国产成人免费| 天堂俺去俺来也www色官网| 无人区码免费观看不卡| 99re在线观看精品视频| 黄片大片在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| av一本久久久久| 亚洲av日韩精品久久久久久密| 亚洲熟女毛片儿| 亚洲精品乱久久久久久| 久久中文字幕人妻熟女| 丝袜在线中文字幕| 制服人妻中文乱码| 男女下面插进去视频免费观看| 怎么达到女性高潮| 在线观看66精品国产| 亚洲成a人片在线一区二区| 亚洲精品美女久久久久99蜜臀| 国产亚洲av高清不卡| www.精华液| 国产精品 国内视频| 天天影视国产精品| 超碰成人久久| 亚洲免费av在线视频| av在线播放免费不卡| 欧美日韩瑟瑟在线播放| 国产精品亚洲av一区麻豆| 女人久久www免费人成看片| 超碰97精品在线观看| av天堂久久9| a在线观看视频网站| 99精品欧美一区二区三区四区| 精品国产超薄肉色丝袜足j| 夫妻午夜视频| 日韩三级视频一区二区三区| 黄网站色视频无遮挡免费观看| 美女扒开内裤让男人捅视频| 精品国产乱码久久久久久男人| 嫩草影视91久久| 丁香六月欧美| 午夜精品久久久久久毛片777| 成年人免费黄色播放视频| 亚洲成人免费电影在线观看| 午夜成年电影在线免费观看| 最近最新中文字幕大全免费视频| 亚洲国产毛片av蜜桃av| 国产成+人综合+亚洲专区| 夜夜夜夜夜久久久久| 成人亚洲精品一区在线观看| 精品电影一区二区在线| 国产99久久九九免费精品| 中文字幕人妻熟女乱码| 99在线人妻在线中文字幕 | av一本久久久久| 国产成人一区二区三区免费视频网站| 久久精品亚洲av国产电影网| a在线观看视频网站| 久久精品成人免费网站| 国产色视频综合| 18禁黄网站禁片午夜丰满| 国产成人影院久久av| 十八禁人妻一区二区| 一级作爱视频免费观看| 国产极品粉嫩免费观看在线| 99久久精品国产亚洲精品| 国产精品av久久久久免费| 日韩成人在线观看一区二区三区| 久久精品亚洲精品国产色婷小说| 亚洲综合色网址| 18禁黄网站禁片午夜丰满| av线在线观看网站| 啦啦啦视频在线资源免费观看| 国产精品免费大片| 成人18禁高潮啪啪吃奶动态图| 国产精品久久电影中文字幕 | 黄色丝袜av网址大全| 亚洲av电影在线进入| 五月开心婷婷网| 99精国产麻豆久久婷婷| 欧美激情 高清一区二区三区| 在线观看免费视频日本深夜| xxxhd国产人妻xxx| 国产精品久久电影中文字幕 | 久久国产精品人妻蜜桃| 窝窝影院91人妻| 亚洲色图综合在线观看| 国产高清videossex| 狂野欧美激情性xxxx| 大香蕉久久成人网| 在线看a的网站| 天堂中文最新版在线下载| 国产免费现黄频在线看| 国产不卡av网站在线观看| 午夜福利一区二区在线看| 国产精品香港三级国产av潘金莲| 一级a爱片免费观看的视频| 美女福利国产在线| 自拍欧美九色日韩亚洲蝌蚪91| 无人区码免费观看不卡| 日本撒尿小便嘘嘘汇集6| 男女下面插进去视频免费观看| 视频区欧美日本亚洲| 天堂√8在线中文| 国产免费av片在线观看野外av| avwww免费| 国产一区在线观看成人免费| 国产在线一区二区三区精| 精品国产超薄肉色丝袜足j| 99久久99久久久精品蜜桃| 成年版毛片免费区| netflix在线观看网站| 精品卡一卡二卡四卡免费| 校园春色视频在线观看| 69精品国产乱码久久久| 日韩人妻精品一区2区三区| 亚洲av电影在线进入| 91精品三级在线观看| 精品国产一区二区三区四区第35| 国产精品影院久久| 动漫黄色视频在线观看| 欧美日韩av久久| 久久精品亚洲精品国产色婷小说| 欧美最黄视频在线播放免费 | 亚洲熟妇熟女久久| 人妻一区二区av| 他把我摸到了高潮在线观看| av在线播放免费不卡| 搡老岳熟女国产| 久久久国产欧美日韩av| 亚洲精品在线观看二区| 99国产极品粉嫩在线观看| 日韩人妻精品一区2区三区| 高清毛片免费观看视频网站 | 一级毛片精品| 最新的欧美精品一区二区| 麻豆国产av国片精品| 夜夜夜夜夜久久久久| 亚洲人成伊人成综合网2020| 黄片小视频在线播放| 亚洲国产毛片av蜜桃av| 老汉色∧v一级毛片| 女人高潮潮喷娇喘18禁视频| 国产99白浆流出| 国产激情欧美一区二区| 男女高潮啪啪啪动态图| 女人被狂操c到高潮| 国产成人啪精品午夜网站| 黄色怎么调成土黄色| 精品高清国产在线一区| 亚洲一区高清亚洲精品| 欧美丝袜亚洲另类 | 精品国产亚洲在线| 亚洲中文字幕日韩| 亚洲精品在线观看二区| 午夜免费观看网址| 18禁国产床啪视频网站| 99在线人妻在线中文字幕 | 精品国产国语对白av| 一二三四在线观看免费中文在| 精品高清国产在线一区| 狂野欧美激情性xxxx| 国产精品美女特级片免费视频播放器 | 久久久久久久久免费视频了| 久久久久精品国产欧美久久久| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩高清在线视频| 一区二区三区国产精品乱码| 一本综合久久免费| 变态另类成人亚洲欧美熟女 | 热99久久久久精品小说推荐| 在线观看一区二区三区激情| 久久亚洲精品不卡| 成年人免费黄色播放视频| 久久精品人人爽人人爽视色| 黄色视频,在线免费观看| 午夜福利免费观看在线| 叶爱在线成人免费视频播放| 18禁美女被吸乳视频| 香蕉丝袜av| 狠狠婷婷综合久久久久久88av| 99精国产麻豆久久婷婷| 丝袜美腿诱惑在线| 一夜夜www| 久久亚洲真实| 久久精品91无色码中文字幕| 国产成人影院久久av| 久久精品亚洲av国产电影网| 69精品国产乱码久久久| 精品无人区乱码1区二区| 国产麻豆69| 精品久久蜜臀av无| 亚洲精品国产一区二区精华液| 欧美精品av麻豆av| 亚洲情色 制服丝袜| 日日爽夜夜爽网站| 人妻丰满熟妇av一区二区三区 | 亚洲欧美一区二区三区久久| 欧美日韩精品网址| 亚洲伊人色综图| 中文字幕精品免费在线观看视频| 老司机亚洲免费影院| 国产人伦9x9x在线观看| 免费日韩欧美在线观看| 国产成人欧美| 国产一区有黄有色的免费视频| 日本欧美视频一区| 久久精品国产清高在天天线| 免费在线观看完整版高清| 国产av精品麻豆| 中亚洲国语对白在线视频| 久久九九热精品免费| 久久久久久久久久久久大奶| 久久久久国内视频| 精品熟女少妇八av免费久了| 欧美大码av| 欧美激情久久久久久爽电影 | 国产欧美日韩一区二区三区在线| 韩国精品一区二区三区| 99精国产麻豆久久婷婷| 日韩大码丰满熟妇| 国产成人免费观看mmmm| 美女 人体艺术 gogo| 熟女少妇亚洲综合色aaa.| 韩国av一区二区三区四区| 黄色 视频免费看| 另类亚洲欧美激情| 久久99一区二区三区| 国产单亲对白刺激| 日日摸夜夜添夜夜添小说| 女性被躁到高潮视频| 人人妻,人人澡人人爽秒播| 亚洲男人天堂网一区| 亚洲av第一区精品v没综合| 如日韩欧美国产精品一区二区三区| 十分钟在线观看高清视频www| 我的亚洲天堂| 亚洲三区欧美一区| 国产在线观看jvid| 80岁老熟妇乱子伦牲交| 啦啦啦免费观看视频1| 99热只有精品国产| 亚洲五月天丁香| 婷婷丁香在线五月| 精品第一国产精品| 国产激情欧美一区二区| 女性被躁到高潮视频| 免费在线观看影片大全网站| 欧美在线黄色| 国产午夜精品久久久久久| 母亲3免费完整高清在线观看| 99热国产这里只有精品6| 成熟少妇高潮喷水视频| 国产成人精品久久二区二区免费| 亚洲成人手机| 成人国产一区最新在线观看| 国产91精品成人一区二区三区| 一本一本久久a久久精品综合妖精| 色94色欧美一区二区| 国产精品 国内视频| 无人区码免费观看不卡| 色尼玛亚洲综合影院| 色综合婷婷激情| 亚洲精品成人av观看孕妇| 久久久久久久国产电影| 两个人免费观看高清视频| 国产成人精品在线电影| 最新的欧美精品一区二区| 亚洲成a人片在线一区二区| 制服人妻中文乱码| 韩国av一区二区三区四区| 女同久久另类99精品国产91| 精品亚洲成a人片在线观看| 精品乱码久久久久久99久播| 黄频高清免费视频| 在线av久久热| 国产成人免费无遮挡视频| 91在线观看av| 国产亚洲精品久久久久久毛片 | 国产成人精品久久二区二区91| 欧美日韩福利视频一区二区| 欧美另类亚洲清纯唯美| 咕卡用的链子| 一区二区三区激情视频| 精品乱码久久久久久99久播| 亚洲五月天丁香| 性少妇av在线| 亚洲色图综合在线观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美+亚洲+日韩+国产| 国产三级黄色录像| 18禁黄网站禁片午夜丰满| 中文字幕高清在线视频| 悠悠久久av| 日韩欧美一区二区三区在线观看 | 久久性视频一级片| 黄片播放在线免费| 精品一区二区三区视频在线观看免费 | 成人亚洲精品一区在线观看| 母亲3免费完整高清在线观看| 青草久久国产| 黑人欧美特级aaaaaa片| 亚洲欧美激情在线| 美女高潮喷水抽搐中文字幕| 久久性视频一级片| √禁漫天堂资源中文www| 搡老熟女国产l中国老女人| 悠悠久久av| 搡老岳熟女国产| 午夜亚洲福利在线播放| 久久久久国产一级毛片高清牌| 香蕉久久夜色| 免费少妇av软件| 久久99一区二区三区| 亚洲精品国产一区二区精华液| 成人永久免费在线观看视频| 国产精品一区二区精品视频观看| tube8黄色片| 久久久国产精品麻豆| 亚洲精品自拍成人| 精品午夜福利视频在线观看一区| 国产成+人综合+亚洲专区| 国产深夜福利视频在线观看| 90打野战视频偷拍视频| 久99久视频精品免费| √禁漫天堂资源中文www| 精品国内亚洲2022精品成人 | 每晚都被弄得嗷嗷叫到高潮| 国产蜜桃级精品一区二区三区 | 亚洲avbb在线观看| 久久精品亚洲精品国产色婷小说| 男女免费视频国产| 首页视频小说图片口味搜索| 亚洲午夜理论影院| 嫩草影视91久久| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品美女久久av网站| 天堂中文最新版在线下载| 精品久久久久久电影网| 久久 成人 亚洲| 亚洲国产欧美日韩在线播放| 黄片小视频在线播放| 中国美女看黄片| 色综合欧美亚洲国产小说| 国产在视频线精品| 人人妻,人人澡人人爽秒播| 国产三级黄色录像| 久久香蕉激情| 人妻丰满熟妇av一区二区三区 | 一区福利在线观看| 日本a在线网址| 久久草成人影院| 一本大道久久a久久精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品粉嫩美女一区| 欧美成狂野欧美在线观看| 久久ye,这里只有精品| 亚洲精品一卡2卡三卡4卡5卡| 日韩精品免费视频一区二区三区| 精品国产乱码久久久久久男人| 可以免费在线观看a视频的电影网站| 91在线观看av| 亚洲国产毛片av蜜桃av| 久久天躁狠狠躁夜夜2o2o| 美女扒开内裤让男人捅视频| 亚洲中文字幕日韩| 欧美黑人欧美精品刺激| 国产伦人伦偷精品视频| 操美女的视频在线观看| 精品久久久精品久久久| 新久久久久国产一级毛片| 国产精品久久久久久人妻精品电影| 国产片内射在线| 精品电影一区二区在线| 最近最新中文字幕大全免费视频| 免费久久久久久久精品成人欧美视频| 精品久久久精品久久久| 久久精品熟女亚洲av麻豆精品| 国产精品成人在线| 一区二区三区国产精品乱码| 美女午夜性视频免费| 超色免费av| 亚洲欧美日韩高清在线视频| 亚洲欧美激情在线| 亚洲片人在线观看| 午夜福利在线免费观看网站| 99精品久久久久人妻精品| 亚洲va日本ⅴa欧美va伊人久久| 少妇被粗大的猛进出69影院| 免费观看a级毛片全部| 欧美乱色亚洲激情| 在线观看免费视频网站a站| 午夜福利在线观看吧| 看免费av毛片| 久久久久久久国产电影| 久久国产精品大桥未久av| 大码成人一级视频| 国产精品一区二区精品视频观看| 欧美激情 高清一区二区三区| 日本wwww免费看| 91九色精品人成在线观看| 亚洲成人免费av在线播放| 一边摸一边抽搐一进一小说 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一区二区三区欧美精品| 在线播放国产精品三级| 国产精品亚洲一级av第二区| 国产免费现黄频在线看| 夫妻午夜视频| 国产一区二区三区视频了| 99热只有精品国产| 国产极品粉嫩免费观看在线| 亚洲午夜精品一区,二区,三区| 人人妻人人爽人人添夜夜欢视频| x7x7x7水蜜桃| 一级黄色大片毛片| www.999成人在线观看| 国产成人一区二区三区免费视频网站| av电影中文网址| 80岁老熟妇乱子伦牲交| 午夜精品久久久久久毛片777| 757午夜福利合集在线观看| 99re6热这里在线精品视频| 少妇的丰满在线观看| 亚洲成人免费电影在线观看| 国产午夜精品久久久久久| 日本wwww免费看| 一边摸一边抽搐一进一出视频| 欧美激情 高清一区二区三区| 精品熟女少妇八av免费久了| 欧美人与性动交α欧美软件| 国产欧美日韩一区二区三| 国产欧美日韩综合在线一区二区| 免费在线观看视频国产中文字幕亚洲| 一区二区日韩欧美中文字幕| 午夜精品久久久久久毛片777| 国产精品 国内视频| av网站在线播放免费| 亚洲欧美激情在线| 女性生殖器流出的白浆| 国产97色在线日韩免费| av免费在线观看网站| 午夜免费观看网址| 国产精品美女特级片免费视频播放器 | 国产极品粉嫩免费观看在线| 亚洲国产精品合色在线| 天堂√8在线中文| 777米奇影视久久| 水蜜桃什么品种好| 91精品三级在线观看| 校园春色视频在线观看| 老司机在亚洲福利影院| 午夜福利免费观看在线| 成人18禁高潮啪啪吃奶动态图| 亚洲成人免费av在线播放| 国产精华一区二区三区| 亚洲午夜理论影院| 精品国产美女av久久久久小说| 人人妻人人澡人人看| 中文亚洲av片在线观看爽 | 亚洲专区字幕在线| 国产高清激情床上av| a级毛片在线看网站| 777久久人妻少妇嫩草av网站| 女人久久www免费人成看片| 久久久久国产一级毛片高清牌| 老熟女久久久| a在线观看视频网站| 亚洲va日本ⅴa欧美va伊人久久| 亚洲熟妇熟女久久| 丰满人妻熟妇乱又伦精品不卡| 自线自在国产av| 嫁个100分男人电影在线观看| 精品国产国语对白av| 久久久精品免费免费高清| 国产成人av教育| 午夜激情av网站| 啦啦啦 在线观看视频| 欧美国产精品一级二级三级| 99久久人妻综合| 香蕉国产在线看| 在线观看66精品国产| 黑丝袜美女国产一区| 国产aⅴ精品一区二区三区波| 日韩欧美国产一区二区入口| 最近最新免费中文字幕在线| 国产亚洲精品一区二区www | 色尼玛亚洲综合影院| 国产一区有黄有色的免费视频| 久9热在线精品视频| 欧美人与性动交α欧美软件| 国产欧美日韩精品亚洲av| 高清在线国产一区| 国产一区在线观看成人免费| 极品人妻少妇av视频| 一区二区三区激情视频| 亚洲av熟女| 久久久久国内视频| 一级毛片女人18水好多| 1024香蕉在线观看| 国产男女内射视频| 午夜精品国产一区二区电影| 中亚洲国语对白在线视频| 一级毛片精品| 欧美老熟妇乱子伦牲交| 99久久综合精品五月天人人| 极品少妇高潮喷水抽搐| 制服诱惑二区| 一边摸一边做爽爽视频免费| 69av精品久久久久久| 美国免费a级毛片| 国产欧美日韩精品亚洲av| 大片电影免费在线观看免费| 一进一出抽搐gif免费好疼 | 亚洲九九香蕉| 国产99久久九九免费精品| 久久天堂一区二区三区四区| 国产男女内射视频| 亚洲成人免费av在线播放| 黄色片一级片一级黄色片| 国产精品久久久久久精品古装| 老鸭窝网址在线观看| 国产在视频线精品| 国产免费现黄频在线看| 高潮久久久久久久久久久不卡| 在线观看免费视频网站a站| 久久久精品国产亚洲av高清涩受| 黑人巨大精品欧美一区二区mp4| 久久久久久久精品吃奶| 国产成+人综合+亚洲专区| 人人妻人人澡人人爽人人夜夜| 精品一区二区三区av网在线观看| 国产不卡一卡二| 亚洲免费av在线视频| 18禁观看日本| 大香蕉久久成人网| 国产精品一区二区精品视频观看| 人人妻人人澡人人爽人人夜夜| 久久人人爽av亚洲精品天堂| 久久亚洲精品不卡| 少妇的丰满在线观看| 日韩欧美国产一区二区入口| 免费观看a级毛片全部| 麻豆乱淫一区二区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品国产精品久久久不卡| 欧美久久黑人一区二区| 一a级毛片在线观看| 久久久久久亚洲精品国产蜜桃av| 999久久久国产精品视频| 在线观看免费午夜福利视频| 国产高清视频在线播放一区| 波多野结衣一区麻豆| 飞空精品影院首页| 精品熟女少妇八av免费久了| 国产成人精品久久二区二区免费| 天堂√8在线中文| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 日本vs欧美在线观看视频| 美女视频免费永久观看网站| 一区在线观看完整版| 国产高清视频在线播放一区| 99精品久久久久人妻精品| 午夜视频精品福利| 在线观看免费午夜福利视频| 少妇粗大呻吟视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久久久久精品古装| 免费看十八禁软件| 国产精品亚洲一级av第二区| 高清毛片免费观看视频网站 | 亚洲精品国产一区二区精华液| 午夜影院日韩av| aaaaa片日本免费| 久久精品国产a三级三级三级| 国产欧美日韩综合在线一区二区| 国产精品久久久久久精品古装| 啪啪无遮挡十八禁网站| 免费日韩欧美在线观看| 在线看a的网站| 亚洲色图av天堂| 精品福利永久在线观看| 亚洲专区国产一区二区| 18在线观看网站| 在线观看www视频免费| 亚洲人成电影观看| 亚洲精品国产区一区二| 国产欧美日韩一区二区精品| 丝袜在线中文字幕| 在线观看免费午夜福利视频| 久久久久久免费高清国产稀缺| 中国美女看黄片| 久久香蕉精品热| 欧美人与性动交α欧美精品济南到| 国产亚洲一区二区精品| 成人黄色视频免费在线看| 日韩免费av在线播放| 色婷婷av一区二区三区视频| 丁香欧美五月| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩一区二区三| 久久久久久久久免费视频了| 99精品欧美一区二区三区四区| 亚洲欧美日韩高清在线视频| 捣出白浆h1v1| 亚洲欧美日韩高清在线视频| 欧美激情极品国产一区二区三区| 看黄色毛片网站| 精品国产乱码久久久久久男人| 不卡av一区二区三区| 99国产精品一区二区蜜桃av | 国产精品一区二区在线观看99| 美国免费a级毛片| 国产成人精品无人区| 老司机在亚洲福利影院|