课题学习(十五)----阅读《测斜仪旋转姿态测量信号处理方法》论文

本文主要是介绍课题学习(十五)----阅读《测斜仪旋转姿态测量信号处理方法》论文,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、 论文内容

1.1 摘要

为准确测量旋转钻井时的钻具姿态,提出了一种新的信号处理方法。测斜仪旋转时,垂直于其旋转轴方向加速度计的输出信号中重力加速度信号分量具有周期性特征,以及非周期性离心加速度分量频率低于重力加速度信号分量频率的特征,用高通或带通滤波将非周期性离心加速度分量分离;对于常规滤波方法无法滤除的通频带内振动信号,增加一个对置加速度计,采取不同隔振措施补偿振动信号,滤除通频带内的振动信号。仿真分析表明,测斜仪旋转时离心力作用所导致的测斜仪加速度信息失真可以被高通或带通滤波器有效滤除;与旋转频率耦合的一个振动信号可以分解为2个频率的振动信号;采取不同隔振措施的对置加速度计,可以对振动信号进行补偿,滤除通频带内的振动信号。台架试验也证实,增加一个对置加速度计,并采取隔振措施棳可以补偿振动信号,滤除通频带内的振动信号。该处理方法为准确计算井斜角和工具面角等钻具姿态参数提供了技术手段。

1.2 加速度计的动态输出信号

   原论文中的测斜仪测斜原理部分就不再赘述,直接从下文开始。
  当加速度计静置时,输出信号基本为重力加速度信号;而当存在振动或因旋转产生离心力时,加速度计的输出信号中便混合了振动或离心加速度信号,即: s f = s g + s v + s e s_f=s_g+s_v+s_e sf=sg+sv+se
  上式中, s g s_g sg为重力加速度信号, s v s_v sv为振动加速度信号, s e s_e se为离心加速度信号。

1.3 重力加速度分量的提取方法

1.3.1 去除离心加速度

  井斜角传感器的三轴加速度计输出信号为:
s x = s g x + s v x + s e x s_x=s_{gx}+s_{vx}+s_{ex} sx=sgx+svx+sex s y = s g y + s v y + s e y s_y=s_{gy}+s_{vy}+s_{ey} sy=sgy+svy+sey s z = s g z + s v z + s e z s_z=s_{gz}+s_{vz}+s_{ez} sz=sgz+svz+sez
   s x , s y , s z 分别是 s_x,s_y,s_z分别是 sx,sy,sz分别是 x , y , z x,y,z x,y,z加速度计的输出信号。
  当井斜角传感器随钻柱旋转时,如图所示,如果z轴平行于钻柱旋转
轴线,则在轴方向不存在离心力,即 e z = 0 e_z=0 ez=0。此时,z轴方向重力分量不随转动改变,如果井斜角固定,其为常值,而z轴方向振动加速度可通过低通滤波滤除,于是有 s ˉ z = s g z \bar{s}_z=s_{gz} sˉz=sgz其中 s ˉ z \bar{s}_z sˉz为z轴加速度计低通滤波后输出信号。因为重力加速度已知,此时以静态测量时g的统计值为准,可以用下式计算井斜角: c o s α = s ˉ z g cos \alpha = \frac{\bar{s}_z}{g} cosα=gsˉz
  设 g x y g_{xy} gxy为重力加速度在xy平面内分量,则有: g x = − g x y c o s θ g_x=-g_{xy}cos\theta gx=gxycosθ g y = − g x y s i n θ g_y=-g_{xy}sin\theta gy=gxysinθ
  旋转钻进时,工具面角因旋转而发生周期变化,因此由上式可知,x,y轴重力加速度分量也随之周期变化,其角频率与旋转角速度有关。因此,该情况下仅采用低通滤波并不能得到真实的重力加速度分量。
  钻柱匀速旋转时,离心加速度为常值,x、y轴重力加速度分量为周期信号,因此可用一个高通滤波器滤去直流成分,将离心加速度分量去除。于是有 :
s x ′ = s v x − s g x y c o s θ s_x' = s_{v_x}-s_{g_{xy}}cos\theta sx=svxsgxycosθ s y ′ = s v y − s g x y s i n θ s_y' = s_{v_y}-s_{g_{xy}}sin\theta sy=svysgxysinθ

1.3.2 去除振动加速度

  不同于钻柱静止或非旋转状态,钻柱旋转时,x、y轴上的振动信号分量耦合了旋转角度的变化。设 v x y v_{xy} vxy为振动加速度在xy平面内分量, v x y v_{xy} vxy相对于高边的夹角为 ϕ \phi ϕ ω v \omega_v ωv为振动角频率,则有: v x = v x y s i n ( w v t + β ) c o s ( θ − ϕ ) v_x=v_{xy}sin(w_vt+\beta)cos(\theta-\phi) vx=vxysin(wvt+β)cos(θϕ) v y = v x y s i n ( w v t + β ) s i n ( θ − ϕ ) v_y=v_{xy}sin(w_vt+\beta)sin(\theta-\phi) vy=vxysin(wvt+β)sin(θϕ)
  上式表明,一个周期振动信号在x、y轴上分量的幅值随钻柱旋转角度呈周期变化。设旋转角速度为 ω r \omega_r ωr,则有 θ = ω r t + θ 0 \theta=\omega_rt+\theta_0 θ=ωrt+θ0其中 θ 0 \theta_0 θ0为t=0时刻的工具面角,可做如下变形 :
在这里插入图片描述
  上式表明,在钻柱旋转时,一个周期振动信号在x、y轴上的分量,耦合了旋转角频率,可分别用2个不同频率、不同相位的半幅值周期振动信号等效 。
  如果振动频率远大于钻柱旋转角速度( ω v ≥ 3 ω r \omega_v≥3\omega_r ωv3ωr),则振动信号所等效分解的2个分量的频率均远离重力加速度分量的变化频率,因此可以采,用合适的滤波器滤除振动信号。
  当 ω v ≈ 2 ω r \omega_v≈2\omega_r ωv2ωr,即振动信号的角频率接近钻柱旋转角频率的2倍时,振动信号中的低频率分量的角频率将接近旋转角速度,会对重力加
速度分量造成严重干扰,而且采用常规滤波方法很难滤除这种干扰。
  当 ω v ≈ ω r \omega_v≈\omega_r ωvωr时,振动信号低频分量接近于0,高频分,接近于 2 ω r 2\omega_r 2ωr,此时振动信号可以被有效滤除。
  由于旋转角速度较低,也就是说x、y轴方向重力加速度分量的变化频率较低(一般不超过4Hz),则 2 ω r + B 2\omega_r+B 2ωr+B以下的振动信号,由于低频分量角频率( ω v − ω r \omega_v-\omega_r ωvωr)与旋转角速度的绝对差值较小,一般难于完全滤除,而且振动信号的频率分量越接近旋转角速度,越难以滤除。此时用滤波后信号计算工具面角或井斜角会有较大误差。
  解决该问题的有效方法是椇增加1个加速度传感器与x或y轴背对放置,并对其采取隔振措施,使其对振动的敏感强度与对称轴有明显不同。例如与轴背对放置第4个加速度计-x,则其输出信号 s − x = s − g x − k s v x + s e x s_{-x}=s_{-g{x}}-ks_{v{x}}+s_{e{x}} sx=sgxksvx+sex,其中k为系数,于是有 s v x = ( s x + s − x − 2 s e x ) / ( 1 − k ) s_{v{x}}=(s_x+s_{-x}-2s_{e{x}})/(1-k) svx=(sx+sx2sex)/(1k)
  高通滤波器可将 s e x s_{e{x}} sex滤除,因此可写成 s v x = ( s x ′ + s − x ′ ) / ( 1 − k ) s_{v{x}}=(s_x'+s_{-x}')/(1-k) svx=(sx+sx)/(1k)
  现在问题的关键是确定系数k。测斜仪不旋转时, s g x s_{g{x}} sgx为常数, s e x s_{e{x}} sex为0,因此k可由下式求得 k = − s x ′ / s − x ′ k=-s_x'/s_{-x}' k=sx/sx
  注意,本论文没有讨论电路以及采样中的白噪声,因为白噪声可以通过低通或者电筒滤波器。所以一般在进行数据处理前,均需要进行低通或者带通滤波器处理。
  此外棳在进行振动信号补偿分析时,假定2个传感器的离心加速度大小相同,方向相反。事实上,由于装配等其他原因,2个传感器的离心加速度不一
定大小相同,方向相反。由于通常情况下,旋转钻进时,钻柱近似匀速旋转,离心加速度近似恒定,可利用高通或带通滤波器滤除,上文给出的方法仍然有效。

1.3 仿真及实验验证

在这里插入图片描述

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

二、往期回顾

课题学习(一)----静态测量
课题学习(二)----倾角和方位角的动态测量方法(基于磁场的测量系统)
课题学习(三)----倾角和方位角的动态测量方法(基于陀螺仪的测量系统)
课题学习(四)----四元数解法
课题学习(五)----阅读论文《抗差自适应滤波的导向钻具动态姿态测量方法》
课题学习(六)----安装误差校准、实验方法
课题学习(七)----粘滑运动的动态算法
课题学习(八)----卡尔曼滤波动态求解倾角、方位角
课题学习(九)----阅读《导向钻井工具姿态动态测量的自适应滤波方法》论文笔记
课题学习(十)----阅读《基于数据融合的近钻头井眼轨迹参数动态测量方法》论文笔记
课题学习(十一)----阅读《Attitude Determination with Magnetometers and Accelerometers to Use in Satellite》
课题学习(十二)----阅读《Extension of a Two-Step Calibration Methodology to Include Nonorthogonal Sensor Axes》
课题学习(十三)----阅读《Calibration of Strapdown Magnetometers in Magnetic Field Domain》论文笔记
课题学习(十四)----三轴加速度计+三轴陀螺仪传感器-ICM20602

这篇关于课题学习(十五)----阅读《测斜仪旋转姿态测量信号处理方法》论文的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/464311

相关文章

Java学习手册之Filter和Listener使用方法

《Java学习手册之Filter和Listener使用方法》:本文主要介绍Java学习手册之Filter和Listener使用方法的相关资料,Filter是一种拦截器,可以在请求到达Servl... 目录一、Filter(过滤器)1. Filter 的工作原理2. Filter 的配置与使用二、Listen

Pandas统计每行数据中的空值的方法示例

《Pandas统计每行数据中的空值的方法示例》处理缺失数据(NaN值)是一个非常常见的问题,本文主要介绍了Pandas统计每行数据中的空值的方法示例,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是空值?为什么要统计空值?准备工作创建示例数据统计每行空值数量进一步分析www.chinasem.cn处

Windows 上如果忘记了 MySQL 密码 重置密码的两种方法

《Windows上如果忘记了MySQL密码重置密码的两种方法》:本文主要介绍Windows上如果忘记了MySQL密码重置密码的两种方法,本文通过两种方法结合实例代码给大家介绍的非常详细,感... 目录方法 1:以跳过权限验证模式启动 mysql 并重置密码方法 2:使用 my.ini 文件的临时配置在 Wi

MySQL重复数据处理的七种高效方法

《MySQL重复数据处理的七种高效方法》你是不是也曾遇到过这样的烦恼:明明系统测试时一切正常,上线后却频频出现重复数据,大批量导数据时,总有那么几条不听话的记录导致整个事务莫名回滚,今天,我就跟大家分... 目录1. 重复数据插入问题分析1.1 问题本质1.2 常见场景图2. 基础解决方案:使用异常捕获3.

最详细安装 PostgreSQL方法及常见问题解决

《最详细安装PostgreSQL方法及常见问题解决》:本文主要介绍最详细安装PostgreSQL方法及常见问题解决,介绍了在Windows系统上安装PostgreSQL及Linux系统上安装Po... 目录一、在 Windows 系统上安装 PostgreSQL1. 下载 PostgreSQL 安装包2.

SQL中redo log 刷⼊磁盘的常见方法

《SQL中redolog刷⼊磁盘的常见方法》本文主要介绍了SQL中redolog刷⼊磁盘的常见方法,将redolog刷入磁盘的方法确保了数据的持久性和一致性,下面就来具体介绍一下,感兴趣的可以了解... 目录Redo Log 刷入磁盘的方法Redo Log 刷入磁盘的过程代码示例(伪代码)在数据库系统中,r

Python实现图片分割的多种方法总结

《Python实现图片分割的多种方法总结》图片分割是图像处理中的一个重要任务,它的目标是将图像划分为多个区域或者对象,本文为大家整理了一些常用的分割方法,大家可以根据需求自行选择... 目录1. 基于传统图像处理的分割方法(1) 使用固定阈值分割图片(2) 自适应阈值分割(3) 使用图像边缘检测分割(4)

Java中Switch Case多个条件处理方法举例

《Java中SwitchCase多个条件处理方法举例》Java中switch语句用于根据变量值执行不同代码块,适用于多个条件的处理,:本文主要介绍Java中SwitchCase多个条件处理的相... 目录前言基本语法处理多个条件示例1:合并相同代码的多个case示例2:通过字符串合并多个case进阶用法使用

Python中__init__方法使用的深度解析

《Python中__init__方法使用的深度解析》在Python的面向对象编程(OOP)体系中,__init__方法如同建造房屋时的奠基仪式——它定义了对象诞生时的初始状态,下面我们就来深入了解下_... 目录一、__init__的基因图谱二、初始化过程的魔法时刻继承链中的初始化顺序self参数的奥秘默认

html5的响应式布局的方法示例详解

《html5的响应式布局的方法示例详解》:本文主要介绍了HTML5中使用媒体查询和Flexbox进行响应式布局的方法,简要介绍了CSSGrid布局的基础知识和如何实现自动换行的网格布局,详细内容请阅读本文,希望能对你有所帮助... 一 使用媒体查询响应式布局        使用的参数@media这是常用的