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

    適用于陸上石油勘探的地震檢波器

    2017-12-18 10:22:53
    石油地球物理勘探 2017年6期
    關鍵詞:動圈檢波器耦合

    魏 繼 東

    (中國石化石油工程地球物理公司勝利分公司,山東東營 257100)

    ·采集技術·

    適用于陸上石油勘探的地震檢波器

    魏 繼 東*

    (中國石化石油工程地球物理公司勝利分公司,山東東營 257100)

    陸上石油勘探中采用過多種基于不同工作原理的檢波器(如動圈式模擬檢波器、MEMS數字檢波器、渦流檢波器、電化學檢波器及光纖檢波器等),而實際地震數據采集中多年來一直是動圈式模擬檢波器(如SM-4、20DX等)占據主導地位。通過分析地震采集階段信號、噪聲的特性以及地震檢波器的性能指標及其與地球物理特征之間的關系,提出了檢波器類型對比試驗的設計原則、采集階段初步評估地震資料的方法。認為在當前去噪能力下(即通過后續(xù)壓噪手段、仍難以將機械屬性的環(huán)境噪聲、原生噪聲及次生噪聲衰減到與電噪聲屬性的本底噪聲相當的水平),對于垂直分量而言,以20DX為代表的動圈式檢波器仍是數據合格、費用合理、施工方便的高效率檢波器; 同時指出了某些新型檢波器存在的不足及今后陸用檢波器的發(fā)展方向。

    檢波器 地震勘探 動圈式檢波器 陸上

    1 引言

    從地質期望角度而言,我們希望識別數千米以下數米尺度的、具有工業(yè)價值的構造或者巖性油氣藏單元。與達成這一目標所需要的反射信息相比,陸上反射波法地震勘探由人工震源激發(fā)的反射波具有主頻低(吸收衰減)、信噪比低(環(huán)境噪聲、原生噪聲、次生噪聲、電噪聲影響)的特點。同時,受技術手段的限制,對由震源激發(fā)的地面地震波機械振動的檢測只能通過具有機電轉換能力的檢波器進行。所以,以上因素決定了陸上石油地震檢波存在三個基本現實,即主頻低、信噪比低、由檢波器間接測量。前二者是地震勘探面臨的“客觀存在”,第三點則反映了人類當前的“主觀能力”。

    自從20世紀30年代推出現今普遍采用的動圈式模擬檢波器以來[1],基于不同的工作原理,涌現了大量不同類型的檢波器,比如MEMS數字檢波器、渦流加速度檢波器、電化學檢波器、光纖檢波器等等。盡管每次隨著不同類型檢波器的出現,均會相應地進行一些試驗、試生產,但是當前野外施工中占據主導地位的,仍然是動圈式模擬檢波器,比如SM-4、20DX等。為什么某些看似具有“指標優(yōu)勢”的檢波器卻始終沒有取代傳統(tǒng)動圈式檢波器的主導地位呢?作者認為,對于這一問題的回答,需要厘清三個方面的問題: ①地震采集階段信號與噪聲的比較特性; ②地震檢波器的指標參數對地震數據的影響; ③檢波器選型試驗的設計與資料比對方法。

    2 采集階段信號與噪聲的特性比較

    在反射波法石油勘探中,信號是指有利于目標地質體識別的有效反射信號,噪聲則是對有效反射信號識別具有干擾作用的非有效信號。噪聲與信號的定義是相對的,在某些情況下,二者可互換[2]。

    對于陸上地震勘探來說,在地震采集(地震信號由震源激發(fā)到被磁帶記錄)階段,影響地震信號振幅、頻率等特征的因素主要包括爆炸子波、大地吸收衰減、組合效應、組內時差、檢波器—大地耦合、檢波器及地震儀的性能、參數等。

    根據噪聲的出現規(guī)律,可將噪聲分為規(guī)則噪聲和隨機噪聲兩類。前者包括面波、折射波、聲波等源自炮點的原生噪聲以及非源自炮點的次生噪聲、具有相干性的環(huán)境噪聲;后者包括風吹草動等引起的頻帶寬、速度多變的環(huán)境噪聲及耦合噪聲、電噪聲等[1]。但若從與檢波器的關系來看,被記錄到磁帶上的噪聲可劃分為以下兩類。

    (1)與檢波器無關的噪聲。指即使檢波器不存在,也獨立客觀存在的噪聲。如表現為地面振動的環(huán)境噪聲、源自炮點的原生噪聲、由震源激發(fā)的次生噪聲等(圖1a)。

    (2)與檢波器有關的噪聲。這種噪聲決定于檢波器、地震儀的性能參數以及檢波器與地表介質之間的耦合關系。可更細致地劃分為: ①機械噪聲,包括檢波器濾波效應、寄生振蕩、耦合效應等,決定于檢波器以及檢波器—大地耦合系統(tǒng)的機械特性; ②電噪聲,主要指地震儀前置放大器等產生的噪聲以及有源檢波器元器件產生的噪聲;這兩類可歸屬為物理噪聲; ③數學噪聲,主要指采樣噪聲以及各種數學運算產生的噪聲,比較微弱(圖1b)。圖1c代表經過地震檢波系統(tǒng)改造后、記錄到磁帶上的數據,是信號與各種噪聲的累加。

    圖1 磁帶數據中信號與噪聲的影響因素(振幅譜)

    有效反射信號強度較噪聲強度的分貝數被稱為信號的動態(tài)范圍。根據筆者調查,山東東營HJ地區(qū)冬季微風夜晚最淺目的層反射波較環(huán)境噪聲的動態(tài)范圍大約為40dB,更深目的層則更小,部分文獻[3,4]也印證了這個量級。同時,震源激發(fā)后會產生次生噪聲,使得有效信號較噪聲的動態(tài)范圍變得更小。在復雜地表區(qū),無論環(huán)境噪聲還是次生噪聲都較普通地表更為強烈。信號較環(huán)境噪聲的動態(tài)范圍經過后續(xù)處理后會有一定程度的提高。

    對于陸上石油勘探而言,由人工震源激發(fā)的反射波信號較噪聲具有下列兩個特點。

    (1)主頻低。油氣藏單元一般埋藏較深,大多在數百米至數千米以下,反射波雙程旅行時多在1s以上。由于吸收、衰減等因素影響,高頻衰減嚴重,最高頻率多在150~200Hz以下;在低降速帶較厚、目的層較深地區(qū),反射波主頻更低。

    (2)信噪比低?;谕瑯拥脑?,反射波振幅在傳播過程中大幅度衰減,環(huán)境噪聲、原生噪聲及次生噪聲的強度變得相對較高。

    Davis[5]認為,接收系統(tǒng)設計中所采用的主要參數均來自對噪聲的研究,而噪聲是由信號來定義的。因此,無論是接收設備的硬件性能指標參數設計,還是震源激發(fā)、觀測系統(tǒng)、地面接收等參數的選擇,均基于對石油勘探中信號與噪聲特性的認識。對信號(尤其是弱小信號)的突顯和和對噪聲的壓制,貫穿地震勘探的每一個環(huán)節(jié),是很多技術、施工手段的出發(fā)點與最終目的[2]。

    3 檢波器指標參數對地震數據的影響

    對地震勘探野外采集來說,受當前技術手段的限制,無法直接測量到震源帶來的大地振動,所以只好將具有機電轉換能力的檢波器置入地面、間接測量大地機械振動,所以檢波器本身的性能指標對于保證地震數據的質量具有重要意義。

    文獻[6]曾經對兩種典型檢波器(動圈式模擬檢波器20DX和MEMS數字檢波器DSU3,表1)的性能指標對地震數據的影響做過討論,本文在此基礎上做更深入的分析、探討。

    表1 DSU3數字檢波器與20DX模擬檢波器的性能比較[7]

    3.1 跟蹤物理量

    根據跟蹤物理量的不同,可將檢波器分為位移型檢波器、速度型檢波器和加速度型檢波器三類。石油勘探中多為速度型或者加速度型。

    速度型檢波器的典型代表是動圈式檢波器,如SM-4或20DX。在MEMS數字檢波器出現之前,絕大多數石油勘探的地震數據都是由動圈式檢波器采集的。相應地,很多處理手段的假設條件以及研究結論,都是建立在速度域數據基礎上的,如雷克子波、隨機噪聲、信噪比計算公式等,都是基于數據的速度域特征[8,9]。

    當前加速度型檢波器的典型代表是MEMS數字檢波器,比如DSU和Vectorseis。在MEMS檢波器被引進中國之初,很多文獻及實際數據處理都將其輸出的數據與動圈式檢波器數據進行直接比較,沒有注意到二者跟蹤物理量的不同,所以均得出了“MEMS檢波器主頻高而信噪比稍低”的結論[7,10-16]。造成這一現象的根本原因是地震信號占據低頻端而加速度數據具有高頻提升的特點,并不是MEMS檢波器的本來屬性。將加速度型檢波器數據較速度型檢波器數據的主頻提高、頻帶拓展歸因于檢波器本身是不正確的,其實本質是跟蹤物理量不同或者量綱不同,與檢波器本身無關。近年來,絕大多數業(yè)內人士趨向認同此觀點,因此對MEMS等加速度檢波器的看法也更加客觀了[6,17-22]。

    文獻[23]表明:在機電比相對較高的前提下,經檢波器反褶積補償動圈式檢波器的低頻損失以后,以上兩類典型檢波器拾取的速度或者加速度在陸上石油勘探的頻帶和能量范圍內相互轉換后具有高度一致性。因為二者的差異僅僅是其跟蹤物理量(速度/加速度)的不同,它們表征的地表振動的物理存在并沒有本質的差異。同時,由檢波器反褶積補償動圈式檢波器的低頻成分或者將MEMS數字檢波器積分為速度后,低頻端(1~10Hz)數據得到了加強,這對于全波形反演、速度分析、深層成像等領域具有重要意義[24-27]。另外,可對以前用動圈式模擬檢波器采集的地震數據先做檢波器反褶積、補償低頻成分,再重新進行后續(xù)處理、解釋,可望獲得新的認識和勘探成果。

    至于到底是速度表征還是加速度表征更有利于石油地震勘探,需要由后續(xù)處理、解釋階段的勘探效果來驗證。

    3.2 絕對動態(tài)范圍

    在AD轉換器、前置放大器等硬件因素確定的情況下,可以采集到的地震信號的動態(tài)范圍是確定的,即電噪聲~AD轉換器最大量程,所有超過這個范圍的信號都不能被識別或者接收(此表述未考慮處理對信噪比的提高作用)。因為該范圍不會隨著輸入信號強度的變化而變化,所以作者將其稱為“絕對”動態(tài)范圍[6]。

    在現今超過120dB絕對動態(tài)范圍的機電背景下,前置放大器為主的電噪聲已經不是主要矛盾,因為與檢波器無關的機械噪聲的強度非常大,是電噪聲的十幾倍甚至幾十倍(圖2顯示中國東部平原地區(qū)冬天夜晚微風天氣時環(huán)境噪聲與電噪聲的對比),并且往往具有相干性(盡管有時以隨機形式表現出來)及空變與時變的特征,很難完全衰減。只有在處理階段有能力大幅度衰減機械噪聲的前提下,降低電噪聲對信噪比的促進作用才能顯現出來,這在現階段是不現實的。所以一味地強調降低電噪聲對提高數據信噪比的幫助并不大,如文獻[28]提到DSU508將本底噪聲降低了-10dB(較DSU428)。就檢波器的機電轉換性能來講,這確實是一個大的進步,但是對實際地震數據信噪比的促進作用非常有限。同時,電噪聲在主要頻段(如DSU的10~200Hz, FDU的3~200Hz)往往接近于白噪,僅靠統(tǒng)計性的疊加效應就可大幅度衰減。比如100次疊加,即可將電噪聲衰減至原強度的1/10以下。

    圖2 環(huán)境噪聲、DSU428電噪聲和新低噪聲模型的功率譜密度(PSD, 1m2·s-4·Hz-1=0dB)

    當然,降低電噪聲會對提高極低頻(1~3Hz及以下)弱信號的機電比有利,但在多次覆蓋的情況下,由普通炸藥或可控震源激發(fā)的極低頻信號可具有較高的機電比,不會對低頻恢復造成太大影響。

    文獻[28]提到了電噪聲的另一特點:越靠近低頻端(1~10Hz),本底噪聲越大,這是因為本底噪聲往往具有1/f規(guī)律(電噪聲強度與頻率成反比)造成的。在這種情況下,因為DSU本身是加速度型檢波器(意味著低頻機械信號較弱),加之電子器件1/f噪聲的存在,就會使DSU系列檢波器的低頻接收能力降低(機電比低)。但是據筆者的研究[23]表明,對于以石油勘探為目標的機械信號,DSU系列是可以滿足要求的。但是對于天然地震中的極低頻以及新低噪聲模型(NLNM)水平的振動,則無法檢測(圖2)。

    同時,對于有源檢波器而言,其輸入到模數轉換器件之前的模擬信號經過了外設電路的改造,在此過程中就產生了電噪聲,進而也就存在一個檢波器自身的絕對動態(tài)范圍。由于制造工藝方面的原因,有源檢波器的動態(tài)范圍很難超過現有地震儀的動態(tài)范圍(大約130dB)。如果有源檢波器的動態(tài)范圍很小,即電噪聲水平很高,就很難輸出高保真的機械信號、忠實地再現地面振動特性。

    比如對于圖3左圖中有源檢波器A而言,其輸出模擬信號的動態(tài)范圍略小于地震儀的動態(tài)范圍,在這種情況下幾乎能夠充分利用地震儀的動態(tài)范圍,是最理想的一種情況;對于圖3右圖中有源檢波器B而言,由于各種原因導致其輸出原始模擬信號的動態(tài)范圍遠小于130dB(電噪聲太強),此時即使將地震儀增益設置為12dB,從地震記錄上也會出現超調現象、讓人誤以為其輸出數據的動態(tài)范圍很大,但是因為其輸出模擬信號的本底電噪聲是固定的,所以對于微弱信號的識別能力并沒有隨著增益由0調整為12dB而提高(盡管此時地震儀的電噪聲水平降低了),數據的實際動態(tài)范圍仍然很小。在無源檢波器情況下能夠識別的微弱信號,在用有源檢波器B接收時,就變得不可識別。 所以,有源檢波器輸出模擬信號的動態(tài)范圍應該在110~130dB,如果太小,就會浪費地震儀的動態(tài)范圍,降低對微弱信號的識別能力。

    圖3 有源檢波器動態(tài)范圍應略小于地震儀動態(tài)范圍

    3.3 相對動態(tài)范圍

    包括20DX在內的多數模擬動圈式檢波器的諧波失真在-60dB左右,而數字檢波器DSU的諧波失真則為-90dB,相應地20DX的動態(tài)范圍為60dB,DSU的動態(tài)范圍為90dB(均指檢波器自身的動態(tài)范圍)。但是應該看到,盡管與前述絕對動態(tài)范圍都屬動態(tài)范圍的范疇,并且單位都是dB,但是由諧波失真定義的檢波器自身動態(tài)范圍與由AD轉換器轉換位數、電噪聲水平所定義的絕對動態(tài)范圍有著很大的差異,其計算公式、決定因素、物理含義均完全不同;如果根據檢波器的動態(tài)范圍60dB遠遠小于地震儀絕對動態(tài)范圍130dB就認為前者限制了后者發(fā)揮作用,是片面的。所以,為了與由AD轉換器轉換位數、電噪聲水平定義的絕對動態(tài)范圍相區(qū)別,作者將由諧波失真定義的動態(tài)范圍稱為相對動態(tài)范圍,因為這個動態(tài)范圍是隨著輸入信號的強度而浮動的,并非固定不變[6]。

    對于石油勘探而言,動圈式檢波器-60dB的諧波畸變即60dB的相對動態(tài)范圍目前是適用的,新型MEMS數字檢波器-90dB的諧波畸變即90dB的相對動態(tài)范圍僅具有實驗室意義,不會表現出明顯的信噪比改善。這一方面是因為相對動態(tài)范圍是隨著輸入信號的幅度而變化的,另外一方面仍然是因為與檢波器無關的機械噪聲的強度太大了[6]。

    3.4 檢波器靈敏度與前放增益

    從機電轉換角度來講,提高信號輸入幅度可采取提高檢波器靈敏度或地震儀前放增益兩種方式。

    針對不同的檢波器,其靈敏度的單位是不同的。比如動圈式速度檢波器的靈敏度單位是V/m/s,MEMS數字加速度檢波器的靈敏度單位是mV/m/s2。在這種情況下要比較靈敏度的高低,一方面要看其最強信號與A/D轉換器最大量程的差距,距離越小,靈敏度越高;另一方面要看電噪聲的水平,電噪聲水平越低而靈敏度越高,保真度(機電比)越高。

    比如,有的檢波器在地震儀適用12dB增益的情況下出現超調現象,并不代表其靈敏度就較其他檢波器更高,靈敏度高低的衡量應該基于相同增益。同時,如果靈敏度的提高是由外設電路實現的,就會不可避免地帶來附加的電噪聲。如果為了充分利用動態(tài)范圍而通過外設電路將靈敏度提高n倍,但隨之而來的是電噪聲提高了等于或者大于n倍,這種靈敏度提升是沒有意義的。很多有源檢波器均存在這個問題。

    所以,基于無源檢波器的靈敏度提高或者地震儀適當加大增益等措施是有意義的,有助于提高高、低頻兩端弱信號的機電比[23],提高信號保真度。從施工便利的角度來講,檢波器串聯(lián)或者提高增益(比如東部強信號地區(qū)用“單點檢波器+12dB增益”,西部弱信號地區(qū)用“單點檢波器+24dB增益”)更具有操作性。

    當然,也不必追求過高的機電比,比如機電比為100/1或1000/1幾乎是無區(qū)別的,此時去噪的主要矛盾是機械噪聲、不再是電噪聲。

    有文獻提出用低靈敏度的高頻部分與高靈敏度檢波器的高頻部分求取匹配濾波算子,然后將進行匹配濾波后的低靈敏度檢波器高頻部分與其低頻部分進行數據重構,以提高數據主頻[29]。這種方式相當于在數據的不同頻段進行不等權組合,已經偏離了地表振動這一源頭數據的頻率屬性,看似提高了主頻,實際沒有物理意義,產生了人為的畸變。

    提高靈敏度的確有利于弱信號的拾取以及數學恢復,但是在野外記錄以及初疊加剖面上難以體現(如單個20DX和一串20DX的初疊加剖面是相似的),須經后續(xù)數學、物理屬性分析才能展現出來。

    3.5 振幅和相位響應

    檢波器的振幅和相位響應主要與檢波器的主頻、阻尼以及寄生振蕩頻率有關。

    動圈式檢波器(比如20DX)與MEMS數字檢波器(比如DSU)具有不同的濾波響應以及各自明確的數學表達式。20DX的主頻為10Hz,DSU的主頻則超過1000Hz。對于石油勘探的頻段而言,前者振幅、相位在低頻段有較大畸變,后者的振幅、相位畸變都很小。但是,因為動圈式檢波器頻率響應是已知的,所以經過“檢波器反褶積”對振幅、相位進行補償后,二者在石油勘探的主要頻段內具有高度一致性[23]。

    動圈式檢波器的寄生振蕩頻率多為其自然頻率的10~15倍以上,所以在多數情況下不必考慮。

    就動圈式檢波器而言,自然頻率太低或太高都不應被提倡。若自然頻率太低(如2.0、4.5Hz),一方面會使得寄生振蕩頻率進入勘探頻段,另一方面會導致費用提高、設備壽命降低;若自然頻率太高(如40、60、100Hz),數據低頻部分會被大大衰減,機電比大大降低,信號畸變嚴重,即使經過檢波器反褶積等處理,仍難以恢復低頻(1~10Hz)段。

    文獻[30]提出用“低頻檢波器的低頻部分+高頻檢波器的高頻部分”進行數據重構,以便拓寬數據頻帶。實際上,二者頻帶上的視覺差異主要是由于振幅譜顯示時采用最大值作為0dB以及靈敏度不同導致的,高頻段歸一化后可看到高頻檢波器缺失低頻部分,而低頻檢波器包含了高頻檢波器的所有信息(圖4)。另外文獻[31]提出用“速度檢波器的低頻部分+加速度檢波器的高頻部分”進行拓頻,因為速度與加速度有著完全不同的物理含義,合并后數據完全偏離了地面振動這一數據源頭,失去了其原有的物理意義。文獻[32,33]曾經嘗試用“高低頻聯(lián)測”的方式進行低頻恢復,即以低頻檢波器的低頻部分作為約束、恢復高頻檢波器的被機械衰減的低頻部分,這從理論上是可行的。但是在實踐中,往往會因為兩種檢波器在空間位置、跟蹤物理量、耦合響應、電噪聲、環(huán)境噪聲等方面的差異,使得校正算子難以保持穩(wěn)定。

    圖4 低頻檢波器包含高頻檢波器信息

    綜上所述,10Hz動圈式檢波器是一個數據、工程意義均較理想的檢波器(須配合檢波器反褶積以恢復低頻)。

    3.6 耦合效應

    在石油地震勘探中,震源被激發(fā)后,反射波傳播到地下反射界面后再反射回地表,帶來地表介質的振動。受當前技術條件限制,無法直接測量大地振動,只能用檢波器置于地表介質間接測量地表的振動。檢波器與大地之間的柔性接觸,構成了一個振動系統(tǒng),被稱為檢波器—大地耦合振動系統(tǒng)。

    借助測耦檢波器技術(專利號:ZL 2015 2 0258826.8),可以測得不同檢波器類型、不同耦合情況下的檢波器—大地耦合響應。圖5為東營HJ地區(qū)普通泥質地表下測得的DSU和20DX的耦合效應。由圖5可見,在低頻部分,兩類檢波器的耦合效應基本一致; 但在160Hz以上,數字檢波器DSU的耦合效應要遠高于模擬檢波器20DX(即放大系數更大)。也就是說,在輸入相同機械振動(環(huán)境噪聲+地震信號)時,數字檢波器的高頻響應會更強,但這是一種高頻畸變,而非高頻提升。

    圖5 據實際試驗數據擬合的兩類檢波器放大系數曲線

    在檢波器—大地耦合響應方面,作者的主要認識包括[34]:

    (1)在檢波器被穩(wěn)固埋置的情況下,用單自由度有阻尼的振動系統(tǒng)可以描述耦合介質具有一定厚度時的檢波器—大地耦合響應;耦合響應的改善,可表現為耦合固有主頻與耦合阻尼比的增加;

    (2)耦合響應降低了地震信號的信噪比和分辨率,通過耦合反褶積可改善數據,提高信號保真度;

    (3)挖坑埋置具有降低高頻段耦合噪聲的作用;

    (4)耦合響應的高頻放大作用以及不一致性,使得檢波器自身優(yōu)秀性能指標的作用降低了,在某些情況下難以轉化為高質量的地球物理數據;

    (5)更小、更輕便、更利于耦合的設計有利于降低檢波器—大地耦合系統(tǒng)的非線性響應。

    3.7 允差

    檢波器各特性參數的實際值與標稱值之間可接受的誤差比率范圍稱為允差。通常在用的多數動圈式模擬檢波器的允差為5%,而允差為2.5%的檢波器則被稱為超級檢波器或者高精度檢波器。模擬數據和實際數據的計算結果[23]證明,在當前地震勘探背景(多次覆蓋、強噪聲、耦合條件差異大且極弱信號識別能力有限)下,高精度/超級檢波器在組合疊加后不會表現為明顯的、可感知的更高質量數據。

    3.8 對電磁干擾的響應與數字化

    文獻[7]認為,應用MEMS技術的數字檢波器DSU不再有任何連接到地震道的電感線圈,所以也就不再受任何電磁干擾信號的影響;但是,實踐證明這個結論并不全面。筆者曾經在一個試驗中發(fā)現部分DSU檢波器也存在某種“高頻干擾”(試驗中采用了36個數字檢波器與36個模擬檢波器進行對比),并且這種高頻干擾的幅度非常大。在總數同為36個的檢波器集合中,存在高頻干擾數字檢波器較模擬檢波器要高得多(DSU:12~13個,20DX:2~3個)。在不同的時間、地點以同樣的方法進行了同樣的試驗,這種高頻干擾卻沒有出現。因此,這種高頻干擾可能是外界干擾與檢波器內部因素相互作用的結果,并且是與電磁干擾而非機械干擾相關的噪聲; 部分野外地震采集中出現的工頻干擾也進一步證實了以上猜測。以上現象進一步說明了DSU檢波器難以對電磁干擾完全“免疫”,并且這種干擾無疑會對地震弱信號的接收產生影響。當然,DSU減少模擬電路比重的做法無疑會有助于減小電磁干擾,但顯然不能完全消除電磁干擾[6]。

    因此,在檢波器內部進行數字化有助于減少天電、工頻干擾,但這并不能帶來數據根本性的改變,因為任何數字信號最初都來源于模擬信號。

    因此筆者認為,在當前去噪能力下,檢波器數字化并不能帶來地震數據根本性的變革,進而體現出較模擬數據更大的優(yōu)勢。數字化的主要意義是工程意義,不是數據意義。藉由檢波器的數字化,進而可實現: ①“小型化”,其目的是衰減耦合效應,降低埋置難度[21,34,35]; ②“無纜化”,旨在通過無線傳輸(檢波器→采集站)避免由沉重的外接電纜引起的耦合響應畸變(外接電纜使得耦合響應不再服從單自由度機械振動系統(tǒng)的數學描述)[21,34,35]; ③“自動化”,其目的在于在以上“三化”基礎上實現檢波器由人工埋置向機械埋置轉變,在確保耦合質量的基礎上,提高施工效率、減少用工總量、節(jié)約生產成本。

    3.9 矢量保真度

    更高的矢量保真度在處理階段擁有優(yōu)勢,但其應用效果方面的結論尚不統(tǒng)一。筆者認為其原因仍是噪聲強度。在信噪比較低地區(qū),矢量保真度提高的優(yōu)勢被噪聲的影響淹沒了[6]。

    3.10 非線性噪聲

    提高制造工藝水平會有效降低非線性噪聲。

    3.11 地震儀因素

    包括前放增益、采樣噪聲、高低截止頻率等,影響較小。

    3.12 組合方式

    檢波器選型問題往往伴隨著組合形式的討論。當前有一種呼聲,認為可以用一種“好的單點檢波器”代替某種“不好的已有檢波器組合”。筆者認為,如果不考慮檢波器濾波效應以及寄生振蕩(前者可以補償[23],后者多超出了地震勘探的主要頻帶范圍),改進檢波器性能主要壓制的是電噪聲,即機電轉換過程中產生的噪聲。而具有相干或者隨機特征的環(huán)境噪聲、原生以及次生噪聲屬于“與檢波器無關的機械噪聲”(圖1),其強度為電噪聲的十幾倍甚至幾十倍(圖2),即使沒有檢波器的存在,它也是客觀存在的,不會因為檢波器性能的改善而得到衰減,必須通過檢波器的方向性效應或者統(tǒng)計性效應以及后續(xù)的處理手段進行衰減,無法通過提高檢波器性能(即降低電噪聲。當前施工中占據主導地位的20DX+FDU采集單元的電噪聲已經非常低)來達到衰減機械噪聲、減少檢波器數量的目的(無法壓制具有空間特征的機械噪聲)。

    所以,如果某種檢波器的電噪聲較20DX+FDU組合稍高或者稍低,數據表現并不會有太大的差異。因為目前的去噪能力有限,去噪的主要矛盾是強度大、頻帶寬、具有相干特征的機械噪聲,而不是由檢波器、地震儀所決定的電噪聲。

    所以,“用某一種指標稍好的檢波器(單個)可代替多個某種指標稍差的檢波器(比如20DX)”是不現實的,是個“偽命題”。即使在某些地區(qū)經由此種檢波器、普通采樣密度的情況下達成了地質目標,也只能說明一個問題:單個20DX+低頻補償可以取得基本等同的接收效果(此類地區(qū)往往次生噪聲較弱)。

    組合還是單點主要取決于組合壓噪能力,參看文獻[36-43]。

    另外,關于組合方式的認識還包括以下幾點。

    (1)10~12個檢波器是當前野外施工的適宜選擇,可兼顧數據(壓噪)與工程(施工方便)的需要[37]。

    (2)單點采集具有頻帶寬、室內自由組合的優(yōu)勢,可以有效地提高地震資料的分辨率。但是,在有些地區(qū)并不能完全摒棄野外組合,特別是在次生干擾非常嚴重、信噪比較低的地區(qū),單純依靠室內手段難以達到有效衰減次生干擾的目的。

    (3)隨著軟硬件技術高速發(fā)展、地震儀帶道能力大幅提升、數據處理能力快速增長,“單點不組合+高采樣密度”是今后發(fā)展的方向。

    綜上所述,地震檢波器更好的性能指標肯定有助于提高數據質量,但是更高的性能指標并不能保證更好的地球物理效果。復雜地質條件導致的地震信號自身的差異、采集階段的信噪比以及處理階段的噪聲衰減能力起決定性作用。電學性能指標的提高必須在原始信噪比較高、去噪能力持續(xù)改進的前提下才能顯現作用。在地震勘探其他環(huán)節(jié)不銜接的情況下,單純追求非常小的失真或者非常大的動態(tài)范圍等電學指標的提高,往往難以表現為地震數據質量的系統(tǒng)性進步[6]。

    4 檢波器的選型與資料對比方法

    如果主要目的在于確定“適于工業(yè)應用的檢波器類型”,按照科學實驗中“唯一變量”的原則,結合勘探現實,應該采取以下步驟。

    (1)研究不同類型檢波器的性能指標參數。優(yōu)良的性能指標是保證數據記錄質量的基礎。

    (2)野外單點對單點試驗,即單個檢波器與單個檢波器比較,二者相距20~30cm是合理的距離。在這個范圍內,耦合響應較小,地震信號、機械噪聲的相似度較高,具有可比性。如果距離太大,可比性就會變差。

    (3)如果想檢驗某一種新型檢波器,應取兩種以上成熟檢波器作為標準檢波器,如20DX、DSU系列。

    (4)將不同類型檢波器(速度型、加速度型)校正到一個域,即將量綱統(tǒng)一為速度或者加速度(過程中主要涉及到檢波器的靈敏度以及微積分計算)。

    (5)根據步驟(1)中提供的檢波器頻率響應校正低頻,參照文獻[23]。

    (6)在單道、單炮、剖面等不同層面上進行信噪比、能量、子波等數學物理數據屬性比較,不要僅根據單炮或者剖面的“面貌”下結論。

    評價一個檢波器的好壞,不可避免地要涉及到評價指標的選擇。當前在野外采集階段評價記錄剖面的時候,常采用頻率掃描等五種方式或指標,但是,對于檢波器選型試驗中的資料對比而言: ①頻率掃描,對于不在同一個域(速度/加速度)的數據來說,頻率掃描的結果對比沒有意義,因為數據基礎不同; ②頻譜分析,現在普遍采用的三種顯示方式(Raw Power,Percent,dB)的情況下,會出現誤導、錯誤,需要根據不同的試驗內容,調整顯示方式,或采用更具有針對性的顯示方式; ③信噪比,從原理上講,“加速度檢波器輸出數據較速度檢波器數據的信噪比低”的現象是由地震信號占據低頻段這一特性決定的,與檢波器類型沒有任何關系,也不是評判一個檢波器好壞的依據;信噪比指標僅僅是一個非量化指標,在數值非常接近的時候,其參考意義有限,不應被作為一個絕對的標準;用加速度、速度檢波器輸出數據進行直接的信噪比比較是不合理的; ④能量,能量對比在不同施工因素(包括工區(qū)、構造、低降速帶、增益、靈敏度、記錄格式等)的情況下沒有可比性,沒有參考意義; ⑤子波,同樣地,在沒有可比性的前提下,采集分析軟件中基于自相關理論的子波分析的表征作用是相對的、有限的。

    (7)得到數據處理特別是地質綜合解釋人員的認可。無論在地震采集、數據處理階段展現多大的優(yōu)勢,提高地質解釋效果才是最終目的。

    5 對幾種典型檢波器的初步認識

    在石油地震勘探中,檢波器的作用是以盡量小的失真產生地面振動單分量或者多分量的電模擬,完整地反映地震波的動力學特征。在這種間接測量過程中,有兩個環(huán)節(jié)上產生了誤差:大地振動與檢波器外殼振動之間的差異——耦合效應; 檢波器外殼振動與地震數據之間的差異——機電效應。所以,檢波器的設計與制造必須最大限度地減小耦合效應帶來的耦合噪聲以及機電效應帶來的電噪聲,盡量真實地記錄大地振動?!皬驮钡孛嬲駝邮钦嬲哪繕耍苷鎸嵱涗浀乇碚駝拥木褪恰昂谩睓z波器。

    評價檢波器的重要標準是看高低頻兩端(從對比的角度來講,高機電比的中頻段應該是一樣的),低頻端與機械濾波響應以及地震儀動態(tài)范圍有關,高頻端與耦合響應以及地震儀的動態(tài)范圍有關。

    基于對陸地石油勘探中信號噪聲的相對特征以及檢波器性能參數方面的研究,獲得以下認識。

    (1)MEMS(如DSU,Vectorseis)數字檢波器是“嚴謹”、“可靠”的檢波器。今后該類型檢波器改進的重點是大幅降低產品成本、深入研究速度域低頻信息強化后的數據處理及石油地質意義、縮小體積以提高耦合能力等方面。該類型檢波器的主要優(yōu)勢在速度域的低頻、而不是加速度域的高頻。

    但需要特別注意的是,以MEMS為代表的各類加速度檢波器輸出的數據較速度檢波器輸出的數據呈現出來的更高主頻是一種“偽高頻”(數學高頻),并不代表其所表征的地表機械振動的主頻真的被提高了(物理高頻),本質是由表征數據的量綱(加速度/速度)差異導致的。比如將1m高的桌子表達為100cm(單位不同),數值增加了,但并不意味著桌子物理高度的增加。這一點在很多處理剖面中被認為是加速度檢波器較速度檢波器的重大進步,是錯誤的。就像認為100cm>1m一樣,是沒有認清數據所對應的量綱以及所代表的物理意義所導致的。

    如果用加速度檢波器接收的剖面質量確有進步的話,其主要原因是極高的空間采樣密度、更合理的觀測系統(tǒng)屬性以及處理手段的改進。

    當然,如果加速度數據較速度數據對地下介質的空間、物性特征具有更好的表征作用,則另當別論。但是,加速度表征不一定意味著一定要采用加速度型檢波器,因為對于成熟的速度/加速度檢波器而言,其數據經由微積分在陸上石油勘探的興趣頻帶(1~200Hz)范圍內是可以互換的。所以,到底是速度表征還是加速度表征更有利于石油勘探,仍有待從地震數據的地球物理屬性角度進行更多的理論分析、試驗驗證,不能經由簡單的生產結果比對進行判斷。

    (2)在當前去噪能力下,20DX等動圈式檢波器串仍然是高效率檢波器:數據合格(合理的機電比),施工方便,費用合理。同時經由檢波器反褶積可以有效補償該類檢波器的低頻衰減,取得與MEMS數字檢波器相當的數據質量[23];也可以改進為測耦檢波器,使檢波器具有獨立測量耦合響應的能力,在加強現場耦合質量監(jiān)控的同時,消除耦合效應、提高高頻弱信號的同相性。

    即使對于高密度采集而言,一串動圈式檢波器的極小面積組合(檢波器分散在1~2m以內)+低頻補償也是較好的選擇(檢波器“開會”會因為小線的干擾產生耦合畸變以及額外的噪聲干擾)。因為目前多數施工單位擁有大量此類檢波器,可以在僅增加少量生產成本的情況下,實現非常高的空間采樣密度。這種采集方式可被稱為“基于動圈式檢波器串和低頻補償的高密度采集”。

    (3)高靈敏度或者高增益是有意義的,但是必須以機電比的提高為前提,檢波器串可較方便地實現這一目的。經由外設電路提高檢波器靈敏度的方式必須以機電比的提高為前提,否則只會降低微弱信號識別能力。當前很多新型檢波器的努力方向都是改進其機電效能,之所以沒有取代20DX在采集中的主導地位,其根本原因是石油勘探中“與檢波器無關的機械噪聲”太強了,由提高靈敏度凸顯出來的微弱信號被強烈的機械噪聲淹沒了。

    (4)低頻檢波器、高頻檢波器、高精度(超級)檢波器不會表現出明顯的數據優(yōu)勢(相對于20DX+低頻補償)。

    6 結束語

    (1)應科學合理地設計檢波器選型試驗,充分做好理論研究、模擬,再進行點試驗,進而做規(guī)模試驗。不宜在未厘清機理或對比不合理的情況下,開展大規(guī)模試驗或采集施工,造成資源的浪費。

    (2)檢波器選型試驗涉及到機電、采集、處理等諸多環(huán)節(jié),宜仔細梳理、系統(tǒng)地看問題,在嚴謹的數學、物理、地質屬性分析基礎上,根據最終的地質解釋效果做出判斷。特別是不能只憑采集階段數據的“面貌”就草率地下結論。

    (3)經由外設電路提高檢波器靈敏度的方式必須以機電比的提高為前提,否則只會降低對微弱信號的識別能力。

    (4)數字化、小型化、無纜化、自動化是檢波器的發(fā)展方向。

    [1] Julien M.Seismic Acquisition from Yesterday to To-morrow.Distinguished Instructor Short Course,Sponsored by SEG and EAGE,2011.

    [2] 魏繼東.石油勘探地震檢波.山東青島:中國海洋大學出版社,2015.

    Wei Jidong.Seismic Sensing in Geophysical Prospecting.Ocean University of China Press,Qingdao,Shandong,2015.

    [3] 李慶忠.地震高分辨勘探中的誤區(qū)與對策.石油地球物理勘探,1997,32(4):751-783.

    Li Qingzhong.Some mistaken concepts in high-resolution seismic exploration and the corresponding countermeasures.OGP,1997,32(4):751-783.

    [4] Crisi P A,Perrin T J.How much wind is enough? SEG Technical Program Expanded Abstracts,2003,22.

    [5] Davis B K.Expert system for designing geophone and source arrays for onshore seismic surveys.SEG Technical Program Expanded Abstracts,1988,7:488-490.

    [6] 魏繼東.地震檢波器性能指標與地球物理效果分析.石油物探,2013,52(3):265-274.

    Wei Jidong.Geophysical requirements for geophone in petroleum seismic prospecting.GPP,2013,52(3):265-274.

    [7] 張丙和,崔樵,裴云廣.新型三分量數字檢波器DSU.石油儀器,2005,19(4):39-44.

    Zhang Binghe,Cui Qiao,Pei Yunguang.New three-component geophone-DSU.Petroleum Instruments,2005,19(4):39-44.

    [8] 張軍華,藏勝濤,周振曉等.地震資料信噪比定量計算及方法比較.石油地球物理勘探,2009,44(4):481-486.

    Zhang Junhua,Zang Shengtao,Zhou Zhenxiao et al.Quantitative computation and comparison of S/N ratio in seismic data.OGP,2009,44(4):481-486.

    [9] Ricker N.The form and nature of seismic waves and the structure of seismograms.Geophysics,1940,5(4):348-366.

    [10] Ronen S,Gibson J,Burnett R et al.Comparison of multi-component data from different MEMS sensors.EAGE 67thConference Exhibition,Madrid,Spain,2005.

    [11] Gibson J and Burnett R.Another Look at MEMS Sensors and Dynamic Range.CSEG Recorder,2005.

    [12] Ronen S,Comeaux L,Cartwright M et al.Comparison between geophones and two MEMS types and repeatability of land data.SEG Technical Program Expanded Abstracts,2005,24.

    [13] Stotter C,Angerer E and Herndler E.Comparison of single sensor 3C MEMS and conventional geophone arrays for deep target.SEG Technical Program Expanded Abstracts,2008,173-177.

    [14] Gibson J,Burnett R,Ronen S et al.MEMS sensors some issues for consideration.Geophysics, 2011,76(3):B79-B87.

    [15] Mougenot D,Cherepovskiy A and Liu J J.MEMS-based accelerometers expectations and practical achievements.First Break,2011,29(2):85-90.

    [16] 陳志德,關昕,李玲等.數字檢波器地震資料高保真寬頻帶處理技術.石油地球物理勘探,2012,47(1):46-55.

    Chen Zhide,Guan Xin,Li Ling et al.Digital detectors seismic data high fidelity wideband processing technology.OGP,2012,47(1):46-55.

    [17] Ronen S,Comeaux L,Cartwright M et al.Comparison between geophones and two MEMS types and repeatability of land data.SEG Technical Program Expanded Abstracts,2005,24.

    [18] Hons M,Stewart R,Lawton D et al.Ground motion through geophones and MEMS accelerometers:sensor comparison in theory,modeling and field data.SEG Technical Program Expanded Abstracts,2007,26.

    [19] Hauer G,Hons M,Stewart R et al.Field data compari-son:3C-2D data acquisition with geophones and accelerometers.SEG Technical Program Expanded Abstracts,2008,27.

    [20] 魏繼東.從一個野外試驗看不同檢波器的差異——以DSU3與20DX為例.地球物理學進展,2013,28(3):1398-1407.

    Wei Jidong.The difference between MEMS accelero-meter and analog geophone showed in a field test.Progress in Geophysics,2013,28(3):1398-1407.

    [21] 于富文,魏繼東,馮玉蘋.檢波器—大地耦合響應對地震數據的影響.地球物理學進展,2017,32(1):339-343.

    Yu Fuwen,Wei Jidong,Feng Yuping.How the ground-geophone coupling affects the seismic data.Progress in Geophysics,2017,32(1):339-343.

    [22] Wei Jidong.Comparing the MEMS accelerometer and the analog geophone.The Leading Edge,2013,32(10):1206-1010.

    [23] 魏繼東.檢波器反褶積對低頻信息的補償作用.石油地球物理勘探,2016,51(2):224-231.

    Wei Jidong.Geophone deconvolution low-frequency compensation for seismic data.OGP,2016,51(2):224-231.

    [24] Plessix R E,Baeten G,de Maag J W et al.Application of acoustic full waveform inversion to a low-frequency large-offset land data set.SEG Technical Program Expanded Abstracts,2010,29:930-934.

    [25] Castagna J P,Sun S,Siegfried R W.Instantaneous spectral analysis:Detection of low-frequency shadows associated with hydrocarbons.The Leading Edge,2003,22(2):120-127.

    [26] Tatalovic R P,Uldall A et al.Event-based low-frequency impedance modeling using well logs and seismic attributes.The Leading Edge,2008,27(5):592-603.

    [27] Kroode F T,Bergler S,Corsten C et al.Broadband seismic data — The importance of low frequencies.Geophysics,2013,78(2):WA3-WA14.

    [28] Lainé J,Mougenot D.A high-sensitivity MEMS-based accelerometer.The Leading Edge,2014,33(11):1234-1242.

    [29] 浙江大學.提高灰?guī)r出露區(qū)地震數據分辨率的分頻匹配方法.發(fā)明專利,申請?zhí)枺?01110152185.4.

    [30] Shen Honglei,Tian Gang,Shi Zhanjie.Partial frequency band match filtering based on high-sensitivity data: method and applications.Applied Geophysics,2013,10(1):15-24.

    [31] 嚴洪瑞,于承業(yè),周志才等.一種寬頻地震信號采集方法.CPS/SEG2014北京國際地球物理會議,148-151.

    Yan Hongrui,Yu Chengye,Zhou Zhicai et al.A wideband seismic data acquisition method.Expended Abstracts of CPS/SEG Beijing 2014 International Geophysical Conference,148-151.

    [32] Zhang Yimeng,Zou Zhihui and Zhou Huawei.Esti-mating and recovering the low-frequency signals in geophone data.SEG Fall Meeting,Las Vegas,2012.

    [33] Chiu S K,Eick P,Davidson M et al.The feasibility and value of low-frequency data collected using co-located 2Hz and 10Hz.SEG Technical Program Expanded Abstracts,2012,31.

    [34] 魏繼東.檢波器—大地耦合系統(tǒng)特性的振動力學解釋、模態(tài)參數識別及其對地震資料的影響與消除.地球物理學進展,2013,28(4):1983-1995.

    Wei Jidong.The vibration mechanics explanation for ground-geophone coupling effect,model identification and its impacting to seismic data and eliminating.Progress in Geophysics,2013,28(4):1983-1995.

    [35] 魏繼東.檢波器—大地耦合響應的野外規(guī)模測量與衰減方法.石油地球物理勘探,2017,52(3):411-417.

    Wei Jidong.Large-scale measuring and correction of geophone-ground coupling effect in field.OGP,2017,52(3):411-417.

    [36] 李慶忠,魏繼東.論檢波器橫向拉開組合的重要性.石油地球物理勘探,2007,43(4):375-382.

    Li Qingzhong,Wei Jidong.Concerning the importance of open horizontal geophone combination.OGP,2007,43(4):375-382.

    [37] 魏繼東.檢波器數量對組合壓噪能力的影響.物探與化探,2011,35(2):238-242,247.

    Wei Jidong.The effects of the number of geophones in array on noise attenuation.Geophysical and Geoche-mical Exploration,2011,35(2):238-242,247.

    [38] 李慶忠,魏繼東.高密度采集中組合效應對高頻截止頻率的影響.石油地球物理勘探,2007,42(4):363-369.

    Li Qingzhong,Wei Jidong.Influence of array effect on cutoff frequency of high frequency in high-density seismic acquisition.OGP,2007,42(4):363-369.

    [39] 魏繼東,丁偉.檢波器野外組合因素對地震資料品質的影響分析.石油物探,2010,49(3):312-318.

    Wei Jidong,Ding Wei.Impact analysis of geophone array factors on seismic data quality.GPP,2010,49(3):312-318.

    [40] 羅仁澤,梁黎明,吳希光等.寬線大組合理論及其在黃土塬地震采集中的應用.天然氣工業(yè),2009,29(2):54-56.

    Luo Renze,Liang Liming,Wu Xiguang et al.The theory of wide line array and its application in loess plateau seismic survey.Natural Gas Industry,2009,29(2):54-56.

    [41] 盧占國,呂景峰,劉新文等.復雜山地寬線大組合優(yōu)化應用.天然氣勘探與開發(fā),2009,32(1):18-25.

    Lu Zhanguo,Lv Jingfeng,Liu Xinwen et al.Complex mountain wide line big combinatorial optimization applications.Natural Gas Exploration and Development,2009,32(1):18-25.

    [42] 王棟,賀振華,孫建庫等.寬線加大基距組合技術在喀什北區(qū)塊復雜山地的應用.石油物探,2010,49(6):606-610.

    Wang Dong,He Zhenhua,Sun Jianku et al.More from the wide line combination technology in Kashgar north block in the hill country of complex applications.GPP,2010,49(6):606-610.

    [43] 朱鵬宇,楊晗,楊海濤等.寬線觀測大組合接收技術在阜康斷裂帶的應用.勘探地球物理進展,2010,33(5):359-362.

    Zhu Pengyu,Yang Han,Yang Haitao et al.Application of wide-line array receiving technique in Fukang Fault Zone.Progress in Exploration Geophysics,2010,33(5):359-362.

    *山東省東營市中國石化石油工程地球物理公司勝利分公司,257100。 Email:563064949@qq.com

    本文于2016年9月13日收到,最終修改稿于2017年9月24日收到。

    1000-7210(2017)06-1127-10

    魏繼東.適用于陸上石油勘探的地震檢波器.石油地球物理勘探,2017,52(6):1127-1136.

    P631

    A

    10.13810/j.cnki.issn.1000-7210.2017.06.001

    (本文編輯:朱漢東)

    魏繼東 高級工程師,1974年生;2008年獲中國海洋大學海洋地球物理學博士學位; 2008~2011年在中國石化勝利油田做博士后研究; 2012年為英國Imperial College London 訪問學者; 發(fā)表論文多篇、專著1部(《 石油勘探地震檢波》); 現在中國石化地球物理公司勝利分公司從事地震數據采集及信號處理方法研究。

    猜你喜歡
    動圈檢波器耦合
    非Lipschitz條件下超前帶跳倒向耦合隨機微分方程的Wong-Zakai逼近
    常規(guī)動圈式地震檢波器測試現狀及漏洞分析
    電子制作(2019年20期)2019-12-04 03:51:52
    檢波器容差對地震信號接收的影響研究
    600kN超大推力電磁振動試驗臺動圈結構的模態(tài)分析與優(yōu)化
    一種井下檢波器測試儀的設計
    電子測試(2018年6期)2018-05-09 07:31:34
    基于高低頻聯(lián)測的常規(guī)檢波器數據低頻振幅和相位同時恢復方法
    基于“殼-固”耦合方法模擬焊接裝配
    大型鑄鍛件(2015年5期)2015-12-16 11:43:20
    動圈式揚聲器低頻驅動下穩(wěn)態(tài)溫度場數值分析研究
    聲學技術(2014年2期)2014-06-21 06:59:06
    求解奇異攝動Volterra積分微分方程的LDG-CFEM耦合方法
    非線性耦合KdV方程組的精確解
    欧美久久黑人一区二区| 深夜精品福利| 最近最新免费中文字幕在线| 大码成人一级视频| 国产精品野战在线观看| 亚洲精品在线美女| 国产精品久久久人人做人人爽| 一夜夜www| 免费在线观看亚洲国产| 波多野结衣巨乳人妻| 国产一级毛片七仙女欲春2 | 叶爱在线成人免费视频播放| 黄色片一级片一级黄色片| 欧美久久黑人一区二区| 精品国产亚洲在线| 欧美激情 高清一区二区三区| 人成视频在线观看免费观看| 亚洲成国产人片在线观看| 波多野结衣巨乳人妻| 12—13女人毛片做爰片一| 91麻豆精品激情在线观看国产| 免费观看人在逋| 久久这里只有精品19| 久久青草综合色| 久久久久久免费高清国产稀缺| 国产主播在线观看一区二区| 久久精品91无色码中文字幕| 国产人伦9x9x在线观看| 黄网站色视频无遮挡免费观看| 男女做爰动态图高潮gif福利片 | 久久精品亚洲精品国产色婷小说| 巨乳人妻的诱惑在线观看| 亚洲九九香蕉| 精品欧美一区二区三区在线| 午夜激情av网站| 97人妻精品一区二区三区麻豆 | 母亲3免费完整高清在线观看| 中文字幕色久视频| 女性被躁到高潮视频| 亚洲美女黄片视频| 午夜影院日韩av| 亚洲第一av免费看| 国产91精品成人一区二区三区| 中文字幕人成人乱码亚洲影| 久久久久久人人人人人| 日韩中文字幕欧美一区二区| 不卡一级毛片| 美女大奶头视频| 午夜免费观看网址| 久久久久久大精品| 国产精品久久久久久人妻精品电影| 两性夫妻黄色片| 成人国语在线视频| xxx96com| 91九色精品人成在线观看| 一二三四在线观看免费中文在| 一进一出抽搐gif免费好疼| 精品久久久久久,| 久久精品国产清高在天天线| 国内毛片毛片毛片毛片毛片| 给我免费播放毛片高清在线观看| xxx96com| 日韩欧美国产一区二区入口| 免费女性裸体啪啪无遮挡网站| 日本撒尿小便嘘嘘汇集6| 热re99久久国产66热| 伊人久久大香线蕉亚洲五| 亚洲精品国产精品久久久不卡| 久久久国产成人免费| 国产成人av激情在线播放| 丝袜在线中文字幕| 国产一区二区三区综合在线观看| 69精品国产乱码久久久| 欧美乱码精品一区二区三区| 国产不卡一卡二| 一级片免费观看大全| 国产精品永久免费网站| 老司机靠b影院| 变态另类丝袜制服| 日本一区二区免费在线视频| 一本综合久久免费| 午夜a级毛片| 久久久久久久久久久久大奶| 国内久久婷婷六月综合欲色啪| 九色亚洲精品在线播放| 婷婷六月久久综合丁香| 老汉色av国产亚洲站长工具| or卡值多少钱| 免费久久久久久久精品成人欧美视频| 成人国产综合亚洲| 午夜福利一区二区在线看| 精品欧美一区二区三区在线| 少妇的丰满在线观看| 又黄又粗又硬又大视频| av欧美777| 亚洲精华国产精华精| 法律面前人人平等表现在哪些方面| 色播在线永久视频| 精品一品国产午夜福利视频| 人成视频在线观看免费观看| 操美女的视频在线观看| 亚洲天堂国产精品一区在线| 91精品国产国语对白视频| 窝窝影院91人妻| 精品卡一卡二卡四卡免费| 搡老熟女国产l中国老女人| 侵犯人妻中文字幕一二三四区| 999精品在线视频| 日韩精品免费视频一区二区三区| 亚洲av电影在线进入| 男人的好看免费观看在线视频 | 成人国产一区最新在线观看| 亚洲,欧美精品.| 欧美日韩乱码在线| 99香蕉大伊视频| 国产精品 欧美亚洲| 亚洲国产日韩欧美精品在线观看 | 日韩国内少妇激情av| 亚洲av成人不卡在线观看播放网| 两个人视频免费观看高清| √禁漫天堂资源中文www| 黑人巨大精品欧美一区二区mp4| 国产精品免费一区二区三区在线| 亚洲人成伊人成综合网2020| 久久久久国产精品人妻aⅴ院| 亚洲av成人不卡在线观看播放网| 成人免费观看视频高清| 午夜精品在线福利| 国产亚洲欧美精品永久| 精品一区二区三区av网在线观看| 熟妇人妻久久中文字幕3abv| 久久久久久免费高清国产稀缺| 久久婷婷人人爽人人干人人爱 | 免费观看精品视频网站| 亚洲精品一卡2卡三卡4卡5卡| 免费看a级黄色片| 欧美激情极品国产一区二区三区| 国产精品久久视频播放| 久久久国产成人精品二区| 国产1区2区3区精品| 欧美成人午夜精品| 婷婷六月久久综合丁香| 久久人妻熟女aⅴ| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久电影中文字幕| 一进一出抽搐gif免费好疼| 久久久精品欧美日韩精品| 久久香蕉激情| 国产精品永久免费网站| 身体一侧抽搐| 88av欧美| 男人操女人黄网站| 欧美 亚洲 国产 日韩一| 美女大奶头视频| 午夜福利在线观看吧| 亚洲情色 制服丝袜| 国产成人精品久久二区二区免费| 麻豆一二三区av精品| 变态另类丝袜制服| 两个人视频免费观看高清| 一区二区三区精品91| aaaaa片日本免费| 夜夜夜夜夜久久久久| 午夜久久久久精精品| 国产亚洲精品综合一区在线观看 | 香蕉国产在线看| 99久久综合精品五月天人人| 精品一区二区三区av网在线观看| 国产免费男女视频| 成人特级黄色片久久久久久久| 国产精品久久久av美女十八| 国产精品,欧美在线| 欧美日韩乱码在线| 国产av又大| 精品第一国产精品| 嫩草影视91久久| 一夜夜www| 久久影院123| 999久久久国产精品视频| 日韩欧美免费精品| 91老司机精品| 男女下面进入的视频免费午夜 | x7x7x7水蜜桃| 亚洲专区国产一区二区| 国产精品日韩av在线免费观看 | 国产一卡二卡三卡精品| 国产麻豆69| 美女国产高潮福利片在线看| 美女高潮到喷水免费观看| 欧美一级毛片孕妇| 国产精品综合久久久久久久免费 | 啦啦啦韩国在线观看视频| 日本五十路高清| 国产精品久久久久久人妻精品电影| 黄片播放在线免费| 国产精品香港三级国产av潘金莲| 精品熟女少妇八av免费久了| 亚洲午夜精品一区,二区,三区| 亚洲成av片中文字幕在线观看| 两人在一起打扑克的视频| 亚洲色图av天堂| 在线天堂中文资源库| 性欧美人与动物交配| 一级黄色大片毛片| 叶爱在线成人免费视频播放| 亚洲精品美女久久久久99蜜臀| 热99re8久久精品国产| 曰老女人黄片| 淫秽高清视频在线观看| 亚洲天堂国产精品一区在线| 女人高潮潮喷娇喘18禁视频| 日韩大码丰满熟妇| 亚洲最大成人中文| 久久久久久久午夜电影| 日本三级黄在线观看| 女人被躁到高潮嗷嗷叫费观| 日本 欧美在线| 欧美久久黑人一区二区| av欧美777| 久99久视频精品免费| 亚洲va日本ⅴa欧美va伊人久久| 国产私拍福利视频在线观看| 久久草成人影院| 精品国产亚洲在线| av视频在线观看入口| 我的亚洲天堂| 国产高清激情床上av| 黄色片一级片一级黄色片| 777久久人妻少妇嫩草av网站| 欧美激情极品国产一区二区三区| 制服人妻中文乱码| 日本欧美视频一区| 午夜精品久久久久久毛片777| 天天一区二区日本电影三级 | 好看av亚洲va欧美ⅴa在| 1024香蕉在线观看| 亚洲欧美激情综合另类| 日日夜夜操网爽| 欧美 亚洲 国产 日韩一| 精品人妻1区二区| 一级a爱片免费观看的视频| 中文字幕人妻丝袜一区二区| 亚洲精华国产精华精| 欧美一级毛片孕妇| 午夜福利影视在线免费观看| 此物有八面人人有两片| 老汉色av国产亚洲站长工具| 国产精品久久视频播放| 国产亚洲欧美精品永久| 欧美性长视频在线观看| 神马国产精品三级电影在线观看 | 欧美日韩精品网址| 国产免费av片在线观看野外av| 亚洲黑人精品在线| 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 国产成人av教育| 亚洲专区国产一区二区| 亚洲av熟女| 成人欧美大片| 久久天堂一区二区三区四区| 亚洲avbb在线观看| 亚洲自拍偷在线| 午夜精品在线福利| 亚洲熟女毛片儿| 香蕉丝袜av| 黄色a级毛片大全视频| 国产精品香港三级国产av潘金莲| 久久久国产成人精品二区| 悠悠久久av| 久久精品成人免费网站| 精品福利观看| 国产精品av久久久久免费| 久久人人爽av亚洲精品天堂| 男人舔女人下体高潮全视频| 最新在线观看一区二区三区| 黄色 视频免费看| 啪啪无遮挡十八禁网站| 少妇熟女aⅴ在线视频| 女生性感内裤真人,穿戴方法视频| 亚洲国产中文字幕在线视频| 女人高潮潮喷娇喘18禁视频| 美女 人体艺术 gogo| 国产亚洲av高清不卡| 日韩欧美免费精品| 看黄色毛片网站| 黄片小视频在线播放| 中文字幕色久视频| 亚洲成人精品中文字幕电影| 精品国产一区二区久久| 亚洲成人久久性| 中亚洲国语对白在线视频| 好看av亚洲va欧美ⅴa在| 99在线人妻在线中文字幕| svipshipincom国产片| 免费看a级黄色片| 在线观看日韩欧美| 国产一区二区三区视频了| 免费人成视频x8x8入口观看| 亚洲va日本ⅴa欧美va伊人久久| or卡值多少钱| 亚洲午夜理论影院| 嫩草影院精品99| 看免费av毛片| 久久精品影院6| 久久人妻av系列| 男女下面进入的视频免费午夜 | 最近最新中文字幕大全电影3 | 人人妻人人爽人人添夜夜欢视频| 最近最新中文字幕大全电影3 | 成人免费观看视频高清| 国产1区2区3区精品| 亚洲久久久国产精品| 性少妇av在线| 亚洲在线自拍视频| 日韩有码中文字幕| 操美女的视频在线观看| 一级黄色大片毛片| 久久中文字幕人妻熟女| 伦理电影免费视频| 亚洲三区欧美一区| 变态另类丝袜制服| 免费在线观看黄色视频的| 可以在线观看毛片的网站| 久久精品aⅴ一区二区三区四区| 97超级碰碰碰精品色视频在线观看| ponron亚洲| 成人手机av| 国产精品免费视频内射| 亚洲精品久久国产高清桃花| 久久香蕉精品热| 欧美在线一区亚洲| 激情在线观看视频在线高清| 99国产精品一区二区蜜桃av| 久久久久国产精品人妻aⅴ院| 免费观看人在逋| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| 午夜免费鲁丝| 欧美激情极品国产一区二区三区| 高潮久久久久久久久久久不卡| 在线免费观看的www视频| 国产高清videossex| 男人舔女人下体高潮全视频| 午夜成年电影在线免费观看| 国产高清激情床上av| 欧美另类亚洲清纯唯美| 午夜a级毛片| 成熟少妇高潮喷水视频| 一区二区日韩欧美中文字幕| 久久伊人香网站| 亚洲精品久久成人aⅴ小说| 老司机福利观看| 国产一区在线观看成人免费| 国产精品久久电影中文字幕| 国产精品久久久久久精品电影 | 亚洲色图综合在线观看| 美女高潮到喷水免费观看| 国产精品电影一区二区三区| 国产欧美日韩精品亚洲av| 亚洲欧洲精品一区二区精品久久久| 男人舔女人下体高潮全视频| 亚洲第一欧美日韩一区二区三区| 亚洲电影在线观看av| 老鸭窝网址在线观看| 亚洲一区二区三区不卡视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av成人一区二区三| 免费看美女性在线毛片视频| 久久精品亚洲熟妇少妇任你| 大型av网站在线播放| 欧美国产精品va在线观看不卡| 91av网站免费观看| 激情在线观看视频在线高清| 视频在线观看一区二区三区| 中文亚洲av片在线观看爽| 午夜免费激情av| 精品国产亚洲在线| 在线观看免费午夜福利视频| 精品人妻1区二区| 十分钟在线观看高清视频www| 精品人妻1区二区| 久久久久久亚洲精品国产蜜桃av| 亚洲国产精品sss在线观看| 十分钟在线观看高清视频www| 日韩成人在线观看一区二区三区| 天天躁夜夜躁狠狠躁躁| 日韩一卡2卡3卡4卡2021年| 日韩精品青青久久久久久| 无人区码免费观看不卡| 久久久国产欧美日韩av| 亚洲avbb在线观看| 午夜福利一区二区在线看| 国产人伦9x9x在线观看| 十八禁网站免费在线| 午夜免费鲁丝| 国产成人精品在线电影| 日韩国内少妇激情av| 高清黄色对白视频在线免费看| 曰老女人黄片| 十八禁人妻一区二区| 动漫黄色视频在线观看| bbb黄色大片| 午夜免费成人在线视频| 91麻豆av在线| 丰满的人妻完整版| 日韩大尺度精品在线看网址 | 99国产精品一区二区三区| 高清黄色对白视频在线免费看| 国语自产精品视频在线第100页| 久久久久久久久久久久大奶| 动漫黄色视频在线观看| 无人区码免费观看不卡| 亚洲国产精品999在线| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费视频网站a站| 亚洲精品国产色婷婷电影| 亚洲男人天堂网一区| 两个人看的免费小视频| 久久香蕉国产精品| 国产精品亚洲美女久久久| 高清黄色对白视频在线免费看| 国产精品免费视频内射| 国产欧美日韩一区二区精品| 亚洲精品一区av在线观看| 性欧美人与动物交配| 男女下面进入的视频免费午夜 | 亚洲国产欧美网| 久久久精品欧美日韩精品| 又黄又粗又硬又大视频| 一级,二级,三级黄色视频| 在线观看www视频免费| 欧美日韩福利视频一区二区| 久久人人精品亚洲av| 久久草成人影院| 日韩欧美一区视频在线观看| 午夜免费鲁丝| 国产精品av久久久久免费| 久久久久国产精品人妻aⅴ院| 久久精品91无色码中文字幕| 悠悠久久av| 性少妇av在线| 亚洲av成人不卡在线观看播放网| 成人手机av| 亚洲熟女毛片儿| 高清黄色对白视频在线免费看| 亚洲午夜理论影院| 十八禁网站免费在线| 伦理电影免费视频| 美女免费视频网站| 亚洲欧美激情综合另类| 99久久久亚洲精品蜜臀av| 欧美黑人欧美精品刺激| 两人在一起打扑克的视频| 精品久久久久久,| 欧美中文综合在线视频| 88av欧美| 91在线观看av| 亚洲五月天丁香| 亚洲电影在线观看av| bbb黄色大片| 亚洲第一av免费看| 久久国产精品人妻蜜桃| 精品久久久久久久毛片微露脸| 色婷婷久久久亚洲欧美| 十八禁人妻一区二区| 日韩免费av在线播放| 国产成人欧美在线观看| 亚洲片人在线观看| 免费观看精品视频网站| 制服丝袜大香蕉在线| 黄色成人免费大全| 操出白浆在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 国产av一区在线观看免费| 长腿黑丝高跟| 欧美 亚洲 国产 日韩一| 麻豆av在线久日| 欧美一区二区精品小视频在线| 操美女的视频在线观看| 国产精品久久久人人做人人爽| 午夜免费鲁丝| 黄频高清免费视频| 动漫黄色视频在线观看| e午夜精品久久久久久久| 两个人免费观看高清视频| 欧美久久黑人一区二区| 极品教师在线免费播放| 女同久久另类99精品国产91| 国产成人av教育| 日韩精品免费视频一区二区三区| 一级毛片女人18水好多| 搡老妇女老女人老熟妇| 精品一区二区三区视频在线观看免费| 夜夜夜夜夜久久久久| 日韩欧美一区二区三区在线观看| 真人做人爱边吃奶动态| av在线天堂中文字幕| 在线观看66精品国产| 日本一区二区免费在线视频| 在线观看一区二区三区| 久久青草综合色| 男女之事视频高清在线观看| 在线观看www视频免费| av片东京热男人的天堂| 亚洲激情在线av| 亚洲国产欧美网| 国产精品 国内视频| 12—13女人毛片做爰片一| 老司机午夜福利在线观看视频| 精品乱码久久久久久99久播| 成人精品一区二区免费| 男女下面插进去视频免费观看| 亚洲 欧美 日韩 在线 免费| 国产区一区二久久| 午夜免费观看网址| 在线观看舔阴道视频| 中亚洲国语对白在线视频| 男人舔女人的私密视频| 成熟少妇高潮喷水视频| 国产黄a三级三级三级人| 欧美乱码精品一区二区三区| 丝袜美腿诱惑在线| 中文字幕精品免费在线观看视频| 51午夜福利影视在线观看| 亚洲一区二区三区不卡视频| 制服丝袜大香蕉在线| 搡老熟女国产l中国老女人| 啦啦啦 在线观看视频| 精品国内亚洲2022精品成人| 国产三级黄色录像| 女人爽到高潮嗷嗷叫在线视频| 搡老妇女老女人老熟妇| 美女国产高潮福利片在线看| 亚洲五月天丁香| 久久精品91蜜桃| 日本一区二区免费在线视频| 最好的美女福利视频网| 亚洲人成电影免费在线| 色哟哟哟哟哟哟| 99香蕉大伊视频| 99国产极品粉嫩在线观看| 亚洲人成伊人成综合网2020| 国产精品影院久久| 夜夜夜夜夜久久久久| 亚洲国产精品999在线| 成熟少妇高潮喷水视频| 国产精品香港三级国产av潘金莲| 国产激情久久老熟女| 国产精品久久久人人做人人爽| 亚洲第一电影网av| 岛国视频午夜一区免费看| 黄色a级毛片大全视频| 国产午夜精品久久久久久| 日日爽夜夜爽网站| 亚洲中文av在线| 国产精品乱码一区二三区的特点 | 国产亚洲精品一区二区www| 制服丝袜大香蕉在线| 搞女人的毛片| 日韩av在线大香蕉| 国产精品1区2区在线观看.| 欧美成人一区二区免费高清观看 | 久久久久九九精品影院| av天堂在线播放| 最新美女视频免费是黄的| 日本免费a在线| 国产99白浆流出| 成人亚洲精品av一区二区| 国产黄a三级三级三级人| 国语自产精品视频在线第100页| 精品无人区乱码1区二区| 国产91精品成人一区二区三区| 亚洲专区字幕在线| 一级作爱视频免费观看| 少妇被粗大的猛进出69影院| 亚洲av成人一区二区三| 久久性视频一级片| 久久狼人影院| 真人做人爱边吃奶动态| 国产成人影院久久av| 精品久久久久久,| 欧美最黄视频在线播放免费| 最近最新免费中文字幕在线| 亚洲精品一区av在线观看| 黄色a级毛片大全视频| 免费人成视频x8x8入口观看| 久久久久亚洲av毛片大全| 黑人欧美特级aaaaaa片| 69av精品久久久久久| 亚洲伊人色综图| 国产人伦9x9x在线观看| 搡老妇女老女人老熟妇| 91字幕亚洲| 九色国产91popny在线| 亚洲一码二码三码区别大吗| 波多野结衣高清无吗| 人人澡人人妻人| 色播亚洲综合网| 国产av在哪里看| 久久青草综合色| 欧美成人一区二区免费高清观看 | 夜夜躁狠狠躁天天躁| 夜夜爽天天搞| 99国产极品粉嫩在线观看| 亚洲欧美日韩高清在线视频| 久久久久久久久久久久大奶| 91成人精品电影| 国内精品久久久久精免费| 亚洲精华国产精华精| 亚洲国产毛片av蜜桃av| 一级毛片高清免费大全| 国产色视频综合| 男女做爰动态图高潮gif福利片 | 国产一区二区在线av高清观看|