CN112379436A - 混合震源激发地震资料静校正量计算方法及装置 - Google Patents
混合震源激发地震资料静校正量计算方法及装置 Download PDFInfo
- Publication number
- CN112379436A CN112379436A CN202011201776.1A CN202011201776A CN112379436A CN 112379436 A CN112379436 A CN 112379436A CN 202011201776 A CN202011201776 A CN 202011201776A CN 112379436 A CN112379436 A CN 112379436A
- Authority
- CN
- China
- Prior art keywords
- shot
- source
- arrival
- seismic
- explosive
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012937 correction Methods 0.000 title claims abstract description 87
- 230000003068 static effect Effects 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000002360 explosive Substances 0.000 claims abstract description 64
- 238000004587 chromatography analysis Methods 0.000 claims abstract description 15
- 230000005284 excitation Effects 0.000 claims description 30
- 238000004590 computer program Methods 0.000 claims description 17
- 230000002159 abnormal effect Effects 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000010304 firing Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 15
- 238000012545 processing Methods 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 12
- 238000003860 storage Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种混合震源激发地震资料静校正量计算方法及装置,所述的方法包括:对炸药震源单炮和可控震源单炮的初至波峰拾取;根据预设炮检距间隔选取预设数量的单炮对的初至波峰,对可控震源单炮的初至波峰进行校正,其中,所述单炮对为对应位置的炸药震源单炮和可控震源单炮;根据校正后的可控震源单炮和炸药震源单炮的初至波峰拾取进行初至层析反演近地表结构模型,得到静校正量,本发明可提高混合震源激发地震资料的静校正量计算精度,从而改善地震资料的成像效果。
Description
技术领域
本发明涉及石油勘探技术领域,具体的讲是一种混合震源激发地震资料静校正量计算方法及装置。
背景技术
层析静校正技术采用了非线性的层析反演方法,理论上可以真实地反演出近地表速度模型,结合野外测量所获得的近地表地层厚度,从而得到资料处理中所需的层析静校正量。在实际应用过程中,层析静校正量的计算主要依赖于初至时间拾取的精度、炮检距范围的选取等因素。但是,在山地、沙漠、戈壁、黄土塬和大型城矿区等复杂地表条件的地震勘探中,为了满足地震原始单炮尽可能均匀分布的采集技术要求,经常需要采用炸药震源和可控震源联合激发的野外采集施工方式(有的区域适合炸药震源施工,有的区域适合可控震源施工)。针对这种混合震源激发采集的地震资料,采用统一的层析反演方法不仅会导致反演精度明显降低,有时甚至会导致反演的近地表速度模型产生畸变,严重影响静校正量计算的精度,进而导致部分区域的构造起伏形态与钻井信息不吻合,地震资料成像效果变差,对后续的精细储层预测、低幅度构造圈闭的落实增加了难度。
发明内容
本发明的一个目的在于提供一种混合震源激发地震资料静校正量计算方法,提高了混合震源激发地震资料的静校正量计算精度,从而改善地震资料的成像效果。本发明的另一个目的在于提供一种混合震源激发地震资料静校正量计算装置。本发明的再一个目的在于提供一种计算机设备。本发明的还一个目的在于提供一种可读介质。
为了达到以上目的,本发明一方面公开了一种混合震源激发地震资料静校正量计算方法,所述的方法包括:
对炸药震源单炮和可控震源单炮的初至波峰拾取;
根据预设炮检距间隔选取预设数量的单炮对的初至波峰,对可控震源单炮的初至波峰进行校正,其中,所述单炮对为对应位置的炸药震源单炮和可控震源单炮;
根据校正后的可控震源单炮和炸药震源单炮的初至波峰拾取进行初至层析反演近地表结构模型,得到静校正量。
优选的,进一步包括预先对可控震源单炮进行最小相位化校正。
优选的,所述预设炮检距的选择范围为200~500米。
优选的,所述对可控震源单炮的初至波峰进行校正具体包括:
确定每个单炮对初至波峰的初至时差;
求取预设数量单炮对的初至时差的平均值;
根据所述平均值对所述可控震源单炮的初至波峰进行校正。
优选的,进一步包括在求取预设数量单炮对的初至时差的平均值,之前:
剔除所有初至时差中的异常值。
优选的,所述预设数量大于等于9。
优选的,所述对应位置的炸药震源单炮和可控震源单炮的单炮对为同一位置或相邻位置的炸药震源单炮和可控震源单炮。
本发明还公开了一种混合震源激发地震资料静校正量计算装置,包括:
初至拾取模块,用于对炸药震源单炮和可控震源单炮的初至波峰拾取;
初至校正模块,用于根据预设炮检距间隔选取预设数量的单炮对的初至波峰,对可控震源单炮的初至波峰进行校正,其中,所述单炮对为对应位置的炸药震源单炮和可控震源单炮;
静校正量求解模块,用于根据校正后的可控震源单炮和炸药震源单炮的初至波峰拾取进行初至层析反演近地表结构模型,得到静校正量。
优选的,所述初至校正模块具体用于确定每个单炮对初至波峰的初至时差;求取预设数量单炮对的初至时差的平均值;根据所述平均值对所述可控震源单炮的初至波峰进行校正。
本发明还公开了一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,
所述处理器执行所述程序时实现如上所述方法。
本发明还公开了一种计算机可读介质,其上存储有计算机程序,
该程序被处理器执行时实现如上所述方法。
本发明通过对混合震源中的炸药震源单炮和可控震源单炮激发后的初至波峰进行拾取,并根据预设炮检距间隔选取预设数量的单炮对的初至波峰以对可控震源单炮的初至波峰进行校正。从而,本发明通过对可控震源单炮的初至波峰时间进行一系列的处理,使可控震源单炮的初至波峰达到与炸药震源单炮的初至波峰时间相同,解决了以往直接利用混合震源激发采集原始单炮进行初至层析反演所导致的静校正量精度低且容易产生畸变的问题,进一步提高了混合震源激发采集地震资料静校正量计算的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明混合震源激发地震资料静校正量计算方法一个具体实施例的流程图;
图2为本发明混合震源激发地震资料静校正量计算方法一个具体实施例S200的流程图;
图3(a)为现有技术中混源激发采集资料常规层析静校正量分布图;
图3(b)为本发明混合震源激发地震资料静校正量计算方法一个具体实施例混源激发采集资料新方法层析静校正量分布图;
图4(a)为现有技术中混源激发采集资料常规层析静校正叠加剖面;
图4(b)为本发明混合震源激发地震资料静校正量计算方法一个具体实施例混源激发采集资料新方法层析静校正叠加剖面;
图5为本发明混合震源激发地震资料静校正量计算装置一个具体实施例的结构图;
图6示出适于用来实现本发明实施例的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
现有技术中,通过炸药震源和可控震源混合震源激发采集得到的地震资料由于不同的震源方式而存在时间偏差的问题。具体的,由于激发方式的改变,导致地表相同或相邻激发点位置的可控震源单炮和炸药震源单炮的原始单炮初至到达时间存在明显差异,可控震源单炮为地面激发、地面接收的单炮记录,炸药震源单炮为地下井底激发、地面接收的单炮记录。由于二者激发点深度的不同,采用混合震源激发得到的地震资料进行层析反演得到的近地表速度模型可能产生畸形,严重影响静校正量的计算精度。为了解决现有技术中存在的问题,本发明提出了一种混合震源激发地震资料静校正量计算方法,利用该方法可以精确计算混合震源激发采集地震资料的静校正量,具有较好的应用效果。
根据本发明的一个方面,本实施例公开了一种混合震源激发地震资料静校正量计算方法。如图1所示,本实施例中,所述方法包括:
S100:对炸药震源单炮和可控震源单炮的初至波峰拾取。
S200:根据预设炮检距间隔选取预设数量的单炮对的初至波峰,对可控震源单炮的初至波峰进行校正,其中,所述单炮对为对应位置的炸药震源单炮和可控震源单炮。
S300:根据校正后的可控震源单炮和炸药震源单炮的初至波峰拾取进行初至层析反演近地表结构模型,得到静校正量。
本发明通过对混合震源中的炸药震源单炮和可控震源单炮激发后的初至波峰进行拾取,并根据预设炮检距间隔选取预设数量的单炮对的初至波峰以对可控震源单炮的初至波峰进行校正。从而,本发明通过对可控震源单炮的初至波峰时间进行一系列的处理,使可控震源单炮的初至波峰达到与炸药震源单炮的初至波峰时间相同,解决了以往直接利用混合震源激发采集原始单炮进行初至层析反演所导致的静校正量精度低且容易产生畸变的问题,进一步提高了混合震源激发采集地震资料静校正量计算的精度。
在优选的实施方式中,所述方法进一步包括预先对可控震源单炮进行最小相位化校正。可以理解的是,由于炸药震源单炮和可控震源单炮极性不同,所以首先需要将可控震源单炮进行最小相位化校正,校正完成后分别对两种震源单炮进行初至波峰拾取,提高地震资源的准确性。
在优选的实施方式中,如图2所示,所述S200中对可控震源单炮的初至波峰进行校正具体包括:
S210:确定每个单炮对初至波峰的初至时差。
S220:求取预设数量单炮对的初至时差的平均值。
S230:根据所述平均值对所述可控震源单炮的初至波峰进行校正。
可以理解的是,由于炸药震源单炮和可控震源单炮激发点深度的不同,需要对其中的一个震源的初至波峰进行时间校正。从而,通过时间校正较好地解决了炸药震源单炮与可控震源单炮的地震资料之间由于激发井深造成的初至到达时间的差异,提高了混合震源激发地震资料的静校正量计算精度,从而改善地震资料的成像效果。
其中,可以理解的是,可以首先对炸药震源单炮和可控震源单炮分别进行初至波峰的拾取,初至拾取结束后统计出工区同一位置(或相邻位置)的炸药震源单炮和可控震源单炮的单炮对,并根据预设炮检距间隔选取预设数量的单炮对的初至波峰。
具体的,可确定每个单炮对初至波峰的初至时差,求取预设数量单炮对的初至时差的平均值,并根据所述平均值对所述可控震源单炮的初至波峰进行校正。例如,在一个具体例子中,第i个单炮对的初至时差计为ki,i=1,2,…n,n为单炮对的预设数量。然后,预设数量单炮对的初至时差的平均值Kn=(k1+k2+......kn)/n。最后,将可控震源单炮的初至波峰时间Ti按照时差值Kn进行校正Ti′=Ti-Kn,以达到与炸药震源单炮初至波峰的时间基本一致的目的。通过以上一系列处理可以消除可控震源单炮与炸药震源单炮之间由于激发井深造成的初至时间差异。最后,把时差校正后的可控震源单炮初至与已有的炸药震源单炮初至合并到一起进行初至层析反演近地表结构模型,获得最终的静校正量用于混合震源激发采集地震资料的基准面静校正处理。
在优选的实施方式中,S200根据预设炮检距间隔选取预设数量的单炮对的初至波峰中,所述预设炮检距的选择范围优选的为200~500米。所述预设数量优选的大于等于9。则挑选出m对有代表性的单炮(m≥9),然后在每对单炮中按一定炮检距间隔(200-500米)挑选出n组相等炮检距的地震道进行初至时差统计(n≥9),求取平均时差值Kn后,将可控震源的初至时间按照时差值Kn进行校正,以达到与炸药震源原始单炮初至时间一致的目的。当然,在其他实施方式中,预设数量和预设炮检距也可以根据实际情况选择其他的数值,本发明对此并不作限定。
在优选的实施方式中,所述方法进一步包括在S220求取预设数量单炮对的初至时差的平均值,之前:
S211:剔除所有初至时差中的异常值。
其中,可以理解的是,通过对预设数量单炮对的初至时差的平均值进行异常值处理,可以防止引入异常值而导致的平均值计算出现误差的问题。例如,在一个具体例子中,可剔除ki≥2K的异常值,其中,K为常值,可根据计算精度或按经验确定。
下面通过一个具体例子来对本发明作进一步的说明。在该具体例子中,首先需要对所有原始单炮进行初至拾取,然后将可控震源原始单炮拾取的初至进行时差校正,使同一位置(或相邻位置)的炸药震源单炮和可控震源单炮初至一致,然后进行统一初至层析反演,获得基准面静校正所需的层析静校正量,具体操作步骤如下:
S1000:分选出炸药震源单炮和可控震源单炮,由于炸药震源单炮和可控震源单炮极性不同,所以首先需要将可控震源单炮进行最小相位化校正,校正完成后分别对两种震源单炮进行初至波峰拾取。
S2000:统计出工区同一位置(或相邻位置)的炸药震源单炮和可控震源单炮(以下称为“单炮对”),并且按照一定间隔、尽可能均匀的原则挑选出m组有代表性的“单炮对”(m≥9)。
S3000:在每组“单炮对”中按一定炮检距间隔(200-500米)挑选出n组相等炮检距的地震道进行初至时差统计(n≥9),每组相等炮检距地震道的初至时差计为ki,ki应该在一个常量K上下变动。
S4000:对每组的初至时差ki进行分析,剔除异常值后(若ki≥2K即为异常值),求取平均值Kn=(k1+k2+......kn)/n。
S5000:将可控震源的初至时间按照时差值Kn进行校正,以达到与炸药震源原始单炮初至时间基本一致。
S6000:把时差校正后的可控震源单炮初至与已有的炸药震源单炮初至合并到一起,统一进行初至层析反演近地表结构模型,获得最终的静校正量,用于混合震源激发采集地震资料的基准面静校正处理。
本发明的方法不仅能够提高静校正量计算精度,避免层析反演建立的近地表结构模型产生畸变,而且能够改善地震资料的叠加成像效果,与常规方法应用效果对比见图3(a)、图3(b)、图4(a)和图4(b)所示。图3(a)是混源激发采集资料常规层析静校正量分布图,图3(b)是混源激发采集资料新方法层析静校正量分布图,图3(b)中炸药震源和可控震源交界处静校正量突变消失,对比可见利用新方法获得的静校正量没有发生畸变,更接近于地表真实情况。图4(a)是混源激发采集资料常规层析静校正叠加剖面,图4(b)是混源激发采集资料新方法层析静校正叠加剖面,图4(b)圈中同相轴的连续性更好,信噪比提升,对比可见利用新方法获得的叠加剖面成像效果改进明显。
综上,该方法通过对可控震源单炮初至进行一系列的处理达到与炸药震源单炮初至相同,解决了以往直接利用混合震源激发采集原始单炮进行初至层析反演所导致的静校正量精度低且容易产生畸变的问题,进一步提高了混合震源激发采集地震资料静校正量计算的精度。
基于相同原理,本实施例还公开了一种混合震源激发地震资料静校正量计算装置。如图5所示,本实施例中,所述装置包括初至拾取模块11、初至校正模块12和静校正量求解模块13。
其中,初至拾取模块11用于对炸药震源单炮和可控震源单炮的初至波峰拾取。
初至校正模块12用于根据预设炮检距间隔选取预设数量的单炮对的初至波峰,对可控震源单炮的初至波峰进行校正,其中,所述单炮对为对应位置的炸药震源单炮和可控震源单炮。
静校正量求解模块13用于根据校正后的可控震源单炮和炸药震源单炮的初至波峰拾取进行初至层析反演近地表结构模型,得到静校正量。
本发明通过对混合震源中的炸药震源单炮和可控震源单炮激发后的初至波峰进行拾取,并根据预设炮检距间隔选取预设数量的单炮对的初至波峰以对可控震源单炮的初至波峰进行校正。从而,本发明通过对可控震源单炮的初至波峰时间进行一系列的处理,使可控震源单炮的初至波峰达到与炸药震源单炮的初至波峰时间相同,解决了以往直接利用混合震源激发采集原始单炮进行初至层析反演所导致的静校正量精度低且容易产生畸变的问题,进一步提高了混合震源激发采集地震资料静校正量计算的精度。
在优选的实施方式中,所述初至校正模块具体用于确定每个单炮对初至波峰的初至时差;求取预设数量单炮对的初至时差的平均值;根据所述平均值对所述可控震源单炮的初至波峰进行校正。
可以理解的是,由于炸药震源单炮和可控震源单炮激发点深度的不同,需要对其中的一个震源的初至波峰进行时间校正。从而,通过时间校正较好地解决了炸药震源单炮与可控震源单炮的地震资料之间由于激发井深造成的初至到达时间的差异,提高了混合震源激发地震资料的静校正量计算精度,从而改善地震资料的成像效果。
其中,可以理解的是,可以首先对炸药震源单炮和可控震源单炮分别进行初至波峰的拾取,初至拾取结束后统计出工区同一位置(或相邻位置)的炸药震源单炮和可控震源单炮的单炮对,并根据预设炮检距间隔选取预设数量的单炮对的初至波峰。
由于该装置解决问题的原理与以上方法类似,因此本装置的实施可以参见方法的实施,在此不再赘述。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备,具体的,计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上所述的由客户端执行的方法,或者,所述处理器执行所述程序时实现如上所述的由服务器执行的方法。
下面参考图6,其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图6所示,计算机设备600包括中央处理单元(CPU)601,其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中,还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605:包括键盘、鼠标等的输入部分606;包括诸如阴极射线管(CRT)、液晶反馈器(LCD)等以及扬声器等的输出部分607;包括硬盘等的存储部分608;以及包括诸如LAN卡,调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器610上,以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地,根据本发明的实施例,上文参考流程图描述的过程可以被实现为计算机软件程序。例如,本发明的实施例包括一种计算机程序产品,其包括有形地包含在机器可读介质上的计算机程序,所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分609从网络上被下载和安装,和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (11)
1.一种混合震源激发地震资料静校正量计算方法,其特征在于,所述的方法包括:
对炸药震源单炮和可控震源单炮的初至波峰拾取;
根据预设炮检距间隔选取预设数量的单炮对的初至波峰,对可控震源单炮的初至波峰进行校正,其中,所述单炮对为对应位置的炸药震源单炮和可控震源单炮;
根据校正后的可控震源单炮和炸药震源单炮的初至波峰拾取进行初至层析反演近地表结构模型,得到静校正量。
2.如权利要求1所述的混合震源激发地震资料静校正量计算方法,其特征在于,进一步包括预先对可控震源单炮进行最小相位化校正。
3.如权利要求1所述的混合震源激发地震资料静校正量计算方法,其特征在于,所述预设炮检距的选择范围为200~500米。
4.如权利要求1所述的混合震源激发地震资料静校正量计算方法,其特征在于,所述对可控震源单炮的初至波峰进行校正具体包括:
确定每个单炮对初至波峰的初至时差;
求取预设数量单炮对的初至时差的平均值;
根据所述平均值对所述可控震源单炮的初至波峰进行校正。
5.如权利要求4所述的混合震源激发地震资料静校正量计算方法,其特征在于,进一步包括在求取预设数量单炮对的初至时差的平均值,之前:
剔除所有初至时差中的异常值。
6.如权利要求1所述的混合震源激发地震资料静校正量计算方法,其特征在于,所述预设数量大于等于9。
7.如权利要求1所述的混合震源激发地震资料静校正量计算方法,其特征在于,所述对应位置的炸药震源单炮和可控震源单炮的单炮对为同一位置或相邻位置的炸药震源单炮和可控震源单炮。
8.一种混合震源激发地震资料静校正量计算装置,其特征在于,包括:
初至拾取模块,用于对炸药震源单炮和可控震源单炮的初至波峰拾取;
初至校正模块,用于根据预设炮检距间隔选取预设数量的单炮对的初至波峰,对可控震源单炮的初至波峰进行校正,其中,所述单炮对为对应位置的炸药震源单炮和可控震源单炮;
静校正量求解模块,用于根据校正后的可控震源单炮和炸药震源单炮的初至波峰拾取进行初至层析反演近地表结构模型,得到静校正量。
9.如权利要求8所述的混合震源激发地震资料静校正量计算装置,其特征在于,所述初至校正模块具体用于确定每个单炮对初至波峰的初至时差;求取预设数量单炮对的初至时差的平均值;根据所述平均值对所述可控震源单炮的初至波峰进行校正。
10.一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,
所述处理器执行所述程序时实现如权利要求1-7任一项所述方法。
11.一种计算机可读介质,其上存储有计算机程序,其特征在于,
该程序被处理器执行时实现如权利要求1-7任一项所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011201776.1A CN112379436A (zh) | 2020-11-02 | 2020-11-02 | 混合震源激发地震资料静校正量计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011201776.1A CN112379436A (zh) | 2020-11-02 | 2020-11-02 | 混合震源激发地震资料静校正量计算方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112379436A true CN112379436A (zh) | 2021-02-19 |
Family
ID=74577838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011201776.1A Pending CN112379436A (zh) | 2020-11-02 | 2020-11-02 | 混合震源激发地震资料静校正量计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112379436A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113514874A (zh) * | 2021-03-30 | 2021-10-19 | 中国石油化工股份有限公司 | 不同类型地震采集系统之间的时差校正装置、方法及系统 |
CN113945980A (zh) * | 2021-09-30 | 2022-01-18 | 国能神东煤炭集团有限责任公司 | 一种初至时间校正方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
US20100157733A1 (en) * | 2008-12-23 | 2010-06-24 | Mark Daly | Compensating seismic data for source variations |
CN102478666A (zh) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | 一种提高可控震源地震记录初至拾取精度的方法 |
CN103777239A (zh) * | 2012-10-17 | 2014-05-07 | 中国石油化工股份有限公司 | 一种提高不同震源资料相位一致性的相位旋转方法 |
CN105319587A (zh) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | 一种基于初至子波提取的可控震源整形滤波方法 |
CN108008447A (zh) * | 2017-10-23 | 2018-05-08 | 中国石油天然气集团公司 | 一种炮集记录处理方法、装置及计算机存储介质 |
-
2020
- 2020-11-02 CN CN202011201776.1A patent/CN112379436A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
US20100157733A1 (en) * | 2008-12-23 | 2010-06-24 | Mark Daly | Compensating seismic data for source variations |
CN102478666A (zh) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | 一种提高可控震源地震记录初至拾取精度的方法 |
CN103777239A (zh) * | 2012-10-17 | 2014-05-07 | 中国石油化工股份有限公司 | 一种提高不同震源资料相位一致性的相位旋转方法 |
CN105319587A (zh) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | 一种基于初至子波提取的可控震源整形滤波方法 |
CN108008447A (zh) * | 2017-10-23 | 2018-05-08 | 中国石油天然气集团公司 | 一种炮集记录处理方法、装置及计算机存储介质 |
Non-Patent Citations (5)
Title |
---|
刘占族;曹孟起;孙军晓;陈院生;晏丰: "煤层气地震勘探中的静校正方法", 《石油地球物理勘探》, vol. 47, no. 03, 30 June 2012 (2012-06-30), pages 425 - 429 * |
刘恋;岳洪波: "可控震源与炸药震源联合采集在煤炭勘探上的应用", 《能源技术与管理》, vol. 45, no. 03, 30 June 2020 (2020-06-30), pages 191 - 193 * |
梁鸿贤: "基于初至信息的可控震源和炸药震源地震资料匹配滤波技术", 《石油物探》, vol. 57, no. 02, 30 April 2018 (2018-04-30), pages 231 - 236 * |
王敏: "匹配滤波技术在江汉平原海相大剖面拼接处理中的应用", 《工程地球物理学报》, vol. 14, no. 02, 30 April 2017 (2017-04-30), pages 165 - 173 * |
闫忠;刘安成: "沼泽地区地震资料处理技术", 《石油地球物理勘探》, vol. 42, no. 1, 31 August 2007 (2007-08-31), pages 40 - 45 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113514874A (zh) * | 2021-03-30 | 2021-10-19 | 中国石油化工股份有限公司 | 不同类型地震采集系统之间的时差校正装置、方法及系统 |
CN113514874B (zh) * | 2021-03-30 | 2024-02-02 | 中国石油化工股份有限公司 | 不同类型地震采集系统之间的时差校正装置、方法及系统 |
CN113945980A (zh) * | 2021-09-30 | 2022-01-18 | 国能神东煤炭集团有限责任公司 | 一种初至时间校正方法及系统 |
CN113945980B (zh) * | 2021-09-30 | 2024-02-09 | 国能神东煤炭集团有限责任公司 | 一种初至时间校正方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112379436A (zh) | 混合震源激发地震资料静校正量计算方法及装置 | |
US11209563B2 (en) | Data optimization method and integral prestack depth migration method | |
CN109765615A (zh) | 一种地层品质因子反演方法及装置 | |
CN104155694A (zh) | 一种反射转换横波共检波点叠加剖面的剩余静校正方法 | |
CN106054252B (zh) | 一种叠前时间偏移的方法及装置 | |
CN111983679B (zh) | 基于沉积参数的相控型白云岩储层地震预测方法和装置 | |
CN112764102B (zh) | 地震资料采集方法及装置 | |
CN113568033A (zh) | 一种三维非规则采样地震采集观测系统设计方法和装置 | |
CN109188515A (zh) | 微地震监测裂缝震源点位置计算方法及系统 | |
CN112946736B (zh) | 一种三维观测系统重建方法及系统 | |
CN111999767B (zh) | 起伏地表的偏移成像方法及装置 | |
CN111337973B (zh) | 地震数据重建方法、系统 | |
CN112099091B (zh) | 地震勘探点数据静校正方法及装置 | |
CN111751880B (zh) | 井震数据联合显示方法及装置 | |
CN104880732B (zh) | 一种十字子集的构建方法及装置 | |
CN113093279B (zh) | 转换波静校正方法及装置 | |
CN113138419B (zh) | 提取下行子波和衰减参数的方法、装置 | |
CN112946745B (zh) | 一种基于引力拟合的转换波剩余静校正方法及系统 | |
CN113156494B (zh) | 面元等效覆盖次数的计算方法及装置 | |
CN111596359B (zh) | 一种倾斜地层宽线数据校正成像方法及装置 | |
CN113296153B (zh) | 轴对称介质各向异性参数确定方法及装置 | |
CN113031071A (zh) | 地震波长波长静校正方法及装置 | |
CN114063165B (zh) | 一种三维地震数据的拼接方法及装置 | |
CN112241023A (zh) | 陆上三维炮点二次定位及定量偏移监控方法和装置 | |
CN116299689A (zh) | 构建近地表速度模型的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |