渤海油田深埋潛山儲層預測關鍵技術研究*
——以渤中19-6潛山為例

周東紅

(中海石油(中國)有限公司天津分公司渤海石油研究院 天津 300452)

渤中凹陷作為渤海油田油氣發現主戰場，經過40多年的勘探，在新近系以及古近系發現了數量眾多的油氣田，儲量規模占總發現的94%，而潛山儲量發現僅占6%，類比鄰區地質條件相似的冀中凹陷、遼河凹陷，潛山儲量發現占比分別達到59.7%和20.3%。隨著中淺層勘探程度的提高，中淺層尋找規模較大的構造圈閉難度越來越高，位于凹陷區的深埋潛山正逐漸成為渤海油田的重點勘探領域之一。

近些年來，中國海油在渤中凹陷開展了一系列的勘探研究工作，明確了位于渤中凹陷西南部的渤中19-6構造區具有發育規模性潛山構造圈閉的條件[1]。但由于該區構造平均埋深大于4 500 m，且處于郯廬斷裂和張蓬斷裂的應力疊加區，斷裂系統極復雜，裂縫型儲層非均質極強；加之該區潛山之上新生界地層中還發育多期火山巖，受多期火山巖屏蔽影響，地震資料品質差，導致渤中19-6深埋潛山在勘探階段的儲層預測面臨巨大挑戰。

針對渤中19-6深埋潛山儲層預測難度大的問題，本文在借鑒國內外類似變質巖潛山地震勘探關鍵技術的基礎上，結合該區資料基礎及地質特點，從凹陷區火山巖下伏地層地震波能量屏蔽補償、潛山風化裂縫帶、內幕裂縫帶儲層地震預測等方面開展研究，綜合利用基于相位一致性分析的微斷裂刻畫技術、基于古地貌約束的疊前彈性參數融合技術、基于曲波變換的內幕大尺度斷裂增強處理技術、基于繞射波分離的小尺度裂縫帶預測技術，提高了變質巖潛山裂縫儲層的預測精度，形成了渤海海域變質巖潛山地震勘探關鍵技術體系，并在渤中19-6大型凝析氣田實際勘探評價過程中取得了良好的應用效果。

1 太古界變質巖潛山風化裂縫帶儲層預測技術

在儲層預測過程中，地震資料品質直接決定了儲層預測成果的預測能力與可靠性。由于渤中19-6構造中深層火山巖發育，構造區約70%面積受高速火山巖屏蔽影響，造成下伏潛山地層反射振幅明顯變弱，將嚴重影響后續儲層預測的精度。目前業界大都采用基于均方根振幅增益控制的方式或者基于已鉆井統計的分頻補償方法進行解決[2-3]。前一種方法破壞了振幅之間的相對能量關系，是一種不保幅的處理手段；后一種方法采用基于已鉆井合成地震記錄統計特征進行補償，不可避免地忽略了由于不同偏移距路徑差異對反射能量的影響，加之受限于研究區已鉆井數量少等因素，其適用性及保幅性均不理想。因此本文采用了基于射線路徑的復雜地質體下伏地層振幅補償技術對火山巖下伏地層反射能量進行補償處理。經過補償處理后，潛山地層真實反射能量得到恢復，潛山內幕斷裂成像更加清楚，提高了特殊地質體覆蓋區地震資料保真度，為后續深埋潛山儲層精細表征奠定了資料基礎。

1.1 基于相位一致性分析的微斷裂刻畫技術

風化帶儲層發育情況和微斷裂及裂縫展布相關，在微斷裂及裂縫的地震刻畫方面，以曲率、螞蟻體等為代表的地震幾何屬性得到廣泛應用。然而，渤中19-6潛山上覆地層為超壓低速泥巖，上、下地層分界面波阻抗差異大，導致地震剖面上產生強振幅反射，這種由強波阻抗差引起的強反射會“淹沒”裂縫對地震子波波形的改造作用，這種情況下利用基于振幅信息的地震不連續性檢測會存在多解性[4]。

為此，引入基于相位一致性分析的微斷裂識別方法。相位一致性是信號中各個頻率成分的相位相似度的一種度量方式，是一個無量綱的量，其本質雖然仍是基于地震信號的不連續性檢測，但與常規方法不同，該方法主要利用相位信息，認為地震信號中不同頻率的傅里葉分量相位最一致的點為特征點，其在時間域具有不連續性特征。因此可以通過相位一致性分析進行地震不連續性檢測，不受地震振幅的橫向變化影響，也無需低通濾波去噪[5]。

對于一維信號f(x)，其傅里葉展開式為

f(x)=∑nAncos(nωx+φn0)=∑nAncos(φn(x))

(1)

式(1)中：An為第n個傅里葉展開分量的振幅值；ω為角頻率；φn0為第n個分量的相位偏移量或初始相位；φn(x)表示位置x處的第n個傅里葉分量的局部相位。

則相位一致性度量函數定義為

(2)

圖1為渤中19-6構造北1塊沿層提取的相位一致性屬性與常規傾角相干屬性的對比。新方法刻畫的裂縫細節更加豐富，集中發育帶展布特征清晰，而常規相干受強振幅影響，未能有效刻畫風化帶裂縫發育規律。同時，圖1a中BZ19-6-4井位于高PC值處，預測裂縫發育，實鉆揭示該井風化帶發育儲層119 m，凈毛比達92%，且測試獲得高產，有效證實了新方法的可靠性。

圖1 相位一致性(a)與傾角相干(b)切片對比 Fig.1 Slice comparison of phase consistency (a)and dip-scanning coherent (b)

1.2 基于古地貌約束的疊前彈性參數融合技術

巖石物理分析表明，風化裂縫帶儲層具有低橫波阻抗、低密度特征，利用疊前彈性參數反演可以預測風化裂縫帶儲層展布。首先以裂縫總孔隙度曲線作為目標曲線，橫波阻抗曲線與密度曲線作為計算曲線，通過參與計算曲線的線性匹配關系模擬目標曲線，并與先驗目標曲線進行誤差分析。基于最小誤差原則，構建基于橫波阻抗和密度雙參數融合的風化帶儲層敏感因子為

F=2.53×10-5SI-0.172 9Rho+0.71

(3)

式(3)中：F為儲層敏感因子；SI為橫波阻抗，kg·m-3·m·s-1；Rho為密度，g·cm-3。

然后，通過稀疏脈沖反演方法得到橫波阻抗及密度參數體，進而得到最終儲層敏感因子體。圖2為疊前反演得到的儲層敏感因子剖面，靠近潛山頂面低敏感因子區為有利儲層發育區，其與潛山內幕裂縫帶能明顯區分開，并且與已鉆井吻合較好。

圖2 疊前儲層敏感因子反演剖面Fig.2 Pre-stack inversion profile of reservoir sensitivity factor

渤中19-6潛山埋深大，地震疊前道集有效入射角范圍小，加之海上勘探階段鉆井少，潛山疊前反演多解性強，單一疊前彈性參數屬性具有一定局限性。從儲層地質成因分析看，變質巖風化帶儲層發育程度及分布特征受地貌形態影響較大，古地貌的高部位風化淋濾作用更強，更有利于儲層發育，這一認識與現有的已鉆井結果相吻合[6-7]。因此，在綜合地質成因、地震響應特征的基礎上，構建了基于古地貌的多屬性融合技術。首先通過印模法分析潛山上覆古近系地層厚度的鏡像關系，識別古地貌單元；然后沿潛山頂面風化帶提取疊前儲層敏感因子反演屬性，并將其與古地貌信息進行線性融合，形成了以疊前橫波阻抗、密度和潛山古地貌參數3種數據聯合的預測思路和計算方法，即

(4)

式(4)中：F為融合后的優質儲層指示屬性；T為古地貌；k為古地貌融合系數，用于調節古地貌約束程度，可結合鉆井在0～5范圍取值。

圖3a為渤中19-6構造新生界沉積之前的古地貌，可見在印支—古近紀，渤中19-6構造整體具有北高南低的特征，其中北塊BZ19-6-47井區和南塊BZ19-6-12Sa井區均為古地貌高部位，其風化淋濾作用更強，更有利于風化帶儲層發育。將上述地質規律融入疊前儲層敏感因子反演結果(圖3b)，得到古地貌約束后的預測結果(圖3c)。從圖中可見，對于鉆井揭示的風化裂縫帶儲層整體發育的BZ19-6-47井區，融合后的新屬性明顯表現出一個整體發育的特征，與鉆井吻合更好。在南塊鉆井揭示古地貌高部位的BZ19-6-1井比BZ19-6-3井儲層更發育，這種差異在融合后的新屬性上表現明顯。總之，融合后的屬性比單一的疊前儲層敏感因子反演結果更符合地質規律，預測精度得到明顯提升。

圖3 渤中19-6構造古地貌(a)及古地貌約束前(b)后(c)疊前儲層預測結果Fig.3 Palaeo-geomorphology of BZ19-6 structure (a) and the pre-stack reservoir prediction results before (b) and after (c) constrained of palaeo-geomorphology

2 太古界變質巖潛山內幕裂縫帶儲層預測技術

渤中19-6太古界潛山自印支期以來經歷了多期構造運動，產生的不同走向的裂縫組成裂縫網絡，為后期油氣聚集提供了良好的儲集場所。內幕裂縫由于埋深較深，鉆探難度大，探井數量少[8]。潛山內幕裂縫的存在與內幕裂隙發育有較大聯系，微觀尺度裂隙發育方向存在一致性，而發育密度較高時，在宏觀上表現為大的裂縫發育特征。

圖4為建立的不同尺度的潛山內幕裂縫的正演模型，其中圖4a為具有不同裂縫寬度的大尺度潛山內幕斷裂模型，圖4c為中、小尺度密集型微斷裂及裂縫模型，圖4b、d分別為對應的正演偏移剖面。正演結果表明:變質巖潛山內幕斷裂或大尺度裂隙帶發育是產生高角度地震反射的主控因素;低頻、高陡、中強連續反射特征，是大尺度內幕斷裂發育區帶的直接指示特征。中尺度密集型和小尺度彌散型裂隙發育會導致地震同相軸連續性變差，具有明顯斷續反射特征，表明其是產生地震散射或繞射的主控因素，散射或繞射波的強弱與中、小尺度裂縫發育程度有較好的一致性。

注：圖中數值為斷裂帶寬度，m。圖4 不同尺度潛山內幕斷裂模型及正演分析Fig.4 Fracture models and forward modeling of buried-hill at different scales

從渤中19-6構造地震剖面來看，潛山內幕地層存在大量的高陡連續反射特征，這種反射特征在某些位置表現為高陡斜交叉反射，與正演分析結論相符。由于潛山內幕埋深大，地震資料品質差，潛山內幕有效反射受其他信號干擾嚴重，潛山裂縫表征存在較大困難。如何增強潛山內幕斷層響應、挖掘繞射波信息是渤中19-6構造當前勘探階段內幕裂縫儲層預測的關鍵。

2.1 基于曲波變換的內幕斷裂增強處理技術

經歷多期構造活動的改造，潛山內幕結構和構造樣式往往復雜多樣，而內幕地震資料一般存在地震照明不均、信噪比低、斷面成像不清等問題，常規斷層增強濾波技術(如各向異性擴散濾波、構造導向濾波)更適用于層狀地層的斷層增強，對于塊狀地層低信噪比條件下的潛山內幕斷層，難以有效提升潛山內幕斷裂成像效果。

本文提出了基于曲波變換的內幕大尺度斷裂反射波增強技術。曲波變換是在小波理論基礎上發展出來的一種適合表征各向異性的多尺度方法。曲波是通過在小波基函數的基礎上添加一個表征方向的參數得到的，所以它不但和小波一樣有局部時頻分析的能力，而且還具有很強的方向選擇和辨識的能力，可以非常有效地表示信號中具有方向性的奇異特性[9-10]。

對地震信號而言，模型正演表明通過曲波變換能夠以較低的冗余度對地震資料進行尺度分解和角度分解，這種分解對本區潛山內幕的資料處理具有重要意義。通過高精度曲波變換，在曲波域進行地震數據的不同尺度分解，并針對不同尺度給出不同的重構系數，主要是通過對反映大尺度斷裂信息的粗尺度數據給予更大權重系數，使重構后的數據能夠突出大斷裂信息。再對重構后數據進行第二次曲波變換，得到不同角度數據體。為了突出高角度的斷層信息，重點對中高角度數據體進行非線性各向異性擴散濾波處理，使其在保持斷裂邊界信息的前提下去除隨機噪音，最后將不同角度的地震數據線性相加，得到內幕大尺度斷裂增強后的數據。

圖5為內幕斷層增強處理前后的剖面對比。在常規地震剖面上，由于內幕資料成像質量差、信噪比低，導致潛山內幕呈雜亂反射，無明顯斷裂響應。傾角、曲率等屬性也難以有效刻畫出其平面展布特征。經過增強處理后的數據，潛山內幕隨機噪音得到良好的壓制，內幕斷層邊界特征和復雜構造特征得到大幅增強，不同斷塊內大尺度斷裂儲層發育差異性明顯，從而有效指示了渤中19-6構造區內幕斷裂發育帶平面展布特征，進而明確了潛山內幕裂縫儲層發育的宏觀規律。

圖5 內幕斷層增強處理前(a)后(b)地震剖面對比Fig.5 Comparison of seismic profiles before (a)and after (b)inner fault enhancement processing

2.2 基于繞射波分離的小尺度裂縫預測技術

裂縫型油氣藏中，構造裂縫系統可以視為多個小尺度地質體的集合，由于其地震信號能量較弱，致使常規的反射波預測方法對裂縫系統的識別并不敏感。前人研究已經證實，繞射波場可以較好地描述小尺度地質體，當地質體的空間尺度接近或小于地震反射數據分辨率時(一般小于λ/4)，從地震數據中分離的繞射波數據能夠精細刻畫小尺度地質體。本文提出了基于主成分分析技術的繞射波提取方法，通常主成分分析技術(PCA)應用在多尺度分析、數據降維壓縮等方面，本文將其引入到反射波與繞射波波場分離中，通過2種波場振幅幅值差別和空間分布差異的特點進行波場分離[11-14]。

基于目標區三維地震數據提取繞射波數據，對比反射剖面(圖6a)和對應的繞射剖面(圖6b)可以得出，在常規地震反射剖面上可以分辨出較為連續的反射同相軸，對于常規的構造解釋工作，只能對較大級別斷層進行解釋，無法對裂縫系統的空間分布進行描述。但是在分離的繞射波剖面中可以看出，進行繞射信息分離后，目標區會保留由于裂縫系統所產生的繞射波。11井附近高角度反射特征并不明顯，但是存在明顯的繞射或散射現象。

圖6 BZ19-6-11井區反射波(a)與繞射波(b)剖面對比Fig.6 Comparison of the profile of the reflection (a) and the diffraction (b) cross the BZ19-6-11

通過在繞射波數據體上提取振幅類屬性對目標區中小尺度裂縫分布進行描述，得到裂縫密度預測結果如圖7所示(紅色為裂縫高密度分布區)。繞射波能量主要分布在2個走滑斷裂相夾持的區域和斷裂相對較發育的一些區域，此區域構造變形強烈，裂縫相對發育。鉆井統計數據顯示繞射波屬性預測成功率達87%，有力地證明了繞射裂縫預測結果能夠較好地描述工區內裂縫的分布規律，同時提取的裂縫信息與構造解釋中的深大斷裂有很好的耦合及伴生關系，為該區構造裂縫發育的地質認識及下一步潛力區劃定提供了重要依據。

圖7 基于繞射數據體的潛山內幕裂縫密度預測結果Fig.7 Fracture density prediction of buried hill interior based on diffraction data volume

3 應用成效及推廣價值

渤中19-6千億方級凝析氣田的發現突破了富油型盆地難以尋找大氣田的認識。該氣田的勘探歷程可分為發現、初步評價、整體評價3個階段，其中在發現階段地球物理研究以尋找規模性構造圈閉為主，初步評價階段地球物理研究以潛山風化裂縫帶儲層預測為主，而整體評價階段地球物理研究則以潛山內幕及風化裂縫帶儲層綜合預測為主。因此，本文提出的深埋潛山儲層預測關鍵技術成功助推了渤中19-6大氣田的發現。

同時，本文針對中深層復雜儲層提出的儲層預測技術預測系列，為渤海類似地質條件的深埋潛山的勘探評價，乃至中國近海海域其他盆地勘探進軍深層，提供了地球物理技術支撐。如：近期利用該套技術預測渤海海域428潛山、曹妃甸12-8構造、渤中8-3構造、渤中22-23等多個構造區深層優質儲層發育，勘探潛力巨大，為下一步勘探目標優選提供了參考。

4 結論

環渤中凹陷深埋變質巖潛山儲層垂向具有明顯分帶性，從潛山頂面向下，依次可分為風化裂縫帶、致密帶、內幕裂縫帶和基巖帶，其中風化裂縫帶和內幕裂縫帶為儲層發育段。綜合利用基于相位一致性分析的微斷裂刻畫技術、基于古地貌約束的疊前彈性參數融合技術、基于曲波變換的內幕大尺度斷裂增強處理技術、基于繞射波分離的小尺度裂縫帶預測技術，提高了變質巖潛山裂縫儲層的預測精度，該套技術體系有效指導了渤中19-6千億方大氣田的勘探評價，對渤海類似地質條件的深埋變質巖潛山的勘探評價具有一定的借鑒意義。