介质密度及Lame系数同时识别的研究

将非均匀介质密度及Ｌａｍｅ系数的识别作为弹性波方程系数反问题，求残差最小，并应用弹性波正演的有限元解及摄动法，使反问题求解转化为二次规划问题，应用Ｌｅｍｋｅ算法求解二次规划问题，并以解为新的初始值重复上述过程，可迭代确定介质的密度和Ｌａｍｅ系数，该方法计算效率较高；可引入介质构造的其他勘探信息，进行联合反演。     

半空间地表反射式弹性波多参数反演

在半空间中及Ｂｏｒｎ近似下，探讨了采用地表发射和地表接收配置下的弹性波方程，得到散射波总场并从中分离出散射波的四类波型转换（Ｐ→Ｐ，Ｐ→Ｓ，Ｓ→Ｓ，Ｓ→Ｐ）．在一个工作频率下，利用这些波型转换得到了多参数（ρ，λ，μ）的具有滤波－反传播形式的反演公式，并对反演公式进行了理论模型验算及分析．     

多锥结合壳流固耦合稳态动力响应分析

针对二维弹性波的反演问题，给出了一个同时重建三个弹性常数的层析成像方法，通过选择不同入射方向和偏振态的入射平面波，实现了三个参数的解耦，得到了相应的加权滤波反投影公式，并对一简单算例进行了数值反演模拟．     

基于多频率分量的叠前AVO多尺度反演方法

为提高频率域反演的分辨率和抗噪能力,减小频率域反演对初始模型的依赖程度,发展了多分量频率域叠前AVO弹性参数反演方法。考虑地震信号在频率域中不同频率成分可实现自动解耦,使得反问题可通过多分量迭代的途径求取。首先,利用贝叶斯公式将多频率分量联系起来,并结合稀疏约束先验提高反演算法抗噪能力。此外,采用超低频平滑模型提高多频率分量AVO反演算法的稳定性和可靠性。然后,模型测试综合验证了该方法的高分辨率特征和横向连续性,获得与理论模型吻合度较高的高分辨率反演结果。最后,实际资料表明该方法能够稳定的从叠前地震数据中获取与测井数据相吻合的高分辨率弹性参数信息,并可有效的运用到储层含油气性质检测中。     

颗粒增强复合材料的界面损伤对有效波速的影响

研究了具有非完好黏接界面的球形粒子对弹性波的散射以及含分布球形增强粒子的复合材料中的弹性波多重散射．首先,在弹簧模型下研究了单个球粒子的弹性波散射问题．然后在统计平均意义下研究了具有随机分布粒子的复合材料中的有效波以及有效波波速和衰减的计算公式．进一步,根据有效波波速和衰减对界面常数的依赖关系,研究了界面损伤对有效波波速和衰减的影响．在此基础上提出了一种依据超声测量数据对界面损伤进行无损估计的方法．针对Sic—Al复合材料进行了数值模拟．计算了不同界面损伤程度情况下的有效波波速和衰减．并利用遗传算法对含不同误差水平的合成实验数据就界面损伤进行了反演计算．数值结果表明对含随机误差的合成实验数据进行反演计算能对界面损伤程度作出近似的估计并且反演算法对不同误差水平具有相对稳定性．     

GRACE时变重力信号反演方法研究的进展和展望

为了解决GRACE时变重力研究全球或区域质量迁移中信号失真问题,学者们多采用信号反演方法,以恢复地表真实质量迁移信息．时变重力信号反演是卫星时变重力研究中最核心的关键技术之一．然而,目前反演方法众多,如单一尺度因子法、变尺度因子法、Forward Modeling、频率域反演方法、空间域最小二乘反演法等,它们之间的共性、差异性及适用性尚不十分清晰,是有待解决的重要科学与技术问题．为此总结了常用的GRACE重力信号反演方法研究的进展与现状,讨论了目前存在的问题,对比了不同反演方法的优缺点、应用领域差异,提供了反演方法选择的思路,并展望了未来重力卫星发展趋势,为卫星时变重力应用研究提供科学参考和依据．     

固体中弹性波衍射层析成像

针对二维弹性波的反演问题，给出了一个同时重建三个弹性常数的层析成像方法，通过选择不同入射方向和偏振态的入射平面波，实现了三个参数的解耦，得到了相应的加权滤波反投影公式，并对一简单算例进行了数值反演模拟。     

利用远场模式反演声波阻尼系数与区域

目的 研究时间调和声波的反散射问题,利用散射波的远场模式反演边界条件中的阻尼系数和边界.方法 应用汉克儿函数的组合来逼近散射波,最终把反散射问题转化为一个最优化问题.结果 给出了一种可行和精确的反演方法 ,并给出了数值例子.结论 给出的反演方法 不仅数学理论简单,而且反演效果也比较好.     

一维弹性波定层逐点反演方法及其精度分析

针对一类一维弹性波反演问题，导出了一种稳定的定层逐点反演方法，推证了该方法的参数解具有整体二阶精度，数值试验表明，该方法是有效而可行的.图1，表1，参7.     

生物软组织黏弹性参数反演方法研究

针对黏弹性参数反演中的模型选择问题,提出应用黏弹性材料的应力松弛曲线作为待反演的参数,不涉及任何黏弹性本构模型,采用有限元仿真与人工神经网络相结合的方法:通过有限元仿真得到不同应力松弛曲线下的剪切波传播速度,然后将仿真得到的数据作为训练样本输入到人工神经网络中进行训练,建立起生物组织黏弹性参数和剪切波传播速度的神经网络输出模型,有效解决了模型选择问题。最后通过粒子群优化算法反演得到组织的黏弹性参数。     

整体波场延拓法基准面校正

提出了一种新的近地表问题解决方法——整体波场延拓法基准面校正(WEDD)。借助地震理论中的互换原理和观测面沉降概念,在基准面校正过程中将炮点和检波点的处理统一起来,并利用炮点-全偏移距域双平方根算子实现波场的单域整体延拓,最后提取目标深度的波场值完成参考面的重新定义。模型试算结果表明了WEDD方法的有效性。     

群桩垂向粘弹性动力分析的拉普拉斯数值反演法

应用担普拉斯数值反演方法分析了群桩垂向粘弹性动力反应。利用弹性和粘弹性间比拟关系及Novak平面应变假设,在拉普拉斯空间获得了单桩和两桩体系垂向粘弹性振动的正确解。根据两桩间两两相互作用关系提出了群桩垂向粘弹性动力分析的方法。最后应用拉普斯数值反演技术就可获得时间域的解。本文的计算结果对群桩分析中所作的假设进行了验证并对弹性解和粘弹性解进行了比较。     

群桩弯曲粘弹性动力分析的LAPLACE数值反演法

应用拉普拉斯数值反演法分析了群桩弯曲粘弹性动力反应。在拉普拉斯空间获得单桩和两桩体系弯曲粘弹性振动的正确解。根据桩间两两相互作用关系给出了群桩弯曲粘弹性动力分析的桩头柔度矩阵。最后应用拉普拉斯数值反演技术得到时间域的解。本文计算结果对群桩弯曲振动分析中所作的假设进行了验证并对群桩的弹性解和粘弹性解作了比较。     

保幅炮域高斯波束偏移

推导基于高斯波束表示的波场反向延拓公式,并结合反褶积成像条件得到保幅的二维炮域高斯波束偏移公式。在高斯波束偏移方法的实现过程中,须通过最速下降法将偏移公式中的二维复值积分简化为一维以提高计算效率。给出炮域高斯波束偏移的具体计算流程。最后,分别对水平层状模型和Marmousi模型进行试算。模型试算结果验证了所提偏移方法的正确性和有效性。     

基于褶积理论的深度域地震资料反演

近年来随着深度域地震资料的普及,为更好地利用这些地震资料来解决地下地质问题,提出了基于褶积理论的深度域地震资料反演。基于褶积理论的地震资料反演,子波与低频模型对反演结果具有十分重要的影响。通过制作人工合成记录与数值模拟研究,得到了深度域子波与时间域子波的关系,为基于褶积理论的深度域反演奠定了基础。将深度域反演应用于南海西部某油田,通过与该区的时间域地震资料反演对比,发现两种反演方法子波与低频模型都有差异。当时间域地震资料的采样间隔约为深度域地震资料采样间隔与层速度之比的两倍时,深度域地震资料的反演结果将优于时间域地震资料反演结果。     

叠前弹性参数反演方法的关键技术

为实现对岩性和隐蔽油气藏的有效勘探,提出一种叠前弹性参数反演方法。该方法以Zoeppritz方程近似式为反演基础。通过比较不同介质条件下的Zoeppritz方程近似式的计算精度,确定适合的方程近似式。反演时,采用无拉伸动校正技术处理模拟的正演地震数据,使用射线追踪法抽取方程计算所需的角度道集。反演测试结果表明,相对于传统的叠后反演方法,该方法可以反演各种弹性参数,且反演结果的横向连续性强,可用于储层预测。     

基于动校正前CRP道集的多波联合AVO反演新方法

叠前AVO反演可以得到地层的弹性信息,继而进行地层岩性和流体类型识别。与速度信息相比,弹性模量信息更有利于储层描述。因此,本文采用叠前AVO反演方法直接反演得到纵波模量反射系数、横波模量反射系数和密度反射系数,联合纵波(PP)资料和转换横波(PS)资料进行AVO反演提高反演结果的精度和稳定性。常规AVO反演方法所采用的数据都是动校正后的CRP道集,动校正处理引入的动校正拉伸效应会影响反演效果,降低反演精度,尤其是密度的反演结果会受到较大影响。本文采用动校正前的CRP道集进行联合AVO反演,在反演过程中引入贝叶斯理论,从而提高反演的稳定性和精度。理论模型试算结果表明,本文提出的反演方法具有更好的稳定性、更高的精度以及更强的抗噪能力。     

一种引入吸收衰减的压缩感知薄储层识别方法

常规基于压缩感知稀疏反演方法是基于频域平稳地震褶积模型进行的,而实际地下介质是粘弹性的,这使得该方法的反演反射系数振幅缺失、位置偏移。针对该问题,本文通过引入时频衰减因子,构建改进的感知矩阵,将稀疏反演从常规的频率域拓展至衰减频率域。反演结果在一定程度上可恢复缺失的振幅,增强弱信号的识别能力。由于噪声影响,以上处理结果中仍存在噪声干扰,故在衰减频率域稀疏反演的基础上,引入平滑的高斯函数,对反演目标函数进一步优化,以有效压制残存的噪声干扰;之后将反演结果与褶积宽频子波可生成高分辨率地震剖面,由此形成了一种引入吸收衰减的压缩感知薄储层识别方法。薄层理论模型及含有河道砂储层的实际地震资料的处理结果表明,本文方法较常规方法,有效地增强了薄层弱信号的振幅及横向连续性,可在保证信噪比的情况下,提高地震资料的分辨率及薄储层的识别能力。     

基于黏弹各向异性理论的地震波场数值模拟与分析

基于Carcione黏弹各向异性理论及具有任意倾斜对称轴的横向各向同性介质(TTI)的本构关系,利用二维三分量高阶交错网格有限差分和完全匹配层边界条件对黏弹TTI介质的地震波场信息进行了数值模拟。研究发现:1改变黏弹TTI介质的弹性参数,对地震波传播造成显著影响。其中,C_(13)主要控制qS1波中质点的偏振方向,进而影响了qS1波波形和反射同相轴的形态;而C_(66)则主要控制qS2波中质点的偏振方向,从而影响了qS2波波形和反射同相轴的形态;2利用层状黏弹TTI介质进行数值模拟,可以很好地观察到拟横波在遇到反射界面时发生的横波再分裂现象;同时由于介质黏弹性的存在,导致地震波出现显著的衰减。这将为进一步研究各向异性介质中复杂的波场信息提供帮助。     

基于 POD 有限元法的粘滞声波方程震源波场重构

叠前逆时偏移是地震勘探中一种流行的地下结构成像方法. 其成像条件需要同时刻的震源波场值与检波器波场值. 这在实际计算中就需要把正演模拟的所有时刻的震源波场数据全部存储下来，存储量需求大. 虽然震源波场重构技术可以降低对于波场数据的存储需求，但会引入额外的计算复杂度. 为解决这个问题，本文提出了POD有限元法，并将其应用于粘滞震源波场重构. 这里的本征正交分解方法( Proper Orthogonal Decomposition, POD)方法是一种降维方法，能够在降低数据量的同时提供足够的计算精度. 数值算例显示，该方法比传统的有限元方法更节省存储空间，能够加快重构速度.