基于匹配滤波的多次波压制方法研究

常规滤波方法对浅海多次波压制不能取得理想效果.通过改变模型表面边界条件,求解波动方程实现多情景地震数值模拟.利用自由边界条件模拟时获得含海平面强反射多次波地震记录,而利用吸收边界条件时得到理想条件下无多次波地震记录.对比分析多次波对有效地震信号的影响得到海平面强反射多次波改变了有效地震信号的振幅、频率和相位.利用归一化方法计算得到压制多次波滤波器算子,并作用于含海平面强反射多次波地震记录.叠加地震剖面对比分析表明,归一化滤波器算子较好的压制了地震记录中的多次波,提高地震资料分辨率.最后,根据多次波产生条件,给出了利用实际区域简化地质模型的地震数值模拟结果计算归一化滤波器算子的方法,从而利用该方法压制实际地震数据中多次波.     

井中垂直双极源体积分方程法三维模拟研究

类比水平方向电偶极源和垂直方向磁偶极源电磁场求解的原理,又利用偶极沿电流方向积分是双极原理,推导出地下均匀半空间中垂直双极源电场和磁场解析计算公式.采用体积分方程法对井中垂直双极源电磁响应进行了三维模拟.它对提高井中垂直双极源圈定储层范围的精度有着非常重要的意义,同时也为井中垂直双极源三维反演问题奠定了基础.     

青藏高原土壤热交换的计算方案及初步分析

从地温满足的热传导方程出发,导出了计算土壤平均和瞬时热通量的计算方案。该计算方案可同时计算出土壤热通量随时间和随深度的变化。它利用整层的地温信息来计算任一层的热通量,这种方案克服了用差分方案进行直接计算的局限性。然后使用中日亚洲季风观测实验期间的地温自动观测站资料和相应的常规观测资料,计算了青藏高原上土壤热通量及其变化,结果表明,不论是常规观测站还是自动观测站,其结果与青藏高原第一次观测实验所用热流板的直接观测结果是相近的,因而这种计算方案是实用而有效的。     

固体—液体层分界面处的P、SV型面波

针对半无限弹性固体—液体层分界面处 Ｐ、 Ｓ Ｖ 型面波不同的数学表达进行了分析研究,认为在有些著作中将它们统称为瑞利波是不够严密的,指出瑞利波和斯通利波是固体—液体层分界面处 Ｐ、 Ｓ Ｖ 型面波的两种极限情况。     

共偏移距可分表示法叠前深度偏移

基于双平方根方程的共偏移距可分表示法叠前深度偏移用于复杂介质成像, 该方法在中点-偏移距坐标中同时向下延拓炮点和检波点波场, 实现方式采用正反傅立叶变换.构造的双平方根方程波场延拓算子能够使波数域变量与空间(速度) 域变量分离, 波数域内进行相移计算, 在空间域对因介质横向变速引起的时移作修正.地震数据偏移不需要逐炮计算, 具有较高的效率.对Marmousi模型数据的偏移成像结果显示, 该方法较好地成像强横向变速介质中的复杂构造.      

波动方程角度域共成像道集

针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.     

湿润气候区无资料站点土壤湿度插补及预报试验

利用前期降水指数法(API)并结合水量平衡方程提出一种土壤湿度的初始化方案,将所建模型用于湿润气候区无资料站点的土壤湿度插补与预报试验,证明其应用的可行性。     

基于稳定逆时传播算子的黏声介质最小二乘逆时偏移

基于GSLS模型黏声介质二阶拟微分方程,采用伪谱法进行数值模拟。针对黏声介质逆时传播过程中产生的高频不稳定问题,提出加入规则化算子对其进行消除的方法,构建了稳定的逆时传播算子。在最小二乘反演的基础上,将黏声介质逆时偏移与最小二乘思路相结合,发展了带有振幅补偿的黏声介质最小二乘逆时偏移( LSRTM)。 Marmousi模型结果表明:相对于常规最小二乘逆时偏移,黏声介质最小二乘逆时偏移校正了地层的黏滞性,得到了更加精确可靠的保幅成像剖面。     

低煤级煤储层游离气含量计算——以准噶尔盆地东南缘为例

我国低煤级煤层气资源量大,约占煤层气资源总量的43.5%。由于对低煤级煤层气赋存特征的认识程度有限,影响了低煤级煤层气的勘探开发。通过对准噶尔盆地南缘低煤级煤储层孔隙与裂隙、吸附特征、含气性等方面的分析,认为该区煤的吸附性能较好,煤中宏观裂隙与显微裂隙发育。相对于中、高煤级煤,该区煤储层大、中孔所占比例较高,为游离气赋存提供了场所。运用气体方程估算了准噶尔盆地东南缘西山窑组B煤组主力煤层中的游离气含量,得出煤层总含气量为2.85~8.94 m³/t,平均为6.12 m³/t。其中游离气占总含气量的2.89%~5.14%,平均3.90%。游离气含量的估算为研究区更加科学合理的进行煤层气勘探开发提供了依据。     

基于改进的EAL边界条件的地震波场数值模拟

在地震波场数值模拟过程中, 人工截断边界会产生数值反射而影响地震波场模拟的精度.应用PML(Perfectly Matched Layer)吸收边界条件是消除人工边界反射的有效手段, 但是该方法计算量较大, 实现较为复杂, 且在某些各向异性介质中不稳定.EAL(Effective Absorbing Layer)边界条件是通过简化PML吸收边界条件的阻尼项推导出来的一种吸收边界条件.EAL的吸收项作用于波动方程的时间导数项, 故EAL边界条件具有计算量小且利于并行实现的特点.应用EAL边界条件能够在保证全波形反演和逆时偏移精度的前提下, 提升计算效率并降低内存消耗.本文对EAL边界条件进行两点改进: (1)引入减速阻尼项; (2)采用高斯型空间加权函数构建阻尼项.正演模拟结果表明, 改进的EAL边界条件在不增加计算量和内存消耗的同时, 有效提升了EAL边界条件的吸收效果且适用于复杂介质和各向异性介质.特别地, 相比于PML吸收边界条件, EAL边界条件在各向异性模型中表现出良好的稳定性.