二维层状介质中衍射波的计算

本文讨论了衍射波的计算,它是文献[5]的继续和发展,是为了应用射线法解决二维地震人机联作工作的一部分。根据波的物理性质,引入了等效炮点、等效检波点和等效界面等概念,将Berryhill导出的衍射波计算公式推广到倾斜界面和楔形结构等情形。这一工作与文献[1—3]中的工作相结合,在实用上解决了二维地震人机联作解释中的正演计算方法问题。     

饱和土中浅埋隧道耦合二维动力响应解答探讨

针对平面波耦合作用下饱和半空间中浅埋隧道二维动力响应的典型问题的解答进行探讨。首先基于Biot饱和孔隙介质理论建立饱和土场方程,并通过中间势函数的代入推导得到饱和土体中的快纵波、慢纵波及横波的Helmholtz势函数方程,随后引入柱坐标系下饱和土体和弹性衬砌的势函数通解,并代入边界条件中,建立以第n组参数为变量的多元一次齐次线性方程组。采用本文的推导思路,可以解决饱和土体中的浅埋隧道在地震作用时P波和S波耦合作用下的响应问题。     

点源二维电场中纵剖面上电位的计算

将总电位ν分解为正常电位．ｕｏ和异常电位ｕ·本文绘出点源二维电场中,异常电位ｕ的傅氏变换Ｕ的变分方程,以及点源二维电场地形影响的异常电位ｕ的傅氏变换Ｕ的积分方程,前者可用有限单元法求解,后者可用边界单元法求解．根据电源点的Ｕ的傅氏反交换,解决了前人尚未解决的直接计算点源二维电场中纵剖面（通过电源点、平行走向的剖面）的电位问题．     

点源二维电场中纵剖面上电位的计算

将总电位ν分解为正常电位．ｕｏ和异常电位ｕ·本文绘出点源二维电场中，异常电位ｕ的傅氏变换Ｕ的变分方程，以及点源二维电场地形影响的异常电位ｕ的傅氏变换Ｕ的积分方程，前者可用有限单元法求解，后者可用边界单元法求解．根据电源点的Ｕ的傅氏反交换，解决了前人尚未解决的直接计算点源二维电场中纵剖面（通过电源点、平行走向的剖面）的电位问题．     

大地电磁Hx型波二维地形改正的方法与效果

用边界无法计算大地电磁Hx型波二维地形影响,并用有限元法计算了起伏地形下水平层状介质的大地电磁响应,用边界元法的结果对畸变后的大地电磁响应进行地形改正,基本恢复了水平地形下层状介质的大地电磁响应的特点．本文给出了物理意义明确、形式简单的视电阻车和阻抗相位地形改正公式．     

二维弹性波的有限元模拟及其初步实践

本文讨论了二维介质中弹性波有限元模拟的一种方法;导出了二维无界空间中集中力点源的理论初动表达式和位移波形表达式,并和相似情况下的有限元结果进行了比较;对用有限元方法算得的几个中小模型的节点位移进行了分析。从所能鉴别的震相的到时、初动符号分布、波形、瑞利波质点运动轨迹等证据看,在得到比较真实的全波理论地震图方面,有限元方法具有很大的潜力。     

二维弹性波的有限元模拟及其初步实践

本文讨论了二维介质中弹性波有限元模拟的一种方法;导出了二维无界空间中集中力点源的理论初动表达式和位移波形表达式,并和相似情况下的有限元结果进行了比较;对用有限元方法算得的几个中小模型的节点位移进行了分析。从所能鉴别的震相的到时、初动符号分布、波形、瑞利波质点运动轨迹等证据看,在得到比较真实的全波理论地震图方面,有限元方法具有很大的潜力。     

流体饱和半空间二维地形三分量弹性波散射间接边界元模拟

无限长局部地形地震波斜入射响应问题称为二维三分量问题,在计算量远小于三维的情况下,一定程度上反映了近地表场地的三维动力响应特征.基于天然土体的成层性及固液两相耦合特性,以层状多孔介质内部移动线荷载(孔隙水压)动力格林函数作为基本解,开展流体饱和半空间二维地形三分量弹性波散射的2.5维间接边界元模拟研究.总场响应由自由波...     

最短路径算法下二维层状介质中多次波追踪

在改进后的最短路径算法（MSPM）中引入分区多步计算技术实现了二维层状起伏介质中的多次透射、反射及转换波波前传播的数值模拟,以及相应的走时和射线路径的跟踪计算.其原理是将二维复杂层状模型按速度界面分成若干个独立的计算区域,采用分步计算技术进行多次波的跟踪计算.基于多次波是通过速度界面的简单入射、透射、反射及转换波按一定规律的不同组合,因此可实施分区多步计算技术.通过某一上、下层界面的透射（或透射转换）波实际上是由上层得到的下行波加上由该界面透射的下行波组成,若为转换波则使用不同的速度模型;而经过某一界面的反射（或反射转换）波实际上是由某层内计算得到的下行（或上行）波再加上由该界面反射的上行（或下行）波组成.这样即可得到分区独立计算,并通过速度界面分步连接达到跟踪多次波的目的.计算结果表明MSPM算法下的分区多步计算技术具有单步SPM算法中的诸多优点,即:算法简单、数值计算稳键、计算精度高、速度快及全球解等,因而是解决多次波跟踪计算行之有效的方法.     

二维复杂层状介质中地震多波走时联合反演成像

采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法.     

层状各向同性多孔介质中井间地震的数值计算

依据Ｂｉｏｔ理论，提出计算理论地震图的相应传播矩阵算法，并针对井间地震震源和检波器的不同配置进行数值计算，结果表明，本文算法高频可达１０００Ｈ２以上．文中算例说明利用震源和检波器相对于储油层的不同位置，可从反射震相记录波形上的差异，判定储油层存在与否．     

层状各向同性多孔介质中井间地震的数值计算

依据Ｂｉｏｔ理论,提出计算理论地震图的相应传播矩阵算法,并针对井间地震震源和检波器的不同配置进行数值计算,结果表明,本文算法高频可达１０００Ｈ２以上．文中算例说明利用震源和检波器相对于储油层的不同位置,可从反射震相记录波形上的差异,判定储油层存在与否．     

井孔声场的本征模及其在分层介质井孔声场计算中的应用

研究无穷大介质中的井孔本征模问题．从有限半径的波导出发,分析半径无限增大的极限情况,证明无穷大介质中的井孔有离散和连续两类本征模,组成正交完备基．进而应用它们分析计算分层介质中井孔的声场,算得孔外分层介质在孔内产生的反射波．     

井孔声场的本征模及其在分层介质井孔声场计算中的应用

研究无穷大介质中的井孔本征模问题．从有限半径的波导出发,分析半径无限增大的极限情况,证明无穷大介质中的井孔有离散和连续两类本征模,组成正交完备基．进而应用它们分析计算分层介质中井孔的声场,算得孔外分层介质在孔内产生的反射波．     

层状介质的声波波动方程反演

基于广义反射透射系数矩阵正演方法 ,讨论了层状介质的声波波动方程反问题 .推导出波数频率域中的雅可比矩阵的解析表达式 ,其计算在正演过程中求出 .采用最小二方法可得到层介质参数 .数值结果表明反演方法的正确有效性 .     

大地电磁二维层状模型参数化反演

本文针对目前二维反演存在求解参数过多和求其偏导数困难的两个问题，在近年来关于连续介质二维大地电磁测深资料快速反演的理论基础上，对二维层状模型反演的目标函数，进行线性近似，采用加权法来近似计算二维模型参数的偏导数。从而实现比较快速的二维层状模型的参数化反演。理论模型的反演试验表明，反演迭代可以较快地以足够的精度拟合视电阻率和相位资料，并收敛接近理论模型。     

各向异性介质中的AVO

分析了横向各向同性和方位各向异性介质的本构关系,由此讨论弹性波在两种各向异性介质中的传播特点,提出可表征这两种介质各向异性程度的广义参数．以此为基础讨论了两种各向异性介质中存在水平界面时的反射系数近似式,将Dely等人推导的横向各向同性介质中的反射系数公式推广到方位各向异性介质的主轴方向上．根据算例讨论修正的Banik和Thomsen的近似式,着重分析两种各向异性介质中的AVO关系及其对实际勘探的影响和指导意义．