CN102278108B - 小口径定向陀螺测斜仪连续测量模式标定方法 - Google Patents
小口径定向陀螺测斜仪连续测量模式标定方法 Download PDFInfo
- Publication number
- CN102278108B CN102278108B CN201110124400.XA CN201110124400A CN102278108B CN 102278108 B CN102278108 B CN 102278108B CN 201110124400 A CN201110124400 A CN 201110124400A CN 102278108 B CN102278108 B CN 102278108B
- Authority
- CN
- China
- Prior art keywords
- omega
- alpha
- axle
- gyro
- speed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Gyroscopes (AREA)
Abstract
本发明提供一种小口径定向陀螺测斜仪连续测量模式标定方法,包括如步骤:步骤1)建立加表输入输出关系数学模型;步骤2)建立陀螺输入输出关系数学模型;步骤3)采用三轴转台通过速率试验及六位置试验对测斜仪进行标定;步骤4)利用步骤1)和2)的数学模型采用计算机辅助程序,对步骤3)所测得36组速率数据和6组位置数据进行处理,计算输出测斜仪陀螺标度因数、安装误差和静态漂移系数。采用本发明标定方法能快速有效的对陀螺测斜仪进行标定,模型估计准确有效,操作简单,能在较短时间内完成所有参数标定,可靠性高。
Description
技术领域
本发明涉及钻井用连续测量陀螺测斜仪的标定方法,属于机电控制技术领域。
背景技术
定向钻井是油气开发的重大技术之一,在国外一些发达国家已经得到应用。定向钻井所用的测斜定向仪,有各种的磁仪表和框架式陀螺仪表。这两种常规测量技术经过数十年的和改进取得了一定的效果,但对小口径钻井测量已经不满足要求。尤其在我国仍大量使用的磁仪表(如磁通门)由于本身受地质和周围铁磁物质环境的影响,测量精度很低,甚至不能使用,一般的框架式陀螺仪表不仅精度低,而且直径大,使用前后的对准,校准也相当繁琐。而钻井用连续测量陀螺测斜仪进行标定是保证其工程质量的关键。现有技术中,对钻井用测斜仪的标定暂时无成型的标定方法,通常凭经验进行标定存在费时且准确度难以确定等问题。
发明内容
针对现有技术存在的上述不足,本发明的目的是提供一种能高效、快速和准确对钻探用小口径定向陀螺测斜仪进行标定的方法。
实现上述目的,本发明采用的技术手段是:一种小口径定向陀螺测斜仪连续测量模式标定方法,采用三轴转台通过六位置试验进行静态标定,包括如下步骤:
1)建立加表静态输入输出关系数学模型:
针对钻探用小口径定向陀螺测斜仪三个石英加速度计,设oxbybzb为惯性组合基体坐标系,与载体坐标系平行;建立石英加速度计输入输出关系数学模型如下:
式中,Fx、Fy、Fz为xb、yb、zb轴的加表输出,axb、ayb、azb为沿IMU坐标系三个轴向的加速度。kx、ky、kz为xb、yb、zb轴加表的标度因数,ax0、ay0、az0为xb、yb、zb轴加表的偏置,k2x、k2y、k2z为xb、yb、zb轴加表的二阶非线性系数(其值均为小量),αxy、αxz、αyx、αyz、αzx、αzy为加表安装误差(其值均为小量);
针对目前项目中陀螺测斜仪上的加表,忽略上式中的小量影响,建立如下简化模型:
各加表的标度因数、偏置、安装误差通过多位置试验进行标定;
2)建立陀螺静态输入输出关系数学模型:
针对钻探用小口径定向陀螺测斜仪的三轴挠性陀螺,设oxbybzb为惯性组合基体坐标系,与载体坐标系平行,陀螺输入输出关系数学模型如下:
式中,ux、uy、uz为IMU坐标系xb、yb、zb轴陀螺的电压输出,kgx、kgy、kgz为各轴陀螺的标度因数,αxy、αxz、αyx、αyz、αzx、αzy为陀螺安装误差,其值均为小量,但考虑到陀螺精度没有加表的高,故此几项系数不能忽略,ωxb、ωyb、ωzb为沿IMU坐标系xb、yb、zb的输入角速率,εx、εy、εz为各轴陀螺静态漂移,它们为:
式中,gxb、gyb、gzb为重力加速度g在IMU坐标系三个轴上的分量;
3)速率试验
试验前,将捷联惯性组合使速率转台转轴与地垂线平行安装方式之一安装在速率台上按表1所示,并使被标定轴与速率台轴重合;按预先设计好的速率测试点使速率台轴以其中某一角速率恒速转动;在速率转台转过某一基准角位置上开始记录陀螺仪输出,而且在一周内若干相邻等角度对称点上记录陀螺仪输出,直到速率转台转完一周;然后,以此角速率作反转试验,同样记录下这些数据;
表1.a速率试验中IMU在速率台上的安装方式
表中,OXYZ为速率台台体坐标系,oxbybzb为IMU坐标系,ω为转台对地速率,ωe为地球自转角速率,φ为当地地理纬度;
利用所得的试验数据可以标定IMU三轴陀螺的标度因数和安装误差;
a)标度因数和安装误差计算:
从各速率试验时陀螺输出数据中剔除转台静止时陀螺仪的输出,然后进行最小二乘法拟合。下面直接给出在速率试验数据,由式(2-1)、(2-2)所得出的各轴陀螺输出值的均值方程,见表2。表中方程用以求各轴陀螺标度因数和安装误差。
表2
表中,ωxb=ω(就是转台教速率角度),ωyb=ω+ωN(教速率角度速率和地球的分量),ωzb=ω+ωu(教速率角度速率和地球的分量)。轴(j=x、y、z)速率试验时,在试验时间内i轴(i=x、y、z)陀螺输出值的均值。
以下以xb轴速率试验为例,介绍由表2中方程计算陀螺标定因数和安装误差的方法。
b)标度因数计算
对xb轴n次速率试验,考虑x轴陀螺,由表2可得如下测试方程:
上式可表示如下:
y=Xb+v
式中,v为测试随机误差向量,
由最小二乘法可以求得
这样就求得了标定因数kgx。
c)安装误差计算
对xb轴速率大小为ω的正、反转试验,考虑y、z轴陀螺,由表2可得如下测试方程:
从上边两组方程可以求得两个安装误差角
用类似的方法对yb、zb轴速率试验时的各轴陀螺输出均值方程进行运算,可以求得另外2个标度因数和四个安装误差角,其中安装误差如下,
d)位置试验
对式(2-2)中的静态漂移系数,通过多位置试验进行标定。
由式(2-1)、(2-2)可以导出
在地球上用力矩反馈法进行静态漂移测试时,上式成为
式中,ωex、ωey、ωez一地速在IMU坐标系xb、yb、zb轴上的分量。
考虑x轴陀螺,n个位置试验可得到如下方程:
对y、z轴陀螺也可得到类似的方程,考虑到测试随机误差,它们可以统一写成如下的形式:
y=Xb+v
其中,y为观测值(n维),X为观测矩阵(n×4维),b为系数向量(n维),v为测试随机误差向量。
由最小二乘法,
若b是m维矢量,则X必须包括m维线性独立矢量,否则此法不可应用。
此时可对每一个位置试验,均可列出一组方程,然后对所有的位置试验列出的方程直接联立求解各个静态漂移系数;
3)采用三轴转台通过六位置试验对测斜仪进行标定:
第一步.将钻探用高精度陀螺测斜仪固定通过工装夹具将其安装在转台上,用水泡或地平仪对测斜仪在三轴速率转台上进行调平,尽量让其中两轴(如X,Y)保持水平,第三轴(如Z)垂直于地表面;
第二步.将Z轴朝天放置,X、Y轴水平。转台绕Z轴逆时针旋转,速率分别为:1°/s,5°/s,10°/s,15°/s,20°/s,25°/s,记录下各速率下陀螺输出的电压数据,加表输出脉冲数,保存为文档;采样时间为5分钟;
第三步.转台绕Z轴顺时针旋转。速率分别为:1°/s,5°/s,10°/s,15°/s,20°/s,25°/s,记录下各速率下陀螺输出数据,保存为文档;采样时间为5分钟;
第四步.将X、Y轴分别朝天,重复第二步,第三步动作,也分别以速率1°/s,5°/s,10°/s,15°/s,20°/s,25°/s正反转,并记录各种状态下陀螺输出数据,保存文档;采样时间为5分钟;
第五步.将陀螺测斜仪X、Y、Z轴分别按如下六位置放置:东、北、天;东、地、北;地、西、北;西、天、北;天、东、北;北、东、地;保存陀螺数据为文本文档,采样时间为5分钟;
六位置测试时,从一个位置到另一个位置,需要转位完成15秒后再保存数据,以保证陀螺数据已稳定到当前状态,测试数据才真实有效,不会影响标定结果;
4)利用步骤1)和2)的数学模型采用计算机辅助程序,对步骤3)所测得36组速率数据和6组位置数据进行处理,计算输出测斜仪陀螺标度因数、安装误差和静态漂移系数。
相比现有技术,本发明具有如下有益效果:(请尽量补充,关系到是否授权的问题!)
本发明设计了一套成型的小口径定向陀螺测斜仪连续测量模式标定方法,能高效快速,准确,非常简便的实现标定全过程;具有可靠性、实用性强,准确性和效率高的特点。
附图说明
图1是三轴转台及测斜仪安装示意图。
图2是钻探用高精度陀螺测斜仪标定时的六位置测试图。
图3是测试流程图(速率和位置试验无先后顺序)。
图4是钻探用高精度陀螺测斜仪的陀螺仪惯性短节结构示意图。
具体实施方式
本发明提供一种小口径定向陀螺测斜仪连续测量模式标定方法,包括如步骤:步骤1)建立加表输入输出关系数学模型;步骤2)建立陀螺输入输出关系数学模型;步骤3)采用三轴转台通过速率试验及六位置试验对测斜仪进行标定;步骤4)利用步骤1)和2)的数学模型采用计算机辅助程序,对步骤3)所测得36组速率数据和6组位置数据进行处理,计算输出测斜仪陀螺标度因数、安装误差和静态漂移系数。采用本发明标定方法能快速有效的对陀螺测斜仪进行标定,模型估计准确有效,操作简单,能在较短时间内完成所有参数标定,可靠性高。
下面结合实施例对本发明作进一步详细说明。
一种钻探用小口径定向陀螺测斜仪,主要由井上控制装置和井下测量设备两部分构成。井上控制装置由一台计算机和控制单元组成,通过电缆连接可实时跟踪显示井下测量设备的反馈信息;井下测量设备由马笼头、扶正器、采集遥测短节、陀螺仪惯性短节、抗压外壳、抗压保温瓶、加重短接、引鞋定向键和定向接头等组成。主要的传感器设备有二个双自由度的小型挠性陀螺和三个石英加速度计。该钻探用高精度陀螺测斜仪可以实现测斜仪的连续测量,大大简化了测试过程,数据自动读取和保存,不需要太多的人为手工操作,计算过程完全由计算机完成。
参见图4,陀螺仪惯性短节(惯性测量组)作为测斜仪的核心部件,由陀螺转子、挠性接头、力矩器、传感器、磁滞电机和密封壳体等元件组成。在陀螺底座上装有一对滚珠轴承,以支撑驱动电机轴的高速旋转,驱动轴一端装着电机的转子,另一端通过挠性接头装着陀螺转子,力矩器线圈骨架和传感器组件用螺钉与底座相连。陀螺底座两端的上、下端盖用电子束焊接并使陀螺密封。所述挠性陀螺的口径小于50mm。
其中,陀螺仪惯性短节包括转动电机47,以及减速器、转动轴和传感器(图中省略减速器、转动轴和传感器),传感器由二个双自由度的小型挠性陀螺41、42和三个互为正交的石英加速度计43、44、45组成,二个双自由度的第一小型挠性陀螺41和第二小型挠性陀螺42互为正交由旋转骨架连接,三个石英加速度计43、44、45互为正交由旋转骨架连接。本发明小口径定向陀螺测斜仪连续测量模式标定方法,适用于该钻探用小口径定向陀螺测斜仪,具体步骤包括:
一.加表标定
1、加表输入输出关系数学模型
针对钻探用小口径定向陀螺测斜仪中的三个石英加速度计,设oxbybzb为惯性组合基体坐标系,与载体坐标系平行。建立石英加速度计输入输出关系数学模型如下:
式中,Fx、Fy、Fz为xb、yb、zb轴的加表输出,axb、ayb、azb为沿IMU坐标系三个轴向的加速度。kx、ky、kz为xb、yb、zb轴加表的标度因数,ax0、ay0、az0为xb、yb、zb轴加表的偏置,k2x、k2y、k2z为xb、yb、zb轴加表的二阶非线性系数(其值均为小量),αxy、αxz、αyx、αyz、αzx、αzy为加表安装误差(其值均为小量)。
针对目前项目中IMU上的加表,建立如下简化模型。
各加表的标度因数、偏置、安装误差可以通过多位置试验进行标定。
2多位置试验
采用三轴转台通过六位置试验进行标定。
2-1试验方法
三轴转台如图1所示。OX、OY、OZ分别为转台的中环轴、内环轴、外环轴。oxyz为内框坐标系,初始时ox、oy轴水平,oz轴垂直朝天。把IMU组件以其坐标系oxbybzb分别与内框坐标系重合的方式安装在转台内框上。
选取如表1所示的六个位置进行试验。在每一个位置,采样数分钟,记录下采样时间和三轴加速度计输出的脉冲数。当某一个位置采数结束后,转动试验台的某一环轴到达另一位置继续上述试验,共进行六个位置的试验。利用这些数据就可以标定三轴加速度计的偏置、标度因数和安装误差。
表1六位置试验各位置及重力加速度在各轴上的分量
表中,H-水平,U-垂直向上,D-垂直向下。
2-2计算
对位置1和位置2试验,将axb=g、ayb=azb=0和axb=-g、ayb=azb=0分别代入(1-2),可得如下两组方程:
由式(1-3.a)、(1-3.b)可以求得
对位置3和4、5和6试验,作类似的处理,可求得其他8个参数。下面给出通过运算得到的结果。
以上是加表的标定。
二.陀螺标定
1、陀螺输入输出关系数学模型
针对钻探用小口径定向陀螺测斜仪中的三轴挠性陀螺,设oxbybzb为惯性组合基体坐标系,与载体坐标系平行。建立陀螺输入输出关系数学模型如下:
式中,ux、uy、uz为IMU坐标系xb、yb、zb轴陀螺的输出,kgx、kgy、kgz为各轴陀螺的标度因数,αxy、αxz、αyx、αyz、αzx、αzy为陀螺安装误差(其值均为小量),ωxb、ωyb、ωzb为沿IMU坐标系xb、yb、zb的输入角速率,εx、εy、εz为各轴陀螺静态漂移,它们为:
式中,gxb、gyb、gzb为重力加速度g在IMU坐标系三个轴上的分量。
标度因数和安装误差可以通过速率试验标定,静态漂移模型系数可以通过多位置试验标定。
2、速率试验
2-1试验方法
为了标定捷联惯性组合中陀螺的标度因数和安装误差,需要进行速率试验。
试验前,将捷联惯性组合按表1所示的安装方式之一安装在速率台上(即,使速率转台转轴与地垂线平行),并使被标定轴与速率台轴重合。按预先设计好的速率测试点(一般不少于11点)使速率台轴以其中某一角速率恒速转动。在速率转台转过某一基准角位置上开始记录陀螺仪输出,而且在一周内若干相邻等角度对称点上记录陀螺仪输出,直到速率转台转完一周(高速时可转若干周)。然后以此角速率作反转试验,同样记录下这些数据。
表1.a速率试验中IMU在速率台上的安装方式
表中,OXYZ为速率台台体坐标系,oxbybzb为IMU坐标系,ω为转台对地速率,ωe为地球自转角速率,φ为当地地理纬度。
对所有速率点重复上述过程。利用所得的试验数据可以标定IMU三轴陀螺的标度因数和安装误差。
2-2标度因数和安装误差计算
要从各速率试验时陀螺输出数据中剔除转台静止时陀螺仪的输出,然后进行最小二乘法拟合。下面直接给出在速率试验数据,由式(2-1)、(2-2)所得出的各轴陀螺输出值的均值方程,见表2。表中方程用以求各轴陀螺标度因数和安装误差。
表2
表中,ωxb=ω(就是转台教速率角度),ωyb=ω+ωN(教速率角度速率和地球的分量),ωab=ω+ωu(教速率角度速率和地球的分量)。轴(j=x、y、z)速率试验时,在试验时间内i轴(i=x、y、z)陀螺输出值的均值。
以下以xb轴速率试验为例,介绍由表2中方程计算陀螺标定因数和安装误差的方法。
2-2-1标度因数计算
对xb轴n次速率试验,考虑x轴陀螺,由表2可得如下测试方程:
上式可表示如下:
y=Xb+v
式中,v为测试随机误差向量,
由最小二乘法可以求得
这样就求得了标定因数kgx。
2-2-2安装误差计算
对xb轴速率大小为ω的正、反转试验,考虑y、z轴陀螺,由表2可得如下测试方程:
从上边两组方程可以求得两个安装误差角
用类似的方法对yb、zb轴速率试验时的各轴陀螺输出均值方程进行运算,可以求得另外2个标度因数和四个安装误差角,其中安装误差如下,
3位置试验
对式(2-2)中的静态漂移系数,通过多位置试验进行标定。
由式(2-1)、(2-2)可以导出
在地球上用力矩反馈法进行静态漂移测试时,上式成为
式中,ωex、ωey、ωez一地速在IMU坐标系xb、yb、zb轴上的分量。
考虑x轴陀螺,n个位置试验可得到如下方程:
对y、z轴陀螺也可得到类似的方程,考虑到测试随机误差,它们可以统一写成如下的形式:
y=Xb+v
其中,y为观测值(n维),X为观测矩阵(n×4维),b为系数向量(n维),v为测试随机误差向量。
由最小二乘法,
若b是m维矢量,则X必须包括m维线性独立矢量,否则此法不可应用。
此时可对每一个位置试验,均可列出一组方程,然后对所有的位置试验列出的方程直接联立求解各个静态漂移系数。
采用本发明方法对所述钻探用小口径定向陀螺测斜仪进行标定,具体步骤包括:
第一步.将钻探用小口径定向陀螺测斜仪固定在三轴速率转台上,尽量让其中两轴(如X,Y)保持水平,第三轴(如Z)垂直于地表面。如图1所示。
第二步.将Z轴朝天放置,X、Y轴水平。转台绕Z轴逆时针旋转,速率分别为:1°/s,5°/s,10°/s,15°/s,20°/s,25°/s,记录下各速率下陀螺输出的电压数据,加表输出脉冲数,保存为文档。采样时间为5分钟。
第三步.转台绕Z轴顺时针旋转。速率分别为:1°/s,5°/s,10°/s,15°/s,20°/s,25°/s,记录下各速率下陀螺输出数据,保存为文档。采样时间为5分钟。
第四步.将X、Y轴分别朝天,重复第二步,第三步动作,也分别以速率1°/s,5°/s,10°/s,15°/s,20°/s,25°/s正反转,并记录各种状态下陀螺输出数据,保存文档。采样时间为5分钟。
第五步.将陀螺测斜仪X、Y、Z轴分别按如下六位置放置:东、北、天;东、地、北;地、西、北;西、天、北;天、东、北;北、东、地;保存陀螺数据为文本文档,采样时间为5分钟。
六位置测试时,从一个位置到另一个位置,需要转位完成15秒后再保存数据,以保证陀螺数据已稳定到当前状态,测试数据才真实有效,不会影响标定结果。
第六步.调用编写好的计算机辅助程序,利用所测得36组速率数据和6组位置数据。计算输出测斜仪陀螺标度因数、安装误差和静态漂移系数。
采用本发明标定方法能快速有效的对陀螺测斜仪进行标定,模型估计准确有效,操作简单,能在较短时间内完成所有参数标定,可靠性高。
Claims (1)
1.小口径定向陀螺测斜仪连续测量模式标定方法,采用三轴转台通过六位置试验进行静态标定,包括如下步骤:
1)建立加表静态输入输出关系数学模型:
针对钻探用小口径定向陀螺测斜仪三个石英加速度计,设oxbybzb为IMU坐标系,与载体坐标系平行;建立石英加速度计输入输出关系数学模型如下:
式中,Fx、Fy、Fz为xb、yb、zb轴的加表输出,axb、ayb、azb为沿IMU坐标系三个轴向的加速度;kx、ky、kz为xb、yb、zb轴加表的标度因数,ax0、ay0、az0为xb、yb、zb轴加表的偏置,k2x、k2y、k2z为xb、yb、zb轴加表的二阶非线性系数,其值均为小量;αxy、αxz、αyx、αyz、αzx、αzy为加表安装误差,其值均为小量;针对目前项目中陀螺测斜仪上的加表,忽略上式中的小量影响,建立如下简化模型:
各加表的标度因数、偏置、安装误差通过多位置试验进行标定;
2)建立陀螺静态输入输出关系数学模型:
针对钻探用小口径定向陀螺测斜仪的三轴挠性陀螺,设oxbybzb为IMU坐标系,与载体坐标系平行,陀螺输入输出关系数学模型如下:
式中,ux、uy、uz为IMU坐标系xb、yb、zb轴陀螺的电压输出,kgx、kgy、kgz为各轴陀螺的标度因数,αxy、αxz、αyx、αyz、αzx、αzy为陀螺安装误差,其值均为小量,但考虑到陀螺精度没有加表的高,故此几项系数不能忽略,ωxb、ωyb、ωzb为沿IMU坐标系xb、yb、zb的输入角速率,εx、εy、εz为各轴陀螺静态漂移,它们为:
式中,gxb、gyb、gzb为重力加速度g在IMU坐标系三个轴上的分量;
3)速率试验
试验前,将捷联惯性组合使速率转台转轴与地垂线平行安装方式之一安装在速率台上按表1所示,并使被标定轴与速率台轴重合;按预先设计好的速率测试点使速率台轴以其中某一角速率恒速转动;在速率转台转过某一基准角位置上开始记录陀螺仪输出,而且在一周内若干相邻等角度对称点上记录陀螺仪输出,直到速率转台转完一周;然后,以此角速率作反转试验,同样记录下这些数据;
表1速率试验中IMU在速率台上的安装方式
表中,OXYZ为速率台台体坐标系,oxbybzb为IMU坐标系,ω为转台对地速率,ωe为地球自转角速率,φ为当地地理纬度;
利用所得的试验数据可以标定IMU三轴陀螺的标度因数和安装误差;
a)标度因数和安装误差计算:
从各速率试验时陀螺输出数据中剔除转台静止时陀螺仪的输出,然后进行最小二乘法拟合;下面直接给出在速率试验数据,由式(2-1)、(2-2)所得出的各轴陀螺输出值的均值方程,见表2,表中方程用以求各轴陀螺标度因数和安装误差;
表2
表中,ωxb=ω,ωyb=ω+ωN,转台对地速率和地球自转角速率的yb轴上分量;
以下以xb轴速率试验为例,介绍由表2中方程计算陀螺标定因数和安装误差的方法;
b)标度因数计算
对xb轴n次速率试验,考虑x轴陀螺,由表2可得如下测试方程:
上式可表示如下:
y=Xb+v
式中,v为测试随机误差向量,
由最小二乘法可以求得
这样就求得了标定因数kgx;
c)安装误差计算
对xb轴速率大小为ω的正、反转试验,考虑y、z轴陀螺,由表2可得如下测试方程:
从上边两组方程可以求得两个安装误差角
用类似的方法对yb、zb轴速率试验时的各轴陀螺输出均值方程进行运算,可以求得另外2个标度因数和四个安装误差角,其中安装误差如下,
d)位置试验
对式(2-2)中的静态漂移系数,通过多位置试验进行标定;
由式(2-1)、(2-2)可以导出
在地球上用力矩反馈法进行静态漂移测试时,上式成为
式中,ωex、ωey、ωez-地速在IMU坐标系xb、yb、zb轴上的分量;
考虑x轴陀螺,n个位置试验可得到如下方程:
对y、z轴陀螺也可得到类似的方程,考虑到测试随机误差,它们可以统一写成如下的形式:
y=Xb+v
其中,y为观测值,n维;X为观测矩阵,nⅹ4维;b为系数向量,n维;v为测试随机误差向量;
由最小二乘法,
若b是m维矢量,则X必须包括m维线性独立矢量,否则此法不可应用;
此时可对每一个位置试验,均可列出一组方程,然后对所有的位置试验列出的方程直接联立求解各个静态漂移系数;
3)采用三轴转台通过六位置试验对测斜仪进行标定:
第一步.将钻探用高精度陀螺测斜仪固定通过工装夹具将其安装在转台上,用水泡或地平仪对测斜仪在三轴速率转台上进行调平,尽量让其中两轴保持水平,第三轴垂直于地表面;
第二步.将Z轴朝天放置,X、Y轴水平;转台绕Z轴逆时针旋转,速率分别为:1°/s,5°/s,10°/s,15°/s,20°/s,25°/s,记录下各速率下陀螺输出的电压数据,加表输出脉冲数,保存为文档;采样时间为5分钟;
第三步.转台绕Z轴顺时针旋转,速率分别为:1°/s,5°/s,10°/s,15°/s,20°/s,25°/s,记录下各速率下陀螺输出数据,保存为文档;采样时间为5分钟;
第四步.将X、Y轴分别朝天,重复第二步,第三步动作,也分别以速率1°/s,5°/s,10°/s,15°/s,20°/s,25°/s正反转,并记录各种状态下陀螺输出数据,保存文档;采样时间为5分钟;
第五步.将陀螺测斜仪X、Y、Z轴分别按如下六位置放置:东、北、天;东、地、北;地、西、北;西、天、北;天、东、北;北、东、地;保存陀螺数据为文本文档,采样时间为5分钟;
六位置测试时,从一个位置到另一个位置,需要转位完成15秒后再保存数据,以保证陀螺数据已稳定到当前状态,测试数据才真实有效,不会影响标定结果;
4)利用步骤1)和2)的数学模型采用计算机辅助程序,对步骤3)所测得36组速率数据和6组位置数据进行处理,计算输出测斜仪陀螺标度因数、安装误差和静态漂移系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110124400.XA CN102278108B (zh) | 2011-05-13 | 2011-05-13 | 小口径定向陀螺测斜仪连续测量模式标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110124400.XA CN102278108B (zh) | 2011-05-13 | 2011-05-13 | 小口径定向陀螺测斜仪连续测量模式标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102278108A CN102278108A (zh) | 2011-12-14 |
CN102278108B true CN102278108B (zh) | 2014-03-26 |
Family
ID=45103796
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110124400.XA Expired - Fee Related CN102278108B (zh) | 2011-05-13 | 2011-05-13 | 小口径定向陀螺测斜仪连续测量模式标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102278108B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102607600B (zh) * | 2012-03-30 | 2014-10-29 | 航天科工惯性技术有限公司 | 滑动式测斜仪简易标校方法 |
CN102889076A (zh) * | 2012-09-14 | 2013-01-23 | 西安思坦仪器股份有限公司 | 陀螺测斜仪标定方法 |
CN104034347B (zh) * | 2013-03-05 | 2018-07-20 | 上海新跃仪表厂 | 一种星用半球谐振陀螺组合指标体系测量方法 |
CN103411625B (zh) * | 2013-07-23 | 2018-08-17 | 西安思坦仪器股份有限公司 | 动力调谐陀螺测斜仪的标定补偿方法及标定补偿模型 |
CN104675385B (zh) * | 2014-09-12 | 2018-03-09 | 北京精密机电控制设备研究所 | 安装误差校正方法及*** |
CN105587310B (zh) * | 2014-11-14 | 2018-11-16 | 中国航空工业第六一八研究所 | 一种陀螺测斜仪内部电磁干扰误差的标定方法 |
CN106884645A (zh) * | 2015-12-16 | 2017-06-23 | 航天科工惯性技术有限公司 | 陀螺测斜仪的标定方法 |
CN110595504A (zh) * | 2019-09-09 | 2019-12-20 | 武汉元生创新科技有限公司 | 一种惯性测量单元的自动标定方法及自动标定*** |
CN112814658B (zh) * | 2021-03-26 | 2023-07-11 | 中国石油化工股份有限公司胜利油田分公司石油工程技术研究院 | 超短半径井柔性测斜***及测量方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2004786C1 (ru) * | 1990-07-10 | 1993-12-15 | нин Лев Николаевич Бел | Инклинометр |
RU2159331C1 (ru) * | 1999-10-05 | 2000-11-20 | Общество с ограниченной ответственностью предприятие "АРКОН" | Способ определения азимута и зенитного угла скважины и гироскопический инклинометр |
RU2269001C1 (ru) * | 2004-05-05 | 2006-01-27 | Федеральное государственное унитарное предприятие "Ижевский механический завод" | Способ измерения траектории скважины по азимуту и двухрежимный бесплатформенный гироскопический инклинометр для его осуществления |
CN1743643A (zh) * | 2005-10-09 | 2006-03-08 | 北京航空航天大学 | 一种用于光纤陀螺油井测斜仪的机械骨架 |
CN1888388A (zh) * | 2006-07-17 | 2007-01-03 | 北京航空航天大学 | 一种适用于全光纤数字测斜仪的初始对准方法 |
CN1888386A (zh) * | 2006-07-17 | 2007-01-03 | 北京航空航天大学 | 适用于全光纤数字测斜仪的捷联惯性组合测量控制装置 |
RU2364718C1 (ru) * | 2008-01-09 | 2009-08-20 | Открытое акционерное общество Научно-производственное предприятие "Научно-исследовательский и проектно-конструкторский институт геофизических исследований геологоразведочных скважин (ОАО НПП "ВНИИГИС") | Установка для градуировки и калибровки инклинометров |
-
2011
- 2011-05-13 CN CN201110124400.XA patent/CN102278108B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2004786C1 (ru) * | 1990-07-10 | 1993-12-15 | нин Лев Николаевич Бел | Инклинометр |
RU2159331C1 (ru) * | 1999-10-05 | 2000-11-20 | Общество с ограниченной ответственностью предприятие "АРКОН" | Способ определения азимута и зенитного угла скважины и гироскопический инклинометр |
RU2269001C1 (ru) * | 2004-05-05 | 2006-01-27 | Федеральное государственное унитарное предприятие "Ижевский механический завод" | Способ измерения траектории скважины по азимуту и двухрежимный бесплатформенный гироскопический инклинометр для его осуществления |
CN1743643A (zh) * | 2005-10-09 | 2006-03-08 | 北京航空航天大学 | 一种用于光纤陀螺油井测斜仪的机械骨架 |
CN1888388A (zh) * | 2006-07-17 | 2007-01-03 | 北京航空航天大学 | 一种适用于全光纤数字测斜仪的初始对准方法 |
CN1888386A (zh) * | 2006-07-17 | 2007-01-03 | 北京航空航天大学 | 适用于全光纤数字测斜仪的捷联惯性组合测量控制装置 |
RU2364718C1 (ru) * | 2008-01-09 | 2009-08-20 | Открытое акционерное общество Научно-производственное предприятие "Научно-исследовательский и проектно-конструкторский институт геофизических исследований геологоразведочных скважин (ОАО НПП "ВНИИГИС") | Установка для градуировки и калибровки инклинометров |
Also Published As
Publication number | Publication date |
---|---|
CN102278108A (zh) | 2011-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102278108B (zh) | 小口径定向陀螺测斜仪连续测量模式标定方法 | |
CN102140913B (zh) | 钻探用小口径定向陀螺测斜仪 | |
CN105180968B (zh) | 一种imu/磁强计安装失准角在线滤波标定方法 | |
CN101290326B (zh) | 石英挠性加速度计测量组件的参数辨识标定方法 | |
CN101514900B (zh) | 一种单轴旋转的捷联惯导***初始对准方法 | |
CN100565115C (zh) | 多位置捷联寻北***方位效应的标定方法 | |
CN1330935C (zh) | 安装误差角与标度因数解耦的微惯性测量单元精确标定方法 | |
CN101975872B (zh) | 石英挠性加速度计组件零位偏置的标定方法 | |
CN105806364B (zh) | 一种矿用回转钻机测斜仪探管的校准方法 | |
CN101493008A (zh) | 基于mems器件的捷联惯性导航陀螺测斜仪 | |
CN104165638B (zh) | 一种双轴旋转惯导***多位置自主标定方法 | |
CN110108300B (zh) | 一种基于卧式三轴转台的imu正六面体标定方法 | |
CN104235618A (zh) | 一种基于mems惯性测量单元的管道测绘及缺陷定位装置及其管道测绘及缺陷定位方法 | |
CN102564461A (zh) | 一种基于双轴转台的光学捷联惯导***的标定方法 | |
CN101629969B (zh) | 一种低精度光纤惯组输出误差的标定补偿和测试方法及装置 | |
CN103852085B (zh) | 一种基于最小二乘拟合的光纤捷联惯导***现场标定方法 | |
CN102889076A (zh) | 陀螺测斜仪标定方法 | |
CN103104251A (zh) | 一种用于小井斜下光纤陀螺测斜仪的方位角和工具面角的精度提高方法 | |
CN100487374C (zh) | 消除动调式陀螺测斜仪惯性器件常值误差的对转定位机构 | |
CN102589573A (zh) | 微型组合导航***中的传感器野外标定方法 | |
CN106017452A (zh) | 双陀螺抗扰动寻北方法 | |
CN106767925A (zh) | 带双轴转位机构的惯导***三位置参数辨识对准方法 | |
CN102174885B (zh) | 一种钻井用高精度点测陀螺测斜仪的标定方法 | |
CN102607594A (zh) | 捷联惯导光纤陀螺***级误差参数现场标定方法 | |
CN105134171B (zh) | 一种两轴光纤陀螺连续测斜***的实现方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140326 Termination date: 20210513 |
|
CF01 | Termination of patent right due to non-payment of annual fee |