考慮夾層影響的滲透率粗化方法
【專利摘要】本發(fā)明公開了一種考慮夾層影響的滲透率粗化方法，包括如下步驟：A、將精細(xì)地質(zhì)模型中的精細(xì)網(wǎng)格與粗化網(wǎng)格對應(yīng)；B、對粗化網(wǎng)格內(nèi)部的砂體進(jìn)行連通性判斷并對連通砂體進(jìn)行編號；C、針對每一個粗化網(wǎng)格，判斷是否有夾層對粗化網(wǎng)格不同方向造成遮擋；D、通過判斷不同粗化網(wǎng)格內(nèi)部夾層是否造成遮擋得到的結(jié)果對滲透率參數(shù)場進(jìn)行具體的粗化。本發(fā)明解決了地質(zhì)模型粗化過程中造成的夾層信息丟失，無法有效將精細(xì)地質(zhì)模型中的滲透率模型有效地過渡到粗化模型中的問題，可以廣泛應(yīng)用于三維地質(zhì)建模領(lǐng)域。
【專利說明】考慮夾層影響的滲透率粗化方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及三維地質(zhì)建模領(lǐng)域，特別是涉及一種考慮夾層影響的滲透率粗化方法。

【背景技術(shù)】
[0002]在油田開發(fā)中后期，夾層的研究越來越精細(xì)，其對油氣運移的影響越來越受到開發(fā)者們的重視，許多學(xué)者對夾層進(jìn)行了不同程度的研究。不同成因的夾層作為一種非滲透性的介質(zhì)，其導(dǎo)致的儲層非均質(zhì)性在不同程度上控制著剩余油的分布。
[0003]夾層引起的滲流屏障和滲流差異直接表現(xiàn)在不同方向的滲透率中，在建立地質(zhì)模型時需要對滲透率及夾層分布進(jìn)行精細(xì)研究才能對其進(jìn)行有效的表征。夾層的分布厚度變化較大，從幾十厘米到幾米、幾十米不等且分布極不規(guī)則。夾層造成的儲層非均質(zhì)性會對油氣運移及油田開發(fā)中后期的剩余油預(yù)測及開采造成較大程度的影響。不同形態(tài)的夾層對粗化網(wǎng)格的滲透率影響不一樣。圖1為四種不同形態(tài)的夾層:圖1(a)中的夾層阻礙了 I方向的滲流通道，粗化為一個網(wǎng)格后，該網(wǎng)格I方向應(yīng)該為非滲透，而J方向具有滲透性；圖1(b)中的傾斜夾層造成了對1、J方向的遮擋，粗化后1、J方向都應(yīng)為非滲透；圖1(c)與圖1(a) 一樣夾層僅造成了 I方向的遮擋；圖1(d)中夾層僅造成了 J方向的遮擋。
[0004]由于精細(xì)模型網(wǎng)格數(shù)量大，而受計算機(jī)硬件的限制，精細(xì)地質(zhì)模型和油藏數(shù)值模擬能夠承受的網(wǎng)格規(guī)模之間還有較大差異，為了合理解決模擬耗時與微機(jī)性能的問題，必須對精細(xì)模型進(jìn)行粗化。隨著微機(jī)性能的穩(wěn)步提升和建模方法的逐步優(yōu)化，儲層地質(zhì)模型越來越精細(xì)，對夾層的描述越來越清楚，將精細(xì)地質(zhì)模型中夾層的信息如何有效的保留到粗化模型中成為研究的難點。
[0005]單純對于滲透率的粗化，國外已經(jīng)有許多成熟的方法及新的研究，國內(nèi)的研究也取得了豐富成果，但對非滲透夾層對滲透率粗化結(jié)果的影響考慮甚少，特別是傾斜狀的夾層，如點壩側(cè)積層，其對油氣的側(cè)向遮擋直接控制了剩余油的分布。
[0006]滲透率不同于孔隙度、飽和度等標(biāo)量屬性，是一種張量屬性，粗化時要考慮到其方向性。滲透率的粗化方法有最簡單的單一算數(shù)平均法、幾何平均法、調(diào)和平均法及加權(quán)平均法，雖然速度快但得到結(jié)果是各向同性的，不能反映出不同方向上滲透性的差異。而利用達(dá)西定律及質(zhì)量守恒方程求解流動方程的粗化方法較多，能夠表現(xiàn)出滲透率的各向異性，但都沒有考慮非滲透夾層的影響:當(dāng)粗化網(wǎng)格的某一方向存在夾層遮擋時，該網(wǎng)格在對應(yīng)方向應(yīng)該為非滲透性的。
[0007]具體參見圖2和圖3，以圖2中1、J、K方向各5個網(wǎng)格粗化為一個網(wǎng)格為例，精細(xì)模型中存在傾斜狀遮擋的夾層，通過序貫高斯模擬方法計算得到精細(xì)網(wǎng)格的滲透率模擬結(jié)果,滲透率分布在50-100毫達(dá)西之間。
[0008]通過成熟商業(yè)化軟件包含的常用滲透率粗化方法對該精細(xì)網(wǎng)格進(jìn)行粗化，圖2中分別是算術(shù)平均法、幾何平均法、Direct1nal averaging方法及Flow-based upscaling方法粗化之后的結(jié)果。由于條件數(shù)據(jù)中存在傾斜夾層的遮擋，如圖2 (b)，精細(xì)網(wǎng)格粗化為一個網(wǎng)格后，該網(wǎng)格在K方向和J方向應(yīng)該為非滲透，而I方向具有滲透性。圖3 (a)為算術(shù)平均法粗化結(jié)果，粗化后結(jié)果為82毫達(dá)西，在50-100毫達(dá)西的正常范圍內(nèi)且不能體現(xiàn)不同方向的滲透性，顯然與實際不符；圖3(b)為幾何平均法粗化結(jié)果，顯示為O毫達(dá)西，結(jié)果與實際能夠?qū)?yīng)，但是得到的模擬結(jié)果為均質(zhì)性的，不同方向的滲透率都為O毫達(dá)西，實際上在I方向應(yīng)該有滲透性，該方法也無法體現(xiàn)滲透率粗化后的非均質(zhì)性；圖3((:)為Direct1nalaveraging方法粗化結(jié)果,該方法雖然能夠產(chǎn)生不同方向的滲透率，但是在這種情況下，粗化結(jié)果在1、J、K三個方向均為82毫達(dá)西，與實際分析結(jié)論相差甚遠(yuǎn)，不能反應(yīng)夾層對滲透率粗化的影響；其中，圖3(d)、圖3(e)和圖3(f)為Flow-based upscaling方法粗化結(jié)果，該方法粗化后1、J、K方向結(jié)果分別為82毫達(dá)西、42毫達(dá)西、34毫達(dá)西，滲透率在J、K方向降低的比較厲害，一定程度上體現(xiàn)了夾層的影響，但是該網(wǎng)格K方向仍然具有滲透性，說明該方法在表征夾層影響上仍然存在缺陷，不能有效刻畫夾層對滲透率模型粗化的影響。


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

[0009]本發(fā)明的目的是為了克服上述【背景技術(shù)】的不足，提供一種考慮夾層影響的滲透率粗化方法，解決了地質(zhì)模型粗化過程中造成的夾層信息丟失，無法有效將精細(xì)地質(zhì)模型中的滲透率模型有效地過渡到粗化模型中的問題。
[0010]本發(fā)明提供的一種考慮夾層影響的滲透率粗化方法，包括如下步驟:A、將精細(xì)地質(zhì)模型中的精細(xì)網(wǎng)格與粗化網(wǎng)格對應(yīng)；B、對粗化網(wǎng)格內(nèi)部的砂體進(jìn)行連通性判斷并對連通砂體進(jìn)行編號；C、針對每一個粗化網(wǎng)格，判斷是否有夾層對粗化網(wǎng)格不同方向造成遮擋；D、通過判斷不同粗化網(wǎng)格內(nèi)部夾層是否造成遮擋得到的結(jié)果對滲透率參數(shù)場進(jìn)行具體的粗化。
[0011]在上述技術(shù)方案中，所述步驟B中，連通性判斷并對連通砂體進(jìn)行編號的過程如下:1)從I方向開始對連通的砂體進(jìn)行編號:從i = O, j = O, k = O開始，以I, J, K依次增加的順序首先找到第一個為砂體的網(wǎng)格，編號為1，之后按照順序依次選擇一個為砂體的網(wǎng)格a，如果在I方向上與其相鄰的網(wǎng)格b為砂體，則將網(wǎng)格b的砂體編號賦值給砂體a ;否則判斷J方向上與其相鄰的網(wǎng)格c是否為砂體，如果是砂體，則將網(wǎng)格c的砂體編號賦值給砂體a ;如果1、J方向都沒有砂體，則按照相同的方法判斷K方向；如果1、J、K三個方向都沒有滿足條件的砂體，則對網(wǎng)格a進(jìn)行新的砂體編號；2)從J方向開始對連通的砂體進(jìn)行編號:判斷J方向兩個相鄰的網(wǎng)格m、η砂體編號是否一致，如果不一致，將編號大的所有網(wǎng)格的砂體編號賦值為該相鄰網(wǎng)格的編號；3)從K方向開始對連通的砂體進(jìn)行編號:判斷K方向兩個相鄰的網(wǎng)格砂體編號是否一致，如果不一致，將編號大的所有網(wǎng)格的砂體編號賦值為該相鄰網(wǎng)格的編號。
[0012]在上述技術(shù)方案中，所述步驟B第I)項中，所述網(wǎng)格b在I方向上的編號為網(wǎng)格a在I方向上的編號減1，所述網(wǎng)格c在J方向上的編號為網(wǎng)格a在J方向上的網(wǎng)格編號減
1
[0013]在上述技術(shù)方案中，所述步驟C中，判斷粗化網(wǎng)格內(nèi)部某一砂體是否穿過了某一方向所在的精細(xì)網(wǎng)格的最大值和最小值分別對應(yīng)的兩個截面，以此來判斷該砂體是否穿過了整個粗化網(wǎng)格。
[0014]在上述技術(shù)方案中，所述步驟C中，同一夾層必須判斷1、J、K三個方向造成的遮擋情況。
[0015]在上述技術(shù)方案中，所述步驟D中，從1、J、K三個方向?qū)Υ只W(wǎng)格的滲透率分別粗化，某一方向夾層對該粗化網(wǎng)格造成遮擋，則該方向的粗化后的滲透率值直接賦值為O ;否則通過滲透率粗化方法對該網(wǎng)格該方向進(jìn)行滲透率粗化。
[0016]在上述技術(shù)方案中，還包括步驟E，所述步驟E的過程如下:單個粗化網(wǎng)格是否被夾層遮擋的判斷結(jié)束后，判斷一個泥巖夾層覆蓋幾個粗化網(wǎng)格，遍歷所有的砂體，記錄砂體所穿過的1、J、K三個方向的粗化網(wǎng)格，并記錄砂體的起始網(wǎng)格，找到所述砂體在某一方向上穿過的最小值所對應(yīng)的粗化網(wǎng)格X及最大值所對應(yīng)的粗化網(wǎng)格Y，如果砂體的起點不在網(wǎng)格X上，結(jié)束點不在網(wǎng)格Y上，找到具有相同性質(zhì)的網(wǎng)格X和網(wǎng)格Y，保證該砂體穿過網(wǎng)格X的最小截面和網(wǎng)格Y的最大截面，記錄網(wǎng)格X和網(wǎng)格Y在該方向上的最大值序號和最小值序號，除去粗化網(wǎng)格中序號在該方向上的小于最小值序號和大于最大值序號的網(wǎng)格，將所述砂體穿過的所有粗化網(wǎng)格的滲透率重新計算，以同樣的判斷方式對另外兩個方向進(jìn)行判斷，以達(dá)到判斷連續(xù)夾層對網(wǎng)格的遮擋情況。
[0017]在上述技術(shù)方案中，還包括步驟F，輸出不同方向的滲透率粗化結(jié)果。
[0018]在上述技術(shù)方案中，所述步驟A中，將精細(xì)地質(zhì)模型中的精細(xì)網(wǎng)格與粗化網(wǎng)格對應(yīng)之前，建立精細(xì)地質(zhì)模型，包括夾層及滲透率模型；并建立相應(yīng)的粗化骨架模型。
[0019]本發(fā)明考慮夾層影響的滲透率粗化方法，具有以下有益效果:通過本發(fā)明得到的滲透率粗化結(jié)果既能夠有效刻畫夾層對粗化后滲透率的影響，也能夠保證粗化后滲透率的非均值性。
[0020]運用本發(fā)明能夠有效利用地質(zhì)上認(rèn)識的夾層，并將夾層信息保留到粗化后的滲透率參數(shù)場中，更加準(zhǔn)確地刻畫流體的滲流路徑，為剩余油的挖潛提供決策依據(jù)。

【專利附圖】

【附圖說明】
[0021]圖1為不同形態(tài)泥巖夾層對滲透率粗化結(jié)果的影響示意圖；
[0022]圖2為示例中通過現(xiàn)有的滲透率粗化方法得到的滲透率模擬結(jié)果及夾層分布情況；
[0023]圖3為常用的滲透率粗化方法粗化后的滲透率模擬結(jié)果；
[0024]圖4為本發(fā)明考慮夾層影響的滲透率粗化方法的流程示意圖；
[0025]圖5為本發(fā)明考慮夾層影響的滲透率粗化方法步驟B中連通砂體編號示意圖；
[0026]圖6為本發(fā)明考慮夾層影響的滲透率粗化方法步驟C中I方向砂體連通性判斷示意圖；
[0027]圖7為本發(fā)明考慮夾層影響的滲透率粗化方法步驟E中一個泥巖夾層覆蓋多個粗化網(wǎng)格時連續(xù)夾層對網(wǎng)格遮擋的判斷示意圖；
[0028]圖8為本發(fā)明考慮夾層影響的滲透率粗化方法中對示例的滲透率粗化結(jié)果；

【具體實施方式】
[0029]下面結(jié)合附圖及實施例對本發(fā)明作進(jìn)一步的詳細(xì)描述，但該實施例不應(yīng)理解為對本發(fā)明的限制。
[0030]圖1至圖3中現(xiàn)有滲透率粗化方法在【背景技術(shù)】中已有描述，在此不再贅述。
[0031]參見圖4，本發(fā)明考慮夾層影響的滲透率粗化方法，包括如下步驟:
[0032]A、建立精細(xì)地質(zhì)模型，包括夾層及滲透率模型；并建立相應(yīng)的粗化骨架模型；將精細(xì)地質(zhì)模型中的精細(xì)網(wǎng)格與粗化網(wǎng)格對應(yīng)。
[0033]B、參見圖5，對粗化網(wǎng)格內(nèi)部的砂體進(jìn)行連通性判斷并對連通砂體進(jìn)行編號，具體過程如下:
[0034]I)從I方向開始對連通的砂體進(jìn)行編號:從i = O, j = O, k = O開始，以I, J, K依次增加的順序首先找到第一個為砂體的網(wǎng)格，編號為1，之后按照順序依次選擇一個為砂體的網(wǎng)格a，如果在I方向上與其相鄰的網(wǎng)格b為砂體(所述網(wǎng)格b在I方向上的編號為網(wǎng)格a在I方向上的編號減I)，則將網(wǎng)格b的砂體編號賦值給砂體a ;否則判斷J方向上與其相鄰的網(wǎng)格c是否為砂體(所述網(wǎng)格c在J方向上的編號為網(wǎng)格a在J方向上的網(wǎng)格編號減I)，如果是砂體，則將網(wǎng)格c的砂體編號賦值給砂體a ;如果1、J方向都沒有砂體，則按照相同的方法判斷K方向；如果1、J、K三個方向都沒有滿足條件的砂體，則對網(wǎng)格a進(jìn)行新的砂體編號，以二維例子為例，得到的結(jié)果如圖5(a)；
[0035]2)從J方向開始對連通的砂體進(jìn)行編號:參見圖5(a)，判斷J方向兩個相鄰的網(wǎng)格m、η砂體編號是否一致,如果不一致，將編號大的所有網(wǎng)格的砂體編號賦值為該相鄰網(wǎng)格的編號，圖5(a)中將砂體編號2賦值為1，得到圖5(b)的結(jié)果，并最后得到圖5(c)的結(jié)果;
[0036]3)從K方向開始對連通的砂體進(jìn)行編號:判斷K方向兩個相鄰的網(wǎng)格砂體編號是否一致，如果不一致，將編號大的所有網(wǎng)格的砂體編號賦值為該相鄰網(wǎng)格的編號。
[0037]C、針對每一個粗化網(wǎng)格，判斷是否有夾層對粗化網(wǎng)格不同方向造成遮擋即判斷連通砂體是否穿過整個粗化網(wǎng)格，由于粗化網(wǎng)格內(nèi)部，同一夾層可能造成粗化網(wǎng)格某一方向的遮擋，也可能造成幾個方向的遮擋，所以同一夾層必須判斷1、J、K三個方向造成的遮擋情況，針對每一個粗化網(wǎng)格，以I方向為例，判斷是否有連通的砂體穿過整個粗化網(wǎng)格(如圖6)。如果有，則該網(wǎng)格在I方向為連通，如果所有砂體在該粗化網(wǎng)格內(nèi)部都不連通，則該粗化網(wǎng)格在I方向被夾層遮擋，具體過程如下:
[0038]I)針對粗化網(wǎng)格內(nèi)部的某一個砂體k，如果該砂體網(wǎng)格數(shù)小于I方向粗化網(wǎng)格包含的精細(xì)網(wǎng)格數(shù)量，說明該砂體規(guī)模過小，不足以穿過整個粗化網(wǎng)格，重新選擇該粗化網(wǎng)格內(nèi)部的另一砂體，否則進(jìn)入下一步；
[0039]2)判斷粗化網(wǎng)格內(nèi)部某一砂體是否穿過了某一方向所在的精細(xì)網(wǎng)格的最大值和最小值分別對應(yīng)的兩個截面，以此來判斷該砂體是否穿過了整個粗化網(wǎng)格，仍以I方向為例，識別出粗化網(wǎng)格中I方向上對應(yīng)的精細(xì)網(wǎng)格最小I值及最大I值所對應(yīng)的兩個截面(分別對應(yīng)圖6中的A面和B面)，如果砂體k在這兩個截面上同時存在，則表示該砂體規(guī)模足夠大且穿過了整個粗化網(wǎng)格，說明該粗化網(wǎng)格在I方向連通，只要I方向連通，其他砂體則不用判斷；如果砂體k在這兩個截面上不同時存在則認(rèn)為該粗化網(wǎng)格在I方向被夾層遮擋，由于另外的砂體可能在該粗化網(wǎng)格內(nèi)部連通，繼續(xù)判斷另一砂體，重新執(zhí)行上一步；
[0040]3)以同樣的方式對J、K方向進(jìn)行判斷；
[0041]依次對每個粗化網(wǎng)格內(nèi)部砂體的連通性進(jìn)行判斷，即能得到夾層對所有粗化網(wǎng)格的遮擋情況。
[0042]D、通過不同粗化網(wǎng)格內(nèi)部夾層是否造成遮擋判斷得到的結(jié)果對滲透率參數(shù)場進(jìn)行具體的粗化:從1、J、K三個方向?qū)Υ只W(wǎng)格的滲透率分別粗化，某一方向夾層對該粗化網(wǎng)格造成遮擋，則該方向的粗化后的滲透率值直接賦值為O ;否則通過滲透率粗化方法對該網(wǎng)格該方向進(jìn)行滲透率粗化，該滲透率粗化方法包括前述的算術(shù)平均法、幾何平均法、Direct1nal averaging 方法及 Flow-based upscaling 方法等。
[0043]E、單個粗化網(wǎng)格是否被夾層遮擋的判斷結(jié)束后，判斷一個泥巖夾層覆蓋幾個粗化網(wǎng)格，如圖7所示(粗線外框表示粗化網(wǎng)格)，圖中包含4個粗化網(wǎng)格，分別為X1、X2、YUY2，單獨判讀時，只有X2網(wǎng)格在I方向具有滲透性；由于圖中的砂體具有連續(xù)性，導(dǎo)致X2、Yl及Y2三個網(wǎng)格整體上會具有滲透性。具體操作如下:遍歷所有的砂體，記錄砂體所穿過的1、J、K三個方向的粗化網(wǎng)格，并記錄砂體的起始網(wǎng)格，以圖7的I方向為例，找到所述砂體在I方向上穿過的最小I值所對應(yīng)的粗化網(wǎng)格(圖7中為Xl及X2網(wǎng)格)及最大I值所對應(yīng)的粗化網(wǎng)格(圖7中為Yl及Y2網(wǎng)格)；找到砂體在I方向的起點及結(jié)束點對應(yīng)的網(wǎng)格，圖中為X2及Y1，保證該砂體穿過網(wǎng)格X2的最小I截面和網(wǎng)格Yl的最大I截面(執(zhí)行步驟C第(2)項)，將所述砂體在I方向上穿過的所有粗化網(wǎng)格I方向的滲透率利用滲透率粗化方法進(jìn)行重新計算。以同樣的判斷方式對另外兩個方向進(jìn)行判斷，以達(dá)到判斷連續(xù)夾層對網(wǎng)格的遮擋情況。
[0044]F、輸出不同方向的滲透率粗化結(jié)果。
[0045]利用本發(fā)明對圖7的模擬結(jié)果粗化后得到的結(jié)果如圖8所示，在K方向和J方向滲透率粗化后結(jié)果均為O毫達(dá)西，I方向粗化后結(jié)果為40毫達(dá)西。最終得到結(jié)果與實際分析情況相符，較好的體現(xiàn)了夾層對滲透率的影響。
[0046]顯然，本領(lǐng)域的技術(shù)人員可以對本發(fā)明進(jìn)行各種改動和變型而不脫離本發(fā)明的精神和范圍。這樣，倘若本發(fā)明的這些修改和變型屬于本發(fā)明權(quán)利要求及其等同技術(shù)的范圍之內(nèi)，則本發(fā)明也意圖包含這些改動和變型在內(nèi)。
[0047]本說明書中未作詳細(xì)描述的內(nèi)容屬于本領(lǐng)域?qū)I(yè)技術(shù)人員公知的現(xiàn)有技術(shù)。
【權(quán)利要求】
1.一種考慮夾層影響的滲透率粗化方法，其特征在于:包括如下步驟: A、將精細(xì)地質(zhì)模型中的精細(xì)網(wǎng)格與粗化網(wǎng)格對應(yīng)； B、對粗化網(wǎng)格內(nèi)部的砂體進(jìn)行連通性判斷并對連通砂體進(jìn)行編號； C、針對每一個粗化網(wǎng)格，判斷是否有夾層對粗化網(wǎng)格不同方向造成遮擋； D、通過判斷不同粗化網(wǎng)格內(nèi)部夾層是否造成遮擋得到的結(jié)果對滲透率參數(shù)場進(jìn)行具體的粗化。
2.根據(jù)權(quán)利要求1所述的考慮夾層影響的滲透率粗化方法，其特征在于:所述步驟B中，連通性判斷并對連通砂體進(jìn)行編號的過程如下: 1)從I方向開始對連通的砂體進(jìn)行編號:從i= O, j = O, k = O開始，以I，J，K依次增加的順序首先找到第一個為砂體的網(wǎng)格，編號為1，之后按照順序依次選擇一個為砂體的網(wǎng)格a，如果在I方向上與其相鄰的網(wǎng)格b為砂體，則將網(wǎng)格b的砂體編號賦值給砂體a ;否則判斷J方向上與其相鄰的網(wǎng)格c是否為砂體，如果是砂體，則將網(wǎng)格c的砂體編號賦值給砂體a ;如果1、J方向都沒有砂體，則按照相同的方法判斷K方向；如果1、J、K三個方向都沒有滿足條件的砂體，則對網(wǎng)格a進(jìn)行新的砂體編號； 2)從J方向開始對連通的砂體進(jìn)行編號:判斷J方向兩個相鄰的網(wǎng)格m、n砂體編號是否一致，如果不一致，將編號大的所有網(wǎng)格的砂體編號賦值為該相鄰網(wǎng)格的編號； 3)從K方向開始對連通的砂體進(jìn)行編號:判斷K方向兩個相鄰的網(wǎng)格砂體編號是否一致，如果不一致，將編號大的所有網(wǎng)格的砂體編號賦值為該相鄰網(wǎng)格的編號。
3.根據(jù)權(quán)利要求2所述的考慮夾層影響的滲透率粗化方法，其特征在于:所述步驟B第I)項中，所述網(wǎng)格b在I方向上的編號為網(wǎng)格a在I方向上的編號減1，所述網(wǎng)格c在J方向上的編號為網(wǎng)格a在J方向上的網(wǎng)格編號減I。
4.根據(jù)權(quán)利要求1至3中任一項所述的考慮夾層影響的滲透率粗化方法，其特征在于:所述步驟C中，判斷粗化網(wǎng)格內(nèi)部某一砂體是否穿過了某一方向所在的精細(xì)網(wǎng)格的最大值和最小值分別對應(yīng)的兩個截面，以此來判斷該砂體是否穿過了整個粗化網(wǎng)格。
5.根據(jù)權(quán)利要求4所述的考慮夾層影響的滲透率粗化方法，其特征在于:所述步驟C中，同一夾層必須判斷1、J、K三個方向造成的遮擋情況。
6.根據(jù)權(quán)利要求1至3中任一項所述的考慮夾層影響的滲透率粗化方法，其特征在于:所述步驟D中，從1、J、K三個方向?qū)Υ只W(wǎng)格的滲透率分別粗化，某一方向夾層對該粗化網(wǎng)格造成遮擋，則該方向的粗化后的滲透率值直接賦值為O ;否則通過滲透率粗化方法對該網(wǎng)格該方向進(jìn)行滲透率粗化。
7.根據(jù)權(quán)利要求1至3中任一項所述的考慮夾層影響的滲透率粗化方法，其特征在于:還包括步驟E，所述步驟E的過程如下:單個粗化網(wǎng)格是否被夾層遮擋的判斷結(jié)束后，判斷一個泥巖夾層覆蓋幾個粗化網(wǎng)格，遍歷所有的砂體，記錄砂體所穿過的1、J、K三個方向的粗化網(wǎng)格，并記錄砂體的起始網(wǎng)格，找到所述砂體在某一方向上穿過的最小值所對應(yīng)的粗化網(wǎng)格X及最大值所對應(yīng)的粗化網(wǎng)格Y，如果砂體的起點不在網(wǎng)格X上，結(jié)束點不在網(wǎng)格Y上，找到具有相同性質(zhì)的網(wǎng)格X和網(wǎng)格Y，保證該砂體穿過網(wǎng)格X的最小截面和網(wǎng)格Y的最大截面，記錄網(wǎng)格X和網(wǎng)格Y在該方向上的最大值序號和最小值序號，除去粗化網(wǎng)格中序號在該方向上的小于最小值序號和大于最大值序號的網(wǎng)格，將所述砂體穿過的所有粗化網(wǎng)格的滲透率重新計算，以同樣的判斷方式對另外兩個方向進(jìn)行判斷，以達(dá)到判斷連續(xù)夾層對網(wǎng)格的遮擋情況。
8.根據(jù)權(quán)利要求1至3中任一項所述的考慮夾層影響的滲透率粗化方法，其特征在于:還包括步驟F，輸出不同方向的滲透率粗化結(jié)果。
9.根據(jù)權(quán)利要求1至3中任一項所述的考慮夾層影響的滲透率粗化方法，其特征在于:所述步驟A中，將精細(xì)地質(zhì)模型中的精細(xì)網(wǎng)格與粗化網(wǎng)格對應(yīng)之前，建立精細(xì)地質(zhì)模型，包括夾層及滲透率模型；并建立相應(yīng)的粗化骨架模型。
【文檔編號】G06T17/05GK104331925SQ201410510718
【公開日】2015年2月4日 申請日期:2014年9月28日 優(yōu)先權(quán)日:2014年9月28日
【發(fā)明者】張昌民, 李少華, 李君 , 喻思羽 申請人:長江大學(xué)