復(fù)雜障礙物區(qū)三維地震觀測系統(tǒng)多方位角變觀設(shè)計(jì)技術(shù)及應(yīng)用

田錦瑞, 亞東菊, 張 昭, 曹秀森, 張曉盼, 牛清華, 邱兆泰

(河北省煤田地質(zhì)局 物測地質(zhì)隊(duì),邢臺 054000)

0 引言

三維地震勘探項(xiàng)目的地震數(shù)據(jù)采集難度,通常與區(qū)內(nèi)的地表復(fù)雜程度存在較大關(guān)系。當(dāng)施工區(qū)域內(nèi)存在湖泊、村莊、工廠等大型障礙區(qū)時,障礙區(qū)下一定深度的地震數(shù)據(jù)將會丟失。這致使區(qū)域中淺層地質(zhì)任務(wù)無法完成,并且地震數(shù)據(jù)采集的施工成本控制將變得極其困難。

目前解決這些問題的主要技術(shù)手段為,增加布置與原有觀測系統(tǒng)同測線方位角的特殊觀測系統(tǒng)。特殊觀測系統(tǒng)設(shè)計(jì)通常首先要依據(jù)障礙區(qū)復(fù)雜程度進(jìn)行分級,在此基礎(chǔ)上將局部的激發(fā)點(diǎn)和接收點(diǎn)的物理位置及組合方式進(jìn)行合理、可行的修改,以彌補(bǔ)因障礙區(qū)造成的目的層地震數(shù)據(jù)缺失。優(yōu)秀的特殊觀測系統(tǒng)設(shè)計(jì)需在覆蓋次數(shù)、炮檢距分布及方位角分布上具備良好的均勻性。

近些年,諸多學(xué)者在特殊觀測系統(tǒng)方面做了許多研究,徐忠華[1]設(shè)計(jì)了Design地震軟件,該軟件具有人工特殊觀測系統(tǒng)設(shè)計(jì)功能,可對實(shí)際覆蓋情況進(jìn)行模擬演示,并及時指導(dǎo)野外生產(chǎn);李飛等[2]提出樹枝型特殊觀測系統(tǒng),在野外現(xiàn)場進(jìn)行實(shí)施中取得了較好的應(yīng)用效果;孫希杰[3]針對大村莊、鐵路及公路等多種障礙物,使用Goody GIS結(jié)合克浪軟件,采用人機(jī)交互方式反復(fù)論證,取得高質(zhì)量的地震原始數(shù)據(jù);劉軍勝[4]利用加減法設(shè)計(jì)思路,通過Omni軟件的觀測系統(tǒng)設(shè)計(jì)模塊,在激發(fā)點(diǎn)障礙區(qū)的特殊觀測系統(tǒng)設(shè)計(jì)應(yīng)用中獲取了高品質(zhì)的地震資料;王樹威[5]應(yīng)用克浪地震采集設(shè)計(jì)系統(tǒng)進(jìn)行特殊觀測系統(tǒng)設(shè)計(jì),在我國東部某復(fù)雜區(qū)的實(shí)際應(yīng)用中取得了良好的野外數(shù)據(jù)資料。

基于此,筆者使用了基于貪心算法的地震數(shù)據(jù),采集特殊觀測系統(tǒng)設(shè)計(jì)優(yōu)化程序[6],結(jié)合在山東省西南部地表復(fù)雜區(qū)的三維地震數(shù)據(jù)采集工作,針對該區(qū)中應(yīng)用不同測線方位角的觀測系統(tǒng)疊加,實(shí)現(xiàn)特殊觀測系統(tǒng)設(shè)計(jì)進(jìn)行了探討,為不同測線方位角的觀測系統(tǒng)疊加,在大型障礙區(qū)的特殊觀測系統(tǒng)設(shè)計(jì)做出了嘗試。

1 系統(tǒng)設(shè)計(jì)優(yōu)化程序簡介

地震數(shù)據(jù)采集特殊觀測系統(tǒng)設(shè)計(jì)優(yōu)化程序是基于貪心算法開展的程序設(shè)計(jì)。

貪心算法是一種以局部最優(yōu)解疊加來近似達(dá)到整體最優(yōu)解的算法。在程序的設(shè)計(jì)過程中,結(jié)合當(dāng)前CMP面元中疊加次數(shù)的缺失情況,選取修補(bǔ)效果最佳的炮檢組合,在不斷修改缺失情況的同時,不斷累加這種最佳組合以達(dá)到最為經(jīng)濟(jì)的變觀方案。

算法的實(shí)現(xiàn)如圖1所示,由CMP面元中缺失的疊加次數(shù)建立模型,并評價所有可用炮檢組合的對于缺失模型的修補(bǔ)程度,選取修補(bǔ)效果最佳的一個組合,將該組合修補(bǔ)后的模型替代原模型,并將該組合添加至最終輸出的方案中,繼續(xù)以上循環(huán)。在該循環(huán)過程中,若選取的最佳組合滿足跳出條件,即選取的組合數(shù)達(dá)到設(shè)定的上限,則修補(bǔ)效果滿足需求,則輸出最終方案,并終止程序。圖1中虛框內(nèi)的流程即為貪心算法在本程序設(shè)計(jì)中的應(yīng)用過程。

圖1 地震數(shù)據(jù)采集變觀設(shè)計(jì)優(yōu)化程序流程圖Fig.1 Flow chart of the variable-view design optimization program for seismic data acquisition

2 方位角觀測系統(tǒng)疊加的理論支撐

當(dāng)前主流的地震數(shù)據(jù)采集設(shè)計(jì)中,整個地震數(shù)據(jù)采集過程通常使用單一測線方位角的觀測系統(tǒng)。單一測線方位角的三維地震數(shù)據(jù)采集系統(tǒng),可以將地震波反射點(diǎn)集中于設(shè)計(jì)的面元中心,對于后期數(shù)據(jù)的準(zhǔn)確性有著重要意義。單一測線方位角的地震數(shù)據(jù)采集系統(tǒng)可以更加高效地施工,使施工成本得到良好控制。

然而,在某些工區(qū)施工過程中,因地震數(shù)據(jù)采集區(qū)域內(nèi)地表?xiàng)l件極為復(fù)雜,大型障礙物(如村莊、河湖)的存在,致使設(shè)計(jì)的接收點(diǎn)與激發(fā)點(diǎn)無法布設(shè),影響障礙物下方的地震數(shù)據(jù)采集,且在原設(shè)計(jì)的測線方位角下進(jìn)行特殊觀測系統(tǒng)設(shè)計(jì)已無法達(dá)成項(xiàng)目任務(wù)。

因此,在實(shí)際工作中我們需要結(jié)合村莊內(nèi)部建筑、道路的實(shí)際特點(diǎn),考慮布設(shè)不同與原設(shè)計(jì)測線方位角的觀測系統(tǒng),沿村莊道路進(jìn)行布設(shè)接收點(diǎn),并在村莊內(nèi)主要道路進(jìn)行采集作業(yè),完成障礙物下方的地震數(shù)據(jù)采集任務(wù)。

實(shí)際施工過程中,由于觀測系統(tǒng)間的測線方位角的不同、激發(fā)點(diǎn)和接收點(diǎn)不符合理論點(diǎn)位,會造成觀測系統(tǒng)面元中的CMP(共中心點(diǎn))位置分散,對于這種情況,數(shù)據(jù)采集質(zhì)量將與面元設(shè)計(jì)大小具有直接關(guān)系。

在地震數(shù)據(jù)采集觀測系統(tǒng)的設(shè)計(jì)過程中,面元的大小是非常重要的[7-10]。面元的大小將直接影響地震數(shù)據(jù)資料的橫向分辨率精度,更對構(gòu)造及斷裂細(xì)節(jié)特征的解釋產(chǎn)生重要的影響,對于各面元疊加時的反射信息的真實(shí)性具有根本性作用[11]。因此,特觀設(shè)計(jì)必須滿足面元大小的要求。鑒于對以上因素的考慮,面元大小必須滿足以下三個方面的要求:

1)滿足橫向分辨率的要求。考慮分辨主要目的層的最小地質(zhì)體的橫向尺度,則需要在地震信號優(yōu)勢頻率的波長內(nèi)要有不少于3個采樣點(diǎn),對應(yīng)的面元就能保證有良好的分辨率,具體公式為式(1)。

b=Vint/(2×fp)

(1)

式中:b為面元尺寸(m);Vint為目的層上覆地層層速度(m/s);fp為目的層主頻(Hz)。

2)滿足最高無混疊頻率的要求?？紤]偏移前最高無混疊頻率的要求,每個傾斜同相軸都有一個偏移前可能的最高無混疊頻率fmax,它依賴于此同相軸的上一層的速度Vint、傾角θ和CMP面元的大小,計(jì)算公式為式(2)。

b=Vint/(4×fmax×sinθ)

(2)

式中:b為面元尺寸(m);Vint為目的層上覆地層層速度(m/s);fmax為最高無混疊頻率(Hz),θ為地層最大傾角(°)。

3)滿足30°繞射收斂要求。

b=Vrms/(4×fmax×sin 30° )

(3)

式中:b為面元尺寸(m);Vrms為均方根速度(m/s);f為最高無混疊頻率(Hz)。

鑒于以上參數(shù),結(jié)合目標(biāo)區(qū)域的相關(guān)參數(shù)得出數(shù)據(jù)如表1所示。

表1 目標(biāo)區(qū)域的面元尺寸分析計(jì)算表Tab. 1 Analysis and calculation of bin size in target area

該區(qū)實(shí)際使用的面元采集參數(shù)為5 m×10 m,依據(jù)論證結(jié)果,在目標(biāo)區(qū)中,當(dāng)選擇面元最小尺寸不大于15.4 m時,則面元內(nèi)數(shù)據(jù)可滿足反射信息正確成像;當(dāng)選擇面元最小尺寸不大于10.26 m時,則可滿足對斷點(diǎn)繞射的充分收斂;當(dāng)選擇面元最小尺寸的不大于60.4 m時,則可滿足最高無混疊頻率的要求[12-13]。

由此可知,在工區(qū)的實(shí)際采集過程中,雖然測線間方向存在夾角且激發(fā)、接收點(diǎn)布設(shè)不符合理論點(diǎn)位,導(dǎo)致面元中CMP位置分散,但其仍在設(shè)計(jì)面元大小范圍內(nèi),對于采集數(shù)據(jù)應(yīng)無顯著影響,即采用不同測線方位角觀測系統(tǒng)疊加采集,不會顯著影響數(shù)據(jù)采集的質(zhì)量。

3 設(shè)計(jì)應(yīng)用范例

為檢驗(yàn)不同測線方位角觀測系統(tǒng)疊加設(shè)計(jì)方案,選擇鮑店煤礦八采區(qū)作為本次應(yīng)用的研究靶區(qū)。

鮑店煤礦行政區(qū)劃歸上隸屬鄒城市管轄,八采區(qū)位于井田西南部,區(qū)內(nèi)有北林廠村、鮑廠村兩個村莊,還分布有在建工地、選煤廠、飼料工廠及廢棄工廠等零星建筑物。

本目標(biāo)區(qū)三維地震勘探觀測系統(tǒng)設(shè)計(jì)為14L×4S×48T×7R×24次線束狀,線距為20 m,道距為10 m,縱向最大炮檢距為240 m,橫向最大炮檢距為190 m,最大非縱炮檢距為306.1 m。

北林廠村對于原設(shè)計(jì)施工的影響最為嚴(yán)重,該村建筑房屋間距較小,村內(nèi)道路呈“井”字形,道路方向與設(shè)計(jì)測線方向夾角5°。村莊內(nèi)設(shè)計(jì)接收點(diǎn)、激發(fā)點(diǎn)均無法正常布設(shè)(圖2)。

對北林廠村整個障礙區(qū)使用最大炮檢距為250 m。將村內(nèi)不能布設(shè)的接收點(diǎn)、激發(fā)點(diǎn)去除后,疊加情況如圖3所示。按照原設(shè)計(jì)的測線方位角建立特殊接收排列并布設(shè)激發(fā)點(diǎn)的傳統(tǒng)特殊觀測系統(tǒng)設(shè)計(jì)方式,經(jīng)過多次嘗試,無法滿足疊加次數(shù)需求。

圖3 北林廠村疊加情況圖Fig.3 Superposition of Beilinchang Village

在實(shí)際采集過程中結(jié)合實(shí)際情況,沿村莊內(nèi)道路布設(shè)接收點(diǎn),依據(jù)村莊內(nèi)道路特點(diǎn)建立不同于原觀測系統(tǒng)測線方位角的觀測系統(tǒng),彌補(bǔ)原設(shè)計(jì)中因激發(fā)點(diǎn)和接收點(diǎn)缺失所造成的疊加次數(shù)減少情況。

將沿村莊道路布設(shè)的接收點(diǎn)定義為一個完整接收排列,應(yīng)用地震數(shù)據(jù)采集特殊觀測系統(tǒng)設(shè)計(jì)優(yōu)化程序進(jìn)行特殊觀測系統(tǒng)設(shè)計(jì)。將程序生成的接收組合與實(shí)際情況相結(jié)合進(jìn)行地震數(shù)據(jù)采集工作,實(shí)際施工中激發(fā)點(diǎn)位、接收點(diǎn)位與實(shí)際疊加覆蓋次數(shù)如圖4、圖5所示。

圖4 北林廠村實(shí)際布設(shè)情況圖Fig.4 Actual layout of Beilinchang Village

圖5 北林廠村實(shí)際疊加情況圖Fig.5 Actual superposition of Beilinchang Village

圖4為該區(qū)域地震數(shù)據(jù)采集所采用的激發(fā)點(diǎn)、接收點(diǎn)的布設(shè)情況。該區(qū)域的特殊觀測系統(tǒng)設(shè)計(jì)所選用的最大炮檢距為250 m,相較于原設(shè)計(jì)增加激發(fā)點(diǎn)140個。

實(shí)際施工中村(障礙物)內(nèi)采用了上面的特殊觀測系統(tǒng),村(障礙物)外采用了原來的觀測系統(tǒng)。圖5為村內(nèi)村外兩種觀測系統(tǒng)疊加后的實(shí)際施工圖,從圖5可以看出,村莊障礙區(qū)域內(nèi)疊加次數(shù)較為均勻,且在北林廠村為主體的障礙區(qū)域內(nèi)面元疊加次數(shù)均已達(dá)到10次左右,無零次疊加區(qū)域,保證了數(shù)據(jù)采集質(zhì)量,達(dá)到了地震數(shù)據(jù)采集的要求。

該項(xiàng)目采集完成后,在地震數(shù)據(jù)處理軟件中重建觀測系統(tǒng),CMP(共中心點(diǎn))與面元的相對關(guān)系見圖6。

圖6 北林廠村區(qū)域內(nèi)CMP與面元的相對關(guān)系圖Fig.6 Actual superposition of Beilinchang Village

圖6(a)面元位置位于北林廠村中,該面元內(nèi)僅有沿村莊內(nèi)道路布設(shè)的特殊方位角觀測系統(tǒng)所產(chǎn)生的CMP點(diǎn)。由于村內(nèi)激發(fā)點(diǎn)與接收點(diǎn)位無法規(guī)則布置,致使CMP位置較為分散。圖6(b)的面元位置位于北林廠村邊,該處面元中同時存在原觀測系統(tǒng)設(shè)計(jì)與特殊方位角觀測系統(tǒng)設(shè)計(jì)的CMP。由于村內(nèi)接收點(diǎn)位無法規(guī)則布置,致使特殊方位角觀測系統(tǒng)產(chǎn)生的CMP位置較為分散。村莊外布設(shè)的原觀測系統(tǒng)由于接收點(diǎn)的布設(shè)相對規(guī)律,所產(chǎn)生的CMP位置較為集中。圖6(c)的面元位置位于北林廠村外,該處面元中的CMP僅為原觀測系統(tǒng)設(shè)計(jì)產(chǎn)生的CMP,由于接收點(diǎn)位布置較規(guī)則,面元中CMP位置較為集中。

對該區(qū)域的地震數(shù)據(jù)進(jìn)行完整處理后,在北林廠村區(qū)域中,分別在村莊中部沿東西、南北兩個方向選取剖面見圖7。

圖7 北林廠村區(qū)域地震剖面圖Fig.7 Regional seismic profile of Beilinchang Village

圖7中3號煤層目的層在地震時間剖面上的反映較為清晰連續(xù),基本完成了該項(xiàng)目的地質(zhì)任務(wù)。其次,特觀前在村莊范圍內(nèi)淺層地震數(shù)據(jù)連續(xù)性較差存在開天窗問題,經(jīng)方位角不同個觀測系統(tǒng)疊加后,村莊內(nèi)外雖有差距,但反射波連續(xù)性較好。

總體來看,不同于原設(shè)計(jì)測線方位角的特觀設(shè)計(jì)與原設(shè)計(jì)的觀測系統(tǒng)交疊應(yīng)用,解決了實(shí)際地震數(shù)據(jù)采集中因障礙區(qū)過大而無法采集的問題,雖然由于障礙區(qū)域內(nèi)接收點(diǎn)與激發(fā)點(diǎn)無法規(guī)則布設(shè),導(dǎo)致特殊方位角觀測系統(tǒng)的CMP位置無法集中于設(shè)計(jì)面元的中心,但面元大小的設(shè)計(jì)有效控制了該情況對于數(shù)據(jù)的影響程度,對于采集數(shù)據(jù)無顯著影響。就應(yīng)用效果而言,雖然在最終處理的時間剖面中特殊方位角觀測系統(tǒng)所參與的剖面質(zhì)量,較原觀測系統(tǒng)設(shè)計(jì)下的剖面仍具有一定差距,但基本可以解決障礙區(qū)域內(nèi)的地質(zhì)任務(wù),尤其在淺層地震剖面的質(zhì)量上,影響較為顯著。

4 結(jié)論

基于貪心算法設(shè)計(jì)的特殊觀測系統(tǒng)設(shè)計(jì),應(yīng)用于不同測線方位角觀測系統(tǒng)疊加得出:

1)三維地震勘探區(qū)內(nèi)大型障礙物影響觀測系統(tǒng)中接收點(diǎn)和激發(fā)點(diǎn)的布設(shè),采用不同測線方位角觀測系統(tǒng)疊加的施工方案,在實(shí)際施工中存在觀測系統(tǒng)面元中CMP位置分散的情況,但面元大小的設(shè)計(jì)可以有效控制該情況對于數(shù)據(jù)的影響程度,對于最終的采集數(shù)據(jù)質(zhì)量無顯著影響。

2)基于貪心算法編寫的特殊觀測系統(tǒng)設(shè)計(jì)優(yōu)化程序,對于不同測線方位角的觀測系統(tǒng)交疊下特殊觀測系統(tǒng)設(shè)計(jì)優(yōu)化計(jì)算合理,設(shè)計(jì)經(jīng)濟(jì)適用,提高了施工效率。經(jīng)實(shí)際應(yīng)用,獲得了大型障礙物下地震反射數(shù)據(jù),保障了采集質(zhì)量。