CN105807317A - 基于切比雪夫伪谱法的各向异性衰减面波模拟方法 - Google Patents
基于切比雪夫伪谱法的各向异性衰减面波模拟方法 Download PDFInfo
- Publication number
- CN105807317A CN105807317A CN201610299398.2A CN201610299398A CN105807317A CN 105807317 A CN105807317 A CN 105807317A CN 201610299398 A CN201610299398 A CN 201610299398A CN 105807317 A CN105807317 A CN 105807317A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- wave
- equation
- sigma
- wave field
- 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
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. analysis, for interpretation, for correction
Abstract
本发明公开基于切比雪夫伪谱法的各向异性衰减面波模拟方法,包括:获取浅层地球物理参数数据、定义观测系统和震源函数;根据自由界面垂向和切向应力连续条件,由波动方程获得自由界面处面波波场更新方程;在垂直和水平方向采用切比雪夫和傅里叶伪谱法计算自由界面处面波波场更新方程的垂向和水平导数;四阶龙格‑库塔离散自由界面处面波波场更新方程的各向异性衰减介质面波方程的时间微分;计算模型水平方向两个边界的波场吸收;将波动方程分解为入射波方程和出射波方程,并得到底边界的波场更新方程,且据以计算当前时刻的面波,以产生面波炮记录和波场快照。通过本发明,以解决现有技术存在的浅层各向异性衰减面波模拟问题。
Description
技术领域
本发明涉及地震勘探的技术领域,尤其涉及一种基于切比雪夫伪谱法的各向异性衰减面波模拟方法。
背景技术
高精度三维地震技术在各个油田得以成功的推进和创新,我国东部老油田,例如,胜利油田、大庆油田等均处于平原地区,而其他主要产区,例如塔河、鄂尔多斯盆地多为沙漠、黄土塬等复杂地表区,探区内资源潜力大,巨厚、疏松的黄土层和浅层沙漠地层不利于有效地震波的产生,且浅地表对地震波的吸收衰减作用强烈,严重影响地震数据处理。
但浅层面波发育,面波波场包含浅地表横波速度、介质对地震波的吸收强弱、各向异性等信息,通过弹性介质的面波已经可实现反演浅层横波速度信息,并可将反演结果应用到静校正流程,但对面波的各向异性衰减研究有限,浅层各向异性衰减的面波模拟研究,可用于建立浅层介质衰减和速度各向异性关系,是我国黄土塬、沙漠地区复杂地质三维地震数据处理技术的理论基础。
发明内容
本发明的主要目的在于提供一种,以解决现有技术存在的浅层各向异性衰减面波模拟的问题。
为解决上述问题,本发明实施例提供一种基于切比雪夫伪谱法的各向异性衰减面波模拟方法,包括:获取浅层地球物理参数数据、定义观测系统和震源函数;根据自由界面垂向和切向应力连续条件,由波动方程获得自由界面处面波波场更新方程;在垂直方向采用切比雪夫伪谱法计算所述自由界面处面波波场更新方程的垂向导数;在水平方向采用傅里叶伪谱法计算所述自由界面处面波波场更新方程的水平导数;所述自由界面处面波波场更新方程的各向异性衰减介质面波方程的时间四阶龙格-库塔离散计算;所述自由界面处面波波场更新方程的模型水平方向的两个边界波场吸收计算;将所述波动方程分解为入射波方程和出射波方程,并将入射波在底边界吸收掉,从而得到底边界的波场更新方程;根据所述底边界的波场更新方程,计算当前时刻的面波,以产生面波炮记录和波场快照;输出所述面波炮记录和波场快照。
根据本发明的技术方案,通过获取浅层地球物理参数数据、定义观测系统和震源函数,再根据自由界面垂向和切向应力连续条件,获得自由界面处面波波场的更新方程,接着,垂直方向采用切比雪夫伪谱法计算导数,水平方向采用傅里叶伪谱法计算导数,各向异性衰减介质面波的时间离散计算,模型水平方向的左右边界吸收计算,而模型的底边界,采用类似于自由边界的方法,将反射到模型内部的地震波场吸收掉,从而得到底界面的波场更新方程,且根据底界面的波场更新方程计算面波波场,并输出面波炮记录和波场快照。如此一来,有效地解决了浅层各向异性衰减面波模拟问题,尤其适用于各向异性衰减介质面波的正演分析工作,还能够快速实现各向异性衰减介质面波模拟,准确记录面波正演结果,并提高地质解释的精度。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的基于切比雪夫伪谱法的各向异性衰减面波模拟方法的流程图;
图2是根据本发明实施例的切比雪夫伪谱法网格剖分的示意图;
图3是根据本发明实施例的一维切比雪夫不均匀网格剖分的示意图;
图4a和图4b分别是根据本发明实施例的勘探区采集的微测井数据的示意图;
图5a和图5b分别根据本发明实施例的面波炮记录的示意图;
具体实施方式
本发明的主要思想在于,基于通过获取浅层地球物理参数数据、定义观测系统和震源函数,再根据自由界面垂向和切向应力连续条件,获得自由界面处面波波场的更新方程,接着,垂直方向采用切比雪夫伪谱法计算导数,水平方向采用傅里叶伪谱法计算导数,各向异性衰减介质面波的时间离散计算,模型水平方向的左右边界吸收计算,而模型的底边界,采用类似于自由边界的方法,将反射到模型内部的地震波场吸收掉,从而得到底界面的波场更新方程,且根据底界面的波场更新方程计算面波波场,并输出面波炮记录和波场快照。如此一来,有效地解决了浅层各向异性衰减面波模拟问题,尤其适用于各向异性衰减介质面波的正演分析工作,还能够快速实现各向异性衰减介质面波模拟,准确记录面波正演结果,并提高地质解释的精度。
为使本发明的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本发明做进一步地详细说明。
根据本发明的实施例,提供了一种基于切比雪夫伪谱法的各向异性衰减面波模拟方法。
图1是根据本发明实施例的基于切比雪夫伪谱法的各向异性衰减面波模拟方法的流程图。
在步骤S102中,获取浅层地球物理参数数据、定义观测系统和震源函数。其中,所述浅层地球物理参数数据包括纵波速度、横波速度、Thomson各向异性参数、纵波品质因子和横波品质因子。所述观测系统包括震源位置、检波器起始位置、检波器个数、检波器间距,定义地震模拟的总时间长度和时间采样率,记录波场快照的时间间隔。所述震源函数包括震源类型和震源主频,并且震源类形还包括爆炸震源和点震源。
进一步来说,震源类型可如式(1)、式(2)所示:
其中,表示体力项,Z表示体力的垂向分量,分别为x和y方向的单位向量,φ表示一标量函数,是关于时间的函数,x、y分别表示相互垂直的二维坐标轴。
在步骤S104中,根据自由界面垂向和切向应力连续条件,由波动方程获得自由界面处面波波场更新方程。
首先,根据各向异性衰减理论,采用标准线性模型表征介质粘弹性,通过粘弹性波动方程中引入的记忆变量表征介质衰减效应,二维直角坐标系下,各向异性衰减介质qP-qSV波动方程,如式(3)、(4)、(5)和(6)所示:
其中,vx、vz为质点在x、z方向的速度分量,σxx、σzz、σxz为应力分量,ρ为介质的密度,fx、fz表示外力;
其中,
其中,表示未松弛的介质弹性系数,cij表示松弛的介质弹性系数, 为应变的时间导数,e1l、e2l、e3l为记忆变量,Muv,v=1,2为广义标准线性模型的松弛函数,v=1对应P波松弛函数,v=2对应SV波松弛函数,X、Y为常数;
其中,为介质的应变与应力松弛时间,为记忆变量的一阶、二阶导数。
然后,根据各向异性衰减波动方程对应的一阶双曲方程,将波动方程在自由界面的波场分解为入射分量和出射分量,从已知波场计算入射波,并依据自由界面σzz=σxz=0,相应的入射波对应的特征值为λ2=-cp,λ4=-cs,则可获得自由界面面波波场更新方程,如式(7)所示:
其中,上标(new)表示重新计算的波场,上标(old)表示最初由波动方程计算的波场。
在步骤S106中,在垂直方向采用切比雪夫伪谱法计算所述自由界面处面波波场更新方程的垂向导数。通过切比雪夫伪谱法计算,式(3)、(4)、(6)可转化为式(8)、(9)、(10),如下所示:
其中,式(8)、式(9)和式(10)实现了垂向导数的切比雪夫网格剖分计算,这种网格在模型边界的网格间距小,模型中间网格间距最大,利于提高自由界面面波模拟精度。
进一步的,切比雪夫不规格网格剖分图如图2所示,在垂直方向的两个边界网格密集,可以提高面波模拟的精度,而且式(7)所示的面波更新方程尤其适用于这种边界。切比雪夫不规格网格剖分的区间是[-1,+1]之间,如图3所示,在区间内由余弦函数定义式(11),如下所示:
通过切比雪夫多项式离散函数f(x)后,可得式(12),如下所示:
并且,函数f(xi)关于x的空间导数如式(13)所示:
其中,bj、dk为其系数,Tj(xi)为切比雪夫函数,Tk(xi)为余弦函数。切比雪夫网格剖分区间为引入扩展函数,并将区间[-1,1]映射到实际物理区间[0,zmax],则在实际物理区间内,函数的空间如式(14)所示:
其中,为区间扩展函数。
在步骤S108中,在水平方向采用傅里叶伪谱法计算所述自由界面处面波波场更新方程的水平导数。在本实施例中,假设某函数为f(x),根据傅里叶变换性质,水平导数计算过程如式(15)所示:
其中,FFT、FFT-1分别为正反傅里叶变换,F(w)为f(x)的傅里叶变换,kx为波数,通过式(15)可计算函数f(x)关于x的空间导数。
在步骤S110中,所述自由界面处面波波场更新方程的各向异性衰减介质面波方程的时间四阶龙格-库塔离散计算。在本实施例中,各向异性衰减介质面波方程的时间离散计算是采用四阶龙格-库塔方法。
首先,将各向异性衰减介质波动方程改写成式(16),如下所示:
其中,矢量V为{vx,vz,σxx,σzz,σxz,e1l,e2l,e3l}T,分别为关于时间、方向x和z的偏导数,A、B表示弹性系数和函数φ1l、φ2l的矩阵,D表示波动方程(9)中的常量矩阵。
然后,各向异性衰减介质波动方程的时间微分表达式变为式(17),如下所示:
并且,式(17)中的参数如式(18)所示:
其中,H1表示Vn点处的斜率;H2表示利用Δ1求得的(Vn+dt/2)点处的斜率;H3表示利用H2求得的(Vn+dt/2)点处的斜率;H4表示利用Δ3求得的(Vn+dt)点处的斜率;M表示速度在空间三个方向的导数算子,Dn表示各向异性衰减方程中常量在开始时刻(第n时刻)的数值,Dn+1/2表示常量在时间段中点(第n+1/2时刻)的数值,Dn+1表示常量在时间段终点(第n+1时刻)的数值。
在步骤S112中,所述自由界面处面波波场更新方程的模型水平方向的两个边界波场吸收计算。在本实施例中,模型水平方向的两个边界波场吸收计算是采用卷积完全匹配层(CPML)吸收边界法,如式(19)所示:
其中,参数a和b分别为边界吸收系数,u为匹配层处的地震波场,ξx(t)为匹配层内线性变换的衰减函数。
在步骤S114中,将所述波动方程分解为入射波方程和出射波方程,并将入射波在底边界吸收掉,从而得到底边界的波场更新方程。其中,所述底边界的波场更新方程如式(20)所示:
其中,cp、cs分别为弹性情况下纵波和横波速度,上标(new)表示重新计算的波场,上标(old)表示最初由波动方程计算的波场。
在步骤S116中,根据所述底边界的波场更新方程,计算当前时刻的面波,以产生面波炮记录和波场快照。在步骤S118中,输出所述面波炮记录和波场快照。
上述已说明了基于切比雪夫伪谱法的衰减各向异性面波模拟方法,以下将提供一些实例来验证上述方法的适用性及准确性。
图4a和图4b是根据本发明实施例的勘探区采集的微测井数据的示意图。利用微测井单道记录,拾取初至,并转化为垂直的t0时间,再与距离绘制时深曲线,然后,利用时深曲线计算、划分速度界面,从而得到近地表速度模型,模型定义为vP=1826m/s,vS=1000m/s,ρ=2300g/cm3,QP=100,QS=60。Thomson各向异性参数可以选择为:ε=0.195,δ=-0.22;模型计算尺寸为1000×1000m,水平方向网格间距dx=1m,垂直方向最大网格间距dz=2m,震源位于Nx=500m和Nz=500m,震源类型采用爆炸震源,子波主频为45Hz。通过本实施例的上述步骤进行模拟后,可获得面波炮记录,如图5a所示。在图5a和图5b中,R表示各向异性衰减介质的面波波场,ux为水平分量,uz为垂向分量。
综上所述,根据本发明的技术方案,通过获取浅层地球物理参数数据、定义观测系统和震源函数,再根据自由界面垂向和切向应力连续条件,获得自由界面处面波波场更新方程,将当前时刻的地震波场更新为各向异性衰减面波的传播方程,接着,垂直方向采用具有不规格网格剖分的切比雪夫伪谱法计算导数,水平方向采用规格网格剖分的傅里叶伪谱法计算导数,各向异性衰减介质面波的时间微分采用四阶龙格-库塔法离散,针对二维模型水平方向的左右边界,采用非分裂CPML吸收较强的边界反射,而模型的底边界,采用类似于自由边界的方法,将反射到模型内部的地震波场吸收掉,从而得到底界面的波场更新方程,以上步骤完成后,叠代计算面波波场,并输出面波炮记录和波场快照。如此一来,有效地解决了浅层各向异性衰减面波模拟问题,尤其适用于各向异性衰减介质面波的正演分析工作,还能够快速实现各向异性衰减介质面波模拟,准确记录面波正演结果,并提高地质解释的精度。
显然,本领域的技术人员应该明白,上述的本发明的基于切比雪夫伪谱法的衰减各向异性面波模拟方法和各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行。这样,本发明不限制于任何特定的硬件和软件结合。存储装置为非易失性存储器,如:ROM/RAM、闪存、磁盘、光盘等。
以上所述仅为本发明的实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (6)
1.一种基于切比雪夫伪谱法的各向异性衰减面波模拟方法,其特征在于,包括以下步骤:
获取浅层地球物理参数数据、定义观测系统和震源函数;
根据自由界面垂向和切向应力连续条件,由波动方程获得自由界面处面波波场更新方程;
在垂直方向采用切比雪夫伪谱法计算所述自由界面处面波波场更新方程的垂向导数;
在水平方向采用傅里叶伪谱法计算所述自由界面处面波波场更新方程的水平导数;
所述自由界面处面波波场更新方程的各向异性衰减介质面波方程的时间四阶龙格-库塔离散计算;
所述自由界面处面波波场更新方程的模型水平方向的两个边界波场吸收计算;
将所述波动方程分解为入射波方程和出射波方程,并将入射波在底边界吸收掉,从而得到底边界的波场更新方程;
根据所述底边界的波场更新方程,计算当前时刻的面波,以产生面波炮记录和波场快照;
输出所述面波炮记录和波场快照。
2.根据权利要求1所述的基于切比雪夫伪谱法的各向异性衰减面波模拟方法,其特征在于,所述浅层地球物理参数数据包括:纵波速度、横波速度、Thomson各向异性参数、纵波品质因子和横波品质因子。
3.根据权利要求1所述的基于切比雪夫伪谱法的各向异性衰减面波模拟方法,其特征在于,所述定义观测系统包括:震源位置、检波器起始位置、检波器个数、检波器间距。
4.根据权利要求1所述的基于切比雪夫伪谱法的各向异性衰减面波模拟方法,其特征在于,所述震源函数包括:震源类型和震源主频。
5.根据权利要求1所述的基于切比雪夫伪谱法的各向异性衰减面波模拟方法,其特征在于,所述自由界面处面波波场更新方程满足如下公式:
其中,vx、vz为质点在x、z方向的速度分量,σxx、σzz、σxz为应力分量,ρ为介质的密度,fx、fz表示外力,表示未松弛的介质弹性系数,e1l、e2l、e3l为记忆变量,分别为衰减介质的纵波和横波速度,上标(new)表示重新计算的波场,上标(old)表示最初由波动方程计算的波场。
6.根据权利要求5所述的基于切比雪夫伪谱法的各向异性衰减面波模拟方法,其特征在于,所述底边界的波场更新方程满足如下公式::
其中,vx、vz为质点在x、z方向的速度分量,σxx、σzz、σxz为应力分量,ρ为介质的密度,fx、fz表示外力,表示未松弛的介质弹性系数,e1l、e2l、e3l为记忆变量,cp、cs分别为弹性情况下纵波和横波速度,上标(new)表示重新计算的波场,上标(old)表示最初由波动方程计算的波场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610299398.2A CN105807317B (zh) | 2016-05-06 | 2016-05-06 | 基于切比雪夫伪谱法的各向异性衰减面波模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610299398.2A CN105807317B (zh) | 2016-05-06 | 2016-05-06 | 基于切比雪夫伪谱法的各向异性衰减面波模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105807317A true CN105807317A (zh) | 2016-07-27 |
CN105807317B CN105807317B (zh) | 2019-04-30 |
Family
ID=56455448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610299398.2A Active CN105807317B (zh) | 2016-05-06 | 2016-05-06 | 基于切比雪夫伪谱法的各向异性衰减面波模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105807317B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109307888A (zh) * | 2017-07-28 | 2019-02-05 | 中国石油化工股份有限公司 | 波场模拟的自由面处理方法及系统 |
CN110609325A (zh) * | 2018-06-15 | 2019-12-24 | 中国石油化工股份有限公司 | 弹性波场数值模拟方法及系统 |
CN110646839A (zh) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | 地震波模拟方法及系统 |
CN111208514A (zh) * | 2020-01-06 | 2020-05-29 | 华侨大学 | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101840001A (zh) * | 2010-02-10 | 2010-09-22 | 中国科学院地质与地球物理研究所 | 一种地质构造三维成像数据的获取方法及装置 |
CN103530451A (zh) * | 2013-09-27 | 2014-01-22 | 中国科学院声学研究所 | 复杂介质弹性波传播模拟的多网格切比雪夫并行谱元法 |
CN105204064A (zh) * | 2015-10-09 | 2015-12-30 | 西南石油大学 | 一种基于优化系数的混合域傅里叶有限差分偏移方法 |
-
2016
- 2016-05-06 CN CN201610299398.2A patent/CN105807317B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101840001A (zh) * | 2010-02-10 | 2010-09-22 | 中国科学院地质与地球物理研究所 | 一种地质构造三维成像数据的获取方法及装置 |
CN103530451A (zh) * | 2013-09-27 | 2014-01-22 | 中国科学院声学研究所 | 复杂介质弹性波传播模拟的多网格切比雪夫并行谱元法 |
CN105204064A (zh) * | 2015-10-09 | 2015-12-30 | 西南石油大学 | 一种基于优化系数的混合域傅里叶有限差分偏移方法 |
Non-Patent Citations (4)
Title |
---|
EKKEHART TESSMER 等: "Elastic wave propagation simulation in the presence of surface topography", 《GEOPHYSICS》 * |
EKKEHART TESSMER 等: "Multi-domain Chebyshev-Fourier Method for the Solution of the Equations of Motion of Dynamic Elasticity", 《JOURNAL OF COMPUTATIONAL PHYSICS》 * |
EKKEHART TESSMER: "3-D seismic modelling of general material anisotropy in the presence of the free surface by a Chebyshev spectral method", 《GEOPHYSICS》 * |
刘鲁波 等: "切比雪夫伪谱法模拟地震波场", 《西北地震学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109307888A (zh) * | 2017-07-28 | 2019-02-05 | 中国石油化工股份有限公司 | 波场模拟的自由面处理方法及系统 |
CN110609325A (zh) * | 2018-06-15 | 2019-12-24 | 中国石油化工股份有限公司 | 弹性波场数值模拟方法及系统 |
CN110646839A (zh) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | 地震波模拟方法及系统 |
CN111208514A (zh) * | 2020-01-06 | 2020-05-29 | 华侨大学 | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 |
CN111208514B (zh) * | 2020-01-06 | 2023-05-05 | 华侨大学 | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105807317B (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ikelle et al. | Introduction to petroleum seismology | |
Zhu et al. | Modeling acoustic wave propagation in heterogeneous attenuating media using decoupled fractional Laplacians | |
Xu et al. | Accurate simulations of pure quasi-P-waves in complex anisotropic media | |
Baig et al. | Microseismic moment tensors: A path to understanding frac growth | |
Askan et al. | Full waveform inversion for seismic velocity and anelastic losses in heterogeneous structures | |
Pan et al. | Love-wave waveform inversion in time domain for shallow shear-wave velocity | |
Liu et al. | A new nonlinear finite fault inversion with three‐dimensional Green's functions: Application to the 1989 Loma Prieta, California, earthquake | |
US20170067323A1 (en) | Multi data reservoir history matching and uncertainty quantification framework | |
Gosselin-Cliche et al. | 3D frequency-domain finite-difference viscoelastic-wave modeling using weighted average 27-point operators with optimal coefficients | |
CN108983285B (zh) | 一种基于矩张量的多种震源波场模拟方法及装置 | |
US20210405233A1 (en) | System and method for real-time passive seismic event localization | |
CN103869362B (zh) | 体曲率获取方法和设备 | |
CN106353797A (zh) | 一种高精度地震正演模拟方法 | |
CN105807317A (zh) | 基于切比雪夫伪谱法的各向异性衰减面波模拟方法 | |
Ren et al. | A hybrid absorbing boundary condition for frequency-domain finite-difference modelling | |
Wang et al. | In situ assessment of soil dynamic parameters for characterizing nonlinear seismic site response using KiK-net vertical array data | |
Boaga et al. | From surface wave inversion to seismic site response prediction: Beyond the 1D approach | |
Homma et al. | A physics-based Monte Carlo earthquake disaster simulation accounting for uncertainty in building structure parameters | |
Taborda et al. | Physics-based ground-motion simulation | |
Narayan et al. | A fourth-order accurate finite-difference program for the simulation of SH-wave propagation in heterogeneous viscoelastic medium | |
Safarshahi et al. | Stochastic finite‐fault modeling of ground motion for the 2010 Rigan earthquake, southeastern Iran | |
CN106501872B (zh) | 一种裂缝储层地应力特征的计算方法及装置 | |
Jie et al. | Time-reverse location of microseismic sources in viscoelastic orthotropic anisotropic medium based on attenuation compensation | |
Raghu Kanth et al. | Stochastic finite fault modeling of subduction Zone Earthquakes in Northeastern India | |
Romero et al. | Seeing through the CO2 plume: Joint inversion-segmentation of the Sleipner 4D seismic data set |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |