時(shí)移地震流動(dòng)單元自動(dòng)追蹤與解釋成果分析

康鯤鵬

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

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

1 國內(nèi)針對(duì)流動(dòng)單元的時(shí)移地震資料解釋現(xiàn)狀

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

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

2.1 巖石物理技術(shù)

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

2.2 采集技術(shù)

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

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

2.3 處理技術(shù)

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

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

3 時(shí)移地震流動(dòng)單元解釋

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

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

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

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

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

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

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

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

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

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

3.2 可視化算法

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

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

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

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

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

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

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

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

3 .4時(shí)移地震流動(dòng)單元解釋

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

1）流動(dòng)單元定義

class CQISSegFlowUnit

{Public：

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

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

void Display（）；//流動(dòng)單元顯示

int InsertVex（StSeedType seed）；//插入解釋點(diǎn)

void DelSeed（bool bMove=true）；//刪除解釋點(diǎn)

void MoveLine（StSeedType seed）；//移動(dòng)解釋面

private：

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

bool checkSam（StSeedType seed）；//檢查解釋點(diǎn)

int InsertVexIn（StSeedType seed）；//在測線剖面上插入解釋點(diǎn)

void InsertMoveLineIn（）；//在測線剖面上插入移動(dòng)解釋面

private://定義成員

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

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

StSeedType m_Seed；//當(dāng)前選準(zhǔn)解釋點(diǎn)

StFlowUnitPara m_FlowUnitPara；//流動(dòng)單元屬性

}

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

struct StSeedType

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

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

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

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

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

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

}

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

struct StFlowUnitPara

{intnSelectLineType；

//流動(dòng)單元類型1：測線剖面流動(dòng)單元；

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

//3：TimeSlice切片流動(dòng)單元。

int nLineType；

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

int nFlowUnitLineType；

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

}

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

4 應(yīng)用實(shí)例

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

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

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

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

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

圖4 時(shí)移地震流動(dòng)單元解釋圖件Fig.4 Flow-unit interpretation maps of time-lapse seismic

5 結(jié)束語

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

[1]王大偉，劉震，趙偉，等.利用時(shí)移地震資料劃分油藏流體流動(dòng)單元的可行性分析[J].地球物理學(xué)報(bào)，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ù)定量研究流動(dòng)單元的方法——以鄂爾多斯盆地W93井區(qū)為例[J].吉林大學(xué)學(xué)報(bào)：地球科學(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油藏流動(dòng)單元的定量劃分[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í)移地震技術(shù)應(yīng)用的可行性研究[J].勘探地球物理進(jìn)展，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.