CN112099082A - 一种共面元共方位角道集的地震回折波走时反演方法 - Google Patents
一种共面元共方位角道集的地震回折波走时反演方法 Download PDFInfo
- Publication number
- CN112099082A CN112099082A CN201910520239.4A CN201910520239A CN112099082A CN 112099082 A CN112099082 A CN 112099082A CN 201910520239 A CN201910520239 A CN 201910520239A CN 112099082 A CN112099082 A CN 112099082A
- Authority
- CN
- China
- Prior art keywords
- velocity
- gather
- travel time
- inflection
- seismic
- 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 48
- 238000004364 calculation method Methods 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 9
- 230000008859 change Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000032798 delamination Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Images
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
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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
技术领域
本发明涉及一种基于共面元共方位角道集地震回折波走时反演的三维起伏地表速度建模方法,属于地震勘探走时反演技术领域。
背景技术
地震勘探走时反演技术利用地震波走时建立地下介质的速度模型。通常分为利用地震初至波走时的初至走时反演方法和利用地震反射波走时的反射波走时反演方法。较反射波而言,初至波的初至走时往往易于识别和拾取,这使得初至走时反演的应用范围更加广泛。
常规的初至走时反演通过最小化给定模型的观测初至走时与实际走时的残差得到最终反演结果,这需要反复迭代,而且每次迭代都要进行大量的正演计算,因而其计算效率较低。在此基础上,通过线性求解射线方程和走时方程重建地下速度结构的初至走时反演方法受到众多学者的青睐,如折射波走时反演和回折波走时反演等。但折射波走时反演需要地下介质存在良好的折射界面,而基于连续介质假设的回折波走时反演更加符合实际情况,故而其具有广阔的应用前景。基于此,更有学者提出了利用回折波走时快速建立速度模型的剥层法,该方法极大地提高了速度建模的效率,但该方法仍利用基于水平地表假设的回折波走时方程,在近地表条件复杂的地区仍难以取得较好的效果。
目前现有的回折波走时速度建模方法均基于地表高程变化平缓和速度横向变化均匀的假设条件,该类方法在近地表条件简单的地区通常可以取得较好的应用效果。然而,该类方法没有考虑到剧烈起伏的地形对反演结果的影响,这使得现有的回折波走时速度建模方法在地形起伏剧烈的地区的应用效果较差。
发明内容
针对现有技术存在的上述缺陷,本发明提出了一种共面元共方位角道集的地震回折波走时反演方法,其解决了现有的回折波走时速度建模方法不适用于地表起伏剧烈的近地表条件复杂地区的问题。
本发明是采用以下的技术方案实现的,本发明所述一种共面元共方位角道集的地震回折波走时反演方法,包括如下步骤:
步骤一:获取回折波走时数据;
步骤二:计算地表速度;
步骤三:获取共面元共方位角道集并编号;
步骤四:抽取当前共面元共方位角道集;
步骤五:将当前共面元共方位角道集内的地震道按照偏移距从小到大进行排序;
步骤六:逐道计算回折波射线参数和速度梯度;利用二分法求解如下方程计算回折波射线参数:
其中,x为偏移距,t为回折波走时,i为每一层编号,n为当前层编号,vs,i-1为第i层炮点端速度,vr,i-1为第i层检波点端速度,pn为回折波射线参数。特别地,当i=1时,vs,i-1即vs,0表示炮点处地表速度,vr,i-1即vr,0表示检波点处地表速度;
速度梯度计算公式如下:
其中,x为偏移距,i为每一层编号,n为当前层编号,vs,i-1为第i层炮点端速度,vr,i-1为第i层检波点端速度,vs,n-1为当前层炮点端速度,vr,n-1为当前层检波点端速度,pn为回折波射线参数,bi为第i层速度梯度,bn为当前层速度梯度。特别地,当i=1时,vs,i-1即vs,0表示炮点处地表速度,vr,i-1即vr,0表示检波点处地表速度;当n=1时,vs,n-1即vs,0表示炮点处地表速度,vr,n-1即vr,0表示检波点处地表速度;
步骤七:逐道计算回折点速度和深度。回折点速度计算公式如下:
其中,n为当前层编号,pn为回折波射线参数,vn为当前层回折点处速度;
回折点深度计算公式如下:
其中,vn为当前层回折点速度,vn-1为当前层回折点所对应的顶界面速度,bn为当前层速度梯度,zn-1为当前层回折点所对应的顶界面的深度,zn为当前层回折点所对应的深度。特别地,当n=1时,vn-1即v0表示面元中点处对应的地表速度,zn-1即z0表示面元中点处对应的深度;
步骤八:判断是否有未经处理的共面元共方位角道集,若有,则将该未经处理的共面元共方位角道集设置为当前共面元共方位角道集,返回步骤六;若无,执行下一步骤;
步骤九:利用所有已计算的回折点速度和深度,通过反距离加权插值算法计算离散速度模型所有网格点速度;
步骤十:输出速度模型。
进一步地,在步骤一所述回折波走时数据中,应包括炮点坐标、炮点井深、检波点坐标、回折波初至走时。
本发明的有益效果是:采用本发明所述的一种共面元共方位角道集的地震回折波走时反演方法,本发明基于炮检点分解原理,利用每一个炮检对的绝对偏移距和绝对走时计算回折波射线参数和速度梯度,其不受相邻道的影响,从而避免了地表起伏对回折波射线参数计算的影响,以及炮检点高程差异大使得回折点深度计算不准等问题,其克服了现有回折波走时速度建模方法难以适用于近地表条件复杂地区的问题;本发明计算简单、易于实现、耗时少,具有很强的适应性,且反演结果可靠性高。
附图说明
图1为本发明的流程图。
图2为本发明测试所用理论速度模型图。
图3为基于水平地表假设回折波走时方程的速度建模方法的反演结果图。
图4为本发明反演结果图。
图5为理论速度模型在Y=625m处X-Z速度切片图。
图6为基于水平地表假设回折波走时方程的速度建模方法的反演结果在Y=625m处X-Z速度切片图。
图7为本发明反演结果在Y=625m处X-Z速度切片图。
具体实施方式
为了使本发明目的、技术方案更加清楚明白,下面结合实施例,对本发明作进一步详细说明。
本发明所述的一种共面元共方位角道集的地震回折波走时反演方法的流程图,如图1所示,包括如下步骤:
步骤一:获取回折波走时数据;
步骤二:计算地表速度;
步骤三:获取共面元共方位角道集并编号;
步骤四:抽取当前共面元共方位角道集;
步骤五:将当前共面元共方位角道集内的地震道按照偏移距从小到大进行排序;
步骤六:逐道计算回折波射线参数和速度梯度;利用二分法求解如下方程计算回折波射线参数:
其中,x为偏移距,t为回折波走时,i为每一层编号,n为当前层编号,vs,i-1为第i层炮点端速度,vr,i-1为第i层检波点端速度,pn为回折波射线参数。特别地,当i=1时,vs,i-1即vs,0表示炮点处地表速度,vr,i-1即vr,0表示检波点处地表速度;
速度梯度计算公式如下:
其中,x为偏移距,i为每一层编号,n为当前层编号,vs,i-1为第i层炮点端速度,vr,i-1为第i层检波点端速度,vs,n-1为当前层炮点端速度,vr,n-1为当前层检波点端速度,pn为回折波射线参数,bi为第i层速度梯度,bn为当前层速度梯度。特别地,当i=1时,vs,i-1即vs,0表示炮点处地表速度,vr,i-1即vr,0表示检波点处地表速度;当n=1时,vs,n-1即vs,0表示炮点处地表速度,vr,n-1即vr,0表示检波点处地表速度;
步骤七:逐道计算回折点速度和深度。回折点速度计算公式如下:
其中,n为当前层编号,pn为回折波射线参数,vn为当前层回折点处速度;
回折点深度计算公式如下:
其中,vn为当前层回折点速度,vn-1为当前层回折点所对应的顶界面速度,bn为当前层速度梯度,zn-1为当前层回折点所对应的顶界面的深度,zn为当前层回折点所对应的深度。特别地,当n=1时,vn-1即v0表示面元中点处对应的地表速度,zn-1即z0表示面元中点处对应的深度;
步骤八:判断是否有未经处理的共面元共方位角道集,若有,则将该未经处理的共面元共方位角道集设置为当前共面元共方位角道集,返回步骤六;若无,执行下一步骤;
步骤九:利用所有已计算的回折点速度和深度,通过反距离加权插值算法计算离散速度模型所有网格点速度;
步骤十:输出速度模型。
进一步地,在步骤一所述回折波走时数据中,应包括炮点坐标、炮点井深、检波点坐标、回折波初至走时。
实施例1:
下面结合具体实施方式,对于本发明进行解释和说明。
为了进一步说明本方法的实现思路及实现过程并证明方法的有效性,用地表起伏且速度横向变化的3D理论模型进行测试,并和基于水平地表假设回折波走时方程的速度建模方法的反演结果进行比较:
S1:本发明的理论速度模型大小为3000m×1250m×600m,理论模型地表起伏且速度横向变化,最大高程差为244m;
S2:在地表布设了390个炮点和2023个检波点;炮间距75m,道间距25m,炮线距125m,接收点线距75m;
S3:采用空间2阶时间10阶精度的规则网格有限差分方法求解二阶常密度声波方程得到模拟的观测地震记录;
S4:从模拟的观测地震记录中拾取回折波走时数据,回折波走时数据包括炮点坐标、炮点井深、检波点坐标、回折波初至走时;
S5:用近道直达波计算检波点速度,通过反距离加权插值算法利用已知检波点速度计算地表速度;
S6:获取共面元共方位角道集并编号,面元大小为50m×75m,方位角为45°,将每个面元按照方位角划分为8个方位;
S7:抽取当前共面元共方位角道集;
S8:将当前共面元共方位角道集内的地震道按照偏移距从小到大进行排序;
S9:利用公式(1)和公式(2)分别逐道计算回折波射线参数和速度梯度;
S10:利用公式(3)和公式(4)分别逐道计算回折点速度和深度;
S11:判断是否有未经处理的共面元共方位角道集,若有,则将该未经处理的共面元共方位角道集设置为当前共面元共方位角道集,返回S8;若无,执行下一步骤;
S12:利用上述所有已计算的回折点速度和深度,通过反距离加权插值算法计算离散速度模型所有点速度;
S13:输出速度模型。
为更好说明本方法的有效性、适用性,截取照明度好、覆盖次数高的部分模型进行解释和说明,所截取的部分理论模型如图2所示,其大小为1000m×1000m×300m,X方向截取1000m-2000m,Y方向截取125m-1125m,Z方向截取0-300m。图3为基于水平地表假设回折波走时方程的速度建模方法的反演结果,图4为本发明反演结果。对比图3和图4,可以看出,基于水平地表假设回折波走时方程的速度建模方法的反演结果则存在明显的速度横向变化,且在纵向上存在速度突变;而本发明的反演结果速度横向变化平缓,纵向随深度线性递增,与理论模型比较相符。为更进一步说明本发明的有益效果,分别抽取X方向速度切片进行比较分析,图5为理论模型在Y=625m处X-Z速度切片,图6为基于水平地表假设回折波走时方程的速度建模方法的反演结果在Y=625m处X-Z速度切片,图7为本发明反演结果在Y=625m处X-Z速度切片,可以看出,基于水平地表假设回折波走时方程的速度建模方法的反演结果在横向上存在明显的条带状分布,且在纵向上一定深度处存在速度突变;而本发明的反演结果速度横向变化平缓,纵向无速度突变,更接近图5所示的理论速度切片。综上所述,本发明反演结果较基于水平地表假设回折波走时方程的速度建模方法的反演结果更接近理论模型,从而表明本发明在近地表起伏较大、速度横向变化的地区,可以取得比基于水平地表假设回折波走时方程的速度建模方法更好的实际应用效果。
以上所述仅为本发明的较佳实例而己,并不以本发明为限制,凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等,均应包含在本发明的专利涵盖范围内。
Claims (5)
1.一种共面元共方位角道集的地震回折波走时反演方法,其特征在于,包括如下步骤:
步骤一:获取回折波走时数据;
步骤二:计算地表速度;
步骤三:获取共面元共方位角道集并编号;
步骤四:抽取当前共面元共方位角道集;
步骤五:将当前共面元共方位角道集内的地震道按照偏移距从小到大进行排序;
步骤六:逐道计算回折波射线参数和速度梯度;
步骤七:逐道计算回折点速度和深度;
步骤八:判断是否有未经处理的共面元共方位角道集,如果有,则将该未经处理的共面元共方位角道集设置为当前共面元共方位角道集,返回步骤五;如果没有,执行下一步骤;
步骤九:利用所有已计算的回折点速度和深度,通过反距离加权插值算法计算离散速度模型所有网格点速度;
步骤十:输出速度模型。
2.根据权利要求1所述的一种共面元共方位角道集的地震回折波走时反演方法,其特征在于,所述步骤二中,地表速度由直达波走时数据计算得到。
3.根据权利要求1所述的一种共面元共方位角道集的地震回折波走时反演方法,其特征在于,所述步骤三中,共面元共方位角道集是指道集内任意两个炮检对连线夹角小于45°且这两个炮检对的中点之间的距离小于检波点线间距的地震道集。
4.根据权利要求1所述的一种共面元共方位角道集的地震回折波走时反演方法,其特征在于,所述步骤六中,回折波射线参数的计算公式如下:
其中,x为偏移距,t为回折波走时,i为每一层编号,n为当前层编号,vs,i-1为第i层炮点端速度,vr,i-1为第i层检波点端速度,pn为回折波射线参数;该公式中只有回折波射线参数pn未知,通过二分法求解该公式,从而得到回折波射线参数pn的值;特别地,当i=1时,vs,i-1即vs,0表示炮点处地表速度,vr,i-1即vr,0表示检波点处地表速度;
速度梯度的计算公式如下:
其中,x为偏移距,i为每一层编号,n为当前层编号,vs,i-1为第i层炮点端速度,vr,i-1为第i层检波点端速度,vs,n-1为当前层炮点端速度,vr,n-1为当前层检波点端速度,pn为回折波射线参数,bi为第i层速度梯度,bn为当前层速度梯度;特别地,当n=1时,vs,n-1即vs,0表示炮点处地表速度,vr,n-1即vr,0表示检波点处地表速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910520239.4A CN112099082B (zh) | 2019-06-17 | 2019-06-17 | 一种共面元共方位角道集的地震回折波走时反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910520239.4A CN112099082B (zh) | 2019-06-17 | 2019-06-17 | 一种共面元共方位角道集的地震回折波走时反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112099082A true CN112099082A (zh) | 2020-12-18 |
CN112099082B CN112099082B (zh) | 2021-08-20 |
Family
ID=73748296
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910520239.4A Expired - Fee Related CN112099082B (zh) | 2019-06-17 | 2019-06-17 | 一种共面元共方位角道集的地震回折波走时反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112099082B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113376689A (zh) * | 2021-03-30 | 2021-09-10 | 中国铁路设计集团有限公司 | 一种考虑层间多次波的弹性反射波走时反演方法 |
CN114488288A (zh) * | 2022-01-25 | 2022-05-13 | 中国海洋大学 | 共偏移距域起伏地表回折波走时速度建模方法 |
CN114879249A (zh) * | 2022-04-13 | 2022-08-09 | 中国海洋大学 | 基于四面体单元走时扰动插值的地震波前走时计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5235555A (en) * | 1992-10-23 | 1993-08-10 | Western Atlas International, Inc. | Diving-wave migration using airy functions |
CN104181592A (zh) * | 2014-08-28 | 2014-12-03 | 中国石油天然气集团公司 | 一种无射线追踪的回折波低阶多项式层析成像方法及装置 |
CN107870361A (zh) * | 2017-10-12 | 2018-04-03 | 南方科技大学 | 一种地震回折波层析成像方法、装置及终端设备 |
CN109490961A (zh) * | 2019-01-09 | 2019-03-19 | 中国海洋大学 | 起伏地表无射线追踪回折波层析成像方法 |
-
2019
- 2019-06-17 CN CN201910520239.4A patent/CN112099082B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5235555A (en) * | 1992-10-23 | 1993-08-10 | Western Atlas International, Inc. | Diving-wave migration using airy functions |
CN104181592A (zh) * | 2014-08-28 | 2014-12-03 | 中国石油天然气集团公司 | 一种无射线追踪的回折波低阶多项式层析成像方法及装置 |
CN107870361A (zh) * | 2017-10-12 | 2018-04-03 | 南方科技大学 | 一种地震回折波层析成像方法、装置及终端设备 |
CN109490961A (zh) * | 2019-01-09 | 2019-03-19 | 中国海洋大学 | 起伏地表无射线追踪回折波层析成像方法 |
Non-Patent Citations (1)
Title |
---|
李凌云等: "浅海海底表层速度建模技术及其应用研究", 《物探与化探》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113376689A (zh) * | 2021-03-30 | 2021-09-10 | 中国铁路设计集团有限公司 | 一种考虑层间多次波的弹性反射波走时反演方法 |
CN113376689B (zh) * | 2021-03-30 | 2022-04-12 | 中国铁路设计集团有限公司 | 一种考虑层间多次波的弹性反射波走时反演方法 |
CN114488288A (zh) * | 2022-01-25 | 2022-05-13 | 中国海洋大学 | 共偏移距域起伏地表回折波走时速度建模方法 |
CN114879249A (zh) * | 2022-04-13 | 2022-08-09 | 中国海洋大学 | 基于四面体单元走时扰动插值的地震波前走时计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112099082B (zh) | 2021-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112099082B (zh) | 一种共面元共方位角道集的地震回折波走时反演方法 | |
CN108196305B (zh) | 一种山地静校正方法 | |
CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
WO2008067588A1 (en) | Method of building a sub surface velocity model | |
CN107219554A (zh) | 陆地地震资料的剩余静校正量的自动获取方法 | |
CN112327362B (zh) | 速度域的海底多次波预测与追踪衰减方法 | |
CN110146924A (zh) | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 | |
CN111856577A (zh) | 一种降低逆时偏移地表炮检距道集计算量的方法 | |
CN110780341B (zh) | 一种各向异性地震成像方法 | |
CN110389377A (zh) | 基于波形互相关系数相乘的微震偏移成像定位方法 | |
US6980482B2 (en) | Seismic survey system | |
CN109490961B (zh) | 起伏地表无射线追踪回折波层析成像方法 | |
CN113960532A (zh) | 一种基于假想源的二次定位计算的微地震定位方法 | |
CN109188527B (zh) | 滩浅海地区快速建立三维近海底速度模型的方法 | |
CN107976714B (zh) | 一种复杂地表分级空间距离加权的超道集计算方法 | |
CN113777654B (zh) | 一种基于伴随状态法初至波走时层析的海水速度建模方法 | |
CN112630825B (zh) | 共炮检距域Beam叠前时间偏移成像方法、***、介质及应用 | |
CN104142513B (zh) | 一种提高地震勘探中层析反演精度的方法 | |
CN104570091B (zh) | 一种获得初至波射线的方法和装置 | |
CN106569278A (zh) | 一种多道相似相干速度谱计算方法 | |
US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
CN112394393B (zh) | 一种crp道集数据体重构的方法 | |
CN108802822A (zh) | 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 | |
CN116774293B (zh) | 一种同相轴自动拾取方法、***、电子设备及介质 | |
CN112255678B (zh) | 一种提高天然气水合物矿体分辨率方法及处理终端 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20210820 |