寬頻寬方位數據同時反演及應用實例

張大海，張 紀，劉 苗，蔣 涔，王慧欣，楊雨松，巢 越

（中國石油化工股份有限公司上海海洋油氣分公司勘探開發研究院，上海 200120）

隨著勘探技術的不斷提高，Y 工區開展了寬頻寬方位的三維地震資料采集處理，雙寬地震資料的信噪比、分辨率以及成像精度均大幅提高。利用資料的寬頻及寬方位特點有助于深部儲層的描述，結合疊前同時反演結果，對深部儲層邊界及縱、橫向不均質性進行精細的識別與刻畫，尋找深部致密砂層中的“甜點”，為后期挖掘西湖凹陷深部儲層油氣潛力奠定基礎。

1 工區概況

東海Y 工區位于西湖凹陷中北部，根據已鉆井資料分析，工區內主要含油氣構造花港組砂巖分布廣泛，砂體厚度穩定，受沉積微相與巖石粒度影響，儲層連續性差，物性變化快，非均質性強。

Y 工區5 口井都有油氣顯示，其中Y1 井最好，H3～H9 氣測顯示共328.25 m/22 層，5 口井分別在H3、H5～H7 及H9 砂層組進行了4 層DST 測試，在Y1 井和Y5 井H3 層測試獲得成果，日產氣分別為20×104m3/d 和1.1×104m3/d。從這5 口井的鉆探來看，Y 地區砂體發育，但物性差異較大，而氣層影響最大的因素是儲層物性。

從鉆井揭示的儲層物性來看，H3 層儲層物性相對較好，但橫向上存在較強的非均質性，各井厚度變化范圍在140～195 m 之間，巖性包含細砂巖，含礫砂巖以及中砂巖，孔隙度變化范圍在2%～11%，Y1、Y4、Y5 井含有油氣層；縱向上H4 及以下砂層組偏致密，為特低孔超低滲儲層，局部發育甜點。因此，有利儲層的預測是Y 地區勘探的關鍵因素，主力氣層H3 整體表現為中低和特低滲儲層，“甜點”相對發育，孔隙度與滲透率相關性較高，因此，尋找深部致密砂層中的“甜點”、減少優質儲層厚度預測誤差是推動Y 工區深部含氣構造評價的關鍵。

2 寬頻寬方位地震數據

寬頻數據對低頻信息的補充使儲層研究中的彈性阻抗及AVO 響應更準確，顯著提高了對巖性-流體的分布及孔隙度的預測可靠性[1-3]。寬方位采集主要有三個優點：寬角度的照明有利于復雜結構下的成像；高覆蓋，高信噪比，多次波干擾得到壓制；豐富的方位角信息使獲得可靠的速度模型成為可能。寬頻寬方位角采集的特點對數據處理也提出了更高的要求。

圖1 分別展示了處理后數據在淺層及深層對頻帶的拓寬結果。斜纜寬頻的采集方式以及有效的三維去鬼波處理方式提高了數據中低頻信號的信噪比，對斷層成像和反演穩定性的提高都有極大幫助，同時有效拓展了數據頻帶，提升了數據分辨率。

圖1 雙寬地震（藍色）與常規地震（紅色）頻譜對比圖Fig.1 Spectrum comparison of double-wide seismic (blue) and conventional seismic (red)

3 寬頻寬方位地震反演

Y 工區雙寬數據疊前同步反演流程和常規地震數據疊前反演流程基本一致，為了更好地挖掘雙寬數據潛力，同時兼顧提高反演精度，在子波的計算及選取和低頻模型的約束試驗等關鍵流程的處理環節采取了不同的處理方法，這些方法的應用直接影響到反演結果的可靠性和有效性[4-6]。Y 工區雙寬數據方位角分布較窄，范圍0°～89°，入射角范圍0°～45°，選取有效范圍為3°～39°，分別分為3°～15°、15°～25°、25°～39°三個部分疊加。

3.1 寬頻子波估算

作為聯系地震與測井紐帶的合成地震記錄，其中子波求取的準確性直接影響到地質層位的標定結果以及反演結果的精度，因而合適子波的求取是決定反演成敗的重要影響因素。因此只有通過子波反演和層位標定多次迭代才能獲取合適的子波以及標定。

寬頻地震帶寬2.5～50 Hz，應用寬頻地震估算振幅譜子波，為了加強疊前地震道集信噪比及最大可能地保留疊前道集信息，通過井上正演分析，分別提取了近角度、中角度和遠角度子波，子波長度則為300 ms，如圖2 為300 ms 寬頻地震子波與104 ms 長度常規地震子波對比。

圖2 分角度疊加體提取子波對比圖Fig.2 Comparison diagram of the extraction of wavelets by an angle-stacking seismic

層位標定選擇了工區內測井優化處理后的曲線來進行，首先根據實際地震資料目標層附近的有效頻帶寬度，選取寬頻數據地震子波，單井的分角度疊加數據標定結果（圖3）顯示，測井曲線上存在的巖性界面與地震剖面上同相軸的對應關系基本一致。

圖3 Y1 井合成地震記錄標定圖Fig.3 Calibration diagram of synthetic seismic records of Well Y1

從圖3 的地震合成記錄標定結果可以看出，主要目的層段正演合成地震與原始地震相關系數較高，地震波組對應關系一致性較好，振幅能量也基本匹配，不同角道集數據標定中原始地震記錄與合成地震記錄的相關系數都比較高，其中近道均超過了0.7；遠角度稍低，在0.4 左右，原因在于遠角度地震數據相對近、中角度本身品質較低。相比以往常規數據的標定相關性，雙寬數據標定中相關系數更高，標定過程的不確定性更小。

近期，《中國矯形外科雜志》編輯部多次接到讀作者的電話和Email，發現有多個網站利用《中國矯形外科雜志》名義非法征稿及騙取有關費用，要求作者將費用匯入指定賬戶等方式騙取作者錢財，侵犯了廣大作者的合法權益?！吨袊C形外科雜志》編輯部在此提醒廣大讀作者，本刊編輯部從未委托任何代理機構為《中國矯形外科雜志》征稿。

3.2 低頻模型搭建

針對Y 研究區儲層橫向變化快的特點，利用測井曲線構建低頻模型方法通常存在不適應性，地層橫向上井控影響太大。因此同時可以考慮引入地震數據中能夠提供的振幅信息和速度信息來改善低頻模型。通過地震處理獲得的地震速度，能較好地反映沉積的變化趨勢。因此，采取多屬性融合的思路建立低頻模型，結合壓實趨勢、構造格架、地震速度、測井曲線、巖石物理關系等多種屬性信息融合建立高精度低頻趨勢模型，為高質量疊前反演打下基礎。結合層位、速度體和測井曲線同時構建初始低頻模型。

地震速度不僅能夠提供低頻信息，而且是疊前地震反演的重要輸入，用于將偏移距域道集轉化到入射角域。原始的地震速度一般都存在局部速度異常、與井速度存在一定差異的問題。應用之前進行了精細的地震速度優化工作，主要包括層速度轉化及空間插值、低通濾波、用井速度校正層速度等，盡量減少速度異常值對速度譜的影響，為低頻模型建立提供合理的約束條件。

低頻模型能夠反映實際的地質趨勢，儲層反演的作用之一是消除子波調諧影響，能真實反映一定厚度的地層特征，但實際操作中反演結果仍然受調諧的影響。為了盡可能降低調諧的影響，優化構建低頻模型，從而保證精確求取反演彈性參數值。 基于上述特點，本次低頻模型建立過程中，借助JASON 軟件包中EMFT 模塊功能，充分融合地震速度、測井以及地質模型信息，通過迭代的方法求取低頻模型，使低頻模型盡可能符合實際的地下沉積特征，為儲層預測反演提供良好的支持，圖4 為研究區速度約束的測井曲線插值縱波阻抗聯井剖面圖，可見低頻趨勢合理，無明顯牛眼現象，利于后續反演研究。

圖4 速度體約束縱波阻抗低頻模型聯井對比剖面圖Fig.4 Well correlation profile velocity volume constrained P-wave impedance low-frequency model

3.3 疊前同時反演

寬頻地震子波提取、精細井震標定、建立以精細構造模型為基礎并以速度體做約束的低頻模型是本次雙寬數據疊前同步反演區別于常規數據反演的特點。圖5 具體顯示了本次疊前同時反演中關鍵流程。

圖5 寬頻地震疊前同時反演的流程圖Fig.5 Flow chart of broadband seismic simultaneous inversion

通過疊前同時反演，可以得到縱橫波阻抗、Vp/Vs、密度及E/LAMDA 等彈性參數?；趲r石物理分析結果，反演得到的Vp/Vs對巖性具有較好區分能力，泥巖表現為中高Vp/Vs值，砂巖表現為中低Vp/Vs值，相比與應用常規地震反演，寬頻地震反演對薄、厚層地質體均有較高分辨率，砂巖、泥巖疊置層次清晰（圖6）。

圖6 疊前反演Vp/Vs 對比聯井剖面Fig.6 Prestack inversion Vp/Vs well correlation profile

4 反演結果討論

本次巖石物理分析遵循了“嵌套”識別的原則，首先進行砂、泥巖巖性的區分，在此基礎上進一步區分物性較好的深部“甜點”儲層。

根據巖性劃分結果，對H3 儲層進行了精細的巖石物理交匯分析。

圖7 顯示，Y 工區H3 主力儲層巖性可用縱橫波速度比參數進行區分，其中泥巖Vp/Vs＞1.68（藍色部分），砂巖Vp/Vs＜1.68（紅色部分）；在砂巖已識別的基礎上，甜點儲層的識別可通過縱波阻抗與E/LAMDA 彈性參數組合進行，縱波阻抗及E/LAMDA雙參數的嵌套識別法與實測結果吻合度較高。

圖7 井上巖石物理參數交會圖Fig.7 Crossplot of petrophysical parameters on the well

根據Y 構造H3 儲層巖石物理分析結果，以Vp/Vs＜1.68 為標準，進行砂體的識別。

圖8 為根據反演結果提取的H3 儲層段Vp/Vs最小值屬性圖，圖9 則是Vp/Vs＜1.68 砂體樣點累計時間厚度圖，圖示H3 砂體全區發育，其中Y1 井以南的研究區中部存在低Vp/Vs區域，同時砂體樣點累計時間厚度也較厚，推測該部位砂體非常發育，同時該部位斷層西側也有較好砂體分布。結合已鉆井分析，Y1、Y2 井位于該砂體，而Y4、Y5 井則位于該砂體以外；工區西北部Y3 井Vp/Vs值較高，同時厚度減小，推測其性質不同于工區中部砂體。

圖8 Y 構造H3 層砂體分布平面圖Fig.8 Distribution map of sand body in H3 layer of Y structure

圖9 Y 構造H3 層砂體時間厚度圖Fig.9 Time thickness map of sand body in H3 layer of Y structure

圖10 為H3 砂體聯井剖面對比圖，通過該段砂體與實測砂體厚度的對比，本次砂體的厚度預測基本符合實測結果（圖中測井H3 色段與背景反演結果對比），其中Y1、Y2、Y3 井與反演結果吻合度較高，而Y5 井稍有差異，分析認為該井H3 段下部砂巖Vp/Vs反演值稍高于門檻值。

圖10 Y 構造H3 砂體聯井對比剖面圖Fig.10 Comparison profile of H3 sand-body joint well constructed in Y structure

在砂體刻畫的基礎上，進一步利用E/LAMDA+AI（縱波阻抗）彈性參數組合刻畫H3 的“甜點”儲層。根據前述甜點儲層的地質特征，在深部致密砂體中尋找物性較好的“甜點”儲層首先是要界定其沉積微相為河道主體部位（心灘、河道砂壩），上述H3 砂體顯示其河道主體部位特征，繼續應用巖石物理分析結論，結合E/LAMDA 高值區+縱波阻抗低值區預測本研究區H3 深部“甜點”儲層的可能分布。

本次定量解釋的依據主要是以Y1 井H3 甜點儲層參數范圍，其中E/LAMDA＞3.3，縱波阻抗范圍為10 000～11 800 kg/cm3·m/s，圖11 為根據上述敏感彈性參數組合解釋得到的H3“甜點”儲層頂面構造圖，圖12 則是H3“甜點”儲層時間厚度圖。由圖可見，甜點儲層主要發育于工區中部大斷裂兩側構造高部位，其中井間的差異性有較好的體現（Y1 井、Y2 井鉆遇低滲氣層，其余三口井并無甜點儲層），另發現Y4井E/LAMDA 值偏高，如用Y1 井的參數進行解釋會增加預測風險，因此“甜點”儲層的預測受井控影響較大，需反復對比以確定適用范圍。

圖11 Y 構造H3 甜點儲層頂面構造圖Fig.11 Top layer map of Y-structure H3 sweetness reservoir

圖12 Y 構造H3 甜點儲層時間厚度圖Fig.12 Time thickness map of Y-structure H3 dessert reservoir

圖13 Y 構造預測甜點儲層空間分布圖Fig.13 Spatial distribution of Y structure prediction sweetness reservoir

根據定量解釋“甜點”儲層與實測的對比結果，H3“甜點”儲層的預測給出了有利目標區，3D 圖顯示了H3“甜點”儲層的空間展布（圖中紅色范圍所示），顯示其潛力較大。

5 結論

（1）通過雙寬數據與常規數據的子波、低頻模型、反演剖面、平面的對比，充分顯示了雙寬數據具備的優勢，以及雙寬數據反演的穩定性及可靠性。

（2）通過對Y 工區雙寬數據進行疊前同時反演，獲得了與井吻合度較高的反演結果，表明寬頻地震數據有效提高了疊前反演質量，標定準確，儲層邊界清晰，儲層內部差異明顯，在此基礎上刻畫了Y 構造H3 主力儲層砂體和“甜點”儲層的展布，結果顯示甜點儲層識別與實測井吻合。

（3）Y 工區雙寬數據在疊前反演中的應用效果比較顯著，疊前反演得到的彈性參數體與井的吻合度較高，儲層的邊界及內部的差異都更清楚，隨著更多雙寬數據的應用，其優勢也將得到更多的驗證，將為西湖凹陷深部儲層的進一步開發提供更多的數據支撐。