钻井布局优化模型.docVIP

钻井布局优化模型 摘要 勘探部门于找矿初期钻井取地质资料，系统勘探时期需在一个区域内按纵横等距的网格点来布置井位。由于钻一口井的费用很高，而当新设计井位与原井位重合（或相当接近）时，不必打井。因此，合理安排井位可以节约钻探费用。 本文将旧井的利用归结为0、1规划模型，从而建立目标函数。由每个旧井到网格节点的距离不超过给定误差W（=0.05单位），可得到约束条件的不等式，以及钻井布局的非线性规划模型。对于问题1，利用MATLAB编程求解，采用全局搜索逐个判断网格平移过程中满足约束条件的旧井个数，得出的结论是最多有4口旧井可使用，编号为2、4、5、10，用Lingo验证其正确；对于问题2，采用全局搜索、坐标变换、矢量旋转判断满足约束条件的旧井个数，得出最多有6口旧井可使用，编号为1、6、7、8、9、11。 关键词：0、1规划模型；非线性规划模型；全局搜索；坐标变换；矢量旋转1.问题描述 勘探部门在某地区找矿。初步勘探时期已零散地在若干位置上钻井，取得了地质资料。进入系统勘探时期后，要在一个区域内按纵横等距的网格点来布置井位，进行“撒网式”全面钻探。由于钻一口井的费用很高，如果新设计的井位与原有井位重合（或相当接近），便可利用旧井的地质资料，不必打这口新井。因此，应该尽量利用旧井，少打新井，以节约钻探费用。比如钻一口新井的费用为500万元，利用旧井资料的费用为10万元，则利用一口旧井就节约费用490万元。 设平面上有个点Pi其坐标为（，），=1，2…表示已经有的个井位。新布置的井位是一个正方形网格N的所有结点（所谓“正方形网格”是指每个格子都是正方形的网格；结点是指纵线和横线的交叉点）。假定每个格子的边长（井位的纵横间距）都是一单位（比如100米）整个网络是可以在平面上任意移动的。若一个已知点Pi与某个网络结点Xi的距离不超过给定误差（=0.05单位），则认为Pi处的旧井资料可以利用，不必在结点Xi处打新井。 为进行辅助决策，勘探部门要求我们研究如下问题： 1）假定网格的横向和纵向是固定的（比如东西向和南北向）并规定两点间的距离为横向距离（横坐标之差绝对值）及纵向距离（纵坐标之差绝对值）的最大值。在平面上平行移动网格N使可利用的旧井数尽可能大。试提供数值计算方法，并对下表1.1的数值例子用计算机进行计算。 2）在欧式（直线距离）距离的误差意义下，平移加旋转，可以转动的情形，给出算法和计算机结果。 表1.1 旧井坐标表 i 1 2 3 4 5 6 7 8 9 10 11 12 Ai 0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50 Bi 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.模型分析 题目要求在网格的横向和纵向固定（比如东西向和南北向）且两点间的距离为横向距离（横坐标之差绝对值）及纵向距离（纵坐标之差绝对值）的最大值的条件下，在平面上平行移动网格N使可利用的旧井数尽可能大。 通过分析目标，知道钻探费用的高低与使用旧井的数目直接相关，因此目标函数是使用旧井的数目最多。针对此，本文建立了一个非线性规划模型。通过进一步分析，我们可以得到非线性规划的约束条件，具体分析如下。 已知当新设计的井位与原有井位重合（或相当接近）时，才可利用旧井的地质资料，不必打这口新井。新布置的井位是一个正方形网格N的所有结点，整个网络可以在平面上任意移动。只有已知井位与某个网络结点的距离不超过给定误差W（=0.05单位）时，该井位才可以被使用。据此，设出旧井位是否可用的0、1变量和网格平移变量，可得出模型的约束条件。 对于问题1，可利用MATLAB编程全局搜索每个旧井是否满足约束条件，来确定网格平移过程中最多可使用旧井个数。得出的结果可以用Lingo加以验证。 对于问题2，要求在欧式（直线距离）距离的误差意义下，平移和旋转网格，使可利用的旧井数最多，此时只能使用MATLAB编程的方法求解。 3.模型假设 （1）假设每个格子的边长（井位的纵横间距）都是一单位（比如100米），整个网络是可以在平面上任意移动的。 （2）假设网格的横向和纵向是固定的（比如东西向和南北向）并规定两点间的距离为横向距离（横坐标之差绝对值）及纵向距离（纵坐标之差绝对值）的最大值。 （3）假设Pi与某个网络结点Xi的距离不超过给定误差（=0.05单位），则认为Pi处的旧井资料可以利用，不必在结点Xi处打新井。 4.符号说明 （1）表示原有井位横坐标（东西方向坐标，） （2）表示原有井位纵坐标（南北方向坐标，） （3）表示网格东向偏移量（） （4）表示网格北向偏移量（） （5）表示0、1变量，0表示第口旧井不被使用，1表示第口旧井被使