简介：常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。

简介：传统上，有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高，但在这个范围之外会产生较强的数值频散；基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度，并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维，形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后，将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中，并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。

简介：储层物性参数的改变伴随着弹性参数某种程度的改变，而这种改变并非简单的线性关系。加上观测信息的缺失、重叠、噪声破坏及模型理想化等原因，致使这种改变关系的获取具有较大的不确定性。针对传统统计岩石物理物性参数反演方法的不足，以利用弹性参数反演储层物性参数为目的，依据贝叶斯反演框架，我们建立了新的储层物性参数目标反演函数。首先，采用兼具确定性与随机性特点的统计岩石物理模型，考虑到不同弹性参数间的精度存在差异，引入权重系数，建立起储层物性参数与弹性参数间的加权统计关系。其次，基于这种加权统计关系，结合马尔科夫链蒙特卡洛随机模拟技术产生储层物性参数、弹性参数随机联合样本空间作为目标函数求解样本空间。最后，建立解的快速求解准则，求取最大后验概率密度对应的储层物性参数取值作为最终解。实际应用表明，该方法具有较高的反演效率，应用前景较好。

简介：本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法（FastFourierTransformMovingAverage,FFT-MA）和逐渐变形算法（GradualDeformationMethod,GDM）得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。

简介：有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程，高阶显式有限差分方法通常用来提高求解精度，已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数，发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明：忽略这些小的差分系数能够近似维持有限差分的精度，但是显著减小了计算量。然后，引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性，采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后，给出了针对均匀和非均匀介质模型的弹性波模拟例子，表明了本文方法的优点。

简介：本文是三分量地震在陆相薄互层地震勘探中应用成功的首例。大庆长垣喇嘛甸油田处于高含水开发后期,储层预测的主要难题是密井网条件下厚度2m以上砂体的边界识别问题,常规单一纵波地震解释存在多解性,引入横波信息可增加地震解释的可靠性。本文根据纵、横波对气藏的响应特征,利用纵、横波联合振幅属性、分频属性对气藏进行了分析,能够很好地反应气藏的边界,体现了横波受流体影响较小的优势。通过测井曲线岩石物理定量的分析了该地区主要是密度、λρ和μρ三种物性参数能够很好地反映储层岩性变化。在以地震数据分频为核心的高精度地震参考标准层解释基础上利用纵、横波联合反演得出的纵、横波速度比、λρ、μρ等属性与测井解释的储层具有很好的一致性,多波属性λρ沿层切片描述了开发小层的砂体平面展布,与井砂岩图具有较高的一致性,刻画了井间砂体边界变化的细节,并指示了剩余油挖潜的有利区域。

简介：针对有限差分数值模拟的频散问题，本文将交错网格技术和紧致差分格式相结合，推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式；对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项，并利用Fourier误差分析方法分析了上述三种差分格式的近似精度；在此基础上，分别采用上述三种差分格式进行了波场数值模拟。结果表明，当差分方程阶数相同时，紧致交错网格差分格式截断误差最小，数值频散最弱，差分精度最高，证实了该方法的有效性。

简介：针对粘弹性裂隙型单斜介质,本文应用高阶交错网格差分法,对其传播快照以及多方位地面三分量记录进行模拟,结果证明了裂隙填充物的性质引起的各向异性与粘弹性性质对单斜介质波场有明显的影响,总结了不同观测方位地震记录呈现出的规律性,这为进一步了解波场在介质中的传播情况,提出更合理的理论模型,提供一种尝试。

简介：孔隙结构特征及类型划分对低渗透储层勘探开发至关重要,基于多重分形理论与核磁共振实验,对东营凹陷南坡沙河街组沙四段（ES4）复杂低渗透砂岩进行孔隙结构研究。首先,根据岩心物性、铸体薄片、压汞等资料所反映的孔隙结构参数差异,将研究区的岩石孔隙结构类型划分了3大类5小类;然后,针对不同类型岩石的核磁共振T2谱进行插值并计算其对数坐标下的一维、三维分形维数以及多重分形谱,并提取多重分形参数奇异性强度α、分布稠密度f（a）,结果显示孔隙结构类型不同,盒维数尤其是多重分形参数值差异明显,孔隙结构好,其a、f（a）偏向高值,以此划分孔隙结构类型与压汞、薄片分析结果基本一致;最后,将该方法应用到核磁共振测井剖面上,应用效果较好,表明多重分形是核磁共振T2谱的一种属性,利用核磁测井T2谱多重分形特征及参数能够连续较好地评价低渗透砂岩孔隙结构类型与预测有效储层。

简介：分数阶S变换（FRST）具有较强的时频聚集性。利用FRST处理地震数据，通过合适的分数阶参数将频率轴旋转到适当位置，即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同，因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体，并利用共频率数据体进行了低频伴影分析，然后提出FRST和盲分离结合的方法，不需要对地震数据的最优分数阶参数进行估计，即可提取识别有效地质特征信息的独立频谱，提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据，并与已知井信息进行比对，验证了其有效性。

简介：分数阶S变换(FRST)具有较强的时频聚集性。利用FRST处理地震数据,通过合适的分数阶参数将频率轴旋转到适当位置,即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同,因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体,并利用共频率数据体进行了低频伴影分析,然后提出FRST和盲分离结合的方法,不需要对地震数据的最优分数阶参数进行估计,即可提取识别有效地质特征信息的独立频谱,提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据,并与已知井信息进行比对,验证了其有效性。

简介：在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正.它们是C-波叠加速度VC2,垂直速度比和有效速度比γ0和γeff,以及各向异性参数χeff.我们将这四种参数作为C波叠加速度模型.C-波速度分析的目的就是确定这种叠加速度模型.C-波叠加速度模型VC2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得.然而错误的传播是C-波反射动校正反演中的严重问题.当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值.中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化.非双曲线分析通过中间偏移距的γ校正量可以产生VC2,若数据不含燥音,其误差小于1%.方法稳健,允许γ启始假定值的误差达20%.该方法也适用垂直非均匀各向异性介质.精度的提高使能够用4分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.