基於trap-3d軟體尋找油氣藏的方法
2023-05-28 07:52:41
基於trap-3d軟體尋找油氣藏的方法
【專利摘要】本發明提供了一種基於TRAP-3D軟體尋找油氣藏的方法，該方法包括：根據三維地震資料和測井資料建立勘探工區的三維巖性及斷層數據立方體；將三維巖性及斷層數據立方體劃分成若干等深度的深度切片，並分別對每個深度切片進行獨立砂體單元劃分；將三維巖性及斷層數據立方體的深度切片順序輸入TRAP-3D軟體進行油氣藏評價。本發明提高了三維圈閉評價的精確性，有利於精確尋找油氣藏，並且可繪出平面上的甜點圖，以及展示縱向不同深度上的油氣圈閉量，並可得到勘探工區內油氣藏的總圈閉量。
【專利說明】基於了以?-30軟體尋找油氣藏的方法
【技術領域】
[0001]本發明涉及油氣勘探【技術領域】，尤其是涉及一種基於軟體尋找油氣藏的方法。
【背景技術】
[0002]適合於油氣聚集、形成油氣藏的場所稱為圈閉(或油氣圈閉目前圈閉概念是尋找石油和天然氣的重要理論。沒有圈閉，油氣無法存在。但是長期以來，對於構造圈閉仍僅限於將構造隆起幅度和構造面積兩項作為其圈閉評價的主要內容。但僅從上述兩項指標難以定量衡量實際圈閉量的大小。而對於斷塊圈閉，通常僅根據斷塊的面積大小對圈閉面積進行估計，而對於該斷塊是否具備圈閉條件這一關鍵問題則缺乏精細研究。
[0003]目前斷層封堵性、構造溢出點以及圈閉量估算研究的商業軟體仍停留在定性分析或通過二維資料分析地下儲層的階段。其常規做法是在二維剖面上做斷層兩側砂體構造疊合圖，據此確定每個砂層的斷層面溢出點，找出各封堵邊界點；將上述封堵邊界點投影至構造等高線平面並連線，以此計算封堵面積和圈閉量。而運用上述常規方法對由許多斷層所形成的複雜斷塊的封堵性進行評價時，僅利用一系列斷層兩側巖性二維對接剖面很難準確分析三維空間中砂層的連通及其總體的封堵情況，從而難以得到精確結果。
【發明內容】
[0004]本發明的目的在於提供一種基於I狀?-30軟體尋找油氣藏的方法，以提高三維圈閉評價的精確性。
[0005]為達到上述目的，本發明提供了一種基於I狀?-30軟體尋找油氣藏的方法，包括以下步驟:
[0006]根據三維地震資料和測井資料建立勘探工區的三維巖性及斷層數據立方體；
[0007]將所述三維巖性及斷層數據立方體劃分成若干等深度的深度切片，並分別對每個深度切片進行獨立砂體單元劃分；
[0008]將所述三維巖性及斷層數據立方體的深度切片順序輸入I狀?-30軟體進行油氣藏評價，其具體包括:
[0009]逐層對每個深度切片內各個獨立砂體單元分別進行圈閉評價，獲取每個深度切片內的各個獨立砂體單元的圈閉量估值；
[0010]根據所述每個深度切片內各個獨立砂體單元的圈閉量估值繪出平面上的甜點圖，以及展示縱向不同深度上的油氣圈閉量，並獲取所述勘探工區內儲藏油氣的總圈閉量估值。
[0011]本發明的評價方法，所述根據三維地震資料和測井資料建立勘探工區的三維巖性及斷層數據立方體，具體包括:
[0012]將從三維地震資料中獲得的構造信息和斷層信息，與從測井資料中獲得的巖性信息相結合建立勘探工區的三維巖性及斷層數據立方體。
[0013]本發明的評價方法，所述分別對每個深度切片進行獨立砂體單元劃分，具體為:
[0014](1)、選擇一個深度切片並將其內的巖性信息和斷層信息分別用相應的特徵參數表不；
[0015](2)、按照預設的連通砂點搜索規則搜索出該深度切片內所有的由若干相互連通的砂點構成的砂點組，其中所述每個砂點組的邊界由泥巖和斷層界定；
[0016](3)、將所述每個砂點組內所有砂點的特徵參數合併，以作為一個獨立砂體單元，並對每個獨立砂體單元分別進行唯一性標識；
[0017](4)、重複上述步驟(1)?(3)，直至逐層順序完成所述三維巖性及斷層數據立方體的所有深度切片的獨立砂體單元劃分。
[0018]本發明的評價方法，所述預設的連通砂點搜索規則包括:
[0019]對於與斷層不相鄰的砂點，按照八方向搜索與該砂點相連通的相鄰砂點；以及，
[0020]對於與斷層相鄰的砂點，按照四方向搜索與該砂點相連通的相鄰砂點。
[0021]本發明的評價方法，在所述逐層對每個深度切片內的各個獨立砂體單元分別進行圈閉評價之前還包括:
[0022]為每個深度切片內的各個獨立砂體單元設置圈閉評價初值。
[0023]本發明的評價方法，在所述為每個深度切片內的各個獨立砂體單元設置圈閉評價初值的過程中，包括:
[0024]如果當前深度切片的某個獨立砂體單元內的砂點與其上一層深度切片的某個獨立砂體單元內的砂點存在上下連通關係，則所述當前深度切片的某個獨立砂體單元直接繼承所述上一層深度切片的某個獨立砂體單元的圈閉評價初值。
[0025]本發明的評價方法，所述逐層對每個深度切片內各個獨立砂體單元分別進行圈閉評價，獲取每個深度切片內的各個獨立砂體單元的圈閉量估值，具體包括:
[0026](1)、選擇一個深度切片，根據圈閉評價更新規則和該深度切片裡的斷層兩側的巖性的對置情況，分別對該深度切片內的各個獨立砂體單元的圈閉評價初值進行更新，以獲得該深度切片內的各個獨立砂體單元的最終圈閉評價值，並將對應的最終圈閉評價值賦予對應獨立砂體單元內的每個砂點；
[0027](2)、根據公式3? =獲取深度切片的圈閉量估值，其中，3?為第』個深度切片的圈閉量估值，為第』個深度切片的第1個獨立砂體單元的空隙容積，
為第』個深度切片的第1個獨立砂體單元的最終圈閉評價值；
[0028](3)、重複上述步驟(1)?(2)直至逐層順序完成所述三維巖性及斷層數據立方體的所有深度切片的圈閉量估值的獲取。
[0029]本發明的評價方法，所述圈閉評價更新規則包括:
[0030]若斷層對側為泥巖，則根據每個深度切片內的各個獨立砂體單元的泥巖的塗抹封堵效應的不同對應調整所述各個獨立砂體單元的圈閉評價值；
[0031]若斷層對側為砂巖，則當每個深度切片內的斷層一側的獨立砂體單元的圈閉評價值高於該斷層對側的獨立砂體單元的圈閉評價值時，按下式計算圈閉評價漏失值:
[0032]= 00馭氺(^2?) /82^1/4
[0033]其中，I私1(為圈閉評價漏失值，00,為斷層兩側獨立砂體單元的圈閉評價值之差，1X2！！為斷層兩側獨立砂體單元的圈閉評價值巖性指數之和，8^1為漏失係數。
[0034]本發明首先把三維地震資料所獲得的構造情況與測井資料獲得的地下巖性的情況結合起來，加上斷層的位置信息，綜合成一個數據立方體；其次將數據立方體進行等深度劃分，劃分出若干等深度的深度切片，然後依次對每個深度切片進行圈閉評價，根據所述每個深度切片內各個獨立砂體單元的圈閉量估值繪出平面上的甜點圖，以及展示縱向不同深度上的油氣圈閉量，並獲取所述勘探工區內儲藏油氣的總圈閉量估值，從而可以將工區儲集油氣的各種圈閉都找出來。由於本發明將數據立方體劃分成若干等深度的深度切片，並依次對每個深度切片進行圈閉評價，從而通過這種更為精細的圈閉評價提高了三維圈閉評價的精確性，有利於精確尋找油氣藏，並且本發明可繪出平面上的甜點圖，以及展示縱向不同深度上的油氣圈閉量，並可得到勘探工區內油氣藏的總圈閉量。
【專利附圖】
【附圖說明】
[0035]此處所說明的附圖用來提供對本發明的進一步理解，構成本申請的一部分，並不構成對本發明的限定。在附圖中:
[0036]圖1為本發明一個實施例的基於1狀？-30軟體尋找油氣藏的方法的流程圖；
[0037]圖2為本發明一個實施例中三維巖性及斷層數據立方體的示意圖；
[0038]圖3為本發明一個實施例中一個深度切片的砂泥巖巖性劃分及斷層分布示意圖；
[0039]圖4為本發明一個實施例中一個深度切片的砂巖聯通單元劃分圖；
[0040]圖5為本發明一個實施例的平面甜點圖；
[0041]圖6為本發明一個實施例中自上而下的各深度切片內各砂體單元圈閉評價方法與評價傳遞方法示意圖；
[0042]圖7為本發明一個實施例中斷層漏失量計算方法示意圖；
[0043]圖8為本發明一個實施例中經I狀?-30軟體處理前後兩個相鄰切片的圈閉評價變化示意圖；
[0044]圖9為本發明一個實施例中經I狀?-30軟體處理前後另外兩個相鄰垂直剖面的圈閉評價變化示意圖；
[0045]圖10為本發明一個實施例中深度域完成圈閉評價後輸出的5條曲線示意圖。
[0046]圖1匕為本發明一個實施例中某工區均勻打採油井12 口的打井前儲量分布示意圖；
[0047]圖116為本發明一個實施例中某工區均勻打採油井12 口的打井後，油氣採完後的剩餘儲量分布示意圖；
[0048]圖123為本發明一個實施例中例油田的三維巖性數據立方體的一個切片的自然伽馬值圖；
[0049]圖126和圖133均為根據圖12&中自然那伽馬值推算的的砂泥巖巖性轉換示意圖；
[0050]圖136為在圖133中插入斷層後的切片示意圖；
[0051]圖14為本發明一個實施例中例油田的一個經過I狀?-30軟體分析後的垂直剖面的圈閉分析剖面圖；
[0052]圖153為本發明一個實施例中例油田按常規斷塊油田常規方法冒估的含油範圍示意圖；
[0053]圖156為本發明一個實施例中例油田經I狀?-30軟體處理後的含油範圍示意圖；
[0054]圖16為本發明一個實施例中IX油田的深度域評價結果輸出的5條統計曲線示意圖。
【具體實施方式】
[0055]為使本發明的目的、技術方案和優點更加清楚明白，下面結合實施例和附圖，對本發明做進一步詳細說明。在此，本發明的示意性實施例及其說明用於解釋本發明，但並不作為對本發明的限定。
[0056]下面結合附圖，對本發明的【具體實施方式】作進一步的詳細說明。
[0057]參考圖1所示，本發明實施例的基於1狀？-30的油氣藏評價方法包括以下步驟:
[0058]步驟3101、根據三維地震資料和測井資料建立勘探工區的三維巖性及斷層數據立方體。即將從三維地震資料中獲得的構造信息和斷層信息，與從測井資料中獲得的巖性信息相結合建立勘探工區的三維巖性及斷層數據立方體。其中，三維巖性及斷層數據立方體由按孔隙度大小區分的泥巖、粉砂巖、中等砂巖和好砂巖四種巖性以及斷層面在三維等間距網格剖分空間中的展布信息組成，斷層應為在三維空間等間距剖分網格中逐點連續展布的緻密(無空缺點)斷層面，例如相應的特徵參數可以為:泥巖用0表示，粉砂巖用1表示，中等砂巖用2表示，好砂巖用3表示，斷層用4表示，要求每個單元具有一個唯一的識別標識(10編號)，如圖4所示。具體來說，根據項目需要，在已建立的構造模型的基礎上，結合測井曲線((?)及聲波曲線等對工區砕削巖進行巖性分類。解釋的結果，以井間插值的方式進行了構造模型的巖性充填，並將沉積相解釋結果作為邊界條件對該模型進行了相控約束。用測井資料加上地震資料進行控制內插，將工區內按自然伽馬值對砂泥巖巖性進行分類。根據60值的大小，把它轉換為好砂巖3，中等砂巖2，粉砂巖1及泥巖0四種巖性。最終，建好的三維巖性及斷層數據立方體，如圖2所示。其中，I'狀?-30即三維圈閉。
[0059]步驟3102、將三維巖性及斷層數據立方體劃分成若干等深度的深度切片，並分別對每個深度切片進行獨立砂體單元劃分。
[0060]如圖3所示，為三維巖性及斷層數據立方體的一個深度切片，其橫坐標為水平方向X，縱坐標為水平方向V。圖中白色為好砂巖，特徵參數為3 ；淺灰色為中等砂巖，特徵參數為2 ；深灰色為粉砂巖，特徵參數為1 ；淺黑色為泥巖，特徵參數為0 ；黑色為斷層，特徵參數為4。其中，分別對每個深度切片進行獨立砂體單元劃分，具體為:
[0061](1)、選擇一個深度切片並將其內的巖性信息和斷層信息分別用相應的特徵參數表不。
[0062](2)、按照預設的連通砂點搜索規則搜索出該深度切片內所有的由若干相互連通的砂點構成的砂點組，其中每個砂點組的邊界由泥巖和斷層界定；該預設的連通砂點搜索規則為:
[0063]對於與斷層不相鄰的砂點，按照八方向(例如0度、45度、90度、135度、180度、225度、270度和315度)搜索與該砂點相連通的相鄰砂點；以及，
[0064]對於與斷層相鄰的砂點，按照四方向(例如0度、90度、180度和270度)搜索與該砂點相連通的相鄰砂點。
[0065](3)、將每個砂點組內所有砂點的特徵參數合併，以作為一個獨立砂體單元，並對每個獨立砂體單元分別進行唯一性標識。
[0066](4)、重複上述步驟(1)?(3)，直至逐層順序完成三維巖性及斷層數據立方體的所有深度切片的獨立砂體單元劃分。
[0067]舉例來說，假設劃分出的若干等深度的深度切片如圖6所示。
[0068](^)切片[為特殊情況(因為它是第一張切片，上面沒有評價控制):圖中若有兩個砂點，由於不知上方是否為泥巖，所以就給以最低評價IV = 16 ；
[0069](^)在切片吧裡，泥巖裡出現一個小砂點，它被上方附的泥巖所覆蓋，因此，就給以最高評價靈=800 ；此外，令四周邊水侵入。
[0070](0)在切片吧中，小砂點變成大砂餅，它們與上方是相連的，於是將上方吧中砂點的評價IV = 800傳遞給下它。
[0071]((1)在切片財中，大砂餅繼續增大，但與邊界相連了，於是邊水侵入，圈閉破壞，評價降為最小，即靈=16。如果它與右方含水砂巖相連，IV也會降為最小。
[0072]步驟31033、逐層對每個深度切片內各個獨立砂體單元分別進行圈閉評價，獲取每個深度切片內的各個獨立砂體單元的圈閉量估值。其中，逐層對每個深度切片內各個獨立砂體單元分別進行圈閉評價，獲取每個深度切片內的各個獨立砂體單元的圈閉量估值，具體包括:
[0073](1)、選擇一個深度切片，根據圈閉評價更新規則和該深度切片裡的斷層兩側的巖性的對置情況，分別對該深度切片內的各個獨立砂體單元的圈閉評價初值進行更新，以獲得該深度切片內的各個獨立砂體單元的最終圈閉評價值，並將對應的最終圈閉評價值賦予對應獨立砂體單元內的每個砂點；其中，圈閉評價更新規則為:
[0074]若斷層對側為泥巖，則根據每個深度切片內的各個獨立砂體單元的泥巖的塗抹封堵效應的不同對應調整各個獨立砂體單元的圈閉評價值；
[0075]若斷層對側為砂巖，則當每個深度切片內的斷層一側的獨立砂體單元的圈閉評價值高於該斷層對側的獨立砂體單元的圈閉評價值時，按下式計算圈閉評價漏失值:
[0076]= 00馭氺(^2?) /82^1/4
[0077]其中[拉漢為圈閉評價漏失值，0017為斷層兩側獨立砂體單元的圈閉評價值之差，1X2！！為斷層兩側獨立砂體單元的圈閉評價值巖性指數之和，3^1為漏失係數，8^1由用戶指定，若32從=5就意味遇到對方5顆砂，才漏到和對方一樣。
[0078]舉例來說，例如圖7所示，101好砂層在八80三處遇斷層，0處斷層對面是泥巖，8處斷層對面是粉砂巖，八處對面是好砂巖，於是可以按照上述圈閉評價更新規則進行更新。
[0079](2)、根據公式3? =獲取深度切片的圈閉量估值，其中，3?為第』個深度切片的圈閉量估值，為第』個深度切片的第1個獨立砂體單元的空隙容積，
為第』個深度切片的第1個獨立砂體單元的最終圈閉評價值。
[0080](3)、重複上述步驟(1)?(2)直至逐層順序完成三維巖性及斷層數據立方體的所有深度切片的圈閉量估值的獲取。
[0081]步驟31031根據每個深度切片內各個獨立砂體單元的圈閉量估值繪出平面上的甜點圖，以及展示縱向不同深度上的油氣圈閉量，並獲取勘探工區內儲藏油氣的總圈閉量估值。各個深度切片的圈閉量估值在深度2方向進行累積後繪出平面的總圈閉量，如圖5所示，它反映著工區平面裡每個點位處地下累計的圈閉量的大小，從該圖可以看到，平面上哪裡圈閉量最多，其含油氣的可能性就大。就可以供布探井做抉擇。
[0082]上述步驟3103^510313即為三維巖性及斷層數據立方體的深度切片順序輸入丁狀?-30軟體進行油氣藏評價過程。此外，本發明實施例中，在上述逐層對每個深度切片內的各個獨立砂體單元分別進行圈閉評價之前還包括:
[0083]為每個深度切片內的各個獨立砂體單元設置圈閉評價初值。並且，在為每個深度切片內的各個獨立砂體單元設置圈閉評價初值的過程中，還包括:
[0084]如果當前深度切片的某個獨立砂體單元內的砂點與其上一層深度切片的某個獨立砂體單元內的砂點存在上下連通關係，則當前深度切片的某個獨立砂體單元直接繼承上一層深度切片的某個獨立砂體單元的圈閉評價初值。
[0085]本發明實施例中，如圖8所示，假設兩個相鄰的切片41與42，圖8左邊分別是相鄰的切片41與42的巖性指數的砂泥巖分布的顯示，圖8右邊是經過I狀?-30軟體進行圈閉分析後，切片41與42的圈閉評價值圖。由右邊兩圖可見由於邊水滲入，也由於斷層兩側砂巖見面，而使圈閉評價值逐步下降。軌跡如箭頭所示，白色箭頭表示邊水通過斷層聯通而入侵，使圈閉評價降低的路徑。
[0086]切片42的封閉含油範圍比切片41更有所減少。
[0087]本發明實施例中，如圖9所示，假設兩個相鄰的切片45與46為複雜巖性分布的複雜斷塊兩條相鄰的切片。圖9左邊分別是相鄰的切片45與46的巖性指數仏的砂泥巖分布的顯示，圖9右邊是經過I狀?-30軟體進行圈閉分析後，切片45與45的圈閉評價值圖，由此可以看到各砂層中的油水分布的情形(淺亮色是水其中，白色為圈閉評價極低的含水區，箭頭所示處為油水界面。
[0088]此外，本發明實施例中，還可以在深度域繪出各深度切片的砂泥巖百分比及每個深度切片的圈閉量的對比圖，如圖10所示。這可更直觀的看到在哪個深度上，哪個地層中圈閉體積最大。
[0089]本發明實施例中，如圖1匕和圖116所示，為打井12 口前後的儲量分布及剩餘未動用儲量分布的示意圖，在上述複雜斷塊模型的工區中均勻打井12 口，每口把所遇油層採幹，結果還留下不能採到的「死油區」以及12 口井未鑽遇的「剩餘油」，從圖116可以看到還有不少未動用儲量可供開發。這些未動用儲量也可以由I狀?-30軟體計算出來。
[0090]下面以例油田為例，對本發明實施例進行舉例說明:
[0091]例油田是一個複雜斷塊油田，本次研究面積約54平方千米，含油層系深度約1-2千米。工區內已打井8 口。見少量薄油層。本次切片分析的深度為1300111-1850111。
[0092]為了提供I狀?-30軟體所需的輸入數據格式，需要準備一個三維巖性及斷層數據立方體。我們在8 口井測井數據的井控條件下作砂層內插，並加上地震解釋的斷層數據，形成三維巖性及斷層數據立方體。假設三維巖性及斷層數據立方體中每個單元點的大小為18111*18111*1111，深度方向按每1111為一片切片進行分析。
[0093](1)根據項目需要，在已建立的構造模型的基礎上，結合測井曲線(⑶)及聲波曲線等對工區砕削巖進行巖性分類。解釋的結果，以井間插值的方式進行了構造模型的巖性充填，並將沉積相解釋結果作為邊界條件對該模型進行了相控約束。用測井資料加上地震資料進行控制內插，將工區內按自然伽馬值對砂泥巖巖性進行分類。如圖123所示，為例油田的某切片自然伽馬⑶值彩色圖。根據⑶值的大小，把它轉換為好砂巖3，砂巖2，粉砂巖1及泥巖0，四種巖性，如圖126和133圖所示，為自然伽馬值的轉換圖。
[0094]下一步再根據地震資料解釋的斷層數據融合進這些切片裡去。
[0095]用通常的地震資料解釋方法，對三維地震數據作過細的斷層解釋，從解釋過的斷層數據中導出斷層文件，再將每個斷層點見坐標插進解釋好的三維巖性及斷層數據立方體。給於每個斷層點標上巖性指數為4。如圖136所示，是砂泥巖劃分的四種分布，右邊是插入斷層(黑色細線)後的情況。此圖符合I'狀?-30程序的輸入要求。
[0096]當每個切片都做完上述工作後，就完成了整個三維巖性及斷層數據立方體的數據準備，就完成了適合於I'狀?-30軟體分析的輸入數據要求。
[0097](11)1^油田實際數據I'狀?-30試算結果
[0098]如圖14所示，對於一條經過I'狀?-30軟體處理的深度剖面，可看到砂層雖然很發育，但是只有少數地方具有圈閉。圖14中右下方幾條長條形黑線是斷層，其它深黑色的是巖性圈閉，圖14左方有一個背斜構造圈閉，如箭頭所示。其中，灰色部分是泥巖，圖中除下方兩根垂向黑線是斷層之外，其它黑色部分均為含油氣的構造圈閉及巖性圈閉，白色為含水砂巖。
[0099]圖153是例油田按斷塊油田常規方法冒估的含油範圍，參考圖156所示，例油田經I'狀?-30軟體處理後的含油範圍與上述冒估的含油範圍基本吻合。由此可見，本申請的I'狀?-30軟體更具有定量的概念，它包含了不少巖性圈閉，及北部的小背斜構造圈閉。
[0100]如圖16所示，軟體計算的按深度統計的各等深度切片中的砂泥巖百分比(曲線1及2)，根據311111(111(16X1計算的每個切片的空隙容積(曲線3)，表達儲集性能的3皿|1」，(曲線4),以及最右邊的曲線5表達著此深度切片中的圈閉量33 =811111 圖16中，雖然深度16300處砂巖百分比高達60 %，但是圈閉量卻不大，根據曲線4及5可看出深度1420-15200這一段砂巖雖然不那麼多，但是圈閉量卻較大，實際開採也證明這正是本油田的主要產層層段。可見I狀?-30軟體給出了較好的定量的結果。
[0101](111)最後例油田的圈閉儲量計算結果如下:
[0102]1)根據資料:平均孔隙度為28.5% ；平均含油飽和度為70% ；砂巖百分比為47.3% ；泥巖百分比為52.7% ；中部主體圈閉面積約10.5如12。
[0103]2)研究範圍總圈閉儲量為333，(即各切片全部加起來)。
[0104]888 = 2】X平均孔隙度^立方單兀體積/800
[0105]I'狀?-30軟體計算結果:研究範圍總圈閉儲量約為0.53父108！:！3。
[0106]本發明實施例首先把三維地震資料所獲得的構造情況與測井資料獲得的地下巖性的情況結合起來，加上斷層的位置信息，綜合成一個數據立方體；其次將數據立方體進行等深度劃分，劃分出若干等深度的深度切片，然後依次對每個深度切片進行圈閉評價，從而可以將工區儲集油氣的各種圈閉都找出來。由於本發明實施例將數據立方體劃分成若干等深度的深度切片，並依次對每個深度切片進行圈閉評價，從而通過這種更為精細的圈閉評價提高了三維圈閉評價的精確性。
[0107]以上所述的具體實施例，對本發明的目的、技術方案和有益效果進行了進一步詳細說明，所應理解的是，以上所述僅為本發明的具體實施例而已，並不用於限定本發明的保護範圍，凡在本發明的精神和原則之內，所做的任何修改、等同替換、改進等，均應包含在本發明的保護範圍之內。
【權利要求】
1.一種基於TRAP-3D軟體尋找油氣藏的方法，其特徵在於，包括以下步驟: 根據三維地震資料和測井資料建立勘探工區的三維巖性及斷層數據立方體； 將所述三維巖性及斷層數據立方體劃分成若干等深度的深度切片，並分別對每個深度切片進行獨立砂體單元劃分； 將所述三維巖性及斷層數據立方體的深度切片順序輸入TRAP-3D軟體進行油氣藏評價，其具體包括: 逐層對每個深度切片內各個獨立砂體單元分別進行圈閉評價，獲取每個深度切片內的各個獨立砂體單元的圈閉量估值； 根據所述每個深度切片內各個獨立砂體單元的圈閉量估值繪出平面上的甜點圖，以及展示縱向不同深度上的油氣圈閉量，並獲取所述勘探工區內儲藏油氣的總圈閉量估值。
2.根據權利要求1所述的方法，其特徵在於，所述根據三維地震資料和測井資料建立勘探工區的三維巖性及斷層數據立方體，具體包括: 將從三維地震資料中獲得的構造信息和斷層信息，與從測井資料中獲得的巖性信息相結合建立勘探工區的三維巖性及斷層數據立方體。
3.根據權利要求1所述的方法，其特徵在於，所述分別對每個深度切片進行獨立砂體單元劃分，具體為: (1)、選擇一個深度切片並將其內的巖性信息和斷層信息分別用相應的特徵參數表示； (2)、按照預設的連通砂點搜索規則搜索出該深度切片內所有的由若干相互連通的砂點構成的砂點組，其中所述每個砂點組的邊界由泥巖和斷層界定； (3)、將所述每個砂點組內所有砂點的特徵參數合併，以作為一個獨立砂體單元，並對每個獨立砂體單元分別進行唯一性標識； (4)、重複上述步驟(I)?(3)，直至逐層順序完成所述三維巖性及斷層數據立方體的所有深度切片的獨立砂體單元劃分。
4.根據權利要求3所述的方法，其特徵在於，所述預設的連通砂點搜索規則包括: 對於與斷層不相鄰的砂點，按照八方向搜索與該砂點相連通的相鄰砂點；以及， 對於與斷層相鄰的砂點，按照四方向搜索與該砂點相連通的相鄰砂點。
5.根據權利要求1所述的方法，其特徵在於，在所述逐層對每個深度切片內的各個獨立砂體單元分別進行圈閉評價之前還包括: 為每個深度切片內的各個獨立砂體單元設置圈閉評價初值。
6.根據權利要求5所述的方法，其特徵在於，在所述為每個深度切片內的各個獨立砂體單元設置圈閉評價初值的過程中，包括: 如果當前深度切片的某個獨立砂體單元內的砂點與其上一層深度切片的某個獨立砂體單元內的砂點存在上下連通關係，則所述當前深度切片的某個獨立砂體單元直接繼承所述上一層深度切片的某個獨立砂體單元的圈閉評價初值。
7.根據權利要求5所述的方法，其特徵在於，所述逐層對每個深度切片內各個獨立砂體單元分別進行圈閉評價，獲取每個深度切片內的各個獨立砂體單元的圈閉量估值，具體包括: (I)、選擇一個深度切片，根據圈閉評價更新規則和該深度切片裡的斷層兩側的巖性的對置情況，分別對該深度切片內的各個獨立砂體單元的圈閉評價初值進行更新，以獲得該深度切片內的各個獨立砂體單元的最終圈閉評價值，並將對應的最終圈閉評價值賦予對應獨立砂體單元內的每個砂點； (2)、根據公式SS^= Sumiindexi^MVJ獲取深度切片的圈閉量估值，其中，SS^為第j個深度切片的圈閉量估值，Indexi為第j個深度切片的第i個獨立砂體單元的空隙容積，MVi為第j個深度切片的第i個獨立砂體單元的最終圈閉評價值； (3)、重複上述步驟(I)?(2)直至逐層順序完成所述三維巖性及斷層數據立方體的所有深度切片的圈閉量估值的獲取。
8.根據權利要求7所述的方法，其特徵在於，所述圈閉評價更新規則包括: 若斷層對側為泥巖，則根據每個深度切片內的各個獨立砂體單元的泥巖的塗抹封堵效應的不同對應調整所述各個獨立砂體單元的圈閉評價值； 若斷層對側為砂巖，則當每個深度切片內的斷層一側的獨立砂體單元的圈閉評價值高於該斷層對側的獨立砂體單元的圈閉評價值時，按下式計算圈閉評價漏失值: LEAK = DDMV* (YXZH)/SEAL/4，其中，LEAK為圈閉評價漏失值，DDMV為斷層兩側獨立砂體單元的圈閉評價值之差，YXZH為斷層兩側獨立砂體單元的圈閉評價值巖性指數之和，SEAL為漏失係數。
【文檔編號】G01V1/28GK104375179SQ201410635467
【公開日】2015年2月25日 申請日期:2014年11月5日 優先權日:2014年11月5日
【發明者】李慶忠, 張立彬, 張曉敏 申請人:中國石油天然氣集團公司, 中國石油集團東方地球物理勘探有限責任公司