• <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在线免费看完整版不卡| 久久狼人影院| av国产精品久久久久影院| 丰满人妻一区二区三区视频av| 高清午夜精品一区二区三区| 欧美成人午夜免费资源| 国产在视频线精品| 夫妻午夜视频| 欧美精品一区二区免费开放| 国产成人freesex在线| 亚洲精品乱码久久久v下载方式| 国产av国产精品国产| 国产黄频视频在线观看| 两个人免费观看高清视频 | 最后的刺客免费高清国语| 成人国产麻豆网| 亚洲国产精品一区二区三区在线| 美女视频免费永久观看网站| 免费播放大片免费观看视频在线观看| 久久鲁丝午夜福利片| 免费黄色在线免费观看| 777米奇影视久久| 大又大粗又爽又黄少妇毛片口| 99热这里只有精品一区| 嫩草影院新地址| 亚洲天堂av无毛| 国产欧美亚洲国产| 丝袜喷水一区| 嘟嘟电影网在线观看| 日本av手机在线免费观看| 亚洲欧洲国产日韩| 日韩欧美 国产精品| 麻豆成人av视频| 十八禁高潮呻吟视频 | 在线 av 中文字幕| 国产亚洲欧美精品永久| 只有这里有精品99| 久久狼人影院| 国产精品久久久久久久久免| 久久精品久久精品一区二区三区| 国产欧美日韩综合在线一区二区 | 久久狼人影院| 国产老妇伦熟女老妇高清| 日日爽夜夜爽网站| 精品国产乱码久久久久久小说| 亚洲在久久综合| 国产高清三级在线| 狂野欧美白嫩少妇大欣赏| 我要看日韩黄色一级片| 99国产精品免费福利视频| 国产精品蜜桃在线观看| 日本av手机在线免费观看| 精品国产一区二区三区久久久樱花| 51国产日韩欧美| 伦理电影大哥的女人| 国产精品国产三级国产专区5o| 欧美国产精品一级二级三级 | 成人国产av品久久久| av播播在线观看一区| 国产在线视频一区二区| videos熟女内射| a 毛片基地| 日韩av免费高清视频| 久久毛片免费看一区二区三区| 欧美丝袜亚洲另类| 性色avwww在线观看| 18禁在线播放成人免费| 深夜a级毛片| 日韩一区二区视频免费看| 黄色毛片三级朝国网站 | 青青草视频在线视频观看| 欧美日韩精品成人综合77777| 精品人妻熟女毛片av久久网站| 国产精品国产三级国产专区5o| 人人澡人人妻人| 亚洲美女视频黄频| 国产免费福利视频在线观看| 插逼视频在线观看| 91精品一卡2卡3卡4卡| 在线观看人妻少妇| 亚洲图色成人| 人人澡人人妻人| 免费久久久久久久精品成人欧美视频 | 青春草亚洲视频在线观看| 街头女战士在线观看网站| 夫妻性生交免费视频一级片| 美女中出高潮动态图| 国内精品宾馆在线| 亚洲自偷自拍三级| 观看av在线不卡| 日韩人妻高清精品专区| 十八禁网站网址无遮挡 | 亚洲国产毛片av蜜桃av| 国产老妇伦熟女老妇高清| 国产精品欧美亚洲77777| 免费观看a级毛片全部| av天堂中文字幕网| 精品久久久噜噜| 91精品国产九色| 国产视频首页在线观看| 亚洲经典国产精华液单| 狠狠精品人妻久久久久久综合| 2021少妇久久久久久久久久久| 成人特级av手机在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲在久久综合| 97在线人人人人妻| 国产黄色视频一区二区在线观看| 中文字幕av电影在线播放| 精品午夜福利在线看| 天堂中文最新版在线下载| a级一级毛片免费在线观看| 亚洲国产欧美日韩在线播放 | 天堂俺去俺来也www色官网| 观看美女的网站| 97超视频在线观看视频| 伊人亚洲综合成人网| 国产探花极品一区二区| 亚洲av成人精品一区久久| 丰满饥渴人妻一区二区三| av在线播放精品| 97超视频在线观看视频| 欧美bdsm另类| 久久久午夜欧美精品| 亚洲精品自拍成人| 91午夜精品亚洲一区二区三区| 午夜福利网站1000一区二区三区| 精品人妻偷拍中文字幕| 久久热精品热| 亚洲精品视频女| 蜜臀久久99精品久久宅男| 搡女人真爽免费视频火全软件| 各种免费的搞黄视频| 五月玫瑰六月丁香| 少妇人妻 视频| 亚洲精品久久久久久婷婷小说| 精品一品国产午夜福利视频| 精品亚洲成a人片在线观看| 欧美xxxx性猛交bbbb| av国产久精品久网站免费入址| 97在线视频观看| 免费不卡的大黄色大毛片视频在线观看| 日韩一本色道免费dvd| 肉色欧美久久久久久久蜜桃| 欧美+日韩+精品| 亚洲精品一二三| av专区在线播放| av在线app专区| 国产精品久久久久久久久免| 国产精品国产三级国产专区5o| 一区二区三区乱码不卡18| 国国产精品蜜臀av免费| 色94色欧美一区二区| 成年女人在线观看亚洲视频| 乱码一卡2卡4卡精品| 日本黄色片子视频| 久久精品久久久久久久性| 男人和女人高潮做爰伦理| 六月丁香七月| 天堂8中文在线网| 视频区图区小说| 亚洲欧美成人综合另类久久久| 欧美日韩精品成人综合77777| 日本av免费视频播放| 亚洲精品日韩在线中文字幕| 亚洲欧美精品专区久久| 边亲边吃奶的免费视频| 男的添女的下面高潮视频| 久久免费观看电影| 欧美精品人与动牲交sv欧美| 26uuu在线亚洲综合色| 晚上一个人看的免费电影| 国产欧美亚洲国产| 一级爰片在线观看| 久久影院123| 色网站视频免费| av在线播放精品| 18禁裸乳无遮挡动漫免费视频| 成人毛片a级毛片在线播放| 国产成人一区二区在线| 欧美最新免费一区二区三区| 亚洲美女黄色视频免费看| 亚洲av.av天堂| 嫩草影院入口| 涩涩av久久男人的天堂| 欧美精品人与动牲交sv欧美| 亚洲欧洲精品一区二区精品久久久 | 国产精品秋霞免费鲁丝片| 久久青草综合色| 99热全是精品| 少妇高潮的动态图| 一本色道久久久久久精品综合| 黄色日韩在线| 在线播放无遮挡| 成年人午夜在线观看视频| 麻豆成人av视频| 精品亚洲成a人片在线观看| 成人二区视频| 午夜福利视频精品| av免费观看日本| 国产精品一区二区性色av| 一级毛片黄色毛片免费观看视频| 中文字幕人妻丝袜制服| 欧美xxⅹ黑人| 我的老师免费观看完整版| 波野结衣二区三区在线| 在线观看免费日韩欧美大片 | 午夜激情久久久久久久| 狂野欧美白嫩少妇大欣赏| 久久久欧美国产精品| 国产精品国产三级专区第一集| 多毛熟女@视频| 久久久久国产精品人妻一区二区| 久久久久久久久久成人| 精品久久久久久电影网| 国产精品一区www在线观看| 人妻 亚洲 视频| 蜜臀久久99精品久久宅男| 最近中文字幕高清免费大全6| 免费看日本二区| 三级经典国产精品| 97超视频在线观看视频| 97精品久久久久久久久久精品| 午夜免费男女啪啪视频观看| 哪个播放器可以免费观看大片| 综合色丁香网| 日日啪夜夜爽| 两个人的视频大全免费| 一级二级三级毛片免费看| 男人狂女人下面高潮的视频| 大陆偷拍与自拍| 最新中文字幕久久久久| 老司机影院成人| 只有这里有精品99| 久久青草综合色| 男人狂女人下面高潮的视频| 一级二级三级毛片免费看| 国产精品久久久久久av不卡| 久久综合国产亚洲精品| 激情五月婷婷亚洲| av一本久久久久| 桃花免费在线播放| 18禁在线播放成人免费| 内地一区二区视频在线| 美女脱内裤让男人舔精品视频| av又黄又爽大尺度在线免费看| 久久久久精品性色| 日韩,欧美,国产一区二区三区| 亚洲美女视频黄频| 国产成人aa在线观看| 成人美女网站在线观看视频| 男女边摸边吃奶| 欧美国产精品一级二级三级 | 国产精品不卡视频一区二区| 国内揄拍国产精品人妻在线| 亚洲精品久久午夜乱码| 丝袜喷水一区| 高清黄色对白视频在线免费看 | 国产精品久久久久久精品电影小说| 欧美变态另类bdsm刘玥| 一区二区av电影网| 亚洲成人一二三区av| 色网站视频免费| √禁漫天堂资源中文www| 黄色视频在线播放观看不卡| www.av在线官网国产| 七月丁香在线播放| 天堂俺去俺来也www色官网| 欧美最新免费一区二区三区| 最近最新中文字幕免费大全7| 婷婷色综合www| 日韩欧美一区视频在线观看 | 免费看日本二区| 亚洲天堂av无毛| 中文字幕制服av| 亚洲国产精品一区二区三区在线| 久久久精品免费免费高清| 久久久久久久国产电影| 久久久a久久爽久久v久久| 日韩欧美一区视频在线观看 | 国产成人精品无人区| 大片电影免费在线观看免费| 十八禁高潮呻吟视频 | 涩涩av久久男人的天堂| 免费播放大片免费观看视频在线观看| 国产精品99久久99久久久不卡 | 五月玫瑰六月丁香| 极品少妇高潮喷水抽搐| 日本午夜av视频| 狂野欧美白嫩少妇大欣赏| 熟女电影av网| 国模一区二区三区四区视频| 国产精品.久久久| av一本久久久久| 亚洲国产毛片av蜜桃av| 国产成人91sexporn| 国产精品久久久久久精品古装| 国产女主播在线喷水免费视频网站| 国精品久久久久久国模美| 18+在线观看网站| 桃花免费在线播放| 久久精品国产a三级三级三级| 亚洲精品久久久久久婷婷小说| 日本欧美国产在线视频| 亚洲av不卡在线观看| 女性被躁到高潮视频| 高清午夜精品一区二区三区| 久久精品国产亚洲av涩爱| 国产高清三级在线| 亚洲,欧美,日韩| 在线看a的网站| 国产 一区精品| 五月玫瑰六月丁香| 在线观看国产h片| 涩涩av久久男人的天堂| 91aial.com中文字幕在线观看| 曰老女人黄片| 日本-黄色视频高清免费观看| 午夜福利视频精品| 精品酒店卫生间| 中国国产av一级| 日韩一区二区视频免费看| 亚洲一级一片aⅴ在线观看| 中文字幕亚洲精品专区| 五月伊人婷婷丁香| 香蕉精品网在线| 2022亚洲国产成人精品| 久久精品国产亚洲网站| 亚洲精品,欧美精品| 亚洲精品成人av观看孕妇| 边亲边吃奶的免费视频| 欧美精品亚洲一区二区| 亚洲激情五月婷婷啪啪| 99热全是精品| 街头女战士在线观看网站| 韩国av在线不卡| 三级国产精品欧美在线观看| 夜夜骑夜夜射夜夜干| 成人亚洲欧美一区二区av| 毛片一级片免费看久久久久| 久久久久久久久久人人人人人人| 欧美日韩国产mv在线观看视频| h日本视频在线播放| 91在线精品国自产拍蜜月| 熟女av电影| 有码 亚洲区| 久久久久久久国产电影| 免费看光身美女| 日韩在线高清观看一区二区三区| 插逼视频在线观看| 国产成人免费观看mmmm| 欧美xxxx性猛交bbbb| 黄色日韩在线| 国产视频首页在线观看| 精品亚洲成国产av| 午夜日本视频在线| 国产高清有码在线观看视频| 欧美精品高潮呻吟av久久| 亚洲真实伦在线观看| 丰满少妇做爰视频| 久热这里只有精品99| 国产一区二区在线观看日韩| 免费人妻精品一区二区三区视频| 男男h啪啪无遮挡| 精品久久久久久久久av| 黄色怎么调成土黄色| 99九九线精品视频在线观看视频| 99热这里只有精品一区| 另类亚洲欧美激情| 美女视频免费永久观看网站| 高清视频免费观看一区二区| 色婷婷av一区二区三区视频| 亚洲国产色片| 国产精品.久久久| 欧美少妇被猛烈插入视频| 一级av片app| 在线天堂最新版资源| 我要看日韩黄色一级片| 十八禁高潮呻吟视频 | 丰满乱子伦码专区| 国产成人aa在线观看| 少妇的逼水好多| 91aial.com中文字幕在线观看| 18禁在线无遮挡免费观看视频| 久久99热6这里只有精品| www.av在线官网国产| 久久久久国产网址| 国产 一区精品| 人妻人人澡人人爽人人| 欧美一级a爱片免费观看看| 99热国产这里只有精品6| 久久午夜综合久久蜜桃| 妹子高潮喷水视频| 国产精品一区二区在线观看99| 亚洲一级一片aⅴ在线观看| 免费在线观看成人毛片| 国产精品一区二区三区四区免费观看| 欧美国产精品一级二级三级 | 国产成人精品无人区| 91久久精品国产一区二区三区| 新久久久久国产一级毛片| 国产精品成人在线| 男人和女人高潮做爰伦理| www.av在线官网国产| 又粗又硬又长又爽又黄的视频| 亚洲国产欧美日韩在线播放 | 麻豆乱淫一区二区| 国产有黄有色有爽视频| 国产欧美另类精品又又久久亚洲欧美| 91精品国产九色| av.在线天堂| 91精品国产国语对白视频| 国产成人一区二区在线| 欧美精品国产亚洲| 国产91av在线免费观看| 亚洲av中文av极速乱| 18禁在线无遮挡免费观看视频| 色5月婷婷丁香| 国内揄拍国产精品人妻在线| 久久久久精品性色| 亚洲在久久综合| 在线观看三级黄色| 国产精品国产三级国产专区5o| 久久久国产精品麻豆| 日韩强制内射视频| 在线观看一区二区三区激情| 亚洲精品日韩在线中文字幕| 一本久久精品| 免费大片黄手机在线观看| xxx大片免费视频| 国产91av在线免费观看| 亚洲第一区二区三区不卡| 国产av国产精品国产| a级片在线免费高清观看视频| 国产男人的电影天堂91| 女的被弄到高潮叫床怎么办| 亚洲精品乱码久久久久久按摩| 欧美 亚洲 国产 日韩一| h日本视频在线播放| 精品少妇久久久久久888优播| 亚洲精华国产精华液的使用体验| 亚洲国产av新网站| 最近手机中文字幕大全| 欧美老熟妇乱子伦牲交| 丰满人妻一区二区三区视频av| 久久久久久久大尺度免费视频| 国产黄频视频在线观看| 日韩 亚洲 欧美在线| 一级毛片黄色毛片免费观看视频| av卡一久久| 亚洲久久久国产精品| 狂野欧美激情性xxxx在线观看| 在线 av 中文字幕| 黄色欧美视频在线观看| 日韩人妻高清精品专区| 老司机影院成人| 精品酒店卫生间| 人人妻人人添人人爽欧美一区卜| 精品视频人人做人人爽| 七月丁香在线播放| 亚洲精品国产av蜜桃| 欧美精品一区二区免费开放| 欧美xxxx性猛交bbbb| 成人毛片60女人毛片免费| 五月天丁香电影| 亚洲电影在线观看av| 丰满少妇做爰视频| 欧美激情极品国产一区二区三区 | 91精品一卡2卡3卡4卡| 国产国拍精品亚洲av在线观看| av黄色大香蕉| 少妇人妻久久综合中文| 久久久久久人妻| 欧美精品一区二区免费开放| 色视频www国产| 一级,二级,三级黄色视频| 我要看黄色一级片免费的| 欧美最新免费一区二区三区| 黄色毛片三级朝国网站 | 国产亚洲最大av| 高清毛片免费看| 五月伊人婷婷丁香| 美女视频免费永久观看网站| 午夜av观看不卡| 五月天丁香电影| 日本黄大片高清| 在现免费观看毛片| 亚洲欧美一区二区三区黑人 | av黄色大香蕉| 国产一区二区在线观看av| 三上悠亚av全集在线观看 | 日本黄色片子视频| 亚洲精品,欧美精品| 麻豆精品久久久久久蜜桃| 成人毛片60女人毛片免费| 亚洲成人手机| 亚洲av二区三区四区| av免费观看日本| 亚洲精品第二区| 日韩一本色道免费dvd| 国产精品久久久久久精品电影小说| 日本免费在线观看一区| 插阴视频在线观看视频| 欧美最新免费一区二区三区| 国产成人a∨麻豆精品| 久久久久国产网址| 啦啦啦中文免费视频观看日本| 欧美人与善性xxx| 欧美日韩视频高清一区二区三区二| 男女无遮挡免费网站观看| 日韩中文字幕视频在线看片| 一级爰片在线观看| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 亚洲怡红院男人天堂| 伊人久久国产一区二区| 桃花免费在线播放| 免费黄频网站在线观看国产| 日日撸夜夜添| 亚洲精品一区蜜桃| 桃花免费在线播放| 日本91视频免费播放| 黑人高潮一二区| 日韩电影二区| 欧美少妇被猛烈插入视频| 国产精品人妻久久久久久| 偷拍熟女少妇极品色| 伊人久久精品亚洲午夜| 日日爽夜夜爽网站| 亚洲av二区三区四区| 能在线免费看毛片的网站| 亚洲电影在线观看av| 久久久久网色| 男女啪啪激烈高潮av片| 三级国产精品片| 亚洲熟女精品中文字幕| 中文字幕制服av| www.色视频.com| 国产日韩欧美在线精品| 有码 亚洲区| 97超碰精品成人国产| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 嫩草影院入口| 中文字幕人妻丝袜制服| 成年人午夜在线观看视频| 日韩一本色道免费dvd| 国产乱来视频区| 亚洲天堂av无毛| 亚洲av在线观看美女高潮| 日本vs欧美在线观看视频 | 国产精品人妻久久久影院| 日韩视频在线欧美| a级毛片免费高清观看在线播放| 久久av网站| 亚洲欧美日韩卡通动漫| 亚洲中文av在线| 日韩免费高清中文字幕av| 亚洲欧美精品自产自拍| av女优亚洲男人天堂| 日日啪夜夜撸| 嘟嘟电影网在线观看| 最黄视频免费看| 涩涩av久久男人的天堂| av免费在线看不卡| 国产免费又黄又爽又色| 成人影院久久| 国产成人精品福利久久| 丰满少妇做爰视频| 午夜福利在线观看免费完整高清在| 一本久久精品| 成年美女黄网站色视频大全免费 | 在线观看www视频免费| 久久精品久久久久久噜噜老黄| 新久久久久国产一级毛片| 国产男人的电影天堂91| 男人添女人高潮全过程视频| 丰满乱子伦码专区| 人妻一区二区av| 久久 成人 亚洲| 三级经典国产精品| 亚洲国产最新在线播放| 亚洲欧美精品自产自拍| 狂野欧美激情性bbbbbb| 91aial.com中文字幕在线观看| 久久久久久久国产电影| 亚洲av成人精品一二三区| 日本91视频免费播放| 国产精品欧美亚洲77777| 高清视频免费观看一区二区| 亚洲色图综合在线观看| 你懂的网址亚洲精品在线观看| 好男人视频免费观看在线| 人体艺术视频欧美日本| 熟女人妻精品中文字幕| 午夜av观看不卡| 乱人伦中国视频| 日本wwww免费看| 一个人免费看片子| 多毛熟女@视频| 精品人妻熟女毛片av久久网站| videossex国产| 伊人亚洲综合成人网| 国产精品久久久久久精品古装| 亚洲国产av新网站| 男女无遮挡免费网站观看| 9色porny在线观看| 亚洲av二区三区四区| av不卡在线播放| 国产精品国产三级专区第一集| 极品教师在线视频| 丰满人妻一区二区三区视频av| 另类亚洲欧美激情| 亚洲欧美日韩卡通动漫| 亚洲电影在线观看av|