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

    檢波器性能參數(shù)現(xiàn)場測定方法

    2021-06-01 12:20:24魏繼東曹國濱
    石油地球物理勘探 2021年3期
    關(guān)鍵詞:檢波器振幅靈敏度

    魏繼東 曹國濱 劉 斌

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

    0 引言

    檢波器是地震勘探中接收地表振動(dòng)的硬件裝置。現(xiàn)今仍廣泛采用的速度型電感動(dòng)圈式檢波器(Electrical Seismometer)于20世紀(jì)30年代被發(fā)明并隨即得到推廣應(yīng)用,此前的接收設(shè)備主要是位移型機(jī)械式檢波器(Mechanical Seismometer)[1-2]。近十幾年來,隨著機(jī)電技術(shù)的發(fā)展,地震數(shù)據(jù)采集施工中出現(xiàn)了以MEMS為代表的加速度型檢波器(Accelerometer)。不論何種類型檢波器,“忠實(shí)(保真)地記錄大地振動(dòng)”是其最主要功能。

    為了表征檢波器接收信號的保真能力(任何檢波器輸出的電信號都是地表機(jī)械振動(dòng)的“模擬與近似”),試用了各種指標(biāo)參數(shù)定量表達(dá)其多方面性能,如靈敏度、自然頻率、阻尼系數(shù)、允差,等等?;诘卣鹂碧降木炔粩嗵岣?,要求檢波器具有盡量高的機(jī)電(轉(zhuǎn)換)參數(shù)保證數(shù)據(jù)畸變盡量小。在地震檢波技術(shù)的發(fā)展過程中,人們發(fā)展了很多室內(nèi)和現(xiàn)場的測試方法標(biāo)定檢波器的性能參數(shù)。其中,實(shí)驗(yàn)室方法受限于設(shè)備、場地等因素,只能檢測少量檢波器,無法及時(shí)、準(zhǔn)確地檢測實(shí)際地震數(shù)據(jù)采集中大量檢波器的多種性能參數(shù),部分重要指標(biāo)甚至無法測量。另一方面,由于沒有厘清檢波器性能指標(biāo)與地震數(shù)據(jù)表現(xiàn)之間正確的因果關(guān)系,當(dāng)前通過比較不同檢波器的頻譜、能量、信噪比、同相軸連續(xù)性等指標(biāo),進(jìn)而評價(jià)檢波器性能參數(shù)的方法是非常粗放甚至錯(cuò)誤的,無法觀測到檢波器性能指標(biāo)的細(xì)部特征,容易產(chǎn)生誤導(dǎo)[3-6]。因此,對于當(dāng)前檢波器性能參數(shù)標(biāo)定方法而言,無論是少量的實(shí)驗(yàn)室檢測還是大規(guī)模的現(xiàn)場試驗(yàn)均存在一定的局限性。施工單位進(jìn)行了大量的、重復(fù)性試驗(yàn),但方法不一,結(jié)論多樣化甚至相互矛盾,未取得穩(wěn)定的、統(tǒng)一的認(rèn)識(shí)與做法,在一定程度上影響了現(xiàn)場采集數(shù)據(jù)質(zhì)量的提高。

    基于以上現(xiàn)狀,在厘清檢波器工作原理與對應(yīng)數(shù)據(jù)特征的基礎(chǔ)上,本文提出一種基于現(xiàn)場試驗(yàn)的估測大量(數(shù)萬道)在用檢波器指標(biāo)參數(shù)的便捷方法,可在占用較小場地、僅用少量單炮數(shù)據(jù)的情況下,在數(shù)分鐘內(nèi)大致估測檢波器的靈敏度、最低可靠頻率、允差等較精細(xì)的指標(biāo),為確保數(shù)據(jù)接收質(zhì)量提供硬件保證。此處之所以說“大致”而非“準(zhǔn)確”,是因?yàn)槭軓?fù)雜現(xiàn)場環(huán)境的影響,估測結(jié)果很難達(dá)到實(shí)驗(yàn)室條件下的測量精度。但因?qū)嶋H數(shù)據(jù)是在現(xiàn)場環(huán)境下采集的,所以該方法可提供一個(gè)從最終“數(shù)據(jù)意義”上評判檢波器是否“可用”的框架性標(biāo)準(zhǔn)。同時(shí),也可用于檢波器研制過程中在缺乏炸藥震源時(shí)指標(biāo)參數(shù)的快速標(biāo)定。

    1 方法原理

    “地震儀+檢波器+地表介質(zhì)”三者的組合可被稱為“地震檢波系統(tǒng)[4]”。從“信號與系統(tǒng)”的角度來講,如果多個(gè)系統(tǒng)的輸入是一致的或高度一致的,其輸出主要反映了系統(tǒng)特性/傳輸函數(shù)的差異。因此,如果同時(shí)給多個(gè)檢波系統(tǒng)輸入一個(gè)高度相似的機(jī)械信號,然后比較其輸出,可觀測到檢波器機(jī)電參數(shù)以及檢波器—大地耦合響應(yīng)的差異。若其中某個(gè)系統(tǒng)的參數(shù)是已知的,可根據(jù)系統(tǒng)其他輸出與已知系統(tǒng)輸出的差異估測系統(tǒng)的其他相應(yīng)參數(shù)(圖1)。

    根據(jù)以上思路,本文提出一種基于現(xiàn)場試驗(yàn)數(shù)據(jù)估測檢波器性能參數(shù)的方法,稱為“微型盒子波(Micro box wave test)”,其原理如下文所述。

    圖1 輸入相同的情況下輸出差異反映系統(tǒng)差異

    1.1 地表臨近質(zhì)點(diǎn)振動(dòng)特征具有高度一致性(輸入分析)

    地震波在地下介質(zhì)中以真速度v沿射線傳播,但在除射線方向以外的任意方向(如地表)觀測時(shí),地震波以視速度v'傳播。將視速度定義為[7]

    (1)

    式中θ為入射角。

    式(1)說明: 當(dāng)入射角θ=0°時(shí),視速度等于真速度; 當(dāng)θ≠0°時(shí),視速度總大于真速度; 當(dāng)入射角θ趨于90°極限時(shí),視速度趨于無窮大,此時(shí)觀測面上各點(diǎn)的擾動(dòng)幾乎同時(shí)到達(dá)。

    圖2顯示低頻檢波器接收到各種常見地震波的視波長譜[8]。從中可見,除次生低速干擾速度稍低外,震源激發(fā)的有效波和干擾波的視速度幾乎都大于100m/s,直到無窮大。因此,若將檢波器布設(shè)在地表1m×1m范圍內(nèi),則地表振動(dòng)(不考慮檢波器因素)的差異主要來源于三個(gè)方面: ①傳播距離; ②大地吸收; ③不具有相干性的環(huán)境噪聲。

    圖2 低頻檢波器接收時(shí)各種地震波的視波長譜[8]

    1.1.1 傳播距離

    陸上震源被激發(fā)后,對于沿地表傳播的地震波而言,1m地表距離產(chǎn)生約3ms時(shí)差(假設(shè)低速帶速度為360m/s)。此時(shí)若選擇頻譜分析時(shí)窗長度為500ms,則相距1m的兩個(gè)質(zhì)點(diǎn)之間的振動(dòng)差異主要在頭、尾兩端(圖3紅色部分),其他部分是相同的(不考慮大地吸收),即最大只有3/500=0.6%是不同的(在頭、尾部振動(dòng)完全不相關(guān)時(shí)); 若分析時(shí)長為3000ms,則最大只有3/3000=0.1%是不同的。

    對于不沿地表傳播的地震波來說,因視速度非常高,故在相同分析時(shí)段內(nèi),其差異非常小。因此,對于陸上勘探激發(fā)的各種地震波及環(huán)境噪聲,1m×1m范圍內(nèi)的地表振動(dòng)具有“高度一致性”。

    圖3 相距1m兩質(zhì)點(diǎn)振動(dòng)記錄在相同時(shí)段具有高度相似性 未考慮大地吸收

    1.1.2 大地吸收

    大地對地震信號具有強(qiáng)烈的吸收衰減作用,地層縱波速度與品質(zhì)因數(shù)之間滿足經(jīng)驗(yàn)公式

    Q≈14×VP2.2

    (2)

    β≈1.949×VP-2.2

    (3)

    式中:Q為品質(zhì)因數(shù)(無單位);β為衰減系數(shù),定義為每傳播一個(gè)波長的振幅衰減量;VP為縱波速度。

    式(2)、式(3)表征地震波衰減的大體規(guī)律,有助于理解大套地層吸收的總趨勢[9]。在此基礎(chǔ)上,根據(jù)工區(qū)地層的VP層速度剖面,可算出高頻信息被衰減的程度[10]。如對于360m/s的低速帶來說,傳播1m距離的累計(jì)衰減率約為0.05dB/Hz[10],據(jù)此可推算100Hz雷克子波從炮檢距“X”m(圖4藍(lán)線)傳播到“X+1”m(圖4綠線)時(shí)振幅衰減的幅度。可見該子波200Hz處大約衰減了-12dB,變?yōu)樵确鹊?/4,這一衰減量有可能導(dǎo)致地表1m×1m范圍內(nèi)質(zhì)點(diǎn)振動(dòng)的較大差異。

    通過比較多個(gè)相距1m的同類檢波器接收的實(shí)際地表振動(dòng)的振幅譜(圖5為其中一個(gè)),發(fā)現(xiàn)距離炮點(diǎn)較近檢波點(diǎn)的高頻成分并沒有比炮檢距增大1m后的檢波點(diǎn)有明顯增強(qiáng)。因?yàn)榈卣鹩涗浭怯裳厮椒较騻鞑サ牟?如面波、初至波、次生干擾波等)以及主要沿垂直方向傳播的地震波(如反射波)構(gòu)成,出現(xiàn)圖5現(xiàn)象的原因可能有以下幾個(gè)方面。

    (1)對于水平傳播的各種地震波來說,受大地吸收影響最大的是沿地表低速帶傳播的初至波(吸收強(qiáng)烈且有高頻成分),對于低頻(吸收弱)面波、次生干擾波等,1m距離產(chǎn)生的吸收差異非常?。?/p>

    圖4 地震波沿地表傳播由大地吸收導(dǎo)致的理論振幅差異

    圖5 相距1m兩檢波器接收地震波的振幅譜

    (2)對于反射波而言,因?yàn)榈徒邓賻У挠绊?,大致垂直地入射到地表,地表相?m的兩個(gè)質(zhì)點(diǎn)對應(yīng)的大地吸收差異非常小[10];

    (3)圖5中炮檢距更小的高頻端振幅沒有明顯偏高的另外一個(gè)原因是檢波器—大地耦合響應(yīng)。大地吸收衰減了高頻成分,但是耦合響應(yīng)會(huì)放大高頻成分,后者的作用超過了前者,使得前者在圖5中未呈現(xiàn)明顯的規(guī)律性。

    (4)或者大套地層的吸收規(guī)律不適于遠(yuǎn)小于波長尺度的吸收量的計(jì)算,尚不確定。

    1.1.3 環(huán)境噪聲

    環(huán)境噪聲具有各態(tài)歷經(jīng)性質(zhì)[7],所以1m×1m范圍內(nèi)差異非常小。

    綜上所述,無論從傳播距離、大地吸收還是環(huán)境噪聲的角度來看,在未植入檢波器的情況下,1m×1m地表范圍內(nèi)由陸上震源激發(fā)的質(zhì)點(diǎn)振動(dòng)在環(huán)境噪聲背景下具有高度一致性,頻率越低一致性越強(qiáng)——因?yàn)橛纱蟮匚諏?dǎo)致的差異越小。

    1.2 地震檢波系統(tǒng)對地震數(shù)據(jù)的影響(系統(tǒng)分析)

    將地表機(jī)械振動(dòng)轉(zhuǎn)換為可用于后續(xù)數(shù)學(xué)計(jì)算的離散數(shù)字過程中,地震檢波系統(tǒng)主要完成了機(jī)電轉(zhuǎn)換與模數(shù)轉(zhuǎn)換兩個(gè)過程,對地震數(shù)據(jù)的主要響應(yīng)可總結(jié)為圖6[4]。地震檢波系統(tǒng)核心要素是檢波器。

    隨著相關(guān)技術(shù)的發(fā)展,陸上石油勘探中出現(xiàn)過各種基于不同工作原理的檢波器,如電感動(dòng)圈式模擬檢波器(SM-24、20DX)、MEMS數(shù)字檢波器(DSUs、VectorSeis)等。表1為SM-24檢波器的參數(shù)表。關(guān)于檢波器指標(biāo)參數(shù)的具體解釋參見文獻(xiàn)[11-12],此處不做詳述。

    文獻(xiàn)[11]表明:靈敏度決定檢波器對機(jī)械信號的電放大能力;自然頻率、阻尼以及電阻決定檢波器的低頻濾波響應(yīng); 允差代表多個(gè)檢波器之間的一致性; 與地震儀或者有源檢波器相關(guān)的本底噪聲對地震數(shù)據(jù)的記錄質(zhì)量有重要影響,主要影響高、低頻兩端的機(jī)電比; 檢波器的物理特征包括直徑、重量、高度等與檢波器—大地耦合響應(yīng)有關(guān)。

    1.2.1 靈敏度

    檢波器靈敏度指的是機(jī)械振動(dòng)轉(zhuǎn)換為電信號后的電壓與其位移、速度或者加速度的比值[14]。不同類型檢波器的靈敏度單位不同,速度型一般為V·m-1·s,加速度型一般為V·m-1·s2或V/g。將SEGY數(shù)據(jù)除以靈敏度或SEGD數(shù)據(jù)乘以標(biāo)定算子(Depreamplifer Scaling Factor)后再除以靈敏度,即可將數(shù)據(jù)統(tǒng)一為速度或加速度量綱。

    圖6 地震檢波系統(tǒng)對機(jī)械振動(dòng)的記錄與改造作用[4]

    1.2.2 低頻濾波響應(yīng)

    不同檢波器具有不同的低頻濾波響應(yīng)[4,6,11]。動(dòng)圈式檢波器低頻濾波響應(yīng)主要受主頻和阻尼系數(shù)兩個(gè)因素影響,主頻以下存在-12dB/OCT的衰減;MEMS檢波器(如DSU3)理論上不存在低頻衰減(實(shí)際上會(huì)由于本底噪聲以及加速度接收的原因而導(dǎo)致較大的極低頻失真)。圖7為兩種典型檢波器(20DX、DSU3)的幅頻響應(yīng)曲線。

    1.2.3 允差

    允差指對某給定檢波器參數(shù)所允許的誤差極限值[13]。允差主要影響數(shù)據(jù)的一致性。由表1可見,主頻、阻尼、電阻及靈敏度等都存在允差概念。

    圖7 兩種典型檢波器的幅頻曲線[5]

    1.2.4 本底噪聲

    檢波系統(tǒng)的本底噪聲(Floor noise)指的是在沒有外界機(jī)械輸入的情況下仍然存在的干擾信號[13]。本底噪聲主要有兩個(gè)來源[14-15]: 一是空氣分子不停的布朗運(yùn)動(dòng)撞擊質(zhì)量塊和彈簧振子引起的布朗熱噪聲; 二是有源檢波器或者地震儀中的電子噪聲。對于無源檢波器而言,電子噪聲主要來源于地震儀前置放大器等元器件; 對于有源檢波器來說,除地震儀噪聲外其自身也會(huì)產(chǎn)生電子噪聲。本文將布朗熱噪聲以及地震儀(包含有源檢波器)產(chǎn)生的電子噪聲均稱為本底噪聲,不再做更細(xì)致區(qū)分。本底噪聲主要影響地震數(shù)據(jù)的絕對動(dòng)態(tài)范圍[6]。同時(shí),由于本底噪聲中的電子噪聲具有1/f的特征(頻率越低強(qiáng)度越大),其強(qiáng)弱對極低頻以及高頻弱信號的記錄有重要影響。表2為DSU1檢波器的部分指標(biāo),其中Noise主要表征了本底噪聲的水平。

    1.2.5 檢波器—大地耦合響應(yīng)

    檢波器的不同材質(zhì)、外形、重量等物理特征會(huì)導(dǎo)致其與大地之間不同的耦合響應(yīng)[3-5]。圖8a是兩個(gè)檢波器相距20cm時(shí)有意差異性埋置時(shí)接收到的振幅譜。在其他因素一致的情況下,二者在高頻端的差異是由于耦合條件不同導(dǎo)致的。

    圖8 兩個(gè)同類檢波器因耦合條件導(dǎo)致的振幅差異(a)以及兩條示意性理論耦合響應(yīng)振幅曲線(b)

    表1 SM24(SENSOR)部分指標(biāo)

    表2 DSU1(SERCEL)部分指標(biāo)

    從振動(dòng)力學(xué)的角度來看,檢波器—大地之間的耦合系統(tǒng)是一個(gè)單自由度機(jī)械振動(dòng)系統(tǒng)。因?yàn)槠胀z波器與陸地地表的耦合諧振頻率多為幾百赫茲,所以低于耦合諧振頻率的頻率成分具有頻率越高、振幅越強(qiáng)的特征(圖8b)。根據(jù)本研究測量,耦合響應(yīng)對地震數(shù)據(jù)的影響多在40~300Hz頻段(40Hz以下耦合響應(yīng)弱,300Hz以上地表振動(dòng)弱)。

    1.3 根據(jù)輸出數(shù)據(jù)反演系統(tǒng)特性(輸出分析)

    根據(jù)以上輸入、系統(tǒng)因素分析,為了估測檢波器的性能指標(biāo),筆者參照文獻(xiàn)[16]中用于調(diào)查干擾波的Box Wave Test方法,設(shè)計(jì)了前述的用于估測檢波器性能指標(biāo)的“微型盒子波法”(圖9)。該方法設(shè)定的檢波器間距較文獻(xiàn)[16]小了約一個(gè)數(shù)量級(1.0m→0.2m),目的是估測檢波器的指標(biāo)參數(shù)。

    圖10為4×4個(gè)檢波器(20DX)按照圖9所示排列接收數(shù)據(jù)的振幅譜。根據(jù)文獻(xiàn)[5],地震數(shù)據(jù)會(huì)受到以下因素的影響:機(jī)械振動(dòng)、機(jī)電轉(zhuǎn)換過程中的濾波響應(yīng)、本底噪聲、耦合響應(yīng)以及電子元器件中數(shù)學(xué)運(yùn)算(如高截濾波)等。但對于不同的頻段而言,其主導(dǎo)因素是不同的,大致可劃分為六個(gè)頻段。就圖10來說,由地表機(jī)械振動(dòng)轉(zhuǎn)換的信號主要影響頻段為2、3、4(綠、青、黃)。這三個(gè)頻段震源激發(fā)的信號較強(qiáng),振幅譜的主要形態(tài)由地表機(jī)械信號決定。同時(shí),每個(gè)檢波器具有不同的低頻濾波響應(yīng)、耦合響應(yīng)(主要影響數(shù)據(jù)的高頻端[17-19]),因此頻段2的主要影響因素是機(jī)械信號+低頻濾波響應(yīng),頻段3的主要影響因素是機(jī)械信號,頻段4則為機(jī)械信號+耦合響應(yīng)。由圖10可見,因震源激發(fā)的地震波視波長很大(?100m),故分布在1m×1m范圍內(nèi)的所有檢波器在頻段2、3的一致性非常高。但超過40Hz以后,由于耦合響應(yīng)的作用,16個(gè)檢波器之間的一致性降低了。

    頻段5(300~400Hz)由于輸入的機(jī)械信號(震源激發(fā)+環(huán)境因素)大大減弱,數(shù)據(jù)中本底噪聲的占比大大增加了,呈現(xiàn)出類似“白噪”的態(tài)勢。

    頻段6的影響因素與頻段5基本相同,但增加了作用明顯的高截濾波,所以其幅度在400Hz(0.8×Nyquist 頻率)以上迅速衰減。

    對于頻段1來說,因?yàn)檎ㄋ幖ぐl(fā)的低頻分量迅速衰減,由機(jī)械振動(dòng)轉(zhuǎn)換而來的電信號非常微弱(但一致性很強(qiáng)); 同時(shí),本底噪聲具有1/f規(guī)律[5,14-15],則該頻段的主導(dǎo)因素變?yōu)楸镜自肼暋S忠蛎總€(gè)檢波器的本底噪聲各不相同,故曲線組的一致性顯著降低。圖11為8個(gè)加速度檢波器DSU3(垂直分量)布設(shè)在1m×1m范圍內(nèi)錄制數(shù)據(jù)的振幅譜。因?yàn)榧铀俣葯z波器具有弱化低頻的作用,本底噪聲的表現(xiàn)更加明顯(紅色部分),低頻端保真度顯著降低。

    圖9 微型盒子波排列示意圖

    這樣,對于接收地表機(jī)械振動(dòng)的地震檢波系統(tǒng)[5]來說,影響地震數(shù)據(jù)(輸出)頻域特征的主要因素包括地表機(jī)械振動(dòng)(輸入)以及低頻濾波響應(yīng)、本底噪聲、檢波器—大地耦合響應(yīng)、允差、靈敏度等。以上因素對地震數(shù)據(jù)的影響可以歸結(jié)為圖12(藍(lán)、紅、綠三種顏色分別代表了不同的檢波器)以及表3(表3中頻帶范圍的劃分?jǐn)?shù)值限于圖10具體情況,并不絕對)。

    另外,若將同類或不同類檢波器埋置到1m×1m狹小面積內(nèi),震源激發(fā)的機(jī)械振動(dòng)高度一致(頻率越低一致性越強(qiáng)),輸出地震數(shù)據(jù)的差異主要反映地震檢波系統(tǒng)特性的差異。同時(shí),由圖12可見20~40Hz是較“特殊”頻段,具有三個(gè)特征: ①大地吸收衰減少、輸入信號一致性強(qiáng); ②地震波能量最強(qiáng)、機(jī)電比最高的頻段,受本底噪聲影響最小(地震子波是鐘型譜); ③既不受低頻濾波響應(yīng)影響,又幾乎不受高頻耦合響應(yīng)的影響。

    圖10 4×4個(gè)檢波器接收同一炮的振幅譜

    圖11 本底噪聲使DSU3檢波器極低頻保真度顯著降低

    圖12 不同性能參數(shù)指標(biāo)對地震數(shù)據(jù)的影響

    在這個(gè)頻段內(nèi),振幅主要決定于靈敏度、振幅響應(yīng)以及二者的變化范圍(允差);如果通過多次疊加消除允差影響、通過頻帶平均消除頻帶差異,則振幅只決定于靈敏度這一單一因素,20~40Hz頻帶內(nèi)振幅均值與檢波器的靈敏度成正比。若不進(jìn)行多次疊加,輸出數(shù)據(jù)振幅的變動(dòng)范圍則反映了靈敏度以及振幅響應(yīng)的變化范圍(無論靈敏度還是振幅響應(yīng)在不同頻率上的分布都不會(huì)保持絕對“平直”),大致相當(dāng)于允差的范圍。

    因此,若將檢波器按照圖9所示設(shè)計(jì)形態(tài)布設(shè),其輸出數(shù)據(jù)的差異主要反映地震檢波系統(tǒng)的差異,這樣就為現(xiàn)場快捷估測檢波器性能參數(shù)提供了可能性。

    表3 地震信號以及地震檢波系統(tǒng)各因素作用的主要頻帶范圍

    2 試驗(yàn)驗(yàn)證

    基于上述原理,對靈敏度、最低可靠頻率以及允差、綜合一致性等指標(biāo)的估測進(jìn)行了試驗(yàn)驗(yàn)證。試驗(yàn)時(shí)檢波器按圖9所示方式擺放,炮檢距大約幾百米,工區(qū)常用激發(fā)因素(如激發(fā)井深為10m、激發(fā)藥量為4kg)、儀器因素(如采樣間隔為1ms、記錄長度為6s、前放增益為12dB)即可。

    2.1 靈敏度

    靈敏度反映了檢波器對機(jī)械信號進(jìn)行“電放大”的能力。根據(jù)前述,可以根據(jù)多個(gè)檢波器累加后20~40Hz頻帶內(nèi)振幅均值標(biāo)定檢波器的靈敏度。

    為了標(biāo)定新型檢波器的靈敏度,可選取已有成熟檢波器作為標(biāo)準(zhǔn)。如測定加速度檢波器的靈敏度時(shí),可采用DSU系列檢波器作為標(biāo)準(zhǔn);測定速度檢波器時(shí),可采用20DX或者其他成熟的速度型檢波器作為標(biāo)準(zhǔn)。因?yàn)橥N類型檢波器(速度/速度,加速度/加速度)之間的靈敏度是簡單的倍數(shù)關(guān)系,不同類型檢波器(速度/加速度)的靈敏度則依賴于頻率,計(jì)算不夠直觀。

    圖13是10種檢波器的測試結(jié)果,其中超級檢波器的靈敏度已知(19.7V·m-1·s,紅色),DSU3的靈敏度已知(4.43V/g,綠色)。通過計(jì)算得知,5Hz動(dòng)圈式檢波器靈敏度約為72V·m-1·s(粉色), 第1~6、10種加速度檢波器的靈敏度約為10V/g(藍(lán)色),均與室內(nèi)試驗(yàn)結(jié)果相符,從而證明了這種方法的正確性。當(dāng)然,非常準(zhǔn)確的測量需在實(shí)驗(yàn)室進(jìn)行,但這種估測結(jié)果對于實(shí)際數(shù)據(jù)采集具有現(xiàn)實(shí)參考意義。

    根據(jù)作者經(jīng)驗(yàn),陸上石油勘探單點(diǎn)采集時(shí)速度型檢波器靈敏度在100~120V·m-1·s、加速度檢波器靈敏度在4~5V/g范圍較為合理,既可以充分利用地震儀的動(dòng)態(tài)范圍,又不會(huì)產(chǎn)生超調(diào)現(xiàn)象。

    作者對此方法進(jìn)行了多炮檢距驗(yàn)證,結(jié)果都是一致的,進(jìn)一步說明了這種方法的穩(wěn)定性與正確性。

    圖13 由微型盒子波估算的檢波器靈敏度

    2.2 最低可靠頻率

    決定檢波器低頻接收能力的主要因素包括檢波器的低頻濾波響應(yīng)(圖7)、本底噪聲(表2)以及靈敏度(主要影響弱信號),可用“最低可靠頻率”這一單一指標(biāo)進(jìn)行直觀衡量。作者將振幅曲線組一致性顯著變差的轉(zhuǎn)折點(diǎn)稱為該型檢波器的“最低可靠頻率LFF-Lowest Faithful Frequency”(圖11紅點(diǎn)),代表在此情況下最低的保真頻率,由機(jī)械信號的強(qiáng)度、本底噪聲、靈敏度以及低頻濾波響應(yīng)共同決定(前三者相同的情況下,低頻濾波響應(yīng)越強(qiáng),低頻記錄能力越差、最低可靠頻率越高)。

    地震儀輸出的數(shù)字信號來源于電信號,而電信號有兩個(gè)來源:由機(jī)械振動(dòng)轉(zhuǎn)換而來的電信號和本底噪聲。本底噪聲屬地震儀和檢波器固有的噪聲,其特征與機(jī)械輸入信號無關(guān)。如果輸入的振動(dòng)信號足夠強(qiáng)、超過本底噪聲,就可以比較保真地記錄振動(dòng)信號;如果輸入振動(dòng)信號與本底噪聲相當(dāng)甚至更弱,本底噪聲的影響就會(huì)顯現(xiàn)。圖14中藍(lán)色環(huán)境噪聲的極低頻逐漸減弱,而粉色本底噪聲逐漸增強(qiáng)(1/f規(guī)律);在本底噪聲強(qiáng)于環(huán)境噪聲的頻段(極低頻段),前者占主導(dǎo)作用,此時(shí)輸出的電信號不再代表地表振動(dòng),而是反映了本底噪聲的強(qiáng)度。部分加速度檢波器數(shù)據(jù)積分為速度后會(huì)產(chǎn)生強(qiáng)烈的極低頻,基于此原因,應(yīng)該將其切除。

    在圖11中,頻段2、3的主要構(gòu)成要素是振動(dòng)信號轉(zhuǎn)換而來的電信號,具有高度的“空間一致性”,所以8個(gè)檢波器高度一致,且越趨向低頻一致性越強(qiáng)(低頻成分視波長更大、吸收更小)。由此一致性可以推知靈敏度、允差以及低頻濾波響應(yīng)都是一致的。但是,頻段1(極低頻,大約2.5Hz以下)中各曲線一致性顯著變差,出現(xiàn)了向上突然抬升的態(tài)勢。因其他因素都是一致的,故極低頻出現(xiàn)差異的唯一原因是“本底噪聲”(表3)。

    本底噪聲具有1/f及隨機(jī)分布的特征[5,14-15],一致性變差的頻率點(diǎn)反映了地震信號由機(jī)械輸入信號為主向本底噪聲為主的轉(zhuǎn)變。

    圖15左圖為10種檢波器(除第9種為8個(gè)以外,其余均為16個(gè)檢波器,記錄時(shí)段0~6s)的低頻振幅譜比較,圖中低頻端振幅譜變粗的頻率點(diǎn)即為對應(yīng)檢波器類型的“最低有效頻率”。由圖15可見,第6種檢波器的低頻接收能力最差,第7種檢波器低頻接收能力最強(qiáng);第9種DSU檢波器的最低有效頻率大約為4Hz(與圖11不同)。另外,隨著輸入信號強(qiáng)度降低,允差、靈敏度等其他因素的作用更明顯,本底噪聲的影響減弱,最低有效頻率就難以標(biāo)定(圖15b,分析時(shí)段5~6s)。因此,最低可靠頻率對比應(yīng)在較強(qiáng)機(jī)械輸入的情況下進(jìn)行。當(dāng)然,經(jīng)多次疊加衰減本底噪聲后最低可靠頻率會(huì)降低。

    圖14 本底噪聲與環(huán)境噪聲的功率譜比較

    圖15 10種檢波器的低頻記錄能力比較(由上到下分別對應(yīng)10種檢波器類型) 左圖對應(yīng)時(shí)段為0~6s,右圖對應(yīng)時(shí)段為5~6s

    2.3 允差

    根據(jù)前述,盡管1m×1m范圍內(nèi)地表機(jī)械振動(dòng)具有高度一致性,但實(shí)際數(shù)據(jù)顯示部分檢波器類型不同檢波器之間的曲線一致性較好、曲線重合度較高(圖16B1、C1),部分一致性則較差(圖16A1、D1)。因?yàn)锳1、D1是加速度檢波器,具有突出視覺高頻的作用(差異大的耦合響應(yīng)作用在高頻,環(huán)境噪聲較震源信號在高頻端也相對更強(qiáng));同時(shí),二者均為有源檢波器,本底噪聲更強(qiáng);這些都是A1、D1一致性差的重要原因。但對于同種檢波器而言,耦合響應(yīng)及本底噪聲應(yīng)極為接近,故圖16中A1、D1一致性差的主要原因是允差。允差對于保持地震數(shù)據(jù)的一致性具有重要影響,但是現(xiàn)場難以測量這一指標(biāo)。在檢波器的性能參數(shù)表中,影響允差的主要參數(shù)包括主頻、阻尼、電阻以及靈敏度等(表1。有源檢波器的本底噪聲也存在允差的概念),多數(shù)在2.5%~5%范圍。

    對于分布在1m×1m范圍內(nèi)的多個(gè)檢波器來說,極低頻容易受到本底噪聲的影響,高頻容易受到耦合響應(yīng)的影響,低、中頻段機(jī)械信號的強(qiáng)度大、機(jī)電比高,作用的主要系統(tǒng)因素有兩個(gè):靈敏度和振幅響應(yīng)允差(暫稱為“合并允差”)。如果把每一組所有檢波器振幅譜的分布范圍作為考核標(biāo)準(zhǔn),就可大致算出該型檢波器合并允差的大小。具體的計(jì)算方法如下。

    (1)計(jì)算每一類型所有檢波器的振幅譜。

    (2)將同一類中所有檢波器各個(gè)頻率振幅譜的“(最大值-最小值)/平均值”作為該頻率的“合并允差”。因?yàn)椤白畲笾?最小值”的范圍代表了該頻率振幅的變動(dòng)范圍,平均值則將不同檢波器的變動(dòng)范圍進(jìn)行了歸一化,實(shí)現(xiàn)了可比性。

    圖17是實(shí)際試驗(yàn)結(jié)果,可見20~40Hz頻段曲線基本平直,反映了檢波器的“合并允差”。其中兩種動(dòng)圈式檢波器(B1—?jiǎng)尤?Hz、C1—?jiǎng)尤?0Hz)的合并允差接近10%,DSU3(D1)大約為18%,而LP-5Hz(A1)則達(dá)到約35%,大約是DSU檢波器的兩倍,顯然太高了,嚴(yán)重影響了數(shù)據(jù)一致性。低頻端的抬升是由于本底噪聲導(dǎo)致的,其轉(zhuǎn)折點(diǎn)對應(yīng)最低可靠頻率,強(qiáng)度大致反映了檢波器本地噪聲的相對強(qiáng)弱。由圖17可知,LP-5Hz、動(dòng)圈10Hz、DSU3、動(dòng)圈5Hz的低頻接收能力依次增強(qiáng)。超過40Hz以上的曲線抬升是由于耦合響應(yīng)不一致導(dǎo)致的。

    因?yàn)槎鄶?shù)反射信號較初至弱很多,所以更容易受到系統(tǒng)因素的影響。圖16下部為4~4.5s的波形圖,可見各個(gè)檢波器之間的一致性較初至(圖16上部)大大下降了,這種不一致使深層反射弱信號的識(shí)別更困難,此不一致性不僅僅來源于合并允差,也跟其他系統(tǒng)因素有關(guān)。為了簡潔、全面地反映不同地震檢波系統(tǒng)的機(jī)電特性,可用綜合一致性對此進(jìn)行表征。

    從傅里葉變換的角度來講,地震檢波系統(tǒng)的所有性能參數(shù)可從每一個(gè)樣點(diǎn)上得到體現(xiàn)。特別地,最小炮檢距的最大值是地震記錄中機(jī)電比最高、輸入一致性最強(qiáng)的樣點(diǎn),其差異最大限度地反映了系統(tǒng)的差異。因此,可用微型盒子波條件下每個(gè)檢波器的最大值(不超調(diào)時(shí))的變動(dòng)范圍量化表征某種類型檢波器的一致性。該指標(biāo)可在一定程度上表示某種類型檢波器的各項(xiàng)指標(biāo)對地震數(shù)據(jù)共同作用后的結(jié)果,與允差類似,但對其指標(biāo)體系的綜合表征作用更強(qiáng),可被稱為綜合一致性。圖18中的綠、黑、粉、藍(lán)曲線分別對應(yīng)圖16中A、B、C、D四種檢波器,其變動(dòng)范圍分別為±20%、±2.5%、±2.5%、±7.5%,綜合反映了檢波系統(tǒng)的機(jī)電參數(shù)浮動(dòng)范圍。如果此范圍過大,弱信號的識(shí)別就難以實(shí)現(xiàn)。

    圖17 20~40Hz的振幅譜相對分布范圍反映了合并允差的大小

    圖18 圖17四種檢波器的綜合一致性對比

    3 具體應(yīng)用中的若干問題

    在具體應(yīng)用中應(yīng)該注意以下幾個(gè)問題。

    3.1 觀測系統(tǒng)

    新型檢波器參數(shù)標(biāo)定與采集作業(yè)用檢波器一致性篩選(挑選不合格檢波器),這兩種應(yīng)用目的宜采用不同的布設(shè)方式,前者如圖19,后者如圖20。

    微型盒子波試驗(yàn)應(yīng)該采用炸藥激發(fā)方式,炮檢距在300~1500m范圍。該距離既不會(huì)導(dǎo)致檢波器超調(diào),又能保證初至波有較強(qiáng)的能量(較高機(jī)電比)。但是,由于安全等各方面的原因,炸藥震源在非采集作業(yè)階段不易獲得,而研制新型檢波器的過程中需要頻繁測量檢波器的性能指標(biāo),這種情況下可以用重錘等在距排列幾十米的位置激發(fā)振動(dòng)作為一個(gè)簡便的替代方法。這種方法對于某些檢測手段不完備的小型公司非常適合。該方法的優(yōu)點(diǎn)是不必采用炸藥震源,隨時(shí)隨地可以試驗(yàn),幾乎沒有成本;缺點(diǎn)是高頻段(約80Hz以上)無法標(biāo)定。高頻段標(biāo)定必須借助炸藥激發(fā)高頻信號或者室內(nèi)檢測設(shè)備來完成。

    當(dāng)前節(jié)點(diǎn)地震儀是采集裝備更新的一個(gè)熱點(diǎn)。

    圖19 新型檢波器參數(shù)標(biāo)定時(shí)布設(shè)方法

    圖20 采集作業(yè)用檢波器一致性篩選時(shí)采用的布設(shè)方法

    因?yàn)槊恳粋€(gè)節(jié)點(diǎn)儀器都是一個(gè)獨(dú)立的機(jī)電單元,所以這種新型設(shè)備需要采用圖19所示的方法來考察其一致性,特別是本底噪聲的一致性。

    如果是采集作業(yè)用檢波器一致性篩選,觀測系統(tǒng)可以采用圖20所示排列方式,炮點(diǎn)位于中心,300m以外布放檢波器,4×4為一組,組與組間隔約5m,便于區(qū)分。如果是檢波器串(如每串6個(gè)),建議每組四串檢波器(四串不至于分布范圍太大,以確保輸入信號的一致性)。以每組為單位,考察組內(nèi)各串檢波器之間的一致性。在能夠接收到較強(qiáng)震源信號的偏移距內(nèi)可將排列擴(kuò)展到數(shù)萬組(每一組的信號輸入是基本一致的)。

    當(dāng)前采集作業(yè)在用的部分檢波器使用年限較長,真實(shí)參數(shù)已經(jīng)大大偏離了標(biāo)稱值,可以采用振幅歸一化 (詳見3.3節(jié))后進(jìn)行各項(xiàng)指標(biāo)的一致性比較,按照少數(shù)服從多數(shù)的原則挑選出明顯數(shù)據(jù)錯(cuò)誤的檢波器。

    3.2 分析時(shí)窗

    輸入信號的高度一致性主要是對較強(qiáng)機(jī)械信號而言的,對隨機(jī)干擾、特別是小波長的隨機(jī)干擾一致性較弱。因此,在選取輸入信號的時(shí)候,應(yīng)選取包含初至波、強(qiáng)反射、折射、強(qiáng)相干噪聲的整個(gè)記錄時(shí)段作為分析時(shí)段,這種情況下機(jī)電比最高、能量最強(qiáng)、頻帶最寬,容易識(shí)別,可更準(zhǔn)確地由檢波器輸出數(shù)據(jù)反推檢波器的各種參數(shù)。

    3.3 振幅歸一化

    不同類型檢波器對比的時(shí)候,往往存在靈敏度、濾波響應(yīng)等參數(shù)不同的情況,此時(shí)應(yīng)該根據(jù)每種檢波器的靈敏度將原始數(shù)據(jù)轉(zhuǎn)換為速度或者加速度量綱后,再進(jìn)行比較。在靈敏度未知或者不準(zhǔn)確的情況下,可以根據(jù)20~40Hz范圍內(nèi)振幅平均值做振幅歸一化來比較檢波器的性能參數(shù)。

    如果檢波器的類型(速度、加速度)不同,要進(jìn)行轉(zhuǎn)換之后再進(jìn)行振幅歸一化[3-5]。因?yàn)椴煌瑱z波系統(tǒng)的性能參數(shù)不同,根據(jù)時(shí)域波形的最大值進(jìn)行歸一化有時(shí)候是不合理的。

    4 關(guān)于檢波器應(yīng)用與研究的幾點(diǎn)看法

    近十幾年來,以MEMS檢波器引進(jìn)國內(nèi)為標(biāo)志,國內(nèi)業(yè)界進(jìn)行了數(shù)十甚至上百種不同型號的檢波器試驗(yàn)與應(yīng)用。但是,即使在此背景下,現(xiàn)場施工仍然以傳統(tǒng)的10Hz動(dòng)圈式檢波器串為主,MEMS數(shù)字檢波器等僅占單點(diǎn)采集總道數(shù)的大約10%[14]。為什么當(dāng)年備受推崇的MEMS數(shù)字檢波器沒有像廠家預(yù)測的那樣[20]被大面積推廣?為什么部分指標(biāo)并不優(yōu)秀的檢波器進(jìn)入了日常采集作業(yè)流程?從技術(shù)應(yīng)用的角度看主要有以下幾個(gè)要素應(yīng)該被重點(diǎn)關(guān)注。

    4.1 速度/加速度

    MEMS加速度數(shù)字檢波器被引進(jìn)十幾年來,很多文獻(xiàn)[20-23]都將其與傳統(tǒng)的速度型動(dòng)圈式檢波器直接對比。因?yàn)樗俣?、加速度是不同的物理量,不同檢波系統(tǒng)的性能參數(shù)也不同,將二者直接對比沒有物理意義。將速度型檢波器與加速度型檢波器數(shù)據(jù)進(jìn)行頻譜、能量、信噪比、子波等指標(biāo)比較并據(jù)此判定優(yōu)劣,沒有反映地震數(shù)據(jù)所代表的物理本質(zhì),僅僅反映了數(shù)學(xué)表象。近年來經(jīng)過有關(guān)文獻(xiàn)的呼吁[3-6],大家逐漸接受了二者應(yīng)該在同一種“域”進(jìn)行比較的觀點(diǎn)。經(jīng)過重新認(rèn)識(shí),MEMS加速度檢波器并沒有表現(xiàn)出引進(jìn)初期期望的大大提高勘探精度的跨越式促進(jìn)作用[20-25](只是過分夸大了MEMS這種優(yōu)質(zhì)檢波器的作用)。

    但在MEMS加速度檢波器可以提高“視覺/數(shù)學(xué)”主頻[3]的啟發(fā)下,國內(nèi)部分單位轉(zhuǎn)向研發(fā)各類新型加速度檢波器,涌現(xiàn)了部分指標(biāo)并不突出甚至性能低劣的檢波器,意圖在于提高數(shù)據(jù)主頻以及地質(zhì)分辨率,并進(jìn)行了大量的試驗(yàn)甚至采集作業(yè),在一定程度上造成了資源的浪費(fèi)。并且,直接將加速度檢波器數(shù)據(jù)輸入處理系統(tǒng)進(jìn)行處理的做法,可能違背了“基本的物理原理”[3,14,20]。文獻(xiàn)[3]初步討論了選擇速度、加速度進(jìn)行處理的時(shí)候應(yīng)該考慮的幾個(gè)因素。文獻(xiàn)[14]認(rèn)為,在進(jìn)行速度、加速度檢波器對比、處理的時(shí)候,必須把加速度檢波器數(shù)據(jù)積分轉(zhuǎn)換為速度數(shù)據(jù),因?yàn)樗俣扔蚋欣诮忉屓藛T拾取準(zhǔn)確的反射時(shí)間(圖21)。對于厚層(圖21左圖,由上而下依次為速度記錄、反射系數(shù)序列、加速度記錄)而言,每一個(gè)反射層對應(yīng)區(qū)分明顯的速度或加速度子波。但對于薄層(圖21右圖)而言,速度雷克子波可區(qū)分兩個(gè)反射層的情況下(圖21右圖上)、加速度子波就難以辨別了(圖21右圖下)。因此,盡管用加速度表征較速度表征的子波主頻由30Hz提高到36Hz,但其時(shí)域分辨能力下降了。

    傳統(tǒng)的速度型動(dòng)圈式檢波器自1937年被廣泛采用以來[1-2],很多處理軟件都默認(rèn)輸入數(shù)據(jù)代表地表振動(dòng)的速度,故世界上大多數(shù)石油公司在將加速度檢波器數(shù)據(jù)輸入處理系統(tǒng)之前將其積分轉(zhuǎn)換為速度并切除低頻畸變后再進(jìn)行處理[3,14,20]。但因?yàn)椴糠旨铀俣葯z波器本底噪聲非常高,積分轉(zhuǎn)換為速度后會(huì)出現(xiàn)畸變的超強(qiáng)低頻分量(一般小于5Hz),所以國內(nèi)部分項(xiàng)目沒有將加速度檢波器數(shù)據(jù)轉(zhuǎn)換為速度數(shù)據(jù),也不愿看到加速度數(shù)據(jù)的“高頻優(yōu)勢”在積分后消失[3],所以就將其直接輸入處理系統(tǒng)。其后果是處理后的剖面視主頻提高了,但是時(shí)間分辨率和信噪比都下降了。

    另外,國內(nèi)部分學(xué)者認(rèn)為采用加速度數(shù)據(jù)可提高縱向分辨率;或者認(rèn)為直接把加速度輸入處理系統(tǒng)后,處理系統(tǒng)可以通過反褶積自動(dòng)消除速度與加速度之間的動(dòng)力學(xué)特征差異,二者的成像效果相當(dāng)、不需預(yù)先進(jìn)行微積分變換。由于作者知識(shí)領(lǐng)域的限制,無法對以上觀點(diǎn)做理論與應(yīng)用方面的深度剖析,希望引起相關(guān)專家的關(guān)注,進(jìn)行進(jìn)一步的驗(yàn)證,避免出現(xiàn)工程層面的原理性失誤。如對多個(gè)區(qū)塊的加速度數(shù)據(jù)直接進(jìn)行處理,會(huì)產(chǎn)生成像扭曲,誤導(dǎo)地質(zhì)專家的更深入認(rèn)識(shí),進(jìn)而降低鉆井成功率。如果經(jīng)過論證,證明加速度處理更有物理意義、地質(zhì)效果更好,就可將以往大量速度檢波器采集的數(shù)據(jù)進(jìn)行微分后重新處理;如果證明速度處理效果更好,那么近年來所有由加速度檢波器施工的區(qū)塊資料應(yīng)該被積分后再重新處理,而不是目前速度與加速度并行、莫衷一是的狀況。

    在實(shí)踐中,有的新型檢波器通過附加微分電路將速度型電信號整形為加速度后輸出,以達(dá)到提高視覺主頻的目的。盡管這種方法理論上對于提高高頻弱信號保真度有一定的作用,但是因?yàn)槲⒎蛛娐返姆€(wěn)定性以及地震勘探目標(biāo)頻帶范圍、附加電噪聲等原因,這種“電微分”的工程意義不大,完全可以通過數(shù)學(xué)微分替代。數(shù)學(xué)微分的效果取決于原始數(shù)據(jù)的“機(jī)電比”[5]。

    圖21 速度域雷克子波(左)較加速度域(右)更有利于提高時(shí)域分辨率

    4.2 低頻/常規(guī)

    低頻信號在克服吸收、非均質(zhì)體散射以及動(dòng)靜校正、疊加、偏移、信噪比等方面有優(yōu)勢,有助于解決復(fù)雜地區(qū)的構(gòu)造問題;深層反射、火成巖反射、鹽下層反射、復(fù)雜逆掩斷裂滑脫面下的反射具有低頻的特性。近年來,隨著專家學(xué)者對于低頻信號(1~10Hz)在石油勘探中的作用認(rèn)識(shí)不斷深化,涌現(xiàn)了多種獲取低頻信號的方法。從接收的角度而言,獲取低頻信號大致有四種方法:①采用低頻無衰減或者衰減非常小的檢波器。這種獲取低頻的方式主要依賴于檢波器的機(jī)械屬性,可稱之為“機(jī)械低頻”(如MEMS數(shù)字檢波器)。②在存在機(jī)械低頻衰減的情況下,根據(jù)檢波器的傳輸函數(shù),經(jīng)由附加電路補(bǔ)償?shù)皖l,這種方式可以稱之“電低頻”(如西方地球物理公司的GAC檢波器)。③在存在機(jī)械低頻衰減的情況下,根據(jù)檢波器的傳輸函數(shù),將損失的低頻通過計(jì)算的方式補(bǔ)償回來[5],這種方法可以稱之為“數(shù)學(xué)低頻”。 ④根據(jù)普通檢波器與低頻檢波器數(shù)據(jù)之間的特征(振幅譜、功率譜)差異,補(bǔ)償普通檢波器的低頻成分[26-27],是以上三種方法的綜合。無論何種實(shí)現(xiàn)形式,它所希望代表的都應(yīng)該是地表機(jī)械振動(dòng)中的低頻分量。

    圖22是四種檢波器貼近埋置后接收同一炮的振幅譜(經(jīng)20~40Hz振幅歸一化)。其中圖22a未做檢波器反褶積[5],可見藍(lán)—綠—青(方法②)—紅(方法①)四種檢波器的低頻振幅依次增強(qiáng),這與理論相符。圖22b的動(dòng)圈式檢波器20DX做了檢波器反褶積Deg[5](方法③)后,與DSU、寬頻檢波器在2Hz以上高度重合; LP-5Hz加速度檢波器存在低頻衰減,故在約7Hz以上才與其余檢波器重合。

    關(guān)于低頻檢波器,本文認(rèn)為:

    (1)對于石油勘探而言,DSU等MEMS數(shù)字檢波器是最直接的選擇,缺點(diǎn)是價(jià)格偏高。DSU-428大約在2Hz以上是可靠的,DSU-508應(yīng)該可以達(dá)到大約1Hz甚至更低。其余很多類別不適于大規(guī)模的石油勘探作業(yè);同時(shí),此類加速度檢波器的低頻優(yōu)勢只有在轉(zhuǎn)換為速度進(jìn)行數(shù)據(jù)處理的時(shí)候才能顯現(xiàn)(限于最低可靠頻率以上)。

    圖22 四種不同類型檢波器的低頻接收能力比較

    (2)采用十幾赫茲的動(dòng)圈式檢波器+反饋電路實(shí)現(xiàn)1~2Hz以上低頻信號接收的情況下,因反饋電路難以保持較高的一致性,故疊加后低頻數(shù)據(jù)的保真度較低。

    (3)10Hz動(dòng)圈式檢波器+檢波器反褶積Deg[5]大約可以實(shí)現(xiàn)1~2Hz以上的數(shù)據(jù)保真,是最經(jīng)濟(jì)的一種方法。更重要的是,Deg可用于以往所有動(dòng)圈式檢波器采集的工區(qū)(已施加低截濾波的除外),可對以往數(shù)據(jù)進(jìn)行重新處理、解釋,發(fā)掘低頻分量的潛力,大大降低采用低頻檢波器重新勘探的成本[5]。有技術(shù)人員認(rèn)為可通過處理系統(tǒng)的反褶積來代替檢波器反褶積[5],這是不現(xiàn)實(shí)的,因?yàn)槎叩乃罁?jù)的物理事實(shí)以及實(shí)現(xiàn)算法均不相同。

    (4)方法④依賴于低頻檢波器的機(jī)械或者電屬性以及數(shù)學(xué)計(jì)算(反褶積),算子不穩(wěn)定,一致性差,操作繁瑣,推廣意義較弱。

    4.3 單點(diǎn)/組合

    “高密度采集”是指空間采樣密度遠(yuǎn)大于常規(guī)施工的采集方式,往往道距較小,主要目的是提高地震資料的信噪比和分辨率。當(dāng)前,高密度采集在很多地區(qū)被大量實(shí)施,且少量施工單位認(rèn)為其必須同時(shí)與單點(diǎn)加速度檢波器相配合。但從分類角度看,高密度本質(zhì)屬“觀測系統(tǒng)”的范疇,單點(diǎn)采集屬“接收”的范疇,加速度檢波器則屬“裝備”的范疇,不能混為一談。部分高密度采集實(shí)例中的剖面品質(zhì)提升實(shí)際是加速度檢波器接收帶來的“視覺/數(shù)學(xué)”高頻優(yōu)勢,如果將其積分為速度后,高頻優(yōu)勢就消失了,時(shí)間分辨率不會(huì)有明顯提高。當(dāng)然,因?yàn)榈谰嗟目s小,空間分辨率有了較大改善。

    文獻(xiàn)[14]認(rèn)為,在低信噪比地區(qū),因?yàn)槿藗內(nèi)匀灰蕾囉趩闻陲@示進(jìn)行質(zhì)量控制,通過組合濾波響應(yīng)衰減低速炮生干擾仍然是人們優(yōu)先的選擇。只有當(dāng)單點(diǎn)檢波器可以在Inline 和Crossline 方向小到足夠無假頻記錄波場并在處理階段有效去噪的時(shí)候,才能夠替代檢波器組合。這在當(dāng)下油價(jià)低迷的現(xiàn)實(shí)條件下很難實(shí)現(xiàn),特別是Crossline方向無法達(dá)到所期望的空間密度。因此,在很多地區(qū),檢波器組合仍是當(dāng)前低成本采集背景下的合理選擇。

    另外,在實(shí)際地震數(shù)據(jù)采集中,將不同檢波器的單點(diǎn)與組合進(jìn)行同時(shí)比對是錯(cuò)誤的,此時(shí)檢波器的濾波效應(yīng)與組合效應(yīng)相互重合,難以得出明晰的結(jié)論。

    當(dāng)前,一些單位熱衷于研發(fā)新型加速度檢波器及低頻檢波器的根本原因,可能是對檢波器性能指標(biāo)與其地球物理表現(xiàn)之間的因果關(guān)系不清楚或商業(yè)利益驅(qū)動(dòng),而對提高地震數(shù)據(jù)保真度并無太大助益。

    除了以上純技術(shù)因素以外,一些更深層面的因素影響了檢波器的應(yīng)用與研究工作。

    (1)部分一線技術(shù)人員缺乏充分理解不同類型檢波器數(shù)據(jù)所代表的物理意涵以及應(yīng)用場景的知識(shí)基礎(chǔ)和工具方法。

    (2)在勘探要素發(fā)生變化的情況下,數(shù)據(jù)評價(jià)沿襲了技術(shù)規(guī)程中的僵化方法、導(dǎo)致了錯(cuò)誤的判斷。檢波器效能應(yīng)該放置到整個(gè)勘探鏈條中進(jìn)行系統(tǒng)評判,而不僅僅是根據(jù)技術(shù)規(guī)程做“現(xiàn)象評價(jià)”[3-4]。

    (3)以“工業(yè)試錯(cuò)”模式開展應(yīng)用研究。在未厘清數(shù)理原理或未做充分點(diǎn)試驗(yàn)情況下,倉促投入大型試驗(yàn)甚至規(guī)模作業(yè),浪費(fèi)了資源、誤導(dǎo)了認(rèn)識(shí)。

    (4)對MEMS等新型檢波器的評價(jià)缺乏獨(dú)立判斷、“只有表揚(yáng)沒有批評”,鮮有“不同的聲音”。這種現(xiàn)象與求真務(wù)實(shí)的科研作風(fēng)是相悖的,背后的方法性、結(jié)構(gòu)性甚至系統(tǒng)性的原因值得反思。

    5 結(jié)束語

    (1)由陸上震源激發(fā)的振動(dòng)在1m×1m范圍內(nèi)的動(dòng)力學(xué)特征(振幅、頻率、相位)高度一致,越趨向低頻一致性越強(qiáng)。

    (2)微型盒子波是一種基于現(xiàn)場試驗(yàn)的、簡便的估測檢波器性能指標(biāo)的方法,可以估測靈敏度、最低可靠頻率、允差以及綜合一致性等,分為常規(guī)采集作業(yè)和小型試驗(yàn)兩種情況。其結(jié)果可以為檢波器提供一個(gè)基于最終數(shù)據(jù)意義的評價(jià)標(biāo)準(zhǔn)。

    (3)在研制新型檢波器的時(shí)候,建議首先進(jìn)行極低成本的微型盒子波試驗(yàn);在確認(rèn)細(xì)部特征確有提高后,再進(jìn)行大規(guī)模現(xiàn)場試驗(yàn)。

    (4)更新失真嚴(yán)重的老舊檢波器是當(dāng)前采集設(shè)備管理的主要任務(wù),不是研制新型檢波器。10Hz動(dòng)圈式檢波器在當(dāng)前處理框架下仍是“高效”檢波器。

    (5)成熟的速度、加速度檢波器對于地表機(jī)械振動(dòng)的表征作用是等效的。到底是速度表征還是加速度表征更有助于提高勘探效果,有待處理、解釋階段的深入系統(tǒng)驗(yàn)證。

    (6)檢波器反褶積[5]可以穩(wěn)定地恢復(fù)1~2Hz以上的低頻信號。低頻恢復(fù)后的地震信號在處理、解釋階段的應(yīng)用有待驗(yàn)證。

    (7)在充分認(rèn)識(shí)檢波器工作原理及其數(shù)據(jù)表現(xiàn)的基礎(chǔ)上,重新制定檢波器選型試驗(yàn)的數(shù)據(jù)評價(jià)標(biāo)準(zhǔn)。

    感謝CGGVeritas公司Shuki Ronen博士、Seismic Oilfield Service公司Robert Heath先生及Sercel公司Jérme Lainé先生的經(jīng)驗(yàn)介紹及觀點(diǎn)分享。

    猜你喜歡
    檢波器振幅靈敏度
    檢波器容差對地震信號接收的影響研究
    導(dǎo)磁環(huán)對LVDT線性度和靈敏度的影響
    一種井下檢波器測試儀的設(shè)計(jì)
    電子測試(2018年6期)2018-05-09 07:31:34
    基于高低頻聯(lián)測的常規(guī)檢波器數(shù)據(jù)低頻振幅和相位同時(shí)恢復(fù)方法
    地下水非穩(wěn)定流的靈敏度分析
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    穿甲爆破彈引信對薄弱目標(biāo)的靈敏度分析
    天堂动漫精品| 成人高潮视频无遮挡免费网站| 日韩大尺度精品在线看网址| 国产精品不卡视频一区二区| 在线观看av片永久免费下载| 久久中文看片网| 亚洲七黄色美女视频| 性色avwww在线观看| 欧美另类亚洲清纯唯美| 亚洲国产精品合色在线| 91av网一区二区| 亚洲无线在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲不卡免费看| 两性午夜刺激爽爽歪歪视频在线观看| 色在线成人网| 久久久精品欧美日韩精品| 亚洲欧美日韩高清专用| 国产伦精品一区二区三区四那| 欧美激情久久久久久爽电影| ponron亚洲| 亚洲精品成人久久久久久| 听说在线观看完整版免费高清| 狂野欧美白嫩少妇大欣赏| 色尼玛亚洲综合影院| 国产毛片a区久久久久| 亚洲成人精品中文字幕电影| 欧美一区二区亚洲| 啦啦啦啦在线视频资源| 日韩精品中文字幕看吧| 日韩精品有码人妻一区| 国产精品久久视频播放| 亚洲人成网站在线播| 亚洲,欧美,日韩| 国产精品久久久久久亚洲av鲁大| 丰满乱子伦码专区| 在线免费观看不下载黄p国产| 人妻久久中文字幕网| 又黄又爽又刺激的免费视频.| 国产成人一区二区在线| 亚洲最大成人中文| 男女做爰动态图高潮gif福利片| 女人十人毛片免费观看3o分钟| 在线观看av片永久免费下载| 亚洲精品色激情综合| 亚洲欧美日韩卡通动漫| 听说在线观看完整版免费高清| 亚洲熟妇熟女久久| 深夜a级毛片| 久久久久久九九精品二区国产| 精品人妻偷拍中文字幕| 人妻夜夜爽99麻豆av| 日日啪夜夜撸| 一级a爱片免费观看的视频| 色噜噜av男人的天堂激情| 国产一区二区在线观看日韩| 丰满的人妻完整版| 久久婷婷人人爽人人干人人爱| 男人舔女人下体高潮全视频| 少妇人妻一区二区三区视频| 搡女人真爽免费视频火全软件 | 成人漫画全彩无遮挡| 国产精品爽爽va在线观看网站| 两个人视频免费观看高清| videossex国产| 99热这里只有是精品在线观看| 国产伦精品一区二区三区四那| 国内少妇人妻偷人精品xxx网站| 亚洲av成人av| 一级黄片播放器| 亚洲无线观看免费| 日本色播在线视频| 欧美在线一区亚洲| 国产在视频线在精品| 在线a可以看的网站| 最新中文字幕久久久久| 在线观看午夜福利视频| 久久久a久久爽久久v久久| 国产午夜福利久久久久久| 波多野结衣高清作品| 午夜福利视频1000在线观看| 免费人成在线观看视频色| 国内精品美女久久久久久| 99久久精品热视频| 国产亚洲av嫩草精品影院| 国产在线精品亚洲第一网站| 成人av一区二区三区在线看| 国产精品爽爽va在线观看网站| 夜夜夜夜夜久久久久| 亚洲精华国产精华液的使用体验 | 欧美性猛交黑人性爽| 日韩 亚洲 欧美在线| 99精品在免费线老司机午夜| 自拍偷自拍亚洲精品老妇| 国产精品精品国产色婷婷| 麻豆精品久久久久久蜜桃| 18+在线观看网站| 日本五十路高清| 我要看日韩黄色一级片| 国产精品伦人一区二区| 色综合色国产| 久久久久久久久久黄片| 午夜爱爱视频在线播放| 国产精品亚洲美女久久久| 国产真实伦视频高清在线观看| 天美传媒精品一区二区| 精品久久久久久久久久久久久| 女同久久另类99精品国产91| 国产一区二区激情短视频| 国产黄色视频一区二区在线观看 | 久久久a久久爽久久v久久| 人妻夜夜爽99麻豆av| 九九热线精品视视频播放| 国产中年淑女户外野战色| 欧美激情国产日韩精品一区| 国产高清激情床上av| 中国美女看黄片| 观看免费一级毛片| 国产三级在线视频| 国模一区二区三区四区视频| 少妇的逼水好多| 村上凉子中文字幕在线| 欧美zozozo另类| 亚洲第一电影网av| 中文字幕av成人在线电影| 亚洲av免费在线观看| 国产淫片久久久久久久久| 赤兔流量卡办理| 2021天堂中文幕一二区在线观| 男女之事视频高清在线观看| 国产精品一区二区免费欧美| 亚洲真实伦在线观看| 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| 欧美+亚洲+日韩+国产| 99九九线精品视频在线观看视频| 亚洲美女搞黄在线观看 | 午夜激情福利司机影院| 婷婷亚洲欧美| 日韩欧美一区二区三区在线观看| 少妇人妻精品综合一区二区 | 国产熟女欧美一区二区| 噜噜噜噜噜久久久久久91| 成人二区视频| 亚洲第一区二区三区不卡| 有码 亚洲区| 99热这里只有是精品在线观看| 99久久九九国产精品国产免费| 99热网站在线观看| 亚洲欧美中文字幕日韩二区| 黄色一级大片看看| 变态另类丝袜制服| 欧洲精品卡2卡3卡4卡5卡区| 大型黄色视频在线免费观看| 日韩欧美免费精品| 精品久久久久久久久久久久久| 国产女主播在线喷水免费视频网站 | 在线观看美女被高潮喷水网站| 亚洲av不卡在线观看| 女同久久另类99精品国产91| 一级毛片久久久久久久久女| 九九爱精品视频在线观看| 精品一区二区三区人妻视频| 亚洲高清免费不卡视频| av卡一久久| 国产精品永久免费网站| 美女被艹到高潮喷水动态| 日韩成人伦理影院| a级一级毛片免费在线观看| 色综合色国产| 自拍偷自拍亚洲精品老妇| 精品久久国产蜜桃| 日韩亚洲欧美综合| 午夜精品一区二区三区免费看| 欧美日韩综合久久久久久| 99久久久亚洲精品蜜臀av| 成人精品一区二区免费| 国产亚洲精品久久久久久毛片| av在线观看视频网站免费| 国产成年人精品一区二区| 久久6这里有精品| 老司机福利观看| 欧美激情久久久久久爽电影| 国产av在哪里看| 国产一区亚洲一区在线观看| 变态另类丝袜制服| 欧美激情在线99| 直男gayav资源| 黄色一级大片看看| 日韩欧美精品v在线| 欧美性感艳星| 久久精品国产99精品国产亚洲性色| 51国产日韩欧美| 精品久久久久久久久久免费视频| 免费av不卡在线播放| 成人av一区二区三区在线看| 亚洲综合色惰| 变态另类丝袜制服| 可以在线观看的亚洲视频| 老师上课跳d突然被开到最大视频| 亚洲av五月六月丁香网| 亚洲乱码一区二区免费版| 国产精品国产三级国产av玫瑰| 国产欧美日韩一区二区精品| 免费看av在线观看网站| 国产伦一二天堂av在线观看| 搡老熟女国产l中国老女人| 亚洲av中文av极速乱| 国产精品爽爽va在线观看网站| 久久午夜亚洲精品久久| 久久精品国产鲁丝片午夜精品| 亚洲中文日韩欧美视频| 国产亚洲精品av在线| 在线a可以看的网站| av国产免费在线观看| 在线观看美女被高潮喷水网站| 亚洲一区高清亚洲精品| 内射极品少妇av片p| 日本色播在线视频| 综合色丁香网| 免费看av在线观看网站| 久久久久久久久久成人| 亚洲乱码一区二区免费版| 午夜亚洲福利在线播放| 蜜桃亚洲精品一区二区三区| 国产精品亚洲美女久久久| 国产精品无大码| 午夜亚洲福利在线播放| 成年av动漫网址| 寂寞人妻少妇视频99o| 欧美性猛交黑人性爽| 久久99热6这里只有精品| 国产激情偷乱视频一区二区| 九九久久精品国产亚洲av麻豆| 国产精品美女特级片免费视频播放器| 亚洲人成网站在线播放欧美日韩| 亚洲aⅴ乱码一区二区在线播放| 99热精品在线国产| 久久久久久伊人网av| 麻豆乱淫一区二区| 欧美日韩一区二区视频在线观看视频在线 | 1000部很黄的大片| 午夜免费激情av| 伊人久久精品亚洲午夜| 综合色av麻豆| 又黄又爽又刺激的免费视频.| 18禁黄网站禁片免费观看直播| 日日摸夜夜添夜夜添av毛片| 午夜a级毛片| 亚洲人与动物交配视频| av.在线天堂| 欧美日韩乱码在线| 赤兔流量卡办理| 久久中文看片网| 晚上一个人看的免费电影| 大型黄色视频在线免费观看| 男人的好看免费观看在线视频| 欧美+日韩+精品| 女生性感内裤真人,穿戴方法视频| 99久久成人亚洲精品观看| av天堂在线播放| 日本一二三区视频观看| 国产精品一区二区三区四区免费观看 | 日本三级黄在线观看| 日韩人妻高清精品专区| 性色avwww在线观看| 日韩,欧美,国产一区二区三区 | 国产私拍福利视频在线观看| 99久国产av精品| 亚洲经典国产精华液单| 综合色丁香网| 91狼人影院| 亚洲熟妇中文字幕五十中出| 欧美成人免费av一区二区三区| 中文字幕av成人在线电影| 国产成人影院久久av| 国产成人freesex在线 | 婷婷精品国产亚洲av| 舔av片在线| 淫妇啪啪啪对白视频| 久久这里只有精品中国| 亚洲久久久久久中文字幕| 男人舔奶头视频| 国产成人影院久久av| 日韩精品中文字幕看吧| а√天堂www在线а√下载| 国产伦精品一区二区三区四那| 色吧在线观看| 亚洲av电影不卡..在线观看| www日本黄色视频网| 三级毛片av免费| av在线蜜桃| 18禁在线无遮挡免费观看视频 | 91久久精品国产一区二区成人| 欧美日本视频| 精品乱码久久久久久99久播| 波多野结衣高清无吗| 最后的刺客免费高清国语| 两性午夜刺激爽爽歪歪视频在线观看| 我要搜黄色片| 国产精品永久免费网站| 九九在线视频观看精品| 欧美3d第一页| 久久久精品大字幕| 午夜福利18| 欧美在线一区亚洲| 久久久久精品国产欧美久久久| 99视频精品全部免费 在线| 国产精品,欧美在线| 精品欧美国产一区二区三| 国产在线男女| 国产 一区精品| 亚洲va在线va天堂va国产| 看片在线看免费视频| 97人妻精品一区二区三区麻豆| 日韩av在线大香蕉| 久久久a久久爽久久v久久| 亚洲精品亚洲一区二区| 日日干狠狠操夜夜爽| 精品一区二区三区av网在线观看| 男女啪啪激烈高潮av片| 欧美激情在线99| 国产色婷婷99| 日日啪夜夜撸| 毛片一级片免费看久久久久| 天堂av国产一区二区熟女人妻| 97超视频在线观看视频| 老师上课跳d突然被开到最大视频| 日本 av在线| 一级av片app| 亚洲精品在线观看二区| 亚洲成a人片在线一区二区| 不卡一级毛片| 简卡轻食公司| 97碰自拍视频| 亚洲成a人片在线一区二区| 中文字幕av成人在线电影| av女优亚洲男人天堂| 亚洲中文字幕日韩| 国产精品嫩草影院av在线观看| 精品不卡国产一区二区三区| 男插女下体视频免费在线播放| 精品福利观看| 最近中文字幕高清免费大全6| 成人欧美大片| 麻豆国产97在线/欧美| 啦啦啦观看免费观看视频高清| 亚洲综合色惰| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 最近最新中文字幕大全电影3| 久久久成人免费电影| 亚洲不卡免费看| 成人国产麻豆网| 精品久久久久久成人av| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品影院6| 老师上课跳d突然被开到最大视频| 亚洲欧美成人综合另类久久久 | 美女 人体艺术 gogo| 久久久久久九九精品二区国产| 好男人在线观看高清免费视频| 插逼视频在线观看| 给我免费播放毛片高清在线观看| av.在线天堂| 亚洲欧美清纯卡通| 国产免费一级a男人的天堂| 久久午夜亚洲精品久久| 国内少妇人妻偷人精品xxx网站| 国产人妻一区二区三区在| 国产免费一级a男人的天堂| www日本黄色视频网| 极品教师在线视频| 欧美bdsm另类| 亚洲图色成人| 国产精品久久久久久精品电影| 又爽又黄a免费视频| 亚洲av.av天堂| 高清毛片免费看| 久久久成人免费电影| 99久久成人亚洲精品观看| 久久久久久久久大av| 久久久a久久爽久久v久久| 亚洲精品色激情综合| 亚洲欧美日韩高清在线视频| 长腿黑丝高跟| 真人做人爱边吃奶动态| 成人漫画全彩无遮挡| 欧美激情国产日韩精品一区| 国产高清有码在线观看视频| 免费av不卡在线播放| 欧美成人一区二区免费高清观看| 国产极品精品免费视频能看的| 国产成人影院久久av| 国产午夜福利久久久久久| 国产欧美日韩一区二区精品| 波多野结衣高清无吗| 夜夜夜夜夜久久久久| 亚洲精品日韩av片在线观看| 日本在线视频免费播放| 亚洲自偷自拍三级| 亚洲中文日韩欧美视频| 搡老妇女老女人老熟妇| 亚洲专区国产一区二区| 亚洲欧美清纯卡通| 精品人妻偷拍中文字幕| 亚洲一区高清亚洲精品| 热99在线观看视频| 久久久久精品国产欧美久久久| 久久久久久久久大av| 国产成人a区在线观看| 亚洲内射少妇av| 在线国产一区二区在线| 久久综合国产亚洲精品| 黄色配什么色好看| 中国美女看黄片| 永久网站在线| 非洲黑人性xxxx精品又粗又长| 欧美日本亚洲视频在线播放| 免费黄网站久久成人精品| 久久久久久伊人网av| 成人美女网站在线观看视频| 久久久久久国产a免费观看| 91在线观看av| 久久综合国产亚洲精品| 97超碰精品成人国产| 欧美精品国产亚洲| 亚洲av.av天堂| 一个人免费在线观看电影| 成人性生交大片免费视频hd| 国产免费男女视频| 男女边吃奶边做爰视频| av在线老鸭窝| eeuss影院久久| 国产中年淑女户外野战色| 亚洲第一区二区三区不卡| 一区二区三区免费毛片| 校园春色视频在线观看| 天天躁夜夜躁狠狠久久av| 男人的好看免费观看在线视频| 少妇高潮的动态图| 深爱激情五月婷婷| 精品国内亚洲2022精品成人| 在线观看一区二区三区| 秋霞在线观看毛片| 国产高清视频在线播放一区| 俺也久久电影网| 色综合色国产| 欧美日韩乱码在线| 亚洲不卡免费看| 亚洲人成网站在线播放欧美日韩| 91在线精品国自产拍蜜月| 国产探花极品一区二区| 亚洲婷婷狠狠爱综合网| 12—13女人毛片做爰片一| 欧美色欧美亚洲另类二区| 3wmmmm亚洲av在线观看| 亚洲美女搞黄在线观看 | 69av精品久久久久久| 亚洲成人av在线免费| 99在线视频只有这里精品首页| 人人妻,人人澡人人爽秒播| 女人被狂操c到高潮| 男人舔女人下体高潮全视频| 黄色欧美视频在线观看| 欧美一区二区亚洲| 小说图片视频综合网站| 欧美xxxx性猛交bbbb| 成人三级黄色视频| 国产在视频线在精品| 色哟哟·www| 精品少妇黑人巨大在线播放 | 免费大片18禁| 日韩精品中文字幕看吧| 一本一本综合久久| 亚洲精品一区av在线观看| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久久亚洲| 高清毛片免费看| 长腿黑丝高跟| 精品少妇黑人巨大在线播放 | 少妇猛男粗大的猛烈进出视频 | 免费电影在线观看免费观看| 午夜免费男女啪啪视频观看 | 久久精品91蜜桃| 国产精品美女特级片免费视频播放器| 日韩中字成人| 免费观看精品视频网站| 激情 狠狠 欧美| 日韩精品中文字幕看吧| 99在线视频只有这里精品首页| 亚洲乱码一区二区免费版| 国产乱人偷精品视频| 久久人人爽人人爽人人片va| 伊人久久精品亚洲午夜| 99久国产av精品国产电影| 亚洲精品一区av在线观看| 国内久久婷婷六月综合欲色啪| 国产亚洲91精品色在线| 不卡视频在线观看欧美| 午夜免费激情av| 亚洲精品在线观看二区| 久久久久精品国产欧美久久久| 激情五月婷婷亚洲| 色94色欧美一区二区| 亚洲av免费高清在线观看| 日韩av不卡免费在线播放| 亚洲人成网站在线播| 九色成人免费人妻av| 伦精品一区二区三区| 亚洲精品国产av成人精品| 一级,二级,三级黄色视频| 亚洲人成网站在线播| 日韩av免费高清视频| 欧美日本中文国产一区发布| 99久久精品国产国产毛片| 国产精品偷伦视频观看了| 亚洲图色成人| 亚洲内射少妇av| 国产欧美亚洲国产| 国产精品久久久久久精品电影小说| 日本与韩国留学比较| 亚洲av在线观看美女高潮| 一级毛片电影观看| av福利片在线观看| 麻豆乱淫一区二区| 最近最新中文字幕免费大全7| 亚洲欧美成人精品一区二区| 18禁在线无遮挡免费观看视频| 91午夜精品亚洲一区二区三区| 亚洲国产成人一精品久久久| 国产精品免费大片| 三级经典国产精品| 亚洲欧美日韩另类电影网站| 亚洲欧美成人综合另类久久久| a级一级毛片免费在线观看| 国产免费视频播放在线视频| 日韩一本色道免费dvd| 亚洲不卡免费看| 精品一区在线观看国产| 人体艺术视频欧美日本| 麻豆乱淫一区二区| 丰满人妻一区二区三区视频av| 日韩在线高清观看一区二区三区| 人妻人人澡人人爽人人| 黄片无遮挡物在线观看| 黄色欧美视频在线观看| 免费观看无遮挡的男女| 中文天堂在线官网| 亚洲人成网站在线播| 欧美另类一区| 最近中文字幕2019免费版| 国产视频内射| 99久久精品热视频| 久久亚洲国产成人精品v| 春色校园在线视频观看| 69精品国产乱码久久久| 在线观看免费视频网站a站| 国产成人a∨麻豆精品| 男人舔奶头视频| 精品亚洲成a人片在线观看| 久久综合国产亚洲精品| 久久99一区二区三区| 狠狠精品人妻久久久久久综合| 欧美另类一区| 99视频精品全部免费 在线| 免费看av在线观看网站| 高清av免费在线| 三级经典国产精品| 少妇人妻精品综合一区二区| 成人国产av品久久久| 亚洲国产精品一区二区三区在线| 亚洲国产成人一精品久久久| a级毛片免费高清观看在线播放| 久久女婷五月综合色啪小说| 国产精品一二三区在线看| 大码成人一级视频| 人妻制服诱惑在线中文字幕| 51国产日韩欧美| 成人国产av品久久久| 国产极品天堂在线| 丁香六月天网| 久久精品久久精品一区二区三区| a级毛片在线看网站| 涩涩av久久男人的天堂| 久久人人爽av亚洲精品天堂| 欧美xxⅹ黑人| 九九爱精品视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 久久久亚洲精品成人影院| 精品亚洲成a人片在线观看| 午夜免费鲁丝| 交换朋友夫妻互换小说| 97在线人人人人妻| 久久99热这里只频精品6学生| 久久免费观看电影| 一个人免费看片子| 国产黄片视频在线免费观看| 亚洲av欧美aⅴ国产| 99久久人妻综合| 成人漫画全彩无遮挡| 我的老师免费观看完整版| 久久精品国产a三级三级三级| 人人澡人人妻人| 国内少妇人妻偷人精品xxx网站| 中文字幕人妻熟人妻熟丝袜美| 久久精品久久精品一区二区三区| 亚洲精品视频女| av在线app专区| 亚洲色图综合在线观看| 国产黄色免费在线视频| 国产成人精品一,二区| 嫩草影院入口| 亚洲熟女精品中文字幕| 日韩一区二区视频免费看|