CN108303730A - 一种工程三维双源面波勘探方法及系统 - Google Patents
一种工程三维双源面波勘探方法及系统 Download PDFInfo
- Publication number
- CN108303730A CN108303730A CN201810008773.2A CN201810008773A CN108303730A CN 108303730 A CN108303730 A CN 108303730A CN 201810008773 A CN201810008773 A CN 201810008773A CN 108303730 A CN108303730 A CN 108303730A
- Authority
- CN
- China
- Prior art keywords
- surface wave
- excitation point
- double source
- seismic
- wave
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000005284 excitation Effects 0.000 claims abstract description 80
- 238000011835 investigation Methods 0.000 claims abstract description 14
- 239000006185 dispersion Substances 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种工程三维双源面波勘探方法及系统,在勘探区域布设震源激发点和检波器,在震源激发点处激发产生地震面波,通过所述检波器采集地震面波,获得地震面波数据;根据所述震源激发点和所述检波器的位置建立所述震源激发点和所述检波器的坐标,获得震源激发点坐标和检波器坐标;根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集;根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。本发明能够通过建立所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集,然后再根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
Description
技术领域
本发明涉及面波勘探领域,特别是涉及一种工程三维双源面波勘探方法及系统。
背景技术
面波勘探,也称弹性波频率测深,是一种浅层地震勘探方法,所谓面波,是一种特殊的地震波,属于一种地滚波,可分为瑞利波(R波)和勒夫波(L波)。面波能量占震源激发总能量的2/3,其中,R波在振动波组中的能量最强,振幅最大,频率最低,容易识别也更易于测量。目前,瑞利面波已广泛用于无损检测和场地评价领域。
由于受观测系统和数据采集技术的局限,二维面波勘探方法采集的地震数据无法直观地从中了解勘探目标体的平面分布特征及空间三维形态。
发明内容
本发明的目的是提供一种能够直观地了解勘探目标体的平面分布特征及空间三维形态的工程三维双源面波勘探方法及系统。
为实现上述目的,本发明提供了如下方案:
一种工程三维双源面波勘探方法,所述勘探方法包括:
在勘探区域布设震源激发点和检波器,在震源激发点处激发产生地震面波,通过所述检波器采集地震面波,获得地震面波数据;
根据所述震源激发点和所述检波器的位置建立所述震源激发点和所述检波器的坐标,获得震源激发点坐标和检波器坐标;
根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集;
根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
可选的,所述在勘探区域布设震源激发点和检波器具体包括:
在所述勘探区域内布设多个所述检波器和多个所述震源激发点,分布多束所述检波器,分布多排震源激发点,每排所述震源激发点的连线与每束所述检波器的连线垂直。
可选的,所述检波器采集地震面波,获得地震面波数据具体包括:
设置地震数据采集窗口,所述地震数据采集窗口的大小为M×N的矩形区域,所述震源激发点设置有炮,设定一排所述震源激发点的连线按照横向排布,一束所述检波器的连线按照竖向排布;
依次激发第i排炮,采集横向第i到第M+i-1排炮,竖向第i到第N+i-1束检波器的区域范围内的地震面波,获取第i地震面波数据,i的取值为1,2,3,...,n;
根据所述第i地震面波数据获取所述地震面波数据。
可选的,所述根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集具体包括:
修订所述检波器和所述震源激发点的坐标;
在所述时间域剔除所述地震面波数据的干扰波,获得有效地震面波数据;
在频率-波数域、频率-速度域分析所述有效地震面波数据,获得所述有效地震面波数据的频散曲线和深度-速度曲线;
根据多条所述频散曲线集成,获得三维双源面波数据集。
可选的,所述根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态具体包括:
获取所述勘探区域内切片测线位置的切片坐标;
根据所述切片坐标对所述三维双源面波数据集切片,获取所述切片测线位置处的地下探测目标的平面分布范围和三维形态。
为了实现上述目的,本发明还提供了如下方案:
一种工程三维双源面波勘探系统,所述勘探系统包括:
地震面波数据获取模块用于获取通过所述检波器采集的地震面波,获得地震面波数据;
坐标建立模块与所述获取模块连接,所述坐标建立模块用于根据所述震源激发点和所述检波器的位置建立所述震源激发点和所述检波器的坐标,获得震源激发点坐标和检波器坐标;
三维双源面波数据集获取模块与所述坐标建立模块连接,所述三维双源面波数据集获取模块用于根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集;
三维形态构建模块与所述三维双源面波数据集获取模块连接,所述三维形态构建模块用于根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供了工程三维双源面波勘探方法及系统,在震源激发点处激发产生地震面波,通过所述检波器采集地震面波,根据所述震源激发点和所述检波器的位置建立所述震源激发点和所述检波器的坐标,再根据地震面波数据获得三维双源面波数据集,根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
工程三维双源面波勘探方法节省了施工成本,明显地简化了现场工作程序,缩短地质勘查周期,提高了野外工作质量和效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的工程三维双源面波勘探方法的流程图;
图2为本发明勘探区域布设震源激发点和检波器的布局图;
图3为本发明根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集的流程图;
图4为本发明提供的工程三维双源面波勘探系统的结构图;
图5为本发明多条所述频散曲线图;
图6为本发明三维双源面波数据图;
图7为本发明地下探测目标的平面分布范围和三维形态图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种能够直观地了解勘探目标体的平面分布特征及空间三维形态的工程三维双源面波勘探方法及系统。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示的一种工程三维双源面波勘探方法的流程图,所述勘探方法包括:
步骤100:在勘探区域布设震源激发点8和检波器7,在震源激发点8处激发产生地震面波,通过所述检波器7采集地震面波,获得地震面波数据,所述检波器检测的为低频信号。
地震面波的激发有两种:
第一种属于人工源激发,在震源激发点8处激发产生地震面波。
第二种直接通过天然源激发地震面波,所述天然源包括海水的潮涌、人类的活动等。
步骤200:根据所述震源激发点8和所述检波器7的位置建立所述震源激发点8和所述检波器7的坐标,获得震源激发点坐标和检波器坐标。
步骤300:根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集。
步骤400:根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
如图2所示的勘探区域布设震源激发点8和检波器7的布局图,所述在勘探区域布设震源激发点8和检波器7具体包括:
在所述勘探区域内布设多个所述检波器7和多个所述震源激发点8,分布多束所述检波器7,分布多排震源激发点8,每排所述震源激发点8的连线与每束所述检波器7的连线垂直。
排列设计采用束线型观测系统,相邻两个所述检波器7之间的距离不大于所述勘探区域的最小地层的厚度。所述检波器7的排列长度大于所述勘探区域勘探目标的深度。
所述检波器7采集地震面波,获得地震面波数据具体包括:
设置地震数据采集窗口,所述地震数据采集窗口的大小为M×N的矩形区域6,所述震源激发点设置有炮,设定一排所述震源激发点8的连线按照横向排布,一束所述检波器7的连线按照竖向排布。
依次激发第i排炮,采集横向第i到第M+i-1排炮,竖向第i到第N+i-1束检波器的区域范围内的地震面波,获取第i地震面波数据,i的取值为1,2,3,...,n;根据所述第i地震面波数据获取所述地震面波数据,直至完全覆盖整个勘探区域的地震数据的采集。
如图3所示的所述步骤300:根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集具体包括:
步骤301:修订所述检波器7和所述震源激发点8的坐标,检查与定义每个检波器7对应的坐标,每个所述震源激发点8对应的坐标。
步骤302:在所述时间域剔除所述地震面波数据的干扰波,对高阶所述地震面波数据和低阶所述地震面波数据有效分离,获得有效地震面波数据,所述检波器7的排列长度越长,所述地震面波数据分离的效果越好。
步骤303:利用FK法在频率-波数域、频率-速度域分析所述有效地震面波数据,获得所述有效地震面波数据的频散曲线和深度-速度曲线。
步骤304:根据多条所述频散曲线集成,获得三维双源面波数据集。
如图5所示的将多条所述频散曲线集成,获得了如图6所示的三维双源面波数据集。
如图7所示的地下探测目标的平面分布范围和三维形态图,所述根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态具体包括:
获取所述勘探区域内切片测线位置的切片坐标;
根据所述切片坐标对所述三维双源面波数据集切片,获取所述切片测线位置处的地下探测目标的平面分布范围和三维形态。
为了实现上述目的,本发明还提供了如下方案:
一种工程三维双源面波勘探系统,在勘探区域布设震源激发点8和检波器7,在震源激发点8处激发产生地震面波,不使用天然源,通过所述检波器7采集地震面波;所述勘探系统包括:
地震面波数据获取模块用于获取通过所述检波器采集的地震面波,获得地震面波数据;
坐标建立模块2与所述获取模块1连接,所述坐标建立模块2用于根据所述震源激发点8和所述检波器7的位置建立所述震源激发点8和所述检波器7的坐标,获得震源激发点坐标和检波器坐标;
三维双源面波数据集获取模块3与所述坐标建立模块2连接,所述三维双源面波数据集获取模块3用于根据所述震源激发点坐标、所述检波器坐标和所述地震波数据获得三维双源面波数据集;
三维形态构建模块4与所述三维双源面波数据集获取模块3连接,所述三维形态构建模块4用于根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (6)
1.一种工程三维双源面波勘探方法,其特征在于,所述勘探方法包括:
在勘探区域布设震源激发点和检波器,在震源激发点处激发产生地震面波,通过所述检波器采集地震面波,获得地震面波数据;
根据所述震源激发点和所述检波器的位置建立所述震源激发点和所述检波器的坐标,获得震源激发点坐标和检波器坐标;
根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集;
根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
2.根据权利要求1所述的一种工程三维双源面波勘探方法,其特征在于,所述在勘探区域布设震源激发点和检波器具体包括:
在所述勘探区域内布设多个所述检波器和多个所述震源激发点,分布多束所述检波器,分布多排震源激发点,每排所述震源激发点的连线与每束所述检波器的连线垂直。
3.根据权利要求1所述的一种工程三维双源面波勘探方法,其特征在于,所述检波器采集地震面波,获得地震面波数据具体包括:
设置地震数据采集窗口,所述地震数据采集窗口的大小为M×N的矩形区域,所述震源激发点设置有炮,设定一排所述震源激发点的连线按照横向排布,一束所述检波器的连线按照竖向排布;
依次激发第i排炮,采集横向第i到第M+i-1排炮,竖向第i到第N+i-1束检波器的区域范围内的地震面波,获取第i地震面波数据,i的取值为1,2,3,...,n;
根据所述第i地震面波数据获取所述地震面波数据。
4.根据权利要求1所述的一种工程三维双源面波勘探方法,其特征在于,所述根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维双源面波数据集具体包括:
修订所述检波器和所述震源激发点的坐标;
在所述时间域剔除所述地震面波数据的干扰波,获得有效地震面波数据;
在频率-波数域、频率-速度域分析所述有效地震面波数据,获得所述有效地震面波数据的频散曲线和深度-速度曲线;
根据多条所述频散曲线集成,获得三维双源面波数据集。
5.根据权利要求1所述的一种工程三维双源面波勘探方法,其特征在于,所述根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态具体包括:
获取所述勘探区域内切片测线位置的切片坐标;
根据所述切片坐标对所述三维双源面波数据切片,获取所述切片测线位置处的地下探测目标的平面分布范围和三维形态。
6.一种工程三维双源面波勘探系统,其特征在于,所述勘探系统包括:
地震面波数据获取模块用于获取通过所述检波器采集的地震面波,获得地震面波数据;
坐标建立模块与所述获取模块连接,所述坐标建立模块用于根据所述震源激发点和所述检波器的位置建立所述震源激发点和所述检波器的坐标,获得震源激发点坐标和检波器坐标;
三维双源面波数据集获取模块与所述坐标建立模块连接,所述三维双源面波数据集获取模块用于根据所述震源激发点坐标、所述检波器坐标和所述地震面波数据获得三维面波数据集;
三维形态构建模块与所述三维双源面波数据集获取模块连接,所述三维形态构建模块用于根据所述三维双源面波数据集获得所述勘探区域不同位置处的地下探测目标的平面分布范围和三维形态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810008773.2A CN108303730A (zh) | 2018-01-04 | 2018-01-04 | 一种工程三维双源面波勘探方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810008773.2A CN108303730A (zh) | 2018-01-04 | 2018-01-04 | 一种工程三维双源面波勘探方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108303730A true CN108303730A (zh) | 2018-07-20 |
Family
ID=62868660
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810008773.2A Pending CN108303730A (zh) | 2018-01-04 | 2018-01-04 | 一种工程三维双源面波勘探方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108303730A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112731519A (zh) * | 2019-10-28 | 2021-04-30 | 中国石油天然气集团有限公司 | 确定颤动时间区间的方法和装置 |
CN113640881A (zh) * | 2021-08-18 | 2021-11-12 | 中南大学 | 多偏移距二维横向高分辨率瞬态面波探测方法 |
CN114517471A (zh) * | 2022-04-20 | 2022-05-20 | 西南交通大学 | 高铁路基有效硬层检测方法、装置、设备及可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6418380B1 (en) * | 1997-10-10 | 2002-07-09 | Compagnie Generale De Geophysique | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration |
CN102288990A (zh) * | 2011-08-03 | 2011-12-21 | 浙江大学 | 超浅层三维地震反射波勘探方法 |
CN106324682A (zh) * | 2016-09-19 | 2017-01-11 | 中国石油化工股份有限公司 | 一种应用于冻土层地区的表层结构调查方法 |
CN107356668A (zh) * | 2016-05-10 | 2017-11-17 | 北京市水利规划设计研究院 | 一种坝体安全鉴定的检测方法 |
-
2018
- 2018-01-04 CN CN201810008773.2A patent/CN108303730A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6418380B1 (en) * | 1997-10-10 | 2002-07-09 | Compagnie Generale De Geophysique | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration |
CN102288990A (zh) * | 2011-08-03 | 2011-12-21 | 浙江大学 | 超浅层三维地震反射波勘探方法 |
CN107356668A (zh) * | 2016-05-10 | 2017-11-17 | 北京市水利规划设计研究院 | 一种坝体安全鉴定的检测方法 |
CN106324682A (zh) * | 2016-09-19 | 2017-01-11 | 中国石油化工股份有限公司 | 一种应用于冻土层地区的表层结构调查方法 |
Non-Patent Citations (1)
Title |
---|
孙秀荣等: "《双源面波数据采集与处理系统及其应用效果》", 《上海国土资源》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112731519A (zh) * | 2019-10-28 | 2021-04-30 | 中国石油天然气集团有限公司 | 确定颤动时间区间的方法和装置 |
CN113640881A (zh) * | 2021-08-18 | 2021-11-12 | 中南大学 | 多偏移距二维横向高分辨率瞬态面波探测方法 |
CN114517471A (zh) * | 2022-04-20 | 2022-05-20 | 西南交通大学 | 高铁路基有效硬层检测方法、装置、设备及可读存储介质 |
CN114517471B (zh) * | 2022-04-20 | 2022-08-09 | 西南交通大学 | 高铁路基有效硬层检测方法、装置、设备及可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102369752B1 (ko) | 지하 구조물의 검출 | |
Parker et al. | Active‐source seismic tomography at the Brady geothermal field, Nevada, with dense nodal and fiber‐optic seismic arrays | |
Bormann et al. | Seismic sources and source parameters | |
US8219320B2 (en) | Seismic attributes for reservoir localization | |
US7590491B2 (en) | Signal integration measure for seismic data | |
CN104280775B (zh) | 一种基于全波形矢量偏移叠加的微地震监测定位方法 | |
Cheraghi et al. | Passive seismic reflection interferometry: A case study from the Aquistore CO2 storage site, Saskatchewan, Canada | |
Schwenk et al. | Surface-wave methods for anomaly detection | |
CN103336297B (zh) | 微破裂向量扫描方法 | |
CN106094029A (zh) | 利用偏移距矢量片地震数据预测储层裂缝的方法 | |
CN107728204A (zh) | 基于叠前纵波各向异性的裂缝预测方法及系统 | |
CN108303730A (zh) | 一种工程三维双源面波勘探方法及系统 | |
CN105467445B (zh) | 一种三维地震观测系统的建立方法及装置 | |
CN107356668A (zh) | 一种坝体安全鉴定的检测方法 | |
CN108919337A (zh) | 城市地下空间地震数据采集系统 | |
CN105445787B (zh) | 一种最优方位子体相干的裂缝预测方法 | |
CN108919351A (zh) | 基于逆时聚焦原理进行观测系统双向聚焦性的评价方法 | |
CN108303741A (zh) | 基于水检积分的地震双检合并鸣震压制方法 | |
Strobbia et al. | Double-array refraction microtremors | |
Taipodia et al. | A review of active and passive MASW techniques | |
Sicking et al. | Forecasting reservoir performance by mapping seismic emissions | |
Mohamed et al. | Near-surface site characterization at Quriyat City, Sultanate of Oman using HVSR and MASW techniques | |
CN104375180B (zh) | 一种地震数据处理方法、装置及系统 | |
Lehujeur et al. | Vertical seismic profiling using double-beamforming processing of nonuniform anthropogenic seismic noise: The case study of Rittershoffen, Upper Rhine Graben, France | |
CN106154354B (zh) | 一种成像道集的应用方法 |
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 |
Application publication date: 20180720 |
|
RJ01 | Rejection of invention patent application after publication |