基于程函方程快速行进法的地震走时层析成像法

地震走时层析成像法是一种利用不同地震台站连续接收到的地震事件,利用地震走时断层扫描反演地球内部速度分布的方法,主要分为正演(射线追踪)和反演两个环节。首先正演利用基于程函方程的快速行进法,实现了对入射波走时场和反射波走时场的计算,接着进行了射线追踪与走时灵敏度矩阵的计算,最后反演采用了子空间算法实现了反射波走时反演,反演成像结果可以分辨出异常体的形状和位置。     

三维地质体对AMT二维反演的影响研究

为探究三维地质体对二维反演结果地影响以及二维反演中极化模式的选择问题,设计了一系列二维、三维地电模型进行正演计算,并对正演结果进行二维反演计算。研究结果表明:相位受到三维畸变的影响较小,视电阻率较大;当三维异常体模型走向延伸较小时,TE模式的反演结果误差较大,TM模式的反演结果相对较好,但是反演得到了虚假的下覆构造,TE+TM联合反演结合了TE和TM模式反演的优点,可以较好地反映异常体的赋存位置及电阻率值;当三维异常体模型走向延伸长度增大时,三维数据的二维反演结果是可靠的,用TM和TE+TM模式反演都比较合理。     

直流电阻率三维正演的预条件共轭梯度法研究

三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。     

井中重力矢量数据联合反演研究

井中重力矢量测量可以近距离感知地下目标体,获得地质体不同方向的重力异常特征,提高地下介质的纵向分辨率.这里首先基于点元法实现了井中重力矢量正演,通过模型正演计算,分析了井中三分量重力异常响应特征,然后基于相关性搜索的黄金分割算法实现了井中重力矢量联合反演,通过模型反演分析了井中重力不同分量数据联合反演的效果,验证了反演...     

基于L-BFGS算法和同时激发震源的频率多尺度全波形反演

全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,庞大的计算量和存储空间需求,限制了全波形反演的发展。在频率多尺度全波形反演中将L BFGS数值优化算法与同时激发震源技术相结合的方法来改善这一现状。首先,对Marmousi模型进行了速度反演：在计算过程中明显发现对计算机内存的占用减少,最终反演结果与实际Marmousi模型的拟合误差为0.095 9,较小;采用10个频带单炮震源正演384炮所需时间约为32 640 s,而采用同时激发震源（384炮）正演一次所需时间仅约为700 s。然后,基于高速楔形体模型进行了抗噪能力研究：原始含噪地震记录信噪比为11.147 3 dB;对反演得到的速度模型进行正演,其地震记录信噪比为22.251 8 dB。最后,基于逆冲断层模型进行了反演速度扰动能力研究,反演得到的最终模型很清晰,与具有速度扰动特性的实际模型非常接近,拟合误差仅为0.036 0。数值模拟试验结果表明：此方法反演精度高,内存开销较小,能够显著提高计算效率,并且具有良好的抗噪能力,能够反演出具有速度扰动特性的介质。     

一种层析成像的正演和反演方法

作者在本文中提出了一种层析成像的正演和反演方法，首先对求解区作网络划分，由给定的正演速度模型在网格点上双线性插值，获得整个网格节点的速度分布，然后，对波阵面作一近似，在此基础上，用最小走时最短路径射线追踪法完成正演，反演则是通过正演的方法实现的，理论模型算例表明，本方法效果较好。     

二度重力人机结合选择法程序

一、功能本程序用于以人机结合选择法解二度重力勘探的反演问题。程序包括正演计算、绘图打印、物质带等效源反演计算三部分。正演计算采用水平多边形棱柱公式,可以有多个地质体,各地质体可以有不同的剩余密度。地形可以是平坦的或起伏的。在每次正演计算前,绘出本次计算所依据的地形地质剖面图。正演计算后,绘出异常剖面图。然后,由用户酌情选     

解重磁反演问题的人机结合选择法

(一)选择法的简要评述选择法是解重磁反演问题的一种老方法。它的做法是根据某一剖面上所观测的异常的特点,假定地下地质体的形状和物性参数;对这个假设的地质体做正演计算,并将计算的异常与观测的异常进行比较;如果它们之间相差很大,就修改原来所假定的地质体模型,     

有限单元法三维电阻率最小二乘反演中存在问题的研究

针对有限单元法在三维电阻率最小二乘反演中关于深部单元的反演能力、反演稳定性和速度问题进行了讨论与改进.首先构造了全新的"体积因子"作为先验信息加入到电阻率反演中,使地质体赋存区域的电阻率变化更符合实际情况,提高了深度较大的地质体的分辨能力.反演计算中对实测电阻率数据和模型正演数据分别取对数,并改进迭代修正量的计算方法,将以前的加减修正量改为与修正量指数之积,这样,不但可以改进电阻率模型出现负值的缺陷,提高反演计算的稳定性,而且还可减少迭代次数;针对反演时间过长的缺陷,将三维全域反演改为三维局域反演,使每次迭代的时间由4 h减少到目前的15 min;对模型的反演计算表明,以上改进是非常有效的.     

直流电阻率测深中二维与三维反演结果的对比与分析

直流电阻率测深中,用二维反演程序对三维地质体进行了反演,并与三维反演结果进行对比和分析。首先,对比二维和三维最小二乘反演在正演模拟算法和先验信息的确定方法;然后,对若干比较典型的模型进行反演实例对比。由于二维程序反演仅考虑单个剖面的电阻率信息,无论是在异常位置、形态及电阻率特性上,其反演精度都比较低;三维电阻率反演综合了多个测深剖面的电阻率信息,其结果与实际模型吻合得非常好。     

层析反演静校正

表层模型层析反演是一种非线性模型反演技术,它利用地震初至波射线的走时和路径反演介质速度结构,不受地表及近地表结构纵横向变化的约束。根据正演初至时间与实际初至时间的误差,修正速度模型,经反复迭代,最终达到要求的误差精度。求取静校正时采用射线法计算炮点到检波点的旅行时,从而得到基准面校正量。     

频率域反射波全波形速度反演

全波形反演不仅利用相位和振幅信息,还利用波形的细节变化,具有刻画模型精确细节的能力.在对稀疏矩阵直接LU分解求解的基础上,采用梯度预处理方法对声波介质速度模型进行了反射波全波形反演.采用误差反向传播算法计算目标函数梯度以及伪Hessian矩阵的对角线元素来做梯度预处理.数值模型的实验结果表明,利用有效的频率段便能反演出分辨率较高的速度结构,用低频反演出的结果作为高频反演的初始模型,减少了解的非唯一性.二维高斯光滑初始模型提供了有利的低频信息,得到较好的反演结果.伪Hessian矩阵的预处理吸收了高斯牛顿法的二次收敛优势,在不增加计算量的前提下,加快收敛速度.      

共聚焦点道集偏移速度建模

偏移速度建模可在多种道集实现.基于等时原理和差异时移分析,提出了共聚焦点道集偏移速度建模方法.偏移速度的更新是利用约束参数迭代反演实现,模型的参数化主要依据实际情况而定.为适应横向变速,选用了垂直速度梯度参数.模型试算表明:偏移速度相对误差在0.25%范围内,反射层深度误差小于10m.     

基于分段波形互相关的井下随采地震数据成像

随采地震勘探是以采煤机为震源的被动地震探测技术，由于采煤机是不断移动的、且其激发的是一种连续信号，因此常规数据处理方法无法直接应用。提出分段波形互相关方法，通过将采煤机产生的数据分段，采用互相关方法提取有效信号的走时，然后利用速度层析成像方法对工作面内部和切眼前方进行速度成像。利用波动方程对含异常体的工作面模型进行随采地震数值模拟，反演得到的成像结果和速度模型基本吻合；实际随采地震数据测试中，利用该方法对采煤机震源数据反演得到采动过程中地震波传播速度的成像结果，实现了对工作面内应力异常变化区域的实时动态监测。研究结果表明:基于分段互相关成像方法能解决采煤机震源信号处理问题，满足随采地震勘探技术实时性和稳定性的要求。      

粘弹介质叠前四参数同步反演及应用

叠前地震资料中包含纵波速度、横波速度、密度和吸收衰减系数等弹性及粘弹性信息,其中吸收衰减属性对储层物性及流体性质尤为敏感.本文从粘弹介质精确Zoeppritz方程出发,基于介质分解理论通过公式推导及近似,建立了包含纵波速度、横波速度、密度和吸收衰减系数的纵波反射系数特征方程,利用贝叶斯反演框架,实现了叠前纵波速度、横波...     

纵横波速度反演在气层识别中的应用

提出一种部分叠加资料的反演方法。部分叠加资料保留了原有的大量信息,在测井及构造层位的约束下,建立最符合地下实际情况的初始模型。通过非线性最优化理论,采用逐道外推技术,能同时反演出每一个点更准确的纵波和横波速度,进而可以求得纵横波速度比及泊松比,为储层预测和储集层内的流体识别提供丰富的弹性参数信息。     

初至波地震层析技术及其在四川复杂地区的应用

近地表速度场的重建是获得准确静校正量的关键。根据实际地震数据的初至波时间建立初始速度模型;采用改进的最短路径射线追踪算法进行正演;利用改进的约束带阻尼联合迭代重建技术(CDSIRT)求解大型、稀疏的方程组来反演近地表的速度场。理论模型和实际资料的反演结果表明,该方法能稳定、快速、准确地重建复杂地区的近地表速度场,建立准确的近地表速度场,从而可计算准确的静校正量,保证复杂地区的构造成像。     

井间地震走时波形层析成像方法

提出了井间地震走时波形顺序反演方法。该方法先利用井间地震走时反演得到速度模型的低频成分 ,然后用井间地震波形反演获得速度模型的高频成分。数值模型试验和实际应用结果表明 ,该方法反演稳健 ,提高了走时成像的分辨率 ,克服了波形成像易于陷入局部极小的缺陷 ,实现了快速高分辨率成像。     

井间地震时频域联合速度和衰减系数同步层析反演

井间地震资料中蕴含丰富的地下储层岩性和物性信息,利用井间地震直达波初至信息,通过层析反演可以得到两井之间的速度剖面。通过对井间地震直达波振幅衰减信息的处理,可以实现井间的吸收衰减层析成像。基于井间地震直达波射线路径的一致性,提出了井间地震时域直达波走时层析反演和频域质心频率衰减层析同步反演方法,增加了约束条件,增强了抗干扰能力。对模型和实际资料的处理结果表明,速度和吸收衰减同步反演方法,提高了速度层析反演精度,得到了对油气更为敏感的衰减属性,增强了井间地震在储层描述和油气预测中的应用能力。