[發(fā)明專(zhuān)利]地層曲面構(gòu)建方法和裝置、設(shè)備及儲(chǔ)存介質(zhì)

在審

• 申請(qǐng)?zhí)枺?/td>: 201910745004.5
• 申請(qǐng)日: 2019-08-13
• 公開(kāi)（公告）號(hào): CN110610539A
• 公開(kāi)（公告）日: 2019-12-24
• 發(fā)明（設(shè)計(jì)）人: 陳敬柱;尹波;史江臣;王天明;梁建平;張守林
• 申請(qǐng)（專(zhuān)利權(quán)）人: 北京中色地科測(cè)繪有限公司
• 主分類(lèi)號(hào): G06T17/05
• 分類(lèi)號(hào): G06T17/05;G06T17/20
• 代理公司: 11751 北京市鼎立東審知識(shí)產(chǎn)權(quán)代理有限公司
• 代理人: 陳佳妹;賈滿(mǎn)意

1.一種地層曲面構(gòu)建方法，其特征在于，包括：

獲取待構(gòu)建的地層曲面的鉆孔數(shù)據(jù)，并基于獲取到的所述鉆孔數(shù)據(jù)構(gòu)建不規(guī)則三角網(wǎng)；其中，所述鉆孔數(shù)據(jù)的個(gè)數(shù)為多個(gè)；

基于所述不規(guī)則三角網(wǎng)，對(duì)多個(gè)所述鉆孔數(shù)據(jù)進(jìn)行綜合預(yù)處理，得到預(yù)處理結(jié)果；其中，所述預(yù)處理結(jié)果包括劃分為n級(jí)的鉆孔數(shù)據(jù)；

對(duì)所述預(yù)處理結(jié)果中的鉆孔數(shù)據(jù)進(jìn)行抽稀處理，得到抽稀處理后的鉆孔數(shù)據(jù)；其中，所述抽稀處理后的鉆孔數(shù)據(jù)參與構(gòu)建所述地層曲面；

利用曲面樣條函數(shù)對(duì)所述抽稀處理后的鉆孔數(shù)據(jù)進(jìn)行擬合，構(gòu)建得到所述地層曲面。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述基于獲取到的所述鉆孔數(shù)據(jù)構(gòu)建不規(guī)則三角網(wǎng)，包括：

將獲取到的所述鉆孔數(shù)據(jù)的坐標(biāo)作為離散點(diǎn)，在多個(gè)所述離散點(diǎn)中選取任意一點(diǎn)作為第一點(diǎn)，并由多個(gè)所述離散點(diǎn)中查找距所述第一點(diǎn)最近的離散點(diǎn)作為第二點(diǎn)，將所述第一點(diǎn)與所述第二點(diǎn)之間的連線作為初始基線；

在所述初始基線右側(cè)運(yùn)用Delaunay法則，由多個(gè)所述離散點(diǎn)中選取第三點(diǎn)，并依次連接所述第一點(diǎn)、所述第二點(diǎn)和所述第三點(diǎn)，生成相應(yīng)的三角形；

由生成的所述三角形中，選取除所述初始基線之外的另外兩邊分別作為重新選定的初始基線，生成相應(yīng)的三角形，直至所重新選定的初始基線的右側(cè)無(wú)所述離散點(diǎn)為止；

其中，所生成的各三角形構(gòu)建成所述不規(guī)則三角網(wǎng)。

3.根據(jù)權(quán)利要求1所述的方法，其特征在于，基于所述不規(guī)則三角網(wǎng)，對(duì)多個(gè)所述鉆孔數(shù)據(jù)進(jìn)行綜合預(yù)處理，包括：

獲取所述不規(guī)則三角網(wǎng)中各三角形的三個(gè)邊長(zhǎng)；

在當(dāng)前獲取的三角形的三個(gè)邊長(zhǎng)均小于預(yù)設(shè)閾值時(shí)，將當(dāng)前獲取的三角形的三個(gè)頂點(diǎn)綜合為一個(gè)點(diǎn)，并將綜合后得到的點(diǎn)作為綜合后的鉆孔數(shù)據(jù)；

其中，所述預(yù)設(shè)閾值表征待構(gòu)建的所述地層曲面在每一級(jí)別下的分辨率；

所述地層曲面的級(jí)別通過(guò)谷歌瓦片地圖算法進(jìn)行切分得到，且所述地層曲面的級(jí)別為n級(jí)，n的取值范圍為：0＜n≤22。

4.根據(jù)權(quán)利要求3所述的方法，其特征在于，對(duì)所述預(yù)處理結(jié)果中的鉆孔數(shù)據(jù)進(jìn)行抽稀處理，包括：

利用谷歌地圖的應(yīng)用程序編程接口將劃分為n級(jí)后的所述鉆孔數(shù)據(jù)的經(jīng)緯度坐標(biāo)換算為谷歌像素索引；

利用所述谷歌像素索引，計(jì)算所述第m級(jí)下鉆孔數(shù)據(jù)的左上角坐標(biāo)(x1，y1)和右下角坐標(biāo)(x2，y2)；

根據(jù)所述左上角坐標(biāo)(x1，y1)和所述右下角坐標(biāo)(x2，y2)生成第n級(jí)下鉆孔數(shù)據(jù)的p行、q列的抽稀格網(wǎng)；

遍歷每個(gè)所述抽稀格網(wǎng)內(nèi)的鉆孔數(shù)據(jù)，在當(dāng)前獲取的抽稀格網(wǎng)內(nèi)存在多個(gè)所述鉆孔數(shù)據(jù)時(shí)，保留距當(dāng)前獲取的抽稀格網(wǎng)中心距離最近的鉆孔數(shù)據(jù)，得到抽稀處理后的鉆孔數(shù)據(jù)；

其中，m的取值范圍為：0≤m≤n-4。

5.根據(jù)權(quán)利要求1至4任一項(xiàng)所述的方法，其特征在于，利用曲面樣條函數(shù)對(duì)所述抽稀處理后的鉆孔數(shù)據(jù)進(jìn)行擬合，構(gòu)建得到所述地層曲面，包括：

讀取根據(jù)平面彈性理論推導(dǎo)出的位移公式：

其中，A_i表征所述地層曲面i點(diǎn)上設(shè)置的垂直位移；Z_A(r)表征在所述地層曲面i點(diǎn)上設(shè)置垂直位移A_i，并限制距i點(diǎn)半徑為R以外位移為0時(shí)，由A_i引起距i點(diǎn)為r處的位移；i為參與構(gòu)建所述地層曲面的鉆孔數(shù)據(jù)中的任一鉆孔數(shù)據(jù)；

其中，r＝ri，x，y為自變量；

利用所述位移公式，對(duì)所述抽稀處理后的鉆孔數(shù)據(jù)進(jìn)行疊加，得到所述曲面樣條函數(shù)：

其中，N為所述抽稀處理后的鉆孔數(shù)據(jù)的個(gè)數(shù)；

讀取所建立的實(shí)測(cè)高程方程：

其中，h_i為實(shí)測(cè)高程，

基于抽稀處理后的鉆孔數(shù)據(jù)，對(duì)所述實(shí)測(cè)高程方程進(jìn)行求解，得到所述曲面樣條函數(shù)中的系數(shù)；

利用確定系數(shù)后的曲線樣條函數(shù)，對(duì)抽稀處理后的鉆孔數(shù)據(jù)進(jìn)行擬合，構(gòu)建得到所述地層曲面。