[發(fā)明專利]一種地震數(shù)據(jù)非均勻間距等值線交互追蹤方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201310184553.2
• 申請(qǐng)日: 2013-05-17
• 公開(kāi)（公告）號(hào): CN104166158B
• 公開(kāi)（公告）日: 2016-11-09
• 發(fā)明（設(shè)計(jì)）人: 詹仕凡;丁建群;冉賢華;陶春峰
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油天然氣集團(tuán)公司;中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司 11127
• 代理人: 湯在彥

1.一種高效便捷的地震數(shù)據(jù)非均勻間距等值線交互追蹤方法，特點(diǎn)是采用以下步驟實(shí)現(xiàn)：

1）在工區(qū)采集任意一層地質(zhì)解釋層位數(shù)據(jù)；

2）采用逐點(diǎn)比較法求出層位數(shù)據(jù)的三維數(shù)據(jù)包圍盒，并以此作為數(shù)據(jù)的有效范圍；

3）在大地坐標(biāo)x,y平面上，采用帶約束條件的狄洛尼不規(guī)則三角網(wǎng)生成算法構(gòu)建層位數(shù)據(jù)不規(guī)則三角網(wǎng)；

4）采用逐點(diǎn)比較法統(tǒng)計(jì)層位數(shù)據(jù)不規(guī)則三角網(wǎng)中所有點(diǎn)的z值范圍，采用尾數(shù)截?cái)喾ㄗ詣?dòng)計(jì)算等值線起始值z(mì)_start和步長(zhǎng)z_step；

5）根據(jù)數(shù)據(jù)范圍確定等值線追蹤的起始值z(mì)_start和步長(zhǎng)z_step；

6）任意連續(xù)點(diǎn)選在層位數(shù)據(jù)生成并顯示在圖上的層位數(shù)據(jù)不規(guī)則三角網(wǎng)上，并根據(jù)二次曲線法依據(jù)這些點(diǎn)定義一條曲線L；

7）采用兩段直線求焦法求取曲線L與層位數(shù)據(jù)不規(guī)則三角網(wǎng)的交點(diǎn)系列SP_cross，并根據(jù)雙線性法計(jì)算交點(diǎn)系列SP_cross中每個(gè)交點(diǎn)p_cross的z值；

8）連接交點(diǎn)系列SP_cross，針對(duì)交點(diǎn)系列SP_cross中的任意兩兩相鄰點(diǎn)p_i、p_i+1采用線性平均插值方法獲得滿足定義的等值線點(diǎn)系列DP_cross；

9）以等值線點(diǎn)系列DP_cross為種子點(diǎn)系列，針對(duì)種子點(diǎn)系列中的每個(gè)種子點(diǎn)dp_cross進(jìn)行如下處理：

采用最小面積法計(jì)算出種子點(diǎn)dp_cross在層位數(shù)據(jù)不規(guī)則三角網(wǎng)中所在的三角形tri_cur，以三角形tri_cur為起點(diǎn)，采用優(yōu)化傳播法生成一條等值線L_contour；

10）在種子點(diǎn)dp_cross處進(jìn)行高程值標(biāo)注，完成追蹤一條等值線；將所述的等值線和追蹤完的所有等值線進(jìn)行合并，以避免重復(fù)的等值線；

11）重復(fù)步驟6）-10），追蹤生成一系列等間距等值線；

12）重復(fù)步驟5）-11），追蹤生成一系列非等間距等值線。

2.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟1）所述的層位數(shù)據(jù)由一系列點(diǎn)d(x,y,z)構(gòu)成的，其中x,y,z表示地質(zhì)解釋層位數(shù)據(jù)的空間三維坐標(biāo)；所述的x,y為大地坐標(biāo)，z為高程坐標(biāo)，對(duì)于任意兩個(gè)數(shù)據(jù)點(diǎn)，其x,y不完全相等。

3.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟2）所述的三維數(shù)據(jù)包圍盒是(x_min,x_max)、(y_min,y_max)、(z_min,z_max)，其中(x_min,x_max)分別是所有數(shù)據(jù)點(diǎn)x坐標(biāo)的最小值、最大值，(y_min,y_max)分別是所有數(shù)據(jù)點(diǎn)y坐標(biāo)的最小值、最大值，(z_min,z_max)分別是所有數(shù)據(jù)點(diǎn)z坐標(biāo)的最小值、最大值。

4.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟9）所述的優(yōu)化傳播法如下：

Step1將采用最小面積法計(jì)算所得到的三角形1tri_cur設(shè)為當(dāng)前三角形tri_cur；

采用下列公式求取等值點(diǎn)c(x,y,z):x=xa+(xb-xa)z-zazb-zay=ya+(yb-ya)z-zazb-zb]]>

式中：給定的高程值z(mì)，

三角形tri_cur的三條邊兩個(gè)頂點(diǎn)為a(x_a,y_a,z_a)和b(x_b,y_b,z_b)，

每個(gè)三角形在三邊上有2個(gè)等值點(diǎn)p₁、p₂；

將與三角形1tri_cur以點(diǎn)p₁所在的邊為相鄰邊的三角形2tri_ncur設(shè)為當(dāng)前三角形tri_cur，點(diǎn)p₁所在的邊為進(jìn)入邊；

Step2在三角形2tri_cur中搜索等值線的離去邊，即找到進(jìn)入相鄰三角形6tri_ecur的邊，假設(shè)該邊為e_next，找到后內(nèi)插該等值點(diǎn)的平面坐標(biāo)，內(nèi)插的方法和Step1中方法相同，只是當(dāng)前邊不參與計(jì)算；

Step3當(dāng)內(nèi)插等值點(diǎn)不是邊的兩個(gè)頂點(diǎn)之一，則找到e_next的相鄰三角形6tri_ecur，設(shè)該三角形6為tri_cur，當(dāng)前離去邊為進(jìn)入邊。

當(dāng)內(nèi)插等值點(diǎn)是邊的兩個(gè)頂點(diǎn)之一時(shí)，需要遍歷所有包含該等值點(diǎn)的三角形，逐個(gè)三角形進(jìn)行Step2操作，找到當(dāng)前點(diǎn)距離最短的等值點(diǎn)所在的離去邊，設(shè)該三角形11為tri_cur，當(dāng)前離去邊為進(jìn)入邊；

Step4重復(fù)Step2-Step3，直到離去邊為邊界邊，此時(shí)等值線為開(kāi)曲線L₁

Step5追蹤另一半等值線，設(shè)與Step1中的三角形1以點(diǎn)p₁所在的邊為相鄰邊的三角形3為當(dāng)前三角形tri_cur，以點(diǎn)p₂所在的邊為進(jìn)入邊，重復(fù)Step2-Step4操作，獲得另一段開(kāi)曲線L₂；

Step6合并開(kāi)曲線L₁和L₂，即得到一條等值線L_contour。