声波方程频率域有限元参数反演

推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的.     

共散射点道集与角道集串级优化叠前偏移速度分析

推导了基于角度域共成像点道集的叠前深度偏移层析速度分析公式,提出一种共散射点（CSP）道集与角道集串级优化叠前偏移的速度分析方法。该方法通过基于CSP道集的叠前时间偏移速度分析获取初始速度,利用基于角度域共成像点道集（ADCIGs）的叠前深度偏移速度分析进行速度更新。实现步骤概括为：将叠前地震数据映射为CSP道集,利用CSP道集叠加速度谱拾取能量团获取均方根（RMS）速度场;通过Dix公式将RMS速度转换为层速度作为层析的初始输入速度,基于ADCIGs实现叠前深度偏移层析速度反演,最终得到高精度的叠前偏移速度场。断层模型和实际资料试算结果验证了该方法的正确性和有效性。     

起伏地表地震波场角度域照明分析

角度域照明可以更加细致地对地下复杂地质体进行照明分析。提出了一种适用于起伏地表条件下的角度域照明分析方法,其主要过程是利用局部倾斜叠加技术分解波场,对起伏地表处采用加权函数判定局部时窗内的振幅包络位置,以区分其位于地表之上或之下;并引入相似系数函数计算反假频局部平面波场,根据局部斜率与平面波场入射角关系,得到局部角度域地震波场;根据帕斯瓦尔等式下的时间-角度域照明关系式,计算得出角度域照明强度。通过起伏地表高速侵入矿体模型试算,得到给定观测系统下角度域照明强度与地下分析点照明随炮点、波场入射角的变化关系。此方法可用于起伏地表条件的观测系统优化、叠前AVA分析等方面。     

岩土体渗流自由面问题的重心插值无网格方法

岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。     

一种基于模糊正域类的模糊属性约简方法

为了实现模糊信息系统属性约简问题,探讨一种基于正域类的模糊属性约简方法。利用模糊二元相似关系建立模糊相似类的方法,进而定义决策属性关于条件属性的模糊正域类,得到正域类属性的重要度。结果表明:一方面,实现模糊信息系统属性约简;另一方面,通过取不同的相似精度,得到不同属性约简集。     

均匀弹性分层介质模型中的同震位移、应变以及应力

地震发生时产生的同震位移、应变以及应力变化,特别是同震应力的变化,在地震触发等问题的研究中有着重要的意义.本文进一步发展了基于均匀弹性水平层状介质,利用广义反射透射系数矩阵和离散波数计算同震位移的方法,使之可以计算相应的应变、应力以及同震库仑应力变化.可适用于多种情况,接收点可以位于地表以及地表以下,震源类型可以是剪切位错源以及拉张位错源.通过与半无限介质的解析解相比较,结果一致,验证了方法的可靠性以及计算精度,可以用于计算地震之后库仑应力变化,为判断余震分布提供参考.在计算同震位移时,使用了梯形积分与Filon积分相结合的积分方式,即提高了同震位移计算的速度,又保证了计算精度,有利于反演问题研究.     

解释位场全张量数据的张量局部波数法及其与常规局部波数法的比较

全张量探测技术以其信息量大、精度高、干扰小等优点在地球物理领域中得到广泛应用.本文提出采用张量局部波数法来进行位场全张量数据的解释,首先给出了张量局部波数的定义,然后推导出利用张量局部波数法进行反演的基本公式.本文方法在进行张量数据反演时无需事先知道场源体的类型(构造指数)即可获得场源体的位置信息,且可根据位置参数对场源体的类型进行估计.通过理论模型证明张量局部波数法可以很好地完成位场全张量数据的反演工作,并将其与常规局部波数法进行对比,证明全张量局部波数法的反演结果更加准确,即使在测点分布不合理的情况下,张量局部波数法仍可以获得准确的结果.最后应用张量局部波数法对美国得克萨斯州实测重力数据进行了反演,其反演结果与已有的研究成果相一致.     

地震散射波的高精度数值模拟与振幅分析

随着地震勘探要求和技术的不断提高,散射波地震勘探方法逐渐引起人们的重视,尤其在构造复杂地区的高精度勘探中,散射波理论具有广泛的应用前景.本文在讨论了现阶段几种散射波数值模拟方法的利弊后,设计了一种基于微扰论的FK域积分法,在散射场的二次震源和空间能量衰减处理两方面做了改进,其计算精度和效率均得到了提高.此外,采用该算法针对散射波动力学进行了正演模拟试算,从扰动介质空间分布和扰动量两个角度定量地分析了散射波振幅动力学方面的特征.在此基础上得到了五个经验常数,并给出了对应的计算公式和物理意义,通过这五个常数可以划分散射波不同形态的动力学特征.     

Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅰ)

Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum) can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and El Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.     

东北地区古生代地层间断

东北地区古生代在南部有一个地层间断面,北部有三个地层间断面：奥陶纪、晚志留世末-早泥盆世和早石炭世晚期-晚石炭世,相应于古亚洲洋构造域发展的不同阶段,结束于早三叠世末.最后的间断面以后,东北南部和北部古生代地层进入同步发展.