五維地震油氣識別方法?

印興耀 張洪學 宗兆云

(1 中國石油大學(華東)地球科學與技術學院 青島 266580)

(2 海洋國家實驗室海洋礦產資源評價與探測技術功能實驗室 青島 266071)

0 引言

基于五維地震數據的地震流體識別源于寬方位地震勘探處理技術的成熟與發展，寬方位地震勘探可以獲得海量的高品質寬方位地震數據，基于炮檢距向量片(Offset vector tile,OVT)處理技術[1?2]得到TB 級高品質“五維”疊前地震道集，進而可進行五維地震數據解釋[3]。寬方位地震勘探是指觀測系統中的橫向與縱向排列的比值大于0.5 的地震勘探方式，相較于傳統的窄方位地震勘探，寬方位地震勘探有很多優勢：寬方位勘探可以增加勘探照明度，獲得較完整的地震波場；寬方位地震可研究振幅隨炮檢距和方位角的變化(Amplitude variation with offset and azimuth,AVOAZ)、地層速度隨方位角的變化(Velocity variation with angle,VVA)，從而增強了識別地下各向異性和流體分布的能力[4?6]。寬方位地震采集在21 世紀初率先在海上得到廣泛應用，之后，在陸上逐漸得到應用[7]。OVT 數據域五維地震資料，是以“片”為單位建立同時包含炮檢距和方位角信息的高品質“五維” (即空間三維坐標＋炮檢距＋方位角)共反射點地震道集。充分考慮五位地震數據中重要的方位角和炮檢距信息，挖掘其中豐富的方位各向異性信息，為五維地震流體識別奠定基礎。

隨著地質勘探目標復雜性及我國能源需求的不斷提高，以很強非均質性和各向異性為特征的復雜裂縫儲層成為我國油氣勘探的重要目標，也是目前面臨的一類最為重要的、最難預測的油氣藏，其預測成為公認的難題之一[8?10]，究其原因是目標儲層的特殊性、地下埋藏條件的復雜性以及相應地球物理數學特征的多解性[11]。各向異性巖石物理是五維地震流體識別的關鍵，疊前各向異性反演算法是實現五維地震流體解釋的有效途徑。五維地震流體識別是在各向異性巖石物理理論指導下建立將與孔隙流體有關的異常特性表征為各向異性因子，并利用五維地震數據實現流體因子反演的過程。

各向異性因子的構建是五維地震流體識別的關鍵，在五維地震巖石物理模型及儲層流體地震響應分析的基礎上，將孔隙及裂縫對流體的異常特性表征為各向異性因子，進而實現復雜儲層含流體類型的判識。Hudson[12]基于線性和硬幣型裂縫假設，提出裂縫流體指示因子。Smith 等[13]率先利用地震波速度相對差異的加權疊加作為儲層油氣識別的流體因子。Schoenberg 等[14]通過線性滑動理論將裂縫法向柔度與切向柔度之比作為另一種裂縫流體因子，目前應用較為廣泛。Goodway等[15]提出的λ-μ-ρ法即利用地下地層的異常拉伸特性識別儲層流體類型。Smith 等[16]提出了流體因子角和交叉角度的兩種流體因子概念。Quakenbush 等[17]提出了泊松比作為流體因子進行油氣識別。Russell 等[18]將參數ρf作為流體因子，隨后在2011年利用孔隙彈性巖石物理理論提出了等效流體模量Kf[19]。Xue 等[20]假設裂縫表面光滑，利用來表征流體因子，其中?N和分別表示代表飽巖石和和干燥巖石的裂縫法向柔度。上述流體是基于不同的假設理論提出，其中只有幾種與各向異性有關的流體因子，然而缺乏孔隙流體理論的支撐，由多孔隙流體理論構造的流體因子沒有考慮各向異性信息的影響。

疊前地震反演是實現五維地震流體識別的有效途徑[21]，其結果直接影響流體因子的可靠性[22]。在疊前地震反演中，各向異性反射系數近似方程是反演框架構建的基礎，如：Rüger等[23?26]給出HTI介質反射透射方程，并進行預測；Downton 等[27?28]利用Rüger 的HTI 介質反射系數近似公式，進而預測各向異性參數；印興耀等[29]利用彈性阻抗反演實現了Gassmann 流體/孔隙項的估算；陳懷震等[30?32]研究HTI 介質各向異性AVO 特征，并進行各向異性參數反演，求取裂縫流體因子；印興耀等[33]基于包含Russell流體因子的兩項彈性阻抗方程實現了Russell 流體因子的直接提取，取得了較好的深層流體識別應用效果；孫瑞瑩等[34]結合測井信息利用隨機反演算法實現流體因子估算；Pan等[35]基于AVOAZ 反演巖石物理反演裂縫流體因子；Zong 等[36]基于HTI 介質提出了預測各向異性參數的新方法；Pan 等[37]對OA 介質儲層進行裂縫參數反演，并預測儲層流體分布。由上述分析可知，目前地震各向異性反演常用的方法是彈性阻抗反演方法和基于貝葉斯框架的AVOAZ 反演方法，其中彈性阻抗反演方法由于其穩定性高應用較為廣泛。

該文結合Gassmann孔隙流體各向異性介質理論、線性滑動理論以及Thomsen 各向異性參數，提出一種新的流體因子類型，即各向異性因子(fani)，經驗證該流體因子在各向同性的情況下，可退化為Gassmann流體項，驗證該流體因子的準確性。在此基礎上，建立各向異性因子與地下介質彈性參數的關系，并利用五維地震數據反演進行流體識別，本方法從巖石物理機制上綜合考慮儲層的各向異性特性以及含流體特性儲層流體的影響，通過實際應用以及流體判識敏感性定量分析發現，各向異性因子(fani)有較高的流體指示敏感性，利用該方法進行五維地震流體識別可取得良好效果。

1 各向異性因子構建

地下巖石是由基質、裂縫和孔隙流體組成的各向異性多孔雙相介質。孔隙流體以及裂縫的存在必然會影響巖石的力學特征，從而使介質巖石表現出一定的彈性異常，以此從各向異性雙相介質理論出發可以更好地研究孔隙流體以及裂縫對介質巖石彈性參數的影響，為敏感流體因子的構建提供更好的指導。關于雙相各向異性介質理論的研究，1951年，Gassmann[38]在一系列假設前提下，推導出表征雙相各向異性介質彈性模量參數的Gassmann方程，該方程表示為

在裂縫參數較小、流體體積模量遠遠小于固體顆粒的假設下，方程(3)可以簡化為

基于長波長假設條件，單組裂縫充填于各項同性介質背景，可被視為有效的HTI 介質。在裂縫法線平行于x軸的情況下，使用線性滑動理論，干燥HTI介質剛度矩陣由方程(5)給出：

其中，λ、μ是圍巖的拉梅常數，裂縫參數為切向柔度?T與法向柔度?N。

將方程(5)代入方程(2)，可以求出

根據方程(6)，分析裂縫以及孔隙流體的影響，按照Russell 對流體因子構建流程[18]，構建各向異性因子為

Bakulin 等[40]基于線性滑動理論研究了多種裂縫組合形態的等效彈性矩陣，并在弱各向異性假設下給出了裂縫參數與Thomsen 各向異性參數的關系：

將方程(6)代入縱橫波速度公式，并結合方程(8)，將方程(7)各向異性因子改寫為如下形式：

為了驗證所建立各向異性因子的準確性，當地下介質不含裂縫時，裂縫柔度?N和?T等于0，即將方程(5)退化為各向同性介質剛度矩陣，代入方程(7)，本文構建的各向異性因子退化為

經驗證，方程(10)與Gassmann流體項相等，驗證了本文構建的各向異性因子的準確性。

為了進一步說明本文構建的各向異性因子在工區流體識別中的優越性，將幾種常用的流體因子(主要是μρ、IP、λρ、Gassmnanf、λ ?μ)對流體區分的指示系數[11]與各向異性因子進行比較，圖1顯示了不同流體因子的指示系數，從其中可以看出各向異性因子fani對流體的區分能力最強。

圖1 流體因子指示系數比較圖Fig.1 Comparison of sensitivity of different fluid factors

2 方位各向異性彈性阻抗反演算法

疊前方位地震反演方法中的方位各向異性彈性阻抗方法具有較好的抗噪性和穩定性，且能夠充分利用方位地震數據的方位各向異性信息，因而在裂縫型儲層預測方面具有廣泛的應用。Rüger[23]推導了HTI 介質的方位反射系數近似方程，使針對裂縫型儲層的疊前方位地震反演成為可能，該反射系數近似公式的另一種以縱橫波模量、密度以及各向異性參數表征的形式為

其中，θ為入射角，φ為方位角，?δ(V)、?ε(V)和?γ(V)為上下介質各向異性參數的差值。

基于方程(11)，根據Connolly[41]推導彈性阻抗的思路，將方位反射系數表示為

不同彈性參數的相對反射系數?M/M、?μ/μ以及?ρ/ρ可以用對數形式來等效：

將方程(12)和方程(13)代入方程(11)可以得到HTI介質方位各向異性彈性阻抗方程。由于方位各向異性彈性阻抗方程的量綱會隨著入射角和方位角的變化出現劇烈變化的不穩定現象，因此為了得到可靠預測地應力的彈性參數和各向異性參數，需要對方位各向異性彈性阻抗做標準化處理[42]，最終整理得到標準化的方位彈性阻抗方程為

其中，

為了能夠從彈性阻抗數據體中反演出穩定的彈性參數和各向異性參數數據體，需要對方程(14)進行線性化處理，即：

利用已知的測井數據以及井旁道方位彈性阻抗值計算公式(16)中各參數對應的常系數。在儲層彈性性質橫向變化較小的情況下，可將常系數看作是恒定不變的值。將疊前方位彈性阻抗反演得到的6個不同方位、不同入射角的彈性阻抗代入公式(16)中，結合各參數對應的常系數，通過求解可得到各采樣點處的彈性參ρ和各向異性參數。由于該方法在進行反演時使用了測井資料進行約束，并結合地震資料的先驗信息，因此得到的彈性參數和各向異性參數較為準確可靠。

3 實際應用

五維地震資料來自中國東部某寬方位地震勘探工區，該區域各向異性特征明顯。對五維地震數據首先進行預處理得到6個方位且每個方位包含3個入射角的共18個分方位部分角度疊加道集，然后選取其中6個疊前方位地震數據進行疊前方位彈性阻抗反演，疊前方位地震數據如圖2所示。從圖2中可以看出由于儲層中存在的裂縫對地震資料的反射振幅產生影響，因此在相同入射角的情況下，不同方位的地震數據之間存在明顯的差異。

根據圖2 五維地震數據，結合該區域井資料以及巖石物理信息作為約束，反演得到不同方位的大中小角度彈性阻抗如圖3 所示。通過基于方位彈性阻抗的各向異性參數反演方法從方位彈性阻抗數據體中提取地層的縱波模量、橫波模量、密度和HTI介質的各向異性參數，進而求取剖面含流體區域的各向異性因子，如圖4 所示，A 井和B 井的有色部分代表流體區，白色部分是干巖區。圖4 中的橘色箭頭所指部分為測井解釋的流體部分，發現該工區儲層中，流體分布區域的各向異性因子在反演數值呈現低值特征，而且各向異性因子數值在巖石發育位置展現出高值特征，該特征與巖石物理分析規律相吻合。同時，含流體裂縫性儲層各向異性流體的估算剖面與A 井和井B 的流體解釋結果具有良好的匹配性。因此，通過各向異性反演各向異性因子可展示儲層流體分布特征，為五維地震定量解釋創建新的理論方法，該方法為復雜裂縫儲層地震流體識別提供了一種新的途徑。

圖2 實際工區不同方位部分角度疊加道集Fig.2 Stack data of different azimuths for reservoir

圖3 實際工區方位各向異性彈性阻抗反演結果Fig.3 Azimuthal anisotropic EI inverted results for reservoir

4 結論

本文結合Gassmann雙相各向異性理論和線性滑動理論推導出的含流體HTI 介質的剛度矩陣，構建了一種新的各向異性因子fani，在弱各向異性假設條件下，建立各向異性因子fani與彈性模量和Thomsen參數的關系。多種流體因子的敏感性比較表明，各向異性因子fani在各向異性儲層比其他流體因子更敏感。最后，充分利用五維實際地震數據中豐富的方位各向異性信息，通過方位彈性阻抗反演實現各向異性因子的求取，分析可知，反演結果與測井解釋資料有很好的吻合，該方法為五維地震流體識別提供一種新的途徑。綜上，為提高儲層流體類型的判識準確度，構建能夠精細反映各向異性儲層油氣分布的流體因子是基礎，選用有效的地震反演算法是關鍵，精度高、穩定性好的各向異性反演算法可以有效避免誤差影響，有助于儲層流體類型的可靠判識，在實際資料應用中，應盡可能選用全方位五維疊前地震數據提高各向異性反演的可靠性。