CN110988993B - 一种偏移成像方法、装置及电子设备 - Google Patents
一种偏移成像方法、装置及电子设备 Download PDFInfo
- Publication number
- CN110988993B CN110988993B CN201911179498.1A CN201911179498A CN110988993B CN 110988993 B CN110988993 B CN 110988993B CN 201911179498 A CN201911179498 A CN 201911179498A CN 110988993 B CN110988993 B CN 110988993B
- Authority
- CN
- China
- Prior art keywords
- imaging
- seismic
- target imaging
- detector
- seismic source
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 495
- 230000001788 irregular Effects 0.000 claims abstract description 133
- 238000000034 method Methods 0.000 claims abstract description 81
- 238000013213 extrapolation Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims description 37
- 238000004590 computer program Methods 0.000 claims description 21
- 230000007246 mechanism Effects 0.000 claims description 4
- 230000005012 migration Effects 0.000 abstract description 8
- 238000013508 migration Methods 0.000 abstract description 8
- 230000006870 function Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 11
- 230000003287 optical effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 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
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010897 surface acoustic wave method Methods 0.000 description 1
- 230000000007 visual effect Effects 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
-
- 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/282—Application of seismic models, synthetic seismograms
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例公开了一种偏移成像方法、装置及电子设备,所述方法包括:基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置,通过间断有限元方法,基于震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,基于震源波场以及检波器波场,对目标成像区域进行偏移成像,以保证成像的准确性。
Description
技术领域
本发明涉及计算机技术领域,尤其涉及一种偏移成像方法、装置及电子设备。
背景技术
地震勘探方法是人类获取地下空间信息的重要手段,其中,地震偏移成像是地震数据处理中的关键步骤。偏移可以使反射波准确归位,具有成像精度高、适应性广的特点,可以直观地展现地下构造的真实形态。
目前,可以使用有限差分方法基于不规则网格剖分,对地形结构复杂、地表起伏较大的成像区域进行偏移成像。
但是,在使用有限差分方法基于不规则网格剖分对成像区域进行偏移成像时,存在以下问题:首先,由于需要引入额外的处理操作(如需选取高精度的数值对波动方程进行求解的处理操作),所以会导致处理复杂度较高、处理效率较低,其次,有限差分方法的数值解只定义在网格中单元的顶点上,所以无法准确的对震源和接收器进行定位,在近地表会产生较强的假象,导致成像不准确。
发明内容
本发明实施例的目的是提供一种偏移成像方法、装置及电子设备,以解决现有技术中在对成像区域进行偏移成像时,存在的数据处理复杂度较高、处理效率较低以及成像不准确的问题。
为解决上述技术问题,本发明实施例是这样实现的:
第一方面,本发明实施例提供的一种偏移成像方法,所述方法包括:
基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定所述震源在不规则网格中的对应位置,以及所述检波器在所述不规则网格中的对应位置;
通过间断有限元方法,基于所述震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场;
在所述震源发射一炮地震波的情况下,获取对应的所述检波器检测到的波形数据;
通过所述间断有限元方法,基于所述检波器检测到的波形数据、所述检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场;
基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像。
可选地,所述震源包括两个及以上震源,所述基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像,包括:
基于所述震源中每个震源对应的所述震源波场,以及与所述每个震源发射的地震波对应的所述检波器波场,确定所述每个震源对应的所述目标成像区域的偏移成像结果;
基于所述每个震源对应的所述目标成像区域的偏移成像结果,对所述目标成像区域进行偏移成像。
可选地,所述目标成像区域中包括成像点,所述基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像,包括:
基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值;
基于所述目标成像区域中的所述成像点的成像值,对所述目标成像区域进行偏移成像。
可选地,所述基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值,包括:
基于预设时间点,获取所述目标成像区域内目标成像点对应的所述震源波场和所述检波器波场,所述目标成像点为所述目标成像区域内的任意一个成像点;
将所述目标成像点对应的所述震源波场和所述检波器波场代入公式
得到所述目标成像点的成像值,其中,x,y分别所述目标成像点在预设水平面内的第一位置坐标和第二位置坐标,z为所述目标成像点的纵向深度坐标,t为所述预设时间点,I(x,y,z)为所述目标成像点的成像值,S(x,y,z,t)为与所述预设时间点对应的所述目标成像点的震源波场,R(x,y,z,t)为与所述预设时间点对应的所述目标成像点的检波器波场;
基于得到所述目标成像点的成像值的方式,确定所述目标成像区域中的其它成像点的成像值。
可选地,所述基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定所述震源在不规则网格中的对应位置,以及所述检波器在所述不规则网格中的对应位置之前,还包括:
基于预设边界确定模型,确定所述目标成像区域的地形构造;
基于所述目标成像区域的地形构造,对所述目标成像区域进行不规则网格剖分,得到所述对目标成像区域进行不规则网格剖分的结果。
第二方面,本发明实施例提供了一种偏移成像装置,所述装置包括:
位置确定模块,用于基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定所述震源在不规则网格中的对应位置,以及所述检波器在所述不规则网格中的对应位置;
第一波场确定模块,用于通过间断有限元方法,基于所述震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场;
数据获取模块,用于在所述震源发射一炮地震波的情况下,获取对应的所述检波器检测到的波形数据;
第二波场确定模块,用于通过所述间断有限元方法,基于所述检波器检测到的波形数据、所述检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场;
成像模块,用于基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像。
可选地,所述震源包括两个及以上震源,所述成像模块,包括:
第一成像单元,用于基于所述震源中每个震源对应的所述震源波场,以及与所述每个震源发射的地震波对应的所述检波器波场,确定所述每个震源对应的所述目标成像区域的偏移成像结果;
第二成像单元,用于基于所述每个震源对应的所述目标成像区域的偏移成像结果,对所述目标成像区域进行偏移成像。
可选地,所述目标成像区域中包括成像点,所述成像模块,包括:
成像值确定单元,用于基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值;
第三成像单元,用于基于所述目标成像区域中的所述成像点的成像值,对所述目标成像区域进行偏移成像。
可选地,所述成像值确定单元,用于:
基于预设时间点,获取所述目标成像区域内目标成像点对应的所述震源波场和所述检波器波场,所述目标成像点为所述目标成像区域内的任意一个成像点;
将所述目标成像点对应的所述震源波场和所述检波器波场代入公式
得到所述目标成像点的成像值,其中,x,y分别所述目标成像点在预设水平面内的第一位置坐标和第二位置坐标,z为所述目标成像点的纵向深度坐标,t为所述预设时间点,I(x,y,z)为所述目标成像点的成像值,S(x,y,z,t)为与所述预设时间点对应的所述目标成像点的震源波场,R(x,y,z,t)为与所述预设时间点对应的所述目标成像点的检波器波场;
基于得到所述目标成像点的成像值的方式,确定所述目标成像区域中的其它成像点的成像值。
可选地,所述装置,还包括:
构造确定模块,用于基于预设边界确定模型,确定所述目标成像区域的地形构造;
不规则剖分模块,用于基于所述目标成像区域的地形构造,对所述目标成像区域进行不规则网格剖分,得到所述对目标成像区域进行不规则网格剖分的结果。
第三方面,本发明实施例提供一种电子设备,包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述计算机程序被所述处理器执行时实现上述第一方面提供的偏移成像方法的步骤。
第四方面,本发明实施例提供一种计算机可读存储介质,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现上述第一方面提供的偏移成像方法的步骤。
由以上本发明实施例提供的技术方案可见,本发明实施例基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置,通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,基于震源波场以及检波器波场,对目标成像区域进行偏移成像。这样,在不需要引入额外处理操作的情况下,提高了数据处理效率,同时,由于间断有限元方法可以使数据解在不规则网格中的任意位置都有定义,所以可以保证震源和检波器的准确定位,提高对目标成像区域的偏移成像的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一种偏移成像方法的流程示意图;
图2为本发明另一种偏移成像方法的流程示意图;
图3为本发明一种对目标成像区域的不规则三角网格剖分的结果的示意图;
图4为本发明一种偏移成像装置的结构示意图;
图5为本发明一种电子设备的结构示意图。
具体实施方式
本发明实施例提供一种偏移成像方法、装置及电子设备。
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
实施例一
如图1所示,本说明书实施例提供一种偏移成像方法,该方法的执行主体可以为服务器,该服务器可以是独立的服务器,也可以是由多个服务器组成的服务器集群。该方法具体可以包括以下步骤:
在S102中,基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置。
其中,目标成像区域可以是野外采集的区域,获取的震源和检波器的位置信息可以为对目标成像区域进行地震波采集时,震源和检波器的实际勘探位置信息,该位置信息可以为通过勘探资料获取的已知数据。
在实施中,地震勘探方法是人类获取地下空间信息的重要手段,其中,地震偏移成像是地震数据处理中的关键步骤。偏移可以使反射波准确归为,绕射波收敛,从而直观地展现地下构造的真实形态。目前,可以使用有限差分方法基于不规则网格剖分,对地形结构复杂、地表起伏较大的成像区域进行偏移成像。
但是,在使用有限差分方法基于不规则网格剖分对成像区域进行偏移成像时,存在以下问题:首先,由于需要引入额外的处理操作(如需选取高精度的数值对波动方程进行求解的处理操作),所以会导致处理复杂度较高、处理效率较低,其次,有限差分方法的数值解只定义在网格中单元的顶点上,所以无法准确的对震源和接收器进行定位,在近地表会产生较强的假象,导致成像不准确。为此,本发明实施例提供一种能够解决上述问题的技术方案,具体可以参见下述内容:
根据实际勘探资料,用户可以获取到多个针对目标成像区域的震源和检波器的位置信息,用户可以将这多个震源和检波器的位置信息输入到服务器中,即服务器可以获取多个针对目标成像区域的震源和检波器的位置信息,如,根据实际勘探资料,针对目标成像区域有278个震源以及480个检波器的位置信息(如位置坐标),用户可以将这278个震源以及480个检波器的位置信息输入到服务器中。
服务器在获取到目标成效区域的不规则网格剖分的结果后,可以将用户输入的这278个震源中每个震源的位置信息与目标成像区域的不规则网格剖分结果进行匹配,以确定每个震源的位置信息在目标成像区域的不规则网格中的对应位置(即确定震源在不规则格网格中的对应位置),同样的,服务器也可以将480个检波器中每个检波器的位置信息与目标成像区域的不规则网格剖分结果进行匹配,以确定每个检波器在目标成像区域的不规则网格中的对应位置。其中,震源(或检波器)在不规则网格中的对应位置可以包括震源(或检波器)的位置信息对应的不规则网格中的网格编号以及对应的网格坐标等。
此外,服务器获取的震源和检波器的位置信息中还可以包括用户根据实际勘探资料输入的震源之间的炮间距、每一个震源的最小偏移距和最大偏移距,以及炮孔径等信息。
在S104中,通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场。
其中,预设速度模型可以通过数学表达式c(x,y,z)来表示,其中,x和y分别为在预设水平面内的第一位置坐标和第二位置坐标,z可以为纵向深度坐标,c(x,y,z)即为由x、y、z确定的位置的波速。
在实施中,在给定的炮点放置震源,服务器可以通过间断有限元方法对震源波场进行不规则网格的正演模拟,以对地下每一点的震源波场进行数值模拟,并得到地下每一点的震源波场。
在S106中,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据。
在实施中,在给定的炮点放置震源,震源可以产生一炮地震波,地震波在被地面反射后,可以到达地面布置的检波器,即在目标成像区域内放置的多个检波器可以接收到由震源发射的地震波的反射波,且检波器可以记录该反射波的波形数据,可以将检波器检测到的波形数据输入到服务器中,即服务器可以获取与震源对应的检波器检测到的波形数据。
在S108中,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场。
在实施中,服务器中可以设置有基于间断有限元方法的处理机制。服务器可以使用上述检波器检测到的波形数据作为边界条件,并通过间断有限元方法,基于检波器在不规则网格中的对应位置和预设速度模型进行逆时外推,得到对应的检波器波场。
此外,上述服务器通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,可以从T1时刻(如第0秒)计算至T2时刻(如第7.996秒),而服务器通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,可以是从T2时刻(即第7.996)秒反向计算至T1时刻(即第0秒)。
在S110中,基于震源波场以及检波器波场,对目标成像区域进行偏移成像。
在实施中,服务器基于得到的震源波场和检波器波场,可以应用震源补充的归一化互相关成像条件,对目标成像区域进行偏移成像。
本发明实施例提供一种偏移成像方法,本发明实施例基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置,通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,基于震源波场以及检波器波场,对目标成像区域进行偏移成像。这样,在不需要引入额外处理操作的情况下,提高了数据处理效率,同时,由于间断有限元方法可以使数据解在不规则网格中的任意位置都有定义,所以可以保证震源和检波器的准确定位,提高对目标成像区域的偏移成像的准确性。
实施例二
如图2所示,本发明实施例提供一种偏移成像方法,该方法的执行主体可以为服务器,该服务器可以是独立的服务器,也可以是由多个服务器组成的服务器集群。该方法具体可以包括以下步骤:
在S202中,基于预设边界确定模型,确定目标成像区域的地形构造。
在实施中,用户可以将目标成像区域的偏移速度值输入到服务器中,以使服务器根据目标成像区域的偏移速度值,确定目标成像区域的地形构造,例如,服务器可以根据地表速度间断面来确定目标成像区域的地表边界位置(即确定目标成像区域的地形构造)。或者,服务器也可以根据实际地表海拔的测量值,确定目标成像区域的地形构造。
用于确定目标成像区域的地形构造的预设边界确定模型可以有多种,上述两种确定目标成像区域的地形构造的预设边界确定模型是一种可选地、可实现的预设边界确定模型,预设边界确定模型可以根据实际应用场景的不同而有所不同,本发明实施例对子不做具体限定。
在S204中,基于目标成像区域的地形构造,对目标成像区域进行不规则网格剖分,得到对目标成像区域进行不规则网格剖分的结果。
在实施中,服务器中可以预先设置有Delaunay三角剖分算法的处理机制,并可以通过Delaunay三角剖分算法对目标成像区域进行不规则三角网格剖分,例如,服务器对目标成像区域进行不规则三角网格剖分后,可以得到如图3所示的,不规则网格剖分的结果。
此外,对目标成像区域进行不规则网格剖分的方法可以是多种多样的,可以根据目标成像区域的地形构造的不同,选用不同的不规则网格剖分的方法,本发明实施例对此不做具体限定。
在S206中,基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置。
在S208中,通过间断有限元方法,基于震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场。
在S210中,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据。
在S212中,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场。
上述S206~S212的具体处理过程可以参见上述实施例一中的S102~S108中的相关内容,在此不再赘述。
在S214中,基于震源中每个震源对应的震源波场,以及与每个震源发射的地震波对应的检波器波场,确定每个震源对应的目标成像区域的偏移成像结果。
在实施中,当震源包括两个及以上震源时,服务器可以获取每个震源对应的震源波场,以及在每个震源发射一炮地震波的情况下,服务器计算得到的对应的检波器波场,在基于每个震源对应的震源波场,以及与每个震源发射的地震波对应的检波器波场,确定每个震源对应的目标成像区域的偏移成像结果。
在实际应用中,上述S214的处理方式可以多种多样,以下再提供一种可选的实现方式,具体可以参见下述步骤一~步骤二。
步骤一,基于震源波场以及检波器波场,确定目标成像区域中的成像点的成像值。
在实施中,可以基于预设时间点,获取目标成像区域内目标成像点对应的震源波场和所述检波器波场,目标成像点为目标成像区域内的任意一个成像点。
然后可以将目标成像点对应的震源波场和所述检波器波场代入公式
得到目标成像点的成像值,其中,x,y分别目标成像点在预设水平面内的第一位置坐标和第二位置坐标,z为目标成像点的纵向深度坐标,t为预设时间点,I(x,y,z)为目标成像点的成像值,S(x,y,z,t)为与预设时间点对应的目标成像点的震源波场,R(x,y,z,t)为与预设时间点对应的目标成像点的检波器波场。
基于得到目标成像点的成像值的方式,确定目标成像区域中的其它成像点的成像值。
步骤二,基于目标成像区域中的成像点的成像值,对目标成像区域进行偏移成像。
在实施中,服务器基于目标成像区域中的成像点的成像值,可以得到每个震源对应的由目标成像区域内的成像点的成像值构成的,对目标成像区域的偏移成像结果(即对目标成像区域进行偏移成像)。
在S216中,基于每个震源对应的目标成像区域的偏移成像结果,对目标成像区域进行偏移成像。
在实施中,服务器可以将每个震源对应的目标成像区域的偏移成像结果进行叠加,以对目标成像区域进行偏移成像。例如,服务器可以将每个震源对应的目标成像区域中每个成像点的成像值进行叠加,以得到目标成像区域内所有成像点的最终的成像值,即对目标成像区域进行偏移成像。
本发明实施例提供一种偏移成像方法,本发明实施例基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置,通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,基于震源波场以及检波器波场,对目标成像区域进行偏移成像。这样,在不需要引入额外处理操作的情况下,提高了数据处理效率,同时,由于间断有限元方法可以使数据解在不规则网格中的任意位置都有定义,所以可以保证震源和检波器的准确定位,提高对目标成像区域的偏移成像的准确性。
实施例三
以上为本说明书实施例提供的偏移成像方法,基于同样的思路,本说明书实施例还提供一种偏移成像装置,如图4所示。
该偏移成像装置包括:位置确定模块401、第一波场确定模块402、数据获取模块403、第二波场确定模块404和成像模块405,中:
位置确定模块401,用于基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定所述震源在不规则网格中的对应位置,以及所述检波器在所述不规则网格中的对应位置;
第一波场确定模块402,用于通过间断有限元方法,基于所述震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场;
数据获取模块403,用于在所述震源发射一炮地震波的情况下,获取对应的所述检波器检测到的波形数据;
第二波场确定模块404,用于通过所述间断有限元方法,基于所述检波器检测到的波形数据、所述检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场;
成像模块405,用于基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像。
本发明实施例中,所述震源包括两个及以上震源,所述成像模块405,包括:
第一成像单元,用于基于所述震源中每个震源对应的所述震源波场,以及与所述每个震源发射的地震波对应的所述检波器波场,确定所述每个震源对应的所述目标成像区域的偏移成像结果;
第二成像单元,用于基于所述每个震源对应的所述目标成像区域的偏移成像结果,对所述目标成像区域进行偏移成像。
本发明实施例中,所述目标成像区域中包括成像点,所述成像模块405,包括:
成像值确定单元,用于基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值;
第三成像单元,用于基于所述目标成像区域中的所述成像点的成像值,对所述目标成像区域进行偏移成像。
本发明实施例中,所述成像值确定单元,用于:
基于预设时间点,获取所述目标成像区域内目标成像点对应的所述震源波场和所述检波器波场,所述目标成像点为所述目标成像区域内的任意一个成像点;
将所述目标成像点对应的所述震源波场和所述检波器波场代入公式
得到所述目标成像点的成像值,其中,x,y分别所述目标成像点在预设水平面内的第一位置坐标和第二位置坐标,z为所述目标成像点的纵向深度坐标,t为所述预设时间点,I(x,y,z)为所述目标成像点的成像值,S(x,y,z,t)为与所述预设时间点对应的所述目标成像点的震源波场,R(x,y,z,t)为与所述预设时间点对应的所述目标成像点的检波器波场;
基于得到所述目标成像点的成像值的方式,确定所述目标成像区域中的其它成像点的成像值。
本发明实施例中,所述装置,还包括:
构造确定模块,用于基于预设边界确定模型,确定所述目标成像区域的地形构造;
不规则剖分模块,用于基于所述目标成像区域的地形构造,对所述目标成像区域进行不规则网格剖分,得到所述对目标成像区域进行不规则网格剖分的结果。
本发明实施例提供一种偏移成像方法,本发明实施例基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置,通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,基于震源波场以及检波器波场,对目标成像区域进行偏移成像。这样,在不需要引入额外处理操作的情况下,提高了数据处理效率,同时,由于间断有限元方法可以使数据解在不规则网格中的任意位置都有定义,所以可以保证震源和检波器的准确定位,提高对目标成像区域的偏移成像的准确性。
实施例四
图5为实现本发明各个实施例的一种电子设备的硬件结构示意图,
该电子设备500包括但不限于:射频单元501、网络模块502、音频输出单元503、输入单元504、传感器505、显示单元506、用户输入单元507、接口单元508、存储器509、处理器510、以及电源511等部件。本领域技术人员可以理解,图5中示出的电子设备结构并不构成对电子设备的限定,电子设备可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。在本发明实施例中,电子设备包括但不限于手机、平板电脑、笔记本电脑、掌上电脑、车载终端、可穿戴设备、以及计步器等。
其中,处理器510,用于基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定所述震源在不规则网格中的对应位置,以及所述检波器在所述不规则网格中的对应位置;
所述处理器510,还用于通过间断有限元方法,基于所述震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场;
所述处理器510,还用于在所述震源发射一炮地震波的情况下,获取对应的所述检波器检测到的波形数据;
所述处理器510,还用于通过所述间断有限元方法,基于所述检波器检测到的波形数据、所述检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场;
所述处理器510,还用于基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像。
此外,所述处理器510,还用于基于所述震源中每个震源对应的所述震源波场,以及与所述每个震源发射的地震波对应的所述检波器波场,确定所述每个震源对应的所述目标成像区域的偏移成像结果;
另外,所述处理器510,还用于基于所述每个震源对应的所述目标成像区域的偏移成像结果,对所述目标成像区域进行偏移成像。
此外,所述处理器510,还用于基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值;
另外,所述处理器510,还用于基于所述目标成像区域中的所述成像点的成像值,对所述目标成像区域进行偏移成像。
此外,所述处理器510,还用于基于预设时间点,获取所述目标成像区域内目标成像点对应的所述震源波场和所述检波器波场,所述目标成像点为所述目标成像区域内的任意一个成像点;
另外,所述处理器510,还用于将所述目标成像点对应的所述震源波场和所述检波器波场代入公式
得到所述目标成像点的成像值,其中,x,y分别所述目标成像点在预设水平面内的第一位置坐标和第二位置坐标,z为所述目标成像点的纵向深度坐标,t为所述预设时间点,I(x,y,z)为所述目标成像点的成像值,S(x,y,z,t)为与所述预设时间点对应的所述目标成像点的震源波场,R(x,y,z,t)为与所述预设时间点对应的所述目标成像点的检波器波场;
此外,所述处理器510,还用于基于得到所述目标成像点的成像值的方式,确定所述目标成像区域中的其它成像点的成像值。
另外,所述处理器510,还用于基于预设边界确定模型,确定所述目标成像区域的地形构造;
此外,所述处理器510,还用于基于所述目标成像区域的地形构造,对所述目标成像区域进行不规则网格剖分,得到所述对目标成像区域进行不规则网格剖分的结果。
本发明实施例提供一种电子设备,本发明实施例基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置,通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,基于震源波场以及检波器波场,对目标成像区域进行偏移成像。这样,在不需要引入额外处理操作的情况下,提高了数据处理效率,同时,由于间断有限元方法可以使数据解在不规则网格中的任意位置都有定义,所以可以保证震源和检波器的准确定位,提高对目标成像区域的偏移成像的准确性。
应理解的是,本发明实施例中,射频单元501可用于收发信息或通话过程中,信号的接收和发送,具体的,将来自基站的下行数据接收后,给处理器510处理;另外,将上行的数据发送给基站。通常,射频单元501包括但不限于天线、至少一个放大器、收发信机、耦合器、低噪声放大器、双工器等。此外,射频单元501还可以通过无线通信系统与网络和其他设备通信。
电子设备通过网络模块502为用户提供了无线的宽带互联网访问,如帮助用户收发电子邮件、浏览网页和访问流式媒体等。
音频输出单元503可以将射频单元501或网络模块502接收的或者在存储器509中存储的音频数据转换成音频信号并且输出为声音。而且,音频输出单元503还可以提供与电子设备500执行的特定功能相关的音频输出(例如,呼叫信号接收声音、消息接收声音等等)。音频输出单元503包括扬声器、蜂鸣器以及受话器等。
输入单元504用于接收音频或视频信号。输入单元504可以包括图形处理器(Graphics Processing Unit,GPU)5041和麦克风5042,图形处理器5041对在视频捕获模式或图像捕获模式中由图像捕获装置(如摄像头)获得的静态图片或视频的图像数据进行处理。处理后的图像帧可以显示在显示单元506上。经图形处理器5041处理后的图像帧可以存储在存储器509(或其它存储介质)中或者经由射频单元501或网络模块502进行发送。麦克风5042可以接收声音,并且能够将这样的声音处理为音频数据。处理后的音频数据可以在电话通话模式的情况下转换为可经由射频单元501发送到移动通信基站的格式输出。
电子设备500还包括至少一种传感器505,比如光传感器、运动传感器以及其他传感器。具体地,光传感器包括环境光传感器及接近传感器,其中,环境光传感器可根据环境光线的明暗来调节显示面板5061的亮度,接近传感器可在电子设备500移动到耳边时,关闭显示面板5061和/或背光。作为运动传感器的一种,加速计传感器可检测各个方向上(一般为三轴)加速度的大小,静止时可检测出重力的大小及方向,可用于识别电子设备姿态(比如横竖屏切换、相关游戏、磁力计姿态校准)、振动识别相关功能(比如计步器、敲击)等;传感器505还可以包括指纹传感器、压力传感器、虹膜传感器、分子传感器、陀螺仪、气压计、湿度计、温度计、红外线传感器等,在此不再赘述。
显示单元506用于显示由用户输入的信息或提供给用户的信息。显示单元506可包括显示面板5061,可以采用液晶显示器(Liquid Crystal Display,LCD)、有机发光二极管(Organic Light-Emitting Diode,OLED)等形式来配置显示面板5061。
用户输入单元507可用于接收输入的数字或字符信息,以及产生与电子设备的用户设置以及功能控制有关的键信号输入。具体地,用户输入单元507包括触控面板5051以及其他输入设备5072。触控面板5051,也称为触摸屏,可收集用户在其上或附近的触摸操作(比如用户使用手指、触笔等任何适合的物体或附件在触控面板5051上或在触控面板5051附近的操作)。触控面板5051可包括触摸检测装置和触摸控制器两个部分。其中,触摸检测装置检测用户的触摸方位,并检测触摸操作带来的信号,将信号传送给触摸控制器;触摸控制器从触摸检测装置上接收触摸信息,并将它转换成触点坐标,再送给处理器510,接收处理器510发来的命令并加以执行。此外,可以采用电阻式、电容式、红外线以及表面声波等多种类型实现触控面板5051。除了触控面板5051,用户输入单元507还可以包括其他输入设备5072。具体地,其他输入设备5072可以包括但不限于物理键盘、功能键(比如音量控制按键、开关按键等)、轨迹球、鼠标、操作杆,在此不再赘述。
进一步的,触控面板5051可覆盖在显示面板5061上,当触控面板5051检测到在其上或附近的触摸操作后,传送给处理器510以确定触摸事件的类型,随后处理器510根据触摸事件的类型在显示面板5061上提供相应的视觉输出。虽然在图5中,触控面板5051与显示面板5061是作为两个独立的部件来实现电子设备的输入和输出功能,但是在某些实施例中,可以将触控面板5051与显示面板5061集成而实现电子设备的输入和输出功能,具体此处不做限定。
接口单元508为外部装置与电子设备500连接的接口。例如,外部装置可以包括有线或无线头戴式耳机端口、外部电源(或电池充电器)端口、有线或无线数据端口、存储卡端口、用于连接具有识别模块的装置的端口、音频输入/输出(I/O)端口、视频I/O端口、耳机端口等等。接口单元508可以用于接收来自外部装置的输入(例如,数据信息、电力等等)并且将接收到的输入传输到电子设备500内的一个或多个元件或者可以用于在电子设备500和外部装置之间传输数据。
存储器509可用于存储软件程序以及各种数据。存储器509可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等;存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外,存储器409可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
处理器510是电子设备的控制中心,利用各种接口和线路连接整个电子设备的各个部分,通过运行或执行存储在存储器509内的软件程序和/或模块,以及调用存储在存储器509内的数据,执行电子设备的各种功能和处理数据,从而对电子设备进行整体监控。处理器510可包括一个或多个处理单元;优选的,处理器510可集成应用处理器和调制解调处理器,其中,应用处理器主要处理操作系统、用户界面和应用程序等,调制解调处理器主要处理无线通信。可以理解的是,上述调制解调处理器也可以不集成到处理器510中。
电子设备500还可以包括给各个部件供电的电源511(比如电池),优选的,电源511可以通过电源管理系统与处理器510逻辑相连,从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。
优选的,本发明实施例还提供一种电子设备,包括处理器510,存储器509,存储在存储器509上并可在所述处理器510上运行的计算机程序,该计算机程序被处理器510执行时实现上述偏移成像方法实施例的各个过程,且能达到相同的技术效果,为避免重复,这里不再赘述。
实施例五
本发明实施例还提供一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述偏移成像方法实施例的各个过程,且能达到相同的技术效果,为避免重复,这里不再赘述。其中,所述的计算机可读存储介质,如只读存储器(Read-OnlyMemory,简称ROM)、随机存取存储器(RandomAccess Memory,简称RAM)、磁碟或者光盘等。
本发明实施例提供一种计算机可读存储介质,基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定震源在不规则网格中的对应位置,以及检波器在不规则网格中的对应位置,通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,在震源发射一炮地震波的情况下,获取对应的检波器检测到的波形数据,通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,基于震源波场以及检波器波场,对目标成像区域进行偏移成像。这样,在不需要引入额外处理操作的情况下,提高了数据处理效率,同时,由于间断有限元方法可以使数据解在不规则网格中的任意位置都有定义,所以可以保证震源和检波器的准确定位,提高对目标成像区域的偏移成像的准确性。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitorymedia),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本发明的实施例可提供为方法、系统或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上所述仅为本发明的实施例而已,并不用于限制本发明。对于本领域技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (8)
1.一种偏移成像方法,其特征在于,包括:
基于预设边界确定模型,确定目标成像区域的地形构造;用户将目标成像区域的偏移速度值输入到服务器中,以使服务器根据目标成像区域的偏移速度值,确定目标成像区域的地形构造;
基于所述目标成像区域的地形构造,对所述目标成像区域进行不规则网格剖分,得到对目标成像区域进行不规则网格剖分的结果;服务器中预先设置有Delaunay三角剖分算法的处理机制,并通过Delaunay三角剖分算法对目标成像区域进行不规则三角网格剖分;
基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定所述震源在不规则网格中的对应位置,以及所述检波器在所述不规则网格中的对应位置;目标成像区域是野外采集的区域,获取的震源和检波器的位置信息为对目标成像区域进行地震波采集时,震源和检波器的实际勘探位置信息;根据实际勘探资料,用户获取到多个针对目标成像区域的震源和检波器的位置信息,用户将这多个震源和检波器的位置信息输入到服务器中,即服务器获取多个针对目标成像区域的震源和检波器的位置信息;震源或检波器在不规则网格中的对应位置包括震源或检波器的位置信息对应的不规则网格中的网格编号以及对应的网格坐标;服务器获取的震源和检波器的位置信息中还包括用户根据实际勘探资料输入的震源之间的炮间距、每一个震源的最小偏移距和最大偏移距,以及炮孔径;
通过间断有限元方法,基于所述震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场;
在所述震源发射一炮地震波的情况下,获取对应的所述检波器检测到的波形数据;在给定的炮点放置震源,震源产生一炮地震波,地震波在被地面反射后,到达地面布置的检波器,即在目标成像区域内放置的多个检波器接收到由震源发射的地震波的反射波,且检波器记录该反射波的波形数据,将检波器检测到的波形数据输入到服务器中,即服务器获取与震源对应的检波器检测到的波形数据;
通过所述间断有限元方法,基于所述检波器检测到的波形数据、所述检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场;服务器通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,从T1时刻计算至T2时刻,而服务器通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,是从T2时刻反向计算至T1时刻;
基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像;
所述震源包括两个及以上震源,所述基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像,包括:
基于所述震源中每个震源对应的所述震源波场,以及与所述每个震源发射的地震波对应的所述检波器波场,确定所述每个震源对应的所述目标成像区域的偏移成像结果;当震源包括两个及以上震源时,服务器获取每个震源对应的震源波场,以及在每个震源发射一炮地震波的情况下,服务器计算得到的对应的检波器波场,在基于每个震源对应的震源波场,以及与每个震源发射的地震波对应的检波器波场,确定每个震源对应的目标成像区域的偏移成像结果;
基于所述每个震源对应的所述目标成像区域的偏移成像结果,对所述目标成像区域进行偏移成像;服务器将每个震源对应的目标成像区域中每个成像点的成像值进行叠加,以得到目标成像区域内所有成像点的最终的成像值,即对目标成像区域进行偏移成像。
2.根据权利要求1所述的方法,其特征在于,所述目标成像区域中包括成像点,所述基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像,包括:
基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值;
基于所述目标成像区域中的所述成像点的成像值,对所述目标成像区域进行偏移成像。
3.根据权利要求2所述的方法,其特征在于,所述基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值,包括:
基于预设时间点,获取所述目标成像区域内目标成像点对应的所述震源波场和所述检波器波场,所述目标成像点为所述目标成像区域内的任意一个成像点;
将所述目标成像点对应的所述震源波场和所述检波器波场代入公式
得到所述目标成像点的成像值,其中,x,y分别为所述目标成像点在预设水平面内的第一位置坐标和第二位置坐标,z为所述目标成像点的纵向深度坐标,t为所述预设时间点,I(x,y,z)为所述目标成像点的成像值,S(x,y,z,t)为与所述预设时间点对应的所述目标成像点的震源波场,R(x,y,z,t)为与所述预设时间点对应的所述目标成像点的检波器波场;
基于得到所述目标成像点的成像值的方式,确定所述目标成像区域中的其它成像点的成像值。
4.一种偏移成像装置,其特征在于,包括:
构造确定模块,用于基于预设边界确定模型,确定目标成像区域的地形构造;用户将目标成像区域的偏移速度值输入到服务器中,以使服务器根据目标成像区域的偏移速度值,确定目标成像区域的地形构造;
不规则剖分模块,用于基于所述目标成像区域的地形构造,对所述目标成像区域进行不规则网格剖分,得到对目标成像区域进行不规则网格剖分的结果;服务器中预先设置有Delaunay三角剖分算法的处理机制,并通过Delaunay三角剖分算法对目标成像区域进行不规则三角网格剖分;
位置确定模块,用于基于对目标成像区域进行不规则网格剖分的结果以及获取的震源和检波器的位置信息,确定所述震源在不规则网格中的对应位置,以及所述检波器在所述不规则网格中的对应位置;目标成像区域是野外采集的区域,获取的震源和检波器的位置信息为对目标成像区域进行地震波采集时,震源和检波器的实际勘探位置信息;根据实际勘探资料,用户获取到多个针对目标成像区域的震源和检波器的位置信息,用户将这多个震源和检波器的位置信息输入到服务器中,即服务器获取多个针对目标成像区域的震源和检波器的位置信息;震源或检波器在不规则网格中的对应位置包括震源或检波器的位置信息对应的不规则网格中的网格编号以及对应的网格坐标;服务器获取的震源和检波器的位置信息中还包括用户根据实际勘探资料输入的震源之间的炮间距、每一个震源的最小偏移距和最大偏移距,以及炮孔径;
第一波场确定模块,用于通过间断有限元方法,基于所述震源在所述不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场;
数据获取模块,用于在所述震源发射一炮地震波的情况下,获取对应的所述检波器检测到的波形数据;在给定的炮点放置震源,震源产生一炮地震波,地震波在被地面反射后,到达地面布置的检波器,即在目标成像区域内放置的多个检波器接收到由震源发射的地震波的反射波,且检波器记录该反射波的波形数据,将检波器检测到的波形数据输入到服务器中,即服务器获取与震源对应的检波器检测到的波形数据;
第二波场确定模块,用于通过所述间断有限元方法,基于所述检波器检测到的波形数据、所述检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场;服务器通过间断有限元方法,基于震源在不规则网格中的对应位置以及预设速度模型进行正演模拟,得到对应的震源波场,从T1时刻计算至T2时刻,而服务器通过间断有限元方法,基于检波器检测到的波形数据、检波器在不规则网格中的对应位置以及预设速度模型进行逆时外推,得到对应的检波器波场,是从T2时刻反向计算至T1时刻;
成像模块,用于基于所述震源波场以及所述检波器波场,对所述目标成像区域进行偏移成像;所述震源包括两个及以上震源,所述成像模块,包括:
第一成像单元,用于基于所述震源中每个震源对应的所述震源波场,以及与所述每个震源发射的地震波对应的所述检波器波场,确定所述每个震源对应的所述目标成像区域的偏移成像结果;当震源包括两个及以上震源时,服务器获取每个震源对应的震源波场,以及在每个震源发射一炮地震波的情况下,服务器计算得到的对应的检波器波场,在基于每个震源对应的震源波场,以及与每个震源发射的地震波对应的检波器波场,确定每个震源对应的目标成像区域的偏移成像结果;
第二成像单元,用于基于所述每个震源对应的所述目标成像区域的偏移成像结果,对所述目标成像区域进行偏移成像;服务器将每个震源对应的目标成像区域中每个成像点的成像值进行叠加,以得到目标成像区域内所有成像点的最终的成像值,即对目标成像区域进行偏移成像。
5.根据权利要求4所述的装置,其特征在于,所述目标成像区域中包括成像点,所述成像模块,包括:
成像值确定单元,用于基于所述震源波场以及所述检波器波场,确定所述目标成像区域中的所述成像点的成像值;
第三成像单元,用于基于所述目标成像区域中的所述成像点的成像值,对所述目标成像区域进行偏移成像。
6.根据权利要求5所述的装置,其特征在于,所述成像值确定单元,用于:
基于预设时间点,获取所述目标成像区域内目标成像点对应的所述震源波场和所述检波器波场,所述目标成像点为所述目标成像区域内的任意一个成像点;
将所述目标成像点对应的所述震源波场和所述检波器波场代入公式
得到所述目标成像点的成像值,其中,x,y分别所述目标成像点在预设水平面内的第一位置坐标和第二位置坐标,z为所述目标成像点的纵向深度坐标,t为所述预设时间点,I(x,y,z)为所述目标成像点的成像值,S(x,y,z,t)为与所述预设时间点对应的所述目标成像点的震源波场,R(x,y,z,t)为与所述预设时间点对应的所述目标成像点的检波器波场;
基于得到所述目标成像点的成像值的方式,确定所述目标成像区域中的其它成像点的成像值。
7.一种电子设备,其特征在于,包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述计算机程序被所述处理器执行时实现如权利要求1至3中任一项所述的偏移成像方法的步骤。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现如权利要求1至3中任一项所述的偏移成像方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911179498.1A CN110988993B (zh) | 2019-11-27 | 2019-11-27 | 一种偏移成像方法、装置及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911179498.1A CN110988993B (zh) | 2019-11-27 | 2019-11-27 | 一种偏移成像方法、装置及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110988993A CN110988993A (zh) | 2020-04-10 |
CN110988993B true CN110988993B (zh) | 2021-01-26 |
Family
ID=70087236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911179498.1A Active CN110988993B (zh) | 2019-11-27 | 2019-11-27 | 一种偏移成像方法、装置及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110988993B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113341456B (zh) * | 2021-07-08 | 2024-02-02 | 清华大学 | 一种基于非规则网格的地震偏移方法、装置及电子设备 |
CN113536638B (zh) * | 2021-07-22 | 2023-09-22 | 北京大学 | 基于间断有限元和交错网格的高精度地震波场模拟方法 |
CN114357942B (zh) * | 2022-03-17 | 2022-06-10 | 南京邮电大学 | 基于间断有限元法的集成电路互连线寄生电容提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9279314B2 (en) * | 2011-08-11 | 2016-03-08 | Conocophillips Company | Heat front capture in thermal recovery simulations of hydrocarbon reservoirs |
CN105388520A (zh) * | 2015-10-22 | 2016-03-09 | 中国石油化工股份有限公司 | 一种地震资料叠前逆时偏移成像方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8289809B2 (en) * | 2008-09-08 | 2012-10-16 | Exxonmobil Upstream Research Company | Common reflection azimuth migration |
US20100118654A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling migration method |
CN102590858B (zh) * | 2011-12-31 | 2014-10-08 | 中国石油集团西北地质研究所 | 基于宽频子波重构的双程波成像方法 |
US9390555B2 (en) * | 2012-11-09 | 2016-07-12 | International Business Machines Corporation | Method to assess the impact of existing fractures and faults for reservoir management |
EP3158367A1 (en) * | 2014-06-17 | 2017-04-26 | Exxonmobil Upstream Research Company | Fast viscoacoustic and viscoelastic full-wavefield inversion |
CN105717539B (zh) * | 2016-01-28 | 2018-01-30 | 中国地质大学(北京) | 一种基于多gpu计算的三维tti介质逆时偏移成像方法 |
CN107798156B (zh) * | 2016-09-02 | 2020-12-11 | 赵建国 | 一种频率域2.5维粘弹性波数值模拟方法及装置 |
CN108072899B (zh) * | 2016-11-10 | 2020-09-15 | 中国石油化工股份有限公司 | 间断Galerkin有限元地震数值模拟算法的自适应实现方法 |
CN108073732A (zh) * | 2016-11-10 | 2018-05-25 | 中国石油化工股份有限公司 | 获得稳定的近最佳匹配层吸收边界条件的方法 |
CN108345032A (zh) * | 2018-01-12 | 2018-07-31 | 中国科学技术大学 | 一种弱照明区域高信噪比偏移成像方法 |
-
2019
- 2019-11-27 CN CN201911179498.1A patent/CN110988993B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9279314B2 (en) * | 2011-08-11 | 2016-03-08 | Conocophillips Company | Heat front capture in thermal recovery simulations of hydrocarbon reservoirs |
CN105388520A (zh) * | 2015-10-22 | 2016-03-09 | 中国石油化工股份有限公司 | 一种地震资料叠前逆时偏移成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110988993A (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110988993B (zh) | 一种偏移成像方法、装置及电子设备 | |
CN107576986B (zh) | 一种确定海底地震反射波旅行时的方法和装置 | |
US9369801B2 (en) | Wireless speaker system with noise cancelation | |
JP2015510113A (ja) | 実際の距離を推定するための較正されたハードウェアセンサ | |
US10963100B2 (en) | Interactive object control method and apparatus, terminal, and computer-readable storage medium | |
CN111986250A (zh) | 物体体积测量方法、装置、测量设备及存储介质 | |
CN109696706B (zh) | 盐构造检测方法、装置及存储介质 | |
CN108828613B (zh) | 一种去除噪点的方法和危化品仓储激光扫描装置 | |
CN107505653A (zh) | 一种确定叠前偏移时间结果的方法和装置 | |
CN112147625A (zh) | 一种标定方法、装置、单目激光测量设备及标定系统 | |
CN107782250A (zh) | 一种深度信息测量方法、装置和移动终端 | |
CN111739093B (zh) | 一种蒸汽管道位移的监测方法、装置及电子设备 | |
CN111964773B (zh) | 一种基于多个交叠线性调频组的振动信号提取方法及装置 | |
CN113341456B (zh) | 一种基于非规则网格的地震偏移方法、装置及电子设备 | |
CN112200130B (zh) | 一种三维目标检测方法、装置及终端设备 | |
CN115265529A (zh) | 目标对象的定位方法、系统、装置、电子设备及存储介质 | |
CN113630712B (zh) | 一种定位方法、装置及设备 | |
CN113009573A (zh) | 地震数据采集方法、装置、计算机设备及存储介质 | |
CN111986251A (zh) | 物体体积测量方法、装置、测量设备及存储介质 | |
CN111897001B (zh) | 波场确定方法与装置 | |
CN111551987A (zh) | 利用变偏移距vsp数据确定岩石弹性性质的方法与装置 | |
CN109375232B (zh) | 一种距离测量方法及装置 | |
CN111238495A (zh) | 一种定位车辆的方法和终端设备 | |
CN111505717A (zh) | 一种三维vsp共面元成像方法与装置 | |
CN110058306B (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 |