简介：通过最新的全球地磁模型——CALS10K.1b,结合CALS3K.4与IGRF11模型,计算并分析了10000BC~1990AD期间中国大陆及邻近地区非偶极子（ND）磁场Z分量的时空变化。为了深入了解ND场的变化,从场源的角度,对2n（n=2-10）极子ND场及其对应的能量进行了分析。结果显示在研究期间ND场的变化可分为3个阶段。在10000BC-2500BC期间,ND场以正值为主并持续了近7500a,在2500BC-1500AD期间强度转弱为以负值为主并持续至1500AD,自此快速增强为以正值为主。东亚地区ND场异常基本在截断阶数（n）为3时即形成,且该异常区已在大陆地区内形成了封闭的圆形区域,这意味着前3阶的ND场占据了总ND场强度的大部分。ND场在核幔边界（CMB）处衰减较快,在地表处则趋于稳定。

简介：全波形反演是一种高精度的地震成像方法，可以对地下介质物性参数模型进行准确的重构。然而在实际应用中，尤其是在三维复杂介质反演中，计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本，提高反演效率。但是使用震源编码技术也带来了两个问题：一方面，参与编码的各个震源之间会产生“串扰噪声”，导致反演结果中出现假象；另一方面，基于震源编码的频率域全波形反演方法周围噪声较为敏感，使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制“串扰噪声”，并基于震源编码技术提出一种频率域自适应全波形反演方法，通过一个与频率相关的自适应选择机制，将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来，在保证反演质量的同时也最大程度地提高了反演效率。

简介：基于Aki-Richards公式和贝叶斯原理,本文发展了利用叠前PP波和PS波资料联合反演P波速度比、S波速度比和密度比的方法。该方法假设参数之间满足正态分布,引入参数协方差矩阵来描述反演参数之间的相关性以提高反演过程的稳定性,并同时使反演的参数序列服从Cauchy分布,引入矩阵Q来描述参数序列的稀疏性以提高反演结果的分辨率。采用本文提出的方法对模型数据和实际多波资料进行反演,结果表明：本文方法正确有效;与传统的单一PP波反演相比,PP波和PS波AVO联合反演具有稳定性更好和反演精度更高等优点。

简介：本文提出了一种新的反演方法:通过采用纵、横波走时数据对(从相同的震源产生的P和S波被同一台站记录)来联合反演纵波速度(Vp)和纵、横波速度比(Vp/Vs),然后单独反演横波速度Vs,在反演过程中同时对地震参数进行定位。该方法不需要假设P和S波的射线路径一致,它是沿着P和S波射线路径计算相对慢度扰动值。该方法直接把Vp/Vs作为一个模型参数,由此能获得比采用从独立反演获得的Vp和Vs计算出Vp/Vs的方法更精确的速度比值。该新方法被应用到反演日本东北地区的壳幔速度及波速比结构的研究中,获得了较好的效果。反演结果表明,在日本东北地区,太平洋俯冲板块为一高Vp,高Vs和低Vp/Vs异常区,而在活火山下方的浅部地幔楔以及背弧深部地区为低Vp,低VS和高Vp/VS异常。虽然这些特征在前人的研究中已经报道过,但与前人的研究结果相比,本次研究所获得的Vp/Vs的空间分布具有较小的分散性,同时,它的分布特征能较好的与地震波速度结构相吻合。

简介：本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。

简介：基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。

简介：在CSEM野外工作中,通常把接收机置于离源3～5倍趋肤深度以外的远区。本文从正演和反演两个角度研究了当源和接收机之间存在导电体时,该导电体对远区中目标体响应的影响。2D有限元正演结果表明导电体主要影响中低频的观测结果,导电体的电导率越低,尺寸越大,影响也越大,并且影响程度随着观测频率的降低而逐渐加大。反演结果表明若不考虑源和勘探区间导电体的存在,反演得到的目标体横向位置会向源的方向移动,带来解释上的误差。在CSEM实际工作中,应该是采用多方位的场源作三维采集和三维反演,以减少源和勘探区间低阻体的影响。更多还原

简介：传统的f-x域经验模态分解法（Empiricalmodedecomposition,EMD）能够有效地对主要由水平同相轴构成的地震记录进行随机噪声衰减。然而,当同相轴倾斜时,f-x域经验模态分解法在衰减随机噪声的同时去除大部分有效信号。本文提出了一种基于f-x域经验模态分解法的改进算法。我们通过局部相似度对所去除的噪声信号中的有效信号进行提取。局部相似度可以用来检测噪声信号中的有效信号点并用来构造一权重算子进行信号提取。新方法与f-x域经验模态分解法、f-x域预测滤波法以及f-x域经验模态分解预测滤波法相比能够在衰减随机噪声的同时保留更多的有用信号。数值模拟实验以及实际地震资料处理结果均表明该方法能更为有效地去噪。

简介：阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以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结构区域内的电场参数估计发挥的重要作用。

简介：

简介：四川盆地及其周缘地区的筇竹寺组和五峰-龙马溪组页岩是目前国内页岩气勘探的主要层位之一,但其地震弹性性质响应规律的区域性特征需要开展相关的实验和理论研究工作予以明确。本研究对干燥状态下的筇竹寺组和五峰-龙马溪组页岩的露头样品进行了超声波速度测试,系统地分析了地震弹性性质随页岩岩石学特征的变化规律。研究结果表明,孔隙度与粘土矿物含量呈正相关、与脆性矿物含量呈负相关;粘土、石英、长石和碳酸盐构成页岩岩石基质,与孔隙共同构成页岩岩石骨架,而干酪根和黄铁矿主要赋存于孔隙中,与页岩骨架的耦合较弱。通过将全部连通性孔隙近似等效于仅存在于骨架粘土矿物之内和采用Gassmann流体替换类似的思路处理干酪根和黄铁矿,可以较为简单地将自相容近似（SCA）理论、微分等效介质模型（DEM）和Gassmann方程组合起来构建研究区页岩的地震岩石物理等效模型。该模型通过关键参数等效孔隙纵横比（采用研究区样品平均值或者由碳酸盐含量进行估算）可以较为准确地预测筇竹寺和五峰-龙马溪组页岩的纵波速度,验证了等效模型的有效性和较广的适用性,可为筇竹寺组和五峰-龙马溪组页岩气储层的测井解释和地震“甜点”预测提供依据。