在位场的区域-剩余分离问题中维纳滤波器的应用(摘要)

在区域物探工作中,重磁位场的区域—剩余分离或分离区域场和局部场的问题,近几十年来,一直被推断解释人员作为一个关键的方法技术课题来研究和探讨的.40年代末至50年代,提出了徒手园滑的各种经验方法和形式不同的各种算术平均法,以及位场的高次导数法;60年代后期至70年代初,提出了所谓的“非线性滤波”法和所谓的“匹配滤波器”;此后,有人给出了一种二维最佳线性数字滤波器的设计原理(王继伦,1977);我们曾提出并在实际中应用了一种分离区域场和局部场的最小均方误差滤波器(1978),这种滤波器实质     

反演的可靠性

虽然最小平方反演是数据分析常用的工具,而非单性问题是不可避免的问题;不过,这个问题通过考查模型响应对参数估算的灵度敏可以进行分析,这种灵敏度分方法还可产生只满足某些分解(或稳定度)准则的极值解。产生这种解的两种非常相似的的方法为“边界弯曲(Edgehog)”法和“最大平方法”,这两种方法是Jackson(1973、1976)提出来的。最小平方反演测定“非单值性度”的这些方法只需要奇异值分解(SVD)解中计算的数据。“边界弯曲”和“最大平方”方法在数学上是相似的,“最大平方”估算在计算上比较简单。这两方法都说明,一种反演的可靠性取决于特殊的误差准则以及与最小平方解有关的雅科毕(Jacobian)矩阵的特征等两个方面。两种相似方法的相似性由合成和实际垂直地震剖面(VSP)的旅行时间反演得到了证明。反演问题的灵敏度分析还提供了解的可靠性的定量测量值。     

用磁异常的谱及其导数求场源埋深和磁矩

随着电算的普及,物探重磁数据的处理已广泛使用了数字滤波技术,旨在提取信息、压低干扰的一类线性滤波,例如:根据信号实际输出与希望输出的方差最小为原则设计的维纳滤波和卡曼滤波;根据输出信号的瞬时功率与干扰的平均功率之比最大为原则设计的最大讯噪比滤波(D.O.诺斯,当干扰为白噪声时称为匹配滤波);根据:①在信号失真最小的同时将干扰(包括随机干扰和规则干扰)压低到一定程度;②在信号保持一定程度失真的同时,将干扰压低到最小为原则设计的最佳线性滤波(B.H.斯特拉霍夫-空间域,王继伦-频率域);还有根据反馈原理,采用Z-变换设计的旨在节省计算时间和减少边缘损失的递归滤波(B.K.布哈塔恰赖阿),等等.上述滤波只有事先知道输入信号和干扰的谱的统计表达式时,才能设计出一个可供实用的数字滤波器.但在实际问题中,这些谱的准确表达式当然是不知道的.常用的处理办法是将场源近似看作是形状规则的物体,只考虑场源埋深h和     

沃希变换在数字滤波中的应用

除众所周知的富氏变换外,滤波处理用其它变换也能实现。本文对使用一维和二维离散沃希交换(DWTS)进行一维和二维数字滤波处理进行了专门研究.特别对最佳滤波器问题(即用一维和二维离散数据其所需计算量明显减少)进行了研究.所谓最佳沃希滤波器是在修正最小均方误差准则的基础上确定的。本文还对被滤波的两种基本类型的信号(周期信号和平稳随机过程)进行了考虑.另在文中给出了两个实例。     

Radon变换的MATLAB实现

Ｒａｄｏｎ变换（包括线性、抛物Ｒａｄｏｎ变换）是地震数据处理中的一种强有力工具。作者在文中阐述了其原理及最小平方算法,并给出了ＭＡＴＬＡＢ语言缩写的源码,试算结果表明了该算法的有效性和程序的正确性,该程序可直接当作工具来使用。     

自适应网格滤波器在消除面波中的应用

根据格点滤波理论,提出了一种自适应最小平方网格算法（简称ＬＭＳ法）进行格点滤波参数的估算,实现此方法的关键在于连续自适应地估计偏相关系数,概述了ＬＭＳ滤波器的主要优点和基本原理,并将其用于理论模型和实际地震记录,其结果表明,ＬＭＳ自适应网格滤波器能非常有效地消除面波,此外ＬＭＳ自适应网格滤波器还能有效地压缩信号和增加信噪比,理论研究与实际资料处理表明,这种新方法具有简单,易行,快速,效果稳定的特点     

用最小二乘法设计电测深滤波器

现有的电测深滤波器还不能满足各种需要。对某些问题,或是在某些情况下,必须自行设计滤波器。有两种设计方法:付里叶变换法和最小二乘法。和付里叶变换法相比,最小二乘法的优点是:计算速度快,计算存储量小,而且不必考虑滤波器的位移。在不能用付里叶变换法的某些情况下,只要有一台袖珍计算器或微型计算机,即可用最小二乘法来设计所需的电测深滤波器。我们参考了O.Koefoed的程序,编制了一个设计电测深滤波器的PC-1500程序。它可用于设计各种装置,不同抽样间距的电测深正、反演滤波器,还可用来设计通用的最小二乘滤波器。     

深拖平缆采集资料中鬼波衰减方法探讨

鬼波衰减是深拖平缆采集资料处理中非常重要的一步,也是实现宽频处理的关键步骤。为此,笔者分析了深拖平缆资料特点,推导了鬼波的时距曲线方程,在此基础上,介绍了线性Radon域最小平方反演方法衰减鬼波的基本原理,该方法首先通过最小平方方法求得无鬼波的Radon域数据,再应用反变换算子中的一次波算子进行反变换,从而实现了对鬼波的压制。采用该理论对模型数据与实际地震资料进行了试处理,结果表明,深拖平缆资料通过本文的鬼波衰减技术处理后,陷频点得到了补偿,有效频带得到了拓宽,达到了宽频处理的效果。     

最小一乘方线性回归计算新方法

最小一乘方线性回归计算新方法李松仕（福建省水文水资源勘测局）１引言线性回归模型的最小一乘法,指的是使下列目标函数达最小的参数估计方法：ｍｉｎ｛Ｆ＝∑ｎｉ＝１｜εｉ｜｝其中,ε为点对回归面（或线）的离差。在统计数学中,称其为“稳健回归”。［１］对于一般...     

提高核磁测井差分谱油气识别方法的研究

对提高核磁共振测井差分谱识刖油气的方法进行了研究.介绍了核磁共振测井求解差分谱的原理,井用NUMAR公司的P型核磁处理软件DPP时P型核磁的原始数据进行了差分谱处理;根据限制最小平方误差法,反演得到T2谱,并对P型核磁测井原始数据进行了处理,与DPP处理的结果进行了对比,效果比较理想.通过比较分析,限制最小平方误差法对于提高低信噪比的油气显示很有意义,由于MAP-Ⅱ算法在低信噪比条件下反演的差分谱偏离真实较远.给实际应用带来较大的不确定性.文章提出的方法为低信噪比的数据处理问题提出了一种可能的解决手段.     

关于振动锤设计理论中几个问题的探讨

在工程地质勘察用的多功能钻机的振动锤的设计中,基本参数的确定是一个关键问题,本文通过对振动锤力学模型和振动方程的分析,探讨了在保证振动锤有一定寿命、有最佳工作效果、受外界因素影响最小的最佳工况时,基本参数(i、ε、λ、R)的组合规律。从振体处于一定位置时弹簧压缩量考虑,探讨了振动锤弹簧设计中,上、下弹簧刚度的分配问题。     

利用非线性模拟预测磁暴的发生

众所周知，磁暴与具有南向方向的星际磁场密切相关。这一结论是通过对线性预测滤波器的研究得出的，研究中采用的资料是Ｄｓｔ和ＩＭＦ数据系统，然而，从一次磁暴到另一次磁暴，该方法都要求对线性滤波器作较大的改动。因此在预测磁暴的活动中，线性滤波器法就显得具有一定的局限性。     

应用数字滤波消除潮位短周期波动对确定海平面变化趋势的影响

由一元线性回归分析的正则方程计算表明,在验潮记录较短的客观条件下,验潮序列中的周期波动对确定海平面变化趋势的影响是显着的.着重讨论采用数字滤波方法消除周期小于4年的短周期波动对确定海平面变化趋势的影响,并设计了最平滤波器,其通带截止频率对应的周期为4年,即应用最平滤波器对月均验潮序列进行低通滤波,基本上可消除周期小于4年的短周期波动对确定海平面变化趋势的影响.     

序贯判别及其应用的一些问题

一、费歇准则下两组线性判别及贝叶斯准则下多组线性判别的简单回顾假定有A、B两类已知总体,样本用P个变量描述。现在要判断未知样本归属A、B中哪一类。我们惯常使用的费歇准则下两组线性判别的做法是将P个变量线性组合成为一个判别函数:(1) K是变量号;λ_K是待定系数,它的确定原则是要使A、B两总体R平均值之差((?)_A-(?)_B)最大,而两总体内部R值的离差平方和最小。即差方比:     

基于Chebyshev走时逼近的三维多次反射射线计算

在利用有限差分等基于网格的数值分析方法解地震波走时所满足的程函方程时,由于速度模型的网格化离散等原因,会使走时在各网格节点之间不具有计算射线路径所要求的光滑性,即走时在邻近网格节点之间不具有连续的一阶导数。因此,直接利用网格节点走时计算射线路径会使最终的射线路径不光滑。为解决这个问题,已有研究者提出了基于B样条插值的逆向梯度方案（法）。然而,在速度发生突变时,B样条逆向梯度法所计算出的射线路径会具有较大的误差。针对这个问题,首先采用适合于解最小零偏差逼近及最佳平方逼近问题的Chebyshev多项式取代B样条对来自于分区多级计算方案的网格节点走时进行最佳逼近,得到在最小平方意义下的最优走时公式;然后采用与B样条逆向梯度法类似的计算过程得到光滑的射线路径。数值实验表明,利用Chebyshev多项式逼近走时可以得到具有很高精度的多次反射射线路径,在多次波偏移成像研究中具有潜在的价值。     

有界域上涉及临界增长含任意特征值的一类共振问题

给出了半线性椭圆方程-△u=λu|u|^2*-2u-f(x,u),λ∈[λ1,λk]；-△u=λu-γ（x,u) h(x),limu→∞γ（x，u）/|u|= ∞ λ∈[λ1,λk]的Dirichent问题的非平凡解的存在性定理等。     

起伏地形条件下二维场分量间的转换与延拓

对于起伏地形条件下直接测得的磁场垂直分量△Z,为了消除地形起伏的影响,便于解释推断,很有必要进行“曲化平”处理.已有许多物探工作者对“曲化平”和场的转换方法做了研究,其中包括冶金部物探公司研究所提出的“曲化平方法”.该方法不但精度较高,而且既能完成向上的“曲化平”计算,也能进行向下的“曲化平”计算.我们在文献的基础上,着重研究了起     

平方慢度射线追踪与非线性地震走时层析成像

重建地下速度分布是一个非线性问题。文章采用平方慢度模型解析计算地震射线传播路径，就模型参数值在网格节点上的情形，可形成一个非负矩阵，这类似单元模型的路径长度矩阵。利用非负矩阵的性质，通过对特殊广义本征问题的分析，引出了一般最小二乘解的“自然权“矩阵，并将正则化解约束在总走时为常数的超平面上。     

计算排序法及其应用

在各种数值计算过程中,数值的排序是一项基本的工作。排序速度的快慢往往直接影响到计算速度。目前各种各样的排序方法不下数十种,但它们都是建立在逐个比较、交换的基础上。常用的线性分类法是一种最原始的排序方法,其排序过程是将全部数据逐个进行比较,找出最小者(或最大者),再找出次最小者,依次进行下去,最后将全部数据按大小顺序排列出来。这种方法原理简单,容易实现。但这种方法的排序时间与数据个数的平方成正比例增长。如有N个数据,就需比较(N-1)·N/2次,同时还     

切面截距加权平均网格化算法

本文介绍了切面截距加权平均网格化算法的方法原理.这种算法能较好地反映随机点的数据的趋势特征,适用广,计算速度也较快.文中通过实例将它与按距离平方倒数加权平均法和局域加权最小二乘拟合法作了比较。