弹性波场和稳恒电流场的Frchet导数计算方法

在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的偏导数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度△Ｇ求得。把任意场源正演问题的解和Ｆｅｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解。文中以稳恒电流场问题为例，给出了直流电法视电阻率偏导数的数值结果。     

井间电阻率成像数值模拟

本文根据稳恒电流场的Ｆｒｅｃｈｅｔ导数解析式，给出一种井间电阻率成像的非线性反演算法。由于该导数完全由格林函数Ｇ及其梯度Ｇ表示，因此格林函数的计算是该方法的重要一环。     

井间电阻率成像中Frechet导数的算法比较

在用有限元法作电阻率反演计算中，求Ｆｒéｃｈｅｔ导数很重要。本文采用２．５维有限元法，对井间电阻率成像中的二极（ＡＭ）装置，分别用标准方法，格林函数法和扰动法计算了Ｆｒéｃｈｅｔ导数数值解，并对三种算法作了比较     

上海张堰铜矿外围地球物理场特征与找矿前景分析

上海张堰铜矿外围地球物理场特征与找矿前景分析杨达忠华东有色地勘局８１４队，镇江，２１２００４ＧｅｏｐｈｙｓｉｃａｌＦｉｅｌｄＣｈａｒａｃｔｅｒｓａｎｄＰｒｏｓｐｅｃｔｉｎｇＡｎａｌｙｓｉｓｉｎＯｕｔ－ｒｉｎｇＲｅｇｉｏｎｏｆＺｈａｎｄｙａｎＣｏｐｐｅ...     

基于波场延拓的多次波压制技术及其改进

水层多次波是海上地震勘探过程中的一个重要干扰因素。为压制水层多次波,本文在倾斜模型条件下采用基本的正向延拓方法去除多次波,取得了较好的压制效果。在此基础上实现了弯曲水底情况下正、反向波场延拓方法压制多次波;经对比发现,反向延拓方法克服了复杂水下介质情况下对格林函数求取的干扰,消除多次波效果更好。     

江苏宁镇地区化探采样点多元场分类结果

江苏宁镇地区化探采样点多元场分类结果何春林江苏省地矿局地质矿产调查研究所，南京，２１１１３５ＣｌａｓｓｉｆｉｃａｔｉｏｎＲｅｓｕｌｔｓｏｆＭｕｌｔｉ－ｅｌｅｍｅｎｔＦｉｅｌｄｏｆＧｅｏｃｈｅｍｉｃａｌＰｒｏｓｐｅｃｔｉｎｇＳａｍｐｌｉｎｇＰｏｉｎｔｓ...     

Chebyshev低通滤波器在位场数据波数域导数换算中的应用

导数换算作为常用的位场数据处理手段具有重要的物理意义。由于实际测量误差和噪声的影响,波数域导数换算存在计算过程不稳定、求导结果精度差的问题。为了降低噪声干扰,提高导数计算精度,提出了在波数域常规导数算子基础上附加Chebyshev低通滤波器的任意阶导数换算方法。通过分析该滤波器的滤波特性,并结合位场异常径向平均功率谱曲线特征,确定了Chebyshev低通滤波法的滤波参数,减少了人为因素对求导结果的干扰。二度体和三度体模型试验表明,与其他3种方法(常规FFT(fast Fourier transformation)法、向上延拓法和ISVD(integrated second vertical derivative)法)相比,该方法计算的导数与理论值的均方根误差最小,求导结果精度较高。在虎林盆地布格重力异常数据处理中,利用该方法计算的垂向二阶导数受噪声干扰小,结果可靠性较高,依此划分出的8条大断裂和11条小断裂在以往研究中均得到证实。     

二维位场向上延拓与向下延拓的样条函数法

本文用样条函数法实现了二维位场的向上延拓和向下延拓，向上延拓用样条函数法计算泊松积分，向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍｔｍ积分方程。模型和实例计算表明，本文方法正确可行且精度高。     

位场垂向高阶导数的Tikhonov正则化迭代法

在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。     

二维横向各向同性介质中的格林函数积分法合成记录

各向异性介质中地震波场计算的研究仍较少。本文根据Ｆ一Ｋ域格林函数的解析解，采用玻恩近似及积分法，计算了二维层状、横向各向同性介质中的合成记录。     

基于格林函数的层间多次波预测

对于海洋地震资料的多次波的预测通常会采用SRME(surface related multiple elimination)方法,但是对于地震资料中的强反射层之间形成的层间多次波,SRME方法则无法实现预测。本文基于惠更斯原理,利用采集接收到的已知地震波场通过格林函数的方法来推测在强反射层之间传播的层间多次波场。首先计算各接收点在强反射层之间传播的格林函数,将接收到的波场(即当前炮集的所有接收点接收到的原始地震数据)与格林函数褶积,将得到的结果相加作为当前接收点的层间多次波的模型。进而再从原始数据中利用自适应减法减去该多次波模型,这样便完成层间多次波的预测和衰减。通过模拟数据和实际数据测试,该方法能够较好地衰减层间多次波能量。     

三维方位各向异性介质数值模拟及波场分析

采用交错网格技术将速度-应力方程中的速度对时间的导数转化为应力对空间的导数, 将弹性波动方程表示为与二阶双曲方程等价的一阶应力-速度公式, 以实现三维三分量地震波场模拟。对获得的波场快照、VSP记录图像和地面记录图像进行波场分析, 发现在波场图像中存在明显的拟P波、拟快横波和慢横波, 还出现了横波分裂、横波分裂盲区及波面三分叉等地震波在方位各向异性介质传播时产生的特殊现象。研究结果表明, 用交错网格三维模拟方法研究方位各向异性介质对方位各向异性介质中地震波传播反演能起到很强的辅助作用。     

三维直流电场积分方程中奇异性的近似处理

直流电场积分方程的核函数是磁并矢格林函数,其数学表达式与电并矢格林函数的数学表达式完全不同。因此,在处理直流电场积分方程的奇异性时不能直接利用文献中针对电并矢格林函数所提出的奇异性消除公式。为了寻求处理磁并矢格林函数奇异性的有效途径,参考文献中针对电并矢格林函数的奇异性消除方法,提出了针对磁并矢格林函数的拟源并矢概念,并求出了当包围奇异点的小邻域为球体、立方体等不同形状时的拟源并矢。如果将这些拟源并矢代入到电场的积分方程中,可以得到只含有正常非奇异积分的数值计算方案。将这个计算方案用于实现关于直流电场的拟解析近似理论,则可以使三维直流电场的快速数值模拟成为可能。     

2.5维各向异性介质中多波波场正演模拟与波场分析

完全3维弹性波数值模拟计算时间长,并且占用庞大的计算资源,这不利于在计算机配置不高的情况下进行科学研究,而二维弹性波数值模拟又达不到三维模拟的精度;同时,当模型、波场空间分布比较复杂时,传统的3维波动方程拟谱法模拟结果比较差.因此,在较高数值精度的一阶应力-速度弹性波动方程的基础上,采用傅氏变换仅计算y方向的偏导数,利用有限差分方法计算x、z方向和时间的偏导数,即利用2.5维数值模拟方法,实现在二维介质中计算三维弹性波场.最后通过数值模拟实现了在各向异性介质中多波波场的数值模拟,验证了2.5维方法是一种高精度、高效率、且能适应复杂模型的正演模拟方法,通过波场分析进一步认识了波在各向异性介质中的传播规律.     

临界滑动场GCSF法分析：以新桥硫铁矿为例

临界滑动场法是边坡稳定性计算新方法,它可迅速准确地计算临界滑动面与临界滑动场以及由各出口对应的局部临界滑动面组成的全局临界滑动场。本文将临界滑动场方法应用于新桥硫铁矿下盘边坡稳定性研究,并针对其他质特性与开采条件作进一步改进,计算产生的边坡全局临界滑动场ＧＣＳＦ可直接用于边坡设计。     

结合有限项多指数级数近似的点源格林函数闭式算法

针对水平层状介质点电流源格林函数,提出了一种快速计算技术。该方法利用点电流源层状介质格林函数谱为单调实函数的特点,采用多指数级数近似格林函数谱,进而借助Somerfeld-Lipschtz积分实现格林函数的闭式计算。计算过程中,格林函数谱以等对数间隔抽样,近似多指数级数的指数常数保持固定并被等对数间隔抽样,近似过程只是指数级数系数谱的计算过程。计算结果表明,采用如下参数能实现点电流源格林函数高精度闭式计算:(1)λmin=10-6,λmax=10,采样点N=128,λ对数抽样;(2)T2min=10-7,T2max=103,T2对数抽样,抽样个数M=64;(3)最优正则化参数α2=10-5,采用最简单模型对解进行约束。     

论早石炭世石燕科13个属的分类位置-地层意义及地理分布

文章对早石炭世石燕科Ｆｕｓｅｌｌａ，Ｕｎｉｓｐｉｒｉｆｅｒ，Ｐｒｏｓｐｉｒａ，ｉｍｂｒｅｘｉａ，Ｆｉｎｏｓｐｉｒｉｆｅｒ，Ｃｅｌｓｉｆｏｍｉｘ，Ｇｒａｎｄｉｓｐｉｒｉｆｅｒ，Ｗａｒｓａｗｉａ，Ｇｙｐｏｓｐｉｒｉｆｅｒ，Ｆｒｅｃｈｅｌｌａ，Ｎｅｏｓｐｉｒｉｆｅｒ，Ｂｒａｃｈｙｔｈｙｒｉｎａ，Ｇｕｉｚｈｏｕｓｐｉｒｉｆｅｒ１３个属的分类位置、相互之间的区别特征和演变规律进行了讨论；在此基础上，进而讨论     

最优化广义离散Shannon奇异核交错网格褶积微分算子地震波场模拟

笔者基于离散Shannon奇异核褶积微分算子(GDSCD)计算地震波速度应力方程的空间导数,推导了一阶GDSCD的具体形式,并提出了优化方法,即在频率域逼近平面波的真实导数,得到了不同半径和采样下限的最优窗函数系数,通过滤波响应分析算子精度,与多种数值方法对比以及模型测试表明,笔者构造的最优化GDSCD模拟地震波具有较高的计算效率和精度.     

二维位场向上延拓和向下延拓的样条函数法

本文用样条函数法实现了二维位场的向上延拓和向下延拓,向上延拓用样条函数法计算泊松积分,向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍ积分方程。模型和实例计算表明,本文方法正确可行且精度高。     

位场各阶垂向导数的ISVD算法及其应用

位场垂向导数广泛应用于位场数据处理与解释当中,目前求取位场各阶垂向导数方法大致可以分为两类：①在频率域进行;②在空间域进行。针对这两类方法存在着稳定性较差的问题,这里介绍了计算位场各阶垂向导数的ISVD算法,并通过对单一地质体模型和地质体组合模型产生的布格重力异常进行处理分析,证明了ISVD算法在计算位场各阶垂向导数时,特别是在计算位场高阶垂向导数时,比在频率域及空间域计算结果具有较好的稳定性,ISVD算法能够有效地识别出浅层小尺度地质体边界。为了验证该方法对实际资料的处理效果,对某区块实际布格重力异常数据应用ISVD算法试算各阶垂向导数,结合地质背景等资料划分出断裂体系。结果证明,应用ISVD算法计算的位场各阶垂向导数,对于确定浅层小尺度地质体边界有较好的效果。