一种快速准确选取可控震源滑动时间参数的方法技术
【技术实现步骤摘要】
一种快速准确选取可控震源滑动时间参数的方法
本专利技术涉及石油勘探工程，是地震勘探中使用可控震源滑动扫描方法生产中快速准确选取可控震源滑动时间参数的方法。
技术介绍
目前，在石油地震勘探中，滑动扫描技术是国际上可控震源高效采集技术中最成熟、最常用的技术之一，该技术已获得了良好的勘探效果和经济效益。滑动扫描是指使用多组可控震源同时进行作业，在该技术应用中，谐波干扰是客观存在的,在处理时需要专门的压制技术，在野外通过设计合理的滑动时间也可以有效减弱谐波干扰。滑动时间是可控震源相邻两次震动时间间隔，滑动时间原则上大于地震记录长度即可，与交替扫描相比，相邻的两次振动(扫描)，突破了第二次扫描必须等第一次扫描记录结束才能开始的限制，压缩了相邻两次扫描的间隔时间，可提高生产效率。滑动时间参数会影响原始资料品质，目前可控震源滑动扫描施工中，确定滑动时间往往是野外应用2套或多套震源，按照滑动扫描方式分别按照不同的滑动时间进行数据采集，室内对上述数据进行分析对比，以确定合适的滑动时间，一般是通过野外大量的点上试验来确定滑动时间参数，工作量大，资料分析周期长，如果施工区域复杂，生产过程中需要考核试验时，将会进一步影响生产效率，同时，确定的滑动时间参数也具有一定的偶然性。
技术实现思路
本专利技术目的在于提供一种可提高生产效率，提高资料品质的快速准确选取可控震源滑动时间参数的方法。本专利技术通过以下步骤实现：1)根据已知的可控震源扫描长度和起止扫描频率，计算出避开二次谐波所要求的最小滑动时间Ti；步骤1)所述的最小滑动时间Ti采用以下公式计算：k(fE-fB)T---(1)]]>其中：T是扫描长度，fB、fE分别是起止扫描频率；k为常数；所述的k为常数，k＝2。2)野外现场采集可控震源原始未相关数据；步骤2)所述的野外现场采集可控震源原始未相关数据是利用可控震源扫描信号S(t)，在震源1和震源2位置，利用相同排列记录可控震源原始未相关数据Raw1(t)和Raw2(t)。3)将原始数据Raw2(t)在时间域向后延迟时间i，得到延迟后的记录Raw2(ti),每次取不同的时间i，获得不同模拟滑动时间i的记录Raw2(ti)；所述的时间i大于步骤1)中估算的最小滑动时间Ti，小于或等于扫描长度T，i取整数，根据实际情况i取3-5个值。4)将由步骤2)得到的Raw1(t)和由步骤3)得到的不同模拟滑动时间i的记录Raw2(ti)分别进行垂直叠加，得到对应模拟滑动时间i的数据Raw(ti)；步骤4)所述的垂直叠加采用以下公式：Raw(ti)＝Raw1(t)+Raw2(ti)(2)；5)将Raw(ti)与可控震源扫描信号S(t)进行互相关，得到模拟滑动时间i的地震记录R(ti)；6)根据对模拟记录R(ti)谐波干扰的大小和范围，选择干扰能量弱和范围小的模拟滑动时间，同时，在记录上读取谐波干扰的起始时间，选择干扰起始时间大于目的层时间的滑动时间，确定可控震源滑动时间参数i。本专利技术实现了可控震源滑动扫描方法施工中滑动时间的快速、准确选择，减少了野外试验工作量，降低了野外试验费用，缩短了试验周期，确保滑动时间参数选取的合理性。附图说明图1是采集可控震源原始未相关记录的观测系统示意图；图2是用于模拟分析的不同位置的相关记录；图3是模拟不同滑动时间的相关记录；图4是模拟不同滑动时间的记录时频图，A)是模拟滑动时间8s的记录时频图；B)是模拟滑动时间10s的记录时频图；C)是模拟滑动时间12s的记录时频图；D)是模拟滑动时间14s的记录时频图。具体实施方式以下结合附图详细说明本专利技术。本专利技术采用以下步骤实现：1)野外通过一条适当长的二维排列接收平面上不同位置的可控震源常规未相关记录，根据已知的可控震源扫描长度T和起止扫描频率fB、fE，由公式(1)计算出避开二次谐波所要求的最小滑动时间Ti；式中的k为常数，k＝2。2)野外现场采集可控震源原始未相关数据步骤2)所述的野外现场采集可控震源原始未相关数据是根据已确定的可控震源参数，利用可控震源扫描信号S(t)，在震源1和震源2位置，利用相同排列记录可控震源原始未相关数据Raw1(t)和Raw2(t)。见图1，其中“△”代表接收排列,“●”代表可控震源位置。3)将原始数据Raw2(t)在时间域向后延迟时间i，该步骤可以在CGG、ProMax、Ω、GeoEast等处理系统上，利用野外静校正应用模块实现。如在GeoEast处理系统上，利用静校正应用模块StApply，在“applicationpotions”选项里面选择“constantcorrection”，然后在“correctionconstant”参数里填写常数校正时间-i，单位为ms,生成SEG-Y格式数据，输出文件名Raw2(ti)。对不同的时间i，重复该步骤，得到不同模拟滑动时间i的记录Raw2(ti)。所述的时间i大于步骤1)中估算的最小滑动时间Ti，小于扫描长度T，i取整数，根据实际情况i取3-5个值。4)将由步骤2)得到的Raw1(t)和由步骤3)得到的不同模拟滑动时间i的记录Raw2(ti)分别进行垂直叠加，得到对应模拟滑动时间i的数据Raw(ti)；步骤4)所述的垂直叠加采用以下公式：Raw(ti)＝Raw1(t)+Raw2(ti)(2)；可在CGG、ProMax、Ω、GeoEast等处理系统上，利用多功能叠加模块将Raw1(t)和Raw2(ti)按同道号进行叠加，形成一个新的炮集。如在GeoEast处理系统上，利用叠加模块VersatiStk，在“inputgathertype”中选择“source”作为输入道集类型，在“headerwordforstack”参数中选择“tracenumber”作为叠加道头字。将Raw1(t)和Raw2(ti)按同道号进行叠加，形成一个新的炮集,生成SEG-Y格式数据，输出文件名Raw(ti)。5)将Raw(ti)与可控震源扫描信号S(t)进行互相关，得到模拟滑动时间i的地震记录R(ti)；可在CGG、ProMax、Ω、GeoEast等处理系统上，利用互相关分析模块，将数据Raw(ti)与扫描信号S(t)互相关，生成互相关后的地震记录R(ti)。如在GeoEast处理系统上，在输入数据Raw(ti)上点击右键选择“CorrelationAnalysis”启动GeoSeismicCorrelation交互相关分析模块，此时数据Raw(ti)已经在交互界面的数据窗口下打开，然后需要在信号窗口里面打开扫描信号S(t)，选择交互界面上端的“Correlation”选项里面的“Crosscorrelation”互相关选项，弹出互相关参数定义窗口，选择填写需要计算的地震数据时窗和信号数据时窗等参数，点击“OK”后生成互相关后的地震记录R(ti)，以SEG-Y格式数据输出。分别对每个i值重复步骤3)至步骤5)，得到不同模拟滑动时间i的地震记录R(ti)；6)分析不同的模拟记录R(ti)即可准确选取可控震源滑动时间参数i；对模拟记录R(ti)进行能量和信噪比定量分析统计，确定记录上谐波干扰的大小和范围，选择谐本文档来自技高网...
【技术保护点】
一种快速准确选取可控震源滑动时间参数的方法，特点是通过以下步骤实现：1）根据已知的可控震源扫描长度和起止扫描频率，计算出避开二次谐波所要求的最小滑动时间Ti；2）野外现场采集可控震源原始未相关数据；步骤2）所述的野外现场采集可控震源原始未相关数据是利用可控震源扫描信号S(t)，在震源1和震源2位置，利用相同排列记录可控震源原始未相关数据Raw1(t)和Raw2(t)；3）将原始数据Raw2(t)在时间域向后延迟时间i，得到延迟后的记录Raw2(ti),每次取不同的时间i，获得不同模拟滑动时间i的记录Raw2(ti)；4）将由步骤2）得到的Raw1(t)和由步骤3）得到的不同模拟滑动时间i的记录Raw2(ti)分别进行垂直叠加，得到对应模拟滑动时间i的数据Raw(ti)；步骤4）所述的垂直叠加采用以下公式：Raw(ti)＝Raw1(t)+Raw2(ti)（2）；5）将Raw(ti)与可控震源扫描信号S(t)进行互相关，得到模拟滑动时间i的地震记录R(ti)；6）根据对模拟记录R(ti)谐波干扰的大小和范围，选择干扰能量弱和范围小的模拟滑动时间，同时，在记录上读取谐波干扰的起始时间，选择干...
【技术特征摘要】
1.一种快速准确选取可控震源滑动时间参数的方法，特点是通过以下步骤实现：1)根据已知的可控震源扫描长度和起止扫描频率，计算出避开二次谐波所要求的最小滑动时间Ti；2)野外现场采集可控震源原始未相关数据；野外现场采集可控震源原始未相关数据是利用可控震源扫描信号S(t)，在震源1和震源2位置，利用相同排列记录可控震源原始未相关数据Raw1(t)和Raw2(t)；3)将原始数据Raw2(t)在时间域向后延迟时间i，得到延迟后的记录Raw2(ti),每次取不同的时间i，获得不同模拟滑动时间i的记录Raw2(ti)；4)将由步骤2)得到的Raw1(t)和由步骤3)得到的不同模拟滑动时间i的记录Raw2(ti)分别进行垂直叠加，得到对应模拟滑动时间i的数据Raw(ti)；垂直叠加采用以下公式：Raw(ti)＝Raw1(...
【专利技术属性】
技术研发人员：王新全，唐东磊，汪长辉，孙梦，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。