本發(fā)明涉及勘探技術(shù)，具體的講是一種地層層序劃分方法及裝置。

背景技術(shù)：

層序地層學(xué)在油氣田勘探開發(fā)階段中已得到廣泛的應(yīng)用，如何對層序進行準確劃分是層序地層學(xué)研究中最基礎(chǔ)、最重要的一個環(huán)節(jié)。準確的層序劃分是進行沉積旋回分析，沉積環(huán)境研究的基礎(chǔ)，對地層的充填模式及沉積相的研究具有重要意義。層序是一套相對整一的、成因上有聯(lián)系的、其頂和底面以不整合面或者與這些不整合面可以對比的整合面為界的地層，是反映自然界周期性海平面升降的結(jié)果。層序反映著沉積旋回的變換，長周期的沉積旋回在地震剖面上有更明顯的響應(yīng)，而短周期的沉積旋回則在測井曲線上有更敏感的特征。因此，要想準確對層序進行劃分，必須結(jié)合地震和測井資料進行綜合分析。

隨著計算機技術(shù)的發(fā)展，層序劃分的方法也越來越多，現(xiàn)有技術(shù)中以小波變換、頻譜屬性趨勢分析和砂泥比曲線分析為主。

小波變換是一種時頻分析方法。利用小波變換劃分層序，由多個不同周期(尺度)沉積旋回疊加的測井曲線，通過小波變換，被分解成各自周期獨立的沉積旋回，以尺度的形式展示出來。研究表明，由多個不同周期沉積旋回疊加的測井曲線，通過考察小波時頻能量圖局部能量團的變化和多種伸縮尺度的周期性震蕩特征，可分析地層的旋回性并與各級層序界面建立對應(yīng)關(guān)系。例如，同一套地層在小波能譜圖上表現(xiàn)為穩(wěn)定的能量團，能量團間的斷點以及小波曲線的變化點對應(yīng)地層中的沉積間歇面或剝蝕面，表征沉積環(huán)境的改變，因此可作為層序界面，這就是小波時頻分析劃分地層的依據(jù)。

頻譜屬性趨勢分析的主要研究方法是以最大熵方法為基礎(chǔ)對測井資料進行最大熵頻譜屬性分析，其實質(zhì)就是根據(jù)研究目的不同把包含有多種頻譜成分的測井曲線利用滑動窗口技術(shù)分解成單一頻譜成分的資料，從而得到頻譜變化屬性(PEFA)曲線，利用PEFA技術(shù)分析測井資料，可以快速識別地層間斷即層序界面。對PEFA曲線進行積分運算得到更有價值的INPEFA曲線，應(yīng)用其拐點識別層序界面和洪泛面。

運用砂泥比曲線劃分層序，根據(jù)測井曲線采用計算機自動劃分高分辨率層序時，首先通過自然伽馬測井曲線計算出泥砂比、砂泥比曲線，然后綜合利用這些曲線求出地層短期基準面變化曲線，最后根據(jù)短期基準面變化曲線，綜合其他信息劃分出不同級別的層序。

現(xiàn)有技術(shù)的層序劃分方法都是對測井資料進行分析，而在測井曲線上只能對短周期高頻層序進行很好的識別，而長周期之間例如二級和三級層序，兩者在測井曲線上的特征差異不大，采用現(xiàn)有技術(shù)難以通過測井曲線區(qū)分是二級還是三級層序。

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

為解決現(xiàn)有技術(shù)中無法準確對長周期低頻層序進行準確識別的問題，本發(fā)明提出一種利用不同尺度識別不同級別層序的方法來進行層序劃分，本發(fā)明實施例提供了一種地層層序劃分方法，包括：

根據(jù)地震剖面的地震反射特征識別二級層序界面、三級層序界面；

將所述二級層序界面、三級層序界面分別定位到所述地震剖面中的一測井的對應(yīng)深度；

在所述二級層序界面的對應(yīng)深度點處對所述測井的測井曲線進行一級重構(gòu)，生成一級重構(gòu)曲線；

在所述三級層序界面的對應(yīng)深度點處對所述一級重構(gòu)曲線進行二級重構(gòu)，生成二級重構(gòu)曲線；

按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成地層層序劃分結(jié)果。

本發(fā)明實施例中，將二級層序界面、三級層序界面分別定位到所述地震剖面中的一測井的對應(yīng)深度包括：

根據(jù)所述地震剖面中的測井的合成記錄的時深關(guān)系分別將所述二級層序界面時間點t₁，t₂，…，t_n對應(yīng)到所述測井的測井曲線上的深度點h₁，h₂，…，h_n；

根據(jù)所述合成記錄的時深關(guān)系將所述三級層序界面的時間點t_n1，t_n2，…，t_nm對應(yīng)到測井曲線上的深度點h_n1，h_n2，…，h_nm；

其中，t_nm為第n個二級層序內(nèi)第m個三級層序界面的時間點，n表示二級層序的次序，m表示三級層序界面的次序，h_nm為第n個二級層序內(nèi)第m個三級層序界面的深度點，n表示二級層序的次序，m表示三級層序界面的次序，m、n∈N₊，N₊為正整數(shù)。

本發(fā)明實施例中，所述的在所述二級層序界面的對應(yīng)深度點處對所述測井的測井曲線進行一級重構(gòu)，生成的一級重構(gòu)曲線為：

S_B(h)＝S_A(h)+S_Amax/8+(S_Amax/8)*(-1)ⁿ⁺¹，h_n<h<h_(n+1)，n∈N₊

其中，S_B(h)為重構(gòu)后的一級重構(gòu)曲線B在各深度點處的測井值，S_A(h)為測井曲線A在各深度點處的原測井值，S_Amax為測井曲線的全局最大值，即原始測井曲線A整個深度段的最大測井值，h_n為二級層序界面對應(yīng)的不同深度點，n表示二級層序的次序，N₊為正整數(shù)。

本發(fā)明實施例中，所述的在所述三級層序界面的對應(yīng)深度點處對所述一級重構(gòu)曲線進行二級重構(gòu)，生成的二級重構(gòu)曲線為：

S_C(h)＝S_B(h)-S_Amax(n)/16+(S_Amax(n)/16)*(-1)^m，h_nm<h<h_n(m+1)，m、n∈N₊

其中，S_C(h)為二級重構(gòu)曲線C在各深度點處的測井值，S_B(h)為重構(gòu)后的一級重構(gòu)曲線B在各深度點處的測井值，S_Amax(n)為曲線A在第n個二級層序內(nèi)部的最大測井值，h_nm為一個二級層序n+1內(nèi)部第m個三級層序界面對應(yīng)的深度點，h_n(m+1)為某一個二級層序n+1內(nèi)部第m+1個三級層序界面對應(yīng)的深度點，n表示二級層序的次序，m表示三級層序界面的次序，N₊為正整數(shù)。

本發(fā)明實施例中，所述的按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成地層層序劃分結(jié)果包括：

按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成小波能量譜系圖和小波系數(shù)曲線圖；

根據(jù)所述小波能量譜系圖和小波系數(shù)曲線圖的間斷點對所述二級重構(gòu)曲線進行層序劃分。

本發(fā)明實施例中，所述按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成地層層序劃分結(jié)果進一步包括：

對所述二級重構(gòu)權(quán)限進行6000尺度的小波分析，生成6000尺度的小波能量譜系圖，根據(jù)該小波能量譜系圖中的間斷點劃分二級層序；

對所述二級重構(gòu)權(quán)限進行1400尺度的小波分析，生成1400尺度的小波能量譜系圖，根據(jù)該小波能量譜系圖中的間斷點劃分三級層序；

對所述二級重構(gòu)權(quán)限進行700尺度的小波分析，生成700尺度的小波能量譜系圖，根據(jù)該小波能量譜系圖中的間斷點劃分四級層序；

對所述二級重構(gòu)權(quán)限進行150尺度的小波分析，生成小波系數(shù)曲線圖，根據(jù)小波系數(shù)曲線圖中的間斷點劃分五級層序。

同時，本發(fā)明還提供一種地層層序劃分裝置，包括：

識別模塊，用于根據(jù)地震剖面的地震反射特征識別二級層序界面、三級層序界面；

定位模塊，用于將所述二級層序界面、三級層序界面分別定位到所述地震剖面中的一測井的對應(yīng)深度；

一級重構(gòu)模塊，用于在所述二級層序界面的對應(yīng)深度點處對所述測井的測井曲線進行一級重構(gòu)，生成一級重構(gòu)曲線；

二級重構(gòu)模塊，在所述三級層序界面的對應(yīng)深度點處對所述一級重構(gòu)曲線進行二級重構(gòu)，生成二級重構(gòu)曲線；

地層劃分模塊，用于按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成地層層序劃分結(jié)果。

本發(fā)明實施例中，定位模塊包括：

二級層序定位單元，根據(jù)所述地震剖面中的測井的合成記錄的時深關(guān)系分別將所述二級層序界面時間點t₁，t₂，…，t_n對應(yīng)到所述測井的測井曲線上的深度點h₁，h₂，…，h_n；

三級層序定位單元，根據(jù)所述合成記錄的時深關(guān)系將所述三級層序界面的時間點t_n1，t_n2，…，t_nm對應(yīng)到測井曲線上的深度點h_n1，h_n2，…，h_nm；

其中，t_nm為第n個二級層序內(nèi)第m個三級層序界面的時間點，n表示二級層序的次序，m表示三級層序界面的次序，h_nm為第n個二級層序內(nèi)第m個三級層序界面的深度點，n表示二級層序的次序，m表示三級層序界面的次序，m、n∈N₊，N₊為正整數(shù)。

本發(fā)明實施例中，所述的地層劃分模塊包括：

小波分析單元，用于按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成小波能量譜系圖和小波系數(shù)曲線圖；

劃分單元，用于根據(jù)所述小波能量譜系圖和小波系數(shù)曲線圖的間斷點對所述二級重構(gòu)曲線進行層序劃分。

本發(fā)明實施例中，小波分析單元進行小波分析包括：

對所述二級重構(gòu)權(quán)限進行6000尺度的小波分析，生成6000尺度的小波能量譜系圖，根據(jù)該小波能量譜系圖中的間斷點劃分二級層序；

對所述二級重構(gòu)權(quán)限進行1400尺度的小波分析，生成1400尺度的小波能量譜系圖，根據(jù)該小波能量譜系圖中的間斷點劃分三級層序；

對所述二級重構(gòu)權(quán)限進行700尺度的小波分析，生成700尺度的小波能量譜系圖，根據(jù)該小波能量譜系圖中的間斷點劃分四級層序；

對所述二級重構(gòu)權(quán)限進行150尺度的小波分析，生成小波系數(shù)曲線圖，根據(jù)小波系數(shù)曲線圖中的間斷點劃分五級層序。

本發(fā)明對長期低頻層序利用地震識別，再通過合成記錄定位到測井上；這樣通過地震對長期層序的準確識別保證了層序曲線中長期低頻層序劃分的結(jié)果準確可靠。并且，本發(fā)明依據(jù)地震對測井曲線經(jīng)過兩次重構(gòu)得到二級重構(gòu)曲線，最后直接對二級重構(gòu)曲線進行不同尺度的小波變換分析，得到不同級別層序的劃分方案，其實現(xiàn)了層序劃分過程中真正的井震結(jié)合。

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

附圖說明

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

圖1為本發(fā)明公開的一種地層層序劃分方法的流程圖；

圖2為本發(fā)明實施例中井的地震剖面圖；

圖3為本發(fā)明實施例中二級重構(gòu)曲線及小波變換能量譜系圖；

圖4本發(fā)明實施例中二級重構(gòu)曲線及小波變換能量譜系圖；

圖5本發(fā)明實施例中二級重構(gòu)曲線及小波變換能量譜系圖；

圖6為二級重構(gòu)曲線及其對應(yīng)的150尺度小波系數(shù)曲線圖；

圖7為本發(fā)明公開的一種地層層序劃分裝置的框圖；

圖8為本發(fā)明實施例中的框圖；

圖9為本發(fā)明實施例中的框圖。

具體實施方式

下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例。基于本發(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護的范圍。

如圖1所示，本發(fā)明提供一種地層層序劃分方法，包括：

步驟S101，根據(jù)地震剖面的地震反射特征識別二級層序界面、三級層序界面；

步驟S102，將所述二級層序界面、三級層序界面分別定位到所述地震剖面中的一測井的對應(yīng)深度；

步驟S103，在所述二級層序界面的對應(yīng)深度點處對所述測井的測井曲線進行一級重構(gòu)，生成一級重構(gòu)曲線；

步驟S104，在所述三級層序界面的對應(yīng)深度點處對所述一級重構(gòu)曲線進行二級重構(gòu)，生成二級重構(gòu)曲線；

步驟S105，按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成地層層序劃分結(jié)果。

本發(fā)明通過重構(gòu)測井曲線，將地震剖面上識別的二、三級層序界面信息附加到測井上，使得其在后續(xù)小波變換中可以被準確識別，解決了原始測井曲線本身難以準確識別長期層序的難題。對重構(gòu)后的測井曲線進行小波分析，通過選取不同尺度的不同小波進行小波變換，依據(jù)得到的小波時頻能量圖和小波系數(shù)曲線圖來分析各個級別的層序旋回。

由于1級層序尺度很大，很容易就在地震上識別；而6級層序尺度很小，是基于巖心進行識別，二者幾乎都不需要以測井曲線為依據(jù)，故本發(fā)明不考慮這兩種層序的識別，本發(fā)明主要針對現(xiàn)有技術(shù)基于測井曲線進行2-5級層序劃分時，無法準確對長周期低頻層序(二級、三級)進行準確識別的問題，提出一種利用不同尺度識別不同級別層序的方法來進行層序劃分。借助地震對長周期旋回的準確識別能力，將其識別結(jié)果轉(zhuǎn)換到測井曲線上，再通過對測井曲線進行分尺度的小波變換分析，逐級準確識別出不同級次的層序旋回。本發(fā)明具體實施時包括如下步驟：

步驟1：依據(jù)地震剖面區(qū)域上明顯的上超下削的地震反射特征，識別出n個二級層序界面，劃分為n+1個二級層序；同時根據(jù)合成記錄的時深關(guān)系將地震上二級層序界面的時間點t₁，t₂，…，t_n對應(yīng)到測井曲線上的深度點h₁，h₂，…，h_n(其中n∈N₊，N₊為正整數(shù))。

步驟2：在上述各二級層序界面對應(yīng)深度點h₁，h₂，…，h_n處對原始測井曲線A進行一級重構(gòu)，得到曲線B，成為一級重構(gòu)曲線。曲線B的重構(gòu)公式為：

S_B(h)＝S_A(h)+S_Amax/8+(S_Amax/8)*(-1)ⁿ⁺¹，h_n<h<h_(n+1)，n∈N₊(N₊為正整數(shù)) (1)

式中：S_B(h)為重構(gòu)后的曲線B在各深度點處的測井值，S_A(h)為曲線A在各深度點處的原測井值，S_Amax為測井曲線的全局最大值，即原始測井曲線A整個深度段的最大測井值，h_n為二級層序界面對應(yīng)的不同深度點，n為對應(yīng)的序號。公式(1)中參數(shù)(S_Amax/8)*(-1)ⁿ⁺¹的作用在于將二級層序界面對應(yīng)的深度點h_n之上的第n個二級層序和之下的第n+1個二級層序分別賦予不同的基準值，相鄰兩個二級層序的基準值相差1/4個S_Amax(依據(jù)小波變換試驗的經(jīng)驗值，相鄰二級層序基準值相差1/4個S_Amax時，在小波能譜圖上可進行明顯區(qū)分)，這樣在后續(xù)進行小波分析時可明顯將不同二級層序區(qū)分開。本發(fā)明實施例中，公式(1)的主要作用是將地震識別的二級層序劃分界限定位到測井曲線上，以準確控制內(nèi)部較短周期層序的劃分。

步驟3：依據(jù)地震剖面局部明顯的上超下截的地震反射特征，識別出三級層序界面；同時根據(jù)合成記錄的時深關(guān)系將地震上三級層序界面的時間點t_n1，t_n2，…，t_nm(t_nm為第n個二級層序內(nèi)第m個三級層序界面的時間點，n表示二級層序的次序，m表示三級層序界面的次序)對應(yīng)到測井曲線上的深度點h_n1，h_n2，…，h_nm(h_nm為第n個二級層序內(nèi)第m個三級層序界面的深度點，n表示二級層序的次序，m表示三級層序界面的次序)(其中m、n∈N₊，N₊為正整數(shù))。

步驟4：在上述三級層序界面對應(yīng)深度點h_n1，h_n2，…，h_nm(其中m、n∈N₊，N₊為正整數(shù))處對曲線B進行二級重構(gòu)，得到曲線C，稱為二級重構(gòu)曲線。曲線C的重構(gòu)公式為：

S_C(h)＝S_B(h)-S_Amax(n)/16+(S_Amax(n)/16)*(-1)^m，h_nm<h<h_n(m+1)，n，m∈N₊ (2)

式中：S_C(h)為曲線C在各深度點處的測井值，S_B(h)為重構(gòu)后的曲線B在各深度點處的測井值，S_Amax(n)為原始測井曲線的局部最大值，是曲線A在第n個二級層序內(nèi)部的最大測井值，h_nm為某一個二級層序n+1內(nèi)部第m個三級層序界面對應(yīng)的深度點，h_n(m+1)為某一個二級層序n+1內(nèi)部第m+1個三級層序界面對應(yīng)的深度點。公式(2)中參數(shù)(S_Amax(n)/16)*(-1)^m+1的含義為在于將某一個二級層序n內(nèi)第m個三級層序界面對應(yīng)的深度點h_nm之上的第m個三級層序和之下的第m+1個三級層序賦予不同的基準值，相鄰兩個三級層序之間的基準值相差1/8個S_Amax(n)(依據(jù)小波變換試驗的經(jīng)驗值，相鄰三級層序基準值相差1/8個S_Amax(n)時，在小波能譜圖上可進行明顯區(qū)分)，這樣在后續(xù)進行小波分析時可明顯將相同二級層序內(nèi)的三級層序區(qū)分開。公式(2)的主要作用是將地震識別的三級層序界面定位到測井曲線上，以準確控制其內(nèi)部較短周期層序的劃分。

步驟5：對重構(gòu)的曲線C進行6000尺度的小波分析，得到小波能量譜系圖，將小波能量譜系圖上的間斷點作為二級層序的分界點，進而對重構(gòu)曲線C中包含的二級層序進行劃分。

步驟6：對重構(gòu)的曲線C進行1400尺度的小波分析，得到小波能量譜系圖，將小波能量譜系圖上的間斷點作為三級層序分界點，進而對重構(gòu)曲線C中包含的三級層序進行劃分。

步驟7：對重構(gòu)的曲線C進行700尺度的小波分析，得到小波能量譜系圖，將小波能量譜系圖上的間斷點作為四級層序分界點，進而對重構(gòu)曲線C中包含的四級層序進行劃分。

步驟8：對重構(gòu)的曲線C進行150尺度的小波分析，得到小波系數(shù)曲線圖，將小波系數(shù)曲線圖上的間斷點作為五級層序的分界點，進而對重構(gòu)曲線C中包含的五級層序進行劃分。

上述本發(fā)明實施例中，進行小波分析中尺度的選擇為經(jīng)驗值，本發(fā)明的實現(xiàn)并不以此為限。

本發(fā)明針對長期低頻層序，依據(jù)地震剖面上的明顯反射特征對其識別，通過合成記錄將其轉(zhuǎn)換到測井曲線上對應(yīng)位置，再結(jié)合測井曲線本身良好的識別短周期高頻層序的能力，對層序進行劃分。具有以下兩個顯著優(yōu)點：

一、對長期低頻層序利用地震識別，再通過合成記錄定位到測井上；這樣通過地震對長期層序的準確識別保證了層序曲線中長期低頻層序劃分的結(jié)果準確可靠。

二、依據(jù)地震對測井曲線經(jīng)過兩次重構(gòu)得到二級重構(gòu)曲線，最后直接對二級重構(gòu)曲線進行不同尺度的小波變換分析，得到不同級別層序的劃分方案，其實現(xiàn)了層序劃分過程中真正的井震結(jié)合。

下面結(jié)合具體的實例對本發(fā)明技術(shù)方案做進一步詳細說明。

在遼河坳陷灘海地區(qū)選取海南洼陷新生界進行研究為例，海南洼陷新生界包括古近系和新近系，其中古近系發(fā)育東一段(d₁)、東二段(d₂)、東三段(d₃)、沙一段(s₁)、沙二段(s₂)以及沙三段(s₃)地層，沙四段(s₄)地層不發(fā)育；新近系發(fā)育館陶組(Ng)和明化鎮(zhèn)組(Nm)地層。進行地層層序劃分包括：

(1)如圖2所示，在地震剖面上依據(jù)區(qū)域和局部的上超下削特征，可在古近系識別出1個二級層序界面(Ng組底)和3個三級層序界面(Ng__S亞段底、d₃段底和s₂段底)。

(2)選取地震剖面上的井well1作為實施例，Well1井的測量深度范圍為190.00-3600.00m。對well1井制作合成記錄，將上述地震上識別出的二、三級層序界面定位到井well1上的對應(yīng)深度，其中1個二級層序界面Ng組底對應(yīng)的深度為1167.00m，3個三級層序界面Ng__S亞段底、E₃d₃段底和E₃s₂段底對應(yīng)的深度分別為721.00m、2230.50m和2971.00m，本發(fā)明實施例中圖2是地震剖面，其左側(cè)值并不是深度值，是時間值；Well1井的單位是深度單位。其中，a：明顯削截，b1、b2：區(qū)域不整合，c1、c2:局部不整合，d：強反射，e:古近系或中生界底。

(3)選取井Well1的AC曲線(本實施例中該井的自然伽馬測井GR、電阻率測井RT等曲線不完整)作為原始測井曲線A，找出其全局最大值為S_Amax＝189us/m，依據(jù)公式(1)及二級層序界面深度點計算得到重構(gòu)的一級重構(gòu)曲線B。

(4)分別在以1167.00m為界的兩個二級層序段內(nèi)找出其局部最大值為S_Amax(1)＝55us/m，S_Amax(2)＝50us/m。依據(jù)公式(2)及三級層序界面深度點計算得到重構(gòu)的二級重構(gòu)曲線C，如圖3中(a)所示，well1井AC曲線二級重構(gòu)曲線。

(5)本發(fā)明實施例中使用Matlab軟件中的小波工具箱對二級重構(gòu)曲線C進行6000尺度的morl小波(在Matlab小波變換工具箱中，經(jīng)過對db小波、sym小波、haar小波以及morl小波進行試驗分析，發(fā)現(xiàn)選擇morl小波進行小波變換的分析結(jié)果最好)變換，得到小波能量譜系圖，如圖3中(b)，曲線C6000尺度的小波變換能量譜系圖。不同層序的沉積環(huán)境不同，反映在小波能量譜系圖上即為不同層序形成的能量團之間存在明顯間斷點，且同一層序內(nèi)部能量線由下至上收斂。因此，從圖3中(b)中可以看出，在標注1處1167.00m處對應(yīng)著Ng組底界的位置，可以很明顯地將能量分為兩部分，分別代表兩個不同的二級層序，即將well1井新生界劃分為兩個二級層序，其與海南洼陷本身的層序地層特征相吻合。而針對原始曲線A同樣進行6000尺度的morl小波變換，如圖3中(c)所示，曲線A6000尺度的小波變換能量譜系圖，可以看出，在原始曲線A的6000尺度的能量譜系圖上，無法在將整體能量團在1167.00m分為兩部分，即難以對二級層序進行準確劃分。

(6)對二級重構(gòu)曲線C進行1400尺度(經(jīng)驗值)的morl小波分析，得到小波能量譜系圖，圖4中(a)為well1井AC曲線二級重構(gòu)曲線，圖4中(b)為曲線C1400尺度小波變換能量譜系圖，可知，3個三級層序的分界點可以清楚識別，即可以對三級層序進行準確劃分；結(jié)合二級層序分界點，可將well1井新生界劃分為5個三級層序，標注1為721.00m，標注2為2230.50m，標注3為2971.00m。

(7)在上述二、三級層序準確劃分的基礎(chǔ)上，開始進行較短周期層序的識別與劃分。對二級重構(gòu)曲線C進行700尺度(經(jīng)驗值)的morl小波分析得到小波能量譜系圖，如圖5所示。圖5中(a)為well1井AC曲線二級重構(gòu)曲線，圖5中(b)為曲線C700尺度小波變換能量譜系圖，根據(jù)圖5中(b)的能量分布，可以識別出4個四級層序分界點；結(jié)合前述二、三級層序分界點，可以將well1井新生界劃分為9個四級層序，標注1表示深度為487.25m，標注2表示1676.50m，標注3表示2683.75m，標注4表示3450.25m。

(8)對于短期層序的識別，由于其頻率高，而能量譜系圖的分辨率較低，因此采用小波系數(shù)曲線圖進行分析。小波系數(shù)曲線的突變處可對應(yīng)各級層序界面、湖泛面以及相轉(zhuǎn)換面，反映沉積環(huán)境的變化。對二級重構(gòu)曲線C進行150尺度(經(jīng)驗值)的morl小波分析，得到小波系數(shù)曲線圖，如圖所示，曲線a為well1井AC曲線二級重構(gòu)曲線，b為曲線C150尺度小波系數(shù)曲線，根據(jù)圖6中的小波曲線變化特征，尋找出10個突變點，結(jié)合前述二、三、四級層序分界點，可以將well1井新生界分為19個五級層序。

綜合以上8步的結(jié)果，可以準確地對海南洼陷well1井新生界進行層序劃分。

本發(fā)明通過在地震上依據(jù)區(qū)域不整合和局部不整合分別對二、三級層序進行識別，再選擇剖面上所含的井，制作合成記錄，將地震上識別的二三級層序分界點定位到測井曲線上對應(yīng)的深度點。在對應(yīng)深度點處通過經(jīng)驗公式對測井曲線進行重構(gòu)，重構(gòu)的目的是為了將二、三級層序界面特征放大，以此在準確劃分長期層序的基礎(chǔ)上控制內(nèi)部短期層序劃分的準確性。通過曲線重構(gòu)得到可以直接用來進行小波變化劃分層序的重構(gòu)曲線，針對不同級別的層序劃分，對重構(gòu)曲線選用不同尺度的小波變換，得到對應(yīng)的小波能量譜系圖或小波系數(shù)曲線圖，通過能譜圖和系數(shù)曲線圖中的間斷點來識別各級層序分界點，最終可對研究區(qū)發(fā)育的層序旋回進行準確的劃分。

同時，本發(fā)明還公開一種地層層序劃分裝置，如圖7所示，該裝置包括：

識別模塊701，用于根據(jù)地震剖面的地震反射特征識別二級層序界面、三級層序界面；

定位模塊702，用于將所述二級層序界面、三級層序界面分別定位到所述地震剖面中的一測井的對應(yīng)深度；

一級重構(gòu)模塊703，用于在所述二級層序界面的對應(yīng)深度點處對所述測井的測井曲線進行一級重構(gòu)，生成一級重構(gòu)曲線；

二級重構(gòu)模塊704，在所述三級層序界面的對應(yīng)深度點處對所述一級重構(gòu)曲線進行二級重構(gòu)，生成二級重構(gòu)曲線；

地層劃分模塊705，用于按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成地層層序劃分結(jié)果。

如圖8所示，本發(fā)明實施例中的定位模塊702包括：

二級層序定位單元7021，根據(jù)所述地震剖面中的測井的合成記錄的時深關(guān)系分別將所述二級層序界面時間點t₁，t₂，…，t_n對應(yīng)到所述測井的測井曲線上的深度點h₁，h₂，…，h_n；

三級層序定位單元7022，根據(jù)所述合成記錄的時深關(guān)系將所述三級層序界面的時間點t_n1，t_n2，…，t_nm對應(yīng)到測井曲線上的深度點h_n1，h_n2，…，h_nm；

其中，t_nm為第n個二級層序內(nèi)第m個三級層序界面的時間點，n表示二級層序的次序，m表示三級層序界面的次序，h_nm為第n個二級層序內(nèi)第m個三級層序界面的深度點，n表示二級層序的次序，m表示三級層序界面的次序，m、n∈N₊，N₊為正整數(shù)。

如圖9所示，地層劃分模塊705包括：

小波分析單元7051，用于按預(yù)設(shè)尺度對所述二級重構(gòu)曲線進行小波分析，生成小波能量譜系圖和小波系數(shù)曲線圖；

劃分單元7052，用于根據(jù)所述小波能量譜系圖和小波系數(shù)曲線圖的間斷點對所述二級重構(gòu)曲線進行層序劃分。

本發(fā)明公開的裝置解決問題的原理與上述方法的實現(xiàn)相似，不再贅述。

本領(lǐng)域內(nèi)的技術(shù)人員應(yīng)明白，本發(fā)明的實施例可提供為方法、系統(tǒng)、或計算機程序產(chǎn)品。因此，本發(fā)明可采用完全硬件實施例、完全軟件實施例、或結(jié)合軟件和硬件方面的實施例的形式。而且，本發(fā)明可采用在一個或多個其中包含有計算機可用程序代碼的計算機可用存儲介質(zhì)(包括但不限于磁盤存儲器、CD-ROM、光學(xué)存儲器等)上實施的計算機程序產(chǎn)品的形式。

本發(fā)明是參照根據(jù)本發(fā)明實施例的方法、設(shè)備(系統(tǒng))、和計算機程序產(chǎn)品的流程圖和/或方框圖來描述的。應(yīng)理解可由計算機程序指令實現(xiàn)流程圖和/或方框圖中的每一流程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結(jié)合?？商峁┻@些計算機程序指令到通用計算機、專用計算機、嵌入式處理機或其他可編程數(shù)據(jù)處理設(shè)備的處理器以產(chǎn)生一個機器，使得通過計算機或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的裝置。

這些計算機程序指令也可存儲在能引導(dǎo)計算機或其他可編程數(shù)據(jù)處理設(shè)備以特定方式工作的計算機可讀存儲器中，使得存儲在該計算機可讀存儲器中的指令產(chǎn)生包括指令裝置的制造品，該指令裝置實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能。

這些計算機程序指令也可裝載到計算機或其他可編程數(shù)據(jù)處理設(shè)備上，使得在計算機或其他可編程設(shè)備上執(zhí)行一系列操作步驟以產(chǎn)生計算機實現(xiàn)的處理，從而在計算機或其他可編程設(shè)備上執(zhí)行的指令提供用于實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的步驟。

本發(fā)明中應(yīng)用了具體實施例對本發(fā)明的原理及實施方式進行了闡述，以上實施例的說明只是用于幫助理解本發(fā)明的方法及其核心思想；同時，對于本領(lǐng)域的一般技術(shù)人員，依據(jù)本發(fā)明的思想，在具體實施方式及應(yīng)用范圍上均會有改變之處，綜上所述，本說明書內(nèi)容不應(yīng)理解為對本發(fā)明的限制。