【技术实现步骤摘要】
本专利技术涉及一种潜山风化带底界面解释方法、装置、介质及设备，属于地球物理油气勘探领域。
技术介绍
1、随着油气勘探开发程度的不断深入，潜山油气藏已经成为国内外油气勘探的一个重要领域。潜山油气藏的研究重点之一是结合储层发育规律，解释潜山储层的顶、底界面，并计算潜山储层厚度。潜山储层自上而下分为风化带和内幕带，通常潜山顶部除构造运动形成构造缝外，在出露地表期间遭受风化淋滤作用改造会形成大量风化缝、收缩缝、溶蚀孔缝等储集空间类型，而越靠近潜山顶部，风化作用越强，储层越发育，油气勘探开发中一般将受构造和风化淋滤作用双重影响的这部分储集层段称为风化带，将风化带以下，基本不受风化作用改造影响、以构造缝为主的储集空间成为内幕带。潜山风化带是潜山储层的最有利发育相带之一，准确解释潜山风化带的顶、底界面，并计算潜山风化带厚度具有重要的研究意义。
2、通常潜山风化带的顶界面在地震剖面上具有较强的反射特征，比较容易解释，而潜山风化带底界面的地震响应往往较弱，并且潜山层段地震资料品质往往较差，这给潜山风化带界底面的解释工作带来较大困难，传统的人工解释耗时耗力，且解释准确度与精度有限，针对潜山风化带底界面解释难题开展方法研究具有重要的理论与应用意义。
技术实现思路
1、针对上述技术问题，本专利技术提供一种基于厚度趋势约束的潜山风化带底界面解释方法，该方法在优化地震资料基础的同时，利用厚度趋势作为层位自动追踪走廊，确保了解释结果的合理性和闭合性，克服了传统的人工解释耗时、耗力、解释准确度与精度有
2、为实现上述目的，本专利技术采取以下技术方案：
3、一种潜山风化带底界面解释方法，包括如下步骤：
4、对地震资料进行滤波优化，压制高角度的断面波，并得到优化后的地震资料；
5、基于潜山古地貌恢复结果与已知井潜山风化带厚度信息计算潜山风化带厚度趋势；
6、基于优化后的地震资料，以已知井点为初始种子点，以潜山风化带厚度趋势为约束，以地震波形相似性为准则，进行种子点外推潜山风化带底界面自动解释。
7、所述的潜山风化带底界面解释方法，优选地，地震资料为一个三维数据体，三维数据体的振幅值为垂直方向的时间坐标t和两个水平方向的空间坐标x、y的函数f(x,y,t)，x∈[0,xm]、y∈[0,ym]、t∈[0,tm]，m表示当前坐标范围的最大值。
8、所述的潜山风化带底界面解释方法，优选地，利用f-k滤波方法对地震资料进行滤波优化。
9、所述的潜山风化带底界面解释方法，优选地，利用f-k滤波方法对地震资料进行滤波优化，压制高角度的断面波，并得到优化后的地震资料，具体步骤如下：
10、首先对三维数据体f(x,y,t)做f-k变换，得到如下结果：
11、
12、其中，kx∈[-kxh,kxh]、ky∈[-kyh,kyh]、ω∈[-ωh,ωh]，kxh、kyh、ωh的取值要求是包含g(kx,ky,ω)的有效信息，即当或者或者时，g1(kx,ky,ω)≈0；
13、压制高角度的断面波信息，求取g1(kx,ky,ω)，使得：
14、
15、其中，θmin为压制高角度的断面波信息的最小倾角；
16、对三维数据体g1(kx,ky,ω)进行f-k反变换，得到优化后的地震资料f1(x,y,t)：
17、
18、所述的潜山风化带底界面解释方法，优选地，基于潜山古地貌恢复结果与已知井潜山风化带厚度信息计算潜山风化带厚度趋势，具体步骤如下：
19、基于潜山古地貌恢复结果提取潜山古地貌相对高度数据b(x,y)，以b(x,y)为约束，以已知井的风化带厚度为硬数据，采用协克里金插值的方法，对已知的风化带厚度数据进行插值，得到潜山风化带厚度趋势q(x,y)，其中，x∈[0,xm]、y∈[0,ym]为三维体的空间坐标范围，为已知井点的空间坐标，下标i为已知井点坐标的索引，nw为已知井的个数。
20、所述的潜山风化带底界面解释方法，优选地，基于优化后的地震资料，以已知井点为初始种子点，以潜山风化带厚度趋势为约束，以地震波形相似性为准则，进行种子点外推潜山风化带底界面自动解释，具体步骤如下：
21、首先，利用已知井的风化带厚度构造层位追踪边界点集合i，公式如下：
22、
23、其中，与为已知井点的潜山风化带顶界与底界的时间位置；种子点外推潜山风化带底界面自动解释的目的就是以已知种子点集合i为基础逐渐外推，直到得到所有空间位置的潜山风化带底界面时间位置，即得到风化带底界集合：
24、
25、其中，ix与iy为三维地震数据体f1(x,y,t)在两个空间维方向的离散索引，nx与ny为两个空间维方向的采样点个数；
26、标记层位追踪边界点集合i中点的空间位置的层位解释状态为已解释，其余空间位置的层位解释状态为未解释，对于层位追踪边界点集合中的每一个点，设该点为(xc,yc,hdown(xc,yc))，其周围8个点的空间维坐标分别为(xc-1,yc-1)、(xc-1,yc)、(xc-1,yc+1)、(xc,yc-1)、(xc,yc+1)、(xc+1,yc-1)、(xc+1,yc)、(xc+1,yc+1)，定义这8个点中层位解释状态为未解释的点为待追踪边界点，为定义(xc,yc)为这8个待追踪边界点的父节点，计算每个待追踪边界点的最优潜山风化带底界面时间位置和外推参考值r；
27、对层位追踪边界点集合i中的所有点的待追踪边界点计算外推参考值r，选取r最小的待追踪边界点，作为生效边界点，并标记其层位解释状态为已解释，移除层位追踪边界点集合i中此生效边界点的父节点，并将次生效边界点加入层位追踪边界点集合i，重复以上过程，直到所有空间位置的层位解释状态均为已解释。
28、本专利技术第二方面提供一种潜山风化带底界面解释装置，包括：
29、第一处理单元，用于对地震资料进行滤波优化，压制高角度的断面波，并得到优化后的地震资料；
30、第二处理单元，用于潜山古地貌恢复结果与已知井潜山风化带厚度信息计算潜山风化带厚度趋势；
31、第三处理单元，用于基于优化后的地震资料，以已知井点为初始种子点，以潜山风化带厚度趋势为约束，以地震波形相似性为准则，进行种子点外推潜山风化带底界面自动解释。
32、本专利技术第三方面提供一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现上述任意一项所述潜山风化带底界面解释方法的步骤。
33、本专利技术第四方面提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任意一项所述潜山风化带底界面解释方法的步骤。
34、本专利技术由于采取以上技术方案，其具有以下优本文档来自技高网...
【技术保护点】
1.一种潜山风化带底界面解释方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的潜山风化带底界面解释方法，其特征在于，地震资料为一个三维数据体，三维数据体的振幅值为垂直方向的时间坐标t和两个水平方向的空间坐标x、y的函数f(x,y,t)，x∈[0,xm]、y∈[0,ym]、t∈[0,tm]，m表示当前坐标范围的最大值。
3.根据权利要求2所述的潜山风化带底界面解释方法，其特征在于，利用f-k滤波方法对地震资料进行滤波优化。
4.根据权利要求3所述的潜山风化带底界面解释方法，其特征在于，利用f-k滤波方法对地震资料进行滤波优化，压制高角度的断面波，并得到优化后的地震资料，具体步骤如下：
5.根据权利要求1所述的潜山风化带底界面解释方法，其特征在于，基于潜山古地貌恢复结果与已知井潜山风化带厚度信息计算潜山风化带厚度趋势，具体步骤如下：
6.根据权利要求1所述的潜山风化带底界面解释方法，其特征在于，基于优化后的地震资料，以已知井点为初始种子点，以潜山风化带厚度趋势为约束，以地震波形相似性为准则，进行种子点外推潜山风化带底
7.一种潜山风化带底界面解释装置，其特征在于，包括：
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-6任意一项所述潜山风化带底界面解释方法的步骤。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-6任意一项所述潜山风化带底界面解释方法的步骤。
...【技术特征摘要】
1.一种潜山风化带底界面解释方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的潜山风化带底界面解释方法，其特征在于，地震资料为一个三维数据体，三维数据体的振幅值为垂直方向的时间坐标t和两个水平方向的空间坐标x、y的函数f(x,y,t)，x∈[0,xm]、y∈[0,ym]、t∈[0,tm]，m表示当前坐标范围的最大值。
3.根据权利要求2所述的潜山风化带底界面解释方法，其特征在于，利用f-k滤波方法对地震资料进行滤波优化。
4.根据权利要求3所述的潜山风化带底界面解释方法，其特征在于，利用f-k滤波方法对地震资料进行滤波优化，压制高角度的断面波，并得到优化后的地震资料，具体步骤如下：
5.根据权利要求1所述的潜山风化带底界面解释方法，其特征在于，基于潜山古地貌恢复结果...
【专利技术属性】
技术研发人员：马淑芳，高云峰，周建楠，聂妍，张显文，杜昕，樊鹏军，肖大坤，
申请(专利权)人：中海石油中国有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。