CN107272084A - 一种创建局部近似匀值引力梯度张量激励的方法 - Google Patents
一种创建局部近似匀值引力梯度张量激励的方法 Download PDFInfo
- Publication number
- CN107272084A CN107272084A CN201710317828.3A CN201710317828A CN107272084A CN 107272084 A CN107272084 A CN 107272084A CN 201710317828 A CN201710317828 A CN 201710317828A CN 107272084 A CN107272084 A CN 107272084A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- munderover
- mtd
- sigma
- 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
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明涉及一种创建局部近似匀值引力梯度张量激励的方法，其技术特点是包括以下步骤：构建高张量局部近似匀值引力梯度张量激励装置；根据所需梯度激励幅值，设置高张量近似匀值引力梯度张量激励的量程、激励步长和测量尺度范围，计算引力梯度激励模型下的空间梯度张量，从而形成量程可控、局部近似匀值的高张量引力梯度激励。本发明通过设计合理，在重力梯度敏感器测量尺度范围内利用匀值引力梯度张量激励对仪器输出进行标定，形成量程可控、局部近似匀值等特性的高张量引力梯度激励，解决了现有重力梯度敏感器标定方法不具备水平重力梯度测量原理中所展现的匀值梯度场特性的问题。
Description
技术领域
本发明属于重力梯度仪技术领域，尤其是一种创建局部近似匀值引力梯度张量激励的方法。
背景技术
重力梯度仪是敏感重力场变化趋势的专用设备，为了实现对重力场变化定量测量的准确性，就必须对重力梯度仪的输出结果进行有效精确的标定。而实现精确标定最有效的方法就是营造匀值、可控的重力梯度场激励，但限于现有旋转加速度计式重力梯度仪的元件水平以及原理方案制约，要想实现对水平局部张量重力梯度仪的有效标定(理论分辨率几十E)，必须创建具备高张量激励且在一定空间尺度内呈现近似匀值特性的质量引力梯度场，以解决对于设备的标定需求。
现阶段定量标定方法主要采用二倍频分辨率等效质量引力梯度张量激励法，即根据万有引力定律，计算出随重力梯度敏感器旋转调制频率相对应各采样点的激励质量引力加速度沿各重梯加速度计敏感轴的有效分量，然后对理论计算得到的一系列数据进行旋转调制频率二倍频处的数据提取，同时对真实激励下测量得到的各采样点输出电信号进行旋转调制频率二倍频处的解调与滤波处理，最后将处理后的输出信号与计算得到的理论有效加速度信息进行比例运算，从而得到重梯加速度计在调整二倍频处的有效分辨率，间接实现了与该分辨能力等效的匀值引力梯度激励标定。但上述标定方法并没有真正营造出符合水平张量重力梯度测量原理中描述的匀值梯度场激励，只体现了集成重力梯度敏感器所用的重梯加速度计具备敏感等效匀值重力梯度激励的测量分辨能力。
发明内容
本发明的目的在于克服现有技术的不足，提供一种创建局部近似匀值引力梯度张量激励的方法，解决现有重力梯度敏感器标定方法不具备水平重力梯度测量原理中所展现的匀值梯度场特性的问题。
本发明解决现有的技术问题是采取以下技术方案实现的：
一种创建局部近似匀值引力梯度张量激励的方法，包括以下步骤：
步骤1、构建高张量局部近似匀值引力梯度张量激励装置；
步骤2、设置高张量近似匀值引力梯度张量激励的量程、激励步长和测量尺度范围，计算引力梯度激励模型下的空间梯度张量，从而形成量程可控、局部近似匀值的高张量引力梯度激励。
所述高张量局部近似匀值引力梯度张量激励装置包括直线运动滑轨、由质量激励单元构成的线型引力质量体和重力梯度敏感器，重力梯度敏感器设置在两个线型引力质量体的中心位置且两个线型引力质量体安装在直线运动滑轨上，两个线型引力质量体的质量中心平面与重力梯度敏感器的旋转中心平面重合。
所述线型引力质量体由钨钼铜合金组成。
所述测量尺度范围为重力梯度敏感器周围的方形区域。
所述步骤2的重力梯度正演递推公式为：
式中， G为万有引力常数；ρ为引力质量体质量密度；ξi和ξi'分别为两侧对称引力质量体8个顶点的x轴坐标值；ηi和ηi'分别为两侧对称引力质量体8个顶点的y轴坐标值；γi和γi'分别为两侧对称引力质量体8个顶点的z轴坐标值；x,y,z为直角坐标系中任意观测点的坐标；为两侧对称引力质量体任意顶点与任意观测点间的距离；为方向系数。
本发明的优点和积极效果是：
本发明通过分析近似匀值引力梯度场的产生机理、分布特性以及对于重力梯度仪标定的实际作用，构建高张量局部近似匀值引力梯度张量激励装置，在重力梯度敏感器测量尺度范围内利用匀值引力梯度张量激励对仪器输出进行标定，形成量程可控、局部近似匀值等特性的高张量引力梯度激励，解决了现有重力梯度敏感器标定方法不具备水平重力梯度测量原理中所展现的匀值梯度场特性的问题。
附图说明
图1为点源质量分布示意图；
图2为在图1中设置观测范围示意图；
图3为Txx、Tyy和Txx-Tyy的计算结果示意图；
图4为本发明构建的高张量局部近似匀值引力梯度张量激励装置结构图；
图5为图4的俯视图；
图6为线型对称质量分布下产生50E激励时的梯度场分布示意图；
图7为线型对称质量分布下产生100E激励时的梯度场分布示意图；
图8为线型对称质量分布下产生150E激励时的梯度场分布示意图；
其中，图4及图5中的1为直线运动滑轨、2为引力质量体(由质量激励单元组成)、3为重力梯度敏感器、4为测量尺度范围。
具体实施方式
以下结合附图对本发明实施例做进一步详述：
一种创建局部近似匀值引力梯度张量激励的方法的设计原理为：
依据点源质量对称分布产生局部匀值线引力梯度场的基础，本发明认为线型质量场平面对称分布同样可产生类似的局部平面匀值引力梯度特性。质量分布如图1所示，两条线型质量在-ξ＜x＜ξ区域产生的引力梯度张量可表示如下：
式中ξ为线型质量体上任意点的x轴坐标值；η，η'分别为两侧对称线型质量体上任意点的y轴坐标值；x为两侧对称线型质量体中间区域任意观测点的x轴坐标值；y为两侧对称线型质量体中间区域任意观测点的y轴坐标值；
对式中各张量进行计算可知，当线质量源在y坐标轴区间取(-∞，+∞)时，Txy和Tyx项均为零，而测量(标定)关心的另外三个张量为：
为分析方便，利用数值积分法对上式进行有限区间积分计算，并且选取一定的观测范围利用图示法进行观察，用以研究对称线型质量分布下的引力梯度场分布特性。假设线密度ρ取1，ξ取2，η取(-10,10)区间，观测范围取以坐标原点为中心边长0.2的正方形区域，如图2所示。Txx、Tyy和Txx-Tyy的计算结果如图3和表1所示。
表1.对称线质量分布下有限测量尺度平面内的各引力梯度张量激励
从计算结果可以看出，当质量源关于观测坐标系y轴成线型对称分布时，以观测坐标系原点为中心的一定观测区域平面内的Txx-Tyy张量近似等于Txx张量且呈现近似匀值的特性，而Tyy张量近似为0。因此，通过对线型质量源对称分布的引力梯度场定量分析，的确在对称中心局部区域可产生既具有一定张量值，又具有近似匀值特性的引力梯度场，故而为后续如何定量创建可用于重力梯度传感器标定所需的近似匀值引力梯度场奠定了基础。
下面结合图4和图5给出的具体实例对本发明进行具体说明。
首先构建高张量局部近似匀值引力梯度张量激励装置，如图4及图5所示。线型引力质量体2的质量中心平面与重力梯度敏感器3的旋转中心平面(图中的XY平面)重合，通过在直线型对称滑轨机构1上对称移动线型引力质量体2，沿线型方向增加、减少引力质量体2或沿对称移动距离方向整体增加、减少线型引力质量体2的列数，来产生具有高张量且在重力梯度敏感器3测量尺度范围4内呈现近似匀值特性的引力梯度张量激励。
然后根据所需梯度激励幅值按照重力梯度正演递推公式计算得到的引力质量体2的线型分布长度，沿直线型对称滑轨机构1移动方向的分布列数，以及与重力梯度敏感器的水平分布间距，配置相应所需的质量激励元，同时控制对称直线滑轨机构1按照得到的相应移动方向进行对称移动，移动定位精度优于1mm。具体方法为：
计算引力梯度激励模型下的空间梯度张量分布公式如下：
式中，而
根据重力梯度敏感器3的标定需求，拟定创建的高张量近似匀值引力梯度张量激励为量程150E，激励步长为50E，划定重力梯度敏感器3的测量尺度范围4为340×340(单位：mm)的方形区域，同时限定所创建的近似匀值引力梯度激励误差不大于激励张量数值的10％为边界条件，假定采用100mm×100mm×100mm的钨钼铜合金(其密度为18×103kg/m3)为引力质量体2的组成单元，利用模型正演递推计算方法搜寻符合创建条件的引力质量体2的对称分布距离、质量分布长度以及梯度场分布特性(如图6至8和表2、3所示)。
表2.线型对称质量分布下梯度仪有效测量尺度内的Txx、Tyy、Txx-yy引力梯度张量激励
表3.线型对称质量分布下梯度仪有效测量尺度内的的Txy引力梯度张量激励
最后，为验证对称引力质量激励下的局部近似匀值梯度张量特性，即验证重力梯度敏感器3的感测值对引力质量体2相对重力梯度敏感器3的水平位置距离在对称中心附近变化不敏感的特性，可以在实现对称引力质量激励时，沿直线型对称滑轨机构1，同向小范围并行移动引力质量体2，以观察重力梯度敏感器感测值的波动变化，实现对局部匀值特性的验证。
需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明包括并不限于具体实施方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。
Claims (5)
1.一种创建局部近似匀值引力梯度张量激励的方法，其特征在于包括以下步骤：
步骤1、构建高张量局部近似匀值引力梯度张量激励装置；
步骤2、设置高张量近似匀值引力梯度张量激励的量程、激励步长和测量尺度范围，计算引力梯度激励模型下的空间梯度张量，从而形成量程可控、局部近似匀值的高张量引力梯度激励。
2.根据权利要求1所述的一种创建局部近似匀值引力梯度张量激励的方法，其特征在于：所述高张量局部近似匀值引力梯度张量激励装置包括直线运动滑轨、由质量激励单元构成的线型引力质量体和重力梯度敏感器，重力梯度敏感器设置在两个线型引力质量体的中心位置且两个线型引力质量体安装在直线运动滑轨上，两个线型引力质量体的质量中心平面与重力梯度敏感器的旋转中心平面重合。
3.根据权利要求2所述的一种创建局部近似匀值引力梯度张量激励的方法，其特征在于：所述线型引力质量体由钨钼铜合金组成。
4.根据权利要求1所述的一种创建局部近似匀值引力梯度张量激励的方法，其特征在于：所述测量尺度范围为重力梯度敏感器周围的方形区域。
5.根据权利要求1所述的一种创建局部近似匀值引力梯度张量激励的方法，其特征在于：所述步骤2计算引力梯度激励模型下的空间梯度张量的公式为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>z</mi>
<mi>x</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>z</mi>
<mi>y</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<mi>G</mi>
<mi>&rho;</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>arctan</mi>
<mfrac>
<mrow>
<msubsup>
<mi>y</mi>
<mi>j</mi>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>z</mi>
<mi>k</mi>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>z</mi>
<mi>k</mi>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>j</mi>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>z</mi>
<mi>k</mi>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>arctan</mi>
<mfrac>
<mrow>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>z</mi>
<mi>k</mi>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>y</mi>
<mi>j</mi>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>j</mi>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mi>arctan</mi>
<mfrac>
<mrow>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mi>j</mi>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>z</mi>
<mi>k</mi>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中， n＝1,2；G为万有引力常数；ρ为引力质量体质量密度；ξi和ξi'分别为两侧对称引力质量体8个顶点的x轴坐标值；ηi和ηi'分别为两侧对称引力质量体8个顶点的y轴坐标值；γi和γi'分别为两侧对称引力质量体8个顶点的z轴坐标值；x,y,z为直角坐标系中任意观测点的坐标；为两侧对称引力质量体任意顶点与任意观测点间的距离；为方向系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710317828.3A CN107272084B (zh) | 2017-05-08 | 2017-05-08 | 一种创建局部近似匀值引力梯度张量激励的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710317828.3A CN107272084B (zh) | 2017-05-08 | 2017-05-08 | 一种创建局部近似匀值引力梯度张量激励的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107272084A true CN107272084A (zh) | 2017-10-20 |
| CN107272084B CN107272084B (zh) | 2019-02-12 |
Family
ID=60074119
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710317828.3A Active CN107272084B (zh) | 2017-05-08 | 2017-05-08 | 一种创建局部近似匀值引力梯度张量激励的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107272084B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109902695A (zh) * | 2019-03-01 | 2019-06-18 | 辽宁工程技术大学 | 一种面向像对直线特征匹配的线特征矫正与提纯方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002103399A1 (en) * | 2001-06-15 | 2002-12-27 | Bhp Billiton Innovation Pty Ltd | Airborne gravity gradiometers |
| US20060156810A1 (en) * | 2005-01-04 | 2006-07-20 | Bell Geospace Inc. | Accelerometer and rate sensor package for gravity gradiometer instruments |
| CN104597520A (zh) * | 2015-01-05 | 2015-05-06 | 中国船舶重工集团公司第七0七研究所 | 一种重力梯度仪重力梯度敏感器输入激励产生方法 |
| CN105717553A (zh) * | 2016-01-29 | 2016-06-29 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
-
2017
- 2017-05-08 CN CN201710317828.3A patent/CN107272084B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002103399A1 (en) * | 2001-06-15 | 2002-12-27 | Bhp Billiton Innovation Pty Ltd | Airborne gravity gradiometers |
| US20060156810A1 (en) * | 2005-01-04 | 2006-07-20 | Bell Geospace Inc. | Accelerometer and rate sensor package for gravity gradiometer instruments |
| CN104597520A (zh) * | 2015-01-05 | 2015-05-06 | 中国船舶重工集团公司第七0七研究所 | 一种重力梯度仪重力梯度敏感器输入激励产生方法 |
| CN105717553A (zh) * | 2016-01-29 | 2016-06-29 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109902695A (zh) * | 2019-03-01 | 2019-06-18 | 辽宁工程技术大学 | 一种面向像对直线特征匹配的线特征矫正与提纯方法 |
| CN109902695B (zh) * | 2019-03-01 | 2022-12-20 | 辽宁工程技术大学 | 一种面向像对直线特征匹配的线特征矫正与提纯方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107272084B (zh) | 2019-02-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105717553B (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
| CN107024198B (zh) | 一种捷联惯导式测波方法及*** | |
| CN106323334B (zh) | 一种基于粒子群优化的磁力计校准方法 | |
| CN103389114B (zh) | 基于g敏感性的三轴MEMS陀螺非正交误差标定方法 | |
| US20200073010A1 (en) | Gravity gradient measurement method and apparatus | |
| CN104459826B (zh) | 旋转加速度计重力梯度仪重力梯度信号仿真方法 | |
| CN103499845B (zh) | 一种利用重力梯度仪测量引力梯度的方法 | |
| CN102628249A (zh) | 全自动惯性传感检测路面平整度装置及检测方法 | |
| CN109709628B (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
| CN106840100A (zh) | 一种数字式倾角传感器及测量方法 | |
| CN109856689B (zh) | 一种超导航磁梯度张量数据抑噪处理方法和*** | |
| CN108458714A (zh) | 一种姿态检测***中不含重力加速度的欧拉角求解方法 | |
| CN106872730A (zh) | 一种加速度测量装置及方法 | |
| CN102175391A (zh) | 红外制导导弹导引头重心位置的测量装置及其测量方法 | |
| CN102878989A (zh) | 一种利用星载线加速度计进行三轴角振动测量的方法 | |
| CN104793257B (zh) | 一种基于高压悬浮的便携式相对重力仪 | |
| CN101975582A (zh) | 一种人员运动轨迹跟踪方法和装置 | |
| CN107272084A (zh) | 一种创建局部近似匀值引力梯度张量激励的方法 | |
| CN106383366A (zh) | 一种旋转簧片式重力梯度仪 | |
| CN103091723A (zh) | 降低重力卫星质心调整误差对地球重力场精度影响的方法 | |
| CN107991049A (zh) | 基于加速度传感器的六自由度振动测试方法及装置 | |
| CN104807443B (zh) | 水平仪 | |
| CN103712634B (zh) | 光纤陀螺振动-磁场交叉耦合度的测量方法 | |
| Xing et al. | A new and high-precision gravity base network in the south of the Tibetan Plateau | |
| CN103969615B (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 |