CN105842731B - 基于波场延拓原理的起伏地表组合震源波场定向方法 - Google Patents
基于波场延拓原理的起伏地表组合震源波场定向方法 Download PDFInfo
- Publication number
- CN105842731B CN105842731B CN201610308447.4A CN201610308447A CN105842731B CN 105842731 B CN105842731 B CN 105842731B CN 201610308447 A CN201610308447 A CN 201610308447A CN 105842731 B CN105842731 B CN 105842731B
- Authority
- CN
- China
- Prior art keywords
- relief surface
- wave field
- elevation
- focus
- domain
- 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
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. analysis, for interpretation, for correction
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).
其中:i是虚数单位,i2=-1,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 CN105842731A (zh) | 2016-08-10 |
CN105842731B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538081B (zh) * | 2020-06-05 | 2021-05-25 | 吉林大学 | 一种地震数据初至波的外推方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
CN102478663A (zh) * | 2010-11-23 | 2012-05-30 | 中国科学院地质与地球物理研究所 | 一种三维地震观测系统偏移噪声获取方法及装置 |
CN102890290A (zh) * | 2012-09-25 | 2013-01-23 | 中国石油天然气股份有限公司 | 一种起伏地表条件下的叠前深度偏移方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009077440A2 (en) * | 2007-12-14 | 2009-06-25 | Shell Internationale Research Maatschappij B.V. | Method of processing data obtained from seismic prospecting |
EP2326971A4 (en) * | 2008-08-11 | 2017-06-14 | Exxonmobil Upstream Research Company | Removal of surface-wave noise in seismic data |
-
2016
- 2016-05-11 CN CN201610308447.4A patent/CN105842731B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
CN102478663A (zh) * | 2010-11-23 | 2012-05-30 | 中国科学院地质与地球物理研究所 | 一种三维地震观测系统偏移噪声获取方法及装置 |
CN102890290A (zh) * | 2012-09-25 | 2013-01-23 | 中国石油天然气股份有限公司 | 一种起伏地表条件下的叠前深度偏移方法 |
Non-Patent Citations (1)
Title |
---|
起伏地表组合震源地震波场定向方法;巩向博 等;《地球物理学报》;20141231;第57卷(第12期);第4150-4156页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105842731A (zh) | 2016-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bauer et al. | Utilizing diffractions in wavefront tomography | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
Gray et al. | True-amplitude Gaussian-beam migration | |
Zeng et al. | An improved vacuum formulation for 2D finite-difference modeling of Rayleigh waves including surface topography and internal discontinuities | |
Veeken | Seismic stratigraphy and depositional facies models | |
US8553499B2 (en) | Migration-based illumination determination for AVA risk assessment | |
Pan et al. | Love-wave waveform inversion in time domain for shallow shear-wave velocity | |
Rawlinson et al. | Inversion of seismic refraction and wide-angle reflection traveltimes for three-dimensional layered crustal structure | |
US8659974B2 (en) | System and method of 3D salt flank VSP imaging with transmitted waves | |
Willis et al. | A novel application of time-reversed acoustics: Salt-dome flank imaging using walkaway VSP surveys | |
Improta et al. | Seismic imaging of complex structures by non-linear traveltime inversion of dense wide-angle data: application to a thrust belt | |
Lin et al. | Accurate diffraction imaging for detecting small-scale geologic discontinuities | |
Wang et al. | Multiparameter TTI tomography of P-wave reflection and VSP data | |
Van Veldhuizen et al. | Acquisition geometry analysis in complex 3D media | |
Mecking et al. | Cavity detection by SH-wave full-waveform inversion—A reflection-focused approach | |
CN107656308B (zh) | 一种基于时间深度扫描的共散射点叠前时间偏移成像方法 | |
Toxopeus et al. | Simulating migrated and inverted seismic data by filtering a geologic model | |
Talukdar et al. | Sub-basalt imaging of hydrocarbon-bearing Mesozoic sediments using ray-trace inversion of first-arrival seismic data and elastic finite-difference full-wave modeling along Sinor–Valod profile of Deccan Syneclise, India | |
Pan et al. | Multichannel analysis of Love waves in a 3D seismic acquisition system | |
Zhu et al. | Recent applications of turning-ray tomography | |
Guo et al. | Topography-dependent eikonal tomography based on the fast-sweeping scheme and the adjoint-state technique | |
Li et al. | Waveform inversion of seismic first arrivals acquired on irregular surface | |
Liu et al. | Tomographic velocity model building of the near surface with velocity-inversion interfaces: A test using the Yilmaz model | |
Tang et al. | Single-well acoustic imaging in anisotropic formations | |
CN108919351A (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 | ||
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 |