CN111999767A - 起伏地表的偏移成像方法及装置 - Google Patents
起伏地表的偏移成像方法及装置 Download PDFInfo
- Publication number
- CN111999767A CN111999767A CN202010703848.6A CN202010703848A CN111999767A CN 111999767 A CN111999767 A CN 111999767A CN 202010703848 A CN202010703848 A CN 202010703848A CN 111999767 A CN111999767 A CN 111999767A
- Authority
- CN
- China
- Prior art keywords
- data
- grid point
- horizontal
- relief
- preset
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 77
- 238000012937 correction Methods 0.000 claims abstract description 58
- 238000000034 method Methods 0.000 claims abstract description 44
- 230000003068 static effect Effects 0.000 claims abstract description 35
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000004590 computer program Methods 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 16
- 239000000523 sample Substances 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 7
- 238000013508 migration Methods 0.000 description 24
- 230000005012 migration Effects 0.000 description 24
- 238000010586 diagram Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 16
- 238000004891 communication Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种起伏地表的偏移成像方法及装置,该方法包括:采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。本发明可以实现在起伏地表上的偏移成像,成像精度高。
Description
技术领域
本发明涉及成像处理技术领域,特别涉及一种起伏地表的偏移成像方法及装置。
背景技术
目前,起伏地表的偏移成像处理,即波场延拓偏移处理,一般要求偏移基准面是水平的,且偏移的零点应在激发和接收的地表。在高差较大的复杂地表区,很难同时满足这些要求,这时可用静校正时移的方法把地震资料校正到一个水平基准面上,由于静校正时移方法只作了垂直方向的时移,没有考虑波的传播方向,因此改变了原时间剖面上绕射波的双曲线性质,偏移后收敛较差,故这种方法是近似的,只能在地表与水平统一基准面的高差较小的地区使用。校正量越大,影响成像精度的程度也越大。因此在山地等起伏地表上,用统一水平基准面条件下的波场延拓偏移,就会影响成像的精度。
发明内容
本发明实施例提出一种起伏地表的偏移成像方法,用以实现在起伏地表上的偏移成像,成像精度高,该方法包括:
采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;
基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;
确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;
基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。
本发明实施例提出一种起伏地表的偏移成像装置,用以实现在起伏地表上的偏移成像,成像精度高,该装置包括:
静校正模块,用于采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;
水平叠加剖面数据获得模块,用于基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;
起伏地表因子确定模块,用于确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;
偏移成像模块,用于基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。
本发明实施例还提出了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述起伏地表的偏移成像方法。
本发明实施例还提出了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述起伏地表的偏移成像方法的计算机程序。
在本发明实施例中,采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。在上述过程中,在对叠加剖面数据进行波长延拓成像时,考虑了预设的起伏地表因子,从而可以得到起伏地表的准确的成像结果,提高了成像精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中起伏地表的偏移成像方法的流程图;
图2为本发明实施例中水平叠加剖面数据中网格点坐标的示意图;
图3为目前传统使用的水平地表波场延拓偏移过程示意图;
图4为本发明实施例中起伏地表波场延拓偏移过程示意图;
图5为本发明实施例中起伏地表两层介质的模型;
图6为通过目前工业界常用的方法进行波场延拓偏移处理偏移结果;
图7为采用本发明实施例提出的方法进行起伏地表波场延拓处理的偏移结果;
图8为本发明实施例中起伏地表的偏移成像装置的示意图;
图9为本发明实施例中计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
在本说明书的描述中,所使用的“包含”、“包括”、“具有”、“含有”等,均为开放性的用语,即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施,其中的步骤顺序不作限定,可根据需要作适当调整。
图1为本发明实施例中起伏地表的偏移成像方法的流程图,如图1所示,该方法包括:
步骤101,采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;
步骤102,基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;
步骤103,确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;
步骤104,基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。
在本发明实施例提出的方法中,在对叠加剖面数据进行波长延拓成像时,考虑了预设的起伏地表因子,从而可以得到起伏地表的准确的成像结果,提高了成像精度。
具体实施时,CMP(共中心点数据道集)基准面是通过预先设置好水平统一基面,CMP(共中心点数据道集)每一个道的炮点、检波点到该水平统一基面的时间为该道炮点、检波点的静校正量。将CMP(共中心点数据道集)每一个道的炮点、检波点静校正量相加除以该CMP道集的道数得到的数值称为静校正的低频静校正量,由该值确定的面叫做CMP基准面。该基准是体现地表起伏的一个面。每一道的炮点、检波点的静校正量减去该低频静校正量的值称为该道炮点、检波点的高频校正量。这样每一道炮、检点的静校正量值分成了低频静校正量与高频静校正量两部分。将原始数据加上炮点、检波点的高频校正量来获得CMP参考面数据,在一实施例中,采用CMP基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据,包括:
采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得炮点校正量和检波点校正量;
根据炮点校正量和检波点校正量,获得CMP参考面数据。
在上述实施例中,根据炮点校正量和检波点校正量,可得到静校正的低频静校正量与高频静校正量,将原始数据加上炮点、检波点的高频校正量来获得CMP参考面数据;这样,在步骤102中,基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据,所述水平叠加剖面数据即是以CMP参考面数据为基础的叠加数据通过加上低频校正量得到的剖面数据。这样的数据就对应了水平统一基准面。
而为了提高起伏地表的成像精度,在进行波场延拓成像时,要考虑起伏地表因子,起伏地表因子表示起伏地表的起伏程度,虽然该发明方法利用的是水平统一基准面数据,但由于有了起伏地表的因子项,水平叠加剖面数据中每个网格点的起伏地表因子不同,在一实施例中,基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,包括:
从水平叠加剖面数据中,提取时间采样间隔,延拓步长,CMP间距、波场速度和每个网格点的波场值;
根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像。
在一实施例中,采用如下公式,根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像:
p(i,j+1,k)=(1/(I+(a-b)T))((I-(a+b)T)(p(i,j+1,k+1)+δ(i,j)p(i,j,k)))-p(i,j,k+1)
(1)
其中,p(i,j+1,k)、p(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)分别为水平叠加剖面数据中网格点坐标(i,j+1,k)、(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)的波场值;i为沿测线方向的网格点坐标,j为沿深度方向的网格点坐标,k为沿时间方向的网格点坐标;
a=DTDτV2/(32Dx 2),DT为时间为采样间隔,Dτ为延拓步长,Dx为CMP间距,V为波场速度;
b=1/6,I=(0,1,0),T=(-1,2,-1);
δ(i,j)为网格点坐标(i,j)的起伏地表因子。
在上述实施例中,若没有起伏地表因子,则为水平地表波场延拓过程。
图2为本发明实施例中水平叠加剖面数据中网格点坐标的示意图,x轴表示测线方向,i表示该方向的坐标;t轴表示时间记录波场方向,k表示该方向的坐标;τ轴表示时间表示的深度方向,是延拓得到的波场,j表示该方向的坐标;r7、r6、1、2表示4个网格点,其对应的波场值分别为p(i,j+1,k)、p(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)。
为了说明起伏地表波场延拓的偏移方法,就有必要看看水平地表波场延拓偏移的过程。图3为目前传统使用的水平地表波场延拓偏移过程示意图,虽然本发明专利与传统的方法对数据的要求是一样的,数据对应的面均为水平统一基准面,但实质偏移过程是不同的,目前传统方法对应的波场延拓偏移公式如下:
p(i,j+1,k)=(1/(I+(a-b)T))((I-(a+b)T)(p(i,j+1,k+1)+p(i,j,k)))-p(i,j,k+1)(2)
公式(2)中,地表起伏因子的值为1,所以它的延拓偏移就是从水平统一基准面开始的。图3中左边第一列r地震记录波场,为已知,下边最后一排全取0。则波场延拓偏移过程为:
a、由网格点r7、r6、1算出网格点2(如图3),然后由网格点r6、r5、2算出网格点3,再由网格点r5、r4、3算出网格点4,……直到算出网格点6为止。
b、由网格点2、1、7算出网格点8,然后由网格点3、2、8算出网格点9,再由网格点4、3、9算出网格点10,……直到算出网格点11为止。
c、类似以上过程,算出网格点13、14、15,再算出网格点17、18等等。
对角线上的全体波场值为偏移的结果。
从水平基准面开始波场延拓偏移,就按上述步骤完成,波场值如图3所示,从τ为0开始,波场全部加进来进行延拓偏移。
而起伏地表波场延拓偏移,波场值如图4所示,图4为本发明实施例中起伏地表波场延拓偏移过程示意图,从τ为0开始,到τ小于CMP参考面的值为止,这一段的偏移延拓中,波场全部加进来进行延拓偏移,但延拓结果的波场为0,也就是由于起伏地表因子控制,该段过程参与延拓偏移的波场值保持了原波场值不动,不参与延拓偏移。只有当延拓到τ大于CMP参考面的值时,该道的r原始波场值才真正加进来,如图4中重新加的网格点r7、r6、r5、r4(第四列),此时的波场是原始的真实波场,这样就实现真了正的波场延拓是从起伏地表CMP参考面进行的,而本发明实施例中,基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据,且考虑了起伏地表因子,这样避免了由CMP参考面到统一基准面校正给波场延拓造成的误差,提高偏移成像的精度。它通过j增加,k减少即τ增加,t减少的过程来推算。从地表τ为0和时间t为tmax时间出发,向τ增加和t减少的方向向地下延拓,得到偏移结果。
图5为本发明实施例中起伏地表两层介质的模型,介质分界面有两个拐点,分别为记A、B。A点时空域的理论坐标为CMP为500,时间为2700ms;B点时空域的理论坐标为CMP为750,时间为3450ms。图6为通过目前工业界常用的方法进行波场延拓偏移处理偏移结果,图7为采用本发明实施例提出的方法进行起伏地表波场延拓处理的偏移结果,通过对比,无论从偏移归位的准确性还绕射波收敛程度,本发明方法处理的剖面要好于目前工业界常用的方法。
综上所述,在本发明实施例提出的方法中,采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。在上述过程中,在对叠加剖面数据进行波长延拓成像时,考虑了预设的起伏地表因子,从而可以得到起伏地表的准确的成像结果,提高了成像精度。
本发明实施例还提出一种起伏地表的偏移成像装置,其原理与起伏地表的偏移成像方法类似,这里不再赘述。
图8为本发明实施例中起伏地表的偏移成像装置的示意图,该装置包括:
静校正模块801,用于采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;
水平叠加剖面数据获得模块802,用于基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;
起伏地表因子确定模块803,用于确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;
偏移成像模块804,用于基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。
在一实施例中,静校正模块801具体用于:
采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得炮点校正量和检波点校正量;
根据炮点校正量和检波点校正量,获得CMP参考面数据。
在一实施例中,偏移成像模块803具体用于:
从水平叠加剖面数据中,提取时间采样间隔,延拓步长,CMP间距、波场速度和每个网格点的波场值;
根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像。
在一实施例中,偏移成像模块803具体用于:
采用如下公式,根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像:p(i,j+1,k)=(1/(I+(a-b)T))((I-(a+b)T)(p(i,j+1,k+1)+δ(i,j)p(i,j,k)))-p(i,j,k+1)
其中,p(i,j+1,k)、p(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)分别为水平叠加剖面数据中网格点坐标(i,j+1,k)、(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)的波场值;i为沿测线方向的网格点坐标,j为沿深度方向的网格点坐标,k为沿时间方向的网格点坐标;
a=DTDτV2/(32Dx 2),DT为时间为采样间隔,Dτ为延拓步长,Dx为CMP间距,V为波场速度;
b=1/6,I=(0,1,0),T=(-1,2,-1);
δ(i,j)为网格点坐标(i,j)的起伏地表因子。
综上所述,在本发明实施例提出的装置中,采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。在上述过程中,在对叠加剖面数据进行波长延拓成像时,考虑了预设的起伏地表因子,从而可以得到起伏地表的准确的成像结果,提高了成像精度。
本申请的实施例还提供一种计算机设备,图9为本发明实施例中计算机设备的示意图,该计算机设备能够实现上述实施例中的起伏地表的偏移成像方法中全部步骤,所述电子设备具体包括如下内容:
处理器(processor)901、存储器(memory)902、通信接口(CommunicationsInterface)903和总线904;
其中,所述处理器901、存储器902、通信接口903通过所述总线904完成相互间的通信;所述通信接口903用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输;
所述处理器901用于调用所述存储器902中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的起伏地表的偏移成像方法中的全部步骤。
本申请的实施例还提供一种计算机可读存储介质,能够实现上述实施例中的起伏地表的偏移成像方法中全部步骤,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的起伏地表的偏移成像方法的全部步骤。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种起伏地表的偏移成像方法,其特征在于,包括:
采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;
基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;
确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;
基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。
2.如权利要求1所述的起伏地表的偏移成像方法,其特征在于,采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据,包括:
采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得炮点校正量和检波点校正量;
根据炮点校正量和检波点校正量,获得CMP参考面数据。
3.如权利要求1所述的起伏地表的偏移成像方法,其特征在于,基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,包括:
从水平叠加剖面数据中,提取时间采样间隔,延拓步长,CMP间距、波场速度和每个网格点的波场值;
根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像。
4.如权利要求3所述的起伏地表的偏移成像方法,其特征在于,采用如下公式,根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像:
p(i,j+1,k)=(1/(I+(a-b)T))((I-(a+b)T)(p(i,j+1,k+1)+δ(i,j)p(i,j,k)))-p(i,j,k+1)
其中,p(i,j+1,k)、p(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)分别为水平叠加剖面数据中网格点坐标(i,j+1,k)、(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)的波场值;i为沿测线方向的网格点坐标,j为沿深度方向的网格点坐标,k为沿时间方向的网格点坐标;
a=DTDτV2/(32Dx 2),DT为时间为采样间隔,Dτ为延拓步长,Dx为CMP间距,V为波场速度;
b=1/6,I=(0,1,0),T=(-1,2,-1);
δ(i,j)为网格点坐标(i,j)的起伏地表因子。
5.一种起伏地表的偏移成像装置,其特征在于,包括:
静校正模块,用于采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得CMP参考面数据;
水平叠加剖面数据获得模块,用于基于CMP参考面数据,对起伏地表的地震数据进行处理,获得水平叠加剖面数据;
起伏地表因子确定模块,用于确定水平叠加剖面数据中每个网格点的预设的起伏地表因子;
偏移成像模块,用于基于水平叠加剖面数据中预设的起伏地表因子,由起伏面开始对水平叠加剖面数据进行波长延拓成像,所述起伏地表因子用于表示起伏地表的起伏程度,每个网格点的预设的起伏地表因子。
6.如权利要求5所述的起伏地表的偏移成像装置,其特征在于,静校正模块具体用于:
采用预设的水平统一基准面,对起伏地表的地震数据进行静校正计算,获得炮点校正量和检波点校正量;
根据炮点校正量和检波点校正量,获得CMP参考面数据。
7.如权利要求5所述的起伏地表的偏移成像装置,其特征在于,偏移成像模块具体用于:
从水平叠加剖面数据中,提取时间采样间隔,延拓步长,CMP间距、波场速度和每个网格点的波场值;
根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像。
8.如权利要求7所述的起伏地表的偏移成像装置,其特征在于,偏移成像模块具体用于:
采用如下公式,根据时间采样间隔,延拓步长,CMP间距、波场速度,以及每个网格点的预设的起伏地表因子,对每个网格点的波场值进行波长延拓成像:p(i,j+1,k)=(1/(I+(a-b)T))((I-(a+b)T)(p(i,j+1,k+1)+δ(i,j)p(i,j,k)))-p(i,j,k+1)
其中,p(i,j+1,k)、p(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)分别为水平叠加剖面数据中网格点坐标(i,j+1,k)、(i,j+1,k+1)、p(i,j,k)、p(i,j,k+1)的波场值;i为沿测线方向的网格点坐标,j为沿深度方向的网格点坐标,k为沿时间方向的网格点坐标;
a=DTDτV2/(32Dx 2),DT为时间为采样间隔,Dτ为延拓步长,Dx为CMP间距,V为波场速度;
b=1/6,I=(0,1,0),T=(-1,2,-1);
δ(i,j)为网格点坐标(i,j)的起伏地表因子。
9.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至4任一项所述方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至4任一项所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010703848.6A CN111999767B (zh) | 2020-07-21 | 2020-07-21 | 起伏地表的偏移成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010703848.6A CN111999767B (zh) | 2020-07-21 | 2020-07-21 | 起伏地表的偏移成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111999767A true CN111999767A (zh) | 2020-11-27 |
CN111999767B CN111999767B (zh) | 2023-09-26 |
Family
ID=73466535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010703848.6A Active CN111999767B (zh) | 2020-07-21 | 2020-07-21 | 起伏地表的偏移成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111999767B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5596547A (en) * | 1995-10-20 | 1997-01-21 | Bancroft; John C. | Prestack seismic migration |
US6201765B1 (en) * | 1995-11-27 | 2001-03-13 | Geco-Prakle (Uk) Limited | Method of monitoring quality of seismic data processing and method of processing vertical seismic profile data |
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
CN101315427A (zh) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | 一种复杂地区地震勘探资料处理的方法及系统 |
CN104570069A (zh) * | 2013-10-12 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种从起伏地表直接成像的地震成像方法和装置 |
CN105717538A (zh) * | 2014-12-02 | 2016-06-29 | 中国石油天然气股份有限公司 | 起伏地表地震数据偏移基准面转换方法及装置 |
CN106338760A (zh) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | 误差补偿的起伏地表偏移方法 |
CN109581494A (zh) * | 2018-10-23 | 2019-04-05 | 中国石油天然气集团有限公司 | 叠前偏移方法及装置 |
-
2020
- 2020-07-21 CN CN202010703848.6A patent/CN111999767B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5596547A (en) * | 1995-10-20 | 1997-01-21 | Bancroft; John C. | Prestack seismic migration |
US6201765B1 (en) * | 1995-11-27 | 2001-03-13 | Geco-Prakle (Uk) Limited | Method of monitoring quality of seismic data processing and method of processing vertical seismic profile data |
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
CN101315427A (zh) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | 一种复杂地区地震勘探资料处理的方法及系统 |
CN104570069A (zh) * | 2013-10-12 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种从起伏地表直接成像的地震成像方法和装置 |
CN105717538A (zh) * | 2014-12-02 | 2016-06-29 | 中国石油天然气股份有限公司 | 起伏地表地震数据偏移基准面转换方法及装置 |
CN106338760A (zh) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | 误差补偿的起伏地表偏移方法 |
CN109581494A (zh) * | 2018-10-23 | 2019-04-05 | 中国石油天然气集团有限公司 | 叠前偏移方法及装置 |
Non-Patent Citations (7)
Title |
---|
刘国峰;刘洪;李博;刘钦;佟小龙;: "起伏地表直接叠前时间偏移", 石油地球物理勘探, no. 02 * |
杨海生等: "基于起伏地表的合成平面波叠前深度偏移", 《石油物探》 * |
杨海生等: "基于起伏地表的合成平面波叠前深度偏移", 《石油物探》, no. 02, 25 March 2011 (2011-03-25) * |
熊晓军等: "起伏地表的自由基准面波场成像方法研究", 《石油天然气学报》 * |
熊晓军等: "起伏地表的自由基准面波场成像方法研究", 《石油天然气学报》, no. 03, 15 June 2008 (2008-06-15) * |
谢万学;李德珍;金德刚;何光明;彭文;陈丹;: "起伏地表叠前成像技术在川东高陡构造工区中的应用", 地球物理学进展, no. 05 * |
郑鸿明, 赵峰, 常玉蓉: "近似层替换的基准面静校正方法", 新疆地质, no. 01 * |
Also Published As
Publication number | Publication date |
---|---|
CN111999767B (zh) | 2023-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031898B (zh) | 数据优化方法及积分法叠前深度偏移方法 | |
CN112099082B (zh) | 一种共面元共方位角道集的地震回折波走时反演方法 | |
CN112379436A (zh) | 混合震源激发地震资料静校正量计算方法及装置 | |
CN111999767A (zh) | 起伏地表的偏移成像方法及装置 | |
CN109581494B (zh) | 叠前偏移方法及装置 | |
CN111596360B (zh) | 一种提高三维初至层析反演精度的方法及装置 | |
CN112099090B (zh) | 地震资料视速度域非一致性长波长静校正方法 | |
CN111337973B (zh) | 地震数据重建方法、系统 | |
CN113866827A (zh) | 一种解释性速度建模地震成像方法、系统、介质和设备 | |
CN110967752B (zh) | 一种起伏地表偏移输入地震数据的静校正方法 | |
CN112379429A (zh) | 地震数据的振幅补偿方法及装置 | |
CN111435172A (zh) | 层析静校正的方法和装置 | |
CN111538077A (zh) | 基于倾角约束的叠前深度偏移方法及装置 | |
CN113031071B (zh) | 地震波长波长静校正方法及装置 | |
CN112379420B (zh) | 高精度弯曲测线叠前时间域成像方法及装置 | |
CN112444876A (zh) | 一种地震数据处理方法及装置 | |
CN114063165B (zh) | 一种三维地震数据的拼接方法及装置 | |
CN113093279B (zh) | 转换波静校正方法及装置 | |
CN113568047B (zh) | 一种叠前成像道集生成方法和装置 | |
CN114002742B (zh) | 一种欧拉高斯束偏移成像方法及装置 | |
CN114966852A (zh) | 双界面控制浅深层速度融合与静校正量计算方法及装置 | |
CN113093279A (zh) | 转换波静校正方法及装置 | |
CN112394389A (zh) | 一种二维地震勘探弯线炮检点确定方法及系统 | |
CN116594058A (zh) | 地震数据的全波形反演方法 | |
CN115542387A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |