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

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

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


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

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

【專利附圖】

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

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