二维地层结构和速度的走时反演

介绍了在倾斜均匀各向同性层状介质情况下，根据地震反射走时反演二维地层结构和地层一地震射线层析成像方法－走时反演。正演模型射线追踪是根据Ｆｅｒｍａｔ原理实现的，即求解满足该原理的非线性方程组得到射线与界面的交点，进而计算相应的走时。反演是先假设一初始模型，用最优化方法使射线追踪走时的残差极小。最后计算了分辨矩阵和信息密度矩阵，以评价反演结果。对有噪情形也进行了反演。     

用遗传算法反演一维声波介质

局部线性化的反演方法不仅在初始模远离真值时收敛速度幔，而且往往易陷入局部极大值中，而遗传算法（ＧＡ）便能解决这一问题，它是一种全局搜索法。强然它同模拟退火法一样利用转移概率进行搜索，但其收敛速度一般比模拟退火法快。本文引进了多个目标函数进行综合评价和实灾变过程，用单道地震道对一维声波介质的速度和密度同时进行了反演，反演了加噪地震道和含薄层介质模型，取得了满意的效果。     

非迭代,非线性二维重力反演计算方法

本文叙述了非迭代、非线性二维重力反演方法的原理、模型计算及应用效果等。与以往诸多的求界面重力反演方法相比,该法的主要特点是：将界面的深度用一幂级数的形式表示,且幂级数系数的计算与界面上下地层的密度差无关。因此,一旦求出了幂级数的系数,由界面的密度差可直接得到界面的深度,且幂级数收敛的速度很快。     

张量CSAMT数据的二维反演

在一单行本中,我们提出了一种ＣＳＡＭＴ正演模拟算法并说明了这种算法在解释野外数据中的应用效果,在本文中,提出了一种用于二维电阻率模型的ＣＳＡＭＴ数据反演算法,这种算法包含一个三维源并可以反演（ＴＥ）模式及（ＴＭ）模式下的数据,用这种算法对假设的数据进行试算,结果表明,不用做平面波校正就能对远区场,过渡区场及近区场进行反演,最终的模型结果显示,原始电阻率模型可被恢复,而不存在任何由上源于效应而产生的     

二维地层结构和速度的走时反演

介绍了在倾斜均匀各向同性层状介质情况下,根据地震反射走时反演二维地层结构和地层速度的地震射线层析成像方法——走时反演。正演模型射线追踪是根据Fermat原理实现的,即求解满足该原理的非线性方程组得到射线与界面的交点,进而计算相应的走时。反演是先假设一初始模型,用最优化方法使射线追踪走时与观察走时的残差极小。最后计算了分辨矩阵和信息密度矩阵,以评价反演结果。对有噪情形也进行了反演。     

基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演

全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。     

生物遗传算法在速度参数反演中的应用

遗传算法是一种新发展的优化算法，具有全局搜索、适应解各类问题的特点。本文应用遗传算法求取最佳叠加速度和反演层速度，不仅可保证精度，还比常规方法减少一倍多的计算时间，获得了较高的计算效率。     

波动方程速度与震源函数联合反演

本文使用遗传算法对均匀弹性介质进行波速和震源函数的联合反演。通过数值计算结果表明,波速和震源函数的理论值与反演值比较吻合。方法简便,有效,并具有一定的抗噪能力。     

一维反散射层速度反演

利用参考速度进行波动方程的坐标变换和玻恩近似,建立了把散射势与速度扰动项相联系的积分方程,通过求解振动项,获得了层速度与界面的反演结果。文中所用的合成记录反演结果表明了本方法运算速度快、精度高,前提条件少。     

用遗传算法反演一维声波介质

局部线性化的反演方法不仅在初始模型远离真值时收敛速度慢,而且往往易陷入局部极大值中,而遗传算法(GA)便能解决这一问题,它是一种全局搜索法。虽然它同模拟退火法一样利用转移概率进行搜索,但其收敛速度一般比模拟退火法快。本文引进了多个目标函数进行综合评价和灾变过程,用单道地震道对一维声波介质的速度和密度同时进行了反演,并反演了加噪地震道和含薄层介质模型,取得了满意的效果。     

反演瞬态瑞雷波频散曲线映射二维横波速度剖面

由于横波速度与岩土力学参数密切相关,因而是工程地质勘察中的关键参数。基于前人的思想,采用等厚薄层权重自适应迭代阻尼最小二乘法自动反演横波速度,提高横波速度的计算精度,进而映射二维横波速度剖面。用上述方法对某高速公路振冲压碾效果进行了无损检测试验,其结果与钻探和开挖验证一致。     

用遗传算法进行瑞利波反演

瑞利波反演属于非线性最优化问题,已有的局部线性化方法(如阻尼最小二乘法)常使解估计陷入目标函数的局部极小值之中,且严重地依赖初始模型.另外,目前被广泛应用的一些简单的瑞利波反演方法也存在明显的缺陷.全局优化法-遗传算法大大放松了对初始模型选择的要求,且不易陷入局部最优解中.本文用已有的简单宜行的反演解释方法获取初始模型,从而确定模型参数的搜索范围,再用遗传算法反演得到最终的介质模型,效果非常理想.     

复杂地质构造波动方程数值模拟

利用全声波方程进行地震波场的数值模拟，虽然可以计算全声波波场值，但是也常常会出现层间多次波，给资料解释造成一定的困难。作者在本文中根据无反射声波方程，采用有限差分数值解法，实现了自激自收时间剖面的制作，避免了层间多次波的产生。本方法不但允许速度有纵向变化，而且允许速度有横向变化，这使得模拟复杂地质构造的地震响应成为可能。实际计算结果表明，该方法是一种行之有效的自激自收时间剖面的计算机数值模拟方法。     

小波变换在波阻抗反演中的应用

根据国内外广泛采用的波阻抗随机反演理论,提出了多道小波波阻抗反演方法,消除了噪音对反演结果的影响。针对地震剖面上不同频带和时段上有不同的信噪比,以及相邻地震道反射波有效成分在波形和能量上有较强的相关性特点,运用多道记录同时进行波阻抗反演,在迭代过程中将实际记录与模型合成记录的残差道进行小波分解,形成残差道分频小波剖面,然后利用Ｋ－Ｌ变换提取分频剖面中的相干部分用于波阻抗变化量的计算,从而避免了噪音参与波阻抗反演。展示了理论模型与实际应用的例子。     

南极和北极地区在全球变化中的作用研究

多次波问题是海洋地震勘探中最突出的问题之一 ,如何有效地压制多次波是解决问题的关键。多次波数据处理的压制方法和技术分为运动学的滤波方法和动力学的波动方程方法两大类。波动方程方法由于几乎不需要先验信息和具有较好的压制效果 ,已经成为多次波压制方法和技术发展的主要趋势。文章首先对波动方程方法进行了较全面的综述 ;然后 ,对波动方程方法中的反馈迭代反演方法做了深入的讨论 ;最后 ,通过理论模型的正反演计算 ,验证了运用波动方程反馈迭代反演法不需要知道震源函数、地下岩层结构和岩性等先验信息。此方法适用于压制多种类型的多次波 ,具有优异的振幅保真性等独特的优点。     

瞬态瑞雷波反演横波的SVD算法及其应用

利用奇异值分解算法对瑞雷波频散曲线进行反演,提高横波速度的反演精度,引入权重矩阵提高数据分辨率,采用自适应修改阻尼因子提高迭代效果并协调分辨率与解的关系.某高速公路软土路基实测资料试算结果表明,用该算法对瞬态瑞雷波频散曲线进行反演,不但具有稳定性好、精度高、分辨能力强的特点,而且能自动分层和反演地层参数,在高速公路软土路基探测等领域中具有广阔的应用前景.     

二维视电阻率断面的快速最小二乘反演

快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方