應用地球物理技術定量評價遼西凹陷沙河街組烴源巖

牛 聰 劉志斌 王彥春 張益明 解吉高 黃 饒

（①中國地質大學（北京）地球物理與信息技術學院，北京100083；②中海油研究總院，北京100028）

應用地球物理技術定量評價遼西凹陷沙河街組烴源巖

牛 聰*①②劉志斌①王彥春②張益明①解吉高①黃 饒①

（①中國地質大學（北京）地球物理與信息技術學院，北京100083；②中海油研究總院，北京100028）

針對遼東灣遼西凹陷沙河街組烴源巖，提出了一種基于地球物理技術的定量評價方法。首先在巖心地化分析數據的基礎上利用測井解釋方法進行總有機碳含量（TOC）預測；接著利用巖石物理分析得到與TOC具有強相關性的地球物理敏感參數，進而擬合兩者的關系方程；然后通過地震反演預測地球物理敏感參數數據體，并由擬合關系將其轉換為TOC數據體；最后按照烴源巖的質量評價指標，預測不同等級烴源巖的空間分布，并與盆地模型參數相結合計算目標凹陷的生烴量。實際應用表明，研究區沙河街組存在一套TOC大多為2%以上、平均厚度可達200m的優質烴源巖，估算其油氣資源量約為26億噸，該方法預測的烴源巖厚度與實測值的相對誤差小于10%。研究結果說明，利用地球物理技術可以實現對烴源巖的豐度分級和定量評價，對類似生烴凹陷的評價具有參考價值。

地球物理 遼西凹陷 沙河街組 烴源巖 定量評價

1 引言

烴源巖評價是含油氣盆地勘探研究的基礎。在過去的研究中，對烴源巖的評價常常是通過對巖樣和巖屑資料的分析形成，其缺點是只有在具有地化資料的情況下才能進行，并且受樣品數量的限制，研究成果僅局限于單點［1］。近年來，人們廣泛利用地球物理技術進行了烴源巖預測并具有廣闊的應用前景［2-8］。前人的研究表明湖相優質烴源巖具有典型的地球物理特征：在地震波響應上為空間可追蹤的低頻連續反射；在測井曲線上表現為高伽馬、高電阻率、高聲波時差、低密度的“三高一低”特征；評價烴源巖品質的總有機碳含量（TOC）、鏡質體反射率（Ro）等指標與密度、縱波阻抗、縱橫波速度比等地球物理彈性參數具有強相關性［9-13］。這些研究雖已經取得了積極的進展，但是在烴源巖豐度分級和定量評價方面的研究仍比較欠缺。

本文基于遼西凹陷的勘探實際，結合地質、沉積、地化及盆地模型等資料，探索了一種烴源巖的地球物理評價方法，確定了研究區內沙河街組烴源巖的發育規模，并對其進行了豐度評價和生烴量估算。

2 研究區概況

遼西凹陷位于渤海東北部海域，是遼東灣地區的一個次級構造單元，總體為北北東—南南西向展布，其東部為遼西凸起，面積約為3830km2，由北、中、南3個洼陷組成。本次研究區位于遼西凹陷北洼內，遼西凸起的西部（圖1）。

前人基于該區構造演化、沉積充填特征及石油地質條件等方面的大量研究，總體認為形成于深湖—半深湖環境的沙河街組和東營組均發育烴源巖，是遼東灣油氣的主要來源［14，15］。但是，由于遼西凹陷位于遼東灣坳陷西緣，其周邊發育綏中、秦皇島、興城等水系，屬于近源快速沉積，對該凹陷沙河街組沙三段烴源巖的分布規模及品質認識不清楚，制約了遼西凹陷資源潛力的評價［16，17］。因此，遼西凹陷沙河街組優質烴源巖的評價與預測，是該凹陷油氣勘探面臨的首要問題。遼西凹陷已有覆蓋研究區的高品質三維地震資料和多口鉆井數據，具備利用地球物理技術評價烴源巖的數據基礎。

圖1 遼西凹陷構造區劃圖

3 研究區烴源巖的特征

3.1 地化特征

目前烴源巖的評價參數主要包括TOC、有機質類型、氯仿瀝青“A”含量、總烴含量等。在實際勘探中也有人提出了有效烴源巖和優質烴源巖的概念，前者常指經過排烴、運移、聚集等環節后仍具有工業價值的烴源巖，而后者為有機質含量高、類型好的有效烴源巖，通常它們的厚度不一定大，但卻具有較高的生烴潛力和排烴強度。

通過分析遼西凹陷及周邊遼中凹陷、遼東凹陷沙河街組的鉆井數據，可知遼東灣地區優質烴源巖主要出現在沙三下、沙三中以及沙一段。目前遼西凹陷內已有13口井鉆遇沙三段，揭示沙三段巖性主要為褐灰色泥巖，厚度較大。前人認為遼西凹陷沙三段烴源巖有較高含量的浮游藻類，有機屑含量較高，腐泥組和殼質組含量普遍大于80%，有機質類型主要為Ⅱ1型，埋深大，因此烴源巖處于成熟階段［16］。研究區內S-1井鉆遇沙三段地層650m，測井解釋泥巖厚度為540m，其中優質烴源巖為228m（TOC＞2%，生烴潛量S1+S2＞10mg／g，氫指數HI＞40mg／g）；S-2井鉆遇沙三段地層750m，測井解釋泥巖厚度為630m，其中優質烴源巖總厚度為340m（TOC＞2%，S1+S2＞10mg／g，HI＞40mg／g）。

本文根據中國近海湖相生烴凹陷評價指標體系［18］（表1）和遼西凹陷的鉆井資料樣品分析，指出TOC＞2%為本區優質烴源巖的主要判別指標。

表1 中國近海湖相烴源巖質量的評價指標（中海油，2015）

3.2 地震相特征

前人認為遼西凹陷沙河街組發育多種沉積相類型，主要有濱湖、半深湖—深湖、三角洲、扇三角洲和濁積扇等沉積體系。其中，凹陷內沙三段主要發育扇三角洲—中、深湖相沉積，中、深湖相面積較大，暗色泥巖分布廣、厚度大。沙三段這套半深湖—深湖相烴源巖表現為平行—亞平行反射結構，具有低頻、連續、強反 射 的 特 點［10，16，19，20］。 圖 2 為 研 究 區 內 過遼西凹陷北洼的一條地震剖面，從中可看出遼西凹陷北洼沙三段（E2s3）烴源巖位于洼陷中心，自下而上整體上表現出低頻、連續、平行的反射特征，其對應一套中、深湖相沉積。

3.3 測井相特征

湖相—中深湖相優質烴源巖的測井曲線表現為高自然伽馬、高聲波時差、高電阻率和低密度的特征。這是因為優質烴源巖往往含有更多的放射性元素鈾和干酪根及孔隙中的液態烴類。其中導電性差的干酪根和油氣的聲波時差一般大于巖石骨架的時差，而干酪根的密度（1.03～1.1g／cm3）和油的密度（0.7～0.98g／cm3）又小于巖石骨架密度（2.3～3.1g／cm3）［21］。表2為遼西凹陷沙河街組烴源巖的測井響應特征對比。測井揭示S-1井在沙三段鉆遇的優質成熟烴源巖（TOC＞2%）與其他品質不佳的泥巖相比（TOC＜1%），測井響應表現為典型的“三高一低”特征。

圖2 遼西北洼E2 s3烴源巖地震反射剖面

表2 遼西凹陷沙河街組烴源巖測井響應特征

4 方法應用

在綜合分析遼西凹陷沙河街組烴源巖地化特征、測井響應特征和地震反射特征的基礎上，結合地化分析和測井解釋等資料，制定了本區優質烴源巖的定量評價流程（圖3）。

圖3 優質烴源巖的定量評價流程

4.1 TOC計算

本文根據已鉆井的實測樣品TOC值，使用Passey經驗公式［22］

計算了靶區內井的TOC曲線。式中：ΔlgR=lg（R／R基線）+0.02（Δt-Δt基線），表示聲波—電阻率曲線幅度差，R基線和Δt基線分別為電阻率與聲波基線值；a=2.297-0.1688LOM，LOM為有機質成熟度的熱變指數。通過分析認為，遼西凹陷沙河街組地層電阻率基線值約為2.5Ω·m，聲波時差基線值約為92μs／ft，熱變指數LOM選8時，測井預測結果與實測值的誤差最小。圖4為S-1井烴源巖TOC實測值與測井計算值的關系圖，由圖可見二者相關性較高，說明由上述參數通過Passey公式建立的計算模型比較可靠。

圖4 S-1井計算TOC值與實測值的關系

4.2 巖石物理分析與敏感參數優選

通過巖石物理分析，能夠明確與TOC之間最為敏感的地球物理參數，進而建立兩者的擬合關系，最終根據關系式對TOC進行估算。對研究區內的S-1和S-2井進行交會分析。由圖5可見，縱橫波速度比為區分沙河街組砂泥巖的敏感參數，砂巖的縱橫波速度比小于2，而泥巖的縱橫波速度比則大于2，且在泥巖段內不同等級的烴源巖表現出不同范圍的縱波阻抗，據此就可形成識別該區泥巖和烴源巖等級的標準（圖5）。

圖5 遼西凹陷沙河街組巖石物理分析

4.3 地震反演TOC數據體

巖石物理與敏感參數優選結果顯示，通過地球物理技術是能夠對遼西凹陷烴源巖的分布范圍和等級進行預測的，即首先通過縱橫波速度比區分砂、泥巖，進而在泥巖之中根據縱波阻抗的范圍識別不同等級的烴源巖。

進一步通過縱波阻抗與TOC的交會分析（圖6），得到研究區內沙河街組TOC與PI（縱波阻抗）之間的計算模型為：TOC=-0.000015×PI+0.154，由此將反演的縱波阻抗轉換為TOC數據體。

需要說明的是，TOC曲線計算、巖石物理分析與敏感參數的優選主要是基于縱橫波速度、密度、電阻率等測井曲線及其組合，各測井曲線的品質和對目標的識別能力將會對后期的預測結果產生影響，因此準確的測井曲線校正工作是應用地球物理技術定量評價烴源巖的關鍵。對比圖4和圖6，可看出利用Passey公式（相關性：96%）和縱波阻抗（相關性：81%）計算的TOC的精度不同，經分析可能是由于緊貼井壁的密度曲線測井方式更容易受到井壁條件、孔內泥漿、孔壁泥餅等外界因素的影響，因此擬合精度較電阻率會低。

圖6 縱波阻抗（PI）與TOC的交會分析

另外受地質沉積環境、埋藏條件、巖性變化等因素的影響，不同凹陷不同層段的烴源巖與圍巖之間的巖石物理特征的差別變化很大，具有不同的TOC敏感識別參數、敏感參數范圍或識別相關性。如筆者在相關凹陷的研究中就發現渤海灣盆地沙南凹陷沙河街組的優質烴源巖的縱波阻抗小于8400g／cm3·m／s，而黃河口凹陷縱波阻抗（擬合相關性：73%）和縱橫波速度比（擬合相關性：82%）都可識別優質烴源巖，但后者的擬合相關性更高，因此黃河口凹陷的敏感參數為縱橫波速度比。同時，烴源巖與圍巖之間巖石物理參數的重疊交叉也會產生一定的多解性，因此在工作中需要針對各凹陷、各層段的特點分別開展分析工作，進而確定相應的敏感參數和計算模型，以期降低多解性。

工區內S-1井和S-2井參與了反演，S-3井為驗證井。圖7為通過疊前同時反演得到的過井vP／vS巖性剖面，其中測井曲線為泥質含量，右側為高值。由圖可見凹陷內沙二段為一套砂巖，沙三中、下段泥巖非常發育，井旁反演結果與實鉆吻合度高（合成記錄標定時相關系數可達0.93）。圖8為TOC數據體剖面，剖面中顏色表示有機質豐度，紅色表示TOC豐度較高，井曲線為TOC曲線，右側為高值。由圖看出沙三中、下段為一套TOC值較高的地層（＞2%），從洼陷邊緣到深部，這套烴源巖的厚度逐漸增大，反映出洼陷深部優質烴源巖更為發育。圖9為提取的沙三中、下段TOC預測平面圖，圖中顯示研究區中東部亮色部分烴源巖TOC大多為2%以上，屬于優質烴源巖發育帶，從凹陷深部到邊緣亮色逐漸減弱，表示地層的泥質含量和TOC逐漸降低，烴源巖等級為中等（1%＜TOC＜2%）或差等（0＜TOC＜1%），沉積相由深湖相—中深湖相轉換為三角洲相。

圖8 遼西凹陷沙三段TOC數據體剖面

4.4 烴源巖定量預測

本文基于Jason反演平臺Body Checking模塊完成了烴源巖的厚度預測。首先針對TOC數據體，根據選定的各級烴源巖的門檻值范圍，得到相應等級烴源巖體的屬性數據體；然后在垂向上累加屬性體的樣本點數，并將其與采樣間隔相乘作為時間厚度；最后通過時間厚度與速度相乘得到厚度值。

圖10為預測研究區優質烴源巖的分布圖，圖中顏色表示預測烴源巖的厚度，顏色越深表示厚度越大。從圖中可看出凹陷中部烴源巖最為發育，通過測量可知該處優質烴源巖的面積近180km2，平均厚度達200m，最厚處可達660m。對比研究區內實鉆井資料（表3），反演預測的烴源巖厚度與實鉆測井解釋結果基本一致，相對誤差較小。

基于本方法研究得到的優質烴源巖的面積、厚度和TOC等數據，結合行業內盆地資源潛力分析方法和靶區的勘探實際，對該區資源量的估算，采用了以下有機碳法

圖10 預測的研究區沙三段優質烴源巖厚度

表3 烴源巖預測與鉆井結果對比

式中：Q為總油氣資源量；H為烴源巖有效厚度（m）；S為烴源巖面積（km2）；ρ為烴源巖密度（g／cm3）；HI為氫指數（mg／g）；TOC為總有機碳含量；Kc為有機碳轉化系數；k為運聚系數。

本研究中，ρ取2.3g／cm3，HI取450mg／g，Kc取1%，k取10%，最終估算研究區沙三段油氣資源量約為2.6×109t。

5 結束語

應用地球物理技術可以對烴源巖的分布范圍進行預測，并能夠根據反演參數的定量關系區分烴源巖的等級，本文所述方案對類似生烴凹陷的評價具有參考價值。遼西凹陷沙河街組烴源巖評價的應用結果表明：

（1）該套優質烴源巖具有高自然伽馬、高聲波時差、高電阻率和低密度的測井響應特征；在地震相上表現為低頻、連續、強反射的地震反射特征。

（2）縱波阻抗參數是反映烴源巖豐度的敏感參數，根據TOC與縱波阻抗間的數學關系，可以通過地震反演實現烴源巖有機質含量的定量評價。

（3）研究區中部沙三段烴源巖TOC大多為2%以上，橫向整體分布穩定，平均厚度達200m。定量預測烴源巖厚度相對誤差小于10%，與實鉆對比吻合度較高，證實了應用地球物理方法定量預測烴源巖技術的有效性和適用性。

［1］ 顧禮敬，徐守余，蘇勁等.利用地震資料預測和評價烴源巖.天然氣地球科學，2011，22（3）：554-560.Gu Lijing，Xu Shouyu，Su Jin et al.Muddy hydrocarbon source rock prediction and evaluation with seismic data.Natural Gas Geoscience，2011，22（3）：554-560.

［2］ 王鵬，鐘建華，張寶權等.虛擬三維地震技術在綏濱拗陷油氣勘探中的應用.石油地球物理勘探，2011，46（增刊1）：97-101.Wang Peng，Zhong Jianhua，Zhang Baoquan et al.Hypothesized 3D seismic application in oil exploration of Suibin Depression.OGP，2011，46（S1）：97-101.

［3］ 王健，石萬忠，舒志國等.富有機質頁巖TOC含量的地球物理定量化預測.石油地球物理勘探，2016，51（3）：596-604.Wang Jian，Shi Wanzhong，Shu Zhiguo et al.TOC content quantitative prediction in organic-rich shale.OGP，2016，51（3）：596-604.

［4］ 李松峰，徐思煌，薛羅等.稀井區烴源巖有機碳的地球物理預測方法——珠江口盆地恩平凹陷恩平組烴源巖勘探實例.石油地球物理勘探，2014，49（2）：369-374.Li Songfeng，Xu Sihuang，Xue Luo et al.Source-rock organic carbon prediction with geophysical approach in the sparsely-drilled area：A case study of Enping Depression，the Pearl Mouth Basin.OGP，2014，49（2）：369-374.

［5］ 許杰，何治亮，董寧等.含氣頁巖有機碳含量地球物理預測.石油地球物理勘探，2013，48（增刊1）：64-68.Xu Jie，He Zhiliang，Dong Ning et al.Total organic carbon content prediction of gas-bearing shale with geophysical methods.OGP，2013，48（S1）：64-68.

［6］ 袁東山，王國斌，湯澤寧等.測井資料評價烴源巖方法及其進展.石油天然氣學報，2009，31（4）：192-203.Yuan Dongshan，Wang Guobin，Tang Zening et al.Methods for evaluation source rocks by well-logging data and its progress.Journal of Oil and Gas Technology，2009，31（4）：192-203.

［7］ 劉軍，汪瑞良，舒譽等.烴源巖TOC地球物理定量預測新技術及在珠江口盆地的應用.成都理工大學學報（自然科學版），2012，39（4）：415-419.Liu Jun，Wang Ruiliang，Shu Yu et al.Geophysical quantitative prediction technology about the total organic carbon in source rocks and application in Pearl River Mouth Basin，China.Journal of Chengdu University of Technology（Science＆ Technology Edition），2012，39（4）：415-419.

［8］ 郭旭升，尹正武，李金磊.海相頁巖含氣量地震定量預測技術及其應用——以四川盆地焦石壩地區為例.石油地球物理勘探，2015，50（1）：144-149.Guo Xusheng，Yin Zhengwu，Li Jinlei.Quantitative seismic prediction of marine shale gas content：a case study in Jiaoshiba area，Sichuan basin.OGP，2015，50（1）：144-149.

［9］ 蘇建龍，屈大鵬，陳超等.疊前地震反演方法對比分析——焦石壩頁巖氣藏勘探實例.石油地球物理勘探，2016，51（3）：581-588.Su Jianlong，Qu Dapeng，Chen Chao et al.Application of prestack inversion technique for shale gas in Jiaoshiba gas field.OGP，2016，51（3）：581-588.

［10］ 吳玉坤，胡明毅，柯嶺等.利用地震相識別優質烴源巖——以遼中凹陷沙三段為例.沉積學報，2013，31（2）：366-373.Wu Yukun，Hu Mingyi，Ke Ling et al.Identification of high-quality source rocks by seismic facies：taking the third member of Shahejie Formation in Liaozhong Depression as an example.Acta Sedimentologica Sinica，2013，31（2）：366-373.

［11］ 王貴文，朱振宇，朱廣宇.烴源巖測井識別與評價方法研究.石油勘探與開發，2002，29（4）：50-52.Wang Guiwen，Zhu Zhenyu，Zhu Guangyu.Logging identification and evaluation of Cambrian-Ordovician source rocks in syneclise of Tarim basin.Petroleum Exploration and Development，2002，29（4）：50-52.

［12］ 陶倩倩，李達，楊希冰等.利用分頻反演技術預測烴源巖.石油地球物理勘探，2015，50（4）：706-713，722.Tao Qianqian，Li Da，Yang Xibing et al.Hydrocarbon source rock prediction with frequency-divided inversion.OGP，2015，50（4）：706-713，722.

［13］ 黃艷輝，劉震，陳婕等.利用地震信息定量預測烴源巖熱成熟度——以瓊東南盆地樂東—陵水凹陷為例.石油地球物理勘探，2013，48（6）：985-994.Huang Yanhui，Liu Zhen，Chen Jie et al.Source rock maturity quantitative prediction using seismic information：a case study in the Ledong-Lingshui Depression，Qiongdongnan Basin.OGP，2013，48（6）：985-994.

［14］ 張功成，梁建設，徐建永等.中國近海潛在富烴凹陷評價方法與烴源巖識別.中國海上油氣，2013，25（1）：13-18.Zhang Gongcheng，Liang Jianshe，Xu Jianyong et al.An evaluation method of potential hydrocarbon-rich sags and their source rock identification.China Offshore Oil and Gas，2013，25（1）：13-18.

［15］ 蘆剛，林暢松，劉景彥等.遼西凹陷沙河街組低凸起帶物源分析.成都理工大學學報（自然科學版），2014，41（2）：171-176.Lu Gang，Lin Changsong，Liu Jingyan et al.Provenance analysis of low uplift of Shahejie Formation，Liaoxi sag，China.Journal of Chengdu University of Technology（Science＆ Technology Edition），2014，41（2）：171-176.

［16］ 梁建設，張功成，苗順德等.遼東灣遼西凹陷沙河街組烴源巖評價及油源研究.沉積學報，2012，30（4）：739-746.Liang Jianshe，Zhang Gongcheng，Miao Shunde et al.Evaluation of Shahejie source rock and oil source research in the Liaoxi Depression，Liaodong Bay，China.Acta Sedimentologica Sinica，2012，30（4）：739-746.

［17］ 朱筱敏，董艷蕾，楊俊生等.遼東灣地區古近系層序地層格架與沉積體系分布.中國科學：D輯，2008，38（增1）：1-10.Zhu Xiaomin，Dong Yanlei，Yang Junsheng et al.Sequence framework and sediment system of Eogene in Liaodong Bay area.Science in China（Series D），2008，38（S1）：1-10.

［18］ 中國海洋石油總公司.Q／HS 1068-2015中國近海湖相生烴凹陷評價指標體系.北京：中國海洋石油總公司企業標準，2015.

［19］ 張寒，朱光有.利用地震和測井信息預測和評價烴源巖——以渤海灣盆地富油凹陷為例.石油勘探與開發，2007，34 （1）：55-60.Zhang Han，Zhu Guangyou.Using seismic and log information to predict and evaluate hydrocarbon source rocks：An example from rich oil depressions in Bohai Bay.Petroleum Exploration and Development，2007，34（1）：55-60.

［20］ 鄧運華，張功成，劉春成等.中國近海兩個油氣帶地質理論與勘探實踐.北京：石油工業出版社，2013.

［21］ 于翔濤.測井技術在烴源巖評價中的應用.長江大學學報（自然科學版），2009，6（2）：198-201.Yu Xiangtao.Application of well logging technology to source rocks evaluation.Journal of Yangtze University（Natural Science Edition），2009，6（2）：198-201.

［22］ 朱振宇，劉洪，李幼銘.ΔlogR技術在烴源巖識別中的應用與分析.地球物理學進展，2003，18（4）：647-649.Zhu Zhenyu，Liu Hong，Li Youming.The analysis and application ofΔlogRmethod in the source rock’s identification.Progress in Geophysics，2003，18（4）：647-649.

P631

A

10.13810／j.cnki.issn.1000-7210.2017.01.018

牛聰，劉志斌，王彥春，張益明，解吉高，黃饒.應用地球物理技術定量評價遼西凹陷沙河街組烴源巖.石油地球物理勘探，2017，52（1）：131-137.

1000-7210（2017）01-0131-07

*北京市朝陽區太陽宮南街6號院A-1003室，100028。Email：niucong＠cnooc.com.cn

本文于2016年1月22日收到，最終修改稿于同年10月29日收到。

（本文編輯：馮杏芝）

牛聰 高級工程師，博士研究生，1980年生；2006年畢業于成都理工大學，獲信號與信息處理專業碩士學位，2013年起在中國地質大學（北京）攻讀地球探測與信息技術專業博士學位；現在中海油研究總院主要從事油氣地球物理勘探技術研究。