还剩3页未读,继续阅读
文本内容:
面向复杂地质的采集观测系统设计近年来,“两宽一高”技术在地震勘探中得到了广泛应用宽宽方向、高密度等技术相结合的方法逐渐取代了传统的采集方法,显著提高了地震资料的质量,降低了采集成本采集观测系统的成本和对地下介质的探测能力受很多参数影响,例如炮距、道距、最大偏移距、纵横比、覆盖次数、面元尺寸和覆盖密度等复杂构造地区,特别是逆掩推覆区,以及含盐丘构造地区的观测系统设计与优化波动方程照明分析方面,谢小碧叠前偏移成像方面,周旭本文在全波形反演方法研究的基础上,利用全波形反演的高分辨率特性,提出了将全波形反演应用于复杂高陡构造区采集观测系统设计的观点试验表明,只有达到全波形反演能够收敛的最小的“最大偏移距”才能保证全波形反演的准确性全波形反演成像技术1自Tarantola1984年提出反演方法以来在地质模型存在的情况下,地震波正演数据可以通过式⑴得到式中——从模型得到的数据,也可以认为是野外地震采集到的数据,这个数据的大小和质量由设定的采集观测系统或者说采集参数所决定;——地质模型,也就是我们常说的地下构造,其中包含我们需要勘探的地质目标一般认为,采集到的地震数据资料越丰富,通过地震资料处理后得到的内容也就越完全,但米集成本也就越高Tarantola认为,可以通过地震数据恢复出地震模型,但我们知道,如果地震数据中包含的信息极少,就算再好的反演方法也不能得到满意的结果目前,几乎所有全波形反演研究人员都把精力集中在全波形反演方法的本身,而没有考虑到数据采集系统对反演结果的影响在已有地震数据的情况下,利用式⑵可以求出地质模型模型速度恢复的准确度直接决定了叠前深度偏移成像的结果对于采集观测系统设计,不考虑模型的体积模量和密度信息,使用简化的标量声波方程⑶进行全波形反演,可以节省大量计算时间3式中将方程加入边界条件和初值条件,离散后会得到相应的差分方程利用有限差分法或者有限元法能够对方程进行求解,并得到不同位置、不同时刻的空间位置的波场用准确模型计算得出的声波压力场定义初始模型波场记录与观测到地震之间的差值为对于得到的地震记录可以通过最小二乘方法进行求解定义全波形反演的目标函数为式中——检波点位置;——总的记录时间;——可视为残差函数的方程⑸的最优化解受炮点位置、检波点位置、记录时间等参数的约束要使目标函数达到最小,就要用初始模型模拟得到的反演速度逐次迭代过程中的一个非常重要的步骤是速度模型的更新,这个步骤的核心是梯度场的计算,为了提高计算效率,这里避免Ferchet矩阵导数的计算,而采用伴随方法计算梯度这种方法利用炮点正传波场与检波点残差的逆时传播波场进行零延迟互相关,逐步累加形成残差函数的梯度伴随方法计算梯度的计算公式如下式中不同时刻炮点正向传播波场;伴随方法计算梯度的过程中使用了炮点正传波场与检波点残差的逆时传播波场的互相关,这个步骤把地震波入射信息与反射信息结合起来,形成了反射全波形反演同样的技术也在叠前深度偏移成像的过程中进行类似的应用,但叠前深度偏移主要利用波形的相位信息,达到同相叠加的目的,最终形成以反射同相轴为反射信息的地震偏移成像结果,其成像分辨率主要受偏移速度模型、地震子波和最大偏移距的影响如图1所示,对于全波形反演方法,由于最终反演得到的是速度模型,而不是反射波同相轴的地震剖面,所以它不受地震波长的限制,其分辨率要比地震偏移成像的分辨率高得多,在此基础上,要获得良好的速度模型分辨率,就要求有比叠前偏移更长的最大偏移距这也是地震全波形反演比常规叠前偏移成像对采集观测系统要求更高的原因之一为此,这里发展了面向全波形反演的地震采集观测系统设计优化研究最大偏移距的影响2由于三维波动正演与全波形反演的计算量巨大,为了阐明道理,我们这里只使用二维模型进行数值试验图2为国内某工区模型,模型深度为3500%模型长度为15000%模型速度范围从1800m/s到3500m/s模型中包含5个层位,并被多个断层所截断,由于断距较大,因此可以认为是o横向速度变化较大的地质结构后面的数值计算结果都基于此模型由于全波形反演对最大偏移距信息比较敏感,因此首先分析不同的最大偏移距对反演模型的影响为了提高计算速度,将模型离散化成25mX25叫在模型地表共布设81炮,炮点范围从2500m至12500叫炮点间距为125叫检波点位于炮点两侧接收,检波点间距为25叫设计最大偏移距分别为500m、1000m、1500m、2000m和2500m五种观测系统方案,由于模型横向范围只有15000m,考虑到满覆盖区域的问题,因此没有设置更大最大偏移距的滚动接收方案使用以上5种观测方案进行正演模拟得到合成记录,我们粗略地认为它们就是野外观测到的地震记录在此我们使用标量声波方程,空间10阶差分精度,时间2阶差分精度,2ms采样间隔,6s记录长度由于模型网格尺寸的限制,为了避免数值频散,震源子波采用9Hz的雷克子波,同时为保证记录范围内波形的完整性,子波的延迟时采用222mso图3中分别为模型8000m处不同偏移距的有限差分正演模拟记录,从地震记录中可以看此随着最大偏移距的增加,单炮记录内的信息逐渐增多,特别是复杂构造断面波和绕射波都可以在单炮记录中出现在500m最大偏移距的情况下,深层的断面波没有出现,一方面是由于接收排列太短没有接收到,另一方面是由于正演模拟算法使用的模型范围过小全波形反演过程中,初始模型的选择会影响反演的收敛速度,很多学者已经进行了相关的研究,本文不做重点分析我们对原始模型使用1000m的平滑半径进行平滑,得到的模型作为初始模型(图4)进行反演可以看到,相对于原始模型(图2),初始模型中高波数成分已经丢失,只含有低波数成分,对于这样的初始模型,只要观测系统中含有足够的信息就可以反演出较为精细的速度模型图5是使用不同最大偏移距观测方案与初始模型(图4)进行全波形反演迭代39次得到的结果反演过程中约束了最小速度1800ni/s与最大速度3500m/s能够看此当最大偏移距为500mo时,经过全波形反演,不但没有收敛到精细的模型,反而破坏了原有低波数模型的准确度,尤其是在模型的中深部产生了类似于波场传播的噪音随着最大偏移距的增加,反演对模型的破坏越来越小,在最大偏移距为2000m时,深部构造逐渐显现,看似噪音,但与正确模型相比,其构造是正确的当偏移距达到2500m时,其结果已经相当令人满意,也没有破坏初始模型,整个过程一直处于收敛状态,这说明此模型在此种条件下,2500m应该是保证全波形反演能够收敛的最小的“最大偏移距”,低于此值,将不能保证全波形反演的收敛这里的数值例子也从侧面证明了很多实际资料很难完成全波形反演任务的原因之一就是观测系统不能满足全波形反演对最大偏移距的要求由于模型横向长度的限制,同时考虑到覆盖次数的均匀,我们只算到了2500m最大偏移距的滚动排列方式为了说明最大偏移距的重要性,我们采用全模型范围激发和接收,其中炮点间隔为125m,检波点间隔为25m,共计121炮,其他正演模拟参数与反演参数同上图6为全模型范围激发与接收39次迭代的全波形反演结果其反演结果质量要远远好于2500m最大偏移距方案,图中断层界面清晰可见,与准确模型之间的差别已经很小,可以用作叠前时间偏移的初始模型为了说明本文提出方法的有效性,我们与目前工业界使用的方法进行对比在面向复杂模型的观测系统设计中,现在最流行的方法就是波动照明分析及逆时偏移成像的观测系统分析对于上述使用的5种观测系统,我们分别从波动照明分析和逆时偏移的分析思路进行计算图7中展示的是不同方案的照明分析结果从结果可以看出随着最大偏移距的增加,不同观测方案对应的照明能量也逐渐增加,对应的观测系统能够探测的深度也在增加,但我们无法判断出哪种观测系统能够满足我们的勘探要求图8是图2所示模型中最下层(L6)对应的不同观测方案下的照明曲线,通过该曲线也能看出,随着最大偏移距的增加,照明能量也在逐渐增强,但同样无法说明哪个方案的观测系统能够满足要求图9中展示了不同最大偏移距方案的逆时偏移结果,其中偏移速度模型使用准确速度模型,随着最大偏移距的逐步增加,逆时偏移结果逐渐变好,由于500m的最大偏移距过小,在偏移剖面中存在大量的噪音,同时断点位置不是很清晰,并且由于逆掩推覆构造的阻挡作用,弱照明区域的成像也不是很好但所有构造都归位到正确的位置当最大偏移距增加到1500m时,逆时偏移成像结果已经很好,断点和断面都比较清晰,下伏岩层成像准确当继续增加最大偏移距时,成像结果继续变好,但变好程度有限这说明,当偏移速度模型准确时,叠前深度偏移成像对数据的要求并不高,利用叠前偏移成像来分析观测系统不是一个有效的手段而实际资料很难获得准确成像的一个主要原因就是很难获得准确的速度模型,所以速度模型的准确建立要比常规叠前偏移成像的要求更高,而全波形反演恰恰是一种很有潜力的速度模型建立手段若想获得良好的地震资料,需要这种面向全波形反演的地震采集观测系统设计与优化的手段基于全波形反演的观测系统设计3地震勘探采集设计中,需要对采集参数进行论证,在对最大偏移距的论证过程中,常常采用分偏移距叠加的思路判断最大偏移距是否能够满足要求,但在对地质目标精度逐渐提高的今天,水平叠加的思路已经远远不能满足要求,应运而生的面向复杂模型的波动照明分析在参数论证上发挥了巨大的作用,但照明分析也存在其自身特有的弱点,它只体现了能流密度,没有充分考虑到反射波动相位信息,所以不能真正地体现偏移成像,这使得波动照明分析很难在精细勘探中发挥作用而叠前深度偏移成像分析没有考虑到速度模型的建立,在准确速度模型的条件下,会使最大偏移距、道距、炮密度等参数的敏感性较差本文提出将全波形反演应用于采集观测系统设计,由于全波形反演的高分辨率特性,使得全波形反演对采集观测系统的要求有了进一步的提高,尤其在复杂高陡构造区,基于全波形反演的观测系统设计将会给采集观测系统设计打开一扇大门,为复杂区观测系统设计提供有力的支撑由于篇幅限制,本文只论证了最大偏移距的影响,其完全可以用于论证分析道距、炮距、炮道密度等与观测系统相关的其他各项参数;在三维情况下,还可以分析观测系统的纵横比、观测方位等参数。