钻井布局问题.doc

钻井布局问题 摘要 本文研究了钻井布局的合理设计，将井和结点抽象成坐标系上的点，依据对网格的方向规定确定两者的相对运动给出算法，可采用软件进行搜索求解或者引入0-1变量表示旧井是否被用建立最优化模型利用Lingo软件求解。 对于第一问，按照题目的规定距离且网格的横向和纵向固定，网格和旧井的相对运动为平移进行约束条件的分析，在平移网格时为了保证既不造成疏漏又要节省运算的前提下确定合适的横纵步长，经分析比较得出0.01为横纵坐标最合适的平移步长，由于网格都是边长相等的正方形，因此步长移动的范围应该不大于1，在满足以上条件下可求出最多可以利用的旧井数为4个，它们的分别为：2，4，5，10号旧井； 对于第二问，在第一问的解析基础上假如网格可以旋转时的分析与求解，因为所有象限的网格大小形状都相同，所以当网格旋转超过时就相当于重复了，因此旋转的范围定为（），规定网格每移动一次后就要旋转一次，这样既保证了查找的精度，至于旋转的弧度步长经分析定位0.01度，求得再旋转情况下可被利用的旧井为1，6，7，8，9，11。 针对以上两问，用画出平移后的旧井的位置图像与平移前的进行比较，直观的验证此次结果的正确性。 对于第三问，我们规定的两点距离为横向距离与纵向距离的最大值进行研究，得出n口旧井都可被利用的横纵坐标x，y的条件为：任意两口旧井的横坐标之差与纵坐标之差的小数部分要满足同时小于0.1。 关键词 0-1变量 最优化模型 搜索求解 lingo 1.问题的重述 设平面上有n个点,其坐标为，表示已有的n个井位。新布置的井位是一个正方形网格N的所有结点（所谓“正方形网格”是指每个格子都是正方形的网格；结点是指纵线和横线的交叉点）。假定每个格子的边长（经纬的纵横间距）都是1单位（都是100米）。整个网格是可以在平面上任意移动的。若一个已知点P与某个网格结点X的距离不超过给定误差（=0.05单位），则认为P处的旧井资料可以利用，不必在结点X处打新井。 研究如下问题： （1） 假定网格的横向和纵向是固定的（比如东西向和南北向），并规定两点间的距离为其横向距离（横坐标之差的绝对值）及纵向距离（纵坐标之差的绝对值）的最大值。在平面上平行移动网格N，使可利用的旧井数尽可能大。是提供数值计算方法，并对下面的数值例子用计算机进行计算。 在欧氏距离误差的意义下，考虑网格的横向纵向不固定（可以旋转）的情形，给出算法及计算结果。 如果有n口旧井，给出判定这些井均可利用的条件和算法（你可以任意选定一种距离）。 数值例子n=12个点的坐标如下表所示： i 1 2 3 4 5 6 7 8 9 10 11 12 a 0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50 b 2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80 2. 问题分析 根据各问中网格方向是否固定，遵照运动相对性原理，对移动形式和算法进行仔细分析。 对于问题一，在网格的横向和纵向固定、两点距离为坐标距离的最大值的前提下进行以下分析： 本题明确要求要以旧井的最大利用率为目的，因此要建立最优化模型，并引入0——1变量进行求解； 平移步长的大小关系到计算时间的长短和搜索求解的误差问题，如果步长太大误差就大，有可能无法最大数量的利用旧井，如果步长太小在实际应用会造成计算上不必要的麻烦； 平移的方法并不是单一的，可以使旧井的的坐标固定，让网格移动，也可以使网格固定平移旧井的坐标，根据与运动的相对性原则可根据两种平移方法分别给出算法； 当平移后使旧井与网格横纵坐标距离均小于0.05时，该旧井可被利用，在这里单位1为实际生活中的一百米； 对于问题二，在网格的横向和纵向不固定、两点距离为欧式距离的前提下进行以下分析： 在这一问中题目明确要求用欧式距离的尺度标准而并非第一问中横纵坐标的距离的标准来衡量旧井与网格的距离，所谓欧式距离就是指横坐标距离的平方与纵坐标距离的平方的和的平方根； 由于本问中坐标轴可以旋转，这就增加了网格与旧井可搜寻的范围，也增加了旧井可利用的机会； 尽管可以旋转，但网格和旧井之间仍然存在运动的相对性，因此可以固定一个时另一个旋转，此外旋转还有顺时针与逆时针两个方向； 为了便于计算并结合实际生活的应用，旋转和平移不能同时进行，可以先平移一步在进行旋转，也可以先旋转一定度数再平移，但是旋转的度数并须要又一个合适的值，与第一问相同，旋转弧度的大小将直接影响求解的该问题的时间和精度； 对应于第三问，我们选择问题一种的距离进行研究，题目给出的12口井的位置具有代表性，因此我们研究满足第一问中4口井的坐标位置，根据有个别到一般的原