[發(fā)明專利]基于無監(jiān)督分類的二維地震數(shù)據(jù)全層位自動追蹤方法

有效

• 申請?zhí)枺?/td>: 201210470551.5
• 申請日: 2012-11-20
• 公開（公告）號: CN102981182A
• 公開（公告）日: 2013-03-20
• 發(fā)明（設(shè)計）人: 陳小二;鄒文;陶正喜;巫盛洪;周晶晶;杜洪;劉璞;巫駿;呂文彪
• 申請（專利權(quán)）人: 中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京銘碩知識產(chǎn)權(quán)代理有限公司 11286
• 代理人: 譚昌馳

1.一種基于無監(jiān)督分類的二維地震數(shù)據(jù)全層位自動追蹤方法，所述方法包括：

讀入二維地震數(shù)據(jù)以及控制層位數(shù)據(jù)；

查找二維地震數(shù)據(jù)波形的極值點；

基于查找的極值點，通過切比雪夫多項式來擬合地震波形，將擬合系數(shù)作為特征向量；

對擬合的地震波形進(jìn)行基于特征向量的無監(jiān)督分類，從而得到二維地震數(shù)據(jù)的全層位自動追蹤結(jié)果。

2.如權(quán)利要求1所述的方法，其中，所述極值點包括極大值或極小值或者查找二維地震數(shù)據(jù)波形的過零點并基于查找的過零點通過切比雪夫多項式來擬合地震波形。

3.如權(quán)利要求1所述的方法，其中，對擬合的地震波形進(jìn)行基于特征向量的無監(jiān)督分類的步驟包括：

根據(jù)特征向量的概率密度函數(shù)建立統(tǒng)計模型：

f(a|c,u,Σ)=Πk∈ϵf(ak|ck,uck,Σck)]]>

=Πk∈ϵ(2π)-nA/2|Σck|-1/2exp{-(ak-uck)′Σ-1ck(ak-uck)/2}---(3),]]>

其中，是特征向量A_k的邊緣分布，c_k∈{1,2,…,n_c}表示分類的類標(biāo)識，總共分為n_c類，表示高斯分布的均值，表示相應(yīng)的方差；

采用最大似然方法對特征向量進(jìn)行聚類分析：

ML=maxθ(f(A|θ))=maxθ(Πk∈ϵf(ak|ck,uck,Σck))---(4),]]>

通過對公式（4）進(jìn)行求導(dǎo)得到：

uc^=1|ϵc|Σk∈ϵcAk---(5),]]>

Σc^=1|ϵc|Σk∈ϵc(Ak-uc^)2---(6),]]>

ck^=argmaxckf(ak|ck,uc,Σc)---(7),]]>

其中，和分別為u，∑和c的估計值；

通過EM最大期望算法不斷迭代得到相應(yīng)的估計值：

uc^=1|ϵ|Σk∈ϵf(c|Ak,uc^,Σc^)Ak---(8)]]>

Σc^=1|ϵ|Σk∈ϵf(c|Ak,uc^,Σc^)(Ak-uc^)2---(9)]]>

通過公式（8）和（9）得到各個聚類團(tuán)的均值和協(xié)方差矩陣；

將各個分類樣本點映射到各自對應(yīng)的聚類團(tuán)中，即將各個極值點分配到不同的類型中，這些類型對應(yīng)于初始得到的各個層位片段所經(jīng)過的層位點，從而得到二維地震數(shù)據(jù)的全層位自動追蹤結(jié)果。