[發(fā)明專利]將測(cè)井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置

有效

• 申請(qǐng)?zhí)枺?/td>: 201410057529.7
• 申請(qǐng)日: 2014-02-19
• 公開(kāi)（公告）號(hào): CN103886129B
• 公開(kāi)（公告）日: 2017-06-13
• 發(fā)明（設(shè)計(jì)）人: 倪廣元;王劍秦;董肇偉
• 申請(qǐng)（專利權(quán)）人: 北京中恒永信科技有限公司
• 主分類號(hào): G06F17/50
• 分類號(hào): G06F17/50
• 代理公司: 北京路浩知識(shí)產(chǎn)權(quán)代理有限公司11002
• 代理人: 李迪

技術(shù)領(lǐng)域

本發(fā)明涉及地質(zhì)屬性建模技術(shù)領(lǐng)域，具體涉及一種將測(cè)井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置。

背景技術(shù)

地質(zhì)體內(nèi)部含有多種反映地層巖性、油氣資源分布等特性的參數(shù)，如孔隙度、滲透率、含有飽和度、沉積相等，地質(zhì)體物性建模就是利用空間插值法，如Kriging、序貫高斯等，求解地質(zhì)體內(nèi)部物理、化學(xué)屬性參數(shù)模型。利用地質(zhì)體物性建模成果，可以有效的發(fā)現(xiàn)儲(chǔ)層，并對(duì)油氣資源分布情況進(jìn)行準(zhǔn)確預(yù)測(cè)，從而指導(dǎo)油田開(kāi)發(fā)。考慮到地層成因，地質(zhì)體中某一區(qū)域的物性，與同一層物性的關(guān)系遠(yuǎn)遠(yuǎn)大于該區(qū)域的上層或下層。為了更精確的表達(dá)地質(zhì)體物性，實(shí)際建模過(guò)程中，在對(duì)地質(zhì)體細(xì)分小層的基礎(chǔ)上，會(huì)進(jìn)一步細(xì)分，對(duì)細(xì)分后的每一層，利用由測(cè)井曲線數(shù)據(jù)得到的關(guān)于該層的物性信息，使用空間插值算法進(jìn)行物性建模。將測(cè)井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型是利用測(cè)井?dāng)?shù)據(jù)，為井穿過(guò)的儲(chǔ)層網(wǎng)格單元賦值。被賦值的網(wǎng)格單元，在進(jìn)行物性建模時(shí)，作為已知量，進(jìn)而確定其他未被測(cè)量過(guò)區(qū)域的屬性值。可見(jiàn)，它是進(jìn)行地質(zhì)體物性建模的基礎(chǔ)。

目前，角點(diǎn)網(wǎng)格是被廣泛使用的描述儲(chǔ)集層的網(wǎng)格模型。作為結(jié)構(gòu)化網(wǎng)格模型，它具有易于實(shí)現(xiàn)的優(yōu)點(diǎn)。高度的靈活性使得它能夠很好的表征儲(chǔ)集層地質(zhì)特征，比如非均質(zhì)性以及斷層。

然而，由于角點(diǎn)網(wǎng)格不使用笛卡爾坐標(biāo)系，在離散測(cè)井?dāng)?shù)據(jù)時(shí)非常復(fù)雜繁瑣。再加上為了滿足對(duì)建模區(qū)域精細(xì)描述的需求，現(xiàn)在的網(wǎng)格劃分都非常細(xì)密。故而，一種快速地將測(cè)井?dāng)?shù)據(jù)離散到儲(chǔ)層模型的方法被迫切的需要。

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

（一）解決的技術(shù)問(wèn)題

針對(duì)現(xiàn)有技術(shù)的不足，本發(fā)明提供了一種將測(cè)井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置，根據(jù)測(cè)井?dāng)?shù)據(jù)值，用盡可能少的時(shí)間為井軌跡所穿過(guò)的網(wǎng)格賦值，從而減少建模人員的等待時(shí)間，提高工作效率。

（二）技術(shù)方案

為實(shí)現(xiàn)以上目的，本發(fā)明通過(guò)以下技術(shù)方案予以實(shí)現(xiàn)：

一種將測(cè)井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法，所述儲(chǔ)層網(wǎng)格模型為角點(diǎn)網(wǎng)格，包括下列步驟：

S1、獲取井軌跡穿越的網(wǎng)格單元；

S2、計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)，并利用所述交點(diǎn)將所述井軌跡劃分為多個(gè)區(qū)塊；

S3、對(duì)落入每一區(qū)塊的所述井軌跡上的測(cè)井曲線屬性值進(jìn)行加權(quán)平均，求得屬性均值，并將所述屬性均值賦給對(duì)應(yīng)的網(wǎng)格單元。

優(yōu)選的，步驟S1包括：

S11、根據(jù)儲(chǔ)層網(wǎng)格模型的網(wǎng)格厚度信息設(shè)置井軌跡中待測(cè)點(diǎn)選取的間隔；

S12、按照所述間隔順次取出待測(cè)點(diǎn)；

S13、找到每一個(gè)待測(cè)點(diǎn)所在的網(wǎng)格單元；

S14、將所述網(wǎng)格單元進(jìn)行整理與合并。

優(yōu)選的，步驟S11包括：

S111、根據(jù)井頭坐標(biāo)，以及儲(chǔ)層網(wǎng)格模型的邊界信息和網(wǎng)格長(zhǎng)寬設(shè)置信息，估算出位于井頭下方的網(wǎng)格單元；

S112、計(jì)算井頭下方網(wǎng)格單元的平均厚度；

S113、結(jié)合井軌跡數(shù)據(jù)的采樣間隔，計(jì)算出待測(cè)點(diǎn)選取的間隔；

其中，步驟S13包括：利用參考索引找出對(duì)應(yīng)網(wǎng)格單元，以網(wǎng)格單元為中心，采用臨近搜索法選取網(wǎng)格單元，對(duì)待測(cè)點(diǎn)和選取的網(wǎng)格單元進(jìn)行包含性測(cè)試；

其中，步驟S14包括：通過(guò)網(wǎng)格單元的索引信息，利用折半查找法，找出未被測(cè)出的網(wǎng)格單元；對(duì)于多個(gè)待測(cè)點(diǎn)落入同一網(wǎng)格單元的情形，只記錄第一個(gè)落入該網(wǎng)格單元的待測(cè)點(diǎn)信息，與最后一個(gè)落入該網(wǎng)格單元的待測(cè)點(diǎn)信息。

優(yōu)選的，步驟S2中計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)包括步驟：

根據(jù)網(wǎng)格單元索引信息確定單元內(nèi)同井軌跡有交點(diǎn)的面；選取面兩端的待測(cè)點(diǎn)，計(jì)算兩待測(cè)點(diǎn)間的井軌跡與網(wǎng)格單元面的交點(diǎn)。

優(yōu)選的，計(jì)算兩待測(cè)點(diǎn)間的井軌跡與網(wǎng)格單元面的交點(diǎn)時(shí)使用快速逼近法，實(shí)現(xiàn)方式是：

選取兩相鄰井軌跡數(shù)據(jù)點(diǎn)組成直線，計(jì)算同網(wǎng)格單元表面的交點(diǎn)，計(jì)算所述交點(diǎn)與數(shù)據(jù)點(diǎn)的距離；

在井軌跡上選取距離相等的數(shù)據(jù)點(diǎn)，并以此數(shù)據(jù)點(diǎn)為基準(zhǔn)，選取與之相鄰的下一個(gè)井軌跡數(shù)據(jù)點(diǎn)，利用重新得到的兩相鄰井軌跡數(shù)據(jù)點(diǎn)，進(jìn)行下一次的求交點(diǎn)運(yùn)算；

當(dāng)交點(diǎn)位于選取的兩相鄰井軌跡數(shù)據(jù)點(diǎn)之間時(shí)，結(jié)束運(yùn)算。

優(yōu)選的，步驟S2中將井軌跡劃分為多個(gè)區(qū)塊時(shí)，將兩個(gè)相鄰的井與網(wǎng)格單元面的交點(diǎn)之間的區(qū)域，作為一個(gè)區(qū)塊，每一區(qū)塊對(duì)應(yīng)儲(chǔ)層網(wǎng)格模型的一個(gè)網(wǎng)格單元。

優(yōu)選的，步驟S3中計(jì)算落入每一區(qū)塊的測(cè)井曲線屬性值均值時(shí)，得到唯一的數(shù)值，并以此數(shù)值作為區(qū)塊對(duì)應(yīng)網(wǎng)格單元的數(shù)值。

本發(fā)明還提供了一種將測(cè)井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的裝置，包括：