地震波场模拟中的褶积微分算子法

褶积微分算子是将偏导数算子的带限傅氏谱进行反傅氏变换得到的,具有对称或反对称和衰减速度快的特点。根据频谱分析理论,选用序列在结束时能光滑趋于零的截断窗口将褶积算子截断,结合考虑微分算子守恒定律,将褶积算子中心值进行修正,得到适合于波场模拟的最优权系数。给出用褶积微分算子模拟声波波场的实现途径与实例。理论分析与模型试算证明了褶积微分算子法的正确性与有效性。     

基于拟柱坐标系的2.5-D声波方程数值模拟

2.5-D 拟柱坐标法是一种新兴的2.5-D 问题的处理方法,该方法最大的优点在于拓宽常规柱坐标法数值模拟的应用范围,能够处理非轴对称模型,而且计算简便,尤其适宜解决大尺度地震勘探耗时多的难题.笔者发展了基于拟柱坐标系采用交错网格高阶有限差分法求解2.5-D拟柱坐标声波方程的数值模拟算法,并分别通过简单的单倾斜界面模型和较复杂的非轴对称模型的数值实例验证该方法的正确性和可行性.     

基于delta波包叠加的时间域深度偏移

delta波包可由高斯波束经傅里叶逆变换得到,是高斯波束在时空域的表达.它最早出现在合成理论地震图的研究中,本文将其应用于偏移领域.通过delta波包叠加表达时间域格林函数,可将高斯波束偏移由频率域转换到时间域,再结合Rayleigh积分和激励时间成像条件,本文给出了基于delta波包叠加的深度偏移算法.该偏移算法可在时间域直接计算,但因包含褶积运算,成像时将耗费大量的计算时间.针对这一问题,本文提出了将褶积简化为乘积的近似公式.近似后的偏移算法,不仅保留了高斯波束偏移的优点,而且计算效率得到显著提升.文中通过两个数值算例验证了上述结论.     

三维起伏地表条件下地震波走时计算的不等距迎风差分法

三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型.     

不同域的局部平面波分解应用与对比

局部平面波分解的计算精度和计算效率对高斯波数偏移存在较大的影响。笔者分析了目前常用的时间域、频率域、频率波数域3种不同域的局部平面波分解方法,给出了3种不同计算方式下局部平面波分解的具体计算公式;并结合局部平面波分解的特征,针对不同域下的局部平面波分解算法采取合理的程序设计思想, 对多组模拟数据应用不同域的局部平面波分解算法进行计算精度和计算效率的对比。通过对比分析可知, 应用频率波数域的局部平面波分解,不但提高了计算的准确度, 而且相对于时间域和频率域局部平面波分解能够提高近60%的计算效率,为下一步进行偏移成像更高效地提供了精确的数据。     

气候变化和人类活动在榆林市荒漠化过程中的相对作用

针对已有研究的不足,提出一个集成时序和空间两方面特征的植被-气候时空关系统计模型,并用来在栅格尺度上定量区分陕西省榆林市荒漠化动态中气候变化和人类活动的相对作用。结果显示,1986-2000年,榆林市荒漠化发展主要发生在南部丘陵沟壑区,逆转主要发生在北部风沙草滩区,这一空间格局特别明显且在研究期内没有发生根本改变。从第1时段(1986-1990年)至第2时段(1991-1995年),荒漠化整体上处于发展状态,且无论是荒漠化发展区还是逆转区,人类活动的作用均占据绝对主导地位,分别为98.7%和101.4%,气候变化的作用甚微;从第2时段(1991-1995年)至第3时段(1996-2000年),荒漠化过程整体上达到稳定态势,在荒漠化发展区,气候变化和人类活动的作用基本相当,分别为46.2%和53.8%,而在荒漠化逆转区,人类活动的作用达到了119.0%,气候变化的作用则为-19.0%。另外,气候变化和人类活动的相对作用有很大的空间异质性。