共查询到20条相似文献，搜索用时 234 毫秒
1.
产生地震动空间效应的主要原因是地质构造的变化会极大地改变复杂场地条件下到达不同观测点的地震波,因此地震波振幅和相位发生明显变化。本文利用有限元和边界元耦合的方法,通过建立雷克子波入射时的非线性响应模型求解复杂场地条件下的波动积分方程,进行复杂近地表构造中地震波场数值模拟,旨在分析复杂场地条件下入射方向不同的地震波的动力响应规律。研究结果显示近地表会极大地改变地震波传播特征,导致地震波传播趋向不一致性的因素是由于复杂场地条件重建了地震波的传播路径,尤其须指出的是非一致性地表位移响应源于软土层中低频子波的传播。在地震响应分析中导致系统不同动力响应的其他原因是主频和地震波速度结构的改变,同时地震波入射角度的改变会引起到达观测点的时间滞后,地表位移响应随着入射角度的增大而增加,地震动的空间效应随着入射角度的减小而愈趋明显。 相似文献
2.
本文针对声-弹耦合介质,为尽可能的减少频率域正演模拟的计算内存,提高计算效率,在一阶非均质位移-应力波动方程的基础上,借助等效交错网格思想并充分考虑密度参数空间变化对地震波传播的影响,推导了声-弹耦合地震波波动方程.在流相介质和固相介质中分别采用非均质情况频率域二阶声压标量波、二阶纯位移控制方程,为保证流、固相介质间地震波能量的稳定传输和有效交换,提出了声-弹耦合界面转换过渡层方法,并详细阐述了过渡层与上下介质空间差分具体耦合方法.在与非均质纯位移波动方程正演结果对比分析的基础上,首先采用各向同性单层流相介质模型进行正演模拟验证了声-弹耦合方程数值模拟中过渡层策略的有效性和准确性,随后又数值模拟了地震波在声-弹耦合介质简单模型和复杂Marmousi2模型中的传播,验证了本文方法稳定性和准确性,同时该方法可以简单的推广到三维情况. 相似文献
3.
偏微分方程作为一种描述客观物理规律的重要工具,可应用于解决大多数科学问题和工程问题.在地球物理学中,求解波动方程是地震波场数值模拟的重要步骤,对于偏移成像和反演起着至关重要的作用.随着勘探目标复杂化,地震波方程的数学描述也更加复杂,现有的经典地震波数值模拟方法存在一些待完善之处,需要一些新型正演算法的应用.物理信息神经网络(Physics-Informed Neural Network,简称PINN)是一种利用神经网络求解偏微分方程近似解的有效手段,通过自动微分将偏微分方程及其定解条件引入损失函数中,在最小化数据驱动误差的同时,使其满足方程的物理约束条件.本文在分析PINN原理并总结其研究现状的基础上,对其在波动方程正演中已有成果进行归纳和分析,展望了其未来可能的研究方向,并对PINN在一维声波方程正演中的应用进行了探究. 相似文献
4.
5.
页岩的强各向异性特征挑战地震波传播数值模拟方法的精度极限, 特别是易引起频散的高频波(>100 Hz)传播的数值模拟.鉴于目前我国页岩气地震勘探主要以常规地震声波资料为主,本文首先介绍了一种VTI介质声波方程的任意偶数阶有限差分数值模拟方法,并讨论其稳定性条件和吸收边界条件.任意偶数阶的差分解可有效提高计算精度,压制数值频散噪声.针对页岩较强的各向异性特征,本文比较了不同模型的声波方程和VTI介质声波方程计算得到的地震响应.数值结果表明,各向异性对地震波的运动学(相位)和动力学(振幅)特性影响作用明显.因此,在页岩气地震勘探资料处理的各个环节必须充分考虑各向异性的影响,采取有别于常规油气勘探的处理流程和技术. 相似文献
6.
传统的人工边界是在声波或弹性波分解为内行和外行分量的基础上设计的.对于较复杂的介质模型（如粘弹介质、弹塑介质、塑性介质模型等）,地震波不能简单地分解为各自满足其运动方程的纵波、横波等成分,且地震波在传播过程中还有衰减现象,因此传统的人工边界已不能有效地用于非弹性的较复杂介质的波场模拟.本文提出了一种适合于更一般的地震模型的波场计算的透射边界条件,并以粘弹介质的一个较简单的二维几何模型为例,通过有限元数值求解这种透射边界和刚性边界下粘弹波的激发和传播,得到地面合成地震记录,证明了这种透射边界用于较复杂介质波场计算是可行的. 相似文献
7.
地震波理论研究进展--介质模型与地震波传播 总被引：17，自引：11，他引：6
介质模型及地震波传播是地震波理论的核心内容之一,本文结合作者目前研究所涉及的内容,阐述了介质模型及地震波传播的相关问题．首先是地震学背景中的地震波传播的主要问题．其次是地震波传播的动力学,其中的主体内容可以由波动方程的建立、波动方程和波动方程的延伸三部分组成．再次是地震波的传播与介质因素和界面因素,以及孔隙介质模型的认识和讨论．最后结合均匀弹性各向同性、粘弹性各向同性、双相介质和各向异性四种介质模型,阐述了地震波传播的主要特性． 相似文献
8.
《地球物理学进展》2016,(2)
地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段,而且是研究地球深部精细构造和地球深部探测的有效工具.射线理论和波动方程理论是地震波场数值模拟的理论基础.射线理论主要刻画地震波在介质中传播的走时场、地震射线等运动学属性;波动方程理论通过求解波动方程来描述地震波在介质中传播的弹性动力学响应(能量衰减、相位特征、偏振属性、以及全波形等).基于波动方程理论的波场数值模拟由于能够引入丰富的波场信息,使得人们对不同介质中地震波的传播过程有了较全面的了解.本文以二维层状均匀介质模型为例,通过射线追踪法和交错网格有限差分法模拟得到的波场快照图、单炮地震记录剖面、合成理论地震图的分析比较,不但对地震波在各向同性层状均匀介质中的传播规律和特点有了深刻的理解和认识,同时又可以相互验证两种不同方法的正确性和有效性. 相似文献
9.
波动方程深度偏移是解决复杂地质体成像的关键技术,基于波动方程的速度建模为其提供更为精确的速度模型.频率域波形反演是目前研究最为广泛的波动方程速度建模方法之一,它推动了波形反演在勘探尺度下的应用.本文通过对频率域波形反演的实现,分析对比了其有效执行过程中与频率相关的影响因素.介绍了时间域的多尺度反演方法在频率域的一种实现方式,对比分析了输入数据的频点带宽和应用的子波频带范围不同时对反演结果的影响.本文通过设计的山地地质模型对频率域波形反演进行了测试和对比,得到的结论为频率域波形反演的有效计算提供了依据和参考. 相似文献
10.
本文基于Kjartansson常Q模型理论,推导了常Q衰减介质中黏声波和黏弹性波的速度-应力方程,并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟,同时引入不分裂的复频移卷积完全匹配层(CPML)吸收边界条件,以消除边界反射.使用基于自适应时间步长记忆方法的中心差分近似时间分数阶导数,与常用的短时记忆方法相比,提高了波动方程的离散化精度和计算效率.通过对比均匀模型下声波的数值解与解析解,验证了算法的精确性,并进一步分析了不同品质因子下地震波的频散及衰减特征.对BP盐丘模型的数值模拟结果可以较好地反映本文数值方法对复杂介质的适应性及频散压制效果. 相似文献
11.
本文首先分析了地震波在黏弹介质的传播规律,基于黏弹介质地震波动方程总结了时变子波振幅谱和相位谱的关系,从而得出结论,准确估计子波相位谱初值和不同时刻的子波振幅谱是实现时变子波准确提取的必要条件.在此基础上,针对传统方法限制子波振幅谱形态且受限于分段平稳假设的问题,提出了一种利用EMD(Empirical Mode Decomposition)和子波振幅谱与相位谱关系的时变子波提取方法,根据子波对数振幅谱光滑连续而反射系数对数振幅谱振荡剧烈的特点,采用EMD方法将不同时刻地震记录的对数振幅谱分解为一组具有不同振荡尺度的模态分量,通过滤除振荡剧烈分量、重构光滑连续分量提取时变子波振幅谱;再应用子波振幅谱和相位谱的关系提取时变子波相位谱,将分别提取的振幅谱和相位谱逐点进行合成,最终实现时变子波的准确提取.本文方法不需要求取Q值,适用于变Q值的情况,具有良好的抗噪性能.数值仿真和叠后实际资料处理结果表明,相比传统的分段提取方法,利用本文方法提取的时变子波准确度更高,研究成果对提高地震资料分辨率具有重要意义. 相似文献
12.
谱分解技术在地震解释领域已得到广泛应用,但常用的谱分解方法存在两方面的不足.一是时间分辨率低,难以对薄层进行刻画;二是在烃类检测中多解性强,难以区分流体类型.为了改善该问题,本文提出一种基于地震复谱分解技术的烃类检测方法.复谱分解是指用一个包含多个不同频率Ricker子波的复子波库对地震道进行分解,从而得到时变子波频率和相位信息的过程.借助稀疏反演技术复谱分解可以获得高分辨率的时频能量谱和时频相位谱.本文首先通过拟合算例验证了复谱分解方法刻画薄层的能力以及求取子波频率和相位的准确性.然后利用基于Kelvin-Voigt模型的黏弹波动方程数值模拟对衰减引起子波相位改变的原因进行了分析.最后通过实际资料应用展示了本文方法在储层预测中的高时间分辨率优势,验证了利用子波相位信息识别气藏的有效性. 相似文献
13.
14.
在讨论小波网络理论方法的基础上,研究了利用地震纪录小波多尺度分解属性资料进行虚井声波时差反演的技术方法。分析了利用地震信号进行小波分解和网络学习、训练的理论方法。研究中发现:对于相邻的地震道,较小一段的相似性比整个地震道的相似性要好。据此,利用小波时-频分析技术方法,可以把相邻道的信息外推到其它地震道上。通过以上综合研究及对实际资料进行反演计算、分析,认为小波网络与人工神经网络相比其网络结构要容易选定,并且收敛速度快。同时,利用地震资料分段时-频分析的相似性较好和小波网络学习、训练及记忆能力较强的特点,可以较好地把井旁道的高、低频信息转换到相邻道上。这样在提高分辨率的同时,又增加了反演结果的真实可靠性 相似文献
15.
Xuebing Zhang Cai Liu Xuan Feng Bonan Li Kexin Li Qin You 《Geophysical Prospecting》2020,68(2):371-381
The widely used wavelets in the context of the matching pursuit are mostly focused on the time–frequency attributes of seismic traces. We propose a new type of wavelet basis based on the classic Ricker wavelet, where the quality factor Q is introduced. We develop a new scheme for seismic trace decomposition by applying the multi-channel orthogonal matching pursuit based on the proposed wavelet basis. Compared with the decomposition by the Ricker wavelets, the proposed method could use fewer wavelets to represent the seismic signal with fewer iterations. Besides, the quality factor of the subsurface media could be extracted from the decomposition results, and the seismic attenuation could be compensated expediently. We test the availability of the proposed methods on both synthetic seismic record and field post-stack data. 相似文献
16.
简述了改进的最佳匹配地震子波的小波函数构造及参数的物理含义;根据地震波的高、低频分量在黏弹介质中传播时被地层吸收的差异,给出了一种在时-频域定性估计地震波衰减特性的方法;分别以改进的最佳匹配地震子波的小波及Morlet小波作为母小波分析地层吸收特性,并比较了两种小波函数刻画地层吸收特性的能力;测试了这种方法对噪声的敏感程度. 将文中提出的方法用于某油田的一段实测地震资料衰减分析,得到的吸收特性剖面能较好地反映油气的空间展布. 相似文献
17.
We present a Gaussian packet migration method based on Gabor frame decomposition and asymptotic propagation of Gaussian packets. A Gaussian packet has both Gaussian‐shaped time–frequency localization and space–direction localization. Its evolution can be obtained by ray tracing and dynamic ray tracing. In this paper, we first briefly review the concept of Gaussian packets. After discussing how initial parameters affect the shape of a Gaussian packet, we then propose two Gabor‐frame‐based Gaussian packet decomposition methods that can sparsely and accurately represent seismic data. One method is the dreamlet–Gaussian packet method. Dreamlets are physical wavelets defined on an observation plane and can represent seismic data efficiently in the local time–frequency space–wavenumber domain. After decomposition, dreamlet coefficients can be easily converted to the corresponding Gaussian packet coefficients. The other method is the Gabor‐frame Gaussian beam method. In this method, a local slant stack, which is widely used in Gaussian beam migration, is combined with the Gabor frame decomposition to obtain uniform sampled horizontal slowness for each local frequency. Based on these decomposition methods, we derive a poststack depth migration method through the summation of the backpropagated Gaussian packets and the application of the imaging condition. To demonstrate the Gaussian packet evolution and migration/imaging in complex models, we show several numerical examples. We first use the evolution of a single Gaussian packet in media with different complexities to show the accuracy of Gaussian packet propagation. Then we test the point source responses in smoothed varying velocity models to show the accuracy of Gaussian packet summation. Finally, using poststack synthetic data sets of a four‐layer model and the two‐dimensional SEG/EAGE model, we demonstrate the validity and accuracy of the migration method. Compared with the more accurate but more time‐consuming one‐way wave‐equation‐based migration, such as beamlet migration, the Gaussian packet method proposed in this paper can correctly image the major structures of the complex model, especially in subsalt areas, with much higher efficiency. This shows the application potential of Gaussian packet migration in complicated areas. 相似文献
18.
联合小波变换与偏振分析自动拾取微地震P波到时 总被引：1，自引：0，他引：1
对微地震P波到时的自动拾取是微地震信号分析和数据处理的主要目标之一。基于小波变换的多尺度分析思想,对微地震信号进行小波处理后的小波系数代替原始信号,应用包含在小波变换系数中的信号偏振信息,提出了联合小波变换与偏振分析自动拾取微地震信号P波到时的方法。通过对嘉阳煤矿监测的实际微地震数据进行小波变换,用多尺度小波分解的各个尺度单支重构信号构成协方差矩阵,求解不同尺度协方差矩阵的最大特征值和次大特征值求取P波到时定位函数,实现P波到时的自动拾取,取得了满意的结果. 相似文献
19.
本文利用第二代小波多尺度分解和快速变换的特点,构造自适应计算网格.对初始计算网格上的数值解进行第二代小波变换,得到数值解对应的小波系数空间.小波系数的大小表示相邻网格上数值变化率,小波系数大的区域网格点上的数值解变化梯度大.当小波系数大于等于预设的阈值时,在小波系数对应的网格点周围插入新的计算网格点,通过阈值可以实现网格的细化,得到多尺度下层层嵌套的细化自适应网格;由有限差分法得到相应网格点的空间导数.比较数值算例得到的波场快照和计算时间,验证了该方法的有效性. 相似文献
20.
震源到接收台站之间的地层响应函数能够反映地下介质信息。对地震波传播过程中的卷积模型进行推导:记录信号是众多震源子波经过时移加权叠加的结果;通过反卷积方法可去除震源子波信息,提取震源到接收台站之间的地层响应函数;地层响应函数中第一个突跳值对应的时间即为P波走时。在河北赤城—张北地区进行人工震源实验,通过反卷积计算得到该地区地层响应函数剖面图,得出P波波速约6 km/s。利用人工震源系统还可以对地下介质波速变化进行长期动态监测,对地震预测具有一定意义。 相似文献