矿井布局问题论文_b09050228(徐力).pdf

钻井布局问题的数学模型 摘要 勘探部门在某地区找矿时，首先进行初步勘探，取几个位置钻井，取得地质资料；然后进行系统 勘探，进行纵横等距的撒网式钻井。显然如果能尽可能多的在系统勘探时利用初步勘探的钻井资料， 就能有效的节约费用。 在不考虑网格方向的情况下，本文首先给出了两个结论，即网格的位置由节点唯一确定，与原始 矿井节点以及单位方格内的矿井映射点有相同的性质。这样就将问题等效为在单位方格内确定网格的 一个节点。要解决这个问题，首先我们提出运用一般的搜索法对网格节点在单位方格内进行遍历（模 型一）。通过对遍历算法进行有效的优化，大量减少了搜索的次数，进而初步计算得到了原井位最多 有 4 个可被利用，并给出了方格节点的坐标为： (0.40  X ,0.50  Y ) X ,Y   i i i i 考虑到搜索算法的复杂度，我们给出了模型二，即在单位方格内通过确定每个矿井映射节点被利 用时节点的区域，来找出方格内被这些区域覆盖次数最高的部分，显然如果将节点放在这部分内，将 会有最多的点被利用，从而也就确定了节点的位置范围。运用 MATLAB 进行计算与判别，得到最多有 4 个可被利用，并求出了网格节点坐标具体的范围： (s X ,t  Y ) 其中 s  (0.37,0.47),t  (0.41,0.51) X ,Y   i i i i 当网格方向可以改变时，我们建立了模型三。考虑到判别条件是欧氏距离，可以将原题简化为一 个圆形进行覆盖，圆的半径为 ε ，再用类比利用模型二进行判断，那么就能相应的找到最优规划。模 型三首先进行了误差分析，根据假设的误差使用夹逼法则，然后，为了减小搜索范围，我们证明了 [ ] 时，最多有 6 个矿井可被利用。 = + ԑ ‒ s , = [ + ԑ ‒ t] 对于第三问的判定算法，我们仍然根据模型三，建立假设模型四。构造出两个极端情况，此时所 有矿井均可被利用。具体算法的见问题三分析步骤。 最后我们对模型四的一个假设进行了检验。虽然这个假设严格的说并不成立，但通过我们用蒙特 卡罗方法进行多次模拟，发现假设成立的概率极高。综上，我们可以先用模型四进行计算，再对结果 进行检验，对极少数不成立的，可以综合特殊情况进行考虑。 关键词 勘探 矿井 遍历算法 蒙特卡罗 数值分析 误差分析 假设检验 一． 问题重述 勘探部门在某地区找矿。初步勘探时期已零散地在若干位置上钻井，取得了地质 料。进入系统勘探时期后，要在一个区域内按纵横等距的网格点来布置井位，进行 “撒网 式”全面钻探。由于钻一口井的费用很高，如果新设计的井位与原有井位重合（或相当接 近），便可利用旧井的地质资料，不必打这口新井。因此，应该尽量利用旧井，少打新井， 以节约钻探费用。比如，钻一口新井的费用为500 万元，利用旧井资料的费用为 10 万元， 则利用一口旧井就节约费用 490 万元。 设平面上有 n 个点 P ，其坐标为(a , b ),i=1,2,…,n，表示已有的n 个井位。新布 i i i 置的井位是一个正方形网格 N 的所有结点（所谓 “正方形网格”是指每个格子都是正方形 的网格；结点是指纵线和横线的交叉点）。假定每个格子的边长（井位的纵横间距）都 是 1 单位（比如 100 米）。整个网格是可以在平面上任意移动的。若一个已知点 P 与某 i 个网格结点 X 的距离不超过给定误差 （＝0.05 单位），则认为 P 处的旧井资料可以利 i