• <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免费久了| 看片在线看免费视频| 午夜精品久久久久久毛片777| 男人舔女人的私密视频| 亚洲五月色婷婷综合| 国产亚洲精品综合一区在线观看 | 黄片小视频在线播放| 变态另类成人亚洲欧美熟女 | 午夜视频精品福利| 每晚都被弄得嗷嗷叫到高潮| 色综合婷婷激情| 亚洲在线自拍视频| 中文字幕色久视频| 亚洲人成电影观看| 狠狠狠狠99中文字幕| 国产91精品成人一区二区三区| 国产精品综合久久久久久久免费 | 久久久国产成人精品二区| 老鸭窝网址在线观看| 亚洲成av片中文字幕在线观看| 动漫黄色视频在线观看| 高潮久久久久久久久久久不卡| 91国产中文字幕| 国产精品久久久久久精品电影 | 欧美国产日韩亚洲一区| 亚洲一码二码三码区别大吗| 女人被躁到高潮嗷嗷叫费观| 国产精品98久久久久久宅男小说| 女人精品久久久久毛片| 国产精品一区二区三区四区久久 | 日韩一卡2卡3卡4卡2021年| 十八禁网站免费在线| 在线观看www视频免费| 老司机午夜十八禁免费视频| 国产黄a三级三级三级人| 久久精品aⅴ一区二区三区四区| 成人三级黄色视频| 成人av一区二区三区在线看| 亚洲 欧美 日韩 在线 免费| 窝窝影院91人妻| 亚洲精品国产色婷婷电影| 久9热在线精品视频| 欧美在线一区亚洲| 91精品三级在线观看| 身体一侧抽搐| 两个人看的免费小视频| 波多野结衣一区麻豆| 在线观看午夜福利视频| 亚洲午夜精品一区,二区,三区| 一级片免费观看大全| 国产一区二区在线av高清观看| 亚洲国产毛片av蜜桃av| 成年版毛片免费区| 亚洲 欧美一区二区三区| 国产在线精品亚洲第一网站| 校园春色视频在线观看| 女人被狂操c到高潮| 真人一进一出gif抽搐免费| 免费搜索国产男女视频| 精品国产乱码久久久久久男人| 90打野战视频偷拍视频| 国产成人精品在线电影| 国产不卡一卡二| 丰满的人妻完整版| 韩国精品一区二区三区| 国产成+人综合+亚洲专区| 中文字幕色久视频| 亚洲国产精品成人综合色| 欧美黄色片欧美黄色片| 夜夜看夜夜爽夜夜摸| 久久久久久久久中文| 亚洲人成电影免费在线| 国产精华一区二区三区| 欧美成人性av电影在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲精品一区av在线观看| 乱人伦中国视频| 777久久人妻少妇嫩草av网站| 91在线观看av| 亚洲午夜理论影院| 精品国内亚洲2022精品成人| 成人18禁高潮啪啪吃奶动态图| 国产麻豆成人av免费视频| 国产精品精品国产色婷婷| 国产主播在线观看一区二区| 中文字幕av电影在线播放| 精品国产一区二区久久| 日韩高清综合在线| 麻豆久久精品国产亚洲av| 久久热在线av| 国产精品精品国产色婷婷| 大香蕉久久成人网| 在线av久久热| 国产亚洲精品综合一区在线观看 | 亚洲精品国产色婷婷电影| 男人的好看免费观看在线视频 | 九色国产91popny在线| 亚洲avbb在线观看| bbb黄色大片| 久久久久久久午夜电影| videosex国产| 天堂√8在线中文| 国产av精品麻豆| 91国产中文字幕| 免费不卡黄色视频| 精品久久蜜臀av无| 波多野结衣巨乳人妻| 日本黄色视频三级网站网址| 午夜福利一区二区在线看| 日本免费a在线| 日韩有码中文字幕| 国产人伦9x9x在线观看| 亚洲免费av在线视频| 午夜精品在线福利| 一区二区三区精品91| 中国美女看黄片| av福利片在线| 免费不卡黄色视频| 亚洲av电影不卡..在线观看| www.熟女人妻精品国产| 日本在线视频免费播放| 欧美日韩亚洲国产一区二区在线观看| 亚洲avbb在线观看| 十分钟在线观看高清视频www| 精品欧美一区二区三区在线| 99久久精品国产亚洲精品| 精品一区二区三区四区五区乱码| 少妇的丰满在线观看| 老司机福利观看| 亚洲欧美精品综合久久99| 91麻豆av在线| 色尼玛亚洲综合影院| 高潮久久久久久久久久久不卡| 亚洲自偷自拍图片 自拍| 美女国产高潮福利片在线看| 日本黄色视频三级网站网址| 丁香欧美五月| 黄色女人牲交| 老汉色∧v一级毛片| 日韩成人在线观看一区二区三区| 欧美黑人精品巨大| 成人三级黄色视频| 国产高清激情床上av| 亚洲精品美女久久av网站| 天堂√8在线中文| 黄片小视频在线播放| 99国产精品一区二区蜜桃av| 搞女人的毛片| 黄色a级毛片大全视频| 男男h啪啪无遮挡| 亚洲五月色婷婷综合| 久久久国产欧美日韩av| 久久久久久大精品| 亚洲人成电影免费在线| 淫妇啪啪啪对白视频| 国产乱人伦免费视频| 啦啦啦免费观看视频1| 极品教师在线免费播放| 国产精品美女特级片免费视频播放器 | 国产日韩一区二区三区精品不卡| 9色porny在线观看| 在线国产一区二区在线| 国产亚洲精品av在线| 老司机在亚洲福利影院| 亚洲 欧美 日韩 在线 免费| 国产99白浆流出| 国产精品影院久久| 男女下面进入的视频免费午夜 | 麻豆一二三区av精品| 国产精品免费视频内射| 国产又爽黄色视频| 国产伦人伦偷精品视频| 欧美色欧美亚洲另类二区 | 午夜福利高清视频| 男人舔女人下体高潮全视频| 亚洲av成人av| 精品国产乱码久久久久久男人| 久久久久久久久久久久大奶| 视频区欧美日本亚洲| 天堂√8在线中文| 91av网站免费观看| 久久久久久久久久久久大奶| 国产色视频综合| 国内精品久久久久久久电影| 非洲黑人性xxxx精品又粗又长| 免费在线观看黄色视频的| 自拍欧美九色日韩亚洲蝌蚪91| 熟女少妇亚洲综合色aaa.| 亚洲精品国产色婷婷电影| 精品高清国产在线一区| 最近最新中文字幕大全电影3 | 国产欧美日韩精品亚洲av| 国产一区二区激情短视频| 人人妻人人爽人人添夜夜欢视频| 亚洲久久久国产精品| 十八禁网站免费在线| 90打野战视频偷拍视频| 免费女性裸体啪啪无遮挡网站| 少妇裸体淫交视频免费看高清 | 大香蕉久久成人网| 亚洲五月天丁香| а√天堂www在线а√下载| 嫁个100分男人电影在线观看| 亚洲avbb在线观看| 精品久久久久久久人妻蜜臀av | 久久精品亚洲精品国产色婷小说| 亚洲色图av天堂| 欧美日韩精品网址| 悠悠久久av| 免费少妇av软件| 国产欧美日韩一区二区三| 久久久久久国产a免费观看| 国产一区二区三区综合在线观看| 国产黄a三级三级三级人| 亚洲男人天堂网一区| 欧美激情极品国产一区二区三区| 精品乱码久久久久久99久播| 欧美老熟妇乱子伦牲交| 亚洲成人国产一区在线观看| 亚洲久久久国产精品| 深夜精品福利| 美女高潮喷水抽搐中文字幕| 亚洲精品美女久久av网站| 51午夜福利影视在线观看| 1024视频免费在线观看| 老司机午夜十八禁免费视频| 九色国产91popny在线| 啦啦啦韩国在线观看视频| 女同久久另类99精品国产91| 日韩欧美一区二区三区在线观看| 老鸭窝网址在线观看| 老司机午夜十八禁免费视频| 久久午夜综合久久蜜桃| 国内精品久久久久精免费| 一边摸一边抽搐一进一出视频| 搡老岳熟女国产| 久久精品人人爽人人爽视色| 高清毛片免费观看视频网站| 色av中文字幕| 国产精品99久久99久久久不卡| 亚洲伊人色综图| 久久久久久国产a免费观看| 在线观看免费日韩欧美大片| 欧美中文日本在线观看视频| 怎么达到女性高潮| 亚洲av五月六月丁香网| 久久中文字幕一级| 久久久久久久精品吃奶| 久久这里只有精品19| 亚洲精品国产一区二区精华液| 9热在线视频观看99| 黄片小视频在线播放| 日韩大码丰满熟妇| www国产在线视频色| 99久久99久久久精品蜜桃| 国产一区二区三区视频了| 色综合亚洲欧美另类图片| 一级,二级,三级黄色视频| 亚洲男人的天堂狠狠| 一本大道久久a久久精品| 色哟哟哟哟哟哟| 日韩欧美一区二区三区在线观看| 国产精品av久久久久免费| 日本一区二区免费在线视频| 久久青草综合色| 亚洲电影在线观看av| 欧美在线黄色| 国产亚洲欧美在线一区二区| 日韩高清综合在线| 成人18禁在线播放| 又黄又粗又硬又大视频| 身体一侧抽搐| 啦啦啦 在线观看视频| 最近最新免费中文字幕在线| 首页视频小说图片口味搜索| 免费一级毛片在线播放高清视频 | 国产不卡一卡二| 亚洲国产精品成人综合色| 少妇的丰满在线观看| 男男h啪啪无遮挡| 精品电影一区二区在线| 男女下面进入的视频免费午夜 | 日本免费一区二区三区高清不卡 | 免费久久久久久久精品成人欧美视频| 波多野结衣巨乳人妻| 少妇粗大呻吟视频| 午夜福利,免费看| 亚洲国产精品成人综合色| 亚洲国产精品成人综合色| 在线av久久热| 国产成人精品久久二区二区免费| 久久久久久免费高清国产稀缺| 国产午夜福利久久久久久| 日韩 欧美 亚洲 中文字幕| 精品久久久久久,| 欧美黑人欧美精品刺激| a在线观看视频网站| 午夜久久久在线观看| 国产精品1区2区在线观看.| av中文乱码字幕在线| 国产精品日韩av在线免费观看 | 亚洲人成网站在线播放欧美日韩| 国产精品一区二区三区四区久久 | 一进一出好大好爽视频| 久久精品成人免费网站| 老司机福利观看| 久久影院123| 侵犯人妻中文字幕一二三四区| 色综合亚洲欧美另类图片| 国产色视频综合| 欧美一级a爱片免费观看看 | 久久香蕉国产精品| 高清毛片免费观看视频网站| 啦啦啦韩国在线观看视频| 黑人巨大精品欧美一区二区蜜桃| av在线天堂中文字幕| 精品熟女少妇八av免费久了| 欧美大码av| 91九色精品人成在线观看| 久久欧美精品欧美久久欧美| 在线观看一区二区三区| 精品国产乱码久久久久久男人| 精品国产一区二区久久| 咕卡用的链子| 婷婷精品国产亚洲av在线| 欧美黑人欧美精品刺激| 日韩高清综合在线| 午夜久久久久精精品| 欧美性长视频在线观看| 波多野结衣高清无吗| 变态另类丝袜制服| 欧美成狂野欧美在线观看| 国产日韩一区二区三区精品不卡| 国产亚洲av高清不卡| 涩涩av久久男人的天堂| 国产av精品麻豆| 精品欧美国产一区二区三| 久久久久精品国产欧美久久久| 在线观看免费视频日本深夜| av欧美777| 99riav亚洲国产免费| 午夜久久久久精精品| 亚洲专区国产一区二区| 中文字幕人成人乱码亚洲影| 满18在线观看网站| 91麻豆av在线| 色婷婷久久久亚洲欧美| 成人三级黄色视频| 久久人妻福利社区极品人妻图片| 久久九九热精品免费| 久久热在线av| 丝袜在线中文字幕| 88av欧美| 国产熟女午夜一区二区三区| 人妻久久中文字幕网| 老司机深夜福利视频在线观看| 999久久久国产精品视频| 十八禁网站免费在线| 日韩高清综合在线| 日本在线视频免费播放| 亚洲中文字幕一区二区三区有码在线看 | 国产高清videossex| 纯流量卡能插随身wifi吗| 日韩欧美三级三区| 亚洲人成电影免费在线| 一级毛片女人18水好多| 天天躁夜夜躁狠狠躁躁| 女生性感内裤真人,穿戴方法视频| 亚洲色图 男人天堂 中文字幕| 黑人欧美特级aaaaaa片| 久久国产精品影院| 亚洲成人国产一区在线观看| 青草久久国产| 亚洲中文日韩欧美视频| 亚洲自拍偷在线| 麻豆av在线久日| 9191精品国产免费久久| 免费高清视频大片| 色综合欧美亚洲国产小说| 在线观看www视频免费| 免费观看人在逋| 久久国产精品人妻蜜桃| 国产精品综合久久久久久久免费 | 两性夫妻黄色片| 精品久久久精品久久久| 亚洲视频免费观看视频| 日本撒尿小便嘘嘘汇集6| 国产亚洲av高清不卡| 制服丝袜大香蕉在线| 九色国产91popny在线| 自线自在国产av| 这个男人来自地球电影免费观看| 久久久久久久久中文| 精品国产美女av久久久久小说| 国内精品久久久久久久电影| 少妇 在线观看| 欧美老熟妇乱子伦牲交| 久久亚洲真实| 大码成人一级视频| 国内精品久久久久精免费| 午夜两性在线视频| 国产真人三级小视频在线观看| 亚洲七黄色美女视频| 亚洲成a人片在线一区二区| 99国产精品99久久久久| 亚洲第一av免费看| 欧美色视频一区免费| 韩国av一区二区三区四区| 在线观看一区二区三区| 日本一区二区免费在线视频| 天天一区二区日本电影三级 | 国产伦人伦偷精品视频| 欧美国产日韩亚洲一区| 亚洲精品中文字幕在线视频| 午夜福利免费观看在线| 人人澡人人妻人| 天天躁夜夜躁狠狠躁躁| 国产高清有码在线观看视频 | 精品国产超薄肉色丝袜足j| 国产成人精品在线电影| 法律面前人人平等表现在哪些方面| 99久久久亚洲精品蜜臀av| 免费女性裸体啪啪无遮挡网站| 中文字幕色久视频| 亚洲 欧美 日韩 在线 免费| 制服丝袜大香蕉在线| 久久午夜综合久久蜜桃| 欧美日韩黄片免| 亚洲人成77777在线视频| 国产精品久久电影中文字幕| 免费在线观看影片大全网站| 国产三级黄色录像| 欧美丝袜亚洲另类 | 1024视频免费在线观看| 色尼玛亚洲综合影院| 两人在一起打扑克的视频| 亚洲熟妇熟女久久| 中文字幕色久视频| 亚洲一区二区三区色噜噜| 满18在线观看网站| 午夜免费鲁丝| 国产精品久久久av美女十八| 国产av又大| 亚洲精品一卡2卡三卡4卡5卡| 成在线人永久免费视频| 亚洲情色 制服丝袜| 亚洲国产看品久久| 最新在线观看一区二区三区| 亚洲五月色婷婷综合| 亚洲精品粉嫩美女一区| 99久久综合精品五月天人人| 精品少妇一区二区三区视频日本电影| 一二三四社区在线视频社区8| 国产成人av教育| 欧美性长视频在线观看| 日韩欧美国产一区二区入口| 久久精品国产亚洲av高清一级| 久久影院123| 91成人精品电影| 69精品国产乱码久久久| 亚洲欧美日韩另类电影网站| 国产精品av久久久久免费| 成人18禁高潮啪啪吃奶动态图| avwww免费| 十分钟在线观看高清视频www| 午夜精品国产一区二区电影| 成人手机av| 亚洲电影在线观看av| 国产精品久久视频播放| 国产精品亚洲美女久久久| 精品乱码久久久久久99久播| 久久精品国产99精品国产亚洲性色 | 成年女人毛片免费观看观看9| av免费在线观看网站| 亚洲av电影在线进入| 久久人人爽av亚洲精品天堂| 亚洲成人精品中文字幕电影| 日本一区二区免费在线视频| 一二三四在线观看免费中文在| 久久久久精品国产欧美久久久| 神马国产精品三级电影在线观看 | 欧美日本视频| 人人妻人人爽人人添夜夜欢视频| 色尼玛亚洲综合影院| 国产xxxxx性猛交| 人人妻人人澡人人看| 少妇粗大呻吟视频| av超薄肉色丝袜交足视频| 国产免费男女视频| 长腿黑丝高跟| 精品午夜福利视频在线观看一区| 身体一侧抽搐| 99国产精品一区二区蜜桃av| 日日摸夜夜添夜夜添小说| 国产主播在线观看一区二区| 人人妻人人澡人人看| 亚洲色图av天堂| 十八禁网站免费在线| 亚洲avbb在线观看| 午夜福利高清视频| 最近最新免费中文字幕在线| 麻豆一二三区av精品| 丰满的人妻完整版| 两人在一起打扑克的视频| 亚洲精品在线美女| 欧美激情高清一区二区三区| av天堂久久9| 亚洲少妇的诱惑av| 一级,二级,三级黄色视频| 精品久久久久久久人妻蜜臀av | 少妇被粗大的猛进出69影院| 久久久久国产一级毛片高清牌| 国产精品影院久久| 欧美成狂野欧美在线观看| www日本在线高清视频| 成人亚洲精品一区在线观看| 亚洲 欧美一区二区三区| 欧美黑人欧美精品刺激| 91老司机精品| 他把我摸到了高潮在线观看| 国产人伦9x9x在线观看| 黄色女人牲交| 国产精品 欧美亚洲| 日本 av在线| 国产成人啪精品午夜网站| 一级a爱片免费观看的视频| 女人被狂操c到高潮| 人人妻人人澡欧美一区二区 | 国产精品1区2区在线观看.| 黄频高清免费视频| 波多野结衣巨乳人妻| 亚洲第一青青草原| 欧美国产日韩亚洲一区| 国产成人系列免费观看| 午夜精品在线福利| 久久精品亚洲精品国产色婷小说| 天堂√8在线中文| 9色porny在线观看| av在线播放免费不卡| 欧美日韩福利视频一区二区| 久久 成人 亚洲| 国产在线观看jvid| 国产成人av教育| 色老头精品视频在线观看| 精品乱码久久久久久99久播| 美女扒开内裤让男人捅视频| 最好的美女福利视频网| 青草久久国产| 女同久久另类99精品国产91| 欧美日韩乱码在线| 热re99久久国产66热| 国产成年人精品一区二区| 法律面前人人平等表现在哪些方面| 女人高潮潮喷娇喘18禁视频| 欧美人与性动交α欧美精品济南到| 国产精品1区2区在线观看.| 大型av网站在线播放| netflix在线观看网站| 国产精品九九99| 久久人人爽av亚洲精品天堂| 精品国产美女av久久久久小说| 自拍欧美九色日韩亚洲蝌蚪91| 男女午夜视频在线观看| 热99re8久久精品国产| 12—13女人毛片做爰片一| 久久性视频一级片| 国语自产精品视频在线第100页| 一夜夜www| 色播在线永久视频| 欧美日韩中文字幕国产精品一区二区三区 | 十八禁人妻一区二区| 9热在线视频观看99| 亚洲精品中文字幕一二三四区| 国产亚洲精品av在线| 女性生殖器流出的白浆| 精品国产乱子伦一区二区三区| 在线永久观看黄色视频| 国产精品 欧美亚洲| 午夜a级毛片| 欧美成人免费av一区二区三区| 久久久国产欧美日韩av| 国产麻豆69| 97碰自拍视频| 久久这里只有精品19| 成人av一区二区三区在线看| 桃红色精品国产亚洲av| 亚洲人成网站在线播放欧美日韩| 久久精品91无色码中文字幕| 手机成人av网站| 亚洲熟妇熟女久久| 国产99久久九九免费精品| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 岛国视频午夜一区免费看| 最近最新中文字幕大全电影3 | 亚洲自偷自拍图片 自拍| 国产午夜精品久久久久久| 夜夜躁狠狠躁天天躁| 欧美黄色片欧美黄色片| 亚洲情色 制服丝袜| 人人妻人人爽人人添夜夜欢视频| 中文字幕久久专区| 国产一区二区三区综合在线观看| 91老司机精品| 国产黄a三级三级三级人| 中文字幕av电影在线播放| 黄色视频不卡| 精品高清国产在线一区| 天天添夜夜摸| 99久久99久久久精品蜜桃| 露出奶头的视频| 真人做人爱边吃奶动态| 露出奶头的视频| 香蕉久久夜色|