將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置制造方法
【專利摘要】本發(fā)明提供一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置，涉及地質(zhì)屬性建?！炯夹g(shù)領(lǐng)域】。所述儲(chǔ)層網(wǎng)格模型為角點(diǎn)網(wǎng)格，該方法包括下列步驟：S1、獲取井軌跡穿越的網(wǎng)格單元；S2、計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)，并利用所述交點(diǎn)將所述井軌跡劃分為多個(gè)區(qū)塊；S3、對(duì)落入每一區(qū)塊的所述井軌跡上的測井曲線屬性值進(jìn)行加權(quán)平均，求得屬性均值，并將所述屬性均值賦給對(duì)應(yīng)的網(wǎng)格單元。本發(fā)明根據(jù)測井?dāng)?shù)據(jù)值，用盡可能少的時(shí)間為井軌跡所穿過的網(wǎng)格賦值，從而減少建模人員的等待時(shí)間，提高了工作效率。
【專利說明】將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地質(zhì)屬性建?！炯夹g(shù)領(lǐng)域】，具體涉及一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置。
【背景技術(shù)】
[0002]地質(zhì)體內(nèi)部含有多種反映地層巖性、油氣資源分布等特性的參數(shù)，如孔隙度、滲透率、含有飽和度、沉積相等，地質(zhì)體物性建模就是利用空間插值法，如Kriging、序貫高斯等，求解地質(zhì)體內(nèi)部物理、化學(xué)屬性參數(shù)模型。利用地質(zhì)體物性建模成果，可以有效的發(fā)現(xiàn)儲(chǔ)層，并對(duì)油氣資源分布情況進(jìn)行準(zhǔn)確預(yù)測，從而指導(dǎo)油田開發(fā)。考慮到地層成因，地質(zhì)體中某一區(qū)域的物性，與同一層物性的關(guān)系遠(yuǎn)遠(yuǎn)大于該區(qū)域的上層或下層。為了更精確的表達(dá)地質(zhì)體物性，實(shí)際建模過程中，在對(duì)地質(zhì)體細(xì)分小層的基礎(chǔ)上，會(huì)進(jìn)一步細(xì)分，對(duì)細(xì)分后的每一層，利用由測井曲線數(shù)據(jù)得到的關(guān)于該層的物性信息，使用空間差值算法進(jìn)行物性建模。將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型是利用測井?dāng)?shù)據(jù)，為井穿過的儲(chǔ)層網(wǎng)格單元賦值。被賦值的網(wǎng)格單元，在進(jìn)行物性建模時(shí)，作為已知量，進(jìn)而確定其他未被測量過區(qū)域的屬性值?？梢?，它是進(jìn)行地質(zhì)體物性建模的基礎(chǔ)。
[0003]目前，角點(diǎn)網(wǎng)格是被廣泛使用的描述儲(chǔ)集層的網(wǎng)格模型。作為結(jié)構(gòu)化網(wǎng)格模型，它具有易于實(shí)現(xiàn)的優(yōu)點(diǎn)。高度的靈活性使得它能夠很好的表征儲(chǔ)集層地質(zhì)特征，比如非均質(zhì)性以及斷層。
[0004]然而，由于角點(diǎn)網(wǎng)格不使用笛卡爾坐標(biāo)系，在離散測井?dāng)?shù)據(jù)時(shí)非常復(fù)雜繁瑣。再加上為了滿足對(duì)建模區(qū)域精細(xì)描述的需求，現(xiàn)在的網(wǎng)格劃分都非常細(xì)密。故而，一種快速地將測井?dāng)?shù)據(jù)離散到儲(chǔ)層模型的方法被迫切的需要。
【發(fā)明內(nèi)容】
[0005](一)解決的技術(shù)問題
[0006]針對(duì)現(xiàn)有技術(shù)的不足，本發(fā)明提供了一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置，根據(jù)測井?dāng)?shù)據(jù)值，用盡可能少的時(shí)間為井軌跡所穿過的網(wǎng)格賦值，從而減少建模人員的等待時(shí)間，提聞工作效率。
[0007](二)技術(shù)方案
[0008]為實(shí)現(xiàn)以上目的，本發(fā)明通過以下技術(shù)方案予以實(shí)現(xiàn):
[0009]一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法，所述儲(chǔ)層網(wǎng)格模型為角點(diǎn)網(wǎng)格，包括下列步驟:
[0010]S1、獲取井軌跡穿越的網(wǎng)格單元；
[0011]S2、計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)，并利用所述交點(diǎn)將所述井軌跡劃分為多個(gè)區(qū)塊；
[0012]S3、對(duì)落入每一區(qū)塊的所述井軌跡上的測井曲線屬性值進(jìn)行加權(quán)平均，求得屬性均值，并將所述屬性均值賦給對(duì)應(yīng)的網(wǎng)格單元。[0013]優(yōu)選的，步驟SI包括:
[0014]S11、根據(jù)儲(chǔ)層網(wǎng)格模型的網(wǎng)格厚度信息設(shè)置井軌跡中待測點(diǎn)選取的間隔；
[0015]S12、按照所述間隔順次取出待測點(diǎn)；
[0016]S13、找到每一個(gè)待測點(diǎn)所在的網(wǎng)格單元；
[0017]S14、將所述網(wǎng)格單元進(jìn)行整理與合并。
[0018]優(yōu)選的，步驟Sll包括:
[0019]S111、根據(jù)井頭坐標(biāo)，以及儲(chǔ)層網(wǎng)格模型的邊界信息和網(wǎng)格長寬設(shè)置信息，估算出位于井頭下方的網(wǎng)格單元；
[0020]SI 12、計(jì)算井頭下方網(wǎng)格單元的平均厚度；
[0021 ] S113、結(jié)合井軌跡數(shù)據(jù)的采樣間隔，計(jì)算出待測點(diǎn)選取的間隔；
[0022]其中，步驟S13包括:利用參考索引找出對(duì)應(yīng)網(wǎng)格單元，以網(wǎng)格單元為中心，采用臨近搜索法選取網(wǎng)格單元，對(duì)待測點(diǎn)和選取的網(wǎng)格單元進(jìn)行包含性測試；
[0023]其中，步驟S14包括:通過網(wǎng)格單元的索引信息，利用折半查找法，找出未被測出的網(wǎng)格單元；對(duì)于多個(gè)待測點(diǎn)落入同一網(wǎng)格單元的情形，只記錄第一個(gè)落入該網(wǎng)格單元的待測點(diǎn)信息，與最后一個(gè)落入該網(wǎng)格單元的待測點(diǎn)信息。
[0024]優(yōu)選的，步驟S2中計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)包括步驟:
[0025]根據(jù)網(wǎng)格單元索引信息確定單元內(nèi)同井軌跡有交點(diǎn)的面；選取面兩端的待測點(diǎn)，計(jì)算兩待測點(diǎn)間的井軌跡與網(wǎng)格單元面的交點(diǎn)。
[0026]優(yōu)選的，計(jì)算兩待測點(diǎn)間的井軌跡與網(wǎng)格單元面的交點(diǎn)時(shí)使用快速逼近法，實(shí)現(xiàn)方式是:
[0027]選取兩相鄰井軌跡數(shù)據(jù)點(diǎn)組成直線，計(jì)算同網(wǎng)格單元表面的交點(diǎn)，計(jì)算所述交點(diǎn)與數(shù)據(jù)點(diǎn)的距離；
[0028]在井軌跡上選取距離相等的數(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)算；
[0029]當(dāng)交點(diǎn)位于選取的兩相鄰井軌跡數(shù)據(jù)點(diǎn)之間時(shí)，結(jié)束運(yùn)算。
[0030]優(yōu)選的，步驟S2中將井軌跡劃分為多個(gè)區(qū)塊時(shí)，將兩個(gè)相鄰的井與網(wǎng)格單元面的交點(diǎn)之間的區(qū)域，作為一個(gè)區(qū)塊，每一區(qū)塊對(duì)應(yīng)儲(chǔ)層網(wǎng)格模型的一個(gè)網(wǎng)格單元。
[0031]優(yōu)選的，步驟S3中計(jì)算落入每一區(qū)塊的測井曲線屬性值均值時(shí)，得到唯一的數(shù)值，并以此數(shù)值作為區(qū)塊對(duì)應(yīng)網(wǎng)格單元的數(shù)值。
[0032]本發(fā)明還提供了一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的裝置，包括:
[0033]網(wǎng)格單元查找模塊，用于找出井軌跡所穿越的網(wǎng)格單元，
[0034]測井?dāng)?shù)據(jù)劃分模塊，用于計(jì)算井軌跡同網(wǎng)格單元面的交點(diǎn)，并對(duì)測井?dāng)?shù)據(jù)劃分區(qū)塊；
[0035]屬性計(jì)算模塊，用于對(duì)落入每一區(qū)塊的所有測井曲線屬性值計(jì)算加權(quán)平均值，并為對(duì)應(yīng)的網(wǎng)格單元賦值；
[0036]交互模塊，對(duì)儲(chǔ)層網(wǎng)格模型和測井?dāng)?shù)據(jù)進(jìn)行訪問。
[0037]優(yōu)選的，所述網(wǎng)格單元查找模塊包括:
[0038]待測點(diǎn)間隔發(fā)生器，用于產(chǎn)生待測點(diǎn)選取間隔；
[0039]待測點(diǎn)抽樣采集器，用于根據(jù)指定間隔選取待測點(diǎn)的抽樣采集；[0040]網(wǎng)格單元搜索裝置，用于準(zhǔn)確、快速判斷待測點(diǎn)所處的網(wǎng)格單元的測試設(shè)備和網(wǎng)格單元搜索；
[0041]網(wǎng)格單元內(nèi)判定器，用于判斷待測點(diǎn)是否位于網(wǎng)格單元內(nèi)；
[0042]其中，所述測井?dāng)?shù)據(jù)劃分模塊包括:
[0043]交點(diǎn)求解器，用于快速計(jì)算井軌跡與網(wǎng)格單元面的交點(diǎn)；
[0044]測井?dāng)?shù)據(jù)分塊器，用于對(duì)測井?dāng)?shù)據(jù)劃分區(qū)塊；
[0045]其中，所述屬性計(jì)算模塊包括:
[0046]屬性求解器，用于計(jì)算落入同一區(qū)塊的測井曲線屬性加權(quán)平均值；
[0047]屬性賦值器，用于將區(qū)塊屬性值設(shè)為對(duì)應(yīng)網(wǎng)格單元屬性值；
[0048]其中，所述交互模塊包括:
[0049]儲(chǔ)層模型訪問子模塊，用于根據(jù)儲(chǔ)層網(wǎng)格模型的三維索引值，獲取儲(chǔ)層網(wǎng)格模型中的特定網(wǎng)格單元的信息，且將指定物性屬性值設(shè)置到對(duì)應(yīng)網(wǎng)格單元中；
[0050]測井?dāng)?shù)據(jù)訪問子模塊，用于提供井軌跡類型數(shù)據(jù)和測井曲線類型數(shù)據(jù)的訪問接□。
[0051]優(yōu)選的，所述測井?dāng)?shù)據(jù)訪問子模塊按照采樣點(diǎn)在對(duì)應(yīng)井軌跡類型數(shù)據(jù)和測井曲線類型數(shù)據(jù)中的索引，獲得采樣點(diǎn)的特定信息；其中所述井軌跡數(shù)據(jù)包含:測量深度信息、大地坐標(biāo)系下的三維坐標(biāo)值和垂直深度；所述測井曲線數(shù)據(jù)包括:測量深度信息、孔隙度值、自然伽馬值和沉積相相位值。
[0052](三)有益效果
[0053]本發(fā)明提供了一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法和裝置，通過獲取井軌跡穿越的網(wǎng)格單元，計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)，并利用所述交點(diǎn)將所述井軌跡劃分為多個(gè)區(qū)塊；進(jìn)而對(duì)落入每一區(qū)塊的所述井軌跡上的測井曲線屬性值進(jìn)行加權(quán)平均，求得屬性均值，并將所述屬性均值賦給對(duì)應(yīng)的網(wǎng)格單元。
[0054]本發(fā)明可以快捷的利用測井?dāng)?shù)據(jù)計(jì)算儲(chǔ)層網(wǎng)格模型中被井穿過的網(wǎng)格單元的屬性值，減少了建模人員的等待時(shí)間，提供工作效率。
[0055]本發(fā)明可以準(zhǔn)確的利用測井?dāng)?shù)據(jù)計(jì)算儲(chǔ)層網(wǎng)格模型中相應(yīng)單元的屬性值，克服了角點(diǎn)網(wǎng)格體系下難以計(jì)算空間位置的缺點(diǎn)，提高了建模結(jié)果的準(zhǔn)確性。大大增加了石油儲(chǔ)層精細(xì)描述系統(tǒng)的可信度，為油井采收、油田二次開發(fā)等工作提供可靠的支撐，從而大幅節(jié)約油田開采成本。
【專利附圖】
【附圖說明】
[0056]為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他的附圖。
[0057]圖1為本發(fā)明實(shí)施例中將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法的流程示意圖；
[0058]圖2為雙線性曲面的示意圖；
[0059]圖3為本發(fā)明實(shí)施例中將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的裝置的結(jié)構(gòu)示意圖；
[0060]圖4為網(wǎng)格單元查找模塊的實(shí)施流程圖；[0061]圖5為計(jì)算角點(diǎn)網(wǎng)格單元平均厚度的示意圖；
[0062]圖6為計(jì)算射線與雙線性曲面交點(diǎn)方法的示意圖；
[0063]圖7為判斷待測點(diǎn)是否位于指定網(wǎng)格單元方法的示意圖；
[0064]圖8為測井?dāng)?shù)據(jù)劃分模塊的實(shí)施流程圖；
[0065]圖9為用快速逼近法求解井軌跡與網(wǎng)格單元面交點(diǎn)過程的示意圖；
[0066]圖10為屬性計(jì)算模塊的實(shí)施流程圖；
[0067]圖11為當(dāng)沒有測井曲線數(shù)據(jù)點(diǎn)落入某一區(qū)塊時(shí)確定網(wǎng)格單元屬性值的方法示意圖。
【具體實(shí)施方式】
[0068]為使本發(fā)明實(shí)施例的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚，下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0069]實(shí)施例1:
[0070]如圖1所示，本發(fā)明實(shí)施例提供了一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的方法，所述儲(chǔ)層網(wǎng)格模型為角點(diǎn)網(wǎng)格，包括下列步驟:
[0071]S1、獲取井軌跡穿越的網(wǎng)格單元；
[0072]S2、計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)，并利用所述交點(diǎn)將所述井軌跡劃分為多個(gè)區(qū)塊；
[0073]S3、對(duì)落入每一區(qū)塊的所述井軌跡上的測井曲線屬性值進(jìn)行加權(quán)平均，求得屬性均值，并將所述屬性均值賦給對(duì)應(yīng)的網(wǎng)格單元。
[0074]本發(fā)明實(shí)施例可以快捷的利用測井?dāng)?shù)據(jù)計(jì)算儲(chǔ)層網(wǎng)格模型中被井穿過的網(wǎng)格單元的屬性值，減少了建模人員的等待時(shí)間，提供工作效率?？梢詼?zhǔn)確的利用測井?dāng)?shù)據(jù)計(jì)算儲(chǔ)層網(wǎng)格模型中相應(yīng)單元的屬性值，克服了角點(diǎn)網(wǎng)格體系下難以計(jì)算空間位置的缺點(diǎn)，提高了建模結(jié)果的準(zhǔn)確性。大大增加了石油儲(chǔ)層精細(xì)描述系統(tǒng)的可信度，為油井采收、油田二次開發(fā)等工作提供可靠的支撐，從而大幅節(jié)約油田開采成本。
[0075]下面對(duì)本發(fā)明實(shí)施例進(jìn)行詳細(xì)的說明:
[0076]優(yōu)選的，步驟SI包括:
[0077]S11、根據(jù)儲(chǔ)層網(wǎng)格模型的網(wǎng)格厚度信息設(shè)置井軌跡中待測點(diǎn)選取的間隔；
[0078]S12、按照所述間隔順次取出待測點(diǎn)；
[0079]S13、找到每一個(gè)待測點(diǎn)所在的網(wǎng)格單元；
[0080]S14、將所述網(wǎng)格單元進(jìn)行整理與合并。
[0081]優(yōu)選的，步驟Sll包括:
[0082]S111、根據(jù)井頭坐標(biāo)，以及儲(chǔ)層網(wǎng)格模型的邊界信息和網(wǎng)格長寬設(shè)置信息，估算出位于井頭下方的網(wǎng)格單元；
[0083]SI 12、計(jì)算井頭下方網(wǎng)格單元的平均厚度；
[0084]S113、結(jié)合井軌跡數(shù)據(jù)的采樣間隔，計(jì)算出待測點(diǎn)選取的間隔。
[0085]優(yōu)選的，步驟S13包括:利用參考索引找出對(duì)應(yīng)網(wǎng)格單元，以網(wǎng)格單元為中心，采用臨近搜索法選取網(wǎng)格單元，對(duì)待測點(diǎn)和選取的網(wǎng)格單元進(jìn)行包含性測試。
[0086]優(yōu)選的，步驟S14包括:
[0087]通過網(wǎng)格單元的索引信息，利用折半查找法，找出未被測出的網(wǎng)格單元；對(duì)于多個(gè)待測點(diǎn)落入同一網(wǎng)格單元的情形，只記錄第一個(gè)落入該網(wǎng)格單元的待測點(diǎn)信息，與最后一個(gè)落入該網(wǎng)格單元的待測點(diǎn)信息。
[0088]優(yōu)選的，步驟S2中計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)包括步驟:
[0089]根據(jù)網(wǎng)格單元索引信息確定單元內(nèi)同井軌跡有交點(diǎn)的面；選取面兩端的待測點(diǎn)，計(jì)算兩待測點(diǎn)間的井軌跡與網(wǎng)格單元面的交點(diǎn)。
[0090]優(yōu)選的，計(jì)算兩待測點(diǎn)間的井軌跡與網(wǎng)格單元面的交點(diǎn)時(shí)使用快速逼近法，實(shí)現(xiàn)方式是:
[0091]選取兩相鄰井軌跡數(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)算。直到交點(diǎn)位于選取的兩相鄰井軌跡數(shù)據(jù)點(diǎn)之間，算法結(jié)束。
[0092]優(yōu)選的，步驟S2中將井軌跡劃分為多個(gè)區(qū)塊時(shí)，將兩個(gè)相鄰的井與網(wǎng)格單元面的交點(diǎn)之間的區(qū)域，作為一個(gè)區(qū)塊，每一區(qū)塊對(duì)應(yīng)儲(chǔ)層網(wǎng)格模型的一個(gè)網(wǎng)格單元。
[0093]優(yōu)選的，步驟S3中計(jì)算落入每一區(qū)塊的測井曲線屬性值均值時(shí)，得到唯一的數(shù)值，并以此數(shù)值作為區(qū)塊對(duì)應(yīng)網(wǎng)格單元的數(shù)值。
[0094]實(shí)施例2:
[0095]如圖3所示，本發(fā)明還提供了一種將測井?dāng)?shù)據(jù)離散到儲(chǔ)層網(wǎng)格模型的裝置，包括:
[0096]網(wǎng)格單元查找模塊，用于找出井軌跡所穿越的網(wǎng)格單元，
[0097]測井?dāng)?shù)據(jù)劃分模塊，用于計(jì)算井軌跡同網(wǎng)格單元面的交點(diǎn)，并對(duì)測井?dāng)?shù)據(jù)劃分區(qū)塊；
[0098]屬性計(jì)算模塊，用于對(duì)落入每一區(qū)塊的所有測井曲線屬性值計(jì)算加權(quán)平均值，并為對(duì)應(yīng)的網(wǎng)格單元賦值；
[0099]交互模塊，對(duì)儲(chǔ)層網(wǎng)格模型和測井?dāng)?shù)據(jù)進(jìn)行訪問。
[0100]優(yōu)選的，所述網(wǎng)格單元查找模塊包括:
[0101]待測點(diǎn)間隔發(fā)生器，用于產(chǎn)生待測點(diǎn)選取間隔；
[0102]待測點(diǎn)抽樣采集器，用于根據(jù)指定間隔選取待測點(diǎn)的抽樣采集；
[0103]網(wǎng)格單元搜索裝置，用于準(zhǔn)確、快速判斷待測點(diǎn)所處的網(wǎng)格單元的測試設(shè)備和網(wǎng)格單元搜索；
[0104]網(wǎng)格單元內(nèi)判定器，用于判斷待測點(diǎn)是否位于網(wǎng)格單元內(nèi)。
[0105]優(yōu)選的，所述測井?dāng)?shù)據(jù)劃分模塊包括:
[0106]交點(diǎn)求解器，用于快速計(jì)算井軌跡與網(wǎng)格單元面的交點(diǎn)；
[0107]測井?dāng)?shù)據(jù)分塊器，用于對(duì)測井?dāng)?shù)據(jù)劃分區(qū)塊。
[0108]優(yōu)選的，所述屬性計(jì)算模塊包括:
[0109]屬性求解器，用于計(jì)算落入同一區(qū)塊的測井曲線屬性加權(quán)平均值；
[0110]屬性賦值器，用于將區(qū)塊屬性值設(shè)為對(duì)應(yīng)網(wǎng)格單元屬性值。
[0111]優(yōu)選的，所述交互模塊包括:[0112]儲(chǔ)層模型訪問子模塊，用于根據(jù)儲(chǔ)層網(wǎng)格模型的三維索引值，獲取儲(chǔ)層網(wǎng)格模型中的特定網(wǎng)格單元的信息，且將指定物性屬性值設(shè)置到對(duì)應(yīng)網(wǎng)格單元中；
[0113]測井?dāng)?shù)據(jù)訪問子模塊，用于提供井軌跡類型數(shù)據(jù)和測井曲線類型數(shù)據(jù)的訪問接□。
[0114]優(yōu)選的，所述測井?dāng)?shù)據(jù)訪問子模塊按照采樣點(diǎn)在對(duì)應(yīng)井軌跡類型數(shù)據(jù)和測井曲線類型數(shù)據(jù)中的索引，獲得采樣點(diǎn)的特定信息；其中所述井軌跡數(shù)據(jù)包含:測量深度信息、大地坐標(biāo)系下的三維坐標(biāo)值和垂直深度；所述測井曲線數(shù)據(jù)包括:測量深度信息、孔隙度值、自然伽馬值和沉積相相位值。
[0115]下面對(duì)各模塊的具體工作過程進(jìn)行描述:
[0116]網(wǎng)格單元查找模塊的實(shí)施流程圖，如圖4所示，包含以下步驟:
[0117]1)獲取待離散測井的井口坐標(biāo)信息，利用儲(chǔ)層模型的邊界信息和建模時(shí)指定的網(wǎng)格在X方向、y方向的增量，獲取參考索引信息。參考索引信息的含義是根據(jù)該索引獲得的儲(chǔ)層模型網(wǎng)格單元位于或近似的位于井口的正下方；
[0118]2)利用參考索引獲取儲(chǔ)層中對(duì)應(yīng)的網(wǎng)格單元，計(jì)算網(wǎng)格單元在z方向上的平均厚度zSpan ;以圖5為例，計(jì)算公式選擇:
[0119]
【權(quán)利要求】
1.一種將測井?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ū)塊； 53、對(duì)落入每一區(qū)塊的所述井軌跡上的測井曲線屬性值進(jìn)行加權(quán)平均，求得屬性均值，并將所述屬性均值賦給對(duì)應(yīng)的網(wǎng)格單元。
2.如權(quán)利要求1所述的方法，其特征在于，步驟SI包括: S11、根據(jù)儲(chǔ)層網(wǎng)格模型的網(wǎng)格厚度信息設(shè)置井軌跡中待測點(diǎn)選取的間隔； S12、按照所述間隔順次取出待測點(diǎn)； S13、找到每一個(gè)待測點(diǎn)所在的網(wǎng)格單元； S14、將所述網(wǎng)格單元進(jìn)行整理與合并。
3.如權(quán)利要求2所述的方法，其特征在于，其中，步驟Sll包括: sill、根據(jù)井頭坐標(biāo)，以及儲(chǔ)層網(wǎng)格模型的邊界信息和網(wǎng)格長寬設(shè)置信息，估算出位于井頭下方的網(wǎng)格單元； S112、計(jì)算井頭下方網(wǎng)格單元的平均厚度； S113、結(jié)合井軌跡數(shù)據(jù)的采樣間隔，計(jì)算出待測點(diǎn)選取的間隔； 其中，步驟S13包括:利用參考索引找出對(duì)應(yīng)網(wǎng)格單元，以網(wǎng)格單元為中心，采用臨近搜索法選取網(wǎng)格單元，對(duì)待測點(diǎn)和選取的網(wǎng)格單元進(jìn)行包含性測試； 其中，步驟S14包括:通過網(wǎng)格單元的索引信息，利用折半查找法，找出未被測出的網(wǎng)格單元；對(duì)于多個(gè)待測點(diǎn)落入同一網(wǎng)格單元的情形，只記錄第一個(gè)落入該網(wǎng)格單元的待測點(diǎn)信息，與最后一個(gè)落入該網(wǎng)格單元的待測點(diǎn)信息。
4.如權(quán)利要求1所述的方法，其特征在于，其特征在于，步驟S2中計(jì)算所述井軌跡同所述網(wǎng)格單元的交點(diǎn)包括步驟: 根據(jù)網(wǎng)格單元索引信息確定單元內(nèi)同井軌跡有交點(diǎn)的面；選取面兩端的待測點(diǎn)，計(jì)算兩待測點(diǎn)間的井軌跡與網(wǎng)格單元面的交點(diǎn)。
5.如權(quán)利要求4所述的方法，其特征在于，計(jì)算兩待測點(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)算。
6.如權(quán)利要求1所述的方法，其特征在于，步驟S2中將井軌跡劃分為多個(gè)區(qū)塊時(shí)，將兩個(gè)相鄰的井與網(wǎng)格單元面的交點(diǎn)之間的區(qū)域，作為一個(gè)區(qū)塊，每一區(qū)塊對(duì)應(yīng)儲(chǔ)層網(wǎng)格模型的一個(gè)網(wǎng)格單元。
7.如權(quán)利要求1所述的方法，其特征在于，步驟S3中計(jì)算落入每一區(qū)塊的測井曲線屬性值均值時(shí)，得到唯一的數(shù)值，并以此數(shù)值作為區(qū)塊對(duì)應(yīng)網(wǎng)格單元的數(shù)值。
8.一種將測井?dāng)?shù)據(jù)離散 到儲(chǔ)層網(wǎng)格模型的裝置，其特征在于，包括:網(wǎng)格單元查找模塊，用于找出井軌跡所穿越的網(wǎng)格單元， 測井?dāng)?shù)據(jù)劃分模塊，用于計(jì)算井軌跡同網(wǎng)格單元面的交點(diǎn)，并對(duì)測井?dāng)?shù)據(jù)劃分區(qū)塊；屬性計(jì)算模塊，用于對(duì)落入每一區(qū)塊的所有測井曲線屬性值計(jì)算加權(quán)平均值，并為對(duì)應(yīng)的網(wǎng)格單元賦值； 交互模塊，對(duì)儲(chǔ)層網(wǎng)格模型和測井?dāng)?shù)據(jù)進(jìn)行訪問。
9.如權(quán)利要求8所述的裝置，其特征在于，其中，所述網(wǎng)格單元查找模塊包括: 待測點(diǎn)間隔發(fā)生器，用于產(chǎn)生待測點(diǎn)選取間隔； 待測點(diǎn)抽樣采集器，用于根據(jù)指定間隔選取待測點(diǎn)的抽樣采集； 網(wǎng)格單元搜索裝置，用于準(zhǔn)確、快速判斷待測點(diǎn)所處的網(wǎng)格單元的測試設(shè)備和網(wǎng)格單元搜索； 網(wǎng)格單元內(nèi)判定器，用于判斷待測點(diǎn)是否位于網(wǎng)格單元內(nèi)； 其中，所述測井?dāng)?shù)據(jù)劃分模塊包括: 交點(diǎn)求解器，用于快速計(jì)算井軌跡與網(wǎng)格單元面的交點(diǎn)； 測井?dāng)?shù)據(jù)分塊器，用于對(duì)測井?dāng)?shù)據(jù)劃分區(qū)塊； 其中，所述屬性計(jì)算模塊包括: 屬性求解器，用于計(jì)算落入同一區(qū)塊的測井曲線屬性加權(quán)平均值； 屬性賦值器，用于將區(qū)塊屬性值設(shè)為對(duì)應(yīng)網(wǎng)格單元屬性值； 其中，所述交互模塊包括: 儲(chǔ)層模型訪問子模塊，用于根據(jù)儲(chǔ)層網(wǎng)格模型的三維索引值，獲取儲(chǔ)層網(wǎng)格模型中的特定網(wǎng)格單元的信息，且將指定物性屬性值設(shè)置到對(duì)應(yīng)網(wǎng)格單元中； 測井?dāng)?shù)據(jù)訪問子模塊，用于提供井軌跡類型數(shù)據(jù)和測井曲線類型數(shù)據(jù)的訪問接口。
10.如權(quán)利要求9所述的裝置，其特征在于，所述測井?dāng)?shù)據(jù)訪問子模塊按照采樣點(diǎn)在對(duì)應(yīng)井軌跡類型數(shù)據(jù)和測井曲線類型數(shù)據(jù)中的索引，獲得采樣點(diǎn)的特定信息；其中所述井軌跡數(shù)據(jù)包含:測量深度信息、大地坐標(biāo)系下的三維坐標(biāo)值和垂直深度；所述測井曲線數(shù)據(jù)包括:測量深度信息、孔隙度值、自然伽馬值和沉積相相位值。
【文檔編號(hào)】G06F17/50GK103886129SQ201410057529
【公開日】2014年6月25日 申請(qǐng)日期:2014年2月19日 優(yōu)先權(quán)日:2014年2月19日
【發(fā)明者】倪廣元, 王劍秦, 董肇偉 申請(qǐng)人:北京中恒永信科技有限公司