一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模方法及系統(tǒng)的制作方法
【專利摘要】本發(fā)明提供一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法及系統(tǒng)，所述方法包括：獲取當(dāng)前砂體的測(cè)井?dāng)?shù)據(jù)；對(duì)測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，得到測(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)；對(duì)廢棄河道數(shù)據(jù)進(jìn)行擬合，得到平面分布的深泓線；根據(jù)深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體；根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體；根據(jù)廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型；對(duì)矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。本發(fā)明的模擬結(jié)果可以刻畫不同規(guī)模層次的非均質(zhì)體，可以快速的對(duì)模型進(jìn)行調(diào)節(jié)。
【專利說明】一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模方法及系統(tǒng)
【技術(shù)領(lǐng)域】
[0001]本發(fā)明關(guān)于油田開采【技術(shù)領(lǐng)域】，特別是關(guān)于儲(chǔ)層的空間建模技術(shù)，具體的講是一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模方法及系統(tǒng)。
【背景技術(shù)】
[0002]三維儲(chǔ)層建模，即建立儲(chǔ)層特征三維分布的數(shù)字化模型，其本質(zhì)是基于三維網(wǎng)格表征儲(chǔ)層特征的分布，其成果是三維數(shù)據(jù)體。這一技術(shù)是上世紀(jì)80年代隨著計(jì)算機(jī)技術(shù)的發(fā)展而發(fā)展起來的?；谟?jì)算機(jī)存儲(chǔ)和顯示技術(shù)，將儲(chǔ)層三維網(wǎng)塊化后，對(duì)各個(gè)網(wǎng)格(grid)賦以各自的儲(chǔ)層參數(shù)值，并按三維空間分布位置存入計(jì)算機(jī)內(nèi)，形成了三維數(shù)據(jù)體，這樣就可以進(jìn)行儲(chǔ)層的三維顯示，可以任意切片和切剖面(不同層位、不同方向剖面)，以及進(jìn)行各種運(yùn)算和分析。對(duì)于一個(gè)給定的儲(chǔ)層，按X，Y，Z三個(gè)方向劃分網(wǎng)格，將儲(chǔ)層在三維空間劃分成眾多網(wǎng)格，即為三維網(wǎng)格化(如圖13所示)。
[0003]在儲(chǔ)層地質(zhì)學(xué)的研究中，對(duì)于儲(chǔ)層在三維空間的表征主要是應(yīng)用地質(zhì)統(tǒng)計(jì)學(xué)的理論和方法，綜合多學(xué)科的信息對(duì)儲(chǔ)層非均質(zhì)的空間分布進(jìn)行確定性和隨機(jī)性的賦值。確定性的途徑主要應(yīng)用各種克里金方法，隨機(jī)性的途徑則主要是應(yīng)用各種高斯模擬(如截?cái)喔咚鼓M、序貫指示模擬)方法。按照隨機(jī)模擬基本單元不同，隨機(jī)建模方法可以分為基于目標(biāo)的方法和基于象元的方法兩類。
[0004]基于象元的方法的模擬基本單元為象元(cell)，即網(wǎng)格化儲(chǔ)層模型中單個(gè)網(wǎng)格(grid)，借助地質(zhì)統(tǒng)計(jì)學(xué)中的變差函數(shù)，以不同變差函數(shù)表征變量的空間變異特征，逐網(wǎng)格分配儲(chǔ)層屬性值。基于象元的模擬方法中，對(duì)于河流相應(yīng)用較多的方法主要為序貫指示模擬。盡管序貫指示模擬應(yīng)用相當(dāng)廣泛，然而由于它是以變差函數(shù)為依據(jù)來反映儲(chǔ)層形態(tài)及空間分布特征，而變差函數(shù)只考慮了兩點(diǎn)相關(guān)性。地質(zhì)現(xiàn)象復(fù)雜性決定了對(duì)其描述應(yīng)該是基于更多點(diǎn)信息的，而變差函數(shù)很難真實(shí)反映儲(chǔ)層地下結(jié)構(gòu)及形態(tài)，難于再現(xiàn)復(fù)雜地質(zhì)的界面信息。
[0005]基于目標(biāo)的方法的模擬基本單元為目標(biāo)物體(即具有離散性質(zhì)的地質(zhì)特征，如沉積相、流動(dòng)單元等)，通過對(duì)目標(biāo)幾何形態(tài)(如長、寬、厚及其之間定量關(guān)系)的定量研究，在建模過程中通過對(duì)多個(gè)網(wǎng)格同時(shí)賦予生目標(biāo)物體而建立儲(chǔ)層的三維模型。但是，依沉積學(xué)原理，沉積砂體的橫向?qū)挾?、厚度受其沉積時(shí)的水動(dòng)力條件制約，兩者之間存在一定的定量統(tǒng)計(jì)關(guān)系，基于目標(biāo)的算法對(duì)沉積因素考慮不夠。
[0006]另外的學(xué)者則對(duì)基于變差函數(shù)的兩點(diǎn)統(tǒng)計(jì)學(xué)儲(chǔ)層建模的改進(jìn)，提出了多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)。應(yīng)用“數(shù)據(jù)事件”代替變差函數(shù)表達(dá)地質(zhì)變量的空間結(jié)構(gòu)性，克服傳統(tǒng)地質(zhì)統(tǒng)計(jì)學(xué)不能再現(xiàn)目標(biāo)幾何形態(tài)的不足。但是多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)作為一種新的建模思路和方法，在再現(xiàn)河道的連續(xù)性方面仍需改進(jìn)，并且對(duì)各種非平穩(wěn)現(xiàn)象(具有趨勢(shì)分布的)的模擬方面仍然存在困難，尤其在對(duì)地質(zhì)體的三維形態(tài)表達(dá)和模擬方面仍需改進(jìn)。
[0007]儲(chǔ)層構(gòu)型建模則多采用序貫指示模擬與人機(jī)結(jié)合再處理的手段，特點(diǎn)是工作量大，且無法表征模型中不確定的因素。也有學(xué)者提出了二次網(wǎng)格加密方法、基于沉積過程的模、層次約束下的隨機(jī)抽樣等隨機(jī)建模方法。
[0008]但各種方法仍然沿用了傳統(tǒng)的基于象元的建模方法。即首先定義網(wǎng)格，之后將構(gòu)型要素作為一種離散屬性定義在網(wǎng)格之上。因此在網(wǎng)格化研究區(qū)時(shí)必須考慮到所定義的網(wǎng)格能予以刻畫最關(guān)鍵的非均質(zhì)體。對(duì)曲流河儲(chǔ)集層構(gòu)型模型而言，雖然點(diǎn)壩內(nèi)部的側(cè)積層厚度較小(0.1m?1.0m級(jí))，但卻具有一定規(guī)模的橫向延伸范圍，是造成點(diǎn)壩砂體剩余油富集的主要因素。在精細(xì)儲(chǔ)層構(gòu)型建模中，為了刻畫此種規(guī)模的非均質(zhì)體，選擇之一便是細(xì)化網(wǎng)格，這樣勢(shì)必會(huì)造成網(wǎng)格節(jié)點(diǎn)個(gè)數(shù)的膨脹。雖然計(jì)算機(jī)有了很大的發(fā)展，對(duì)于一個(gè)縱向、橫向上劃分的網(wǎng)格精度還是有一定的限制，并不能采取任意細(xì)化的網(wǎng)格。例如，對(duì)于一個(gè)IkmX IkmX 50m的模型，如果使用5mX 5mX Im的建模網(wǎng)格，網(wǎng)格總數(shù)將達(dá)到200 X IO4個(gè)。即使這樣的網(wǎng)格依然無法刻畫垂向上小于0.5m的非均值體。因此，在建模階段采用基于網(wǎng)格的方法導(dǎo)致網(wǎng)格數(shù)目急劇增加，增加了計(jì)算時(shí)間，也限制了隨機(jī)模擬實(shí)現(xiàn)的個(gè)數(shù)。雖然已有隨機(jī)建模算法對(duì)高級(jí)別構(gòu)型要素，如單一河道砂體的三維建模有較好的效果，但難于滿足砂體內(nèi)部構(gòu)型建模解剖與建模的需要。尤其對(duì)于微小級(jí)次的構(gòu)型單元，人工繪制進(jìn)行三維建模的工作量巨大，效率低下。

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

[0009]為了克服現(xiàn)有技術(shù)存在的上述問題，本發(fā)明提供了一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法及系統(tǒng)，對(duì)模擬地質(zhì)體的定義采用了基于實(shí)數(shù)集的點(diǎn)、線、面所限定的體，因而模型在模擬過程中不受網(wǎng)格尺寸的限制，使得模擬結(jié)果可以刻畫不同規(guī)模層次的非均質(zhì)體，可以快速的對(duì)模型進(jìn)行調(diào)節(jié)，在模擬過程中只是對(duì)地質(zhì)體的空間位置和形態(tài)的定義參數(shù)進(jìn)行調(diào)節(jié)。
[0010]本發(fā)明的目的之一是，提供一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法，包括:獲取當(dāng)前砂體的測(cè)井?dāng)?shù)據(jù)；對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)；對(duì)所述的廢棄河道數(shù)據(jù)進(jìn)行擬合，得到平面分布的深泓線；根據(jù)所述的深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體；根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體；根據(jù)所述廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型；對(duì)所述的矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。
[0011]本發(fā)明的目的之一是，提供了一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)，包括:測(cè)井?dāng)?shù)據(jù)獲取裝置，用于獲取當(dāng)前砂體的測(cè)井?dāng)?shù)據(jù)；測(cè)井?dāng)?shù)據(jù)處理裝置，用于對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)；擬合裝置，用于對(duì)所述的廢棄河道數(shù)據(jù)進(jìn)行擬合，得到平面分布的深泓線；三維幾何體建立裝置，用于根據(jù)所述的深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體；空間分布體建立裝置，用于根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體；矢量化模型確定裝置，用于根據(jù)所述廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型；網(wǎng)格化裝置，用于對(duì)所述的矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。
[0012]本發(fā)明的有益效果在于，提供了一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法及系統(tǒng)，與現(xiàn)有技術(shù)相比，該方法提出的主要目的是在儲(chǔ)層構(gòu)型建模中描述不同級(jí)別的非均質(zhì)體以及準(zhǔn)確刻畫各級(jí)地質(zhì)體的復(fù)雜空間形態(tài)。在模擬階段的地質(zhì)體的存儲(chǔ)借鑒了矢量圖形存儲(chǔ)的思路，即對(duì)模擬地質(zhì)體的定義采用了基于實(shí)數(shù)集的點(diǎn)、線、面所限定的體，因而模型在模擬過程中不受網(wǎng)格尺寸的限制，使得模擬結(jié)果可以刻畫不同規(guī)模層次的非均質(zhì)體。應(yīng)用基于矢量體存儲(chǔ)方法的另外一個(gè)優(yōu)點(diǎn)是可以快速的對(duì)模型進(jìn)行調(diào)節(jié)，因?yàn)樵谀M過程中只是對(duì)地質(zhì)體的空間位置和形態(tài)的定義參數(shù)進(jìn)行調(diào)節(jié)，因而這樣的調(diào)節(jié)方式較基于網(wǎng)格的調(diào)節(jié)更便捷，從而使結(jié)果既忠實(shí)條件數(shù)據(jù)點(diǎn)又易再現(xiàn)模擬目標(biāo)地質(zhì)體的幾何形態(tài)。
[0013]為讓本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂，下文特舉較佳實(shí)施例，并配合所附圖式，作詳細(xì)說明如下。
【專利附圖】

【附圖說明】
[0014]為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他的附圖。
[0015]圖1為本發(fā)明實(shí)施例提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法的流程圖；
[0016]圖2為圖1中的步驟S104的具體流程圖；
[0017]圖3為圖1中的步驟S105的具體流程圖；
[0018]圖4為圖1中的步驟S106的具體流程圖；
[0019]圖5為圖1中的步驟S107的具體流程圖；
[0020]圖6為圖5中的步驟S503的具體流程圖；
[0021]圖7為本發(fā)明實(shí)施例提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)的示意圖；
[0022]圖8為本發(fā)明提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)中三維幾何體建立裝置400的結(jié)構(gòu)框圖；
[0023]圖9為本發(fā)明提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)中空間分布體建立裝置500的結(jié)構(gòu)框圖；
[0024]圖10為本發(fā)明提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)中矢量化模型確定裝置600的結(jié)構(gòu)框圖；
[0025]圖11為本發(fā)明提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)中網(wǎng)格化裝置700的結(jié)構(gòu)框圖；
[0026]圖12為本發(fā)明提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)中數(shù)值填充模塊703的結(jié)構(gòu)框圖；
[0027]圖13為現(xiàn)有技術(shù)中的三維網(wǎng)格化示意圖；
[0028]圖14為曲流河的河曲變化形態(tài)實(shí)施方式一的示意圖；
[0029]圖15為曲流河的河曲變化形態(tài)實(shí)施方式二的示意圖；
[0030]圖16為曲流河的河曲變化形態(tài)實(shí)施方式三的示意圖；
[0031]圖17為廢棄河道的空間矢量表征方法得到的廢棄河道剖面形態(tài)示意圖；
[0032]圖18為廢棄河道的空間矢量表征方法得到的廢棄河道剖面沿著深泓線的移動(dòng)的三維幾何體示意圖；
[0033]圖19為側(cè)積層剖面的空間矢量表達(dá)方法對(duì)應(yīng)的理論反曲曲線形態(tài)示意圖；
[0034]圖20為側(cè)積層剖面的空間矢量表達(dá)方法對(duì)應(yīng)的平移反曲曲線所圈定的側(cè)積層剖面示意圖；
[0035]圖21為將側(cè)積層剖面對(duì)應(yīng)的兩個(gè)反曲曲線沿兩個(gè)圓相交的圓弧線進(jìn)行平行移動(dòng)側(cè)積層的空間分布體示 意圖；
[0036]圖22為將點(diǎn)壩砂體數(shù)據(jù)填充到圖18的廢棄河道三維幾何體的區(qū)域中得到的第一填充三維幾何體的示意圖；
[0037]圖23為矢量化的砂體儲(chǔ)層構(gòu)型模型的剖面圖；
[0038]圖24為依據(jù)水平密度填充后得到的矢量化的砂體儲(chǔ)層構(gòu)型模型的剖面圖；
[0039]圖25為網(wǎng)格化后的單個(gè)點(diǎn)壩砂體構(gòu)型模擬示意圖。
【具體實(shí)施方式】
[0040]下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例。基于本發(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0041]本發(fā)明提出一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法，圖1為該方法的具體流程圖，由圖1可知，所述的方法包括:
[0042]SlOl:獲取當(dāng)前砂體的測(cè)井?dāng)?shù)據(jù)，在具體的實(shí)施方式中，即獲取每口井在當(dāng)前建模層位的測(cè)井?dāng)?shù)據(jù)。
[0043]S102:對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)。在具體的實(shí)施方式中，對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，依據(jù)其形態(tài)確定砂體的構(gòu)型類型，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)。
[0044]S103:對(duì)所述的廢棄河道數(shù)據(jù)進(jìn)行擬合，得到平面分布的深泓線。
[0045]雖然自然界的廢棄河道形態(tài)呈現(xiàn)出各種形態(tài)，但是用數(shù)學(xué)方法去刻畫它則主要考慮到其共性，通過設(shè)置不同的參數(shù)則可以描述不同的河曲形態(tài)。本發(fā)明中對(duì)于廢棄河道的矢量表達(dá)是采用了 2部分參數(shù)。第I部分是河曲控制點(diǎn)，共為7個(gè)，可用p0、pl、……p6來表示。第2部分是河曲形態(tài)參數(shù)，可用α來表示。
[0046]對(duì)于深泓線的空間形態(tài)表征是采用了關(guān)鍵點(diǎn)控制下的三次樣條內(nèi)插法來進(jìn)行表達(dá)。對(duì)于一個(gè)河曲采用了7個(gè)關(guān)鍵控制點(diǎn)(？0，？1，...，？6)和1個(gè)形態(tài)參數(shù)角α (如圖14所示)。不同的河曲形態(tài)是通過設(shè)定形態(tài)控制關(guān)鍵點(diǎn)和形態(tài)參數(shù)角進(jìn)行差值來確定其空間各個(gè)點(diǎn)的位置，圖14和圖15中是使用不同參數(shù)得到的兩個(gè)不同形態(tài)的河曲形態(tài)例子。圖14、圖15、圖16中的曲線(即p4、p3、p2、p1、ρ6組成的曲線)即為深泓線。
[0047]S104:根據(jù)所述的深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體。廢棄河道的三維幾何體是通過沿著深泓線平移廢棄河道的剖面來表征的。廢棄河道的剖面形態(tài)中，最主要一點(diǎn)就是設(shè)定深泓線與剖面的交點(diǎn)，即河道剖面最大厚度點(diǎn)(圖17中的a3)的相對(duì)位置。a3即為深泓線上的點(diǎn)。圖2為步驟S104的具體流程圖，由圖2可知，該步驟具體包括:[0048]S201:從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取砂體厚度；
[0049]S202:將所述深泓線上的任一點(diǎn)作為數(shù)據(jù)點(diǎn)。如圖17所示，選取的數(shù)據(jù)點(diǎn)為a3，砂體厚度即為圖17中的H，即|a2a3|。
[0050]S203:根據(jù)所述的數(shù)據(jù)點(diǎn)以及曲河流點(diǎn)壩砂體經(jīng)驗(yàn)公式確定河道寬度。當(dāng)其寬深比確定后，側(cè)積層的傾角就已經(jīng)給定，而廢棄河道深度H應(yīng)該滿足寬深比的要求。此處提及的曲河流點(diǎn)壩砂體經(jīng)驗(yàn)公式包括:
[0051]①由點(diǎn)壩砂體厚度(hs)經(jīng)壓實(shí)校正求得原始砂體厚度，這一厚度也近似等于河道滿岸深度(hr)。 [0052]②依據(jù)滿岸深度由經(jīng)驗(yàn)公式lgWfl.541g hr+0.83，計(jì)算得到滿岸寬度(Wr)。
[0053]③利用公式Wm=7.44ffrL01計(jì)算出曲流帶寬度(Wm)，并依公式Wd=6.5894ffr+0.0677計(jì)算點(diǎn)壩跨度(Wd)；
[0054]④根據(jù)經(jīng)驗(yàn)規(guī)模Wl=2/3Wr計(jì)算出側(cè)積層水平寬度(W1)，并使用公式β =33exp (-0.1ffZ/HZ)計(jì)算側(cè)積層傾角(β )。
[0055]S204:根據(jù)所述的數(shù)據(jù)點(diǎn)、側(cè)積層水平寬度、河道寬度以及砂體厚度確定廢棄河道的剖面。如圖17所示，數(shù)據(jù)點(diǎn)為a3，根據(jù)側(cè)積層水平寬度、河道寬度可以確定出a1、a4，根據(jù)砂體厚度可以確定出a2，弧線aia2則可以使用以點(diǎn)&3為中心點(diǎn)的1/4橢圓的方法來刻畫其形態(tài)，其中長軸和短軸分別為Iap3I和|a2a3|。同理對(duì)于弧線a2a4也是使用以點(diǎn)a3為中心點(diǎn)的1/4橢圓的方法來刻畫其形態(tài)，其中長軸和短軸分別為|a3a4|和|a2a3|，如此即可確定出廢棄河道的剖面(如圖17所示)。
[0056]S205:遍歷所述深泓線上的點(diǎn)，得到多個(gè)廢棄河道的剖面。深泓線上的所有點(diǎn)皆可作為數(shù)據(jù)點(diǎn)構(gòu)建剖面，如此得到多個(gè)剖面。
[0057]S206:多個(gè)廢棄河道的剖面之間圈定的空間體積即為廢棄河道的三維幾何體。沿所模擬的河道深泓線在不同的位置計(jì)算河道空間的剖面形態(tài)后，各個(gè)面之間所圈定的空間體積即為廢棄河道的三維幾何體，如圖18所示。在實(shí)際模擬中，河道深泓線的參數(shù)以及剖面形態(tài)等參數(shù)都可以依據(jù)地質(zhì)家對(duì)于河曲形態(tài)的理解和研究區(qū)的具體問題進(jìn)行設(shè)置。
[0058]由圖1可知，該方法還包括:
[0059]S105:根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體。圖3為步驟S105的具體流程圖，由圖3可知，該步驟具體包括:
[0060]S301:根據(jù)反曲曲線公式確定在實(shí)數(shù)范圍內(nèi)的反曲曲線形態(tài)。
[0061]曲流河構(gòu)型模式中另一個(gè)重要的構(gòu)型要素則是側(cè)積層。在本發(fā)明中對(duì)于側(cè)積層剖面形態(tài)描述采用了反曲曲線。反曲曲線的數(shù)學(xué)表達(dá)為

I
[0062]J ~—

1 + e
[0063]其在實(shí)數(shù)范圍內(nèi)的形狀為圖19所示。
[0064]S302:對(duì)所述的反曲曲線形態(tài)進(jìn)行坐標(biāo)平移或坐標(biāo)旋轉(zhuǎn)，得到側(cè)積層剖面。在實(shí)際模擬中，可采用坐標(biāo)平移或者坐標(biāo)旋轉(zhuǎn)的方法生成形態(tài)各異的側(cè)積層剖面形態(tài)，如圖20所示，是對(duì)一個(gè)反曲曲線進(jìn)行平移所圈定的剖面，這個(gè)剖面構(gòu)成了側(cè)積層剖面。側(cè)積層剖面包括著兩個(gè)反曲曲線。
[0065]S303:從所述的側(cè)積層數(shù)據(jù)中獲取側(cè)積層中心厚度、側(cè)積層水平跨度。[0066]S304:設(shè)置兩個(gè)圓心距為側(cè)積層中心厚度、直徑為側(cè)積層水平跨度的圓，如圖21中所示的兩個(gè)圓，兩個(gè)圓的圓心的距離Ic1C2I等于側(cè)積層中心的厚度，兩個(gè)圓的直徑為側(cè)積層水平跨度，而延展的距離則由圖21中的Iap2I和|a2a3|來控制，則兩個(gè)圓可以形成如圖21所示的兩個(gè)相交的圓弧線。
[0067]S305:將所述側(cè)積層剖面對(duì)應(yīng)的兩個(gè)反曲曲線沿兩個(gè)圓相交的圓弧線進(jìn)行平行移動(dòng)；
[0068]S306:平行移動(dòng)所圈定的空間即為側(cè)積層的空間分布體。如圖21所示。模擬過程中是通過給定圖21中的|aia2|和|a2a3|的分布函數(shù)從而控制側(cè)積層側(cè)向延伸的距離和延伸的方式。將兩個(gè)反曲曲線沿著所定義好的圓弧線進(jìn)行平行移動(dòng)所圈定的空間體積即是側(cè)積層的空間分布體積。
[0069]如上所述，本發(fā)明雖然模擬仍然采用了傳統(tǒng)的基于目標(biāo)體的模擬，但是并不設(shè)定模型的網(wǎng)格尺寸，而是直接在模型內(nèi)用面和體的方法生成模擬地質(zhì)體。模型只存儲(chǔ)地質(zhì)體的空間分布形態(tài)參數(shù)而不存儲(chǔ)每個(gè)網(wǎng)格上的目標(biāo)體的屬性。故而必須用數(shù)學(xué)公式或形態(tài)參數(shù)刻畫模型中的基本地質(zhì)體。雖然地下的地質(zhì)體形態(tài)是極其復(fù)雜的，但是這些復(fù)雜的地質(zhì)體都可以用一些簡(jiǎn)單形狀比如帶狀、球狀、橢球狀，楔形等經(jīng)過一定的修正和調(diào)整來刻畫。其中地質(zhì)學(xué)家的儲(chǔ)層構(gòu)型概念地質(zhì)模式是定義模擬基本單元的基礎(chǔ)。故而本發(fā)明以現(xiàn)在的構(gòu)型地質(zhì)研究作為基礎(chǔ)，以構(gòu)型建模算法的開發(fā)為目標(biāo)，對(duì)各類儲(chǔ)層構(gòu)型要素進(jìn)行三維形態(tài)刻畫和表征。
[0070]由圖1可知，該方法還包括:
[0071]S106:根根據(jù)所述廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型。圖4為步驟S106的具體流程圖，由圖4可知，該步驟具體包括:
[0072]S401:將所述的點(diǎn)壩砂體數(shù)據(jù)填充到所述廢棄河道三維幾何體所圈定的區(qū)域中，得到第一填充三維幾何體，如圖22所示，將點(diǎn)壩砂體數(shù)據(jù)填充到圖18的廢棄河道三維幾何體的區(qū)域中得到第一填充三維幾何體。
[0073]S402:從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層在井身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度；
[0074]S403:根據(jù)所述側(cè)積層在井身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度，將所述側(cè)積層的空間分布體填充到第一填充三維幾何體，得到第二填充幾何體。側(cè)積層在測(cè)井上的位置用于表征此口井鉆遇側(cè)積層的位置(如圖23示例中的LI和L2)，即必須于井點(diǎn)的構(gòu)型要素解釋類型相吻合，角度用于表征側(cè)積層發(fā)育的角度(見圖23中示例)。如圖23所示，將側(cè)積層的空間分布體分別填充到井A、井B，與點(diǎn)壩砂體表面相交的點(diǎn)為a、C。
[0075]S404:從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層的水平密度，水平密度用于表征側(cè)積層發(fā)育的多少。
[0076]S405:根據(jù)所述的水平密度將將所述側(cè)積層的空間分布體填充到所述第二填充幾何體，得到矢量化的砂體儲(chǔ)層構(gòu)型模型。如圖24所示，將側(cè)積層的空間分布體分別填充到點(diǎn)壩砂體中，與點(diǎn)壩砂體表面的交點(diǎn)為b、d。圖25即為網(wǎng)格化后的單個(gè)點(diǎn)壩砂體構(gòu)型模擬示意圖。
[0077]上述步驟利用所定義的基本構(gòu)型要素地質(zhì)體進(jìn)行隨機(jī)模擬。地下儲(chǔ)層構(gòu)型表征的核心是井間預(yù)測(cè)，而預(yù)測(cè)的關(guān)鍵是預(yù)知對(duì)象的分布規(guī)律或模式。雖然模型的調(diào)節(jié)只是對(duì)各個(gè)地質(zhì)體的表征參數(shù)進(jìn)行調(diào)節(jié)，但是約束依然是各級(jí)儲(chǔ)層構(gòu)型定量地質(zhì)模式。如各個(gè)構(gòu)型要素的展布規(guī)模和相互之間的切割和疊置關(guān)系。國內(nèi)外進(jìn)行過不少的儲(chǔ)層地質(zhì)知識(shí)庫的研究，初步建立了垂向信息與側(cè)向規(guī)模的關(guān)系，涉及了河流相、三角洲相、扇三角洲相、深水濁積相等，取得了豐富的研究成果，這些研究都是進(jìn)行隨機(jī)模擬的基礎(chǔ)。故而本發(fā)明將現(xiàn)在的構(gòu)型地質(zhì)定量模式研究作為基礎(chǔ)，以構(gòu)型隨機(jī)模擬算法的開發(fā)為目標(biāo)，對(duì)各類儲(chǔ)層構(gòu)型要素的三維空間隨機(jī)模擬方法進(jìn)行研究，最終得到基于矢量表達(dá)的儲(chǔ)層構(gòu)型模型。
[0078]在定義好各個(gè)構(gòu)型要素之后則是進(jìn)行基于目標(biāo)的隨機(jī)模擬。與傳統(tǒng)的基于目標(biāo)體的模擬方法一樣，認(rèn)為目標(biāo)中心點(diǎn)位置符合齊次泊松點(diǎn)過程，而在無井區(qū)，模擬采用非齊次的泊松點(diǎn)過程，從而滿足了井間與井點(diǎn)分布具有差異的要求。模擬仍然采用迭代的算法，并且在模擬過程的必須受到各個(gè)目標(biāo)體的疊置切割模式和井?dāng)?shù)據(jù)的約束。
[0079]地質(zhì)構(gòu)型分析所得到的單井構(gòu)型解釋是作為模擬階段的硬數(shù)據(jù)，最終的構(gòu)型模型必須滿足井?dāng)?shù)據(jù)的約束。本發(fā)明在模擬程序中采用了優(yōu)先擬合井點(diǎn)數(shù)據(jù)的做法。首先通過調(diào)整各個(gè)地質(zhì)體的形態(tài)參數(shù)使其滿足井點(diǎn)約束然。之后在選定的地質(zhì)模式約束下進(jìn)行空間其它構(gòu)型要素的投放以滿足目標(biāo)函數(shù)。在點(diǎn)壩砂體內(nèi)部模擬側(cè)積層為例(如圖24)，依據(jù)井?dāng)?shù)據(jù)模擬出了側(cè)積層a和C。而依據(jù)側(cè)積層間距等定量構(gòu)型模式信息可以模擬出側(cè)積層b和d。
[0080]由圖1可知，該方法還包括:
[0081]S107:對(duì)所述的矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。圖5為步驟S107的具體流程圖，由圖5可知，該步驟包括:
[0082]S501:根據(jù)所述的側(cè)積層數(shù)據(jù)確定網(wǎng)格尺寸，即網(wǎng)格必須滿足刻畫側(cè)積層的最小厚度。在具體的實(shí)施方式中，如側(cè)積層數(shù)據(jù)為0.8米厚，則確定的網(wǎng)格尺寸可以是0.5米。
[0083]S502:根據(jù)所述的網(wǎng)格尺寸對(duì)所述矢量化的砂體儲(chǔ)層構(gòu)型模型建立模擬區(qū)的網(wǎng)格；
[0084]S503:對(duì)所述模擬區(qū)的網(wǎng)格根據(jù)構(gòu)型要素類型進(jìn)行數(shù)值填充，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。圖6為步驟S503的具體流程圖，由圖6可知，該步驟具體包括:
[0085]S601:判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為廢棄河道；
[0086]S602:當(dāng)判斷為是時(shí)，將所述的網(wǎng)格填充為廢棄河道數(shù)據(jù)；
[0087]S603:否則，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為側(cè)積層；
[0088]S604:當(dāng)判斷為是時(shí)，將所述的網(wǎng)格填充為側(cè)積層數(shù)據(jù)；
[0089]S605:否則，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為點(diǎn)壩砂體；
[0090]S606:當(dāng)判斷為是時(shí)，將所述的網(wǎng)格填充為點(diǎn)壩砂體數(shù)據(jù)；
[0091]S607:否則，確定所述網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型的比例。當(dāng)網(wǎng)格中包括多種構(gòu)型要素類型時(shí)，分別確定各個(gè)要素所占的比例，如:網(wǎng)格中既有點(diǎn)壩砂體，又有側(cè)積層，則確定二者的比例，為2:5。
[0092]S608:獲取所述比例中的最大值。
[0093]S609:確定所述最大值對(duì)應(yīng)的構(gòu)型要素類型；
[0094]S610:將所述的網(wǎng)格填充為最大值對(duì)應(yīng)的構(gòu)型要素類型對(duì)應(yīng)的數(shù)據(jù)。
[0095]如上的實(shí)施例中，網(wǎng)格中既有點(diǎn)壩砂體，又有側(cè)積層，二者的比例為2:5，則比例最大之為5，對(duì)應(yīng)的構(gòu)型要素類型為側(cè)積層，因此該網(wǎng)格填充的數(shù)據(jù)為側(cè)積層數(shù)據(jù)。如圖25為單個(gè)點(diǎn)壩砂體構(gòu)型模型網(wǎng)格化結(jié)果。
[0096]由于后期物性參數(shù)的建模方法仍然是基于網(wǎng)格的，而且后期的數(shù)值模擬也是基于離散網(wǎng)格數(shù)據(jù)體，故而模擬的最終結(jié)果應(yīng)該是一個(gè)基于網(wǎng)格的模型。所以必須將矢量存儲(chǔ)的空間地質(zhì)體進(jìn)行離散化，將各個(gè)儲(chǔ)層構(gòu)型屬性賦于每一個(gè)網(wǎng)格。由于基于矢量存儲(chǔ)的地質(zhì)模型不受空間網(wǎng)格大小的限制，并且不再考慮模擬速度的要求，此階段的網(wǎng)格定義可以采取任意細(xì)化以滿足刻畫不同規(guī)模的非均質(zhì)地質(zhì)體。故而如何對(duì)基于矢量表達(dá)的模型進(jìn)行網(wǎng)格化處理也是本發(fā)明的內(nèi)容之一。由于此階段不早模擬，沒有迭代等計(jì)算，所以可以選擇能夠刻畫最小規(guī)模的構(gòu)型要素的模型網(wǎng)格進(jìn)行離散化，如圖25為單個(gè)點(diǎn)壩砂體構(gòu)型模型網(wǎng)格化結(jié)果。
[0097]如上即為本發(fā)明提供的方法，本發(fā)明還提出一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)，圖7為該系統(tǒng)的結(jié)構(gòu)示意圖，由圖7可知，所述的系統(tǒng)包括:
[0098]測(cè)井?dāng)?shù)據(jù)獲取裝置100，用于獲取當(dāng)前砂體的測(cè)井?dāng)?shù)據(jù)，在具體的實(shí)施方式中，SP獲取每口井在當(dāng)前建模層位的測(cè)井?dāng)?shù)據(jù)。
[0099]測(cè)井?dāng)?shù)據(jù)處理裝置200，用于對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)。在具體的實(shí)施方式中，對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，依據(jù)其形態(tài)確定砂體的構(gòu)型類型，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)。此處的數(shù)字處理具體包括:依據(jù)測(cè)井相進(jìn)行。
[0100]擬合裝置300，用于對(duì)所述的廢棄河道數(shù)據(jù)進(jìn)行擬合，得到平面分布的深泓線。
[0101]雖然自然界的廢棄河道形態(tài)呈現(xiàn)出各種形態(tài)，但是用數(shù)學(xué)方法去刻畫它則主要考慮到其共性，通過設(shè)置不同的參數(shù)則可以描述不同的河曲形態(tài)。本發(fā)明中對(duì)于廢棄河道的矢量表達(dá)是采用了 2部分參數(shù)。第I部分是河`曲控制點(diǎn)，共為7個(gè)，可用p0、pl、……p6來表示。第2部分是河曲形態(tài)參數(shù)，可用α來表示。
[0102]對(duì)于深泓線的空間形態(tài)表征是采用了關(guān)鍵點(diǎn)控制下的三次樣條內(nèi)插法來進(jìn)行表達(dá)。對(duì)于一個(gè)河曲采用了7個(gè)關(guān)鍵控制點(diǎn)(？0，？1，...，？6)和1個(gè)形態(tài)參數(shù)角α (如圖14所示)。不同的河曲形態(tài)是通過設(shè)定形態(tài)控制關(guān)鍵點(diǎn)和形態(tài)參數(shù)角進(jìn)行差值來確定其空間各個(gè)點(diǎn)的位置，圖14和圖15中是使用不同參數(shù)得到的兩個(gè)不同形態(tài)的河曲形態(tài)例子。圖14、圖15、圖16中的曲線(即p4、p3、p2、p1、ρ6組成的曲線)即為深泓線。
[0103]三維幾何體建立裝置400，用于根據(jù)所述的深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體。廢棄河道的三維幾何體是通過沿著深泓線平移廢棄河道的剖面來表征的。廢棄河道的剖面形態(tài)中，最主要一點(diǎn)就是設(shè)定深泓線與剖面的交點(diǎn)，即河道剖面最大厚度點(diǎn)(圖17中的a3)的相對(duì)位置。a3即為深泓線上的點(diǎn)。圖8為三維幾何體建立裝置400的結(jié)構(gòu)框圖，由圖8可知，三維幾何體建立裝置具體包括:
[0104]砂體厚度獲取模塊401，用于從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取砂體厚度；
[0105]數(shù)據(jù)點(diǎn)選取模塊402，用于將所述深泓線上的任一點(diǎn)作為數(shù)據(jù)點(diǎn)。如圖17所示，選取的數(shù)據(jù)點(diǎn)為a3，砂體厚度即為圖17中的H，即|a2a3|。
[0106]河道寬度確定模塊403，用于根據(jù)所述的數(shù)據(jù)點(diǎn)以及曲河流點(diǎn)壩砂體經(jīng)驗(yàn)公式確定凹岸水平寬度以及河道寬度。
[0107]當(dāng)其寬深比確定后，側(cè)積層的傾角就已經(jīng)給定，而廢棄河道深度H應(yīng)該滿足寬深比的要求。此處提及的曲河流點(diǎn)壩砂體經(jīng)驗(yàn)公式包括:
[0108]①由點(diǎn)壩砂體厚度(hs)經(jīng)壓實(shí)校正求得原始砂體厚度，這一厚度也近似等于河道滿岸深度(hr)。
[0109]②依據(jù)滿岸深度由經(jīng)驗(yàn)公式lgWfl.541g hr+0.83，計(jì)算得到滿岸寬度(Wr)。
[0110]③利用公式Wm=7.44ffrL01計(jì)算出曲流帶寬度(Wm)，并依公式Wd=6.5894ffr+0.0677計(jì)算點(diǎn)壩跨度(Wd)；
[0111]④根據(jù)經(jīng)驗(yàn)規(guī)模Wl=2/3Wr計(jì)算出側(cè)積層水平寬度(W1)，并使用公式β =33exp (-0.1ffZ/HZ)計(jì)算側(cè)積層傾角(β )。
[0112]廢棄河道剖面確定模塊404，用于根據(jù)所述的數(shù)據(jù)點(diǎn)、側(cè)積層水平寬度、河道寬度以及砂體厚度確定廢棄河道的剖面。如圖17所示，數(shù)據(jù)點(diǎn)為a3，根據(jù)側(cè)積層水平寬度、河道寬度可以確定出apa4,根據(jù)砂體厚度可以確定出a2,弧線則可以使用以點(diǎn)a3為中心點(diǎn)的1/4橢圓的方法來刻畫其形態(tài)，其中長軸和短軸分別為Iap3I和|a2a3|。同理對(duì)于弧線a2a4也是使用以點(diǎn)a3為中心點(diǎn)的1/4橢圓的方法來刻畫其形態(tài)，其中長軸和短軸分別為
a3a4|和|a2a3|， 如此即可確定出廢棄河道的剖面(如圖17所示)。
[0113]數(shù)據(jù)點(diǎn)遍歷模塊405，用于遍歷所述深泓線上的點(diǎn)，得到多個(gè)廢棄河道的剖面。深泓線上的所有點(diǎn)皆可作為數(shù)據(jù)點(diǎn)構(gòu)建剖面，如此得到多個(gè)剖面。
[0114]三維幾何體確定模塊406，多個(gè)廢棄河道的剖面之間圈定的空間體積即為廢棄河道的三維幾何體。沿所模擬的河道深泓線在不同的位置計(jì)算河道空間的剖面形態(tài)后，各個(gè)面之間所圈定的空間體積即為廢棄河道的三維幾何體，如圖18所示。在實(shí)際模擬中，河道深泓線的參數(shù)以及剖面形態(tài)等參數(shù)都可以依據(jù)地質(zhì)家對(duì)于河曲形態(tài)的理解和研究區(qū)的具體問題進(jìn)行設(shè)置。
[0115]由圖7可知，該系統(tǒng)還包括:
[0116]空間分布體建立裝置500，用于根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體。圖9為空間分布體建立裝置500的結(jié)構(gòu)框圖，由圖9可知，空間分布體建立裝置具體包括:
[0117]反曲曲線形態(tài)確定模塊501，用于根據(jù)反曲曲線公式確定在實(shí)數(shù)范圍內(nèi)的反曲曲線形態(tài)。
[0118]曲流河構(gòu)型模式中另一個(gè)重要的構(gòu)型要素則是側(cè)積層。在本發(fā)明中對(duì)于側(cè)積層剖面形態(tài)描述采用了反曲曲線。反曲曲線的數(shù)學(xué)表達(dá)為
[0119]y=^~7

1 + e
[0120]其在實(shí)數(shù)范圍內(nèi)的形狀為圖19所示。
[0121]側(cè)積層剖面確定模塊502，用于對(duì)所述的反曲曲線形態(tài)進(jìn)行坐標(biāo)平移或坐標(biāo)旋轉(zhuǎn)，得到側(cè)積層剖面。在實(shí)際模擬中，可采用坐標(biāo)平移或者坐標(biāo)旋轉(zhuǎn)的方法生成形態(tài)各異的側(cè)積層剖面形態(tài)，如圖20所示，是對(duì)一個(gè)反曲曲線進(jìn)行平移所圈定的剖面，這個(gè)剖面構(gòu)成了側(cè)積層剖面。側(cè)積層剖面包括著兩個(gè)反曲曲線。
[0122]水平跨度獲取模塊503，用于從所述的側(cè)積層數(shù)據(jù)中獲取側(cè)積層中心厚度、側(cè)積層水平跨度。
[0123]設(shè)置模塊504，用于設(shè)置兩個(gè)圓心距為側(cè)積層中心厚度、直徑為側(cè)積層水平跨度的圓，
[0124]如圖21中所示的兩個(gè)圓，兩個(gè)圓的圓心的距離Ic1C2I等于側(cè)積層中心的厚度，兩個(gè)圓的直徑為側(cè)積層水平跨度，而延展的距離則由圖21中的|aia2|和|a2a3|來控制，則兩個(gè)圓可以形成如圖21所示的兩個(gè)相交的圓弧線。
[0125]平行移動(dòng)模塊505，用于將所述側(cè)積層剖面對(duì)應(yīng)的兩個(gè)反曲曲線沿兩個(gè)圓相交的圓弧線進(jìn)行平行移動(dòng)；
[0126]空間分布體確定模塊506，平行移動(dòng)所圈定的空間即為側(cè)積層的空間分布體。如圖21所示。模擬過程中是通過給定圖21中的|aia2|和|a2a3|的分布函數(shù)從而控制側(cè)積層側(cè)向延伸的距離和延伸的方式。將兩個(gè)反曲曲線沿著所定義好的圓弧線進(jìn)行平行移動(dòng)所圈定的空間體積即是側(cè)積層的空間分布體積。
[0127]如上所述，本發(fā)明雖然模擬仍然采用了傳統(tǒng)的基于目標(biāo)體的模擬，但是并不設(shè)定模型的網(wǎng)格尺寸，而是直接在模型內(nèi)用面和體的方法生成模擬地質(zhì)體。模型只存儲(chǔ)地質(zhì)體的空間分布形態(tài)參數(shù)而不存儲(chǔ)每個(gè)網(wǎng)格上的目標(biāo)體的屬性。故而必須用數(shù)學(xué)公式或形態(tài)參數(shù)刻畫模型中的基本地質(zhì)體。雖然地下的地質(zhì)體形態(tài)是極其復(fù)雜的，但是這些復(fù)雜的地質(zhì)體都可以用一些簡(jiǎn)單形狀比如帶狀、球狀、橢球狀，楔形等經(jīng)過一定的修正和調(diào)整來刻畫。其中地質(zhì)學(xué)家的儲(chǔ)層構(gòu)型概念地質(zhì)模式是定義模擬基本單元的基礎(chǔ)。故而本發(fā)明以現(xiàn)在的構(gòu)型地質(zhì)研究作為基礎(chǔ)，以構(gòu)型建模算法的開發(fā)為目標(biāo)，對(duì)各類儲(chǔ)層構(gòu)型要素進(jìn)行三維形態(tài)刻畫和表征。由圖7可知，該系統(tǒng)還包括:
[0128]矢量化模型確定裝置600，用于根根據(jù)所述廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型。圖10為本發(fā)明提供的一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)中矢量化模型確定裝置600的結(jié)構(gòu)框圖，由圖10可知，矢量化模型確定裝置具體包括:
[0129]第一填充模塊601，用于將所述的點(diǎn)壩砂體數(shù)據(jù)填充到所述廢棄河道三維幾何體的所圈定的區(qū)域中，得到第一填充三維幾何體，如圖22所示，將點(diǎn)壩砂體數(shù)據(jù)填充到圖18的廢棄河道三維幾何體的區(qū)域中得到第一填充三維幾何體。
[0130]角度獲取模塊602，用于從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層在井身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度；
[0131]第二填充模塊603，用于根據(jù)所述側(cè)積層在井身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度，將所述側(cè)積層的空間分布體填充到第一填充三維幾何體，得到第二填充幾何體。側(cè)積層在測(cè)井上的位置用于表征此口井鉆遇側(cè)積層的位置(如圖23示例中的LI和L2)，即必須于井點(diǎn)的構(gòu)型要素解釋類型相吻合，角度用于表征側(cè)積層發(fā)育的角度(見圖23中示例)。如圖23所示，將側(cè)積層的空間分布體分別填充到井A、井B，與點(diǎn)壩砂體表面相交的點(diǎn)為a、C。
[0132]水平密度獲取模塊604，用于從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層的水平密度，水平密度用于表征在點(diǎn)壩砂體內(nèi)部非井點(diǎn)控制處側(cè)積層發(fā)育的可能性。
[0133]矢量化模型確定模塊605，用于根據(jù)所述的水平密度將將所述側(cè)積層的空間分布體填充到所述第二填充幾何體，得到矢量化的砂體儲(chǔ)層構(gòu)型模型。如圖24所示，將側(cè)積層的空間分布體分別填充到點(diǎn)壩砂體中，與點(diǎn)壩砂體表面的交點(diǎn)為b、d。圖25即為網(wǎng)格化后的單個(gè)點(diǎn)壩砂體構(gòu)型模擬示意圖。[0134]上述步驟利用所定義的基本構(gòu)型要素地質(zhì)體進(jìn)行隨機(jī)模擬。地下儲(chǔ)層構(gòu)型表征的核心是井間預(yù)測(cè)，而預(yù)測(cè)的關(guān)鍵是預(yù)知對(duì)象的分布規(guī)律或模式。雖然模型的調(diào)節(jié)只是對(duì)各個(gè)地質(zhì)體的表征參數(shù)進(jìn)行調(diào)節(jié)，但是約束依然是各級(jí)儲(chǔ)層構(gòu)型定量地質(zhì)模式。如各個(gè)構(gòu)型要素的展布規(guī)模和相互之間的切割和疊置關(guān)系。國內(nèi)外進(jìn)行過不少的儲(chǔ)層地質(zhì)知識(shí)庫的研究，初步建立了垂向信息與側(cè)向規(guī)模的關(guān)系，涉及了河流相、三角洲相、扇三角洲相、深水濁積相等，取得了豐富的研究成果，這些研究都是進(jìn)行隨機(jī)模擬的基礎(chǔ)。故而本發(fā)明將現(xiàn)在的構(gòu)型地質(zhì)定量模式研究作為基礎(chǔ)，以構(gòu)型隨機(jī)模擬算法的開發(fā)為目標(biāo)，對(duì)各類儲(chǔ)層構(gòu)型要素的三維空間隨機(jī)模擬方法進(jìn)行研究，最終得到基于矢量表達(dá)的儲(chǔ)層構(gòu)型模型。
[0135]在定義好各個(gè)構(gòu)型要素之后則是進(jìn)行基于目標(biāo)的隨機(jī)模擬。與傳統(tǒng)的基于目標(biāo)體的模擬方法一樣，認(rèn)為目標(biāo)中心點(diǎn)位置符合齊次泊松點(diǎn)過程，而在無井區(qū)，模擬采用非齊次的泊松點(diǎn)過程，從而滿足了井間與井點(diǎn)分布具有差異的要求。模擬仍然采用迭代的算法，并且在模擬過程的必須受到各個(gè)目標(biāo)體的疊置切割模式和井?dāng)?shù)據(jù)的約束。
[0136]地質(zhì)構(gòu)型分析所得到的單井構(gòu)型解釋是作為模擬階段的硬數(shù)據(jù)，最終的構(gòu)型模型必須滿足井?dāng)?shù)據(jù)的約束。本發(fā)明在模擬程序中采用了優(yōu)先擬合井點(diǎn)數(shù)據(jù)的做法。首先通過調(diào)整各個(gè)地質(zhì)體的形態(tài)參數(shù)使其滿足井點(diǎn)約束然。之后在選定的地質(zhì)模式約束下進(jìn)行空間其它構(gòu)型要素的投放以滿足目標(biāo)函數(shù)。在點(diǎn)壩砂體內(nèi)部模擬側(cè)積層為例(如圖24)，依據(jù)井?dāng)?shù)據(jù)模擬出了側(cè)積層a和C。而依據(jù)側(cè)積層間距等定量構(gòu)型模式信息可以模擬出側(cè)積層b和d。
[0137]由圖7可知，該系統(tǒng)還包括:
[0138]網(wǎng)格化裝置700，用于對(duì)所述的矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。圖11為網(wǎng)格化裝置700的結(jié)構(gòu)框圖，由圖11可知，網(wǎng)格化裝置包括:
[0139]網(wǎng)格尺寸確定模塊701，用于根據(jù)所述的側(cè)積層數(shù)據(jù)確定網(wǎng)格尺寸。在具體的實(shí)施方式中，如側(cè)積層數(shù)據(jù)為垂向上最厚處為0.8米，則確定的網(wǎng)格尺寸為0.5米。
[0140]網(wǎng)格建立模塊702，用于根據(jù)所述的網(wǎng)格尺寸對(duì)所述矢量化的砂體儲(chǔ)層構(gòu)型模型建立模擬區(qū)的網(wǎng)格；
[0141]數(shù)值填充模塊703，用于對(duì)所述模擬區(qū)的網(wǎng)格根據(jù)構(gòu)型要素類型進(jìn)行數(shù)值填充，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。圖12為數(shù)值填充模塊703的結(jié)構(gòu)框圖，由圖12可知，數(shù)值填充模塊具體包括:
[0142]第一判斷單元7031，用于判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為廢棄河道；
[0143]第一填充單元7032，用于當(dāng)所述的第一判斷單元判斷為是時(shí)，將所述的網(wǎng)格填充為廢棄河道數(shù)據(jù)；
[0144]第二判斷單元7033，用于當(dāng)所述的第一判斷單元判斷為否時(shí)，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為側(cè)積層；
[0145]第二填充單元7034，用于當(dāng)所述的第二判斷單元判斷為是時(shí)，將所述的網(wǎng)格填充為側(cè)積層數(shù)據(jù)；
[0146]第三判斷單元7035，用于當(dāng)所述的第二判斷單元判斷為否時(shí)，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為點(diǎn)壩砂體；[0147]第三填充單元7036，用于當(dāng)所述的第三判斷單元判斷為是時(shí)，將所述的網(wǎng)格填充為點(diǎn)壩砂體數(shù)據(jù)；
[0148]比例確定單元7037，用于當(dāng)所述第三判斷單元判斷為否時(shí)，確定所述網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型的比例。當(dāng)網(wǎng)格中包括多種構(gòu)型要素類型時(shí)，分別確定各個(gè)要素所占的比例，如:網(wǎng)格中既有點(diǎn)壩砂體，又有側(cè)積層，則確定二者的比例，為2:5。
[0149]最大值獲取單元7038，用于獲取所述比例中的最大值。
[0150]要素類型確定單元7039，用于確定所述最大值對(duì)應(yīng)的構(gòu)型要素類型；
[0151]第四填充單元7040，用于將所述的網(wǎng)格填充為最大值對(duì)應(yīng)的構(gòu)型要素類型對(duì)應(yīng)的數(shù)據(jù)。
[0152]如上的實(shí)施例中，網(wǎng)格中既有點(diǎn)壩砂體，又有側(cè)積層，二者的比例為2:5，則比例最大之為5，對(duì)應(yīng)的構(gòu)型要素類型為側(cè)積層，因此該網(wǎng)格填充的數(shù)據(jù)為側(cè)積層數(shù)據(jù)。如圖25為單個(gè)點(diǎn)壩砂體構(gòu)型模型網(wǎng)格化結(jié)果。
[0153]由于后期物性參數(shù)的建模方法仍然是基于網(wǎng)格的，而且后期的數(shù)值模擬也是基于離散網(wǎng)格數(shù)據(jù)體，故而模擬的最終結(jié)果應(yīng)該是一個(gè)基于網(wǎng)格的模型。所以必須將矢量存儲(chǔ)的空間地質(zhì)體進(jìn)行離散化，將各個(gè)儲(chǔ)層構(gòu)型屬性賦于每一個(gè)網(wǎng)格。由于基于矢量存儲(chǔ)的地質(zhì)模型不受空間網(wǎng)格大小的限制，并且不再考慮模擬速度的要求，此階段的網(wǎng)格定義可以采取任意細(xì)化以滿足刻畫不同規(guī)模的非均質(zhì)地質(zhì)體。故而如何對(duì)基于矢量表達(dá)的模型進(jìn)行網(wǎng)格化處理也是本發(fā)明的內(nèi)容之一。由于此階段不早模擬，沒有迭代等計(jì)算，所以可以選擇能夠刻畫最小規(guī)模的構(gòu)型要素的模型網(wǎng)格進(jìn)行離散化，如圖25為單個(gè)點(diǎn)壩砂體構(gòu)型模型網(wǎng)格化結(jié)果。
[0154]下面結(jié)合具體的實(shí)施例，詳細(xì)介紹本發(fā)明的技術(shù)方案。研究區(qū)井距最大約380m，最小約25m，平均井距IOOm左右。研究區(qū)目的層位為館陶組上段，其細(xì)分為5個(gè)小層，其中本文選取S3小層為例。S3小層發(fā)育點(diǎn)壩、天然堤、決口扇等微相。巖性自下而上為中、細(xì)粒度砂巖漸變?yōu)榧?xì)、粉砂巖至泥巖具典型正韻律特征，底部有沖刷面，河流二元結(jié)構(gòu)特征明顯。
[0155](I)讀取工區(qū)85 口井的自然電位、自然伽馬、微電極曲線等測(cè)井曲線數(shù)據(jù)、分層數(shù)據(jù)。
[0156](2)依據(jù)測(cè)井曲線特征，確定S3小層各井的儲(chǔ)層構(gòu)型類型。如曲流河河道砂體中的泥質(zhì)側(cè)積層的巖性主要包括泥巖、頁巖、粉砂質(zhì)泥巖及部分泥質(zhì)粉砂巖。巖電標(biāo)定結(jié)果顯示研究區(qū)內(nèi)泥質(zhì)側(cè)積層較薄，一般為0.2?0.Sm，測(cè)井曲線表現(xiàn)出微電極曲線明顯回返，幅度差減小，自然伽馬曲線見回返，自然電位曲線輕微回返或者不明顯的回返。過工區(qū)的一條沉積相剖面連井對(duì)比分析。通過以上地質(zhì)分析所解釋的單井各級(jí)構(gòu)型要素將作為構(gòu)型模擬的條件數(shù)據(jù)。
[0157](3)對(duì)于本工區(qū)，根據(jù)經(jīng)驗(yàn)公式和井?dāng)?shù)據(jù)，推斷單河道河流滿岸深度為7.2m，平均河流滿岸寬度162m，單一側(cè)積體水平寬度約為106m左右，側(cè)積層的傾角約為5?10°。
[0158](4)得到本工區(qū)中使用本發(fā)明得到的單個(gè)點(diǎn)壩砂體構(gòu)型模擬示意圖，模型較好的再現(xiàn)了點(diǎn)壩內(nèi)部及點(diǎn)壩內(nèi)部側(cè)積層的空間分布特征。
[0159]綜上所述，本發(fā)明提供了一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法及系統(tǒng)，采用的基于空間矢量的儲(chǔ)層構(gòu)型建模方法可以有效解決現(xiàn)有技術(shù)存在的問題。該方案提出的主要目的是在儲(chǔ)層構(gòu)型建模中描述不同級(jí)別的非均質(zhì)體以及準(zhǔn)確刻畫各級(jí)地質(zhì)體的復(fù)雜空間形態(tài)。在模擬階段的地質(zhì)體的存儲(chǔ)借鑒了矢量圖形存儲(chǔ)的思路，即對(duì)模擬地質(zhì)體的定義采用了基于實(shí)數(shù)集的點(diǎn)、線、面所限定的體。因而模型在模擬過程中不受網(wǎng)格尺寸的限制，使得模擬結(jié)果可以刻畫不同規(guī)模層次的非均質(zhì)體。應(yīng)用基于矢量體存儲(chǔ)方法的另外一個(gè)優(yōu)點(diǎn)是可以快速的對(duì)模型進(jìn)行調(diào)節(jié)，因?yàn)樵谀M過程中只是對(duì)地質(zhì)體的空間位置和形態(tài)的定義參數(shù)進(jìn)行調(diào)節(jié)，因而這樣的調(diào)節(jié)方式較基于網(wǎng)格的調(diào)節(jié)更便捷，從而使結(jié)果既忠實(shí)條件數(shù)據(jù)點(diǎn)又易再現(xiàn)模擬目標(biāo)地質(zhì)體的幾何形態(tài)。
[0160]該方案依據(jù)地質(zhì)家對(duì)研究工區(qū)的地質(zhì)模式定義模型中的構(gòu)型要素。用構(gòu)型要素形態(tài)參數(shù)和數(shù)學(xué)關(guān)系來描述其空間形態(tài)及彼此之間的相對(duì)位置關(guān)系。
[0161]該方案利用所定義的構(gòu)型要素進(jìn)行隨機(jī)模擬。雖然仍然采用了傳統(tǒng)的基于目標(biāo)體的模擬，但是并不設(shè)定模型的網(wǎng)格尺寸，而是直接在模型內(nèi)用面和體的方法生成構(gòu)型模型。模型只存儲(chǔ)構(gòu)型要素空間分布形態(tài)參數(shù)。并且模型的條件化過程僅對(duì)構(gòu)型要素的形態(tài)表征參數(shù)和位置關(guān)系參數(shù)進(jìn)行調(diào)節(jié)，這樣避免了基于網(wǎng)格的方法中網(wǎng)格個(gè)數(shù)增加帶來的計(jì)算時(shí)間挑戰(zhàn)。
[0162]該方案對(duì)模擬的結(jié)果進(jìn)行網(wǎng)格化。由于不再考慮模擬速度的要求，此階段可以采取任意細(xì)化的網(wǎng)格以滿足刻畫所關(guān)注的關(guān)鍵構(gòu)型要素，如點(diǎn)壩砂體中的側(cè)積層。
[0163]本領(lǐng)域普通技術(shù)人員可以理解實(shí)現(xiàn)上述實(shí)施例方法中的全部或部分流程，可以通過計(jì)算機(jī)程序來指令相關(guān)的硬件來完成，所述的程序可存儲(chǔ)于一般計(jì)算機(jī)可讀取存儲(chǔ)介質(zhì)中，該程序在執(zhí)行時(shí)，可包括如上述各方法的實(shí)施例的流程。其中，所述的存儲(chǔ)介質(zhì)可為磁碟、光盤、只讀存儲(chǔ)記憶體(Read-Only Memory, ROM)或隨機(jī)存儲(chǔ)記憶體(Random AccessMemory, RAM)等。
[0164]本領(lǐng)域技術(shù)人員還可以了解到本發(fā)明實(shí)施例列出的各種功能是通過硬件還是軟件來實(shí)現(xiàn)取決于特定的應(yīng)用和整個(gè)系統(tǒng)的設(shè)計(jì)要求。本領(lǐng)域技術(shù)人員可以對(duì)于每種特定的應(yīng)用，可以使用各種方法實(shí)現(xiàn)所述的功能，但這種實(shí)現(xiàn)不應(yīng)被理解為超出本發(fā)明實(shí)施例保護(hù)的范圍。
[0165]本發(fā)明中應(yīng)用了具體實(shí)施例對(duì)本發(fā)明的原理及實(shí)施方式進(jìn)行了闡述，以上實(shí)施例的說明只是用于幫助理解本發(fā)明的方法及其核心思想；同時(shí)，對(duì)于本領(lǐng)域的一般技術(shù)人員，依據(jù)本發(fā)明的思想，在【具體實(shí)施方式】及應(yīng)用范圍上均會(huì)有改變之處，綜上所述，本說明書內(nèi)容不應(yīng)理解為對(duì)本發(fā)明的限制。
【權(quán)利要求】
1.一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的方法，其特征是，所述的方法具體包括: 獲取當(dāng)前砂體的測(cè)井?dāng)?shù)據(jù)； 對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)； 對(duì)所述的廢棄河道數(shù)據(jù)進(jìn)行擬合，得到平面分布的深泓線； 根據(jù)所述的深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體； 根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體； 根據(jù)所述廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型； 對(duì)所述的矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。
2.根據(jù)權(quán)利要求1所述的方法，其特征是，根據(jù)所述的深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體具體包括: 從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取砂體厚度； 將所述深泓線上的任一點(diǎn)作為數(shù)據(jù)點(diǎn)； 根據(jù)所述的數(shù)據(jù)點(diǎn)以及曲河流點(diǎn)壩砂體經(jīng)驗(yàn)公式確定側(cè)積層水平寬度以及河道寬 度； 根據(jù)所述的數(shù)據(jù)點(diǎn)、側(cè)積層水平寬度、河道寬度以及砂體厚度確定廢棄河道的剖面； 遍歷所述深泓線上的點(diǎn)，得到多個(gè)廢棄河道的剖面； 多個(gè)廢棄河道的剖面之間圈定的空間體積即為廢棄河道的三維幾何體。
3.根據(jù)權(quán)利要求2所述的方法，其特征是，根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體具體包括: 根據(jù)反曲曲線公式確定在實(shí)數(shù)范圍內(nèi)的反曲曲線形態(tài)； 對(duì)所述的反曲曲線形態(tài)進(jìn)行坐標(biāo)平移或坐標(biāo)旋轉(zhuǎn)，得到側(cè)積層剖面； 從所述的側(cè)積層數(shù)據(jù)中獲取側(cè)積層中心厚度、側(cè)積層水平跨度； 設(shè)置兩個(gè)圓心距為側(cè)積層中心厚度、直徑為側(cè)積層水平跨度的圓； 將所述側(cè)積層剖面對(duì)應(yīng)的兩個(gè)反曲曲線沿兩個(gè)圓相交的圓弧線進(jìn)行平行移動(dòng)； 平行移動(dòng)側(cè)積層剖面所圈定的空間即為側(cè)積層的空間分布體。
4.根據(jù)權(quán)利要求3所述的方法，其特征是，根據(jù)所述廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型網(wǎng)格具體包括: 將所述的點(diǎn)壩砂體數(shù)據(jù)填充到所述廢棄河道三維幾何體所圈定的區(qū)域中，得到第一填充三維幾何體； 從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層在井身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度； 根據(jù)所述側(cè)積層在井身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度，將所述側(cè)積層的空間分布體填充到第一填充三維幾何體，得到第二填充幾何體； 從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層的水平密度； 根據(jù)所述的水平密度將將所述側(cè)積層的空間分布體填充到所述第二填充幾何體，得到矢量化的砂體儲(chǔ)層構(gòu)型模型。
5.根據(jù)權(quán)利要求4所述的方法，其特征是，對(duì)所述的矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型具體包括: 根據(jù)所述的側(cè)積層數(shù)據(jù)確定網(wǎng)格尺寸； 根據(jù)所述的網(wǎng)格尺寸對(duì)所述矢量化的砂體儲(chǔ)層構(gòu)型模型建立模擬區(qū)的網(wǎng)格； 對(duì)所述模擬區(qū)的網(wǎng)格根據(jù)構(gòu)型要素類型進(jìn)行數(shù)值填充，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。
6.根據(jù)權(quán)利要求5所述的方法，其特征是，對(duì)所述模擬區(qū)的網(wǎng)格根據(jù)構(gòu)型要素類型進(jìn)行數(shù)值填充具體包括: 判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為廢棄河道； 當(dāng)判斷為是時(shí)，將所述的網(wǎng)格填充為廢棄河道數(shù)據(jù)； 否則，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為側(cè)積層； 當(dāng)判斷為是時(shí)，將所述的網(wǎng)格填充為側(cè)積層數(shù)據(jù)； 否則，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為點(diǎn)壩砂體； 當(dāng)判斷為是時(shí) ，將所述的網(wǎng)格填充為點(diǎn)壩砂體數(shù)據(jù)； 否則，確定所述網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型的比例； 獲取所述比例中的最大值； 確定所述最大值對(duì)應(yīng)的構(gòu)型要素類型； 將所述的網(wǎng)格填充為最大值對(duì)應(yīng)的構(gòu)型要素類型對(duì)應(yīng)的數(shù)據(jù)。
7.一種基于空間矢量的砂體儲(chǔ)層構(gòu)型建模的系統(tǒng)，其特征是，所述的系統(tǒng)具體包括: 測(cè)井?dāng)?shù)據(jù)獲取裝置，用于獲取當(dāng)前砂體的測(cè)井?dāng)?shù)據(jù)； 測(cè)井?dāng)?shù)據(jù)處理裝置，用于對(duì)所述的測(cè)井?dāng)?shù)據(jù)進(jìn)行數(shù)字處理，得到側(cè)積層數(shù)據(jù)、廢棄河道數(shù)據(jù)以及點(diǎn)壩砂體數(shù)據(jù)； 擬合裝置，用于對(duì)所述的廢棄河道數(shù)據(jù)進(jìn)行擬合，得到平面分布的深泓線； 三維幾何體建立裝置，用于根據(jù)所述的深泓線以及點(diǎn)壩砂體數(shù)據(jù)建立廢棄河道三維幾何體； 空間分布體建立裝置，用于根據(jù)反曲曲線公式、側(cè)積層數(shù)據(jù)確定側(cè)積層的空間分布體； 矢量化模型確定裝置，用于根據(jù)所述廢棄河道三維幾何體、側(cè)積層的空間分布體以及點(diǎn)壩砂體數(shù)據(jù)確定矢量化的砂體儲(chǔ)層構(gòu)型模型； 網(wǎng)格化裝置，用于對(duì)所述的矢量化的砂體儲(chǔ)層構(gòu)型模型進(jìn)行網(wǎng)格化，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。
8.根據(jù)權(quán)利要求7所述的系統(tǒng)，其特征是，所述的三維幾何體建立裝置具體包括: 砂體厚度獲取模塊，用于從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取砂體厚度； 數(shù)據(jù)點(diǎn)選取模塊，用于將所述深泓線上的任一點(diǎn)作為數(shù)據(jù)點(diǎn)； 河道寬度確定模塊，用于根據(jù)所述的數(shù)據(jù)點(diǎn)以及曲河流點(diǎn)壩砂體經(jīng)驗(yàn)公式確定側(cè)積層水平寬度、河道寬度； 廢棄河道剖面確定模塊，用于根據(jù)所述的數(shù)據(jù)點(diǎn)、側(cè)積層水平寬度、河道寬度以及砂體厚度確定廢棄河道的剖面； 數(shù)據(jù)點(diǎn)遍歷模塊，用于遍歷所述深泓線上的點(diǎn)，得到多個(gè)廢棄河道的剖面； 三維幾何體確定模塊，多個(gè)廢棄河道的剖面之間圈定的空間體積即為廢棄河道的三維幾何體。
9.根據(jù)權(quán)利要求8所述的系統(tǒng)，其特征是，所述的空間分布體建立裝置具體包括: 反曲曲線形態(tài)確定模塊，用于根據(jù)反曲曲線公式確定在實(shí)數(shù)范圍內(nèi)的反曲曲線形態(tài)； 側(cè)積層剖面確定模塊，用于對(duì)所述的反曲曲線形態(tài)進(jìn)行坐標(biāo)平移或坐標(biāo)旋轉(zhuǎn)，得到側(cè)積層劑面； 水平跨度獲取模塊，用于從所述的側(cè)積層數(shù)據(jù)中獲取側(cè)積層中心厚度、側(cè)積層水平跨度； 設(shè)置模塊，用于設(shè)置兩個(gè)圓心距為側(cè)積層中心厚度、直徑為側(cè)積層水平跨度的圓；平行移動(dòng)模塊，用于將所述側(cè)積層剖面對(duì)應(yīng)的兩個(gè)反曲曲線沿兩個(gè)圓相交的圓弧線進(jìn)行平行移動(dòng)； 空間分布體確定模塊，平行移動(dòng)側(cè)積層剖面所圈定的空間即為側(cè)積層的空間分布體。
10.根據(jù)權(quán)利要求9所述的系統(tǒng)，其特征是，所述的矢量化模型確定裝置具體包括: 第一填充模塊，用于將所述的點(diǎn)壩砂體數(shù)據(jù)填充到所述廢棄河道三維幾何體所圈定的區(qū)域中，得到第一填充三維幾何體； 角度獲取模塊，用于從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層在井身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度； 第二填充模塊，用于根據(jù)所述側(cè)積層在井`身上的位置、側(cè)積層的發(fā)育角度以及側(cè)積層的水平密度，將所述側(cè)積層的空間分布體填充到第一填充三維幾何體，得到第二填充幾何體； 水平密度獲取模塊，用于從所述的點(diǎn)壩砂體數(shù)據(jù)中獲取側(cè)積層的水平密度； 矢量化模型確定模塊，用于根據(jù)所述的水平密度將將所述側(cè)積層的空間分布體填充到所述第二填充幾何體，得到矢量化的砂體儲(chǔ)層構(gòu)型模型。
11.根據(jù)權(quán)利要求10所述的系統(tǒng)，其特征是，所述的網(wǎng)格化裝置具體包括: 網(wǎng)格尺寸確定模塊，用于根據(jù)所述的側(cè)積層數(shù)據(jù)確定網(wǎng)格尺寸； 網(wǎng)格建立模塊，用于根據(jù)所述的網(wǎng)格尺寸對(duì)所述矢量化的砂體儲(chǔ)層構(gòu)型模型建立模擬區(qū)的網(wǎng)格； 數(shù)值填充模塊，用于對(duì)所述模擬區(qū)的網(wǎng)格根據(jù)構(gòu)型要素類型進(jìn)行數(shù)值填充，得到網(wǎng)格化的砂體儲(chǔ)層構(gòu)型模型。
12.根據(jù)權(quán)利要求11所述的系統(tǒng)，其特征是，所述的數(shù)值填充模塊具體包括: 第一判斷單元，用于判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為廢棄河道； 第一填充單元，用于當(dāng)所述的第一判斷單元判斷為是時(shí)，將所述的網(wǎng)格填充為廢棄河道數(shù)據(jù)； 第二判斷單元，用于當(dāng)所述的第一判斷單元判斷為否時(shí)，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為側(cè)積層； 第二填充單元，用于當(dāng)所述的第二判斷單元判斷為是時(shí)，將所述的網(wǎng)格填充為側(cè)積層數(shù)據(jù)； 第三判斷單元，用于當(dāng)所述的第二判斷單元判斷為否時(shí)，判斷所述模擬區(qū)的網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型是否僅為點(diǎn)壩砂體； 第三填充單元，用于當(dāng)所述的第三判斷單元判斷為是時(shí)，將所述的網(wǎng)格填充為點(diǎn)壩砂體數(shù)據(jù)； 比例確定單元，用于當(dāng)所述第三判斷單元判斷為否時(shí)，確定所述網(wǎng)格對(duì)應(yīng)的構(gòu)型要素類型的比例； 最大值獲取單元，用于獲取所述比例中的最大值； 要素類型確定單元，用于確定所述最大值對(duì)應(yīng)的構(gòu)型要素類型； 第四填充單元，用于將 所述的網(wǎng)格填充為最大值對(duì)應(yīng)的構(gòu)型要素類型對(duì)應(yīng)的數(shù)據(jù)。
【文檔編號(hào)】G06F17/50GK103699751SQ201310744975
【公開日】2014年4月2日 申請(qǐng)日期:2013年12月30日 優(yōu)先權(quán)日:2013年12月30日
【發(fā)明者】李宇鵬, 吳勝和 申請(qǐng)人:中國石油大學(xué)(北京)