彭玉林 吳佳樂 劉宜文 陳景亮 鄒梅 侯奇 高雪江

摘? 要：準噶爾盆地是新疆油田公司地震勘探的主戰(zhàn)場，其復雜多變的近地表條件給地震勘探帶了很大困難。建立精確的近地表模型，掌握近地表結構變化規(guī)律，對地震勘探3大環(huán)節(jié)（采集、處理、解釋）的研究至關重要。在分析DEM數(shù)字高程數(shù)據(jù)特征基礎上，研發(fā)DEM數(shù)字高程校正方法，實現(xiàn)與實測地表高程的無縫拼接，填補了沒有實測高程區(qū)域的空缺，為精確建立近地表結構模型奠定了基礎。在此基礎上，利用校正后的DEM數(shù)字高程建立盆地級高精度近地表結構模型，為后續(xù)指導新探區(qū)表層調(diào)查點的布設、地震采集激發(fā)井深的精確設計等發(fā)揮了重要作用。

關鍵字：準噶爾;表層結構; 區(qū)域建模;靜校正;表層部署

分層建模法是近地表結構建模常用方法之一，相比利用地震單炮初至時間折射和層析反演方法而言，分層法具工作量小、占用計算機時少、建模效率高等特點，是目前盆地級區(qū)域建模的主要技術方法[1-4]。該方法是在獲取目標區(qū)域地表高程數(shù)據(jù)前提下，通過區(qū)域內(nèi)表層調(diào)查點（小折射、微測井）的垂向厚度、速度等成果信息進行空間插值，得到區(qū)域內(nèi)每個物理點的表層結構模型[5-8]。以往地震勘探測網(wǎng)間存在不同密度的地表高程數(shù)據(jù)空白區(qū)，導致區(qū)域模型整體建立存在一定問題。為彌補地表高程缺失問題，以往直接利用地表衛(wèi)片影像中的DEM數(shù)字高程代替區(qū)域測網(wǎng)中的地震實測高程，由于DEM數(shù)字高程與地震實測高程之間存在誤差，直接使用DEM數(shù)字高程不能滿足區(qū)域近地表結構精細建模需求。因此，本文提出一種在高精度數(shù)字高程的基礎上，利用實測高程網(wǎng)格化過渡填充方式進行校正，獲得盆地區(qū)域精細地表高程的方法。本方法通過設置過渡區(qū)域和網(wǎng)格化方式，保證數(shù)字模型高程和實測物理高程無痕過渡銜接，確保DEM數(shù)字高程與實測高程之間的閉合，為盆地級近地表模型構建奠定了數(shù)據(jù)基礎。

1? 研究思路

本次研究應用的準噶爾盆地DEM數(shù)字高程基礎數(shù)據(jù)，來源于ASTER GDEM V2版本的DEM數(shù)字高程數(shù)據(jù)，經(jīng)投影轉換、格式調(diào)整和重新網(wǎng)格化后生成20×20 m平面網(wǎng)格的高程節(jié)點，是盆地統(tǒng)一高程的初始網(wǎng)格數(shù)據(jù)。由于DEM數(shù)字高程生產(chǎn)過程中，原始數(shù)據(jù)采集受原始資料、采點設備及采集過程中人為因素影響，導致原始數(shù)據(jù)存在誤差，且DEM建立時進行內(nèi)插計算與建模過程會產(chǎn)生內(nèi)插建模誤差，該誤差不僅與內(nèi)插算法相關，同時與原始數(shù)據(jù)分布密度相關。另投影轉化和平面網(wǎng)格化過程也存在一定誤差[9-12]，這些誤差的累積，最終導致在含有地震測量數(shù)據(jù)的區(qū)域內(nèi)，相同網(wǎng)格內(nèi)的DEM數(shù)字高程與地震實測高程存在最大約20 m的閉合誤差，導致以此高程數(shù)據(jù)為基礎建立的盆地區(qū)域近地表模型發(fā)生一定程度突變（圖1），影響模型的整體應用。為保證盆地區(qū)域近地表結構建模精度，需對盆地初始DEM數(shù)字高程進行以實測地震高程為參考點的平面校正，消除兩套地表高程數(shù)據(jù)的閉合差。

1.1? 網(wǎng)格化過渡填充方式對DEM數(shù)字高程進行校正

數(shù)據(jù)網(wǎng)格化是將空間上不均勻分布的數(shù)據(jù)，按一定方法，如滑動平均法、克里格法或其他適當數(shù)值推算、插值等方法，計算規(guī)則網(wǎng)格中的代表值（趨勢值）的過程。數(shù)據(jù)網(wǎng)格化的基本功能是在遵循所研究變量的空間變化趨勢的基礎上，將空間上分散的數(shù)值換算成規(guī)則分布的網(wǎng)格數(shù)值，彌補空白網(wǎng)格的數(shù)值，為不同變量的綜合分析及對比提供統(tǒng)一的空間結構[12-14]。

為在準噶爾盆地范圍內(nèi)獲得精細、可靠、連續(xù)、全覆蓋的高程基準數(shù)據(jù)，需將采集的物理點高程和DEM數(shù)字高程數(shù)據(jù)相結合，形成統(tǒng)一高程網(wǎng)格數(shù)據(jù)?？紤]到采集區(qū)塊內(nèi)物理點高程為實測高程數(shù)據(jù)，可信度高，因此，以物理點高程作為DEM數(shù)字高程校正的基準數(shù)據(jù)。在物理點可控范圍內(nèi)，對DEM數(shù)字高程進行校正，保證校正后的DEM數(shù)字高程與物理點高程一致，并在缺少物理點高程控制的區(qū)域內(nèi)使用DEM數(shù)字高程數(shù)據(jù)進行網(wǎng)格化填充，達到提高DEM數(shù)字高程精度的目的。本次以自主研發(fā)的集盆地數(shù)據(jù)管理、應用、質控于一體的《近地表結構解釋系統(tǒng)》為平臺，開展DEM數(shù)字高程校正研究。該系統(tǒng)能方便快捷地獲取盆地所有區(qū)域部署采集的物理點高程和初始的DEM數(shù)字高程（20×20 m平面網(wǎng)格），并參考DEM網(wǎng)格邊長等參數(shù)對地震測量成果進行網(wǎng)格化，生成與初始DEM數(shù)字高程網(wǎng)格坐標一致的中間實測物理點高程網(wǎng)格平面數(shù)據(jù)體。在此網(wǎng)格基礎上，去除沒有實測物理點網(wǎng)格數(shù)據(jù)覆蓋的空白區(qū)域，獲得物理點實測高程與初始DEM數(shù)字高程誤差平面數(shù)據(jù)體。DEM數(shù)字高程校正的目的就是利用此誤差數(shù)據(jù)體對初始DEM數(shù)字高程進行平面校正，整體校正基本可分為以下3種：

地震實測高程點包絡線以外區(qū)域? 據(jù)預先指定的過渡區(qū)域尺度，在中間高程網(wǎng)格的空值區(qū)域劃分標記出過渡區(qū)域。將過渡區(qū)域外的空值網(wǎng)格點坐標代入全盆地DEM數(shù)字高程網(wǎng)格，提取DEM數(shù)據(jù)并替換對應的網(wǎng)格點數(shù)值。此時，中間高程網(wǎng)格點可分為3個部分：有效物理點控制范圍內(nèi)的有值網(wǎng)格點（W區(qū)）;標記為過渡區(qū)域的空值網(wǎng)格點（Q區(qū)）;過渡區(qū)域外按DEM數(shù)值填充的有值網(wǎng)格點（M區(qū)）。并將在平面網(wǎng)格中實測高程與中間網(wǎng)格DEM數(shù)字高程之間的誤差標記為[Di]，按上述區(qū)域劃分，可知W與M區(qū)[Di]均為已知，其中M區(qū)域[Di]為零，Q區(qū)[Di]需通過插值計算。為避免過渡區(qū)域內(nèi)出現(xiàn)“臺階”似突變現(xiàn)象，采用距離加權平均方式對Q區(qū)網(wǎng)格點[Di]值進行空間插值。

以給定搜索半徑內(nèi)與目標點距離最近的W區(qū)3個控制點及M區(qū)3個控制點為基點，在三角網(wǎng)格中計算任意一網(wǎng)格空值點的高程誤差。假設與目標點J最近的6個控制點分別為[i]、[i+1]、[i+2]、[i+3]、[i+4]、[i+5]（圖2左），此6個點與相鄰網(wǎng)格點的高程誤差分別為[Di]、[Di+1]、[Di+2]、[Di+3]、[Di+4]、[Di+5]，則J點的誤差[Di]可表示為：

其中，[Ki]為[i]點對[j]點的權系數(shù)，與距離相關，可表示為：

上式中[xi]為[i]點與[j]點的直線距離[13-14]。由此，可得到Q區(qū)域中所有網(wǎng)格點的高程誤差值，與已知W區(qū)和M區(qū)中高程誤差值相結合，即得到全區(qū)域的高程誤差值。

地震實測高程點包絡線以內(nèi)區(qū)域? 劃分為有效物理點控制范圍內(nèi)的有值網(wǎng)格點（W區(qū)）與標記為過渡區(qū)域的空值網(wǎng)格點（Q區(qū)），Q區(qū)的[Di]需利用W區(qū)的[Di]進行空間插值方式計算。

由于盆地地震實測物理點近4 000萬個，全部用于空間插值的計算量巨大。盆地中同樣為實測地表高程的表層調(diào)查點均按一定密度間隔分布在地震測網(wǎng)中，總數(shù)據(jù)點為76萬個，遠少于實測地震數(shù)據(jù)點，用其作為插值基準點將大幅減小計算數(shù)據(jù)量。因此，在實測高程點包絡線以內(nèi)區(qū)域，可利用平面網(wǎng)格中表層控制點的實測高程與初始DEM數(shù)字高程之間的誤差值進行空間插值，獲得完整的高程誤差平面圖，實現(xiàn)DEM數(shù)字高程的整體校正。

以給定搜索半徑內(nèi)與目標點距離最近的3個控制點為基點，在三角網(wǎng)格中計算任意一網(wǎng)格空值點的高程誤差。假設與目標點J最近的3個控制點分別為[i]、[i+1]、[i+2]（圖2右），此3個點與相鄰炮-檢點的高程誤差分別為[Di]、[Di+1]、[Di+2]，則[j]點的誤差[Di]表示為：

其中，[Ki]為[i]點對[j]點的權系數(shù)，與距離相關，可表示為：

上式中[xi]為[i]點與[j]點的直線距離[15]。由此，可得到Q區(qū)域中所有網(wǎng)格點的高程誤差值，與已知W區(qū)中的高程誤差值相結合，即得到全區(qū)域內(nèi)高程誤差值。

DEM數(shù)字高程校正? 當?shù)玫絒j]點高程誤差后，可由下式得到[j]點經(jīng)校正之后的高程值[Sj]：

上式中[Sj]為[j]點需整體整合的DEM數(shù)字高程，[RSi]為[j]點實際DEM數(shù)字高程值。同理，可得到實測物理點校正后的DEM數(shù)字高程值，這樣就完成了用物理點實測高程對DEM數(shù)字高程的整體校正。這種高程校正由于在平面上受表層調(diào)查點密度和插值半徑約束，只利用實測物理點可靠的低頻成分對DEM數(shù)字高程的低頻成分進行校正，保留了DEM數(shù)字高程可信的高頻成分。因此，這種高程校正具良好的適應性，校正結果具實測物理點高程與DEM數(shù)字高程的可靠低、高頻成分。

1.2? 雙緩存技術提高DEM高程校正效率

由前述可知，平面區(qū)域上每個物理網(wǎng)格點的高程誤差計算是以網(wǎng)格數(shù)據(jù)的形式進行的。由于全盆地物理點數(shù)據(jù)量非常大，僅從《近地表結構解釋系統(tǒng)》加載整個探區(qū)物理點高程數(shù)據(jù)會一次性向內(nèi)存提交近10G數(shù)據(jù)量，加上全探區(qū)高程數(shù)據(jù)網(wǎng)格化流程需占用大量系統(tǒng)資源，故在數(shù)字高程校正計算時所申請使用的內(nèi)存很容易超出系統(tǒng)內(nèi)存容量，導致客戶端程序崩潰。本次采用以時間換空間的解決思路，不僅在加載全探區(qū)采集項目物理點數(shù)據(jù)時使用本地緩存文件代替內(nèi)存空間來存儲預處理過的高程數(shù)據(jù)，還在全探區(qū)高程網(wǎng)格生成和提取過程中采取動態(tài)分塊和本地緩存相結合的方式，避免校正計算過程中出現(xiàn)一次占用過多系統(tǒng)資源，導致程序無法正常工作。這樣處理雖會增多加載、釋放分塊緩存數(shù)據(jù)的時間，但可解決全探區(qū)大數(shù)據(jù)量無法順利計算的問題，同時可支持更精細的探區(qū)級網(wǎng)格數(shù)據(jù)生成與提取[16-17]。

2? 應用效果

2.1? DEM數(shù)字高程校正的效果

通過本文所述技術的綜合應用，解決了以往地震測網(wǎng)中存在實測高程空白、數(shù)字高程與地震實測高程存在閉合差等問題，為實現(xiàn)盆地級高精度近地表模型的整體建立提供支撐。圖3為校正前后地震實測高程與DEM數(shù)字高程間的誤差對比。如圖所示，在校正前DEM數(shù)字高程與物理點實測高程之間存在明顯不閉合問題，校正后的DEM數(shù)字高程與實測高程閉合良好（誤差小于2 m），所得盆地區(qū)域級DEM數(shù)字高程精度可靠。在以此數(shù)據(jù)為基礎建立的盆地級近地表結構模型基礎上，提取盆地級近地表結構剖面（圖4），相比圖1來看，當DEM數(shù)字高程精度提高后，之前在地震測網(wǎng)空白區(qū)存在的模型閉合差問題得到有效解決，模型精度進一步提高。

2.2? 指導表層調(diào)查點的部署

在獲取盆地級高精度DEM數(shù)字高程后，進行區(qū)域表層調(diào)查點采集部署時，改變以往測網(wǎng)建模方式為區(qū)域整體模型構建方式（（圖5左-右）。以往方式只能沿二維測網(wǎng)對表層結構進行縱橫向變化分析，不能對重要目標區(qū)進行任意線的連井表層結構分析，導致難以對目標區(qū)表層結構變化進行針對性的表層調(diào)查部署。進行區(qū)域模型構建后，能精細地分析任意區(qū)域的表層結構變化，針對性地部署表層調(diào)查點。針對探井A、B間的結構變化關系，改變以往均勻部署調(diào)查點的方式，可在重要結構突變區(qū)域進行針對性部署，在提高部署效果的同時有效降低了表層調(diào)查點采集成本[18-20]。

2.3? 精準設計激發(fā)井深

地震采集過程中，井炮激發(fā)井深的選取非常重要，關系到子波的頻帶寬度、反射信號能量。恰當?shù)木钗恢每蓽p少震源產(chǎn)生的次生干擾、虛反射等。通過對準噶爾盆地地震勘探多年來井炮激發(fā)井深的認知，最佳激發(fā)深度應選擇低降速帶界面之下的高速層，盡量靠近高速層頂界面，這實際上對地震采集的井深設計精度提出了很高要求。

以往激發(fā)井深設計是依據(jù)工區(qū)內(nèi)歷史表層調(diào)查資料進行近地表結構模型建立，確定每個激發(fā)點位置的井深，由于缺少實際的地表高程，地表高程的變化又與低降速層界面的變化具不完全相關性，激發(fā)位置是否在高速層頂界面之下就很難保證，其誤差大小依賴于高程精度。目前獲取精確盆地DEM數(shù)字高程，建立能準確反映表層結構空間垂向變化的近地表模型，可為精準設計激發(fā)井深奠定基礎。圖6為依據(jù)本文方法校正后的DEM數(shù)字高程建模前后所得近地表模型設計的激發(fā)深度試驗單炮效果對比。受空間地表高程精度影響，前后設計的激發(fā)深度相差5 m，但單炮效果相差較大。圖6左為在高速頂界面之上3 m處激發(fā)，激發(fā)介質速度1 600 m/s;圖6右為在高速層頂界面之下2 m處激發(fā)，激發(fā)介質速度為2 000 m/s。從圖中可看到，高速層中激發(fā)的單炮從能量、信噪比、頻率等特征均優(yōu)于高速層之上激發(fā)的單炮效果。

3? 結論及認識

（1） 實際資料應用表明，本文所述盆地DEM數(shù)字高程校正方法具良好的穩(wěn)定性和精度。此方法結合盆地地震實測物理點并以盆地表層調(diào)查點實測高程為基準網(wǎng)格，對以往直接獲取的初始DEM數(shù)字高程進行了低頻誤差校正，較好地提高了建模所需高程數(shù)據(jù)的精度，為盆地級近地表模型構建奠定了堅實的數(shù)據(jù)基礎。

（2） DEM數(shù)字高程校正技術的實現(xiàn)，為地震采集激發(fā)井深的精確設計提供了基礎，確保了接近高速層頂界面位置激發(fā)，單炮記錄優(yōu)質品率明顯提高。

（3） 本次研究和推廣應用的DEM數(shù)字高程校正是以現(xiàn)有的初始DEM高程數(shù)據(jù)體、盆地地震測網(wǎng)和表層調(diào)查點數(shù)據(jù)為基準，當后續(xù)隨精度進一步提升后，仍需進行重新校正，以保證數(shù)據(jù)更新的及時性。

參考文獻

[1]? ? 張恒超，李學聰. 模型法、擴展廣義互換法（EGRM）、沙丘曲線法靜校正對比研究及在ZG沙漠資料處理中的應用[J]. 地球物理學進展，2010，25（6）：2193-2198.

[2]? ? 王進海，熊民生. 模型法靜校正在黃土塬地區(qū)的應用[J]. 石油地球物理勘探，1995（S1）：48-53+185.

[3]? ? 夏忠謀，何曉冬. 模型法靜校正在準噶爾盆地腹部地區(qū)的應用[J]. 中國西部油氣地質，2005（2）：203-206.

[4]? ? 李衛(wèi)忠，Saleh A.Al-Maghlouth， Adnan Hashem，等. 基于模型的復雜地表區(qū)靜校正方法[J]. 石油地球物理勘，2010，45（6）：797-801+936+783.

[5]? ? 張立雙. 模型法靜校正技術研究與應用[J]. 化工管理，2015（35）：194.

[6]? ? 梁承敏，薛為平，邵雨. 綜合建模靜校正技術在低幅度構造研究中的應用[C] //CPS /SEG北京 2003國際地球物理會議暨展覽， 北京， 2003.

[7]? ? 李衛(wèi)忠.復雜地表地區(qū)近地表模型與靜校正[J].油氣地質與采收率，2007，14（3） ： 73-77.

[8]? ? 劉立平.準噶爾盆地南緣山前帶靜校正技術研究[J].中國石油大學勝利學院學報， 2016，30（5）：6-7.

[9]? ? 張甲，張俊濤. DEM誤差與精度分析研究[J].中國科技縱橫，? ? ? ? 2014（20）：225.

[10]? 胡鵬，黃雪蓮，吳艷蘭，等. DEM若干理論問題思考[J]. 哈爾濱工業(yè)大學學報， 2006 （12）：2143-2147.

[11]? 王耀革，王志偉，朱長青. DEM誤差的空間自相關特征分析[J]. 武漢大學學報（信息科學版），2008（12）：1259-1262.

[12]? 田瑞云，王玉寬，傅斌，等. 基于DEM的地形單元多樣性指數(shù)及其算法[J]. 地理科學進展， 2013，32（1）：121-129.

[13]? 張偉，覃慶炎，簡興祥. 自然鄰點插值算法及其在二維不規(guī)則數(shù)據(jù)網(wǎng)格化中的應用[J]. 物探化探計算技術， 2011，33（3）：291-295+? 228.

[14]? 張春玲，王振峰，李宏. 基于最優(yōu)插值法的Argo數(shù)據(jù)網(wǎng)格化實驗[J]. 海洋測繪，2012， 32（3）：29-31.

[15]? 汪隆六. 三角形線性插值在區(qū)域重力測量數(shù)據(jù)網(wǎng)格化中的應用[J]. 物化探計算技術，1986（2）：141-146.

[16]? 陳柏生，杜吉祥. 用于超級節(jié)點網(wǎng)絡的雙緩存技術[J].華中科技大學學報（自然科學版）， 2011 （1）：24-27.

[17]? 劉婷婷，潘彤陳，彬茹.利用雙緩存技術提高繪圖速率[J]. 電腦編程技巧與維護， 2012（18） ：22，47.

[18]? 夏竹，張少華，王學軍.中國西部復雜地表地區(qū)的近地表特征與表層結構探討[J]. 石油地球物理勘探， 2003， 38（4）： 414-424.

[19]? 劉治凡，毛海波，邵雨，等. 復雜地表區(qū)基準面和靜校正方法的選擇[J]. 石油物探， 2003 ， 42（2）： 240-246.

[20]? 楊德義，彭蘇萍，常鎖亮.地震勘探常規(guī)低速帶調(diào)查方法應用中的問題探討[J].物探與化探， 2003， 27（5）：387-390.