CN103792572A - 利用地震记录中的面波信息调查表层结构的方法 - Google Patents

利用地震记录中的面波信息调查表层结构的方法 Download PDF

Info

Publication number
CN103792572A
CN103792572A CN201210418186.3A CN201210418186A CN103792572A CN 103792572 A CN103792572 A CN 103792572A CN 201210418186 A CN201210418186 A CN 201210418186A CN 103792572 A CN103792572 A CN 103792572A
Authority
CN
China
Prior art keywords
frequency
wave
data
surface structure
ground roll
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
Application number
CN201210418186.3A
Other languages
English (en)
Inventor
冯刚
段卫星
宁鹏鹏
卢湘鹏
刘斌
张光德
何京国
张在武
田根海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petrochemical Corp Shengli Offshore Drilling Co Geophysics Exploration And Develop
China Petroleum and Chemical Corp
Original Assignee
China Petrochemical Corp Shengli Offshore Drilling Co Geophysics Exploration And Develop
China Petroleum and Chemical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petrochemical Corp Shengli Offshore Drilling Co Geophysics Exploration And Develop, China Petroleum and Chemical Corp filed Critical China Petrochemical Corp Shengli Offshore Drilling Co Geophysics Exploration And Develop
Priority to CN201210418186.3A priority Critical patent/CN103792572A/zh
Publication of CN103792572A publication Critical patent/CN103792572A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种利用地震记录中的面波信息调查表层结构的方法,该利用地震记录中的面波信息调查表层结构的方法包括:步骤1,输入原始数据,进行原始数据处理;步骤2,将步骤1中处理后的数据在频率波数域进行面波波形的提取,获得基阶组分的面波频散数据;步骤3,将步骤2中的得到的所述面波频散数据在距离频率域数据叠加成频散文件;以及步骤4,根据步骤3中在距离频率域所得到的频散文件进行反演拟合,得到地层模型综合成果图。该利用地震记录中的面波信息调查表层结构的方法克服传统表层结构方法对地层速度和波阻抗的严格要求,获得了准确的浅表层厚度和速度资料。

Description

利用地震记录中的面波信息调查表层结构的方法
技术领域
 本发明涉及地震勘探领域,特别是涉及到一种利用地震记录中的面波信息调查表层结构的方法。
背景技术
地震勘探近地表结构调查中常用的方法各有优缺点,岩性取心比较可靠,但成本高,且对地层具有破坏性;折射波法和反射波法对于波阻抗差异较小的地质体界面反映较弱,不容易分辨,特别是折射波法要求下层介质速度一定要大于上层介质速度,如果地层存在低速夹层、速度倒转,折射法将无能为力。面波勘探法是一种新型的浅层地震勘探方法,不受上述物性条件制约,仅要求各层介质间存在横波速度差异,因此在地层划分方面有较好的分辨能力。面波法除了能弥补传统方法的不足,还具有快速,经济,易采集处理、对环境无破坏,抗干扰强以及进行连续表层探测等诸多优点。
面波分为拉夫波和瑞雷波,与拉夫波相比,瑞雷波在振动波组中能量最强、振幅最大、频率最低,容易识别也易于测量,因此面波勘探一般指的是瑞雷波勘探。以下面波勘探指的是瑞雷波勘探。
迄今为止,虽然国内外对于利用瑞雷波调查近地层结构开展了大量的研究,发表了许多论文,但关于石油勘探中面波应用的公开文献数量较少,从事这方面研究的力量也不多。从这些公开发表的文献上看,目前国内石油勘探领域对利用瑞雷面波调查表层结构的研究主要可分为两条路线。一条路线是直接应用工程勘探中瞬态瑞雷波勘探技术调查地表结构,如米晓利等人的研究(米晓利,王欣等.瞬态瑞雷波技术在地震勘探方面的应用研究.物探装备,2002.12(2):95~100)。这种方式目前已经取得了一定的实际效果,但需要开展专门的野外工作,增加了石油地震勘探的投资。另一条路线是借鉴瞬态瑞雷波勘探的思路,利用单炮记录中的面波信息来调查地表结构,如裴江云、席拥军、周熙襄、王振国等人的研究(周熙襄,王振国等. 利用瑞利面波调查表层结构的可行性探讨.石油地球物理勘探,2004.39(02):181~186,122~248)。这种方式充分挖掘单炮记录数据,不需要专门投资,把常规生产中的干扰波转化成有用信息。如果研究成功,将会具有广阔的应用前景和较高的经济效益。但目前该路线主要处于理论探索阶段,还有一些理论与技术上的问题等待解决。
因此,亟待研制一套利用地震记录中的面波提高表层结构调查精度的流程与方法。
发明内容
本发明的目的是提供一种应用于地震勘探的新的浅表层结构的利用地震记录中的面波信息调查表层结构的方法,该方法的建立基于地震单炮面波信息调查表层结构的处理流程,克服传统表层结构方法对地层速度和波阻抗的严格要求,获得了准确的浅表层厚度和速度资料。   
本发明的目的可通过如下技术措施来实现:利用地震记录中的面波信息调查表层结构的方法,该利用地震记录中的面波信息调查表层结构的方法包括:步骤1,输入原始数据,进行原始数据处理;步骤2,将步骤1中处理后的数据在频率波数域进行面波波形的提取,获得基阶组分的面波频散数据;步骤3,将步骤2中的得到的所述面波频散数据在距离频率域数据叠加成频散文件;以及步骤4,根据步骤3中在距离频率域所得到的频散文件进行反演拟合,得到地层模型综合成果图。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,在进行原始数据处理时,整理地震采集获得的单端激发的共炮点等道距单炮记录,包括核定采集距离参数、识别和清理干扰波型、观察频谱特征、设定频率波数转换的频段上限。
在步骤2中,进行面波波形的提取时,在显出的频率波数谱图形上,人工选择、追踪相应于提取波型的幅度峰,自动找出该频率的幅度峰脊,读取相速度,圈出幅度峰的范围,获得该面波频散数据。
在步骤3中,通过该频散文件观察排列道全部叠加的该面波频散数据或偏移距由近到远的每两道间的所述面波频散数据,并将该频散文件存盘,
在步骤3中,将该频散文件转入深度速度域反演页面作数值反演。
在步骤4中,将在距离频率域所得到频散数据进行反演拟合,拟合满意后将反演模型参数置入该频散文件,得出近地表层数、厚度及剪切波速等模型参数。
面波作为地震资料中传统意义的干扰波,实际上在地层介质中传播时,携带所经过介质丰富的地质信息,将这些反映地层属性的有用信息提取出来,可解决地质问题,尤其是浅表层地质问题。面波勘探利用面波在非均匀介质中的频散特性,得出面波与表层结构(厚度与速度)之间的关系。本发明中的利用地震记录中的面波信息调查表层结构的方法,将面波法引入到地震勘探低(降)速带结构的调查中,建立一种新的地震勘探表层调查新技术,对于进行近地表结构研究具有很高的利用价值。与常规表层结构调查方法相比,本发明所提供的面波法不但能避开基于折射波理论方法的困境和微测井方法带来的高昂成本,还能利用生产单炮记录中的面波建立高密度的表层结构剖面,从而解决小折射和微测井测点间隔较稀的问题。
 
附图说明
图1为本发明的利用地震记录中的面波信息调查表层结构的方法的流程图;
图2为本发明的一具体实施例中单边激发单炮解编记录的示意图;
图3为在时域距清理单炮记录上干扰波的显示图;
图4为在频率波数域提取基阶面波波形时的显示图;
图5为在距离频率域相关叠加获得的频散数据显示图;
图6为反演地层模型求得的表层综合成果图。
 
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合所附图式,作详细说明如下。
本发明的利用地震记录中的面波信息调查表层结构的方法的发明原理为,设面波的频率为fi,相邻检波器记录的面波的时间差为△t或相位差为△Φ,则相邻道△X长度内面波的传播速度为:
在满足空间采样定理的条件下,测量范围N △x内平均波速为:
Figure 2012104181863100002DEST_PATH_IMAGE002
Figure 2012104181863100002DEST_PATH_IMAGE004
在同一测点测量出一系列频率fi的VRi值,就可以得到一条VR--f曲线,即所谓的频散曲线。
面波法勘探实质上是根据面波传播的频散特性,利用人工震源激发产生多种频率成分的面波,寻找出波速随频率的变化关系,最终确定出地表岩土的面波速度随场点坐标(x,z)的变化关系,以解决浅表层结构勘探等问题。
如图1所示,图1为本发明的利用地震记录中的面波信息调查表层结构的方法的流程图。
步骤1,输入原始数据,进行原始数据处理。整理地震采集获得的单端激发的共炮点等道距单炮记录,包括核定采集距离参数、识别和清理干扰波型、观察频谱特征、设定频率波数转换的频段上限等。图2为本发明的一具体实施例中单边激发单炮解编记录的示意图,即步骤1中输入的原始数据示意图,根据该图可进行数据整理。如图3,图3为时域距清理单炮记录上干扰波的显示图。从图3可以看出干扰波清除后的单炮记录形态。流程进入到步骤2。
步骤2,将步骤1中整理后的数据在频率波数域进行面波波形的提取,获得基阶组分的面波频散数据。处理过程是在显出的频率波数谱图形上,人工选择、追踪相应于提取波型的幅度峰,自动找出该频率的幅度峰脊,读取相速度,圈出幅度峰的范围,将此频率波数范围提供给下步距离频率域叠加页面。如图4所示,图4为在频率波数域提取基阶面波波形时的显示图。从图4可以看出白色阴影圈出的频率波数范围连续稳定,包含了面波的基阶组分,其中的白点显示出幅度峰脊,对应着不同频率的基阶面波频散相速。流程进入到步骤3。
步骤3,将步骤2中的选定的频散数据在距离频率域数据叠加成一个频散文件,提供的处理功能容许全面观察排列道全部叠加的频散数据或偏移距由近到远的每两道间的频散数据,可以将数据存盘,也可以直接转入深度速度域反演页面作数值反演。如图5所示,图5为在距离频率域相关叠加获得的频散数据显示图。从图5可以看出全部远近道互相关叠加求得的频散数据,可反映浅层近地表波速参数,但是并不能反映深层波速参数。。流程进入到步骤4。
步骤4、根据步骤3中在距离频率域所得到频散文件进行反演拟合,得到地层模型综合成果图。将在距离频率域所得到频散数据进行反演拟合,拟合满意后将反演模型参数置入频散数据文件,得出近地表层数、厚度及剪切波速等模型参数。也就是说,距离频率域得到的频散数据或频散文件是由步骤2获得,在步骤3进行转化。步骤2得到频散数据是提取面波基阶组分后获得,在距离频率域进行转化获得频散文件以进行步骤4的反演。如图6所示,图6为本发明的一实施例中反演地层模型求得的表层综合成果图。从图6可以看出连续1号线点为面波频散数据点,2号线点为模型正演频散数据点,3号线为地层波速断面,左下方为地层模型参数及模型拟合度列表。列表包含:频散数据点、反演地层模型参数、地层波速断面、模型正演频散数据点及模型拟合度。流程结束。

Claims (6)

1.利用地震记录中的面波信息调查表层结构的方法,其特征在于,该利用地震记录中的面波信息调查表层结构的方法包括:
步骤1,输入原始数据,进行原始数据处理;
步骤2,将步骤1中处理后的数据在频率波数域进行面波波形的提取,获得基阶组分的面波频散数据;
步骤3,将步骤2中的得到的所述面波频散数据在距离频率域数据叠加成频散文件;以及
步骤4,根据步骤3中在距离频率域所得到的频散文件进行反演拟合,得到地层模型综合成果图。
2.根据权利要求1所述的利用地震记录中的面波信息调查表层结构的方法,其特征在于,在步骤1中,在进行原始数据处理时,整理地震采集获得的单端激发的共炮点等道距单炮记录,包括核定采集距离参数、识别和清理干扰波型、观察频谱特征、设定频率波数转换的频段上限。
3.根据权利要求1所述的利用地震记录中的面波信息调查表层结构的方法,其特征在于,在步骤2中,进行面波波形的提取时,在显出的频率波数谱图形上,人工选择、追踪相应于提取波型的幅度峰,自动找出该频率的幅度峰脊,读取相速度,圈出幅度峰的范围,获得该面波频散数据。
4.根据权利要求1所述的利用地震记录中的面波信息调查表层结构的方法,其特征在于,在步骤3中,通过该频散文件观察排列道全部叠加的该面波频散数据或偏移距由近到远的每两道间的所述面波频散数据,并将该频散文件存盘。
5.根据权利要求1所述的利用地震记录中的面波信息调查表层结构的方法,其特征在于,在步骤3中,将该频散文件转入深度速度域反演页面作数值反演。
6.根据权利要求1所述的利用地震记录中的面波信息调查表层结构的方法,其特征在于,在步骤4中,将在距离频率域所得到频散数据进行反演拟合,拟合满意后将反演模型参数置入该频散文件,得出近地表层数、厚度及剪切波速等模型参数。
CN201210418186.3A 2012-10-29 2012-10-29 利用地震记录中的面波信息调查表层结构的方法 Pending CN103792572A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210418186.3A CN103792572A (zh) 2012-10-29 2012-10-29 利用地震记录中的面波信息调查表层结构的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210418186.3A CN103792572A (zh) 2012-10-29 2012-10-29 利用地震记录中的面波信息调查表层结构的方法

Publications (1)

Publication Number Publication Date
CN103792572A true CN103792572A (zh) 2014-05-14

Family

ID=50668423

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210418186.3A Pending CN103792572A (zh) 2012-10-29 2012-10-29 利用地震记录中的面波信息调查表层结构的方法

Country Status (1)

Country Link
CN (1) CN103792572A (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105572731A (zh) * 2015-12-21 2016-05-11 中国石油天然气集团公司 一种确定勘探区表层结构信息的方法和装置
CN106094020A (zh) * 2016-05-31 2016-11-09 中国石油天然气集团公司 一种地震反演方法及装置
CN107015271A (zh) * 2017-02-28 2017-08-04 中国石油天然气集团公司 表层结构调查方法和装置
CN107561589A (zh) * 2017-10-25 2018-01-09 中国石油化工股份有限公司 一种近地表横波层速度模型建立方法
CN108181650A (zh) * 2017-12-19 2018-06-19 中国石油天然气集团公司 一种数据处理方法和装置
CN108572390A (zh) * 2018-05-22 2018-09-25 武汉市市政建设集团有限公司 利用表面波谱扰动预测浅部洞穴的探测方法
WO2019071515A1 (zh) * 2017-10-12 2019-04-18 南方科技大学 面波勘探方法及终端设备
CN112285775A (zh) * 2020-10-15 2021-01-29 中国科学院地质与地球物理研究所 一种海底浅部沉积层剪切波速度探测装置及方法
CN112748459A (zh) * 2019-10-29 2021-05-04 中国石油天然气集团有限公司 频散曲线自动拾取方法及装置
CN113640881A (zh) * 2021-08-18 2021-11-12 中南大学 多偏移距二维横向高分辨率瞬态面波探测方法
CN113805234A (zh) * 2021-10-13 2021-12-17 四川省冶金地质勘查院 在被动源地震数据中增强面波的处理方法
CN114063159A (zh) * 2020-08-10 2022-02-18 中国石油化工股份有限公司 一种地震面波速度确定方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2510016A1 (en) * 2001-12-21 2003-07-10 Universite De Sherbrooke Method and algorithm for using surface waves
CN101907727A (zh) * 2010-08-17 2010-12-08 中国科学院地质与地球物理研究所 一种面波多分量转换波静校正方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2510016A1 (en) * 2001-12-21 2003-07-10 Universite De Sherbrooke Method and algorithm for using surface waves
CN101907727A (zh) * 2010-08-17 2010-12-08 中国科学院地质与地球物理研究所 一种面波多分量转换波静校正方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
周熙襄等: "利用瑞利面波调查表层结构的可行性探讨", 《石油地球物理勘探》 *
徐贵来: "地震瑞利面波测深应注意的问题", 《铀矿地质》 *
徐贵来: "瑞雷面波测深主要影响因素分析及应用实例", 《全国成矿理论与深部找矿新方法及勘查开发关键技术交流研讨会》 *
颜俊华等: "利用面波信息调查表层结构的方法", 《新疆石油地质》 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105572731A (zh) * 2015-12-21 2016-05-11 中国石油天然气集团公司 一种确定勘探区表层结构信息的方法和装置
CN106094020A (zh) * 2016-05-31 2016-11-09 中国石油天然气集团公司 一种地震反演方法及装置
CN107015271A (zh) * 2017-02-28 2017-08-04 中国石油天然气集团公司 表层结构调查方法和装置
CN107015271B (zh) * 2017-02-28 2019-02-15 中国石油天然气集团公司 表层结构调查方法和装置
CN109923440A (zh) * 2017-10-12 2019-06-21 南方科技大学 面波勘探方法及终端设备
WO2019071515A1 (zh) * 2017-10-12 2019-04-18 南方科技大学 面波勘探方法及终端设备
CN107561589B (zh) * 2017-10-25 2019-04-30 中国石油化工股份有限公司 一种近地表横波层速度模型建立方法
CN107561589A (zh) * 2017-10-25 2018-01-09 中国石油化工股份有限公司 一种近地表横波层速度模型建立方法
CN108181650A (zh) * 2017-12-19 2018-06-19 中国石油天然气集团公司 一种数据处理方法和装置
CN108572390A (zh) * 2018-05-22 2018-09-25 武汉市市政建设集团有限公司 利用表面波谱扰动预测浅部洞穴的探测方法
CN112748459A (zh) * 2019-10-29 2021-05-04 中国石油天然气集团有限公司 频散曲线自动拾取方法及装置
CN112748459B (zh) * 2019-10-29 2024-05-28 中国石油天然气集团有限公司 频散曲线自动拾取方法及装置
CN114063159A (zh) * 2020-08-10 2022-02-18 中国石油化工股份有限公司 一种地震面波速度确定方法及装置
CN112285775A (zh) * 2020-10-15 2021-01-29 中国科学院地质与地球物理研究所 一种海底浅部沉积层剪切波速度探测装置及方法
CN113640881A (zh) * 2021-08-18 2021-11-12 中南大学 多偏移距二维横向高分辨率瞬态面波探测方法
CN113640881B (zh) * 2021-08-18 2022-06-28 中南大学 多偏移距二维横向高分辨率瞬态面波探测方法
CN113805234A (zh) * 2021-10-13 2021-12-17 四川省冶金地质勘查院 在被动源地震数据中增强面波的处理方法
CN113805234B (zh) * 2021-10-13 2023-12-12 四川省冶金地质勘查院 在被动源地震数据中增强面波的处理方法

Similar Documents

Publication Publication Date Title
CN103792572A (zh) 利用地震记录中的面波信息调查表层结构的方法
Cheng et al. A new passive seismic method based on seismic interferometry and multichannel analysis of surface waves
CN101329405B (zh) 一种简单的多参数地震反演方法
Mi et al. Near-surface imaging from traffic-induced surface waves with dense linear arrays: An application in the urban area of Hangzhou, China
CN101551463B (zh) 三维观测系统噪声压制估算方法
Pan et al. Delineating shallow S‐wave velocity structure using multiple ambient‐noise surface‐wave methods: An example from Western Junggar, China
CN102176054B (zh) 近地表综合信息处理解释方法
CN102176053A (zh) 提升波动方程叠前深度偏移成像效果的方法
CN102109616A (zh) 一种沉积旋回约束的波阻抗反演方法
Zhou et al. Migration velocity analysis and prestack migration of common-transmitter GPR data
CN103758511A (zh) 一种井下逆时偏移成像识别隐蔽储层的方法及装置
Hayashi et al. CMP spatial autocorrelation analysis of multichannel passive surface-wave data
CN112505749A (zh) 一种基于线形台阵多次覆盖的微动数据采集方法
de Ridder et al. Ambient seismic noise tomography at Ekofisk
Zeng et al. High‐resolution shallow structure at brady hot springs using ambient noise tomography (ANT) on a trenched distributed acoustic sensing (DAS) array
Martorana et al. Assessment of a geological model by surface wave analyses
Nardone et al. Shallow velocity structure and site effects at Mt. Vesuvius, Italy, from HVSR and array measurements of ambient vibrations
Dong’an et al. Full wave seismic exploration technology
Su et al. High-resolution seismic processing technique with broadband, wide-azimuth, and high-density seismic data—A case study of thin-sand reservoirs in eastern China
RU2386984C1 (ru) Способ поиска углеводородов
Pandey et al. Developing a distributed acoustic sensing seismic land streamer: Concept and validation
CN112764100B (zh) 储层油气范围检测方法及装置
Jianhua* et al. Applied research of 3D DAS-VSP in Eastern China
Al Hasani et al. Experiences with distributed acoustic sensing using both straight and helically wound fibers in surface-deployed cables—A case history in Groningen, The Netherlands
Bezuidenhout et al. Rayleigh group velocity extraction from ambient seismic noise to map the south Eastern Cape Karoo region, South Africa

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20140514