主权项 |
一种地震属性聚类方法,其特征在于,包括:输入预先聚类别个数k及待聚类的地震属性数据{x<sub>i</sub>,i=1,...,N}共N个;在N个地震属性数据中随机选取其中k个,作为k个类别中心的初始点{m<sub>p</sub>,p=1,...,k};为每一个地震属性数据到k个类别中心的距离d(x<sub>i</sub>,m<sub>p</sub>)估计下界l(x<sub>i</sub>,m<sub>p</sub>),初始化均为0;对每一个地震属性数据计算d(x<sub>i</sub>,m<sub>p</sub>),选取距离最近的类别中心代表的类别作为该地震属性数据的所属类别,并设置地震属性数据所属类别中心<img file="FDA0000970803810000011.GIF" wi="642" he="63" />对每一个地震属性数据设置到所属类别中心的距离估计值变量u(x<sub>i</sub>)=min<sub>p=1,...,k</sub>d(x<sub>i</sub>,m<sub>p</sub>);循环执行如下步骤直至类别中心收敛或循环迭代至设定次数,输出每个地震属性数据所属类别标号:对所有类别中心计算d(m<sub>q</sub>,m<sub>p</sub>),q=1,...,k,并设<img file="FDA0000970803810000012.GIF" wi="584" he="117" />确认所有地震属性数据u(x<sub>i</sub>)≥s(c(x<sub>i</sub>));对于所有同时满足条件m<sub>p</sub>≠class(x<sub>i</sub>)、u(x<sub>i</sub>)>l(x<sub>i</sub>,m<sub>p</sub>)及<img file="FDA0000970803810000013.GIF" wi="434" he="119" />的地震属性数据x<sub>i</sub>和类别中心m<sub>p</sub>:如果标志r(x<sub>i</sub>)为真,计算d(x<sub>i</sub>,c(x<sub>i</sub>))并置r(x<sub>i</sub>)为假,否则赋值d(x<sub>i</sub>,c(x<sub>i</sub>))=u(x<sub>i</sub>);如果d(x<sub>i</sub>,c(x<sub>i</sub>))>l(x<sub>i</sub>,m<sub>r</sub>),m<sub>r</sub>≠c(x<sub>i</sub>)或者<img file="FDA0000970803810000014.GIF" wi="573" he="119" />计算d(x<sub>i</sub>,m<sub>p</sub>);如果d(x<sub>i</sub>,m<sub>p</sub>)<d(x<sub>i</sub>,c(x<sub>i</sub>)),赋值c(x<sub>i</sub>)=m<sub>p</sub>;求取每个类别所属地震属性数据的平均值mean(m<sub>p</sub>);更新每个地震属性数据到所属类别中心的距离下界:l(x<sub>i</sub>,c(x<sub>i</sub>))=max{0,l(x<sub>i</sub>,c(x<sub>i</sub>)‑d(c(x<sub>i</sub>),mean(c(x<sub>i</sub>)))};更新每个地震属性数据相关的上下界估计值:u(x<sub>i</sub>)=u(x<sub>i</sub>)+d(c(x<sub>i</sub>),mean(c(x<sub>i</sub>))),置r(x<sub>i</sub>)为真;对k个类别中心进行更新m<sub>p</sub>=mean(m<sub>p</sub>),使用每个类别所占有的地震属性数据的均值作为新的类别中心;输入待聚类的地震属性数据之前,还包括对每一个地震属性数据进行高斯归一化处理;所述高斯归一化处理按如下公式进行:<maths num="0001"><math><![CDATA[<mrow><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>μ</mi></mrow><mi>σ</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000970803810000021.GIF" wi="245" he="102" /></maths>其中,{x<sub>i</sub>,i=1,2,...,N}为地震属性数据,{y<sub>i</sub>,i=1,2,...,N}为高斯归一化后的地震属性数据,<img file="FDA0000970803810000022.GIF" wi="766" he="128" /> |