节点文献

各向异性及双相介质弹性波正演模拟与波场特征分析

Forward Simulation of Elastic Wave and Analysis of Wave Field Characteristics in Anisotropic & Two-phase Medium

【作者】 苏云

【导师】 李录明; 罗省贤;

【作者基本信息】 成都理工大学 , 地球探测与信息技术, 2010, 博士

【摘要】 对地震波传播的研究是地震学、测震学、地震勘探与探测和地震工程学科的基础,因此受到国内外地震学界的广泛重视。一方面地球物理学者们已认识到从地幔到地壳上部的沉积岩都存在不同程度的各向异性;另一方面传统的地震波理论将地下介质近似为单相介质,但储层实际是双相介质或多相介质,当岩石孔隙较大或者孔隙中流体的弹性模量、密度与岩石骨架的弹性模型、密度差值较时,大采用基于单相介质的弹性波理论就会有较大的误差,因此基于各向异性介质及多相介质的地震波理论研究,对于越来越复杂的油气储藏勘探具有重要的理论及实际意义。本文深入研究了各向异性及双相介质中的地震波理论,分别从弹性波方程的建立、弹性波方程的数值求解、相速度及群速度特征方程、弹性波在弹性界面上散射方面来探讨了各向异性及双相介质中的波场特征以及弹性波传播机理。论文首先在分析、研究国内外关于各向异性介质理论、双相介质理论及其相关的正演模拟技术研究成果的基础上,从柯西方程、虎克定律和运动微分方程出发论证了各向异性介质、粘弹性各向异性介质、双相各向同性介质及双相各向异性介质以位移为变量的弹性波方程。其次引入了目前比较流行的交错网格高阶有限差分方法,先对有限差分方法的边界条件、频散、差分格式以及震源加载方式进行了讨论,结合前人的研究成果,选取完全匹配层法处理人工边界、采用通量校正传输方法消除数值频散及选择胀缩震源加载在速度项上,再利用交错网格有限差分方法求解与二阶偏微分方程等价的一阶速度-应力方程,数值模拟各种介质中波场特征,并对不同地质模型的多波波场进行分析,通过上述分析进一步认识了弹性波在各向异性介质和双相介质中的传播规律,这对多波信息的识别具有一定的指导意义。由于弹性波的频散和衰减以及弹性波在分界面处的能量分配规律是地震波传播理论中很重要的问题之一,因此本文从弹性波方程出发,导出相速度和群速度及其特征方程、弹性波在自由界面和弹性界面能量分配方程等。在具体研究中结合各向异性介质的克里斯托夫方程,首次分析了汤姆森参数对相速度随相位角、群速度随群角的影响情况;结合基于BISQ理论的双相介质的克里斯托夫方程,详细地研究了弹性波的相速度、逆品质因子与频率、孔隙度、渗透率、饱和度及波的传播方向之间的关系。联合基于Biot理论的双相各向同性介质的弹性波方程和4个边界条件证明了牟氏方程,讨论了储层顶界面、储层内分界面以及储层底界面处的AVA特征曲线,并与基于单相介质的佐普里兹方程计算结果进行了对比,在此基础上并分析了储层内孔隙度的变化对AVA特征曲线的影响。通过上述分析进一步认识了弹性波衰减特性以及储层参数对弹性波的影响程度,这为多波联合叠前AVA反演岩性参数和储层参数及流体识别奠定了基础。本文以河南油田泌阳地区某口井测井曲线为储层参数,基于Biot理论建立了无耗散情形时的双相层状模型,并采用高阶交错网格有限差分技术模拟了层状介质中的波场记录,分析了孔隙度的改变、孔隙内流体成份差异对共炮点记录的影响,并且对相应的地质模型分界面上的AVA特征曲线及各模型AVA正演记录进行了分析。综上所述,本文通过正演模拟方法分析研究了各向异性介质及双相介质中的弹性波多波波场特征,所取得的成果对于研究地震波传播理论和认识弹性波场特征具有参考意义,对复杂介质的地震勘探与开发、油气横向预测等均具有实际应用价值。

【Abstract】 The worldwide experts in seismology attaches great importance to the research ofseismic wave propagation ,which is the basis for seismology, seismological, seismicprospecting and exploration and earthquake engineering .On the one hand,the scholarsof the geophysics have recognized that the sedimentary rock from the mantle to theupper crust have varying degrees of anisotropy;on the other hand, the traditionaltheory of seismic wave treat the subsurface medium as the single-phase mediumapproximatively,and actually the reservoir is the two-phase or multi-phasemedium.When the porosity of the rock and the elastic modulus and density of thefluids in the pore are large,the seismic wave theory based on single-phase media willhave large error.Thus the theoretical studies of seismic wave based on anisotropicmedia and multi-phase media has important theoretical and practical significance forthe complex exploration of oil and gas reserves.This thesis discusses the seismic wave theory of anisotropic and two-phasemedium in detail.Specifically, this thesis discuss and analysis wave fieldcharacteristics and mechanism of seismic wave propagation in the anisotropic andtwo-phase medium by establishing、solving in the way of numerical method andextending the wave equation.First,this thesis overview the research status quo of anisotropic medium theory、two-phase medium theory and numerical simulation technology. Unifyingpredecessor’s researches, we derived the anisotropic media、viscoelastic anisotropicmedia、two-phase isotropic media and two-phase anisotropic medium’s the secondorder partical differential wave equation with a displacement variable by using theCauchy equations、Hooke’s Law and differential equations of motion. Second thispaper introduces and discuses the staggered grid high-order finite-difference method,including the boundary conditions、dispersion、differential format of the finitedifference method and the way of loading souce.Unifying predecessor’s researches, weuse PML to deal with the artificial boundaries、use FCT to eliminate numericaldispersion and add the source on the velocity item. By solving the first-order velocitystressequations which is equal with the second order partial differential equations using staggered grid finite difference method、simulating the propagation of seismicwave of various mediums and analysising the multi-wave field of different geologicalmodels, we have understood better the law of the wave propagation in anisotropicmedia and two-phase media, which has some signigicances for distinguishing themulti-wave information.This thesis has extended the phase velocity、group velocity and energydistribution equation in the free and flexibility interface from elastic wave equation,because seismic wave dispersion、attenuation and the energy distribution law ofseismic wave propagation are very important issues of the theory of seismic wave.This part of the work include: we have derived Chrisotoffel equations of theanisotropic media and analyzed the phase velocity and group velocity with the changeof the phase and group angle when the Thomson parameter had different values; wehave derived the Chrisotoffel equation of the two-phase medium based on BISQtheory and studied the relationship between the phase velocity of elastic waves、Qfactor and frequency、porosity、permeability、saturation and the direction of wavepropagation; Jointing the wave equation of two-phase isotropic medium based on theBiot thory and four boundary conditions, we have proved the Mou’s equation anddiscussed the AVA characteristic curves of the top reservoir interface、the interfacewithin the reservoir and the bottom reservoir interface, which is compared with theresult calculated using the single-phase medium’s Zoeppritz equation, at the same timewe have analysed the AVA curves changes with the reservoir porosity.Through theanalysis,we further understand the elastic wave attenuation characteristics and thedegree of influence about elastic wave caused by reservoir parameters,which isfoundation for inversing lithology parameters using multi-wave joint prestack AVAand identifying the reservoir parameters and fluid parameters.Finally,this thesis has established layer model without dissipation based on Biotthory using reservoir parameters of a certain well in the region of Miyan of Henanoilfield、have simulated the wave field records using high-order staggered-grid finitedifference technique and have analysed the impact of common shot records caused bythe change of porosity and AVA characteristic curves of each interface.Overall, this paper analysis the wave field characteristics and seismic wavepropagation characteristics of the anisotropic medium and two-phase medium usingforward simulation, the research results have a reference value for studying the seismicwave propagation and understanding the seismic wave propagation,and have practicalapplication value for the seismic exploration and development in the complex media、forecasting the oil and gas,etc.

节点文献中: 

本文链接的文献网络图示:

本文的引文网络