CN113534255B - 一种任意形态不连续面自适应表达的方法 - Google Patents
一种任意形态不连续面自适应表达的方法 Download PDFInfo
- Publication number
- CN113534255B CN113534255B CN202110768611.0A CN202110768611A CN113534255B CN 113534255 B CN113534255 B CN 113534255B CN 202110768611 A CN202110768611 A CN 202110768611A CN 113534255 B CN113534255 B CN 113534255B
- Authority
- CN
- China
- Prior art keywords
- medium
- constitutive relation
- solid
- discontinuous
- free surface
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000012937 correction Methods 0.000 claims abstract description 12
- 239000006185 dispersion Substances 0.000 claims abstract description 7
- 239000007787 solid Substances 0.000 claims description 34
- 239000007788 liquid Substances 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000001788 irregular Effects 0.000 claims description 4
- 238000004088 simulation Methods 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 4
- 238000006073 displacement reaction Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 10
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 241001344923 Aulorhynchidae Species 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011438 discrete method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
-
- 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
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring or predicting earthquakes
-
- 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/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
-
- 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/64—Geostructures, e.g. in 3D data cubes
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/641—Continuity of geobodies
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种任意形态不连续面自适应表达的方法,包括如下步骤:导入正演的模型,包括纵横波速度以及密度,对于各向异性介质,导入各向异性参数,并根据模型参数设置空间步长和时间步长;对自由表面进行阶梯化离散;对自由表面上不同点进行本构关系修正;当模型地表起伏时,利用修正后的本构关系更改模型的第一层参数,这样在不对原始的有限差分法代码进行修改的情况下就可以引入自由表面的影响;将修正后的本构关系带入到位移应力方程中经过一系列的运算可在各异性介质的情况下引入自由表面的影响。本发明所述方法可以准确地对不连续面附近的波场进行数值模拟,有利于对此类地震数据的信息进行提取和分析。
Description
技术领域
本发明属于地震勘探技术领域,尤其涉及一种任意形态不连续面自适应表达的方法。
背景技术
在地震学中,通常把地震传播过程中突然发生变化的地方成为不连续面。根据不连续面两侧介质物理的物理性质(固体,液体和气体),可以将不连续面划分为以下四类:1)固体/固体不连续面;2)固体/液体不连续面;3)固体自由表面(陆地表面),也称固体/真空不连续面;4)液体自由表面(海平面),也称液体/真空不连续面。这类不连续面附近均会发生复杂的波转换现象和强能量的界面波生成。例如,固体/液体不连续面往往用于描述海底界面或者地球内外核边界,在其附近会产生一种沿该界面传播的Scholte wave,以及发生P-S波转换,会产生强振幅面波(Rayleigh wave或love wave),其频散特征在近地表面地震勘探中有重要的意义和用途;液体自由表面用于描述海洋表面,是海洋勘探中多次波产生的原因,其可以被视为噪音进行除去(多次波去噪技术),也可以用于照明角度和区域补偿(多次波成像技术)。
在复杂的介质地震波数值模拟中,往往会遇到以上不连续面的组合,对于这些不连续面的正确数值建模表达会直接影响数值模拟的准确性和精度,特别是对于地表观测地震数据和海底地震仪观测地震数据,其检波器和地震仪直接位于不连续面上或附附近,因此受不连续面的影响更大。能否正确的对不连续面附近的波长进行数值模拟会直接对这类地震数据中提取和分析信息产生影响。
此外,这类不连续面还会破坏一些数值离散方法的完整性,例如有限差分法,伪谱法等基于偏微分方程强解形式的算法。目前工业界和科研界最常用的地震数值模拟技术是有限差分法,其适用范围为参数均匀或光滑变化介质,直接应用于含强阻抗对比的不连续面时会导致计算不稳定,数值假象产生和计算精度降低,进而影响最终模拟结果的准确性。
目前,一般采用直接离散化边界条件,以差分形式来近似边界条件中的控件导数,以固体自由表面的实现为例,有中心差分显示近似格式(Alterman和Karal,1968)、单边显示近似格式(Vidale和Clayton,1986)以及复合近似等格式(Ilan等,1975;Ilan和Loewenthal,1976;Lan和Zhang,2011);另外,在一阶速度-应力方程表达形式下对应力张量和速度矢量分别操作,以固体自由表面的实现为例,有真空格式(Zahradnik等,1993;Bohlen和Saenger,2006)、应力镜像法(Levander,1988;Graves,1996;Kristek等,2002)和一些介质平均方法(Mittet,2002;Xu等,2007);更有,基于Mimetic算子的离散方法,例如dela Puente等(2014),Shragge和Tapley(2017),Konuk和Shragge(2020),Shthi等(2021);除以上三种方法,还有有限元类方法,包含谱yuan法(Komatitsch和Tromp,1999;Komatitsch等,1999;Tromp等,2008;Peter等,2011),这类基于偏微分弱解形式的数值模拟方法,可参考现有的普元发开源程序包:SPECFEM2D和SPECFEM3D。
但是,现有技术存在以下问题:1)应用于高泊松比介质的正演模拟会产生计算不稳定的现象,并且最高只能到达二阶的计算精度。在截面起伏情况下其计算稳定性还会受到截面最大曲率的影响,曲率越大越容易不稳定;2)计算精度低,且受泊松比变化影响;在截面起伏情况下,截面两侧物理场所需的关于界面的对称操作变得及其复杂;3)算法实现复杂,需额外开辟内存控件来存储集合拓扑变量,且计算稳定性还会收到截面最大曲率的影响,曲率越大越容易不稳定;4)相比较于基于有限差分法的数值模拟技术,算法理论和实现复杂,计算量大,内存消耗费多。
发明内容
针对上述背景技术中提到的技术缺陷,为了解决上述题述问题,本发明提供了一种任意形态不连续面自适应表达的方法,其特征在于,包括如下步骤:
S1:导入初始正演模型,包括纵横波的速度和密度,设置震源子波的类型和主频的参数,并根据各模型的参数设置空间步长和时间步长;
S2:根据模型参数选取时间步长和空间步长,设置差分的阶数和震源子波;
S3:当模型表面起伏时,对模型的自由表面进行阶梯化离散,将非规则的起伏地表变成规则的阶梯状网格;
S4:运用介质参数修正法对自由表面上非水平的点进行本构关系和密度的修正;
S5:利用修正后的本构关系来替换自由表面处原来的本构关系,在不对原始的有限差分法代码进行修改的情况下,即可引入自由表面的影响;
S6:对于各向异性介质,将修正后的本构关系带入到有限差分方法的代码中,即可在各向异性介质的情况下引入自由表面的影响。
进一步的,所述自由表面包括液体自由表面和固体自由表面。
进一步的,所述不连续面包括固体/液体不连续面和液体/真空不连续面。
进一步的,对于各异性介质,需要导入各向异性参数。
进一步的,所述异性介质适用于固体自由表面部分。
本发明所述通过介质参数修正实现任意形态不连续面自适应表达的方法,可直接继承与应用于现有的工业界和学术界主流的速度-应力交错网格有限分数值程序,精确贴合不连续面的空间位置,没有传统方法中的版网格偏移误差;通过简单的修正不连续面附近网格点的介质参数来实现,且该修正过程只需一次,计算高效、操作简单;可应用高阶有限差分算子,实现与基于偏微分方程弱解类精确实现方法相当的计算精度;适用于各向异性介质附近的不连续面,且具有泊松比自适应的特点。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面对实施例描述中所需要使用的附图作一个简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的介质参数修正实现任意形态不连续面自适应表达的方法的流程图;
图2为本发明实施例的自由表面离散化示意图;
图3为本发明实施例的固体与空气的位置关系示意图,图a)指空气在自由表面左侧,图b)指空气在自由表面前面,图c)空气在自由表面右侧,图d)指空气在自由表面后面,图e)指空气在自由表面上面。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种任意形态不连续面自适应表达的方法,其特征在于,包括如下步骤:
S1:导入初始正演模型,包括纵横波的速度和密度,设置震源子波的类型和主频的参数,并根据各模型的参数设置空间步长和时间步长;
S2:根据模型参数选取时间步长和空间步长,设置差分的阶数和震源子波;
S3:当模型表面起伏时,对模型的自由表面进行阶梯化离散,将非规则的起伏地表变成规则的阶梯状网格;
S4:运用介质参数修正法对自由表面上非水平的点进行本构关系和密度的修正;
S5:利用修正后的本构关系来替换自由表面处原来的本构关系,在不对原始的有限差分法代码进行修改的情况下,即可引入自由表面的影响,所述影响包括但不限于自由表面引起的地震波的反射;
S6:对于各向异性介质,将修正后的本构关系带入到有限差分方法的代码中,即可在各向异性介质的情况下引入自由表面的影响。
进一步的,所述自由表面包括液体自由表面和固体自由表面;所述不连续面包括固体/液体不连续面和液体/真空不连续面;所述异性介质适用于固体自由表面部分。
具体地,对固体/液体不连续面假设固体介质1对应的本构关系为固体介质2对应的本构关系为根据平均介质思想,可将两个连接的介质等效成一个等效平均介质,其本构关系为同理,对于固体/液体不连续面,假设液体介质1对应的本构关系为液体介质2对应的本构关系为根据平均介质思想,可将两个连接的介质等效成一个等效平均介质,其本构关系为
上述式中:τ指应力张量矩阵;M1指弹性介质1,M2指弹性介质2,A代表平均介质;E指等效的弹性系数矩阵;λ1、λ2和μ1、μ2是拉梅常数;ε为应变张量矩阵。
进一步的,对于固体自由表面,例如陆地表面,即固体/真空不连续面,假设上层真空介质的本构关系为其中根据λ1→0,μ1=0,下层固体介质的本构关系为则等效平均介质的本构关系为λ1→0,μ1=0;
进一步的,对于地下各向异性介质情况,即固体/真空不连续面,假设上层真空介质的本构关系为其中对于下层各向异性介质,其本构关系为则等效平均介质的本构关系为其中
上述公式中c11,c13,c33,c44为各向异性介质的弹性参数。
上述式中:τ指应力张量矩阵;C代表弹性系数;M1指弹性介质1,M2指弹性介质2,A代表平均介质;E指等效的弹性系数矩阵;λ1、λ2和μ1、μ2是拉梅常数;ε为应变张量矩阵。
具体的,在步骤S1中,对于各向异性介质,需要导入各向异性参数c11,c13,c33,c44。同时需要确定差分的阶数、震源类型和震源主频等参数,根据模型参数设置空间步长和时间步长。
进一步的,当模型的地表面起伏时,需要对模型的自由面进行阶梯化离散,离散如图1所述,将非规则的欺负地面变成规则的阶梯状网格,如图2,图2中0L、0R指网格点的外角点,IL、IR代表网格点的内角点,图2显示了起伏地表的实现方式;对自由表面网格点进行划分,并对不同网格点根据表1进行本构关系修正。
在步骤S4中,运用介质参数法对自由表面上非水平的点进行本构关系修正,修正关系如下表所示,其中VL代表自由表面左边为空气的情况,VR代表自由表面右边为空气的情况,VB代表自由表面后边为空气的情况。H代表自由表面上方为空气的情况,具体可参考图3,图3中A指空气,S指地下介质。
表1:起伏地表自由表面处本构关系
上述表中:τij指应力张量的第ij分量;εij指应变张量的第ij分量;λ和μ指拉梅常数;ρ指介质的密度;α=(2μ(λ+μ)/(λ+2μ)),β=(μλ)/(λ+2μ)。
进一步的,利用修正后的本构关系,只需更改模型的第一层参数,即自由表面处的参数,这样在不对原始的有限差分法代码进行修改的情况下就可以引入自由表面的影响。在自由表面处的本构关系如下所示:
τzz|z=0=0,
τxy|z=0=τyx|z=0=μεyx=μεxy,
τyz|z=0=τzy|z=0=0,
τxz|z=0=τzx|z=0=0,
ui是位移的第i个分量,τij指应力张量的第ij分量;εij指应变张量的第ij分量;λ和μ指拉梅常数;ρ指介质的密度。
进一步的,对于各向异性介质,以VTI介质为例,VTI为一种各向异性介质,其修正后的本构关系如下式所示:
修正后的本构关系带入到位移应力方程中可得:
由此即可在各向异性介质的情况下引入自由表面的影响。
本发明所述通过介质参数修正实现任意形态不连续面自适应表达的方法,可直接继承与应用于现有的工业界和学术界主流的速度-应力交错网格有限分数值程序,精确贴合不连续面的空间位置,没有传统方法中的版网格偏移误差;通过简单的修正不连续面附近网格点的介质参数来实现,且该修正过程只需一次,计算高效、操作简单;可应用高阶有限差分算子,实现与基于偏微分方程弱解类精确实现方法相当的计算精度;适用于各向异性介质附近的不连续面,且具有泊松比自适应的特点。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例和落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (3)
1.一种任意形态不连续面自适应表达的方法,其特征在于,包括如下步骤:
S1:导入初始正演模型,包括纵横波的速度和密度,设置震源子波的类型和主频的参数,并根据各模型的参数设置空间步长和时间步长;
S2:根据模型参数选取时间步长和空间步长,设置差分的阶数和震源子波;
S3:当模型表面起伏时,对模型的自由表面进行阶梯化离散,将非规则的起伏地表变成规则的阶梯状网格;
S4:运用介质参数修正法对自由表面上非水平的点进行本构关系和密度的修正;
S5:利用修正后的本构关系来替换自由表面处原来的本构关系,在不对原始的有限差分法代码进行修改的情况下,即可引入自由表面的影响;
S6:对于各向异性介质,将修正后的本构关系带入到有限差分方法的代码中,即可在各向异性介质的情况下引入自由表面的影响;
所述不连续面包括固体/液体不连续面和液体/真空不连续面;
对固体/液体不连续面假设固体介质1对应的本构关系为固体介质2对应的本构关系为根据平均介质思想,可将两个连接的介质等效成一个等效平均介质,其本构关系为
同理,对于固体/液体不连续面,假设液体介质1对应的本构关系为液体介质2对应的本构关系为根据平均介质思想,可将两个连接的介质等效成一个等效平均介质,其本构关系为
固体/真空不连续面,假设上层真空介质的本构关系为其中根据λ1→0,μ1=0,下层固体介质的本构关系为则等效平均介质的本构关系为λ1→0,μ1=0,固体/真空不连续面,假设上层真空介质的本构关系为其中对于下层各向异性介质,其本构关系为则等效平均介质的本构关系为其中c33 —→0,c44 —=0;
上述式中:τ指应力张量矩阵;c代表弹性系数;M1指弹性介质1,M2指弹性介质2,A代表平均介质;E指等效的弹性系数矩阵;λ1、λ2和μ1、μ2是拉梅常数;ε为应变张量矩阵;
对于各向异性介质,需要导入各向异性参数c11,c13,c33,c44。
2.根据权利要求1所述通过介质参数修正实现任意形态不连续面自适应表达的方法,其特征在于:所述自由表面包括液体自由表面和固体自由表面。
3.根据权利要求1所述通过介质参数修正实现任意形态不连续面自适应表达的方法,其特征在于:所述异性介质适用于固体自由表面部分。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110768611.0A CN113534255B (zh) | 2021-07-07 | 2021-07-07 | 一种任意形态不连续面自适应表达的方法 |
US17/393,317 US20230020158A1 (en) | 2021-07-07 | 2021-08-03 | Media parameter-modified method for realizing an adaptive expression of an arbitrary discontinuous surface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110768611.0A CN113534255B (zh) | 2021-07-07 | 2021-07-07 | 一种任意形态不连续面自适应表达的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113534255A CN113534255A (zh) | 2021-10-22 |
CN113534255B true CN113534255B (zh) | 2024-07-16 |
Family
ID=78097967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110768611.0A Active CN113534255B (zh) | 2021-07-07 | 2021-07-07 | 一种任意形态不连续面自适应表达的方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20230020158A1 (zh) |
CN (1) | CN113534255B (zh) |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2329043B (en) * | 1997-09-05 | 2000-04-26 | Geco As | Method of determining the response caused by model alterations in seismic simulations |
EP2164694B1 (en) * | 2007-07-02 | 2012-05-30 | MAGMA Giessereitechnologie GmbH | Method for describing the statistical orientation distribution of particles in a simulation of a mould filling process. |
BRPI1012351B1 (pt) * | 2009-02-16 | 2019-11-05 | Maersk Olie & Gas | método para prever propriedades de uma formação anisotrópica na crosta terrestre |
CN102749643B (zh) * | 2011-04-22 | 2015-06-03 | 中国石油天然气股份有限公司 | 一种面波地震记录的频散响应获取方法及其装置 |
GB2491346A (en) * | 2011-05-23 | 2012-12-05 | Imp Innovations Ltd | Three dimensional full-wavefield seismic tomography for use in mining or in extraction of metalliferous mineral and/or diamond deposits |
CN105676280A (zh) * | 2016-01-21 | 2016-06-15 | 中国矿业大学(北京) | 基于旋转交错网格的双相介质地质数据获取方法和装置 |
US20180210101A1 (en) * | 2016-12-30 | 2018-07-26 | China Petroleum & Chemical Corporation | System and method for seismic inversion |
CN108037532B (zh) * | 2017-12-06 | 2020-05-12 | 中国科学院地质与地球物理研究所 | 用于标定合成地震记录的方法、装置、系统和计算机可读介质 |
CN112424646A (zh) * | 2018-06-10 | 2021-02-26 | 吉奥奎斯特系统公司 | 地震数据解释系统 |
CN111781635B (zh) * | 2019-04-04 | 2023-02-24 | 中国石油天然气集团有限公司 | 海底四分量弹性波高斯束深度偏移方法和装置 |
CN112433246A (zh) * | 2020-11-04 | 2021-03-02 | 郝新源 | 一种地表炮检距道集获取方法及系统 |
-
2021
- 2021-07-07 CN CN202110768611.0A patent/CN113534255B/zh active Active
- 2021-08-03 US US17/393,317 patent/US20230020158A1/en active Pending
Non-Patent Citations (1)
Title |
---|
弹性波有限差分模拟中自由表面的自适应表达;曹健;《石油物探》;20180725;第1-3节 * |
Also Published As
Publication number | Publication date |
---|---|
CN113534255A (zh) | 2021-10-22 |
US20230020158A1 (en) | 2023-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107894612B (zh) | 一种q吸收衰减补偿的声波阻抗反演方法及系统 | |
CN106033124B (zh) | 一种基于随机最优化的多震源粘声最小二乘逆时偏移方法 | |
CA2920008A1 (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
CN111025387B (zh) | 一种页岩储层的叠前地震多参数反演方法 | |
CN110542928A (zh) | 基于vti各向异性传播矩阵的地震响应模拟方法 | |
CN109212599A (zh) | 一种地震数据的叠前同步反演方法 | |
CN113031068A (zh) | 一种基于反射系数精确式的基追踪叠前地震反演方法 | |
CN111505714B (zh) | 基于岩石物理约束的弹性波直接包络反演方法 | |
CN114895351A (zh) | 模拟地震波在任意不连续界面处传播的介质模型化的方法及装置 | |
CN116859460B (zh) | 一种适用于弹性波近场反射的海底物性参数反演方法 | |
Jiang et al. | Full waveform inversion based on inversion network reparameterized velocity | |
CN113534255B (zh) | 一种任意形态不连续面自适应表达的方法 | |
CN111273349B (zh) | 一种用于海底浅部沉积层的横波速度提取方法及处理终端 | |
Zhang et al. | Attenuation compensation for wavefield‐separation‐based least‐squares reverse time migration in viscoelastic media | |
CN111273346B (zh) | 去除沉积背景的方法、装置、计算机设备及可读存储介质 | |
WO2018067014A1 (en) | Seismic modeling | |
US8437219B2 (en) | Correcting an acoustic simulation for elastic effects | |
Chapman et al. | Elastic corrections to acoustic finite-difference simulations | |
CN115079263A (zh) | 高精度提取体积模量的地震叠前反演方法 | |
Li et al. | Seismic acoustic impedance inversion using total variation with overlapping group sparsity | |
CN112099079B (zh) | 自适应分频串联反射率反演方法及系统 | |
CN113985480B (zh) | 一种基于角度校正的avo反演方法及装置 | |
CN117388921B (zh) | 弹性参数的叠前反演方法、装置和电子设备 | |
CN117607954B (zh) | 一种基于指数变密度模型的自由参数与基底同时反演方法 | |
CN112684505B (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 |