一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

极区惯导编排中地球近似模型的适用性分析     

林秀秀  卞鸿巍  马恒  王荣颖 《测绘学报》2019,48(3):303-312

极区基于地球第二近似椭球模型横向坐标系编排复杂,编排本身存在原理性误差,而地球第一近似球体模型虽然存在几何误差,但在极区惯导编排简单,且相对于中低纬度地区,第一近似几何误差较小。惯导极区编排地球模型的选用问题随着极区惯导的深入研究逐渐凸显。针对这一问题,在定义适用于极区的地球近似模型的基础上,推导建立了基于地球第一近似模型几何误差引起的惯导系统误差方程,通过仿真分析对比地球近似模型的惯导误差,分析地球近似模型的适用性。仿真结果表明,地球第一近似模型几何误差会随惯导运行时间发散,而第二近似模型的编排误差随时间变化不大,当惯导运行时间长则地球第一近似模型不再适用。  相似文献   

基于近似真地表浮动面叠前深度偏移成像技术应用研究     

窦强峰  罗勇  杨晓海  谭佳 《物探与化探》2022,(2):444-450

常规基于水平地表假设的CMP浮动面处理方法不适用于准噶尔前陆冲断带复杂构造区成像。选择近似真地表浮动面并统一作为时间域和深度域处理面,从叠前时间域预处理包括优选近地表底界面,地表一致性近似真地表静校正,地表一致性双平方根动校正求取速度及剩余静校正,使时间域预处理到叠前深度偏移保持流程、参数的一致性。该基准面的选择使波场走时更接近实际传播路径,成像精度更高,地质构造准确。  相似文献   

高温高压下黄铁矿热力学性质的第一性原理研究   总被引：1，自引：0，他引：1  

刘善琪  李永兵  石耀霖 《岩石学报》2018,34(6):1813-1818

黄铁矿是自然界中分布最为广泛的硫化物矿物,同时也是重要的造矿矿物,在金属矿床、沉积岩、变质岩、花岗岩、基性-超基性岩浆岩、以及地幔岩中都有大量出现。因此,研究黄铁矿在不同温度压力下的热力学性质可以为深入探讨与黄铁矿有关的成岩、成矿、成藏问题提供有用的矿物学依据。本文利用基于密度泛函微扰理论的第一性原理方法,采用准谐近似计算了黄铁矿在高温高压下的热力学性质。我们计算的黄铁矿的晶格常数、零压下的体积模量及其对压力的导数与前人的实验及理论计算结果吻合得很好,零压下等压热容和熵随温度的变化与实验结果有很好的一致性。尤其是,本文计算了直至2500K、100GPa的高温高压下黄铁矿的等温体积模量、热膨胀系数、热容和熵等热力学性质。这为在有硫参与的情况下,人们开展下地壳-岩石圈地幔深度的地球动力学模拟和建立地球物理模型提供了有用的信息。  相似文献   

快速近似主成分分析算法   总被引：2，自引：0，他引：2  

朱述龙 《遥感学报》1999,3(1):43-47

通过分析现有主成分分析算法的不足,研究了如何利用小波包算法实现快速近似主成分分析算法的问题,并对两种算法的复杂度进行了比较。实验结果表明：提出的快速近似主成分分析算法在精度和速度两个方面都具有明显的优势。  相似文献   

单三角锁近似平差计算的虚拟附合导线法     

姜晨光  闫桂峰  王世周 《测绘与空间地理信息》1999,22(1):6-7

本文指出了单三角锁近似平差计算传统算法的局限性;提出利用虚拟附合导线的形式,解决单三角锁计算问题的新算法;给出了新算法的详细计算过程。  相似文献   

冰川均衡调整重力与径向位移近似关系的不确定性     

贾路路  汪汉胜  相龙伟 《地球科学》2014,39(7):905-914

根据不同地幔粘滞度的冰川均衡调整(glacial isostatic adjustment, GIA)模型, 研究了地球内部各个圈层对GIA粘性重力扰动速率的贡献, 检验了粘性重力扰动速率与径向位移速率的近似关系及其是否独立于地幔粘滞度, 同时利用绝对重力和GPS(global positioning system)径向位移数据从实测角度对Wahr的近似关系进行比较和验证.结果表明: 岩石圈对GIA重力扰动速率和大地水准面异常速率的贡献都超过了86%, 而岩石圈以下5个圈层的总贡献不大于14%;利用近似关系, 由重力信号转换的径向位移速率与有限元模拟的结果相对差异大约为15%, 且相对差异的大小不依赖于地幔粘滞度的变化; 根据北美绝对重力和GPS径向位移数据得到实测的粘性重力-径向位移比值为0.141±0.014 μGal/mm, 与Wahr的理论值(0.154 μGal/mm)非常接近, 相对差异仅为9.2%.因此, 定量给出了粘性重力-径向位移近似关系的不确定性为9.2%~15.0%, 为利用此近似关系分离GIA和现今地表质量变化粘弹信号的不确定性估计提供了重要参考.   相似文献   

万米AUV下潜深度估算研究          下载免费PDF全文

汪明星  常津铖  任翀  陈朝晖 《海洋技术学报》2022,(3):56-65

AUV是深渊观测的重要平台,万米级AUV是国际上深渊观测的研究热点。传统对于深度的测算主要依据EOS-80方法,使用压强数据进行深度估算,但由于深海中水层密度、重力加速度变化等因素影响,实际下潜深度有待进一步研究。本文以2021年青岛海洋科学与技术试点国家实验室研制的万米级AUV样机在马里亚纳海沟开展的下潜试验为例,通...  相似文献   

高斯牛顿迭代法解算非线性Bursa-Wolf模型的精度分析     

邓兴升  孙虹虹  汤仲安 《测绘科学》2014,(5):93-95

在参数估计中,非线性模型直接精密解算缺乏高效的方法,线性近似存在模型误差,而线性化取高次项导致模型复杂不具实用性。研究表明经典边角网可以不考虑线性近似的模型误差问题,本文以任意旋转角的非线性Bursa-Wolf模型参数解算为例,以不存在模型误差的直接严密解为参照对比,采用线性近似模型的高斯牛顿迭代方法解算非线性模型。试验结果显示,线性化取一次项虽然存在模型误差,但高斯牛顿迭代能以指定精度收敛,可获得更优于非线性严密直接解的精度,该发现对非线性模型解算的研究具有参考价值。  相似文献   

非线性阻尼吸振器应用于主系统受迫振动的一阶近似解   总被引：1，自引：0，他引：1  

陈载思  王树勇等 《广东海洋大学学报》2000,20(4):65-68

研究利用非线性阻尼吸振器应用于主系统受迫振动的情形，运用伽辽金法求出系统运动和微分方程的近似解。  相似文献