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

一、 论文内容

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

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/224742.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

渲染(iOS渲染过程解析)

渲染 渲染原理 一个硬核硬件科普视频 CPU和GPU CPU(Central Processing Unit):现代计算机整个系统的运算核心、控制核心,适合串行计算。GPU(Graphics Processing Unit):可进行绘图运算工作的…

系列四、过滤器简介

一、简介 1.1、概述 过滤器作为JavaWEB的三大组件(Servlet程序、Filter过滤器、Listener监听器),它的主要功能是用来拦截请求的,当客户端要访问某个资源时,先来到配置好的过滤器,过滤器可以在用户访问某个…

Docker架构、镜像操作和容器操作

一、docker基本管理和概念 1、概念 docker:开源的应用容器引擎。基于go语言开发的。运行在Linux系统中的开源的轻量级的“虚拟机” docker的容器技术可用在一台主机上轻松到达为任何应用创建一个轻量级到的,可移植的,自给自足的容器 dock…

基于remix+metamask+ganache的智能合约部署调用

在我们部署合约时为了让它更接近真实区块链去中心化体验,我们需要调用小狐狸(Metamask)来进行真实交易,而metamask里没有内置虚拟测试币,我们需要进行调用Ganache来添加带有虚拟测试币的账号。以上就是三者的关系&…

编程实战:类C语法的编译型脚本解释器(十)编译表达式

系列入口:编程实战:类C语法的编译型脚本解释器(九)编译语句 本文介绍表达式的编译。 一、代码概览 表达式的编译就是不断获取下一个标识符,直到遇到不属于表达式的东西。 完整代码如下: Expression* GetExp…

Java+Swing: (jframe自定义图标和居中显示) 整理1

package com.test;import javax.swing.*; import java.awt.*; import java.net.URL;/*** Author:xiexu* Date:2023/12/3 19:13*/ public class JframeTest {JFrame jFrame;JButton jButton;public JframeTest() {// 容器组件(jframe, jpanel,…

深圳冬季穿搭造型

深圳冬季穿搭造型 今天是2023年11月17日,北方在11月初就降下了大雪,那是我的老家,无比怀念。 而在深圳,冬天是体会不到那份冷冽的寒冷的,所以在深圳,每天的晚上,我都会出去散步,昨…

Kubernetes Service控制器详解以及切换为ipvs代理模式

文章目录 一、Service 存在的意义二、Pod与Service的关系三、Service定义与创建四、Service三种常用类型五、Service代理模式六、切换Service代理模式七、service总体工作流程八、kube-proxy ipvs和iptables的异同九、Service DNS名称 一、Service 存在的意义 service的引入主…

推荐10款App安全测试工具

移动互联网时代,我们的生活和工作深受 App 影响。伴随移动 App 的广泛应用,App 安全日益重要。本文介绍了 App 开发可能用到的安全测试工具。 当今, 全球移动用户大约超过37亿。 Google Play 上大约有 220 万个 App, 苹果App Sto…

二极管是什么

二极管 电子元器件百科 文章目录 二极管前言一、二极管是什么二、二极管的类别三、二极管的应用实例四、二极管的作用原理总结前言 二极管是一种重要的电子器件,通过其整流行为和管理方向的特性,可以在电路中实现电流控制和电压整流等功能。 一、二极管是什么 二极管是一种…

idea__SpringBoot微服务03——yaml(新注解)(新的依赖)

yaml 一、数据格式二、注入配置文件(yaml注入)(新注解ConfigurationProperties)三、注入配置文件(properties注入)(新注解PropertySource)四、yaml配置文件占位符${}五 、yaml跟prop…

【虚拟机】Docker基础 【二】【数据卷和挂载本地目录】

2.2.数据卷 容器是隔离环境,容器内程序的文件、配置、运行时产生的容器都在容器内部,我们要读写容器内的文件非常不方便。大家思考几个问题: 如果要升级MySQL版本,需要销毁旧容器,那么数据岂不是跟着被销毁了&#x…

智能优化算法应用:基于蜜獾算法无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用:基于蜜獾算法无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用:基于蜜獾算法无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.蜜獾算法4.实验参数设定5.算法结果6.参考文献7.MATLAB…

深入理解 new 操作符:创建对象的秘密武器(下)

🤍 前端开发工程师(主业)、技术博主(副业)、已过CET6 🍨 阿珊和她的猫_CSDN个人主页 🕠 牛客高级专题作者、在牛客打造高质量专栏《前端面试必备》 🍚 蓝桥云课签约作者、已在蓝桥云…

SQL Server 数据库,创建触发器避免数据被更改

5.4触发器 触发器是一种特殊类型的存储过程,当表中的数据发生更新时将自动调用,以响应INSERT、 UPDATE 或DELETE 语句。 5.4.1什么是触发器 1.触发器的概念 触发器是在对表进行插入、更新或删除操作时自动执行的存储过程,触发器通常用于强…

Elastcsearch:通过 Serverless 提供更多服务

作者:Ken Exner 人们使用 Elasticsearch 解决最大数据挑战的方式一直令我们感到惊讶。 从超过 40 亿次下载、70,000 次提交、1,800 名贡献者以及我们全球社区的反馈中可以清楚地看出这一点。 Elastic 在广泛的用例中发挥的作用促使我们简化复杂性,让搜索…

9.基于SpringBoot3+I18N实现国际化

1. 新建资源文件 在resources目录下新建目录i18n, 然后 新建messages_en.properties文件 user.login.erroraccount or password error!新建messages_zh_CN.properties文件 user.login.error帐户或密码错误!2. 新建LocaleConfig.java文件 Configurati…

gpt3、gpt2与gpt1区别

参考:深度学习:GPT1、GPT2、GPT-3_HanZee的博客-CSDN博客 Zero-shot Learning / One-shot Learning-CSDN博客 Zero-shot(零次学习)简介-CSDN博客 GPT-2 模型由多层单向transformer的解码器部分构成,本质上是自回归模型…

Hazelcast分布式内存网格(IMDG)基本使用,使用Hazelcast做分布式内存缓存

文章目录 一、Hazelcast简介1、Hazelcast概述2、Hazelcast之IMDG3、数据分区 二、Hazelcast配置1、maven坐标2、集群搭建(1)组播自动搭建 3、客户端4、集群分组5、其他配置 三、Hazelcast分布式数据结构1、IMap2、IQueue:队列3、MultiMap4、I…

MySQL和MongoDB简介以及它们之间的区别

本文主要介绍MySQL和MongoDB的简介以及它们之间的区别。 目录 MySQL简介MySQL的优缺点MySQL的应用场景MongoDB简介MongoDB的优缺点MongoDB的应用场景MySQL和MongoDB的区别 MySQL简介 MySQL是一种开源的关系型数据库管理系统,是世界上最流行的数据库之一。它支持多用…