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

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

一、 论文内容

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

相关文章

Linux换行符的使用方法详解

《Linux换行符的使用方法详解》本文介绍了Linux中常用的换行符LF及其在文件中的表示,展示了如何使用sed命令替换换行符,并列举了与换行符处理相关的Linux命令,通过代码讲解的非常详细,需要的... 目录简介检测文件中的换行符使用 cat -A 查看换行符使用 od -c 检查字符换行符格式转换将

SpringBoot实现数据库读写分离的3种方法小结

《SpringBoot实现数据库读写分离的3种方法小结》为了提高系统的读写性能和可用性,读写分离是一种经典的数据库架构模式,在SpringBoot应用中,有多种方式可以实现数据库读写分离,本文将介绍三... 目录一、数据库读写分离概述二、方案一:基于AbstractRoutingDataSource实现动态

Java中的String.valueOf()和toString()方法区别小结

《Java中的String.valueOf()和toString()方法区别小结》字符串操作是开发者日常编程任务中不可或缺的一部分,转换为字符串是一种常见需求,其中最常见的就是String.value... 目录String.valueOf()方法方法定义方法实现使用示例使用场景toString()方法方法

Java中List的contains()方法的使用小结

《Java中List的contains()方法的使用小结》List的contains()方法用于检查列表中是否包含指定的元素,借助equals()方法进行判断,下面就来介绍Java中List的c... 目录详细展开1. 方法签名2. 工作原理3. 使用示例4. 注意事项总结结论:List 的 contain

macOS无效Launchpad图标轻松删除的4 种实用方法

《macOS无效Launchpad图标轻松删除的4种实用方法》mac中不在appstore上下载的应用经常在删除后它的图标还残留在launchpad中,并且长按图标也不会出现删除符号,下面解决这个问... 在 MACOS 上,Launchpad(也就是「启动台」)是一个便捷的 App 启动工具。但有时候,应

SpringBoot日志配置SLF4J和Logback的方法实现

《SpringBoot日志配置SLF4J和Logback的方法实现》日志记录是不可或缺的一部分,本文主要介绍了SpringBoot日志配置SLF4J和Logback的方法实现,文中通过示例代码介绍的非... 目录一、前言二、案例一:初识日志三、案例二:使用Lombok输出日志四、案例三:配置Logback一

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很

mysql出现ERROR 2003 (HY000): Can‘t connect to MySQL server on ‘localhost‘ (10061)的解决方法

《mysql出现ERROR2003(HY000):Can‘tconnecttoMySQLserveron‘localhost‘(10061)的解决方法》本文主要介绍了mysql出现... 目录前言:第一步:第二步:第三步:总结:前言:当你想通过命令窗口想打开mysql时候发现提http://www.cpp

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T

MySQL INSERT语句实现当记录不存在时插入的几种方法

《MySQLINSERT语句实现当记录不存在时插入的几种方法》MySQL的INSERT语句是用于向数据库表中插入新记录的关键命令,下面:本文主要介绍MySQLINSERT语句实现当记录不存在时... 目录使用 INSERT IGNORE使用 ON DUPLICATE KEY UPDATE使用 REPLACE