含两种不相混流体的饱和孔隙介质的波场模拟

文中对含两种不相混流体的饱和孔隙介质模型进行了波场模拟,该模型基于封闭型系统假设,考虑了流相与固相的相对运动和孔隙率、饱和度(毛细管压力)的松弛机制,可以较好的模拟波场的数值衰减。与目前常用的等效流体方法和基于开敞系统假设的非饱和孔隙介质模型相比,更符合勘探阶段的实际情况。前人对该介质模型平面波的速度和衰减做了一定的研究,但对整个波场的计算研究还未见报道。本文推导了该模型包含毛细管压力和孔隙度松弛机制的波动方程,并利用有限元的方法进行了波场模拟,并对波场特征进行了分析。数值模拟结果表明,在地震频段,非湿相位移波场中慢波p3较为清晰;毛细管压力(饱和度)和孔隙度的松弛效应对非湿相流体位移有较大影响,随松弛系数的增大,位移减小。     

黏弹性双相介质参数反演的同伦方法

本文以介质的动态时域响应（位移、速度、加速度）作为反演的依据,对黏弹性双相介质的材料参数进行了反分析研究.根据附加条件,并考虑到由于测量误差和测量数据的不充分造成的反演问题的不适定性,将黏弹性双相介质的参数反演问题转化为非线性泛函的极小值问题,然后应用大范围收敛的同伦方法求出使泛函极小的根作为反问题的解.最后,以二维半空间黏弹性双相介质模型为例进行了数值反演分析,数值结果表明本文方法在参数的反演过程中目标函数稳定收敛,且具有一定的抗“噪声”干扰性能;与利用速度和加速度信息反演的结果相比,利用位移信息反演的结果具有较高的精度和抗噪声能力.     

基于经验曲波变换的面波压制方法（英文）

鉴于面波严重影响地震中深层有效反射,而传统的曲波变换面波压制方法无法根据数据实现自适应提取曲波基,影响了去噪效果,本文提出了一种经验曲波变换压制面波的方法。该方法不仅能将地震数据分为多个尺度和方向,而且不同于传统的曲波变换,能根据数据本身的信息,自适应提取曲波基,区分面波与有效波的频谱支集。由于一次有效波与面波能量十分接近,导致面波与有效波曲波系数出现少许重叠,结合奇异值分解的方法,最终在曲波域有效实现了面波与有效信号的信噪分离。模拟数据及实际资料处理结果表明,该方法比传统的曲波变换面波压制方法效果更好,自动化程度更高,是一种先进有效的相对保幅去噪技术。     

地震波形多尺度反演的一点讨论

基于在反演过程中对初始模型依赖性强、易陷入局部极值等问题,本文引入小波分析,提出多尺度地震波形反演方法,从而将参数反演问题转化到小波域中重要系数优化问题。利用多尺度之间的内在联系及小波域中重要系数的稀疏性,有将改进了局部极值、计算量等问题。并对几种多尺度反演策略进行了比较讨论。基于波动方程正演及褶积模型的两种反演方法的数值实便结果显示了本方法良好的效果。     

SH波场中完全匹配层吸收边界研究(英文)

无限域中的波动方程数值模拟往往需要稳定有效的吸收边界来消除人为边界截断所引起的虚假反射.本文首先写出了全匹配层(PML)内SH波的波动方程推导结果,并给出了方程的Crank-Nicolson计算格式与其中窄间导数2阶,6阶,10阶精度的有限差分算法以及伪谱法算法.然后设计了均匀各向同性介质模型和分层溶洞模型并引入图像处理中的信噪比(SNR)概念来定量研究边界吸收效果同PML宽度、不同精度算法的关系.数值结果表明当匹配层宽度比较薄时可以用低精度的有限差格式来获得比较好的吸收效果,当匹配层比较宽时,采用高精度的算法可以获得很好的吸收效果.最后对"反射系数"进行了讨论,指出"反射系数"的不足和文中用SNR来定量衡量吸收边界     

周期成层Patchy模型中纵波的频散和衰减研究

本文用孔隙介质力学的方法来研究纵波在介观周期成层Patchy模型的传播问题.相对于White方法,该方法少了不同介质分界面处的流体压力不等的假设,推导过程更为严密;和Dutta的双相方程解耦方法比较,文中方法直接求解双相方程,形式上更为直观.当纵波通过周期成层孔隙模型时,在低频段用孔隙弹性力学得到纵波频散和衰减结果与用White公式得到的结果符合得很好;周期成层的Patchy孔隙模型由于其空间排列的周期性使其成为孔隙弹性声子晶体,在其频散和衰减曲线上的高频段会出现多个禁带和通带,这使得速度不会随频率单调递增,衰减峰也不只一个.对于地震频段,随着渗透率降低模型的衰减峰向低频移动,这和实际观测结果一致,和传统Biot理论预测结果正相反;随着含气量上升纵波衰减峰值先上升后下降,其最大值出现在含气量0.1左右,这和相关实验结果一致.研究结果表明该模型可以定性解释观测到的地震衰减结果.     

White球状Patchy模型中纵波传播研究

在球坐标系下用直接求解孔隙弹性方程的方法计算了介观尺度下空间周期排列的White球状Patchy模型中纵波传播问题.首先对纵波的衰减和频散进行了计算,并引入了物理学上声子晶体原理来解释高频时纵波在White球状模型中传播的异常现象.在含水饱和度和速度关系的研究中发现,在低频段用等效流体理论和Gassmann理论估计流体Patchy饱和岩石中的纵波速度完全能够满足当前地震勘探的要求.随后的具有相同含气饱和度但有不同周期的Patchy模型研究结果表明,随着空间周期变大,低频的纵波频散变得明显,纵波衰减峰频率向低频移动,但峰值几乎不变.最后,对单元外层含水中心含油的White球状Patchy模型和中心含气White球状Patchy模型进行研究、对比,发现孔隙流体流动对孔隙介质中的纵波频散、衰减影响显著.另外,在具体数值求解过程中用缩减方程组规模的方法解决了线性方程组严重病态得不到正确结果的问题.     

基于学习型超完备字典的地震数据去噪(英文)

基于变换基函数的方法,是地震去噪处理中最常用的技术之一,它利用地震数据在某种基函数变换域内的稀疏性和可分离性来达到剔除噪声的目的。但传统的做法是事先选定一组固定的变换基并在对应域内进行处理,其效果往往并不十分令人满意。为了探索新的改进方法,我们引入学习型超完备冗余字典,即根据地震模型数据进行学习和训练,以寻求最优的稀疏表示字典,而不是只选用固定的变换基。本文在字典学习中融入全变差最小化策略以压制伪吉布斯现象。我们选用离散傅里叶变换作为初始变换,并以随机噪声为例,对单一的全局变换、未经学习的超完备冗余字典和学习型超完备冗余字典做了比较。结果表明,利用经过训练的超完备冗余字典,在对地震数据进行稀疏表示的同时,也达到了有效去除噪声的目的,可视性和信噪比都得到了明显提高。我们也比较了均匀和不均匀字典子块的效果,结果表明,不均匀的字典子块更利于地震数据去噪。     

基于VSP直达波资料的粘弹介质参数遗传算法反演

粘弹性参数变得越来越重要,其反演算法也逐渐成为众多研究者的研究热点。而遗传算法是一种随机、自适应、启发式的算法, 具有很好的鲁棒性和全局收敛性, 本文基于VSP直达波方程,引入了遗传算法来进行粘弹参数反演, 首先将频率域直达波方程表示为复速度的函数,然后通过遗传算法反演出复速度。而复速度和品质因子又是复速度的函数,从而便可很容易的得出。但若直接反演复速度, 反演参数太多, 不容易实现, 所以又将复速度表示成参数C0和C∞的函数,以减少反演参数数量。最后给出了理论模型实验,以证明该算法的有效性。