CN104678435A - 一种提取Rayleigh面波频散曲线的方法 - Google Patents
一种提取Rayleigh面波频散曲线的方法 Download PDFInfo
- Publication number
- CN104678435A CN104678435A CN201410580207.0A CN201410580207A CN104678435A CN 104678435 A CN104678435 A CN 104678435A CN 201410580207 A CN201410580207 A CN 201410580207A CN 104678435 A CN104678435 A CN 104678435A
- Authority
- CN
- China
- Prior art keywords
- curve
- real part
- frequency
- dispersion curve
- correlation
- 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
- 239000006185 dispersion Substances 0.000 title claims description 17
- 238000000034 method Methods 0.000 title claims description 13
- 238000001228 spectrum Methods 0.000 claims description 14
- 239000000284 extract Substances 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- 238000005314 correlation function Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 5
- 238000004613 tight binding model Methods 0.000 claims description 5
- 206010038743 Restlessness Diseases 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 206010044565 Tremor Diseases 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 claims description 2
- 230000002087 whitening effect Effects 0.000 claims description 2
- 238000000605 extraction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
一种提取Rayleigh面波频散曲线的方法。其包括以下步骤：a、对相距为R的地震检波器对记录的噪声互相关信号进行傅里叶变换，得到频谱W(f)；b、对得到的频谱做Gauss窄带滤波得到WG(f)；c、提取频谱WG(f)的实部曲线并找出其零点，记fn为曲线上的第n个零点处的频率；d、找出第一类零阶贝塞尔函数的根记为Zn；e、根据互相关谱实部曲线与第一类零阶贝塞尔函数的关系求单频点处的Rayleigh面波相速度值V(f)=2πfnR/Zn；f、将算得的单频点相速度连接绘制为一条平滑曲线，做为Rayleigh面波的频散曲线。通过贝塞尔函数与互相关谱的关系，有效的提高了频散曲线的计算效率。
Description
技术领域
本发明涉及工程地震勘探或地球深部构造探测领域中的Rayleigh面波探测方法的关键技术，即通过反演地下横波速度结构时Rayleigh面波的频散曲线的提取，我利用地震记录互相关谱实部曲线和第一类零阶贝塞尔函数的关系，提取出了Rayleigh面波的频散曲线。
背景技术
Rayleigh面波是一种沿地球自由表面传播的弹性波,其沿地表传播的速度与波及深度内的介质弹性参数有关,包括介质的密度、纵波和横波速度，而主要的影响参数是介质的横波速度。Rayleigh波在层状介质中具有频散特征, 其传播速度将随着频率的改变而发生变化，根据Rayleigh面波的这一特征,可以求得地层不同深度的弹性参数。
频散曲线的提取Rayleigh面波勘探中的关键步骤，简单来说就是求取对应某频率下瑞雷面波的相速度问题。传统的提取频散曲线方法的思想为，对不同检波器或地震台站实际接收到的Rayleigh面波信号做fourier变换得到其相位，计算不同检波器或台站处的相位差，再根据检波器或台站间的相对距离即可求出某频点的相速度，假设p1台站处的Rayleigh面波信号为R1(t)，p2台站处的Rayleigh面波信号为R2(t)，分别做傅里叶变换的到R1和R2的相位谱R1(f,φ)和R2(f,φ),若两台站间的距离为S，则某个频率处的相速度的可由下式算得：
其中n为周期的倍数，计算时须对n进行估计，即相位解缠。
目前的计算方法通过估计一个速度v的区间范围，从而确定n的可能值，最后带入上式算出相速度，然而在估计值不合理的情况下，n值无法准确得到，即意味着相位的旅行周期无法确定，故计算出的相速度会存在偏差。
发明内容
本发明主要是解决现有计算频散曲线方法中无法准确确定相位旅行周期n的问题，通过互相关谱实部与第一类零阶贝塞尔函数的对应关系计算相速度，这样就可以避开对相位旅行周期n的估计，提高计算结果的准确性。
仍然假设两个台站记录的Rayleigh面波信号为R1(t)和R2(t)，由于信号为离散采样，则他们的互相关函数可表示为：
得到互相关函数后计算其频谱C( QUOTE )，
根据Aki(1957)文献中42式:
该式指出噪声互相关谱和第一类零阶贝塞尔函数的对应关系，实际计算中考虑数据处理的非线性效应，我们通过互相关谱的实部曲线和贝塞尔函数的根的对应关系计算相速度：
通过该式计算出各个频率 QUOTE 处的相速度C( QUOTE )后，即可得到Rayleigh面波的频散曲线。
附图说明：
图1是本专利提取频散曲线流程图；
图2是本专利实施效果图。
具体实施方式：
我们基于背景噪声的方法采集数据，然后进行数据预处理并计算互相关得到台站间的经验格林函数，然后根据互相关谱实部与贝塞尔函数的对应关系求出Rayleigh面波的相速度，绘制频散曲线(图1，图2)。
1数据采集：
在研究区域布设台站对或台阵，进行长时背景噪声观测，根据目标探测深度选取合适的台站间距和采样率。
2数据预处理：
对得到的背景噪声信号进行谱白化、滤波等预处理过程，然后将各个台站的微动信号分段，分段时长根据探测目标深度选取。
3 互相关运算：
将分段后的信号按相同时间段内两两进行互相关运算，然后将相同台站对的不同时间段互相关运算的结果进行叠加，叠加后的互相关函数即可视为台站间的经验格林函数。
4 提取谱实部曲线：
对经验格林函数进行傅里叶变换得到互相关谱，然后提取出谱实部曲线。
5 计算相速度：
实际数据中计算得到的实部曲线与第一类零阶贝塞尔函数振幅有差异，但是实部曲线的零点与贝塞尔函数的根仍然一一对应，因此我们将零点处对应频率的相速度值一一求出。
6 绘制频散曲线：
将得到的各个零点处对应频率的相速度值绘制在V-f平面坐标内，连接各个频点的相速度值，即得到Rayleigh面波的频散曲线。
Claims (2)
1.一种提取频散曲线的方法，其特征在于，所述方法包括：
a数据采集：
在研究区域布设台站对或台阵，进行长时背景噪声观测，根据目标探测深度选取合适的台站间距和采样率；
b数据预处理：
对得到的背景噪声信号进行谱白化、滤波等预处理过程，然后将各个台站的微动信号分段，分段时长根据探测目标深度选取；
c 互相关运算：
将分段后的信号按相同时间段内两两进行互相关运算，然后将相同台站对的不同时间段互相关运算的结果进行叠加，叠加后的互相关函数即可视为台站间的经验格林函数；
d 提取谱实部曲线：
对经验格林函数进行傅里叶变换得到互相关谱，然后提取出谱实部曲线；
e 计算相速度：
实际数据中计算得到的实部曲线与第一类零阶贝塞尔函数振幅有差异，但是实部曲线的零点与贝塞尔函数的根仍然一一对应，因此我们将零点处对应频率的相速度值一一求出；
f 绘制频散曲线：
将得到的各个零点处对应频率的相速度值绘制在V-f平面坐标内，连接各个频点的相速度值，即得到Rayleigh面波的频散曲线。
2.根据权利要求1所述的Rayleigh面波频散曲线提取方法，其特征在于：所述方法通过计算互相关函数谱的实部曲线，根据实部曲线零点与零阶贝塞尔函数根的对应关系计算相速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410580207.0A CN104678435A (zh) | 2014-10-27 | 2014-10-27 | 一种提取Rayleigh面波频散曲线的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410580207.0A CN104678435A (zh) | 2014-10-27 | 2014-10-27 | 一种提取Rayleigh面波频散曲线的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104678435A true CN104678435A (zh) | 2015-06-03 |
Family
ID=53313755
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410580207.0A Pending CN104678435A (zh) | 2014-10-27 | 2014-10-27 | 一种提取Rayleigh面波频散曲线的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104678435A (zh) |
Cited By (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106934183A (zh) * | 2017-04-21 | 2017-07-07 | 王兵 | 频散曲线确定方法和装置，及纵横波速度确定方法和装置 |
| CN107092030A (zh) * | 2017-05-12 | 2017-08-25 | 武汉港湾工程质量检测有限公司 | 通过瑞雷波探测软黏土包的方法 |
| CN108318918A (zh) * | 2017-12-29 | 2018-07-24 | 福建省建筑设计研究院有限公司 | 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用 |
| CN108572390A (zh) * | 2018-05-22 | 2018-09-25 | 武汉市市政建设集团有限公司 | 利用表面波谱扰动预测浅部洞穴的探测方法 |
| CN109188504A (zh) * | 2018-08-08 | 2019-01-11 | 武汉市工程科学技术研究院 | 一种地下管廊探测方法 |
| CN109477904A (zh) * | 2016-06-22 | 2019-03-15 | 休斯敦大学系统 | 地震或声波频散的非线性信号比较和高分辨率度量 |
| CN109471172A (zh) * | 2018-12-26 | 2019-03-15 | 中国科学院地球化学研究所 | 一种基于同相轴形态差异的面波提纯方法及装置 |
| CN109521090A (zh) * | 2018-10-18 | 2019-03-26 | 天津大学 | 一种激光无损表征薄膜杨氏模量的优化方法 |
| WO2019071515A1 (zh) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | 面波勘探方法及终端设备 |
| CN109763831A (zh) * | 2018-12-13 | 2019-05-17 | 武汉市工程科学技术研究院 | 隧道盾构施工检测方法 |
| CN110058299A (zh) * | 2018-09-14 | 2019-07-26 | 南方科技大学 | 地震定位方法、装置及终端设备 |
| CN110426741A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
| WO2020029015A1 (zh) * | 2018-08-06 | 2020-02-13 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
| CN111045076A (zh) * | 2019-12-10 | 2020-04-21 | 核工业北京地质研究院 | 多模式瑞雷波频散曲线并行联合反演方法 |
| CN111290017A (zh) * | 2020-03-04 | 2020-06-16 | 南方科技大学 | 一种震电波场联合提取瑞雷波频散特征的面波勘探方法 |
| CN111766625A (zh) * | 2020-07-06 | 2020-10-13 | 中国科学技术大学 | 一种基于深度学习的地震背景噪声频散曲线的提取方法 |
| CN112051610A (zh) * | 2020-10-21 | 2020-12-08 | 中国地质大学（北京） | 一种矢量场多模式面波频散计算方法及系统 |
| CN112748459A (zh) * | 2019-10-29 | 2021-05-04 | 中国石油天然气集团有限公司 | 频散曲线自动拾取方法及装置 |
| CN112861721A (zh) * | 2021-02-09 | 2021-05-28 | 南方科技大学 | 一种自动提取背景噪声频散曲线的方法及装置 |
| CN112882098A (zh) * | 2021-01-20 | 2021-06-01 | 长江勘测规划设计研究有限责任公司 | 一种天然源面波速度计算及地形校正方法 |
| CN112946753A (zh) * | 2019-12-11 | 2021-06-11 | 中国石油天然气集团有限公司 | 近地表地层结构的分析系统及方法 |
| CN112987090A (zh) * | 2019-12-02 | 2021-06-18 | 中国石油天然气集团有限公司 | 面波频散曲线拾取方法及装置 |
| CN114185093A (zh) * | 2021-12-07 | 2022-03-15 | 中国石油大学(北京) | 一种基于瑞雷面波反演的近地表速度模型建立方法及装置 |
| CN116400406A (zh) * | 2023-04-21 | 2023-07-07 | 中国地震局地球物理研究所 | 一种基于阵列的被动源多模式面波频散曲线提取方法 |
| CN116540298A (zh) * | 2023-07-05 | 2023-08-04 | 中国科学院地质与地球物理研究所 | 一种海洋漏能振型面波频散谱计算方法、系统和电子设备 |
| CN116577829A (zh) * | 2023-05-15 | 2023-08-11 | 中国矿业大学(北京) | 一种基于背景噪音频散曲线自动化提取方法 |
| US11754744B2 (en) | 2020-03-04 | 2023-09-12 | Southern University Of Science And Technology | Surface wave prospecting method for jointly extracting Rayleigh wave frequency dispersion characteristics by seismoelectric field |
| CN118483743A (zh) * | 2024-05-06 | 2024-08-13 | 南方科技大学 | 一种基于应变场提取面波频散曲线的方法 |
-
2014
- 2014-10-27 CN CN201410580207.0A patent/CN104678435A/zh active Pending
Cited By (45)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11209566B2 (en) | 2016-06-22 | 2021-12-28 | University Of Houston System | Nonlinear signal comparison and high-resolution measurement of seismic or acoustic wave dispersion |
| CN109477904A (zh) * | 2016-06-22 | 2019-03-15 | 休斯敦大学系统 | 地震或声波频散的非线性信号比较和高分辨率度量 |
| CN106934183A (zh) * | 2017-04-21 | 2017-07-07 | 王兵 | 频散曲线确定方法和装置，及纵横波速度确定方法和装置 |
| CN107092030B (zh) * | 2017-05-12 | 2019-04-12 | 武汉港湾工程质量检测有限公司 | 通过瑞雷波探测软黏土包的方法 |
| CN107092030A (zh) * | 2017-05-12 | 2017-08-25 | 武汉港湾工程质量检测有限公司 | 通过瑞雷波探测软黏土包的方法 |
| CN109923440A (zh) * | 2017-10-12 | 2019-06-21 | 南方科技大学 | 面波勘探方法及终端设备 |
| WO2019071515A1 (zh) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | 面波勘探方法及终端设备 |
| CN108318918A (zh) * | 2017-12-29 | 2018-07-24 | 福建省建筑设计研究院有限公司 | 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用 |
| CN108572390A (zh) * | 2018-05-22 | 2018-09-25 | 武汉市市政建设集团有限公司 | 利用表面波谱扰动预测浅部洞穴的探测方法 |
| US11137511B2 (en) | 2018-08-06 | 2021-10-05 | Southern University Of Science And Technology | Active source surface wave prospecting method, surface wave exploration device and computer-readable storage medium |
| CN111164462B (zh) * | 2018-08-06 | 2022-05-06 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
| CN111164462A (zh) * | 2018-08-06 | 2020-05-15 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
| WO2020029015A1 (zh) * | 2018-08-06 | 2020-02-13 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
| CN109188504A (zh) * | 2018-08-08 | 2019-01-11 | 武汉市工程科学技术研究院 | 一种地下管廊探测方法 |
| CN109188504B (zh) * | 2018-08-08 | 2020-04-17 | 武汉市工程科学技术研究院 | 一种地下管廊探测方法 |
| CN110058299A (zh) * | 2018-09-14 | 2019-07-26 | 南方科技大学 | 地震定位方法、装置及终端设备 |
| CN109521090A (zh) * | 2018-10-18 | 2019-03-26 | 天津大学 | 一种激光无损表征薄膜杨氏模量的优化方法 |
| CN109763831A (zh) * | 2018-12-13 | 2019-05-17 | 武汉市工程科学技术研究院 | 隧道盾构施工检测方法 |
| CN109471172A (zh) * | 2018-12-26 | 2019-03-15 | 中国科学院地球化学研究所 | 一种基于同相轴形态差异的面波提纯方法及装置 |
| CN110426741A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
| CN110426741B (zh) * | 2019-08-02 | 2021-07-16 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
| CN112748459B (zh) * | 2019-10-29 | 2024-05-28 | 中国石油天然气集团有限公司 | 频散曲线自动拾取方法及装置 |
| CN112748459A (zh) * | 2019-10-29 | 2021-05-04 | 中国石油天然气集团有限公司 | 频散曲线自动拾取方法及装置 |
| CN112987090A (zh) * | 2019-12-02 | 2021-06-18 | 中国石油天然气集团有限公司 | 面波频散曲线拾取方法及装置 |
| CN112987090B (zh) * | 2019-12-02 | 2024-05-28 | 中国石油天然气集团有限公司 | 面波频散曲线拾取方法及装置 |
| CN111045076A (zh) * | 2019-12-10 | 2020-04-21 | 核工业北京地质研究院 | 多模式瑞雷波频散曲线并行联合反演方法 |
| CN112946753A (zh) * | 2019-12-11 | 2021-06-11 | 中国石油天然气集团有限公司 | 近地表地层结构的分析系统及方法 |
| CN111290017A (zh) * | 2020-03-04 | 2020-06-16 | 南方科技大学 | 一种震电波场联合提取瑞雷波频散特征的面波勘探方法 |
| US11754744B2 (en) | 2020-03-04 | 2023-09-12 | Southern University Of Science And Technology | Surface wave prospecting method for jointly extracting Rayleigh wave frequency dispersion characteristics by seismoelectric field |
| CN111766625B (zh) * | 2020-07-06 | 2021-12-14 | 中国科学技术大学 | 一种基于深度学习的地震背景噪声频散曲线的提取方法 |
| CN111766625A (zh) * | 2020-07-06 | 2020-10-13 | 中国科学技术大学 | 一种基于深度学习的地震背景噪声频散曲线的提取方法 |
| CN112051610A (zh) * | 2020-10-21 | 2020-12-08 | 中国地质大学（北京） | 一种矢量场多模式面波频散计算方法及系统 |
| CN112882098A (zh) * | 2021-01-20 | 2021-06-01 | 长江勘测规划设计研究有限责任公司 | 一种天然源面波速度计算及地形校正方法 |
| CN112882098B (zh) * | 2021-01-20 | 2022-03-11 | 长江勘测规划设计研究有限责任公司 | 一种天然源面波速度计算及地形校正方法 |
| CN112861721A (zh) * | 2021-02-09 | 2021-05-28 | 南方科技大学 | 一种自动提取背景噪声频散曲线的方法及装置 |
| CN112861721B (zh) * | 2021-02-09 | 2024-05-07 | 南方科技大学 | 一种自动提取背景噪声频散曲线的方法及装置 |
| CN114185093B (zh) * | 2021-12-07 | 2023-05-12 | 中国石油大学(北京) | 一种基于瑞雷面波反演的近地表速度模型建立方法及装置 |
| CN114185093A (zh) * | 2021-12-07 | 2022-03-15 | 中国石油大学(北京) | 一种基于瑞雷面波反演的近地表速度模型建立方法及装置 |
| CN116400406A (zh) * | 2023-04-21 | 2023-07-07 | 中国地震局地球物理研究所 | 一种基于阵列的被动源多模式面波频散曲线提取方法 |
| CN116400406B (zh) * | 2023-04-21 | 2023-12-19 | 中国地震局地球物理研究所 | 一种基于阵列的被动源多模式面波频散曲线提取方法 |
| CN116577829A (zh) * | 2023-05-15 | 2023-08-11 | 中国矿业大学(北京) | 一种基于背景噪音频散曲线自动化提取方法 |
| CN116577829B (zh) * | 2023-05-15 | 2024-06-04 | 中国矿业大学(北京) | 一种基于背景噪音频散曲线自动化提取方法 |
| CN116540298A (zh) * | 2023-07-05 | 2023-08-04 | 中国科学院地质与地球物理研究所 | 一种海洋漏能振型面波频散谱计算方法、系统和电子设备 |
| CN116540298B (zh) * | 2023-07-05 | 2023-09-12 | 中国科学院地质与地球物理研究所 | 一种海洋漏能振型面波频散谱计算方法、系统和电子设备 |
| CN118483743A (zh) * | 2024-05-06 | 2024-08-13 | 南方科技大学 | 一种基于应变场提取面波频散曲线的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104678435A (zh) | 一种提取Rayleigh面波频散曲线的方法 | |
| JP6945895B2 (ja) | 表面波探査方法および端末デバイス | |
| WO2020029015A1 (zh) | 一种人工源面波勘探方法、面波勘探装置及终端设备 | |
| CN102866421B (zh) | 识别小断距断点的散射波叠前成像方法 | |
| RU2694621C1 (ru) | Способ и устройство для обработки сейсмических данных | |
| CN102937721B (zh) | 利用初至波走时的有限频层析成像方法 | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN104237945B (zh) | 一种地震资料自适应高分辨处理方法 | |
| CN107894613B (zh) | 弹性波矢量成像方法、装置、存储介质及设备 | |
| Liu et al. | Extracting seismic attenuation coefficients from cross-correlations of ambient noise at linear triplets of stations | |
| CN102590862B (zh) | 补偿吸收衰减的叠前时间偏移方法 | |
| CN103984011A (zh) | 一种动态q补偿偏移方法 | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN104199088B (zh) | 一种提取入射角道集的方法及系统 | |
| CN107450103B (zh) | 一种基于边界积分反算子的鬼波压制方法 | |
| CN115993641A (zh) | 一种提取被动源面波频散曲线的方法 | |
| CN106257309B (zh) | 叠后地震数据体处理方法及装置 | |
| CN101825722B (zh) | 一种鲁棒的地震信号瞬时频率的估计方法 | |
| CN104280774B (zh) | 一种单频地震散射噪声的定量分析方法 | |
| CN106324702A (zh) | 一种地震干涉法成像观测系统设计的定量评价方法 | |
| Guan et al. | Improving the retrieval of high-frequency surface waves using convolution-based three-station interferometry for dense linear arrays | |
| CN117665930A (zh) | 一种基于面波频散反传与叠加的震源定位方法及系统 | |
| EP3232234B1 (en) | Method and apparatus for estimating surface wave coda using time-reversal experiments | |
| CN112764100B (zh) | 储层油气范围检测方法及装置 | |
| CN102890288B (zh) | 一种地震波层速度反演方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150603 |
|
| WD01 | Invention patent application deemed withdrawn after publication |