CN111948706A - 正交各向异性介质地震成像方法及装置 - Google Patents
正交各向异性介质地震成像方法及装置 Download PDFInfo
- Publication number
- CN111948706A CN111948706A CN201910408859.9A CN201910408859A CN111948706A CN 111948706 A CN111948706 A CN 111948706A CN 201910408859 A CN201910408859 A CN 201910408859A CN 111948706 A CN111948706 A CN 111948706A
- Authority
- CN
- China
- Prior art keywords
- medium
- orthotropic
- ray tracing
- seismic
- seismic imaging
- 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 100
- 238000000034 method Methods 0.000 claims abstract description 76
- 230000005012 migration Effects 0.000 claims abstract description 51
- 238000013508 migration Methods 0.000 claims abstract description 51
- 238000012545 processing Methods 0.000 claims abstract description 28
- 230000006870 function Effects 0.000 claims description 33
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 9
- 230000008569 process Effects 0.000 abstract description 8
- 238000011161 development Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 8
- 230000009466 transformation Effects 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005755 formation reaction Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000000750 progressive effect Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- RKTYLMNFRDHKIL-UHFFFAOYSA-N copper;5,10,15,20-tetraphenylporphyrin-22,24-diide Chemical compound [Cu+2].C1=CC(C(=C2C=CC([N-]2)=C(C=2C=CC=CC=2)C=2C=CC(N=2)=C(C=2C=CC=CC=2)C2=CC=C3[N-]2)C=2C=CC=CC=2)=NC1=C3C1=CC=CC=C1 RKTYLMNFRDHKIL-UHFFFAOYSA-N 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 239000002023 wood Substances 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
-
- 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
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
技术领域
本申请涉及地震成像技术领域,具体涉及一种正交各向异性介质地震成像方法及装置。
背景技术
与直接数值求解波动方程不同(有限差分、有限元),射线理论是一类波动理论在高频渐进近似下的射线级数方法。射线理论不仅可以求取射线路径、旅行时,还可以计算近似的地震波振幅以及地震记录,其要求地下的介质参数相对平滑且可存在起伏地层界面。射线理论可分为运动学和动力学两类,前者主要计算射线路径、波前以及旅行时,后者主要估算波场的振幅、位移矢量和地震记录。基于射线追踪的积分法叠前深度偏移相对与波动方程偏移具有计算效率高、成像频带宽、输入输出灵活、硬件门槛低等优点,一直是实际生产中的主力方法。积分法偏移只是利用地震资料中的可信度相对较高的相位信息,存在波动方程偏移无法比拟的稳定输出。实际生产中的叠前深度偏移项目都需要做积分法叠前深度偏移,尤其对宽方位OVT数据,由于其不再是炮集数据,无法采用波动方程方法进行炮集数据偏移,只能采用积分法完成叠前深度偏移。
地下介质的各向异性普遍存在。实际地球介质引起地震各向异性的因素很多、成因很复杂。许多地球物理学家和地震学家通过对地震波在地球介质中的传播现象进行观测,对地震波在各向异性介质中的传播规律和形成机理方面做了大量的研究工作,认识到地球介质存在各向异性;发现地壳上的大多数沉积岩展现出地震各向异性特征。各向异性介质中的地震模拟和实际地震资料中的处理技术有了一定进步,尽管各向异性介质相对复杂,但在地震成像效果方面有了较大的提高。不同实际地下介质模型需要不同的各向异性介质进行近似刻画。最简单的各向异性介质是横向各向同性介质(TI),该介质仅有单一的对称界面,当界面为水平时,则为横向各向同性介质(VTI),主要表示为沉积地层环境下介质。当界面为垂直时,为垂向各向同性介质(HTI),主要表征为垂向发育的裂缝介质。当地层出现倾角时,则对称界面与水平面存在一定夹角,为倾斜横向各向同性(TTI)介质。在构造运动或外力的作用下,由于地层的脆性会产生不同尺度的裂缝,进而成为地下流体存储和运移的主要途径。对于含有裂缝的沉积地层则无法用横向各向同性(TI)介质来描述,需要正交各向异性(ORT)介质,但当前尚未有一种能够同时保证正交各向异性介质的射线追踪的稳定性和准确性的地震成像方式。
因此,亟需设计一种正交各向异性介质地震成像方式,以实现对裂缝型储层及勘探开发潜力巨大的非常规储层的更为准确的刻画。
发明内容
针对现有技术中的问题,本申请提供一种正交各向异性介质地震成像方法及装置,能够有效实现正交各向异性介质射线追踪,追踪过程高效、可靠且准确性高,进而能够更加准确且稳定的实现针对裂缝型储层及勘探开发潜力巨大的非常规储层的定量刻画。
为解决上述技术问题,本申请提供以下技术方案:
第一方面,本申请提供一种正交各向异性介质地震成像方法,包括:
根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数;
基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式;
应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
进一步地,所述应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,包括:
应用龙格库塔方法和射线的初始值求解所述运动学射线追踪系统和动力学射线追踪系统,得到所述正交各向异性介质对应的地震波旅行时场信息。
进一步地,所述根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像,包括:
根据所述正交各向异性介质对应的地震波旅行时场信息,基于积分法叠前深度偏移方法对所述正交各向异性介质进行偏移处理,得到所述正交各向异性介质的地震成像。
第二方面,本申请提供一种正交各向异性介质地震成像装置,包括:
汉密尔顿函数确定模块,用于根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数;
射线追踪模块,用于基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式;
地震成像模块,用于应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
进一步地,所述地震成像模块包括:
地震波旅行时场信息获取单元,用于应用龙格库塔方法和射线的初始值求解所述运动学射线追踪系统和动力学射线追踪系统,得到所述正交各向异性介质对应的地震波旅行时场信息。
进一步地,所述地震成像模块包括:
偏移处理单元,用于根据所述正交各向异性介质对应的地震波旅行时场信息,基于积分法叠前深度偏移方法对所述正交各向异性介质进行偏移处理,得到所述正交各向异性介质的地震成像。
第三方面,本申请提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现所述正交各向异性介质地震成像方法的步骤。
第四方面,本申请提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现所述的正交各向异性介质地震成像方法的步骤。
由上述技术方案可知,本申请提供一种正交各向异性介质地震成像方法及装置,根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数;基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式;应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像,能够有效实现正交各向异性介质射线追踪,追踪过程高效、可靠且准确性高,进而能够更加准确且稳定的实现针对裂缝型储层及勘探开发潜力巨大的非常规储层的地震成像,并能够根据更加准确的地震成像结果进行更为可靠的石油勘探。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例中的正交各向异性介质地震成像系统的结构示意图。
图2为本申请实施例中的正交各向异性介质地震成像方法的流程示意图。
图3为本申请实施例中的正交各向异性介质地震成像方法中步骤300的流程示意图。
图4为本申请具体举例中的正交各向异性介质地震成像方法的流程示意图。
图5为本申请实施例中的正交各向异性介质地震成像装置的结构示意图。
图6为本申请实施例中的正交各向异性介质地震成像装置中地震成像模块30的结构示意图。
图7为本申请实施例中的电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
考虑到当前尚未有一种能够同时保证正交各向异性介质的射线追踪的稳定性和准确性的地震成像方式,本申请提供一种正交各向异性介质地震成像方法及装置,通过利用正交各向异性介质相速度方式,定义了新的含有以慢度和坐标位置为参变量的汉密尔顿函数。根据运动学射线追踪原理,利用汉密尔顿函数推导笛卡尔坐标系与相速度慢度的时间一阶导数具体表达式,给出了正交各向异性介质下的运动学射线追踪方程;根据动力学射线追踪原理,利用汉密尔顿函数推导了波前正交坐标系中的射线中心坐标系转换矩阵的一阶导数具体表达,给出了正交各向异性介质下的动力学射线追踪方程;利用龙格库塔方法和射线的初始值求解计算了正交各向异性介质模型的旅行时场信息,进而可用常规的积分法叠前深度偏移方法完成正交各向异性介质的偏移实现。本申请后续实施例给出的正交各向异性介质的运动学和动力学射线追踪方程的准确、简洁表达式,可以实现裂缝型储层等复杂、强各向异性介质的地震成像。
基于上述内容,本申请还提供一种正交各向异性介质地震成像系统,正交各向异性介质地震成像系统包含有一正交各向异性介质地震成像装置,参见图1,该正交各向异性介质地震成像装置可以为一种服务器1,该服务器1可以与至少一个客户端设备2通信连接,服务器1还可以与石油勘探现场设备通信连接。其中,用户可以通过登录客户端设备2的网站或者APP等,与服务器1建立连接。服务器1可以执行所述正交各向异性介质地震成像的过程,具体根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数;基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式;应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
可以理解的是,所述客户端设备2可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中,所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中,进行正交各向异性介质地震成像的部分可以在如上述内容所述的服务器1侧执行,即,如图1所示的架构,也可以所有的操作都在所述客户端设备2中完成。具体可以根据所述客户端设备2的处理能力,以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备2中完成,所述客户端设备2还可以包括处理器,该处理中用于执行所述正交各向异性介质地震成像方法。
上述的客户端设备可以具有通信模块(即通信单元),可以与远程的服务器进行通信连接,实现与所述服务器的数据传输。所述服务器可以包括任务调度中心一侧的服务器,其他的实施场景中也可以包括中间平台的服务器,例如与任务调度中心服务器有通信链接的第三方服务器平台的服务器。所述的服务器可以包括单台计算机设备,也可以包括多个服务器组成的服务器集群,或者分布式装置的服务器结构。
所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信,包括在本申请提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然,所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol,远程过程调用协议)、REST协议(Representational State Transfer,表述性状态转移协议)等。
为了能够有效实现正交各向异性介质射线追踪,追踪过程高效、可靠且准确性高,进而能够更加准确且稳定的实现针对裂缝型储层及勘探开发潜力巨大的非常规储层的定量刻画,本申请提供一种执行主体为正交各向异性介质地震成像装置的正交各向异性介质地震成像方法的实施例,参见图2,所述正交各向异性介质地震成像方法具体包含有如下内容:
步骤100:根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数。
各向异性介质(anisotropic medium)是物理性质具有方向特性的介质。各向异性意味着不同方向的不同性质。这种介质的物理或机械性能(吸光度,折射率,电导率,拉伸强度等)沿着不同的轴测量时产生差异。各向异性的一个例子是通过偏振器的光,还有例如木材,一个方向比另一方向更容易被折断。地震学家根据晶体对称性的分类体系,按照地下介质中波动物理可实现的对称性,将实际地球介质各向异性基本对称性分为10类,主要表现为弹性系数矩阵的差异,分别为:三斜各向异性介质、单斜各向异性介质、正交各向异性介质(ORT介质)、三方各向异性介质Ⅰ、三方各向异性介质Ⅱ、四方各向异性介质Ⅰ、四方各向异性介质Ⅱ、六方各向异性介质(TI介质)、立方各向异性介质和各向同性介质。
步骤200:基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式。
在步骤200中,所述正交各向异性介质地震成像装置利用定义的汉密尔顿函数求取笛卡尔坐标系下的时间一阶导数具体表达式,即运动学射线追踪系统,T为时间;所述正交各向异性介质地震成像装置利用定义的汉密尔顿函数提出了波前正交坐标系中新的波前正交坐标系转换矩阵的一阶导数具体表达,即动力学射线追踪系统,其中,Q、P为转换矩阵,M=1,2。
步骤300:应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
在步骤300中,所述正交各向异性介质地震成像装置利用正交各向异性介质运动学和动力学射线追踪系统,求解正交各向异性介质下的地震波旅行时场信息,实现积分法叠前深度偏移。
从上述描述可知,本申请实施例提供的正交各向异性介质地震成像方法,是一种高效、高精度的正交各向异性介质射线追踪方法,提出的正交各向异性介质运动学射线追踪方程显示表达式更加简明,计算更加高效;提出了精确的正交各向异性介质动力学射线追踪方程可实现正确的旁轴旅行时计算,进而实现正交各向异性介质的积分法叠前深度偏移。
在本申请的正交各向异性介质地震成像方法的一个实施例中,参见图3,所述正交各向异性介质地震成像方法中的步骤300具体包含有如下内容:
步骤301:应用龙格库塔方法和射线的初始值求解所述运动学射线追踪系统和动力学射线追踪系统,得到所述正交各向异性介质对应的地震波旅行时场信息。
步骤302:根据所述正交各向异性介质对应的地震波旅行时场信息,基于积分法叠前深度偏移方法对所述正交各向异性介质进行偏移处理,得到所述正交各向异性介质的地震成像。
为进一步说明本方案,本申请还提供一种正交各向异性介质地震成像方法的具体应用实例,具体包含有如下内容:
1)定义了正交各向异性介质新的Hamiltonian函数,具体形式如下所示
2)笛卡尔坐标系坐标函数和慢度矢量与时间的一阶导数关系式为下式所示,
在定义的Hamiltonian函数的基础上,利用上式即可求取正交各向异性介质运动学射线追踪方程:
其中,A=1+2ε1,B=1+2ε2,C=1+2δ1,D=1+2δ2,E=1+2δ3,F=ε1-δ1,G=ε2-δ2。
3)利用定义的Hamiltonian函数,可以计算出正交各向异性介质波前正交坐标系下的动力学射线系统,即旁轴射线坐标系统
其中
其中,
通过龙格库塔方法求解以上正交各向异性介质中拟声学情况下的运动学和动力学射线追踪方程即可完成三维情况下射线追踪,获得地下介质旅行时场信息。利用旅行时场信息,最终实现正交各向异性介质的地震成像。
从上述描述可知,本申请应用实例提供的正交各向异性介质地震成像方法,能够推导的正交各向异性介质运动学射线追踪方程更简洁、高效;推导的正交各向异性介质动力学射线追踪方程更精确;可以基于推导的射线追踪方程完成正交各向异性介质的叠前深度偏移方法。
基于上述内容,本申请还提供一种正交各向异性介质地震成像方法的具体举例,参见图4,所述正交各向异性介质地震成像方法具体包含有如下内容:
首先要初始化正交各向异性介质的6个物理参数,Vp0、ε1、δ1、ε2、δ2和δ3,这些参数可以通过测井数据计算获取,也可以通过地震层析及地震全波形反演方法求取。为了验证本方法的准确性,给出了如下的正交各向异性介质值,Vp0=2000m/s,ε1=0.2,δ1=0.06,ε2=0.12,δ2=0.05,δ3=0.1。
三维模型的大小为nx=ny=1001,nz=801,dx=dy=dz=10m,在1-8秒间每隔0.5秒设置雷克子波的脉冲数据,主频为30Hz,炮点位置在地面中心处,检波点位置在地面,空间间隔为40m。在炮点和检波点位置进行正交各向异性介质射线追踪计算,获取地下介质的旅行时信息。
具体射线步骤是在出射点位置给出初始的射线方向利用正交各向异性介质运动学射线追踪方程和动力学运动方程计算射线的空间坐标及转换矩阵的一阶时间导数利用龙格库塔算法计算出射线出射的空间位置和转换矩阵(和),依次完成时间传播中的射线追踪并记录下地下介质的空间旅行时场信息。最后完成不同接收点处的正交各向异性介质射线追踪并记录旅行时场信息。在此基础上,利用常规积分法叠前深度偏移方法,对地震脉冲数据进行偏移。通过定性分析偏移结果,可以看出本文旅行时计算是准确的。
为定量分析偏移结果的准确性,定量计算了地震脉冲响应在正交各向异性介质解析偏移空间的波前时间曲线,通过对比可以发现解析解与正交各向异性介质积分法叠前深度偏移的结果一致,没有偏差出现,定量说明了本申请提供的正交各向异性介质地震成像方法的准确性。
为进一步说明本申请方法的准确性,将本申请的正交各向异性介质地震成像方法退化为横向各向异性介质(VTI)方法并与已有的方法进行对比。给出了如下VTI各向异性介质参数,Vp0=2000m/s,ε1=ε2=0.2,δ1=δ2=0.06,δ3=0。三维模型的大小为nx=ny=1001,nz=801,dx=dy=dz=10m,在1-8秒间每隔0.5秒设置雷克子波的脉冲数据,主频为30Hz,炮点位置在地面中心处,检波点位置在地面,空间间隔为40m,然后在炮点和检波点位置进行VTI各向异性介质射线追踪计算,获取地下介质的旅行时信息。然后利用常规叠前深度偏移方法计算偏移结果。正交各向异性介质退化到VTI介质后即是本申请方法计算的X-line、Crossline和Z切片的偏移结果与不同时刻解析波前曲线做了对比分析之后可知,本申请方法具有广泛的适用范围且准确无误,是一种高效、稳定和高精度的复杂各向异性介质地震成像方法。
也就是说,本申请实例的目的在于提供一种正交各向异性介质稳定、高效射线追踪方法,在定于新的汉密尔顿函数的基础上,定量提出了笛卡尔坐标系中坐标与相速度慢度的时间一阶导数具体表达式,波前正交坐标系中的射线坐标系转换矩阵的一阶导数具体表达,通过求取正交各向异性介质旅行时场信息进行积分法叠前深度偏移。具体可用于裂缝性、非常规页岩储层的正交各向异性介质成像,同时也适用于横向各向同性(VTI)介质成像;同时进行精确的正交各向异性介质运动学和动力学射线追踪计算;可以实现正交各向异性介质旅行时场的精确计算和积分法叠前深度偏移。
为了能够有效实现正交各向异性介质射线追踪,追踪过程高效、可靠且准确性高,进而能够更加准确且稳定的实现针对裂缝型储层及勘探开发潜力巨大的非常规储层的定量刻画,本申请提供一种用于实现所述正交各向异性介质地震成像方法中全部或部分内容的正交各向异性介质地震成像装置的实施例,参见图5,所述正交各向异性介质地震成像装置具体包含有如下内容:
汉密尔顿函数确定模块10,用于根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数。
射线追踪模块20,用于基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式。
地震成像模块30,用于应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
在本申请的正交各向异性介质地震成像装置的一个实施例中,参见图6,所述正交各向异性介质地震成像装置中的地震成像模块30具体包含有如下内容:
地震波旅行时场信息获取单元31,用于应用龙格库塔方法和射线的初始值求解所述运动学射线追踪系统和动力学射线追踪系统,得到所述正交各向异性介质对应的地震波旅行时场信息。
偏移处理单元32,用于根据所述正交各向异性介质对应的地震波旅行时场信息,基于积分法叠前深度偏移方法对所述正交各向异性介质进行偏移处理,得到所述正交各向异性介质的地震成像。
从上述描述可知,本申请应用实例提供的正交各向异性介质地震成像装置,能够推导的正交各向异性介质运动学射线追踪方程更简洁、高效;推导的正交各向异性介质动力学射线追踪方程更精确;可以基于推导的射线追踪方程完成正交各向异性介质的叠前深度偏移方法。
从硬件层面来看,本申请的实施例还提供能够实现上述实施例中的正交各向异性介质地震成像方法中全部步骤的一种电子设备的具体实施方式,参见图7,所述电子设备具体包括如下内容:
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604;
其中,所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信;所述通信接口603用于实现正交各向异性介质地震成像装置、客户终端、现场设备以及其他参与机构之间的信息传输;
所述处理器601用于调用所述存储器602中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的正交各向异性介质地震成像方法中的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
步骤100:根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数。
步骤200:基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式。
步骤300:应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
从上述描述可知,本申请实施例提供的电子设备,能够推导的正交各向异性介质运动学射线追踪方程更简洁、高效;推导的正交各向异性介质动力学射线追踪方程更精确;可以基于推导的射线追踪方程完成正交各向异性介质的叠前深度偏移方法。
本申请的实施例还提供能够实现上述实施例中的正交各向异性介质地震成像方法中全部步骤的一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的正交各向异性介质地震成像方法的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
步骤100:根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数。
步骤200:基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式。
步骤300:应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
从上述描述可知,本申请实施例提供的计算机可读存储介质,能够推导的正交各向异性介质运动学射线追踪方程更简洁、高效;推导的正交各向异性介质动力学射线追踪方程更精确;可以基于推导的射线追踪方程完成正交各向异性介质的叠前深度偏移方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
本领域技术人员应明白,本说明书的实施例可提供为方法、系统或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书的实施例而已,并不用于限制本说明书实施例。对于本领域技术人员来说,本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本说明书实施例的权利要求范围之内。
Claims (10)
1.一种正交各向异性介质地震成像方法,其特征在于,包括:
根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数;
基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式;
应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
3.根据权利要求1所述的正交各向异性介质地震成像方法,其特征在于,所述应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,包括:
应用龙格库塔方法和射线的初始值求解所述运动学射线追踪系统和动力学射线追踪系统,得到所述正交各向异性介质对应的地震波旅行时场信息。
4.根据权利要求1所述的正交各向异性介质地震成像方法,其特征在于,所述根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像,包括:
根据所述正交各向异性介质对应的地震波旅行时场信息,基于积分法叠前深度偏移方法对所述正交各向异性介质进行偏移处理,得到所述正交各向异性介质的地震成像。
5.一种正交各向异性介质地震成像装置,其特征在于,包括:
汉密尔顿函数确定模块,用于根据正交各向异性介质的相速度特征确定对应的汉密尔顿函数;
射线追踪模块,用于基于运动学射线追踪原理,利用所述汉密尔顿函数分别确定所述正交各向异性介质对应的运动学射线追踪系统和动力学射线追踪系统,其中,所述运动学射线追踪系统为笛卡尔坐标系与慢度的时间一阶导数表达式,所述动力学射线追踪系统为波前正交坐标系对应的射线中心坐标系转换矩阵的一阶导数表达式;
地震成像模块,用于应用所述运动学射线追踪系统和动力学射线追踪系统得到所述正交各向异性介质对应的地震波旅行时场信息,并根据该地震波旅行时场信息对所述正交各向异性介质进行偏移处理以实现所述正交各向异性介质的地震成像。
7.根据权利要求5所述的正交各向异性介质地震成像装置,其特征在于,所述地震成像模块包括:
地震波旅行时场信息获取单元,用于应用龙格库塔方法和射线的初始值求解所述运动学射线追踪系统和动力学射线追踪系统,得到所述正交各向异性介质对应的地震波旅行时场信息。
8.根据权利要求5所述的正交各向异性介质地震成像装置,其特征在于,所述地震成像模块包括:
偏移处理单元,用于根据所述正交各向异性介质对应的地震波旅行时场信息,基于积分法叠前深度偏移方法对所述正交各向异性介质进行偏移处理,得到所述正交各向异性介质的地震成像。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至4任一项所述正交各向异性介质地震成像方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至4任一项所述的正交各向异性介质地震成像方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910408859.9A CN111948706B (zh) | 2019-05-16 | 2019-05-16 | 正交各向异性介质地震成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910408859.9A CN111948706B (zh) | 2019-05-16 | 2019-05-16 | 正交各向异性介质地震成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111948706A true CN111948706A (zh) | 2020-11-17 |
CN111948706B CN111948706B (zh) | 2023-06-30 |
Family
ID=73336679
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910408859.9A Active CN111948706B (zh) | 2019-05-16 | 2019-05-16 | 正交各向异性介质地震成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111948706B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
US20100309749A1 (en) * | 2009-06-03 | 2010-12-09 | Terra Nova Sciences Llc | Methods and systems for multicomponent time-lapse seismic measurement to calculate time strains and a system for verifying and calibrating a geomechanical reservoir simulator response |
US20140188393A1 (en) * | 2012-12-27 | 2014-07-03 | King Abdullah University of Science and Technology (KAUST) | Efficient wavefield extrapolation in anisotropic media |
US20140200813A1 (en) * | 2013-01-11 | 2014-07-17 | Cgg Services Sa | Systems and methods for seismic data processing using kinematic analysis of source-receive migration adcigs |
WO2014117284A2 (en) * | 2013-01-30 | 2014-08-07 | Shravan Hanasoge | Wave propagation and imaging method |
WO2015042386A1 (en) * | 2013-09-20 | 2015-03-26 | Westerngeco Llc | Eikonal solver for quasi p-waves in anisotropic media |
CN106547023A (zh) * | 2017-01-16 | 2017-03-29 | 青岛海洋地质研究所 | 一种精度高、计算稳定的复杂介质地震波场延拓方法 |
CN107290780A (zh) * | 2017-08-03 | 2017-10-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 射线方程获取方法、高斯束计算方法及叠前深度偏移方法 |
CN109116413A (zh) * | 2018-07-30 | 2019-01-01 | 中国石油化工股份有限公司 | 成像域立体层析速度反演方法 |
-
2019
- 2019-05-16 CN CN201910408859.9A patent/CN111948706B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
US20100309749A1 (en) * | 2009-06-03 | 2010-12-09 | Terra Nova Sciences Llc | Methods and systems for multicomponent time-lapse seismic measurement to calculate time strains and a system for verifying and calibrating a geomechanical reservoir simulator response |
US20140188393A1 (en) * | 2012-12-27 | 2014-07-03 | King Abdullah University of Science and Technology (KAUST) | Efficient wavefield extrapolation in anisotropic media |
US20140200813A1 (en) * | 2013-01-11 | 2014-07-17 | Cgg Services Sa | Systems and methods for seismic data processing using kinematic analysis of source-receive migration adcigs |
WO2014117284A2 (en) * | 2013-01-30 | 2014-08-07 | Shravan Hanasoge | Wave propagation and imaging method |
WO2015042386A1 (en) * | 2013-09-20 | 2015-03-26 | Westerngeco Llc | Eikonal solver for quasi p-waves in anisotropic media |
CN106547023A (zh) * | 2017-01-16 | 2017-03-29 | 青岛海洋地质研究所 | 一种精度高、计算稳定的复杂介质地震波场延拓方法 |
CN107290780A (zh) * | 2017-08-03 | 2017-10-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 射线方程获取方法、高斯束计算方法及叠前深度偏移方法 |
CN109116413A (zh) * | 2018-07-30 | 2019-01-01 | 中国石油化工股份有限公司 | 成像域立体层析速度反演方法 |
Non-Patent Citations (2)
Title |
---|
VÁCLAV VAVRYČUK: "Velocity, attenuation, and quality factor in anisotropic viscoelastic media: A perturbation approach", 《GEOPHYSICS》 * |
刘强等: "各向异性介质共炮域高斯束偏移", 《石油地球物理勘探》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111948706B (zh) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20140129194A1 (en) | Methods of three-dimensional potential field modeling and inversion for layered earth models | |
CN108983285B (zh) | 一种基于矩张量的多种震源波场模拟方法及装置 | |
US10495768B2 (en) | Method of operating a data-processing system for the simulation of the acoustic wave propagation in the transversely isotropic media comprising an hydrocarbon reservoir | |
BR112016016088B1 (pt) | Método implementado por computador de migração de dados sísmicos | |
JP2020523589A (ja) | リバースタイムマイグレーションからの角度ドメイン共通画像ギャザーのモデリング | |
CN106646611B (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
CN108415073B (zh) | 角度域逆散射偏移成像方法及装置 | |
Stovas et al. | Application of perturbation theory to a P-wave eikonal equation in orthorhombic media | |
Le Bouteiller et al. | A discontinuous Galerkin fast-sweeping eikonal solver for fast and accurate traveltime computation in 3D tilted anisotropic media | |
NO20150821A1 (en) | Efficient Wavefield Extrapolation in Anisotropic Media | |
CN103869362A (zh) | 体曲率获取方法和设备 | |
WO2022232572A1 (en) | Method and system for high resolution least-squares reverse time migration | |
CN111781635B (zh) | 海底四分量弹性波高斯束深度偏移方法和装置 | |
Gao et al. | Parallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach | |
CN104122584A (zh) | 根据地震数据确定方向性的方法及装置 | |
CA2953100A1 (en) | System and method for rock property estimation of subsurface geologic volumes | |
US20220283329A1 (en) | Method and system for faster seismic imaging using machine learning | |
Poppeliers et al. | The relative importance of assumed infrasound source terms and effects of atmospheric models on the linear inversion of infrasound time series at the source physics experiment | |
CN108693557A (zh) | 一种地震数据处理方法及装置 | |
Gulley et al. | A numerical approach for modelling fault-zone trapped waves | |
Sun et al. | Full-azimuth anisotropic prestack time migration in the local-angle domain and its applications on fracture detection | |
CN111948706B (zh) | 正交各向异性介质地震成像方法及装置 | |
CN109001804B (zh) | 一种基于三维地震数据确定有效应力的方法、装置及系统 | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
US20090043549A1 (en) | Methods, apparatus, and products for seismic ray tracing |
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 |