应用Vista55处理地震数据流程
一、 二层水平介质模型
模型基本参数:单边放炮,每炮24道接收,共12炮,道间距25m,炮间距50m(2个道间距),偏移距250m(10个道间距)。
采样率2000微秒,每道采样点1000个。
反射界面深度800m,上层介质速度2500m/s,下层介质速度3000m/s。
11 数据的输入
首先选择File/New Project新建一个Project,如下图:
按住 不放,出现按钮组合 ,选择 新建一个二维数据集,然后按住 不放,在出现的按钮组合 中选择 ,在d出的OpenDialog中选择该数据文件Msgy,单击OK即可将该数据加入到数据集中,如图所示:
界面上部显示了该数据集的基本信息。可以单击 查看该数据集,如下图:
界面上部显示了该地震图,下部是各道的炮间距,状态栏显示的是鼠标处的道号和振幅信息。
该模型是由射线追踪模拟出的理想二层水平介质模型,不需要做什么预处理,可以直接进行下面的实质性处理。在做实质性处理之前,必须给数据建立观测系统,并将观测系统相关信息写入道头,以便进行实质性处理。
12 建立观测系统
在Data List窗口的数据集M中点击 ,或者选择Interactive/Geometry Window Display,在d出的对话框中选择M数据集,即可出现观测系统界面,默认出现的是设置炮检关系及炮点坐标界面,在第一行中填入相应得增量,如下图:
主要参数增量为炮点增量2个站点(桩号),首尾检波器桩号也相应增2,炮点坐标增量为2个桩的长度50m
按住 不放,在出现的组合按钮 中选择增加单行按钮 增加一行,并设置相应值作为初始值,如下图:
设第一个炮点位于第1个站点,坐标为0m,因此第一炮的第一个检波器位于第11个站点,最后一个检波器站点位于第34个站点。
然后按住 不放,在出现的组合按钮 中选择增加多行按钮 ,在d出的增加炮点对话框中填入剩下的炮数,如下所示:
单击OK即可填充,填充完后如下图:
设置完炮检关系及炮点坐标后,点击 设置检波器坐标,其基本方法与上面是一致的,
检波器站点增量为1,坐标增量为25m,初始设置为:第1个检波器,即第1炮的第1个位于11号桩,坐标为250m,然后填充剩下的检波器个数(即填充到最后一炮的最后一个检波器,位于56号桩)即可。填充后如下图:
然后点击 并点击 计算CMP及检波器叠加次数和炮检距离,以便写入道头相应的位置,完成后续的处理。可以点击 查看CMP及检波器叠加次数,并可以此判断建立的观测系统是否正确。如下图:
可以看出,CMP最大覆盖次数为6次,根据已知的模型参数信息,由公式可以算出理论上最大的覆盖次数
N=241/(22)=6
因此上图所示说明建立的观测系统是正确的。可以按住保存按钮 不放,在出现的按钮组合 中选择 将相关的信息写入道头,并还可以选择 将炮检关系及他们的坐标写入文本文件,以便以后查看。
加入观测系统后(主要是将相关的信息写入道头),便可以对数据的输入按所要的处理模块选择不同的排列方式,这对于后面的实质性处理是非常重要的。
13 速度分析
可以选取若干个CMP道集进行速度分析(可以选择速度谱法和常速度叠加法(CVS)),以便获得最佳的叠加速度,为随后的动校正提供速度。
速度分析首先要得到速度谱,常速度叠加图和一个最佳炮检距道集,其工作流如下所示:
首先要选择做速度分析的道集,输入数据选择M,#Bins From Center Bin填5,即抽取以选择的CMP为中心,左右各5道共11个CMP道集。右键→Data Input Control→Data/header Selection,出现
表示目前没有定义可选的数据集,单击 建立可选的CMP道集,如下图:
点击 增加一行,定义要选择的CMP道集号,如下图所示:
即选择了CMP道集号为30的CMP道集,因此通过VelZone的道集为CMP号为25,26,27,28,29,30,31,32,33,34,35共11个CMP道集。
设置常速度叠加法的速度扫描范围为100,5000,100,速度谱法的速度扫描范围为100,5000,100,如下图所示:
设置完毕,点击 执行上面的工作流,执行完后得到三个输出的数据集,分别是速度谱,一个最佳炮检距道集和常速度叠加图。
选择Interactive/Velocity Tools/Interactive Velocity Analysis进行速度分析,在d出的对话框中选择速度谱,一个最佳炮检距道集和常速度叠加数据集,如下图:
单击OK后即可进行速度分析,拾取最佳叠加速度,见下图:
可以单击 拾取速度,拾取速度可以根据三个图谱拾取最佳叠加速度,见下图:
拾取到速度后,点击保存按钮 保存为速度文件,为动校正准备。
14 动校正
动校正是将CMP道集中不同炮检距的各道校正为共中心点的自激自收道。其工作流命令如下:
数据输入选择数据集M,数据排列方式为CMP_NO,如下图:
双击动校正命令图标,设置动校正速度文件(即速度分析中得到的速度文件),如下图:
设置完毕点击 执行工作流,即可得到动校正后的结果。见下图:
15 水平叠加
水平叠加将动校正后同一个CMP道集的各道叠加为一道。其工作流命令如下:
双击Input命令图标,数据输入选择动校正后的数据,数据排列方式选择CMP_NO,如下图:
双击CmpStr命令图标,选择叠加方式,如下图:
OK确定,单击 执行工作流即可得到水平叠加剖面图。水平叠加剖面图如下:
局部放大,见下图:
可以看出反射界面大致位于双程时640ms处,可以由模型参数算出理论反射界面的双程时时间:t=2800/2500=064s=640ms,可以看出得到的水平叠加剖面基本是准确的。
二、 Marmousi模型
模型基本信息:单边放炮,96道接收,共240炮,道间距25m,炮间距25m,偏移距200m(8个道间距)。注:该模型炮点位于检波器排列之后。
采样率4000微秒,每道采样750个点。
基本的处理方法与上面处理模型的方法是差不多的。
21 数据的输入
数据的输入与上面处理模型时是一样的,就不再重复了,下图是加入Marmousi模型数据:
不知该模型数据已作过滤波处理还是模拟时的原因,数据反射波能量很强,深部的反射也很强,反射信息很丰富,信噪比较高,不需作其他的处理就可以做动校正和水平叠加了。之前仍需给数据加观测系统。
22 建立观测系统
建立观测系统的方法与前面也是一样的,也不再重复了,仅就个别参数作些说明。下图是炮检关系和炮点坐标:
炮点是位于检波器排列之后的,即
第1个检波器......第96个检波器←炮点
因此和上面二层水平介质模型的炮检关系有一点小的差别,只要正确认请炮检的位置关系和准确确定炮检坐标,其实是一样的。
下图是检波器坐标:
同样点击 和 计算CMP和检波器的覆盖次数,点击 显示,见下图:
CMP最大覆盖次数:N=961/(21)=48次,因此上面所建的观测系统是正确的。
按住保存按钮 不放,在出现的按钮组合中 中选择 将相关的信息写入道头,并还可以选择 将炮检关系及他们的坐标写入文本文件,以便以后查看。
23 速度分析
选取几个CMP道集进行速度分析,可以选取CMP200,300,400,500,600这6个CMP道集进行分析,获得这几个CMP道集的速度谱,常速度叠加图和一个最佳炮检距道集,以便拾取这几个CMP的速度。
工作流如下:
双击VelZone图标,填入参数,如下图:
右键点击VelZone图标→Data Input Control→Data/header Selection并选择CMP道集,如下图所示:
并设置常速叠加和速度谱扫描范围,设为100,5000,100,设置完毕后,点击 执行工作流。执行完毕即可得到这几个CMP道集的速度谱,常速度叠加图和一个最佳炮检距道集。选择Interactive/Velocity Tools/Interactive Velocity Analysis进行速度分析,在d出的对话框中选择速度谱,一个最佳炮检距道集和常速度叠加数据集,如下图:
点击OK后即可对各CMP道集进行速度分析,拾取各CMP速度,见下图:
可以通过左右移动 来选择不同的CMP道集,拾取不同的CMP的速度。拾取完后,点击保存按钮保存为速度文件,为后面的动校正提供速度信息。因模型中部结构比较复杂,反射信息比较复杂,很难准确拾取到反射层的速度,造成速度不准确或不充分,对后面的处理有一定的影响。
24 动校正
得到速度信息后,可以对按CMP排列的地震数据进行动校正。工作流如下:
动校正后浅层且大炮检距的地震道会发生动校正拉伸现象,导致浅层的同相轴遭到严重损害,为处理方便,根据动校正拉伸后频率下降,对动校正后的数据设置一个带通滤波器,滤掉因拉伸而导致的畸变。
双击Input图标,设置输入数据及其排列方式,如下所示:
双击Nmo命令图标,设置动校正速度文件选择为速度分析中得到的速度文件,可以点击按钮 查看速度信息,见下图:
点击 查看各CMP的速度信息,见下图:
带通滤波器设置如下:
设置完毕后点击 执行工作流,执行完毕后得到动校正后的结果。见下图:
25 水平叠加
水平叠加将动校正后同一个CMP的各道叠加为一道,以加强有效信号,消除干扰,提高信噪比。水平叠加工作流见下图:
双击各命令图标并设置完参数后,执行工作流即可得到水平叠加剖面,见下图:
利用已给的Marmousi模型速度(已知的是z-层速度,应将其转化为t-Rms)文件,对Marmousi模型做动校正并水平叠加后也得到一个水平叠加剖面,但是由于CMP号对应不上,效果也没有改善多少,也有可能是方法本身的问题。水平叠加剖面见下图:
对数据做反褶积压缩地震子波后再做动校正和水平叠加,叠加剖面的局部有所改善,见下图:
因为Marmousi模型结构比较复杂,反射信息很复杂,仅仅通过动校正和水平叠加得到各CMP自激自收的信号,这是很难得到地下反射界面的信息的,成像效果也不可能很好,这是模型本身的结构和方法本身所决定的。
(一)合成地震记录的原理及制作过程
合成地震记录即是将测井资料经人工合成转换为地震记录(地震道)。据第五章的讨论,我们知道地震记录x(t)是地震子波w(t)与地层反射系数序列r(t)的褶积,即
x(t)=w(t)r(t) (7-1-1)
因此,理论合成地震记录制作的处理流程为
反射波地震勘探原理和资料解释
其实现的过程可以示意地表示于图7-1-1,图中(b)可是由测井资料得到的地下地层速度模型。由地层速度模型通过下列公式算出反射系数序列(如果测井资料中有密度(ρi)资料,则下式中的vi→Zi=vi·ρi):
反射波地震勘探原理和资料解释
图7-1-1 地震记录的形成及其反演的原理
反射系数序列与子波褶积便得图中(g)的合成地震记录(地震道)。其中子波的选择则是通过设计若干个地震子波,与测井反射系数曲线褶积得到若干个地震道,把这若干个地震道与井旁实际地震道进行对比,选用符合程度最好的作为要选用的子波。这种从一口井的资料计算出一道合成地震记录的过程也叫一维地震模型。在实际地震解释中用得更多的是二维地震模型,它的实施过程如下:
(1)从实际地震剖面出发,提出一个地层岩性或油气分布的初步解释方案。
(2)把这个解释方案作为地质模型,将模型的参数输入计算机,计算出合成地震剖面。
(3)比较合成地震剖面和实际地震剖面,看两者是否符合。如果两者符合,说明解释方案正确;否则,就要修改解释方案,并重新进行(2)、(3)项工作,直到解释结果和实际情况能满意地符合为止。
(二)合成地震记录的应用
(1)与实际地震记录对比,能更好地解释地震波的层位,用于地震层位的标定。
(2)利用包括和不包括多次波的合成记录与实际记录对比,能更好地分辨出多次波。
(3)利用多道合成地震记录制作合成地震剖面图来确定和寻找地层和岩性油气藏,以及详细的含油气砂层变化等。
合成地震记录是用声波测井或垂直地震剖面资料经过人工合成转换成的地震记录(地震道)它是地震模型技术中应用非常广泛的一种,也是层位标定、油藏描述等工作的基础,是把地质模型转化为地震信息的中间媒介 是联合高分辨率的测井信息与区域性的地震信息的桥梁,其精度直接影响到地质层位的准确标定目前油气勘探工作越来越向隐蔽性油气藏发展,目标尺度越来越小,对合成地震记录提出了更高的要求[1]由于合成记录的制作过程中存在诸多的制约因素,合成地震记录与实际地震剖面往往并不能完全一致本文将着重从原理上分析这些制约因素,并总结相应的技术对策 1、合成地震记录的制作原理 合成记录的制作是一个简化的一维正演的过程,合成记录F(t)是地震子波S(t)与反射系数R(t)褶积的结果 F(t) =S(t)×R(t) 合成地震记录制作的一般流程是:由速度和密度测井曲线计算得到反射系数,将反射系数与提取的地震子波进行褶积得到初始合成地震记录根据较精确的速度场对初始合成地震记录进行校正,再与井旁地震道匹配调整,得到最终合成地震记录 2、合成记录的制作过程中的质量控制 合成地震记录制作主要包括计算反射系数、提取地震子波和匹配调整时深关系三个环节,如何对这三个环节进行有效的质量控制成为决定最终精度的关键,在实际制作过程中依照下列方法进行相应的质量控制[2]
由于崎岖海底的存在,横跨海底界面强烈的侧向速度变化使得下伏地层随着海底起伏,构造形态严重畸变,根本不能反映构造的真实面貌,严重影响了下伏地层的地震成像。近年来,针对崎岖海底采用多种手段,从实验室正演模拟、采集参数的试验、处理方法及时深转换等做了大量的一系列的攻关,进一步揭示了深水崎岖海底区地震波传播的本质特征,及崎岖海底对地震波的影响机理和成像畸变的因素;通过对崎岖海底区地震处理的攻关,特别是对崎岖海底区绕射多次波的压制改善了地震资料的品质;通过对层替换技术、波场延拓技术、叠前深度偏移处理等多种方法进行了处理试验,确定了叠前深度偏移对崎岖海底的处理流程,解决了由崎岖海底造成的构造畸变问题。
同时,在长电缆大偏移距条件下,有些常规处理技术已不能应用,如以双曲线反射走时为基础的动校正,速度分析和水平叠加以及压制多次波的方法。近来国际上速度分析的研究可归纳为三个方面:一是叠前速度分析方法向非双曲线反射走时方程为基础发展,二是改善层速度的计算方法,三是偏移速度分析方法发展迅速,这与叠前深度偏移的兴起有关,主要是层析成像方法。
(一)已有地震资料分析
深水区多年度陆续采集了部分地震资料,有些资料由于年代久远已无法利用,为此已有地震资料分析主要有目的地针对1979年和1997年采集的地震资料进行分析。主要针对噪音分析、主要干扰波类型、多次波发育分布情况等方面进行分析。
1噪音分析
噪音分析主要是评估涌浪噪音的分布频带和固有噪音的主要频带,我们采用的分析手段主要是FK分析和频谱分析。涌浪噪音主要是低频噪音,其频带主要集中在10Hz以下。固有噪音的频带主要集中在30~65Hz之间,其主要噪音源是地震采集船的螺旋桨转动。
2主要干扰波类型、多次波发育分布分析
主要干扰波类型、多次波发育分布分析主要是评估干扰波类型、多次波发育的主要频带。采用的分析手段主要是FK分析和频谱分析。干扰波的主要类型是线性干扰。线性干扰波的主要频带分布集中在20Hz以下。
多次波主要表现为海底多次等长周期多次波,其频带分布与一次波极其相似,主要能量集中在30~60Hz之间,能量较一次有效反射强,掩盖了有效的一次波反射,并等时重复出现。其次,多次波还表现为崎岖海底区的绕射多次波。由于存在崎岖不平的海底,海底多次波在地震剖面上的反映也不一致,海底较平时,由于多次波和正常地层速度上的差异,可以通过Tau-P域去多次等传统的方法来消除,但海底崎岖造成海底的角度很大的斜层,这种很强的海底斜层产生的多次波,由于其速度和下覆地层没有太大的差别,就很难通过常规的方法来消除,使得地震剖面的中深层横向能量很不均匀,造成偏移剖面划弧现象(图5-1)。
这些多次波不但严重干扰了凹陷内有效反射,造成凹陷内地震资料信噪比极低,而且对基底反射也产生较强的干扰,严重影响了该地区地震资料的地质解释和研究。因此,压制和消除多次波成为深水地震资料采集和处理的重点。
通过分析,复杂海底与地下结构是影响该地区资料品质的主要因素。深水地震资料具有以下特点:海底构造复杂,水深变化剧烈,侧面反射以及斜坡带内能量反射很弱;噪音以低频干扰、中深层高频干扰、异常值为主;浅层的频带较中深层宽,中深层信噪比和分辨率低;多次波干扰以深层海底及长周期多次波为主,能量强,存在散射多次波;崎岖基底引起的中深层速度拾取复杂。
图5-1 崎岖海底区强绕射多次波
(二)处理技术方法
根据对原始资料的分析,对工区地质情况的调查,结合地质任务和处理要求,采取的处理对策为:SRME、高精度拉冬和LIFT技术组合多次波衰减技术;通过确定性子波处理和沿海底构造处理的串联组合多道反褶积技术压制延续相位;针对信噪比很低的斜坡带,采用频谱整形技术提高该区域资料的信噪比;针对凹陷内随机噪音严重的地方,采用多域去噪技术提高信噪比;进行高精度速度分析,构造复杂部位加密控制点,对目标区前后剖面进行认真对比,反复迭代以提高速度分析的正确性及合理性;利用叠前深度偏移解决该区崎岖海底及高陡构造成像问题。
1多次波衰减技术
衰减多次波是本次地震数据处理的重点和难点之一。虽然压制多次波的方法有很多,但没有一个能在所有条件下除去所有的多次反射波。
针对工区多次波具有的特征,经过多次试验,采取了SRME(海底多次波衰减)、高精度拉冬及LIFT多域组合多次波衰减技术,通过三步法对多次波进行逐步压制,并取得了非常理想的效果。
在海上地震勘探中,近道多次波是最难处理的相干噪音之一,特别是在浅层气的影响下,近道多次波更是难以压制。常规衰减近道多次波的技术是预测反褶积组合内切除,该技术简单有效,但在衰减多次波的同时,有效信号也被去掉了,破坏了道集的完整性,给后续处理带来一定的困扰。
本次我们研发了一种有效衰减近道多次波的LIFT技术,该技术是根据AVO原理模拟有效信号,通过局部时窗进行信噪分离。实践证明,该技术既能有效衰减近道多次波,又能很好地保留有效信号,为后续处理打下坚实的基础。
2串联组合反褶积技术
海洋地震资料因为采集时气q沉放离海面有一定的深度,所以在气q因压力爆炸后的瞬间会产生较大的气泡升至海面,再加上涌浪的干扰,期间的信号因此来回摆动,所以经检波器接收到的信号中就产生了延续性的相位。此次采集所产生的延续相位在浅水和深水区的深层均有较强的体现,有的甚至覆盖有效信号。所以针对严重的延续相位,采用了确定性子波反褶积和多道反褶积串联组合来压制,并取得了较为理想的效果(图5-2)。相比统计性子波反褶积,确定性子波反褶积更有针对性,且有效地保护了浅层信号、频率振幅特征。
图5-2 串联反褶积效果图
3频谱整形技术
针对斜坡带与基底低信噪比区域,在迭代前采用了频谱整形技术来提高信噪比(图5-3)。
图5-3 频谱整形效果对比图
4多域去噪技术
斜坡、凹陷等处能量反射很弱,造成信噪比很低,采取多域去噪技术提高信噪比。多域去噪方法是利用信号和噪声在不同域的差异,将干扰波与有效波的差异最大化,分别在炮域、共偏移距域,运用拟三维FXY滤波、线性干扰消除等技术提高地震资料的信噪比(图5-4)。
图5-4 多域地震资料信噪比对比图
5高精度速度分析技术
在常规数据处理方法中,速度分析普遍采用相干性度量法。这种方法没有考虑相近或干涉同相轴、剩余静校、非双曲型时差等有关的噪声以及其他非随机噪声的影响,因此影响了时间和速度分辨率。本次处理采用的是最新研制的相位相关统计方法,这种方法的优点是在时间和速度分辨率上比常规方法更高更可靠,更有助于对小幅度构造的分析和解释。
时间分辨率的检测:合成CDP道集中,两组同相轴中间隔均为30ms,从图5-5中可以看出,相位相关统计速度谱与常规速度谱相比,时间分辨率明显提高。
速度分辨率的检测:使用同一时间而采用不同速度的两个同相轴,速度差由大到小不断变化,观察速度谱中能量团,直至分不开为止。从图5-6中可以看到,当常规速度谱中能量团已分不清楚时,使用本方法,同一深度的两个能量团可清楚地分开,特别在深部,效果更为明显。
图5-5 两种速度谱分辨率对比
图5-6 两种速度谱能量团对比
6叠前深度偏移(PSDM)成像技术
崎岖海底最核心的问题就是:由于存在崎岖不平的海底,横跨海底界面的侧向速度强烈变化,使得地震射线路径复杂,时距曲线为非双曲线,常规处理方法中的CMP道集不再是共反射点道集,叠加剖面不再是零偏移距剖面,造成下伏地层的成像差及构造形态的严重畸变。陈礼、葛勇等人利用理论模型讨论了用常规时间偏移、叠后深度偏移及叠前深度偏移技术解决深水崎岖海底地震成像问题的有效性。通过对深水模型数据各种偏移结果的对比分析认为,常规时间偏移和叠后深度偏移均不能解决崎岖海底地区地震成像问题,而叠前深度偏移是解决这一问题的有效方法。
叠前深度偏移技术通常用来实现复杂构造准确偏移成像,解决复杂地质问题。对于地下深度成像而言,最棘手的问题不是偏移方法,而是地下速度模型的建立。深度偏移是一个迭代过程,是一个不断建立模型、试验模型、运行偏移,根据成像修正模型的反复过程。叠前深度偏移对地下形态基本不作假设,速度深度模型直接用叠前资料建立,地下速度纵横向均可变化,CMP道集考虑非双曲效应。由此得到的数据体不但能提高信噪比、使空间归位正确,同时还能直接得到地质上合理的深度成像数据体,供地质解释之用,很显然是解决崎岖海底造成构造畸变的一个较好的办法。图5-7为过LW3-1构造的叠前时间偏移和叠前深度偏移的剖面,对比可以看出,时间偏移剖面LW3-1构造区周边构造倾角大,下伏地层结构成像精度较低,结构不清楚,构造形态严重畸变,而叠前深度偏移在纵横向分辨率、振幅的相对保持、对复杂构造的能量聚焦以及构造形态方面都获得比较明显的改善,能够满足地质解释的要求。
图5-7 三维叠前时间偏移剖面和叠前深度偏移剖面对比图
经过多次试验和论证,我们选择崎岖海底发育并可能有火山岩发育的白云6-1构造区资料进行叠前深度偏移试验。
图5-8(上)为04EC2458深度偏移成果剖面,图5-8(下)为该测线最终偏移时间剖面,从深度和时间剖面的对比来看,深度剖面保持了原有的分辨率和信噪比,剖面面貌比较自然,大部分地区海底崎岖影响基本消除,随海底崎岖起伏的同相轴基本上被拉平,反映了地下真实的构造形态,但局部地区(图5-8左部)仍然还存在上下地震反射起伏共鸣的现象,说明海底崎岖的影响仍未消除,分析这些局部海底崎岖影响仍未消除的情况,可以发现,这些不理想的情况的上方海底为一些较浅的海沟,仔细观察,可以发现这些较浅的海沟里充填了较厚的沉积物(图5-9),通过速度分析,发现这些沉积物层速度很低,大约1670m/s,比1500m/s稍高,但比起隆起上的地层速度1820m/s要低得多,如此低速的沉积物可能是一些晚期沉积的淤泥。
图5-8 04EC2458叠前深度偏移剖面和最终时间偏移剖面对比
图5-9 海底崎岖速度分析
通过计算,如果海底有一个400ms的海沟,可能造成下伏地震相位的下拉达75ms,而如果海沟之中有充填物250ms,则可造成下伏地震相位下拉达25ms,由此可见,下伏地层构造形态的畸变不仅是海底崎岖本身的影响,海沟中充填物的厚度也有相当的影响,而在这次的叠前深度偏移中没有考虑这个因素,所以在一些局部地方尚需要进一步改善。
7叠后LIFT提高信噪比处理技术
由于本工区构造复杂,浅、中、深层信噪比和频率成分差异很大,我们采用LIFT去噪处理技术,有效提高处理成果的质量(图5-10)。
图5-10 LIFT技术信噪比对比图
1叠前 AVO属性处理
(1)AVO速度调整和层速度场的建立
除了振幅处理不当会造成 AVO 分析陷阱外,如果速度分析处理不当也可能造成技术陷阱,同时也会影响到资料的品质和横向分辨率。虽然常规处理中已经产生了一个比较准确的速度场,但这个速度场由于考虑到多方面因素的影响,并不能完全适应AVO属性处理。
因此,在AVO 速度分析中,首先对信噪比低的剖面段继续采用细致的常速扫描与谱点加密的方法,保证速度拾取的精确程度 (图5-17),从而获得精度较高的叠加速度场。由于本区地层为倾斜地层,因此叠加速度场与均方根速度场之间存在以下关系:
VRMS=VSTK×cosθ
图5-17 精细速度调整
式中:VRMS为均方根 速度;VSTK为 叠 加 速度;θ为地层倾角。
根据实际资料,本工区地层倾角在10°~15°之间,cosθ约在097~099之间,因此我们可以利用以上的公式将叠加速度场近似转换为均方根速度场。
对叠加速度场进行平滑 (图5-18),然后将平滑后的叠加速度场转化为层速度,结合地质层位解释对层速度场进行调整 (图5-19),最后利用该速度场完成初次AVO 处理并对处理结果进行分析,根据实际处理效果情况再次调整速度场,完成最终 AVO速度场的建立 (图5-20)。
图5-18 叠加速度场
图5-19 均方根速度场转化的层速度场
图5-20 排2井三维 AVO速度场建立流程
(2)AVO入射角度的试算与选择
入射角=tan-1(offset/2 H)
式中:offset代表偏移距;H 代表目的层深度。
从共 中 心 点 面 元 道 集 内 的 高 差 变 化 情 况 看, 过 排 2 井 道 集 (x:313980,y:4982130)高程在 289~2922m 之间,最 大 高 差 为 32m;工 区 南 部 农 田 区 道 集 (x:309000,y:4974410)高程在2978~3025m之间,最大高差为47m;工区南部水库区道集 (x:309000,y:4963400)高程在3194~3064m之间,最大高差为13m。在完成近地表校正后,共中心点面元道集内的高差对 AVO 入射角计算的影响很小,可忽略不计。
排2井三维地层埋深差异较大,同一地层最大埋深差达1500m,因此取500~2750m为目的层深度范围。根据现有资料,排2井油层深度为1014m,经分析,该层对应地震资料最大偏移距为1070m,根据 AVO 入射角计算公式,满足 AVO 分析的最大入射角为43°;排2井石炭系深度为1374m,经分析该层对应地震资料最大偏移距为1400m,根据AVO入射角计算公式,满足 AVO分析的最大入射角为45°;三维工区内石炭系最大地层埋深为2750m,按最大偏移距2870m 计算,满足 AVO分析的最大入射角为46°。综合分析,试处理中确定的最大入射角为45°。
为了保证 AVO 属性分析的效果,对 AVO 属性分析的角度进行试算,分别试算了0°~30°、0°~35°、0°~40°、0°~45°(图5-21~图5-24)。最终 确 定 0°~35°作 为 本 区 AVO属性分析试处理的入射角。
图5-21 入射角0°~45°碳氢检测剖面
图5-22 入射角0°~40°碳氢检测剖面
图5-23 入射角0°~35°碳氢检测剖面
图5-24 入射角0°~30°碳氢检测剖面
(3)AVO属性体处理
应用分选出的0°~35°角度限制道集,结合全区层速度场,利用Shuey近似公式进行AVO属性体的叠加处理,得到 AVO属性体。
另外也对另一种近似方程——Richards方程,抽取过排2井的纵线进行了试处理,并取得了一定的效果 (图5-25,图5-26)。
图5-25 过排2井(排8井)纵线λ剖面
图5-26 过排2井(排8井)纵线μ剖面
(4)P波数据体的后续处理
为了满足后续反演处理和综合分析的要求,还要对P波数据体进行后续的处理工作,主要是应用 Omega处理系统的STOLT偏移方法对 P波数据体进行叠后偏移处理;应用零相位反褶积、蓝色滤波提高P波数据体的分辨率;应用三维 RNA提高P波数据体的信噪比 (图5-27,图5-28)。
图5-27 P波偏移剖面
图5-28 提频去噪后P波偏移剖面
对于P波属性数据体和成果数据体,尤其是过井线,进行了详细的对比分析,认为P波数据体与成果数据体基本相当,从合成记录对比 (图5-29,图5-30)上看,两者的频率、相位、能量都一致,同样都可以很好地反映地质现象,但在细节方面,P波剖面振幅的强弱关系反映更明显,保幅性更好 (图5-31~图5-46)。
图5-29 过排2井成果剖面标定图
图5-30 过排2井P波剖面标定图
图5-31 过排2井P波剖面
图5-32 过排2井纵线成果剖面
图5-33 过排8井P波剖面
图5-34 过排8井纵线成果剖面
(5)AVO属性体的归位处理
由于 AVO属性处理是直接应用叠前道集进行属性体的叠加,因此 AVO 属性体的归位一直是一个难题。
图5-35 过排201井纵线成果剖面
图5-36 过排201井P波剖面
图5-37 过排208井纵线成果剖面
图5-38 过排208井P波剖面
图5-39 过排9井纵线成果剖面
图5-40 过排9井P波剖面
图5-41 过排12井纵线成果剖面
图5-42 过排12井P波剖面
图5-43 过排16井纵线成果剖面
图5-44 过排16井P波剖面
图5-45 过排17井纵线成果剖面
图5-46 过排17井P波剖面
在本次处理中,把P波数据体和梯度 G 数据体分别进行偏移,再将偏移后的 P、G数据体进行相应运算,获得偏移归位后的碳氢检测、拟波松比等属性体,彻底解决了这个难题,并取得了很好的效果 (图5-47,图5-48)。
图5-47 偏移前过排2井-排8井纵线碳氢检测剖面
图5-48 偏移归位后过排2井-排8井纵线碳氢检测剖面
(6)AVO属性处理效果分析
AVO叠前属性处理取得了较好的处理效果。在过排2井附近的道集上 (图5-49)可以看到明显的 AVO正异常现象。
从过排2井 (排8井)P波剖面、碳氢检测剖面、拟泊松比剖面上可以看到,油井与负相位砂体吻合很好 (图5-50~图2-52)。
图5-49 过排2井道集 AVO正异常显示
图5-50 过排2井-排8井纵线P波剖面
图5-51 过排2井-排8井纵线碳氢检测剖面
图5-52 过排2井-排8井纵线拟泊松比剖面
而过排201井 (排204井)、排203井、排208井P波剖面、碳氢检测剖面上,没有明显的反映 (图5-53~图5-58)。
图5-53 过排201井-排204井纵线P波剖面
图5-54 过排201井-排204井纵线碳氢检测剖面
图5-55 过排203井纵线P波剖面
图5-56 过排203井纵线碳氢检测剖面
图5-57 过排208井纵线P波剖面
图5-58 过排208井纵线碳氢检测剖面
排9井、排12井、排16井和排17井在碳氢检测剖面等属性剖面上均没有明显 AVO反映 (图5-59~图5-66),这与实钻结果也是吻合的。
图5-59 过排9井纵线P波剖面
图5-60 过排9井纵线碳氢检测剖面
图5-61 过排12井纵线P波剖面
图5-62 过排12井纵线碳氢检测剖面
图5-63 过排16井纵线P波剖面
图5-64 过排16井纵线碳氢检测剖面
图5-65 过排17井纵线P波剖面
图5-66 过排17井纵线碳氢检测剖面
2叠后属性处理
当储层物性和充填在储层中的流体性质发生变化时,会造成地震反射系数、传播速度、振幅、频率等多种属性的变化。这些变化表现为波形、能量、频率、相位等一系列基于运动学、动力学的地震属性的变化。地震属性比地震剖面在检测储层或流体性质变化方面敏感得多,并且许多地震属性都是非线性的,它将增加预测的准确性。鉴于本区目前勘探存在的困难,有必要开展叠后属性处理工作,提高勘探的成功率。
(1)精细标定及构造解释
精细构造解释是进行属性提取工作的基础,只有如此才能保证所提取的地震属性能够准确反映所研究目的层段或储层的特征。需要做好以下三个方面的工作:
极性判断:首先进行正演分析对比法。采取正极性子波和负极性子波分别进行排2井自激自收正演,可以看到正极性子波正演结果中储层附近波组特征表现为上弱波峰,下强波峰,中间夹一个强波谷的特点 (图5-67)。而表现负极性子波正演结果中储层附近波组为两个相对弱的波峰夹一个相对强波谷的特征 (图5-68)。对照过排2井的地震剖面,可以发现地震剖面上储层处的地震响应特征与正极性子波正演结果一致 (图5-69)。其次采用能量判识方法——选择排2井靠近塔西河组下部的一套较厚的含砾细砂岩,厚度13m,地震资料可以分辨其顶底。其顶底分别对应地震的波峰和波谷。从反射系数曲线上可以看到顶部反射系数大于底部反射系数 (图5-70),所以顶部反射在地震资料反射中对应能量应该大于底部反射能量。从地震资料读取该反射层附近的能量,可以看到波谷能量最大在-1300附近,而波峰能量最大达到5000左右 (图5-71),从对应关系上看,波峰顶应该对应含砾细砂岩的顶,这只有在使用正极性子波条件下才能达到该条件,使合成记录道和地震道相对应。再其次采用正负子波标定对比法——从排2井正负子波合成记录对比标定剖面上 (图5-72)可以看到正极性子波与负极性子波在目的层段标定效果都不错,但在11~12s处正极性子波合成记录波组与地震波组更加匹配。最后采用多口井综合标定法——采用多口井标定对比,发现正极性子波标定结果与实际地震道对应效果良好。综合以上四种方法,判定该区地震资料为正极性。
图5-67 排2井正极性子波正演结果
图5-68 排2井负极性子波正演结果
层位及储层标定:本次研究对车排子地区已钻探井均进行了合成地震记录标定,标定采用如下原则——以井点附近地震道提取子波、利用VSP做为时深关系指导、以塔西河组及沙湾组底部反射为标志层,在此基础上进行细微的调整。通过标定认识到塔西河组底界以及沙湾组底界为连续强振幅波谷反射同相轴,全区可追踪,对应地震反射层为TN1t、TN1s,沙一段1砂组底部为较强连续振幅反射同相轴,全区基本可追踪,对应地震反层为 TN1s1。为方便层位解释,采用了波峰反射的解释作为控制层位 (图5-73)。
图5-69 排2井地震剖面正极性子波标定结果
图5-70 排2井塔西河组下部含砾细砂岩顶底反射系数对比
图5-71 排2井塔西河组下部含砾细砂岩顶底反射能量对比
图5-72 排2井正负极性子波标定结果对比
图5-73 排2井区标准层、控制层位、储层标定结果
由于排2井油层只有39m 厚,因此地震资料对其分辨能力及其在地震剖面上的对应关系需要精细标定。从声波时差曲线上计算得到,排2井油层砂岩平均速度2120m/s,泥岩平均速度2450~2600m/s,在油层顶为正反射系数,油层底为负反射系数,因地震资料为正极性资料,故而波谷对应油层顶界,波峰对应油层底界。由于地震资料视主频70Hz,以1/4波长产生调谐波为最大分辨率,最大分辨厚度应为9m。对于39m 的砂层无法区分,但由于排2井油层发育在泥岩段中,理论情况下只有砂体顶界面会产生地震反射,并且三维地震资料有效频宽大至10~110Hz,从而进一步提高了垂向分辨率,使得排2井砂层在地震剖面上有响应。因此得到结论:①排2井区三维地震资料,在沙一段1砂组对应的地震反射中,强振幅波谷反映了砂岩存在,并对应砂岩顶界,强波谷的横向变化反映了砂层的横向变化;②砂层顶界对应强振幅波谷,下部较强振幅波峰与强波谷相连,波峰与强波谷之间反映了多个砂层存在,砂层总厚度较大。
断层及层位解释:研究区内构造解释采用断裂和控制层位同时进行解释的方案,并对整个三维工区地震资料进行了解释。主要利用了时间切片、相干体、三维可视化等多种技术,理顺了断裂结构,落实了构造 (图5-74);编制了车排子地区排2井三维区塔西河组底界、沙一段1砂组底界、沙湾组底界、白垩系底界等4层构造图 (图5-75)。
图5-74 排2三维断裂系统图
图5-75 排2三维四层构造图
(2)地震属性提取技术研究
地震属性的提取方式有剖面提取属性和层面提取属性两种。剖面提取属性可以获得研究目标的纵向信息以及点与点之间横向变化情况。沿层提取属性获得的是各类属性沿界面横向变化的信息,常用来预测薄储层和与断层有关的隐蔽油气藏。各种不同属性分类都有对应的地质意义 (表5-4),用来指导工作中采用合理的属性提取方法。
表5-4 地震属性分类表
另外还有其他常用属性:
方差体属性:利用地震数据中相邻道之间地震信号的相关性,通过计算样点的方差值,揭示数据体中的不连续信息。其作用在于进行断层 、岩性识别 (大时窗利于分析大断层,小时窗利于分析岩性体、小断层)。
地震波吸收衰减:该现象是由岩石基质的固有黏d性,包括颗粒之间和裂隙表面的内摩擦损耗、孔隙岩石内液体相对流动、局部饱和效应以及几何漫射等引起的。影响地震波吸收衰减的主要因素有岩石性质、岩石孔隙度和孔隙内流体成分等。当地震波在地下传播时,随着离震源的距离的增加,能量逐渐衰减。而一些特定的因素可能加速能量的吸收,如天然气的存在能引起高频段的地震波异常高的吸收率。在一定的时窗内,分析地震能量的吸收,作为频率的函数测量能量衰减的速度,可以检测储层的变化。对于裂缝性油藏,裂缝、溶孔以及含油气性都会引起储层的孔隙度、饱和度、层速度和地震振幅频率等属性的变化,从而引起地震吸收系数的变化。因此,可以利用地震能量吸收分析预测裂缝储层的发育情况。
工作中结合工区储层特征提取敏感地震属性共6种:以排2井钻遇油层为例 (图5-76),振幅类属性提取了均方根振幅、累加负振幅、平均波谷振幅;复地震道属性提取了瞬时频率、瞬时相位;频能统计类提取了弧线长度。另外提取了方差体、地震波吸收衰减等两种属性。
从地震资料与属性对比图 (图5-77)中可以看出,所提取的属性异常边界与地震资料同相轴波形、能量变化点相对应,可以说,所提取的属性是能够反映所研究的地质目标的。
图5-76 沿排2油层多种属性平面图
图5-77 沿排2油层多种属性异常边界与地震同相轴边界对比图
图5-78 沿排2油层平均波谷振幅属性分频段平面图
为更加深入研究频率、速度谱信息,还采用了分频段属性分析、时频分析、速度谱分析等技术手段。
分频段属性分析:从地震资料有效频段中按 10~30Hz、30~50Hz、50~70Hz、70~90Hz、90~110Hz共5种频段范围分别进行地 震 资 料 6 种 敏 感 属 性 提 取 试 验(图5-78)。通过对比分析,认识到该地震资料 50~70Hz是最佳的属性频段。由于50~70Hz地震资料的理论分辨厚度范围为8~11m,而工 区 中 钻 遇 油 层 厚 度 范 围 为2~5m,因此,所提取的各种属性中体现的异常并不反映砂层厚度概念,仅是地层物性、含流体性等变化的反映。
时频分析技术:为拓展和提高该区纵向(时间轴方向)频率属性研究的深度,针对有利地震资料频段50~70Hz范围,开展此项技术研究,期望能够发现油层段在时间轴方向存在有规律异常。工作中采用以下工作流程——首先从三维地震资料中抽取过井的二维测线,然后对其进行50~70Hz带通滤波,再对滤波结果求取瞬时频率属性,最后抽取过井点CDP处瞬时频率值与时间交汇得到成果图件,如图5-79。通过分析钻遇油层井油层位置频率特征及相邻井对应段的频率特征认识到时频变化没有规律,因此,该技术不能应用于该区储层研究。
速度谱分析:速度谱资料往往在油气分布处有异常反映,可以凭借该特征辅助判定油气的存在。为研究该区速度谱特征,在原有高精度速度谱基础上,针对目的层段,缩小速度扫描时窗,提高速度谱变化精细程度,期望能够发现有利规律。但速度谱资料在研究目的层段没有针对油层或可能储层的速度异常现象,因此该技术不适用于该区(图5-80)。
图5-79 排8及附近井储层时频分析图
图5-80 排2、排201井点处速度谱图
3多井约束反演处理
(1)合成记录标定与子波求取
反演过程中的合成记录不同于层位解释时的合成记录,它的标定要求细节更加精细,合成记录道中每个同相轴都有地震道同相轴相对应,这样才能将地质和地震精确对应起来。
合成记录标定的过程是反射系数与子波褶积的过程,子波求取的过程是合成记录与反射系数反褶积的过程,两者是正反运算的有机整体。合成记录与地震子波是影响反演处理过程中的时深关系、初始波阻抗模型与波阻抗反演结果是否准确的重要因素,而一个高质量的合成记录与地震子波的获得是一个循环反馈过程:利用实际地震资料多道记录自相关统计的方法,在一个经验时深关系 (排2井 VSP速度)的控制下,先利用一个主频70Hz初始标准雷克子波 (实际地震资料的主频为70Hz)作最初的合成记录道,将此合成记录道与井旁地震道对比,做测井曲线与实际地震资料之间的时深关系校正,在校正后合成记录上选目的层段的合适的时窗提取子波,并用此子波重作合成记录,校正时深关系,如此反复,直到合成记录与实际地震资料在能量、相位、频率等方面都匹配程度很高时,认为所得到的合成记录与所提取的地震子波是合适的。
为保证合成记录标定的可靠程度,在完成单井标定后,提取标定速度与本区 VSP速度进行比较,从对比图 (如图5-81)中可以看到所有井的标定速度与 VSP速度一致性良好,这说明标定是可靠的。单井标定完成后,为保证标定结果在横向上一致,还需要进行多井横向标定。图5-82中可以 看 到 标 定结 果 在 横 向 上是一致 的,特别是排2井、排8井钻遇的储层情况与实际地质情况一致。
图5-81 各井标定速度与本区 VSP速度对比图
图5-82 排2、排8井连井标定剖面
(2)地质模型的建立
建立尽可能接近实际地层沉积情况的波阻抗约束模型,是减少反演最终结果多解性的十分重要的环节。建立波阻抗模型的过程实际上就是把地震界面信息与测井波阻抗信息正确结合起来的过程。地震资料包含着区域的构造信息,控制模型的横向变化;测井资料包含丰富的高、低频信息,控制模型的纵向阻抗变化关系,为波阻抗界面间的地层赋予合适的波阻抗信息。声波、密度测井资料在纵向上详细揭示了岩层的波阻抗变化细节,三维地震资料则在三维空间内记录了波阻抗界面的地震反射。测井资料在三维地震地质反射界面内合理内插外推的结合,为精确地反演出地层波阻抗数据提供了有效的先验约束模型。
地下沉积体的空间接触关系是十分复杂的,计算机无法一次确定各个层位之间的拓扑关系,因此建立地质框架是通过地质框架结构表按沉积体的沉积顺序,从下往上逐层定义各层与其他层的接触关系。由于本区存在着上超这种现象,因此在模型的建立过程中,必须在地质框架结构表中定义出来。通过合理的定义上下层位的接触关系,使建立的初始波阻抗模型 (图5-83)符合实际的地下沉积模式,沙一段1砂组表现为从北向南逐渐加厚的特征,2、3砂组则表现为基本厚度一致的特征。
图5-83 排2三维反演初始波阻抗模型图
(3)反演参数选择
针对不同地区的资料特点选择适合该区的反演参数是反演项目的质量保证。结合对该区基础资料及地质特征的认识,对稀疏脉冲反演中的多个敏感参数进行试验和选择,特别针对λ、子波影响、频带补偿、色标范围调试等四个方面。根据稀疏脉冲反演的目标函数可知,地震反射系数的稀疏和合成记录与原始地震道的残差最小这两项是相互矛盾的,这是由于在算法上,它遵循以下原则:λ值小,强调反射系数之和最小,即强调稀疏性,稀疏脉冲反演剖面细节少,分辨率低,残差大;否则反之。但是λ值太大,过分强调地震残差最小,一味地使合成记录与原始地震道吻合,结果使一些噪音也加到了反演剖面中,同时由于忽略了反射系数的稀疏,使得反演结果失去了波阻抗纵向变化的低频背景。因此,在反演参数调试中很重要的一步就是寻找一个合适的λ值,使得反演剖面既保持细节又不损失低频背景,这个工作是通过对井旁边合成记录与原始地震道吻合程度的控制来完成的。λ值可以用Jason软件反演的质量控制工具来确定。据此选定本次反演应用的λ值为16。
该区反演面积比较大,实际地震子波受施工因素及实际地层物性特征的影响在能量、相位、频率等特征上会有微小差异,针对这种情况,采用空间特征变化的空变子波进行反演 (图5-84),这有助于对地下地层特征进行正确反演,使得到的反演结果更加接近地下岩层的真实地质特征。保证空变子波在有效频带范围内基本稳定而略有差异,满足了该区反演的实际需要。
图5-84 排2三维空变子波图
针对本区储层,尤其是油层厚度薄的情况,采取了合理的高频补偿,补偿示意图见图5-85,这使得储层分辨率得到合理的提高。
图5-85 排2三维反演高频补偿示意图
色标调试是正确反映储层的关键步骤之一。本次反演采用了剖面色标调试和三维立体色标调试两种方法。剖面色标调试采用将油层顶底投在剖面上,调整色标,逐渐使色标变化范围与厚度一致,同时注重储层横向变化 (图5-86),最终得到色标调试结果。三维立体色标调试采用三维立体镂空方法,将钻遇油层范围镂空出,调整色标范围,使油层范围与实际钻探范围一致,并记录色标变化点 (图5-87)。结合两种色标调试方法,再精细调整,最终得到合理色标范围,并将油层颜色调整为醒目的黄红色。
图5-86 排2三维反演结果剖面色标调试图
图5-87 排2反演结果三维立体色标调试图
(4)稀疏脉冲反演处理
上述合理的时深关系、准确的层位断层数据、校正过的测井数据、空变子波及高精度的三维约束模型等是下一步反演处理的数据基础。在反演处理时首先选择多条二维连井骨架剖面进行了大量、反复的试验,采用严格的质量控制,检查并适当调整反演参数,最大程度地保障反演结果的可靠性。
考虑到反演出的波阻抗数据体仍然相对缺乏高频、低频信息,我们对其做了高低频信息补偿。将前面生成的含有丰富高频和低频信息的初始模型数据体与所得到的带陷阱的反演波阻抗结果做匹配合并,补偿其缺乏的频率成分。
约束地震反演过程,是所用测井数据、钻井、试采数据、构造层位解释数据、地震数据等各种数据紧密结合反演,并根据地质储层变化情况不断加深认识、反复修正,逐步完善反演结果的过程。每反演出一次结果,处理、解释人员就结合在一起,对效果进行反复对比、分析,根据掌握的地质和各井钻探,钻采资料提出下一次反演处理应改进的问题和措施,如此反复循环处理。通过以上处理技术和质量控制手段,得到最终反演数据体。
(5)反演效果分析
稀疏脉冲反演是测井约束地震反演技术中最为可靠的技术,在目前的储层描述与评价中得到了广泛的应用。该技术成功地将地震资料与高频丰富的测井资料相结合,充分发挥了地震在平面上连续采集、测井在纵向上分辨率高的优势,使点与面达到和谐的统一,把用于构造解释的常规地震资料的界面型剖面转换成可与钻井资料直接对比的岩层型测井剖面,给储层的追踪、描述以及预测工作带来了方便。其反演结果与地震资料所具有的振幅、频率、相位等特征都有较好的对应关系。
纵观反演结果,其具有以下显著特点:常规地震剖面,其波峰、波谷的极值点对应地层的分界面,是界面型剖面;而测井反演处理的资料,其波峰、波谷对应的是岩层,是岩层型界面,实质是层速度剖面。
反演结果如何,可以通过以下两点分析:
1)井点处反演的结果与井的吻合程度。反演结果是否与实钻井吻合,可以通过参与反演的井和未参与反演的井加以验证。反演结果与排2井、排8井等的钻井结果吻合的很好 (图5-88)。
2)反演结果符合地质变化规律。可以从反演资料同一层系地层波阻抗的变化是否均匀,反演结果的沉积模式是否与地质规律吻合等进行验证。如排201-排204连井反演剖面(图5-89)上,排201井钻遇沙湾组I砂层,排204井没钻到该砂体,这与实际地质情况是吻合的。
图5-88 过排2-排8井反演剖面
图5-89 过排201-排2-排204井反演剖面
以上就是关于vista 5.5地震处理软件的使用手册,以及这个软件的详细介绍(软件优势、技术思路、技术特点)。谢谢了!全部的内容,包括:vista 5.5地震处理软件的使用手册,以及这个软件的详细介绍(软件优势、技术思路、技术特点)。谢谢了!、合成地震记录、合成地震记录的概念,等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)