矿井布局问题论文

1、.：.；钻井规划问题的数学模型摘要勘探部门在某地域找矿时，首先进展初步勘探，取几个位置钻井，获得地质资料；然后进展系统勘探，进展纵横等距的撒网式钻井。显然假设能尽能够多的在系统勘探时利用初步勘探的钻井资料，就能有效的节约费用。 在不思索网格方向的情况下，本文首先给出了两个结论，即网格的位置由节点独一确定，与原始矿井节点以及单位方格内的矿井映射点有一样的性质。这样就将问题等效为在单位方格内确定网格的一个节点。要处理这个问题，首先我们提出运用普通的搜索法对网格节点在单位方格内进展遍历模型一。经过对遍历算法进展有效的优化，大量减少了搜索的次数，进而初步计算得到了原井位最多有4个可被利用，并给出了方格

2、节点的坐标为：思索到搜索算法的复杂度，我们给出了模型二，即在单位方格内经过确定每个矿井映射节点被利用时节点的区域，来找出方格内被这些区域覆盖次数最高的部分，显然假设将节点放在这部分内，将会有最多的点被利用，从而也就确定了节点的位置范围。运用MATLAB进展计算与判别，得到最多有4个可被利用，并求出了网格节点坐标详细的范围：当网格方向可以改动时，我们建立了模型三。思索到判别条件是欧氏间隔 ，可以将原题简化为一个圆形进展覆盖，圆的半径为，再用类比利用模型二进展判别，那么就能相应的找到最优规划。模型三首先进展了误差分析，根据假设的误差运用夹逼法那么，然后，为了减小搜索范围，我们证明了Xi=ai+-s

3、,Yi=bi+-t时，最多有6个矿井可被利用。对于第三问的断定算法，我们依然根据模型三，建立假设模型四。构造出两个极端情况，此时一切矿井均可被利用。详细算法的见问题三分析步骤。 最后我们对模型四的一个假设进展了检验。虽然这个假设严厉的说并不成立，但经过我们用蒙特卡罗方法进展多次模拟，发现假设成立的概率极高。综上，我们可以先用模型四进展计算，再对结果进展检验，对极少数不成立的，可以综合特殊情况进展思索。关键词 勘探 矿井 遍历算法 蒙特卡罗 数值分析 误差分析 假设检验问题重述勘探部门在某地域找矿。初步勘探时期已零散地在假设干位置上钻井，获得了地质资料。进入系统勘探时期后，要在一个区域内按纵横等

4、距的网格点来布置井位，进展“撒网式全面钻探。由于钻一口井的费用很高，假设新设计的井位与原有井位重合或相当接近，便可利用旧井的地质资料，不用打这口新井。因此，应该尽量利用旧井，少打新井，以节约钻探费用。比如，钻一口新井的费用为500万元，利用旧井资料的费用为10万元，那么利用一口旧井就节约费用490万元。 设平面上有n个点Pi，其坐标为(ai, bi),i=1,2,n，表示已有的n个井位。新布置的井位是一个正方形网格N的一切结点所谓“正方形网格是指每个格子都是正方形的网格；结点是指纵线和横线的交叉点。假定每个格子的边长井位的纵横间距都是1单位比如100米。整个网格是可以在平面上恣意挪动的。假设一

5、个知点Pi与某个网格结点Xi的间隔 不超越给定误差0.05单位，那么以为Pi处的旧井资料可以利用，不用在结点Xi处打新井。为进展辅助决策，勘探部门要求我们研讨如下问题： 1)假定网格的横向和纵向是固定的比如东西向和南北向，并规定两点间的间隔 为其横向间隔 横坐标之差绝对值及纵向间隔 纵坐标之差绝对值的最大值。在平面上平行挪动网格N，使可利用的旧井数尽能够大。试提供数值计算方法，并对下面的数值例子用计算机进展计算。 2)在欧氏间隔 的误差意义下，思索网格的横向和纵向不固定可以旋转的情形，给出算法及计算结果。 3)假设有n口旧井，给出断定这些井均可利用的条件和算法他可以恣意选定一种间隔 。 n=1

6、2个点的坐标如下表所示： i123456789101112ai0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50 bi2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80 符号阐明Pi(ai,bi) 原始勘测中恣意一点Pi(ai，,bi) 将旧井移至单位正方形里面时候的新井Pif可利用的最大钻井数n 正方形的候选位置的数目 旧井离新表格节点的间隔 Qs,t 网格的一个节点ui 判别旧井节点能否在表格范围内的标志位Xi ，Yi a,b等原始节点需求转化至单位正方形

7、内部时候需求挪动的整数大小。 模型假设1. 初步勘探时，所取的点比较分散，不存在两个点的资料被系统勘探时同一个点利用的情况。2. 系统勘探所取的点数比初步勘探时取的点多。3. 初步勘探所取的点在系统勘探的勘探范围之内。4. 假定每个格子的边长井位的纵横间距都是1单位比如100米。整个网格是可以在平面上恣意挪动的。5. 假设一个知点Pi与某个网格结点Xi的间隔 不超越给定误差=0.05单位，那么以为Pi处的旧井资料可以利用，不用在结点Xi处打新井。 问题的分析及阐明思索到勘探部门在某地域找矿时是分为两步进展的，即：第一步，初步勘探，获得地质资料。第二步，系统勘探，进展“撒网式全面钻探。那么，为了

8、节约钻探费用，我们自然希望充分利用第一步的数据，来减少勘探次数。对实践问题进展分析，可以以为初步勘探时钻井的位置可对应与二维坐标系中的点，成为初始点。系统勘探时的全面钻探可以为是二维坐标系中横坐标与纵坐标间距相等的点为节点所构成的网格。这样，问题的本质就转化为如何定位网格，使尽能够多的初始点位于网格节点的误差范围内。为了更好的表述问题，下面给出两条结论：结论一：假定网格的横向和纵向是固定的，那么只需确定其中一个节点，就可确定整个网格。证明：很显然，当网格节点中的一个点的坐标确定时，其他点的坐标可由如下公式确定：结论二：对于坐标系里的每一个点，可以定义其映射点，映射点与原来的点对网格有一样的位置

9、关系。证明：意为对向下取整。由于网格横向和纵向的单位都是1，所以当点位于网格中一个节点的误差范围内时，那么其映射点必然位于的误差范围内。反之假设点不在任何一个节点的误差范围内，那么映射点也不在任何一个节点的误差范围内。模型的建立模型一 遍历法思索曾经给定点矿井以及网格的横向与纵向固定，那么问题就变为如何确定出网格的最优平移位置。由结论一，定出网格的位置只需确定出一个节点，同时，由于网格以1为单位，那么在单位方格内，必然有且只需一个节点，这样我们在四个点所构成的方块单位方格内对节点进展搜索，就能将网格的全部能够的平移情况进展遍历。思索到标题中所给数据准确到0.01，那么可取0.01为两个坐标的步

10、出息展搜索。对于搜索到的每一个坐标，作为网格的节点，来确定整个网格，再计算究竟有多少初始点能进入网格节点的误差范围，最后进展比较，选出包容初始点最多的网格作为处理方案。搜索法从实际上是可行的，然而这种方法的计算量往往比较大。如上例，要搜索的节点个数为即循环次数：10010010000次而每一次循环都要对12个点依次进展判别，这必然非常复杂，当实践的数据量比较大，数据精度要求比较高时，这种直接搜索的方法计算起来就非常耗时，甚至是不可行的。为了减少循环次数，必需想方法对算法进展优化。这里我们思索对于每一个初始点，由结论二，首先在单位方格内找到该点的映射点，这样就将一切点都映射到了单位方格内。然后进

11、展搜索，在映射点附近搜索网格节点，使映射点位于网格节点的误差范围内，再对节点构成的网格进展分析、比较，选出其中的最优。经过这样的过程，循环的次数为：1111121452次这样就大大减少了循环次数，优化了搜索算法，使之耗时短，使方案变得可行。 运用优化后的搜索算法，对所给数据进展计算，得到最多可以有4口井不用打而可以利用初步勘探的资料。可以利用的初始点为： (1.41,3.50) (3.37,3.51) (3.40,5.50) (8.38,4.50)网格坐标为：模型二 框图分析法思索到遍历节点法毕竟是一种复杂度极高的算法，当数据量大的时候就会有严重缺陷，因此这里给出一种简单易行而又非常准确的算法

12、：方框图分析法。 首先，根据定理二，确定出初始点在单位方格内的坐标即映射点，变换前后点的坐标如下：表1：坐标映射表n123456789101112初始点ai0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50 bi2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80 变换点ai0.50 0.41 0.00 0.37 0.40 0.72 0.72 0.43 0.57 0.38 0.98 0.50 bi0.00 0.50 0.50 0.51 0.50 0.00 0.2

13、4 0.10 0.01 0.50 0.41 0.80 运用matlab软件将变换前的点画在一坐标系下，如图1：点图：1同时将映射之后的点画于一个坐标系下面得到如图2：映射点图：2将一切的旧井节点平移至单位正方形可以得到Q=x,y|0=x=1,0=yPiai,bi ai=ai-ai, bi=bi-bi在上述变换后，问题1大致等价于用一个变长为2的正方形去覆盖尽能够多的Pi,正方形的中心就是网格的一个节点所在位置。如图-3所示：察看可得，正方形有n2个候选位置。假设正方形左边经过Pi，右边经过Pj，那么aj【ai, ai+2*】,bj【bi-2*,bi】代入数据得：aj【ai, ai+0.1】,b

14、j【bi-0.1,bi】1根据上表1可知：ai=【0.5,0.41,0.00,0.37,0.40,0.72,0.72,0.43,0.57,0.38,0.98,0.50】 bi=【0.00,0.50,0.50,0.51,0.50,0.00,0.24,0.10,0.01,0.50,0.41,0.80】;将以上数据代入1式得：aj0.5,0.6,0.41,0.51,0.00,0.10,0.37,0.47,0.40,0.50,0.72,0.82,0.72,0.82,0.43,0.53,0.57,0.67,0.38,0.48,0.98,1.08,0.50,0.60; bj-0.1,0.00,0.40,0

15、.50,0.40,0.50,0.41,0.51,0.40,0.50,-0.10,0.00,0.14,0.24,0.00,0.10, -0.09,0.01,0.40,0.50,0.31,0.41,0.70,0.80;经过将上述数据运用matlab数学软件分析可以得到，满足条件的旧井数目为4，方格节点的纵横坐标为0.42,0.46，亦即0.42+xi,0.46+yi, xi，yiZ模型三 动态分析法网格平行挪动，对网格上面的任一点Pi，当网格挪动整数个单位时，Pi相对于最近的网格节点的距离不变， Pi在网格上平移整数个单位时，其相对于最近的网格节点的间隔 不变。那么如今把一切旧井节点平移至正方形：

16、Q=(x,y)|0=x=1,0=y=1;设网格的一个节点为(s,t)，0=s,t设Pi:ui=1否那么ui=0；那么以上可以归纳为：2等价于以某一角度为步长旋转网格，一个半径为的圆覆盖更多的Pi，圆心是网格的一个节点位置。旋转方式：如下图可旋转的网格图4：图4Pi可利用当且仅当3最多可以利用旧井位为max1nui,其中，ui=10，。代入数据，即(s+xi-ai)2+t+yi-bin=0.0025xi,yi 化简该算式得：(-s-xi+ai)2+-t-yi+bin a=0.5 1.41 3.00 3.37 3.40 4.72 5.43 7.57 8.38 8.98 9.50; b=2.00 3

17、.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80; a=0.5 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50; plot(a,b,r) plot(a,b,*) hold on grid mkdir chapter_1999 cd chapter_1999 a=0.5 0.41 0.00 0.37 0.40 0.72 0.72 0.43 0.57 0.38 0.98 0.50; b=0.00 0.50 0.50 0.51 0.50 0.00 0.24 0.10 0.01 0.

18、50 0.41 0.80; plot(a,b,*) grid on plot(a,b,*) grid on a=0.5 0.41 0.00 0.37 0.40 0.72 0.72 0.43 0.57 0.38 0.98 0.50; b=0.00 0.50 0.50 0.51 0.50 0.00 0.24 0.10 0.01 0.50 0.41 0.80; c=a=0.37&a d=b=0.41&b avg_x=(0.37+0.47)/2avg_x = 0.4200 avg_y=(0.51+0.41)/2avg_y = 0.4600第二问：C言语代码：#include#include#define kesi 0.05int main()double a12=0.5,1.41,3.00,3.37,3.40,4.72,5.43,7.57,8.38,8.98,9.50;double b12=2.00,3.50,1.50,3.51,5.50,2.00,6.24,4.10,2.01,4.50,3.41,0.80;double c12,d12,u_112,u_212,s,t;int i=0,flag=0,k=0,j=0,p=0,x10;while(i12)ci=fabs(ai-floor(ai);di=fa