CN108376187A - 一种海域流动点外部扰动引力垂向分量的无奇异计算方法 - Google Patents
一种海域流动点外部扰动引力垂向分量的无奇异计算方法 Download PDFInfo
- Publication number
- CN108376187A CN108376187A CN201810054341.5A CN201810054341A CN108376187A CN 108376187 A CN108376187 A CN 108376187A CN 201810054341 A CN201810054341 A CN 201810054341A CN 108376187 A CN108376187 A CN 108376187A
- Authority
- CN
- China
- Prior art keywords
- external disturbance
- vertical component
- marine site
- disturbance gravitation
- point
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种海域流动点外部扰动引力垂向分量的无奇异计算方法,其主要技术特点是:利用局部积分域恒等式变换实现海域格网点外部扰动引力垂向分量的无奇异计算;联合局部积分域恒等式变换和内插方法实现海域非格网点外部扰动引力垂向分量的无奇异计算。本发明针对表层法计算海域流动点外部扰动引力垂向分量存在奇异性的固有缺陷,通过引入局部积分域恒等式变换和非网格点内插方法,建立海域流动点外部扰动引力垂向分量的无奇异计算模型,能够有效获取海域流动点外部扰动引力垂向分量,通过实验验证,该方法实现了海域流动点外部扰动引力垂向分量的高精度计算,满足全海域和全高度段对外部扰动引力垂向分量快速赋值的实际需求。
Description
技术领域
本发明属于海洋重力技术领域,尤其是一种海域流动点外部扰动引力垂向分量的无奇异计算方法。
背景技术
海域外部扰动引力垂向分量计算在大地测量和空间科学技术研究、人造卫星及各类航天器飞行保障中具有非常重要的应用价值。表层法模型具有核函数的数学结构简单,输入数据保障要求容易得到满足等特点,在计算海域流动点外部扰动引力垂向分量中得到广泛应用,但是,表层法存在超低空计算的积分奇异性问题。
发明内容
本发明的目的在于克服现有技术的不足,提出一种海域流动点外部扰动引力垂向分量的无奇异计算方法,解决基于表层法计算海域流动点外部扰动引力垂向分量存在的奇异性问题。
本发明解决其技术问题是采取以下技术方案实现的:
一种海域流动点外部扰动引力垂向分量的无奇异计算方法,包括以下步骤:
步骤1、利用局部积分域恒等式变换实现海域格网点外部扰动引力垂向分量的无奇异计算;
步骤2、联合局部积分域恒等式变换和内插方法实现海域非格网点外部扰动引力垂向分量的无奇异计算。
所述步骤1局部积分域恒等式变换公式为:
式中,Tp为外部空间计算点P处的扰动位;ψ为计算点和流动点之间的球面角距;l为计算点和流动点之间的空间距离;γ为海面平均重力;μ称为广义面密度;R为地球椭球平均半径;ψ0为数据区块的积分半径。
步骤2的实现方法为:当计算点为非网格点时,利用周围的4个网格点计算结果内插出计算点上的扰动引力,设计算点距4个网格点的距离分别为:s1、s2、s3和s4,由局部积分域恒等式变换公式计算得到4个网格点外部扰动引力垂向分量影响量分别为:δg1、δg2、δg3和δg4,然后由下式内插得到计算点扰动引力分量的对应部分:
本发明的优点和积极效果是:
本发明针对表层法计算海域流动点外部扰动引力垂向分量存在奇异性的固有缺陷,通过引入局部积分域恒等式变换和非网格点内插方法,建立海域流动点外部扰动引力垂向分量的无奇异计算模型,能够有效获取海域流动点外部扰动引力垂向分量,通过实验验证,该方法实现了海域流动点外部扰动引力垂向分量的高精度计算,满足全海域和全高度段对外部扰动引力垂向分量快速赋值的实际需求。
具体实施方式
以下通过实施例对本发明做进一步详述。
一种海域流动点外部扰动引力垂向分量的无奇异计算方法,包括以下步骤:
步骤1、利用局部积分域恒等式变换实现海域格网点外部扰动引力垂向分量的无奇异计算。具体实现方法为:
利用表层法计算地球外部扰动引力垂向分量的积分式为:
其中:
r=R+h
式中,Tp为外部空间计算点P处的扰动位;和分别代表计算点和球面流动点的地心向径、地心纬度和地心经度;ψ为计算点和流动点之间的球面角距;l为计算点和流动点之间的空间距离;Δg为海面观测重力异常;T为海面点扰动位;ζ为海面高程异常;γ为海面平均重力;μ称为广义面密度;R为地球椭球平均半径;h为计算高度(大地高);dσ为单位球积分面积元。
由公式(1)得知,当计算高度h→0和球面角距ψ→0,即计算点与数据点重合时,公式(1)的数值积分将出现不同程度的奇异性问题。为了解决公式(1)的积分奇异性问题,对公式(1)进行局部恒等积分式变换,得到如下公式:
假设该数据区块的积分半径为ψ0,公式(2)可换为:
将公式(3)右侧第二项积分成解析式,得到局部积分域恒等变换式:
不难看出,公式(4)不再存在积分奇异性问题。
步骤2、联合局部积分域恒等式变换和内插方法实现海域非格网点外部扰动引力垂向分量的无奇异计算。具体实现过程如下:
公式(4)是以数据网格点作为计算点为前提条件建立起来的。当计算点为非网格点时,恒等变换公式(4)要求的积分区域对称性假设不再严格成立,故它们不能直接用于非数据网格点的计算。此时,可利用计算点(非网格点)周围的4个网格点计算结果内插出计算点上的扰动引力。设计算点距4个网格点的距离分别为:s1、s2、s3和s4,由公式(4)给出的局部积分域恒等式变换计算得到4个网格点外部扰动引力垂向分量影响量分别为:δg1、δg2、δg3和δg4,则可由下式内插得到计算点扰动引力分量的对应部分:
下面结合实例说明本发明的效果。
采用2160阶次的全球位模型EGM2008作为标准场开展数值计算检验,由其模拟产生外部扰动引力三分量的基准值,同时计算1°×1°、20′×20′、5′×5′和2′×2′四组地面重力异常和高程异常的伪观测量,作为确定外部空间扰动引力场的基础数据。试验区涵盖了重力场变化比较剧烈的马里亚纳海沟,区域范围为:λ:90°E~190°E,计算某一点扰动引力采用的四种分辨率基础数据覆盖范围为:1°×1°→30°×30°;20′×20′→10°×10°;5′×5′→4°×4°;2′×2′→2°×2°;参考场的阶次取为N=36。表1列出了四组地面重力异常和高程异常计算量的统计结果。
表1 地面重力异常和高程异常统计结果
为了考察本发明提出的无奇异表层法计算模型对超低空扰动引力垂向分量的逼近效果,这里特别选取位于马里亚纳海沟的两个点作为计算点,分别采用不同的改化模型计算不同高度的扰动引力垂向分量,并将其与相对应的位模型计算基准值作对比分析。其中,一个计算点选取与2′×2′数据网格点重合,具体位置为另一个计算点选取为非数据网格点,具体位置为
表2 P1点和P2点在各高度上扰动引力垂向分量值(单位:10-5m˙s-2)
对于格网点P1,利用格网点无奇异算法(公式(4))计算外部扰动引力垂向分量,对于非格网点P2,利用格网点无奇异算法(公式(4))先计算该点周围的4个网格点外部扰动引力垂向分量,然后利用公式(5)内插出非格网点P2的外部扰动引力垂向分量。为了比对分析本发明算法的先进性,引入传统表层法(即公式(1))计算结果作为参考。表3和表4分别列出了两个计算点在不同高度上不同赋值模型计算结果与对应基准值的互比结果。
表3 P1点不同赋值模型计算结果与对应基准值的互比结果(单位:10-5m˙s-2)
高度/km | 0 | 0.02 | 1 | 3 | 5 | 10 | 30 | 50 | 100 | 200 | 300 |
传统模型 | -110.51 | -109.59 | -67.73 | -21.63 | -8.59 | -2.17 | -0.20 | -0.01 | 0.10 | 0.12 | 0.09 |
本发明 | -0.15 | -0.15 | -0.07 | -0.08 | -0.07 | -0.03 | 0.03 | 0.00 | -0.13 | -0.13 | -0.07 |
表4 P2点不同赋值模型计算结果与对应基准值的互比结果(单位:10-5m˙s-2)
高度/km | 0 | 0.02 | 1 | 3 | 5 | 10 | 30 | 50 | 100 | 200 | 300 |
传统模型 | -104.45 | -102.76 | -34.45 | -1.51 | -0.10 | -0.02 | 0.06 | 0.10 | 0.16 | 0.17 | 0.13 |
本发明 | -1.60 | -1.60 | -1.41 | -1.27 | -1.19 | -1.01 | -0.59 | -0.44 | -0.36 | -0.19 | -0.09 |
由表3和表4可知,本发明提出的基于表层法的扰动引力垂向分量无奇异计算模型,在重力场变化比较剧烈的海沟区,该模型的计算误差也不超过2mGal,表明能有效克服了传统表层法积分模型固有的奇异性问题,较好地满足了全海域和全高度段对局部扰动引力垂向分量快速赋值的实际需求,具有良好的应用前景。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (3)
1.一种海域流动点外部扰动引力垂向分量的无奇异计算方法,其特征在于包括以下步骤:
步骤1、利用局部积分域恒等式变换实现海域格网点外部扰动引力垂向分量的无奇异计算;
步骤2、联合局部积分域恒等式变换和内插方法实现海域非格网点外部扰动引力垂向分量的无奇异计算。
2.根据权利要求1所述的一种海域流动点外部扰动引力垂向分量的无奇异计算方法,其特征在于:所述步骤1局部积分域恒等式变换公式为:
式中,Tp为外部空间计算点P处的扰动位;ψ为计算点和流动点之间的球面角距;l为计算点和流动点之间的空间距离;γ为海面平均重力;μ称为广义面密度;R为地球椭球平均半径;ψ0为数据区块的积分半径。
3.根据权利要求1所述的一种海域流动点外部扰动引力垂向分量的无奇异计算方法,其特征在于:步骤2的实现方法为:当计算点为非网格点时,设计算点距4个网格点的距离分别为:s1、s2、s3和s4,由局部积分域恒等式变换公式计算得到4个网格点外部扰动引力垂向分量影响量分别为:δg1、δg2、δg3和δg4,然后由下式内插得到计算点扰动引力分量的对应部分:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810054341.5A CN108376187B (zh) | 2018-01-19 | 2018-01-19 | 一种海域流动点外部扰动引力垂向分量的无奇异计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810054341.5A CN108376187B (zh) | 2018-01-19 | 2018-01-19 | 一种海域流动点外部扰动引力垂向分量的无奇异计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108376187A true CN108376187A (zh) | 2018-08-07 |
CN108376187B CN108376187B (zh) | 2021-09-10 |
Family
ID=63016604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810054341.5A Active CN108376187B (zh) | 2018-01-19 | 2018-01-19 | 一种海域流动点外部扰动引力垂向分量的无奇异计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108376187B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965128A (zh) * | 2021-02-08 | 2021-06-15 | 中国人民解放军92859部队 | 一种无奇异性顾及局域保障条件计算外部重力异常的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120296616A1 (en) * | 2011-05-16 | 2012-11-22 | Masaya Tsunoda | Three-dimensional fluid simulation method |
CN104035138A (zh) * | 2014-04-25 | 2014-09-10 | 西安测绘研究所 | 一种全球及局部海洋扰动重力的精确快速计算方法 |
CN104750983A (zh) * | 2015-03-20 | 2015-07-01 | 中国人民解放军信息工程大学 | 一种空间分层网格扰动引力场模型构建与扰动引力快速确定方法 |
US20150219455A1 (en) * | 2012-08-08 | 2015-08-06 | Hexagon Technology Center Gmbh | Inclination sensor |
CN104834320A (zh) * | 2015-03-23 | 2015-08-12 | 中国人民解放军信息工程大学 | 一种空间分层扰动引力场网格模型快速构建方法 |
CN106646647A (zh) * | 2016-11-28 | 2017-05-10 | 中国人民解放军92859部队 | 基于带限航空矢量重力确定大地水准面的一步积分直接法 |
CN106646644A (zh) * | 2016-11-28 | 2017-05-10 | 中国人民解放军92859部队 | 基于带限航空矢量重力确定大地水准面的两步积分反解法 |
CN107480338A (zh) * | 2017-07-18 | 2017-12-15 | 哈尔滨工业大学深圳研究生院 | 基于时空积分域处理时域边界元奇异积分的方法 |
-
2018
- 2018-01-19 CN CN201810054341.5A patent/CN108376187B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120296616A1 (en) * | 2011-05-16 | 2012-11-22 | Masaya Tsunoda | Three-dimensional fluid simulation method |
US20150219455A1 (en) * | 2012-08-08 | 2015-08-06 | Hexagon Technology Center Gmbh | Inclination sensor |
CN104035138A (zh) * | 2014-04-25 | 2014-09-10 | 西安测绘研究所 | 一种全球及局部海洋扰动重力的精确快速计算方法 |
CN104750983A (zh) * | 2015-03-20 | 2015-07-01 | 中国人民解放军信息工程大学 | 一种空间分层网格扰动引力场模型构建与扰动引力快速确定方法 |
CN104834320A (zh) * | 2015-03-23 | 2015-08-12 | 中国人民解放军信息工程大学 | 一种空间分层扰动引力场网格模型快速构建方法 |
CN106646647A (zh) * | 2016-11-28 | 2017-05-10 | 中国人民解放军92859部队 | 基于带限航空矢量重力确定大地水准面的一步积分直接法 |
CN106646644A (zh) * | 2016-11-28 | 2017-05-10 | 中国人民解放军92859部队 | 基于带限航空矢量重力确定大地水准面的两步积分反解法 |
CN107480338A (zh) * | 2017-07-18 | 2017-12-15 | 哈尔滨工业大学深圳研究生院 | 基于时空积分域处理时域边界元奇异积分的方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965128A (zh) * | 2021-02-08 | 2021-06-15 | 中国人民解放军92859部队 | 一种无奇异性顾及局域保障条件计算外部重力异常的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108376187B (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108010103B (zh) | 复杂河道地形快速精细生成方法 | |
CN102722909B (zh) | 一种基于自适应尺度dem的流水线拓扑网络动态模拟方法 | |
CN106646644B (zh) | 基于带限航空矢量重力确定大地水准面的两步积分反解法 | |
CN104392113B (zh) | 一种近海海面冷空气大风风速的估算方法 | |
CN108416082A (zh) | 一种海域流动点外部扰动引力水平分量的无奇异计算方法 | |
CN115688607A (zh) | 基于多通道频谱组合的带状跨海似大地水准面精化方法 | |
CN108376187A (zh) | 一种海域流动点外部扰动引力垂向分量的无奇异计算方法 | |
Hello et al. | Combined use of sensitivity information and observations to improve meteorological forecasts: A feasibility study applied to the'Christmas storm'case | |
CN110287620B (zh) | 适用于地表观测面的球坐标系密度界面正演方法及系统 | |
Altaf et al. | Efficient identification of uncertain parameters in a large‐scale tidal model of the European continental shelf by proper orthogonal decomposition | |
CN113985490B (zh) | 利用地形和地壳密度数据进行地表重力仿真的方法及装置 | |
Zhang et al. | Watershed characteristics extraction and subsequent terrain analysis based on digital elevation model in flat region | |
CN113885101B (zh) | 一种基于椭球模型构建重力梯度基准图方法 | |
Chu et al. | Japan Sea thermohaline structure and circulation. Part II: A variational P-vector method | |
Nyoka et al. | Computation of Gravity Field Functionals with a localized level Ellipsoid | |
Yamazaki et al. | Non‐hydrostatic atmospheric cut cell model on a block‐structured mesh | |
Singh et al. | Development of geoid model-A case study on western India | |
Huang et al. | Prediction of design typhoon wind speeds and profiles using refined typhoon wind field model | |
CN106646647B (zh) | 基于带限航空矢量重力确定大地水准面的一步积分直接法 | |
CN111506871A (zh) | 基于频域特性构建海域垂线偏差模型的径向基函数格网法 | |
CN106646648B (zh) | 基于带限航空矢量重力确定大地水准面的两步积分直接法 | |
Khedr | Sea level modeling and realization of lowest astronomical tide (LAT) for Alexandria Harbor, Egypt | |
Kikuchi et al. | A study of topographic multiplier considering the effect of complex terrains and tropical cyclones | |
Sīle et al. | Applying Numerical Weather Prediction Models to the Production of New European Wind Atlas: Sensitivity studies of the wind climate to the planetary boundary layer parametrization | |
CN104112180B (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 |