首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 328 毫秒
1.
SNRVEI地球模型连续分布的简正模及其意义   总被引:5,自引:1,他引:5       下载免费PDF全文
基于Peltier等建立的粘弹地球表面负荷形变理论,运用简正模分析法,实现了SNRVEI粘弹地球模型表面负荷勒夫数从拉氏域至时间域的变换.数值结果印证了韩大仲关于真实地球存在“连续分布的简正模”的结论,同时还表明在保证较高计算精度前提下,为便于计算,将地球核慢边界以上部分划分成较少个物性均匀层是可行的.本文的方法和结果经过数值检验,可用于冰后相对海平面变化、地球自转非潮汐加速度变化等地表负荷均衡动态响应观测量的研究.  相似文献   

2.
利用谱元法计算SNREI地球的表面负荷变形   总被引:1,自引:0,他引:1       下载免费PDF全文
谱元法作为一种基于变分原理的数值方法,已经成为地球动力学数值模拟的重要手段.本文旨在对Martinec(2000)正演模拟表面负荷问题的谱-有限元方法进行改进.简单地介绍了谱元法计算负荷边值问题的相关理论,并推导了中性分层液体中的平衡方程以及其对应的弱形式解.在地球半径方向上,采取高阶样条函数作为负荷解的基函数,并利用高斯-勒让德方法进行数值积分.均质地球模型数值试验表明,相较于Martinec的方法,本文方法可以有效地加快数值解的收敛速度.最后,利用PREM模型计算了SNREI地球的负荷勒夫数.数值计算结果表明本文方法的收敛性和稳定性都非常好.本文方法获得的负荷勒夫数与采用传统龙格-库塔积分方法获得的结果相对误差在0.01%的量级.  相似文献   

3.
地球表面质量负荷的静态响应   总被引:7,自引:2,他引:7       下载免费PDF全文
本文给出了单位点质量负荷作用在球状成层地球模型上的解。对较新的G-D1066A地球模型求出了n直到10000阶的负荷勒夫数,并利用这组负荷勒夫数计算了格林函数的值。简述了如何利用格林函数计算地球对表面任何负荷的响应。讨论了负荷潮研究在地球物理和海洋学方面的一些应用。  相似文献   

4.
时萍  孙文科 《地震学报》1988,10(2):183-191
本文较系统地研究和讨论了海潮负荷效应计算的精度问题。建立了总的及单项的误差估算模型。对中国大陆M2波重力、倾斜和应变的负荷计算都给出了精度估值。结果表明,海潮负荷效应计算精度明显地受地球模型和海潮图精度的限制,尤其是地球模型选取不同,对结果影响更大。越靠近沿海,其相对精度越低。并对目前的负荷研究、改正精度等进行了分析和评价。   相似文献   

5.
根据卫-卫跟踪观测技术的测量原理,基于能量守恒法建立了一种新的双星相互跟踪和三星相互跟踪的卫星观测方程. 通过数值模拟,采用预处理共轭梯度法恢复120阶地球重力场. 模拟结果表明:第一,双星相互跟踪恢复地球重力场的精度和美国喷气动力实验室公布的EIGEN GRACE02S的结果相符合;第二,三星相互跟踪恢复地球重力场的精度较双星提高约2倍.  相似文献   

6.
SNREI地球对表面负荷和引潮力的形变响应   总被引:5,自引:2,他引:5       下载免费PDF全文
基于PREM模型,利用非自转、球型分层、各向同性、理想弹性(SNREI)地球的形变理论,讨论了地球在不同驱动力作用下的形变特征.采用地球位移场方程的4阶Runge Kutta数值积分方法,解算了在表面负荷和日月引潮力作用下地球表面和内部形变和扰动位,并给出了地球表面的负荷Love数和体潮Love数.结果表明在固体内核中的形变很小,液核中低阶(n<10)负荷位移随半径的变化非常复杂.当负荷阶数超过10时,地核中的形变和扰动位都很小,地球的响应主要表现为弹性地幔中的径向位移,且随深度增加急剧减弱,负荷阶数越高这种衰减的速度越快.SNREI地球的地表负荷Love数和体潮Love数与信号频率的依赖关系很弱.在计算体潮Love数的过程中,采用了SNREI地球的运动方程,同时考虑了由于地球自转和椭率引起的核幔边界附加压力,这一近似处理方法获得的结果能很好地符合地球表面重力潮汐实际观测结果.  相似文献   

7.
研究固体潮应力场时,需要考虑负荷潮应力场的影响,特别是在沿海地区,在地球的表层负荷潮应力场的大小甚至接近或超过固体潮应力场.根据Alterman方法,本文给出计算球状径向不均匀弹性地球模型对点源负荷的应力响应(负荷潮应力格林函数)的方法;作为计算实例,计算出GB地球模型的负荷潮应力格林函数.只要给出全球海潮模型,在球面上完成海潮潮高与负荷潮应力格林函数的褶积,即可计算出海潮负荷在地球模型表面上和内部任意点产生的负荷潮应力.  相似文献   

8.
不同地球模型对负荷潮汐改正的影响   总被引:5,自引:0,他引:5       下载免费PDF全文
负荷潮汐改正与地壳上地幔结构有密切关系,主要原因是负荷潮汐取决于高阶勒夫数。利用较新的1066A地球模型,对负荷勒夫数进行了重新推算,并依据负荷勒夫数组成重力、倾斜、应变及位移格林函数。发现在近区用1066A模型和经典的G-B模型求出的格林函数值差别较大。据此,推估出不同地球模型对负荷潮汐改正的影响,其中远区估算按解析方法进行,近区由数值法进行。  相似文献   

9.
弹性波边界元法正演模拟   总被引:10,自引:4,他引:10       下载免费PDF全文
弹性波边界元地震模型方法(BEESM),实现了二维和三维问题的纵、横波及转换波的同时模拟,并且能模拟任意复杂构造的地震声波正演模型.根据地震模型的特点,本文发展了数值积分计算与矩阵消元同步进行的块状高斯消元法;用解析法处理奇异积分;用无限元法处理边界吸收问题;采用单元长度随介质速度和计算频率变化的变单元算法,及自动剖分单元等技术,提高了计算精度,节省了内存,缩短了计算时间.  相似文献   

10.
横向非均匀地球负荷问题的CLFE有限元算法的有效性   总被引:1,自引:2,他引:1       下载免费PDF全文
本文设计了三个数值算例分别模拟径向分层黏弹地球对二阶球谐负荷、均厚圆盘负荷和全球末次冰期真实冰负荷的响应.利用耦合位扰动Laplace方程的三维有限元法(CLFE),计算出地球表面的位移、大地水准面的变化、冰后海平面变化和这些量的现今变化速率,并将其结果与传统谱方法的结果进行对比.结果表明,两种方法的结果总体上表现出良好的一致性,因此,CLFE方法被证实是正确的、有效的、能被用于研究地幔横向非均匀性的影响.值得注意的是,新的方法完全没有横向非均匀小扰动假设的限制.  相似文献   

11.
利用远震P波波形反演渤海地震的震源参数   总被引:6,自引:1,他引:6       下载免费PDF全文
本文研究用理论地震图反演地震震源参数的方法,利用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.
信号的瞬时参数与正交小波基   总被引:12,自引:3,他引:9       下载免费PDF全文
对于可表成u(t)=a(t)cosθ(t)的信号,称a(t)为其瞬时振幅,θ(t)为瞬时位相,为瞬时频率.在许多实际问题中,瞬时振幅、位相、频率等瞬时参数的提取有重要意义.但已知的一些瞬时参数计算方法均各有其局限性.本文构造了一种频谱支集紧且有解析表达式的小波,利用它得到一种计算信号瞬时参数的能抗噪声干扰的快速、局部算法.理论分析和数值模拟的结果表明,该方法是实用而有效的,可望能满足许多实际问题的需要.  相似文献   

13.
为模拟震源的动力学破裂过程,本文讨论一种求解剪切裂纹动力学扩展问题的有限差分法。在研究二维反平面破裂的基础上,我们把所用方法推广到三维问题。研究了许多简单震源模型,通过对比同一问题的解析解和其它数值解,对方法的正确性和计算精度进行了检验,结果表明我们的数值方法是可行的。最后,作为初步应用的例子,我们研究了非均匀断层的自发破裂问题。  相似文献   

14.
根据起伏地形上观测的异常场分量(沿任一方向的分量)换算出垂直与水平分量,本文提供了二维情况下的转换方法。其算法是:(1)通过共形映照,将w平面上的地形曲线L的上部映照为z平面的上半平面,将曲线L映照成z平面上的实轴。(2)求解边界为实轴的Hilbert边值问题。 通过理论模型的计算,验证了该方法的正确性。  相似文献   

15.
利用重力资料反演三维密度界面的质面系数法   总被引:9,自引:2,他引:9       下载免费PDF全文
文中给出一深度为h的水平矩形质面在地面产生的重力公式,对比了此公式算出的重力异常与三维长方体公式算出的重力异常,得出:当三维地质体的厚度与其埋深之比小于0.2-0.4时,二者的重力效应几乎相等,它们之间最大偏差在1-3%以内。用给出的重力公式进一步导出了重力反演三维密度界面的质面系数法。用理想模式检验,结果表明,本文给出的方法收敛快,反演误差较小,计算效果较好。最后用一实例,验证了方法的适用性。  相似文献   

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.
三分量地震波场的近似解析离散模拟技术   总被引:3,自引:0,他引:3       下载免费PDF全文
首先将二维各向异性介质中地震波动二阶偏微分方程组降为关于时间t的一组一阶偏微分波动方程组,然后利用Taylor展式和插值逼近的方法构制了一种与有限元法、反射率法、射线追踪法、有限差分法等传统方法不同的近似解析离散方法(NADM).对双层各向同性介质以及横向各向同性介质中的三分量地震波场进行了模拟,并与传统的有限差分法作了比较,结果表明NADM算法可行.  相似文献   

19.
东亚地区地磁非偶极场长期变化的分析   总被引:1,自引:3,他引:1       下载免费PDF全文
本文利用东亚及其邻近地区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)。因此,磁电勘探法不具备在地面上寻找此种矿产的物理前提。 本文还对倾斜极化球体和柱体的磁电勘探异常进行了计算,给出了主剖面上理论曲线的异常分布规律,指出了磁电异常的数量级。同时,还对井中磁电异常的分布规律给出了计算资料。为了对比,也对相应条件下的自然电场异常做了计算,指出了电场异常与磁场异常之间的不同特点及其对实际找矿的意义。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号