本申請(qǐng)涉及油氣勘探技術(shù)領(lǐng)域，特別涉及一種沉積旋回劃分的方法和裝置。

背景技術(shù)：

在進(jìn)行油氣勘探的過程中，在對(duì)目標(biāo)地層進(jìn)行沉積旋回劃分時(shí)，往往需要先對(duì)沉積旋回進(jìn)行分析。其中，上述的沉積旋回指的是自然界周期性海平面升降的結(jié)果。具體地：在時(shí)間上具有周期性，在成因上具有相似性，在地質(zhì)上表現(xiàn)出旋回性，在地球物理上表現(xiàn)出一定響應(yīng)特征。例如，地震波場(chǎng)中地震波頻譜特征的規(guī)律性變化。因此，通過對(duì)沉積旋回的分析，可以進(jìn)一步實(shí)現(xiàn)對(duì)于地震相及沉積相的分析、層序地層研究乃至對(duì)儲(chǔ)層的預(yù)測(cè)。

目前，為了劃分沉積旋回曲線，常用的方法是根據(jù)測(cè)井資料獲得沉積旋回曲線?；蛘咴谌鄙贉y(cè)井資料的地方，根據(jù)地震資料獲取沉積旋回曲線。其中，所述根據(jù)地震資料獲取沉積旋回曲線的方法，通常是根據(jù)計(jì)算出曲線的屬性，根據(jù)曲線的趨勢(shì)，人為主觀地進(jìn)行折線化，進(jìn)而得到沉積旋回曲線。應(yīng)用這種沉積旋回曲線劃分方法，具體實(shí)施時(shí)，往往會(huì)存在誤差大、效率低，只能定性或半定量分析的技術(shù)問題。

針對(duì)上述問題，目前尚未提出有效的解決方案。

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

本申請(qǐng)實(shí)施例提供了一種沉積旋回劃分的方法和裝置，以解決現(xiàn)有沉積旋回劃分方法中劃分沉積旋回曲線時(shí)存在的誤差大、效率低的技術(shù)問題。

本申請(qǐng)實(shí)施例提供了一種沉積旋回劃分的方法，包括：

對(duì)地震資料進(jìn)行時(shí)頻分析，得到時(shí)頻數(shù)據(jù)體；

根據(jù)所述時(shí)頻數(shù)據(jù)體，得到峰值頻率屬性曲線；

對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線；

對(duì)所述低頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第一沉積旋回曲線，對(duì)所述高頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第二沉積旋回曲線；

根據(jù)所述第一沉積旋回曲線和所述第二沉積旋回曲線，進(jìn)行沉積旋回劃分。

在一個(gè)實(shí)施方式中，對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線，包括：

根據(jù)待研究的沉積旋回級(jí)別，確定中心頻率；

根據(jù)所述中心頻率和通帶帶寬，設(shè)計(jì)帶通濾波器；

通過所述帶通濾波器對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線。

在一個(gè)實(shí)施方式中，根據(jù)所述中心頻率和通帶帶寬，設(shè)計(jì)帶通濾波器，包括：

以所述中心頻率與所述通帶帶寬的一半之間的差值作為低通頻率，以所述中心頻率與所述通帶帶寬的一半的相加值作為高通頻率；

根據(jù)所述低通頻率和所述高通頻率，設(shè)計(jì)所述帶通濾波器。

在一個(gè)實(shí)施方式中，對(duì)所述低頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第一沉積旋回曲線，對(duì)所述高頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第二沉積旋回曲線，包括：

根據(jù)所述低頻屬性的帶通頻率分量曲線，得到低頻屬性最值序列；

對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的低頻屬性最值序列；

根據(jù)所述處理后的低頻屬性最值序列，通過線性插值得到所述第一沉積旋回曲線；

根據(jù)所述高頻屬性的帶通頻率分量曲線，得到高頻屬性最值序列；

對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的高頻屬性最值序列；

根據(jù)所述處理后的高頻屬性最值序列，通過線性插值得到所述第二沉積旋回曲線。

在一個(gè)實(shí)施方式中，對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的低頻屬性最值序列，包括：

根據(jù)所述待研究的沉積旋回級(jí)別，確定第一旋回閾值；

根據(jù)所述低頻屬性最值序列，確定低頻屬性最值序列的平均值；

按照低頻屬性最值序列的排列順序，分別以所述低頻屬性最值序列中兩個(gè)相鄰的最值作為低頻處理子單元，分別對(duì)所述低頻屬性最值序列中各個(gè)低頻處理子單元的兩個(gè)最值進(jìn)行作差，得到對(duì)應(yīng)的低頻處理子單元的差值的絕對(duì)值；

根據(jù)所述低頻處理子單元的差值的絕對(duì)值、所述低頻屬性最值序列的平均值和所述第一旋回閾值，對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的低頻屬性最值序列。

在一個(gè)實(shí)施方式中，根據(jù)所述低頻處理子單元的差值的絕對(duì)值、所述低頻屬性最值序列的平均值和所述第一旋回閾值，對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，包括：

在所述低頻處理子單元的差值的絕對(duì)值大于等于所述第一旋回閾值的情況下，則保留所述低頻處理子單元中的兩個(gè)最值和所述兩個(gè)最值所分別對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述低頻處理子單元的差值的絕對(duì)值小于所述第一旋回閾值，且所述低頻處理子單元中的兩個(gè)最值都小于等于所述低頻屬性最值序列的平均值的情況下，則保留所述低頻處理子單元中的兩個(gè)最值中較小的最值和所述較小的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述低頻處理子單元的差值的絕對(duì)值小于所述第一旋回閾值，且所述低頻處理子單元中的兩個(gè)最值都大于所述低頻屬性最值序列的平均值的情況下，則保留所述低頻處理子單元中的兩個(gè)最值中較大的最值和所述較大的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述低頻處理子單元的差值的絕對(duì)值小于所述第一旋回閾值，且所述低頻處理子單元中的兩個(gè)最值中一個(gè)較大最值大于所述低頻屬性最值序列的平均值，一個(gè)較小最值小于所述低頻屬性最值序列的平均值的情況下，則保留所述低頻處理子單中的兩個(gè)最值中與所述低頻屬性最值序列的平均值差值的絕對(duì)值較大的最值和對(duì)應(yīng)的時(shí)間坐標(biāo)。

在一個(gè)實(shí)施方式中，對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的高頻屬性最值序列，包括：

根據(jù)所述待研究的沉積旋回級(jí)別，確定第二旋回閾值；

根據(jù)所述高頻屬性最值序列，確定高頻屬性最值序列的平均值；

按照高頻屬性最值序列的排列順序，分別以所述高頻屬性最值序列中兩個(gè)相鄰的最值作為高頻處理子單元，分別對(duì)所述高頻屬性最值序列中各個(gè)高頻處理子單元的兩個(gè)最值進(jìn)行作差，得到對(duì)應(yīng)的高頻處理子單元的差值的絕對(duì)值；

根據(jù)所述高頻處理子單元的差值的絕對(duì)值、所述高頻屬性最值序列的平均值和所述第二旋回閾值，對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的高頻屬性最值序列。

在一個(gè)實(shí)施方式中，根據(jù)所述高頻處理子單元的差值的絕對(duì)值、所述高頻屬性最值序列的平均值和所述第二旋回閾值，對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，包括：

在所述高頻處理子單元的差值的絕對(duì)值大于等于所述第二旋回閾值的情況下，則保留所述高頻處理子單元中的兩個(gè)最值和所述兩個(gè)最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述高頻處理子單元的差值的絕對(duì)值小于所述第二旋回閾值，且所述高頻處理子單元中的兩個(gè)最值都小于等于所述高頻屬性最值序列的平均值的情況下，則保留所述高頻處理子單元中的兩個(gè)最值中較小的最值和所述較小的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述高頻處理子單元的差值的絕對(duì)值小于所述第二旋回閾值，且所述高頻處理子單元中的兩個(gè)最值都大于所述高頻屬性最值序列的平均值的情況下，則保留所述高頻處理子單元中的兩個(gè)最值中較大的最值和所述較大的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述高頻處理子單元的差值的絕對(duì)值小于所述第二旋回閾值，且所述高頻處理子單元中的兩個(gè)最值中一個(gè)較大最值大于所述高頻屬性最值序列的平均值，一個(gè)較小最值小于所述高頻屬性最值序列的平均值的情況下，則保留所述高頻處理子單中的兩個(gè)最值中與所述高頻屬性最值序列的平均值差值的絕對(duì)值較大的最值和對(duì)應(yīng)的時(shí)間坐標(biāo)。

在一個(gè)實(shí)施方式中，對(duì)所述低頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第一沉積旋回曲線，對(duì)所述高頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第二沉積旋回曲線之后，所述方法還包括：

對(duì)所述第一沉積旋回曲線進(jìn)行多次折線化處理，直到處理后的第一旋回曲線符合預(yù)設(shè)要求，將符合預(yù)設(shè)要求的第一沉積旋回曲線作為所述第一沉積旋回曲線，

對(duì)所述第二沉積旋回曲線進(jìn)行多次折線化處理，直到處理后的第二旋回曲線符合預(yù)設(shè)要求，將符合預(yù)設(shè)要求的第二沉積旋回曲線作為所述第二沉積旋回曲線。

基于相同的發(fā)明構(gòu)思，本申請(qǐng)實(shí)施例還提供了一種沉積旋回劃分的裝置，包括：

第一獲取模塊，用于對(duì)地震資料進(jìn)行時(shí)頻分析，得到時(shí)頻數(shù)據(jù)體；

第二獲取模塊，用于根據(jù)所述時(shí)頻數(shù)據(jù)體，得到峰值頻率屬性曲線；

分頻濾波模塊，用于對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線；

折線化處理模塊，用于對(duì)所述低頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第一沉積旋回曲線，對(duì)所述高頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第二沉積旋回曲線；

沉積旋回劃分模塊，用于根據(jù)所述第一沉積旋回曲線和所述第二沉積旋回曲線，進(jìn)行沉積旋回劃分。

在一個(gè)實(shí)施方式中，所述分頻濾波模塊包括：

確定濾波器單元，用于根據(jù)待研究的沉積旋回級(jí)別，確定中心頻率，并根據(jù)所述中心頻率和通帶帶寬，確定帶通濾波器；

分頻濾波單元，用于通過所述帶通濾波器對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線。

在一個(gè)實(shí)施方式中，所述折線化處理模塊包括：

低頻最值序列獲取單元，用于根據(jù)所述低頻屬性的帶通頻率分量曲線，得到低頻屬性最值序列；

低頻最值序列處理單元，用于對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的低頻屬性最值序列；

低頻折線化處理單元，用于根據(jù)所述處理后的低頻屬性最值序列，通過線性插值得到所述第一沉積旋回曲線；

高頻最值序列獲取單元，用于根據(jù)所述高頻屬性的帶通頻率分量曲線，得到高頻屬性最值序列；

高頻最值序列處理單元，用于對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的高頻屬性最值序列；

高頻折線化處理單元，用于根據(jù)所述處理后的高頻屬性最值序列，通過線性插值得到所述第二沉積旋回曲線。

在本申請(qǐng)實(shí)施例中，先通過濾波器進(jìn)行分頻處理，再對(duì)分頻后的曲線分別進(jìn)行折線化處理，得到相應(yīng)的沉積旋回曲線。通過折線化處理的方式解決了現(xiàn)有沉積旋回劃分方法中進(jìn)行沉積旋回曲線劃分時(shí)存在的誤差大、效率低，只能定性或半定量分析的技術(shù)問題，達(dá)到了充分利用地震資料，快速、準(zhǔn)確進(jìn)行沉積旋回劃分的技術(shù)效果。

附圖說明

為了更清楚地說明本申請(qǐng)實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本申請(qǐng)中記載的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1是根據(jù)本申請(qǐng)實(shí)施例的沉積旋回劃分的方法的處理流程圖；

圖2是根據(jù)本申請(qǐng)實(shí)施例的沉積旋回劃分的裝置的組成結(jié)構(gòu)圖；

圖3是應(yīng)用本申請(qǐng)實(shí)施例提供沉積旋回劃分的方法/裝置的過程示意圖；

圖4是應(yīng)用本申請(qǐng)實(shí)施例提供沉積旋回劃分的方法/裝置處理的試驗(yàn)區(qū)地震剖面圖；

圖5是應(yīng)用本申請(qǐng)實(shí)施例提供沉積旋回劃分的方法/裝置獲得的單道數(shù)據(jù)、單道時(shí)頻分析和峰值頻率屬性示意圖；

圖6是應(yīng)用本申請(qǐng)實(shí)施例提供沉積旋回劃分的方法/裝置獲得的分頻濾波后的時(shí)頻屬性圖；

圖7是應(yīng)用本申請(qǐng)實(shí)施例提供沉積旋回劃分的方法/裝置獲得的沉積旋回曲線。

具體實(shí)施方式

為了使本技術(shù)領(lǐng)域的人員更好地理解本申請(qǐng)中的技術(shù)方案，下面將結(jié)合本申請(qǐng)實(shí)施例中的附圖，對(duì)本申請(qǐng)實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本申請(qǐng)一部分實(shí)施例，而不是全部的實(shí)施例。基于本申請(qǐng)中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都應(yīng)當(dāng)屬于本申請(qǐng)保護(hù)的范圍。

考慮到現(xiàn)有的沉積旋回劃分方法，具體實(shí)施時(shí)，通常需要先劃分相應(yīng)的沉積旋回曲線。劃分沉積旋回曲線一般都是根據(jù)計(jì)算出的頻率屬性曲線，根據(jù)曲線的趨勢(shì)，人為主觀地進(jìn)行折線化。因此，導(dǎo)致現(xiàn)有沉積旋回曲線劃分方法中劃分沉積旋回曲線時(shí)往往存在誤差大、效率低，只能定性或半定量分析的技術(shù)問題。針對(duì)產(chǎn)生上述技術(shù)問題的根本原因，本申請(qǐng)考慮可以先根據(jù)地震資料獲得峰值頻率屬性曲線，再通過濾波器對(duì)該曲線進(jìn)行分頻濾波，對(duì)濾波后得到曲線分別進(jìn)行相應(yīng)的折線化處理，得到對(duì)應(yīng)的沉積旋回曲線。從而解決了現(xiàn)有的沉積旋回曲線劃分方法中劃分沉積旋回曲線時(shí)存在的誤差大、效率低，只能定性或半定量分析的技術(shù)問題，實(shí)現(xiàn)了充分利用地震資料，快速、準(zhǔn)確進(jìn)行沉積旋回劃分的技術(shù)效果，達(dá)到了定量分析沉積旋回的目的。

基于上述思考思路，本申請(qǐng)?zhí)峁┝艘环N沉積旋回劃分的方法。請(qǐng)參閱圖1。本申請(qǐng)?zhí)峁┑某练e旋回劃分的方法，具體可以包括以下步驟。

步驟101：對(duì)地震資料進(jìn)行時(shí)頻分析，得到時(shí)頻數(shù)據(jù)體。

在一個(gè)實(shí)施方式中，為了根據(jù)地震資料獲得時(shí)頻數(shù)據(jù)，具體可以是：對(duì)疊后地震數(shù)據(jù)逐道進(jìn)行廣義S變換，得到時(shí)頻數(shù)據(jù)體。需要說明的是，本實(shí)施方式中所提到的廣義S變換只是為了更好的說明本申請(qǐng)實(shí)施例中的時(shí)頻分析，所列舉的一個(gè)示意說明。除通過廣義S變換進(jìn)行時(shí)頻分析外，還可以通過短時(shí)傅里葉變化、小波變化和S變換等方法進(jìn)行時(shí)頻分析，獲取時(shí)頻數(shù)據(jù)體。具體實(shí)施時(shí)，可以根據(jù)具體情況，按照具體要求，靈活選擇合適方法對(duì)地震數(shù)據(jù)進(jìn)行時(shí)頻分析，獲取對(duì)應(yīng)的時(shí)頻數(shù)據(jù)體。對(duì)此，本申請(qǐng)不作限定。在本申請(qǐng)所涉及的實(shí)施方式中，主要通過廣義S變換對(duì)地震數(shù)據(jù)進(jìn)行時(shí)頻分析。需要說明的是，在本實(shí)施方式中，從時(shí)間域地震數(shù)據(jù)體到時(shí)頻域數(shù)據(jù)體可以是一個(gè)升維的過程。

步驟102：根據(jù)所述時(shí)頻數(shù)據(jù)體，得到峰值頻率屬性曲線。

在本實(shí)施方式中，得到瞬時(shí)譜峰值頻率屬性數(shù)據(jù)體可以等同于得到所述峰值頻率屬性曲線。即當(dāng)?shù)玫剿矔r(shí)譜峰值頻率屬性數(shù)據(jù)體時(shí)，可以通過簡單的方法，獲得相應(yīng)的峰值頻率屬性曲線。例如，可以根據(jù)瞬時(shí)譜峰值頻率屬性數(shù)據(jù)體繪制得到峰值頻率屬性曲線。也可以通過計(jì)算軟件在計(jì)算得到瞬時(shí)譜峰值頻率屬性數(shù)據(jù)體的同時(shí)，自動(dòng)生成相應(yīng)的峰值頻率屬性曲線。對(duì)此，本申請(qǐng)不作限定。

在一個(gè)實(shí)施方式中，為了根據(jù)時(shí)頻數(shù)據(jù)體獲得瞬時(shí)譜峰值頻率屬性數(shù)據(jù)體，即峰值頻率屬性曲線?？梢裕合雀鶕?jù)時(shí)頻數(shù)據(jù)體獲取瞬時(shí)譜峰值頻率屬性，其中，所述瞬時(shí)譜峰值頻率屬性可以是時(shí)間方向上單點(diǎn)振幅譜最大值所對(duì)應(yīng)的頻率。因此，可以通過對(duì)單點(diǎn)振幅譜序列進(jìn)行比較求取瞬時(shí)譜峰值頻率屬性，進(jìn)而根據(jù)求得的瞬時(shí)譜峰值頻率屬性得到瞬時(shí)譜峰值頻率屬性數(shù)據(jù)體，即峰值頻率屬性曲線。需要說明的是，在本實(shí)施方式中，相對(duì)的，從時(shí)頻域數(shù)據(jù)體到瞬時(shí)譜峰值頻率數(shù)據(jù)體是一個(gè)降維的過程。

步驟103：對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線。

在一個(gè)實(shí)施方式中，為了得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線，具體實(shí)施時(shí)，可以按照以下步驟執(zhí)行：

S1：根據(jù)待研究的沉積旋回級(jí)別，確定中心頻率。

在本實(shí)施方式中，需要說明的是，所述旋回級(jí)別可以是根據(jù)使用者的需求而定的，有時(shí)候需要研究第一旋回，則設(shè)定相應(yīng)的第一旋回級(jí)別；有時(shí)候需要研究第二旋回，則設(shè)定相應(yīng)的第二旋回級(jí)別。對(duì)此，本申請(qǐng)不作限定。此外，還需要說明的是，中心頻率可以依據(jù)旋回級(jí)別確定。例如，第一旋回選擇的中心頻率較低，第二旋回選擇中心頻率較高。在具體施工中，通常應(yīng)用地震資料來研究旋回往往最多能達(dá)到三級(jí)。比如一級(jí)旋回選擇的中心頻率為10HZ，二級(jí)旋回選擇的中心頻率為15HZ，以此類推。具體如何根據(jù)旋回級(jí)別確定中心頻率，可以根據(jù)具體實(shí)施情況和具體要求，靈活設(shè)定。對(duì)此，本申請(qǐng)不作限定。

S2：根據(jù)所述中心頻率和通帶帶寬，設(shè)計(jì)帶通濾波器。

在一個(gè)實(shí)施方式中，為了設(shè)計(jì)符合要求的帶通濾波器進(jìn)行分頻濾波，具體可以：以所述中心頻率與所述通帶帶寬的一半之間的差值作為低通頻率，以所述中心頻率與所述通帶帶寬的一半的相加值作為高通頻率；根據(jù)所述低通頻率和所述高通頻率，設(shè)計(jì)所述帶通濾波器。

在一個(gè)實(shí)施方式中，所述濾波器為巴特沃斯帶通濾波器。當(dāng)然還可以使用其他的濾波器。具體實(shí)施時(shí)，可以根據(jù)具體情況選擇合適的濾波器。對(duì)此，本申請(qǐng)不作限定。

S3：通過所述帶通濾波器對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線。

在本實(shí)施方式中，需要說明的是，所述高頻和所述低頻在此都是相對(duì)的概念。并不特指某個(gè)具體的頻率數(shù)值范圍。結(jié)合本實(shí)施方式，可以舉例說明：當(dāng)通過濾波器分頻濾波得到兩條帶通頻率分量曲線時(shí)，可以以其中中心頻率相對(duì)較高的一條帶通頻率分量曲線為高頻屬性的帶通頻率分量曲線；可以以其中中心頻率相對(duì)較低的一條帶通頻率分量曲線為低頻屬性的帶通頻率分量曲線。對(duì)于具體的頻率范圍，本申請(qǐng)?jiān)诖瞬蛔飨薅ā?/p>

需要說明的是，在本申請(qǐng)實(shí)施例中，所述第一旋回具體指的是低級(jí)旋回，所述第二旋回具體指的是高級(jí)旋回。相應(yīng)的，所述第一沉積旋回曲線具體指的是低級(jí)沉積旋回曲線，所述第二沉積旋回曲線具體指的是高級(jí)沉積旋回曲線。所述第一旋回級(jí)別具體指的是低級(jí)旋回級(jí)別，所述第二旋回級(jí)別具體指的是高級(jí)旋回級(jí)別。所述第一旋回閾值具體指的是低級(jí)旋回閾值，所述第二旋回閾值具體指的是高級(jí)旋回閾值。

步驟104：對(duì)所述低頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第一沉積旋回曲線，對(duì)所述高頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第二沉積旋回曲線。

在一個(gè)實(shí)施方式中，為了根據(jù)低頻屬性的帶通頻率分量曲線獲得第一沉積旋回曲線，為了根據(jù)高頻屬性的帶通頻率分量曲線獲得第二沉積旋回曲線，具體實(shí)施時(shí)可以按照以下步驟執(zhí)行：

S1：根據(jù)所述低頻屬性的帶通頻率分量曲線，得到低頻屬性最值序列。

在一個(gè)實(shí)施方式中，為了獲得低頻屬性最值序列，具體實(shí)施時(shí)，可以是：沿所述低頻屬性的帶通頻率分量曲線，每次從低頻屬性的帶通頻率分量曲線上第一預(yù)設(shè)個(gè)數(shù)的點(diǎn)中取出最大值和最小值；記錄取出的最大值和最小值所對(duì)應(yīng)的時(shí)間坐標(biāo)；將所取出的多組最大值和最小值按照時(shí)間坐標(biāo)表示的時(shí)間順序排列，即可得到低頻屬性最值序列。在本實(shí)施方式中，需要說明的是第一預(yù)設(shè)個(gè)數(shù)的設(shè)置可以是用1除以中心頻率得到的商乘以1000所得到的數(shù)值。例如，一條低頻屬性的帶通頻率分量曲線上按時(shí)間順序排列有如下的點(diǎn)：1、3、4、6、7、9、1、4、7、10、5、6、3、9、8、5、3、6、2、1。假設(shè)計(jì)算得到的第一預(yù)設(shè)個(gè)數(shù)的數(shù)值為4，則每次取出4個(gè)點(diǎn)，從4個(gè)點(diǎn)中選出最大值和最小值。第一次，從1、3、4、6這4個(gè)點(diǎn)中選出最大值6和最小值1，并記錄所選的6和1對(duì)應(yīng)的時(shí)間坐標(biāo)；第二次，從7、9、1、4這4個(gè)點(diǎn)中選出最大值9和最小值1，并記錄這里的9和1對(duì)應(yīng)的時(shí)間坐標(biāo)；第三次，從7、10、5、6這4個(gè)點(diǎn)中選出最大值10和最小值5，并記錄這里的10和5所對(duì)應(yīng)的時(shí)間坐標(biāo)；第四次，從3、9、8、5這4個(gè)點(diǎn)中選出最大值9和最小值3，并記錄這里的9和3所對(duì)應(yīng)的時(shí)間坐標(biāo)；第五次，從3、6、2、1這4個(gè)點(diǎn)中選出最大值6和最小值1，并記錄這里6和1所對(duì)應(yīng)的時(shí)間坐標(biāo)。根據(jù)選出的各個(gè)點(diǎn)的時(shí)間坐標(biāo)，將選出的各個(gè)最大值點(diǎn)或最小值點(diǎn)按照各個(gè)點(diǎn)所對(duì)應(yīng)的時(shí)間坐標(biāo)表示的時(shí)間順序進(jìn)行排列，得到低頻屬性最值序列：1、6、9、1、10、5、3、9、6、1。

S2：對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的低頻屬性最值序列。

在一個(gè)實(shí)施方式中，為了獲得符合要求的處理后的低頻屬性最值序列，需要對(duì)低頻屬性最值序列進(jìn)行數(shù)值處理，具體實(shí)施可以按照以下步驟執(zhí)行：

S2-1：根據(jù)所述待研究的沉積旋回級(jí)別，確定第一旋回閾值。

在一個(gè)實(shí)施方式中，所述第一旋回閾值一般相對(duì)較大，例如，可以將第一旋回閾值設(shè)置為1。當(dāng)然，這里設(shè)為1只是為了更好地說明本申請(qǐng)實(shí)施例。具體實(shí)施時(shí)，還可以設(shè)置為與1接近的其他數(shù)值，可以根據(jù)具體情況設(shè)置合適的數(shù)值作為第一旋回閾值。對(duì)此，本申請(qǐng)不作限定。

S2-2：根據(jù)所述低頻屬性最值序列，確定低頻屬性最值序列的平均值。

在一個(gè)實(shí)施方式中，例如低頻屬性最值序列為：1、6、9、1、10、5、3、9、6、1，對(duì)所述序列中各點(diǎn)求平均值為5.1，則5.1即為該低頻屬性最值序列的平均值。

S2-3：按照低頻屬性最值序列的排列順序，分別以所述低頻屬性最值序列中兩個(gè)相鄰的最值作為低頻處理子單元，分別對(duì)所述低頻屬性最值序列中各個(gè)低頻處理子單元的兩個(gè)最值進(jìn)行作差，得到對(duì)應(yīng)的低頻處理子單元的差值的絕對(duì)值。

在一個(gè)實(shí)施方式中，為了對(duì)低頻屬性最值序列進(jìn)行數(shù)值處理，需要先按照低頻屬性最值序列的排列順序，依次取兩個(gè)相鄰的最值作為一個(gè)低頻處理子單元，分別計(jì)算低頻處理子單元的差值的絕對(duì)值。例如，對(duì)于低頻屬性最值序列1、6、9、1、10、5、3、9、6、1。可以按照該低頻屬性最值序列的排列順序，依次取出1和6作為第一個(gè)低頻處理子單元，取出9和1作為第二個(gè)低頻處理子單元，取出10和5作為第三個(gè)低頻處理子單元，取出3和9作為第四個(gè)低頻處理子單元，取出6和1作為第五個(gè)低頻處理子單元。分別計(jì)算各個(gè)低頻處理子單元中兩個(gè)數(shù)值的差值的絕對(duì)值，得到第一個(gè)低頻處理子單元的差值的絕對(duì)值為5，第二個(gè)低頻處理子單元的差值的絕對(duì)值為8，第三個(gè)低頻處理子單元的差值的絕對(duì)值為5，第四個(gè)低頻處理子單元的差值的絕對(duì)值為6，第五個(gè)低頻處理子單元的差值的絕對(duì)值為5。

S2-4：根據(jù)所述低頻處理子單元的差值的絕對(duì)值、所述低頻屬性最值序列的平均值和所述第一旋回閾值，對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的低頻屬性最值序列。

在一個(gè)實(shí)施方式中，為了按要求根據(jù)所述低頻處理子單元的差值的絕對(duì)值、所述低頻屬性最值序列的平均值和所述第一旋回閾值，對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，具體可以按照以下方法執(zhí)行：

在所述低頻處理子單元的差值的絕對(duì)值大于等于所述第一旋回閾值的情況下，則保留所述低頻處理子單元中的兩個(gè)最值和所述兩個(gè)最值所分別對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述低頻處理子單元的差值的絕對(duì)值小于所述第一旋回閾值，且所述低頻處理子單元中的兩個(gè)最值都小于等于所述低頻屬性最值序列的平均值的情況下，則保留所述低頻處理子單元中的兩個(gè)最值中較小的最值和所述較小的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述低頻處理子單元的差值的絕對(duì)值小于所述第一旋回閾值，且所述低頻處理子單元中的兩個(gè)最值都大于所述低頻屬性最值序列的平均值的情況下，則保留所述低頻處理子單元中的兩個(gè)最值中較大的最值和所述較大的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述低頻處理子單元的差值的絕對(duì)值小于所述第一旋回閾值，且所述低頻處理子單元中的兩個(gè)最值中一個(gè)較大最值大于所述低頻屬性最值序列的平均值，一個(gè)較小最值小于所述低頻屬性最值序列的平均值的情況下，則保留所述低頻處理子單中的兩個(gè)最值中與所述高頻屬性最值序列的平均值差值的絕對(duì)值較大的最值和對(duì)應(yīng)的時(shí)間坐標(biāo)。

舉例而言，對(duì)于低頻屬性最值序列1、6、9、1、10、5、3、9、6、1，經(jīng)過計(jì)算得到第一個(gè)低頻處理子單元的差值的絕對(duì)值為5，第二個(gè)低頻處理子單元的差值的絕對(duì)值為8，第三個(gè)低頻處理子單元的差值的絕對(duì)值為5，第四個(gè)低頻處理子單元的差值的絕對(duì)值為6，第五個(gè)低頻處理子單元的差值的絕對(duì)值為5。此外，該低頻屬性最值序列的平均值為5.1，為了方便進(jìn)行示例說明，此處第一旋回閾值假設(shè)為6。由于第二、第四個(gè)低頻處理子單元的差值的分別為8和6均大于等于第一旋回閾值，則保留第二、第四個(gè)低頻處理子單元中的各個(gè)最值，即保留9、1和3、9，以及這4個(gè)最值對(duì)應(yīng)的時(shí)間坐標(biāo)。由于第一、第三、第五個(gè)低頻處理子單元的差值的絕對(duì)值都為5，小于第一旋回閾值6，且上述各個(gè)低頻處理子單元中的兩個(gè)最值均符合其中一個(gè)較大最值大于所述低頻屬性最值序列的平均值，一個(gè)較小最值小于所述低頻屬性最值序列的平均值的情況。例如，第一個(gè)低頻處理子單元中1小于5.1，而6大于5.1，同樣第三個(gè)低頻處理子單元中，10大于5.1，而5小于5.1，第五個(gè)低頻處理子單元中，6大于5.1，而1小于5.1。因此，對(duì)上述各個(gè)低頻處理子單元分別保留各個(gè)低頻處理子單中的兩個(gè)最值中與所述高頻屬性最值序列的平均值差值的絕對(duì)值較大的最值和對(duì)應(yīng)的時(shí)間坐標(biāo)。例如對(duì)于第一個(gè)低頻處理子單元中最值1和5.1差值的絕對(duì)值為4.9，最值6和5.1的差值的絕對(duì)值為0.9，由于4.9大于0.9，因此第一個(gè)低頻處理子單元保留最值1和對(duì)應(yīng)的時(shí)間坐標(biāo)。同理第三個(gè)低頻處理子單元保留最值10和對(duì)應(yīng)的時(shí)間坐標(biāo)，第五個(gè)低頻處理子單元保留最值1和對(duì)應(yīng)的時(shí)間坐標(biāo)。再將最后保留的最值按照所對(duì)應(yīng)的時(shí)間坐標(biāo)所表示的時(shí)間順序排列，得到處理后的低頻屬性最值序列：1、9、1、10、3、9、1。

S3：根據(jù)所述處理后的低頻屬性最值序列，通過線性插值得到所述第一沉積旋回曲線。

在一個(gè)實(shí)施方式中，為了對(duì)處理后的低頻屬性最值序列進(jìn)行線性插值，進(jìn)而得到第一旋回曲線，具體可以按照以下步驟執(zhí)行：

S3-1：按照處理后的低頻屬性最值序列的排列順序，從所述處理后的低頻屬性最值序列中選取連續(xù)的三個(gè)最值作為一個(gè)低頻折線處理組，得到多個(gè)低頻折線處理組，其中，相鄰的兩個(gè)低頻折線處理組共用一個(gè)最值；

S3-2：對(duì)所述多個(gè)低頻折線處理組中的各個(gè)低頻折線處理組分別按照以下步驟執(zhí)行，得到各個(gè)處理后的低頻折線組：

S3-2-1：確定所述低頻折線處理組中三個(gè)最值中的最大值，并將所述最大值賦值為100，所述最大值的時(shí)間坐標(biāo)保持不變；

S3-2-2：確定所述低頻折線處理組中的三個(gè)最值中的最小值，并將所述最小值賦值為0，所述最小值的時(shí)間坐標(biāo)保持不變；

S3-2-3：確定所述低頻折線處理組中的三個(gè)最值中的中間值，并去掉所述中間值；

S3-3：根據(jù)所處理后的低頻折線組，通過線性插值得到所述第一沉積旋回曲線。

舉例而言，處理后的低頻屬性最值序列為：1、4、9、1、3、4、7、9、8。按照處理后的低頻屬性最值序列的排列順序，依次從處理后的低頻屬性最值序列中取出3個(gè)最值得到一組低頻折線處理組，取出1、4、9作為第一低頻折線處理組，取出9、1、3作為第二低頻折線處理組，取出3、4、7作為第三低頻折線處理組，取出7、9、8作為第四低頻折線處理組。需要說明的是，這里相鄰的兩個(gè)低頻折線處理組共用一個(gè)最值。例如，第一低頻折線處理組和第二低頻折線處理組共用最值9，第二低頻折線處理組與第三低頻折線處理組共用最值3。參照規(guī)則，對(duì)上述四組低頻折線處理組分別進(jìn)行如下處理：將第一低頻折線處理組中的數(shù)值最大值9賦值為100，但保留該9的時(shí)間坐標(biāo)，即將原最值9改為100，而100的時(shí)間坐標(biāo)還是原最值9的時(shí)間坐標(biāo)。類似，將第一低頻折線處理組中的數(shù)值最小值賦值為1，并保留原最值1的時(shí)間坐標(biāo)作為該組數(shù)值為0的時(shí)間坐標(biāo)。此外，在第一低頻折線處理組中，最值4的數(shù)值大小既不是該組中的最大值，也不是該組中的最小值，是該組中的中間值，因此將最值4直接刪除。得到處理后的第一低頻折線處理組為0、100，其中0和100按照原來對(duì)應(yīng)最值的時(shí)間坐標(biāo)所表示的時(shí)間順序排列。參照上述對(duì)第一低頻折線處理組處理方法，可以得到處理后的第二低頻折線處理組為100、0，處理后的第三低頻折線處理組為0、100，處理后的第四低頻折線處理組為0、100。進(jìn)而可以根據(jù)處理后的第一低頻折線組、處理后的第二低頻折線組、處理后的第三低頻折線處理組和處理后的第四低頻折線處理組，對(duì)上述各個(gè)處理后的低頻折線處理組中賦值后的值根據(jù)時(shí)間坐標(biāo)位置進(jìn)行線性插值，即可得到對(duì)應(yīng)的第一沉積旋回曲線。

S4：根據(jù)所述高頻屬性的帶通頻率分量曲線，得到高頻屬性最值序列。

在一個(gè)實(shí)施方式中，為了獲得高頻屬性最值序列，具體實(shí)施時(shí)，可以是：沿所述高頻屬性的帶通頻率分量曲線，每次從高頻屬性的帶通頻率分量曲線上第一預(yù)設(shè)個(gè)數(shù)的點(diǎn)中取出最大值和最小值；記錄取出的最大值和最小值所對(duì)應(yīng)的時(shí)間坐標(biāo)；將所取出的多組最大值和最小值按照時(shí)間坐標(biāo)表示的時(shí)間順序排列，即可得到高頻屬性最值序列。在本實(shí)施方式中，需要說明的是第一預(yù)設(shè)個(gè)數(shù)的設(shè)置可以是用1除以中心頻率得到的商乘以1000所得到的數(shù)值。具體實(shí)施可以參照低頻屬性最值序列的獲取。在此，本申請(qǐng)不再贅述。

S5：對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的高頻屬性最值序列。

在一個(gè)實(shí)施方式中，為了獲得符合要求的處理后的高頻屬性最值序列，需要對(duì)高頻屬性最值序列進(jìn)行數(shù)值處理，具體實(shí)施可以按照以下步驟執(zhí)行：

S5-1：根據(jù)所述待研究的沉積旋回級(jí)別，確定第二旋回閾值。

在一個(gè)實(shí)施方式中，所述第二旋回閾值一般相對(duì)較小，例如，可以將第二旋回閾值設(shè)置為0.5。當(dāng)然，這里設(shè)為0.5只是為了更好地說明本申請(qǐng)實(shí)施例。具體實(shí)施時(shí)，還可以設(shè)置為與0.5接近的其他數(shù)值，可以根據(jù)具體情況設(shè)置合適的數(shù)值作為第二旋回閾值。對(duì)此，本申請(qǐng)不作限定。

S5-2：根據(jù)所述高頻屬性最值序列，確定高頻屬性最值序列的平均值。

在本實(shí)施方式中，確定高頻屬性最值序列的平均值的方法和確定低頻屬性最值序列的平均值的方法類似，具體可以參照確定低頻屬性最值序列的平均值的方法。在此，本申請(qǐng)不再贅述。

S5-3：按照高頻屬性最值序列的排列順序，分別以所述高頻屬性最值序列中兩個(gè)相鄰的最值作為高頻處理子單元，分別對(duì)所述高頻屬性最值序列中各個(gè)高頻處理子單元的兩個(gè)最值進(jìn)行作差，得到對(duì)應(yīng)的高頻處理子單元的差值的絕對(duì)值。

在一個(gè)實(shí)施方式中，為了對(duì)高頻屬性最值序列進(jìn)行數(shù)值處理，需要先按照高頻屬性最值序列的排列順序，依次取兩個(gè)相鄰的最值作為一個(gè)高頻處理子單元，分別計(jì)算高頻處理子單元的差值的絕對(duì)值。具體實(shí)施與計(jì)算低頻處理子單元的差值的絕對(duì)值類似，可以參照計(jì)算低頻處理子單元的差值的絕對(duì)值的過程。在此，本申請(qǐng)不再贅述。

S5-4：根據(jù)所述高頻處理子單元的差值的絕對(duì)值、所述高頻屬性最值序列的平均值和所述第二旋回閾值，對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的高頻屬性最值序列。

在一個(gè)實(shí)施方式中，為了按要求根據(jù)所述高頻處理子單元的差值的絕對(duì)值、所述高頻屬性最值序列的平均值和所述第二旋回閾值，對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，具體可以按照以下方法執(zhí)行：

在所述高頻處理子單元的差值的絕對(duì)值大于等于所述第二旋回閾值的情況下，則保留所述高頻處理子單元中的兩個(gè)最值和所述兩個(gè)最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述高頻處理子單元的差值的絕對(duì)值小于所述第二旋回閾值，且所述高頻處理子單元中的兩個(gè)最值都小于等于所述高頻屬性最值序列的平均值的情況下，則保留所述高頻處理子單元中的兩個(gè)最值中較小的最值和所述較小的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述高頻處理子單元的差值的絕對(duì)值小于所述第二旋回閾值，且所述高頻處理子單元中的兩個(gè)最值都大于所述高頻屬性最值序列的平均值的情況下，則保留所述高頻處理子單元中的兩個(gè)最值中較大的最值和所述較大的最值所對(duì)應(yīng)的時(shí)間坐標(biāo)；

在所述高頻處理子單元的差值的絕對(duì)值小于所述第二旋回閾值，且所述高頻處理子單元中的兩個(gè)最值中一個(gè)較大最值大于所述高頻屬性最值序列的平均值，一個(gè)較小最值小于所述高頻屬性最值序列的平均值的情況下，則保留所述低頻處理子單中的兩個(gè)最值中與所述高頻屬性最值序列的平均值差值的絕對(duì)值較大的最值和對(duì)應(yīng)的時(shí)間坐標(biāo)。

具體實(shí)施時(shí)，對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理與對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理類似，可以參照對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理的過程。在此，本申請(qǐng)不再贅述。

S6：根據(jù)所述處理后的高頻屬性最值序列，通過線性插值得到所述第二沉積旋回曲線。

在一個(gè)實(shí)施方式中，為了對(duì)處理后的高頻屬性最值序列進(jìn)行線性插值，進(jìn)而得到第二旋回曲線，具體可以按照以下步驟執(zhí)行：

S6-1：按照處理后的高頻屬性最值序列的排列順序，從所述處理后的高頻屬性最值序列中選取連續(xù)的三個(gè)最值作為一個(gè)高頻折線處理組，得到多個(gè)高頻折線處理組，其中，相鄰的兩個(gè)高頻折線處理組共用一個(gè)最值；

S6-2：對(duì)所述多個(gè)高頻折線處理組中的各個(gè)高頻折線處理組分別按照以下步驟執(zhí)行，得到各個(gè)處理后的高頻折線組：

S6-2-1：確定所述高頻折線處理組中三個(gè)最值中的最大值，并將所述最大值賦值為100，所述最大值的時(shí)間坐標(biāo)保持不變；

S6-2-2：確定所述高頻折線處理組中的三個(gè)最值中的最小值，并將所述最小值賦值為0，所述最小值的時(shí)間坐標(biāo)保持不變；

S6-2-3：確定所述高頻折線處理組中的三個(gè)最值中的中間值，并去掉所述中間值；

S6-3：根據(jù)所處理后的高頻折線組，通過線性插值得到所述第二沉積旋回曲線。

具體實(shí)施時(shí)，根據(jù)所述處理后的高頻屬性最值序列，通過線性插值得到所述第二沉積旋回曲線的具體處理過程與根據(jù)低頻屬性最值序列求得第一沉積旋回曲線類似，可以參照求解第一沉積旋回曲線的過程。在此，本申請(qǐng)不再贅述。

步驟105：根據(jù)所述第一沉積旋回曲線和所述第二沉積旋回曲線，進(jìn)行沉積旋回劃分。

在一個(gè)實(shí)施方式中，可以根據(jù)沉積旋回劃分的結(jié)果，進(jìn)行地質(zhì)勘探。其中，所述根據(jù)沉積旋回劃分的結(jié)果，進(jìn)行地質(zhì)勘探，可以包括：根據(jù)沉積旋回劃分的結(jié)果，對(duì)地震相及沉積相進(jìn)行分析，分析層序地層或者預(yù)測(cè)儲(chǔ)層。具體地，可以根據(jù)沉積旋回劃分的結(jié)果，進(jìn)行井位部署。需要說明的是，上述所列舉幾種利用沉積旋回劃分結(jié)果進(jìn)行地質(zhì)勘探的方式只是為了更好地說明本申請(qǐng)實(shí)施例，具體實(shí)施時(shí)，也可以根據(jù)具體情況或者施工要求，利用所述沉積旋回劃分的結(jié)果進(jìn)行其他相應(yīng)的地質(zhì)勘探工作。對(duì)此，本申請(qǐng)不作限定。

在本申請(qǐng)實(shí)施例中，由于先通過濾波器對(duì)峰值頻率屬性曲線進(jìn)行分頻處理，再對(duì)分頻后的曲線分別進(jìn)行折線化處理，得到相應(yīng)的沉積旋回曲線，解決了現(xiàn)有沉積旋回劃分方法中沉積旋回曲線劃分時(shí)存在的誤差大、效率低，只能定性或半定量分析的技術(shù)問題，實(shí)現(xiàn)了充分利用地震資料，快速、準(zhǔn)確劃分沉積旋回的技術(shù)效果，達(dá)到了定量分析沉積旋回的目的，進(jìn)而可以為井位部署提供部署依據(jù)。

在一個(gè)實(shí)施方式中，為了獲得更加準(zhǔn)確的第一沉積旋回曲線和第二沉積旋回曲線，在根據(jù)所述處理后的低頻屬性最值序列，通過線性插值得到所述第一沉積旋回曲線，根據(jù)所述處理后的高頻屬性最值序列，通過線性插值得到所述第二沉積旋回曲線之后，所述方法還可以包括：重復(fù)線性差值處理，對(duì)所述第一沉積旋回曲線進(jìn)行多次折線化處理，直到處理后的第一旋回曲線符合預(yù)設(shè)要求，將符合預(yù)設(shè)要求的第一沉積旋回曲線作為所述第一沉積旋回曲線，對(duì)所述第二沉積旋回曲線進(jìn)行多次折線化處理，直到處理后的第二旋回曲線符合預(yù)設(shè)要求，將符合預(yù)設(shè)要求的第二沉積旋回曲線作為所述第二沉積旋回曲線。進(jìn)而可以提高最后得到的第一沉積旋回曲線和第二沉積旋回曲線的準(zhǔn)確度，進(jìn)一步減少誤差。

在一個(gè)實(shí)施方式，可以進(jìn)一步擴(kuò)展，可以根據(jù)具體要求或者施工情況，對(duì)3級(jí)旋回級(jí)別進(jìn)行分析。由于是分析對(duì)象有3個(gè)旋回級(jí)別，對(duì)應(yīng)的在設(shè)計(jì)濾波器進(jìn)行分頻時(shí)需要考慮到3個(gè)旋回級(jí)別，設(shè)計(jì)濾波器時(shí)，分別設(shè)定一級(jí)旋回中心頻率和二級(jí)旋回中心頻率，根據(jù)一級(jí)旋回中心頻率和二級(jí)旋回中心頻率設(shè)計(jì)所相應(yīng)的分頻濾波器，通過該濾波器得到3條帶通濾波分量曲線，再對(duì)這3條帶通濾波分量曲線分別進(jìn)行相應(yīng)的數(shù)值處理和折線化處理，具體可以參照上述根據(jù)低頻帶通濾波分量曲線獲得第一沉積旋回曲線的過程進(jìn)行處理，進(jìn)而可以得到對(duì)應(yīng)的3條沉積旋回曲線，實(shí)現(xiàn)對(duì)3級(jí)旋回級(jí)別的分析處理。當(dāng)然對(duì)于1級(jí)旋回級(jí)別的分析處理，得到1條沉積旋回曲線，可以參照上述處理2級(jí)旋回級(jí)別，得到2條沉積旋回曲線和處理3級(jí)旋回級(jí)別，得到3條沉積旋回曲線的情況。在此，本申請(qǐng)不再贅述。

基于同一發(fā)明構(gòu)思，本發(fā)明實(shí)施例中還提供了一種沉積旋回劃分的裝置，如下面的實(shí)施例所述。由于裝置解決問題的原理與沉積旋回劃分的方法相似，因此沉積旋回劃分的裝置的實(shí)施可以參見沉積旋回劃分的方法的實(shí)施，重復(fù)之處不再贅述。以下所使用的，術(shù)語“單元”或者“模塊”可以實(shí)現(xiàn)預(yù)定功能的軟件和/或硬件的組合。盡管以下實(shí)施例所描述的裝置較佳地以軟件來實(shí)現(xiàn)，但是硬件，或者軟件和硬件的組合的實(shí)現(xiàn)也是可能并被構(gòu)想的。請(qǐng)參閱圖2，是本發(fā)明實(shí)施例的沉積旋回劃分的裝置的一種組成結(jié)構(gòu)圖，該裝置可以包括：第一獲取模塊201、第二獲取模塊202、分頻濾波模塊203、折線化處理模塊204和沉積旋回劃分模塊205，下面對(duì)該結(jié)構(gòu)進(jìn)行具體說明。

第一獲取模塊201，用于對(duì)地震資料進(jìn)行時(shí)頻分析，得到時(shí)頻數(shù)據(jù)體；

第二獲取模塊202，用于根據(jù)所述時(shí)頻數(shù)據(jù)體，得到峰值頻率屬性曲線；

分頻濾波模塊203，用于對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線；

折線化處理模塊204，用于對(duì)所述低頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第一沉積旋回曲線，對(duì)所述高頻屬性的帶通頻率分量曲線進(jìn)行折線化處理，得到第二沉積旋回曲線；

沉積旋回劃分模塊205，用于根據(jù)所述第一沉積旋回曲線和所述第二沉積旋回曲線，進(jìn)行沉積旋回劃分。

在一個(gè)實(shí)施方式中，為了按照旋回級(jí)別，進(jìn)行分頻濾波，得到對(duì)應(yīng)的多條帶通濾波分量曲線，所述分頻濾波模塊204可以包括：

確定濾波器單元，用于根據(jù)待研究的沉積旋回級(jí)別，確定中心頻率，并根據(jù)所述中心頻率和通帶帶寬，確定帶通濾波器；

分頻濾波單元，用于通過所述帶通濾波器對(duì)所述峰值頻率屬性曲線進(jìn)行分頻濾波，得到低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線。

在一個(gè)實(shí)施方式中，為了對(duì)帶通頻率分量曲線分別進(jìn)行相應(yīng)的數(shù)值處理和折線化處理，得到對(duì)應(yīng)的沉積旋回曲線，所述折線化處理模塊205可以包括：

低頻最值序列獲取單元，用于根據(jù)所述低頻屬性的帶通頻率分量曲線，得到低頻屬性最值序列；

低頻最值序列處理單元，用于對(duì)所述低頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的低頻屬性最值序列；

低頻折線化處理單元，用于根據(jù)所述處理后的低頻屬性最值序列，通過線性插值得到所述第一沉積旋回曲線；

高頻最值序列獲取單元，用于根據(jù)所述高頻屬性的帶通頻率分量曲線，得到高頻屬性最值序列；

高頻最值序列處理單元，用于對(duì)所述高頻屬性最值序列進(jìn)行數(shù)值處理，得到處理后的高頻屬性最值序列；

高頻折線化處理單元，用于根據(jù)所述處理后的高頻屬性最值序列，通過線性插值得到所述第二沉積旋回曲線。

在一個(gè)實(shí)施方式中，為了獲得更加準(zhǔn)確的第一沉積旋回曲線和第二沉積旋回曲線，在通過折線化處理模塊205獲得所述第一沉積旋回曲線和所述第二沉積旋回曲線之后，還可以通過折線化處理模塊205對(duì)所述第一沉積旋回曲線和所述第二沉積旋回曲線分別進(jìn)行多次的折線化處理，直到處理后的第二旋回曲線和第一沉積旋回曲線符合預(yù)設(shè)要求，將符合預(yù)設(shè)要求的第二沉積旋回曲線和第一沉積旋回曲線分別作為所述第二沉積旋回曲線和第一沉積旋回曲線。

在一個(gè)實(shí)施方式，可以進(jìn)一步擴(kuò)展，可以根據(jù)具體要求或者施工情況，對(duì)3級(jí)旋回級(jí)別進(jìn)行分析。由于是分析對(duì)象有3個(gè)旋回級(jí)別，分頻濾波模塊204中的確定濾波器單元在執(zhí)行具體操作時(shí)，考慮到3個(gè)旋回級(jí)別，確定相對(duì)應(yīng)的濾波器：先分別設(shè)定一級(jí)旋回中心頻率和二級(jí)旋回中心頻率，再根據(jù)一級(jí)旋回中心頻率和二級(jí)旋回中心頻率設(shè)計(jì)所相應(yīng)的分頻濾波器，進(jìn)而分頻濾波模塊204中的分頻濾波單元可以通過該濾波器得到3條對(duì)應(yīng)的帶通濾波分量曲線，再通過折線化處理模塊205這3條帶通濾波分量曲線分別進(jìn)行相應(yīng)的數(shù)值處理和折線化處理，具體可以參照上述說明。在此，本申請(qǐng)不再贅述。當(dāng)然對(duì)于1級(jí)旋回級(jí)別的分析處理，得到1條沉積旋回曲線，可以參照上述處理2級(jí)旋回級(jí)別，得到2條沉積旋回曲線和處理3級(jí)旋回級(jí)別，得到3條沉積旋回曲線的情況。在此，本申請(qǐng)同樣也不再贅述。

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

需要說明的是，上述實(shí)施方式闡明的系統(tǒng)、裝置、模塊或單元，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有某種功能的產(chǎn)品來實(shí)現(xiàn)。為了描述的方便，在本說明書中，描述以上裝置時(shí)以功能分為各種單元分別描述。當(dāng)然，在實(shí)施本申請(qǐng)時(shí)可以把各單元的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)。

此外，在本說明書中，諸如第一和第二這樣的形容詞僅可以用于將一個(gè)元素或動(dòng)作與另一元素或動(dòng)作進(jìn)行區(qū)分，而不必要求或暗示任何實(shí)際的這種關(guān)系或順序。在環(huán)境允許的情況下，參照元素或部件或步驟(等)不應(yīng)解釋為局限于僅元素、部件、或步驟中的一個(gè)，而可以是元素、部件、或步驟中的一個(gè)或多個(gè)等。

從以上的描述中，可以看出，本申請(qǐng)實(shí)施例提供的沉積旋回劃分的方法和裝置由于先通過濾波器對(duì)峰值頻率屬性曲線進(jìn)行分頻處理，再對(duì)分頻后的曲線分別進(jìn)行折線化處理，得到相應(yīng)的沉積旋回曲線，因此解決了現(xiàn)有方法中沉積旋回曲線劃分時(shí)存在的誤差大、效率低，只能定性或半定量分析的技術(shù)問題；由于本方法和裝置是根據(jù)地震資料獲取峰值頻率屬性曲線，再根據(jù)峰值頻率屬性曲線進(jìn)行后續(xù)的處理，因此充分地利用了地震資料，解決了在測(cè)井資料缺少的地區(qū)無法根據(jù)測(cè)井資料沉積旋回曲線劃分的技術(shù)問題，達(dá)到了擴(kuò)展適用范圍的技術(shù)效果，進(jìn)而可以實(shí)現(xiàn)沉積旋回劃分，為井位部署提供依據(jù)；又由于本方法和裝置通過對(duì)分頻過濾的曲線分別進(jìn)行相對(duì)應(yīng)的數(shù)值處理，再進(jìn)行相應(yīng)的折線化處理，提高了沉積旋回曲線劃分的準(zhǔn)確性，進(jìn)而實(shí)現(xiàn)了快速、準(zhǔn)確確定沉積相旋回曲線，進(jìn)而進(jìn)行沉積旋回劃分的技術(shù)效果。

在一個(gè)具體的實(shí)施例中，應(yīng)用本申請(qǐng)?zhí)峁┑某练e旋回劃分的方法/裝置對(duì)某地進(jìn)行具體的沉積旋回的研究。具體可以按照以下方法執(zhí)行。其中，具體的施工流程可以參閱圖3。

請(qǐng)先參閱圖4，待處理的是某地實(shí)際時(shí)間域純波地震資料，其中，采樣率為4ms。

步驟1：對(duì)時(shí)間域的地震數(shù)據(jù)體進(jìn)行廣義S變換時(shí)頻分析，得到基于廣義S變換時(shí)頻分析的時(shí)頻域數(shù)據(jù)體，即時(shí)頻數(shù)據(jù)體。其中原始地震剖面可以參閱圖4所示。選擇過該地震剖面中某一點(diǎn)的單道地震數(shù)據(jù)參閱圖5中的左圖所示，對(duì)該單道數(shù)據(jù)進(jìn)行廣義S變換時(shí)頻分析結(jié)果可以參閱圖5中的中間圖所示。

步驟2：對(duì)圖5中的中間圖所示的單道時(shí)頻析結(jié)果，求取時(shí)間方向上單點(diǎn)振幅譜最大值所對(duì)應(yīng)的頻率，即瞬時(shí)譜峰值頻率屬性。具體結(jié)果可以參閱圖5中的右圖所示。

步驟3：根據(jù)所需要研究的沉積旋回級(jí)別(Ta～Tb),確定所要求取的沉積旋回曲線的中心頻率(F＝1/T)。由此設(shè)計(jì)帶通濾波器，使得低通頻率和高通頻率分別為：fa＝F-M和fb＝F+M(M一般取5～10)；對(duì)峰值頻率屬性曲線進(jìn)行帶通濾波，得到相應(yīng)頻率成分的帶通頻率分量曲線：低頻屬性的帶通頻率分量曲線CF1和高頻屬性的帶通頻率分量曲線CF2。具體可以參閱圖6。

步驟4：分別對(duì)低頻屬性的帶通頻率分量曲線和高頻屬性的帶通頻率分量曲線進(jìn)行折線化，具體可以根據(jù)所設(shè)計(jì)的帶通濾波器中心頻率F，對(duì)帶通濾波后的帶通頻率分量曲線CF上每間隔T個(gè)樣點(diǎn)上分別求取最大值A(chǔ)max和最小值A(chǔ)min并記錄其對(duì)應(yīng)的時(shí)間坐標(biāo)C_Amax及C_Amin。按照時(shí)間坐標(biāo)順序排列所求取的最值序列。求取最值序列的平均值M_mean。根據(jù)待研究的沉積旋回級(jí)別(Ta～Tb)確定閾值，一般地，第一旋回閾值設(shè)置較大如1.0左右，第二旋回閾值設(shè)置較小如0.5左右。對(duì)最值序列中相鄰兩個(gè)最值進(jìn)行做差，保留插值大于等于閾值的最值序列。對(duì)于差值小于閾值的最值做如下處理：如果相鄰兩最值小于平均值M_mean，則保留較小的最值及對(duì)應(yīng)的時(shí)間坐標(biāo)；反之，如果相鄰兩最值大于平均值M_mean，則保留較大的最值及對(duì)應(yīng)的時(shí)間坐標(biāo)；此外，如果相鄰兩最值一個(gè)大于平均值M_mean，一個(gè)小于平均值M_mean，則保留與平均值M_mean相差較大的最值及對(duì)應(yīng)的時(shí)間坐標(biāo)。對(duì)經(jīng)過上述處理后的最值序列，每選擇連續(xù)3個(gè)作為一組，其中下一組要包含上一組的最后一個(gè)值，然后確定這三個(gè)值中的最大值和最小值并且將最大值賦值為100，最小值賦值為0。對(duì)賦值后的值根據(jù)坐標(biāo)位置進(jìn)行線性插值，便得到了所要求取的沉積旋回曲線：第一沉積旋回曲線和第二沉積旋回曲線。具體結(jié)果可以參閱圖7所示。

通過應(yīng)用本申請(qǐng)?zhí)岢龅某练e旋回劃分的方法/裝置根據(jù)試驗(yàn)區(qū)內(nèi)的一種無井地區(qū)疊后地震資料進(jìn)行沉積旋回曲線的確定。由于是針對(duì)無井地區(qū)，所以通過工區(qū)內(nèi)待鉆井點(diǎn)的測(cè)井解釋結(jié)果來驗(yàn)證其在該井點(diǎn)周圍沉積旋回劃分。通過比較，驗(yàn)證了本申請(qǐng)?zhí)峁┑某练e旋回劃分的方法和裝置解決了現(xiàn)有方法中存在的沉積旋回曲線劃分誤差大、效率低，只能定性或半定量分析的技術(shù)問題，實(shí)現(xiàn)了充分利用地震資料，快速、準(zhǔn)確確定沉積相旋回曲線的技術(shù)效果，達(dá)到了定量分析沉積旋回的目的。

盡管本申請(qǐng)內(nèi)容中提到不同的沉積旋回劃分的方法或裝置，但是，本申請(qǐng)并不局限于必須是行業(yè)標(biāo)準(zhǔn)或?qū)嵤├枋龅那闆r等，某些行業(yè)標(biāo)準(zhǔn)或者使用自定義方式或?qū)嵤├枋龅膶?shí)施基礎(chǔ)上略加修改后的實(shí)施方案也可以實(shí)現(xiàn)上述實(shí)施例相同、等同或相近、或變形后可預(yù)料的實(shí)施效果。應(yīng)用這些修改或變形后的數(shù)據(jù)獲取、處理、輸出、判斷方式等的實(shí)施例，仍然可以屬于本申請(qǐng)的可選實(shí)施方案范圍之內(nèi)。

雖然本申請(qǐng)?zhí)峁┝巳鐚?shí)施例或流程圖所述的方法操作步驟，但基于常規(guī)或者無創(chuàng)造性的手段可以包括更多或者更少的操作步驟。實(shí)施例中列舉的步驟順序僅僅為眾多步驟執(zhí)行順序中的一種方式，不代表唯一的執(zhí)行順序。在實(shí)際中的裝置或客戶端產(chǎn)品執(zhí)行時(shí)，可以按照實(shí)施例或者附圖所示的方法順序執(zhí)行或者并行執(zhí)行(例如并行處理器或者多線程處理的環(huán)境，甚至為分布式數(shù)據(jù)處理環(huán)境)。術(shù)語“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過程、方法、產(chǎn)品或者設(shè)備不僅包括那些要素，而且還包括沒有明確列出的其他要素，或者是還包括為這種過程、方法、產(chǎn)品或者設(shè)備所固有的要素。在沒有更多限制的情況下，并不排除在包括所述要素的過程、方法、產(chǎn)品或者設(shè)備中還存在另外的相同或等同要素。

上述實(shí)施例闡明的裝置或模塊等，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有某種功能的產(chǎn)品來實(shí)現(xiàn)。為了描述的方便，描述以上裝置時(shí)以功能分為各種模塊分別描述。當(dāng)然，在實(shí)施本申請(qǐng)時(shí)可以把各模塊的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)，也可以將實(shí)現(xiàn)同一功能的模塊由多個(gè)子模塊的組合實(shí)現(xiàn)等。以上所描述的裝置實(shí)施例僅僅是示意性的，例如，所述模塊的劃分，僅僅為一種邏輯功能劃分，實(shí)際實(shí)現(xiàn)時(shí)可以有另外的劃分方式，例如多個(gè)模塊或組件可以結(jié)合或者可以集成到另一個(gè)系統(tǒng)，或一些特征可以忽略，或不執(zhí)行。

本領(lǐng)域技術(shù)人員也知道，除了以純計(jì)算機(jī)可讀程序代碼方式實(shí)現(xiàn)控制器以外，完全可以通過將方法步驟進(jìn)行邏輯編程來使得控制器以邏輯門、開關(guān)、專用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來實(shí)現(xiàn)相同功能。因此這種控制器可以被認(rèn)為是一種硬件部件，而對(duì)其內(nèi)部包括的用于實(shí)現(xiàn)各種功能的裝置也可以視為硬件部件內(nèi)的結(jié)構(gòu)?；蛘呱踔粒梢詫⒂糜趯?shí)現(xiàn)各種功能的裝置視為既可以是實(shí)現(xiàn)方法的軟件模塊又可以是硬件部件內(nèi)的結(jié)構(gòu)。

本申請(qǐng)可以在由計(jì)算機(jī)執(zhí)行的計(jì)算機(jī)可執(zhí)行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類型的例程、程序、對(duì)象、組件、數(shù)據(jù)結(jié)構(gòu)、類等等。也可以在分布式計(jì)算環(huán)境中實(shí)踐本申請(qǐng)，在這些分布式計(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ù)人員可以清楚地了解到本申請(qǐng)可借助軟件加必需的通用硬件平臺(tái)的方式來實(shí)現(xiàn)?；谶@樣的理解，本申請(qǐng)的技術(shù)方案本質(zhì)上或者說對(duì)現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，該計(jì)算機(jī)軟件產(chǎn)品可以存儲(chǔ)在存儲(chǔ)介質(zhì)中，如ROM/RAM、磁碟、光盤等，包括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī)，移動(dòng)終端，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請(qǐng)各個(gè)實(shí)施例或者實(shí)施例的某些部分所述的方法。

本說明書中的各個(gè)實(shí)施例采用遞進(jìn)的方式描述，各個(gè)實(shí)施例之間相同或相似的部分互相參見即可，每個(gè)實(shí)施例重點(diǎn)說明的都是與其他實(shí)施例的不同之處。本申請(qǐng)可用于眾多通用或?qū)Ｓ玫挠?jì)算機(jī)系統(tǒng)環(huán)境或配置中。例如：個(gè)人計(jì)算機(jī)、服務(wù)器計(jì)算機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的電子設(shè)備、網(wǎng)絡(luò)PC、小型計(jì)算機(jī)、大型計(jì)算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布式計(jì)算環(huán)境等等。

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