共查询到19条相似文献，搜索用时 156 毫秒
1.
二维层状介质点电源边界单元解 总被引：1，自引：0，他引：1
本文用付氏变换将点源二维地电断面的三维边值问题变成二维边值问题,然后讨论了边界元法在解决层状介质中存在二维柱体时点源场视电阻率计算问题中的应用。由于选取的基本解避免了对水平层状界面的积分计算,从而大大减少了边界剖分节点数。在计算中设法略去了对边界上位场法向导数的求解,使未知量大为减少,从而便利了求解二维边值问题,并由付氏反变换求得三维空间中的位。本方法所需计算机内存少,原始数据的输入简单,且有较高的计算精度。 相似文献
2.
3.
论文使用高分辨有限体积法模拟二维声波在多层介质中的传播现象.该方法需要求解双曲型偏微分方程的黎曼解,数值解在强速度变化的分界面上满足数值通量守恒,可以有效地计算数值解的不连续性.根据数值实验给出的波场快照和接收点上的波形变化,可以清晰地观察到声波在介质分界面上的多次反射波,说明该方法不但可以得到高分辨率的波场快照,而且... 相似文献
4.
5.
根据瞬变电磁场扩散方程与波动方程间存在的数学对应形式,将已知瞬变电磁场数据转换为虚拟波场数据,在此基础上,提出了用克希霍夫积分法将所转换的波场从地面向地下反向延拓的方法.在延拓的数值计算中,为提高计算效率,采用三维边界元技术,把边界积分分解为诸三角单元积分的积分之和, 建立了曲面延拓方程式;通过地面测点的波场值,求出地下某一点的波场值,实现了地表为曲面的向下延拓成像计算.通过对理论模型计算和实际资料处理,证明了该方法可以增强瞬变电磁法识别地下电性分界面的能力,使瞬变电磁法对地下目标体的三维精细探测成为可能. 相似文献
6.
一般工程的地震反应分析中地震波假设为竖直向上入射的体波,场地自由波场可简化为成层半空间模型的一维波动问题,通过数值方法获得出平面自由场。在数值计算方法中以显式有限差分法为代表的时间域方法是当前地震波模拟的主流方法,其精度通常受到计算方法、模型尺寸、网格大小、边界设定等诸多条件的影响,数值算法稳定性还受到许多计算参数的限制。为了改进时域数值方法中存在的这些缺点和不足,本文结合分层均匀介质中的理论地震图合成方法提出了一种新的时域数值算法。其基本思想是将计算节点直接设置在波阻抗界面处,各个层内不再设置计算节点,然后根据地震波在波阻抗界面处的反射、透射关系及其在相邻界面之间的传播时间建立计算界面节点振动的数学表达式,据此计算得到各层界面节点的振动,并通过对自由面节点进行适当处理得到成层半空间出平面自由场。界面节点振动可看作为入射波通过波阻抗界面处产生的次一级子波源,本文称之为界面子波,因此把这一计算方法称为界面子波算法。数值算例表明该方法具有很高的精度和很快的计算效率。理论上还可将这一方法推广至二维波动情形,即地震波倾斜入射条件下成层半空间出平面自由场的计算。 相似文献
7.
8.
矩形网格有限差分法在地震波传播数值模拟方面具有计算速度快的显著优势,但该方法在处理复杂边界问题上存在着效率低的严重缺陷.本文针对分段光滑曲线边界定义了尖点处的一种正则导数,给出了矩形网格情形分段光滑曲线网格边界点法向导数的一种插值计算方法.采用矩形网格有限差分法对复杂边界地球介质模型进行地震波场数值模拟,并采用波场系列快照技术揭示地震波在起伏地表和复杂介质中的传播规律.模拟结果表明：法向导数插值计算方法为矩形网格有限差分法处理复杂边界提供了有效途径,采用波场系列快照技术可以清晰地展现地震波在反射界面的反射和透射规律、在尖点的绕射规律以及在自由表面的直达波和多次反射规律. 相似文献
9.
本文将小波有限元法引入到流体饱和多孔隙介质二维波动方程的正演模拟中,以二维Daubechies小波的尺度函数代替多项式函数作为插值函数,构造二维张量积小波单元.引入一类特征函数解决了Daubechies小波没有显式解析表达式所带来的基函数积分值计算问题,并推导出计算分数节点上Daubechies小波函数值的递推公式,从而构造出由小波系数空间到波场位移空间的快速小波变换.数值模拟结果表明该方法是有效的. 相似文献
10.
不规则界面介质中散射波场的理论研究由Aki-Larner(1970)开始,自此就这一问题发展出一系列的模拟方法,主要有:纯数值法(有限元、有限差分、谱元等)、高频射线法、边界法等.边界法是基于边界积分方程的一种半解析的方法,对于界面的散射问题,它优于其他方法. 相似文献
11.
12.
13.
14.
Tong JiangShozo NishiokaHideaki NagasakaEiich Kuribayashi 《Soil Dynamics and Earthquake Engineering》1993,12(3)
The reflectivity method for complicated irregular formations has been extended to three-dimensional irregular (layered) topography. The Aki-Larner technique is generalized to solve the integral equations for 3-D boundary conditions, and propagator matrices are enlarged to express the total wave-field. The formulation in terms of the propagator matrix has been applied to horizontal layers included in the irregularity. The model can simulate more complex, layered, topography. Frequency responses on the surface of 2-D and 3-D layered soft basins are calculated for incident plane waves. They compare favorably with the results of another boundary method. Transfer function versus dimensionless frequency has been also calculated for the 2-D layered soft basin. Significant differences between the layered model and the homogeneous model have been observed; therefore it is necessary to consider the complicated irregular formations at the actual site. 相似文献
15.
16.
用平面波延拓方程进行地震数据的叠前速度反演 总被引：1，自引：1，他引：1
本文讨论地震勘探数据的叠前速度反演方法及其在海洋地震勘探数据上的反演试验.反演主要的计算步骤是:1.采用Fourier-Hankel变换把球面波分解为平面谐波;2.用平面谐波的延拓方程将上行波与下行波同时向下延拓,并计算每一层底部的反射系数和下一层的波阻抗;3.用最小二乘法从波阻抗中确定该层的声波速度.重复第2步与第3步,直到某一预定深度时结束.通过反演试验,对地震振幅比例的改变,子波变形,以及第1层速度和密度的误差对反演方法的稳定性及其精度的影响进行了分析.还通过实际海洋地震勘探数据的反演试验,对这一方法在地震勘探中的应用前景作了论述. 相似文献
17.
提出了针对轴对称二维井间模型的一维、二维联合反演方法.该方法基于纵向成层背景地层的Green函数,以一维反演结果作为联合反演的迭代初始值,同时反演一维背景地层和二维异常剖面的电导率.采用递推技术计算Green函数的系数,可以很方便地同时得到Green函数对空间坐标的积分和对背景层电导率的微分,大大加快了计算雅可比矩阵〖WTHX〗M〖WTBZ〗的速度并使计算方便、准确.反演实例说明了该方法的有效性.反演实例还显示,背景层电导率较二维成像剖面的电导率收敛速度要快,测量数据的精度对二维剖面电导率成像质量的影响大于对背景层电导率分布的影响. 相似文献
18.
本文发展了积分方程法用于层状介质中三维不均匀体的大地电磁模拟算法（简称MT）.分为二个步骤:第一步,异常体用等效的散射电流代替,通过层状介质中的格林函数,建立以散射电流为未知参数的积分方程;第二步,把求得的散射电流乘上相应的格林函数,即得地面上的二次电磁场,由此而计算出各种MT响应. 文中采用了数值滤波与插值、群变换以及格林矩阵带状化三个方面的数值处理方法,提高了计算效率.通过与已发表的三维MT计算结果对比及格林函数互易性检验,表明了该算法的正确性.在此基础上,进行了数值模拟,初步讨论了三维MT曲线的畸变特点. 相似文献