黃土塬區(qū)致密儲(chǔ)集層模型地震正反演模擬

王大興，張盟勃，楊文敬，蔡克漢，高利東，朱軍（. 中國(guó)石油長(zhǎng)慶油田公司勘探開(kāi)發(fā)研究院，西安 7008；2. 低滲透油氣田勘探開(kāi)發(fā)國(guó)家工程實(shí)驗(yàn)室，西安 7008；3. 中國(guó)石油長(zhǎng)慶油田公司勘探部，西安 7008）

黃土塬區(qū)致密儲(chǔ)集層模型地震正反演模擬

王大興1, 2，張盟勃1, 2，楊文敬2, 3，蔡克漢1，高利東1，朱軍1
（1. 中國(guó)石油長(zhǎng)慶油田公司勘探開(kāi)發(fā)研究院，西安 710018；2. 低滲透油氣田勘探開(kāi)發(fā)國(guó)家工程實(shí)驗(yàn)室，西安 710018；3. 中國(guó)石油長(zhǎng)慶油田公司勘探部，西安 710018）

為明確鄂爾多斯盆地黃土塬區(qū)地表地震波場(chǎng)傳播規(guī)律和致密儲(chǔ)集層地震響應(yīng)特征，通過(guò)建立真實(shí)地表?xiàng)l件下儲(chǔ)集層的地質(zhì)-地球物理模型，開(kāi)展了全彈性地震波動(dòng)方程正演模擬和疊前地震彈性參數(shù)反演研究。模擬地震數(shù)據(jù)和實(shí)際地震數(shù)據(jù)同步處理分析表明：起伏疏松地表是導(dǎo)致地震靜校正問(wèn)題和干擾波發(fā)育的主要原因，近道約束的層析靜校正方法可以解決黃土塬地震靜校正問(wèn)題并提高其精度，疊前反演橫波波阻抗識(shí)別砂體效果明顯，泊松比可以被用于識(shí)別含油儲(chǔ)集層。地震正反演模擬研究為致密儲(chǔ)集層地震預(yù)測(cè)提供了理論和實(shí)驗(yàn)依據(jù)，并在鄂爾多斯盆地黃土塬區(qū)致密油勘探開(kāi)發(fā)中見(jiàn)到良好的應(yīng)用效果。圖11表3參12

黃土塬地表；致密儲(chǔ)集層；波動(dòng)方程正演模擬；近道約束層析靜校正；疊前彈性參數(shù)反演

引用：王大興, 張盟勃, 楊文敬, 等. 黃土塬區(qū)致密儲(chǔ)集層模型地震正反演模擬[J]. 石油勘探與開(kāi)發(fā), 2017, 44(2): 243-251.

WANG Daxing, ZHANG Mengbo, YANG Wenjing, et al. Seismic forward and inverse simulation in a tight reservoir model of loess plateau region[J]. Petroleum Exploration and Development, 2017, 44(2): 243-251.

0 引言

鄂爾多斯盆地中生界三疊系延長(zhǎng)組長(zhǎng) 7段致密油資源豐富，平面上分布于湖盆沉積中部，縱向上位于延長(zhǎng)組中部[1-2]，受地表?xiàng)l件和儲(chǔ)集層特征影響，致密油地震勘探面臨極大挑戰(zhàn)。首先，致密油目標(biāo)區(qū)地表為溝壑縱橫、起伏劇烈的黃土山地，導(dǎo)致地震干擾波非常發(fā)育、靜校正問(wèn)題極其嚴(yán)重，地震成像難度大；其次，長(zhǎng) 7段致密儲(chǔ)集層物性差、非均質(zhì)性強(qiáng)，與圍巖地震特性差異小，“甜點(diǎn)”識(shí)別困難。因此，弄清復(fù)雜地表地震波的傳播規(guī)律和致密儲(chǔ)集層的地球物理參數(shù)特征是黃土塬地震成像處理和致密儲(chǔ)集層預(yù)測(cè)的關(guān)鍵，開(kāi)展地質(zhì)-地球物理模型正反演模擬是揭示地震波傳播規(guī)律和驗(yàn)證地球物理方法有效性的重要手段。目前，波動(dòng)方程法是地震波場(chǎng)正演模擬的重要方法，國(guó)內(nèi)外很多學(xué)者對(duì)正演模擬的理論方法、波場(chǎng)的傳播規(guī)律做了大量研究[3]，但多為基于理論模型的算法研究，對(duì)于實(shí)際復(fù)雜地表模型的地震波傳播規(guī)律研究較少。在應(yīng)用研究方面，王雪秋、張華等將重點(diǎn)放在了復(fù)雜地表地質(zhì)條件的地震波數(shù)值模擬上，討論了起伏地表所面臨的諸多問(wèn)題[4-5]；裴正林等進(jìn)行了復(fù)雜地表復(fù)雜構(gòu)造模型的彈性波方程正演模擬研究，明確了起伏地形可造成不同炮點(diǎn)的波場(chǎng)特征差異及靜校正問(wèn)題[6]；王錫文、秦廣勝等采用優(yōu)化的波動(dòng)方程有限差分法正演模擬了復(fù)雜地區(qū)地震采集觀測(cè)系統(tǒng)，為地震采集設(shè)計(jì)提供了科學(xué)依據(jù)[7]；蘭鋒等應(yīng)用黏彈性波動(dòng)方程模擬地震低頻陰影，驗(yàn)證了低 Q（儲(chǔ)集層品質(zhì)因子）值儲(chǔ)集層的地震波衰減是產(chǎn)生低頻陰影現(xiàn)象的重要原因，為油水預(yù)測(cè)提供了理論依據(jù)[8]；王大興、張盟勃等通過(guò)對(duì)低滲透致密砂巖楔狀模型進(jìn)行波動(dòng)方程正演及疊前彈性參數(shù)反演研究，優(yōu)選了表征致密儲(chǔ)集層敏感性的彈性參數(shù)，明確了地震彈性參數(shù)對(duì)儲(chǔ)集層的分辨能力[9]。

本研究針對(duì)黃土塬區(qū)實(shí)際地表結(jié)構(gòu)、致密儲(chǔ)集層特征，建立了全剖面地質(zhì)-地球物理模型，基于完全彈性波動(dòng)方程理論，采用高階有限差分算法進(jìn)行模型正演模擬；通過(guò)野外地震單炮和正演模擬單炮記錄對(duì)比分析，明確了復(fù)雜地表結(jié)構(gòu)下地震波場(chǎng)特征和傳播規(guī)律；通過(guò)實(shí)際、模擬單炮同步處理和反演，對(duì)處理疊加剖面和反演結(jié)果進(jìn)行對(duì)比分析，驗(yàn)證了本文黃土塬地區(qū)關(guān)鍵處理技術(shù)和儲(chǔ)集層預(yù)測(cè)方法的可行性。

1 黃土塬區(qū)致密儲(chǔ)集層地質(zhì)-地球物理模型的建立

通過(guò)地表結(jié)構(gòu)、巖性及速度等參數(shù)調(diào)查測(cè)量，建立起伏地表結(jié)構(gòu)模型；分析致密儲(chǔ)集層地球物理參數(shù)特征，建立地下儲(chǔ)集層模型，進(jìn)而建立黃土塬地表?xiàng)l件下致密儲(chǔ)集層全剖面地質(zhì)-地球物理參數(shù)模型。建模的原則是：①參考的實(shí)際地震測(cè)線所在位置能反映典型的黃土塬地貌特征，所過(guò)井能代表致密儲(chǔ)集層的普遍特征；②地表結(jié)構(gòu)（高程、速度及密度）和儲(chǔ)集層參數(shù)（厚度、速度、密度）完全依據(jù)實(shí)際測(cè)量和鉆井測(cè)試數(shù)據(jù)；③地層格架嚴(yán)格采用鉆井地質(zhì)分層資料，地下儲(chǔ)集層模型參數(shù)在地質(zhì)分層的約束下從井點(diǎn)處內(nèi)插外推。

1.1 地表結(jié)構(gòu)參數(shù)模型

地表模型的建立嚴(yán)格參照鄂爾多斯盆地南部新安邊地區(qū)HCX測(cè)線近地表調(diào)查測(cè)量數(shù)據(jù)，圖1為地表結(jié)構(gòu)參數(shù)模型，可以看出，黃土塬地表高差變化大，高程范圍-1 800～-1 400 m，低降速帶厚度10～300 m，地表可分為 4層結(jié)構(gòu)，從上至下分別為風(fēng)化層、低速層、降速層及高速層，巖性分別為干黃土、黃土、膠泥和砂巖。真實(shí)黃土塬表層速度連續(xù)變化，為方便建模，各層縱波速度分別取值600 m/s、800 m/s、1 500 m/s 及2 500 m/s。

圖1 黃土塬地區(qū)地表結(jié)構(gòu)參數(shù)模型

1.2 致密油區(qū)地層、儲(chǔ)集層參數(shù)模型

本文所建致密油區(qū)地層和儲(chǔ)集層參數(shù)模型是在研究區(qū)新生界、白堊系、侏羅系及三疊系大套地質(zhì)層序的約束下，由經(jīng)過(guò)HCX測(cè)線的7口井（H223、H231、H232、X71、S544、S581、S582井）測(cè)井參數(shù)內(nèi)插外推得到。圖2為結(jié)合地表結(jié)構(gòu)模型建立的致密油區(qū)地層、儲(chǔ)集層地質(zhì)-地球物理模型，模型長(zhǎng)度為60 000 m，深度超過(guò)2 500 m，主要目的層為長(zhǎng)7段，模型網(wǎng)格大小為10 m（縱向）×2 m（橫向），地層和儲(chǔ)集層具體參數(shù)見(jiàn)表1、表2。

圖2 黃土塬區(qū)致密油區(qū)地層地質(zhì)-地球物理模型

表1 致密油區(qū)主要地層地球物理參數(shù)表

表2 致密油區(qū)主要巖性地球物理參數(shù)表

2 復(fù)雜介質(zhì)波動(dòng)方程數(shù)值模擬

2.1 二維各向同性介質(zhì)彈性模擬方程

在各向同性彈性波數(shù)值模擬時(shí)，對(duì)以位移表達(dá)的波動(dòng)方程進(jìn)行降階處理，得到速度-應(yīng)力表示的基本方程：

2.2 復(fù)雜介質(zhì)彈性波方程交錯(cuò)網(wǎng)格高階有限差分解法

傳統(tǒng)的有限差分方法在求解波動(dòng)方程時(shí)，會(huì)產(chǎn)生數(shù)值頻散（或稱網(wǎng)格頻散），導(dǎo)致模擬結(jié)果分辨率降低。為了提高有限差分的精度，減少數(shù)值模擬中頻散的發(fā)生，需要保留高階項(xiàng)，但這將會(huì)大大增加計(jì)算量。因此，針對(duì)黃土塬區(qū)起伏疏松的地表?xiàng)l件，采用交錯(cuò)網(wǎng)格高階有限差分解法。

2.2.1 有限差分法的精度分析

圖3中藍(lán)色曲線自下而上分別為K取1～10，即差分精度為 2～20階的曲線，紅色直線對(duì)應(yīng)于虛譜差分算子。由圖 3可以看出，一階導(dǎo)數(shù)的不同長(zhǎng)度中心差分算子在1/l值為0.5時(shí)，振幅譜為零，即中心差分算子精度不能達(dá)到虛譜差分算子的精度；而交錯(cuò)網(wǎng)格差分算子的振幅譜在接近1/l值為0.5時(shí)，略偏離于虛譜差分算子的譜線，但隨著差分算子長(zhǎng)度的增加，這種偏差越來(lái)越小?？梢?jiàn)，交錯(cuò)網(wǎng)格差分算子的精度可以達(dá)到虛譜差分算子的精度。取二階近似時(shí)，規(guī)則網(wǎng)格和交錯(cuò)網(wǎng)格上的近似誤差系數(shù)比為 0.25（(1/24)/(1/6)），四階近似時(shí)近似誤差系數(shù)比為0.141（(3/640)/(1/30)）；p（p為偶數(shù)）階近似時(shí)有，顯然，隨著階數(shù)的增加，交錯(cuò)網(wǎng)格的局部精度迅速提高。2.2.2 規(guī)則網(wǎng)格和交錯(cuò)網(wǎng)格上一階導(dǎo)數(shù)有限差分系數(shù)

的收斂性分析

從圖 4可以看出，隨著階數(shù)的不斷增高，規(guī)則網(wǎng)格上的差分系數(shù)（或權(quán)重）收斂于2，…）；而在交錯(cuò)網(wǎng)格上收斂于（m=1，2，…)，即規(guī)則網(wǎng)格上的差分系數(shù)收斂速度近似為1/m，而在交錯(cuò)網(wǎng)格上收斂速度近似為（1/m）2?？梢?jiàn)，交錯(cuò)網(wǎng)格上有限差分的權(quán)重衰減很快。因此，交錯(cuò)網(wǎng)格法的差分算子局部特性更好，收斂性好。

圖3 中心差分算子和交錯(cuò)網(wǎng)格差分算子分別與虛譜差分算子振幅譜對(duì)比圖

圖4 一階導(dǎo)數(shù)的20點(diǎn)交錯(cuò)差分算子與相應(yīng)虛譜差分算子對(duì)比

2.3 地質(zhì)-地球物理模型數(shù)值模擬

對(duì)前面建立的地質(zhì)-地球物理模型采用交錯(cuò)網(wǎng)格四階有限差分法進(jìn)行了波動(dòng)方程模擬，觀測(cè)系統(tǒng)完全按照野外實(shí)際地震采集參數(shù)設(shè)計(jì)。觀測(cè)系統(tǒng)具體參數(shù)如下：道檢距25 m，炮檢距50 m，最大炮檢距3 800 m，接收道數(shù)305道，覆蓋次數(shù)76次，點(diǎn)震源中間激發(fā)兩邊接收，35 Hz雷克子波激發(fā)，采樣間隔1 ms，所得單炮正演模擬結(jié)果及與實(shí)際地震記錄對(duì)比見(jiàn)圖5。

由于模擬中沒(méi)有考慮表層吸收衰減，也沒(méi)有加隨機(jī)噪音，所以正演模擬單炮（見(jiàn)圖 5a）波場(chǎng)特征較野外實(shí)際單炮（見(jiàn)圖5b）清晰，但基本能反映黃土塬地表地震單炮的共性：受地表高程和低降速帶的影響，黃土塬區(qū)單炮初至波存在不同程度的扭曲，在黃土層較厚地區(qū)，初至波有分層現(xiàn)象，直達(dá)波速度約為 800m/s，初至波速度約為1 500～2 500 m/s，靜校正問(wèn)題特別嚴(yán)重。另外單炮上面波比較發(fā)育，其他噪音有線性干擾、多次折射干擾等，從正演模擬單炮速度分析和建立的近地表速度模型對(duì)比認(rèn)為，單炮上的“帽子”是沿著淺層（干黃土、濕黃土層）水平傳播到檢波點(diǎn)的直達(dá)波，折射波分別是膠泥頂部和高速砂巖頂部所產(chǎn)生；面波是膠泥頂部的滑行波，主頻為8～12 Hz，速度700～1 200 m/s，存在頻散效應(yīng)；近炮點(diǎn)面波和有效波相消干涉，致使在面波形成的三角區(qū)域內(nèi)資料信噪比極低，資料的信噪比整體較低。

圖5 正演模擬單炮與野外實(shí)際單炮地震記錄

3 地震處理、解釋和儲(chǔ)集層預(yù)測(cè)方法驗(yàn)證分析

3.1 黃土塬區(qū)地震處理

上文分析認(rèn)為地表靜校正和去噪是黃土塬地震成像需要解決的關(guān)鍵問(wèn)題[10]。其中靜校正的難點(diǎn)在于地形起伏較大，黃土塬地表速度分層不明顯，折射界面不穩(wěn)定，初至波存在串層現(xiàn)象。理論上層析反演可以模擬任意復(fù)雜的近地表速度分布，但在沒(méi)有小折射、微測(cè)井約束的情況下，對(duì)于近地表淺層的刻畫(huà)比較粗糙。因此靜校正方法的選擇是地震成像的關(guān)鍵[11]。研究中認(rèn)為初至波的近道信息可以近似等價(jià)于小折射資料，將近道的信息充分利用起來(lái)，相當(dāng)于每一炮都做了小折射，將這些等價(jià)的小折射波信息進(jìn)行層析反演，就可以得到淺層高精度的速度模型，將反演的淺層模型作為約束權(quán)重場(chǎng)再進(jìn)行全初至波的層析反演，來(lái)提高整個(gè)地表模型的反演精度，因此提出了近道約束層析靜校正方法，其方法原理如下。

在層析反演中，M條射線N個(gè)未知數(shù)建立的層析方程組可表示為：

L個(gè)約束條件建立的約束方程組可表示為：

旅行時(shí)和約束方程構(gòu)建成的聯(lián)合方程組可表示為：

拉格朗日最優(yōu)約束的目標(biāo)函數(shù)可表示為：

圖6 正演模擬單炮層析反演的地表模型

圖7 正演單炮靜校正效果與實(shí)際單炮靜校正效果對(duì)比

圖8 正演模擬剖面（a）和實(shí)際地震剖面（b）對(duì)比

3.2 致密儲(chǔ)集層預(yù)測(cè)

如前所述，長(zhǎng) 7段致密儲(chǔ)集層非常致密，導(dǎo)致儲(chǔ)集層和非儲(chǔ)集層地球物理差異很小[12]，從測(cè)井彈性參數(shù)交會(huì)圖上看（見(jiàn)圖9），在縱波阻抗上，泥巖和砂巖儲(chǔ)集層相互疊置，含油砂巖、含水砂巖及干砂巖在波阻抗上完全重疊，常規(guī)縱波預(yù)測(cè)難度較大。橫波則可以很好地區(qū)分巖性，砂巖橫波阻抗大于 0.6×1010g/(m2·s)；泊松比基本可以區(qū)分含油砂巖和干砂巖，干砂巖的泊松比最大，含油砂巖泊松比最?。ㄒ话阈∮?0.23），含水砂巖居于兩者中間，這些地球物理參數(shù)特征是致密儲(chǔ)集層預(yù)測(cè)的基礎(chǔ)。

基于以上巖石物理分析結(jié)果，對(duì)正演模擬的數(shù)據(jù)和實(shí)際地震數(shù)據(jù)進(jìn)行了疊前彈性反演研究，分別求取表征巖性、含油性參數(shù)剖面進(jìn)行對(duì)比，驗(yàn)證兩個(gè)數(shù)據(jù)體反演結(jié)果的一致性，同時(shí)驗(yàn)證利用疊前反演資料進(jìn)行致密油“甜點(diǎn)”預(yù)測(cè)的可行性。圖10和圖11分別是地震正演模擬和實(shí)際地震數(shù)據(jù)的疊前反演剖面，可以看出，在同樣的數(shù)值（色標(biāo)）下，兩組數(shù)據(jù)的反演結(jié)果反映儲(chǔ)集層砂巖、油層的變化趨勢(shì)基本一致，自西向東，砂體逐漸尖滅，含油性橫向非均質(zhì)性強(qiáng)。對(duì)比反演參數(shù)剖面和試油結(jié)果（見(jiàn)表3），橫波阻抗越大，泊松比越小，代表砂巖厚度越大，油層越厚，產(chǎn)量越高，與完鉆結(jié)果基本一致，證實(shí)了利用疊前反演進(jìn)行致密油“甜點(diǎn)”預(yù)測(cè)是可行的。

圖9 縱、橫波波阻抗與泊松比交會(huì)圖

近年來(lái)，針對(duì)鄂爾多斯盆地黃土塬地表，應(yīng)用近道約束層析靜校正、多域多方法噪音壓制為核心的關(guān)鍵地震處理技術(shù)，最終處理剖面上中生界目的層主頻達(dá)40 Hz，古生界目的層地震主頻達(dá)35 Hz，井震吻合程度為 85%，保真度較高，提高了黃土塬地震成像質(zhì)量。在致密儲(chǔ)集層預(yù)測(cè)中，以疊前彈性參數(shù)反演為主，輔以疊后多屬性分析技術(shù)，提高了致密儲(chǔ)集層預(yù)測(cè)的精度，在實(shí)際致密油勘探開(kāi)發(fā)時(shí)，使致密儲(chǔ)集層預(yù)測(cè)符合率達(dá)77.3%，含油性預(yù)測(cè)符合率達(dá)71.4%。

圖10 模擬數(shù)據(jù)疊前反演的橫波阻抗（a）和泊松比（b）剖面

圖11 實(shí)際數(shù)據(jù)疊前反演的橫波阻抗（a）和泊松比（b）剖面

表3 鉆井儲(chǔ)集層及試油參數(shù)

4 結(jié)論

基于真實(shí)地表地質(zhì)-地球物理模型的波動(dòng)方程正演模擬是了解地震波場(chǎng)傳播規(guī)律的直觀方法，全彈性地震波動(dòng)方程的四階有效差分法能解決模擬中的頻散和邊界問(wèn)題，具有收斂快、精度高的特點(diǎn)；模擬的單炮記錄和實(shí)際單炮記錄具有很好的可對(duì)比性，波場(chǎng)分析表明，起伏疏松的黃土塬地表結(jié)構(gòu)使得單炮記錄的干擾波非常發(fā)育，復(fù)雜地表和低降速帶的變化是導(dǎo)致靜校正問(wèn)題的主要原因。

模擬單炮和實(shí)際單炮同步處理解釋是驗(yàn)證地震處理和儲(chǔ)集層預(yù)測(cè)方法有效性的重要手段，正反演數(shù)據(jù)處理表明，利用大炮檢距近道初至波約束的層析靜校正方法反演的表層速度模型精度較高，能夠解決復(fù)雜地表?xiàng)l件的靜校正問(wèn)題；多域多方法聯(lián)合去噪方法提高了資料的信噪比，疊前同時(shí)反演的橫波阻抗、泊松比參數(shù)可以表征致密儲(chǔ)集層巖性及含油性，并取得了實(shí)際應(yīng)用效果。

符號(hào)注釋：

ak——網(wǎng)格差分系數(shù)，無(wú)因次；A——Jacob矩陣，無(wú)因次；A1，A2——彈性波方程系數(shù)矩陣，無(wú)因次；An——?dú)w一化振幅，無(wú)因次；B——（M+K）×N維矩陣，無(wú)因次；cp——p階近似時(shí)近似誤差系數(shù)比，無(wú)因次；C——由約束條件形成的K×N維矩陣，無(wú)因次；Co——算子系數(shù)，無(wú)因次；fi——第i個(gè)網(wǎng)格單元約束和無(wú)約束反演結(jié)果殘差，m/s；F——約束條件下的旅行時(shí)差矩陣，s；h——空間采樣間隔，m；i——層析反演中網(wǎng)格單元序號(hào)；H——所有約束條件下的旅行時(shí)殘差矩陣，s；j——層析反演中用于約束的網(wǎng)格單元序號(hào)；k——算子長(zhǎng)度變量，取值為1—K，無(wú)因次；K——算子長(zhǎng)度的一半，無(wú)因次；l——網(wǎng)格點(diǎn)數(shù)；m——階數(shù)；M——射線條數(shù)，即道數(shù)；N——模型網(wǎng)格數(shù)；p——偶數(shù)差分階數(shù)，無(wú)因次； is~——第i個(gè)網(wǎng)格約束的反演速度，m/s；isΔ ——第i個(gè)網(wǎng)格的慢度修正量；ΔS——慢度修正量，無(wú)因次；t——時(shí)間，s；Δti——第i個(gè)網(wǎng)格的地震道實(shí)際初至?xí)r間與模型正演初至?xí)r間之差，s；ΔT——實(shí)際和正演走時(shí)時(shí)差矩陣，s；vx，vz——沿x，z方向的速度，m/s；β——約束系數(shù)，無(wú)因次；λ——介質(zhì)的拉梅常數(shù)，Pa；μ——介質(zhì)的剪切模量，Pa；ρ——密度，g/cm3；τxx，τzz——沿x，z方向的正應(yīng)力，Pa；τzx——切向應(yīng)力，Pa；ω——角頻率，rad/s。

[1] 楊華, 付金華. 超低滲透油藏勘探理論與技術(shù)[M]. 北京: 石油工業(yè)出版社, 2012. YANG Hua, FU Jinhua. Exploration theory and technology of ultra low permeability reservoir[M]. Beijing: Petroleum Industry Press, 2012.

[2] 姚涇利, 鄧秀琴, 趙彥德, 等. 鄂爾多斯盆地延長(zhǎng)組致密油特征[J]. 石油勘探與開(kāi)發(fā), 2013, 40(2): 150-158. YAO Jingli, DENG Xiuqin, ZHAO Yande, et al. Characteristics of tight oil in Triassic Yanchang Formation, Ordos Basin[J]. Petroleum Exploration and Development, 2013, 40(2): 150-158.

[3] 牟永光, 裴正林. 三維復(fù)雜介質(zhì)地震數(shù)值模擬[M]. 北京: 石油工業(yè)出版社, 2005. MU Yongguang, PEI Zhenglin. Three dimensional seismic numerical modeling complex medium[M]. Beijing: Petroleum Industry Press, 2005.

[4] 王雪秋, 孫建國(guó), 張文志. 復(fù)雜地表地質(zhì)條件下地震波數(shù)值模擬綜述[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2005, 35(增刊): 12-18. WANG Xueqiu, SUN Jianguo, ZHANG Wenzhi. The state of the art in numerical modeling including surface topograph[J]. Journal of Jilin University (Earth Science Edition), 2005, 35(Supp.): 12-18.

[5] 張華, 李振春, 韓文功. 起伏地表?xiàng)l件下地震波數(shù)值模擬方法綜述[J]. 勘探地球物理進(jìn)展, 2007, 30(5): 334-339. ZHANG Hua, LI Zhenchun, HAN Wengong. Summary of seismic wave numerical simulation method under the condition of undulating terrain[J]. Progress in Exploration Geophysics, 2007, 30(5): 334-339.

[6] 裴正林, 何光明, 謝芳. 復(fù)雜地表復(fù)雜構(gòu)造模型的彈性波方程正演模擬[J]. 石油地球物理勘探, 2010, 45(6): 807-818. PEI Zhenglin, HE Guangming, XIE Fang. Elastic wave equation forward modeling for complex surface and complex structure model[J]. Petroleum Geophysical Exploration, 2010, 45(6): 807-818.

[7] 王錫文, 秦廣勝, 趙衛(wèi)鋒, 等. 正演模擬技術(shù)在地震采集設(shè)計(jì)中的應(yīng)用[J]. 地球物理學(xué)進(jìn)展, 2012, 27(2): 642-650. WANG Xiweng, QIN Guangsheng, ZHAO Weifeng, et al. The application of forward modeling technique in seismic acquisition design[J]. Progress in Geophysics, 2012, 27(2): 642-650.

[8] 蘭鋒, 李明, 李艷東, 等. 應(yīng)用黏彈性波動(dòng)方程模擬地震低頻陰影[J]. 石油地球物理勘探, 2015, 50(1): 84-90. LAN Feng, LI Ming, LI Yandong, et al. Low-frequency shadow numerical simulation with visco-elastic wave equation[J]. Petroleum Geophysical Exploration, 2015, 50(1): 84-90.

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

[10] 王大興, 張杰, 高利東. 非縱地震勘探技術(shù)及其在鄂爾多斯盆地的應(yīng)用[J]. 石油地球物理勘探(物探技術(shù)研討會(huì)?？?, 2012: 133-136. WANG Daxing, ZHANG Jie, GAO Lidong. Off-line seismic exploration technology and application in the Ordos Basin[J]. Petroleum Geophysical Exploration(Special Issue), 2012: 133-136.

[11] 陳娟, 胡建, 王永剛, 等. 鄂爾多斯盆地黃土塬區(qū)地震資料處理方法及應(yīng)用[J]. 石油天然氣學(xué)報(bào), 2012, 34(5): 65-68. CHEN Juan, HU Jian, WANG Yonggang, et al. Seismic data processing method and its application in loess plateau Ordos Basin[J]. Journal of Petroleum and Natural Gas, 2012, 34(5): 65-68.

[12] 王大興, 王永剛. 致密儲(chǔ)層巖石力學(xué)參數(shù)預(yù)測(cè)方法研究及應(yīng)用[R].珀斯: 第三屆國(guó)際巖石物理研討會(huì), 2015. WANG Daxing, WANG Yonggang. Study and application of rock mechanics parameter prediction method for tight reservoir[R]. Perth: 3rd International Workshop on Rock Physics, 2015.

（編輯 黃昌武）

Seismic forward and inverse simulation in a tight reservoir model of loess plateau region

WANG Daxing1, 2, ZHANG Mengbo1, 2, YANG Wenjing2, 3, CAI Kehan1, GAO Lidong1, ZHU Jun1
(1. Exploration and Development Research Institute of PetroChina Changqing Oilfield Company, Xi’an 710018, China; 2. National Engineering Laboratory for Exploration and Development of Low-Permeability Oil & Gas Fields, Xi’an 710018, China; 3. Exploration Department of PetroChina Changqing Oilfield Company, Xi’an 710018, China)

To find out the seismic wave field propagation principles in loess plateau near surface of the Ordos Basin and the seismic response characteristics of tight oil reservoirs, this study established a geological-geophysical model under the real conditions of ground surface of loess plateau, and launched full elastic seismic wave equation forward modeling and pre-stack elastic seismic inversion study. Comparison of modeling and real seismic data shows that, the loose and wavy loess plateau surface is the main reason for causing the problems of seismic static correction and interference wave. Tomographic static correction method with the constraint of traces near shot point can effectively solve the problem of seismic static correction in the loess plateau and enhance seismic interpretation accuracy, S-wave impedance obtained from pre-stack seismic inversion can identify sandstone effectively, and Poisson's ratio can identify oil-bearing reservoirs. The seismic forward and inverse simulation and rock physical analysis provide a solid theoretical and experimental basis for the seismic prediction of tight oil reservoir, and have worked well in the oil exploration and development in the loess plateau of the Ordos Basin.

loess plateau surface; tight oil reservoir; wave equation forward simulation; tomographic inversion static correction; pre-stack elastic inversion

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

TE<122.2 class="emphasis_bold">122.2 文獻(xiàn)標(biāo)識(shí)碼：A122.2

A

1000-0747(2017)02-0243-09

10.11698/PED.2017.02.08

王大興（1963-），男，甘肅武威人，博士，中國(guó)石油長(zhǎng)慶油田公司教授級(jí)高級(jí)工程師，主要從事石油地球物理方法研究和地震儲(chǔ)集層預(yù)測(cè)工作。地址：陜西省西安市未央?yún)^(qū)，中國(guó)石油長(zhǎng)慶油田公司勘探開(kāi)發(fā)研究院，郵政編碼：710018。E-mail：wdx1_cq@petrochina.com.cn

2016-05-27

2017-02-17