适用于层状介质的分频构形反演方法技术

本发明专利技术公开了一种适用于层状介质的分频构形反演方法，包括：利用克里金插值建立低频段阻抗体，利用有色反演将地震数据转化为中频阻抗体，利用基于地震波形相似度控制的构形建模以及随机模拟方法，即构形反演方法，建立高频阻抗体，将低频、中频、高频阻抗建模或反演的方法割裂开，再将三部分阻抗体合成一体，获得高分辨率的阻抗体，实现薄互层储层预测。通过本发明专利技术的技术方案，有效解决了高频段会由于井少造成插值过程中的随机性和不稳定性，提高了反演结果的稳定性和分辨率，提升了反演效果。

Inversion Method of Fractional Configuration Suitable for Layered Media

The invention discloses a frequency division configuration inversion method suitable for layered media, which includes: using Kriging interpolation to establish low frequency impedance body, using colored inversion to transform seismic data into intermediate frequency impedance body, using configuration modeling based on similarity control of seismic waveform and random simulation method, namely configuration inversion method. High-frequency impedance volume is established, and the modeling or inversion methods of low, medium and high frequency impedance are split, then the three parts of impedance volume are synthesized to obtain high-resolution impedance volume, so as to realize thin interbed reservoir prediction. The technical scheme of the invention effectively solves the randomness and instability in the interpolation process caused by fewer wells in the high frequency band, improves the stability and resolution of the inversion results, and improves the inversion effect.

【技术实现步骤摘要】
适用于层状介质的分频构形反演方法
本专利技术涉及地震学
，尤其涉及一种适用于层状介质的分频构形反演方法。
技术介绍
地质统计学反演原理，低频部分通过测井低频数据统计变差，获得低频变差函数，中频部分应用常规确定性反演方法(稀疏脉冲反演)得到波阻抗体，以了解储层的大致分布，并用于求取中频变差函数。高频部分从井点出发获得高频变差，通过三部分变差函数合成，通过随机模拟产生井间波阻抗，将全频带的波阻抗转换成反射系数并与确定性反演方法求得的子波进行褶积产生合成地震道，通过反复迭代直至合成地震道与原始地震道达到一定程度的匹配，反演获得多个等概率的波阻抗数据体。常规地质统计学存在问题，第一，统计合理的变差函数需对已知井的个数及分布要求严格，要求井数较多且分布均匀，对于横向变化快或者岩性变化复杂地区预测精度偏低。第二，对于高频随机模拟部分，常规地质统计需要随机模拟高达几百赫兹频率范围(地震高截频不到100Hz)。100Hz以上频率信息变差函数拟合难度大，效果不理想，反演结果随机性强，计算效率低。
技术实现思路
针对上述问题中的至少之一，本专利技术提供了一种适用于层状介质的分频构形反演方法，利用分频构建的思想，将低、中、中高、超高频四个部分阻抗建模或反演的方法割裂开，每一个部分都利用不同建模和反演方法进行计算，提高了反演结果的稳定，获得高分辨率的阻抗体，实现薄互层储层预测。利用地震波形相似程度插值代替变差函数统计的随机模拟，增强了高频信息中一部分结果的稳定性，反演结果平面上更具有地震相特征，连续性比常规地质统计学反演效果好。为实现上述目的，本专利技术提供了一种适用于层状介质的分频构形反演方法，包括：步骤S101，对真井的井曲线进行滤波得到低频阻抗，统计低频变差函数，并利用克里金插值法进行插值建模，获得低频阻抗体；步骤S102，利用有色反演法获得中频阻抗体，并将所述中频阻抗体与所述低频阻抗体合成为绝对波阻抗体；步骤S103，抽取所述真井的井旁道的井阻抗曲线并滤波获得高频阻抗以形成井集，在地震数据网格中抽取伪井的井旁道及序号生成对伪井进行反演的随机路径；步骤S104，按照所述随机路径的顺序并利用所述井集进行插值建模，获得待插点伪井的中高频信息；步骤S105，按照所述随机路径顺序对所述伪井进行超高频的随机模拟，以合成记录误差残差最小为约束建立贝叶斯判别规则，对采用MCMC算法进行寻优获得的多个等概率结果进行平均以获得超高频信息；步骤S106，以所述随机路径的顺序，将所有所述伪井的中高频信息和超高频信息合成后加入所述井集；步骤S107，以步骤S104和步骤S105方法对所述地震数据网格中的剩余道进行反演；步骤S108，以所述绝对波阻抗体的阻抗信息和高频阻抗的阻抗信息分别与真井对应频段匹配并组合为全频带阻抗，对所述地震数据网格进行逐道迭代形成高分辨率反演体；其中，所述低频频率范围为0≤f＜10Hz，所述中频为10≤f＜80Hz，所述高频包括所述中高频和所述超高频，所述中高频为80≤f＜200Hz，所述超高频为200≤f＜500Hz。在上述技术方案中，优选地，所述步骤S104中，利用待插点伪井的井旁道和所述真井的井旁道相似度和距离为权重，利用所述井集进行中高频的插值建模，获得待插点伪井的中高频信息。在上述技术方案中，优选地，所述步骤S105中按照所述随机路径顺序对所述伪井进行超高频的随机模拟具体包括：按所述随机路径顺序依次模拟，通过概率分析统计测井数据中超高频的均值和方差对超高频的随机模拟进行指导，同时利用地震平面变差分析指导高频的平面延伸的模拟。在上述技术方案中，优选地，所述步骤S106中，以所述随机路径的顺序将先生成的伪井的中高频信息和超高频信息合成并加入到所述井集中，将所述随机路径中所有路径点迭代反演完成，以完成所述井集的扩充。与现有技术相比，本专利技术的有益效果为：利用分频构建的思想，将低、中、中高、超高频四个部分阻抗建模或反演的方法割裂开，每一个部分都利用不同建模和反演方法进行计算，提高了反演结果的稳定，获得高分辨率的阻抗体，实现薄互层储层预测。利用地震波形相似程度插值代替变差函数统计的随机模拟，增强了高频信息中一部分结果的稳定性，反演结果平面上更具有地震相特征，连续性比常规地质统计学反演效果好。附图说明图1为本专利技术一种实施例公开的适用于层状介质的分频构形反演方法的流程图；图2为本专利技术一种实施例公开的适用于层状介质的分频构形反演方法的流程示意框图；图3为本专利技术一种实施例公开的分频构形反演的频谱构建形式示意图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本专利技术保护的范围。下面结合附图对本专利技术做进一步的详细描述：如图1至图3所示，根据本专利技术提供的一种适用于层状介质的分频构形反演方法，包括：步骤S101，对真井的井曲线进行带通滤波得到低频阻抗(0～10Hz)，统计低频变差函数，并利用克里金插值法进行插值建模，获得低频阻抗体；步骤S102，确定性反演部分，利用有色反演法获得中频阻抗体(10～80Hz)，并将中频阻抗体与低频阻抗体合成为绝对波阻抗体(频带范围0Hz～80Hz)；步骤S103，数据准备，抽取真井的井旁道的井阻抗曲线并滤波获得高频阻抗(频带范围80～500Hz)以形成井集，在地震数据网格(50*50)中抽取伪井的井旁道及序号，生成对伪井进行反演的随机路径；步骤S104，伪井构形插值，即中高频的插值建模，实现方法为：利用待插点伪井的井旁道和真井的井旁道相似度和距离为权重，利用井集进行中高频(频带范围80～200Hz)的插值建模，获得待插点伪井的中高频信息；步骤S105，伪井超高频(频带范围200～500Hz)的随机模拟：按照随机路径顺序对伪井进行超高频的随机模拟，以合成记录误差残差最小为约束，建立贝叶斯判别规则，采用MCMC(MarkovChainMonteCarlo，马尔科夫链蒙特卡洛)算法进行寻优，获得多个等概率结果进行平均，获得超高频信息；步骤S106，将伪井的中高频信息和超高频信息合成后，加入井集，以随机路径的顺序，将先生成的伪井加入真井的井集，影响下一个路径点伪井的插值结果，所有路径点反演完成后，井集扩充完毕；步骤S107，以步骤S104和步骤S105方法对地震数据网格中的剩余道进行反演；步骤S108，以0～80Hz绝对波阻抗体的阻抗信息和80～500Hz的高频阻抗的阻抗信息分别与真井对应频段的方差和均值进行匹配，并组合为0～500Hz的全频带阻抗，对地震数据网格进行逐道迭代，形成高分辨率反演体；其中，低频频率范围为0≤f＜10Hz，中频为10≤f＜80Hz，高频包括中高频和超高频，中高频为80≤f＜200Hz，超高频为200≤f＜500Hz。在上述实施例中，其分频统计建模和分频反演的思想，区别于常规地质统计学中分频统计变差、合成一体反演的思想，具体体现为：在步骤S101中利用克里金插值建立低频段阻抗体，步骤S102中利用有色反演将地震数据转化为中频阻抗体，步本文档来自技高网...

【技术保护点】
1.一种适用于层状介质的分频构形反演方法，其特征在于，包括：步骤S101，对真井的井曲线进行滤波得到低频阻抗，统计低频变差函数，并利用克里金插值法进行插值建模，获得低频阻抗体；步骤S102，利用有色反演法获得中频阻抗体，并将所述中频阻抗体与所述低频阻抗体合成为绝对波阻抗体；步骤S103，抽取所述真井的井旁道的井阻抗曲线并滤波获得高频阻抗以形成井集，在地震数据网格中抽取伪井的井旁道及序号生成对伪井进行反演的随机路径；步骤S104，按照所述随机路径的顺序并利用所述井集进行插值建模，获得待插点伪井的中高频信息；步骤S105，按照所述随机路径顺序对所述伪井进行超高频的随机模拟，以合成记录误差残差最小为约束建立贝叶斯判别规则，对采用MCMC算法进行寻优获得的多个等概率结果进行平均以获得超高频信息；步骤S106，以所述随机路径的顺序，将所有所述伪井的中高频信息和超高频信息合成后加入所述井集；步骤S107，以步骤S104和步骤S105方法对所述地震数据网格中的剩余道进行反演；步骤S108，以所述绝对波阻抗体的阻抗信息和高频阻抗的阻抗信息分别与真井对应频段匹配并组合为全频带阻抗，对所述地震数据网格进行逐道迭代形成高分辨率反演体；其中，所述低频频率范围为0≤f＜10Hz，所述中频为10≤f＜80Hz，所述高频包括所述中高频和所述超高频，所述中高频为80≤f＜200Hz，所述超高频为200≤f＜500Hz。...

【技术特征摘要】
1.一种适用于层状介质的分频构形反演方法，其特征在于，包括：步骤S101，对真井的井曲线进行滤波得到低频阻抗，统计低频变差函数，并利用克里金插值法进行插值建模，获得低频阻抗体；步骤S102，利用有色反演法获得中频阻抗体，并将所述中频阻抗体与所述低频阻抗体合成为绝对波阻抗体；步骤S103，抽取所述真井的井旁道的井阻抗曲线并滤波获得高频阻抗以形成井集，在地震数据网格中抽取伪井的井旁道及序号生成对伪井进行反演的随机路径；步骤S104，按照所述随机路径的顺序并利用所述井集进行插值建模，获得待插点伪井的中高频信息；步骤S105，按照所述随机路径顺序对所述伪井进行超高频的随机模拟，以合成记录误差残差最小为约束建立贝叶斯判别规则，对采用MCMC算法进行寻优获得的多个等概率结果进行平均以获得超高频信息；步骤S106，以所述随机路径的顺序，将所有所述伪井的中高频信息和超高频信息合成后加入所述井集；步骤S107，以步骤S104和步骤S105方法对所述地震数据网格中的剩余道进行反演；步骤S108，以所述绝对波阻抗体的阻抗信息和高频阻抗的阻抗信息分别与真井对应频段匹配并组合为全频带阻抗，对...

【专利技术属性】
技术研发人员：刘力辉，
申请(专利权)人：北京诺克斯达石油科技有限公司，
类型：发明
国别省市：北京,11