储层提前或滞后的矢量入靶井眼轨道优化

2016-09-15罗朝东孙巧雷王旭东夏成宇王大勇

西部探矿工程 2016年4期

关键词：斜角井眼斜率

罗朝东，孙巧雷，王旭东，夏成宇，王大勇，杨　峰

（1.中石化西南石油工程有限公司钻井工程研究院，四川德阳618000；2.长江大学机械工程学院，湖北荆州434023）

储层提前或滞后的矢量入靶井眼轨道优化

罗朝东*1，孙巧雷2，王旭东1，夏成宇2，王大勇1，杨峰1

（1.中石化西南石油工程有限公司钻井工程研究院，四川德阳618000；2.长江大学机械工程学院，湖北荆州434023）

在钻定向井的过程中，按设计的具有井斜和方位变化的轨道钻进，若储层提前或滞后，就需要重新设计储层提前或滞后的轨道。针对这个问题，在满足现场工具造斜能力范围、靶点井斜控制范围和入靶点精度等情况下，建立了储层提前或滞后的井眼轨道优化的数学模型。从C#调用Matlab优化工具箱中求解多约束非线性多元函数最小值，从而建立储层提前或滞后的井眼轨道优化模块，计算储层终点改变后不同垂深各段的井眼曲率和各段的长度，指导现场设计钻具组合选取及判断是否需要上提钻杆。

储层提前滞后；轨道；优化设计；数学模型

在定向钻井［1］过程中，钻井的轨道是先期设计的，但是储层的位置往往具有不确定性。储层位置的不确性是指储层可能在水平方向和垂直方向上存在着一定的误差，即在钻进的过程中可能出现储层提前或者滞后的情形，因此需要实时监测实钻轨道，并在钻进的过程中，根据储层的位置及靶区形势［2-3］进行相应的钻进方案调整，对井眼的轨迹进行控制、对钻进的轨道进行优化设计［4-5］。

井眼轨迹控制是水平井钻井技术的关键环节［6］，国内外众多学者对钻井过程中井眼轨迹控制［7-9］的问题进行了探究，然而针对储层提前或滞后轨道优化的问题，并未建立起轨道优化的数学模型，笔者通过理论分析，建立了储层提前或者滞后轨道的钻井轨道目标函数和优化的计算模型，并通过具体算例进行了分析计算。

1　储层提前或滞后轨道优化理论模型建立

1.1建立坐标系与优化要求

在钻井过程中，储层的分布特点［10］是在水平方向尺寸大、垂直方向尺寸小，储层的提前或滞后主要考虑垂直方向的误差。在进行井眼轨道设计时，假定三段式井眼轨迹（造斜+微增+造斜）如图1所示，同时假定储层的位置和井眼方向是明确的，以P面水平方向位移为X，垂深为Z，建立二维坐标系O-XZ，由当前井底位置A和井眼方向是确定的，假定A点坐标为（XA，YA，ZA）、井斜角为αA。图中A为钻头现在所处位置，DM为储层提前或滞后当前靶点到储层顶界的垂深，MF为储层顶界到储存中间位置的距离（即储层厚度的一半），AC为优化轨道，S1为稳斜段，造斜率为Kα1；S2为微增斜段，微造斜率为Kα2；S3为造斜段，造斜率为Kα3；B、G、C处对应的井斜角依次为α1、α2、α3。

图1　轨道设计示意图

对于优化后设计的轨迹必须满足以下要求：

（1）对到达终点C的要求：C点到储层中间位置的距离小于微小量（即||F′F≤b1，如||F′F≤0.05m），确保接近储层中部。

（2）到达终点C的井斜角控制范围为（αDmin，αDmax）。

（3）优化轨道的各段造斜率在各段现场工具的造斜能力范围内（即Kαmin1≤Kα1≤Kαmax1，Kαmin3≤Kα3≤Kαmax3Kαmin2≤Kα2≤Kαmax2）

1.2建立目标函数

假设轨道的井斜变化率为Kαi，井身长为Si，则井眼轨迹的三维方程为：

通过对函数F′F（F′F为设计轨迹上的C点至靶点间的距离）的优化计算，求得在F′F足够小（如F′F＜0.05m）时的井斜变化率和段长等参数，进而可得到所求设计轨道。建立的目标函数是：寻找待钻井轨道长度S的最小值，即：

2　计算模型的建立

目前常见的三段式轨道形状包括造斜稳方位（Kα1≠0，Kφ=0）、微造斜段（Kα2≠0，Kφ=0）和造斜段（Kα3≠0，Kφ=0），下面将通过上述三段式轨道建立待钻井轨道最优化数学模型。

假设稳斜段长为S1，微造斜段长为S2，造斜段长为S3，造斜率为Kα1，微造斜率为Kα2，造斜率为Kα3。则C点的坐标应为：

其中：

根据上述的设计准则，储层提前或滞后轨道设计的最优化数学模型（记为W）：

其中：

F′F=ZF-ZF'

b1为F′与F点间的距离精度控制常数，H0为储层厚度。通过对对上述模型（W）的优化计算，可得到满足所有约束条件下最优轨道的轨道参数。

3　实例计算

根据上述的优化模型，讨论储层改变后轨迹的优化。

3.1地层参数

四川盆地川西坳陷新场地区某井的井身数据为：井底井深为2725m，现井底A点坐标为E为23.5m，N 为-77m，垂深2206m，井斜为60°，方位为168°；由于储层提前导致储层顶界变为垂深2306m，储层厚度10m，优化设计待钻井眼轨迹。现场提供工具AB段造斜率Kα1增斜变化范围为（2°/30m，9°/30m）或降斜变化范围为（-5°/30m，-2°/30m），BG段的Kα2微增变化范围（2°/30m，3°/30m），GC段的Kα3造斜率（3°/30m，9°/30m），F′F精度为0.05m。靶点C处井斜角变化范围为（89.5°，90.5°），由（7）建立的6变量S（S1，S2，S3，R1，R2，R3）优化模型为：