小波变换多尺度地震波形反演

完全非线性地震波形反演问题是石油地球物理勘探领域中一个非常重要而又难度很大的问题。本文提出了多尺度地震波形反演的小波变换方法,对于一维非线性地震波形反演问题,将该方法和已有的简单迭代法及多重网格法作了比较,数值实验结果表明,本方法效果较好。     

小波变换多尺度地震波形反演

完全非线性地震波形反演问题是石油地球物理勘探领域中一个非常重要而又难度很大的问题。本文提出了多尺度地震波形反演的小波变换方法,对于一维非线性地震波形反演问题,将该方法和已有的简单迭代法及多重网格法作了比较,数值实验结果表明,本方法效果较好。     

中国重力异常的小波变换与多尺度分析

阐述了小波变换及多尺度分析方法的原理,利用二维小波多尺度分析方法对中国布格重力异常进行了分解,并从其中划分出有意义的剩余异常．模型结果证明方法有效．由于本文方法可以将重力异常分解成各种几何尺度意义下的不同成分,这为重力解释和研究地壳的结构提供更多更新的依据．     

重力异常和DEM的小波变换与多尺度相关分析

为深入认识高程与重力异常相关性的内在机制,利用小波变换信号分解特性,研究了重力异常和相应数字地形模型（DEM）在不同小波尺度上的统计相关性,由相关系数、乘常数、加常数的变化特点,得出传统意义上重力异常与高程的相关,实际上表现为二者小波变换小尺度细节的相关,这种相关不是简单孤立的关系,而是包含物理意义的平均相关。     

中国重力异常的小波变换与多尺度分析

阐述了小波变换及多尺度分析方法的原理，利用二维小波多尺度分析方法对中国布格重力异常进行了分解，并从其中划分出有意义的剩余异常．模型结果证明方法有效．由于本文方法可以将重力异常分解成各种几何尺度意义下的不同成分，这为重力解释和研究地壳的结构提供更多更新的依据．     

基于多尺度有限元法的频率域航空电磁三维正演研究

正演模拟是电磁数据反演的基础, 其计算速度与精度一直是制约电磁反演的两大核心问题.在三维电磁正反演中, 传统方法通过加密网格或增加插值基函数阶数提高计算精度, 但由此也降低了计算效率, 制约了三维电磁反演的实用化.因此, 如何实现大尺度模型高精度快速正演是目前电磁三维正反演中亟需解决的问题.本文将多尺度有限元法应用到麦克斯韦方程求解中.我们首先在粗网格尺度上构建满足局部特性微分算子的多尺度基函数, 进而在粗网格尺度上对原问题进行求解, 通过建立粗细两套网格间场的映射关系, 在未知数较少的粗网格上实现电磁问题求解之后, 利用粗细两套网格间场的映射关系获取细网格上电磁场响应, 由此可以在保证计算精度前提下快速获取不同尺度电磁场正演响应, 计算速度得到很大提高.此外, 本文还基于八叉树思想进行网格优化, 进一步改善三维正演效率.我们通过对典型地电结构进行多尺度有限元正演模拟并与传统有限元结果对比验证算法的有效性.最后, 我们通过模拟加拿大Voisey's Bay卵形体镍铜硫化矿区航空电磁响应以检验本文算法模拟地下复杂异常体的能力.

     

基于单程波正演模拟的多尺度地震资料特征分析

多尺度地震资料在油藏地球物理中发挥着重要的作用,通过开展单程波正演模拟并对模拟结果进行偏移,有利于深入剖析多尺度地震资料之间的联系与差异.正演模拟表明,地面地震资料、VSP资料和井间地震资料三种尺度的地震资料都是相同地下地质体对应于该观测尺度下的地震波响应特征,差异在于由于观测方式的不同,资料所包含的波场信息丰富程度不一样,反映地质体信息的范围不一样,在分辨能力上具有一定的差异性.地面地震资料具有较强的横向连续性但纵向分辨率低,VSP资料具备比较精确的时深关系,井间地震资料分辨率高但是覆盖范围小,对三种多尺度地震资料开展综合研究有利于更好的揭示地下地质信息.     

基于GPU并行的重力、重力梯度三维正演快速计算及反演策略

利用NVIDIA CUDA编程平台,实现了基于GPU并行的重力、重力梯度三维快速正演计算方法.采用当前在重力数据约束反演或联合反演中流行的物性模型(密度大小不同、规则排列的长方体单元)作为地下剖分单元,对任意三维复杂模型体均可用很多物性模型进行组合近似,利用解析方法计算出所有物性模型在计算点的异常值并累加求和,得到整个模型体在某一计算点引起的重力(或重力梯度)值.针对精细的复杂模型体产生的问题,采用GPU并行计算技术,主要包括线程有效索引与优化的并行归约技术进行高效计算.在显卡型号为NVIDIA Quadro 2000相对于单线程CPU程序,重力和重力梯度Uxx、Uxy正演计算可以分别达到60与50倍的加速.本文还讨论了GPU并行计算在两种反演方法中的策略,为快速三维反演技术提供了借鉴.     

基于面向对象思想的正演反演统一框架结构

运用面向对象思想分析地球物理正演与反演问题,提出了一种基于有限单元正演与正则化非线性最小二乘反演的统一框架结构.该框架采用稀疏线性系统求解作为正演的基础类;分别建立一维、二维、三维的几何模型类,用于完成模型的剖分与有限单元钢度矩阵的生成,其中与单元物质属性相关的计算采用纯虚函数的形式,在具体解问题中实现;以正则化非线性最小二乘最优化算法作为反演基础类,反演中有关模型的正演数据计算及Fréchet矩阵求解等采用纯虚函数的形式,在具体问题中进行计算.框架结构的建立提高了代码的重用性与程序开发速度;对地质单元采用统一的空间剖分,有利于开展联合反演程序设计.     

基于波动方程正演模拟的多尺度地震资料反射特征分析

地面地震、VSP和井间地震等多尺度地震资料在油气勘探开发中发挥着重要的作用,但不同的观测尺度和资料采集处理方式限制了多种资料之间的有效匹配和综合应用.通过对同一地质体模型采用波动方程开展正演模拟分析,进一步开展逆时偏移有利于更好地建立起多尺度地震资料之间的联系与差异.分析表明,地面地震、VSP和井间地震等不同尺度下的观测资料都能在对应的范围内准确反映地质体的地震反射特征,而观测尺度和观测方式的不同则导致多尺度地震资料在分辨率、波场信息丰富程度和观测范围等方面存在差异;正演模拟地震道的时频分析进一步表明多尺度地震资料在时间域和频率域都存在着较好的联系,特别是井间地震和VSP等井中地震资料能够提供更为合理、可靠的高频信息,在提高地面地震资料分辨率方面具有明显的优势,为开展多尺度地震资料的综合研究与应用提供了基础.     

台湾地壳构造的三维重力正演研究

用三维重力正演方法处理了台湾地区的布格重力异常和地质资料.根据由浅而深分层正演(剥皮)得出的剩余重力异常,发现某些断层的深度并不大;而另一些断层或构造则越向下越清楚.根据正演的结果,还讨论了台湾所处的地球动力学环境.      

河北省及邻区重力异常多尺度小波分析

静态重力场变化可为地壳深部结构研究提供重要基础数据.本文利用小波多尺度分解方法,对河北及邻近区域(110°～120°E,36°～42°N)的布格重力异常、Airy均衡重力异常数据资料进行了处理分析,得到了不同尺度下的分解结果.结果表明:张渤带、山西带在中上、中下地壳存在一定程度的异常,破坏性地震主要发生在重力异常正负差...     

用代数重建法作正演和反演模型计算

本文主要是软件技术开发工作,结合地震学的物理模型[3],给出不同的“数据”组,作正演计算(其结果相当于实测数据),然后用代数重建法(ART)作反演计算,进行数字图象重建,获得了较好的计算结果。     

基于自适应有限元正演的大地电磁法三维反演算法研究

目前大地电磁法三维反演大多基于有限差分正演.反演使用规则六面体网格无法有效模拟复杂地形, 同时正反演网格同套存在严重影响反演可靠性的问题.针对上述两个问题, 本文利用自适应有限元算法, 开发了基于非规则六面体的高精度大地电磁三维正演方法; 反演中, 使用独立的反演网格和正演网格来提高正演计算精度和反演可靠性.这一思路既保证了正演响应和灵敏度矩阵计算的精确性, 又降低了因反演参数过多造成的不唯一性.最后, 通过算例验证了正演算法的精确性和反演算法对地形处理的有效性.

     

基于反演基质矿物模量和多约束条件的双相介质AVO正演方法（英文）

基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。     

横向均匀介质中地震走时的正演和反演

横向均匀介质可由若干水平分层足够好地逼近,本文引入零厚度层来描述成层介质,从而建立了回折波、反射波、首波和直达波统一的时距关系以及地震走时随模型参数变化的公式。本文详细讨论了水平分层介质里地震走时曲线的特征,合成地震走时的试验表明了由地震走时,用阻尼最小二乘法可以成功地反演速度模型参数,所反演的模型含有速度梯度层,也可以含有低速层。     

青藏高原三维地壳结构的重力正演模拟

青藏高原地壳结构的重力学研究主要集中于莫霍面深度反演及单个剖面的重力正演,但是,在亚东-格尔木与下察隅-共和之间没有地震测线的区域,地壳结构的研究受到限制.为了探讨从青藏高原中部到东部的地壳结构特征,本文采用三维重力正演的方法,对青藏高原三维地壳结构进行了重力正演模拟.研究结果显示,地壳各层界面起伏较大,表明高原地壳变形强烈.研究区域中,藏南喜马拉雅地体界面相对较浅,处于界面陡变带,往北进入拉萨地体,界面相对变深.正演剖面中,大的界面起伏与大的断裂构造有关,莫霍面的陡变也暗示着莫霍面存在错断现象.东构造结莫霍面有抬升现象,在玛多附近莫霍面有下凹现象.地壳分层正演显示,32 km以下地壳产生的重力异常占总异常的80%以上,低速层及高速体产生的重力异常最大可占到总异常的10%左右.     

基于小波变换尺度相关性的轮廓数据去噪

针对机械零件CT切片轮廓数据存在大量噪声的特点,在总结小波变换方法的基础上,提出了基于曲率小波变换的多尺度分析方法,利用尺度间相关性达到去噪和特征点检测目的,不仅解决了单一尺度无法有效去除切片轮廓存在的冗余数据和噪声,而且由于线特征剧烈波动而导致轮廓特征点检测困难,拟合精度差等的问题也得到有效地解决。     

地球物理学中的电磁场正演与反演

本文在近年来众多的地球物理研究者的研究基础上,总结了当前地电磁模型正反演已有成果,定量分析了各种主要正反演的性能测试,指出不同正反演法的优点、缺点以及应用范围局限,提出了各种方法的发展趋势以及当前计算地球物理领域的核心内容,指出了计算地球物理领域的数值模拟发展方向.     

基于卫星重力的青藏高原东缘多尺度时变重力场研究

利用经过去相关滤波处理的GRACE时变重力场模型获得了青藏高原东缘2003—2012年的卫星年重力变化图像,并针对该区域近年发生的三次特大地震,结合震前及震后月重力场变化图像,分析与强震有关的卫星重力场变化特征。从区域年重力变化图像可以看出,三次大震均发生在年重力变化较低的时段内,震前小幅值变化可能是地震发生的中短期前兆;从汶川地震和玉树地震发生前后的月重力场变化图像可以发现,发震前后断层附近的重力变化模式发生变化,这可能印证了震后位场变化恢复理论;从汶川地震前后的龙门山断层附近点上的周重力变化趋势可以明显发现,汶川地震发生（第20周）后近9周的时间,断层东西侧呈现了相反的重力变化特征,这可能是对震后壳幔物质调整过程的反映。