海洋CSEM的共中心点域快速反演

海洋可控源电磁法(CSEM)对海底高阻体的反映比较灵敏,可用于天然气水合物探测资料的定性解释和反演研究.海洋CSEM资料的共中心点(CMP)域转换方式,可在横向上较好地分辨高阻储层.本文提出在CMP域实现一维频率域海洋可控源电磁资料高斯-牛顿反演算法.鉴于一维反演是解释地球物理资料的基础,较于二维和三维反演方法有着更高的计算效率和更低的硬件要求,将二维模型的响应在CMP域单元内表达为一维模型的响应,进而运用一维高斯-牛顿反演解释二维海洋CSEM资料.模型数据试算表明,海洋CSEM的CMP域反演速度较快,能够实现二维CSEM资料的反演解释.     

一维大地电磁信赖域反演法研究

信赖域算法是一种线性化的最优化算法,该算法在每次迭代时强制性地要求新的步长不超过某一控制量,这等价于在以当前迭代点为中心的一个邻域内对一个近似于原问题的近似模型求极值,从而提高了算法的稳定性.将信赖域算法引入到层状介质大地电磁数据的反演,得到了一维层状介质的信赖域反演法.数值试验结果发现该算法收敛性强,收敛速度快.抗噪...     

基于平均导数优化方法的VTI介质频率空间域正演

本文提出了一种新的基于平均导数优化方法(average-derivative optimal method,简称ADM)的二维VTI介质qP波波动方程频率空间域二阶9点格式,这种新算法将二维VTI介质qP波波动方程中中心空间导数项的差分近似表示为正交方向上3个网格点的加权平均形式.通过最小二乘优化方法求取空间导数项和加速度项的加权优化系数从而使数值频散达到极小化,每个波长所需要的网格点数在1%的误差范围内仅为3.57个网格点数,而VTI介质常规9点差分格式在相同的误差范围内则需要约12个网格点数,新方法的计算精度明显提高.复杂BP2007 2D VTI海洋标准模型数值模拟结果也验证了本文VTI介质9点ADM算法的有效性和准确性.     

交叉偶极各向异性快速处理新方法

在交叉偶极子横波各向异性反演中,通常先利用波形反演法得到各向异性目标函数,再采用二维的全局寻优方法来求取地层方位横波各向异性参数,但对于低各向异性地层,存在搜索效率低、低精度解的问题.本文结合Alford旋转法和波形反演法,以两个分别沿着方位轴和各向异性参数轴的一维搜索代替了原有的二维搜索方法,这种一维搜索方法使用布伦特(Brent)算法有效地实现对各向异性参数的求取.使用该方法处理了现场测井数据,我们获得了和原先的二维搜索方法相似的结果.新方法的优点是很大程度上提高了数据处理的效率和稳定性.     

频率域航空电磁法2.5维正反演研究（英文）

频率域航空电磁法具有经济、高效和适应能力强等特点,广泛应用于矿产普查、地质填图、地下水资源调查等领域。目前,实测数据处理技术以电阻率深度转换成像技术以及一维反演方法为主。然而,对二、三维复杂地电模型,一维反演往往难以达到满意的结果。虽然三维正反演能较好的解决该问题,但其计算效率难以满足海量实测数据的计算要求。因此,本文在前人研究的基础上开展了三维源二维地电模型(2.5维)的频率域航空电磁法正反演算法研究。将总场分解为一次场和二次场来消除源的奇异性,一次场在均匀全空间或层状介质中求得,二次场用等参有限元方法计算得到,利用大规模系数矩阵并行直接求解器计算线性方程组,有效提高了计算效率。在正演的基础上,实现了基于奇异值分解(SVD)的阻尼最小二乘反演算法,通过"拟正演"和互换定理来计算雅克比矩阵。通过理论模型计算和实测数据很好的验证了本文方法在频率域航空电磁数据正反演应用中的准确性和有效性。     

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.     

电偶源频率电磁测深二维阻抗视电阻率计算的源效应校正法

研讨了频率域电磁法中不同源装置的大地电磁测深、线源频率电磁测深和偶极源频率电磁测深阻抗视电阻率的源效应影响特征。在唯象分析的基础上，提出了几种电磁测深法阻抗视电阻率的相互换算法──源效应校正法（大地电磁测深二维ＴＥ极化视电阻率和其它两种电磁法的赤道装置二维阻抗视电阻车）。模型试验表明，利用这一源效应校正法可以由大地电磁二维视电阻率近似地计算出线源频率电磁二维阻抗视电阻率。这一方法被尝试应用于由线源频率电磁二维阻抗视电阻率估算偶极源频率电磁二维阻抗视电阻率。     

基于改进NAD方法的频率域声波逆时偏移

逆时偏移成像具有高分辨率的特点,在频率域实施能提升偏移成像效率.数值求解波动方程是逆时偏移成像的关键环节,数值算法的优劣直接决定着偏移过程的计算效率与成像结果的质量.本文在频率域近似解析离散化(NAD)方法的基础上通过优化网格差分模板系数构造改进NAD方法,保证数值离散精度的同时可缓解阻抗矩阵的病态程度,以此提高频率域...     

二维时间空间域和频率空间域声波全波形速度反演方法的对比研究

本文将二维时间空间域和频率空间域声波全波形速度反演方法分别应用到Marmousi模型,进行数值试验.两种方法均采用相同的观测系统和其他的参数,理论模型的数值试验结果证实了:使用较多的计算集群的CPU进行二维频率空间域直接法声波全波形反演时,其加速有限(正演数值模拟的计算量主要用于稀疏矩阵的LU分解,炮点计算波场时为线性关系).二维时间空间域声波全波形反演计算时更灵活,多炮同时计算时,可以多倍提高其计算效率;二维声波全波形速度反演时,直接法求解频率空间域的计算速度远快于时间空间域,所需要的计算机内存也比时间空间域少.二维声波全波形速度反演时,相比较于时间空间域的方法,频率空间域直接法声波全波形反演具有计算速度快和节省计算机内存需求的优势.     

可视化滤波器在加格达奇地震台震相资料分析中的应用

根据实际工作需要,为确保数字地震资料震相分析的准确性,基于MATLAB GUI设计了可视化滤波器。实现相应的数字低通、高通、带通、带阻滤波器,界面友好,具有可实用性。半年来,该滤波器在加格达奇地震台利用率高,发挥了较大作用。     

用传播矩阵法研究层状交错地层中的多分量感应测井

为了实现交错沉积等复杂环境中的电磁场数值模拟,本文在常规横向同性模型的基础上引入了电导率主轴坐标系相对地层坐标系的层理方位角和倾角,建立了交错地层模型.并利用传播矩阵法建立了一维层状交错地层模型中的多分量感应测井仪器响应的正演模拟算法.首先将频率-波数域中的电磁场分解为上行和下行模式波,给出了任意朝向的磁偶极子在无限大地层中模式波的解析解.进一步通过引入地层界面上的透射、局部反射以及广义反射系数矩阵,推导了一维层状地层中的模式波表达式.在此基础上,利用二维Gauss-Legendre积分实现了Fourier逆变换,得到了可用于多分量感应测井模拟的频率-空间域磁场并矢格林函数.最后,通过多个数值模拟结果考察了井眼倾角、层理方位角和倾角变化对多分量感应测井响应的影响.     

频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用

海上倾斜缆采集技术具有多样的陷波特征,通过去鬼波处理可获得宽频数据.针对海水面波浪起伏及缆深误差引起的鬼波延迟时间估计误差以及崎岖海底和目的层深度变化使得鬼波和一次反射波的振幅差异系数随偏移距的变化而难以给定一个固定值的问题,本文推导出频率慢度域中鬼波滤波算子以及自适应迭代反演求解上行波算法,该鬼波滤波算子与不同水平慢度对应的鬼波和一次反射波的振幅差异系数以及鬼波延迟时间有关.并基于计算出的理论下行波与实际下行波之间的平方误差最小理论实现自适应反演迭代最优计算该振幅差异系数和鬼波延迟时间.合成的及某海上采集的倾斜缆数据去鬼波处理结果表明,频率慢度域自适应迭代反演算法能较好地去除海上变深度缆鬼波,能达到拓宽地震记录频带目的.     

瞬变电磁法2.5维有限元数值模拟研究

本文对瞬变电磁法2.5维有限元正演进行了研究.从频率域麦克斯韦方程组出发,经过傅里叶变换推导出了走向y方向频率域电磁场响应的变分问题,然后运用频率域和时间域的转换公式求解出时间域瞬变电磁场的解.在求解频率域电磁场响应时为提高精度采用了基于二次插值的高阶有限元的算法,即单元网格插值为二次函数,同时推导出了经有限元离散后的泛函问题;在求解时间域电磁响应采用了正余弦变换的数字滤波算法.通过基本模型的正演,验证了算法的可行性.同时,也对比了基于G-S变换的线性有限元算法的数值结果,结果表明,本文采用的算法精度更高,层状模型最大延迟采样时间提高到了100 ms以上.     

2.5维起伏地表条件下频率域地空电磁正演模拟

频率域地空电磁探测方法是指在地面布设人工场源,在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势,尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大,因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算,分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题,利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分,源和异常体附近网格加密处理,保证计算精度,远离目标区域网格逐渐稀疏,模拟无穷远边界,降低对计算资源的要求.在单元内进行插值,将有限元方程组变换为线性方程组,采用总场算法,利用具有一定面积的伪δ函数表达源电流分布,源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应,再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解,验证了本文算法的精度,同时本文还对地空电磁场与地面电磁场的响应特性进行了对比.     

频率域航空电磁数据的加权横向约束反演

传统的一维反演技术已经被广泛应用于航空电磁数据解释中.然而,利用单点水平层状介质模型模拟地下复杂地电结构有时会遇到困难.突出表现在反演参数的横向不连续性,即使相邻测点的反演结果也会出现突变.本文针对航空电磁直升机吊舱系统可进行密集采样,相邻测点地下电性结构应具有某种程度连续性的特点,研究航空电磁数据横向约束反演理论,并提出参数加权约束方法.首先阐述频率域航空电磁正演和加权横向约束反演理论,着重介绍这种拟二维反演方法的基本原理和实施步骤,以及将该方法成功应用于频率域航空电磁数据反演处理的方法技术.最后,通过对理论和实测数据反演处理,并与传统的一维反演结果进行对比,验证加权横向约束反演方法的有效性.     

基于形态学灰度重构的乳腺微钙化点提取

针对乳腺X线摄片,提出了一种计算机辅助检测乳腺微钙化点的方法:首先利用灰度形态重构完成穹顶提取,除去乳腺背景结构同时又能很好地保留微钙化点;其次运用拉普拉斯算子在空间域内对图像进行锐化以增强微钙化点;然后在频率域内利用巴特沃思高通滤波器对图像进行高通滤波;最后利用二维最大熵阈值分割对图像二值化,从而检测出微钙化点并在图像上标识出来。计算机辅助诊断的方法不仅避免了医生因大量阅片引起疲劳而造成误诊,起到了协助医生诊断的作用,而且提高了诊断的效率和准确率。     

Walsh列率域中多维分形模型与GIS环境下地球物理信号处理

地球物理信号通常在多个尺度段表现尺度不变性,这些不变性起因于不同的地质、地球物理或成矿过程的自相似性. 利用这种在多个尺度段的尺度不变性可以设计多维分形滤波器,滤波所得信号表征了具尺度不变性的地质地球物理或成矿过程,可以用于成矿预测或环境评价. 本文研究了Walsh变换列率空间地球物理信号的列率功率谱密度与列率之间的分形与多维分形关系, 试验证实了大洋钻探、石油以及煤系地层地球物理测井资料在Walsh域的多维分形性质,提出了用于分解地球物理场,提取有用信号并用于矿产资源勘探或环境评价的多维分形W-A模型. 利用波列率域中的多维分形关系构造了W A图解(W-A Plot). 借助W-A图解可以确定最小平方误差（LS）意义下Walsh功率谱变化的不同自相似性的频率分界点,从而用于设计W-A分形滤波器. 这种滤波器可将地球物理场分解成具有不同自相似性的局部场,并且保留原场的各向异性结构. 与通常使用的基于Fourier变换的滤波技术相比,W-A模型具有许多优点：W-A适用于检测地球物理信号中的突变、线性、环状、局部与纹理结构等弱信号. 同时,由于Walsh变换中只有简单的变号(加法与减法),其计算速度远快于建立在复数乘法之上的Fourier变换,所以W-A计算速度远快于Fourier域的滤波方法,可以用于地球物理信号的现场实时处理. 用加拿大Nova-cotia省西南地区的布格重力异常进行了W-A方法的试算,处理结果反映了地质、矿产分布规律,能够很好地进行矿产预测.     

共炮检距高斯波束偏移的一种快速实现算法

与共炮高斯波束偏移相比,共炮检距高斯波束偏移具有直接抽取炮检距域共成像点道集的优势.过去,共炮检距高斯波束偏移以损失成像精度的代价采用最速下降法来降低积分的维数,从而提高计算效率.但经过最速下降近似简化的偏移公式仍是频率域的,需要在每个频点进行计算.为此,本文提出一种快速实现算法来避免采用最速下降法.本文通过分析一个水平层状速度模型的偏移过程和Marmousi速度模型的成像结果来检验不同插值方法对快速实现算法的成像精度和计算效率的影响,并建议采用二维三次卷积插值方法.同时本文在Marmousi速度模型下验证了快速实现算法相对于最速下降法在成像精度和计算效率上的优势.此外,本文将采用二维三次卷积插值的快速实现算法应用于Sigsbee2A模型并获得了清晰的盐下图像.     

快速精确的二维频率空间域弹性波数值模拟

二维频率空间域的数值模拟方法具有以下的优势:多炮模拟时,计算成本比时间域方法低;无累计误差;在地震反演中处理多震源模拟时,只需要有限的几个频率就可以得到好的反演结果.差分离散化形成的稀疏系数矩阵,需要求解一个巨大规模的线性方程组,最大瓶颈是需要海量的计算机内存,导致计算量庞大.本文在前人研究的基础上,采用嵌套剖分网格排序法,极大限度减少对计算机内存的需求,从而减少了计算量.针对弹性波数值模拟的特征,提出二维频率空间域弹性波多炮模拟的快速计算流程.数值模拟试验证明使用嵌套剖分排序法的弹性波多炮数值模拟比压缩存储法具有节省存储量、计算效率高等优势,为后续的二维频率空间域弹性波全波形反演奠定了很好的基础.     

用传播矩阵法研究层状正交各向异性地层中多分量感应测井响应

本文采用传播矩阵技术研究并建立了层状正交各向异性地层中多分量感应测井响应的有效算法.首先通过Fourier变换将频率空间域中的Maxwell方程组求解问题转化为频率波数域中关于电磁场水平分量常微分方程组的定解问题.利用该方程组系数矩阵的本征值和归一化本征向量将电磁场分解成上行波和下行波模式的组合,推导出均匀正交各向异性介质中由任意方向磁偶极子产生的电磁波模式解析表达式;在此基础上,利用叠加原理和边界条件研究了电磁波在层状正交各向异性地层中的反射和透射,给出各个界面上的广义反射系数和不同地层中电磁波振幅的递推公式,进而得到电磁波模式的解析解.为了有效确定频率空间域中的电磁场,采用二维Patterson自适应求积算法结合有限连分式展开技术计算傅氏逆变换.最后通过数值模拟结果证明了该算法的有效性,考察了不同各向异性系数、不同井眼倾角以及仪器长度和工作频率变化等情况下的多分量感应测井响应特征.