Forward modeling of elastic wave propagation in porous media has great importance for understanding and interpreting the influences of rock properties on characteristics of seismic wavefield. However,the finite-differ...Forward modeling of elastic wave propagation in porous media has great importance for understanding and interpreting the influences of rock properties on characteristics of seismic wavefield. However,the finite-difference forward-modeling method is usually implemented with global spatial grid-size and time-step; it consumes large amounts of computational cost when small-scaled oil/gas-bearing structures or large velocity-contrast exist underground. To overcome this handicap,combined with variable grid-size and time-step,this paper developed a staggered-grid finite-difference scheme for elastic wave modeling in porous media. Variable finite-difference coefficients and wavefield interpolation were used to realize the transition of wave propagation between regions of different grid-size. The accuracy and efficiency of the algorithm were shown by numerical examples. The proposed method is advanced with low computational cost in elastic wave simulation for heterogeneous oil/gas reservoirs.展开更多
Terrain characteristics can be accurately represented in spectrum space. Terrain spectra can quantitatively reflect the effect of topographic dynamic forcing on the atmosphere. In wavelength space, topographic spectra...Terrain characteristics can be accurately represented in spectrum space. Terrain spectra can quantitatively reflect the effect of topographic dynamic forcing on the atmosphere. In wavelength space, topographic spectral energy decreases with decreasing wavelength, in spite of several departures. This relationship is approximated by an exponential function. A power law relationship between the terrain height spectra and wavelength is fitted by the least-squares method, and the fitting slope is associated with grid-size selection for mesoscale models. The monotonicity of grid size is investigated, and it is strictly proved that grid size increases with increasing fitting exponent, indicating that the universal grid size is determined by the minimum fitting exponent. An example of landslide-prone areas in western Sichuan is given, and the universal grid spacing of 4.1 km is shown to be a requirement to resolve 90% of terrain height variance for mesoscale models, without resorting to the parameterization of subgrid-scale terrain variance. Comparison among results of different simulations shows that the simulations estimate the observed precipitation well when using a resolution of 4.1 km or finer. Although the main flow patterns are similar, finer grids produce more complex patterns that show divergence zones, convergence zones and vortices. Horizontal grid size significantly affects the vertical structure of the convective boundary layer. Stronger vertical wind components are simulated for finer grid resolutions. In particular, noticeable sinking airflows over mountains are captured for those model configurations.展开更多
基金supported by the National Basic Research Program of China (No. 2013CB228604)the National Science and Technology Major Project (No. 2011ZX05030-004-002,2011ZX05019-003)the National Natural Science Foundation (No. 41004050)
文摘Forward modeling of elastic wave propagation in porous media has great importance for understanding and interpreting the influences of rock properties on characteristics of seismic wavefield. However,the finite-difference forward-modeling method is usually implemented with global spatial grid-size and time-step; it consumes large amounts of computational cost when small-scaled oil/gas-bearing structures or large velocity-contrast exist underground. To overcome this handicap,combined with variable grid-size and time-step,this paper developed a staggered-grid finite-difference scheme for elastic wave modeling in porous media. Variable finite-difference coefficients and wavefield interpolation were used to realize the transition of wave propagation between regions of different grid-size. The accuracy and efficiency of the algorithm were shown by numerical examples. The proposed method is advanced with low computational cost in elastic wave simulation for heterogeneous oil/gas reservoirs.
基金supported by the Key Research Program of the Chinese Academy of Sciences (Grant No. KZZD-EW-05-01)the special grant (Grant No. 41375052) from the National Natural Science Foundation of Chinafunded by an open project of the State Key Laboratory of Severe Weather (Grant No. 2013LASW-A06)
文摘Terrain characteristics can be accurately represented in spectrum space. Terrain spectra can quantitatively reflect the effect of topographic dynamic forcing on the atmosphere. In wavelength space, topographic spectral energy decreases with decreasing wavelength, in spite of several departures. This relationship is approximated by an exponential function. A power law relationship between the terrain height spectra and wavelength is fitted by the least-squares method, and the fitting slope is associated with grid-size selection for mesoscale models. The monotonicity of grid size is investigated, and it is strictly proved that grid size increases with increasing fitting exponent, indicating that the universal grid size is determined by the minimum fitting exponent. An example of landslide-prone areas in western Sichuan is given, and the universal grid spacing of 4.1 km is shown to be a requirement to resolve 90% of terrain height variance for mesoscale models, without resorting to the parameterization of subgrid-scale terrain variance. Comparison among results of different simulations shows that the simulations estimate the observed precipitation well when using a resolution of 4.1 km or finer. Although the main flow patterns are similar, finer grids produce more complex patterns that show divergence zones, convergence zones and vortices. Horizontal grid size significantly affects the vertical structure of the convective boundary layer. Stronger vertical wind components are simulated for finer grid resolutions. In particular, noticeable sinking airflows over mountains are captured for those model configurations.
