CN104101897A - 一种利用矢量合成纵波进行地震勘探的方法 - Google Patents

一种利用矢量合成纵波进行地震勘探的方法 Download PDF

Info

Publication number
CN104101897A
CN104101897A CN201410329213.9A CN201410329213A CN104101897A CN 104101897 A CN104101897 A CN 104101897A CN 201410329213 A CN201410329213 A CN 201410329213A CN 104101897 A CN104101897 A CN 104101897A
Authority
CN
China
Prior art keywords
wave
seismic
vector
compressional wave
reflection
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
CN201410329213.9A
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.)
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Original Assignee
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201410329213.9A priority Critical patent/CN104101897A/zh
Publication of CN104101897A publication Critical patent/CN104101897A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种利用矢量合成纵波进行地震勘探的方法,向待勘探区发射地震纵波后,利用埋设于地震纵波接收点的检波器在垂直地面方向和平行地面方向上分别对反射地震纵波的垂直矢量纵波和水平矢量纵波进行接收,根据所述反射地震纵波的出射方向与铅垂线或水平面的夹角将所述反射地震纵波的垂直矢量纵波和水平矢量纵波进行矢量合成得到所述反射地震纵波,通过分析合成的所述反射地震纵波的信息进行勘探。本发明利用矢量分解和合成的原理,通过对非垂直出射的反射地震纵波进行多角度的接收、合成,可以获得真正的反射地震波出射时的信息,提高有效波的能量和信噪比,进而获得有效的反射地震波信息并进行进一步地勘探。

Description

一种利用矢量合成纵波进行地震勘探的方法
技术领域
本发明属于地球勘探物理学领域,更具体地讲,涉及一种利用矢量合成纵波进行地震勘探的方法。
背景技术
地震勘探是利用人工方法产生并发射地震波,地震波向地下传播过程中遇到不同岩层的分界面时,部分地震波会反射回地面,在地面利用检波器收集这些反射回来的地震波,并对所收集的反射波信息进行研究以查明地下地质目标的一种勘探方法。
目前,在石油、天然气的地震勘探中,仍然以纵波地震勘探方法为主。纵波地震勘探是利用收集的纵波信息进行勘探的方法,但现在使用的纵波地震勘探方法有一个假设前提,即在地表存在低降速地层。具体来说,地震纵波在地下的传播过程中遵循斯奈尔(Snell)定律,即当地震纵波从地下高速地层传播到地表低降速地层时,根据斯奈尔定律,地震纵波将垂直出射,因此利用一个垂直地表埋置的检波器就可以很好的接收到垂直出射的地震纵波。
但是,在碳酸盐岩出露区等其它区域,由于地表并不存在低降速地层,或者低降速地层的地震波传播速度与地下高速地层的地震波传播速度差别不大,因而导致地震纵波在地表出射时并不是垂直地出射,而是与地面呈一定角度地出射。此时如果仍然采用一个垂直地表埋置的检波器,就不能很好地接收到与地面呈一定角度出射的反射地震纵波,进而无法有效地进行反射波信息的分析和进一步的勘探。
发明内容
针对现有技术中的不足,本发明旨在解决上述技术问题中的一个或多个。
本发明的目的在于提供一种能够解决地震勘探中地震纵波非垂直出射所造成反射的地震纵波接收不全的问题的一种地震勘探方法,从而获得真正的地震纵波出射时的信息并进行有效分析和勘探。
为了实现上述目的,本发明提供了一种利用矢量合成纵波进行地震勘探的方法,向待勘探区发射地震纵波后,利用埋设于地震纵波接收点的检波器在垂直地面方向和平行地面方向上分别对反射地震纵波的垂直矢量纵波和水平矢量纵波进行接收,根据所述反射地震纵波的出射方向与铅垂线或水平面的夹角将所述反射地震纵波的垂直矢量纵波和水平矢量纵波进行矢量合成得到所述反射地震纵波,通过分析合成的所述反射地震纵波的信息进行勘探。
根据本发明的利用矢量合成纵波进行地震勘探的方法的一个实施例,所述方法具体包括以下步骤:
a、对待勘探区进行近地表结构调查,建立近地表速度结构模型,并获得地震纵波在近地表地层中的传播速度v0、地震纵波在地下勘探深度范围内地层中传播的均方根速度v和反射地震纵波的双程旅行时间t0
b、通过震源向所述待勘探区发射地震纵波,选取多个地震纵波接收点并分别埋设多分量检波器,每个多分量检波器将接收在垂直地面方向和平行地面方向上分别接收反射地震纵波的垂直矢量纵波和水平矢量纵波;
c、利用下式1计算出射到第i个多分量检波器的反射地震纵波的出射角θi,其中,1≤i≤n且i为正整数,n为多分量检波器的数量,
式中, θ i = arctg x i v 0 t 0 v 2 ——式1,
θi为出射到第i个多分量检波器的反射地震纵波的出射方向与铅垂线之间的夹角,单位是°且0°<θi<90°;
xi为从震源到第i个多分量检波器的距离,单位是m;
v0为地震纵波在近地表地层中的传播速度,单位是m/s;
v为地震纵波在地下勘探深度范围内地层中传播的均方根速度,单位是m/s;
t0为反射地震纵波的双程旅行时间。单位是s;
d、利用步骤c计算得到的θi和下式2矢量合成出射到第i个多分量检波器的反射地震纵波,
Zi=Xi*cosθi+Yi*sinθi——式2,
式中,Zi为出射到第i个多分量检波器的反射地震纵波;
Xi为第i个多分量检波器在垂直地面方向上接收到的反射地震纵波的垂直矢量纵波;
Yi为第i个多分量检波器在平行地面方向上接收到的反射地震纵波的水平矢量纵波;
e、分析步骤d矢量合成的出射到各多分量检波器的反射地震纵波信息进行勘探。
根据本发明的利用矢量合成纵波进行地震勘探的方法的一个实施例,所述方法适用于反射地震纵波非垂直出射的待勘探区。
本发明利用矢量分解和合成的原理,通过对非垂直出射的反射地震纵波进行多角度的接收、合成,可以获得真正的反射地震纵波出射时的信息,提高有效波的能量和信噪比,进而获得有效的反射地震纵波信息并进行进一步地勘探。本发明适用于所有涉及反射的地震纵波非垂直出射的勘探区,具有广泛的应用前景。
具体实施方式
在下文中,将详细描述本发明的利用矢量合成纵波进行地震勘探的方法。
本发明具体是根据矢量分解、合成的原理,对于非垂直出射的反射地震纵波,先通过接收其垂直地面方向和平行地面方向上的矢量地震纵波,再通过计算非垂直出射的反射地震纵波的出射角并继而合成得到该非垂直出射的反射地震纵波。具体来说,可以在地震波接收点利用多分量检波器分别对垂直地面方向和平行地面方向的矢量地震纵波进行接收,并利用近地表的相关模型计算反射地震纵波的出射角,再根据该出射角值对垂直和水平方向接收到的矢量地震纵波进行矢量合成。
根据本发明的示例性实施例,所述利用矢量合成纵波进行地震勘探的方法具体为:
向待勘探区发射地震纵波后,利用埋设于地震纵波接收点的检波器在垂直地面方向和平行地面方向上分别对反射地震纵波的垂直矢量纵波和水平矢量纵波进行接收,根据所述反射地震纵波的出射方向与铅垂线或水平面的夹角将所述反射地震纵波的垂直矢量纵波和水平矢量纵波进行矢量合成得到所述反射地震纵波,通过分析合成的所述反射地震纵波的信息进行勘探。
需要说明的是,其中涉及的地震纵波是指利用人工方法产生的地震纵波,反射地震纵波是指地震波发射到地下之后反射回地表的地震纵波,垂直矢量纵波是指反射地震纵波在垂直地面方向上的矢量纵波,水平矢量纵波是指反射地震纵波在平行地面方向上的矢量纵波。
根据本发明的一个实施例,上述利用矢量合成纵波进行地震勘探的方法具体可以包括以下步骤:
步骤a:
首先,对待勘探区进行近地表结构调查,建立近地表速度结构模型,并计算地震纵波在近地表地层的传播速度v0,同时利用已有的地震资料获得地震纵波在地下勘探深度范围内地层中传播的均方根速度v和反射地震纵波的双程旅行时间t0
地震纵波在近地表地层的传播速度v0也称为地表层速度或近地表速度,地震纵波在地下勘探深度范围内地层中传播的均方根速度v也称为地下底层均方根速度,反射纵波的双程旅行时间t0,即地震纵波从地面传到地下目标地层后再从目标地层传回地面的时间。获取上述两个速度值和一个时间值是为了后续计算反射波的出射角,其中,近地表结构调查是直接通过测量地层厚度和传播时间,由速度=厚度÷传播时间的公式计算速度,这是目前能获得最准确的近地表速度和厚度信息的方法,而均方根速度则是对地下主要地层的传播速度先平方、再平均、然后再开方获得;时间t0可以由检波器接收的地震纵波资料上直接读取。上述近地表速度结构模型的建立、各速度值的计算和时间值的获得均可以采用本领域的常规方法获得,本发明不限于此。
由于在反射波的传播过程中,当反射波遇到速度不一样的介质界面时,其传播的角度会发生变化,根据透射波遵循斯奈尔定律:因此可以利用上、下地层的速度计算出反射波出射的角度。但一般情况下,近地表地层的速度较低,而其下面地层的速度较高,两者之间差异较大,因此对反射波的出射角的影响较大,建立高精度的近地表速度结构模型更利于计算反射波的出射角。
步骤b:
然后,通过震源向待勘探区发射地震纵波,选取多个地震波接收点并分别埋设多分量检波器。其中,每个多分量检波器分别在垂直地面方向和平行地面方向上接受反射地震纵波的垂直矢量纵波和水平矢量纵波。具体地,上述多分量检波器为地震纵波接收器,并且应至少包括在垂直地面方向上的接收器和在平行地面方向上的接收器。
步骤c:
之后,利用下式1计算出射到第i个多分量检波器的反射纵波的出射角θi,1≤i≤n且i为正整数,n为多分量检波器的数量,——式1。
式中,θi为出射到第i个多分量检波器的反射地震纵波的出射方向与铅垂线之间的夹角,单位是°且0°<θi<90°;xi为从震源到第i个多分量检波器的距离,单位是m;v0为地震纵波在地表地层中的传播速度,单位是m/s;v为地震纵波在勘探深度范围内地层中传播的均方根速度,单位是m/s;t0为反射地震纵波的双程旅行时间,单位是s。
在步骤c中,利用上述参数可以计算得到出射到每个多分量检波器的反射地震纵波的出射角。具体地,v0、v、t0可以在步骤a中得到,xi可以通过具体测量得到。
步骤d:
再利用步骤c计算得到的θi和下式2矢量合成出射到第i个多分量检波器的反射地震纵波,
Zi=Xi*cosθi+Yi*sinθi——式2
式中,Zi为出射到第i个多分量检波器的反射地震纵波;Xi为第i个多分量检波器在垂直地面方向上接收到的反射地震纵波的垂直矢量纵波;Yi为第i个多分量检波器在平行地面方向上接收到的反射地震纵波的水平矢量纵波。
在步骤d中,利用出射到每个多分量检波器的反射地震纵波的出射角以及相应的多分量检波器接收到的反射地震纵波的垂直矢量纵波和水平矢量纵波就可以合成计算出射到每个多分量检波器的反射地震纵波。
步骤e:
最后,分析步骤d矢量合成得到的出射到各多分量检波器的反射地震纵波信息进行勘探。其中,具体的分析步骤和勘探步骤可以为本领域常规的方式,在此不进行过多的描述,并且本发明也不对具体的分析、勘探方法进行限定。
本发明的上述方法适用于所有的反射地震纵波非垂直出射的待勘探区,尤其可广泛应用于地表的地震纵波传播速度较高的地震勘探区。
综上所述,本发明解决了反射地震纵波非垂直出射时,垂直地表的地震纵波接收器(检波器)不能较好地接收出射的地震纵波的难题,本发明的方法有利于提高资料信噪比,进而提高地震勘探的精度。
尽管上面已经通过结合示例性实施例描述了本发明,但是本领域技术人员应该清楚,在不脱离权利要求所限定的精神和范围的情况下,可对本发明的示例性实施例进行各种修改和改变。

Claims (3)

1.一种利用矢量合成纵波进行地震勘探的方法,其特征在于,向待勘探区发射地震纵波后,利用埋设于地震纵波接收点的检波器在垂直地面方向和平行地面方向上分别对反射地震纵波的垂直矢量纵波和水平矢量纵波进行接收,根据所述反射地震纵波的出射方向与铅垂线或水平面的夹角将所述反射地震纵波的垂直矢量纵波和水平矢量纵波进行矢量合成得到所述反射地震纵波,通过分析合成的所述反射地震纵波的信息进行勘探。
2.根据权利要求1所述的利用矢量合成纵波进行地震勘探的方法,其特征在于,所述方法具体包括以下步骤:
a、对待勘探区进行近地表结构调查,建立近地表速度结构模型,并获得地震纵波在近地表地层中的传播速度v0、地震纵波在地下勘探深度范围内地层中传播的均方根速度v和反射地震纵波的双程旅行时间t0
b、通过震源向所述待勘探区发射地震纵波,选取多个地震纵波接收点并分别埋设多分量检波器,每个多分量检波器分别在垂直地面方向和平行地面方向上接受反射地震纵波的垂直矢量纵波和水平矢量纵波;
c、利用下式1计算出射到第i个多分量检波器的反射地震纵波的出射角θi,其中,1≤i≤n且i为正整数,n为多分量检波器的数量,
θ i = arctg x i v 0 t 0 v 2 ——式1,
式中,θi为出射到第i个多分量检波器的反射地震纵波的出射方向与铅垂线之间的夹角,单位是°且0°<θi<90°;
xi为从震源到第i个多分量检波器的距离,单位是m;
v0为地震纵波在近地表地层中的传播速度,单位是m/s;
v为地震纵波在地下勘探深度范围内地层中传播的均方根速度,单位是m/s;
t0为反射地震纵波的双程旅行时间,单位是s;
d、利用步骤c计算得到的θi和下式2矢量合成出射到第i个多分量检波器的反射地震纵波,
Zi=Xi*cosθi+Yi*sinθi——式2,
式中,Zi为出射到第i个多分量检波器的反射地震纵波;
Xi为第i个多分量检波器在垂直地面方向上接收到的反射地震纵波的垂直矢量纵波;
Yi为第i个多分量检波器在平行地面方向上接收到的反射地震纵波的水平矢量纵波;
e、分析步骤d矢量合成得到的出射到各多分量检波器的反射地震纵波信息进行勘探。
3.根据权利要求1或2所述的利用矢量合成纵波进行地震勘探的方法,其特征在于,所述方法适用于反射地震纵波非垂直出射的待勘探区。
CN201410329213.9A 2014-07-11 2014-07-11 一种利用矢量合成纵波进行地震勘探的方法 Pending CN104101897A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410329213.9A CN104101897A (zh) 2014-07-11 2014-07-11 一种利用矢量合成纵波进行地震勘探的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410329213.9A CN104101897A (zh) 2014-07-11 2014-07-11 一种利用矢量合成纵波进行地震勘探的方法

Publications (1)

Publication Number Publication Date
CN104101897A true CN104101897A (zh) 2014-10-15

Family

ID=51670191

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410329213.9A Pending CN104101897A (zh) 2014-07-11 2014-07-11 一种利用矢量合成纵波进行地震勘探的方法

Country Status (1)

Country Link
CN (1) CN104101897A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109358361A (zh) * 2018-10-18 2019-02-19 河海大学 一种提取地震射线矢量的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4766574A (en) * 1987-03-31 1988-08-23 Amoco Corporation Method for depth imaging multicomponent seismic data
CN103513277A (zh) * 2013-09-27 2014-01-15 中国石油天然气股份有限公司 一种地震地层裂隙裂缝密度反演方法及系统
CN103592683A (zh) * 2013-10-17 2014-02-19 中国石油天然气集团公司 一种用于转换波地震资料的地表一致性振幅补偿方法及装置
CN103645497A (zh) * 2013-11-08 2014-03-19 中国石油大学(北京) 一种基于出射角的多分量波场分离方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4766574A (en) * 1987-03-31 1988-08-23 Amoco Corporation Method for depth imaging multicomponent seismic data
CN103513277A (zh) * 2013-09-27 2014-01-15 中国石油天然气股份有限公司 一种地震地层裂隙裂缝密度反演方法及系统
CN103592683A (zh) * 2013-10-17 2014-02-19 中国石油天然气集团公司 一种用于转换波地震资料的地表一致性振幅补偿方法及装置
CN103645497A (zh) * 2013-11-08 2014-03-19 中国石油大学(北京) 一种基于出射角的多分量波场分离方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡立新 等: ""井中近源三分量波场调查试验与数据分析"", 《石油物探》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109358361A (zh) * 2018-10-18 2019-02-19 河海大学 一种提取地震射线矢量的方法

Similar Documents

Publication Publication Date Title
CN107526101B (zh) 一种获取地震反射波的采集和处理方法
CN102033242B (zh) 一种深层倾斜裂缝储层地震振幅预测方法
CN103645503B (zh) 一种三维时间域照明分析及振幅补偿方法
CN102759746B (zh) 一种变偏移距垂直地震剖面数据反演各向异性参数方法
CN103630934B (zh) 一种确定转换波检波点大的横波静校正量的方法
CN104614765A (zh) 一种增强地震波激发照明的设计方法
CN102213769A (zh) 一种利用三维垂直地震剖面资料确定各向异性参数的方法
CN105093292A (zh) 一种地震成像的数据处理方法和装置
CN101738636B (zh) 一种三维vsp高斯束法多波联合偏移成像方法
CN104570072A (zh) 一种粘弹性介质中的球面pp波反射系数建模方法
CN101609165B (zh) 一种基于起伏地表的折射波静校正方法
CN104678434A (zh) 一种预测储层裂缝发育参数的方法
Xu et al. The horizontal-to-vertical spectral ratio and its applications
CN104316965A (zh) 一种裂缝方位和强度的预测方法及系统
CN106646661A (zh) 一种矿床水文地质综合勘查系统
CN109669223A (zh) 利用微测井获得浅地表声波测井与合成地震记录的方法
Kaklamanos et al. Introduction to the special section on advances in site response estimation
Ning et al. High-frequency surface-wave imaging from traffic-induced noise by selecting in-line sources
CN103675900A (zh) 一种确定转换波叠前时间偏移最佳速度剖面的方法
Tran et al. An assessment of surface wave techniques at the Texas A&M national geotechnical experimentation site
CN103558637A (zh) 基于三分量传感器的远探测方法
Murphy et al. Detailed P-and S-wave velocity models along the LARSE II transect, southern California
CN103472482A (zh) 基于基因排序体系的多波地震资料时间域匹配方法及系统
CN104101897A (zh) 一种利用矢量合成纵波进行地震勘探的方法
CN101609166B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20141015

RJ01 Rejection of invention patent application after publication