2D各向同性和各向异性介质中高阶交错网格有限差分波场模拟

来源 :长安大学 | 被引量 : 3次 | 上传用户:yyxu123
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在地震勘探中,习惯的做法是通过弹性波场分解技术把P波场和S波场分离出来以达到成像的目的。可是这个目标有时很难达到,特别是当速度场较为复杂的时候。有一种有效的多分量分析和模拟方法是直接求解弹性波动方程而得到纯P波场和纯S波场,得到的方程称为分离的弹性波动方程。在本文中,我们首先采用高阶交错网格有限差分法对比了两种这样的波动方程:一阶(速度-应力)和二阶(位移-应力)分离的弹性波动方程,与一阶(速度-应力)和二阶(位移-应力)完全(混合)的弹性波动方程。然后比较了得到的波场快照、共炮点道集的地震剖面和单独的合成地震图。模拟结果表明得到的结果是等价的,不管是一阶或二阶分离的弹性波方程都能获得纯P波场和纯S波场。将纯P波场和纯S波场叠加得到的混合波场与用对应的一阶或二阶完全弹性波动方程计算出的混合波场结果相同。求解这种混合方程的计算代价稍微低于分离方程。分离的方程主要用于实现分离的P和S波场,且可以用来测试多分量波场分离程序的有效性。二阶分离的弹性波方程提供了纯P波场和纯S波场的位移信息,因此是一个很好的选择。其次,我们还使用高阶交错网格有限差分法模拟各向异性介质中的混合波场,可以得到包含准纵波和准横波的混合波场。为了从各向异性介质中分离qP波和qSV波,使用各向同性介质中的散度和旋度算子是没有效的。我们要通过求解Christoffel方程得到的伪求导算子来分离纯qP波和qSV波。虽然伪求导算子可以得到纯qP波和qSV波,但它需要大量内存的计算,且实际模型的模拟比较复杂。
其他文献
当前电力青工思想状况的调查与思考刘志云跨世纪一代青年的思想状况如何,在很大程度上决定着中华民族的前途和命运。国有大中型企业的青工,其思想状况又如何呢?为此,山西电建二公
带着一份渴望与期待,怀揣一份激情和喜悦。今天,我们欢聚一堂,共同庆祝安徽省黄梅戏剧院建院60华诞。六十年,一个甲子的回眸,荣辱兴衰、风雨兼程,寄托了几代黄梅戏工作者的梦
赣中铁矿田是我国南方重要的铁矿类型。本文依托参与的江西赣中铁矿田资源远景调查及后续矿山接替资源找铁矿专项行动为研究基础,通过对研究区原有地质物探等综合资料学习研
本文主要研究电离层的时间变化性对解算GPS系统硬件延迟的影响。基于GPS系统基本的观测量(伪距和载波相位)可以得到电离层的总电子含量,这个总电子含量中包含有GPS卫星和接收
王铎在五十岁之后,写了很多草书精品,诸如《赠张抱一草书诗卷》、《杜甫诗卷》、《秋兴八首》等。而《王屋山图诗卷》无疑是其这一时期的代表作。此卷为1647年(清顺治四年)王铎56岁时所作,纸本,墨迹,纵26cm,横512.5cm,计87行,凡490字,以一幅王屋图山水起首,实在是不可多见。  从起首的王屋图可以看出,王铎笔下的山水,与当时影响巨大的松江派枯瘦的山水有很大不同。当时他曾旗帜鲜明地表达对松
地震资料反褶积、波阻抗反演是提高分辨率和储层预测识别的关键技术。反褶积的基本原理,是利用采集到的地震数据和已知的子波经过反褶积得到反射系数;波阻抗反演则是由反射系
志留纪末,加里东运动将上扬子海盆抬升为陆,广遭剥蚀,直至中石炭世早期地壳下降,由鄂西向西海侵,为一套半封闭的潮上~潮间带沉积,晚石炭世上升为陆遭受剥蚀,残厚各地不一。QLX地区处
地球外辐射带充满了能量高达MeV的相对论电子,给航天器的安全运行带来了极大的威胁。同步轨道区域运行着几百颗通讯、导航及气象卫星,这个区域的相对论电子始终是人们研究的重
本文收集整理了1973年以来全球329个M≥7.0浅源地震(H≤60km)的余震序列资料(M≥4.0,1年时间),统计了全部序列的最大余震震级M1、主震-最大余震震级差△M、主震与最大余震的震中
应该说,“实践是检验真理的唯一标准”是马克思主义的基本原则。但在我们这个以马克思主义为指导思想的国家,这条原则竟然有过如此复杂、曲折的遭遇,这不由使人再三深思。 历