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

    青海瑪多MS7.4地震前后秒采樣地電場(chǎng)動(dòng)態(tài)變化

    2022-02-23 12:44:00席繼樓趙家騮高尚華王曉蕾李國(guó)佑張興王濤
    地球物理學(xué)報(bào) 2022年2期
    關(guān)鍵詞:大武瑪多階躍

    席繼樓,趙家騮,高尚華,王曉蕾,李國(guó)佑,張興,王濤

    1 中國(guó)地震局地震預(yù)測(cè)研究所,北京 100036 2 中國(guó)地震局地震預(yù)測(cè)重點(diǎn)實(shí)驗(yàn)室,北京 100036 3 青海省地震局,西寧 810001

    0 引言

    據(jù)中國(guó)地震臺(tái)網(wǎng)中心測(cè)定,北京時(shí)間2021年5月22日2時(shí)4分29秒,在青海省果洛藏族自治州瑪多縣黃河鄉(xiāng)(E98.34°,N34.59°)發(fā)生了7.4級(jí)地震,震源深度為17 km(王龍等,2021).瑪多MS7.4地震是我國(guó)大陸地區(qū)繼2017年九寨溝MS7.0地震之后,也是2008年汶川MS8.0地震以來,震級(jí)最高的一次地震.

    地電場(chǎng)是重要地球物理場(chǎng)之一,以地球表層的天然電場(chǎng)及其隨時(shí)間變化為主要觀測(cè)對(duì)象,以地下介質(zhì)場(chǎng)源產(chǎn)生的地電場(chǎng)及其動(dòng)態(tài)變化過程為重要研究目標(biāo).多年來開展的理論分析、數(shù)值模擬和試驗(yàn)研究等研究結(jié)果表明,在強(qiáng)地震活動(dòng)過程中,地電場(chǎng)的動(dòng)態(tài)變化,主要與地應(yīng)力和地應(yīng)變作用下,地下介質(zhì)的構(gòu)造斷層活動(dòng)、裂隙發(fā)育和電解質(zhì)運(yùn)移和滲流等物理過程,以及地下介質(zhì)電磁學(xué)屬性及其動(dòng)態(tài)變化等相關(guān)聯(lián)(孫正江和王華俊,1984;Huang and Ikeya,1999;Huang, 2002;錢復(fù)業(yè)和趙玉林,2005;蘇巍等,2006;黃清華和林玉峰,2010;徐志鋒和吳小平,2010;湯吉等,2010;Ren et al.,2010;張丹等,2013;譚大城等,2013,2014;Huang et al.,2015;Ren et al.,2016;王軍等,2016;席繼樓等,2018b;王宇等,2020).

    21世紀(jì)以來,基于最新建設(shè)完成的數(shù)字化地電場(chǎng)觀測(cè)臺(tái)網(wǎng),我國(guó)地震工作者對(duì)發(fā)生在中國(guó)大陸及周邊地區(qū)的印尼蘇門答臘MS8.7、尼泊爾博克拉MS8.1、汶川MS8.0、玉樹MS7.1、蘆山MS7.0、九寨溝MS7.0、長(zhǎng)寧MS6.0等多次中強(qiáng)以上地震前后的區(qū)域性地電場(chǎng)異常變化,開展了大量數(shù)據(jù)分析、信息識(shí)別和異常特征提取等研究工作,逐步積累了一系列非常寶貴的地電場(chǎng)觀測(cè)數(shù)據(jù)及震例信息資料(張學(xué)民等,2006;湯吉等,2010;Huang, 2011;田山等,2012;馬欽忠等,2013;席繼樓等,2016a,2018a,2020).

    “十二五”期間,“中國(guó)地震背景場(chǎng)探測(cè)項(xiàng)目”實(shí)施前后,部分臺(tái)站新建和改建了可觀測(cè)和存儲(chǔ)秒采樣數(shù)據(jù)的地電場(chǎng)觀測(cè)系統(tǒng),主要用于在該臺(tái)站日常運(yùn)維檢測(cè)和異常變化落實(shí)過程中,能夠準(zhǔn)實(shí)時(shí)評(píng)判和校驗(yàn)地電場(chǎng)觀測(cè)數(shù)據(jù)變化的可靠性.研究顯示,相對(duì)于分鐘采樣地電場(chǎng)觀測(cè),在這種基于秒采樣的地電場(chǎng)觀測(cè)數(shù)據(jù)中,包含更為豐富的較短周期地電場(chǎng)變化,能夠比較完整的記錄和再現(xiàn)地震活動(dòng)過程中地電場(chǎng)時(shí)變特性,并且可在較寬頻率范圍(DC~0.5 Hz)內(nèi)更精確反映地電場(chǎng)的動(dòng)態(tài)變化信息(席繼樓等,2016b,2020).

    大武地震臺(tái)是青海省地震局所屬的專業(yè)地震臺(tái)站,位于青海省果洛藏族自治州瑪沁縣大武鎮(zhèn)(E100.22°、N34.49°),地處甘肅、青海和四川三省的交界地段,庫(kù)瑪斷裂帶東段北側(cè),巴顏喀拉山和阿尼瑪卿山之間,距離瑪多MS7.4地震震中約170 km,如圖1所示.同時(shí),大武地震臺(tái)的地電場(chǎng)觀測(cè)系統(tǒng)始建于2007年,并于2013年進(jìn)行了系統(tǒng)性技術(shù)改造,是我國(guó)大陸地區(qū)第一個(gè)改造和建設(shè)完成的可以實(shí)現(xiàn)秒采樣數(shù)據(jù)觀測(cè)和存儲(chǔ)的地電場(chǎng)臺(tái)站.

    本文主要針對(duì)瑪多地震前后,在大武地震臺(tái)觀測(cè)記錄的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)資料,利用傳統(tǒng)的統(tǒng)計(jì)分析和信號(hào)處理方法,就其動(dòng)態(tài)變化特征及其可能的物理機(jī)制開展分析和討論,以期獲取地震活動(dòng)過程中,地電場(chǎng)的場(chǎng)源模型及其可能變化機(jī)理方面的進(jìn)一步深入認(rèn)識(shí).

    圖1 瑪多MS7.4地震震中位置圖1)Fig.1 The epicenter location about Madoi MS7.4 earthquake

    1 大武地震臺(tái)地電場(chǎng)觀測(cè)

    地電場(chǎng)觀測(cè)遵循如下基本觀測(cè)原理:即在地表選定兩個(gè)以上觀測(cè)方位,每個(gè)方位按照?qǐng)D2所示方法,布設(shè)測(cè)量電極M和N,觀測(cè)M和N電極之間的電位差ΔUMN及電極間距LMN,并按照公式(1)所示方法,計(jì)算在該測(cè)量方位的地電場(chǎng)強(qiáng)度EMN,單位記為mV·km-1.

    (1)

    為準(zhǔn)確分析和檢驗(yàn)地電場(chǎng)觀測(cè)數(shù)據(jù)的可靠性及觀測(cè)精度,檢測(cè)和識(shí)別電磁環(huán)境干擾、觀測(cè)場(chǎng)地、系統(tǒng)故障等人為因素產(chǎn)生的非天然電場(chǎng)變化,大武地震臺(tái)的地電場(chǎng)觀測(cè)系統(tǒng)采用了“雙L”型測(cè)量裝置布設(shè)方式,沿著臺(tái)站圍墻的內(nèi)側(cè)布設(shè)和埋設(shè)觀測(cè)外線路,如圖3所示.其中,兩個(gè)正交測(cè)向(實(shí)線所示)約為北東(NE)和北西(NW)方位,第三測(cè)向(虛線所示)接近北南(NS)方位,每個(gè)測(cè)向均布設(shè)長(zhǎng)、短兩種極距的獨(dú)立裝置,長(zhǎng)極距裝置用下標(biāo)L表示、短極距裝置用下標(biāo)S表示,共組合成為6組獨(dú)立的測(cè)道(ch0~ch5).

    圖2 地電場(chǎng)觀測(cè)基本原理示意圖Fig.2 The basic principle of geoelectric field observation

    圖3 大武地震臺(tái)地電場(chǎng)觀測(cè)系統(tǒng)布設(shè)示意圖Fig.3 The layout method of the geoelectric field observation system at Dawu seismic station

    圖3所示的大武地震臺(tái)地電場(chǎng)觀測(cè)系統(tǒng)中,測(cè)量電極采用Pb-PbCl2固體不極化電極,電極之間100%濃度鹽水電位差不大于1.0 mV,地下埋設(shè)深度不小于3 m.測(cè)量外線路采用高絕緣、高防腐性能的鎧裝多芯銅電纜,地下埋設(shè)深度不小于0.5 m.觀測(cè)儀器采用新型網(wǎng)絡(luò)化地電場(chǎng)觀測(cè)儀器“ZD9A-2B地電場(chǎng)儀”,配置6個(gè)高精度數(shù)據(jù)采集通道,依次與圖3所示ch0~ch5測(cè)道配接,主要技術(shù)性能包括:①測(cè)量準(zhǔn)確度優(yōu)于(0.1%讀數(shù)+0.02%滿度);②測(cè)量分辨力優(yōu)于10 μV;③頻帶范圍覆蓋DC~0.1 Hz;④測(cè)量范圍≥±1000.00 mV;⑤動(dòng)態(tài)范圍≥100 dB;⑥工頻共模抑制比≥150 dB;⑦工頻串模抑制比≥100 dB;⑧最大電源功率損耗≤5 W;⑨工作環(huán)境溫度范圍-30~40 ℃(實(shí)驗(yàn)室溫度測(cè)試試驗(yàn)).

    由于秒采樣地電場(chǎng)觀測(cè)的數(shù)據(jù)量比較大,為了全面反映該地電場(chǎng)動(dòng)態(tài)變化全過程,圖4給出了2020年1月—2021年8月期間,大武地震臺(tái)地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)的日均值曲線和日標(biāo)準(zhǔn)偏差曲線(以NES、NWS、NSS三個(gè)測(cè)道為例).其中,圖4a、圖4b和圖4c分別為三個(gè)測(cè)向短極距測(cè)道的日均值曲線,圖4d、圖4e和圖4f分別為對(duì)應(yīng)測(cè)道的日標(biāo)準(zhǔn)偏差曲線.從圖4可以看到,在2020年4月至2021年8月期間,大武地震臺(tái)的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)中出現(xiàn)了比較顯著性的異常變化.

    圖4 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)統(tǒng)計(jì)分析曲線(2020-01-01—2021-08-22)(a) NES日均值; (b) NWS日均值; (c) NSS日均值; (d) NES標(biāo)準(zhǔn)偏差; (e) NWS標(biāo)準(zhǔn)偏差; (f) NSS標(biāo)準(zhǔn)偏差.Fig.4 Statistical analysis curve of second sampling geoelectric field data at Dawu seismic station (from Jan 1, 2020 to Aug 22, 2021)(a) NES daily mean; (b) NWS daily mean; (c) NSS daily mean; (d) NES standard deviation; (e) NWS standard deviation; (f) NSS standard deviation.

    2 地電場(chǎng)觀測(cè)數(shù)據(jù)時(shí)序分析

    2.1 時(shí)序分析方法

    一般來說,地表觀測(cè)的地電場(chǎng)變化主要可分為三類:第一類,信號(hào)源與地電場(chǎng)臺(tái)站的直線距離遠(yuǎn)遠(yuǎn)大于測(cè)區(qū)觀測(cè)裝置系統(tǒng)的幾何尺寸,該地電場(chǎng)變化可視為“遠(yuǎn)源”場(chǎng),具有較強(qiáng)的均勻性變化特征,如外空源產(chǎn)生的地電日變化、地電暴等背景變化;第二類,信號(hào)源與地電場(chǎng)臺(tái)站的直線距離比較接近測(cè)區(qū)觀測(cè)裝置系統(tǒng)的幾何尺寸,該地電場(chǎng)變化可視為“近源”場(chǎng),具有一定的非均勻性和局部性特征,如近場(chǎng)地電磁環(huán)境干擾變化;第三類,當(dāng)信號(hào)源與地電場(chǎng)臺(tái)站的距離介于“遠(yuǎn)源”和“近源”之間時(shí),如果該地電場(chǎng)的時(shí)空變化無明顯局部性特征,則可近似視為“準(zhǔn)遠(yuǎn)源”場(chǎng).

    分析和研究顯示,與地震活動(dòng)性有關(guān)的地電場(chǎng)變化,由于其信號(hào)源距離觀測(cè)臺(tái)站的距離比測(cè)區(qū)尺寸大得多,一般具有“遠(yuǎn)源”場(chǎng)或“準(zhǔn)遠(yuǎn)源”場(chǎng)變化特征(馬欽忠,2008;馬欽忠等,2011).具體表現(xiàn)在,在圖3中所標(biāo)示各測(cè)道記錄到的地電場(chǎng)變化動(dòng)態(tài)過程,應(yīng)具有較強(qiáng)同步性,且同測(cè)向的不同測(cè)道地電場(chǎng)變化具有較強(qiáng)的一致性,或近似一致性.

    本文主要依據(jù)上述理論研究結(jié)果,以及地電場(chǎng)場(chǎng)源特性的分析和識(shí)別方法,基于圖4所示的統(tǒng)計(jì)分析數(shù)據(jù)資料,對(duì)大武地震臺(tái)的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)的變化特征及可能變化機(jī)理,進(jìn)行分析和討論.

    2.2 震前動(dòng)態(tài)變化

    通過逐日查閱大武地震臺(tái)地電場(chǎng)觀測(cè)的工作日志記錄,在排除雷電、磁暴和磁擾等空間因素、降雨等氣象因素、電磁環(huán)境干擾因素、以及觀測(cè)系統(tǒng)故障等影響的情況下,瑪多MS7.4地震之前大武地震臺(tái)的地電場(chǎng)異常變化主要表現(xiàn)為單向階躍變化和雙向擾動(dòng)變化.

    2.2.1 階躍變化

    階躍變化是典型的地電場(chǎng)動(dòng)態(tài)變化之一,在多次大地震之前有過記錄(馬欽忠,2008;馬欽忠等,2011;席繼樓等,2016a,2018a,2020).圖4所示瑪多地震前后的多發(fā)性地電場(chǎng)階躍變化,其時(shí)序變化過程具有一定的顯著性.

    為進(jìn)一步分析和討論該階躍變化的客觀性和可靠性,首先對(duì)大武地震臺(tái)地電場(chǎng)觀測(cè)數(shù)據(jù)及臺(tái)站工作日志等資料,進(jìn)行了逐日分析和核準(zhǔn),明確了圖4中各大幅度階躍變化的發(fā)生時(shí)間,以及電磁環(huán)境干擾、氣象因素以及場(chǎng)地因素等相關(guān)信息,如表1所示.在此基礎(chǔ)上,主要選取具有代表性的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)資料,如圖5—7所示,對(duì)這種階躍變化的時(shí)序特征進(jìn)行量化分析.

    圖5為2020年11月16日秒采樣地電場(chǎng)日變化曲線,即圖4所示多發(fā)性大幅度階躍變化之一.從圖5可以看到,在該階躍變化的前后,均有一個(gè)相對(duì)緩慢的過渡性變化過程,同時(shí)伴隨著小幅度的擾動(dòng)變化,其變化特征主要表現(xiàn)在兩個(gè)方面:①各測(cè)道的該地電場(chǎng)變化,具有較強(qiáng)的同步性和相關(guān)性;②同測(cè)向、不同測(cè)道地電場(chǎng)躍變幅度有一定的差別,比值約為0.3~0.5左右.

    圖6為2020年11月12日地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)的日變化曲線,其中在正常背景變化曲線之上,疊加了一段小幅度階躍變化,持續(xù)時(shí)間約8 h.該階躍變化的基本特征為:①各測(cè)道的該地電場(chǎng)變化,具有較強(qiáng)的同步性和相關(guān)性;②同測(cè)向、不同測(cè)道地電場(chǎng)變化幅度基本一致.

    圖7為2020年6月28日地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)的日變化曲線,其中在正常的背景變化曲線之上,短時(shí)間疊加了一段小幅度的階躍變化,持續(xù)時(shí)間約40 min,其基本特征與圖6所示的階躍變化基本相似.

    圖5 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)大幅度階躍變化曲線(2020-11-16)(a) NEL日變化; (b) NWL日變化; (c) NSL日變化; (d) NES日變化; (e) NWS日變化; (f) NSS日變化.Fig.5 Large amplitude step change curve of second sampling geoelectric field data at Dawu seismic station (Nov 16, 2020)(a) NEL daily variation; (b) NWL daily variation; (c) NSL daily variation; (d) NES daily variation; (e) NWS daily variation; (f) NSS daily variation.

    圖6 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)小幅度階躍變化曲線(2020-11-12)(a) NEL日變化; (b) NWL日變化; (c) NSL日變化; (d) NES日變化; (e) NWS日變化; (f) NSS日變化.Fig.6 Small amplitude step change curve of second sampling geoelectric field data at Dawu seismic station (Nov 12, 2020)(a) NEL daily variation; (b) NWL daily variation; (c) NSL daily variation; (d) NES daily variation; (e) NWS daily variation; (f) NSS daily variation.

    圖7 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)疊加性脈沖變化曲線(2020-06-28)(a) NEL日變化;(b) NWL日變化;(c) NSL日變化;(d) NES日變化;(e) NWS日變化;(f) NSS日變化.Fig.7 Superimposed pulse variation curve of second sampling geoelectric field data at Dawu seismic station (Jun 28, 2020)(a) NEL daily variation; (b) NWL daily variation; (c) NSL daily variation; (d) NES daily variation; (e) NWS daily variation; (f) NSS daily variation.

    依據(jù)大武地震臺(tái)的工作日志記載,2020年10月10日—2021年1月31日期間,臺(tái)站內(nèi)部曾實(shí)施過打鉆作業(yè)兩處,該打鉆作業(yè)位置在B2電極附近,距離東圍墻直線距離約15 m左右,計(jì)劃在此開展鉆孔傾斜觀測(cè),如表1所示.由于在該時(shí)段的地電場(chǎng)觀測(cè)數(shù)據(jù)中,未發(fā)現(xiàn)局部性或頻發(fā)性干擾變化,初步分析認(rèn)為,大武地震臺(tái)的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)中出現(xiàn)的這種顯著性階躍變化,與臺(tái)站內(nèi)部實(shí)施的打鉆作業(yè)過程無明確關(guān)聯(lián)性.

    表1 大武地震臺(tái)地電場(chǎng)大幅度階躍變化的影響因素調(diào)查Table 1 Investigation on influencing factors about large step changes of the geoelectric field at Dawu seismic station

    綜合分析顯示,圖5、圖6和圖7所示的階躍變化過程,有一個(gè)共同的特點(diǎn),就是該變化均具有單向性和變化同步性,以及比較典型的“遠(yuǎn)源”場(chǎng)或“準(zhǔn)遠(yuǎn)源”場(chǎng)變化特性.其波形特征主要表現(xiàn)為:在正常背景變化基礎(chǔ)上,不定期向某一個(gè)方向突然增大,并且在持續(xù)一段時(shí)間以后瞬間減小,隨后逐步恢復(fù)平穩(wěn).分析認(rèn)為,這種單向階躍變化,可能與較強(qiáng)地應(yīng)力和地應(yīng)變持續(xù)作用下,地下流體運(yùn)移過程和滲透作用的突發(fā)性變化有關(guān).

    2.2.2 擾動(dòng)變化

    擾動(dòng)變化也是瑪多MS4.7地震前出現(xiàn)的比較典型地電場(chǎng)動(dòng)態(tài)變化.從2020年6月份開始,在大武地震臺(tái)地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)中,不定期疊加了如圖8所示的較高頻率突發(fā)性擾動(dòng)變化,該變化的出現(xiàn)時(shí)間、持續(xù)時(shí)間段以及變化幅度均具有一定不確定性.但總體來看,各個(gè)測(cè)道地電場(chǎng)變化同步性,以及同測(cè)向、不同測(cè)道地電場(chǎng)變化一致性均比較強(qiáng),變化幅度與當(dāng)天日變化幅度大約在同一數(shù)量級(jí),且與電磁環(huán)境、場(chǎng)地條件及氣象變化等影響因素?zé)o直接關(guān)聯(lián),表現(xiàn)出較強(qiáng)的“遠(yuǎn)源”場(chǎng)變化特性.

    另外,從2021年5月2日開始,各測(cè)道同步出現(xiàn)如圖9所示另一類擾動(dòng)變化,其主要變化特征為:①變化過程持續(xù)性,2021年5月2日至5月24日期間,該擾動(dòng)性變化一直存在,中間未曾間斷過;②脈沖寬度漸變性,由最大約200多秒漸變至約2 s左右;③正反交替性,正反兩個(gè)方向脈沖變化交替出現(xiàn);④變化間隔漸變性,每一組正反向脈沖變化的時(shí)間間隔,經(jīng)歷了由最大約710 s左右到最小約270 s漸變過程;⑤擾動(dòng)變化同步性,各個(gè)測(cè)道擾動(dòng)變化具有較強(qiáng)同步性,且每個(gè)測(cè)道的正反向脈沖幅度基本保持穩(wěn)定.

    針對(duì)圖9b所示NWL測(cè)道與其他測(cè)道的變化不一致現(xiàn)象,2021年5月8日至5月9日期間,大武地震臺(tái)對(duì)觀測(cè)外線路進(jìn)行了全面檢查,并且將連接“ZD9A-2B地電場(chǎng)儀”長(zhǎng)、短極距測(cè)道的信號(hào)線進(jìn)行互換測(cè)試,以檢測(cè)觀測(cè)系統(tǒng)可能存在的相關(guān)問題.測(cè)試結(jié)果顯示,圖9b所示的變化差異性與觀測(cè)系統(tǒng)的工作狀態(tài)無關(guān).通過對(duì)臺(tái)站周邊的環(huán)境因素的仔細(xì)考察,由于A1電極附近的地形地貌比較復(fù)雜(斜坡、斷崖式陡坎等),目前尚不排除由此產(chǎn)生的可能影響.

    綜合分析圖8和圖9所示地電場(chǎng)擾動(dòng)變化,其共同特點(diǎn)為同步性、相關(guān)性和雙向性,各測(cè)道地電場(chǎng)數(shù)據(jù)均在正反兩個(gè)方向同步出現(xiàn)交替性或非交替性的脈沖變化.其中,每一組脈沖變化均表現(xiàn)為突然增大,并在持續(xù)一段時(shí)間以后瞬間截止和快速恢復(fù),基本特征與圖7所示的小幅度階躍變化比較類似.分析認(rèn)為,這種雙向擾動(dòng)變化,可能與較強(qiáng)地應(yīng)力和地應(yīng)變的持續(xù)作用下,地下流體運(yùn)移過程和滲透作用的震顫性變化有關(guān).

    2.3 同震變化

    2021年5月22日瑪多MS7.4地震發(fā)震期間,在大武地震臺(tái)的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)中,出現(xiàn)了圖10所示的同震變化,基本特征為:各個(gè)測(cè)道的地電場(chǎng)變化具有較好的同步性,但變化形態(tài)略具差異性,同測(cè)向、各測(cè)道數(shù)據(jù)變化的峰值時(shí)刻、峰峰值大小等特性也不盡相同.由于大武地震臺(tái)的“ZD9A-2B地電場(chǎng)儀”觀測(cè)不到大于1 Hz及以上較高頻率地電場(chǎng)變化,圖10所示各測(cè)道變化形態(tài)的差異性,可能與高頻信號(hào)欠采樣產(chǎn)生的頻率混疊效應(yīng)有關(guān).

    圖8 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)高頻擾動(dòng)變化曲線(2020-09-06)(a) NEL日變化; (b) NWL日變化; (c) NSL日變化; (d) NES日變化; (e) NWS日變化; (f) NSS日變化.Fig.8 High frequency disturbance variation curve of second sampling geoelectric field data at Dawu seismic station (Sep 6, 2020)(a) NEL daily variation; (b) NWL daily variation; (c) NSL daily variation; (d) NES daily variation; (e) NWS daily variation; (f) NSS daily variation.

    圖9 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)持續(xù)性擾動(dòng)變化曲線(2021-05-03)(a) NEL日變化; (b) NWL日變化; (c) NSL日變化; (d) NES日變化; (e) NWS日變化;(f) NSS日變化.Fig.9 Continuous disturbance variation curve of second sampling geoelectric field data at Dawu seismic station (May 3, 2021)(a) NEL daily variation; (b) NWL daily variation; (c) NSL daily variation; (d) NES daily variation; (e) NWS daily variation; (f) NSS daily variation.

    圖10 大武地震臺(tái)秒采樣地電場(chǎng)同震變化曲線(2021-05-22)(a) NEL同震變化; (b) NWL同震變化; (c) NSL同震變化; (d) NES同震變化; (e) NWS同震變化; (f) NSS同震變化.Fig.10 Co-seismic variation curve of second sampling geoelectric field data at Dawu seismic station (May 22, 2021)(a) NEL co-seismic variation; (b) NWL co-seismic variation; (c) NSL co-seismic variation; (d) NES co-seismic variation; (e) NWS co-seismic variation; (f) NSS co-seismic variation.

    圖11 大武地震臺(tái)秒采樣地電(磁)場(chǎng)同震變化對(duì)比分析曲線(a) N地震波形; (b) E地震波形; (c) Z地震波形; (d) H地磁場(chǎng)變化; (e) Z地磁場(chǎng)變化; (f) D地磁場(chǎng)變化; (g) 合成地電場(chǎng)變化.Fig.11 Co-seismic variation comparative analysis curve of second sampling geoelectric (magnetic) field data at Dawu seismic station(a) N seismic waveform; (b) E seismic waveform; (c) Z seismic waveform; (d) H geomagnetic field variation; (e) Z geomagnetic field variation; (f) D geomagnetic field variation; (g) Synthetic geoelectric field variation.

    圖12 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)典型震后變化曲線(2021-08-08)(a) NEL日變化; (b) NWL日變化; (c) NSL日變化; (d) NES日變化; (e) NWS日變化; (f) NSS日變化.Fig.12 Typical post-earthquake variation curve of second sampling geoelectric field data at Dawu seismic station (Aug 8, 2021)(a) NEL daily variation; (b) NWL daily variation; (c) NSL daily variation; (d) NES daily variation; (e) NWS daily variation; (f) NSS daily variation.

    圖13 大武地震臺(tái)地電場(chǎng)秒采樣數(shù)據(jù)變化功率譜密度(2020-01-01—2021-08-22)(a) NEL功率譜密度; (b) NWL功率譜密度; (c) NSL功率譜密度; (d) NES功率譜密度; (e) NWS功率譜密度; (f) NSS功率譜密度.Fig.13 Power spectral density of second sampling geoelectric field data at Dawu seismic station (from Jan 1, 2020 to Aug 22, 2021)(a) NEL power spectral density; (b) NWL power spectral density; (c) NSL power spectral density; (d) NES power spectral density; (e) NWS power spectral density; (f) NSS power spectral density.

    為進(jìn)一步分析圖10所示地電場(chǎng)同震變化的基本特性,首先對(duì)該同震變化數(shù)據(jù)進(jìn)行了矢量合成(以NES和NWS兩個(gè)正交測(cè)道為例),計(jì)算結(jié)果如圖11g所示.與此同時(shí),仔細(xì)分析和核對(duì)大武地震臺(tái)“GM4磁通門磁力儀”記錄的地磁場(chǎng)秒采樣觀測(cè)數(shù)據(jù),截取了相關(guān)時(shí)間段的H、Z、D三個(gè)分量的相對(duì)地磁變化,如圖11d—f所示.特別是,專門提取了大武地震臺(tái)“CTS-1E寬頻帶地震計(jì)”記錄到的瑪多地震的波形信息(1 Hz抽樣),如圖11a—c所示.

    在圖11中,為完整呈現(xiàn)和分析地電場(chǎng)的同震變化,分別截取了瑪多MS7.4主震前后的500 s測(cè)震、地電場(chǎng)和地磁場(chǎng)同步觀測(cè)數(shù)據(jù)資料.其中,GM4磁通門磁力儀和CTS-1E寬頻帶地震計(jì)的時(shí)間信息,已經(jīng)利用全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System,GNSS)進(jìn)行授時(shí),與世界時(shí)(Greenwich Mean Time,GMT)自動(dòng)同步;ZD9A-2B地電場(chǎng)儀的時(shí)間信息,則利用地電場(chǎng)變化與地磁場(chǎng)變化的時(shí)序相關(guān)性,按照GM4磁通門磁力儀時(shí)間信息進(jìn)行了校對(duì)和同步.

    對(duì)比分析圖11所示的同震變化曲線,其基本特征為:①地電場(chǎng)變化和地磁場(chǎng)變化具有較強(qiáng)的時(shí)序相關(guān)性;②地電場(chǎng)和地磁場(chǎng)的主要變化過程與直達(dá)橫波(Sg波)基本同步;③與直達(dá)縱波(Pg波)相關(guān)的地電場(chǎng)變化幅度比較微弱;④地電(磁)場(chǎng)變化尾波衰減均比較快,未識(shí)別到與面波相關(guān)的變化過程.

    由于圖11顯示的地電(磁)場(chǎng)變化形態(tài)與地震波形具有較大差異,分析認(rèn)為,圖10所示該頻段(DC~0.5 Hz)地電場(chǎng)的同震變化,主要與強(qiáng)震動(dòng)過程中,在橫波(S波)作用下,臺(tái)站及周邊地區(qū)地下流體振蕩性運(yùn)移及滲流作用有關(guān),與地表介質(zhì)的彈性振動(dòng)過程的關(guān)聯(lián)性不是太大.

    2.4 震后變化

    據(jù)中國(guó)地震臺(tái)網(wǎng)中心測(cè)定,2021年8月13日12時(shí)21分35秒,在青海省果洛州瑪多縣(E97.54°,N34.58°)又發(fā)生了一次MS5.8地震,震源深度8 km,距離瑪多MS7.4地震震中的直線距離約73 km,如圖1所示.由于兩次地震的可能關(guān)聯(lián)性,在這里特別對(duì)大武地震臺(tái)2021年5月份至8月份的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù),進(jìn)行了補(bǔ)充分析和討論.

    逐日分析結(jié)果顯示,在瑪多MS7.4地震的震后時(shí)間段,大武地震臺(tái)的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù),主要經(jīng)歷了如下變化過程:①2021年5月22日瑪多MS7.4地震主震過后,圖9所示持續(xù)性擾動(dòng)變化于2021年5月24日開始大幅度減小,并在5月27日之后基本恢復(fù)平穩(wěn);②在2021年5月27日—7月15日期間,在排除降雨、雷電、磁暴和磁擾等影響因素的情況下,該地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)未發(fā)現(xiàn)明顯異常變化;③從2021年7月16日開始,該地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)中,又一次不定期出現(xiàn)多種形態(tài)的突發(fā)性脈沖變化和擾動(dòng)變化;④從2021年7月26日開始,該地電場(chǎng)觀測(cè)數(shù)據(jù)中再次出現(xiàn)如圖9所示持續(xù)性擾動(dòng)變化過程;⑤2021年8月7以后,在持續(xù)性擾動(dòng)變化波形中,每天不定期疊加了如圖12所示的較高頻率、較大幅度的隨機(jī)性擾動(dòng)變化.

    從上述分析結(jié)果可以看到,在2021年8月13日發(fā)生的瑪多MS5.8地震前后,大武地震臺(tái)的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)出現(xiàn)了與瑪多MS7.4地震前后比較類似的動(dòng)態(tài)變化過程,顯示了在該發(fā)震區(qū)域的不同地震前后,地電場(chǎng)的異常特征及演化過程可能具有的重現(xiàn)性.

    3 地電場(chǎng)觀測(cè)數(shù)據(jù)頻譜分析

    結(jié)合前面的時(shí)序分析結(jié)果,基于開源計(jì)算機(jī)編程軟件Numerical Python的快速傅里葉變換(Fast Fourier Transform,F(xiàn)FT)數(shù)字信號(hào)處理算法,對(duì)大武地震臺(tái)的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)進(jìn)行逐日頻譜分析.同時(shí),利用周期圖方法,對(duì)每一天的不大于1小時(shí)周期的地電場(chǎng)功率譜密度(Power Spectral Density,PSD)進(jìn)行計(jì)算、分析和討論.

    主要分析過程如下:①針對(duì)每一天的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù),以7200 s周期長(zhǎng)度加漢寧窗(Hanning)函數(shù);②采用單邊快速傅里葉變換(FFT)算法,分段計(jì)算加窗函數(shù)之后的地電場(chǎng)離散傅里葉變化數(shù)據(jù)(Discrete Fourier Transform,DFT),并進(jìn)一步計(jì)算其功率譜密度(PSD);③將每天加窗函數(shù)以后分段得到的多組功率譜密度,逐頻點(diǎn)計(jì)算其算術(shù)平均值;④將計(jì)算得到的功率譜密度(PSD)數(shù)值換算為分貝值(dB)表示;⑤利用Matplotlib Python繪圖函數(shù),針對(duì)上述分析和計(jì)算結(jié)果,繪制地電場(chǎng)功率譜密度的時(shí)頻圖像.

    圖13為2020年1月1日至2021年8月22日(橫坐標(biāo)每段標(biāo)識(shí)的間隔時(shí)間為四個(gè)月)期間,大武地震臺(tái)地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)功率譜密度計(jì)算結(jié)果,頻率范圍約為0.0002~0.5 Hz(對(duì)數(shù)坐標(biāo)).其中,2020年4月8日—5月30日期間數(shù)據(jù)空缺(未及時(shí)收集),為了繪圖方便,用2020年4月7日數(shù)據(jù)進(jìn)行填充.

    分析圖13所示地電場(chǎng)功率譜密度數(shù)據(jù),其基本變化特征為:①各測(cè)道的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)的功率譜密度計(jì)算結(jié)果,具有較強(qiáng)一致性;②從2020年4月份開始,地電場(chǎng)功率譜密度較之前逐步抬升,最高可達(dá)50~60 dB左右,并經(jīng)歷了“抬升-回落-平穩(wěn)-再抬升”變化過程,該變化過程與時(shí)間序列分析結(jié)果基本對(duì)應(yīng);③從瑪多MS7.4地震發(fā)生之后,至瑪多MS5.8地震發(fā)生之前,地電場(chǎng)功率譜密度再一次經(jīng)歷了“(抬升)-回落-平穩(wěn)-再抬升”的變化過程,該變化過程與瑪多MS7.4地震之前的變化過程基本類似.

    從上述功率譜密度的變化特征可以看到,在瑪多MS7.4地震前后,大武地震臺(tái)秒采樣地電場(chǎng)動(dòng)態(tài)變化過程,在“時(shí)間域”和“頻率域”分析結(jié)果具有比較強(qiáng)的一致性,進(jìn)一步檢驗(yàn)了該地電場(chǎng)動(dòng)態(tài)變化的客觀性和可靠性.與此同時(shí),圖13所示的分析結(jié)果,也從“頻率域”再一次佐證,在該發(fā)震區(qū)域的不同地震前后,地電場(chǎng)動(dòng)態(tài)變化特征及演化過程可能具有的重現(xiàn)性.

    4 結(jié)論與討論

    地電場(chǎng)是重要的地球物理場(chǎng),本文利用傳統(tǒng)的時(shí)間序列統(tǒng)計(jì)分析方法和FFT信號(hào)處理方法,從“時(shí)間域”和“頻率域”兩個(gè)方面,對(duì)瑪多MS7.4地震前后,距離震中約170 km的青海省果洛州大武地震臺(tái)記錄到的秒采樣地電場(chǎng)動(dòng)態(tài)變化的基本特征,進(jìn)行了多方面的分析和討論.

    研究結(jié)果顯示:(1)瑪多地震前后,大武地震臺(tái)記錄到了比較典型的地電場(chǎng)震前異常變化和同震變化,該變化過程在不同的觀測(cè)裝置下,具有較強(qiáng)的同步性和相關(guān)性;(2)在瑪多地震發(fā)生之前,大武地震臺(tái)的地電場(chǎng)異常變化,主要表現(xiàn)為不定期交叉和交疊出現(xiàn)的單向階躍變化和雙向擾動(dòng)變化,其中雙向擾動(dòng)變化延續(xù)更長(zhǎng)時(shí)間段,并在不同階段,表現(xiàn)為突發(fā)性或持續(xù)性變化特征;(3)瑪多地震的地電場(chǎng)同震變化,其主要變化過程與橫波(S波)相關(guān),而與縱波(P波)相關(guān)的變化幅度比較微弱,未識(shí)別到與面波相關(guān)的波形記錄;(4)在瑪多地震前后,大武地震臺(tái)的秒采樣地電場(chǎng)動(dòng)態(tài)變化過程,在“時(shí)間域”和“頻率域”具有較強(qiáng)一致性.

    綜合分析認(rèn)為,在排除電磁暴、雷電等空間電磁信號(hào)變化的情況下,在瑪多MS7.4地震前后,大武地震臺(tái)秒采樣(DC~0.5 Hz頻段)地電場(chǎng)動(dòng)態(tài)變化過程,主要與較強(qiáng)地應(yīng)力和地應(yīng)變作用下,地下流體在構(gòu)造斷層、孔隙及裂隙中的運(yùn)移過程和滲流作用的突發(fā)性、振蕩性和震顫性變化有關(guān).其中:①單向階躍變化,可能與大區(qū)域較強(qiáng)應(yīng)力和地應(yīng)變持續(xù)作用下,地下流體的突發(fā)性變化有關(guān),并且該突發(fā)性變化可能不定期源自震源構(gòu)造體、以及臺(tái)站及周邊區(qū)域構(gòu)造斷層帶,從而產(chǎn)生“遠(yuǎn)源”或“準(zhǔn)遠(yuǎn)源”地電場(chǎng)階躍變化;②雙向擾動(dòng)變化,可能與裂隙發(fā)育過程、以及非震構(gòu)造斷層帶的蠕滑過程中,地下流體震顫性變化有關(guān),該震顫性變化可能源自震源構(gòu)造體、以及震源臨近區(qū)域的斷層構(gòu)造帶等,從而產(chǎn)生“遠(yuǎn)源”或“準(zhǔn)遠(yuǎn)源”地電場(chǎng)擾動(dòng)性變化;③地電場(chǎng)同震變化則主要與橫波(S波)到達(dá)臺(tái)站及周邊區(qū)域時(shí),該區(qū)域的地下流體在構(gòu)造斷層及孔隙、裂隙中的振蕩性變化過程有關(guān).

    青海省大武地震臺(tái)位于庫(kù)瑪斷裂帶東段北側(cè),巴顏喀拉山和阿尼瑪卿山之間,與瑪多MS7.4地震的發(fā)震構(gòu)造(昆侖山口—江錯(cuò)斷裂),同屬于巴顏喀拉塊體邊界構(gòu)造帶.由于這種地理位置和地質(zhì)構(gòu)造的特殊性,瑪多地震活動(dòng)過程中,在較強(qiáng)地應(yīng)力和地應(yīng)變作用下,地下流體在發(fā)震構(gòu)造體和臺(tái)站周邊構(gòu)造斷層中發(fā)生運(yùn)移和滲流的可能性均比較大,可以較好的支持上述分析和討論過程及結(jié)論.特別是,當(dāng)?shù)叵铝黧w在臺(tái)站周邊所在構(gòu)造斷層中突發(fā)性運(yùn)移和滲流時(shí),由此產(chǎn)生的地電場(chǎng)變化過程,在高導(dǎo)介質(zhì)的近距離傳播時(shí),由地面記錄到的相關(guān)動(dòng)態(tài)變化可能具有一定的顯著性.從這個(gè)角度來看,本文所討論的瑪多MS7.4地震前后,大武地震臺(tái)秒采樣地電場(chǎng)動(dòng)態(tài)變化過程及變化特征,以及所揭示的可能物理機(jī)制和變化機(jī)理等,對(duì)于地電場(chǎng)觀測(cè)臺(tái)站勘址、建設(shè)及觀測(cè)數(shù)據(jù)分析和應(yīng)用,以及地電場(chǎng)觀測(cè)技術(shù)系統(tǒng)的進(jìn)一步研究和發(fā)展,均具有一定的啟迪和借鑒意義.

    總體而言,在瑪多MS7.4地震前后,大武地震臺(tái)記錄到的地電場(chǎng)震前變化和同震變化,是我國(guó)在7.0級(jí)以上的大地震前后,首次近距離觀測(cè)記錄的秒采樣地電場(chǎng)動(dòng)態(tài)變化信息,相關(guān)數(shù)據(jù)資料彌足珍貴.本文通過對(duì)該數(shù)據(jù)及其動(dòng)態(tài)變化信息的綜合性分析和討論,在地電場(chǎng)場(chǎng)源特性及變化機(jī)理等方面,獲得了一些比較明確的認(rèn)識(shí),積累了一定的數(shù)據(jù)、資料和經(jīng)驗(yàn).該研究結(jié)果,可望在地震監(jiān)測(cè)、預(yù)測(cè)研究方面發(fā)揮積極的參考、促進(jìn)和推進(jìn)作用.

    致謝感謝上海市地震局馬欽忠研究員、甘肅省地震局譚大誠(chéng)研究員的多方面交流和指導(dǎo),感謝云南省地震局李文新高級(jí)工程師,重慶市地震局陳敏高級(jí)工程師、董蕾工程師,四川省地震局顏曉曄高級(jí)工程師,甘肅省地震局牛延平高級(jí)工程師、張彩艷高級(jí)工程師等提供的不同區(qū)域和時(shí)段的地電場(chǎng)秒采樣觀測(cè)數(shù)據(jù)資料,感謝中國(guó)地震局地震預(yù)測(cè)研究所王斌高級(jí)工程師提供的測(cè)震數(shù)據(jù)資料,感謝審稿專家及編輯部的全方位指導(dǎo)意見和建議.

    猜你喜歡
    大武瑪多階躍
    另一種時(shí)間觀
    讀者(2023年3期)2023-02-28 08:24:28
    音頻式自媒體品牌課:微課開發(fā)新選擇——以北大武志紅的心理學(xué)課為例
    基于階躍雙包層光纖的螺旋型光纖傳感器
    2021年瑪多MS7.4地震的深部構(gòu)造背景
    2021年5月22日青海瑪多MS 7.4地震總結(jié)
    撤離
    探討單位階躍信號(hào)的教學(xué)
    從《大武》“樂”看戲劇教化人心之能效
    戲曲研究(2017年4期)2017-05-31 07:39:38
    果洛州瑪沁縣大武鎮(zhèn)商品豬旋毛蟲感染情況調(diào)查報(bào)告
    一種階躍函數(shù)在矩形時(shí)間窗口頻域特性的分析方法
    亚洲电影在线观看av| 不卡视频在线观看欧美| 男人狂女人下面高潮的视频| 国产免费视频播放在线视频 | 777米奇影视久久| 国产极品天堂在线| 亚洲经典国产精华液单| 可以在线观看毛片的网站| 国产乱来视频区| 最后的刺客免费高清国语| 日本爱情动作片www.在线观看| 日韩制服骚丝袜av| 亚洲成人精品中文字幕电影| 一区二区三区高清视频在线| 精品少妇黑人巨大在线播放| 夜夜爽夜夜爽视频| 国产一级毛片在线| 久久精品夜色国产| 亚洲综合精品二区| 日韩欧美精品免费久久| 亚洲欧美日韩卡通动漫| 免费观看性生交大片5| 91久久精品国产一区二区三区| 国产成人精品婷婷| 亚洲一级一片aⅴ在线观看| 搡女人真爽免费视频火全软件| 亚洲国产精品sss在线观看| 亚洲av国产av综合av卡| 日韩av不卡免费在线播放| 国产白丝娇喘喷水9色精品| 中文资源天堂在线| 亚洲av成人精品一区久久| 亚洲久久久久久中文字幕| 免费看av在线观看网站| 秋霞在线观看毛片| 禁无遮挡网站| 亚洲综合精品二区| 91久久精品国产一区二区成人| 久久久久久国产a免费观看| 精品一区二区免费观看| 久久久午夜欧美精品| 高清视频免费观看一区二区 | 久久久亚洲精品成人影院| 91久久精品电影网| 欧美区成人在线视频| 国产一区二区三区综合在线观看 | 国产淫片久久久久久久久| 国产av国产精品国产| 青春草亚洲视频在线观看| av一本久久久久| 激情 狠狠 欧美| 亚洲av福利一区| 亚洲熟妇中文字幕五十中出| 亚洲综合色惰| 亚洲av成人精品一二三区| 亚洲欧美成人综合另类久久久| 99久久精品热视频| 国产不卡一卡二| 日韩欧美精品v在线| 亚洲成人精品中文字幕电影| 少妇人妻一区二区三区视频| 日韩人妻高清精品专区| 亚洲精品乱码久久久久久按摩| 男的添女的下面高潮视频| 成人亚洲精品一区在线观看 | av网站免费在线观看视频 | 91av网一区二区| 国产精品福利在线免费观看| 日本与韩国留学比较| 亚洲av不卡在线观看| 国产亚洲精品久久久com| 欧美成人一区二区免费高清观看| 欧美zozozo另类| 男人爽女人下面视频在线观看| 秋霞在线观看毛片| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品一,二区| 一级爰片在线观看| 又爽又黄a免费视频| 男女那种视频在线观看| 国产在视频线在精品| 99久久人妻综合| 久久久久久久国产电影| 精品久久久久久成人av| av黄色大香蕉| 午夜精品一区二区三区免费看| 亚洲在线观看片| 国产激情偷乱视频一区二区| 精品一区二区三卡| 免费观看精品视频网站| 国产伦精品一区二区三区四那| 久久久精品免费免费高清| 久久久精品免费免费高清| 免费av观看视频| 久久久久久久久久成人| www.av在线官网国产| 成人亚洲精品一区在线观看 | 亚洲精品成人久久久久久| 免费观看a级毛片全部| 日韩伦理黄色片| 国产伦精品一区二区三区视频9| 热99在线观看视频| 亚州av有码| 男女边吃奶边做爰视频| 舔av片在线| 免费观看无遮挡的男女| 午夜福利在线观看吧| 亚洲精品国产av成人精品| 国产精品一区www在线观看| 免费在线观看成人毛片| 蜜桃亚洲精品一区二区三区| 99热这里只有精品一区| 亚洲人成网站高清观看| 欧美高清成人免费视频www| 熟妇人妻久久中文字幕3abv| 成人亚洲精品一区在线观看 | 男女国产视频网站| 欧美性感艳星| 性插视频无遮挡在线免费观看| 亚洲四区av| 久久精品夜色国产| 国产v大片淫在线免费观看| av线在线观看网站| 一级爰片在线观看| av在线亚洲专区| 国产高清不卡午夜福利| 亚洲精品视频女| 男女下面进入的视频免费午夜| 九色成人免费人妻av| 一个人免费在线观看电影| 成人无遮挡网站| 最近2019中文字幕mv第一页| 亚洲av电影不卡..在线观看| 国产在线一区二区三区精| 丝袜喷水一区| 精品久久久久久久人妻蜜臀av| 国产黄片美女视频| 午夜福利网站1000一区二区三区| 日本色播在线视频| 97超视频在线观看视频| 色综合站精品国产| 一个人观看的视频www高清免费观看| 成年人午夜在线观看视频 | 成人午夜精彩视频在线观看| 国产伦在线观看视频一区| 少妇高潮的动态图| 国产爱豆传媒在线观看| av国产免费在线观看| 精品久久久久久电影网| 精品一区在线观看国产| 欧美zozozo另类| 2021少妇久久久久久久久久久| 秋霞在线观看毛片| 丝瓜视频免费看黄片| 联通29元200g的流量卡| 日韩在线高清观看一区二区三区| 免费无遮挡裸体视频| 1000部很黄的大片| 欧美日韩在线观看h| 亚洲一级一片aⅴ在线观看| 久久久精品欧美日韩精品| 69人妻影院| 中文乱码字字幕精品一区二区三区 | 久久久精品94久久精品| 一区二区三区四区激情视频| 美女主播在线视频| 日韩伦理黄色片| 久久久精品欧美日韩精品| 97热精品久久久久久| 国产精品一二三区在线看| 精品亚洲乱码少妇综合久久| 看黄色毛片网站| 黄色日韩在线| 日韩av不卡免费在线播放| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品色激情综合| 日本av手机在线免费观看| 久久久久网色| 亚洲图色成人| 日日摸夜夜添夜夜添av毛片| 肉色欧美久久久久久久蜜桃 | 午夜福利高清视频| 看非洲黑人一级黄片| 亚洲欧洲日产国产| 精品久久久久久久末码| 日韩一区二区三区影片| 中文在线观看免费www的网站| 久久久久久久久久成人| 日本黄大片高清| 精品久久久久久久久亚洲| 国产黄a三级三级三级人| 国产av国产精品国产| 成人性生交大片免费视频hd| 91精品国产九色| 精品久久久久久久人妻蜜臀av| 女人久久www免费人成看片| 男人舔奶头视频| 最近最新中文字幕免费大全7| 中国国产av一级| 国产精品久久视频播放| 熟女人妻精品中文字幕| 尤物成人国产欧美一区二区三区| 小蜜桃在线观看免费完整版高清| 美女脱内裤让男人舔精品视频| 国产精品福利在线免费观看| 国产伦精品一区二区三区视频9| 听说在线观看完整版免费高清| 国模一区二区三区四区视频| 一级黄片播放器| 十八禁国产超污无遮挡网站| 一个人看的www免费观看视频| 日韩精品青青久久久久久| 亚洲av免费在线观看| 国产高清国产精品国产三级 | 亚洲,欧美,日韩| 国产永久视频网站| 亚洲综合色惰| 国产精品女同一区二区软件| av女优亚洲男人天堂| av黄色大香蕉| 中文资源天堂在线| 国产免费福利视频在线观看| 日本wwww免费看| 综合色丁香网| 亚洲欧美中文字幕日韩二区| 直男gayav资源| 天堂av国产一区二区熟女人妻| 国产男女超爽视频在线观看| 哪个播放器可以免费观看大片| 色哟哟·www| 国产在视频线精品| 欧美日韩亚洲高清精品| 久久久久久久大尺度免费视频| 观看免费一级毛片| 国产成人一区二区在线| 国产乱人视频| 亚洲av.av天堂| 免费看美女性在线毛片视频| 大香蕉97超碰在线| 精品99又大又爽又粗少妇毛片| 在线播放无遮挡| 男人舔奶头视频| 在线免费观看的www视频| 69人妻影院| 韩国高清视频一区二区三区| 夜夜爽夜夜爽视频| 色网站视频免费| 3wmmmm亚洲av在线观看| 欧美日韩视频高清一区二区三区二| 国产精品美女特级片免费视频播放器| 国产精品麻豆人妻色哟哟久久 | 麻豆国产97在线/欧美| 欧美成人一区二区免费高清观看| 特大巨黑吊av在线直播| 亚洲人成网站在线观看播放| 成人毛片60女人毛片免费| 美女国产视频在线观看| 色播亚洲综合网| 亚洲精品456在线播放app| 国产精品不卡视频一区二区| 国内少妇人妻偷人精品xxx网站| 国产精品人妻久久久影院| 插阴视频在线观看视频| 免费黄网站久久成人精品| 亚洲最大成人中文| 男女啪啪激烈高潮av片| 国产精品日韩av在线免费观看| 中文乱码字字幕精品一区二区三区 | 成年av动漫网址| 少妇熟女aⅴ在线视频| av在线观看视频网站免费| 我要看日韩黄色一级片| 男插女下体视频免费在线播放| av线在线观看网站| 大香蕉久久网| av免费观看日本| 亚洲精品久久午夜乱码| 久久人人爽人人片av| 国产在线男女| 国产真实伦视频高清在线观看| 久久久久久久午夜电影| 亚洲国产色片| 久久久a久久爽久久v久久| 最近中文字幕2019免费版| 日韩欧美精品v在线| 视频中文字幕在线观看| 男人爽女人下面视频在线观看| 青春草视频在线免费观看| 欧美三级亚洲精品| 国产一区二区三区av在线| 97在线视频观看| 好男人在线观看高清免费视频| 久久久久久久亚洲中文字幕| 美女脱内裤让男人舔精品视频| 日韩不卡一区二区三区视频在线| 久久精品久久精品一区二区三区| 国产精品蜜桃在线观看| 久久亚洲国产成人精品v| 亚洲精品乱码久久久久久按摩| 七月丁香在线播放| 久久久久九九精品影院| 国产免费又黄又爽又色| 肉色欧美久久久久久久蜜桃 | 草草在线视频免费看| 18禁在线无遮挡免费观看视频| 成人亚洲精品av一区二区| 免费大片黄手机在线观看| 精华霜和精华液先用哪个| 91久久精品国产一区二区三区| 麻豆乱淫一区二区| 天堂av国产一区二区熟女人妻| 国内精品美女久久久久久| 色吧在线观看| 一级毛片我不卡| 真实男女啪啪啪动态图| 色网站视频免费| 美女黄网站色视频| 在现免费观看毛片| 男女那种视频在线观看| 99视频精品全部免费 在线| 亚州av有码| 亚洲欧美中文字幕日韩二区| 啦啦啦韩国在线观看视频| 久久热精品热| 久久久久久国产a免费观看| 水蜜桃什么品种好| 国产伦在线观看视频一区| 国产真实伦视频高清在线观看| av在线天堂中文字幕| 婷婷色av中文字幕| 色综合色国产| 亚洲国产精品专区欧美| 色播亚洲综合网| 欧美高清性xxxxhd video| 欧美日韩综合久久久久久| 我的女老师完整版在线观看| av国产久精品久网站免费入址| 夜夜爽夜夜爽视频| 成人毛片a级毛片在线播放| 久久久国产一区二区| 久久精品夜色国产| 欧美+日韩+精品| 三级国产精品欧美在线观看| 搞女人的毛片| 又爽又黄无遮挡网站| 日韩伦理黄色片| 人体艺术视频欧美日本| 亚洲av中文字字幕乱码综合| 国产在线一区二区三区精| 成人一区二区视频在线观看| 老司机影院毛片| 美女主播在线视频| 嫩草影院入口| 久久久久久久亚洲中文字幕| 亚洲av中文av极速乱| 国产黄色小视频在线观看| 熟妇人妻不卡中文字幕| 久久久久久久久久人人人人人人| 搡女人真爽免费视频火全软件| 亚洲精品影视一区二区三区av| 又黄又爽又刺激的免费视频.| 国产免费福利视频在线观看| 91久久精品国产一区二区三区| 婷婷色av中文字幕| 少妇高潮的动态图| 五月伊人婷婷丁香| videossex国产| 中文字幕久久专区| 在线观看av片永久免费下载| 最近视频中文字幕2019在线8| 建设人人有责人人尽责人人享有的 | 免费黄网站久久成人精品| 国产伦在线观看视频一区| 最新中文字幕久久久久| 最近手机中文字幕大全| 欧美97在线视频| 国产精品久久久久久久电影| 一边亲一边摸免费视频| 韩国高清视频一区二区三区| 亚洲经典国产精华液单| 国产亚洲91精品色在线| 国产成人aa在线观看| 尾随美女入室| 日本wwww免费看| 麻豆成人午夜福利视频| 亚洲无线观看免费| 国产伦精品一区二区三区视频9| 在线免费十八禁| 人妻一区二区av| 亚洲怡红院男人天堂| 亚洲精品日本国产第一区| 麻豆成人av视频| eeuss影院久久| 男人和女人高潮做爰伦理| 中文资源天堂在线| 男人舔女人下体高潮全视频| 亚洲最大成人手机在线| 少妇的逼水好多| 精品国产三级普通话版| 在线免费十八禁| 亚洲性久久影院| 丰满人妻一区二区三区视频av| 色视频www国产| 在线天堂最新版资源| 乱人视频在线观看| 亚洲国产精品国产精品| 免费观看av网站的网址| 99久国产av精品国产电影| 亚洲av成人精品一区久久| 日本与韩国留学比较| 一级毛片我不卡| 麻豆成人av视频| 国产欧美另类精品又又久久亚洲欧美| 51国产日韩欧美| 欧美不卡视频在线免费观看| 又粗又硬又长又爽又黄的视频| 天堂俺去俺来也www色官网 | 美女主播在线视频| 国产精品久久视频播放| 成年女人在线观看亚洲视频 | 久久国产乱子免费精品| 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区| 午夜激情久久久久久久| 欧美日韩亚洲高清精品| 人妻夜夜爽99麻豆av| 国产高清有码在线观看视频| 又黄又爽又刺激的免费视频.| 日本与韩国留学比较| 午夜免费激情av| 日韩成人伦理影院| 久久久久久久久大av| 国产麻豆成人av免费视频| h日本视频在线播放| 少妇的逼好多水| 亚洲精品国产av蜜桃| 边亲边吃奶的免费视频| 欧美zozozo另类| 欧美一级a爱片免费观看看| 日本午夜av视频| 久久久久精品性色| 十八禁国产超污无遮挡网站| 一区二区三区高清视频在线| 国产精品女同一区二区软件| 亚洲精品第二区| 国产精品一区二区在线观看99 | 欧美日韩国产mv在线观看视频 | 青青草视频在线视频观看| 十八禁国产超污无遮挡网站| 美女大奶头视频| 人人妻人人看人人澡| 九九久久精品国产亚洲av麻豆| 少妇人妻精品综合一区二区| 老司机影院毛片| 免费看美女性在线毛片视频| 成人美女网站在线观看视频| 五月伊人婷婷丁香| 国产人妻一区二区三区在| 人妻制服诱惑在线中文字幕| 九色成人免费人妻av| 九九爱精品视频在线观看| 亚洲精品亚洲一区二区| 亚洲va在线va天堂va国产| 一个人看视频在线观看www免费| 久久人人爽人人爽人人片va| 免费黄网站久久成人精品| 在线a可以看的网站| 一级二级三级毛片免费看| 免费少妇av软件| 26uuu在线亚洲综合色| 国产美女午夜福利| 大陆偷拍与自拍| 免费看美女性在线毛片视频| 亚洲在线观看片| 日韩中字成人| 97精品久久久久久久久久精品| 在线a可以看的网站| 超碰av人人做人人爽久久| 欧美97在线视频| 免费av毛片视频| 久久草成人影院| freevideosex欧美| 日本欧美国产在线视频| 国产毛片a区久久久久| 色综合站精品国产| 亚洲在线观看片| 我的老师免费观看完整版| 六月丁香七月| 久久久久免费精品人妻一区二区| 99久久精品热视频| 国产精品伦人一区二区| 男的添女的下面高潮视频| 天天躁日日操中文字幕| 中文字幕av在线有码专区| 日韩伦理黄色片| 久久午夜福利片| 亚洲av中文字字幕乱码综合| 精品人妻偷拍中文字幕| 春色校园在线视频观看| 亚洲最大成人中文| 亚洲美女视频黄频| 2021少妇久久久久久久久久久| 日韩电影二区| 久久精品人妻少妇| 国产黄片视频在线免费观看| 80岁老熟妇乱子伦牲交| 黑人高潮一二区| 精品人妻一区二区三区麻豆| 91av网一区二区| 最近的中文字幕免费完整| 国产熟女欧美一区二区| 一级毛片久久久久久久久女| 国产在线男女| 国产亚洲午夜精品一区二区久久 | 欧美日韩亚洲高清精品| 午夜日本视频在线| 久久久久精品性色| 精品一区二区三区人妻视频| 一区二区三区免费毛片| 国产精品.久久久| 啦啦啦中文免费视频观看日本| av在线亚洲专区| 最近视频中文字幕2019在线8| 国产一区有黄有色的免费视频 | 亚洲aⅴ乱码一区二区在线播放| 亚洲性久久影院| 男女国产视频网站| 国产欧美日韩精品一区二区| 久久人人爽人人爽人人片va| 观看免费一级毛片| av免费在线看不卡| 高清毛片免费看| 国内揄拍国产精品人妻在线| 欧美3d第一页| av在线观看视频网站免费| 水蜜桃什么品种好| 国产黄片视频在线免费观看| 久久精品久久精品一区二区三区| 久久久久免费精品人妻一区二区| 最近视频中文字幕2019在线8| 国产在线一区二区三区精| 26uuu在线亚洲综合色| 国产一级毛片在线| 国产真实伦视频高清在线观看| 久久精品国产亚洲网站| 亚洲精品亚洲一区二区| 简卡轻食公司| 国产女主播在线喷水免费视频网站 | 亚洲人与动物交配视频| 精品一区二区三区人妻视频| 国产精品嫩草影院av在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲在久久综合| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 久久久精品94久久精品| 久久这里有精品视频免费| 亚洲18禁久久av| 国产成人精品一,二区| 精品午夜福利在线看| 又黄又爽又刺激的免费视频.| 国产成人免费观看mmmm| 亚洲天堂国产精品一区在线| 久久这里只有精品中国| 国产精品蜜桃在线观看| 黄片无遮挡物在线观看| 久久久久性生活片| 亚洲av中文av极速乱| 三级国产精品片| 免费高清在线观看视频在线观看| 国产美女午夜福利| 亚洲激情五月婷婷啪啪| 久久精品夜色国产| 亚洲激情五月婷婷啪啪| 日韩 亚洲 欧美在线| 天天躁夜夜躁狠狠久久av| 亚洲成人久久爱视频| 亚洲精品日韩在线中文字幕| 国产人妻一区二区三区在| 午夜福利网站1000一区二区三区| 国产片特级美女逼逼视频| 国内少妇人妻偷人精品xxx网站| 国产淫片久久久久久久久| 国产永久视频网站| 最近2019中文字幕mv第一页| 在线观看av片永久免费下载| 人妻系列 视频| 国产真实伦视频高清在线观看| 亚洲av不卡在线观看| 国产 一区精品| 如何舔出高潮| 久久精品久久久久久噜噜老黄| 成人二区视频| 哪个播放器可以免费观看大片| 日韩一区二区三区影片| 国产真实伦视频高清在线观看| 干丝袜人妻中文字幕| 亚洲精品,欧美精品| 中文字幕制服av| 99久久精品热视频| 久久久a久久爽久久v久久| 国产国拍精品亚洲av在线观看| 久久久久久久久中文| 久久99热这里只频精品6学生| 99久久九九国产精品国产免费| 国产精品久久视频播放| 全区人妻精品视频| 免费观看的影片在线观看| 成人毛片60女人毛片免费| 韩国av在线不卡| 欧美日韩精品成人综合77777| 亚洲精品乱码久久久v下载方式| av免费在线看不卡| 亚洲三级黄色毛片| 97在线视频观看|