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

    鄂爾多斯南緣地區(qū)重力變化場源特征

    2021-10-13 14:31:24張永奇韓美濤鄭增記曹建平
    地震地磁觀測與研究 2021年3期
    關(guān)鍵詞:場源重力場歐拉

    張永奇 韓美濤 鄭增記 曹建平

    1)中國西安 710068 陜西省地震局2)中國北京 100029 中國地震局地質(zhì)研究所地震動力學(xué)國家重點實驗室

    0 引言

    地震孕育過程中產(chǎn)生的地殼變形和震源介質(zhì)變化在地表會引起重力變化,通過開展連續(xù)重力觀測或定期相對重力復(fù)測,有可能捕捉到地震的這些前兆信息。重力場具有明確的地球物理含義,屬于構(gòu)造運動的一種表現(xiàn)形式,一般分為“靜態(tài)”重力場和“動態(tài)”重力場?!办o態(tài)”重力場是指,基于地球重力模型或?qū)崪y重力數(shù)據(jù),獲取布格重力及均衡重力異常信息,用來研究地質(zhì)構(gòu)造特征、地殼物性密度變化及地殼均衡狀態(tài)等地質(zhì)時間尺度的地學(xué)問題(陳石等,2011a,2011b,2014);“動態(tài)”重力場是指,在年時間尺度內(nèi),利用流動作業(yè)方式獲取的重力異常變化信息,主要用來認識強震孕育機理、捕捉地震前兆信息和判定地震危險區(qū)(祝意青等,2009,2017,2020;申重陽等2010,2011;陳石等,2015;Chen et al,2016a)。

    當(dāng)前主要基于“場”模式開展地震監(jiān)測預(yù)測,而應(yīng)用地表觀測的重力場變化信息來開展異常源的研究,即“以場求源”問題的研究仍不夠深入。眾所周知,地表觀測的重力場變化是測量誤差、地表水影響、測點三維形變、地下介質(zhì)密度變化及物質(zhì)遷移等信息的綜合反映(陳運泰等,1980;陳石等,2011c)。準(zhǔn)確分離并識別這些因素對重力場變化信息的影響,是比較困難的事情,原因如下:①深部介質(zhì)變化的先驗信息及參數(shù)難以觀測或獲?。虎诂F(xiàn)有重力儀器觀測精度偏低,信噪比不高。此外,在重力場變化的相關(guān)研究中,經(jīng)常采用大地測量學(xué)方法進行數(shù)據(jù)計算和解釋,間接造成了上述困難無法克服。為此,嘗試將三維歐拉反褶積方法引入重力場變化分析。該方法是Reid 等(1990)結(jié)合Thompson(1982)在歐拉齊次方程的研究基礎(chǔ)上建立起來的,其最大優(yōu)勢是,在先驗信息缺乏情況下,仍能較好完成位場數(shù)據(jù)反演與解釋工作(范美寧,2006)。

    文中對鄂爾多斯南緣地區(qū)EGM2008 重力場模型數(shù)據(jù)及2014—2017 年相對重力觀測數(shù)據(jù)進行精細化處理。在此基礎(chǔ)上,以布格重力異常、1 年尺度差分重力場和3 年尺度累積重力場為研究對象,分析上述重力場異常變化的場源特征。設(shè)計一個與實際重力觀測相似的理論模型模擬地表重力場變化,同時利用三維歐拉反褶積方法對場源參數(shù)進行反演,優(yōu)選與模型接近的構(gòu)造指數(shù)及滑動窗口大小等參數(shù),對相對重力觀測網(wǎng)數(shù)據(jù)進行歐拉反演,并利用水平梯度濾波法對發(fā)散的反演結(jié)果進行優(yōu)化,結(jié)合地震目錄資料,對獲取的三維場源信息進行分析與解釋。

    1 測區(qū)概況與數(shù)據(jù)處理

    鄂爾多斯南緣地區(qū)位于我國青藏塊體與華北塊體交匯區(qū)域,地質(zhì)構(gòu)造形式多樣,動力學(xué)環(huán)境復(fù)雜。同時,該地區(qū)位于我國重要的汾渭地震帶和南北地震帶北段,歷史上曾發(fā)生多次6 級以上大震,屬于我國地震活動重點監(jiān)視區(qū)。該區(qū)發(fā)育多條活動斷裂,其中鄂爾多斯南緣斷裂、渭河斷裂及秦嶺北緣斷裂控制著該區(qū)主要構(gòu)造變形演化格局,區(qū)內(nèi)其他斷裂進一步把鄂爾多斯南緣地區(qū)分隔為不同樣式的地質(zhì)構(gòu)造單元,單條或多條斷裂邊界或交匯處成為地震孕育的主要場所。為了監(jiān)測鄂爾多斯南緣地區(qū)地震活動潛在危險性,尤其是渭河斷陷盆地內(nèi)主要斷裂的地球物理場異常變化,陜西省地震局于20 世紀(jì)80 年代建立了相對重力觀測網(wǎng),每年觀測1 期。2008 年汶川8.0 級地震后,依托“震后恢復(fù)重建”等項目,對原有相對重力觀測網(wǎng)進行優(yōu)化改造,改造后的重力測網(wǎng)覆蓋陜西全境,加上“陸態(tài)網(wǎng)絡(luò)”一、二期建立的重力站點,重力站點總計達200 余座,點間距基本在20—30 km,站點分布見圖1。這些重力站點覆蓋渭河盆地及鄰區(qū)大部分?jǐn)嗔?,且分布相對合理,為監(jiān)測斷裂活動引起的重力變化奠定了堅實的基礎(chǔ)(曹建平等,2014;張永奇等,2014)。

    圖1 鄂爾多斯南緣地區(qū)重力點及斷裂分布Fig.1 Distribution of gravity points and faults in the southern margin of the Ordos

    基于重力資料開展鄂爾多斯南緣地區(qū)的研究,以往集中在數(shù)據(jù)精細化處理、重力場變化及重力網(wǎng)的映震能力分析等方面,針對深部場源特征的研究相對較少。因此,嘗試基于EGM2008 重力場模型數(shù)據(jù)及相對重力數(shù)據(jù),利用歐拉反褶積方法,對該區(qū)進行場源特征反演。反演前將相對重力觀測數(shù)據(jù)進行平差計算,具體流程為:①采用中國地震局實用化攻關(guān)推廣軟件——LGADJ,對多期重力觀測資料進行獨立平差計算,解算重力儀一次項周期誤差等參數(shù)的區(qū)域適定解,計算各臺儀器觀測精度,合理確定其先驗方差;②以四川廣元、西安子午、甘肅天水等絕對測點為起算基準(zhǔn),采用經(jīng)典平差方法進行平差計算,對其中因周圍環(huán)境改變而出現(xiàn)異常的個別測點進行改正或刪除;③采用克里金插值算法擬合離散重力場,獲得不同時間尺度的重力場。

    2 方法原理

    歐拉反褶積方法不依賴場源密度參數(shù),僅與場源位置的幾何參數(shù)相關(guān),使得反演問題得以簡化(Reid et al,1990)。通過求解歐拉方程,可以確定場源位置的三維空間位置,并可對不同類型地質(zhì)體給出參數(shù)估計,公式如下

    構(gòu)造指數(shù)N反映了位場異常強度隨場源深度變化的衰減率,一般與場源幾何結(jié)構(gòu)有關(guān)。N的選取是場源位置確定的關(guān)鍵。根據(jù)不同構(gòu)造特征,針對重力場異常數(shù)據(jù),N的范圍一般取0—2 較為合適,其中斷層、巖脈、臺階等構(gòu)造指數(shù)皆不同(范美寧,2006;高倩等,2015;劉芳等,2016;翟麗娜等,2020)。三維重力歐拉反褶積方法的反演流程如下:①對實測離散重力位場數(shù)據(jù)進行網(wǎng)格化;②計算重力位場異常在x、y、z三個方向的導(dǎo)數(shù);③選擇適當(dāng)?shù)臉?gòu)造指數(shù)N及窗口大小。根據(jù)異常體性質(zhì)選擇構(gòu)造指數(shù)N,窗口大小根據(jù)異常體大小進行選擇,一般應(yīng)覆蓋異常范圍,不得低于3 倍測點間距。根據(jù)精度要求確定允許誤差;④解方程組,求出場源位置參數(shù)及背景值;⑤根據(jù)解的匯聚位置、深度調(diào)整參數(shù)。如果對反演結(jié)果滿意,結(jié)束計算,否則重新選擇參數(shù)重復(fù)步驟②、③,直至對結(jié)果滿意為止(高倩等,2015;劉芳等,2016;翟麗娜等,2020)。

    3 模型反演

    歐拉反褶積方法不需要已知場源先驗信息,即可對場源參數(shù)進行反演與解釋。反演計算時,應(yīng)選擇合適的構(gòu)造指數(shù)N,并求出異常在x、y、z三個方向的導(dǎo)數(shù),即可組成線性方程組進行求解(范美寧,2006)。在實際應(yīng)用中,地質(zhì)情況復(fù)雜,場源形態(tài)也不確定,所以如何選擇構(gòu)造指數(shù),需要對不同模型進行反復(fù)計算來確定(魯寶亮等,2009)。此外,滑動窗口大小的選取,一般需考慮測點間距、異常形態(tài)特征等,在一個窗口內(nèi)應(yīng)有足夠多的測點滿足方程組的求解,并要覆蓋一個場源異常形態(tài)。對于反演結(jié)果的評價,需要通過求得場源位置參數(shù)的水平誤差和垂直誤差大小,及反演場源位置解的三維空間聚集程度和與地質(zhì)認識是否符合等判斷,即場源位置越集中,其精度及收斂性越好,反之解的收斂性差(高倩等,2015;劉芳等,2016)。

    基于對歐拉反褶積法的認識,根據(jù)鄂爾多斯地塊南緣地區(qū)現(xiàn)有重力觀測網(wǎng)形態(tài)及規(guī)模,設(shè)計一個相對符合實際觀測情況的理論模型。以3 個直面正六面體組合模型為研究模型,假設(shè)其分別位于上地殼、中地殼和下地殼。其中模型1(Model1)位于上地殼,單元尺寸為20 km × 20 km × 20 km,中心埋深為5.5 km;模型2(Model2)位于中地殼,單元尺寸為10 km×10 km×10 km,中心埋深為11.5 km;模型3(Model3)位于下地殼,單元尺寸為30 km×30 km×30 km,中心埋深為22.5 km。3 個模型密度變化統(tǒng)一取為3.0 g·cm-3,正演計算模型在地表重力測網(wǎng)的理論異常值,用于模擬觀測的重力變化信號和優(yōu)選反演參數(shù),模型異常見圖2,矩形框表示場源模型在地表的水平投影。地表重力測網(wǎng)點間距設(shè)計為5 km,測網(wǎng)長度在X、Y方向均為200 km、深度方向為50 km,3 個不同深度的場源模型參數(shù)見表1,可見3個模型的尺度特征不同,埋深也不同,正演結(jié)果顯示理論異常范圍在(0—120)×10-8m·s-2。

    圖2 不同深度場源組合模型及正演異常Fig.2 Combination models of different depth field sources and forward modeling anomalies

    表1 理論場源模型參數(shù)Table 1 The theoretical field source model parameters

    基于上述理論測試結(jié)果,實際測試時設(shè)計多個參數(shù)組合,其中構(gòu)造指數(shù)分別選擇0.5、1、1.5、2,滑動窗口分別選擇4—8 倍的測點距離。通過對多次試驗結(jié)果進行對比分析,認為選擇構(gòu)造指數(shù)N=1,滑動窗口大小在4 個測點距離時,反演結(jié)果的聚集程度更符合模型的幾何結(jié)構(gòu),組合模型的反演深度誤差范圍限制在10%以內(nèi)。圖3 給出計算所得一組最佳反演結(jié)果,圖中彩色圓圈為反演的場源參數(shù),共獲得467 組解,黑色矩形框為理論組合模型位置,可見聚集程度相對較好,場源特征明顯。其中模型2(Model2)和模型3(Model3)反演結(jié)果最好,模型1(Model1)反演結(jié)果稍差。

    圖3 基于歐拉反褶積的反演結(jié)果空間分布Fig.3 Spatial distribution of inversion results based on Euler deconvolution

    4 歐拉反演解優(yōu)化

    歐拉反演的基本條件是模型滿足歐拉齊次方程的假設(shè)條件,此時反演結(jié)果應(yīng)該是準(zhǔn)確的。然而實際上地質(zhì)體的位置和構(gòu)成情況復(fù)雜,無法用任何一種模型準(zhǔn)確表達,這就會導(dǎo)致反演解出現(xiàn)發(fā)散現(xiàn)象(姚長利等,2004;管志寧等,2005;石磊等,2008),因此需要尋找一種減弱或消除反演解發(fā)散的方法,研究表明水平梯度濾波法是一種比較有效的方法。在利用歐拉反褶積方法開展反演研究時,當(dāng)窗口滑動到異常體邊緣時,梯度值變小甚至趨近于零,此時反演容易出現(xiàn)發(fā)散解。理論上,水平梯度決定了場源的水平位置,構(gòu)造指數(shù)決定了場源的深度,因此構(gòu)造指數(shù)越大,反演的深度越深,結(jié)果將更加發(fā)散(姚長利等,2004;鄭秋月等,2021),使用水平梯度濾波法即可克服上述問題。

    設(shè)Hi=為研究區(qū)任意測點i處的重、磁異常水平梯度模,則該區(qū)域內(nèi)重磁異常的水平梯度模為n為研究區(qū)網(wǎng)格點總數(shù)),以Hi≥kHA為濾波準(zhǔn)則(k為常數(shù))。若k值太小,將無法剔除發(fā)散解區(qū)域;若k值太大,會導(dǎo)致有效場源信息丟失。一般,k值在[0.5,1.5]范圍內(nèi)選擇(姚長利等,2004)。

    水平梯度濾波法優(yōu)勢在于:①可剔除水平梯度較小的區(qū)域網(wǎng)格點,選擇水平梯度模相對較大的區(qū)域參與反演;②可去除因測網(wǎng)不均勻插值帶來的虛假網(wǎng)格點信息(姚長利等,2004;管志寧等,2005;鄭秋月等,2021)。文中在計算場源位置時,采用水平梯度濾波法,對傳統(tǒng)的歐拉反演結(jié)果予以優(yōu)化,其中k值取為1。

    5 重力場變化及場源特征反演與解釋

    5.1 鄂爾多斯南緣地區(qū)“靜態(tài)”重力場變化特征

    地震孕育發(fā)生是區(qū)域構(gòu)造應(yīng)力積累與釋放的過程,震源位置和震源機制是巖石圈動力學(xué)過程的反映,通過地震可以了解地下深部的地質(zhì)構(gòu)造和應(yīng)力狀態(tài),同時可為分析重力異常變化與地震分布的關(guān)系提供便利。為此,統(tǒng)計鄂爾多斯南緣地區(qū)歷史地震目錄,繪制歷史地震分布圖,見圖4。

    圖4 鄂爾多斯南緣地區(qū)地震分布Fig.4 Earthquake distribution in the southern margin of the Ordos area

    重力位場異常可以表征地質(zhì)構(gòu)造運動過程,常用于劃分構(gòu)造單元、識別斷裂形態(tài)和盆地邊界(馬宗晉等,2006)。一般認為,自由空氣重力異??梢员碚鞯叵挛镔|(zhì)盈虧程度;布格重力異??梢苑从成畈繕?gòu)造特征,被稱為構(gòu)造運動的“活化石”(陳石等,2011b)。在分析鄂爾多斯南緣地區(qū)“動態(tài)”重力場變化之前,分析“靜態(tài)”重力場變化特征?;贓GM2008 地球重力模型,計算自由空氣重力異常和布格重力異常(Sandwell et al,2009;Pavlis et al,2012),結(jié)果見圖5。

    圖5 鄂爾多斯地塊南緣地區(qū)重力異常(a)自由空氣異常;(b)布格重力異常Fig.5 Gravity anomaly in the southern margin of the Ordos block

    圖5(a)為自由空氣重力異常,異常范圍基本在(-140—110)×10-5m·s-2。該異常與地形具有較強的相關(guān)性,在高海拔的山地,重力異常表現(xiàn)為正值,在低海拔的平原和盆地,重力異常表現(xiàn)為負值,自由空氣異?;痉从沉嗽搮^(qū)構(gòu)造單元的大體輪廓。圖5(b)為布格重力異常,研究區(qū)整體表現(xiàn)為負值異常,范圍在(-240—-80)×10-5m·s-2。研究區(qū)西北端靠近青藏高原東北緣地區(qū),其布格重力異常呈現(xiàn)極小負值,向SE,重力異常呈逐漸增大趨勢,異常整體形態(tài)與斷裂及構(gòu)造區(qū)分布具有一定對應(yīng)關(guān)系。結(jié)合圖4 所示的地震分布,發(fā)現(xiàn)大部分地震震中與斷裂帶發(fā)育位置具有較好的一致性,但有些地震與斷裂之間的對應(yīng)關(guān)系不明顯。此外,重力異常的變化特征與地震分布也具有一定相關(guān)性。進一步說明,重力異常變化、斷裂活動與地震孕育發(fā)生具有較好的協(xié)調(diào)性和一致性。

    5.2 “靜態(tài)”重力場歐拉反演與解釋

    為了全面掌握鄂爾多斯南緣地區(qū)重力異常變化的場源特征,對該區(qū)布格重力異常進行歐拉反演計算。反演前,對布格重力異常進行坐標(biāo)轉(zhuǎn)換,將大地坐標(biāo)轉(zhuǎn)換為平面直角坐標(biāo),采用克里格插值算法進行插值,并采用高斯濾波進行平滑。在此基礎(chǔ)上,對布格重力異常數(shù)據(jù)進行歐拉反演,策略如下:①構(gòu)造指數(shù)取1.5;②X、Y方向點距均設(shè)置為10 km;③水平梯度濾波倍數(shù)k=0.5;④滑動窗口設(shè)為4—8 倍點距。由此獲得基于布格重力異常的鄂爾多斯南緣地區(qū)場源特征,見圖6。從“靜態(tài)場”角度分析鄂爾多斯南緣地區(qū)場源特征,由圖6 可知,在4—8 倍的歐拉反演結(jié)果中,場源分布的水平位置和深度基本相同,聚集性相對較好的區(qū)域集中在斷裂邊界帶附近,尤其是研究區(qū)西部的海原—六盤山地區(qū)以及南北地震帶北部地區(qū),異常源聚集性較好且與斷裂走向基本一致,呈NS 向分布。渭河盆地內(nèi)部場源呈EW 向分布特征,主要沿斷裂走向展布。在鄂爾多斯地塊內(nèi)部出現(xiàn)平?jīng)觥獩艽ê蛻c城—富縣地區(qū)2 個場源特征較顯著區(qū)域。秦嶺造山帶場源集中在勉—略大斷裂、月河斷裂和山陽斷裂附近。從場源深度看,0—30 km 深度的場源占整個反演結(jié)果的90%以上,表明該區(qū)布格重力異常變化主要由淺部異常源引起。結(jié)合圖4 可知,地震分布的三維位置與基于布格重力異常計算的場源位置吻合程度較高。

    圖6 基于布格重力異常的鄂爾多斯南緣地區(qū)場源特征Fig.6 Field source characteristics in the southern margin of Ordos based on Bouguer gravity anomaly

    5.3 鄂爾多斯南緣地區(qū)“動態(tài)”重力場變化特征

    選取研究區(qū)2014—2017 年相對重力數(shù)據(jù),計算“動態(tài)”重力異常變化,結(jié)果見圖7。其中(a)—(c)圖分別表示2014年8月—2015年8月、2015年8月—2016年7月、2016年7月—2017 年7 月的1 年尺度重力差分變化結(jié)果,(d)圖表示2014 年8 月—2017 年7 月3 年尺度的重力累積變化結(jié)果。圖中灰色線表示斷裂帶,藍色等值線表示負變化,紅色等值線表示正變化,黑色等值線表示“0”值線,等值線間距10×10-8m·s-2。因部分地區(qū)重力站點分布不均,采用克里格插值算法,對平差后重力數(shù)據(jù)進行插值及高斯濾波,等值線覆蓋范圍與測網(wǎng)范圍基本一致。在測點稀疏地區(qū),重力場變化的可信度偏低,不做深入分析。

    由圖7(a)可見,2014 年8 月—2015 年8 月,該區(qū)重力變化基本以西安為界,以西大部分地區(qū)以負值變化為主,以東地區(qū)以正值變化為主。在成縣—太白斷裂與秦嶺北緣斷裂交匯處出現(xiàn)局部重力正值變化,量值約20×10-8m·s-2。另一個重力正值變化區(qū)位于海原—六盤山斷裂帶以東固原、慶城一帶,量值在30×10-8m·s-2。以渭河斷裂與長安—臨潼斷裂交界為中心出現(xiàn)不規(guī)則重力四象限分布特征,整體來看,重力場展布特征受主要斷裂控制。由圖7(b)可見,2015 年8 月—2016 年7 月,研究區(qū)重力以負值變化為主,在西秦嶺北緣斷裂、六盤山—寶雞斷裂帶附近出現(xiàn)重力變化梯度帶,差異變化達70×10-8m·s-2。其他地區(qū)重力場變化相對平穩(wěn),變化量基本在(-40— -10)×10-8m·s-2。由圖7(c)可見,2016 年7 月—2017 年7 月重力場變化相對復(fù)雜。與圖6(a)相比,重力場出現(xiàn)反向變化,基本以西安為界,以西地區(qū)重力表現(xiàn)為大范圍正值變化,而東部則為負值變化。在六盤山—寶雞斷裂附近不僅存在重力變化梯度帶,還具有四象限分布特征。在秦嶺中部及鄂爾多斯東南部地區(qū)也出現(xiàn)了重力變化梯度帶。由圖7(d)可見,2014 年8 月—2017年7 月,研究區(qū)重力變化相對平穩(wěn),且以負值變化為主,變化量基本在(-40— -10)×10-8m·s-2,僅在海原—六盤山斷裂以西地區(qū)出現(xiàn)局部正值變化,變化量最大為50×10-8m·s-2,可見3 年尺度重力累積變化不明顯,未出現(xiàn)大震、中強震的孕育特征。

    由圖7 可見,1 年尺度的重力場變化在某些區(qū)域存在一定繼承性,但大部分地區(qū)的重力變化存在明顯差異。海原—六盤山—隴縣斷裂帶附近重力變化梯度帶位置在不同時間尺度上基本一致,但是變化量存在差異,重力變化等值線的展布和斷裂分布具有一定對應(yīng)關(guān)系。研究表明:半年尺度,甚至1 年尺度的重力變化依然存在正負波動的變化特征。這種波動可能由近地表局部影響、觀測誤差所導(dǎo)致,目前尚無法有效分離各種場源引起的重力變化。理論上,若為深部場源引起的區(qū)域重力變化,其在不同時間尺度內(nèi)應(yīng)具有相似性和一致性,此區(qū)域即為比較可靠的異常源位置,也是物質(zhì)匯集、應(yīng)力積聚區(qū)域,未來發(fā)生地震的可能性相對較大。

    5.4 “動態(tài)”重力場歐拉反演結(jié)果與解釋

    基于常規(guī)歐拉反褶積方法,對不同時間尺度的重力變化進行反演,尋找聚集性好、場源重復(fù)性好的區(qū)域。此類區(qū)域往往物質(zhì)積累相對集中、構(gòu)造活動頻繁。實際反演時設(shè)定網(wǎng)格化間距為10 km,構(gòu)造指數(shù)為1,窗口大小為5—8 倍測點間距。反演完成后,對反演結(jié)果進行統(tǒng)計,發(fā)現(xiàn)反演解個數(shù)存在差異,具體結(jié)果見表2。

    表2 改正前后基于不同滑動窗口獲得的反演結(jié)果個數(shù)Table 2 The number of inversion results obtained based on different sliding windows before and after correction

    對不同反演結(jié)果進行對比分析,將6 倍測點間距獲得的結(jié)果作為最優(yōu)解來進行分析,具體反演結(jié)果見圖8(a)。圖中圓圈表示反演的場源位置,圓圈大小與場源深度相對應(yīng),采用紅、黃、綠、藍4 種顏色表示2014—2017 年不同時間段的反演結(jié)果。由圖7 可知,從引起重力場變化的不同時間尺度結(jié)果來看,場源位置相對分散,無明顯規(guī)律性。總體而言,研究區(qū)邊緣地區(qū)場源深度反演結(jié)果相對較深,內(nèi)部區(qū)域場源則相對較淺。在不同時間尺度下,海原—六盤山—隴縣斷裂帶、成縣—太白斷裂、華山山前斷裂以及韓城斷裂附近地區(qū)場源位置集聚性較高,而鄂爾多斯地塊南部及秦嶺造山帶南部場源位置相對分散。作為一種地震前異常信號,場源位置的分散型分布特征表明,鄂爾多斯南緣地區(qū)近年來不具備中強地震孕育與發(fā)生的條件,同時說明,2014—2017 年重力場變化屬于一種隨機的、正常的波動變化,地震緊迫性相對較弱。此外,統(tǒng)計得到不同深度的反演解個數(shù),具體結(jié)果見表3。

    表3 改正前后不同深度的反演結(jié)果個數(shù)Table 3 The number of inversion results at different depths before and after correction

    由表3 可知,2014—2017 年,發(fā)生重力變化的場源深度集中在0—50 km,其中場源深度為5—30 km 的反演結(jié)果約占60%。對比發(fā)現(xiàn),2014 年8 月—2015 年8 月,場源深度0—10 km 的反演解數(shù)量較少,深度50—100 km 的反演解數(shù)量較多,說明該期重力場變化受深部物質(zhì)影響較大。在該時間范圍之前曾于2013 年7 月22 日發(fā)生岷縣—漳縣MS6.6 地震,該地震距研究區(qū)西部僅150 km,因此重力變化可能受到地震遠程應(yīng)力調(diào)整及物質(zhì)運移的影響。其他3 期的反演結(jié)果顯示,大部分場源深度在30 km 以上,與該區(qū)震源深度基本吻合。

    利用常規(guī)歐拉反演方法可獲得比較合理的場源特征,但反演結(jié)果發(fā)散性嚴(yán)重,采用水平梯度濾波法對反演結(jié)果進行優(yōu)化,發(fā)現(xiàn)歐拉反演結(jié)果更加聚集收斂,發(fā)散現(xiàn)象顯著減弱,場源特征更加明顯,具體結(jié)果見圖8(b),與圖8(a)對比可知,研究區(qū)邊緣大部分區(qū)域反演結(jié)果數(shù)量明顯減少。這是因為,水平梯度法刪除了測點稀疏造成的不真實反演結(jié)果??傮w來看,研究區(qū)西北部的海原—六盤山地區(qū)出現(xiàn)場源叢集現(xiàn)象,鄂爾多斯地塊內(nèi)部場源集聚性較好地區(qū)分布在慶城一帶。渭河斷陷盆地內(nèi)部反演結(jié)果集聚性較好區(qū)域基本分布在主要斷裂附近,且場源展布與斷裂走向基本一致。秦嶺造山帶的場源特征相對分散,且場源位置相對較深。同時,統(tǒng)計水平梯度濾波反演結(jié)果,與常規(guī)歐拉反演結(jié)果(表2,表3)進行對比,可知:從5—8 倍測點間距滑動窗口結(jié)果看,采用水平梯度濾波法效果較好,反演解數(shù)量較常規(guī)歐拉反演結(jié)果至少減少50%,剔除了大部分不真實的場源信息;采用水平梯度濾波法處理后,對場源深度反演結(jié)果改正效果顯著,進一步表明相對重力觀測結(jié)果主要反映深度在30 km 以內(nèi)的異常源信息。

    圖8 鄂爾多斯南緣地區(qū)2014—2017 年重力變化場源特征(a)常規(guī)歐拉反演結(jié)果;(b)水平梯度濾波結(jié)果Fig.8 Field source characteristics of gravity changes in the southern margin of Ordos from 2014 to 2017

    5.5 重力變化場源特征與地震的關(guān)系

    鄂爾多斯南緣地區(qū)在地質(zhì)構(gòu)造劃分上主要包含鄂爾多斯地塊南部、渭河斷陷盆地及秦嶺造山帶,在地震帶劃分上主要包含汾渭地震帶的渭河地震帶及南北地震帶北部。該研究區(qū)歷史上曾發(fā)生多次6 級以上強震,具備強震發(fā)生的歷史背景和構(gòu)造條件。據(jù)歷史地震記錄和震源機制解結(jié)果,鄂爾多斯南緣地區(qū)震源深度相對較淺,多在30 km 以內(nèi)(王衛(wèi)東,2009;杜建軍,2016)。對圖4、圖6、圖8 及表2 進行細致分析,發(fā)現(xiàn)基于布格重力異常數(shù)據(jù)計算的“靜態(tài)”場源位置和基于相對重力觀測數(shù)據(jù)計算的“動態(tài)”場源位置,均與研究區(qū)歷史地震分布位置具有高度一致性,表明場源集聚性較好的地區(qū)有利于地震的孕育與發(fā)生(劉芳等,2016;鄭秋月等,2021)。本研究結(jié)果表明,海原—六盤山斷裂帶、隴縣—寶雞斷裂帶、渭河斷裂帶西段以及華山山前斷裂東段場源積聚性較好,未來發(fā)生地震的可能性較高。其他構(gòu)造區(qū)異常源呈現(xiàn)相對性和隨機性,地震危險性相對較弱。此外,據(jù)歐拉反演深度,重力場源集中在30 km 以內(nèi),位于中上地殼,與王衛(wèi)東(2009)、杜建軍(2016)等研究所得震源深度一致。

    6 結(jié)論

    作為我國東、西部不同構(gòu)造單元的連接部位,鄂爾多斯南緣地區(qū)不僅是物質(zhì)流通的重要通道,也是青藏塊體和華北塊體相互作用的關(guān)鍵部位。受2 種不同構(gòu)造應(yīng)力的影響,該區(qū)山地隆升,盆地下陷,斷裂發(fā)育,地震頻發(fā),歷史上發(fā)生多次6 級以上強震,甚至8 級巨震。因此,鄂爾多斯南緣地區(qū)的地震危險性歷來是地學(xué)研究的重點關(guān)注問題之一。近年來,隨著地震觀測手段的日益豐富、觀測技術(shù)的不斷提高,對孕震機理、發(fā)震構(gòu)造、前兆信號的研究日趨成熟,并取得長足進展(陳石等,2011c,2014,2015;祝意青等,2017,2020),如基于重力觀測資料開展地震監(jiān)測預(yù)測,取得豐碩成果(祝意青等,2009;申重陽等,2010;Chen et al,2016a,2016b)。但是,大多數(shù)工作是基于重力異?!皥觥遍_展的,針對重力異?!霸础遍_展的工作相對較少,主要原因是缺少足夠的先驗信息。為此,本研究嘗試采用在地球物理勘探工作中應(yīng)用廣泛的三維反褶積方法,進行重力變化場源特征研究,取得以下認識。

    (1)利用歐拉反褶積方法進行反演,關(guān)鍵在于選擇合適的構(gòu)造指數(shù)、滑動窗口大小等參數(shù)。通過構(gòu)建模型,正演得到基本可以反映真實重力變化的模型參數(shù)。通過不同實驗方案進行測試,發(fā)現(xiàn)當(dāng)構(gòu)造指數(shù)為1,滑動窗口為4—8 倍點位間距時,反演效果最佳。

    (2)反演結(jié)果表明:海原—六盤山—隴縣斷裂、隴縣—寶雞斷裂、渭河斷裂和華山山前斷裂附近,在不同時間尺度下場源特征一致性較好。結(jié)合圖4,可知以上地區(qū)也是中小地震頻發(fā)區(qū)域,即更易于應(yīng)力積累、能量積聚而發(fā)生地震的位置。由圖8、表3 可知,場源反演深度集中在5—30 km 范圍內(nèi),屬中上地殼。據(jù)地震震源機制解,可知該深度范圍也是鄂爾多斯南緣地區(qū)的主要震源深度。

    (3)歐拉反褶積方法適用于反演重力場變化異常的場源信息,不同的場源位置及深度揭示了不同的構(gòu)造信息,多期重力變化反演三維場源特征一致性較好,表明該區(qū)存在物質(zhì)積聚和應(yīng)力增強;反之,則該區(qū)地震危險性較弱。圖8(a)顯示,在不同時間尺度下,場源位置相對分散,表明大部分地區(qū)的重力變化是隨機的物質(zhì)運動現(xiàn)象,屬于較短期的地殼活動現(xiàn)象或者與觀測誤差不確定性有關(guān)。但圖8(b)顯示,經(jīng)水平梯度濾波的歐拉反演結(jié)果聚集性和收斂性更強,場源特征更加明顯。因此,在實際反演時,需要對歐拉反演結(jié)果做進一步濾波處理。

    需要說明的是,選擇2014—2017 年相對重力觀測數(shù)據(jù)開展研究,主要是因為,2013 年陜西重力測網(wǎng)進行了優(yōu)化改造,增建大量站點,而2018 年進行了二次優(yōu)化改造,刪減了大量站點。因此,陜西重力測網(wǎng)在該時間范圍內(nèi)點位密度相對較高,分布較均勻。此外,本研究在建立正演模型時,僅考慮了直立六面體組合,實際地質(zhì)構(gòu)造更復(fù)雜,而文中測試模型相對簡單。盡管如此,本研究結(jié)果仍能驗證歐拉反褶積方法的可靠性和適用性。希望通過本研究,為鄂爾多斯南緣地區(qū)地震危險性研判提供一種新的研究思路和技術(shù)途徑。

    文中流動重力數(shù)據(jù)由國家重力臺網(wǎng)中心提供,撰寫過程中多次與云南省地震局鄭秋月工程師進行討論,匿名審稿專家及編輯部老師對本文提出的修改意見與建議,對文章質(zhì)量的提升大有益處,在此一并表示感謝!

    猜你喜歡
    場源重力場歐拉
    例談求解疊加電場的電場強度的策略
    歐拉閃電貓
    汽車觀察(2022年12期)2023-01-17 02:20:42
    基于深度展開ISTA網(wǎng)絡(luò)的混合源定位方法
    信號處理(2022年10期)2022-11-16 00:50:56
    歐拉魔盒
    哈哈畫報(2022年1期)2022-04-19 11:27:20
    精致背后的野性 歐拉好貓GT
    車迷(2022年1期)2022-03-29 00:50:26
    基于矩陣差分的遠場和近場混合源定位方法
    基于空間分布的重力場持續(xù)適配能力評估方法
    歐拉的疑惑
    衛(wèi)星測量重力場能力仿真分析
    一種識別位場場源的混合小波方法
    男插女下体视频免费在线播放| 日韩免费av在线播放| 日本熟妇午夜| 曰老女人黄片| 日韩 欧美 亚洲 中文字幕| 男女床上黄色一级片免费看| 日韩欧美免费精品| 欧美黑人巨大hd| 麻豆成人午夜福利视频| 欧美极品一区二区三区四区| 国产精品一区二区免费欧美| 成人永久免费在线观看视频| 老汉色∧v一级毛片| 色在线成人网| 欧美av亚洲av综合av国产av| 精品乱码久久久久久99久播| 亚洲成人中文字幕在线播放| 人人妻人人澡欧美一区二区| 丰满的人妻完整版| 噜噜噜噜噜久久久久久91| 亚洲欧美一区二区三区黑人| 黄片大片在线免费观看| 狂野欧美白嫩少妇大欣赏| 性欧美人与动物交配| 99国产精品一区二区蜜桃av| 老司机午夜十八禁免费视频| 久久久久性生活片| 最新在线观看一区二区三区| www.熟女人妻精品国产| 香蕉久久夜色| 日本黄色视频三级网站网址| 人妻夜夜爽99麻豆av| 国内少妇人妻偷人精品xxx网站 | 国产伦精品一区二区三区视频9 | 97碰自拍视频| 国产成人系列免费观看| 国产91精品成人一区二区三区| 舔av片在线| 久久中文字幕一级| 欧美成人一区二区免费高清观看 | 又黄又爽又免费观看的视频| 国产亚洲精品久久久久久毛片| 亚洲欧美精品综合一区二区三区| 不卡一级毛片| 日本五十路高清| 国产精品av视频在线免费观看| 国产精品,欧美在线| 午夜福利免费观看在线| 成人av在线播放网站| 亚洲在线观看片| 999久久久精品免费观看国产| 日韩欧美三级三区| 国产精品美女特级片免费视频播放器 | 99视频精品全部免费 在线 | 国产精品综合久久久久久久免费| 欧美精品啪啪一区二区三区| 我的老师免费观看完整版| 国产精品香港三级国产av潘金莲| 久久久国产精品麻豆| 嫁个100分男人电影在线观看| 精品熟女少妇八av免费久了| 黄色丝袜av网址大全| 怎么达到女性高潮| 特级一级黄色大片| 亚洲欧美精品综合一区二区三区| 999久久久精品免费观看国产| 精品久久久久久,| 国产精品九九99| 久久精品91无色码中文字幕| 国产欧美日韩一区二区精品| 女人被狂操c到高潮| 精品国产乱子伦一区二区三区| 老熟妇乱子伦视频在线观看| ponron亚洲| 色吧在线观看| 最近最新中文字幕大全免费视频| 狂野欧美白嫩少妇大欣赏| 国内久久婷婷六月综合欲色啪| 国产av麻豆久久久久久久| 黄色日韩在线| 嫩草影视91久久| 岛国视频午夜一区免费看| 精品久久久久久久毛片微露脸| 久久久久国内视频| 欧美色欧美亚洲另类二区| 午夜精品一区二区三区免费看| 一个人看的www免费观看视频| 人人妻,人人澡人人爽秒播| 亚洲av成人av| av视频在线观看入口| 在线免费观看不下载黄p国产 | 日韩欧美国产一区二区入口| 丰满的人妻完整版| 日本免费a在线| 一进一出抽搐动态| 最新美女视频免费是黄的| 国产毛片a区久久久久| 国产av在哪里看| 久久久国产精品麻豆| 欧美色视频一区免费| 欧美国产日韩亚洲一区| 中文字幕高清在线视频| 国产精品国产高清国产av| 欧美日本亚洲视频在线播放| 哪里可以看免费的av片| 日韩高清综合在线| 日本三级黄在线观看| netflix在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 国产蜜桃级精品一区二区三区| 狠狠狠狠99中文字幕| 日韩有码中文字幕| 淫妇啪啪啪对白视频| 日本成人三级电影网站| www日本黄色视频网| 国产亚洲av高清不卡| 久久久久亚洲av毛片大全| 91麻豆av在线| 国产精品电影一区二区三区| 色老头精品视频在线观看| 亚洲美女黄片视频| 叶爱在线成人免费视频播放| 精品一区二区三区视频在线 | 国产欧美日韩精品一区二区| 在线十欧美十亚洲十日本专区| 一个人看的www免费观看视频| 村上凉子中文字幕在线| 九色成人免费人妻av| 国产精品久久视频播放| 欧美成人免费av一区二区三区| 欧美乱色亚洲激情| 久久欧美精品欧美久久欧美| 亚洲18禁久久av| 欧美在线一区亚洲| 人妻久久中文字幕网| 国产成+人综合+亚洲专区| 久久久久久九九精品二区国产| 日韩 欧美 亚洲 中文字幕| 露出奶头的视频| 1024手机看黄色片| 99热6这里只有精品| 日韩欧美在线二视频| 曰老女人黄片| 久久久国产欧美日韩av| 欧美日韩中文字幕国产精品一区二区三区| 真人一进一出gif抽搐免费| 一本精品99久久精品77| 99视频精品全部免费 在线 | 俺也久久电影网| av黄色大香蕉| 在线观看免费视频日本深夜| 男女做爰动态图高潮gif福利片| 一夜夜www| 亚洲一区高清亚洲精品| 亚洲无线观看免费| 亚洲精品一区av在线观看| 精品久久久久久久毛片微露脸| 国产成人一区二区三区免费视频网站| 国产欧美日韩精品一区二区| 国产精品综合久久久久久久免费| 国产av不卡久久| 久久久国产成人免费| 久久久久久久久免费视频了| 欧美三级亚洲精品| 久久天堂一区二区三区四区| 岛国在线免费视频观看| 中出人妻视频一区二区| 亚洲欧美日韩东京热| 女同久久另类99精品国产91| 国产男靠女视频免费网站| 免费电影在线观看免费观看| 国产精品爽爽va在线观看网站| 美女 人体艺术 gogo| 麻豆一二三区av精品| 999久久久国产精品视频| 琪琪午夜伦伦电影理论片6080| 成人av在线播放网站| 真人做人爱边吃奶动态| 悠悠久久av| 国产精品九九99| x7x7x7水蜜桃| 韩国av一区二区三区四区| 久久久国产欧美日韩av| 男女之事视频高清在线观看| 欧美黄色淫秽网站| 欧美精品啪啪一区二区三区| 亚洲乱码一区二区免费版| 亚洲18禁久久av| 嫁个100分男人电影在线观看| 亚洲国产中文字幕在线视频| 免费无遮挡裸体视频| 日韩欧美一区二区三区在线观看| 欧美激情在线99| 少妇的丰满在线观看| 国产精品乱码一区二三区的特点| 国产真人三级小视频在线观看| 日韩av在线大香蕉| 蜜桃久久精品国产亚洲av| 一二三四在线观看免费中文在| 国产精品一及| 久久久久国内视频| 淫妇啪啪啪对白视频| 成人av一区二区三区在线看| 亚洲精品国产精品久久久不卡| 无人区码免费观看不卡| 看黄色毛片网站| 神马国产精品三级电影在线观看| 露出奶头的视频| 欧美中文日本在线观看视频| 欧美激情久久久久久爽电影| 久久亚洲精品不卡| a在线观看视频网站| 国产在线精品亚洲第一网站| 国产一区二区三区视频了| 亚洲国产欧美人成| 亚洲人成伊人成综合网2020| 91在线观看av| 欧美av亚洲av综合av国产av| 日本a在线网址| 日韩欧美在线二视频| 久久精品91蜜桃| 日韩三级视频一区二区三区| 在线观看免费午夜福利视频| 国产亚洲精品综合一区在线观看| 国产男靠女视频免费网站| 两个人视频免费观看高清| 老司机福利观看| 国产亚洲av高清不卡| 一二三四社区在线视频社区8| 国产亚洲精品久久久com| 国产探花在线观看一区二区| 成人国产一区最新在线观看| 国产精品av久久久久免费| 国产激情久久老熟女| 久久国产精品影院| 一二三四社区在线视频社区8| 老司机午夜十八禁免费视频| 不卡av一区二区三区| 色噜噜av男人的天堂激情| 性色avwww在线观看| 国产激情偷乱视频一区二区| 国产av在哪里看| 久久中文字幕一级| 成年女人毛片免费观看观看9| 网址你懂的国产日韩在线| 18禁国产床啪视频网站| 真人一进一出gif抽搐免费| 日日摸夜夜添夜夜添小说| 大型黄色视频在线免费观看| 一个人免费在线观看的高清视频| 欧美日韩福利视频一区二区| 国产成人精品久久二区二区免费| 日本五十路高清| 91在线精品国自产拍蜜月 | 亚洲第一欧美日韩一区二区三区| 亚洲最大成人中文| 久久婷婷人人爽人人干人人爱| 久久午夜亚洲精品久久| 美女被艹到高潮喷水动态| 人人妻,人人澡人人爽秒播| 岛国在线免费视频观看| 天堂影院成人在线观看| 国产午夜精品论理片| 国内揄拍国产精品人妻在线| 日本免费a在线| 久久天躁狠狠躁夜夜2o2o| 日本在线视频免费播放| 久久久久国产一级毛片高清牌| 国产伦一二天堂av在线观看| 亚洲成人精品中文字幕电影| 制服丝袜大香蕉在线| 91麻豆av在线| av福利片在线观看| 亚洲av成人av| 性欧美人与动物交配| 日本成人三级电影网站| 国产精品九九99| 久久久久国产一级毛片高清牌| 老熟妇乱子伦视频在线观看| 久久婷婷人人爽人人干人人爱| 亚洲国产色片| 亚洲一区二区三区不卡视频| 男女下面进入的视频免费午夜| 一本一本综合久久| 91九色精品人成在线观看| 午夜a级毛片| 国产欧美日韩一区二区精品| 亚洲国产欧美网| 国产不卡一卡二| 又大又爽又粗| 国内精品久久久久精免费| 国产三级中文精品| 叶爱在线成人免费视频播放| 午夜福利成人在线免费观看| 精品国产超薄肉色丝袜足j| 91麻豆精品激情在线观看国产| 好男人在线观看高清免费视频| 伦理电影免费视频| 欧美三级亚洲精品| 极品教师在线免费播放| 女人被狂操c到高潮| 成年女人看的毛片在线观看| 90打野战视频偷拍视频| 精品日产1卡2卡| 亚洲国产高清在线一区二区三| 亚洲,欧美精品.| 国产麻豆成人av免费视频| 91在线观看av| 欧美日本亚洲视频在线播放| a级毛片在线看网站| 国产成人啪精品午夜网站| 不卡一级毛片| 女生性感内裤真人,穿戴方法视频| 老司机午夜十八禁免费视频| 757午夜福利合集在线观看| 他把我摸到了高潮在线观看| 国产精品自产拍在线观看55亚洲| 精品午夜福利视频在线观看一区| 婷婷亚洲欧美| 亚洲av免费在线观看| 在线免费观看不下载黄p国产 | 老熟妇乱子伦视频在线观看| 天堂动漫精品| 最新在线观看一区二区三区| 亚洲av五月六月丁香网| 19禁男女啪啪无遮挡网站| 亚洲第一欧美日韩一区二区三区| 身体一侧抽搐| 久久久久久久精品吃奶| 可以在线观看毛片的网站| avwww免费| 狂野欧美白嫩少妇大欣赏| 久久国产精品影院| 欧美性猛交╳xxx乱大交人| 成人18禁在线播放| 久久欧美精品欧美久久欧美| 一本久久中文字幕| 99国产精品一区二区三区| 亚洲av免费在线观看| 美女扒开内裤让男人捅视频| 久久久久久九九精品二区国产| 精品久久久久久久末码| 欧美乱码精品一区二区三区| 日本精品一区二区三区蜜桃| 亚洲av电影不卡..在线观看| 成人av一区二区三区在线看| 欧美黑人巨大hd| 国产成+人综合+亚洲专区| 国产欧美日韩精品亚洲av| 老司机福利观看| 久久久久久九九精品二区国产| 午夜福利在线观看吧| 久久久久免费精品人妻一区二区| 深夜精品福利| 亚洲av电影不卡..在线观看| 日韩欧美国产一区二区入口| 观看美女的网站| 精品国内亚洲2022精品成人| 国内揄拍国产精品人妻在线| 免费在线观看成人毛片| 国产一区二区三区视频了| 夜夜爽天天搞| 国产1区2区3区精品| 欧美午夜高清在线| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 嫩草影院入口| netflix在线观看网站| 长腿黑丝高跟| 黄色丝袜av网址大全| 一级a爱片免费观看的视频| 国产美女午夜福利| 小说图片视频综合网站| 欧美性猛交黑人性爽| 午夜福利成人在线免费观看| 人妻夜夜爽99麻豆av| 嫩草影院入口| 特级一级黄色大片| 国产亚洲欧美在线一区二区| 高清在线国产一区| 国产成人系列免费观看| 香蕉国产在线看| 成年女人毛片免费观看观看9| 久久久久久人人人人人| www.www免费av| 天天一区二区日本电影三级| 成人特级av手机在线观看| 男女那种视频在线观看| 巨乳人妻的诱惑在线观看| 三级毛片av免费| 动漫黄色视频在线观看| 99国产精品一区二区蜜桃av| 亚洲无线观看免费| 国产成人福利小说| 欧美大码av| 少妇丰满av| 可以在线观看的亚洲视频| 成年版毛片免费区| 香蕉丝袜av| 国产伦精品一区二区三区四那| 18禁国产床啪视频网站| 精品99又大又爽又粗少妇毛片 | 男插女下体视频免费在线播放| 激情在线观看视频在线高清| 欧美日韩乱码在线| 午夜日韩欧美国产| 午夜免费观看网址| 亚洲av电影不卡..在线观看| 色在线成人网| 亚洲成a人片在线一区二区| ponron亚洲| 欧美乱妇无乱码| 99视频精品全部免费 在线 | 熟女电影av网| 日韩av在线大香蕉| 国产精品亚洲av一区麻豆| 免费在线观看影片大全网站| 男插女下体视频免费在线播放| 老司机在亚洲福利影院| 亚洲av第一区精品v没综合| 中文资源天堂在线| 变态另类丝袜制服| 亚洲在线观看片| 久久精品国产清高在天天线| 欧美中文综合在线视频| 美女免费视频网站| 国产精品亚洲av一区麻豆| 午夜免费成人在线视频| 一进一出抽搐动态| 两个人的视频大全免费| 男插女下体视频免费在线播放| av福利片在线观看| 免费搜索国产男女视频| 少妇熟女aⅴ在线视频| 亚洲va日本ⅴa欧美va伊人久久| 九九在线视频观看精品| 老鸭窝网址在线观看| 欧美激情久久久久久爽电影| 偷拍熟女少妇极品色| 波多野结衣巨乳人妻| 日韩三级视频一区二区三区| 亚洲av成人精品一区久久| 91麻豆精品激情在线观看国产| 久久伊人香网站| 国产精品影院久久| bbb黄色大片| 好男人电影高清在线观看| 亚洲人成电影免费在线| 999精品在线视频| 久久久国产成人免费| 国内精品美女久久久久久| 成年人黄色毛片网站| 亚洲中文字幕日韩| 国产精品久久视频播放| 国产精品99久久99久久久不卡| 欧美黑人巨大hd| 成年版毛片免费区| 久久亚洲真实| 日韩精品中文字幕看吧| 狠狠狠狠99中文字幕| 亚洲欧美日韩东京热| 岛国在线免费视频观看| www日本在线高清视频| 搡老妇女老女人老熟妇| 99国产精品一区二区三区| 亚洲av第一区精品v没综合| 老司机深夜福利视频在线观看| 少妇裸体淫交视频免费看高清| 国产成+人综合+亚洲专区| 制服丝袜大香蕉在线| 一个人看视频在线观看www免费 | 一个人观看的视频www高清免费观看 | 日本熟妇午夜| 国产黄色小视频在线观看| 搞女人的毛片| 91av网站免费观看| 一级毛片高清免费大全| 欧美日韩中文字幕国产精品一区二区三区| 免费一级毛片在线播放高清视频| 国产三级黄色录像| 在线国产一区二区在线| 成在线人永久免费视频| 亚洲专区字幕在线| 听说在线观看完整版免费高清| 国产亚洲精品一区二区www| 在线观看美女被高潮喷水网站 | 这个男人来自地球电影免费观看| 又紧又爽又黄一区二区| 亚洲五月婷婷丁香| 国产亚洲av高清不卡| 欧美黑人巨大hd| 欧美日韩精品网址| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩瑟瑟在线播放| 亚洲黑人精品在线| 亚洲最大成人中文| 51午夜福利影视在线观看| 搡老熟女国产l中国老女人| 中文字幕熟女人妻在线| 黑人操中国人逼视频| 免费观看人在逋| 日韩精品中文字幕看吧| 久久精品91无色码中文字幕| 国产精品亚洲一级av第二区| 嫁个100分男人电影在线观看| 999久久久精品免费观看国产| 亚洲人成电影免费在线| 此物有八面人人有两片| 色视频www国产| 在线视频色国产色| 欧美中文综合在线视频| 91麻豆精品激情在线观看国产| 精品久久久久久久久久免费视频| 亚洲熟妇熟女久久| 欧美高清成人免费视频www| 欧美日韩乱码在线| 免费无遮挡裸体视频| 99riav亚洲国产免费| 999精品在线视频| 狂野欧美白嫩少妇大欣赏| 久久久国产欧美日韩av| 亚洲avbb在线观看| avwww免费| 久久伊人香网站| 日本精品一区二区三区蜜桃| 看片在线看免费视频| 最近视频中文字幕2019在线8| 91在线观看av| 亚洲精品中文字幕一二三四区| 网址你懂的国产日韩在线| 俺也久久电影网| 国内少妇人妻偷人精品xxx网站 | 亚洲国产精品sss在线观看| 熟女少妇亚洲综合色aaa.| 99久久精品热视频| 精品国产乱子伦一区二区三区| 亚洲精品粉嫩美女一区| 午夜久久久久精精品| 美女高潮喷水抽搐中文字幕| 丁香六月欧美| 日本免费a在线| 婷婷精品国产亚洲av| 国产亚洲精品久久久com| 成人国产一区最新在线观看| 国产成人欧美在线观看| 黄色女人牲交| 一级毛片女人18水好多| 最近视频中文字幕2019在线8| 国产毛片a区久久久久| 男插女下体视频免费在线播放| 一进一出抽搐gif免费好疼| 亚洲成a人片在线一区二区| 精品国产三级普通话版| 欧美一级a爱片免费观看看| 久久久国产精品麻豆| 不卡av一区二区三区| 99riav亚洲国产免费| 日本免费一区二区三区高清不卡| 最新在线观看一区二区三区| 国产高清激情床上av| 国产成人av激情在线播放| 噜噜噜噜噜久久久久久91| 成人国产综合亚洲| 免费观看的影片在线观看| 人人妻,人人澡人人爽秒播| 精品久久久久久成人av| 国产高清视频在线观看网站| 在线国产一区二区在线| 免费搜索国产男女视频| 我的老师免费观看完整版| 日本黄色片子视频| 在线永久观看黄色视频| 国产激情偷乱视频一区二区| 精品免费久久久久久久清纯| 亚洲熟女毛片儿| 99久久国产精品久久久| 午夜福利欧美成人| 淫妇啪啪啪对白视频| 久久中文字幕人妻熟女| 亚洲va日本ⅴa欧美va伊人久久| 黄片大片在线免费观看| 嫩草影视91久久| 国产精品久久电影中文字幕| 深夜精品福利| 亚洲专区字幕在线| 成人无遮挡网站| 好男人电影高清在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲精品456在线播放app | 久久精品国产99精品国产亚洲性色| 国产精品一及| 欧美高清成人免费视频www| 母亲3免费完整高清在线观看| 一进一出抽搐gif免费好疼| 美女被艹到高潮喷水动态| 精品国产乱码久久久久久男人| 国产亚洲av高清不卡| 亚洲欧美激情综合另类| 久久亚洲精品不卡| 日日摸夜夜添夜夜添小说| 国产成人精品久久二区二区91| 999精品在线视频| 久久草成人影院| 天堂动漫精品| 天天一区二区日本电影三级| 欧美成狂野欧美在线观看| 亚洲国产精品sss在线观看| 淫妇啪啪啪对白视频| 黄频高清免费视频| 国内精品美女久久久久久| 国产单亲对白刺激| 中文字幕人妻丝袜一区二区| 一个人看视频在线观看www免费 | 亚洲avbb在线观看| 国产精品九九99|