煤岩双相介质弹性波场数值模拟-地球探测与信息技术专业论文.docx

摘 要随着地震工程和能源地震勘探的深入发展，为了较精确的描述地下含流体多孔介质，引入双相介质理论。Biot(1956a、b，1962a、b)奠定了双相介质波动传播理论的基础。在 此理论基础之上对煤岩介质中弹性波场进行数值模拟研究，了解地下煤岩介质的地质构 造特征以及孔隙度大小变化等情况。这对有效地防治煤矿突水、瓦斯突出等自然灾害具 有重要的意义。本文基于双相介质理论建立理论模型，讨论了双相介质的本构关系，通过给出的双 相介质中弹性波的传播方程求解出了弹性波的传播速度，为具体的数值模拟奠定了理论 基础。对数值模拟方法原理有限差分法作了详细阐述和讨论，重点介绍了二维弹性波动 方程有限差分法，对有限差分的实现作了系统研究，对交错网格差分给出了算法流程图。在对弹性波场进行数值模拟时，采用透射吸收边界条件和零相位子波激励源，给出了适 用于二维各向同性介质中的空间四阶、时间二阶有限差分的稳定性条件，用基于位移方 程的二阶中心差分法获得基于Biot理论的二维双相介质中弹性波动方程的离散数值解， 利用Matlab语言编程，分别对几个具有代表性的地层模型进行了弹性波场正演模拟，获 得层状介质模型中固、液双相介质的水平、竖直方向上的质点位移振荡图。正演结果表 明，弹性波在均匀层状煤岩介质中传播时随着孔隙率增大振幅明显变小，在孔隙率减小 的介质中传播时振幅明显变大，弹性波在含有煤岩的层状介质中传播时随着孔隙率的变 化振幅也发生明显的变化，这在弹性波质点位移振荡图中可以直观的看出。以上工作表明，煤岩双相介质弹性波场数值模拟方法对了解地下煤岩介质的地质构 造特征以及孔隙度大小变化等情况具有现实意义和理论价值。也验证了弹性波探测技术 应用于煤矿地质工作中的有效性和可行性。关键词：弹性波双相介质理论有限差分法吸收边界条件正演模拟AbstractAstheearthquakeengineeringandenergyseismicexplorationarefurtherdeveloped,thetheoryofbi-phasemediumisintroducedformore accurate descriptionofthe undergroundporousmediumcontainingfluid.Biot(1956a,b,1962a,b)hasestablishedthefoundationofbi-phasemediumwavepropagationtheory.Basedonthetheory,theelasticwavefieldincoal-rockmediumissimulated,thestatusoftheundergroundcoal-rockmedium’s geologicaltectoniccharacteristicsandthechangesofporositysizeisinvestigated.Ithas importantsignificancetoeffectivelycontrolthewaterandgasoutburstandothernaturaldisastersinmines.Inthispaper,thetheoreticalmodelisestablishedonthebasisofbi-phasemedium,thebi-phasemediumconstitutiverelationshipisdiscussed,andtheelasticwavepropagationvelocityiscalculatedthroughtheequationofelasticwavepropagationinbi-phasemedium,whichlayatheoreticalfoundationforthe specific numerical simulation. Moreover, the principle ofnumericalsimulation,finitedifference,is described and discussed. The finite differencemethodusedin2-delasticwaveequationisemphaticallyintroduced,therealizationofthefinitedifferenceissystematicallyresearched,aswellasthealgorithmofstaggeredgridsdifferenceisgivenout.Inthesimulationofwavefield,thetransmissiveabsorbingboundaryconditionsandexcitingsou