碳酸鹽巖縫洞雕刻方法及裝置制造方法
【專利摘要】本發(fā)明實(shí)施例提供了一種碳酸鹽巖縫洞雕刻方法及裝置，其中，該方法包括：將碳酸鹽巖的工區(qū)平面圖劃分為多個縫洞帶，并將每個縫洞帶劃分為多個縫洞單元，所述縫洞帶為不規(guī)則的閉合區(qū)域，所述縫洞單元是不規(guī)則的閉合區(qū)域；根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體；計(jì)算每個縫洞帶中每個縫洞單元內(nèi)每個縫洞體的體積和儲量；對于每個縫洞單元，將該縫洞單元內(nèi)縫洞體的體積累加，得到該縫洞單元的體積，將該縫洞單元內(nèi)縫洞體的儲量累加，得到該縫洞單元的儲量。該方案可以一次性將縫洞帶內(nèi)的縫洞體全部雕刻出來，并量化縫洞體的體積和儲量。
【專利說明】碳酸鹽巖縫洞雕刻方法及裝置

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地球物理勘探地震資料計(jì)算機(jī)識別【技術(shù)領(lǐng)域】，特別涉及一種碳酸鹽巖 縫洞雕刻方法及裝置。

【背景技術(shù)】
[0002] 在地球物理勘探處理和解釋過程中，對碳酸鹽巖縫洞進(jìn)行雕刻、量化統(tǒng)計(jì)并形成 報(bào)表以滿足勘探生產(chǎn)開發(fā)的需要，是碳酸鹽巖區(qū)域必須解決的問題。
[0003] 只有將碳酸鹽巖縫洞雕刻出來，統(tǒng)計(jì)分布范圍并量化其體積、儲量才能進(jìn)行下一 步的勘探生產(chǎn)和開發(fā)。但是，目前針對碳酸鹽巖這樣的特殊地質(zhì)目標(biāo)，常規(guī)解釋軟件無法滿 足量化縫洞的體積、儲量等信息需要，無法實(shí)現(xiàn)工業(yè)化生產(chǎn)。


【發(fā)明內(nèi)容】

[0004] 本發(fā)明實(shí)施例提供了一種碳酸鹽巖縫洞雕刻方法，以滿足量化縫洞的體積、儲量 的需要，該方法包括：將碳酸鹽巖的工區(qū)平面圖劃分為多個縫洞帶，并將每個縫洞帶劃分為 多個縫洞單元，所述縫洞帶為不規(guī)則的閉合區(qū)域，所述縫洞單元是不規(guī)則的閉合區(qū)域；根據(jù) 預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體；計(jì)算每個縫洞帶中每 個縫洞單元內(nèi)每個縫洞體的體積和儲量；對于每個縫洞單元，將該縫洞單元內(nèi)縫洞體的體 積累加，得到該縫洞單元的體積，將該縫洞單元內(nèi)縫洞體的儲量累加，得到該縫洞單元的儲 量。
[0005] 在一個實(shí)施例中，根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi) 的縫洞體，包括：對于每個縫洞帶，獲取該縫洞帶的線號個數(shù)、道號個數(shù)和樣點(diǎn)個數(shù)；將該 縫洞帶內(nèi)每個種子點(diǎn)的瞬時振幅值存儲在三維數(shù)組中，并為瞬時振幅值在所述預(yù)設(shè)縫洞振 幅值范圍內(nèi)的種子點(diǎn)設(shè)置第一特征值，為瞬時振幅值在所述預(yù)設(shè)縫洞振幅值范圍外的種子 點(diǎn)設(shè)置第二特征值；以該縫洞帶內(nèi)第一個具有第一特征值的種子點(diǎn)為中心定義一個縫洞 體；在該種子點(diǎn)的上方、下方、左方、右方、前方和后方與該種子點(diǎn)聯(lián)通的種子點(diǎn)中，將具有 第一特征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，并以該確定出的種子點(diǎn)為中 心，從與該確定出的種子點(diǎn)上方、下方、左方、右方、前方和后方聯(lián)通的種子點(diǎn)中，將具有第 一特征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，直至與該確定出的種子點(diǎn)聯(lián)通的 種子點(diǎn)中，不包括具有第一特征值的種子點(diǎn)。
[0006] 在一個實(shí)施例中，計(jì)算每個縫洞帶中每個縫洞單元內(nèi)縫洞體的體積，包括：根據(jù)工 區(qū)的線線間距、道號間距和采樣間隔，計(jì)算縫洞體內(nèi)每個種子點(diǎn)的體積；累加縫洞體內(nèi)每個 種子點(diǎn)的體積，得到縫洞體的體積。
[0007] 在一個實(shí)施例中，計(jì)算每個縫洞帶中每個縫洞單元內(nèi)縫洞體的儲量，包括：根據(jù) 縫洞體內(nèi)每個種子點(diǎn)的體積和縫洞體內(nèi)每個種子點(diǎn)的孔隙度數(shù)據(jù)值，計(jì)算縫洞體的有效體 積；根據(jù)含油飽和度、密度和縫洞體的有效體積，計(jì)算縫洞體的儲量。
[0008] 在一個實(shí)施例中，根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi) 的縫洞體，包括：根據(jù)用戶輸入的預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi) 的縫洞體；所述方法還包括：確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體之前，對所述工 區(qū)的地震資料數(shù)據(jù)進(jìn)行瞬時振幅變換，得到瞬時振幅數(shù)據(jù)，所述瞬時振幅數(shù)據(jù)為用戶確定 所述預(yù)設(shè)縫洞振幅值范圍提供依據(jù)。
[0009] 在一個實(shí)施例中，還包括：輸出外部文件，所述外部文件包括每個縫洞體的體積、 儲量、位置、面積和厚度。
[0010] 本發(fā)明實(shí)施例提供了一種碳酸鹽巖縫洞雕刻裝置，以滿足量化縫洞的體積、儲量 的需要。該裝置包括：劃分模塊，用于將碳酸鹽巖的工區(qū)平面圖劃分為多個縫洞帶，并將每 個縫洞帶劃分為多個縫洞單元，所述縫洞帶為不規(guī)則的閉合區(qū)域，所述縫洞單元是不規(guī)則 的閉合區(qū)域；確定模塊，用于根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi) 的縫洞體；計(jì)算模塊，用于計(jì)算每個縫洞帶中每個縫洞單元內(nèi)每個縫洞體的體積和儲量； 累加模塊，用于對于每個縫洞單元，將該縫洞單元內(nèi)縫洞體的體積累加，得到該縫洞單元的 體積，將該縫洞單元內(nèi)縫洞體的儲量累加，得到該縫洞單元的儲量。
[0011] 在一個實(shí)施例中，所述確定模塊，包括：獲取單元，用于對于每個縫洞帶，獲取該縫 洞帶的線號個數(shù)、道號個數(shù)和樣點(diǎn)個數(shù)；設(shè)置單元，用于將該縫洞帶內(nèi)每個種子點(diǎn)的瞬時振 幅值存儲在三維數(shù)組中，并為瞬時振幅值在所述預(yù)設(shè)縫洞振幅值范圍內(nèi)的種子點(diǎn)設(shè)置第一 特征值，為瞬時振幅值在所述預(yù)設(shè)縫洞振幅值范圍外的種子點(diǎn)設(shè)置第二特征值；定義單元， 用于以該縫洞帶內(nèi)第一個具有第一特征值的種子點(diǎn)為中心定義一個縫洞體；雕刻單元，用 于在該種子點(diǎn)的上方、下方、左方、右方、前方和后方與該種子點(diǎn)聯(lián)通的種子點(diǎn)中，將具有第 一特征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，并以該確定出的種子點(diǎn)為中心， 從與該確定出的種子點(diǎn)上方、下方、左方、右方、前方和后方聯(lián)通的種子點(diǎn)中，將具有第一特 征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，直至與該確定出的種子點(diǎn)聯(lián)通的種子 點(diǎn)中，不包括具有第一特征值的種子點(diǎn)。
[0012] 在一個實(shí)施例中，所述計(jì)算模塊，包括：第一計(jì)算單元，用于根據(jù)工區(qū)的線線間距、 道號間距和采樣間隔，計(jì)算縫洞體內(nèi)每個種子點(diǎn)的體積；第二計(jì)算單元，用于累加縫洞體內(nèi) 每個種子點(diǎn)的體積，得到縫洞體的體積。
[0013] 在一個實(shí)施例中，所述計(jì)算模塊，還包括：第三計(jì)算單元，用于根據(jù)縫洞體內(nèi)每個 種子點(diǎn)的體積和縫洞體內(nèi)每個種子點(diǎn)的孔隙度數(shù)據(jù)值，計(jì)算縫洞體的有效體積；第四計(jì)算 單元，用于根據(jù)含油飽和度、密度和縫洞體的有效體積，計(jì)算縫洞體的儲量。
[0014] 在一個實(shí)施例中，所述確定模塊，具體用于根據(jù)用戶輸入的預(yù)設(shè)縫洞振幅值范圍， 確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體；所述裝置還包括：變換模塊，用于根據(jù)用戶 輸入的預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體之前，對所述工 區(qū)的地震資料數(shù)據(jù)進(jìn)行瞬時振幅變換，得到瞬時振幅數(shù)據(jù)，所述瞬時振幅數(shù)據(jù)為用戶確定 所述預(yù)設(shè)縫洞振幅值范圍提供依據(jù)。
[0015] 在一個實(shí)施例中，還包括：輸出模塊，用于輸出外部文件，所述外部文件的內(nèi)容包 括每個縫洞體的體積、儲量、位置、面積和厚度。
[0016] 在本發(fā)明實(shí)施例中，通過將碳酸鹽巖的工區(qū)平面圖劃分為多個縫洞帶，并將每個 縫洞帶劃分為多個縫洞單元，然后，根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞 單元內(nèi)的縫洞體，并計(jì)算每個縫洞單元內(nèi)縫洞體的體積和儲量，最后得到每個縫洞單元的 體積和儲量，實(shí)現(xiàn)了一次性將縫洞帶內(nèi)的縫洞體全部雕刻出來，并量化縫洞體的體積和儲 量，使得可以對縫洞體進(jìn)行快速、便捷地雕刻，并使解釋人員對縫洞體能夠有一個更加直觀 的認(rèn)識，以滿足下一步勘探生產(chǎn)的需求。

【專利附圖】

【附圖說明】
[0017] 此處所說明的附圖用來提供對本發(fā)明的進(jìn)一步理解，構(gòu)成本申請的一部分，并不 構(gòu)成對本發(fā)明的限定。在附圖中：
[0018] 圖1是本發(fā)明實(shí)施例提供的一種碳酸鹽巖縫洞雕刻方法的流程圖；
[0019] 圖2是本發(fā)明實(shí)施例提供的一種縫洞帶劃分示意圖；
[0020] 圖3是本發(fā)明實(shí)施例提供的一種縫洞單元劃分示意圖；
[0021] 圖4是本發(fā)明實(shí)施例提供的一種縫洞的聯(lián)通示意圖；
[0022] 圖5是本發(fā)明實(shí)施例提供的一種縫洞的三維顯示效果圖；
[0023] 圖6是本發(fā)明實(shí)施例提供的一種碳酸鹽巖縫洞雕刻裝置的結(jié)構(gòu)框圖。

【具體實(shí)施方式】
[0024] 為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下面結(jié)合實(shí)施方式和附圖，對 本發(fā)明做進(jìn)一步詳細(xì)說明。在此，本發(fā)明的示意性實(shí)施方式及其說明用于解釋本發(fā)明，但并 不作為對本發(fā)明的限定。
[0025] 在本發(fā)明實(shí)施例中，提供了一種碳酸鹽巖縫洞雕刻方法，如圖1所示，該方法包 括：
[0026] 步驟101 :將碳酸鹽巖的工區(qū)平面圖劃分為多個縫洞帶，并將每個縫洞帶劃分為 多個縫洞單元，所述縫洞帶為不規(guī)則的閉合區(qū)域，所述縫洞單元是不規(guī)則的閉合區(qū)域；
[0027] 步驟102 :根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞 體；
[0028] 步驟103 :計(jì)算每個縫洞帶中每個縫洞單元內(nèi)每個縫洞體的體積和儲量；
[0029] 步驟104 :對于每個縫洞單元，將該縫洞單元內(nèi)縫洞體的體積累加，得到該縫洞單 元的體積，將該縫洞單元內(nèi)縫洞體的儲量累加，得到該縫洞單元的儲量。
[0030] 由圖1所示的流程可知，在本發(fā)明實(shí)施例中，通過將碳酸鹽巖的工區(qū)平面圖劃分 為多個縫洞帶，并將每個縫洞帶劃分為多個縫洞單元，然后，根據(jù)預(yù)設(shè)縫洞振幅值范圍，確 定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體，并計(jì)算每個縫洞單元內(nèi)縫洞體的體積和儲量， 最后得到每個縫洞單元的體積和儲量，實(shí)現(xiàn)了一次性將縫洞帶內(nèi)的縫洞體全部雕刻出來， 并量化縫洞體的體積和儲量，使得可以對縫洞體進(jìn)行快速、便捷地雕刻，并使解釋人員對縫 洞體能夠有一個更加直觀的認(rèn)識，以滿足下一步勘探生產(chǎn)的需求。
[0031] 具體實(shí)施時，可以通過以下方式將工區(qū)平面圖劃分為多個縫洞帶，并將每個縫洞 帶劃分為多個縫洞單元：例如，解釋用戶在工區(qū)平面圖上用鼠標(biāo)圈出閉合縫洞帶，如圖2所 示，劃分為縫洞帶1和縫洞帶2,這些縫洞帶用來圈定需要雕刻縫洞的工區(qū)邊界，表現(xiàn)為多 個閉合的多邊形，每個縫洞帶的名字唯一；使用同樣方法定義縫洞單元，用鼠標(biāo)在縫洞帶內(nèi) 圈出縫洞單元，如圖3所示，縫洞單元指縫洞帶內(nèi)需要進(jìn)一步統(tǒng)計(jì)縫洞信息的邊界范圍。
[0032] 具體實(shí)施時，可以通過以下步驟來實(shí)現(xiàn)根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞 帶中每個縫洞單元內(nèi)的縫洞體：
[0033] 首先，對于每個縫洞帶，獲取該縫洞帶的線號個數(shù)、道號個數(shù)和樣點(diǎn)個數(shù)；例如， 搜索某個指定的縫洞帶所在的工區(qū)線道號最小范圍（line_min，line_max，cmp_min，cmp_ max)，搜索目的層所在的采樣索引范圍（sample_min，sample_max)，則該縫洞帶的線號個 數(shù)、道號個數(shù)、樣點(diǎn)個數(shù)分別為：
[0034] line_count = line_max - line_min+l ；
[0035] cmp_count = cmp_max - cmp_min+l ;
[0036] sample_count = sample_max - sample_min+l〇
[0037] 其次，將該縫洞帶內(nèi)每個種子點(diǎn)的瞬時振幅值存儲在三維數(shù)組中，并為瞬時振幅 值在所述預(yù)設(shè)縫洞振幅值范圍內(nèi)的種子點(diǎn)設(shè)置第一特征值，為瞬時振幅值在所述預(yù)設(shè)縫洞 振幅值范圍外的種子點(diǎn)設(shè)置第二特征值；例如，定義bool類型三維數(shù)組，以存放該縫洞帶 內(nèi)每個種子點(diǎn)的瞬時振幅特征值，并將特征數(shù)組中每個特征點(diǎn)初始化為false :
[0038] bool flag[line_count, cmp_count, sample_count]
[0039] 對每一個flag元素賦值為false ;
[0040] 依次讀取該縫洞帶范圍內(nèi)每個種子點(diǎn)的瞬時振幅數(shù)據(jù)，根據(jù)瞬時振幅值是否位于 用戶定義的縫洞振幅范圍（即上述預(yù)設(shè)縫洞振幅值范圍）內(nèi)，對振幅特征數(shù)組進(jìn)行賦值，假 設(shè)某個樣點(diǎn)線號為line_no,道號為cmp_no,樣點(diǎn)號為sample_no,那么該點(diǎn)在所加載數(shù)據(jù) 范圍內(nèi)的線索引、道索引、樣點(diǎn)索引分別為：
[0041] line_idx = line_no - line_min ；
[0042] cmp_idx = cmp_no - cmp_min ;
[0043] sample_idx = sample_no - sample_min ;
[0044] 該位置的特征值根據(jù)用戶定義的縫洞振幅范圍確定，如果該樣點(diǎn)的瞬時振幅值在 所述預(yù)設(shè)縫洞振幅值范圍內(nèi)，則為該樣點(diǎn)設(shè)置的特征值為true (即第一特征值），否則，為 該樣點(diǎn)設(shè)置的特征值為false (第二特征值）。
[0045] 再次，以該縫洞帶內(nèi)第一個具有第一特征值的種子點(diǎn)為中心定義一個縫洞體；例 如，分別對該縫洞帶內(nèi)每一個flag特征點(diǎn)進(jìn)行判斷，如果該樣點(diǎn)的特征值為false，則忽略 該特征點(diǎn)，否則定義一個縫洞體，并對位于（line_i，cmp_i，sample_i)位置的種子特征點(diǎn) 執(zhí)行雕刻操作。
[0046] 最后，對具有第一特征值的種子點(diǎn)執(zhí)行以下步驟進(jìn)行縫洞的雕刻：在該種子點(diǎn)的 上方、下方、左方、右方、前方和后方與該種子點(diǎn)聯(lián)通的種子點(diǎn)中，將具有第一特征值的種子 點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，并以該確定出的種子點(diǎn)為中心，從與該確定出的 種子點(diǎn)上方、下方、左方、右方、前方和后方聯(lián)通的種子點(diǎn)中，將具有第一特征值的種子點(diǎn)確 定為屬于所述一個縫洞體的種子點(diǎn)，直至與該確定出的種子點(diǎn)聯(lián)通的種子點(diǎn)中，不包括具 有第一特征值的種子點(diǎn)。例如，首先將flag(line_i，cmp_i，sample」）的特征值賦值為 false,防止該點(diǎn)被再次分析，并將種子特征點(diǎn)（line」，cmp_i，sample」）的索引保存到縫 洞體中，取出位于下列位置種子點(diǎn)的特征值：
[0047] (line」，cmp_i，sample_i_l)
[0048] (line_i, cmp_i, sample_i+l)
[0049] (line_i, cmp_i-l, sample_i)
[0050] (line」，cmp_i+l，sample」）
[0051] (line_i_l，cmp_i，sample」）
[0052] (line_i+l，cmp_i，sample」）
[0053] 上面點(diǎn)分別表示種子特征點(diǎn)（line_i，cmp_i，sample_i)上、下、左、右、前、后六個 點(diǎn)，即與種子特征點(diǎn)（line」，cmp_i, sample」）上方、下方、左方、右方、前方和后方聯(lián)通的 種子點(diǎn)（如圖4所示），將這些種子點(diǎn)中特征值為true的種子點(diǎn)，確定為同一個縫洞體，將 這些種子點(diǎn)中特征值為false的種子點(diǎn)，確定為不是同一個縫洞體，并以確定出的種子點(diǎn) 為中心，從與該確定出的種子點(diǎn)聯(lián)通的種子點(diǎn)中，將具有第一特征值（即true)的種子點(diǎn)確 定為屬于一個縫洞體的種子點(diǎn)，直至與該確定出的種子點(diǎn)聯(lián)通的種子點(diǎn)中，不包括具有第 一特征值的種子點(diǎn)，該一個縫洞體雕刻完畢。繼續(xù)判斷縫洞帶內(nèi)其他種子點(diǎn)的特征值，并以 具有第一特征值的種子點(diǎn)為中心定義另一個縫洞體，對該具有第一特征值的種子點(diǎn)進(jìn)行縫 洞雕刻的操作。
[0054] 具體的，根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體， 包括：根據(jù)用戶輸入的預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體， 為了便于用戶給出預(yù)設(shè)縫洞振幅值范圍，在本實(shí)施例中，根據(jù)用戶輸入的預(yù)設(shè)縫洞振幅值 范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體之前，還包括：對所述工區(qū)的地震資料數(shù) 據(jù)進(jìn)行瞬時振幅變換，得到瞬時振幅數(shù)據(jù)，所述瞬時振幅數(shù)據(jù)為用戶確定所述預(yù)設(shè)縫洞振 幅值范圍提供依據(jù)。因?yàn)榭p洞在地震資料中的表現(xiàn)是一個個的最大或最小振幅范圍，即縫 洞的振幅特征為大于某個振幅值和小于某個振幅值兩個區(qū)間，這樣的振幅值范圍很難讓用 戶指定所雕刻的振幅值范圍，所以在做縫洞雕刻之前，對地震資料做瞬時振幅變換。
[0055] 首先，對地震數(shù)據(jù)進(jìn)行Hilbert變換，假設(shè)地震振幅值為s (t):
[0056] .s (t) = Hjs! =/;(t)xs(t)= i s(r)/7(t-r)dr = - f , J一 π x- τ
[0057] 則瞬時振幅計(jì)算公式如下：
[0058] ^(t) = ^/s；(i) + s:(t) '
[0059] 對地震資料做了瞬時振幅變換后，其最大、最小振幅區(qū)間都變成了某個最大振幅 區(qū)間，這樣的振幅區(qū)間方便用戶指定某個振幅范圍（即預(yù)設(shè)縫洞振幅值范圍），比如大于某 個振幅值的樣點(diǎn)特征。
[0060] 在雕刻出縫洞體之后，還需要量化縫洞體的體積，在本實(shí)施例中，計(jì)算每個縫洞帶 中每個縫洞單元內(nèi)縫洞體的體積，包括：根據(jù)縫工區(qū)的線線間距、道號間距和采樣間隔，計(jì) 算縫洞體內(nèi)每個種子點(diǎn)的體積；累加縫洞體內(nèi)每個種子點(diǎn)的體積，得到縫洞體的體積。具 體的，根據(jù)工區(qū)線線間距（line_distance)道號間距（cmp_distance)采樣間隔（sampel_ interval),計(jì)算每個有效點(diǎn)（即縫洞體內(nèi)每個種子點(diǎn)）的體積其中l(wèi)ine_distance, cmp_distance單位為米，如果是深度域類型數(shù)據(jù)sample_interval單位為米，而時間域類 型數(shù)據(jù)則需要將時間間隔變換到深度域：
[0061] ? [深度域：sample-interval 1 .va/77/；/e _ interval = , - - (Jl、t 丨川域：sample _ interval x velocity I
[0062] 其中velocity指當(dāng)前目地層的速度，于某個特定目的層velocity可以認(rèn)為是一 個常量，該參數(shù)由用戶指定。
[0063] Vj = line_distanceXcmp_distanceXsample_distance
[0064] 那么每個縫洞的體積則為：V = 縫洞體內(nèi)種子點(diǎn)的個數(shù)。在雕刻出縫洞體之 后，還需要量化縫洞體的體積，在本實(shí)施例中，計(jì)算每個縫洞帶中每個縫洞單元內(nèi)縫洞體的 儲量，包括：根據(jù)縫洞體內(nèi)每個種子點(diǎn)的體積和縫洞體內(nèi)每個種子點(diǎn)的孔隙度數(shù)據(jù)值，計(jì)算 縫洞體的有效體積；根據(jù)含油飽和度、密度和縫洞體的有效體積，計(jì)算縫洞體的儲量。例如， 對縫洞體內(nèi)的每個種子點(diǎn)讀取相應(yīng)位置的孔隙度數(shù)據(jù)值P〇ri，該參數(shù)可以由外部作為孔隙 度數(shù)據(jù)體輸入，根據(jù)孔隙度數(shù)據(jù)值計(jì)算縫洞體的有效體積：V valid = / ViXpoiV
[0065] 那么縫洞體的儲量則為：R = VvalidX含油飽和度X密度，其中，含油飽和度和密 度可以通過測井獲得，由用戶輸入即可。
[0066] 最后，為了將雕刻結(jié)果呈現(xiàn)給用戶，并便于供其他軟件打開，在本實(shí)施例中，上述 碳酸鹽巖縫洞雕刻方法，還包括：輸出外部文件，所述外部文件包括每個縫洞體的體積、儲 量、位置、面積和厚度，該外部文件可以是報(bào)表的形式，如下表1所示。其中，縫洞體的位置 是縫洞體中心的種子點(diǎn)所在的位置，縫洞體的面積是縫洞體從上向下投影的面積。此外，還 可以使用三維可視化顯示雕刻結(jié)果（如圖5所示），以更加直觀的方式觀察縫洞的位置、體 積等信息。
[0067]

【權(quán)利要求】
1. 一種碳酸鹽巖縫洞雕刻方法，其特征在于，包括： 將碳酸鹽巖的工區(qū)平面圖劃分為多個縫洞帶，并將每個縫洞帶劃分為多個縫洞單元， 所述縫洞帶為不規(guī)則的閉合區(qū)域，所述縫洞單元是不規(guī)則的閉合區(qū)域； 根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體； 計(jì)算每個縫洞帶中每個縫洞單元內(nèi)每個縫洞體的體積和儲量； 對于每個縫洞單元，將該縫洞單元內(nèi)縫洞體的體積累加，得到該縫洞單元的體積，將該 縫洞單元內(nèi)縫洞體的儲量累加，得到該縫洞單元的儲量。
2. 如權(quán)利要求1所述碳酸鹽巖縫洞雕刻方法，其特征在于，根據(jù)預(yù)設(shè)縫洞振幅值范圍， 確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體，包括： 對于每個縫洞帶，獲取該縫洞帶的線號個數(shù)、道號個數(shù)和樣點(diǎn)個數(shù)； 將該縫洞帶內(nèi)每個種子點(diǎn)的瞬時振幅值存儲在三維數(shù)組中，并為瞬時振幅值在所述預(yù) 設(shè)縫洞振幅值范圍內(nèi)的種子點(diǎn)設(shè)置第一特征值，為瞬時振幅值在所述預(yù)設(shè)縫洞振幅值范圍 外的種子點(diǎn)設(shè)置第二特征值； 以該縫洞帶內(nèi)第一個具有第一特征值的種子點(diǎn)為中心定義一個縫洞體； 在該種子點(diǎn)的上方、下方、左方、右方、前方和后方與該種子點(diǎn)聯(lián)通的種子點(diǎn)中，將具 有第一特征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，并以該確定出的種子點(diǎn)為中 心，從與該確定出的種子點(diǎn)上方、下方、左方、右方、前方和后方聯(lián)通的種子點(diǎn)中，將具有第 一特征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，直至與該確定出的種子點(diǎn)聯(lián)通的 種子點(diǎn)中，不包括具有第一特征值的種子點(diǎn)。
3. 如權(quán)利要求2所述碳酸鹽巖縫洞雕刻方法，其特征在于，計(jì)算每個縫洞帶中每個縫 洞單元內(nèi)縫洞體的體積，包括： 根據(jù)工區(qū)的線線間距、道號間距和采樣間隔，計(jì)算縫洞體內(nèi)每個種子點(diǎn)的體積； 累加縫洞體內(nèi)每個種子點(diǎn)的體積，得到縫洞體的體積。
4. 如權(quán)利要求3所述碳酸鹽巖縫洞雕刻方法，其特征在于，計(jì)算每個縫洞帶中每個縫 洞單元內(nèi)縫洞體的儲量，包括： 根據(jù)縫洞體內(nèi)每個種子點(diǎn)的體積和縫洞體內(nèi)每個種子點(diǎn)的孔隙度數(shù)據(jù)值，計(jì)算縫洞體 的有效體積； 根據(jù)含油飽和度、密度和縫洞體的有效體積，計(jì)算縫洞體的儲量。
5. 如權(quán)利要求1至4中任一項(xiàng)所述碳酸鹽巖縫洞雕刻方法，其特征在于，根據(jù)預(yù)設(shè)縫洞 振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體，包括：根據(jù)用戶輸入的預(yù)設(shè)縫洞 振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體； 所述方法還包括： 確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞體之前，對所述工區(qū)的地震資料數(shù)據(jù)進(jìn)行瞬 時振幅變換，得到瞬時振幅數(shù)據(jù)，所述瞬時振幅數(shù)據(jù)為用戶確定所述預(yù)設(shè)縫洞振幅值范圍 提供依據(jù)。
6. 如權(quán)利要求1至4中任一項(xiàng)所述碳酸鹽巖縫洞雕刻方法，其特征在于，還包括： 輸出外部文件，所述外部文件包括每個縫洞體的體積、儲量、位置、面積和厚度。
7. -種碳酸鹽巖縫洞雕刻裝置，其特征在于，包括： 劃分模塊，用于將碳酸鹽巖的工區(qū)平面圖劃分為多個縫洞帶，并將每個縫洞帶劃分為 多個縫洞單元，所述縫洞帶為不規(guī)則的閉合區(qū)域，所述縫洞單元是不規(guī)則的閉合區(qū)域； 確定模塊，用于根據(jù)預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的縫洞 體； 計(jì)算模塊，用于計(jì)算每個縫洞帶中每個縫洞單元內(nèi)每個縫洞體的體積和儲量； 累加模塊，用于對于每個縫洞單元，將該縫洞單元內(nèi)縫洞體的體積累加，得到該縫洞單 元的體積，將該縫洞單元內(nèi)縫洞體的儲量累加，得到該縫洞單元的儲量。
8. 如權(quán)利要求7所述碳酸鹽巖縫洞雕刻裝置，其特征在于，所述確定模塊，包括： 獲取單元，用于對于每個縫洞帶，獲取該縫洞帶的線號個數(shù)、道號個數(shù)和樣點(diǎn)個數(shù)； 設(shè)置單元，用于將該縫洞帶內(nèi)每個種子點(diǎn)的瞬時振幅值存儲在三維數(shù)組中，并為瞬時 振幅值在所述預(yù)設(shè)縫洞振幅值范圍內(nèi)的種子點(diǎn)設(shè)置第一特征值，為瞬時振幅值在所述預(yù)設(shè) 縫洞振幅值范圍外的種子點(diǎn)設(shè)置第二特征值； 定義單元，用于以該縫洞帶內(nèi)第一個具有第一特征值的種子點(diǎn)為中心定義一個縫洞 體； 雕刻單元，用于在該種子點(diǎn)的上方、下方、左方、右方、前方和后方與該種子點(diǎn)聯(lián)通的種 子點(diǎn)中，將具有第一特征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，并以該確定出 的種子點(diǎn)為中心，從與該確定出的種子點(diǎn)上方、下方、左方、右方、前方和后方聯(lián)通的種子點(diǎn) 中，將具有第一特征值的種子點(diǎn)確定為屬于所述一個縫洞體的種子點(diǎn)，直至與該確定出的 種子點(diǎn)聯(lián)通的種子點(diǎn)中，不包括具有第一特征值的種子點(diǎn)。
9. 如權(quán)利要求8所述碳酸鹽巖縫洞雕刻裝置，其特征在于，所述計(jì)算模塊，包括： 第一計(jì)算單元，用于根據(jù)工區(qū)的線線間距、道號間距和采樣間隔，計(jì)算縫洞體內(nèi)每個種 子點(diǎn)的體積； 第二計(jì)算單元，用于累加縫洞體內(nèi)每個種子點(diǎn)的體積，得到縫洞體的體積。
10. 如權(quán)利要求9所述碳酸鹽巖縫洞雕刻裝置，其特征在于，所述計(jì)算模塊，還包括： 第三計(jì)算單元，用于根據(jù)縫洞體內(nèi)每個種子點(diǎn)的體積和縫洞體內(nèi)每個種子點(diǎn)的孔隙度 數(shù)據(jù)值，計(jì)算縫洞體的有效體積； 第四計(jì)算單元，用于根據(jù)含油飽和度、密度和縫洞體的有效體積，計(jì)算縫洞體的儲量。
11. 如權(quán)利要求7至10中任一項(xiàng)所述碳酸鹽巖縫洞雕刻裝置，其特征在于，所述確定模 塊，具體用于根據(jù)用戶輸入的預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單元內(nèi)的 縫洞體； 所述裝置還包括： 變換模塊，用于根據(jù)用戶輸入的預(yù)設(shè)縫洞振幅值范圍，確定每個縫洞帶中每個縫洞單 元內(nèi)的縫洞體之前，對所述工區(qū)的地震資料數(shù)據(jù)進(jìn)行瞬時振幅變換，得到瞬時振幅數(shù)據(jù)，所 述瞬時振幅數(shù)據(jù)為用戶確定所述預(yù)設(shè)縫洞振幅值范圍提供依據(jù)。
12. 如權(quán)利要求7至10中任一項(xiàng)所述碳酸鹽巖縫洞雕刻裝置，其特征在于，還包括： 輸出模塊，用于輸出外部文件，所述外部文件的內(nèi)容包括每個縫洞體的體積、儲量、位 置、面積和厚度。
【文檔編號】G01V1/28GK104155689SQ201410427787
【公開日】2014年11月19日 申請日期:2014年8月27日 優(yōu)先權(quán)日:2014年8月27日
【發(fā)明者】詹仕凡, 張瑋, 陶春峰, 李磊, 劉永雷, 呂東, 但光劍, 趙佳瑜, 萬忠宏 申請人:中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司