常Q滞弹性介质地震波动数值模拟——时域本构优化逼近          下载免费PDF全文

谢志南  郑永路  章旭斌 《地球物理学报》2018,61(10):4007-4020

在地震波动模拟中计入常Q滞弹性阻尼,可有效降低模拟波形的误差.就时域有限差分和有限元模拟而言,常基于广义标准线性体建立阻尼介质的时域本构逼近.广义标准线性体由若干标准线性体并联得到,增加标准线性体个数能有效提高模拟精度,但计算量及计算存储将成倍增长.目前尚未有普适的标准线性体个数优化取值方案.本文基于广义标准线性体参数的非线性优化拟合方法,详细分析了时域本构逼近误差的影响因素,清楚揭示了逼近误差仅取决于频带宽度,与频带上下限取值无关这一特性,阐明了构建具有普适性标准线性体个数优化取值方案的可行性.论证了波形模拟精度主要取决于波传播距离与模拟波长的比值以及标准线性体的个数取值.综合考虑上述两个控制因素,结合在波动正反演问题中广为采纳的波形时频误差衡量准则,对不同Q值介质给出了标准线性体个数优化取值表.进一步,本文提出采用不同个数标准线性体以近似不同Q值的阻尼介质时域本构,解决了以往波动数值模拟中统一采用相同个数标准线性体而导致的计算量及计算存储浪费或模拟精度低下等问题,并基于数值实验验证了这一方法的精度.本文工作对推进滞弹性介质波动数值模拟及其在全波形反演问题中的应用具有理论价值和实践意义.  相似文献   

巴颜喀拉块体东缘形变及九寨沟地震孕震环境数值分析   总被引：6，自引：3，他引：3       下载免费PDF全文

庞亚瑾  程惠红  张怀  石耀霖 《地球物理学报》2017,60(10):4046-4055

2017年8月8日四川发生九寨沟M7.0地震,是继2008年汶川M8.0地震后发生在巴颜喀拉块体东部的又一强震.现今GPS速度观测数据显示,2008年汶川地震前后的1999-2007年和2011-2016年两个时间段内巴颜喀拉块体东部地表速度场存在明显的差异.本文以实际GPS速度观测资料为约束,构建三维有限元地球动力学模型,分别计算分析了两个时段内震源区及周边现今地壳形变、弹性应变能和应力积累特征,进一步探讨汶川地震的发生对九寨沟地区变形及应力的影响.数值模拟结果显示,汶川地震之后(2011-2016年)巴颜喀拉块体东部的地壳形变、弹性应变能积累及应力积累速率均明显大于震前,增加量值达1.5-3倍;九寨沟地震发震断裂上库仑应力增长率在1999-2007年约为0.7 kPa·a~(-1),2011-2016年间增至1.2 kPa·a~(-1).上述结果表明,现今巴颜喀拉块体东部地壳应力积累过程有利于左旋走滑型九寨沟地震的发生,汶川地震的发生调整了区域应力状态,加速了九寨沟地震的孕育过程.  相似文献   

基于非结构化网格的最小二乘逆时偏移   总被引：1，自引：1，他引：0       下载免费PDF全文

杨凯  张剑锋 《地球物理学报》2017,60(3):1053-1061

相对于传统的逆时偏移,最小二乘逆时偏移具有更高的成像质量,这种改善是通过迭代反演来获得的,另外其精度与效率依赖于求解波动方程算法的精度与效率.本文给出了基于非结构化网格的最小二乘逆时偏移的方法,该方法能够充分结合最小二乘算法与非结构化网格精细刻画地下界面以及随速度自适应剖分的优点;并采用带补偿的拉普拉斯滤波算法,来消除梯度计算中的低频噪声,从而加速目标函数的收敛速度.通过简单倾斜模型以及复杂Marmousi模型测试,显示了该方法的有效性和潜力.  相似文献   

基于平面波加速的VTI介质最小二乘逆时偏移   总被引：1，自引：1，他引：0       下载免费PDF全文

李振春  黄金强  黄建平  朱峰 《地球物理学报》2017,60(1):240-257

地震各向异性集中表现为速度各向异性,势必影响地震波运动学特征.传统声波逆时偏移(RTM)和最小二乘逆时偏移(LSRTM)没有考虑介质各向异性特征,导致反射波不能正确归位、同相轴出现扭曲及寻优速度慢或不收敛等,VTI介质逆时偏移(VTI-RTM)矫正了声波成像的不足,但仍存在低频干扰严重、中深部成像不佳、振幅保持差等缺陷.为此,本文首先实现了VTI介质最小二乘逆时偏移(VTI-LSRTM)方法,为了节省I/O及内存需求并提高效率,进一步引入平面波编码技术,提出了一种基于平面波加速的VTI介质最小二乘逆时偏移(VTI-PLSRTM)策略.在此基础上开展了简单模型及复杂Marmousi模型成像试验,并与标准逆时偏移剖面对比表明:本方法能够校正各向异性造成的相位畸变,且在迭代中自动压制串扰及低频噪声、补偿中深部能量,是一种兼具质量与效率的保幅成像策略;对速度误差的敏感性测试说明该方法需要相对正确的偏移速度及Thomsen参数模型.  相似文献   

基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移   总被引：2，自引：2，他引：0       下载免费PDF全文

黄金强  李振春 《地球物理学报》2017,60(2):704-721

近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制(η0)和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明:①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略.  相似文献   

鲜水河断裂带北西段跨断层基线变化机理探讨   总被引：1，自引：1，他引：0       下载免费PDF全文

李玉江  陈连旺  刘少峰  季灵运  占伟 《地球物理学报》2017,60(2):554-563

2013年芦山M_s7.0地震发生后,鲜水河断裂北西段的侏倭、虚墟等跨断层测点基线测项均出现不同程度的异常变化,异常核实结果认为基线异常属实.该异常变化是否与芦山地震的发生存在一定关联性,值得深入研究.本文首先分析芦山地震发生前后鲜水河断裂北西段不同基线长度的变化特征.基于鲜水河断裂带及邻区三维非线性黏弹性有限元模型,以GPS观测资料、同震静态滑移量作为约束,通过开展多组数值模拟实验,探讨下地壳不同流变特征下,芦山地震发生前后不同阶段,不同跨断层基线长度的时序动态变化特征.初步研究结果表明,1)侏倭、虚墟基线测点原始观测资料均表现出地震发生时基线长度减小,发生后短时间内反向快速恢复及增加的协同变化;2)在以GPS观测数据为约束的块体间差异性运动的动力学边界条件下,鲜水河断裂表现为左旋走滑运动;3)模拟结果显示,地震发生时测点间距离迅速减小,表现为断层右旋的特征,这与实际观测资料反映的结果一致;4)鲜水河断裂带两侧地块的下地壳黏滞系数分别取10~(18)Pa·s、10~(19)Pa·s时,由于黏弹性松弛效应的影响,造成基线长度在地震发生后短时间尺度内快速增加.对比分析认为,地震后短时间内实测资料反映的基线长度快速增加的特征可能是黏弹性松弛效应与构造作用共同作用的结果.  相似文献   

水下多路液压快速接头的夹头结构分析与优化     

陈再玉  王豪  肖易萍  吴志星  肖德明  任广新  陈家旺 《海洋工程》2021,39(6):143-151

深海油气田水下生产系统装备之间的精准连接是国际公认的技术难题,对接一般通过深海连接器完成。弹性夹头是水下连接器的重要组成部分,其性能关系到海洋油气生产系统的安全性和可靠性。弹性夹头的主要作用是在水下机器人的辅助下连接和锁紧连接器的固定端和移动端。采用ABAQUS软件的显式动力学分析模块计算夹头夹紧心轴这一过程中的应力、反作用力状态。三维模型仿真结果表明,锁紧过程中弹性夹头的最大应力不超过430 MPa,小于所选取的材料屈服强度480 MPa。采用正交试验和方差分析的方法研究夹头的长度、夹头与心轴的间隙等几个关键尺寸从而分析影响夹头应力和连接所需推力的最显著因素。对夹头的结构进行优化设计,同时也为夹头心轴结构的设计提供了理论依据。最后,通过试验研究对样机进行测试,测试结果表明,夹紧动作所需的最大推力在700 N以下。  相似文献   

基于滑动窗口导纳技术反演西太平洋中部岩石圈的有效弹性厚度     

凌子龙  赵俐红  彭祎辉  吴招才  曲彦丞 《海洋地质与第四纪地质》2021,41(1):138-146

西太平洋中部地区是西太平洋板块边缘沟-弧-盆体系构造演化的关键区域,其地质特征与构造演化一直是地学家关注的焦点问题之一。开展岩石圈有效弹性厚度的研究对于认识该区域的形成演化具有重要的科学意义。本文采用滑动窗口导纳技术,并在挠曲模型中考虑了表面荷载和内部荷载同时存在的情况,计算得到该区域的岩石圈有效弹性厚度(Te)。计算结果显示,研究区的Te值整体上为0~50 km,其变化基本上与构造单元相吻合,且与主要的构造边界密切相关。除海底火山具有相对较小的Te值(15~20 km)外,太平洋板块整体上具有较强的岩石圈强度(25~30 km)。马里亚纳海沟和菲律宾海沟的岩石圈强度从外隆起到海沟方向表现为明显的减弱,表明岩石圈由外隆起向海沟发生了弱化。帕里西维拉海盆西部相较于东部具有较弱的岩石圈强度,这可能与海盆的非对称扩张有关。卡罗琳板块的岩石圈整体上表现为相对均一的低Te值特征(<15 km)。欧里皮克海隆、卡罗琳海岭和索罗尔海槽的Te值为3 km,这可能是强烈的火山作用所导致的结果。  相似文献   

青藏高原挠曲均衡重力异常特征及其地质意义          下载免费PDF全文

张星宇  杜劲松  陈超  梁青 《地质科技通报》2023,42(2):223-233

自新生代印度板块的块持续碰撞与俯冲作用下,青藏高原经历了快速隆升与复杂的岩石圈改造过程,但高原现今的垂向动力学机制和地壳形变特征仍然存在争议。基于非均一有效弹性厚度的挠曲模型,利用地形和地球重力场模型数据,计算了青藏高原及邻区的挠曲均衡重力异常。结果显示,青藏高原的均衡重力异常在-120~90 mGal之间,高原中部为明显的正异常特征,边缘为显著的均衡负异常。极小值出现在青藏高原西北部及其相邻的帕米尔高原,极大值则出现在与之紧邻的喜马拉雅块体西北部。此外,在青藏高原北面和东面,塔里木盆地和四川盆地显示出大片的均衡正异常。这些特征说明青藏高原及邻区地壳现今处于非均衡的状态,在板块碰撞挤压作用下,老的块体地壳整体发生抬升,导致了均衡正异常特征;而年轻的造山区域,地壳形变主要表现为地表抬升与下地壳强烈增厚,形成了均衡负异常。在高原中部和北部,均衡调整方向与地壳垂向运动趋势相一致;但在高原南面（喜马拉雅块体）和东面（四川盆地）,均衡调整方向与地表形变观测结果相反。这说明印度板块碰撞与俯冲仍然控制着青藏高原南部、东部及其相邻块体的地壳形变过程,然而在更北的区域,地壳正通过均衡调整恢复均衡状态。   相似文献   

弹性波CT技术在大足石刻岩体破碎带探测中的应用          下载免费PDF全文

孔得天  李高  郑旭辉  石金山 《CT理论与应用研究》2018,27(1):35-44

大足石刻是中国现存的重要石质文物之一,具有重要的历史和文化价值,其中有的石刻造像曾被列为国家石质文物保护“一号工程”;最近,文物保护部门又启动了对危岩体的加固和防渗工程。为确定大足石刻宝顶山大佛湾岩体中破碎带发育情况,应用跨孔弹性波CT技术对指定位置的工程岩体进行了探测。CT探测得到的岩体弹性波速成像清楚地揭示了两个钻孔之间岩体中破碎带的规模和产状,与岩体所在的区域地质构造应力场作用协调一致,证明了探测成像结果的可靠性,为大佛湾岩体的渗流和稳定性评估提供了依据。   相似文献