共查询到18条相似文献，搜索用时 343 毫秒
1.
位场向下与向上延拓之间存在固有的内在联系,向上延拓解算具有稳定可靠的优良特性,本文据此提出了借助向上延拓信息实现航空重力向下延拓稳定解算的两种方法,分别建立了点对点向下解析延拓模型和最小二乘向下解析延拓模型.其核心思想是,依据泰勒级数展开模型,将位场向下延拓解算过程转换为向上延拓计算和垂向偏导数解算两个步骤,通过第一步的处理有效抑制数据观测噪声对解算结果的干扰,通过第二步的处理成功实现向下延拓反问题的稳定解算,较好地解决了向下延拓解算固有的不适定性问题.分析研究了两种解析延拓模型的计算精度及适用条件,利用超高阶位模型EGM2008建立的模拟标准场数据对两种模型解算结果的合理性和有效性进行了数值验证,证明本文新方法实用易行,具有较高的应用价值. 相似文献
2.
重力场向下延拓能够突出局部和浅部的异常信息,分离叠加的异常特征.但是向下延拓通常具有过程不稳定、下延深度小、结果不准确等问题.针对向下延拓所存在的不足,本文利用重力场及其垂向一阶导数,基于辛普森(Simpson)求积公式,推导出重力场向下延拓米尔尼(Milne)公式.将本文向下延拓方法应用于模型数据,向下延拓模型结果及误差曲线表明,相对于向下延拓快速傅里叶变换(FFT)法和积分迭代法,向下延拓Milne法的深度更大,相对误差更小;相对模型值,向下延拓Milne法能够获得稳定且准确的结果.对加拿大乃查科(Nechako)盆地地区实测航空重力数据进行本文方法向下延拓验证,处理结果表明,相对于实测异常,本文方法向下延拓结果能够很好还原实测数据,并且在进一步向下延拓中反映原始异常的趋势,增强局部和细小异常信息. 相似文献
3.
建立高精度高分辨率区域地球重力场模型是大地测量学科的重要科学目标之一,可为固体地球物理、地球动力学、地震学以及资源勘探等领域提供基础信息,也是当前大地测量现代化发展实现GNSS测定海拔高程的重大需求.提出了一种基于矩谐分析的区域重力场建模方法,在局部直角坐标系中求解地球引力位的Laplace方程,推导了扰动位、重力异常、重力扰动、大地水准面差距和垂线偏差的矩谐展开式以及矩谐系数阶方差和阶误差的计算公式,给出了利用重力观测值进行矩谐分析求解矩谐系数的数学模型和算法,为抑制矩谐分析中存在的周期延拓边界效应,提出了扩展计算区域范围的策略.利用EGM2008重力位模型生成模拟地面重力和航空重力观测值,加入标准差为2 mGal的高斯白噪声,分别设计基于地面重力数据和基于航空重力数据的模拟数值试验,对矩谐分析建模方法进行了验证与分析.数值结果表明基于矩谐分析构建的区域重力场模型可靠且可达很高的精度,由地面和航空重力数据计算的2.5′×2.5′大地水准面精度分别达到1 cm和1.4 cm,4 km飞行高度处航空重力观测值的向下延拓误差仅为3.1mGal,矩谐分析可为区域重力场精细结构的逼近提供一种新的选择. 相似文献
4.
针对位场向下延拓的不适定性,我们将位场向下延拓视为向上延拓的反问题,提出以位场最小曲率作为约束条件来求解稳定的下延位场.我们将剖面位场向上延拓表达式用傅里叶矩阵的形式表示,以矩阵乘法形式给出延拓的表达式,同时向待反演的下延位场引入最小曲率约束,得到向下延拓的最小曲率解,并利用正交变换给出了更为简洁的频率域解.随后,利用Kronecker积将上述全部结果拓展至三维位场,给出了三维位场向下延拓的最小曲率解.此外,我们将位场数据的填充、扩充问题与向下延拓问题统筹考虑,提出一种新的向下延拓迭代格式,该算法面向实际资料处理需求、无须预扩充或填补数据.下延迭代时,对原始数据直接向下延拓,而空白部分利用上一次下延位场估计的上延值替代其空白值并对其向下延拓,直至获得最小曲率约束下稳定的向下延拓结果.同时,我们也讨论了利用改进L曲线和广义交叉验证(GCV)计算正则参数最优估计的问题.对理论模型和实际航空重力资料进行了向下延拓检验,处理结果表明位场向下延拓的最小曲率方法解能满足实际位场资料对向下延拓处理的需求,具有较高的下延精度. 相似文献
5.
顾及远区影响的向下延拓实用算法 总被引：2，自引：0，他引：2
向下延拓是Stokes边值解算中必不可少的环节.为了研究向下延拓的远区影响,本文在泊松实用公式的基础上给出了顾及远区影响的向下延拓算法的推导过程.对4°×4°实验区内5'×5'分辨率模型扰动重力的向下延拓结果表明:对于地面重力数据的向下延拓,积分半径为1°、2°时,顾及远区影响的算法可将精度分别提高0.36 m Gal、0.17 m Gal;对于4000 m航高的重力数据的向下延拓,积分半径为1°、2°时,该算法可将精度分别提高1.40 m Gal、0.67 m Gal.延拓高度越高,积分半径越小,顾及远区影响的算法对精度的改善效果越明显. 相似文献
6.
在大地坐标系下利用Forsberg局部扰动位协方差模型（即DPM模型）导出了实用的局部重力异常协方差模型（即GAC模型）和局部扰动重力协方差模型（即GDC模型）,两种模型形式完全一致.针对GAC模型,提出了两种模型参数的拟合方法,即按照泊松积分向上延拓获得的不同高度数据进行拟合以及按照测量区域的平面实测数据进行拟合,通过某地区的实测数据检验得出两种参数拟合方法得到的参数值相差不大,这种差别在向上和向下延拓过程中的影响可以忽略.依据本文的算例,GAC模型作为配置法的协方差模型用于延拓时,其向上延拓的精度在1.8×10－5 m·s-2左右,向下延拓的精度在5×10－5 m·s-2左右,完全满足局部重力场在中等山区的延拓要求.通过对不同高度下GAC模型用于延拓效果的对比,可以得出基于GAC模型的延拓精度随着高度的增加而衰减,在满足测量精度要求下其最大向下延拓高度约为7 km.总体而言,本文推导的GAC模型能够很好地利用地形数据,较好地满足了航空重力测量在局部重力场的延拓要求. 相似文献
7.
8.
联合不同类型重力测量数据确定地球重力场模型的迭代法 总被引：1，自引：0，他引：1
不同的重力测量数据包含了不同波段的地球重力场信息,因此要恢复更高精度的地球重力场模型,就必须对不同类型的重力测量数据进行联合处理.以地面重力异常Δg为例,推导了利用迭代法联合不同类型重力测量数据反演地球重力场模型的基本原理公式,并给出了其具体实现步骤,接着采用全球的重力异常Δg数据和扰动位T数据,基于迭代法对卫星重力梯度SGG数据解算的重力场模型进行了进一步的精化.结果表明,初始的卫星重力梯度SGG模型和经过全球重力异常Δg数据精化后的模型,其对应的累计大地水准面误差分别达到1.128cm和0.048cm、累计重力异常误差分别达到0.416mGal和0.018mGal的精度;在经过全球扰动位T数据进一步精化后的模型,其对应的累计大地水准面误差达到0.043cm、累计重力异常误差达到0.016mGal的精度. 相似文献
9.
我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性. 相似文献
10.
将水平观测面上的实测位场值，垂直投影至下部的延拓水平面上，作为该水平面上的位场初始值. 根据该水平面上的初始值，用快速傅里叶变换(FFT)的方法向上延拓计算观测面上的位场值. 用观测面上的实测值与计算值的差值，对延拓面上的位场值进行校正. 如此反复迭代，直至观测面上的实测值与计算值的差值小到可以忽略. 这种空间域的迭代法原理简单，不用解线性代数方程组，有较高的计算速度和良好的延拓效果. 本文用迭代法对模型数据和实际数据进行向下延拓，对比了迭代法与常规的FFT法在位场向下延拓中的效果，迭代法显著优于FFT法. 相似文献
11.
全张量测量技术是在空中或海上用加载了多个加速度计的移动平台技术测量位场的五个独立分量.各张量分量包含不同方向的地下地质体信息,水平张量分量T_(xx)、T_(yy)、T_(xy)、T_(xz)、T_(yz)通常用于识别和映射与地质构造或地层变化有关的测量区域中的目标,垂直张量分量Tzz用于估计深度.然而,这些分量传统上是彼此分开解释,经常遇到错失关键信息的风险.本文所用全张量欧拉反褶积是在单独z方向的欧拉反演基础上发展而来的,它融合了重力异常垂直分量以及其三个方向导数、水平分量以及其三个方向导数.全张量数据信息得以有效应用的同时,欧拉反褶积结果也比常规欧拉反褶积结果更加收敛.最后,结合美国墨西哥湾地区实测航空FTG数据,用重力梯度张量数据进行联合欧拉三维反演研究,有效的识别岩盖的边界信息,划分岩盖范围,为进一步研究盖层底下深部复杂地质情况提供可靠的解释结果. 相似文献
12.
JIANG Tao LI JianCheng DANG YaMin ZHANG ChuanYin WANG ZhengTao KE BaoGui 《中国科学:地球科学(英文版)》2014,57(7):1637-1644
Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy,and can provide fundamental information for geophysics,geodynamics,seismology,and mineral exploration.Rectangular harmonic analysis(RHA)is proposed for regional gravity field modeling in this paper.By solving the Laplace’s equation of gravitational potential in local Cartesian coordinate system,the rectangular harmonic expansions of disturbing potential,gravity anomaly,gravity disturbance,geoid undulation and deflection of the vertical are derived,and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients(RHC).We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations.In order to reduce the edge effects caused by periodic continuation in RHA,we propose the strategy of extending the size of computation domain.The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal.The accuracy of the 2.5′×2.5′geoid undulations computed from ground and airborne gravity data is 1 and 1.4cm,respectively.The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 mGal.Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model,which may be a new option for the representation of the fine structure of regional gravity field. 相似文献
13.
捷联式航空重力测量系统与平台式系统相比具有体积小、重量轻、功耗低等许多优点,近些年来取得了显著的研究进展.本文给出了捷联式航空重力测量的两种算法模型:捷联式惯性标量重力测量(SISG)和旋转不变式标量重力测量(RISG)模型,并对其误差模型作了初步讨论.利用我国首套捷联式航空重力仪SGA-WZ01在某海域的部分试验数据,对两种算法模型进行了比较分析,表明其差值之标准差对于200s的滤波长度小于0.5mGal.同时,利用两组重复测线数据估算了不同滤波尺度下的两种算法的内符合精度,表明SISG算法略优于RISG算法.对于200s和300s的滤波长度,SISG的内符合精度分别为1.06mGal和0.80mGal. 相似文献
14.
融合多源数据的高精度、高分辨率的局部重力场建模是物理大地测量学的前沿和热点问题.本文研究了基于径向基函数融合多源数据的局部重力场建模方法,利用Monte-Carlo方差分量估计实现了不同类型的观测数据的合理定权,引入了最小标准差法确定基函数的适宜网络,分析了地形因素对于基函数网络确定及局部重力场建模精度的影响.以泊松小波基函数为构造基函数,结合残差地形模型,融合实测的陆地重力异常、船载重力异常及航空重力扰动数据构建了局部区域陆海统一的似大地水准面模型.研究结果表明:引入残差地形模型平滑了地形质量引入的高频扰动信号,简化了基函数的网络设计;并提高了重力似大地水准面的精度,平原地区其精度提高了4mm,地形起伏较大的山区其精度提高了约5cm.总体而言,基于"三步法"构建的局部重力似大地水准面在荷兰、比利时及德国相关区域,其精度分别达到1.12cm、2.80cm以及2.92cm. 相似文献
15.
Hasan Yildiz Rene Forsberg Carl Christian Tscherning Daniel Steinhage Graeme Eagles Johannes Bouman 《Studia Geophysica et Geodaetica》2017,61(1):53-68
An airborne gravity campaign was carried out at the Dome-C survey area in East Antarctica between the 17th and 22nd of January 2013, in order to provide data for an experiment to validate GOCE satellite gravity gradients. After typical filtering for airborne gravity data, the cross-over error statistics for the few crossing points are 11.3 mGal root mean square (rms) error, corresponding to an rms line error of 8.0 mGal. This number is relatively large due to the rough flight conditions, short lines and field handling procedures used. Comparison of the airborne gravity data with GOCE RL4 spherical harmonic models confirmed the quality of the airborne data and that they contain more high-frequency signal than the global models. First, the airborne gravity data were upward continued to GOCE altitude to predict gravity gradients in the local North-East-Up reference frame. In this step, the least squares collocation using the ITGGRACE2010S field to degree and order 90 as reference field, which is subtracted from both the airborne gravity and GOCE gravity gradients, was applied. Then, the predicted gradients were rotated to the gradiometer reference frame using level 1 attitude quaternion data. The validation with the airborne gravity data was limited to the accurate gradient anomalies (TXX, TYY, TZZ and TXZ) where the long-wavelength information of the GOCE gradients has been replaced with GOCO03s signal to avoid contamination with GOCE gradient errors at these wavelengths. The comparison shows standard deviations between the predicted and GOCE gradient anomalies TXX, TYY, TZZ and TXZ of 9.9, 11.5, 11.6 and 10.4 mE, respectively. A more precise airborne gravity survey of the southern polar gap which is not observed by GOCE would thus provide gradient predictions at a better accuracy, complementing the GOCE coverage in this region. 相似文献
16.
Determination of gravity anomaly at sea level from inversion of satellite gravity gradiometric data 总被引：1，自引：0，他引：1
Gravity gradients can be used to determine the local gravity field of the Earth. This paper investigates downward continuation of all elements of the disturbing gravitational tensor at satellite level using the second-order partial derivatives of the extended Stokes formula in the local-north oriented frame to determine the gravity anomaly at sea level. It considers the inversion of each gradient separately as well as their joint inversion. Numerical studies show that the gradients Tzz, Txx, Tyy and Txz have similar capability of being continued downward to sea level in the presence of white noise, while the gradient Tyz is considerably worse than the others. The bias-corrected joint inversion process shows the possibility of recovering the gravity anomaly with 1 mGal accuracy. Variance component estimation is also tested to update the observation weights in the joint inversion. 相似文献
17.
本文在法方程层面融合GOCE卫星的Vxx、Vyy、Vzz和Vxz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型. 相似文献
18.
Numerical Studies on the Harmonic Downward Continuation of Band-Limited Airborne Gravity 总被引：1，自引：0，他引：1
In this paper, the numerical stability and efficiency of methods of harmonic downward continuation from flying altitudes are treated for sampled gravity field data. The problem is first formulated in its continuous form, i.e. as the inverse solution of the spherical Dirichlet problem, and is then approximated by Gaussian quadrature to yield a finite system of linear equations. The numerical stability of this system is investigated for both error-free gravity data and for the noisy and band-limited gravity measurements usually obtained from airborne gravity surveys. It can be shown that the system becomes ill-conditioned, once the ratio between flying altitude and data sampling rate exceeds a certain limit. It can also be shown that noisy measurements tend to generate a solution that is practically useless, long before the system becomes ill-conditioned. Therefore, instead of treating the general solution of the discrete downward continuation problem, the more modest question is studied, for which range of flying altitudes and sampling rates, the numerical solution of the discrete linear system can be considered as practically useful. Practically useful will be defined heuristically as of sufficient accuracy and stability to satisfy the requirements of the user. The question will be investigated for the specific application of geoid computation from gravity data sampled at flying altitudes. In this case, a stable solution with a standard deviation of a few centimeters is required. Typical flight parameters are heights of 2–6 km, a minimum half-wavelength resolution of 2 km, and data noise between 0.5 and 1.5 mGal. Different methods of geoid determination, different solution techniques for the resulting systems of linear equations, and different minimization principles will be compared. As a result operational parameters will be defined which, for a given noise level, will result in a geoid accuracy of a few centimeters for the estimated band-limited gravity field spectrum. 相似文献