CN109581525A - 旋转加速度计式重力梯度敏感器原始采样频率的选择方法 - Google Patents
旋转加速度计式重力梯度敏感器原始采样频率的选择方法 Download PDFInfo
- Publication number
- CN109581525A CN109581525A CN201811406519.4A CN201811406519A CN109581525A CN 109581525 A CN109581525 A CN 109581525A CN 201811406519 A CN201811406519 A CN 201811406519A CN 109581525 A CN109581525 A CN 109581525A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- gravity gradient
- gradient sensor
- accelerometer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,其技术特点在于:包括以下步骤:步骤1、对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号进行离散采样;步骤2、对步骤1得到的一组信号通过抽样的方法得到多组不同采样频率的系统输出数字信号;步骤3、对步骤2得到的多组信号分别计算其在旋转频率二倍频的功率谱密度值;步骤4、对步骤3中多组信号的采样频率倒数和对应的功率谱密度值利用最小二乘的方法进行一阶线性拟合,得到一次项系数a1和常值项系数a0;步骤5、确定最佳原始采样频率。本发明能够输出信号的二倍频附近白噪声能量最小,且硬件实现代价最低。
Description
技术领域
本发明属于重力梯度敏感器技术领域,涉及重力梯度敏感器原始采样频率的选择方法,尤其是一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法。
背景技术
旋转加速度计式重力梯度仪是对地球表面微小重力梯度变化进行连续测量的仪器。如图2所示,作为核心敏感器的重力梯度测量组件是由安装在旋转圆盘上的四个等距离成对安装的石英挠性加速度计组成。每个加速度计的检测质心到圆盘中心距离相等,加速度计敏感轴正切于圆,两组加速度计正交安装,且每组加速度计敏感轴反向安置。其中,1号加速度计与3号相反,2号与4号相反。实际工作中,圆盘以恒定的角速率ω旋转,对重力梯度信号进行调制,此时四个加速度计组合输出为:
(a1+a3)-(a2+a4)=2R[(Γxx-Γyy)sin(ωt+θ)-2Γxycos(ωt+θ)]
其中,ω为圆盘旋转速率;R为加速度计中心到圆盘中心的距离;θ为初始相位。最后,加速度计组合信号通过2ω频率解调和低通滤波得到最终的重力梯度张量Γxx-Γyy和Γxy。
重力梯度信号十分微弱,在近地表的最大变化幅值仅为几百E左右(1E=10-9s-2),这就对加速度计组合输出信号在调制频率附近的本底噪声提出了很高的要求。系统设计的目标在于降低加速度计组合输出信号调制频率附近的本底噪声,提高系统测量分辨率。对系统信号进行Allan方差分析,分析结果表明在调制频率附近的噪声模式主要有白噪声和粉红噪声两种,二者的功率谱密度函数如下表示:
式中Sw(f)表示白噪声的功率谱密度幅值,N表示白噪声幅值,Sp(f)表示粉红噪声随频率的功率谱密度幅值,B表示粉红噪声幅值。
帕塞瓦尔定理表明,信号时域总能量等于频域总能量。对离散信号X则可表示为:
式中,E表示信号X的总能量,SX(f)表示信号X的功率谱密度函数,fs表示信号X的采样频率。如果信号X是白噪声模式序列,强度为N,则上式可表示为:
E=N·fs
由此可知,在信号总能量不变的条件下,提高离散信号的采样频率,可对数线性降低白噪声强度。由于白噪声在频域上是均匀分布的,即可降低调制频率点附近的白噪声强度,而提高采样频率对调制频率附近的粉红噪声没有影响。但增加系统原始采样频率不仅对精密ADC提出了更高速率的要求,并且占用大量的处理器资源,还会显著增加数据量,加重数据传输负荷。
因此,需建立一种原始采样频率参数优选方法,在降低白噪声对测量信号干扰的同时最大程度上减小数据传输压力。
发明内容
本发明的目的在于克服现有技术的不足,提出一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,用于降低重力梯度信号调制频率附近加速度计白噪声对测量信号的影响,以提高测量精度。
本发明解决其现实问题是采取以下技术方案实现的:
一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,包括以下步骤:
步骤1、对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号进行离散采样;
步骤2、对步骤1得到的一组信号通过抽样的方法得到多组不同采样频率的系统输出数字信号;
步骤3、对步骤2得到的多组信号分别计算其在旋转频率二倍频的功率谱密度值;
步骤4、对步骤3中多组信号的采样频率倒数和对应的功率谱密度值利用最小二乘的方法进行一阶线性拟合,得到一次项系数a1和常值项系数a0;
步骤5、确定最佳原始采样频率。
而且,所述步骤1的具体方法为:以采样频率为fs0的方式在静态下对正常旋转工作的重力梯度敏感器中加速度计组合输出信号进行离散采样,得到四只加速度计输出组合信号Aout0;
而且,所述步骤2的具体方法为:通过重采样或抽样的方法,将采样频率为fs0的序列Aout0重采样或抽样为一组或多组的信号Aout1、Aout2、Aout3……,其对应的采样频率为fs1、fs2、fs3……;
而且,所述步骤3的具体方法为:将采样频率分别为fs0、fs1、fs2、fs3……的加速度计组合输出信号进行功率谱密度分析,得到对应频率下加速度计组合信号在旋转频率二倍频附近的本底噪声功率谱密度幅值S0、S1、S2、S3……;
而且,所述步骤4的具体方法为:通过最小二乘的方法将采样频率的倒数序列[1/fs0 1/fs1 1/fs2 1/fs3···]和其对应的功率谱密度幅值序列[S0 S1 S2 S3···]使用最小二乘的方法进行一阶线性拟合,得到一次项系数a1和常值项系数a0,其中a1在数值上表示在采样频率为1Hz时加速度计组合输出信号在旋转频率二倍频附近白噪声的功率谱密度幅值,a0表示加速度计组合输出信号在旋转频率二倍频附近粉红噪声的功率谱密度幅值;
而且,所述步骤5的具体方法为:假定系统输出信号在旋转频率二倍频附近白噪声与粉红噪声能量之比为1/100,需通过调整采样频率fs使与a0项的比值为1/100,即满足要求,认为该采样频率fs为重力梯度敏感器最佳原始采样频率。
本发明的优点和有益效果:
本发明通过建立重力梯度敏感器输出信号白噪声和粉红噪声的模型计算出系统最佳原始采样频率,使得在该采样频率下系统输出信号的二倍频附近白噪声能量最小,且硬件实现代价最低。
附图说明
图1是本发明的重力梯度敏感器原始采样频率选择方法流程图;
图2是本发明的旋转加速度计式重力梯度敏感器原理示意图;
图3是本发明的以100Hz为采样频率的重力梯度敏感器加速度计组合信号输出离散信号时域图;
图4(a)-图4(b)是本发明的利用抽样方法得到的不同采样频率的离散信号时域图;其中,图4(a):采样频率50Hz、图4(b):采样频率25Hz;
图5(a)-图5(c)是本发明的不同频率采样的离散信号的功率谱密度分布图;其中,图5(a):采样频率100Hz;图5(b):采样频率50Hz;图5(c):采样频率25Hz。
具体实施方式
以下结合附图对本发明实施例作进一步详述:
一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,如图1所示,包括以下步骤:
步骤1、以采样频率为fs0的方式在静态下对正常旋转工作的重力梯度敏感器中加速度计组合输出信号进行离散采样,得到四只加速度计输出组合信号Aout0;
在本实施例中,对旋转加速度计式重力梯度敏感器中加速度计组合输出信号进行离散采样,以敏感器原始采样频率100Hz为例,得到相应的离散信号Aout0,其时域图如图3所示。
步骤2、通过重采样或抽样的方法,将采样频率为fs0的序列Aout0重采样或抽样为一组或多组的信号Aout1、Aout2、Aout3……,其对应的采样频率为fs1、fs2、fs3……;
在本实施例中,利用抽样的方法将离散信号Aout0的采样率改为50Hz、25Hz,其对应的离散信号Aout1、Aout2的时域图如图4(a)和图4(b)所示。
步骤3、将采样频率分别为fs0、fs1、fs2、fs3……的加速度计组合输出信号进行功率谱密度分析,得到对应频率下加速度计组合信号在旋转频率二倍频附近的本底噪声功率谱密度幅值S0、S1、S2、S3……;
在本实施例中,分别对采样频率为100Hz、50Hz、25Hz的离散信号Aout0、Aout1、Aout2做离散功率谱密度分析,得到其功率谱密度分布如图5(a)-图5(c)所示,由此得到不同采样频率下加速度计组合信号在旋转频率二倍频附近的本底噪声功率谱密度幅值分别为
步骤4、通过最小二乘的方法将采样频率的倒数序列[1/fs0 1/fs1 1/fs2 1/fs3···]和其对应的功率谱密度幅值序列[S0 S1 S2 S3···]使用最小二乘的方法进行一阶线性拟合,得到一次项系数a1和常值项系数a0,其中a1在数值上表示在采样频率为1Hz时加速度计组合输出信号在旋转频率二倍频附近白噪声的功率谱密度幅值,a0表示加速度计组合输出信号在旋转频率二倍频附近粉红噪声的功率谱密度幅值;
在本实施例中,通过最小二乘的方法将采样频率的倒数序列和其对应的功率谱密度幅值序列y=[0.0012 0.0017 0.0026]使用最小二乘的方法进行一阶线性拟合,得到一次项系数a1=0.0464和常值项系数a0=0.0007。
步骤5、假定系统输出信号在旋转频率二倍频附近白噪声与粉红噪声能量之比为1/100,需通过调整采样频率fs使与a0项的比值为1/100,即满足要求,认为该采样频率fs为重力梯度敏感器最佳原始采样频率。
在本实施例中,利用公式计算得到该重力梯度敏感器最佳原始采样频率为6628.6Hz。
需要强调的是,本发明所述实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (6)
1.一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,其特征在于:包括以下步骤:
步骤1、对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号进行离散采样;
步骤2、对步骤1得到的一组信号通过抽样的方法得到多组不同采样频率的系统输出数字信号;
步骤3、对步骤2得到的多组信号分别计算其在旋转频率二倍频的功率谱密度值;
步骤4、对步骤3中多组信号的采样频率倒数和对应的功率谱密度值利用最小二乘的方法进行一阶线性拟合,得到一次项系数和常值项系数;
步骤5、确定最佳原始采样频率。
2.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,其特征在于:所述步骤1的具体方法为:以采样频率为fs0的方式在静态下对正常旋转工作的重力梯度敏感器中加速度计组合输出信号进行离散采样,得到四只加速度计输出组合信号Aout0。
3.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,其特征在于:所述步骤2的具体方法为:通过重采样或抽样的方法,将采样频率为fs0的序列Aout0重采样或抽样为一组或多组的信号Aout1、Aout2、Aout3……,其对应的采样频率为fs1、fs2、fs3……。
4.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,其特征在于:所述步骤3的具体方法为:将采样频率分别为fs0、fs1、fs2、fs3……的加速度计组合输出信号进行功率谱密度分析,得到对应频率下加速度计组合信号在旋转频率二倍频附近的本底噪声功率谱密度幅值S0、S1、S2、S3……。
5.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,其特征在于:所述步骤4的具体方法为:通过最小二乘的方法将采样频率的倒数序列[1/fs0 1/fs1 1/fs2 1/fs3···]和其对应的功率谱密度幅值序列[S0 S1 S2 S3···]使用最小二乘的方法进行一阶线性拟合,得到一次项系数a1和常值项系数a0,其中a1在数值上表示在采样频率为1Hz时加速度计组合输出信号在旋转频率二倍频附近白噪声的功率谱密度幅值,a0表示加速度计组合输出信号在旋转频率二倍频附近粉红噪声的功率谱密度幅值。
6.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器原始采样频率的选择方法,其特征在于:所述步骤5的具体方法为:假定系统输出信号在旋转频率二倍频附近白噪声与粉红噪声能量之比为1/100,需通过调整采样频率fs使与a0项的比值为1/100,即满足要求,认为该采样频率fs为重力梯度敏感器最佳原始采样频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811406519.4A CN109581525B (zh) | 2018-11-23 | 2018-11-23 | 旋转加速度计式重力梯度敏感器原始采样频率的选择方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811406519.4A CN109581525B (zh) | 2018-11-23 | 2018-11-23 | 旋转加速度计式重力梯度敏感器原始采样频率的选择方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109581525A true CN109581525A (zh) | 2019-04-05 |
CN109581525B CN109581525B (zh) | 2020-06-19 |
Family
ID=65924247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811406519.4A Active CN109581525B (zh) | 2018-11-23 | 2018-11-23 | 旋转加速度计式重力梯度敏感器原始采样频率的选择方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109581525B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111830589A (zh) * | 2020-06-18 | 2020-10-27 | 中国船舶重工集团公司第七0七研究所 | 旋转加速度计重力梯度敏感器冗余测量信号离线融合方法 |
CN113447994A (zh) * | 2021-06-25 | 2021-09-28 | 中国船舶重工集团公司第七0七研究所 | 一种基于旋转加速度计原理的重力梯度全张量测量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0735522A1 (en) * | 1995-03-27 | 1996-10-02 | Hewlett-Packard Company | Disk data recording and reproducing device |
CN102680948A (zh) * | 2012-05-15 | 2012-09-19 | 东南大学 | 一种线性调频信号调频率和起始频率估计方法 |
CN105044795A (zh) * | 2015-05-29 | 2015-11-11 | 中国船舶重工集团公司第七0七研究所 | 利用光纤陀螺减小旋转加速度计重力梯度仪测量误差的方法 |
US20170275989A1 (en) * | 2014-09-10 | 2017-09-28 | Fracture ID, Inc. | Apparatus and method using measurements taken while drilling cement to obtain absolute values of mechanical rock properties along a borehole |
CN108287372A (zh) * | 2017-12-21 | 2018-07-17 | 中国船舶重工集团公司第七0七研究所 | 一种基于惯性技术的重力梯度敏感器转速选择方法 |
-
2018
- 2018-11-23 CN CN201811406519.4A patent/CN109581525B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0735522A1 (en) * | 1995-03-27 | 1996-10-02 | Hewlett-Packard Company | Disk data recording and reproducing device |
CN102680948A (zh) * | 2012-05-15 | 2012-09-19 | 东南大学 | 一种线性调频信号调频率和起始频率估计方法 |
US20170275989A1 (en) * | 2014-09-10 | 2017-09-28 | Fracture ID, Inc. | Apparatus and method using measurements taken while drilling cement to obtain absolute values of mechanical rock properties along a borehole |
CN105044795A (zh) * | 2015-05-29 | 2015-11-11 | 中国船舶重工集团公司第七0七研究所 | 利用光纤陀螺减小旋转加速度计重力梯度仪测量误差的方法 |
CN108287372A (zh) * | 2017-12-21 | 2018-07-17 | 中国船舶重工集团公司第七0七研究所 | 一种基于惯性技术的重力梯度敏感器转速选择方法 |
Non-Patent Citations (1)
Title |
---|
杨 晔 等: ""旋转加速度计式重力梯度仪输出解调与滤波"", 《中国惯性技术学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111830589A (zh) * | 2020-06-18 | 2020-10-27 | 中国船舶重工集团公司第七0七研究所 | 旋转加速度计重力梯度敏感器冗余测量信号离线融合方法 |
CN111830589B (zh) * | 2020-06-18 | 2022-07-26 | 中国船舶重工集团公司第七0七研究所 | 旋转加速度计重力梯度敏感器冗余测量信号离线融合方法 |
CN113447994A (zh) * | 2021-06-25 | 2021-09-28 | 中国船舶重工集团公司第七0七研究所 | 一种基于旋转加速度计原理的重力梯度全张量测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109581525B (zh) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109856689B (zh) | 一种超导航磁梯度张量数据抑噪处理方法和系统 | |
CN105717553B (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
CN109581524A (zh) | 一种旋转加速度计式重力梯度敏感器动态测量解调方法 | |
CN104459826B (zh) | 旋转加速度计重力梯度仪重力梯度信号仿真方法 | |
Peng et al. | A new type of tri-axial accelerometers with high dynamic range MEMS for earthquake early warning | |
CN109581525A (zh) | 旋转加速度计式重力梯度敏感器原始采样频率的选择方法 | |
JPH11506198A (ja) | 低次グラディオメータを用いて高次グラディオメータを得るための方法及び装置 | |
CN101251556A (zh) | 正弦信号四参数检测方法和虚拟仪器信号检测装置 | |
WO2020140378A1 (zh) | 一种旋转加速度计重力梯度仪运动误差事后补偿方法 | |
CN105044798A (zh) | 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法 | |
CN108287372B (zh) | 一种基于惯性技术的重力梯度敏感器转速选择方法 | |
CN108535544A (zh) | 一种基于正交相位解调技术的高精度数字相位测量方法 | |
CN107153224A (zh) | 检波器动态性能综合测试与评价方法 | |
CN104914485B (zh) | 加速度计在1g输入激励下的分辨率测试方法 | |
CN115236738A (zh) | 一种基于汉克尔变换的勒夫波频散提取方法 | |
Hutt et al. | Guidelines for standardized testing of broadband seismometers and accelerometers | |
CN106383366A (zh) | 一种旋转簧片式重力梯度仪 | |
CN108931824A (zh) | 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法 | |
Webb et al. | High-resolution seafloor absolute pressure gauge measurements using a better counting method | |
CN206146426U (zh) | 一种提供稳定地磁场环境的装置 | |
CN105674956A (zh) | 一种风电机组塔体倾斜度测量系统 | |
CN104484545B (zh) | 一种基于方差贡献率的泄流结构振动响应动态融合方法 | |
CN109085654B (zh) | 一种旋转加速度计重力梯度仪数字建模仿真方法 | |
Byrne | Instrument noise in seismometers | |
CN205749941U (zh) | 重力加速度的非水平测量系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |