本發(fā)明涉及一種識別方法與裝置，尤其涉及一種角度彈性參數(shù)巖性識別方法及裝置，屬于石油地球勘探技術(shù)領(lǐng)域。

背景技術(shù)：

近年來我國油氣勘探的重點(diǎn)逐漸向巖性油氣藏勘探轉(zhuǎn)移。不同于構(gòu)造油氣藏，這些新型儲(chǔ)層受構(gòu)造和儲(chǔ)層非均質(zhì)性的影響，成藏條件復(fù)雜，識別難度大，投資風(fēng)險(xiǎn)高。因此，開發(fā)一套高精度的巖性識別地球物理解決方案顯得尤為迫切。

常規(guī)利用疊前地震資料進(jìn)行巖性識別的技術(shù)主要有疊前地震彈性參數(shù)反演方法。彈性參數(shù)反演基于疊前共偏移距道集或角度部分疊加道集，采用最大似然估計(jì)、約束稀疏脈沖等反演算法，反演出包含巖性信息的彈性參數(shù)體。

由于彈性參數(shù)本身物理意義的限制，其對巖性的敏感性具有一定限度，巖性區(qū)分能力并不能滿足如今精細(xì)化油氣勘探的要求且常規(guī)方法并沒有充分利用地震數(shù)據(jù)的角度信息以及多波地震中的轉(zhuǎn)換橫波地震信息。事實(shí)上，地震數(shù)據(jù)的角度信息以及轉(zhuǎn)換橫波地震信息均包含有豐富的巖性信息。因此，造成了目前的彈性參數(shù)巖性識別方法普遍存在敏感性不強(qiáng)的問題。

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

為了解決上述技術(shù)問題，本發(fā)明的目的在于提供一種角度彈性參數(shù)巖性識別方法及裝置，以克服常規(guī)巖性識別方法采用的彈性參數(shù)敏感性不足，造成巖性難以區(qū)分的問題。

為了實(shí)現(xiàn)上述技術(shù)目的，本發(fā)明提供了角度彈性參數(shù)巖性識別方法，該角度彈性參數(shù)巖性識別方法包括以下步驟：

在常規(guī)轉(zhuǎn)換橫波彈性阻抗方程的基礎(chǔ)上，得到與縱波彈性阻抗同量級的轉(zhuǎn)換橫波彈性阻抗方程；

根據(jù)所述轉(zhuǎn)換橫波彈性阻抗方程以及縱波彈性阻抗方程，得到轉(zhuǎn)換橫波彈性阻抗偽測井曲線和縱波彈性阻抗偽測井曲線；

利用轉(zhuǎn)換橫波彈性阻抗偽井曲線替代橫波阻抗測井曲線、縱波彈性阻抗偽井曲線替代縱波阻抗測井曲線，帶入常規(guī)彈性參數(shù)計(jì)算公式，得到某一角度范圍內(nèi)彈性參數(shù)偽測井曲線；

利用所述彈性參數(shù)偽測井曲線結(jié)合巖性測井曲線進(jìn)行敏感性對比，選取敏感性最強(qiáng)的角度彈性參數(shù)對應(yīng)的彈性參數(shù)為目標(biāo)角度彈性參數(shù)，對應(yīng)的角度為最佳角度；

輸入常規(guī)多波角度部分疊加地震道集；

對地震道集進(jìn)行層位進(jìn)行解釋；

采用常規(guī)約束稀疏脈沖反演算法進(jìn)行轉(zhuǎn)換橫波彈性阻抗反演，得到帶有角度信息的轉(zhuǎn)換橫波彈性阻抗數(shù)據(jù)體、縱波彈性阻抗數(shù)據(jù)體；

根據(jù)測井曲線目標(biāo)角度彈性參數(shù)及最佳角度分析結(jié)果，帶入轉(zhuǎn)換橫波彈性阻抗、縱波彈性阻抗反演數(shù)據(jù)體到相應(yīng)彈性參數(shù)計(jì)算公式中，得到角度彈性參數(shù)數(shù)據(jù)體；

根據(jù)角度彈性參數(shù)數(shù)據(jù)體巖性取值范圍，進(jìn)行巖性識別，完成所述角度彈性參數(shù)巖性的識別。

在本發(fā)明的角度彈性參數(shù)巖性識別方法中，優(yōu)選地，與縱波彈性阻抗同量級的轉(zhuǎn)換橫波彈性阻抗方程按照如下步驟獲得：

引入目的層平均橫波速度、密度，對常規(guī)轉(zhuǎn)換橫波彈性阻抗方程進(jìn)行標(biāo)準(zhǔn)化，得到標(biāo)準(zhǔn)化后的轉(zhuǎn)換橫波彈性阻抗方程：

其中：

θ表示入射角，單位為度；

V_p為縱波速度，單位為米/秒；

V_s為橫波速度，單位為米/秒；

ρ為密度，單位為千克/立方米；

V_s0為目的層段橫波速度，單位為米/秒；

ρ₀為密度的平均值，單位為千克/立方米。

在本發(fā)明的角度彈性參數(shù)巖性識別方法中，優(yōu)選地，常規(guī)多波角度部分疊加地震道集包括縱波角度部分疊加地震道集、轉(zhuǎn)換橫波角度部分疊加地震道集，其中，要求疊加道集角度與所述最佳角度相同。

在本發(fā)明的角度彈性參數(shù)巖性識別方法中，優(yōu)選地，轉(zhuǎn)換橫波彈性阻抗數(shù)據(jù)體SEI(θ)，與常規(guī)縱波彈性阻抗數(shù)據(jù)體EI(θ)，通過組合得到角度彈性參數(shù)數(shù)據(jù)體：P_ang(θ)＝f(SEI(θ),EI(θ))。

在本發(fā)明的角度彈性參數(shù)巖性識別方法中，優(yōu)選地，P_ang(θ)＝f(SEI(θ),EI(θ))中，常規(guī)彈性參數(shù)P的構(gòu)建關(guān)系表述為縱波阻抗I_p與橫波阻抗I_s的函數(shù)P＝f(I_p,I_s)。

在本發(fā)明的角度彈性參數(shù)巖性識別方法中，優(yōu)選地，最佳角度按照如下步驟確定：

利用縱波、橫波、密度測井曲線，根據(jù)轉(zhuǎn)換橫波彈性阻抗方程以及常規(guī)縱波彈性阻抗方程，分別采用不同角度計(jì)算角度彈性參數(shù)偽井曲線，結(jié)合巖性測井曲線進(jìn)行敏感性對比，選取敏感性最強(qiáng)的角度彈性參數(shù)對應(yīng)的角度為最佳角度。

在本發(fā)明的角度彈性參數(shù)巖性識別方法中，函數(shù)表達(dá)式P_ang(θ)＝f(SEI(θ),EI(θ))中的函數(shù)關(guān)系按照常規(guī)彈性參數(shù)的構(gòu)建關(guān)系來建立。常規(guī)彈性參數(shù)P的構(gòu)建關(guān)系可表述為縱波阻抗I_p與橫波阻抗I_s的函數(shù)P＝f(I_p,I_s)。角度彈性參數(shù)的求取采用同樣的函數(shù)關(guān)系，用帶有角度信息的縱波彈性阻抗EI(θ)代替縱波阻抗I_p、轉(zhuǎn)換橫波阻抗SEI(θ)替代橫波阻抗I_s。

例如，常規(guī)彈性參數(shù)中對砂巖和泥巖較為敏感的泊松比參數(shù)可表述為則角度泊松比參數(shù)為

在本發(fā)明的角度彈性參數(shù)巖性識別方法中，獲得某一角度范圍內(nèi)彈性參數(shù)偽測井曲線時(shí)，本領(lǐng)域技術(shù)人員可以根據(jù)需要確定具體的角度范圍。

本發(fā)明還提供了一種角度彈性參數(shù)巖性識別裝置，該角度彈性參數(shù)巖性識別裝置包括：

轉(zhuǎn)換模塊，用于將常規(guī)轉(zhuǎn)換橫波彈性阻抗方程轉(zhuǎn)換為與縱波彈性阻抗同量級的轉(zhuǎn)換橫波彈性阻抗方程；

測井曲線輸入模塊，用于根據(jù)本發(fā)明提出的轉(zhuǎn)換橫波彈性阻抗方程以及縱波彈性阻抗方程，得到轉(zhuǎn)換橫波彈性阻抗偽測井曲線和縱波彈性阻抗偽測井曲線；

角度彈性參數(shù)偽測井曲線計(jì)算模塊，用于利用轉(zhuǎn)換橫波彈性阻抗偽井曲線替代橫波阻抗測井曲線、縱波彈性阻抗偽井曲線替代縱波阻抗測井曲線，帶入常規(guī)彈性參數(shù)計(jì)算公式，得到某一角度范圍內(nèi)彈性參數(shù)偽測井曲線；

最佳角度模塊，用于結(jié)合巖性測井曲線進(jìn)行敏感性對比，選取敏感性最強(qiáng)的角度彈性參數(shù)對應(yīng)的彈性參數(shù)為目標(biāo)角度彈性參數(shù)，對應(yīng)的角度為最佳角度；

多波地震數(shù)據(jù)輸入模塊，用于輸入常規(guī)多波角度部分疊加地震道集；

解釋模塊，用于對地震道集進(jìn)行層位進(jìn)行解釋；

反演模塊，用于采用常規(guī)約束稀疏脈沖反演算法進(jìn)行轉(zhuǎn)換橫波彈性阻抗反演，得到帶有角度信息的轉(zhuǎn)換橫波彈性阻抗數(shù)據(jù)體、縱波彈性阻抗數(shù)據(jù)體；

角度彈性參數(shù)數(shù)據(jù)體模塊，用于根據(jù)測井曲線目標(biāo)角度彈性參數(shù)及最佳角度分析結(jié)果，將轉(zhuǎn)換橫波彈性阻抗、縱波彈性阻抗反演數(shù)據(jù)體帶入到相應(yīng)彈性參數(shù)計(jì)算公式中，得到角度彈性參數(shù)數(shù)據(jù)體；

識別模塊，用于根據(jù)角度彈性參數(shù)數(shù)據(jù)體巖性取值范圍，進(jìn)行巖性識別。

根據(jù)本發(fā)明的具體實(shí)施方式，輸入的測井曲線需要利用常規(guī)測井校正技術(shù)進(jìn)行校正，確保所有曲線能準(zhǔn)確地反映地層信息。

在本發(fā)明的角度彈性參數(shù)巖性識別裝置中，優(yōu)選地，與縱波彈性阻抗同量級的轉(zhuǎn)換橫波彈性阻抗方程按照如下步驟獲得：

引入目的層平均橫波速度、密度，對常規(guī)轉(zhuǎn)換橫波彈性阻抗方程進(jìn)行標(biāo)準(zhǔn)化，得到標(biāo)準(zhǔn)化后的轉(zhuǎn)換橫波彈性阻抗方程：

其中：

θ表示入射角，單位為度；

V_p為縱波速度，單位為米/秒；

V_s為橫波速度，單位為米/秒；

ρ為密度，單位為千克/立方米；

V_s0為目的層段橫波速度，單位為米/秒；

ρ₀為密度的平均值，單位為千克/立方米。

在本發(fā)明的角度彈性參數(shù)巖性識別裝置中，優(yōu)選地，常規(guī)多波角度部分疊加地震道集包括縱波角度部分疊加地震道集、轉(zhuǎn)換橫波角度部分疊加地震道集，其中，要求疊加道集角度與所述最佳角度相同。

在本發(fā)明的角度彈性參數(shù)巖性識別裝置中，優(yōu)選地，轉(zhuǎn)換橫波彈性阻抗數(shù)據(jù)體SEI(θ)，與常規(guī)縱波彈性阻抗數(shù)據(jù)體EI(θ)，通過組合得到角度彈性參數(shù)數(shù)據(jù)體：P_ang(θ)＝f(SEI(θ),EI(θ))。

在本發(fā)明的角度彈性參數(shù)巖性識別裝置中，優(yōu)選地，P_ang(θ)＝f(SEI(θ),EI(θ))中，常規(guī)彈性參數(shù)P的構(gòu)建關(guān)系表述為縱波阻抗I_p與橫波阻抗I_s的函數(shù)P＝f(I_p,I_s)。

在本發(fā)明的角度彈性參數(shù)巖性識別裝置中，優(yōu)選地，最佳角度按照如下步驟進(jìn)行：

利用縱波、橫波、密度測井曲線，根據(jù)轉(zhuǎn)換橫波彈性阻抗方程以及常規(guī)縱波彈性阻抗方程，分別采用不同角度計(jì)算角度彈性參數(shù)偽井曲線，結(jié)合巖性測井曲線進(jìn)行敏感性對比，選取敏感性最強(qiáng)的角度彈性參數(shù)對應(yīng)的角度為最佳角度。

本發(fā)明提供的角度彈性參數(shù)巖性識別方法和裝置，針對常規(guī)彈性參數(shù)受本身物理意義的限制且?guī)r性區(qū)分能力并不能滿足現(xiàn)如今精細(xì)化油氣勘探的要求的問題，考慮到地震數(shù)據(jù)的角度信息以及轉(zhuǎn)換橫波地震信息均包含有豐富的巖性信息。

本發(fā)明提供的角度彈性參數(shù)巖性識別方法和裝置，利用多波地震資料巖性信息豐富的特點(diǎn)，提出新的轉(zhuǎn)換橫波彈性阻抗方程，從多波地震彈性阻抗反演出發(fā)，利用縱波、轉(zhuǎn)換橫波彈性阻抗構(gòu)建新的巖性識別參數(shù)，具有較高的提取精度與易用性。

附圖說明

圖1為實(shí)施例中的角度彈性參數(shù)巖性識別裝置圖；

圖2為實(shí)施例的角度彈性參數(shù)巖性識別方法的流程圖；

圖3為實(shí)施例灰?guī)r蓋成、泥巖儲(chǔ)層的縱波速度、橫波速度、密度曲線；

圖4為得到的縱波阻抗、橫波阻抗曲線；

圖5為得到的常規(guī)10°、20°、30°轉(zhuǎn)換橫波彈性阻抗曲線；

圖6為得到的本發(fā)明10°、20°、30°轉(zhuǎn)換橫波彈性阻抗曲線；

圖7為實(shí)施例轉(zhuǎn)換橫波彈性阻抗敏感性隨入射角變化情況；

圖8為常規(guī)彈性參數(shù)與角度彈性參數(shù)敏感性對比曲線；

圖9為某工區(qū)實(shí)際測井曲線與角度彈性參數(shù)敏感性對比曲線；

圖10為某實(shí)際工區(qū)角度彈性參數(shù)反演剖面。

具體實(shí)施方式

為了對本發(fā)明的技術(shù)特征、目的和有益效果有更加清楚的理解，現(xiàn)對本發(fā)明的技術(shù)方案進(jìn)行以下詳細(xì)說明，但不能理解為對本發(fā)明的可實(shí)施范圍的限定。

實(shí)施例1

本實(shí)施例提供了一種角度彈性參數(shù)巖性識別裝置，該裝置的結(jié)構(gòu)如圖1所示，該裝置具體包括：

轉(zhuǎn)換模塊，用于將常規(guī)轉(zhuǎn)換橫波彈性阻抗方程轉(zhuǎn)換為與縱波彈性阻抗同量級的轉(zhuǎn)換橫波彈性阻抗方程；

測井曲線輸入模塊，用于根據(jù)本發(fā)明提出的轉(zhuǎn)換橫波彈性阻抗方程以及縱波彈性阻抗方程，得到轉(zhuǎn)換橫波彈性阻抗偽測井曲線和縱波彈性阻抗偽測井曲線；

角度彈性參數(shù)偽測井曲線計(jì)算模塊，用于利用轉(zhuǎn)換橫波彈性阻抗偽井曲線替代橫波阻抗測井曲線、縱波彈性阻抗偽井曲線替代縱波阻抗測井曲線，帶入常規(guī)彈性參數(shù)計(jì)算公式，得到某一角度范圍內(nèi)彈性參數(shù)偽測井曲線；

最佳角度模塊，用于結(jié)合巖性測井曲線進(jìn)行敏感性對比，選取敏感性最強(qiáng)的角度彈性參數(shù)對應(yīng)的彈性參數(shù)為目標(biāo)角度彈性參數(shù)，對應(yīng)的角度為最佳角度；

多波地震數(shù)據(jù)輸入模塊，用于輸入常規(guī)多波角度部分疊加地震道集；

解釋模塊，用于對地震道集進(jìn)行層位進(jìn)行解釋；

反演模塊，用于采用常規(guī)約束稀疏脈沖反演算法進(jìn)行轉(zhuǎn)換橫波彈性阻抗反演，得到帶有角度信息的轉(zhuǎn)換橫波彈性阻抗數(shù)據(jù)體、縱波彈性阻抗數(shù)據(jù)體；

角度彈性參數(shù)數(shù)據(jù)體模塊，用于根據(jù)測井曲線目標(biāo)角度彈性參數(shù)及最佳角度分析結(jié)果，將轉(zhuǎn)換橫波彈性阻抗、縱波彈性阻抗反演數(shù)據(jù)體帶入到相應(yīng)彈性參數(shù)計(jì)算公式中，得到角度彈性參數(shù)數(shù)據(jù)體；

識別模塊，用于根據(jù)角度彈性參數(shù)數(shù)據(jù)體巖性取值范圍，進(jìn)行巖性識別。

本實(shí)施例提供了一種角度彈性參數(shù)巖性識別方法，該識別方法的工藝流程如圖2所示，具體包括以下步驟：

步驟101，在常規(guī)轉(zhuǎn)換橫波彈性阻抗方程的基礎(chǔ)上，通過數(shù)學(xué)推導(dǎo)，得到與縱波彈性阻抗相同量級的新轉(zhuǎn)換橫波彈性阻抗方程；

常規(guī)轉(zhuǎn)換橫波彈性阻抗方程的一般形式為：

SEI(θ)＝V_s^m(θ)ρ^n(θ)

其中：

其中，θ表示入射角、V_p為縱波速度、V_s為橫波速度、ρ為密度。

該方程的不足在于轉(zhuǎn)換橫波彈性阻抗值的數(shù)量級隨著入射角度的變化非常巨大，不利于不同角度轉(zhuǎn)換橫波彈性阻抗間的比較，尤其是不利于同縱波彈性阻抗間的組合使用。圖3所示為灰?guī)r、白云巖組成的三層模型的縱波速度、橫波速度、密度曲線，計(jì)算的縱波阻抗、橫波阻抗曲線如圖4所示。根據(jù)圖3中所示模型曲線，利用常規(guī)轉(zhuǎn)換橫波彈性阻抗方程分別計(jì)算入射角度為10°、20°、30°時(shí)的轉(zhuǎn)換橫波彈性阻抗曲線，如圖5所示。由圖5可以看到，10°、20°、30°的數(shù)量級分別為10的2次方、4次方、5次方，不同入射角度轉(zhuǎn)換橫波彈性阻抗間數(shù)量級差別較大，且與橫波阻抗數(shù)量級不同，不利于組合利用。

本發(fā)明引入目的層平均橫波速度、密度，對常規(guī)轉(zhuǎn)換橫波彈性阻抗方程進(jìn)行標(biāo)準(zhǔn)化，得到標(biāo)準(zhǔn)化后的轉(zhuǎn)換橫波彈性阻抗方程：

其中V_s0、ρ₀分別表示目的層段橫波速度以及密度的平均值。

利用轉(zhuǎn)換橫波彈性阻抗方程，采用圖3所示三層模型曲線，計(jì)算轉(zhuǎn)換橫波彈性阻抗曲線如圖6所示。由圖6可以看到，不同入射角度轉(zhuǎn)換橫波彈性阻抗在與橫波阻抗在同一數(shù)量級，利于組合利用。

本實(shí)施例中，縱波彈性阻抗方程采用常規(guī)彈性阻抗反演技術(shù)所采用的方程。

步驟102，輸入縱波速度、橫波速度、密度測井曲線。輸入的測井曲線需要利用常規(guī)測井校正技術(shù)進(jìn)行校正，確保所有曲線能準(zhǔn)確地反映地層信息。

步驟103，根據(jù)步驟101所述的轉(zhuǎn)換橫波、縱波彈性阻抗方程以及步驟102輸入的測井曲線計(jì)算轉(zhuǎn)換橫波、縱波彈性阻抗偽測井曲線。

步驟104，利用轉(zhuǎn)換橫波彈性阻抗偽井曲線替代橫波阻抗測井曲線、縱波彈性阻抗偽井曲線替代縱波阻抗測井曲線，帶入常規(guī)彈性參數(shù)計(jì)算公式，得到某一角度范圍內(nèi)的彈性參數(shù)偽井曲線。結(jié)合巖性測井解釋結(jié)果進(jìn)行敏感性對比，選取敏感性最強(qiáng)的角度彈性參數(shù)為目標(biāo)角度彈性參數(shù)，對應(yīng)的角度為最佳角度θ。

本實(shí)施例利用目標(biāo)巖性與圍巖間彈性參數(shù)取值的相對差異作為巖性敏感性評價(jià)參數(shù)，選取敏感性最強(qiáng)對應(yīng)的入射角為最佳角度。如圖7所示為圖3中模型白云巖和灰?guī)r的轉(zhuǎn)換橫波彈性阻抗相對差異隨入射角的變化情況?？梢钥吹睫D(zhuǎn)換橫波彈性阻抗對巖性的敏感性并不是隨著入射角度的增加而單調(diào)遞增的，當(dāng)入射角大于一定角度時(shí)(本實(shí)施例中為35°)，敏感性反而降低，因此實(shí)際操作過程中需要根據(jù)實(shí)際測井曲線分析多波角度彈性參數(shù)敏感性隨入射角的情況，選取最佳的入射角度。

本發(fā)明實(shí)施例選取敏感性最強(qiáng)對應(yīng)的多波角度彈性參數(shù)為最佳的巖性識別參數(shù)。如圖8所示為圖3模型常規(guī)彈性參數(shù)與角度彈性參數(shù)敏感性對比情況?？梢钥吹桨自茙r層段常規(guī)縱波速/橫波速度參數(shù)、泊松比參數(shù)以及拉梅參數(shù)/剪切摸量相對灰?guī)r的差異均小于相應(yīng)的35°多波角度彈性參數(shù)，而其中35°多波角度拉梅參數(shù)/剪切模量的敏感性又最強(qiáng)。

本實(shí)施例利用不同巖性間彈性參數(shù)取值的相對差異作為巖性敏感性評價(jià)參數(shù)，選取敏感性最強(qiáng)對應(yīng)的入射角為最佳角度。

步驟105，本實(shí)施例要求輸入常規(guī)多波角度部分疊加地震道集，包括縱波角度部分疊加地震道集、轉(zhuǎn)換橫波角度部分疊加地震道集。要求疊加道集角度與步驟104選取的最佳角度相同。

步驟106，對地震道集進(jìn)行層位進(jìn)行解釋，精細(xì)的地層信息解釋結(jié)果能有效保證后續(xù)約束稀疏脈沖反演的精度。

步驟107，結(jié)合上述步驟，采用常規(guī)約束稀疏脈沖反演算法進(jìn)行彈性阻抗反演。本步驟中除了采用的轉(zhuǎn)換橫波彈性阻抗方程與常規(guī)方法不同外，其他均相同。

步驟108，根據(jù)步驟107約束稀疏脈沖反演結(jié)果，得到帶有角度信息的轉(zhuǎn)換橫波彈性阻抗以及縱波彈性阻抗數(shù)據(jù)體。

步驟109，根據(jù)步驟104中巖性敏感性分析結(jié)果，帶入轉(zhuǎn)換橫波彈性阻抗SEI(θ)、縱波彈性阻抗反演數(shù)據(jù)體到相應(yīng)彈性參數(shù)計(jì)算公式P_ang(θ)＝f(SEI(θ),EI(θ))中，得到角度彈性參數(shù)數(shù)據(jù)體P_ang(θ)。

步驟110，根據(jù)步驟109的計(jì)算結(jié)果，按照不同巖性地層的角度彈性參數(shù)取值范圍，進(jìn)行巖性識別。

圖9所示為某地區(qū)實(shí)際測井曲線、計(jì)算得到的常規(guī)及角度拉梅參數(shù)/剪切模量曲線。由圖9可以看到，在白云巖儲(chǔ)層段(圖中方框所示)角度拉梅參數(shù)/剪切模量(虛線所示)較常規(guī)縱波速度、橫波速度、密度以及拉梅參數(shù)/剪切模量(實(shí)線)幅度較大，與灰?guī)r圍巖區(qū)分明顯，白云巖儲(chǔ)層的角度拉梅參數(shù)/剪切模量值均小于1。因此可以利用角度拉梅參數(shù)/剪切模量來進(jìn)行巖性識別。

圖10為某地區(qū)角度拉梅參數(shù)/剪切模量值反演結(jié)果，可以看到反演剖面中值低于1的層段與測井巖性解釋結(jié)果(橢圓框所示為白云巖儲(chǔ)層)吻合較好，因此可以利用拉梅參數(shù)/剪切模量值來識別有利巖性。

以上實(shí)施例說明，本發(fā)明的角度彈性參數(shù)巖性識別方法和裝置具有較高的提取精度與易用性。

上述實(shí)施例闡明的系統(tǒng)、裝置、模塊或單元，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有某種功能的產(chǎn)品來實(shí)現(xiàn)。

為了描述的方便，描述以上裝置時(shí)以功能分為各種單元分別描述。當(dāng)然，在實(shí)施本申請時(shí)可以把各單元的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)。

通過以上的實(shí)施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請可借助軟件加必需的通用硬件平臺(tái)的方式來實(shí)現(xiàn)?；谶@樣的理解，本申請的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，在一個(gè)典型的配置中，計(jì)算設(shè)備包括一個(gè)或多個(gè)處理器(CPU)、輸入/輸出接口、網(wǎng)絡(luò)接口和內(nèi)存。該計(jì)算機(jī)軟件產(chǎn)品可以包括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī)，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請各個(gè)實(shí)施例或者實(shí)施例的某些部分所述的方法。該計(jì)算機(jī)軟件產(chǎn)品可以存儲(chǔ)在內(nèi)存中，內(nèi)存可能包括計(jì)算機(jī)可讀介質(zhì)中的非永久性存儲(chǔ)器，隨機(jī)存取存儲(chǔ)器(RAM)和/或非易失性內(nèi)存等形式，如只讀存儲(chǔ)器(ROM)或閃存(flash RAM)。內(nèi)存是計(jì)算機(jī)可讀介質(zhì)的示例。計(jì)算機(jī)可讀介質(zhì)包括永久性和非永久性、可移動(dòng)和非可移動(dòng)媒體可以由任何方法或技術(shù)來實(shí)現(xiàn)信息存儲(chǔ)。信息可以是計(jì)算機(jī)可讀指令、數(shù)據(jù)結(jié)構(gòu)、程序的模塊或其他數(shù)據(jù)。計(jì)算機(jī)的存儲(chǔ)介質(zhì)的例子包括，但不限于相變內(nèi)存(PRAM)、靜態(tài)隨機(jī)存取存儲(chǔ)器(SRAM)、動(dòng)態(tài)隨機(jī)存取存儲(chǔ)器(DRAM)、其他類型的隨機(jī)存取存儲(chǔ)器(RAM)、只讀存儲(chǔ)器(ROM)、電可擦除可編程只讀存儲(chǔ)器(EEPROM)、快閃記憶體或其他內(nèi)存技術(shù)、只讀光盤只讀存儲(chǔ)器(CD-ROM)、數(shù)字多功能光盤(DVD)或其他光學(xué)存儲(chǔ)、磁盒式磁帶，磁帶磁磁盤存儲(chǔ)或其他磁性存儲(chǔ)設(shè)備或任何其他非傳輸介質(zhì)，可用于存儲(chǔ)可以被計(jì)算設(shè)備訪問的信息。按照本文中的界定，計(jì)算機(jī)可讀介質(zhì)不包括短暫電腦可讀媒體(transitory media)，如調(diào)制的數(shù)據(jù)信號和載波。

本說明書中的各個(gè)實(shí)施例均采用遞進(jìn)的方式描述，各個(gè)實(shí)施例之間相同相似的部分互相參見即可，每個(gè)實(shí)施例重點(diǎn)說明的都是與其他實(shí)施例的不同之處。尤其，對于系統(tǒng)實(shí)施例而言，由于其基本相似于方法實(shí)施例，所以描述的比較簡單，相關(guān)之處參見方法實(shí)施例的部分說明即可。

本申請可用于眾多通用或?qū)Ｓ玫挠?jì)算機(jī)系統(tǒng)環(huán)境或配置中。例如：個(gè)人計(jì)算機(jī)、服務(wù)器計(jì)算機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的消費(fèi)電子設(shè)備、網(wǎng)絡(luò)PC、小型計(jì)算機(jī)、大型計(jì)算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布式計(jì)算環(huán)境等等。

本申請可以在由計(jì)算機(jī)執(zhí)行的計(jì)算機(jī)可執(zhí)行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類型的例程、程序、對象、組件、數(shù)據(jù)結(jié)構(gòu)等等。也可以在分布式計(jì)算環(huán)境中實(shí)踐本申請，在這些分布式計(jì)算環(huán)境中，由通過通信網(wǎng)絡(luò)而被連接的遠(yuǎn)程處理設(shè)備來執(zhí)行任務(wù)。在分布式計(jì)算環(huán)境中，程序模塊可以位于包括存儲(chǔ)設(shè)備在內(nèi)的本地和遠(yuǎn)程計(jì)算機(jī)存儲(chǔ)介質(zhì)中。

雖然通過實(shí)施例描繪了本申請，本領(lǐng)域普通技術(shù)人員知道，本申請有許多變形和變化而不脫離本申請的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請的精神。