【摘 要】
：
本文将分数阶导数应用到粘滞介质理论中的Constant-Q模型中得到一种以空间分数阶偏微分方程描述的模型，并对此模型的性质做了分析。本文分别将推广的伪谱法和双向G算法应用于求解空间分数阶导数，进而求解粘滞声波空间分数阶波动方程，并对多个模型进行正演与分析，得到一些有益的结论。
【机 构】
：
中国石油大学(北京)油气资源与探测国家重点实验室 北京102249;中国石油大学(北京)CNPC物探重点实验室 北京102249
论文部分内容阅读
本文将分数阶导数应用到粘滞介质理论中的Constant-Q模型中得到一种以空间分数阶偏微分方程描述的模型，并对此模型的性质做了分析。本文分别将推广的伪谱法和双向G算法应用于求解空间分数阶导数，进而求解粘滞声波空间分数阶波动方程，并对多个模型进行正演与分析，得到一些有益的结论。
其他文献
本文提出一种基于主成分分析的频率域FWI算法，利用对残差矩阵的主成分实现了数据降维和FWI的提速，而且避免了串扰噪声的产生。根据主成分分析和震源编码技术在不同频率尺度的互补性，提出了一种联合的高效FWI算法。
本文研究以S震相做为参考震相，进行SS及其前至震相的提取，并由动力射线追踪方法对所选用的地震及相应的震相进行射线追踪，得到了可以成像的区域，用汪荣江的QSEIS程序进行波形拟合，以得到与地震记录所反映结构相对应的LAB深度分布及相应的速度结构。
本文提出用WEPIF法估算叠前纵波和转换横波Q值.假设震源子波可以用具有4个待定参数的常相位子波来渐进逼近，借助于粘弹介质中的单程波传播理论，推导出了叠前道集中各层反射波的EPIF和Q值的关系.进而，为提高瞬时频率的估计精度和抗燥性能，可以在小波域中高精度的计算瞬时频率.WEPIF法估计叠前纵波和转换横波Q值，计算结果较稳定，精度相对较高.
本文在前人的研究基础上，系统的推导了相位编码方法的理论公式，并将该算法进行了模型试算。为了验证算法的正确性及有效性，本文将该算法对复杂程度不同的模型进行了偏移成像试处理，以测试该方法对复杂构造模型的偏移成像能力。
本文应用L1和L2两种范数准则构建目标函数，并分别采用CG(共轭梯度法)和L-BFGS(拟牛顿法的一种)两种反演算法来研究频率域粘滞声波波动方程反演问题。本文充分考虑了地下介质的衰减特性，以粘滞声波方程为基础同时反演出速度和Q值的信息。又将L1范数和L-BFGS算法相结合提高了反演精度。
核磁共振(NMR)微观响应特征对于复杂孔隙结构的研究具有重要意义。论文研究了核磁共振微观孔隙数值模拟的随机行走法基本原理及其实现步骤；针对理想孔隙模型，利用随机行走法模拟了磁化强度衰减，验证了方法的正确性；计算了不同孔隙尺寸及形状、不同流体的核磁共振微观响应，研究了扩散半径、行走者数目等参数对模拟结果的影响，为下一步开展基于像素的数字岩石数值模拟准备了条件。
本文利用中国地震局华北科学探测台站的远震数据成功的反演出斐济-汤加地区的410km间断面的起伏情况，所得数据前人利用近震所做的数据吻合较好，侧向的深度变化规律与矿物学相变模型所预测的规律有一定的一致性。410km间断面异常升降的区域反映了斐济-汤加俯冲带区域结构的复杂性，同时说明了俯冲带的现象并不能完全用矿物学和热动力学上的模型来完全解释。
本文通过对Watters研究小行星撞击对地幔对流影响中所用模型的改进，使其更具有真实性和可信度。本文将小行星撞击效果视为两种等效效果的叠加：一是热异常，二是陨石地形的回升。本文采用了分层常粘度，并同时考虑了固定和自由滑移两种底部边界情况，对上地幔对流进行了数值模拟，结果表明：直径1000~2500km的陨石坑会对上地幔对流产生较强的影响，包括一定程度地幔柱的偏移，但不会导致对流格局的重组，而是否会
本文提出了一种基于反射波相位的界面方位判定方法，指出该方法是基于单极激发、偶极接收的测井仪，同时介绍了该界面方位判定方法的具体判定过程。该方法目前在数值模拟中已经得到了验证，并提出利用反射横波的相位进行界面方位判定，从理论上来说也是可能的。
本文设计了两个数值模型（模型一、模型二）对煤层厚度与地震属性关系进行研究。通过对所设计的模型的分析对比，得到以下结论：(1)尽管煤层厚度较薄、沉积特征有差异(煤层合并分叉等)，但是其初至对应于煤层的顶界面，可以利用初至确定煤层的顶界面；（2）煤层厚度与属性值之间基本呈线性关系；（3）可以通过地震属性确定煤层内部构造突变位置及煤体类型变化部位。