等效交错网格高阶有限差分法标量波波场模拟     

段沛然  李青阳  赵志强  李振春 《地球物理学进展》2019,34(3):1032-1040

密度是岩石物性参数的重要组成部分,对储层评价、岩性解释和油藏描述等具有重要作用,因此实际资料应用中含有密度信息的正演建模是必不可少的.作为波动理论应用最广泛的实际资料处理技术——逆时偏移和全波形反演,其正演建模一般采用经典的二阶标量波方程,且不考虑密度的空间变化.而一阶速度-应力控制方程通过交错网格有限差分法正演建模,其精度高于二阶方程但计算成本过高,很难在三维实际资料中应用.鉴于非均质正演在实际生产中所面临问题,本文借用交错网格的思想,充分考虑空间变密度对波传播的影响,提出基于等效交错网格的高阶有限差分建模算法并应用于非均质标量波正演.文中从数学上给出一阶方程和二阶方程在数值模拟中的等价性证明,并分析了震源、边界条件和稳定性.最后通过简单层状模型进行数值测试,对比不同建模方法以验证新方法的准确性;利用高阶精度方法对Sigsbee-2a模型做正演测试,验证本文提出的方法的稳定性.  相似文献   

莺歌海盆地黄流组一段限定-非限定朵体构型及其演化     

姜平  李明军  周伟  李明 《地球物理学进展》2019,34(6):2451-2458

浅海陆架扇是当前海洋沉积研究的前沿.基于三维地震资料、钻测井资料,利用地震相分析和沉积动力学方法,研究限定朵体、非限定朵体的构型、沉积演化及其动力学机制.研究结果表明:(1)陆架盆地浊流沉积体系发育水道—堤岸复合体、限定朵体和非限定朵体3类沉积单元.(2)浊流沉积体系受陆架盆地地形的影响,早期发育限定朵体,晚期发育水道—堤岸复合体、非限定朵体.(3)浊流在陆架盆地坡脚发生水跃,浊流动力学特征(浊流流速、Fr数)的变化——高速、超临界流在水跃后转换为低速、亚临界流,并伴随着浊流活动的侵蚀和沉积活动,导致了限定朵体向非限定朵体的转化.朵体构型及其演化研究对该地区优质储层预测具有重要指导意义.  相似文献   

PCA算法在CRP地震道集优化处理中的应用研究     

周家雄  马光克  隋波  徐伟  张万龙  曾勇坚  卢回忆 《地球物理学进展》2019,34(4):1577-1583

  相似文献   

黏弹性介质中瑞雷波有限差分数值模拟与波场分析     

王常波  田坤  刘立彬  李国磊  张学涛 《CT理论与应用研究》2019,28(2):153-166

瑞雷面波经常被用来反演地表浅层横波速度，受到越来越广泛的关注。对瑞雷波的研究一般都基于完全弹性介质，而实际地层更接近黏弹性介质，对黏弹性介质中的瑞雷面波进行模拟更具实际意义。本文采用广义标准线性体模型来描述黏弹性介质，并采用交错网格有限差分法对考虑水平自由表面的黏弹介质进行正演模拟，再与弹性介质中的结果进行对比分析。首先采用非线性最优化算法根据期望常数品质因子直接求取松弛时间来拟合常Q模型，并给出广义标准线性固体的具体算例，实施自由表面条件时采用声学-弹性边界近似法，通过剪切模量不变来考虑自由表面上、下横向应力保持连续的条件。对于非自由表面，采用非分裂的多轴卷积完全匹配层来吸收波场。然后对几种典型的数值模型进行正演模拟计算，数值解与解析解的对比验证了本文方法的准确性与有效性，正演结果的对比表明波场尤其是面波频散会受黏弹性影响，因此有必要在面波勘探中考虑黏弹性因素。   相似文献   

福建及其邻区背景噪声面波群速度层析成像     

张丽娜  罗艳  陈智勇  谢志招  周蓝捷  曾祥方 《地震》2018,38(3):134-143

基于福建及邻区108个宽频带地震台站2016年6月到7月两个月垂直和水平分量波形连续记录, 利用相位加权叠加算法提高信噪比, 计算得到108个台站对的高质量经验格林函数。 对所获取的经验格林函数, 采用时频分析的方法在1~20 s频段内量取了大量高质量的Rayleigh波和Love波群速度频散数据。 在此基础上, 采用基于射线追踪的二维层析成像方法反演得到了福建及其邻区1~20 s的Rayleigh波和Love波群速度分布。 分辨率测试结果表明群速度分布的分辨率能达50 km。 成像结果显示1~10 s的群速度分布与地壳中上部地质特征有很好的一致性, 区内的福州盆地和漳州盆地在浅层结构中表现出明显的低速异常。 长周期的群速度则揭示了漳州西北的高地热区内中下地壳低速体, 政和-大埔断裂两侧的速度差异, 表明其可能是一个深大断裂, 并呈现明显的东西差异。  相似文献   

波粒相互作用导致环电流质子沉降的卫星共轭观测          下载免费PDF全文

王杰  袁志刚  余雄东  薛祖祥 《地球物理学报》2020,63(6):2131-2140

波粒相互作用是环电流损失的重要机制之一,但波粒相互作用导致的环电流离子沉降而损失迄今为止缺乏直接的观测证据.基于磁层及电离层卫星的协同观测,本文报道了发生在2015年9月7日,由电磁离子回旋波(EMIC波)导致环电流质子沉降的共轭观测事件.在等离子体层的内边界,Van Allen Probe B卫星观测到,存在EMIC波的区域和不存在EMIC波的区域相比,离子通量的投掷角分布的各向异性变弱.我们将Van Allen Probe B卫星沿着磁力线投影到电离层高度,同时在该投影区域内DMSP 16卫星在亚极光区域观测到环电流质子沉降.而且,通过从理论上计算质子弹跳平均扩散系数,我们进一步证实观测的EMIC波确实能将环电流质子散射到损失锥中.本文的研究工作为EMIC波导致环电流质子沉降提供了直接的观测证据,揭示了环电流衰减的重要物理机制:EMIC波将环电流质子散射到损失锥中,从而沉降到低高度大气层中而损失.  相似文献   

一种模拟非均匀介质中弹性波传播新的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方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

利用远震研究四川盆地及其周缘瑞利面波相速度和方位各向异性          下载免费PDF全文

王小龙  雷兴林  马胜利  张烨  郭志  闻学泽 《地球物理学报》2020,63(2):445-459

基于四川盆地及周边的245个宽频带台站2010年9月—2014年9月期间的远震记录,提取双台路径瑞利面波相速度频散资料,反演得到四川盆地20～120s的高分辨率瑞利面波相速度及各向异性空间分布.在丰富区域地球物理基础数据的同时,结合已有研究成果对地壳上地幔变形耦合进行探讨,结果表明短周期(20～30s)的相速度分布与四川盆地的地质构造特征相吻合,作为川滇地块、松潘—甘孜地块和四川盆地之间的边界——龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用;松潘—甘孜地块,特别是川滇地块中下地壳普遍表现为明显的低速异常,表明中下地壳相对软弱;而四川盆地的中下地壳整体呈现相对高速,表明四川盆地具有相对坚硬的中下地壳.研究区域东南角接近北扬子地块与南扬子地块的缝合部位,呈现高速异常.四川盆地南部和东南邻区不同周期均具有较强的各向异性,且快波方向较为一致,反映这些地区不同深度变形耦合较好.四川盆地西部、北部及东北部邻区,不同周期的各向异性快波方向变化较大,不同深度变形耦合较差.这些特征与绕喜马拉雅东构造结的物质流动被扬子地块的高速地壳阻挡的宏观认识基本一致.  相似文献   

基于激波捕捉类Boussinesq模型对岛礁地形上波浪传播的数值模拟   总被引：2，自引：2，他引：0  

刘雨诗  刘维杰  赵西增  陶金波 《海洋工程》2019,37(5):48-56

为了探究激波捕捉类Boussinesq模型在模拟岛礁地形上规则波和不规则波传播的可行性,采用基于完全非线性Boussinesq方程并具有激波捕捉能力的数值模型Funwave-TVD对规则波和不规则波在岛礁地形上的传播进行了数值模拟,通过与试验数据对比,分析模型中空间步长的影响,验证模型在模拟波高、平均水位分布以及波谱空间演变的能力,结果表明:采用合适的空间步长,模型能较好地模拟规则波和不规则波在岛礁地形上的传播和演化过程。对于规则波,较小的空间步长可改善破碎点处波高峰值的预测,并能更好地预测波浪破碎后波高的空间分布,相比结合经验破碎的Boussinesq模型,Funwave-TVD能更好地模拟规则波在岛礁地形上的破碎,以及破碎以后行进涌波的再生成过程;对于不规则波,Funwave-TVD总体而言能较好地模拟涌浪有效波高、次重力波的生成及空间演化和平均水位,但会低估礁坪上次重力波波高,较粗的空间步长也会低估礁坪上涌浪有效波高。  相似文献   

汶川M_W7.9地震起始破裂断层几何结构          下载免费PDF全文

宫猛  徐锡伟  李康 《地球物理学报》2020,63(3):1224-1234

本文收集使用紫坪铺水库台网记录到的汶川地震主震P波波形资料,利用P波反投影叠加法获取了2008年5月12日汶川M_W7.9地震起始破裂的时空演化过程.通过分析本次大地震起始破裂阶段(0～1s)破裂点在三维空间内的分布特征,确定了本次大地震起始破裂位置及起始破裂断层几何结构模型.得到以下结果:汶川地震起始破裂点位于31.013±0.002°N、103.392±0.002°E,震源深度为8.2±0.4km,发震时刻为2008年5月12日14∶27∶58.80±0.4.汶川地震起始破裂的最佳断层面走向为NE48°,倾向NW35°,起始阶段破裂的深度范围为地下7.5～9km.  相似文献