CN103823241B - 一种在炮检距域中求取折射剩余静校正量的方法 - Google Patents
一种在炮检距域中求取折射剩余静校正量的方法 Download PDFInfo
- Publication number
- CN103823241B CN103823241B CN201210461796.1A CN201210461796A CN103823241B CN 103823241 B CN103823241 B CN 103823241B CN 201210461796 A CN201210461796 A CN 201210461796A CN 103823241 B CN103823241 B CN 103823241B
- Authority
- CN
- China
- Prior art keywords
- value
- static correction
- residual
- point
- residual static
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是地球物理勘探中利用初至时间在基础静校正量的基础上求取折射剩余静校正量的方法。用基础静校正量计算所有炮检距的激发点或接收点在多个炮检距中的不同的剩余静校正量,计算分解炮检对剩余静校正量,计算所有炮检对剩余校正量的均方根误差;计算一个激发点和接收点在共炮检距域中的最终剩余静校正量,用激发点和接收点的剩余静校正量到已知的基础静校正量上,重复迭代直到均方根误差值达到给定范围完成共炮检距域中折射剩余静校正。本发明避免了常规在共激发点域、共接收点域和共中心点域计算剩余静校正量的不足,提高了短波长静校正的精度,改善长波长静校正的精度,适应二维勘探和三维勘探。
Description
技术领域
本发明涉及地球物理勘探静校正方法,是一种利用初至时间,在基础静校正量的基础上求取折射剩余静校正量的方法。
背景技术
地震勘探中静校正处理主要是解决地震波在近地表介质旅行时间变化影响反射波叠加成像问题。初至波反映近地表地层变化,利用初至波的折射静校正方法或层析静校正方法可以反演出相对精确的近地表模型,解决大部分长波长和中长波长的静校正问题,在解决短波长静校正方面,还存在明显不足。短波长静校正直接影响地震剖面成像的质量,是非常重要的。为了进一步提高短波长静校正量的精度,大都采取了求取剩余静校正量的方法。从使用数据的类型和原理上来分,常规求取剩余静校正量的方法有两个:一个是利用反射波信息和最优叠加反射波同相轴的原理求取地震道的剩余静校正量,再分解该地震道的剩余静校正量,得到激发点和接收点的剩余静校正量。该方法需要精度较高的叠加速度场,而得到精度较高的叠加速度场本身就是一个比较困难的问题,尤其在复杂的近地表区域,静校正量严重影响叠加速度的分析结果。叠加速度分析又是一个交互迭代过程,所以用该方法计算剩余静校正量的速度较慢;另一个方法是利用初至信息或折射波信息,在共激发点域、共接收点域和共中点域中,利用拟合初至波成一条光滑的曲线原理来求取激发点和接收点的剩余静校正量。由于使用的数据仅仅是初至和静校正量,所以该方法计算速度较快。但是该方法要求初至是连续的,在应用静校正量后,相邻地震道的初至不能出现较大跳动。为了避免相同炮检距的初至来自不同的层或折射界面倾斜的情况,在计算中要求对初至进行区分左右支,对每支的初至进行多项式拟求取剩余静校正量。该方法是在共激发点域、共接收点域或共中点域中进行的,比较适合二维勘探,在三维勘探中效果不明显。
发明内容
本发明目的在于提供一种利用初至信息,适合二维勘探和三维勘探,进一步提高短波长静校正量的精度,快速求取在炮检距域中求取折射剩余静校正量的方法。
本发明通过以下步骤实现:
1)采集地震数据,抽取共炮检距道集G:
所述的抽取共炮检距道集是在炮检距D和浮动范围d以内,抽取所有初至和相应炮检对的中点坐标,形成一个共炮检距道集G。
2)应用基础静校正量:
所述的应用基础静校正量是在共炮检距道集G中,把已知激发点和接收点基础静校正量的和,应用到初至上,形成时间t1。
3)计算炮检对剩余静校正量:
所述的计算炮检对剩余静校正量是在炮检距道集G中,以一个炮检对的中点坐标为圆心,以平滑半径R形成一个圆,计算道集G中所有炮检对中点坐标在该圆中的时间t1的平均值,形成时间t2,时间t1与时间t2时差就是该炮检对的剩余静校正量dt;
所述的时间t2利用以下公式计算:
式中:t1i是指在道集G中,其它炮检对中点坐标到当前炮检对中点坐标的距离小于平滑半径R的时间t1;
N是相应的炮检对的总个数;i是序号;
所述的平滑半径R取最大炮检距的一半;
所述的平滑半径R取1000-2000m;
4)分解炮检对剩余校正量:
所述的分解炮检对剩余校正量是:按照激发点的比例系数a和接收点的比例系数,分解炮检对的剩余校正量dt到相应的激发点ds和接收点dr中。
步骤4)中所述比例系数a取值范围0至1;相应的接收点的比例系数是1减a,取值范围是1至0。
步骤4)中所述激发点ds=dt*a;接收点dr=dt*(1-a)。
5)采用另一个炮检距D+2d,重复迭代步骤1)到4),直到计算完所有炮检距;得到激发点或接收点在多个炮检距中的不同的剩余静校正量;
6)利用以下公式计算所有炮检对剩余校正量dt的均方根误差rms:
式中:N是步骤1)到5)中所有道集G中所有炮检对的总个数;i是序号;
7)利用以下公式计算一个激发点dsi在共炮检距域中的最终剩余静校正量ts:
式中:Ns是同一个激发点在步骤1)到4)中产生剩余静校正量的总个数;i是序号。
利用以下公式计算一个接收点在共炮检距域中的最终剩余静校正量tr:
式中:Nr是同一个接收点在步骤1)到4)中产生剩余静校正量的总个数;i是序号;
8)用激发点和接收点的剩余静校正量到已知的基础静校正量上,重复迭代步骤1)到7),直到前后两次迭代的rms差值达到给定范围,完成在共炮检距域中的折射剩余静校正量。
本发明在共炮检距域中计算剩余静校正量,由于使用坐标计算距离和确定位置,避免了常规在共激发点域、共接收点域和共中心点域计算剩余静校正量的不足,所以同时适应二维勘探和三维勘探。在该域求取剩余静校正量应用后,再在共激发点或共接收点或共中心点域计算剩余静校正量,结果在后三个域中的剩余静校正量较小,可以忽略不计,所以在共炮检距域中计算完剩余静校正量后,可以不需要在其它域中计算剩余静校正量了。并且在炮检距域中计算剩余静校正量,不仅提高了短波长静校正的精度,而且还会改善长波长静校正的精度。
具体实施方式
在共炮检距域中求取折射剩余静校正量的方法,以道距10米、计算炮检距200m到2000m为例,有以下步骤实现:
1)采集地震数据,抽取共炮检距道集:给出一个炮检距(激发点到接收点之间的距离)D=200和浮动范围d=10,抽取所有炮检距在范围[190,210)中的初至picker、和相应炮检对(一个激发点和一个接收点)的中点坐标(x,y),形成一个共炮检距道集G。
2)应用基础静校正量:在道集G中,把已知激发点和接收点基础静校正量的和,应用到初至上,形成时间t1。
3)计算炮检对剩余静校正量:在道集G中,以一个炮检对的中点坐标(x,y)为圆心,以给出的平滑半径R(如1000m)为半径形成一个圆,计算道集G中所有炮检对中点坐标在该圆中的t1的平均值,形成时间t2,时差dt(t2-t1)就是该炮检对的剩余静校正量。利用公式(1)计算时间t2:
公式(1)中t1是指在道集G中,其它炮检对中点坐标到当前炮检对中点坐标的距离小于给定半径R的t1,N是相应的炮检对的总个数。
4)分解炮检对剩余校正量:给出激发点的比例系数和接收点的比例系数,分解炮检对的剩余校正dt到相应的激发点ds和接收点dr中。
步骤4)中所述激发点的比例系数是指用户给定的系数a(如0.5),取值范围[0,1];相应的接收点的比例系数是1-a,取值范围是[1,0]。
步骤4)中所述激发点ds=dt*a;接收点dr=dt*(1-a)。
5)重复下一个炮检距:给出下一个炮检距D+2d,重复迭代步骤1)到4)。直到计算完所有炮检距。
6)计算dt的均方根误差。利用公式(2)计算均方根误差rms:
公式(2)中N是步骤1)到5)中所有道集G中所有炮检对的总个数。
7)计算激发点和接收点的剩余静校正量:
利用公式(3)计算一个激发点在共炮检距域中的最终剩余静校正量:
公式(3)中,ds是步骤4)产生的激发点剩余静校正量,N是该量的总个数。
利用公式(4)计算一个接收点在共炮检距域中的最终剩余静校正量:
公式(4)中,dr是步骤4)产生的接收点剩余静校正量,N是该量的总个数。
8)应用激发点和接收点的剩余静校正量到已知的基础静校正量上,重复迭代步骤1)到7)。直到迭代的次数达到了用户给定的次数(4次)或前后两次迭代的rms差值达到了用户给定的范围(如0.1),就完成了在共炮检距域中的折射剩余静校正量的工作。
Claims (7)
1.一种在炮检距域中求取折射剩余静校正量的方法,特点是通过以下步骤实现:
1)采集地震数据,抽取共炮检距道集G:
所述的抽取共炮检距道集G是在炮检距D和浮动范围d以内,抽取所有初至和相应炮检对的中点坐标生成;
2)应用基础静校正量:
3)计算炮检对剩余静校正量:
所述的计算炮检对剩余静校正量是在炮检距道集G中,以一个炮检对的中点坐标为圆心,以平滑半径R形成一个圆,计算道集G中所有炮检对中点坐标在该圆中的时间t1的平均值,形成时间t2,时间t1与时间t2时差就是该炮检对的剩余静校正量dt;
所述的时间t2利用以下公式计算:
式中:t1i是指在道集G中,其它炮检对中点坐标到当前炮检对中点坐标的距离小于平滑半径R的时间t1;
N是相应的炮检对的总个数;i是序号;
4)分解炮检对剩余校正量:
5)采用另一个炮检距重复迭代步骤1)到4),直到计算完所有炮检距;得到激发点或接收点在多个炮检距中的不同的剩余静校正量;
所述的另一个炮检距是炮检距D加两倍的浮动范围d;
6)利用以下公式计算所有炮检对剩余校正量dt的均方根误差rms:
式中:N是步骤1)到5)中所有道集G中所有炮检对的总个数;i是序号;
7)利用以下公式计算一个激发点在共炮检距域中的最终剩余静校正量ts:
式中:Ns是同一个激发点在步骤1)到4)中产生剩余静校正量的总个数;dsi是激发点剩余静校正量;i是序号;
利用以下公式计算一个接收点在共炮检距域中的最终剩余静校正量tr:
式中:Nr是同一个接收点在步骤1)到4)中产生剩余静校正量的总个数;dri是接收点剩余静校正量;i是序号;
8)应用激发点和接收点的剩余静校正量叠加到已知的基础静校正量上,重复迭代步骤1)到7),直到前后两次迭代的rms差值达到给定范围,完成在共炮检距域中的折射剩余静校正量。
2.根据权利要求1的方法,特点是步骤2)所述的应用基础静校正量是在共炮检距道集G中,把已知激发点和接收点基础静校正量的和,应用到初至上,形成时间t1。
3.根据权利要求1的方法,特点是步骤3)所述的平滑半径R取最大炮检距的一半。
4.根据权利要求1的方法,特点是步骤3)所述的平滑半径R取1000-2000m。
5.根据权利要求1的方法,特点是步骤4)所述的分解炮检对剩余校正量是:按照激发点的比例系数a和接收点的比例系数,分解炮检对的剩余校正量dt到相应的激发点剩余静校正量ds和接收点剩余静校正量dr中。
6.根据权利要求5的方法,特点是所述比例系数a取值范围0至1;相应的接收点的比例系数是1减a,取值范围是1至0。
7.根据权利要求5的方法,特点是所述激发点剩余静校正量ds=dt*a;接收点剩余静校正量dr=dt*(1-a)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210461796.1A CN103823241B (zh) | 2012-11-16 | 2012-11-16 | 一种在炮检距域中求取折射剩余静校正量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210461796.1A CN103823241B (zh) | 2012-11-16 | 2012-11-16 | 一种在炮检距域中求取折射剩余静校正量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103823241A CN103823241A (zh) | 2014-05-28 |
CN103823241B true CN103823241B (zh) | 2016-08-10 |
Family
ID=50758389
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210461796.1A Active CN103823241B (zh) | 2012-11-16 | 2012-11-16 | 一种在炮检距域中求取折射剩余静校正量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103823241B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104122587B (zh) * | 2014-06-30 | 2016-10-19 | 中国石油天然气集团公司 | 一种基于炮检距域的异常初至识别方法以及*** |
CN104197777B (zh) * | 2014-08-18 | 2017-02-08 | 中国人民解放军63853部队 | 火炮回转半径测量方法 |
CN105911592B (zh) * | 2016-06-30 | 2017-11-10 | 中国石油化工股份有限公司 | 一种实钻约束的三维地震资料长波长静校正方法 |
CN107918154B (zh) * | 2016-10-10 | 2019-06-04 | 中国石油化工股份有限公司 | 一种地表一致性的基准面静校正量的组合方法 |
CN111624655B (zh) * | 2019-02-27 | 2023-02-07 | 中国石油天然气集团有限公司 | 初至波剩余静校正量的确定方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2254146B (en) * | 1991-03-08 | 1994-12-07 | Geco As | Determining residual static shift in seismic data |
US6188964B1 (en) * | 1999-09-14 | 2001-02-13 | Ut-Battelle, Llc | Method for using global optimization to the estimation of surface-consistent residual statics |
CN1365008A (zh) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
CN102692648A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种基于大炮初至的折射波剩余静校正方法 |
-
2012
- 2012-11-16 CN CN201210461796.1A patent/CN103823241B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2254146B (en) * | 1991-03-08 | 1994-12-07 | Geco As | Determining residual static shift in seismic data |
US6188964B1 (en) * | 1999-09-14 | 2001-02-13 | Ut-Battelle, Llc | Method for using global optimization to the estimation of surface-consistent residual statics |
CN1365008A (zh) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
CN102692648A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种基于大炮初至的折射波剩余静校正方法 |
Non-Patent Citations (2)
Title |
---|
交互迭代静校正方法;王彦春 等;《石油物探》;19980630;第37卷(第2期);第63-70页 * |
多域迭代折射静校正及其应用;王新红 等;《石油物探》;20011231;第40卷(第4期);第56-61页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103823241A (zh) | 2014-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103823241B (zh) | 一种在炮检距域中求取折射剩余静校正量的方法 | |
US7630865B2 (en) | Complex analysis of kinematics for non-hyperbolic moveout corrections | |
Dong et al. | Theory and practice of refraction interferometry | |
Zhang et al. | Common-image gathers in the incident phase-angle domain from reverse time migration in 2D elastic VTI media | |
US20110096627A1 (en) | Method and system for seismic imaging and earth modeling using beam tomography | |
CN104656142A (zh) | 一种利用垂直地震剖面与测井联合的地震层位标定方法 | |
CN108363101A (zh) | 一种斜井井间地震高斯束叠前深度偏移成像方法 | |
NO180026B (no) | Fremgangsmåte for korreksjon av amplitudeverdier ved seismiske målinger | |
CN104483705A (zh) | 一种三维剩余静校正方法 | |
Purdy | The correction for the travel time effects of seafloor topography in the interpretation of marine seismic data | |
CN106199704A (zh) | 一种三维三分量海底电缆地震资料速度建模方法 | |
WO2015053659A1 (ru) | Способ получения априорного годографа для выполнения литолого стратиграфической привязки | |
Shen et al. | Improving subsalt imaging with illumination-based weighting of RTM 3D angle gathers | |
AU2004286575B2 (en) | Kirchhoff prestack time migration method for PS waves | |
Gong et al. | Combined migration velocity model-building and its application in tunnel seismic prediction | |
CN101963674A (zh) | 共转换点抽道集方法 | |
CN106932821B (zh) | 地震层析反演中的一种目标射线追踪方法 | |
CN106338760B (zh) | 误差补偿的起伏地表偏移方法 | |
CN105510958A (zh) | 一种适用复杂介质的三维vsp观测***设计方法 | |
Cormier et al. | Calculation of strong ground motion due to an extended earthquake source in a laterally varying structure | |
CN104101901A (zh) | 一种转换波弯曲射线保幅各向异性叠前时间偏移方法 | |
Cantillo et al. | Ocean bottom nodes (OBN) repeatability and 4D | |
Alshuhail et al. | Near-surface characterization using vertical array seismic data from smart DAS upholes | |
Yang et al. | Least-squares reverse time migration with velocity errors | |
CN105954798A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |