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

    地面微地震有源噪聲自動(dòng)識(shí)別與匹配壓制方法

    2018-04-09 05:53:44吳國(guó)忱尚新民芮擁軍崔慶輝
    石油地球物理勘探 2018年2期
    關(guān)鍵詞:時(shí)窗傳播速度噪聲源

    刁 瑞 吳國(guó)忱 尚新民 芮擁軍 崔慶輝

    (①中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島 266555; ②中國(guó)石化勝利油田分公司物探研究院,山東東營(yíng) 257022)

    1 引言

    由于頁巖氣、頁巖油和煤層氣等非常規(guī)油氣儲(chǔ)層具有低孔隙度、低滲透率的物性特征,決定了其具有較大的開采難度[1,2]。非常規(guī)儲(chǔ)層水力壓裂改造技術(shù)成為提高油氣采收率的主要手段,特別是針對(duì)頁巖氣、頁巖油和煤層氣等非常規(guī)油氣資源,致密儲(chǔ)層的水力壓裂改造具有十分重要的作用和意義,微地震監(jiān)測(cè)技術(shù)則被用于評(píng)價(jià)壓裂效果[3,4]。地面陣列式微地震監(jiān)測(cè)技術(shù)采用多條測(cè)線數(shù)千個(gè)檢波器排列。基于星型排列的地面微地震監(jiān)測(cè)采用超過1000道、6000~24000個(gè)檢波器,測(cè)線排列長(zhǎng)度為2~10km。在地表進(jìn)行微地震監(jiān)測(cè),具有監(jiān)測(cè)范圍廣、數(shù)據(jù)采集量大等特點(diǎn)。相對(duì)井中微地震監(jiān)測(cè)技術(shù)而言,地面微地震數(shù)據(jù)具有噪聲干擾嚴(yán)重、信噪比低、事件識(shí)別難度大的特點(diǎn)。劉玉海等[5]利用相鄰道信號(hào)之間的相關(guān)性,提出了一種基于隨機(jī)信號(hào)統(tǒng)計(jì)理論的互相關(guān)檢測(cè)方法。張旭亮等[6]利用K-L變換從多道微地震記錄中提取微地震信號(hào),去除隨機(jī)噪聲和相干噪聲。趙翠霞等[7]基于正余弦加權(quán)逼近法實(shí)現(xiàn)了工頻干擾的自動(dòng)識(shí)別與壓制,有效提高了地震資料的品質(zhì)。刁瑞等[8,9]將改進(jìn)S變換引入微地震監(jiān)測(cè)數(shù)據(jù)的去噪處理,有效分離了微地震信號(hào)分量與噪聲干擾分量; 并基于陣列式微地震數(shù)據(jù)的優(yōu)點(diǎn),利用互相關(guān)的盲源分離去噪方法壓制隨機(jī)噪聲干擾。賈瑞生等[10]基于經(jīng)驗(yàn)?zāi)B(tài)分解及獨(dú)立成分分析提出了一種微地震信號(hào)降噪方法。李會(huì)儉等[11]將多尺度形態(tài)學(xué)理論應(yīng)用于弱信號(hào)分析和識(shí)別,通過分析不同尺度下的信號(hào)特征,估計(jì)并檢測(cè)出微弱信號(hào)和噪聲。宋維琪等[12]研究了地面微地震資料τ-p變換去噪方法,增強(qiáng)了地面微地震事件的辨識(shí)度。崔慶輝等[13]針對(duì)微地震數(shù)據(jù)信噪比低的特點(diǎn),聯(lián)合應(yīng)用多種噪聲壓制技術(shù),明顯提高了微地震信號(hào)的信噪比。

    在進(jìn)行地面微地震監(jiān)測(cè)采集過程中,由于地表各種不利因素的影響,微地震數(shù)據(jù)中強(qiáng)能量的地面有源噪聲非常多,包括井場(chǎng)噪聲、鉆機(jī)噪聲、建筑工地噪聲、車輛噪聲、風(fēng)噪、人步行和物體墜落噪聲等。壓裂井場(chǎng)噪聲以典型的面波方式傳播,向兩側(cè)逐漸衰減。車輛噪聲是典型的寬頻瞬變?cè)肼?,噪聲振幅具有由小變大、再變小的特征。人步行或物體墜落是典型的窄頻穩(wěn)定噪聲源,人步行引起的噪聲幾乎為固頻信號(hào),形成單一頻率的尖脈沖[14]。儲(chǔ)層壓裂改造產(chǎn)生的震源能量相對(duì)較弱,各類地面有源噪聲的能量強(qiáng),嚴(yán)重影響微地震數(shù)據(jù)的信噪比,微地震事件湮沒在噪聲中,增加了微地震事件的識(shí)別難度。必須通過針對(duì)性的地面有源噪聲壓制處理手段,提高微地震數(shù)據(jù)的信噪比,以保證微地震震源點(diǎn)的定位精度。在充分分析微地震噪聲特征基礎(chǔ)上,根據(jù)地面有源噪聲與微地震事件在能量、頻率、傳播速度及源位置不同的特點(diǎn),提出了一種地面微地震有源噪聲自動(dòng)識(shí)別與匹配壓制方法。在長(zhǎng)短時(shí)窗能量比和微地震量板聯(lián)合自動(dòng)識(shí)別有源噪聲的基礎(chǔ)上,進(jìn)行噪聲源位置和傳播速度三維最優(yōu)并行搜索,確定有源噪聲的位置坐標(biāo)和噪聲傳播速度后,根據(jù)噪聲標(biāo)準(zhǔn)道和自適應(yīng)匹配算子,對(duì)地面有源噪聲進(jìn)行自適應(yīng)壓制處理。

    2 有源噪聲特征分析

    對(duì)非常規(guī)油氣藏進(jìn)行儲(chǔ)層壓裂改造,水力壓裂造成的巖石破裂類似于一個(gè)震源,而震源能量相對(duì)較弱,微地震信號(hào)從震源點(diǎn)向上傳播過程中,能量又被地層強(qiáng)烈吸收衰減,造成微地震事件湮沒在噪聲干擾里。只有通過分析地面有源噪聲的特征,才能針對(duì)性地壓制噪聲,達(dá)到提高微地震數(shù)據(jù)品質(zhì)的目的。微地震事件與地面有源噪聲的特征差別主要包括以下幾個(gè)方面。

    (1)源位置差異。壓裂段位于地下儲(chǔ)層中,微地震事件由震源點(diǎn)從下向上傳播;而噪聲源位于地面,噪聲干擾沿地面或高速層頂界面?zhèn)鞑?。由于微地震事件和噪聲干擾的傳播路徑不同,到達(dá)不同站點(diǎn)的時(shí)差存在差異,微地震監(jiān)測(cè)數(shù)據(jù)中微地震事件的同相軸形態(tài)與噪聲干擾的同相軸形態(tài)不同,即不同站點(diǎn)之間的旅行時(shí)差存在特征差異。

    (2)傳播速度差異。壓裂造成的震源點(diǎn)由深層往地面?zhèn)鞑?,地下介質(zhì)的傳播速度較高,而地面噪聲源主要是通過地面或沿高速層頂界面滑行,傳播速度為低降速帶或高速層速度,巖石破裂地震波和有源噪聲的傳播速度存在明顯的差異,即微地震數(shù)據(jù)中的同相軸斜率或視速度不同。

    (3)能量差異。從整體上分析,微地震事件和地面有源噪聲的能量差異明顯,儲(chǔ)層壓裂造成的巖石破裂能量相對(duì)較弱,地面有源噪聲的能量相對(duì)較強(qiáng)。另外,地面站點(diǎn)在壓裂段上方的地面均勻分布,壓裂震源點(diǎn)到達(dá)地面站點(diǎn)時(shí)能量基本相當(dāng)或差別不大。地面有源噪聲在傳播過程中衰減劇烈,距離較近的地面站點(diǎn)中噪聲能量較強(qiáng),而距離較遠(yuǎn)的地面站點(diǎn)中噪聲能量較弱,即地面有源噪聲在不同距離的站點(diǎn)中能量差異明顯。

    (4)頻率差異。由于檢波器布設(shè)于地面,地震波經(jīng)歷了地層的吸收衰減作用,地面微地震監(jiān)測(cè)數(shù)據(jù)的頻率范圍是5~100Hz。車輛、建筑工地等噪聲具有寬頻特征,人步行、鉆機(jī)干擾、工頻干擾等具有窄頻或單頻特征,其中工頻干擾的頻率約為50Hz。

    圖1為地震波傳播路徑及模擬旅行時(shí)記錄,可見震源點(diǎn)位于地下儲(chǔ)層附近,微地震信號(hào)從下向上傳播,地面站點(diǎn)均勻分布,微地震事件具有雙曲線特征,并且地層速度相對(duì)較高,站點(diǎn)之間的旅行時(shí)差相對(duì)較小。噪聲源位于地面,有源噪聲沿高速層頂界面?zhèn)鞑ィ?dāng)噪聲源位于所有站點(diǎn)一側(cè)時(shí),地面有源噪聲具有線性特征(圖1b); 當(dāng)噪聲源位于站點(diǎn)中部時(shí),地面有源噪聲具有雙曲線特征,并多次重復(fù)出現(xiàn)(圖1c)。由于正演模擬過程中沒有考慮吸收衰減,因此不存在能量差異,但從旅行時(shí)和傳播速度兩個(gè)方面的差異特征,可有效識(shí)別微地震事件和地面有源噪聲。

    3 方法原理

    根據(jù)微地震事件與地面有源噪聲在旅行時(shí)、傳播速度、能量和頻率方面的差異,提出了微地震有源噪聲自動(dòng)識(shí)別與匹配壓制方法,對(duì)微地震數(shù)據(jù)中的有源噪聲進(jìn)行自適應(yīng)壓制。該技術(shù)的主要思路為:首先,在微地震數(shù)據(jù)中地面有源噪聲發(fā)育的情況下,采用長(zhǎng)短時(shí)窗能量比方法和微地震量板自動(dòng)識(shí)別有源噪聲,得到微地震數(shù)據(jù)中有源噪聲的到達(dá)時(shí)刻;其次,根據(jù)疊加能量最強(qiáng)的原則,沿地面進(jìn)行三維最優(yōu)并行搜索,同時(shí)確定有源噪聲的位置坐標(biāo)和噪聲傳播速度;然后,根據(jù)有源噪聲的位置坐標(biāo)和傳播速度,計(jì)算每一個(gè)微地震站點(diǎn)的動(dòng)校正量,并對(duì)微地震數(shù)據(jù)進(jìn)行動(dòng)校正處理;最后,根據(jù)有源噪聲標(biāo)準(zhǔn)道和自適應(yīng)匹配算子,對(duì)有源噪聲進(jìn)行自適應(yīng)壓制處理,并進(jìn)行反動(dòng)校正處理,可得到高信噪比的微地震數(shù)據(jù)。該技術(shù)主要包括三方面內(nèi)容: ①長(zhǎng)短時(shí)窗能量比與微地震量板聯(lián)合的有源噪聲自動(dòng)識(shí)別; ②源位置和傳播速度三維最優(yōu)并行搜索; ③基于匹配算子的有源噪聲自適應(yīng)壓制。

    3.1 長(zhǎng)短時(shí)窗能量比與微地震量板聯(lián)合自動(dòng)識(shí)別

    微地震數(shù)據(jù)初至拾取方法主要有長(zhǎng)短時(shí)窗能量比方法(STA/LTA)、偏振分析及AIC方法、初至拾取SLPEA算法、基于匹配追蹤的事件識(shí)別方法等[15,16]。劉勁松等[17]利用改進(jìn)的高階累計(jì)量的偏度和峰度計(jì)算方法識(shí)別微地震事件。王鵬等[18]利用Renyi熵值表示微地震數(shù)據(jù)的時(shí)頻稀疏程度,建立以低熵值的道數(shù)為判別閾值的目標(biāo)函數(shù),在對(duì)低信噪比數(shù)據(jù)處理中取得了較好效果。長(zhǎng)短時(shí)窗能量比方法根據(jù)環(huán)境噪聲振幅弱、頻帶寬,而地震信號(hào)振幅強(qiáng)、頻帶窄的特點(diǎn),通過計(jì)算一長(zhǎng)一短兩個(gè)滑動(dòng)時(shí)窗內(nèi)地震記錄特征函數(shù)的平均值之比作為拾取地震波初至到時(shí)的依據(jù)。短時(shí)窗能量值的變化往往比長(zhǎng)時(shí)窗能量值的變化快,用短時(shí)窗能量平均值與長(zhǎng)時(shí)窗能量平均值進(jìn)行對(duì)比,當(dāng)存在能量變化的情況下,長(zhǎng)短時(shí)窗能量比值就會(huì)明顯增大。

    僅僅依賴地震信號(hào)與環(huán)境噪聲的差異進(jìn)行初至拾取是比較困難的,長(zhǎng)短時(shí)窗能量比方法無法區(qū)分微地震事件和噪聲干擾,應(yīng)該綜合考慮多種特征差異自動(dòng)識(shí)別噪聲干擾。在此基于能量、傳播速度和源位置不同的特點(diǎn),采用長(zhǎng)短時(shí)窗能量比與微地震量板聯(lián)合的方法,自動(dòng)識(shí)別微地震數(shù)據(jù)中的有源噪聲。長(zhǎng)短時(shí)窗能量比方法從能量差異的方面進(jìn)行判別,微地震量板從旅行時(shí)差異或時(shí)距關(guān)系的方面進(jìn)行進(jìn)一步識(shí)別,最終實(shí)現(xiàn)有源噪聲的自動(dòng)識(shí)別。

    相對(duì)于微地震事件而言,地面有源噪聲能量相對(duì)較強(qiáng),可利用長(zhǎng)短時(shí)窗能量比方法確定微地震數(shù)據(jù)中噪聲源分布范圍。在利用能量比法進(jìn)行計(jì)算時(shí),由于初至點(diǎn)以上少數(shù)樣點(diǎn)的幅值有可能接近零,為了避免這種情況的發(fā)生,同時(shí)提高初至拾取的穩(wěn)定性,對(duì)前后時(shí)窗的能量分別加一個(gè)穩(wěn)定因子。穩(wěn)定能量比方法的數(shù)學(xué)表達(dá)式為

    (1)

    式中:f(t)為微地震數(shù)據(jù);T1為短時(shí)窗開始時(shí)刻;T0為短時(shí)窗與長(zhǎng)時(shí)窗交叉時(shí)刻;T2為長(zhǎng)時(shí)窗結(jié)束時(shí)刻;A為振幅因子;α為穩(wěn)定系數(shù)(可據(jù)不同數(shù)據(jù)調(diào)整)。設(shè)置一個(gè)自動(dòng)識(shí)別能量閾值θ,當(dāng)長(zhǎng)短時(shí)窗能量比值大于θ時(shí),則該時(shí)刻就存在有源噪聲,當(dāng)長(zhǎng)短時(shí)窗能量比值小于θ時(shí),則不存在有源噪聲。

    長(zhǎng)短時(shí)窗能量比方法識(shí)別的信號(hào)中可能包含微地震事件,需要通過微地震量板做進(jìn)一步識(shí)別。微地震量板即為震源點(diǎn)到所有接收點(diǎn)的走時(shí)形態(tài),儲(chǔ)層壓裂產(chǎn)生的震源點(diǎn)位于射孔點(diǎn)周圍,以射孔點(diǎn)為震源點(diǎn),在建立速度模型和已知站點(diǎn)分布的基礎(chǔ)上,通過正演模擬方法獲得震源點(diǎn)的正演記錄,其中壓裂震源點(diǎn)的走時(shí)形態(tài)與微地震量板基本一致。計(jì)算識(shí)別信號(hào)與微地震量板之間的均方差值φ,通過差值φ與微地震量板因子δ識(shí)別信號(hào),其數(shù)學(xué)表達(dá)式為

    (2)

    Y-8HF井鉆遇砂礫巖油藏,位于東營(yíng)凹陷北部陡坡帶東翼,主力含油層系為沙四段砂礫巖體,是常溫、常壓、低孔隙度、特低滲油藏。對(duì)砂礫巖體致密油藏進(jìn)行儲(chǔ)層壓裂改造,Y-8HF井共13段壓裂,壓裂方向近為南北向,微地震監(jiān)測(cè)站點(diǎn)位于壓裂段上方。圖2a為兩個(gè)壓裂段的微地震量板,是根據(jù)兩個(gè)壓裂段的射孔點(diǎn)分別計(jì)算得到的,藍(lán)色線是第2壓裂段,紅色線是第5壓裂段; 圖2b為微地震量板與實(shí)際監(jiān)測(cè)數(shù)據(jù)的對(duì)比,微地震量板與實(shí)際數(shù)據(jù)的旅行時(shí)形態(tài)基本一致(即φ<δ),該識(shí)別信號(hào)在微地震事件集合內(nèi),最終確定為微地震事件,避免了微地震事件的誤拾和有源噪聲的漏拾。

    圖2 不同壓裂段微地震量板(a)及與實(shí)際地震數(shù)據(jù)(b)的對(duì)比

    3.2 源位置和傳播速度三維最優(yōu)并行搜索

    噪聲源分布于地面或近地表,因此可以沿地面的大地坐標(biāo)X和Y方向進(jìn)行噪聲源位置的搜索;有源噪聲沿高速頂界面?zhèn)鞑ィ梢源_定傳播速度v的范圍?;谝陨蟽煞矫娴囊蛩兀梢赃M(jìn)行噪聲源位置(X和Y)和噪聲傳播速度(v)的同時(shí)最優(yōu)搜索,并引入并行搜素機(jī)制,大幅提高運(yùn)算效率。

    首先,利用不同噪聲源位置和傳播速度對(duì)自動(dòng)識(shí)別的有源噪聲做動(dòng)校正處理,消除傳播距離造成的時(shí)差;其次,計(jì)算動(dòng)校正后數(shù)據(jù)疊加能量,可得三維疊加能量譜,計(jì)算中采用并行算法可大幅提高效率;最后,求取疊加能量譜最大值,該最大值對(duì)應(yīng)的位置和傳播速度即為最優(yōu)搜索結(jié)果。

    沿大地坐標(biāo)X方向平均剖分為M個(gè)網(wǎng)格,用i表示網(wǎng)格序號(hào)(i=0,1,2,…,M); 沿大地坐標(biāo)Y方向?qū)⒌乇砥骄史譃镹個(gè)網(wǎng)格,用j表示網(wǎng)格序號(hào)(j=0,1,2,…,N); 將速度掃描范圍均分為Z個(gè)速度值,用k表示速度掃描序號(hào)(k=0,1,2,…,Z)。最優(yōu)并行搜素疊加能量數(shù)學(xué)表達(dá)式為

    (3)

    式中:S為地面檢波器站點(diǎn)總數(shù);fp(t)為第p個(gè)站點(diǎn)數(shù)據(jù);li,j,p為第i×j個(gè)網(wǎng)格點(diǎn)到第p個(gè)站點(diǎn)傳播距離;li,j,min為第i×j個(gè)網(wǎng)格點(diǎn)到所有站點(diǎn)最小距離;vk為第k個(gè)掃描速度值。對(duì)動(dòng)校正后的數(shù)據(jù)求取疊加能量,當(dāng)Ei,j,k值最大時(shí),所對(duì)應(yīng)的第i×j個(gè)網(wǎng)格點(diǎn)和掃描速度值vk即為最優(yōu)的噪聲源位置和傳播速度,最優(yōu)噪聲源網(wǎng)格點(diǎn)記為ib×jb,最優(yōu)噪聲傳播速度記為vbest。

    3.3 基于匹配算子的有源噪聲自適應(yīng)壓制

    利用確定的最優(yōu)噪聲源位置和傳播速度對(duì)微地震數(shù)據(jù)進(jìn)行動(dòng)校正處理,動(dòng)校正后有源噪聲同相軸表現(xiàn)為水平同相軸[19-21],將各道數(shù)據(jù)疊加取均值可得到噪聲標(biāo)準(zhǔn)道[22],在引入自適應(yīng)匹配算子的基礎(chǔ)上,對(duì)所有站點(diǎn)的有源噪聲進(jìn)行自適應(yīng)壓制。

    有源噪聲標(biāo)準(zhǔn)道n(t)的數(shù)學(xué)表達(dá)式為

    (4)

    式中:lib,jb,p為噪聲源所在網(wǎng)格到第p個(gè)站點(diǎn)的距離;lib,jb,min為噪聲源所在網(wǎng)格到所有站點(diǎn)中的最小距離;vbest為噪聲傳播速度。

    以有源噪聲標(biāo)準(zhǔn)道n(t)為基準(zhǔn),通過維納方程和Levinson算法進(jìn)行自適應(yīng)匹配,求得第p個(gè)站點(diǎn)的自適應(yīng)匹配算子hp(t),然后從微地震數(shù)據(jù)中減掉有源噪聲,得到噪聲壓制后數(shù)據(jù)

    (5)

    4 模型試算及資料試處理

    應(yīng)用微地震噪聲壓制方法,能有效提高微地震數(shù)據(jù)的品質(zhì)及微地震事件的定位精度[23]。 通過正演模擬數(shù)據(jù)驗(yàn)證方法的可行性,在模擬過程中采用Y-8HF井的實(shí)際微地震監(jiān)測(cè)觀測(cè)系統(tǒng),該井儲(chǔ)層為致密砂礫巖體[24]。圖3a為Y-8HF井的地面微地震監(jiān)測(cè)站點(diǎn)分布示意圖,藍(lán)色線為地面臺(tái)站分布,共計(jì)約900個(gè)站點(diǎn),南北方向分布范圍是1700m,東西方向分布范圍600m,黃色線為井軌跡,水平段位于靶點(diǎn)A和靶點(diǎn)B之間,由南向北共計(jì)13段壓裂;圖3b為地面噪聲源分布示意圖,噪聲源位于地面站點(diǎn)的東南方,坐標(biāo)為(5000m,1000m),距離檢波器站點(diǎn)約為2~3km。

    根據(jù)壓裂井區(qū)的地層特征,高速層速度采用2000m/s,對(duì)噪聲源的傳播特征進(jìn)行模擬,得到有源噪聲正演模擬結(jié)果(圖4a);將正演模擬數(shù)據(jù)與實(shí)際微地震數(shù)據(jù)(圖4b)進(jìn)行對(duì)比,可見兩者的噪聲具有相同特征。由于噪聲源位于地面站點(diǎn)的東南方,且距離較遠(yuǎn),地面有源噪聲具有線性特征,并且噪聲源具有一定的規(guī)律性,有源噪聲在微地震數(shù)據(jù)中重復(fù)出現(xiàn),間隔時(shí)間約為500~600ms。

    基于所有微地震站點(diǎn)數(shù)據(jù)的疊加能量最大原則, 沿地面有源噪聲分布范圍進(jìn)行三維最優(yōu)并行搜索,在搜索噪聲傳播速度的同時(shí),計(jì)算有源噪聲分布范圍內(nèi)每一個(gè)位置點(diǎn)到達(dá)微地震站點(diǎn)的波至?xí)r間,根據(jù)對(duì)應(yīng)的波至?xí)r間對(duì)微地震數(shù)據(jù)進(jìn)行動(dòng)校正處理,并將動(dòng)校正處理后的所有道數(shù)據(jù)疊加,以求取每一個(gè)位置點(diǎn)的微地震數(shù)據(jù)疊加能量值,圖5為三維最優(yōu)并行搜索能量譜示意圖。

    圖3 地面微地震監(jiān)測(cè)站點(diǎn)分布(a)及地面噪聲源分布(b)示意圖

    圖4 有源噪聲正演模擬數(shù)據(jù)(a)與實(shí)際微地震數(shù)據(jù)(b)對(duì)比

    當(dāng)某點(diǎn)的疊加能量值達(dá)到最大時(shí),說明該點(diǎn)和對(duì)應(yīng)的噪聲傳播速度能最合理地消除由于傳播距離不同造成的正常時(shí)差,實(shí)現(xiàn)最優(yōu)動(dòng)校正處理,該點(diǎn)即是有源噪聲的位置, 其對(duì)應(yīng)的傳播速度即是有源噪聲的傳播速度。圖6a為東西方向能量譜,圖6b為南北方向能量譜,圖6c為速度方向能量譜,通過三維最優(yōu)搜索得到噪聲源位置為(5000m,1000m),噪聲傳播速度為2000m/s,與模型數(shù)據(jù)一致,驗(yàn)證了該方法的有效性。速度方向能量譜的聚焦性較差,顯示有源噪聲的疑似位置呈西北—東南方向分布,這是由于觀測(cè)系統(tǒng)的南北方向排列相對(duì)較長(zhǎng),而東西方向排列相對(duì)較短,并且該有源噪聲位于排列的東南方向,距離較遠(yuǎn),因此能量譜聚焦性相對(duì)較差。在進(jìn)行觀測(cè)系統(tǒng)設(shè)計(jì)時(shí),應(yīng)盡量將檢波器分散布設(shè),以使各個(gè)方位角均勻分布,利于提高定位精度。

    圖5 三維最優(yōu)并行搜索能量譜

    根據(jù)地面有源噪聲位置坐標(biāo)、噪聲傳播速度及微地震監(jiān)測(cè)觀測(cè)系統(tǒng)坐標(biāo),計(jì)算每一個(gè)微地震站點(diǎn)的動(dòng)校正量。將動(dòng)校正量代入微地震數(shù)據(jù)中,對(duì)微地震數(shù)據(jù)進(jìn)行動(dòng)校正處理,消除由于不同微地震監(jiān)測(cè)檢波器之間傳播距離不同造成的正常時(shí)差。圖7為動(dòng)校正處理前、后正演模擬數(shù)據(jù),可見動(dòng)校正處理后地面有源噪聲同相軸基本是一條水平直線,所有道數(shù)據(jù)疊加能量即為最大值,說明三維最優(yōu)搜索確定的噪聲源位置和噪聲傳播速度是合理的。

    為了進(jìn)一步驗(yàn)證方法的可行性和有效性,對(duì)實(shí)際監(jiān)測(cè)微地震數(shù)據(jù)進(jìn)行噪聲壓制處理。圖8a為噪聲壓制前微地震數(shù)據(jù), 其中包括正演模擬微地震事件和地面建筑工地有源噪聲,該噪聲源位于排列的一側(cè),有源噪聲具有線性特征; 圖8b為地面有源噪聲壓制后微地震數(shù)據(jù),從圖中可以看出,地面有源噪聲自適應(yīng)壓制技術(shù)可以有效地去除微地震數(shù)據(jù)中的地面有源噪聲,并不會(huì)損失有效信號(hào),噪聲壓制效果較好。圖9a為實(shí)際監(jiān)測(cè)微地震數(shù)據(jù),其中包含周期性的地面有源噪聲,該噪聲源位于排列的中部,有源噪聲具有雙曲特征,隨著傳播距離的增加,噪聲干擾明顯減弱,距離噪聲源較遠(yuǎn)的臺(tái)站幾乎接收不到噪聲干擾; 圖9b為噪聲干擾壓制后數(shù)據(jù),地面有源噪聲得到了較好的壓制,數(shù)據(jù)信噪比明顯提升,有利于后續(xù)的微地震事件識(shí)別和精確定位。

    圖6 東西方向(a)、南北方向(b)及速度方向(c)能量譜

    圖7 動(dòng)校正處理前(a)、后(b)的模擬數(shù)據(jù)對(duì)比

    圖8 模擬地面有源噪聲壓制前(a)、后(b)效果對(duì)比藍(lán)色箭頭指示微地震事件,紅色箭頭指示有源噪聲

    圖9 實(shí)際地面有源噪聲壓制前(a)、后(b)效果對(duì)比紅色箭頭指示有源噪聲

    4 結(jié)論與展望

    微地震數(shù)據(jù)中地面有源噪聲較多,而儲(chǔ)層壓裂改造產(chǎn)生的震源能量相對(duì)較弱,微地震事件湮沒于噪聲中。在充分分析微地震噪聲特征基礎(chǔ)上,根據(jù)地面有源噪聲與微地震事件在能量、頻率、傳播速度及源位置等方面的差異,提出了一種地面微地震有源噪聲自動(dòng)識(shí)別與匹配壓制方法,獲得了高信噪比微地震數(shù)據(jù)。正演模擬數(shù)據(jù)和實(shí)際微地震數(shù)據(jù)去噪效果驗(yàn)證了該方法的實(shí)用性和有效性。

    在微地震數(shù)據(jù)的有源噪聲自動(dòng)識(shí)別過程中,采用長(zhǎng)短時(shí)窗能量比與微地震量板相聯(lián)合的方法。雖然長(zhǎng)短時(shí)窗能量比方法無法確定強(qiáng)能量信號(hào)是有源噪聲還是微地震事件,但以微地震量板為準(zhǔn)則,進(jìn)行交互判斷和識(shí)別,避免了微地震事件的誤拾和有源噪聲的漏拾。通過三維最優(yōu)搜索同時(shí)確定地面噪聲源的坐標(biāo)位置和傳播速度,并引入并行搜索機(jī)制以大幅提高運(yùn)算效率。根據(jù)所獲地面噪聲源坐標(biāo)位置和傳播速度,在動(dòng)?;A(chǔ)上,通過自適應(yīng)匹配算子對(duì)有源噪聲進(jìn)行壓制處理,實(shí)現(xiàn)了有源噪聲的自適應(yīng)壓制。

    隨著頁巖油氣、煤層氣等非常規(guī)資源的不斷開發(fā)和儲(chǔ)層壓裂改造技術(shù)的應(yīng)用,對(duì)微地震監(jiān)測(cè)技術(shù)的需求越來越多。國(guó)內(nèi)微地震監(jiān)測(cè)技術(shù)也取得迅速發(fā)展,如高靈敏度監(jiān)測(cè)儀器研發(fā)、觀測(cè)系統(tǒng)針對(duì)性設(shè)計(jì)、震源機(jī)制反演、三維速度模型優(yōu)化校正、精確定位和壓裂裂縫綜合解釋等,但還需持續(xù)深入研究和探索。微地震監(jiān)測(cè)技術(shù)仍將是一個(gè)研究熱點(diǎn),特別是地面與井中聯(lián)合監(jiān)測(cè)、地面臺(tái)站長(zhǎng)期監(jiān)測(cè)、有纜地震與無纜節(jié)點(diǎn)聯(lián)合采集等技術(shù),具有廣闊的應(yīng)用前景。

    感謝中國(guó)石化勝利油田分公司物探研究院微地震項(xiàng)目組的鼎力支持!

    [1]劉振武,撒利明,楊曉等.頁巖氣勘探開發(fā)對(duì)地球物理技術(shù)的需求.石油地球物理勘探,2011,46(5):810-820.

    Liu Zhenwu,Sa Liming,Yang Xiao et al.Needs of geophysical technologies for shale gas exploration.OGP,2011,46(5):810-820.

    [2]刁瑞,吳國(guó)忱,尚新民等.三維地震與地面微地震聯(lián)合校正方法.地球物理學(xué)報(bào),2017,60(1):283-292.

    Diao Rui,Wu Guochen,Shang Xinmin et al.Joint correction method based on 3D seismic and surface microseismic data.Chinese Journal of Geophysics,2017,60(1):283-292.

    [3]孫龍德,撒利明,董世泰.中國(guó)未來油氣新領(lǐng)域與物探技術(shù)對(duì)策.石油地球物理勘探,2013,48(2):317-324.

    Sun Longde,Sa Liming,Dong Shitai.New challenges for the future hydrocarbon in China and geophysical technology strategy.OGP,2013,48(2):317-324.

    [4]刁瑞,吳國(guó)忱,崔慶輝等.地面陣列式微地震監(jiān)測(cè)關(guān)鍵技術(shù)研究.巖性油氣藏,2017,29(1):104-109.

    Diao Rui,Wu Guochen,Cui Qinghui et al.Key techniques for surface array microseismic monitoring.Lithologic Reservoirs,2017,29(1):104-109.

    [5]劉玉海,尹成,潘樹林等.基于互相關(guān)函數(shù)法的地面微地震信號(hào)檢測(cè)研究.石油物探,2012,51(6):633-637.

    Liu Yuhai,Yin Chen,Pan Shulin et al.Ground microseismic signal detection based on cross-correlation function.GPP,2012,51(6):633-637.

    [6]張旭亮,桂志先,王鵬等.基于K-L變換的微地震資料去噪的分析及應(yīng)用.工程地球物理學(xué)報(bào),2013,10(1):81-84.

    Zhang Xuliang,Gui Zhixian,Wang Peng et al.The analysis and application of micro-seismic data denoising based onK-Ltransform.Chinese Journal of Engineering Geophysics,2013,10(1):81-84.

    [7]趙翠霞,尚新民,刁瑞等.基于正余弦加權(quán)逼近法的50Hz工業(yè)干擾壓制方法.油氣地球物理,2015,13(4):14-17.

    Zhao Cuixia,Shang Xinmin,Diao Rui et al.Method to suppress 50Hz industrial noise based on sine and cosine weighted approximation.Petroleum Geophysics,2015,13(4):14-17.

    [8]刁瑞,單聯(lián)瑜,尚新民等.微地震監(jiān)測(cè)數(shù)據(jù)時(shí)頻域去噪方法.物探與化探,2015,39(1):112-117.

    Diao Rui,Shan Lianyu,Shang Xinmin et al.The denoising method for microseismic monitoring data in time-frequency domain.Geophysical and Geochemical Exploration,2015,39(1):112-117.

    [9]刁瑞,吳國(guó)忱,尚新民等.地面陣列式微地震數(shù)據(jù)盲源分離去噪方法.物探與化探,2017,41(3):521-526.

    Diao Rui,Wu Guochen,Shang Xinmin et al.The blind separation denoising method for surface array micro-seismic data.Geophysical and Geochemical Exploration,2017,41(3):521-526.

    [10]賈瑞生,趙同彬,孫紅梅等.基于經(jīng)驗(yàn)?zāi)B(tài)分解及獨(dú)立成分分析的微震信號(hào)降噪方法.地球物理學(xué)報(bào),2015,58(3):1013-1023.

    Jia Ruisheng,Zhao Tongbin,Sun Hongmei et al.Microseismic signal denoising method based on empirical mode decomposition and independent component analysis.Chinese Journal of Geophysics,2015,58(3):1013-1023.

    [11]李會(huì)儉,王潤(rùn)秋,曹思遠(yuǎn)等.利用多尺度形態(tài)學(xué)識(shí)別微地震監(jiān)測(cè)中的弱信號(hào).石油地球物理勘探,2015,50(6):1105-1111.

    Li Huijian,Wang Runqiu,Cao Siyuan et al.Weak signal identification in microseismic monitoring with multi-scale morphology.OGP,2015,50(6):1105-1111.

    [12]宋維琪,劉太偉.地面微地震資料τ-p變換噪聲壓制.石油地球物理勘探,2015,50(1):48-53.

    Song Weiqi,Liu Taiwei.Surface microseismic noise suppression withτ-ptransform.OGP,2015,50(1):48-53.

    [13]崔慶輝,尹成,刁瑞等.地面微地震監(jiān)測(cè)數(shù)據(jù)處理難點(diǎn)及對(duì)策.油氣藏評(píng)價(jià)與開發(fā),2017,7(1):7-13.

    Cui Qinghui,Yin Cheng,Diao Rui et al.Difficulties and countermeasure of surface microseismic monitoring data processing.Reservoir Evaluation and Development,2017,7(1):7-13.

    [14]趙爭(zhēng)光,秦月霜,楊瑞召.地面微地震監(jiān)測(cè)致密砂巖儲(chǔ)層水力裂縫.地球物理學(xué)進(jìn)展,2014,29(5):2136-2139.

    Zhao Zhengguang,Qin Yushuang,Yang Zhaorui.Hydraulic fracture mapping for a tight sands reservoir by surface based microseismic monitoring.Progress in Geophysics,2014,29(5):2136-2139.

    [15]譚玉陽,于靜,馮剛等.微地震事件初至拾取SLPEA算法.地球物理學(xué)報(bào),2016,59(1):185-196.

    Tan Yuyang,Yu Jing,Feng Gang et al.Arrival picking of microseismic events using the SLPEA algorithm.Chinese Journal of Geophysics,2016,59(1):185-196.

    [16]曹俊海,顧漢明,尚新民.基于局部相關(guān)譜約束的多道匹配追蹤算法識(shí)別微地震信號(hào).石油地球物理勘探,2017,54(4):704-714.

    Cao Junhai,Gu Hanming,Shang Xinmin.Microseismic signal identification with multichannel matching pursuit based on local coherence spectrum constraint.OGP,2017,54(4):704-714.

    [17]劉勁松,王赟,姚振興.微地震信號(hào)到時(shí)自動(dòng)拾取方法.地球物理學(xué)報(bào),2013,56(5):1660-1666.

    Liu Jinsong,Wang Yun,Yao Zhenxing.On micro-seismic first arrival identification:A case study.Chinese Journal of Geophysics, 2013,56(5):1660-1666.

    [18]王鵬,常旭,王一博.基于時(shí)頻稀疏性分析法的低信噪比微震事件識(shí)別與恢復(fù).地球物理學(xué)報(bào),2014,57(8): 2656-2665.

    Wang Peng,Chang Xu,Wang Yibo.Automatic event detection and event recovery in low SNR microseismic signals based on time-frequency sparseness.Chinese Journal of Geophysics,2014,57(8):2656-2665.

    [19]刁瑞,胡曉婷,崔慶輝等.微地震監(jiān)測(cè)中等效速度定位精度分析及應(yīng)用.油氣藏評(píng)價(jià)與開發(fā),2015,5(1):49-53.

    Diao Rui,Hu Xiaoting,Cui Qinghui et al.Precision analysis and application of equivalent velocity in micro-seismic monitor.Reservoir Evaluation and Development,2015,5(1):49-53.

    [20]王維波,周瑤琪,春蘭.地面微地震監(jiān)測(cè)SET震源定位特性研究.中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,36(5):45-54.

    Wang Weibo,Zhou Yaoqi,Chun Lan.Characteristics of source localization by seismic emission tomography for surface based on microseismic monitoring.Journal

    of China University of Petroleum,2012,36(5):45-54.

    [21]譚玉陽,何川,張洪亮等.基于初至旅行時(shí)差的微地震速度模型反演.石油地球物理勘探,2015,50(1):54-59.

    Tan Yuyang,He Chuang,Zhang Hongliang et al.Microseismic velocity model inversion based on moveouts of first arrivals.OGP,2015,50(1):54-59.

    [22]高少武,趙波,祝樹云等.余弦函數(shù)自適應(yīng)法識(shí)別與消除單頻干擾.石油地球物理勘探,2011,46(1):64-69.

    Gao Shaowu,Zhao Bo,Zhu Shuyun et al.Identification and suppression of single frequency interference by using cosine function adaptive method.OGP,2011,46(1):64-69.

    [23]王納申,張譯丹,黃家旋等.用微地震技術(shù)評(píng)價(jià)姬源油田體積壓裂的效果.物探與化探,2017,41(1):165-170.

    Wang Nashen,Zhang Yidan,Huang Jiaxuan et al.The application of micro-seismic monitoring and evaluation technology to evolution of volume fracturing effect in the Jiyuan oilfield.Geophysical and Geochemical Exploration,2017,41(1): 165-170.

    [24]張?jiān)沏y,劉海寧,李紅梅等.應(yīng)用微地震監(jiān)測(cè)數(shù)據(jù)估算儲(chǔ)層壓裂改造體積.石油地球物理勘探,2017,52(2):309-314.

    Zhang Yunyin,Liu Haining,Li Hongmei et al.Reservoir fracturing volume estimation with microseismic monitoring data.OGP,2017,52(2):309-314.

    猜你喜歡
    時(shí)窗傳播速度噪聲源
    代謝綜合征患者臂踝脈搏波傳播速度與頸動(dòng)脈粥樣硬化的關(guān)系
    GRAPES-GFS模式2 m溫度預(yù)報(bào)的最優(yōu)時(shí)窗滑動(dòng)訂正方法
    氣象科技(2021年5期)2021-11-02 01:36:28
    一種基于改進(jìn)時(shí)窗法的爆炸沖擊波檢測(cè)方法
    汽車后視鏡-A柱區(qū)域氣動(dòng)噪聲源特征識(shí)別
    汽車工程(2021年12期)2021-03-08 02:34:18
    新雷
    一類廣義canmassa—Holm方程的無限傳播速度與漸近行為
    一種基于相位增量隨機(jī)化的寬帶噪聲源產(chǎn)生技術(shù)
    不同介入時(shí)窗和療程對(duì)高壓氧治療急性脊髓損傷的影響研究
    UWB搜救生物雷達(dá)分段時(shí)窗探測(cè)系統(tǒng)的實(shí)現(xiàn)
    傳媒全球化語境下的媒介話語批判
    五月玫瑰六月丁香| 视频中文字幕在线观看| 国产精华一区二区三区| 国产久久久一区二区三区| 国产亚洲精品av在线| 国产高潮美女av| 亚洲av成人精品一二三区| 日本免费a在线| 秋霞在线观看毛片| 小蜜桃在线观看免费完整版高清| 1000部很黄的大片| 国产成人免费观看mmmm| 你懂的网址亚洲精品在线观看 | 国产精品美女特级片免费视频播放器| 国产综合懂色| 高清午夜精品一区二区三区| 最近最新中文字幕免费大全7| 女人十人毛片免费观看3o分钟| 国产精品国产三级国产专区5o | 日本爱情动作片www.在线观看| 午夜激情欧美在线| 在线播放无遮挡| 午夜福利高清视频| 国产精品日韩av在线免费观看| 1024手机看黄色片| 国产精品麻豆人妻色哟哟久久 | av免费在线看不卡| 男插女下体视频免费在线播放| 久久精品国产亚洲av天美| 少妇熟女aⅴ在线视频| 亚洲欧美精品自产自拍| 久久这里只有精品中国| av卡一久久| 免费黄网站久久成人精品| 国产成人freesex在线| 亚洲欧美精品专区久久| 综合色丁香网| 成人亚洲精品av一区二区| 久久久久免费精品人妻一区二区| 国产亚洲5aaaaa淫片| 97在线视频观看| 亚洲国产精品成人久久小说| 欧美潮喷喷水| 日韩av在线免费看完整版不卡| 国产一区二区在线av高清观看| 久久久色成人| 亚洲av中文字字幕乱码综合| 汤姆久久久久久久影院中文字幕 | 亚洲欧美日韩卡通动漫| 免费看日本二区| 国产伦一二天堂av在线观看| 看非洲黑人一级黄片| 校园人妻丝袜中文字幕| 别揉我奶头 嗯啊视频| 国产黄片美女视频| 婷婷六月久久综合丁香| 亚洲高清免费不卡视频| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 亚洲欧洲国产日韩| 午夜老司机福利剧场| 亚洲成人久久爱视频| 啦啦啦观看免费观看视频高清| 日本黄大片高清| 国产一区亚洲一区在线观看| av女优亚洲男人天堂| 免费观看精品视频网站| 午夜福利高清视频| 97在线视频观看| av在线播放精品| 婷婷色av中文字幕| 一个人看的www免费观看视频| 草草在线视频免费看| 毛片一级片免费看久久久久| 亚洲国产精品国产精品| 大话2 男鬼变身卡| 国产视频首页在线观看| 又爽又黄无遮挡网站| 特级一级黄色大片| 男人舔女人下体高潮全视频| 97热精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 精品不卡国产一区二区三区| 99在线人妻在线中文字幕| 日韩 亚洲 欧美在线| 小蜜桃在线观看免费完整版高清| 91aial.com中文字幕在线观看| 天美传媒精品一区二区| 内地一区二区视频在线| 亚洲第一区二区三区不卡| 亚洲熟妇中文字幕五十中出| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品av在线| 国产一区二区三区av在线| 久久99热这里只频精品6学生 | 日产精品乱码卡一卡2卡三| 我要看日韩黄色一级片| 国产91av在线免费观看| 熟女电影av网| 日韩欧美三级三区| 日本免费a在线| 国产成人aa在线观看| 亚洲人成网站在线播| 联通29元200g的流量卡| 天天一区二区日本电影三级| 久久久久久久久大av| 亚洲av免费高清在线观看| 精品国内亚洲2022精品成人| 村上凉子中文字幕在线| 久久精品久久精品一区二区三区| 国产又黄又爽又无遮挡在线| 波多野结衣高清无吗| 日韩成人av中文字幕在线观看| 国产v大片淫在线免费观看| 国产精品国产高清国产av| 九九久久精品国产亚洲av麻豆| 国产在视频线精品| 日韩在线高清观看一区二区三区| 黄片wwwwww| 精品久久久久久久末码| 99久久精品一区二区三区| 亚洲精品日韩在线中文字幕| 欧美激情在线99| 午夜久久久久精精品| 国产极品精品免费视频能看的| 国产成人福利小说| 久热久热在线精品观看| 亚洲中文字幕一区二区三区有码在线看| 国产国拍精品亚洲av在线观看| 亚洲人与动物交配视频| 九草在线视频观看| 亚洲国产欧洲综合997久久,| 国产视频内射| 最近中文字幕2019免费版| 又粗又硬又长又爽又黄的视频| 国产美女午夜福利| 成年版毛片免费区| 有码 亚洲区| 日日干狠狠操夜夜爽| 中文字幕免费在线视频6| 小说图片视频综合网站| 久久久久久久午夜电影| 午夜精品国产一区二区电影 | 老女人水多毛片| 日韩欧美 国产精品| 国产高清国产精品国产三级 | av免费在线看不卡| 非洲黑人性xxxx精品又粗又长| 国产精品99久久久久久久久| 日韩欧美三级三区| 欧美丝袜亚洲另类| 色网站视频免费| 床上黄色一级片| 亚洲人与动物交配视频| 免费看光身美女| 搡女人真爽免费视频火全软件| 久久久色成人| 久久亚洲国产成人精品v| 黄片wwwwww| 午夜a级毛片| 国产黄片视频在线免费观看| 日韩欧美国产在线观看| 神马国产精品三级电影在线观看| 久久精品久久精品一区二区三区| 老女人水多毛片| 亚洲精品日韩在线中文字幕| 成人无遮挡网站| 啦啦啦韩国在线观看视频| 在线观看av片永久免费下载| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久久丰满| 亚洲欧美精品自产自拍| 69人妻影院| 老师上课跳d突然被开到最大视频| 纵有疾风起免费观看全集完整版 | 午夜免费激情av| 国产伦精品一区二区三区四那| 欧美日韩在线观看h| 日韩高清综合在线| 伦精品一区二区三区| 久久久久久久久大av| 51国产日韩欧美| 成年版毛片免费区| 日本猛色少妇xxxxx猛交久久| 嘟嘟电影网在线观看| 美女脱内裤让男人舔精品视频| 国产亚洲午夜精品一区二区久久 | 免费在线观看成人毛片| 精品不卡国产一区二区三区| 日韩成人av中文字幕在线观看| 老女人水多毛片| 只有这里有精品99| 久久午夜福利片| 亚洲av中文字字幕乱码综合| 亚洲综合色惰| 毛片女人毛片| 夜夜爽夜夜爽视频| 成人美女网站在线观看视频| 能在线免费观看的黄片| 精品99又大又爽又粗少妇毛片| 欧美一级a爱片免费观看看| 色综合色国产| 91久久精品电影网| 色综合亚洲欧美另类图片| 青青草视频在线视频观看| 久久久久网色| 欧美一级a爱片免费观看看| av福利片在线观看| 国产精品久久久久久精品电影| 超碰av人人做人人爽久久| 国产在视频线精品| 亚洲av不卡在线观看| 午夜视频国产福利| 三级国产精品欧美在线观看| 日韩大片免费观看网站 | 国产片特级美女逼逼视频| 麻豆成人av视频| 色吧在线观看| 中文字幕亚洲精品专区| 精品人妻偷拍中文字幕| 国语对白做爰xxxⅹ性视频网站| 日韩欧美在线乱码| 搡女人真爽免费视频火全软件| 熟女人妻精品中文字幕| 精品国产露脸久久av麻豆 | 午夜激情欧美在线| 亚洲欧美日韩无卡精品| 男女下面进入的视频免费午夜| 免费播放大片免费观看视频在线观看 | av又黄又爽大尺度在线免费看 | 在线免费观看不下载黄p国产| 日韩高清综合在线| 国产精品永久免费网站| 亚洲国产精品合色在线| 五月伊人婷婷丁香| 国产视频首页在线观看| 日韩成人av中文字幕在线观看| 日韩亚洲欧美综合| 午夜福利成人在线免费观看| 成人一区二区视频在线观看| a级一级毛片免费在线观看| 特级一级黄色大片| 国产精品国产三级国产av玫瑰| 国产黄色视频一区二区在线观看 | 国语对白做爰xxxⅹ性视频网站| 国产激情偷乱视频一区二区| 高清日韩中文字幕在线| 极品教师在线视频| 在线观看av片永久免费下载| 亚洲va在线va天堂va国产| 日韩在线高清观看一区二区三区| 国产亚洲91精品色在线| 国产在线一区二区三区精 | 亚洲欧美日韩高清专用| 亚洲在久久综合| 亚洲国产精品专区欧美| 婷婷色av中文字幕| 91午夜精品亚洲一区二区三区| 好男人视频免费观看在线| 人人妻人人看人人澡| 亚洲欧美成人综合另类久久久 | 黄片wwwwww| 亚洲av免费高清在线观看| 国产午夜精品论理片| 七月丁香在线播放| 欧美极品一区二区三区四区| 久久久午夜欧美精品| 午夜福利在线观看吧| 欧美日本视频| 亚洲av成人av| 国产av一区在线观看免费| 久久久午夜欧美精品| 男插女下体视频免费在线播放| 色5月婷婷丁香| 在线免费观看的www视频| 亚洲一区高清亚洲精品| 久久热精品热| 2021少妇久久久久久久久久久| 熟妇人妻久久中文字幕3abv| 99在线人妻在线中文字幕| 成人亚洲精品av一区二区| 中文天堂在线官网| 亚洲精品国产av成人精品| 欧美人与善性xxx| 九色成人免费人妻av| 99热这里只有精品一区| 亚洲自拍偷在线| 国产片特级美女逼逼视频| 一个人看的www免费观看视频| 久久久久久大精品| 国产精品麻豆人妻色哟哟久久 | 久久6这里有精品| 看片在线看免费视频| 亚洲伊人久久精品综合 | 国内精品美女久久久久久| av在线播放精品| 可以在线观看毛片的网站| 日本wwww免费看| 建设人人有责人人尽责人人享有的 | 村上凉子中文字幕在线| 国产成人91sexporn| 欧美性猛交╳xxx乱大交人| 日本爱情动作片www.在线观看| 日韩成人av中文字幕在线观看| 亚洲欧洲日产国产| 日本免费a在线| 中文在线观看免费www的网站| 91aial.com中文字幕在线观看| av女优亚洲男人天堂| 人体艺术视频欧美日本| 岛国在线免费视频观看| 国产伦在线观看视频一区| АⅤ资源中文在线天堂| 亚洲在线自拍视频| av免费在线看不卡| 久久精品国产亚洲av天美| 国产三级在线视频| 三级毛片av免费| 小蜜桃在线观看免费完整版高清| 日韩在线高清观看一区二区三区| 一夜夜www| 国产69精品久久久久777片| 国产成人一区二区在线| 精品国产三级普通话版| 久久精品影院6| 熟女人妻精品中文字幕| 久久久久久久久中文| 日本免费一区二区三区高清不卡| 九九久久精品国产亚洲av麻豆| 国产一区二区在线av高清观看| 毛片女人毛片| 久久精品国产亚洲av天美| 女人久久www免费人成看片 | 欧美成人一区二区免费高清观看| 蜜桃久久精品国产亚洲av| 色网站视频免费| 久久99热这里只有精品18| 日本免费a在线| 人妻系列 视频| 99热精品在线国产| 日本免费在线观看一区| av免费观看日本| 久久人人爽人人爽人人片va| av视频在线观看入口| 亚洲国产精品成人综合色| eeuss影院久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩av在线免费看完整版不卡| 亚洲综合色惰| 99热这里只有精品一区| 久久亚洲精品不卡| 亚洲人成网站在线播| 91aial.com中文字幕在线观看| 秋霞在线观看毛片| 久久久午夜欧美精品| 日韩 亚洲 欧美在线| 大话2 男鬼变身卡| 亚洲av免费高清在线观看| or卡值多少钱| 中文字幕av在线有码专区| 又粗又硬又长又爽又黄的视频| 真实男女啪啪啪动态图| 免费不卡的大黄色大毛片视频在线观看 | 99久久人妻综合| 精品久久久久久久久亚洲| 国产亚洲一区二区精品| av播播在线观看一区| 少妇人妻精品综合一区二区| 男女啪啪激烈高潮av片| 99久国产av精品国产电影| 国产精品久久久久久久电影| 能在线免费看毛片的网站| h日本视频在线播放| 两个人视频免费观看高清| 伊人久久精品亚洲午夜| 一区二区三区四区激情视频| 亚洲国产精品合色在线| 91精品伊人久久大香线蕉| 九草在线视频观看| 视频中文字幕在线观看| 亚洲最大成人中文| 久久6这里有精品| 国产伦在线观看视频一区| 听说在线观看完整版免费高清| kizo精华| 观看美女的网站| 特级一级黄色大片| 亚洲aⅴ乱码一区二区在线播放| 国产探花极品一区二区| 两个人视频免费观看高清| 久久鲁丝午夜福利片| 少妇裸体淫交视频免费看高清| 色5月婷婷丁香| 在线观看美女被高潮喷水网站| 亚洲性久久影院| av在线天堂中文字幕| 日韩成人伦理影院| 国产黄色小视频在线观看| 欧美人与善性xxx| 精品免费久久久久久久清纯| 黄色一级大片看看| 好男人视频免费观看在线| 一级黄色大片毛片| 亚洲精品国产av成人精品| 亚洲国产高清在线一区二区三| 男人的好看免费观看在线视频| 夜夜看夜夜爽夜夜摸| 1024手机看黄色片| 汤姆久久久久久久影院中文字幕 | 国产亚洲av片在线观看秒播厂 | 成人漫画全彩无遮挡| 国产成人一区二区在线| 色播亚洲综合网| 国产黄片美女视频| 国产在视频线精品| 国产免费一级a男人的天堂| 国产片特级美女逼逼视频| 男女啪啪激烈高潮av片| 男人狂女人下面高潮的视频| 午夜a级毛片| 色尼玛亚洲综合影院| 国产成人免费观看mmmm| 国产探花在线观看一区二区| 九九热线精品视视频播放| 亚洲精华国产精华液的使用体验| 狠狠狠狠99中文字幕| 欧美不卡视频在线免费观看| 亚洲国产最新在线播放| 人妻夜夜爽99麻豆av| 熟女人妻精品中文字幕| 日韩欧美精品v在线| 女的被弄到高潮叫床怎么办| 国产私拍福利视频在线观看| 精品国内亚洲2022精品成人| 欧美性感艳星| 久久久久久久久久久免费av| 男人的好看免费观看在线视频| 精品久久久久久久久亚洲| 免费看av在线观看网站| 国内揄拍国产精品人妻在线| 国产 一区 欧美 日韩| 99在线人妻在线中文字幕| 国产又色又爽无遮挡免| 搡老妇女老女人老熟妇| 免费av不卡在线播放| 国产精品熟女久久久久浪| 在线观看av片永久免费下载| 91精品国产九色| 色综合亚洲欧美另类图片| av在线老鸭窝| 国产高清国产精品国产三级 | 国产精品熟女久久久久浪| 男女啪啪激烈高潮av片| 国产精品蜜桃在线观看| 久久久久久久久大av| 一级爰片在线观看| 综合色av麻豆| 人妻制服诱惑在线中文字幕| 春色校园在线视频观看| 天天躁夜夜躁狠狠久久av| 超碰97精品在线观看| 男女下面进入的视频免费午夜| 丝袜喷水一区| 国产亚洲最大av| 国产精品人妻久久久久久| 久久99热这里只频精品6学生 | 久久6这里有精品| 亚洲中文字幕日韩| 少妇人妻精品综合一区二区| 26uuu在线亚洲综合色| 国产色爽女视频免费观看| 丰满乱子伦码专区| 成人高潮视频无遮挡免费网站| 国语对白做爰xxxⅹ性视频网站| ponron亚洲| 欧美3d第一页| 免费电影在线观看免费观看| 精品久久久久久久人妻蜜臀av| 日本av手机在线免费观看| 老司机影院毛片| 国产精品一区www在线观看| 2021天堂中文幕一二区在线观| 内地一区二区视频在线| 一边亲一边摸免费视频| 伦精品一区二区三区| av卡一久久| 身体一侧抽搐| 久久精品久久久久久久性| 国内揄拍国产精品人妻在线| 久久精品熟女亚洲av麻豆精品 | 亚洲最大成人手机在线| 欧美又色又爽又黄视频| 欧美成人午夜免费资源| 国产一区二区在线av高清观看| 国产探花极品一区二区| 久久国内精品自在自线图片| 神马国产精品三级电影在线观看| 精品久久久久久成人av| 亚洲av成人av| 日本熟妇午夜| 99久国产av精品国产电影| 自拍偷自拍亚洲精品老妇| 高清毛片免费看| 国产高清视频在线观看网站| 国产免费男女视频| 99久国产av精品国产电影| 人人妻人人澡欧美一区二区| 乱码一卡2卡4卡精品| 久久精品综合一区二区三区| 色综合色国产| 亚洲欧美日韩高清专用| 美女大奶头视频| videos熟女内射| 日本猛色少妇xxxxx猛交久久| 国产精品99久久久久久久久| 日韩欧美精品免费久久| 婷婷色麻豆天堂久久 | 久久精品综合一区二区三区| 网址你懂的国产日韩在线| 精品欧美国产一区二区三| 国产真实伦视频高清在线观看| 天堂影院成人在线观看| 99九九线精品视频在线观看视频| 最近中文字幕2019免费版| 久久人人爽人人片av| 久久久成人免费电影| 久久这里有精品视频免费| 69av精品久久久久久| 国内精品美女久久久久久| 蜜桃亚洲精品一区二区三区| 黑人高潮一二区| 尾随美女入室| 国产亚洲av嫩草精品影院| 亚洲真实伦在线观看| 久久精品国产鲁丝片午夜精品| 国产精品三级大全| 亚洲色图av天堂| 亚洲成人久久爱视频| 热99在线观看视频| 国产成人午夜福利电影在线观看| 天天躁日日操中文字幕| 久久欧美精品欧美久久欧美| 亚洲国产精品国产精品| 嫩草影院精品99| 三级男女做爰猛烈吃奶摸视频| 少妇熟女欧美另类| 亚洲av.av天堂| 哪个播放器可以免费观看大片| 老师上课跳d突然被开到最大视频| 九九热线精品视视频播放| 在线观看66精品国产| 午夜激情欧美在线| 国产精品一及| 国产免费男女视频| 国产淫片久久久久久久久| 国产一区二区在线观看日韩| 一级毛片aaaaaa免费看小| 99久久无色码亚洲精品果冻| 日韩av在线免费看完整版不卡| 亚洲国产精品国产精品| 精品久久国产蜜桃| 老司机影院毛片| 建设人人有责人人尽责人人享有的 | 欧美成人一区二区免费高清观看| 国产av在哪里看| 老师上课跳d突然被开到最大视频| 国产精品女同一区二区软件| 身体一侧抽搐| 九九爱精品视频在线观看| 欧美性猛交黑人性爽| 九九在线视频观看精品| 成年免费大片在线观看| 亚洲国产欧洲综合997久久,| 少妇熟女aⅴ在线视频| 在线免费十八禁| 深夜a级毛片| 中文在线观看免费www的网站| 国产黄a三级三级三级人| 赤兔流量卡办理| a级毛片免费高清观看在线播放| 国产精华一区二区三区| 男女那种视频在线观看| 欧美一区二区亚洲| 色网站视频免费| 一级黄片播放器| eeuss影院久久| 又黄又爽又刺激的免费视频.| 成人特级av手机在线观看| 在线播放无遮挡| 国产精品久久视频播放| ponron亚洲| 午夜福利成人在线免费观看| 国产精品一区二区三区四区久久| 日韩在线高清观看一区二区三区| 久久6这里有精品| 免费观看在线日韩| 麻豆精品久久久久久蜜桃| 免费看a级黄色片| 午夜福利成人在线免费观看| 久久久久久久久大av| 美女被艹到高潮喷水动态| 偷拍熟女少妇极品色| 色综合色国产| 99久国产av精品国产电影| 婷婷色麻豆天堂久久 | 免费观看在线日韩| 最近的中文字幕免费完整| 亚洲自拍偷在线| 国产乱人偷精品视频| 有码 亚洲区| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 国产av在哪里看| 国产在线一区二区三区精 | 最近的中文字幕免费完整|