多项式微商法圈定异常.pdfVIP

第 37卷 第 5期 成都理工大学学报 (自然科学版) Vo1．37No．5 2010年 1O月 JOURNALOFCHENGDUUNIVERSITYOFTECHNOLOGY (Science TechnologyEdition) Oct．2010 [文章编号]167卜9727(2010)05—0535—05 多 项 式 微 商 法 圈 定 异 常 蒋开明 顾 民 葛 良全 (成都理工大学应用核技术与 自动化工程学院，成都 610059) [摘要]在处理勘查数据过程中，重点是发现和确定数据异常区域，而异常区域的边界应该是 变化幅度最大的点构成的边界。利用二元多项式微商法，找到所有峰顶，对异常数据峰顶的陡 峭度和峰值幅度两方面进行综合判断，从而圈定异常峰，提取反映特性变化的峰的边缘，从而 圈定异常区域。文 中对峰值 的算法进行 了理论推导和分析 ，并给 出了异常圈定的处理步骤 。 实际处理数据的结果表明，该方法是有效和准确的。利用文 中方法圈定 出的异常面积 比采用 标准方差与等值线圈定异常方法小，表征 了真正的异常边界。 [关键词]多项式微商法；异常峰；异常边界；航空勘探 [分类号]P631 [文献标识码]A 在勘查数据异常处理过程 中，通过各种方法 1 异常圈定原理 得到剩余值后 ，合理确定剩余值异常下限至关重 要n]。目前最常用的异常下限值确定方法是平均 曲面上的峰可 以看成是二元函数极大值 ，因 值加相应倍数的标准差 ，然后在等值线 中显示 出 此可以通过求取二元 函数 的极大值 来寻找峰 来。这种方法的异常下限值仅从数值本身大小考 位 。 虑 ，同时，等值线显示出来 的是高于异常下限值的 要圈定异常的剩余值 ，用二元多项式对数据 区域 ，异常区域的边界不能代表真正的异常区域 进行拟合 ，然后求取极大值得到峰位 。具体做法 边界 ，异常区域的边界应该是变化幅度最大的点 为：首先用 (2m+1)×(2m+1)二维窗滑动取数 构成的边界 。从表面图观察异常，异常往往表现 据 ，对窗内数据进行二元多项式最IllS．乘拟合 ，计 为比较陡峭、峰值 比较大 的峰，因此，从峰高值和 算出窗 口中心位置的导数和偏导数，滑动窗 口，计 峰的陡峭程度两个方面综合判断是否是异常更为 算下一窗 口中心的导数和偏导数 。当得到所有数 合理。当异常峰 圈定后 ，通过提取 峰周 围变化幅 据的导数和偏导数后，最后寻找极大值 ；从拟合精 度最大的点构成峰的边缘 ，从而实现异常区域 的 度和计算的复杂度考虑 ，本文选择 5×5窗 口，用 圈定 。 三阶多项式进行滑动拟合 。 探讨利用二元多项式微商法 ，找到所有峰顶 ， 二元三阶多项式进行拟合的公式为 然后对每个峰从陡峭度和峰值幅度两方面进行判 f(x，3，)=C。+c1+c2+c332+c4xyq-c5Y 断 ，从而圈定异常峰 ，然后提取反映特性变化的峰 +f6 。+c7z。Y+C8xy。+f9 。 (1) 的边缘，从而圈定异常区域 。 其导数 以及偏导数分别为 [收稿 日期]2010—01—28 [基金项 目]国家 自然科学基金资助项 目 [作者简介]蒋开明(1965一)，男，工程师，研究方向：自动化 系统与可靠性、信息处理技术，E-mail：jjkkmm@163．com。 · 536 ·