高阶交错网格有限差分弹性波场模拟的精度分析

交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明：时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。     

沁水盆地南部中高阶煤高压甲烷吸附行为

煤层气吸附作用是发生在煤基质内表面的物理过程,而煤岩复杂孔裂隙网络为高压甲烷吸附提供了丰富的空间。开展沁水盆地南部高阶煤30℃高压甲烷等温吸附实验,结合煤岩煤质参数与孔隙特征参数,通过改进的D-R模型分析了煤岩性质、孔隙特征与吸附参数的相关性。煤岩性质对最大吸附能力和吸附热参数的影响是多因素叠加的综合效应,而最大吸附能力与微孔体积,吸附体积校正参数与大中孔比表面积呈较好的正相关性,表明甲烷分子在煤基质内表面会根据孔径尺度大小呈现不同的吸附方式。据此提出高压甲烷在煤基质微孔中呈紧密堆积状态而在大中孔中呈多层分子堆叠状态的新认识,为进一步研究煤层气吸附机理提供了新的思路。     

基于f-x域时频非凸正则化低秩矩阵近似的共偏移距道集去噪方法

随机噪声压制是地震数据处理的关键环节,而时频稀疏低秩近似算法逐道处理地震数据过程中无法利用信号的道间相干性。为此,将时频稀疏低秩近似与f-x域去噪结合,提出一种f-x域时频非凸正则化低秩矩阵近似算法。该算法对f-x域中每一单频分量作时频分解后,再对时频系数矩阵作低秩矩阵近似计算,能够利用信号和噪声的时频谱差异实现非平稳信号去噪处理。与共炮点道集和共中心点道集相比,共偏移距道集具有平缓甚至接近水平的同相轴结构,基本满足f-x域去噪的线性同相轴假设前提,建议将所提算法应用于共偏移距道集去噪处理。通过数值模拟和实际地震数据试算,证明本文方法能够有效压制随机噪声,同时保持有效信号不被损害。     

基于新的差分结构的时-空域高阶有限差分波动方程数值模拟方法

传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性.     

一种模拟非均匀介质中弹性波传播新的k-space方法

传统的伪谱（PS）方法,采用傅里叶变换（FT）计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分（FD）近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.

     

一个近似的VTI介质声波方程

声学近似是提高各向异性介质准纵波数值模拟和偏移处理计算速度的有效方法之一.为了获得VTI(具有垂直对称轴的横向各向同性)介质声波方程,根据VTI介质的精确相速度表达式,用待定系数法简化其中的根号项,得到一个近似相速度公式,由此公式通过反Fourier变换推导出一个时间二阶、空间四阶偏微分方程.定量相速度计算表明,当Thomsen参数ε=δ(即相应介质为椭圆各向异性介质)时,由该方程所确定的相速度是精确的;当ε≠δ时,该方程所确定的相速度随ε和δ差值的增大逐渐增加.该方程所确定的最大相速度百分比相对误差当ε=0.1、δ=0.2时为0.13%,当ε=0.8、δ=0.3时为-1.65%.有限差分数值模拟算例表明该方程是一个纯准纵波方程,其数值模拟波场快照中没有准横波.     

双线性插值近似网格的栅格数据投影变换

在栅格数据投影变换常用算法的基础上,介绍一种新的基于双线性内插近似网格算法的栅格数据投影变换,并详细描述其算法过程,认真分析其特点,并证明该算法在大范围栅格数据投影变换中有显著优点。     

利用极化角点指数进行数字曲线快速多边形近似

针对曲线结构特征,引入极化角点指数作为特征角点的候选依据,提出了一种与曲线起始点位置无关的基于角点检测的多边形近似算法。实验证明,此方法快速有效,结果稳定,噪声抑制效果较好。     

基于改进的Farid高阶统计量模型隐写检测

在分析经典Farid高阶统计量模型基础上,采用离散Haar小波对其线性预测模型相关性区域进行定位描述,并改进其偏向不一致缺陷,通过数学推导提出R-Farid线性预测模型;并提出了R-Farid高阶统计量模型。实验结果显示,其具有更好的隐写检测效果。