[發(fā)明專利]一種地質(zhì)體變屬性插值方法及系統(tǒng)

有效

• 申請?zhí)枺?/td>: 201811159460.3
• 申請日: 2018-09-30
• 公開（公告）號: CN110968930B
• 公開（公告）日: 2023-03-21
• 發(fā)明（設計）人: 孟祥賓
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06T19/00
• 代理公司: 北京思創(chuàng)畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種地質(zhì)體變屬性插值方法，其特征在于，包括：

步驟1：加載在設定工區(qū)范圍內(nèi)已解釋的地質(zhì)體成果數(shù)據(jù)；

步驟2：根據(jù)地質(zhì)體的控制點數(shù)據(jù)，構(gòu)建地質(zhì)體的表面網(wǎng)格模型；

步驟3：以塊體幾何表面三角形網(wǎng)格數(shù)據(jù)作為限定條件，對限定條件進行規(guī)則化處理，獲得規(guī)則的限定點、限定線段、限定面，獲得規(guī)則化表面網(wǎng)格模型；

步驟4：采用三維限定DELAUNAY剖分算法對所述規(guī)則化表面網(wǎng)格模型進行剖分，獲得具有拓撲一致性的四面體網(wǎng)格；

步驟5：利用基于梯度的變屬性插值算法對步驟4得到的四面體網(wǎng)格進行屬性插值，獲得插值結(jié)果；

利用幾何造型方法，構(gòu)建復雜地質(zhì)體的表面網(wǎng)格模型，幾何造型方法包括：二維限定Delaunay剖分、曲面求交、曲面縫合和復雜地質(zhì)體快速追蹤法；

步驟3具體包括：

步驟301：根據(jù)所述三角形網(wǎng)格數(shù)據(jù)，設置約束線段集和離散點集；

步驟302：將所述約束線段集包含的所有線段的端點記入所述離散點集；

步驟303：選擇所述約束線段集中的任意兩條線段，判斷所述兩條線段是否相交，若兩條線段相交且共線，則將共線的部分記入所述約束線段集，且在所述約束線段集中刪除所述兩條線段，若兩條線段僅相交不共線，則將交點記入所述離散點集；

步驟304：若三角形網(wǎng)格不滿足規(guī)則化的要求，則增加新點，判斷所述新點是否屬于所述離散點集，若不屬于所述離散點集，則將所述新點記入所述離散點集；判斷所述約束線段集中的線段是否包含所述新點，若包含所述新點，則以所述新點作為端點，將包含所述新點的線段劃分為兩條新線段，將所述兩條新線段記入所述約束線段集，且刪除包含所述新點的線段；

步驟305：重復執(zhí)行步驟303-304，直到所有三角形網(wǎng)格滿足規(guī)則化的要求為止，獲得規(guī)則的限定點、限定線段、限定面，獲得由規(guī)則的限定點、限定線段、限定面組成的規(guī)則化表面網(wǎng)格模型。

2.根據(jù)權(quán)利要求1所述的變屬性插值方法，其特征在于，根據(jù)三角形的面積和三角形的最小角判定是否滿足規(guī)則化的要求。

3.根據(jù)權(quán)利要求1所述的變屬性插值方法，其特征在于，步驟4具體包括：

步驟401：采用三維限定DELAUNAY剖分算法對所述規(guī)則化表面網(wǎng)格模型進行剖分獲得四面體網(wǎng)格；

步驟402：檢查限定條件三角形網(wǎng)格數(shù)據(jù)中每個點、每條邊、每個三角形是否在剖分結(jié)果中，若不在剖分結(jié)果中，則將交點分別加入到當前三角形網(wǎng)格和四面體網(wǎng)格；

步驟403：采用變網(wǎng)格方法對所述四面體網(wǎng)格進行優(yōu)化，獲得具有拓撲一致性的四面體網(wǎng)格。

4.根據(jù)權(quán)利要求3所述的變屬性插值方法，其特征在于，所述變網(wǎng)格方法為根據(jù)地質(zhì)體構(gòu)造形態(tài)增加輔助點來構(gòu)造四面體網(wǎng)格。

5.根據(jù)權(quán)利要求1所述的變屬性插值方法，其特征在于，采用基于梯度的變屬性插值算法對復雜地質(zhì)體進行屬性插值，所述變屬性插值算法使用的公式為：

SampleValue[i] ＝ mBaseValue +mGX*gx +mGY*gy + mGZ*gz (1)

其中，SampleValue[i]代表第i個四面體的屬性插值結(jié)果，mBaseValue代表待插四面體所在塊體的基礎屬性值，mGX、mGY、mGZ分別代表X、Y、Z坐標方向的屬性變化梯度，gx、gy、gz分別代表X、Y、Z坐標方向的坐標增量。