CN102759749B - 一种瑞雷面波速度分析方法 - Google Patents
一种瑞雷面波速度分析方法 Download PDFInfo
- Publication number
- CN102759749B CN102759749B CN201210231886.1A CN201210231886A CN102759749B CN 102759749 B CN102759749 B CN 102759749B CN 201210231886 A CN201210231886 A CN 201210231886A CN 102759749 B CN102759749 B CN 102759749B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- frequency
- rayleigh surface
- surface wave
- penetration depth
- 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.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种瑞雷面波速度分析方法,步骤为:一,将含有m道,每道有n个采样点,时间采样率为的单炮瑞雷面波地震记录读取到二维数组g中;二,对g进行二维傅里叶变换,得到f-k谱;三,建立穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系;四,依据瑞雷面波穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,将f-k谱中对应(k,f)坐标的能量映射到h-v谱中对应(v,h)坐标的位置,选定h-v谱输出的坐标范围并成图;五,基于h-v谱,从浅至深依次拾取h-v谱中的强能量团,并提取对应的(v,h)坐标值,信息提取简单、方便和精度较高的优点。
Description
技术领域
本发明属于地震瑞雷面波速度参数提取技术领域,特别涉及一种瑞雷面波速度分析方法。
背景技术
瑞雷面波在研究地球内部结构、近地表各种工程勘查问题中,尤其在近地表地层结构与属性特征获取方面,是一种非常有效的地震勘探方法。由于瑞雷面波存在频散特性,现有方法一方面主要是基于f-k谱来提取瑞雷面波频散曲线,根据频散曲线特性进而可提取相应的地层结构、构造等信息;同时,由于瑞雷面波与横波速度密切相关,因此,处理瑞雷面波的方法另一方面是通过提取瑞雷面波的频散曲线来反演横波速度,进而获取地层结构、构造等信息来达到勘探的目的。其实,瑞雷面波中隐含的信息是非常丰富的,现有方法存在两方面的缺陷:一方面是处理后提取的信息有限,体现的地质信息过于单一、不够直观,直接解决实际地质问题的能力有限;另一方面,瑞雷面波穿透深度概念比较模糊,未能直接建立瑞雷面波隐含的地质信息与深度的关系,从而影响瑞雷面波解决实际地质问题的能力。
发明内容
为了克服上述现有技术的不足,本发明的目的在于提出了一种瑞雷面波速度分析方法,该方法具有瑞雷面波速度信息提取简单、方便,瑞雷面波速度与深度对应直观、明了,直接解决实际地质问题和精度较高的优点。
为了实现上述步骤,本发明采用的技术方案是:一种瑞雷面波速度分析方法,包括如下步骤:
第一步,将含有m道,每道有n个采样点,时间采样率为 的单炮瑞雷面波地震记录读取到二维数组g中;
第二步,对g进行二维傅里叶变换,得到f-k谱,二维离散傅里叶正变换为:
其中,k为波数,为波数采样率,j 0=1,2,…,m,为波数采样序号;f为频率,为频率采样率,l 0=1,2,…,n,为频率采样序号;j=1,2,…,m,为道序号;l=1,2,…,n,为采样点序号;为道间距,为最小偏移距,x为偏移距;为时间采样率,t为采样时间; 为虚数单位;
第三步,建立穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,具体如下:
瑞雷面波的穿透深度(h)与波长(λ)和泊松比(σ)有关,其关系式为:
则穿透深度h与速度v的关系为:
其中,C为与泊松比有关的深度校正因子,通常,近地表地层的泊松比范围在0.2~0.35之间,因此,深度校正因子C取值范围一般为0.625~0.75;
第四步,依据瑞雷面波穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,将f-k谱中对应(k,f)坐标的能量映射到h-v谱中对应(v,h)坐标的位置,同时选定h-v谱输出的坐标范围并成图;
第五步,基于h-v谱,从浅至深依次拾取h-v谱中的强能量团,并提取对应的(v,h)坐标值,便实现了瑞雷面波速度分析。
本发明的有益效果是:
与现有技术相比,本发明通过二维傅立叶变换获取瑞雷面波的f-k谱,依据瑞雷面波穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,将f-k谱转换成深度-速度(h-v)谱,并进行瑞雷面波速度分析,依据瑞雷面波穿透深度(h)与速度(v)的关系获取地层结构、构造等信息,直接解决实际地质问题,同时可获取瑞雷面波速度,提供其它处理技术参数。
附图说明
图1为本发明实际瑞雷面波单炮地震记录及与之对应的f-k谱,其中图1(a)为瑞雷面波单炮地震记录;图1(b)为瑞雷面波的f-k谱。
图2为本发明泊松比σ与深度校正因子C之间的关系曲线图。
图3为本发明实际瑞雷面波地震记录的h-v谱及提取的深度-速度(h-v)曲线,其中图3(a)为瑞雷面波地震记录的h-v谱;图3(b)为提取的深度-速度(h-v)曲线。
具体实施方法
下面结合附图对本发明进一步详细说明。
一种瑞雷面波速度分析方法,包括如下步骤:
第一步,将含有m道,每道有n个采样点,时间采样率为的单炮瑞雷面波地震记录读取到二维数组g中,参见表1、图1、图2、图3;
第二步,对g进行二维傅里叶变换,得到f-k谱,二维离散傅里叶正变换为:
其中,k为波数,为波数采样率,j 0=1,2,…,m,为波数采样序号;f为频率,为频率采样率,l 0=1,2,…,n,为频率采样序号;j=1,2,…,m,为道序号;l=1,2,…,n,为采样点序号;为道间距,为最小偏移距,x为偏移距;为时间采样率,t为采样时间; 为虚数单位;
第三步,建立穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,具体如下:
瑞雷面波的穿透深度(h)与波长(λ)和泊松比(σ)有关,其关系式为:
则穿透深度h与速度v的关系为:
其中,C为与泊松比有关的深度校正因子,通常,近地表地层的泊松比范围在0.2~0.35之间,因此,深度校正因子C取值范围一般为0.625~0.75,参见表1, 图2;
第四步,依据瑞雷面波穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,将f-k谱中对应(k,f)坐标的能量映射到h-v谱中对应(v,h)坐标的位置,同时选定h-v谱输出的坐标范围并成图,参见图1、图2、图3(a);
第五步,基于h-v谱,从浅至深依次拾取h-v谱中的强能量团,并提取对应的(v,h)坐标值,便实现了瑞雷面波速度分析,参见图3。
表1为本发明泊松比σ、波长λ与瑞雷面波穿透深度h之间的关系。当地层的泊松比σ为0.1时,瑞雷面波穿透深度h是其波长λ的0.55倍;当地层的泊松比σ为0.15时,瑞雷面波穿透深度h是其波长λ的0.575倍;当地层的泊松比σ为0.20时,瑞雷面波穿透深度h是其波长λ的0.625倍,等等。深度校正因子C的取值需根据实际情况而定,由于近地表地层的泊松比σ范围在0.20~0.35之间,因此,深度校正因子C取值范围一般为0.625~0.75,即穿透深度h是波长λ的0.625~0.75。
图1为本发明实际瑞雷面波单炮地震记录及与之对应的f-k谱。其中图1(a)为瑞雷面波单炮地震记录,横坐标为道号,纵坐标为瑞雷面波旅行时,也称为采样时间,图1(a)表明瑞雷面波的波场特征频散现象明显,能量较强、成扫把状分布;图1(b)为瑞雷面波的f-k谱,横坐标为波数k,纵坐标为频率f,图1(b)表明波谱特征明显,能量较强,且集中,频率f范围为4.2-9Hz,波数k范围为0.009-0.04 m-1,但是该谱与其它物性参数之间的关系表现不直接。
图2为本发明泊松比σ与深度校正因子C之间的关系曲线图,横坐标为泊松比σ,纵坐标为深度校正因子C,图中标号“+”为与泊松比σ值相对应的深度校正因子C值。
图3为本发明实际瑞雷面波地震记录的h-v谱及提取的深度-速度(h-v)曲线。其中图3(a)为瑞雷面波地震记录的h-v谱,横坐标为速度v,纵坐标为深度h,图3(a)建立了速度v与深度h的直接关系,波谱特征明显,能量集中,存在速度v随深度h增加而增大的特征,且瑞雷面波穿透深度h范围为16-76m,速度v范围为220-540m/s,在26-29m存在一物性相对均匀的地层;图3(b)为提取的深度-速度(h-v)曲线,横坐标为速度v,纵坐标为深度h,图中标号“○”为在h-v谱中拾取同一深度强能量团对应的(v,h)坐标值,图3(b)表明了速度v随深度h的衰减趋势和瑞雷面波穿透深度h与速度v的对应关系。
实施例
将一含有100道,每道1000个采样点,时间采样率为4ms的实际单炮面波地震记录为例说明本实例的实施步骤:
第一步,将含有100道,每道有1000个采样点,时间采样率为4ms的单炮瑞雷面波地震记录读取到二维数组g中;
第二步,对g进行二维傅里叶变换,得到f-k谱,二维离散傅里叶正变换为:
其中,k为波数,为波数采样率,j 0=1,2,…,100,为波数采样序号;f为频率,为频率采样率,l 0=1,2,…,1000,为频率采样序号;j=1,2,…,100,为道序号;l=1,2,…,1000,为采样点序号;为道间距,为最小偏移距,x为偏移距;为时间采样率,t为采样时间; 为虚数单位;
第三步,建立穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,具体如下:
取与泊松比有关的深度校正因子C为0.65,则其关系式为:
则穿透深度(h)与频率(f)的关系为:
第四步,依据瑞雷面波穿透深度(h)、速度(v)与频率(f)、波数(k)、波长(λ)、泊松比(σ)之间的关系,将f-k谱中对应(k,f)坐标的能量映射到h-v谱中对应(v,h)坐标的位置,同时选定h-v谱输出的坐标范围并成图,其中穿透深度(h)的范围为:0~-100m(取地层沿地表向下为负),速度(v)的范围为:200~800m/s;
第五步,基于h-v谱,从浅至深依次拾取h-v谱中的强能量团,并提取对应的(v,h)坐标值,便实现了瑞雷面波速度分析。
Claims (2)
1.一种瑞雷面波速度分析方法,其特征在于,包括如下步骤:
第一步,将含有m道,每道有n个采样点,时间采样率为Δt的单炮瑞雷面波地震记录读取到二维数组g中;
第二步,对g进行二维傅里叶变换,得到f-k谱,二维离散傅里叶正变换为:
k=j0·Δk
f=l0·Δf
x=x0+(j-1)·Δx
t=l·Δt
其中,k为波数,Δk为波数采样率,j0=1,2,…,m,为波数采样序号;f为频率,Δf为频率采样率,l0=1,2,…,n,为频率采样序号;j=1,2,…,m,为道序号;l=1,2,…,n,为采样点序号;Δx为道间距,x0为最小偏移距,x为偏移距;Δt为时间采样率,t为采样时间;为虚数单位;
第三步,建立穿透深度h、速度v与频率f、波数k、波长λ、泊松比σ之间的关系,具体如下:
瑞雷面波的穿透深度h与波长λ和泊松比σ有关,其关系式为:
h=Cλ
则穿透深度h与速度v的关系为:
其中,C为与泊松比有关的深度校正因子,通常,近地表地层的泊松比范围在0.2~0.35之间,因此,深度校正因子取值范围一般为0.625~0.75;
第四步,依据瑞雷面波穿透深度h、速度v与频率f、波数k、波长λ、泊松比σ之间的关系,将f-k谱中对应(k,f)坐标的能量映射到h-v谱中对应(v,h)坐标的位置,同时选定h-v谱输出的坐标范围并成图;
第五步,基于h-v谱,从浅至深依次拾取h-v谱中的强能量团,并提取对应的(v,h)坐标值,便实现了瑞雷面波速度分析。
2.根据权利要求1所述的一种瑞雷面波速度分析方法,其特征在于,包括如下步骤:
将一含有100道,每道1000个采样点,时间采样率为4ms的实际单炮面波地震记录读取到二维数组g中;
第二步,对g进行二维傅里叶变换,得到f-k谱,二维离散傅里叶正变换为:
k=j0·Δk
f=l0·Δf
x=x0+(j-1)·Δx
t=l·Δt
其中,k为波数,Δk为波数采样率,j0=1,2,…,100,为波数采样序号;f为频率,Δf为频率采样率,l0=1,2,…,1000,为频率采样序号;j=1,2,…,100,为道序号;l=1,2,…,1000,为采样点序号;Δx为道间距,x0为最小偏移距,x为偏移距;Δt为时间采样率,t为采样时间;为虚数单位;
第三步,建立穿透深度h、速度v与频率f、波数k、波长λ、泊松比σ之间的关系,具体如下:
取与泊松比有关的深度校正因子C为0.65,则其关系式为:
h=0.65λ
则穿透深度h与频率f的关系为:
第四步,依据瑞雷面波穿透深度h、速度v与频率f、波数k、波长λ、泊松比σ之间的关系,将f-k谱中对应(k,f)坐标的能量映射到h-v谱中对应(v,h)坐标的位置,同时选定h-v谱输出的坐标范围并成图,其中穿透深度h的范围为:0~-100m取地层沿地表向下为负,速度v的范围为:200~800m/s;
第五步,基于h-v谱,从浅至深依次拾取h-v谱中的强能量团,并提取对应的(v,h)坐标值,便实现了瑞雷面波速度分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210231886.1A CN102759749B (zh) | 2012-07-06 | 2012-07-06 | 一种瑞雷面波速度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210231886.1A CN102759749B (zh) | 2012-07-06 | 2012-07-06 | 一种瑞雷面波速度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102759749A CN102759749A (zh) | 2012-10-31 |
CN102759749B true CN102759749B (zh) | 2014-11-19 |
Family
ID=47054258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210231886.1A Expired - Fee Related CN102759749B (zh) | 2012-07-06 | 2012-07-06 | 一种瑞雷面波速度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102759749B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108572390B (zh) * | 2018-05-22 | 2019-10-25 | 武汉市市政建设集团有限公司 | 利用表面波谱扰动预测浅部洞穴的探测方法 |
CN114415234B (zh) * | 2022-01-24 | 2023-05-12 | 西南交通大学 | 基于主动源面波频散和h/v确定浅地表横波速度的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6151275A (en) * | 1998-09-11 | 2000-11-21 | Pgs Tensor, Inc. | Method of dual wavefield reinforcement |
CN101907727B (zh) * | 2010-08-17 | 2012-05-30 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
-
2012
- 2012-07-06 CN CN201210231886.1A patent/CN102759749B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6151275A (en) * | 1998-09-11 | 2000-11-21 | Pgs Tensor, Inc. | Method of dual wavefield reinforcement |
CN101907727B (zh) * | 2010-08-17 | 2012-05-30 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
Non-Patent Citations (2)
Title |
---|
杜正涛等.瑞雷面波勘探技术在第四系分层方面的应用.《物探与化探》.1999,第23卷(第4期),277-282. * |
瑞雷面波勘探技术在第四系分层方面的应用;杜正涛等;《物探与化探》;19990831;第23卷(第4期);277-282 * |
Also Published As
Publication number | Publication date |
---|---|
CN102759749A (zh) | 2012-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jin et al. | The global boundary stratotype section and point (GSSP) for the base of Changhsingian stage (Upper Permian) | |
CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
CN100429529C (zh) | 利用全波列、偶极横波测井资料确定气层的方法 | |
CN103399346B (zh) | 一种井震联合初始波阻抗建模方法 | |
CN103424777B (zh) | 一种提高地震成像分辨率的方法 | |
CN103698807B (zh) | 梯状二维宽线观测系统设计方法 | |
CN104330827B (zh) | 一种地表模型静校正的处理方法 | |
CN106646601A (zh) | 多信息联合约束的浅中深层三维q体建立方法 | |
CN104865614A (zh) | 基于可变骨架参数的复杂储层流体识别方法 | |
CN107255831A (zh) | 一种叠前频散属性的提取方法 | |
CN104280772A (zh) | 一种井中微地震震相识别方法 | |
CN110824563A (zh) | 一种基于Xgboost算法的储层岩性预测方法 | |
CN107356965A (zh) | 基于加权叠加噪音压制策略的反射系数反演储层预测方法 | |
CN102749646B (zh) | 一种瑞雷面波深度-频率分析方法 | |
CN102759749B (zh) | 一种瑞雷面波速度分析方法 | |
CN102866425A (zh) | 一种基于盲源地震信号平稳叠加模型的盲分离方法 | |
CN106094044A (zh) | 一种多道瞬变电磁法(mtem)虚拟波场提取装置与方法 | |
CN105093318A (zh) | 一种自适应波动方程波场延拓静校正方法 | |
CN105866837A (zh) | 一种海上非重复性时移地震多层系互均衡处理方法 | |
CN105676292A (zh) | 一种基于二维曲波变换的三维地震数据去噪方法 | |
CN108375789A (zh) | 联合采集地震数据的同步匹配方法 | |
CN107356963A (zh) | 一种数据驱动的自适应的地震信号相干体属性分析方法 | |
CN102778691B (zh) | 一种计算检波器组内静校正时差的方法 | |
CN104459789A (zh) | 一种地震波初至拾取的方法 | |
CN104297798A (zh) | 巨厚黄土塬区多井微测井调查方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141119 Termination date: 20150706 |
|
EXPY | Termination of patent right or utility model |