排序方式: 共有70条查询结果,搜索用时 15 毫秒
51.
传统有限差分系数是通过泰勒级数展开求取的,这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散,本文发展了一种新的优化差分系数方法:首先将泰勒级数展开与多点采样方法结合应用于空间频散关系,基于最大范数建立直观有效的优化目标函数,采用Remez算法求解该目标函数,从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性(VTI)介质中的声波和弹性波方程.另外,本文将二维混合吸收边界条件推广到三维VTI介质中,用于吸收人工截断边界反射;基于各向异性特征,合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题,本文波场外推过程中采用图形处理器(GPU)取代传统的中央处理器(CPU).数值精度分析表明,相比较于传统的泰勒级数展开方法,优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率,混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果. 相似文献
52.
Two seismic modelling approaches, that is, two-dimensional pre-stack elastic finite-difference and one-dimensional convolution methods, are compared in a modelling exercise over the fluid-flow simulation model of a producing deep-water turbidite sandstone reservoir in the West of Shetland Basin. If the appropriate parameterization for one-dimensional convolution is used, the differences in three-dimensional and four-dimensional seismic responses from the two methods are negligible. The key parameters to ensure an accurate seismic response are a representative wavelet, the distribution of common-depth points and their associated angles of incidence. Conventional seismic images generated by the one-dimensional convolutional model suffer from lack of continuity because it only accounts for vertical resolution. After application of a lateral resolution function, the convolutional and finite-difference seismic images are very similar. Although transmission effects, internal multiples and P-to-S conversions are not included in our convolutional modelling, the subtle differences between images from the two methods indicates that such effects are of secondary nature in our study. A quantitative comparison of the (normalized root-mean-square) amplitude attributes and waveform kinematics indicates that the finite-difference approach does not offer any tangible benefit in our target-oriented seismic modelling case study, and the potential errors from one-dimensional convolution modelling are comparatively much smaller than the production-induced time-lapse changes. 相似文献
53.
弹性波矢量波场逆时偏移可以综合利用纵横波场信息,对地下空间进行清晰成像,且对于成像介质没有角度限制,可以对复杂构造进行更清晰的成像.而弹性波逆时偏移中最重要的就是求解波动方程的算法,其直接影响成像的精度以及效率.本文引进电力系统谐波分析中常用的余弦组合窗函数,并通过一种新的优化算法得到了改进的余弦组合窗函数从而得到优化后的有限差分算子.并将此算子应用于解耦方程的矢量波场分离算法从而提高了成像精度.数值测试表明基于新算法的逆时偏移的成像精度和清晰度得到了明显的提高. 相似文献
54.
55.
地震波传播有限差分模拟的人工边界问题 总被引:6,自引:3,他引:6
在用有限差分法模拟地震波传播中,差分网格人工边界的反射问题是至关的。本文简要介绍了几种消除人工边界反射的方法及其主要优缺点,联立这些不同的吸收边界条件于实际的有限差分计算可以给出好的结果。 相似文献
56.
针对不同陆面模式对土壤湿度方程求解方法以及对土壤分层结构的差异, 本文选取了三类陆面模式(CLM, CABLE和ECMWF陆面模式)作比较研究。为了避免不同模式参数化方案引起的上边界误差, 上边界采用固定蒸发率、 入渗率和固定表层土壤湿度三类边界条件。土壤分层采用101层(细网格)和11层(粗网格)两种, 并考虑土壤性质沿深度变化。结果表明: 当土壤性质均匀时, 求解的差别主要在第三类边界条件下CLM 求出的水分入渗速度比其它两种快; 改用粗网格后由于土壤深层厚度加大无法与细网格得出的土壤湿度廓线相重合。当土壤性质非均匀时, 模拟结果间差别加大, 只有ECMWF模式模拟的土壤湿度廓线是严格连续的。对于模式和上边界的不同组合, 粗、 细网格模拟结果间均方根偏差不一致。一般而言, CABLE模拟的偏差除第一类条件较小外, 其它都是最大的。第二三类边界条件引起的偏差较大, 第一类最小。上述结果提示我们, 在比较不同陆面模式以及用观测资料来检验模拟结果时应充分考虑土壤分层及土壤性质非均匀性的可能影响。 相似文献
57.
The staggered grid finite-difference method is a powerful tool in seismology and is commonly used to study earthquake source
dynamics. In the staggered grid finite-difference method stress and particle velocity components are calculated at different
grid points, and a faulting problem is a mixed boundary problem, therefore different implementations of fault boundary conditions
have been proposed. Viriuex and Madariaga (1982) chose the shear stress grid as the fault surface, however, this method has
several problems: (1) Fault slip leakage outside the fault, and (2) the stress bump beyond the crack tip caused by S waves
is not well resolved. Madariaga et al. (1998) solved the latter problem via thick fault implementation, but the former problem remains and causes a new issue;
displacement discontinuity across the slip is not well modeled because of the artificial thickness of the fault. In the present
study we improve the implementation of the fault boundary conditions in the staggered grid finite-difference method by using
a fictitious surface to satisfy the fault boundary conditions. In our implementation, velocity (or displacement) grids are
set on the fault plane, stress grids are shifted half grid spacing from the fault and stress on the fictitious surface in
the rupture zone is given such that the interpolated stress on the fault is equal to the frictional stress. Within the area
which does not rupture, stress on the fictitious surface is given a condition of no discontinuity of the velocity (or displacement).
Fault normal displacement (or velocity) is given such that the normal stress on the fault is continuous across the fault.
Artificial viscous damping is introduced on the fault to avoid vibration caused by onset of the slip. Our implementation has
five advantages over previous versions: (1) No leakage of the slip prior to rupture and (2) a zero thickness fault, (3) stress
on the fault is reliably calculated, (4) our implementation is suitable for the study of fault constitutive laws, as slip
is defined as the difference between displacement on the plane of z = + 0 and that of z = − 0, and (5) cessation of slip is achieved correctly. 相似文献
58.
本文提出一种结合广义射线和有限差分理论而形成的混合方法(以下简称混合方法).在介质均匀区(包括平行介质分层区)应用广义射线理论,在介质不均匀区应用有限差分理论.在这种结合过程中,把均匀区广义射线的结果作为不均匀区有限差分区波场连续的输入.对于均匀介质模型,混合方法理论地震图和广义射线理论地震图、广义反射透射系数矩阵及离散波数法地震图的对比,取得了满意的结果,在数值上证明了该方法的正确性.计算了1990年8月5日Hakone地震前震在日本Ashigara山谷盲测点的响应,理论地震图和实际资料的比较同样取得好的结果.这种方法可用于地震勘探、工程地震、强地面运动的研究. 相似文献
59.
Both finite-element and finite-difference numerical models are applied to simulate storm surges and associated currents generated
by tropical cyclones that struck the coast of Andhra Pradesh, located on the east coast of India. During a cyclone, the total
water level at any location on the coast is made up of the storm surge, surge–wind wave interaction and the tide. The advanced
circulation two-dimensional depth-integrated (ADCIRC-2DDI) model based on finite-element formulation and the two-dimensional
finite-difference model of storm surges developed at IIT Delhi, hereafter referred as IITD storm surge model, are used. These
models are driven by astronomical tides at the open ocean boundary and cyclonic asymmetric winds over the surface of the computational
domain. Comparison of model simulated sea-surface elevations with coarse and finer spatial resolutions suggests that the grid
resolution near the coast is very crucial for accurate determination of the surges in addition to the local bathymetry. The
model underpredicts surges, and the peak surge location shifts more to the right of the landfall as the spatial resolution
of the model becomes coarser. The numerical experiments also demonstrate that the ADCIRC model is robust over the IITD storm
surge model for surge computations as the coastline is better represented in the former. 相似文献
60.
The numerical error associated with finite-difference simulation of wave propagation in discontinuous media consists of two
components. The first component is a higher-order error that leads to grid dispersion; it can be controlled by higher-order
methods. The second component results from misalignment between numerical grids and material interfaces. We provide an explicit
estimate of the interface misalignment error for the second order in time and space staggered finite-difference scheme applied
to the acoustic wave equation. Our analysis, confirmed by numerical experiments, demonstrates that the interface error results
in a first-order time shift proportional to the distance between the interface and computational grids. A 2D experiment shows
that the interface error cannot be suppressed by higher-order methods and indicates that our 1D analysis gives a good prediction
about the behavior of the numerical solution in higher dimensions.
相似文献