亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        信陽臺(tái)地磁觀測(cè)數(shù)據(jù)干擾分析及觀測(cè)質(zhì)量改善

        2021-12-16 07:35:38陳賢成萬里黃恩賢
        高原地震 2021年3期
        關(guān)鍵詞:信陽臺(tái)站臺(tái)階

        陳賢,成萬里,黃恩賢

        (河南省地震局,河南 鄭州 450000)

        0 引 言

        在任何一種地震觀測(cè)手段的觀測(cè)數(shù)據(jù)中,都不可避免的會(huì)出現(xiàn)各種干擾現(xiàn)象,例如自然干擾、環(huán)境干擾、人為干擾、儀器故障等干擾[1-2]。地磁觀測(cè)也是如此,對(duì)觀測(cè)數(shù)據(jù)產(chǎn)生的干擾種類繁多,這些干擾數(shù)據(jù)的曲線形態(tài)主要表現(xiàn)為記錄曲線的單點(diǎn)突跳、頻繁的單點(diǎn)或多點(diǎn)突跳、一定時(shí)段內(nèi)數(shù)據(jù)抬升或下降(臺(tái)階)。造成這種記錄曲線異常的主要原因有以下幾種:磁暴、高壓直流輸電、車輛、建設(shè)施工、城際輕軌、地鐵、背景噪聲、人為因素以及儀器內(nèi)部干擾等[3]。只有充分認(rèn)識(shí)地磁場(chǎng)的干擾因素,并將干擾因素進(jìn)行分類總結(jié),正確地做好預(yù)處理工作,保障好觀測(cè)資料質(zhì)量,才能更好地分析地磁場(chǎng)的異常變化,為地震監(jiān)測(cè)與預(yù)報(bào)提供真實(shí)準(zhǔn)確的信息支持[4-5]。

        基于地磁數(shù)據(jù)跟蹤分析軟件,以信陽地震臺(tái)(以下簡(jiǎn)稱“信陽臺(tái)”)近年來地磁觀測(cè)資料為例,分析近年來受到的各種干擾因素影響,并提供典型干擾數(shù)據(jù)圖例,總結(jié)干擾特征并給出合理的降低干擾影響的措施,以便提高地震臺(tái)站人員識(shí)別各種干擾的能力和數(shù)據(jù)預(yù)處理的水平。

        1 觀測(cè)場(chǎng)地背景

        信陽地震臺(tái)周邊觀測(cè)地質(zhì)構(gòu)造較為復(fù)雜,信陽臺(tái)地理坐標(biāo)為:32.13°N,114.02°E。位于信陽市正西方向8 km南灣轄區(qū),賢山北西向5°約1 km處,距南灣水庫(kù)約1.5 km。大地構(gòu)造上,該區(qū)隸屬秦嶺褶皺系,地震臺(tái)站附近發(fā)育有桐柏—商城斷裂帶、龜山—梅山斷裂帶等。地貌上,地震臺(tái)站位于桐柏山東麓低山丘陵區(qū)與山前沖積平原過渡帶。依據(jù)地貌類型,從西南到東北可劃為桐柏—大別山低山區(qū)、南灣—東雙河丘陵區(qū)、平昌關(guān)—信陽低丘陵平原區(qū)。本研究區(qū)域涉及華北斷塊、秦嶺斷褶系和揚(yáng)子斷塊區(qū)三個(gè)大地構(gòu)造單元。該區(qū)域在喜馬拉雅運(yùn)動(dòng)以前漫長(zhǎng)的地質(zhì)歷史中,主要處于近南北向的水平擠壓應(yīng)力場(chǎng)中,使區(qū)域構(gòu)造線大致呈北西西向展布。區(qū)域主要活動(dòng)構(gòu)造有:華北斷塊區(qū)南緣斷裂帶、欒川—魯山—漯河—太和斷裂帶、白河—谷城斷裂帶、青峰斷裂帶、新安—郟縣斷裂帶、許昌—鄢陵斷裂帶、蘭田—朱陽關(guān)—南陽—信陽斷裂帶、封門口—五指嶺斷裂帶。該區(qū)域地震構(gòu)造環(huán)境中對(duì)本區(qū)和鄰區(qū)有發(fā)震影響的主要構(gòu)造有:華北斷塊區(qū)南緣斷裂帶、蘭田—朱陽關(guān)—南陽—信陽斷裂帶、欒川—魯山—漯河—太和斷裂帶、丹鳳—西峽斷裂帶(圖1)。

        圖1 信陽臺(tái)附近地質(zhì)構(gòu)造

        信陽臺(tái)GM4地磁儀器位于臺(tái)站院內(nèi)磁房?jī)?nèi),GM4探頭距離臺(tái)站院墻外約為45 m,院墻外時(shí)常有過往車輛干擾、高壓直流輸電干擾等電磁干擾、施工人員在磁房附近施工、砍樹、挖樹等人為干擾。

        2 觀測(cè)資料干擾性分析

        利用地震前兆地磁數(shù)據(jù)處理軟件和地磁數(shù)據(jù)跟蹤分析軟件,通過分析2016年以來信陽臺(tái)GM4型磁通門磁力儀觀測(cè)數(shù)據(jù)受到多種干擾情況,認(rèn)為GM4地磁儀受多種干擾影響,其影響持續(xù)時(shí)間、數(shù)據(jù)變化幅值、干擾變化形態(tài)如表1)。將干擾類型分為內(nèi)觀測(cè)環(huán)境干擾和外觀測(cè)環(huán)境干擾。內(nèi)觀測(cè)環(huán)境干擾主要為觀測(cè)室內(nèi)溫濕度變化儀器本身背景噪聲的影響;外觀測(cè)環(huán)境干擾主要為高壓直流輸電和臺(tái)站附近車輛干擾;內(nèi)、外觀測(cè)環(huán)境干擾都能通過及時(shí)有效的措施來降低其干擾影響,從而改善觀測(cè)資料質(zhì)量。

        表1 信陽臺(tái)GM4地磁儀受多種干擾特征及干擾形態(tài)對(duì)比

        2.1 內(nèi)觀測(cè)環(huán)境干擾分析

        2017年12月,信陽臺(tái)工作人員通過數(shù)據(jù)跟蹤分析軟件對(duì)地磁GM4數(shù)據(jù)進(jìn)行多臺(tái)對(duì)比,發(fā)現(xiàn)Z分量數(shù)據(jù)自2017年11月13日儀器重設(shè)補(bǔ)償值開始,與省內(nèi)其它地磁臺(tái)站同類型儀器數(shù)據(jù)相比變化不同步,存在明顯反常現(xiàn)象。表現(xiàn)為與觀測(cè)室溫度變化呈現(xiàn)一定相關(guān)性,即溫度上升,觀測(cè)數(shù)據(jù)出現(xiàn)下降趨勢(shì);溫度下降,觀測(cè)數(shù)據(jù)呈現(xiàn)上升趨勢(shì);溫度恒定,觀測(cè)數(shù)據(jù)基本穩(wěn)定(圖2)。

        由圖2可以看出信陽臺(tái)GM4觀測(cè)數(shù)據(jù)Z分量隨溫度變化較為反常。為核實(shí)此反常變化,分別對(duì)GM4儀器主機(jī)、探頭、模擬盒、通信系統(tǒng)、供電線路、UPS、供電電壓、避雷設(shè)施等進(jìn)行了檢查。同時(shí)對(duì)觀測(cè)房房頂、墻裙作防水、隔熱、磁房?jī)?nèi)部設(shè)施進(jìn)行檢查,且臺(tái)站周邊也無任何施工干擾,無明顯異常,判定觀測(cè)數(shù)據(jù)Z分量隨溫度異常變化為觀測(cè)室溫濕度原因。信陽臺(tái)工作人員在GM4地磁房門口位置架設(shè)了探頭、模擬盒進(jìn)行對(duì)比觀測(cè),從同步對(duì)比觀測(cè)中可以直觀地看出,Z分量數(shù)據(jù)漂移與溫度變化相關(guān)性大。因此著重分析了兩個(gè)階段:(1)2~6月溫度顯著上升,數(shù)據(jù)出現(xiàn)漂移下降趨勢(shì),下降幅度約40 nT,兩者相關(guān)系數(shù)R=0.69(圖3a);(2)9~12月溫度顯著下降,數(shù)據(jù)出現(xiàn)漂移上升趨勢(shì),上升幅度約120 nT,兩者相關(guān)系數(shù)R=0.98(圖3b)。

        圖2 信陽臺(tái)GM4觀測(cè)數(shù)據(jù)隨溫度變化曲線

        圖3 信陽臺(tái)GM4觀測(cè)數(shù)據(jù)Z分量與溫度觀測(cè)曲線

        儀器故障對(duì)地磁觀測(cè)數(shù)據(jù)造成的干擾是不可避免的、也是突發(fā)性的。儀器工作時(shí)間比較久之后,儀器內(nèi)部會(huì)出現(xiàn)一些存儲(chǔ)軟件故障,對(duì)地磁觀測(cè)數(shù)據(jù)也會(huì)產(chǎn)生干擾,并且這種突發(fā)情況下產(chǎn)生的干擾難以明確,一般情況下只能采取排除法檢查故障原因。例如信陽臺(tái)FHD數(shù)采CF存儲(chǔ)卡故障,無法正常儲(chǔ)存數(shù)據(jù),導(dǎo)致原始數(shù)據(jù)被壓制,甚至出現(xiàn)錯(cuò)誤數(shù)據(jù),這種情況預(yù)處理只能將原始錯(cuò)誤數(shù)據(jù)刪除,氣溫、氣壓這些觀測(cè)環(huán)境突變也會(huì)對(duì)儀器的觀測(cè)有影響,甚至影響儀器內(nèi)部元件,導(dǎo)致對(duì)地磁觀測(cè)造成干擾。

        2.2 外觀測(cè)環(huán)境干擾特征

        高壓直流輸電干擾。在部分臺(tái)站的地磁觀測(cè)中,當(dāng)受到地電觀測(cè)、高頻干擾、地鐵干擾時(shí),觀測(cè)數(shù)據(jù)的表現(xiàn)形態(tài)為周期性規(guī)律的方波或大量毛刺[6]。這與高壓直流輸電干擾有著明顯的區(qū)別,仔細(xì)觀測(cè)數(shù)據(jù)變化形態(tài)就可以很好判定為非高壓直流輸電故障干擾。但當(dāng)?shù)卮庞^測(cè)受到臺(tái)站周圍的人為活動(dòng)、車輛運(yùn)行、基建工程等干擾時(shí),其干擾后的觀測(cè)數(shù)據(jù)變化形態(tài)表現(xiàn)大多為尖峰、緩變、臺(tái)階等,這和高壓直流輸電故障干擾表現(xiàn)形態(tài)大致相同,這時(shí)就需要我們?nèi)ヅ卸ǎ瑢⑷藶楫a(chǎn)生的干擾與高壓直流輸電干擾區(qū)分開來。判定區(qū)分的方法為:(1)通常人為、車輛運(yùn)行、基建工程等環(huán)境干擾時(shí),總場(chǎng)F、水平分量 H、垂直分量Z、偏角D會(huì)出時(shí)出現(xiàn)尖峰、緩變、臺(tái)階。而高壓直流輸電故障時(shí),主要是對(duì)地磁的垂直分量Z產(chǎn)生影響。當(dāng)對(duì)垂直分量Z影響較大時(shí),才會(huì)對(duì)水平分量H,偏角D產(chǎn)生一定影響,且垂直分量明顯大于水平分量。(2)人為、車輛運(yùn)行、基建工程等環(huán)境干擾只對(duì)單臺(tái)產(chǎn)生影響,而高壓直流輸電故障通常會(huì)對(duì)線路周邊的多個(gè)臺(tái)站產(chǎn)生影響,多個(gè)臺(tái)站在基本相同的時(shí)間內(nèi)出現(xiàn)相同的干擾形態(tài)。(3)位于輸電線路兩側(cè)的地磁觀測(cè)臺(tái)站,由于地磁觀測(cè)參數(shù)各方向分量的差異,垂直分量Z方向與受到干擾的臺(tái)階方向是相反的,同側(cè)臺(tái)站臺(tái)階方向則是相同的,比如在受到干擾時(shí),??h、洛陽、代縣臺(tái)階向上,信陽、武漢、應(yīng)城、蒙城臺(tái)階則是向下的。(4)高壓直流輸電故障干擾幅度,隨著臺(tái)站距線路距離變大而變小,并與不平衡電流的大小成比,同一條線路的不同臺(tái)站的干擾幅度成比值,此比值為常數(shù),比如信陽臺(tái)干擾幅度為潢川的臺(tái)的1.5倍,為蒙城臺(tái)的2倍減去0.1~ 0.2 nT。

        高壓直流輸電故障對(duì)地磁觀測(cè)產(chǎn)生影響時(shí),其形態(tài)主要表現(xiàn)為尖峰、緩變、臺(tái)階,臺(tái)階又可分為單個(gè)臺(tái)階、多重臺(tái)階、緩變臺(tái)階、跨日臺(tái)階。以信陽臺(tái)為例分別對(duì)上述形態(tài)做出描述(通過對(duì)比原始圖和預(yù)處理后的圖,可以更清晰地看出數(shù)據(jù)受干擾后的變化)。2020年12月3日高壓直流輸電干擾,由于臺(tái)階當(dāng)天沒有結(jié)束,所以稱之為跨天臺(tái)階。(1)尖峰:所謂尖峰即是單點(diǎn)脈沖,01時(shí)10分表現(xiàn)為尖峰。(2)單個(gè)臺(tái)階:突然變化且變化量大小相同的方波,05時(shí)49分至06時(shí)45分表現(xiàn)為單個(gè)臺(tái)階。(3)多重臺(tái)階:在一個(gè)變化的臺(tái)階中再次出現(xiàn)新的變化量的臺(tái)階,11時(shí)30分至14時(shí)10分表現(xiàn)為多重臺(tái)階。(4)緩變臺(tái)階:在長(zhǎng)時(shí)間內(nèi)緩慢上升或下降的臺(tái)階,16時(shí)05分至18時(shí)58分的臺(tái)階中為緩慢上升所以稱其為緩變臺(tái)階。臺(tái)站附近車輛、施工干擾。

        2017年以來信陽臺(tái)觀測(cè)環(huán)境受周邊住戶居民過往車輛、防雷改造、優(yōu)化改造、標(biāo)準(zhǔn)化改造等干擾影響較為明顯,施工期間零線與接地線混為一團(tuán),造成觀測(cè)區(qū)附近有電流干擾,漏電電流引起Z分量觀測(cè)曲線出現(xiàn)密集型臺(tái)階變化,在全天記錄曲線上分布較為均勻,每間隔一定時(shí)間出現(xiàn)一組明顯臺(tái)階,持續(xù)時(shí)間為5 min,臺(tái)階幅度大小基本在0.5~2.5 nT。由于分布較為密集,磁通門磁力儀FGM01在記錄秒數(shù)據(jù)在預(yù)處理過程中十分繁瑣。當(dāng)磁房附近的車輛快速通過時(shí),這種干擾對(duì)磁房附近的磁場(chǎng)影響是短時(shí)的、突然性的,往往這種干擾對(duì)數(shù)據(jù)會(huì)產(chǎn)生比較明顯的突跳、尖峰;當(dāng)磁房附近停放車輛時(shí),這種干擾對(duì)磁房附近的磁場(chǎng)的影響會(huì)比較長(zhǎng),可以根據(jù)停放時(shí)間記錄下產(chǎn)生干擾的時(shí)間節(jié)點(diǎn),往往這種干擾對(duì)數(shù)據(jù)會(huì)產(chǎn)生明顯的臺(tái)階,臺(tái)階量的大小與停放車輛對(duì)大地磁場(chǎng)的影響量,臺(tái)階區(qū)間就是停放的時(shí)間節(jié)點(diǎn)。總之車輛干擾對(duì)地磁觀測(cè)分量干擾的幅度不盡相同,開始、結(jié)束時(shí)間相同,干擾時(shí)間與車輛通過或停放的時(shí)間具有一定的對(duì)應(yīng)性(圖4)。對(duì)于車輛干擾一般情況下需要對(duì)產(chǎn)生的突跳和尖峰做刪除預(yù)處理,對(duì)產(chǎn)生的臺(tái)階核實(shí)影響的時(shí)間節(jié)點(diǎn),然后做臺(tái)階預(yù)處理。

        圖4 信陽臺(tái)FGM01受高壓直流和車輛施工干擾記錄曲線(2020年12月14日)

        3 觀測(cè)資料干擾應(yīng)對(duì)措施分析

        針對(duì)上文所示內(nèi)觀測(cè)環(huán)境干擾和外觀測(cè)環(huán)境干擾,可采取以下措施,防止對(duì)地磁觀測(cè)產(chǎn)生的影響。

        (1)內(nèi)觀測(cè)環(huán)境干擾。此類干擾是偶然發(fā)生的,一般情況下當(dāng)儀器內(nèi)部元器件故障時(shí),觀測(cè)數(shù)據(jù)會(huì)產(chǎn)生缺數(shù)或者突跳,在分?jǐn)?shù)據(jù)和秒數(shù)據(jù)同時(shí)出現(xiàn);而當(dāng)觀測(cè)室溫濕度變化幅度較大時(shí),觀測(cè)數(shù)據(jù)會(huì)出現(xiàn)漂移等現(xiàn)象。針對(duì)此類干擾需地震臺(tái)站工作人員定期對(duì)儀器工作環(huán)境和工作狀態(tài)進(jìn)行周、月、季檢查。如每周進(jìn)行一次周檢查,每月進(jìn)行一次月檢查,每個(gè)季度進(jìn)行一次大范圍內(nèi)巡檢,檢查內(nèi)容包括更換觀測(cè)室干燥劑、除濕劑等。每次檢查后做好登記簿,觀測(cè)儀器工作環(huán)境和工作狀態(tài)情況,保持在觀測(cè)室溫濕度維持一定范圍內(nèi),檢查儀器運(yùn)行情況。在進(jìn)行臺(tái)站標(biāo)準(zhǔn)化、優(yōu)化改造等大型施工后,也需及時(shí)檢查儀器連接線與接地線的連接方式,以免儀器連接線與接地線混接,導(dǎo)致產(chǎn)生一些干擾電流,影響儀器觀測(cè)資料質(zhì)量??傊疁p少內(nèi)觀測(cè)環(huán)境背景噪聲,才能降低儀器受工作環(huán)境和工作狀態(tài)的影響。

        (2)外觀測(cè)環(huán)境干擾。此類干擾是偶然發(fā)生的,一般情況下當(dāng)儀器附件有車輛通過或者人員施工時(shí),觀測(cè)數(shù)據(jù)會(huì)產(chǎn)生臺(tái)階或者突跳,在分鐘值數(shù)據(jù)和秒數(shù)據(jù)同時(shí)出現(xiàn);而當(dāng)高壓直流輸電對(duì)觀測(cè)數(shù)據(jù)產(chǎn)生影響時(shí),觀測(cè)數(shù)據(jù)會(huì)根據(jù)輸電影響時(shí)間對(duì)應(yīng)出現(xiàn)規(guī)則的臺(tái)階或緩變臺(tái)階現(xiàn)象。針對(duì)此類干擾因觀測(cè)區(qū)在一定時(shí)間內(nèi)收到不規(guī)律電磁干擾,主要包括高壓直流輸電干擾和觀測(cè)室周邊過往車輛干擾。對(duì)于高壓直流干擾,地震臺(tái)站工作人員首先要對(duì)比其他臺(tái)站或者查詢高壓直流網(wǎng)站分清是否為高壓直流干擾。高壓直流干擾處理觀測(cè)數(shù)據(jù)的方法主要過程如下:①剔除尖峰選用數(shù)據(jù)處理軟件中的“剔除尖峰”;②去除臺(tái)階選用數(shù)據(jù)處理軟件中的“去除臺(tái)階―錄入時(shí)間和改正量”;③對(duì)于緩慢變化時(shí)間較短的變化數(shù)據(jù),也可以用“剔除尖峰”方法處理;④緩慢變化時(shí)間較長(zhǎng)的,選用數(shù)據(jù)處理軟件中的“去除緩變―錄入時(shí)間和改正量”。對(duì)于觀測(cè)室周邊過往車輛干擾,應(yīng)在地震臺(tái)站觀測(cè)室附近增加警示標(biāo)志,劃定地磁觀測(cè)區(qū)范圍,盡量避免過往車輛、攜帶鐵磁性物質(zhì)的人員靠近觀測(cè)室,以降低周邊環(huán)境對(duì)觀測(cè)資料質(zhì)量造成的影響。

        4 結(jié) 論

        基于地磁數(shù)據(jù)跟蹤分析軟件對(duì)信陽臺(tái)地磁觀測(cè)資料受內(nèi)外觀測(cè)環(huán)境干擾因素分析,可以得出以下結(jié)論:

        (1)信陽臺(tái)地磁觀測(cè)儀器受觀測(cè)室溫濕度變化影響較大,而且觀測(cè)數(shù)據(jù)變現(xiàn)為較長(zhǎng)期存在,在日常觀測(cè)中要時(shí)刻保持觀測(cè)室溫濕度保持在一定范圍內(nèi)。

        (2)地磁儀器靈敏度較高,對(duì)于觀測(cè)時(shí)間較長(zhǎng)的地磁儀器,要長(zhǎng)期檢查主備機(jī)工作狀態(tài),特別是地磁數(shù)據(jù)存儲(chǔ)卡CF卡、SD卡,需保證存儲(chǔ)卡對(duì)數(shù)據(jù)采集狀態(tài)。

        (3)信陽臺(tái)地磁觀測(cè)環(huán)境受國(guó)內(nèi)多條高壓直流輸電干擾影響,地震臺(tái)站工作人員需在日常數(shù)據(jù)預(yù)處理過程中,核對(duì)每日觀測(cè)數(shù)據(jù)受高壓直流影響情況,準(zhǔn)確判斷出受干擾情況,正常預(yù)處理受高壓直流影響的地磁資料。

        (4)隨著中國(guó)經(jīng)濟(jì)全面快速的發(fā)展,對(duì)于周邊過往車輛干擾,必須在觀測(cè)室圍墻懸掛地震觀測(cè)標(biāo)志,盡量避免車輛距離觀測(cè)室30 m。受到車輛不同形態(tài)的干擾,采用多種方法進(jìn)行觀測(cè)數(shù)據(jù)的預(yù)處理,利用多臺(tái)對(duì)比、差值及相關(guān)的方法進(jìn)行檢驗(yàn),從而保證地磁觀測(cè)數(shù)據(jù)的真實(shí)、可靠、連續(xù)、穩(wěn)定,更好地服務(wù)于地震預(yù)報(bào)實(shí)踐和科學(xué)研究。

        猜你喜歡
        信陽臺(tái)站臺(tái)階
        中國(guó)科學(xué)院野外臺(tái)站檔案工作回顧
        氣象基層臺(tái)站建設(shè)
        西藏科技(2021年12期)2022-01-17 08:46:38
        戰(zhàn)“疫”大考中的信陽答卷
        繡繡信陽八大景
        青年歌聲(2018年8期)2018-10-22 01:16:20
        繡繡信陽八大景
        青年歌聲(2018年2期)2018-10-20 02:02:56
        走在除法的臺(tái)階上
        信陽茶魂
        電影故事(2016年5期)2016-06-15 20:27:30
        基層臺(tái)站綜合觀測(cè)業(yè)務(wù)管理之我見
        西藏科技(2015年6期)2015-09-26 12:12:13
        臺(tái)階
        小說月刊(2015年11期)2015-04-23 08:47:31
        77級(jí)臺(tái)階
        一本加勒比hezyo无码视频| 亚洲精品无码久久久久久| 偷看农村妇女牲交| 国产性生交xxxxx免费| 亚洲三级香港三级久久| 亚洲色图视频在线播放| 91久久精品一二三区色| 国产丝袜爆操在线观看| 久久精品中文字幕女同免费| 久久婷婷五月综合色丁香| 无人高清电视剧在线观看| 三级网址在线| 精品国产亚欧无码久久久| 国产精品黄色在线观看| 99国产精品久久一区二区三区 | 国产呦精品系列在线播放| 亚洲国产精品500在线观看| 国产精品天堂在线观看| 夜夜高潮夜夜爽夜夜爱爱一区 | 大岛优香中文av在线字幕| 国产自拍视频免费在线观看| 亚洲高清在线天堂精品| 艳妇臀荡乳欲伦69调教视频| 天美传媒一区二区| 正在播放一区| 亚洲免费成年女性毛视频| 在线观看二区视频网站二区 | 久久精品国产亚洲av麻豆色欲 | 人妻精品无码一区二区三区| 一区二区三无码| 国产精品自产拍av在线| 亚洲综合网国产精品一区| 狼人香蕉香蕉在线28 - 百度| 天堂а√在线最新版中文| 日韩欧美国产亚洲中文| 亚洲精品国产熟女久久久| 国产成人精品亚洲日本在线观看| 曰本女人与公拘交酡免费视频| 91国在线啪精品一区| 亚洲国产91精品一区二区| 久久亚洲av成人无码电影 |