本公開(kāi)涉及地震勘探和開(kāi)發(fā)領(lǐng)域，更具體地，涉及一種確定地層品質(zhì)因子Q的值的方法及系統(tǒng)。

背景技術(shù)：

在地震勘探和開(kāi)發(fā)領(lǐng)域，地層品質(zhì)因子Q是表征地層介質(zhì)吸收衰減作用的常用參量，Q值在烴類(lèi)檢測(cè)、油氣儲(chǔ)層預(yù)測(cè)等工作中起著非常重要的作用。

常規(guī)的Q值估算一般基于井間地震、VSP(垂直地震剖面)、疊前CMP(共中心點(diǎn))道集以及疊后資料獲取等原始數(shù)據(jù)來(lái)實(shí)現(xiàn)。用于確定地層品質(zhì)因子Q的值方法一般可分為三類(lèi)：頻率域法、時(shí)頻率域法和時(shí)間域法。實(shí)際操作中，由于采集的原始地震數(shù)據(jù)大多數(shù)是時(shí)間信號(hào)，所以使用時(shí)間域信號(hào)來(lái)估算Q值是最直接的方法，同時(shí)由于時(shí)間域法中所采用的原始數(shù)據(jù)最為可靠，所以有可能得到最為精確的Q值。但是，目前所采用的諸如振幅衰減法、上升時(shí)間法等的時(shí)間域法，其所獲得的Q值的精確度并不理想。

因此，發(fā)明人認(rèn)為，開(kāi)發(fā)一種能得到具有高精確度的Q值的方法，對(duì)于提高地震資料品質(zhì)、油氣地球物理屬性反演和儲(chǔ)層描述等是非常有必要的。

公開(kāi)于本公開(kāi)背景技術(shù)部分的信息僅僅旨在加深對(duì)本公開(kāi)的一般背景技術(shù)的理解，而不應(yīng)當(dāng)被視為承認(rèn)或以任何形式暗示該信息構(gòu)成已為本領(lǐng)域技術(shù)人員所公知的現(xiàn)有技術(shù)。

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

本公開(kāi)提供了一種能直接使用時(shí)間域信號(hào)確定具有高精確度的地層品質(zhì) 因子Q的值的方法。本公開(kāi)還提供了相應(yīng)的系統(tǒng)。

根據(jù)本公開(kāi)的一方面，提出了一種基于時(shí)間域信號(hào)確定地層品質(zhì)因子Q的值的方法，該方法可以包括：選取衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)；基于所選取的s(t)和r(t)確定下列公式中的待定系數(shù)c和n：其中，為衰減后的地震子波信號(hào)r(t)的解析信號(hào)，為衰減前的地震子波信號(hào)s(t)的解析信號(hào)，為的一階導(dǎo)數(shù)；以及根據(jù)確定地層品質(zhì)因子Q的值，其中Δt是衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)之間的走時(shí)差。

根據(jù)本公開(kāi)的另一方面，提出了一種基于時(shí)間域信號(hào)確定地層品質(zhì)因子Q的值的系統(tǒng)，該系統(tǒng)可以包括：選取單元，用于選取衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)；待定系數(shù)確定單元，用于基于所選取的s(t)和r(t)確定下列公式中的待定系數(shù)c和n：其中，為衰減后的地震子波信號(hào)r(t)的解析信號(hào)，為衰減前的地震子波信號(hào)s(t)的解析信號(hào)，為的一階導(dǎo)數(shù)；以及Q值確定單元，用于根據(jù)確定地層品質(zhì)因子Q的值，其中Δt是衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)之間的走時(shí)差。

通過(guò)應(yīng)用本公開(kāi)，可直接使用時(shí)間域信號(hào)來(lái)確定地層品質(zhì)因子Q的值，并且經(jīng)驗(yàn)證所確定的Q的值具有較高的精確度。

附圖說(shuō)明

通過(guò)結(jié)合附圖對(duì)本公開(kāi)的示例性實(shí)施例進(jìn)行更詳細(xì)的描述，本公開(kāi)的上述以及其它目的、特征和優(yōu)勢(shì)將變得更加明顯，其中，在本公開(kāi)的示例性實(shí)施例中，相同的參考標(biāo)號(hào)通常代表相同部件。

圖1示出了根據(jù)本公開(kāi)的示例性實(shí)施例的確定地層品質(zhì)因子Q的值的方法的流程圖。

圖2示出了在本公開(kāi)的一個(gè)示例性實(shí)施例中選取的衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)的波形示意圖。

圖3示出了基于圖2所示的地震子波得到的衰減前的地震子波解析信號(hào)其一階導(dǎo)數(shù)以及相對(duì)應(yīng)的衰減后的地震子波解析信號(hào)的示意圖。

圖4(A)示出了設(shè)定的測(cè)試模型的地層速度與地層深度的關(guān)系示意圖。

圖4(B)示出了與圖4(A)相同的測(cè)試模型中設(shè)定的測(cè)試模型的地層品質(zhì)因子Q的值與地層深度的關(guān)系示意圖。

圖5示出了與圖4(A)相同的測(cè)試模型中正演得到的零偏VSP下行波記錄示意圖。

圖6(A)示出了將根據(jù)本公開(kāi)的一個(gè)實(shí)施例應(yīng)用于圖4和圖5所示的測(cè)試模型，反演得到的地層品質(zhì)因子Q的值與模型設(shè)定的Q值的比較結(jié)果示意圖。

圖6(B)是基于圖6(A)所示結(jié)果的地層品質(zhì)因子Q值與模型設(shè)定的Q值的相對(duì)誤差曲線圖。

圖7示出了某一實(shí)際應(yīng)用場(chǎng)景中VSP下行波記錄示意圖。

圖8(A)示出了在該實(shí)際應(yīng)用場(chǎng)景中得到的地層速度與道號(hào)的關(guān)系示意圖。

圖8(B)示出了將本公開(kāi)的一個(gè)示例性實(shí)施例應(yīng)用于該實(shí)際應(yīng)用場(chǎng)景得到的地層品質(zhì)因子Q值與道號(hào)的關(guān)系示意圖。

圖9(A)示出了在該實(shí)際應(yīng)用場(chǎng)景中進(jìn)行反Q補(bǔ)償前的過(guò)井剖面示意圖。

圖9(B)示出了在該實(shí)際應(yīng)用場(chǎng)景中使用根據(jù)本公開(kāi)的一個(gè)示例性實(shí)施例所確定的Q值進(jìn)行反Q補(bǔ)償后的過(guò)井剖面示意圖。

具體實(shí)施方式

下面將參照附圖更詳細(xì)地描述本公開(kāi)的優(yōu)選實(shí)施例。雖然附圖中顯示了 本公開(kāi)的優(yōu)選實(shí)施例，然而應(yīng)該理解，可以以各種形式實(shí)現(xiàn)本公開(kāi)而不應(yīng)被這里闡述的實(shí)施例所限制。相反，提供這些實(shí)施例是為了使本公開(kāi)更加透徹和完整，并且能夠?qū)⒈竟_(kāi)的范圍完整地傳達(dá)給本領(lǐng)域的技術(shù)人員。

本公開(kāi)公開(kāi)了一種基于時(shí)間域信號(hào)確定地層品質(zhì)因子Q的值的方法，該方法包括：選取衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)；基于所選取的s(t)和r(t)，確定公式中的待定系數(shù)c和n，其中，為衰減后的地震子波信號(hào)r(t)的解析信號(hào)，為衰減前的地震子波信號(hào)s(t)的解析信號(hào)，為的一階導(dǎo)數(shù)；根據(jù)確定地層品質(zhì)因子Q的值，其中Δt是衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)之間的走時(shí)差。

實(shí)施例1

圖1示出了根據(jù)本公開(kāi)的示例性實(shí)施例的確定地層品質(zhì)因子Q的值的方法的流程圖。在該實(shí)施例中，該方法可以包括以下步驟：

步驟101，選取衰減前的地震子波信號(hào)s(t)和與該s(t)相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)。

在一個(gè)示例中，可以從初至波中選取衰減前的地震子波信號(hào)s(t)。初至波可較好地反應(yīng)地震波的信息，相比于非初至波所受到的干擾更小。本領(lǐng)域技術(shù)人員可以采用任意已知的方法來(lái)選取衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)。。

步驟102，基于所選取的s(t)和r(t)確定公式中的待定系數(shù)c和n，其中，為衰減后的地震子波信號(hào)r(t)的解析信號(hào)，為衰減前的地震子波信號(hào)s(t)的解析信號(hào)，為的一階導(dǎo)數(shù)。

對(duì)于給定的連續(xù)時(shí)間信號(hào)x(t)，其解析信號(hào)記為即：

其中，H[x(t)]為希爾伯特變換，是x(t)和1/πt的褶積，即：

$H\[x\left(t\right)\]=x\left(t\right)*\frac{1}{\mathrm{\π}t}=\frac{1}{\mathrm{\π}}{\∫}_{-\mathrm{\∞}}^{\mathrm{\∞}}\frac{x\left(\mathrm{\τ}\right)}{t-\mathrm{\τ}}d\mathrm{\τ}=\frac{1}{\mathrm{\π}}{\∫}_{-\mathrm{\∞}}^{\mathrm{\∞}}\frac{x(t-\mathrm{\τ})}{\mathrm{\τ}}d\mathrm{\τ}.$

對(duì)于衰減前的地震子波信號(hào)s(t)，其解析信號(hào)為：

根據(jù)Fourier變換定理，衰減前的地震子波信號(hào)s(t)的頻率域信號(hào)可以表示為S(f)，其解析信號(hào)的頻率域信號(hào)可以表示為

對(duì)于衰減后的地震子波信號(hào)r(t)，其解析信號(hào)為：

根據(jù)Fourier變換定理，衰減后的地震子波信號(hào)r(t)的頻率域信號(hào)可以表示為R(f)，其解析信號(hào)的頻率域信號(hào)可以表示為

根據(jù)Futterman衰減模型，由于上述s(t)和r(t)是互相對(duì)應(yīng)的衰減前后的地震子波信號(hào)，則有：

R(f)＝cS(f)e^-kf，

其中c是與頻率無(wú)關(guān)的衰減因子(如幾何擴(kuò)散、反射透射損失等)。

當(dāng)kf＜＜1時(shí)，將上述指數(shù)衰減項(xiàng)按Tayer級(jí)數(shù)展開(kāi)，并取一階近似，則有：

根據(jù)Fourier反變換定理，將上式從頻率域轉(zhuǎn)換到時(shí)間域，則有：

其中，c和m為衰減系數(shù)，為地層品質(zhì)因子Q，進(jìn)而可推導(dǎo)出$Q=-i\frac{\mathrm{\Δ}t}{2}\frac{c}{n}.$

在一個(gè)示例中，可以采用諸如最小二乘法的數(shù)值擬合來(lái)確定待定系數(shù)c和n。本領(lǐng)域技術(shù)人員應(yīng)理解，確定待定系數(shù)c和n并不限于以上所限定的方式，而是可以采用本領(lǐng)域技術(shù)人員已知的任何其他方法來(lái)確定。

步驟103，根據(jù)確定地層品質(zhì)因子Q的值，其中Δt是衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)之間的走時(shí)差。

可以采用本領(lǐng)域技術(shù)人員已知的任何方法得到衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)之間的走時(shí)差Δt。基于所確定的待定系數(shù)c和n，并根據(jù)可確定地層品質(zhì)因子Q的值。

進(jìn)一步地，根據(jù)本公開(kāi)，可以判斷上述步驟所確定的地層品質(zhì)因子Q的值是否滿足預(yù)定約束條件，如果不滿足，則重新選取衰減前的地震子波信號(hào)s(t)(例如，初至波中的另一個(gè)子波)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)，并基于重新選取的s(t)和r(t)來(lái)確定相應(yīng)的待定系數(shù)c和n以及地層品質(zhì)因子Q的值，直至所確定的地層品質(zhì)因子Q的值滿足所述預(yù)定約束條件。其中，預(yù)定約束可包括下面的至少一者：在預(yù)定深度范圍內(nèi)的地層速度和所確定的Q值之間的相關(guān)度γ1是否大于第一預(yù)定閾值；使用所確定的Q值對(duì)過(guò)井線進(jìn)行反Q補(bǔ)償，補(bǔ)償后的井旁道與井上合成記錄之間的相關(guān)度γ2_after與補(bǔ)償前的井旁道與井上合成記錄之間的相關(guān)度γ2_bef之間的差值是否大于第二預(yù)定閾值。

例如，在預(yù)定深度范圍內(nèi)的地層速度和所確定的Q的值之間的相關(guān)度γ1可表示為：

$\mathrm{\γ}1=\frac{\underset{H}{\∫}Q\left(h\right)V\left(h\right)dh}{\sqrt{\underset{H}{\∫}Q\left(h\right)dh}\sqrt{\underset{H}{\∫}{V}^2\left(h\right)dh}},$

其中，h為地層深度，Q(h)為地層深度為h時(shí)的地層品質(zhì)因子，V(h)為地層深度為h時(shí)的地層速度，H為積分區(qū)間(例如預(yù)定深度范圍)。例如，當(dāng)γ1大于第一預(yù)定閾值(例如0.6)時(shí)，可判斷當(dāng)前所確定的地層品質(zhì)因子Q的值滿足該預(yù)定約束條件。

例如，使用所確定的Q值對(duì)過(guò)井線進(jìn)行反Q補(bǔ)償，井旁道與井上合成記錄的相關(guān)度γ2可表示為：

$\mathrm{\γ}2=\frac{\underset{T}{\∫}x\left(t\right)y\left(t\right)dt}{\sqrt{\underset{T}{\∫}{x}^2\left(t\right)dt}\sqrt{\underset{T}{\∫}{y}^2\left(t\right)dt}},$

其中，x(t)為井旁道，y(t)為井上合成記錄，t為時(shí)間變量，T為積分區(qū)間。例如，當(dāng)補(bǔ)償后的得到的γ2_after減去補(bǔ)償前得到的γ2_bef的差值大于第二預(yù)定閾值(例如0.2)時(shí)，可判斷當(dāng)前所確定的地層品質(zhì)因子Q的值滿足該預(yù)定約束條件。

本領(lǐng)域技術(shù)人員可根據(jù)需要設(shè)定第一預(yù)定閾值和/或第二預(yù)定閾值。

應(yīng)用示例

以下給出了本公開(kāi)的實(shí)施例的若干應(yīng)用示例，用于從不同角度對(duì)本公開(kāi)的技術(shù)效果予以介紹。本領(lǐng)域技術(shù)人員應(yīng)該理解應(yīng)用示例僅為了便于理解本公開(kāi)，其中的數(shù)值及其他細(xì)節(jié)僅為示例性的，而非意在限制本公開(kāi)。

圖2和圖3示出了用于測(cè)試某模型在一定地層深度的Q值的相關(guān)波形示意圖。該模型中該地層深度的地層品質(zhì)因子Q被設(shè)定為50。圖2中實(shí)線所示波形為從初至波中選取的衰減前的地震子波信號(hào)s(t)，虛線所示波形為衰減后的地震子波信號(hào)r(t)。圖3示出了基于圖2所示的地震子波得到的衰減前的地震子波解析信號(hào)的實(shí)部(曲線2)、其一階導(dǎo)數(shù)的實(shí)部(曲線2)以及衰減后的地震子波解析信號(hào)的實(shí)部(曲線3)的示意圖。從圖3中可看出，曲線3可以近似看作是曲線2和曲線1的線性組合，即采用最小二乘法可得到c＝0.5、n＝-i×1.482×10^-4。根據(jù)本公開(kāi)可確定相應(yīng)的地層品質(zhì)因子Q的值為50.60，與模型設(shè)定值50的相對(duì)誤差為1.2％，從而可以看出根據(jù)本公開(kāi)確定的Q值具有較高的精確度。

圖4、圖5和圖6示出了針對(duì)層狀模型的測(cè)試結(jié)果。圖4(A)示出了設(shè)定的測(cè)試模型的地層速度與地層深度的關(guān)系示意圖，其中橫坐標(biāo)表示地層深度，縱坐標(biāo)表示地層速度；圖4(B)示出了該測(cè)試模型中設(shè)定的地層品質(zhì)因 子Q的值與地層深度的關(guān)系示意圖，其中橫坐標(biāo)表示地層深度，縱坐標(biāo)表示地層品質(zhì)因子Q的值。可以看出，該模型共分為六層，Q值分布在區(qū)間[50,150]的范圍內(nèi)，與速度走向趨勢(shì)一致。圖5示出了該測(cè)試模型中正演得到零偏VSP下行波記錄示意圖，檢波間距為10m，這里只考慮與頻率有關(guān)的地層吸收衰減及地層界面反射透射損失，可以看出，隨著傳播距離的增加，地震波的能量逐漸衰減，延續(xù)時(shí)間增加，分辨率降低。

圖6(A)示出了將根據(jù)本公開(kāi)的一個(gè)實(shí)施例應(yīng)用于圖4和圖5所示的測(cè)試模型，反演得到的地層品質(zhì)因子Q的值與模型設(shè)定的Q值的比較結(jié)果示意圖，其中，選取的地震數(shù)據(jù)為相鄰VSP接收道記錄，實(shí)線為模型設(shè)定的Q值，虛線為采用本公開(kāi)得到的Q值，可以看到根據(jù)本公開(kāi)反演得到的Q值能較好地反映地層品質(zhì)因子Q的走向，與設(shè)定的Q值的變化趨勢(shì)基本保持一致。圖6(B)是基于圖6(A)所示結(jié)果的地層品質(zhì)因子Q值與模型設(shè)定的Q值的相對(duì)誤差曲線圖，本公開(kāi)得到的地層品質(zhì)因子Q值與模型設(shè)定的Q值的誤差小于5％，表明本公開(kāi)在理論上具有較好的精度。

圖7、圖8和圖9展示了將本公開(kāi)中的一個(gè)實(shí)施例應(yīng)用于某一實(shí)際應(yīng)用場(chǎng)景的情況。圖7示出了實(shí)際VSP下行波記錄示意圖，該VSP下行波包含初至波、多次波和反射波等，樣品采樣間隔為1ms，檢波間距為10m，同時(shí)選取各道的初至波作為衰減前的地震子波信號(hào)。本示例中選取間隔5道的記錄來(lái)確定地層品質(zhì)因子Q的值。

在本示例中，判斷所確定的地層品質(zhì)因子是否可靠的約束條件有兩個(gè)：1，所確定的Q值曲線與速度曲線的相關(guān)度；2，使用所確定的Q值對(duì)過(guò)井線進(jìn)行反Q補(bǔ)償，補(bǔ)償前后井旁道與井上合成記錄的相關(guān)度的改善程度。如果所確定的Q值不滿足其中某一個(gè)約束條件，則調(diào)整用于截取初至波的長(zhǎng)度參數(shù)，重新計(jì)算Q值。

圖8(A)示出了在該實(shí)際應(yīng)用場(chǎng)景中得到的地層速度與道號(hào)的關(guān)系示意圖，圖8(B)示出了將本公開(kāi)的一個(gè)示例性實(shí)施例應(yīng)用于該實(shí)際應(yīng)用場(chǎng)景得 到的地層品質(zhì)因子Q值與道號(hào)的關(guān)系示意圖?？梢钥吹?，兩條曲線的趨勢(shì)基本走勢(shì)一致，且都隨深度的增加呈遞增趨勢(shì)，圖8(A)所示的Q值曲線和圖8(B)所示的地層速度曲線的相關(guān)度γ1可達(dá)到0.78。結(jié)果表明，采用本公開(kāi)得到的地層品質(zhì)因子Q的精確度較高。

圖8(B)中示出的地層品質(zhì)因子Q的值大概分布區(qū)間為[80,300]，可利用該Q值與地層速度V擬合出該工區(qū)的經(jīng)驗(yàn)關(guān)系為：Q＝13.0V^2.10。可將該關(guān)系應(yīng)用于整個(gè)工區(qū)，并用于三維資料的反Q補(bǔ)償處理。

圖9(A)示出了在該實(shí)際應(yīng)用場(chǎng)景中進(jìn)行反Q補(bǔ)償前的過(guò)井剖面示意圖(中間10道平直同相軸為井上合成記錄)，圖9(B)示出了在該實(shí)際應(yīng)用場(chǎng)景中使用根據(jù)本公開(kāi)的一個(gè)示例性實(shí)施例所確定的Q值進(jìn)行反Q補(bǔ)償后的過(guò)井剖面示意圖，目的層為1700ms附近。對(duì)比圖9(A)和圖9(B)可以看到，使用根據(jù)本公開(kāi)中的示例性實(shí)施例所確定的Q值進(jìn)行反Q補(bǔ)償后，記錄的復(fù)合波分離，同相軸增多(見(jiàn)圖中標(biāo)注1和2)，一次有效波能量得到較好的抬升(見(jiàn)標(biāo)注3)，相位得到校正，與井上合成記錄的匹配度提高(見(jiàn)標(biāo)注4)，井旁道與井上合成記錄相關(guān)度γ2從反Q補(bǔ)償前的γ2_bef＝0.4235提高到補(bǔ)償后的γ2_after＝0.7324。結(jié)果表明，采用本公開(kāi)得到的地層品質(zhì)因子Q具有較好的可信度和準(zhǔn)確度。

本領(lǐng)域技術(shù)人員應(yīng)理解，上面對(duì)本公開(kāi)的實(shí)施例的描述的目的僅為了示例性地說(shuō)明本公開(kāi)的實(shí)施例的有益效果，并不意在將本公開(kāi)的實(shí)施例限制于所給出的任何示例。

實(shí)施例2

在該實(shí)施例中，提供了一種確定地層品質(zhì)因子Q的值的系統(tǒng)，該系統(tǒng)可以包括：選取單元，用于選取衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)的單元；待定系數(shù)確定單元，基于所選取的s(t)和r(t)，用于確定公式中的待定系數(shù)c和n的單元，其中為衰減后的地震子波解析信號(hào)，為衰減前的地震子波解析信號(hào)，為的一階導(dǎo)數(shù)，c和n為待定系數(shù)；Q值確定單元，用于根據(jù)確定地層品質(zhì)因子 Q的值，其中Δt是衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)之間的走時(shí)差。

例如，該系統(tǒng)還可包括判斷單元，該判斷單元可用于判斷所確定的地層品質(zhì)因子Q的值是否滿足預(yù)定約束條件，如果不滿足，則可重新選取衰減前的地震子波信號(hào)s(t)和相對(duì)應(yīng)的衰減后的地震子波信號(hào)r(t)，并基于重新選取的s(t)和r(t)確定相應(yīng)的待定系數(shù)c和n以及地層品質(zhì)因子Q的值，直至所確定的地層品質(zhì)因子Q的值滿足所述預(yù)定約束條件。

例如，預(yù)定約束條件可包括下面的至少一者：在預(yù)定深度范圍內(nèi)的地層速度和所確定的Q的值之間的相關(guān)度γ1是否大于第一預(yù)定閾值；使用所確定的Q值對(duì)過(guò)井線進(jìn)行反Q補(bǔ)償，補(bǔ)償后的井旁道與井上合成記錄之間的相關(guān)度γ2_after與補(bǔ)償前的井旁道與井上合成記錄之間的相關(guān)度γ2_bef之間的差值是否大于第二預(yù)定閾值。

例如，在預(yù)定深度范圍內(nèi)的地層速度和所確定的Q值之間的相關(guān)度γ1可表示為：

其中h為地層深度，Q(h)為地層深度為h時(shí)的地層品質(zhì)因子，V(h)為地層深度為h時(shí)的地層速度，H為積分區(qū)間。

例如，井旁道與井上合成記錄之間的相關(guān)度γ2可表示為：

$\mathrm{\γ}2=\frac{\underset{T}{\∫}x\left(t\right)y\left(t\right)dt}{\sqrt{\underset{T}{\∫}{x}^2\left(t\right)dt}\sqrt{\underset{T}{\∫}{y}^2\left(t\right)dt}},$

其中，x(t)為井旁道，y(t)為井上合成記錄，t為時(shí)間變量，T為積分區(qū)間。

例如，可以通過(guò)諸如最小二乘法等的數(shù)值擬合來(lái)確定待定系數(shù)c和n。

例如，可從初至波中選取衰減前的地震子波信號(hào)s(t)。

本領(lǐng)域技術(shù)人員應(yīng)理解，上面對(duì)本公開(kāi)的實(shí)施例的描述的目的僅為了示例性地說(shuō)明本公開(kāi)的實(shí)施例的有益效果，并不意在將本公開(kāi)的實(shí)施例限制于 所給出的任何示例。

本公開(kāi)可以是系統(tǒng)、方法和/或計(jì)算機(jī)程序產(chǎn)品。計(jì)算機(jī)程序產(chǎn)品可以包括計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，其上載有用于使處理器實(shí)現(xiàn)本公開(kāi)的各個(gè)方面的計(jì)算機(jī)可讀程序指令。

計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)可以是可以保持和存儲(chǔ)由指令執(zhí)行設(shè)備使用的指令的有形設(shè)備。計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)例如可以是――但不限于――電存儲(chǔ)設(shè)備、磁存儲(chǔ)設(shè)備、光存儲(chǔ)設(shè)備、電磁存儲(chǔ)設(shè)備、半導(dǎo)體存儲(chǔ)設(shè)備或者上述的任意合適的組合。計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)的更具體的例子(非窮舉的列表)包括：便攜式計(jì)算機(jī)盤(pán)、硬盤(pán)、隨機(jī)存取存儲(chǔ)器(RAM)、只讀存儲(chǔ)器(ROM)、可擦式可編程只讀存儲(chǔ)器(EPROM或閃存)、靜態(tài)隨機(jī)存取存儲(chǔ)器(SRAM)、便攜式壓縮盤(pán)只讀存儲(chǔ)器(CD-ROM)、數(shù)字多功能盤(pán)(DVD)、記憶棒、軟盤(pán)、機(jī)械編碼設(shè)備、例如其上存儲(chǔ)有指令的打孔卡或凹槽內(nèi)凸起結(jié)構(gòu)、以及上述的任意合適的組合。這里所使用的計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)不被解釋為瞬時(shí)信號(hào)本身，諸如無(wú)線電波或者其他自由傳播的電磁波、通過(guò)波導(dǎo)或其他傳輸媒介傳播的電磁波(例如，通過(guò)光纖電纜的光脈沖)、或者通過(guò)電線傳輸?shù)碾娦盘?hào)。

這里所描述的計(jì)算機(jī)可讀程序指令可以從計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)下載到各個(gè)計(jì)算/處理設(shè)備，或者通過(guò)網(wǎng)絡(luò)、例如因特網(wǎng)、局域網(wǎng)、廣域網(wǎng)和/或無(wú)線網(wǎng)下載到外部計(jì)算機(jī)或外部存儲(chǔ)設(shè)備。網(wǎng)絡(luò)可以包括銅傳輸電纜、光纖傳輸、無(wú)線傳輸、路由器、防火墻、交換機(jī)、網(wǎng)關(guān)計(jì)算機(jī)和/或邊緣服務(wù)器。每個(gè)計(jì)算/處理設(shè)備中的網(wǎng)絡(luò)適配卡或者網(wǎng)絡(luò)接口從網(wǎng)絡(luò)接收計(jì)算機(jī)可讀程序指令，并轉(zhuǎn)發(fā)該計(jì)算機(jī)可讀程序指令，以供存儲(chǔ)在各個(gè)計(jì)算/處理設(shè)備中的計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中。

用于執(zhí)行本公開(kāi)操作的計(jì)算機(jī)程序指令可以是匯編指令、指令集架構(gòu)(ISA)指令、機(jī)器指令、機(jī)器相關(guān)指令、微代碼、固件指令、狀態(tài)設(shè)置數(shù)據(jù)、或者以一種或多種編程語(yǔ)言的任意組合編寫(xiě)的源代碼或目標(biāo)代碼，所述編程 語(yǔ)言包括面向?qū)ο蟮木幊陶Z(yǔ)言—諸如Smalltalk、C++等，以及常規(guī)的過(guò)程式編程語(yǔ)言—諸如“C”語(yǔ)言或類(lèi)似的編程語(yǔ)言。計(jì)算機(jī)可讀程序指令可以完全地在用戶計(jì)算機(jī)上執(zhí)行、部分地在用戶計(jì)算機(jī)上執(zhí)行、作為一個(gè)獨(dú)立的軟件包執(zhí)行、部分在用戶計(jì)算機(jī)上部分在遠(yuǎn)程計(jì)算機(jī)上執(zhí)行、或者完全在遠(yuǎn)程計(jì)算機(jī)或服務(wù)器上執(zhí)行。在涉及遠(yuǎn)程計(jì)算機(jī)的情形中，遠(yuǎn)程計(jì)算機(jī)可以通過(guò)任意種類(lèi)的網(wǎng)絡(luò)—包括局域網(wǎng)(LAN)或廣域網(wǎng)(WAN)—連接到用戶計(jì)算機(jī)，或者，可以連接到外部計(jì)算機(jī)(例如利用因特網(wǎng)服務(wù)提供商來(lái)通過(guò)因特網(wǎng)連接)。在一些實(shí)施例中，通過(guò)利用計(jì)算機(jī)可讀程序指令的狀態(tài)信息來(lái)個(gè)性化定制電子電路，例如可編程邏輯電路、現(xiàn)場(chǎng)可編程門(mén)陣列(FPGA)或可編程邏輯陣列(PLA)，該電子電路可以執(zhí)行計(jì)算機(jī)可讀程序指令，從而實(shí)現(xiàn)本公開(kāi)的各個(gè)方面。

這里參照根據(jù)本公開(kāi)實(shí)施例的方法、裝置(系統(tǒng))和計(jì)算機(jī)程序產(chǎn)品的流程圖和/或框圖描述了本公開(kāi)的各個(gè)方面。應(yīng)當(dāng)理解，流程圖和/或框圖的每個(gè)方框以及流程圖和/或框圖中各方框的組合，都可以由計(jì)算機(jī)可讀程序指令實(shí)現(xiàn)。

這些計(jì)算機(jī)可讀程序指令可以提供給通用計(jì)算機(jī)、專(zhuān)用計(jì)算機(jī)或其它可編程數(shù)據(jù)處理裝置的處理器，從而生產(chǎn)出一種機(jī)器，使得這些指令在通過(guò)計(jì)算機(jī)或其它可編程數(shù)據(jù)處理裝置的處理器執(zhí)行時(shí)，產(chǎn)生了實(shí)現(xiàn)流程圖和/或框圖中的一個(gè)或多個(gè)方框中規(guī)定的功能/動(dòng)作的裝置。也可以把這些計(jì)算機(jī)可讀程序指令存儲(chǔ)在計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中，這些指令使得計(jì)算機(jī)、可編程數(shù)據(jù)處理裝置和/或其他設(shè)備以特定方式工作，從而，存儲(chǔ)有指令的計(jì)算機(jī)可讀介質(zhì)則包括一個(gè)制造品，其包括實(shí)現(xiàn)流程圖和/或框圖中的一個(gè)或多個(gè)方框中規(guī)定的功能/動(dòng)作的各個(gè)方面的指令。

也可以把計(jì)算機(jī)可讀程序指令加載到計(jì)算機(jī)、其它可編程數(shù)據(jù)處理裝置、或其它設(shè)備上，使得在計(jì)算機(jī)、其它可編程數(shù)據(jù)處理裝置或其它設(shè)備上執(zhí)行一系列操作步驟，以產(chǎn)生計(jì)算機(jī)實(shí)現(xiàn)的過(guò)程，從而使得在計(jì)算機(jī)、其它可編 程數(shù)據(jù)處理裝置、或其它設(shè)備上執(zhí)行的指令實(shí)現(xiàn)流程圖和/或框圖中的一個(gè)或多個(gè)方框中規(guī)定的功能/動(dòng)作。

附圖中的流程圖和框圖顯示了根據(jù)本公開(kāi)的多個(gè)實(shí)施例的系統(tǒng)、方法和計(jì)算機(jī)程序產(chǎn)品的可能實(shí)現(xiàn)的體系架構(gòu)、功能和操作。在這點(diǎn)上，流程圖或框圖中的每個(gè)方框可以代表一個(gè)模塊、程序段或指令的一部分，所述模塊、程序段或指令的一部分包含一個(gè)或多個(gè)用于實(shí)現(xiàn)規(guī)定的邏輯功能的可執(zhí)行指令。在有些作為替換的實(shí)現(xiàn)中，方框中所標(biāo)注的功能也可以以不同于附圖中所標(biāo)注的順序發(fā)生。例如，兩個(gè)連續(xù)的方框?qū)嶋H上可以基本并行地執(zhí)行，它們有時(shí)也可以按相反的順序執(zhí)行，這依所涉及的功能而定。也要注意的是，框圖和/或流程圖中的每個(gè)方框、以及框圖和/或流程圖中的方框的組合，可以用執(zhí)行規(guī)定的功能或動(dòng)作的專(zhuān)用的基于硬件的系統(tǒng)來(lái)實(shí)現(xiàn)，或者可以用專(zhuān)用硬件與計(jì)算機(jī)指令的組合來(lái)實(shí)現(xiàn)。

以上已經(jīng)描述了本公開(kāi)的各實(shí)施例，上述說(shuō)明是示例性的，并非窮盡性的，并且也不限于所披露的各實(shí)施例。在不偏離所說(shuō)明的各實(shí)施例的范圍和精神的情況下，對(duì)于本技術(shù)領(lǐng)域的普通技術(shù)人員來(lái)說(shuō)許多修改和變更都是顯而易見(jiàn)的。本文中所用術(shù)語(yǔ)的選擇，旨在最好地解釋各實(shí)施例的原理、實(shí)際應(yīng)用或?qū)κ袌?chǎng)中的技術(shù)的技術(shù)改進(jìn)，或者使本技術(shù)領(lǐng)域的其它普通技術(shù)人員能理解本文披露的各實(shí)施例。