Stability and Grid Dispersion of the P-SV 4th-Order Staggered-Grid Finite-Difference Schemes

Stability and grid dispersion in the P-SV 4th-order in space, 2nd-order in time, displacement-stress staggered-grid finite-difference scheme is investigated in the case of a homogeneous unbounded medium. All results, however, also apply to the velocity-stress and displacement- velocity-stress finite-difference schemes. Independent stability conditions for the P and S waves are obtained by exact separation of equations for the two types of waves. Since the S-wave group velocity can differ from the actual velocity as much as 5% for the sampling ratio 1/5, commonly used in numerical modelling, the sampling of the minimum S wavelength by 6 grid spacings (with the velocity difference not larger than 2.5%) is recommended. Grid dispersion is strongest for a wave propagating in a direction of a coordinate axis and weakest for a wave propagating along a plane diagonal. Grid dispersion in the 4 ^th -order scheme for the sampling ratios s = 1/5 and s = 1/6 is smaller than grid dispersion in the 2 ^nd -order scheme for s = 1/10 and s = 1/12, respectively.     

地震波有限差分模拟综述

本文从有限差分法数值模拟技术的各个方面对地震波有限差分模拟的发展和现状进行了论述.波场的数值模拟技术是认识地震波传播规律,检验各种处理方法正确性的重要工具,地震波的数值模拟是地震波传播规律研究的必要手段,贯穿于地震资料的采集、处理、解释的整个过程中.有限差分法数值模拟技术相对于射线方法具有更高的精度,同时比有限元方法计算量小,因此在实际应用中占很重要的地位.     

高精度频率域弹性波方程有限差分方法及波场模拟

有限差分方法是波场数值模拟的一个重要方法,但常规的有限差分法本身存在着数值频散问题,会降低波场模拟的精度与分辨率,为了克服常规差分算子的数值频散,本文采用25点优化差分算子,再根据最优化理论求取的优化系数,建立了频率空间域中弹性波波动方程的差分格式;为了消除边界反射,引入最佳匹配层,构造了各向同性介质中弹性波方程在不同边界和角点处的边界条件. 最后由弹性波波动方程和边界条件,通过频率域有限差分法,分别利用不同震源对弹性波在均匀各向同性介质、层状介质及凹陷模型中的传播过程进行了数值正演模拟,得到了单频波波场、时间切片和共炮点道集,为下一步的研究工作（如成像、反演）提供了研究基础.     

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E-polarization induction problems

Introduction The algorithm of the finite difference method (FDM) for solving the 3-D induction problem used by CHEN (1985) is based on the concept suggested by Brewitt-Taylor and Weaver (1976), which assumes the conductivity to be a smoothly-varying function of position. After the comparisons between many model calculations carried out by the author and other methods showed that the vertical components estimated by our method are acceptable in general cases. However, the reliability of the…     

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E -polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.     

Tube wave modeling by the finite-difference method with varying grid spacing

A finite-difference approach of aP-SV modeling scheme is applied to compute seismic wave propagation in heterogeneous isotropic media, including fluid-filled boreholes. The discrete formulation of the equation of motion requires the definition of the material parameters at the grid points of the numerical mesh. The grid spacing is chosen as coarse as possible with respect to the accurate representation of the shortest wavelength. If we assume frequencies lower than 250 Hz then the grid spacing is usually chosen in the range of a few meters. One encounters difficulties because of the large-scale difference between the grid spacing and the size of the borehole, usually several centimeters.These difficulties can be overcome by a grid refinement technique. This technique provides the construction of grids with varying grid spacing. The grid spacing in the vicinity of the borehole is chosen such that the borehole is properly represented. An example demonstrates the accuracy of this technique by comparisons with other methods. Unlike many analytical methods, the FD method can handle complex subsurface geometries. Further numerical examples of walk-awayVSP configurations show tube wave propagation within fluid-filled boreholes of realistic diameters.     

Poroelastic finite-difference modeling for ultrasonic waves in digital porous cores

Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with small-scale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference (FD) method of Biot’s poroelastic equations, with an arbitrary even-order (2L) accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer (CPML) absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultrasonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thicknesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Q _sc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments.     

完全匹配层在数值模拟孔隙介质声波传播中的应用

The nonsplitting perfectly matched layer （NPML） absorbing boundary condition （ABC） was first provided by Wang and Tang （2003） for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot＇s equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber （DW） method.     

完全匹配层吸收边界在孔隙介质弹性波模拟中的应用

模拟弹性波在孔隙介质中传播,需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上,首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先,本文通过均匀孔隙模型的数值解与解析解的对比,验证所提出的数值方法的正确性;然后,本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能,将完全匹配层与廖氏和阻尼吸收边界进行了对比,研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明,在孔隙介质中,完全匹配层作为吸收边界能十分有效地吸收衰减外行波,无论对体波还是面波,是一种高效边界吸收算法.     

频率域波动方程的双参数识别方法

利用反射波地震记录恢复波动方程的两个系数以研究地球内部构造，是当前应用地球物理研究的热点之一，本文在波动方程反演方法的理论基础上，导出了频率域识别介质的密度波速两个重要参数的新算法，从而把当前的波动方程单参数反演双参数反演方法，通过理论模型及实际测井资料的反演试算对比，均藜得良好效果，在算法中还省去了Ｇｒｅｅｎ函数的繁杂计算，节省了计算量，增加实用性。     

弹性波场数值模拟的隐式差分多重网格算法

为了得到稳定的弹 性波数值模拟,而不得不选择隐式差分方程;为了提高解的精度,又不得不增加节点数目, 但同时也降低了隐式迭代求解的收敛速度. 为此,本文使用隐式差分的多重网格算法进行弹 性波数值模拟,多重网格算法通过粗网格收敛较快的迭代过程求出近似解,以近似解为初值 使用细网格进行精确的迭代求解,从而加速了隐式迭代求解的过程,能够以较高计算速度、 精度、稳定性完成弹性波传播过程的数值模拟.     

弹性波数值模拟的延迟边界方法

在地震波场的波动方程数值模拟中,由于计算量的限制,必须加入人为的边界,使模拟计算可以在一定的空间范围内进行. 由于边界节点上的波场值不能像模拟区域内部的节点一样使用中心差分来计算,使其计算精度大大降低,从而产生边界反射. 为了消除边界反射,本文提出了延迟边界方法,根据弹性波在传播方向上等距离质点的等相位延迟性质和振幅衰减特性,由内部波场的时空分布,推算出边界波场的相位延迟的大小和振幅衰减系数,从而提高边界节点上的波场值计算精度,消除边界反射的产生.     

VTI介质准P波旋转交错有限差分数值模拟

本文采用旋转交错网格差分格式对VTI(垂直对称轴的横向各向同性)介质准P波一阶应力-速度方程进行数值模拟。并在PML边界条件和稳定性条件下得出Marmousi等复杂模型的高精度波场快照和地震记录,分析了各向异性对地震波的影响。数值结果表明:旋转交错网格有限差分能获得高精度的地震模拟数据,PML边界有较好的吸收效果。     

傅里叶有限差分法三维波动方程正演模拟

傅里叶有限差分(FFD)法兼有相位屏法和隐式有限差分法二者的优势,能够处理复杂地质构造中的波传播问题,但在三维情形下,算子的双向分裂会引起明显的方位各向异性误差.本文用Fourier变换计算双向分裂过程中的高阶交叉项,消除了方位各向异性误差.该方法充分利用了FFD法在双域实现的算法结构,明显减少了由于引入误差校正所带来的计算量.将该方法应用于修改后的三维French模型的地震正演问题,并将得到的叠后记录、单炮记录同全波有限差分法的模拟结果进行对比,结果证实了该方法对一次反射波具有较高的模拟精度,在内存需求和计算效率方面则具有更大的优势.     

起伏地表三维电阻率法有限差分数值模拟

有限差分法是地球物理数值模拟中最常用的方法之一,为了研究起伏地表对三维电阻率法的影响问题,本文基于网格变换方法将起伏地表映射成水平地表,推导出映射后坐标系下稳定电流场基本方程.采用有限差分法对映射后的控制方程进行数值模拟,通过对三维山谷、山脊地形进行数值模拟,分析不同测量装置的地形影响,结果表明:起伏地表对稳定电流场的影响很大.     

Fracture of a Liquefied Crack and the Physics of Rayleigh Waves

The standard free-surface boundary conditions for in-plane crack dynamics are shown to be identical to the conditions for crack dynamics on a liquefied crack. The surfaces of both the free and liquefied cracks do not separate during faulting and hence the static normal stress is not relaxed by the faulting. A crack with either free or liquid boundary conditions deforms in the transverse direction during slip. It follows that both the free and liquefied cracks may represent solutions to the heat-flow paradox. As an application of the proof, we derive a physical understanding of the properties of harmonic Rayleigh waves on a uniform elastic half-space without solving a cubic equation.     

矿井电阻率成像技术的现状与展望

主要介绍国内外关于矿井电阻率层析成像技术在矿井生产中应用的相关资料,综合论述了矿井电阻率成像技术发展现状,介绍了矿井电阻率层析成像技术的特殊问题、井下勘探主要方法和数值模拟技术,并提出了作者对矿井电阻率成像技术发展趋势的展望．     

Boundary conditions of the dynamic numerical model for the Chinese mainland lithosphere

The numerical dynamic model of Chinese mainland lithosphere's stress and strain field was constructed with elasto-viscous creep constitutive relation. The most recent data of the stress field of Chinese mainland and the horizontal movement velocity of the crust blocks were used as constraint conditions. The values of the boundary force were computed by trial-and-error method. The effect of the Qinghai-Xizang Plateau's excess potential energy to the movement of Chinese mainland was studied also in this model. The results of the numerical computing show that, recent rapid raising of the Qinghai-Xizang Plateau and the generation of normal faults in the southern part of the plateau resulted from the convergence of Indian Plate to the Eurasian Plate, and also from the rapid convective thinning of the lower lithosphere. Horizontal extension of the Qinghai-Xizang Plateau is the main dynamic factor to form the present tectonic framework of the Chinese mainland. The compressive loads on the eastern boundary of the model were mainly applied by the subduction of the Pacific Plate. The compression from the Philippine Plate maybe slight.     

低速带的波场数值模拟及影响因素初步分析

低速带具有“自由表面、低速度、高吸收”的特征,能引起反射波能量减弱,产生很强的频散,造成近地表波场畸变,严重影响地震勘探的精度。而地震采集过程又不可避免遇到以下低速带问题:①低速带存在与否会对地震波传播造成什么影响;②近地表变化大,怎么才能得到精确的低速带参数;③存在低速带时,激发井深、最佳的药量(与子波频率有关)、最佳的耦合(与低速带速度相关)等采集参数会对采集资料的品质产生什么影响。本文通过数值模拟,集中对低速带与地震波激发的井深、速度、子波主频等核心采集参数的关系进行分析研究,同时在低速带背景下对数值模拟的网格、边界等问题也进行了讨论。模拟得到的波场特征及低速带参数影响的机理分析,对做好野外地震波激发,获取高品质的地震资料,具有理论意义和指导作用。