CN114966856B 基于多频带地震资料的碳封存场址优 选方法、系统和设备 （中国科学院地质与地球物理研究所）

申请号202210919676.5所号基于多频带地震资料的碳封存场址优选方算稳定地层波阻抗数据体与波阻抗反演数据体2获取已知井位的原始地震资料和测井资料，进基于所述扩频模拟子波，构建目标层位的等时地层格架基于所述三维叠后地震资料和扩频模拟子波，进行多尺度分解获得地震资料多尺度基于所述扩频模拟子波和标准化测井数据，进行井基于所述测井资料、提频后的地震资料和时深转换根据所述断溶体储层解释模型，计算二氧化碳注入其中，s(o)表示傅氏变换后的地震记录频谱，o()表示傅氏变换后的地震子波频谱，e()表示傅氏变换后的反射系数频谱，w表示角频率；将所述地震记录褶积模型在频率域的表示转化为地震记其中，Ins()为地震记录线性系统频域表示，Ino()为地震子波线性系统频域表示，S(t)-i(t)+E(t)3通过低通滤波器将复赛谱中的子波复赛谱序列和反射设子波o(t)的最大相位分量为u(t)、最o(t)=u(t)·v(t)27(t)=i(t)+j(t)+i(-t)+(-t)(t)(-t)为地震子波最小相位分量v(t)所对应的最大相位函数的复赛谱，i(-t)为地震子波最大相位分量u(t)所对应的基于所述振幅谱的复赛谱确定一组具有相同振幅谱的n(r,y)-u(r,y)设置变量l为分解的Bimf层数，l的初h(x,y)为：h(x,)=ri-1(x,y)计算关于h(x,y)的局部最大值矩阵和局通过设置观测矩阵获取观测矩阵中的最小值，进而获4求得NUX为最大值包络面，将过度公式B中的局部最大值矩阵Max(r,y)替换为局部最Bimfi(x,y)=h(x,)-m-1(x,y)ni(x,)=ri-1(x,y)-Bimfi(x,)基于每个已知井位的测井资料中的声波时差曲线和密度曲线做乘积运算获取波阻抗以目标层段地震主频为依据构建雷克子波，将雷克子波与反射系数曲线褶积计算后，将每个钻井井位井眼处的标志层深度数据与标志层的震记录与井旁地震道扩频模拟子波的相关性，当波形相关性高于预设的第一相关阈值时，将所述地震资料多尺度Bimf分量以尺度由大到小的顺序，5每添加一个Bimf分量获得更新后的有效地震资料计算一次更新后的合成地震记录与将所述第二时深转化关系作为最终测得的时基于所述扩频模拟子波，计算待判别地震道波形以所述稳定的目标函数作为初始模型的输入，通过马尔科夫链蒙特卡罗方法MCMC和6将采样区域进行纵向和横向平移，遍历计算所有采样区对所述三维方差属性体进行切片，获取方差属性数据在平面上的分设未知地下波阻抗模型参数m为维空间向量m=fm,…,mg}；通过非线性函数核G建立未知地下波阻抗模型参数和k维数通过一次地震采集过程获得叠前深度域偏移地震剖面a'，7设等价的深度域模型概率分布服从高斯分布，有8基于虚拟井井位平面坐标与地震数据平面坐标的相关关系频波阻抗反演数据的一一对应关系，进而将虚拟井井旁低频波阻抗反演数据对虚拟井赋基于波阻抗反演数据体提取虚拟井井位波阻抗数据，根据地层格架限定的时窗范围资料获取模块，配置为获取已知井位的原始地震资预处理模块，配置为基于所述测井资料，进行剔等时地层格架模型获取模块，配置为基于所述扩Bimf分量获取模块，配置为基于所述三维叠后地震资料和扩频模拟波阻抗反演数据体获取模块，配置为基于所述未定地层波阻抗数据体获取模块，配置为基于波阻抗异常数据体获取模块，配置为计算所述稳断溶体波阻抗数据体获取模块，配置为通过去除三维方差属性体中低于平均值的区模型解释模块，配置为比对所述井旁地质解释结果和断根据所述断溶体储层解释模型，计算二氧化碳注入述处理器执行以实现权利要求1-7任一项所述的基于多频带地震资料的碳封存场址优选方9令，所述计算机指令用于被所述计算机执行以实现权利要求1-7任一项所述的基于多频带[0020]比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、[0024]n(r,y)-u(r,y)x为矩阵行坐标，y为矩阵列坐标，v为叠后三维地震资料数据的值，ro为矩阵h(x,y)为：[0027]h(x,)=ri-1(x,y)[0029]在变量矩阵h(x,y)中以预设尺寸的观测矩阵进行观测，提取观测矩阵中的元素;[0046]ni(x,)=ri-1(x,y)-Bimfi(x,)Ti-1表示第l-1层的余量，第l-1层的余量作为第l层计算[0052]基于每个已知井位的测井资料中的声波时差曲线和密度曲线做乘积运算获取波[0056]其中，Ti,表示声波测井标志层深[0058]每添加一个Bimf分量获得更新后的地震资料计算一次更新后的合成地震记录与关性出现下降趋势时，取第二相关性峰值时的更新后的合成地震记录与更新后的地震资[0080]图1是本发明实施例中基于多频带地震资料的碳封存场址优选方法的流程示意区面积约27km2的三维叠后地震资料数据，信号记录双程旅行时为4s，采样点时间间隔为[0107]其中，s(o)表示傅氏变换后的地震记录频谱，o()表示傅氏变换后的地震子波频[0110]其中，Ins()为地震记录线性系统频域表示，lno()为地震子波线性系统频域表[0114]通过低通滤波器将复赛谱中的子波复赛谱序列和反射系数复H(f)和A(f)为待求的关于频率f的多项式，用于拟合地[0122]27(t)=i(t)+j(t)+i(-t)+(-t)(t)(-t)为地震子波最小相位分量V(t)所对应的最小相位函数的复赛谱，i(-t)为地震子波最大相位分量QUOTE的结果，高分辨率意味着需要高频带的地震数据，但是现在有效频带的分布范围为20-。[0131]n(r,y)-u(r,y)x为矩阵行坐标，y为矩阵列坐标，v为叠后三维地震资料数据的值，ro为矩阵h(x,y)为：[0134]h(x,)=ri-1(x,y)[0136]在变量矩阵h(x,y)中以预设尺寸的观测矩阵进行观测，提取观测矩阵中的元素[0147]求得NUX为最大值包络面，将过度公式B中的局部最大值矩阵Max(r,y)替换为局[0151]Bimfi(x,y)=h(x,)-m-1(x,y)[0153]ni(x,)=ri-1(x,y)-Bimfi(x,)Ti-1表示第l-1层的余量，第l-1层的余量作为第l层BimfBimf分量则表征了信号的内部结构特征和[0164]基于每个已知井位的测井资料中的声波时差曲线和密度曲线做乘积运算获取波[0168]其中，Ti,表示声波测井标志层深度对应的地震资料的双程旅行时；Tia为声波时[0170]每添加一个Bimf分量获得更新后的有效地震资料计算一次更新后的合成地震记[0187]设三维叠后地震资料中的各采样点数据为ss，其中p表示地震测网线号，a代[0194]设未知地下波阻抗模型参数m为维空间向量m=fm,…,[0204]对所述正演模型进行泰勒展开并略去二次以上的高阶项，获得预测数据简略表尝试对待反演模型参数进行量化约束以控制解空间大小，以期得到稳定精确的反演结果，[0211]通过一次地震采集过程获[0212]假设反演深度域模型参数有先验概率分布P(M=m)，根据贝叶斯公式有概率分布解[0218]n为噪音且服从高斯分布均值为0，cn为一等值的对角矩阵且对角线元素为数据[0220]对于反演深度域模型参数先验概率分布P(M-m)，设置初始模型，且有[0226]加入随机噪音的概率分布变形中的分母P(D=d)依赖于数据采集和处理，设为常[0236]随机逆反演的基本公式的迭代公式为：m*=mf+围针对所有虚拟井波阻抗数据的插值计算，设定计算区域为地质格架内工区的全部范围，[0243]比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、[0246]比对测井解释结果与断溶体结构-特征值模型的特征值能量数据，划定特征值大其可以同时(并行)执行或以颠倒的次序执行，这些简单的变化都在本发明的保护范围之[0249]本发明第二实施例的基于多频带地震资料的碳