李 翔, 尹 陳, 吳永宏, 熊曉軍, 劉 陽(yáng), 簡(jiǎn)世凱

(1.油氣藏地質(zhì)及開發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室 (成都理工大學(xué))，成都 610059;2.中國(guó)石油川慶鉆探工程有限公司 地球物理勘探公司，成都 610213 3.中國(guó)石油西南油氣田分公司 重慶氣礦，重慶 400021)

基于Tsvankin參數(shù)頁(yè)巖儲(chǔ)層各向異性強(qiáng)度預(yù)測(cè)

李 翔1, 尹 陳2, 吳永宏3, 熊曉軍1, 劉 陽(yáng)1, 簡(jiǎn)世凱1

(1.油氣藏地質(zhì)及開發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室 (成都理工大學(xué))，成都 610059;2.中國(guó)石油川慶鉆探工程有限公司 地球物理勘探公司，成都 610213 3.中國(guó)石油西南油氣田分公司 重慶氣礦，重慶 400021)

在頁(yè)巖氣勘探開發(fā)過(guò)程中，各向異性強(qiáng)弱與儲(chǔ)層品質(zhì)、地震波速度、巖石脆性、巖石彈性參數(shù)等都有關(guān)系，所以對(duì)頁(yè)巖各向異性特征的定量描述具有重要意義。綜合考慮頁(yè)巖頁(yè)片狀層理大量發(fā)育，基于VTI介質(zhì)模型對(duì)頁(yè)巖進(jìn)行模擬，通過(guò)對(duì)Thomsen參數(shù)簡(jiǎn)化而利用Tsvankin參數(shù)進(jìn)行頁(yè)巖各向異性表征，然后以疊前偏移結(jié)果為基準(zhǔn)，修正均方根速度使得疊前偏移后的CRP道集層拉平效果明顯，迭代計(jì)算出使得疊前偏移效果最優(yōu)時(shí)的各向異性強(qiáng)度的Tsvankin參數(shù)值，并利用Tsvankin參數(shù)對(duì)疊前道集剩余時(shí)差進(jìn)行校正，最后輸出Tsvankin參數(shù)數(shù)據(jù)體，分析Tsvankin參數(shù)與頁(yè)巖層脆性發(fā)育及儲(chǔ)層性質(zhì)的關(guān)系。將該方法應(yīng)用于四川長(zhǎng)寧“寧201”井區(qū)的下志留統(tǒng)龍馬溪組頁(yè)巖，綜合地震資料分析發(fā)現(xiàn)Tsvankin參數(shù)譜值域范圍與頁(yè)巖脆性概率分布趨勢(shì)正相關(guān)，能夠預(yù)測(cè)頁(yè)巖儲(chǔ)層各向異性強(qiáng)度并能反映頁(yè)巖地層脆性強(qiáng)弱，在“甜點(diǎn)”預(yù)測(cè)中有一定的指導(dǎo)意義。

Tsvankin參數(shù)；頁(yè)巖；各向異性；脆性

頁(yè)巖氣賦存于具有超低孔隙度的頁(yè)巖中，是一種在海相、陸相及過(guò)渡相地層均有發(fā)育的很重要的非常規(guī)資源。頁(yè)巖由于其構(gòu)成礦物顆粒在成分、形狀及大小上均有差異，加之結(jié)構(gòu)上大量發(fā)育層理及片理，所以頁(yè)巖儲(chǔ)層各向異性特征十分明顯。四川盆地作為中國(guó)最重要的頁(yè)巖氣勘探區(qū)域，探究各向異性特征對(duì)頁(yè)巖儲(chǔ)層“甜點(diǎn)”預(yù)測(cè)及后期開發(fā)意義重大。

為了探究頁(yè)巖儲(chǔ)層各向異性發(fā)育情況，鄧?yán)^新等[1]在實(shí)驗(yàn)室超聲條件下進(jìn)行試驗(yàn)，探究了泥、頁(yè)巖縱橫波速度及各向異性參數(shù)隨壓力的變化過(guò)程；Carl等[2]探究了Thomsen參數(shù)的適用性并結(jié)合SEM資料研究了頁(yè)巖微觀物性對(duì)其各向異性特征的影響；王倩等[3]根據(jù)巖石物理實(shí)驗(yàn)探究了頁(yè)巖彈性模量、泊松比與層理面的關(guān)系，設(shè)計(jì)了一種頁(yè)巖各向異性指數(shù)；謝劍勇等[4]基于巖樣測(cè)試研究了頁(yè)巖各向異性參數(shù)間的關(guān)系并對(duì)影響各向異性的因素進(jìn)行分析。這些研究證明了頁(yè)巖由于其層理發(fā)育可利用橫向各向同性(TI)介質(zhì)進(jìn)行模擬，在巖石物理領(lǐng)域，通過(guò)Thomsen參數(shù)評(píng)價(jià)頁(yè)巖各向異性特征具有較好的適用性[5]。

以上學(xué)者從實(shí)驗(yàn)的角度對(duì)頁(yè)巖各向異性特征進(jìn)行研究；然而在頁(yè)巖氣勘探過(guò)程中，巖心資料珍貴而數(shù)量極少，實(shí)驗(yàn)分析較難以滿足生產(chǎn)需求。如何從理論計(jì)算的角度進(jìn)行各向異性參數(shù)預(yù)測(cè)是一個(gè)亟待解決的問(wèn)題。張廣智等[6]基于測(cè)井資料，采用巖石物理建模方法，反演計(jì)算了碳酸鹽巖地區(qū)橫波速度及各向異性參數(shù)；孫偉家等[7]從正演的角度對(duì)VTI介質(zhì)強(qiáng)各向異性特征進(jìn)行數(shù)值模擬。而如何充分利用地震資料進(jìn)行各向異性參數(shù)預(yù)測(cè)的研究成果較少。本文基于頁(yè)巖工區(qū)疊前地震資料，設(shè)計(jì)了一套針對(duì)頁(yè)巖儲(chǔ)層的Tsvankin各向異性參數(shù)計(jì)算方法。該參數(shù)能夠?qū)?yè)巖工區(qū)各向異性強(qiáng)度進(jìn)行定量描述，從而對(duì)勘探開發(fā)有一定的指導(dǎo)意義[8-16]。

1 方法原理

1.1 各向異性強(qiáng)度計(jì)算方法

頁(yè)巖由粒徑lt;0.003 9 mm的細(xì)粒碎屑、黏土及有機(jī)質(zhì)組成，頁(yè)巖中廣泛發(fā)育頁(yè)狀及薄片狀的層理，是一種典型的各向異性介質(zhì)。我們利用表征薄層旋回沉積的VTI介質(zhì)來(lái)模擬頁(yè)巖儲(chǔ)層頁(yè)片狀層理，VTI介質(zhì)的彈性剛度矩陣為

(1)

VTI介質(zhì)及HTI介質(zhì)的彈性剛度矩陣均含5個(gè)獨(dú)立的彈性參數(shù)。為便于理論研究，L.Thomsen[17]于1986年提出了一套新的彈性參數(shù)表達(dá)式，將5個(gè)彈性參數(shù)簡(jiǎn)化為3個(gè)，定義如下

(2)

其中：ε值與介質(zhì)的縱波各向異性正相關(guān)，當(dāng)ε=0時(shí)，介質(zhì)的縱波無(wú)各向異性;δ是連接vP0和vP90(0和90表示入射角為0°和90°)之間的一種隨入射角變化而變化的過(guò)渡參數(shù)。

綜合考慮基于常規(guī)聲波地震資料的縱波勘探方法，T.Alkhalifah等[18]重新定義了一個(gè)有效的各向異性參數(shù)“η”(內(nèi)含Thomsen參數(shù))，將彈性參數(shù)進(jìn)一步簡(jiǎn)化為

(3)

在常規(guī)各向同性介質(zhì)速度分析過(guò)程中，反射波時(shí)距曲線滿足雙曲線特征。TI介質(zhì)中，時(shí)距曲線不再滿足簡(jiǎn)單的雙曲線規(guī)律，由T.Alkhalifah等[18]重新定義的時(shí)距方程

(4)

式中：x為偏移距；t0為零偏移距的雙程旅行時(shí)；vNM0為入射角為0°時(shí)的動(dòng)校正速度。(4)式充分地描述了長(zhǎng)距離傳播的縱波正常時(shí)差與零入射角時(shí)的動(dòng)校正速度vNM0及各向異性參數(shù)η的關(guān)系。方程右端第三項(xiàng)不同于常規(guī)的雙曲線時(shí)距方程，其主要控制時(shí)差方程中的非雙曲線形態(tài)，若該式中η=0則該方程與常規(guī)雙曲線時(shí)距方程一致。分母中增加x校正因子，是為了將泰勒展開式中的截?cái)嗾`差進(jìn)行彌補(bǔ)。由(4)式能夠看出，當(dāng)介質(zhì)各向異性較強(qiáng)且偏移距較大時(shí)，反射波時(shí)距關(guān)系會(huì)偏離各向同性雙曲線關(guān)系。

由此，計(jì)算各向異性η參數(shù)可以通過(guò)進(jìn)行速度分析得到均方根速度，并依此速度對(duì)地震數(shù)據(jù)進(jìn)行疊前時(shí)間偏移。由于存在各向異性段的道集，其同相軸形態(tài)受(4)式中的η值影響，形態(tài)不會(huì)是標(biāo)準(zhǔn)的雙曲線；因而在常規(guī)疊前偏移后其層拉平效果不明顯。當(dāng)偏移迭代計(jì)算結(jié)果變化不大且層拉平效果不明顯時(shí)，基于(4)式引入η參數(shù)重新計(jì)算地震波時(shí)差，通過(guò)修改η值在CMP道集動(dòng)校拉平效果最優(yōu)時(shí)輸出此時(shí)的各向異性η參數(shù)作為該道的各向異性強(qiáng)度檢測(cè)值。

1.2 各向異性參數(shù)計(jì)算流程

綜合考慮上述方法，本文設(shè)計(jì)了一套針對(duì)頁(yè)巖工區(qū)疊前各向異性強(qiáng)度計(jì)算流程，其具體步驟如下。

步驟1： 基于頁(yè)巖氣地區(qū)三維疊前地震數(shù)據(jù)，對(duì)數(shù)據(jù)進(jìn)行疊前道集保邊去噪處理。通過(guò)引入最小標(biāo)準(zhǔn)差準(zhǔn)則，在壓制資料噪聲的同時(shí)能很好地保留邊緣等細(xì)節(jié)信息。在實(shí)際資料處理過(guò)程中發(fā)現(xiàn)該步驟對(duì)計(jì)算結(jié)果影響顯著，基于保邊去噪后品質(zhì)較好的地震數(shù)據(jù)能計(jì)算得到更加穩(wěn)定的η參數(shù)。

步驟2： 對(duì)疊前道集進(jìn)行速度分析，計(jì)算RMS速度(均方根速度)，基于得到的RMS速度進(jìn)行疊前時(shí)間偏移，并觀察偏移后結(jié)果的道集拉平效果，可反復(fù)迭代計(jì)算至偏移后疊前道集位置變化不大時(shí)結(jié)束計(jì)算，作為計(jì)算η參數(shù)的基礎(chǔ)。

步驟3： 給定初始的各向異性參數(shù)η場(chǎng)，將η值代入(4)式重新計(jì)算地震波走時(shí)。利用計(jì)算得到的非雙曲線形態(tài)地震波走時(shí)對(duì)疊前道集剩余時(shí)差進(jìn)行修正，觀察修正后疊前道集層拉平效果，當(dāng)拉平效果滿足要求時(shí)輸出。迭代計(jì)算完成整塊三維疊前數(shù)據(jù)得到三維η數(shù)據(jù)體，作為工區(qū)各向異性強(qiáng)度預(yù)測(cè)結(jié)果。

圖1為上述步驟的技術(shù)流程。

圖1 各向異性參數(shù)η場(chǎng)計(jì)算流程Fig.1 The flowchart of calculation of anisotropic parameter field

通過(guò)以上3個(gè)步驟計(jì)算得到頁(yè)巖工區(qū)的各向異性參數(shù)η數(shù)據(jù)體。由此對(duì)η數(shù)據(jù)體進(jìn)行分析，繪制研究區(qū)目的層η平面展布圖，可結(jié)合地震資料分析各向異性強(qiáng)度與地震屬性的聯(lián)系，并基于η參數(shù)進(jìn)一步實(shí)現(xiàn)各向異性疊前偏移。

2 應(yīng)用實(shí)例

四川長(zhǎng)寧地區(qū)頁(yè)巖氣資源豐富，是中國(guó)石油天然氣股份有限公司頁(yè)巖氣建產(chǎn)核心區(qū)域之一[19]。其中寧201-H1井日產(chǎn)氣gt;0.15×106m3，其勘探目的層為下志留統(tǒng)龍馬溪組頁(yè)巖。

將本文提出的各向異性強(qiáng)度預(yù)測(cè)方法應(yīng)用在長(zhǎng)寧的“寧201”井區(qū)三維疊前地震數(shù)據(jù)中，基于常規(guī)疊前時(shí)間偏移結(jié)果為約束，迭代計(jì)算使得層拉平效果最優(yōu)時(shí)的各向異性參數(shù)譜，并根據(jù)得到的η譜對(duì)疊前數(shù)據(jù)剩余時(shí)差進(jìn)行校正。工區(qū)實(shí)踐如圖2所示，圖2-A中計(jì)算所得的η參數(shù)譜能量越集中且偏離中線越遠(yuǎn)表示各向異性發(fā)育強(qiáng)度越大，圖2-B為未經(jīng)各向異性參數(shù)校正的疊前道集，圖2-C為經(jīng)各向異性校正后的疊前道集。對(duì)比發(fā)現(xiàn)：經(jīng)過(guò)各向異性參數(shù)校正后，淺層遠(yuǎn)偏移距處的層拉平效果改善明顯。與此同時(shí)，圖2-A中拾取的η參數(shù)譜在1.4 s處偏離中線位置，且能量較強(qiáng)，對(duì)應(yīng)各向異性參數(shù)校正前后的CRP道集剖面；如圖2-B和圖2-C中紅色箭頭所指，疊前道集同相軸的層拉平效果有一定程度的改善，說(shuō)明各向異性參數(shù)對(duì)CRP道集剩余時(shí)差校正有積極作用，計(jì)算得到的各向異性參數(shù)數(shù)據(jù)體對(duì)后期的疊前偏移意義同樣重大。

圖2 各向異性參數(shù)對(duì)疊前CRP道集校正結(jié)果Fig.2 The results of anisotropic parameter for the pre-stack gathers(A)各向異性參數(shù)譜剖面; (B)各向異性參數(shù)校正前的疊前道集; (C)各向異性參數(shù)校正后的疊前道集

對(duì)三維工區(qū)計(jì)算其各向異性參數(shù)體，針對(duì)目的層繪制平面圖(圖3)。我們發(fā)現(xiàn)“寧201”井處于各向異性強(qiáng)度中高值段(圖中紅黃色段)，與實(shí)鉆高產(chǎn)氣結(jié)論吻合，說(shuō)明各向異性發(fā)育強(qiáng)度與頁(yè)巖儲(chǔ)層品質(zhì)有關(guān)，可以一定程度上指導(dǎo)儲(chǔ)層“甜點(diǎn)”預(yù)測(cè)。龍馬溪組頁(yè)巖脆性概率呈西南高東北低的趨勢(shì)，這與各向異性參數(shù)數(shù)據(jù)強(qiáng)弱趨勢(shì)范圍完全一致，說(shuō)明各向異性強(qiáng)度與頁(yè)巖脆性發(fā)育存在一定程度的正相關(guān)關(guān)系，圖中紅黃色段是脆性概率較大區(qū)域，是有利的“甜點(diǎn)”區(qū)。

圖3 沿目的層向下50 ms各向異性參數(shù)平面圖Fig.3 The layout of anisotropy parameter downward for 50 ms along the target horizon

工區(qū)正東Inline 476線(圖3中黑線)疊后地震剖面如圖4，目的層段龍馬溪組(圖中藍(lán)線)地震同相軸連續(xù)性較差，可以觀察到小斷裂較為發(fā)育，說(shuō)明此處地層普遍具有各向異性，與根據(jù)地震剩余時(shí)差計(jì)算所得的高各向異性參數(shù)結(jié)論吻合，該處頁(yè)巖層物性較差，為非有效儲(chǔ)層段。將Inline 476線投影在各向異性參數(shù)平面圖上，可以看到其處于強(qiáng)各向異性發(fā)育段(圖3中白色區(qū)域高值區(qū))，證明各向異性參數(shù)對(duì)評(píng)價(jià)頁(yè)巖儲(chǔ)層物性有一定意義。

圖4 目的層段龍馬溪組頁(yè)巖Inline 476線地震剖面Fig.4 The seismic profile along Inline 476 for Longmaxi Formation

3 結(jié) 論

a.頁(yè)巖中廣泛發(fā)育各向異性，利用Tsvankin各向異性參數(shù)可以獨(dú)立地表征縱波在頁(yè)巖中傳播時(shí)各向異性特征強(qiáng)弱，通過(guò)迭代可以得到一個(gè)最優(yōu)各向異性參數(shù)。該參數(shù)能夠進(jìn)一步校正常規(guī)疊前偏移后的剩余時(shí)差，為精細(xì)刻畫頁(yè)巖斷裂及儲(chǔ)層性質(zhì)提供依據(jù)。

b.得到的各向異性參數(shù)數(shù)據(jù)體可以為頁(yè)巖工區(qū)“甜點(diǎn)”預(yù)測(cè)提供一定的依據(jù)。將該方法應(yīng)用到長(zhǎng)寧地區(qū)，分析認(rèn)為“寧201”井位于各向異性參數(shù)中高值區(qū)域，較優(yōu)的儲(chǔ)層發(fā)育著中等強(qiáng)度的各向異性特征。龍馬溪組頁(yè)巖脆性概率呈西南高東北低的趨勢(shì)，與各向異性參數(shù)的場(chǎng)數(shù)據(jù)強(qiáng)弱趨勢(shì)一致并呈現(xiàn)正相關(guān)關(guān)系。在頁(yè)巖氣勘探中各向異性較強(qiáng)的區(qū)域其巖石脆性較強(qiáng)，是有利的“甜點(diǎn)”區(qū)。

c.在工區(qū)正東部各向異性場(chǎng)很強(qiáng)，呈片狀分布，結(jié)合地震資料分析得出：高各向異性可能由斷裂廣泛發(fā)育引起，地震同相軸連續(xù)性較差的地區(qū)，是勘探中應(yīng)當(dāng)回避的非“甜點(diǎn)”區(qū)。

[1] 鄧?yán)^新,史謌,劉瑞珣,等. 泥巖、頁(yè)巖聲速各向異性及其影響因素分析[J].地球物理學(xué)報(bào),2004,47(5):862-868.

Deng J X, Shi G, Liu R X,etal. Analysis of the velocity anisotropy and its affection factors in shale and mudstone[J]. Chinese Journal of Geophysics, 2004, 47(5): 862-868. (in Chinese)

[2] Carl H, Chandra S. Elastic anisotropy of shales[J]. The Leading Edge, 2011, 30(3): 324-331.

[3] 王倩,王鵬,項(xiàng)德貴,等.頁(yè)巖力學(xué)參數(shù)各向異性研究[J].天然氣工業(yè),2012,32(12):62-65.

Wang Q, Wang P, Xiang D G,etal. Anisotropy property of mechanical parameters of shales[J]. Natural Gas Industry, 2012, 32(12): 62-65. (in Chinese)

[4] 謝劍勇,魏建新,狄?guī)妥?等.頁(yè)巖各向異性參數(shù)間關(guān)系及其影響因素分析[J].石油地球物理勘探,2015,50(6):1141-1145.

Xie J Y, Wei J X, Di B R,etal. Correlation of shale anisotropy parameters and their influence factors[J]. Oil Geophysical Prospecting, 2015, 50(6):1141-1145. (in Chinese)

[5] Mavko G, Mukerji T, Dvorkin J. The Rock Physics Handbook: Tools for Seismic Analysis of Porous Media[M]. Cambridge: Cambridge University Press, 2009.

[6] 張廣智,陳懷震,王琪,等.基于碳酸鹽巖裂縫巖石物理模型的橫波速度和各向異性參數(shù)預(yù)測(cè)[J].地球物理學(xué)報(bào),2013,56(5):1707-1715.

Zhang G Z, Chen H Z, Wang Q,etal. Prediction for in-situ formation stress of shale based on rock physics equivalent model[J]. Chinese Journal of Geophysics, 2013, 56(5): 1707-1715. (in Chinese)

[7] 孫偉家,符力耘,管西竹,等.頁(yè)巖氣地震勘探中頁(yè)巖各向異性的地震模擬研究[J].地球物理學(xué)報(bào),2013,56(3):961-970.

Sun W J, Fu L Y, Guan X Z,etal. A study on anisotropy of shale using seismic forward modeling in shale gas exploration[J]. Chinese Journal of Geophysics, 2013, 56(3): 961-970. (in Chinese)

[8] Meadows M, Coen S. Exact inversion of plane layered isotropic and anisotropic elastic media by the state-space approach [J]. Geophysics, 1986, 51(11): 2031-2050.

[9] Meadows M A, Abriel W L. 3-D poststack phase-shift migration in transversely isotropic media[J]. Seg Technical Program Expanded Abstracts, 1994, 13(1): 1205.

[10] Sena A G. Seismic travel-time equations for azimuthally anisotropic and isotropic media: Estimation of interval elastic properties [J]. Geophysics, 1991, 56(12): 2090-2101.

[11] 孫長(zhǎng)贊,王建民,馮文霞,等.各向異性疊前時(shí)間偏移技術(shù)在大慶探區(qū)的應(yīng)用[J].石油地球物理勘探, 2010, 45(增刊1): 71-73.

Sun C Z, Wang J M, Feng W X,etal. Application of anisotropic pre-stack time migration technique in Daqing exploration area[J]. Oil Geophysical Prospecting, 2010, 45(S1): 71-73. (in Chinese)

[12] 萬(wàn)志超,滕吉文,張秉明.各向異性介質(zhì)中地震波速度分析的研究現(xiàn)狀[J].地球物理學(xué)進(jìn)展,1997,12(3):35-44.

Wan Z C, Teng J W, Zhang B M. Velocity analysis in anisotropic media[J]. Chinese Journal of Geophysics, 1997, 12(3): 35-44. (in Chinese)

[13] Johnston J E, Christensen N I. Seismic anisotropy of shales[J]. Journal of Geophysical Research Solid Earth, 1995, 100(B4): 5991-6003.

[14] Sayers C M. The effect of kerogen on the elastic anisotropy of organic-rich shales[J]. Geophysics, 2013, 78: 65-74.

[15] Alkhalifah T. Velocity analysis for transversely isotropic media[J]. Seg Technical Program Expanded Abstracts, 1994, 60(9): 1550-1566.

[16] 王康寧,孫贊東,侯昕曄.基于傅里葉級(jí)數(shù)展開的縱波方位各向異性裂縫預(yù)測(cè)[J].石油物探,2015,54(6):755-761.

Wang K N, Sun Z D, Hou X Y. Fracture prediction by P-wave azimuthal anisotropy based on Fourier series decomposition[J]. Geophysical Prospecting For Petroleum, 2015, 54(6): 755-761. (in Chinese)

[17] Thomsen L. Weak elastic anisotropy [J]. Geophysics, 1986, 51(10): 1954-1966.

[18] Alkhalifah T, Tsvankin I. Velocity analysis for transversely isotropic media[J]. Geophysics, 1995, 60: 1550-1566.

[19] 蘭潔.西南油氣田四川長(zhǎng)寧區(qū)塊頁(yè)巖氣產(chǎn)量超6000×104m3[J].天然氣與石油,2014,32(5):40.

Lan J. Southwest oil and gas fields Sichuan Chang-ning block shale gas production over 6000×104m3[J]. Natural Gas and Oil, 2014, 32(5): 40. (in Chinese)

PredictionmethodforshalereservoiranisotropybasedontheTsvankinparameter

LI Xiang1, YIN Chen2, WU Yonghong3, XIONG Xiaojun1, LIU Yang1, JIAN Shikai1

1.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,ChengduUniversityofTechnology,Chengdu610059,China;2.GeophysicalExplorationCompany,ChuanqingDrillingEngineeringCo.Ltd.,CNPC,Chengdu610213,China;3.ChongqingGasField,SouthwestOilamp;GasFieldCompary,CNPC,Chongqing420021,China

The anisotropy of rock is related to seismic velocity, brittleness, rock elastic parameters and so on. Quantitative description of anisotropy for shale is very important in the early exploration and post-development process of shale gas. VTI media model is used to simulate the shale depending on the constraint with pre-stack migration results and iteratively calculate the anisotropic strength Tsvankin parameter, and the Tsvankin parameter is used to correct the remaining time difference in pre-stack gather. Finally, the output volume of the Tsvankin parameter data is obtained. The anisotropy parameter from Longmaxi Formation shale is used to discuss the significance of shale anisotropy of reservoir exploration and to describe the shale brittleness. It reveals that the Tsvankin parameter of the Longmaxi Formation shale has positive correlation with probability distribution of shale brittleness. It also can be used to predict the intensity of shale anisotropy and reflect the strong and weak of formation brittleness.

Tsvankin parameter; shale; anisotropy; brittleness

P631.4; TE132

A

10.3969/j.issn.1671-9727.2017.06.11

1671-9727(2017)06-0738-06

2016-11-30。

國(guó)家科技重大專項(xiàng)(2016ZX05023004-001-003); 國(guó)家自然科學(xué)基金項(xiàng)目(41274130)。

李翔(1993-)，男，碩士研究生，研究方向：地球物理勘探, E-mail：lixiang_cdut@163.com。

熊曉軍(1980-)，男，博士，教授，研究方向：地球物理勘探, E-mail:xiongxiaojun07@cdut.cn。