CN105842731A - 基于波场延拓原理的起伏地表组合震源波场定向方法 - Google Patents
基于波场延拓原理的起伏地表组合震源波场定向方法 Download PDFInfo
- Publication number
- CN105842731A CN105842731A CN201610308447.4A CN201610308447A CN105842731A CN 105842731 A CN105842731 A CN 105842731A CN 201610308447 A CN201610308447 A CN 201610308447A CN 105842731 A CN105842731 A CN 105842731A
- Authority
- CN
- China
- Prior art keywords
- wave field
- relief surface
- seismic source
- elevation
- time
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 230000005284 excitation Effects 0.000 claims abstract description 8
- 238000013213 extrapolation Methods 0.000 claims description 5
- 230000000644 propagated effect Effects 0.000 claims description 4
- 239000002352 surface water Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000006185 dispersion 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/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
技术领域
本发明专利属于勘探地震领域,涉及一种地震勘探组合震源的定向方法,尤其是在复杂起伏地表条件下的地震波场定向。
背景技术:
地震勘探方法是探测地下地质体构造的一种地球物理勘探方法,其重要的环节之一是地震数据的采集。而数据采集质量的好坏与震源激发有密切关系,野外勘探时常用组合震源的方式以提高震源激发的能量和地下目标体反射信号的信噪比。
目前组合震源可考虑以下几种情况:水平地表条件下简单线性组合震源可用组合原理的方向因子公式来确定震源的激发延时;倾角固定的倾斜地表条件下有类似的方向因子公式;复杂起伏地表条件下可用旋转坐标方法的几何路径计算激发延时。如巩向博等,地球物理学报2014,“起伏地表组合震源地震波场定向方法”。然而,对于非均匀复杂起伏地表情况,通过几何关系计算的组合震源激发延时有一定的误差。
发明内容:
本发明针对现有技术的不足,提供了一种基于波场延拓原理的起伏地表条件下组合震源激发定向波场的方法。
为了实现本发明的上述目的,基于波场延拓原理的起伏地表组合震源定向波场方法包括以下步骤:
a、建立起伏地表地球物理参数模型,将模型参数离散化至每个网格节点;
b、记录此模型的地表水平位置与相应高程坐标,沿起伏地表设置N个水平间隔一致的点震源S1,S2,…,Sn,分别记录其高程坐标为zS1,zS2,…zSn;
c、将震源的高程进行排序,取其中的最大高程,作为水平基准面高程zH=max(zS1,zS2,…zSn),max是取最大值函数;
d、从模型最左网格至最右网格,将起伏地表模型向上填充等效介质至水平基准面,等效介质物理参数为起伏地表表层位置处的地球物理参数;
e、将起伏地表的N个地震点震源投影至在水平基准面上,记录起伏地表至水平基准面高程ΔzS1=zH-zS1,ΔzS2=zH-zS2,…,ΔzSn=zH-zSn;
f、根据水平地表组合震源的方向因子公式计算相邻震源的延时时差参数τ;
g、按照延时参数依次激发N个点震源,得到N个时间域波场,分别为u1(x1,zH,t),u2(x2,zH,τ+t),…,u2(x2,zH,(n-1)*τ+t);
h、将此N个时间域波场分别对时间和水平坐标两个维度进行傅里叶变换,得到N个频率波数域波场,分别为U1(kx,zH,ω),U2(kx,zH,ω),…,Un(kx,zH,ω),其中kx是水平波数;
i、根据频率波数域的下行波动方程,将N个频率波数域波场以自激自收方式正向外推(向下延拓)至实际起伏地表高程,结果分别为U1(kx,zS1,ω),U2(kx,zS2,ω),…,Un(kx,zSn,ω)
j、将延拓后N个频率波数域波场叠加,得到频率域起伏地表组合震源定向波场记录;
k、将上步结果进行时间和水平坐标的二维反傅里叶变换,得到时间域定向传播的波场。
有益效果:经模拟试验,本发明公开的起伏地表组合震源波场定向的方法,计算准确,有效得到定向激发的地震波场,并改善了由于地表非均匀性引起的波场计算误差,在实际应用中更具有普遍性。从图中明显看出,本发明通过a~k步骤能够得到任意方向角的地震波场,且起伏地表对定向传播的地震波场的影响很小。
附图说明
图1复杂起伏地表地质模型图
图2组合震源定向波场在0.3秒处的波场快照图。
2a是方向角θ为30度,2b是方向角θ为0度,2c是方向角θ为-30度
具体实施方式
下面结合附图和实例对本发明进一步的详细说明。
一种起伏地表条件下组合震源激发定向波场的方法,在图1所示的非均匀起伏地表模型情况下,包括如下步骤:
a、建立起伏地表地球物理参数模型,将模型参数离散化至每个网格节点;
b、记录此模型的地表水平位置与相应高程坐标,沿起伏地表设置N个水平间隔一致的点震源S1,S2,…,Sn,分别记录其高程坐标为zS1,zS2,…zSn,设震源间距为d;
c、将震源的高程进行排序,取其中的最大高程,作为水平基准面高程zH=max(zS1,zS2,…zSn),max是取最大值函数;
d、从模型最左网格至最右网格,将起伏地表模型向上填充等效介质至水平基准面,等效介质物理参数为起伏地表表层位置处的地球物理参数;
e、将起伏地表的N个地震点震源投影至在水平基准面上,记录起伏地表至水平基准面高程ΔzS1=zH-zS1,ΔzS2=zH-zS2,…,ΔzSn=zH-zSn;
f、设计地震波传播的方向角为θ,根据水平地表组合震源的方向因子公式计算相邻震源的相位延迟β,根据相位延迟与延时关系得到震源的延时时差参数τ=β/ω,其中ω是角频率;
g、按照延时延时依次激发N个点震源,得到N个时间域波场,分别为u1(x1,zH,t),u2(x2,zH,τ+t),…,u2(x2,zH,(n-1)*τ+t);
h、将此N个时间域波场分别对时间和水平坐标两个维度进行傅里叶变换,得到N个频率波数域波场,分别为U1(kx,zH,ω),U2(kx,zH,ω),…,Un(kx,zH,ω),其中kx是水平波数;
i、根据频率波数域的下行波动方程,将N个频率波数域波场以自激自收方式正向外推(向下延拓)至实际起伏地表高程,结果分别为U1(kx,zS1,ω),U2(kx,zS2,ω),…,Un(kx,zSn,ω),其中波场延拓公式为
U(kx,zs1,ω)=-U(kx,zH,ω)exp(ikzΔzs1),
U(kx,zs2,ω)=-U(kx,zH,ω)exp(ikzΔzs2),
…,
U(kx,zsn,ω)=-U(kx,zH,ω)exp(ikzΔzsn).
其中kz是垂直波数,符合频散关系v代表地下介质速度,是空间变量的函数v(x,z);
j、将延拓后N个频率波数域波场叠加,得到频率域起伏地表组合震源定向波场记录;
k、将上步结果进行时间和水平坐标的二维反傅里叶变换,得到时间域定向传播的波场。
Claims (1)
1.一种基于波场延拓原理的起伏地表组合震源波场定向方法,其特征在于,包括以下步骤:
a、建立起伏地表地球物理参数模型,将模型参数离散化至每个网格节点;
b、记录此模型的地表水平位置与相应高程坐标,沿起伏地表设置N个水平间隔一致的点震源S1,S2,…,Sn,分别记录其高程坐标为zS1,zS2,…zSn;
c、将震源的高程进行排序,取其中的最大高程,作为水平基准面高程zH=max(zS1,zS2,…zSn),max是取最大值函数;
d、从模型最左网格至最右网格,将起伏地表模型向上填充等效介质至水平基准面,等效介质物理参数为起伏地表表层位置处的地球物理参数;
e、将起伏地表的N个地震点震源投影至在水平基准面上,记录起伏地表至水平基准面高程ΔzS1=zH-zS1,ΔzS2=zH-zS2,…,ΔzSn=zH-zSn;
f、根据水平地表组合震源的方向因子公式计算相邻震源的延时时差参数τ;
g、按照延时参数依次激发N个点震源,得到N个时间域波场,分别为u1(x1,zH,t),u2(x2,zH,τ+t),…,u2(x2,zH,(n-1)*τ+t);
h、将此N个时间域波场分别对时间和水平坐标两个维度进行傅里叶变换,得到N个频率波数域波场,分别为U1(kx,zH,ω),U2(kx,zH,ω),…,Un(kx,zH,ω),其中kx是水平波数;
i、根据频率波数域的下行波动方程,将N个频率波数域波场以自激自收方式正向外推至实际起伏地表高程,结果分别为U1(kx,zS1,ω),U2(kx,zS2,ω),…,Un(kx,zSn,ω)
j、将延拓后N个频率波数域波场叠加,得到频率域起伏地表组合震源定向波场记录;
k、将j步骤结果进行时间和水平坐标的二维反傅里叶变换,得到时间域定向传播的波场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610308447.4A CN105842731B (zh) | 2016-05-11 | 2016-05-11 | 基于波场延拓原理的起伏地表组合震源波场定向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610308447.4A CN105842731B (zh) | 2016-05-11 | 2016-05-11 | 基于波场延拓原理的起伏地表组合震源波场定向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105842731A true CN105842731A (zh) | 2016-08-10 |
CN105842731B CN105842731B (zh) | 2018-02-23 |
Family
ID=56592179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610308447.4A Expired - Fee Related CN105842731B (zh) | 2016-05-11 | 2016-05-11 | 基于波场延拓原理的起伏地表组合震源波场定向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105842731B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538081A (zh) * | 2020-06-05 | 2020-08-14 | 吉林大学 | 一种地震数据初至波的外推方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
US20100299071A1 (en) * | 2007-12-14 | 2010-11-25 | Kiyashchenko Denis | Method of processing data obtained from seismic prospecting |
US20110069581A1 (en) * | 2008-08-11 | 2011-03-24 | Christine E Krohn | Removal of Surface-Wave Noise In Seismic Data |
CN102478663A (zh) * | 2010-11-23 | 2012-05-30 | 中国科学院地质与地球物理研究所 | 一种三维地震观测系统偏移噪声获取方法及装置 |
CN102890290A (zh) * | 2012-09-25 | 2013-01-23 | 中国石油天然气股份有限公司 | 一种起伏地表条件下的叠前深度偏移方法 |
-
2016
- 2016-05-11 CN CN201610308447.4A patent/CN105842731B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
US20100299071A1 (en) * | 2007-12-14 | 2010-11-25 | Kiyashchenko Denis | Method of processing data obtained from seismic prospecting |
US20110069581A1 (en) * | 2008-08-11 | 2011-03-24 | Christine E Krohn | Removal of Surface-Wave Noise In Seismic Data |
CN102478663A (zh) * | 2010-11-23 | 2012-05-30 | 中国科学院地质与地球物理研究所 | 一种三维地震观测系统偏移噪声获取方法及装置 |
CN102890290A (zh) * | 2012-09-25 | 2013-01-23 | 中国石油天然气股份有限公司 | 一种起伏地表条件下的叠前深度偏移方法 |
Non-Patent Citations (1)
Title |
---|
巩向博 等: "起伏地表组合震源地震波场定向方法", 《地球物理学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538081A (zh) * | 2020-06-05 | 2020-08-14 | 吉林大学 | 一种地震数据初至波的外推方法 |
CN111538081B (zh) * | 2020-06-05 | 2021-05-25 | 吉林大学 | 一种地震数据初至波的外推方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105842731B (zh) | 2018-02-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
Bauer et al. | Utilizing diffractions in wavefront tomography | |
AU2012220584B2 (en) | Sensitivity kernel-based migration velocity analysis in 3D anisotropic media | |
CN105182408B (zh) | 一种合成地震记录的制作方法和装置 | |
CN101036069B (zh) | 振幅保存叠前偏移方法 | |
Martin et al. | High-resolution teleseismic body-wave tomography beneath SE Romania—I. Implications for three-dimensional versus one-dimensional crustal correction strategies with a new crustal velocity model | |
CN102841375A (zh) | 一种复杂条件下基于角度域共成像点道集的层析速度反演方法 | |
EP2548052B1 (en) | System and method of 3d salt flank vsp imaging with transmitted waves | |
AU2004220765B2 (en) | Method for stable estimation of anisotropic parameters for P-wave prestack imaging | |
Improta et al. | Seismic imaging of complex structures by non-linear traveltime inversion of dense wide-angle data: application to a thrust belt | |
Jarillo Michel et al. | Waveform inversion for microseismic velocity analysis and event location in VTI media | |
CN107656308B (zh) | 一种基于时间深度扫描的共散射点叠前时间偏移成像方法 | |
Guo et al. | Topography-dependent eikonal tomography based on the fast-sweeping scheme and the adjoint-state technique | |
Pan et al. | Multichannel analysis of Love waves in a 3D seismic acquisition system | |
Zhu et al. | Recent applications of turning-ray tomography | |
Smithyman et al. | Waveform tomography of field vibroseis data using an approximate 2D geometry leads to improved velocity models | |
Tang et al. | Single-well acoustic imaging in anisotropic formations | |
Cao et al. | Full-wave directional illumination analysis in the frequency domain | |
CN105842731B (zh) | 基于波场延拓原理的起伏地表组合震源波场定向方法 | |
De Matteis et al. | 1-D P-velocity models of Mt. Vesuvius volcano from the inversion of TomoVes96 first arrival time data | |
Bruno | High-resolution seismic imaging in complex environments: A comparison among common-reflection-surface stack, common-midpoint stack, and prestack depth migration at the Ilva-Bagnoli brownfield site, Campi Flegrei, Italy | |
Bosch et al. | Full-waveform inversion of intensity-focused seismic data | |
Barnes et al. | Diving wave tomography: A robust method for velocity estimation in a foothills geological context | |
He et al. | The effect of strong near-surface scattering on the quality of seismic imaging | |
Feighner | Three-dimensional seismic imaging of the Ryepatch geothermal reservoir |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180223 Termination date: 20200511 |