热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函数法数值模拟,详细研究两个热耗散衰减波的频散和衰减特征,着重分析热导率、热膨胀系数及比热的变化对波速和衰减的影响.研究表明热导率作为主要参数决定了波速与衰减的临界变化,热膨胀系数对波速和衰减的幅度有明显影响,比热则兼顾了前两个热弹系数的影响特征.最后,利用热弹性动力学频率域的二阶格林函数进行波场快照数值模拟,展示热弹性介质中纵波、横波和热波的传播行为.     

基于改进反射系数近似方程的纵横波阻抗同步反演

基于地震波反射系数近似公式的叠前反演是油气勘探的重要工具.本文在已有研究的基础上,推导了一个改进的射线参数域地震纵波反射系数近似方程.该方程建立了地震纵波反射系数与纵波阻抗和横波阻抗的非线性关系,在中、小角度的范围内较现有的反射系数线性近似公式精度更高.另外,由于该方程仅包含纵波和横波阻抗反射系数项,因此基于新方程的反演能够有效地降低同步反演纵波速度、横波速度、密度三个参数的不适定性.在此基础上,结合广义线性反演法(GLI)理论和贝叶斯理论,相应地发展了一种叠前地震同步反演方法.模型测试和实际资料的应用表明,基于新方程的反演方法能够利用有限角度(偏移距)的数据稳定地反演纵波和横波阻抗,由于在反演过程中,不需要假设纵横波速度为常数,因此该方法还能有效地提高反演结果的精度.     

双相各向异性介质中弹性波传播特征研究

随着地震工程和能源地震勘探的深入发展,人们所遇到的地下介质愈来愈复杂．常规的各向异性介质理论或双相各向同性介质理论难以精确描述含流体的各向异性介质,如裂缝性气藏、含水页岩等．本文以Biot双相各向异性介质理论为基础,利用弹性平面波方程,推导出了任意双相各向异性介质中弹性波的Christoffel方程．根据Christoffel方程,计算并分析了频率对双相横向各向同性介质中弹性波的相速度、衰减、双相振幅比和偏振特征的影响．结果表明,在4类波(快纵波、慢纵波、快横波和慢横波)中,频率对慢纵波影响最大;当耗散很大时,快纵波、快横波和慢横波的流固相振幅比值近似为1．对偏振特征分析的结果表明,在双相各向异性介质中,弹性波的固相位移偏振方向与流相位移偏振方向将不再保持同向或反向,而是呈不同大小的夹角．      

各向同性介质弹性阻抗的张量表示

作为声波阻抗的广义形式,当前地震勘探领域所应用的弹性阻抗描述了在地震纵波非垂直入射的情况下地层对反射纵波的阻抗特征,因此能够反映地层的横波速度信息.但是,根据其理论公式的推导,传统的弹性阻抗实质是反射波的一个衍生标量属性,且存在诸如弱阻抗差、小入射角和相邻地层纵/横波速度比恒定等假设条件.因此在相应的反射系数重构、地震反演与解释过程中不可避免地存在误差.通过分析各向同性弹性介质中速度场与应力场之间的关系,给出了弹性阻抗张量的精确公式表征.弹性阻抗张量能够表示每个应力分量与速度分量间的关系,因此完备地描述了介质的弹性力学特性.文章通过理论模型和实际数据应用证明了导出的弹性阻抗张量表达式较经典的弹性阻抗公式有更高的精度和适用范围.     

基于散度和旋度纵横波分离方法的改进

纵、横波的分离是多波多分量地震资料处理中很重要的一步,其分离结果直接影响到后续数据处理的质量.各向同性介质中纵波为无旋场,横波为无散场,因此可以在频率-波数域利用散度和旋度算子对地震记录进行纵、横波分离,但是此处理过程必须知道地表处的纵、横波速度.本文给出了一种估算地表纵、横波速度的方法,可以在纵、横波速度值未知的情况下,将其估算出来.针对弹性波场进行散度和旋度运算时,纵、横波的相位和振幅比发生改变的问题,本文给出了相位和纵、横波振幅比的校正方法.     

基于Low rank近似的黏弹性介质衰减补偿微地震逆时定位与裂缝成像

真实地下介质具有黏弹性,地震波在传播过程中会发生耗散与频散.忽视黏弹性介质的吸收衰减效应,逆时延拓过程中地震波将会出现振幅减弱、相位失真等现象,无法准确定位震源真实位置,因此需要对黏弹性介质中传播的波场进行衰减补偿,并通过采用合适的成像算子对微地震震源进行定位与裂缝成像.本文基于耗散与频散解耦的分数阶黏弹性波动方程模拟波场,采用low rank分解近似混合域算子,分离衰减相关项并反转耗散项符号,并在补偿的衰减项波场的波数域中进行低通滤波,压制噪声的影响;使用优化后的成像算子进行微地震震源定位,并通过分离散射波场,对散射波进行逆时反传寻找裂缝.数值实验证明,本文方法通过low rank近似有效提高了计算效率,衰减补偿算子在滤波器约束下能够稳定地补偿反向延拓的波场,优化后的成像算子能够在压制随机噪声的同时进一步提高计算效率和定位分辨率.     

层状弹性介质的波阻抗反演

讨论层状弹性介质的纵波阻抗和横波阻抗的同时反演问题，由纵波线源激发产生的弹性波场，通过倾斜叠加可以得到其平面波分解，分解出来的平面波场对应于以不同倾角斜入射激发所产生的波场，它们满足由二维弹性波基本方程派生出的一维偏微分方程组。利用在地面同时接收的位移速度的水平和垂直分量，可以由每个特定角度斜入射产生的平面波场所满足的一维方程组反演层状弹性介质的两个参数，在已知密度情况下，可以反演纵波阻抗和横波阻     

横向各向同性多孔介质中的地震波传播

基于各向异性多孔介质中的广义Ｂｉｏｔ理论，导出了横向各向同性多孔介质中波传播的特征方程．指出在多孔介质中有４种类型的频散和耗散波传播：准纵波ＱＰ１（快纵波）、准纵波ＱＰ２（慢纵波）、准横波ＱＳＶ和横渡ＳＨ．文中给出了４种波速度的解析表达式．数值计算频率曲线和衰减曲线与Ｓｃｈｍｉｔｔ（１９８９）用均值处理得到的结果类似．还给出了波传播过程中３种类型准体波之间的耦合系数（或称转换系数）．     

叠前纵波和转换波地震资料Q值提取及反Q滤波(英文)

纵波和转换波联合的多波地震勘探技术是解决复杂油气勘探的有效技术,提高转换波的分辨率是其关键问题之一。影响转换波分辨率的主要原因是地层对地震波的吸收,有效地计算地层Q值、消除地层吸收对转换波传播的影响,是提高转换波分辨率的关键。本文提出了从叠前转换波道集中估算横波Q值的方法,并利用沿射线路径的波场延拓,将一种稳定有效的反Q滤波方法应用到叠前共炮点纵波和转换波道集的衰减补偿中。模型资料结果表明,本文提出的估算转换横波Q值的方法精度较高;模型资料和实际资料的吸收补偿结果表明,此稳定全反Q滤波能有效地提高叠前纵波和转换波资料的分辨率。     

层析成象中的数学方法问题及其应用进展(Ⅱ)

三、基于波动方程反演的散射层析成象方法地震波的传播是服从弹性波方程的,但是要严格的从弹性波方程反演来研究散射层析成象是远远超出当前数学发展水平的,其初步探讨可见[13]。熟知弹性介质中存在两种波:压缩波(纵波)与剪切波(横波),如介质参数的变化率较小,则压缩波场u近似地满足变系数的声波方程如下     

怀来、 延庆及其周边地区地壳 介质各向异性研究

提出了利用人工爆破P波走时反演地壳介质方位各向异性参数的方法. 在假定介质是弱各向异性介质的情况下, 使用扰动理论得到了线性化的反演公式, 其中待反演的弱各向异性参数是P波走时的线性函数. 如果在反演公式中参考走时取相同震中距接收点的P波平均走时, 那么所获得的弱各向异性参数与参考介质速度的选取无关. 反演得到的弱各向异性参数可以看作是不同震中距和不同深度范围内介质的等效弱各向异性参数. 等效弱各向异性参数在一定程度上反映了不同深度范围内水平方向相速度随方位的变化. 这种变化可能是不同时期构造应力作用的结果. 2007年中国地震局在首都圈怀来地区实施了一次大吨位人工爆破实验, 以爆破点为中心, 布设了高密度的地震观测台网和台阵. 台站相对于爆破点具有360°的全方位覆盖, 所得到的地震记录数据为研究怀来、 延庆地区地壳介质P波方位各向异性提供了必要条件. 我们通过走时反演获得了与水平方位相关的弱各向异性参数, 并对弱各向异性参数进行坐标变换, 得到了能够直观描述岩石弱各向异性的具有水平对称轴的横向各向同性介质, 给出了对应的3个独立弱各向异性参数及其对称轴方位, 讨论了介质各向异性与构造应力场的关系. 结果表明该地区地壳介质存在明显的方位各向异性, 其最大值约为4.6%.      

各向异性介质qP波传播描述I：伪纯模式波动方程

地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征．各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用．为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程．本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程．文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势．     

反卷积在人工震源数据处理中的应用

震源到接收台站之间的地层响应函数能够反映地下介质信息。对地震波传播过程中的卷积模型进行推导:记录信号是众多震源子波经过时移加权叠加的结果;通过反卷积方法可去除震源子波信息,提取震源到接收台站之间的地层响应函数;地层响应函数中第一个突跳值对应的时间即为P波走时。在河北赤城—张北地区进行人工震源实验,通过反卷积计算得到该地区地层响应函数剖面图,得出P波波速约6 km/s。利用人工震源系统还可以对地下介质波速变化进行长期动态监测,对地震预测具有一定意义。     

山西临汾流动地震观测台址P／S波速度结构测试

用单孔检层法测定山西临汾流动地震观测台址两个钻孔的P／S波速度结构,为研究地震观测台址速度结构对地震波接收效果的影响提供基础资料。与常规P／S波测井不同的是,在P波速度测试中采用高频24道水听器串共炮点接收,有效降低单点多炮接收因重复性问题带来的零时误差,提高P波速度测定的精度。该方法可为类似工程场址的高精度P／S波速度测试提供新途径。     

含裂隙双相各向异性介质中的地震波传播

基于各向异性介质和含裂隙双相介质中的地震波传播理论,导出了具有四阶对称轴的含裂隙双相各向异性介质的本构关系与波传播的运动方程．指出在含裂隙双相各向异性介质中可能传播６种类型的准纵横波,即快纵波ＱＰ１,慢纵波ＱＰ２,两个分裂的快横波ＱＳＶ１,ＱＳＨ１和两个分裂的慢横波ＱＳＶ２,ＱＳＨ２．并以平面波传播为例作了进一步分析     

Propagation characteristics of converted refracted wave and its application in static correction of converted wave

Three-component seismic exploration through P-wave source and three-component geophone is an effective technique used in complicated reservoir exploration. In three-component seismic exploration data processing,one of the difficulties is static correction of converted wave. This paper analyzes propagation characteristics of non-converted and converted refracted waves,and discovers a favor-able condition for the formation of converted refracted wave,i.e. the velocity of overlaying medium S wave is much lower than that of underlying medium S wave. In addition,the paper proposes the static correction method of converted wave based on PPS converted refracted wave,and processes the real three-component seismic data with better results of static correction of converted wave.     

Natural gas origins of large and medium-scale gas fields in China sedimentary basins

The integrated study of the geological and seismic reflection data from the drilling area of CCSD has discovered that the density and the P-wave velocity of orthogneiss are almost the same as that of the paragneiss in the area; but the orthogneiss and the paragneiss hold different reflection attributes. The strong seismic reflector packes coinciding spatially with the paragneiss suites have implied that the paragneiss buried in the metamorphic crust itself can cause bone-like seismic reflector sets. The P-wave velocity of paragneiss shows little apparent difference with that of the orthogneiss; but its transverse wave velocity is lower, with the Vp/Vs ratios being high. The paragneiss has partially inherited the layering structures and textures of the protolithe of sedimentary rocks, hence shows strong heterogeneity and anisotropy, that is why the paragneiss are able to produce the bone-like reflectors in the upper crust. The low transverse wave velocity of paragneiss often means weak shear resistance, which will further cause cracks or fractures in the rock, consequentially increase its porosity and permeability during tectonic movements, and form the paragneiss reservoirs of low-permeability zones for gases uplifted from the deeper crust. Because the paragneiss in the crustal metamorphic basement can cause the seismic reflectors, seismic reflection sections are able to provide information about the paragneiss under certain prerequisites.     

黏声介质平面波有限差分叠前深度偏移及在莺歌海盆地的应用

莺琼盆地诸探区中存在底劈现象,深层气源产生的气体沿底劈产生的裂隙通道向上漫溢.漫溢过程中,一是充填在遇到的砂体中形成气藏;二是弥漫在上溢通道中,使得通道中的纵波速度发生变化,进而纵波波阻抗差异变小,反射变弱.另外,通道中气体的存在,会加强地层的吸收衰减,使得地震波振幅变弱、高频成分损失导致同相轴分辨率降低.利用OBC数据进行多波地震勘探和利用黏声介质的叠前深度偏移都是改善模糊区成像质量的重要方法技术.为此,本文提出用黏声介质平面波有限差分法叠前深度偏移成像方法改善气体充填区域的成像质量.黏声介质成像目的是补偿地震波的吸收衰减;平面波偏移成像目的是适应海上单炮数量巨大,提高波动方程叠前深度偏移成像的效率;有限差分法叠前偏移的目的是适应该区浅层气分布局域性极强、Q值的空间变化大的情况.在莺歌海某探区的实际数据上的黏声介质平面波有限差分叠前深度偏移试验证明,本方法是改善模糊区成像质量的较为有效的途径.     

二维弹性随机介质中的波场特征

通过波动方程的交错网格有限差分正演,模拟了地震波在二维弹性随机介质中的传播及其自激自收时间记录.本文将理论记录剖面分割成三个不同的时间区段,在三个时间区段上分别计算并提取剖面的统计特征.运用正交试验设计和分析方法,本文全面研究了混合型弹性随机介质的五个统计特征（自相关长度、粗糙度因子、扰动标准差、纵波背景速度、纵横波速度比）及其交互作用与各种波场特征之间的关联,并得出了若干结论.