1.一种用于砂泥岩薄互层砂体厚度预测的方法,包括:
获取90°相位地震数据,并从所述90°相位地震数据中提取90°相位子波;
基于所述90°相位子波和单砂体地质模型绘制目标砂体地震反射的过零点时间间隔与目标砂体时间厚度的关系模板;
利用目标砂体的空间展布形态确定目标砂体地震反射的过零点时间间隔;
基于目标砂体地震反射的过零点时间间隔和所述关系模板计算目标砂体的厚度。
2.根据权利要求1所述的方法,其特征在于,建立所述关系模板的步骤进一步包括:
利用所述90°相位子波与所述单砂体地质模型的反射系数褶积合成目标砂体的地震记录;
拾取每个地震记录的两个过零点时间,并计算每个地震记录的过零点时间间隔;
以所述单砂体地质模型的砂体时间厚度为横坐标、对应的过零点时间间隔为纵坐标,绘制过零点时间间隔与砂体时间厚度关系曲线。
3.根据权利要求2所述的方法,其特征在于,所述单砂体地质模型为楔形模型,其最小时间厚度为零,最大时间厚度不小于地震子波视周期的一半,顶面反射系数和底面反射系数符号相反、幅度相等。
4.根据权利要求1-3中任一项所述的方法,其特征在于,确定目标砂体地震反射的过零点时间间隔进一步包括:
利用地震切片解释技术确定目标砂体的空间展布形态;
以目标砂体所在的地震层位为时窗中心,以地震子波的周期为时窗长度,确定地震切片分析时窗;
在所述分析时窗内,以采样间隔为增量,从浅到深提取水平地震切片;
基于目标砂体的空间展布形态,在所述地震切片中寻找不包含目标砂体反射影像的两个地震切片,所述两个地震切片所对应的时间为目标砂体的两个过零点时间,这两个过零点时间之差即为目标砂体地震反射的过零点时间间隔。
5.根据权利要求4所述的方法,其特征在于,基于目标砂体地震反射的过 零点时间间隔和所述关系模板计算目标砂体的厚度的步骤进一步包括:
基于目标砂体地震反射的过零点时间间隔和所述关系模板确定砂体时间厚度;
将所述砂体时间厚度的一半与砂体速度相乘以得到目标砂体的厚度。
6.一种用于砂泥岩薄互层砂体厚度预测的系统,包括:
子波提取模块,获取90°相位地震数据,并从所述90°相位地震数据中提取90°相位子波;
关系模板建立模块,基于所述90°相位子波和单砂体地质模型绘制目标砂体地震反射的过零点时间间隔与目标砂体时间厚度的关系模板;
过零点时间间隔计算模块,利用目标砂体的空间展布形态确定目标砂体地震反射的过零点时间间隔;
目标砂体厚度计算模块,基于目标砂体地震反射的过零点时间间隔和所述关系模板计算目标砂体的厚度。
7.根据权利要求6所述的系统,其特征在于,所述关系模板建立模块包括:
单砂体地震记录合成单元,利用所述90°相位子波与所述单砂体地质模型的反射系数褶积合成目标砂体的地震记录;
地震记录过零点时间间隔计算单元,拾取每个地震记录的两个过零点时间,并计算每个地震记录的过零点时间间隔;
关系模板绘制单元,以所述单砂体地质模型的砂体时间厚度为横坐标、对应的过零点时间间隔为纵坐标,绘制过零点时间间隔与砂体时间厚度关系曲线。
8.根据权利要求7所述的系统,其特征在于,所述单砂体地震记录合成单元中采用的单砂体地质模型为楔形模型,其最小时间厚度为零,最大时间厚度不小于地震子波视周期的一半,顶面反射系数和底面反射系数符号相反、幅度相等。
9.根据权利要求6-8中任一项所述的系统,其特征在于,所述过零点时间间隔计算模块包括:
目标砂体空间展布形态确定单元,利用地震切片解释技术确定目标砂体的空间展布形态;
地震切片分析时窗计算单元,以目标砂体所在的地震层位为时窗中心,以地震子波的周期为时窗长度,确定地震切片分析时窗;
地震切片提取单元,在所述分析时窗内,以采样间隔为增量,从浅到深提取 水平地震切片;
目标砂体地震反射过零点时间间隔获取模块,基于目标砂体的空间展布形态,在所述地震切片中寻找不包含目标砂体反射影像的两个地震切片,所述两个地震切片所对应的时间为目标砂体的两个过零点时间,这两个过零点时间之差即为目标砂体地震反射的过零点时间间隔。
10.根据权利要求9所述的系统,其特征在于,所述目标砂体厚度计算模块包括:
砂体时间厚度计算单元,基于目标砂体地震反射的过零点时间间隔和所述关系模板确定砂体时间厚度;
目标砂体厚度计算单元,将所述砂体时间厚度的一半与砂体速度相乘以得到目标砂体的厚度。