[發(fā)明專利]基于中心擴(kuò)散加傾角屬性的三維地震層位自動(dòng)追蹤方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201210299046.9
• 申請(qǐng)日: 2012-08-22
• 公開（公告）號(hào): CN102819040A
• 公開（公告）日: 2012-12-12
• 發(fā)明（設(shè)計(jì)）人: 姚興苗;楊析儒;胡光岷
• 申請(qǐng)（專利權(quán)）人: 電子科技大學(xué)
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 成都行之專利代理事務(wù)所(普通合伙) 51220
• 代理人: 溫利平

1.一種基于中心擴(kuò)散加傾角屬性的三維地震層位自動(dòng)追蹤方法，其特征在于：包括如下步驟：

步驟一、調(diào)整種子點(diǎn)位置：

如果種子點(diǎn)位置不在局部振幅最大的位置則對(duì)種子點(diǎn)位置按如下方法進(jìn)行調(diào)整：在種子點(diǎn)的時(shí)窗范圍內(nèi)尋找離種子點(diǎn)最近的極大值點(diǎn)，如果找到極大值點(diǎn)，則將種子點(diǎn)移動(dòng)到極大值點(diǎn)的位置；如果未找到極大值點(diǎn)，則種子點(diǎn)位置不動(dòng)；

步驟二、數(shù)據(jù)初始化準(zhǔn)備：

首先，計(jì)算出整個(gè)三維地震數(shù)據(jù)體的梯度向量；其次，將各初始的種子點(diǎn)的位置放入隊(duì)列queue中存放，并用一個(gè)數(shù)組array來標(biāo)識(shí)每個(gè)種子點(diǎn)是否已經(jīng)進(jìn)行過相關(guān)追蹤；

步驟三、中心擴(kuò)散追蹤：

1）計(jì)算種子點(diǎn)的傾角屬性：

利用步驟二計(jì)算出的梯度向量，計(jì)算出數(shù)據(jù)體中每個(gè)種子點(diǎn)的局部梯度方差矩陣，再求出該矩陣的特征向量，將特征向量投影到追蹤平面上，得到特征向量在追蹤平面的投影角度；

2）確定映射點(diǎn)：根據(jù)第1）步確定的投影角度算出映射點(diǎn)的位置；

3）確定追蹤點(diǎn)：

在第2）步確定的映射點(diǎn)的時(shí)窗范圍內(nèi)尋找振幅極大值點(diǎn)，在每個(gè)振幅極大值點(diǎn)周圍取一段波形數(shù)據(jù)，和種子點(diǎn)的波形數(shù)據(jù)進(jìn)行相關(guān)系數(shù)計(jì)算，然后找出相關(guān)系數(shù)最大的映射點(diǎn)，并計(jì)算該映射點(diǎn)和種子點(diǎn)的振幅比；判斷振幅比是否大于設(shè)定的門限值，如果是，則放棄該映射點(diǎn)，并將該映射點(diǎn)所在的種子道標(biāo)記為已追蹤；如果否，則將該映射點(diǎn)確定為追蹤點(diǎn)。

2.根據(jù)權(quán)利要求1所述的基于中心擴(kuò)散加傾角屬性的三維地震層位自動(dòng)追蹤方法，其特征在于：所述用數(shù)組array來標(biāo)識(shí)每個(gè)種子點(diǎn)是否已經(jīng)進(jìn)行過相關(guān)追蹤的方法是：初始化為0，表示沒有進(jìn)行過相關(guān)追蹤；如果進(jìn)行過相關(guān)追蹤，則標(biāo)識(shí)為1，表示進(jìn)行過相關(guān)追蹤。

3.根據(jù)權(quán)利要求1所述的基于中心擴(kuò)散加傾角屬性的三維地震層位自動(dòng)追蹤方法，其特征在于：根據(jù)第1）步確定的投影角度計(jì)算映射點(diǎn)位置的方法為：

首先按如下公式計(jì)算映射點(diǎn)和種子點(diǎn)相距的位置點(diǎn)數(shù)d：

d=h*tan(α)

其中，h表示種子道和相鄰道點(diǎn)的位置距離，α為特征向量在追蹤平面的投影角度；則，映射點(diǎn)的位置為a₀±d，其中，a₀為種子點(diǎn)的位置。

4.根據(jù)權(quán)利要求1所述的基于中心擴(kuò)散加傾角屬性的三維地震層位自動(dòng)追蹤方法，其特征在于：所述相關(guān)系數(shù)的計(jì)算公式為：

r=E([X-E(X)][Y-E(Y)])D(X)*D(Y)]]>

其中，X和Y分別是每個(gè)映射點(diǎn)的振幅值，E(X)和E(Y)分別是每段波形的平均振幅值，D(X)和D(Y)分別表示X和Y向量的方差。