本发明的目的:克服
现有技术的不足,提出一种基于神经网络的钻井测斜 仪方位角误差补偿方法,在整个测量范围内实现对钻井测斜仪方位角误差的补 偿,输出准确和高精度的钻井测斜仪方位角。
本发明的技术解决方案是:一种基于神经网络的钻井测斜仪方位角误差补 偿方法,其特征在于:利用神经网络的非线性映射功能,选择三输入单输出的 深层前向神经网络,建立钻井测斜仪输出的倾斜角、工具面角和初步计算的方 位角同真实方位角间的映射关系,主要实现步骤有:
(1)根据钻井测斜仪方位角输出
信号的特点确定神经网络的结构模型,建 立神经网络模型时以钻井测斜仪输出的倾斜角、工具面角和初步计算的方位角 做为神经网络的输入量,真实方位角为输出量来构建三输入单输出的神经网络 系统,即选择由
输入层、隐层1、隐层2和
输出层组成的四层结构的前向神经 网络;
神经网络理论证明可以用一个三层前向网络来逼近任意一个非线性映射, 但对于比较复杂的函数用深层前向网络时会取得更好的效果。因此,本发明选 择四层(即输入层、隐层1、隐层2、输出层)结构的前向神经网络。对于钻 井测斜仪来说,钻井测斜仪的倾斜角、方位角和工具面角对方位角的计算都有 影响,会产生方位角计算误差。因此,选择初步计算的的倾斜角、方位角和工 具面角作为神经网络的输入,修正补偿后的方位角作为神经网络的输出,即神 经网络输入层
节点数为3,输出层节点数为1,隐层1节点数可选为4~8,隐 层2节点数可选为4~8。
(2)获取神经网络的训练样本,模拟钻井测斜仪的倾斜角、工具面角和方 位角是在三轴位置转台上进行的,通过三轴位置转台绕其三个转动轴的转动来 模拟钻井测斜仪的倾斜角、工具面角和方位角,同时采集测斜仪输出的原始数 据并计算出倾斜角和初步计算方位角,形成神经网络的学习样本。
采用
人工神经网络解决实际问题时,主要的工作是收集样本数据。本发明 中学习样本的获取是在三轴位置转台上进行的。通过三轴位置转台绕其三个转 动轴的转动来模拟钻井测斜仪的倾斜角、工具面角和方位角。测试点选取规则 为:
倾斜角测试点的选取规则是在0~75度范围内选取,在接近于零时取样点密, 角度接近75度时时测试取样点间隔增大,如在10度以内取的间隔小,大于10 度以10度的间隔选取;方位角的测试点的选取原则是在0~360度范围内均匀 选取,取样间隔根据精度和实验时间的要求选取;工具面角的测试点的选取原 则是在0~360度范围内均匀选取,取样间隔根据精度和实验时间的要求选取。
(3)进行神经网络训练,得到最优的神经网络模型参数。
采用步骤(1)得到的神经网络模型和步骤(2)得到的学习样本训练神经 网络并得到最优权值。
(4)根据倾斜角、初步计算的方位角、工具面角来计算真实方位角。
本发明的原理是:本发明针对采用陀螺罗经法的钻井测斜仪,惯性测量部 分由一个双轴动力调谐速率陀螺、两个石英加速度计和相应的
电子线路组成。 其中陀螺用来敏感地球
角速度分量、加速度用来敏感重力分量。通过陀螺和加 速度的输出值来计算钻井测斜仪也就是井管的方位角、倾斜角和工具面角。在 表示井管的方位角、倾斜角和工具面角时采用地理
坐标系XYZ(东北天)和测 斜仪坐标系xyz,其中x轴、y轴是两轴陀螺和两轴加速度计的敏感轴。测斜仪 坐标系xyz同测量井管的
姿态是一致的,因此xyz坐标系相对XYZ坐标系的姿 态便是井管的姿态,由此得到井管的轨迹。xyz坐标系相对XYZ坐标系的关系 如图3所示,XYZ是东北天地理坐标系坐标轴,xyz测斜仪坐标系坐标轴。x1y1z1 和x2y2z2是坐标系转换过程中的坐标系。图3中A为方位角、I为倾斜角、T为 工具面角,A&、I&、T&分别为相应的角速度矢量。坐标系旋转的顺序为:先按Z 轴顺
时针旋转到坐标系x1y1z1,,旋转角为方位角A,然后按坐标系x1y1z1的y1 轴旋转倾斜角I到坐标系x2y2z2,最后按x2y2z2坐标系的z2轴旋转测斜仪到坐标 系xyz
地球自转角速度和
重力加速度在东北天坐标系下的投影分量为:
ω=[0 Ωcosφ Ωsinφ]T (1)
a=[0 0 g]T (2)
式中:
Ω是地球自转角速度,φ是测量点的纬度,g是测量点的重力加速度。
经过坐标系旋转后,在测斜仪坐标系的角速度和加速度分量为:
由(1)、(2)、(3)、(4)式可算得工具面角T、倾斜角I、方位角A为:
以目前加速度的精度,由其计算的工具面角T、倾斜角I的误差比较小。方 位角的计算式中包括了工具面角T、倾斜角I和在x轴y轴上的角速度分量,方 位角的计算误差是它们共同作用的结果。当陀螺有误差和系统有安装误差时, 不同姿态下引起测斜仪方位角计算误差不同,且有很好的重复性。因此,事先 通过试验找到其误差的模型,并对测斜仪的方位角进行补偿,可以提高方位角 的精度。从式(7)可见当倾斜角I接近于90度附近时系统计算误差会比较大, 因此在采用陀螺罗盘方式时对倾斜角I进行限制,一般取倾斜角I小于75度。
本发明与现有技术相比的优点在于:
(1)本发明克服了建表法中把倾斜角、初步计算的方位角、工具面角做为 输入变量时插值算法精度下降的缺点,实现了钻井测斜仪倾斜角、初步计算的 方位角、工具面角到方位角的映射,达到对计算方位角误差补偿的效果。
(2)倾斜角测试点的分段选择方法充分反映了钻井测斜仪方位角计算误 差的特点,形成的学习样本更有针对性,即加快了神经网络的训练速度,又提 高了其计算精度。在达到相同精度的情况下,应用神经网络方法可比查表法采 取更少的测试点数,从而提高了效率。
附图说明
图1为本发明的神经网络的示意图;
图2为本发明获取神经网络学习样本的
流程图;
图3为本发明的测斜仪坐标系与地理系间的位置关系图。
具体实施方式:
针对以上提出的技术方案,采取以下几个步骤实现本发明:
(1)建立神经网络模型
本发明取钻井测斜仪的倾斜角、工具面角和初步计算的方位角作为神经网 络的输入变量,真实方位角作为神经网络的输出变量,选取四层(即输入层、 隐层1、隐层2、输出层)结构的前向神经网络。神经网络输入层节点数为3, 输出层节点数为1。考虑到钻井测斜仪的倾斜与方位对方位角的影响比较复杂, 隐层1节点数选为8,隐层2节点数选为8,神经网络结构如图1所示,其中 X1、X2、X3是神经网络的输入变量分别对应钻井测斜仪输出的倾斜角、工具 面角和初步计算的方位角,Y对应校正后的方位角,是神经网络的输出变量。
(2)确定钻井测斜仪的倾斜角和方位角的测试点
钻井测斜仪测试点的选取即要
覆盖钻井测斜仪的测量范围又要反映出钻井 测斜仪输出数据的特点。倾斜角的测试点在0~75度范围内取样,取样点为3 度、5度、10度、15度、25度、35度、45度、55度、65度、75度。方位 角的测试点在0~360度范围内均匀取样,取样间隔视精度要求而定。由于神经 网络具有泛化能力,取样间隔可比查表法放宽,本发明中取样间隔为30度。 工具面角的测试点在0~360度范围内均匀取样,取样间隔为30度。
(3)获取神经网络学习样本
在(2)所确定测试点处采集钻井测斜仪输出的倾斜角、初步计算方位角和 工具面角作为学习样本的输入值,三轴位置转台模拟的方位角作为学习样本的 期望输出值。具体操作是先把倾斜角、方位角、工具面角转到零位,然后顺序 转动倾斜角、方位角、工具面角,转遍全部测试点,同时采集测斜仪输出数据, 形成覆盖整个测量范围的神经网络的学习样本。学习样本获取流程如图2所示, 测试点个数Ni(i=1,2,3)分别表示倾斜角、方位角、工具面角的测试点个数, 可根据系统精度要求确定。
(4)进行神经网络训练,得到最优的神经网络模型参数
在前面(1)、(2)、(3)步确定的神经网络模型和学习样本的基础上,采 用BP算法对网络进行训练。在训练过程中,首先给出一组模型参数,以此参 数计算神经网络的输出,再对神经网络的输出和真实值进行比较得到计算误差。 然后根据误差,按BP算法
修改模型参数,使网络朝向能正确响应的方向不断 发展,直到网络的输出与期望的输出之差在允许的范围之内。此时的模型参数 便是最优的网络模型参数,也就是神经元间的连接权值及神经元的
阈值。权值 修改过程如下。
定义误差函数为
其中tp是神经网络的期望输出值,yp是神 经网络计算值。按
去修改神经元间连接权值,最终达到ep最小,确 定此时参数为最优参数。
(5)根据倾斜角、工具面角和初步计算的方位角来计算真实方位角 由
式可见,方位角的计算式中包括了工 具面角T、倾斜角I和在x轴y轴上的角速度分量,方位角的计算误差是它们共 同作用的结果。当陀螺有误差和系统有安装误差时,不同姿态下引起测斜仪方 位角计算误差不同,且有很好的重复性。
在第(4)步中得到了最优参数,也就是建立了初步计算出的工具面角、倾 斜角、方位角与真实方位角间的非线性映射关系,即向神经网络输入初步计算 出的工具面角、倾斜角、方位角,神经网络便会输出真实方位角。
在计算中,前一层神经元的输出作为下一层神经元的输入。设神经网络输 入信号:初步计算出的工具面角、倾斜角、方位角分别为x1,x2,x3,
输出信号为y。 下面是由x1,x2,x3求出y的过程。
第一隐层神经元输入为:
上式中n是输入神经元的个数,在此n=3,m1是输入神经元的个数,在此 m1=8,w1i,j为第一隐层神经元i与输入神经元间的连接权值,θ1j为第一隐层神 经元的阈值。
第一隐层神经元输出为:
O1i=f(I1i) (9)
第二隐层中共有m2个神经元,各神经元的输入是:
上式中m2=8,w2i,j为第二隐层神经元与第一隐层神经元间的连接权值,θ2j 为第二隐层神经元的阈值。
第二隐层神经元输出为:
O2i=f(I2i) (11)
取输出神经元的阈值为零,取线性函数为输出层神经元的激发函数,则输 出神经元的输出(也就是整个网络的输出)为:
上式中vi为输出神经元与第二隐层神经元的连接权值,由式(12)求出的y 便是补偿后真实的方位角。