浅析钻井布局问题的数学模型

钻井布局问题的数学模型

摘要

勘探部门在某地区找矿时，首先进行初步勘探，取几个位置钻井，取得地质资料；然后进行系统勘探，进行纵横等距的撒网式钻井。显然假如能尽可能多的在系统勘探时利用初步勘探的钻井资料，就能有效的节约费用。

在不考虑网格方向的情况下，本文首先给出了两个结论，即网格的位置由节点唯一确定，与原始矿井节点以及单位方格内的矿井映射点有相同的性质。如此就将问题等效为在单位方格内确定网格的一个节点。要解决那个问题，首先我们提出运用一般的搜索法对网格节点在单位方格内进行遍历（模型一）。通过对遍历算法进行有效的优化，大量减少了搜索的次数，进而初步计算得到了原井位最多有4个可被利用，并给出了方格节点的坐标为：

Z ∈++i i i i Y X Y X , )50.0,40.0(

考虑到搜索算法的复杂度，我们给出了模型二，即在单位方格内通过确定每个矿井映射节点被利用时节点的区域，来找出方格内被这些区域覆盖次数最高的部分，显然假如将节点放在这部分内，将会有最多的点被利用，从而也就确定了节点的位置范围。运用MATLAB 进行计算与判不，得到最多有4个可被利用，并求出了网格节点坐标具体的范围：

Z ∈∈∈++i i Y ,X )51.0,41.0(),47.0,37.0( 其中 ),s (t s Y t X i i

当网格方向能够改变时，我们建立了模型三。考虑到判不条件是欧氏距离，能够将原题简化为一个圆形进行覆盖，圆的半径为ε，再用类比利用模型二进行推断，那么就能相应的找到最优规划。模型三首先进行了误差分析，依照假设的误差使用夹逼法则，然后，为了减小搜索范围，我们证明了

时，最多有6个矿井可被利用。

关于第三问的判定算法，我们仍然依照模型三，建立假设模型四。构造出两个极端情况，现在所有矿井均可被利用。具体算法的见问题三分析步骤。

最后我们对模型四的一个假设进行了检验。尽管那个假设严格的讲并不成立，但通过我们用蒙特卡罗方法进行多次模拟，发觉假设成立的概率极高。综上，我们能够先用模型四进行计算，再对结果进行检验，对极少数不成立的，能够综合专门情况进行考虑。

关键词

勘探矿井遍历算法蒙特卡罗数值分析误差分析假设检验

一．问题重述

勘探部门在某地区找矿。初步勘探时期已零散地在若干位置上钻井，取得了地质资料。进入系统勘探时期后，要在一个区域内按

纵横等距的网格点来布置井位，进行“撒网式”全面钻探。由于钻

一口井的费用专门高，假如新设计的井位与原有井位重合（或相当

接近），便可利用旧井的地质资料，不必打这口新井。因此，应该

尽量利用旧井，少打新井，以节约钻探费用。比如，钻一口新井的

费用为500万元，利用旧井资料的费用为10万元，则利用一口旧井

就节约费用490万元。

设平面上有n个点P i，其坐标为(a i,b i),i=1,2,…,n，表示已有的n个井位。新布置的井位是一个正方形网格N的所有结点（所谓“正方形网格”是指每个格子差不多上正方形的网格；结点是指纵线和横线的交叉点）。假定每个格子的边长（井位的纵横间距）差不多上1单位（比如100米）。整个网格是能够在平面上任意移动的。若一个已知点P i与某个网格结点X i的距离不超过给定误差 （＝0.05单位），则认为P i处的旧井资料能够利用，不必在结点X i 处打新井。

为进行辅助决策，勘探部门要求我们研究如下问题：

1)假定网格的横向和纵向是固定的（比如东西向和南北向），并规定两点间的距离为其横向距离（横坐标之差绝对值）及纵向距离（纵坐标之差绝对值）的最大值。在平面上平行移动网格N，使可利用的旧井数尽可能大。试提供数值计算方法，并对下面的数值例子用计算机进行计算。

2)在欧氏距离的误差意义下，考虑网格的横向和纵向不固定（能够旋转）的情形，给出算法及计算结果。

3)假如有n口旧井，给出判定这些井均可利用的条件和算法（你能够任意选定一种距离）。

n=12个点的坐标如下表所示：

二．符号讲明

(原始勘测中任意一点

将旧井移至单位正方形里面时候的新井

f 可利用的最大钻井数

n 正方形的候选位置的数目

旧井离新表格节点的距离

Q（s,t）网格的一个节点

推断旧井节点是否在表格范围内的标志位

a,b等原始节点需要转化至单位正方形内部时候需要移动的整数大小。

三．模型假设

1. 初步勘探时，所取的点比较分散，不存在两个点的资料被系统

勘探时同一个点利用的情况。

2. 系统勘探所取的点数比初步勘探时取的点多。

3. 初步勘探所取的点在系统勘探的勘探范围之内。

4. 假定每个格子的边长（井位的纵横间距）差不多上1单位（比

如100米）。整个网格是能够在平面上任意移动的。

5. 若一个已知点P i与某个网格结点X i的距离不超过给定误差ε

（=0.05单位），则认为P i处的旧井资料能够利用，不必在结点X i处打新井。

四．问题的分析及讲明

考虑到勘探部门在某地区找矿时是分为两步进行的，即：第一步，初步勘探，取得地质资料。第二步，系统勘探，进行“撒网式”全面钻探。那么，为了节约钻探费用，我们自然希望充分利用第一步的数据，来减少勘探次数。

对实际问题进行分析，能够认为初步勘探时钻井的位置可对应与二

维坐标系中的点，成为初始点。系统勘探时的全面钻探可认为是二维坐标系中横坐标与纵坐标间距相等的点为节点所构成的网格。如此，问题的实质就转化为如何定位网格，使尽可能多的初始点位于网格节点的误差范围内。

为了更好的表述问题，下面给出两条结论：

结论一：假定网格的横向和纵向是固定的，那么只要确定其中一个节点，就可确定整个网格。

证明：专门显然，当网格节点中的一个点的坐标)s,(t 确定时，其他点的坐标可由如下公式确定：

Z ∈++i i i y y t x s ,x ),(i

结论二：关于坐标系里的每一个点),(i i b a ，能够定义其映射点⎣⎦⎣⎦),(i i i i b b a a --，映射点与原来的点对网格有相同的位置关系。 证明：⎣⎦i a 意为对i a 向下取整。由于网格横向和纵向的单位差不多上

1，因此当点),(i i b a 位于网格中一个节点),s (t 的误差范围内时，则其映射点必定位于⎣⎦⎣⎦),(i i i i b b a a --的误差范围内。反之假如点),a (i i b 不在任何一个节点的误差范围内，那么映射点也不在任何一个节点的误差范围内。

五． 模型的建立