频率域可变网格数值算法研究

来源 :中国地球物理学会第二十八届年会 | 被引量 : 0次 | 上传用户:jdbc_jiang
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  本文通过在频率域使用可变网格数值算法计算射线参数和单频波场,发现对于单频波场的计算,网格越细越好。频率越高,使用的网格应越细;对于射线参数的计算,网格并不是越细越好。频率越低,使用的网格应越粗。此外,要想有效地压制数值频散,网格的选取不仅要满足频散关系,还要使网格间距的变化尽可能光滑连续。
其他文献
本文研究利用有效孔隙度在泥质含量高的地层中进行流体替换的方法。通过该方法的原理分析结果发现,利用有效孔隙度进行流体替换的改进方法解决了常规方法中存在的泥岩束缚水问题,使得在低孔隙的泥质砂岩地区得到的弹性参数差异更明显。本文指出在很多常规方法无法解释的观测地区,可以利用这种新方法得到的结果与测井数据以及三维地震数据,特别是四维地震数据进行匹配、解释。
本文在收集大洋洲(澳大利亚)地区卫星重力与磁力数据的基础上,对其进行常规处理,处理手段包括异常分离、磁力化极、延拓、水平垂直方向导数等方法,得到了该地区区域布格重力异常不同高度的向上延拓、不同方向导数异常及磁异常化极等常规处理结果。研究并编程实现了位场归一化差分方法,并利用该方法对某模型数据及澳大利亚大陆重力数据进行处理。
本文介绍了Q值提取式的推导过程,通过设计多个地层Q值模型,分析了谱比法、质心频率法和对数谱统计属性法三种方法提取Q值,在无噪情况下的精度,在含噪情况下的稳定性和可靠性。
本文出了一种走时层析与时域全波形联合反演的方法,利用走时层析快速、高效地获取质量较好的初始速度场,然后基于时域全波形反演实现高精度的刻画细节的叠前速度场。模型试算表明,由于采用较好的初始速度场作为输入,使得该联合反演具有较高的反演效率和反演精度。
本文首先通过不同参数的2D和3D理论模型试验总结了水平总梯度、解析信号振幅、Tilt导数、Tilt导数水平梯度、Theta图法、归一化标准差等目前常用的位场线性特征增强方法的应用条件和优缺点,得出提取效果最好且受场源深度影响较小的方法。另外,通过将不同的延拓高度作为不同尺度对边缘检测结果进行多尺度叠加,这样根据其变化规律不仅能够反映线性构造的走向、延伸及连续性还能显示出不同深度层次的变化情况和断裂
本文介绍敏感核的准确计算是实现波形层析成像的关键。本文在二维近似解析中心差分方法的基础上,提出了三维近似解析中心差分方法来求解波动方程,进而计算敏感核。
本文在频率域用全波形反演方法基于Stokes方程做双参数的同时反演研究,反演了Stokes方程中速度、粘滞系数两个参数,反演中密度设为已知。结果表明了算法对初始模型的依赖性不强,以及频率域全波形方法在Stokes方程双参数反演中的可行性。
本文通过对2D拟合数据叠加剖面数据和实际数据的计算,结果证明FGST比短时傅里叶变换计算精度高,由于它基于傅里叶变换产生的谱优于小波变换,比普通的S变换速度更快。而FGST所求得的相位剩余与不同的相位值转换有关。可以对不容易察觉的断裂做出敏感的反应,在对地层不连续性成图上有着重要的意义。
本文认为上下左右行波分解成像条件可以有效的解决假象问题,但需要额外的计算时间。这些因素导致逆时偏移比较昂贵。基于此,本文提出了仅在目标成像区域内进行源波场重建与上下左右行波分解的逆时偏移方案,来达到降低逆时偏移成本的目的。
本文研究了变步长高阶有限差分叠前逆时偏移成像方法,在保证偏移成像质量的前提下,可以有效地提高计算效率,减少存储空间。通过模型试算及实际应用表明,变步长高阶有限差分逆时偏移有效地解决了近地表低速资料的逆时偏移成像精度问题,与全区域都采用小空间步长和时间步长的常规算法相比,既提高了计算效率,又保证了波场计算的精度和偏移成像的效果,是一种简便可行的方法。采用双倍变步长算法计算效率就可以提高2~3倍。另外