【技术实现步骤摘要】
本专利技术涉及地球物理勘探和开发，尤其涉及一种基于速度低频趋势建模的地震重构方法。
技术介绍
1、由于深洼带储层速度和密度受到压实作用的影响，随着埋深增大、压实作用增强，直接影响到三类储层的速度和密度变化特征，也即影响到波阻抗变化特征。受压实影响，平面上不同沉积相砂泥岩波阻抗存在差异性，测井曲线垂直向上砂泥岩波阻抗差随着埋深的增大逐渐减小，使得同一类储层不仅导致测井岩性波阻抗交叉重叠严重，而且波阻抗差异性减小产生弱反射降低了地震识别能力，压实作用引起的同一类储层波阻抗差异变化，将直接影响到地震记录特征的变化，使得从现有地震记录中提取三类有效储层对应的地震属性变化的规律性不强，难以优选出合适的地震属性来描述三类储层。
2、常用的地层去压实校正方法有压实率恢复法、压实系数恢复法、去压实恢复法、地层原始厚度恢复图版法、层序地层分析法、镜质组最大反射率恢复法、地球物理法等。压实系数恢复法是最常用的方法，其中压实系数是反演盆地埋藏过程的重要参数，考虑压实作用引起的同一类储层波阻抗差异变化，将直接影响到地震记录特征的变化，也是准确获取正常压实段泥岩压实系数的前提，波阻抗的低频(时间域)或低波数(深度域)分量随深度变化代表了地层波阻抗的压实趋势。为此，通过提取波阻抗的低频或低波数分量，可以求得时间域或深度域地层波阻抗压实系数。
技术实现思路
1、鉴于上述问题，提出了本专利技术以便提供克服上述问题或者至少部分地解决上述问题的一种基于速度低频趋势建模的地震重构方法。
2、根据
3、步骤s1：将测井资料转化为时间域，对测井波阻抗进行分频处理，选取与趋势相符的低频分量为测井低频分量数据体；
4、步骤s2：基于地震速度数据体及测井密度与速度的拟合关系计算地震波阻抗低频分量，利用所述测井低频分量数据体对所述地震波阻抗低频分量进行校正，得到校正后的三维空间地震波阻抗低频分量数据体；
5、步骤s3：根据校正后的三维空间地震波阻抗低频分量数据体确定目的层段对应的压实趋势线，计算目的层去压实后的波阻抗值，计算波阻抗去压实系数；
6、步骤s4：建立去压实前后的波阻抗关系、反射系数数学关系，建立多道地震记录与地震子波和多道波阻抗差分的数学关系；
7、步骤s5：在波阻抗低频分量趋势约束情形下，建立去压实后的波阻抗反演目标函数，采用非线性反演方法得到压实校正后的地震记录。
8、可选的，所述步骤s1：将测井资料转化为时间域，对测井波阻抗进行分频处理，选取与趋势相符的低频分量为测井低频分量数据体具体包括：
9、测井数据计算时间域测井波阻抗数据体，采用广义s变换时频谱分析，确定目的层井段代表储层压实趋势的波阻抗低频分量范围；
10、采用广义s反变换对测井波阻抗低频分量进行重构，获得测井时间域波阻抗低频分量。
11、可选的，所述步骤s2：基于地震速度数据体及测井密度与速度的拟合关系计算地震波阻抗低频分量，利用测井低频分量对地震波阻抗分量进行校正，得到校正后的三维空间地震波阻抗低频分量数据体具体包括：
12、基于地震偏移速度数据体、解释的主要目的层数据，采用层控多约束层速度反演mcvi方法计算地震层速度；
13、进行目的层测井速度和测井密度交汇图分析，获得密度与速度的拟合公式，并应用于地震层速度换算出地层密度，计算地震波阻抗体；
14、利用广义s变换获取低频波阻抗体；
15、利用测井波阻抗低频分量进行校正，得到校正后的三维空间地震波阻抗低频分量数据体。
16、可选的，所述步骤s3：根据校正后的三维空间地震波阻抗低频分量数据体确定目的层段对应的压实趋势线，计算目的层去压实后的波阻抗值，计算波阻抗去压实系数具体包括：
17、根据低频波阻抗值是因为压实作用随深度变化特征，确定目的层段对应的压实趋势线，计算目的层去压实后的波阻抗值；
18、利用去压实前后的波阻抗值计算波阻抗压实系数，得到目的层段测井和三维波阻抗压实系数随深度变化曲线。
19、可选的，所述步骤s3中计算井旁地震波阻抗压实系数与测井波阻抗去压实系数的差异，结合解释层位，采用克里金插值校正技术，建立层控地震波阻抗压实系数误差数据体；
20、对地震波阻抗压实系数进行校正，得到校正后的三维空间地震波阻抗压实系数数据体。
21、可选的，所述步骤s4：建立去压实前后的波阻抗关系、反射系数数学关系，建立多道地震记录与地震子波和多道波阻抗差分的数学关系具体包括：
22、根据反射系数与波阻抗之间的非线性关系和差分算子，建立反射系数、多道波阻抗差分和地震记录的关系。
23、可选的，所述步骤s5：在波阻抗低频分量趋势约束情形下，建立去压实后的波阻抗反演目标函数，采用非线性反演方法得到压实校正后的地震记录具体包括：
24、建立低频约束目标函数，采用非线性反演方法得到压实校正后的波阻抗，得到去压实校正后的反射系数；
25、利用井旁资料提取的地震子波，通过褶积得到去压实校正后的地震记录。
26、可选的，所述步骤s5中的非线性反演方法中采用了谱投影梯度算法算法求解到压实校正后的波阻抗。
27、可选的，所述步骤s5中建立去压实后的波阻抗反演目标函数具体包括：多道波阻抗反演问题表示为
28、
29、其中，w为子波时间序列组成的行向量，xlow为波阻抗低频分量，从测井波阻抗或地震速度谱获得；l为低通滤波算子；y为已知的多道地震记录，为所反演的反射系数时间序列具有稀疏性和横向光滑性；为合成地震记录与实际地震记录吻合约束；xlow＝lx为波阻抗低频分量趋势约束；
30、可选的，所述步骤s5中建立去压实后的波阻抗反演目标函数还包括：去压实后的波阻抗反演问题表示为
31、
32、式中，c为压实系数组成的矩阵；
33、反演求得去压实后的波阻抗半对数xc；
34、利用差分算子计算去压实后的多道反射系数rc；
35、利用井旁地震纪录提取的子波计算去压实后的地震记录yc。
36、本专利技术提供的一种基于速度低频趋势建模的地震重构方法包括：步骤s1：将测井资料转化为时间域，对测井波阻抗进行分频处理，选取与趋势相符的低频分量为测井低频分量数据体；步骤s2：基于地震速度数据体及测井密度与速度的拟合关系计算地震波阻抗低频分量，利用所述测井低频分量数据体对所述地震波阻抗低频分量进行校正，得到校正后的三维空间地震波阻抗低频分量数据体；步骤s3：根据校正后的三维空间地震波阻抗低频分量数据体确定目的层段对应的压实趋势线，计算目的层去压实后的波阻抗值，计算波阻抗去压实系数；步骤s4：建立去压实前后的波阻抗关系、反射系数数学关系，建立多道地震记录与地震子波和多道波阻抗差分的数学关系；步骤s5：在本文档来自技高网...
【技术保护点】
1.一种基于速度低频趋势建模的地震重构方法，其特征在于，所述重构方法包括：
2.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S1：将测井资料转化为时间域，对测井波阻抗进行分频处理，选取与趋势相符的低频分量为测井低频分量数据体具体包括：
3.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S2：基于地震速度数据体及测井密度与速度的拟合关系计算地震波阻抗低频分量，利用测井低频分量对地震波阻抗分量进行校正，得到校正后的三维空间地震波阻抗低频分量数据体具体包括：
4.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S3：根据校正后的三维空间地震波阻抗低频分量数据体确定目的层段对应的压实趋势线，计算目的层去压实后的波阻抗值，计算波阻抗去压实系数具体包括：
5.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S3中计算井旁地震波阻抗压实系数与测井波阻抗去压实系数的差异，结合解释层位，采用克里金插值校正技术，建立层控
6.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S4：建立去压实前后的波阻抗关系、反射系数数学关系，建立多道地震记录与地震子波和多道波阻抗差分的数学关系具体包括：
7.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S5：在波阻抗低频分量趋势约束情形下，建立去压实后的波阻抗反演目标函数，采用非线性反演方法得到压实校正后的地震记录具体包括：
8.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S5中的非线性反演方法中采用了谱投影梯度算法算法求解到压实校正后的波阻抗。
9.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S5中建立去压实后的波阻抗反演目标函数具体包括：
10.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤S5中建立去压实后的波阻抗反演目标函数还包括：
...【技术特征摘要】
1.一种基于速度低频趋势建模的地震重构方法，其特征在于，所述重构方法包括：
2.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤s1：将测井资料转化为时间域，对测井波阻抗进行分频处理，选取与趋势相符的低频分量为测井低频分量数据体具体包括：
3.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤s2：基于地震速度数据体及测井密度与速度的拟合关系计算地震波阻抗低频分量，利用测井低频分量对地震波阻抗分量进行校正，得到校正后的三维空间地震波阻抗低频分量数据体具体包括：
4.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤s3：根据校正后的三维空间地震波阻抗低频分量数据体确定目的层段对应的压实趋势线，计算目的层去压实后的波阻抗值，计算波阻抗去压实系数具体包括：
5.根据权利要求1所述的一种基于速度低频趋势建模的地震重构方法，其特征在于，所述步骤s3中计算井旁地震波阻抗压实系数与测井波阻抗去压实系数的差异，结合解释层位，...
【专利技术属性】
技术研发人员：时秀朋，谭绍泉，陈学国，尹克敏，王千军，张建华，钱焕菊，李守济，魏敏，冯国志，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。