康鯤鵬

（商丘師范學(xué)院計算機科學(xué)系，河南商丘476000）

時移地震是利用多次采集的時移地震資料，結(jié)合巖石物理和生產(chǎn)動態(tài)數(shù)據(jù)，通過提取和綜合分析各種敏感屬性，分析解釋油氣藏動態(tài)，并對油藏進行模擬處理分析，以確定剩余油的分布范圍，指導(dǎo)開發(fā)井的部署和調(diào)整，提高采收率。因此，時移地震方法廣泛應(yīng)用于油藏監(jiān)測、基礎(chǔ)測量和監(jiān)測測量。流動單元是從宏觀到微觀的不同級次上的、空間三維連續(xù)分布的、具有相似的影響流體流動的巖石特征和流體本身滲流特征的儲集巖體；同一流動單元內(nèi)各處相對均質(zhì)，不同流動單元之間儲、滲能力差異顯著，且有較好的隔擋界面或滲流屏障。其發(fā)育特征和空間分布狀況受沉積作用、構(gòu)造作用和成巖作用等的綜合控制。對流動單元的研究有很多種方法，如精細沉積學(xué)研究方法、FZI劃分法等，近年來利用時移地震資料進行流動單元解釋也得到發(fā)展。

1 國內(nèi)針對流動單元的時移地震資料解釋現(xiàn)狀

時移地震技術(shù)是一門綜合地質(zhì)勘探技術(shù)，包含諸如巖石物理技術(shù)、采集技術(shù)和處理解釋技術(shù)等，目前，我國在這方面的理論研究比較多，王大偉等[1]通過對地震資料的解釋分析了油層內(nèi)部滲流的影響，為劃分油藏流體流動單元提供依據(jù)。桑淑云[2]針對SZ油田，從時移地震差異合理性分析、地震敏感屬性、井上動態(tài)分析、流動單元劃分、時移地震約束下的油藏數(shù)值模擬等方面著手，經(jīng)過分析對比，在地震可分辨的基礎(chǔ)上對油田主體部分的小層的水驅(qū)效果及剩余油分布進行了分析；他們在分析過程中用到的軟件多為國外專業(yè)或者是類似MATLAB這樣的通用軟件。國內(nèi)基于對時移地震資料解釋分析的軟件很少，而對流動單元的時移地震資料解釋軟件研究目前更是處于起步階段。

2 時移地震技術(shù)分析

2.1 巖石物理技術(shù)

對于時移地震技術(shù)而言，巖石物理技術(shù)是連接地震響應(yīng)和油藏參數(shù)（如壓力、流體飽和度等）的橋梁。人們注意到近幾年來實際的時移地震響應(yīng)變化往往大于實驗室中觀測到的由流體以及壓力變化引起的響應(yīng)變化。另外，通過油藏模型的正演發(fā)現(xiàn)，當(dāng)油藏的厚度和調(diào)諧厚度相近時引起的變化比速度的變化大得多。這主要是因為速度的變化導(dǎo)致調(diào)諧位置的變化，這與地震的主頻有關(guān)。當(dāng)兩次地震響應(yīng)相減時，差異就比較大。因此，在進行時移地震的可行性分析時，必須考慮油藏本身的變化，以及油藏的其他參數(shù)，如厚度和地震主頻。時移地震可行性研究綜合考慮了巖石物理、地震資料的品質(zhì)及油藏的若干信息和條件。在風(fēng)險圖中劃分出高、中、低風(fēng)險區(qū)。

2.2 采集技術(shù)

近幾年來，時移地震技術(shù)多應(yīng)用在重復(fù)采集的三維區(qū)中，采集技術(shù)方面的進展相對緩慢。通常把在同一地區(qū)不同時期采集的數(shù)據(jù)稱為繼承性數(shù)據(jù)（Legacy Data）。這種數(shù)據(jù)往往在采集時沒有考慮到時移地震技術(shù)的應(yīng)用，成本低。另外一種采集方式為目的性重復(fù)采集，在已有的三維區(qū)為了時移地震技術(shù)的應(yīng)用進行另一個時間的采集，這種數(shù)據(jù)的成本次之，但其重復(fù)性相對比較好。如果把檢波器永久固定在同一個位置，在不同時間進行地震波激發(fā)，稱之為永久性檢波器的采集，其成本最高。

在采集設(shè)計上，目的性重復(fù)采集主要要求檢波器以及炮點盡量保持一致。GPS技術(shù)的發(fā)展使這種一致性比較容易達到。在采集方面，Western Geco近年推出的Q采集系統(tǒng)采用高密度網(wǎng)格數(shù)據(jù)接收，可將響應(yīng)的重復(fù)性和位置的重復(fù)性問題變成相同位置取舍檢波器數(shù)據(jù)進行處理的問題。

2.3 處理技術(shù)

均衡處理技術(shù)是在完成常規(guī)處理后為了提高可重復(fù)性，對不同時間的地震響應(yīng)進行匹配，使油藏以外的差異性變小，而油藏位置上的差異性比較大。具有相同面元的數(shù)據(jù)體就可以進行均衡處理，該處理有時很難消除處理過程中的差異，很多情況下要對原始數(shù)據(jù)體重復(fù)進行處理。這種重新處理的過程一般采用常規(guī)的處理流程就可完成，但其處理參數(shù)和原則與三維處理方法不同，主要包括：

1）噪聲的獨立消除，即不同數(shù)據(jù)體的噪聲由于采集條件不同，去噪必須是相對獨立的；2）相同的疊加速度和偏移距；3）相同的網(wǎng)格生成；4）相同的偏移速度。

3 時移地震流動單元解釋

3.1 解釋中的可視化技術(shù)

可視化即科學(xué)計算可視化，它是計算機應(yīng)用中的一個重要領(lǐng)域，可把函數(shù)值計算或?qū)嶒灥玫降拇罅繑?shù)據(jù)表現(xiàn)為人的視覺可以感受的圖像。三維空間數(shù)據(jù)場的顯示是實現(xiàn)科學(xué)計算可視化的核心?？梢暬夹g(shù)是將計算機圖形學(xué)、圖像處理技術(shù)、虛擬現(xiàn)實技術(shù)相結(jié)合，將難以理解的抽象的數(shù)據(jù)場轉(zhuǎn)換為直觀的圖形圖像信息。將N維域中的一組離散的數(shù)據(jù)映射為彩色象素值，并能進行交互處理的理論、方法和技術(shù)。

一般來說，三維空間數(shù)據(jù)場是連續(xù)的，而數(shù)值計算結(jié)果或測量所的數(shù)據(jù)則是離散的，是采集連續(xù)的三維場的結(jié)果。體繪制技術(shù)就是要將這一種三維空間樣本直接轉(zhuǎn)換為屏幕上的二維圖像，盡可能準確地重現(xiàn)原始的三維數(shù)據(jù)場。屏幕上的二維圖像決定于幀緩存中對應(yīng)與每一個象素點的光亮度值，這也是一個二維的離散數(shù)據(jù)場。因此，體繪制技術(shù)的實質(zhì)是將離散的三維空間數(shù)據(jù)場轉(zhuǎn)化為離散的二維數(shù)據(jù)。

將離散的三維空間數(shù)據(jù)場轉(zhuǎn)化為離散的二維數(shù)據(jù)點陣，首先必須進行三維空間數(shù)據(jù)場的重新采樣。其次，應(yīng)該考慮三維空間中每一個數(shù)據(jù)對二維圖像的貢獻，因而必須實現(xiàn)圖像的合成。所以，體繪制技術(shù)的實現(xiàn)是一個三維離散數(shù)據(jù)場的重新采樣的圖象合成的過程。實現(xiàn)重新采樣從理論上說應(yīng)有以下步驟：

1）選擇適當(dāng)?shù)闹貥?gòu)函數(shù)，對離散的三維數(shù)據(jù)場進行三維卷積運算，重構(gòu)連續(xù)的三維數(shù)據(jù)場；2）對連續(xù)的三維數(shù)據(jù)場根據(jù)給定的觀察方向進行幾何變換；3）由于屏幕上采樣點的分辨率是已知的，由此可計算出被采樣信號的Nyquist頻率極限，采用低通濾波函數(shù)去掉高于這一極限的頻率；4）對濾波后的函數(shù)進行重新采樣。

由于進行三維卷積運算比較費時，因此可采用離散方法實現(xiàn)。上述提到的物體空間為序和圖象空間為序2種不同的體繪制算法，只是實現(xiàn)重新采樣的不同方法。

繪制技術(shù)是科學(xué)計算可視化中的核心技術(shù)。體可視化的實現(xiàn)方法包括面繪制和體繪制2種。其中，面繪制的過程包括等值面轉(zhuǎn)化，可見部分繪制和陰影部分繪制；而體繪制過程包括圖像的可見部分繪制和陰影部分的繪制。體數(shù)據(jù)是一個三維連續(xù)的幾何模型，通常代表一個幾何化或參數(shù)化定義的多邊形或平面顯示列表，這個幾何模型被掃描轉(zhuǎn)化為一組體數(shù)據(jù)集，使體數(shù)據(jù)集可視，則原物可能被直接投影或減至二維象素空間，并且存儲為一個光柵圖形于一個幀緩沖區(qū)中。

體數(shù)據(jù)可能先以一定的方法轉(zhuǎn)化為幾何物體，如等值面處理、等值線處理、面抽象處理或依邊界處理。然后，幾何物體通過常規(guī)幾何繪制投影到屏幕。體繪制是直接將原物顯示而沒有任何體數(shù)據(jù)中介轉(zhuǎn)化到表面表示；體繪制采用帶有透明度的體元微粒模型繪制；它可以再現(xiàn)細節(jié)，獲得高質(zhì)量的圖形。因面體繪制計算量大，對硬件的要求也比較高。面繪制對原物可視化的基本技術(shù)，先將原物轉(zhuǎn)化為一個平面表示，然后用常用的計算機圖形技術(shù)投影平面到屏幕；一般情況下，面繪制采用確定表面的數(shù)據(jù)，且表面可能是實在的并且是可視的；這是一種有損繪制，但是效率比較高。體繪制和面繪制這兩種方法是相互補充的，無論采用體繪制，面繪制還是混合繪制方法，主要依據(jù)特定的可視化任務(wù)。

體繪制-面繪制技術(shù)是可視化實現(xiàn)的重要方法，它可以揭示數(shù)據(jù)場內(nèi)部的細微與難于定義的結(jié)構(gòu)，可以顯示數(shù)據(jù)場的整體特征。但因為其計算量大，與數(shù)據(jù)場的規(guī)模成正比以及隨機存取大量的體數(shù)據(jù)而需占用巨大的帶寬，難于滿足實時交互需要。

可視化技術(shù)可以直接從地震數(shù)據(jù)中看到蘊含地質(zhì)現(xiàn)象和規(guī)律，無需經(jīng)傳統(tǒng)的地質(zhì)解釋就可直觀、全方位地看到地層界面的起伏、斷層的形態(tài)，甚至可直接觀察到沉積體系的空間展布，大大提高地質(zhì)解釋的效率和可靠性。

3.2 可視化算法

三維可視化技術(shù)的思想是用直觀的圖形輸出代替枯燥的數(shù)據(jù)輸出，借助人類強大的視覺及形象思維能力，對數(shù)據(jù)進行本質(zhì)上的理解，從而洞察、發(fā)現(xiàn)數(shù)據(jù)中隱藏的現(xiàn)象和規(guī)律，為獲取深層次信息提供了強有力的手段，同時也極大提高了工作效率?？梢暬暮诵氖侨S數(shù)據(jù)場的可視化[10]，可視化算法可分為2大類。第一類是基于面的體繪制算法，它首先由三維空間數(shù)據(jù)場構(gòu)造出中間幾何圖元，最常見的中間幾何圖元就是等值面。在三維地震數(shù)據(jù)體的面繪制中，最常用的就是進行各種切片顯示。第二類算法并不構(gòu)造中間幾何圖元，而是直接由三維數(shù)據(jù)場產(chǎn)生屏幕上的二維圖像，稱為直接體繪制算法，這種算法能產(chǎn)生三維數(shù)據(jù)場的整體圖像，包括每一個細節(jié)，并具有圖象質(zhì)量高、便于并行處理的優(yōu)點。

三維地震已經(jīng)成為目前油氣勘探開發(fā)的主要方法，三維地震數(shù)據(jù)體中蘊含著地下巖層的豐富信息，以前由于沒有合適的三維顯示設(shè)備，只利用一條條地震剖面以及水平切片來顯示三維數(shù)據(jù)，地質(zhì)解釋人員依據(jù)這些二維圖像來推測、想象地下地層的空間形態(tài)和結(jié)構(gòu)。對于龐大的三維數(shù)據(jù)體，即使將每一條剖面和水平切片都顯示出來，一個三維地質(zhì)體的完整信息也是分散在各個獨立的二維圖像中，很難綜合起來進行聯(lián)想，觀察也不直觀，不能從三維的角度去觀察和分析地質(zhì)體的空間形態(tài)。更何況在實際工作中不可能對每條剖面和切片都進行觀察，常常是只對部分剖面和水平切片進行解釋，只利用了一小部分信息，大部分信息沒得到利用，其結(jié)果是很難了解地層結(jié)構(gòu)的三維細節(jié)，不可避免地漏失了大量的小油氣田。因此，采用三維可視化技術(shù)及虛擬現(xiàn)實技術(shù)實現(xiàn)三維數(shù)據(jù)三維解釋。

對于三維地震數(shù)據(jù)場，目前常用的切片式可視化方法帶有明顯的三維數(shù)據(jù)二維解釋的痕跡，所表達的信息是片面的、孤立的，難以反映原始數(shù)據(jù)的全貌，丟失了寶貴的三維細節(jié)，仍然沒有擺脫二維解釋方法的局限性，本質(zhì)上仍然是一種三維數(shù)據(jù)二維顯示的方法，沒有體現(xiàn)三維數(shù)據(jù)場可視化的真正含義，無法直接看到反射面的形態(tài)以及相互之間的接觸關(guān)系。

與傳統(tǒng)剖面解釋方法完全不同，常規(guī)的三維解釋是通過對每一條地震剖面上的每個層位、每條斷層拾取后，再通過三維空間的組合來完成的。三維體可視化解釋是通過對來自于地下界面的地震反射率數(shù)據(jù)體采用各種不同的透明度參數(shù)在三維空間內(nèi)直接解釋地層的構(gòu)造、巖性及沉積特點。這種三維立體掃描和追蹤技術(shù)可使解釋人員快速選定目標(biāo)，結(jié)合精細的鉆井標(biāo)定，可幫助解釋人員準確快速的描述各種復(fù)雜的地質(zhì)現(xiàn)象。三維可視化是根據(jù)數(shù)據(jù)體的透明度屬性，假定地下界面的反射率是地下界面的原始、真正的三維模型，本質(zhì)上講，它是由三維空間中的構(gòu)造、地層及振幅屬性綜合組成的。無論是做三維區(qū)域分析，還是做特定前景目標(biāo)評價（包括流體界面識別），都可以通過這種‘進去看’的方式快速完成。在基于三維像素的立體可視化中，每個數(shù)據(jù)樣點都被轉(zhuǎn)換成為一個三維像素（其大小近似面元間距和采樣間隔的三維像素）。每一個三維像素具有與原三維數(shù)據(jù)母體相對應(yīng)的數(shù)值，一個三色（紅、綠、藍）值以及一個暗度變量，該變量用來調(diào)整數(shù)據(jù)體的透明度。這樣，每一個地震道被轉(zhuǎn)換成為一個三維像素柱。圖1為多屬性可視化。

圖1 多屬性體可視化Fig.1 Multi-attribute visualization

3.3 時移地震流動單元解釋參數(shù)設(shè)置

本系統(tǒng)的時移地震流動單元解釋中需要對解釋點、線、面的顏色、大小等參數(shù)[3-4]進行恰當(dāng)?shù)脑O(shè)置，以便在形成圖件的時候更利于觀察。參數(shù)設(shè)置對話框如圖2所示。

圖2 時移地震流動單元解釋參數(shù)設(shè)置Fig.2 Flow-unit parameter setting of time-lapse seismic

3 .4時移地震流動單元解釋

圖2所示的流動單元點、線、面都是需要處理的對象，采用以下數(shù)據(jù)結(jié)構(gòu)進行描述，定義流動單元類CQISSegFlowUnit。

1）流動單元定義

class CQISSegFlowUnit

{Public：

CQISSegFlowUnit（）；//構(gòu)造函數(shù)

virtual～CQISSegFlowUnit（）；//析構(gòu)函數(shù)public://定義方法

void Display（）；//流動單元顯示

int InsertVex（StSeedType seed）；//插入解釋點

void DelSeed（bool bMove=true）；//刪除解釋點

void MoveLine（StSeedType seed）；//移動解釋面

private：

bool isSameSeed（StSeedType ArraySeed，StSeedType Seed）；//判斷兩解釋點是否相同

bool checkSam（StSeedType seed）；//檢查解釋點

int InsertVexIn（StSeedType seed）；//在測線剖面上插入解釋點

void InsertMoveLineIn（）；//在測線剖面上插入移動解釋面

private://定義成員

std::list＜std::list＜StSeedType＞＞::iterator pCurLine；//當(dāng)前解釋面

std::list＜StSeedType＞::iterator pCurLineSeed；//當(dāng)前行的解釋點

StSeedType m_Seed；//當(dāng)前選準解釋點

StFlowUnitPara m_FlowUnitPara；//流動單元屬性

}

2）解釋點結(jié)構(gòu)體

struct StSeedType

{int CrosslineNum；//解釋點坐標(biāo)Xline

int TimeNum；//解釋點坐標(biāo)TimeNum

int InlineNum；//解釋點坐標(biāo)InLine

float x；//解釋點顯示坐標(biāo)x

float y；//解釋點顯示坐標(biāo)y

float z；//解釋點顯示坐標(biāo)z

}

3）流動單元屬性結(jié)構(gòu)體定義

struct StFlowUnitPara

{intnSelectLineType；

//流動單元類型1：測線剖面流動單元；

//2：聯(lián)絡(luò)測線剖流動單元；

//3：TimeSlice切片流動單元。

int nLineType；

//流動單元走向：1：測線剖面；2：聯(lián)絡(luò)

int nFlowUnitLineType；

//映射相鄰剖面解釋TRUE：映射；

}

在剖面或切片上解釋流動單元時需要利用增加解釋點的方法實現(xiàn)[5]。首先要判斷增加的解釋點是否在當(dāng)前剖面或者切片上，如果在，就分側(cè)面或剖面上2種情況進行解釋；否則什么也不做，直接返回。具體算法流程圖如圖3所示。

4 應(yīng)用實例

利用本系統(tǒng)定義的數(shù)據(jù)結(jié)構(gòu)在剖面和側(cè)面上解釋用到的流動單元，取得了良好的效果，對時移地震數(shù)據(jù)體“東方-海底電纜-中間成果”主測線（360～500）、聯(lián)絡(luò)測線（800～1 100）、時間（450～750 ms），采用本系統(tǒng)的算法進行流動單元解釋。若從時移地震數(shù)據(jù)體中Section0面開始流動單元解釋，解釋方法如下：

1）把Section0作為引導(dǎo)面，在面上解釋流動單元線。

根據(jù)Section0引導(dǎo)面信息，沿流動單元走向添加解釋點（如圖4（a）），Section0引導(dǎo)面解釋完畢。如果流動單元全部解釋完畢，則解釋結(jié)束，否則轉(zhuǎn)入2）；

2）在數(shù)據(jù)體中拖動Section0到Section1，此時如果：①選擇流動單元解釋信息自動映射，系統(tǒng)自動把Section0上的流動單元面投影到Section1上，在Section0與Section1之間通過插值繪出流動單元體，根據(jù)Section1上同一流動單元的形狀修改Section1上的流動單元解釋點（如圖4（b）～圖4（d））。如果解釋流動單元全部解釋完畢，則解釋結(jié)束，否則以Section1作為新的引導(dǎo)面Section0，轉(zhuǎn)入2）。②選擇流動單元解釋信息不自動映射，把Section1作為Section0，轉(zhuǎn)入1）繼續(xù)。

圖3 流動單元解釋算法流程圖Fig.3 Flow chart of flow-unit algorithm

圖4 時移地震流動單元解釋圖件Fig.4 Flow-unit interpretation maps of time-lapse seismic

5 結(jié)束語

本軟件第二版已開發(fā)完成，主要功能達到設(shè)計要求。軟件已通過檢測，并取得國家版權(quán)局計算機軟件著作權(quán)登記證書。軟件已在多家單位實現(xiàn)銷售，用戶一致認為在Windows環(huán)境下運行的本軟件，技術(shù)新穎，實用性強；軟件易學(xué)易用、易于推廣；軟件界面友好，性能穩(wěn)定。軟件中提供了豐富的油藏綜合分析方法，特別是基于地質(zhì)約束的儲層預(yù)測方法、典型相關(guān)屬性優(yōu)化和儲層預(yù)測方法等更是其他軟件所沒有，有助于提高油藏預(yù)測準確度，減少勘探風(fēng)險，必將在我國的油氣勘探生產(chǎn)中發(fā)揮重要作用。

[1]王大偉，劉震，趙偉，等.利用時移地震資料劃分油藏流體流動單元的可行性分析[J].地球物理學(xué)報，2007，50（2）：46-48.

WANG Da-wei，LIU Zhen，ZHAO Wei，et al.Feasibility analysis on division of flow units using time-lapse seismic data[J].Chinese Journal of Geophysics，2007，50（2）：46-48.

[2]陳小宏，牟永光.四維地震油藏監(jiān)測技術(shù)及其應(yīng)用[J].石油地球物理勘探，1998（6）：198-201.

CHEN Xiao-hong，MOU Yong-guang.Four-dimensional seismic reservoir monitoring technique and its application[J].Oil Geophysical Prospecting，1998（6）：198-201.

[3]康立明，任戰(zhàn)利.多參數(shù)定量研究流動單元的方法——以鄂爾多斯盆地W93井區(qū)為例[J].吉林大學(xué)學(xué)報：地球科學(xué)版，2008，38（5）：282-287.

KANG Li-ming，REN Zhan-li.Study on multi-parameters discrimination method for flow units—taking W93 wellblock in ordos basin as an example[J].Journal of Jilin University：Earth Science Edition，2008，38（5）：282-287.

[4]史成恩，解偉，孫衛(wèi)，等.靖安油田盤古梁長6油藏流動單元的定量劃分[J].石油與天然氣地質(zhì)，2006，27（2）:239-243.

SHI Cheng-en，XIE Wei，SUN Wei.Quantitative division of flow units in Panguliang Chang 6 reservoir，Jing’an oilfield[J].Oil&Gas Geology，2006，27（2）：239-243.

[5]趙偉.中國海上時移地震技術(shù)應(yīng)用的可行性研究[J].勘探地球物理進展，2003，26（1）:33-37.

ZHAO Wei.Feasibility study on time-lapse seismic offshore china[J].Progress in Exploration Geophysics，2003，26（1）33-37.

[6]Smith T M，Sondergeld C H，Chandra S R.Gassmann fluid substitutions:A tutorial[J].Geophysics，2003，68（2）:430-440.