角點網格是一種結構化網格類型,網格位置能用i,j,k定義,並且單元網格的長寬大小可變,垂向連接頂底網格點的網格面可以是傾斜的,能夠更加精確地描述斷層兩翼的深度變化流體分布和流體滲流特征。
假如角點網格在X,Y,Z方向上划分的數量分別為NI,NJ,NK.那么XY平面上的網格數量為NI*NJ。三維網格數量為NI*NJ*NZ。
ECLIPSE定義角點網格的由兩個關鍵字,即COORD和ZCORN.
(1)先用COORD定義線,一共有(NI+1)*(NJ+1)條線,其中N表示XY平面網格數目,一條線由兩個點來定,一個點有X、Y、Z三個數據,所以緊跟在COORD后面的數據一共有(NI+1)*(NJ+1)*2*3 個!注意:定義線的順序是一條線一條線的定義,方向是從左到右,再從左到右。
例如以3*3*3的網格系統,COORD就定義16條線。在eclipse中就有16行數據(見附錄文件)

(2)然后用ZCORN定義每個網格的點,一個網絡8個點嘛,所以緊跟在ZCORN后面的數據應該有N*8個!,這8個點全是Z值。注意:定義
網格角點的順序都是從左到右,再從左到右從頂到底!
如圖1中的3*3*3=27個網格,第一層網格有9個網格,則首先定義網格編號為(1,1,1),(2,1,1), (3,1,1)的上邊的Z值,然后定義網格編號為(1,1,1),(2,1,1), (3,1,1)的下邊的Z值;然后才是第二排網格和第三排網格。平面網格數量為3*3,每排(I)有六個數,共有6排。具體如下:1001,1001,1002,1002,1003,1003;//網格編號為(1,1,1)左上側,(1,1,1)右上側,(2,1,1)左上側,(2,1,1)右上側;(3,1,1)左上側,(3,1,1)
右上側,
1001,1001,1002,1002,1003,1003; //網格編號為(1,1,1)左下側,(1,1,1)右下側,(2,1,1)左下側,(2,1,1)右下側;(3,1,1)左下側,(3,1,1)
右下側,
1004 1004 1005 1005 1006 1006 //網格編號為(1,2,1)左上側,(1,2,1)右上側,(2,2,1)左上側,(2,2,1)右上側;(3,2,1)左上側,(3,2,1)
右上側,
1004 1004 1005 1005 1006 1006 //網格編號為(1,2,1)左下側,(1,2,1)右下側,(2,2,1)左下側,(2,2,1)右下側;(3,2,1)左下側,(3,2,1)
右下側,
1007 1007 1008 1008 1009 1009 //網格編號為(1,3,1)左上側,(1,3,1)右上側,(2,3,1)左上側,(2,3,1)右上側;(3,3,1)左上側,(3,3,1)
右上側,
1007 1007 1008 1008 1009 1009//網格編號為(1,3,1)左下側,(1,3,1)右下側,(2,3,1)左下側,(2,3,1)右下側;(3,3,1)左下側,(3,3,1)
右下側,
排完第一層網格的頂面,排完第一層網格的底面, 排完第2層網格的頂面,排完第2層網格的底面, 排完第3層網格的頂面,排完第3層網格的底面。
顯然,數據點為36*6=216.
附錄A:3*3*3網格系統,最初采用矩形網格來定義,具體的參數如下:



來源:http://www.wendangku.net/doc/33cbea8b680203d8ce2f242b-3.html
------------------------------------------------------------------------------------------------------------------------------------------------------------------------
來源:http://www.yanfengliu.com/2016/02/02/%e8%a7%92%e7%82%b9%e7%bd%91%e6%a0%bc%e8%af%a6%e8%a7%a3/
邏輯結構上屬於m×n×L個子單元的規則拓撲結構模型。其中,X方向邏輯上有(m+1)條線,剖分成m個單元;在Y方向上有(n+1)條線,剖分成n個單元;在深度Z方向上有(L+1)跳線,剖分成L個單元,每個單元屬於不規則六面體。各子單元室一個不規則六面體,定義其空間坐標的8個節點有頂、底兩個控制面積每個節點z值定義。
假如角點網格X,Y,Z方向上划分的數量分別為NI,NJ,NK,那么XY方向上的網格數量為NI×NJ。三維網格數量為NI×NJ×NK。
Eclipse定義角點網格有兩個關鍵詞,即COORD和ZCORN.
(1)先用COORD定義線,一共有(NI+1)×(NJ+1)條線,其中NI、NJ表示XY方向網格數目,一條線由兩個點來定,一個點有X、Y、Z三個數據,所以緊跟在COORD后面的數據一共有(NI+1)×(NJ+1)×2×3個。定義順序是一條一條先的定義,方向是從左到右,再從左到右。
(2)然后用ZCORD定義每個網格的點,一個網格8個點,所以緊跟在COORD后面的數據應該有NI×NJ×8個,這個8個點全部是Z值,定義網格角點的順序是從左到右,再從左到右,從頂到底。
假設一個網格的一條邊C、D在線AB上,C、D的z值為zc,zd,則C、D的空間坐標(xc,yc)、(xd,yd)由A(xa,ya,za)與B(xb,yb,zb)與兩個z值定義。
xc=xa+(zc-za)*(xb-xa)/(zb-za)
yc=ya+(zc-za)*(yb-ya)/(zb-za)
xd=xa+(zd-za)*(xb-xa)/(zb-dz)
yd=ya+(zd-za)*(yb-ya)/(zb-za)
即由頂部A點與底部B點線性內插所得。
3.2 角點網格讀取
首先讀取COORD的內容,這是(NI+1)×(NJ+1)條縱向上的線,每個線有2個端點組成,一個在頂部,一個在底部、控制了這條線上不同網格的xy的計算,把這些數據存儲到Line[NI+1, NJ+1, 2][3]的數組內。
然后讀取ZCORD的數據,它是每個網格的8個定點,最外層循環是z方向,次級是y方向,最內循環是x方向。
