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

    基于地震活動性和震源機制解研究渭河-運城盆地主要斷裂帶的特征及應力場分布

    2022-06-10 01:56:20余占洋沈旭章鄭文俊劉旭宙
    地震地質(zhì) 2022年2期
    關(guān)鍵詞:渭河應力場震源

    余占洋 沈旭章* 梁 浩 鄭文俊 劉旭宙

    1)中山大學地球科學與工程學院,廣東省地球動力作用與地質(zhì)災害重點實驗室,廣州 510275

    2)南方海洋科學與工程廣東省實驗室(珠海),珠海 519082

    3)甘肅省地震局,蘭州 730000

    0 引言

    本文利用近年來積累的大量震相報告資料,使用雙差定位法對渭河-運城盆地的地震事件進行精定位; 然后根據(jù)P波初動資料,使用格點嘗試法反演該區(qū)域內(nèi)部ML≥2地震的震源機制解; 并以反演的震源機制解結(jié)果為基礎(chǔ),應用MSATSI方法反演研究區(qū)域內(nèi)的應力場。通過截取不同的垂直斷層剖面,結(jié)合該地區(qū)已有的活動斷裂、構(gòu)造變形等研究結(jié)果,分析區(qū)域內(nèi)部的斷層結(jié)構(gòu)、區(qū)域應力場特征及地震發(fā)生的動力機制,探討區(qū)域現(xiàn)今構(gòu)造變形的模式,進而為渭河-運城盆地防震減災工作提供堅實基礎(chǔ)和重要參考。

    圖 1 渭河-運城盆地主要斷裂分布圖(圖中斷裂根據(jù)國家地震局“鄂爾多斯周緣活動斷裂系”課題組(1988)修改)Fig. 1 Distribution of main faults in Weihe-Yuncheng Basin(The fault in the figure was modified according to the Research Group of Active Fault System around Ordos,1988,issued by China Earthquake Administration).不同顏色的圓表示2009年1月1日—2018年6月1日該區(qū)域ML≥3的不同震級和深度的地震事件。F1秦嶺北緣斷裂; F2渭河斷裂; F3口鎮(zhèn)-關(guān)山斷裂; F4華山山前斷裂; F5韓城斷裂; F6雙泉-臨猗斷裂; F7中條山山前斷裂; F8羅云山山前斷裂; F9霍山山前斷裂; F10太谷斷裂; F11交城斷裂

    1 數(shù)據(jù)和方法

    1.1 觀測資料

    本研究選取了中國地震臺網(wǎng)中心“全國地震編目系統(tǒng)”記錄的2009年1月1日—2018年6月1日在渭河-運城盆地(圖 1)發(fā)生的14381個地震事件的數(shù)據(jù)進行重定位。在進行重定位前,需要對震相數(shù)據(jù)進行篩選,剔除震源參數(shù)及走時信息記錄不完整、震相記錄少于4個和記錄的臺站數(shù)量少于2個的地震事件,從而確保重定位的精確性。經(jīng)篩選,剩余地震事件11856個,共計443074個震相數(shù)據(jù),其中P波到時數(shù)據(jù)103874條,S波到時數(shù)據(jù)339200條。

    1.2 地震精定位方法

    本研究使用雙差定位法對區(qū)域內(nèi)的地震事件進行重定位。雙差定位法是一種相對地震定位方法,是在主事件地震定位法的基礎(chǔ)上發(fā)展起來的,將地震叢集中的地震事件兩兩組對,根據(jù)走時差的觀測值與理論值的殘差(即雙差)確定震源參數(shù)(Waldhauseretal.,2000)。雙差定位法的一項基本要求是地震事件對之間的距離遠小于事件到臺站間的距離和波傳播路徑上速度不均勻體的線性尺度,這樣才能使2個地震到同一個臺站的走時差只由這2個地震事件之間的相對位置和地震波的波速決定,能夠有效減小由于對地殼結(jié)構(gòu)了解不夠精細而引起的誤差,使重定位的結(jié)果更加準確(朱艾斕等,2005; 陳九輝等,2009; 盛書中等,2016)。

    圖 2 各區(qū)域的地殼速度模型Fig. 2 Crustal velocity model for each region.

    重定位反演的基本方程為

    WGm=Wd

    (1)

    式中,W為對每個方程進行加權(quán)的對角矩陣,G為包含偏導數(shù)的M×4N階(M為雙差的觀測數(shù),N為地震數(shù))矩陣,m是長度為4N的矢量,含有待定的震源參數(shù)變化(其形式為[Δx,Δy,Δz,ΔT]T),d是包含雙差的數(shù)據(jù)矢量(Waldhauseretal.,2000)。

    1.3 震源機制解反演方法

    本次研究使用格點嘗試法反演震源機制解。Aki(1966)提出采用大量小地震的P波初動資料推測應力場方向的綜合震源機制解法。在此基礎(chǔ)上,許忠淮等(1983)進一步提出了格點嘗試法。其原理為: 按震源的雙力偶點源模式,由震源坐標架在地平坐標架中的3個角度(Az,ih,R)來確定震源機制解。格點嘗試法以一定步長(如5°×5°×5°等)改變Az、ih和R的角度,觀察哪組(Az,ih,R)所對應的2個理論正交節(jié)面的矛盾符號比最小。

    定義矛盾符號比ψ=矛盾符號數(shù)/符號總數(shù),介于極小值ψmin和ψmin+5%之間所有可能的P、B、T軸的分布范圍稱為離散區(qū),對離散區(qū)內(nèi)所有可能的解求平均,計算得出的平均解即為地震事件的震源機制解結(jié)果。

    1.4 應力場反演方法

    本研究使用MSATSI程序包進行應力場反演。Michael(1984,1987)提出了自助線性應力反演(LSIB)方法,這是一種較為常用的利用震源機制解進行應力場反演的方法。Hardebeck等(2006)進一步發(fā)展了LSIB方法,并引入阻尼最小二乘法,得到了區(qū)域尺度的應力張量阻尼反演方法(SATSI)。此后,Martinez-Garzon等(2014)改進了SATSI算法,提出了MSATSI算法。

    在應力場反演中,首先將研究區(qū)域劃分成1°×1°的網(wǎng)格。為使相鄰網(wǎng)格的應力張量差異最小化,我們使用阻尼最小二乘法,該方法可有效減小人為劃分網(wǎng)格而產(chǎn)生的應力場的異常偏轉(zhuǎn),能更好地凸顯整體應力偏轉(zhuǎn)特征(崔華偉等,2019)。根據(jù)MSATSI軟件包提供的方法設(shè)置了一系列阻尼系數(shù),得到了數(shù)據(jù)擬合殘差與應力場反演模型長度的折中曲線,如圖 3 所示。折中曲線的拐點即為最佳阻尼系數(shù),本研究的最優(yōu)阻尼系數(shù)e=1.2。反演時每個網(wǎng)格中至少有2個震源機制解參與反演。反演應力場參數(shù)的置信水平默認范圍為68%~95%,本研究將置信水平設(shè)置為95%。重采樣次數(shù)的默認范圍是1000~5000次(Martinez-Garzonetal.,2014),本次研究對所有數(shù)據(jù)進行2000次抽樣以確定應力場參數(shù)的不確定度。

    圖 3 數(shù)據(jù)擬合殘差與應力場反演模型長度的折中曲線Fig. 3 Trade-off curve between model length and data fitting misfit in the stress field inversion model.藍色十字為最優(yōu)阻尼系數(shù),空心圓圈表示反演中采用的阻尼參數(shù)所得出的結(jié)果

    2 研究結(jié)果

    2.1 地震事件重定位結(jié)果

    依照前文設(shè)置的參數(shù)及地殼速度模型,應用雙差定位法對研究區(qū)域內(nèi)11856個地震事件進行了重定位。圖 4 顯示了定位前后的地震事件以及定位過程中使用的52個臺站的分布情況。圖4a 為定位前地震事件分布圖,圖4b 為定位后地震事件分布圖。結(jié)果顯示,定位后地震事件分布的離散程度降低,多數(shù)地震集中分布在太谷斷裂(F10)、交城斷裂(F11)、霍山山前斷裂(F9)和羅云山山前斷裂(F8)等斷裂帶附近; 另外,山西地塹帶內(nèi)部及周邊地震活動頻繁,但震級普遍較小。

    重定位后,所有地震的走時均方根殘差由0.618s降為0.255s。對比重定位前、后的震源深度直方圖(圖 5)可以發(fā)現(xiàn),重定位前震源深度集中分布在5~10km,重定位后震源深度的分布發(fā)生了明顯改變,地震集中分布在10~25km范圍內(nèi),總體上集中分布在20km以內(nèi),小部分地震發(fā)生在25~35km深度處。從圖 4 中可以看出,盆地內(nèi)部的震源深度相對較淺,以5~15km為主,向盆地兩端震源深度有加深的趨勢,最深約達30km,這與前人的研究結(jié)果也較為一致(蔡妍等,2014)。

    圖 4 研究區(qū)域內(nèi)重定位前(a)、后(b)地震事件的分布以及定位所用臺站的分布圖Fig. 4 The distribution of seismic events and stations before(a)and after(b)relocation in the research area.

    圖 5 研究區(qū)域內(nèi)重定位前(a)、后(b)地震事件的震源深度直方圖Fig. 5 Distribution of focal depth before(a)and(b)after relocation.

    2.2 震源機制結(jié)果

    本研究使用許忠淮等(1983)提出的利用P波初動資料反演震源機制解的格點嘗試法,使用從中國地震臺網(wǎng)中心獲取的震相報告的P波初動資料,計算了渭河-運城盆地內(nèi)346個ML≥2地震的震源機制解。

    根據(jù)Zoback(1992)給出的世界應力圖劃分準則(表1)將震源機制解劃分正斷型(NF)、正走滑型(NS)、逆沖型(TF)、逆走滑型(TS)、走滑型(SS)、過渡型(U)等6種類型的震源機制解。結(jié)果表明: 走滑斷層類型地震有69個,占19.9%; 正斷層及正走滑斷層類型有145個,占41.9%; 逆斷層和逆走滑斷層類型有36個,占10.4%; 過渡型有96個,占27.8%。

    表 1 震源機制解的類型劃分表Table1 Categories of tectonic stress regime for focal mechanism

    圖 6 研究區(qū)域內(nèi)的震源機制解分布圖Fig. 6 The distribution of focal mechanism solutions.

    從圖 6 中可以看出,渭河-運城盆地內(nèi)的正斷型和走滑型地震事件分布廣泛,占全部地震事件的60%以上,且大多集中在太谷斷裂(F10)、交城斷裂(F11)、霍山山前斷裂(F9)和羅云山山前斷裂(F8)附近,表明渭河-運城盆地的控盆斷裂及內(nèi)部斷裂的運動方式以拉張運動為主,兼有一定的剪切運動。

    圖 7 利用不同方法得到的震源機制解Fig. 7 Focal mechanism solutions obtained by different methods.a 格點嘗試法的結(jié)果; b Hash方法的結(jié)果

    表 2 利用不同方法得到的震源機制解Table2 Focal mechanism solutions of earthquakes obtained by different methods

    圖 8 利用不同方法得到的震源機制解Fig. 8 Focal mechanism solutions obtained by different methods.a—c為本文結(jié)果,分別對應高陵地震、河津地震、原平地震; d—f為前人的研究結(jié)果,分別對應高陵地震、河津地震、原平地震

    表 3 通過不同方法得到的MS>4.0地震的震源機制解Table3 Focal mechanism solutions of earthquakes above MS4.0 obtained by different methods

    2.3 應力場反演結(jié)果

    圖 9 研究區(qū)的應力場反演結(jié)果Fig. 9 Results of stress field inversion in the research area.繪出了每個網(wǎng)格點的最大水平應力方向,紅色表示正斷層作用體系,綠色表示走滑斷層體系,網(wǎng)格的顏色表示R值的大小,顏色越深則R值越大

    3 討論

    渭河-運城盆地的地震重定位結(jié)果顯示,地震活動多沿著現(xiàn)今地表斷裂的走向發(fā)生,其震源深度的空間分布指示斷裂在盆地深部的幾何結(jié)構(gòu)和空間位置。為厘定斷裂在盆地深部的展布特征,本研究以20km寬的廊帶為約束,分別對A、B、C區(qū)域做與現(xiàn)今斷裂交切的深度剖面進行分析,共做8條剖面。

    圖 10 A區(qū)域的剖面位置及地震事件剖面圖Fig. 10 Sectional location map and seismic event profiles of region A.a 剖面分布圖; b—d AA′、BB′和CC′剖面

    A區(qū)域(渭河斷陷帶)內(nèi)發(fā)育有秦嶺北緣斷裂(F1)、渭河斷裂(F2)和口鎮(zhèn)-關(guān)山斷裂(F3)3條主要斷裂,在此區(qū)域截取了AA′、BB′和CC′ 3條跨越秦嶺山系、渭河地塹和鄂爾多斯高原且與主要斷裂交切的NNE向剖面,如圖 10 所示。剖面圖中顯示了MS≥3.5地震的震源機制解,其他區(qū)域剖面中所顯示的也是MS≥3.5地震的震源機制解,下文不再說明。AA′剖面跨越秦嶺北緣斷裂(F1)西段和渭河斷裂(F2)。F1西段全長110km,走向近EW—NWW,斷裂帶寬數(shù)米至數(shù)十米不等,地表露頭和探槽揭露為高角度N傾的正斷層(王明秋,2011; 白相東,2018); F2全長50km,西段走向近EW,東段走向NEE,沿斷裂走向發(fā)育明顯的斷層陡坎,鉆孔揭露為高角度正斷層(田勤虎等,2014)。在深部,重定位地震事件多在F1左盤附近,呈線性延伸至地下10km處,斷裂傾角近直立,與地表露頭所示的斷面產(chǎn)狀一致。地表露頭揭露F2產(chǎn)狀近直立,但因地震數(shù)據(jù)較少,其在盆地深部產(chǎn)狀仍不明確,推測F2為與F1近平行的N傾斷裂。BB′剖面穿過秦嶺北緣斷裂(F1)東段、渭河斷裂和口鎮(zhèn)-關(guān)山斷裂。F1東段整體走向近EW,全長51km,地表露頭顯示斷裂下盤為花崗巖,上盤為黃土,與黃土接觸的斷面上發(fā)育擦痕,指示F1東段為具走滑性質(zhì)的左旋正斷層; F2的性質(zhì)變化不大; F3整體走向近EW,總長約100km,傾向S,傾角為40°~70°,斷裂沿線發(fā)育的階地位錯、黃土陡坎與地裂縫等活動標志多集中分布于斷裂西段(米豐收等,1993)。深部的地震事件基本分布在深度10km以內(nèi),同時可以發(fā)現(xiàn)口鎮(zhèn)-關(guān)山斷裂向深部延伸傾角變緩且具有一定的走滑性質(zhì)。CC′剖面跨過渭河斷裂(F2)和口鎮(zhèn)-關(guān)山斷裂(F3),震源機制解揭示渭河斷裂的斷層性質(zhì)為正斷兼走滑??傮w而言,渭河地塹剖面的地震活動性相對較弱,區(qū)域內(nèi)部的斷層以正斷為主,兼有一定走滑性質(zhì)。

    圖 11 B區(qū)域的剖面位置及地震事件剖面圖Fig. 11 Sectional location map and seismic event profiles of region B.a 剖面分布圖; b、c DD′和EE′剖面

    B區(qū)域內(nèi)發(fā)育有華山山前斷裂(F4)、韓城斷裂(F5)、雙泉-臨猗斷裂(F6)、中條山山前斷裂(F7)和羅云山山前斷裂(F8)5條主要斷裂,在此區(qū)域截取了DD′和EE′ 2條剖面,如圖 11 所示。DD′剖面穿過雙泉-臨猗斷裂、中條山山前斷裂和韓城斷裂,其中,中條山山前斷裂和韓城斷裂是圍限侯馬-運城盆地的主要斷裂,雙泉-臨猗斷裂則是侯馬-運城盆地內(nèi)部峨嵋臺地隆起的南緣斷裂。在地表,中條山山前斷裂位于運城盆地與中條山之間,全長約130km,走向NEE—NE,斷面傾向NNW—NW,傾角為60°~80°,為高角度正斷層(司蘇沛等,2014; 郭春杉等,2019); 韓城斷裂整體走向NE,長約120km,為高角度正斷層(扈桂讓等,2017); 雙泉-臨猗斷裂走向NE,斷面總體S傾但出露不全,斷裂下盤為黃土臺地,上盤為運城盆地沖積平原,高差明顯,指示斷層為正斷層。在深部,地震事件主要沿中條山山前斷裂和韓城斷裂分布,中條山山前斷裂向盆地的分支斷裂亦有揭露,與錯斷中條山山前洪積扇陡坎的斷裂有較好的對應(田建梅等,2013)。中條山山前斷裂和韓城斷裂均顯示出高角度的正斷性質(zhì),與地表產(chǎn)狀對應,而雙泉-臨猗斷裂地震事件較少,活動性較弱。EE′剖面主要跨越羅云山山前斷裂北端。羅云山山前斷裂是臨汾盆地西側(cè)的控盆斷裂,全長145km,總體走向NNE,其南段折向NWW后與韓城斷裂相連,為高角度正斷層(孫昌斌等,2011)。羅云山山前斷裂的地震事件主要分布在深度10~20km之間,震源深度較大,大致勾勒出高角度的斷裂面,與地表產(chǎn)狀相對應。與A區(qū)域相比,B區(qū)域的地震活動性明顯增強。

    C區(qū)域(山西地塹系南部)內(nèi)發(fā)育了霍山山前斷裂(F9)、太谷斷裂(F10)和交城斷裂(F11)3條主要斷裂,在此區(qū)域截取了FF′、GG′和HH′ 3條剖面,如圖 12 所示。斷裂調(diào)查顯示,霍山山前斷裂為臨汾盆地東北側(cè)的邊界斷裂,總體走向NNE,運動性質(zhì)以正斷為主,兼具右旋走滑性質(zhì)(畢麗思等,2017); 交城斷裂是晉中盆地西界的控盆斷裂,走向NE,長約125km,為SE傾正斷層,傾角為50°~70°(孫昌斌等,2012); 太谷斷裂是晉中盆地東邊界的控盆斷裂,走向與交城斷裂大致平行,全長約100km,為NW傾正斷層。FF′剖面穿過霍山山前斷裂南端,地震活動性更強,震源深度分布廣,有少量地震事件的震源深度已超過30km,斷裂帶呈現(xiàn)出高角度正斷層的特點,與地表產(chǎn)狀相匹配。GG′剖面穿過太谷斷裂和交城斷裂,可以看出地震事件基本分布在2條斷裂之間,指示盆地內(nèi)部的地震活動較為頻繁。HH′剖面跨越太谷斷裂和交城斷裂的北端,地震事件的震源深度主要集中在10~30km之間,且2條斷裂的傾角很大,接近直立。與A、B區(qū)域相比,C區(qū)域的地震活動性最強。

    圖 12 C區(qū)域的剖面位置及地震事件剖面圖Fig. 12 Sectional location map and seismic event profiles of region C.a 剖面分布圖; b—d FF′、GG′、HH′剖面

    4 結(jié)論

    本研究選取了近年來渭河-運城盆地的震相報告資料,采用雙差定位法對地震位置進行重定位,得到了8106個地震事件的重定位結(jié)果; 利用格點嘗試法得到了346個MS≥2地震的震源機制解,并基于震源機制解結(jié)果通過應力張量阻尼反演方法使用MSATSI軟件包反演獲得該區(qū)域的應力場,得出了以下主要結(jié)論:

    (1)渭河-運城盆地的地震事件重定位之后相較于定位前的分布更集中,大部分地震在斷裂帶附近密集分布; 根據(jù)截取的剖面可知,研究區(qū)域內(nèi)大部分斷裂的傾角較高,有的接近直立狀態(tài),如太谷斷裂和交城斷裂等。

    (2)研究區(qū)域內(nèi)的震源機制解類型以走滑型和正斷型為主,斷裂帶的構(gòu)造運動性質(zhì)以走滑和正斷為主。

    致謝審稿專家對本文提出了中肯建議; 國家地震科學數(shù)據(jù)共享中心提供了地震震相報告; 德國地球科學研究中心(GFZ)提供了反演程序MSATSI; Waldhauser共享了雙差定位程序; 萬永革老師提供了計算不同震源機制的最小空間旋轉(zhuǎn)角的程序; 文中圖件使用GMT(Wesseletal.,2019)繪制。在此一并表示感謝!

    猜你喜歡
    渭河應力場震源
    退耕還林工程對渭河洪峰流量的影響
    震源的高返利起步
    鋁合金多層多道窄間隙TIG焊接頭應力場研究
    焊接(2016年9期)2016-02-27 13:05:22
    渭河故道恢復濕地生態(tài)的引水方案創(chuàng)新
    可控震源地震在張掖盆地南緣逆沖斷裂構(gòu)造勘探中的應用
    同步可控震源地震采集技術(shù)新進展
    考慮斷裂破碎帶的丹江口庫區(qū)地應力場與水壓應力場耦合反演及地震預測
    1996—2009年渭河干流氮素污染特征
    基于位移相關(guān)法的重復壓裂裂縫尖端應力場研究
    斷塊油氣田(2014年5期)2014-03-11 15:33:49
    岸坡應力場及卸荷帶劃分量化指標研究
    成人三级做爰电影| 亚洲片人在线观看| 成人特级黄色片久久久久久久| 69av精品久久久久久| ponron亚洲| 免费看十八禁软件| 国产精品av久久久久免费| 岛国视频午夜一区免费看| 韩国精品一区二区三区| 亚洲自拍偷在线| avwww免费| 波多野结衣高清无吗| 亚洲国产中文字幕在线视频| 美女 人体艺术 gogo| 制服丝袜大香蕉在线| 免费在线观看完整版高清| 91老司机精品| 亚洲情色 制服丝袜| 一级毛片女人18水好多| 欧美日韩瑟瑟在线播放| 亚洲人成网站在线播放欧美日韩| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩一区二区三| 久久久久久大精品| 这个男人来自地球电影免费观看| 国产精品亚洲av一区麻豆| av天堂久久9| 黑人操中国人逼视频| 亚洲一区二区三区不卡视频| 欧美成人午夜精品| 少妇粗大呻吟视频| 国产在线观看jvid| 精品乱码久久久久久99久播| 18禁美女被吸乳视频| av福利片在线| 久久午夜综合久久蜜桃| 日韩国内少妇激情av| 日本黄色视频三级网站网址| 啦啦啦免费观看视频1| 国产黄a三级三级三级人| 色精品久久人妻99蜜桃| 精品国产一区二区久久| 久久热在线av| 国产成+人综合+亚洲专区| 国产精品综合久久久久久久免费 | 19禁男女啪啪无遮挡网站| 老鸭窝网址在线观看| 老司机午夜十八禁免费视频| 午夜免费鲁丝| 手机成人av网站| 老汉色∧v一级毛片| 99re在线观看精品视频| 国产成年人精品一区二区| 日韩欧美免费精品| 国产精华一区二区三区| 亚洲免费av在线视频| xxx96com| 欧美乱码精品一区二区三区| 国产精品99久久99久久久不卡| 亚洲全国av大片| 久久亚洲精品不卡| 久久久久久大精品| 久久久久九九精品影院| 国产麻豆69| 制服人妻中文乱码| 精品国产国语对白av| 国产亚洲av高清不卡| 自线自在国产av| 亚洲精品久久成人aⅴ小说| 琪琪午夜伦伦电影理论片6080| a级毛片在线看网站| 久久久精品国产亚洲av高清涩受| 在线播放国产精品三级| 久久热在线av| 国产麻豆成人av免费视频| 精品久久久久久久毛片微露脸| 久久香蕉国产精品| 18禁美女被吸乳视频| 黄色片一级片一级黄色片| 丰满人妻熟妇乱又伦精品不卡| 成熟少妇高潮喷水视频| 免费搜索国产男女视频| 国内精品久久久久久久电影| e午夜精品久久久久久久| 精品日产1卡2卡| 亚洲第一电影网av| 久久午夜综合久久蜜桃| 可以在线观看毛片的网站| 少妇熟女aⅴ在线视频| 欧美亚洲日本最大视频资源| 免费少妇av软件| 久久久精品欧美日韩精品| 好男人在线观看高清免费视频 | 成人手机av| 好男人在线观看高清免费视频 | 免费在线观看亚洲国产| 精品久久蜜臀av无| 亚洲情色 制服丝袜| 亚洲av成人一区二区三| 亚洲一区二区三区不卡视频| 国产精品免费视频内射| 精品国产亚洲在线| 亚洲少妇的诱惑av| avwww免费| 99国产精品免费福利视频| 欧美日韩中文字幕国产精品一区二区三区 | 久久婷婷成人综合色麻豆| 免费看a级黄色片| 国产欧美日韩一区二区三| 国产精品爽爽va在线观看网站 | 91麻豆精品激情在线观看国产| 亚洲熟妇中文字幕五十中出| 日韩欧美免费精品| 国产亚洲av嫩草精品影院| 两人在一起打扑克的视频| 91大片在线观看| 久久精品亚洲精品国产色婷小说| 成人国产一区最新在线观看| 亚洲av第一区精品v没综合| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区视频在线观看免费| 90打野战视频偷拍视频| 欧美色欧美亚洲另类二区 | 亚洲欧美一区二区三区黑人| 天堂影院成人在线观看| 国产一区二区三区视频了| 亚洲成av片中文字幕在线观看| 18美女黄网站色大片免费观看| 成人18禁在线播放| 老司机靠b影院| 亚洲久久久国产精品| 好男人在线观看高清免费视频 | 亚洲一区二区三区色噜噜| 变态另类丝袜制服| 午夜老司机福利片| 成人国产综合亚洲| 亚洲久久久国产精品| 国产极品粉嫩免费观看在线| 岛国在线观看网站| АⅤ资源中文在线天堂| 日本一区二区免费在线视频| 欧美乱色亚洲激情| av欧美777| 欧美成人一区二区免费高清观看 | 欧美日韩中文字幕国产精品一区二区三区 | 一级毛片高清免费大全| 午夜福利欧美成人| 在线观看日韩欧美| 欧美激情久久久久久爽电影 | 亚洲 欧美一区二区三区| 日韩精品青青久久久久久| 久久中文字幕一级| 天堂影院成人在线观看| 国产精品综合久久久久久久免费 | 脱女人内裤的视频| 啦啦啦免费观看视频1| 1024视频免费在线观看| 亚洲av成人一区二区三| 亚洲视频免费观看视频| bbb黄色大片| 又黄又粗又硬又大视频| 久久天堂一区二区三区四区| 欧美人与性动交α欧美精品济南到| 一区二区三区国产精品乱码| 日本 av在线| 午夜福利视频1000在线观看 | 免费观看人在逋| 多毛熟女@视频| 国产亚洲av高清不卡| 男人舔女人下体高潮全视频| 男女床上黄色一级片免费看| 久久久久久久午夜电影| 极品教师在线免费播放| www国产在线视频色| 日韩大码丰满熟妇| 97超级碰碰碰精品色视频在线观看| 成人永久免费在线观看视频| 妹子高潮喷水视频| 搡老岳熟女国产| 老汉色av国产亚洲站长工具| 国产高清有码在线观看视频 | 国产高清videossex| 欧美性长视频在线观看| 国产高清有码在线观看视频 | 久久久精品欧美日韩精品| 国产国语露脸激情在线看| 可以在线观看毛片的网站| 精品国产亚洲在线| 久久天堂一区二区三区四区| 国内毛片毛片毛片毛片毛片| 成人免费观看视频高清| 欧美成人免费av一区二区三区| 国产精品爽爽va在线观看网站 | 欧美一级a爱片免费观看看 | 欧美+亚洲+日韩+国产| 伦理电影免费视频| 欧美另类亚洲清纯唯美| 成年女人毛片免费观看观看9| 午夜老司机福利片| 久久久久国产一级毛片高清牌| 欧美色欧美亚洲另类二区 | 色哟哟哟哟哟哟| 国产一级毛片七仙女欲春2 | 亚洲国产高清在线一区二区三 | 男女下面进入的视频免费午夜 | 国产野战对白在线观看| xxx96com| 亚洲人成网站在线播放欧美日韩| 欧美日韩瑟瑟在线播放| √禁漫天堂资源中文www| 久久精品人人爽人人爽视色| 欧美久久黑人一区二区| 久久伊人香网站| 成人国语在线视频| 午夜福利在线观看吧| 午夜免费观看网址| 女人被躁到高潮嗷嗷叫费观| 亚洲色图av天堂| 欧美日本视频| 国产精品电影一区二区三区| 国产精品一区二区精品视频观看| 国产黄a三级三级三级人| 免费不卡黄色视频| 精品不卡国产一区二区三区| 成人亚洲精品一区在线观看| 丰满人妻熟妇乱又伦精品不卡| 999精品在线视频| 久久香蕉国产精品| 69av精品久久久久久| 男女床上黄色一级片免费看| 曰老女人黄片| 久久欧美精品欧美久久欧美| 久久精品aⅴ一区二区三区四区| 无遮挡黄片免费观看| 男女下面进入的视频免费午夜 | 性欧美人与动物交配| 亚洲九九香蕉| 亚洲aⅴ乱码一区二区在线播放 | 欧美成人一区二区免费高清观看 | 亚洲男人天堂网一区| 99香蕉大伊视频| 中文字幕精品免费在线观看视频| 99国产精品免费福利视频| av片东京热男人的天堂| 在线观看www视频免费| 国产一区二区三区在线臀色熟女| 国产精品一区二区在线不卡| 国产在线观看jvid| 亚洲人成电影免费在线| 亚洲精品一卡2卡三卡4卡5卡| 女生性感内裤真人,穿戴方法视频| 亚洲一区高清亚洲精品| 黄色 视频免费看| 亚洲性夜色夜夜综合| 亚洲avbb在线观看| 国产精品久久久久久亚洲av鲁大| a在线观看视频网站| 午夜日韩欧美国产| 99热只有精品国产| 婷婷丁香在线五月| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 精品午夜福利视频在线观看一区| 国产精品综合久久久久久久免费 | 窝窝影院91人妻| 巨乳人妻的诱惑在线观看| 亚洲成av人片免费观看| 欧美精品啪啪一区二区三区| 女生性感内裤真人,穿戴方法视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品在线美女| cao死你这个sao货| 亚洲国产日韩欧美精品在线观看 | 电影成人av| 丁香六月欧美| 国产精品久久久人人做人人爽| 18禁黄网站禁片午夜丰满| 欧美不卡视频在线免费观看 | 午夜福利一区二区在线看| 在线观看日韩欧美| 国产精品亚洲一级av第二区| 搡老妇女老女人老熟妇| 麻豆一二三区av精品| 一区福利在线观看| 在线av久久热| 男人舔女人的私密视频| 久久性视频一级片| 97人妻天天添夜夜摸| 深夜精品福利| svipshipincom国产片| 黑人巨大精品欧美一区二区蜜桃| 成年人黄色毛片网站| 久久精品亚洲精品国产色婷小说| 精品国产一区二区三区四区第35| 麻豆av在线久日| 法律面前人人平等表现在哪些方面| 99国产精品一区二区三区| 丁香欧美五月| 1024视频免费在线观看| 黄色a级毛片大全视频| 亚洲三区欧美一区| 亚洲熟妇熟女久久| 精品无人区乱码1区二区| 久久久久久久久中文| 999精品在线视频| 亚洲精品美女久久av网站| 男人操女人黄网站| 国产精品亚洲美女久久久| 亚洲欧美日韩无卡精品| 母亲3免费完整高清在线观看| 国产激情久久老熟女| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 高清毛片免费观看视频网站| 一级毛片精品| 97人妻精品一区二区三区麻豆 | 国产成人精品久久二区二区91| 91精品国产国语对白视频| 精品国产美女av久久久久小说| 久久草成人影院| 日韩精品青青久久久久久| 男人操女人黄网站| 999精品在线视频| 亚洲人成电影观看| 国产亚洲av高清不卡| 中文字幕另类日韩欧美亚洲嫩草| 亚洲午夜理论影院| 久久久国产精品麻豆| 欧美一级a爱片免费观看看 | 国产又色又爽无遮挡免费看| 国产单亲对白刺激| 久久人人精品亚洲av| av福利片在线| 12—13女人毛片做爰片一| 黄片大片在线免费观看| 亚洲成人免费电影在线观看| 久久久久久久久久久久大奶| 精品福利观看| 日本 欧美在线| 99国产精品一区二区三区| 成在线人永久免费视频| 亚洲国产精品合色在线| 伊人久久大香线蕉亚洲五| 国产精品乱码一区二三区的特点 | 日日干狠狠操夜夜爽| 啦啦啦韩国在线观看视频| 国产精品影院久久| 精品久久久久久久毛片微露脸| 国产精品一区二区在线不卡| 欧美日本亚洲视频在线播放| 国产高清有码在线观看视频 | 亚洲精品在线观看二区| 色精品久久人妻99蜜桃| 韩国av一区二区三区四区| 99精品欧美一区二区三区四区| 精品久久蜜臀av无| 在线观看免费视频日本深夜| 亚洲中文av在线| 久久午夜综合久久蜜桃| 国产xxxxx性猛交| 69精品国产乱码久久久| 午夜激情av网站| 香蕉丝袜av| 国产成人免费无遮挡视频| 婷婷六月久久综合丁香| 免费不卡黄色视频| 欧美一区二区精品小视频在线| 51午夜福利影视在线观看| 国产99久久九九免费精品| 在线观看免费视频网站a站| 亚洲人成电影免费在线| 国产亚洲欧美精品永久| 亚洲精品中文字幕在线视频| 亚洲精品中文字幕一二三四区| 成人欧美大片| 亚洲成人国产一区在线观看| 可以免费在线观看a视频的电影网站| 一边摸一边做爽爽视频免费| 久久久久久大精品| svipshipincom国产片| 成人av一区二区三区在线看| 国产亚洲欧美在线一区二区| 精品久久久久久成人av| 亚洲中文字幕一区二区三区有码在线看 | 精品福利观看| 久久久久久国产a免费观看| 欧美成人免费av一区二区三区| 999久久久精品免费观看国产| 中文字幕另类日韩欧美亚洲嫩草| xxx96com| 国产又爽黄色视频| 中文字幕人妻熟女乱码| 一边摸一边抽搐一进一出视频| 长腿黑丝高跟| 成年女人毛片免费观看观看9| 成人特级黄色片久久久久久久| 亚洲av电影在线进入| 男男h啪啪无遮挡| 久久精品人人爽人人爽视色| av电影中文网址| 欧美激情高清一区二区三区| 在线十欧美十亚洲十日本专区| 国产亚洲精品久久久久5区| 91九色精品人成在线观看| 成在线人永久免费视频| 在线十欧美十亚洲十日本专区| 亚洲av熟女| 成人国产综合亚洲| 又大又爽又粗| 老司机靠b影院| 久久香蕉激情| 亚洲伊人色综图| av片东京热男人的天堂| 欧美在线黄色| 黄色丝袜av网址大全| 久久中文看片网| 午夜免费鲁丝| 久久久国产精品麻豆| 亚洲国产日韩欧美精品在线观看 | 久久久精品国产亚洲av高清涩受| 99国产极品粉嫩在线观看| 日本撒尿小便嘘嘘汇集6| 最好的美女福利视频网| 亚洲欧美精品综合久久99| 精品一品国产午夜福利视频| 午夜久久久久精精品| 欧美黄色片欧美黄色片| 国产精品野战在线观看| 香蕉丝袜av| 又黄又爽又免费观看的视频| 丝袜人妻中文字幕| 国产高清视频在线播放一区| 久久久精品国产亚洲av高清涩受| 国产精品久久久av美女十八| 精品久久久久久成人av| 精品乱码久久久久久99久播| 久久 成人 亚洲| 亚洲人成77777在线视频| 欧美乱码精品一区二区三区| 欧美成人免费av一区二区三区| 美女午夜性视频免费| 国产精品 欧美亚洲| 久久中文看片网| 欧美最黄视频在线播放免费| 嫩草影院精品99| 熟妇人妻久久中文字幕3abv| 国产精品 国内视频| 亚洲三区欧美一区| 国产麻豆成人av免费视频| 欧美日本亚洲视频在线播放| 国产精品久久电影中文字幕| 深夜精品福利| 色综合婷婷激情| 一夜夜www| av免费在线观看网站| 一本久久中文字幕| 18禁观看日本| 大码成人一级视频| 亚洲视频免费观看视频| 免费一级毛片在线播放高清视频 | 国产精品99久久99久久久不卡| 嫁个100分男人电影在线观看| 麻豆久久精品国产亚洲av| 久久久精品国产亚洲av高清涩受| 巨乳人妻的诱惑在线观看| 国产精品永久免费网站| 亚洲熟妇熟女久久| 无限看片的www在线观看| 99国产精品一区二区蜜桃av| 日韩中文字幕欧美一区二区| 午夜老司机福利片| 巨乳人妻的诱惑在线观看| 天堂影院成人在线观看| 19禁男女啪啪无遮挡网站| 搡老熟女国产l中国老女人| 高清黄色对白视频在线免费看| 国产视频一区二区在线看| 中文字幕人妻熟女乱码| 成人三级黄色视频| 国产人伦9x9x在线观看| 美女高潮到喷水免费观看| 成人永久免费在线观看视频| 一边摸一边抽搐一进一小说| 999久久久国产精品视频| 桃色一区二区三区在线观看| 午夜福利一区二区在线看| 中文字幕久久专区| 免费女性裸体啪啪无遮挡网站| 97人妻天天添夜夜摸| tocl精华| 在线视频色国产色| 九色国产91popny在线| 亚洲精品国产精品久久久不卡| 欧美乱码精品一区二区三区| 1024视频免费在线观看| 999久久久国产精品视频| 老司机午夜福利在线观看视频| 少妇被粗大的猛进出69影院| 麻豆一二三区av精品| 日韩欧美免费精品| 黑人巨大精品欧美一区二区mp4| 99在线视频只有这里精品首页| 亚洲伊人色综图| 欧美中文综合在线视频| www.999成人在线观看| 好看av亚洲va欧美ⅴa在| 母亲3免费完整高清在线观看| 国产成人欧美| 老司机午夜福利在线观看视频| 久久精品国产清高在天天线| 黄色毛片三级朝国网站| 香蕉丝袜av| av免费在线观看网站| 好男人在线观看高清免费视频 | 国产精品免费一区二区三区在线| 久久国产亚洲av麻豆专区| 国产亚洲精品久久久久5区| 曰老女人黄片| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲成人精品中文字幕电影| 国产高清videossex| 欧美日本中文国产一区发布| 韩国av一区二区三区四区| 亚洲黑人精品在线| 国产午夜福利久久久久久| 非洲黑人性xxxx精品又粗又长| 国产又爽黄色视频| 国产精品一区二区免费欧美| 亚洲无线在线观看| 久久久久国产精品人妻aⅴ院| 欧美成人一区二区免费高清观看 | 狠狠狠狠99中文字幕| 亚洲欧美激情在线| 亚洲人成电影免费在线| 日本欧美视频一区| 久久久久久久久免费视频了| 国产精品自产拍在线观看55亚洲| 一边摸一边抽搐一进一出视频| 欧美乱妇无乱码| av中文乱码字幕在线| 国产精品 欧美亚洲| 欧美日本中文国产一区发布| 欧美绝顶高潮抽搐喷水| 久久精品成人免费网站| 在线播放国产精品三级| 黄色视频不卡| 女性被躁到高潮视频| 丝袜人妻中文字幕| 国产精品久久久久久亚洲av鲁大| 欧美另类亚洲清纯唯美| 日韩三级视频一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 一二三四在线观看免费中文在| 69精品国产乱码久久久| 美女大奶头视频| 国产精品一区二区免费欧美| 丝袜美腿诱惑在线| 久久国产精品影院| 操出白浆在线播放| 日韩精品免费视频一区二区三区| 亚洲色图综合在线观看| 大型av网站在线播放| 国产一区二区三区视频了| 国产区一区二久久| 欧美黄色片欧美黄色片| 久久精品亚洲精品国产色婷小说| 桃红色精品国产亚洲av| 日本五十路高清| 欧美中文综合在线视频| 国产亚洲欧美在线一区二区| 欧美日韩精品网址| 久99久视频精品免费| 国产xxxxx性猛交| 中文字幕久久专区| 亚洲一码二码三码区别大吗| 成人av一区二区三区在线看| 人妻久久中文字幕网| 99国产精品免费福利视频| 无限看片的www在线观看| 免费看十八禁软件| 亚洲av成人不卡在线观看播放网| 成人18禁在线播放| 国产成人精品在线电影| 亚洲一区中文字幕在线| 一级片免费观看大全| 亚洲av五月六月丁香网| 国产欧美日韩综合在线一区二区| 久久亚洲精品不卡| 真人做人爱边吃奶动态| 夜夜爽天天搞| 久久人人97超碰香蕉20202| 欧美国产日韩亚洲一区| 国产精品av久久久久免费| 亚洲精品中文字幕在线视频| 国产激情久久老熟女| 天天添夜夜摸| 亚洲电影在线观看av| 可以免费在线观看a视频的电影网站| 日韩三级视频一区二区三区| 在线观看66精品国产| 久久久国产成人精品二区| 长腿黑丝高跟| 亚洲精品一卡2卡三卡4卡5卡| 亚洲免费av在线视频| 亚洲精品美女久久av网站| 成年版毛片免费区| 日韩一卡2卡3卡4卡2021年| 亚洲,欧美精品.| 亚洲一区中文字幕在线| 精品无人区乱码1区二区| 久久精品国产亚洲av香蕉五月| 9191精品国产免费久久| 国产精华一区二区三区| 男女下面进入的视频免费午夜 | 九色亚洲精品在线播放|