CN111983721A - 一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统 - Google Patents
一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统 Download PDFInfo
- Publication number
- CN111983721A CN111983721A CN202010867940.6A CN202010867940A CN111983721A CN 111983721 A CN111983721 A CN 111983721A CN 202010867940 A CN202010867940 A CN 202010867940A CN 111983721 A CN111983721 A CN 111983721A
- Authority
- CN
- China
- Prior art keywords
- data
- sandstone
- determining
- wave impedance
- mud
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Abstract
本发明涉及一种砂岩铀矿“泥‑砂‑泥”地质结构识别方法及系统。该方法包括对三维地震数据进行预处理,确定三维地震纯波数据和三维地震成果数据;根据三维地震成果数据和测井数据确定目的层的层位和构造数据;根据测井数据确定砂岩和泥岩的波阻抗范围;根据三维地震纯波数据、三维地震成果数据和测井数据确定目的层的波阻抗值;根据砂岩和泥岩的波阻抗范围以及波阻抗值确定砂岩和泥岩的分布范围;根据砂岩的分布范围确定有效砂岩;根据泥岩的分布范围确定有效泥岩;根据有效砂岩和有效泥岩的分布位置识别砂岩型铀矿有利成矿“泥‑砂‑泥”地层结构;本发明能够快速识别地下的有利成矿的“泥‑砂‑泥”地质结构,加快砂岩铀矿的勘探周期。
Description
技术领域
本发明涉及砂岩型铀矿勘查领域,特别是涉及一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统。
背景技术
在砂岩型铀矿的勘查中,砂岩是砂岩型铀矿的赋存空间,上覆和下伏的泥岩层可以成为良好的隔水层,为砂岩型铀矿的成矿及后期开采提供了条件,因此如何快速准确识别砂岩铀矿有利成矿的“泥-砂-泥”地质结构,对于砂岩型铀矿成矿环境的研究非常重要。
而现有技术中还不能快速识别地下的有利成矿的“泥-砂-泥”地质结构,因此,需要一种砂岩铀矿“泥-砂-泥”地质结构识别方法,为砂岩铀矿有利成矿地层结构的勘查提供技术方法,加快砂岩铀矿的勘探周期。
发明内容
本发明的目的是提供一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统,能够快速识别地下的有利成矿的“泥-砂-泥”地质结构,加快砂岩铀矿的勘探周期。
为实现上述目的,本发明提供了如下方案:
一种砂岩铀矿“泥-砂-泥”地质结构识别方法,包括:
获取待勘查区的三维地震数据;
对所述三维地震数据进行预处理,确定三维地震纯波数据和三维地震成果数据;所述三维地震纯波数据用于地震波阻抗反演;所述三维地震成果数据用于地震资料解释;所述预处理包括解编、观测系统定义、静校正、去噪、振幅补偿、反褶积、速度分析、剩余静校正、叠加和偏移处理;
获取所述待勘查区的施工测井的测井数据;所述测井数据包括井径、电阻率、声波时差、电流、密度、自然电位、定量伽马和岩性;
根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据;
根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围;
根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值;
根据所述砂岩的波阻抗值范围、所述泥岩的波阻抗范围以及所述目的层的波阻抗值确定砂岩的分布范围和泥岩的分布范围;
根据所述砂岩的分布范围确定有效砂岩;所述有效砂岩为厚度20-60m,横向展布面积大于0.01km2的砂岩;
根据所述泥岩的分布范围确定有效泥岩;所述有效泥岩为厚度大于5m,横向展布面积大于0.01km2的泥岩;
根据所述有效砂岩和所述有效泥岩的分布位置识别所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构;所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构为所述有效砂岩上部和下部均覆有所述有效泥岩。
可选的,所述根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据,具体包括:
对所述测井数据进行合成地震记录,得到所述目的层的层位的标定;
对所述三维地震成果数据进行层位追踪,以及根据所述三维地震成果数据的同向轴特征识别断裂位置,确定所述目的层的层位和构造数据。
可选的,所述根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围,具体包括:
根据所述测井数据,采用公式波阻抗=密度/声波时差确定地层波阻抗值;
对所述测井数据进行岩石物理分析,确定所述目的层的砂岩和泥岩的波阻抗差异;
根据所述目的层的砂岩和泥岩的波阻抗差异确定所述砂岩的波阻抗值范围和所述泥岩的波阻抗范围。
可选的,所述根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值,具体包括:
利用所述三维地震纯波数据和所述测井数据合成地震记录;
从所述三维地震纯波数据提取地震子波;
根据所述三维地震成果数据和所述测井数据建立反演初始模型;
利用所述三维地震纯波数据、所述地震子波以及所述反演初始模型进行地震波阻抗反演,确定所述目的层的波阻抗值。
一种砂岩铀矿“泥-砂-泥”地质结构识别系统,包括:
三维地震数据获取模块,用于获取待勘查区的三维地震数据;
预处理模块,用于对所述三维地震数据进行预处理,确定三维地震纯波数据和三维地震成果数据;所述三维地震纯波数据用于地震波阻抗反演;所述三维地震成果数据用于地震资料解释;所述预处理包括解编、观测系统定义、静校正、去噪、振幅补偿、反褶积、速度分析、剩余静校正、叠加和偏移处理;
测井数据获取模块,用于获取所述待勘查区的施工测井的测井数据;所述测井数据包括井径、电阻率、声波时差、电流、密度、自然电位、定量伽马和岩性;
层位和构造数据确定模块,用于根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据;
波阻抗范围确定模块,用于根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围;
目的层的波阻抗值确定模块,用于根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值;
分布范围确定模块,用于根据所述砂岩的波阻抗值范围、所述泥岩的波阻抗范围以及所述目的层的波阻抗值确定砂岩的分布范围和泥岩的分布范围;
有效砂岩确定模块,用于根据所述砂岩的分布范围确定有效砂岩;所述有效砂岩为厚度20-60m,横向展布面积大于0.01km2的砂岩;
有效泥岩确定模块,用于根据所述泥岩的分布范围确定有效泥岩;所述有效泥岩为厚度大于5m,横向展布面积大于0.01km2的泥岩;
地层结构识别模块,用于根据所述有效砂岩和所述有效泥岩的分布位置识别所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构;所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构为所述有效砂岩上部和下部均覆有所述有效泥岩。
可选的,所述层位和构造数据确定模块具体包括:
层位的标定确定单元,用于对所述测井数据进行合成地震记录,得到所述目的层的层位的标定;
层位和构造数据确定单元,用于对所述三维地震成果数据进行层位追踪,以及根据所述三维地震成果数据的同向轴特征识别断裂位置,确定所述目的层的层位和构造数据。
可选的,所述波阻抗范围确定模块具体包括:
地层波阻抗值确定单元,用于根据所述测井数据,采用公式波阻抗=密度/声波时差确定地层波阻抗值;
波阻抗差异确定单元,用于对所述测井数据进行岩石物理分析,确定所述目的层的砂岩和泥岩的波阻抗差异;
波阻抗范围确定单元,用于根据所述目的层的砂岩和泥岩的波阻抗差异确定所述砂岩的波阻抗值范围和所述泥岩的波阻抗范围。
可选的,所述目的层的波阻抗值确定模块具体包括:
地震记录确定单元,用于利用所述三维地震纯波数据和所述测井数据合成地震记录;
地震子波提取单元,用于从所述三维地震纯波数据提取地震子波;
反演初始模型建立单元,用于根据所述三维地震成果数据和所述测井数据建立反演初始模型;
目的层的波阻抗值确定单元,用于利用所述三维地震纯波数据、所述地震子波以及所述反演初始模型进行地震波阻抗反演,确定所述目的层的波阻抗值。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明所提供的一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统,综合利用三维地震数据和测井数据,通过技术方法的实施,形成一种砂岩铀矿有利成矿“泥-砂-泥”地质结构的识别方法,通过开展地震波阻抗反演获取目的层的波阻抗分布情况,根据砂岩与泥岩的波阻抗差异,从地震波阻抗反演结果上快速识别出砂岩和泥岩的分布情况,再利用有利成矿地层结构的特征进行判定,最终筛选出砂岩铀矿有利成矿“泥-砂-泥”地层结构,加快砂岩铀矿的勘探周期。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所提供的一种砂岩铀矿“泥-砂-泥”地质结构识别方法流程示意图;
图2为本发明所提供的一种砂岩铀矿“泥-砂-泥”地质结构识别系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统,能够快速识别地下的有利成矿的“泥-砂-泥”地质结构,加快砂岩铀矿的勘探周期。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的一种砂岩铀矿“泥-砂-泥”地质结构识别方法流程示意图,如图1所示,本发明所提供的一种砂岩铀矿“泥-砂-泥”地质结构识别方法,包括:
S101,获取待勘查区的三维地震数据。其中,三维地震数据的采集范围,应根据待勘查区的范围确定,满覆盖区域应大于待勘查区的范围。
三维地震数据采集是通过地震仪采集地震野外实测数据,观测系统面元应≤10m×10m,覆盖次数不低于48次,最大偏移距根据勘探深度确定,计算方式为:最大偏移距≥勘探深度×1.5。检波器的主频不高于10Hz,震源可选用可控震源或炸药震源。
S102,对所述三维地震数据进行预处理,确定三维地震纯波数据和三维地震成果数据;所述三维地震纯波数据用于地震波阻抗反演;所述三维地震成果数据用于地震资料解释;所述预处理包括解编、观测系统定义、静校正、去噪、振幅补偿、反褶积、速度分析、剩余静校正、叠加和偏移处理。具体利用地震数据处理软件进行预处理。
S103,获取所述待勘查区的施工测井的测井数据;所述测井数据包括井径、电阻率、声波时差、电流、密度、自然电位、定量伽马和岩性。测井数据的采集过程是使用测井仪的参数探头向井下滑动,每隔0.05m测定一个数据。岩性是地质人员依据钻井取芯,通过观察描述各段岩芯的岩性,获取各地层的岩性。
S104,根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据。
S104具体包括:
对所述测井数据进行合成地震记录,得到所述目的层的层位的标定。具体利用地震数据解释软件。
对所述三维地震成果数据进行层位追踪,以及根据所述三维地震成果数据的同向轴特征识别断裂位置,确定所述目的层的层位和构造数据。
S105,根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围;
S105具体包括:
根据所述测井数据,采用公式波阻抗=密度/声波时差确定地层波阻抗值。
对所述测井数据进行岩石物理分析,确定所述目的层的砂岩和泥岩的波阻抗差异。岩石物理分析是通过波阻抗和其他参数的交会实现,交会图的横坐标轴为波阻抗,纵坐标轴为定量伽玛或电阻率等其他参数,点位颜色为岩性,将目的层相应的参数的测井数据投影到交会图上,根据点位颜色可以划分出砂岩区和泥岩区,从横坐标轴上分别读取砂岩区和泥岩区的波阻抗值范围。
根据所述目的层的砂岩和泥岩的波阻抗差异确定所述砂岩的波阻抗值范围和所述泥岩的波阻抗范围。
S106,根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值。
S106具体包括:
利用所述三维地震纯波数据和所述测井数据合成地震记录。具体利用地震波阻抗反演软件。
从所述三维地震纯波数据提取地震子波。
根据所述三维地震成果数据和所述测井数据建立反演初始模型。
利用所述三维地震纯波数据、所述地震子波以及所述反演初始模型进行地震波阻抗反演,确定所述目的层的波阻抗值。
S107,根据所述砂岩的波阻抗值范围、所述泥岩的波阻抗范围以及所述目的层的波阻抗值确定砂岩的分布范围和泥岩的分布范围。
对目的层的波阻抗反演结果的每个波阻抗数据进行判断,当该数值处于砂岩的波阻抗值范围内时,定义该位置为砂岩,当该数据处于泥岩的波阻抗值范围内时,定义该位置为泥岩,当该数据既不处于砂岩的波阻抗值范围内又不处于泥岩的波阻抗值范围内时,定义该位置为空值,最终获取目的层的砂岩和泥岩的分布范围。
S108,根据所述砂岩的分布范围确定有效砂岩;所述有效砂岩为厚度20-60m,横向展布面积大于0.01km2的砂岩。
S109,根据所述泥岩的分布范围确定有效泥岩;所述有效泥岩为厚度大于5m,横向展布面积大于0.01km2的泥岩。
具体的确定过程为:
砂岩的纵向厚度的计算方法是:首先读取砂岩在纵向上的数值个数(Num),用该数值个数乘以数据的采样间隔Sample,则获得砂岩的时间厚度H_time,该时间厚度为地震波的双程旅行时,用该双程旅行时除2并乘以砂岩速度Velocity,即可得到该砂岩的实际厚度H_depth。砂岩速度Velocity通过测井数据计算出,为相应地层的砂岩声波时差的倒数。举以下实例,假设Number=30,Sample=1ms,Velocity=2000m/s,则有H_time=Number×Sample=30ms,则该砂岩的实际厚度H_depth=H_time/2×Velocity=30ms/2×2000m/s=30m。
泥岩的纵向厚度的计算方法是:首先读取泥岩在纵向上的数值个数(Num),用该数值个数乘以数据的采样间隔Sample,则获得泥岩的时间厚度H_time,该时间厚度为地震波的双程旅行时,用该双程旅行时除2并乘以泥岩速度Velocity,即可得到该泥岩的实际厚度H_depth。泥岩速度Velocity通过测井数据计算出,为相应地层的泥岩声波时差的倒数。举以下实例,假设Number=20,Sample=1ms,Velocity=2200m/s,则有H_time=Number×Sample=20ms,则该砂岩的实际厚度H_depth=H_time/2×Velocity=20ms/2×2200m/s=22m。
砂岩的横向分布面积的计算方法是:首先依次读取各个纵向上的数值的横向分布个数,假设个数分别为N1、N2、N3…,读取其中的最大值Number_max=SUM(N1,N2,N3…),用该值乘以三维地震勘探的面元DX,即可计算出砂岩的横向分布面积X=Number×DX。举以下实例,假设Number_max=200,DX=5米×5米,则可以计算出假设的砂岩的横向分布面积X=Number×DX=200×5米×5米=5000平方米=0.5km2。
泥岩的横向分布面积的计算方法是:首先依次读取各个纵向上的数值的横向分布个数,假设个数分别为N1、N2、N3…,读取其中的最大值Number_max=SUM(N1,N2,N3…),用该值乘以三维地震勘探的面元DX,即可计算出泥岩的横向分布面积X=Number×DX。举以下实例,假设Number_max=100,DX=5米×5米,则可以计算出假设的砂岩的横向分布面积X=Number×DX=100×5米×5米=2500平方米=0.25km2。
S110,根据所述有效砂岩和所述有效泥岩的分布位置识别所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构;所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构为所述有效砂岩上部和下部均覆有所述有效泥岩。
图2为本发明所提供的一种砂岩铀矿“泥-砂-泥”地质结构识别系统结构示意图,如图2所示,本发明所提供的一种砂岩铀矿“泥-砂-泥”地质结构识别系统,包括:三维地震数据获取模块201、预处理模块202、测井数据获取模块203、层位和构造数据确定模块204、波阻抗范围确定模块205、目的层的波阻抗值确定模块206、分布范围确定模块207、有效砂岩确定模块208、有效泥岩确定模块209和地层结构识别模块210。
三维地震数据获取模块201用于获取待勘查区的三维地震数据。
预处理模块202用于对所述三维地震数据进行预处理,确定三维地震纯波数据和三维地震成果数据;所述三维地震纯波数据用于地震波阻抗反演;所述三维地震成果数据用于地震资料解释;所述预处理包括解编、观测系统定义、静校正、去噪、振幅补偿、反褶积、速度分析、剩余静校正、叠加和偏移处理;
测井数据获取模块203用于获取所述待勘查区的施工测井的测井数据;所述测井数据包括井径、电阻率、声波时差、电流、密度、自然电位、定量伽马和岩性。
层位和构造数据确定模块204用于根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据。
波阻抗范围确定模块205用于根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围。
目的层的波阻抗值确定模块206用于根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值。
分布范围确定模块207用于根据所述砂岩的波阻抗值范围、所述泥岩的波阻抗范围以及所述目的层的波阻抗值确定砂岩的分布范围和泥岩的分布范围。
有效砂岩确定模块208用于根据所述砂岩的分布范围确定有效砂岩;所述有效砂岩为厚度20-60m,横向展布面积大于0.01km2的砂岩。
有效泥岩确定模块209用于根据所述泥岩的分布范围确定有效泥岩;所述有效泥岩为厚度大于5m,横向展布面积大于0.01km2的泥岩。
地层结构识别模块210用于根据所述有效砂岩和所述有效泥岩的分布位置识别所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构;所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构为所述有效砂岩上部和下部均覆有所述有效泥岩。
所述层位和构造数据确定模块204具体包括:层位的标定确定单元和层位和构造数据确定单元。
层位的标定确定单元用于对所述测井数据进行合成地震记录,得到所述目的层的层位的标定。
层位和构造数据确定单元用于对所述三维地震成果数据进行层位追踪,以及根据所述三维地震成果数据的同向轴特征识别断裂位置,确定所述目的层的层位和构造数据。
所述波阻抗范围确定模块205具体包括:地层波阻抗值确定单元、波阻抗差异确定单元和波阻抗范围确定单元。
地层波阻抗值确定单元用于根据所述测井数据,采用公式波阻抗=密度/声波时差确定地层波阻抗值。
波阻抗差异确定单元用于对所述测井数据进行岩石物理分析,确定所述目的层的砂岩和泥岩的波阻抗差异。
波阻抗范围确定单元用于根据所述目的层的砂岩和泥岩的波阻抗差异确定所述砂岩的波阻抗值范围和所述泥岩的波阻抗范围。
所述目的层的波阻抗值确定模块206具体包括:地震记录确定单元、地震子波提取单元、反演初始模型建立单元和目的层的波阻抗值确定单元。
地震记录确定单元用于利用所述三维地震纯波数据和所述测井数据合成地震记录。
地震子波提取单元用于从所述三维地震纯波数据提取地震子波。
反演初始模型建立单元用于根据所述三维地震成果数据和所述测井数据建立反演初始模型。
目的层的波阻抗值确定单元用于利用所述三维地震纯波数据、所述地震子波以及所述反演初始模型进行地震波阻抗反演,确定所述目的层的波阻抗值。本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种砂岩铀矿“泥-砂-泥”地质结构识别方法,其特征在于,包括:
获取待勘查区的三维地震数据;
对所述三维地震数据进行预处理,确定三维地震纯波数据和三维地震成果数据;所述三维地震纯波数据用于地震波阻抗反演;所述三维地震成果数据用于地震资料解释;所述预处理包括解编、观测系统定义、静校正、去噪、振幅补偿、反褶积、速度分析、剩余静校正、叠加和偏移处理;
获取所述待勘查区的施工测井的测井数据;所述测井数据包括井径、电阻率、声波时差、电流、密度、自然电位、定量伽马和岩性;
根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据;
根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围;
根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值;
根据所述砂岩的波阻抗值范围、所述泥岩的波阻抗范围以及所述目的层的波阻抗值确定砂岩的分布范围和泥岩的分布范围;
根据所述砂岩的分布范围确定有效砂岩;所述有效砂岩为厚度20-60m,横向展布面积大于0.01km2的砂岩;
根据所述泥岩的分布范围确定有效泥岩;所述有效泥岩为厚度大于5m,横向展布面积大于0.01km2的泥岩;
根据所述有效砂岩和所述有效泥岩的分布位置识别所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构;所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构为所述有效砂岩上部和下部均覆有所述有效泥岩。
2.根据权利要求1所述的一种砂岩铀矿“泥-砂-泥”地质结构识别方法,其特征在于,所述根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据,具体包括:
对所述测井数据进行合成地震记录,得到所述目的层的层位的标定;
对所述三维地震成果数据进行层位追踪,以及根据所述三维地震成果数据的同向轴特征识别断裂位置,确定所述目的层的层位和构造数据。
3.根据权利要求1所述的一种砂岩铀矿“泥-砂-泥”地质结构识别方法,其特征在于,所述根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围,具体包括:
根据所述测井数据,采用公式波阻抗=密度/声波时差确定地层波阻抗值;
对所述测井数据进行岩石物理分析,确定所述目的层的砂岩和泥岩的波阻抗差异;
根据所述目的层的砂岩和泥岩的波阻抗差异确定所述砂岩的波阻抗值范围和所述泥岩的波阻抗范围。
4.根据权利要求1所述的一种砂岩铀矿“泥-砂-泥”地质结构识别方法,其特征在于,所述根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值,具体包括:
利用所述三维地震纯波数据和所述测井数据合成地震记录;
从所述三维地震纯波数据提取地震子波;
根据所述三维地震成果数据和所述测井数据建立反演初始模型;
利用所述三维地震纯波数据、所述地震子波以及所述反演初始模型进行地震波阻抗反演,确定所述目的层的波阻抗值。
5.一种砂岩铀矿“泥-砂-泥”地质结构识别系统,其特征在于,包括:
三维地震数据获取模块,用于获取待勘查区的三维地震数据;
预处理模块,用于对所述三维地震数据进行预处理,确定三维地震纯波数据和三维地震成果数据;所述三维地震纯波数据用于地震波阻抗反演;所述三维地震成果数据用于地震资料解释;所述预处理包括解编、观测系统定义、静校正、去噪、振幅补偿、反褶积、速度分析、剩余静校正、叠加和偏移处理;
测井数据获取模块,用于获取所述待勘查区的施工测井的测井数据;所述测井数据包括井径、电阻率、声波时差、电流、密度、自然电位、定量伽马和岩性;
层位和构造数据确定模块,用于根据所述三维地震成果数据和所述测井数据确定目的层的层位和构造数据;
波阻抗范围确定模块,用于根据所述测井数据确定砂岩的波阻抗值范围和泥岩的波阻抗范围;
目的层的波阻抗值确定模块,用于根据所述三维地震纯波数据、三维地震成果数据和所述测井数据确定所述目的层的波阻抗值;
分布范围确定模块,用于根据所述砂岩的波阻抗值范围、所述泥岩的波阻抗范围以及所述目的层的波阻抗值确定砂岩的分布范围和泥岩的分布范围;
有效砂岩确定模块,用于根据所述砂岩的分布范围确定有效砂岩;所述有效砂岩为厚度20-60m,横向展布面积大于0.01km2的砂岩;
有效泥岩确定模块,用于根据所述泥岩的分布范围确定有效泥岩;所述有效泥岩为厚度大于5m,横向展布面积大于0.01km2的泥岩;
地层结构识别模块,用于根据所述有效砂岩和所述有效泥岩的分布位置识别所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构;所述砂岩型铀矿有利成矿“泥-砂-泥”地层结构为所述有效砂岩上部和下部均覆有所述有效泥岩。
6.根据权利要求5所述的一种砂岩铀矿“泥-砂-泥”地质结构识别系统,其特征在于,所述层位和构造数据确定模块具体包括:
层位的标定确定单元,用于对所述测井数据进行合成地震记录,得到所述目的层的层位的标定;
层位和构造数据确定单元,用于对所述三维地震成果数据进行层位追踪,以及根据所述三维地震成果数据的同向轴特征识别断裂位置,确定所述目的层的层位和构造数据。
7.根据权利要求5所述的一种砂岩铀矿“泥-砂-泥”地质结构识别系统,其特征在于,所述波阻抗范围确定模块具体包括:
地层波阻抗值确定单元,用于根据所述测井数据,采用公式波阻抗=密度/声波时差确定地层波阻抗值;
波阻抗差异确定单元,用于对所述测井数据进行岩石物理分析,确定所述目的层的砂岩和泥岩的波阻抗差异;
波阻抗范围确定单元,用于根据所述目的层的砂岩和泥岩的波阻抗差异确定所述砂岩的波阻抗值范围和所述泥岩的波阻抗范围。
8.根据权利要求5所述的一种砂岩铀矿“泥-砂-泥”地质结构识别系统,其特征在于,所述目的层的波阻抗值确定模块具体包括:
地震记录确定单元,用于利用所述三维地震纯波数据和所述测井数据合成地震记录;
地震子波提取单元,用于从所述三维地震纯波数据提取地震子波;
反演初始模型建立单元,用于根据所述三维地震成果数据和所述测井数据建立反演初始模型;
目的层的波阻抗值确定单元,用于利用所述三维地震纯波数据、所述地震子波以及所述反演初始模型进行地震波阻抗反演,确定所述目的层的波阻抗值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010867940.6A CN111983721A (zh) | 2020-08-26 | 2020-08-26 | 一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010867940.6A CN111983721A (zh) | 2020-08-26 | 2020-08-26 | 一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111983721A true CN111983721A (zh) | 2020-11-24 |
Family
ID=73443480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010867940.6A Pending CN111983721A (zh) | 2020-08-26 | 2020-08-26 | 一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111983721A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640873A (zh) * | 2021-08-18 | 2021-11-12 | 核工业二0八大队 | 砂岩型铀矿叠前地震预测方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454685A (zh) * | 2013-08-09 | 2013-12-18 | 中国石油天然气股份有限公司 | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 |
CN103954995A (zh) * | 2014-04-22 | 2014-07-30 | 核工业北京地质研究院 | 一种砂岩型铀矿勘查中的砂体识别方法 |
CN105353407A (zh) * | 2015-10-28 | 2016-02-24 | 中国石油化工股份有限公司 | 一种叠后地震波阻抗反演方法 |
WO2016135505A1 (en) * | 2015-02-26 | 2016-09-01 | Ireson Richard | A method of producing and utilising high resolution impedance logs derived from vsp data for use in assessing an oilfield subterranean formation |
CN107576982A (zh) * | 2017-09-01 | 2018-01-12 | 核工业北京地质研究院 | 一种砂岩型铀矿地震综合预测方法 |
WO2019062655A1 (zh) * | 2017-09-27 | 2019-04-04 | 中国石油天然气股份有限公司 | 薄夹层的确定方法和装置 |
CN111257926A (zh) * | 2018-12-03 | 2020-06-09 | 核工业二0八大队 | 一种利用老旧地震资料进行古河谷铀储层预测的方法 |
CN111352172A (zh) * | 2020-02-21 | 2020-06-30 | 吉林大学 | 一种用井震联合法获取铀异常在砂体中空间分布位置的方法 |
-
2020
- 2020-08-26 CN CN202010867940.6A patent/CN111983721A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454685A (zh) * | 2013-08-09 | 2013-12-18 | 中国石油天然气股份有限公司 | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 |
CN103954995A (zh) * | 2014-04-22 | 2014-07-30 | 核工业北京地质研究院 | 一种砂岩型铀矿勘查中的砂体识别方法 |
WO2016135505A1 (en) * | 2015-02-26 | 2016-09-01 | Ireson Richard | A method of producing and utilising high resolution impedance logs derived from vsp data for use in assessing an oilfield subterranean formation |
CN105353407A (zh) * | 2015-10-28 | 2016-02-24 | 中国石油化工股份有限公司 | 一种叠后地震波阻抗反演方法 |
CN107576982A (zh) * | 2017-09-01 | 2018-01-12 | 核工业北京地质研究院 | 一种砂岩型铀矿地震综合预测方法 |
WO2019062655A1 (zh) * | 2017-09-27 | 2019-04-04 | 中国石油天然气股份有限公司 | 薄夹层的确定方法和装置 |
CN111257926A (zh) * | 2018-12-03 | 2020-06-09 | 核工业二0八大队 | 一种利用老旧地震资料进行古河谷铀储层预测的方法 |
CN111352172A (zh) * | 2020-02-21 | 2020-06-30 | 吉林大学 | 一种用井震联合法获取铀异常在砂体中空间分布位置的方法 |
Non-Patent Citations (3)
Title |
---|
吴曲波 等: "《砂岩型铀矿地震勘探中若干问题研究》", 《地球物理学进展》, vol. 32, no. 5, pages 2097 - 2105 * |
李子伟等: "地震勘探方法在伊犁盆地南缘砂岩型铀矿勘查中的应用研究", 《铀矿地质》, vol. 32, no. 04, pages 236 - 240 * |
李子伟等: "砂岩型铀矿波阻抗反演砂体识别技术研究", 《世界核地质科学》, vol. 34, no. 04 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640873A (zh) * | 2021-08-18 | 2021-11-12 | 核工业二0八大队 | 砂岩型铀矿叠前地震预测方法及装置 |
CN113640873B (zh) * | 2021-08-18 | 2024-04-12 | 核工业二0八大队 | 砂岩型铀矿叠前地震预测方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Suzuki et al. | Case studies of electrical and electromagnetic methods applied to mapping active faults beneath the thick quaternary | |
Tronicke et al. | Crosshole traveltime tomography using particle swarm optimization: A near-surface field example | |
Cueto et al. | Karst‐induced sinkhole detection using an integrated geophysical survey: a case study along the Riyadh Metro Line 3 (Saudi Arabia) | |
CN104714254B (zh) | 一种基于瞬变电磁法探测煤矿积水采空区的方法和装置 | |
Ahmadi et al. | High-resolution 2D seismic imaging and forward modeling of a polymetallic sulfide deposit at Garpenberg, central Sweden | |
CN109799540A (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
Gammaldi et al. | Fluid accumulation zone by seismic attributes and amplitude versus offset analysis at Solfatara Volcano, Campi Flegrei, Italy | |
Wen et al. | Review of geophysical exploration on mined-out areas and water abundance | |
CN111852467B (zh) | 一种砂岩铀矿矿体延伸范围的圈定方法及系统 | |
Strahser et al. | Polarisation and slowness of seismoelectric signals: a case study | |
WO2013057484A2 (en) | Seismic survey data collection and processing | |
CN111983721A (zh) | 一种砂岩铀矿“泥-砂-泥”地质结构识别方法及系统 | |
Rapetsoa et al. | Cost‐effective in‐mine seismic experiments to image platinum deposits and associated geological structures at Maseve platinum mine, South Africa | |
Sagar et al. | Electrical Resistivity Tomography in Geotechnical Engineering Applications | |
Rapetsoa et al. | Seismic solutions utilizing existing in-mine infrastructure for mineral exploration: A case study from Maseve platinum mine, South Africa | |
Maufroy et al. | Travel time inversion from ground level to gallery: protocol for the characterization of P‐wave seismic signature in a fractured‐porous Urgonian platform at hectometric scale | |
Onyishi et al. | Source parameter imaging and Euler deconvolution of aeromagnetic anomalies over parts of the middle Benue trough, Nigeria | |
CN112731557A (zh) | 一种预测砂岩铀矿成矿有利区的综合地球物理方法 | |
Stucchi et al. | High‐resolution reflection seismic survey at the Patigno landslide, Northern Apennines, Italy | |
Bergman et al. | High-resolution reflection seismic imaging of the upper crust at Laxemar, southeastern Sweden | |
JP2002055172A (ja) | 地盤内空洞の調査方法 | |
Botter et al. | Seismic attribute analysis of a fault zone in the Thebe field, Northwest shelf, Australia | |
Stober et al. | Geophysical Methods, Exploration and Analysis | |
Bhattacharya | Application of geophysical techniques in groundwater management | |
Amalina et al. | Application of 2-D Resistivity Imaging and Seismic Refraction method in identifying the structural geological contact of sedimentary lithologies |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |