夏紅敏，鄭小文，黃捍東，肖炎輝，董金超，廖俊

(1.中石化石油勘探開(kāi)發(fā)研究院，北京 100083； 2.中海油研究總院，北京 100028；3.中國(guó)石油大學(xué)(北京) 地球物理學(xué)院，北京 102249)

0 引言

低幅度構(gòu)造是指構(gòu)造比較平緩、閉合幅度較低的一類(lèi)地質(zhì)體，一般構(gòu)造幅度10 m左右。含油后，這種低幅度構(gòu)造及圈閉的幅度小，面積也小，但有可能形成富集高產(chǎn)的“小而肥”油藏，因此成為勘探的重要方向。近年來(lái)，隨著勘探程度的逐漸提高，規(guī)模小、構(gòu)造幅度低且被斷層復(fù)雜化的油氣藏所占的比例逐年上升，儲(chǔ)量品位日趨下降，油田可采儲(chǔ)量與采出資源量之間的矛盾日益尖銳，低幅度構(gòu)造的勘探逐漸受到重視[1-5]。

目前已有一些針對(duì)低幅度構(gòu)造的識(shí)別方法及成圖技術(shù)，如均值濾波法和變速成圖法等，但是他們都有些不足，要么無(wú)法準(zhǔn)確地建立速度變化大的三維速度場(chǎng)，導(dǎo)致畸變、計(jì)算量大、成本高，要么過(guò)分依賴(lài)于構(gòu)造背景數(shù)據(jù)，從而出現(xiàn)低幅構(gòu)造“假象”，影響低幅構(gòu)造識(shí)別效果。傳統(tǒng)方法低幅度構(gòu)造識(shí)別難度大，主要是構(gòu)造背景頻段多，影響復(fù)雜低幅度構(gòu)造的判斷。本文采用去低頻化的方法，提出利用趨勢(shì)分解法和小波變換法對(duì)構(gòu)造走時(shí)數(shù)據(jù)進(jìn)行細(xì)節(jié)與趨勢(shì)(或低頻背景)分解，使斜坡背景或復(fù)雜構(gòu)造背景趨于水平化。從而可以弱化區(qū)域構(gòu)造形態(tài)的影響，消除不同構(gòu)造背景對(duì)低幅構(gòu)造解釋的干擾，突出局部構(gòu)造的細(xì)節(jié)變化，使局部構(gòu)造的形態(tài)更為直觀、清晰[6-15]。

1 基本原理

本文利用兩類(lèi)技術(shù)從不同的角度對(duì)構(gòu)造解釋數(shù)據(jù)進(jìn)行分解，進(jìn)而將構(gòu)造解釋數(shù)據(jù)分為構(gòu)造大背景和低幅構(gòu)造細(xì)節(jié)，最終削弱構(gòu)造大背景影響，增強(qiáng)局部微構(gòu)造特征，使低幅構(gòu)造得以彰顯。這兩種技術(shù)手段分別是對(duì)三維地震構(gòu)造解釋數(shù)據(jù)進(jìn)行處理的趨勢(shì)分解構(gòu)造成圖技術(shù)和對(duì)二維構(gòu)造平面數(shù)據(jù)進(jìn)行處理的小波分解構(gòu)造成圖技術(shù)，圖1為技術(shù)路線(xiàn)圖。

圖1 低幅度構(gòu)造識(shí)別技術(shù)路線(xiàn)Fig.1 Roadmap of low-amplitude structure recognition technology

1.1 趨勢(shì)分解法原理

趨勢(shì)分解法是在三維地震解釋構(gòu)造走時(shí)數(shù)據(jù)體上，根據(jù)構(gòu)造上下起伏及走勢(shì)搜索拐點(diǎn)位置，再根據(jù)拐點(diǎn)位置的構(gòu)造走時(shí)信息插值得到構(gòu)造趨勢(shì)面，最后用原來(lái)的構(gòu)造數(shù)據(jù)減去構(gòu)造趨勢(shì)面，即可獲得低幅度構(gòu)造的細(xì)節(jié)。趨勢(shì)分解構(gòu)造成圖技術(shù)主要包括4個(gè)步驟：①趨勢(shì)化預(yù)處理，即對(duì)構(gòu)造層位數(shù)據(jù)預(yù)處理，濾掉隨機(jī)因子對(duì)趨勢(shì)面的影響，有效控制構(gòu)造趨勢(shì)面；②構(gòu)造起伏拐點(diǎn)搜索，構(gòu)造起伏拐點(diǎn)是指構(gòu)造由逐漸下降轉(zhuǎn)為逐漸上升的轉(zhuǎn)折點(diǎn)，往往對(duì)應(yīng)于低幅構(gòu)造的圈閉溢出點(diǎn)，是勾勒低頻構(gòu)造背景的關(guān)鍵控制點(diǎn)，在成圖過(guò)程中可以根據(jù)構(gòu)造的幅度自適應(yīng)調(diào)整搜索拐點(diǎn)半徑尋找局部極值點(diǎn)來(lái)確定拐點(diǎn)；③構(gòu)造低頻背景勾勒，針對(duì)構(gòu)造起伏拐點(diǎn)，利用反距離加權(quán)法得到低頻背景；④去低頻化構(gòu)造成圖，將構(gòu)造解釋與低頻背景疊合，提取低幅構(gòu)造數(shù)據(jù)，形成低幅度構(gòu)造平面圖。

從圖2中可以看出，構(gòu)造起伏拐點(diǎn)的搜索及搜索半徑對(duì)構(gòu)造勾勒至關(guān)重要，它既控制低頻構(gòu)造背景，又通過(guò)自適應(yīng)調(diào)整搜索拐點(diǎn)半徑控制拐點(diǎn)的數(shù)量，使構(gòu)造的低頻化處理更加合理、適當(dāng)。圖3是地震剖面中的構(gòu)造起伏拐點(diǎn)位置及搜索半徑。在實(shí)現(xiàn)過(guò)程中，基于低幅構(gòu)造大小自適應(yīng)調(diào)整搜索拐點(diǎn)半徑。利用反距離加權(quán)法插值的方法勾勒低頻構(gòu)造背景，可保證低頻構(gòu)造背景的平滑和準(zhǔn)確性。

圖2 構(gòu)造拐點(diǎn)示意Fig.2 Schematic diagram of construction inflection point

圖3 構(gòu)造起伏拐點(diǎn)及搜索半徑示意Fig.3 Construction of undulation inflection point and search radius map

通過(guò)趨勢(shì)分解法進(jìn)行構(gòu)造成圖，最大的優(yōu)勢(shì)就在于不用構(gòu)造速度場(chǎng)，通過(guò)拐點(diǎn)的有效搜索和確定，使得構(gòu)造解釋既包含低幅度構(gòu)造趨勢(shì)又突出局部微小起伏。

1.2 小波分解法原理

小波變換是一種多尺度變換分析方法，它繼承和發(fā)展了短時(shí)傅里葉變換局部化的思想，同時(shí)又克服了窗口大小不隨頻率變化等缺點(diǎn)，能夠提供一個(gè)隨頻率改變的“時(shí)間—頻率”窗口，它的特點(diǎn)是通過(guò)變換能夠充分突出問(wèn)題某些方面的特征，能對(duì)時(shí)間(空間)頻率的局部化分析，通過(guò)伸縮平移運(yùn)算對(duì)信號(hào)(函數(shù))逐步進(jìn)行多尺度細(xì)化，最終達(dá)到高頻處頻率細(xì)分，低頻處頻率粗分，能自動(dòng)適應(yīng)時(shí)頻信號(hào)分析的要求，從而可聚焦到信號(hào)的任意細(xì)節(jié)，進(jìn)而獲得低頻構(gòu)造背景[6-7]。

在本文中，我們用到的二維小波正變換公式為：

(1)

二維構(gòu)造走時(shí)解釋平面數(shù)據(jù)經(jīng)過(guò)二維小波多尺度分析后，選取大尺度數(shù)據(jù)進(jìn)行二維小波逆變換即得到構(gòu)造大背景。其二維小波逆變換的公式為：

(2)

(3)

其中：Cg為小波函數(shù)的傅里葉變換，Gg為g的傅里葉變換。

小波分頻進(jìn)行構(gòu)造識(shí)別是通過(guò)人工選取合適的尺度參數(shù)，利用二維小波變換技術(shù)對(duì)構(gòu)造解釋平面數(shù)據(jù)進(jìn)行多尺度分析獲得低頻(多尺度)構(gòu)造背景，進(jìn)而將背景與局部細(xì)節(jié)(低幅度構(gòu)造)分離，達(dá)到凸顯低幅度細(xì)節(jié)的目的。充分利用小波分解的時(shí)頻特征，對(duì)不同頻率成分進(jìn)行分解和突出，這樣將整個(gè)構(gòu)造圖聚焦到微小起伏等細(xì)節(jié)上來(lái)，得到能凸顯細(xì)節(jié)變化的去低頻化低幅度構(gòu)造圖[16-19]。

1.3 去低頻化構(gòu)造成圖

利用去低頻化思路，采用趨勢(shì)分解法和小波分頻法進(jìn)行低幅度構(gòu)造成圖。值得注意的是，在趨勢(shì)分解法實(shí)現(xiàn)過(guò)程中，先進(jìn)行趨勢(shì)化預(yù)處理，濾掉隨機(jī)因子干擾趨勢(shì)面的影響。然后進(jìn)行構(gòu)造低頻背景勾勒，通過(guò)搜索構(gòu)造起伏拐點(diǎn)，利用反距離加權(quán)法得到低頻背景。隨后在構(gòu)造大背景預(yù)測(cè)圖中，重新調(diào)整成圖參數(shù)，進(jìn)一步提高構(gòu)造的大尺度連續(xù)性。為了提高構(gòu)造成圖的細(xì)節(jié)變化，可利用小波分頻的原理。優(yōu)選出小波基函數(shù)，對(duì)構(gòu)造數(shù)據(jù)進(jìn)行二維小波變換，得到二維尺度時(shí)頻圖。再通過(guò)小波逆變換，得到基于小波分解的低幅度構(gòu)造平面圖。

圖4為兩種方法的局部對(duì)比，從圖中可知，兩種方法體現(xiàn)構(gòu)造宏觀趨勢(shì)一致，只是局部有差異，因此在實(shí)現(xiàn)過(guò)程中，將二者有效結(jié)合，利用這兩種方法，可以消除大構(gòu)造背景影響，展現(xiàn)低幅構(gòu)造細(xì)節(jié)。

a—趨勢(shì)分解低頻背景;b—小波分解低頻背景;c—趨勢(shì)分解低幅構(gòu)造細(xì)節(jié);d—小波分解低幅構(gòu)造細(xì)節(jié)a—low frequency background of trend decomposition;b—low frequency background of wavelet decomposition;c—low amplitude strunctural details of trend decomposition;d—low amplitude strunctural details of wavelet decomposition圖4 兩種構(gòu)造成圖法對(duì)比Fig.4 Comparison map of two construction methods

2 實(shí)例應(yīng)用及效果分析

本次選取的研究區(qū)位于E盆地，目的層為白堊系M組砂巖油層，主要發(fā)育低幅度構(gòu)造和構(gòu)造—巖性復(fù)合圈閉，如圖5所示，其構(gòu)造幅度不超過(guò)10～15 m，由于構(gòu)造幅度小于地質(zhì)構(gòu)造解釋精度，加上儲(chǔ)層薄的特點(diǎn)，常規(guī)構(gòu)造方法難以精細(xì)刻畫(huà)低幅度構(gòu)造圈閉中的油氣分布。

圖5 E盆地北部示意Fig.5 Slope zone of a basin in Ecuador

經(jīng)過(guò)上述去低頻化的構(gòu)造分析，利用趨勢(shì)分解和小波分解法消除大構(gòu)造背景影響，展現(xiàn)低幅構(gòu)造細(xì)節(jié)。兩種方法對(duì)比，各有千秋。趨勢(shì)分解法在處理不規(guī)則或有缺失的數(shù)據(jù)時(shí)有優(yōu)勢(shì)，但成圖的效果依賴(lài)于搜索半徑，稍遜色于尺度連續(xù)的小波變換；二維小波分解法則需要數(shù)據(jù)完整，否則需要進(jìn)行數(shù)據(jù)外推，一定程度上影響低幅構(gòu)造識(shí)別精度。圖6、7顯示的是低幅度構(gòu)造解釋分析的局部效果，并將平面圖與剖面圖對(duì)比分析，從地震剖面中，可以清晰識(shí)別4個(gè)低幅構(gòu)造(圖中標(biāo)注①、②、③、④)，其在平面上也可清晰展現(xiàn)，平面、剖面對(duì)應(yīng)良好，成圖效果理想。

圖6 趨勢(shì)分解法細(xì)節(jié)構(gòu)造Fig.6 Detail construction map of trend decomposition method

圖8為研究區(qū)M1層低幅構(gòu)造平面對(duì)比，從圖中不難看出：兩種方法預(yù)測(cè)結(jié)果相似性較大，趨勢(shì)分解法在工區(qū)西部，幅度較小的低幅構(gòu)造識(shí)別上效果更精細(xì)，小波分解法可以輔助低幅構(gòu)造的解釋。但趨勢(shì)分解法成圖效果依賴(lài)搜索半徑，使得構(gòu)造特征不連續(xù)，稍遜色于尺度連續(xù)的小波變換。經(jīng)過(guò)小波變換后的大尺度時(shí)頻分析，保持住構(gòu)造低頻背景，從而保證構(gòu)造連續(xù)性。既有構(gòu)造背景趨勢(shì)，又具有局部聚焦性。整體上看，利用小波分解法消除了大構(gòu)造背景影響，展現(xiàn)了低幅構(gòu)造細(xì)節(jié)。將兩者有機(jī)結(jié)合，通過(guò)與單一方法對(duì)比，結(jié)合后得到的構(gòu)造圖不僅消除了大構(gòu)造背景的影響，而且可以展現(xiàn)低幅構(gòu)造細(xì)節(jié)?？梢杂行?shí)現(xiàn)對(duì)研究區(qū)低幅度構(gòu)造的識(shí)別和構(gòu)造成圖分析。

圖7 小波分解法細(xì)節(jié)構(gòu)造Fig.7 Detail construction map of wavelet decomposition

a—M1層原始局部構(gòu)造平面;b—基于趨勢(shì)分解法M1層局部構(gòu)造平面;c—基于小波分解法M1層局部構(gòu)造平面a—original local tectonic plane of M1 layer;b—local tectonic plane of M1 layer based on trend decomposition method;c—local tectonic plane of M1 layer based on wavelet decomposition method圖8 研究區(qū)M1層低幅度局部構(gòu)造平面對(duì)比Fig.8 Contrast map of low-amplitude structural plane of M1 layer in the study area

3 結(jié)論

本文利用趨勢(shì)分解法和小波分解法聯(lián)合的低幅度構(gòu)造識(shí)別及成圖方法，既增強(qiáng)了局部微構(gòu)造特征，又兼顧了低幅度構(gòu)造宏觀背景。在E盆地北部中的應(yīng)用中效果明顯，給后續(xù)的地質(zhì)解釋和儲(chǔ)層描述提供了依據(jù)。通過(guò)分析研究，得出以下結(jié)論：

1)由于低幅構(gòu)造尺度大小不一，本文充分利用二維小波變換尺度的冗余特性，實(shí)現(xiàn)了對(duì)構(gòu)造解釋成果的二維圖像多尺度分解；

2)為避免小波變換在構(gòu)造解釋數(shù)據(jù)缺失的情況下，邊緣效應(yīng)嚴(yán)重的問(wèn)題，本文使用趨勢(shì)分解法，實(shí)現(xiàn)了去低頻化的多尺度處理，并在較大程度上克服了因數(shù)據(jù)點(diǎn)缺失帶來(lái)的邊緣效應(yīng)；

3)針對(duì)目前低幅構(gòu)造成圖技術(shù)的不足，本文從不同角度采用趨勢(shì)分解法與二維小波分解法分別對(duì)低幅構(gòu)造進(jìn)行成圖，實(shí)現(xiàn)區(qū)域背景與局部低幅度構(gòu)造的有效分離，增強(qiáng)局部微構(gòu)造特征，使其在構(gòu)造平面圖上得以凸顯。