共查询到20条相似文献，搜索用时 40 毫秒
1.
为使物理模拟实验效果与实际探测情形更为接近,在室外自然条件下建立中尺度土壤石油污染实验模型,油污染区扩展深度超过1 m,采用实地探测中常用的500 MHz雷达天线进行长期定时探测.通过实测雷达图像特征、土壤含水量含油量分析,并对比前人开展的小尺度室内模拟试验结果,综合评价探地雷达对油污染区的探测效果.研究表明探地雷达探测图像异常特征与污染区扩散阶段密切相关:包气带内油污染区会引起振幅增强;毛细带的油污染区则表现为水位面反射轴附近清晰可辨的高幅异常区,且水位面反射轴呈下凹状;随扩散过程持续进行,异常区下移与水位面反射轴相交,并产生水平扩张.当污染土含油饱和度大于20%时,可通过雷达图像异常区圈定污染范围;当污染土含油饱和度大于15%时,可通过频谱图出现低频响应的位置圈定污染区水平范围.中尺度实验结果与室内小尺度模拟结果具有一致性,可作为油污染区雷达图像异常的解译依据. 相似文献
2.
In this paper, we present an improvement to our previously published nearly analytic discrete method (NADM) to solve acoustic
and elastic wave equations. We compare the numerical errors of the improved NADM with the original NADM and also the fourth-order
Lax–Wendroff correction and present examples of three-component wave fields in 2D transversely isotropic media with strong
velocity contrasts. Comparing with the original NADM, we find that the improved method requires significantly less storage
space and can increase the time accuracy from second order of the original NADM to fourth order, while the space accuracy
remains the same as that of the original one. Theoretical analyses and numerical results suggest that our improved NADM is
suitable for large-scale numerical modelling as it can effectively suppress numerical dispersion and source-generated noises
caused by discretising wave equations when too-coarse grids are used. 相似文献
3.
4.
合成三维横向非均匀介质远震体波接收函数的Maslov方法 总被引：5，自引：3，他引：2
三维横向非均匀介质远震体波接收函数的合成算法对于三维介质远震体波接收函数研究至关重要．本文发展了基于Ｍａｓｌｏｖ理论的横向非均匀介质接收函数合成算法,针对远震体波接收函数计算的特点,利用遗传算法完成三维动力学射线追踪,并采用递归算法组织射线的追踪过程,从而在保证追踪精度的同时提高了射线追踪的效率,并可完全避免传统射线追踪过程中的射线编码问题．正演计算结果表明,此法是一种高效、实用的三维横向非均匀介质接收函数计算方法． 相似文献
5.
震源定位是微地震监测技术要解决的主要问题.目前,井下微地震监测多采用走时拟合法计算震源位置.常规方法受到环境噪声、初至拾取误差、速度模型误差等因素的影响,定位结果存在一定误差.为了提高定位精度,本文提出了一种改进的基于网格搜索的微地震震源定位方法.本文方法根据P波的偏振特征参数计算概率密度函数求取震源方位角,并采用改进的目标函数和搜索算法计算震源的径向距离和深度.模型数据和实际资料的处理结果表明,本文方法具有较强的抗噪性,计算得到的震源方位角更加接近真实值;与常规目标函数相比,本文方法采用的目标函数具有更好的收敛性,其定位结果受初至拾取误差和速度模型误差的影响更小;本文提出的搜索算法能够消除由于错误拾取造成的观测到时中的异常值对定位结果的影响. 相似文献
6.
复杂速度模型的地震交切定位方法 总被引：3，自引：0，他引：3
交切法是最基本的地震定位方法之一,该定位方法具有稳定性强、速度快的优点。由于速度模型假设为均匀或水平均匀介质,与实际介质有较大偏差,因此定位精度较低,特别是震源深度。为克服传统交切法的缺点,我们对其进行了改进。在改进的方法中,震源轨迹不再假设为圆形或双曲线形,而是利用最小走时树射线追踪方法进行精确计算。数值模型计算表明,本文提出的改进方法,可用于复杂速度模型的地震定位。 相似文献
7.
Patrick N.J. Rasolofosaon 《Geophysical Prospecting》2010,58(4):637-655
Seismic anisotropy in geological media is now widely accepted. Parametrizations and explicit approximations for the velocities in such media, considered as purely elastic and moderately anisotropic, are now standards and have even been extended to arbitrary types of anisotropy. In the case of attenuating media, some authors have also recently published different parametrizations and velocity and attenuation approximations in viscoelastic anisotropic media of particular symmetry type (e.g., transversely isotropic or orthorhombic). This paper extends such work to media of arbitrary anisotropy type, that is to say to triclinic media. In the case of homogeneous waves and using the so‐called ‘correspondence principle’, it is shown that the viscoelastic equations (for the phase velocities, phase slownesses, moduli, wavenumbers, etc.) are formally identical to the corresponding purely elastic equations available in the literature provided that all the corresponding quantities are complex (except the unit vector in the propagation direction that remains real). In contrast to previous work, the new parametrization uses complex anisotropy parameters and constitutes a simple extension to viscoelastic media of previous work dealing with non‐attenuating elastic media of arbitrary anisotropy type. We make the link between these new complex anisotropy parameters and measurable parameters, as well as with previously published anisotropy parameters, demonstrating the usefulness of the new parametrization. We compute the explicit complete directional dependence of the exact and of the approximate (first and higher‐order perturbation) complex phase velocities of the three body waves (qP, qS1 and qS2). The exact equations are successfully compared with the ultrasonic phase velocities and phase attenuations of the three body waves measured in a strongly attenuating water‐saturated sample of Vosges sandstone exhibiting moderate velocity anisotropy but very strong attenuation anisotropy. The approximate formulas are checked on experimental data. Compared to the exact solutions, the errors observed on the first‐order approximate velocities are small (<1%) for qP‐waves and moderate (<10%) for qS‐waves. The corresponding errors on the quality factor Q are moderate (<6%) for qP‐waves but critically large (up to 160%) for the qS‐waves. The use of higher‐order approximations substantially improves the accuracy, for instance typical maximum relative errors do not exceed 0.06% on all the velocities and 0.6% on all the quality factors Q, for third‐order approximations. All the results obtained on other rock samples confirm the results obtained on this rock. The simplicity of the derivations and the generality of the results are striking and particularly convenient for practical applications. 相似文献
8.
由于地下介质对地震波能量有强烈的吸收作用,降低了地震资料的分辨率以及地震资料处理、解释的精度,因此选择合适的参数对地层吸收衰减情况进行有效描述是对地下构造进行高分辨率、高精度处理及解释的有效途径。目前品质因子Q是对岩石弹性参数进行有效描述的重要参数之一。用品质因子Q来衡量地震波能量在介质中传播的吸收衰减情况,是目前提高地震资料分辨率的有效方法,同时也是提高含油气地层解释精度的有效途径。近些年,许多学者提出了多种计算Q值的方法,其中谱比法是在实际中应用最为广泛的Q值估计方法。本文在前人方法的基础上,对谱比法进行了改进,并结合走时层析方法反演Q值,充分利用谱比法在精度、稳定性等方面的优势以及走时层析方法计算效率高的优点,不仅可以提高反演的精度而且能够保证反演的稳定性,提高计算效率。通过模型试算证明谱比法衰减层析方法能有效估计Q值变化情况,具有较好的发展前景。 相似文献
9.
小生境遗传算法在地震CT中的应用研究 总被引：3，自引：1，他引：2
针对含单相饱和流体的孔隙介质,综合费马原理、等效孔隙介质理论、Radon变换等,提出了基于地震CT技术求解速度、孔隙度和岩石密度的小生境多参数反演算法．以此为基础,经过改进,形成了物性参数的隔离小生境反演算法．理论分析与实例计算结果表明,该方法简化了波在层状介质中的传播方程,具有一定的实用价值,可以有效解决井间测井地层问题,获得较高精度的地层物性参数;而且模型简单,有精度高、多峰优化等特点,可用于地下层状介质中油气的寻找,满足油气预测与勘探开发的要求． 相似文献
10.
Modeling and RTM for arbitrary-order pure acoustic wave equation in VTI media using normalized pseudo-analytical method 
下载免费PDF全文
下载免费PDF全文 The conventional pseudo-acoustic wave equations(PWEs) in vertical transversely isotropic(VTI)media may generate SV-wave artifacts and propagation instabilities when anisotropy parameters cannot satisfy the pseudo-acoustic assumption. One solution to these issues is to use pure acoustic anisotropic wave equations, which can produce stable and pure P-wave responses without any SVwave pollutions. The commonly used pure acoustic wave equations(PAWEs) in VTI media are mainly derived from the decoupled P-SV dispersion relation based on first-order Taylor-series expansion(TE), thus they will suffer from accuracy loss in strongly anisotropic media. In this paper, we adopt arbitrary-order TE to expand the square root term in Alkhalifah's accurate acoustic VTI dispersion relation and solve the corresponding PAWE using the normalized pseudoanalytical method(NPAM) based on optimized pseudodifferential operator. Our analysis of phase velocity errors indicates that the accuracy of our new expression is perfectly acceptable for majority anisotropy parameters. The effectiveness of our proposed scheme also can be demonstrated by several numerical examples and reverse-time migration(RTM) result. 相似文献
11.
本文在孔隙弹性介质理论的指导下,基于入射角AVO近似方程推导了包含Russell流体项的两项AVO近似方程和相应的弹性阻抗方程,通过分析可知其精度符合反演要求.在贝叶斯理论框架下,建立了包含正则化约束的弹性阻抗反演方法,在此基础上直接提取Russell流体项.该方法可在缺少大角度叠前地震资料的情况下进行叠前直接反演得到流体因子,减少传统方法带来的累积误差.模型试算表明,该方法具有较好的准确度和稳定性.实际工区应用取得了良好的效果,表明该方法有实用性. 相似文献
12.
Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data. 相似文献
13.
在地震波动模拟中计入常Q滞弹性阻尼,可有效降低模拟波形的误差.就时域有限差分和有限元模拟而言,常基于广义标准线性体建立阻尼介质的时域本构逼近.广义标准线性体由若干标准线性体并联得到,增加标准线性体个数能有效提高模拟精度,但计算量及计算存储将成倍增长.目前尚未有普适的标准线性体个数优化取值方案.本文基于广义标准线性体参数的非线性优化拟合方法,详细分析了时域本构逼近误差的影响因素,清楚揭示了逼近误差仅取决于频带宽度,与频带上下限取值无关这一特性,阐明了构建具有普适性标准线性体个数优化取值方案的可行性.论证了波形模拟精度主要取决于波传播距离与模拟波长的比值以及标准线性体的个数取值.综合考虑上述两个控制因素,结合在波动正反演问题中广为采纳的波形时频误差衡量准则,对不同Q值介质给出了标准线性体个数优化取值表.进一步,本文提出采用不同个数标准线性体以近似不同Q值的阻尼介质时域本构,解决了以往波动数值模拟中统一采用相同个数标准线性体而导致的计算量及计算存储浪费或模拟精度低下等问题,并基于数值实验验证了这一方法的精度.本文工作对推进滞弹性介质波动数值模拟及其在全波形反演问题中的应用具有理论价值和实践意义. 相似文献
14.
Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement. 相似文献
15.
We use the slowness-azimuth station correction （SASC） method to improve the location accuracy of the Wenchuan aftershocks recorded by the Nagqu and Hotan seismic arrays. The results show that the standard deviations of back-azimuth and slowness errors of Wenchuan aftershocks recorded by the Nagqu array decreased by 32% and 58% respectively after correction. The decrease is 38 % and 71% for the Hotan array. After the correction, the location accuracy of all Wenchuan aftershocks recorded by the Nagqu array is improved. For the Hotan array, the accuracy is improved in the slowness estimation for 78 % of aftershocks and in back-azimuth estimation for all aftershocks. 相似文献
16.
Numerical simulation of seismic wavefields in TTI media using the rotated staggered-grid compact finite-difference scheme 下载免费PDF全文
下载免费PDF全文 Numerical simulation in transverse isotropic media with tilted symmetry axis(TTI) using the standard staggered-grid finite-difference scheme(SSG)results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme(RSG) is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme(CSG) is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme(RSGC). The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG. 相似文献
17.
Luděk Klimeš 《Pure and Applied Geophysics》1996,148(3-4):539-563
A new computational scheme for calculating the first-arrival travel times on a rectangular grid of points is proposed. The new proposed method is of second-order accuracy. This means that the error of the calculated travel time is proportional to the second power of the grid spacing. The method should be sufficiently accurate for all applications in smooth seismic models. On the other hand, the method is not, in its present form, proposed for models with structural interfaces which make the method unstable and generate travel-time errors of the first order. Equations are also presented for the appropriate evaluation of the errors of calculated travel times to check their accuracy, and the proposed method is compared with other numerical methods. The method is developed, described and demonstrated in 2-D, but may also be extended to 3-D models and to general models with structural interfaces. 相似文献
18.
Manuel Stark Tobias Heckmann Livia Piermattei Fabian Dremel Andreas Kaiser Patrick Machowski Florian Haas Michael Becht 《地球表面变化过程与地形》2021,46(10):2019-2043
Uncrewed aerial systems (UAS), combined with structure-from-motion photogrammetry, has already proven to be very powerful for a wide range of geoscience applications and different types of UAS are used for scientific and commercial purposes. However, the impact of the UAS used on the accuracy of the point clouds derived is not fully understood, especially for the quantitative analysis of geomorphic changes in complex terrain. Therefore, in this study, we aim to quantify the magnitude of systematic and random error in digital elevation models derived from four commonly used UAS (XR6/Sony α6000, Inspire 2/X4s, Phantom 4 Pro+, Mavic Pro) following different flight patterns. The vertical error of each elevation model is evaluated through comparison with 156 GNSS reference points and the normal distribution and spatial correlation of errors are analysed. Differences in mean errors (−0.4 to −1.8 cm) for the XR6, Inspire 2 and Phantom 4 Pro are significant but not relevant for most geomorphological applications. The Mavic Pro shows lower accuracies with mean errors up to 4.3 cm, thus showing a higher influence of random errors. QQ plots revealed a deviation of errors from a normal distribution in almost all data. All UAS data except Mavic Pro exhibit a pure nugget semivariogram, suggesting spatially uncorrelated errors. Compared to the other UAS, the Mavic Pro data show trends (i.e. differences increase with distance across the survey—doming) and the range of semivariances is 10 times greater. The lower accuracy of Mavic Pro can be attributed to the lower GSD at the same flight altitude and most likely, the rolling shutter sensor has an effect on the accuracy of the camera calibration. Overall, our study shows that accuracies depend highly on the chosen data sampling strategy and that the survey design used here is not suitable for calibrating all types of UAS camera equally. 相似文献
19.
基于地震波反射系数近似公式的叠前反演是油气勘探的重要工具.本文在已有研究的基础上,推导了一个改进的射线参数域地震纵波反射系数近似方程.该方程建立了地震纵波反射系数与纵波阻抗和横波阻抗的非线性关系,在中、小角度的范围内较现有的反射系数线性近似公式精度更高.另外,由于该方程仅包含纵波和横波阻抗反射系数项,因此基于新方程的反演能够有效地降低同步反演纵波速度、横波速度、密度三个参数的不适定性.在此基础上,结合广义线性反演法(GLI)理论和贝叶斯理论,相应地发展了一种叠前地震同步反演方法.模型测试和实际资料的应用表明,基于新方程的反演方法能够利用有限角度(偏移距)的数据稳定地反演纵波和横波阻抗,由于在反演过程中,不需要假设纵横波速度为常数,因此该方法还能有效地提高反演结果的精度. 相似文献
20.
地震定位是地震监测与减灾研究重要基础.基于均匀或横向均匀介质模型,利用震源轨迹确定震源位置的交切法具有稳健和效率高的优点,但定位精度较低,特别是震源深度.为提高震源定位精度,我们提出适用于三维复杂速度模型的地震定位交切法.将地壳速度模型由均匀或横向均匀介质模型扩展为三维复杂速度模型;均匀或横向均匀介质模型对应的原假设为球面或双曲面的震源轨迹通过最小走时树射线追踪技术予以确定.确定震源位置的震源轨迹以到时差作为约束条件;将震源定位于震源轨迹交汇最密集的点处,即总的到时差残差(RDT)最小的点处.定位结果的不确定性可通过RDT值较小节点的空间分布予以定性表示.考察了准确速度模型、扰动速度模型、扰动观测到时及地震在台网外等4种情况下改进方法的地震定位效果,结果表明改进的交切法可用于三维复杂速度模型的地震定位;综合利用P波与S波的到时差信息,可明显改善震源位置约束;使用多条震源轨迹进行定位,有助于减少由随机因素导致的定位误差. 相似文献