简介:本文通过地震正演模型和面波振幅数据对比分析了天然地震、爆炸,与核试验数据,以探索改进辨别天然地震与爆炸的方法。所提方法是基于波数域和频率域的双积分变换解。用于研究的地震事件都发生在北朝鲜,其中包括2001年6月26日(39.212°N,125.383°E)和2001年10月30日(38.748°N,125.267°E)所记录的爆炸数据,2006年10月9日(41.275°N,129.095°E)所记录的一次核试验数据、以及2002年4月14日(39.207°N,125.686°E)和2002年6月7日(38.703°N,125.638°E)所记录的两次天然地震。基于这些数据通过地震波分析与理论模拟表明爆炸引起的地震与天然地震有着不同的波型特征。爆炸引起的信号特征是P波的能量较S波强,在爆炸记录上0.05-0.5Hz频率之间Rg波清晰呈现,而在然地震记录上没有。这是由于爆炸地震记录上P波是优势波,与SH成分发生了耦合。
简介:本文通过人工变换T2分布和建立管-球模型模拟法研究含水合物地层渗透率与水合物含量之间的关系。首先,在渗透率的模拟试验中,我们改变了束缚水与可动水的比例、总孔隙度以及与之关联的T2分布。试验结果表明,相对渗透率与水合物含量之间的关系受到这些因素的制约。随后,我们用管-球模型表示水合物生长的孔隙空间,并把水合物的生长过程看成是向孔隙空间随机扔小球的过程。在此过程中,采用两种方法计算渗透率,一是Schlumberger'sT2公式(即SDR模型),二是Darcy定律与Poiseuille流动方程相结合的方法。前人的实验研究表明,在一定的水合物含量范围内,渗透率基本保持不变。以此为参考,我们将计算结果与之进行比较。我们发现,采用SDR模型时,渗透率的数值模拟曲线与Masuda模型N=15时的结果相近。而采用Darcy定律时,渗透率模拟值较高,但与实验结果的趋势相一致,都会出现渗透率的平直阶段。尤其,当水合物晶体在孔隙体内优先生成时,优先的概率越高,渗透率的平直范围越大。
简介:针对祁连山冻土区DK-4井孔含水合物岩层构建岩石物理模型,分别采用K-T方程模型方法和区分填充模式的等效介质模型方法(模式I和模式II)。K-T方程主要模拟地震波在两相介质中传播,基于弹性模量计算速度;而等效介质模型主要依据对水合物地层介质的两种假设:模式I是将水合物作为孔隙填充物的一部分,模式II是将水合物作为岩石骨架的一部分。首先根据粉砂岩层段的测井数据提取了水合物地层骨架的物性参数,包括纵波速度、横波速度、密度、体积模量和剪切模量。然后依据水合物地层各主要成分的物性参数,建立了基于K-T方程的岩石物理模型和区分填充模式的等效介质模型。将两类模型的速度曲线分别与实际地层数据进行了对比:由K-T方程建立的岩石物理模型,其理论计算的速度偏离实际值;而区分填充模式的等效介质模型,其理论计算的速度符合实际值,并且填充模式II模型的速度曲线比填充模式I模型更接近实际地层情况。采用区分填充模式的等效介质模型,能够更好地实现对祁连山冻土区水合物粉砂岩地层的模拟。
简介:根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震(OBS)共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。
简介:天然气水合物沉积层在常规叠后地震剖面上,一般表现为似海底反射(BSR)特征,常作为水合物识别的重要标志,但由于地震数据分辨率的限制与多解性的存在,类似的反射特征不一定是水合物的表现。本文先从南中国海神狐海域地质条件分析人手,综述研究区域天然气水合物发育的构造、沉积与运移环境;其次以加权方程与Zoeppritz公式的正演模拟为基础,讨论天然气水合物与BSR特征的关系;然后选取经过三口已钻获天然气水合物站位的地震测线,提取角度域共成像道集用于BSR的AVA特征研究;最后依据叠前地震弹性参数同步反演思路,对水合物重点赋存区域开展精细预测研究,获取了研究区域含水合物沉积层稳定带的分布特征与饱和度数。该方法可作为采用多道地震数据估算天然气水合物含量的一种有效技术。
简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:阵列声波信号是典型的非线性、非平稳信号,Hilbert~Huang变换(HHT)是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解(EMD)和对瞬时频率的求解,可以获得声波信号的时一频谱。其关键技术就是进行经验模态分解,任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征,在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程,并列举了一些声波测井波列实例,说明了该方法的有效性。
简介:阵列声波信号是典型的非线性、非平稳信号,Hilbert-Huang变换(HHT)是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解(EMD)和对瞬时频率的求解,可以获得声波信号的时-频谱。其关键技术就是进行经验模态分解,任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征,在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程,并列举了一些声波测井波列实例,说明了该方法的有效性。
简介:阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期(超过100秒),逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。
简介:随着全张量重力梯度(FTG)测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU(GraphicsProcessingUnit)和预处理分解技术(Preconditioningmethods)在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛(SSOR)技术与不完全乔列斯基分解共轭梯度算法(ICCG)提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。