CN101598811B - 一种二维垂直地震剖面数据计算炮点静校正的方法 - Google Patents
一种二维垂直地震剖面数据计算炮点静校正的方法 Download PDFInfo
- Publication number
- CN101598811B CN101598811B CN2008101143789A CN200810114378A CN101598811B CN 101598811 B CN101598811 B CN 101598811B CN 2008101143789 A CN2008101143789 A CN 2008101143789A CN 200810114378 A CN200810114378 A CN 200810114378A CN 101598811 B CN101598811 B CN 101598811B
- Authority
- CN
- China
- Prior art keywords
- data
- vertical seismic
- static correction
- vsp
- 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.)
- Active
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是地球物理勘探二维垂直地震剖面数据计算炮点静校正的方法。对零井源距垂直地震剖面地震数据进行初至拾取,求取观测点深度处的层速度,拾取二维垂直地震剖面的纵波直达波初至,对二维垂直地震剖面原始数据进行三分量偏振分析,旋转处理,然后反演计算二维垂直地震剖面数据每一炮的静校正量,将每一道地震数据点的振幅值移动完成静校正。本发明不需要进行野外近地表调查,并且克服了野外调查点稀疏,所求取的静校正精度不够的不利之处,克服了野外近地表调查方法调查深度较浅的缺点,提高了垂直地震剖面的处理质量。
Description
技术领域
本发明涉及地球物理勘探技术,具体是一种二维垂直地震剖面数据计算炮点静校正的方法。
背景技术
常规的垂直地震剖面(沿设计的一条或多条过井的炮线进行地面激发,用井中检波器进行记录的垂直地震剖面方法,也称作walkaway VSP,简称为WVSP)数据是在地面单点激发,井下多点接收得到的地震数据;为了扩大垂直地震剖面数据的覆盖范围与提高覆盖次数,设计了二维垂直地震剖面(WVSP)技术,该技术是在地面沿设计的过观测井的炮线上连续多点激发,而在井下多点接收得到的地震数据。在二维垂直地震剖面数据采集时,由于不同炮点激发位置的地表起伏与地下地质结构的差异,导致不同炮之间接收的地震数据有静校正差异,在成像时各炮数据不能同相叠加,降低了最终的数据成果质量。通常的做法是利用近地表调查方法来消除这种静校正差异,即调查得到调查点的高程、风化层厚度与相应的速度,然后用统一的高程与速度进行替换,得到替换前后的时差作为静校正量,在炮点与调查点不一致时采用插值的方法估算得到静校正量。但是这存在两个问题:一是近地表调查得到的数据对应的深度相比二维垂直地震剖面的接收点深度浅很多;二是近地表调查点的间距比二维垂直地震剖面的炮点间距大很多,不能每个炮点都有准确的调查结果。这就导致近地表调查计算炮点静校正的方法不能满足二维垂直地震剖面技术的需要。二维垂直地震剖面技术没有相适应的炮点静校正技术来提高二维垂直地震剖面的处理质量。
发明内容
本发明目的是提供一种提高垂直地震剖面的处理质量的二维垂直地震剖面数据计算炮点静校正的方法。
本发明具体是采用了以下技术步骤:
1)采集地震数据;
步骤1)所述的采集地震数据是用在地面炮线的不同位置分别人工激发地震纵波,在井中不同深度放置三分量地震波接收器接收地震波场记录三分量二维地震波波场与常规的零井源距垂直地震剖面数据。
步骤1)所述的地震波接收器两两相互垂直。
2)对零井源距垂直地震剖面(VSP)地震数据进行初至拾取,得到时间-深度数据对(t-d),然后利用以下公式求取观测点深度处的层速度:
式中:Vi为观测点深度处的层速度,X为激发点与观测井之间的水平距离,即偏移距;H为从地面面起算的观测点深度;T为拾取的直达纵波的初至时间;
步骤2)所述的层速度计算如果没有零井源距数据,则用二维垂直地震剖面中的小井源距数据。
3)根据观测点层速度计算每一层沿垂直方向传播的纵波层速,反演出Vn,Vn=(v1,v2,…,vn),分层数为n;
4)对二维垂直地震剖面数据进行常规的三分量偏振分析,旋转处理,然后拾取二维垂直地震剖面的纵波直达波初至;
5)对于n层的模型,定义每一层的各向异性系数为:
6)对二维垂直地震剖面采集到的三分量原始数据进行常规的三分量偏振分析,旋转处理,然后反演计算二维垂直地震剖面数据每一炮的静校正量:
步骤6)所述的计算得到二维垂直地震剖面静校正量利用以下公式:
式中:M为垂直地震剖面接收点个数,为第j个垂直地震剖面初至时距曲线、第m个垂直地震剖面接收点实际观测的初至时间,Tj,m(V,αn)为由速度模型求出的第j个垂直地震剖面初至时距曲线、第m个垂直地震剖面接收点的理论初至时间,statj为第j个炮点的静校正量。
7)将每一道地震数据第i点的振幅值移动到第(i-stat)点上,完成静校正。
本发明不需要进行野外近地表调查,并且克服了野外调查点稀疏,所求取的静校正精度不够的不利之处,克服了野外近地表调查方法调查深度较浅(最多几十米)的缺点,提高了垂直地震剖面的处理质量。
附图说明:
图1是三分量(X、Y、Z)原始数据(单炮点点)的波场图;
图2是常规零井源距数据处理得到的初至与速度分层图;
图3是计算得到的各炮点静校正量图;
图4(上)是静校正前的共检波点纪录波场图;
图4(下)是静校正后的共检波点纪录波场图;
图5(上)是经过进行静校正前的成像剖面图。
图5(下)是经过进行静校正后的成像剖面图。
具体实施方式
野外通过机械可控震源在地面设计的过井炮线上人工激发地震波,在井中放置三分量地震波串接收器接受地震波场信号,再由电缆传到地面记录仪器记录三分量垂直地震剖面地震波波场,分别得到零井源距垂直地震剖面和二维垂直地震剖面数据;
对采集到的三分量零井源距垂直地震剖面原始地震数据的垂直分量(图1)进行初至拾取,得到激发炮点在井下所有检波器记录数据的地震波初至时间Ti,利用公式(1)求取每一个井下观测点层速度,然后对层速度进行地质分层得到层速度模型;
对二维垂直地震剖面(WVSP)采集到的三分量原始数据进行常规的三分量检波器偏振分析,检波器旋转校正处理,可以得到直达地震波的分量数据,然后在该数据上逐炮拾取二维垂直地震剖面的纵波直达波初至,得到所有地面激发炮点在井下所有检波器记录数据的地震波初至时间Tj,m(第j个垂直地震剖面初至时距曲线、第m个垂直地震剖面接收点实际观测的初至时间),对于图1所示的多层介质模型,根据S.I.Lapin(1978)提出的公式(3)可以求出射线参数p:
其中,L是炮点的井源距,N是速度模型的层数,hn为模型第n层的厚度。通过(4)式求出垂直地震剖面初至旅行时:
对于有J个炮点的二维垂直地震剖面记录,我们实际观测的垂直地震剖面初至时距曲线有J个,则可以构建如下的目标函数R(Vn,An):
其中,M为垂直地震剖面接收点个数,为第j个垂直地震剖面初至时距曲线、第m个垂直地震剖面接收点实际观测的初至时间,Tj,m(Vn,αn)为由速度模型求出的第j个垂直地震剖面初至时距曲线、第m个垂直地震剖面接收点的理论初至时间,statj为第j个炮点的静校正量。
由(5)式,令R(Vn,An)→min,通过阻尼最小二乘法(Levenberg-Marquardt算法)我们可求得每一个炮点的静校正量,并同时得到与该静校正量匹配的沿垂直方向传播的纵波层速度Vn=(v1,v2,…,vn)和各向异性系数An=(a1,a2,…,an)。
通过公式(5)计算得到二维垂直地震剖面数据每一炮的静校正量stat,计算得到的每一炮的炮点静校正如图3所示。对每一炮数据应用该静校正量,即将每一道地震数据第i点的振幅值移动到第(i-stat)点上。完成静校正的应用,静校正应用前后的共检波点数据对比如图4所示。
对应用静校正前后的数据分别进行常规处理(波场分离,反褶积,成像等),最后分别得到成像后的数据,如图5所示。
本发明实施例用采集到的常规零井源距垂直地震剖面和二维垂直地震剖面,进行本发明的炮点静校正后,成像剖面上扭曲的同相轴得到了明显改善,通过图4和图5的对比可见,本发明静校正后的剖面质量明显提高。
Claims (3)
1.一种二维垂直地震剖面数据计算炮点静校正的方法,其特征在于采用以下技术步骤:
1)采集地震数据,在地面过井炮线的不同位置分别人工激发地震纵波,在井中不同深度放置三分量地震波接收器接收地震波场记录三分量二维地震波波场与常规的零井源距垂直地震剖面数据;
2)对零井源距垂直地震剖面(VSP)地震数据进行初至拾取,得到时间-深度数据对(t-d),然后利用以下公式求取观测点深度处的层速度:
式中:Vi为观测点深度处的层速度,X为激发点与观测井之间的水平距离,即偏移距;H为从地面起算的观测点深度;T为拾取的直达纵波的初至时间;
3)根据观测点层速度计算每一层沿垂直方向传播的纵波层速,反演出Vn,Vn=(v1,v2,…,vn),分层数为n;
4)对二维垂直地震剖面数据进行常规的三分量偏振分析,旋转处理,然后拾取二维垂直地震剖面的纵波直达波初至;
5)对于n层的模型,定义每一层的各向异性系数为:
6)对二维垂直地震剖面采集到的三分量原始数据进行常规的三分量偏振分析,旋转处理,然后反演计算二维垂直地震剖面数据每一炮的静校正量:
计算利用以下公式:
式中:M为垂直地震剖面接收点个数,为第j个垂直地震剖面初至时距曲线、第m个垂直地震剖面接收点实际观测的初至时间,Tj,m(Vn,αn)为由速度模型求出的第j个垂直地震剖面初至时距曲线、第m个垂直地震剖面接收点的理论初至时间,statj为第j个炮点的静校正量;
7)将每一道地震数据第i点的振幅值移动到第i-stat点上,完成静校正。
2.根据权利要求1所述的二维垂直地震剖面数据计算炮点静校正的方法,其特征在于步骤1)所述的地震波接收器两两相互垂直。
3.根据权利要求1所述的二维垂直地震剖面数据计算炮点静校正的方法,其特征在于步骤2)所述的层速度计算如果没有零井源距数据,则用二维垂直地震剖面中的小井源距数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101143789A CN101598811B (zh) | 2008-06-04 | 2008-06-04 | 一种二维垂直地震剖面数据计算炮点静校正的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101143789A CN101598811B (zh) | 2008-06-04 | 2008-06-04 | 一种二维垂直地震剖面数据计算炮点静校正的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101598811A CN101598811A (zh) | 2009-12-09 |
CN101598811B true CN101598811B (zh) | 2011-04-20 |
Family
ID=41420314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101143789A Active CN101598811B (zh) | 2008-06-04 | 2008-06-04 | 一种二维垂直地震剖面数据计算炮点静校正的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101598811B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102213769A (zh) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 |
CN102466821B (zh) * | 2010-11-18 | 2014-01-15 | 中国石油天然气集团公司 | 一种利用局部爆炸反射面地震剖面评价静校正的方法 |
CN102478666B (zh) * | 2010-11-23 | 2013-11-13 | 中国石油天然气集团公司 | 一种提高可控震源地震记录初至拾取精度的方法 |
CN103105624B (zh) * | 2011-11-11 | 2015-07-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于数据库技术的纵横波时差定位方法 |
CN104570119B (zh) * | 2013-10-29 | 2017-11-17 | 中国石油化工股份有限公司 | 一种三维垂直地震剖面反射波拉伸校正方法 |
CN107843928A (zh) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | 一种基于网格搜索算法的炮点纠偏方法 |
CN109143358B (zh) * | 2018-08-31 | 2019-09-13 | 中国石油大学(北京) | 一种获得深层、超深层碎屑岩地层压力结构的方法及装置 |
CN109991662A (zh) * | 2019-05-15 | 2019-07-09 | 中油奥博(成都)科技有限公司 | 浅部地层二维或三维弹性参数测量和计算的装置及方法 |
CN112305601B (zh) * | 2019-07-31 | 2024-05-17 | 中国石油化工股份有限公司 | 基于层析反演的地震勘探采集激发井的设计方法 |
CN111650645B (zh) * | 2020-06-16 | 2022-05-20 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp弯线校正处理方法和装置 |
CN112462426B (zh) * | 2020-11-02 | 2024-05-28 | 中国石油天然气集团有限公司 | 横波矢量静校正方法及装置 |
CN113917539B (zh) * | 2021-10-15 | 2023-07-14 | 中国地质科学院地球物理地球化学勘查研究所 | 火山岩覆盖区地震数据叠前道集处理方法、系统和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1308240A (zh) * | 2001-01-21 | 2001-08-15 | 周熙襄 | 用初至折射波对地震勘探资料进行短波长静校正处理的方法 |
CN1365008A (zh) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
CN1523373A (zh) * | 2003-02-20 | 2004-08-25 | 中国石油集团东方地球物理勘探有限责 | 上行透射转换波计算横波静校正的方法 |
CN1580817A (zh) * | 2003-08-01 | 2005-02-16 | 中国石油天然气集团公司 | 分步静校正计算方法 |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
-
2008
- 2008-06-04 CN CN2008101143789A patent/CN101598811B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1365008A (zh) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
CN1308240A (zh) * | 2001-01-21 | 2001-08-15 | 周熙襄 | 用初至折射波对地震勘探资料进行短波长静校正处理的方法 |
CN1523373A (zh) * | 2003-02-20 | 2004-08-25 | 中国石油集团东方地球物理勘探有限责 | 上行透射转换波计算横波静校正的方法 |
CN1580817A (zh) * | 2003-08-01 | 2005-02-16 | 中国石油天然气集团公司 | 分步静校正计算方法 |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
Non-Patent Citations (1)
Title |
---|
王宝彬等.利用VSP资料估算反射界面倾角方法综述.《勘探地球物理进展》.2008,第31卷(第1期),期刊第16-25页. * |
Also Published As
Publication number | Publication date |
---|---|
CN101598811A (zh) | 2009-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101598811B (zh) | 一种二维垂直地震剖面数据计算炮点静校正的方法 | |
US10577926B2 (en) | Detecting sub-terranean structures | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
US7359283B2 (en) | System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers | |
CN102213769A (zh) | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 | |
CN101551463B (zh) | 三维观测系统噪声压制估算方法 | |
CN102053261B (zh) | 一种地震数据处理方法 | |
CN102540250B (zh) | 基于方位保真角度域成像的裂缝型油气储层地震探测方法 | |
CN101354444B (zh) | 一种确定地层岩性和孔隙流体的方法 | |
CN103630934B (zh) | 一种确定转换波检波点大的横波静校正量的方法 | |
CN103026265A (zh) | 分开波型的地震采集方法 | |
CN102375154A (zh) | 基于宽方位三维地震的裂缝参数确定方法 | |
CN102590862A (zh) | 补偿吸收衰减的叠前时间偏移方法 | |
CN102053263A (zh) | 调查表层结构的方法 | |
CN102053260B (zh) | 获得地震纵波的方位速度的方法及处理地震数据的方法 | |
Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
CN102053262B (zh) | 获得地震转换波的方位速度的方法及处理地震数据的方法 | |
CN103149588A (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN101609163B (zh) | 基于波动理论的多尺度地震资料联合成像方法 | |
CN102798888B (zh) | 一种利用非零井源距数据计算纵横波速度比的方法 | |
CN101900830B (zh) | 一种地震散射p-p波成像方法 | |
Orlando et al. | P and PS data to reduce the uncertainty in the reconstruction of near-surface alluvial deposits (Case study—Central Italy) | |
Holzschuh | 13GA–EG1 Eucla-Gawler seismic survey–acquisition and processing of the western Gawler craton section | |
Raji | Reflection Enhancement Processing in Ultra-Shallow Seismic Data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |