长方体ΔT场及其梯度场无解析奇点理论表达式

首次明确指出文献中长方体ΔT场及其梯度场理论表达式在上半无源空间存在某些点的场值无法计算的解析“奇点”问题 .通过详细的理论推导 ,深入分析了出现解析“奇点”的原因 ,并在此基础上导出新的长方体ΔT场及其梯度场在上半无源空间无解析“奇点”的理论表达式 .对比模型正演计算结果表明 ,新导出的ΔT场及其梯度场无解析“奇点”理论表达式的正确     

起伏地形条件下长方体磁场无解析奇点表达式（英文）

已有的均匀磁化长方体磁场计算方法大都假设观测点位于上半无源空间,对于起伏地形条件而言,这些方法可能存在解析"奇点"。为此,本文基于地磁场基本理论,采用变量替换的积分方法,导出了改进的磁场表达式,详细讨论并有效解决了整个无源空间的所有奇点问题。相比前人的方法,其积分过程更自然、简单,最后的积分结果形式更加统一,并且不需要坐标变换即可求出无源区任意点处的磁场值,从而简化了正演过程。对比模型试验表明,新导出的磁场无解析"奇点"理论表达式是正确的,并能适应地形起伏的情况。     

基于BTTB矩阵的快速高精度三维磁场正演

本文改进了一种快速、高精度空间域三维正演算法,用来计算地下场源在水平观测面产生的磁异常ΔT场及其梯度场,以解决传统空间域正演计算效率低的问题.算法采用长方体对场源区域进行剖分,观测点与场源剖分单元体中心点在水平面上的投影重合.改进的算法具有以下三个特点:(1)采用无解析奇点的解析解公式计算磁异常,保证计算精度.(2)通...     

三度体（均质模型）位场波谱的正演计算

本文给出任意指向的均质直线段、多边形面和多面体的重磁异常谱的解析表达式。利用它们可以进行不规则三度体的重磁异常谱的数值计算。此外,尚导出任意指向的斜平行六面体重磁异常谱的解析表达式。这些表达式结构简单,易于计算,用作位场的正演反演计算都十分方便。     

均匀半空间瞬变电磁场直接时域响应数值分析

近源时域电磁场具有信号强、探测深度大和精度高等优点,但传统瞬变电磁场理论中偶极子近似在近源区会引起较大误差,推导瞬变电磁场直接时域解析式是解决这一问题的关键.本文在点电荷微元假设下通过时域格林函数,采用分离变量等方法推导出了上半空间一次有源波动场和反射波的时域解析式和下半空间二次无源波动场的时域解析式,结合均匀半空间瞬变电磁场的边界条件给出了均匀半空间瞬变电磁场的直接时域解析式,进而利用第一型曲线积分,通过沿回线源叠加推导出圆回线源在瞬变电磁场中的直接时域解析式.然后在半空间表面上,与传统的电偶极源假设下的表达式作了比较.数值结果表明两者在远源区的计算结果相差甚微,而近源区则存在很大误差.本文利用真正点元（点电荷）严密推导给出的均匀半空间表面上瞬变电磁场的直接时域解析式适用于全场区探测,克服了偶极子假设下只适用远场区的不足,为瞬变电磁法的进一步发展和实际勘探提供了新的理论基础.     

多学科前兆场联合反演孕震区力学参数的探讨

根据地震前兆位移场的解析表达式，进一步给出各种应变，倾斜，地下水和地电阻率的前兆场的解析表达式，讨论了根据多种前兆场观测资料反演潜在震源参数的可能性。由于在这些表达式中具有相同的未知数F，x,y,hα，β，充分利用上述各种前兆场观测数据，联合反演潜在震源参数应数应该是可行的，但是，对于前兆台站来说，由于受局部条件和仪器漂移等因素的影响，观测资料往往与面上观测结果不一致。在实际反演计算时，还有许多技术问题需要解决。     

二维位场在有源空间的解析延拓

位场的向下延拓对重、磁等位场资料的处理与解释具有重要意义,但现有的延拓方法在位场向下接近场源时会产生剧烈震荡,计算过程很不稳定,而且理论上都不能越过场源.本文在复数域里定义了一个具有位场和几何双重特性的新函数PFGR,提出了一套新的延拓方法,在理论上实现了位场在包含场源的下半空间里的解析延拓,经过模型试算验证,表明该方法可以延拓越过场源,准确确定场源奇点,计算过程稳定,对测量中的随机噪声也具有较好的抗干扰能力.     

无压含不层非稳定流诱发地壳形变的定量计算

基于无压含水层迟后疏干的Ｂｏｕｌｔｏｎ理论，导出地下水无压水层完整井非稳定流单井排放诱发的地表位移的解析表达式，讨论了公式的计算方案，给出了模拟计算实例。     

两相饱和介质半无限空间位错位移场

弹性位错理论是研究地震震源正演和反演的重要工具之一. 单相介质位错理论求得的地面变形场是稳定的、不变的，而实际地震后的地面变形场是变化的. 本文利用集中力作用下无限空间中两相饱和介质Green函数的柱函数表达式，通过弹性运动方程积分解的Helmholtz方程的Hansen矢量变换，再叠加自由表面影响场的方法求得两相饱和介质半无限空间位错位移场. 本文结果结合Nur的DD模式，在进水排水的假说下能解释地震后地面变形场的变化.     

用传递矩阵法求解多层层状弹性半空间轴对称动力问题的奇…

本文利用Ｌａｐｌａｃｅ－Ｈａｎｋｅｌ联合变换及传递矩阵技术，导出了多层层状弹性半空间轴对称动力问题在层间完全接触情况下奇异解的一般解析表达式。本文方法概念清晰、公式简洁、易于应用。文中给出了计算实例，所得结果与文献的计算结果吻合较好。     

地震孕育体源流变模型（一）-位移场

为了研究地震孕育过程，在弹性包体理论的基础上建立了流变介质的体源包体模型理论．本文基于弹性包体模型的位移场的解析解，运用流变力学理论中的对应原理，导出流变体介质中，包体在任意一点的3个方向分别产生的位移的粘弹性解析解.      

地震前兆形变场研究

假设在均匀弹性介质半空间中，存在一单一力源F^→，通过求解弹性静力学方程，得到地球表面位移解3个分量（u、v、w)的解析表达式。讨论了理论垂直形变场的基本特征。力的大小是决定形变幅度的基本因素，力越大，形变幅度亦越大，反之，则越小。力源深度h越大，形变幅度越小，影响范围越大，反之，形变幅度大，影响范围小。力源倾角为正时（向下），形变中心区下降，倾角为负（向上）时，形变中心区隆起。力的方位角是决定形变场图像形态的重要因素。利用大同地区（1984－1987年）和张北地区（1984－1991年）水准测量资料反演这二个地区震前力学参数，结果表明，力源中心与实际震中位置相差小于40km，力的方位角近于水平，理论形变图与实测形变图基本相似。     

论稳定电流场的入射、反射、透射间的关系

根据分界面上稳定电流场的边界条件，解稳定电流场Ｌａｐｌａｃｅ方程，给出椭球形分界面入射场、反射场及透射场的表达式，并导出圆柱形、球形分界面上入射角、反射角及透射角的关系式后，指出，只在少数特殊情况下三者才相等，对一般形状的大多数光滑连续分界面情况，它们是不相等的．其值取决于分界面的具体形状．     

用传递矩阵法求解多层层状弹性半空间轴对称动力问题的奇异解

本文利用laplace—Hankel联合变换及传递矩阵技术,导出了多层层状弹性半空间轴对称动力问题在层间完全接触情况下奇异解的一般解析表达式。本文方法概念清晰、公式简洁、易于应用。文中给出了计算实例,所得结果与文献的计算结果吻合较好。     

点源两层大地三维地电模型视电阻率边界元解

本文用两层均匀大地人工点源作用下的空间电位分布作为基本解，导出了下伏基岩中陷含三维地质体情况下，关于电位的边界积分方程组，用边界单元法解积分方程求得地表电位，从而求得了一定装置下的视电阻率。本文方法的独特优点是积分只须在异常体边界面上进行了测线可以沿任意方向。用联剖法计算了均匀半空间下含球形异常体的ρｓ曲线，所得结果与其他学者用解析方法所求得的结果一致。用中间梯度法计算下层基岩中含有一长方体，测线     

磁力测深公式的推导及应用

本文从地球物理场出发，通过严密的数学物理推导，提出了磁力测深公式，并导出了著名的重力测深公式的更一般表达式。为了得到更为精确的解，本文又提出了应用最优化理论的空间域迭代算法，并根据地球物理模型对于高斯优化法进行了改进，在算法上有所创新，使其在计算机上实现时算法简单，计算速度快。文中还列举了一些模型的实验结果，对磁力测深公式及迭代算法的有效性及可靠性进行了论证。用本文方法对实际资料的处理结果是令人满意的，证明该方法有一定的实用性。     

三度体（物性随深度变化模型）位场波谱的正演计算

本文将均质的任意二维、三维物体位场的波谱解析表达式的研究成果推广到变密度、变磁化强度的更一般的情形。对密度差随深度呈指数函数衰减或线性变化的模型,获得了任意倾斜多边形质量面、斜平行六面体以及一般的多面体等形体的重力谱的解析表达式。它们的结构与均质体相应表达式一样简单,易于计算。以上结果表明,在很一般的条件下,位场波谱具有指数函数和的形式。     

任意震源力作用下多层粘弹性半无限地基平面问题的解析解

本文利用Ｆｏｕｒｉｅｒ变换及矩阵传递技术，导出了在任意震源力作用下多层粘弹性半无限地基平面问题的一般解析表达式。本文的方法概念清晰、公式简便实用。文中给出了计算实例，所得结果与现有文献的数值计算结果一致。     

InSAR干涉形变场远程模拟系统实现与地震地表形变场模拟算例分析

利用InSAR(Interferometric SAR)干涉测量技术可以获得地表形变场(视线向, LOS), 将其与弹性半空间中一定断层模型模拟计算出的地表形变场进行正演分析, 以便获取发震断层的几何学和运动学参数, 是一种典型的前向模拟模式。 该计算模式由于模拟计算程序调试、 安装过程的复杂性, 使得在每台微机上安装此类程序不仅费时费力, 且浪费大量的计算资源。 文中首先介绍了实现InSAR干涉形变场模拟模型由单机推广到Internet/Intranet上进行远程计算最终建立“InSAR干涉形变场远程模拟系统”的过程, 并基于此系统平台对1997年玛尼M_S7.9地震的InSAR地表形变场进行了模拟试算、分析等。     

球面水平应变场无偏计算的实现与滤波

GNSS观测资料为研究地壳应变场提供了必要的依托,但如何较为完整或恰当地反映各种地壳应变信息仍存在缺陷.从连续变化的角度出发,在ITRF参考框架下水平运动场的基础上,借助于多核函数获得水平运动场的解析式,据此,利用球面上无偏差应变计算式给出应变场结果,再利用多核函数法进行空间滤波,最后给出了利用本方法进行应变计算的实际算例.该方法有如下特点:①适用于各种空间尺度的地壳应变场描述,并且全面而又客观;②根据不同的需求和实际资料情况可提取不同频谱及其以内的应变信息;③可对应变参量进行严密的精度评定;④由于获得了解析式,故可在此基础上进行更深入的分析与计算.