双平方根波动方程偏移速度分析

传统的剩余校正（RMO）偏移速度分析方法基于走时原理,在陡倾角和欠照明地区,因为不能得到充分的角度域信息而失效.本文将展示一种基于波场延拓理论的偏移速度分析方法,即波动方程偏移速度分析（WEMVA）.这种方法先利用成像优化方法获得剩余成像,再利用剩余成像反演剩余速度.此类方法继承了波动方程偏移方法的优点和缺点.波动方程偏移速度分析是一种线性反演方法,它要求对Born近似的展开序列作一阶截断.高阶部分的丢失必然带来巨大的截断误差,因此剩余成像必须也进行线性化,以适应大速度扰动和大延拓步长.因此,在此类算法中,剩余成像的获取和线性化是偏移速度分析的关键.在叠前偏移算子中,因为双平方根算子的数学表达式更为简洁,所以本文基于对波动方程偏移速度分析初步讨论,并通过模型验证其原理.     

提高浅层地震反射静校正精度的一种方法与数值模拟

浅层地震反射法是一种常用的勘探方法.在浅层地震资料处理中,静校正的精度直接影响速度反演的结果和叠加剖面的质量,在地形平缓时,固定基准面静校正可以满足勘探精度的要求,但在复杂地形条件下,其存在较大误差,即使采用浮动基准面,仍会由于地表一致性假设而残余静校正量,不能消除地形起伏引起的影响,为了提高浅层地震反射静校正的精度必须在常规静校正后进行一次剩余静校正,本文给出起伏地形条件下,滑动基准面(过共中心点的水平面即为该共中心点的滑动基准面)的剩余静校正量,该校正量与炮检距、反射层埋深、地层波速以及炮点和接收点高程有关,适用于单一介质和层状介质情况,本文通过对典型地形起伏的3个水平均匀层状介质理论地质模型的速度谱计算和分析,阐明在复杂地形条件下,应用本文提出的剩余静校正方法可以消除地形起伏的影响,提高静校正精度,在此基础上做动校正可以得到高质量的水平叠加剖面.     

非对称走时克希霍夫叠前时间偏移方法在塔河地区的应用

针对塔河油田“复式”成藏组合所具有的油气圈闭类型复杂、储层厚度薄、构造幅度小、岩性横向变化大和非均质性强等特点,采用非对称走时克希霍夫叠前时间偏移方法进行数据处理.与常规叠后时间偏移和弯曲射线克希霍夫叠前时间偏移相比,非对称走时克希霍夫叠前时间偏移成果剖面绕射波归位合理,断层、断点清晰,不整合面成像及奥陶系内幕碳酸盐缝洞体反射结构更加清楚,横向分辨率进一步提高.     

川滇及邻区Pn波速度结构和各向异性研究

川滇地区是中国大陆构造变形及地震活动最强烈的区域之一.本文利用从国际地震中心(ISC)、中国地震台网以及云南、四川、重庆三个省级地震台网的地震观测报告中严格挑选出的53673个Pn波到时数据,用Pn波层析成像方法反演了川滇地区的Pn波速度结构和各向异性.结果表明该地区上地幔顶部速度结构的横向变化与现代构造运动有明显的关联:盆地地区上地幔顶部速度明显高于其周边区域,四川盆地尤为突出,其上地幔顶部P波速度是整个研究区域的高值;上地幔顶部P波速度低值区与研究区现代火山活动区域(如腾冲)及高热泉活动区域(如康定)有明显的相关性.强震大都发生在Pn波速低值区或Pn波速低值边界区.Pn波速度低速带与狭义的南北地震带相吻合,反映出上地幔顶部的流变性对构造应力的传递起着重要作用.Pn波各向异性表明:位于川青块体与四川盆地交界的龙门山地区,快波方向与该区域绝对运动方向不一致,而是基本沿龙门山断裂走向,这表明该地区的各向异性主要受龙门山断裂对Moho面切割变形的影响.而在其他区域的Pn波各向异性方向大体与当地块体的绝对运动方向一致,说明其各向异性主要是由板块运动产生的软流圈变形引起.     

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

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

利用远震研究四川盆地及其周缘瑞利面波相速度和方位各向异性

基于四川盆地及周边的245个宽频带台站2010年9月-2014年9月期间的远震记录,提取双台路径瑞利面波相速度频散资料,反演得到四川盆地20~120 s的高分辨率瑞利面波相速度及各向异性空间分布.在丰富区域地球物理基础数据的同时,结合已有研究成果对地壳上地幔变形耦合进行探讨,结果表明短周期（20~30 s）的相速度分布与四川盆地的地质构造特征相吻合,作为川滇地块、松潘-甘孜地块和四川盆地之间的边界——龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用;松潘-甘孜地块,特别是川滇地块中下地壳普遍表现为明显的低速异常,表明中下地壳相对软弱;而四川盆地的中下地壳整体呈现相对高速,表明四川盆地具有相对坚硬的中下地壳.研究区域东南角接近北扬子地块与南扬子地块的缝合部位,呈现高速异常.四川盆地南部和东南邻区不同周期均具有较强的各向异性,且快波方向较为一致,反映这些地区不同深度变形耦合较好.四川盆地西部、北部及东北部邻区,不同周期的各向异性快波方向变化较大,不同深度变形耦合较差.这些特征与绕喜马拉雅东构造结的物质流动被扬子地块的高速地壳阻挡的宏观认识基本一致.

     

鄂尔多斯西缘北段的地壳结构和块体间变形关系

本文对喜马拉雅计划二期部分台站的远震波形数据进行接收函数提取,利用接收函数共转换点叠加方法获得阿拉善地块、鄂尔多斯地块以及银川—河套盆地下方0~80 km深度的速度间断面结构.结果表明：鄂尔多斯地块成层性好,地壳厚度为38~42 km,康拉德界面为18~22 km,阿拉善地区的Moho面深度为38~45 km.河套盆地地壳厚度约52 km,银川断陷盆地和贺兰山下方的Moho面最深为~55 km.鄂尔多斯西缘构造边界下方Moho面变化明显,且黄河断裂为深大断裂直接切割莫霍界面.根据本文的间断面成像结果我们进一步确定阿拉善地块与鄂尔多斯地块分属不同的大地构造单元.与此同时,我们推测贺兰山以西70~80 km范围内和鄂尔多斯地块西缘北段存在地壳增厚变形的可能.

     

青藏高原东南缘Moho面速度密度跃变研究

青藏高原东南缘地下深部结构的研究对了解青藏高原的变形机制和动力学过程具有重要意义.本文利用四川、云南固定台站记录到的远震波形资料,首先采用接收函数H-k叠加方法获得青藏高原东南缘台站下方的地壳厚度和波速比.进而利用接收函数一次转换波和多次波幅度信息确定了青藏高原东南缘Moho面上的S波速度和密度跃变.研究结果表明：研究区由南到北地壳厚度逐渐增加,从永德、沧源、孟连地区的33 km左右增至巴塘地区的69.7 km左右,厚度变化了近乎37 km.四川盆地和松潘甘孜块体南部的姑咱地区具有高泊松比、速度密度跃变较小特征,表明这两个地区含有较多铁镁物质.腾冲地区、龙门山西侧的汶川地区、四川盆地西南缘的沐川地区以及则木河断裂的石门坎至东川地区同属于高泊松比、速度密度跃变较大,显示这些地区壳内存在部分熔融.

     

集成奇异谱分析和自回归滑动平均预测日本近海海平面变化

本文从日本沿岸选取了28个验潮站及联测的GPS站,利用奇异谱分析（Singular Spectrum Analysis,SSA）和SSA+自回归滑动平均（Auto Regression Moving Average,ARMA）方法预测了2014—2018年的近海海平面变化和地壳垂直变化.并用同时段的验潮及GPS的实际测量值进行验证,结果显示,SSA+ARMA预测的相对海平面精度为0.0357~0.0607 m,地壳垂直运动的精度为0.0049~0.0077 m,绝对海平面的精度为0.0433~0.0683 m,且三者SSA+ARMA的预测结果均优于只用SSA预测的结果.在此基础上本文利用SSA+ARMA预测了日本沿岸2019—2023年的近海绝对海平面变化,结果显示,2019—2023年的平均海面高较往年（2014—2018）升高0.0353 m,2003—2023年绝对海平面的变化率为0.0039 m·a^-1,预测结果较为理想.

     

利用爆破事件分析不同速度模型下常规定位方法定位结果

地震定位对速度模型的依赖性很强。四川地区地形复杂,常规工作中可选取多种速度模型进行定位。川西龙门山断裂带为东南部四川盆地和西北部青藏高原东部山区的明显分界线,近年在此断裂带上发生多次较大地震。对发生在该断裂带附近的6个爆破事件和15个天然地震重新定位,并对比结果。研究表明,相同台站包围情况下,川滇3D速度模型稳定性最好,但对浅表爆破不太准确。相比HypoSat(一维速度模型)组合,台站分布对Hypo2000(一维速度模型)和Hypo2000(赵珠速度模型)组合的定位结果影响较大。