二维不规则界面介质局域离散波数边界积分法的研究

来源 :北京大学 | 被引量 : 0次 | 上传用户:wei71
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
不规则界面介质中散射波场的理论研究由Aki-Larner(1970)开始,自此就这一问题发展出一系列的模拟方法,主要有:纯数值法(有限元、有限差分、谱元等)、高频射线法、边界法等。边界法是基于边界积分方程的一种半解析的方法,对于界面的散射问题,它优于其它方法。从Wong(1975)首次将边界积分方程应用于模拟不规则界面介质中散射波场起,边界法得到了迅速发展,发展出各具特点多种计算方法,按所求的未知量性质分类主要有:以位移和应力为未知量的直接积分法,以界面分布的力源为未知量的间接积分法,以基函数的系数为未知量的基函数法等。按所应用的Green函数的性质可以分为空间域法和离散波数域法。边界法一个共同的缺点是需要求解大型满矢线性方程,因此计算量往往巨大,这使得一些问题如高频或大规模问题,模拟计算较难进行,因此设法减少边界法的计算量是一非常值得研究的问题,这正是本文的主要目的和任务。 现有边界法的改进是在离散波数域中进行的。有如下几种:一Kawase法,它通过应用半空间Green函数,将地表分布的未知量限定在地表起伏部分,从而减少计算量,但是半空间Green函数的普适性较差;二Campillo法,它的思想与Kawase相似,它通过使用层状Green函数,将模拟过程中的未知量由分布于全计算区域内限定为分布在界面的不规则部分上,但是因为层状Green的获得需要数值计算,所以计算量的减少并不是特别显著;三是Bouchon法,它通过设定阈值,将未知量系数矩阵中小于阈值的量定为零,将原来满矢线形方程变为有限带宽,从而减少计算量,但这种处理会影响计算结果的精确度。与上面三种方法不同,本文提出另外一种提高计算效率的新的改进方法“局域离散波数边界积分法”。它的改进基于Bouchon—Campolli离散波数法,首先将计算区域中的未知量分为对应于界面水平部分的力源和起伏部分的力源,充分利用Fourier变换,以小的计算量求出两者之间的数学关系,通过简单的代数运算,将未知量限定为起伏界面上分布力源,从而减少计算量。理论及模型验证,以上方法保留了原Bouchon—Campillo离散波数法的计算精确度,而计算效率至少提高一个量级。 在文章的第四第五章,将“局域离散波数边界积分法”(简称局域离散波数法)应用于理论和实际地震波场的模拟中,研究了不同尺度地形的频率响应特征,得到地形尺度与频率响应的似共振的响应关系;接下来模拟分析了多层介质中SH波和P-SV波的激发传播特征。对于实际地震波场的模拟主要进行了两部分,一是分析1999台湾集集地震中的地形效应,二是讨论张北地震在北京市引起的强地面运动,分析了北京凹陷第四纪沉积盖层和第三纪沉积盖层的地点放大效应。震中张北县距离北京市200多公里,在这么长范围内,利用通常的数值计算法(如有限差分等)及边界法,将震源与凹陷盆地放入同一个计算区域内模拟强地面运动,都会因存储空间巨大和计算量巨大,较难实施。通常多利用分区域连接法计算,但分区边界会代入一定的误差,计算不够精确。而“局域离散波数法”可以达到较远震源与凹陷盆地的同区域模拟,避免分区域模拟引入的误差。
其他文献
目前,中学正在进行新课程改革,要求将学习的主动性交给学生,让学生成为教学的主体。要怎么做才能实现这一目标呢?文言文教学从来都是一个难题,五四运动以来的白话文运动把我们带入了一个新的时代,对于现在的孩子来说通过文言文和古人沟通、交流是一个不可思议的事情,如果能够在文言文的教学上取得突破,对于中学语文教学来说意义重大。  一、引导学生反思历史  从中国人信奉“一日为师,终身为父”,西方人则“吾爱吾师,
重力场是反映地球介质密度变化和各种环境(固体地球潮汐、内部热流、固体和液体之间的质量交换及迁移、表面负荷和地震构造运动等)下动力学特征的最基本和最直接的物理量。重
任何一种风尚的兴起,在于其皇家所规定的、可以维护统治、抚平天下的制度,北宋对画院评判所制定的标准、绘画风貌和行事法则,无不是体现了制度的准绳和皇家的意志,显示出皇家
求学路漫漫,不知不觉间,已过了十几载,而我对它的感情,是又爱又恨的……  从小便知道,学习是迈向成功的一条最好途径,没有学问,便难在这个社会上立足,更不谈实现个人的理想。我想通过学习实现我的理想,但同时又害怕自己做得不够好,落于人后。这种情感使我时刻不敢放松,甚至神经紧绷,刻苦学习。都说熬过漫长黑夜就是一片蓝天,我单纯地想只要迈过高中这道槛,后面的就是阳光大道。有时单纯是一种美,也是一种动力。  
波动方程正演模拟是研究地震波传播规律的主要手段,是探测地球内部结构的基础。目前,大多数关于地震波传播的研究工作,都是将地球近似为一个理想弹性体。但是,地震波在传播过程中
北京航科文化传媒有限公司是北京海鹰科技情报研究所的全资子公司,是致力探索发展“军民融合”道路而诞生的年轻生力军。多年来,航科公司始终秉承“国家利益高于一切”的企业
本文通过定期观测和遥测台网采集了大量监测数据,并建立了地质灾害监测地理信息系统。开发了枇杷坪地质灾害监测地理信息系统,建立了滑坡三维地层模型并将其可视化。论文取
首先介绍了本论文的一些基础准备工作,包括V5-2000仪器时间序列的读取、噪音处理以及基本参数的计算。V5-2000仪器附有自己的一套数据处理软件,但它不提供读取原始时间序列的软
地震体波走时层析成像是壳幔精细速度结构成像的重要途径。基于高频近似射线理论的地震波走时反演是建立在一定的模型参数化和正演射线追踪基础之上的。不同的成像方法主要区
磁暴作为磁层物理一个重要的动力学过程,其研究一直引起人们的广泛关注。磁暴期间,地磁活动加剧,全球对流增强,同时,大量电离层上行离子注入等离子体片中,这些离子在相应机制作用下