CN107144873A - 一种砂岩型铀矿三维地震数据观测方法 - Google Patents
一种砂岩型铀矿三维地震数据观测方法 Download PDFInfo
- Publication number
- CN107144873A CN107144873A CN201710236881.0A CN201710236881A CN107144873A CN 107144873 A CN107144873 A CN 107144873A CN 201710236881 A CN201710236881 A CN 201710236881A CN 107144873 A CN107144873 A CN 107144873A
- Authority
- CN
- China
- Prior art keywords
- mrow
- seismic
- sandstone
- msub
- respect
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种砂岩型铀矿的三维地震数据观测方法,特别是涉及一种用于探测砂岩型铀矿的三维地震数据采集和接收方法。包括以下步骤:(1)确定观测系统参数,需要确定的砂岩型铀矿三维地震数据采集观测系统的关键参数包括面元尺寸、覆盖次数、最大炮检距、方位参数、横向滚动距、接收线距、最大非纵距、最大偏移距;(2)确定地震信号激发技术;(3)接收技术,设置砂岩型铀矿地震勘探的检波器和其他仪器的接收参数;(4)部署检波器检波器、激发震源;(5)激发和接收地震信号。
Description
技术领域
本发明涉及一种砂岩型铀矿的三维地震数据观测方法,特别是涉及一种用于探测砂岩型铀矿的三维地震数据采集和接收方法。
背景技术
三维地震勘探技术运用于石油、煤炭的勘探和开发已日臻成熟,然而,目前尚未见关于砂岩型铀矿三维地震勘探技术开发和应用研究的报道。砂岩型铀矿的成矿因素和赋矿层位有别于石油、煤炭资源,正因这些特点,也导致砂岩型铀矿三维地震勘探存在如下难点:①砂岩型铀矿埋藏深度浅;②砂岩型铀矿与围岩的物性差异小;③目标层砂泥薄互层多。如果教条地使用石油煤炭系统的三维地震数据观测方法进行砂岩型铀矿地震数据采集,将导致采集的数据缺失浅层砂岩型铀矿有效信息,难以有效地识别目标砂体等不足。因此,需开发一种针对砂岩型铀矿的三维地震数据观测方法,提升勘探效果。
发明内容
本发明的目的在于提供一种砂岩型铀矿三维地震数据观测方法,使用该方法能够有效、经济地采集到包含砂岩型铀矿弱信息的三维地震原始数据,为后期砂岩型铀矿精细勘查提供技术和资料保障。
实现本发明目的的技术方案:
一种砂岩型铀矿三维地震数据观测方法,包括以下步骤:
(1)确定观测系统参数
需要确定的砂岩型铀矿三维地震数据采集观测系统的关键参数包括面元尺寸、覆盖次数、最大炮检距、方位参数、横向滚动距、接收线距、最大非纵距、最大偏移距;
(2)确定地震信号激发技术
地震信号激发技术主要是设置砂岩型铀矿地震勘探中的震源激发参数,涉及两种地震信号激发震源:炸药震源和可控机械震源;
(3)接收技术
设置砂岩型铀矿地震勘探的检波器和其他仪器的接收参数;
(4)部署检波器检波器、激发震源
根据步骤(1)中设计的观测参数部署检波器、激发震源的位置;
砂岩型铀矿三维地震观测系统部署原则是:
(a)检波器的布设方向垂直于地质构造延伸方向,激发震源的布设方向平行于构造延伸方向;
(b)检波器的布设方向和激发震源点的布设方向是垂直的;
(5)激发和接收地震信号
根据步骤(2)中确定的地震信号激发源的参数,使用步骤(2)中确定的地震信号激发源激发地震信号;
根据步骤(3)中确定的检波器和其他仪器的接收参数,启动地震数据采集系统接收信号。
进一步的,如上所述的一种砂岩型铀矿三维地震数据观测方法,
步骤(1)中,面元尺寸的确定方法如下:
A、满足横向分辨率的要求
定义在地震信号谱密度曲线上与最大值对应的频率为地震信号优势频率;在地震信号优势频率的波长内设置2-3个采样点,从而控制对应的面元满足系统需要的分辨率,具体公式如下:
x=Vint/(2*Fp)
式中:x-确定所求面元尺寸的第一上限值;Vint-目的层上覆地层层速度;Fp-目的层反射波主频,表示由由不同频率的波动叠合在一起构成的地震波中能量最强的波动所对应的频率成份;
公式中Vint、Fp根据工区地质资料确定;
B、满足最高无混叠频率的要求
Fmax表示每个倾斜地震反射波同相轴的一个偏移前可能的最高无混叠频率,它根据目的层反射波均方根速度Vrms、倾角θ、面元尺寸确定,具体公式如下:
y=Vrms/(4Fmax*Sinθ)
式中:y-确定所求面元尺寸的第二上限值;Vrms-目的层反射波均方根速度,Fmax-最高无混叠频率,θ-目的层倾角,表示每个倾斜地震反射波同相轴与水平方向的角度;
公式中Vrms、Fmax、θ根据工区地质资料确定;
C、通过上述A、B两步计算得到面元尺寸b的取值上限:
若x≤y,则b≤x;
若y≤x,则b≤y;
D、确定面元尺寸b的取值下限
确定面元尺寸b的取值下限为5m;
E、确定面元尺寸b的取值范围
若y≤x,5m≤b≤y;
若x≤y,5m≤b≤x。
F、在面元尺寸b的取值范围内选取数值确定面元尺寸。
进一步的,如上所述的一种砂岩型铀矿三维地震数据观测方法,
砂岩型铀矿地震勘探中的面元尺寸选择为以下几种之一:
面元尺寸是一组两个数,第一个数的含义是面元的纵向距离,其数值的两倍对应检波点距;第二个数的含义是面元的横向距离,其数值的两倍对应炮点距;
5m*5m,对应10m检波点距*10m炮点距;
7.5m*7.5m,对应15m检波点距*15m炮点距;
5*10m,对应10m检波点距*20m炮点距;
10*5m,对应20m检波点距*10m炮点距。
进一步的,如上所述的一种砂岩型铀矿三维地震数据观测方法,
(一)步骤(1)中,覆盖次数的确定方法如下:
首先使用克朗软件噪声分析模块估算出以往地震剖面中目的层信号的信噪比(S/N)raw,再根据新资料期望达到的(S/N)required,按照以下公式计算覆盖次数Foldrequired:
公式中,Foldold表示以往地震资料的覆盖次数;
根据要解决的技术问题限定在砂岩型铀矿地震勘探中新资料期望达到的(S/N)required:
如要解决地下地质体的构造问题,则(S/N)required≥1;
如要解决砂岩型铀矿勘查中砂体预测问题,则(S/N)required≥2;
(二)步骤(1)中,最大炮检距的确定方法如下:
最大炮检距用X表示,根据a、b、c三个条件确定其取值范围;
a.最大炮检距X在0.9倍至1.2倍主要目标深度之间;
b.动校拉伸率D与最大炮检距X的关系如下:
式中:T0-目的层双程反射时间,V-目标层地震均方根速度;
根据最大炮检距X满足动校拉伸率D小于12.5%的要求,确定最大炮检距X≥5VT0;
c.最大炮检距满足速度分析要求,根据如下公式确定:
式中:P-速度分析精度,根据系统设定确定,P<6%;
(三)步骤(1)中,方位角参数的确定方法如下:
方位角参数=横向最大偏移距/纵向最大偏移距,在砂岩型铀矿地震勘探领域方位角参数取值在0.9-1之间;
横向最大偏移距为地震勘探的常用参数,横向最大偏移距=((接收线距*(接收线数-1)/2)+1/2*接收线距-1/2*炮点距);
纵向最大偏移距为地震勘探的常用参数,纵向最大偏移距=(检波点距*(单线接收道数-1)/2);
综上,方位角=((接收线距*(接收线数-1)/2)+1/2*接收线距-1/2*炮点距)/(检波点距*(单线接收道数-1)/2);
根据上述公式,通过设置接收线距、接收线数、炮点距、检波点距、单线接收道数的取值,使得方位角参数落入0.9-1范围;
(四)接收线距的确定方法如下:
接收线距≤垂直入射时的菲涅尔带半径:
式中:R-菲涅尔带半径,根据工区地质资料确定;
(五)横向滚动距的确定方法如下:
横向滚动距=接收线距;
(六)最大横向偏移距的确定方法如下:
最大横向偏移距Ymax满足:
式中:Va-目标层反射波平均速度,根据工区地质资料确定;
(七)最大最小偏移距的确定方法如下:
最浅目标层的埋深Xm根据工区地质资料确定,最大最小偏移距小于Xm。
进一步的,如上所述的一种砂岩型铀矿三维地震数据观测方法,
最大炮检距在800-940m之间,覆盖次数在64~128之间,接收线距在60~80m之间。
进一步的,如上所述的一种砂岩型铀矿三维地震数据观测方法,
接收线数为16线,单线接收道数为128道,面元为5m*10m,覆盖次数为128次,检波点距为10m,炮点距为20m,接收线距为80m,炮线距为40m,纵向最大偏移距为635m,最大横向偏移距为630m,最大最小偏移距为89m,最大炮检距为894m,横向滚动距离为80m,横纵比为0.992。
进一步的,如上所述的一种砂岩型铀矿三维地震数据观测方法,
接收线数为16线,单线接收道数为64道,面元为10m*5m,覆盖次数为128次,检波点距为20m,炮点距为20m,接收线距为80m,炮线距为40m,纵向最大偏移距为630m,最大横向偏移距为635m,最大最小偏移距为89m,最大炮检距为894m,横向滚动距离为80m,横纵比为0.992。
进一步的,如上所述的一种砂岩型铀矿三维地震数据观测方法,
①步骤(2)中,地震信号激发震源采用炸药震源时:
激发炮井的深度选择在高速层的埋藏深度以下3m,炸药选用高密度高爆速炸药,药量选用1.5kg~2kg;
其中,采用锤击多次叠加小折射法对高速层埋藏深度进行探测,锤击小折射法的观测系统参数为:
仪器使用工程地震仪;采样间隔0.25ms;记录长度0.5s;记录格式为SEG-2;激发源为20磅大锤,采用大于10次的能量叠加的方式激发;观测方式为相遇观测;接收道数24道;检波器使用主频为10Hz~40Hz的单只纵波检波器;最小偏移距2~6m;24道检波器之间的距离设置为:1m,1m,2m,2m,3m,5m,5m,8m,8m,10m,10m,10m,10m,10m,8m,8m,5m,5m,3m,2m,2m,1m,1m;使用截距时间法对采集的数据进行处理和解释;
②步骤(2)中,地震信号激发震源采用可控机械震源时:
反射波信噪比≥1时,可控震源台次使用1台*2次至1台*4次;
反射波信噪比<1时,可控震源台次使用2台*2次;
扫描频率使用宽频激发,低频低于6Hz,高频大于100Hz,扫描时间8s~14s,线性扫频激发;
步骤(3)接收技术中,设置砂岩型铀矿地震勘探的接收参数如下:
①选用主频≤10Hz的压电式加速度检波器,采用4~10只检波器点式组合的方式进行地震数据接收,并埋置于地下20~30cm的硬实土壤中;
②选用真24位地震数据采集系统进行地震信号接收,仪器的采样间隔为0.5ms-1ms,记录长度为3s-4s,前放增益为0-12dB。
本发明的有益技术效果在于:首次提出了一套系统的砂岩型铀矿三维地震数据观测方法。主要包括观测系统参数设置、激发技术和接收技术。
利用此砂岩型铀矿三维地震数据观测方法,可收集一套有效反应砂岩型铀矿弱信息的三维地震数据,通过此数据精确解释铀矿成矿环境要素信息,指明研究区铀矿找矿方向。
具体发明内容
下面结合具体实施例对本发明技术方案进行详细说明。
本发明一种砂岩型铀矿三维地震数据观测方法,包括以下步骤:
(1)确定观测系统参数
需要确定的砂岩型铀矿三维地震数据采集观测系统的关键参数包括面元尺寸、覆盖次数、最大炮检距、方位参数、横向滚动距、接收线距、最大非纵距、最大偏移距;
步骤(1)中,面元尺寸的确定方法如下:
A、满足横向分辨率的要求
定义在地震信号谱密度曲线上与最大值对应的频率为地震信号优势频率;在地震信号优势频率的波长内设置2-3个采样点,从而控制对应的面元满足系统需要的分辨率,具体公式如下:
x=Vint/(2*Fp)
式中:x-确定所求面元尺寸的第一上限值;Vint-目的层上覆地层层速度;Fp-目的层反射波主频,表示由由不同频率的波动叠合在一起构成的地震波中能量最强的波动所对应的频率成份;
公式中Vint、Fp根据工区地质资料确定;
B、满足最高无混叠频率的要求
Fmax表示每个倾斜地震反射波同相轴的一个偏移前可能的最高无混叠频率,它根据目的层反射波均方根速度Vrms、倾角θ、面元尺寸确定,具体公式如下:
y=Vrms/(4Fmax*Sinθ)
式中:y-确定所求面元尺寸的第二上限值;Vrms-目的层反射波均方根速度,Fmax-最高无混叠频率,θ-目的层倾角,表示每个倾斜地震反射波同相轴与水平方向的角度;
公式中Vrms、Fmax、θ根据工区地质资料确定;
C、通过上述A、B两步计算得到面元尺寸b的取值上限:
若x≤y,则b≤x;
若y≤x,则b≤y;
D、确定面元尺寸b的取值下限
确定面元尺寸b的取值下限为5m;
E、确定面元尺寸b的取值范围
若y≤x,5m≤b≤y;
若x≤y,5m≤b≤x。
F、在面元尺寸b的取值范围内选取数值确定面元尺寸。
砂岩型铀矿地震勘探中的面元尺寸选择为以下几种之一:
面元尺寸是一组两个数,第一个数的含义是面元的纵向距离,其数值的两倍对应检波点距;第二个数的含义是面元的横向距离,其数值的两倍对应炮点距;
5m*5m,对应10m检波点距*10m炮点距;
7.5m*7.5m,对应15m检波点距*15m炮点距;
5*10m,对应10m检波点距*20m炮点距;
10*5m,对应20m检波点距*10m炮点距。
(一)步骤(1)中,覆盖次数的确定方法如下:
首先使用克朗软件噪声分析模块估算出以往地震剖面中目的层信号的信噪比(S/N)raw,再根据新资料期望达到的(S/N)required,按照以下公式计算覆盖次数Foldrequired:
公式中,Foldold表示以往地震资料的覆盖次数;
根据要解决的技术问题限定在砂岩型铀矿地震勘探中新资料期望达到的(S/N)required:
如要解决地下地质体的构造问题,则(S/N)required≥1;
如要解决砂岩型铀矿勘查中砂体预测问题,则(S/N)required≥2;
(二)步骤(1)中,最大炮检距的确定方法如下:
最大炮检距用X表示,根据a、b、c三个条件确定其取值范围;
a.最大炮检距X在0.9倍至1.2倍主要目标深度之间;
b.动校拉伸率D与最大炮检距X的关系如下:
式中:T0-目的层双程反射时间,V-目标层地震均方根速度;
根据最大炮检距X满足动校拉伸率D小于12.5%的要求,确定最大炮检距X≥5VT0;
c.最大炮检距满足速度分析要求,根据如下公式确定:
式中:P-速度分析精度,根据系统设定确定,P<6%;
(三)步骤(1)中,方位角参数的确定方法如下:
方位角参数=横向最大偏移距/纵向最大偏移距,在砂岩型铀矿地震勘探领域方位角参数取值在0.9-1之间;
横向最大偏移距为地震勘探的常用参数,横向最大偏移距=((接收线距*(接收线数-1)/2)+1/2*接收线距-1/2*炮点距);
纵向最大偏移距为地震勘探的常用参数,纵向最大偏移距=(检波点距*(单线接收道数-1)/2);
综上,方位角=((接收线距*(接收线数-1)/2)+1/2*接收线距-1/2*炮点距)/(检波点距*(单线接收道数-1)/2);
根据上述公式,通过设置接收线距、接收线数、炮点距、检波点距、单线接收道数的取值,使得方位角参数落入0.9-1范围;
(四)接收线距的确定方法如下:
接收线距≤垂直入射时的菲涅尔带半径:
式中:R-菲涅尔带半径,根据工区地质资料确定;
(五)横向滚动距的确定方法如下:
横向滚动距=接收线距;
(六)最大横向偏移距的确定方法如下:
最大横向偏移距Ymax满足:
式中:Va-目标层反射波平均速度,根据工区地质资料确定;
(七)最大最小偏移距的确定方法如下:
最浅目标层的埋深Xm根据工区地质资料确定,最大最小偏移距小于Xm。
(2)确定地震信号激发技术
地震信号激发技术主要是设置砂岩型铀矿地震勘探中的震源激发参数,涉及两种地震信号激发震源:炸药震源和可控机械震源;
①步骤(2)中,地震信号激发震源采用炸药震源时:
激发炮井的深度选择在高速层的埋藏深度以下3m,炸药选用高密度高爆速炸药,药量选用1.5kg~2kg;
其中,采用锤击多次叠加小折射法对高速层埋藏深度进行探测,锤击小折射法的观测系统参数为:
仪器使用工程地震仪;采样间隔0.25ms;记录长度0.5s;记录格式为SEG-2;激发源为20磅大锤,采用大于10次的能量叠加的方式激发;观测方式为相遇观测;接收道数24道;检波器使用主频为10Hz~40Hz的单只纵波检波器;最小偏移距2~6m;24道检波器之间的距离设置为:1m,1m,2m,2m,3m,5m,5m,8m,8m,10m,10m,10m,10m,10m,8m,8m,5m,5m,3m,2m,2m,1m,1m;使用截距时间法对采集的数据进行处理和解释;
②步骤(2)中,地震信号激发震源采用可控机械震源时:
反射波信噪比≥1时,可控震源台次使用1台*2次至1台*4次;
反射波信噪比<1时,可控震源台次使用2台*2次;
扫描频率使用宽频激发,低频低于6Hz,高频大于100Hz,扫描时间8s~14s,线性扫频激发;
步骤(3)接收技术中,设置砂岩型铀矿地震勘探的接收参数如下:
①选用主频≤10Hz的压电式加速度检波器,采用4~10只检波器点式组合的方式进行地震数据接收,并埋置于地下20~30cm的硬实土壤中;
②选用真24位地震数据采集系统进行地震信号接收,仪器的采样间隔为0.5ms‐1ms,记录长度为3s‐4s,前放增益为0‐12dB。
(3)接收技术
设置砂岩型铀矿地震勘探的检波器和其他仪器的接收参数;
(4)部署检波器检波器、激发震源
根据步骤(1)中设计的观测参数部署检波器、激发震源的位置;
砂岩型铀矿三维地震观测系统部署原则是:
(a)检波器的布设方向垂直于地质构造延伸方向,激发震源的布设方向平行于构造延伸方向;
(b)检波器的布设方向和激发震源点的布设方向是垂直的;
(5)激发和接收地震信号
根据步骤(2)中确定的地震信号激发源的参数,使用步骤(2)中确定的地震信号激发源激发地震信号;
根据步骤(3)中确定的检波器和其他仪器的接收参数,启动地震数据采集系统接收信号。
最大炮检距在800-940m之间,覆盖次数在64~128之间,接收线距在60~80m之间。接收线数为16线,单线接收道数为128道,面元为5m*10m,覆盖次数为128次,检波点距为10m,炮点距为20m,接收线距为80m,炮线距为40m,纵向最大偏移距为635m,最大横向偏移距为630m,最大最小偏移距为89m,最大炮检距为894m,横向滚动距离为80m,横纵比为0.992。或者接收线数为16线,单线接收道数为64道,面元为10m*5m,覆盖次数为128次,检波点距为20m,炮点距为20m,接收线距为80m,炮线距为40m,纵向最大偏移距为630m,最大横向偏移距为635m,最大最小偏移距为89m,最大炮检距为894m,横向滚动距离为80m,横纵比为0.992。
Claims (8)
1.一种砂岩型铀矿三维地震数据观测方法,其特征在于,包括以下步骤:
(1)确定观测系统参数
需要确定的砂岩型铀矿三维地震数据采集观测系统的关键参数包括面元尺寸、覆盖次数、最大炮检距、方位参数、横向滚动距、接收线距、最大非纵距、最大偏移距;
(2)确定地震信号激发技术
地震信号激发技术主要是设置砂岩型铀矿地震勘探中的震源激发参数,涉及两种地震信号激发震源:炸药震源和可控机械震源;
(3)接收技术
设置砂岩型铀矿地震勘探的检波器和其他仪器的接收参数;
(4)部署检波器检波器、激发震源
根据步骤(1)中设计的观测参数部署检波器、激发震源的位置;
砂岩型铀矿三维地震观测系统部署原则是:
(a)检波器的布设方向垂直于地质构造延伸方向,激发震源的布设方向平行于构造延伸方向;
(b)检波器的布设方向和激发震源点的布设方向是垂直的;
(5)激发和接收地震信号
根据步骤(2)中确定的地震信号激发源的参数,使用步骤(2)中确定的地震信号激发源激发地震信号;
根据步骤(3)中确定的检波器和其他仪器的接收参数,启动地震数据采集系统接收信号。
2.如权利要求1所述的一种砂岩型铀矿三维地震数据观测方法,其特征在于:
步骤(1)中,面元尺寸的确定方法如下:
A、满足横向分辨率的要求
定义在地震信号谱密度曲线上与最大值对应的频率为地震信号优势频率;在地震信号优势频率的波长内设置2-3个采样点,从而控制对应的面元满足系统需要的分辨率,具体公式如下:
x=Vint/(2*Fp)
式中:x-确定所求面元尺寸的第一上限值;Vint-目的层上覆地层层速度;Fp-目的层反射波主频,表示由由不同频率的波动叠合在一起构成的地震波中能量最强的波动所对应的频率成份;
公式中Vint、Fp根据工区地质资料确定;
B、满足最高无混叠频率的要求
Fmax表示每个倾斜地震反射波同相轴的一个偏移前可能的最高无混叠频率,它根据目的层反射波均方根速度Vrms、倾角θ、面元尺寸确定,具体公式如下:
y=Vrms/(4Fmax*Sinθ)
式中:y-确定所求面元尺寸的第二上限值;Vrms-目的层反射波均方根速度,Fmax-最高无混叠频率,θ-目的层倾角,表示每个倾斜地震反射波同相轴与水平方向的角度;
公式中Vrms、Fmax、θ根据工区地质资料确定;
C、通过上述A、B两步计算得到面元尺寸b的取值上限:
若x≤y,则b≤x;
若y≤x,则b≤y;
D、确定面元尺寸b的取值下限
确定面元尺寸b的取值下限为5m;
E、确定面元尺寸b的取值范围
若y≤x,5m≤b≤y;
若x≤y,5m≤b≤x。
F、在面元尺寸b的取值范围内选取数值确定面元尺寸。
3.如权利要求2所述的一种砂岩型铀矿三维地震数据观测方法,其特征在于:砂岩型铀矿地震勘探中的面元尺寸选择为以下几种之一:
面元尺寸是一组两个数,第一个数的含义是面元的纵向距离,其数值的两倍对应检波点距;第二个数的含义是面元的横向距离,其数值的两倍对应炮点距;
5m*5m,对应10m检波点距*10m炮点距;
7.5m*7.5m,对应15m检波点距*15m炮点距;
5*10m,对应10m检波点距*20m炮点距;
10*5m,对应20m检波点距*10m炮点距。
4.如权利要求1所述的一种砂岩型铀矿三维地震数据观测方法,其特征在于:
(一)步骤(1)中,覆盖次数的确定方法如下:
首先使用克朗软件噪声分析模块估算出以往地震剖面中目的层信号的信噪比(S/N)raw,再根据新资料期望达到的(S/N)required,按照以下公式计算覆盖次数Foldrequired:
<mrow>
<msub>
<mi>Fold</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>q</mi>
<mi>u</mi>
<mi>i</mi>
<mi>r</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>/</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>q</mi>
<mi>u</mi>
<mi>i</mi>
<mi>r</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>/</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>w</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&times;</mo>
<msub>
<mi>Fold</mi>
<mrow>
<mi>o</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
公式中,Foldold表示以往地震资料的覆盖次数;
根据要解决的技术问题限定在砂岩型铀矿地震勘探中新资料期望达到的(S/N)required:
如要解决地下地质体的构造问题,则(S/N)required≥1;
如要解决砂岩型铀矿勘查中砂体预测问题,则(S/N)required≥2;
(二)步骤(1)中,最大炮检距的确定方法如下:
最大炮检距用X表示,根据a、b、c三个条件确定其取值范围;
a.最大炮检距X在0.9倍至1.2倍主要目标深度之间;
b.动校拉伸率D与最大炮检距X的关系如下:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>T</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中:T0-目的层双程反射时间,V-目标层地震均方根速度;
根据最大炮检距X满足动校拉伸率D小于12.5%的要求,确定最大炮检距X≥5VT0;
c.最大炮检距满足速度分析要求,根据如下公式确定:
<mrow>
<mi>X</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mi>O</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mi>P</mi>
</msub>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</msqrt>
</mrow>
2
式中:P-速度分析精度,根据系统设定确定,P<6%;
(三)步骤(1)中,方位角参数的确定方法如下:
方位角参数=横向最大偏移距/纵向最大偏移距,在砂岩型铀矿地震勘探领域方位角参数取值在0.9-1之间;
横向最大偏移距为地震勘探的常用参数,横向最大偏移距=((接收线距*(接收线数-1)/2)+1/2*接收线距-1/2*炮点距);
纵向最大偏移距为地震勘探的常用参数,纵向最大偏移距=(检波点距*(单线接收道数-1)/2);
综上,方位角=((接收线距*(接收线数-1)/2)+1/2*接收线距-1/2*炮点距)/(检波点距*(单线接收道数-1)/2);
根据上述公式,通过设置接收线距、接收线数、炮点距、检波点距、单线接收道数的取值,使得方位角参数落入0.9-1范围;
(四)接收线距的确定方法如下:
接收线距≤垂直入射时的菲涅尔带半径:
<mrow>
<mi>R</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>V</mi>
<mrow>
<mn>4</mn>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
式中:R-菲涅尔带半径,根据工区地质资料确定;
(五)横向滚动距的确定方法如下:
横向滚动距=接收线距;
(六)最大横向偏移距的确定方法如下:
最大横向偏移距Ymax满足:
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mfrac>
<msub>
<mi>V</mi>
<mi>a</mi>
</msub>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
</mrow>
</mfrac>
<msqrt>
<mfrac>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
式中:Va-目标层反射波平均速度,根据工区地质资料确定;
(七)最大最小偏移距的确定方法如下:
最浅目标层的埋深Xm根据工区地质资料确定,最大最小偏移距小于Xm。
5.如权利要求1所述的一种砂岩型铀矿三维地震数据观测方法,其特征在于:
最大炮检距在800-940m之间,覆盖次数在64~128之间,接收线距在60~80m之间。
6.如权利要求1所述的一种砂岩型铀矿三维地震数据观测方法,其特征在于:
接收线数为16线,单线接收道数为128道,面元为5m*10m,覆盖次数为128次,检波点距为10m,炮点距为20m,接收线距为80m,炮线距为40m,纵向最大偏移距为635m,最大横向偏移距为630m,最大最小偏移距为89m,最大炮检距为894m,横向滚动距离为80m,横纵比为0.992。
7.如权利要求1所述的一种砂岩型铀矿三维地震数据观测方法,其特征在于:
接收线数为16线,单线接收道数为64道,面元为10m*5m,覆盖次数为128次,检波点距为20m,炮点距为20m,接收线距为80m,炮线距为40m,纵向最大偏移距为630m,最大横向偏移距为635m,最大最小偏移距为89m,最大炮检距为894m,横向滚动距离为80m,横纵比为0.992。
8.如权利要求1所述的一种砂岩型铀矿三维地震数据观测方法,其特征在于:
①步骤(2)中,地震信号激发震源采用炸药震源时:
激发炮井的深度选择在高速层的埋藏深度以下3m,炸药选用高密度高爆速炸药,药量选用1.5kg~2kg;
其中,采用锤击多次叠加小折射法对高速层埋藏深度进行探测,锤击小折射法的观测系统参数为:
仪器使用工程地震仪;采样间隔0.25ms;记录长度0.5s;记录格式为SEG-2;激发源为20磅大锤,采用大于10次的能量叠加的方式激发;观测方式为相遇观测;接收道数24道;检波器使用主频为10Hz~40Hz的单只纵波检波器;最小偏移距2~6m;24道检波器之间的距离设置为:1m,1m,2m,2m,3m,5m,5m,8m,8m,10m,10m,10m,10m,10m,8m,8m,5m,5m,3m,2m,2m,1m,1m;使用截距时间法对采集的数据进行处理和解释;
②步骤(2)中,地震信号激发震源采用可控机械震源时:
反射波信噪比≥1时,可控震源台次使用1台*2次至1台*4次;
反射波信噪比<1时,可控震源台次使用2台*2次;
扫描频率使用宽频激发,低频低于6Hz,高频大于100Hz,扫描时间8s~14s,线性扫频激发;
步骤(3)接收技术中,设置砂岩型铀矿地震勘探的接收参数如下:
①选用主频≤10Hz的压电式加速度检波器,采用4~10只检波器点式组合的方式进行地震数据接收,并埋置于地下20~30cm的硬实土壤中;
②选用真24位地震数据采集系统进行地震信号接收,仪器的采样间隔为0.5ms‐1ms,记录长度为3s‐4s,前放增益为0‐12dB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710236881.0A CN107144873A (zh) | 2017-04-12 | 2017-04-12 | 一种砂岩型铀矿三维地震数据观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710236881.0A CN107144873A (zh) | 2017-04-12 | 2017-04-12 | 一种砂岩型铀矿三维地震数据观测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107144873A true CN107144873A (zh) | 2017-09-08 |
Family
ID=59773829
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710236881.0A Pending CN107144873A (zh) | 2017-04-12 | 2017-04-12 | 一种砂岩型铀矿三维地震数据观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107144873A (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108897041A (zh) * | 2018-08-16 | 2018-11-27 | 中国石油天然气股份有限公司 | 一种铀矿富集区的预测方法和装置 |
CN111190240A (zh) * | 2020-01-13 | 2020-05-22 | 吉林大学 | 基于三维地震解释的砂岩型铀矿成矿构造要素的提取方法 |
CN111239815A (zh) * | 2020-01-20 | 2020-06-05 | 吉林大学 | 基于三维地震属性的砂岩型铀储层成矿沉积要素提取方法 |
CN111239801A (zh) * | 2020-04-03 | 2020-06-05 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据采集方法、装置、设备及存储介质 |
CN112255672A (zh) * | 2020-09-24 | 2021-01-22 | 中国石油天然气股份有限公司 | 一种地震数据采集覆盖次数优化方法 |
CN112346114A (zh) * | 2020-10-29 | 2021-02-09 | 中国石油天然气集团有限公司 | 油气地震数据采集接收道数的确定方法及装置 |
CN112764093A (zh) * | 2019-10-21 | 2021-05-07 | 中国石油天然气集团有限公司 | 三维地震数据采集方法及系统 |
CN112782768A (zh) * | 2019-11-04 | 2021-05-11 | 中国石油天然气集团有限公司 | 一种地震勘探激发参数的试验方法及装置 |
CN113109864A (zh) * | 2020-01-09 | 2021-07-13 | 中国石油天然气集团有限公司 | 地震勘探炮点变观方法及装置 |
CN113156512A (zh) * | 2020-01-07 | 2021-07-23 | 中国石油天然气集团有限公司 | 适合可控震源激发的三维地震观测方法、设备及系统 |
CN113514886A (zh) * | 2021-07-22 | 2021-10-19 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法 |
CN113687410A (zh) * | 2021-08-19 | 2021-11-23 | 吉林大学 | 一种金属矿地震勘探数据采集方法 |
WO2023000257A1 (zh) * | 2021-07-22 | 2023-01-26 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050157589A1 (en) * | 2004-01-20 | 2005-07-21 | Andreas Laake | Survey design using earth observation data |
CN102213770A (zh) * | 2010-04-01 | 2011-10-12 | 罗仁泽 | 一种地震资料低信噪比地区组合面积激发接收方法 |
CN103901465A (zh) * | 2014-03-21 | 2014-07-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 全息三维地震勘探观测系统的设计方法 |
CN104142514A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种三维地震观测系统定量设计方法 |
CN104536041A (zh) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | 一种地震观测系统参数的优化方法 |
CN104932005A (zh) * | 2015-06-19 | 2015-09-23 | 中石化石油工程技术服务有限公司 | 一种三维观测系统布设的方法及对应的三维观测系统 |
CN105223630A (zh) * | 2014-06-26 | 2016-01-06 | 中石化石油工程地球物理有限公司胜利分公司 | 基于地质模型的全方位观测系统参数论证方法 |
-
2017
- 2017-04-12 CN CN201710236881.0A patent/CN107144873A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050157589A1 (en) * | 2004-01-20 | 2005-07-21 | Andreas Laake | Survey design using earth observation data |
CN102213770A (zh) * | 2010-04-01 | 2011-10-12 | 罗仁泽 | 一种地震资料低信噪比地区组合面积激发接收方法 |
CN104142514A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种三维地震观测系统定量设计方法 |
CN103901465A (zh) * | 2014-03-21 | 2014-07-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 全息三维地震勘探观测系统的设计方法 |
CN105223630A (zh) * | 2014-06-26 | 2016-01-06 | 中石化石油工程地球物理有限公司胜利分公司 | 基于地质模型的全方位观测系统参数论证方法 |
CN104536041A (zh) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | 一种地震观测系统参数的优化方法 |
CN104932005A (zh) * | 2015-06-19 | 2015-09-23 | 中石化石油工程技术服务有限公司 | 一种三维观测系统布设的方法及对应的三维观测系统 |
Non-Patent Citations (1)
Title |
---|
陆基孟等: "《地震勘探原理》", 28 February 2009, 中国石油大学出版社 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108897041A (zh) * | 2018-08-16 | 2018-11-27 | 中国石油天然气股份有限公司 | 一种铀矿富集区的预测方法和装置 |
CN112764093B (zh) * | 2019-10-21 | 2023-08-22 | 中国石油天然气集团有限公司 | 三维地震数据采集方法及系统 |
CN112764093A (zh) * | 2019-10-21 | 2021-05-07 | 中国石油天然气集团有限公司 | 三维地震数据采集方法及系统 |
CN112782768A (zh) * | 2019-11-04 | 2021-05-11 | 中国石油天然气集团有限公司 | 一种地震勘探激发参数的试验方法及装置 |
CN113156512A (zh) * | 2020-01-07 | 2021-07-23 | 中国石油天然气集团有限公司 | 适合可控震源激发的三维地震观测方法、设备及系统 |
CN113109864B (zh) * | 2020-01-09 | 2023-07-25 | 中国石油天然气集团有限公司 | 地震勘探炮点变观方法及装置 |
CN113109864A (zh) * | 2020-01-09 | 2021-07-13 | 中国石油天然气集团有限公司 | 地震勘探炮点变观方法及装置 |
CN111190240A (zh) * | 2020-01-13 | 2020-05-22 | 吉林大学 | 基于三维地震解释的砂岩型铀矿成矿构造要素的提取方法 |
CN111190240B (zh) * | 2020-01-13 | 2021-09-07 | 吉林大学 | 基于三维地震解释的砂岩型铀矿成矿构造要素的提取方法 |
CN111239815A (zh) * | 2020-01-20 | 2020-06-05 | 吉林大学 | 基于三维地震属性的砂岩型铀储层成矿沉积要素提取方法 |
CN111239815B (zh) * | 2020-01-20 | 2021-07-13 | 吉林大学 | 基于三维地震属性的砂岩型铀储层成矿沉积要素提取方法 |
CN111239801A (zh) * | 2020-04-03 | 2020-06-05 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据采集方法、装置、设备及存储介质 |
CN112255672A (zh) * | 2020-09-24 | 2021-01-22 | 中国石油天然气股份有限公司 | 一种地震数据采集覆盖次数优化方法 |
CN112346114A (zh) * | 2020-10-29 | 2021-02-09 | 中国石油天然气集团有限公司 | 油气地震数据采集接收道数的确定方法及装置 |
CN113514886B (zh) * | 2021-07-22 | 2021-12-10 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法 |
WO2023000257A1 (zh) * | 2021-07-22 | 2023-01-26 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法 |
CN113514886A (zh) * | 2021-07-22 | 2021-10-19 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法 |
CN113687410A (zh) * | 2021-08-19 | 2021-11-23 | 吉林大学 | 一种金属矿地震勘探数据采集方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107144873A (zh) | 一种砂岩型铀矿三维地震数据观测方法 | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
Sloan et al. | Detecting clandestine tunnels using near-surface seismic techniques | |
Beres et al. | Mapping the architecture of glaciofluvial sediments with three-dimensional georadar | |
Malehmir et al. | Seismic methods in mineral exploration and mine planning: A general overview of past and present case histories and a look into the future | |
She et al. | Shallow crustal structure of the middle‐lower Yangtze River region in eastern China from surface‐wave tomography of a large volume airgun‐shot experiment | |
CN102288990B (zh) | 超浅层三维地震反射波勘探方法 | |
CN102495434A (zh) | 地下工程超前地质预报的方法 | |
CN102455439A (zh) | 基于克希霍夫积分法的绕射波场分离方法 | |
Steer et al. | Deep structure of the continental lithosphere in an unextended orogen: An explosive‐source seismic reflection profile in the Urals (Urals Seismic Experiment and Integrated Studies (URSEIS 1995)) | |
CN102062869B (zh) | 一种地震勘探观测系统设计方法 | |
Ridwan et al. | Development of an engineering bedrock map beneath Jakarta based on microtremor array measurements | |
CN102053261A (zh) | 一种地震数据处理方法 | |
Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
CN107356668A (zh) | 一种坝体安全鉴定的检测方法 | |
Baglari et al. | A state-of-the-art review of passive MASW survey for subsurface profiling | |
CN106324682A (zh) | 一种应用于冻土层地区的表层结构调查方法 | |
Igboekwe et al. | Investigation into the weathering layer using uphole method of seismic refraction | |
Arosio et al. | Seismic characterization of an abandoned mine site | |
CN102520441B (zh) | 一种井巷定向探测地震勘探数据采集方法 | |
Cunningham et al. | Interplay of late Cenozoic siliciclastic supply and carbonate response on the southeast Florida platform | |
Gorstein et al. | Combination of HVSR and MASW methods to obtain shear wave velocity model of subsurface in Israel | |
Maunde et al. | Seismic refraction investigation of fracture zones and bedrock configuration for geohydrologic and geotechnical studies in part of Nigeria’s Capital City, Abuja | |
Mohamed Adel et al. | Definition of soil characteristics and ground response at the northwestern part of the Gulf of Suez, Egypt | |
Carvalho et al. | Characterization of an intraplate seismogenic zone using geophysical and borehole data: the Vila Franca de Xira Fault, Portugal |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170908 |