二维完全非均匀介质中位场格林函数的数值解

利用数值模式匹配理论,对具有轴对称的任意二维非均匀介质中位场的格林函数给出数值解。应用这一数值解,对复杂介质环境中双侧向测井响应进行了高效的数值分析。由于采用了半解析-半数值的混合方法,较之二维有限元分析具有很高的效率,节省大量计算时间。     

轴对称二维任意非均匀介质中偏轴点源激励的位场数值解

对具有任意多平面分层,且每一平面分层中又具有任意多径向柱面分层的复杂非均匀介质中偏轴点源激励的位场分布进行了高效数值分析．对于这种三维非均匀位场分布,应用数值本征模展开描述其径向的基模和高阶模,用集合反射矩阵和集合透射矩阵匹配各平面界面的边界条件,各平面分层中的位场则用递推算法获得．偏轴点源位场的求解,为电法测井中高分辨薄五层探测提供了理论分析基础．     

轴对称二维任意非均匀介质中偏轴点源激励的位场数值解

对具有任意多平面分层，且每一平面分层中又具有任意多径向柱面分层的复杂非均匀介质中偏轴点源激励的位场分布进行了高效数值分析．对于这种三维非均匀位场分布，应用数值本征模展开描述其径向的基模和高阶模，用集合反射矩阵和集合透射矩阵匹配各平面界面的边界条件，各平面分层中的位场则用递推算法获得．偏轴点源位场的求解，为电法测井中高分辨薄五层探测提供了理论分析基础．     

横向各向同性介质二维格林函数的F-K算法

采用算子表示法,在Ｆ-Ｋ域内将横向各向同性介质中的波动方程进行算子分解,然后利用留数定理求出格林函数的解析解。并给出了一个算例,得出时空域中ｑＰ波、ｑ（ＳＶ）波和ＱｓＨ波的格林函数剖面。     

横向各向同性介质二维格林函数的F-K算法

采用算子表示法，在Ｆ-Ｋ域内将横向各向同性介质中的波动方程进行算子分解，然后利用留数定理求出格林函数的解析解。并给出了一个算例，得出时空域中ｑＰ波、ｑ（ＳＶ）波和Ｑ_ｓＨ波的格林函数剖面。     

任意曲线上二维重磁位场转换的B样条函数法

重磁位场转换在重磁资料解释中是十分重要和必要的。随着物探仪器的改进和资料解释精度要求的提高,人们对位场转换方法的研究越来越深入。本文提出一种原理简单、计算精度高的新方法--B样条函数法,应用其插值、求导、求积的性质,以解决曲线上二维重磁位场转换的两类问题。一是Z_α分量向H_α分量的转换及向上延拓,二是水平和垂向导数的换算及磁位的计算。 样条函数为分段多项式函数。B样条函数是用σ函数逼近的一类样条函数,其节点     

横向非均匀介质远震体波接收函数的波场特征

系统研究了台站下方存在倾斜界面时,远震体波接收函数的波场特征及多方位接收函数对台站下方介质横向非均匀特征的反映能力．我们利用合成三维横向非均匀介质接收函数的Maslov方法,具体模拟了台站下方存在倾斜界面时的接收函数响应．结果表明,当台站下方存在倾向一致的倾斜界面时,远震体波接收函数的径向分量和切向分量分别随震源方位角的变化呈现对称和反对称的规律性变化．利用不同震源方位角接收函数径向分量和切向分量的变化规律,可以估计台站下方界面的倾向和倾角．当台站下方各界面的倾向不同时,随方位角的变化,接收函数只能直观给出界面的整体倾向．实际观测数据的分析结果表明,对于利用单个台站接收函数研究台站下方介质的横向非均匀特征来说,简单倾斜界面可以看作较好的一级近似．     

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

三维多值走时地震波场重建及格林函数计算

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

电测井Green函数的Taylor展开

本文证明了在分区均匀的介质中,电测井的Green函数是各区电导率的解析函数,并给出了它的Taylor展开式。如果背景值选择得适当,这个展开式是收敛的。     

均匀半空间导体的并矢格林函数解析式

导出了低频时谐电磁场中均匀半空间导体的并矢格林函数解析式,它由初等函数和修正贝塞尔函数的有限形状所表示．该解析式可用Kelvin函数的近似多项式计算,可显著提高并矢格林函数的计算速度．用"电磁分析方法验证模型"的第15号问题作为验证模型,结果表明,计算值和测量值吻合良好．     

电测井Green函数的Taylor展开

本文证明了在分区均匀的介质中,电测井的Green函数是各区电导率的解析函数,并给出了它的Taylor展开式。如果背景值选择得适当,这个展开式是收敛的。     

渗透率场敏感系数的数值计算

由渗流微分方程定解问题，利用格林互易定理从理论上导出了渗透率场敏感系数 的计算公式，并借助于数值积分和差分方法给出了渗透率场敏感系数的离散形式. 由Peacem an方程建立了井压渗透率场敏感系数与网格压力渗透率场敏感系数的关系. 对理论模型利用 三维不均匀非稳定渗流场的压强数值解计算了井压对渗透率场的敏感系数，并用敏感系数的 直接计算方法进行了验证，结果表明本文的计算方法与直接计算方法相吻合.     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。