首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
核-幔边界的动力学背景   总被引:1,自引:1,他引:1  
根据传播矩阵方法,并把由联合反演得到的同时满足长波地形起伏、板块运动速度、重力位异常资料以及地震层析先验知识的全地幔三维异常密度作为载荷,以求取核-幔边界的动力学背景.计算结果显示:1.所求的核-幔边界起伏图像与Hager等根据格林函数方法所求得的核-幔边界的起伏在全球范围内基本相符.2.核-幔边界处的环型场仅在数量上降为地表处的环型场的1/8左右,而极型场较地表处的极型场的流动图像有显著变化,数值也增为地表处极型场的3倍左右.  相似文献   

2.
根据传播矩阵方法,并把由联合反演得到的同时满足长波地形起伏、板块运动速度、重力位异常资料以及地震层析先验知识的全地幔三维异常密度作为载荷,以求取核-幔边界的动力学背景.计算结果显示:1.所求的核-幔边界起伏图像与Hager等根据格林函数方法所求得的核-幔边界的起伏在全球范围内基本相符.2.核-幔边界处的环型场仅在数量上降为地表处的环型场的1/8左右,而极型场较地表处的极型场的流动图像有显著变化,数值也增为地表处极型场的3倍左右.  相似文献   

3.
非均匀介质热蠕变流动的数值求解   总被引:5,自引:1,他引:4       下载免费PDF全文
针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格-粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov Galerkin方法,求解能量方程时采用了流线迎风Petrov Galerkin方法,网格-粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格-粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性.  相似文献   

4.
在前人工作的基础上对分形插值方法作了详细的探讨,给出了分形插值函数的显式表达方式.在量纲分析的基础上给出了垂直比例因子的局部显式表达式,旨在提高地震道插值重建的精度及突出局部信息,并从单道地震图的角度分析其在地震道插值重建中的应用效果.研究了垂直比例因子的变化对分形插值精度的影响.数值实验表明,随着垂直比例因子的增大,分形垂直的误差逐渐增大,二者之间呈显出指数增长的趋势.该法克服了随机分形插值方法必须进行多步迭代的弱点,提高了计算效率.通过对理论地震道插值重建的分析,说明了本文分形插值方法的高精度和高效率.本文提出的显式分形插值方法既能够突出地震道数据的局部信息,又较好地保持了地震道数据的总体变化趋势.  相似文献   

5.
为了提高多变量复杂结构抗震可靠性设计计算速度和精度,本文采用了一种为较理想的计算方法-响应表面法,文中将结构的材料强度和荷载分布视为某特定分布的随机变量,通过少数几个点的随机抽样计算,得到了近似的极限状态函数g(x)我们将这个函数称为响应表面,然后通过二阶矩方法进行了可靠指标计算,得以了满意的结果,更重要的是在计算响应表面g(x)的过程中,本文采用了逐步求精的插值方法,将g(x)在失效设计点附近区  相似文献   

6.
层析反演是速度建模中最重要的方法之一,结合偏移成像在成像域进行走时层析速度反演是当前比较成熟有效且广泛应用的技术.本文从高斯束偏移成像条件出发,在波动方程的一阶Born近似和Rytov近似下,推导了成像域走时扰动与速度扰动的线性关系,建立了成像域走时层析方程及其显式表达的层析核函数.该核函数的本质是有限频层析核函数,利用该核函数替换常规射线层析核函数可以明显提高层析反演精度.该核函数的计算关键是背景波场格林函数的计算,本文利用高斯束传播算子计算格林函数进而得到走时层析核函数,实现方式灵活高效且计算精度较高.基于高斯束传播算子的偏移成像与层析成像相结合进行深度域建模迭代,体现了速度建模与偏移成像一体化的思想.数值计算及实际数据应用证明了基于高斯束传播算子的成像域走时层析方法的有效性.  相似文献   

7.
粘弹性层状地基动力柔度系数的半解析解   总被引:1,自引:0,他引:1  
本文以三次样条函数和傅氏级数作为求解空间的基函数,应用基于复阻尼理论的Lagrange方程求计算粘弹性层状半无限地基在表面单位简谐力作用下的动力柔度阵。这种半解析半离散方法,将二维问题降为一维数值问题处理,同时利用了样条函数良好的插值性能和傅氏级数的正交性,大幅度节省计算工作量和内存要求,便于利用微机实现,并具有较高的计算精度,从而为进一步建立复杂地基的动力刚度阵和研究结构与地基相互作用问题提供一个简便有效的数值计算工具。  相似文献   

8.
动态随机有限元法   总被引:12,自引:0,他引:12  
本提出分析了不确定性结构随机动力响应的一种新的数据方法--动态随机有限元法。经的特点是引入反映结构振动特征的动态形函数,中引出了相应的计算列式。算例表明,本方法的计算效率及数值精度均优于一般的随机有了限元法,具有应用前景。  相似文献   

9.
时间域瞬变电磁法中心方式全程视电阻率的数值计算   总被引:32,自引:9,他引:32       下载免费PDF全文
给出了一种时间域瞬变电磁法视电阻率的数值计算方法,利用该方法可以容易地求出中心方式的全程视电阻率. 根据中心方式磁场垂直分量时间变化率Bz/T的核函数Y′(Z)的表现特征,以参数Z把整个瞬变过程分为早期阶段(Z>1.6)、早期到晚期的转折点(Z=1.6)和晚期阶段(Z<1.6). 首先分别得到早期视电阻率和晚期视电阻率的精确值,然后通过转折点构成一条完整的全程视电阻率曲线. 虽然磁场垂直分量Bz的核函数Y(Z)是参数Z的单值函数,但同样存在一个从早期到晚期的转折点Z=1.6,转换点两边仍然可以得到一条早期曲线和一条晚期曲线. 在数值计算中,当迭代步长ΔZ<0.005Z时,视电阻率的相对误差小于0.5髎. 理论模型和实际数据计算表明,与早期和晚期近似值比较,全程视电阻率具有更高的精度和分辨率.  相似文献   

10.
提出了一种新的数值方法——基本结构有限元法. 从经典的伽辽金法(Galerkin method)出发,首先定义了基本结构插值基函数,在此基础上建立了基本结构方程,然后用有限元法进行进一步的具体实现. 该算法被成功应用到大地电磁测深一维连续介质正演计算中. 采用了6种不同的单元插值基函数进行计算,并对计算结果进行了比较和分析.  相似文献   

11.
如何有效协调地震波数值模拟过程中的计算精度和计算效率,一直是许多地球物理工作者密切关注且努力寻求解决的问题.这一问题求解的好坏,将直接影响到人们正确认识地下构造的地球物理响应.本文从离散Shannon奇异核理论出发,推导了基于Shannon奇异核的交错网格褶积微分算子,分析并求出了影响算子精度窗函数的最优化参数.通过将...  相似文献   

12.
层析反演的核心是敏感核函数的求取和旅行时反演方程的建立.在射线理论下,射线追踪技术是计算敏感核函数和旅行时的主要方法,因此初至层析的实现需要有一种可以准确追踪初至波路径的射线追踪方法.本文推导了抛物旅行时插值射线追踪公式,提出了基于抛物旅行时插值的最短路径射线追踪方法.理论模型和实际资料处理结果表明,利用该方法实现初至层析速度建模可行.  相似文献   

13.
基于WNAD方法的非一致网格算法及其弹性波场模拟   总被引:3,自引:2,他引:1       下载免费PDF全文
加权近似解析离散化(WNAD) 方法是近年发展的一种在粗网格步长条件下能有效压制数值频散的数值模拟技术. 在地震勘探的实际应用中, 不是所有情况都适合使用空间大网格步长. 为适应波场模拟的实际需要, 本文给出了求解波动方程的非一致网格上的WNAD算法. 这种方法在低速区、介质复杂区域使用细网格, 在其他区域采用粗网格计算. 在网格过渡区域, 根据近似解析离散化方法的特点, 采用了新的插值公式, 使用较少的网格点得到较高的插值精度. 数值算例表明, 非一致网格上的WNAD方法能够有效压制数值频散, 显著减少计算内存需求量和计算时间, 进一步提高了地震波场的数值模拟效率.  相似文献   

14.
大地电磁自适应正则化反演算法   总被引:36,自引:11,他引:36       下载免费PDF全文
针对大地电磁正则化反演中正则化因子的选取困难问题提出了自适应正则化反演算法(Adaptive Regularized Inversion Algorithm, ARIA). 在该算法中, ①提出了一种新的数据方差处理方法:数据方差规范化,使得数据方差的大小只对数据的拟合发生影响,不对数据目标函数和模型约束目标函数的权重产生影响,从而减少了正则化因子取值的影响因素;②提出了粗糙度核矩阵的概念,并给出了由基本结构插值基函数计算粗糙度核矩阵的公式,使得模型目标函数的构建更为简便、直接;③根据数据目标函数、模型约束目标函数和正则化因子之间的关系,提出了两种正则化因子自适应调节方法. 本文详细阐述了最平缓模型约束下的大地电磁一维连续介质反演的ARIA实现,以几个算例的分析比较来说明ARIA的有效性.  相似文献   

15.
采用代数重建技术重建三维物体,提出了一种基于splatting原理的ART权值计算方法,以增强重建图像的空间分辨率.该算法基于三维锥束扫描方式,采用双三次样条局部基函数作为插值核,并利用其足印函数对物体进行插值来计算像素的权值.与常量插值相比,这种方法对物体的灰度值有更精确的近似,在相同的实验条件下可重建出分辨率高的物体.给出了适用于三维锥束重建的splatting算法的具体实现,并对实验数据进行了实物重建及误差分析.  相似文献   

16.
大地电磁测深二维正演中辅助场的新算法   总被引:6,自引:2,他引:4  
马为  陈小斌  赵国泽 《地震地质》2008,30(2):525-533
文中提出一种计算大地电磁二维正演辅助场的新方法。在主场满足线性插值基函数的前提下,通过构建二次插值基函数来提高辅助场计算中数值求导的精度。通过与线性插值函数法的计算结果进行对比,证明该算法在基本不增加计算量的同时可以大幅度提高大地电磁二维正演精度。文中还研究了网格对辅助场计算的影响,结果表明二次插值法较线性插值法对网格的要求更为宽松。此外还以均匀半空间模型为例讨论了地表网格纵向间距与模型之间的关系,这对正演计算中的网格设计有参考意义  相似文献   

17.
长距离GPS基线及卫星定轨计算精度   总被引:1,自引:0,他引:1  
利用6个GPS固定观测占3天的观测资料,采用伯尔尼3.5版软件用多种方法计算了长距离GPS基线边及GPS卫星轨道参数,结果表明,对800~3000km的GPS基线,利用精密星历计算,边长及其经纬度分量重复测量精度要达10^-8~10^-9,同精密层历相比,2天轨道弧段的定轨精度为1~4m(坐标中误差)显然,采用IGS(国际地球动力学GPS服务处)计算中心的精密星历可以满足地球动力学及地震预报研究对  相似文献   

18.
蕴震系统的Kolmogorov熵及其可预报性   总被引:3,自引:0,他引:3  
本文简单地回顾了从热力学熵到拓扑熵的熵概念发展历史以及熵理论在地震预报中的应用,着重介绍了Kolmogorov熵的理论背景、计算方法、估计方法以及它与系统动力学行为的可预报性问题的关系。以京津唐地区为例,对时间域计算了K熵的唐山大震前后的变化情况,得出大震前系统存在K熵下降的趋势,震后又回升。以多台的宽间域综合方法计算得到K熵大于零,对于10^-5~10^-6的观测精度其可报期限为4-5年。  相似文献   

19.
本文研究了在高层建筑抗震分析中的阻尼计算问题,在考虑土壤-结构相互作用的情况下,研究了等效模态阻尼比的五种计算方法,振型加权法、应变能加权法、动能加权法、传函数法和比例阻尼法,文中输入EL-Centro地震记录,对一个十层楼房的实例进行了数值计算,并对各种方法所得到的阻尼和响应值进行分析和比较,最后对这些方法作出了评价。  相似文献   

20.
分形插值地震数据重建方法研究   总被引:10,自引:3,他引:7       下载免费PDF全文
对分形插值方法作了较详细的探讨,给出了分形插值函数的显式表达方式,同时给出了垂直比例因子的局部显式表达式,旨在提高地震道插值重建的精度及突出局部信息,并从单道地震图的角度分析其在地震道插值重建中的应用效果.利用该方法对理论模型和济阳坳陷实际地震台站资料进行了重建处理,结果表明,分形插值重建的地震道是原始地震道的良好近似,缺失道的振幅和相位都得到了很好的恢复.该法克服了随机分形插值方法必须进行多步迭代的弱点,提高了计算效率.通过对单道地震图插值重建结果的分析,说明了本文分形插值方法具有较高的精度和较高的效率,有深入研究的潜力.本文提出的显式分形插值方法既能够突出地震道数据的局部信息,又较好地保持了地震道数据的总体变化趋势.  相似文献   

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

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