舒夢(mèng)珵， 張峰， 鄒佳儒， 霍守東*， 羅鳴， 張萬(wàn)棟，李文拓， 穆盛強(qiáng)， 梁瑤

1 中國(guó)石油大學(xué)(北京)地球物理學(xué)院， 北京 102249 2 中國(guó)科學(xué)院地質(zhì)與地球物理研究所， 北京 100029 3 中海石油(中國(guó))有限公司海南分公司， 海南?？?570312 4 中海石油(中國(guó))有限公司湛江分公司， 廣東湛江 524057 5 中國(guó)地質(zhì)科學(xué)院地質(zhì)研究所， 自然資源部深地動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室， 北京 100037

0 引言

南海資源豐富，勘探的開發(fā)前景廣闊，北部南海鶯歌海盆地和瓊東南盆地超過(guò)75%的地區(qū)高溫超壓地層發(fā)育，高溫高壓天然氣資源量大約有4萬(wàn)億m3的資源量,是我國(guó)南海典型的高溫高壓盆地，也是南海大氣區(qū)建設(shè)和國(guó)家南海能源基地建設(shè)重點(diǎn)區(qū)域(陳多福等，2004; 李緒深等，2017).鉆井是超高溫高壓油氣勘探開發(fā)的關(guān)鍵環(huán)節(jié)，是南海資源開發(fā)不可缺少的技術(shù).但由于南海地區(qū)高溫超壓地層的發(fā)育，在鉆井過(guò)程中，受關(guān)鍵層位、尤其是薄弱層影響，較易導(dǎo)致井漏、溢流等復(fù)雜事故，部分井甚至無(wú)法實(shí)現(xiàn)鉆探目標(biāo)的精準(zhǔn)卡層，會(huì)大大增加鉆井風(fēng)險(xiǎn)，以致鉆井事故頻發(fā)，給海上超高溫高壓油氣勘探開發(fā)帶來(lái)巨大安全隱患.

目前，通過(guò)地震數(shù)據(jù)進(jìn)行關(guān)鍵層位及薄弱層識(shí)別是最為有效的手段.但由于常規(guī)地震數(shù)據(jù)分辨率不足，在隨鉆的卡層和風(fēng)險(xiǎn)評(píng)估中無(wú)法精準(zhǔn)識(shí)別特殊地層的分布，亟需高分辨率處理技術(shù)來(lái)提高地震數(shù)據(jù)的縱向分辨率，從而提高地震數(shù)據(jù)識(shí)別風(fēng)險(xiǎn)的能力(李慶忠， 1993).常規(guī)的地震資料高分辨率方法主要有反褶積(章珂等，1999)、時(shí)變譜白化(范小冬等，1995)、反Q濾波(姚振興等，2003；王珺等，2008)和基于反演方法的拓頻處理(曹孟起等，2003).雖然上述方法均取得了實(shí)際應(yīng)用效果，但都存在著一些不足.反褶積法地震數(shù)據(jù)并不能完全滿足地層反射系數(shù)序列的統(tǒng)計(jì)規(guī)律和子波相位的前提假設(shè)，對(duì)高頻噪聲敏感，因此會(huì)影響反褶積方法的應(yīng)用效果.譜白化法沒(méi)有考慮子波在傳播過(guò)程中的時(shí)變特性，得到的地震記錄結(jié)果通常無(wú)法保持原始剖面的反射能量相對(duì)關(guān)系，而且對(duì)于能量弱的同相軸效果較差.反Q濾波法需要求出地下地質(zhì)體的吸收衰減Q因子模型，但實(shí)際應(yīng)用過(guò)程中Q因子很難求取準(zhǔn)確，極大地限制了該方法在實(shí)際資料處理中的應(yīng)用.拓展頻帶方法雖然在原理上可以保持地震數(shù)據(jù)有效的高頻成分，恢復(fù)的薄層信息相對(duì)保真，但其穩(wěn)定性和可靠性受到參數(shù)的影響比較大.常規(guī)的譜反演方法通常在L2范數(shù)的約束條件下進(jìn)行求解，提高地震數(shù)據(jù)分辨率效果有限.針對(duì)上述不足，本文選取基于譜反演的方法進(jìn)行地震數(shù)據(jù)高分辨率處理，通過(guò)利用反射系數(shù)的奇偶分解原理(劉建輝等，2018)及偶分量對(duì)于薄層分辨能力更高等優(yōu)點(diǎn)，通過(guò)反演方法求解目標(biāo)函數(shù)(劉萬(wàn)金，2013)，得到分辨率更高的反射系數(shù).

近年來(lái)迅速發(fā)展的壓縮感知理論(Donoho, 2006; Baraniuk, 2007; Candes and Wakin,2008)給信號(hào)帶來(lái)了新的采集與處理方法.在壓縮感知問(wèn)題中，可以通過(guò)求解一個(gè)范數(shù)的稀疏約束優(yōu)化問(wèn)題重建稀疏信號(hào)，從而達(dá)到保幅、保真提高信號(hào)高頻信息的目的.但該問(wèn)題是一個(gè)典型的NP難題.有學(xué)者為了解決這一問(wèn)題，提出了一系列的迭代閾值類算法.Blumensath等(2007)提出了硬閾值迭代算法.Mallat和Zhang(1992)提出了匹配追蹤(Matching Pursuit，MP)算法，通過(guò)從字典矩陣中選擇與信號(hào)最為匹配的原子構(gòu)建稀疏逼近，求出信號(hào)殘差，然后通過(guò)迭代選擇與殘差最為匹配的原子，直到殘差滿足要求，但該算法收斂需要很多次迭代.Tropp 和Gilbert(2007)針對(duì)匹配追蹤算法的不足，提出了正交匹配追蹤算法，該算法能夠確保殘差與已選擇原子之間的正交性，使得在有限次迭代中達(dá)到收斂的效果.Blumensath和Davies(2008)提出了基于貪婪算法的方向追蹤法，該算法通過(guò)選擇合適的方向和步長(zhǎng)進(jìn)行搜索最優(yōu)解.Needell和Tropp(2009)提出了一種新的迭代重建算法——壓縮采樣匹配追蹤算法(Compressive Sampling Matching Pursuit, CoSaMP)，該算法在達(dá)到最優(yōu)化算法的重建精度的同時(shí)大大降低了算法的計(jì)算復(fù)雜度以及計(jì)算空間.Foucart(2011)在正交匹配追蹤算法和壓縮采樣匹配追蹤算法的基礎(chǔ)之上，提出了一種新的求解欠定方程的迭代算法——硬閾值追蹤算法(Hard Thresholding Pursuit, HTP)，該算法有效提升了計(jì)算效率.該類算法還有正則化的正交匹配追蹤算法(Regularized Orthogonal Matching Pursuit, ROMP)、改進(jìn)的后退型OOMP(Improved backward optimized OMP, IBOOMP)等.本文提出將基于壓縮感知理論應(yīng)用于譜反演算法中以實(shí)現(xiàn)提高地震數(shù)據(jù)分辨率的目的.通過(guò)對(duì)上述稀疏約束優(yōu)化問(wèn)題求解方法對(duì)比，最終本文選用硬閾值法求解方法直接進(jìn)行壓縮感知L0范數(shù)求解，該方法不但可以保持尺度獨(dú)立，而且還可大大提高算法及信號(hào)重建的穩(wěn)定性.

同時(shí)，本文提出了一種穩(wěn)健的壓縮感知譜反演算法及技術(shù)流程：首先利用薄層匹配追蹤算法求取初始反射系數(shù)；在此基礎(chǔ)上，進(jìn)行基于L0范數(shù)壓縮感知理論的疊后稀疏反演，得到最終反射系數(shù)模型；最后基于整形正則化進(jìn)行子波提取，并進(jìn)行小波分解以及高頻子波替換，以達(dá)到保幅、保真地拓展高頻，大幅提高地震數(shù)據(jù)分辨率的目的.將本文方法應(yīng)用于隨鉆過(guò)程中，不但能夠相對(duì)保幅、保真地提高地震數(shù)據(jù)的分辨率，而且還可精細(xì)刻畫超高溫高壓關(guān)鍵層位及薄弱層的地震響應(yīng)特征.

1 基本原理

本文提出了一套基于壓縮感知L0范數(shù)譜反演算法提高地震數(shù)據(jù)分辨率方法和實(shí)現(xiàn)流程，主要分為如下三個(gè)步驟：

(1)基于薄層匹配追蹤算法求取地震數(shù)據(jù)初始反射系數(shù).

(2)在步驟(1)所獲得初始反射系數(shù)的基礎(chǔ)上，進(jìn)行基于壓縮感知L0范數(shù)譜反演計(jì)算，反演求得最終反射系數(shù).

(3)基于整形正則化進(jìn)行子波提取，并進(jìn)行小波分解及高頻子波替換，最終實(shí)現(xiàn)提高分辨率的目的.

1.1 薄層匹配追蹤

薄層匹配追蹤有單層假設(shè)和雙層假設(shè)，為了能夠提高地震數(shù)據(jù)的分辨能力，本文采用雙層假設(shè)匹配追蹤算法.假設(shè)兩個(gè)單反射層位于時(shí)間ta和tb，反射系數(shù)為a和b，對(duì)應(yīng)的初始子波分別為wa和wb，那么，目標(biāo)函數(shù)obj表達(dá)為：

obj=min∑(d-awa-bwb)2,

(1)

其中，d為地震記錄.對(duì)該方程求極值:

(2)

線性方程的解為：

(3)

在初始計(jì)算時(shí)可以設(shè)置零反射系數(shù)或隨機(jī)初始反射系數(shù)模型，通過(guò)薄層匹配追蹤，可以得到初始的反射系數(shù)模型，如圖1所示.

圖1 (a)輸入地震剖面；(b)初始反射系數(shù)

1.2 基于L0范數(shù)壓縮感知的疊后稀疏反演

在獲得初始反射系數(shù)后，文中在求解反射系數(shù)的過(guò)程中引入L0范數(shù)約束，目標(biāo)函數(shù)為

obj=∑(d-a*w)2+λ‖a‖0(a≠0),

(4)

其中，d為疊后地震數(shù)據(jù)，a為反射系數(shù)，w為子波褶積算子矩陣，λ為L(zhǎng)0范數(shù)的權(quán)重，‖a‖0為反射系數(shù)L0范數(shù)，表示非零值的個(gè)數(shù).由于(4)式中L0范數(shù)約束的優(yōu)化問(wèn)題求解為非凸的，實(shí)際計(jì)算較為困難，因此我們采用硬閾值算法(Li et al., 2011)直接求解L0范數(shù)約束的優(yōu)化問(wèn)題.

為對(duì)比不同范數(shù)約束下目標(biāo)函數(shù)求解的效果，我們將(4)中目標(biāo)函數(shù)的L0范數(shù)約束換為L(zhǎng)1范數(shù)約束及L2范數(shù)約束，可得到L1范數(shù)約束及L2范數(shù)約束的優(yōu)化問(wèn)題.采用一道地震數(shù)據(jù)為例，得到三種范數(shù)約束下優(yōu)化問(wèn)題求解的反射系數(shù)反演結(jié)果如圖2所示.對(duì)比結(jié)果表明，相較于L1范數(shù)約束和L2范數(shù)約束，L0范數(shù)結(jié)果更加稀疏，信噪比更高，能夠得到頻帶更寬的地震反射系數(shù)模型.為進(jìn)一步提頻處理提供較好的結(jié)果，隨后將匹配追蹤結(jié)果作為初始輸入進(jìn)行基于L0范數(shù)壓縮感知的疊后稀疏反演，相關(guān)測(cè)試結(jié)果(圖3)表明，相較于匹配追蹤的結(jié)果，通過(guò)L0范數(shù)疊后稀疏反演(圖3c)結(jié)果可以得到更加準(zhǔn)確的反射系數(shù).

圖2 (a) 輸入地震數(shù)據(jù)反射系數(shù)； (b) L0范數(shù)反演； (c) L1范數(shù)反演； (d) L2范數(shù)反演

圖3 (a) 輸入地震剖面； (b) 匹配追蹤結(jié)果； (c) 稀疏反演結(jié)果

1.3 連續(xù)小波變換-Morlet小波

Morlet小波ψ(t)可視為高斯包絡(luò)調(diào)幅下的復(fù)正弦函數(shù)(Morlet et al., 1982)：

(5)

其中σ為Morlet子波的調(diào)制頻率，通常情況下取值需大于5.33.Morlet子波ψ(t)的傅里葉變換為：

(6)

典型的Morlet小波時(shí)間域波形及其振幅譜如圖4所示.隨后采用地震道的傅里葉變換D(ω)和反射系數(shù)傅里葉變換R(ω)優(yōu)化得到Morlet子波的最優(yōu)尺度，可通過(guò)最小化式(7)的優(yōu)化目標(biāo)函數(shù)得到最優(yōu)尺度aopt：

圖4 (a) 時(shí)間域Morlet小波； (b) 頻率域Morlet小波

(7)

(8)

最終利用高頻Morlet小波或者其他高頻子波ψhigh(t)對(duì)ψ(t)進(jìn)行替換并重構(gòu)(如圖5所示即為一個(gè)Morlet子波及一個(gè)高頻子波)，可自然地拓展低、高頻能量：

圖5 (a) 低頻Morlet小波示意圖； (b) 替換用Morlet高頻小波示意圖

(9)

應(yīng)用上述步驟對(duì)圖6a所示數(shù)據(jù)進(jìn)行處理，得到最終寬頻處理結(jié)果如圖6b所示，可明顯提高地震剖面的分辨率.

圖6 提高分辨率處理前后對(duì)比

基于以上三個(gè)步驟，本文形成了基于壓縮感知L0范數(shù)譜反演提高分辨率技術(shù)流程(圖7)，在技術(shù)流程中，為了能夠保證子波計(jì)算過(guò)程中的穩(wěn)定性，采用整形正則化子波提取方法(Huo，2015)，以保證在反射系數(shù)較小區(qū)域或信噪比較低的地震數(shù)據(jù)情況下可提高求取子波的穩(wěn)定性.

圖7 基于壓縮感知L0范數(shù)譜反演提高分辨率技術(shù)流程

2 基于壓縮感知譜反演算法實(shí)際數(shù)據(jù)試算效果

通過(guò)實(shí)驗(yàn)數(shù)據(jù)測(cè)試，合成地震記錄標(biāo)定結(jié)果(圖8)顯示，文中提出的提高分辨方法能夠相對(duì)保幅、保真的恢復(fù)原始數(shù)據(jù)中的高頻信息.原始實(shí)驗(yàn)數(shù)據(jù)應(yīng)用26 Hz雷克子波較好標(biāo)定，對(duì)地層厚度識(shí)別能力有限，提高分辨率結(jié)果能夠應(yīng)用48 Hz雷克子波進(jìn)行很好標(biāo)定，地層厚度識(shí)別能力提高一倍.同時(shí)，拓頻后地震數(shù)據(jù)的波組特征與合成記錄吻合較好(如圖8中黑色實(shí)線及黑色虛線方框范圍顯示)，很好地保留了原始地震數(shù)據(jù)波組強(qiáng)弱對(duì)應(yīng)關(guān)系，并且相較于原始地震數(shù)據(jù)，能夠突出更為精細(xì)的小尺度沉積地層分布規(guī)律.

圖8 提高分辨率前、后井震標(biāo)定對(duì)比圖

另一方面，從測(cè)試區(qū)的沿T31層位均方根屬性(RMS)圖可以看出，提高分辨率處理結(jié)果(圖9b)對(duì)于地層、斷層和水道體的刻畫與原始地震數(shù)據(jù)基本一致，同時(shí)可反映出更多地質(zhì)體的細(xì)節(jié)信息(如圖9b中黑色圓圈顯示區(qū)域).由此說(shuō)明，文中方法在實(shí)現(xiàn)了相對(duì)保幅、保真提高分辨率目的的基礎(chǔ)之上，還能對(duì)地震數(shù)據(jù)中的小尺度地質(zhì)信息進(jìn)行精細(xì)刻畫.

圖9 提高分辨率前(a)、后(b)沿T31層段RMS屬性切片對(duì)比圖

3 隨鉆中實(shí)際資料分析

本文在基于壓縮感知譜反演算法的基礎(chǔ)上，將其應(yīng)用于隨鉆過(guò)程中的精準(zhǔn)卡層識(shí)別，并形成針對(duì)薄弱層精準(zhǔn)卡層識(shí)別及薄儲(chǔ)層精細(xì)深度矯正兩項(xiàng)技術(shù)流程.

3.1 薄弱層精準(zhǔn)卡層識(shí)別

文章選取井C進(jìn)行實(shí)際數(shù)據(jù)試算，井C位于鶯歌海凹陷斜坡帶南段，主要目的層為黃流組一段和二段，該區(qū)存在薄弱層發(fā)育，在鉆井過(guò)程中造成鉆井液漏失的復(fù)雜情況.

原始數(shù)據(jù)顯示(圖10)，該區(qū)薄弱層主要位于區(qū)域地震反射T31(藍(lán)色層)與強(qiáng)波峰(綠色層)位置中間，為復(fù)合波谷或中強(qiáng)波峰反射特征，由于研究區(qū)地震數(shù)據(jù)頻率較低，常規(guī)數(shù)據(jù)處理方法得到的處理成果不能有效地對(duì)超高溫高壓薄弱層進(jìn)行單獨(dú)識(shí)別.

圖10 原始地震數(shù)據(jù)薄弱層空間分布示意圖

利用文中方法對(duì)現(xiàn)有地震數(shù)據(jù)進(jìn)行保幅、保真高分辨率處理，如圖11所示通過(guò)井C合成記錄標(biāo)定后的井旁道分析，原始數(shù)據(jù)復(fù)合波谷處薄層及其空間分布范圍得到有效分離，鉆井證實(shí)該套薄層為薄弱層，該處發(fā)生鉆井液漏失的鉆井復(fù)雜情況.

圖11 提高分辨率前、后井的合成地震記錄對(duì)比圖

文章通過(guò)研究區(qū)多井縱波阻抗和密度曲線的交會(huì)，進(jìn)一步開展巖石物理分析(圖12).通過(guò)巖性的標(biāo)定，確定了圖中紅色方框內(nèi)為砂體阻抗及密度分布范圍.進(jìn)一步，將漏失段多井測(cè)量的阻抗和密度值在交會(huì)圖中標(biāo)出，范圍為圖中紫色方框區(qū)域，揭示了試驗(yàn)區(qū)薄弱漏失層波阻抗的數(shù)值分布范圍(阻抗值9300～11400 m·s-1·g·cm-3).基于此結(jié)果可以半定量的刻畫出工區(qū)內(nèi)實(shí)驗(yàn)井薄弱層漏失的空間分布范圍.

圖12 三口實(shí)驗(yàn)井波阻抗-密度測(cè)井曲線交會(huì)圖

拓頻前后反演結(jié)果對(duì)比(圖13)表明，拓頻后反演結(jié)果能夠很好地反映出薄弱層的反射特征，與已知鉆井巖性信息吻合；拓頻后反演數(shù)據(jù)縱向分辨率提高，橫向連續(xù)性增強(qiáng)，通過(guò)反演的波阻抗數(shù)據(jù)標(biāo)定，大大提高了砂體薄弱層及漏失區(qū)域的可追蹤性及可識(shí)別性.

圖13 提高分辨率前(a)、后(b)過(guò)井反演剖面對(duì)比

通過(guò)上述的分析，文中總結(jié)了高精度薄弱層識(shí)別卡層的技術(shù)流程如圖14所示.

圖14 高精度薄弱層識(shí)別技術(shù)流程

3.2 薄儲(chǔ)層精細(xì)深度矯正

文中選取鶯歌海盆地樂(lè)東區(qū)中央底辟帶的測(cè)試數(shù)據(jù)，進(jìn)行薄儲(chǔ)層精細(xì)刻畫及深度矯正實(shí)際數(shù)據(jù)計(jì)算，測(cè)試目標(biāo)區(qū)主要目的層為黃流組一段，發(fā)育多套薄砂體沉積，現(xiàn)有地震數(shù)據(jù)難以精確識(shí)別砂體分布.井D為該區(qū)一口探井，在隨鉆過(guò)程中對(duì)砂體深度的精準(zhǔn)卡層預(yù)測(cè)，可以有效指導(dǎo)施工方案調(diào)整，為鉆井安全實(shí)施提供依據(jù).

原始數(shù)據(jù)顯示(圖15a)該區(qū)砂體儲(chǔ)層主要位于區(qū)域地震反射T30(黃色層)下覆，為復(fù)合波反射特征.由于研究區(qū)地震數(shù)據(jù)頻率較低，常規(guī)數(shù)據(jù)處理方法得到的成果剖面無(wú)法對(duì)目標(biāo)井黃流組一段(N1hl1-B至N1hl1-C)中三套主力砂體進(jìn)行刻畫.而提高分辨率數(shù)據(jù)結(jié)果(圖15b)可以清楚的突顯出復(fù)合波組中薄砂體的空間分布.

圖15 提高分辨率前(a)、后(b)目標(biāo)井D過(guò)井地震偏移剖面對(duì)比圖

基于提高分辨率數(shù)據(jù)進(jìn)行地震反演(圖16)可以看到，原始數(shù)據(jù)的反演結(jié)果無(wú)法對(duì)井D的三套主力砂體有效刻畫(圖16a，暖色為砂體波阻抗特征)，而基于本文提高分辨率數(shù)據(jù)得到的高精度反演結(jié)果能夠清楚地反映三套砂體的空間分布及相互間的交接關(guān)系(圖16b，暖色為砂體波阻抗特征)，其中可識(shí)別砂體的最小厚度為二號(hào)砂體13 m.

圖16 提高分辨率前(a)、后(b)目標(biāo)井D過(guò)井波阻抗反演剖面對(duì)比圖

另外，在T30反射層的上覆，與原始地震數(shù)據(jù)相比，提高分辨率之后的數(shù)據(jù)能夠很好的突出鶯歌海組二段N2hgh2-A層和N2hgh2-B層的反射特征，從剖面上看(圖17b)，這兩層在提高分辨率剖面上表現(xiàn)為中弱振幅連續(xù)反射特征，全區(qū)可追蹤.由于原始數(shù)據(jù)在T30反射層前無(wú)法進(jìn)行深度標(biāo)定卡層，在鉆遇T30時(shí)，鉆遇地層深度與預(yù)測(cè)深度的誤差為20 m，而通過(guò)提高分辨率數(shù)據(jù)，可以在鉆遇T30前，通過(guò)N2hgh2-A層和N2hgh2-B層進(jìn)行深度標(biāo)定，在鉆遇T30儲(chǔ)層頂時(shí)，鉆遇地層深度與預(yù)測(cè)深度的誤差僅為8 m，預(yù)測(cè)誤差為2.2‰(圖17c).

圖17 提高分辨率前(a)、后(b)—(c)目標(biāo)井D井深度識(shí)別誤差對(duì)比圖

通過(guò)隨鉆過(guò)程中對(duì)砂體深度的精準(zhǔn)卡層預(yù)測(cè)的應(yīng)用，文章提出了高精度薄儲(chǔ)層識(shí)別及精準(zhǔn)卡層技術(shù)流程(圖18)，實(shí)現(xiàn)了隨鉆過(guò)程中的精準(zhǔn)卡層及薄儲(chǔ)層識(shí)別，為鉆井過(guò)程中的風(fēng)險(xiǎn)評(píng)估和鉆井方案的調(diào)整提供了可靠依據(jù).

圖18 高精度薄儲(chǔ)層識(shí)別及精準(zhǔn)地層識(shí)別技術(shù)流程

4 結(jié)論與展望

本文系統(tǒng)的論述了目前主流的提高分辨率處理方法和壓縮感知算法的研究現(xiàn)狀及方法優(yōu)勢(shì)，提出了基于壓縮感知譜反演算法提高分辨率技術(shù)及計(jì)算流程，并將該技術(shù)應(yīng)用于隨鉆過(guò)程精準(zhǔn)卡層及風(fēng)險(xiǎn)層評(píng)估的應(yīng)用中，取得了較好效果.

文中提出的基于壓縮感知L0范數(shù)譜反演算法，可以較好的實(shí)現(xiàn)反射系數(shù)的求解，對(duì)弱反射系數(shù)有更好的保護(hù)性.同時(shí)，基于該算法的計(jì)算流程，可以有效的保證L0范數(shù)求解反射系數(shù)過(guò)程中的計(jì)算穩(wěn)定性，提高計(jì)算迭代效率，實(shí)現(xiàn)相對(duì)保幅、保真的提高地震數(shù)據(jù)分辨率的目的.

該算法在超高溫高壓隨鉆過(guò)程的實(shí)際數(shù)據(jù)計(jì)算中，通過(guò)對(duì)實(shí)際地震數(shù)據(jù)保幅、保真高分辨率處理，有效的識(shí)別了隨鉆過(guò)程中鉆遇的關(guān)鍵層位(尤其是薄弱層)的地震響應(yīng)特征及空間分布規(guī)律.并且，文中針對(duì)薄弱層和薄儲(chǔ)層精準(zhǔn)卡層的情況提出了針對(duì)性的計(jì)算流程，在實(shí)際數(shù)據(jù)應(yīng)用中取得了較好的應(yīng)用效果.

該方法在隨鉆過(guò)程中的應(yīng)用，可以有效避免因?yàn)楸∪鯇与y以識(shí)別及關(guān)鍵層位卡層不準(zhǔn)而引起的鉆井事故和經(jīng)濟(jì)損失，為南海西部油田超高溫高壓鉆井的安全生產(chǎn)及鉆井方案的調(diào)整提供保障.

致謝感謝審稿人和編輯對(duì)本文提出的寶貴修改建議.