基于惯性导航系统的海底管道轨迹描述技术研究
上传者:马继光|上传时间:2015-05-05|密次下载
基于惯性导航系统的海底管道轨迹描述技术研究
384 第十五届中国海洋(岸)工程学术讨论会论文集
基于惯性导航系统的海底管道轨迹描述技术研究
史永晋1,2
(1.中国石油大学(华东),山东 东营257061;2.胜利石油管理局 钻井工艺研究院,山东 东营 257017)
摘要:由于海底地质容易发生体貌变化,难以精确提供浅海海底管道的实际轨迹和实时地理信息。因此,无法确定管道内缺陷点的精确位置,导致维修设备难以对其进行精确的定位维修。主要利用INS(惯性导航系统)测录及姿态解算系统来完成海底管道轨迹描述和缺陷定位,通过分析管道航迹的经纬度定位得到管道内壁缺陷与大地坐标间的精确位置。该系统算法通过实验室试验进行了误差分析和校正,最终在海上油田输油管道内进行了工程应用,得到的数据轨迹与实际轨迹满足误差范围要求。计算分析和实验证明,采用惯性导航系统技术来解决缺陷定位和管道轨迹问题是准确和可行的。 关键词:缺陷定位;海底管道;轨迹分析;定位技术
随着近浅海石油开发规模的不断扩大和发展,尤其是油气混输技术的日渐成熟和应用,利用管道输送油、气,具有安全可靠、成本低、效率高的优势,因此海底管道在海洋油气开发中得到广泛应用[1]。定期对在役海底管道进行位貌检验和缺陷检测是确保管道安全、高效运行的可靠手段[2]。但是,我国浅海海底管道的敷设受装备、工艺、经验等因素的制约,敷设时难以精确测定记录其实际轨迹。另外,我国滩海交界地带海洋环境、浅层工程地质、海底地貌条件十分复杂,海床极不稳定,易遭受冲蚀,海底管道运行过程由于海底表层地基不稳定、海流冲淘及海上意外事故等原因,易发生位貌变化,难以建立详尽实时的地理信息系统。检测装置如果只记录自身周向坐标与行走距离,难以进行精确大地定位,不能为管线维修装置提供缺陷的精确位置[3]。因此国内外对海底管道内的定位技术十分重视,进行了大量的研究工作。其中,胜利油田钻井工艺研究院承担的国家“863”项目“海底管道内爬行器及其检测技术”和“海底管道内爬行器工程样机的研制”就是针对海底管道内壁缺陷检测定位的一项高新课题。海底管道内缺陷定位技术是课题的关键技术之一。我们经过多方调研技术比较,决定采用INS( Inertial Navigation System,惯性导航系统)技术来解决该问题。
1系统组成
在海底管道检测器中,主要用INS(惯性导航系统)测录及其姿态解算系统来完成海底管道内缺陷定位功能。它可实现海底管道爬行检测器的INS测录和离线INS数据处理,从而得到航迹的经纬度定位。
工程中多采用捷联惯性导航系统[4],捷联惯性导航系统是基于惯性测量元件和导航计算机构成的多功能自主式导航系统。捷联惯性导航系统通过惯性敏感元件(陀螺仪和加速度计)测量载体相对惯性空间的角速度和加速度,实时输出载体当前的角速度和加速度。捷联惯性导航系统具有不依赖外界信息,隐蔽性好等特点。但是由于捷联惯性导航系统中的陀螺的误差存在随着时间积累而逐渐增大的缺点,对于载体运行时间较长而且又需要高精度定位要求的系统来说,是无法克服的弱点。当前解决此类问题的方法多是采用组合导航的方式,即惯性导航系统+GPS导航系统。但是,海底管道的运行环境有其特殊性。海底管道是双层管道,两个管道间是聚氨脂泡沫保温层,有的管道外部还有混凝土保护层,并且深埋在2-3 m深的海底[5]。海水、海泥双层管道和聚氨脂泡沫保温层组成了复杂的运行环境,该环境具有较强的电磁屏蔽性,单纯的GPS定位的方法无法采用[6-8]。这是对本课题组最大的挑战。
INS测录及姿态解算系统,包括测量、记录、存储及其离线解算四个方面的工作。硬件整体系统结构如图1所示。
内容需要下载文档才能查看
图1 INS测录及姿态解算系统硬件结构图
INS测录及姿态解算系统安装在INS测录监视舱中,为数据处理系统计算航迹的经纬度提供载体的三个角速度( 横滚角/俯仰角/航向角)和X/Y/Z 三个方向的加速度和时间等。该系统由微型惯性测量装置、嵌入
第十五届中国海洋(岸)工程学术讨论会论文集 385
式计算机、存储模块和信息接口几部分组成。INS解算系统功能由分析计算机中的数据处理软件完成。它根据捷联惯性导航系统的姿态解算原理及误差修正的方法,进行载体运行轨迹的姿态解算,输出载体的姿态及经纬度的运行轨迹。
通过计算机编程可实现INS解算系统的友好人机界面。解算的处理过程中涉及大量的矩阵推导和积分求解过程,为此采用了后台调用强大的数据处理工具—Matlab的先进设计方案。在捷联惯性导航系统中,陀螺和加速度计直接固联在运载体上,直接测量运载体的角运动和线运动。
2 工作原理及算法分析
2.1 工作原理
惯性导航系统的工作原理是建立在牛顿经典力学(一个物体如果没有外力作用,将保持匀速直线运动)的基础上的。物体的加速度正比于作用在物体上的外力。如果能够测量得到加速度,那么通过加速度对时间的连续数学积分就可计算得到物体的速度和位置的变化。捷联式惯性导航系统是一种利用安装在运载体上的陀螺仪和加速度计来测定运载体位置的一个系统。通过陀螺仪和加速度计的测量数据,可以确定运载体在惯性参考坐标系中的运动,同时也能够计算出运载体在惯性参考坐标系中的位置。 2.2 算法分析
捷联惯导姿态的解算方法有欧拉角法、方向余弦法、四元素法, 而四元数法比其他方法明显优越就是在于它有计算量小、计算精度高、可避免奇异性,因此本项目运用的是四元素法[9-10]。捷联惯导姿态的算法中主要设计两个坐标系:地理坐标系Oxnynzn和载体坐标系Oxbybzb。地理坐标系又称导航坐标系,即常用的NED坐标系。载体坐标系与载体固连,其原点在载体的质心O,Oyb轴沿目标载体纵轴且指向前方,Oxb轴沿目标载体立轴指向上方,Ozb轴沿目标载体横轴指向右。载体坐标系相对于地理坐标系的姿态角由航向角、俯仰角、横滚角来表示。载体纵轴轴向在水平面的投影与北向基准线之间的夹角α称为航向角,北向偏东为正。载体纵轴与水平面之间的夹角β称为俯仰角;载体横轴与水平面之间的夹角θ称为横滚角[11]。
内容需要下载文档才能查看
图2 地理坐标系和载体坐标系示意
载体坐标系到地理坐标系的姿态矩T阵为[12]:
O
→xnynznOZ
α
O
n
xb1yb1zb1
→O
β
O
xb1
xb2yb2zb2→Oxbybzb
Oxb2
θ
(1) (2)
载体坐标系到地理坐标系的姿态矩阵T 为
xb T11 xn
= Tn yTy=nbb 21
T31 zb zn
T12T22
T32
T13 xb
y T23 b T33 zb
其中,
内容需要下载文档才能查看
四元素法用微分方程求解方向余弦矩阵,实现坐标系的转换。
Q=qo+q1i+q2j+q3k
(3) (4)
386 第十五届中国海洋(岸)工程学术讨论会论文集
式中: qj (j=0,1,2,3,)为实数,称它为四元素, i2+j2+k2= 1,ij=k,jk=i,ki=j. 四元素微分方程的标量形式写成矩阵形式为:
0-ω1-ω2
1 ω10ω3
Q=
2 ω2-ω30
ω3ω2-ω1
其中,ω1,ω2,ω3—载体的旋转角速度矢量。其中Q的微分方程的初值问题为:
Q=
1
Q W,Q(t2
-ω3 q0
q -ω2 1 ω1 q2
0 q3
)=Q0
(5)
(6)
其中t0为运动的初始时刻,Q0为初始对准获得的四元素。
在(6)求解过程中,通常采用四阶龙格库塔法进行计算。
1
Qj+1=Qj+(k1+k2+k3+k4)dt
6
式中,dt为计算周期,k1,k2,k3,k4是四阶龙格库塔系数,可表示为
K1=0.5Q(t) W(t)
(7)
K2=0.5(Q(t)+K1dt/2) W(t+dt/2)
K3=0.5(Q(t)+K2dt/2) W(t+dt/2)K4=0.5(Q(t)+K3dt/2) W(t+dt)
(8)
(9)
(10)
解算出Q后,可由下式求姿态矩阵
222 q02(q1q2 q0q3)2(q1q3+q0q2) +q12 q2 q3
222
T= 2(q1q2 q0q3)q02(q1q3+q0q2) q12+q2 q3
222 2(q1q3 q0q2)2(q2q3+q0q1)q0 q12 q2+q3
通过姿态矩阵T可将加速度计测量的比力f转换成为f
bp
f
p
=Tf
b
2.3 轨迹计算
得到转换后的加速度后,可进行速度的即时修正。于是,惯性导航的基本方程为
式中,ie
(11)
为地球自转角速度在当地地理坐标系中的表示,ep为地理坐标系相对于地球固定坐标系的转
Vep=(fp (2ie+ep))×Vep+g
动角速度,g为当地重力矢量,它由万有引力重力影响和由地球转动引起的向心加速度组成。
Vep=[VN
VE
VD]表示相对于地球在当地地理坐标系的速度。
T
VN
(纬度)
(R0+h)
VsecL
(经度) λ=E
(R0+h)L=
(12) (13)
2.4 误差分析
如果没有外界信息的标定和补偿,管线误差随时间积累,因此需要对误差进行分析[13-14]。 惯性导航系统的主要误差源,一般都考虑以下几个方面: 1)惯性仪表的安装误差和标度因子误差。
2)陀螺的漂移和加度计的零位误差。这部分包括常值漂移和随机两部分。 3)初始条件误差、包括导航参数和姿态航向的初始误差。
4)计算误差,主要考虑姿态航向系统的计算误差,也就是数学平台的计算误差。 5)飞行器的角运动所引起的动态误差。
对于我们采用IMU(惯性测量传感器)来进行导航解算,主要考虑的误差是1、2、3、5项。由于管道检测器在管道内行驶,速度很慢,近似为匀速运动,陀螺和加速度计引起的误差对爬行检测器的解算速度会产生很大的影响,由此也引起较大位置解算误差。
第十五届中国海洋(岸)工程学术讨论会论文集 387
下面根据误差公式,计算陀螺的随机游走,加速度计的误差,以及初始对准误差引起的系统误差。设取用陀螺漂移为0.05°/s,加速度计误差为1 mg,计算陀螺漂移引起的系统误差。
设三个方向的陀螺漂移分别为:εE,εN,εD,ωie为地球自转角速度,
ωs=g/R为舒勒频率, 为纬度。
2.4.1 计算陀螺的随机游走引起的系统误差
纬度误差:
ω11
δ (t)=2s2sinωiet sinωst)εE
ωs ωieωieωs
ωs2ωiesin 11sin
+ωtωtcoscos)2 εN sie222
ωωωωωsieie sie
ωs2cos ωiecos cos
coscosωωtt εDies2222
()ωωωωωωsieie iesie
λ(t)= tg ωie
(1 cosωiet)
ωietg
ωtωt(coscos) ies εE
ωs2 ωie2
(14)
经度误差:
sec ω ωs2tg sin (s2 ωie2(cos )2)
sinωst sinωiet tcos εN 2222
()()ωωωωωω iesieiesie ωs2sin ωie2sin
sinωiet sinωst tsin εD2222
ωie(ωs ωie) ωie(ωs ωie)
(15)
内容需要下载文档才能查看经度
内容需要下载文档才能查看差
纬度差
图3 陀螺引起的经度差和纬度差
图3分别为1个小时随机游走的经纬度差。可见经度差是随时间而积累的,纬度差是周期性变化的。
一个小时的经度差误差约为330 120 m。一个小时的经纬度差约为326 910 m。这样的精度不能达满足我们的需要。
2.4.2 加速度计的漂移引起的系统误差
E, N为东向和北向的加速度计漂移。引起的纬度和经度误差分别为:
δ (t)=
N
(1 cosωst)) g
(16)
Esec( )
(17) (1 cosωst)
g
图4为加速度计误差引起的经纬度误差。可见是周期性变化的。一个小时后,加速度计引起的系统误差:经度方向为160 m,纬度方向为1 215 m。
δλ(t)=
388
内容需要下载文档才能查看-3
第十五届中国海洋(岸)工程学术讨论会论文集
-3
经度差
内容需要下载文档才能查看纬度差
图4 加速度计引起的纬度误差和经度误差
2.4.3 初始条件的误差
首先考虑初始条件误差φN(0),φE(0),φD(0),δVN(0),δVE(0),δ (0)为三个姿态角初始误差,北向和东向的速度初始误差,纬度初始误差。设初始姿态角误差为[0.1 0.1 0.1];纬度差0.001度。
δVN(t)=δVN(0)cosωst+
g
(ωsinωst ωiesinωiet)φE(0)+
ωs ωies
(16)
gωiesin gωcos
cosωiet cosωst]φN(0)+2ie2[cosωst cosωiet]φD(0)22[ωs ωieωs ωie+
R
(ωie2ωssinωst ωs2ωiesinωiet)δ (0)2
ωs ωie
2
δ (t)=
δVNω
sinωst+2s2(cosωiet cosωst)φE(0)+Rωsωs ωie
2
(17)
ωs2ωiesin 11
(sinωiet sinωst)φN(0)+22
ωs ωieωieωsωscos
(ωiesinωst ωssinωiet)φD(0)+
ωs2 ωie2
2
1
(ωs2cosωiet ωie2cosωst)δ (0)2
ωs ωie
δV(0)sec sinωsttg δλ(t)=E+(ωs2sinωiet ωieωssinωst)φE(0)+ωs ωieRωs
sec (ωs2 ωie2(cos 2)) ωs2tg sin ωcostcosωiet cos φN(0)+ s2222
ωs ωieωs ωie ωie2sin ωie2sin
ωω costcostsin 2 φD(0)+ies2
ωs2 ωie2 ωs ωie
ωstg
(ωsinωiet ωiesinωst)δ (0)+δλ(0)
ωs ωies
经度
内容需要下载文档才能查看差
(18)
内容需要下载文档才能查看-3
-3
纬度差
秒秒
图5 初始误差引起的系统经纬度误差
由以上的分析可知陀螺的漂移误差影响最大。
3 实际测量及数据处理算例
3.1校正测试
将INS系统绕400 m操场一周,采用IMU YH5100,用DGPS所采集的数据作为真实轨迹。图6反映了所算出的速度已经发散;从图7可以看出,姿态角的解算还在可接受的范围;图8中所示的校正前解算出的位置比较。图9所示的是采用了我们新设计的解算方案。可见,虽然与DGPS轨迹还没有吻合,但是已基本达到了航迹计算的要求。
下载文档
热门试卷
- 2016年四川省内江市中考化学试卷
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
- 山东省滨州市三校2017届第一学期阶段测试初三英语试题
- 四川省成都七中2017届高三一诊模拟考试文科综合试卷
- 2017届普通高等学校招生全国统一考试模拟试题(附答案)
- 重庆市永川中学高2017级上期12月月考语文试题
- 江西宜春三中2017届高三第一学期第二次月考文科综合试题
- 内蒙古赤峰二中2017届高三上学期第三次月考英语试题
- 2017年六年级(上)数学期末考试卷
- 2017人教版小学英语三年级上期末笔试题
- 江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
- 重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
- 江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
- 江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
- 山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
- 【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
- 四川省简阳市阳安中学2016年12月高二月考英语试卷
- 四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
- 安徽省滁州中学2016—2017学年度第一学期12月月考高三英语试卷
- 山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
- 福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
- 甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷
网友关注
- 2018北京公务员考试行测每日一练判断推理练习题11.23
- 2018北京公务员考试行测每日一练判断推理练习题11.01
- 2018北京公务员考试行测每日一练判断推理练习题11.27
- 2018北京公务员考试行测每日一练判断推理练习题答案11.01
- 2017年北京公务员考试行测真题——判断推理题(八)
- 2017年北京公务员考试行测真题——判断推理题(五)
- 2016北京公务员考试行测真题——数量关系题(四)
- 行测题库:行测每日一练判断推理练习题12.04
- 2018北京公务员考试行测每日一练判断推理练习题答案11.10
- 2017年北京公务员考试行测真题——判断推理题(一)
- 2018北京公务员考试行测每日一练言语理解练习题11.22
- 公务员面试热点模拟题:电商逼死了实体店?
- 2018北京公务员考试:行测每日一练言语理解练习题10.11
- 2018北京公务员考试行测每日一练言语理解练习题答案11.13
- 2018北京公务员考试行测每日一练判断推理练习题11.06
- 2017年北京公务员考试行测真题——判断推理题(四)
- 2018北京公务员考试行测每日一练判断推理练习题答案11.23
- 2017年北京公务员考试行测真题——判断推理题(二)
- 2018北京公务员考试行测每日一练言语理解练习题11.13
- 2018北京公务员面试模拟题:把道德修养当做人生必修课
- 2018北京公务员考试行测每日一练资料分析练习题11.07
- 2018北京公务员考试行测每日一练判断推理练习题答案11.06
- 2018北京公务员考试行测每日一练数量关系练习题答案11.03
- 面试模拟题:“状元”的未来
- 2017年北京公务员考试行测真题——数量关系题(一)
- 2018北京公务员考试行测每日一练数量关系练习题答案11.28
- 2017年北京公务员考试行测真题——数量关系题(二)
- 2018北京公务员考试行测演练厅之生活常识模拟题
- 面试模拟题:环境保护与GDP的爱恨情仇
- 2018北京公务员考试行测每日一练资料分析练习题11.08
网友关注视频
- 第8课 对称剪纸_第一课时(二等奖)(沪书画版二年级上册)_T3784187
- 第五单元 民族艺术的瑰宝_16. 形形色色的民族乐器_第一课时(岭南版六年级上册)_T1406126
- 沪教版八年级下次数学练习册21.4(2)无理方程P19
- 苏教版二年级下册数学《认识东、南、西、北》
- 沪教版牛津小学英语(深圳用)五年级下册 Unit 1
- 19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T502436
- 苏科版数学七年级下册7.2《探索平行线的性质》
- 七年级英语下册 上海牛津版 Unit3
- 苏科版数学八年级下册9.2《中心对称和中心对称图形》
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 2
- 冀教版小学数学二年级下册第二周第2课时《我们的测量》宝丰街小学庞志荣
- 外研版英语七年级下册module3 unit2第二课时
- 第19课 我喜欢的鸟_第一课时(二等奖)(人美杨永善版二年级下册)_T644386
- 人教版二年级下册数学
- 冀教版小学数学二年级下册第二单元《租船问题》
- 二年级下册数学第一课
- 8.练习八_第一课时(特等奖)(苏教版三年级上册)_T142692
- 【部编】人教版语文七年级下册《过松源晨炊漆公店(其五)》优质课教学视频+PPT课件+教案,辽宁省
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 8
- 北师大版小学数学四年级下册第15课小数乘小数一
- 《空中课堂》二年级下册 数学第一单元第1课时
- 沪教版牛津小学英语(深圳用) 五年级下册 Unit 12
- 苏科版数学 八年级下册 第八章第二节 可能性的大小
- 河南省名校课堂七年级下册英语第一课(2020年2月10日)
- 每天日常投篮练习第一天森哥打卡上脚 Nike PG 2 如何调整运球跳投手感?
- 沪教版八年级下册数学练习册21.3(3)分式方程P17
- 沪教版八年级下册数学练习册21.4(1)无理方程P18
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
- 第12章 圆锥曲线_12.7 抛物线的标准方程_第一课时(特等奖)(沪教版高二下册)_T274713
精品推荐
- 2016-2017学年高一语文人教版必修一+模块学业水平检测试题(含答案)
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
分类导航
- 互联网
- 电脑基础知识
- 计算机软件及应用
- 计算机硬件及网络
- 计算机应用/办公自动化
- .NET
- 数据结构与算法
- Java
- SEO
- C/C++资料
- linux/Unix相关
- 手机开发
- UML理论/建模
- 并行计算/云计算
- 嵌入式开发
- windows相关
- 软件工程
- 管理信息系统
- 开发文档
- 图形图像
- 网络与通信
- 网络信息安全
- 电子支付
- Labview
- matlab
- 网络资源
- Python
- Delphi/Perl
- 评测
- Flash/Flex
- CSS/Script
- 计算机原理
- PHP资料
- 数据挖掘与模式识别
- Web服务
- 数据库
- Visual Basic
- 电子商务
- 服务器
- 搜索引擎优化
- 存储
- 架构
- 行业软件
- 人工智能
- 计算机辅助设计
- 多媒体
- 软件测试
- 计算机硬件与维护
- 网站策划/UE
- 网页设计/UI
- 网吧管理