基于近似真地表浮动面叠前深度偏移成像技术应用研究
Application of PSDM imaging technology based on floating datum of approximate true surface
责任编辑: 叶佩
收稿日期: 2020-12-29 修回日期: 2021-11-2
Received: 2020-12-29 Revised: 2021-11-2
作者简介 About authors
窦强峰(1986-),男,工程师,主要从事地震资料处理方法研究工作。Email:
常规基于水平地表假设的CMP浮动面处理方法不适用于准噶尔前陆冲断带复杂构造区成像。选择近似真地表浮动面并统一作为时间域和深度域处理面,从叠前时间域预处理包括优选近地表底界面,地表一致性近似真地表静校正,地表一致性双平方根动校正求取速度及剩余静校正,使时间域预处理到叠前深度偏移保持流程、参数的一致性。该基准面的选择使波场走时更接近实际传播路径,成像精度更高,地质构造准确。
关键词:
The conventional CMP floating datum processing method based on horizontal surface assumption is not suitable for the imaging of the complex tectonic area in the Southern Junggar Thrust Belt.This paper selects the floating datum of approximate true surface as the processing datum in both the time domain and the depth domain.The pre-stack time-domain preprocessing includes selecting optimal the near-surface bottom boundary,the static correction for surface consistency of approximate true surface,surface-consistent DSR NMO velocity analysis,and residual static correction,in order to maintain the consistency of the process and parameters from the time-domain preprocessing to pre-stack depth migration (PSDM).The selection of the floating datum makes the wavefield travel time closer to the actual travel path and realizes higher imaging precision and accurate geological structure.
Keywords:
本文引用格式
窦强峰, 罗勇, 杨晓海, 谭佳.
DOU Qiang-Feng, LUO Yong, YANG Xiao-Hai, TAN Jia.
0 引言
准噶尔前陆冲断带位于准噶尔盆地天山南北麓,是近地表纵横向结构复杂和地下构造复杂的“双复杂”地区。地表主要为山地、冲沟,河流等,地形变化剧烈,表层介质复杂;地下多为山前推覆带发育背斜和成排成带的逆冲断块。从20世纪80年代开始,历经数十年的攻关研究,新疆油田在该区发现了齐古油田等重要油气,2019年高探1井获得重大突破,日产油气千余吨,打破国内陆上碎屑岩单井日产量最高记录,证明该区地区油气资源丰富,勘探潜力巨大。这主要得益于基于近似真地表浮动面叠前深度偏移成像等技术的不断发展及应用。
常规基于CMP浮动基准面叠前偏移处理技术,在山前复杂构造地区不再适用,主要表现在几个方面:①数据处理面问题,时间域处理CMP浮动面是静校正量平滑所得的一个时间量,并没有具体物理意义[1],而常规叠前深度偏移面主要是CMP浮动面的时间量乘以基准面静校正计算时的替换速度,认为CMP浮动面时间量均由替换速度引起,忽略低降速带对CMP浮动面的贡献,无疑就加入了算法引起的误差[2];②静校正问题,常规两步法静校正不具有地表一致性的特点,同一个CMP道集中的炮检点使用同一个基准面,但是不同CMP道集之间的浮动面是不同的,同时受地形及低降速带速度变化的影响,同一CMP道集中的高频静校正量差异大,很难满足最小静校正误差的准则[3-4];③动校正问题,传统单平方根算子动校正方法不具有地表一致性特点,随着CMP的位置变化,相同炮检点动校时移量也随之变化[5],同时采用对称的旅行时方程而存在较大的误差,导致地下构造成像的畸变,破坏同相轴的连续性。前人为此做了大量工作,主要偏向于偏移算法研究[6⇓⇓⇓⇓⇓⇓-13],1978年Stolt最先提出f-k域波动方程叠前偏移方法[14],1979年Yilmaz提出双平方根算子的扩展相移叠前偏移方法[15],1996年Popovic,刘文革等给出中点—半偏移距域的基于分布傅里叶的DSR叠前深度偏移公式[16-17],2002年Li等提出基于起伏地表的山地地震数据偏移方法[18],2009年Miao提出基于起伏地表的转换波双平方根动校公式和速度分析方法[19],2016年秦宁等发展了一种起伏地表直接成像的高斯束叠前深度偏移技术[20]。以往研究中,叠前预处理阶段依然采用基于CMP浮动面的常规处理方法,未考虑双平方根动校正等处理方法。基于以上分析,本文研究提出一种基于近似真地表浮动面的叠前深度偏移成像处理方法,可以有效提高山前复杂构造成像精度。
1 近似真地表浮动面静校正
图1所示为本研究的主要技术思路,包括统一时间域和深度域处理的浮动面、地表一致性近似真地表浮动面静校正、地表一致性双平方根动校正及剩余静校正、近似真地表叠前深度偏移等。
图1
这里提到的近似真地表静校正是指将地震资料从实际地表高程校正到地表高程小平滑面上,用来消除风化层厚度和速度对地震资料的影响,相对常规静校正,减小对实际波场走时的改造。
1) 近似真地表浮动面选择:准噶尔冲断带地区一般选择海拔1 500 m为固定基准面,替换速度在 2 000~2 500 m/s,利用静校正量平滑计算得到CMP浮动面,该面是按一定空间范围将炮、检点静校正进行非线性平滑的一个时间平滑量。通过替换速度将该面换算到深度域高程面(图2白线所示),可以看出该面是高于实际的地表高程面的,常规做法会以该CMP浮动面作为叠前深度偏移面,不难发现与实际炮、检点位置有较大误差。本次研究采用实际CMP高程的小平滑面作为数据处理浮动面(图2绿线所示),即近似真地表浮动面,CMP高程平滑半径据前人研究的经验,一般选择1/4~1/2排列长度,具体参数按不同地区实际试验结果确定。可以看出该面接近真实地表,和地表高程变化趋势相同,高程在实际地表高程面附近。从叠前时间域预处理到叠前深度偏移,炮、检点都校正到该面上,保证时间域和深度域浮动面的统一。
图2
图2
不同浮动面及底界深度域示意
Fig.2
The schematic diagram of different floating datum and surface bottom in depth domain
3) 地表一致性静校正量计算:高波数的道间时差不能通过常规成像及速度分析技术恢复速度的高频扰动,必须在成像处理前靠静校正及剩余静校正消除道间时差的高频扰动部分[21]。其原因是速度的高频扰动对旅行时没有影响,基于旅行时的速度估计无法反演小尺度速度变化。本文采用回转波层析成像方法建立近地表模型,该方法是建立一个速度随深度逐渐递增的速度模型,在这样一个模型下产生的初至波为回转波,利用其正演走时和实际初至走时之差,结合回转波路径进行速度反演。该方法可以稳定地反演出近地表速度异常体,且可以通过最大穿透深度确定近地表底界面,为后续深度偏移提供合理的浅表层速度模型,适合于复杂介质的成像。在模型迭代反演的过程中,实际初至的时间和反演模型所计算的初至时间,在理论上二者的数值差(此处定义为δt)等于0,就认为反演的模型和实际地表完全一致。但是在实际生产中,由于反演模型相对实际地表是一个低频趋势,无法真正刻画速度的高频变化,所以δt不可能等于0,我们将该值分解到炮检点分别得到两个短波长值(炮点为δts,检波点为δtr)。那么如图3所示,最终近似真地表浮动面的静校正值就应该是:
图3
式中:STs、STr分别为炮点、检波点静校正量;Vs、Vr分别为炮点、检波点处近地表速度;Vc为替换速度;Hs1、Hr1分别为炮点、检波点实际高程到近似真地表底界面之间的厚度;Hs2、Hr2分别为炮点、检波点近似真地表浮动面到近似真地表底界面之间的厚度。
图4
图4
不同静校正方法道集对比
a—静校正前CMP道集;b—近似真地表浮动面静校正后CMP道集;c—常规两步法静校正后CMP道集
Fig.4
The comparation among different statics CMP gathers
a—CMP gather before statics;b—CMP gather after rugged topography statics;c—CMP gather after conventional statics
2 地表一致性双平方根动校正
常规动校正是基于水平基准面或者CMP浮动面进行对称性的单平方根动校正,其中水平基准面动校正由于静校正量较大,对波长改造严重,动校不准确。CMP浮动面动校正在山前资料中,不同CMP道集相同炮检点动校正时移不同,是非地表一致性的,在均方根速度分析及剩余静校正中会引入较大误差。复杂构造地区,炮点走时和检波点走时很明显有较大差异,常规单平方根动校正无法从理论本身解决该问题而引入一定误差,这种误差会致使同相轴无法完全聚焦,产生构造畸变。
近似真地表浮动面的双平方根动校正,是在本文上述提到的地表一致性静校正基础上,在该浮动面上分别针对炮、检点开展非对称走时的一种动校正方法。图5中,ts、tr是炮、检点各自走时,a、b分别表示炮、检点位置相对于CMP点的高程差,t0表示CMP点处双程旅行时,x为炮检距,v为均方根速度。那么,CMP点处动校正方程可表述为:
图5
从上述方程中不难看出,炮检点相对于CMP点高程差异所造成的高频时差,是通过动校正的方式解决。图6a、6b分别是不同动校正方法的道集对比,其中图6a是在图4c所示数据基础上进行常规单平方根动校正后道集,图6b是在图4b所示数据基础上进行双平方根动校正后道集。不难发现,双平方根动校正相对常规单平方根动校正,道集成像质量明显改善。图7a所示为在近似真地表浮动面静校正基础上进行单平方根算子动校正的叠加剖面及道集,图7b所示为在近似真地表浮动面静校正基础上进行双平方根算子动校正的叠加剖面及道集,可以看出双平方根动校正叠加在构造主体部位成像合理(图7b),而单平方根动校叠加(图7a)明显存在随浮动面的扭动,且两翼成像质量差,说明近似真地表浮动面叠前预处理必须采用双平根动校正才能准确成像。同时还需在双平方根动校正道集求取相应剩余静校正,为叠前深度偏移提供偏前预处理数据。
图6
图6
不同动校正方法道集对比
a—常规静校正单平方根校正道集;b—近似真地表浮动面静校正平方根动校正道集
Fig.6
The comparation between different NMO CMP gathers
a—CMP gather after SSR NMO;b—CMP gather after DSR NMO
图7
图7
不同动校正方法叠加对比
a—近似真地表浮动面单平方根算子动校正剖面及道集;b—近似真地表浮动面双平方根算子动校叠加剖面及道集
Fig.7
The comparation between different NMO CMP gathers and stacks
a—CMP gather and stack after SSR NMO from rugged topography;b—CMP gather and stack after DSR NMO from ruged topography
3 近似真地表浮动面叠前深度偏移
为了节省篇幅,本文把具体偏移理论算法论述一并略去,仅展示本次研究方法的正确性和有效性。与常规叠前深度偏移不同之处是:
1)叠前时间域预处理和叠前深度偏移所采用的浮动面(偏移面)都是完全统一的,即近似真地表浮动面,更接近实际地表面;
2)近地表底界选取贴近地表高程,静校正量时差相对更小,对波场改造更小;
3)双平方根动校正可以求取更精确的均方根速度场,为叠前深度偏移提供初始速度,同时对应剩余静校正量更适应起伏地表。
图8为不同浮动面叠前深度偏移剖面和层速度叠合对比剖面,从图中箭头所指范围显而易见,近似真地表浮动面偏移同相轴聚焦,成像品质更好,浅层偏移归位合理,翼部高陡成像更加清楚。
图8
图8
不同浮动面叠前深度偏移剖面和层速度叠合对比
a—常规CMP浮动面叠前深度偏移剖面;b—近似真地表浮动面叠前深度偏移剖面
Fig.8
The comparation between PSDM stack velocity overlays on different topography
a—PSDM on RG topography;b—PSDM on rugged topography
4 结论及认识
1)基于近似真地表浮动面叠前深度偏移成像处理流程是一套完整技术序列,包括了地表一致性近似真地表静校正、地表一致性双平方根动校正及速度分析、浮动面剩余静校正、叠前深度偏移成像等,实现了时间域处理与深度域处理的一致性,提高了偏移成像精度,在准噶尔前陆冲断带山前复杂构造区取得良好效果。
2)常规基于水平地表假设CMP浮动面处理技术,在山前起伏地表成像应用条件及效果都受到一定条件限制,而基于近似真地表处理方法从叠前浮动面(偏移面)的选择、偏前时间域数据处理、近地表速度融合建模等关键技术环节,相对于常规处理都更接近实际波场传播规律,为进一步提高成像精度提供技术基础,更加适应于山前地表结构、地下构造复杂的“双复杂”地区的地震成像。
参考文献
关于浮动基准面与起伏地表面的讨论
[J]. ,
Discussion on floating datum and rugged topography
[J]. ,
准南复杂山地探区基准面静校正方法与质控策略
[J]. ,
Strategy of static correction in complicated mountainous area on the south margin of Junggar Basin
[J]. ,
最小静校正误差浮动基准面方法
[J]. ,
A new method for determination of floating datum in static corrections
[J]. ,
真地表动校叠加技术
[J]. ,
True surface dynamic correction and stacking techinque
[J]. ,
起伏地表叠前成像技术在川东高陡构造工区中的应用
[J]. ,
Application of pre-stack seismic imaging from rugged topography in complex structure survey of eastern Sichuan Basin
[J]. ,
双平方根方程三维叠前深度偏移
[J]. ,
Double square root equation 3D pre-stack depth migration
[J]. ,
窄方位地震数据双平方根方程偏移方法探讨
[J]. ,
Double square root migration methods of narrow azimuth seismic data
[J]. ,
双平方根算子模型的f-k方法
[J]. ,
Modeling of double square root operators in f-k domain
[J]. ,
SSR与DSR组合的波动方程速度建模方法及应用
[J]. ,
SSR and DSR wave equation velocity modeling technique and its application
[J]. ,
基于双平方根方程的保幅地震偏移
[J]. ,
Amplitude preserved seismic migration based on double square root
[J]. ,
山前带地震数据成像处理流程探讨
[J]. ,
Discussion on the imaging processing workflow for foothill seismic data
[J]. ,
复杂地表复杂地下地区地震成像技术研究
[J]. ,
A study on seismic imaging techniques in complex surface and subsurface areas
[J]. ,
复杂地表条件下波动方程叠前深度成像
[J]. ,
Pre-stack wave equation depth migration for irregular topography
[J]. ,
Prestack partial migration
[J]. ,DOI:10.1190/1.1441064 URL [本文引用: 1]
Prestack migration by split step DSR
[J]. ,DOI:10.1190/1.1444065 URL [本文引用: 1]
非零偏移距DSR叠前深度偏移
[J]. ,
Prestack depth migration by DSR in non-zero offset domain
[J]. ,
Surface consistent full wave imaging from topography for complex structures
[C]// ,
起伏地表高斯束叠前深度偏移技术及应用
[C]// ,
Gaussian beam pre-stack depth migration and application based on the rugged surface
[C]// ,
准噶尔盆地玛湖地区高密度地震勘探中静校正问题解决方法
[J]. ,
The solution to high dencity seismic statics correction in Mahu area of Junggar Basin
[J]. ,
面向起伏地表偏移成像的表层静校正方法
[J]. ,
A new static correction method for the migration from rugged topography
[J]. ,
山前带地震勘探策略与成像处理方法
[J]. ,
Seismic exploration strategy and image processing in mountain areas
[J]. ,
山前带地震勘探策略与成像处理方法
[J]. ,
Seismic exploration strategy and image processing in mountain areas
[J]. ,
天山南北复杂构造成像技术进展及应用效果
[J]. ,
Progeres and application of imaging technology for complex structure in the southern and northern Tianshan Mountains
[J]. ,
时空移动成像条件及偏移速度分析
[J]. ,
Space-time-shift imaging condition and migration velocity analysis
[J]. ,
速度误差和地震噪声对最小二乘逆时偏移的影响分析
[J]. ,
An analysis of the influence of velocity error and seismic noise on least squares reverse-time migration
[J]. ,
/
〈 | 〉 |