基于神经网络的随机地震反演方法

针对随机地震反演中存在的两个主要问题,随机实现含有噪声和难以从大量随机实现中挖掘有效信息,提出了一种基于神经网络的随机地震反演方法.通过对多组随机实现及其正演地震数据的计算,构建了基于序贯高斯模拟的训练集.这也为应用神经网络求解地球物理反问题,提供了一种有效建立训练集的方法.较之传统的神经网络反演,这种训练集不仅保证了学习样本具有多样性,同时还引入了空间相关性.数值模拟结果表明,该方法只需要通过单层前馈神经网络,就可以比较有效的解决一个500个阻抗参数的反演问题.     

粘弹各向异性介质中地震波场模拟与特征

通过引入记忆变量,可以避免粘弹性应力-应变关系中的褶积运算,使波场数值模拟易于实现.通过伪谱法对粘弹各向异性介质中的qP波、qS波数值模拟,结合理论分析,研究了粘弹各向异性介质中速度各向异性和衰减各向异性.衰减各向异性要比速度各向异性更为显著,并且qS波比qP波的衰减各向异性明显.粘弹各向异性介质中,粘弹性对波的影响主要在于波的衰减,各向异性主要影响波前面形状.     

可控力源激发EDA介质中横波分裂现象正演模拟研究

在研究了各向异性介质中波动方程的正演模拟现状的基础上,对波动理论进行了总结,并采用FCT有限差分法对EDA介质中的波动方程进行了数值模拟,得到了EDA介质中横波分裂的波场快照,并对波场特性进行了分析。     

针对致密砂岩气储层复杂孔隙结构的岩石物理模型及其应用

针对致密砂岩气储层复杂的微观孔隙结构进行岩石物理建模,在模型中比较了单一孔隙纵横比、双孔隙模型两种表征孔隙结构的表征方式。岩石物理正演分析表明,两种孔隙结构模型均可解释致密砂岩复杂的速度-孔隙度关系。岩石物理反演结果表明,双孔隙模型对测井横波速度的预测精度更高,说明该模型更适用于表征研究区致密砂岩的孔隙结构,反演的软孔比例参数能够反映地层中孔隙结构的非均匀分布。应用双孔隙模型计算致密砂岩地层岩石骨架的弹性模量,与Krief及Pride等传统经验公式相比,该方法考虑了岩石骨架模量与矿物基质、孔隙度和孔隙结构等微观物性因素的关系,理论上更具有严谨性。对致密砂岩骨架模量计算结果的分析表明,少量微裂隙的存在即能够显著影响致密砂岩骨架的弹性性质,同时孔隙空间中的球形孔隙是致密气的主要赋存空间。并且,通过致密砂岩骨架弹性模量,进一步计算了可用于地层评价的Biot系数等岩石物理参数。致密砂岩骨架模量的预测结果可为Gassmann流体替换理论、BISQ孔隙弹性介质理论等岩石物理方法提供关键参数。     

基于模拟退火粒子群优化算法的裂缝型储层各向异性参数地震反演

发育垂直定向排列裂缝的地下岩石可等效为具有水平对称轴的横向各向同性（horizontal transverse isotropic,HTI）介质。针对HTI介质模型,本文研究了裂缝型储层的各向异性参数地震振幅随方位角变化（amplitude variations with azimuth,AVAZ）的反演方法。首先,在地震AVAZ反演流程中,提出采用模拟退火粒子群优化算法实现裂缝型储层各向异性参数反演。之后,通过理论模型测试,验证了基于模拟退火粒子群优化算法的地震AVAZ反演的有效性。最后,将反演方法应用于四川盆地龙马溪组页岩气储层实际方位地震数据;在反演之前先利用傅里叶级数方法估计裂缝方位并对实际数据进行方位校正,以提供更准确的输入数据;通过计算得到的P波、S波各向异性参数可用于评价裂缝发育程度。反演结果表明,研究区域构造顶部裂缝较发育,与地质基本理论一致,验证了反演方法的合理性。     

遗传算法和LM算法联合反演瑞雷波相速度

上世纪中叶以来,微震探查法等基于面波相速度频散测定以及其反演的地下构造勘探方法,在地壳研究、工程地质勘探等方面得到了广泛的应用.本研究分析了遗传算法（GA）应用在面波频散反演当中出现的问题,采用遗传算法和Levenberg-Marquardt 算法（LM法）的联合运用法——以遗传算法反演出来的大局的搜索结果为初始点,再进行LM法搜索,由此克服两者的缺点,实现高精度的反演.     

双相介质中地震波的频率-空间域数值模拟

本文利用优化的25点频率-空间域有限差分算法对基于BISQ模型双相各向同性介质中的地震波进行了数值模拟.通过与经典的Biot模型理论模拟结果进行对比,分析了Biot流动（宏观流体流动）和Squirt流动（微观流体流动）耦合作用对地震波在孔隙介质中传播特性的影响.数值模拟在地震频段进行,结果显示：在理想相界和黏滞相界情况下,Squirt流动机制都比Biot流动机制产生了更大的速度频散和能量衰减.其中,在Biot流动和Squirt流动耦合作用下的快P波的速度和振幅小于仅考虑Biot流动影响下快P波速度和振幅,而且慢P波的衰减也更加强烈.本文还研究了地震波在双层双相各向同性介质分界面处的反射和透射特征,双相介质中波的反射与透射现象类似于单相介质的情况.模拟结果表明,利用优化25点频率-空间域有限差分法模拟双相孔隙介质中的地震波场是可行的,这为开展双相孔隙介质全波形反演问题的研究提供了可能.     

交错网格任意阶导数有限差分格式及差分系数推导

交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。     

非弹性层状介质地震波频变AVO响应模拟及分析

以非弹性层状介质为模型,基于广义传播矩阵理论计算地震波频变反射系数,算法中同时考虑了与频散和衰减有关的地层岩性因素,以及与薄互层有关的地层结构因素.实现了岩石物理模型、反射系数这两个计算过程的"无缝"连接,精确考虑了由复数弹性模量表示的地层非弹性因素,也为在同一反射模型中考虑源于不同物理机制的频散与衰减提供了方法.数值模拟结果验证了算法的有效性和稳定性,计算结果表明,非弹性薄层的反射振幅随频率先增加后减小,不一定表现常规"低频亮点"异常;同时,薄互层条件下的频散与衰减使得地震反射波的频谱以及AVO特征呈现与频率相关的复杂变化.本文完善了频变AVO算法,为含油气储层频变AVO响应的模拟和分析提供了方法.