本發(fā)明涉及盆地模擬和油氣成藏過(guò)程模擬技術(shù)領(lǐng)域，特別是勘探程度較高的盆地級(jí)、凹陷級(jí)不同層次油氣系統(tǒng)的一種基于角點(diǎn)網(wǎng)格的高精度油氣運(yùn)移模擬方法。

背景技術(shù)：

隨著油氣藏勘探開發(fā)的不斷深入，對(duì)基于高精度地質(zhì)模型下的油氣成藏模擬提出了更高的要求。傳統(tǒng)的直角正交網(wǎng)格模型在油氣數(shù)值模擬過(guò)程中應(yīng)用較為廣泛，其特點(diǎn)是計(jì)算速度快，其所有網(wǎng)格的長(zhǎng)寬相等，垂向連接網(wǎng)格面垂直。但是在有斷層情況下，這種方式就不能建立精細(xì)的地質(zhì)模型。角點(diǎn)網(wǎng)格模型克服了正交網(wǎng)格模型的缺點(diǎn)，可以根據(jù)控制區(qū)域的長(zhǎng)度和寬度靈活調(diào)整各油層的網(wǎng)格步長(zhǎng)，從而使斷層附近的網(wǎng)格邊線與斷層的走向重合。該技術(shù)能夠更精準(zhǔn)的完成油氣運(yùn)移模擬，為油藏開發(fā)提供更好的分析工具。更精準(zhǔn)的角點(diǎn)網(wǎng)格模型帶來(lái)了更精確的模擬刻畫，同時(shí)也模擬算法實(shí)現(xiàn)帶來(lái)了更大的難度。目前的各類模擬軟件，僅能夠?qū)崿F(xiàn)傳統(tǒng)正交網(wǎng)格模型的模擬方法，如何在更精準(zhǔn)的角點(diǎn)網(wǎng)格模型下，處理好構(gòu)造作用、沉積作用、地?zé)嶙饔?、生烴作用、排烴作用及運(yùn)烴作用、聚烴作用之間的關(guān)系，將該技術(shù)運(yùn)用到油藏開發(fā)中，成為了目前亟待解決的問(wèn)題。

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

本發(fā)明的目的是為了解決以上技術(shù)問(wèn)題，提供一種基于角點(diǎn)網(wǎng)格的高精度油氣運(yùn)移模擬方法。

本發(fā)明的目的可通過(guò)如下技術(shù)措施來(lái)實(shí)現(xiàn)：

一種基于角點(diǎn)網(wǎng)格的高精度油氣運(yùn)移模擬方法，包括：

步驟1，建立基于角點(diǎn)網(wǎng)格的地質(zhì)模型和屬性模型，進(jìn)行前期的熱史、生排烴模擬，形成模擬成果數(shù)據(jù)；

步驟2，首先計(jì)算每個(gè)單元體的中心坐標(biāo)，加密體處理；其次是進(jìn)行運(yùn)移量分配評(píng)估；第三是進(jìn)行油氣運(yùn)移計(jì)算；第四是進(jìn)行初次運(yùn)移、查找通道，將當(dāng)前通道的油氣運(yùn)移達(dá)到體平衡；

其中：在步驟1中，對(duì)于地質(zhì)模型和屬性模型，對(duì)于研究工區(qū)的每套對(duì)應(yīng)地層、斷層都要建立，網(wǎng)格數(shù)統(tǒng)一；斷層要根據(jù)實(shí)際情況進(jìn)行封堵性設(shè)置，屬性模型信息至少應(yīng)包括巖性、孔隙度、滲透率；在步驟2中，對(duì)于局部加密的每個(gè)加密體，每個(gè)小網(wǎng)格要進(jìn)行滲透率、孔隙度處理，然后統(tǒng)計(jì)到大網(wǎng)格三個(gè)方向的滲透率和綜合孔隙度。

上述方案進(jìn)一步包括：

進(jìn)行前期的熱史、生排烴模擬，要形成標(biāo)準(zhǔn)格式的溫度文件、超壓文件、生排烴文件、油氣流體勢(shì)文件。

在步驟2中，通過(guò)達(dá)西定律進(jìn)行運(yùn)移量分配評(píng)估，計(jì)算每個(gè)單元格的油粘度、油密度、液體深度壓力、XYZ三個(gè)方向的流體速度。

在步驟2中，進(jìn)行油氣運(yùn)移計(jì)算，根據(jù)地質(zhì)體屬性計(jì)算每個(gè)網(wǎng)格能容納的烴量，計(jì)算每個(gè)網(wǎng)格向相鄰網(wǎng)格的運(yùn)移情況。

在步驟2中，單元格間運(yùn)移，要對(duì)于每個(gè)單元格的各個(gè)方向，根據(jù)巖性計(jì)算網(wǎng)格間接觸面積、排替壓力，進(jìn)行斷層封堵性判斷。

本發(fā)明中的基于角點(diǎn)網(wǎng)格的高精度油氣運(yùn)移模擬方法，采用了能夠較為準(zhǔn)確描述裂縫空間展布和特征的角點(diǎn)網(wǎng)格模型，同時(shí)運(yùn)用了動(dòng)態(tài)的體平衡技術(shù)，從而提高了運(yùn)移模擬精度，為勘探開發(fā)提供有效依據(jù)。

附圖說(shuō)明

圖1為本發(fā)明的一種具體實(shí)施例的流程簡(jiǎn)圖；

圖2是本發(fā)明的角點(diǎn)網(wǎng)格模型的應(yīng)用效果圖；

圖3是本發(fā)明的油氣運(yùn)移模擬應(yīng)用效果圖。

具體實(shí)施方式

為使本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂，通過(guò)典型實(shí)施例，作詳細(xì)說(shuō)明如下。

基于角點(diǎn)網(wǎng)格的高精度油氣運(yùn)移模擬方法包括：

步驟1，建立基于角點(diǎn)網(wǎng)格的地質(zhì)模型和屬性模型，進(jìn)行運(yùn)移前的熱史、生烴、排烴等模擬，逐期次獲取每個(gè)期次的各類成果數(shù)據(jù)（包括體數(shù)據(jù)、斷層信息、生排烴信息、溫度信息、超壓信息、油氣流體勢(shì)等）和屬性數(shù)據(jù)（包括巖性、孔隙度、滲透率等）；步驟2，對(duì)每個(gè)期次進(jìn)行油運(yùn)聚模擬、加密體處理、運(yùn)移量分配、單元格間運(yùn)移等運(yùn)算處理，最后輸出油氣運(yùn)移強(qiáng)度、運(yùn)移方向和運(yùn)移量。

在步驟1中，要首先建立基于角點(diǎn)網(wǎng)格的地質(zhì)模型和屬性模型，對(duì)于地質(zhì)模型和屬性模型，對(duì)于研究工區(qū)的每套對(duì)應(yīng)地層、斷層都要建立，網(wǎng)格數(shù)統(tǒng)一，局部可進(jìn)行加密處理。斷層要根據(jù)實(shí)際情況進(jìn)行封堵性設(shè)置，屬性信息至少應(yīng)包括巖性、孔隙度、滲透率等。對(duì)建好的地質(zhì)模型和屬性模型，進(jìn)行前期的熱史、生排烴模擬，形成標(biāo)準(zhǔn)格式的溫度文件、超壓文件、生排烴文件、油氣流體勢(shì)文件，以便于在油氣運(yùn)移模擬的屬性提取。

在步驟2中，首先計(jì)算每個(gè)單元體的中心坐標(biāo)，對(duì)每個(gè)加密體的小網(wǎng)格處理，然后統(tǒng)計(jì)到大網(wǎng)格三個(gè)方向的滲透率和綜合孔隙度；其次是通過(guò)達(dá)西定律進(jìn)行運(yùn)移量分配評(píng)估，計(jì)算每個(gè)單元格的油粘度、油密度、液體深度壓力、XYZ三個(gè)方向的流體速度；第三是進(jìn)行油氣運(yùn)移，根據(jù)地質(zhì)體屬性計(jì)算每個(gè)網(wǎng)格能容納的烴量，計(jì)算每個(gè)網(wǎng)格可以向相鄰網(wǎng)格的運(yùn)移情況；第四是進(jìn)行初次運(yùn)移、查找通道，將當(dāng)前通道的油氣運(yùn)移達(dá)到體平衡。

下面結(jié)合附圖1對(duì)做進(jìn)一步說(shuō)明。

在步驟101，通過(guò)收集建模所需各項(xiàng)基礎(chǔ)數(shù)據(jù)，應(yīng)用模型建立工具，建立基于角點(diǎn)網(wǎng)格的地質(zhì)模型和屬性模型。流程進(jìn)入到步驟102。

在步驟102，通過(guò)熱史、生排烴方法等，進(jìn)行運(yùn)移模擬前的模擬處理，形成各期次的成果文件。流程進(jìn)入到步驟103。

在步驟103，根據(jù)期次循環(huán)處理以下各步驟。獲取每個(gè)期次的各類成果數(shù)據(jù)和屬性數(shù)據(jù)。成果數(shù)據(jù)包括體數(shù)據(jù)、斷層信息、生排烴信息、溫度信息、超壓信息、油氣流體勢(shì)等；屬性數(shù)據(jù)包括巖性、孔隙度、滲透率等。流程進(jìn)入到步驟104。

在步驟104，對(duì)于有局部加密的數(shù)據(jù)體，要進(jìn)行加密體處理。對(duì)每個(gè)加密體，獲取加密網(wǎng)格邏輯單元起止序號(hào)，計(jì)算XYZ方向被加密網(wǎng)格數(shù)，提取加密體屬性（巖性、孔隙度、滲透率等），對(duì)每個(gè)加密網(wǎng)格進(jìn)行處理，統(tǒng)計(jì)到大網(wǎng)格三個(gè)方向的滲透率及綜合孔隙度。流程進(jìn)入到步驟105。

在步驟105，進(jìn)行運(yùn)移量分配估算，對(duì)每個(gè)單元格分別計(jì)算油粘度、油密度、流體深度壓力、XYZ方向的流體速度。流程進(jìn)入到步驟106。

在步驟106，計(jì)算油氣運(yùn)移。根據(jù)地質(zhì)體屬性計(jì)算網(wǎng)格能容納的烴量，計(jì)算每個(gè)網(wǎng)格向相鄰單元格的運(yùn)移情況，然后進(jìn)行初次運(yùn)移、查找通道，最終將當(dāng)前通道中的油氣運(yùn)移到體平衡，輸出該期次的油氣運(yùn)移強(qiáng)度、運(yùn)移量、運(yùn)移方向文件等。

其中，在單元格間運(yùn)移時(shí)，首先獲取角點(diǎn)網(wǎng)格和斷層文件，其次依次計(jì)算向上和運(yùn)移的其他方向，根據(jù)巖性，計(jì)算網(wǎng)格間接觸面積、排替壓力，判斷斷層封堵性，得到單元格間的運(yùn)移關(guān)系。

圖2是車西地區(qū)基于角點(diǎn)網(wǎng)格的現(xiàn)階段三維地質(zhì)模型效果圖；

圖3是車西地區(qū)應(yīng)用該運(yùn)移模擬方法形成的Ng沉積末時(shí)期運(yùn)移后各地層疊加的油氣分布強(qiáng)度效果圖。

流程結(jié)束。