利用地震资料反演地层的碳酸盐含量

基于碳酸盐含量与地层速度、密度之间的关系, 在井资料约束下, 使用人工神经网络方法反演高分辨率地震资料所反映的地层碳酸盐含量, 并应用于南海北部陆坡ODP184航次1146和1148孔区, 取得较好效果.方法的关键是从井旁地震道中提取多种属性, 利用逐步回归法, 确定6种与碳酸盐含量相关性最好的地震属性, 分别是平均频率、道积分绝对振幅、主频、时间、道微分瞬时振幅和瞬时频率, 然后进行地层碳酸盐含量反演.反演结果相对于岩心分析的碳酸盐含量的误差大多在±5%之内, 较为精确地揭示了地震地层剖面上碳酸盐含量的分布.      

单斜介质中方位NMO速度Thomsen参数反演方法研究

采用启发式共轭梯度法,即随机爬山法 + 共轭梯度法,利用单斜介质中P波方位NMO速度椭圆轴向偏转角度接近于零这一特性,简化P波方位NMO速度公式,并利用多方位P波NMO速度,反演出某一初始CMP观测线与自然坐标系之间的夹角,作为进一步进行Thomsen各向异性参数反演的基础. 根据各向异性介质中方位NMO速度与Thomsen参数之间的关系,建立了利用三种波的多方位NMO速度及垂直传播速度反演单层单斜各向异性介质Thomsen各向异性参数的目标函数. 对计算的理论值添加具有一定标准差的正态分布的随机噪声,用以模拟实际观测存在的误差,通过对加噪后的数据进行多次反演的误差分析,表明了所建立的目标函数及选用的反演方法是有效可行的,而且相对稳定.     

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

改进的等效半空间法及瑞雷波频散曲线反演

层状模型的理论频散曲线表明,层状模型表面瑞雷波同一频率对应的相速度具有多值特点,这就是所谓的瑞雷波多阶性. 面波的多阶性导致了实测频散曲线的复杂性,当地层中存在软弱夹层时这一问题尤为突出,从而给实测频散曲线的解释带来了困难. 本文从工程实用的角度出发,根据等效半空间理论,提出了一种计算理论频散曲线的新算法——改进的等效半空间法,避开了面波多阶性这一复杂问题. 由此方法计算得到的面波相速度并不同于某个具体阶数的面波,而是对多阶面波的综合反映. 按照改进的等效半空间法编制了拟合反演程序,并在工程中进行了运用,取得较好的效果.     

电阻率二维神经网络反演

由于非线性特性地球物理反演一直以来都是一个比较困难的问题. 近十年来,非线性反演方法如人工神经网络、遗传算法在地球物理数据解释中得到越来越多的应用,但目前基本仍限于一维反演问题. 对于二维反问题,反演参数较多,神经网络反演运用较少. 本文利用BP神经网络优化方法,实现了电阻率二维非线性反演. 与传统线性化的迭代反演比较,神经网络反演能够克服传统方法的不足、获得更好的反演结果.     

利用GPS观测的时间序列资料反演地壳地幔黏性结构

地球介质的黏度是地球动力学研究中的重要参数. 在考虑震后震源区介质物性变化的情况下,本文提出利用GPS观测的时间序列资料反演地壳、地幔黏度的方法. 计算中,将遗传算法与Marquardt法相结合,利用遗传算法搜索结果作为Marquardt法反演的初始模型,运用微扰法形成Jacobian矩阵. 根据余震分布及地震波速度变化确定震源区的范围及介质的物性变化情况,利用黏弹性有限单元模型（Maxwell体）及台湾集集地震后GPS观测的时间序列数据,反演了台湾地区地壳、地幔的黏性系数. 反演计算及数值实验表明,该方法反演黏度的效率高,稳定性好. 反演结果初步显示台湾地区的地壳、地幔的黏度分别为12×1018Pa·s和36×1019 Pa·s.     

2004年3月18日高纬磁层顶多重通量管事件分析

2004-03-18 23:10～23:50 UT期间,“双星（Double Star）”探测一号卫星（TC 1）在向阳面磁层顶高纬晨侧由内向外穿越磁层顶,其时TC_1的GSM坐标为 (75RE, -55RE, -54RE), RE为地球半径.穿越过程中TC_1观测到了8个通量管和1个磁通量传输事件（FTEs）.在此期间Cluster星簇位于向阳面太阳风内,其GSM坐标为(180RE, -31RE, -62RE),其4颗卫星监测到行星际磁场（IMF）的BZ分量持续南向,BＹ有较大的负值.本文的研究表明：TC_1观测到的前7个通量管具有准周期重现性,周期大约是1～4 min,明显小于以前所观测到的FTEs的平均周期(8～11 min);所有的通量管都具有较强的核心场.本文分别使用最小方差分析法（MVA）和Grad_Shafranov反演方法（GSR）对通量管的轴向进行了分析和对比,发现所有的通量管主轴基本沿晨昏向,结果显示GSR方法在轴向分析上比MVA优越.本文使用GSR方法对通量管的磁场结构进行了分析,恢复出了通量管的磁场在卫星穿越面的结构图;此外,本文还对这次多重通量管事件进行了deHoffmann Teller(HT)分析,结果表明,所有通量管大致朝南极方向运动,均来源于向日面低纬区域.这说明它们可能起源于向日面低纬区,由该区的磁场分量重联产生.     

求解地震静校正问题的双尺度反演方法

通过对地震静校正问题的分析,综合线性反演方法的计算速度快和非线性反演方法寻找全局最优解能力强的优点,提出了求解地震静校正问题的双尺度反演方法.在大尺度下采用非线性反演方法为小尺度下的线性反演计算初始模型,利用小尺度下的线性反演方法获得精细模型.对反演参数的选取进行了讨论,使双尺度反演算法具有自适应的特点.理论模型和实际资料的计算表明该方法计算结果精度高、计算速度快.     

叠后地震属性分析在油气田勘探开发中的应用

地震属性分析技术一直是地震特殊处理和解释的主要研究内容.随着油气勘探开发的发展,地震属性分析技术已经成为油藏地球物理研究的核心内容,是勘探地震与开发地震之间纽带.本文针对鄂尔多斯盆地的低幅度构造、低孔隙、低渗透率、致密性隐蔽油气藏的特点,综合应用相干数据体分析、地震相自动分类定性识别砂体厚度、地震振幅属性分析、频谱分解、多井约束的储层叠后反演等叠后属性分析技术,探索了一套适合该区油气特征的储层横向预测及油气识别模式.为该区油气勘探开发,储量计算提供可靠依据.同时也为隐蔽性油气藏的勘探开发积累了经验.     

各向异性介质中的弹性阻抗及其反演

地震反演已成为油藏描述中的重要组成部分.绝大多数的常规地震反演是叠后地震数据体上进行,很少考虑各向异性存在的情况.随着勘探开发的发展,地震各向异性和叠前地震波阻抗反演引起了人们极大关注.本文在各向同性介质中弹性阻抗研究基础上,推导出了各向异性介质中的弹性阻抗方程,提出了地震各向异性介质中用弹性阻抗进行储层参数描述的技术路线和框架,并对反演过程中存在的问题进行了有益探讨.