发明名称 适用于小角度井斜角测量的陀螺测斜仪姿态测量解算方法
摘要 本发明涉及陀螺测斜仪姿态提取技术,具体是一种适用于小角度井斜角测量的陀螺测斜仪姿态测量解算方法。本发明解决了现有陀螺测斜仪在小井斜角测量时因姿态提取方法导致方位角和横滚角测量精度差的问题。适用于小角度井斜角测量的陀螺测斜仪姿态测量解算方法,该方法是采用如下步骤实现的:1)组成陀螺测斜仪;2)求取姿态变换矩阵;3)提取俯仰角;4)若俯仰角是小角度井斜角,解算出横滚角;5)若俯仰角是大角度井斜角,则按大井斜状态完成解算。本发明所述的适用于小角度井斜角测量的陀螺测斜仪姿态测量解算方法彻底解决了现有陀螺测斜仪在小井斜角测量时因姿态提取方法导致方位角和横滚角测量精度差的问题。
申请公布号 CN102536207B 申请公布日期 2014.03.26
申请号 CN201110454317.9 申请日期 2011.12.30
申请人 中北大学 发明人 张晓明;龙达峰;李杰;刘俊;崔星;刘喆;黄建林;尚剑宇;范玉宝;李永慧
分类号 E21B47/022(2012.01)I 主分类号 E21B47/022(2012.01)I
代理机构 太原科卫专利事务所(普通合伙) 14100 代理人 朱源
主权项 1.一种适用于小角度井斜角测量的陀螺测斜仪姿态测量解算方法,其特征在于:该方法是采用如下步骤实现的:1)采用三轴陀螺仪、三轴加速度计和相应的电子线路组成陀螺测斜仪,并确定陀螺测斜仪的载体坐标系和导航坐标系;2)选取陀螺测斜仪的工作方式,并根据陀螺测斜仪的工作方式选取相应的姿态变换矩阵求解方法,求取姿态变换矩阵;3)通过求取的姿态变换矩阵提取陀螺测斜仪的姿态所对应的俯仰角,根据俯仰角和井斜角之间的关系计算出井斜角,并判断井斜角是否为小角度井斜角;4)若井斜角是小角度井斜角,则利用安装在陀螺测斜仪纵轴方向的陀螺仪作为辅助测量解算出陀螺测斜仪的姿态所对应的横滚角,并利用陀螺测斜仪的姿态所对应的方位角和横滚角之间的关系提取方位角,以此完成陀螺测斜仪的三个姿态角的解算;若井斜角是大角度井斜角,则按大井斜状态完成陀螺测斜仪的三个姿态角的解算;所述步骤1)中,陀螺测斜仪采用北东地<img file="694454DEST_PATH_IMAGE001.GIF" wi="36" he="18" />地理坐标系为导航坐标系;所述步骤2)中,陀螺测斜仪的工作方式选取静态测量方式或随钻测量方式;当陀螺测斜仪的工作方式选取静态测量方式时,姿态变换矩阵求解方法采用kalman滤波初始对准方法;当陀螺测斜仪的工作方式选取随钻测量方式时,姿态变换矩阵求解方法采用kalman滤波初始对准方法和四子样等效旋转矢量算法;姿态变换矩阵为:<img file="549278DEST_PATH_IMAGE002.GIF" wi="517" he="84" />;上式中:<img file="329015DEST_PATH_IMAGE003.GIF" wi="14" he="20" />为俯仰角,<img file="661907DEST_PATH_IMAGE004.GIF" wi="17" he="18" />为方位角,<img file="413963DEST_PATH_IMAGE005.GIF" wi="14" he="18" />为横滚角;所述步骤3)中,俯仰角的提取公式如下:<img file="236425DEST_PATH_IMAGE006.GIF" wi="100" he="25" />;式中:<img file="441142DEST_PATH_IMAGE003.GIF" wi="14" he="20" />为俯仰角,<img file="577725DEST_PATH_IMAGE007.GIF" wi="25" he="25" />为<img file="246604DEST_PATH_IMAGE008.GIF" wi="38" he="23" />;根据俯仰角和井斜角的互余关系计算出井斜角;所述步骤4)中,若井斜角是小角度井斜角,横滚角的积分解算公式如下:<img file="177651DEST_PATH_IMAGE009.GIF" wi="125" he="36" />;上式中:<img file="869663DEST_PATH_IMAGE010.GIF" wi="14" he="33" />为横滚角的变化率,<img file="541428DEST_PATH_IMAGE011.GIF" wi="28" he="26" />为陀螺测斜仪载体坐标系相对于惯性系转动的角速率分量,<img file="64814DEST_PATH_IMAGE012.GIF" wi="22" he="25" />为地球自转速率,<img file="432341DEST_PATH_IMAGE013.GIF" wi="16" he="18" />为地球纬度,<img file="611650DEST_PATH_IMAGE011.GIF" wi="28" he="26" />为陀螺测斜仪载体坐标系相对于地理坐标系转动的角速度分量;积分计算以上一次大井斜角所解算的横滚角<img file="886773DEST_PATH_IMAGE005.GIF" wi="14" he="18" />作为积分初值;方位角的提取公式如下:<img file="467927DEST_PATH_IMAGE014.GIF" wi="165" he="46" />;上式中:<img file="740777DEST_PATH_IMAGE004.GIF" wi="17" he="18" />为方位角,<img file="469698DEST_PATH_IMAGE005.GIF" wi="14" he="18" />为横滚角,<img file="486196DEST_PATH_IMAGE015.GIF" wi="25" he="25" />为<img file="921856DEST_PATH_IMAGE016.GIF" wi="165" he="21" />,<img file="693503DEST_PATH_IMAGE017.GIF" wi="25" he="25" />为<img file="581825DEST_PATH_IMAGE018.GIF" wi="165" he="21" />,<img file="416662DEST_PATH_IMAGE019.GIF" wi="26" he="25" />为<img file="706829DEST_PATH_IMAGE020.GIF" wi="162" he="21" />,<img file="587060DEST_PATH_IMAGE021.GIF" wi="25" he="25" />为<img file="962678DEST_PATH_IMAGE022.GIF" wi="165" he="21" />;若井斜角是大角度井斜角,方位角和横滚角的提取公式如下:<img file="383295DEST_PATH_IMAGE023.GIF" wi="121" he="101" />;上式中,<img file="793548DEST_PATH_IMAGE004.GIF" wi="17" he="18" />为方位角,<img file="844680DEST_PATH_IMAGE005.GIF" wi="14" he="18" />为横滚角,<img file="769911DEST_PATH_IMAGE024.GIF" wi="25" he="25" />为<img file="931902DEST_PATH_IMAGE025.GIF" wi="61" he="21" />,<img file="931082DEST_PATH_IMAGE026.GIF" wi="24" he="25" />为<img file="215433DEST_PATH_IMAGE027.GIF" wi="64" he="20" />,<img file="562713DEST_PATH_IMAGE028.GIF" wi="25" he="25" />为<img file="528395DEST_PATH_IMAGE029.GIF" wi="58" he="21" />,<img file="647661DEST_PATH_IMAGE030.GIF" wi="25" he="25" />为<img file="837334DEST_PATH_IMAGE031.GIF" wi="61" he="20" />。
地址 030051 山西省太原市尖草坪区学院路3号