一种可控震源自适应地表一致性反褶积方法技术

本发明专利技术涉及石油地球物理勘探数据在处理技术，是一种快速有效的可控震源自适应地表一致性反褶积方法。具体步骤是：使用常规的可控震源激发并记录地震波数据；采用可控震源地震记录并在选定的时窗进行吸收分析，获得炮点或检波点在ｔ↓［１］～ｔ↓［２］时窗内Ｑ吸收模型道的振幅谱；对炮的扫描信号做傅氏变换，得到扫描相关子波的振幅谱；求出反褶积因子ｈ（ｔ）；得到浅、中、深层波形明显改善的可控震源和炸药震源激发的成像剖面。本发明专利技术可实现提高分辨率和相位转换，具有地表一致性统计的自适应能力，对于含有较大近地表吸收影响的可控震源初至波得以明显改善，使地震资料的处理质量明显提高。

A method of adaptive surface coherence deconvolution for vibroseis

The present invention relates to the processing technique of petroleum geophysical prospecting data, and is a fast and effective vibroseis adaptive surface consistency deconvolution method. Specific steps are: using vibroseis conventional excitation and seismic wave data; using vibroseis seismic record and the window in the selected absorption analysis, get shot or receiver in t: 1 ~ t, the lower amplitude absorption window Q model 2 the spectrum of gun sweep; scanning signal Fourier transform, amplitude scanning correlation wavelet spectrum; calculate deconvolution factor H (T); imaging profiles are generated deep and shallow wave significantly improved vibroseis and explosive source. The invention can realize high resolution and phase transformation, adaptive ability with surface consistent statistics, for vibroseis containing larger near surface absorption effect to wave in order to improve the quality of seismic data processing, which significantly improved.

【技术实现步骤摘要】

本专利技术涉及石油地球物理勘探技术，具体是一种对地震资料处理的可 控震源自适应地表一致性反褶积的方法。
技术介绍
地震勘探是采用人工激发地震波勘测地下石油、天然气及煤田等有关 的地壳结构和地层岩性方法的总称。由于地下介质的力学性质不同，地震 波在传播中，会引起地震波的频率、振幅、相位、速度及波场类型等参数 的空间变化，通过在地表或井下采用高精度记录设备和相应观测方法可以 获得这些特征参数变化的数据，并经地震数据处理技术的去伪存真和成像， 最终转换为地质学家(包括油气田及煤田开发工程师)可以识别的图像。 地质学家利用这些地震信息不仅可以了解地层构造的几何形态和断裂大小 及分布，而且可以了解古地理和古沉积相信息，进而预测储层与油气藏(煤 田)的关系。地震勘探一般可以分为地震数据采集、地震数据处理和地质解释三个 大的环节。地震数据采集是通过人工激发地震波，利用地面接收仪器记录 由地下传播回来的地震波所引起的地面振动信号的过程。地震资料处理是 根据地震波的传播理论，利用计算机等设备和相应的处理软件，对野外采 集的地震数据进行加工处理，以获得能反映地下地质结构的地震信息。地 质解释是基于地质理论和地震勘探理论的结合，将地震信息转换为地质认 识，最终获得地质目标的认识，寻找有利含油气圈闭。自70年代初可控震源用于地震勘探以来，它以其灵活、经济、环保等优越性而得到了充分肯定。可控震源是一种利用机械振动激发地震波的人 工震源装置，它具有人为控制振动信号的扫描时间和扫描频率的特点，相 比炸药震源更加灵活、高效。然而，可控震源相关子波通常被认为是零相 位子波，多年来对于可控震源采集数据处理主要是采用纯相位滤波方法来进行零相位至小相位的子波转换，这种方法对于含近地表和大地Q吸收的 实际地震数据存在一定问题。实际可控震源相关子波在大地Q吸收衰减的 作用下不再是零相位，而是混合相位，并且随着Q吸收的增大，其影响也 随之增加。为此，本文研究提出了，该方法的主要特点是通过计算近地表和大地吸收衰减重构可控震源吸收子波，并通过子波反褶积达到消除近地表和大地吸收衰减引起混合 相位子波问题，从而达到提高可控震源数据的分辨率和相位转换的目的， 同时该方法具有相对保持储层振幅、频率、相位和波形信息的能力。
技术实现思路
本专利技术针对目前可控震源地震采集数据在处理过程中存在的问题，提 出了一种快速有效的可控震源自适应地表一致性反褶积方法。本专利技术提供以下可控震源自适应地表一致性反褶积的方法，具体步骤是1) 使用常规的可控震源激发并记录地震波数据；2) 采用可控震源地震记录并在选定的时窗进行吸收分析，获得炮点或 检波点在t t2时窗内Q吸收模型道的振幅谱；所述的时窗是为了减少规则干扰波影响所述的分析是对时窗内地震数据做傅氏变换，采用频率域振幅中值滤波 消除噪音。上述的频率域振幅中值滤波即将某一炮点或检波点各个地震道的傅氏变 换结果按采样点排列后取中值，获得一个炮点或检波点模型道的振幅谱。3)对炮的扫描信号做傅氏变换，得到扫描相关子波的振幅谱，然后用仁 t2时窗内Q吸收模型道的振幅谱除以扫描相关子波的振幅谱可以求得Q 吸收的振幅谱Q(f),根据以下福特曼模型用公式可以求出Q吸收滤波因子，把扫描信号^)与Q吸收滤波因子A(t)褶积可求得Q吸雄)=exp-§1/1收扫描信号s'W = ^)=M(0，式中s(O为含Q吸收的扫描信号；^)为可控震源扫描信号；A0为Q吸收滤波因子；将含Q吸收的扫描信号s (o与可控震源扫描信号W)相关得到模拟Q吸 收的可控震源相关子波《々)'(0^(0，式中《力)为模拟Q吸收的可控震源相关子波；s(O为含Q吸收的扫 描信号；W)为可控震源扫描信号；；为相关符号；4) 用以下公式求出反褶积因子h(t), f = J )*/^)-^)]2 ，式中; ha)为棘的反褶棚子；将h(t)应用于对应道集，即将反褶积因子h(t)与对应道集的实际地震数据进行褶积，完成在Q吸收条件下的地表一致性可控震源反褶积和相位 转换；5) 采用上述方法可以得到消除可控震源和炸药震源激发产生的时间和相位差异后，浅、中、深层波形连续性和一致性明显改善的可控震源和炸 药震源激发的成像剖面。本专利技术能自动提取地表一致性Q吸收量，能模拟地表一致性Q吸收可 控震源相关子波，并采用自适应反褶积因子长度和自适应白噪。①在Q吸收条件下，可实现提高分辨率和相位转换的目的，使炸药震 源和可控震源数据能得到较好的闭合；②具有地表一致性统计的自适应能 力，可实现可控震源的地表一致性反褶积处理，使可控震源地震数据具有 统计意义下的叠前相对保持振幅和波形的能力；③能进一步拓宽可控震源 数据频带的范围，使可控震源勘探数据可以进一步提高分辨率；④对于含 有较大近地表吸收影响的可控震源初至波，通过处理得以明显改善，从而 解决了可控震源的折射波静校正的拾取问题。本专利技术使地震资料的处理质量明显提高，处理的地震资料成果提供的 井位钻井成功率明显提高。 附图说明图1是无吸收的可控震源扫描信号子波和频谱； 图2是含Q吸收的可控震源扫描信号子波和频谱； 图3是实际可控震源炮集数据； 图4是不同震源激发的叠加剖面；图5是自适应可控震源地表一致性反褶积处理过程图； 图6是可控震源反褶积处理前后的理论炮集分析； 图7是可控震源反褶积处理前后的实际炮集数据； 图8是可控震源反褶积处理前后叠加剖面；具体实施例方式本专利技术提供以下可控震源自适应地表一致性反褶积的方法，具体包含 以下步骤1、 使用可控震源激发地震波；使用可控震源激发地震波并用仪器记录，可控震源机械振动产生的扫 描信号如图l所示。其中，图l (a)为不含Q吸收可控震源扫描信号，(b) 是相关子波，(c)是子波的频谱。图2 (a)为含大地Q吸收条件下可控震 源扫描信号，(b)是相关子波，(c)是子波的频谱。比较图1和图2中的 子波差异可以看出不含Q吸收条件下的子波是零相位子波，而在Q吸收条 件下的子波变为混合相位子波。从频谱中也可以看出Q吸收作用的影响。图3是实际含Q吸收条件下的可控震源激发的实际炮集数据，从图中 可以明显看出由于可控震源的混合相位子波产生的初至前的子波旁瓣延 续，从而造成可控震源数据初至不清。图4是同一条炸药震源和可控震源 激发的地震剖面，从图中可以看出两种不同形式震源激发的地震数据间存 在明显的闭合差。2、 求地表一致性Q吸收模型道振幅谱；对步骤1得到的可控震源地震记录，选择L t2的时窗进行Q吸收分析， 选取时窗原则类似其他反褶积方法，尽量减少规则干扰波影响。在本例中 科选择图3中1000-2000毫秒作为分析时窗，在仁 t2的时窗内，对时窗内 地震数据做傅氏变换，对于时窗内地震数据包含的噪音可采用频率域振幅 中值滤波消除(频率域振幅中值滤波即将某一炮点或检波点各个地震道的 傅氏变换结果按采样点排列后取中值，即可获得一个炮点或检波点模型道的振幅谱)。通过以上处理，即可获得某一炮点或检波点在仁 t2时窗内Q 吸收模型道的振幅谱。3、求取Q吸收模拟子波；通过步骤2获得某一炮点或检波点在仁 t2时窗内Q吸收模型道的振幅谱后，对该炮的扫描信号做傅氏变换，得到扫描相关子波的振幅谱(扫描信号通常记录在地震数据的辅助道内)。然后用ti 本文档来自技高网...

【技术保护点】
可控震源自适应地表一致性反褶积的方法，其特征在于具体步骤是：　１）使用常规的可控震源激发并记录地震波数据；　２）采用可控震源地震记录并在选定的时窗进行吸收分析，获得炮点或检波点在ｔ↓［１］～ｔ↓［２］时窗内Ｑ吸收模型道的振幅谱； 　３）对炮的扫描信号做傅氏变换，得到扫描相关子波的振幅谱，然后用ｔ↓［１］～ｔ↓［２］时窗内Ｑ吸收模型道的振幅谱除以扫描相关子波的振幅谱可以求得Ｑ吸收的振幅谱Ｑ（ｆ），根据以下福特曼模型用公式可以求出Ｑ吸收滤波因子Ａ（ｔ）＝ｅｘｐ［－ π／Ｑ｜ｆ｜］，把扫描信号ｓ（ｔ）与Ｑ吸收滤波因子Ａ（ｔ）褶积可求得Ｑ吸收扫描信号：ｓ′（ｔ）＝ｓ（ｔ）＊Ａ（ｔ），　式中：ｓ′（ｔ）为含Ｑ吸收的扫描信号；ｓ（ｔ）为可控震源扫描信号；Ａ（ｔ）为Ｑ吸收滤波因子；　将含Ｑ吸收的扫描 信号ｓ′（ｔ）与可控震源扫描信号ｓ（ｔ）相关得到模拟Ｑ吸收的可控震源相关子波Ｒ′↓［ｓ↓［ｉ］］（ｔ）＝ｔ′（ｔ）＊ｓ（ｔ），　式中：Ｒ′↓［ｓ↓［ｉ］］（ｔ）为模拟Ｑ吸收的可控震源相关子波；ｓ′（ｔ）为含Ｑ吸收的扫描信号；ｓ（ｔ）为 可控震源扫描信号；＊为相关符号；　４）用以下公式求出反褶积因子ｈ（ｔ），ε＝∑［Ｒ′↓［ｓ↓［ｉ］］（ｔ）＊ｈ（ｔ）－ｂ（ｔ）］↑［２］，式中：＊＊＊；ｈ（ｔ）为待求的反褶积因子；　将ｈ（ｔ）应用于对应道集，即将反褶积因子ｈ（ｔ ）与对应道集的实际地震数据进行褶积，完成在Ｑ吸收条件下的地表一致性可控震源反褶积和相位转换；　５）采用上述方法可以得到消除可控震源和炸药震源激发产生的时间和相位差异后，浅、中、深层波形连续性和一致性明显改善的可控震源和炸药震源激发的成 像剖面。...

【技术特征摘要】
1、可控震源自适应地表一致性反褶积的方法，其特征在于具体步骤是1)使用常规的可控震源激发并记录地震波数据；2)采用可控震源地震记录并在选定的时窗进行吸收分析，获得炮点或检波点在t1～t2时窗内Q吸收模型道的振幅谱；3)对炮的扫描信号做傅氏变换，得到扫描相关子波的振幅谱，然后用t1～t2时窗内Q吸收模型道的振幅谱除以扫描相关子波的振幅谱可以求得Q吸收的振幅谱Q(f)，根据以下福特曼模型用公式可以求出Q吸收滤波因子<maths id=math0001 num=0001 ><math><![CDATA[ <mrow><mi>A</mi><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac> <mi>&pi;</mi> <mi>Q</mi></mfrac><mo>|</mo><mi>f</mi><mo>|</mo><mo>]</mo><mo>,</mo> </mrow>]]></math> id=icf0001 file=A2008101036960002C1.tif wi=35 he=11 top= 104 left = 22 img-content=drawing img-format=tif orientation=portrait inline=yes/></maths>把扫描信号s(t)与Q吸收滤波因子A(t)褶积可求得Q吸收扫描信号s′(t)＝s(t)*A(t)，式中s′(t)为含Q吸收的扫描信号；s(t)为可控震源扫描信号；A(t)为Q吸收滤波因子；将含Q吸收的扫描信号s′(t)与可控震源扫描信号s(t)相关得到模拟Q吸收的可控震源相关子波<maths id=math0002 num=0002 ><math><![CDATA[ <mrow><msubsup> <mi>R</mi> <msub><mi>s</mi><mi>i</mi> </msub> <mo>&prime;</mo></msubsup><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><mo>=</mo><msup> <mi>s</mi> <mo>&prime;</mo></msup><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><mover> <mo>*</mo> <mo>&OverBar;</mo></mover><mi>s</mi><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><mo>,</mo> </mrow>]]></math> id=icf0002 file=A2008101036960002C2.tif wi=31 he=6 top= 164 left = 75 img-c...

【专利技术属性】
技术研发人员：凌云，高军，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：13[中国|河北]