基于慢度法的小波变换多尺度叠前弹性波反演方法

针对常规的非线性反演方法存在稳定性较差及易陷入局部极值的问题,采用小波变换将原问题分解在不同的尺度上进行迭代反演,并以大尺度搜索为起点,逐步缩小尺度,最终降至目标函数的原始尺度,同时通过界约束以增加解的稳定性.对理论模型和实际资料分别进行的叠前弹性波反演结果表明,所提方法能有效提高反演的稳定性及精度.     

叠前弹性参数反演方法的关键技术

为实现对岩性和隐蔽油气藏的有效勘探,提出一种叠前弹性参数反演方法。该方法以Zoeppritz方程近似式为反演基础。通过比较不同介质条件下的Zoeppritz方程近似式的计算精度,确定适合的方程近似式。反演时,采用无拉伸动校正技术处理模拟的正演地震数据,使用射线追踪法抽取方程计算所需的角度道集。反演测试结果表明,相对于传统的叠后反演方法,该方法可以反演各种弹性参数,且反演结果的横向连续性强,可用于储层预测。     

VTI介质三维VSP多波联合叠前AVA岩性参数反演方法及应用

三维三分量(3D 3C) VSP属于地下地震勘探方法,三维三分量VSP地震资料合有深度和丰富的地层岩性信息,有地面地震不可比拟的优势.为了发挥和应用三维三分量VSP资料的优势,提高油气勘探的成功率,需要反演出准确的岩性及弹性参数,用于精确预测储层、识别流体性质.以各向异性介质(VTI)弹性波传播理论为基础,研究多波精确振幅特征方程(AVA),并在此基础上建立了三维各向异性介质VSP多波联合叠前AVA反演岩性参数理论和方法.该方法利用对3D 3C VSP资料处理的NMO道集,直接反演地层纵波速度、横波速度、密度及2个各向异性系数,在此基础上可进一步计算得到地层的弹性参数,这些参数对岩性及流体识别有重要作用.这套纵横波联合叠前反演方法适应直井及斜井3D 3C VSP资料,并经实际斜井3D 3C VSP资料处理,得到了与实际储层吻合的反演结果,证明了该方法的正确性和有效性.     

基于叠前弹性参数的剩余油描述方法及应用

中后期老油田开发如何精准描述剩余油在井间的分布成为亟待解决的难题.文章基于长垣油田典型区块,利用其高密度地震资料,开展基于叠前弹性参数反演的剩余油预测应用技术研究.建立了密度、体积模量和剪切模量等弹性参数与含油饱和度的关系,采用Zeop-pritz公式为核心的弹性参数反演算法,得到了密度、体积模量和剪切模量等弹性参数,实现了L1井区含油饱和度分布预测.经由后验井验证,含油饱和度描述准确率达75％以上.剩余油整体分布相对零散,呈甜豆状分布,在河道相砂体边部含油饱和度较高,指导措施调整后取得了较好的挖潜效果.     

利用叠前弹性参数反演进行储层含气性预测

常规的叠后波阻抗反演没有考虑地震反射振幅随入射角变化的情况,反演结果只有声波阻抗一种参数,信息单一,很难进行储层预测。基于AVO原理发展起来的叠前反演技术可以同时预测地层的纵波速度、横波速度和密度。虽然这些与岩性和流体有关的岩石物理信息比波阻抗更为丰富,但其他一些重要的岩石物理参数(如拉梅系数和剪切模量)只能由这3个参数间接计算得到,这样不可避免地引入了额外的计算误差,并且降低了信噪比。为了减小这种误差,利用Zoeppritz方程的Gray近似公式来拟合振幅随入射角的变化特征,采用直接反演出拉梅系数、剪切模量和密度的叠前弹性参数反演方法,大大提高了反演结果的准确性。使用模型进行验证,并将该方法应用于中国西部某地区的含气储层预测,取得了较好的效果。实际资料的应用表明,该方法能用于储层含气性预测。     

基于混合遗传算法的叠前随机反演方法

针对常规随机反演方法计算效率低的问题,提出一种基于混合遗传算法的叠前随机反演方法。该方法充分利用测井资料中的高频信息,并以地震数据作为约束,首先通过快速傅里叶滑动平均(fast Fourier transform-moving average,FFT-MA)谱模拟算法进行随机模拟得到基于地质统计学的初始模型信息,随后结合提出的混合遗传算法对模拟结果进行快速优化,得到符合地下地质结构的反演剖面,实现对叠前弹性参数的反演。混合遗传算法避免了一般遗传算法常见问题,如收敛速度慢以及产生"早熟"现象,与模拟退火相结合能够快速收敛达到全局最优,保证了反演精度。数值试验结果表明,融入混合遗传算法的叠前随机反演方法,在充分利用叠前信息的同时可以保证反演结果有效收敛,并且与模型数据吻合较好,与传统的叠前反演方法相比具有较高的分辨率,在储层识别和油藏描述中起到了重要作用。     

半空间弹性波强散射双参数非线性反演

研究了半空间内弹性波强散射源的逆问题,利用矩量法进行离散处理,并利用正则化原则避免求解病态矩阵方程,基于DBI(Distorted-Born Iterative)方法,提出了一种非线性反演算法．和Bom方法不同的是：该算法在每一步迭代过程中不仅只修正总波场,同时也修正格林函数．数值模拟计算结果证明了本方法的有效性。     

基于叠前反演流体识别与频散AVO反演的研究

由于纵波不能在流体中传播,叠后地震反演得到纵波波阻抗是无法识别储层中的流体。基于横波在流体中传播特点,利用叠前数据丰富的运动学与动力学信息获得弹性参数进行储层流体识别是可行的,但受到地震反问题固有的“病态解”等不确定性的影响,在进行叠前地震反演时会产生较为严重的多解性与不确定性,从而影响了后续的储层岩性与流体识别的精度。基于以上问题本论文以储层埋藏较深、钻井资料较少的D2地区为研究目标,采用测井岩石物理敏感分析、敏感流体因子反演和频散AVO反演三种技术综合来提高储层预测的精度,实际资料的应用效果表明采用多种技术方法来进行储层的流体识别是可行的。     

层状半空间弹性参数非线性迭代反演

对层状半空间弹性参数的反演这一具有较强的工程应用价值的问题进行了研究。先通过引入共轭波函数给出了泛函导数的表达式,然后构造一均方误差函娄把反演问题转化为寻求最优解,并通过共轭梯度法导出了非线性迭代公式,最后进行了数值计算。数值计算结果表明,该公式具有较好的收敛性和稳定性。     

基于地震叠前反演的零散砂体岩体力学场建模

针对大庆油田外围扶杨低渗透油层分布零散, 岩体力学性质变化快, 储层改造难的问题, 在大量岩体力学原位模型测定的基础上, 研究纵横波速度之间的关系。从叠前弹性参数反演得到纵横波速度体, 利用岩石物理关系式得到井间的岩体力学参数场, 运用多元回归分析方法得到井点弹性参数值。结合井点与井间弹性参数值, 采用地质统计学方法建立非均质岩体力学三维场, 为水力压裂参数优化设计提供地质依据, 达到增产的效果。     

基于叠前振幅随角度变化反演的储层流体识别方法

从理论模型入手,首先研究流体参数变化对地震弹性参数及地震记录的影响,然后用胜利油田某区块的实际资料进行叠前振幅随角度变化(AVA)弹性反演.在反演基础上,以振幅随炮检距变化(AVO)分析为基础,提取反映流体特征的不同属性,在剖面与平面上指示出储层流体的特征.反演结果与实际钻井资料吻合率较高,说明叠前反演与AVO属性分析结合是预测储层流体性质的一种有效方法.     

Nonlinear optimization method of radar target polarization enhancement

n the basis of nonlinear programming model and solution of radar target polarization enhancement,a novel efficient iterative algorithm is presented. The optimization effieiency and precision are improved. The perfor-mance of algorithm is analyzed and verified by experiments.     

基于最优Born近似的叠前深度偏移方法

波动方程法中的叠前深度偏移是实现复杂构造和岩性地震成像的关键技术。其中，广义屏法是近几年发展起来的一种新方法，它基于散射理论、屏近似和一般Born近似，具有较高的精度和效率，是条件稳定的。为了提高算法的稳定性，提出了基于最优Born近似的叠前深度偏移方法。即对散射波场成像所使用的介质速度横向扰动公式不做任何近似，并且采用与两个延拓深度层中间位置处的地震波场最接近的波场代替层间的平均波场，而不是用上一延拓层处的地震波场代替层间的平均波场。模型试算结果表明，最优Born近似法对复杂地质体具有较强的适应性，在保证计算效率和成像精度的前提下，明显提高了算法的稳定性。     

基于最优Born近似的叠前深度偏移方法

波动方程法中的叠前深度偏移是实现复杂构造和岩性地震成像的关键技术。其中 ,广义屏法是近几年发展起来的一种新方法 ,它基于散射理论、屏近似和一般Born近似 ,具有较高的精度和效率 ,是条件稳定的。为了提高算法的稳定性 ,提出了基于最优Born近似的叠前深度偏移方法。即对散射波场成像所使用的介质速度横向扰动公式不做任何近似 ,并且采用与两个延拓深度层中间位置处的地震波场最接近的波场代替层间的平均波场 ,而不是用上一延拓层处的地震波场代替层间的平均波场。模型试算结果表明 ,最优Born近似法对复杂地质体具有较强的适应性 ,在保证计算效率和成像精度的前提下 ,明显提高了算法的稳定性。     

基于算子优化的Fourier有限差分法保幅叠前深度偏移

从声波方程的单程波分解出发,对传统波动方程叠前深度偏移的传播算子和边界条件进行优化,推导出基于优化算子的保幅边界条件和3个新偏移算子。通过应用新算子和边界条件,断层和倾斜层处的成像效果有明显改善,并且叠前深度偏移的质量也得到提高。模型试算和实际资料处理结果表明,使用优化的边界条件和偏移算子,能够提高横向速度变化情况下的成像精度。     

二维井间电磁场的非线性反演方法

采用高精度的非线性反演方法，对轴对称二维井间模型的电导率分布进行了成像。该成像方法既考虑到雅可比矩阵的线性项，也考虑了其偏微分项，因而反演精度高且稳定性强。采用矩阵求逆技术计算了纵向成层背景地层Green函数的系数，可以很方便地得到Green函数对空间坐标的积分，从而使雅可比矩阵的计算更方便、准确。反演实例表明，井间电磁成像技术更适宜于对相对高电导率异常区域进行成像。由于井间电磁成像技术对异常区块的垂向探测范围大于径向探测范围，测量数据中包含的垂向信息多于径向信息，因而垂向分辨率高于径向分辨率。     

二维井间电磁场的非线性反演方法

采用高精度的非线性反演方法,对轴对称二维井间模型的电导率分布进行了成像。该成像方法既考虑到雅可比矩阵的线性项,也考虑了其偏微分项,因而反演精度高且稳定性强。采用矩阵求逆技术计算了纵向成层背景地层Green函数的系数,可以很方便地得到Green函数对空间坐标的积分,从而使雅可比矩阵的计算更方便、准确。反演实例表明,井间电磁成像技术更适宜于对相对高电导率异常区域进行成像。由于井间电磁成像技术对异常区块的垂向探测范围大于径向探测范围,测量数据中包含的垂向信息多于径向信息,因而垂向分辨率高于径向分辨率。     

Applying a dual optimization method to quantify carbon fluxes: recent progress in carbon flux inversion

The widely performed Bayesian synthesis inversion method （BSIM） utilizes prior carbon flux and atmospheric carbon dioxide observations to optimize the unknown flux. The prior flux is usually computed from ecological models with large biases. The BSIM is useful in solving the problem of insufficient data, but it will increase the inaccuracies in the estimates caused by the biased prior flux. In this study, we propose a dual optimization method （DOM） to introduce a set of scaling factors as new state variables to correct for the prior flux according to information on plant functional types. The DOM estimates the scaling factors and carbon flux simultaneously by minimizing the cost function. The statistical properties of the DOM, which compare favorably with the BSIM, are provided in this article. We tested the DOM through simulation experiments which represent a true ecosystem. The results, according to the root mean squared error, show that the DOM has a higher accuracy than the BSIM in flux estimates.     

高精度叠前时间域速度分析方法研究及应用

以水平层状为假设的常规速度分析和倾角时差校正都不能很好地解决倾斜地层反射点弥散问题,而高精度共散射点道集速度分析方法则可以很好地解决。研究优于常规方法的椭圆展开共反射点叠加速度分析方法,用它获得的速度场作为高精度共散射点道集速度分析方法的初始输入。同时,研究等价偏移距单平方根旅行时方程的推导,以及共散射点道集的提取过程及基于共散射点道集的速度分析方法与叠前时间偏移。通过两种方法的组合,从模型和实际资料的试算可以看出速度谱精度得到提高,获取得到高精度速度场,并且获取叠前时间偏移剖面的效率得到提高。