076] A :構(gòu)建大范圍三角形，包含所有的集合ST中的所有散點(diǎn)，大三角形三個(gè)頂點(diǎn)不包 含在集合ST中；
[0077] B :從集合ST中任意選取一點(diǎn)，構(gòu)建第一個(gè)Delaunay三角剖分，共包含三個(gè)三角 形；
[0078] C :依次計(jì)算集合ST中尚未計(jì)算的點(diǎn)，找出使其位于相鄰兩三角形外接圓的公共 部分中的兩個(gè)相鄰三角形，并刪除相鄰三角形的公共邊，并連接該點(diǎn)與剩余兩點(diǎn)，形成新的 Delaunay三角剖分；
[0079] D :刪除含大范圍三角形頂點(diǎn)的三角形，即可得到最后的Delaunay三角剖分；
[0080] 同理依照以上方法，構(gòu)建地底三角形TRIb，以及每個(gè)大層的三角網(wǎng)剖分；
[0081] 步驟4、計(jì)算標(biāo)準(zhǔn)矩形區(qū)域G單元網(wǎng)格的計(jì)算線與網(wǎng)格線：
[0082] 計(jì)算點(diǎn)P1在三角網(wǎng)TRIT，TRIb上的投影點(diǎn)：由于點(diǎn)集P內(nèi)所有點(diǎn)均分布標(biāo)準(zhǔn)矩形 區(qū)域G內(nèi)，所以在三角網(wǎng)TRI T，TRIb上均能找到該點(diǎn)對(duì)應(yīng)的投影點(diǎn)，連接該點(diǎn)的兩投影點(diǎn)得 到線I1 ;如果點(diǎn)P1為網(wǎng)格計(jì)算點(diǎn)，則I1為計(jì)算線；如果點(diǎn)P1為網(wǎng)格節(jié)點(diǎn)，則I 1為輔助線；根 據(jù)三角形所確定的平面方程，計(jì)算出點(diǎn)P1在三角網(wǎng)TRIt、TRI b上的投影點(diǎn)的Z值分別為Zt、 Zb ;
[0083] 則當(dāng)前標(biāo)準(zhǔn)工程范圍表達(dá)成nxXny條計(jì)算線和（n x+l) X (ny+l)條輔助線，由網(wǎng)格 節(jié)點(diǎn)和網(wǎng)格計(jì)算點(diǎn)對(duì)應(yīng)關(guān)系可知每條計(jì)算線對(duì)應(yīng)四條輔助線；
[0084] 步驟5、等距劃分輔助線與計(jì)算線：
[0085] 以
為基準(zhǔn)點(diǎn)，以Uz為單位長(zhǎng)，平均劃分輔助線，稱劃分點(diǎn)為輔助點(diǎn)，則輔助 點(diǎn)Z值分布為Z = [n2, {n}，II1] Xuz，其中
' η為區(qū)間[η2， nj之間的自然數(shù)集合；
[0086] 以Z = 0為基準(zhǔn)點(diǎn)，以Uz為單位長(zhǎng)，平均劃分計(jì)算線，稱劃分點(diǎn)為計(jì)算點(diǎn)，則計(jì)算點(diǎn) Z值分布為：Z = [m2, Imhm1] XUz，其中
m為區(qū)間[1?,!?]之間的自然數(shù) 集合；
[0087] 如圖3所示，由等距劃分輔助線與計(jì)算線機(jī)制可知一個(gè)計(jì)算點(diǎn)對(duì)應(yīng)周圍八個(gè)輔助 點(diǎn)，對(duì)應(yīng)的輔助點(diǎn)構(gòu)建成該計(jì)算點(diǎn)所在的單元體；
[0088] 步驟6、計(jì)算離散點(diǎn)地質(zhì)屬性，確定每個(gè)離散點(diǎn)所在的單元體，并確認(rèn)單元體的地 質(zhì)屬性；
[0089] 把采集到的鉆孔空間位置信息、鉆孔的地層分布情況、地質(zhì)背景參數(shù)、地質(zhì)物探信 息、鉆孔綜合分層資料、離散點(diǎn)空間位置作為已知信息，通過(guò)模糊、概率、隨機(jī)過(guò)程的不確定 性方法，對(duì)離散點(diǎn)進(jìn)行度量計(jì)算，得到每種地質(zhì)屬性在該離散點(diǎn)的概率分布情況；
[0090] 離散點(diǎn)的地質(zhì)屬性計(jì)算過(guò)程為：
[0091] A :設(shè)鉆孔資料中共有m種地質(zhì)屬性，則其概率分布為集合L : {11，12, 13,…，lm}， 其中
[0092] B :取集合L中的最大值，記下標(biāo)為i，則Li>Lj，其中j = 1，2,…，i-1，i+1，...m ;則 第i種地質(zhì)屬性為該離散點(diǎn)的地質(zhì)屬性；
[0093] C :取集合L中的最大值，記下標(biāo)為il, i2,…，ik共k個(gè)，且Lil = Li2 =…= Lik>Lj,其中j古il, i2,…，ik ;則該離散點(diǎn)處于地質(zhì)體的分界處,記該點(diǎn)地質(zhì)屬性為第 il, i2,…，ik中的任意一種地質(zhì)屬性即可。
[0094] 步驟7、建立計(jì)算點(diǎn)與輔助點(diǎn)之間的對(duì)應(yīng)關(guān)系；
[0095] A :對(duì)計(jì)算線與輔助線進(jìn)行編號(hào)，編號(hào)規(guī)則：[列號(hào)，行號(hào)]，記為[I，J]
[0096] B :建立計(jì)算線與輔助線的對(duì)應(yīng)關(guān)系，編號(hào)為[I，J]的計(jì)算線對(duì)應(yīng)的輔助線的編號(hào) 為：[I，J]，[1+1，J]，[1+1，J+1]，[I，J+1];
[0097] C :依次對(duì)每條的輔助線上輔助點(diǎn)，從上至下以自然數(shù)順序進(jìn)行編號(hào)；依次對(duì)每條 的計(jì)算線上計(jì)算點(diǎn)，從上至下以自然數(shù)順序進(jìn)行編號(hào)；
[0098] D :從上至下，依次找出每個(gè)計(jì)算點(diǎn)對(duì)應(yīng)的輔助點(diǎn)；其對(duì)應(yīng)的原則為：輔助點(diǎn)Z值分 布位于[計(jì)算點(diǎn)Z值-Uz，計(jì)算點(diǎn)Z值+uj區(qū)間中；
[0099] 步驟8、構(gòu)建三維地質(zhì)柵格單元體：
[0100] A :-個(gè)計(jì)算點(diǎn)對(duì)應(yīng)周圍八個(gè)輔助點(diǎn)，對(duì)應(yīng)的輔助點(diǎn)構(gòu)建成該計(jì)算點(diǎn)所在的單元 體；
[0101] B :對(duì)該單元體的面進(jìn)行編號(hào)，統(tǒng)計(jì)該單元體的所有面的類型，將面分為堅(jiān)直面與 水平面，統(tǒng)計(jì)每個(gè)面包含的點(diǎn)編號(hào)，標(biāo)記面在網(wǎng)格劃分中的下標(biāo)，將所得統(tǒng)計(jì)信息即該單元 體合并計(jì)算信息；
[0102] C :記該計(jì)算點(diǎn)的地質(zhì)屬性為該計(jì)算點(diǎn)所在單元體的地質(zhì)屬性；
[0103] 如圖4所示，則將當(dāng)前工程范圍表達(dá)為由多個(gè)包含地質(zhì)屬性、合并計(jì)算信息、位置 信息的單元體集合的三維地質(zhì)柵格體T，其中單元體Tu k E T ;
[0104] 步驟9、先列后行，根據(jù)計(jì)算順序和單元體位置信息，將單元體分類；通過(guò)構(gòu)建單 元體，每一條計(jì)算線上計(jì)算點(diǎn)對(duì)應(yīng)的單元體形成一條單元體列，與計(jì)算線的行數(shù)與列數(shù)一 致，如圖5所示；
[0105] 以單元體T1, jik為例，其位置信息[i，j，k]，i為所在列順序數(shù)i e {1，2,……，nj， j為所在行順序數(shù)j e {1，2,……，ny}，K為堅(jiān)直方向順序數(shù)；
[0106] 類型1 :不與其他任何單元體相交；
[0107] 類型2 :與一個(gè)單元體相交；
[0108] 類型3 :與兩個(gè)單元體相交；
[0109] 類型4 :同時(shí)與三個(gè)單元體相交；
[0110] 步驟10、根據(jù)單元體的計(jì)算信息、位置信息將單元體合并；合并地質(zhì)屬性相同的 單元體，并得到每種地質(zhì)的分布實(shí)體：
[0111] 將不與其他任何單元體相交的單元體合并：定單元體TtJiO 1。為起始地質(zhì)體，且該地 質(zhì)體的地質(zhì)屬性與單元體τ。^。一致；
[0112] 實(shí)施例2
[0113] 本實(shí)施例的一種三維地質(zhì)柵格體模型構(gòu)建方法，與實(shí)施例1基本相同，不同之處 在于合并的單元體類型不一樣：
[0114] 將與一個(gè)單元體相交的單元體合并：如該單元體與相交單元體的地質(zhì)屬性不同， 不合并；反之，找出相交單元體的所在地質(zhì)體，從地質(zhì)體中刪除該單元體的相交表面以及其 所在地質(zhì)數(shù)據(jù)，并把其余5個(gè)面加入至相交單元體的所在地質(zhì)體中；該單元體與相交單元 體合并至同一個(gè)地質(zhì)體中。
[0115] 實(shí)施例3
[0116] 本實(shí)施例的一種三維地質(zhì)柵格體模型構(gòu)建方法，與實(shí)施例1基本相同，不同之處 在于合并的單元體類型不一樣：
[0117] 將與兩個(gè)單元體相交的單元體合并：如該單元體與兩相交單元體的地質(zhì)屬性均不 同，均不合并；如只與其中一相交單元體地質(zhì)屬性相同，則按照實(shí)施例2方式進(jìn)行合并即 可；如三單元體地質(zhì)屬性相同，把單元體加入至相鄰單元體所在地質(zhì)體中，且需合并相交單 元體所在的地質(zhì)體；
[0118] 實(shí)施例4
[0119] 本實(shí)施例的一種三維地質(zhì)柵格體模型構(gòu)建方法，與實(shí)施例1基本相同，不同之處 在于合并的單元體類型不一樣：
[0120] 將與三個(gè)單元體相交的單元體合并：如該單元體與兩相交單元體的地質(zhì)屬性均不 同，均不合并；如只與其中一相交單元體地質(zhì)屬性相同，則按照實(shí)施例3方式進(jìn)行合并即 可；如四單元體地質(zhì)屬性相同，合并四單元體成新的單元體。
[0121] 以上示意性地對(duì)本發(fā)明創(chuàng)造及其實(shí)施方式進(jìn)行了描述，實(shí)際的技術(shù)方案并不局限 于此。所以，如果本領(lǐng)域的普通技術(shù)人員受其啟示，在不脫離本創(chuàng)造宗旨的情況下，不經(jīng)創(chuàng) 造性的設(shè)計(jì)出與該技術(shù)方案相似的實(shí)施方式，均應(yīng)屬于本專利的保護(hù)范圍。
【主權(quán)項(xiàng)】
1. 一種Η維地質(zhì)柵格體模型構(gòu)建方法，方法流程如下： 步驟1、確定標(biāo)準(zhǔn)工程平面范圍：已知標(biāo)準(zhǔn)工程平面范圍為矩形區(qū)域[(XI，Χ2)，（yi， 72)]、柵格體的單位為；口枯，11^。，得到標(biāo)準(zhǔn)矩形區(qū)域6[山;。，義