大地坐标系应变张量表达及其与地心直角坐标系的相互转换

推导并给出了不同正交曲线坐标系应变张量和转动张量的普适表达和在旋转椭球坐标系下的应变张量表达,给出了不同正交曲线坐标系之间应变张量转换的普适表达式,以及在大地坐标系与地心直角坐标系这两种坐标系之间应变张量矩阵相互转换的具体表达,可供实际研究工作应用。     

曲线坐标系下的块体运动与应变模型研究

本文建立了顾及地球扁率和局部切标架随点变化特性的椭球坐标系下的刚体运动模型和块体运动与应变模型,以及球坐标系下顾及局部切标架随点变化特性的严密的块体运动与应变模型,分析了球坐标系下块体运动与应变模型及椭球坐标系下的块体运动与应变模型间的差异;通过计算具体讨论了地球扁率和曲线坐标系的局部切标架随点变化特性对欧拉矢量与应变张量的影响.结果表明：地球扁率对刚体欧拉矢量和应变参数的影响甚小,具体计算时可以不予考虑,但曲线坐标系的局部切标架随点变化特性对两者的影响较大,在建模过程中需要顾及,常用的Savage模型需要修正.     

利用GPS地形变资料在大地坐标系中计算应变

根据正交曲线坐标中应变张量的定义，得到了用大地坐标的椭球面上计算应变的方法。由此可直接利用ＧＰＳ或其他３维大地测量技术得到的位移或其他形变量计算地应变。这一方法克服了利用平面坐标作计算的限制，无投影误差影响。本文给出了用首都圈ＧＰＳ地形变监测网１９９５￣１９９６年椭球面上的位移求得的应变量。     

球坐标系下三维大地电磁正演研究

大地电磁正演理论研究热点一直以来主要集中在如何提高计算效率和精度,但在剖面足够长、探测深度足够大的情况下,传统的笛卡尔坐标系数值模拟方式难以准确拟合地球曲率形态.本文研究了基于球坐标系的三维大地电磁正演,推导了交错网格有限差分三维正演公式,与一维解析解和三维标准模型测试对比,验证了正演算法的正确性.通过理论模型计算,对比分析球坐标和笛卡尔坐标系正演结果表明：球坐标系模拟更合理,避免了传统笛卡尔坐标拉伸投影所引入的误差,可代替目前的笛卡尔坐标模拟方法.基于球坐标和笛卡尔坐标系的三维大地电磁正演响应值随着频率变低差异越明显.球坐标和笛卡尔坐标计算结果差异度与频率、模型结构和电阻率有关.本文模型计算结果在数万秒周期处已出现接近10%的差异,对于较大尺度的长周期大地电磁,地球曲率的影响不能忽略.     

由多个震源机制解分析甘肃及边邻地区应力场特征

基于震源断层面解的空间取向和断层滑动方向,由相应力轴张量在地理坐标系中的表达式,通过坐标旋转、代数平均可求得多个震源机制解的P、B、T轴参数的平均应力场。据甘肃及边邻地区的分区活动特征,将研究地区分为8个小区域,由各个区域内多个地震震源机制解,采用力轴张量计算法,得到了各个区应力张量的定量分析结果。     

球坐标系中图形单元应变与旋转张量及其误差解算

利用泰勒级数展开和弹性力学的几何方程,推导了球坐标系中由GPS位移数据解算图形单元应变和旋转张量的解析公式. 通过线性化处理,并利用误差传播定律,详细推导了图形单元应变和旋转张量的误差公式. 利用川滇地区最新的GPS测站位移速率数据,采用图形单元法解算了该区的面应变率及最大剪切应变率分布,并对其进行了初步的分析. 同时阐述了图形单元应变数学模型的局限性,讨论了图形单元法计算应变的意义以及图形单元的选择问题,进一步分析了应变计算的定权方法,还讨论了GPS观测网图形单元的地心半径差异与应变的关系问题.      

应变固体潮不变量理论值计算及其应用

由球坐标系下的应变固体潮理论值所构成应变张量矩阵推导并给出了应变固体潮不变量理论表达式,并采用摩尔圆定理证明了公式的正确性。在此基础上,利用应变不变量与坐标系选择无关的性质,又利用与目前所使用的钻孔应变观测系统的4条测线呈均匀分布的特点,推导并给出了5组计算不变量的观测组合公式,并以应变固体潮理论值取代观测值,按观测组合公式计算出了5组应变固体潮不变量理论值。结果发现,不变量的5组数值几乎相同,仅在小数点后2位有差异,证明了公式的正确性与不变量的唯一性。然后根据由莫尔圆定理给出的5组观测组合公式对漳州钻孔4分量应变观测资料进行了不变量计算,并对计算结果进行了相关与回归分析。最后,对不变量理论值在钻孔应变观测数据处理中的应用作了较详细的介绍,考虑到摩尔圆定理不适用于钻孔应变,又给出了采用加衬模型计算不变量的方法与观测组合公式,并对有关问题进行了讨论。     

全张量重力梯度数据滤波处理(英文)

在石油和金属矿勘探中,相对于重力数据,重力梯度张量数据含有高频的信号成分,能更好的描述小的异常特征。然而,全张量重力梯度仪测量值中含有高频随机噪声。从高频信号成分中分离出噪声将是处理重力梯度张量数据的一个挑战。本文在拉普拉斯方程约束条件下推导了重力梯度张量的笛卡尔方程和位场的表达式,然后应用笛卡尔方程通过最优线性反演方法拟合测量的重力梯度张量值。从而去除测量值中的噪声。通过模型实验,证明了这种方法不仅能很好的去除高频的随机噪声,而且能增强被噪声淹没的弱异常信号。与传统的低通滤波方法相比,避免了通过牺牲分辨率来达到去除噪声的缺点。最后将该方法应用到Bell Geospace在Vinton Dome测得的Air-FTG梯度张量数据中,并取得了很好的效果。     

关于手持GPS测定活断层位置中求解坐标系统转换参数的问题

测定活断层的坐标、绘制活断层空间位置分布图,是城市活断层探测的一个重要目的。手持GPS可以方便、快捷的测定目标点的WGS-84坐标,测量精度也满足1∶1万地形图的要求。本文介绍如何求解系统转换参数,将WGS-84坐标转换为我国目前使用的是1954年北京坐标系或1980年国家大地坐标系。进而可以利用GIS技术建立活断层信息管理数据库系统。     

利用多个震源机制解求东大别地区平均应力场

利用地震的震源机制解资料,得到应力张量在地理坐标系下的表达式,进而可计算平均应力张量.通过求平均应力张量的本征向量,即可得到其主轴方向,并由此推断区域应力场方向.利用东大别地区219次地震的震源机制解资料,比较了不同相似程度和不同起算震级资料得到的主轴方向,认为该方法计算结果非常稳定,主轴方位角的误差小于5°,倾角的误...     

Expression of strain and rotation tensor in geodetic coordinates

In this paper, unit moving trihedron is first constructed for a point on the surface of a revolution ellipsoid. Via translation, the origin of the trihedron coincides with that of Cartesian coordinates established at the center of the ellipsoid, and then through two coordinate rotations, the trihedron completely coincides with the Cartesian coordi-nates. Transformation formulae between the moving trihedron and unit Cartesian coordinate frameworks as well as transformation of point displacement between two unit coordinate frameworks are presented. Based on the above transformation formulae between two different coordinate frameworks, due to the fact that the displacement and moving trihedron of the point are both functions of the geodetic coordinates, components in the corresponding axis for differential of displacement vector and geodetic curves arc differential at the point in geodetic system can be obtained through complicated derivation. Displacement gradient matrix at the point in geodetic system is also given. Finally, expressions of strain and rotation tensor in geodetic coordinates are presented. Geometric meanings of the rotation tensor are explained in detail. The intrinsic relationship between strain tensors of sphere and ellip-soid are also discussed.     

An optimized method to transform the Cartesian to geodetic coordinates on a triaxial ellipsoid

Studia Geophysica et Geodaetica - A general triaxial ellipsoid is suitable to represent the reference surface of the celestial bodies. The transformation from the Cartesian to geodetic coordinates...     

Increasing numerical efficiency of iterative solution for total least-squares in datum transformations

Cartesian coordinate transformation between two erroneous coordinate systems is considered within the Errors-In-Variables (EIV) model. The adjustment of this model is usually called the total Least-Squares (LS). There are many iterative algorithms given in geodetic literature for this adjustment. They give equivalent results for the same example and for the same user-defined convergence error tolerance. However, their convergence speed and stability are affected adversely if the coefficient matrix of the normal equations in the iterative solution is ill-conditioned. The well-known numerical techniques, such as regularization, shifting-scaling of the variables in the model, etc., for fixing this problem are not applied easily to the complicated equations of these algorithms. The EIV model for coordinate transformations can be considered as the nonlinear Gauss-Helmert (GH) model. The (weighted) standard LS adjustment of the iteratively linearized GH model yields the (weighted) total LS solution. It is uncomplicated to use the above-mentioned numerical techniques in this LS adjustment procedure. In this contribution, it is shown how properly diminished coordinate systems can be used in the iterative solution of this adjustment. Although its equations are mainly studied herein for 3D similarity transformation with differential rotations, they can be derived for other kinds of coordinate transformations as shown in the study. The convergence properties of the algorithms established based on the LS adjustment of the GH model are studied considering numerical examples. These examples show that using the diminished coordinates for both systems increases the numerical efficiency of the iterative solution for total LS in geodetic datum transformation: the corresponding algorithm working with the diminished coordinates converges much faster with an error of at least 10^-5 times smaller than the one working with the original coordinates.     

Two modified algorithms to transform Cartesian to geodetic coordinates on a triaxial ellipsoid

The paper presents modified Lin and Wang??s (1995) and Hedgley??s (1976) algorithms to tackle the problem of transforming Cartesian to geodetic coordinates on a triaxial ellipsoid. Originally, the methods were developed for an ellipsoid of revolution but due to their universality, they may be adapted to the more complicated problem stated on a triaxial ellipsoid what is in fact done in this work. Two modified methods are compared to the vector method recently introduced by Feltens. The modified methods turn out to be more accurate and faster than the algorithm presented by Feltens.     

利用地球重力位模型计算重力和重力梯度

高阶高精度地球重力场模型具有广泛的用途。本文利用地球重力位模型计算重力和重力梯度在应用中很有实用阶值,同时也是计算重力场其它量的关键。利用伪局部笛卡尔坐标与球坐标的关系计算了重力与重力梯度在伪局部笛卡尔坐标系下的分量;利用张量变换的原理给出了已知重力与重力梯度在某一坐标系下的分量求它们在另一坐标系下分量的方法,并具体给出了重力与重力梯度在局部笛卡儿坐标系下的分量计算公式,同时还给出计算重力场五参量与垂线偏差的计算公式,本研究推进了地球重力场的可视化进程。     

椭球坐标系下多震相地震射线追踪

地震射线追踪方法技术在地震学领域有着较为广泛的应用,然而大多数算法建立在直角坐标系或球坐标系下,实际地球并非完美的球体,而是两极略扁的椭球体,因此,球坐标系下计算结果与真实情况存在一定误差.传统的做法一般是在球坐标系下进行计算,而后进行椭球校正.本文提出了一种直接在椭球体模型中采用分区多步最短路径算法进行多震相地震射线追踪的方法技术,实现了椭球坐标系下多震相地震波射线路径追踪和走时计算.与解析解的对比表明:该算法具有较高的计算精度,适用于任意形状的椭球体,且不需要进行额外的走时校正.数值模拟结果表明,计算所得P波和PcP反射波的走时与AK135走时表的误差小于0.1 s.当震中距较大时,使用球对称模型和椭球体模型计算所得的走时差异显著,说明采用椭球坐标系的必要性.     

Closed-form transformation between geodetic and ellipsoidal coordinates

We present formulas for direct closed-form transformation between geodetic coordinates (φ, λ, h) and ellipsoidal coordinates (β, λ, u) for any oblate ellipsoid of revolution. These will be useful for those dealing with ellipsoidal representations of the Earth’s gravity field or other oblate ellipsoidal figures. The numerical stability of the transformations for near-polar and near-equatorial regions is also considered.     

Computation of strain and rotation tensor as well as their uncertainties for small arrays in spherical coordinate system

Based on Taylor series expansion and strain components expressions of elastic mechanics, we derive formulae of strain and rotation tensor for small arrays in spherical coordinates system. By linearization process of the formulae, we also derive expressions of strain components and Euler vector uncertainties respectively for subnets using the law of error propagation. Taking GPS velocity field in Sichuan-Yunnan area as an example, we compute dilation rate and maximum shear strain rate field using the above procedure, and their characteristics are preliminarily car- ried on. Limits of the strain model for small array are also discussed. We make detailed explanations on small array method and the choice of small arrays. How to set weights of GPS observations are further discussed. Moreover relationship between strain and radius of GPS subnets is also analyzed.