钻井布局研究 数学建模.docVIP

钻井布局问题研究 摘要 本文主要研究了钻井布局过程中使可利用旧井位最大化的问题，即如何移动规划中的正方行网格（边长为1）使满足与网格结点的距离不超过=0.05个单位的旧井的个数最多。文中先引入了0-1变量，旧井可利用（与结点距离不超过0.05）为1，不可利用为0主要进行了平行移动（不可旋转，只可横向、纵向移动）和自由移动（可旋转）的两方面研究。在进行平行移动的研究中两点间的距离为其横向距离（横坐标之差的绝对值）及纵向距离（纵坐标之差的绝对值）的最大值。自由移动的研究是在欧氏距离误差的意义下进行的。在解决平移问题的过程中根据运动的相对性，文中将网格的移动转换成了旧井的整体移动。 对于问题一，然后假设旧井横向移动了，纵向移动了，用取整法将旧井移动后与其最近的结点坐标表示为。画出树形图将旧井位坐标、移动后旧井位坐标、结点坐标之间、给定误差之间的关系直观化后，根据给定误差确定横向、纵向移动步长为0.01。移动范围不超过1。建立最优化模型，用Matlab搜索求解并画出点阵模型和用Lingo全局求解求出在平行移动的情况下可被利用的旧井最多了4个，它们分别为：、、、。 对于问题二，网格除在纵向和横向方向移动之外，还进行旋转，我们把原旧井坐标的其中一个作为坐标原点进行顺时针转动，即网格为逆时针转动，根据条件我们确定旋转步长为1度，旋转范围（0，），分析旧井点坐标，移动距离、旋转角度、移动后井点坐标、结点坐标的关系，建立最优化模型，再利用Metlab软件编写程序，用Matlab搜索求解并画出点阵模型，其能利用的旧井数量为6口；分别为，，，，，。 关键词：0-1变量 取整 最优化模型 Matlab搜索求解 Lingo全局求解 1问题重述 在平面上有n个井位，坐标为。现要重新布置井位，要求把新井布置在一个每个格子的边长都是1个单位（100米）的正方形网格N的所有结点（纵线和横线的交叉点）上。如果旧井P距网格结点的距离不超过=0.05个单位，则认为旧井可以利用，不必在此结点上打新井。整个网格在平面上可以任意移动。 研究以下问题： （1）假定网格的横向和纵向固定的，并规定两点间的距离为其横坐标之差的绝对值及纵坐标之差的绝对值的最大值，在平面上平行移动网格N，使可利用的旧井数尽可能大。对表中提供的12个旧井进行分析。为了方便计算我们可以把网格看做固定的，把旧井P的整个面进行平行移动。 （2）在欧氏距离误差的意义下，考虑网格的横向纵向不固定（可以旋转）的情形，给出算法及计算结果。 为了方便问题研究，我们用Matlab软件绘出移动网格前旧井的点阵模型： P12（9.50，0.80）P1（0.50，2.00）P11（8.98，3.41）P10（8.38，4.50）P P12（9.50，0.80） P1（0.50，2.00） P11（8.98，3.41） P10（8.38，4.50） P9（7.57，2.01） P8（5.43，4.10） P6（4.72，2.00） P7（4.72，6.24） P5（3.40，5.50） P4（3.37，3.51） P2（1.41，3.50） P3（3.00，1.50） 2问题分析 问题一分析 因为表中给出的所有点坐标都在第一象限，所以我们在第一象限研究问题。我们可以用直接Matlab编程求出结果后，使结果直观化，再建立模型。用Lingo全局求解，得到具体结果。因为最后结果考虑的是距离问题，根据运动的相对性。所以移动网格和整体移动井位本质上是相同的。因此我们可以考虑移动井位平面，这样我们只需要先引入变量x(横向移动距离)和y(纵向移动距离)。再引入0-1变量（距离小于0.05为1，大于0.05为0），为了使几个点的关系直观化，我们画出树图进行分析： P(i) P(i) a(i) b(i) a(i)+x b(i)+y [b(i)+y+0.5] [b(i)+y+0.5] [a(i)+x+0.5] [a(i)+x+0.5] 纵向平移Y 横向平移X 大于0.05 小于等于0.05 小于等于0.05 大于0.05 P(i)不可利用 P(i) 可 利 用 P(i)不可利用 因为问题的本质是将移动后的井位和网格的结点进行横坐标比较和纵坐标比较，而正方形网格边长为1，即所有结点坐标都是整数，所以我们可以把问题看为原井位移动后的横坐标、纵坐标与其横坐标、纵坐标进行（四舍五入）取整后的值进行比较。因为给定误差是0.05，所以为了节省时间和确保没有疏漏我们确定移动步长为0.01。 问题二分析： 问题二的研究可分两步进行，先横向、纵向移动，再旋转一定角度，因此问题二只要在问题一的基础上考虑到旋转角度，并考虑旋转角度对坐标的影响。 3假设条件 (1)问题一中的网格横向纵向绝对固定，误差远远小于0.05。 (2)在平行移动中，误差