伊萬順 徐 海 都小芳 趙 夏 岳震國 黃 衛(wèi)
(1.中國石油大學CNPC物探重點實驗室,北京102249;2.中國地質(zhì)大學 能源學院,北京100004;3.東方地球物理公司,河北 涿州072751)
傳統(tǒng)意義上的復雜地區(qū)三維地震資料采集的炮點偏移設計,通常是通過手工進行偏點的偏移,僅適合低復雜度與偏點數(shù)量較小的情況,效率低下而且準確度低[1,2]。近年來隨著油氣勘探進一步深入,極復雜地區(qū)(如城區(qū)、海陸過渡帶等)成為勘探的熱點,手工炮點的偏移設計遠不能滿足生產(chǎn)的需要,如某些工業(yè)區(qū)與城區(qū)內(nèi)的線束,其偏點(偏移的炮點)數(shù)量一般占理論炮點的30%~50%,采用人工方法進行偏移設計相當困難。第一,影響設計與施工效率,增加野外施工的成本,不能保證偏點炮的有效性;第二,不能保證有效炮數(shù)量的最大化以及覆蓋次數(shù)的最優(yōu)化;第三,不能保證后期資料提交的準確性[3,4]。本文根據(jù)實例,著重研究實現(xiàn)了一種以衛(wèi)星照片、地形圖、勘察測量數(shù)據(jù)等地理信息為基礎(chǔ),進行地震勘探部署設計、現(xiàn)場采集觀測設計的計算機自動優(yōu)化技術(shù)和方法,為評估野外采集方法、降低勘探成本、避免失誤風險、優(yōu)化地震資料采集觀測效果提供更為科學的分析技術(shù)和手段。
以國外X地區(qū)的海陸過渡帶三維地震勘探項目為例。該項目規(guī)模大,陸地部分地表復雜,分布著密集的管線群、高速公路、石油化工廠、污水處理廠、飛機場、居民區(qū)等。地震采集存在大量空點(不能放炮的理論設計點),其中,某些線束偏移炮點(以下簡稱為偏點)數(shù)量高達50%。針對這種地面條件,如何既保證炮點布設的合理性,又保證炮點的可實施性,最大限度地滿足覆蓋次數(shù)均勻與施工效率的要求,是目前野外地震采集中所面臨的最大挑戰(zhàn)。
X區(qū)設計規(guī)則為:橫向偏移距≤750m,縱向偏移距≤800m;其觀測系統(tǒng)設計為12線24炮,炮線點距100m×50m,檢波線點距200m×50 m,最大偏移距為6km,滿覆蓋為360次。
根據(jù)偏點規(guī)則,嘗試了以下偏移方法與匹配方法。以炮線為參考系定義方向,平行于檢波線方向的炮線稱為炮線的橫向,垂直于檢波線方向的炮線稱為炮線的縱向。
1.2.1 偏移方法與分析
a.先橫后縱的串聯(lián)查找(圖1)。平行于檢波線方向的橫向匹配,然后再從垂直于檢波線的炮線方向進行查找;縱向距離取絕對值的最小距離,橫向采用最優(yōu)距離。
圖1 先橫向后縱向查找Fig.1 Inline-cross line search
b.橫向-橫縱并聯(lián)查找(圖2)。橫向查找完畢后,再由橫向與縱向的方向同時查找(藍色框內(nèi)查找);縱向距離取絕對值的最小距離,橫向采用最優(yōu)距離。
c.橫縱的并聯(lián)查找(圖3)。即橫向與縱向同時查找,只要滿足偏移條件,就認為匹配的偏點是有效的。
圖2 先橫向再橫縱向查找Fig.2 Inline-inline &cross line search
圖3 橫縱向同時查找Fig.3 Inline &cross line search
根據(jù)以上偏點的查找思路,對這3種查找方法應用實際數(shù)據(jù)進行了試驗分析(表1),結(jié)論如下。
表1 3種匹配方法的對比試驗結(jié)果Table 1 The result of the comparison of three methods
第一種與第二種偏移查找方法皆符合地球物理偏點思路,能夠最大化地滿足空點區(qū)域CMP疊加剖面的覆蓋次數(shù);但是第一種方法找到的匹配點數(shù)要比第二種方法少,原因是僅考慮了各個方向的單向查找。
第三種與前兩種相比,偏點的查找數(shù)量最多(與第二種方法相比,實際處理數(shù)據(jù)多出10~30個點,視具體情況而定),但是縱向偏移的炮點太多,不完全符合地球物理學家的偏點思路。
1.2.2 匹配方法與分析
假設在橫向匹配過程中偏移距離不能超過250m,抽象模型算法與分析如下。
a.最小絕對值法。在橫向匹配過程中,當空點與偏移點之間距離的絕對值最小時,認為此偏點與空點的匹配形式是最佳的(圖4);但其中僅有3個點能夠匹配,其余2個點已經(jīng)超過偏移距離250m的限制。最小絕對值的匹配方式:紅色點代表可實施設計點,藍色點代表偏點,紫色點代表空點。
b.最優(yōu)距離法。在橫向匹配過程中,不以距離最小絕對值為判定標準,而是以怎樣最大化地增加空點區(qū)域覆蓋次數(shù)為目的的偏點匹配形式(圖5),其中僅有5個點皆能夠匹配,沒有超過偏移距離250m的限制。最優(yōu)距離匹配方式:紅色點代表可實施設計點,藍色點代表偏點,紫色點代表空點。
圖4 最小絕對值匹配方法示意圖Fig.4 Least absolute value match
圖5 最優(yōu)距離匹配方法示意圖Fig.5 Optimization distance match
通過模型的理論分析以及實際數(shù)據(jù)處理后,發(fā)現(xiàn)第二種匹配方法最佳(表1),找到匹配點較多,符合地球物理學家的偏點思路。
通過方法分析,選取了最優(yōu)距離與橫向-橫縱并聯(lián)查找方法組合,進行了編程實現(xiàn)[5-7]。
設SPS炮點文件中設計炮點文件為aij(i=0,1,2,…,m1;j=0,1,2,…,n1);不存在的設計點為bji(i=0,1,2,…,m2;j=0,1,2,…,n2),存在的偏點為cji(i=0,1,2,…,m3;j=0,1,2,…,n3)。其中bji與cji為經(jīng)過二分法遞歸排序得到的結(jié)果,此時令cji中的j為一重循環(huán),bji中的j為二重循環(huán),實現(xiàn)沿橫向查找偏點所匹配的空點??拯c的橫向查找方式以max(abs(cji-bij))為準,橫縱向以 min(abs(cji-bij))距離為準,實現(xiàn)偏點的匹配。
運行發(fā)現(xiàn)該方法配合,查找最合理,雖然找到的匹配點沒有橫縱并聯(lián)查找的個數(shù)多,但這種查找方式更符合地球物理學家面元覆蓋次數(shù)均一化的設計思路(圖6)。
基于衛(wèi)星遙感圖片可將其作為踏勘依據(jù),設計不同地形條件下的激發(fā)點,組織質(zhì)控與測量人員提前踏勘,對激發(fā)點進行可行性分析,按照計劃去核實炮點的可實施性,對不合適的炮點還需要二次偏移。在商業(yè)區(qū)、居住區(qū)及城鎮(zhèn)等大面積復雜障礙區(qū),要根據(jù)觀測系統(tǒng)的特點以及地形地物特點,通過模擬分析,在施工前需提前確定偏移炮點的分布(如偏移方向與偏移距離等),通過室內(nèi)與野外密切結(jié)合,實時地使室內(nèi)與野外設計達到最優(yōu)化。
圖6 偏點優(yōu)化流程Fig.6 Offset shot optimization flow
以橫跨城區(qū)與工業(yè)區(qū)的某線束為例,其地表分布著密集的管線群、高速公路、石油化工廠、污水處理廠、居民區(qū),地表情況極其復雜(圖7)。為滿足目標區(qū)各目的層覆蓋次數(shù)以及最大限度地降低空炮率,需要布設大量的偏移炮點。該線理論施工(假定地表地形規(guī)則,沒有障礙物)設計炮點4 723炮(圖8)?,F(xiàn)場地球物理專家根據(jù)地表地形條件實際設計了5 306炮點(圖9),在非優(yōu)化條件下,比理論多設計了583炮(5 306-4 723=583)。若按此方案激發(fā),在考慮覆蓋次數(shù)均勻的情況下,不但不能保證偏點的有效性,而且還會產(chǎn)生大量的無效炮,影響施工效率,增加施工成本。
據(jù)實際設計,在城區(qū)與工業(yè)區(qū)的無障礙物處可施工的設計點(正點)應為2 839個(圖10中粉紅色炮點),障礙物處的空點和對應有效偏點應為1 884個(圖10中藍色炮點)。
圖7 工區(qū)地面衛(wèi)星圖Fig.7 High resolution satellite photo
據(jù)偏移規(guī)則,借助偏移優(yōu)化算法的自動匹配,對踏勘設計炮點進行優(yōu)化處理與分析后得到如下結(jié)果:在踏勘初始設計方案中有89個空點(圖11中藍色炮點)不能在所有可利用的設計偏點中找到其對應的偏移炮點,672個設計偏移點找不到與其對應的障礙物空點(圖11中粉紅色炮點);經(jīng)優(yōu)化處理后被廢除,最終優(yōu)化處理后的有效偏移炮點個數(shù)為1 795個(表2,表3)。
圖8 理論設計炮點分布圖Fig.8 Theoretical designed shot points distribution
圖9 踏勘初始設計炮點分布圖Fig.9 The scouting designed shot points distribution
圖10 優(yōu)化后實際炮點分布圖Fig.10 The shot points distribution after optimizing正點:粉紅色點;空點:藍色點
圖11 廢棄偏點與不匹配空點分布圖Fig.11 Distribution of the abandoned offset shot and unmatched points
表2 理論點、偏點、正點、空點的關(guān)系Table 2 The relation between theoretical,offset,enforceable theoretical and skipped points
表3 設計點、正點、偏點、空點統(tǒng)計關(guān)系Table 3 The Statistical relation between theoretical,offset,enforceable theoretical and skipped points
通過對偏點優(yōu)化設計方法的討論以及國外X復雜過渡帶三維地震勘探項目炮點施工設計的應用,得到以下結(jié)論。
a.在匹配數(shù)量以及設計觀念上,經(jīng)過試驗分析可知,基于橫-橫縱向查找,并以橫向最優(yōu)距離、縱向最小距離的匹配方法是最優(yōu)化的。
b.通過偏點優(yōu)化程序?qū)碗s區(qū)施工方案進行反復設計、不斷優(yōu)化,最大限度地提高了覆蓋次數(shù)。
c.將數(shù)量巨大的偏點設計從手工設計中解放出來,設計人員可隨意盡可能多地布設可以放炮的偏點,不需考慮偏移規(guī)則,處理后的多余偏點可通過程序自動剔除,提高了質(zhì)控人員偏點設計的工作效率。
d.降低了因不合理偏移點而產(chǎn)生的無效炮數(shù),節(jié)約了成本,提高了野外作業(yè)的工作效率。
e.解決了偏點與空點的相互匹配問題,降低了資料整理人員的工作量,提高了工作效率。
[1]鄔倫,劉瑜,張晶,等.地理信息系統(tǒng)——原理、方法和應用[M].北京:科學出版社,2001.
[2]Beasley C J,Vermeer G J O.3-D Seismic Survey De-sign[M].Tulsa:Soc of Exploration Geophysicists,2002.
[3]Telford W M,Geldart L P.Applied Geophysics Second Edition[M].New York:Cambridge University Press,2004:186-190.
[4]易遠元,雷鳴,安俊.基于地理信息的三維地震資料采集設計[J].科技經(jīng)濟市場,2009(10):29-32.
[5]莊道川,曾慶才,張亞斌,等.超大偏移距地震資料處理方法[J].西南石油學院學報,2003,25(2):11-14.
[6]盧開澄.單目標、多目標與整數(shù)規(guī)劃[M].北京:清華大學出版社,1997.
[7]袁亞湘,孫文瑜.最優(yōu)化理論與方法[M].北京:科學出版社,1997.