共查询到20条相似文献,搜索用时 328 毫秒
1.
2.
谱元法作为一种基于变分原理的数值方法,已经成为地球动力学数值模拟的重要手段.本文旨在对Martinec(2000)正演模拟表面负荷问题的谱-有限元方法进行改进.简单地介绍了谱元法计算负荷边值问题的相关理论,并推导了中性分层液体中的平衡方程以及其对应的弱形式解.在地球半径方向上,采取高阶样条函数作为负荷解的基函数,并利用高斯-勒让德方法进行数值积分.均质地球模型数值试验表明,相较于Martinec的方法,本文方法可以有效地加快数值解的收敛速度.最后,利用PREM模型计算了SNREI地球的负荷勒夫数.数值计算结果表明本文方法的收敛性和稳定性都非常好.本文方法获得的负荷勒夫数与采用传统龙格-库塔积分方法获得的结果相对误差在0.01%的量级. 相似文献
3.
本文给出了单位点质量负荷作用在球状成层地球模型上的解。对较新的G-D1066A地球模型求出了n直到10000阶的负荷勒夫数,并利用这组负荷勒夫数计算了格林函数的值。简述了如何利用格林函数计算地球对表面任何负荷的响应。讨论了负荷潮研究在地球物理和海洋学方面的一些应用。 相似文献
4.
5.
6.
基于PREM模型,利用非自转、球型分层、各向同性、理想弹性(SNREI)地球的形变理论,讨论了地球在不同驱动力作用下的形变特征.采用地球位移场方程的4阶Runge Kutta数值积分方法,解算了在表面负荷和日月引潮力作用下地球表面和内部形变和扰动位,并给出了地球表面的负荷Love数和体潮Love数.结果表明在固体内核中的形变很小,液核中低阶(n<10)负荷位移随半径的变化非常复杂.当负荷阶数超过10时,地核中的形变和扰动位都很小,地球的响应主要表现为弹性地幔中的径向位移,且随深度增加急剧减弱,负荷阶数越高这种衰减的速度越快.SNREI地球的地表负荷Love数和体潮Love数与信号频率的依赖关系很弱.在计算体潮Love数的过程中,采用了SNREI地球的运动方程,同时考虑了由于地球自转和椭率引起的核幔边界附加压力,这一近似处理方法获得的结果能很好地符合地球表面重力潮汐实际观测结果. 相似文献
7.
研究固体潮应力场时,需要考虑负荷潮应力场的影响,特别是在沿海地区,在地球的表层负荷潮应力场的大小甚至接近或超过固体潮应力场.根据Alterman方法,本文给出计算球状径向不均匀弹性地球模型对点源负荷的应力响应(负荷潮应力格林函数)的方法;作为计算实例,计算出GB地球模型的负荷潮应力格林函数.只要给出全球海潮模型,在球面上完成海潮潮高与负荷潮应力格林函数的褶积,即可计算出海潮负荷在地球模型表面上和内部任意点产生的负荷潮应力. 相似文献
8.
9.
10.
本文设计了三个数值算例分别模拟径向分层黏弹地球对二阶球谐负荷、均厚圆盘负荷和全球末次冰期真实冰负荷的响应.利用耦合位扰动Laplace方程的三维有限元法(CLFE),计算出地球表面的位移、大地水准面的变化、冰后海平面变化和这些量的现今变化速率,并将其结果与传统谱方法的结果进行对比.结果表明,两种方法的结果总体上表现出良好的一致性,因此,CLFE方法被证实是正确的、有效的、能被用于研究地幔横向非均匀性的影响.值得注意的是,新的方法完全没有横向非均匀小扰动假设的限制. 相似文献
11.
本文研究用理论地震图反演地震震源参数的方法,利用16个WWSSN台记录的渤海地震(1969年7月18日,M=7.4)远震P波波形,用我们的反演方法重新测定了该地震的震源参数。文中以剪切位错源的理论地震图与实际观测记录波形的相关程度为判据,测得该地震的震源参数为:地震矩3.9×1027达因·厘米;震源破裂持续时间6秒;位错面两个可能的选择解是:(1)θs=207°,δ=87°NW,λ=-159°,右旋;(2)θs=298.6°,δ=69°NE,λ=-3.6°,左旋。其中θS、δ、λ分别是震源位错面的方位角、倾角和错滑角;震源深度为25公里。数值运算的结果还表明,本文提供的测定震源参数的方法实际效能良好。 文中还与用P波初动符号方法测定的震源参数结果进行了对比分析。 相似文献
12.
13.
14.
15.
16.
频域激发极化法中体极化与面极化内在联系的探讨 总被引:2,自引:0,他引:2
矿化的体极化介质,其激发极化效应来自其中所含电子导电矿物颗粒的表面极化效应。基于上述认识,本文从理论上推导出含球状电子导体的体极化介质复电阻率的“微观”数学表达式。将它与体极化介质复电阻率的“宏观”表达式--Cole-Cole模型相比较,得出了:1.Cole-Cole模型中的特征参量与球状电子导体的表面极化系数、含量以及导体半径等的综合关系式;2.电子导体表面极化系数在频域中的一般表达式。这表达式同目前在实验室中观测到的结果是一致的。 相似文献
17.
共炮检距二维弹性波地震剖面的偏移方法 总被引:1,自引:0,他引:1
在均匀各向同性介质的分界面上入射纵波时,除产生纵波反射P-P波外,还会产生转换横波的反射P-SV波。同样,当入射波为SV横波时,在分界面上会同时产生SV-SV反射波和SV-P反射转换波,对这种二维弹性波地震记录如何进行偏移是正在探讨的问题。为了能够处理实际地震数据,本文提出在共炮检距剖面上进行二维弹性波记录的偏移方法。具体做法是,首先用弹性波方程将两分量的地震记录分解为P-P(或SV-SV)波和P-SV(或SV-P)波,然后用我们已导出的适合于处理共炮检距剖面的叠前偏移方法对它们进行处理。用模拟数据进行的试验说明方法正确,效果良好。 相似文献
18.
19.
本文利用东亚及其邻近地区50个地磁台站1955-1977年地磁三要素(D、H、Z)的实测年平均值算出的非偶极场值X、Y、Z,选用“按距离加权最小二乘曲面拟合”方法,计算、绘制、分析了东亚地区地磁非偶极场及其长期变化的时空分布特征,并与国际参考场(IGRF)模式进行了粗略的对比。所得的主要结果是:1.1955.0-1965.0年代,东亚地区Z分量的增长速度在迅速减小;1965.0年代以后,Z值在逐年下降。2.从整体上看,东亚地区地磁非偶极场虽然形态与IGRF给出的结果基本吻合,但长期变化有明显差异。1965.0年代的IGRF长期变化模式与实测结果偏差甚大。为了使IGRF模式更好地拟合地磁长期变化处于转折时期的复杂地区的情况,看来时间间隔应更短些,例如5年为好,采用现在通用的方法是不适宜的。 相似文献
20.
本文从麦氏方程出发,运用电流磁场的等效理论阐明了在磁电勘探中,当矿体的二次极化电流具有垂直对称轴分布时,在地表面及以上空间中任一点的磁异常场(H)均等于零。因而,具有任何垂直对称状电化学矿产模式(石油、天然气、地下水及金属矿床等)的磁电勘探法异常场线积分均不存在。此外,由于对各种形状和任何极化倾角的矿体情况,地表面上均无电流密度垂直分量(in)存在,故地表面上异常场线积分均为零(∮L H·dL=0)。因此,磁电勘探法不具备在地面上寻找此种矿产的物理前提。 本文还对倾斜极化球体和柱体的磁电勘探异常进行了计算,给出了主剖面上理论曲线的异常分布规律,指出了磁电异常的数量级。同时,还对井中磁电异常的分布规律给出了计算资料。为了对比,也对相应条件下的自然电场异常做了计算,指出了电场异常与磁场异常之间的不同特点及其对实际找矿的意义。 相似文献