3.MORPAS-信息综合与有利度计算.ppt

MORPAS3 0信息综合与有利度计算 陈建国中国地质大学 武汉 jgchen 建模的基本方法 知识驱动 无须训练区 点 专家系统模型参数由专家确定Fuzzylogic Dempster Shaferbelieftheory PROSPECTOR数据驱动 必须有训练区 点 统计方法根据训练区计算确定模型参数证据权法 特征分析 逻辑回归 神经网络 分形方法 矿产资源区域评价数学模型 矿产资源评价强调了数学模型的应用 所涉及的大致有以下几种数学模型 统计模型 特别是多元统计模型 如因子分析 对应分析 判别分析 回归分析 典型相关分析等 利用这类模型 可对用于评价的数据从各个角度进行统计分析 除对连续型变量进行分析外 多半还有其离散分支 纯代数模型 如特征分析 主成分分析 聚类分析 趋势分析 逻辑信息法等这类模型仅立足于矩阵元算 一般无统计能力 其他数学模型 如立足于图论 模糊集 数字滤波 信息论等数学分支的数学模型 这类模型有时也含有统计概念 专家系统模型 其数学支持主要是逻辑数学模型 用于系统的知识表达与推理过程 其中也会使用到广泛的数学模型 如统计模型 模糊集 人工神经网络等 矿产资源潜力评价指南 地调局2001 5附件三 矿产资源区域评价数学模型 矿产资源评价需要解决三个问题 1 评价的数学模型要立足于数据的空间相关性 2 系统要能直接调用和操作空间图形 3 评价系统要以数据管理为核心 实现数据驱动 GIS的应用 使后面两个问题的解决成为可能 矿产资源潜力评价指南 地调局2001 5附件三 建模的基本步骤 数据处理 原始数据 图件 综合 预测结果 权重计算 优化 模型 信息综合方法 证据权法信息量法BP神经网络ART神经网络特征分析法 证据权法 证据权法 WeightsofEvidence 是一种定量评价方法 它最初用于医疗诊断 20世纪80年代末 加拿大地调所FritsAgterberg和GraemeBonham Carter将该方法引入到GIS支持下的矿产资源潜力评价中 其基本过程是 将每一种地学信息视为成矿预测的一个证据因子 通过分析 计算出每一个证据因子对成矿贡献的权重值 最终对各证据因子加权求和 得出成矿有利度值 从而对矿产远景区进行定位预测 证据权法 假设研究区被划分成面积相等的 个单元 其中有 个单元为有矿单元 则随机选取一个单元有矿的概率是 T则先验几率 O D 1 D T D 对于任意一个证据因子 其权重定义为 为因子存在区的单元数 为因子不存在区的单元数 表示有矿 表示无矿 分别为证据因子存在区和不存在区的权重值 即成矿关联度 对于原始数据缺失区域权重值为 定义关联度显著性指标C为 C W W 证据权法 证据权重法要求各证据因子之间相对于矿点分布满足条件独立 对于 个证据因子 若它们都关于矿点条件独立 则研究区内任一 单元为有矿的后验几率为 由此可得出研究区内任一 单元为有矿单元的后验概率为 最后根据后验概率圈出找矿远景区 证据权法 工作步骤 各类图件证据层 面积计算 点 区分析点 线分析 矿点密度分析 权重计算W W 根据C值进行优化 证据加权求和 先验概率 成矿有利度 后验概率 证据权法 工作步骤 1 权重计算 研究区 某类地层 落在该地层上的的矿点数不落在该地层的矿点数出露该地层的面积未出露该地层的面积 证据权法 工作步骤 1 权重计算 15070 7 0 5 1 4ln 1 4 0 3325030 3 0 5 0 6ln 0 6 0 51Total10010 2 1 类别面积矿点数矿点 面积权重 证据权法 工作步骤 1 权重计算 1MissingData 25 3 0 026030 3 0 60 0 5ln 0 5 0 6931540 4 0 15 2 7ln 2 7 0 98Total7510 1 2 3 类别面积矿点数矿点 面积权重 证据权法 工作步骤 2 权重优化 显著度曲线C 证据存在 证据不存在 断层BUFFER距离或物化探异常范围 证据权法 工作步骤 3 证据加权求和 图层加权 有利度图 最小有利度 最大有利度 证据权法 工作步骤 3 证据加权求和 证据层1证据层2加权和A110 33 0 0 0 33B120 33 0 69 0 36C130 33 0 98 1 31E22 0 51 0 69 1 20 A B C E 2 1 1 3 证据层1 证据层2 图层综合 2 信息量法 信息量计算法也属于BAYES统计分析方法 其实质是用信息量的大小来评价地质因素 标志与研究对象的关系密切程度 该方法应用于区域矿产预测 是由 维索科奥斯特罗夫斯卡娅 及 恰金 先后提出的 其基本过程是 将每一种地学信息视为成矿预测的一个证据因子 通过分析 计算出每一个证据因子对成矿贡献的权重值 最终对各证据因子加权求和 得出成矿有利度值 从而对矿产远景区进行定位预测 概述 信息量法 ln 其中 为 标志j状态提供事件 有矿 发生的信息量 为 标志j状态存在条件下事件 实现的概率 为事件 发生的概率 根据概率乘法定理 上式可变为 ln 具体计算时 总体概率用样本频率来估计 ln ln 其中 为具有标志值 的含矿单元数 为研究区中含矿单元总数 为具有标志值 的单元数 为研究区单元总数 公式 信息量法 计算步骤是 1 选择参与运算的变量集与研究对象2 计算各地质因素 找矿标志所提供的找矿信息量 定量地评价各地质因素和标志对指导找矿的作用 2 计算每个单元中各标志信息量的总和 其大小反映了该单元相对的找矿意义 用以评价找矿远景区进行预测 特征分析法 特征分析 Botbol 1971 是一种多元统计分析方法 在矿产资源定位预测中 常采用它来圈定预测远景区 原理 通过研究模型单元的控矿变量特征 查明变量之间的内在联系 确定各个地质变量的成矿和找矿意义 建立起某种类型矿产资源体的成矿有利度类比模型 然后将模型应用到预测区 将预测单元与模型单元的各种特征进行类比 用它们的相似程度表示预测单元的成矿有利度 并据此圈定出有利成矿的远景区 特征分析方法要求自变量必须是二态或三态变量 特征分析法 设有n个取样单元 已知样品 m个地质变量 第j个变量在第i个单元上的取值为xij yi为第i个单元的成矿有利度 设它们之间满足如下线性关系 Y XA其中aj j 1 2 m 就是变量的权系数 它反映了第j个变量的重要性 数学模型 特征分析法 关键问题是如何确定变量权aj 这通常有三种方法 1 平方和法2 乘积矩阵主分量法求乘积矩阵的特征值与特征向量 其中最大特征值所对应的特征向量作为变量权系数3 概率矩阵主分量法 数学模型 BP人工神经网络 黑箱 有导师引导学习 训练 BP神经网络 BP神经网络的基本结构与学习算法一个典型的BP神经网络由三层构成 即 输入层 隐含层 或者称中间层 由一层或多层组成 和输出层 各阶层之间实行全互连接方式 分导师知识学习训练和模式识别决策两个过程 导师知识学习训练过程归结起来又分为 模式顺传播 误差逆传播 导师知识记忆训练 学习收敛4个步骤 模式顺传播误差逆传播网络学习训练全局误差迭代收敛 顺传过程 网络初始化 给各连接权 Wij Wkp 及阀值 bk p 赋 1 1 间随机值 选取或建立用于网络训练的示教学习模式对Ai xi1 xi2 xim Yi yi1 yi2 yip 提供给网络 11111 000示教模式对111100 750示教模式对211010 750示教模式对310110 750示教模式对401110 750示教模式对511000 500示教模式对610010 500示教模式对701100 500示教模式对800110 500示教模式对910000 250示教模式对1001000 250示教模式对1100100 250示教模式对1200010 250示教模式对1300000 000示教模式对14 ART神经网络 竞争型网络 ART神经网络模型 1 竞争型网络与模式分类竞争型神经网络是一种以无导师示教方式进行训练的网络 只给网络提供一些学习样本 而不提供相应的希望输出模式 网络通过自身的训练 自动对输入样本进行分类 竞争型神经网络经常作为基本的网络形式 构成其它一些具有自组织能力的网络 具有代表性的有 自适应共振理论 AdaptiveResonanceTheory ART 自组织特征映射 Self OrganizingFeatureMap SOM 对向传播 CounterPropagation CP 网络等 二 主要模型与基本算法 2 竞争型网络结构及学习算法 二 主要模型与基本算法 3 竞争型神经网络常用分类算法 二 主要模型与基本算法 4 ART神经网络模型自适应共振 AdaptiveResonanceTheory 简写为ART 是一种无导师学习网络 由美国Boston大学的格罗斯伯格 StephenGrossberg 1976年提出来的 实际上是一种依据特征参数对样本归类的网络分类器这一理论充分利用了生物神经细胞之间自兴奋与侧抑制的动力学原理 让输入模式通过网络双向连接权的识别与比较 达到共振来完成对自身的记忆 并以同样的方法实现网络回想 当提供给网络回想的是一个网络中已输出向量Bk b1k b2k bmk 输出中已记忆的 或者是与已记忆的模式十分相似的模式时 网络会把这个模式回想出来 提出正确的分类 如果提供给网络回想的是一个网络记忆中不存在的模式 则网络将在不影响已有记忆的前提下 将这一模式记忆下来 并将分配给一个新的分类单元作为这一记忆模式的分类标志 二 主要模型与基本算法 5 ART神经网络模型 二 主要模型与基本算法 二 主要模型与基本算法 二 主要模型与基本算法 总结 证据权法的优点 充分考虑了某一地学因子有利及不利成矿两个方面 信息量法只考虑了有