The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid colu...In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.展开更多
Patterns of ground movement and pore water pressure variation are obtained through a case study using a finite element method. With the progress of excavation, ground subsidence, ground inclination and horizontal disp...Patterns of ground movement and pore water pressure variation are obtained through a case study using a finite element method. With the progress of excavation, ground subsidence, ground inclination and horizontal displacement accelerates. Along the striking direction, a subsidence basin is formed on the ground surface induced by underground mining. The maximum subsidence is around 5.41m. The ratio of ground subsidence to the thickness of the coal seam is 1.08. The maximum inclination is 11.5 mm/m. The maximum horizontal displacement is 2.15 mm/m. At the time the coal has been excavated, the maximum pore water pressure reaches 25 kPa. In order to improve protection of structures lo- cated over the area with underground mining, the variation of additional stresses of a railway bridge induced by ground surface deformation is analyzed. The main effect of underground mining on the railway bridge is the tensile stress and the maximum value reaches as high as 4.29 MPa, which is greater than the concrete tensile strength.展开更多
Combined with highway construction, the analysis on the relationship betweentunnel construction and coal resource exploitation was processed, which was based onthe research of rational exploitation depth of coal.3D FE...Combined with highway construction, the analysis on the relationship betweentunnel construction and coal resource exploitation was processed, which was based onthe research of rational exploitation depth of coal.3D FEM numerical analysis for tunnelexcavation was carried out according to engineering geological features of coal measurestrata in the project area.Based on the analysis of displacement and stress of the surroundingrock in the tunnel after excavation, the characteristics for displacement andstress of the tunnel support structure were analyzed when the underlying coal bed wasexploited with sublevel and full caving method.In addition, combined with the related codeand standard, the economic and safe prohibiting exploited depth of the underlying coalbed was proposed, so that a scientific basis for tunnel construction of coal measure strataand reasonable exploitation of the mineral resources in complex geological conditions canbe offered.展开更多
The Zipingpu concrete-faced rockfill dam(CFRD)experienced strong ground motion from the 2008 Wenchuan earthquake.Separation between concrete face slabs and the cushion layer was observed after the earthquake.The separ...The Zipingpu concrete-faced rockfill dam(CFRD)experienced strong ground motion from the 2008 Wenchuan earthquake.Separation between concrete face slabs and the cushion layer was observed after the earthquake.The separation voids under the stage III slabs make up 55%of the total area of the stage III slabs.The observed maximum height of the separation voids was nearly 23 cm at the top of the stage III slabs.Separation voids were also observed locally below the top of stage II slabs near the left abutment,with a maximum height of 7 cm.In this study,a static and dynamic elasto-plastic finite element analysis on Zipingpu CFRD was conducted to capture the separation during the Wenchuan earthquake.The rockfill materials were described using a state-dependent elasto-plastic model that considered particle breakage.The model parameters of rockfill materials were obtained from feedback analysis.The numerical results were largely consistent with the field measurements during construction and after the Wenchuan earthquake.A three-dimensional state-dependent elasto-plastic model that can trace the separation and re-contact of a soil-structure interface was employed to investigate the interaction between concrete face slabs and a cushion layer.The analysis showed the distribution of separation voids observed in the Zipingpu CFRD has a close relationship to the water level and slab dislocations at the time of the earthquake.The phenomenon of the separation from the Wenchuan earthquake was successfully captured by the proposed numerical procedure.展开更多
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
基金supported by NSFC(No.41174118)one of the major state S&T special projects(No.2008ZX05020-004)+1 种基金a Postdoctoral Fellowship of China(No.2013M530106)China Scholarship Council(No.2010644006)
文摘In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.
基金Project 50574088 supported by the Nation Natural Science Foundation of China
文摘Patterns of ground movement and pore water pressure variation are obtained through a case study using a finite element method. With the progress of excavation, ground subsidence, ground inclination and horizontal displacement accelerates. Along the striking direction, a subsidence basin is formed on the ground surface induced by underground mining. The maximum subsidence is around 5.41m. The ratio of ground subsidence to the thickness of the coal seam is 1.08. The maximum inclination is 11.5 mm/m. The maximum horizontal displacement is 2.15 mm/m. At the time the coal has been excavated, the maximum pore water pressure reaches 25 kPa. In order to improve protection of structures lo- cated over the area with underground mining, the variation of additional stresses of a railway bridge induced by ground surface deformation is analyzed. The main effect of underground mining on the railway bridge is the tensile stress and the maximum value reaches as high as 4.29 MPa, which is greater than the concrete tensile strength.
基金Supported by the National Natural Science Foundation Special Originality Innovation Research Colony of China(50621403)
文摘Combined with highway construction, the analysis on the relationship betweentunnel construction and coal resource exploitation was processed, which was based onthe research of rational exploitation depth of coal.3D FEM numerical analysis for tunnelexcavation was carried out according to engineering geological features of coal measurestrata in the project area.Based on the analysis of displacement and stress of the surroundingrock in the tunnel after excavation, the characteristics for displacement andstress of the tunnel support structure were analyzed when the underlying coal bed wasexploited with sublevel and full caving method.In addition, combined with the related codeand standard, the economic and safe prohibiting exploited depth of the underlying coalbed was proposed, so that a scientific basis for tunnel construction of coal measure strataand reasonable exploitation of the mineral resources in complex geological conditions canbe offered.
基金supported by the National Natural Science Foundation of China(Grant Nos.51279025,51379028&91215301)the Program for New Century Excellent Talents in University(Grant No.NCET-12-0083)China’s Post-doctoral Science Fund(Grant No.2015M580225)
文摘The Zipingpu concrete-faced rockfill dam(CFRD)experienced strong ground motion from the 2008 Wenchuan earthquake.Separation between concrete face slabs and the cushion layer was observed after the earthquake.The separation voids under the stage III slabs make up 55%of the total area of the stage III slabs.The observed maximum height of the separation voids was nearly 23 cm at the top of the stage III slabs.Separation voids were also observed locally below the top of stage II slabs near the left abutment,with a maximum height of 7 cm.In this study,a static and dynamic elasto-plastic finite element analysis on Zipingpu CFRD was conducted to capture the separation during the Wenchuan earthquake.The rockfill materials were described using a state-dependent elasto-plastic model that considered particle breakage.The model parameters of rockfill materials were obtained from feedback analysis.The numerical results were largely consistent with the field measurements during construction and after the Wenchuan earthquake.A three-dimensional state-dependent elasto-plastic model that can trace the separation and re-contact of a soil-structure interface was employed to investigate the interaction between concrete face slabs and a cushion layer.The analysis showed the distribution of separation voids observed in the Zipingpu CFRD has a close relationship to the water level and slab dislocations at the time of the earthquake.The phenomenon of the separation from the Wenchuan earthquake was successfully captured by the proposed numerical procedure.