[發明專利]面向地質鉆孔的虛擬地層推斷方法及裝置

在審

• 申請號: 202110180713.0
• 申請日: 2021-02-08
• 公開（公告）號: CN112837413A
• 公開（公告）日: 2021-05-25
• 發明（設計）人: 沈言根;李安波;黃鍵初;閭國年
• 申請（專利權）人: 南京師范大學
• 主分類號: G06T17/05
• 分類號: G06T17/05;G06T17/20;G06N5/04
• 代理公司: 南京蘇高專利商標事務所(普通合伙) 32204
• 代理人: 馮艷芬

1.一種面向地質鉆孔的虛擬地層推斷方法，其特征在于該方法包括：

(1)讀取鉆孔數據、鉆孔地層信息和約束邊界矢量數據，生成鉆孔集合D、地層集合S、約束邊界集合Boundary和地層厚度集合H；

(2)基于約束邊界集合Boundary，將鉆孔集合D劃分為邊界鉆孔集合BD和非邊界鉆孔集合FBD，并根據地層集合S確定需推斷地層的地層編號區間C；

(3)從地層編號區間C中獲取任一地層編號c，構建推斷鉆孔集合ND與非推斷鉆孔集合ID；

(4)基于鉆孔集合ID，構建尖滅鉆孔集合AD，合并集合ID和集合AD為鉆孔集合UD，并根據集合UD創建三角網M；

(5)基于三角網M、邊界鉆孔集合BD和非邊界鉆孔集合FBD，獲取推斷鉆孔集合ND中所有鉆孔缺失的虛擬底層，并采用反距離加權插值計算虛擬地層的厚度；

(6)循環執行步驟(3)-(5)，直至地層編號區間C被遍歷完，完成所有鉆孔虛擬地層的推斷，更新地層信息，保存含有虛擬地層的鉆孔數據。

2.根據權利要求1所述的面向地質鉆孔的虛擬地層推斷方法，其特征在于：步驟(1)具體包括：

(1-1)讀取鉆孔數據到鉆孔集合D＝{d_i|i＝1，2，…,DN}；其中，i表示鉆孔序號，DN表示鉆孔數量；

(1-2)讀取所有鉆孔的地層信息，存入地層集合S＝{s_ij|i＝1,2,…,DN,j＝1,2,…,SN_i}，其中，s_ij表示鉆孔d_i的第j個地層，SN_i表示鉆孔d_i的地層數量；

(1-3)讀取約束邊界矢量數據到約束邊界集合Boundary；

(1-4)計算每一地層的厚度，存入地層厚度集合H＝{s_ij.h|i＝1,2,…,DN,j＝1,2,…,SN_i},其中，s_ij.h＝s_ij.et-s_ij.eb表示s_ij的地層厚度，s_ij.et為s_ij的層頂深度，s_ij.eb為s_ij的層底深度。

3.根據權利要求1所述的面向地質鉆孔的虛擬地層推斷方法，其特征在于：步驟(2)具體包括：

(2-1)獲取鉆孔集合D內每一鉆孔的坐標，判斷鉆孔是否位于約束邊界集合Boundary中的約束邊界上，若是則標記鉆孔為邊界鉆孔，并存入邊界鉆孔集合BD中，否則標記為非邊界鉆孔，并存入非邊界鉆孔集合FBD中；

(2-2)從地層集合S中獲取鉆孔集合D內所有鉆孔的末尾地層的地層編號，得到地層編號的最小值α與最大值β，構建需推斷地層的地層編號區間C＝[α+1，β],區間內取正整數。

4.根據權利要求1所述的面向地質鉆孔的虛擬地層推斷方法，其特征在于：步驟(3)具體包括：

(3-1)從地層編號區間C中獲取任一地層編號c；

(3-2)從地層集合S中獲取鉆孔集合D內每一鉆孔的末尾地層的地層編號，若該地層編號小于c，則代表該地層編號對應的鉆孔為需推斷鉆孔，將其保存到需推斷鉆孔集合ND，否則代表該地層編號對應的鉆孔為非推斷鉆孔，將其保存到非推斷鉆孔集合ID中；

(3-3)若非推斷鉆孔集合ID中的非推斷鉆孔數量LN≥3，則返回執行步驟(3-1)，否則執行步驟(4)。

5.根據權利要求1所述的面向地質鉆孔的虛擬地層推斷方法，其特征在于：步驟(4)具體包括：

(4-1)根據非推斷鉆孔集合ID，采用1/2尖滅規則構建尖滅鉆孔集合AD；

(4-2)合并非推斷鉆孔集合ID與尖滅鉆孔集合AD為鉆孔集合UD；

(4-3)根據鉆孔集合UD內的鉆孔坐標信息，創建Delaunary三角網M＝{m_t|t＝1,2,…,TM}，其中，m_t表示三角網內第t個三角形，TM表示三角網內三角形數量。