一维垂直各向异性介质频率域海洋可控源电磁资料反演方法

本文提出了一维垂直各向异性(VTI)介质倾斜偶极源频率域海洋可控源电磁(CSEM)资料高斯-牛顿反演方法.在电阻率各向异性介质水平偶极源和垂直偶极源海洋CSEM正演算法的基础上,利用欧拉旋转方法,实现了各向异性介质倾斜偶极源海洋CSEM正演算法.海洋可控源电磁场关于地下介质横向电阻率(ρ_h)和垂向电阻率(ρ_v)的偏导数(即灵敏度矩阵)是解析计算的,结合垂直各向异性介质横向电阻率与垂向电阻率的关系,将各向异性率融入到正则化因子选择中,实现了正则化因子的自适应选择.理论模型合成数据和实测资料反演算例表明,我们提出的反演方法能够较准确的重构海底围岩和基岩的各向异性电阻率以及高阻薄层的埋藏深度、厚度和垂向电阻率.     

青藏高原东南缘基于背景噪声的Rayleigh面波方位各向异性研究

青藏高原东南缘地区是现今地壳形变和地震活动最强烈的地区之一,也是研究青藏高原现今变形机制和构造演化规律的重要区域.本研究使用云南区域地震台网的55个宽频带地震台站连续地震背景噪声数据,采用双台站互相关方法获得Rayleigh(瑞利)面波经验格林函数,提取相速度频散曲线,反演得到云南地区周期5~34s范围内方位各向异性分布图像.反演结果揭示:短周期(5~12s)Rayleigh面波快波优势方向与区域断裂走向有很好的一致性,快波方向随着断裂走向的变化而变化.周期16~26s快波优势方向与反映上地壳特性的5~12s图像总体图像相似,但细节略有不同.其中,滇中块体内易门断裂和滇中块体内东侧的普渡河断裂附近,各向异性快波方向从NS向NW方向旋转;易门断裂以西呈NW向.这反映了青藏高原物质东流和川滇块体受到青藏块体的南东向挤压作用.周期30~34s范围的各向异性,滇缅泰块体和印支块体,快波优势方向为NS和NNW向;而在滇中块体内部,各向异性快波方向呈顺时针旋转变化,可能与青藏高原物质向东逃逸有关.本文还开展了与体波各向异性的对比分析,通过与近震S波分裂、Pms转换波分裂和远震SKS、PKS和SKKS(以后简称为XKS)分裂的对比研究,发现随着周期的增大,得到的快波优势方向与XKS剪切波快波偏振方向趋向一致,与地壳快剪切波偏振方向呈一定夹角.本研究认为,青藏高原东南缘地区壳幔各向异性具有不同的特征和形成机制.     

不同围压和流体饱和状态下致密砂岩弹性各向异性特征

为研究致密砂岩声波速度及其各向异性随围压的变化规律以及不同流体饱和状态下的弹性各向异性特征,钻取了不同方向的岩心并在实验室超声波频率下对致密砂岩的声学特性进行了测量,分别给出干燥和饱和水状态下,不同方向样品纵横波速度、刚性系数以及各向异性系数随围压的变化规律,并对实验结果进行了分析讨论.实验结果表明致密砂岩纵横波速度、纵横波速度比以及刚性系数均随围压增加而增加,但其在不同饱和状态下的变化率却截然不同;纵横波速度比、各向异性系数在饱和水状态下变化规律不明显,表明孔隙流体的存在对于岩石物理性质有着非常重要的影响.这方面的实验工作不但对于考察不同流体性质对致密岩石弹性各向异性影响是必要的,而且有助于致密砂岩油水和气层的识别.     

基于截断牛顿法的VTI介质声波多参数全波形反演

不同类别参数间的相互耦合使多参数地震全波形反演的非线性程度显著增加,地震波速度与各向异性参数取值数量级的巨大差异也会使反演问题的性态变差.合理使用Hessian逆算子可以减弱这两类问题对反演的影响,提高多参数反演的精度,而截断牛顿法是一种可以比较准确地估计Hessian逆算子的优化方法.本文采用截断牛顿法在时间域进行了VTI介质的声波双参数同时反演的研究.不同模型的反演试验表明,在VTI介质声波双参数同时反演中,截断牛顿法比有限内存BFGS(Limited-memory Broyden-Fletcher-Goldfarb-Shanno,L-BFGS)法能更准确地估计Hessian逆算子,进而较好地平衡两类不同参数的同时更新,得到了比较精确的反演结果.     

一种裂缝流体因子的提出及应用

为了解决各向异性下的流体识别问题,将纵波各向异性裂缝预测以及Russell的流体因子融合到直角坐标系中,提出了一种能够同时检测裂缝发育情况以及流体性质的新的裂缝流体因子(Factor of Fluid-filled Fracture,FFF),并通过一组岩性参数检验了裂缝流体因子在裂缝预测及流体识别中的有效性.在理论研究的基础上,选取松辽盆地某地区的火成岩裂缝及流体识别研究为应用实例.通过与测井流体及裂缝信息的对比验证,裂缝流体因子能够较为准确地预测研究区裂缝和流体的分布情况,且裂缝流体因子在单井上的计算结果与单井含气饱和度吻合度较高.此外,根据实际应用效果,指出裂缝流体因子在应用中的局限性:裂缝流体因子在平面成图时受地层厚度影响较大,且无法预测裂缝方向.     

低孔隙度泥页岩应力依赖的各向异性裂纹演化特性研究

本文在实验室对所获取的东营地区层理发育的低孔隙度页岩和泥岩的各向异性裂纹演化特性进行了研究,获得了各向同性条件下泥页岩的力学与超声波响应特性,分析了应力幅度对于页岩声波速度和各向异性的影响.主要结论包括:(1)泥页岩在循环载荷下存在滞后效应,表明其经历了去压实或油气产生导致的超压;(2)泥岩和页岩具有不同程度的各向异性,随着各向同性压力的增高微裂隙逐渐闭合,样品的各向异性程度减弱;(3)分析了岩石韧度和裂纹损伤参数随压力的变化特征,相比泥岩,页岩各向异性程度更高,随压力变化更明显,其裂纹导致的附加各向异性更强;(4)分析了各向异性岩石的动态弹性模量特征,由于软裂隙空间的闭合,动态弹性模量在低压条件下都随着围压的增加有硬化趋势.     

多测线变偏移距VSP地震各向异性反演

本文从一般弱各向异性介质参数反演中得到了使用两条相互正交的变偏VSP测线计算弱各向异性(WA)参数的反演公式. 如果仅仅使用qP波, 则可以确定9个独立的WA参数, 这9个WA参数可以完全地描述井中接收点在两个剖面内介质的各向异性性质. 通过对走时曲线进行最小二乘条件下的三次样条光滑, 可以获得慢度矢量的垂直分量和水平分量. 如果介质是横向非均匀介质, 则水平分量仅仅被用作反演时的约束条件. 为了获得偏振矢量, 本文引入质心计算方法, 该方法计算简单、 稳定, 而且不存在解的奇异问题. 在寻找与一般各向异性介质最接近的高对称性正交各向异性介质和TTI各向异性介质时, 使用qP波各向异性坐标变换方法和最小二乘求解方法, 得到了与一般各向异性介质最接近的正交各向异性和TTI各向异性参数及其对称轴方向参数的计算公式. 使用这些方法, 对瓜哇海地区布设的两条相互正交的变偏VSP测线数据进行各向异性反演, 获得了井中10个接收点处介质的WA参数. 数值计算和实际资料反演表明, 本文所使用的反演方法能够准确地得到VSP井中接收点处介质的WA参数, 这是地震勘探中研究地壳介质各向异性性质最直接和最可靠的方法.      

鄂尔多斯块体及周缘瑞利面波相速度和方位各向异性研究

<正>鄂尔多斯块体被认为是华北克拉通自中生代破坏后残留的稳定地块,内部矿产资源丰富,活动微弱,而周缘地区新生代受到强烈改造,构造运动强烈,地震频发。由于鄂尔多斯块体处于中国东西部过渡区域,受到中国东西部双重动力作用的影响,其深部动力学问题一直是地球科学家研究的热点。     

海南东北部地区地壳地震各向异性特征初步研究

基于海南省地震台网2000~2013年的区域地震波形数据,用剪切波分裂系统分析方法(SAM)获得了海南琼东北部地区"九五"数字台网中2个台站的剪切波分裂参数。结果表明,快剪切波偏振优势方向代表了原地最大主压应力方向。七星岭台NE方向的快剪切波偏振优势方向与区域水平主压应力场方向不一致,与NE走向的断裂一致,体现了局部构造和局部应力场的复杂性;青山岭NNE向的快剪切波偏振优势方向揭示了NNE走向断裂的构造意义。同时,本研究证实,位于活动断裂上或几条活动断裂交汇部位的台站的快剪切波偏振优势方向与对所选用的小地震起控制作用的活动断裂走向一致,而快剪切波偏振优势方向较为离散则反映了该区域复杂的断裂构造和应力分布特征。     

基于各向异性的自适应拟合推估及其在缺失数据拟合中的应用

随机信号协方差函数的拟合和确定是拟合推估的关键。在常规拟合推估方法中, 通常认为随机信号具有各向同性, 而事实上各向异性现象却更为普遍。另外, 在拟合过程中, 很难保证信号协方差与观测噪声协方差的先验方差因子一致。基于此, 本文提出了基于变异函数的各向异性自适应拟合推估方法, 并将其应用于InSAR监测缺失数据的拟合。