In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem. 相似文献
The factors affecting permeability change under repeated mining of coal seams are important study aspects that need to be explored. This study combined various stress variation characteristics of protective seam mining and simplified the stress path of repeated mining in protective seam mines. Based on the results from the bespoke gas flow and displacement testing apparatus, seepage tests for simulated repetitive mining were carried out. The results simulated the actual behavior very well. With any drastic increase in the mining influence, the axial deviation stress in the stress path increased, and the greater the difference in coal permeability during the unloading and stress recovery stage, the more substantial the increase in permeability. The change in coal permeability was significantly influenced by the severity of simulated repeated mining cycles. When the mining stress exceeded a critical value, the permeability of the coal sample increased with the increase in the number of loading and unloading cycles, but the reverse was true when the mining stress was lower than the critical value. The effective sensitivity of seepage to the applied stress decreased with an increase in the number of stress cycles. With a decrease in the deviation stress, that is, with lower severity of mining influence, the effective sensitivity of coal seepage to stress gradually decreased.
Natural Resources Research - In the past few decades, a variety of data-driven predictive modeling techniques has led to a dramatic advancement in mineral prospectivity mapping (MPM). The random... 相似文献