首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 156 毫秒
1.
电磁波法探测地下异常体的数值模型和实验模型   总被引:1,自引:0,他引:1       下载免费PDF全文
本文研究均匀有耗媒质全空间及半空间中三维异常体散射的数值模型和实验模型。用矩量法解体积分方程,导出自小块积分的简单准确的计算公式,说明非自小块积分以及索末菲积分的简化计算方法。给出两个利用金属镜象板的水池测试装置,说明与测试方法有关的若干主要问题。计算结果与测量结果基本吻合。  相似文献   

2.
利用积分方程法计算双侧向测井仪的井眼校正曲线   总被引:5,自引:0,他引:5       下载免费PDF全文
对电测井积分方程中的各积分项进行了理论推导,利用任意闭合曲面对曲面外一点所张立体角为零的特点将电势积分项由对柱体侧面的积分转化为对柱体底面的积分,从而使在纵向采用线性插值条件下给出了电势积分项的解析解.对电流积分项进行了化简,使之由曲面积分化为一维数值积分.利用所得结果计算了双侧向测井仪的井眼校正曲线.  相似文献   

3.
文中推导出单轴各向异性媒质半空间上方垂直磁偶极源电磁场的精确解析表达式. 应用圆柱波函数的球面波展开式和超几何函数理论,场分量中的索末菲尔德型积分被表示成快速、绝对收敛的球面波函数系展开式;展开系数是以物性参数为复宗量的勒让德多项式. 该展开式数学物理意义明显,并且不受场点和源点的位置、媒质的物性参数和频率等条件的限制. 利用本文的结果可十分方便地计算和分析任意场点处的电磁场分布.  相似文献   

4.
用散射、叠加方法推导出两层大地的并矢格林函数。使用这些函数,含三维异常体的二层大地的边值问题转变为积分方程。使用矩量法,可解此积分方程。 使用偶极-偶极装置计算了激发极化和电阻率响应的几个数值结果,并对这一算法作了几方面检验。  相似文献   

5.
区域地震范围的宽频带理论地震图算法研究   总被引:1,自引:2,他引:1       下载免费PDF全文
本文研究一种计算区域地震范围宽频带理论地震图的快速算法.对于层状地球介质模型,使用广义反射、透射系数矩阵和离散波数积分方法计算理论地震图.为了减少计算时间,改进了Filon波数积分方法,并用它来计算波数积分.数值模拟表明,在计算震中距为300km的理论地震图时,计算速度提高了1倍.  相似文献   

6.
将空气及地面以下两层视为三层有耗散的线性媒质,求出第三层中的水平偶极子激发的电磁场在地表处的积分形式解,进行数值计算,分析其频谱特征。结果表明电磁辐射异常的优势频率为0~2Hz,与实际观测资料相吻合,在超低频或甚低频,有比较明显的电磁场异常。并讨论其对地震电磁辐射监测的指导意义。  相似文献   

7.
本文研究并建立了一种模拟各向异性地层中三维电性异常体电磁响应的积分方程算法.首先讨论了并矢Green函数及其相关积分的计算,将水平层状各向异性地层中的电场并矢Green函数分解成含有奇异项的直达波与非奇异的来自各个层界面的反射和透射波两个部分,再应用等效体积单元和表面积分技术对积分方程的奇异核进行离散化处理以便提高离散方程的精度.然后为了节省计算机内存以及计算时间,引入基于Krylov子空间的迭代算法求解积分方程的离散化矩阵方程.最后通过与现有文献中的结果作对比从而检验了所述算法的有效性,并结合具体算例考察分析了地层的各向异性对三维电性异常体电磁响应的影响特征和规律.  相似文献   

8.
水平n层导电介质中时谐垂直电偶极子的电磁场   总被引:5,自引:5,他引:0  
对于包含n层线性、均匀、各向同性的水平成层导电介质,本文应用汉克尔变换及其逆变换推导求得了位于顶层导电媒质中的时谐垂直电偶极子在该层产生的电磁场的表达通式.通过这些表达式能够很方便地求得包含任意多层导体半空间中的电磁场.文中最后计算了导电半体空间和二层导电介质中时谐垂直电偶极子产生的电磁场,验正了本文所得结果的正确性。  相似文献   

9.
二维线性流变体波的有限元模拟   总被引:4,自引:7,他引:4       下载免费PDF全文
用加里津方法推导了二维线性流变体内波传播的有限元方程,方程内包含了边界项,使得处理人工边界影响有较好的结果.同时考虑了其它一些因素,例如,吸收介质、爆炸源、自激自收等.实例表明,作者在仿真性方面的努力是成功的,计算的垂直地震剖面理论地震图质量比较好,结果也显示了地表低速层对地表水平记录剖面的干扰.  相似文献   

10.
为了实现大地电磁测深三维高效快速正演算法,本文研究了大地电磁测深积分方程三维正演模拟的理论与方法.首先详细地推导了Fredholm型积分方程;然后给出积分方程的离散方法、并矢格林函数积分的求解方法以及大型矩阵方程的存储与求解方法;最后,通过理论地电模型的计算检验了方法和程序编制的正确性,并展示了几种理论地电模型的数值计算结果,分析其异常特征.算例结果表明,该正演算法在CPU主频为1.81GHz、内存为1GB的计算机上能正确地模拟多个异常体的异常特征;地表电性不均匀体会使地表视电阻率曲线产生十分严重的畸变;TE模式相对TM模式纵向分辨率比较好,而TM模式的横向分辨率较好;随着异常体埋深的增加,反映异常的分辨能力进一步降低.  相似文献   

11.
An analytical approach is developed to study the dynamic response of a flexible plate on single-layered saturated soil. The analysis is based on Biot's two-phased theory of poroelasticity and also on the classical thin-plate theory. First, the governing differential equations for saturated soil are solved by the use of Hankel transform. The general solutions of the skeleton displacements, stresses, and pore pressures, derived in the transformed domain, are subsequently incorporated into the imposed boundary conditions, which leads to a set of dual integral equations describing the corresponding mixed boundary value problem. These governing integral equations are finally reduced to the Fredholm integral equations of the second kind and solved by standard numerical procedures. The accuracy of the present solution is validated via comparisons with existing solutions for an ideal elastic half-space. Furthermore, some numerical results are presented to show the influences of the layer depth, the plate flexibility, and the soil porosity on the dynamic compliances.  相似文献   

12.
区间B样条小波有限元GPR模拟双相随机混凝土介质   总被引:2,自引:2,他引:0       下载免费PDF全文
冯德山  王珣 《地球物理学报》2016,59(8):3098-3109
基于可分离小波理论,由一维区间B样条小波尺度函数的张量积构造二维B样条小波基,并将它作为GPR波动方程求解的插值函数,通过引入转换矩阵,实现小波系数空间与雷达电磁场之间的转换.应用Galerkin算法,推导了二维区间B样条小波有限元GPR波动方程离散格式,求出了2阶1尺度与2阶2尺度BSWI尺度函数的积分值及联系系数,给出了该算法的详细求解过程.编制了BSWI的Matlab模拟程序,应用该程序对两个典型实例进行了正演,结果表明:BSWI能采用较少的单元达到与FEM相似的精度,而BSWI算法尺度提升能提高解的精度,但耗时会急剧增加.最后,将BSWI算法应用于双相随机混凝土模型,说明随机介质模型理论能灵活、有效地描述实际混凝土介质的分布,正演剖面与实测剖面特征更相符,能更真实地模拟雷达波的传播过程,可为提高GPR的探测效果和解释准确性提供理论基础.  相似文献   

13.
Three-dimensional time-harmonic response of a poroelastic half space subjected to an arbitrary buried loading is investigated. The analysis starts with the field equations in cylindrical coordinates based on Biot's general theory of poroelasticity. General solutions for the displacements are first derived using the Fourier expansions and Hankel integral transform with respect to the circumferential and radial coordinates, respectively. The transformed-domain solutions are obtained in explicit form. The physical-domain displacements and stress components are then obtained numerically by inverse integral transform. Comparisons illustrating the accuracy of the developed approach are made with existing solutions for an elastic half space, which is reduced directly from the general solution developed in the paper. Numerical results are presented for the displacements of a saturated soil subjected to a horizontal internal excitation.  相似文献   

14.
Satellite missions CHAMP and GRACE dedicated to global mapping of the Earth’s gravity field yield accurate satellite-to-satellite tracking (SST) data used for recovery of global geopotential models usually in a form of a finite set of Stokes’s coefficients. The US-German Gravity Recovery And Climate Experiment (GRACE) yields SST data in both the high-low and low-low mode. Observed satellite positions and changes in the intersatellite range can be inverted through the Newtonian equation of motion into values of the unknown geopotential. The geopotential is usually approximated in observation equations by a truncated harmonic series with unknown coefficients. An alternative approach based on integral inversion of the SST data of type GRACE into discrete values of the geopotential at a geocentric sphere is discussed in this article. In this approach, observation equations have a form of Green’s surface integrals with scalar-valued integral kernels. Despite their higher complexity, the kernel functions exhibit features typical for other integral kernels used in geodesy for inversion of gravity field data. The two approaches are discussed and compared based on their relative advantages and intended applications. The combination of heterogeneous gravity data through integral equations is also outlined in the article. panovak@kma.zcu.cz  相似文献   

15.
在文献[1]中我们利用复变函数论的方法只解决了二维的曲面上的位场转换问题.在本文中,我们利用单层位势理论将它化为一个奇异积分方程求解的问题,并且提出了用奇点分离法来解此奇异积分方程,这对二维和三维情况都适用.通过理论模型的计算,结果的精度是令人满意的.  相似文献   

16.
This paper addresses the horizontal vibration of a rigid disk embedded in a poroelastic half-space in contact with a fluid half-space using the poroelastic theory of potentials. The solution of this problem is expressed in terms of dual integral equations that are converted into Fredholm integral equations of the second kind and solved numerically. Selected numerical results for the horizontal dynamic impedance coefficient are examined based on different poroelastic materials, embedment depths, and excitation frequencies; furthermore, the results are analyzed for the cases in which there is and is no fluid overlying the poroelastic medium to examine the effect of fluid. The results of this study are helpful for designing a foundation embedded in the seabed due to dynamic horizontal forces.  相似文献   

17.
Summary Green's theorem on harmonic functions makes it possible to determine the integral relationship between the harmonic function and its derivative with respect to the normal on a closed Lyapunov surface. The conditions of solvability are given by Fredholm's theory of integral equations. The solution for a sphere was presented by Molodenskii[3] and the general solution with the help of Molodenskii's parameter k by Ostach[4]. The present paper indicates a possibility of solving this problem with the help of a system of linear algebraic equations, a simplified modification of the Ostach-Molodenskii solution and, finally, a method, based on Eremeev's solution of the fundamental integral equation[5].  相似文献   

18.
Integral abutment bridges (IABs) are jointless structures without bearings or expansion joints which require minimum or zero maintenance. The barrier to the application of long‐span integral abutment bridges is the interaction of the abutment with the backfill soil during the thermal expansion and contraction of the bridge deck, that is, serviceability, or when the bridge is subjected to dynamic loads, such as earthquakes. The interaction of the bridge with the backfill leads to settlements and ratcheting of the soil behind the abutment and, as a result, the soil pressures acting on the abutment build up in the long term. This paper provides a solution for the aforementioned challenges by introducing a novel isolator that is a compressible inclusion of reused tyre‐derived aggregates placed between the bridge abutment and the backfill. The compressibility of typical tyre‐derived aggregates was measured by laboratory tests, and the compressible inclusion was designed accordingly. The compressible inclusion was then applied to a typical integral frame abutment model, which was subjected to static and dynamic loads representing in‐service and seismic loads correspondingly. The response of both the conventional and the isolated abutment was assessed based on the settlements of the backfill, the soil pressures and the actions of the abutment. The study of the isolated abutment showed that the achieved decoupling of the abutment from the backfill soil results in significant reductions of the settlements of the backfill and of the pressures acting on the abutment. Hence, the proposed research enables extending the length limits of integral frame bridges subjected to earthquake excitations. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

19.
本文比较详细地讨论和分析了均匀介质中零偏移距散射曲面拓扑奇性和相应回波解析奇性之间的关系和分类.其基础是衍射理论的Kirchhoff积分,原理是通过曲面几何函数的解析奇性和距离函数的拓扑奇性之间的联系.基本的数学工具是Morse理论和突变理论.  相似文献   

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

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