CN104459826B - 旋转加速度计重力梯度仪重力梯度信号仿真方法 - Google Patents
旋转加速度计重力梯度仪重力梯度信号仿真方法 Download PDFInfo
- Publication number
- CN104459826B CN104459826B CN201410610089.3A CN201410610089A CN104459826B CN 104459826 B CN104459826 B CN 104459826B CN 201410610089 A CN201410610089 A CN 201410610089A CN 104459826 B CN104459826 B CN 104459826B
- Authority
- CN
- China
- Prior art keywords
- signal
- accelerometer
- analog
- demodulation
- gravity
- 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.)
- Active
Links
Landscapes
- Navigation (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种旋转加速度计重力梯度仪重力梯度信号仿真方法,根据万有引力公式计算出圆盘上四只加速度计感应到的引力加速度计信号,四组加速度信号分别与各种误差数据相加,计算机把计算得到的含有误差的加速度数据发送给模数转换器,然后经过电流转电压信号转换、电压放大和加减操作,对加法信号和减法信号分别进行1倍频解调和固定频率信号解调,根据解调结果判断是否调整加速度计标度因数,如果解调结果满足要求,则对仿真***总输出信号进行2倍频梯度解调,得到重力梯度值。通过本发明,可以详细研究旋转加速度计重力梯度仪的信号特征、参数调整、误差补偿以及信号处理方法。
Description
技术领域
本发明涉及一种旋转加速度计重力梯度仪重力梯度信号仿真方法,尤其涉及一种基于精密仪器模拟真实***环境的重力梯度信号仿真方法。
背景技术
地球重力场反映了地球内部物质分布和组成信息,通过精确测量地球重力场可以反演出物质分布。高精度重力梯度测量对于空间科学、地球科学、地质科学的发展以及在惯性导航等方面具有非常重要的意义,同时重力梯度测量被认为是资源勘探的最有效的手段之一,对于基础地质调查、基础地质研究、大地测量、油气矿藏资源勘查等领域具有重要的应用前景,航空重力梯度测量对山区、无人区、海岛和沿海大陆架部分的基础数据获取具有重要应用价值。
在直线加速度场中可以通过测量两点之间的加速度之差测量重力的变化,就是重力梯度,重力梯度反映的是重力加速度的空间微分,能够反映场源体的细节,具有比重力更高的分辨率。采用差分形式的重力梯度仪不受重力仪运动速度及加速度的影响,在实际的重力梯度仪***中,梯度信号非常微弱,这对重力梯度敏感器、信号调理、噪声抑制等技术提出了严峻挑战。
本世纪初,澳大利亚BHP公司与美国Lockheed Martin公司联合研制了Falcon航空重力梯度仪,并将其安装在Cessna Grand Caravan飞机上进行了一系列能源勘探实验,取得了很好的效果,并投入了商业运行。目前,国际上正在研制的重力梯度测量***主要包括旋转加速度计重力梯度测量***、超导重力梯度测量***、冷原子重力梯度测量***等,而真正成功投入勘探应用的重力梯度测量***只有基于美国Lockheed Martin公司研制的旋转加速度计重力全张量重力梯度仪,还有在此基础上发展起来的澳大利亚BHP公司研制的Falcon部分张量重力梯度仪和Bell Geospace(现在的Lockheed Martin)研制的Air-FTGTM与ARKeX公司研制的FTGeX全张量重力梯度仪。目前Falcon精度达到了8.1E,Air-FTGTM的探测精度达到8E,FTGeX的敏感度达到7E/√Hz。
发明内容
技术问题:本发明提出一种利用高精密仪表实现模拟真实***环境,可以帮助研制人员分析旋转加速度计重力梯度仪***误差源、研究滤波算法,有助于为旋转加速度计重力梯度仪电路设计提供验证测试平台的旋转加速度计重力梯度仪重力梯度信号仿真方法。
技术方案:本发明的旋转加速度计重力梯度仪重力梯度信号仿真方法,包括以下步骤:
1)对仿真***参数进行初始化。
2)将四只高精度加速度计均匀正交放置在旋转圆盘上,相对两个加速度计的输入轴方向相反,相邻两个加速度计的输入轴相互垂直,且四个加速度计的输入轴方向正切于圆盘边沿,将圆盘坐标系记为OXYZ,原点为O,圆盘圆心位于坐标系原点,圆盘旋转轴垂直于圆盘平面。
3)将待检测物体置于步骤2)确定的圆盘坐标系中,根据万有引力公式,计算出检测物体对4只加速度计所产生的引力加速度信号,得到4组理想加速度信号,然后分别与实际物理环境下的加速度计固有噪声信号相加,得到4组加速度计信号。
4)将所述步骤3)中得到的4组加速度计信号,根据加速度计数学模型公式,分别计算得到4组实际加速度计输出的电流信号。
5)将所述步骤4)中得到的4组实际加速度计的电流信号经过高精度数模转换,产生4路对应实际加速度计信号的模拟电流信号,然后4路模拟电流信号分别经过电流转电压放大,产生对应的4路放大的模拟电压信号。将两两相对的加速度计的模拟电压信号经过低噪声加法放大进行信号相加操作,各自得到一路加法模拟信号。
6)将所述步骤5)得到的两路加法模拟信号经过双通道电压模数转换,使模拟信号转换为数字信号,得到的数字信号分别设为D12和D34,对数字信号D12和D34分别进行1倍频数字信号解调,各自得到一组解调结果。
7)将所述步骤6)中得到的两组解调结果进行如下分析判断:
若数字信号D12的1倍频解调结果不为零,则调整加速度计A1的标度因数,然后返回步骤4);
若数字信号D34的1倍频解调结果不为零,则调整加速度计A4的标度因数,然后返回步骤4);
若数字信号D12和D34的1倍频解调结果均为零,则进入步骤8)。
8)将所述步骤5)中输出的两路加法模拟信号经减法放大得到一路减法信号,将减法信号经过单通道电压模数转换,使模拟信号转换为数字信号,输出的数字信号记为E,对数字信号E进行频率解调,得到一组解调数据。
9)将所述步骤8)中得到的一组解调结果进行分析,当解调数据的值为零时,直接进入步骤10),否则同时调整加速度计A3、A4的标度因数,并保证二者相等后返回步骤4)。
10)对所述步骤8)中得到的数字信号E进行2倍频梯度解调,得到检测物体对圆盘中心处产生的重力梯度值。
本发明方法的优选方案中,步骤1)中进行初始化的仿真***参数为:加速度计数据文件,加速度计标度因数和零位偏置,加速度计检测质量中心到圆盘中心的距离,圆盘旋转角频率,***放大增益,初始相位角,信号采样率。
本发明是利用精密仪器实现旋转加速度计重力梯度信号仿真的方法,在精密仪器高精度基础上,通过加入各种误差源数据,分析重力梯度信号特征,找出误差源影响其信号的规律,对旋转加速度计重力梯度仪***的设计起着重要的理论和技术支撑作用。
有益效果:本发明与现有技术相比,具有以下优点:
本发明是利用精密仪表模拟重力梯度信号仿真的一种方法,通过该仿真***,在理想的加速度数据中加入各种误差数据,如标度因数误差、加速度计噪声、旋转角加速度误差、加速度计安装误差等,对含有噪声的梯度信号进行分析,得出误差源影响梯度信号的规律。重力梯度信号仿真方法具有参数可调、操作便捷、针对性和实用性强等特点,是一种接近实际物理***环境的分析方法。
重力梯度信号非常微弱,在实际物理***环境中较难提取出来,通过本发明的重力梯度信号仿真方法,可以实现重力梯度信号全方位模拟。实际加速度计数据与各种误差数据相加经过高精度数模转换、加法放大和减法放大电路,输出含有重力梯度信息的重力梯度信号,然后经过高精度模数转换,将模拟信号转为数字信号,最后对该数字信号进行分析,根据误差解调结果对加速度计标度因数进行反馈调整,最终解调出高精度重力梯度信号;通过本仿真方法,可以模拟再现旋转加速度计重力梯度仪的噪声和误差机理,有助于验证电路设计的正确性,方便验证旋转加速度计重力梯度仪的数据处理方法和误差补偿算法,对旋转加速度计重力梯度仪***的设计有重要的实际价值。。
附图说明
图1为旋转加速度计重力梯度仪原理图。
图2为仿真***工作示意图。
图3为梯度信号动态解调示意图。
图4为仿真***工作流程图。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步地说明。
旋转加速度计重力梯度仪测量原理如图1所示,四个加速度计相互正交的安装在旋转圆盘上,加速度计敏感轴正切于圆盘边沿,相对的两只加速度计敏感轴方向相反。图2是仿真***工作示意图,图3是梯度信号动态解调示意图,图4是仿真***工作流程图,具体流程步骤如下:
1)计算机对所有参数进行初始化,参数主要有:需要进行测试仿真的加速度计数据文件、加速度计标度因数和零位偏置、加速度计检测质量中心到圆盘中心的距离、圆盘旋转角频率、***放大增益、初始相位角、信号采样率等,设加速度计A1、A2、A3、A4标度因数分别设为K11、K12、K13、K14,零位偏置分别设为K01、K02、K03、K04,加速度计检测质量中心到圆盘中心的距离设为R,圆盘旋转角频率设为ω,***放大增益设为K,初始相位角设为θ,信号采样率设为fs,以上参数的初始化根据测试条件和实验要求进行设定,其中四只加速度计的标度因数设置为不一致,零位偏置设置后,保证一级放大器不饱和,***放大增益为仿真***真实信号放大倍数。
2)将四只高精度加速度计均匀正交放置在圆盘上,加速度计敏感轴正切于圆盘边沿,相对的两只加速度计的敏感轴方向相反,相邻两只加速度计的敏感轴方向相互垂直,以圆盘中心为原点,建立圆盘坐标系,定为东北天地理坐标系,记为OXYZ,坐标系原点为O,圆盘的旋转轴垂直于圆盘平面。
3)将待检测物体置于步骤2)所设的坐标系中,根据万有引力公式,计算出检测物体对4只加速度计所产生的引力加速度信号,得到4组理想加速度信号,然后分别与固定频率信号、实际物理环境下加速度计固有噪声信号相加,其中固定频率信号设为δsinωst,δ为固定频率信号幅值,ωs为固定频率信号频率。
4)将所述步骤3)中得到4组加速度计信号,根据加速度计数学模型公式,分别计算得到4组实际加速度计输出的电流信号。
5)计算机将所述步骤4)中得到的4组实际加速度电流信号发出,送给高精度数模转换器,产生4路对应实际加速度信号的模拟电流信号,输出的4路模拟电流信号分别设为A1、A2、A3和A4,其中信号A1模拟的是加速度计A1实际输出的电流信号,其他信号,以此类推。然后4路模拟电流信号A1、A2、A3和A4分别经过电流转电压放大器,产生对应的4路放大的模拟电压信号,输出的4路模拟电压信号分别设为V1、V2、V3和V4。信号V1和V2经过加法放大器1#,实现信号相加操作,加法放大器1#的输出设为V12,信号V3和V4经过加法放大器2#,实现信号相加操作,加法放大器2#的输出设为V34。
6)将所述步骤5)得到的两路加法模拟信号V12和V34经过双通道电压模数转换器,使模拟信号转换为数字信号,得到的数字信号分别设为D12和D34,计算机读取D12和D34数据,然后对D12和D34进行1倍频数字信号解调,分别得到两组解调结果。
7)将所述步骤6)中得到的两组解调结果分别进行分析,若信号D12的1倍频解调结果不为零,则调整加速度计A1标度因数K11;若信号D34的1倍频解调结果不为零,则调整加速度计A4标度因数K14,然后返回步骤4),否则进入步骤8),直到所述步骤6)中得到的两组解调结果为零为止。
8)将所述步骤5)中输出的两路加法模拟信号V12和V34经减法放大器得到一路减法信号,减法放大的输出信号经过单通道电压模数转换器,使模拟信号转换为数字信号,输出的数字信号设为E,计算机读取E数据,然后对数字信号E在频率ωs处进行固定频率解调,得到一组解调数据。
9)将所述步骤8)中得到的一组解调结果进行分析,当解调结果的值不为零时,同时调整加速度计A3、A4的标度因数K13和K14,并保证K13=K14,然后返回步骤4),否则进入步骤10),直到所述步骤8)中得到的解调结果为零为止。
10)对所述步骤8)中得到的数字信号E进行2倍频梯度解调,得到检测物体对圆盘中心处产生的重力梯度值。
11)计算机对步骤10)中得到的重力梯度值进行波形显示和数据保存。
利用本本发明方法的仿真***,可以实现各种误差源分析,验证电路设计的正确性,实现误差补偿研究,并可实现数据处理方法研究。
应理解上述实施例仅用于说明本发明技术方案的具体实施方式,而不用于限制本发明的范围。在阅读了本发明之后,本领域技术人员对本发明的各种等同形式的修改和替换均落于本申请权利要求所限定的保护范围。
Claims (2)
1.一种旋转加速度计重力梯度仪重力梯度信号仿真方法,其特征在于,该方法包括以下步骤:
1)对仿真***参数进行初始化;
2)将四只高精度加速度计A1,A2,A3和A4均匀正交放置在旋转圆盘上,A1与A2正对,A3与A4正对,相对两个加速度计的输入轴方向相反,相邻两个加速度计的输入轴相互垂直,且四个加速度计的输入轴方向正切于圆盘边沿,将圆盘坐标系记为OXYZ,原点为O,圆盘圆心位于坐标系原点,圆盘旋转轴垂直于圆盘平面;
3)将待检测物体置于所述步骤2)确定的圆盘坐标系中,根据万有引力公式,计算出检测物体对4只加速度计所产生的引力加速度信号,得到4组理想加速度信号,然后分别与实际物理环境下的加速度计固有噪声信号相加,得到4组加速度计信号;
4)将所述步骤3)中得到的4组加速度计信号,根据加速度计数学模型公式,分别计算得到4组实际加速度计输出的电流信号;
5)将所述步骤4)中得到的4组实际加速度计的电流信号经过高精度数模转换,产生4路对应所述加速度计信号的模拟电流信号,然后所述4路模拟电流信号分别经过电流转电压放大,产生对应的4路放大的模拟电压信号;将两两相对的加速度计的模拟电压信号经过低噪声加法放大进行信号相加操作,各自得到一路加法模拟信号;
6)将所述步骤5)得到的两路加法模拟信号经过双通道电压模数转换,使模拟信号转换为数字信号,得到的数字信号分别设为D12和D34,对所述数字信号D12和D34分别进行1倍频数字信号解调,各自得到一组解调结果;
7)将所述步骤6)中得到的两组解调结果进行如下分析判断:
若数字信号D12的1倍频解调结果不为零,则调整加速度计A1的标度因数,然后返回步骤4);
若数字信号D34的1倍频解调结果不为零,则调整加速度计A4的标度因数,然后返回步骤4);
若数字信号D12和D34的1倍频解调结果均为零,则进入步骤8);
8)将所述步骤5)中输出的两路加法模拟信号经减法放大得到一路减法信号,将所述减法信号经过单通道电压模数转换,使模拟信号转换为数字信号,输出的数字信号记为E,对数字信号E进行频率解调,得到一组解调数据;
9)将所述步骤8)中得到的一组解调结果进行分析,当解调数据的值为零时,直接进入步骤10),否则同时调整加速度计A3、A4的标度因数,并使二者相等后返回步骤4);
10)对所述步骤8)中得到的数字信号E进行2倍频梯度解调,得到检测物体对圆盘中心处产生的重力梯度值。
2.根据权利要求1所述的旋转加速度计重力梯度仪重力梯度信号仿真方法,其特征在于,所述步骤1)中进行初始化的仿真***参数为:加速度计数据文件,加速度计标度因数和零位偏置,加速度计检测质量中心到圆盘中心的距离,圆盘旋转角频率,***放大增益,初始相位角,信号采样率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410610089.3A CN104459826B (zh) | 2014-11-03 | 2014-11-03 | 旋转加速度计重力梯度仪重力梯度信号仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410610089.3A CN104459826B (zh) | 2014-11-03 | 2014-11-03 | 旋转加速度计重力梯度仪重力梯度信号仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104459826A CN104459826A (zh) | 2015-03-25 |
CN104459826B true CN104459826B (zh) | 2017-09-26 |
Family
ID=52906170
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410610089.3A Active CN104459826B (zh) | 2014-11-03 | 2014-11-03 | 旋转加速度计重力梯度仪重力梯度信号仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104459826B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044795A (zh) * | 2015-05-29 | 2015-11-11 | 中国船舶重工集团公司第七0七研究所 | 利用光纤陀螺减小旋转加速度计重力梯度仪测量误差的方法 |
CN104898176B (zh) * | 2015-06-10 | 2017-10-20 | 东南大学 | 一种旋转加速度计重力梯度仪重力梯度解调方法 |
CN105044798A (zh) * | 2015-06-29 | 2015-11-11 | 东南大学 | 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法 |
CN105717553B (zh) * | 2016-01-29 | 2018-01-16 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
CN106353705B (zh) * | 2016-09-23 | 2019-02-05 | 沈阳东软医疗***有限公司 | 一种梯度信号输出的方法及装置 |
CN107643548B (zh) * | 2017-09-04 | 2019-03-12 | 东南大学 | 一种多频率重力梯度激励信号产生方法 |
CN109085654B (zh) * | 2018-06-11 | 2020-01-07 | 东南大学 | 一种旋转加速度计重力梯度仪数字建模仿真方法 |
CN110927822B (zh) * | 2019-12-03 | 2021-05-18 | 吉林大学 | 一种评估希尔伯特变换算法获得磁梯度张量准确性的方法 |
CN113447994B (zh) * | 2021-06-25 | 2022-04-12 | 中国船舶重工集团公司第七0七研究所 | 一种基于旋转加速度计原理的重力梯度全张量测量方法 |
CN113885099B (zh) * | 2021-09-28 | 2024-02-27 | 中国船舶重工集团公司第七0七研究所 | 重力梯度仪加速度计标度因数不一致性动态实时估计方法 |
CN114089441B (zh) * | 2021-11-26 | 2022-08-30 | 华中科技大学 | 一种重力梯度仪测量***数值仿真方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5357802A (en) * | 1993-05-07 | 1994-10-25 | Textron, Incorporated | Rotating accelerometer gradiometer |
US5922951A (en) * | 1997-06-11 | 1999-07-13 | The Broken Hill Proprietary Company Ltd. | Gravity gradiometer |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007038820A1 (en) * | 2005-10-06 | 2007-04-12 | Technological Resources Pty Limited | Gravity gradiometer |
CA2690893A1 (en) * | 2006-11-20 | 2008-05-29 | Technological Resources Pty Limited | A gravity gradiometer |
CN102043068A (zh) * | 2010-11-05 | 2011-05-04 | 中国船舶重工集团公司第七○七研究所 | 一种标度因数在线可调的高分辨率加速度计 |
CN102053282A (zh) * | 2010-11-09 | 2011-05-11 | 北京奥地探测仪器公司 | 井中三分量磁力梯度测量装置 |
CN103472493B (zh) * | 2013-09-06 | 2015-11-25 | 北京航天控制仪器研究所 | 一种用于重力梯度测量的微弱信号检测电路 |
-
2014
- 2014-11-03 CN CN201410610089.3A patent/CN104459826B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5357802A (en) * | 1993-05-07 | 1994-10-25 | Textron, Incorporated | Rotating accelerometer gradiometer |
US5922951A (en) * | 1997-06-11 | 1999-07-13 | The Broken Hill Proprietary Company Ltd. | Gravity gradiometer |
Non-Patent Citations (2)
Title |
---|
旋转加速度计重力梯度仪误差分析;李海兵;《中国惯性技术学报》;20091031;第17卷(第5期);第1-4部分,图1 * |
旋转加速度计重力梯度仪重力梯度信号仿真;蔡体菁;《物探与化探》;20151231;第39卷;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104459826A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104459826B (zh) | 旋转加速度计重力梯度仪重力梯度信号仿真方法 | |
Munschy et al. | Scalar, vector, tensor magnetic anomalies: Measurement or computation? | |
CN106997035A (zh) | 一种基于磁梯度不变量的磁梯度计校正方法 | |
CN102494699B (zh) | 捷联式航空重力仪测量参数置信度评估方法 | |
US8359920B2 (en) | Gravity sensing instrument | |
CN103499845B (zh) | 一种利用重力梯度仪测量引力梯度的方法 | |
CN104898176A (zh) | 一种旋转加速度计重力梯度仪重力梯度解调方法 | |
CN113885098B (zh) | 一种重力敏感器低频频率响应误差在线建模及补偿方法 | |
van Leeuwen | BHP develops airborne gravity gradiometer for mineral exploration | |
CN106772488A (zh) | 一种航空重力测量gps后处理方法 | |
Bian et al. | Experimental aeromagnetic survey using a rotary-wing aircraft system: A case study in Heizhugou, Sichuan, China | |
Hammer et al. | Exploring for stratigraphic traps with gravity gradients | |
Yin et al. | Calculation of the magnetic gradient tensor from total magnetic anomaly field based on regularized method in frequency domain | |
Schiffler et al. | Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data | |
Han et al. | A predictive model for the streamwise velocity in the near‐neutral atmospheric surface layer | |
CN109085654A (zh) | 一种旋转加速度计重力梯度仪数字建模仿真方法 | |
CN108931824A (zh) | 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法 | |
CN102621506B (zh) | 岩矿石标本磁参数的总场磁力仪4方位测定方法 | |
Zhang et al. | A component compensation method for magnetic interferential field | |
Georgievich et al. | Multielement triaxial fluxgate gradiometer with variable base | |
Sui et al. | Error analysis and correction of a downhole rotating magnetic full-tensor gradiometer | |
Sonnenfeld et al. | Comparing E field changes aloft to lightning mapping data | |
Yu et al. | Generation of a multifrequency gravitational-gradient excitation and extraction of the gravitational gradient with a rotating accelerometer gravity gradiometer | |
Rahimi et al. | Improving the calibration process of inertial measurement unit for marine applications | |
Díaz Michelena et al. | Anisotropic MagnetoResistance (AMR) instrument to study the Martian magnetic environment from the surface: Expected scientific return |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |