首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 171 毫秒
1.
埋地目标体矢量电磁散射的一种快速正演算法   总被引:20,自引:4,他引:16       下载免费PDF全文
利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。  相似文献   

2.
三维多值走时地震波场重建及格林函数计算   总被引:7,自引:0,他引:7       下载免费PDF全文
针对三维复杂介质地震波传播数值计算中出现的多值走时情况,阐述了地震波场重建及格林函数计算中的困难,提出一种在相空间拉格朗日流形上的波场重建及格林函数计算方法.本文算法应用于三维叠前深度偏移处理流程时,可高效地获得均匀网格点上地震波多值走时及振幅的数值计算结果.文中还阐明两类数值计算过程判据,以在多值走时区域及计算的全区域中,控制射线追踪过程中的射线密度及格林函数计算精度.算例验证了本文方法的有效性.  相似文献   

3.
近场地震动格林函数的解析法与数值法对比研究   总被引:3,自引:2,他引:3       下载免费PDF全文
分析了在均匀弹性介质中,分别用解析方法和三维波动有限元数值方法计算无限全空间近场地震动的格林函数的理论与方法。将震源处理分为:(1)计算单一位错点源单位脉冲引起的格林函数;(2)用有限断层模型将断面划分为若干子源并有时间延迟,叠加所有子源的格林函数。本以1994年美国北岭地震为例,用解析法和数值方法分别计算了LV3、PCD、MCN三个场地的格林函数,并进行了对比分析。  相似文献   

4.
一种新的三维大地电磁积分方程正演方法   总被引:4,自引:4,他引:0       下载免费PDF全文
采用规则六面体单元和并矢Green函数奇异积分等效积分技术,已有的大地电磁积分正演方法具有不能有效模拟地下复杂地质体和计算精度偏低的缺点.本文提出了一种新的三维大地电磁积分方程正演技术,即采用四面体单元、解析的并矢Green函数奇异积分表达式,达到既能模拟地下复杂异常体,又能有效提高已有积分方程法计算精度的目的.首先,采用四面体网格技术离散地下复杂异常体,获得四面体单元上的大地电磁积分方程.然后,利用针对四面体单元开发的新的奇异值积分的解析表达式,准确计算线性方程中的并矢Green函数的奇异积分,从而获得精确的线性方程.借助于PARDISO高性能并行直接求解器,实现了三维大地电磁问题的高精度求解.最后,基于国际标准3D-1模型和六棱柱模型,通过与其他方法结果的对比分析,验证了本文方法的正确性、处理高电导率对比度的能力(1000:1)和处理复杂模型的能力.  相似文献   

5.
瞬变电磁场时域格林函数解   总被引:3,自引:3,他引:0       下载免费PDF全文
近源时间域电磁场具有信号强、探测深度大、精度高等优点,但传统勘探电磁场理论中偶极子近似在近源会引起较大误差,导致这一优势的发挥受到了制约.开展直接时间域电磁场解析式研究,是解决这一问题的途径之一.本文提出在点电荷微元假设下,引入时域格林函数,求取瞬变电磁场时间域解析解.采用积分运算法,把电磁场阻尼波动方程的求解问题转化为求其格林函数积分形式解的问题;建立辅助路径解决奇点问题,利用复分析中的约当引理、留数定理和广义函数等理论和方法,推导计算出时间域格林函数的时空四重广义积分.得到达朗贝尔方程的直接时域格林函数精确解析式,与传统方法“比拟”出的公式具有相同的形式,验证了本文推导的时域格林函数解析公式的正确性;推导出扩散方程的直接时间域解析解.通过与时变点电荷源时间域的电磁响应近似表达式进行对比,得出本文所推导的公式计算精度较高的结论;建立了全空间回线源瞬变电磁场问题的直接时间域求解公式.为解决全场区瞬变电磁场精细探测直接时域解析问题提供了基础理论.  相似文献   

6.
格林函数的奇异性处理   总被引:1,自引:0,他引:1       下载免费PDF全文
格林函数在解决电磁场和直流电场问题的积分方程法和边界单元法等方法中有着广泛的应用.在应用格林函数时会遇到格林函数奇异性的问题,即其在源区域的内部,当场点和源点重合时,会碰到对格林函数为无穷大,积分为奇异积分.对于不同问题中遇到的不同形式的格林函数奇异性问题的各种处理方法进行了分析和评述.对于标量格林函数奇异性问题处理方法有:挖去法、级数展开法、绕开法和解析法等.对于并矢格林函数奇异性问题处理方法有:分量处理法、源并矢法和拟源并矢法等处理方法.通过对实际工作中所遇到的新方法新问题中的不同形式的格林函数的研究,提供了奇异性问题的处理的方法和途径.  相似文献   

7.
本研究提出一种计算同震垂线偏差变化的方法,给出关于四种独立震源的格林函数,并将其应用于重力卫星GRACE观测数据的解析研究. 方法:针对球形地球模型,依据目前为止作者提出的球形地球位错理论的基本研究思路,计算四种独立震源产生的位错Love数并组成相应的格林函数.  相似文献   

8.
植被电磁散射的半空间模型研究   总被引:1,自引:0,他引:1       下载免费PDF全文
本文给出了地面植被电磁散射的半空间模型研究方法.在以往的相关文献中,均采用自由空间格林函数求解单个叶片散射体的散射场,本文利用半空间并矢格林函数求出了单个散射体的散射场,然后运用Monte Carlo方法模拟生成地面植被层,得到了在半空间下植被层的电磁散射特性,计算出了单、双站雷达散射截面,并与自由空间下的散射场做了相互对比.结果表明,在半空间格林函数下得到的散射场更为适用于描述地面植被的电磁散射特性.  相似文献   

9.
地球表面质量负荷的静态响应   总被引:7,自引:2,他引:7       下载免费PDF全文
本文给出了单位点质量负荷作用在球状成层地球模型上的解。对较新的G-D1066A地球模型求出了n直到10000阶的负荷勒夫数,并利用这组负荷勒夫数计算了格林函数的值。简述了如何利用格林函数计算地球对表面任何负荷的响应。讨论了负荷潮研究在地球物理和海洋学方面的一些应用。  相似文献   

10.
不同地球模型对负荷潮汐改正的影响   总被引:5,自引:0,他引:5       下载免费PDF全文
负荷潮汐改正与地壳上地幔结构有密切关系,主要原因是负荷潮汐取决于高阶勒夫数。利用较新的1066A地球模型,对负荷勒夫数进行了重新推算,并依据负荷勒夫数组成重力、倾斜、应变及位移格林函数。发现在近区用1066A模型和经典的G-B模型求出的格林函数值差别较大。据此,推估出不同地球模型对负荷潮汐改正的影响,其中远区估算按解析方法进行,近区由数值法进行。  相似文献   

11.
The theory of elastic wave scattering is a fundamental concept in the study of elastic dynamics and wave motion,and the wave function expansion technique has been widely used in many subjects.To supply the essential tools for solving wave scattering problems induced by an eccentric source or multi-sources as well as multi-scatters,a whole-space transform formula of cylindrical wave functions is presented and its applicability to some simple cases is demonstrated in this study.The transforms of wave functions in cylindrical coordinates can be classifi ed into two basic types: interior transform and exterior transform,and the existing Graf’s addition theorem is only suitable for the former.By performing a new replacement between the two coordinates,the exterior transform formula is fi rst deduced.It is then combined with Graf’s addition theorem to establish a whole-space transform formula.By using the whole-space transform formula,the scattering solutions by the sources outside and inside a cylindrical cavity are constructed as examples of its application.The effectiveness and advantages of the whole-space transform formula is illustrated by comparison with the approximate model based on a large cycle method.The whole-space transform formula presented herein can be used to perform the transform between two different cylindrical coordinates in the whole space.In addition,its concept and principle are universal and can be further extended to establish the coordinate transform formula of wave functions in other coordinate systems.  相似文献   

12.
对已经提出的四种经常使用的包线函数模型进行了评价,认为对于评价单峰点包线函数,除了要求形式简单以外,其优劣可以用函数对峰点位置和高峰隆起程度的控制能力作为标准。借助结构动力试验中阻尼比的确定方法,提出用平稳度指数来评价单峰点包线函数平缓程度。据此对四种包线函数的平稳度进行了分析,还对某些包线函数提出了一些改进。通过分析指出如果把“金井清谱”和“Penzien谱”等模型,看作是时间域中包线,这两个模型的调节能力都很强,而且比较方便,其中“Penzien谱”模型可以用来表示具有两个峰点的包线。文中的比较分析结果可以用于地震动包线函数选择和模拟。  相似文献   

13.
转换点位置的计算是转换波资料处理中的一个关键问题. 本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法. 研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法. 结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置.  相似文献   

14.
利用全自动多功能三轴仪进行了粉土场地的液化后变形试验,探讨了粉粒含量对液化后应力应变关系的影响.试验时施加动加载使试样达到设定的液化程度,然后立即施加单调荷载,以模拟地震现场发生强震时初始液化到液化后变形的整个过程.研究发现,不同粉粒含量粉土场地发生液化后变形得到的液化后应力应变关系可采用推导的同一理论关系式表示,其系数的差异表明了粉粒含量对粉土液化后变形的影响.验证结果表明该关系式能较好地模拟粉粒含量对粉土场地液化后应力-应变的影响,为同一地区或相似地区粉粒含量变化的粉土场地液化评估提供依据.  相似文献   

15.
散射问题中柱面坐标波函数的全域变换公式   总被引:1,自引:0,他引:1  
为满足偏心源和多体散射理论研究发展的需求,阐明了柱面坐标下波函数变换的基本类型,指出了现有变换公式的局限性,提出了全域变换公式,给出了应用实例并讨论了应用前景。研究表明,不同柱面坐标下波函数变换可以分为内域和外域两种基本类型,现有Graf加法公式仅适用于内域变换问题。提出的新公式突破了现有公式的限制,可完成弹性波散射问题中不同圆柱坐标系下波函数全空间域上的变换。应用新公式,构造了圆形空穴的外线源散射和内线源散射的解答,通过实例与大圆弧近似方法进行了对比,表明了新公式的优势和有效性。波函数变换基本类型的划分、全域变换概念的引入和公式构建的原理具有普适性,可以推广到其它坐标系下波函数的坐标变换公式的建立中。  相似文献   

16.
Introduction Up to now,simulation of non-stationary ground motion processes in engineering usually is bygenerating stationary ground motion processes,then making them non-stationary by an enve-lope function(Amin and Ang,1968);and the envelope function is usually calculated based on sta-tistical analysis of ground motion records,and the stationary ground motion processes are simu-lated using filtering method,spectral representation method,or time domain method,etc.Because the envelope function …  相似文献   

17.
Most seismic reflection imaging methods are confronted with the difficulty of accurately knowing input velocity information. To eliminate this, we develop a special prestack depth migration technique which avoids the necessity of constructing a macro-velocity model. It is based upon the weighted Kirchhoff-type migration formula expressed in terms of model-independent stacking velocity and arrival angle. This formula is applied to synthetic sub-basaltic data. Numerical results show that the method can be used to successfully image beneath basalts.  相似文献   

18.
We give an overview of the methods designed for reconstructing close-to-harmonic functions from the magnetic field ΔT. The formula of Yu.P. Tafeev is refined. It is shown that this refined formula directly leads to the relation derived by V.M. Gordin and his colleagues that allows isolating the harmonic component in the function ΔT. V.N. Strakhov’s linearized representation of the function ΔT is immediately derived from the main approximate Tafeev formula for Q ΔT. The experience of using Strakhov’s ΔS function in the interpretation of the magnetic anomaly ΔT generated by the Krivoi Rog structure is described. It is noted that the problem of reconstructing the corresponding harmonic functions from the data of magnetic and gravity surveys has much in common. The specific features of measuring the magnetic field H and magnetic induction B in the material media are considered, and the physical interpretation of these fields is presented.  相似文献   

19.
地震信号瞬时频率的估算   总被引:6,自引:1,他引:5       下载免费PDF全文
本文比较了四类计算瞬时频率的方法.基于相位差分的方法对噪声敏感;基于零交叉点的方法是针对周期信号提出的,对非周期信号效果不理想,通过插值改善估算结果;基于FIR微分器的方法有三种,基于两点的FIR方法与标准公式较为接近,但会引入半个采样间隔的时移;基于三点的FIR方法可以消除时移,但计算结果会出现负频尖峰;Claerbout方法不仅会引入半个采样间隔的时移,而且估算的频率值可能超过Nyqusit频率.本文在Claerbout公式的基础上导出了一种估算瞬时频率的公式,该公式消除了时移,且负频尖峰也得到抑制;基于时频表示一阶矩的方法的计算结果较为平滑,同时具有一定的抗噪能力.最后用本文导出的公式和时频表示一阶矩两种方法计算了南海东北部一条地震剖面的瞬时频率,并将其结果与Seismic Unix中相关函数的计算结果做了比较,表明本文导出公式的计算结果较为理想.  相似文献   

20.
基于对数目标函数的跨孔雷达频域波形反演   总被引:2,自引:1,他引:1       下载免费PDF全文
波形反演在探地雷达领域的应用已有十余年历史,但绝大部分算例属于时间域波形反演.频率域波形反演由于能够灵活地选择迭代频率并可以使用不同类型的目标函数,因而更加多样化.本文的频率域波形反演基于时间域有限差分(FDTD)法,采用对数目标函数,可在每一次迭代过程中同时或者单独反演介电常数和电导率.文中详细推导了频率域波形反演的理论公式,给出对数目标函数下的梯度表达式,并使用离散傅氏变换(DFT)实现数据的时频变换,能够有效地减少大模型反演的内存需求.在后向残场源的时频域转换过程中,提出仅使用以当前频点为中心的一个窄带数据,可以消除高频无用信号的干扰,获得可靠的反演结果.为加速收敛,采用每迭代十次则反演频率跳跃一定频带宽度的反演策略.实验证明适当的频率跳跃能够在不降低分辨率的基础上有效地提高反演效率.通过两组不同情形下合成数据反演的分析对比,证明基于对数目标函数的波形反演结果准确可靠.最后,将该方法应用到一组实际数据,得到较好的反演结果.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号