应用计算机“推导”位场水平导数公式的方法及程序

根据台劳级数展开式，将窗口内网格结点位场数据在计算点展开，且对应点两两相减，由最小二乘法求得窗口内网格结点的水平导数算子，将此算子与观测的位场数据褶积，就可求得位场的水平导数。     

用B样条函数计算非等距点重磁异常垂向二次导数

位场垂向二次导数是基于拉氏方程求得的。1969年B.K.Bhattacharyya实现了二元三次样条函数计算位场导数的方法。本文采用了非等距结点B样条函数求导的方法,它具有灵活统一的优点。取作为二次微商计算式。由于用曲线插值求导取代曲面插值求导,程序更便于编写。在计算机上分别用三种模型作了试算,效果很好,精度比其它方法(包括上述二元三次样条函数法)提高了一个级次。     

位场高阶导数的波数域迭代法

位场导数换算是位场数据处理和解释的基本方法,但其不稳性限制了高阶导数在实际资料处理中的应用。针对于此,提出了位场高阶导数的波数域迭代法。在导数算子中引入一个包含求导阶次和高阶导数算子本身相关的低通滤波器,然后通过迭代法进行逐次逼近,进而得到高阶导数波数域迭代法的递推公式;同时证明了该方法的收敛性,对比分析了迭代法与常规方法导数算子的滤波特性。在高阶导数的模型检验和实例应用中,迭代法与常规方法的对比分析结果表明,迭代法具有稳定性强和计算精度高的优点。     

位场垂向导数的稳定算法

垂向导数广泛应用于位场数据处理和解释中,常规的空间域求导公式多以台劳级数为基础,这不可避免会引入截断误差,且不便于求取高阶导数。另一种应用较多的波数域算法。其求导算子在高波数成份产生振荡失真,很不稳定,在实际应用中很难取得理想结果。这里提出了一种计算位场数据垂向导数的稳定算法,在波数域求垂向积分,在空间域计算二个二阶水平方向导数,借助拉普拉斯方程,求取位场异常的各阶垂向导数。算法只在求垂向积分时引入了一次付里叶变换,其它计算全部在空间域中进行,从而有效地避免了吉布斯效应的重复累加传递,有效降低了常规波数域求导算子在高波数成份产生的振荡失真。通过理论模型对比,证明了这里所提出的算法较常规的波数域方法具有良好的稳定性,特别是在计算高阶垂向导数时有明显优势,能有效地提高与垂向导数有关的位场定量解释方法的准确性。     

二重网格方法用于曲面位场转换

本文介绍了二重网格方法及在求解曲面位场转换的积分方程中的应用。给出该积分方程的二重网格方法的具体算法和二重网格迭代算子的结构图。并提出了三点注意事项。     

地震波场模拟中的褶积微分算子法

褶积微分算子是将偏导数算子的带限傅氏谱进行反傅氏变换得到的,具有对称或反对称和衰减速度快的特点。根据频谱分析理论,选用序列在结束时能光滑趋于零的截断窗口将褶积算子截断,结合考虑微分算子守恒定律,将褶积算子中心值进行修正,得到适合于波场模拟的最优权系数。给出用褶积微分算子模拟声波波场的实现途径与实例。理论分析与模型试算证明了褶积微分算子法的正确性与有效性。     

基于窗口数据的重力相关成像

位场相关成像可以快速估计场源的空间位置。为了减小叠加异常对相关成像的影响,现有的成像方法通常是采用剩余异常或异常垂向导数成像。笔者通过加入一个数据截取窗口,计算窗口内数据的相关系数,以减少相邻场源间的影响,从而达到提高成像分辨率的目的;并通过相关系数极值的大小选取窗口,从理论上使成像结果更加可靠。通过重力模型试验,验证本方法能够提高成像分辨率。将该法应用于葫芦岛市黑鱼沟的某矿洞的深度计算中,取得较好效果。     

利用重力归一化总梯度及相位法研究断裂构造

用波数域的位场转换和快速傅立叶变换理论研究重力归一化总梯度和归一化相位;在波数域向下延拓和导数计算的滤波算子中,分别引入圆滑滤波因子,抑制了向下延拓计算中对噪声干扰的放大作用,增加了计算的稳定性;利用改进的方法,计算并阐述了4种断裂构造GH场等值线和相位曲线的变化特征;别列兹金重力归一化总梯度法中极大值理论在研究断裂构造中具有局限性,对于无限延伸断裂构造,应根据GH场等值线的走向、形状的特点,结合归一化相位曲线的转折点确定断裂空间位置.     

利用相关系数法确定三种位场垂向导数换算的滤波参数

垂向导数被广泛应用于位场数据的处理和解释中。在波数域中常规垂向导数算子受噪声的影响较大。通常解决方法是在垂向导数算子中引入低通滤波器,来压制高频噪声的干扰。Butterworth低通滤波、Hanning窗滤波以及向上延拓滤波是三种效果较好的低通滤波,但存在滤波参数较难选择的问题。针对这个问题,本文利用相关系数法确定各方法的滤波参数,避免了人为因素对求导结果的影响。模型试验和松辽盆地镇赉地区实际数据处理结果表明,使用相关系数法确定滤波参数后能快速获得较好的垂向导数结果,其中Butterworth低通滤波的导数结果相比其他两种方法的效果要好。     

基于泛克里格方法的位场扩边处理

在频率域进行位场数据处理与转换时,为了减弱边界效应以及Gibbs效应,需要对原始位场数据进行扩边处理,使得网格数据在x、y方向上的点数均为2的幂次方,并且扩边后的边界值均相等。目前比较常用的扩边算法是余弦扩边方法,为了寻求更高的处理精度,这里将泛克里格插值方法引入位场扩边处理中,在频率域计算了位场数据的一阶垂向导数,通过理论模型对比了该方法与传统的余弦扩边以及其他常用网格化方法的频率域处理精度,结果表明:泛克里格方法能很好地实现位场连续且平滑的扩边,并且与其他扩边方法相比有着更高的精度。     

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

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

井中磁测资料预处理与弱信号识别

井中磁测仪器距离场源近,磁测资料受到钻孔打穿磁性体产生的内磁场影响,往往存在强干扰压制有用信号,因此研究井中磁测资料预处理与弱信号识别方法具有重要意义。笔者将地面重磁资料边界识别中倾斜角总水平导数和解析信号振幅方法改进后应用于井中磁测资料预处理中,提出倾斜角垂向导数和均值归一解析信号振幅方法。通过模型分析证明了两种方法都可识别强弱不同的井旁磁性体位置,反映边界信息,且给出了计算均值归一解析信号振幅过程中窗口选择的依据。将两种方法应用于江苏某工区ZK002孔磁测资料预处理中,得出倾斜角垂向导数方法应用于实测资料中会放大干扰,而均值归一解析信号振幅方法增强了该钻孔弱异常信息,且压制了干扰信号,处理结果与磁化率测井和钻孔资料相对应。     

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

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

一种地震共反射点道集数据的叠前反演方法

叠前地震资料含有地层的纵波、横波速度和密度等信息。利用叠前反演获得隐藏在地震数据中的这些基本参数后,即可揭示大量岩性及孔隙流体性质的信息。这里推导了平面波在层状弹性介质中传播的正演算子,提出了一种基于Brent正交搜索方向组的叠前三参数反演方法．该方法不需要求解庞大而复杂的导数矩阵。通过自适应退火因子和罚函数来处理约束条件,提高了算法的稳定性。将K—L变换引入到方向置换过程,有效防止了搜索方向组的线性相关。经理论模型和油田实际数据的反演结果表明,该反演方法是一种利用叠前地震数据进行储层预测的有效手段。     

重磁数据稳定向下延拓的水平导数迭代法

向下延拓是重磁数据处理的常用手段,能有效地区分叠加异常、增强浅部异常,但现有向下延拓算法的计算大多是不稳定的,且易受噪声的干扰,往往造成异常形态的畸变.提出一种基于水平导数和向上延拓联合迭代的向下延拓算法,由于向上延拓和水平导数的计算是稳定的,因此该向下延拓方法可有效地增强结果的准确性和稳定性.理论模型试验表明该方法的向下延拓结果比Fourier变换计算结果更加稳定、准确,且受噪音干扰小.将该方法应用于实际数据的处理,结果显示该方法能稳定和准确地完成异常的向下延拓任务,且有效地增强了浅部局部异常.      

A new ionospheric tomographic algorithm — constrained multiplicative algebraic reconstruction technique (CMART)

For the limitation of the conventional multiplicative algebraic reconstruction technique (MART), a constrained MART (CMART) is proposed in this paper. In the new tomographic algorithm, a popular two-dimensional multi-point finite difference approximation of the second order Laplacian operator is used to smooth the electron density field. The feasibility and superiority of the new method are demonstrated by using the numerical simulation experiment. Finally, the CMART is used to reconstruct the regional electron density field by using the actual GNSS data under geomagnetic quiet and disturbed days. The available ionosonde data from Beijing station further validates the superiority of the new method.     

估计迹长概率分布函数的新方法及其应用

本文提出一种估计迹长概率分布函数的新方法。该方法的优点是不需要有节理在测窗上出露长度的数据，有效的减少了野外工作量