本發(fā)明涉及地球物理勘探數(shù)據(jù)處理技術(shù)領(lǐng)域，特別是涉及到一種適用于測(cè)井約束反演的測(cè)井縱波阻抗去壓實(shí)處理方法。

背景技術(shù)：

測(cè)井約束反演是目前儲(chǔ)層預(yù)測(cè)中的一種重要方法，應(yīng)用范圍較為廣泛。在濁積巖、砂礫巖體和灘壩砂巖等巖性、構(gòu)造-巖性油藏勘探開(kāi)發(fā)中發(fā)揮了重要的作用。通常意義下的測(cè)井約束反演所得到的結(jié)果是地震的縱波阻抗反演體，具有明確的物理意義，在儲(chǔ)層預(yù)測(cè)中，通過(guò)結(jié)合儲(chǔ)層在測(cè)井縱波阻抗曲線上的分布特征來(lái)指導(dǎo)整個(gè)工區(qū)的儲(chǔ)層預(yù)測(cè)。但大量的實(shí)際應(yīng)用表明測(cè)井約束反演所得到的縱波阻抗反演體在儲(chǔ)層預(yù)測(cè)中有時(shí)存在著較大多解性，主要表現(xiàn)為儲(chǔ)層與圍巖在縱波阻抗分布區(qū)間上存在較大的重疊區(qū)間，難以有效區(qū)分。造成這一現(xiàn)象的根本原因在于，在原始的測(cè)井縱波阻抗曲線上，儲(chǔ)層與圍巖之間本身就存在著較大的重疊區(qū)間，因而在測(cè)井約束下所得到的地震縱波阻抗反演體也必然存在這一問(wèn)題。針對(duì)這一問(wèn)題，現(xiàn)有的解決方法是采取聲波曲線重構(gòu)的方法來(lái)放大儲(chǔ)層與圍巖之間的速度差異，具體做法是在原始的聲波曲線中引入自然伽馬、電阻率、自然電位等此類對(duì)儲(chǔ)層響應(yīng)敏感的測(cè)井曲線，通過(guò)數(shù)學(xué)擬合公式，將這些原始測(cè)井曲線重新擬合為一條新的曲線，稱為“偽聲波”曲線，再與密度曲線相乘，得到一條“偽縱波阻抗”曲線，從而參與到最終的測(cè)井約束反演中。分析認(rèn)為，這種聲波曲線重構(gòu)方法所得到的“偽聲波”曲線雖然能放大儲(chǔ)層與圍巖之間的速度差異，但采取這種“偽聲波”曲線最終得到的縱波阻抗反演體的預(yù)測(cè)結(jié)果并不準(zhǔn)確，原因在于采取這種重構(gòu)方式得到的“偽聲波”曲線會(huì)嚴(yán)重的影響到由此而最終計(jì)算得到的反射系數(shù)曲線，這種曲線處理方式會(huì)明顯改變測(cè)井的反射系數(shù)序列，致使其與子波褶積時(shí)所得到的合成記錄道與原始地震道之間的匹配性變差而影響到最終的反演結(jié)果。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的是解決了以上技術(shù)問(wèn)題，提供一種適用于測(cè)井約束反演的測(cè)井縱波阻抗去壓實(shí)處理方法，通過(guò)對(duì)原始測(cè)井的縱波阻抗曲線進(jìn)行合理的處理，從而有效地放大儲(chǔ)層與圍巖之間縱波阻抗差異，使得最終的縱波阻抗反演體能夠準(zhǔn)確地預(yù)測(cè)儲(chǔ)層的分布情況。

本發(fā)明的目的通過(guò)如下的技術(shù)方案來(lái)實(shí)現(xiàn)：

本適用于測(cè)井約束反演的測(cè)井縱波阻抗去壓實(shí)處理方法包括：

步驟1，將測(cè)井縱波阻抗曲線與巖性曲線對(duì)應(yīng)，分巖性提取出各巖性的測(cè)井縱波阻抗曲線；

步驟2，計(jì)算各巖性的縱波阻抗非壓實(shí)分量；

步驟3，剔除步驟2中計(jì)算得到的各巖性縱波阻抗非壓實(shí)分量，用各巖性的縱波阻抗值減去各自巖性的縱波阻抗非壓實(shí)分量，得到剔除縱波阻抗非壓實(shí)分量后的各巖性縱波阻抗；

步驟4，分巖性分別擬合計(jì)算出剔除了縱波阻抗非壓實(shí)分量的各巖性縱波阻抗的壓實(shí)趨勢(shì)，并將其剔除掉。即利用剔除了縱波阻抗非壓實(shí)分量的各巖性縱波阻抗減去各自巖性的縱波阻抗壓實(shí)趨勢(shì)，得到去壓實(shí)趨勢(shì)后的各巖性縱波阻抗；

步驟5，將步驟2中計(jì)算得到的各巖性縱波阻抗非壓實(shí)分量回歸到步驟4中得到的去壓實(shí)趨勢(shì)后的各巖性縱波阻抗中，從而得到各巖性的去壓實(shí)縱波阻抗；

步驟6，將步驟5中得到的各巖性去壓實(shí)趨勢(shì)后的縱波阻抗分別進(jìn)行各自巖性的縱波阻抗基線賦值，得到基線賦值后的各巖性縱波阻抗；步驟7，將基線賦值后的各巖性縱波阻抗進(jìn)行合并，得到一條去壓實(shí)趨勢(shì)處理后的測(cè)井縱波阻抗曲線，并參與到后續(xù)的縱波阻抗反演中。

上述技術(shù)方案進(jìn)一步說(shuō)明：

在步驟1中，通過(guò)巖性曲線與測(cè)井的縱波阻抗曲線之間的對(duì)應(yīng)關(guān)系，將每種巖性的測(cè)井縱波阻抗單獨(dú)提取出來(lái)，從而將連續(xù)型的測(cè)井組波阻抗曲線分解成多個(gè)離散型的單巖性測(cè)井縱波阻抗曲線。

在步驟2中，各巖性的縱波阻抗值隨著埋深的增加會(huì)顯著增大，但不同巖性之間影響縱波阻抗變化的因素并不完全相同。以砂、泥巖為例，泥巖縱波阻抗的變化主要受壓實(shí)的影響，非壓實(shí)因素的影響較小，可忽略不計(jì)，而砂巖的縱波阻抗的變化除了受壓實(shí)因素的影響外，還受到含流體性、孔隙性等非壓實(shí)因素的影響，將這部分影響因素統(tǒng)稱為縱波阻抗非壓實(shí)分量，縱波阻抗非壓實(shí)分量可由巖石物理實(shí)驗(yàn)室模擬計(jì)算得出。

在步驟3中，用各巖性的縱波阻抗值減去各自巖性的縱波阻抗非壓實(shí)分量后，各巖性的縱波阻抗值隨著埋深的變化的影響因素只剩下壓實(shí)因素。剔除縱波阻抗非壓實(shí)分量的目的是為了后續(xù)能夠更準(zhǔn)確的計(jì)算出各巖性的縱波阻抗壓實(shí)趨勢(shì)。在縱波阻抗非壓實(shí)分量的剔除過(guò)程中，因?yàn)槟鄮r的縱波阻抗非壓實(shí)分量較小，可以忽略不計(jì)，但砂巖、灰?guī)r等巖性的縱波阻抗非壓實(shí)分量需要剔除。

在步驟4中，對(duì)各巖性縱波阻抗的壓實(shí)趨勢(shì)進(jìn)行擬合計(jì)算。泥巖的縱波阻抗壓實(shí)趨勢(shì)，可直接由泥巖的縱波阻抗曲線進(jìn)行擬合計(jì)算，而砂巖、灰?guī)r等巖性的縱波阻抗壓實(shí)趨勢(shì)，只有在剔除了縱波阻抗非壓實(shí)分量之后才能進(jìn)行擬合計(jì)算。在擬合計(jì)算出各巖性的縱波阻抗壓實(shí)趨勢(shì)并將其剔除后，就得到了去壓實(shí)趨勢(shì)后的各巖性縱波阻抗曲線，此時(shí)需要注意的是，泥巖僅剔除了泥巖的縱波阻抗壓實(shí)趨勢(shì)，而砂巖、灰?guī)r等巖性不僅剔除了它們的縱波阻抗壓實(shí)趨勢(shì)，在此之前還剔除了它們的縱波阻抗非壓實(shí)分量。

在步驟5中，需要將各巖性的縱波阻抗非壓實(shí)分量進(jìn)行回歸，其原因在于，泥巖的去壓實(shí)趨勢(shì)縱波阻抗僅僅是剔除了它的縱波阻抗壓實(shí)趨勢(shì)，與之相對(duì)應(yīng)，砂巖、灰?guī)r等巖性也只能是剔除了它們的縱波阻抗壓實(shí)趨勢(shì)，因而需要將步驟2中計(jì)算得到的砂巖、灰?guī)r等巖性的縱波阻抗非壓實(shí)分量進(jìn)行回歸，這樣得到的各巖性的去壓實(shí)縱波阻抗曲線相對(duì)于各自巖性的原始縱波阻抗曲線而言，均只剔除了各自巖性的縱波阻抗壓實(shí)趨勢(shì)。

在步驟5中，處理得到的各巖性的去壓實(shí)縱波阻抗具有以下特征：處理之后的各巖性縱波阻抗的值域分布范圍較為集中，并且由于消除了縱波阻抗壓實(shí)趨勢(shì)，同一巖性的縱波阻抗在淺層、中層、深層都比較接近。

在步驟6中，各巖性的測(cè)井縱波阻抗基線是從原始的測(cè)井縱波阻抗曲線中選取埋深相對(duì)較淺的層段中統(tǒng)計(jì)出來(lái)的各巖性縱波阻抗平均值，因?yàn)樵诼裆钕鄬?duì)較淺的層段中，各巖性的縱波阻抗受到壓實(shí)的影響相對(duì)較小，各巖性之間，如泥巖、砂巖、灰?guī)r等巖性在淺層埋深段中疊合區(qū)間較小，縱波阻抗平均值存在較大的差異。

在步驟6中，在選取埋深相對(duì)較淺的層段時(shí)，應(yīng)遵循以下原則：優(yōu)選目的層系的相對(duì)淺層部分，若沉積相帶變化不大，可取得較好的效果，如果沉積相帶差異較大，則可考慮取目的層系之上的相對(duì)淺層層系；在目的層系之上的相對(duì)淺層層系選擇中，優(yōu)先選擇與目的層系具有相似沉積背景的層系。遵循以上原則有利于經(jīng)過(guò)處理而最終得到的縱波阻抗曲線所計(jì)算的反射系數(shù)序列與子波褶積所得到的合成地震記錄道到原始地震之間具有最佳的匹配性。

在步驟7中，基線賦值后的各巖性縱波阻抗在深度域上都是單一巖性的離散縱波阻抗值，需要進(jìn)行合并，而合并后得到的是一條消除壓實(shí)影響的測(cè)井縱波阻抗曲線，各巖性之間的縱波阻抗具有較大的差異，因而在后續(xù)的測(cè)井約束縱波阻抗反演中能夠取得較好的效果。

本發(fā)明的適用于測(cè)井約束反演的測(cè)井縱波阻抗曲線去壓實(shí)處理方法能夠合理地放大儲(chǔ)層與圍巖之間的縱波阻抗差異，減少了儲(chǔ)層與圍巖在原始測(cè)井縱波阻抗曲線上因縱波阻抗值疊合區(qū)間較大而造成反演結(jié)果的多解性，使得最終的縱波阻抗反演結(jié)果能夠準(zhǔn)確地預(yù)測(cè)儲(chǔ)層的分布情況。

附圖說(shuō)明

圖1為本發(fā)明的具體實(shí)施例的流程圖；

圖2是實(shí)測(cè)井的原始縱波阻抗(A)曲線圖；

圖3是實(shí)測(cè)井原始縱波阻抗(A)的概率密度分布圖；

圖4是從實(shí)測(cè)井原始縱波阻抗(A)分巖性剝離出來(lái)的砂巖縱波阻抗(I)曲線和泥巖縱波阻抗(a)曲線圖；

圖5是砂巖縱波阻抗非壓實(shí)分量圖；

圖6是剔除了砂巖縱波阻抗非壓實(shí)分量后的砂巖縱波阻抗(II)曲線圖；

圖7是砂巖縱波阻抗(II)與埋深（垂深）交匯圖；

圖8是砂巖縱波阻抗(II)壓實(shí)趨勢(shì)擬合計(jì)算圖；

圖9是泥巖縱波阻抗(a)與埋深（垂深）交匯圖；

圖10是泥巖縱波阻抗(a)壓實(shí)趨勢(shì)擬合計(jì)算圖；

圖11是砂巖縱波阻抗(II)剔除砂巖縱波阻抗壓實(shí)趨勢(shì)后的砂巖縱波阻抗(III)曲線圖；

圖12是泥巖縱波阻抗(a)剔除泥巖縱波阻抗壓實(shí)趨勢(shì)后的泥巖縱波阻抗(b)曲線圖；

圖13是將砂巖縱波阻抗非壓實(shí)分量回歸到砂巖縱波阻抗(III)而得到的砂巖縱波阻抗(IV)曲線圖；

圖14是工區(qū)內(nèi)實(shí)測(cè)井統(tǒng)計(jì)的砂、泥巖縱波阻抗基線；

圖15是將砂巖縱波阻抗曲線(IV)進(jìn)行砂巖縱波阻抗基線賦值而得到的砂巖縱波阻抗(V)曲線圖；

圖16是將泥巖縱波阻抗(b)進(jìn)行泥巖縱波阻抗基線賦值而得到的泥巖縱波阻抗(c)曲線圖；

圖17是砂巖縱波阻抗(V)曲線與泥巖縱波阻抗(c)曲線合并而得到了測(cè)井縱波阻抗去壓實(shí)處理后的縱波阻抗(B)曲線圖；

圖18是去壓實(shí)處理后的縱波阻抗(B)曲線的概率密度分布圖。

具體實(shí)施方式

為使本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂，下文特舉出較佳實(shí)施例，并配合所附圖式，作詳細(xì)說(shuō)明如下。

如圖1所示，圖1為本發(fā)明的適用于測(cè)井約束反演的測(cè)井縱波阻抗去壓實(shí)處理方法的流程圖。在步驟101，選取實(shí)鉆井的測(cè)井縱波阻抗曲線。流程進(jìn)入到步驟102。

在步驟102，從實(shí)鉆井的測(cè)井縱波阻抗曲線中按照巖性的不同，分別將各自巖性的測(cè)井縱波阻抗提取出來(lái)，得到各自巖性的測(cè)井縱波阻抗曲線。流程進(jìn)入到步驟103。

圖2是實(shí)測(cè)井原始縱波阻抗(A)曲線，在顯示的井段（2220-2450米）鉆遇了砂巖和泥巖兩種巖性；圖4是實(shí)測(cè)井原始縱波阻抗(A)分巖性剝離出來(lái)的砂巖縱波阻抗(I)曲線和泥巖縱波阻抗(a)曲線圖。

在步驟103，分巖性計(jì)算各自巖性的縱波阻抗非壓實(shí)分量。流程進(jìn)入到步驟104。

圖5是砂巖的縱波阻抗非壓實(shí)分量，主要包含孔隙、流體等因素的影響，由實(shí)驗(yàn)室模擬計(jì)算。泥巖的縱波阻抗非壓實(shí)分量較小，可以忽略不計(jì)。

在步驟104，剔除各巖性縱波阻抗非壓實(shí)分量，即用步驟2中得到的各巖性縱波阻抗減去其縱波阻抗非壓實(shí)分量。流程進(jìn)入到步驟105。

在步驟105，分巖性計(jì)算各巖性的縱波阻抗壓實(shí)分量。流程進(jìn)入到步驟106。

圖7為砂巖的縱波阻抗與埋深（垂深）的交匯圖，圖8為在砂巖的縱波阻抗（II）與埋深（垂深）交匯分析的基礎(chǔ)上，進(jìn)行擬合計(jì)算得到砂巖縱波阻抗（II）的壓實(shí)趨勢(shì)；圖9為泥巖縱波阻抗與埋深（垂深）的交匯圖，圖10為在泥巖的縱波阻抗（a）與埋深（垂深）交匯分析的基礎(chǔ)上，進(jìn)行擬合計(jì)算得到泥巖縱波阻抗（a）的壓實(shí)趨勢(shì)。

在步驟106，剔除各巖性縱波阻抗壓實(shí)趨勢(shì)。流程進(jìn)入步驟107。

圖11為砂巖縱波阻抗（II）減去砂巖縱波阻抗（II）的壓實(shí)趨勢(shì)而得到的砂巖縱波阻抗（III）；圖12為泥巖縱波阻抗（a）減去泥巖縱波阻抗（a）的壓實(shí)趨勢(shì)而得到的泥巖縱波阻抗（b）。

在步驟107，分巖性進(jìn)行各巖性縱波阻抗非壓實(shí)分量回歸。砂巖需要進(jìn)行縱波阻抗非壓實(shí)分量回歸，而泥巖不需要。流程進(jìn)入步驟108。

圖13為砂巖縱波阻抗（III）加上砂巖縱波阻抗非壓實(shí)分量而得到的砂巖縱波阻抗（IV）。

在步驟108，分巖性進(jìn)行各巖性的縱波阻抗基線賦值。各巖性的縱波阻抗基線是從工區(qū)內(nèi)實(shí)鉆井的測(cè)井縱波阻抗曲線中選取埋深相對(duì)較淺的層段中統(tǒng)計(jì)出來(lái)的各巖性縱波阻抗平均值，因?yàn)樵诼裆钕鄬?duì)較淺的層段中，各巖性的縱波阻抗受到壓實(shí)的影響相對(duì)較小，各巖性之間，如泥巖、砂巖等巖性的縱波阻抗在淺層埋深段中疊合區(qū)間較小，縱波阻抗平均值存在較大的差異。在選取埋深相對(duì)較淺的層段時(shí)，應(yīng)遵循以下原則：優(yōu)選目的層系的相對(duì)淺層部分，若沉積相帶變化不大，可取得較好的效果，如果沉積相帶差異較大，則可考慮取目的層系之上的相對(duì)淺層層系；在目的層系之上的相對(duì)淺層層系選擇中，優(yōu)先選擇與目的層系具有相似沉積背景的層系。遵循以上原則有利于經(jīng)過(guò)處理而最終得到的縱波阻抗曲線所計(jì)算的反射系數(shù)序列與子波褶積所得到的合成地震記錄道與原始地震之間具有最佳的匹配性。流程進(jìn)入步驟109.

圖14為工區(qū)內(nèi)所有實(shí)鉆井的測(cè)井縱波阻抗在測(cè)井深度為1450-1850米井段所統(tǒng)計(jì)得到的砂巖和泥巖的縱波阻抗基線，這兩條縱波阻抗基線分別代表著砂、泥巖在該埋深段（1450-1850米）的原始測(cè)井縱波阻抗的平均值。圖15是將砂巖縱波阻抗曲線(IV)進(jìn)行砂巖縱波阻抗基線賦值而得到的砂巖縱波阻抗(V)曲線，即對(duì)砂巖縱波阻抗（IV）進(jìn)行整體漂移，將其平均值漂移到砂巖縱波阻抗基線上。圖16是將泥巖縱波阻抗(b)進(jìn)行泥巖縱波阻抗基線賦值而得到的泥巖縱波阻抗(c)曲線，即對(duì)泥巖縱波阻抗(b)進(jìn)行整體漂移，將其平均值漂移到泥巖縱波阻抗基線上。

在步驟109，將基線賦值后的各巖性縱波阻抗曲線進(jìn)行合并，從而得到了最終的經(jīng)過(guò)測(cè)井縱波阻抗去壓實(shí)處理后的縱波阻抗曲線。流程進(jìn)入步驟110。

圖17是砂巖縱波阻抗(V)曲線與泥巖縱波阻抗(c)曲線合并而得到的縱波阻抗(B)曲線，該曲線即原始測(cè)井縱波阻抗經(jīng)去壓實(shí)處理后而最終得到的縱波阻抗曲線；

在步驟110中，利用處理最終得到的縱波阻抗（B）曲線進(jìn)行測(cè)井約束反演。

圖18是縱波阻抗(B)曲線的概率密度分布圖，與實(shí)測(cè)井原始縱波阻抗(A)的概率密度分布圖（圖3）相比，雖然砂、泥巖之間的縱波阻抗仍存在著一定的疊合區(qū)間，但其主體已經(jīng)分離，均在著較大的區(qū)分，能夠滿足儲(chǔ)層預(yù)測(cè)的需求，可在測(cè)井約束反演中取得理想的效果。