熊晶璇 賀振華 劉 鴻 趙容容 楊衛(wèi)寧 張恩嘉 王光銀 唐 虎

1.成都理工大學(xué) 2.中國石油天然氣集團公司東方地球物理勘探有限公司3.中國石油天然氣集團公司山地地震技術(shù)試驗基地 4.中國石油西南油氣田公司勘探事業(yè)部

0 引言

前陸盆地發(fā)育在穩(wěn)定克拉通與線性收縮造山帶的結(jié)合部位，有著鮮明的構(gòu)造—沉積特征，油氣成藏條件優(yōu)越[1-4]。眾多學(xué)者論述了前陸盆地是重要的含油氣單元[5-8]。近年來，塔西南、準(zhǔn)噶爾南緣，川西南等前陸盆地油氣勘探取得重要進展[9-13]，也證實了該領(lǐng)域巨大的勘探潛力。前陸盆地油氣資源豐富，但其地下構(gòu)造復(fù)雜，多存在高陡構(gòu)造上下盤展布差異大，斷裂發(fā)育，地層傾角變化大，構(gòu)造倒轉(zhuǎn)[1-2,6-7]，目標(biāo)地質(zhì)體埋藏深、尺度小等現(xiàn)象[14-17]，在速度模型上表現(xiàn)為速度縱橫向變化大，普遍存在速度倒轉(zhuǎn)和局部高頻等現(xiàn)象、山地復(fù)雜構(gòu)造的深度域速度模型建立和更新等問題成為業(yè)界難題[17]。因此，針對性研究適合復(fù)雜山地油氣勘探的深度域速度建模技術(shù)具有重要意義。

初始層速度模型建立是深度域速度建模的第一步。在山地復(fù)雜構(gòu)造區(qū)域，因地下地質(zhì)情況復(fù)雜造成速度場復(fù)雜，使深度域速度模型的迭代更新更依賴于初始層速度模型[18]，高質(zhì)量的初始層速度模型既可加快速度更新的收斂，又能進一步提高偏移成像的效果和效率。經(jīng)典Dix反演是常用的初始層速度模型建立方法，它建立在水平層狀或小傾角地層與速度橫向無變化假設(shè)的理論基礎(chǔ)上[19]，但山地復(fù)雜構(gòu)造時距曲線方程多為非雙曲線型，經(jīng)典Dix反演抗干擾能力弱，提供的層速度變得不可靠[20]。為了獲得可靠的層速度，Hubral等學(xué)者對經(jīng)典Dix公式進行了研究[21-23]，Clapp等學(xué)者針對Dix 公式假設(shè)[24-25]，提出了改進的約束Dix反演層速度的求取方法。Koren 和Ravve提出了約束層速度反演 CVI 方法[26]，進一步提高速度模型的橫向連續(xù)性[27-29]。但如前所述，山地復(fù)雜構(gòu)造區(qū)速度場復(fù)雜，CVI方法作為常用方法，通過單一固定的影響域?qū)λ芯植繕?gòu)造進行約束，反演結(jié)果精確性的風(fēng)險，甚至?xí)刖植空`差， 并在反演中將誤差傳遞放大，影響結(jié)果的正確性。除此之外，常規(guī)算法中對速度的離散數(shù)據(jù)進行網(wǎng)格化時，多基于速度體橫向變化平緩，在山地復(fù)雜地質(zhì)區(qū)，因存在構(gòu)造倒轉(zhuǎn)和局部異常地質(zhì)體帶來的速度倒轉(zhuǎn)和局部速度高頻等現(xiàn)象，常規(guī)的插值平滑方法會降低速度模型對倒轉(zhuǎn)構(gòu)造和小規(guī)模地質(zhì)體的靈敏性，增加速度更新迭代收斂的次數(shù)，甚至?xí)硭俣饶Ｐ偷呐で?，影響疊前深度偏移的成像質(zhì)量。

本文通過形成構(gòu)造約束的速度趨勢函數(shù)和基于徑向基函數(shù)薄板理論的三維空間速度網(wǎng)格重構(gòu)函數(shù)，在常規(guī)約束速度反演研究的基礎(chǔ)上，形成了山地復(fù)雜構(gòu)造區(qū)初始層速度改進反演方法，該方法在復(fù)雜構(gòu)造區(qū)域反演得到的層速度場更符合地質(zhì)規(guī)律，反演效果優(yōu)于常規(guī)方法。該方法研究成果有助于推進我國復(fù)雜構(gòu)造區(qū)油氣勘探開發(fā)。

1 方法原理

常規(guī)約束速度反演方法的目標(biāo)函數(shù)由均方根速度擬合、指數(shù)漸進邊界速度趨勢函數(shù)擬合（以下簡稱趨勢函數(shù)擬合）和防振蕩阻尼擬合三部分組成[26,30]。其中，均方根速度擬合是均方根速度對反演結(jié)果的貢獻。防振蕩阻尼擬合用來壓制界面上下縱向速度梯度的過大跳躍。速度趨勢函數(shù)擬合用于改善Dix反演帶來的高頻振蕩，影響甚至決定了反演結(jié)果的合理性和穩(wěn)定性。Ravve[31]等提出兩個粗網(wǎng)格點間用指數(shù)形式的趨勢速度填充，對改善經(jīng)典Dix反演時出現(xiàn)的高頻振蕩，具有重要的貢獻意義。速度趨勢函數(shù)即指數(shù)漸進邊界速度趨勢函數(shù)（Exponential, Asymptotically Bounded Model，簡稱 EAB 函數(shù)），其數(shù)學(xué)表達式為[31-33]：

Ravve[31]等進一步提出一定影響半徑內(nèi)的粗網(wǎng)格上所有的垂直函數(shù)都相互影響，影響權(quán)重與距離成反比。影響域具有一定的全局性，在構(gòu)造平緩區(qū)域，反距離加權(quán)的影響權(quán)重也具有統(tǒng)計學(xué)意義，可以較好地解決速度反演的橫向連續(xù)性問題。影響權(quán)重和半徑的數(shù)學(xué)表達為[31-33]：

式中R表示影響域的半徑，m；dij表示兩道的距離，m；i表示RMS垂直函數(shù)的序號；j表示速度趨勢函數(shù)的序號；表示某給定的值，ξ表示預(yù)設(shè)的極小值，一般為10－2。

加入影響半徑和權(quán)重后趨勢函數(shù)變?yōu)閇32-33]：

式中，內(nèi)層求和表示對第i個垂直速度函數(shù)的所有Ni+1個縱向網(wǎng)格點的累加，外層求和表示橫向網(wǎng)格點i為中心影響半徑內(nèi)的全部Mj個垂直速度函數(shù)的累加。wVn表示影響半徑內(nèi)各道對中心點的影響權(quán)重系數(shù)。

1.1 構(gòu)造約束的趨勢函數(shù)

由式（2）、（3）可知，影響半徑和權(quán)重是保證EAB 函數(shù)計算穩(wěn)定性和合理性的重要因子。在構(gòu)造平緩區(qū)域，速度模型橫向變化平緩，固定影響半徑的引入能提高計算穩(wěn)定性，使單道反演擴展為半徑范圍內(nèi)的多點信息聯(lián)合反演，起到很好的約束效果。但在構(gòu)造復(fù)雜區(qū)域，保證速度橫向連續(xù)的基礎(chǔ)上，要剔除異常高頻振蕩，還要保留地質(zhì)體引起的速度突變和倒轉(zhuǎn)，單個固定的影響半徑存在較大局限性，使用大尺度影響半徑滿足平緩構(gòu)造剔除速度異常的需求，在構(gòu)造變化位置速度模型的精度受損；采用小尺度影響半徑保證構(gòu)造部位的速度精度，平緩構(gòu)造的速度異常未被有效剔除，并且在構(gòu)造復(fù)雜區(qū)域，采用反距離加權(quán)函數(shù)作為影響權(quán)重，很可能為當(dāng)前分析點引入較大誤差，以上問題均會影響成像效果，如圖1所示。

圖1 使用不同影響半徑反演的速度模型的單道速度曲線和偏移成像道集圖

經(jīng)過研究，利用時間偏移剖面的道間相關(guān)系數(shù)，反映地質(zhì)構(gòu)造的橫向連續(xù)和變化特征，在反演方程中自動形成影響半徑，使用自動形成的不同影響半徑匹配不同局部構(gòu)造的規(guī)模和形態(tài)進行反演約束，即構(gòu)造約束的速度趨勢函數(shù)，能提高速度約束反演精度[34-35]，盡可能降低橫向速度振蕩，并提高層速度場對小規(guī)模地質(zhì)體的響應(yīng)。求取道間相關(guān)系數(shù)的方法很多，如歸一化互相算法、基于幾何結(jié)構(gòu)張量的相關(guān)算法、多屬性復(fù)相關(guān)算法等。在山地復(fù)雜地質(zhì)構(gòu)造區(qū)域，資料信噪比常不理想，過于簡單的相關(guān)算法很難達到理想效果，過于復(fù)雜的算法在硬件成本、人力成本和時間成本上很難滿足生產(chǎn)需求，并且資料信噪比低易導(dǎo)致求解困難。因此，通過對比研究，選擇對多道互相關(guān)算法進行改進，分析道采用模型道替代，根據(jù)目的層形成時窗，在時窗內(nèi)進行相關(guān)系數(shù)的計算并取最優(yōu)系數(shù)，保證相關(guān)系數(shù)計算的穩(wěn)定性，公式為 ：

式中ei表示兩地震道的相關(guān)系數(shù)，無量綱；trji表示構(gòu)建當(dāng)前道的模型道的地震道數(shù)據(jù)；tr2i表示除當(dāng)前道的其他地震道數(shù)據(jù)；表示時窗內(nèi)的地震道數(shù)據(jù)的平均值；m表示時窗內(nèi)的縱向樣點數(shù)，n表示納入模型道計算的地震道數(shù)。

由前面分析可知，當(dāng)ei大于預(yù)設(shè)的接受系數(shù)，說明道間相關(guān)系數(shù)大，地震道地質(zhì)特征相似，可累加該處速度信息進行反演，此時影響半徑為納入模型道計算的道數(shù)的累加。因此，構(gòu)造約束的速度趨勢函數(shù)的目標(biāo)方程可寫為：

式（5）是一個最優(yōu)化問題，牛頓法和擬牛頓法是求解最優(yōu)化問題最常用的數(shù)值方法，由于在計算中不需要精確地計算Hessian矩陣，又因L-BFGS（Limited-memory BFGS）算法內(nèi)存占用少，計算復(fù)雜度低，穩(wěn)定性較好，所以本文一并使用L-BFGS算法來求解最優(yōu)值。

1.2 保持構(gòu)造信息的速度體重構(gòu)算法

上述計算是在給定空間拾取點上沿著時間方向離散進行的，數(shù)據(jù)離散化可以減少計算成本、提高計算速度，并且簡化邏輯，降低模型過擬合的風(fēng)險，同時離散化后的特征對異常數(shù)據(jù)也有很強的魯棒性。但初始層速度反演，最終還是要面臨將離散的速度樣點重建為速度體，以供后續(xù)速度更新和偏移成像步驟使用的問題。常規(guī)方法對此討論不多，常使用普通的插值和平滑函數(shù)解決此問題，但無論是一維還是多維插值平滑多基于水平和垂直方向。在構(gòu)造平緩區(qū)域普通的插值和平滑函數(shù)不僅能降低反演求解的多解性，保證反演算法的穩(wěn)定，還能剔除速度模型的異常。但在構(gòu)造復(fù)雜區(qū)域，特別是大傾角構(gòu)造區(qū)域和局部小規(guī)模異常地質(zhì)體區(qū)域，常規(guī)插值平滑不僅會模糊掉小尺度地質(zhì)體和構(gòu)造倒轉(zhuǎn)處的速度響應(yīng)信息，甚至?xí)硭俣饶Ｐ偷呐で?。并且，在山地?fù)雜構(gòu)造區(qū)域因資料信噪比等原因，速度拾取點的空間展布并不一定是均勻給出，這對速度模型網(wǎng)格連續(xù)化帶來更大挑戰(zhàn)。利用徑向基函數(shù)薄板理論形成保持構(gòu)造信息的三維空間速度體重構(gòu)算法可有效解決此問題。

徑向基函數(shù)具有形式簡單，對支撐域和邊界無過多要求，僅需以分析點為中心的子域并集覆蓋整體域就能得到精確結(jié)果。保持構(gòu)造信息的速度體重構(gòu)算法將分布在空間上的節(jié)點構(gòu)造徑向基函數(shù)， 將具有相同時間的層速度在物理模型上類似地形成一個帶有自由邊界的平面薄板，在此平面薄板上的速度信息重構(gòu)可以看作是下面問題s(x,y) 的求解。

薄板樣條（TPS）函數(shù)[36]是離散數(shù)據(jù)插值方法中插值精度非常高的方法[37-38]。Duchon從約束優(yōu)化的角度出發(fā)，以泛函的觀點提出薄板樣條（Thin-plate spline）理論，由其滿足插值條件，并取彎曲能最小的結(jié)論可知[38]，滿足這一問題的s(x,y)具有下面的形式 ：

結(jié)合式（7）、（8），并做下面記號，可以得到與式（7）、（8）等價的方程組（10）：

其中：

圖2給出了一個單層保持構(gòu)造信息的速度體重構(gòu)算法例子，圖 2-a 給出了空間離散節(jié)點的分布及約束條件，插值節(jié)點數(shù)為 1 466 個，空間分布非均勻，約束條件的變化趨勢梯度較大。圖2-b、c顯示了插值結(jié)果及與插值條件的疊合比較，從圖2-b、c中可以看出插值結(jié)果很好地反映了插值條件的變化趨勢，速度模型與插值條件吻合，同時符合地質(zhì)規(guī)律的變化，結(jié)果合理平滑精細(xì)，這正是后續(xù)深度域速度層析反演所需要的。

圖2 構(gòu)造約束的空間三維體插值平滑圖

1.3 復(fù)雜構(gòu)造約束層速度反演方法

山地復(fù)雜構(gòu)造區(qū)的約束層速度反演改進算法從算法分解來看，求解公式（5）得到各網(wǎng)格點的趨勢速度，將其帶入指數(shù)漸進趨勢函數(shù)公式中[31-33]：

形成構(gòu)造約束的速度趨勢函數(shù)擬合：

2 理論模型應(yīng)用效果

采用合成理論模型對本文提出的初始層速度改進反演方法（以下簡稱本文方法）進行應(yīng)用。

理論模型模擬速度倒轉(zhuǎn)情況縱向長度為1 s，橫向長度為1 200 m，共設(shè)置了四層速度，第一層速度值為1 500 m/s，第二層速度值為3 100 m/s，第三層速度為2 000 m/s，第四層速度為4 000 m/s，其中第二層速度和第三層速度發(fā)生倒轉(zhuǎn)，為重點測試本文形成的自適應(yīng)速度趨勢擬合算法和保持構(gòu)造信息的速度體重構(gòu)算法的穩(wěn)定性，在此測試中采用小阻尼權(quán)，減少影響。

圖3-a～d對比可知，本文方法使用了多種參數(shù)都能較好的還原理論模型，模型的形態(tài)和速度都吻合較好，特別是在速度倒轉(zhuǎn)時表現(xiàn)了較好的穩(wěn)定性，速度倒轉(zhuǎn)層不僅速度誤差小而且未出現(xiàn)任何異常高頻，從圖3-e單道曲線可知，速度模型的形態(tài)、速度和深度非常接近真實模型，在速度倒轉(zhuǎn)處最大誤差小于3%，除去速度轉(zhuǎn)換、正則化、網(wǎng)格化、反演迭代平滑等精度損失，本文方法很好地還原了真實模型，證明方法的正確性和較高精度。

圖3 理論模型一和不同參數(shù)的反演速度模型圖

理論模型二模擬了山地復(fù)雜地質(zhì)情況，不僅存在地層尖滅，還存在小尺度局部地質(zhì)體。在速度模型上表現(xiàn)為橫向速度突變、速度倒轉(zhuǎn)和局部高頻速度團，這是速度反演的大挑戰(zhàn)。模型縱向長度為1 s，橫向長度為1 200 m，共設(shè)置了六層速度，分別是1 500 m/s、1 800 m/s、2 000 m/s、2 800 m/s、3 200 m/s、4 000 m/s，在第二層設(shè)置了地層尖滅引起的橫向速度突變，在第三層設(shè)置了小尺度地質(zhì)體引起的速度倒轉(zhuǎn)，第四層設(shè)置了小尺度地質(zhì)體引起的局部速度高頻。采用了不同的參數(shù)測試本文方法在復(fù)雜地質(zhì)環(huán)境下的可靠性和穩(wěn)定性。

由圖 4-a ～ c 對比可知，對于復(fù)雜地質(zhì)速度環(huán)境，本文方法能很好地還原理論模型，從模型的形態(tài)、速度、深度均能與理論模型較好吻合，特別是在地層尖滅和小目標(biāo)地質(zhì)體的速度反演上表現(xiàn)出較高的穩(wěn)定性。僅從此理論模型看，本文基于構(gòu)造的速度趨勢函數(shù)算法與三維復(fù)相關(guān)自適應(yīng)算法[35]的效果類似。從反演速度模型的單道曲線也可以看出除去速度轉(zhuǎn)換、正則化、網(wǎng)格化、反演迭代平滑等帶來的精度損失，速度曲線的形態(tài)、速度和深度非常接近真實模型，證明本方法的正確性和在復(fù)雜地質(zhì)和速度條件下的較高精度。

圖4 理論模型二和不同參數(shù)的反演速度模型和單道速度曲線圖

3 實際數(shù)據(jù)的應(yīng)用效果

為進一步驗證本文方法在實際地震數(shù)據(jù)上的反演效果，將本方法應(yīng)用于深丘地區(qū)弧形褶皺帶交匯處復(fù)雜構(gòu)造區(qū)的地震資料處理上，并通過反演速度模型和偏移成像效果與常規(guī)速度約束反演方法進行對比。實際資料剖面長度16 km，縱向長度8.0 s，給定初始的均方根速度場空間分析點共32道，反演得到的速度模型網(wǎng)格重建為大小650 m×1 250 m。由于實際資料的信噪比不高，為保證疊前深度偏移效果，常規(guī)方法選擇1 km的固定影響半徑。采用兩種方法反演得到的初始深度域?qū)铀俣饶Ｐ鸵妶D5。

圖5 反演得到的初始層速度模型圖

圖5可知，常規(guī)方法采用了較大尺度平滑，但地質(zhì)信息也被模糊掉。采用本文方法得到的速度模型較好地保留了地質(zhì)信息，模型速度變化趨勢與地質(zhì)認(rèn)識一致。

山地構(gòu)造復(fù)雜區(qū)，地層傾角陡，地層速度變化大，如果速度模型精度不高很容易導(dǎo)致構(gòu)造翼部成像不準(zhǔn)或者在逆掩構(gòu)造頂部成像效果差，甚至不成像。疊前深度偏移成像質(zhì)量是檢驗速度模型質(zhì)量的有效手段，將兩套速度場利用積分法進行疊前深度偏移，得到偏移成像結(jié)果。

從圖6、7實際資料的疊前深度偏移結(jié)果可知，本文方法形成的速度模型成像效果更佳，共成像點道集在構(gòu)造部位同相軸更清晰，偏移剖面構(gòu)造成像不論在構(gòu)造翼部還是構(gòu)造頂部成像都更準(zhǔn)確，細(xì)節(jié)更清晰，與地質(zhì)規(guī)律更一致；后續(xù)深度域速度更新時也僅需幾次迭代更新，就可以完成最終深度域速度模型的建立。

圖6 深度偏移部分成像道集對比圖

圖7 疊前深度深度成像剖面對比圖

4 結(jié)論

1）本文在對常規(guī)約束速度反演方法充分研究的基礎(chǔ)上，利用疊前時間偏移剖面上的構(gòu)造信息使速度趨勢函數(shù)求解的影響半徑和權(quán)重隨構(gòu)造特征變化，并利用徑向基函數(shù)薄板理論形成保持構(gòu)造信息的三維空間速度體重構(gòu)算法，最終形成適用于復(fù)雜地質(zhì)區(qū)域的初始層速度改進約束反演方程組。

2）在構(gòu)造復(fù)雜區(qū)，本文方法在改善速度模型縱橫向連續(xù)性的同時，更精細(xì)地描述地質(zhì)特征的速度結(jié)構(gòu)，更好保留小規(guī)模局部地質(zhì)體的速度特征。經(jīng)理論模型測試和實際資料應(yīng)用，本文方法在山地復(fù)雜地質(zhì)區(qū)域有效，且精度較高，得到的層速度模型更有利于深度域速度建模的后續(xù)步驟。

3）本文方法中地震數(shù)據(jù)體道間相關(guān)性的計算非常關(guān)鍵，如果資料信噪比極低，相關(guān)性的計算精度會受到影響。