基于蒙特卡洛方法的AVO属性定量分析研究

AVO属性分析可对碎屑岩沉积环境下的岩性流体性质进行识别,但由于存在一个大范围的岩性导制类似的AVO效应,预测只能定性分析。这里利用测井曲线进行趋势分析,可以得出不同岩性参数的概率密度函数,基于蒙特卡洛方法生成概率统计模型,按照不同深,能够得到砂岩不同流体的概率分布图。并且,对叠前地震数据提取目的层计算得到泊松比属性图,选取不同区域对模型得到的概率分布图进行校准;用贝叶斯公式计算不同流体的概率值,可以达到定量分析的目的。     

自适应基质矿物体积模量提取和流体替换验证(英文)

岩石物理学研究中Gassmann方程被广泛应用于预测岩石中的地震波速度,由于输入的基质矿物体积模量参数不准确,极大的影响预测结果的可靠性,特别是复杂基质矿物组合的碳酸盐岩储层。因此本文结合Russell流体因子和Gassmann-Boit-Geertsma方程计算式,通过引入干岩石骨架泊松比,提出了一种基质矿物体积模量提取方法,能够自适应反演岩石基质矿物等效体积模量,提高流体影响预测的可靠性,通过实际资料流体替换验证,该方法的预测结果是可靠的,并且计算效率高、适应性强。     

弹性波场和稳恒电流场的Frechet导数计算方法

在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的编数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度ΔＧ求得。把任意场源正演问题的解和Ｆｒｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解     

一种基于逆时偏移的高效照明补偿成像方法

对地震波偏移成像而言,由于观测采集系统不规则、不完备采样以及地下复杂构造对波场传播的影响,在成像点处通常存在照明不均匀的现象,对偏移成像结果进行照明补偿是提高地质体成像精度的有效方法.本文给出了一种基于逆时偏移实现的高效照明补偿成像方法.在高频近似下,本文首先定义了衡量地下照明均匀程度的雅克比矩阵,并以此为基础给出照明补偿算子.通过波场的边界积分表达式,进一步将照明补偿算子利用外推波场进行表达,避免了照明补偿算子中高频格林函数的显式计算,以此构造出高效的照明补偿成像算法.合成数据和实际数据的数值实验结果验证了本文提出方法的正确性.     

弹性波场和稳恒电流场的Frchet导数计算方法

在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的偏导数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度△Ｇ求得。把任意场源正演问题的解和Ｆｅｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解。文中以稳恒电流场问题为例，给出了直流电法视电阻率偏导数的数值结果。     

伪谱法运算效率的提高及井问波场的正演模拟

伪谱法求解波动方程的常规算法,运算量大,运算时间长。笔者综合提高伪谱法计算精度和效率的方法,改进了常规计算流程里计算傅氏变换和反变换中需要多次使用计算的常量,把这些常量提到时间循环体外,从而避免很多重复运算,减少了运算次数,使运算效率提高了一倍以上,并结合目前隐蔽圈闭在油气勘探中的重要性,构建了隐蔽圈闭中一种特殊的构造,即透镜体模型,模拟了高速和低速透镜体模型中的井间波场特征,对比分析了2种模型的合成井间地震记录、时间切片及振幅特征,与理论分析的结果能够很好的相吻合。     

伪谱法运算效率的提高及井间波场的正演模拟

伪谱法求解波动方程的常规算法,运算量大,运算时间长。笔者综合提高伪谱法计算精度和效率的方法,改进了常规计算流程里计算傅氏变换和反变换中需要多次使用计算的常量,把这些常量提到时间循环体外,从而避免很多重复运算,减少了运算次数,使运算效率提高了一倍以上,并结合目前隐蔽圈闭在油气勘探中的重要性,构建了隐蔽圈闭中一种特殊的构造,即透镜体模型,模拟了高速和低速透镜体模型中的井间波场特征,对比分析了2种模型的合成井间地震记录、时间切片及振幅特征,与理论分析的结果能够很好的相吻合。     

基于射线参数估计的快速 2D Kirchhoff叠前深度偏移

叠前深度偏移是油气勘探领域解决复杂地震构造成像问题的主要手段.这里介绍一种基于射线参数估计的快速2D Kirchhoff叠前深度偏移方法.该方法根据炮点和接收点的相对位置,对入射波的射线参数进行先验约束,辅以设定振幅阈值,仅对主要反射事件进行偏移等,显著减小了射线追踪计算,较传统Kirchhoff叠前深度偏移显著提高了偏移成像效率.数值模拟计算结果表明,该算法能正确对复杂构造进行快速成像.这里介绍的快速Kirchhoff叠前深度偏移方法,可用于野外现场质量监控与精确偏移成像前的速度分析与建模.     

混合法计算三维非均匀介质中的电磁场

介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。     

基于分数阶傅里叶变换的混合Cadzow滤波法

常规频率域SVD滤波法在随机噪声压制处理时,容易造成弯曲同相轴损伤。针对这一问题,本文提出一种基于分数阶傅里叶变换的混合Cadzow滤波法压制三维地震数据中随机噪声的方法。首先应用分数阶傅里叶变换,将地震数据变换到时频面,再依据Eigenimage滤波法与Cadzow滤波法建立混合的高维Hankel矩阵,然后对其运行奇异值分解,最后通过秩约化的方法来压制随机噪声。此方法用于四川某地区的地震数据处理,结果表明该方法可有效地去除随机噪声,保护有效信号,提高叠后地震数据的信噪比。