鄂爾多斯盆地致密氣儲集層孔隙流體地震預測方法

楊華，王大興，張盟勃，王永剛，劉力輝，張盟黎

（1. 中國石油天然氣集團公司，北京 100007；2. 低滲透油氣田勘探開發(fā)國家工程實驗室，西安 710018；3. 中國石油長慶油田公司勘探開發(fā)研究院，西安 710018；4. 成都理工大學，成都 610059）

鄂爾多斯盆地致密氣儲集層孔隙流體地震預測方法

楊華1，王大興2,3，張盟勃2,3，王永剛3，劉力輝4，張盟黎3

（1. 中國石油天然氣集團公司，北京 100007；2. 低滲透油氣田勘探開發(fā)國家工程實驗室，西安 710018；3. 中國石油長慶油田公司勘探開發(fā)研究院，西安 710018；4. 成都理工大學，成都 610059）

基于鄂爾多斯盆地上古生界二疊系下石盒子組盒 8段致密儲集層巖石物理分析結果，考慮體積模量對孔隙流體敏感的特點，提出了致密儲集層巖相控制下利用壓縮系數(shù)檢測流體性質的方法。該方法首先是進行巖相識別，利用疊前反演的縱橫波交會技術計算致密砂巖相的分布；其次是在巖相控制下（排除泥質含量影響），利用可穩(wěn)定反演的波阻抗和速度計算壓縮系數(shù)；最后利用氣水層壓縮系數(shù)的差異進行儲集層孔隙流體性質判識。蘇里格致密氣田勘探開發(fā)生產(chǎn)實踐表明：應用壓縮系數(shù)進行孔隙流體性質識別的效果較好，能有效預測氣水層的分布。圖 9表2參19

鄂爾多斯盆地；二疊系下石盒子組；致密儲集層；彈性參數(shù)；壓縮系數(shù)；孔隙流體預測

0 引言

鄂爾多斯盆地致密儲集層與圍巖波阻抗差異小，氣層厚度一般較薄，常規(guī)地震預測難度大，利用單一地球物理參數(shù)特別是地震疊后參數(shù)無法預測儲集層。地球物理理論技術發(fā)展經(jīng)歷了聲學和波動學理論階段、彈性和黏彈性理論階段以及各向異性巖石力學理論階段，相應儲集層預測的地球物理參數(shù)也由地震疊后的“三瞬”（瞬時振幅、瞬時頻率、瞬時相位）屬性、波速和波阻抗等，發(fā)展為疊前反演的泊松比、拉梅系數(shù)、黏滯系數(shù)、吸收衰減系數(shù)（地震品質因子Q）、彈性阻抗和彈性模量等彈性參數(shù)，以及壓縮系數(shù)、脆性指數(shù)、兩向應力差、抗壓強度等巖石力學參數(shù)[1-2]。近年來，隨著常規(guī)油氣資源的發(fā)現(xiàn)日益減少，勘探和開發(fā)非常規(guī)油氣資源成為了必然選擇，其中非常規(guī)油氣資源的甜點預測成為了研究熱點，國內非常規(guī)儲集層甜點區(qū)評價研究中的甜點包括地質甜點和工程甜點。地質甜點著眼于烴源巖、儲集層與裂縫等綜合評價；工程甜點著眼于巖石可壓性、地應力各向異性綜合評價[1]。由于致密儲集層的巖石力學特征對其有利區(qū)預測和水平井甜點靶區(qū)優(yōu)選意義重大，因此，巖石力學參數(shù)成為致密儲集層預測的關鍵。

按照孔隙彈性介質理論[2]，彈性三模量（體積模量、切變模量、彈性模量）和泊松比等是儲集層研究重要的彈性參數(shù)，這些參數(shù)中孔隙流體的體積模量對儲集層特別是氣層較為敏感。儲集層骨架一定時（比如都為砂巖），體積模量主要受孔隙流體體積模量影響，因此，應用儲集層的體積模量進行氣水識別也相對有效。儲集層體積模量的倒數(shù)即為儲集層壓縮系數(shù)（以下簡稱壓縮系數(shù)），可反映孔隙流體氣液兩相的差異，是評價儲集層流體特性的關鍵指標[3]。以往壓縮系數(shù)預測主要應用測井巖石物理分析和偶極子縱橫波測井資料計算，其準確性受地震密度反演穩(wěn)定性及精度的影響較大，利用地震疊前資料預測壓縮系數(shù)方面的研究與應用仍處于探索階段。

鄂爾多斯盆地上古生界二疊系下石盒子組以河流—三角洲相沉積為主，發(fā)育水下分流河道和辮狀河道微相，砂體致密、氣層厚度薄，為典型的致密砂巖氣[3-4]。本文對鄂爾多斯盆地蘇里格氣田下石盒子組盒 8段致密氣儲集層的巖石力學特征、地震預測方法進行了探索，提出利用多波地震資料或高品質的縱波疊前地震資料進行相控預測壓縮系數(shù)的方法，該方法要求先排除泥質含量影響，其次在致密砂巖相計算過程中，密度參數(shù)不直接參與計算，而用穩(wěn)定反演的波阻抗和速度參數(shù)來計算壓縮系數(shù)，然后將壓縮系數(shù)用于致密氣儲集層的氣水識別。

1 壓縮系數(shù)的地震預測方法

根據(jù)流體混合（氣、水）模型，一般情況下純氣和純水的體積模量相差約2.5 GPa[3]。鄂爾多斯盆地蘇里格氣田盒 8段沉積期大型辮狀河道砂體沉積呈廣覆式分布，普遍含氣，儲集層非均質性強。蘇里格氣田開發(fā)中區(qū)平均儲量豐度為0.44×108m3/km2，高含氣區(qū)豐度平均為 0.88×108m3/km2，低含氣區(qū)豐度平均為0.28×108m3/km2[5]，含氣飽和度變化大。在蘇里格氣田盒 8段致密砂巖中，相對高孔滲砂巖含氣，部分低孔滲砂巖含水，氣水分異不明顯，以毛細管水和滯留水為主，其余為少量原狀束縛水（原始地層狀態(tài)下難流動，壓裂后產(chǎn)水）；部分特低孔滲砂巖為不含氣水的干砂巖，成為巖性遮擋層。通過對歷年43口井偶極聲波測井資料重新處理校正，測井巖石物理分析表明，蘇里格氣田致密砂巖含流體（氣、水）層的地球物理屬性差異相對較小[6]，就含氣砂巖和含水砂巖而言，其橫波速度基本無變化，縱波速度差異在 100～400 m/s，縱橫波速度比差異在 0.1～0.2，泊松比差異在 0.06～0.10。在區(qū)分巖性的基礎上，彈性“三模量”中的體積模量預測流體的物理意義比較明確，檢測窗口也相對較大（見表 1），致密含氣砂巖的體積模量與氣水同層砂巖的體積模量相差0.9～2.1 GPa。

表1 蘇里格氣田上古生界地層測井和巖石力學參數(shù)統(tǒng)計表

當儲集層孔隙內充滿氣時，其縱波能量衰減顯著增大[6]，泊松比降低，體積模量顯著減小導致其壓縮系數(shù)升高，因為氣層的壓縮系數(shù)要比相同巖性相同孔隙度的水層大。在蘇里格氣田盒 8段儲集層中高含氣飽和度（含氣飽和度大于 60%）段（見圖 1），壓縮系數(shù)對含氣飽和度的變化最敏感。蘇里格氣田西部盒 8段19口井34個層點統(tǒng)計表明（見圖1），氣層壓縮系數(shù)大于1.68 GPa-1，而干層大多小于1.0 GPa-1，氣水同層大多介于兩者之間。因此，由體積模量分析出發(fā)，計算壓縮系數(shù)和泊松比等巖石力學參數(shù)，以及縱橫波聯(lián)合反演及交會解釋是預測氣水分布的有效方法。

圖1 不同類型儲集層壓縮系數(shù)與含水飽和度關系圖（19口井34個層點）

1.1 巖相巖性預測思路

利用實際地震資料計算泥質含量和彈性參數(shù)時，經(jīng)常用到反演的縱、橫波速度（波阻抗）等，由于反演結果存在一定誤差，因此很可能導致計算的彈性參數(shù)誤差較大。針對此，提出了用巖相控制預測巖石力學參數(shù)的思想，具體做法是：①測井上綜合若干相關測井曲線特征劃分巖相，將地層劃分為砂巖、泥巖等不同巖性；②用坐標旋轉的方法根據(jù)井上巖相曲線計算出地震巖相，使得無井區(qū)也能有巖性信息；③在巖相控制下計算泥質含量、彈性模量等參數(shù)，例如，頁巖的泥質含量為 100%，在計算時可以把頁巖定義為非儲集層的泥巖，這種類型的巖相不需用公式來計算泥質含量，直接給定數(shù)值為100%，或者不同的巖相若擬合出不同的公式，用不同的公式來分巖相計算參數(shù)。

1.2 壓縮系數(shù)的計算方法

壓縮系數(shù)反映了巖石的總體可壓縮性，定義為體積應變與流體靜壓力之比，與抗壓強度意義正好相反，其計算無需考慮泥質含量，壓縮系數(shù)還可以定義為體積模量的倒數(shù)：

體積模量的計算公式為：

由于密度獲取受制于地震疊前資料品質和反演方法的精度，密度值誤差較大，而且求解的穩(wěn)定性也差，因此，一般利用波阻抗和速度這些可以穩(wěn)定反演的參數(shù)來代替密度，計算公式如下：

雖然壓縮系數(shù)的計算無需考慮泥質含量，但它和抗壓強度為同一性質的參數(shù)，為明確其巖石物理意義，也需要在巖相控制下在砂巖區(qū)計算，具體流程見圖2。

1.3 壓縮系數(shù)識別氣層的邊界條件分析

根據(jù)蘇里格氣田盒8段、山1段儲集層的西區(qū)12口井31個數(shù)據(jù)點、中區(qū)9口井26個數(shù)據(jù)點和東區(qū)9口井26個數(shù)據(jù)點的壓縮系數(shù)與泊松比交會關系，將其儲集層流體類型的參數(shù)界限值分析結果分別列于表2。以東區(qū)圖版分析為例（見圖 3），氣層壓縮系數(shù)大于1.64 GPa-1，泊松比小于 0.18；差氣層壓縮系數(shù)小于1.64 GPa-1，泊松比大于0.18；干層壓縮系數(shù)小于1.12 GPa-1，泊松比大于0.20；水層、氣水同層壓縮系數(shù)小于1.25 GPa-1，泊松比大于0.24。

圖2 相控壓縮系數(shù)計算流程

表2 蘇里格氣田不同類型儲集層泊松比、壓縮系數(shù)參數(shù)界限

圖3 蘇里格氣田東區(qū)儲集層壓縮系數(shù)與泊松比交會圖（9口井26個數(shù)據(jù)點）

1.4 巖相識別方法

鄂爾多斯盆地不論是常規(guī)儲集層還是非常規(guī)致密儲集層，沉積相（帶）決定巖（性）相，巖相決定物性（相），物性（相）決定含油氣性和甜點分布。致密儲集層物性預測難度極大，特別是對于滲透率約為0.1×10-3μm2的儲集層，地震難以預測，但綜合孔隙度與飽和度描述的致密儲集層物性（相）與非儲集層就有一定的差異。劉力輝等[7]把在地震上能識別的物性（相）定義為地震物相，即與儲集層物性相關的、地震彈性參數(shù)可分辨的屬性類別，根據(jù)巖石的微觀結構，綜合巖性（沉積成因）、物性（壓實與成巖成因）和地震彈性參數(shù)的可識別性來劃分。文獻[7]詳述了地震物相的計算方法，計算巖相比計算物相簡單，只需考慮巖性，不用考慮巖石的物性信息。

1.4.1 測井劃分巖相

綜合相關測井曲線信息，先在測井資料上人為解釋出幾類巖相，如圖 4中解釋的黃色為砂巖相，黑色為泥巖、油頁巖和煤層等，統(tǒng)一歸為非砂巖相。

圖4 測井巖相劃分示意圖

1.4.2 坐標旋轉計算地震巖相體

坐標旋轉技術在數(shù)學上就是一種映射，即

由公式（4），可將兩個參數(shù)線性組合形成一個新參數(shù)，通過尋找合適的旋轉角度，使生成的新參數(shù)對地震物相有最大的區(qū)分能力。目前的坐標轉換技術是一種線性組合技術，是針對地震信息同時包含巖性和儲集層物性信息、單一彈性參數(shù)難以直接準確描述而提出的技術，其技術思路為：從地震疊前反演的井旁道抽取彈性波阻抗曲線，比如大小角度射線彈性波阻抗，通過旋轉一定角度，生成一個新的坐標系統(tǒng)，在此坐標系統(tǒng)下將兩個不同的彈性參數(shù)，通過線性組合生成一個新的復合參數(shù)，以便能最佳地區(qū)分各類地震巖相體。

由圖5可見，單一的彈性參數(shù)識別巖相能力有限，通過坐標旋轉形成新的復合參數(shù)體后可有效區(qū)分不同巖相體（見圖6）。坐標旋轉要求兩個參數(shù)的量綱一致，因此，選擇大小角度射線彈性波阻抗[8-9]或大小角度縱波波阻抗和橫波波阻抗來進行。其次，用于地震巖相計算的彈性波阻抗反演方法最好選用模型約束少的反演方法，使反演出的彈性參數(shù)有反映地質體邊界的能力[9-10]。這樣要求的目的是能夠用地震巖相體的地層切片信息，通過平面地震地貌和沉積模式對比[10-11]，驗證地震巖相體計算的合理性。

圖5 測井大小角度彈性波阻抗交會圖

圖6 坐標旋轉后交會圖

根據(jù)新參數(shù)的旋轉公式，用地震彈性參數(shù)和旋轉后劃分巖相的新門檻值，計算出地震巖相體（見圖 7）。

圖7 蘇里格氣田地震疊前反演縱波波阻抗（a）、橫波波阻抗（b）、泊松比（c）、巖相（d）和壓縮系數(shù)剖面（e）

2 應用實例

在諸多的巖石力學參數(shù)中，壓縮系數(shù)是研究致密儲集層很重要的參數(shù)[12]。實際應用結果與鉆井資料以及地質情況對比表明，在巖相控制下計算的壓縮系數(shù)能夠指示致密儲集層的甜點發(fā)育區(qū)。壓縮系數(shù)對含氣飽和度的敏感性明顯優(yōu)于泊松比、彈性模量和衰減等參數(shù)[13]，壓縮系數(shù)是刻畫致密含氣儲集層甜點的重要參數(shù)。在評價蘇里格氣田盒8段致密砂巖含氣性時，壓縮系數(shù)按含氣層、差氣層、含水層、干層順序逐漸變小，因此，壓縮系數(shù)在蘇里格氣田盒8段氣層的研究中具有顯著意義。

2.1 工區(qū)概況

蘇里格氣田為典型大型巖性氣藏，二疊系下石盒子組盒8段分為上、下兩段，儲集層以盒8下段水下分流河道沉積的巖屑質石英砂巖為主，盒 8段砂巖的孔隙度為4%～14%，平均孔隙度為8.8%；滲透率主要為（0.05～5.00）×10-3μm2，平均滲透率為 0.872×10-3μm2[3]，為典型致密氣田[14-16]。但在大面積致密儲集層背景下局部發(fā)育相對高滲儲集層，以水下分流河道中部砂巖為主，巖性多為粗砂巖，分布不連續(xù)，是儲量和產(chǎn)能的甜點區(qū)。儲集空間以殘余粒間孔、溶蝕孔、晶間孔為主，甜點區(qū)為中高含氣飽和度（55%～85%），發(fā)育微米級和亞微米—納米級兩大孔喉體系，儲集層滲流能力與喉道大小關系分析表明，占喉道總量不足25%的大孔喉對滲透率的貢獻值達到了 80%。因此，尋找相對高孔滲區(qū)和氣水分布預測是制約氣田有效開發(fā)的關鍵[17-18]。測井巖石物理分析表明，壓縮系數(shù)對含氣飽和度的變化最敏感，可顯著反映氣液兩相差異，壓縮系數(shù)可由疊前反演縱橫波彈性參數(shù)求取，是預測致密儲集層甜點的重要參數(shù)。

2.2 應用實例

2.2.1 鉆后方法檢驗

蘇里格地區(qū)全數(shù)字多波地震勘探始于2003年[3,19]，截至2010年共完成全數(shù)字三分量三維地震103 km2、全數(shù)字三分量二維地震5 360 km。采用單點三分量數(shù)字檢波器接收，單深井激發(fā)，覆蓋次數(shù)100，觀測系統(tǒng)為 6035-45-10-45-6035。由于在低降速帶較?。ㄐ∮?0 m）的沙漠區(qū)勘探，3個分偏移距疊加的覆蓋次數(shù)達到30次以上，資料信噪比高[18]，因此通過縱橫波聯(lián)合反演獲得了高品質的縱橫波速度和波阻抗剖面。按圖2壓縮系數(shù)的計算流程，在地震巖相預測剖面的基礎上，計算壓縮系數(shù)，結果如圖7e所示。

圖7剖面有8口鉆井，S1、S2和S8井試氣井口產(chǎn)量分別 14.7×104、4.5×104和 9.8×104m3/d，為高產(chǎn)工業(yè)氣井；S3、S5和S7井試氣井口產(chǎn)量僅為0.8×104、1.4×104和 1.5×104m3/d，為低產(chǎn)氣井；而S4和 S6井則為產(chǎn)水氣井，S4井產(chǎn)氣3.1×104m3/d、產(chǎn)水4.5 m3/d，S6井產(chǎn)氣2.0×104m3/d、產(chǎn)水7.5 m3/d。

鉆后檢驗結果表明：高含氣儲集層表現(xiàn)為高橫波波阻抗（紅—黃色）、低泊松比（紅—黃色）和高壓縮系數(shù)（紅—黃色），如S1、S2和S8井井旁反演剖面；差含氣儲集層表現(xiàn)為高橫波波阻抗（紅—黃色）、中—低泊松比（黃色）和中壓縮系數(shù)（黃色），如 S3、S5和S7井井旁反演剖面；氣水同層表現(xiàn)為高橫波波阻抗（紅—黃色）、中泊松比（黃色）和低壓縮系數(shù)（綠—藍色），如產(chǎn)水量大的 S4、S6井井旁反演剖面的中高壓縮系數(shù)層太薄，因而該井壓裂后氣水同出，壓縮系數(shù)預測結果與實際試氣結果符合較好。

2.2.2 鉆前流體檢測

探井E01井位于蘇里格氣田北部，鄰井鉆探、試氣資料表明，本區(qū)盒 8段砂巖儲集層較發(fā)育，含氣普遍，是勘探的有利目標區(qū)。同時因沉積相帶變化而導致砂體變薄或儲集層物性變差、氣層厚度變薄，且鄰井不同程度產(chǎn)水，因而存在一定地質風險。地震資料顯示（見圖8），該井處盒8段疊前反演橫波波阻抗高（紅色）、泊松比低（紅—黃色），表明該井處砂體發(fā)育（預測砂層厚30 m）且存在氣層，相控壓縮系數(shù)剖面顯示盒8段下部氣層發(fā)育、盒8段上部含水。實鉆該井盒8段鉆遇砂巖35.9 m，測井解釋盒8段下部氣層10.1 m、氣水同層2.5 m，盒8段上部含氣水層15 m，試氣獲無阻流量10.119 5×104m3/d、產(chǎn)水7.5 m3/d，地震預測與實鉆結果吻合。

圖8 蘇里格氣田E01井地震預測橫波波阻抗（a）、泊松比（b）、巖相（c）、壓縮系數(shù)（d）和盒8段綜合測井圖（e）

近期部署完鉆的 E02井地震預測剖面顯示（見圖9），該井處盒8段疊前反演橫波波阻抗高（紅色）、泊松比低（紅—黃色），表明該井處砂體發(fā)育（預測砂層厚30 m）且存在氣層，相控壓縮系數(shù)剖面顯示盒 8段下部氣層發(fā)育、盒 8段上部含水。實鉆該井盒8段鉆遇砂巖32.3 m，測井解釋盒8段下部氣層4.8 m，盒8段上部氣層5.3 m（待試），地震預測與實鉆結果吻合。

2.3 壓縮系數(shù)流體識別的應用條件

基于地震資料估算儲集層壓縮系數(shù)進行孔隙流體特征（含氣性）預測是一個新的探索，在地質應用上首先要考慮儲集層壓縮系數(shù)與流體飽和度之間是否存有一定關系，這需要做大量的巖石物理分析；其次，要研究巖石物理建模的流體替換分析結論，對含氣性的敏感程度進行評價，此外，在地球物理上要考慮儲集層與非儲集層區(qū)分的可能性和可信度；最后，要有滿足疊前彈性反演的高品質地震分偏移距疊加數(shù)據(jù)，最佳應用條件是工區(qū)具備地震三分量勘探獲得的高品質縱橫波數(shù)據(jù)。上述應用范圍和適應性是利用疊前地震資料預測致密儲集層孔隙流體性質的必要條件，既要具備儲集層地質應用條件，又要具備地震疊前反演彈性參數(shù)的高品質資料基礎。

圖9 蘇里格氣田E02井地震預測橫波波阻抗（a）、泊松比（b）和壓縮系數(shù)（c）剖面

3 結論

通過蘇里格氣田下石盒子組盒 8段致密儲集層大量的巖石物理分析，提出了相控壓縮系數(shù)預測孔隙流體特征的方法，即在地層巖相識別的基礎上，利用孔隙氣液兩相壓縮系數(shù)的差異進行儲集層流體性質識別；利用地震疊前資料可穩(wěn)定反演的波阻抗和速度參與計算，提高了預測的穩(wěn)定性和精度。該方法應用于蘇里格氣田下石盒子組盒 8段致密儲集層氣層、氣水同層、水層識別，與實鉆結果吻合較好，并在蘇里格致密氣田甜點預測、流體識別、實際井位部署中取得了較好的預測效果。壓縮系數(shù)的獲取受制于原始地震資料的質量和地震分偏移距疊加成像的品質，需要引起研究者的足夠重視。

符號注釋：

Cc——壓縮系數(shù)，GPa-1；K——體積模量，GPa；m1，m2，n1，n2——坐標變換系數(shù)，無因次；vp——縱波速度，m/s；vs——橫波速度，m/s；x，y——原坐標系；X，Y——新坐標系；Zp——縱波波阻抗，kg/（m2·s）；Zs——橫波波阻抗，kg/（m2·s）；ρ——密度，kg/m3。

[1] 鄒才能, 張國生, 楊智, 等. 非常規(guī)油氣概念、特征、潛力及技術：兼論非常規(guī)油氣地質學[J]. 石油勘探與開發(fā), 2013, 40(4)：385-399.ZOU Caineng, ZHANG Guosheng, YANG Zhi, et al. Geological concepts, characteristics, resource potential and key techniques of unconventional hydrocarbon： On unconventional petroleum geology[J]. Petroleum Exploration and Development, 2013, 40(4)：385-399.

[2] 黃緒德, 李明. 孔隙彈性與油氣預測[M]. 北京： 石油工業(yè)出版社,2012.HUANG Xude, LI Ming. Pore elasticity and oil and gas prediction[M]. Beijing： Petroleum Industry Press, 2012.

[3] 楊華, 王喜雙, 王大興, 等. 蘇里格氣田多波地震勘探關鍵技術[M]. 北京： 石油工業(yè)出版社, 2013.YANG Hua, WANG Xishuang, WANG Daxing, et al. Critical technology for multi-wave seismic exploration in the Sulige Gasfield[M]. Beijing： Petroleum Industry Press, 2013.

[4] 段長江, 羅順社, 周彪. 蘇里格地區(qū)中二疊統(tǒng)下石盒子組盒 8段沉積相研究[J]. 沉積與特提斯地質, 2012, 32(1)： 55-60.DUAN Changjiang, LUO Shunshe, ZHOU Biao. Sedimentary facies in the eighth member of the Middle Permian Lower Shihezi Formation in the Sulige region[J]. Sedimentary Geology and Tethyan Geology, 2012, 32(1)： 55-60.

[5] 盧濤, 劉艷俠, 武力超, 等. 鄂爾多斯盆地蘇里格氣田致密砂巖氣藏穩(wěn)產(chǎn)難點與對策[J]. 天然氣工業(yè), 2015, 35(6)： 43-52.LU Tao, LIU Yanxia, WU Lichao, et al. Challenges to and countermeasures for the production stabilization of tight sandstone gas reservoirs of the Sulige Gasfield, Ordos Basin[J]. Natural Gas Industry, 2015, 35(6)： 43-52.

[6] 王大興, 辛可鋒, 李幼銘, 等. 地層條件下砂巖含水飽和度對波速及衰減影響的實驗研究[J]. 地球物理學報, 2006, 49(3)：908-914.WANG Daxing, XIN Kefeng, LI Youming, et al. An experimental study of influence of water saturation on velocity and attenuation in sandstone under stratum conditions[J]. Geophysics, 2006, 49(3)：908-914.

[7] 劉力輝, 李建海, 劉玉霞. 地震物相分析方法與“甜點”預測[J].石油物探, 2013, 52(4)： 432-437.LIU Lihui, LI Jianhai, LIU Yuxia. Seismic reservoir property facies analysis and sweet spot prediction[J]. Geophysical Prospecting for Petroleum, 2013, 52(4)： 432-437.

[8] 馬勁風. 地震勘探中廣義彈性阻抗的正反演[J]. 地球物理學報,2003, 46(1)： 118-124.MA Jinfeng. Forward modeling and inversion method of generalized elastic impedance in seismic exploration[J]. Geophysics, 2003, 46(1)：118-124.

[9] 劉力輝, 王緒本. 一種改進的射線彈性阻抗公式及彈性參數(shù)反演[J]. 石油物探, 2011, 50(4)： 155-159.LIU Lihui, WANG Xuben. Elastic impedance inversion based on an improved ray-path elastic impedance formula[J]. Geophysical Prospecting for Petroleum, 2011, 50(4)： 155-159.

[10] 劉力輝, 王緒本, 陳雪菲．地震地貌學在L區(qū)濁積砂體預測中的應用[J]. 石油物探, 2011, 52(2)： 99-102.LIU Lihui, WANG Xuben, CHEN Xuefei. Application of seismic geomorphology on tuibidite in L area[J]. Geophysical Prospecting for Petroleum, 2011, 52(2)： 99-102.

[11] 劉保國, 劉力輝．實用地震沉積學在沉積相分析中的應用[J]. 石油物探, 2008, 47(3)： 266-271.LIU Baoguo, LIU Lihui. Application of practical seismic sedimentary in depositional facies analysis[J]. Geophysical Prospecting for Petroleum, 2008, 47(3)： 266-271.

[12] 王大興, 劉力輝, 王永剛. 致密儲層巖石力學參數(shù)預測方法研究及應用[R]. 上海： 第三屆地震巖石物理學術研討會, 2014.WANG Daxing, LIU Lihui, WANG Yonggang. The study and application of the prediction method for the prediction method of the tight reservoir rock mechanics parameter[R]. Shanghai： 3rd Seismic Rock Physics Academic Seminar, 2014.

[13] 楊雙定. 縱橫波參數(shù)提取與偶極橫波老井資料重新處理[R]. 西安：中國石油集團測井有限公司長慶事業(yè)部, 2013.YANG Shuangding. Primary wave and shear wave parameters extraction and dipole shear wave data of old well to treatment[R].Xi’an： Changqing Bureau, Well Logging Cooperation, CNPC, 2013.

[14] 楊華, 付金華, 劉新社, 等. 蘇里格大型致密砂巖氣藏形成條件及勘探技術[J]. 石油學報, 2012, 33(S1)： 28-35.YANG Hua, FU Jinhua, LIU Xinshe, et al. Formation conditions and exploration technology of large-scale tight sandstone gas reservoir in Sulige[J]. Acta Petrolei Sinica, 2012, 33(S1)： 28-35.

[15] 楊華, 席勝利, 魏新善, 等. 蘇里格地區(qū)天然氣勘探潛力分析[J].天然氣工業(yè), 2006, 26(12)： 45-48.YANG Hua, XI Shengli, WEI Xinshan, et al. Analysis on gas exploration potential in Sulige area of the Ordos Basin[J]. Natural Gas Industry, 2006, 26(12)： 45-48.

[16] 楊華, 劉新社, 楊勇. 鄂爾多斯盆地致密氣勘探開發(fā)形勢與未來發(fā)展展望[J]. 中國工程科學, 2012, 14(6)： 40-48.YANG Hua, LIU Xinshe, YANG Yong. Status and prospects of tight gas exploration and development in the Ordos Basin[J]. Engineering Sciences, 2012, 14(6)： 40-48.

[17] 羅順社, 彭宇慧, 魏新善, 等. 蘇里格氣田致密砂巖氣水相滲曲線特征與分類[J]. 西安石油大學學報(自然科學版), 2015, 30(6)：55-61.LUO Shunshe, PENG Yuhui, WEI Xinshan, et al. Characteristics and classification of gas-water relative permeability curves of tight sandstone reservoirs in Sulige gas field[J]. Journal of Xi’an Shiyou University (Natural Science Edition), 2015, 30(6)： 55-61.

[18] 郝騫, 張吉, 路中奇, 等. 蘇里格氣田西區(qū)蘇 48區(qū)塊儲層微觀特征及其對儲層的影響作用[J]. 天然氣地球科學, 2012, 23(6)：1011-1018.HAO Qian, ZHANG Ji, LU Zhongqi, et al. Microscopic characteristic and its influences on reservoirs from the Su 48 Block in the West Sulige gas field[J]. Natural Gas Geoscience, 2012, 23(6)：1011-1018.

[19] 王大興, 張盟勃, 陳娟, 等. 多波地震技術在致密砂巖氣藏勘探中的應用[J]. 石油地球物理勘探, 2014, 49(5)： 946-953, 1005.WANG Daxing, ZHANG Mengbo, CHEN Juan, et al. Application of multi-wave seismic in tight sandstone gas reservoir exploration[J].Oil Geophysical Prospecting, 2014, 49(5)： 946-953, 1005.

（編輯 黃昌武）

Seismic prediction method of pore fluid in tight gas reservoirs,Ordos Basin, NW China

YANG Hua1, WANG Daxing2,3, ZHANG Mengbo2,3, WANG Yonggang3, LIU Lihui4, ZHANG Mengli3
(1. China National Petroleum Corporation, Beijing 100007, China; 2. National Engineering Laboratory for Exploration and Development of Low-Permeability Oil & Gas Fields, Xi’an 710018, China; 3. Exploration and Development Research Institute of PetroChina Changqing Oilfield Company, Xi’an 710018, China; 4. Chengdu University of Technology, Chengdu 610059, China)

Based on the rock physical analysis of tight gas reservoirs in He8 Member of Permian Lower Shihezi Formation of Ordos Basin, and considering that bulk modulus is sensitive to pore fluid, this study proposes a fluid property detection method by compression coefficient of tight sandstone reservoirs under the constraint of reservoir lithofacies. In this method, lithofacies is identified first by calculating distribution of tight sandstone facies with the cross plot of Vp and Vs obtained from pre-stack seismic inversion; secondly, the compression coefficient is calculated by P-wave impendance and velocity from stable pre-stack seismic inversion with the restriction of lithofacies (excluding the influence of clay content); and finally, pore fluid properties are determined using the differences of compression coefficients in gas and water layers. Its application in tight gas exploration and development in Sulige gas field of Ordos Basin shows that this pore fluid prediction method by calculating compression coefficient can effectively and efficiently delineate the distribution of gas-bearing and water-bearing sandstone.

Ordos Basin; Permian Lower Shihezi Formation; tight reservoir; elastic parameters; compression coefficient; pore fluid prediction

國家自然科學基金項目“非常規(guī)油氣勘探與開發(fā)地球物理基礎理論與方法研究”（41390451）；國家科技重大專項“鄂爾多斯盆地大型低滲透巖性地層油氣藏開發(fā)示范工程”（2016ZX05050）

TE122.2

A

1000-0747(2017)04-0513-08

10.11698/PED.2017.04.04

楊華, 王大興, 張盟勃, 等. 鄂爾多斯盆地致密氣儲集層孔隙流體地震預測方法[J]. 石油勘探與開發(fā), 2017,44(4)： 513-520.

YANG Hua, WANG Daxing, ZHANG Mengbo, et al. Seismic prediction method of pore fluid in tight gas reservoirs, Ordos Basin,NW China[J]. Petroleum Exploration and Development, 2017, 44(4)： 513-520.

楊華（1963-），男，山東菏澤人，博士，中國石油天然氣集團公司教授級高級工程師，主要從事石油地質綜合研究及油氣勘探開發(fā)管理工作。地址：北京市東城區(qū)東直門北大街 9號，中國石油天然氣集團公司，郵政編碼：100007。E-mail：yang-hua@cnpc.com.cn

聯(lián)系作者簡介：王大興（1963-），男，甘肅武威人，博士，中國石油長慶油田公司教授級高級工程師，主要從事石油地球物理方法研究和地震儲集層預測工作。地址：陜西省西安市未央?yún)^(qū)長慶興隆園小區(qū)，長慶油田公司勘探開發(fā)研究院，郵政編碼：710018。E-mail：wdx1_cq@petrochina.com.cn

2016-04-08

2017-06 -09