跨孔地震CT中的逐次线性化方法

非线性问题的逐次线性化方法的每个线性化过程可通过最小二乘ＱＲ分解算法求解．为了精确地描述和控制逐次线性化迭代过程，引入迭代阻尼系数．首先根据模型拟合表最小的准则确定满足最佳一次线性化效果的最佳初始迭代阻尼系数；然后随遂次线性化迭代而减小迭代阻尼系数，并加入松弛因子对每次迭代输出的波速修正量进行控制；最后根据混沌理论，采用Ｌｙａｐｕｎｏｖ指数和模型拟合差的变化，描述第ｋ次迭代所处的状态，在解估计的方差急增之前及时停止迭代，获得可信度和分辨率最佳的解估计．模型计算与实例证明方法有效．     

地震勘探时空域双曲滤波因子的特点与提高信噪比的能力

为改善强随机噪声背景中地震勘探资料信噪比（S/N）,提出了混沌振子算法;该算法可实现对未知同相轴的确定. 本文的时空域双曲滤波（hyperbolic time\|distance relation filter,HTDF）是上述算法的后续处理,即处理得到存在的弱同相轴.HTDF的理论基础是最小平方滤波和混沌振子算法. 为说明滤波过程与效果,对于包括一个反射同相轴的地震记录进行了多道滤波因子计算. 结果表明,与一个输出道相应的双曲滤波因子组具有复杂的时频域结构,突出的特点是振幅谱“波动式陷频”表现;弥补或消除这个不足的方法是与输出道相应的多个滤波因子之间存在的振幅谱“互补”性.利用相同随机噪声背景下的共炮点理论记录,经与域双曲速度滤波、褶积滤波处理方法比较,在振幅谱、子波初至、零相位性、S/N改善等方面,HTDF的滤波效果优于另外两种.最后探讨了该滤波技术的机理,包括同相轴内各子波间的相似性,以及对组内因子实施的多次有效相加处理.     

重力异常和重力梯度联合辅助导航算法及仿真

鉴于潜艇的惯性导航系统误差随时间积累的缺陷,重力无源辅助导航近几年成为水下潜艇辅助导航研究的热点和前沿问题.采用多模型自适应Kalman滤波并行算法对重力异常和重力梯度联合辅助导航,通过对沿航迹重力异常和重力梯度信息因子的实时分配,加权处理得到潜艇位置的最优实时估计.选取不同特征区域的2条航线进行仿真试验,对算法的匹配...     

一种改进的低纬度磁场化极方法——变频双向阻尼因子法

本文在考察分析原有的低纬度化极阻尼因子的基础上,提出了对原有阻尼因子的改造方法,即提出变频双向阻尼因子,改造后的阻尼因子不仅考虑了对振幅的压制,并兼顾了对相位谱的改造,同时还增加了频率控制因子.通过模型数据和实际资料检验,表明在低纬度地区,特别是在极低纬度地区,变频双向阻尼方法的化极结果能更有效地突出地质信息.     

抗差自适应Kalman滤波在动态变形监测中的应用研究

抗差自适应Kalman滤波算法中,抗差等价权矩阵和自适应因子的计算,要求观测信息具有多余观测量且准确可靠,但在动态变形监测应用中,通常滤波观测值仅为三维坐标且存在较强噪声和粗差的影响。为此,先对该算法中的自适应因子和抗差等价权矩阵的计算进行研究和改进,然后计算了某高速公路边坡的GPS动态监测数据。结果表明,抗差自适应Kalman滤波能够有效地抵制动态变形监测中观测值异常的影响。     

一种稳定的波场延拓反Q滤波算法

反Q滤波算法中稳定性是关键。针对地下介质为层状Q模型结构,本文提出一种稳定的反Q滤波算法,基于波场延拓原理,对于每个常Q层反Q滤波算法分两步完成:(1)利用一种稳定的波场延拓算法,将地表波场记录直接延拓到当前层上;(2)在当前层内进行常Q反滤波。在步骤(1)中,利用稳定的波场延拓算法并结合稳定因子,将地表波场记录直接延拓到当前层上,避免了上覆层反Q滤波误差的累积。在步骤(2)中,对于当前层内的反Q滤波,通过对该层地震信号的Gabor谱分析,拾取时变的限幅频率,并确定相应于该频率值的限幅增益值来实现反Q滤波算法的稳定性。最后通过理论模型与实际数据验证了算法的优越性。     

变阻尼共轭梯度算法及其性能分析

为了提高反演的分辨率和计算效率,本文在传统阻尼共轭梯度法的基础上,提出了变阻尼共轭梯度算法.由于在最小二乘线性反演过程中,多数情况下都要计算偏导数矩阵,而偏导数矩阵列向量的长度大小决定了解向量在对应分量方向上前进的快慢,变阻尼共轭梯度算法的提出正是利用了偏导数矩阵的这一特点.从理论上讲,它要优于传统的固定阻尼共轭梯度法.最后通过计算实例证明了该算法计算精度高,稳定性好,收敛速度快.因此采用变阻尼共轭梯度算法进行地球物理反演是切实可行的.     

逆子波域消除多次波方法研究

SRMA（与表面相关多次波的衰减）算法包含预测和相减两步.相减算法中,当多次波与反射波同相轴相交时,如何有效减去多次波、保留反射波,是面临的主要问题.通过分析非正交性对滤波器（逆子波）的影响,可以证明：逆子波因非正交性产生的误差呈近似的高斯分布.在此基础上,本文提出了在逆子波域（单道自适应相减滤波的滤波算子的集合）,利用其误差的概率分布特征,对逆子波进行估计,用逆子波的估计对逆子波进行校正来消除多次波的方法.其步骤为：首先用SRMA方法预测出表面多次波,并对每一单炮进行单道自适应相减,得到逆子波,形成逆子波域;其次,在逆子波域采用中值滤波,提取接近真实逆子波的逆子波估计;第三,在逆子波域用逆子波估计对畸变的逆子波进行校正;最后采用校正后的逆子波来衰减多次波.通过简单模型和SMARRT模型的测试,该方法不仅能够有效减去多次波,而且在相交的区域,能够保持反射波同相轴的连续性并恢复其正确的振幅.     

用一种新滤波函数进行CT图像的局部重建

本文阐述了用一种新型滤波函数进行CT图象局部重建的原理,并给出了这种重建的二次误差估计.本文选用一种新的滤波函数,它的旁瓣衰减很快,经证明在可能选取的滤波函数中,它的衰减是最快的.新滤波函数由于空域的迅速衰减性,因此具有良好的局部重建性质.同时保持该滤波函数的其它特性,能够有效地抑制噪声,补偿混迭效应,并缓解Gibbs效应.应用新滤波函数,忽略远离计算点的数据,从而直接用CBP实现图像局部重建.在计算机上对模拟和实测数据进行局部重建的实验结果表明本文建议的图像局部重建算法简便快速,重建图象Gibbs效应减小并与用全部数据的CBP重建图像比较仍具有较高的空间及密度分辨率.     

改进的阻尼法低纬度磁异常化极方法

中低纬度地区的磁异常资料进行化极处理能减少或消除斜磁化影响,对提高磁测资料的利用程度有重要意义.本文在总结分析原有的低纬度磁异常化极过程的基础上,改进了原有化极方法的滤波过程.改进后的阻尼化极方法,采用将振幅谱和相位谱分开进行滤波的过程,过程中只需要考虑对化极因子的振幅谱进行压制即可.通过公式分析表明该方法更加合理,减少了人为因素的影响;经单个模型和组合模型试验表明,此化极方法准确、稳定、可靠,可以有效地突出低纬度地区的异常体信息.

     

基于浮点遗传算法的近震定位方法

文中介绍了一种适用于基于浮点遗传算法的地震定位问题方法, 在该算法中选用了浮点基因表示方法、 参数和基因取值范围一致来避免编码解码, 提出了算法的动态结束条件及处理方法。 通过实际运算检验证明使用该方法来实现地震定位可以避免地震定位中的不稳定性和对初值的依赖性, 从而保证地震定位的可靠性和稳定性, 在时效方面也达到了相当的水平, 是一个具有实际使用价值的算法。     

Noncyclic variation and quiet day daily variation

Summary Several methods for the determination of a linear trend and the periodic variation in geomagnetic hourly values are compared. Where the true model for the periodic variation is not known or difficult to determine the recommended method (where there exists a noncyclic variation which may be regarded as a linear trend) is that presented by Chapman and Miller and denoted here by method B. Where a detailed analysis of the data is desired the method presented here (and denoted as method C) would in general give more reliable estimates of the parameters, once the correct model is chosen, than methods A and B. If only 24 hourly values are to be used method A would give estimates of the parameters which are almost as good as those for method C.     

用广义射线和有限差分计算近场理论地震图的混合方法

本文提出一种结合广义射线和有限差分理论而形成的混合方法（以下简称混合方法）．在介质均匀区（包括平行介质分层区）应用广义射线理论，在介质不均匀区应用有限差分理论．在这种结合过程中，把均匀区广义射线的结果作为不均匀区有限差分区波场连续的输入．对于均匀介质模型，混合方法理论地震图和广义射线理论地震图、广义反射透射系数矩阵及离散波数法地震图的对比，取得了满意的结果，在数值上证明了该方法的正确性．计算了１９９０年８月５日Ｈａｋｏｎｅ地震前震在日本Ａｓｈｉｇａｒａ山谷盲测点的响应，理论地震图和实际资料的比较同样取得好的结果．这种方法可用于地震勘探、工程地震、强地面运动的研究．     

均匀倾斜层状地层中地震波球面扩散补偿方法研究

在分析地层深度、炮检距以及地层倾角变化对球面扩散补偿因子影响的基础上,总结了不同因素对补偿因子影响的趋势,以及资料处理中常规球面扩散补偿方法的应用条件．实际应用效果表明,该方法能更加合理的保证不同深度地震信息能量的均衡程度,提高资料的可信度。     

多导体柱微波成像的连续编码遗传算法

提出一种能同时重建多个二维良导体目标外形轮廓的新方法——连续编码遗传算法．将目标的横截面轮廓近似表达为三角级数形式,由边界条件得到一积分方程组,在此基础上将成像问题转化为约束优化问题,级数的各项系数为待优化量．积分方程组为约束条件,目标函数定义为实际测量的散射场与反演过程中得到的散射场之间的相对误差函数．采用连续编码遗传算法求解,优化过程通过选择、交叉、变异等遗传操作的选代而实现,待优化目标函数进行线性变换并采用模拟退火原理确定目标后表达为适应度函数,采用联赛选择与比例选择相结合的选择机制和单点交叉方式,变异操作通过对基因施加微小随机扰动实现．上一代中适应度最高的个体直接保留．数值模拟反演实验验证了方法的有效性．与其他反演方法相比,本法具有简单、通用、鲁棒性强等特点．     

多导体柱微波成像的连续编码遗传算法

提出一种能同时重建多个二维良导体目标外形轮廓的新方法——连续编码遗传算法．将目标的横截面轮廓近似表达为三角级数形式,由边界条件得到一积分方程组,在此基础上将成像问题转化为约束优化问题,级数的各项系数为待优化量．积分方程组为约束条件,目标函数定义为实际测量的散射场与反演过程中得到的散射场之间的相对误差函数．采用连续编码遗传算法求解,优化过程通过选择、交叉、变异等遗传操作的选代而实现,待优化目标函数进行线性变换并采用模拟退火原理确定目标后表达为适应度函数,采用联赛选择与比例选择相结合的选择机制和单点交叉方式,变异操作通过对基因施加微小随机扰动实现．上一代中适应度最高的个体直接保留．数值模拟反演实验验证了方法的有效性．与其他反演方法相比,本法具有简单、通用、鲁棒性强等特点．     

井地有限线源三维电阻率反演研究

有限线源的三维电阻率反演,是地学探测研究的重要课题之一.本文主要讨论了以下几个方面的内容:(1)结合反演方程式改进的共轭梯度算法;(2)基于扰动法的线源Jacobi矩阵的近似计算方法及光滑系数矩阵在各个方向上的光滑因子的计算公式;(3)通过分析阻尼系数与修改量校正量对反演结果的影响,提出了利用校正系数对模型修正量进行校正;(4)较为系统的讨论了阻尼系数λ对反演分辨能力的影响,较好的改善了三维电阻率反演中的电性异常体重心的"上漂"现象.数值实验表明,改进后的共轭梯度法反演对初始模型的依赖程度较小,反演能够稳定收敛,对模型的空间位置分辨率较高.     

起伏地表QR径向基函数有限差分及其在弹性波逆时偏移中的应用

弹性波逆时偏移不受倾角和偏移孔径的限制, 能够实现任意复杂构造的高精度多波成像, 是目前最精确的多分量资料偏移成像方法之一.逆时偏移算法的核心是波场延拓, 传统波场延拓以水平基准面为边界条件, 基于固定采样步长进行规则网格剖分, 采用阶梯近似法处理起伏地表和复杂构造界面时会产生台阶散射, 严重影响起伏地表复杂构造的成像精度.基于无网格节点模型, 定量分析了弹性波模拟中径向基函数有限差分法的频散关系和稳定性条件.基于此, 提出一种基于QR径向基函数的高精度有限差分方法, 并提出一种优化的起伏地表自适应节点剖分方法, 推导了精确的无网格自由边界条件和弹性波无网格混合吸收边界条件, 形成了新的基于无网格的起伏地表弹性波数值模拟方法.此外, 本文将此无网格径向基函数有限差分方法应用于精确的纵横波场矢量分解公式, 实现了起伏地表弹性波逆时偏移成像.通过对高斯山丘模型, 起伏凹陷模型和起伏地表Marmousi-2模型进行数值试算, 验证了本文方法的有效性和可行性.

     

地震波动强度变化的数学建模分析与仿真

传统利用灰色关联分析方法对地震波动强度变化进行数学建模分析与仿真时,对地震波动强度变化的数列进行仿真分析时,忽略了地震波动强度的时间属性对结果的影响,导致分析结果准确性较低。本论述提出新的地震波动强度变化数学建模分析与仿真方法,通过地震波动强度序列的经验分布确定门限自回归模型的门限值,依据该门限值、AIC最小准则以及最小残差平方等方法获取地震波动强度序列的门限自回归模型,分析自回归模型的极限环和振荡的属性特点,得到地震波动强度变化的初步数值模拟结果。本论述构建了基于均生函数的地震波动强度序列的数学模型,通过均生函数数学建模方法拟合地震波动强度时间序列,依据时间序列基于双评分准则选取拟合周期,实现地震波动强度的数值仿真。实验结果表明,所提方法对地震波动强度变化模型具有较高的准确性和稳定性。     

井间电磁法综述

地面电磁法在实际工作中在很大程度上受到了探测深度和分辨率的限制.为了克服这些缺点,出现了井间电磁法.井间电磁法指的是在两个(或多个)钻孔中分别发射或接收电磁波信号,利用电磁波信号进行成像并探测井间物理性质的地球物理方法.由于发射机和接收机可以分别放置在很深钻孔中,其具有大透距、大探测深度的特点,因而广泛应用于工程环境物探、矿产勘查、石油勘探等中.针对不同的应用,产生了各具特点的一些特殊方法,包括井间无线电波成像、跨孔雷达、井间电磁成像.井间无线电波成像仪目前只测量电场强度数据,工作频率低,一般是单频的电磁波,频率范围通常在1 kHz至10 MHz.由于缺少走时数据修正射线路径,井间无线电波成像主要是进行基于直射线追踪的衰减层析成像.井间无线电波成像既可用于工程与环境地球物理也可用于找矿.跨孔雷达是钻孔雷达的一种探测方式,用高频电磁脉冲探测两个井孔间介电常数和电导率的变化.跨孔雷达层析成像也叫地质雷达CT,既可进行走时成像,还可进行衰减成像.一般来说,地质雷达CT的电磁波工作频率较高,中心频率通常在10 MHz和1 GHz之间,因此在分辨率指标上占有优势,跨孔雷达主要用于工程与环境地球物...