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

    立體氣槍陣列延遲激發(fā)震源特性及在淺海區(qū)OBS探測(cè)中的應(yīng)用

    2016-07-28 09:19:07吳志強(qiáng)郝天珧唐松華張訓(xùn)華鄭彥鵬肖國林趙維娜
    地球物理學(xué)報(bào) 2016年7期

    吳志強(qiáng), 郝天珧, 唐松華, 張訓(xùn)華, 鄭彥鵬, 肖國林, 趙維娜

    1 國土資源部海洋油氣資源與環(huán)境地質(zhì)重點(diǎn)實(shí)驗(yàn)室, 青島 266071 2 青島海洋地質(zhì)研究所, 青島 266071 3 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室, 青島 266000 4 中國科學(xué)院地質(zhì)與地球物理研究所, 北京 100029 5 上海海洋石油局第一海洋地質(zhì)調(diào)查大隊(duì), 上海 201200 6 國家海洋局第一海洋研究所, 青島 266061 7 中國海洋大學(xué), 青島 266063

    ?

    立體氣槍陣列延遲激發(fā)震源特性及在淺海區(qū)OBS探測(cè)中的應(yīng)用

    吳志強(qiáng)1,2,3, 郝天珧4, 唐松華5, 張訓(xùn)華1,2,3, 鄭彥鵬6, 肖國林1,2,3, 趙維娜7,2,3

    1 國土資源部海洋油氣資源與環(huán)境地質(zhì)重點(diǎn)實(shí)驗(yàn)室, 青島2660712 青島海洋地質(zhì)研究所, 青島2660713 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室, 青島2660004 中國科學(xué)院地質(zhì)與地球物理研究所, 北京1000295 上海海洋石油局第一海洋地質(zhì)調(diào)查大隊(duì), 上海2012006 國家海洋局第一海洋研究所, 青島2660617 中國海洋大學(xué), 青島266063

    摘要利用中、小容量氣槍組成的立體氣槍陣列延遲激發(fā)震源和海底地震儀(OBS)在我國北部淺海海域開展了人工地震深部地球物理探測(cè)試驗(yàn).基于水深條件和壓制水體虛反射、提升低頻能量的需要,使氣槍震源有足夠的輸出能量和高品質(zhì)子波特性,研究了立體氣槍陣列延遲激發(fā)震源工作機(jī)理,經(jīng)遠(yuǎn)場(chǎng)子波理論模擬優(yōu)選了組合參數(shù)并進(jìn)行了海上試驗(yàn)工作.結(jié)果表明,中、小容量氣槍組成的立體氣槍陣列延遲激發(fā)震源,適應(yīng)了淺水海域的激發(fā)環(huán)境,降低了由虛反射造成的局部陷波和干擾作用,有效地改善了OBS信號(hào)的品質(zhì),獲得了Ps,Pg,PmP,Pn 等多種震相.創(chuàng)新了由中、小容量氣槍組成的立體氣槍陣列延遲激發(fā)震源在淺海區(qū)OBS探測(cè)中的應(yīng)用,也填補(bǔ)了南黃海海域深地震探測(cè)數(shù)據(jù)的空白,為南黃海、渤海深部地殼結(jié)構(gòu)研究及含油氣盆地形成演化研究提供了重要的基礎(chǔ)資料.

    關(guān)鍵詞立體氣槍陣列; 遠(yuǎn)場(chǎng)子波; OBS深部探測(cè); 虛反射壓制; 淺海區(qū)

    1引言

    氣槍震源具有激發(fā)地震信號(hào)一致性好,定位精確高、重復(fù)性好、綠色環(huán)保和成本低的優(yōu)勢(shì),成為了海洋地震勘探的主流震源,國際上95%以上的海洋主動(dòng)源地震勘探采用了氣槍震源(唐杰等,2009).眾所周知,氣槍或氣槍陣列的容量越高,激發(fā)的地震波能量越大、主頻越低,探測(cè)的深度越大,但地震記錄的分辨率較低;反之,則激發(fā)的地震波能量低、主頻高,探測(cè)的深度小,地震記錄的分辨率高.因此,在淺部高分辨率探測(cè)中,基本上使用小容量氣槍或氣槍陣列;在海洋石油勘探中,一般采用中、小容量氣槍組成氣槍陣列震源,其探測(cè)深度可達(dá)數(shù)千米.

    海洋深部地震探測(cè)要求氣槍震源激發(fā)的地震信號(hào)主頻低、能量大,傳播距離在100km以上、探測(cè)深度應(yīng)在30 km以上.為了保障氣槍震源具有足夠的輸出能量,一般使用4~6條大容量氣槍組成氣槍陣列震源(羅桂純等,2006;林建民等,2008;丘學(xué)林等,2007;趙明輝等,2008;支鵬遙,2012).國內(nèi)外有過許多以大容量氣槍陣列作為震源的深部地震探測(cè)試驗(yàn),如美國的“LARSE”計(jì)劃(Nazareth,1994),新西蘭的“TAICRUST”計(jì)劃(McIntosh,1996),我國的南海深部探測(cè)項(xiàng)目(趙明輝等,2004; 吳振利等,2008;呂川川等,2011;丘學(xué)林等,2012;)和渤海深部探測(cè)項(xiàng)目(支鵬遙等,2012),均獲得了深達(dá)莫霍面的有效反射與折射震相,為深部地殼結(jié)構(gòu)研究打下了堅(jiān)實(shí)基礎(chǔ).

    到目前為止,OBS深部探測(cè)均使用的大容量平面氣槍陣列震源,即將組成氣槍陣列震源的所有氣槍均沉放在海平面之下相同深度下,它具有操作方便、排列簡(jiǎn)單和能實(shí)現(xiàn)子陣列最大能量同時(shí)疊加等優(yōu)點(diǎn).但是,大容量平面氣槍陣列震源也存在較大的缺點(diǎn),由于大容量氣槍的氣泡半徑大,氣槍必須沉放在大于氣泡半徑的深度,才能保證其最佳的輸出能量;氣槍的沉放深度大,氣槍到海平面之間的水體虛反射(鬼波)改變了氣槍震源激發(fā)的地震子波特征,地震信號(hào)的低頻分量和高頻分量受到了不同程度的抑制,降低了地震數(shù)據(jù)的頻帶寬度,也造成了氣槍陣列震源的遠(yuǎn)場(chǎng)子波的低頻部分振幅振蕩跳躍,頻譜曲線呈“鋸齒”狀分布,降低了原始地震資料品質(zhì)(Cambois等, 2009).另一方面,大容量的單槍及簡(jiǎn)單的陣列組合,在抑制氣泡方面能力有限,而氣泡抑制程度決定了原始地震資料的品質(zhì)(陳浩林等,2008).

    位于我國北部海域的南黃海、渤海水深只有15~40 m,在進(jìn)行OBS深部探測(cè)施工中,為保證氣槍震源的最佳輸出能量,將大容量氣槍沉放較深而與海底距離較小時(shí),氣槍瞬間釋放的2000psi高壓氣體造成了海底泥沙泛起,隨氣泡收縮進(jìn)入氣槍槍膛,引起氣槍的故障率大幅度上升,造成資料采集工作效率低下和成本上升.因此,在淺水區(qū)進(jìn)行OBS探測(cè),使用大容量氣槍陣列震源存在資料品質(zhì)差和成本增加等問題.其次,在淺水區(qū)大容量氣槍震源激發(fā)的地震波,在海底與海平面之間形成了強(qiáng)能量周期震蕩干擾,與折射/反射震相混合在一起,極大地降低了地震資料的品質(zhì).

    中、小容量氣槍的氣泡半徑較小,較淺的沉放深度就能發(fā)揮其最大輸出能量,適應(yīng)水深較淺的勘探環(huán)境,把大容量氣槍改成多條中、小容量氣槍有機(jī)相干組合,在總?cè)萘肯喈?dāng)?shù)那闆r下,不但輸出能量增加,還能能削弱氣泡振蕩,提高初泡比(楊懷春等,2004),因此,采用多條中、小容量氣槍相干組合成陣列震源,不但能改善激發(fā)子波的品質(zhì),同時(shí)也能達(dá)到與大容量氣槍相當(dāng)?shù)妮敵瞿芰?但是,中、小容量氣槍激發(fā)的地震波主頻高,低頻部分的能量弱,常規(guī)的組合方式對(duì)低頻能量起到了壓制作用.因此,使用與水深條件相適宜的中、小容量的氣槍進(jìn)行OBS深部地震探測(cè),不能機(jī)械地照搬常規(guī)的組合與激發(fā)方式,必須對(duì)其進(jìn)行改進(jìn),提高低頻輸出能量,改善遠(yuǎn)場(chǎng)子波頻譜特性,以形成遠(yuǎn)場(chǎng)子波品質(zhì)優(yōu)、輸出能量強(qiáng)、穿透深度大的氣槍陣列震源,滿足淺水區(qū)深部OBS探測(cè)的需要.

    針對(duì)南黃海、渤海OBS深部探測(cè)目標(biāo)和存在的施工難題,經(jīng)過對(duì)氣槍性能的分析研究,確定以使用與水深條件相適宜的中、小容量的氣槍組成氣槍陣列震源,達(dá)到提高施工效率和降低故障率的要求.同時(shí),采用立體氣槍陣列延遲激發(fā)的技術(shù)方法,著力提高氣槍陣列震源低頻輸出能量.為了達(dá)到上述目標(biāo),開展了立體氣槍陣列延遲激發(fā)技術(shù)的研究并對(duì)震源遠(yuǎn)場(chǎng)子波特征進(jìn)行了理論模擬計(jì)算和分析.在此基礎(chǔ)上,設(shè)計(jì)了遠(yuǎn)場(chǎng)子波品質(zhì)優(yōu)、低頻輸出能量強(qiáng)和穿透深度大的立體氣槍陣列延遲激發(fā)震源,并應(yīng)用到渤海、南黃海OBS深部地質(zhì)探測(cè),獲得了品質(zhì)良好的Ps,Pg,PmP,Pn 等多種震相,達(dá)到了預(yù)定目標(biāo);同時(shí)形成了針對(duì)淺水區(qū)OBS深部探測(cè)資料采集的立體氣槍陣列延遲激發(fā)震源設(shè)計(jì)技術(shù)方法,為淺水海域的OBS深部探測(cè)資料采集,在氣槍陣列震源設(shè)計(jì)方面起到了示范作用,也為深水海域OBS深部探測(cè)和油氣地震資料采集提供了借鑒.

    2立體氣槍陣列延遲激發(fā)震源激發(fā)特性

    2.1基本原理

    氣槍陣列震源(簡(jiǎn)稱槍陣,下同)沉放深度是影響其遠(yuǎn)場(chǎng)子波品質(zhì)的關(guān)鍵因素,理論模擬和實(shí)踐表明(何漢漪,2001;陳浩林等,2008;趙明輝等,2008; 李緒宣等,2009),為了獲得最佳的激發(fā)子波能量,需將槍陣沉放在略大于氣泡半徑的同一深度上,稱為平面(常規(guī))槍陣.此時(shí)激發(fā)的地震波在向下傳播的同時(shí)也向上傳播,上行波到達(dá)海面被反射后向下傳播,與正常下行波之間交互干涉形成了震源水體虛反射(簡(jiǎn)稱:震源鬼波),并引起了陷波作用,限制了地震資料的頻帶寬度.槍陣沉放深度小,氣槍外界水的壓力小,氣槍激發(fā)后高壓氣體將以更快的速度釋放,激發(fā)子波的頻率得到了提高,氣泡干擾得到壓制;但是,此時(shí)一部分能量變成了海水破碎能,從而大大降低了氣體振蕩的能量,受震源鬼波的陷波作用,低頻信號(hào)能量受到壓制,勘探深度受到了影響.槍陣沉放深度大,氣槍外界水的壓力變大,氣槍激發(fā)后高壓氣體釋放速度的變緩,激發(fā)子波的頻率降低,氣泡的震蕩加劇,從而大大降低了槍陣的波泡比,受震源鬼波陷波作用,低頻信號(hào)能量得到釋放,高頻信號(hào)能量受到壓制,雖然地震波的穿透深度得到了提高,但是地震波頻帶變窄、分辨率降低,原始資料的品質(zhì)下降.

    因此,平面槍列難以在勘探深度和原始資料品質(zhì)兩方面做到完全兼顧和統(tǒng)一,必須改進(jìn)槍陣的組合與激發(fā)方式,以獲得深部高品質(zhì)的地震信號(hào).

    借鑒陸地地震勘探中井中延遲激發(fā)技術(shù)壓制虛反射的經(jīng)驗(yàn)(Moldoveanu,2000;高銀波等,2006;譚紹泉,2003;趙殿棟等,2001;張會(huì)星等,2001),在海洋深水區(qū)地震勘探中,將平面槍陣改造成將陣列中的子陣列沉到放不同深度的多層震源,通過氣槍控制器按設(shè)定的時(shí)間延遲激發(fā),即從最上層子陣列開始依次激發(fā)各子陣列,使直達(dá)波同相疊加而虛反射非同相疊加,改善了地震子波的頻譜,提高了激發(fā)子波的信噪比(Cambois等, 2009;全海燕等,2011;趙仁永等,2011).Cambois等(Cambois等,2009)介紹了多層氣槍陣列震源的技術(shù)方法在海洋深水區(qū)油氣地震勘探中的應(yīng)用,指出其在提升低頻分量、拓寬頻帶和提高分辨率的作用;全海燕等(全海燕等,2011)通過對(duì)理論模擬和實(shí)際數(shù)據(jù)對(duì)比分析研究認(rèn)為,多層氣槍陣列震源能夠改善遠(yuǎn)場(chǎng)子波的頻譜品質(zhì),提高它的低頻段、高頻段能量和降低震源鬼波產(chǎn)生的陷波作用.

    2.2震源鬼波對(duì)氣槍子波的改造作用

    氣槍在水中激發(fā)時(shí),將氣槍中的高壓空氣瞬間釋放到水中,迅速形成球形的氣泡,由于氣泡內(nèi)壓力大于周圍水體壓力,導(dǎo)致氣泡迅速膨脹形成壓力脈沖,即氣槍的主脈沖,它向四周傳播形成地震波.同時(shí),部分向上傳播的信號(hào)經(jīng)海平面反射之后向下傳播并到達(dá)檢波點(diǎn),形成了震源鬼波.由于海水與空氣界面的反射系數(shù)為負(fù),因此震源鬼波與近場(chǎng)信號(hào)極性相反、到達(dá)接收點(diǎn)的時(shí)間滯后于下行地震波.存在到達(dá)時(shí)間延遲的震源鬼波與近場(chǎng)子波疊加,形成遠(yuǎn)場(chǎng)子波(圖1).

    圖1 氣槍遠(yuǎn)場(chǎng)子波形成示意圖Fig.1 Schematic diagram for the formation of air-gun far-field wavelet

    遠(yuǎn)場(chǎng)子波信號(hào)可表示為

    (1)

    式中,x(t)為近場(chǎng)子波信號(hào),Δt為震源鬼波延遲到達(dá)時(shí)間,x(t+Δt)為震源鬼波信號(hào),α為海水與空氣界面的反射系數(shù),約等于-1.0,則式(1)可改寫為

    (2)

    圖2為容量150in3的G槍在沉放6 m、壓力2000psi條件下的近場(chǎng)子波波形,它反映了氣槍激發(fā)后氣泡膨脹、壓縮、再膨脹、再壓縮循環(huán)往復(fù)的震蕩過程,由于振蕩過程中氣體分子摩擦消耗能量,所以每一次形成的脈沖振幅就越來越小,直至最后氣泡浮出水面破裂震蕩結(jié)束.

    圖2 容量為150in3的G槍近場(chǎng)子波Fig.2 The near field wavelet of G gun with volume 150in3

    圖3為該槍的遠(yuǎn)場(chǎng)子波波形,由震源鬼波造成的負(fù)極性子波脈沖延遲疊加,遠(yuǎn)場(chǎng)子波的正負(fù)脈沖震蕩并隨時(shí)間的增加近指數(shù)衰減的特征,與最小相位子波波形相近,說明震源鬼波改善了氣槍子波波形,使之更接近地震勘探的理論子波,遠(yuǎn)場(chǎng)子波的氣泡震蕩比近場(chǎng)子波小,表明鬼波有壓制氣泡作用的積極意義.

    圖3 容量為150in3的G槍沉放深度6 m時(shí)的遠(yuǎn)場(chǎng)子波Fig.3 The far field wavelet of G gun with volume 150in3

    對(duì)式(2)進(jìn)行傅氏變換得到震源的頻譜特征:

    (3)

    式中,f為頻率,(1-e-i2π·fΔt)可視為一次反射的濾波因子,即虛反射濾波器,設(shè)為H(f).假設(shè)震源與檢波器的連線與海面垂直,則

    (4)

    式中hg為氣槍陣列的沉放深度,c為海水速度.則濾波器的振幅譜為

    (5)

    由(5)式可得,當(dāng)頻率f=nc/2hg(n=1,2,3…)時(shí),濾波器的振幅譜受到震源鬼波的壓制最大,稱f=nc/2hg為陷波頻率(圖4),其大小由海水中聲波傳播速度c和槍陣沉放深度hg決定.

    圖4 容量為150in3的G槍沉放深度6 m時(shí)的 遠(yuǎn)場(chǎng)子波頻譜Fig.4 The far field wavelet spectrum of G gun with volume 150in3

    2.3立體氣槍陣列延遲激發(fā)震源頻率特性研究

    常規(guī)的氣槍陣列震源由包含著一定數(shù)量的單槍、相干槍的子陣列組成,所有的氣槍均沉放在針對(duì)勘探目標(biāo)所選擇的沉放深度下,氣槍的激發(fā)時(shí)間要求嚴(yán)格一致,由此它形成了一個(gè)疊合的下行波,也形成了一個(gè)疊合的上行波(圖5上),并且上行波與下行波具有相同的能量.因此,由上行波而產(chǎn)生的震源鬼波能量強(qiáng)、干擾大.

    圖5 常規(guī)槍陣(上)與順序激發(fā)的多層槍陣(下)產(chǎn)生下行波(實(shí)線)和 水體虛反射(虛線)的效應(yīng)圖(據(jù)Cambois等,2009)Fig.5 The downgoing wave(up) and water ghost(down) effect by conventional and sequence fired multi-levels air-gun array source (after Cambois et al.2009)

    Cambois等(Cambois等,2009)提出的多層氣槍陣列震源的工作原理如圖5(下)所示,將各子陣列沉放在不同的深度層上,從最上層子陣開始順序地延遲激發(fā)各層子陣列,延遲的時(shí)間是上層子陣列激發(fā)的下行波波前到達(dá)下一層子陣列的走時(shí),這樣在保證下行波波前同時(shí)疊加、能量不變的同時(shí),到達(dá)海平面的上行波能量不能同時(shí)疊加而受到削弱,降低了水體虛反射效應(yīng).與陸地炸藥井中延遲激發(fā)相比,激發(fā)介質(zhì)海水的物性相對(duì)均一,聲波速度基本恒定,而且子陣列沉放深度相對(duì)穩(wěn)定,其變化可以忽略不計(jì),精確的氣槍觸發(fā)控制完全可以做到下行波前同相疊加.該技術(shù)實(shí)現(xiàn)起來相對(duì)簡(jiǎn)單,只需要對(duì)現(xiàn)在的氣槍陣列的激發(fā)方式進(jìn)行小的改進(jìn).

    組合方式與容量相同的子陣列組成的多層氣槍陣列震源的遠(yuǎn)場(chǎng)子波信號(hào)可表示為

    y(t)=nx(t)-x(t+Δt1)-x(t+Δt2)…

    -x(t+Δtn),

    (6)

    式中,n為子陣列數(shù)量,x(t)為單子陣的近場(chǎng)子波信號(hào),x(t+Δtn)為第n個(gè)子陣列的虛反射信號(hào),Δt1、Δt2、Δtn為由式(4)求得的各子陣虛反射延遲到達(dá)時(shí)間,對(duì)式(6)做傅里葉變換可以得到多層槍陣的遠(yuǎn)場(chǎng)子波頻譜為

    Y(f)=(n-ei2πfΔt1-ei2πfΔt2…-ei2πfΔtn)F(f)

    =H(f)*F(f)

    (7)

    式中F(f)為子陣的近場(chǎng)子波頻譜,H(f)為多層槍陣鬼波濾波器,它的振幅譜為

    -4sin2(πfΔtn)]1/2

    (8)

    當(dāng)各子陣列的沉放深度相同時(shí),即為常規(guī)的平面氣槍陣列震源,則其鬼波濾波器的振幅為

    (9)

    對(duì)比式(8)和式(9)可以看出,多層氣槍陣列震源延遲激發(fā)方式,分散和降低了震源鬼波對(duì)遠(yuǎn)場(chǎng)子波的陷波作用,陷波頻率點(diǎn)分散且陷波作用大幅度降低,低頻和高頻段能量都得到了釋放;平面氣槍陣列震源陷波點(diǎn)集中、陷波作用強(qiáng).

    但是,多層氣槍陣列延遲激發(fā)震源在淺水區(qū)應(yīng)用還存在不能克服的困難,首先,氣槍的觸發(fā)同步精度一般在0~1.5 ms,為了發(fā)揮多層氣槍陣列延遲激發(fā)震源的優(yōu)勢(shì),需要將子陣列沉放深度差異加大到6 m(延遲激發(fā)時(shí)間4 ms)以上,這在深水采集環(huán)境中是不成問題的,但在淺水區(qū)由于水深條件的限制,幾乎不可能實(shí)現(xiàn)的;其次,沉放深度差異的加大,帶來了氣泡壓制作用的減弱和震源方位角的復(fù)雜化,造成了地震資料品質(zhì)的降低(Cambois等,2009).

    2.4立體氣槍陣列延遲激發(fā)震源組合設(shè)計(jì)方法及遠(yuǎn)場(chǎng)子波模擬分析

    基于以上分析,確定選用觸發(fā)同步精度高的中、小容量氣槍組成沉放深度差異≤3 m的立體氣槍陣列,在順序延遲激發(fā)時(shí)間設(shè)計(jì)上,不完全拘泥將上層子陣列激發(fā)地震波到達(dá)下一層子陣列的時(shí)間作為延遲激發(fā)時(shí)間,而是通過遠(yuǎn)場(chǎng)子波模擬計(jì)算和對(duì)其技術(shù)指標(biāo)的總體性能評(píng)價(jià)與優(yōu)選后,確定立體氣槍陣列組合與延遲激發(fā)的設(shè)計(jì)方案.

    通過對(duì)各種型號(hào)氣槍的性能對(duì)比分析,選擇可靠性強(qiáng)、觸發(fā)精度高、輸出能量大的G槍,作為淺海區(qū)OBS深部地質(zhì)探測(cè)的氣槍震源,在對(duì)不同的立體氣槍陣列組合與延遲激發(fā)時(shí)間設(shè)計(jì)方案進(jìn)行遠(yuǎn)場(chǎng)子波理論模擬分析的基礎(chǔ)上,針對(duì)勘探目標(biāo)需求,以拓展低頻、拓寬頻帶,提高遠(yuǎn)場(chǎng)子波特性和增加激發(fā)地震波穿透能力為目標(biāo),確定了由48條不同容量的G槍組成的4個(gè)子陣列、總?cè)萘繛?060in3的立體氣槍陣列組合,其中最大單槍容量為380in3,最小單槍容量40in3,按氣泡半徑計(jì)算公式(何漢漪,2001)算得最大氣泡半徑為5 m左右,滿足了淺海區(qū)OBS探測(cè)的震源激發(fā)的需求.

    在此基礎(chǔ)上進(jìn)行遠(yuǎn)場(chǎng)子波模擬,分析不同立體氣槍陣列延遲激發(fā)組合方案的遠(yuǎn)場(chǎng)子波特征.模擬環(huán)境參數(shù)為:氣槍工作壓力為2000psi,海水密度為1030 kg·m-3,海水溫度15 ℃,海水聲波速度為1500.0 m·s-1,海平面反射系數(shù)為-1.0,觀測(cè)點(diǎn)距離氣槍陣列震源中心9000 m.首先,通過改變子陣列的沉放深度組合和延遲激發(fā)時(shí)間,模擬相應(yīng)的遠(yuǎn)場(chǎng)子波,從中優(yōu)選出子波頻譜光滑、有效頻帶寬初泡比大的5組陣列組合(圖6);然后,進(jìn)一步分析遠(yuǎn)場(chǎng)子波波形和初峰值、峰峰值、陷波效應(yīng)等特征,優(yōu)選出最優(yōu)多層氣槍陣列組合.

    表1為不同立體氣槍陣列延遲激發(fā)震源與平面氣槍震源遠(yuǎn)場(chǎng)子波性能對(duì)比.平面氣槍陣列震源產(chǎn)生的續(xù)至波谷與海平面虛反射振幅疊合作用,造成其虛反射振幅高于初峰值,形成了強(qiáng)烈的陷波效應(yīng).因此它的頻帶較窄、低頻段能量相對(duì)較低、頻譜曲線鋸齒狀缺口較多,但海平面虛反射對(duì)壓制氣槍的氣泡效應(yīng)是有利的,因此它具有最高的波泡比;而立體氣槍陣列震源由于各子陣列的沉放深度不一,海平面虛反射分散疊加,虛反射振幅低于初峰值,其中四層氣槍陣列震源的虛反射最低、雙層氣槍陣列震源的虛反射次之;與平面氣槍陣列震源相比,立體氣槍陣列頻譜曲線相對(duì)光滑,低頻段能量較強(qiáng)、頻帶寬.

    圖7(左)為總?cè)萘繛?060 in3的S4-13-9-7-11-D1.5組合(A線)與其相同組合的總?cè)萘?000 in3(B線)氣槍陣列震源遠(yuǎn)場(chǎng)子波波形對(duì)比圖,這兩種槍陣均有較高的峰峰值和波泡比,虛反射效應(yīng)相對(duì)較低.其中,總?cè)萘?000 in3的組合中以4條1500 in3的Bolt氣槍替代較小容量的氣槍,由于在淺海區(qū)氣槍的沉放深度受到了限制,不能完全發(fā)揮大容量氣槍的能量?jī)?yōu)勢(shì),造成該槍陣的初峰值低于容量6060 in3組合的立體氣槍陣列震源.在頻譜對(duì)比結(jié)果(圖7右)上,容量6060 in3組合的低頻能量高于容量9000 in3的組合,頻譜曲線相對(duì)光滑,陷波效應(yīng)低,兼顧了獲得沉積層震相的要求.

    圖6 6組槍陣組合示意圖(上)及其遠(yuǎn)場(chǎng)子波(中)和頻譜(下)對(duì)比Fig.6 Schematic diagram for six air-gun array patterns (upper) and far field waveforms (middle) and spectrum (lower) simulation

    圖7 兩種氣槍陣列震源遠(yuǎn)場(chǎng)子波波形(左)和頻譜(右)對(duì)比圖Fig.7 The contrast of air-gun array far field wavelet (left) and spectrum (right)

    方式類型S4-8-8-8-8-D0.0S4-13-9-7-11-D1.5S4-14-10-8-12-D1.5S4-13-10-10-13-D2.0S4-11-8-8-11-D2.0S4-10-7-7-10-D2.0?組合方式平面、同時(shí)激發(fā)四層、1.5ms延遲激發(fā)四層、1.5ms延遲激發(fā)雙層、2ms延遲激發(fā)雙層、2ms延遲激發(fā)雙層、2ms延遲激發(fā)沉放模式8-8-8-813-9-7-1114-10-8-1213-10-10-1311-8-8-1110-7-7-10初峰值(bar·m)71.190.289.290.789.889.8虛反射強(qiáng)度(bar·m)77.646.445.854.153.453.4峰峰值(bar·m)148.7136.6135.0144.8143.2143.2波泡比16.815.313.012.014.414.4頻帶寬度(Hz)7~584.5~654.5~656~555~625~62低頻能量不突出強(qiáng)、突出強(qiáng)、突出較高較高較高頻譜曲線不光滑光滑光滑較光滑較粗糙較粗糙陷波頻率92Hz不明顯不明顯不明顯不明顯不明顯

    3應(yīng)用效果

    2013年8月,在國家自然科學(xué)基金“黃海及鄰區(qū)殼幔結(jié)構(gòu)及深淺構(gòu)造關(guān)系的綜合地球物理研究”(41210005)、國土資源部地質(zhì)調(diào)查專項(xiàng)和國家863計(jì)劃“海陸聯(lián)合深部地球物理探測(cè)關(guān)鍵技術(shù)研究”(2009AA093401)的支持下,青島海洋地質(zhì)研究所、中國科學(xué)院地質(zhì)與地球物理研究所和國家海洋局第一海洋研究所聯(lián)合開展了南黃?!綎|半島—渤海海陸聯(lián)合深部地震探測(cè).海區(qū)采用氣槍陣列震源密點(diǎn)激發(fā)(激發(fā)間距125 m)、大間距OBS接收(OBS臺(tái)站間距6 km),陸地地震臺(tái)站密點(diǎn)接收(臺(tái)站間距1 km),兩端各進(jìn)行一個(gè)大藥量放炮激發(fā).

    本次海上OBS深部探測(cè)首次采用立體氣槍陣列延遲激發(fā)震源技術(shù),考慮到渤海和南黃海海域的水深條件,對(duì)不同組合方式震源的遠(yuǎn)場(chǎng)子波模擬對(duì)比后,在南黃海選擇了容量6060 in3的S4-13-9-7-11-D1.4氣槍陣列組合方式;鑒于渤海海域探測(cè)區(qū)的水深不大于20 m,氣槍沉放后距海底的距離不能太小,以保證氣槍的工作安全,故采用了S4-11-8-8-11-D2.0組合方式.兩種組合方式的遠(yuǎn)場(chǎng)子波特征如圖6所示,性能對(duì)比如表1.

    圖8上為2011年在渤海采用6條1500 in3Bolt氣槍組成的容量9000 in3的平面氣槍陣列采集的OBS記錄,圖8下為2013年在相同位置,采用6060 in3的S4-11-8-8-11-D2.0組合采集的OBS記錄.可以看出,采用了新的震源激發(fā)技術(shù)后,采集的OBS數(shù)據(jù)質(zhì)量明顯得到提高,總體具有振幅能量強(qiáng)、震相豐富,視頻率相對(duì)較高,較為連續(xù)可靠的優(yōu)勢(shì)波組特征,OBS記錄信噪比高,氣槍信號(hào)有效傳播距離達(dá)150 km,大于平面氣槍陣列的傳播距離.

    圖8 平面槍陣(上)與立體槍陣延遲激發(fā)(下)OBS記錄剖面對(duì)比Fig.8 Comparison of OBS profile with conventional air-gun array (upper) and Tridimensional air-gun array (lower)

    南黃海海域OBS深部地震探測(cè)是首次進(jìn)行,通過對(duì)分布南黃海測(cè)線上的全部39個(gè)臺(tái)站記錄的震相特征和分布情況分析(表2)可以看出,本次探測(cè)獲得大量有價(jià)值的廣角反射/折射地震震相.其中包括:直達(dá)水波震相Pw,陸相沉積層折射震相Ps1,海相沉積層折射震相Ps2,地殼內(nèi)的折射波震相Pg、莫霍面的反射震相PmP和上地幔的折射震相Pn,填補(bǔ)了南黃海海域深部數(shù)據(jù)的空白,所得到的OBS記錄均獲得了延續(xù)長(zhǎng)度較大的震相,為OBS速度剖面反演奠定了良好的基礎(chǔ).

    圖9為位于千里巖隆起之上的OBS13臺(tái)站折合顯示剖面(折合速度6000 m·s-1,下同).緊隨直達(dá)水波震相之后是一組連續(xù)可追蹤、視速度高達(dá)6 m·s-1的折射震相Pb,這是該構(gòu)造區(qū)內(nèi)南黃海陸相盆地下的高速折射波震相,本文將其定義為Pb震相,它相當(dāng)于地殼內(nèi)的折射波震相Pg,因其發(fā)育在隆起區(qū)內(nèi),為了便于后期的走時(shí)擬合,故將其單獨(dú)定義.該震相在臺(tái)站左支2~60 km范圍內(nèi)可連續(xù)追蹤,其中在偏移距大約40 km處,其走時(shí)由0.4 s快速增加至0.6 s,結(jié)合侯方輝等(侯方輝等,2012)對(duì)該區(qū)域多道地震的解釋結(jié)果,推測(cè)該處應(yīng)為山東半島膠萊盆地在海區(qū)的延伸邊界,走時(shí)的突然增加代表該區(qū)域沉積層增厚,發(fā)育了沉積盆地;Pb震相在右半支可以追蹤到約35 km,剖面上顯示Pb震相在右半支明顯滯后于左半支,走時(shí)的大幅增加表明向SE方向已進(jìn)入沉積坳陷-北部坳陷;在右半支100~140 km處觀測(cè)到視速度高達(dá)8 m·s-1的PmP折射震相,由于臺(tái)站位于淺部地震波強(qiáng)反射屏蔽區(qū),接收的震相振幅相對(duì)較弱,信噪比相對(duì)較低.

    表2 南黃海OBS震相分布統(tǒng)計(jì)表

    圖10為位于南黃海盆地北部坳陷OBS19臺(tái)站,由于坳陷內(nèi)沉積物厚度較大,綜合地震剖面上可以清楚的觀察到Ps震相,左半支追蹤到偏移距約25 km處(圖11上),右半支追蹤到偏移距10 km處(圖11下),折合走時(shí)可達(dá)2~3 s,視速度約4.5 km·s-1;Pg震相在左半支25~55 km以及右半支10~25 km可以連續(xù)追蹤,起伏不大,說明該處地殼頂面較為平緩;在左半支的100~150 km處以及右半支的110~150 km處可以清楚觀察到莫霍面上的反射震相PmP,其視速度約8 km·s-1.

    圖9 13臺(tái)站OBS水聽器分量地震記錄剖面(折合速度為6 km·s-1)Fig.9 The hydrophone component seismic record profile at Station No.13 (reduced velocity is 6 km·s-1)

    圖10 OBS19站位初步處理剖面Fig.10 The preliminary processing seismic record profile at Station No.19 (reduced velocity is 6 km·s-1)

    圖11為位于中部隆起的OBS37站位初步處理剖面,站位的北支(圖11上)和南支(圖11下)均可以識(shí)別來自沉積層的廣角反射震相,并能在偏移距(Offset)65~100 km范圍內(nèi)識(shí)別到PmP震相;并由于站位位于中部隆起與北部坳陷的交界南側(cè),其北支偏移距范圍內(nèi)的沉積層底界深度變化較大,其所能見到的Pg震相延續(xù)范圍只有20 km左右;南支偏移距范圍內(nèi)的沉積層埋藏深度變化不大,Pg震相相對(duì)穩(wěn)定且延續(xù)距離在90 km以上.

    圖11 OBS37站位初步處理剖面Fig.11 The preliminary processing seismic record profile at Station No.37 (reduced velocity is 6 km·s-1)

    通過對(duì)南黃海OBS數(shù)據(jù)的速度結(jié)構(gòu)模擬等工作,獲得了2D速度結(jié)構(gòu)模型(祁江豪,2005),同時(shí)結(jié)合區(qū)域地質(zhì)、多道地震、重力及磁力資料進(jìn)行綜合分析、對(duì)比、驗(yàn)證后認(rèn)為,千里巖隆起上變質(zhì)巖推覆體推覆于海相沉積層之上,是蘇魯造山帶在海區(qū)的延伸標(biāo)志,變質(zhì)巖體速度約為6.3 km·s-1;北部坳陷自上而下分布陸相沉積層和海相沉積層,陸相沉積層速度約1.7~4 km·s-1,海相沉積層速度約5.2~6 km·s-1,厚度在3~8 km,再次驗(yàn)證了多道地震(吳志強(qiáng)等,2015)和重力、磁力反演推測(cè)結(jié)果(黃松等,2010).

    4結(jié)論及討論

    理論分析認(rèn)為,中、小容量的氣槍組成的立體氣槍陣列延遲激發(fā)震源,采用將氣槍子陣列沉放在不同的深度并由淺到深順序地激發(fā)的工作方法,能夠有效地衰減了由上行波反射引起的鬼波作用,拓寬了地震頻帶,低頻能量得到了提升,提高了震源穿透能力.這種工作方法可以有效避免大容量氣槍陣列在淺海區(qū)應(yīng)用的限制,既可以獲得深達(dá)莫霍面的反射/折射震相,也可以兼顧探測(cè)深部沉積層分布特征,有效地提高了OBS資料品質(zhì).

    在深水區(qū)海洋油氣地震勘探中,多層氣槍陣列延遲激發(fā)技術(shù)獲得了成功,有效地拓寬了地震頻帶,提高地震成像精度,為反演處理提供高品質(zhì)的基礎(chǔ)資料(Cambois等,2009;趙仁永等,2011).對(duì)多層氣槍陣列延遲激發(fā)技術(shù)進(jìn)行了改進(jìn),降低了子陣列沉放深度差異,采用了靈活的延遲激發(fā)時(shí)間,形成了立體氣槍陣列延遲激發(fā)震源設(shè)計(jì)技術(shù)方法,使之完全適應(yīng)淺水區(qū)苛刻的采集工作環(huán)境.該工作方法在OBS深部探測(cè)中還屬于首次應(yīng)用,且已有的海洋OBS深部探測(cè)絕大部分采用的是大容量氣槍組成的超大容量氣槍陣列震源(總?cè)萘恳话愠^6000 in3),以保證有較好的低頻性能和較大的穿透深度.本次中、小容量的氣槍組成的立體氣槍陣列延遲激發(fā)震源工作方法在我國淺海區(qū)(渤海、黃海)的應(yīng)用結(jié)果表明,該工作方法在淺海環(huán)境下也能激發(fā)出傳播距離長(zhǎng)、穿透深度深的地震波,并能有效壓制震源鬼波效應(yīng),可以得到高品質(zhì)的OBS數(shù)據(jù),既適用于巖石圈深部結(jié)構(gòu)、物質(zhì)組成和狀態(tài)的探測(cè)研究,也適用于中-深目標(biāo)層的多道地震探測(cè),為海洋地震勘探提供一種較理想的人工震源.

    立體氣槍陣列延遲激發(fā)的重要作用是壓制震源的虛反射效應(yīng).但是,震源遠(yuǎn)場(chǎng)子波的特征也會(huì)受到新組合方式的影響.首先,震源鬼波是氣泡脈沖的有效衰減器,立體氣槍震源會(huì)得到比常規(guī)震源低的波泡比(PTB);其次,由于個(gè)別子陣列沉放深度加大,在高靜水壓力下激發(fā),降低氣泡幅度、加大了震蕩周期,導(dǎo)致如圖6所示的部分高頻頻率能量的減少;第三,子陣列沉放深度差異造成氣泡幅度的差異,意味著氣泡脈沖衰減性能的降低.因此,與常規(guī)的氣槍陣列震源相比,遠(yuǎn)場(chǎng)子波模擬分析顯得更加重要,需要進(jìn)行細(xì)致的理論模擬分析,優(yōu)選高性能的氣槍陣列震源組合方式,盡量避免它所帶來的副作用.

    致謝感謝審稿專家和編輯部的大力支持.

    References

    Chen H L, Quan H Y, Yu G P, et al. 2008a. Summary of airgun source theory and technology(1). EGP, 18(4): 211-217.

    Chen H L, Quan H Y, Yu G P, et al. 2008b. Summary of airgun source theory and technology(2). EGP, 18(5): 300-308, 312. Gao Y B, Zhang Y. 2006. Forward modeling on delayed blast in the desert area of Junggar basin. GPP, 45(6): 642-645.

    Guillaume Cambois, Andrew Long, Gregg Parkes, et al. 2009. Multi-level airgun array: a simple and effective way to enhance the low frequency content of marine seismic data. SEG Houston 2009 International Exposition and Annual Meeting: 152-156.

    He H Y. 2001. High Resolution Marine Seismic Survey techniques and its application. China. Beijing: Geological Publishing House. Hou F H, Li R H, Zhang X H, et al. 2012. Jiaoliao basin extending into South Yellow Sea: evidence form a new seismic profile in South Yellow Sea. Marine Geology Frontiers (in Chinese), 28(3): 12-16.

    Li X X, Wen S L, Gu H M, et al. 2009. A numerical simulation of wavelets from offshore air-gun array seismic source. China Offshore Oil and Gas, 21(4): 215-220.

    Lin J M, Wang B S, Ge H K, et al. 2008. Study on large volume airgun source characteristics and seismic phase analysis. Chinese Journal of Geophysics, 51(1): 206-212. Luo G C, Wang B S, Ge H K, et al. 2006. Progress in earth′s deep structures exploration by airgun source. Progress in Geophysics (in Chinese), 21(2): 400-407. Lü C C, Hao T Y, Qiu X L, et al. 2011. Astudy on the deep structure of the northern part of southwest sub-basin from ocean bottom seismic data, South China Sea. Chinese J. Geophys. (in Chinese), 54(12): 3129-3138.

    McIntosh K, Nakamura Y, Wang T K, et al. 2005. Crustal-scale seismic profiles across Taiwan and the western Philippine Sea. Tectonophysics, 401: 23-54.

    Moldoveanu N. 2000. Vertical source array in marine seismic exploration.∥ 70th Annual International Meeting, SEG, Expanded Abstracts, 53-56. Nazareth J J, Clayton R W. 2003. Crustal structure of the Borderland-Continent Transition Zone of southern California adjacent to Los Angeles. JGR, 108(B8): 2404.

    Qi J H. 2015. Research on the Crustal Velocity Structure in the Southern Yellow Sea[D]. Beijing: China University of Geosciences. Qiu X L, Chen Y, Zhu R X, et al. 2007. The application of large volume airgun sources to the onshore-offshore seismic surveys: Implication from the experimental results in northern South China Sea. Chinese Science Bulletin (in Chinese), 52(4): 463-469.

    Qiu X L, Zhao M H, Xu H L, et al. 2012. Important processes of deep seismic surveys in the South China Sea: Retrospection and expectation. Journal of Tropical Oceanography, 31(3): 1-9.

    Quan H Y, Chen X H, Wei X B, et al. 2011. An analysis of the time-lapse explosion of airgun array in marine seismic. OGP, 46(4): 513-516.

    Tan S Q. 2003. The delay stack explosion technique and its applications. GPP, 42(4): 427-433.

    Tang J, Wang B S, Ge H K, et al. 2009. Study of experiment and simulation of large volume air-gun in deep structures exploration. Earthquake Research in China, 25(1): 1-10.

    Wu Z L, Yuan A G, Li J B, et al. 2008. New progress of deep crust sounding in the mid-northern South China Sea using ocean bottom seismometers. South China Journal of Seismology, 28(1): 21-28.

    Wu Z Q, Liu L H, Xiao G L, et al. 2015. Progress and enlightenment of integrated geophysics exploration of marine residual basin in the South Yellow Sea. Progress in Geophys., 30(5): 1692-1705.

    Yang H C, Gao J S. 2004. The air-gun source characteristics research for marine seismic exploration. Geophysical Prospecting for Petroleum, 43(4): 323-326. Zhang H X, Ning S N, Yang F, et al. 2001. The theoretical study of time-lapse explode methodand its application. Coal Geology of China, 13(4): 47-49.

    Zhao D D, Tan S Q, Zhang Q H, et al. 2001. Development and application of special source inseismic exploration. OGP, 36(4): 383-389.

    Zhao M H, Qiu X L, Ye C M, et al. 2004. Analysis on deep crustal structure along the onshore-offshore seismic profile across the Binhai(littoral) fault zone in northeastern South Chine Sea. 2004. Chinese J. Geophys. (in Chinese), 47(5): 845-852.

    Zhao M H, Qiu XL, Xia S H, et al. 2008. Large volume air-gun sources and its seismic waveform characters. Chinese J. Geophys. (in Chinese), 51(2): 558-565.

    Zhao R Y, Zhang Z B, XuanY H. 2011. Application of over/under streamer and over/under source seismic acquisition in the Pearl River Mouth Basin. OGP, 46(4): 517-521.

    Zhi P Y, Liu B H, Hua Q F, et al. 2012. Ocean bottom seismometer experiment and its preliminary results in Bohai Sea. Advances in Earth Science, 27(7): 769-777.

    附中文參考文獻(xiàn)

    陳浩林, 全海燕, 於國平等. 2008. 氣槍震源理論與技術(shù)綜述(上). 物探裝備, 18(4): 211-217.

    陳浩林, 全海燕, 於國平等. 2008. 氣槍震源理論與技術(shù)綜述(下). 物探裝備, 18(5): 300-308, 312.

    高銀波, 張研. 2006. 準(zhǔn)噶爾盆地沙漠區(qū)延遲爆炸激發(fā)正演模擬分析. 石油物探, 45(6): 642-645. 何漢漪. 2001. 海上高分辨率地震技術(shù)及其應(yīng)用. 北京;地質(zhì)出版社. 侯方輝, 李日輝, 張訓(xùn)華等. 2012. 膠萊盆地向南黃海延伸—來自南黃海地震剖面的新證據(jù). 海洋地質(zhì)前言, 28(3): 12-16.

    黃松, 郝天珧, 徐亞等. 2010. 南黃海殘留盆地宏觀分布特征研究. 地球物理學(xué)報(bào), 53(6): 1344-1353.

    李緒宣, 溫書亮, 顧漢明等. 2009. 海上氣槍陣列震源子波數(shù)值模擬研究. 中國海上油氣, 21(4): 215-220.

    羅桂純, 王寶善, 葛洪魁等. 2006. 氣槍震源在地球深部結(jié)構(gòu)探測(cè)中的應(yīng)用研究進(jìn)展. 地球物理學(xué)進(jìn)展, 21(2): 400-407.

    林建民、王寶善、葛洪魁等. 2008. 大容量氣槍震源特征及地震波傳播的震相分析. 地球物理學(xué)報(bào), 51(1): 206-212.

    呂川川, 郝天珧, 丘學(xué)林等. 2011, 南海西南次海盆北緣海底地震儀測(cè)線深部地殼結(jié)構(gòu)研究. 地球物理學(xué)報(bào), 54(12): 3129-3138.

    祁江豪. 2015. 南黃海地區(qū)地殼速度結(jié)構(gòu)研究[D]. 北京: 中國地質(zhì)大學(xué)(北京).

    丘學(xué)林, 陳龥, 朱日祥等. 2007. 大容量氣槍震源在海陸聯(lián)測(cè)中的應(yīng)用: 南海北部試驗(yàn)結(jié)果分析. 科學(xué)通報(bào), 52(4): 463-469.

    丘學(xué)林, 趙明輝, 徐輝龍等. 2012. 南海深地震探測(cè)的重要科學(xué)進(jìn)程: 回顧和展望. 熱帶海洋學(xué)報(bào), 31(3): 1-9.

    全海燕, 陳小宏, 韋秀波等. 2011. 氣槍陣列延遲激發(fā)技術(shù)探討. 石油地球物理勘探, 46(4): 513-516.

    吳振利, 阮愛國, 李家彪等. 2008. 南海中北部地殼深部結(jié)構(gòu)探測(cè)新進(jìn)展. 華南地震, 28(1): 21-28.

    吳志強(qiáng), 劉麗華, 肖國林等. 2015. 南黃海海相殘留盆地綜合地球物理調(diào)查進(jìn)展與啟示. 地球物理學(xué)進(jìn)展, 30(5): 1692-1705.

    譚紹泉. 2003. 震源延遲疊加技術(shù)及應(yīng)用效果. 石油物探, 42(4): 427-433.

    唐杰、王寶善、葛洪魁等. 2009. 大容量氣槍震源的實(shí)驗(yàn)與模擬研究. 中國地震, 25(1): 1-10.

    楊懷春, 高生軍. 2004. 海洋地震勘探中空氣槍震源激發(fā)特征研究. 石油物探, 43(4): 323-326.

    張會(huì)星, 寧書年, 楊峰等. 2001. 時(shí)延爆炸理論及應(yīng)用研究. 中國煤田地質(zhì), 13(4): 47-49.

    趙殿棟, 譚紹泉, 張慶淮等. 2001. 地震勘探中特殊震源的研制與應(yīng)用. 石油地球物理勘探, 36(4): 383-389.

    趙明輝, 丘學(xué)林, 葉春明等. 2004. 南海東北部海陸深地震聯(lián)測(cè)與濱海斷裂帶兩側(cè)地殼結(jié)構(gòu)分析. 地球物理學(xué)報(bào), 47(5): 845-852.

    趙明輝, 丘學(xué)林, 夏少紅等. 2008. 大容量氣槍震源及其波形特征. 地球物理學(xué)報(bào), 51(2): 558-565.

    趙仁永, 張振波, 軒義華. 2011. 上下源、上下纜地震采集技術(shù)在珠江口盆地的應(yīng)用. 石油地球物理勘探, 46(4): 517-521.

    支鵬遙, 劉保華, 華清峰等. 2012. 渤海海底地震儀探測(cè)試驗(yàn)及初步成果. 地球科學(xué)進(jìn)展, 27(7): 769-777.

    (本文編輯劉少華)

    基金項(xiàng)目國家自然科學(xué)基金項(xiàng)目(41210005,90814011),國家“863”計(jì)劃項(xiàng)目(2009AA093401),國土資源部地質(zhì)調(diào)查專項(xiàng)(GZH200800503,GZH201200509)聯(lián)合資助.

    作者簡(jiǎn)介吳志強(qiáng),男,1964年生,博士,研究員,主要從事海洋油氣地球物理研究工作.

    doi:10.6038/cjg20160722 中圖分類號(hào)P631, P738

    收稿日期2015-08-18,2016-03-09收修定稿

    Tridimensional air-gun array with delay fired source signal characteristics and the application in OBS exploration in shallow seas

    WU Zhi-Qiang1,2,3, HAO Tian-Yao4, TANG Song-Hua5, ZHANG Xun-Hua1,2,3, ZHENG Yan-Peng6,XIAO Guo-Lin1,2,3, ZHAO Wei-Na7,2,3

    1QingdaoInstituteofMarineGeology,MinistryofLandandResources,Qingdao266071,China2TheKeyLaboratoryofMarineHydrocarbonResourceandEnvironmentalGeology,MinistryofLandandResources,Qingdao266071,China3FunctionlaboratoryforMarineMineralResourceGeologyandExploration,QingdaoNationalOceanographyLaboratory,Qingdao266000,China4InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China5FirstMarineGeologicalParty,ShanghaiOffshorePetroleumBureau,Shanghai201200,China6FirstInstituteofOceanography,StateOceamicAdministrationPeople′sRepubicofChina,Qingdao266061,China7OceanUniversityofChina,Qingdao266063,China

    AbstractThe Bohai Sea is in the North China block and the South Yellow Sea is located in the Yangtze block. These seas are rich in oil and gas resources, related with block convergence and plate movement. The aim is to carry out seismic deep OBS detection in the Southern Yellow Sea and the Bohai Sea to get the seismic reflection and refraction information of lithosphere and the Moho surface for the study on the crust-mantle structure characteristics, the contact relation of the North China block and the Yangtze block and the tectonic effect of the plate movement and its control on the basin formation, evolution and hydrocarbon accumulation. The deep crust OBS exploration in the deep sea uses a conventional air-gun array which is made of several sub-arrays each containing a number of large volume air-guns. But in the shallow sea,it is difficult to play an effective role due to the gun depth limit and serious seafloor ghost reflection and source water ghost. To meet the needs for attenuation of water ghost reflection and enhancement of the low frequency content of seismic wavelet, based on the water depth condition and analysis of the principle of multi-level air-gun array and delay fired source applied in oil and gas exploration in the deep sea, we develop a method to delay the sub-array fired time for the tridimensional air-gun array source composed of a number of smaller volume air-guns. Compared with the multi-levels air-gun array,we apply a smaller depth difference and flexible sub-array sinking for sub air-gun array in order to adapt to the harsh environment for shooting the air-gun in the shallow water, thus we can get better effect for bubble suppression and ghost attenuation. A pattern of flexible delay fired time is designed for improving to output energy of the source and the quality of the wavelet. We design different combination schemes for tridimensional air-gun array and delay fired source and model its far-field wavelets. The amplitude-frequence distribution of our designed tridimensional air-gun array with delay fired source is more smooth than that of a conventional one, and more energy can be extended in the low frequency bands. Based on water depth conditions in the OBS exploration area,we applied respectively different tridimensional air-gun array and delay fired source pattern in the Bohai Sea and the Southern Yellow Sea. Compared with the deep crust OBS exploration using conventional large volume air-gun array in the Bohai Sea in 2012,the quality of the OBS data accquired by using new source is improved obviously, the seismic signal has a stronger amplitude, richer phase,higher resolution, better continuity and higher S/N ratio. The maximum propagation distance of source signal is over 150 km, which is more than that of the conventional large volume air-gun array. The deep crust OBS exploration which is the first one in the South Yellow Sea indicates the OBSs can record rich seismic phases including Ps,Pg,PmP and Pn both in the uplift areas where distributed the high velocity layers which had shielding effect on the seismic signals and the depression areas with large deposition thickness. The OBSs record quality is higher and the effective propagation distance of the air-gun signal is over 140 km.The deep crust OBS exploration test using tridimensional air-gun array with delay fired source technique shows that this method can generate seismic signals with long propagation distance and deep penetration in the shallow water, which can effectively suppress the source ghost effect and improve OBSs data quality. The method provides an ideal artificial source for marine seismic exploration, which can help study deep lithospheric structure, composition and status, and multichannel seismic detection for deep target layers. An important function for the tridimensional air-gun array and delay fired source technique is to attenuate the source ghost and to enhance low frequency energy. But the air-gun far field signals are also affected by the new design. Firstly, the source ghost is an effective attenuator of bubble pulses,the tridimensional air-gun array and delay fired source technique deliver a lower peak-to-bubble ratio (PTB) than a conventional air-gun array. Secondly, since the guns are towed deeper, they fire under higher hydrostatic pressure, which would reduce bubble sizes and oscillation periods and result in decrease of the high frequency contents. Thirdly, the difference of the bubble sizes caused by different array depths leads to less flexibility in array design and bubble pulse attenuation. Therefore, the far field wavelet modeling is the key to ensure an effective tridimensional air-gun array with delay fired source design. To avoid the side effects, we must carefully carry out the simulation and analysis of far-field wavelets to develop an optimal design for tridimensional air-gun array and delay fired source.

    KeywordsTridimensional air-gun array source; Delay fired source; Far-field wavelet simulation; OBS depth exploration; Source ghost suppression; Shallow sea

    吳志強(qiáng), 郝天珧, 唐松華等. 2016. 立體氣槍陣列延遲激發(fā)震源特性及在淺海區(qū)OBS探測(cè)中的應(yīng)用.地球物理學(xué)報(bào),59(7):2573-2586,doi:10.6038/cjg20160722.Wu Z Q, Hao T Y, Tang S H, et al. 2016. Tridimensional air-gun array with delay fired source signal characteristics and the application in OBS exploration in shallow seas. Chinese J. Geophys. (in Chinese),59(7):2573-2586,doi:10.6038/cjg20160722.

    超色免费av| 久久中文字幕人妻熟女| 日韩成人在线观看一区二区三区| 极品少妇高潮喷水抽搐| 色在线成人网| 一a级毛片在线观看| 亚洲精品国产精品久久久不卡| 18在线观看网站| 777米奇影视久久| 久久 成人 亚洲| 精品一区二区三区视频在线观看免费 | 日韩三级视频一区二区三区| 免费人成视频x8x8入口观看| 欧洲精品卡2卡3卡4卡5卡区| 麻豆成人av在线观看| 99re在线观看精品视频| 黄色视频,在线免费观看| 日韩有码中文字幕| 色播在线永久视频| 国产日韩欧美亚洲二区| 国产伦人伦偷精品视频| 成年女人毛片免费观看观看9 | 午夜福利,免费看| 1024视频免费在线观看| 人妻久久中文字幕网| 久久精品亚洲精品国产色婷小说| 国产成人欧美| 国产xxxxx性猛交| 在线观看一区二区三区激情| 国产不卡一卡二| 亚洲成av片中文字幕在线观看| 丝袜美腿诱惑在线| 亚洲五月天丁香| 91国产中文字幕| 制服诱惑二区| 日本vs欧美在线观看视频| 99久久精品国产亚洲精品| 国产熟女午夜一区二区三区| 成年人黄色毛片网站| 视频区图区小说| 国产成人av教育| 成年女人毛片免费观看观看9 | 亚洲色图av天堂| 黄色成人免费大全| 亚洲精品在线美女| 精品久久蜜臀av无| 91精品国产国语对白视频| 黄色 视频免费看| 欧美黄色片欧美黄色片| 99精品久久久久人妻精品| 精品久久久久久久久久免费视频 | 美女视频免费永久观看网站| 亚洲av日韩精品久久久久久密| 亚洲人成电影免费在线| 啦啦啦免费观看视频1| 天天影视国产精品| 国产99久久九九免费精品| 亚洲色图av天堂| 一级a爱片免费观看的视频| 人人妻人人爽人人添夜夜欢视频| 久久久久久久精品吃奶| x7x7x7水蜜桃| 在线av久久热| 国产精品永久免费网站| 国产成人精品久久二区二区免费| 水蜜桃什么品种好| 精品一区二区三卡| 亚洲九九香蕉| 在线观看午夜福利视频| 精品一区二区三区四区五区乱码| 精品一区二区三区四区五区乱码| 法律面前人人平等表现在哪些方面| 久久久久精品国产欧美久久久| 国产97色在线日韩免费| 久久久久国产一级毛片高清牌| 国产精品永久免费网站| 99精品在免费线老司机午夜| 久久 成人 亚洲| 久久午夜综合久久蜜桃| 日日夜夜操网爽| 色婷婷久久久亚洲欧美| 欧美日韩av久久| 亚洲精品中文字幕在线视频| 精品人妻1区二区| 精品乱码久久久久久99久播| 久久亚洲真实| av国产精品久久久久影院| ponron亚洲| 日韩熟女老妇一区二区性免费视频| 18禁黄网站禁片午夜丰满| 看黄色毛片网站| 午夜免费成人在线视频| 国产亚洲精品久久久久5区| 久久久久国产一级毛片高清牌| 久久久久国产精品人妻aⅴ院 | 欧美丝袜亚洲另类 | 精品国产乱子伦一区二区三区| 超碰97精品在线观看| 十八禁网站免费在线| 久久香蕉激情| 少妇猛男粗大的猛烈进出视频| 国产深夜福利视频在线观看| 亚洲一码二码三码区别大吗| 久久热在线av| 新久久久久国产一级毛片| 国产亚洲精品久久久久5区| 成人三级做爰电影| 热re99久久精品国产66热6| 国产精品二区激情视频| 国产成人免费观看mmmm| 国产精品98久久久久久宅男小说| 国产一区二区激情短视频| 国产国语露脸激情在线看| 看黄色毛片网站| 欧美av亚洲av综合av国产av| 亚洲一区二区三区欧美精品| 亚洲 欧美一区二区三区| 大型黄色视频在线免费观看| 在线看a的网站| 国产精品久久电影中文字幕 | 男女下面插进去视频免费观看| 午夜成年电影在线免费观看| 亚洲国产欧美日韩在线播放| 露出奶头的视频| 欧美丝袜亚洲另类 | 亚洲国产精品合色在线| 高清视频免费观看一区二区| 色婷婷av一区二区三区视频| 欧美激情久久久久久爽电影 | 一区福利在线观看| 后天国语完整版免费观看| 啦啦啦视频在线资源免费观看| 国产精品影院久久| 亚洲人成电影观看| 国产成人精品久久二区二区免费| 在线观看一区二区三区激情| 女性生殖器流出的白浆| 亚洲国产毛片av蜜桃av| 热99国产精品久久久久久7| 美女国产高潮福利片在线看| 亚洲成人免费电影在线观看| 国产成人欧美| 欧美 日韩 精品 国产| 午夜福利欧美成人| 男女床上黄色一级片免费看| 精品福利永久在线观看| 又大又爽又粗| 欧美另类亚洲清纯唯美| 欧美大码av| 亚洲成人免费电影在线观看| 美女扒开内裤让男人捅视频| 黑人巨大精品欧美一区二区mp4| 国产欧美日韩精品亚洲av| 色婷婷av一区二区三区视频| 激情视频va一区二区三区| 美女高潮到喷水免费观看| 色婷婷av一区二区三区视频| 精品国产亚洲在线| 欧美激情高清一区二区三区| 黄片播放在线免费| 丁香欧美五月| 1024视频免费在线观看| 精品久久久精品久久久| 高清视频免费观看一区二区| 一级片'在线观看视频| 成人特级黄色片久久久久久久| 亚洲中文av在线| 首页视频小说图片口味搜索| 国产99白浆流出| 99精国产麻豆久久婷婷| 国产精品一区二区在线不卡| 午夜精品国产一区二区电影| 国产精品 国内视频| 久久久久久久久免费视频了| 99国产极品粉嫩在线观看| 人成视频在线观看免费观看| 黑人欧美特级aaaaaa片| 99re在线观看精品视频| 如日韩欧美国产精品一区二区三区| 在线国产一区二区在线| 亚洲精品一二三| 女警被强在线播放| 久久香蕉精品热| 久久久国产精品麻豆| 老司机福利观看| 丰满饥渴人妻一区二区三| 亚洲av电影在线进入| 校园春色视频在线观看| av天堂在线播放| 色精品久久人妻99蜜桃| 久久精品亚洲精品国产色婷小说| 黄色视频不卡| 久久午夜综合久久蜜桃| 国产xxxxx性猛交| 国产av又大| 一本综合久久免费| 国产精品二区激情视频| 大码成人一级视频| 乱人伦中国视频| 亚洲精品自拍成人| av电影中文网址| 久久中文看片网| 在线观看日韩欧美| 黄色a级毛片大全视频| 啦啦啦在线免费观看视频4| 在线观看免费午夜福利视频| svipshipincom国产片| 99国产精品一区二区三区| 日韩欧美国产一区二区入口| 亚洲片人在线观看| 女人久久www免费人成看片| 欧美日韩一级在线毛片| 久久久久久久久免费视频了| 国产亚洲精品久久久久久毛片 | 国产精华一区二区三区| 国产aⅴ精品一区二区三区波| 欧美人与性动交α欧美软件| 天天躁夜夜躁狠狠躁躁| 80岁老熟妇乱子伦牲交| 性色av乱码一区二区三区2| 日韩有码中文字幕| 最近最新中文字幕大全电影3 | 一边摸一边抽搐一进一出视频| 自线自在国产av| 亚洲男人天堂网一区| 中国美女看黄片| 色尼玛亚洲综合影院| 国产乱人伦免费视频| 精品人妻在线不人妻| 亚洲aⅴ乱码一区二区在线播放 | 黄色片一级片一级黄色片| 久久久久精品国产欧美久久久| 午夜视频精品福利| 亚洲男人天堂网一区| 99久久国产精品久久久| 亚洲av成人一区二区三| 日韩免费av在线播放| 中文字幕另类日韩欧美亚洲嫩草| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久人人人人人| 狠狠婷婷综合久久久久久88av| 亚洲第一欧美日韩一区二区三区| 女同久久另类99精品国产91| 亚洲美女黄片视频| 国产男女超爽视频在线观看| 免费女性裸体啪啪无遮挡网站| 在线观看免费日韩欧美大片| 国产在线精品亚洲第一网站| 丁香六月欧美| 亚洲精品成人av观看孕妇| 韩国精品一区二区三区| 国产精品二区激情视频| 午夜福利在线观看吧| 亚洲成国产人片在线观看| 欧美 亚洲 国产 日韩一| 动漫黄色视频在线观看| 19禁男女啪啪无遮挡网站| 精品国产一区二区久久| 丁香六月欧美| 村上凉子中文字幕在线| 亚洲精品一二三| 亚洲精品粉嫩美女一区| 国产精品偷伦视频观看了| 欧美日韩视频精品一区| 国产国语露脸激情在线看| 亚洲成a人片在线一区二区| 久久精品亚洲精品国产色婷小说| 每晚都被弄得嗷嗷叫到高潮| 男人舔女人的私密视频| 国产精品九九99| 午夜福利在线观看吧| 99热国产这里只有精品6| 自拍欧美九色日韩亚洲蝌蚪91| 婷婷成人精品国产| 夜夜爽天天搞| 一级毛片精品| 午夜福利乱码中文字幕| cao死你这个sao货| 视频在线观看一区二区三区| 黑人欧美特级aaaaaa片| 日韩有码中文字幕| 欧美成人免费av一区二区三区 | 18禁观看日本| 亚洲精品美女久久久久99蜜臀| 亚洲av熟女| 欧美黄色淫秽网站| 亚洲午夜理论影院| 丝瓜视频免费看黄片| 欧美黑人欧美精品刺激| 亚洲黑人精品在线| av天堂久久9| 成年人午夜在线观看视频| 亚洲国产欧美日韩在线播放| 欧美激情极品国产一区二区三区| 国产成人精品无人区| 伦理电影免费视频| 亚洲成人手机| 一边摸一边抽搐一进一出视频| 欧美黄色片欧美黄色片| 麻豆av在线久日| svipshipincom国产片| 国产成人欧美| 99热只有精品国产| 一级作爱视频免费观看| 91麻豆av在线| 欧美精品亚洲一区二区| 精品人妻熟女毛片av久久网站| 免费人成视频x8x8入口观看| 首页视频小说图片口味搜索| 精品福利观看| 成人国产一区最新在线观看| 亚洲aⅴ乱码一区二区在线播放 | 视频在线观看一区二区三区| 精品少妇久久久久久888优播| 美女视频免费永久观看网站| 国产不卡一卡二| 国产熟女午夜一区二区三区| 国产精品久久电影中文字幕 | 在线永久观看黄色视频| 99国产精品一区二区蜜桃av | 久久精品国产清高在天天线| 少妇裸体淫交视频免费看高清 | 99精国产麻豆久久婷婷| 一区在线观看完整版| 法律面前人人平等表现在哪些方面| 国产精品一区二区精品视频观看| 亚洲男人天堂网一区| 日本一区二区免费在线视频| 91在线观看av| 在线视频色国产色| 精品人妻熟女毛片av久久网站| 免费在线观看黄色视频的| 成年人午夜在线观看视频| 女人被躁到高潮嗷嗷叫费观| 精品亚洲成国产av| 久久精品aⅴ一区二区三区四区| 香蕉国产在线看| 久久精品亚洲精品国产色婷小说| 国产蜜桃级精品一区二区三区 | 欧美国产精品va在线观看不卡| 十八禁网站免费在线| 欧美人与性动交α欧美精品济南到| 看片在线看免费视频| 欧美精品高潮呻吟av久久| 亚洲七黄色美女视频| 国产人伦9x9x在线观看| 亚洲熟女毛片儿| 日韩欧美一区视频在线观看| a级毛片在线看网站| 91大片在线观看| 午夜免费成人在线视频| 精品免费久久久久久久清纯 | 亚洲第一av免费看| 1024视频免费在线观看| 高潮久久久久久久久久久不卡| 午夜精品国产一区二区电影| 极品教师在线免费播放| 免费av中文字幕在线| 动漫黄色视频在线观看| 香蕉丝袜av| 91老司机精品| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 日韩 欧美 亚洲 中文字幕| 狂野欧美激情性xxxx| 免费黄频网站在线观看国产| 老司机福利观看| 国产一区在线观看成人免费| 国产成人啪精品午夜网站| 国产高清激情床上av| 法律面前人人平等表现在哪些方面| 在线观看免费视频网站a站| 国产激情久久老熟女| 国产真人三级小视频在线观看| 老司机午夜福利在线观看视频| 在线观看一区二区三区激情| 国产精品久久电影中文字幕 | 国产成人欧美在线观看 | 久久人人爽av亚洲精品天堂| 天堂动漫精品| 国产1区2区3区精品| 久久久国产精品麻豆| 王馨瑶露胸无遮挡在线观看| 少妇裸体淫交视频免费看高清 | 妹子高潮喷水视频| 国产欧美日韩一区二区精品| 日韩制服丝袜自拍偷拍| 在线av久久热| 国产高清激情床上av| 黑人猛操日本美女一级片| 人妻 亚洲 视频| 欧美国产精品一级二级三级| 免费在线观看完整版高清| 我的亚洲天堂| 成人免费观看视频高清| 一本综合久久免费| 国产精品久久久久久人妻精品电影| 777久久人妻少妇嫩草av网站| 亚洲精品自拍成人| 欧美亚洲日本最大视频资源| 欧美日韩一级在线毛片| 国产男靠女视频免费网站| 性少妇av在线| 午夜影院日韩av| 欧美一级毛片孕妇| 国产精品国产av在线观看| 丰满迷人的少妇在线观看| 男女床上黄色一级片免费看| 在线观看免费高清a一片| 露出奶头的视频| 成人永久免费在线观看视频| 满18在线观看网站| 亚洲欧美日韩另类电影网站| 两性夫妻黄色片| 老熟妇仑乱视频hdxx| 视频在线观看一区二区三区| a级毛片黄视频| 亚洲黑人精品在线| 中文亚洲av片在线观看爽 | 乱人伦中国视频| 国产成人一区二区三区免费视频网站| av片东京热男人的天堂| 国产xxxxx性猛交| 国产精品 欧美亚洲| 他把我摸到了高潮在线观看| 丰满迷人的少妇在线观看| 日本黄色日本黄色录像| 色综合婷婷激情| 欧美日韩乱码在线| 国产精品98久久久久久宅男小说| 欧美成狂野欧美在线观看| 12—13女人毛片做爰片一| 精品国产一区二区久久| 国产又爽黄色视频| 制服人妻中文乱码| 99在线人妻在线中文字幕 | 欧美在线一区亚洲| 丰满饥渴人妻一区二区三| 在线观看日韩欧美| 免费黄频网站在线观看国产| 国产精品欧美亚洲77777| 91麻豆精品激情在线观看国产 | 飞空精品影院首页| 国产区一区二久久| 80岁老熟妇乱子伦牲交| 脱女人内裤的视频| 中文字幕人妻丝袜制服| 免费观看人在逋| 王馨瑶露胸无遮挡在线观看| 99国产精品一区二区蜜桃av | 热99re8久久精品国产| 午夜福利视频在线观看免费| 中出人妻视频一区二区| 他把我摸到了高潮在线观看| 夜夜躁狠狠躁天天躁| 亚洲av日韩精品久久久久久密| 国产亚洲精品一区二区www | 午夜精品国产一区二区电影| 国产成人欧美在线观看 | 一级毛片精品| 男女之事视频高清在线观看| 久热爱精品视频在线9| 午夜精品在线福利| 涩涩av久久男人的天堂| 国内久久婷婷六月综合欲色啪| 人妻一区二区av| 色综合婷婷激情| 亚洲av日韩在线播放| 女同久久另类99精品国产91| 一进一出抽搐动态| 一级a爱片免费观看的视频| 国产片内射在线| 国产欧美亚洲国产| 不卡av一区二区三区| 又大又爽又粗| 嫩草影视91久久| 老司机福利观看| 日韩欧美三级三区| 99久久99久久久精品蜜桃| 国产成人精品无人区| 国产亚洲av高清不卡| 精品熟女少妇八av免费久了| 日本欧美视频一区| 两个人看的免费小视频| 成在线人永久免费视频| 老司机在亚洲福利影院| 国产欧美日韩一区二区精品| 国产一区有黄有色的免费视频| 欧美精品av麻豆av| 亚洲中文av在线| 老司机亚洲免费影院| 大香蕉久久成人网| 新久久久久国产一级毛片| 美女扒开内裤让男人捅视频| 亚洲熟妇熟女久久| 欧美精品啪啪一区二区三区| 久久人妻福利社区极品人妻图片| 黑人操中国人逼视频| 亚洲性夜色夜夜综合| 午夜91福利影院| 一个人免费在线观看的高清视频| 啦啦啦视频在线资源免费观看| 亚洲精品国产精品久久久不卡| 午夜免费鲁丝| 欧美黑人欧美精品刺激| 伦理电影免费视频| 久热这里只有精品99| 亚洲综合色网址| 99久久人妻综合| 很黄的视频免费| 日韩欧美国产一区二区入口| 国产精品秋霞免费鲁丝片| 嫩草影视91久久| 国产成人一区二区三区免费视频网站| 大型黄色视频在线免费观看| 日本一区二区免费在线视频| 成人国产一区最新在线观看| 久久久精品国产亚洲av高清涩受| 好男人电影高清在线观看| 国产国语露脸激情在线看| 嫩草影视91久久| x7x7x7水蜜桃| 极品少妇高潮喷水抽搐| 性少妇av在线| 午夜福利视频在线观看免费| 久久精品国产99精品国产亚洲性色 | 午夜影院日韩av| 国产高清国产精品国产三级| 亚洲中文av在线| 国产乱人伦免费视频| 免费在线观看黄色视频的| av在线播放免费不卡| 免费在线观看日本一区| 在线看a的网站| 少妇粗大呻吟视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美日韩精品亚洲av| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久视频播放| 午夜两性在线视频| 国产精品自产拍在线观看55亚洲 | 欧美乱色亚洲激情| 俄罗斯特黄特色一大片| 欧美精品啪啪一区二区三区| 天堂中文最新版在线下载| 精品福利观看| 欧美成狂野欧美在线观看| 男人舔女人的私密视频| 久久久久久亚洲精品国产蜜桃av| 巨乳人妻的诱惑在线观看| 欧美另类亚洲清纯唯美| 宅男免费午夜| 中文字幕av电影在线播放| 成在线人永久免费视频| 高清欧美精品videossex| 久久久久久久久久久久大奶| 少妇 在线观看| 久久九九热精品免费| 999久久久精品免费观看国产| 国产欧美日韩一区二区精品| 乱人伦中国视频| 校园春色视频在线观看| 男女午夜视频在线观看| 国产亚洲欧美精品永久| 18禁国产床啪视频网站| 人妻丰满熟妇av一区二区三区 | 中出人妻视频一区二区| 日韩欧美一区二区三区在线观看 | 成人国产一区最新在线观看| 成人国语在线视频| 国产欧美日韩一区二区三区在线| 亚洲欧美激情综合另类| 国产xxxxx性猛交| 久久ye,这里只有精品| 美女高潮喷水抽搐中文字幕| 看片在线看免费视频| videosex国产| 激情在线观看视频在线高清 | 亚洲午夜精品一区,二区,三区| 免费观看人在逋| 成年人午夜在线观看视频| 大型av网站在线播放| 亚洲精品一卡2卡三卡4卡5卡| 欧美亚洲 丝袜 人妻 在线| 久久久久久免费高清国产稀缺| 两性夫妻黄色片| 三上悠亚av全集在线观看| 精品久久久久久电影网| 国产成人系列免费观看| 免费av中文字幕在线| 男女床上黄色一级片免费看| 人人妻人人添人人爽欧美一区卜| 成人亚洲精品一区在线观看| 欧美性长视频在线观看| 麻豆av在线久日| 12—13女人毛片做爰片一| 99热网站在线观看| 天天影视国产精品| 国产人伦9x9x在线观看| 国产一区二区三区综合在线观看| 亚洲av成人一区二区三| 亚洲一码二码三码区别大吗| 又紧又爽又黄一区二区| 久久人人爽av亚洲精品天堂| 一本一本久久a久久精品综合妖精| 脱女人内裤的视频| 国产精品一区二区在线观看99| 不卡av一区二区三区| 夜夜爽天天搞| 精品亚洲成a人片在线观看| 脱女人内裤的视频| 久久久久国内视频| 免费黄频网站在线观看国产| 一级毛片精品| 黄片大片在线免费观看|