一种砂岩型铀矿地震综合预测方法技术

技术编号:17031404 阅读:22 留言:0更新日期:2018-01-13 18:32
本发明专利技术属于铀矿地震综合预测方法技术领域,具体涉及一种砂岩型铀矿地震综合预测方法。本发明专利技术使用地震波阻抗反演法反算地下砂体分布情况,得到与铀成矿有直接关系的目标地层含砂率信息,同时利用钻孔数据统计的转换关系或经验转换关系将目标层波阻抗信息转换为与成矿作用有关的孔隙度信息,最后利用三维地震数据体提取与铀成矿环境相关的敏感地震属性信息,通过这三类信息综合分析研究预测工区的铀成矿潜力,能够达到快速、有效地预测砂岩型铀矿目标储层发育范围的目的。

【技术实现步骤摘要】
一种砂岩型铀矿地震综合预测方法
本专利技术属于铀矿地震综合预测方法
,具体涉及一种砂岩型铀矿地震综合预测方法。
技术介绍
砂岩型铀矿预测评价工作是砂岩型铀矿勘查中一项关键任务。常规地质预测技术主要利用地质调查、钻探资料进行预测,缺点在于预测的成本较高、周期较长;使用地球化学勘查、放射性物探、电磁法等地球物理方法进行预测评价,其结果往往达不到精度要求、效果不佳;而地震勘探技术对于砂岩型铀矿成矿环境有较好的探测精度,使用地震技术对砂岩型铀矿目标储层进行综合预测具有较好的前景。
技术实现思路
本专利技术要解决的技术问题是提供一种砂岩型铀矿地震综合预测方法,该方法能够快速、有效地预测砂岩型铀矿目标储层的平面展布范围。为了实现这一目的,本专利技术采取的技术方案是:一种砂岩型铀矿地震综合预测方法,该方法包括以下步骤:(1)在研究区,采集一套砂岩型铀矿三维地震纯波数据;采集野外三维地震数据,对野外三维地震数据依次进行处理以得到三维地震纯波数据;(2)使用基于模型反演法对步骤(1)的三维地震纯波数据进行反演计算,得到三维波阻抗数据体,进而计算三维岩性数据体;(3)求取目标层位砂体的含砂率分布数据,并绘制成等值线图;利用步骤(2)得到的三维岩性数据体,利用geoview软件的tracemath模块计算目标地层的含砂率数据,并绘制平面等值线图;(4)求取孔隙度数据与波阻抗数据的关系①在研究区测井曲线中有孔隙度数据、声波数据、密度数据三类测井数据的情况下:使用excel软件或geoview软件的crossplot模块进行孔隙度和波阻抗的交会分析;分析时,设置x横轴为波阻抗数据,波阻抗数据=声波和密度数据的乘积;y纵轴为孔隙度数据;然后通过excel软件或geoview软件的crossplot模块的线性拟合工具拟合出波阻抗转孔隙度的转换式y=ax+b,y为孔隙度数据,x为波阻抗数据,a,b为需要拟合的参数;②在研究区无孔隙度测井数据的情况下:使用下述的公式(1),将波阻抗转换为孔隙度数据:其中,AC固体骨架表示岩石固体骨架的声波时差,AC流体表示岩石孔隙间流体的声波时差,IMP为波阻抗数据,POR为孔隙度数据;(5)利用步骤(4)求取的波阻抗转孔隙度的转换关系,将波阻抗数据转换为孔隙度数据,得到目标地层孔隙度数据分布,并绘制等值线图;(6)利用三维地震纯波数据,提取地震均方根振幅、瞬时相位、弧长3种地震属性,对此三种地震属性进行聚类分析,基于钻孔资料建立约束条件,回归拟合出一种新的地震属性组合;(7)交会分析上述目标地层含砂率数据、孔隙度分布数据以及地震属性组合分布数据,综合预测砂岩型铀矿成矿有利地段;综合预测砂岩型铀矿成矿有利地段是指对目标地层含砂率分布图中数值大于X的区域、孔隙度值大于Y的区域、地震属性组合值大于Z的区域进行交会分析研究,据此再进行成矿有利地段的综合预测。进一步的,如上所述的一种砂岩型铀矿地震综合预测方法,步骤(1)中,通过地震仪采集野外三维地震数据,对野外三维地震数据进行处理的过程是依次进行静校正、去噪、振幅补偿、反褶积、动校叠加、偏移处理。进一步的,如上所述的一种砂岩型铀矿地震综合预测方法,步骤(2)中,具体步骤如下:①初始模型的建立:首先收集研究区的钻井数据,对其中的声波曲线和密度曲线进行平滑处理和标准化处理;②利用geoview软件的STRATA模块建立反演初始模型;③波阻抗反演计算:对步骤(1)得到的三维地震纯波数据进行反演计算,得到三维波阻抗数据体;④确定波阻抗转岩性的门槛值:每个研究区的波阻抗转岩性的门槛值不尽相同,门槛值的确定基于工区岩石物性参数分析;⑤岩性数据计算:基于步骤(2)④得到的波阻抗转岩性的门槛值,将步骤(2)③得到的三维波阻抗数据体转换为三维岩性数据体,三维岩性数据体的形式为LITH(x,y,t),其中x表示三维地震数据的Inline号,y表示三维地震数据的Xline号,t表示时间;LITH(x,y,t)数据体的取值为0或1,0表示泥岩,1表示砂岩。进一步的,如上所述的一种砂岩型铀矿地震综合预测方法,步骤(2)①中,平滑处理时采用3点或5点平滑处理;密度数据的3点或5点滑动平均处理分别如以下公式所示:●密度3点滑动平均公式:dden(i)=(dden(i-1)+dden(i)+dden(i+1))/3其中,di代表某个采样点的密度值,di-1为该采样点的前一个采样点的密度值,di+1为该采样点的后一个采样点的密度值;●密度5点滑动平均公式:dden(i)=(dden(i-2)+dden(i-1)+dden(i)+dden(i+1)+dden(i+2))/5其中,di代表某采样点的密度值,di-2为该采样点的前两个采样点的密度值,di-1为该采样点的前一个采样点的密度值,di+1为该采样点的后一个采样点的密度值,di+2为该采样点的后两个采样点的密度值;声波数据的3点或5点滑动平均处理分别如以下公式所示:●声波3点滑动平均公式:dson(i)=(dson(i-1)+dson(i)+dson(i+1))/3其中,di代表某个采样点的声波值,di-1为该采样点的前一个采样点的声波值,di+1为该采样点的后一个采样点的声波值;●声波5点滑动平均公式:dson(i)=(dson(i-2)+dson(i-1)+dson(i)+dson(i+1)+dson(i+2))/5其中,di代表某采样点的声波值,di-2为该采样点的前两个采样点的声波值,di-1为该采样点的前一个采样点的声波值,di+1为该采样点的后一个采样点的声波值,di+2为该采样点的后两个采样点的声波值;步骤(2)①中,标准化处理时使用geoview软件的loggingnomalize模块进行处理,将用于建立初始模型的井的声波和密度数据标准化至同一值域范围;步骤(2)②中,建立反演初始模型时:在使用geoview软件的STRATA模块进行反演计算时,反演参数的设置是:10-15Hz高截频;叠代次数10~20次;采样率为1ms~2ms;最大阻抗变化范围为25%~50%;预白化率为1%;运算块大小为1ms~2ms,该运算块大小与采样率相同;比例因子为1;反演类型为多道反演,Inline方向10~20道,Xline方向10~20道;步骤(2)④中,岩石物性参数分析利用geoview软件的crossplot模块。进一步的,如上所述的一种砂岩型铀矿地震综合预测方法,步骤(3)中,具体步骤如下:①基于三维地震纯波数据输入Landmark地震解释软件,将解释得到的三维地层数据导入geoview软件中,将步骤(2)的三维岩性数据体导入geoview软件;②确定目标层范围,利用geoview软件的tracemath模块在目标层顶底范围内,进行循环统计计算,得到地层含砂率分布数据RATIO(x,y);③使用“tracemath”模块需要编写代码,具体代码如下:据此可得到目标地层的含砂率RATIO(x,y)分布数据;④对步骤(3)②中的含砂率RATIO(x,y)数据,进行等值线成图。进一步的,如上所述的一种砂岩型铀矿地震综合预测方法,步骤(5)中,具体步骤如下:①利用步骤(4)得到的转换关系式,利用geoview软件的tracemath模块计算本文档来自技高网...

【技术保护点】
一种砂岩型铀矿地震综合预测方法,其特征在于:该方法包括以下步骤:(1)在研究区,采集一套砂岩型铀矿三维地震纯波数据;采集野外三维地震数据,对野外三维地震数据依次进行处理以得到三维地震纯波数据;(2)使用基于模型反演法对步骤(1)的三维地震纯波数据进行反演计算,得到三维波阻抗数据体,进而计算三维岩性数据体;(3)求取目标层位砂体的含砂率分布数据,并绘制成等值线图;利用步骤(2)得到的三维岩性数据体,利用geoview软件的trace math模块计算目标地层的含砂率数据,并绘制平面等值线图;(4)求取孔隙度数据与波阻抗数据的关系①在研究区测井曲线中有孔隙度数据、声波数据、密度数据三类测井数据的情况下:使用excel软件或geoview软件的cross plot模块进行孔隙度和波阻抗的交会分析;分析时,设置x横轴为波阻抗数据,波阻抗数据=声波和密度数据的乘积;y纵轴为孔隙度数据;然后通过excel软件或geoview软件的cross plot模块的线性拟合工具拟合出波阻抗转孔隙度的转换式y=ax+b,y为孔隙度数据,x为波阻抗数据,a,b为需要拟合的参数;②在研究区无孔隙度测井数据的情况下:使用下述的公式(1),将波阻抗转换为孔隙度数据:...

【技术特征摘要】
1.一种砂岩型铀矿地震综合预测方法,其特征在于:该方法包括以下步骤:(1)在研究区,采集一套砂岩型铀矿三维地震纯波数据;采集野外三维地震数据,对野外三维地震数据依次进行处理以得到三维地震纯波数据;(2)使用基于模型反演法对步骤(1)的三维地震纯波数据进行反演计算,得到三维波阻抗数据体,进而计算三维岩性数据体;(3)求取目标层位砂体的含砂率分布数据,并绘制成等值线图;利用步骤(2)得到的三维岩性数据体,利用geoview软件的tracemath模块计算目标地层的含砂率数据,并绘制平面等值线图;(4)求取孔隙度数据与波阻抗数据的关系①在研究区测井曲线中有孔隙度数据、声波数据、密度数据三类测井数据的情况下:使用excel软件或geoview软件的crossplot模块进行孔隙度和波阻抗的交会分析;分析时,设置x横轴为波阻抗数据,波阻抗数据=声波和密度数据的乘积;y纵轴为孔隙度数据;然后通过excel软件或geoview软件的crossplot模块的线性拟合工具拟合出波阻抗转孔隙度的转换式y=ax+b,y为孔隙度数据,x为波阻抗数据,a,b为需要拟合的参数;②在研究区无孔隙度测井数据的情况下:使用下述的公式(1),将波阻抗转换为孔隙度数据:其中,AC固体骨架表示岩石固体骨架的声波时差,AC流体表示岩石孔隙间流体的声波时差,IMP为波阻抗数据,POR为孔隙度数据;(5)利用步骤(4)求取的波阻抗转孔隙度的转换关系,将波阻抗数据转换为孔隙度数据,得到目标地层孔隙度数据分布,并绘制等值线图;(6)利用三维地震纯波数据,提取地震均方根振幅、瞬时相位、弧长3种地震属性,对此三种地震属性进行聚类分析,基于钻孔资料建立约束条件,回归拟合出一种新的地震属性组合;(7)交会分析上述目标地层含砂率数据、孔隙度分布数据以及地震属性组合分布数据,综合预测砂岩型铀矿成矿有利地段;综合预测砂岩型铀矿成矿有利地段是指对目标地层含砂率分布图中数值大于X的区域、孔隙度值大于Y的区域、地震属性组合值大于Z的区域进行交会分析研究,据此再进行成矿有利地段的综合预测。2.如权利要求1所述的一种砂岩型铀矿地震综合预测方法,其特征在于:步骤(1)中,通过地震仪采集野外三维地震数据,对野外三维地震数据进行处理的过程是依次进行静校正、去噪、振幅补偿、反褶积、动校叠加、偏移处理。3.如权利要求1所述的一种砂岩型铀矿地震综合预测方法,其特征在于:步骤(2)中,具体步骤如下:①初始模型的建立:首先收集研究区的钻井数据,对其中的声波曲线和密度曲线进行平滑处理和标准化处理;②利用geoview软件的STRATA模块建立反演初始模型;③波阻抗反演计算:对步骤(1)得到的三维地震纯波数据进行反演计算,得到三维波阻抗数据体;④确定波阻抗转岩性的门槛值:每个研究区的波阻抗转岩性的门槛值不尽相同,门槛值的确定基于工区岩石物性参数分析;⑤岩性数据计算:基于步骤(2)④得到的波阻抗转岩性的门槛值,将步骤(2)③得到的三维波阻抗数据体转换为三维岩性数据体,三维岩性数据体的形式为LITH(x,y,t),其中x表示三维地震数据的Inline号,y表示三维地震数据的Xline号,t表示时间;LITH(x,y,t)数据体的取值为0或1,0表示泥岩,1表示砂岩。4.如权利要求3所述的一种砂岩型铀矿地震综合预测方法,其特征在于:步骤(2)①中,平滑处理时采用3点或5点平滑处理;密度数据的3点或5点滑动平均处理分别如以下公式所示:●密度3点滑动平均公式:dden(i)=(dden(i-1)+dden(i)+dden(i+1))/3其中,di代表某个采样点的密度值,di-1为该采样点的前一个采样点的密度值,di+1为该采样点的后一个采样点的密度值;●密度5点滑动平均公式:dden(i)=(dden(i-2)+dden(i-1)+dden(i)+dden(i+1)+dden(i+2))/5其中,di代表某采样点的密度值,di-2为该采样点的前两个采样点的密度值,di-1为该采样点的前一个采样点的密度值,di+1为该采样点的后一个采样点的密度值,di+2为该采样点的后两个采样点的密度值;声波数据的3点或5点滑动平均处理分别如以下公式所示:●声波3点滑动平均公式:dson(i)=(dson(i-1)+dson(i)+dson(i+1))/3其中,di代表某个采样点的声波值,di-1为该采样点的前一个采样点的声波值,di+1为该采样点的后一个采样点的声波值;●声波5点滑动平均公式:dson(i)=(dson(i-2)+dson(i-1)+dson(i)+dson(i+1)+dson(i+2))/5其中,di代表某采样点的声波值,di-2为该采样点的前两个采样点的声波值,di-1为该采样点的前一个采样点的声波值,di+1为该采样点的后一个采样点的声波值,di+2为该采样点的后两个采样点的声波值;步骤(2)①中,标准化处理时使用geoview软件的loggingnomalize模块进行处理,将用于建立初始模型的井的声波和密度数据标准化至同一值域范围;步骤(2)②中,建立反演初始模型时:在使用geoview软件的STRATA模块进行反演计算时,反演参数的设置是:10-15Hz高截频;叠代次数10~20次;采样率为1ms~2ms;最大阻抗变化范围为25%~50%;预白化率为1%;运算块大小为1ms~2ms,该运算块大小与采样率相同;比例因子为1;反演类型为多道反演,Inline方向10~20道,Xline方向10~20道;步骤(2)④中,岩石物性参数分析利用geoview软件的crossplot模块。5.如权利要求1所述的一种砂岩型铀矿地震综合预测方法,其特征在于:步骤(3)中,具体步骤如下:①基于三维地震纯波数据输入Landmark地震解释软件,将解释得到的三维地层数据导入geoview软件中,将步骤(2)的三维岩性数据体导入geoview软件;②确定目标层范围,利用geoview软件的tracemath模块在目标层顶底范围内,进行循环统计计算,得到地层含砂率分布数据RATIO(x,y);③使用“tracemath”模块需要编写代码,具体代码如下:据此可得到目标地层的含砂率RATIO(x,y)分布数据;④对步骤(3)②中的含砂率RATIO(x,y)数据,进行等值线成图。6.如权利要求1所述的一种砂岩型铀矿地震综合预测方法,其特征在于:步骤(5)中,具体步骤如下:①利用步骤(4)得到的转换关系式,利用geoview软件的tracemath模块计算得到孔隙度三维数据体POR(x,y,t),x表示三维地震数据的Inline号,y表示三维地震数据的Xline号,t表示时间;②确定目标层范围,利用geoview软件的tracemath模块在目标层顶底范围内,进行循环统计计算,得到目标地层平均孔隙度分布数据POR_AVA(x,y),x表示三维地震数据的Inline号,y表示三维地震数据的Xline号;使用tracemath模块编写的代码如下:③据此可得到目标地层的孔隙度POR_AVA(x,y)分布数据;④对步骤(5)③中得到的孔隙度POR_AVA(x,y)数据,进行等值线成图。7.如权利要求1所述的一种砂岩型铀矿地震综合预测方法,其特征在于:步骤(6)中,将步骤(1)三维地震纯波数据输入至Landmark软件中,确定目标储层的位置,在此位置上下各20ms厚度内提取地震均方根振幅、瞬时相...

【专利技术属性】
技术研发人员:吴曲波黄伟传
申请(专利权)人:核工业北京地质研究院
类型:发明
国别省市:北京,11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1