共查询到20条相似文献,搜索用时 31 毫秒
1.
本文提出一种适于层状介质内点源模型的近场理论地震图算法,采用的位移表达式由Has kell矩阵的分解形式和离散波数积分组成。文中还出一起计算若干不同深度震源的近场理论地震图时的简便算法。这样,由本文的点源模型模拟复杂地震震源过程时可免去大量的审复计算。为检验算法功能,文中还给出了与文献[12]的数值计算的对比结果。 相似文献
2.
我国地震地电场观测时的漏电干扰问题会影响到观测质量和地震监测效能。针对这一实际问题,本文提出了点源干扰源定位方法。该方法首先建立地表点源对地电场观测影响模型,根据模型计算点源对地电场观测的干扰幅度,然后求解干扰源相对于中心电极的位置。将定位方法应用于高邮地震台和汉王地震台实际干扰源排查工作中,定位结果与实际干扰源位置基本一致,测区外干扰源定位误差小于2%,测区内干扰源定位误差小于20 m,并对影响计算误差的因素进行了分析。两个台站的实际应用表明:本文提出的定位方法为地电观测实践中查找点源漏电干扰源提供解决方案,能够有效地定位点源干扰源位置,提高干扰排查效率。 相似文献
3.
4.
二维层状介质点电源边界单元解 总被引:1,自引:0,他引:1
本文用付氏变换将点源二维地电断面的三维边值问题变成二维边值问题,然后讨论了边界元法在解决层状介质中存在二维柱体时点源场视电阻率计算问题中的应用。由于选取的基本解避免了对水平层状界面的积分计算,从而大大减少了边界剖分节点数。在计算中设法略去了对边界上位场法向导数的求解,使未知量大为减少,从而便利了求解二维边值问题,并由付氏反变换求得三维空间中的位。本方法所需计算机内存少,原始数据的输入简单,且有较高的计算精度。 相似文献
5.
从分层均匀介质中地震波传播的基本理论出发, 参考已有的利用广义反射透射系数矩阵方法计算分层均匀介质中点源引起的地面运动的方法和程序, 考虑了任意几何特征的有限震源及有限震源的震源机制随时间和空间发生变化的可能性, 并考虑了并行计算的发展方向, 本文对点源情况下的计算流程进行了改进, 重新设计编写了计算程序GRTMatSyn, 使其不但适应于点源的情形, 也适应于任意几何形状的、 有限的、 震源机制随时间和空间变化的有限震源的情形. 为了验证该程序的可靠性和计算效率, 设计了必要的有限震源模型和观测点分布, 计算了地面运动, 并与已有的被广泛认可的计算有限平面断层引起的地面运动的程序CompSyn的计算结果进行了对比分析. 结果表明, GRTMatSyn的计算结果可靠、 计算效率较高. 相似文献
6.
本文从一阶方程组形式的波动方程出发,发展了一种计算二维不均匀介质中点源P-SV波响应的近似方法。该方法通过引入线分布的应力作为震源,利用二维有限差分方法计算出线源响应,然后再经过波形校正和几何扩散校正得出相应的近似点源响应。通过把波形和振幅与精确解比较表明,该方法具有较好的精度。由于有限差分方法对于介质中速度和密度的分布没有特殊要求,另一方面,本文所给出的震源可以适用于位错点源、爆炸源或集中力源,因此上述方法十分适合于研究横向不均匀介质中的近场强地运动、爆炸振动或地震勘探等问题。 相似文献
7.
点源三维地电场的求解是一大型数值计算问题.本文用有限差分方法求解,最后形成一个线性方程组Ax=b,这里A是大型稀疏的带状对称矩阵.解大型稀疏方程组的完全Cholesky分解直接算法,一般要求巨大的机器内存来存储系数矩阵A,而且计算速度极慢.因此引入不完全Cholesky共轭梯度(ICCG)算法及按行索引的稀疏存储模式,充分利用A的稀疏性,使得计算速度大大提高,而内存要求则大大减少,因此ICCG算法是地电三维正演的强有力工具. 相似文献
8.
分析G-R-闵公式后得出,它是从点源幅射场导出的,与地震宏观场不相应。因在近场震源不能视作点源。在分析过程中,导出相当于点源、线源或面源以及复杂源的幅射场的深度公式。对实际的地震,它的源类型一无所知,所以不预作假定,将源指标几何扩散率n作为待定参数,导出了一个物理意义明确而又普遍化的震源深度公式。此公式的诸解法中,以计算方法准确、精度高,作图法有直观的优点,但准确性差、精度低,图算法只作获取粗略值和考察数据均匀性用。计算了9个8级以上巨震、1个71/2大震,和两个M_L=3的有感地震的震源深度,经对比结果很好,利用计算得出震源类型,结合宏观场研究了这些地震的震源几何学。并且首次在国际上给出宏观地震震源深度值的标准误差。 相似文献
9.
10.
11.
12.
13.
研究固体潮应力场时,需要考虑负荷潮应力场的影响,特别是在沿海地区,在地球的表层负荷潮应力场的大小甚至接近或超过固体潮应力场.根据Alterman方法,本文给出计算球状径向不均匀弹性地球模型对点源负荷的应力响应(负荷潮应力格林函数)的方法;作为计算实例,计算出GB地球模型的负荷潮应力格林函数.只要给出全球海潮模型,在球面上完成海潮潮高与负荷潮应力格林函数的褶积,即可计算出海潮负荷在地球模型表面上和内部任意点产生的负荷潮应力. 相似文献
14.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7%左右. 相似文献
15.
非规则采样快速傅里叶变换(NFFT)主要用于快速计算非规则采样数据的频谱及重建.该方法为非规则采样数据频谱重建技术的核心算法.在实现NFFT算法时,高速度和高精度计算是其应用的前提和关键.本文针对二维NFFT计算效率,应用表驱动思路进行改进,将Gauss褶积算子由矩形改进为椭圆以减少计算量,将e指数计算改进为乘法以加快计算速度,并建表解决NFFT算法在地震资料处理中的应用问题.本文同时给出了非规则采样地震数据NFFT谱重建方法.最后本文给出算例验证提出方法的计算速度和精度,和非规则采样地震资料重建结果. 相似文献
16.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7%左右. 相似文献
17.
18.
井地电法是深部地质矿产勘查的重要方法之一,研究高效率、高精度、适用于起伏地形等复杂条件下的正演模拟算法具有理论和实际意义.本文研究了点源场井地电位测量的三维有限元正演模拟,采用与井轴一致倾斜的计算区域,使剖分单元随地形起伏,用预处理共轭梯度法求解线性方程组,利用行压缩存储和改进的行压缩存储方式来存储刚度矩阵和预处理矩阵,实现了起伏地形和倾斜井情形下的正演计算.此外,本文利用仿射坐标变换技术,给出了平行五面体单元精确的单元积分公式,这种处理技术同样可推广应用于其他平行多面体单元的单元分析,相比于采用等参变换和高斯数值积分,能极大地提高计算效率.文中构建了一些模型算例,其计算结果验证了程序的可靠性和高效性. 相似文献
19.