频率波数域预测和减去法压制多次波

研究了压制与海面有关的地震勘探多次反射波的两种方法。两种方法均以波动方程为基础。当地震子波已知时，可在频率波数域中通过子波反褶积方法消除与海水表面有关的多次波；在未知子波的情况下，可依据消除后波场能量最小为准则，求出一滤波因子，用它去标定估计的多次波，然后从地震记录中减去预测的多次波，达到压制它们的目的。为了降低维数，x、y的二维空间坐标用柱坐标的径向距离代替，其波数通过汉克尔变换求得。     

均匀横向各向同性介质中P波及SV波的射线规律

常规的反射数据处理都假定介质为均匀且各向同位。研究表明，各向异性的存在会严重影响处理和解释的结果。在已有知识的基础上，研究了二维均匀横向各向同性介质中的射线传播规律并描述了速度各向异性波场的特征。对于准P波，各向异性参数ε，δ的绝对值决定波前面的形状，而ε与δ之差决定射线密度；对于准SV波，σ则是主要的控制因素。     

湘鄂西—渝东褶皱带多重滑脱构造变形特征与构造模型

湘鄂西—渝东褶皱带为多层系、多类型天然气大型聚集区,揭示其结构与成因将为拓展该区勘探新领域提供重要依据,为其它陆内强构造变形区的油气评价提供重要参考。本文利用两条拼接的跨研究区的地震反射剖面,以断层相关褶皱几何学、运动学分析方法为指导,对湘鄂西—渝东褶皱带进行精细构造解析。研究结果表明,湘鄂西—渝东褶皱带构造变形主要与4套区域滑脱层有关,具有多层次拆离、多级别滑脱、多方向冲断及多样式组合的特征。纵向上由深及浅依次以低幅度构造、双重构造和冲褶构造为主要特征,空间上自南东向北西总体为由基底滑脱向盖层滑脱、由单层滑脱向多层滑脱的陆内递进式变形。湘鄂西—渝东地区存在多套“源—储—盖”组合,基于多重滑脱构造变形系统及“相对稳定”保存单元概念,认为未来常规天然气勘探可重点关注寒武系盐下领域,非常规页岩气需继续加强对寒武系、志留系、二叠系构造弱变形区的评价。     

叠前共炮点道集的奇性反演研究

讨论了地震波场的Fourier积分算子奇性反演理论与方法,得出常数背景假设条件下二维共炮点道集奇性反演的解析表达式．构造了两个理论地震模型,分别进行二维有限差分正演与奇性反演的数值计算将二维Fourier积分算子奇性反演方法用于野外的二维地震勘探叠前共炮点记录．在勘探范围内,当介质波速变化相对背景波速符合小扰动假设时,理论模型与实际地震数据的计算结果表明这种奇性反演方法的有效性和实用性     

各向异性介质中波动方程有限元法模拟及其稳定性

利用六面体单元和三线性插值函数，推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程．给出用时域有限差分求解该有限元方程时，时间、空间采样率和弹性参数应满足的稳定性条件．提出避免每一时间步都求解一大型线性代数方程组的方法．在开始求波场前，先求质量矩阵的逆的因子形式，以后每时间步只需作刚度矩阵的逆的因子形式与向量的乘积运算，可大大提高有限元素法在微机上进行波动方程数值模拟的运算速度．结合两个数值模拟算例，说明该方法的可行性，并与伪谱法作了计算速度的比较，有限元法的计算速度高于伪谱法的计算速度．     

各向异性介质相速度、群速度与克利斯托费尔方程

本文详细回顾了克利斯托费尔方程(Christoffel)的推导.在此基础上,提炼出时空域和频率域克利斯托费尔方程,前者可作正演计算介质的相速度,后者可计算介质的群速度.并计算了横观各向同性介质的相速度和群速度.最后对无损耗各向异性介质的相速度和群速度的关系进行了简单讨论.     

准噶尔盆地西北缘前陆冲断带二叠纪逆冲断裂活动的沉积响应

笔者从160口单井剖面、15条联井剖面及5个层位平面渐次展开了准噶尔西北缘前陆冲断带二叠纪逆冲断裂活动的沉积响应研究。各种成因扇体(冲积扇、水下扇、扇三角洲)在佳木河组和上乌尔禾组地层中最为发育，其沉积分布严格受不同时期活动的同生断裂控制。由早二叠世佳木河期至晚二叠世上乌尔禾期，扇体在克拉玛依-夏子街地区叠置程度较好，由盆缘向盆内迁移、推进明显；在车排子中拐地区叠置程度相对较差，局部略呈向盆缘老山收缩、后退趋势，表明红山嘴车排子断裂带与克拉玛依-百口泉断阶带、乌尔禾-夏子街断褶带逆冲活动的强度及地域迁移性有明显差异。二叠纪总体为一由弱到强的前展式推覆冲断及渐进式扇体迁移模式，即随同生控扇断裂由老山向盆缘的前展式推覆活动，相应地，各构造带扇体逐渐由盆缘向盆地方向发育、推进，扇体面积不断扩大，显示出明显的渐进式迁移的沉积响应，两者形成良好的耦合性。进而引入“活动性指数”的概念与方法，对前陆冲断带同生断裂的冲断活动强度进行了定量化统计分析，据二叠纪冲断推覆事件的地层、沉积标识，识别出3个逆冲推覆幕、7个逆冲推覆事件，划分出P1j-P1f、P2x—P2w、P3w三套构造层序，并指出在P1f、P2x及P2 3w各时期，由红山嘴-车排子→克拉玛依-百口泉→乌尔禾-夏子街构造带，同生断裂活动性趋于增强；在乌尔禾-夏子街断褶带，由早二叠世佳木河期到中、晚二叠世乌尔禾期，冲断推覆强度逐渐增大。     

双谱拾取时差参数场

研究Siliqi和Meur等提出的高密度拾取均方根速度和非椭圆率参数η的方法。导出了以最大偏移距对应的剩余时差dtn和零偏移距时间τ0为变量的时移双曲线方程,采用相似度扫描拾取dtn以及τ0,然后再把它们换算成均方根速度和非椭圆率参数。在理论数据处理过程中,发现时移双曲线扫描的精度比非双曲线方法有了一定的提高,特别是η值的误差减小了很多,这说明时移双曲线扫描的优越性。     

二维横向各向同性介质中人为边界反射的吸收——差分法弹性波场模拟

本文对有限差分法弹性波场三分量正演模拟中遇到的人为边界反射问题进行了研究,分别获得2-D横向各向同性介质中U_x,U_z,与U_y位移分量模拟时的边界吸收条件.从吸收边界上反射系数与入射角之间的关系可以说明,本文提出的边界吸收条件的吸收能力比较强.     

2021年日本福岛县冲地震的震源破裂过程分析基于采用经验格林函数方法的波形反演

利用中小震作为经验格林函数,选取0.2~2.0 Hz频段的强震数据进行波形反演,获得了2021年福岛县冲地震的破裂过程。结果表明:该地震的破裂主要集中在断层面距离震源约25 km的区域内,沿震源向东北延伸约5 km,向西南延伸约20 km;在该区域内,识别出两个滑动量集中的区域,均分布在距离震源约15 km西南侧,主要滑动量集中区域最大滑动量约3.2 m,深度基本与震源一致;次要滑动量集中区域略比主要滑动量集中区域浅约18 km。该震源模型对应矩震级M_w7.3,破裂速度为2.4 km/s。通过选择不同的中小震组合进行波形反演,结果对该震源模型特性无显著影响,揭示了该震源模型的稳健性;基于该震源模型合成反演分析中未使用的强震观测台站强震动,获得的合成波形与观测波形有很好的相关性,充分证明了该震源模型时空特征的可靠性。