松遼盆地葡南油田西緣薄砂體儲層地震預(yù)測

尹兵祥， 邱長星， 劉洪濤

(1.中國石油大學(xué)(華東) 地球科學(xué)與技術(shù)學(xué)院,青島 266580；2.大慶油田有限公司 第七采油廠,大慶 163517)

?

松遼盆地葡南油田西緣薄砂體儲層地震預(yù)測

尹兵祥1， 邱長星1， 劉洪濤2

(1.中國石油大學(xué)(華東) 地球科學(xué)與技術(shù)學(xué)院,青島 266580；2.大慶油田有限公司 第七采油廠,大慶 163517)

葡南油田開發(fā)層屬姚一段葡萄花油層，儲層為三角洲前緣亞相的水下分流河道砂體、席狀砂沉積和湖相濱湖砂沉積，砂體儲層窄小，橫向變化快，連續(xù)性差，地震預(yù)測是必不可少的方法。通過建立精細(xì)地質(zhì)模型，進(jìn)行多井約束波阻抗反演，采用多屬性儲層預(yù)測技術(shù)，結(jié)合測井解釋結(jié)果，對葡南油田西部邊緣葡331-1區(qū)塊的儲層橫向分布進(jìn)行了預(yù)測，預(yù)測結(jié)果符合砂體分布特征，井點(diǎn)符合率高，對于認(rèn)識該區(qū)油氣分布規(guī)律和開發(fā)井位部署提供了直接的技術(shù)支持，對相鄰地區(qū)的勘探開發(fā)具有參考意義。

葡南油田； 葡萄花油層； 薄砂體儲層； 地震預(yù)測

0 引言

由于地震資料在橫向上具有較強(qiáng)的連續(xù)性，且地震資料中包含著豐富的儲層物性信息，相比之下，鉆井、測井資料則因?yàn)闄M向稀疏而在井間對比和預(yù)測中不具有優(yōu)勢，因此地震技術(shù)在儲層預(yù)測中發(fā)揮著重要的作用[1-6]。作為儲層預(yù)測的核心技術(shù)之一的井約束波阻抗地震反演基于構(gòu)造地質(zhì)模型，以井點(diǎn)阻抗為約束條件，從地震波數(shù)據(jù)求取地層波速或波阻抗數(shù)據(jù)，進(jìn)而根據(jù)層間或是層內(nèi)波速和/或波阻抗差異進(jìn)行不同巖性的區(qū)分[5-10]。將反演阻抗和波速甚至其他類型的數(shù)據(jù)體類比于地震屬性，將其與原始地震數(shù)據(jù)的屬性類型共同進(jìn)行儲層預(yù)測，可以進(jìn)一步提高儲層地震預(yù)測的效果，更準(zhǔn)確地描述儲層分布狀態(tài)和油藏儲量規(guī)模[10-12]。

1 研究區(qū)概況

葡南油田在構(gòu)造上位于松遼盆地大慶長垣南部的葡萄花構(gòu)造北部，主要開發(fā)目的層為下白堊統(tǒng)姚家組一段葡一組油層組，該層組為一套典型的砂泥巖薄互層，厚度約45 m～55 m，自上而下分為11個(gè)小層，儲層為三角洲前緣亞相的水下分流河道砂體、席狀砂沉積和湖相濱湖砂沉積，單層最大厚度在1 m～2 m之間。圖1為葡154-28井的單井綜合圖，葡I油層組厚度為55.4 m，其中砂巖厚度為11.7 m，共7層，單層最大2.8 m。砂體橫向變化快，規(guī)律性差，整個(gè)層組則表現(xiàn)為疊合連片的狀態(tài)[13]。葡331-1區(qū)塊位于葡南油田西緣，已于2011年提交控制儲量，面積為7.3 km2，控制儲量58×104t，表現(xiàn)出一定的資源潛力。構(gòu)造位置相對較低的葡148-16井PI1、PI3、PI4號層為水層，PI5-PI10號層為高含水低產(chǎn)工業(yè)油層，呈現(xiàn)明顯的上水下油情況，表明該區(qū)油水關(guān)系較為復(fù)雜。同時(shí)該區(qū)儲層橫向變化快，連續(xù)性差，單純依靠井間對比難以滿足儲層預(yù)測的需要，必須借助于地震數(shù)據(jù)開展儲層預(yù)測工作[14-16]。

圖2 葡331-1區(qū)塊東西向典型地震剖面

圖1 葡158-28單井綜合圖

圖2為研究區(qū)內(nèi)一條東西向地震剖面及目的層構(gòu)造解釋，顯示構(gòu)造趨勢為東高西低，發(fā)育有若干正斷層，最大落差超過一個(gè)地震子波周期，不同斷層的組合將目的層切割成地壘、地塹、斷階等斷塊，局部也具有微幅度構(gòu)造的顯示，為油氣聚集提供了不同的場所和條件。圖3為研究區(qū)葡I油層組頂面構(gòu)造圖，顯示總體構(gòu)造形態(tài)為東南高、西北低，被若干條斷層切割。最顯著的構(gòu)造高點(diǎn)在研究區(qū)東南角，界面海拔深度小于900 m；在東部南北向大斷層以西的斷層上升盤有兩個(gè)斷塊高點(diǎn)，具有形成油氣藏的構(gòu)造條件，是現(xiàn)有井網(wǎng)外擴(kuò)的有利位置。

2 井約束波阻抗反演

2.1 波阻抗反演的基本原理

井約束地震反演采用地震資料構(gòu)造解釋成果建立地質(zhì)模型，通過合成地震記錄實(shí)現(xiàn)層位的準(zhǔn)確標(biāo)定，根據(jù)相鄰地震道之間的相關(guān)性質(zhì)，利用神經(jīng)網(wǎng)絡(luò)方法自適應(yīng)地外推地震子波[3]，井旁道以測井資料為約束條件，實(shí)現(xiàn)阻抗等數(shù)據(jù)的計(jì)算，并在直接反演的中頻相對阻抗數(shù)據(jù)的基礎(chǔ)上，補(bǔ)充測井?dāng)?shù)據(jù)所包含的低頻和高頻信息，拓展阻抗數(shù)據(jù)體的頻寬，以便提高分辨率并得出巖層絕對波阻抗。由于砂體儲層一般比圍/泥巖阻抗值大，因此有利于巖性識別，地震反演已經(jīng)成為儲層預(yù)測的重要技術(shù)[8-9]。圖4為測井約束地震反演的具體流程。

圖3 葡331-1區(qū)塊葡I油層組頂面構(gòu)造圖

圖4 測井約束地震反演方法流程

波阻抗反演的原理為：給定多層地質(zhì)模型中各層的波速、密度、厚度等參數(shù)分別為：

v(i)，ρ(i)，d(i)，i=1,2,3,…,N

各層中地震波垂直傳播的時(shí)間為τ(i)=2d(i)/v(i)，目的層界面的反射時(shí)間為其上各層內(nèi)時(shí)間的累加即

(1)

由此可以將地震記錄描述為

i=1,2,3,…,M

(2)

其中：s(i)為實(shí)際地震信號；r為地層界面的反射系數(shù)；w為地震子波；M為樣點(diǎn)總個(gè)數(shù)；i為地震記錄內(nèi)的樣點(diǎn)序號。式(2)的矩陣表示形式為式(3)。

S=W·R

(3)

式中：S=(s(1)s(2)…s(m))T；

R=(r(1)r(2)…r(n))T；

W為M×N階子波矩陣：

(4)

地層反射系數(shù)與其對數(shù)波阻抗的關(guān)系用矩陣表示為式(5)。

R=D·L

(5)

其中：R=[r(1)r(2) …r(N)]T；

L=[l(1)l(2) …l(N)]T；

D為N行，N+1列的系數(shù)矩陣；l(i)為地層模型中各層波阻抗的對數(shù)值，將式(5)帶入式(3)得式(6)。

S=WDL

(6)

設(shè)實(shí)際測量的地震記錄為：T=[t(1)t(2)…t(M)]T，模型道S與實(shí)際地震道T之間的差為E=T-S，而誤差能量可以表示為式(7)。

J=ETE=(T-WDL)T(T-WDL)

(7)

這樣目標(biāo)函數(shù)J使待求波阻抗與實(shí)際地震觀測資料發(fā)生直接聯(lián)系，使J達(dá)到最小的物理含義是：尋求一個(gè)最佳的地質(zhì)模型，使由此模型計(jì)算的合成地震數(shù)據(jù)與實(shí)際觀測資料的誤差能量最小。

式(7)的最小平方解為式(8)。

L=(DTWTWD)-1DTWTT

(8)

由于式(8)無法引入約束條件，所以實(shí)際地震記錄中的噪音會嚴(yán)重地影響反演結(jié)果，可以考慮用共軛梯度法求解來避免此問題[5]，該方法的主要優(yōu)點(diǎn)是①算法精確、穩(wěn)定；②不需做矩陣反演計(jì)算，避免了大矩陣處理中的病態(tài)問題；③具有比較強(qiáng)的抗干擾能力；④求解過程中能容易地執(zhí)行約束條件。因此井約束波阻抗反演通常采用該方法，通過多次迭代和修改參數(shù)、逐次逼近地層模型來求取地層的波阻抗[6-7]。

2.2 地震子波估算

地震子波估算是反演的關(guān)鍵步驟之一，合理的子波是高精度反演的前提[9]。獲取子波的方式一般有標(biāo)準(zhǔn)理論子波、統(tǒng)計(jì)估算子波和地震提取子波等三種。從地震資料中提取的子波與地震數(shù)據(jù)匹配良好，可以很好地實(shí)現(xiàn)合成記錄與井旁道的相關(guān)性。本次反演在各井提取到的子波都非常接近于主頻為50 Hz的Ricker子波，因此直接采用了50 Hz的Ricker子波。

2.3 地層框架模型的建立

建立地層框架模型就是定義各地層界面及斷層面之間的相互關(guān)系，構(gòu)建地層空間格架，為將測井曲線數(shù)據(jù)通過內(nèi)插外推形成三維數(shù)據(jù)體提供約束和控制，使其橫向內(nèi)插外推沿層進(jìn)行，而不是以等深為原則，這樣得到的測井曲線類型的數(shù)據(jù)體才更合理。同時(shí)該框架模型也對地震數(shù)據(jù)的反演計(jì)算具有模型控制作用。本次解釋的目的層頂、層底反射的時(shí)間厚度約為30 ms～40 ms，只有兩個(gè)界面，為保證目的層反演的效果，將頂界面、底界面分別向上和向下推移10 ms，分別形成另外的兩個(gè)層用作模型的頂、底界面。

2.4 高低頻分量補(bǔ)充

相對于地震反演計(jì)算直接得到的阻抗數(shù)據(jù)，井點(diǎn)的測井曲線計(jì)算的阻抗數(shù)據(jù)具有顯著的分辨率優(yōu)勢，而且其垂向的緩慢變化趨勢也是相對清楚的，從各井的鉆井阻抗內(nèi)插外推得到的測井阻抗數(shù)據(jù)體，也同樣具有很寬的頻帶范圍和很高的垂向分辨率[11]。

地震反演阻抗數(shù)據(jù)的高低頻分量補(bǔ)充就是將其與測井阻抗數(shù)據(jù)體按道進(jìn)行合并，將地震數(shù)據(jù)缺失的低頻和高頻成分從測井阻抗數(shù)據(jù)體中提取合并到地震反演阻抗數(shù)據(jù)體中來[5,9]。根據(jù)原始地震資料及直接反演阻抗的頻譜范圍，確定低頻補(bǔ)償?shù)拈撝禐?0 Hz，高頻補(bǔ)償?shù)念l段確定為90 Hz ～220 Hz，高、低頻率補(bǔ)償后就得到能反映地層波阻抗的數(shù)據(jù)體。

2.5 三維地震反演效果分析

反演效果主要從以下幾個(gè)方面分析：①阻抗數(shù)據(jù)的縱橫向分布變化是否符合地質(zhì)規(guī)律；②井旁道與鉆井?dāng)?shù)據(jù)之間的吻合程度；③阻抗數(shù)據(jù)比原始地震數(shù)據(jù)的分辨率要有明顯的提高[9]。由于反演過程中的多個(gè)環(huán)節(jié)都加強(qiáng)了質(zhì)量控制，保證了反演阻抗數(shù)據(jù)具有較高的分辨率，能夠較好地反映目的層巖性變化特征。從圖5可以看出，目的層相對于上、下地層是一套高阻抗的地層，這與實(shí)際地質(zhì)情況是符合的。

3 砂巖儲層地震預(yù)測

3.1 多屬性儲層地震預(yù)測

儲層地震預(yù)測以地震資料為主，綜合地質(zhì)、測井、試井等多種資料分析研究油氣儲層的空間分布、巖性物性、流體性質(zhì)等油氣藏特征[10]。為此需要分析待預(yù)測的地質(zhì)目標(biāo)數(shù)據(jù)與各種地震屬性之間的相關(guān)關(guān)系，不僅從原始地震數(shù)據(jù)提取屬性，而且從反演成果數(shù)據(jù)提取屬性，幾種典型屬性如圖6所示，由于地震資料是多種地質(zhì)因素的綜合效應(yīng)，單個(gè)屬性難以有效反應(yīng)儲層分布。

圖7為多屬性儲層地震預(yù)測的技術(shù)思路，首先是將井點(diǎn)的儲層厚度數(shù)據(jù)與多種地震(包括反演阻抗)屬性進(jìn)行相關(guān)分析，確定建立預(yù)測模型所需的屬性類型，選取的原則首先是相關(guān)系數(shù)的大小，其次是選取的屬性應(yīng)該相互之間不相關(guān)，屬性的數(shù)量不宜太多，此次預(yù)測選擇了3種不同的屬性，根據(jù)所選擇的屬性類型，采用最小二乘法確定預(yù)測模型各屬性的系數(shù)，得到本次儲層預(yù)測的模型關(guān)系為：儲層厚度=1.82+0.316×10-6×阻抗+0.418×10-4×均方根振幅-0.264×能量半衰時(shí)，這樣就可以根據(jù)該模型做預(yù)測了，但是初次預(yù)測結(jié)果在井點(diǎn)的預(yù)測值與真實(shí)值難免存在誤差，將各井點(diǎn)的偏差分別求出，進(jìn)行網(wǎng)格化對應(yīng)到各個(gè)地震道點(diǎn)，然后將該偏差從初步預(yù)測數(shù)據(jù)中對應(yīng)相減，得到最后的預(yù)測結(jié)果。除上述根據(jù)屬性回歸模型的預(yù)測方法外，還可以采用地質(zhì)統(tǒng)計(jì)學(xué)和協(xié)克里金的方法進(jìn)行預(yù)測。由于地震數(shù)據(jù)的最小有效單元為一個(gè)同相軸，厚度更小的薄層預(yù)測可信度降低，因此這里最多只能做到層組的預(yù)測。

3.2 儲層厚度預(yù)測結(jié)果

采用前述儲層地震預(yù)測技術(shù)，我們對葡331-1區(qū)塊的砂巖厚度進(jìn)行了預(yù)測。圖8為整個(gè)葡I油層組的儲層厚度預(yù)測，按照0.3 m誤差為限，井點(diǎn)符合率大于80 %。從中可以看出，研究區(qū)葡I油層組儲層厚度在0 m～9 m之間，分布比較零散，符合該區(qū)單個(gè)砂體儲層面積小、厚度薄的特征，在多層砂體疊合的情況下，仍然難有明顯的分布規(guī)律性，而是呈現(xiàn)出眾多小的厚度中心。盡管砂體厚度不大，但是由于該區(qū)域油氣地質(zhì)條件好，油源充分，在構(gòu)造條件的配合下，砂體成藏幾率非常大，加上埋藏深度只有1 000 m 左右，即使比較小的含油砂體也是具有開發(fā)價(jià)值的。

圖5 反演阻抗數(shù)據(jù)體的剖面顯示

圖6 典型屬性示例

圖7 多屬性儲層預(yù)測技術(shù)示意圖

表1 葡331-1區(qū)塊潛力區(qū)有利目標(biāo)評價(jià)

3.3 有利目標(biāo)優(yōu)選

綜合構(gòu)造、儲層、油氣顯示等因素，在葡南油田西緣區(qū)塊開發(fā)井網(wǎng)以外的潛力區(qū)，篩選了3個(gè)局部的有利目標(biāo)區(qū)(如圖8中點(diǎn)線區(qū)域所示)，評價(jià)的相關(guān)參數(shù)及評價(jià)結(jié)果列于表1，同時(shí)給出了葡I油層及分層組預(yù)測的儲層厚度數(shù)據(jù)。

圖8 葡331-1區(qū)塊儲層地震預(yù)測

4 結(jié)論

針對砂巖儲層與泥巖在地震反射波屬性及波阻抗等方面的差異，通過三維地震資料和鉆井、測井資料在模型框架下的有機(jī)結(jié)合，利用井約束波阻抗反演和多屬性分析技術(shù)預(yù)測了目的層的薄互層砂巖的橫向分布，深化了對油藏特征和規(guī)律的認(rèn)識，有助于增儲潛力區(qū)塊優(yōu)選和油藏地質(zhì)評價(jià)。地震反演效果與軟件的方法技術(shù)和資料的品質(zhì)好壞有著極為重要的關(guān)系，也受井位數(shù)量及分布情況等因素的影響。實(shí)踐中應(yīng)根據(jù)研究區(qū)的地質(zhì)特征及具體資料情況，選擇合適的反演方法及流程和參數(shù)，在單一阻抗數(shù)據(jù)預(yù)測效果不佳的情況下，可以考慮綜合波速、密度等多種數(shù)據(jù)類型進(jìn)行綜合預(yù)測。

[1] 樂友喜, 楊麗. 儲層地震預(yù)測基礎(chǔ)理論方法研究[J].天然氣地球科學(xué),2009,20(4): 563-570. YUE Y X, YANG L. Basic theories and methods of seismic reservoir prediction. Natural gas geoscience, 2009,20(4): 563-570.(In Chinese)

[2] WANG Y H.Sparseness-constrained least-squares inversion: application to seismic wave reconstruction[J]. Geophysics, 2003,68(5):1633-1638.

[3] COOKE D J.Techniques for improving seismic attribute versus well-log property cross plots with synthetic and real data examples from Prudhoe Bay Field[C].Annual Meeting Abstracts, Society of Exploration Geophysicists, 1999:69-72.

[4] PEARSON R,HART B.Convergence of 3-D seismic attribute based reservoir property prediction and geologic interpretation as a risk reduction tool: A case study from a Permian intraslope basin[C]. Annual Meeting Abstracts, Society of Exploration Geophysicists,1999:896-899.

[5] 郭朝斌,楊小波, 陳紅岳，等. 約束稀疏脈沖反演在儲層預(yù)測中的應(yīng)用[J].石油物探,2006,45(4):397-400. GUO CH B, YANG X B, CHEN H Y, et al. Constrained sparse pulse research in north of Haitongji depression. Geophysical Prospecting for Petroleum, 2006,45(4):397-400.(In Chinese)

[6] 何胡軍, 王秋語, 程會明. 基于匹配追蹤算法子波分解技術(shù)在薄互層儲層預(yù)測中的應(yīng)用[J]. 物探化探計(jì)算技術(shù), 2010, 32(6):641-644. HE H J, WANG Q Y, CHENG H M. The application of wavelet decomposition technique based on matching pursuit algorithm in thin interbedded reservoir prediction. Computing Techniques for Geophysical and Geochemical Exploration, 2010, 32(6):641-644.(In Chinese)

[7] 劉丹, 徐 偉. 隨機(jī)反演在陸豐13-1油田儲層預(yù)測中的應(yīng)用[J]. 物探化探計(jì)算技術(shù), 2012, 34(3):331-335. LIU D, XU W. Application of stochastic inversion in the reservoir prediction of Lufeng 13-1 oil field. Computing Techniques for Geophysical and Geochemical Exploration, 2012, 34(3):331-335.(In Chinese)

[8] 韓翀, 臧殿光, 李建華.地質(zhì)統(tǒng)計(jì)學(xué)反演在四川盆地L地區(qū)儲層預(yù)測中的應(yīng)用[J]. 物探化探計(jì)算技術(shù), 2010, 32(3):310-315. HAN CH, ZANG D G, LI J H. The application of geostatistical inversion to reservoir prediction of l area in Sichuan basin. Computing Techniques for Geophysical and Geochemical Exploration, 2010, 32(3):310-315.(In Chinese)

[9] 隋淑玲，唐軍，蔣宇冰，等. 常用地震反演方法技術(shù)特點(diǎn)與適用條件[J]. 油氣地質(zhì)與采收率, 2012, 19(4):38-41. SUI SH L, TANG J, JIANG Y B, et al. Technical characteristics and applicable conditions analysis on seismic inversion traditional methods. Petroleum Geology and Recovery Efficiency, 2012, 19(4):38-41.(In Chinese)

[10]郭華軍，劉慶成．地震屬性技術(shù)的歷史、現(xiàn)狀及發(fā)展趨勢勢[J]. 物探與化探, 2008,32(1):19-22. GUO H J, LIU Q CH. The discussion of earthquake attribute technology's history, present situation and development tendency. Geophysical & Geochemical Exploration, 2008,32(1):19-22.(In Chinese)

[11]楊飛，趙政權(quán). 利用地震屬性進(jìn)行砂巖儲層橫向預(yù)測[J]. 物探化探計(jì)算技術(shù), 2012, 34(4):415-418. YANG F, ZHAO ZH Q. The application of seismic attributes analysis to sandstone reservoir lateral prediction. Computing Techniques for Geophysical and Geochemical Exploration, 2012, 34(4):415-418.(In Chinese)

[12]張國棟, 劉 萱, 田麗花,等. 綜合應(yīng)用地震屬性與地震反演進(jìn)行儲層描述[J]. 石油地球物理勘探, 2010, 45(增刊1):137-144. ZHANG G D, LIU X, TIAN L H, et al. Integrated application of seismic attribute and seismic inversion in reservoir characterization. Oil Geophysical Prospecting, 2010, 45(Supplement 1): 137-144.(In Chinese)

[13]趙長勛, 姜巖, 趙峰華,等. 葡南油田薄層河道砂巖體分布規(guī)律及預(yù)測[J]. 油氣地質(zhì)與采收率, 2009, 16(3):4-6+10. ZHAO CH X, JIANG Y, ZHAO F H, et al. Distribution and prediction of the thin underwater distributary channel sandbodies in Punan oilfield. Petroleum Geology and Recovery Efficiency, 2009, 16(3):4-6+10.(In Chinese)

[14]劉洪濤. 松遼盆地葡南油田砂體預(yù)測與評價(jià)[J].大慶石油地質(zhì)與開發(fā), 2010, 29(3):158-162. LIU H T. Prediction and evaluation of the sandbodies in Punan oilfield of Songliao basin. Petroleum Geology and Oilfield Development in Daqing, 2010, 29(3):158-162.(In Chinese)

[15]侯君, 王金榮, 王長生,等. 薄差儲層預(yù)測技術(shù)在葡南油田的應(yīng)用[J]. 西南石油學(xué)院學(xué)報(bào), 2006, 28(3):53-56. HOU J, WANG J R, WANG CH SH, et al. Application of the thin & poor reservoir predicted technology to the Punan oilfield. Journal of Southwest Petroleum Institute, 2006, 28(3):53-56.(In Chinese)

[16]劉小平, 楊曉蘭, 劉彥君. 大慶油田葡南地區(qū)葡萄花油層薄層砂巖儲層預(yù)測[J]. 石油物探, 2007, 46(2):162-165. LIU X P, YANG X L, LIU Y J. Thin sand reservoir prediction of Putaohua oil layers in Punan area of Daqing oilfield. Geophysical Prospecting for Petroleum, 2007, 46(2):162-165.(In Chinese)

Seismic prediction of thin sand body reservoir in west margin of Punan oilfield in Songliao basin

YIN Bing-xiang1,QIU Chang-xing1,LIU Hong-tao2

(1.School of geosciences, China University of Petroleum, Qingdao 266580,China;2.Daqing Oil Field Co., Ltd. No.7 Oil Extraction Factory, Daqing,163517,China)

In Punan oilfield, the development layer is Putaohua reservoir of the first member of Yaojia formation, the reservoir include sub-channel branch sand, sheet sand, and offshore sand, which are formed in Delta front sub-face.The reservoir drastic variation laterally results in poor continuity. So that, seismic technique is indispensable in reservoir prediction. In the reservoir prediction of Pu331-1 Block which located at the west margin of Punan oilfield, multiple technique or data are integrated together such as the establishment of fine geological model, well constrained impedance inversion, multiple seismic attribute analysis and logging interpretation result. The predicted reservoir distribution meets the state in this area and fit the data at well site. The study provides data for distribution of oil and gas directly, which is used to locating exploitation well. The technique is also helpful to exploration and exploitation of petroleum in adjacent area.

Punan oilfield; Putaohua oil formation; thin sand reservoir; seismic prediction

2014-08-29 改回日期：2015-03-18

國家重大科學(xué)工程建設(shè)項(xiàng)目計(jì)劃(2011ZX05009-002-3)

尹兵祥(1970-)，男，博士，副教授，從事儲層地震預(yù)測和綜合地球物理教學(xué)與研究工作,E-mail：Yinbx@upc.edu.cn。

1001-1749(2015)04-0494-08

P 631.4

A

10.3969/j.issn.1001-1749.2015.04.13