利用τ-p变换技术实现多波波场分离

通过研究多波地震资料空间三个分量的波场特征及在τ-p域的特性,采用τ-p变换技术,编制了波场分离程序;选取适当的空间和时间τ-p参数,对模型和实际多波地震资料进行波场分离,取得了满意的效果。     

去假频高分辨率τ-p变换法波场分离

τ-p变换是消除相干噪声、波场分离的重要方法之一,但常规τ-p变换存在算子假频和端点效应,分辨率太低,无法直接使用。本文提出一种新的非迭代高分辨率τ-p变换方法,可同时压制假频和端点效应,提高了变换域的分辨率。通过理论模型记录与实际资料的试验表明了该方法的有效性。     

τ-p变换的两种实现方法

τ-p变换是一种经典的函数投影变换方法,在地震数据处理中有着广泛的应用,可以根据地震数据射线参数的差异,实现信噪分离、地震道插值、平面波分解等。但是,由于τ-p变换的精度和分辨率受到相应数学反问题的限制,在保证较高变换精度的前提下,τ-p变换计算速度的加快,以及相应滤波算子的设计等都值得研究。本文提出了基于径向道变换的τ-p变换方法和基于斜率分解的τ-p变换方法,阐明了两种方法的基本原理并与常规τ-p变换方法进行对比。基于径向道变换的τ-p变换方法利用快速Fourier变换和径向道变换,能有效减少τ-p变换的耗时;基于斜率分解的τ-p变换方法运用斜率分解算法,能进行高分辨率τ-p变换,并且提供冗余变换维度（τ-x-p域）,使滤波器算子的设计更灵活。数值实验结果表明,两种实现方法分别在计算速度和重构精度上优于传统τ-p变换算法。通过面波噪声压制的实际数据处理和比较,证明本文提出的两种τ-p变换实现方法可以为实际处理提供更加有效和灵活的实施方案。     

非零偏VSP多波波场分离方法

这里研究了中值滤波法在波场分离中的应用,分析了中值滤波中参数的选择对分离效果的影响,并应用中值滤波法分离上行波、下行波.同时,对VSP中纵波、横波的偏振方向进行了分析,利用上行波、下行波振幅矢量关系对VSP记录进行纵波、横波分离,得到良好的分离效果.理论模型处理结果表明,应用中值滤波分离上行波、下行波,再利用偏振关系分...     

基于τ-p变换的地震道重构

基于τ-p变换原理,研究一种基于τ-p变换的地震数据重构方法,运用该方法对待重构道用零道替代,通过τ-p正反变换及一定的迭代处理实现地震道数据重构。理论模型与实际资料应用表明,基于τ-p变换的地震道重构方法效果明显,重构数据稳健。     

基于双相介质弹性波波场分离的逆时偏移方法

随着地震勘探的发展,以单相介质模型为基础的声波逆时偏移和弹性波逆时偏移很难比较精确地描述地下介质和满足复杂的工程要求,尤其是在处理地下含流体多孔介质时,所以需要引入双相介质模型。首先基于Biot理论双相介质一阶速度-应力弹性波方程,利用散度场和旋度场的波场分离理论将固相和流相中的纵横波场解耦分离,并由波场快照图得出固相与流相中的纵波波场与横波波场均完成分离,而快纵波与慢纵波相互耦合均存在于纵波波场中。采用高阶交错网格有限差分方法,结合PML吸收边界条件,推导了具有波场分离的双相介质弹性波正传与反传的偏移延拓算子,结合互相关成像条件,可以得到双相介质弹性波固相与流相中的各分量偏移成像剖面。最后对模型使用该双相介质弹性波波场分离的逆时偏移方法进行偏移成像试算并取得成功,效果较常规双相介质弹性波逆时偏移成像效果有了明显提升。     

基于稀疏离散τ-p变换的叠后地震道外推

这里研究的是一种基于稀疏离散τ-p变换的地震道缺失重建方法。该方法根据叠后零偏移距剖面,在局部时窗内可以看作是一系列线性同相轴的组合,使用稀疏离散τ-p和预条件双共轭梯度算法进行地震道外推,使空间方向的缺失道得到恢复。τ-p变换由于信息不足,有限的孔径和离散等因素,导致在τ-p域的结果存在假象、不准确,而在τ-p域进行稀疏对于重建缺失的信息是非常有用和必要的。同时,还可使用迭代的预条件双共轭梯度算法进行计算。     

基于散度和旋度的弹性波波场分离数值模拟方法

弹性波波场分离数值模拟方法是研究混合波场中纯纵波和纯横波波场传播规律的一种重要手段,通常采用算子分裂来实现,这不但增加了计算量,而且还浪费了大量计算时间.因此提出了基于散度和旋度的弹性波波场分离数值模拟新方法,仅依赖于交错网格来计算出散度场(纯纵波)和旋度场(纯横波),从而提高计算效率.在人工截断边界处采用完全匹配层吸收边界条件(PML).数值模拟结果表明,该方法能够从混合波场中准确地分离出纯纵横波波场,且精度高,边界吸收效果好,计算速度快,计算过程稳定,可用于指导实际多波多分量地震资料处理,同时指出本文算法对相互耦合的波型(如面波)无法实现准确分离.     

面向海上地震资料的τ-p域保幅处理方法

AVO分析和反演需要地震资料处理时最大限度地保持振幅相对关系,然而常规处理中球面扩散补偿、拉伸畸变和多次波均会影响振幅相对关系,海上地震资料尤为严重。提出一套面向海上地震资料的τ-p域保幅处理方法。该方法通过τ-p域无需球面扩散补偿且多次波周期性增强的特点,减少球面扩散补偿与多次波对振幅相对关系的破坏,通过τ-p域拉伸因子非时空变换的特点减少动校拉伸产生的频率和能量畸变。最后通过挪威海上实际资料的处理验证了方法的合理性与可行性。     

2.5维各向异性介质中多波波场正演模拟与波场分析

完全3维弹性波数值模拟计算时间长,并且占用庞大的计算资源,这不利于在计算机配置不高的情况下进行科学研究,而二维弹性波数值模拟又达不到三维模拟的精度;同时,当模型、波场空间分布比较复杂时,传统的3维波动方程拟谱法模拟结果比较差.因此,在较高数值精度的一阶应力-速度弹性波动方程的基础上,采用傅氏变换仅计算y方向的偏导数,利用有限差分方法计算x、z方向和时间的偏导数,即利用2.5维数值模拟方法,实现在二维介质中计算三维弹性波场.最后通过数值模拟实现了在各向异性介质中多波波场的数值模拟,验证了2.5维方法是一种高精度、高效率、且能适应复杂模型的正演模拟方法,通过波场分析进一步认识了波在各向异性介质中的传播规律.     

基于修正伪Helmholtz分解算子的横向各向同性介质纵横波波场分离

纵横波波场分离是弹性波逆时偏移成像中的关键步骤,可以有效消除纵横波串扰和成像畸变,提高成像精度。在各向异性介质中,常用的波场分离法有非平稳滤波器或者低秩近似的方法,但由于使用多次傅里叶变换,导致波场分离计算成本较高。借鉴基于Helmhotz算子构建各向同性介质纵横波解耦方程的思路,提出消除振幅畸变的修正伪Helmholtz分解算子,通过待定系数方法求解横向各向同性介质中P波和S波2种模式的表达式,并将其变换为一阶解耦的伪弹性波方程,实现时空域纵横波波场解耦。通过简单模型的波场分离测试,得到分离后的P波与S波波场,验证了该方法的波场分离的有效性。此外,将解耦得到的矢量P波与矢量S波应用到弹性逆时偏移之中,利用矢量点积互相关成像条件得到清晰的弹性波逆时偏移成像结果,说明本方法在复杂介质中有着较好的适用性,同时可以有效地应用到VTI介质弹性逆时偏移之中。     

基于抛物线方程的坐标拉伸τ-p变换的P-P、P-SV波分离

分离P-P波和P-SV波最常用的是τ-p变换法, 但是在τ-p域中, P-P波和P-SV波常有重叠部分, 不能很好分离。地震勘探中, 当炮检距小于界面深度时, P-P波和P-SV波的时距曲线方程都可近似为抛物线方程。针对这种情况, 提出了一种基于抛物线方程的坐标拉伸τ-p变换的P-P、P-SV 波分离方法, 该方法首先对地震记录沿空间方向上进行坐标拉伸, 再进行τ-p正变换, 在τ-p域中分离P-P波和P-SV波, 其次对P-P波和P-SV波分别进行τ-p反变换, 实现P-P波和P-SV波的分离。与常规τ-p变换法对比分析表明, 利用该方法能很好地将合成地震记录中的P-P波和P-SV波分离。     

合肥市浅层地温场分布特征及影响因素分析

通过对合肥市15个地温监测孔的长期观测,分析了该市浅层地温场的分布特征及其影响因素。研究结果表明:合肥市浅层地温场整体呈西低东高分布特征,区内恒温带平均温度在17.4℃左右,恒温带深度在12～32m左右,地温梯度在2.3～4.1℃之间。浅层地温场主要受地质构造控制,同时受地层岩性及地下水补径排条件影响。200m以浅岩性以泥岩、砂岩为主,胶结程度差,结构松散,且地下水贫乏,浅层地温能的开发利用适宜采用地埋管型地源热泵。     

河南省城市浅层地温场特征及影响因素分析

本文在大量钻孔测温资料的基础上,系统分析了河南省城市浅层地温场分布特征,分析了不同地貌类型城市恒温带特征。全区地下水恒温带深度平均深度24.8m,温度一般15.5℃~17.5℃;冲积平原区松散层恒温带深度最浅、温度最高,内陆河谷盆地区松散层恒温带深度最深、温度最低。近山前地带基岩浅埋区,地温梯度低;沿深大断裂带和构造隆(凸)起区,地温梯度高;济源—商丘断裂的新乡—延津段、内黄凸起和通许凸起地温梯度高。通过分析地温增温率特征和地温恢复能力,得出颗粒越粗地温恢复能力K值较大,富水性越强、水力坡度越大K值越大。对影响浅层地温场的多种因素的系统研究表明,该区浅层地温场受城市、人类活动、地下水流场、地下水埋深、构造、地下水补给、排泄等因素影响明显。     

τ-P变换及向下延拓方法在声纳浮标折射地震测深速度反演中的应用

地震折射波场Φ（ｘ，ｔ）通过τ－ｐ变换得到倾斜迭加场φ（ｐ，τ）。从φ（ｐ，τ）中获得τ（ｐ）函数，并将其向下延拓到深度—慢度剖面，从而反演出横向均匀的速度结构。该技术能充分利用折射波场的丰富信息，在迭代反演中不仅能简单快速地得到反演结果，而且所得解具有良好的稳定性，较少受主观因素影响。将该方法应用于海洋声纳浮标地震折射测深数据处理，得到较为理想的结果     

雄安新区地热水的化学场特征及影响因素分析

【研究目的】雄安新区地热水化学场特征及其影响因素对了解区域地热资源聚敛机制具有重要意义。【研究方法】在区内5个主要构造单元（容城凸起、牛驼镇凸起、霸县凹陷、保定凹陷和高阳低凸起）开展地热水水化学采样和测试分析工作。【研究结果】研究区内主要分布有砂岩热储和岩溶热储,水化学类型以Cl-Na型和Cl·HCO₃-Na为主。从浅入深,地热水TDS值呈增大趋势。自西向东,地热水TDS含量表现为增大的趋势,指示着地下水流向。受断裂带导水作用影响,部分深部雾迷山组与馆陶组地热水发生混合现象。雾迷山组岩溶热储地热水的变质系数和脱硫系数低,热储封闭性较好,处于相对还原状态,而馆陶组和明化镇组砂岩热储地热水封闭性相对较差。【结论】区内地下水化学场主要受水岩相互作用程度、断裂、储层封闭性等因素影响。研究成果对于认识雄安新区地热资源、科学合理开发利用地热能,推进北方地区冬季清洁取暖具有重要意义。     

压制截断效应的Radon变换在波场分离中的应用

由于离散运算,Radon变换方法中存在截断效应,传统的最小二乘约束反演并不能有效抑制这种效应。介绍了一种频域空间稀疏约束算法,在反演迭代过程中,根据前一次迭代的结果,通过贝叶斯原理将加权矩阵与前一次迭代的结果联系起来,得到新的加权矩阵;然后求解这个加权矩阵方程,得到频率域的稀疏解。比较了用最小二乘反演和压制截断效应的Radon域变换的效果,后者效果明显优于前者。模拟了一个5层的水平地层的地质模型,用这种压制截断效应的Radon变换进行了VSP资料的上下行波场的分离,得到了较好的效果。     

基于小波变换方法的济南市区泉水影响因素研究

济南市区泉水的喷涌主要受大气降水和地下水开采等因素的影响,为查明降水量和地下水开采对济南市区泉水的影响,选取济南市区、济南西郊和济南东郊的水源地十年开采量数据,以及同期市区泉水位动态观测数据和降水量数据,采用交叉小波变换、小波相干谱和多元小波相干谱的方法进行分析。结果表明：①市区泉水位与降水量、地下水开采量均存在约1 a的频域周期,市区泉水位滞后于降水量133.22 d,市区泉水位与济南西郊开采量和济南东郊开采量的响应时滞分别为125.43 d和83.85 d,济南市区泉水与济南西郊地下水、济南东郊地下水存在水力联系;②济南市区泉水位与降水量的平均小波相干值（AWC）为0.58,市区泉水位与济南西郊开采、市区开采和济南东郊开采的AWC分别为0.47、0.40和0.32,PASC分别为13.90%、16.81%和10.09%,在研究期开采条件下,降水量和岩溶地下水开采量对市区泉水位的影响明显,市区开采对市区泉水位的影响最大,济南东郊开采对市区泉水位的影响最小;③在两种影响因素的作用下,降水量和济南西郊开采量、降水量和市区开采量两种组合可以作为市区泉水位变化的最大影响因素;在三种影响因素的作用下,降水量和济南西郊开采量、市区开采量组合可以成为市区泉水位变化的最大影响因素。     

位场希尔伯特变换的理论及应用

本从位场理论出发，推导了连续变量条件下任意形体位场的希尔伯特变换，导出了几种规则二度体磁场的希尔伯特变换解析表达式。在总结磁场垂直分量、水平分量及其希尔伯特变换之间的关系的基础上。讨论了由离散希尔伯特变换值反演地质体形状、确定座标原点，曲线切点座标的方法原理。为实现利用位场资料用计算机连续、自动反演产状参数提供了一种新的可行途径，基于大量模型实验结果，在辽东湾油气远景评价项目中，应用希尔伯特变换方法反演了磁性岩深度，结合其它地质、物探资料的综合解释，对判定生油沉积厚度，追索断裂、构造单元走向、圈定含油气远景区，提供了有意义的信息。     

基于多模式分离的S变换与小波变换提取面波频散曲线对比分析

为了对比分析小波变换和S变换在计算瑞雷波频散曲线时的精度,通过设计两层介质和六层介质模型,采用交错网格高阶有限差分方法对两个模型进行正演模拟,利用F-K法分离出时间—空间域不同模态的瑞雷波记录,再采用小波变换和S变换计算出不同模态的瑞雷波频散曲线,分别与理论值进行对比分析,最后利用两种方法对实际资料进行处理。结果表明:利用小波变换计算得到的频散曲线在低频段误差相对较大,S变换计算得到的频散曲线精度总体高于小波变换。