论文部分内容阅读
波动互易关系是同一介质或结构的两种波动状态之间的关系,常用于天然地震模拟与震源反演、地球物理测井、地震勘探等地学问题中,也可用于工业和医学的超声无损检测领域。前人建立了几种点源在均匀介质中的互易关系,但对于其在流固分区均匀介质中激发波场的互易关系,以及其他类型点源在均匀介质中激发波场的互易关系,研究得还不充分,这使得互易关系在工程中的应用受到阻碍。本文推导流固分层介质和孔隙介质中几类点源两两之间激发波场的互易关系,并采用数值实验阐释这些互易关系,论述其在波场计算方面的应用意义。地震矩张量是与互易性相关联、在地震学中有重要意义的物理量,但它是针对弹性介质推导的。考虑到震源位于孔隙岩石层,本文通过推导获得了孔隙介质位错辐射波场的等效地震矩张量表达式,计算并分析了波场的特性。论文在以下几个方面取得了进展:给出了声波测井中常用的偶极和四极声源的等效体力,在此基础上获得了流固分层介质中点源激发波场的互易关系,包括集中力、爆炸源、双力偶源、偶极源和四极源。采用三维时域有限差分(FDTD)算法,针对由弹性地层与流体井孔组合的介质模型,编程计算了声源分别位于井内和井外时的波场,通过数值结果阐释了互易关系。结果表明:无论发射源和接收点均位于井内流体中,还是分别处于井内流体或井外地层中,只要发射源和接收信号在类型和指向这两方面满足本文给出的互易条件,原状态和互易状态下的计算波形曲线就重合。数值实验还表明:当原状态与互易状态的发射源和接收器都位于井内时,互易物理量的井外界面微弱反射波的波形也重合。特别地,通过记录单极源激发声场在井轴上的水平位移(或记录水平偶极源激发声场在井轴上的声压),可以有效抑制接收信号中的井孔导波,加强对井外界面反射波的提取效果,这对于当前正在研制的井中远探测成像仪的设计有指导意义。论文通过算例阐明了互易关系在检验数值算法方面的适用性和局限性,指出:应用互易关系可在复杂模型中检验算法和源的施加的正确性,但不能检验界面处理的正确性;应用互易关系,可以简化波场模拟的计算量。基于孔隙介质Biot方程组推导了孔隙介质位错的等效地震矩张量,发现孔隙介质中的位错需要通过两个独立的地震矩张量描述,它们分别作用于孔隙单元整体和孔隙流体上。整体地震矩张量具有传统的弹性介质地震矩张量的意义,而流体地震矩张量则是孔隙介质特有的,引起二者的因素包括固相位错和经断层面向周围岩石的流体注入。对于张开位错,需两个地震矩共同描述;对于剪切位错,流体地震矩为零,剪切位错等效为孔隙单元整体上的双力偶;对于注入源,两个地震矩同为各向同性张量,且对角线元素比值由介质参数决定。给出了孔隙介质地震矩张量辐射位移场表达式,发现流体地震矩仅影响辐射纵波,不影响辐射横波。针对整体和流体地震矩同为单位张量的计算表明,流体地震矩仅对辐射慢纵波有较大影响,因此在计算无限大介质远距离辐射波场时流体地震矩可以忽略,但当计算近距离波场时或计算远距离长时间的流体扩散效应时,流体地震矩不可忽略。进一步的计算表明,对于非均匀孔隙介质,当源距离界面很近时,流体地震矩激发的慢纵波在衰减前就达到界面并转化为快纵波和横波,因此可以对远距离波场产生明显影响。考虑震电效应的计算还表明,流体地震矩对于辐射电磁波有较强影响。得到了孔隙介质中几类点源激发波场的互易关系,包括整体集中力、流体集中力、注入源、整体地震矩、流体地震矩和爆炸源。给出了流体-孔隙介质分层结构中体力形式和点源波场互易关系。采用孔隙地层井孔模型下三维FDTD数值实验阐释了互易关系。通过井内单极源和井外爆炸源的互易实验证明了简单地采用弹性介质互易关系代替孔隙介质互易关系会造成互易物理量波形不重合。在考虑孔隙介质震电效应的情况下,基于Pride震电耦合波控制方程组获得了孔隙介质中震源与电偶极子和磁偶极子激发波场的互易关系。通过格林函数法模拟无限大均匀孔隙介质中集中力和电偶极子的辐射震电耦合波场并阐释互易关系,发现如果互易关系中忽略孔隙介质特有的渗流位移项,将导致互易物理量波形在辐射电磁波部分不重合。采用解析方法模拟水平分层结构中地表集中力与电偶极子激发的震电波场,针对反射波阐释互易关系,互易物理量波形完全重合,这表明在地震勘探中震电和电震方法可以在一定程度上相互替代。计算了地下双力偶源与电偶极子、磁偶极子激发波场,其互易物理量波形完全重合。最后给出了流体-孔隙介质分层结构中震电波场互易关系。这一工作对于震电测井、震电勘探的波场分析、对于地震同震电磁场的模拟具有重要意义。